Methods and systems, including computer programs encoded on a computer storage medium can be used for an integrated methodology that can be used by a computing system to automate processes for generating, and updating (e.g., in real-time), subsurface reservoir models. The methodology and automated approaches employ technologies relating to machine learning and artificial intelligence (AI) to process seismic data and information relating to seismic facies.

A method for modeling a subsurface volume includes receiving a plurality of ordered seismic images including representations of objects in the subsurface volume, generating flow fields based on a difference between individual images of the plurality of ordered seismic images, and identifying the objects in the seismic images based on the flow fields and the plurality of ordered seismic images.

Well information may define subsurface configuration of different wells. Marker information defining marker positions within the wells may be obtained. A dissimilarity matrix for the wells may generated, with the element values of the dissimilarity matrix determined based on comparison of corresponding subsurface configuration of the wells. A gated dissimilarity matrix may be generated from the dissimilarity matrix based on the marker positions within the wells. The elements values of the gated dissimilarity matrix corresponding to one set of marker positions and not corresponding to the other set of marker positions may be changed. Correlation between the wells may be determined based on the gated dissimilarity matrix such that correlation exists between a marker position in one well and a marker position in another well.

Disclosed herein are geologic modeling methods and systems employing function-based representations of horizons intersected by partial faults. An illustrative method embodiment includes: (a) obtaining a seismic image volume; (b) identifying a horizon within the seismic image volume, said horizon being intersected by a partial fault; (c) deriving a function-based representation of the horizon, the representation being continuous except across the partial fault; (e) constructing a watertight subsurface model using the function-based representation; (f) assigning petrophysical parameter values to compartments of the watertight subsurface model; and, optionally, (g) storing or displaying the watertight subsurface model.

Systems and methods develop a three-dimensional model of a subterranean formation based on vertical seismic profiles at a plurality of well locations. This approach can include receiving seismic data for the subterranean formation including the vertical seismic profiles; for each vertical seismic profile, injecting a ground force into the vertical seismic profile to provide a reference trace at depth zero to estimate energy loss in each receiver providing data in the vertical seismic profile and estimating time and depth variant quality factors for the well location associated with the vertical seismic profile based on the seismic profile; estimating quality factors for points within a three-dimensional volume representing the subterranean formation by interpolating between the time and depth variant quality factors for the location associated with each vertical seismic profile; and combining estimated quality factors to generate a three-dimensional quality factor model of the three-dimensional volume representing the subterranean formation.

A method for generating a geophysical image of a subsurface region includes defining a computational sub-volume for the geophysical image including a predetermined number of seismic traces of a plurality of seismic traces and a predetermined number of samples per each one of the plurality of seismic traces, generating a data matrix corresponding to a first sub-volume of the subsurface region based on the defined computational sub-volume, the data matrix comprising the predetermined number of samples for the predetermined number of traces of a portion of a seismic dataset corresponding to the first sub-volume. The method also includes performing a singular value decomposition and estimating a coherence between the predetermined number of traces of the data matrix by performing a weighted sum of a variance of a plurality of right singular vectors of the data matrix, and assigning the estimated coherence to a location in the geophysical image.

A computer executable algorithm adapted to propagate a boundary surface of a seed that is placed within a region of interest of a visual representation of a 3D seismic data so as to follow a natural contour of said region of interest, wherein said algorithm is executable to: (i) generate at least one attribute volume comprising at least on attribute derivable from said 3D seismic data set; (ii) generate at least one characteristic parameter for a plurality of candidate events of said 3D seismic data within a predefined gate region located forward of said propagating boundary surface; (iii) generate and assign a probability characteristic for said plurality of candidate events based on said at least one attribute volume and said at least one characteristic parameter; and propagate said boundary surface towards and incorporating any one of said plurality of candidate events that fulfils an acceptance criteria of said probability characteristic so as to generate a surface along the natural contour of said region of interest.

Geologic modeling methods and systems may use design-space to design-space mapping to facilitate simulation grid generation for multiple interpretations of a subsurface region. As one example, one or more embodiments of a geologic modeling method may comprise: obtaining first and second geologic models having different structural interpretations of a subsurface region; mapping each of the geologic models to associated design space models representing an unfaulted subsurface region; determining a design-to-design space mapping from the first design space model to the second design space model; using said mapping to copy parameter values from the first design space model to the second of the design space model; gridding each of the design space models to obtain design space meshes; partitioning cells in the first and second design space meshes along faults; reverse mapping the partitioned design space meshes to the physical space to obtain first and second physical space simulation meshes.

A system can include a processor; memory operatively coupled to the processor; and processor-executable instructions stored in the memory to instruct the system to: receive a marker on a well log for a well in a geographic region; and iteratively propagate the marker automatically to a plurality of well logs for other wells in the geographic region.

Systems and methods to correct misties across multiple 2D seismic surveys using a correction solution calculated based only on the intersecting points between different surveys.