The presently disclosed seismic acquisition technique employs a receiver array and a processing methodology that are designed to attenuate the naturally occurring seismic background noise recorded along with the seismic data during the acquisition. The approach leverages the knowledge that naturally occurring seismic background noise moves with a slower phase velocity than the seismic signals used for imaging and inversion and, in some embodiments, may arrive from particular preferred directions. The disclosed technique comprises two steps: 1) determining from the naturally occurring seismic background noise in the preliminary seismic data a range of phase velocities and amplitudes that contain primarily noise and the degree to which that noise needs to be attenuated, and 2) designing an acquisition and processing method to attenuate that noise relative to the desired signal.

Certain implementations of a three-dimensional elastic frequency domain iterative solver for full waveform inversion can be implemented as a method in which frequency domain numerical simulation of elastic waves is per formed in three-dimensional (3D) media.

Certain implementations of a three-dimensional elastic frequency domain iterative solver for full waveform inversion can be implemented as a method in which frequency domain numerical simulation of elastic waves is per formed in three-dimensional (3D) media.

A method for exploring a geological formation estimates a multiple-inclusive source wavelet and a reflectivity model by performing an inversion minimizing the difference between the shot records and modeled shots. The modeled shots, that include primaries and multiples, are obtained by propagating up-going and down-going wavefields through an input velocity model, which is not updated during the inversion. The shot records are converted into a geological product suitable for hydrocarbon exploration in the geological formation using the multiple-inclusive source wavelet and the reflectivity model.

Described herein are implementations of various technologies for a method for seismic data processing. The method may receive seismic data for a region of interest. The seismic data may be acquired in a seismic survey. The method may determine sparse seismic data by selecting shot points in the acquired seismic data using statistical sampling. The method may determine simulated seismic data based on an earth model for the region of interest, a reflection model for the region of interest, and the selected shot points. The method may determine an objective function that represents a mismatch between the sparse seismic data and the simulated seismic data. The method may update the reflection model using the objective function.

Systems and methods for simulating an effect of multiple seismic sources concurrently on a geologic formation are provided. Data is read from a seismic source file that describes at least two seismic source types. The concurrent propagation of acoustic energy from the at least two seismic sources types through the geologic formation is modeled. A seismic output file is then generated.

A method for seismic imaging includes receiving a multi-shot seismic data set that was collected using one or more streamers having recorders configured to detect seismic waves that propagate through a subterranean domain. The method also includes partitioning the multi-shot seismic data set into windows including a source dimension. The method also includes defining one or more first basis functions that describe the windows of the multi-shot seismic data set. The method also includes generating a model that describes a decomposition of the multi-shot seismic data set using the one or more first basis functions. The method also includes defining one or more second basis functions that describe a selected output data. The method also includes combining the one or more second basis functions with the model to produce a result for a source side wavefield and a receiver side wavefield.

A method for seismic processing includes steps of seismic signal forward propagation and seismic data back propagation. The subsurface medium image is created after correlating and summarizing forward and backward propagation results. To address migration footprint and noise due to the incomplete data acquisition aperture and migration approximation in the migration operator, the iteration inversion strategy incorporates tensor flow calculated from seismic image. A regularization operator based on structure tensor of image is applied to seismic image inversion.

Systems and methods of detecting marine seismic survey parameters are provided. A data processing system can obtain seismic data from seismic data acquisition units disposed on a seabed responsive to an acoustic signal propagated from an acoustic source through a water column. The data processing system can determine from the seismic data, a direct arrival time for the acoustic signal at each of the plurality of seismic data acquisition units, and can obtain an estimated depth value of each of the plurality of seismic data acquisition units and an estimated water column transit velocity of the acoustic signal. The data processing system can apply a depth model and a water column transit velocity model to the estimated depth value and to the estimated water column transit velocity determine an updated depth value and an updated water column transit velocity for each of the plurality of seismic data acquisition units.

A multi-stage FWI workflow uses multiple-contaminated FWI models to predict surface-related multiples. A method embodying the present technological advancement, can include: using data with free surface multiples as input into FWI; generating a subsurface model by performing FWI with the free-surface boundary condition imposed on top of the subsurface model; using inverted model from FWI to predict multiples; removing predicted multiples from the measured data; using the multiple-free data as input into FWI with absorbing boundary conditions imposed on top of the subsurface model; and preparing a multiple free data set for use in conventional seismic data processing.