A surface-consistent refraction analysis method to automatically derive near surface corrections for seismic data processing. The method uses concepts from surface-consistent analysis applied to refracted arrivals. The method includes the use of CMP-offset-azimuth binning, evaluation of mean travel time and standard deviation for each bin, rejection of anomalous first break (FB) picks, derivation of CMP-based travel time-offset functions, conversion to velocity-depth functions, evaluation of long wavelength statics and calculation of surface-consistent residual statics through waveform cross-correlation. Residual time lags are evaluated in multiple CMP-offset-azimuth bins by similarity analysis with a pilot trace for all the other traces in the gather where the correlation window is centered at the refracted arrival. The similarity analysis may take the form of computerized cross-correlation, or other criteria such as semblance. The residuals are then used to build a system of linear equations that is simultaneously inverted for surface-consistent shot and receiver time shift corrections plus a possible subsurface residual term. All the steps are completely automated and require a fraction of the time needed for conventional near surface analysis.

A surface-consistent refraction analysis method to automatically derive near surface corrections for seismic data processing. The method uses concepts from surface-consistent analysis applied to refracted arrivals. The method includes the use of CMP-offset-azimuth binning, evaluation of mean travel time and standard deviation for each bin, rejection of anomalous first break (FB) picks, derivation of CMP-based travel time-offset functions, conversion to velocity-depth functions, evaluation of long wavelength statics and calculation of surface-consistent residual statics through waveform cross-correlation. Residual time lags are evaluated in multiple CMP-offset-azimuth bins by similarity analysis with a pilot trace for all the other traces in the gather where the correlation window is centered at the refracted arrival. The similarity analysis may take the form of computerized cross-correlation, or other criteria such as semblance. The residuals are then used to build a system of linear equations that is simultaneously inverted for surface-consistent shot and receiver time shift corrections plus a possible subsurface residual term. All the steps are completely automated and require a fraction of the time needed for conventional near surface analysis.

An example system for noise removal in distributed acoustic sensing data may include a distributed acoustic sensing (DAS) data collection system and an information handling system coupled thereto. The information handling system may receive seismic information from the DAS data collection system. The seismic information may include seismic traces associated with a plurality of depths in the wellbore. The information handling system may also generate a noise pilot trace by stacking one or more of the seismic traces, and subtract the noise pilot trace from the seismic information received from the DAS data collection system.

A system may obtain orthogonal pairs of acoustic waveforms obtained by a rotating downhole acoustic tool. The downhole tool may include an acoustic transmitter and an acoustic receiver. A computing device communicatively coupled to the downhole tool that instructs the transmitter to generate the various waveforms detected by the acoustic receiver and store the various waveforms in a buffer. Each waveform of a subset of the various waveforms oriented in a similar direction may be stacked to generate a composite waveform with reduced noise. The computing system may determine an orthogonal pair of waveforms based on the composite waveform and remaining waveforms.

Methods, apparatuses, and systems are disclosed that reduce seismic noise. In one embodiment, a method of processing seismic data includes accessing seismic data representative of a plurality of seismic input traces acquired by one or more seismic sensors. The method also includes stacking the plurality of seismic input traces into a stacked trace. The method also includes generating, utilizing at least one processor unit, a function of similarity between at least two of the plurality of seismic input traces. The method also includes scaling at least one of the seismic input traces or the stacked trace with the function of similarity.

A method may include obtaining various seismic traces for a geological region of interest. The method may further include determining an offset attribute and an azimuthal attribute. The method may further include determining, using the offset attribute and the azimuthal attribute, a virtual trace bin for the geological region of interest. The method may further include generating a virtual trace using a subset of the seismic traces and corresponding to the virtual trace bin. The method may further include generating a velocity model for the geological region of interest using a virtual shot gather including the virtual trace and various virtual traces. A respective virtual trace among the virtual traces may correspond to a respective virtual trace bin among various virtual trace bins. The method may further include generating a seismic image of the geological region of interest using the velocity model.

A method of performing notch compensation and a system to perform notch compensation for a first seismic streamer are described. The method includes disposing the first seismic streamer at a first depth, where the seismic streamer includes a first set of sensors to receive reflections resulting from a seismic source, the reflections indicating a notch at a frequency. The method also includes disposing a second seismic streamer at a second depth, the second depth being less than the first depth and the second seismic streamer including a second set of sensors to receive reflections resulting from the seismic source. The method further includes processing the reflections received by the first set of sensors and the second set of sensors together to derive the match filter, and applying the match filter to the reflections received by the first set of sensors of the first seismic streamer to compensate for the notch.

A method of performing notch compensation and a system to perform notch compensation for a first seismic streamer are described. The method includes disposing the first seismic streamer at a first depth, where the seismic streamer includes a first set of sensors to receive reflections resulting from a seismic source, the reflections indicating a notch at a frequency. The method also includes disposing a second seismic streamer at a second depth, the second depth being less than the first depth and the second seismic streamer including a second set of sensors to receive reflections resulting from the seismic source. The method further includes processing the reflections received by the first set of sensors and the second set of sensors together to derive the match filter, and applying the match filter to the reflections received by the first set of sensors of the first seismic streamer to compensate for the notch.

A computer implemented method for denoising a set of seismic datasets, specifically belonging to different 3D subsets of a 4D survey the method including: (a) receiving a baseline and a monitor seismic dataset which were acquired by surveying over the same subsurface formation over different periods of time; (b) cross-equalizing the monitor seismic dataset to match to the baseline seismic dataset in terms of amplitude, frequency, phase and timing of events; (c) computing an initial 4D difference between the monitor and baseline seismic datasets; (d) formulating a common noise template from the initial 4D difference; (e) de-noising the baseline and monitor seismic datasets, independently, using the common noise template in a curvelet domain; (f) updating the initial 4D difference to form an updated 4D difference, which reflects de-noised baseline and monitor datasets from step (e); and iterating the steps (d) through (F) until the updated 4D difference satisfies a predetermined criteria.

A method and apparatus of locating subsurface geohazards in a geographical area that includes: receiving a plurality of seismic trace signals in the geographical area based on a shot gather from a seismic shot source; isolating and stacking the plurality of seismic trace signals to generate a windowed trace signal associated with refraction traces from the seismic shot source; transforming the windowed trace signal to a frequency domain; calculating a low frequency to high frequency ratio for the transformed trace signal; outputting the calculated ratio to a two-dimensional array representing the geographical area at a source location and at a mean receiver location; repeating the steps for a plurality of other shot gathers in the geographical area; and multiplying each source location ratio with one or more mean receiver location ratios on the two-dimensional array to generate a final frequency ratio map.