A method and apparatus for separating the unknown contributions of two or more sources from a commonly acquired wave field including the determination of a wavenumber dependent model which reconstructs the wave field of the sources independently below a frequency set by the slowest physical propagation velocity, and applying an inversion based on the model to the commonly acquired wave field to separate the contributions.

A method is described for removing the surface ghost from and/or separating wave field data and/or for estimating redatuming operators of the wave field data by effective use of a transform that that relies on the non-uniform distribution of distances with respect to a reference surface or of tuned-source radiation directions of sources and or the non-uniform distribution of receivers with respect to a reference surface to partition or map the wave field from at least two different cones in the transformed domain and using the contribution of sources and or receivers inside at least one of the at least two different cones to estimate a first wave field of interest, a second separated or ghost wave field and/or redatuming operator.

The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for suppressing noises in seismic data. One computer-implemented method includes receiving, at a data processing apparatus, a set of seismic data associated with a subsurface region; flattening, by the data processing apparatus, the set of seismic data according to an identified seismic event; dividing, by the data processing apparatus, the set of seismic data into a plurality of spatial windows; randomizing, by the data processing apparatus, the set of seismic data according to a random sequential order; filtering, by the data processing apparatus, the randomized seismic data; and reorganizing, by the data processing apparatus, the filtered seismic data according to a pre-randomization order.

Methods and systems that attenuate noise in seismic data. In one aspect, noise attenuation methods iteratively generate a low-speed noise model of low-speed noise recorded in the seismic data. The seismic data is arranged into a sparse seismic-data matrix. Low-speed noise refers to noise that propagates at speeds less the speed of sound in water. The low-speed noise model includes the low-speed noise and includes interpolated low-speed noise that approximates portions of the low-speed noise typically affected by spatial aliasing and streamer surface irregularities. The low-speed noise model may be subtracted from the sparse seismic-data matrix to obtain a noise-attenuated sparse seismic-data matrix.

A method for wavefield separation of sonic data is provided. The method comprises estimating direct phases of waveforms of sonic data observed with two or more sensors by using cross-correlation of waveform traces at adjacent sensor locations, removing the direct phases from the observed waveforms, and extracting event signals from the waveforms after removing the direct phases.

Seismic shot gather data is received from a computer data store for processing. The received seismic shot gather data is separated into downgoing and upgoing wavefields, a time-frequency-wavenumber (t-f-k) three-dimensional (3D) data cube comprising multiple time-frequency (t-f) slices is formed. The downgoing wavefields are wavelet transformed from a time (t) domain to a t-f domain and the upgoing wavefields are wavelet transformed from the t domain to the t-f domain. A wavelet cross-correlation is performed between the downgoing wavefields in the t-f domain and the upgoing wavefields in a t-f-k domain to generate wavelet cross-correlated data. Soft-threshold filtering if performed for each t-f slice of the t-f-k 3D data cube. An inverse wavelet transform is performed to bring wavelet cross-correlated data from the t-f-k domain to a time-receiver (t-x) domain. All seismic shots of the received seismic shot gather data are looped over and the wavelet cross-correlated data is stacked as a virtual source gather.

A technique includes receiving sensor data; sorting the data into a gather representation that corresponds to a plurality of shots of an energy source; and determining a signal cone based at least in part on at least one characteristic of the gather representation. The technique includes processing the sensor data in a processor-based machine to attenuate noise to generate data representing a signal based at least in part on the determined signal cone and the gather representation.

A system and method are provided for filtering scatterer noise energy from land seismic waves using three dimensional (3D) iterative filtering in a cross-spread source-receiver pattern. The system and method isolate scatterer noise energy associated with a scatterer based on a filtering process performed using a scatterer referential time delay. Then, the isolated scatterer noise energy can be subtracted from the surface wave data. This process can be repeated for each scatterer in the covered area to remove each of their contributions, and can be performed, for example, after a preliminary filtering of the surface wave data using 3D fk filtering or the like.

Waves from cement bond logging with a sonic logging-while-drilling tool (LWD-CBL) are often contaminated with tool waves and may yield biased CBL amplitudes. The disclosed LWD-CBL wave processing corrects the first echo amplitudes of LWD-CBL before calculating the BI. The LWD-CBL wave processing calculates a tool wave amplitude and a phase angle difference as the difference of the phases between the tool waves and casing waves. The tool waves are then used to correct the LWD-CBL casing wave amplitude and remove errors introduced from tool waves. In conjunction with the sets of operations described, the LWD-CBL wave processing also include array preprocessing operations. Array preprocessing may employ variation of bandpass filtering and frequency-wavenumber (F-K) filtering operations to suppress tool wave.

Waves from cement bond logging with a sonic logging-while-drilling tool (LWD-CBL) are often contaminated with tool waves and may yield biased CBL amplitudes. The disclosed LWD-CBL wave processing corrects the first echo amplitudes of LWD-CBL before calculating the BI. The LWD-CBL wave processing calculates a tool wave amplitude and a phase angle difference as the difference of the phases between the tool waves and casing waves. The tool waves are then used to correct the LWD-CBL casing wave amplitude and remove errors introduced from tool waves. In conjunction with the sets of operations described, the LWD-CBL wave processing also includes array preprocessing operations. Array preprocessing may employ variation of bandpass filtering and frequency-wavenumber (F-K) filtering operations to suppress tool waves.