Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Reactive Transport Modelling of CO2 Storage in Saline Aquifers to Elucidate Fundamental Processes, Trapping Mechanisms, and Sequestration Partitioning

Journal Article · · Geological Storage of Carbon Dioxide
The ultimate fate of CO{sub 2} injected into saline aquifers for environmental isolation is governed by three interdependent yet conceptually distinct processes: CO{sub 2} 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 cap-rock 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 modeling of CO{sub 2} storage in a shale-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 10-year prograde (active-injection) and retrograde (post-injection) regimes. Intra-aquifer permeability structure controls the path of immiscible CO{sub 2} 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 cap rock: the degree of solubility trapping decreases with increasing formation-water salinity, while 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 CO{sub 2} 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 cap-rock integrity. We have identified four distinct mechanisms: dawsonite [NaAlCO{sub 3}(OH){sub 2}] cementation occurs throughout the intra-aquifer plume, while calcite-group carbonates [principally, (Fe,Mg,Ca)CO{sub 3}] precipitate via disparate processes along lateral and upper plume margins, and by yet another process within inter-bedded and cap-rock 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 CO{sub 2} 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 CO{sub 2} 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 cap-rock 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 CO{sub 2} 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 cap-rock 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 cap rock and, therefore, containment of the immiscible plume and solubility trapped CO{sub 2}.
Research Organization:
Lawrence Livermore National Lab., Livermore, CA (US)
Sponsoring Organization:
US Department of Energy (US)
DOE Contract Number:
W-7405-ENG-48
OSTI ID:
15015125
Report Number(s):
UCRL-JRNL-205627
Journal Information:
Geological Storage of Carbon Dioxide, Journal Name: Geological Storage of Carbon Dioxide
Country of Publication:
United States
Language:
English

References (27)

Reactive Transport Modeling of Geohazards Associated with CO2 Injection for EOR and Geologic Sequestration conference April 2013
Thermodynamics of hydrothermal systems at elevated temperatures and pressures journal September 1969
The Viscosity of Carbon Dioxide journal January 1998
Critical phenomena in hydrothermal systems; state, thermodynamic, electrostatic, and transport properties of H 2 O in the critical region journal June 1991
Reactive transport modeling of plug-flow reactor experiments: quartz and tuff dissolution at 240°C journal August 1998
Reactive Transport Modeling of Geologic CO2 Sequestration at Sleipner conference September 2007
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 journal January 1992
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 journal January 1988
Hydrothermal experiments and thermo-kinetic modelling of water-sandstone interactions journal October 1994
Theory of Hydrothermal Systems journal May 1984
rates of feldspar dissolution at pH 3–7 with 0–8 m M oxalic acid journal October 1996
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 journal November 1996
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 journal August 1992
Processes at the magnesium-bearing carbonates/solution interface. II. kinetics and mechanism of magnesite dissolution. journal July 1999
Dissolution Rates of Calcite (104) Obtained by Scanning Force Microscopy: Microtopography-Based Dissolution Kinetics on Surfaces with Anisotropic Step Velocities journal March 1998
The effect of malonate on the dissolution kinetics of albite, quartz, and microcline as a function of pH at 70°C journal January 1995
Dependence of albite dissolution kinetics on ph and time at 25°c and 70°c journal November 1986
A parametric model for constitutive properties governing multiphase flow in porous media journal April 1987
The dissolution kinetics of quartz as a function of pH and time at 70°C journal January 1988
Seismic monitoring of CO 2 injected into a marine acquifer conference March 2012
Correlating quartz dissolution kinetics in pure water from 25 to 625°C journal June 1994
Chapter 5. Dissolution and Precipitation Kinetics of Sheet Silicates book December 1995
The Ladbroke Grove-Katnook Carbon Dioxide Natural LaboratoryA Recent CO2 Accumulation in a Lithic Sandstone Reservoir book January 2003
Precipitation and dissolution rate constants for cristobalite from 150 to 300°C journal January 1995
Frio Brine Sequestration Pilot in the Texas Gulf Coast book January 2003
Muscovite dissolution kinetics as a function of pH and time at 70°C journal July 1989
The Long-Term Fate of CO2 Injected into an Aquifer book January 2003