Abstract

Carbon capture and storage has become a practice to reduce the greenhouse effect of carbon dioxide (CO2) on the global climate. Recent studies have generated increasing concerns about CO2 leakage from underground structures. This has called for more research on CH4–CO2 swapping in natural gas hydrate (NGH) reservoirs to lock CO2 in a solid state in underground structures. Because the CH4–CO2 swapping is too slow to be efficient, this study proposes to use geothermal energy to accelerate the process. This paper presents a technical feasibility analysis of using geothermal energy to assist CH4–CO2 swapping for simultaneously storing CO2 in NGH reservoirs and producing the dissociated natural gas. Mathematical models were developed to compute heat transfer from geothermal zones to NGH reservoirs. A case study was carried out using the data from an NGH reservoir in the Shenhu area, Northern South China Sea. The result of the case study indicates that heat conduction dictates the heat transfer process when the heat convection flow rate is less than 0.01 m3/s over a heat-releasing borehole length of 2,000 m. Heat convection can significantly accelerate the heat transfer inside the gas hydrate reservoir. The 15°C (designed gas hydrate dissociation temperature in the studied case) heat front will propagate to the upper and lower boundaries of the gas hydrate reservoir (39 ft or 12 m) in 220 days by heat conduction only. This time can be shortened to 140 days with the aid of a fluid convection rate of 0.005 m3/s. Geothermal heating can significantly increase the initial productivity of wells in heated gas hydrate reservoirs in CO2 swapping processes. When the gas hydrate reservoir is heated from 6 to 16°C, the fold of increase is expected to exceed five in the studied case. This study shows that CH4–CO2 swapping process using geothermal stimulation is a promising method for producing natural gas and locking CO2 permanently in NGH reservoirs. Further studies should first focus on investigations of the effect of CO2-hydrate formation on the CO2 mass transfer inside reservoirs.

1. Introduction

The increasing trend of carbon dioxide (CO2) level in the atmosphere is a great concern to humankind because of its global warming effect. Many carbon capture and storage projects have been initiated in the past 5 years to place CO2 in underground structures such as depleted oil reservoirs. Recent studies indicated a high risk of leakage of CO2 through CO2-exposed oil wells [1]. The leakage can be due to the loss of sealing integrity of wellbore cement resulting from CO2–cement interaction over time [2]. According to the computer simulation by Zhang et al. [3], it should take hundreds of years for CO2 to fully penetrate a 3.5 cm-thick cement sheath through the pore space of cement to reach well casing. But it should take less than a year for CO2 to fully penetrate a 3.5 cm-thick cement sheath through cracks in the cement. Such cracks in cement sheath were evidenced in careful studies such as that by Duguid et al. [1]. Because the sealing capacity deterioration of wellbore cement sheath is not avoidable, it is logical to “blame” the nature of the high mobility of CO2 in its supercritical fluid state in the underground storage reservoirs. It is also logical to store the CO2 in its solid state in low-temperature underground structures such as natural gas hydrate (NGH) reservoirs. This process is called CH4–CO2 swapping in previous studies such as that by Cha et al. [4].

A vast amount (1015–1018 m3 at the standard condition) of natural gas is stored in the NGH in the permafrost and offshore sediments [57]. The amount of natural gas in the NGH deposits is believed to be greater than that in coal, oil, and conventional gas reservoirs combined [8]. The natural gas in the NGH reservoirs is considered the major source of relatively clean energy to support the future development of the world economy [911].

Most of the previous studies on NGH reservoirs have focused on how to produce natural gas from the NGH reservoirs cost-effectively. Three methods have been identified to harvest the natural gas in NGH reservoirs, namely (1) depressurization, (2) thermal stimulation, and (3) the use of hydrate inhibitors. Using the depressurization method, the pressure in the hydrate-bearing zone is lowered below the hydration dissociation pressure at the reservoir temperature [12, 13]. Because NGH is an essential portion of the matrix structure in the reservoir, dissociation of NGH releases water and causes destabilization of the matrix structure. As a result, formation sand production and even wellbore collapse frequently occur during gas production when depressurization is used as a production scheme. With the thermal stimulation method, the temperature of the hydrate-bearing zone is raised above the hydration–dissociation temperature by circulating hot water from the surface [1416]. This method is very costly due to the power consumption in heating the water. Hydrate inhibitors for gas production from NGH reservoirs include salts and alcohol-type chemicals. They can shift the P–T equilibrium by competing with the NGH for guest and host molecules [17, 18]. This method is also costly due to the use of a large number of expensive inhibitors, especially alcohol-type chemicals. A combination of these methods has been tested to facilitate long-term gas production [1921]. None of these methods has been proven economical in offshore operations, although some operations are run at a marginal profit in onshore NGH reservoirs in permafrost zones.

