Methods and systems for forming a three-dimensional (“3D”) seismic image of a subterranean region of interest is disclosed. The method includes obtaining a seismic dataset a seismic trace for each of a plurality of pairs of one source and one receiver location and obtaining a 3D travel-time cube for each source location and each receiver location. The method further includes dividing the seismic dataset into a plurality of seismic subsets composed of set of source locations, set of receiver locations a seismic trace for each pair of source and receiver location and the 3D travel-time cube for each source for each receiver location. The method still further includes transmitting, to a random-access memory block of a computer processing unit the seismic subset, and forming a seismic partial image based on the seismic subset, and determining the 3D seismic image based on a combination of the seismic partial images.

Disclosed are methods, systems, and computer-readable medium to perform operations including: receiving an input dataset that represents partial spatial information of an area of interest; providing the input dataset to a spatial context generator, wherein the spatial context generator comprises a machine learning model trained to generate, based on the partial spatial information, contextual spatial information for the area of interest; and using the spatial context generator to generate, based on the partial spatial information, at least one output dataset associated with the area of interest, where each output dataset comprises simulated contextual spatial information for the area of interest.

Hybrid seismic inversion methods and apparatuses perform wave equation inversion and stochastic inversion to generate one or more final models for the reservoir characterization of the survey region. A method may include retrieving seismic data using seismic data recording sensors; storing the seismic data in the database; retrieving well data using the well born sensor in the wellbore; storing the seismic data in the database; storing geology integration information and one or more background models in the database; retrieving the seismic data and processing the seismic data to mitigate the seismic data for a seismic hybrid inversion; and performing the seismic hybrid inversion including performing wave equation inversion and stochastic inversion to generate the one or more final models for the reservoir characterization of the survey region.

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.

A method may include obtaining seismic data for a geological region of interest. The method may further include determining various velocity semblance values for the geological region of interest using a time window and the seismic data. The method may further include determining, automatically by a computer processor, one or more stacking velocities for the geological region of interest using a traced path based on the velocity semblance values and a path tracing algorithm. The path tracing algorithm may recursively determine an accumulated amplitude array based on the velocity semblance values. The path tracing algorithm may further determine the traced path among the velocity semblance values and the accumulated amplitude array, the traced path corresponding to the one or more stacking velocities. The method may further include generating a velocity model of the geological region of interest using the one or more stacking velocities.

A fracture model for a hydraulic fracture in a wellbore can be updated and calibrated. Information about a microseismic event can be received from a sensor that is monitoring a subterranean formation. The information can be received subsequent to a fracking fluid being introduced into the formation. An observed geometry of a hydraulic fracture can be determined based on the information and a predicted geometry of the fracture can be determined based on properties of the fracking fluid and a fracture model. The fracture model can be updated using the information about the microseismic event where it is determined that an uncertainty value of the observed geometry does not exceed a pre-set maximum. The uncertainty value can be based on the predicted geometry of the hydraulic fracture.

The disclosed technology provides solutions for identifying noise in seismic profile data sets. In some aspects, a process of the disclosed technology includes steps for receiving wellbore data including seismic measurements, processing the wellbore data to generate a seismic input image including visual representations of the one or more seismic measurements, and processing the seismic input image to identify a noise region in the seismic input image. Systems and machine-readable media are also provided.

In one embodiment, a reverse time migration module is configured with a visco-acoustic wave equation for media with heterogeneous attenuation solved using a pseudo-analytical method. Seismic data is obtained for a zone of interest, and a model is created for the zone of interest. The model has spatial variability in velocity and quality factor. Pseudo-analytic Q-compensating reverse time migration (PA-Q-RTM) is performed using the reverse time migration module and the model for the zone of interest to obtain PA-Q-RTM seismic data. In another embodiment, a system for processing seismic data includes a reverse time migration module configured with a visco-acoustic wave equation that is solved from the pseudo-analytical method. The system operates to obtain seismic data for a zone of interest and perform PA-Q-RTM using a model for the zone of interest to obtain PA-Q-RTM seismic data for the zone of interest. Other embodiments and features are also disclosed.

Disclosed are systems and methods for numerically simulating seismic-wave propagation in tilted transversely isotropic (TTI) media, using an adaptive Lebedev staggered grid. In various embodiments, the adaptive grid includes multiple horizontal zones having different associated grid spacings, which may be determined based on a vertical wave-velocity model. The numerical simulation may involve iteratively solving a set of finite-difference equations including finite-difference coefficients that vary spatially depending on the grid spacing. Additional embodiments and features are described.

A global objective function is initialized to an initial value. A particular model simulation process is executed using prepared input data. A mismatch value is computed by using a local function to compare an output of the particular model simulation process to corresponding input data for the particular model simulation process. Model objects associated with the particular model simulation process are sent to another model simulation process. An optimization process is executed to predict new values for input data to reduce the computed mismatch value.