THE KINETICS OF FLUID-ROCK INTERACTION IN CO2 STORAGE SITES: CONSTRAINTS FROM MEASUREMENTS ON NATURAL CO2 CHARGED WATERS
The progress of individual reactions, inferred from changes in groundwater chemistry is modelled using mass balance techniques. Mineral modes, in conjunction with published surface area measurements and flow rates estimated from hydraulic head measurements, are then used to quantify the kinetics of mineral dissolution.
Maximum estimated dissolution rates for plagioclase and K-feldspar are 2x10-14 and 4x10-16 mol·m-2·s-1, respectively. Fluid ion-activity products are close to equilibrium (e.g. ΔGr for plagioclase between 2 and 10 kJ/mol) and lie in the region in which mineral surface reaction rates show a strong dependence on ΔGr. Local variation in ΔGr is attributed to the injection and disassociation of CO2 which initially depresses mineral saturation in the fluid, promoting feldspar dissolution. With progressive flow through the aquifer feldspar hydrolysis reactions consumed H+ and liberate solutes to solution which increase mineral saturation in the fluid and rates slow as a consequence.
The measured plagioclase dissolution rates at low ΔGr of 2x10-14 mol·m-2·s-1 would be compatible with far-from-equilibrium rates of ~1x10-13 mol·m-2·s-1 as observed in some experimental studies. This suggests that the discrepancy between field and laboratory reactions rates may in part be explained by the differences in the thermodynamic state of natural and experimental fluids, with field-scale reactions occurring close to equilibrium whereas most laboratory experiments are run far-from-equilibrium.