A method, system and non-transitory computer-readable medium for forming a seismic image of a geological structure are provided. After obtaining seismic wave data including a plurality of seismic wave traces at a first region of the geological structure, a predicted time dispersion error of an actual time dispersion error that results from a use of a finite difference approximation in calculating predicted seismic wave data at a second region of the geological structure as if a seismic wave propagates from the first region to the second region of the geological structure, is calculated. A corrected predicted seismic wave data at the second region of the geological structure is calculated by applying the finite difference approximation to the seismic wave data at the first region of the geological structure compensated with the predicted time dispersion error. A seismic image of the second region of the geological structure is generated using the corrected predicted seismic wave data, such that the actual time dispersion error is negated by the predicted time dispersion error.

Embodiments provide for a method that utilizes the azimuthally spaced receivers of a sonic logging tool. Signals from monopole and dipole sources are reflected from the geologic interfaces and recorded by arrays of receivers of the same tool. For the incident P-waves from the monopole source, phase arrival times for the azimuthal receivers are compensated for stacking using properties of wave propagation in the borehole, and for the incident SH-waves from the dipole source, signs of waveforms for the receivers are changed for specified azimuths.

Disclosed are methods of marine 3D seismic data acquisition that do not require compensation for winds and currents.

Implementations for monitoring seismic data recorded in a marine survey of a subterranean formation for coverage gaps are described herein. Implementations include computing Fresnel sum operators for Fresnel zones of the subterranean formation based on a Kirchhoff migration impulse response at horizons of a representative plane layer model of a survey area of the subterranean formation. Implementations also include computing an acceptability map of the survey area based on the Fresnel sum operators. The acceptability map reveals coverage gaps in the survey area. Geoscientist may use the acceptability map to infill seismic data in areas of the survey area that correspond to the coverage gaps.

Seismic data may provide valuable information with regard to the description such as the location and/or change of hydrocarbon deposits within a subsurface region of the Earth. The present disclosure generally discusses techniques that may be used by a computing system to interpolate or deblend data utilizing a projection on convex sets (POCS) interpolation algorithm. The utilized POCS interpolation algorithm operates in parallel for frequency of a set of frequencies of a seismic frequency spectrum.

Methods and systems, including computer programs encoded on a computer storage medium can be used for identifying primary-wave (P-wave) and secondary-wave (S-wave) characteristics of an underground formation by separating P-wave and S-wave modes of seismic data generated by applying a seismic source to a subterranean region of a geological area. Particle motion vectors of a P-wave are parallel to a propagation vector of the P-wave, whereas particle motion vectors of an S-wave are perpendicular to a propagation vector of the S-wave. The parallel and perpendicular relationship between the motion and propagation vectors of the respective P- and S-waves provide a basis for separating P- and S-wave components from a wavefield. The separation methodology extracts P-wave components and S-wave components from the wavefield based on an estimated angle between propagation vectors and wave motion vectors for the wavefield.

A method can include receiving seismic survey data of a subsurface environment from a seismic survey that includes a source arrangement of sources that is spatially denser than a receiver arrangement of receivers; processing the seismic survey data using the principle of reciprocity for performing interpolation across the receivers to generate processed seismic survey data; and generating an image of at least a portion of the subsurface environment using the processed seismic survey data.

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.

Recorded seismic data includes seismic traces having respective source orientation angles, where the source orientation angles represent deviations in seismic source orientation relative to an inline survey direction. A plurality of designature operators corresponding to respective designature orientation angles within a defined set of designature orientation angles may be generated. For a given member of the defined set of designature orientation angles, a corresponding designature operator may be applied to the recorded seismic data to generate designatured seismic data for the given designature orientation angle. For a given seismic trace having a given source orientation angle, the designatured seismic data may be interpolated to generate a designatured version of the given seismic trace. The results may be stored in a tangible, computer-readable medium.

Computational systems and methods for interpolating a pressure wavefield based on pressure wavefield and particle motion wavefield data are disclosed. The pressure and particle motion wavefields are sampled at a sample rate that is less than a lower bound for sample rates typically used to interpolate the pressure wavefield from the pressure wavefield samples alone. The particle motion wavefield can an acceleration wavefield or a time derivative of a velocity wavefield.