A computer-implemented method and computing system apparatus programmed to perform operations of the computer-implemented method for obtaining a subsurface stack image, subsurface angle gathers, and a subsurface velocity model over an entire survey region having high velocity contrast geo-bodies. Particularly, user inputs, input velocity models, and surface-seismic data are obtained by fixed source and receiver pairs and then used by the computer program product embedded within the computing system apparatus to minimize the number of iterations, required to obtain a final velocity model, a final stack image, and final angle gathers wherein their flatness deviation is equal to, or less than, a user-defined flatness value. Therefore, the attributes developed by said computer-implemented method and system can help solve the imaging problem of sub high velocity contrast geo-bodies like subsalt, or salt overhung deep mini basins.

A system is described for determining a likelihood of a type of fluid in a subterranean reservoir. The system may include a processor and a non-transitory computer-readable medium that includes instructions executable by the processor to cause the processor to perform various operations. The processor may receive pre-stack seismic data having seismically-acquired data elements for geometric locations in a subterranean reservoir. The processor may determine, using the pre-stack seismic data, input features for each geometric location and may execute a trained model on the input features for determining a likelihood of a type of fluid in the subterranean reservoir and for determining a list of features affecting the likelihood. The processor may subsequently output the likelihood and the list of features.

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 system and method can be used for to calibrating time-lapse seismic volumes by cross-migration rescaling and reorientation for use in determining optimal wellbore placement or production in a subsurface environment. Certain aspects include methods for cross-migration of data sets processed using different migration techniques. Pre-processing of the data sets, optimization of rescaling and reorientation, and identification of adjustment parameters associated with minimum global error can be used to achieve a time-dependent formation data set that addresses error in all input data sets.

A computer-implemented method and computing system apparatus programmed to perform operations of the computer-implemented method for obtaining a subsurface stack image, subsurface angle gathers, and a subsurface velocity model over an entire survey region having high velocity contrast geo-bodies. Particularly, user inputs, input velocity models, and surface-seismic data are obtained by fixed source and receiver pairs and then used by the computer program product embedded within the computing system apparatus to minimize the number of iterations, required to obtain a final velocity model, a final stack image, and final angle gathers wherein their flatness deviation is equal to, or less than, a user-defined flatness value. Therefore, the attributes developed by said computer-implemented method and system can help solve the imaging problem of sub high velocity contrast geo-bodies like subsalt, or salt overhung deep mini basins.

A horizontal component of marine seismic survey data from an ocean bottom seismic survey can be migrated using a primary wave velocity model. The horizontal component can comprise a shear converted wave. An image of a subsurface location can based on the migration can be produced. Migrating the horizontal component can comprise wave-equation migrating the horizontal component, where the horizontal component is input as both a source wavefield and a receiver wavefield.

A method is described for seismic imaging including receiving a pre-stack seismic dataset and an earth model at one or more computer processors; performing least-squares reverse time migration of the pre-stack seismic dataset using the earth model to create a digital seismic image, wherein the least-squares reverse time migration includes wave-equation forward modeling based on an asymptotic expression for reflection in a subsurface Kirchhoff integral; and generating a display of the digital seismic image on a graphical user interface.

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.

Determination of an impulse response at a subsurface image level can include extrapolation of an up-going pressure wavefield to a subsurface image level, extrapolation of a down-going velocity wavefield to the subsurface image level, and determination of the impulse response at the subsurface image level from a hypothetical seismic source by spectral division of the extrapolated up-going pressure wavefield by the extrapolated down-going velocity wavefield.

A method of analysing seismic data to detect possible hydrocarbons includes determining a set of data tiles from a seismic data cube of seismic data and testing each data tile in the set of data tiles to determine whether it corresponds to a possible fluid contact.