Saudi Arabian Oil Company

Trofimuk Institute of Petroleum Geology and Geophysics

Trofikmuk Institute of Petroleum Geology and Geophysics

Dmitry Neklyudov

Andrey Bakulin

Ilya Silvestrov

Maxim Dmitriev

Maxim Protasov

Kirill Gadylshin

Mikhail Belonosov

Vladimir Tcheverda

Vladimir Cheverda

Disclosed are methods, systems, and computer-readable medium to perform operations including: receiving prestack single sensor seismic data; representing traveltime moveout of the prestack single sensor seismic data locally as a second-order curve; calculating, using the prestack single sensor seismic data, local kinematic parameters that define the second-order curve; and performing, based in part on the local kinematic parameters, wavefield transformation on the single sensor seismic data to generate enhanced prestack single sensor seismic data.

Enhancement of single sensor seismic data

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.

Enhancement of seismic data

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; calculating a difference between the first time-frequency spectrum and the second time-frequency spectrum to generate a noise estimate in the first seismic trace; constructing, based on (i) the noise estimate, (ii) the first time-frequency spectrum, and (iii) the second time-frequency spectrum, a time-frequency mask (TFM); and using the constructed TFM to generate a third time-frequency spectrum of an output trace that corresponds to the first and second seismic traces.

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; calculating a difference between the first time-frequency spectrum and the second time-frequency spectrum to generate a noise estimate in the first seismic trace; characterizing the initial noise estimate as White Gaussian Noise (WGN); calculating, based on the characterization of the initial noise estimate, a third time-frequency spectrum of a refined noise estimate; constructing, based on the first time-frequency spectrum, the second time-frequency spectrum, and the third time-frequency spectrum, a time-frequency mask (TFM); and using the constructed TFM to generate a fourth time-frequency spectrum of an output trace that corresponds to the first and second seismic traces.

Methods, systems, and apparatus, including computer programs encoded on a computer storage medium, for determining a velocity model for a geological region. In one aspect, a method comprises: obtaining a current velocity model for the geological region; obtaining pre-stack and post-stack seismic data characterizing the geological region; and for each of a plurality of iterations: identifying a plurality of reflection events from the post-stack seismic data and the current velocity model; determining a respective observed travel time for each of the plurality of reflection events, comprising, for each reflection event, determining the respective observed travel time for the reflection event based at least in part on kinematic features derived from a respective seismic trace included in the pre-stack seismic data; and updating the current velocity model based at least in part on the observed travel times of the plurality of reflection events.

Seismic velocity modeling

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.

Providing seismic images of the subsurface using enhancement of pre-stack seismic data

The disclosure provides systems and methods to enhance pre-stack data for seismic data analysis by: sorting the reflection seismic data acquired from cross-spread gathers into sets of data sections; performing data enhancement on the sets of data sections to generate enhanced traces by: (i) applying forward normal-moveout (NMO) corrections such that arrival times of primary reflection events become more flat, (ii) estimating beamforming parameters including a nonlinear traveltime surface and a summation aperture, (iii) generating enhanced traces that combine contributions from original traces in the sets of data sections, and (iv) applying inverse NMO corrections to the enhanced traces such that temporal rearrangements due to the forward NMO corrections are undone.

Systems and methods to enhance 3-D prestack seismic data based on non-linear beamforming in the cross-spread domain

Certain implementations of a three-dimensional elastic frequency domain iterative solver for full waveform inversion can be implemented as a method in which frequency domain numerical simulation of elastic waves is performed in three-dimensional (3D) media.

Growing community of inventors

Novosibirsk, Russia

Dmitry Neklyudov