Predictive modeling of CO2 sequestration in deep saline sandstone reservoirs: Impacts of geochemical kinetics

Victor N. Balashov, George D. Guthrie, J. Alexandra Hakala, Christina L. Lopano, J. Donald Rimstidt, Susan L. Brantley

Research output: Contribution to journalArticle

40 Citations (Scopus)

Abstract

One idea for mitigating the increase in fossil-fuel generated CO2 in the atmosphere is to inject CO2 into subsurface saline sandstone reservoirs. To decide whether to try such sequestration at a globally significant scale will require the ability to predict the fate of injected CO2. Thus, models are needed to predict the rates and extents of subsurface rock-water-gas interactions. Several reactive transport models for CO2 sequestration created in the last decade predicted sequestration in sandstone reservoirs of ~17 to ~90kg CO2 m-3. To build confidence in such models, a baseline problem including rock+water chemistry is proposed as the basis for future modeling so that both the models and the parameterizations can be compared systematically. In addition, a reactive diffusion model is used to investigate the fate of injected supercritical CO2 fluid in the proposed baseline reservoir+brine system. In the baseline problem, injected CO2 is redistributed from the supercritical (SC) free phase by dissolution into pore brine and by formation of carbonates in the sandstone. The numerical transport model incorporates a full kinetic description of mineral-water reactions under the assumption that transport is by diffusion only. Sensitivity tests were also run to understand which mineral kinetics reactions are important for CO2 trapping.The diffusion transport model shows that for the first ~20years (20a) after CO2 diffusion initiates, CO2 is mostly consumed by dissolution into the brine to form CO2,aq (solubility trapping). From 20 to 200a, both solubility and mineral trapping are important as calcite precipitation is driven by dissolution of oligoclase. From 200 to 1000a, mineral trapping is the most important sequestration mechanism, as smectite dissolves and calcite precipitates. Beyond 2000a most trapping is due to formation of aqueous HCO3-. Ninety-seven percent of the maximum CO2 sequestration, 34.5kg CO2 per m3 of sandstone, is attained by 4000a even though the system does not achieve chemical equilibrium until ~25,000a. This maximum represents about 20% CO2 dissolved as CO2,aq, 50% dissolved as HCO3,aq-, and 30% precipitated as calcite. The extent of sequestration as HCO3- at equilibrium can be calculated from equilibrium thermodynamics and is roughly equivalent to the amount of Na+ in the initial sandstone in a soluble mineral (here, oligoclase). Similarly, the extent of trapping in calcite is determined by the amount of Ca2+ in the initial oligoclase and smectite. Sensitivity analyses show that the rate of CO2 sequestration is sensitive to the mineral-water reaction kinetic constants between approximately 10 and 4000a. The sensitivity of CO2 sequestration to the rate constants decreases in magnitude respectively from oligoclase to albite to smectite.

Original languageEnglish (US)
Pages (from-to)41-56
Number of pages16
JournalApplied Geochemistry
Volume30
DOIs
StatePublished - Mar 1 2013

Fingerprint

Sandstone
carbon sequestration
trapping
sandstone
Minerals
Calcium Carbonate
Calcite
kinetics
Kinetics
modeling
calcite
smectite
brine
Mineral Waters
Dissolution
dissolution
reaction kinetics
mineral
Reaction kinetics
Water

All Science Journal Classification (ASJC) codes

  • Environmental Chemistry
  • Pollution
  • Geochemistry and Petrology

Cite this

