Reactive transport modeling in porous media

Reactive transport modeling in porous media refers to the creation of computer models integrating chemical reaction with transport of fluids through the Earth's crust. Such models predict the distribution in space and time of the chemical reactions that occur along a flowpath. Reactive transport modeling in general can refer to many other processes, including reactive flow of chemicals through tanks, reactors, or membranes; particles and species in the atmosphere; gases exiting a smokestack; and migrating magma.

Overview

Reactive transport models are constructed to understand the composition of natural waters; the origin of economic mineral deposits; the formation and dissolution of rocks and minerals in geologic formations in response to injection of industrial wastes, steam, or carbon dioxide; and the generation of acidic waters and leaching of metals from mine wastes. They are often relied upon to predict the migration of contaminant plumes; the mobility of radionuclides in waste repositories; and the biodegradation of chemicals in landfills. When applied to the study of contaminants in the environments, they are known as fate and transport models. ^[1]

Development of reactive transport modeling

Modern reactive transport modeling has arisen from several separate schools of thought. ^[2] Hydrologists primarily concerned with the physical nature of mass transport assumed relatively simple reaction formulations, such as linear distribution coefficients or linear decay terms, which could be added to the advection-dispersion equation. By assuming linear, equilibrium sorption, for example, the advection-dispersion equation can be modified by a simple retardation factor and solved analytically. Such analytical solutions are limited to relatively simple flow systems and reactions.

Geochemical models, on the other hand, have been developed to provide thermodynamic descriptions of multicomponent systems without regard to transport. Reaction path models were created, for instance, to describe the sequence of chemical reactions resulting from chemical weathering or hydrothermal alteration in batch systems, in terms of the overall reaction progress. By adopting the reference frame of a packet of fluid and treating reaction progress as travel time (or distance along a flowpath), however, a batch reaction path model could be thought of as describing advective transport through an aquifer. ^[3]

The most sophisticated multi-component reactive transport models consider both reaction and transport. ^{[4] [5]} Early studies developed the theoretical basis of reactive transport models, and the numerical tools necessary to solve them, and applied them to problems of reactive contaminant transport ^[6] and flow through reacting hydrothermal systems. ^[7]

Reactive transport models have found increased application in recent years with improvements in the power of personal computers and modeling software. ^{[5] [8]}

Processes considered in reactive transport models

Reactive transport models couple a large number chemical reactions with mass transport. Certain applications, such as geothermal energy production and ore deposit modeling, require the additional calculation of heat transfer. In modeling carbon sequestration and hydraulic fracturing, moreover, it may be necessary to describe rock deformation resulting from mineral growth or abnormally high fluid pressure. Description of transport through the unsaturated zone and multiphase flow modeling, as applied to transport of petroleum and natural gas; non-aqueous phase liquids (DNAPL or LNAPL); and supercritical carbon dioxide requires increasingly complex models which are prone to considerable uncertainty.

In many cases the processes simulated in reactive transport models are highly related. Mineral dissolution and precipitation, for example, can affect the porosity and permeability of the domain, which in turn affect the flow field and groundwater velocity. Heat transport greatly affects the viscosity of water and its ability to flow. Below are many of the physical and chemical processes which can be simulated with reactive transport models.

Geochemical reactions:

• Acid-base reactions
• Aqueous complexation
• Mineral dissolution and precipitation
• Reduction and oxidation (redox) reactions, including those catalyzed by enzymes, surfaces, and microorganisms
• Sorption, ion exchange, and surface complexation
• Gas dissolution and exsolution
• Stable isotope fractionation
• Radioactive decay

Mass Transport:

• Advection
• Molecular scale diffusion
• Hydrodynamic dispersion
• Colloid-facilitated transport

Heat transport:

• Advection
• Conduction
• Convection

Medium deformation:

• Compression or expansion of the domain
• Fracture formation

Solving reactive transport models

Some of the simplest reactive transport problems can be solved analytically. Where equilibrium sorption is described by a linear distribution coefficient, for example, the sorbing solute's velocity is retarded relative to that of a nonreactive tracer; the relative velocities can be described with a retardation factor. Analytical solutions are exact solutions of the governing equations.

