A method for improved prediction and enhancement of hydrocarbon recovery from ultra-tight/unconventional reservoirs for both the primary production and any subsequent solvent huff‘n’puff periods based on facilitating the diffusion process may include steps of defining one or more initial properties of a reservoir and integrating characterization data of the reservoir; defining a wellbore trajectory for each of at least one well and one or more parameters associated with a completion/reservoir stimulation design; specifying operating conditions for a current development cycle; performing diffusion-based dynamic fracture/reservoir simulation for calculating hydrocarbon recovery and efficiency of a hydrocarbon process; and; determining whether to commence or continue enhanced oil recovery (EOR) or enhanced gas recovery (EGR) cycles.

Systems and methods are disclosed for generating subsurface feature prediction probability distributions from a subsurface feature as a function of position in a subsurface volume of interest. For example, a computer-implemented method may include: obtaining subsurface data and well data, generating subsurface feature values, generating subsurface feature realizations, generating subsurface feature realization uncertainty values, generating subsurface parameter values, generating subsurface parameter realizations, generating subsurface feature prediction probability distributions, generating a first representation of likelihoods of the subsurface feature values, and displaying the first representation.

The invention notably relates to a computer-implemented method of geomodelling. The method comprises providing a pseudo-stratigraphic grid. The pseudo-stratigraphic grid represents a reservoir and has pillars. Each pillar includes respective cells. Each cell has a respective stratigraphic layering index. The method then comprises providing a surface. The surface has a first region and a second region. The second region is complementary to the first region. The method also comprises, for each first pillar intercepted by the first region, determining a respective first stratigraphic layering value based on the relative position of the surface in the first pillar. The method also comprises, for each second pillar intercepted by the second region, determining a respective second stratigraphic layering value by interpolating and/or extrapolating first stratigraphic layering values. This provides an improved solution of geomodeling.

A DAS VSP technique is used to determine the induced fracture height and fracture density of an induced fracture region. The DAS VSP technique obtains pre-hydraulic fracturing DAS VSP survey time-lapse data to establish a baseline reference for the direct acoustic wave travel time. The DAS VSP technique obtains one or more time-lapse data corresponding to the subsequent monitor surveys conducted after each hydraulic fracturing stage along the well. Forward modeling is used to determine a theoretical acoustic wave travel time difference. The forward modeling uses seismic anisotropy to describe the behavior of seismic waves traveling through the induced fracture regions. An inversion scheme is then used to invert for the induced fracture height and the fracture density using the forward modeling. The two extracted induced fracture characteristics may then be used to determine optimal hydraulic fracturing parameters.

Hybrid 3D geocellular grids are generated to represent a subset of a natural fracture network (“NFN”) directly in the simulation, while the remainder of the NFN is approximated by a multi-continuum formulation. The resulting output is a 3D geocellular grid that possesses a higher level of mesh resolution in those areas surrounding the first fracture subsets, and lower mesh resolution in the areas of the second fracture subset.

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.

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.

A method of predicting hydraulic fracturing and associated risks of hydraulic fracturing operation is proposed. The methods use the mathematical simulation which allow to predict the geometry of a hydraulic fracture and location of fluids, propping agents (proppant), fibers and other materials therein. Reconsidering the fracturing design allows to remedy the possible risks (overflush, screen-out, bridging, gel contamination, temperature effects).

Systems and methods of the present disclosure are directed to reservoir simulation modeling using upon rock compaction tables derived from physical pore compressibility tests. The illustrative methods transform rock mechanics-based pore compressibility tests into compliant rock compaction tables for reservoir simulators using Dimensionless Stress to Pore Pressure Conversion, to thereby transfer geomechanical changes due to confining stress into expressions of geomechanical changes due to pore pressure.

One or more methods for validating reservoir simulation models. At least one of the methods include determining one or more time segments of fluid recovery of a reservoir by analyzing a production history of the reservoir; running a simulation model, for the first time segment, to generate one or more drainage volumes; generating, for a first time segment, a plurality of grid regions along one or more no-flow boundaries of the one or more drainage volumes; generating, for the first time segment, a plurality of sector models corresponding to the plurality of grid regions; and performing, for the first time segment, a history matching process corresponding to a time phase simultaneously on each of the plurality of sector models to generate, for each of the sector models, a history matching output.