Next Article in Journal
Feasibility Study of a Heating, Cooling and Domestic Hot Water System Combining a Photovoltaic-Thermal System and a Ground Source Heat Pump
Previous Article in Journal
Coupled Effects of Moisture Content and Inherent Clay Minerals on the Cohesive Strength of Remodelled Coal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation Study on the Effect of Fracturing Technology on the Production Efficiency of Natural Gas Hydrate

1
Engineering College, Jilin University, Changchun 130026, China
2
Key Laboratory of Drilling and Exploitation Technology in Complex Conditions, Ministry of Land and Resources, Changchun 130026, China
*
Authors to whom correspondence should be addressed.
Energies 2017, 10(8), 1241; https://doi.org/10.3390/en10081241
Submission received: 14 July 2017 / Revised: 15 August 2017 / Accepted: 16 August 2017 / Published: 21 August 2017
(This article belongs to the Section L: Energy Sources)

Abstract

:
Natural gas hydrate (NGH) concentrations hold large reserves of relatively pure unconventional natural gases, consisting mainly of methane. Depressurization is emerging as the optimum conversion technology for converting NGH in its reservoir to its constituent water and natural gas. NGH concentrations commonly have a pore fill of over 80%, which means that NGH is a low-permeability reservoir, as NGH has displaced water in terms of porosity. Fracturing technology (fracking) is a technology employed for increasing permeability-dependent production, and has been proven in conventional and tight oil and gas reservoirs. In this work, we carried out numerical simulations to investigate the effects on depressurization efficiency of a variably-fractured NGH reservoir, to make a first order assessment of fracking efficiency. We performed calculations for the variations in original NGH saturation, pressure distribution, CH4 gas production rate, and cumulative production under different fracturing conditions. Our results show that the rate of the pressure drop within the NGH-saturated host strata increases with increased fracturing. The CH4 gas production rate and cumulative production are greatly improved with fracturing. Crack quantity and spacing per volume have a significant effect on the improvement of NGH conversion efficiencies. Possibly most important, we identified an optimum fracking value beyond which further fracking is not required.

1. Introduction

Natural gas hydrates (NGH) are white, solid, crystalline compounds, composed of light hydrocarbons, CO2, H2S, and water, under certain temperature and pressure conditions. Among them, methane hydrate is the most widely-distributed in nature [1,2,3,4,5]. Marine gas hydrate is mainly developed in hydrate stability zones (HSZ) below a 500-m water depth in the open ocean, where temperatures and pressures are suitable for its spontaneous formation. About 90% of the world's oceans present temperature and pressure conditions that are suitable for the formation of HSZ, which provide a good environmental basis for NGH development [6,7]. When there is enough methane, and water enters the pores of undercompacted seabed sediments, it is possible to form NGH [7].
NGH has a low environmental risk, and exhibits the efficient sequestration of natural gas from both biogenic and thermogenic sources. Very clean natural gas can be produced from NGH concentrations, especially from sandy turbidites, from which it is already known in the industry how to produce conventional hydrocarbons. It is estimated that reserves are approximately 2 × 1016 m3, equivalent to double the world's proven, conventional, total carbon. World estimates for gas from NGH reserves of NGH in sands are >40,000 Tcf (1 Tcf = 1 × 1013 ft3 = 283.17 × 109 m3) [8,9,10,11,12,13]. NGH should be converted in situ to its constituent gas and water. A number of conversion methods exist, but early production testing and modeling indicate that depressurization will be the ideal method to use [14,15,16,17,18].
In 2007, samples from the SH2, SH3, and SH7 sites were successfully drilled from the Shenhu area in the northern South China Sea (Figure 1) [19,20,21]. Taking the SH7 site as an example, the hydrate layers were 18–34 m thick, and the sediment porosity, permeability, and NGH pore fill were 33–48%, 7.5 × 10−14 m2, and 20–44%, respectively [19,21].
In 2013, a large number of gas hydrate samples were obtained from the northeastern continental slope of the South China Sea. The NGH host sediments in this area comprised sandy layers of turbidite sediments. Natural gas hydrate saturation of porosity ranged from 45% to 100%; the amount of gas seemed to be very large. The NGH was massive, layered, tumor-like veins, dispersed in all combinations of silty clay and clastic limestone [22]. In 2013, Nippon and Mie Prefecture NGH were produced using depressurization. The trial production achieved good results [23,24,25,26]; however, commercial exploitation lacks depth of study [23,26]. In recent years, TOUGH+HYDRATE v1.0 models have been widely used in gas hydrate simulations. Taking the Shenhu area, South China Sea, as an example, Li [19] used this model to evaluate gas production potential by depressurization and thermal stimulation from the SH7 site. Hu [27] and Jin [28] used this model to analyze sensitive parameters (porosity, permeability, NGH pore fill, bottom-hole pressure, and thermal stimulation intensity) via depressurization and thermal stimulation from the SH2 site. Su [29] used this model to evaluate the effects of thermal stimulations on gas production from the SH2 site. Zhang [30] used this model to preliminarily estimate gas production potential by depressurization and thermal stimulation from the SH2 site. The results showed that cumulative CH4 volume in the Shenhu area was not large, and NGH exploitation was seriously inhibited, mainly because of the low porosity and poor permeability.
The primary repository of NGH is in sandy and silty marine turbidites. Clay and clastic limestone and secondary permeability related to structure also host NGH in some areas of the South China Sea, but their viability as producible gas resources is unlikely in their natural state. The NGH-enriched stratum encountered in the well, which we model in this paper, has a low porosity and poor permeability, which hinder the transfer of heat and pressure between the production wells and the stratum, and reduces the production efficiency of the decomposition of gas and the continuous dissociation of NGH. The problem of how to efficiently and economically exploit natural gas hydrate in low-permeability marine sediment reservoirs is an important issue. Being able to produce NGH from marine sediments, in which the grain size distribution hosting NGH is less producible than well-sorted sand, will greatly enlarge the potential use of NGH as a gas resource.
In 1947, the world's first successful fracturing well was built in the state of Kansas, the United States [31,32]. Since the 1970s, the technology has expanded to include low-permeability oil and gas shale, as well as tight sand reservoirs [33]. Multi-stage fracturing was perfected and has been used to produce gas and oil from shale reservoirs since about 2003 in the United States [34]. Fracturing technology can be divided into vertical and horizontal fracturing (Figure 2), according to the different characteristics of the formations [35].
Increasing the permeability of the hydrate layer can improve heat and mass transfer efficiency, increase the gas migration channels in the hydrate layer, accelerate the hydrate dissociation and discharge, and increase the gas production rate and the cumulative production of CH4 [4,13,28,30,36].
The application of fracturing technology to the exploitation of gas hydrate can effectively increase reservoir permeability and enhance the efficiency of gas hydrate production. The combination of fracturing technology and existing methods to exploit natural gas hydrates would solve some of the problems of existing methods [17,27,37,38]. Fracturing technology is the primary mechanical solution for increasing porosity, interconnecting artificial cracks and natural cracks, and forming a mutual interconnected crack network in NGH concentrations, in order to enhance the conductivity of NGH concentrations. Depressurization and fracturing combine to exploit NGH [27,37,38], which increases the rate of the pressure drop within the NGH-saturated host strata and promotes NGH decomposition in cracked areas, in addition to being conducive to the discharge of methane. Thermal stimulation and fracturing combine to exploit NGH [27,37], as hot fluid enters NGH concentrations through the crack network and directly heats NGH to increase the mining radius and reduce heat loss. For CO2 replacement [17], CO2 enters NGH concentrations through the crack networks. The large contact area between CO2 and NGH increases replacement efficiency by CO2.
Of course, fracturing is also inadequate in NGH exploitation. The fracturing effect is not easy to control, sometimes resulting in a failure in fracturing. In addition, there is a certain risk that fracturing can cause formation instability [39,40].
Oceanic NGH mining may exacerbate global greenhouse effects and worsen the marine ecological environment, resulting in a series of environmental effects, and also probably lead to submarine landslides, seabed collapse, and other geological disasters [39,40]. There is a certain amount of risk in the application of fracturing technology to oceanic NGH mining, as it will further reduce the stability of submarine formations. Taking the trial production in China and Japan as an example [23,25], in order to avoid risks and ensure safety, researchers monitored seabed deformation, reservoir stability, and in situ methane leakage during the trial production. The results proved that in situ monitoring is an effective way to reduce the risks of NGH mining. Similarly, if a reasonable implementation and monitoring plan is adopted, the risk of fracturing oceanic NGH can be reduced.
In this paper, TOUGH+HYDRATE numerical simulation software was used to simulate and compare NGH reservoirs before and after fracturing. The influence of the fracturing parameters on gas hydrate production efficiency was studied.

