Geological storage of CO2 in saline aquifers—A review of the experience from existing storage operations

The experience from CO₂ injection at pilot projects (Frio, Ketzin, Nagaoka, US Regional Partnerships) and existing commercial operations (Sleipner, Snøhvit, In Salah, acid-gas injection) demonstrates that CO₂ geological storage in saline aquifers is technologically feasible. Monitoring and verification technologies have been tested and demonstrated to detect and track the CO₂ plume in different subsurface geological environments. By the end of 2008, approximately 20 Mt of CO₂ had been successfully injected into saline aquifers by existing operations. Currently, the highest injection rate and total storage volume for a single storage operation are approximately 1 Mt CO₂/year and 25 Mt, respectively. If carbon capture and storage (CCS) is to be an effective option for decreasing greenhouse gas emissions, commercial-scale storage operations will require orders of magnitude larger storage capacity than accessed by the existing sites. As a result, new demonstration projects will need to develop and test injection strategies that consider multiple injection wells and the optimisation of the usage of storage space. To accelerate large-scale CCS deployment, demonstration projects should be selected that can be readily employed for commercial use; i.e. projects that fully integrate the capture, transport and storage processes at an industrial emissions source.     

Large-scale impact of CO2 storage in deep saline aquifers: A sensitivity study on pressure response in stratified systems

Large volumes of CO₂ captured from carbon emitters (such as coal-fired power plants) may be stored in deep saline aquifers as a means of mitigating climate change. Storing these additional fluids may cause pressure changes and displacement of native brines, affecting subsurface volumes that can be significantly larger than the CO₂ plume itself. This study aimed at determining the three-dimensional region of influence during/after injection of CO₂ and evaluating the possible implications for shallow groundwater resources, with particular focus on the effects of interlayer communication through low-permeability seals. To address these issues quantitatively, we conducted numerical simulations that provide a basic understanding of the large-scale flow and pressure conditions in response to industrial-scale CO₂ injection into a laterally open saline aquifer. The model domain included an idealized multilayered groundwater system, with a sequence of aquifers and aquitards (sealing units) extending from the deep saline storage formation to the uppermost freshwater aquifer. Both the local CO₂-brine flow around the single injection site and the single-phase water flow (with salinity changes) in the region away from the CO₂ plume were simulated. Our simulation results indicate considerable pressure buildup in the storage formation more than 100 km away from the injection zone, whereas the lateral distance migration of brine is rather small. In the vertical direction, the pressure perturbation from CO₂ storage may reach shallow groundwater resources only if the deep storage formation communicates with the shallow aquifers through sealing units of relatively high permeabilities (higher than 10^?18 m²). Vertical brine migration through a sequence of layers into shallow groundwater bodies is extremely unlikely. Overall, large-scale pressure changes appear to be of more concern to groundwater resources than changes in water quality caused by the migration of displaced saline water.     

Assessment of basin-scale hydrologic impacts of CO2 sequestration,Illinois basin

Idealized, basin-scale sharp-interface models of CO₂ injection were constructed for the Illinois basin. Porosity and permeability were decreased with depth within the Mount Simon Formation. Eau Claire confining unit porosity and permeability were kept fixed. We used 726 injection wells located near 42 power plants to deliver 80 million metric tons of CO₂/year. After 100 years of continuous injection, deviatoric fluid pressures varied between 5.6 and 18 MPa across central and southern part of the Illinois basin. Maximum deviatoric pressure reached about 50% of lithostatic levels to the south. The pressure disturbance (>0.03 MPa) propagated 10–25 km away from the injection wells resulting in significant well–well pressure interference. These findings are consistent with single-phase analytical solutions of injection. The radial footprint of the CO₂ plume at each well was only 0.5–2 km after 100 years of injection. Net lateral brine displacement was insignificant due to increasing radial distance from injection well and leakage across the Eau Claire confining unit. On geologic time scales CO₂ would migrate northward at a rate of about 6 m/1000 years. Because of paleo-seismic events in this region (M5.5–M7.5), care should be taken to avoid high pore pressures in the southern Illinois basin.     

Hydraulic characterisation of the Stuttgart formation at the pilot test site for CO2 storage, Ketzin, Germany

