Influence of Reservoir Stimulation on Marine Gas Hydrate Conversion Efficiency in Different Accumulation Conditions

In this paper, we used a method of combining reservoir stimulation technique (RST) with depressurization to investigate the conversion efficiency of marine natural gas hydrate (NGH) reservoirs in the Shenhu area, on the northern slope of the South China Sea, which differ in intrinsic permeability and initial NGH saturation conditions. We also analyzed the influence of the variable-stimulation effect on marine NGH conversion efficiency in different accumulation conditions, providing a reference scheme for improving the NGH conversion efficiency in the Shenhu area. In this work, we performed calculations for the variations in CH4 production rate and cumulative volume of CH4 in different initial NGH saturation, intrinsic permeability, and stimulation effect conditions. Variance analysis and range analysis methods were used to analyze the significance of these key factors and their interaction. Furthermore, we investigated the sensitivity of stimulation effect on NGH conversion efficiency. The simulation results showed that the stimulation effect has a significant influence on NGH conversion efficiency, and the influence of interaction between these three factors was not obvious. Possibly most importantly, we clarified that the sparsely fractured networks (N = 3) had a better effect for higher NGH conversion efficiency under higher saturation conditions. For lower permeability cases, the influence between sparsely (N = 3) and densely (N = 5) fractured networks were similar.


Background
Natural gas hydrates (NGH) are white and pale yellow, solid, ice-like cage type crystalline compounds [1]; formed by small-molecule gases such as light hydrocarbons, carbon dioxide, and water under low-temperature and high-pressure conditions.These are also known as "combustible ice" [2].NGH is an alternate energy resource with great reserves [3], very clean natural gas can be produced from NGH deposits, especially from sandy turbidites, from which conventional hydrocarbons can be produced [4].It is estimated that natural gas from NGH in sands are >40,000 Tcf [5,6].NGH should be converted in situ to its constituent gas and water [7].However, it is a potential geohazard as hydrates can trigger submarine landslides [8].In addition, hydrate dissociation releases methane in to the atmosphere which is a powerful greenhouse gas.Therefore, efficient and safe exploitation methods are the focus of NGH study.A number of conversion methods exist [9,10], but early production testing and modeling indicate that depressurization will be the ideal method to use [11,12].
The NGH enriched area of the Shenhu area is located on the northern slope of the South China Sea, between Xisha Trough and Dongsha Islands (Figure 1) [13,14].The NGH enrichment area is characterized by thickness of sedimentary formations, high content of organic matter, and high sedimentation; the thermogenic gas from deep formations and microbial and thermogenic gases from the shallow formations provided a sufficient gas source for the formation of NGH [15,16].In this area, the water depth is 1000-1700 m [17].The temperature of the ocean floor, heat flow, geothermal gradient and organic matter contents are 3.3-3.7 • C, 74-78 mW/m 2 , 43-67.7 • C/km and 0.46-1.4%,respectively, which satisfied the favorable geological and thermodynamic conditions of NGH formation and stability [14,15,18,19].
Energies 2018, 11, x 2 of 16 can be produced [4].It is estimated that natural gas from NGH in sands are >40,000 Tcf [5,6].NGH should be converted in situ to its constituent gas and water [7].However, it is a potential geohazard as hydrates can trigger submarine landslides [8].In addition, hydrate dissociation releases methane in to the atmosphere which is a powerful greenhouse gas.Therefore, efficient and safe exploitation methods are the focus of NGH study.A number of conversion methods exist [9,10], but early production testing and modeling indicate that depressurization will be the ideal method to use [11,12].
The NGH enriched area of the Shenhu area is located on the northern slope of the South China Sea, between Xisha Trough and Dongsha Islands (Figure 1) [13,14].The NGH enrichment area is characterized by thickness of sedimentary formations, high content of organic matter, and high sedimentation; the thermogenic gas from deep formations and microbial and thermogenic gases from the shallow formations provided a sufficient gas source for the formation of NGH [15,16].In this area, the water depth is 1000-1700 m [17].The temperature of the ocean floor, heat flow, geothermal gradient and organic matter contents are 3.3-3.7 °C, 74-78 mW/m 2 , 43-67.7 °C/km and 0.46-1.4%,respectively, which satisfied the favorable geological and thermodynamic conditions of NGH formation and stability [14,15,18,19].

