A method includes receiving seismic data related to a subterranean formation. The method also includes receiving a selection of a property of the subterranean formation that is permitted to vary during a simulation of a model of the subterranean formation. The method also includes simulating fluid flow in the model of the subterranean formation based at least partially on the seismic data, and the selected property. The method also includes generating an updated model based at least partially upon a result of simulating the fluid flow in the model.

The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for models the accumulation and migration of hydrocarbons. One computer-implemented method includes: identifying one of grid cells neighboring an accumulation as a recent back-filled grid cell; setting an oil phase potential of identified grid cell as an accumulation potential of the accumulation; comparing oil phase potentials of grid cells neighboring the recent back-filled grid cell with the accumulation potential of the accumulation, where the oil phase potential of each of the grid cells neighboring the recent back-filled grid cell is stored as a key in the node corresponding to the respective grid cell; selecting one of the grid cells neighboring the accumulation as a next back-filled grid cell; and updating the accumulation potential of the accumulation based on the oil phase potential of the selected grid cell.

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.

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.

The invention discloses a dolomite reservoir prediction method and system based on well and seismic combination, and storage medium. The method steps include: obtaining the dolomite index characteristic curve through well log sensitivity analysis, and distinguishing the dolomite and limestone according to the difference in their response range; after the artificial intelligence deep learning is performed on the dolomite index characteristic curve of the drilling area, the dolomite index characteristic curve of the virtual drilling area is obtained; according to the dolomite index characteristic curve of the drilling area and the virtual drilling area, the post-stack seismic data is used for inversion to obtain the distribution and development status of the dolomite reservoir in the test area. The invention effectively distinguishes the dolomite and limestone through the dolomite index characteristic curve, and accurately predicts the distribution and development status of the dolomite reservoir in the test area with less wells.

An apparatus comprising a data acquisition tool including NMR sensors, a data acquisition processor communicatively coupled with the NMR sensors, and a first memory storing instructions that cause the data acquisition processor to perform operations comprising acquiring data of earth formation fluid, varying at least one of a magnetic field gradient and an inter-echo time, and acquiring additional data. The apparatus further comprises a data processing unit comprising a second memory storing instructions that cause the data processor to perform operations comprising receiving data acquired by the data acquisition tool, constructing a mathematical model of the data, conducting a first inversion of the mathematical model to obtain a first set of NMR responses, performing a forward model of the first set of NMR responses obtained from the first inversion, conducting a second inversion to obtain a second set of NMR responses, and determining earth formation fluid properties.

The present disclosure relates to a method for determination of a real subsoil geological formation. In at least one embodiment, the method includes receiving a model representing the real subsoil, determining a first fluvial geological formation in said model using parametric surfaces, determining a subsequent fluvial geological formation as a deformation of the first fluvial geological formation using parametric surfaces, and subtracting the first fluvial geological formation from the subsequent fluvial geological formation to create a new geological formation named point bar formation.

A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.

A parallel-processing hydrocarbon (HC) migration and accumulation methodology is applied to basin data to determine migration pathways and traps for high-resolution petroleum system modeling. HC is determined in parallel to have been expelled in source rocks associated with a plurality of grid cells divided into one or more subdomains. Potential trap peaks are identified within the plurality of grid cells. An invasion percolation (IP) process is performed until the HC stops migrating upon arrival to the plurality of trap peaks. A determination is made as to whether the grid cells containing HC contains an excess volume of HC. An accumulation process is performed to model the filling of the HC at a trap associated with the identified potential trap peaks. The trap boundary cell list is updated in parallel together with an HC potential value. Trap filling terminates when excess HC is depleted or a spill point is reached.