CH4–CO2 swapping is a process where one CO2 molecule replaces one CH4 molecule without destroying the hydrate structure with little water production. The process reduces the possibility of matrix collapse and stratum failure [4]. It is, therefore, more competitive compared with other gas-extraction methods. The process can lock the CO2 in the reservoir permanently with a low probability of leaking into the atmosphere. However, the swapping efficiency for CH4 production is low due to the mass transfer barriers caused by the formation of CO2-hydrates without external energy added to delay the formation of CO2-hydrates [22].

The recent work of Mahmood and Guo [23] provides a new idea for solving this problem. They proposed a geothermal method for producing natural gas from marine gas hydrate reservoirs. The method uses geothermal energy to accelerate gas hydrate decomposition during the depressurization process. Under the reservoir temperature controlled in a certain range, the thermal energy can accelerate the dissociation of the CH4-hydrates, delay the formation of CO2-hydrates, and increase the diffusion coefficient of CO2. The combination of these three effects is expected to create more diffusion channels to promote CO2 penetration. The low-level depressurization allows for gas production with manageable sand production and wellbore collapse.

Mahmood and Guo’s [23] work considered heat transfer due to conduction only and neglected the effect of gas flow on temperature change inside the reservoir. This left a gap between the reliable heat transfer modeling and well productivity modeling. This gap is filled in the current study using a numerical solution of the diffusion–convection equation to better simulate the heat transfer in the CH4–CO2 swapping processes.

2. Heat Transfer Modeling

The heat transfer from a geothermal zone to an NGH reservoir takes two steps: (1) heat transfer from the geothermal zone to a heat-releasing borehole in the NGH reservoir, and (2) heat transfer from the heat-releasing borehole into the NGH reservoir.

Fu et al. [24] proposed a y-shaped wellbore configuration and presented an analytical model for heat transfer from a geothermal zone to an NGH reservoir using two horizontal boreholes. The horizontal borehole at the bottom of the geothermal zone is used to receive heat from the geothermal zone, and the horizontal borehole at the top of the NGH reservoir is used for leasing heat into the NGH reservoir. The heat-receiving borehole is designed horizontally to increase heat-receiving efficiency through adequate fluid contact/retention time for heat transfer. An alternative means of improving the heat-receiving efficiency is to use a deeper vertical borehole if drilling technology (thermal-stable drilling fluids) permits. Fu et al.’s [24] system design was modified in this study to replace the horizontal borehole with a vertical borehole in the geothermal zone, assuming it is drillable with the existing drilling fluids. The modified well structure is illustrated in Figure 1. The CO2 in tank (1) is transferred by compressor (2) to injection well (3), injected through pipe (4) to the bottom of the wellbore in geothermal zone (5), where the CO2 is heated by geothermal energy. The heated CO2 in the heat-receiving hole (6) is guided to the heat-receiving hole (7) placed in the NGH reservoir (8). The CO2 is partially released to the NGH reservoir, and the remaining CO2 returns to the surface through well tubing (9). The amount of CO2 released to the NGH reservoir is controlled by the wellhead choke (10) of the gas-producing well (11). The NGH reservoir is heated by the CO2 stream through heat conduction and heat convection. The released natural gas from the NGH is collected by the gas wellbore (12) and produced through the annulus (13), and processed by the fluid separator (14).

Fu et al.’s [24] study with water as the working fluid demonstrated that the deliverable temperature of the fluid to the heat-receiving borehole is controlled by mainly three factors, including as follows:(i)Pipe and annular insulation.(ii)Fluid flow rate.(iii)Depth of the geothermal zone.

