Surface wave coda in seismic data recorded with a data acquisition system over an underground formation is estimated using time-reversal experiments. First time-reversal experiments use a first time-reversal mirror including a target source and one or more other sources to obtain estimates of surface waves traveling from other receivers to a target receiver. Second time-reversal experiments obtain a coda estimate for a surface wave traveling from the target source to the target receiver using a second time-reversal mirror including the target receiver and the other receivers.

A method includes receiving information for a subsurface region; based at least in part on the information, identifying sub-regions within the subsurface region; assigning individual identified sub-regions a dimensionality of a plurality of different dimensionalities that correspond to a plurality of different models; via a model-based computational framework, generating at least one result for at least one of the individual identified sub-regions based at least in part on at least one assigned dimensionality; and consolidating the at least one result for multiple sub-regions.

A data acquisition program, which includes core, image log, microseismic, DAS, DTS, and pressure data, is described. This program can be used in conjunction with a variety of techniques to accurately monitor and conduct well stimulation.

A method and apparatus for imaging seismic data includes obtaining an initial model of a subsurface formation, wherein the model includes a plurality of nodes that form at least part of a grid; an initial dip value for the nodes; and a set of origin coordinates for each of the nodes; performing bottom-up ray tracing for each node in the model, resulting in a set of arrival coordinates for each node; identifying a plurality of gathers from the seismic data; for each gather: calculating a set of midpoint coordinates; defining a midpoint vicinity surrounding the set of midpoint coordinates; identifying the nodes having arrival coordinates within the midpoint vicinity; and estimating a unique aperture for each of the gathers based on the respective origin coordinates; storing the estimated apertures in a table; and generating a subsurface volume or image with subsurface reflectors determined with apertures of the respective gathers.

Apparatus, computer instructions and method for deghosting seismic data related to a subsurface of a body of water. The method includes inputting data recorded by detectors that are towed by a vessel, the data being associated with waves travelling from the subsurface to the detectors; applying a migration procedure to the data to determine a first image of the subsurface; applying a mirror migration procedure to the data to determine a second image of the subsurface; joint deconvoluting the first image and the second image for deghosting a reflectivity of the subsurface; and generating a final image of the subsurface based on the deghosted reflectivity of the joint deconvoluting step.

A computer-implemented method for determining a subsurface feature, includes: determining a first velocity model based on an initial velocity model; generating a second velocity model based on measured seismic data at one or more first frequencies, the first velocity model, and a full waveform inversion (FWI); and in response to the second velocity model not satisfying a preset condition, performing a seismic forward simulation on the second velocity model to generate simulated seismic data at one or more second frequencies lower than the one or more first frequencies; updating the first velocity model based on the simulated seismic data at the one or more second frequencies; and updating the second velocity model based on the measured seismic data at the one or more first frequencies, the updated first velocity model, and the FWI, to determine the subsurface feature.

Method for characterizing subterranean formation is described. One method involves simulating a poroelastic pressure response of known fracture geometry utilizing a geomechanical model to generate a simulated poroelastic pressure response. Compiling a database of simulated poroelastic pressure responses. Measuring a poroelastic pressure response of the subterranean formation during a hydraulic fracturing operation to generate a measured poroelastic pressure response. Identifying a closest simulated poroelastic pressure response in the library of simulated poroelastic pressure response. Estimating a geometrical parameter of a fracture or fractures in the subterranean formation based on the closest simulated poroelastic pressure response.

A method may include obtaining a P-wave velocity model and velocity ratio data regarding a geological region of interest. The method may further include generating, based on the P-wave velocity model and the velocity ratio data, an initial S-wave velocity model regarding the geological region of interest. The method may further include determining various velocity boundaries within the initial S-wave velocity model using a trained model. The method may further include updating the initial S-wave velocity model using the velocity boundaries, an automatically-selected cross-correlation lag value based on various seismic migration gathers, and a migration-velocity analysis to produce an updated S-wave velocity model. The method further includes generating a combined velocity model for the geological region of interest using the updated S-wave velocity model and the P-wave velocity model.

A method of estimating lateral positioning uncertainties of a seismic image is performed at a computer system. The computer system receives a velocity model, the velocity model including a plurality of base velocity values used for generating the seismic image, each base velocity value having a low limit and a high limit. The computer system derives a plurality of lateral velocity gradient uncertainties from the velocity model and generates multiple lateral velocity gradient profiles, each lateral velocity gradient profile including a random realization of the plurality of lateral velocity gradient uncertainties. The computer system calculates perturbation raypaths originating from a surface point of the seismic image based on the velocity model and the multiple lateral velocity gradient profiles and estimates a lateral positioning uncertainty for a target location at a predefined depth of the seismic image based on lateral distributions of the perturbation raypaths at the predefined depth.

A method includes the steps of receiving a wavefield generated by reflections in a subsurface region and recorded by a plurality of seismic receivers and compensating the recorded wavefield for amplitude attenuation. The method further includes modelling a propagation of a source wavefield forward in time, from an initial time-state to a final time-state through an earth model that is representative of the subsurface region, wherein the modelling includes phase and amplitude effects of attenuation and modelling a propagation of the compensated recorded wavefield backward in time from a final time-state to an earlier time-state through the earth model, wherein the subsurface region has an absorption characteristic that dampens the recorded wavefield wherein the modelling includes phase and amplitude effects of attenuation.