A solution for multicomponent reactive transport under equilibrium and kinetic reactions

[1]Analysis of the space!time evolution of reactive solutes in porous systems is complex owing to the presence of different types of chemical reactions. The complete description of a reactive transport scenario entails calculating the spatial and temporal distribution of species concentrations and reaction rates. Here we develop an exact explicit expression for the space!time distribution of reaction rates for a scenario where the geochemical system can be described by an arbitrary number of equilibrium (fast) reactions and one kinetic (slow) reaction, in the absence of non!constant!activity immobile species. The key result is that the equilibrium reaction rate is the sum of two terms representing the availability of reactants. One term involves diffusion and dispersion and represents the contribution of mixing. The other term includes the contribution of the kinetic reaction. The approach also yields the local concentrations of all dissolved species. Yet the latter are not needed for the direct computation of equilibrium reaction rates. We illustrate the approach by means of a simple reactive transport scenario, involving a common ion effect in the presence of a kinetic and an equilibrium reaction leading to precipitation and dissolution processes within a one!dimensional fully saturated porous medium. The example highlights the highly nonlinear and nonmonotonic response of the system to the controlling input parameters. Citation:Sanchez!Vila, X., L. D. Donado, A. Guadagnini, and J. Carrera (2010), A solution for multicomponent reactive transport under equilibrium and kinetic reactions,Water Resour. Res.,46, W07539, doi:10.1029/2009WR008439.

[1]  C. Steefel,et al.  Reactive transport modeling: An essential tool and a new research approach for the Earth sciences , 2005 .

[2]  Jeanne M. VanBriesen,et al.  Multicomponent transport with coupled geochemical and microbiological reactions: model description and example simulations , 1998 .

[3]  C. Ayora,et al.  On the behavior of approaches to simulate reactive transport. , 2001, Journal of contaminant hydrology.

[4]  T. Skoulikidis,et al.  Inversion of Marble Sulfation. Reconversion of Gypsum Films into Calcite on the Surfaces of Monuments and Statues , 1984 .

[5]  Jacob Kielland.  Individual Activity Coefficients of Ions in Aqueous Solutions , 1937 .

[6]  Brian Berkowitz,et al.  Carbonate dissolution and precipitation in coastal environments: Laboratory analysis and theoretical consideration , 2004 .

[7]  D. Langmuir.  Aqueous Environmental Geochemistry , 1997 .

[8]  Jesús Carrera,et al.  Multicomponent reactive transport in multicontinuum media , 2009 .

[9]  J. Rubin,et al.  Transport of reacting solutes in porous media: Relation between mathematical nature of problem formulation and chemical nature of reactions , 1983 .

[10]  X. Sanchez‐Vila,et al.  An analytical approach to transient homovalent cation exchange problems , 2009 .

[11]  C. Steefel,et al.  Approaches to modeling of reactive transport in porous media , 1996 .

[12]  C. Steefel,et al.  A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution rea , 1994 .

[13]  A. Lasaga.  Chemical kinetics of water‐rock interactions , 1984 .

[14]  X. Sanchez‐Vila,et al.  Conditional Probability Density Functions of Concentrations for Mixing-Controlled Reactive Transport in Heterogeneous Aquifers , 2009 .

[15]  Brian Berkowitz,et al.  On Characterization of Anomalous Dispersion in Porous and Fractured Media , 1995 .

[16]  E. C. Childs.  Dynamics of fluids in Porous Media , 1973 .

[17]  Rainer Helmig,et al.  Numerical simulation of biodegradation controlled by transverse mixing , 1999 .

[18]  J. Carrera,et al.  Reactive transport modeling of calcite dissolution in the fresh-salt water mixing zone , 2005 .

[19]  Peter Knabner,et al.  A new numerical reduction scheme for fully coupled multicomponent transport‐reaction problems in porous media , 2005 .

[20]  Jesús Carrera,et al.  A mixing ratios‐based formulation for multicomponent reactive transport , 2007 .

[21]  Vladimir Cvetkovic,et al.  Transport of reactive tracers in rock fractures , 1999, Journal of Fluid Mechanics.

[22]  Peter Dietrich,et al.  Finiteness of steady state plumes , 2005 .

[23]  Marc Massot,et al.  Entropic structure of multicomponent reactive flows with partial equilibrium reduced chemistry , 2004 .

[24]  T. P. Clement,et al.  Modeling Multispecies Reactive Transport in Ground Water , 1998 .

[25]  S. Gorelick,et al.  Multiple‐Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore‐Scale Heterogeneity , 1995 .

[26]  Gour-Tsyh Yeh,et al.  A Model for Simulating Transport of Reactive Multispecies Components: Model Development and Demonstration , 1991 .

[27]  M. Dentz,et al.  Modeling non‐Fickian transport in geological formations as a continuous time random walk , 2006 .

[28]  W. Sanford,et al.  Simulation of calcite dissolution and porosity changes in saltwater mixing zones in coastal aquifers , 1989 .

[29]  X. Sanchez‐Vila,et al.  On matrix diffusion: formulations, solution methods and qualitative effects , 1998 .

[30]  Henning Prommer,et al.  Effects of hydrodynamic dispersion on plume lengths for instantaneous bimolecular reactions , 2004 .

[31]  Jesús Carrera,et al.  CHEPROO: A Fortran 90 object-oriented module to solve chemical processes in Earth Science models , 2009, Comput. Geosci..

[32]  D. Benson,et al.  Application of a fractional advection‐dispersion equation , 2000 .

[33]  R. Knapp.  Spatial and temporal scales of local equilibrium in dynamic fluid-rock systems , 1989 .

[34]  X. Sanchez‐Vila,et al.  Application of a mixing-ratios based formulation to model mixing-driven dissolution experiments , 2009 .

[35]  Alberto Guadagnini,et al.  A procedure for the solution of multicomponent reactive transport problems , 2005 .

[36]  Jesús Carrera,et al.  A mathematical formulation for reactive transport that eliminates mineral concentrations , 1998 .

[37]  Alberto Guadagnini,et al.  Reaction rates and effective parameters in stratified aquifers , 2008 .

[38]  Carl I. Steefel,et al.  Chapter 2. APPROACHES TO MODELING OF REACTIVE TRANSPORT IN POROUS MEDIA , 1996 .

[39]  C. Steefel,et al.  A new kinetic approach to modeling water-rock interaction: The role of nucleation, precursors, and Ostwald ripening , 1990 .

[40]  Jesús Carrera,et al.  A formulation for decoupling components in reactive transport problems , 2004 .

[41]  Kathryn L. Nagy,et al.  Chemical weathering rate laws and global geochemical cycles , 1994 .

[42]  W. Sanford,et al.  Ground-water control of evaporite deposition , 1990 .

[43]  Albert J. Valocchi,et al.  Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state , 2007 .

[44]  Wolfgang Kinzelbach,et al.  Numerical Modeling of Natural and Enhanced Denitrification Processes in Aquifers , 1991 .

[45]  Timothy R. Ginn,et al.  Fractional advection‐dispersion equation: A classical mass balance with convolution‐Fickian Flux , 2000 .

[46]  David W. Blowes,et al.  Multicomponent reactive transport modeling in variably saturated porous media using a generalized formulation for kinetically controlled reactions , 2002 .

[47]  M. Dentz,et al.  Transport‐controlled reaction rates under local non‐equilibrium conditions , 2007 .

[48]  Peter C. Lichtner,et al.  Continuum formulation of multicomponent-multiphase reactive transport , 1996 .