Gas hydrates in the Shenhu area are mainly hosted in clayey silt sediments, which have the relatively high irreducible fluid saturation and gas entry pressure. And then, they will have an impact on gas production from hydrate-bearing clayey silt sediments, which was evaluated by the numerical simulations of SH2 site in Shenhu area in this paper. The results showed that, with the increase in irreducible water saturation and irreducible gas saturation, the amount of water production and gas production was obviously reduced. When the irreducible water saturation increased from 0.10 to 0.50, the cumulative CH4 production volume decreased from 1668799 m3 to 1536262 m3, and the cumulative water production volume dropped from 620304 m3 to 564797 m3, respectively. When the irreducible gas saturation increased from 0.01 to 0.05, the cumulative CH4 production volume dropped from 1812522 m3 to 1622121 m3, and the cumulative water production volume dropped from 672088 m3 to 600617 m3, respectively. In addition, the capillary pressure increased obviously with the increase in gas entry pressure, but the effect on gas production was small and the effect on water production could be negligible. In conclusion, irreducible water and gas saturation had an important effect on the gas production from gas hydrate, whereas the effects of gas entry pressure could be ignored.
National Natural Science Foundation of China418721855189091441902321519913641. Introduction
Natural gas hydrate (NGH) is an ice-like crystalline compound formed by water molecules and natural gas molecules, such as methane, ethane, and propane, under high pressure and low temperature [1]. NGH is considered one of the potentially clean energies in future, and it is widely distributed in marine sediments and the permafrost regions [2, 3]. It has attracted more and more attentions from researchers all over the world. According to the preliminary estimates, the total amount of energy hosted in the NGH reservoir is two times that of conventional oil and gas resources [4]. However, because of the characteristics of unconsolidated and permeable overburden layer, submarine hydrate decomposition may cause geological hazards and climate change [5–7].
Different from the production of oil and natural gas resources, gas hydrate should be decomposed into gas and water at first, which is an endothermic phase transition. There are four main methods for hydrate production, including depressurization, thermal stimulation, CO2-CH4 replacement, and inhibitor injection [8–10]. In addition, Chinese scientists have proposed a new type of production method called solid fluidization, which transfer natural gas hydrate through mechanical crushing fluidization into a gas-liquid-solid multiphase pipeline [11]. Among them, depressurization method was successfully applied to the offshore gas hydrate production tests in Japan and China. In 2013, Japan’s first marine hydrate production lasted for six days, and the total gas volume reached 1.195×105 m3 [12]. In 2017, two production wells were drilled in Japan’s second offshore marine hydrate production test; the gas production of the first production well reached 3.5×104 m3 in 12 days and the total gas production of the second production well was up to 2×105 m3 for 24 days [13]. In July 9, 2017, China successfully completed its first offshore hydrate production test in the Shenhu area, northern South China Sea, which was the first time in the world to produce gas from clayey silt NGH reservoir [14]. The total volume of gas production was up to 3.09×105 m3 in 60 days [13, 14]. In May 2017, China successfully produced 81 m3 of gas in the northern South China Sea by solid fluidization [11]. From February 17 to March 30, 2020, China’s second offshore hydrate production test achieved a total gas production of 1.4986×106 m3 by using a horizontal production well, with an average daily gas production of 3.57×104 m3 [15]. Although these four offshore methane hydrate productions have made great progress, it is far from being commercialized and economical.
China conducted three hydrate drilling expeditions in 2007 (GMGS-1, Guangzhou Marine Geological Survey), 2015 (GMGS-3), and 2016 (GMGS-4) in the Shenhu area, and gas hydrate samples were successfully obtained in these three expeditions [16]. In GMGS-1, gas hydrates were found at sites SH2, SH3, and SH7 [17, 18], and this area is the typical Class 2 hydrate reservoir, which have a mobile water zone under hydrate-bearing layer (HBL). Based on the analysis of core samples, geophysics, and geochemistry, a series of important data, such as hydrate saturation, porosity of sediments, reservoir temperature, and pressure, was obtained. By depressurization experiments to estimate gas hydrate saturation in site SH2, it was concluded that the maximum hydrate saturation could get up to 0.48 [19]. The particle size composition of hydrate-bearing sediments in site SH2 is that silts (0.04–0.063 mm) accounts for 73.94%, clay (<0.04 mm) accounts for 24.68%, and sand (>0.063 mm) accounts for 1.4%, which means gas hydrates are hosted in clayey silt sediments [20]. And the results of X-ray diffraction pattern indicate that the clay minerals in the hydrate layer are mainly montmorillonite (illite/smectite mixed layer), illite, chlorite, and kaolinite [21]. The features of small particle size and high clay content would cause the multiphase flow in the process of gas production more complicated, which mainly includes three aspects. The first is increasing irreducible water saturation. Water will be adsorbed on the surface of sediment particles because of its hydrophilicity [22], and some of the irreducible water is formed by capillarity; in addition, clay will absorb a large amount of water due to strongly hydrophilic [23]. The second is increasing irreducible gas saturation. Previous studies about shale gas have shown that the gas sorption capacity could be influenced by clay mineral composition [24, 25]. The sorption capacity of clays for methane increases in the order illite < chlorite < kaolinite < illite/smectite mixed layer ≪ montmorillonite [25], which means that the clay mineral composition will lead to the different sorption capacities of methane for the different hydrate-bearing clayey silt sediments. In addition, the absorbed gas will decrease with the increase in water saturation [26, 27]. The third is increasing gas entry pressure. Gas entry pressure depends on the radius of the largest interconnected pathway through the sediment matrix, the water wettability of the sediments, and the interfacial tension between the fluid phases, so the small particle diameter and high clay content could increase the gas entry pressure [28]. Those changes caused by the small particle size and high clay content could change the gas and water relative permeability and capillary pressure, which would further affect gas production from HBL.
Based on the data obtained from hydrate drilling expeditions in the South China Sea, a lot of numerical simulations were carried out. Su et al. [29] studied the gas production behaviors of the NGH reservoir in site SH3 by depressurization method in a vertical well, and it was found that the gas production could not reach the industrial production according to the relative criterion and the absolute criterion. Li et al. [30] used a single horizontal well in the middle of the HBL to produce gas from the NGH reservoir in site SH7 and found that the gas production rate could reach 3.46×104 m3/d when the horizontal well length was 1000 m [30]. Jin et al. [31] investigated the enhancement of hot water injection on gas production from the NGH reservoir in site SH2 and found that the gas production increased with hot water temperature and injection pressure. Sun et al. [32] studied the geomechanical response of the NGH reservoir during gas production under different conditions including hydrate saturation heterogeneity, hydrate formation permeability, and gas formation permeability and found that the largest seafloor subsidence was only about 0.14 m. Konno et al. [12] simulated Japan’s first marine hydrate production through numerical simulation. The simulation results showed that the gas production gradually increased with the expansion of the decomposition area, and the increase in hydrate formation permeability increased the ratio of gas to water. Yang et al. [33] studied gas production by depressurization combined with warm water flooding, and the results showed that this method had the characteristics of high efficiency, low cost, and high safety. Ruan et al. [34] investigated the influence of factors such as initial hydrate saturation, permeability reduction index, absolute/relative permeability, and intrinsic porosity on hydrate production, the authors also pointed out the importance of relative permeability model and permeability decay index for numerical simulation. Ma et al. [35] found the increase of water permeability reduction index from 2.5 to 4.5 could make the total gas production decreases from 1.93×106m3 to 1.34×106m3.
Although there are many articles studying the effects of various factors on gas hydrate production through numerical simulation, little research has studied the effects of hydrogeological parameters such as irreducible water saturation, irreducible gas saturation, and gas entry pressure on gas hydrate production. The main purpose of this study is to evaluate the effects of different irreducible fluid saturation and gas entry pressure on gas and water production, gas and water relative permeability, and capillary pressure.
2. Simulation Model2.1. Model Parameters
In this study, Tough+Hydrate was used to conduct the numerical simulations. This software could model the nonisothermal gas release, phase behavior, and flow of fluids and heat in complex geologic media [36]. The schematic of the simulated NGH reservoir and well design in site SH2 in this study is shown in Figure 1. The NGH reservoir in site SH2 could be divided into three layers, including permeable overburden layer, HBL, and permeable underlying layer. The thickness of HBL is assumed to be 40 m, and its intrinsic permeability kH, porosity, and average hydrate saturation are set to 10 mD [19], 38%, and 0.40 [36, 37], respectively. The thicknesses of the permeable overburden layer and the permeable underlying layer are all set to 30 m, which could satisfy the heat and mass exchange in the process of 10 years’ gas production from the NGH reservoir. Other physical characteristics of overburden layer and underlying layer, such as porosity and permeability, are the same as the HBL [19, 29, 36]. The production well is in the center of this cylindrical hydrate reservoir with a radius of 0.1 m and a special no-flow boundary, which is called the Neumann boundary, is applied at rmax=100m [37]. The upper and lower boundaries are assumed to be constant temperature and pressure boundary conditions, and these boundaries are permeable. The hydrate used in this numerical simulation is methane hydrate. The salinity of pore water is 3.05% [32]. Other properties and conditions of the numerical model used in this study are listed in Table 1.
Schematic of the simulated NGH reservoir and well design in site SH2.
Main hydrate reservoir properties and conditions in site SH2.
Overburden layer thickness, ∆ZO
30 m
Underlying layer thickness, ∆ZU
30 m
HBL thickness, ∆ZH
40 m
Completion interval (LP)
20 m (in the middle of HBL)
Borehole radius
0.1 m
Initial pressure, PB (at base of HBL)
14.97 MPa
Initial temperature, TB (at base of HBL)
14.87 °C
Gas composition
CH4
Porosity, Φ (all formations)
38%
Initial saturation in the HBL
SW=0.60, SH=0.40
Intrinsic permeability, kx=ky=kz (all formations)
10 mD
Water salinity (mass fraction)
3.05%
Grain density, ρR (all formations)
2600 kg/m3
Dry thermal conductivity, kdry (all formations)
1.0 W/m/°C
Wet thermal conductivity, kwet (all formations)
3.1 W/m/°C
Composite thermal conductivity model
kθ=kdry+SA+SHkwet−kdry+ϕSIλI
Capillary pressure model
Pcap=−P0S∗−1/λ−11−λS∗=SW−SirW/SmxW−SirW
λ
0.45
Relative permeability model
krW=SW−SirW/1−SirWnkrG=SG−SirG/1−SirGnG
n
3.572 [38]
nG
3.572 [38]
As shown in Figure 1, without confining impermeable overburden layer, a large volume of undesirable water will flow into production well, thereby reducing gas production and the effectiveness of depressurization. Meanwhile, the lack of a confining impermeable overburden layer will cause gas leak into the ocean through the overburden layer. In order to prevent gas leaking into the ocean and undesirable water production, the perforated interval does not completely pass through the whole HBL [29, 38]. In this study, the perforation interval has a total length of Lp=20 m, which is located at the middle of HBL [35]. Previous studies have proven that depressurization is the most effective and economical method of hydrate production [39–41], so the method of depressurization is used in this study. The production well has a constant pressure of 3 MPa, which is above the quadruple point of methane gas hydrate. This pressure will accelerate the decomposition of hydrate and avoid the formation of ice [42, 43]. The entire production lasts 3650 days.
The values of irreducible water saturation SirW, irreducible gas saturation SirG, and gas entry pressure P0 will change with the particle size, clay content, and hydrate saturation. However, the change rule of relevant data with geological parameters is not given in the previous research through experiments. Therefore, in this study, those parameters are assumed to be fixed values to study the influence of these parameters on hydrate production. The relative permeability model used in this simulation is as follows [44]:
(1)KrW=SW−SirW1−SirWn,KrG=SG−SirG1−SirGnG,KrH=0,where KrW, KrG, and KrH are water relative permeability, gas relative permeability, and hydrate relative permeability, respectively, SW and SG are water saturation and gas saturation, respectively. Hydrate exists as a solid in pores, so its relative permeability is 0.
The capillary pressure function used is as follows [45]:
(2)Pcap=−P0S∗−1/λ−11−λ,S∗=SW−SirWSmxW−SirW,−Pmax≤Pcap≤0,where Pcap is capillary pressure, S∗ is scaled saturation, and Pmax is maximum capillary pressure.
In these two models, SirW, SirG, and P0 are key parameters which could affect the multiphase flow by affecting fluid relative permeability and capillary pressure. In the cases with different SirW, SirW is set to 0.10, 0.30, and 0.50, respectively, while SirG and P0 are set to 0.05 and 0.1 MPa, respectively; In the cases with different SirG, SirG is set to 0.01, 0.03 and 0.05, respectively, while SirW and P0 are set to 0.30 and 0.1 MPa, respectively. In the cases with different P0, P0 is set to 0.01, 0.10, and 0.50 MPa, respectively, while SirW and SirG are set to 0.30 and 0.05, respectively. The parameters in this paper are shown in Table 2.
Simulation parameters in different cases.
Simulation ID
Parameters
Bottom hole pressure (MPa)
Irreducible water saturation
Irreducible gas saturation
Gas entry pressure (MPa)
1
Irreducible water saturation
3
0.1
0.05
0.1
2
0.3
3
0.5
4
Irreducible gas saturation
3
0.3
0.01
0.1
5
0.03
6
0.05
7
Gas entry pressure
3
0.3
0.05
0.01
8
0.1
9
0.5
2.2. Initial and Boundary Conditions
Initial pressure and temperature at the bottom of HBL are assumed to 14.97 °C and 14.87 MPa, respectively, which are also used by other papers [37]. The overburden layer above the HBL is permeable, and the pore water in the stratum could be exchanged with seawater, which means pore water pressure could be determined by hydrostatic pressure formula. So, the following formula is used to calculate the initial pore pressure of whole model:
(3)Ppw=PB+ρswgΔz×10−6,where Ppw is the hydrostatic pore water pressure in MPa, PB is the pressure at the bottom of the HBL, and its value is 14.87 MPa, and g is the gravity acceleration. ∆z is the distance in m from the bottom of HBL, and it is negative above HBL and positive below HBL. ρsw is the average pore water density, assumed to be 1035 kg/m3 [30]. The temperature gradient is set to 4.695 °C/100 m. According to the hydrostatic pressure formula and temperature gradient, the temperature and pressure of the upper and lower boundaries could be calculated to initialize the temperature and pressure of the whole simulation domain through the self-balancing function of software.
3. Results and Discussions3.1. Spatial Distributions of Physical Properties in Reservoir
Figure 2 shows the evolution of the pressure distribution over time within the entire reservoir in the case with SirW=0.30, SirG=0.05, and P0=0.1 MPa under the production pressure of 3 MPa. It could be seen from Figure 2 that pressure distribution has the following characteristics: (1) in the first ten days of production, the water near the production well flows into the production well under the pressure gradient, and the effect of the depressurization is only within 10 meters; as shown in Figure 2(a), due to the large range of the low-pressure zone below 12 MPa in the initial stage, a large number of hydrates are decomposed, so that the rate of gas production at this stage is high; (2) there is an obvious pressure drop around the production well, and the range of depressurization increases with time; however, because of the low permeability of HBL, the maximum range of pressure drop is within 30 m; (3) it could be seen from the figure that the influence range of the low-pressure zone below 12 MPa increases between 10–365 d, which is due to the gradual extraction of fluid from the reservoir resulting in the expansion of low-pressure zone; but after 365 d, more and more water from the overburden and underlying layers enters the production well, which leads to the increase in water production rate and has a negative effect on depressurization, so the range of the low-pressure area decreases.
Evolution of the pressure distribution P on the 10th day, 365th day, 1825th day, and 3650th day in the case with SirW=0.30, SirG=0.05, and P0=0.1 MPa.
The location of decomposed hydrate zone could be analyzed by hydrate saturation distribution. The following conclusion could be concluded from Figure 3: (1) during the 3650 days of production, the main range of hydrate decomposition area is limited to 20 meters; the existence of hydrates reduces the effective permeability of the low permeability clayey slit sediments, which affects the range of the depressurization and increases the time of the decomposed gas flows to the production well; (2) the decomposition of hydrates is also accompanied by the formation of secondary hydrates; this phenomenon is mainly due to the fact that the gas generated by hydrate decomposition enters the undecomposed area by diffusion or convection, thereby causing hydrate reformation phenomenon; (3) in the initial stage, the hydrate decomposition zone is mainly concentrated around the production well; and (4) on the 365th day, hydrate decomposition occurs at both the upper and lower parts of the HBL, and the lower part of the HBL decomposes rapidly, which is mainly caused by the high-temperature fluid of the underlying and overburden layers flowing into the production well, and the temperature of the underlying layer is higher than that of the overburden layer [30].
Evolution of the hydrate saturation distribution SH on the 10th day, 365th day, 1825th day, and 3650th day in the case with SirW=0.30, SirG=0.05, and P0=0.1MPa.
Figure 4 shows the evolution of the gas saturation distribution SG on the 10th day, 365th day, 1825th day, and 3650th day in the case with SirW=0.30, SirG=0.05, and P0=0.1MPa. The following characteristics could be concluded from Figure 4: (1) the gas saturation in the whole process of production is almost below 0.2, and the gas saturation gradually decreases with the production; (2) in the early stage of production, gas is mainly distributed around the production well; but when the water from the underlying and overburden layers flows into the production well, the gas is mainly distributed on the lower decomposition edge; and (3) it could be seen from Figure 4(a) and Figure 4(b) that the range of gas is larger than hydrate decomposition area, which confirms the reason for the formation of secondary hydrate.
Evolution of the gas saturation distribution SG on the 10th day, 365th day, 1825th day, and 3650th day in the case with SirW=0.30, SirG=0.05, and P0=0.1MPa.
Figure 5 shows the evolution of the water saturation distribution SW on the 10th day, 365th day, 1825th day, and 3650th day in the case with SirW=0.30, SirG=0.05, and P0=0.1MPa. It can be seen from the partial enlarged figure in Figure 5(a) that with the increase in the distance from the production well, the SW in the area near the production well decreases first and then increases. The formation of the area with low water saturation is due to the expansion of the gas produced by hydrate decomposition under the action of pressure reduction, thus displacing the water in the pores. As shown in Figure 5(b), water from the overburden and the underlying layers begins to invade the hydrate decomposed area around the production well. Because the temperature of the water in the pores of the underlying layer is higher than that of the water in the overburden layer, hydrate decomposition in the lower part of the HBL is faster than that in the upper part. As shown in Figures 5(c) and 5(d), with the development of gas production from the hydrate reservoir, the phenomenon of water intrusion becomes more and more obvious.
Evolution of the water saturation distribution SW on the 10th day, 365th day, 1825th day, and 3650th day in the case with SirW=0.30, SirG=0.05, and P0=0.1MPa.
3.2. The Effects of Irreducible Water Saturation
Irreducible water saturation SirW is one of the fundamental parameters in characterizing reservoirs and also determined as the key factor in reservoir calculation and fluid type identification. High clay content and small particle size could increase SirW. Therefore, the SirW in these three cases are assumed to be 0.10, 0.30, and 0.50, respectively, to evaluate the effect of SirW on gas production from NGH reservoir. In these cases, the SirG is assumed to be 0.05, and the value of the P0 is assumed to be 0.1 MPa.
Figure 6 shows the variation of relative permeability of gas and water (KrG and KrW) with water saturation in different cases with different SirW. According to the number of coexisting movable phases in the flow process, it could be divided into the gas single-phase flow region (A1, A2, and A3), gas-liquid two-phase flow region (B1, B2, and B3), and water single-phase flow region (C). When the water saturation is less than SirW, the water is adsorbed on the surface of sediment particles or trapped in the pore throat, so it cannot flow. When the water saturation is greater than SirW, the water starts to flow. With the increase in water saturation, the water relative permeability increases gradually. As the increase in water saturation will occupy the pore space, the gas saturation will decrease, so the gas relative permeability decreases. With the increase in SirW, the region of the gas single-phase flow becomes larger and the two-phase flow region decreases. This is because when the water saturation is low, all water is adsorbed by sediment particles or trapped by small pores or throats, thus becoming irreducible water. Meanwhile, the increase in SirW could increase the proportion of immovable water, thus reducing the relative water permeability. The relative permeability curves of gas and water become steeper with the increase in SirW, mainly because the irreducible water adsorbed by clay particles makes the effective pore diameter smaller.
Relationship between gas, water relative permeability, and water saturation in different cases with different irreducible water saturation SirW.
As shown in Figure 7, with the decrease in water saturation and increase in gas saturation, capillary pressure increases gradually. When the water saturation decreases to a certain extent, the capillary pressure reaches the maximum value, and its maximum value remains unchanged as the water saturation continues to decrease to SirW. The irreducible water occupies the pore space, which makes the pore of the movable water occupied decreases and the slope of capillary pressure increases. In the hydrate-bearing sediments, the particle diameter, clay content (especially montmorillonite content), particle gradation, and consolidation degree of the sediment could change the irreducible water content, thus affecting the change of capillary pressure.
Relationship between capillary pressure and water or gas saturation in different irreducible water saturation cases.
Figure 8 compares total gas production rate QPT and total gas production VPT in three different cases with different SirW. It could be clearly seen from Figure 8 that QPT is obviously declining with the increase in SirW. At the end of 3650 days’ production, the total gas production of three different SirW cases is 1668799, 1622121, and 1536262 m3, respectively. When SirW decreases from 0.50 to 0.10, the total gas production increases by 8.63%. It could be easily drawn that the whole production could be divided into three stages according to the difference in QPT between different SirW cases. In the initial stage, the difference in QPT between different SirW cases is very small, and then in the following stage, the differences in QPT between different SirW gradually increase. However, in the last stage of production, the difference in QPT between different SirW cases gradually decreases. In the initial stage, gas produced by hydrate decomposition could flow into production well more easily as it is close to the production well, which means the gas production rate is less affected by SirW. As production progresses, gas needs to migrate longer distances to reach the production well, so the impact of gas relative permeability on gas production from hydrate reservoir increases. And the gas relative permeability decreases with the increase in SirW, which leads to the decrease in gas production rate in the second stage. In the final stage, the decrease in pressure gradient slows the decomposition rate of hydrate, resulting in lower gas saturation. As could be seen from Figure 6, when gas saturation is low, the difference between gas relative permeability becomes smaller, so the difference of gas production rates between different cases at the end of gas production from hydrate reservoir becomes smaller.
Evolution of total gas production rate QPT and total gas production VPT in three different cases with different SirW.
It could be seen from Figure 9 that, as the range of pressure reduction increases gradually, more water flows into the production well, so that the water production rate increases gradually. In the early stage of gas production from the hydrate reservoir, the increase rate of water production rate is faster, and then the increase rate gradually slows down. As shown in Figure 9, the increase in SirW will cause the decrease in QW. At the end of 3650 days’ production, the total water production of three different SirW cases is 620304, 600617, and 564797 m3, respectively. When the irreducible water saturation decreases from 0.50 to 0.10, the water production increases by 9.83%. On the one hand, the water flow into the production well decreases due to the decrease in the movable water proportion in the pore, and on the other hand, the water relative permeability decreases due to the increase in SirW.
Evolution of total water production rate QW and total water production VW in three different cases with different SirW.
Figure 10 shows the comparison of gas relative permeability KrG, water relative permeability KrW, and capillary pressure Pcap distribution between different cases with different SirW on the 3650th day. As seen in Figure 10, all the shapes of KrG, KrW and Pcap distribution are like a tail, which is consistent with the gas distribution in Figure 4(d). It could be seen from Figures 10(a)–10(c) that the difference of KrG distribution between different SirW cases is very small, which may be the main reason of the small difference of QPT between different SirW cases in the last stage of production. And it is also could be deduced that KrG decreases as SirW increases in the second stage of production in Figure 8. KrG near the well is higher than other locations, which is mainly due to the fact that more gas accumulate near the perforated section under the influence of depressurization. The low permeability of clayey silt sediments reduces the depressurization effect and fluid flow rate in the hydrate reservoir, resulting in low hydrate decomposition rate and low gas saturation. As shown in Figure 6, low gas saturation could reduce the KrG. Therefore, KrG in Figure 10 is very small, and its value is mainly distributed between 0 and 0.0002, and low KrG is the root cause of low gas production rate in the later period of production. It could be seen from Figures 10(d) to 10(f) that KrW drops significantly with the increase in SirW, which results to the decline of QW with the increase in SirW in Figure 9. The increase in SirW leads to the decrease in movable water, so the KrW decreases greatly. KrW near the well is lower than other locations, and this is mainly ascribed that more gas and less water accumulate near the perforated section. In addition, KrW is mainly distributed between 0.45 and 0.95, much larger than KrG. The higher KrW leads to an increase in water production and a decrease in production efficiency. It is worth noting that Pcap, in Figures 10(g)–10(i), increases moderately with the increase in SirW. This may be due to the fact that the increase in SirW reduces the effective porosity and narrows the fluid flow path, which leads to an increase in Pcap.
Comparison of gas relative permeability KrG, water relative permeability KrW, and capillary distribution Pcap between three cases with different irreducible water saturation SirW on the 3650th day.
3.3. The Effects of Irreducible Gas Saturation
Due to the characteristics of higher clay content, small particle diameter, and more small pore throats, clayey silt sediments may contain more irreducible gas. To study the effects of the increase in SirG on gas production from hydrate reservoir, the values of SirG in three cases are assumed to be 0.01, 0.03, and 0.05. In these cases, SirW is assumed to be 0.30, and the value of the P0 is assumed to be 0.1 MPa.
In the water relative permeability equation used in this paper, the change of SirG has no effect on the relative permeability of water. As shown in Figure 11, because the gas produced by hydrate decomposition is trapped in the reservoir and could not be a part of movable fluids, increasing SirG will decrease the gas relative permeability. When the saturation of water in porous media is greater than 0.6, the relative permeability of gas drops below 0.13.
Relationship between gas, water relative permeability, and water saturation under different irreducible gas saturation.
As shown in Figure 12, when the values of SirG are equal to 0.01, 0.03, and 0.05, QPT decreases from 1238 m3/d, 1109 m3/d, and 995 m3/d on the 1st day to 330 m3/d, 324 m3/d, and 319 m3/d on the 3650th day, respectively. When the SirG decreases from 0.05 to 0.01, the total gas production in 3650 days increases from 1622121 m3 to 1812522 m3, increasing by about 11.74%. The increase in SirG leads to the decrease in gas production, which is caused by the fact that the increase in SirG causes more free gas produced by hydrate decomposition to remain in the hydrate reservoir. In addition, as the production progresses, the difference of gas production rate between different SirG gradually decreases. There are two reasons for this phenomenon. One is that with the expansion of the pressure reduction range, the pressure gradient inside the reservoir decreases, resulting in the slowing down of the gas flow rate. The second is that the pressure drop of the hydrate decomposition front becomes smaller as the hydrate decomposition front is gradually far away from the production well, so that the decomposition driving force is smaller. Further, this leads to a decrease in gas saturation in the hydrate reservoir, which reduces the difference in KrG between different SirG as shown in Figure 11.
Evolution of total gas production rate QPT and total gas production VPT in three different cases with different SirG.
Figure 13 shows a clear comparison between three different cases with different SirG in QW and VW. As shown in the figure, SirG also has a significant impact on water production. When the SirG decreases from 0.05 to 0.01, the total water production in 3650 days increases from 600617 m3 to 672088 m3, increasing by about 11.90%. The increase in SirG causes a significant decrease in QW, which is due to the fact that more gas occupies the pores of hydrate-bearing sediments, thus reducing the reservoir’s effective porosity and permeability.
Evolution of total water production rate QW and total water production VW in three different cases with different SirG.
Figure 14 shows the evolution of KrG, KrW, and Pcap distribution at different SirG cases on the 3650th day. As the Figure 14 shows, (1) on the 3650th day, the difference of KrG between three different cases is very small, leading to a very little difference in gas production rate in the later stage of production, which is due to the low gas saturation makes smaller difference in KrG; (2) the increase in SirG reduces KrW, which is due to that the increase in SirG increases gas saturation in the pores of sediments, thus reducing water saturation and KrW; (3) the increase in SirG makes capillary pressure increase significantly, which is because the irreducible gas occupies pore space, thus making pore throat smaller and increasing capillary pressure.
Comparison of gas relative permeability KrG, water relative permeability KrW, and capillary distribution Pcap between there cases with different irreducible gas saturation SirG on the 3650th day.
3.4. The Effects of Gas Entry Pressure
Figure 15 shows the relationship between capillary pressure and water or gas saturation in different gas entry pressure (P0) cases. Gas entry pressure represents the pressure of gas starting to enter the porous medium, which corresponds to the capillary pressure at the maximum throat. With the increase in P0, the range of maximum capillary pressure increases, and the capillary pressure in the same water saturation condition increases. When the gas entry pressure is 0.5 MPa, the capillary pressure increases rapidly to the maximum value with the increase in water saturation, which means the diameter of pore throat is smaller. When the gas entry pressure is 0.01 MPa and 0.1 MPa, the capillary pressure rises slowly and then rapidly to 1 MPa, in which the slowly rising stage represents the larger pore throat diameter.
Relationship between capillary pressure and water or gas saturation in different gas entry pressure cases.
Figure 16 compares QPT and VPT in three different cases with different gas entry pressure P0. It could be seen from the figure that the increase in P0 could slightly increase gas production. When the P0 increases from 0.01 MPa to 0.1 MPa, the gas production in 3650 days is only reduced by 0.46%, which could be ignored; when the gas entry pressure increases from 0.1 MPa to 0.5 MPa, the gas production in 3650 days is only reduced by 2.36%. The gas produced by hydrate decomposition will occupy the original space of the hydrate preferentially. Then, the gas will expand under the effect of depressurization. When the gas pressure exceeds the sum of water pressure in the pore throat and the capillary pressure, the water will drive the water in the capillary tube to flow. In this paper, the gas expansion effect is obvious due to the large depressurization amplitude in the hydrate reservoir, which is easy to overcome the flow resistance caused by capillary pressure. At the same time, the gas entry pressure has no effect on the gas relative permeability and reservoir permeability in this paper. Therefore, the gas entry pressure has no obvious effect on gas production. Figure 17 shows the evolution of total water production rate QW and total water production VW in three different cases with different gas entry pressure P0. The gas entry pressure has little effect on the water production, which is because the gas entry pressure has little effect on the water relative permeability in this study.
Evolution of total gas production rate QPT and total gas production VPT in three different cases with different gas entry pressure P0.
Evolution of total water production rate QW and total water production VW in three different cases with different gas entry pressure P0.
Figures 18 and 19 show the evolution of gas relative permeability KrG, water relative permeability KrW, and capillary pressure Pcap distribution in different SirG cases on the 3650th day. It could be clearly seen that P0 reduces gas relative permeability on the 3650th day. However, the change of P0 does not change the effective space for fluid flow as irreducible water or gas saturation, so it has little effect on gas production. Meanwhile, as shown in Figures 18(d)–18(f), P0 has little effect on the relative permeability of water on the 3650th day, which is the main reason that the gas entry pressure has no effect on the water production. However, the increase in gas entry pressure will significantly increase the capillary pressure in hydrate-bearing sediments, causing the amount of gas production to slightly decline.
Comparison of gas relative permeability KrG and water relative permeability KrW distribution in three different cases with different gas entry pressure P0 on the 3650th day.
Comparison of capillary pressure Pcap distribution in three different cases with different gas entry pressure P0 on the 3650th day.
4. Conclusions
The effects of hydrogeological parameters such as irreducible water saturation, irreducible gas saturation, and gas entry pressure on gas hydrate production were studied in this paper. The main conclusions could be drawn as follows.
When the irreducible water saturation decreases from 0.10 to 0.50, the cumulative gas production decreases from 1668799 m3 to 1536262 m3, and the total water production drops from 620304 m3 to 564797 m3, respectively. The increase in the irreducible water saturation could occupy the pore space, so it will reduce the relative permeability of gas and water. At the same time, it will make the effective pore throat diameter smaller and increase the capillary pressure, so it will reduce the gas production and water production
When the SirG decreases from 0.05 to 0.01, the total gas production in 3650 days increases from 1622121 m3 to 1812522 m3, increasing by about 11.74%. When the SirG decreases from 0.05 to 0.01, the total water production in 3650 days increases from 600617 m3 to 672088 m3, increasing by about 11.90%. When the irreducible gas saturation increases, a large number of gas will be trapped in the pores, and the increase in irreducible gas saturation will also affect the relative permeability of the gas, so the gas production would decrease with the increase in irreducible gas saturation. The increase in the irreducible gas saturation will reduce the effective pore space and the effective permeability of the reservoir, so the water production will decrease
The change of gas entry pressure has a slight effect on gas production and has no effect on water production. The increase in gas entry pressure has no effect on the relative permeability of water and gas, so its effect on water and gas production is not obvious. However, the gas entry pressure will change the capillary pressure distribution in the reservoir, which has a slight impact on the gas production in this paper
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study has been supported by the National Natural Science Foundation of China (Grant No. 51991364, Grant No. 41902321, Grant No. 51890914, and Grant No. 41872185).
UddinM.CoombeD.LawD.GunterB.Numerical studies of gas hydrate formation and decomposition in a geological reservoir2008130310.1115/1.29569782-s2.0-57249108611ChattiI.DelahayeA.FournaisonL.PetitetJ. P.Benefits and drawbacks of clathrate hydrates: a review of their areas of interest2005469–101333134310.1016/j.enconman.2004.06.0322-s2.0-15844389145XuW.RuppelC.Predicting the occurrence, distribution, and evolution of methane gas hydrate in porous marine sediments1999104B35081509510.1029/1998JB900092LiX. S.XuC. G.ZhangY.RuanX. K.LiG.WangY.Investigation into gas production from natural gas hydrate: a review201617228632210.1016/j.apenergy.2016.03.1012-s2.0-84962345979NingF.YuY.KjelstrupS.VlugtT. J. H.GlavatskiyK.Mechanical properties of clathrate hydrates: status and perspectives2012556779679510.1039/c2ee03435b2-s2.0-84860362572HatzikiriakosS. G.EnglezosP.The relationship between global warming and methane gas hydrates in the earth199348233963396910.1016/0009-2509(93)80375-Z2-s2.0-0027790608KvenvoldenK. A.Potential effects of gas hydrate on human welfare19999673420342610.1073/pnas.96.7.34202-s2.0-003361655710097052LeeJ.ParkS.SungW.An experimental study on the productivity of dissociated gas from gas hydrate by depressurization scheme201051122510251510.1016/j.enconman.2010.05.0152-s2.0-77955467465WangY.FengJ. C.LiX. S.ZhangY.LiG.Large scale experimental evaluation to methane hydrate dissociation below quadruple point in sandy sediment201616237238110.1016/j.apenergy.2015.10.0992-s2.0-84945966531SunY.LiB.GuoW.LüX.ZhangY.LiK.WangP.JinG.JiaR.QuL.Comparative analysis of the production trial and numerical simulations of gas production from multilayer hydrate deposits in the Qilian Mountain permafrost20142145646610.1016/j.jngse.2014.09.0052-s2.0-84907997010ZhouS.LiQ.ChenW.Research on the solid fluidization well testing and production fo shallow non-diagenetic natural gas hydrate in the deep water area201729418KonnoY.FujiiT.SatoA.AkamineK.NaikiM.MasudaY.YamamotoK.NagaoJ.Key findings of the world’s first offshore methane hydrate production test off the coast of Japan: toward future commercial production20173132607261610.1021/acs.energyfuels.6b031432-s2.0-85017554850ChenL.FengY.OkajimaJ.KomiyaA.MaruyamaS.Production behavior and numerical analysis for 2017 methane hydrate extraction test of Shenhu, South China Sea201853556610.1016/j.jngse.2018.02.0292-s2.0-85042793036LiJ.YeJ.QinX.QiuH.WuN.The first offshore natural gas hydrate production test in South China Sea20181151610.31035/cg2018003People’s dailyMining combustible ice is expected to speed up again2020http://www.cgs.gov.cn/xwl/ddyw/202003/t20200331_504547.htmWuS.WangJ.On the China’s successful gas production test from marine gas hydrate reservoirs20186312810.1360/N972017-00645YangS.ZhangH.WuN.SuX.SchultheissP.HollandM.ZhangG.LiangJ.LuJ.RoseK.High concentration hydrate in disseminated forms obtained in Shenhu Area, North Slope of South China Sea2008Proceedings of the 6th International Conference on Gas Hydrates (ICGH 2008)July 2008Vancouver, BC, CanadaWangX.HutchinsonD. R.WuS.YangS.GuoY.Elevated gas hydrate saturation within silt and silty clay sediments in the Shenhu area, South China Sea2011116B510.1029/2010JB0079442-s2.0-79958029447SuZ.CaoY.WuN.HeY.Numerical analysis on gas production efficiency from hydrate deposits by thermal stimulation: application to the Shenhu area, South China Sea20114229431310.3390/en40202942-s2.0-79951614946ChenF.ZhouY.SuX.LiuG.LuH.WangJ.The variation of grain size in hydrate fomation and its relationship with hydrate saturation in Shenhu area, South China Sea20112011595100LuH.ChenH.ChenF.LiaoZ.Mineralogy characteristics of gas hydrate borehole sediments in Shenhu area, South China Sea200920092839LiuD.YueX. A.WangL.HouJ.ZhangJ. H.Model study on the irreducible water saturation by centrifuge experiments2006531–2778210.1016/j.petrol.2006.03.0252-s2.0-33746547426LiJ.LiX.WangX.LiY.WuK.ShiJ.YangL.FengD.ZhangT.YuP.Water distribution characteristic and effect on methane adsorption capacity in shale clay201615913515410.1016/j.coal.2016.03.0122-s2.0-84964378331ChengA. L.HuangW. L.Selective adsorption of hydrocarbon gases on clays and organic matter200435441342310.1016/j.orggeochem.2004.01.0072-s2.0-1842484066JiL.ZhangT.MillikenK. L.QuJ.ZhangX.Experimental investigation of main controls to methane adsorption in clay-rich rocks201227122533254510.1016/j.apgeochem.2012.08.0272-s2.0-84868532642HildenbrandA.KroossB. M.BuschA.GaschnitzR.Evolution of methane sorption capacity of coal seams as a function of burial history – a case study from the Campine Basin, NE Belgium200666317920310.1016/j.coal.2005.07.0062-s2.0-32644445572KroossB. M.van BergenF.GensterblumY.SiemonsN.PagnierH. J. M.DavidP.High-pressure methane and carbon dioxide adsorption on dry and moisture- equilibrated Pennsylvanian coals2002512699210.1016/S0166-5162(02)00078-22-s2.0-0036640880HanF.BuschA.van WageningenN.YangJ.LiuZ.KroossB. M.Experimental study of gas and water transport processes in the inter-cleat (matrix) system of coal: anthracite from Qinshui Basin, China201081212813810.1016/j.coal.2009.11.0012-s2.0-73649141264SuZ.HeY.WuN.ZhangK.MoridisG. J.Evaluation on gas production potential from laminar hydrate deposits in Shenhu area of South China Sea through depressurization using vertical wells201286-87879810.1016/j.petrol.2012.03.0082-s2.0-84860345107LiG.MoridisG. J.ZhangK.LiX. S.Evaluation of gas production potential from marine gas hydrate deposits in Shenhu area of South China Sea201024116018603310.1021/ef100930m2-s2.0-78649565940JinG.XuT.XinX.WeiM.LiuC.Numerical evaluation of the methane production from unconfined gas hydrate- bearing sediment by thermal stimulation and depressurization in Shenhu area, South China Sea20163349750810.1016/j.jngse.2016.05.0472-s2.0-84990068706SunJ.ZhangL.NingF.LeiH.LiuT.HuG.LuH.LuJ.LiuC.JiangG.LiangJ.WuN.Production potential and stability of hydrate-bearing sediments at the site GMGS3-W19 in the South China Sea: a preliminary feasibility study20178644747310.1016/j.marpetgeo.2017.05.0372-s2.0-85020734880YangH.BaiY. H.LiQ. P.The simulation of gas production from oceanic gas hydrate reservoir by the combination of ocean surface warm water flooding with depressurization20122851287129510.1007/s10409-012-0095-62-s2.0-84870378385RuanX.SongY.LiangH.YangM.DouB.Numerical simulation of the gas production behavior of hydrate dissociation by depressurization in hydrate-bearing porous medium20122631681169410.1021/ef201299p2-s2.0-84863338199MaX. L.LiB.JiaR.GuoW.SunY. H.Effect of water permeability reduction index on gas production from hydrate-bearing clayey-silt sediments by depressurization2018KnE Materials Science180189SuZ.HuangL.WuN.YangS.Effect of thermal stimulation on gas production from hydrate deposits in Shenhu area of the South China Sea201356460161010.1007/s11430-013-4587-42-s2.0-84876060992SuZ.MoridisG. J.ZhangK.WuN.A huff-and-puff production of gas hydrate deposits in Shenhu area of South China Sea through a vertical well201286-87546110.1016/j.petrol.2012.03.0202-s2.0-84860326998HuangL.SuZ.WuN. Y.Evaluation on the gas production potential of different lithological hydrate accumulations in marine environment20159178279810.1016/j.energy.2015.08.0922-s2.0-84946040522LiX.-S.WangY.LiG.ZhangY.Experimental investigations into gas production behaviors from methane hydrate with different methods in a cubic hydrate simulator20122621124113410.1021/ef201641v2-s2.0-84857394969FengJ. C.WangY.LiX.-S.Hydrate dissociation induced by depressurization in conjunction with warm brine stimulation in cubic hydrate simulator with silica sand201617418119110.1016/j.apenergy.2016.04.0902-s2.0-84964509977MoridisG. J.CollettT. S.BoswellR.KuriharaM.ReaganM. T.Toward production from gas hydrates: assessment of resources, technology and potentialSPE Unconventional Reservoirs Conference. Society of Petroleum EngineersFebuary 2008Keystone, CO, USA10.2118/114163-msKonnoY.UchiumiT.OyamaH.JinY.NagaoJ.MasudaY.OuchiH.Dissociation behavior of methane hydrate in sandy porous media below the quadruple point20122674310432010.1021/ef300628c2-s2.0-84863905930ZhaoJ.ZhuZ.SongY.LiuW.ZhangY.WangD.Analyzing the process of gas production for natural gas hydrate using depressurization201514212513410.1016/j.apenergy.2014.12.0712-s2.0-84921060382MoridisG. J.KowalskyM. B.PruessK.2008Earth Sciences Division, Lawrence Berkeley National Laboratory10.2172/927149van GenuchtenM. T.A closed-form equation for predicting the hydraulic conductivity of unsaturated soils198044589289810.2136/sssaj1980.03615995004400050002x2-s2.0-0019057216