Modeling fractured and faulted regions: Local grid refinement methods for implicit solvers

In order to give a full 3D description of reservoir fault zones, we use a non-overlapping, non-matching local grid refinement (LGR) technique in the multi-phase, multi-component hydrocarbon migration simulator Som (Secondary Oil Migration). The resulting small spatial scales severely restrict the time steps if the mass transport equations are solved explicitly. Thus, an implicit formulation is derived where the Jacobian terms are analytical expressions based on binary mixture thermodynamics. We will present details of the implementation of this method in the LGR framework using Newton-GMRES and a multigrid-like, two-level iterative solver.

[1]  Magne S. Espedal,et al.  Parallelization of a Compositional Reservoir Simulator , 2000 .

[2]  Yousef Saad,et al.  Hybrid Krylov Methods for Nonlinear Systems of Equations , 1990, SIAM J. Sci. Comput..

[3]  J. W. Watts,et al.  A Compositional Formulation of the Pressure and Saturation Equations , 1986 .

[4]  Richard E. Ewing,et al.  Numerical treatment of multiphase flows in porous media : proceedings of the international workshop, held at Beijing, China, 2-6 August 1999 , 2000 .

[5]  Magne S. Espedal,et al.  Implicit treatment of molar mass equations in secondary oil migration , 2002 .

[6]  Tom Manzocchi,et al.  The representation of two phase fault-rock properties in flow simulation models , 2002, Petroleum Geoscience.

[7]  Rune Teigland,et al.  On some variational acceleration techniques and related methods for local refinement , 1998 .

[8]  William Gropp,et al.  Skjellum using mpi: portable parallel programming with the message-passing interface , 1994 .

[9]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[10]  H. Reme,et al.  Use of local grid refinement and a Galerkin technique to study secondary migration in fractured and faulted regions , 1999 .

[11]  Zhangxin Chen,et al.  Analysis of a Compositional Model for Fluid Flow in Porous Media , 2000, SIAM J. Appl. Math..

[12]  Eva Farkas,et al.  General Purpose Compositional Model , 1985 .

[13]  Louis J. Durlofsky,et al.  Representation of Fault Zone Permeability in Reservoir Flow Models , 2001 .

[14]  Magne S. Espedal,et al.  Numerical Simulation of Compositional Fluid Flow in Porous Media , 2000 .

[15]  Zhangxin Chen.  Formulations and Numerical Methods of the Black Oil Model in Porous Media , 2000, SIAM J. Numer. Anal..

[16]  Ivar Aavatsmark,et al.  Discretization on Unstructured Grids for Inhomogeneous, Anisotropic Media. Part I: Derivation of the Methods , 1998, SIAM J. Sci. Comput..

[17]  Ivar Aavatsmark,et al.  Discretization on Unstructured Grids For Inhomogeneous, Anisotropic Media. Part II: Discussion And Numerical Results , 1998, SIAM J. Sci. Comput..

[18]  M. Edwards.  Elimination of Adaptive Grid Interface Errors in the Discrete Cell Centered Pressure Equation , 1996 .

[19]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[20]  E. Reiso,et al.  Control-Volume Discretization Method for Quadrilateral Grids with Faults and Local Refinements , 2001 .

[21]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[22]  Ivar Aavatsmark,et al.  Discretization on Non-Orthogonal, Quadrilateral Grids for Inhomogeneous, Anisotropic Media , 1996 .

[23]  John A. Trangenstein,et al.  Mathematical structure of the black-oil model for petroleum reservoir simulation , 1989 .

[24]  Magne S. Espedal,et al.  Implicit Treatment of Compositional Flow , 2004 .

[25]  S. McCormick,et al.  The fast adaptive composite grid (FAC) method for elliptic equation , 1986 .

[26]  Hilde Reme,et al.  Parallelization of a Compositional Simulator with a Galerkin Coarse/Fine Method , 1999, Euro-Par.