2. Simulation Preparation

2.1. Numerical Model and Simulation Parameters

2.1.1. Numerical Simulation Code

The simulator model used in this work was TOUGH+HYDRATE v1.0 by Moridis from the Lawrence Berkeley National Laboratory (Berkeley, CA, USA). The model can simulate the formation and decomposition of natural gas hydrate, phase equilibrium, seepage, and heat and mass transfer processes under complex conditions and non-isothermal conditions. In addition, the model can simulate production from natural CH4-hydrate deposits in the subsurface (i.e., in permafrost and deep ocean sediments) as well as laboratory experiments of hydrate dissociation/formation in porous/fractured media [41]. TOUGH+HYDRATE v1.0 can also simulate the formation and decomposition of gas hydrate under equilibrium and kinetics. The model includes four phases (liquid, gas, hydrate, and ice), and four components (water, methane, hydrate, and water-soluble inhibitors, such as salt, alcohol, etc.). By constructing fractured-porous media, the gas hydrate under different geological conditions in permafrost regions and from deep-sea reservoirs can be simulated using the methods of depressurization, heating injection, and injection inhibition. TOUGH+HYDRATE v1.0 is a numerical simulator developed by Moridis at the Lawrence Berkeley National Laboratory, USA. It is the first member of TOUGH+, and the successor to TOUGH2.

2.1.2. System Parameters and Initialization of the Model

The geologic system used in this study was according to the SH7 site in the Shenhu area, South China Sea. The hydrate samples from the SH7 site were almost pure methane hydrate (99.2%) in NGH concentrations [19,21]. Therefore, only methane hydrate was simulated. The system parameters and the initial conditions of the simulation are shown in Table 1. The main parameters in the simulation were derived from previous literature on gas hydrate reservoirs in the region [19,21,27].
The water depth of the simulated gas hydrate reservoir was 1108 m, and the hydrate-bearing layer (HBL) was located in the area of 155–177 m below the seabed, with a thickness of 22 m. The gas hydrate reserve temperature was T0 = 14.15 °C, the pressure was P0 = 13.83 MPa, the saturation was SH = 44%, the porosity was 38%, and the permeability was 7.5 × 10−14 m2 (75 mD).

2.2. Design of the Production Well and Fracturing Crack

2.2.1. Production Well Design

