Practical implementation of the fractional flow approach to multi-phase flow simulation

Fractional flow formulations of the multi-phase flow equations exhibit several attractive attributes for numerical simulations. The governing equations are a saturation equation having an advection diffusion form, for which characteristic methods are suited, and a global pressure equation whose form is elliptic. The fractional flow approach to the governing equations is compared with other approaches and the implication of equation form for numerical methods discussed. The fractional flow equations are solved with a modified method of characteristics for the saturation equation and a finite element method for the pressure equation. An iterative algorithm for determination of the general boundary conditions is implemented. Comparisons are made with a numerical method based on the twopressure formulation of the governing equations. While the fractional flow approach is attractive for model problems, the performance of numerical methods based on these equations is relatively poor when the method is applied to general boundary conditions. We expect similar difficulties with the fractional flow approach for more general problems involving heterogenous material properties and multiple spatial dimensions. q 1999 Elsevier Science Ltd. All rights reserved.

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

[2]  Chin-Fu Tsang,et al.  A perspective on the numerical solution of convection‐dominated transport problems: A price to pay for the easy way out , 1992 .

[3]  M. Espedal,et al.  Characteristic, local grid refinement techniques for reservoir flow problems , 1992 .

[4]  P. Milly,et al.  A mass-conservative procedure for time-stepping in models of unsaturated flow , 1985 .

[5]  Richard W. Healy,et al.  A finite‐volume Eulerian‐Lagrangian Localized Adjoint Method for solution of the advection‐dispersion equation , 1993 .

[6]  Richard E. Ewing,et al.  Incorporation of Mixed Finite Element Methods in Compositional Simulation for Reduction of Numerical Dispersion , 1983 .

[7]  R. Ewing,et al.  Characteristic adaptive subdomain methods for reservoir flow problems , 1990 .

[8]  M. Espedal,et al.  On the numerical solution of non-linear reservoir flow models with gravity , 1995 .

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

[10]  J. J. Douglas,et al.  Finite Difference Methods for Two-Phase Incompressible Flow in Porous Media , 1983 .

[11]  P. Forsyth,et al.  Variable spatial and temporal weighting schemes for use in multi-phase compositional problems , 1996 .

[12]  Larry C. Young,et al.  A study of spatial approximations for simulating fluid displacements in petroleum reservoirs , 1984 .

[13]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

[14]  H. J. Morel-Seytoux,et al.  Two-Phase Flows in Porous Media , 1973 .

[15]  Jonathan F. Sykes,et al.  Compositional simulation of groundwater contamination by organic compounds: 2. Model applications , 1993 .

[16]  H. J. Welge,et al.  A Simplified Method for Computing Oil Recovery by Gas or Water Drive , 1952 .

[17]  G. Moridis,et al.  Secondary Water Recovery by Air Injection: 1. The Concept and the Mathematical and Numerical Model , 1991 .

[18]  Peter A. Forsyth,et al.  A positivity preserving method for simulation of steam injection for NAPL site remediation , 1993 .

[19]  S. E. Buckley,et al.  Mechanism of Fluid Displacement in Sands , 1942 .

[20]  Michel Vauclin,et al.  Experimental and numerical analysis of two-phase infiltration in a partially saturated soil , 1986 .

[21]  Richard E. Ewing,et al.  Mixed finite element approximation of phase velocities in compositional reservoir simulation , 1984 .

[22]  Bernhard A. Schrefler,et al.  A FULLY COUPLED MODEL FOR WATER FLOW AND AIRFLOW IN DEFORMABLE POROUS MEDIA , 1993 .

[23]  Randel Haverkamp,et al.  A Comparison of Numerical Simulation Models For One-Dimensional Infiltration1 , 1977 .

[24]  H. Morel‐Seytoux,et al.  A Two-Phase Numerical Model for Prediction of Infiltration: Case of an Impervious Bottom , 1985 .

[25]  Emil O. Frind,et al.  Two‐phase flow in heterogeneous porous media: 1. Model development , 1991 .

[26]  Emil O. Frind,et al.  Comparative error analysis in finite element formulations of the advection-dispersion equation , 1985 .

[27]  G. Chavent.  Mathematical models and finite elements for reservoir simulation , 1986 .

[28]  Magne S. Espedal,et al.  Macrodispersion for two-phase, immiscible flow in porous media , 1994 .

[29]  Magne S. Espedal,et al.  Continuous-time finite element analysis of multiphase flow in groundwater hydrology , 1994 .

[30]  K. Pruess,et al.  TOUGH User's Guide , 1987 .

[31]  Magne S. Espedal,et al.  Multiphase Flow Simulation with Various Boundary Conditions , 1994 .

[32]  Alain Bourgeat,et al.  TWO-PHASE FLOW IN HETEROGENEOUS POROUS MEDIA , 1991 .

[33]  Jonathan F. Sykes,et al.  Compositional simulation of groundwater contamination by organic compounds: 1. Model development and verification , 1993 .

[34]  Jonathan F. Sykes,et al.  Modeling the transport of volatile organics in variably saturated media , 1989 .

[35]  Richard E. Ewing,et al.  Simulation of multiphase flows in porous media , 1991 .

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

[37]  Philip John Binning,et al.  A mass conservative numerical solution for two‐phase flow in porous media with application to unsaturated flow , 1992 .

[38]  T. F. Russell,et al.  Modeling of multiphase multicontaminant transport in the subsurface , 1995 .

[39]  Richard E. Ewing,et al.  Eulerian-Lagrangian Localized Adjoint Methods for a Nonlinear Advection-Diffusion Equation , 1994 .

[40]  R. Ewing,et al.  Characteristics Petrov-Galerkin subdomain methods for two-phase immiscible flow , 1987 .

[41]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[42]  George F. Pinder,et al.  On the Simulation of Nonaqueous Phase Organic Compounds in the Subsurface , 1986 .

[43]  K. Rathfelder,et al.  Mass balance errors in modeling two-phase immiscible flows: causes and remedies , 1993 .

[44]  Magnus Wangen,et al.  Vertical migration of hydrocarbons modelled with fractional flow theory , 1993 .