CO<Subscript>2</Subscript> leakage through an abandoned well: problem-oriented benchmarks

The efficiency and sustainability of carbon dioxide (CO₂) storage in deep geological formations crucially depends on the integrity of the overlying cap-rocks. Existing oil and gas wells, which penetrate the formations, are potential leakage pathways. This problem has been discussed in the literature, and a number of investigations using semi-analytical mathematical approaches have been carried out by other authors to quantify leakage rates. The semi-analytical results are based on a number of simplifying assumptions. Thus, it is of great interest to assess the influence of these assumptions. We use a numerical model to compare the results with those of the semi-analytical model. Then we ease the simplifying restrictions and include more complex thermodynamic processes including sub- and supercritical fluid properties of CO₂ and non-isothermal as well as compositional effects. The aim is to set up problem-oriented benchmark examples that allow a comparison of different modeling approaches to the problem of CO₂ leakage.     

Dynamics of CO<Subscript>2</Subscript> fluxes and concentrations during a shallow subsurface CO<Subscript>2</Subscript> release

A field facility located in Bozeman, Montana provides the opportunity to test methods to detect, locate, and quantify potential CO₂ leakage from geologic storage sites. From 9 July to 7 August 2008, 0.3 t CO₂ day⁻¹ were injected from a 100-m long, ~2.5-m deep horizontal well. Repeated measurements of soil CO₂ fluxes on a grid characterized the spatio-temporal evolution of the surface leakage signal and quantified the surface leakage rate. Infrared CO₂ concentration sensors installed in the soil at 30 cm depth at 0–10 m from the well and at 4 cm above the ground at 0 and 5 m from the well recorded surface breakthrough of CO₂ leakage and migration of CO₂ leakage through the soil. Temporal variations in CO₂ concentrations were correlated with atmospheric and soil temperature, wind speed, atmospheric pressure, rainfall, and CO₂ injection rate.     

Simulation of CO<Subscript>2</Subscript> hydrate formation in cold aquifers: nonequilibrium approach

In this paper, we focus on the geological storage of CO₂ in reservoirs with zones that are cold enough to facilitate CO₂ hydrate formation at local pressures. A 2D hydro-chemical mechanical model which has five layers (three layers with aquifers and two layers with cap rock in which we introduced two fractures) is created. We apply a reactive transport reservoir simulator, RetrasoCodeBright (RCB), in which hydrate is treated as a pseudo mineral. Following the recent modifications to account for hydrate dynamics in the code through a kinetic approach (Kvamme et al., Proceedings of the 7th International Conference on Gas Hydrates (ICGH 2011), 2011b), we have further improved the simulator to implement the nonequilibrium thermodynamic calculations. In the present study, we spot the light on the hydrate formation effects on porosity in different regions, as well as on the flow pattern. These simulations are based on classical relationships between porosity and permeability, but the outline of ongoing modifications is presented as well. A critical question in such systems is whether hydrate formation can contribute to stabilizing the storage, given that hydrates are pore filling and cannot be stable toward mineral surfaces. The implications of hydrate formation on the geo-mechanical properties of the model reservoir are other aspects addressed in this study.     

A stochastic inversion workflow for monitoring the distribution of CO<Subscript>2</Subscript> injected into deep saline aquifers

We present a two-step stochastic inversion approach for monitoring the distribution of CO₂ injected into deep saline aquifers for the typical scenario of one single injection well and a database comprising a common suite of well logs as well as time-lapse vertical seismic profiling (VSP) data. In the first step, we compute several sets of stochastic models of the elastic properties using conventional sequential Gaussian co-simulations (SGCS) representing the considered reservoir before CO₂ injection. All realizations within a set of models are then iteratively combined using a modified gradual deformation algorithm aiming at reducing the mismatch between the observed and simulated VSP data. In the second step, these optimal static models then serve as input for a history matching approach using the same modified gradual deformation algorithm for minimizing the mismatch between the observed and simulated VSP data following the injection of CO₂. At each gradual deformation step, the injection and migration of CO₂ is simulated and the corresponding seismic traces are computed and compared with the observed ones. The proposed stochastic inversion approach has been tested for a realistic, and arguably particularly challenging, synthetic case study mimicking the geological environment of a potential CO₂ injection site in the Cambrian-Ordivician sedimentary sequence of the St. Lawrence platform in Southern Québec. The results demonstrate that the proposed two-step reservoir characterization approach is capable of adequately resolving and monitoring the distribution of the injected CO₂. This finds its expression in optimized models of P- and S-wave velocities, density, and porosity, which, compared to conventional stochastic reservoir models, exhibit a significantly improved structural similarity with regard to the corresponding reference models. The proposed approach is therefore expected to allow for an optimal injection forecast by using a quantitative assimilation of all available data from the appraisal stage of a CO₂ injection site.     

Feasibility of CO<Subscript>2</Subscript> migration detection using pressure and CO<Subscript>2</Subscript> saturation monitoring above an imperfect primary seal of a geologic CO<Subscript>2</Subscript> storage formation: a numerical investigation

A numerical model was developed to investigate the potential to detect fluid migration in a (homogeneous, isotropic, with constant pressure lateral boundaries) porous and permeable interval overlying an imperfect primary seal of a geologic CO₂ storage formation. The seal imperfection was modeled as a single higher-permeability zone in an otherwise low-permeability seal, with the center of that zone offset from the CO₂ injection well by 1400 m. Pressure response resulting from fluid migration through the high-permeability zone was detectable up to 1650 m from the centroid of that zone at the base of the monitored interval after 30 years of CO₂ injection (detection limit = 0.1 MPa pressure increase); no pressure response was detectable at the top of the monitored interval at the same point in time. CO₂ saturation response could be up to 774 m from the center of the high-permeability zone at the bottom of the monitored interval, and 1103 m at the top (saturation detection limit = 0.01). More than 6% of the injected CO₂, by mass, migrated out of primary containment after 130 years of site performance (including 30 years of active injection) in the case where the zone of seal imperfection had a moderately high permeability (10^??17 m² or 0.01 mD). Free-phase CO₂ saturation monitoring at the top of the overlying interval provides favorable spatial coverage for detecting fluid migration across the primary seal. Improved sensitivity of detection for pressure perturbation will benefit time of detection above an imperfect seal.     

A generalized cubic equation of state with application to pure CO<Subscript>2</Subscript> injection in aquifers

A generalized cubic equation of state is given. The Peng-Robinson and the Soave-Redlich-Kwong equations are special cases of this equation. The generalized equation of state is precisely as simple and computationally efficient as these classical equations. Through comparison with the Span-Wagner equation for CO ₂, we obtain an improved density accuracy in predefined temperature-pressure domains. The generalized equation is then verified through two relevant examples of CO ₂ injection and migration. Comparisons are made with other standard cubic EOS in order to show the range of solutions obtained with less accurate EOS.     

LiFeSi<Subscript>2</Subscript>O<Subscript>6</Subscript> and NaFeSi<Subscript>2</Subscript>O<Subscript>6</Subscript> at low temperatures: an infrared spectroscopic study

 Synthetic aegirine LiFeSi₂O₆ and NaFeSi₂O₆ were characterized using infrared spectroscopy in the frequency range 50–2000 cm⁻¹, and at temperatures between 20 and 300 K. For the C2/c phase of LiFeSi₂O₆, 25 of the 27 predicted infrared bands and 26 of 30 predicted Raman bands are recorded at room temperature. NaFeSi₂O₆ (with symmetry C2/c) shows 25 infrared and 26 Raman bands. On cooling, the C2/c–P2₁/c structural phase transition of LiFeSi₂O₆ is characterized by the appearance of 13 additional recorded peaks. This observation indicates the enlargement of the unit cell at the transition point. The appearance of an extra band near 688 cm⁻¹ in the monoclinic P2₁/c phase, which is due to the Si–O–Si vibration in the Si₂O₆ chains, indicates that there are two non-equivalent Si sites with different Si–O bond lengths. Most significant spectral changes appear in the far-infrared region, where Li–O and Fe–O vibrations are mainly located. Infrared bands between 300 and 330 cm⁻¹ show unusually dramatic changes at temperatures far below the transition. Compared with the infrared data of NaFeSi₂O₆ measured at low temperatures, the change in LiFeSi₂O₆ is interpreted as the consequence of mode crossing in the frequency region. A generalized Landau theory was used to analyze the order parameter of the C2/c–P2₁/c phase transition, and the results suggest that the transition is close to tricritical. Received: 21 January 2002 / Accepted: 22 July 2002     

Accessing effects and signals of leakage from a CO<Subscript>2</Subscript> reservoir to a shallow freshwater aquifer by reactive transport modelling

