Next Article in Journal
Differentiation Study of the Damage Characteristics of Rock Cultural Heritage Sites Due to the Sulfate Weathering Process
Previous Article in Journal
Edible Oils from Selected Unconventional Sources—A Comprehensive Review of Fatty Acid Composition and Phytochemicals Content
Previous Article in Special Issue
Capture of CO2 Using Mixed Amines and Solvent Regeneration in a Lab-Scale Continuous Bubble-Column Scrubber
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Potential Benefits of Horizontal Wells for CO2 Injection to Enhance Storage Security and Reduce Leakage Risks

by
Marcos Vitor Barbosa Machado
1,*,
Mojdeh Delshad
2,* and
Kamy Sepehrnoori
2
1
PETROBRAS, Rio de Janeiro 20231-030, RJ, Brazil
2
Hildebrand Department of Petroleum and Geosystems Engineering, The University of Texas at Austin, Austin, TX 78712, USA
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(23), 12830; https://doi.org/10.3390/app132312830
Submission received: 12 October 2023 / Revised: 17 November 2023 / Accepted: 19 November 2023 / Published: 29 November 2023

Abstract

:

Featured Application

This work suggests that horizontal injectors are a safe way to keep CO2 trapped in the storage site and reduce/prevent any unexpected event associated with leakage through the caprock since the horizontal well minimizes or avoids CO2 contact with the caprock.

Abstract

This study used numerical simulations of CO2 storage to identify the benefits of horizontal wells for geological carbon storage, such as enhancing CO2 trapped in porous media due to relative permeability and capillary hysteresis. Two injection schemes were tested: one using a vertical injector and the other employing a horizontal well. The results revealed two main findings. Firstly, the horizontal injection well effectively prevented or minimized CO2 penetration into the caprock across various sensitivity scenarios and over a thousand years of CO2 redistribution. Secondly, horizontal wells provided a safe approach to trapping CO2, increasing its entrapment as a residual phase by up to 19% within the storage site. This, in turn, reduced or prevented any unexpected events associated with CO2 leakage through the caprock. Additionally, the paper proposes a practical method for designing the optimal length of a horizontal well. This method considers a combination of two parameters: the additional CO2 that can be trapped using a horizontal well and the gravity number. In the case of the reservoir model of this study, a horizontal branch with a length of 2000 m was found to be the most effective design in enhancing CO2 entrapment and reducing CO2 buoyancy.

1. Introduction

Scenarios considered by international entities such as the International Agency (IEA) point out the complexity of actions needed to achieve goals established for reducing emissions of greenhouse gases in the coming decades, simultaneously with the desired transformation in energy generation sources. The report “World Energy Outlook 2022” [1] emphasizes that multiple technologies and energy sources will play an essential joint role in providing energy resources for the planet in a more sustainable scenario. Among the various options for reducing CO2 emissions, Carbon Capture and Storage (CCS) presents itself as a technology with significant potential to reduce CO2 emissions in a (probable) scenario of continued use of fossil fuels in the coming decades.
A CCS project involves CO2 capture from high-emission industries and injecting it into geological formations such as aquifers and depleted hydrocarbon reservoirs. One of the most critical barriers to long-term and large-volume CO2 storage in geological formations is the proof of safe and reliable storage. For that, CO2 can take advantage of different trapping mechanisms in porous media [2,3,4,5]; for example, free CO2 migration is controlled by the structural and stratigraphic trapping exerted by the caprock during the short term encompassing the injection time, known as a primary mechanism. During this stage, the caprock plays an essential role in the security of the CCS operation due to CO2 buoyancy which can move it up to the surface or the seafloor in the case of offshore storage sites. In the Frio CO2 field demonstration project conducted in the U.S. [6], CO2 was injected in a deeper zone below several well-known shale seals. However, some authors [7] are currently reviewing the need for a caprock to control the plume rise because of the following reasons: (i) there is no prescriptive regulation concerning geologic seals; and (ii) the concept of the composite confining system is introduced, as a set of discontinuous barriers which create long and tortuous paths that attenuate mobile CO2 rise.
During the mid- and long term, part of mobile CO2 will be dissolved in water (solubility trapping), as time passes, especially in low-salinity brine under high-pressure and low-temperature conditions [8,9]. Saturation changes caused by the rising plume can lead to more CO2 being trapped as a residual phase due to relative permeability and capillary hysteresis [10,11,12]. Beyond that, in few cases, some CO2 can be trapped as minerals because of the resulting pH of the brine and the mineralogy of the rock [13,14]. Therefore, those additional mechanisms increase storage security, when the buoyant CO2 is immobile in the pore space or no longer exists as a free phase (generally as super-critical CO2).
Modeling CO2 retention in porous media depends on the rock properties, such as porosity and permeability, and intrinsic parameters related to each trapping mechanism mentioned beforehand [15]. In this study, the entrapment of CO2 in porous media will be rigorously modeled to allow a fair evaluation between two well geometries to inject CO2: vertical and horizontal wells. The hypothesis to be investigated is related to CO2 spreading around the well. With vertical wells, CO2 is concentrated nearby with a limited spread due to the dominance of gravitational forces. However, with horizontal wells, CO2 is broadly spread in porous media, developing a larger swept area, and it is thus likely to become trapped in that area due to hysteresis as a residual phase. Besides that, there is a possibility of contacting poor CO2-saturated brine, enhancing the solubility and ionic entrapment. This analysis will be addressed during this comparison, as well as the better equilibrium between gravity and viscous forces provided by horizontal wells due to the main contribution of vertical permeability for the flow, reducing the buoyancy CO2 flux.
The idea that the mass of residually trapped CO2 is enhanced, when the plume migrates in the lateral direction was observed numerically by Han et al. [16], who investigated different levels of vertical and horizontal permeabilities in a model with a source point of injection. In the conclusion section of this work, the authors postulated that horizontal wells can be a likely means to increase the lateral extension of CO2 plumes, taking advantage of the residual trapping mechanism. However, they did not simulate cases with horizontal wells, as is in the goals of our current study. In addition to that, our paper aims to expand the study by Okwen et al. [17]. They compared vertical and horizontal wells to store CO2 in saline aquifers but did not consider the effect of hysteresis as a trapping mechanism, because their investigation was focused only on the injection time (50 years), where the system only underwent the drainage path. In other words, our work intends to evaluate the long-term CO2 storage, modeling all of the trapping mechanisms that can occur.
Therefore, the innovative aspects of our work are summarized below:
  • Providing technical arguments in favor of using horizontal wells in low- and mid-permeability targets (average less than 100–200 mD), considering the greater mass of CO2 trapped as a residual fluid and the better equilibrium between gravity and viscous forces to control the CO2 buoyancy, evaluating long-term storage (thousands of years of redistribution);
  • Proposing a practical way to design the horizontal well length to maximize the CO2 entrapment and make the flow less gravity-driven;
  • Evaluating the impact of both strategies (vertical and horizontal injectors) on the risk of CO2 penetration into the caprock, considering adsorption, diffusion, and the Darcy reactive flow.
The horizontal wells are applied in actual operating CCS projects worldwide. Sleipner in the North Sea and In Salah in the Sahara Desert are examples of storage sites with horizontal injectors [18]. Therefore, the broader impact of this research is intrinsically related to better control of the CO2 plume redistribution, which has the following advantages:
  • Preventing an unexpected and undesired interaction with the caprock, mainly related to the risk of CO2 leakage due to fault/fracture activation or the creation of new fractures when the injection exceeds the minimum rock strength [19,20]. In addition, related to the geochemical reactions among the CO2−brine−mineral system, which can affect the integrity of caprock and its sealing capacity with alteration in its petrophysical properties induced by mineral changes [21];
  • Better control and monitoring of the CO2 plume propagation to the CCS project evaluation with data from the injection or observation/monitoring wells, besides geophysical data, such as time-lapse seismic, vertical-seismic profiles (VSPs), and micro-seismic [5];
  • Providing insights into the CO2 storage in depleted oil/gas reservoirs where it is possible to use existing horizontal and multilateral wells.
The following conditions are outside the scope of this paper:
  • Consideration of impurities and free water content in CO2 stream injected;
  • Modeling the dry-out effect due to water vaporization with CO2 injection or another injectivity issue, which can be found in Machado et al. [22];
  • CO2 leakage through wells with poor cement jobs, as pointed out by Gholami et al. [23], can be the most important reason behind the migration and leakage;
  • Drilling design planning and stability concerns for horizontal well construction;
  • Investigation of solutions with horizontal branches to inject CO2 and produce brine to relieve pressure buildup [24,25];
  • Economic evaluation of vertical and horizontal wells. This specific evaluation would depend on factors such as the environment (onshore or offshore), well depth and length, and location of the operation.

