A Fully Mass Conservative Numerical Method for Multiphase Flow in Fractured Porous Reservoirs

Numerical methods that are mass conservative, computationally stable, and efficient are essential for simulating multiphase flow in fractured porous reservoirs. In this paper, a fully mass conservative numerical method that meets these requirements is proposed and analyzed. Unlike in the implicit pr...

Full description

Saved in:
Bibliographic Details
Published in:Transport in porous media 2021-09, Vol.139 (2), p.171-184
Main Authors: Cai, Hailiang, Li, Peichao, Feng, Meng, Hao, Youzhi, Lu, Detang, Xian, Yuxi
Format: Article
Language:English
Subjects:
Capillary pressure
Civil Engineering
Classical and Continuum Physics
Computational efficiency
Earth and Environmental Science
Earth Sciences
Geotechnical Engineering & Applied Earth Sciences
Hydrogeology
Hydrology/Water Resources
Industrial Chemistry/Chemical Engineering
Initial conditions
Iterative methods
Multiphase flow
Numerical analysis
Numerical methods
Permeability
Reservoirs
Robustness (mathematics)
Saturation
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Numerical methods that are mass conservative, computationally stable, and efficient are essential for simulating multiphase flow in fractured porous reservoirs. In this paper, a fully mass conservative numerical method that meets these requirements is proposed and analyzed. Unlike in the implicit pressure and explicit saturation (IMPES) method, in this method, the pressure and saturation equations are solved sequentially in each iteration step. In this framework, the calculation of the saturation-related parameters no longer depends on the initial conditions of the current time step, but on the results of the current and previous iterations. Through this treatment, the discontinuities of the saturation and capillary pressure in the pre- and post-two time steps can be overcome to achieve fully mass conservation. Two numerical examples are designed to verify the robustness and efficiency of the presented method. The numerical results indicate that this new method is fully mass conservative, computationally more stable, and more efficient than the IMPES method. Furthermore, the proposed method is suitable not only for reservoirs with homogeneous permeability distributions but also for reservoirs with heterogeneous permeability distributions, which greatly improves the applicability of the new method.
ISSN:0169-3913
1573-1634
DOI:10.1007/s11242-021-01636-9