Modeling Reservoir Geometry With Irregular Grids

Abstract

The widely used Cartesian coordinate grid has some disadvantages in the description of boundaries, faults and discontinuities. In addition, a five-point scheme can cause significant grid orientation effects. A ninepoint scheme reduces this effect but makes the treatment of boundaries and heterogeneities more difficult.Orthogonal curvilinear coordinate systems have improved the modeling of reservoir shape and flow geometry. Mathematically they are based on a coordinate transformation and discretization of the transformed equations in the usual manner. Due to the lack of orthogonality additional mixed derivative terms are introduced which are difficult to discretize and therefore are usually neglected. However, because grid lines have to be coordinate lines, strictly orthogonal systems are not flexible enough to describe reservoirs of complicated shape. This paper describes a practical method for using irregular or locally irregular grids in reservoir simulation with the advantages of flexible approximation of reservoir geometry, simple treatment of boundary conditions and reduced grid orientation effects.Finite difference equations are set up by the so-called balance method. This method uses an integral formulation of the reservoir model equations equivalent to the commonly used differential equations. Integrating over grid blocks results in "balance equations" for each References and illustrations at end of paper 37 block. This can be done for various types of networks formed by triangles, convex quadrangles, polar meshes, curvilinear or locally refined grids. HEINRICHS 1 has proved consistence for this discretization scheme. The mesh can be refined locally. Well locations can be selected as mesh points to ensure that they are situated in the center of grid blocks. For triangular grids, the more isotropic distribution of grid points diminishes the orientation effect significantly.Numerical examples are presented comparing the proposed difference scheme with a nine-point Cartesian scheme . The performance of the method is illustrated by symmetry elements and complex simulation problems.

Grid Orientation Revisited: Near-well, Early-time Effects and Solution Coupling Methods

“…Other gridding methods have been proposed with the express purpose of reducing GOE such as Pruess and Bodvarsson (1983) and Heinemann et al (1991). Currently, irregular grids are often suggested for handling GOE in the literature since they can seemingly reduce directionally preferred transport directions.…”

q-Families of CVD(MPFA) Schemes on General Elements: Numerical Convergence and the Maximum Principle

“…In general a full tensor equation arises whenever the computational grid is non-aligned with the principal axes of the local tensor field. A full tensor can occur when representing cross bedding, modeling any anisotropic medium that is non-aligned with the computational grid [10], using non K-orthogonal [11,12] or unstructured grids and when upscaling rock properties from fine scale diagonal tensor simulation to the grid block scale [13]. Consequently, the standard diagonal tensor simulator will suffer from an inconsistent O(1) error in flux when applied to cases involving these major features.…”

Kinetic simulation of methane hydrate formation and dissociation in porous media

“…Finite volume difference method (FVDM) is used in spatial discretization of the governing equations (Heinemann et al, 1989;Pruess and Moridis, 1999). The discretized grid blocks are arranged with the standard method described by .…”