Numerical simulation of 1-D oil and water displacements in petroleum reservoirs using the correction procedure via reconstruction (CPR) method

The focus of this work is to investigate and to apply, for the first time, a very high-resolution correction procedure via reconstruction (CPR) numerical discretization technique for the hyperbolic saturation equation that describes 1-D oil-water displacement through heterogeneous porous media. The...

Full description

Saved in:
Bibliographic Details
Published in:Computational geosciences 2020-02, Vol.24 (1), p.1-15
Main Authors: Galindez-Ramirez, G., Carvalho, D. K. E., Lyra, P. R. M.
Format: Article
Language:English
Subjects:
Accuracy
Coefficients
Computer simulation
Earth and Environmental Science
Earth Sciences
Galerkin method
Geotechnical Engineering & Applied Earth Sciences
Hydrogeology
Mathematical Modeling and Industrial Mathematics
Mathematical models
Multiphase flow
Original Paper
Oscillations
Petroleum
Polynomials
Porous media
Reconstruction
Runge-Kutta method
Saturation
Soil Science & Conservation
Two phase flow
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:The focus of this work is to investigate and to apply, for the first time, a very high-resolution correction procedure via reconstruction (CPR) numerical discretization technique for the hyperbolic saturation equation that describes 1-D oil-water displacement through heterogeneous porous media. The CPR method can achieve high-order accuracy via a compact stencil consisting of the current cell and its immediate neighbors; in addition, the CPR recovers simplified versions of nodal discontinuous Galerkin (NDG), spectral volume (SV), and spectral difference (SD) methods using an adequate polynomial reconstruction function, whose coefficients are preprocessed and stored. Indeed the CPR versions of NDG and SV/SD are highly efficient. In order to suppress numerical oscillations near shocks, which are typical from higher-order schemes, a hierarchical MLP (multi-dimensional limiting process) is used in the reconstruction stage. The integration in time is carried out using a third-order RK (Runge-Kutta) method. A number of 1-D two-phase flow benchmark problems are solved, to verify the accuracy, efficiency, and shock-capturing capability of the adopted methodology.
ISSN:1420-0597
1573-1499
DOI:10.1007/s10596-019-09884-x