A Numerical Model of Tracer Transport in a Non-isothermal Two-Phase Flow System for $${\text{ CO}}_2$$ Geological Storage Characterization

Abstract

For the purpose of characterizing geologically stored CO 2 including its phase partitioning and migration in deep saline formations, different types of tracers are being developed. Such tracers can be injected with CO 2 or water, and their partitioning and/or reactive transfer from one phase to another can give information on the interactions between the two fluid phases and the development of their interfacial area. Kinetic rock-water interactions and geochemical reactions during two-phase flow of CO 2 and brine have been incorporated in numerical simulators (e.g., Xu et al., TOUGHREACT User's Guide: A Simulation Program for Non-isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geologic Media. LBNL Report 55460, V.1.2., Berkeley, CA, 2004). However, chemical equilibrium between the fluid phases is typically assumed, and multi-component, multiphase, non-isothermal codes for CO 2 -brine systems that incorporate kinetic mass transfer of tracers between the two fluid phases are not readily available. New models or further developments of existing models are therefore needed to provide the capability for interpreting the signals of novel tracers, including tracers with kinetic/time-dependent interface transfer. This paper presents such new numerical model of tracer transport in a non-isothermal two-phase flow system. The model consists of five different governing equations describing liquid phase (aqueous) flow, gas (CO 2 ) flow, heat transport and the movement of the tracers within the two phases, as well as allowing kinetic transport of the tracers between the two phases. A finite element method is adopted for the spatial discretization and a finite difference approach is used for temporal discretization. Some special technologies and solution strategies are adopted for increasing the convergence, ensuring the numerical stability and eliminating non-physical oscillations. The new numerical model is validated against the code TOUGH2/ECO2N as well as some analytical/semi-analytical solutions. Good agreement between the simulated and analytical results indicates that the model has capability to simulate two-phase flow and

An existence result for nonisothermal immiscible incompressible 2‐phase flow in heterogeneous porous media

“…In this work, a coupled reservoir 2‐phase flow model is described, which accounts for varying reservoir temperature to capture flow physics accurately. Although considerable progress has been made in the computational simulation of 2‐phase problems under nonisothermal conditions (see, eg, some related works() and the refrences therein), to the best knowledge of the authors, the mathematical analysis of such coupled models under nonisothermal conditions is still incomplete. The previous results on the existence of a weak solution of a simplified system describing nonisothermal 2‐phase flow in porous media were obtained in Bocharov and Monakhov() and then revisited in Monakhov .…”

Effect of the Water-Air Coupling on the Stability of Rainfall-Induced Landslides Using a Coupled Infiltration and Hydromechanical Model

“…The water saturation and pore-air pressure are set as the basic variables, and the solution of Equations ( 3) and ( 4) is replaced by that of Equations ( 9) and ( 4). The iterative method for solving the two-phase flow equations was proven to have good performance of convergence and numerical stability [33].…”

Modelling the Unsteady Flow of Water into a Partly Saturated Soil

“…[28], who used the well-known pseudo-transient continuation approach to solve the nonlinear transient water infiltration problem, The authors recommended stability diagrams for the exact solution of the Richards equation. In geo-environment applications, Bunsri et al [29] solved the Richards's equation with advection-dispersion and solute transport equations by the Galerkin technique.…”