2. Petrophysical Modeling for Sandstone and Shale

This study considered two sandstone geological models: homogeneous and heterogeneous saline aquifers (Figure 1). The caprock was modeled as a water-saturated shale, assuming homogeneous properties in both cases. The homogeneous model (Figure 1A) was 1000 m in width and 1200 m in thickness, of which 1000 m was the shale caprock. The gridblock size was 10 × 10 m2 horizontally and 1 m thick.
Figure 1B shows the heterogeneous model with 7000 m in width and 2500 m in thickness (including the thickness of the shale caprock of about 1500 m). This latter model represented part of an actual sandstone aquifer, and it was built with seismic and well data from three wells drilled in that area, located at the Campos basin, Rio de Janeiro state in Brazil. This sector model had a pore volume of about 160 million m3; the gridblock size was 100 × 100 m2 horizontally and 5 m thick. A grid refinement was performed around the injector to reduce the gridblock size from 100 m to 25 m, using the workflow proposed by Machado et al. [26].
Table 1 summarizes the key petrophysical properties assigned to both models.
The relative permeability curves for CO2 and brine in sandstone and shale were obtained by Bennion and Bachu [27]. The CO2—brine capillary pressure curves were obtained from a J-function fitted to the data from Abdoulghafour et al. [28] for the sandstone and from Bennion and Bachu [29] for the shale. As the mean ratio k/φ was the same in both models, the capillary pressure curve resulting from the J-function is also the same. Figure 2 and Figure 3 show the final drainage curves used in the simulation for the aquifer and the caprock, respectively. The imbibition curves will be generated according to the hysteresis model to be detailed in the next section.
The capillary pressure in the shale caprock exhibited a threshold of 318 kPa, as depicted in Figure 3. This threshold represented the minimum pressure required to exceed the capillary forces and enable free super-critical CO2 migration into the caprock.

3. Modeling CO2 Entrapment for Sandstone and Shale Formations

The sandstone saline aquifer and its shale caprock in the aforementioned models were numerically simulated using CMG-GEM [30]. CMG-GEM is known for its versatility and robustness in predicting and analyzing various thermodynamic properties and fluid behavior for CCS projects [2,3,11,22,26]. It is based on the discretization of component conservation and energy balance equations in space and time using finite volume and finite difference methods [31,32].
One of the key strengths of the CMG-GEM model lies in its accurate prediction of CO2 phase behavior for wide ranges of pressure and temperature. It can effectively capture the transition of CO2 from gaseous to liquid states, as well as its behavior in supercritical conditions, which is crucial for understanding and studying CO2 storage. Overall, the CMG-GEM simulation model enables accurate prediction and analysis of CO2 behavior and the following trapping mechanisms:
  • CO2 solubility in brine according to the method by Li and Nghiem [33] based on Henry’s law. This model is based on Henry’s constant calculation according to Equation (1) as a function of pressure and temperature. The effect of salt on the gas solubility in the aqueous phase is modeled by the salting-out coefficient [34].
    l n H i = l n H i * + v i ¯ R T p p *
    where:
    • H i is Henry’s constant at current pressure (p) and temperature (T);
    • H i * is Henry’s constant at reference pressure (p*) and temperature (T);
    • v i ¯ is the partial molar volume at infinite dilution;
    • R is universal gas constant;
    • I is species dissolved in water (CO2(aq) in this work).
The   H i * constant was computed considering the initial conditions of the models in Table 2.
  • Solubility trapping in brine can be enhanced by diffusion. To model this effect, the diffusion coefficient (D) for super-critical CO2 in brine is applied to compute the effective CO2 diffusion (Deff) considering tortuosity τ [35]:
    D e f f = D τ
    In shales, where diffusion is a relevant mechanism due to lower permeability, the tortuosity value ranges approximately from 40 to 70 [36]. An effective diffusion equal to 2.8 × 10−7 cm2/s [37] was used in the simulations.
  • Chemical trapping by adsorption of CO2 in the gas phase was modeled using a Langmuir model [38], which is a widely accepted isotherm adsorption equation [39,40]. With data for shales [41,42] and sandstone [43], two isotherms were matched to model the adsorption. Table 3 summarizes Langmuir parameters (BCO2 and ωCO2,max) obtained according to Equation (3), that is, the extended Langmuir isotherm for multicomponent adsorption [44,45]:
    ω C O 2 = ω C O 2 ,   m a x B C O 2 y C O 2 ,   g p 1 + p j B j y j ,   g
    where
    • BCO2 is the parameter for Langmuir isotherm relation;
    • BCO2 is the moles of adsorbed CO2 per unit mass of rock;
    • ωCO2,max is the maximum moles of adsorbed CO2 per unit mass of rock;
    • YCO2,g is the molar fraction of adsorbed CO2 in the gas phase.
      Table 3. Parameters of Langmuir isotherm to model the CO2 adsorption.
      Table 3. Parameters of Langmuir isotherm to model the CO2 adsorption.
      SandstoneShale
      B C O 2 0.00017 kPa−10.00017 kPa−1
      ω C O 2 ,   m a x 0.07 gmol/kg0.18 gmol/kg
  • The residual CO2 trapping due to the relative permeability and capillarity hysteresis with the saturation changes was modeled with the maximum gas trapped (Sgt) converted to the Land’s constant (C) [46] in the two-phase Carlson’s model [47], as recommended by Jarrell et al. [48], according to Equation (4):
    S g t = S g   m a x 1 + C S g   m a x
    where Sg max is the maximum gas saturation.
