CO<Subscript>2</Subscript> injection into saline carbonate aquifer formations I: laboratory investigation

Although there are a number of mathematical modeling studies for carbon dioxide (CO₂) injection into aquifer formations, experimental studies are limited and most studies focus on injection into sandstone reservoirs as opposed to carbonate ones. This study presents the results of computerized tomography (CT) monitored laboratory experiments to analyze permeability and porosity changes as well as to characterize relevant chemical reactions associated with injection and storage of CO₂ in carbonate formations. CT monitored experiments are designed to model fast near well bore flow and slow reservoir flows. Highly heterogeneous cores drilled from a carbonate aquifer formation located in South East Turkey were used during the experiments. Porosity changes along the core plugs and the corresponding permeability changes are reported for different CO₂ injection rates and different salt concentrations of formation water. It was observed that either a permeability increase or a permeability reduction can be obtained. The trend of change in rock properties is very case dependent because it is related to distribution of pores, brine composition and thermodynamic conditions. As the salt concentration decreases, porosity and the permeability decreases are less pronounced. Calcite deposition is mainly influenced by orientation, with horizontal flow resulting in larger calcite deposition compared to vertical flow.     

Linear Stability Analysis of Double-Diffusive Convection in Porous Media,with Application to Geological Storage of CO<Subscript>2</Subscript>

Onset of double-diffusive buoyancy-driven flow resulted from vertical temperature and concentration gradients in a horizontal layer of a saturated and homogenous porous medium is investigated using amplification factor theory. After injection of CO₂ into a deep saline aquifer, the density of the brine saturated with CO₂ increases slightly. This increase in density induces natural convection. The effect of geothermal gradient is also considered in this work as a second incentive for convection and the double-diffusion convection was studied. Linear stability analysis is used to predict the inception of instabilities and initial wavelength of the convective instabilities. The analysis presented is applied to acid gas injection (as an analogue for CO₂ storage) into saline aquifers in the Alberta basin. It is found that the geothermal gradient does not have significant effect on the onset of convection for these aquifers. It is shown that the geothermal effects on the onset of natural convection are negligible as compared to the solutal effects induced by dissolution and diffusion of CO₂ in deep saline aquifers. Therefore, the linear stability analysis and the long-term numerical simulation of CO₂ sequestration into such saline aquifers may be assumed to be isothermal in terms of natural convection occurrence.     

Numerical Simulations of the Thermal Impact of Supercritical CO<Subscript>2</Subscript> Injection on Chemical Reactivity in a Carbonate Saline Reservoir

Geological sequestration of CO₂ offers a promising solution for reducing net emissions of greenhouse gases into the atmosphere. This emerging technology must make it possible to inject CO₂ into deep saline aquifers or oil- and gas-depleted reservoirs in the supercritical state (P > 7.4MPa and T > 31.1°C) to achieve a higher density and therefore occupy less volume underground. Previous experimental and numerical simulations have demonstrated that massive CO₂ injection in saline reservoirs causes a major disequilibrium of the physical and geochemical characteristics of the host aquifer. The near-well injection zone seems to constitute an underground hydrogeological system particularly impacted by supercritical CO₂ injection and the most sensitive area, where chemical phenomena (e.g. mineral dissolution/precipitation) can have a major impact on the porosity and permeability. Furthermore, these phenomena are highly sensitive to temperature. This study, based on numerical multi-phase simulations, investigates thermal effects during CO₂ injection into a deep carbonate formation. Different thermal processes and their influence on the chemical and mineral reactivity of the saline reservoir are discussed. This study underlines both the minor effects of intrinsic thermal and thermodynamic processes on mineral reactivity in carbonate aquifers, and the influence of anthropic thermal processes (e.g. injection temperature) on the carbonates’ behaviour.     

Characteristics of Salt-Precipitation and the Associated Pressure Build-Up during CO<Subscript>2</Subscript> Storage in Saline Aquifers

