A method for reducing effects of free surface multiple reflections from seismic signal measurements. The measurements result from seismic energy imparted into the Earth's subsurface from collocated measurements related to pressure and vertical component of motion in response to the imparted seismic energy. The method includes entering as input to a computer the measurements related to pressure and vertical component of motion. In the computer, a down-going component of the measurements is determined. An impulse response of the Earth in the absence of a free surface from the down-going component is determined.

Computing device, computer instructions and method process input seismic data d recorded in a first domain by seismic receivers that travel in water, the input seismic data d including pressure and particle motion measurements, including up-going and down-going wave-fields. A model p is generated in a second domain by solving an inverse problem for the input seismic data d, wherein applying an L transform to the model p describes the input data d. An L′ transform, which is different from the L transform, is then applied to the model p to obtain an output seismic data in the first domain, the output seismic data having a characteristic imparted by the transform L′. The characteristic is related to pressure wave-fields and/or particle motion wave-fields interpolated at positions in-between the input seismic receivers. An image of the surveyed subsurface is generated based on the output seismic dataset.

A method for determining seismic sensor depths in a body of water includes accepting as input to a computer measurements of seismic signals made by a plurality of seismic sensors disposed in a body of water. A depth increment and a range of sensor depths for correlation of signals from each of the plurality of seismic sensors is defined. In the computer, the input seismic measurements are extrapolated to each depth increment in the range. A depth of each seismic sensor is determined by correlating the seismic signal measurements with depth-extrapolated measurements of the seismic signal measurements.

A DUnet engine produces a processed image of seismic data acquired over an underground formation. The DUnet engine includes: a contractive path that performs multilayer convolutions and contraction to extract a code from the seismic data input to the DUnet, an expansive path configured to perform multilayer convolutions and expansion of the code, using features provided by the contractive path through skip connections, and a model level that performs multilayer convolutions on outputs of the contractive path and expansive paths to produce the processed image and/or an image that is a difference between the processed image and the seismic data. A fraction of the seismic data may be selected for training the DUnet engine using an anchor method that automatically extends an initial seismic data subset, based on similarity measurements. A reweighting layer may further combine inputs received from layers of the DUnet model to preserve signal amplitude trend.

A sensor system includes a triplet element including a first hydrophone, a second hydrophone, and a third hydrophone configured to receive an incoming signal at a first phase, a second phase, and a third phase, respectively, the first to third hydrophones extending along a first direction, and a processor configured to determine an incidence direction of the incoming signal, and to dynamically generate a cardioid null in the incidence direction to reject the incoming signal based on the incoming signal at the first to third phases.

A method for directional Q compensation of seismic data may comprise calculating angle-dependent subsurface travel times; applying directional Q compensation to the prestack seismic data to obtain Q-compensated data in time-space domain, wherein the directional Q compensation is based on the angle-dependent subsurface travel times; and using the Q-compensated data to generate an image of the subsurface. Directional Q compensation may comprise determining an angle-dependent forward E operator and an angle-dependent adjoint E* operator using the angle-dependent subsurface travel times; and applying a sparse inversion algorithm using the angle-dependent operators to obtain a model of Q-compensated data. The angle-dependent operators may be preconditioned by introducing ghost and source effects in a wavelet matrix and a transpose of the wavelet matrix, respectively, such that applying a sparse inversion algorithm using the preconditioned angle-dependent operators is used to obtain a model of Q-compensated, deghosted data without source effects.

A method for reducing inline directivity of an air-gun source signature by optimizing spatial distribution of air-guns is provided according to the present application, which relates to a field of design and optimization of an air-gun source. An evaluation standard in the air-gun distribution in an air-gun array direction is proposed. By a combination optimization along both the inline and depth directions, a design scheme having evidently broader effective bandwidth and effective take-off angle width than a design scheme of a conventional source is obtained, with which the directivity of the air-gun source signature can be reduced.

Processes and systems for deblending blended seismic data with attenuated source signatures and source ghost are described. Processes and systems compute blended upgoing pressure wavefield based on blended pressure wavefield and blended vertical velocity wavefield recorded in a marine survey of a subterranean formation. Downgoing vertical velocity wavefield is computed based on near-field pressure wavefields generated by source elements of sources activated in the marine survey. Deblended wavefield is computed based on the blended upgoing pressure wavefield and the downgoing vertical velocity source wavefield. The deblended wavefield may be used to generate an image of the subterranean formation with the source signatures and source ghosts contained in the blended pressure wavefield and blended vertical velocity wavefield.

Techniques are disclosed relating to geophysical surveying and data processing using streamers at different depths. In one embodiment, a method includes obtaining geophysical data specific to a geophysical formation that includes a first set of data representative of a first particle motion signal and a first pressure signal recorded using sensors at a first depth and a second set of data representative of a second particle motion signal and a second pressure signal recorded using sensors at a second, greater depth. In this embodiment, the method includes generating first and second wave separation equations from the first and second sets of data and determining a cross-line wavenumber value that reconciles the first and second equations. The determined cross-line wavenumber may be used to separate measured up-going and down-going wavefields.

A method for seismic data processing can include obtaining seismic data acquired based upon trigger times and not based upon positions of triggered source elements. The seismic data can include near-continuously recorded seismic data in split records. The split records can be spliced together into a single near-continuous record to produce a trace with seismic data from a single acquired line. The seismic data can be processed by performing a spatial shift for each of a number of time samples to correct for motion of a number of seismic receivers.