STREAMLINE-BASED MATHEMATICAL MODEL FOR CO2 MISCIBLE FLOODING

According to the research theory of improved black oil simulator, a practical mathematical model for C02 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance fast, streamline method is employed to solve this model as an alternative to traditional finite difference methods. Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method (BEM), an explicit total variation diminishing (TVD) method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less information need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an effective tool for early stage screening andmiscible oil field.reservoir dynamic management of the CO2 miscible oil field.     

Laboratory Experiments on Environmental Friendly Means to Improve Coalbed Methane Production by Carbon Dioxide/Flue Gas Injection

Scaled in situ laboratory core flooding experiments with CO₂, N₂ and flue gas were carried out on coal in an experimental high P,T device. These experiments will be able to give an insight into the design of the injection system, management, control of the operations and the efficiency of an ECBM project. Although the experience gained by the oil industry represents a valuable starting point, several problems are still to be studied and solved before CO₂ improved deep coalbed methane production may be operationally feasible. These are all related to the heterogeneous nature of the pore structure of coal, and in particular to the presence of fractures. More specifically, a number of questions need to be addressed, e.g. what are the conditions under which the fluid in the micro pores of the coal is displaced by the CO₂ in the presence of competitive adsorption; what is the role of compositional heterogeneity and fracture anisotropy of coal for the injection design and the efficiency of the sequestration in relation to the swelling and shrinkage characteristics of coal; how does the mobile and the immobile water in the coal affect the exchange process. These questions can be answered by means of downscaled laboratory experiments that are capable of accurately describing the coupled process of multiphase flow, competitive adsorption and geo-mechanics. The laboratory conditions have been simulated to match pressure and temperature at depths of 800 to 1,000 m. Under those conditions the injected CO₂ remains supercritical. Upto now, the results show that dewatering will be an essential step for successful ECBM combined with a CO₂ sequestration process.     

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.     

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.     

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.     

Adsorption and strain: The CO2-induced swelling of coal

Enhanced coal bed methane recovery (ECBM) consists in injecting carbon dioxide in coal bed methane reservoirs in order to facilitate the recovery of the methane. The injected carbon dioxide gets adsorbed at the surface of the coal pores, which causes the coal to swell. This swelling in confined conditions leads to a closure of the coal reservoir cleat system, which hinders further injection. In this work we provide a comprehensive framework to calculate the macroscopic strains induced by adsorption in a porous medium from the molecular level. Using a thermodynamic approach we extend the realm of poromechanics to surface energy and surface stress. We then focus on how the surface stress is modified by adsorption and on how to estimate adsorption behavior with molecular simulations. The developed framework is here applied to the specific case of the swelling of CO₂-injected coal, although it is relevant to any problem in which adsorption in a porous medium causes strains.     

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.     

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.     

Coreflooding Studies to Investigate the Potential of Carbonated Water Injection as an Injection Strategy for Improved Oil Recovery and CO<Subscript>2</Subscript> Storage

Carbonated water injection (CWI) is a CO₂-augmented water injection strategy that leads to increased oil recovery with added advantage of safe storage of CO₂ in oil reservoirs. In CWI, CO₂ is used efficiently (compared to conventional CO₂ injection) and hence it is particularly attractive for reservoirs with limited access to large quantities of CO₂, e.g. offshore reservoirs or reservoirs far from large sources of CO₂. We present the results of a series of CWI coreflood experiments using water-wet and mixed-wet Clashach sandstone cores and a reservoir core with light oil (n-decane), refined viscous oil and a stock-tank crude oil. The experiments were carried out to assess the performance of CWI and to quantify the level of additional oil recovery and CO₂ storage under various experimental conditions. We show that the ultimate oil recovery by CWI is higher than the conventional water flooding in both secondary and tertiary recovery methods. Oil swelling as a result of CO₂ diffusion into the oil and the subsequent oil viscosity reduction and coalescence of the isolated oil ganglia are amongst the main mechanisms of oil recovery by CWI that were observed through the visualisation experiments in high-pressure glass micromodels. There was also evidence of a change in the rock wettability that could also influence the oil recovery. The coreflood test results also reveal that the CWI performance is influenced by oil viscosity, core wettability and the brine salinity. Higher oil recovery was obtained with the mixed-wet core than the water-wet core, with light oil than with the viscous oil and low salinity carbonated brine than high-salinity carbonated brine. At the end of the flooding period, an encouraging amount of the injected CO₂ was stored in the brine and the remaining oil in the form of stable dissolved CO₂. The experimental results clearly demonstrate the potential of CWI for improving oil recovery as compared with the conventional water flooding (secondary recovery) or as a water-based EOR (enhanced oil recovery) method for watered out reservoirs.     

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.     