Mitigation and control of borehole pressure at the bottom of an injection well is directly related to the effective management of well injectivity during geologic carbon sequestration activity. Researchers have generally accepted the idea that high rates of CO₂ injection into low permeability strata results in increased bottom-hole pressure in a well. However, the results of this study suggested that this is not always the case, due to the occurrence of localized salt precipitation adjacent to the injection well. A series of numerical simulations indicated that in some cases, a low rate of CO₂ injection into high permeability formation induced greater pressure build-up. This occurred because of the different types of salt precipitation pattern controlled by buoyancy-driven CO₂ plume migration. The first type is non-localized salt precipitation, which is characterized by uniform salt precipitation within the dry-out zone. The second type, localized salt precipitation, is characterized by an abnormally high level of salt precipitation at the dry-out front. This localized salt precipitation acts as a barrier that hampers the propagation of both CO₂ and pressure to the far field as well as counter-flowing brine migration toward the injection well. These dynamic processes caused a drastic pressure build-up in the well, which decreased injectivity. By modeling a series of test cases, it was found that low-rate CO₂ injection into high permeability formation was likely to cause localized salt precipitation. Sensitivity studies revealed that brine salinity linearly affected the level of salt precipitation, and that vertical permeability enhanced the buoyancy effect which increased the growth of the salt barrier. The porosity also affected both the level of localized salt precipitation and dry-out zone extension depending on injection rates. High temperature injected CO₂ promoted the vertical movement of the CO₂ plume, which accelerated localized salt precipitation, but at the same time caused a decrease in the density of the injected CO₂. The combination of these two effects eventually decreased bottomhole pressure. Considering the injectivity degradation, a method is proposed for decreasing the pressure build-up and increasing injectivity by assigning a ‘skin zone’ that represents a local region with a transmissivity different from that of the surrounding aquifer.     

Effect of Well Orientation (Vertical vs. Horizontal) and Well Length on the Injection of CO<Subscript>2</Subscript> in Deep Saline Aquifers

Simulations of CO₂ injection into confined saline aquifers were conducted for both vertical and horizontal injection wells. The metrics used in quantifying the performances of different injection scenarios included changes in pressure near the injection well, mass of CO₂ dissolved into brine (solubility trapping), and storage efficiency, all evaluated with an assumed injection period of 50 years. Metrics were quantified as functions of well length, well orientation, CO₂ injection rate, and formation anisotropy (ratio of vertical to horizontal conductivity). When equal well lengths are compared, there is not a significant difference between the predicted performances of horizontal and vertical wells. However, the length of a horizontal well may exceed the length of a vertical well because the length of the horizontal well is not constrained to the vertical thickness of the geologic formation. Simulations show that, as the length of the horizontal well is allowed to increase, the geologic formation can receive a significantly higher CO₂ injection rate without exceeding a maximum allowable pressure. This result is observed in both isotropic and anisotropic formations, and suggests that horizontal wells may be advantageous under pressure-limited conditions. However, the use of horizontal wells does not significantly improve the storage efficiency, and under strongly anisotropic conditions, a vertical well provides higher storage efficiency than a horizontal well. We conclude that horizontal wells may be preferable if the goal is to sequester a large amount of CO₂ in a short period of time, but do not offer a significant advantage in terms of long-term capacity of a potential repository.     

An Experimental Study of CO<Subscript>2</Subscript> Exsolution and Relative Permeability Measurements During CO<Subscript>2</Subscript> Saturated Water Depressurization

Dissolution of CO₂ into brine is an important and favorable trapping mechanism for geologic storage of CO₂. There are scenarios, however, where dissolved CO₂ may migrate out of the storage reservoir. Under these conditions, CO₂ will exsolve from solution during depressurization of the brine, leading to the formation of separate phase CO₂. For example, a CO₂ sequestration system with a brine-permeable caprock may be favored to allow for pressure relief in the sequestration reservoir. In this case, CO₂-rich brine may be transported upwards along a pressure gradient caused by CO₂ injection. Here we conduct an experimental study of CO₂ exsolution to observe the behavior of exsolved gas under a wide range of depressurization. Exsolution experiments in highly permeable Berea sandstones and low permeability Mount Simon sandstones are presented. Using X-ray CT scanning, the evolution of gas phase CO₂ and its spatial distribution is observed. In addition, we measure relative permeability for exsolved CO₂ and water in sandstone rocks based on mass balances and continuous observation of the pressure drop across the core from 12.41 to 2.76 MPa. The results show that the minimum CO₂ saturation at which the exsolved CO₂ phase mobilization occurs is from 11.7 to 15.5%. Exsolved CO₂ is distributed uniformly in homogeneous rock samples with no statistical correlation between porosity and CO₂ saturation observed. No gravitational redistribution of exsolved CO₂ was observed after depressurization, even in the high permeability core. Significant differences exist between the exsolved CO₂ and water relative permeabilities, compared to relative permeabilities derived from steady-state drainage relative permeability measurements in the same cores. Specifically, very low CO₂ and water relative permeabilities are measured in the exsolution experiments, even when the CO₂ saturation is as high as 40%. The large relative permeability reduction in both the water and CO₂ phases is hypothesized to result from the presence of disconnected gas bubbles in this two-phase flow system. This feature is also thought to be favorable for storage security after CO₂ injection.     

