A method for generating seismic attribute gathers, the method including: computing, with a computer, seismic images with a field dataset; generating, with a computer, synthetic data corresponding to the seismic images; computing, with a computer, an attribute volume by applying an expectation method to the synthetic data; mapping, with a computer, the attribute volume to the seismic images; and generating, with a computer, seismic attribute gathers by stacking the seismic images mapped to the attribute volume.

Computing device, computer instructions and method for correcting an image, of a surveyed surface, due to a free-surface reflection. The method includes calculating a free surface reflection operator for a seismic source displaced in water based on a position of the source, and an air-water interface datum; receiving recorded seismic data d recorded with seismic sensors (r), wherein the recorded seismic data is associated with a pressure and/or a particle motion produced by a seismic wave in earth; correcting the recorded seismic data d based on the free surface reflection operator to obtain transformed seismic data; and generating an image of the surveyed subsurface, based on the transformed seismic data, wherein the image is indicative of various layers of the earth. The free surface reflection operator varies while a source signal is being emitted by the source.

Predicting elastic parameters of a subsurface includes modelling changes in the shear modulus and changes in the bulk modulus of the subsurface as a combination of a host medium shear modulus and host medium bulk modulus and a plurality of inclusion shear moduli and inclusion bulk moduli. Each inclusion shear modulus and inclusion bulk modulus associated with a unique inclusion geometry. An inclusion-based rock physical model is used to solve the models for changes in shear modulus and changes in bulk modulus to predict an effective shear modulus of the subsurface and an effective bulk modulus of the subsurface.

Computing device, computer instructions and method for estimating a broadband wavelet associated with a given seismic data set. The method includes receiving broadband seismic data; constructing and populating a misfit function; calculating the broadband wavelet based on the misfit function and the broadband seismic data; and estimating physical reservoir properties of a surveyed subsurface based on the broadband wavelet. The broadband wavelet is constrained, through the misfit function, by (1) an amplitude only long wavelet, and (2) an amplitude and phase short wavelet. The amplitude and phase short wavelet is shorter in time than the amplitude only long wavelet.

A method, including: obtaining, with a computer, an initial geophysical model; modeling, with a computer, a forward wavefield based on the initial geophysical model with wave equations including a second order z-derivative in a rotated coordinate system that accounts for a tilted transverse isotropic (TTI) medium; modeling, with a computer, an adjoint wavefield with adjoint wave equations including a second order z-derivative in a rotated coordinate system that accounts for a tilted transverse isotropic (TTI) medium, wherein the wave equations and the adjoint wave equations include relaxation terms accounting for anelasticity of earth in an update of a primary variable and an evolution relationship for the relaxation terms; and obtaining, with a computer, a gradient of a cost function based on a combination of a model of the forward wavefield and a model of the adjoint wavefield.

A tangible, non-transitory computer-readable medium configured to store instructions executable by a processor of an electronic device to access a beam migration image of a subsurface target. In addition, the computer-readable medium is configured to store instructions executable by a processor of an electronic device to determine a decomposition criteria based on at least one of subsurface dip inclinations, subsurface dip azimuths, or a combination thereof. Further, the computer-readable medium is configured to store instructions executable by a processor of an electronic device to decompose the beam migration image into a plurality of partial images according to the decomposition criteria to provide various views of the subsurface target. The plurality of partial images are usable by seismic interpreters in exploration for hydrocarbons within the subsurface target.

Method for acquiring, at reduced acquisition cost, seismic data using simultaneous, field-encoded sources in the field (702), and then constructing pseudo source-records (703) that better meet the requirements for using additional simultaneous computer-encoded sourcing for computer simulations or forward modeling (706) as part of (707) iterative FWI (Full Wavefield Inversion) or RTM (Reverse Time Migration), with additional reduction in computational costs. By better meeting the requirements of simultaneous sourcing for FWI or RTM (701), artifacts and crosstalk are reduced in the output. The method can be used for marine streamer acquisition and other non-fixed spread geometries to acquire both positive and negative offsets and to mitigate the “missing data” problem for simultaneous-source FWI. It can also be used for land data to overcome issues with moving spreads and long continuous records.

A system and method for analyzing geologic features including fluid estimation and lithology discrimination may include the steps of identifying areas of interest on a seismic horizon, computing statistical data ranges for the seismic amplitudes within the areas of interest, and analyzing the geologic features based on the amplitude variation with offset (AVO) or angle (AVA) curves including the statistical data ranges.

A method for obtaining zero-offset and near zero offset seismic data from a marine survey, with separation of direct arrival information and reflectivity information, the method including: modeling a direct arrival estimate at a passive near-field hydrophone array by using a notional source separation on active near-field hydrophone data; generating reflection data for the passive near-field hydrophone array by subtraction of the modeled direct wave from data recorded by the passive near-field hydrophone array; generating near zero-offset reflectivity traces by stacking the reflection data for the passive near-field hydrophone array on a string-by-string basis or on a combination of strings basis; generating reflectivity information at the active near-field hydrophone array by subtracting the direct arrival estimate modeled using the notional source separation from the active near-field hydrophone data; and generating an estimate of zero-offset reflectivity traces by calculating a cross-correlation between the between the reflectivity information at the active near-field hydrophone array and the near zero-offset traces and performing an optimized stacking with summation weights based on coefficients of the cross-correlation.

A method including selecting a forward model based on a modeled well structure and including a single modeled acoustic source located in a modeled wellbore and a plurality of modeled acoustic sensors located in a modeled source area, simulating an acoustic signal generated by the single modeled acoustic source and received by each modeled acoustic sensor, calculating phases of the simulated acoustic signals received at each modeled acoustic sensor, obtaining with a principle of reciprocity a plurality of modeled acoustic sources in the modeled source area and a single modeled acoustic sensor in the modeled wellbore, calculating phase delays of the simulated acoustic signals between each modeled acoustic source and the single modeled acoustic sensor, detecting acoustic signals generated by a flow of fluid using acoustic sensors in a wellbore, and processing the acoustic signals using the phase delays to generate a flow likelihood map.