Competitive Methane Desorption by Supercritical CO2 Injection in Coal

A large diameter (∼70 mm) dry coal sample was used to study the competitive displacement of CH₄ by injection of supercritical CO₂, and CO₂–CH₄ counter-diffusion in coal matrix. During the test, a staged loading procedure, which allows the calibration of the key reservoir modelling parameters in a sequential and progressive manner, was employed. The core-flooding test was history matched using an Enhanced Coalbed Methane (ECBM) simulator, in which Fick’s Law for mixed gas diffusion and the extended Langmuir equations are implemented. The system pressure rise during the two loading stages and the CO₂ breakthrough time in the final production stage were matched by using the pair of constant sorption times (9 and 3.2 days) for CH₄ and CO₂, respectively. The corresponding diffusion coefficients for CH₄ and CO₂ were estimated to be 1.6 ×  10⁻¹² and 4.6 ×  10⁻¹² m²/s, respectively. Comparison was made with published gas diffusion coefficients for dry ground samples (ranging from < 0.063 to ∼3 mm) of the same coal at relatively low pressures (< 4 MPa). The CO₂/CH₄ gas diffusion coefficient ratio was well within the reported range (2–3), whereas the CH₄ diffusion coefficient obtained from history matching of the core-flooding test is approximately 15 times smaller than that arrived by curve-fitting the measured sorption uptake rate using a unipore diffusion model. The calibrated model prediction of the effluent gas composition was in good agreement with the test data for CO₂ mole fraction of up to 20%.     

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.     

Comments on “Abrupt-Interface Solution for Carbon dioxide Injection into Porous Media” by M. Dentz and D. Tartakovsky

Several numerical simulations were conducted to compare the results with analytical solutions given by “Abrupt-Interface Solution for Carbon Dioxide Injection into Porous Media” by M. Dentz and D. Tartakovsky, “Injection and Storage of CO₂ in Deep Saline Aquifers: Analytical Solution for CO₂ Plume Evolution During Injection” by J. M. Nordbotten et al.     

Front Tracking Using a Hybrid Analytical Finite Element Approach for Two-Phase Flow Applied to Supercritical CO<Subscript>2</Subscript> Replacing Brine in a Heterogeneous Reservoir and Caprock

Predicting fluid replacement by two-phase flow in heterogeneous porous media is of importance for issues such as supercritical CO₂ sequestration, the integrity of caprocks and the operation of oil water/brine systems. When considering coupled process modelling, the location of the interface is of importance as most of the significant interaction between processes will be happening there. Modelling two-phase flow using grid based techniques presents a problem as the fluid–fluid interface location is approximated across the scale of the discretisation. Adaptive grid methods allow the discretisation to follow the interface through the model, but are computationally expensive and make coupling to other processes (thermal, mechanical and chemical) complicated due to the constant alteration in grid size and effects thereof. Interface tracking methods have been developed that apply sophisticated reconstruction algorithms based on either the ratio of volumes of a fluid in an element (Volume of Fluid Methods) or the advective velocity of the interface throughout the modelling regime (Level set method). In this article, we present an “Analytical Front Tracking” method where a generic analytical solution for two-phase flow is used to “add information” to a finite element model. The location of the front within individual geometrical elements is predicted using the saturation values in the elements and the velocity field of the element. This removes the necessity for grid adaptation, and reduces the need for assumptions as to the shape of the interface as this is predicted by the analytical solution. The method is verified against a standard benchmark solution and then applied to the case of CO₂ pooling and forcing its way into a heterogeneous caprock, replacing hot brine and eventually breaking through. Finally the method is applied to simulate supercritical CO₂ injected into a brine saturated heterogeneous reservoir rock leading to significant viscous fingering and developement of preferential flow paths. The results are compared with to a finite volume simulation.     

