A three-dimensional analog simulation test system for gas-liquid countercurrent in an abandoned mine goaf includes a three-dimensional analog simulation device configured to simulate an analog environment of gas-liquid countercurrent in an abandoned mine goaf, a gas supply system, an automatic water pressure control system and a hydraulic loading system. The three-dimensional analog simulation device includes a reaction frame and a three-dimensional analog simulation chamber. The reaction frame includes a beam, a base and a column for connecting the beam and the base. The three-dimensional analog simulation chamber is arranged inside the reaction frame. The test system can be used to simulate the gas-liquid countercurrent in an abandoned mine goaf so as to study the evolution of water accumulation and gas enrichment in long-term abandonment of closed mines and the migration evolution rules of gas-liquid two-phase in abandoned mines with high gas density.

Systems and methods for determining a numerical age for new geological events within a new scheme by ordering relations between geological events within a new scheme and/or within a new scheme and a preexisting scheme into a preferred hierarchy, dynamically excluding lower relations in the preferred hierarchy that conflict with higher relations due to irreconcilable ages of the relations, and using the ordered relations remaining in the preferred hierarchy to determine a numerical age for the new geological events within the new scheme.

The present invention relates to a method for determining hydrocarbon production for a reservoir. The method comprises determining a projector matrix based on a Jacobian matrix function of the gridded model, then splitting the Jacobian matrix into subsets of consecutive lines. For each subset of consecutive lines, creating a respective square matrix based on said subset. A determining eigenvectors and respective eigenvalues associated with the respective square matrix and then determining relevant eigenvectors having respective eigenvalues below a predetermined threshold. The projector is determined as a concatenation of the extended eigenvectors ordered according to multiple criteria: the respective order value of the subset; and the respective eigenvalue of the relevant eigenvector.

Method for mapping a 3D grid or mesh from a faulted subsurface domain to a continuous design domain, wherein the grid may be used to represent a discrete model of a subsurface material property (such as permeability) to use, for example, in a reservoir simulator. The mapping is geometry-based, not physics-based. The mapping is determined by an iterative optimization procedure designed to penalize deformation of tessellated mesh cells (703) in the design domain compared to their geometric quality in the faulted domain (701), but subject to stitching constraints (702) appearing as a penalty term or Lagrange multiplier term in the optimization objective function to influence the final mesh to co-locate pairs of points identified on opposite sides of a fault as having been located together before the fault occurred.

Techniques for evaluating hydrocarbon reserves using tool response models are provided. An example method includes measuring a first fluid distribution of a first formation proximate to a first wellbore and measuring a second fluid distribution of a second formation proximate to a second wellbore. The method further includes generating a first tool response model for the first formation based at least in part on the first fluid distribution and generating a second tool response model for the second formation based at least in part on the second fluid distribution. The method further includes comparing results of the first tool response model to results of the second tool response model to determine a fluid distribution difference between the first formation and the second formation and implementing a drilling command to alter drilling of one of the first and second wellbores based at least in part on the fluid distribution difference.

A method to measure fracture length and geometry/complexity from pressure decay and diffusion and near wellbore conductivity measurements with far-field conductivity estimates.

Simulation model cells of a reservoir having gas condensate production from layers of a subsurface hydrocarbon reservoir are determined by the Total Dynamic Productivity Index (TDPI) method. The determinations are based on rock and fluid physics, including non-Darcy flow and pseudo-pressure integral, which are responsible for pressure drop near wellbores in the gas condensate wells. A three-dimensional (3D) grid matrix is linked with a well placement optimization algorithm to target high gas energy spots and increase productivity, efficiency and recovery from the gas condensate wells.

The invention simulates flows in a geological reservoir having a heterogeneous pore size. From laboratory measurements on samples taken in the geological reservoir, pore size distribution classes are determined and a triple-porosity model representative of each class is determined. The flow simulator according to the invention implements the triple-porosity model, a thermodynamic equation of state accounting for an equivalent dimension of the pores of the small-size medium, fluid exchanges exclusively between the large-pore and small-pore media and between the small-pore and fracture media, and the capillary pressure as a function of the saturation in a small-pore medium.

Systems and methods for estimating a likelihood of an object element in a given position in a subsurface volume of interest are disclosed. Exemplary implementations may: obtain target subsurface data from the subsurface volume of interest; obtain an object element set corresponding to the subsurface volume of interest; generate correlation values as a function of position in the subsurface volume of interest by applying the object filters to the target subsurface data; and generate object element likelihood values by applying the object templates to positions in the subsurface volume of interest corresponding to the correlation values.

A system and method for modeling a geological structure may include, in an initial model, computing a first function for a geological structure including a first set of iso-surfaces. A processor may detect if the first set of iso-surfaces intersect a set of geological markers within a threshold proximity. If not, the initial model may be corrected using an induced mesh having an increased cell resolution compare to the initial model for computing a second function for the geological structure including a second set of iso-surfaces that intersect the geological markers within the threshold proximity. A processor may insert the second set of iso-surfaces into a second model to locally increase its resolution relative to the initial model by dividing cells in the second model along the second set of iso-surfaces. For each new geological structure, the above steps may be repeated using the second model as the initial model.