Methods for updating a physical properties model of a subsurface region that combine advantages of FWI and tomography into a joint inversion scheme are provided. One method comprises obtaining measured seismic data; generating simulated seismic data using an initial model; computing at least one of an FWI gradient and a tomography gradient; minimizing a joint objective function E, wherein the objective function E is based on a combination of the FWI gradient and the tomography gradient or preconditioning of the FWI gradient or the tomography gradient; generating a final model based on the minimized joint objective function E; and using the final model to generate a subsurface image. The joint objective function E may be a function of one or both a FWI objective function C.sub.FWI and a tomography objective function C.sub.Tomo. The joint objective function E may be defined as a weighted sum of C.sub.FWI and C.sub.Tomo.

A method and apparatus for generating a high-resolution seismic image, including extracting a reflectivity distribution from a geological model; utilizing the reflectivity distribution to label features of the model; generating forward-modeled data from the model; migrating the forward-modeled data to create a migrated image; and training a deep neural network with the labeled synthetic geological model and the migrated image to create a reflectivity prediction network. A method and apparatus includes: selecting a first subset of the field data; applying a low-pass filter to the first subset to generate a first filtered dataset; migrating the first filtered dataset to create a first migrated image; applying a high-pass filter to the first subset to generate a second filtered dataset; migrating the second filtered dataset to create a second migrated image; and training a deep neural network to predict a target distribution of high-frequency signal.

Systems and methods for seismic imaging of a subterranean geological formation include receiving parameter data representing one or more parameters of a seismic survey, the seismic data specifying an incident angle and an azimuth angle for each trace of the seismic survey; determining a relationship between the incident angle and the azimuth angle for each trace and a location in a spiral coordinate system, and generating a weighting function for applying a weight value to each trace seismic data based on the incident angle and the azimuth angle associated with each trace; and determining a residual moveout value of the seismic data for each location in the spiral coordinate system by applying the weighting function to each; and generating a seismic image representing the residual moveout value of the seismic data for each location in the spiral coordinate system.

A method and system for locating a reflector in a formation. The method may comprise broadcasting a sonic waveform as a shear formation body wave or a compressional formation body wave into the formation, recording a reflected wave from a reflector with the one or more receivers as dipole data by the dipole receiver and quadrupole data by the quadrupole receiver, and processing the dipole data and the quadrupole data with an information handling system to determine a location of the reflector from the borehole sonic logging tool. The system may comprise a borehole sonic logging tool and an information handling system. The borehole sonic logging tool may comprise one or more transmitters configured to transmit a sonic waveform into a formation and one or more receivers configured to record a reflected wave as a dipole receiver for dipole data and a quadrupole receiver for quadrupole data.

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.

Methods and apparatuses for seismic data processing perform a least-squares reverse time migration method in which surface-attribute-independent coefficients for the surface attribute gathers are demigrated to reduce the computational cost.

Methods and systems for adjusting migration gather are disclosed. A method for adjusting migration gather may include generating a first migration gather based on acquired seismic data. The method may also include generating a compensation volume based on reference data produced utilizing a reference seismic model. The method may further include applying the compensation volume to adjust the first migration gather.

A system and method for reducing migration distortions in migrated images of the Earth's subsurface. Recorded seismic data may be migrated, using a migration velocity model, to generate a migration image comprising distortions. Synthetic seismic data may be generated, using the migration velocity model, for a grid of scattered points. The synthetic seismic data may be migrated, using the migration velocity model, to generate impulse responses for the scattered points. The impulse responses are used as point spread functions (PSFs) which approximates the blurring operator, e.g., the Hessian operator. An optimal reflectivity model may be selected using image-domain least-squares migration (LSM), based on the PSFs, with a regularization of the difference between the migration image and a reflectivity model and a total variation (TV) regularization of the reflectivity model. An image of the optimal reflectivity model may be generated that has reduced migration distortions compared to the original migration image.

A method for applying separated wave-field imaging onshore (1) by artificially creating up-going and down-going fields and (2) by using these fields in a migration algorithm. If there are any surface multiples in the data, the resulting image created using the migration algorithm will be distorted by the unknown free-surface reflection coefficient. In fact, the surface multiples may be generated with a complex series of reflection coefficients. The distortions found in the resulting image created using the migration algorithm are then removed.

A method for imaging one or more dipping structures is provided. The method comprises providing a given velocity model, calculating travel time of a seismic wave using the velocity model, estimating coherency of event signals of array data that are time-corrected for a trial reflector, weighting a waveform sample corresponding to the travel time based on the coherency, and mapping the weighted waveform sample.