Hydrates in the Shenhu Area
In 2007, NGH samples were successfully drilled from the SH2, SH3, and SH7 sites of GMGS-1 research area in the Shenhu area, in the northern South China Sea (Figure 1) [14].The drilling results showed that the top of the hydrate-bearing layers (HBL) are located at 115-229 m below the ocean floor [14,15].In 2015, a total of 23 wells across 19 sites, were drilled in GMGS-3 research area in the Shenhu area in the northern South China Sea (Figure 1) [15,20,21].The well logging data showed that all of these stations have NGH.Hydrate samples were collected in four of these wells.The hydrate layers are 13-70 m thick.The NGH saturation varies from 13% to 53% in these wells.High saturations of up to 75% were observed sporadically in some of these wells [15,20].The exploration results of GMGS-1 and GMGS-3 showed that the NGH enrichment area in the Shenhu area holds large NGH reserves; however, NGH saturations and permeability of reservoirs show clear differences both in horizontal and vertical directions, and strong heterogeneity [20].Table 1 showed the characteristics of the NGH deposits with greater exploitation potential.

Hydrates in the Shenhu Area
In 2007, NGH samples were successfully drilled from the SH2, SH3, and SH7 sites of GMGS-1 research area in the Shenhu area, in the northern South China Sea (Figure 1) [14].The drilling results showed that the top of the hydrate-bearing layers (HBL) are located at 115-229 m below the ocean floor [14,15].In 2015, a total of 23 wells across 19 sites, were drilled in GMGS-3 research area in the Shenhu area in the northern South China Sea (Figure 1) [15,20,21].The well logging data showed that all of these stations have NGH.Hydrate samples were collected in four of these wells.The hydrate layers are 13-70 m thick.The NGH saturation varies from 13% to 53% in these wells.High saturations of up to 75% were observed sporadically in some of these wells [15,20].The exploration results of GMGS-1 and GMGS-3 showed that the NGH enrichment area in the Shenhu area holds large NGH reserves; however, NGH saturations and permeability of reservoirs show clear differences both in horizontal and vertical directions, and strong heterogeneity [20].Table 1 showed the characteristics of the NGH deposits with greater exploitation potential.NGHs are mainly stored in sandy and silty marine sediments.Clay and clastic limestone and secondary permeability related to structure also host NGH in some areas of the South China Sea [18].The NGH-enriched reservoirs in the Shenhu area have poor permeability, which hinders the transfer of heat and pressure between the production wells and the reservoir, and reduces the conversion efficiency of the continuous dissociation of NGH [14,22,23].The problem of how to efficiently and economically exploit natural gas hydrates in low-permeability marine sediment reservoirs is an important issue.Being able to exploit NGHs from marine sediments, especially in low-permeability silt and clay sediments with a less producible capacity, will greatly expand the potential use of NGH as a gas resource [23].
The reservoir stimulation technique (RST) is to stimulate the reservoir through hydraulic fracturing method [24].RST is aimed at fracturing more main cracks and multi-level secondary cracks to form a fractured network system, creating communication between main cracks and natural cracks at the same time [25].'Breaking up' the effective reservoir to increase the contact area between cracked walls and the reservoir, and reduce the seepage distance from effective reservoir to cracks, would greatly improve the permeability of the reservoir [26,27].RST has expanded to include low-permeability oil and gas shale, as well as tight sand reservoirs [28].The depressurization method and RST combine to exploit NGHs, which will increase the transmission rate of depressurization within NGH deposits and promote NGH decomposition in fractured zones, in addition to being conducive to the discharge of methane [23,[29][30][31].
The formation of NGH in host sediment pore space results in a higher bulk modulus and increased mechanical strength [4,6].In an ideal case, a hydrate deposit probably had a sufficiently brittle response to fracturing; our model would provide a base case with which actual testing can be compared in order to assess the likelihood of artificial fracturing when inducing additional permeability in semi-consolidated marine sediments, which, without NGH, would be expected to respond in a more mechanically-ductile manner [23].

Objective
The main objective of this study is to investigate the influence of RST on the conversion efficiency of NGHs under different accumulation conditions such as initial NGH saturation, intrinsic permeability, and variable-stimulation effect of the depressurization method.We provide a reference program for increasing conversion efficiency of NGH accumulations from the Shenhu area in the northern South China Sea, in different permeability and saturation conditions, especially in low-permeability and high-saturation conditions.

Numerical Simulation Code
The simulator model used in this work was TOUGH+HYDRATE v1.0, a numerical simulator developed by Moridis from the Lawrence Berkeley National Laboratory (Berkeley, CA, USA), which is the first iteration of TOUGH+, and the successor to TOUGH2.The model can simulate the formation and decomposition of natural gas hydrates, 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 CH 4 -hydrate deposits in the subsurface (that is, in permafrost and deep ocean sediments) as well as laboratory experiments of hydrate dissociation and formation in porous and fractured media, using the methods of depressurization, heating injection, and injection inhibition [32].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.).In recent years, TOUGH+HYDRATE v1.0 models have been widely used in NGH simulations.Li et al. [14] used this model to evaluate NGH conversion potential through depressurization and thermal stimulation from the SH7 site.Su et al. [22] used depressurization and thermal stimulation method to analyze NGH conversion efficiency from the SH2 site.Chen et al. [23] used this model to investigate the effect of fracturing technology on the production efficiency of NGHs by depressurization from the SH7 site.

