Finite Difference Simulation of Geologically Complex Reservoirs With Tensor Permeabilities

Abstract

Summary: The gridblock permeabilities used in reservoir simulation are commonly determined through the upscaling of a fine scale geostatistical reservoir description. Though it is well established that per- meabilities 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 con- tinuous approach for the pressure equation by use of a method analogous to that of previous researchers (Edwards and Rogers1; Aavatsmark et al.2), consider methods for upwinding in multiphase flow problems, and additionally discuss some relevant implemen- tation 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 sys- tems. 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 capture the effects of the underlying geology accurately. Finally, the significant errors that can be incurred through inaccurate representation of the full per- meability tensor are demonstrated for several cases.