The paper presents an approach for the interpretation of hydraulic tests of a CO₂ storage reservoir. The sandstone reservoir is characterised by a fluviatile channel structure embedded in a low-permeability matrix. Pumping tests were carried out in three wells, with simultaneous pressure monitoring in each well.The hydraulic parameters (permeability and storativity) and the boundary configurations were calibrated using three different approaches: (i) parameter calibration and type curve interpretation for single-hole tests, (ii) calibration of the entire build-up phase for cross-hole tests, and (iii) calibration of the initial pressure response for cross-hole pumping tests. In addition, the arrival time of the pressure response was determined and provides additional information about the pathways of hydraulic connection.The measured pumping test permeabilities of the formation were much lower than those measured on the cores, which is very unusual. The pumping test permeabilities are mainly between 50 mD and 100 mD (millidarcy), while core samples show a mean aquifer permeability of 500–1100 mD. Based on this it was concluded that some kind of continuous low-permeability structure exists, which was supported by core material. Three possible aquifer configurations were considered. The first and second were derived from traditional pumping test analysis and were conceptualised using flow boundaries. Each of the analyses provides a different result. A method was developed in which these differences were resolved by interpreting the pressure response with respect to its spatial and temporal sensitivity. This solution lead to a third configuration which was mainly based on spatially-variable permeabilities. Taking into account the pumping test results, the geological background and the behaviour of injected CO₂, we consider only the third configuration to be realistic. The results are in good agreement with modelled CO₂ arrival times and pressure history.     

Carbon dioxide storage capacity in uneconomic coal beds in Alberta,Canada: Methodology,potential and site identification

Methodology is presented for a first-order regional-scale estimation of CO₂ storage capacity in coals under sub-critical conditions, which is subsequently applied to Cretaceous-Tertiary coal beds in Alberta, Canada. Regions suitable for CO₂ storage have been defined on the basis of groundwater depth and CO₂ phase at in situ conditions. The theoretical CO₂ storage capacity was estimated on the basis of CO₂ adsorption isotherms measured on coal samples, and it varies between ∼20 kt CO₂/km² and 1260 kt CO₂/km², for a total of approximately 20 Gt CO₂. This represents the theoretical storage capacity limit that would be attained if there would be no other gases present in the coals or they would be 100% replaced by CO₂, and if all the coals will be accessed by CO₂. A recovery factor of less than 100% and a completion factor less than 50% reduce the theoretical storage capacity to an effective storage capacity of only 6.4 Gt CO₂. Not all the effective CO₂ storage capacity will be utilized because it is uneconomic to build the necessary infrastructure for areas with low storage capacity per unit surface. Assuming that the economic threshold to develop the necessary infrastructure is 200 kt CO₂/km², then the CO₂ storage capacity in coal beds in Alberta is greatly reduced further to a practical capacity of only ∼800 Mt CO₂.     

Geoelectrical methods for monitoring geological CO2 storage: First results from cross-hole and surface–downhole measurements from the CO2SINK test site at Ketzin (Germany)

The feasibility of monitoring CO₂ migration in a saline aquifer at a depth of about 650 m with cross-hole and surface–downhole electrical resistivity tomography (ERT) is investigated at the CO₂SINK test site close to Ketzin (Germany). The permanent vertical electrical resistivity array (VERA) consists of 45 electrodes (15 in the injection well Ktzi201 and 15 in each of the two observation wells Ktzi200 and Ktzi202), successfully placed on the electrically insulated casings, in the depth range of about 590–740 m with a spacing of about 10 m. The three Ketzin wells are arranged as perpendicular triangle with distances of 50 and 100 m.First synthetic modelling studies indicate an increase of the electrical resistivity of about 200% caused by CO₂ injection, corresponding to a bulk CO₂ saturation of 50%, which is in good agreement with laboratory studies. Finite difference inversion of field data delivers three-dimensional resistivity distributions between the wells which are consistent with the reservoir modelling studies.To increase the limited observation area provided by the cross-hole measurements, additional surface–downhole measurements were deployed. A main CO₂ migration in SE–NW direction is deduced from surface to downhole resistivity experiments.The first cross-hole time-lapse results show that the resolution and the coverage of the electrode array in the Ketzin setting are sufficient to resolve the expected resistivity changes on the characteristic length scale of the electrode array. Significant resistivity changes could be measured, however, detailed information on the CO₂ plume could not be resolved yet by VERA under the existing geological circumstances.     

