Communications on Applied Mathematics and Computation ›› 2023, Vol. 5 ›› Issue (2): 946-964.doi: 10.1007/s42967-022-00200-6

Previous Articles    

Numerical Approach of a Coupled Pressure-Saturation Model Describing Oil-Water Flow in Porous Media

Paula Luna1, Arturo Hidalgo2   

  1. 1 ISIS Neutron and Muon Source, Rutherford Appleton Laboratory, Didcot OX11,0QX, UK;
    2 Departamento de Ingeniería Geológica y Minera, Escuela Técnica Superior de Ingenieros de Minas y Energía, Center for Computational Simulation, Universidad Politécnica de Madrid, Ríos Rosas, 21, 28003 Madrid, Spain
  • Received:2021-09-10 Revised:2022-06-14 Online:2023-06-20 Published:2023-05-26
  • Contact: Arturo Hidalgo, arturo.hidalgo@upm.es;Paula Luna, paula.luna.dapica@gmail.com E-mail:arturo.hidalgo@upm.es;paula.luna.dapica@gmail.com
  • Supported by:
    This work has been carried out under the project:Development of a numerical multiphase flow tool for applications to petroleum industry, with the Spanish oil company REPSOL.

Abstract: Two-phase flow in porous media is a very active field of research, due to its important applications in groundwater pollution, CO2 sequestration, or oil and gas production from petroleum reservoirs, just to name a few of them. Fractional flow equations, which make use of Darcy's law, for describing the movement of two immiscible fluids in a porous medium, are among the most relevant mathematical models in reservoir simulation. This work aims to solve a fractional flow model formed by an elliptic equation, representing the spatial distribution of the pressure, and a hyperbolic equation describing the space-time evolution of water saturation. The numerical solution of the elliptic part is obtained using a finite-element (FE) scheme, while the hyperbolic equation is solved by means of two different numerical approaches, both in the finite-volume (FV) framework. One is based on a monotonic upstream-centered scheme for conservation laws (MUSCL)-Hancock scheme, whereas the other makes use of a weighted essentially non-oscillatory (ENO) reconstruction. In both cases, a first-order centered (FORCE)-α numerical scheme is applied for intercell flux reconstruction, which constitutes a new contribution in the field of fractional flow models describing oil-water movement. A relevant feature of this work is the study of the effect of the parameter α on the numerical solution of the models considered. We also show that, in the FORCE-α method, when the parameter α increases, the errors diminish and the order of accuracy is more properly attained, as verified using a manufactured solution technique.

Key words: Two-phase flow, Reservoir simulation, Porous media, FORCE-α, Finite volume

CLC Number: