Element balance formulation in reactive compositional flow and transport with parameterization technique

Abstract

We present a novel nonlinear formulation for modeling reactive-compositional flow and transport in the presence of complex phase behavior due to a combination of thermodynamic and chemical equilibria in multi-phase systems. We apply this formulation to model precipitation/dissolution of minerals in reactive multiphase flow in subsurface reservoirs. The proposed formulation is based on the consistent element balance reduction of the molar (overall composition) formulation. To predict the complex phase behavior in such systems, we include the chemical equilibrium constraints to the multiphase multicomponent negative flash calculations and solve the thermodynamic phase and chemical phase equilibria simultaneously. In this solution, the phase equilibrium is represented by the partition coefficients, whereas the chemical equilibrium reaction is represented by the activity coefficients model. This provides a generic treatment of chemical and thermodynamic equilibrium within the successive substitution loop of multiphase flash to accommodate chemical equilibrium reactions (precipitation and dissolution reactions). Equilibrium Rate Annihilation matrix allows us to reduce the governing component conservation equations to element conservation equations, while the coupling between chemical and thermodynamic equilibrium is captured by a simultaneous solution of modified multiphase flash equations. The element balance equation written in terms of overall component mole fractions is modified and defined in terms of element mole fractions. Therefore, the primary set of governing equations are the element balance equations and the kinetic equations. This element composition of the mixture serves as an input to the modified multiphase flash computations, whereas the output is fractions of components in each phase, including solids. The nonlinear element-based governing equations are solved with the modified version of the operator-based linearization (OBL) approach where the governing equations are formulated in terms of space-and state-dependent parameters constrained by the solution of the extended multiphase flash. The element balance molar formulation along with the modified multiphase flash has been tested in a simple transport model with dissolution and precipitation reactions. The simulation of multidimensional problems of practical interest has been performed using the adaptive OBL technique. This is the first time when a robust multiphase multicomponent flash based on element fractions is coupled with an element balance-based compositional formulation and tested for multidimensional problems of practical interest. The proposed technique improves both robustness and performance of complex chemical models. simulation being the main tool in this process. Therefore, there are continuous efforts to improve the performance of the reservoir simulators and also implement a higher degree of physics into them in order to capture detailed subsurface processes with greater accuracy and without loss of computational efficiency. Reser...

An Advanced Discrete Fracture Methodology for Fast, Robust, and Accurate Simulation of Energy Production From Complex Fracture Networks

“…The proposed fracture network processing framework has been fully integrated with the open‐source Delft Advanced Research Terra Simulator (DARTS). DARTS is a scalable parallel simulation framework, which has been successfully applied for modeling of energy transition applications, including hydrocarbon (Khait & Voskov, 2018a; Lyu et al., 2021), geothermal (Khait & Voskov, 2018b; Wang et al., 2020) and CO 2 sequestration (Kala & Voskov, 2020; Lyu et al., 2021) cases. The ongoing effort to include fully coupled geomechanical modeling into DARTS allowed us to directly address induced seismicity problems (Novikov et al., 2021) which is another challenge in energy transition usually directly relevant to fracture networks.…”

Simulation of Foam-Assisted CO2 Storage in Saline Aquifers

“…The accuracy of this thermodynamic model vs. experimental results has been validated in Morshuis (2019). Foam-assisted CO 2 storage simulations for a brine aquifer are performed with the Delft Advanced Research Terra Simulator (DARTS) which is capable of modeling complex flow and transport related to various energy applications (Khait and Voskov, 2017;Kala and Voskov, 2020;Wang et al, 2020).…”

Simulation of Foam Enhanced-Oil-Recovery Processes Using Operator-Based Linearization Approach

“…To choose a proper resolution, a sensitivity study is required. A coarser parameterization space, corresponding to a smaller OBL resolution, can provide satisfying results, with higher efficiency compared with some academic or legacy simulators in previous tests Kala and Voskov 2020;Wang et al 2020). The proper implementation of the OBL approach can speed up Jacobian assembly by a factor of 14 and even faster on a graphics processing unit architecture (Khait et al 2020), which saves computational costs significantly.…”

“…The OBL approach is applied to improve the performance of numerical simulation. The Delft Advanced Research Terra Simulator, which is capable of modeling complex flow and transport related to various energy applications Kala and Voskov 2020;Wang et al 2020), is used to perform all simulation cases. First, the wet-foam model, as shown in Appendix A, is examined, where oil changes only the mobility of full-strength foam in the low-quality regime and has no direct effect on the high-quality regime.…”