CCS scenarios optimization by spatial multi-criteria analysis: Application to multiple source sink matching in Hebei province

A method, based on spatial analysis of the different criteria to be taken into consideration for building scenarios of CO₂ capture and storage (CCS), has been developed and applied to real case studies in the Hebei province. Totally 88 point sources (42 from power sector, 9 from iron and steel, 18 from cement, 16 from ammonia, and 3 from oil refinery) are estimated and their total emission amounts to 231.7 MtCO₂/year with power, iron and steel, cement, ammonia and oil refinery sharing 59.13%, 25.03%, 11.44%, 3.5%, and 0.91%, respectively. Storage opportunities can be found in Hebei province, characterised by a strong tectonic subsidence during the Tertiary, with several kilometres of accumulated clastic sediments. Carbon storage potential for 25 hydrocarbon fields selected from the Huabei complex is estimated as 215 MtCO₂ with optimistic assumption that all recovered hydrocarbon could be replaced by an equivalent volume of CO₂ at reservoir conditions. Storage potential for aquifers in the Miocene Guantao formation is estimated as 747 MtCO₂ if closed aquifer assumed or 371 MtCO₂ if open aquifer and single highly permeable horizon assumed. Due to poor knowledge on deep hydrogeology and to pressure increase in aquifer, injecting very high rates requested by the major CO₂ sources (>10 MtCO₂/year) is the main challenge, therefore piezometry and discharge must be carefully controlled. A source sink matching model using ArcGIS software is designed to find the least-cost pathway and to estimate transport route and cost accounting for the additional costs of pipeline construction due to landform and land use. Source sink matching results show that only 15–25% of the emissions estimated for the 88 sources can be sequestrated into the hydrocarbon fields and the aquifers if assuming sinks should be able to accommodate at least 15 years of the emissions of a given source.     

Coupled reservoir-geomechanical analysis of CO2 injection and ground deformations at In Salah,Algeria

In Salah Gas Project in Algeria has been injecting 0.5–1 million tonnes CO₂ per year over the past 5 years into a water-filled strata at a depth of about 1800–1900 m. Unlike most CO₂ storage sites, the permeability of the storage formation is relatively low and comparatively thin with a thickness of about 20 m. To ensure adequate CO₂ flow-rates across the low-permeability sand-face, the In Salah Gas Project decided to use long-reach (about 1–1.5 km) horizontal injection wells. In an ongoing research project we use field data and coupled reservoir-geomechanical numerical modeling to assess the effectiveness of this approach and to investigate monitoring techniques to evaluate the performance of a CO₂ injection operation in relatively low-permeability formations. Among the field data used are ground surface deformations evaluated from recently acquired satellite-based inferrometry (InSAR). The InSAR data shows a surface uplift on the order of 5 mm per year above active CO₂ injection wells and the uplift pattern extends several km from the injection wells. In this paper we use the observed surface uplift to constrain our coupled reservoir-geomechanical model and conduct sensitivity studies to investigate potential causes and mechanisms of the observed uplift. The results of our analysis indicate that most of the observed uplift magnitude can be explained by pressure-induced, poro-elastic expansion of the 20-m-thick injection zone, but there could also be a significant contribution from pressure-induced deformations within a 100-m-thick zone of shaly sands immediately above the injection zone.     

CO2SINK—From site characterisation and risk assessment to monitoring and verification: One year of operational experience with the field laboratory for CO2 storage at Ketzin, Germany

