A least one seismic attribute is determined for each voxel of the seismic volume. A first horizon is selected for mapping and a sparse global grid is generated which includes the horizon, at least one constraint point identifying the horizon, and a number of points having a depth in the seismic volume. A value of at least one seismic attribute is determined for each point and their depths are adjusted based on the value of the seismic attribute. A map of the horizon can be generated based on the adjusted depths. Multiple local grids can be generated based on the sparse global grid, and the depths of the local grid points adjusted to generate a map of the horizon at voxel level resolution. The seismic volume can be mapped into multiple horizons, where previously mapped horizons can function as constraints on the sparse global grid.

The present invention proposes a method to represent seismic fault damage zones and fracture density in the geological models of reservoirs in a simple, agile and automated way, so that it can be easily replicated by geologists in any production design. It was developed as a group of workflows, inserted in the commercial software Petrel, widely used in the company for the 3D numerical modeling of reservoirs.

Methods, systems, and non-transitory computer readable media for identifying type curve regions as a function of position in a region of interest are disclosed. Exemplary implementations may include: obtaining a spatial clustering model from the non-transitory storage medium; obtaining well data from the non-transitory storage medium; obtaining production parameter data from the non-transitory storage medium; and delineating each of the type curve regions in the region of interest by applying the spatial clustering model to the well data and the production parameter data.

A method can include receiving data associated with a geologic environment; based on at least a portion of the data, estimating relationships for multiple properties of the geologic environment; and based at least in part on the relationships, performing simultaneous joint inversion for at least one property of the geologic environment.

System and methods for automated seismic processing are provided. Historical seismic project data associated with one or more historical seismic projects is obtained from a data store. The historical seismic project data is transformed into seismic workflow model data. At least one seismic workflow model is generated using the seismic workflow model data. Responsive to receiving seismic data for a new seismic project, an optimized workflow for processing the received seismic data is determined based on the at least one generated seismic workflow model. Geophysical parameters for processing the seismic data with the optimized workflow are selected. The seismic data for the new seismic project is processed using the optimized workflow and the selected geophysical parameters.

A backpropagation enabled model is trained for estimating time-lapse property changes of a subsurface volume. Synthetic models of the subsurface volume are generated, with pre-determined property changes before and after a time lapse. These models are used to compute baseline-monitor pairs of synthetic seismic traces before and after the time lapse, wherein the baseline synthetic traces are computed from the synthetic model before the time lapse and the monitor synthetic traces are computed from the synthetic model after the time lapse. A ground truth 4D attribute characterizing the time-lapse property changes in the synthetic models is defined, and a backpropagation enabled model is trained by feeding the baseline-monitor pairs of synthetic seismic traces and the corresponding ground truth 4D attribute. The thus obtained trained backpropagation enabled model can be used to estimate time-lapse property changes of the actual subsurface Earth volume from actual baseline-monitor pairs of seismic traces.

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 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.

A method and system are described for characterizing fractures in subsurface regions, using the fracture characterizations in subsurface models, and using the subsurface models in hydrocarbon operations. In the methods and systems, one or more zones are identified in a subsurface model for fracture characterization; a multi-layer model is created for the selected zone; and macromechanical geological loads are applied and simulated to generate fractures for the zone.

A method, including performing, with a computer, up/down separation of geophysical data, which produces an approximate up-going wavefield and an approximate down-going wavefield; creating an areal source based at least in part on the down-going wavefield; and performing, with a computer, a full wavefield inversion process with the areal source, and an objective function measuring a misfit between modeled up-going wavefields and recorded up-going wavefields, wherein the full wavefield inversion process generates a final subsurface physical property model.