Pressure Buildup During Supercritical Carbon Dioxide Injection From a Partially Penetrating Borehole into Gas Reservoirs

Injecting CO₂ into a subsurface formation causes a buildup of pressure in the vicinity of the injection well. While a large injection rate can reduce the cost associated with injection, an indefinitely large injection rate can result in excessive formation damage. To obtain an optimal injection rate without exceeding the safe pressure limits, one will like to have some knowledge of the transient pressure buildup characteristics resulting from a particular injection rate. While elaborate numerical simulations can provide reliable pressure buildup predictions, they require extensive knowledge about the formation, which is normally not available at the start of an injection process. To alleviate this problem, using some simplifying assumptions, we have developed a solution to predict the transient buildup of pressure resulting from injection of supercritical carbon dioxide from a partially penetrating well into a gas reservoir. The solution in space and time is first obtained in the Fourier–Laplace transform space, and then inverted back into real space (in cylindrical coordinates) and time. We use the solution to study pressure transient characteristics for different formation permeabilities and anisotropy ratios. Results obtained using the solution compared well with those from numerical simulations.     

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.     

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.     

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.     

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.     

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.     

Effect of Vertical Heterogeneity on Long-Term Migration of CO<Subscript>2</Subscript> in Saline Formations

For deep injection of CO₂ in thick saline formations, the movements of both the free gas phase and dissolved CO₂ are sensitive to variations in vertical permeability. A simple model for vertical heterogeneity was studied, consisting of a random distribution of horizontal impermeable barriers with a given overall volume fraction and distribution of lengths. Analytical results were obtained for the distribution of values for the permeability, and compared to numerical simulations of deep CO₂ injection and convection in heterogeneous formations, using multiple realizations for the permeability distribution. It is shown that for a formation of thickness H, the breakthrough times in two dimensions for deep injection scale as H ² for moderate injection rates. In comparison to heterogeneous shale distributions, a homogeneous medium with equivalent effective vertical permeability has a longer breakthrough time for deep injection, and a longer onset time for convection.     

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.     

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.     

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.     

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.     

Modeling of CO<Subscript>2</Subscript> Leakage up Through an Abandoned Well from Deep Saline Aquifer to Shallow Fresh Groundwaters

This article presents a numerical modeling application using the code TOUGHREACT of a leakage scenario occurring during a CO₂ geological storage performed in the Jurassic Dogger formation in the Paris Basin. This geological formation has been intensively used for geothermal purposes and is now under consideration as a site for the French national program of reducing greenhouse gas emissions and CO₂ geological storage. Albian sandstone, situated above the Dogger limestone is a major strategic potable water aquifer; the impacts of leaking CO₂ due to potential integrity failure have, therefore, to be investigated. The present case–study illustrates both the capacity and the limitations of numerical tools to address such a critical issue. The physical and chemical processes simulated in this study have been restricted to: (i) supercritical CO₂ injection and storage within the Dogger reservoir aquifer, (ii) CO₂ upwards migration through the leakage zone represented as a 1D vertical porous medium to simulate the cement–rock formation interface in the abandoned well, and (iii) impacts on the Albian aquifer water quality in terms of chemical composition and the mineral phases representative of the porous rock by estimating fluid–rock interactions in both aquifers. Because of CPU time and memory constraints, approximation and simplification regarding the geometry of the geological structure, the mineralogical assemblages and the injection period (up to 5 years) have been applied to the system, resulting in limited analysis of the estimated impacts. The CO₂ migration rate and the quantity of CO₂ arriving as free gas and dissolving, firstly in the storage water and secondly in the water of the overlying aquifer, are calculated. CO₂ dissolution into the Dogger aquifer induces a pH drop from about 7.3 to 4.9 limited by calcite dissolution buffering. Glauconite present in the Albian aquifer also dissolves, causing an increase of the silicon and aluminum in solution and triggering the precipitation of kaolinite and quartz around the intrusion point. A sensitivity analysis of the leakage rate according to the location of the leaky well and the variability of the petro-physical properties of the reservoir, the leaky well zone and the Albian aquifers is also provided.     

Analytical Model for CO2 Injection into Brine Aquifers-Containing Residual CH4

During CO₂ injection into brine aquifers-containing residual and/or dissolved CH₄, three distinct regions develop: (1) a single-phase, dry-out region around the well-bore filled with pure supercritical CO₂; (2) a two-phase, two-component system containing CO₂ and brine; and (3) a two-phase, two-component system containing CH₄, and brine. This article extends an existing analytical solution, for pressure buildup during CO₂ injection into brine aquifers, by incorporating dissolved and/or residual CH₄. In this way, the solution additionally accounts for partial miscibility of the CO₂?CCH₄?Cbrine system and the relative permeability hysteresis associated with historic imbibition of brine and current drainage due to CO₂ injection and CH₄ bank development. Comparison of the analytical solution results with commercial simulator, CMG-GEM, shows excellent agreement among a range of different scenarios. The presence of residual CH₄ in a brine aquifer summons two competing phenomena, (1) reduction in relative permeability (phase interference), which increases pressure buildup by reducing total mobility, and (2) increase in bulk compressibility which decreases pressure buildup of the system. If initial CH₄ is dissolved (no free CH₄), these effects are not as important as they are in the residual gas scenario. Relative permeability hysteresis increased the CH₄ bank length (compared to non-hysteretic relative permeability), which led to further reduction in pressure buildup. The nature of relative permeability functions controls whether residual CH₄ is beneficial or disadvantageous to CO₂ storage capacity and injectivity in a candid brine aquifer.     

Investigation of hydrodynamic instability of CO<Subscript>2</Subscript> injection into an aquifer

The two-dimensional problem of supercritical carbon dioxide injection into an aquifer is solved. Shocks and rarefaction waves propagating in a sequence from an injection well into the formation are described within the framework of a complete nonisothermal model of flows in a porous medium. In the approximation of isothermal immiscible water and carbon dioxide flow the hydrodynamic stability of the leading displacement front is investigated for various reservoir pressures and temperatures. The parameters of unstable fronts are determined using a sufficient instability condition formulated in analytic form. The approximate analytic results are supported by the direct numerical simulation of CO₂ injection using the complete model in which thermal effects and phase transitions are taken into account.     

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.     

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.     

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.     

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.     

Quantification of Density-Driven Natural Convection for Dissolution Mechanism in CO<Subscript>2</Subscript> Sequestration

Dissolution of CO₂ into brine causes the density of the mixture to increase. The density gradient induces natural convection in the liquid phase, which is a favorable process of practical interest for CO₂ storage. Correct estimation of the dissolution rate is important because the time scale for dissolution corresponds to the time scale over which free phase CO₂ has a chance to leak out. However, for this estimation, the challenging simulation on the basis of convection–diffusion equation must be done. In this study, pseudo-diffusion coefficient is introduced which accounts for the rate of mass transferring by both convection and diffusion mechanisms. Experimental tests in fluid continuum and porous media were performed to measure the real rate of dissolution of CO₂ into water during the time. The pseudo diffusion coefficient of CO₂ into water was evaluated by the theory of pressure decay and this coefficient is used as a key parameter to quantify the natural convection and its effect on mass transfer of CO₂. For each experiment, fraction of ultimate dissolution is calculated from measured pressure data and the results are compared with predicted values from analytical solution. Measured CO₂ mass transfer rate from experiments are in reasonable agreement with values calculated from diffusion equation performed on the basis of pseudo-diffusion coefficient. It is suggested that solving diffusion equation with pseudo diffusion coefficient herein could be used as a simple and rapid tool to calculate the rate of mass transfer of CO₂ in CCS projects.     

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.