System Parameters and Initialization of the Model
The geologic system used in this study was according to the drilling results of GMGS-1 and GMGS-3 in the Shenhu area, the northern South China Sea.The hydrate samples from this area were dominated by methane hydrate, and in some areas, they were almost pure methane hydrate (99.2%) [14].Therefore, only methane hydrate was simulated.The system parameters, part of the initial conditions and mathematical models of the simulation are shown in Table 2.The main parameters in the simulation were derived from the previous literature on NGH reservoirs in this area [14,20,23].

Production Well Design
In this paper, a single vertical well which had a 6 m height production interval was located in the middle of the NGH deposit with a radius of r w = 0.1 m, and the simulation system was cylindrical (Figure 2A).The production interval design referred to Su [22] and Li's [35] studies.Setting the production well in the middle of the NGH deposit was to reduce natural gas spillage through overburden (OB) or underburden (UB).

Production Well Design
In this paper, a single vertical well which had a 6 m height production interval was located in the middle of the NGH deposit with a radius of rw = 0.1 m, and the simulation system was cylindrical (Figure 2A).The production interval design referred to Su [22] and Li's [35] studies.Setting the production well in the middle of the NGH deposit was to reduce natural gas spillage through overburden (OB) or underburden (UB).

Cracks Stimulation Design
Reservoir stimulation is a complex process.The fracturing cracks morphology were affected by a lot of parameters such as the physical and mechanical properties of the formation, and the stress distribution of target formation [23].Therefore, in this work, in order to simplify the physical model, the fractured network system was simplified into multiple horizontal cracks, with different crack densities to represent the different stimulation effects.The larger crack quantity and smaller crack spacing represents better stimulation effects with a densely fractured network in the fracturing zone.The opposite conditions represent poor stimulation effects, with a sparsely fractured network in the fracturing zone.According to the different stimulation effect, the cracks were divided into spacing categories, Δl = 3, 2, 1 m, for three, four, and five cracks respectively, and these cracks were uniformly