Influence of Permeability Heterogeneity on Miscible CO<Subscript>2</Subscript> Flooding Efficiency in Sandstone Reservoirs: An Experimental Investigation

In this study, we systematically investigate the effect of core-scale heterogeneity on the performance of miscible CO₂ flooding under various injection modes (secondary and tertiary). Manufactured heterogeneous core plugs are used to simulate vertical and horizontal heterogeneity that may be present in a reservoir. A sample with vertical heterogeneity (i.e. a layered sample) is constructed using two axially cut half plugs each with a distinctly different permeability value. In these samples, the permeability ratio (PR) defines the ratio between the permeabilities of adjacent half plugs. Horizontal heterogeneity (i.e. a composite sample) is introduced by stacking two or three short cylindrical core segments each with a different permeability value. Our special sample construction techniques have also enabled us to investigate the effect of permeability ratio and crossflow in layered samples and axial arrangement of core segments in composite samples on the ultimate recovery of the floods. Core flooding experiments are conducted with an n-Decane–brine–CO₂ system at a pore pressure of 17.2 MPa and a temperature of 343 K. At this temperature, the minimum miscibility pressure of CO₂ with n-Decane is 12.6–12.7 MPa so it is expected that at 17.2 MPa CO₂ is fully miscible with n-Decane. The results obtained for both the composite and layered samples indicate that CO₂ injection would achieve the highest recovery factor (RF) when performed under the secondary mode (e.g. layered: 79.00%, composite: 89.83%) compared with the tertiary mode (e.g. layered: 73.2%, composite: 86.2%). This may be attributed to the effect of water shielding which impedes the access of the injected CO₂ to the residual oil under the tertiary injection mode. It is also found that the oil recovery from a layered sample decreases noticeably with an increase in the PR as higher PR makes the displacement more uneven due to CO₂ channelling. The RFs of 93.4, 87.89, 77.9 and 69.8% correspond to PRs of 1, 2.5, 5, and 12.5, respectively. In addition, for the layered samples, crossflow was found to have an important role during the recovery process; however, due to excessive channelling, this effect tends to diminish as PR increases. Compared with the layered heterogeneity, the effect of composite heterogeneity on the RF seems to be very subtle as the RF is found to be almost independent from the permeability sequence along the length of a composite sample. This outcome may have been caused by the small diameter of the plugs resulting in invariable 1-D floods.     

Modeling Gas Transport in the Shallow Subsurface During the ZERT CO<Subscript>2</Subscript> Release Test

We used the multiphase and multicomponent TOUGH2/EOS7CA model to carry out predictive simulations of CO₂ injection into the shallow subsurface of an agricultural field in Bozeman, Montana. The purpose of the simulations was to inform the choice of CO₂ injection rate and design of monitoring and detection activities for a CO₂ release experiment. The release experiment configuration consists of a long horizontal well (70 m) installed at a depth of approximately 2.5 m into which CO₂ is injected to mimic leakage from a geologic carbon sequestration site through a linear feature such as a fault. We estimated the permeability of the soil and cobble layers present at the site by manual inversion of measurements of soil CO₂ flux from a vertical-well CO₂ release. Based on these estimated permeability values, predictive simulations for the horizontal well showed that CO₂ injection just below the water table creates an effective gas-flow pathway through the saturated zone up to the unsaturated zone. Once in the unsaturated zone, CO₂ spreads out laterally within the cobble layer, where liquid saturation is relatively low. CO₂ also migrates upward into the soil layer through the capillary barrier and seeps out at the ground surface. The simulations predicted a breakthrough time of approximately two days for the 100kg d⁻¹ injection rate, which also produced a flux within the range desired for testing detection and monitoring approaches. The seepage area produced by the model was approximately five meters wide above the horizontal well, compatible with the detection and monitoring methods tested. For a given flow rate, gas-phase diffusion of CO₂ tends to dominate over advection near the ground surface, where the CO₂ concentration gradient is large, while advection dominates deeper in the system.     