In this paper, gas-hydrate reservoirs were exploited using a vertical well, and we investigated the depressurization efficiency effects of a variably-fractured NGH reservoir to make a first order assessment of fracking efficiency. The simulation system was cylindrical. As shown in Figure 3, the production well was designed according to Su [29] and Li [42]. The well was located at the center of the cylindrical simulation system, with a well radius of rw = 0.1 m. The production interval was located in the middle of gas hydrate reservoir, at a height of 6 m. The production interval was set in the middle of the hydrate layer; because of pores and cracks in the upper and lower layers, natural gas is likely to overflow from the cover layer. The upper and lower hydrate layers temporarily seal and reduce the escape of natural gas through the caprock to a certain extent.

2.2.2. Fracturing Crack Design

The formation of NGH in host sediment pore space results in a higher bulk modulus and increased mechanical strength. In an ideal case, where a hydrate deposit had a sufficient brittle response to hydraulic fracturing, our model would provide a base case with which actual testing can be compared in order to assess the likelihood of artificial fracturing of inducing additional permeability in semi-consolidated marine sediments, which, without NGH, would be expected to respond in a more mechanically-ductile manner.
This paper addresses vertical well fracturing cracks. When fracturing a target formation, the key parameters of cracks are not only affected by the stress distribution of the formation, but are also closely related to the physical and mechanical properties of the rock and jet parameters [43]. Therefore, in this paper, a reasonable simplification of crack fracturing was performed:
  • To consider the influence of the main fracture, secondary cracks after fracturing were ignored;
  • Cracks only formed in the horizontal direction;
  • There was only water in the cracks, without hydrate or broken rock particles.
As shown in Figure 4, the crack length was Lf = 40 m. The cracks were divided into spacing categories, Δl = 1 m, 2 m, 3 m, and 5 m, and crack height, h = 10 mm, respectively, for one, three, and five cracks.
The permeability of the porous medium varied according to the porosity. The porosity and permeability have the following relationships [44,45,46]:
k k 0 = F S = ( 0 ) n
k k 0 = F S = ( c 0 c ) n
where k0 is the formation permeability, k is the formation permeability after the porosity change, 0 is the formation porosity, is the porosity of the formation after the change, and c is a non-zero critical porosity. In Equation (1), n is 2 or 3; in Equation (2), n is 10 or more. After calculations, the permeability of the cracks was 5.2 × 10−13 m2 when h = 10 mm. Table 2 shows the cracks parameters.

3. Simulation Results and Discussion

As shown in Table 2, we carried out numerical simulations using two group parameters to analyze the effects on the depressurization efficiency of crack quantity and other parameters. A crack quantity of 0 referred to the hydrate layer before fracturing. When the crack quantity was equal to 1, there was one crack in the hydrate layer. When the crack quantity was equal to 3, there were three cracks in the hydrate layers. When the crack quantity was 1, the crack was located in the center of the production interval. When the crack quantities were 3 and 5, the middle crack was located in the center of the production interval. The distributions of the remaining cracks are shown in Table 3.

3.1. Effect of Crack Quantity on the Efficiency of Natural Gas Hydrate Exploited by Depressurization

To study the effect of crack quantity on depressurization efficiency, four kinds of crack parameters were used in the first group. The spatial distributions of SH, P, and CH4 gas production rates and cumulative volumes of the natural gas hydrate layer were compared under different fracturing conditions.

3.1.1. Spatial Distributions of SH

Figure 5 shows the distribution of SH after one month and one year of depressurization under four conditions. Figure 5A shows the SH distribution after one month of exploitation. As the crack quantity increased, NGH dissociation area increased and NGH dissociated faster in the fracturing area. In the early stage of dissociation, the fracturing cracks (cracks) promoted NGH dissociation; Figure 5B represents SH distribution after one year of exploitation. According to Figure 5B, the disadvantages of depressurization are obvious; NGH dissociation was seriously inhibited, mainly due to insufficient energy supply, low heat efficiency and mass transfer property, and gas migration channels. Figure 5B 1-1 and 1-2 show that secondary natural gas hydrates appear at the decomposition front. Because of the large amount of CH4 in the production wells, as well as the endothermic reaction of hydrate decomposition, the decomposition front increases in pressure, decreases in temperature, and forms secondary natural gas hydrates. Figure 5B demonstrates that, when the NGH layer contains three or five cracks, the hydrate decomposition rate was promoted. In comparison with Figure 5B 1-3 and 1-4, the rate of NGH dissociation increased with increasing the crack quantity and decreasing the crack spacings. Therefore, the spacings of the cracks have a certain effect on hydrate efficiency.

3.1.2. Spatial Distributions of P

Figure 6 shows the distribution of P after one month and one year by depressurization under four conditions. When the temperature was 14 °C, the equilibrium pressure of methane hydrate was approximately 11 MPa [47]. When the pressure was lower than the equilibrium pressure, the hydrate began to decompose. Figure 6A shows that the rate of the pressure drop within the NGH-saturated host strata increased upon increasing the fracturing cracks in the initial stage. Figure 6B shows that the rate of the pressure drop within the NGH-saturated host strata decreased gradually. Moreover, the rate of the pressure drop within the NGH-saturated host strata, with fracturing cracks, was higher than the original NGH strata. Comparing the radius of the pressure drop under the four conditions, the rate of the pressure drop and the CH4 collection within the NGH-saturated host strata increased with fracturing cracks, and the cracks avoided the formation of secondary hydrates at the decomposition front. In the undecomposed region of natural gas hydrate, the transfer rate of the pressure drop first increased in the presence of fracturing cracks, and then decreased with increasing crack quantities.

3.1.3. Production Rate and Cumulative Volume of CH4

Figure 7 shows the CH4 production rate and the cumulative volume curve after one year of depressurization under the four conditions. The CH4 production rate first increased and then decreased, reaching its maximum in 40 days. At that time, the CH4 gas production rates of 1-3 and 1-4 were significantly greater than that of 1-1. This is because, in the initial stage of hydrate dissolution, the porosity and permeability of the NGH-saturated host strata with fracturing cracks was higher than that of the original NGH strata, so the NGH could be decomposed faster. Thus, CH4 was rapidly discharged.
In the early stages of exploitation, the CH4 gas production rate of 1-2 was higher than that of 1-1. NGH dissociation was faster from a single fracture, and the CH4 migration rate was small. At this time, one crack did not promote hydrate decomposition. Comparing the CH4 gas production rates of 1-3, 1-4, and 1-1 showed that the rate of CH4 production increased with increasing crack quantities. As the crack quantities increased, the porosity and permeability of the fractured area increased, which enlarged the diameter of the gas migration channels for CH4 migration into the production wells, and promoted NGH decomposition. However, compared with the CH4 production rates of 1-3 and 1-4, a crack quantity of greater than three weakened the hydrate dissociation efficiency rate.
As exploitation progressed (after 190 days), the 1-2 fracturing effect became apparent and the CH4 gas production rate increased and exceeded that of 1-1. This finding was observed because fracturing can promote a pressure drop in the internal hydrate layer and promote the decomposition of the hydrate as exploitation progressed. At the same time, one crack can effectively transport the decomposed CH4 into the production well because of the lower rate of natural gas hydrate dissociation, thereby reducing the pressure of the decomposition front and increasing the CH4 production rate. After 300 days, the CH4 production rate tended to stabilize, and the effects on the depressurization efficiency of the fractured NGH reservoir was obvious. At the same time, the CH4 production rates increased by 10.48%, 19.85%, and 25.62%, respectively, and the cumulative volume of CH4 increased by 0.67%, 19.07%, and 25.87%, respectively, in comparison with 1-2, 1-3, 1-4, and 1-1. Increasing the number of cracks greatly improved the natural gas hydrate decomposition efficiency by depressurization. Comparing the three quantities of cracks, the 1-3 CH4 production rate increased by 8.51%, and the cumulative volume of CH4 increased by 18.27% over 1-2. When comparing 1-4 with 1-3, the CH4 production rate increased by 4.81% and the cumulative volume of CH4 increased by 5.71%. Overall, as the number of fractured cracks increased from one to three to five, the efficiency of NGH dissociation first increased and then decreased.

3.2. Effect of Crack Spacing on the Efficiency of Natural Gas Hydrate Exploited by Depressurization

From the distribution of SH in Figure 5, the cracks promoted NGH dissociation by depressurization, and the dissociation rates varied according to the crack spacing conditions. Therefore, when the crack spacing was the appropriate value, the rate of NGH dissociation was the highest in that fracturing region.

3.2.1. CH4 Production Rate and Cumulative Volume

Figure 8 shows the CH4 production rate and cumulative volume after five years of depressurization exploitation under four crack spacing conditions. The CH4 production rate did not have a large disparity at the initial stage of exploitation, and CH4 gas production rates gradually decreased with the progression of exploitation. The CH4 production rate of the hydrate layer, Δl =3 m, was considerably greater than those of the other three conditions at the same time. By observing the cumulative volume curve of CH4, under the conditions of Δl = 1 m, Δl = 2 m, and Δl = 5 m, the cumulative volume of CH4 was similar. However, the cumulative volume of CH4 was higher under the condition of Δl = 3 m, and the cumulative volume of CH4 was increased by 21.14%, 26.19%, and 17.78% for the conditions of Δl = 1 m, Δl = 2 m and Δl = 5 m, respectively. Comparing Δl = 3 m and the original NGH layer in Figure 7, the cumulative volume of CH4 increased by 43.49% in one year. The results show that the fracturing effect has the most obvious influence on the NGH production and the production efficiency is the highest when the crack spacing is Δl = 3 m. To determine why the CH4 production rates and cumulative volumes were considerably greater than those of the other three conditions, the SH distributions of these four conditions were analyzed.

3.2.2. Spatial Distribution of SH

Figure 9 shows the distribution of SH after one year and five years of depressurization under four crack spacing conditions. One year after the exploitation, the smaller crack spacing of 2-1 led to the complete dissociation of NGH around the crack areas; NGH dissociation in a non-fracturing area was slower. With 2-2, we can see that the crack spacing was moderate and NGH around the fracturing area dissociated faster. The distances among the three cracks in 2-4 were larger, and the upper and lower cracks were closer to the upper and lower boundaries of the NGH layer, which promoted the dissociation of the upper and lower NGH. Because of the large distance between adjacent cracks, the NGH dissociation between cracks had no obvious effect. In 2-3, the NGH dissociation among the cracks was promoted by the interactions between two adjacent cracks. Since the upper and lower cracks in 2-3 are closer to the upper and lower boundaries of the hydrate layer than those in 2-2, the cracks in this area appear larger and their role in NGH decomposition is more obvious. It can be seen from Figure 9B that, after five years of exploitation, the cracks in 2-1 and 2-2 are more concentrated, and the rate of NGH dissociation is greater in the fracturing area. NGH in the fracturing area was limited by the heat compensation mechanism; the temperature of the cracking area gradually reduced and formed secondary hydrates, which inhibited NGH dissociation to a certain extent. When Δl = 3 m, the three cracks were distributed in the upper, middle, and lower regions of the NGH layer, respectively. The cracks not only promoted NGH dissociation around fracturing regions, but also promoted the dissociation of the upper and lower NGH layers. The NGH dissociation range could increase and the concentration of the NGH dissociation region could be avoided, which would prevent a local temperature decrease and a pressure increase, which reduces the inhibition of the hydrate dissociation via the thermal compensation mechanism. Meanwhile, the cracks in 2-4 were far apart and the NGH dissociation rate was small. This also revealed that the CH4 production rate and cumulative volume were much higher than those of the other three cases, where the crack spacing was Δl = 3 m.