Balashov, Victor N. ; Guthrie, George D. ; Hakala, J. Alexandra ; Lopano, Christina L. ; Rimstidt, J. Donald ; Brantley, Susan L. / Predictive modeling of CO2 sequestration in deep saline sandstone reservoirs : Impacts of geochemical kinetics. In: Applied Geochemistry. 2013 ; Vol. 30. pp. 41-56.
@article{4306fae31d45435ea009a83a8dc69856,
title = "Predictive modeling of CO2 sequestration in deep saline sandstone reservoirs: Impacts of geochemical kinetics",
abstract = "One idea for mitigating the increase in fossil-fuel generated CO2 in the atmosphere is to inject CO2 into subsurface saline sandstone reservoirs. To decide whether to try such sequestration at a globally significant scale will require the ability to predict the fate of injected CO2. Thus, models are needed to predict the rates and extents of subsurface rock-water-gas interactions. Several reactive transport models for CO2 sequestration created in the last decade predicted sequestration in sandstone reservoirs of ~17 to ~90kg CO2 m-3. To build confidence in such models, a baseline problem including rock+water chemistry is proposed as the basis for future modeling so that both the models and the parameterizations can be compared systematically. In addition, a reactive diffusion model is used to investigate the fate of injected supercritical CO2 fluid in the proposed baseline reservoir+brine system. In the baseline problem, injected CO2 is redistributed from the supercritical (SC) free phase by dissolution into pore brine and by formation of carbonates in the sandstone. The numerical transport model incorporates a full kinetic description of mineral-water reactions under the assumption that transport is by diffusion only. Sensitivity tests were also run to understand which mineral kinetics reactions are important for CO2 trapping.The diffusion transport model shows that for the first ~20years (20a) after CO2 diffusion initiates, CO2 is mostly consumed by dissolution into the brine to form CO2,aq (solubility trapping). From 20 to 200a, both solubility and mineral trapping are important as calcite precipitation is driven by dissolution of oligoclase. From 200 to 1000a, mineral trapping is the most important sequestration mechanism, as smectite dissolves and calcite precipitates. Beyond 2000a most trapping is due to formation of aqueous HCO3-. Ninety-seven percent of the maximum CO2 sequestration, 34.5kg CO2 per m3 of sandstone, is attained by 4000a even though the system does not achieve chemical equilibrium until ~25,000a. This maximum represents about 20{\%} CO2 dissolved as CO2,aq, 50{\%} dissolved as HCO3,aq-, and 30{\%} precipitated as calcite. The extent of sequestration as HCO3- at equilibrium can be calculated from equilibrium thermodynamics and is roughly equivalent to the amount of Na+ in the initial sandstone in a soluble mineral (here, oligoclase). Similarly, the extent of trapping in calcite is determined by the amount of Ca2+ in the initial oligoclase and smectite. Sensitivity analyses show that the rate of CO2 sequestration is sensitive to the mineral-water reaction kinetic constants between approximately 10 and 4000a. The sensitivity of CO2 sequestration to the rate constants decreases in magnitude respectively from oligoclase to albite to smectite.",
author = "Balashov, {Victor N.} and Guthrie, {George D.} and Hakala, {J. Alexandra} and Lopano, {Christina L.} and Rimstidt, {J. Donald} and Brantley, {Susan L.}",
year = "2013",
month = "3",
day = "1",
doi = "10.1016/j.apgeochem.2012.08.016",
language = "English (US)",
volume = "30",
pages = "41--56",
journal = "Applied Geochemistry",
issn = "0883-2927",
publisher = "Elsevier Limited",

}

Predictive modeling of CO2 sequestration in deep saline sandstone reservoirs : Impacts of geochemical kinetics. / Balashov, Victor N.; Guthrie, George D.; Hakala, J. Alexandra; Lopano, Christina L.; Rimstidt, J. Donald; Brantley, Susan L.

In: Applied Geochemistry, Vol. 30, 01.03.2013, p. 41-56.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Predictive modeling of CO2 sequestration in deep saline sandstone reservoirs

T2 - Impacts of geochemical kinetics

AU - Balashov, Victor N.

AU - Guthrie, George D.

AU - Hakala, J. Alexandra

AU - Lopano, Christina L.

AU - Rimstidt, J. Donald

AU - Brantley, Susan L.

PY - 2013/3/1

Y1 - 2013/3/1