Burnside and Naylor [49] obtained a Sgt distribution for sandstone, shale, and carbonates based on more than 30 published coreflood data for CO2 and brine. The mean value for each lithology is presented in Table 4.
  • Ionic trapping due to acidic water reactions for bicarbonate and carbonate ions generation was modeled using kinetic parameters from the PHREEQC database [50,51].
    OH + H+ = H2O
    CO2 + H2O = H+ + HCO3
    CO32− + H+ = HCO3
    The synthetic water composition is given in Table 5 for the homogeneous case [52] and the heterogeneous case (Personal communication with PETROBRAS. 2023. Rio de Janeiro-RJ, Brazil).
    Table 5. Ionic composition of formation brine.
    Table 5. Ionic composition of formation brine.
    IonsMolality Concentration (gmol/kgw)
    HomogenousHeterogeneous
    pH8.86.8
    Ca2+0.04750.027
    Na+0.32341.272
    Cl0.5511.3
    Mg2+0.00840.028
    K+0.00710.013
    SO42−0.0220.000085
  • The following mineralization reactions with primary minerals were modeled using kinetic parameters from PHREEQC for Transition State Theory (TST)-derived rate laws:
    Quartz [SiO2] = SiO2 (aq)
    Kaolinite [Al2Si2O5(OH)4] + 6.0 H+ = 5.0 H2O + 2 Al3+ + 2 SiO2 (aq)
    Calcite [CaCO3] + H+ = Ca2+ + HCO3
    Illite [K0.6Mg0.25Al2.3Si3.5O10(OH)2 + 11.2 H2O = 3.5 H4SiO4 + 2.3 Al(OH)4 + 0.6 K+ + 0.25 Mg2+ + 1.2 H+
    Albite [NaAlSi3O8] + 4 H+ = 3 SiO2 (aq) + Al3+ + Na+ + 2 H2O
    Anorthite [CaAl2Si2O8] + 8 H2O = 2 H4SiO4 + 2 Al(OH)4 + Ca2+
    Chlorite [Mg5Al2Si3O10(OH)8] + 16.0 H+ = 5.0 Mg2+ + 2.0 Al3+ + 3.0 H4SiO4 + 6.0 H2O
    Dolomite [CaMg(CO3)2] = 2 CO32− + Mg2+ + Ca2+
    Pyrite [FeS2] + 2 H+ = 2 HS + Fe2+
    The initial mineral compositions for the shale caprock [52] and for the sandstone [53] are summarized in Table 6. Illite is the most abundant clay mineral in the shale matrix, occurring in the sandstone matrix. The long-term exposure of clay minerals to Sc-CO2 (super-critical CO2) can generate strong CO2 adsorption [54], justifying the modeling of this phenomenon to increase CO2 storage capacity.
    Table 6. Mineral composition for the sandstone aquifer and the shale caprock.
    Table 6. Mineral composition for the sandstone aquifer and the shale caprock.
    SandstoneSiliceous Shale
    MineralNormalized Volume
    Fraction
    Normalized Volume
    Fraction
    quartz0.6140.277
    kaolinite0.0210.015
    calcite0.0200.073
    illite0.0110.374
    albite0.2100.063
    anorthite0.0760.115
    chlorite0.0480.038
    dolomite0.0000.021
    pyrite0.0000.023
Permeability alteration due to mineral precipitation or dissolution was computed by applying the Kozeny−Carman equation with an exponent value of 3, recommended by Zeidouni et al. [55]:
k k = k n / r f
where
r f = r f n φ n φ k 3 1 φ k 1 φ n 2 ,
where the resistance factor rf is modeled by the Kozeny−Carman equation or the power law relationship, and kn and kk refer to permeability at previous (n) and current (k) timesteps, respectively. The porosity, φ, in Equation (18) is calculated as follows:
φ = 1 + c f p p * φ * j = 1 n N j ρ m , j N j 0 ρ m , j ,
where φ* is the reference porosity for the case without mineral precipitation/dissolution, Nj is the total moles of mineral j per bulk volume at the current time, N j 0   is the total moles of mineral j per bulk volume at the initial time, ρm,j is the mineral molar density, cf is the rock compressibility, and p* is the reference pressure.

4. Potential Changes in the Caprock Integrity with a Vertical Injector

Figure 4 brings up the motivation for this study. It compares the injection of 7000 tonnes of CO2 (3.5% of the initial pore volume) over 30 years through a vertical well perforated at the bottom of the aquifer (Figure 4A) and two cases with horizontal wells (400 m and 800 m in length indicated by Lw; Figure 4B,C), simulated in the homogeneous aquifer model. The horizontal wells were located at the same depth of the shallowest perforation of the vertical well at 1170 m. In this figure, the CO2 in terms of global mole fraction, after 3000 years of redistribution, contacted and penetrated, mainly by diffusion, the caprock, when the injection well was vertical. On the other hand, the longer the horizontal well, the further away the CO2 plume was from the caprock.
From Figure 4, in the vertical case, the upward movement of CO2 was more prominent, leading to contact with the caprock, which had much lower permeability. This resulted in the CO2 plume spreading extensively at the interface between the caprock and the host sandstone formation.
In contrast, in the horizontal well case, the CO2 movement exhibited two distinct peaks, representing the toe and heel of the horizontal well. These peaks indicated areas of higher injectivity compared to other sections of the well, primarily due to the larger open area available for fluid flow. The presence of these peaks in the horizontal well configuration suggests that the CO2 distribution and injectivity were not uniform along the entire length of the well.
A more significant CO2 spreading with a horizontal well allowed more CO2 to be trapped due to hysteresis. Figure 5 in the study illustrates the CO2 inventory after a period of 3000 years of redistribution. The figure utilizes the following terminology to depict the different forms of the CO2 inventory:
  • % free Sc-CO2: free CO2 as a super-critical fluid;
  • % CO2 adsorbed: CO2 adsorbed on the rock by chemical trapping;
  • % CO2 in water: dissolved CO2 in aqueous phase, CO2 (aq);
  • % residual CO2: CO2 trapped as a residual gas due to relative permeability and capillarity hysteresis.
Figure 5 highlights the smaller percentage of free CO2 as a supercritical fluid (1.6% of free Sc-CO2 in purple) and the greater amount as a residual fluid (37% of residual CO2 in yellow) provided when a horizontal injector well was used (Figure 5B). The free Sc-CO2 increased to 13% with a vertical well, and the residual CO2 reduced to 24% (Figure 5A). The CO2 dissolution in the brine (% CO2 in water in blue) was roughly the same in both cases (vertical and horizontal wells), in agreement with the conclusions from Okwen et al. [17] for short-term storage evaluations.
The higher amount of CO2 adsorbed (4.1% of CO2 adsorbed in green) in Figure 5A with a vertical injector is due to the CO2 penetration in the caprock, which contained 22% of the mass injected. As discussed before, adsorption was much more relevant in the shale rock compared to in the sandstone aquifer minerals, and this was confirmed by Figure 6, with the inventory for each region (caprock and aquifer). Figure 6B shows that almost 11% of CO2 was adsorbed in the shale and that most of the CO2 in the caprock penetrated it by diffusion since CO2 dissolved in the brine was the highest amount (41.8%). Since there was no invasion of CO2 to caprock for the case with the horizontal well, the adsorbed quantity was expected to be lower, as shown in Figure 5.
The mass of mineralized CO2 was omitted from the previous graphs, because it was only 0.0005% of the total injected after 3000 years. It was mainly related to the calcite (83% in moles), dolomite (8% in moles), and pyrite (9% in moles) precipitation in the caprock with the vertical injector where CO2 contacted the caprock. There was no evidence of mineralization inside the aquifer since mainly non-reactive minerals constitute the sandstone matrix. Thus, those reactions are more relevant for modeling the caprock.
To compute how horizontal wells can make the CO2 flow less gravity-driven compared to vertical wells, the concept of the gravity number introduced by Zhou et al. [56] was applied, specified by Equation (20):
M N g v 1 + M
where
N g v = Δ ρ g L k v H u μ b r i n e , which is the characteristic time ratio for fluid to flow in the transverse direction due to gravity;
M is the mobility ratio;
Δ ρ is the density difference between the brine and CO2;
g is the acceleration due to gravity;
L is the length of the aquifer;
H is the aquifer thickness;
kv is the average vertical permeability;
u is the Darcy velocity (=interstitial velocity × φ);
μ b r i n e is the brine viscosity.
The higher the gravity number (Equation (20)) is, the more gravity-dominant the flow becomes. In our comparison, a value of 200 was calculated for a vertical well, whereas it was 100 for a horizontal well with 800 m of length, reinforcing the buoyancy observation in Figure 4.
The invasion depth in caprock due to diffusion reached over 26 m after 3000 years using a vertical well (based on the values of CO2 in the aqueous phase (CO2(aq)). The CO2 gas saturation plume (CO2(g)) in Figure 7 was only 17 m thick, highlighting that the diffusive flow was more relevant than Darcy flow in the caprock, as concluded by Currenti et al. [57] during the numerical investigation of CO2 leakage scenarios in shale caprock. Higher capillary pressure, critical water saturation typical of the shale (Figure 3), and lower permeability all inhibited the CO2 convective flow through the caprock.
This computed upward CO2 velocity of 26 m/3000 years had the same magnitude compared with other numerical studies, such as 3 m/100,000 years [36], 8 m/100,000 years [58], and 1 to 77 m/1000 years [57]. The numbers are unequal due to differences in the models’ properties and grid resolution and should be used only for qualitative comparison. In our case, a vertical refinement by a factor of 10 times was tested 50 m under and above the interface between sandstone and shale rocks using gridblocks of a 0.1 m thickness. With the finer gridblocks, the penetration was reduced by only 15% from the estimated penetration (22 m).
Those results show that CO2 penetration in the caprock is possible (mainly by diffusion) and can generate leakage to the surface or other non-target zones when the migration is associated with flow through fault/fracture potentially activated by the pressure buildup [57]. Figure 8 compares the CO2 propagation in terms of saturation and aqueous mole fraction when an open fracture with 10 mD of intrinsic permeability and an aperture of 0.00014 m was placed inside the caprock [57]. This represented a single natural fracture with a length of 1600 m which traversed the entire thickness of the caprock (it did not intersect the sandstone). It was characterized by the fixed intrinsic permeability and aperture for numerical flow simulation through it. As the dynamic and geomechanical behavior of the fracture (such as fracture propagation) was not considered in the modeling process, there was no requirement to specify the mechanical properties. The flow in the fracture was modeled using the Embedded Discrete Fracture Model (EDFM) [59,60,61,62] according to the procedure proposed by Machado et al. [63], which validated the EDFM application to reactive transport simulation in fractured media. In the current case, the dissolved CO2 in brine diffused into caprock and rose through the fracture due to its higher permeability than that of the shale caprock. Therefore, the penetration of the CO2(aq) (70 m) was also greater than by CO2(g) (17 m), as observed in Figure 7, making the leakage more likely when there were active fractures or faults in the shale caprock. To prevent or minimize this risk, the application of horizontal wells will be investigated in a heterogeneous and more realistic aquifer model in the next section.

5. Sensitivity Analysis with a Horizontal Injector

In the previous section, some advantages emerged regarding horizontal wells for a CCS project. In summary, they are shown as follows:
  • the wider lateral extension of the CO2 plume in comparison to the vertical rise;
  • more CO2 trapped as a residual fluid as it spreads in the porous media;
  • less gravity-dominant flow.
Consequently, the risk of contacting the caprock can be reduced, preventing leakage by diffusion or flow through geological features in the shale, such as fractures or faults.
To evaluate those points better, a thorough sensitivity analysis will be conducted in the heterogeneous aquifer model to compare the horizontal well performance and a vertical injector. As our goal mainly focuses on the flow inside the aquifer and preventing flow in the caprock, heterogeneity is more critical in the sandstone layers. Table 1 provides the base case properties for the heterogeneous aquifer, keeping the caprock permeability homogeneous. The pore volume injected (PVI) will be 6.5% over 25 years.

5.1. Sensitivity to the Horizontal Well Length (Lw)

Figure 9 brings the intuitive conclusion from the previous discussion on the lateral spreading of CO2 when a horizontal well was used. In this case, this lateral spreading became wider for a longer horizontal well. From Lw = 2000 m, after 3000 years of redistribution, there was no CO2 in contact with the caprock in terms of global mole fraction, which accounted for the total CO2 concentration in both brine and gas phases. Hence, the case with Lw = 2000 m was further considered in the base case with a horizontal well. This case showed no invasion of CO2 into the caprock, while the equivalent vertical injector had 2% brine with dissolved CO2 diffused into the caprock.

5.2. Sensitivity to the Injection Rate

The base case injected 6.5% PV at 350,000 metric tonnes/year as a typical field project injection rate. The CO2 mass injected was kept the same for these cases (~6.5% PVI), but the injection rate was varied. Figure 10 shows that the CO2 distribution was insensitive to the injection rate for long-duration evaluation since the total mass injected was the same and the long-term entrapment is a function of petrophysical and trapping parameters, which remained unchanged.

5.3. Sensitivity to the Horizontal Permeability

The permeability of the aquifer was increased by a factor of 10, keeping the ratio of kv/kh equal to 0.1. In both cases, the CO2 plume reached the caprock (Figure 11), emphasizing the main advantage of horizontal wells over vertical wells is in low- to mid-permeability reservoirs.

5.4. Sensitivity to the kv/kh Ratio

Assuming the properties of the base case (Table 1), the ratio of vertical to horizontal permeability was changed between 0.01 and 0.6. The results in Figure 12 confirm the lower risk of contacting the shale caprock when a horizontal well was used with a broad range of kv/kh ratios.

5.5. Sensitivity to the Natural Water Flow

Aquifers typically receive groundwater recharge, generating a natural flux in porous media. For the specific target represented by the heterogeneous base model, the natural water flow velocity was assumed to be 0.0015 ft/d [64]. To model this effect, two open boundary conditions were imposed on the lateral edges of the aquifer, coupled with an analytical function to calculate the water influx to the aquifer using the Carter−Tracy method [65]. The left edge was modeled as the source, and the right edge was modeled as the sink, as represented in Figure 13. Calibrating the analytical aquifer permeability, it was possible to calculate a similar velocity of 0.0015 ft/d (yellow region in Figure 13).
Despite the low velocity, this flux affected the CO2 plume propagation. Increased plume spreading makes the monitoring more challenging and impossible to prevent contact with the caprock even with a horizontal well. Table 7 shows the lateral spread of the CO2(g) plume and the thickness inside the shale from Figure 14 and the percentage of brine mass with CO2(aq) inside the caprock from Figure 15.
It is clear that the natural water flux contributed to expanding the CO2 lateral extension, limiting its vertical migration, e.g., reducing the leakage risk but making it more difficult for plume monitoring with observation wells due to its spread over thousands of meters and with seismic methods due to its small thickness (the thickness of less than 10 m was challenging due to the seismic vertical resolution).
However, this contact with the shale overburden using the vertical well case was still more pronounced compared to that using the horizontal well, which further increased the risk of flowing via existing or activated fractures and faults in the caprock.

5.6. Sensitivity to the Time of Redistribution

All of the previous analyses show that horizontal wells were more resilient than vertical ones in preventing contact with the caprock and consequently avoided CO2 leakage through paths that existed or could be created inside the shale. Figure 16B expands the timeframe of the base case over 20,000 years and confirms no contact with the caprock when a horizontal well with Lw = 2000 m was considered in the base case. The equivalent plot for the vertical well in Figure 16A shows that the CO2 diffusion into the caprock reached 3% of the mass of the brine with CO2(aq) at 20,000 years, against 2% at 3000 years.

6. Horizontal Well Design

The study identified two primary findings associated with the use of horizontal wells in CO2 storage:
  • Horizontal injection wells effectively prevent or minimize CO2 penetration into the caprock across various sensitivity scenarios. This finding is closely linked with the buoyancy effect, which can be modeled by the gravity number introduced in Equation (20). Consequently, the gravity number is one of the variables used to determine the optimal length of horizontal wells;
  • Horizontal wells offer a safer approach to trapping CO2 by increasing its entrapment as a residual phase compared to the use of vertical wells as injectors. In this context, the additional CO2 trapped as a residual phase becomes a crucial variable.
Therefore, the ideal length for a horizontal well would be one that maximizes CO2 entrapment while minimizing buoyancy effects.
To confirm our hypothesis that the wider the lateral CO2 spreading, the more the residual entrapment, Figure 17 brings the evolution of the residual CO2 trapped by hysteresis for the cases considered in Figure 9. Figure 17 shows that a longer Lw increased residual saturations as the plume contacted new permeability regions of the aquifer. However, the residual CO2 mass developed a plateau after 1500 years for all these cases, which represented the optimum value. For instance, from Lw = 2000 m, there was no additional entrapment for longer horizontal branches. Therefore, Lw = 2000 m seems to be an optimum length for the entrapment.
Figure 18 combines two parameters plotted against the length of the horizontal well:
  • On the left: the difference between the CO2 saturation trapped as a residual phase (∆CO2 residual) obtained with the vertical well minus the one with the horizontal well after 1500 years, corresponding to the time of plateau in Figure 17;
  • On the right: the gravity numbers according to Equation (20) for different well lengths.
Due to the permeability heterogeneity, a well length of Lw = 2000 m had more access to higher permeability regions (i.e., concentrated in the middle of this model). Hence, it developed higher z-transverse velocities (u in Equation (20)), causing a lower gravity number for the well length of 2000 m (Figure 18).
This combination allows setting an optimum horizontal well length to enhance the CO2 entrapment and to control the rise of the CO2 plume, preventing its contact with the caprock and mitigating the associated risks. In this case, the optimum considering both purposes was Lw = 2000 m, which is technically possible according to actual CCS field projects [18,25].
Figure 18. The difference between the residual trapped CO2 for vertical and horizontal wells (axis on the left) and the gravity number (axis on the right) against the well length after 1500 years.
Figure 18. The difference between the residual trapped CO2 for vertical and horizontal wells (axis on the left) and the gravity number (axis on the right) against the well length after 1500 years.
Applsci 13 12830 g018

7. Conclusions

The main conclusions are summarized below:
  • The CO2 intrusion into the caprock was mainly dominated by diffusion of the dissolved CO2 in the brine. This can become more critical, if this solution accesses active faults and fractures of the shale caprock;
  • The impact of mineral precipitation in the caprock and its integrity was negligible, with only an insignificant precipitation of about 0.0005% of the total injected CO2 after 3000 years. Therefore, considering reservoir parameters such as mineralogy, pH, injection rate, pressure, and temperature, there was a small risk for caprock integrity due to mineralization/dissolution;
  • Horizontal wells were more effective in controlling the CO2 buoyancy due to two mechanisms:
    enhancement of the CO2 entrapment as a residual phase by up to 19% due to relative permeability/capillary pressure hysteresis where the more spread plume contacted new portions of the aquifer pore volume;
    the better balance between gravity and viscous forces provided by horizontal wells compared to by vertical wells, with up to an 18% reduction in this balance compared to the vertical well case. The conclusion was based on simulations of homogeneous and heterogeneous sandstone saline aquifers over thousands of years.
  • A sensitivity test showed the following results:
    hypothesis was mainly valid for low- and mid-permeability aquifers (<200 mD);
    for long-term evaluation, there was no significant impact on the injection rate for the same CO2 mass injected;
    the natural flux in aquifers can affect the plume propagation, triggering its contact with the caprock even with horizontal wells; and
    these conclusions regarding horizontal wells persisted over tens of thousands of years.
  • A practical method was proposed to design the optimum length for horizontal wells, combining the maximum CO2 saturation trapped as the residual phase and the gravity number.
These results present an exciting opportunity to further enhance CO2 entrapment in CCS projects. The study’s findings provide a basis for exploring alternative project options. These alternatives may involve optimizing well geometry, considering uncertainties in trapping properties and exploring different geological realizations.
By incorporating well geometry optimization techniques and accounting for uncertainties, it is possible to improve CO2 entrapment and storage efficiency. This approach allows for a more comprehensive assessment of potential benefits and risks associated and opens up promising perspectives for advancing CCS projects.

Author Contributions

Conceptualization, M.V.B.M., M.D. and K.S.; methodology, M.V.B.M.; software, M.V.B.M.; validation, M.V.B.M., M.D. and K.S.; formal analysis, M.V.B.M., M.D. and K.S.; investigation, M.V.B.M.; resources, M.V.B.M. and M.D.; data curation, M.V.B.M. and M.D.; writing—original draft preparation, M.V.B.M.; writing—review and editing, M.D. and K.S.; visualization, M.V.B.M.; supervision, M.D. and K.S.; project administration, M.D. and K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used in this study are provided within the article. However, the models used are not accessible to the public due to PETROBRAS’ policy.

Acknowledgments

The authors would like to thank PETROBRAS for authorizing the use of models and Chowdhury Mamun for thorough English revision.

Conflicts of Interest

Author Marcos Vitor Barbosa Machado is employed by the company PETROBRAS. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Nomenclature

B C O 2 parameter for Langmuir isotherm relation, kPa−1;
CLand’s constant, dimensionless;
cfrock compressibility, kPa−1;
Ddiffusion coefficient, cm2/s;
Deffeffective diffusion coefficient, cm2/s;
gacceleration due to gravity, m/s2;
Haquifer thickness, m;
H i Henry’s constant at the current pressure (p) and temperature (T), dimensionless;
H i * Henry’s constant at the reference pressure (p*) and temperature (T), dimensionless;
JLeverett J-function, dimensionless;
k or khaverage horizontal permeability, mD [9.869 × 10−16 m2];
k v average vertical permeability, mD [9.869 × 10−16 m2];
krlrelative permeability, dimensionless;
Llength of the aquifer, m;
Lwhorizontal well length, m;
Mmobility ratio, dimensionless;
Njthe total moles of mineral j, gmol/m3;
N g v characteristic time ratio for fluid to flow in the transverse direction due to gravity, dimensionless;
ppressure, kPa;
PcCO2−brine capillary pressure, kPa;
Runiversal gas constant, 8.314 kPa·L/mol·K;
rfresistance factor, dimensionless;
Sgttrapped gas saturation, dimensionless;
Sg maxmaximum gas saturation, dimensionless;
Ttemperature, °C;
uthe Darcy velocity (real velocity × φ), m/s;
y C O 2 ,   g molar fraction of adsorbed CO2 in the gas phase, dimensionless;
Zglobal mole fraction, dimensionless.
Greek symbols
φrock porosity, fraction;
μ b r i n e brine viscosity, cP [10−3 Pa.s];
ρmmineral molar density, gmol/m3;
ρdensity, kg/m3;
τtortuosity, dimensionless;
v i ¯ partial molar volume at infinite dilution, L/mol;
ω C O 2 moles of adsorbed CO2 per unit mass of rock, gmole/kg of rock;
ω C O 2 ,   m a x maximum moles of adsorbed CO2 per unit mass of rock, gmole/kg of rock.

References

  1. Birol, D.F. World Energy Outlook 2022; IEA Publications: Paris, France, 2022. [Google Scholar]
  2. Nghiem, L.; Shrivastava, V.; Kohse, B.; Hassam, M.; Yang, C. Simulation of Trapping Processes for CO2 Storage in Saline Aquifers. In Proceedings of the Canadian International Petroleum Conference, Calgary, AB, Canada, 16–18 June 2009; Petroleum Society of Canada: Calgary, AB, Canada, 2009. [Google Scholar]
  3. Han, W.S.; McPherson, B.J.; Lichtner, P.C.; Wang, F.P. Evaluation of Trapping Mechanisms in Geologic CO2 Sequestration: Case Study of SACROC Northern Platform, a 35-Year CO2 Injection Site. Am. J. Sci. 2010, 310, 282–324. [Google Scholar] [CrossRef]
  4. Delshad, M.; Kong, X.; Tavakoli, R.; Hosseini, S.A.; Wheeler, M.F. Modeling and Simulation of Carbon Sequestration at Cranfield Incorporating New Physical Models. Int. J. Greenh. Gas Control 2013, 18, 463–473. [Google Scholar] [CrossRef]
  5. Rackley, S.A. Carbon Capture and Storage, 2nd ed.; Butterworth-Heinemann: Cambridge, MA, USA, 2017; ISBN 978-0-12-812041-5. [Google Scholar]
  6. Hovorka, S. Optimization of Geological Environments for Carbon Dioxide Disposal in Saline Aquifers in the United States (Part One); University of Texas: Austin, TX, USA, 2008; p. 990445. [Google Scholar]
  7. Bump, A.P.; Bakhshian, S.; Ni, H.; Hovorka, S.D.; Olariu, M.I.; Dunlap, D.; Hosseini, S.A.; Meckel, T.A. Composite Confining Systems: Rethinking Geologic Seals for Permanent CO2 Sequestration. Int. J. Greenh. Gas Control 2023, 126, 103908. [Google Scholar] [CrossRef]
  8. Duan, Z.; Sun, R. An Improved Model Calculating CO2 Solubility in Pure Water and Aqueous NaCl Solutions from 273 to 533 K and from 0 to 2000 Bar. Chem. Geol. 2003, 193, 257–271. [Google Scholar] [CrossRef]
  9. Portier, S.; Rochelle, C. Modelling CO2 Solubility in Pure Water and NaCl-Type Waters from 0 to 300 °C and from 1 to 300 Bar. Chem. Geol. 2005, 217, 187–199. [Google Scholar] [CrossRef]
  10. Spiteri, E.J.; Juanes, R.; Blunt, M.J.; Orr, F.M. Relative Permeability Hysteresis: Trapping Models and Application to Geological CO2 Sequestration. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dallas, TX, USA, 9 October 2005; p. SPE-96448-MS. [Google Scholar]
  11. Nghiem, L.; Yang, C.; Shrivastava, V.; Kohse, B.; Hassam, M.; Card, C. Risk Mitigation through the Optimization of Residual Gas and Solubility Trapping for CO2 Storage in Saline Aquifers. Energy Procedia 2009, 1, 3015–3022. [Google Scholar] [CrossRef]
  12. Qi, R.; Laforce, T.; Blunt, M. Design of Carbon Dioxide Storage in Aquifers. Int. J. Greenh. Gas Control 2009, 3, 195–205. [Google Scholar] [CrossRef]
  13. Xu, T.; Yue, G.; Wang, F.; Liu, N. Using Natural CO2 Reservoir to Constrain Geochemical Models for CO2 Geological Sequestration. Appl. Geochem. 2014, 43, 22–34. [Google Scholar] [CrossRef]
  14. Gunter, W.D.; Bachu, S.; Benson, S. The Role of Hydrogeological and Geochemical Trapping in Sedimentary Basins for Secure Geological Storage of Carbon Dioxide. Geol. Soc. Lond. Spéc. Publ. 2004, 233, 129–145. [Google Scholar] [CrossRef]
  15. Han, W.S.; Kim, K.-Y.; Esser, R.P.; Park, E.; McPherson, B.J. Sensitivity Study of Simulation Parameters Controlling CO2 Trapping Mechanisms in Saline Formations. Transp. Porous. Med. 2011, 90, 807–829. [Google Scholar] [CrossRef]
  16. Han, W.S.; Lee, S.-Y.; Lu, C.; McPherson, B.J. Effects of Permeability on CO2 Trapping Mechanisms and Buoyancy-Driven CO2 Migration in Saline Formations. Water Resour. Res. 2010, 46, W07510. [Google Scholar] [CrossRef]
  17. Okwen, R.; Stewart, M.; Cunningham, J. Effect of Well Orientation (Vertical vs. Horizontal) and Well Length on the Injection of CO2 in Deep Saline Aquifers. Transp. Porous. Med. 2011, 90, 219–232. [Google Scholar] [CrossRef]
  18. Eiken, O.; Ringrose, P.; Hermanrud, C.; Nazarian, B.; Torp, T.A.; Høier, L. Lessons Learned from 14 Years of CCS Operations: Sleipner, In Salah and Snøhvit. Energy Procedia 2011, 4, 5541–5548. [Google Scholar] [CrossRef]
  19. Esposito, A.; Benson, S.M. Evaluation and Development of Options for Remediation of CO2 Leakage into Groundwater Aquifers from Geologic Carbon Storage. Int. J. Greenh. Gas Control 2012, 7, 62–73. [Google Scholar] [CrossRef]
  20. Vialle, S.; Druhan, J.L.; Maher, K. Multi-Phase Flow Simulation of CO2 Leakage through a Fractured Caprock in Response to Mitigation Strategies. Int. J. Greenh. Gas Control 2016, 44, 11–25. [Google Scholar] [CrossRef]
  21. Sharma, S.; Agrawal, V.; McGrath, S.; Hakala, J.A.; Lopano, C.; Goodman, A. Geochemical Controls on CO2 Interactions with Deep Subsurface Shales: Implications for Geologic Carbon Sequestration. Environ. Sci. Processes Impacts 2021, 23, 1278–1300. [Google Scholar] [CrossRef] [PubMed]
  22. Machado, M.V.B.; Delshad, M.; Sepehrnoori, K. Injectivity Assessment for CCS Field-Scale Projects with Considerations of Salt Deposition, Mineral Dissolution, Fines Migration, Hydrate Formation, and Non-Darcy Flow. Fuel 2023, 353, 129148. [Google Scholar] [CrossRef]
  23. Gholami, R.; Raza, A.; Iglauer, S. Leakage Risk Assessment of a CO2 Storage Site: A Review. Earth-Sci. Rev. 2021, 223, 103849. [Google Scholar] [CrossRef]
  24. Buscheck, T.A.; Sun, Y.; Chen, M.; Hao, Y.; Wolery, T.J.; Bourcier, W.L.; Court, B.; Celia, M.A.; Julio Friedmann, S.; Aines, R.D. Active CO2 Reservoir Management for Carbon Storage: Analysis of Operational Strategies to Relieve Pressure Buildup and Improve Injectivity. Int. J. Greenh. Gas Control 2012, 6, 230–245. [Google Scholar] [CrossRef]
  25. Kim, M.; Kwon, S.; Ji, M.; Shin, H.; Min, B. Multi-Lateral Horizontal Well with Dual-Tubing System to Improve CO2 Storage Security and Reduce CCS Cost. Appl. Energy 2023, 330, 120368. [Google Scholar] [CrossRef]
  26. Machado, M.V.B.; Delshad, M.; Sepehrnoori, K. A Practical and Innovative Workflow to Support the Numerical Simulation of CO2 Storage in Large Field-Scale Models. SPE Reserv. Eval. Eng. 2023, 26, 1541–1552. [Google Scholar] [CrossRef]
  27. Bennion, D.B.; Bachu, S. Drainage and Imbibition Relative Permeability Relationships for Supercritical CO2/Brine and H2S/Brine Systems in Intergranular Sandstone, Carbonate, Shale, and Anhydrite Rocks. SPE Reserv. Eval. Eng. 2008, 11, 487–496. [Google Scholar] [CrossRef]
  28. Abdoulghafour, H.; Sarmadivaleh, M.; Hauge, L.P.; Fernø, M.; Iglauer, S. Capillary Pressure Characteristics of CO2-Brine-Sandstone Systems. Int. J. Greenh. Gas Control 2020, 94, 102876. [Google Scholar] [CrossRef]
  29. Bennion, D.B.; Bachu, S. Permeability and Relative Permeability Measurements at Reservoir Conditions for CO2-Water Systems in Ultralow-Permeability Confining Caprocks. In Proceedings of the EUROPEC/EAGE Conference and Exhibition, London, UK, 11 June 2007; p. SPE-106995-MS. [Google Scholar]
  30. CMG. GEM Compositional & Unconventional Simulator, (version 2022.10); Windows, CMG: Calgary, AB, Canada, 2022. [Google Scholar]
  31. Balhoff, M. An Introduction to Multiphase, Multicomponent Reservoir Simulation; Developments in Petroleum Science; Elsevier: Amsterdam, The Netherlands, 2022; ISBN 978-0-323-99235-0. [Google Scholar]
  32. Machado, M.V.B. Numerical Petroleum Reservoir Modeling: Integrated Simulation Practice, 1st ed.; PETROBRAS: Rio de Janeiro, Brazil, 2023; ISBN 9786588763070. (In Portuguese) [Google Scholar]
  33. Li, Y.-K.; Nghiem, L.X. Phase Equilibria of Oil, Gas and Water/Brine Mixtures from a Cubic Equation of State and Henry’s Law. Can. J. Chem. Eng. 1986, 64, 486–496. [Google Scholar] [CrossRef]
  34. Bakker, R.J. Package FLUIDS 1. Computer Programs for Analysis of Fluid Inclusion Data and for Modelling Bulk Fluid Properties. Chem. Geol. 2003, 194, 3–23. [Google Scholar] [CrossRef]
  35. Rezk, M.G.; Foroozesh, J.; Abdulrahman, A.; Gholinezhad, J. CO2 Diffusion and Dispersion in Porous Media: Review of Advances in Experimental Measurements and Mathematical Models. Energy Fuels 2022, 36, 133–155. [Google Scholar] [CrossRef]
  36. Busch, A.; Alles, S.; Gensterblum, Y.; Prinz, D.; Dewhurst, D.N.; Raven, M.D.; Stanjek, H.; Krooss, B.M. Carbon Dioxide Storage Potential of Shales. Int. J. Greenh. Gas Control 2008, 2, 297–308. [Google Scholar] [CrossRef]
  37. Montegrossi, G.; Cantucci, B.; Piochi, M.; Fusi, L.; Misnan, M.S.; Rashidi, M.R.A.; Abu Bakar, Z.A.; Tuan Harith, Z.Z.; Bahri, N.H.S.; Hashim, N. CO2 Reaction-Diffusion Experiments in Shales and Carbonates. Minerals 2022, 13, 56. [Google Scholar] [CrossRef]
  38. Langmuir, I. The Constitution and Fundamental Properties of Solids and Liquids. J. Frankl. Inst. 1917, 183, 102–105. [Google Scholar] [CrossRef]
  39. Xie, W.; Wang, M.; Wang, H. Adsorption Characteristics of CH4 and CO2 in Shale at High Pressure and Temperature. ACS Omega 2021, 6, 18527–18536. [Google Scholar] [CrossRef]
  40. Wang, J.; Samara, H.; Jaeger, P.; Ko, V.; Rodgers, D.; Ryan, D. Investigation for CO2 Adsorption and Wettability of Reservoir Rocks. Energy Fuels 2022, 36, 1626–1634. [Google Scholar] [CrossRef]
  41. Ambrose, R.J.; Hartman, R.C.; Akkutlu, I.Y. Multi-Component Sorbed-Phase Considerations for Shale Gas-in-Place Calculations. In Proceedings of the SPE Production and Operations Symposium, Oklahoma City, OK, USA, 27 March 2011; p. SPE-141416-MS. [Google Scholar]
  42. Hartman, R.C.; Ambrose, R.J.; Akkutlu, I.Y.; Clarkson, C.R. Shale Gas-in-Place Calculations Part II—Multi-Component Gas Adsorption Effects. In Proceedings of the North American Unconventional Gas Conference and Exhibition, The Woodlands, TX, USA, 14 June 2011; p. SPE-144097-MS. [Google Scholar]
  43. Ding, J.; Yan, C.; Wang, G.; He, Y.; Zhao, R. Competitive Adsorption Between CO2 and CH4 in Tight Sandstone and Its Influence on CO2-Injection Enhanced Gas Recovery (EGR). Int. J. Greenh. Gas Control 2022, 113, 103530. [Google Scholar] [CrossRef]
  44. Arri, L.E.; Yee, D.; Morgan, W.D.; Jeansonne, M.W. Modeling Coalbed Methane Production with Binary Gas Sorption. In Proceedings of the SPE Rocky Mountain Regional Meeting, Casper, WY, USA, 18 May 1992; p. SPE-24363-MS. [Google Scholar]
  45. Hall, F.E.; Zhou, C.; Gasem, K.A.M.; Robinson, R.L.; Yee, D. Adsorption of Pure Methane, Nitrogen, and Carbon Dioxide and Their Binary Mixtures on Wet Fruitland Coal. In Proceedings of the SPE Eastern Regional Meeting, Charleston, WV, USA, 8 November 1994; p. SPE-29194-MS. [Google Scholar]
  46. Land, C.S. Calculation of Imbibition Relative Permeability for Two- and Three-Phase Flow from Rock Properties. Soc. Pet. Eng. J. 1968, 8, 149–156. [Google Scholar] [CrossRef]
  47. Carlson, F.M. Simulation of Relative Permeability Hysteresis to the Nonwetting Phase. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 4 October 1981; p. SPE-10157-MS. [Google Scholar]
  48. Jarrell, P.M.; Fox, C.; Stein, M.; Webb, S. Practical Aspects of CO2 Flooding; SPE Monograph Series; Henry L. Doherty Memorial Fund of AIME; Society of Petroleum Engineers: Richardson, TX, USA, 2002; ISBN 978-1-55563-096-6. [Google Scholar]
  49. Burnside, N.M.; Naylor, M. Review and Implications of Relative Permeability of CO2/Brine Systems and Residual Trapping of CO2. Int. J. Greenh. Gas Control 2014, 23, 1–11. [Google Scholar] [CrossRef]
  50. Parkhurst, D.L.; Thorstenson, D.C.; Plummer, L.N. PHREEQE: A Computer Program for Geochemical Calculations; U.S. Geological Survey: Denver, CO, USA, 1980.
  51. Parkhurst, D.L.; Appelo, C.A.J. Description of Input and Examples for PHREEQC Version 3: A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations. In U.S. Geological Survey Techniques and Methods, Book 6; U.S. Geological Survey: Denver, CO, USA, 2013; p. 497. [Google Scholar]
  52. Zeng, L.; Vialle, S.; Ennis-King, J.; Esteban, L.; Sarmadivaleh, M.; Sarout, J.; Dautriat, J.; Giwelli, A.; Xie, Q. Role of Geochemical Reactions on Caprock Integrity during Underground Hydrogen Storage. J. Energy Storage 2023, 65, 107414. [Google Scholar] [CrossRef]
  53. Xiao, Y.; Xu, T.; Pruess, K. The Effects of Gas-Fluid-Rock Interactions on CO2 Injection and Storage: Insights from Reactive Transport Modeling. Energy Procedia 2009, 1, 1783–1790. [Google Scholar] [CrossRef]
  54. Wan, J.; Tokunaga, T.K.; Ashby, P.D.; Kim, Y.; Voltolini, M.; Gilbert, B.; DePaolo, D.J. Supercritical CO2 Uptake by Nonswelling Phyllosilicates. Proc. Natl. Acad. Sci. USA 2018, 115, 873–878. [Google Scholar] [CrossRef]
  55. Zeidouni, M.; Pooladi-Darvish, M.; Keith, D. Analytical Solution to Evaluate Salt Precipitation during CO2 Injection in Saline Aquifers. Int. J. Greenh. Gas Control 2009, 3, 600–611. [Google Scholar] [CrossRef]
  56. Zhou, D.; Fayers, F.J.; Orr, F.M. Scaling of Multiphase Flow in Simple Heterogeneous Porous Media. In Proceedings of the SPE/DOE Improved Oil Recovery Symposium, Tulsa, OK, USA, 17 April 1994; p. SPE-27833-MS. [Google Scholar]
  57. Currenti, G.; Cantucci, B.; Montegrossi, G.; Napoli, R.; Misnan, M.S.; Rashidi, M.R.A.; Abu Bakar, Z.A.; Harith, Z.Z.T.; Bahri, N.H.S.; Hashim, N. CO2 Leakage Scenarios in Shale Overburden. Minerals 2023, 13, 1016. [Google Scholar] [CrossRef]
  58. Kampman, N.; Busch, A.; Bertier, P.; Snippe, J.; Hangx, S.; Pipich, V.; Di, Z.; Rother, G.; Harrington, J.F.; Evans, J.P.; et al. Observational Evidence Confirms Modelling of the Long-Term Integrity of CO2-Reservoir Caprocks. Nat. Commun. 2016, 7, 12268. [Google Scholar] [CrossRef] [PubMed]
  59. Li, L.; Lee, S.H. Efficient Field-Scale Simulation of Black Oil in a Naturally Fractured Reservoir through Discrete Fracture Networks and Homogenized Media. SPE Reserv. Eval. Eng. 2008, 11, 750–758. [Google Scholar] [CrossRef]
  60. Moinfar, A.; Varavei, A.; Sepehrnoori, K.; Johns, R.T. Development of an Efficient Embedded Discrete Fracture Model for 3D Compositional Reservoir Simulation in Fractured Reservoirs. SPE J. 2014, 19, 289–303. [Google Scholar] [CrossRef]
  61. Cavalcante Filho, J.S.d.A.; Shakiba, M.; Moinfar, A.; Sepehrnoori, K. Implementation of a Preprocessor for Embedded Discrete Fracture Modeling in an IMPEC Compositional Reservoir Simulator. In Proceedings of the SPE Reservoir Simulation Symposium, Houston, TX, USA, 23 February 2015; p. D012S021R013. [Google Scholar]
  62. Sepehrnoori, K.; Xu, Y.; Yu, W. Embedded Discrete Fracture Modeling and Application in Reservoir Simulation; Developments in Petroleum Science; Elsevier: Amsterdam, The Netherlands, 2020; ISBN 978-0-12-821872-3. [Google Scholar]
  63. Machado, M.V.B.; Delshad, M.; Sepehrnoori, K. A Computationally Efficient Approach to Model Reactive Transport during CO2 Storage in Naturally Fractured Saline Aquifers. SSRN 2023, 4520437. [Google Scholar] [CrossRef]
  64. Correia, L.C.; Alves, M.D.G.; Silva Júnior, G.C.D. Estimativa Da Recarga de Água Subterrânea Utilizando o Método WTF Na Porção Continental Da Bacia Sedimentar de Campos, Rio de Janeiro, Brasil. R. Águas Subter. 2021, 35, 1–15. [Google Scholar] [CrossRef]
  65. Carter, R.D.; Tracy, G.W. An Improved Method for Calculating Water Influx. Trans. AIME 1960, 219, 415–417. [Google Scholar] [CrossRef]
Figure 1. Models studied in this work: (A) a homogeneous sandstone aquifer; (B) a heterogeneous sandstone aquifer. Coordinates are in meters. Blue areas in both images represent the shale caprock.
Figure 1. Models studied in this work: (A) a homogeneous sandstone aquifer; (B) a heterogeneous sandstone aquifer. Coordinates are in meters. Blue areas in both images represent the shale caprock.
Applsci 13 12830 g001
Figure 2. Drainage relative permeability (left) and capillary pressure (right) for the sandstone.
Figure 2. Drainage relative permeability (left) and capillary pressure (right) for the sandstone.
Applsci 13 12830 g002
Figure 3. Drainage relative permeability (left) and capillary pressure (right) for the shale caprock.
Figure 3. Drainage relative permeability (left) and capillary pressure (right) for the shale caprock.
Applsci 13 12830 g003
Figure 4. Cross-sectional view of CO2 redistribution after 3000 years through a vertical well (A), a horizontal well with Lw = 400 m (B), and a horizontal well with Lw = 800 m (C). The dashed black line marks the top of the sandstone aquifer.
Figure 4. Cross-sectional view of CO2 redistribution after 3000 years through a vertical well (A), a horizontal well with Lw = 400 m (B), and a horizontal well with Lw = 800 m (C). The dashed black line marks the top of the sandstone aquifer.
Applsci 13 12830 g004
Figure 5. CO2 inventory after 3000 years through a vertical well (A) and a horizontal well with Lw = 800 m (B).
Figure 5. CO2 inventory after 3000 years through a vertical well (A) and a horizontal well with Lw = 800 m (B).
Applsci 13 12830 g005
Figure 6. CO2 inventory after 3000 years in the sandstone aquifer (A) and in the shale caprock (B) with a vertical injector.
Figure 6. CO2 inventory after 3000 years in the sandstone aquifer (A) and in the shale caprock (B) with a vertical injector.
Applsci 13 12830 g006
Figure 7. Profiles along the vertical well after 3000 years of CO2 saturation (A) and aqueous CO2 mole fraction (B), inside the caprock.
Figure 7. Profiles along the vertical well after 3000 years of CO2 saturation (A) and aqueous CO2 mole fraction (B), inside the caprock.
Applsci 13 12830 g007
Figure 8. Cross-sectional view of CO2 redistribution, injected through a vertical well, after 3000 years represented by CO2 saturation (A) and CO2 aqueous mole fraction (B). The dashed black line marks an open natural fracture in the caprock.
Figure 8. Cross-sectional view of CO2 redistribution, injected through a vertical well, after 3000 years represented by CO2 saturation (A) and CO2 aqueous mole fraction (B). The dashed black line marks an open natural fracture in the caprock.
Applsci 13 12830 g008
Figure 9. Cross-sectional views of CO2 global mole fraction after 3000 years for different well geometries. From left to right: vertical injector and horizontal wells with Lw = 1000, 2000, 3000, and 4000 m. The dashed black line marks the boundary between the aquifer and the caprock; the coordinates are in meters.
Figure 9. Cross-sectional views of CO2 global mole fraction after 3000 years for different well geometries. From left to right: vertical injector and horizontal wells with Lw = 1000, 2000, 3000, and 4000 m. The dashed black line marks the boundary between the aquifer and the caprock; the coordinates are in meters.
Applsci 13 12830 g009
Figure 10. Cross-sectional views of CO2 global mole fraction after 3000 years for vertical and horizontal wells for the base case at 350,000 tonnes/year (A) and the injection rate reduced by a factor of 2 (B). The light blue zone is the caprock; the coordinates are in meters.
Figure 10. Cross-sectional views of CO2 global mole fraction after 3000 years for vertical and horizontal wells for the base case at 350,000 tonnes/year (A) and the injection rate reduced by a factor of 2 (B). The light blue zone is the caprock; the coordinates are in meters.
Applsci 13 12830 g010
Figure 11. Cross-sectional views of CO2 global mole fraction after 3000 years comparing vertical and horizontal wells for the base case with average horizontal permeability of 140 mD (A) and with permeability increase by a factor of 10 (B). The light blue zone is the caprock; the coordinates are in meters.
Figure 11. Cross-sectional views of CO2 global mole fraction after 3000 years comparing vertical and horizontal wells for the base case with average horizontal permeability of 140 mD (A) and with permeability increase by a factor of 10 (B). The light blue zone is the caprock; the coordinates are in meters.
Applsci 13 12830 g011
Figure 12. Cross-sectional views of the CO2 global mole fraction (Z) after 3000 years with a vertical injector (A) and a horizontal injector of Lw = 2000 m (B) for different ratios of kv/kh. The light blue zone is the caprock; the coordinates are in meters.
Figure 12. Cross-sectional views of the CO2 global mole fraction (Z) after 3000 years with a vertical injector (A) and a horizontal injector of Lw = 2000 m (B) for different ratios of kv/kh. The light blue zone is the caprock; the coordinates are in meters.
Applsci 13 12830 g012
Figure 13. Modeling a natural water influx into the aquifer using open lateral boundaries. The colors represent the water velocity in ft/day; positive values are the influx, and negative values are the efflux. The caprock is in yellow, because the velocity is zero.
Figure 13. Modeling a natural water influx into the aquifer using open lateral boundaries. The colors represent the water velocity in ft/day; positive values are the influx, and negative values are the efflux. The caprock is in yellow, because the velocity is zero.
Applsci 13 12830 g013
Figure 14. CO2 saturation filtered only inside the caprock region after 3000 years considering a vertical injector (A), a vertical injector and water influx (B), and a horizontal injector with water influx (C).
Figure 14. CO2 saturation filtered only inside the caprock region after 3000 years considering a vertical injector (A), a vertical injector and water influx (B), and a horizontal injector with water influx (C).
Applsci 13 12830 g014
Figure 15. Cross-sectional views of the CO2 global mole fraction after 3000 years for the base case without water influx (A) and with water influx (B). The dashed white line marks the boundary of the aquifer and caprock; the coordinates are in meters.
Figure 15. Cross-sectional views of the CO2 global mole fraction after 3000 years for the base case without water influx (A) and with water influx (B). The dashed white line marks the boundary of the aquifer and caprock; the coordinates are in meters.
Applsci 13 12830 g015
Figure 16. Cross-sectional views of the CO2 global mole fraction after 20,000 years with a vertical well (A) and a horizontal well with Lw = 2000 m (B). The light blue zone is the caprock; the coordinates are in meters.
Figure 16. Cross-sectional views of the CO2 global mole fraction after 20,000 years with a vertical well (A) and a horizontal well with Lw = 2000 m (B). The light blue zone is the caprock; the coordinates are in meters.
Applsci 13 12830 g016
Figure 17. The masses of residual CO2 trapped by relative permeability and capillary pressure hysteresis for different well geometries over 3000 years.
Figure 17. The masses of residual CO2 trapped by relative permeability and capillary pressure hysteresis for different well geometries over 3000 years.
Applsci 13 12830 g017
Table 1. Summary of the main petrophysical properties used in the geological models.
Table 1. Summary of the main petrophysical properties used in the geological models.
Homogeneous ModelHeterogeneous Model
SandstoneShaleSandstoneShale
porosity (φ)0.150.100.21 (mean)0.10
permeability (k)100 mD0.001 mD140 mD (mean)0.001 mD
kv/kh ratio0.10.10.10.1
pore compressibility5.8 × 10−7 kPa−15 × 10−8 kPa−15.8 × 10−7 kPa−15 × 10−8 kPa−1
relative permeabilityFigure 2Figure 3Figure 2Figure 3
capillary pressureFigure 2Figure 3Figure 2Figure 3
Table 2. Summary of the initial conditions of the reservoir models.
Table 2. Summary of the initial conditions of the reservoir models.
Homogeneous Saline AquiferHeterogeneous Saline Aquifer
Initial pressure11,800 kPa at 1000 m7159 kPa at 730 m
Temperature80 °C70 °C
Salinity50,000 ppm70,000 ppm
H i * 5.66 × 1056.44 × 105
Table 4. Maximum gas trapped due to hysteresis of relative permeabilities and capillarity.
Table 4. Maximum gas trapped due to hysteresis of relative permeabilities and capillarity.
SandstoneShale
Sgt0.250.35
Table 7. Lateral CO2 plume extension in the aquifer, thickness in the caprock and the diffused mass of brine mass saturated with CO2.
Table 7. Lateral CO2 plume extension in the aquifer, thickness in the caprock and the diffused mass of brine mass saturated with CO2.
Aquifer CO2 Plume Spread (m)Caprock CO2 Plume Thickness (m)Brine in the Caprock
(% in Mass)
vertical injector2980–5510 36−62.0−1.1
horizontal injector with a 2000 m length0–1440 0–60.0–0.26
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Machado, M.V.B.; Delshad, M.; Sepehrnoori, K. Potential Benefits of Horizontal Wells for CO2 Injection to Enhance Storage Security and Reduce Leakage Risks. Appl. Sci. 2023, 13, 12830. https://doi.org/10.3390/app132312830

AMA Style

Machado MVB, Delshad M, Sepehrnoori K. Potential Benefits of Horizontal Wells for CO2 Injection to Enhance Storage Security and Reduce Leakage Risks. Applied Sciences. 2023; 13(23):12830. https://doi.org/10.3390/app132312830

Chicago/Turabian Style

Machado, Marcos Vitor Barbosa, Mojdeh Delshad, and Kamy Sepehrnoori. 2023. "Potential Benefits of Horizontal Wells for CO2 Injection to Enhance Storage Security and Reduce Leakage Risks" Applied Sciences 13, no. 23: 12830. https://doi.org/10.3390/app132312830

APA Style

Machado, M. V. B., Delshad, M., & Sepehrnoori, K. (2023). Potential Benefits of Horizontal Wells for CO2 Injection to Enhance Storage Security and Reduce Leakage Risks. Applied Sciences, 13(23), 12830. https://doi.org/10.3390/app132312830

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop