A system provides seismic images of the subsurface by enhancing pre-stack seismic data. The system obtains seismic data comprising a plurality of seismic traces that are generated by measuring reflections of seismic waves emitted into a geological formation. The system sorts seismic data into at least one multidimensional gather comprising a data domain. The system determines local kinematical attributes of a seismic trace. The system forms an ensemble of seismic traces, each representing a reference point. The system applies local moveout corrections to each seismic trace of the ensemble. The system applies residual statics and phase corrections for each seismic trace that is corrected by the local moveout corrections. The system sums the seismic traces of the ensemble to obtain an output seismic trace having an increased signal-to-noise ratio (SNR) relative to the seismic trace that represents the reference point for the ensemble of seismic traces.

Techniques for determining a drill bit location includes identifying a plurality of acoustic energy signals received at a plurality of sets of acoustic receivers from a passive acoustic energy source that is part of a wellbore drilling system; processing the plurality of acoustic energy signals; determining a location of a drill bit of the wellbore drilling system based on the processed plurality of acoustic signals; and updating a geo-steering path of the drill bit based on the determined location of the drill bit.

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.

This specification describes workflows for, but is not limited to, performing full waveform inversion (FWI) to build high resolution velocity models to improve the accuracy of seismic imaging of a subterranean formation. This specification describes processes to automatically edit and enhance S/N quality of seismic data (such as land seismic data) to prepare the datasets for FWI. The methods for automatic corrections and pre-processing include: automatic iterative surface-consistent residual statics calculation, automatic rejection of anomalous traces (such as dead traces), and the automatic correction of surface-consistent amplitude anomalies (such as by scalar or deconvolution approaches). The operations include automatic “muting” of noise before first arrivals.

One method interpolates simulated seismic data of a coarse spatial sampling to a finer spatial sampling using a neural network. The neural network is previously trained using a set of simulated seismic data with the finer spatial sampling and a subset thereof with the coarse spatial sampling. The data is simulated using an image of the explored underground formation generated using real seismic data. The seismic dataset resulting from simulation and interpolation is used for denoising the seismic data acquired over the underground formation. Another method demigrates seismic data at a sparse density and then increases density by interpolating traces using a neural network.

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.

Methods, systems, and computer-readable medium to perform operations including: generating a first time-frequency spectrum of a first seismic trace from an original seismic dataset; generating a second time-frequency spectrum of a second seismic trace from an enhanced seismic dataset, where the second seismic trace corresponds to the first seismic trace; and re-combining an amplitude spectrum of the first time-frequency spectrum and a phase spectrum of the second time-frequency spectrum to generate a third time-frequency spectrum of an output trace that corresponds to the first and second seismic traces.

An apparatus, method, and system for determining body wave slowness from guided borehole waves. The method includes selecting a target axial resolution based on the size of a receiver array, obtaining a plurality of waveform data sets corresponding to a target formation zone and each acquired at a different shot position, computing a slowness-frequency 2D dispersion semblance map for each waveform data set, stacking the slowness-frequency 2D dispersion semblance maps to generate a stacked 2D semblance map, and determining a body wave slowness from the extracted dispersion curve. The method may also include generating a self-adaptive weighting function based on a dispersion model and the extracted dispersion curve, fitting the weighted dispersion curve and the dispersion model to determine a body wave slowness that minimizes the misfit between the weighted dispersion curve and the dispersion model. The method can be applied to both frequency-domain and time-domain processing.

Methods for full waveform inversion (FWI) in the midpoint-offset domain include using a computer system to sort seismic traces into common midpoint-offset bins (XYO bins). For each XYO bin, a linear moveout correction is applied to a collection of seismic traces within the XYO bin. The collection of seismic traces is stacked to form a pilot trace. The computer system determines a surface-consistent residual static correction for each seismic trace. The computer system determines that the surface-consistent residual static correction for each seismic trace is less than a threshold. Responsive to the determining that the surface-consistent residual static correction is less than the threshold, the computer system stacks the collection of seismic traces to provide the pilot trace. The computer system groups the pilot traces for the XYO bins into a set of virtual shot gathers. The computer system performs one-dimensional FWI based on each virtual shot gather.

A method of enhancing seismic images includes receiving a seismic gather. The seismic gather includes a plurality of seismic traces. A feature trace is generated based on the plurality of seismic traces in the seismic gather. For each of the plurality of seismic traces in the seismic gather, a correlation trace is generated based on that seismic trace and the feature trace, the correlation trace is modified using an activation function, and an enhanced trace is generated by multiplying that seismic trace with the modified correlation trace.