Analytical Solution to the Riemann Problem of Three-Phase Flow in Porous Media

In this paper we study one-dimensional three-phase flow through porous media of immiscible, incompressible fluids. The model uses the common multiphase flow extension of Darcy’s equation, and does not include gravity and capillarity effects. Under these conditions, the mathematical problem reduces to a 2 × 2 system of conservation laws whose essential features are: (1) the system is strictly hyperbolic; (2) both characteristic fields are nongenuinely nonlinear, with single, connected inflection loci. These properties, which are natural extensions of the two-phase flow model, ensure that the solution is physically sensible. We present the complete analytical solution to the Riemann problem (constant initial and injected states) in detail, and describe the characteristic waves that may arise, concluding that only nine combinations of rarefactions, shocks and rarefaction-shocks are possible. We demonstrate that assuming the saturation paths of the solution are straight lines may result in inaccurate predictions for some realistic systems. Efficient algorithms for computing the exact solution are also given, making the analytical developments presented here readily applicable to interpretation of lab displacement experiments, and implementation of streamline simulators.

[1]  Tai-Ping Liu,et al.  The Riemann problem for general systems of conservation laws , 1975 .

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

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

[4]  M. J. Oak,et al.  Three-phase relative permeability of Berea sandstone , 1990 .

[5]  H. L. Stone Estimation of Three-Phase Relative Permeability And Residual Oil Data , 1973 .

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

[7]  S. Axler Linear Algebra Done Right , 1995, Undergraduate Texts in Mathematics.

[8]  Morris Muskat,et al.  Physical principles of oil production , 1949 .

[9]  L. Baker Three-Phase Relative Permeability Correlations , 1988 .

[10]  Joshua Kiddy K. Asamoah,et al.  Fractal–fractional age-structure study of omicron SARS-CoV-2 variant transmission dynamics , 2022, Partial Differential Equations in Applied Mathematics.

[11]  Cass T. Miller,et al.  Multiphase flow and transport modeling in heterogeneous porous media: challenges and approaches , 1998 .

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

[13]  R. LeVeque Numerical methods for conservation laws , 1990 .

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

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

[16]  Jérôme Jaffré,et al.  Estimation of relative permeabilities in three-phase flow in porous media , 1999 .

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

[18]  H. L. Stone Probability Model for Estimating Three-Phase Relative Permeability , 1970 .

[19]  R. Aris First-order partial differential equations , 1987 .

[20]  Martin J. Blunt,et al.  A 3D Field-Scale Streamline-Based Reservoir Simulator , 1997 .

[21]  G. I. Barenblatt Scaling: Self-similarity and intermediate asymptotics , 1996 .

[22]  G. I. Barenblatt,et al.  Theory of Fluid Flows Through Natural Rocks , 1990 .

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

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

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

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

[27]  L. A. Rapoport,et al.  Mechanism of Water Flooding in the Presence of Free Gas , 1956 .

[28]  Barbara Lee Keyfitz,et al.  Nonlinear evolution equations that change type , 1990 .

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

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

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

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

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

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

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

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

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

[38]  John A. Trangenstein Three-Phase Flow with Gravity , 1988 .

[39]  J. Bear Dynamics of Fluids in Porous Media , 1975 .

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

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

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

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

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

[45]  W. Brent Lindquist Current progress in hyperbolic systems, Riemann problems and computations : proceedings of the AMS-IMS-SIAM joint summer research conference held July 16-22, 1988 with support from the National Science Foundation and the Office of Naval Research , 1989 .