Guo and Zhang [25] presented an analytical model for heat transfer from a horizontal borehole deep into an NGH reservoir by heat conduction only. The model does not consider heat convection by the fluid flow into the NGH reservoir. The heat transfer process in porous media dominated by heat conduction and heat convection is governed by the following convection–diffusion equation [26]:where T is the temperature in °C, x is the distance in m, t is time in s, ϕ is porosity in fraction, u is the convective velocity in m/s, and α is the thermal diffusion coefficient defined by the following:where k is the thermal conductivity of the solid system in W/(m K), is the density of the solid system in kg/m3, and cp is the heat capacity of the solid system in J/(kg K).

The convective velocity is formulated as follows:where is the inlet fluid flow rate (injection rate) in m3/s, is the outlet fluid flow rate (return rate) in m3/s, x is the radial distance in m, and L is the length of the heat-releasing hole in m.

The diffusion–convection Equation (1) with coefficients defined by Equations (2) and (3) describes heat transfer into the reservoir due to the combined thermal diffusion coefficient of the solid and fluid when α is evaluated considering the contribution of both rock and pore fluid, and u takes the velocity of the injected fluid. This equation is believed adequate for simulating the total heat transfer into the gas hydrate reservoir.

Equation (1) can be solved by the finite difference method (FDM) with the following initial and boundary conditions:where T0 is the initial reservoir temperature in °C, rw is the radius of the heat-releasing borehole in m, Tw is the temperature in the heat-releasing borehole in °C, is the special derivative of the temperature outside the external boundary in °C/m, and is the special derivative of the temperature inside the external boundary in °C/m. The profile of the temperature in the heat-releasing borehole (Tw) is given by the analytical model presented by Fu et al. [24]. It is understood that the fluid temperature at any point of the heat-releasing borehole will change with injection time. However, a steady temperature profile is assumed in the whole process. It was noted by Wu et al. [27] that during the flow of CO2 in vertical and horizontal wellbores, the temperature will be affected by the Joule–Thomson effect and friction heat. However, these effects are not considered in Fu et al.’s [24] analytical model and thus are assumed to be negligible in this study.

With an explicit scheme of formulation of FDM, the discretization of Equation (1) takes the form ofwhere n is the ordinal of a time spot, i is the ordinal of a grid block, Δx is the length of a grid block, and Δt is the length of a timestep, i.e.,

The numerical solution of Equation (7) gives the dynamic distribution of temperature inside the NGH reservoir as follows:

3. Gas Well Productivity Modeling

The productivity of the gas production wellbore can be predicted by the analytical model that was originally proposed by Joshi [28]. This model was modified by Guo [29] to include the effects of non-Darcy flow and real gas pseudo-pressure. The model takes the following form:wherewhere qg is the gas production rate in Mscf/day, kH is the horizontal permeability in mD, kV is the vertical permeability in mD, h is the thickness of the NGH reservoir in ft, pe is the reservoir pressure in psi, m(pe) is the real gas pseudo-pressure at pe, pwf is the flowing bottom hole pressures in psi, m(pwf) is the real gas pseudo-pressure at pwf, Ti is the reservoir temperature in °R, L is the length of the horizontal wellbore in ft, s is the Darcy skin factor, D is the non-Darcy coefficient in day/Mscf, and reH is the equivalent radius of the drainage area of the horizontal well in ft. The real gas pseudo-pressure is defined as follows and is computed by the spreadsheet developed by Guo and Ghalambor [30]:where p is pressure in psi, pb is base pressure in psi, is gas viscosity in cp, and z is gas compressibility factor.

Because heating the NGH reservoir will increase real gas pseudo-pressure proportionally, the heating effect should not be considered Equation (13). Applying Equation (13) to both nonheated and heated reservoir conditions gives:where qgNH is the gas production rate of the well in the nonheated NGH reservoir in Mscf/day, peNH is the driving pressure of the nonheated NGH reservoir in psi, pwfNH is the flowing bottom hole pressure in the nonheated NGH reservoir in psi, qgH is the gas production rate of well in the heated NGH reservoir in Mscf/day, peH is the driving pressure of the heated NGH reservoir in psi, and pwfH is the flowing bottom hole pressure in the heated NGH reservoir in psi. The driving pressure is defined as the NGH dissociation pressure at a given reservoir temperature. Dividing Equation (18) by Equation (17) gives the following relation:where FOI is the fold of increase in well productivity due to reservoir heating in the swapping process.

