Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

[1]  Joachim Hoffmann,et al.  Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems—the MoMaS benchmark case , 2010 .

[2]  J. Allison,et al.  MINTEQA2/PRODEFA2, a geochemical assessment model for environmental systems: Version 3. 0 user's manual , 1991 .

[3]  P. Gouze,et al.  Chemical reactions and porosity changes during sedimentary diagenesis , 2002 .

[4]  K. P. Saripalli,et al.  Changes in Hydrologic Properties of Aquifer Media Due to Chemical Reactions: A Review , 2001 .

[5]  Thomas Kalbacher,et al.  Reactive transport codes for subsurface environmental simulation , 2015, Computational Geosciences.

[6]  Karsten Pruess,et al.  TOUGHREACT - A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: Applications to geothermal injectivity and CO2 geological sequestration , 2006, Comput. Geosci..

[7]  D. Pellegrini,et al.  Coupled modeling of cement/claystone interactions and radionuclide migration. , 2004, Journal of contaminant hydrology.

[8]  Satish Karra,et al.  PFLOTRAN User Manual A Massively Parallel Reactive Flow and Transport Model for Describing Surface and Subsurface Processes , 2015 .

[9]  Karsten Pruess,et al.  TOUGHREACT Version 2.0: A simulator for subsurface reactive transport under non-isothermal multiphase flow conditions , 2011, Comput. Geosci..

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

[11]  Diederik Jacques,et al.  Modeling Coupled Hydrologic and Chemical Processes: Long‐Term Uranium Transport following Phosphorus Fertilization , 2008 .

[12]  Georg Kosakowski,et al.  A class of analytical solutions for multidimensional multispecies diffusive transport coupled with precipitation‐dissolution reactions and porosity changes , 2012 .

[13]  Christophe Tournassat,et al.  Influence of reaction kinetics and mesh refinement on the numerical modelling of concrete/clay interactions , 2009 .

[14]  Nebo Jovanovic,et al.  Modelling the long-term effect of irrigation with gypsiferous water on soil and water resources , 1999 .

[15]  J. Quirk,et al.  Permeability of porous solids , 1961 .

[16]  Kerry T.B. MacQuarrie,et al.  Solution of the MoMaS reactive transport benchmark with MIN3P—model formulation and simulation results , 2010 .

[17]  Olaf Kolditz,et al.  A geochemical transport model for thermo‐hydro‐chemical (THC) coupled processes with saline water , 2011 .

[18]  Vincent Lagneau,et al.  Operator-splitting-based reactive transport models in strong feedback of porosity change: The contribution of analytical solutions for accuracy validation and estimator improvement. , 2010, Journal of contaminant hydrology.

[19]  Karsten Pruess,et al.  TOUGHREACT User's Guide: A Simulation Program for Non-isothermal Multiphase Reactive geochemical Transport in Variable Saturated Geologic Media , 2004 .

[20]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[21]  Y. Le Gallo,et al.  Coupled reaction-flow modeling of diagenetic changes in reservoir permeability, porosity and mineral compositions , 1998 .

[22]  D. Jacques,et al.  Conceptual model analysis of interaction at a concrete–Boom Clay interface , 2014 .

[23]  G. Houben.  Iron oxide incrustations in wells. Part 1: genesis, mineralogy and geochemistry , 2003 .

[24]  Christoph Clauser,et al.  Numerical simulation of pore space clogging in geothermal reservoirs by precipitation of anhydrite , 2005 .

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

[26]  Diederik Jacques,et al.  Modelling coupled water flow, solute transport and geochemical reactions affecting heavy metal migration in a podzol soil , 2008 .

[27]  S. Churakov,et al.  Exact analytical solutions for a diffusion problem coupled with a precipitation‐dissolution reaction and feedback of porosity change , 2011 .

[28]  W. D. Carrier.  Goodbye, Hazen; Hello, Kozeny-Carman , 2003 .

[29]  Stephanie Rinck-Pfeiffer,et al.  Interrelationships between biological, chemical, and physical processes as an analog to clogging in aquifer storage and recovery (ASR) wells. , 2000 .