CO₂ geological storage is a transitional technology for the mitigation of climate change. In the vicinity of potential CO₂ reservoirs in Hungary, protected freshwater aquifers used for drinking water supplies exist. Effects of disaster events of CO₂ escape and brine displacement to one of these aquifers have been studied by kinetic 1D reactive transport modelling in PHREEQC. Besides verifying that ion concentrations in the freshwater may increase up to drinking water limit values in both scenarios (CO₂ or brine leakage), total porosity of the rock is estimated. Pore volume is expected to increase at the entry point of CO₂ and to decrease at further distances, whereas it shows minor increase along the flow path for the effect of brine inflow. Additionally, electrical conductivity of water is estimated and suggested to be the best parameter to measure for cost-effective monitoring of both worst-case leakage scenarios.     

Preparation of carbon molecular sieves and its impregnation with Co and Ni for CO<Subscript>2</Subscript>/N<Subscript>2</Subscript> separation

The carbon molecular sieves (CMSs) prepared by carbonaceous materials as precursors are effective in CO₂/N₂ separation. However, selectivity of these materials is too low, since hydrocarbon cracking for developing the desired microporosity in carbonaceous materials has not been done effectively. Hence, in this study, cobalt and nickel impregnation on the precursor was conducted to introduce catalysts for hydrocarbon cracking. Cobalt and nickel impregnation, carbonization under N₂ atmosphere, and chemical vapor deposition (CVD) by benzene were conducted on the extruded mixtures of activated carbon and coal tar pitch under different conditions to prepare CMSs. The best CMS prepared by carbon deposition on the cobalt-impregnated samples exhibited CO₂ adsorption capacity of 54.79 mg/g and uptake ratio of 28.9 at 0 °C and 1 bar. In terms of CO₂ adsorption capacity and uptake ratio, CMSs prepared by carbon deposition on non-impregnated and cobalt-impregnated samples presented the best results, respectively. As benzene concentration and CVD time increased, equilibrium adsorption capacity of CO₂ decreased, and uptake ratio increased. Cobalt was found to be the best catalyst for benzene cracking in the CVD process.     

CO<Subscript>2</Subscript> mineralization of natural wollastonite into porous silica and CaCO<Subscript>3</Subscript> powders promoted via membrane electrolysis

CO₂ is a greenhouse gas, whose emissions threaten the existence of human beings. Its inherently safe sequestration can be performed via CO₂ mineralization, which is relatively slow under natural conditions. In this work, an energy-saving membrane electrolysis technique was proposed for accelerating the CO₂ mineralization of wollastonite into SiO₂ and CaCO₃ products. The electrolysis process involved splitting NH₄Cl into HCl and NH₃·H₂O via hydrogen oxidation and water reduction at the anode and cathode of the electrolytic system, respectively. In contrast to the chlor-alkali electrolysis, this method did not involve Cl^? oxidation and the standard potential of the anode was reduced. Additionally, NH₄Cl was used as the electrolyte instead of NaCl; as a result, the generation of NH₃·H₂O instead of NaOH occurred in the catholyte and the cathodic pH dramatically decreased, thus reducing the cathodic potential for hydrogen evolution. The observed changes led to a 73.5% decrease in the energy consumption. Moreover, after the process of CO₂ mineralization was optimized, SiO₂ with a specific surface area of 221.8 m² g^?1 and CaCO₃ with a purity of 99.9% were obtained.     

Experimental study of pure and mixtures of CO<Subscript>2</Subscript> and CH<Subscript>4</Subscript> adsorption on modified carbon nanotubes

In this work 3-[2-(2-aminoethylamino)ethylamino]propyl trimethoxysilane (TRI) was employed to functionalize MWCNT containing hydroxyl groups (OH-MWCNT), and the XRD, FTIR, TGA and CHNS elemental analysis techniques were used to characterize the resulted adsorbents. The characterization results for amine-MWCNT showed amine groups effectively attached to the surface of the MWCNT. The equilibrium adsorption capacity of pure CO₂ and CH₄ and their binary mixture on the pristine MWCNT, OH-MWCNT and amine-MWCNT was measured through a set of equilibrium adsorption experiments at 303.2 and 318.2 K. Capacities of all three types of adsorbents for CO₂ adsorption were higher than those for methane adsorption. Also, amine-MWCNT demonstrated better performance on CO₂ adsorption than the other two adsorbents, especially at low partial pressures. The capacity of amine-MWCNT for pure CO₂ adsorption was 2.5 and 4 times as much as those for pristine MWCNT and OH-MWCNT, respectively, at the temperature of 303.2 K and the pressure of 0.2 bar. The binary adsorption experiment revealed that CO₂/CH₄ selectivity for pristine MWCNT and amine-MWCNT in all molar fractions of CO₂ is about 1.77 and 7, respectively.     