For NGH reservoirs containing no free gas in the initial conditions (Class 1 W), no gas should flow before the reservoir is depressurized to NGH dissociation pressure. Therefore, the peNH and peH are defined as the NGH dissociation pressures at nonheated and heated reservoir conditions. Their values can be determined from the NGH equilibrium curves based on temperature. The pwfNH and pwfH should be selected between the dissociation pressures of CH4-hydrate and CO2-hydrate in the swapping process.

4. Case Study

The Shenhu NGH reservoir in the middle of the North Continental Slope of the South China Sea is composed of clayey silt with permeabilities between 1.5 and 7.4 mD [31]. The NGH saturation ranges from 11.7% to 34% [32]. Table 1 provides basic data relevant to well productivity analysis [33, 34].

Table 2 presents a summary of model parameter values for the well structure shown in Figure 1. The result of Fu et al. [24] interwellbore heat transfer model predicted fluid temperatures of 40°C at the inlet and 28°C at the outlet of the heat-releasing borehole.

Table 3 provides basic data relevant to the heat transfer inside the NGH reservoir. Figure 2 presents Equation (12)-calculated temperature profiles for return flow rates of 100% (no fluid flowing into the reservoir). It indicates that the temperature should increase to 6°C at a distance of 20 m from the wellbore in 6 months. Figure 3 shows model-calculated temperature profiles for a return flow rate of 50% (50% of the injected fluid flows into the reservoir). It implies that the temperature should increase to 6.5°C at a distance of 20 m from the wellbore in 6 months. Figure 4 provides a plot of model-calculated temperature profiles for a return flow rate of 0% (100% of the injected fluid flows into the reservoir). It illustrates that the temperature should increase to 10.5°C at a distance of 20 m from the wellbore in 6 months. Figure 5 shows Equation (12)-calculated effect of fluid flowing into the reservoir on temperature profiles for an injection time of 6 months. It indicates that for a fluid convection flow rate under 0.01 m3/s, heat conduction is the main mechanism that controls heat transfer inside the reservoir.

Figure 6 presents the calculated propagation front of temperature 15°C as a function of fluid circulation time for two fluid convection levels. It shows that the rate of front propagation decreases with time, which is expected for a radial heat-transfer system. The curves imply that the 15°C front will propagate to the upper and lower boundaries of the gas hydrate reservoir (39 ft or 12 m) in 140 days for a 50% return rate of the injected fluid and in 220 days for the pure heat-conduction condition. If this temperature of 15°C is used for dissociating NGH, the heat convection can significantly accelerate heat transfer and thus the CH4–CO2 swapping process.

Figure 7 shows hydrate-phase equilibrium curves for CH4 and CO2. Point A represents the initial reservoir condition (2,053 psi and 279°K). The line BC shows the desired flowing bottom hole pressures and temperatures in geothermal-heated conditions. Coordinates of some points on the gas hydrate dissociation line and the line BC are shown in Table 4. Figure 8 plots the fold of increase given by Equation (19), which indicates that a nonlinear relationship exists between the FOI and heated reservoir temperature. The FOI reaches greater than five at the heated reservoir temperature of about 16°C.

5. Discussion

The accuracy of the result presented in this feasibility analysis is subject to error due to the assumptions made in mathematical modeling. Three issues are discussed here.

First, the mathematical model should overestimate the heat transfer efficiency due to omitted effects: (1) the heat transfer model does not consider the loss of heat to the produced gas, (2) the reservoir temperature drops due to NGH depressurization [13] are not considered, (3) heating of the dissociated gas will cause the gas pressure to increase, which is not considered, (4) Class 1 W hydrate reservoir where there is no free gas in the reservoir in the initial condition, and (5) the formation of CO2 hydrates should reduce CO2 mass transfer, which was not studied in this work.

Second, the well productivity model assumes a Class 1 W hydrate reservoir. This should result in underestimated well productivity by the model if the model is applied to other types of gas hydrate reservoirs where free gas exists in the initial condition. Also, it is understood that the well productivity model assumes that the driving pressure is the hydrate dissociation pressure at the external flow boundary of the dissociated region of the reservoir. Because the boundary distance is time-dependent and controlled by heat transfer efficiency, the well productivity should also be fluid-circulation time-dependent.

6. Conclusions