Pressure Buildup During CO<Subscript>2</Subscript> Injection into a Closed Brine Aquifer

CO₂ injected into porous formations is accommodated by reduction in the volume of the formation fluid and enlargement of the pore space, through compression of the formation fluids and rock material, respectively. A critical issue is how the resulting pressure buildup will affect the mechanical integrity of the host formation and caprock. Building on an existing approximate solution for formations of infinite radial extent, this article presents an explicit approximate solution for estimating pressure buildup due to injection of CO₂ into closed brine aquifers of finite radial extent. The analysis is also applicable for injection into a formation containing multiple wells, in which each well acts as if it were in a quasi-circular closed region. The approximate solution is validated by comparison with vertically averaged results obtained using TOUGH2 with ECO2N (where many of the simplifying assumptions are relaxed), and is shown to be very accurate over wide ranges of the relevant parameter space. The resulting equations for the pressure distribution are explicit, and can be easily implemented within spreadsheet software for estimating CO₂ injection capacity.     

Investigation of CO<Subscript>2</Subscript> Plume Behavior for a Large-Scale Pilot Test of Geologic Carbon Storage in a Saline Formation

The hydrodynamic behavior of carbon dioxide (CO₂) injected into a deep saline formation is investigated, focusing on trapping mechanisms that lead to CO₂ plume stabilization. A numerical model of the subsurface at a proposed power plant with CO₂ capture is developed to simulate a planned pilot test, in which 1,000,000 metric tons of CO₂ is injected over a 4-year period, and the subsequent evolution of the CO₂ plume for hundreds of years. Key measures are plume migration distance and the time evolution of the partitioning of CO₂ between dissolved, immobile free-phase, and mobile free-phase forms. Model results indicate that the injected CO₂ plume is effectively immobilized at 25 years. At that time, 38% of the CO₂ is in dissolved form, 59% is immobile free phase, and 3% is mobile free phase. The plume footprint is roughly elliptical, and extends much farther up-dip of the injection well than down-dip. The pressure increase extends far beyond the plume footprint, but the pressure response decreases rapidly with distance from the injection well, and decays rapidly in time once injection ceases. Sensitivity studies that were carried out to investigate the effect of poorly constrained model parameters permeability, permeability anisotropy, and residual CO₂ saturation indicate that small changes in properties can have a large impact on plume evolution, causing significant trade-offs between different trapping mechanisms.     

Injection and Storage of CO<Subscript>2</Subscript> in Deep Saline Aquifers: Analytical Solution for CO<Subscript>2</Subscript> Plume Evolution During Injection

Injection of fluids into deep saline aquifers is practiced in several industrial activities, and is being considered as part of a possible mitigation strategy to reduce anthropogenic emissions of carbon dioxide into the atmosphere. Injection of CO₂ into deep saline aquifers involves CO₂ as a supercritical fluid that is less dense and less viscous than the resident formation water. These fluid properties lead to gravity override and possible viscous fingering. With relatively mild assumptions regarding fluid properties and displacement patterns, an analytical solution may be derived to describe the space–time evolution of the CO₂ plume. The solution uses arguments of energy minimization, and reduces to a simple radial form of the Buckley–Leverett solution for conditions of viscous domination. In order to test the applicability of the analytical solution to the CO₂ injection problem, we consider a wide range of subsurface conditions, characteristic of sedimentary basins around the world, that are expected to apply to possible CO₂ injection scenarios. For comparison, we run numerical simulations with an industry standard simulator, and show that the new analytical solution matches a full numerical solution for the entire range of CO₂ injection scenarios considered. The analytical solution provides a tool to estimate practical quantities associated with CO₂ injection, including maximum spatial extent of a plume and the shape of the overriding less-dense CO₂ front.     

Effects of CO<Subscript>2</Subscript> Compressibility on CO<Subscript>2</Subscript> Storage in Deep Saline Aquifers

