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.

A method of determining a shear-wave attenuated vertical component vertical seismic profile (VSP) dataset is disclosed. The method includes, obtaining a multi-component VSP dataset, including a vertical and a horizontal component, transforming the vertical component into a vertical spectrum and the horizontal component into a horizontal spectrum, and designing a band-pass filter based, at least in part, on an energetic signal of the horizontal spectrum. The method further includes determining a muted vertical amplitude spectrum by applying the pass-band filter to an amplitude spectrum of the vertical spectrum, determining an estimated noise model based on the muted vertical amplitude spectrum and the vertical spectrum; and determining the shear-wave attenuated vertical component VSP dataset by adaptively subtracting the estimated noise model from the vertical component of the multi-component VSP dataset. A system including a seismic source, a plurality of seismic receivers, and a seismic processor for executing the method is disclosed.

A system for seismic imaging of a subterranean geological formation, the system includes a receiver configured to obtain seismic data comprising a data volume representing a post-stacked image. The system includes a filtering module configured to: apply frequency-wavenumber (F-K) filter to the data volume extract a negative-dip structure image and apply the frequency-wavenumber (F-K) filter to the data volume extract a positive-dip structure image. The system includes a diffraction rendering module configured to: multiply the positive-dip structure image with the negative-dip structure image and generate a diffraction-enhanced seismic image representing a geological formation of the data volume.

A simultaneous sources seismic acquisition method is described that introduces notch diversity to improve separating the unknown contributions of one or more sources from a commonly acquired set of wavefield signals while still allowing for optimal reconstruction properties in certain diamond-shaped regions. In particular, notch diversity is obtained by heteroscale encoding.

Systems and methods for generating seismic images of subterranean features including: receiving raw seismic data of a subterranean formation; selecting a portion of the raw seismic data; transforming the selected portion of the raw seismic data from a first domain to a second domain; generating soft constraint data corresponding to the selected portion of the raw seismic data; calculating at least one weight using the generated soft constraint data; generating a weighted transformed data set by applying at least one weight to the transformed selected portion of the raw seismic data; selecting at least one data point of the generated weighted transformed data set; and removing the selected at least one data point from the weighted transformed data set to generate revised seismic data.

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.

Disclosed are methods, systems, and computer-readable medium for a full seismic wavefield de-aliasing workflow. To achieve the de-aliasing, the workflow employs a four-dimension (4D) anti-leakage anti-aliasing regularization algorithm. The workflow involves application of successive de-aliasing steps while restricting computations only to the significant spatial dimensions. In areas of strong elastic property variation in the near-surface, the benefit of de-aliasing the full wavefield is both significant and demonstrable. In addition to achieving de-aliased sampling of the full wavefield, the workflow reduces the complexity of both the computational and geophysical aspects of the problem of de-aliasing full wavefields.

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.

A system for seismic imaging of a subterranean geological formation, the system includes a receiver configured to obtain seismic data comprising a data volume representing a post-stacked image. The system includes a filtering module configured to: apply frequency-wavenumber (F-K) filter to the data volume extract a negative-dip structure image and apply the frequency-wavenumber (F-K) filter to the data volume extract a positive-dip structure image. The system includes a diffraction rendering module configured to: multiply the positive-dip structure image with the negative-dip structure image and generate a diffraction-enhanced seismic image representing a geological formation of the data volume.

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.