4. Conclusions

In this paper, an NGH reservoir was simulated using TOUGH+HYDRATE v1.0 via the depressurization method, and the effects on production efficiency of a variably-fractured NGH reservoir were compared. Analyses on the change in distribution of SH, P, as well as the CH4 production rate and the cumulative volume under different crack quantities, crack heights, and crack spacings are as follows:
  • The exploitation of a fractured NGH reservoir using the depressurization method increased the transfer rate of the pressure drop to the interior NGH layer, and adjacent cracks promoted hydrate decomposition.
  • The exploitation of a fractured NGH reservoir using the depressurization method increased CH4 production rates by a maximum of 25.62%, and CH4 production rates increased with increasing crack quantity, although the growth range was reduced.
  • The exploitation of NGH reservoirs with different crack spacings (Δl = 1 m, 2 m, 3 m and 5 m) using the depressurization method initially increased, and then decreased the CH4 production rate and cumulative volume, which indicated that Δl = 3 m was the most favorable crack spacing for NGH exploitation, and the cumulative volume of CH4 increased by 43.49%.
In addition, the increase in productivity in points 2 and 3 was suggested by the numerical simulation, but has not been validated in the field.
In conclusion, the fracturing effect can effectively improve the exploitation efficiency of natural gas hydrate reservoirs. When the crack quantity was three and the crack spacing was Δl = 3 m, the natural gas hydrate exploitation efficiency was the most obvious.

Acknowledgments

This study has been supported by National Natural Science Foundation of China (Grant No.41672361, Grant No. 41502343 and Grant No. 51474112).

Author Contributions

Lin Yang analyzed the data and wrote the manuscript. Lin Yang, Yong Chen and Xitong Li performed the simulations. Rui Jia, Lin Yang and Chen Chen conceived the original ideas. All authors discussed the results and commented on the manuscript. Chen Chen, Youhong Sun and Wei Guo directed the overall project.

Conflicts of Interest

The authors declare no conflicts of interest.

Notation

zposition of HBL below ocean surface(m)
ZHHBL thickness(m)
H1Depth of HBL below seafloor(m)
H2Depth of seafloor(m)
GThermal gradient below seafloor(°C/m)
P0Initial pressure (at base of HBL)(MPa)
ΔPProduction pressure(MPa)
PcapCapillary pressure(MPa)
P01Atmosphere pressure(Pa)
T0Initial temperature (at base of HBL)(°C)
kx, ky, kzIntrinsic permeability(m2)
kcPermeability of fracturing cracks (h2 = 10 mm)(m2)
krAAqueous relative permeability(m2)
krGGas relative permeability(m2)
KdryDry thermal conductivity(W/(kg·°C))
KwetWet thermal conductivity(W/(kg·°C))
KΘThermal conductivity(W/(kg·°C))
ΦPorosity
ρRGrain density(kg/m3)
SHSaturation of natural gas hydrate
SASaturation of aqueous
rRadius(m)
XsSalinity
λVan Genuchten exponent—Table 1
hCrack height(mm)
LfCrack length(m)
ΔlCrack spacing(m)

Subscripts and Superscripts

AAqueous phase
BBase of HBL
capCapillary
GGas phase
HBLHydrate-bear layer
irAIrreducible aqueous phase
irGIrreducible gas
nPermeability reduction exponent—Table 1
nGGas permeability reduction exponent—Table 1
OBOverburden
UBUnderburden