The CO₂SINK pilot project at Ketzin is aimed at a better understanding of geological CO₂ storage operation in a saline aquifer. The reservoir consists of fluvial deposits with average permeability ranging between 50 and 100 mDarcy. The main focus of CO₂SINK is developing and testing of monitoring and verification technologies. All wells, one for injection and two for observation, are equipped with smart casings (sensors behind casing, facing the rocks) containing a Distributed Temperature Sensing (DTS) and electrodes for Electrical Resistivity Tomography (ERT). The in-hole Gas Membrane Sensors (GMS) observed the arrival of tracers and CO₂ with high temporal resolution. Geophysical monitoring includes Moving Source Profiling (MSP), Vertical Seismic Profiling (VSP), crosshole, star and 4-D seismic experiments. Numerical models are benchmarked via the monitoring results indicating a sufficient match between observation and prediction, at least for the arrival of CO₂ at the first observation well. Downhole samples of brine showed changes in the fluid composition and biocenosis. First monitoring results indicate anisotropic flow of CO₂ coinciding with the “on-time” arrival of CO₂ at observation well one (Ktzi 200) and the later arrival at observation well two (Ktzi 202). A risk assessment was performed prior to the start of injection. After one year of operations about 18,000 t of CO₂ were injected safely.     

Constraints on the in situ density of CO2 within the Utsira formation from time-lapse seafloor gravity measurements

At Sleipner, CO₂ is being separated from natural gas and injected into an underground saline aquifer for environmental purposes. Uncertainty in the aquifer temperature leads to uncertainty in the in situ density of CO₂. In this study, gravity measurements were made over the injection site in 2002 and 2005 on top of 30 concrete benchmarks on the seafloor in order to constrain the in situ CO₂ density. The gravity measurements have a repeatability of 4.3 μGal for 2003 and 3.5 μGal for 2005. The resulting time-lapse uncertainty is 5.3 μGal. Unexpected benchmark motions due to local sediment scouring contribute to the uncertainty. Forward gravity models are calculated based on both 3D seismic data and reservoir simulation models. The time-lapse gravity observations best fit a high temperature forward model based on the time-lapse 3D seismics, suggesting that the average in situ CO₂ density is about to 530 kg/m³. Uncertainty in determining the average density is estimated to be ±65 kg/m³ (95% confidence), however, this does not include uncertainties in the modeling. Additional seismic surveys and future gravity measurements will put better constraints on the CO₂ density and continue to map out the CO₂ flow.     

Assessing the economic feasibility of regional deep saline aquifer CO2 injection and storage: A geomechanics-based workflow applied to the Rose Run sandstone in Eastern Ohio, USA

Typical top-down regional assessments of CO₂ storage feasibility are sufficient for determining the maximum volumetric capacity of deep saline aquifers. However, they do not reflect the regional economic feasibility of storage. This is controlled, in part, by the number and type of injection wells that are necessary to achieve regional CO₂ storage goals. In contrast, the geomechanics-based assessment workflow that we present in this paper follows a bottom-up approach for evaluating regional deep saline aquifer CO₂ storage feasibility. The CO₂ storage capacity of an aquifer is a function of its porous volume as well as its CO₂ injectivity. For a saline aquifer to be considered feasible in this assessment it must be able to store a specified amount of CO₂ at a reasonable cost per ton of CO₂. The proposed assessment workflow has seven steps that include (1) defining the storage project and goals, (2) characterizing the geology and developing a geomechanical model of the aquifer, (3) constructing 3D aquifer models, (4) simulating CO₂ injection, (5,6) evaluating CO₂ injection and storage feasibility (with and without injection well stimulation), and (7) determining whether it is economically feasible to proceed with the storage project. The workflow was applied to a case study of the Rose Run sandstone aquifer in the Eastern Ohio River Valley, USA. We found that it is feasible in this region to inject 113 Mt CO₂/year for 30 years at an associated well cost of less than US $1.31/t CO₂, but only if injectivity enhancement techniques such as hydraulic fracturing and injection induced micro-seismicity are implemented.     

3D geomechanical modelling for CO2 geologic storage in the Dogger carbonates of the Paris Basin

