Reactive transport modelling of CO2 storage in saline aquifers to elucidate fundamental processes, trapping mechanisms and sequestration partitioning

Abstract The ultimate fate of CO2 injected into saline aquifers for environmental isolation is governed by three interdependent yet conceptually distinct processes: CO2 migration as a buoyant immiscible fluid phase, direct chemical interaction of this rising plume with ambient saline waters, and its indirect chemical interaction with aquifer and caprock minerals through the aqueous wetting phase. Each process is directly linked to a corresponding trapping mechanism: immiscible plume migration to hydrodynamic trapping, plume-water interaction to solubility trapping, and plume-mineral interaction to mineral trapping. In this study, reactive transport modelling of CO2 storage in a shele-capped sandstone aquifer at Sleipner has elucidated and established key parametric dependencies of these fundamental processes, the associated trapping mechanisms, and sequestration partitioning among them during consecutive ten-year prograde (active-injection) and retrograde (post-injection) regimes. Intra-aquifer permeability structure controls the path of immiscible CO2 migration, thereby establishing the spatial framework of plume-aquifer interaction and the potential effectiveness of solubility and mineral trapping. Inter-bedded thin shales, which occur at Sleipner, retard vertical and promote lateral plume migration, thereby significantly expanding this framework and enhancing this potential. Actual efficacy of these trapping mechanisms is determined by compositional characteristics of the aquifer and caprock: the degree of solubility trapping decreases with increasing formation-water salinity, whereas that of mineral trapping is proportional to the bulk concentration of carbonate-forming elements, principally Fe, Mg, Ca, Na and Al. In the near-field environment of Sleipner-like settings, 80–85% by mass of injected CO2 remains and migrates as an immiscible fluid phase, 15–20% dissolves into formation waters, and less than 1% precipitates as carbonate minerals. This partitioning defines the relative effectiveness of hydrodynamic, solubility, and mineral trapping on a mass basis. Seemingly inconsequential, mineral trapping has enormous strategic significance: it maintains injectivity, delineates the storage volume, and improves caprock integrity. Four distinct mechanisms have been identified: dawsonite [NaAlCO3(OH)2] cementation occurs throughout the intra-aquifer plume, while calcite-group carbonates [principally (Fe,Mg,Ca)CO3] precipitate via disparate processes along lateral and upper plume margins, and by yet another process within inter-bedded and caprock shales. The coupled mineral dissolution/precipitation reaction associated with each mechanism reduces local porosity and permeability. For Sleipner-like settings, the magnitude of such reduction for dawsonite cementation is near negligible; hence, this process effectively maintains initial CO2 injectivity. Of similarly small magnitude is the reduction associated with formation of carbonate rind along upper and lateral plume boundaries; these processes effectively delineate the CO2 storage volume, and for saline aquifers anomalously rich in Fe-Mg-Ca may partially self-seal the plume. Porosity and permeability reduction is most extreme within shales, because their clay-rich mineralogy defines bulk Fe-Mg concentrations much greater than those of saline aquifers. In the basal caprock shale of our models, these reductions amount to 4.5 and 13%, respectively, after the prograde regime. During the retrograde phase, residual saturation of immiscible CO2 maintains the prograde extent of solubility trapping while continuously enhancing that of mineral trapping. At the close of our 20-year simulations, initial porosity and permeability of the basal caprock shale have been reduced by 8 and 22%, respectively. Extrapolating to hypothetical complete consumption of Fe-Mg-bearing shale minerals (here 10 vol.% Mg-chlorite) yields an ultimate reduction of about 52 and 90%, respectively, after 130 years. Hence, the most crucial strategic impact of mineral trapping in Sleipner-like settings: it continuously improves hydrodynamic seal integrity of the caprock and, therefore, containment of the immiscible plume and solubility-trapped CO2.

[1]  A. Blum,et al.  Chapter 7. FELDSPAR DISSOLUTION KINETICS , 1995 .

[2]  J. Nitao,et al.  Reactive Transport Modeling of Geologic CO2 Sequestration at Sleipner , 2003 .

[3]  P. Bergmo,et al.  The Long-Term Fate of CO2 Injected into an Aquifer , 2003 .

[4]  W. Wakeham,et al.  The Viscosity of Carbon Dioxide , 1998 .

[5]  J. C. Baker,et al.  Continental-scale magmatic carbon-dioxide seepage recorded by Dawsonite in the Bowen-Gunnedah-Sydney Basin System, eastern Australia , 1995 .

[6]  Lincoln Paterson,et al.  Rate of Dissolution Due to Convective Mixing in the Underground Storage of Carbon Dioxide , 2003 .

[7]  K. Knauss,et al.  The effect of malonate on the dissolution kinetics of albite, quartz, and microcline as a function of pH at 70°C , 1995 .

[8]  H. Barnes,et al.  Precipitation and dissolution rate constants for cristobalite from 150 to 300°C , 1995 .

[9]  R. Garrels,et al.  Solutions, Minerals and Equilibria , 1965 .