Enhanced CO<Subscript>2</Subscript> Storage and Sequestration in Deep Saline Aquifers by Nanoparticles: Commingled Disposal of Depleted Uranium and CO<Subscript>2</Subscript>

Geological storage of anthropogenic CO₂ emissions in deep saline aquifers has recently received tremendous attention in the scientific literature. Injected buoyant CO₂ accumulates at the top part of the aquifer under a sealing cap rock. Potential buoyant movement of CO₂ has caused some concern that the high-pressure CO₂ could breach the seal rock. However, CO₂ will diffuse into the brine underneath and generate a slightly denser fluid that may induce instability and convective mixing. Onset times of instability and convective mixing performance depend on the physical properties of the rock and fluids, such as permeability and density contrast. We present the novel idea of adding nanoparticles (NPs) to injected CO₂ to increase density contrast between the CO₂-rich brine and the underlying resident brine and, consequently, decrease onset time of instability and increase convective mixing. The analyses show that 0.001 volume fraction of NPs added to the CO₂ stream shortens onset time of mixing by approximately 80% and increases convective mixing by 50%. If it thus originally takes 5 years for the overlying CO₂ to start convective mixing, by adding NPs, onset time of mixing reduces to 1 year, and after initiation of convective mixing, mixing improves by 50%. A reduction of the CO₂ leakage risk ensues. In addition to other metallic NPs, use of processed depleted uranium oxide (DU) as the NPs is also proposed. DU-NPs are potentially stable and might be safely commingled with CO₂ to store in saline aquifers.     

Phase Behaviour,Fluid Properties and Recovery Efficiency of Immiscible and Miscible Condensate Displacements by SCCO<Subscript>2</Subscript> Injection: Experimental Investigation

This paper presents a quantitative investigation of the interfacial tension dependent relative permeability (IFT-DRP) and displacement efficiency of supercritical CO₂ injection into gas-condensate reservoirs. A high-pressure high-temperature experimental laboratory was established to simulate reservoir conditions and to perform relative permeability measurements on sandstone cores at a constant reservoir temperature of 95°C and displacement velocity of 10 cm/h. This investigation covers immiscible displacements (1100 and 2100 psi), near-miscible displacement (3000 psi) and miscible displacements (4500 and 5900 psi). The coreflooding results demonstrated that displacement pressure is a key factor governing the attainment of optimum sweep efficiency. The ultimate condensate recovery increased by almost threefold when CO₂ was injected at near-miscible conditions (i.e., 23.40% ultimate recovery at 1100 psi compared to 69.70% at 3000 psi). Miscible flooding was found to give the optimum condensate recovery (9% extra ultimate recovery compared to near-miscible injection). Besides improving the ultimate recovery, miscible floods provided better mobility ratios and delayed gas breakthrough (0.62 PV BT at 5900 psi compared to 0.21 PV BT at 1100 psi). In addition to the elimination of IFT forces in miscible displacements, favourable ratios of fluid properties and phase behaviour relationships between the SCCO₂ and condensate were believed to be the driving force for the improved recovery as they provided a stabilising effect on the displacement front and stimulated swelling of the condensate volume. This paper incorporates the theoretical aspects of phase behaviour and fluid properties that largely affect the microscopic displacement efficiency and serves as a practical guideline for operators to aid their project designs and enhance their recovery capabilities.     

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.     

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.