Cracks Stimulation Design
Reservoir stimulation is a complex process.The fracturing cracks morphology were affected by a lot of parameters such as the physical and mechanical properties of the formation, and the stress distribution of target formation [23].Therefore, in this work, in order to simplify the physical model, Energies 2018, 11, 339 6 of 16 the fractured network system was simplified into multiple horizontal cracks, with different crack densities to represent the different stimulation effects.The larger crack quantity and smaller crack spacing represents better stimulation effects with a densely fractured network in the fracturing zone.The opposite conditions represent poor stimulation effects, with a sparsely fractured network in the fracturing zone.According to the different stimulation effect, the cracks were divided into spacing categories, ∆l = 3, 2, 1 m, for three, four, and five cracks respectively, and these cracks were uniformly distributed through the production well, which increased the communication between the production well and the NGH deposit.The parameters of cracks were showed in Table 3.Furthermore, the crack length was L f = 40 m, and crack height was h = 10 mm, as shown in Figure 2B.The permeability of the cracks varied according to the porosity.The porosity and permeability have the following relationships [36][37][38] where k 0 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.

Analysis of the Factors' Significance and Influence Rules on Gas Production Efficiency
In this work, we considered the interaction of various factors and used the whole simulation experiments method (L 27 (3) 13 ) to analyze the significance of intrinsic permeability k, initial NGH saturation S H0 , and the stimulation effect (represented by cracks quantity N) and their interaction on NGH conversion efficiency of NGH deposits in the Shenhu area.Because these factors and their interactions both would affect NGH conversion efficiency, it was inaccurate to use a single variable approach to describe the influence rules of these three factors on NGH conversion efficiency; thus, a whole simulation experiment program which contained the interactions between these factors was necessary.In this work, using CH 4 cumulative volume Q cv to represent NGH conversion efficiency, the depressurization exploitation time was set to five years, because the CH 4 production rates were stabilized at this time.
Table 4 showed the levels of these three key factors, and Table 5 showed the 3 3 whole simulation experiments program and the simulation experiments results.As shown in Table 5, for example, the k × S H0 , k × N and S H0 × N meant the interaction column between intrinsic permeability k and initial NGH saturation S H0 , intrinsic permeability k and crack quantity N, and initial NGH saturation S H0 and crack quantity N, respectively, and the 9th, 10th, 12th, and 13th columns were used for error analysis, which were vacant columns and not written into Table 5 in this work.As shown in Table 5, x was the average value of the same level of each factor.As shown in Table 6, analysis of variance was applied to analyze the significance and sensitivity of these three factors on NGH conversion efficiency.In these two tables, ) where, SS T was total sum of the square, C was the correction parameter, n was the test quantity and n = 27.In Equations ( 4) and ( 5), SS i were the sum of the square of each factor and SS e was the error sum of the square, while the number of repetitions of each factor K i = 9.
Energies 2018, 11, 339 8 of 16 df j=k,SH0,N = x j − 1, x j = 3 (7) where, df T is the total degree of freedom, df e is the degree of freedom for error, df j,m are the degrees of freedom for factors j and m, x j is the level of each factor and MS is mean square.F 0.05(2,8) and F 0.01 (2,8) are derived from the standard F(f 1 ,f 2 ) table.
By comparing F i and F 0.05(2,8) and F 0.01(2,8) , the impact of factor i was significant when F i > F 0.05(2,8) , and had a more significant influence when F i > F 0.01 (2,8) .As shown in Table 6, F k > F SH0 > F 0.05(2,8) > F 0.01(2,8) , which meant intrinsic permeability k and initial NGH saturation S H0 had a more significant influence on CH 4 cumulative volume.Furthermore, when F 0.01(2,8) > F N > F 0.05(2,8) , crack quantity N had a significant influence on CH 4 cumulative volume, which was smaller than that of k and S H0 , and the significance of interaction between k, S H0 , and N were small.The results showed that the impact of intrinsic permeability, initial NGH saturation, and the stimulation effect on NGH conversion efficiency were significant; however, the interaction had no significant effect.In order to identify the influence rules of each factor for NGH conversion efficiency, multiple comparisons of the three factors were used.As shown in Tables 7-9, by comparing x cv of each level, the results showed, for intrinsic permeability k, NGH conversion efficiency was substantially increased while increasing k; however, the growth rate decreased.For initial NGH saturation S H0 , a lower S H0 led to a higher NGH conversion efficiency.For crack quantity N, a densely fractured network had a higher NGH conversion efficiency, but the impact of a sparsely fractured network was slightly less than that of a densely fractured network.By comparing the difference between the maximum and minimum x cv of each level, it was showed that the influence on NGH conversion efficiency was increased by about four times, in comparison with S H0 and k.

