Reactive transport modeling to assess geochemical monitoring for detection of CO2 intrusion into shallow aquifers

Abstract The hypothesis is tested if changes in electric conductivity of groundwater (EC) in response to gaseous CO 2 intrusion are sufficient to be detected using probe measurements and geophysical electromagnetic measurements, e.g. airborne electromagnetic measurements. Virtual reactive scenario modelling is used to simulate the effects of the presence of calcite, CO 2 intrusion rates, depth of the aquifer formation, initial salinity of groundwater and CO 2 intrusion time on changes in EC. In all simulations, EC rises rapidly in response to CO 2 intrusion, however in different magnitudes. When calcite is present, EC changes are strong (+1.11 mS/cm after 24 hours of CO 2 intrusion) mainly due to calcite dissolution, whereas in aquifers without calcite changes are very low (+0.02 mS/cm after 24 hours) and close to the resolution range of probes. Increased depth (250 m / 500 m), i.e. higher temperature and pressure, and higher intrusion rates (up to full saturation) result in stronger rises in EC (+5.08 mS/cm in 500 m depth and 100 % saturation), and initial salinity has a negligible influence on changes in EC. Temporally limited CO 2 intrusion leads to EC values close to pre- CO 2 -intrusion-levels in the long-term. Measurement resolution of commercial EC probes is sufficient to detect CO 2 intrusion in almost all cases. In terms of geophysical electromagnetic measurements, applications in the field of monitoring saltwater-freshwater interfaces indicate a sufficient measurement resolution to detect changes in all simulations. However, practical limitations are expected due to the dependence of measurement resolutions on the applied measurement devices and site-specific geological settings.

[1]  Jens Birkholzer,et al.  On mobilization of lead and arsenic in groundwater in response to CO2 leakage from deep geological storage , 2009 .

[2]  Yue Hao,et al.  Geochemical detection of carbon dioxide in dilute aquifers , 2009, Geochemical transactions.

[3]  David R. Cole,et al.  Potential environmental issues of CO2 storage in deep saline aquifers: Geochemical results from the Frio-I Brine Pilot test, Texas, USA , 2009 .

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

[5]  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 .

[6]  James W. Ball,et al.  User's manual for WATEQ4F, with revised thermodynamic data base and text cases for calculating speciation of major, trace, and redox elements in natural waters , 1991 .

[7]  P. Wentzell,et al.  Models for Conductance Measurements in Quality Assurance of Water Analysis , 1994 .

[8]  S. Brantley,et al.  Chemical weathering rates of silicate minerals , 1995 .

[9]  Using helicopter electromagnetic data to predict groundwater quality in fractured crystalline bedrock in a semi-arid region, Northeast Brazil , 2010 .

[10]  Zhenhao Duan,et al.  An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533 K and from 0 to 2000 bar , 2003 .

[11]  Yousif K. Kharaka,et al.  A Compilation of Rate Parameters of Water-Mineral Interaction Kinetics for Application to Geochemical Modeling , 2004 .

[12]  Richard Hammack,et al.  Near-surface monitoring for the ZERT shallow CO2 injection project , 2009 .

[13]  A. Lasaga Kinetic theory in the earth sciences , 1998 .