This paper presents a technical feasibility analysis of using geothermal energy to assist CH4–CO2 swapping for simultaneously storing CO2 in NGH reservoirs and producing the dissociated natural gas. Mathematical models were developed to compute heat transfer from geothermal zones to NGH reservoirs. A case study was carried out using the data from an NGH reservoir in the Shenhu area, Northern South China Sea. The following conclusions are drawn:(i)For the studied gas hydrate reservoir, mathematical modeling of heat conduction and convection shows that heat conduction dictates the heat transfer process when the heat convection flow rate is less than 0.01 m3/s over a heat-releasing borehole length of 2,000 m.(ii)Heat convection accelerates heat transfer inside the gas hydrate reservoir. The 15°C (designed gas hydrate dissociation temperature in the studied case) heat front will propagate to the upper and lower boundaries of the gas hydrate reservoir (39 ft or 12 m) in 220 days by heat conduction only. This time can be shortened to 140 days with the aid of a fluid convection rate of 0.005 m3/s.(iii)Geothermal heating can significantly increase the initial productivity of wells in heated gas hydrate reservoirs in CO2 swapping processes. The rate of a fold of increase increases nonlinearly with heated reservoir temperature. When the gas hydrate reservoir is heated from 6 to 16°C, the fold of increase is expected to exceed five in the studied case.(iv)The accuracy of the result presented in this work is subject to the error due to the assumptions made in mathematical modeling. Future research work should consider heat loss to the produced gas, temperature drops due to NGH depressurization, gas pressure change due to heating, initial free gas in the reservoir, and the effect of formed CO2 hydrates on CO2 mass transfer. The well productivity model should be modified to consider the effects of the initial free gas and the dynamic flow boundary due to the dynamic heat transfer process.

In summary, the CH4–CO2 swapping process using geothermal stimulation is a promising method for producing natural gas and locking CO2 permanently in NGH reservoirs. Further studies should first focus on investigations of the effect of CO2-hydrate formation on the CO2 mass transfer inside reservoirs.

Nomenclature

c:Specific heat of reservoir rock, J/kg °C
Cp:Heat capacity of the fluid inside the wellbore, J/kg °C
dc:Inner diameter of insulation-cement sheath, ft
dp:Inner diameter of pipe, ft
Dc:Outer diameter of insulation-cement sheath, ft
Dp:Outer diameter of pipe, ft
D:Non-Darcy coefficient in day/Mscf
G:Geothermal gradient, °C/ft
h:Reservoir thickness, ft
Iani:Anisotropy index
kH:Horizontal permeability, mD
kv:Vertical permeability, mD
Kc:Thermal conductivity of the insulation-cement, W/(m °C)
Kp:Thermal conductivity of insulation pipe, W/(m °C)
L:Length of the horizontal well, ft
:Fluid flow rate, m3/s
p:Base pressure, psia
pe:Reservoir pressure, psia
pwf:Wellbore flowing pressure, psia
peH:Driving pressure of the heated NGH reservoir, psia
peNH:Driving pressure of the nonheated NGH reservoir, psia
pwfH:Flowing bottom hole pressure in the heated NGH reservoir, psia
pwfNH:Flowing bottom hole pressure in the nonheated NGH reservoir, psia
qg:Well production rate, Mscf/day
qgNH:Gas production rate of well in the nonheated NGH reservoir, Mscf/day
qgH:Gas production rate of well in the heated NGH reservoir, Mscf/day
Qf:Fluid circulation rate, m3/s
Qin:Inlet fluid flow rate, m3/s
Qout:Outlet fluid flow rate, m3/s
rw:Wellbore radius of laterals, ft
reH:Radius of drainage area, ft
s:Darcy skin factor
t:Time, s
tc:Thickness of cement sheath, ft
T:Reservoir temperature, °C
Ti:Initial reservoir temperature, °R
Tf:Temperature of the injected fluid, °C
Tin:Inlet temperature, °C
Tout:Outlet temperature, °C
u:Convective velocity, m/s
x:Radial distance, m
z:Gas compressibility factor
:Thermal diffusion coefficient
ϕ:Porosity
ρ:Density of reservoir rock, kg/m3
γg:Gas-specific gravity
:Gas viscosity, cp.

Data Availability

The data were generated during the study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by BIRD through the project “Safe, sustainable, and resilient development of offshore reservoirs and natural gas upgrading through innovative science and technology: Gulf of Mexico—Mediterranean,” award number EC-19 Fossil Energy.