Range Analysis Method
As shown in Table 6, the sensitivity of the stimulation effect (represented by crack quantity N) on the cumulative volume of CH 4 was significant.In order to investigate the influence of the stimulation effect on NGH conversion efficiency using range analysis method, we compared the values of CH 4 cumulative volume Q cv under different stimulation conditions as shown in Tables 10 and 11.
Where x cv was the average value of CH 4 cumulative volume Q cv , R cv was the range of CH 4 cumulative volume Q cv , and r cv was the rate of change between maximum x cv and minimum x cv .Value of R cv was calculated by subtracting the minimum x from the maximum x, with the following expression, R = x max − x min (12) r = R/x min (13) As shown in Table 10, r cv was decreased with increasing S H0 , and the r cv for case S H0 = 0.3 was the largest.This was because a higher S H0 had a lower effective permeability, and stimulation effect was more obvious in lower effective permeability condition, which led to a less difference of Q cv .The results showed that the sensitivity of variable stimulation effect on NGH conversion efficiency is significant in low-saturation conditions.As shown in Table 11, r cv for case k = 7.5 mD was much bigger than that for cases k = 40 mD and 75 mD, and r cv for cases k = 40 mD and k = 75 mD were similar.This was because, RST had a greater improvement on effective permeability under low permeability conditions, and there were enough seepage channels for the discharge of methane in higher permeability conditions.The results showed that the sensitivity of the variable stimulation effect on NGH conversion efficiency was significant in low-permeability conditions.
The influence of intrinsic permeability on the sensitivity of stimulation effects on NGH conversion efficiency was bigger than that of initial NGH saturation, in comparison with the r cv of Tables 10 and 11. Figure 3 showed the CH 4 production rate (Q pr ) and the cumulative volume (Q cv ) curves under low-permeability conditions for different values of N. As shown in Q pr curves of Figure 3A, in the early stage of exploitation, the Q pr changed greatly, then tended to stabilize after about 200 days.The Q pr for cases N = 3 and N = 5 were similar and much bigger than that for case N = 4.As shown in Q cv curves, the Q cv for cases N = 3 and N = 5 were similar and much bigger than that for case N = 4, and the Q cv increased by 52.4%, in comparison with cases N = 3, N = 5, and N = 4.Our previous research results indicated that the crack spacing ∆l and crack quantity N are the significant parameters for NGH exploitation, and the effect of crack spacing on gas production efficiency is greater than that of crack quantity [23].Crack spacing ∆l = 3 m is the most favorable crack spacing, and the gas production efficiency is increased with increasing crack quantity N [23].The effective permeability is bigger under low-saturation conditions such that crack spacing ∆l = 3 m had a better improvement effect for gas production efficiency.Therefore, the Q pr for cases N = 3 and N = 5 were similar and much bigger than that for case N = 4.
In Figure 3B, S H0 = 0.4.As shown in Q pr curves of Figure 3B, the differences of Q pr were small under three crack quantity conditions, and the Q pr for case N = 5 was the largest.As shown in Q cv curves, the Q cv for case N = 5 was the largest, and increased by 26.4%, in comparison with cases N = 3 and N = 5.
In Figure 3C, S H0 = 0.5.As shown in Q pr curves of Figure 3C, in the early stage of exploitation, the Q pr for case N = 5 was the biggest, and the Q pr for case N = 4 was bigger than that for case N = 3.As exploitation progressed, the Q pr for case N = 4 decreased and trended to stability, and the Q pr for cases N = 3 and N = 4 were similar.As shown in Q cv curves, the Q cv for case N = 5 was much bigger than that for cases N = 3 and N = 4, and case N = 4 was slightly bigger than case N = 3.This was because densely fractured networks had more cracks in the fracturing zone, with a high-density fractured network, which had a better improvement effect for the permeability of NGH deposits.In comparison with cases N = 5 and N = 3, the Q cv increased by 30.6%.
These results showed that, under low-permeability conditions, the influence of densely fractured networks on NGH conversion efficiency was most significant.However, the influences were similar for densely and sparsely fractured networks in low-permeability and low-saturation cases.
because densely fractured networks had more cracks in the fracturing zone, with a high-density fractured network, which had a better improvement effect for the permeability of NGH deposits.In comparison with cases N = 5 and N = 3, the Qcv increased by 30.6%.
These results showed that, under low-permeability conditions, the influence of densely fractured networks on NGH conversion efficiency was most significant.However, the influences were similar for densely and sparsely fractured networks in low-permeability and low-saturation cases.Figures 3C and 4 showed the CH 4 production rate Q pr and the cumulative volume Q cv curves under high-saturation (S H0 = 0.5) and different N conditions.By comparing the Q pr curves of Figures 3C  and 4A,B, under k = 7.5 mD and 40 mD conditions, the Q pr for case N = 5 was bigger than that for cases N = 3 and 4 during the whole exploitation process, and Q pr for cases N = 3 and 4 were similar.In high-permeability cases, the Q pr and Q cv for case N = 3 were the largest.
In Figure 4A, k = 40 mD.As shown in Q cv curves, the Q cv for case N = 5 was much bigger than that for cases N = 3 and N = 4, and the Q cv for cases N = 3 and N = 4 were similar, and Q cv increased by 14.1%, in comparison with cases N = 4 and N = 5.In Figure 4B, k = 75 mD.As shown in Q cv curves, in the early stage, the Q cv of three cases were similar, and then, the Q cv for case N = 3 was the biggest, and the Q cv for case N = 5 was bigger than that for case N = 4 at later time.The Q cv increased by 21.3%, in comparison with cases N = 3 and N = 4.This was because intrinsic permeability was the most important factor affecting gas production efficiency.Under high-permeability conditions, the decomposition front was close to the production well in the early stage, the effects of crack spacing and crack quantity on gas production efficiency were small in this time.As exploitation progressed, hydrate decomposition area was gradually removed from the production well, the effect of RST on gas production efficiency became significant.Therefore, production rate and cumulative volume were higher for N = 3 compared to N = 4 and N = 5 at later time.However, in low-permeability cases, densely fractured networks had more cracks in the fracturing zone, which had a better improvement effect for the permeability of NGH deposits.
These results showed that, under high-saturation conditions, the influence of densely fractured networks on NGH conversion efficiency was most significant.However, the influence of sparsely fractured networks was better in high-permeability and high-saturation cases.
were higher for N = 3 compared to N = 4 and N = 5 at later time.However, in low-permeability cases, densely fractured networks had more cracks in the fracturing zone, which had a better improvement effect for the permeability of NGH deposits.
These results showed that, under high-saturation conditions, the influence of densely fractured networks on NGH conversion efficiency was most significant.However, the influence of sparsely fractured networks was better in high-permeability and high-saturation cases.