Impacts of mineralogical compositions on different trapping mechanisms during long-term CO<Subscript>2</Subscript> storage in deep saline aquifers

Deep saline aquifers in sedimentary basins are considered to have the greatest potential for CO₂ geological storage in order to reduce carbon emissions. CO₂ injected into a saline sandstone aquifer tends to migrate upwards toward the caprock because the density of the supercritical CO₂ phase is lower than that of formation water. The accumulated CO₂ in the upper portions of the reservoir gradually dissolves into brine, lowers pH and changes the aqueous complexation, whereby induces mineral alteration. In turn, the mineralogical composition could impose significant effects on the evolution of solution, further on the mineralized CO₂. The high density of aqueous phase will then move downward due to gravity, give rise to “convective mixing,” which facilitate the transformation of CO₂ from the supercritical phase to the aqueous phase and then to the solid phase. In order to determine the impacts of mineralogical compositions on trapping amounts in different mechanisms for CO₂ geological storage, a 2D radial model was developed. The mineralogical composition for the base case was taken from a deep saline formation of the Ordos Basin, China. Three additional models with varying mineralogical compositions were carried out. Results indicate that the mineralogical composition had very obvious effects on different CO₂ trapping mechanisms. Specific to our cases, the dissolution of chlorite provided Mg²⁺ and Fe²⁺ for the formation of secondary carbonate minerals (ankerite, siderite and magnesite). When chlorite was absent in the saline aquifer, the dominant secondary carbon sequestration mineral was dawsonite, and the amount of CO₂ mineral trapping increased with an increase in the concentration of chlorite. After 3000 years, 69.08, 76.93, 83.52 and 87.24 % of the injected CO₂ can be trapped in the solid (mineral) phase, 16.05, 11.86, 8.82 and 6.99 % in the aqueous phase, and 14.87, 11.21, 7.66 and 5.77 % in the gas phase for Case 1 through 4, respectively.     

CO<Subscript>2</Subscript> fluid inclusions in Jack Hills zircons

The discovery of Hadean to Paleoarchean zircons in a metaconglomerate from Jack Hills, Western Australia, has catalyzed intensive study of these zircons and their mineral inclusions, as they represent unique geochemical archives that can be used to unravel the geological evolution of early Earth. Here, we report the occurrence and physical properties of previously undetected CO₂ inclusions that were identified in 3.36–3.47 Ga and 3.80–4.13 Ga zircon grains by confocal micro-Raman spectroscopy. Minimum P–T conditions of zircon formation were determined from the highest density of the inclusions, determined from the density-dependence of the Fermi diad splitting in the Raman spectrum and Ti-in-zircon thermometry. For both age periods, the CO₂ densities and Ti-in-zircon temperatures correspond to high-grade metamorphic conditions (≥5 to ≥7 kbar/~670 to 770 °C) that are typical of mid-crustal regional metamorphism throughout Earth’s history. In addition, fully enclosed, highly disordered graphitic carbon inclusions were identified in two zircon grains from the older population that also contained CO₂ inclusions. Transmission electron microscopy on one of these inclusions revealed that carbon forms a thin amorphous film on the inclusion wall, whereas the rest of the volume was probably occupied by CO₂ prior to analysis. This indicates a close relationship between CO₂ and the reduced carbon inclusions and, in particular that the carbon precipitated from a CO₂-rich fluid, which is inconsistent with the recently proposed biogenic origin of carbon inclusions found in Hadean zircons from Jack Hills.     

Post-injection trapping of mobile CO<Subscript><Emphasis Type="Bold">2</Emphasis></Subscript> in deep aquifers: Assessing the importance of model and parameter uncertainties