The injection of supercritical CO₂ in deep saline aquifers leads to the formation of a CO₂ plume that tends to float above the formation brine. As pressure builds up, CO₂ properties, i.e. density and viscosity, can vary significantly. Current analytical solutions do not account for CO₂ compressibility. In this article, we investigate numerically and analytically the effect of this variability on the position of the interface between the CO₂-rich phase and the formation brine. We introduce a correction to account for CO₂ compressibility (density variations) and viscosity variations in current analytical solutions. We find that the error in the interface position caused by neglecting CO₂ compressibility is relatively small when viscous forces dominate. However, it can become significant when gravity forces dominate, which is likely to occur at late times of injection.     

Impact of CO<Subscript>2</Subscript> Injection on Flow Behavior of Coalbed Methane Reservoirs

On the basis of observations at four enhanced coalbed methane (ECBM)/CO₂ sequestration pilots, a laboratory-scale study was conducted to understand the flow behavior of coal in a methane/CO₂ environment. Sorption-induced volumetric strain was first measured by flooding fresh coal samples with adsorptive gases (methane and CO₂). In order to replicate the CO₂–ECBM process, CO₂ was then injected into a methane-saturated core to measure the incremental “swelling.” As a separate effort, the permeability of a coal core, held under triaxial stress, was measured using methane. This was followed by CO₂ flooding to replace the methane. In order to best replicate the conditions in situ, the core was held under uniaxial strain, that is, no horizontal strain was permitted during CO₂ flooding. Instead, the horizontal stress was adjusted to ensure zero strain. The results showed that the relative strain ratio for CO₂/methane was between 2 and 3.5. The measured volumetric strains were also fitted using a Langmuir-type model, thus enabling calculation of the strain at any gas pressure and using the analytical permeability models. For permeability work, effort was made to increase the horizontal stress to achieve the desired zero horizontal strain condition expected under in situ condition, but this became impossible because the “excess” stress required to maintain this condition was very large, resulting in sample failure. Finally, when CO₂ was introduced and horizontal strain was permitted, permeability reduction was an order of magnitude greater, suggesting that the “excess” stress would have reduced it significantly further. The positive finding of the work was that the “excess” stresses associated with injection of CO₂ are large. The excess stresses generated might be sufficient to cause microfracturing and increased permeability, and improved injectivity. Also, there might be a weakening effect resulting from repeated CO₂ injection, as has been found to be the case with thermal cycling of rocks.     

Numerical Simulation Studies of the Long-term Evolution of a CO<Subscript>2</Subscript> Plume in a Saline Aquifer with a Sloping Caprock

We have used the TOUGH2-MP/ECO2N code to perform numerical simulation studies of the long-term behavior of CO₂ stored in an aquifer with a sloping caprock. This problem is of great practical interest, and is very challenging due to the importance of multi-scale processes. We find that the mechanism of plume advance is different from what is seen in a forced immiscible displacement, such as gas injection into a water-saturated medium. Instead of pushing the water forward, the plume advances because the vertical pressure gradients within the plume are smaller than hydrostatic, causing the groundwater column to collapse ahead of the plume tip. Increased resistance to vertical flow of aqueous phase in anisotropic media leads to reduced speed of up-dip plume advancement. Vertical equilibrium models that ignore effects of vertical flow will overpredict the speed of plume advancement. The CO₂ plume becomes thinner as it advances, but the speed of advancement remains constant over the entire simulation period of up to 400 years, with migration distances of more than 80 km. Our simulations include dissolution of CO₂ into the aqueous phase and associated density increase, and molecular diffusion. However, no convection develops in the aqueous phase because it is suppressed by the relatively coarse (sub-) horizontal gridding required in a regional-scale model. A first crude sub-grid-scale model was developed to represent convective enhancement of CO₂ dissolution. This process is found to greatly reduce the thickness of the CO₂ plume, but, for the parameters used in our simulations, does not affect the speed of plume advancement.     

Modeling of Rocks and Cement Alteration due to CO<Subscript>2</Subscript> Injection in an Exploited Gas Reservoir