Analysis under Other Accumulation Conditions
Figure 5 showed the CH4 production rate Qpr and the cumulative volume Qcv curves under other accumulation conditions.As showed in Qpr curves of Figure 5, in the early stage of exploitation, the Qpr for case N = 5 was smaller than that for case N = 3.However, as exploitation progressed, the Qpr for case N = 5 reached the stable value earlier, and the stable value for case N = 5 was larger than that for case N = 3.By comparing the Qpr curves of Figure 5, the higher the permeability, the lower the saturation, the earlier the Qpr for case N = 5 exceeded case N = 3.This was because, under higher permeability and lower saturation conditions, there were more seepage channels to facilitate the discharge of methane to the production well, and the NGH dissociation area was closed to the production well in the early stage of exploitation, therefore, the improvement effect of densely fractured networks for NGH deposit permeability was not obvious.As exploitation progressed, the improvement effect of densely fractured networks was increased with increasing the distance between the production well and decomposition front.Additionally, in higher permeability and lower saturation cases, the NGH dissociation rates were faster.

Analysis under Other Accumulation Conditions
Figure 5 showed the CH 4 production rate Q pr and the cumulative volume Q cv curves under other accumulation conditions.As showed in Q pr curves of Figure 5, in the early stage of exploitation, the Q pr for case N = 5 was smaller than that for case N = 3.However, as exploitation progressed, the Q pr for case N = 5 reached the stable value earlier, and the stable value for case N = 5 was larger than that for case N = 3.By comparing the Q pr curves of Figure 5, the higher the permeability, the lower the saturation, the earlier the Q pr for case N = 5 exceeded case N = 3.This was because, under higher permeability and lower saturation conditions, there were more seepage channels to facilitate the discharge of methane to the production well, and the NGH dissociation area was closed to the production well in the early stage of exploitation, therefore, the improvement effect of densely fractured networks for NGH deposit permeability was not obvious.As exploitation progressed, the improvement effect of densely fractured networks was increased with increasing the distance between the production well and decomposition front.Additionally, in higher permeability and lower saturation cases, the NGH dissociation rates were faster.
In Figure 5A, k = 40 mD, S H0 = 0.3.As shown in Q cv curves, the Q cv for cases N = 3 and N = 5 were similar and bigger than that for case N = 4 at five years, and the Q cv increased by about 13.3%, in comparison with cases N = 3, N = 5, and N = 4. Furthermore, as exploitation progressed, the Q cv for cases N = 5 was bigger than that for case N = 3.In Figure 5B, S H0 = 0.4.As shown in Q cv curves, the Q cv for case N = 3 was slightly bigger than that for case N = 5, and the Q cv increased by 33.8%, in comparison with cases N = 4 and N = 3.In Figure 5C, S H0 = 0.3.As shown in Q cv curves, the Q cv for case N = 5 was bigger than that for cases N = 3 and N = 4, and the Q cv increased by 13.8% in comparison with cases N = 5 and N = 4.In Figure 5D, S H0 = 0.4.As shown in Q cv curves, the Q cv for cases N = 3 and 5 were much bigger than that for case N = 4, and the Q cv for case N = 5 was the largest.The Q cv increased by 18.4%, in comparison with cases N = 4 and N = 5.This was because crack spacing ∆l = 3 m is the most favorable crack spacing, and the effect of crack spacing on gas production efficiency is greater than that for crack quantity.Therefore, the production rate for N = 3 was bigger that for N = 4.Under low-saturation cases, the effective permeability was higher than that for high-saturation cases, there were more heat and mass transfer channels in low-saturation conditions.Therefore, a higher effective permeability led to a higher gas production efficiency.
These results showed that the influence of densely fractured networks on NGH conversion efficiency was most significant under these accumulation conditions.
than that for crack quantity.Therefore, the production rate for N = 3 was bigger that for N = 4.Under low-saturation cases, the effective permeability was higher than that for high-saturation cases, there were more heat and mass transfer channels in low-saturation conditions.Therefore, a higher effective permeability led to a higher gas production efficiency.
These results showed that the influence of densely fractured networks on NGH conversion efficiency was most significant under these accumulation conditions.

