Finite Difference Simulation of Geologically Complex Reservoirs With Tensor Permeabilities

Abstract

The grid block permeabilities used in reservoir simulation are commonly determined via the upscaling of a fine scale geostatistical reservoir description. Though it is well established that permeabilities computed in this manner are in general full tensor quantities, most finite difference reservoir simulators still treat permeability as a diagonal tensor. In this paper, we implement a capability to handle full tensor permeabilities in a general purpose finite difference simulator and apply this capability to the modeling of several complex geological systems. We formulate a flux continuous approach for the pressure equation using a method analogous to that of previous researchers (Edwards and Rogers; Aavatsmark et al.), consider methods for upwinding in multiphase flow problems, and additionally discuss some relevant implementation and reservoir characterization issues. The accuracy of the finite difference formulation, assessed through comparisons to an accurate finite element approach, is shown to be generally good, particularly for immiscible displacements in heterogeneous systems. The formulation is then applied to the simulation of upscaled descriptions of several geologically complex reservoirs involving crossbedding and extensive fracturing. The method performs quite well for these systems and is shown to accurately capture the effects of the underlying geology. Finally, the significant errors which can be incurred through inaccurate representation of the full permeability tensor are demonstrated for several cases. Introduction Recent advances in reservoir characterization permit the construction of realistic, highly detailed, heterogeneous reservoir descriptions. Such models typically contain far too many grid blocks to simulate directly and therefore require some type of upscaling before they can be used for reservoir simulation. The most important of the upscaled rock properties, for purposes of flow simulation, is the absolute permeability. Accurate procedures for the scale up of permeability generate full tensor permeabilities on the coarse scale, even in cases where the underlying fine scale permeability description is isotropic. Therefore, simulation models generated through scale up of complex reservoir descriptions will in general be characterized by full tensor permeabilities. For many models, however, the off-diagonal components of the effective (or equivalent grid block) permeability tensors can be expected to be small relative to the diagonal components and can generally be ignored. This will typically be the case, for example, when the fine scale permeability is correlated along the coordinate directions (e.g., strictly layered systems) or is correlated nearly isotropically. However, for other types of systems the cross terms of the permeability tensor can be expected to be quite significant. These include formations containing complex crossbedding, dipping layers not aligned with the coordinate system, or extensive fracturing. In these cases the upscaled simulation model will contain full tensor permeabilities with significant off-diagonal terms, which must be accommodated by the simulator if reservoir performance is to be predicted accurately. Our intent here is to develop an approach for the modeling of complex geological systems using a general purpose reservoir simulator. Toward that end, we first develop and implement a formulation for full tensor permeability models, applicable for curvilinear grids, into a general purpose finite difference reservoir simulator. Second, we apply this formulation to the simulation of flow through a variety of upscaled geologic descriptions. P. 253^

Hierarchical modeling of flow in naturally fractured formations with multiple length scales

“…The grid block effective permeabilities are used with a finite difference simulator to compute flow through the fracture and rock system. The simulator is enhanced to use a control-volume finite difference formulation for general tensor permeability input (i.e., nine-point stencil for 2-D and 27-point stencil for 3-D; the stencil denotes the grid points involved in the computation of a grid variable) [Aavatsmark et al, 1994;Edwards and Rogers, 1994;Lee et al, 1998Lee et al, , 1999].…”

Grid design and numerical modeling of multiphase flow in complex reservoirs using orthogonal collocation schemes

“…(3) comprise cell edge-based quantities. Several new schemes have been proposed and developed for full tensor simulation in recent years, which included the works of some investigators [4,8,10,12,17,21,22]. These schemes all employ nine-point operators for the pressure equation on logically Cartesian grids in two dimensions.…”

Limitations to upscaling of groundwater flow models dominated by surface water interaction

“…An important aspect of coarsening is that K c becomes often anisotropic, even if the original hydraulic permeability K is isotropic. Recently, methods became available that cope with the off-diagonal components in a tensor [Aavatsmark et al, 1996;Lee et al, 1998;Anderman et al, 2002].…”