[10]  Jack C. Parker,et al.  A parametric model for constitutive properties governing multiphase flow in porous media , 1987 .

[11]  K. Knauss,et al.  The dissolution kinetics of quartz as a function of pH and time at 70°C , 1988 .

[12]  A. Blum,et al.  Feldspar dissolution kinetics , 1995 .

[13]  R. A. Robie,et al.  Thermodynamic properties of minerals and related substances at 298.15 K and 1 bar (10[5] pascals) pressure and at higher temperatures , 1995 .

[14]  Guntram Jordan,et al.  Dissolution Rates of Calcite (104) Obtained by Scanning Force Microscopy: Microtopography-Based Dissolution Kinetics on Surfaces with Anisotropic Step Velocities , 1998 .

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

[16]  E. Oelkers,et al.  Calculation of the thermodynamic properties of aqueous species at high pressures and temperatures. Effective electrostatic radii, dissociation constants and standard partial molal properties to 1000 °C and 5 kbar , 1992 .

[17]  H. Helgeson,et al.  Thermodynamics of hydrothermal systems at elevated temperatures and pressures , 1969 .

[18]  N. P. Fofonoff,et al.  Algorithms for Computation of Fundamental Properties of Seawater. Endorsed by Unesco/SCOR/ICES/IAPSO Joint Panel on Oceanographic Tables and Standards and SCOR Working Group 51. Unesco Technical Papers in Marine Science, No. 44. , 1983 .

[19]  Ola Eiken,et al.  Seismic Monitoring of CO2 Injected Into a Marine Acquifer , 2000 .

[20]  Oleg S. Pokrovsky,et al.  Processes at the magnesium-bearing carbonates/solution interface. II. kinetics and mechanism of magnesite dissolution. , 1999 .

[21]  Carl I. Steefel,et al.  Reactive Transport Modeling of Geologic CO{sub 2} Sequestration in Saline Aquifers: The Influence of Intra-Aquifer Shales and the Relative Effectiveness of Structural, Solubility, and Mineral Trapping During Prograde and Retrograde Sequestration , 2001 .

[22]  J. Tester,et al.  Correlating quartz dissolution kinetics in pure water from 25 to 625°C , 1994 .

[23]  Maxwell N. Watson,et al.  The Ladbroke Grove-Katnook carbon dioxide natural laboratory: A recent CO2 accumulation in a lithic sandstone reservoir , 2004 .

[24]  J. Drever,et al.  rates of feldspar dissolution at pH 3–7 with 0–8 m M oxalic acid , 1996 .

[25]  K. L. Nagy,et al.  Chapter 5. DISSOLUTION AND PRECIPITATION KINETICS OF SHEET SILICATES , 1995 .

[26]  E. C. Beutner Slaty cleavage and related strain in Martinsburg Slate, Delaware Water Gap, New Jersey , 1978 .

[27]  Denis Norton,et al.  Theory of Hydrothermal Systems , 1984 .

[28]  James W. Johnson,et al.  Critical phenomena in hydrothermal systems; state, thermodynamic, electrostatic, and transport properties of H 2 O in the critical region , 1991 .

[29]  Susan D. Hovorka,et al.  Frio brine sequestration pilot in the Texas Gulf Coast , 2003 .

[30]  Kevin G. Knauss,et al.  Reactive transport modeling of plug-flow reactor experiments: quartz and tuff dissolution at 240°C , 1998 .

[31]  W. Wagner,et al.  A New Equation of State for Carbon Dioxide Covering the Fluid Region from the Triple‐Point Temperature to 1100 K at Pressures up to 800 MPa , 1996 .

[32]  H. Helgeson,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures; revised equations of state for the standard partial molal properties of ions and electrolytes , 1988 .

[33]  Adrian E. Scheidegger,et al.  The physics of flow through porous media , 1957 .

[34]  E. Oelkers,et al.  SUPCRT92: a software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000 ° C , 1992 .

[35]  H. Helgeson,et al.  Summary and critique of the thermodynamic properties of rock forming minerals , 1978 .

[36]  J. H. Stewart Origin of Sedimentary Rocks , 1974 .

[37]  K. Knauss,et al.  Muscovite dissolution kinetics as a function of pH and time at 70°C , 1989 .

[38]  Mark H. Holtz,et al.  Residual Gas Saturation to Aquifer Influx: A Calculation Method for 3-D Computer Reservoir Model Construction , 2002 .

[39]  Lincoln Paterson,et al.  Role of Convective Mixing in the Long-Term Storage of Carbon Dioxide in Deep Saline Formations , 2005 .

[40]  B. Fritz,et al.  Hydrothermal experiments and thermo-kinetic modelling of water-sandstone interactions☆ , 1994 .

[41]  K. Knauss,et al.  Dependence of albite dissolution kinetics on ph and time at 25°c and 70°c , 1986 .

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

[43]  Kathryn L. Nagy,et al.  Dissolution and precipitation kinetics of sheet silicates , 1995 .