Conclusions
In this paper, NGH deposits in the Shenhu area, on the northern slope of the South China Sea was simulated using TOUGH+HYDRATE v1.0 via RST and the single vertical well depressurization method.Based on the simulation results, the following conclusions were drawn: (1) Combining RST and the single vertical well depressurization method to exploit NGH deposits under different intrinsic permeability and initial NGH saturation conditions, the sensitivity of stimulation effects on NGH conversion efficiency was significant.Furthermore, the sensitivity of intrinsic permeability was larger than that for initial NGH saturation, and the influence of the interaction between these three factors was not obvious.

Conclusions
In this paper, NGH deposits in the Shenhu area, on the northern slope of the South China Sea was simulated using TOUGH+HYDRATE v1.0 via RST and the single vertical well depressurization method.Based on the simulation results, the following conclusions were drawn: (1) Combining RST and the single vertical well depressurization method to exploit NGH deposits under different intrinsic permeability and initial NGH saturation conditions, the sensitivity of stimulation effects on NGH conversion efficiency was significant.Furthermore, the sensitivity of intrinsic permeability was larger than that for initial NGH saturation, and the influence of the interaction between these three factors was not obvious.(2) For the stimulated NGH deposits, NGH conversion efficiency was substantially increased with increasing intrinsic permeability.However, the growth rate decreased, and a lower NGH saturation led to a higher NGH conversion efficiency.The influence on NGH conversion efficiency was increased by about four times in comparison with initial NGH saturation and intrinsic permeability.(3) The sensitivity of the variable stimulation effect on NGH conversion efficiency decreased with increasing initial NGH saturation and intrinsic permeability, respectively, and the sensitivity was most significant under lower intrinsic permeability condition.The influence of intrinsic permeability on the sensitivity of stimulation effect on NGH conversion efficiency was bigger than that for initial NGH saturation.(4) The stimulation effects required for a higher NGH conversion efficiency were different under different accumulation conditions.For sparsely fractured networks, the influence was significant under higher permeability and saturation conditions.Furthermore, under lower permeability and saturation conditions, the influence between sparsely and densely fractured networks were similar.For other accumulation cases, dense fracture networks had a significant influence.
It should be stressed here the conclusions above are based on purely numerical simulations.With the development of the reservoir stimulation technique, it would probably be applied to exploit marine gas hydrates and improve conversion efficiency.This will greatly expand the potential use of NGH as a gas resource.Of course, it still requires experimental verification, which is currently under consideration.

Figure 1 .
Figure 1.Location of the research area and drilling sites in the Shenhu area, on the northern slope of the South China Sea.

Figure 1 .
Figure 1.Location of the research area and drilling sites in the Shenhu area, on the northern slope of the South China Sea.

Figure 2 .
Figure 2. (A) Design of the production well, and (B) diagram of the horizontal stimulation cracks.

Figure 2 .
Figure 2. (A) Design of the production well, and (B) diagram of the horizontal stimulation cracks.

Table 1 .
The characteristics of the NGH deposits with greater exploitation potential.

Table 2 .
Production trial properties and mathematical models.

Table 3 .
Parameters of cracks.

Table 4 .
The simulation experiments factor level table.

Table 6 .
Analysis of variance of Q cv .

Table 7 .
Multiple comparison of factor k.

Table 8 .
Multiple comparison of factor S H0 .

Table 9 .
Multiple comparison of factor N.

Table 10 .
Range analysis of Q cv for varying N in different S H0 conditions.

Table 11 .
Range analysis of Q cv of variably N in different k conditions.