Complex reactive transport problems are more commonly solved numerically. In this case, the governing equations are approximated so that they can be solved by computer algorithms. The governing equations, including both reaction and transport terms, can be solved simultaneously using a one-step or global implicit simulator. This technique is straightforward conceptually, but computationally very difficult. ^[9]

Instead of solving all the relevant equations together, the transport and chemical reaction equations can be solved separately. Operator splitting, as this technique is known, uses appropriate numerical techniques to solve the reaction and transport equations at each time step. ^[1] Various methods exist, including the sequential non-iterative approach (SNIA), Strang splitting, and sequential iterative approach (SIA). ^[9] Since the reaction and transport terms are handled separately, separate programs for batch reaction and transport can be linked together. Cross-linkable re-entrant software objects designed for this purpose readily enable construction of reactive transport models of any flow configuration. ^{[10] [11]}

Challenges

Reactive transport modeling requires input from numerous fields, including hydrology, geochemistry and biogeochemistry, microbiology, soil physics, and fluid dynamics. ^[2] The numerical formulation and solution of reactive transport problems can be especially difficult due to errors arising in the coupling process, beyond those inherent to the individual processes. Valocchi and Malmstead (1992), for example, reported on the potential errors arising from the operator splitting technique. ^[12]

Even in the absence of numerical difficulties, the general lack of knowledge available to practitioners creates uncertainty. Field sites are typically heterogeneous, both physically and chemically, and sampling is often sparse. The prevailing assumption of Fickian dispersion is often inadequate. Equilibrium constants and kinetic rate laws for relevant reactions are often poorly known. The complexity of many processes requires expertise in one or more of the aforementioned fields. Many processes, such as long-term nuclear waste storage, cannot be experimentally verified; reactive transport problems can only attempt to predict such long-term behavior. The current descriptions of multi-phase flow and mechanical deformation processes are still being developed.

Software programs in common use

• PFLOTRAN*
• ChemPlugin
• MIN3P
• CHESS, ^[13] HYTEC
• CrunchFlow ^{[7] [14]}
• The Geochemist's Workbench ^{[1] [15]}
• HYDROGEOCHEM ^[16]
• THMC
• PHREEQC, ^{[17] [18]}
• PHAST ^[19]
• Reaktoro ^[20]
• TOUGHREACT ^[21]
• OpenGeoSys
• PHT3D
• PNBRNS ^[5]
• HP1 / HP2

See also

• Chemical thermodynamics
• Chemical kinetics
• Geochemistry
• Geomicrobiology
• Hydrogeology
• Groundwater model
• Geochemical modeling
• Reservoir simulation
• Chemical process modeling
• Chemical transport model

References

