Reactive transport in porous media: a review of recent mathematical efforts in modeling geochemical reactions in petroleum subsurface reservoirs

Abstract

The rapid advancements in the computational abilities of numerical simulations have attracted researchers to work on the area of reactive transport in porous media to improve the hydrocarbon production processes from mature reservoirs. In the hydrology community, reactive transport is well developed where the main research focuses on studying the movement of groundwater and contaminants in aquifers, and quantifying the effect of chemical reactions between the rocks and water. Recently, great efforts have been made to adapt similar models for petroleum applications where multiphase flow is experienced in the subsurface reservoirs. In such systems, thermodynamic and chemical equilibrium is key in establishing an accurate description of the states of the fluids existing in the reservoir. This paper presents a detailed and comprehensive review on the concepts of geochemical modeling, and how it can be mathematically adapted to petroleum multiphase flow problems in porous media. We introduce key physical concepts outlining the treatment of chemical reactions in experimental trials and then explain in detail how a network of chemical reactions can be modeled mathematically for numerical simulation applications. The steps of characterizing the physical behavior of the fluid flow—through a set of governing equations by either natural or molar variables formulations, and the methodology to simplify and incorporate the numerical algorithms into an existing reservoir simulation scheme are shown as well. We finally present two numerical cases from the literature to highlight the key variations between the different variable formulations and comment on the advantages and disadvantages of each approach.

Carbonated Water Injection for Enhanced Oil Recovery and CO2 Geosequestration in Different CO2 Repositories: A Review of Physicochemical Processes and Recent Advances

“…A schematic illustration of the pore and the reservoir scale levels of CO 2 mineral dissolution and precipitation (A) with pressure and (B) without pressure. Reproduced with permission from ref . Copyright 2021 Springer Nature.…”

Coupling of Rigorous Multiphase Flash with Advanced Linearization Schemes for Accurate Compositional Simulation

“…To achieve thermodynamic equilibrium, the fugacity of each component in each phase should be equal. A generic equation could represent a three-phase system at equilibrium by, (1) where f i is the fugacity of the vapor (V), liquid (L), or aqueous (Aq) phase, and i is the index of the component. The equation indicates that when three phases are present, their respective composition will keep changing until thermodynamic equilibrium is achieved.…”

An emulation-based approach for interrogating reactive transport models

“…Unfortunately, due to the computational expense of many modern multi-component RTMs (e.g. Abd and Abushaikha, 2021;Seigneur et al, 2021;Gharasoo et al, 2022), it is normally impractical to perform sensitivity analyses in more than a few dimensions and it is up to the investigator to use their knowledge of the system to choose which sensitivity analyses are necessary to explore (Steefel et al, 2005b). Ideally, we would be able to systematically perform sensitivity analyses over many model parameters, considering how model outputs vary as a function of multiple input parameters simultaneously (i.e.…”