A Finite Volume Framework for the Fully Implicit Thermal-Hydro-Mechanical-Compositional Modeling in Geo-Energy Applications

Abstract

The role of Thermal-Hydro-Mechanical-Compositional analysis in the development of geo-energy resources has been amplified in recent years. As an example, challenges such as wellbore stability, land subsidence and induced seismicity highlight the necessity for comprehensive geomechanical evaluations which are then coupled with thermo-hydrodynamical processes within the reservoir. Numerical simulations of the coupled thermo-poromechanical processes provide a general-purpose tool capable of performing these evaluations at both continuum laboratory and field scales. However, efficient integration of the coupled system of fluid mass, energy and momentum conservation equations poses multiple numerical and implementation difficulties, such as combining different numerical methods on staggered grids and associated limitations on admissible grids. This paper introduces a new fully-implicit scheme of the Finite Volume Method (FVM) for modeling thermal compositional flow in thermo-poroelastic rocks. The scheme uses the gradient-based variant of coupled multi-point approximations of fluid mass, momentum, heat convection and conduction fluxes, which are derived from their respective local balances. The novelty of the scheme is that it incorporates temperature into the approximation of these fluxes. Consequently, the approximation of displacement gradients depends on temperatures, while the approximation of temperature itself is derived from the balance of heat conduction fluxes. At the same time, we utilize a single-point upstream weighting for the temperature-dependent terms in heat convection fluxes. The resulting scheme respects the local balance of fluxes in the presence of temperature gradients. Besides, it also supports star-shaped and various boundary conditions. Overall, the scheme represents a unified FVM-based approach for the integration of all conservation laws relevant to geo-energy applications on a cell-centered collocated grid. Furthermore, the implemented two-stage block-partitioned preconditioning strategy enables the efficient solution of obtained linear systems. The proposed modeling framework has been implemented in an open-source Delft Advanced Research Terra Simulator (DARTS). Moreover, the flexibility regarding compositional fluid properties is reinforced by the Operator-Based Linearization (OBL) technique incorporated into DARTS. The proposed modeling framework has undergone rigorous validation in convergence study, and comparisons against established analytical and numerical solutions. The framework covers advanced physical phenomena including thermal expansion and contraction, porosity dependent on pressure, temperature and strain, and multiphase flow with phase changes and chemical alterations. The framework capabilities and the performance of the preconditioning strategy have been assessed in the mechanical extension of the 10th SPE Comparative study (SPE10) model.