Development of a reactive transport solver in MATLAB Reservoir Simulation Toolbox using the fully-implicit sequential iterative approach
Abstract view|89|times PDF download|31|times
Abstract
Different operations dealing with the subsurface, such as subsurface CO2 disposal, hazardous waste disposal, geothermal energy extraction, underground hydrogen storage, etc., can change the fluid/flow system underground. The injection of fluids with thermodynamic and chemical properties different from those of the reservoir fluid can trigger a series of chemical reactions, which may affect the fluid and/or rock properties. Depending on the system under study, these changes may be advantageous or unfavorable. Reactive transport modeling is a choice for investigating how these changes can alter the system. In this study, a reactive transport solver is developed in the MATLAB Reservoir Simulation Toolbox using the sequential fully-implicit approach. The developed reactive transport solver is illustrated using reactions and geometries using reactions and geometries relevant for assessing the sealing capacity of a fractured caprock of a deep saline aquifer used for underground CO2 disposal, and the limitations and advantages of the approach are stated. Moreover, the results of the simulation for two fracture models, the discrete fracture matrix and embedded discrete fracture matrix models, are compared. The simulations demonstrate that hydrogen ion concentration or pH is the primary parameter affecting the extent of dissolution, while the other aqueous species concentrations are less influential. It is also shown that at higher flow rates, dissolution substantially occurs in the vicinity of the main fracture, along the flow direction, while at lower flow rates, because the injected fluid becomes fully buffered closer to the inlet, dissolution only occurs in the vicinity of the inlet over the course of the simulation. Applying the discrete fracture matrix and embedded discrete fracture matrix models to one of the scenarios demonstrates that both yield equivalent results.
Document Type: Original article
Cited as: Moslehi, S., Fazeli, H., Doster, F., Kord, S. Development of a reactive transport solver in MATLAB Reservoir Simulation Toolbox using the fully-implicit sequential iterative approach. Advances in Geo-Energy Research, 2025, 16(2): 114-130. https://doi.org/10.46690/ager.2025.05.04
Keywords
Full Text:
PDFReferences
Abd, A.S., Abushaikha, A.S. Reactive transport in porous media: A review of recent mathematical efforts in modeling geochemical reactions in petroleum subsurface reservoirs. SN Applied Sciences, 2021, 3: 401.
Ahusborde, E., El Ossmani, M. A sequential approach for numerical simulation of two-phase multicomponent flow with reactive transport in porous media. Mathematics and Computers in Simulation, 2017, 137: 71-89.
Ahusborde, E., El Ossmani, M., Moulay, M.I. A fully implicit finite volume scheme for single phase flow with reactive transport in porous media. Mathematics and Computers in Simulation, 2019, 164: 3-23.
Ahusborde, E., Kern, M., Vostrikov, V. Numerical simulation of two-phase multicomponent flow with reactive transport in porous media: Application to geological sequestration of CO2. ESAIM: Proceedings and Surveys, 2015, 50: 21-39.
Andrews, E., Hyman, J., Sweeney, M., et al. Fracture intensity impacts on reaction front propagation and mineral weathering in three-dimensional fractured media. Water Resources Research, 2023, 59(2): e2022WR032121.
Applegate, D., Appleyard, P. Capability for hydrogeochemical modelling within discrete fracture networks. Energies, 2022, 15: 6199.
Banshoya, S.I., Berre, I., Keilegavlen, E. Simulation of reactive transport in fractured porous media. Transport in Porous Media, 2023, 2: 643-667.
Barnes, D.A., Bacon, D.H., Kelley, S.R. Geological sequestration of carbon dioxide in the cambrian mount simon sandstone: Regional storage capacity, site characterization, and large-scale injection feasibility, michigan basin. Environmental Geosciences, 2009, 16: 163-183.
Bear, J. Dynamics of Fluids in Porous Media. New York, USA, Dover Publications, 2013.
Berre, I., Doster, F., Keilegavlen, E. Flow in fractured porous media: A review of conceptual models and discretization approaches. Transport in Porous Media, 2019, 130: 215-236.
Billings, M.P. Structural Geology. New Jersey, USA, PrenticeHall, 1972.
Cao, H. Development of techniques for general purpose simulators. Stanford, Stanford University, 2002.
Chen, Y., Ma, G., Wang, H. The simulation of thermo-hydro-chemical coupled heat extraction process in fractured geothermal reservoir. Applied Thermal Engineering, 2018, 143: 859-870.
Colin, M., Xavier, R., Halvor, N., et al. Robust chemical solver for fully-implicit simulations. ArXiv Preprint ArXiv:1806.03010, 2018.
de Hoop, S., Jones, E., Voskov, D. Accurate geothermal and chemical dissolution simulation using adaptive mesh refinement on generic unstructured grids. Advances in Water Resources, 2021, 154: 103977.
de Hoop, S., Voskov, D., Bertotti, G. Studying the effects of heterogeneity on dissolution processes using operator based linearization and high-resolution lidar data. Paper Presented at ECMOR XVII, Edinburgh, United Kingdom, 14-17 September, 2020.
Deng, H., Fitts, J.P., Crandall, D., et al. Alterations of fractures in carbonate rocks by CO2-acidified brines. Environmental Science & Technology, 2015, 49: 10226-10234.
Elgendy, A., Pizzolato, A., Maniglio, M., et al. Reactive transport modelling of H2 storage in depleted gas fields: An approach to implement biogeochemical reactions in a compositional reservoir simulator. Paper SPE 214434 Presented at the SPE EuropEC-Europe Energy Conference featured at the 84th EAGE Annual Conference & Exhibition, Vienna, Austria, 5-8 June, 2023.
Fahs, M., Carrayrou, J., Younes, A., et al. On the efficiency of the direct substitution approach for reactive transport problems in porous media. Water, Air, and Soil Pollution, 2008, 193: 299-308.
Fan, Y., Durlofsky, L.J., Tchelepi, H.A. Numerical simulation of the in-situ upgrading of oil shale. SPE Journal, 2010, 15: 368-381.
Fan, Y., Durlofsky, L.J., Tchelepi, H.A. A fully-coupled flow-reactive-transport formulation based on element conservation, with application to CO2 storage simulations. Advances in Water Resources, 2012, 42: 47-61.
Fang, Y., Yeh, G.T., Burgos, W.D. A general paradigm to model reaction-based biogeochemical processes in batch systems. Water Resources Research, 2003, 39(4): 1083.
Fitts, J.P., Peters, C.A. Caprock fracture dissolution and CO2 leakage. Reviews in Mineralogy and Geochemistry, 2013, 77: 459-479.
Franc, J., Møyner, O., Tchelepi, H.A. Coupling-strength criteria for sequential implicit formulations. Journal of Computational Physics, 2023, 492: 112413.
Gholami, R. Hydrogen storage in geological porous media: Solubility, mineral trapping, H2S generation and salt precipitation. Journal of Energy Storage, 2023, 59: 106576.
Gilmore, T., Bonneville, A., Sullivan, C., et al. Characterization and design of the FutureGen 2.0 carbon storage site. International Journal of Greenhouse Gas Control, 2016, 53: 1-10.
Gong, B., Zhang, Y., Fan, Y., et al. A novel approach to model enhanced coal bed methane recovery with discrete fracture characterizations in a geochemical simulator. Journal of Petroleum Science and Engineering, 2014, 124: 198-208.
Gylling, B., Trinchero, P., Molinero, J., et al. A DFN-based high performance computing approach to the simulation of radionuclide transport in mineralogically heterogeneous fractured rocks. Paper H51C-1479 Presented at the AGU Fall Meeting Abstracts, San Francisco, California, 12-16 December, 2016.
Hammond, G., Lichtner, P., Lu, C. Subsurface multiphase flow and multicomponent reactive transport modeling using high-performance computing. Journal of Physics: Conference Series, 2007, 78(1): 012025.
Hammond, G.E., Valocchi, A.J., Lichtner, P.C. Application of Jacobian-free Newton-Krylov with physics-based preconditioning to biogeochemical transport. Advances in Water Resources, 2005, 28: 359-376.
Hemme, C., Van Berk, W. Hydrogeochemical modeling to identify potential risks of underground hydrogen storage in depleted gas fields. Applied Sciences, 2018, 8: 2282.
Hu, M., Rutqvist, J. Multi-scale coupled processes modeling of fractures as porous, interfacial and granular systems from rock images with the numerical manifold method. Rock Mechanics and Rock Engineering, 2022, 55: 3041-3059.
Jiang, C., Wang, X., Pu, S., et al. Incipient karst generation in jointed layered carbonates: Insights from three-dimensional hydro-chemical simulations. Journal of Hydrology, 2022, 610: 127831.
Jiang, C., Wang, X., Pu, S., et al. Karst genesis and wormhole formation in carbonate joint networks: A comparison between 3D and 2D modeling. Journal of Hydrology, 2023, 619: 129303.
Jiang, Y. Techniques for modeling complex reservoirs and advanced wells. Stanford, Stanford University, 2007.
Kampman, N., Bertier, P., Busch, A., et al. Validating reactive transport models of CO2-brine-rock reactions in caprocks using observations from a natural CO2 reservoir. Energy Procedia, 2017, 114: 4902-4916.
Karimi-Fard, M., Durlofsky, L.J., Aziz, K. An efficient discrete-fracture model applicable for general-purpose reservoir simulators. SPE Journal, 2004, 9: 227-236.
Keller, S.J. Analyses of subsurface brines of Indiana. Bloomington, Indiana, US Geological & Water Survey, 1983.
Krogstad, S., Lie, K.A., Møyner, O., et al. MRST-AD an open-source framework for rapid prototyping and evaluation of reservoir simulation problems. Paper SPE 173317 Presented at SPE Reservoir Simulation Conference, Houston Texas, 23-25 February, 2015.
Labotka, D.M., Panno, S.V., Locke, R.A., et al. Isotopic and geochemical characterization of fossil brines of the Cambrian Mt. Simon sandstone and Ironton-Galesville formation from the Illinois Basin, USA. Geochimica et Cosmochimica Acta, 2015, 165: 342-360.
Lasaga, A.C. Chemical kinetics of water-rock interactions. Journal of Geophysical Research: Solid Earth, 1984, 89: 4009-4025.
Lee, K.J., Moridis, G.J., Ehlig-Economides, C.A. Compositional simulation of hydrocarbon recovery from oil shale reservoirs with diverse initial saturations of fluid phases by various thermal processes. Energy Exploration & Exploitation, 2017, 35: 172-193.
Li, J., Tomin, P., Tchelepi, H. Sequential fully implicit newton method for flow and transport with natural black-oil formulation. Computational Geosciences, 2023, 27: 485-498.
Li, S., Kang, Z., Feng, X., et al. Three-dimensional hydrochemical model for dissolutional growth of fractures in karst aquifers. Water Resources Research, 2020, 56: e2019WR025631.
Lie, K.A. An Introduction to Reservoir Simulation Using MATLAB/GNU Octave: User Guide for the MATLAB Reservoir Simulation Toolbox (MRST). Cambridge, UK, Cambridge University Press, 2019.
Lie, K.A., Moyner, O. Advanced Modelling with the MATLAB Reservoir Simulation Toolbox. Cambridge, UK, Cambridge University Press, 2021.
Liu, P., Yao, J., Couples, G.D., et al. Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks. Chemical Engineering Science, 2017, 172: 143-157.
Liu, X., Ormond, A., Bartko, K., et al. A geochemical reaction-transport simulator for matrix acidizing analysis and design. Journal of Petroleum Science and Engineering, 1997, 17: 181-196.
Luhmann, A.J., Kong, X., Tutolo, B.M., et al. Experimental dissolution of dolomite by CO2-charged brine at 100 °C and 150 bar: Evolution of porosity, permeability, and reactive surface area. Chemical Geology, 2014, 380: 145-160.
Machado, M.V.B., Delshad, M., Sepehrnoori, K. Modeling self-sealing sechanisms in fractured carbonates induced by CO2 injection in saline aquifers. ACS Omega, 2023, 8: 48925.
Mayer, K.U., Frind, E.O., Blowes, D.W. Multicomponent reactive transport modeling in variably saturated porous media using a generalized formulation for kinetically controlled reactions. Water Resources Research, 2002, 38: 13-1-13-21.
Mehnert, E., Damico, J.R., Grigsby, N.P., et al. Geologic carbon sequestration in the Illinois basin: Numerical modeling to evaluate potential impacts. Champaign, Illinois, Illinois State Geological Survey, 2019.
Moinfar, A., Varavei, A., Sepehrnoori, K., et al. Development of an efficient embedded discrete fracture model for 3D compositional reservoir simulation in fractured reservoirs. SPE Journal, 2014, 19: 289-303.
Molson, J., Aubertin, M., Bussière, B. Reactive transport modelling of acid mine drainage within discretely fractured porous media: Plume evolution from a surface source zone. Environmental Modelling & Software, 2012, 38: 259-270.
Morel, F.M., Hering, J.G. Principles and Applications of Aquatic Chemistry. New York, USA, John Wiley & Sons, 1993.
Nghiem, L., Sammon, P., Grabenstetter, J., et al. Modeling CO2 storage in aquifers with a fully-coupled geochemical EOS compositional simulator. Paper SPE 89474 Presented at SPE Improved Oil Recovery Conference, Tulsa, Oklahoma, 17-21 April, 2004.
Nghiemw, L., Shrivastavar, A., Kohse, B., et al. Simulation of CO2 EOR and sequestration processes with a geochemical EOS compositional simulator. Paper PETSOC2004-051 Presented at Canadian International Petroleum Conference, Calgary, Alberta, 8-10 June, 2004.
Palandri, J.L., Kharaka, Y.K. A compilation of rate parameters of water-mineral interaction kinetics for application to geochemical modeling. Reston, Viriginia, US Geological Survey, 2004.
Palkovic, M.J. Depositional characterization of the Eau Claire formation at the Illinois Basin-Decatur Project: Facies, mineralogy and geochemistry. Urbana, Illinois, University of Illinois at Urbana-Champaign, 2015.
Parkhurst, D.L., Appelo, C. User’s guide to PHREEQC (Version 2): A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Denver, Colorado, US Geological Survey, 1999.
Pruess, K. TOUGH2-A general-purpose numerical simulator for multiphase fluid and heat flow. Berkeley, California, Lawrence Berkeley Laboratory, 1991.
Pruess, K., Garcia, J., Kovscek, T., et al. Intercomparison of numerical simulation codes for geologic disposal of CO2. Berkeley, California, Lawrence Berkeley National Laboratory, 2002.
Saaltink, M.W., Carrera, J., Ayora, C. On the behavior of approaches to simulate reactive transport. Journal of Contaminant Hydrology, 2001, 48: 213-235.
Salih, H.H., Dastgheib, S.A. Treatment of a hypersaline brine, extracted from a potential CO2 sequestration site, and an industrial wastewater by membrane distillation and forward osmosis. Chemical Engineering Journal, 2017, 325: 415-423.
Seigneur, N., Socié, A., Mayer, K.U. A compositional global implicit approach for modeling coupled multicomponent reactive transport. Water Resources Research, 2023, 59: e2021WR031774.
Shaik, R., Tomin, P., Voskov, D. Modeling of near-well matrix acidization. Paper Presented at 43rd Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, 12-14 February, 2018.
Sherman, T., Sole-Mari, G., Hyman, J., et al. Characterizing reactive transport behavior in a three-dimensional discrete fracture network. Transport in Porous Media, 2023, 146: 307-327.
Shewchuk, J.R. Triangle: Engineering a 2D auality mesh generator and Delaunay triangulator, Workshop on applied computational geometry, Paper Presented at Applied Computational Geometry Towards Geometric Engineering, Philadelphia, PA, 27-28 May, 1996.
Song, G., Song, X., Ji, J., et al. Evolution of fracture aperture and thermal productivity influenced by chemical reaction in enhanced geothermal system. Renewable Energy, 2022a, 186: 126-142.
Song, G., Song, X., Xu, F., et al. Numerical parametric investigation of thermal extraction from the enhanced geothermal system based on the thermal-hydraulic-chemical coupling model. Journal of Cleaner Production, 2022b, 352: 131609.
Steefel, C.I., Hu, M. Reactive transport modeling of mineral precipitation and carbon trapping in discrete fracture networks. Water Resources Research, 2022, 58:e2022WR032321.
Steefel, C.I., Lasaga, A.C. A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution reactions with application to reactive flow in single phase hydrothermal systems. American Journal of Science, 1994, 294:529-592.
Stein, E., Frederick, J.M., Hammond, G.E., et al. Modeling Coupled Reactive Flow Processes in Fractured Crystalline Rock. Albuquerque, New Mexico, Sandia National Laboratory, 2017.
Thibeau, S., Nghiem, L.X., Ohkuma, H. A modeling study of the role of selected minerals in enhancing CO2 mineralization during CO2 aquifer storage. Paper SPE 109739 Presented at SPE Annual Technical Conference and Exhibition Anaheim, California, 11-14 November, 2007.
Vu, P.T., Ni, C., Li, W., et al. Particle-based workflow for modeling uncertainty of reactive transport in 3D discrete fracture networks. Water, 2019a, 11(12):2502.
Vu, P.T., Ni, C., Li, W., et al. Particle tracking approach to model chemical reaction transport in 3D discrete fracture networks. Geophysical Research Abstracts, 2019b, 21:EGU2019-6678.
Wang, T., Sun, Z., Sun, H., et al. Development of a scalable parallel compositional simulator for thermo-hydromechanical coupling in fractured rocks using an embedded discrete fracture model. SPE Journal, 2024, 29(5):2545-2565.
Xu, T., Pruess, K. Modeling multiphase non-isothermal fluid flow and reactive geochemical transport in variably saturated fractured rocks: 1. Methodology. American Journal of Science, 2001, 301:16-33.
Xu, T., Sonnenthal, E., Spycher, N., et al. TOUGHREACT user’s guide: A simulation program for non-isothermal multiphase reactive geochemical transport in variable saturated geologic media. Berkeley, USA, Lawrence Berkeley National Laboratory, 2004.
Yan, B., Mi, L., Wang, Y., et al. Multi-porosity multi-physics compositional simulation for gas storage and transport in highly heterogeneous shales. Journal of Petroleum Science and Engineering, 2018, 160:498-509.
Yeh, G., Tripathi, V. A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components. Water Resources Research, 1989, 25:93-108.
Zhang, W., Han, D., Wang, B., et al. Thermal-hydraulic-mechanical-chemical modeling and simulation of an enhanced geothermal system based on the framework of extended finite element methods-embedded discrete fracture model. Journal of Cleaner Production, 2023, 415: 137630.
DOI: https://doi.org/10.46690/ager.2025.05.04
Refbacks
- There are currently no refbacks.
Copyright (c) 2025 The Author(s)
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.