Multi-Stage Preconditioners for Thermal-Compositional-Reactive Flow in Porous Media

by   Matthias A. Cremon, et al.

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe convergence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Peclet number) to diffusion dominated (low Peclet number). The number of linear iterations is reduced by 40-85 to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.


A constrained pressure-temperature residual (CPTR) method for non-isothermal multiphase flow in porous media

In petroleum reservoir simulation, the standard preconditioner, the Cons...

A block preconditioner for non-isothermal flow in porous media

In petroleum reservoir simulation, the industry standard preconditioner,...

Constrained Pressure-Temperature Residual (CPTR) Preconditioner Performance for Large-Scale Thermal CO2 Injection Simulation

This work studies the performance of a novel preconditioner, designed fo...

An Upwind Generalized Finite Difference Method (GFDM) for Analysis of Heat and Mass Transfer in Porous Media

In this paper, an upwind meshless GFDM is presented for the first time a...

An Improved Multi-Stage Preconditioner on GPUs for Compositional Reservoir Simulation

The compositional model is often used to describe multicomponent multiph...

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

This paper focuses on developing a reduction-based algebraic multigrid m...

A boundary-oriented reduced Schwarz domain decomposition technique for parametric advection-diffusion problems

We present in this paper the results of a research motivated by the need...

Please sign up or login with your details

Forgot password? Click here to reset