The injection of CO₂ in exploited natural gas reservoirs as a means to reduce greenhouse gas (GHG) emissions is highly attractive as it takes place in well-known geological structures of proven integrity with respect to gas leakage. The injection of a reactive gas such as CO₂ puts emphasis on the possible alteration of reservoir and caprock formations and especially of the wells’ cement sheaths induced by the modification of chemical equilibria. Such studies are important for injectivity assurance, wellbore integrity, and risk assessment required for CO₂ sequestration site qualification. Within a R&D project funded by Eni, we set up a numerical model to investigate the rock–cement alterations driven by the injection of CO₂ into a depleted sweet natural gas pool. The simulations are performed with the TOUGHREACT simulator (Xu et al. in Comput Geosci 32:145–165, 2006) coupled to the TMGAS EOS module (Battistelli and Marcolini in Int J Greenh Gas Control 3:481–493, 2009) developed for the TOUGH2 family of reservoir simulators (Pruess et al. in TOUGH2 User’s Guide, Version 2.0, 1999). On the basis of field data, the system is considered in isothermal (50°C) and isobaric (128.5 bar) conditions. The effects of the evolving reservoir gas composition are taken into account before, during, and after CO₂ injection. Fully water-saturated conditions were assumed for the cement sheath and caprock domains. The gas phase does not flow by advection from the reservoir into the interacting domains so that molecular diffusion in the aqueous phase is the most important process controlling the mass transport occurring in the system under study.     

CO2 Plume Evolution and Pressure Buildup of Large-scale CO2 Injection into Saline Aquifers in Sanzhao Depression, Songliao Basin, China

Large-scale injection of carbon dioxide (CO₂) into saline aquifers in sedimentary basins is a promising approach to mitigate global climate change. Songliao Basin, a large continental clastic sedimentary basin in northeastern China, is one of the great potential candidate sites for future CO₂ storage in China. In this paper, a three-dimensional CO₂ storage model was built to evaluate the CO₂ plume evolution and pressure buildup of large-scale CO₂ injection into the saline aquifers in the Sanzhao Depression of the Songliao Basin. CO₂ was injected into the aquifers through five wells, each with an annual injection rate of 3 Mt over 50?years. The results show that the clastic Yaojia formation at the depth between 900 and 1,600 m with thickness of 150 m might be the favorable layer to store a considerable amount of CO₂, and the overlying Nenjiang formation could ensure long-term CO₂ containment. The relative low permeability of the upper part of the Yaojia formation seems to play a role of a secondary seal on carbon storage. In current injection scenario, CO₂ plume migrates into the formations in the southeast of the depression, which could have potential risk of polluting the freshwater. Therefore, the injection site should stay far away from the southeast of the depression. Moreover, it is very crucial to investigate the permeability distribution of the Yaojia formation because it significantly dominates the CO₂ plume migration. After only 6?months of injection, the pressure buildup at each injection site is affected by pressure interference from neighboring sites. The maximum pressure buildup in the formations is 7.8?MPa after 50?years of injection, and it can even reach 10.5?MPa when the injection layers are with lower permeability. The maximum pressure buildup at the bottom of the Nenjiang formation is 6.7?MPa. The gradient of maximum limited formation pressure is about 18?MPa/km, which might cause fractures to open in the formations of the Sanzhao Depression. Continuous injection of CO₂ for 50?years may not cause damage to the caprock even when the lower permeability occurred in the upper part of the Yaojia formation. The safety of CO₂ storage will be enhanced if the upper part of the storage formation has lower permeability than the lower part.     

Analytical Solution Incorporating History-Dependent Processes for Quick Assessment of Capillary Trapping During CO<Subscript>2</Subscript> Geological Storage

Carbon storage in saline formations is considered as a promising option to ensure the necessary decrease of CO₂ anthropogenic emissions. Its industrial development in those formations is above all conditioned by its safety demonstration. Assessing the evolution of trapped and mobile CO₂ across time is essential in the perspective of reducing leakage risks. In this work, we focus on residual trapping phenomenon occurring during the wetting of the injected CO₂ plume. History dependent effects are of first importance when dealing with capillary trapping. We then apply the classical fractional flow theory (Buckley–Leverett type model) and include trapping and hysteresis models; we derive an analytical solution for the temporal evolution of saturation profile and of CO₂ trapped quantity when injecting water after the gas injection (“artificial imbibition”). The comparison to numerical simulations for different configurations shows satisfactory match and justifies, in the case of industrial CO₂ storage, the assumptions of incompressible flow with no consideration of capillary pressure. The obtained analytical solution allows the quick assessment of both the quantity and the location of mobile gas left during imbibition.     

