A front-tracking method for the simulation of three-phase flow in porous media

Abstract Under certain physically reasonable assumptions, three-phase flow of immiscible, incompressible fluids can be described by a 2×2 nongenuinely nonlinear, hyperbolic system. We combine analytical solutions to the corresponding Riemann problem with an efficient front-tracking method to study Cauchy and initial-boundary value problems. Unlike finite difference methods, the front-tracking method treats all waves as discontinuities by evolving shocks exactly and approximating rarefactions by small entropy-violating discontinuities. This way, the method can track individual waves and give very accurate (or even exact) resolution of discontinuities. We demonstrate the applicability of the method through several numerical examples, including a streamline simulation of a water-alternating-gas (WAG) injection process in a three-dimensional, heterogeneous, shallow-marine formation.

[1]  P. Lax.  Hyperbolic systems of conservation laws II , 1957 .

[2]  C. Dafermos.  Polygonal approximations of solutions of the initial value problem for a conservation law , 1972 .

[3]  Filtration of a three-phase liquid (Buckley-Leverett model) , 1972 .

[4]  Tai-Ping Liu.  The Riemann problem for general 2×2 conservation laws , 1974 .

[5]  Gary A. Pope,et al.  The Application of Fractional Flow Theory to Enhanced Oil Recovery , 1980 .

[6]  F. Helfferich,et al.  Theory of Multicomponent, Multiphase Displacement in Porous Media , 1981 .

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

[8]  John A. Trangenstein,et al.  Conservation laws of mixed type describing three-phase flow in porous media , 1986 .

[9]  Michael Shearer,et al.  Loss of strict hyperbolicity of the Buckley-Leverett equations for three phase flow in a porous medium , 1987 .

[10]  F. Brezzi,et al.  Mathematical models and finite elements for reservoir simulation : single phase, multiphase and multicomponent flows through porous media , 1988 .

[11]  D. W. Pollock.  Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[12]  L. Holden,et al.  A NUMERICAL METHOD FOR FIRST ORDER NONLINEAR SCALAR CONSERVATION LAWS IN ONE-DIMENSION , 1988 .

[13]  John A. Trangenstein,et al.  Loss of real characteristics for models of three-phase flow in a porous medium , 1989 .

[14]  F. Fayers.  Extension of Stone's Method 1 and conditions for real characteristics in three-phase flow , 1989 .

[15]  Lars Holden,et al.  On the strict hyperbolicity of the Buckley-Leverett equations for three-phase flow in a porous medium , 1989 .

[16]  Nils Henrik Risebro,et al.  Front Tracking Applied to a Nonstrictly Hyperbolic System of Conservation Laws , 1991, SIAM J. Sci. Comput..

[17]  A. Bressan.  Global solutions of systems of conservation laws by wave-front tracking , 1992 .

[18]  A. H. Falls,et al.  Features of three-component, three-phase displacement in porous media , 1992 .

[19]  A. H. Falls,et al.  Theory of Three-Component, Three-Phase Displacement in Porous Media , 1992 .

[20]  N. Risebro.  A front-tracking alternative to the random choice method , 1993 .

[21]  Aslak Tveito,et al.  A Conservative Front Tracking Scheme for 1D Hyperbolic Conservation Laws , 1993 .

[22]  Abraham S. Grader,et al.  Simulation of three-phase displacement experiments , 1996 .

[23]  L. Lake,et al.  Enhanced Oil Recovery , 2017 .

[24]  Martin J. Blunt,et al.  Theoretical Analysis of Three Phase Flow Experiments in Porous Media , 1996 .

[25]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[26]  A. Bressan,et al.  Uniqueness of Weak Solutions to Systems of Conservation Laws , 1997 .

[27]  F. Fayers,et al.  Solutions to the Three-Phase Buckley-Leverett Problem , 1996 .

[28]  F. Fayers,et al.  Mathematical Properties of Three-Phase Flow Equations , 1997 .

[29]  B. Piccoli,et al.  Well-posedness of the Cauchy problem for × systems of conservation laws , 2000 .

[30]  Fabio Ancona,et al.  A Note on the Riemann Problem for General n × n Conservation Laws , 2001 .

[31]  Alexander Kurganov,et al.  Semidiscrete Central-Upwind Schemes for Hyperbolic Conservation Laws and Hamilton-Jacobi Equations , 2001, SIAM J. Sci. Comput..

[32]  Michael Andrew Christie,et al.  Tenth SPE Comparative Solution Project: a comparison of upscaling techniques , 2001 .

[33]  Dan Marchesin,et al.  Wave Structure in WAG Recovery , 2001 .

[34]  Inga Berre,et al.  A Streamline Front Tracking Method for Two- and Three-Phase Flow Including Capillary Forces , 2002 .

[35]  M. Blunt,et al.  Elliptic Regions and Stable Solutions for Three-Phase flow in Porous Media , 2002 .

[36]  H. Holden,et al.  Front Tracking for Hyperbolic Conservation Laws , 2002 .

[37]  Kevin Zumbrun,et al.  Capillary instability in models for three-phase flow , 2002 .

[38]  T. Patzek,et al.  Relative Permeabilities for Strictly Hyperbolic Models of Three-Phase Flow in Porous Media , 2004 .

[39]  Ruben Juanes,et al.  A Front-Tracking Method for Hyperbolic Three-Phase Models , 2004 .

[40]  Ruben Juanes,et al.  Analytical Solution to the Riemann Problem of Three-Phase Flow in Porous Media , 2004 .

[41]  Ruben Juanes,et al.  A Front-Tracking Method for Efficient Simulation of Miscible Gas Injection Processes , 2005 .

[42]  P. Lax.  Hyperbolic systems of conservation laws , 2006 .