On the selection of primary variables in numerical formulation for modeling multiphase flow in porous media.

Selecting the proper primary variables is a critical step in efficiently modeling the highly nonlinear problem of multiphase subsurface flow in a heterogeneous porous-fractured media. Current simulation and ground modeling techniques consist of (1) spatial discretization of mass and/or heat conservation equations using finite difference or finite element methods; (2) fully implicit time discretization; (3) solving the nonlinear, discrete algebraic equations using a Newton iterative scheme. Previous modeling efforts indicate that the choice of primary variables for a Newton iteration not only impacts computational performance of a numerical code, but may also determine the feasibility of a numerical modeling study in many field applications. This paper presents an analysis and general recommendations for selecting primary variables in simulating multiphase, subsurface flow for one-active phase (Richards' equation), two-phase (gas and liquid) and three-phase (gas, water and nonaqueous phase liquid or NAPL) conditions. In many cases, a dynamic variable switching or variable substitution scheme may have to be used in order to achieve optimal numerical performance and robustness. The selection of primary variables depends in general on the sensitivity of the system of equations to the variables selected at given phase and flow conditions. We will present a series of numerical tests and large-scale field simulation examples, including modeling one (active)-phase, two-phase and three-phase flow problems in multi-dimensional, porous-fractured subsurface systems.

[1]  J. Istok.  Groundwater Modeling by the Finite Element Method , 1989 .

[2]  Yu-Shu Wu,et al.  A consistent approach for applying numerical boundary conditions for multiphase subsurface flow , 1996 .

[3]  J. Parker,et al.  Correction to “A parametric model for constitutive properties governing multiphase flow in porous media” by J. C. Parker, R. J. Lenhard, and T. Kuppusamy , 1987 .

[4]  Peter A. Forsyth,et al.  Three‐dimensional modelling of steam flush for DNAPL site remediation , 1994 .

[5]  Karsten Pruess,et al.  Software qualification of selected TOUGH2 modules , 1996 .

[6]  T. Narasimhan,et al.  AN INTEGRATED FINITE DIFFERENCE METHOD FOR ANALYZING FLUID FLOW IN POROUS MEDIA , 1976 .

[7]  K. H. Coats,et al.  Simulation of three-dimensional, two-phase flow in oil and gas reservoirs , 1967 .

[8]  Hubert J. Morel-Seytoux,et al.  A Two‐Phase Numerical Model for Prediction of Infiltration: Applications to a Semi‐Infinite Soil Column , 1985 .

[9]  Jack C. Parker,et al.  A parametric model for constitutive properties governing multiphase flow in porous media , 1987 .

[10]  W. E. Wilson,et al.  Conceptual hydrologic model of flow in the unsaturated zone, Yucca Mountain, Nevada , 1984 .

[11]  Peter A. Forsyth,et al.  Numerical simulation of gas venting for NAPL site remediation , 1991 .

[12]  K. Pruess,et al.  TOUGH2-A General-Purpose Numerical Simulator for Multiphase Fluid and Heat Flow , 1991 .

[13]  Karsten Pruess,et al.  Numerical modeling of steam injection for the removal of nonaqueous phase liquids from the subsurface: 2. Code validation and application , 1992 .

[14]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[15]  D. W. Peaceman.  Fundamentals of numerical reservoir simulation , 1977 .

[16]  L. K. Kuiper,et al.  Numerical simulation of NAPL flow in the subsurface , 1998 .

[17]  J. Parker,et al.  An efficient finite element method for modeling multiphase flow , 1989 .

[18]  Yu-Shu Wu,et al.  A three-dimensional multiphase flow model for assesing NAPL contamination in porous and fractured media, 1. Formulation , 1994, Journal of Contaminant Hydrology.

[19]  Khalid Aziz,et al.  Use of Voronoi Grid in Reservoir Simulation , 1994 .

[20]  H. H. Rachford,et al.  Numerical Calculation of Multidimensional Miscible Displacement , 1962 .

[21]  Yu-Shu Wu,et al.  A three-dimensional multiphase flow model for assessing NAPL contamination in porous and fractured media, 2. Porous medium simulation examples , 1994 .

[22]  Hydrogeology of the unsaturated zone, North Ramp area of the Exploratory Studies Facility, Yucca Mountain, Nevada , 1999 .

[23]  R D Hoak,et al.  Water Resources , 1961, Nature.

[24]  Three-Dimensional Geothermal Reservoir Simulation , 1978 .

[25]  Peter A. Forsyth,et al.  Robust numerical methods for saturated-unsaturated flow with dry initial conditions in heterogeneous media , 1996 .

[26]  Marge Kars,et al.  Resources for Research , 2001, Lancet.

[27]  Charles R. Faust,et al.  Transport of Immiscible Fluids Within and Below the Unsaturated Zone: A Numerical Model , 1985 .

[28]  Philip John Binning,et al.  Practical implementation of the fractional flow approach to multi-phase flow simulation , 1999 .

[29]  P. Huyakorn,et al.  Techniques for Making Finite Elements Competitve in Modeling Flow in Variably Saturated Porous Media , 1984 .

[30]  Yu-Shu Wu,et al.  A multiple-porosity method for simulation of naturally fractured petroleum reservoirs , 1988 .

[31]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[32]  P. A. Forsytha,et al.  Nonlinear iteration methods for nonequilibrium multiphase subsurface flow , 1998 .

[33]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[34]  P. Witherspoon,et al.  Numerical modeling of steam injection for the removal of nonaqueous phase liquids from the subsurface. 1. Numerical formulation , 1992 .

[35]  R. L. Cooley.  Some new procedures for numerical solution of variably saturated flow problems , 1983 .

[36]  Peter A. Forsyth.  Simulation of nonaqueous phase groundwater contamination , 1988 .

[37]  George F. Pinder,et al.  A Multiphase Approach to the Modeling of Porous Media Contamination by Organic Compounds: 1. Equation Development , 1985 .

[38]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .

[39]  Gs Bodvarsson,et al.  DEVELOPMENT AND CALIBRATION OF THE THREE-DIMENSIONAL SITE-SCALE UNSATURATED ZONE MODEL OF YUCCA MOUNTAIN, NEVADA , 1996 .

[40]  Gudmundur S. Bodvarsson,et al.  A site-scale model for fluid and heat flow in the unsaturated zone of Yucca Mountain, Nevada , 1999 .

[41]  P. Berg.  Long-term simulation of water movement in soils using mass-conserving procedures , 1999 .

[42]  Howard B. Bradley,et al.  Petroleum engineering handbook , 1987 .

[43]  Gudmundur S. Bodvarsson,et al.  Preliminary development of the LBL/USGS three-dimensional site-scale model of Yucca Mountain, Nevada , 1995 .

[44]  Yu-Shu Wu,et al.  A vertical equilibrium model for assessing nonaqueous phase liquid contamination and remediation of groundwater systems , 1994 .

[45]  Peter A. Forsyth,et al.  A Control Volume Finite Element Approach to NAPL Groundwater Contamination , 1991, SIAM J. Sci. Comput..

[46]  George F. Pinder,et al.  A Multiphase Approach to the Modeling of Porous Media Contamination by Organic Compounds: 2. Numerical Simulation , 1985 .

[47]  G. W. Thomas,et al.  Principles of Hydrocarbon Reservoir Simulation , 1981 .

[48]  Yu-Shu Wu,et al.  A virtual node method for handling well bore boundary conditions in modeling multiphase flow in porous and fractured media. , 2000 .

[49]  R. B. Scott,et al.  Preliminary geologic map of Yucca Mountain, Nye County, Nevada, with geologic sections , 1984 .

[50]  FINITE ELEMENT ANALYSIS OF LAMINAR/ TURBULENT FLOW IN POROUS AND FRACTURED MEDIA , 1985 .

[51]  C. S. Simmons,et al.  Recent Studies of Flow and Transport in the Vadose Zone , 1991 .

[53]  T. N. Narasimhan,et al.  A PRACTICAL METHOD FOR MODELING FLUID AND HEAT FLOW IN FRACTURED POROUS MEDIA , 1985 .

[54]  Jim Douglas,et al.  A Method for Calculating Multi-Dimensional Immiscible Displacement , 1959 .