References

  1. Sloan, E.D.; Koh, C.A. Clathrate Hydrates of Natural Gases, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  2. Sun, Y.H.; Jia, R.; Guo, W.; Zhang, Y.Q.; Zhu, Y.H.; Li, B.; Li, K. Design and experiment study of the steam mining system for natural gas hydrates. Energy Fuels 2012, 26, 7280–7287. [Google Scholar] [CrossRef]
  3. Su, Z.; Cao, Y.C.; Yang, R.; Zhang, K.N.; Wu, N.Y. Feasibility of gas production from hydrate reservoir considering heat conduction: Taking Shenhu area in the South China Sea as an example. Geoscience 2011, 25, 608–616. [Google Scholar]
  4. Qorbani, K.; Kvamme, B.; Kuznetsova, T. Using a reactive transport simulator to simulate CH4 production from bear island basin in the Barents Sea utilizing the depressurization method. Energies 2017, 10, 187. [Google Scholar] [CrossRef]
  5. Katagiri, J.; Yoneda, J.; Tenma, N. Multiobjective optimization of the partical aspect ratio for gravel pack in a methane-hydrate reservoir using pore scale simulation. J. Nat. Gas Sci. Eng. 2016, 35, 920–927. [Google Scholar] [CrossRef]
  6. Max, M.D.; Johnson, A.H. Exploration and Production of Oceanic Natural Gas Hydrate: Critical Factors for Commercialization; Springer International Publishing AG: Basel, Switzerland, 2016; p. 405. [Google Scholar]
  7. Sun, X.J.; Chen, Y.F.; Li, L.D.; Cui, Q.; Li, Q.P. Triaxial compression test on synthetic core sample with simulated hydrate-bearing sediments. Pet. Drill. Tech. 2012, 40, 52–57. [Google Scholar]
  8. Max, M.D.; Johnson, A.H. Diagenetic methane hydrate formation in permafrost: A new gas play? In Proceedings of the OTC Arctic Technology Conference, Houston, TX, USA, 7–9 February 2011. [Google Scholar]
  9. Johnson, A.H. Gas Hydrate In: GEA, 2011: The Global Energy Assessment. IIASA, Laxenburg; Austria and Cambridge University Press: Cambridge, UK; New York, NY, USA, 2012; pp. 35–43. [Google Scholar]
  10. Mas, M.D.; Johnson, A.H. Could gas hydrate in fine-grained sediments be a precursor for some shale gas deposits? Pet. Geosci. 2012, 18, 231–238. [Google Scholar]
  11. Max, M.D.; Johnson, A.H. Methane Hydrate/Clathrate Conversion. In Clean Hydrocarbon Fuel Conversion Technology; Khan, M.R., Ed.; Woodhead Publishing Series in Energy No. 19; Woodhead Publishing Ltd.: Cambridge, UK, 2011; pp. 413–434. ISBN 1 84569 727 8. [Google Scholar]
  12. Liang, Y.P.; Li, X.S.; Li, B. Assessment of Gas Production Potential from Hydrate Reservoir in Qilian Mountain Permafrost Using Five-Spot Horizontal Well System. Energies 2015, 8, 10796–10817. [Google Scholar] [CrossRef]
  13. Ruan, X.; Li, X.S.; Xu, C.G. Numerical investigation of the production behavior of methane hydrate under depressurization conditions combined with well-wall heating. Energies 2017, 10, 161. [Google Scholar] [CrossRef]
  14. Wood, D.A. Gas hydrate research advances steadily on multiple fronts: A collection of published research (2009–2015). J. Nat. Gas Sci. Eng. 2015, 24, A1–A8. [Google Scholar] [CrossRef]
  15. Wang, B.; Huo, P.; Luo, T.; Fan, Z.; Liu, F.L.; Xiao, B.; Yang, M.J.; Zhao, J.F.; Song, Y.C. Analysis of the Physical Properties of Hydrate Sediments Recovered from the Pearl River Mouth Basin in the South China Sea: Preliminary Investigation for Gas Hydrate Exploitation. Energies 2017, 10, 531. [Google Scholar] [CrossRef]
  16. Zhao, J.; Zhang, L.; Chen, X.; Zhang, Y.; Liu, Y.; Song, Y. Combined replacement and depressurization methane hydrate recovery method. Energy Explor. Exploit. 2016, 34, 129–139. [Google Scholar] [CrossRef]
  17. Birkedal, K.A.; Hauge, L.P.; Graue, A.; Ersland, G. Transport Mechanisms for CO2-CH4 Exchange and Safe CO2 Storage in Hydrate-Bearing Sandstone. Energies 2015, 8, 4073–4095. [Google Scholar] [CrossRef]
  18. Voronov, V.P.; Gorodetskii, E.E.; Muratov, A.R. Study of methane replacement in hydrates by carbon dioxide in a cyclic process. J. Nat. Gas Sci. Eng. 2014, 21, 1107–1112. [Google Scholar] [CrossRef]
  19. Li, G.; Li, X.S.; Zhang, K.N.; Li, B.; Zhang, Y. Effects of impermeable boundaries on gas production from hydrate accumulations in the Shenhu area of the south China sea. Energies 2013, 6, 4078–4096. [Google Scholar] [CrossRef]
  20. Sun, Z.X.; Xin, Y.; Sun, Q.; Ma, R.L.; Zhang, J.G.; Lv, S.H.; Cai, M.Y.; Wang, H.X. Numerical simulation of the depressurization process of a natural gas hydrate reservoir: An attempt at optimization of field operational factors with multiple wells in a real 3D geological model. Energies 2016, 9, 714. [Google Scholar] [CrossRef]
  21. Li, G.; Moridis, G.J.; Zhang, K.; Li, X.S. Evaluation of gas production potential from marine gas hydrate deposits in Shenhu area of South China Sea. Energy Fuels 2010, 24, 6018–6033. [Google Scholar] [CrossRef]
  22. Liang, J.Q.; Wang, H.B.; Su, X.; Fu, S.Y.; Wang, L.F.; Guo, Y.Q.; Chen, F.; Shang, J.J. Natural gas hydrate formation conditions and the associated controlling factors in the northern slope of the South China Sea. Nat. Gas Ind. 2014, 34, 128–135. [Google Scholar]
  23. Yamamoto, K. Methane hydrate offshore production test in the eastern Nankai trough; a milestone on the path to real energy resource. In Proceedings of the 8th International Conference on Gas Hydrates (ICGH8), Beijing, China, 28 July–1 August 2014. [Google Scholar]
  24. Kida, M.; Suzuki, K.; Kawamura, T.; Oyama, H.; Nagao, J. Characteristics of natural gas hydrates occurring in pore-spaces of marine sediments collected from the Eastern Nankai trough, off Japan. Energy Fuels 2009, 23, 5580–5586. [Google Scholar]
  25. Kanno, T.; Takekoshi, M.; Wang, X.; Chee, S.; Fukuhara, M.; Osawa, O.; Yamamoto, K.; Fujii, T.; Takayama, T.; Suzuki, K. Temperature measurement of gas hydrate dissociation during the world-first offshore production test. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 5–8 May 2014. [Google Scholar]
  26. Yamamoto, K. Overview and introduction: Pressure core-sampling and analyses in the 2012–2013 MH21 offshore test of gas production from methane hydrates in the eastern Nankai Trough. Mar. Pet. Geol. 2015, 66, 296–309. [Google Scholar] [CrossRef]
  27. Hu, L.T.; Zhang, K.N.; Gao, T. Numerical studies of gas hydrate production from gas hydrate zong using heat injection and depressurization in Shenhu area, the South China Sea. Geoscience 2011, 25, 675–682. [Google Scholar]
  28. Jin, G.R.; Xu, T.F.; Liu, X.; Xin, X. Optimization of gas production from hydrate deposit using joint depressurization and thermal stimulation. J. Central South Univ. (Sci. Technol.) 2015, 46, 1534–1543. [Google Scholar]
  29. Su, Z.; Li, H.; Wu, N.Y.; Yang, S.X. Effect of thermal stimulation on gas production from hydrate deposits in Shenhu area of the South China Sea. Earth Sci. 2013, 56, 601–610. [Google Scholar] [CrossRef]
  30. Zhang, K.N.; Moridis, G.J.; Wu, N.Y.; Li, X.; Reagan, M.T. Evaluation of alternative horizontal well designs for gas production from hydrate deposits in the Shenhu Area, South China Sea. In Proceedings of the International Oil and Gas Conference and Exhibition in China, Beijing, China, 8–10 June 2010. [Google Scholar]
  31. Wu, Q.; Xu, Y.; Liu, Y.Z.; Ding, Y.H.; Wang, X.Q.; Wang, T.F. The current situation of stimulated reservoir volume for shale in U.S. and its inspiration to China. Oil Drill. Product. Technol. 2011, 33, 1–7. [Google Scholar]
  32. Lancaster, D.E.; Holditch, S.A.; Mcketta, S.F.; Hill, R.E.; Guidry, F.K.; Jochen, J.E. Resercoir evaluation, completion techniques, and recent reesuls from Barnett shale development in the fort worth basin. In Proceedings of the SPE Annual Technical Conference and Exhibition, Washington, DC, USA, 4–7 October 1992. [Google Scholar]
  33. Li, Y.; Cao, G. Development technology for low-permeability sandstone reservoirs in Shengli Oil field. Pet. Explor. Dev. 2005, 32, 123–126. [Google Scholar]
  34. Wu, Q.; Xu, Y.; Wang, T.F.; Wang, X.Q. The revolution of reservoir stimulation: An introduction of volume fracturing. Nat. Gas Ind. 2011, 31, 7–12. [Google Scholar]
  35. Zhang, Y.J.; Li, Z.W.; Guo, L.L.; Gao, P.; Jin, X.P.; Xu, T.F. Electricity generation from enhance geothermal systems by oilfield produced water circulating through reservoir stimulated by staged fracturing technology for horizontal wells: A case study in Xujiaweizi area in Daqing Oilfield, China. Energy 2014, 78, 788–805. [Google Scholar] [CrossRef]
  36. Sun, Y.H.; Li, B.; Wei, G. Comparative Analysis of a Production Trial and Numerical Simulations of Gas Production from Multilayer Hydrate Deposits in the Qilian Mountain Permafrost. J. Nat. Gas Sci. Eng. 2014, 21, 456–466. [Google Scholar] [CrossRef]
  37. Song, Y.; Zhang, L.; Lv, Q.; Yang, M.; Ling, Z.; Zhao, J. Assessment of gas production from natural gas hydrate using depressurization, thermal stimulation and combined method. RSC Adv. 2016, 6, 47357–47367. [Google Scholar] [CrossRef]
  38. Qorbani, K.; Kvamme, B. Non-equilibrium simulation of CH4 production from gas hydrate reservoirs through the depressurization method. J. Nat. Gas Sci. Eng. 2016, 35, 1544–1554. [Google Scholar] [CrossRef]
  39. Kim, A.R.; Cho, G.C.; Lee, J.Y.; Kim, S. Numerical simulation on geomechanical stability during gas hydrate production by depressurization. In Proceedings of the 11th International Symposium on Cold Regions Development Studies (IACORDS), Incheon, Korea, 18–20 May 2016. [Google Scholar]
  40. Goto, S.; Matsubayashi, O.; Nagakubo, S. Simulation of gas-hydrate dissociation caused by repeated tectonic uplift events. J. Geophys. Res. Solid Earth 2016, 5, 3200–3219. [Google Scholar] [CrossRef]
  41. Moridis, G.J.; Kowalsky, M.B.; Pruess, K. TOUGH+HYDRATE v1.0 User’s Manual: A Code for the Simulation of System Behavior in Hydrate-Bearing Geologic Media; Lawrence Berkeley National Laboratory: Berkeley, CA, USA, 2008.
  42. Li, G.; Li, X.S.; Chen, Q.; Chen, Z.Y. Numerical simulation of gas production from gas hydrate zone in Shenhu are, South China Sea. Acta Chim. Sin. 2010, 68, 1083–1092. [Google Scholar]
  43. Zhao, Y.Z.; Qu, L.Z.; Wang, X.Z.; Cheng, Y.F.; Shen, H.C. Simulation experiment on prolongation law of hydraulic fracture for different lithologic formations. J. China Univ. Pet. 2007, 31, 63–66. [Google Scholar]
  44. Phillips, O.M. Flow and Reactions in Permeable in Permeable Rocks; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  45. Verma, A.; Pruess, K. Thermohydrologic conditions and silica redistribution near high-level nuclear wastes emplaced in saturated geological formations. J. Geophys. Res. 1988, 93, 1159–1173. [Google Scholar]
  46. Xu, T.; Ontoy, Y.; Molling, P.; Spycher, N.; Parini, M.; Pruess, K. Reactive transport Modeling of injection well scaling and acidizing at Tiwi field, Philippines. Geothermics 2004, 33, 477–491. [Google Scholar] [CrossRef]
  47. Wang, S.H.; Song, H.B.; Yan, W. The change of external conditions effects on the phase equilibrium curve of gas hydrate and the thickness of hydrate stability zone. Prog. Geophys. 2005, 20, 761–768. [Google Scholar]