1. Bethke, C.M., 2008, Geochemical and Biogeochemical Reaction Modeling. Cambridge University Press, 547 pp.
2. Steefel, C.I., D.J. DePaolo, and P.C. Lichtner, 2005, Reactive transport modeling: An essential tool and a new research approach for the Earth sciences. Earth and Planetary Science Letters240, 539-558
3. Lichtner, P.C., 1988, The quasi-stationary state approximation to coupled mass transport and fluid-rock interaction in a porous medium. Geochemica et Cosmochemica Acta52, 143-165
4. Lichtner, P.C., 1985, Contiuum model for simultaneous chemical reactions and mass transport in hydrothermal systems. Geochemica et Cosmochemica Acta49, 779-800
5. Gharasoo, M., Centler, F., Regnier, P., Harms, H., Thullner, M., 2012. A reactive transport modelling approach to simulate biogeochemical processes in pore structures with pore-scale heterogeneities. Environmental Modelling & Software 30, 102-114.
6. Yeh, G.T. and V.S. Tripathi, 1989, A critical evaluation of recent developments in hydrogeochemical transport models of reactive multi-chemical components. Water Resources Research25, 93-108
7. Steefel, C.I. and A.C. Lasaga, 1994, 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 Science294, 529-592
8. Brady, P.V. and C.M. Bethke, 2000, Beyond the Kd approach. Ground Water38, 321-322
9. Steefel, C.I. and K.T.B. MacQuarrie, 1996, Approaches to modeling of reactive transport in porous media. In P.C. Lichtner, C.I. Steefel, and E.H. Oelkers (eds.), 1996, Reactive Transport in Porous Media. Reviews in Mineralogy34, 85-129
10. "Aqueous Solutions launches ChemPlugin". Researchpark.illinois.edu. Research Park at the University of Illinois. Retrieved 7 June 2013.
11. "ChemPlugin.ORG". Aqueous Solutions LLC. Retrieved 3 May 2013.
12. Valocchi, A.J. and M. Malmstead, 1992, Accuracy of operator splitting for advection-dispersion-reaction problems. Water Resources Research28(5), 1471-1476
13. van der Lee, J., and L. De Windt, 2000, CHESS, another speciation and complexation computer code. Technical Report no. LHM/RD/93/39, Ecole des Mines de Paris, Fontainebleau
14. Steefel, C.I., 2001, GIMRT, Version 1.2: Software for modeling multicomponent, multidimensional reactive transport, User's Guide. Report UCRL-MA-143182, Lawrence Livermore National Laboratory, Livermore, California.
15. Bethke, C.M., B. Farrell, and M. Sharifi, 2021, The Geochemist's Workbench Release 15 (five volumes). Aqueous Solutions LLC, Champaign, IL
16. Cheng, H.P. and G.T. Yeh, 1998, Development of a three-dimensional model of subsurface flow, heat transfer, and reactive chemical transport: 3DHYDROGEOCHEM. Journal of Contaminant Hydrology34, 47-83
17. Parkhurst, D.L., 1995, User's Guide to PHREEQC, a computer model for speciation, reaction-path, advective-transport and inverse geochemical calculations. US Geological Survey Water-Resources Investigations Report 95-4227.
18. Parkhurst, D.L. and C.A.J. Appelo, 1999, User's Guide to PHREEQC (version 2), a computer program for speciation, batch-reaction, one-dimensional transport and inverse geochemical calculations. US Geological Survey Water-Resources Investigations Report 99-4259.
19. Parkhurst, D.L., Kipp, K.L., and Charlton, S.R., 2010, PHAST Version 2—A program for simulating groundwater flow, solute transport, and multicomponent geochemical reactions: U.S. Geological Survey Techniques and Methods 6–A35, 235 p.
20. Leal, A. M. M. et al. (2018). Reaktoro, a unified, open-source framework for modeling chemically reactive systems, https://github.com/reaktoro/reaktoro
21. Xu, T., E.L. Sonnenthal, N. Spycher and K. Pruess, 2004, TOUGHREACT user's guide: A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media. Report LBNL-55460, Lawrence Berkeley National Laboratory, Berkeley, California.

Further reading

• Appelo, C.A.J. and D. Postma, 2005, Geochemistry, Groundwater, and Pollution. Taylor & Francis, 683 pp. ISBN   978-0415364287
• Bethke, C.M., 2008, Geochemical and Biogeochemical Reaction Modeling. Cambridge University Press, 547 pp. ISBN   978-0521875547
• Lichtner, P.C., C.I. Steefel, and E.H. Oelkers (eds.), 1996, Reactive Transport in Porous Media. Reviews in Mineralogy34, 438 pp. ISBN   978-0-939950-42-3
• Merkel, B.J., B. Planer-Friedrich, and D.K. Nordstrom, 2008, Groundwater Geochemistry: A Practical Guide to Modeling of Natural and Contaminated Aquatic Systems. Springer, 242 pp. ISBN   978-3540746676
• Zhang, F., G.T. Yeh, and J.C. Parker (eds.), 2012, Groundwater Reactive Transport Models. Behtham Publishers, 254 pp. ISBN   9781608053063
• Zhu, C. and G. Anderson, 2002, Environmental Applications of Geochemical Modeling. Cambridge University Press, 300 pp. ISBN   978-0521005777