CO₂ injection into a depleted hydrocarbon field or aquifer may give rise to a variety of coupled physical and chemical processes. During CO₂ injection, the increase in pore pressure can induce reservoir expansion. As a result the in situ stress field may change in and around the reservoir. The geomechanical behaviour induced by oil production followed by CO₂ injections into an oil field reservoir in the Paris Basin has been numerically modelled. This paper deals with an evaluation of the induced deformations and in situ stress changes, and their potential effects on faults, using a 3D geomechanical model. The geomechanical analysis of the reservoir–caprock system was carried out as a feasibility study using pressure information in a “one way” coupling, where pressures issued from reservoir simulations were integrated as input for a geomechanical model. The results show that under specific assumptions the mechanical effects of CO₂ injection do not affect the mechanical stability of the reservoir–caprock system. The ground vertical movement at the surface ranges from ?2 mm during oil production to +2.5 mm during CO₂ injection. Furthermore, the changes in in situ stresses predicted under specific assumptions by geomechanical modelling are not significant enough to jeopardize the mechanical stability of the reservoir and caprock. The stress changes issued from the 3D geomechanical modelling are also combined with a Mohr–Coulomb analysis to determine the fault slip tendency. By integrating the stress changes issued from the geomechanical modelling into the fault stability analysis, the critical pore pressure for fault reactivation is higher than calculated for the fault stability analysis considering constant horizontal stresses.     

Screening CO2 storage options in The Netherlands

This paper describes the development and application of a methodology to screen and rank Dutch reservoirs suitable for long-term large scale CO₂ storage. The screening focuses on off- and on-shore individual aquifers, gas and oil fields. In total 176 storage reservoirs have been taken into consideration: 138 gas fields, 4 oil fields and 34 aquifers, with a total theoretical storage potential of about 3200 Mt CO₂. The reservoirs are screened according to three criteria: potential storage capacity, storage costs and effort needed to manage risk. Due to the large number of reservoirs, which limits the possibility to use any pair-wise comparison method (e.g. Multi-Criteria programs such as Bosda or Naiade), a spreadsheet tool was designed to provide an assessment of each of the criteria through an evaluation of the fields present in the database and a set of scores provided by a (inter)national panel of experts. The assessment is sufficiently simple and allows others to review it, re-do it or expand it. The results of the methodology show that plausible comparisons of prospective sites with limited characterization data are possible.     

Pathways towards large-scale implementation of CO2 capture and storage: A case study for the Netherlands

We sketch four possible pathways how carbon dioxide capture and storage (CCS) (r)evolution may occur in the Netherlands, after which the implications in terms of CO₂ stored and avoided, costs and infrastructural requirements are quantified. CCS may play a significant role in decarbonising the Dutch energy and industrial sector, which currently emits nearly 100 Mt CO₂/year. We found that 15 Mt CO₂ could be avoided annually by 2020, provided some of the larger gas fields that become available the coming decade could be used for CO₂ storage. Halfway this century, the mitigation potential of CCS in the power sector, industry and transport fuel production is estimated at maximally 80–110 Mt CO₂/year, of which 60–80 Mt CO₂/year may be avoided at costs between 15 and 40 €/t CO₂, including transport and storage. Avoiding 30–60 Mt CO₂/year by means of CCS is considered realistic given the storage potential represented by Dutch gas fields, although it requires planning to assure that domestic storage capacity could be used for CO₂ storage. In an aggressive climate policy, avoiding another 50 Mt CO₂/year may be possible provided that nearly all capture opportunities that occur are taken. Storing such large amounts of CO₂ would only be possible if the Groningen gas field or large reservoirs in the British or Norwegian part of the North Sea will become available.     

Modelling of CO2 arrival time at Ketzin – Part I

The injection of CO₂ at the Ketzin storage site and the chemical detection of its arrival in the observation well allowed testing of different numerical simulation codes. ECLIPSE 100 (E100, black-oil simulator), ECLIPSE 300 (E300, compositional CO₂STORE) and MUFTE-UG were used for predictive modelling applying a constant injection rate of 1 kg s^?1 CO₂ and for a history match applying the actual variable injection rate which ranged from 0 to 0.7 kg s^?1 and averaged 0.23 kg s^?1. The geological model applied, is based on all available geophysical and geological information and has been the same for all programs.The results of the constant injection regime show a good agreement among the programs with a discrepancy of 21–33% for the CO₂ arrival times. However, it is determined from the comparison of the cumulative mass of CO₂ at the time of CO₂ arrival that the injection regime is an important factor for the accurate prediction of CO₂ migration within a saline aquifer. Comparing the actual variable injection regime with the simulations applying a constant injection rate the results are relatively inaccurate.Regarding the actual variable injection regime, which was evaluated using all three simulators, the computational results show a good agreement with the data actually measured at the first observation well. Here, the calculated arrival times exceeded the actual ones by 8.1% (E100), 9.2% (E300) and 17.7% (MUFTE-UG).It can be concluded that irrespective of the deviations of the simulations, due to combinations of different codes and slight differences in input parameters, all three programs are well equipped to give a reliable estimate of the arrival of CO₂. Deviations in the results mainly occur due to different input data and grid size choices done by the different modelling teams working independently of each other. Deviations of the simulations results compared to the actual CO₂ arrival time result from uncertainties in the implementation of the geological model, which was set up based on well log data and analogue studies.     