Figure 1. Location of the research area and drilling sites in the Shenhu area, north slope of the South China Sea.
Figure 1. Location of the research area and drilling sites in the Shenhu area, north slope of the South China Sea.
Energies 10 01241 g001
Figure 2. Vertical versus horizontal fracture stimulation.
Figure 2. Vertical versus horizontal fracture stimulation.
Energies 10 01241 g002
Figure 3. Production well design and schematic of the marine hydrate deposit at the SH7 site.
Figure 3. Production well design and schematic of the marine hydrate deposit at the SH7 site.
Energies 10 01241 g003
Figure 4. Fracturing cracks diagram.
Figure 4. Fracturing cracks diagram.
Energies 10 01241 g004
Figure 5. (A) Hydrate distribution at one month; (B) Hydrate distribution at one year.
Figure 5. (A) Hydrate distribution at one month; (B) Hydrate distribution at one year.
Energies 10 01241 g005aEnergies 10 01241 g005b
Figure 6. (A) Pressure distribution at one month; (B) Hydrate distribution at one year.
Figure 6. (A) Pressure distribution at one month; (B) Hydrate distribution at one year.
Energies 10 01241 g006
Figure 7. Production rate and cumulative volume of CH4.
Figure 7. Production rate and cumulative volume of CH4.
Energies 10 01241 g007
Figure 8. The production rate and the cumulative volume curve of CH4 in the gas hydrate layer with three cracks under different crack spacings (Δl = 1 m, 2 m, 3 m, 5 m) were exploited by 5 MPa in five years.
Figure 8. The production rate and the cumulative volume curve of CH4 in the gas hydrate layer with three cracks under different crack spacings (Δl = 1 m, 2 m, 3 m, 5 m) were exploited by 5 MPa in five years.
Energies 10 01241 g008
Figure 9. The hydrate distribution in the gas hydrate layer with three cracks under different crack spacings (Δl = 1 m, 2 m, 3 m, 5 m) exploited by 5 MPa. (A) Hydrate distribution at one year; (B) Hydrate distribution at five years.
Figure 9. The hydrate distribution in the gas hydrate layer with three cracks under different crack spacings (Δl = 1 m, 2 m, 3 m, 5 m) exploited by 5 MPa. (A) Hydrate distribution at one year; (B) Hydrate distribution at five years.
Energies 10 01241 g009
Table 1. Production trial properties.
Table 1. Production trial properties.
ParameterValue
Initial pressure P0 (at base of HBL)13.83 MPa
Initial temperature T0 (at base of HBL)14.15 °C
Depth of seafloor1108 m
Thermal gradient0.0433 °C/m
HBL thickness ZH22 m
Depth of HBL H1155–177 m
Initial saturation in the HBLSH = 0.44; SA = 0.56
Gas composition100% CH4
Porosity Φ0.38
Water salinity (mass fraction) Xs0.0305
Intrinsic permeability kx = ky = kz7.5 × 10−14 m2
Permeability of fracture fracturing5.2 × 10−13 m2
Grain density ρR2600 kg/m3
Dry thermal conductivity Kdry1.0 W/(kg·°C)
Wet thermal conductivity Kwet3.1 W/(kg·°C)
Production pressure ΔP5 MPa
Composite thermal conductivity model K θ = K dry + ( S A + S H )
( K wet K dry ) + α S I K I
Capillary pressure model P cap = P 01 [ ( S * ) 1 / λ 1 ] 1 λ
S * = ( S A S irA ) / ( S mxA S irA ) .
SirA0.29
λ0.45
P01105 Pa
Relative permeability model K rA = ( S * ) n
K rG = ( S G * ) n G
S A * = ( S A S irA ) / ( 1 S irA )
S G * = ( S G S irG ) / ( 1 S irG )
EPM #2 model
N3.572
nG3.572
SirA0.30
SirG0.05
Table 2. Parameters of cracks and original formations.
Table 2. Parameters of cracks and original formations.
ParameterValue of CracksValue of Original Formation
Crack quantity1, 3, 50
Crack spacing1 m, 2 m, 3 m, 5 m0
Crack permeability5.2 × 10−13 m2 (10 mm)7.5 × 10−14 m2
Table 3. Simulation number and crack variable settings.
Table 3. Simulation number and crack variable settings.
GroupSimulation NumberCrack QuantityCrack Height/mmCrack Spacing/m
11-1000
1-21100
1-33102
1-45101
22-13101
2-23102
2-33103
2-43105

Share and Cite

MDPI and ACS Style

Chen, C.; Yang, L.; Jia, R.; Sun, Y.; Guo, W.; Chen, Y.; Li, X. Simulation Study on the Effect of Fracturing Technology on the Production Efficiency of Natural Gas Hydrate. Energies 2017, 10, 1241. https://doi.org/10.3390/en10081241

AMA Style

Chen C, Yang L, Jia R, Sun Y, Guo W, Chen Y, Li X. Simulation Study on the Effect of Fracturing Technology on the Production Efficiency of Natural Gas Hydrate. Energies. 2017; 10(8):1241. https://doi.org/10.3390/en10081241

Chicago/Turabian Style

Chen, Chen, Lin Yang, Rui Jia, Youhong Sun, Wei Guo, Yong Chen, and Xitong Li. 2017. "Simulation Study on the Effect of Fracturing Technology on the Production Efficiency of Natural Gas Hydrate" Energies 10, no. 8: 1241. https://doi.org/10.3390/en10081241

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