Coupled Vadose Zone and Atmospheric Surface-Layer Transport of Carbon Dioxide from Geologic Carbon Sequestration Sites

Geologic CO 2 sequestration is being considered as a way to offset fossil fuel–related CO 2 emissions to reduce the rate of increase of atmospheric CO 2 concentrations. The accumulation of vast quantities of injected CO 2 in geologic sequestration sites may entail health and environmental risks from potential leakage and seepage of CO 2 into the near-surface environment. We are developing and applying a coupled subsurface and atmospheric surface-layer modeling capability built within the framework of the integral finite difference reservoir simulator TOUGH2. The overall purpose of the modeling studies is to predict CO 2 concentration distributions under a variety of seepage scenarios and geologic, hydrologic, and atmospheric conditions. These concentration distributions will provide the basis for determining aboveground and near-surface instrumentation needs for CO 2 sequestration monitoring and verification, as well as for assessing health, safety, and environmental risks. A key feature of CO 2 is its large density (ρ = 1.8 kg m −3 ) relative to air (ρ = 1.2 kg m −3 ), a property that may allow small leaks to cause concentrations in air above the occupational exposure limit of 4% in low-lying and enclosed areas such as valleys and basements where dilution rates are low. The approach we take to coupled modeling involves development of T2CA, a TOUGH2 module for modeling the multicomponent transport of water, brine, CO 2 , gas tracer, and air in the subsurface. For the atmospheric surface-layer advection and dispersion, we use a logarithmic vertical velocity profile to specify constant time-averaged ambient winds, and atmospheric dispersion approaches to model mixing due to eddies and turbulence. Initial simulations with the coupled model suggest that atmospheric dispersion quickly dilutes diffuse CO 2 seepage fluxes to negligible concentrations, and that rainfall infiltration can cause CO 2 to return to the subsurface as a dissolved component in infiltrating rainwater.

[1]  F. Gifford,et al.  Use of routine meteorological observations for estimating atmospheric dispersion , 1961 .

[2]  R. Stull An Introduction to Boundary Layer Meteorology , 1988 .

[3]  Joan M. Ogden,et al.  Carbon sequestration research and development , 1999 .

[4]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[5]  K. Pruess,et al.  TOUGH2 User's Guide Version 2 , 1999 .

[6]  Rex Britter,et al.  ATMOSPHERIC DISPERSION OF DENSE GASES , 1989 .

[7]  S. Arya Air Pollution Meteorology and Dispersion , 1998 .

[8]  Donald Golder,et al.  Relations among stability parameters in the surface layer , 1972 .

[9]  K. W. Steinberg,et al.  Overview of Petroleum Environmental Research Forum (PERF) dense gas dispersion modeling project , 2001 .

[10]  Stefan Bachu,et al.  Aquifer disposal of CO2: Hydrodynamic and mineral trapping , 1994 .

[11]  Curtis M. Oldenburg,et al.  On Leakage and Seepage from Geologic Carbon Sequestration Sites: Unsaturated Zone Attenuation , 2003 .

[12]  Curtis M. Oldenburg,et al.  Near-surface monitoring strategies for geologic carbon dioxide storage verification , 2003 .

[13]  C. Oldenburg,et al.  Simulation of propagating fronts in geothermal reservoirs with the implicit Leonard total variation diminishing scheme , 2000 .

[14]  Curtis M. Oldenburg,et al.  EOS7R: Radionuclide transport for TOUGH2 , 1995 .

[15]  Dennis D. Baldocchi,et al.  Modeling CO2 and water vapor exchange of a temperate broadleaved forest across hourly to decadal time scales , 2001 .