N2 - One idea for mitigating the increase in fossil-fuel generated CO2 in the atmosphere is to inject CO2 into subsurface saline sandstone reservoirs. To decide whether to try such sequestration at a globally significant scale will require the ability to predict the fate of injected CO2. Thus, models are needed to predict the rates and extents of subsurface rock-water-gas interactions. Several reactive transport models for CO2 sequestration created in the last decade predicted sequestration in sandstone reservoirs of ~17 to ~90kg CO2 m-3. To build confidence in such models, a baseline problem including rock+water chemistry is proposed as the basis for future modeling so that both the models and the parameterizations can be compared systematically. In addition, a reactive diffusion model is used to investigate the fate of injected supercritical CO2 fluid in the proposed baseline reservoir+brine system. In the baseline problem, injected CO2 is redistributed from the supercritical (SC) free phase by dissolution into pore brine and by formation of carbonates in the sandstone. The numerical transport model incorporates a full kinetic description of mineral-water reactions under the assumption that transport is by diffusion only. Sensitivity tests were also run to understand which mineral kinetics reactions are important for CO2 trapping.The diffusion transport model shows that for the first ~20years (20a) after CO2 diffusion initiates, CO2 is mostly consumed by dissolution into the brine to form CO2,aq (solubility trapping). From 20 to 200a, both solubility and mineral trapping are important as calcite precipitation is driven by dissolution of oligoclase. From 200 to 1000a, mineral trapping is the most important sequestration mechanism, as smectite dissolves and calcite precipitates. Beyond 2000a most trapping is due to formation of aqueous HCO3-. Ninety-seven percent of the maximum CO2 sequestration, 34.5kg CO2 per m3 of sandstone, is attained by 4000a even though the system does not achieve chemical equilibrium until ~25,000a. This maximum represents about 20% CO2 dissolved as CO2,aq, 50% dissolved as HCO3,aq-, and 30% precipitated as calcite. The extent of sequestration as HCO3- at equilibrium can be calculated from equilibrium thermodynamics and is roughly equivalent to the amount of Na+ in the initial sandstone in a soluble mineral (here, oligoclase). Similarly, the extent of trapping in calcite is determined by the amount of Ca2+ in the initial oligoclase and smectite. Sensitivity analyses show that the rate of CO2 sequestration is sensitive to the mineral-water reaction kinetic constants between approximately 10 and 4000a. The sensitivity of CO2 sequestration to the rate constants decreases in magnitude respectively from oligoclase to albite to smectite.

AB - One idea for mitigating the increase in fossil-fuel generated CO2 in the atmosphere is to inject CO2 into subsurface saline sandstone reservoirs. To decide whether to try such sequestration at a globally significant scale will require the ability to predict the fate of injected CO2. Thus, models are needed to predict the rates and extents of subsurface rock-water-gas interactions. Several reactive transport models for CO2 sequestration created in the last decade predicted sequestration in sandstone reservoirs of ~17 to ~90kg CO2 m-3. To build confidence in such models, a baseline problem including rock+water chemistry is proposed as the basis for future modeling so that both the models and the parameterizations can be compared systematically. In addition, a reactive diffusion model is used to investigate the fate of injected supercritical CO2 fluid in the proposed baseline reservoir+brine system. In the baseline problem, injected CO2 is redistributed from the supercritical (SC) free phase by dissolution into pore brine and by formation of carbonates in the sandstone. The numerical transport model incorporates a full kinetic description of mineral-water reactions under the assumption that transport is by diffusion only. Sensitivity tests were also run to understand which mineral kinetics reactions are important for CO2 trapping.The diffusion transport model shows that for the first ~20years (20a) after CO2 diffusion initiates, CO2 is mostly consumed by dissolution into the brine to form CO2,aq (solubility trapping). From 20 to 200a, both solubility and mineral trapping are important as calcite precipitation is driven by dissolution of oligoclase. From 200 to 1000a, mineral trapping is the most important sequestration mechanism, as smectite dissolves and calcite precipitates. Beyond 2000a most trapping is due to formation of aqueous HCO3-. Ninety-seven percent of the maximum CO2 sequestration, 34.5kg CO2 per m3 of sandstone, is attained by 4000a even though the system does not achieve chemical equilibrium until ~25,000a. This maximum represents about 20% CO2 dissolved as CO2,aq, 50% dissolved as HCO3,aq-, and 30% precipitated as calcite. The extent of sequestration as HCO3- at equilibrium can be calculated from equilibrium thermodynamics and is roughly equivalent to the amount of Na+ in the initial sandstone in a soluble mineral (here, oligoclase). Similarly, the extent of trapping in calcite is determined by the amount of Ca2+ in the initial oligoclase and smectite. Sensitivity analyses show that the rate of CO2 sequestration is sensitive to the mineral-water reaction kinetic constants between approximately 10 and 4000a. The sensitivity of CO2 sequestration to the rate constants decreases in magnitude respectively from oligoclase to albite to smectite.

UR - http://www.scopus.com/inward/record.url?scp=84875395417&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84875395417&partnerID=8YFLogxK

U2 - 10.1016/j.apgeochem.2012.08.016

DO - 10.1016/j.apgeochem.2012.08.016

M3 - Article

AN - SCOPUS:84875395417

VL - 30

SP - 41

EP - 56

JO - Applied Geochemistry

JF - Applied Geochemistry

SN - 0883-2927

ER -