Biofilm enhanced geologic sequestration of supercritical CO2

In order to develop subsurface CO₂ storage as a viable engineered mechanism to reduce the emission of CO₂ into the atmosphere, any potential leakage of injected supercritical CO₂ (SC-CO₂) from the deep subsurface to the atmosphere must be reduced. Here, we investigate the utility of biofilms, which are microorganism assemblages firmly attached to a surface, as a means of reducing the permeability of deep subsurface porous geological matrices under high pressure and in the presence of SC-CO₂, using a unique high pressure (8.9 MPa), moderate temperature (32 °C) flow reactor containing 40 millidarcy Berea sandstone cores. The flow reactor containing the sandstone core was inoculated with the biofilm forming organism Shewanella fridgidimarina. Electron microscopy of the rock core revealed substantial biofilm growth and accumulation under high-pressure conditions in the rock pore space which caused >95% reduction in core permeability. Permeability increased only slightly in response to SC-CO₂ challenges of up to 71 h and starvation for up to 363 h in length. Viable population assays of microorganisms in the effluent indicated survival of the cells following SC-CO₂ challenges and starvation, although S. fridgidimarina was succeeded by Bacillus mojavensis and Citrobacter sp. which were native in the core. These observations suggest that engineered biofilm barriers may be used to enhance the geologic sequestration of atmospheric CO₂.     

Geomechanical analysis of the Naylor Field,Otway Basin,Australia: Implications for CO2 injection and storage

A geomechanical assessment of the Naylor Field, Otway Basin, Australia has been undertaken to investigate the possible geomechanical effects of CO₂ injection and storage. The study aims to evaluate the geomechanical behaviour of the caprock/reservoir system and to estimate the risk of fault reactivation. The stress regime in the onshore Victorian Otway Basin is inferred to be strike–slip if the maximum horizontal stress is calculated using frictional limits and DITF (drilling induced tensile fracture) occurrence, or normal if maximum horizontal stress is based on analysis of dipole sonic log data. The NW–SE maximum horizontal stress orientation (142°N) determined from a resistivity image log is broadly consistent with previous estimates and confirms a NW–SE maximum horizontal stress orientation for the Otway Basin.An analytical geomechanical solution is used to describe stress changes in the subsurface of the Naylor Field. The computed reservoir stress path for the Naylor Field is then incorporated into fault reactivation analysis to estimate the minimum pore pressure increase required to cause fault reactivation (ΔP_p).The highest reactivation propensity (for critically-oriented faults) ranges from an estimated pore pressure increase (ΔP_p) of 1 MPa to 15.7 MPa (estimated pore pressure of 18.5–33.2 MPa) depending on assumptions made about maximum horizontal stress magnitude, fault strength, reservoir stress path and Biot's coefficient. The critical pore pressure changes for known faults at Naylor Field range from an estimated pore pressure increase (ΔP_p) of 2 MPa to 17 MPa (estimated pore pressure of 19.5–34.5 MPa).     

Numerical investigation concerning the impact of CO2 geologic storage on regional groundwater flow