Dissolution of CO2 From Leaking Fractures in Saline Formations

Leakage of CO₂ through fractures in saline formations will increase the CO₂??brine interface and promote CO₂ dissolution. We use a 2D, finite difference MATLAB model to simulate dissolution rates from a vertical fracture, with CO₂ flowing through it, in a secondary storage formation. The instigation of convection currents increases dissolution rates leading to higher dissolution in higher Rayleigh number systems. Comparison of our results with fracture flow rates shows that for typical fracture apertures dissolution from a fracture is small relative to the amount of CO₂ flowing through the fracture. Temporal and spatial variations in fracture permeability may reduce fracture flow rates and increase the relative amount of CO₂ dissolved from the fracture compared to the CO₂ flowing through the fracture. Further work on CO₂ dissolution in relation to fracture heterogeneity, flow of CO₂ within fractures and the interaction of multiple fractures will improve our ability to predict CO₂ dissolution rates for site characterisation.     

Sensitivity Study of Simulation Parameters Controlling CO<Subscript>2</Subscript> Trapping Mechanisms in Saline Formations

The primary purpose of this study is to understand quantitative characteristics of mobile, residual, and dissolved CO₂ trapping mechanisms within ranges of systematic variations in different geologic and hydrologic parameters. For this purpose, we conducted an extensive suite of numerical simulations to evaluate the sensitivities included in these parameters. We generated two-dimensional numerical models representing subsurface porous media with various permutations of vertical and horizontal permeability (k _v and k _h), porosity (f{\phi}), maximum residual CO₂ saturation (S_gr^max{S_{\rm gr}^{\max}}), and brine density (ρ _br). Simulation results indicate that residual CO₂ trapping increases proportionally to k_v, k_h, S_gr^max{k_{\rm v}, k_{\rm h}, S_{\rm gr}^{\max}} and ρ _br but is inversely proportional to f.{\phi.} In addition, the amount of dissolution-trapped CO₂ increases with k _v and k _h, but does not vary with f{\phi } , and decreases with S_gr^max{S_{\rm gr}^{\max}} and ρ _br. Additionally, the distance of buoyancy-driven CO₂ migration increases proportionally to k _v and ρ _br only and is inversely proportional to k_h, f{k_{\rm h}, \phi } , and S_gr^max{S_{\rm gr}^{\max}} . These complex behaviors occur because the chosen sensitivity parameters perturb the distances of vertical and horizontal CO₂ plume migration, pore volume size, and fraction of trapped CO₂ in both pores and formation fluids. Finally, in an effort to characterize complex relationships among residual CO₂ trapping and buoyancy-driven CO₂ migration, we quantified three characteristic zones. Zone I, expressing the variations of S_gr^max{S_{\rm gr}^{\max}} and k _h, represents the optimized conditions for geologic CO₂ sequestration. Zone II, showing the variation of f{\phi} , would be preferred for secure CO₂ sequestration since CO₂ has less potential to escape from the target formation. In zone III, both residual CO₂ trapping and buoyancy-driven migration distance increase with k _v and ρ _br.     

Matched Boundary Extrapolation Solutions for CO<Subscript>2</Subscript> Well-Injection into a Saline Aquifer

The injection of supercritical CO₂ through wells into deep brine reservoirs is a topic of interest for geologic carbon sequestration. The injected CO₂ is predominantly immiscible with the brine and its low density relative to brine leads to strong buoyancy effects. The displacement of brine by CO₂ in general is a multidimensional, complex nonlinear problem that requires numerical methods to solve. The approximations of vertical equilibrium and complete gravity segregation (sharp interface) have been introduced to reduce the complexity and dimensionality of the problem. Furthermore, for the radial displacement process considered here, the problem can be formulated in terms of a similarity variable that reduces spatial and temporal dependencies to a single variable. However, the resulting ordinary differential equation is still nonlinear and exact solutions are not available. The existing analytical solutions are approximations limited to certain parameter ranges that become inaccurate over a large portion of the parameter space. Here, I use a matched boundary extrapolation method to provide much greater accuracy for analytical/semi-analytical approximations over the full parameter range.