Predicting the fate of the injected CO₂ is crucial for the safety of carbon storage operations in deep saline aquifers: especially the evolution of the position, the spreading and the quantity of the mobile CO₂ plume during and after the injection has to be understood to prevent any loss of containment. Fluid flow modelling is challenging not only given the uncertainties on subsurface formation intrinsic properties (parameter uncertainty) but also on the modelling choices/assumptions for representing and numerically implementing the processes occurring when CO₂ displaces the native brine (model uncertainty). Sensitivity analysis is needed to identify the group of factors which contributes the most to the uncertainties in the predictions. In this paper, we present an approach for assessing the importance of model and parameter uncertainties regarding post-injection trapping of mobile CO₂. This approach includes the representation of input parameters, the choice of relevant simulation outputs, the assessment of the mobile plume evolution with a flow simulator and the importance ranking for input parameters. A variance-based sensitivity analysis is proposed, associated with the ACOSSO-like meta-modelling technique to tackle the issues linked with the computational burden posed by the use of long-running simulations and with the different types of uncertainties to be accounted for (model and parameter). The approach is tested on a potential site for CO₂ storage in the Paris basin (France) representative of a project in preliminary stage of development. The approach provides physically sound outcomes despite the challenging context of the case study. In addition, these outcomes appear very helpful for prioritizing the future characterisation efforts and monitoring requirements, and for simplifying the modelling exercise.     

Analytical model for screening potential CO<Subscript>2</Subscript> repositories

Assessing potential repositories for geologic sequestration of carbon dioxide using numerical models can be complicated, costly, and time-consuming, especially when faced with the challenge of selecting a repository from a multitude of potential repositories. This paper presents a set of simple analytical equations (model), based on the work of previous researchers, that could be used to evaluate the suitability of candidate repositories for subsurface sequestration of carbon dioxide. We considered the injection of carbon dioxide at a constant rate into a confined saline aquifer via a fully perforated vertical injection well. The validity of the analytical model was assessed via comparison with the TOUGH2 numerical model. The metrics used in comparing the two models include (1) spatial variations in formation pressure and (2) vertically integrated brine saturation profile. The analytical model and TOUGH2 show excellent agreement in their results when similar input conditions and assumptions are applied in both. The analytical model neglects capillary pressure and the pressure dependence of fluid properties. However, simulations in TOUGH2 indicate that little error is introduced by these simplifications. Sensitivity studies indicate that the agreement between the analytical model and TOUGH2 depends strongly on (1) the residual brine saturation, (2) the difference in density between carbon dioxide and resident brine (buoyancy), and (3) the relationship between relative permeability and brine saturation. The results achieved suggest that the analytical model is valid when the relationship between relative permeability and brine saturation is linear or quasi-linear and when the irreducible saturation of brine is zero or very small.     

The Role of CaCO<Subscript>3</Subscript> Reactions in the Contemporary Oceanic CO<Subscript>2</Subscript> Cycle

The present analysis adjusts previous estimates of global ocean CaCO₃ production rates substantially upward, to 133 × 10¹² mol yr^?1 plankton production and 42 × 10¹² mol yr^?1 shelf benthos production. The plankton adjustment is consistent with recent satellite-based estimates; the benthos adjustment includes primarily an upward adjustment of CaCO₃ production on so-called carbonate-poor sedimentary shelves and secondarily pays greater attention to high CaCO₃ mass (calcimass) and turnover of shelf communities on temperate and polar shelves. Estimated CaCO₃ sediment accumulation rates remain about the same as they have been for some years: ~20 × 10¹² mol yr^?1 on shelves and 11 × 10¹² mol yr^?1 in the deep ocean. The differences between production and accumulation of calcareous materials call for dissolution of ~22 × 10¹² mol yr^?1 (~50 %) of shelf benthonic carbonate production and 122 × 10¹² mol yr^?1 (>90 %) of planktonic production. Most CaCO₃ production, whether planktonic or benthonic, is assumed to take place in water depths of <100 m, while most dissolution is assumed to occur below this depth. The molar ratio of CO₂ release to CaCO₃ precipitation (CO₂↑/CaCO₃↓) is <1.0 and varies with depth. This ratio, Ψ, is presently about 0.66 in surface seawater and 0.85 in ocean waters deeper than about 1000 m. The net flux of CO₂ associated with CaCO₃ reactions in the global ocean in late preindustrial time is estimated to be an apparent influx from the atmosphere to the ocean, of +7 × 10¹² mol C yr^?1, at a time scale of 10²–10³ years. The CaCO₃-mediated influx of CO₂ is approximately offset by CO₂ release from organic C oxidation in the water column. Continuing ocean acidification will have effects on CaCO₃ and organic C metabolic responses to the oceanic inorganic C cycle, although those responses remain poorly quantified.