Large-scale storage of carbon dioxide in saline aquifers may cause considerable pressure perturbation and brine migration in deep rock formations, which may have a significant influence on the regional groundwater system. With the help of parallel computing techniques, we conducted a comprehensive, large-scale numerical simulation of CO₂ geologic storage that predicts not only CO₂ migration, but also its impact on regional groundwater flow. As a case study, a hypothetical industrial-scale CO₂ injection in Tokyo Bay, which is surrounded by the most heavily industrialized area in Japan, was considered, and the impact of CO₂ injection on near-surface aquifers was investigated, assuming relatively high seal-layer permeability (higher than 10 microdarcy). A regional hydrogeological model with an area of about 60 km × 70 km around Tokyo Bay was discretized into about 10 million gridblocks. To solve the high-resolution model efficiently, we used a parallelized multiphase flow simulator TOUGH2-MP/ECO2N on a world-class high performance supercomputer in Japan, the Earth Simulator. In this simulation, CO₂ was injected into a storage aquifer at about 1 km depth under Tokyo Bay from 10 wells, at a total rate of 10 million tons/year for 100 years. Through the model, we can examine regional groundwater pressure buildup and groundwater migration to the land surface. The results suggest that even if containment of CO₂ plume is ensured, pressure buildup on the order of a few bars can occur in the shallow confined aquifers over extensive regions, including urban inlands.     

Integration of counter-current relative permeability in the simulation of CO2 injection into saline aquifers

Carbon dioxide (CO₂) injection into saline aquifers is one of the promising options to sequester large amounts of CO₂ in geological formations. During as well as after injection of CO₂ into an aquifer, CO₂ migrates towards the top of the formation due to density differences between the formation brine and the injected CO₂. The time scales of CO₂ migration towards the top of an aquifer and the fraction of CO₂ that is trapped as residual gas depends strongly on the driving forces that are acting on the injected CO₂.When CO₂ migrates to the top of an aquifer, brine may be displaced downwards in a counter-current flow setting particularly during the injection period. A majority of the published work on counter-current flow settings have reported significant reductions in the associated relative permeability functions as compared to co-current measurements. However, this phenomenon has not yet been considered in the simulation of CO₂ storage into saline aquifers.In this paper we study the impact of changes in mobility for the two-phase brine/CO₂ system as a result of transitions between co- and counter-current flow settings. We have included this effect in a simulator and studied the impact of the related mobility reduction on the saturation distribution and residual saturation of CO₂ in aquifers over relevant time scales. We demonstrate that the reduction in relative permeability in the vertical direction changes the plume migration pattern and has an impact on the amount of gas that is trapped as a function of time. This is to our best knowledge the first attempt to integrate counter-current relative permeability into the simulation of injection and subsequent migration of CO₂ in aquifers. The results and analysis presented in this paper are directly relevant to all ongoing activities related to the design of large-scale CO₂ storage in saline aquifers.     

Analytical solution for estimating storage efficiency of geologic sequestration of CO2

During injection of carbon dioxide (CO₂) into deep saline aquifers, the available pore volume of the aquifer may be used inefficiently, thereby decreasing the effective capacity of the repository for CO₂ storage. Storage efficiency is the fraction of the available pore space that is utilized for CO₂ storage, or, in other words, it is the ratio between the volume of stored CO₂ and the maximum available pore volume. In this note, we derive and present simple analytical expressions for estimating CO₂ storage efficiency under the scenario of a constant-rate injection of CO₂ into a confined, homogeneous, isotropic, saline aquifer. The expressions for storage efficiency are derived from models developed previously by other researchers describing the shape of the CO₂-brine interface. The storage efficiency of CO₂ is found to depend on three dimensionless groups, namely: (1) the residual saturation of brine after displacement by CO₂; (2) the ratio of CO₂ mobility to brine mobility; (3) a dimensionless group (which we call a “gravity factor”) that quantifies the importance of CO₂ buoyancy relative to CO₂ injection rate. In the particular case of negligible residual brine saturation and negligible buoyancy effects, the storage efficiency is approximately equal to the ratio of the CO₂ viscosity to the brine viscosity. Storage efficiency decreases as the gravity factor increases, because the buoyancy of the CO₂ causes it to occupy a thin layer at the top of the confined formation, while leaving the lower part of the aquifer under-utilized. Estimates of storage efficiency from our simple analytical expressions are in reasonable agreement with values calculated from simulations performed with more complicated multi-phase-flow simulation software. Therefore, we suggest that the analytical expressions presented herein could be used as a simple and rapid tool to screen the technical or economic feasibility of a proposed CO₂ injection scenario.