The Effects of the Length and Conductivity of Artiﬁcial Fracture on Gas Production from a Class 3 Hydrate Reservoir

: Natural gas hydrate is considered as a potential energy resource. To develop technologies for the exploitation of natural gas hydrate, several ﬁeld gas production tests have been carried out in permafrost and continental slope sediments. However, the gas production rates in these tests were still limited, and the low permeability of the hydrate-bearing sediments is identiﬁed as one of the crucial factors. Artiﬁcial fracturing is proposed to promote gas production rate by improving reservoir permeability. In this research, numerical studies about the effect of fracture length and ﬂuid conductivity on production performance were carried out on an artiﬁcially fractured Class 3 hydrate reservoir (where the single hydrate zone is surrounded by an overlaying and underlying hydrate-free zone), in which the equivalent conductivity method was applied to depict the artiﬁcial fracture. The results show that artiﬁcial fracture can enhance gas production by offering an extra ﬂuid ﬂow channel for the migration of gas released from hydrate dissociation. The effect of fracture length on production is closely related to the time frame of production, and gas production improvement by enlarging the fracture length is observed after a certain production duration. Through the production process, secondary hydrate formation is absent in the fracture, and the high conductivity in the fracture is maintained. The results indicate that the increase in fracture conductivity has a limited effect on enhancing gas production.


Introduction
Natural gas hydrates, commonly occurring in offshore sediments and terrestrial permafrost with relatively low temperature and high pressure, are non-stoichiometric solid compounds in which small guest gas molecules (mainly methane) are captured in the water cages [1,2]. Although accurate resource estimation of the natural gas hydrate is still challenging, its considerable energy potential has been validated even based on the conservative evaluation results [3][4][5]. The research of natural gas hydrate has entered the stage of development in exploitation technology after several production trials in terrestrial and marine hydrate reservoirs, including Mallik in Canada [6], Alaskan North Slope in USA [7], Eastern Nankai Trough in Japan [8,9], and the South China Sea [10,11]. In the pilot production conducted in Eastern Nankai Trough of Japan in 2017, two production tests were carried out by depressurization method; however, due to the failure of sand management, the first production lasted for 12 days, with 41,000 m 3 of methane gas produced, which is way lower than the requirement for commercial exploitation of natural gas hydrate [9]. From 10 May to 9 July 2017, China performed its first offshore natural gas hydrate production test in the South China Sea using the formation fluid extraction method, achieving the production duration of 60 days, with a total gas production volume of 3.09 × 10 5 m 3 [10] but at a relatively low production rate (~5000 m 3 /day). As to the rough estimation, about 0.5 million m 3 /day is a marginal value for a profitable production of methane gas from marine gas hydrate reservoir [1]. Obviously, the gas production levels reached in the previous production tests are all way lower than the marginal value. There are quite some factors that might be associated with the low production level. One of the main reasons is that the low efficiencies in both depressurization spreading in the hydrate reservoir and gas releasing from reservoir to production well are crucial, attributing to the presence of hydrate can result in a sharp decline in effective permeability of hydratebearing sediments [12,13]. Increasing the permeability of the hydrate-bearing sediments by reservoir stimulation techniques might also benefit promoting the gas production rate, as evinced by its effect in shale gas production.
In consideration of the significant effect of hydraulic fracturing in unconventional oil and gas productions for production enhancement [14], the technique of artificial fracture is proposed for gas hydrate exploitation [15], and its feasibility has been studied theoretically and experimentally [16][17][18]. Ito et al. [19] successfully created a hydraulic fracture in unconsolidated sediments for producing methane from hydrate in lab, and the developed fracture was found to occur along the interface between the sand and mud layers. Ito and Narita [20] demonstrated that the propagation behaviors of fractures in unconsolidated sand induced by tension were governed by the compressive stresses, which were quite different from those in hard rocks. In the experimental study carried out in a triaxial pressure cell, Konno et al. [21] found that the permeability of the methane hydrate-bearing sediments increased after fracturing and would be maintained even after fracture closed. Too et al. [16,17] identified the challenges of hydraulic fracturing in sandy sediments and artificially created fractures in the synthetic sediment made in the lab, illustrating the prospect of applying in the natural hydrate-bearing sediments. During the field trial of gas hydrate production in the South China Sea, micro-fractures were induced to increase the permeability of the reservoir sediments to enhance the production performance [10]. Therefore, creating an artificial fracture in a hydrate reservoir is feasible in terms of field engineering and could be a prospective method to increase the gas production rate in a natural gas hydrate reservoir.
Numerical simulations were carried out with the purpose of evaluating the efficiency of artificial fracture. Chen et al. [22] evaluated the effect of fracturing technology on the production efficiency of hydrate reservoir, showing the gas production rate was greatly enhanced with fracturing, but the effect of fracture length was not investigated. Recently, Feng et al. [23] considered the fractured zone around the production well as a highpermeability elliptic zone (e.g., the semi-major axis is 50 m, and the semi-mini axis is 20 m) and showed that it can enhance the gas production efficiently. However, a penny-shaped fracture is usually induced in the hydrate-bearing sediments for the vertical well [16,17,21], and it is commonly described by fracture length and conductivity [24]. Hence, it is necessary to find out how fracture length and conductivity affect gas production from a fractured hydrate reservoir. Since depressurization is regarded as the most promising way in field production [25,26], it was adopted as the production method in this study.
According to the classification of natural gas hydrate deposits [26,27], Class 3 hydrate deposit is characterized by the single hydrate zone that is surrounded by an overlaying and underlying hydrate-free zone. For Class 2 and Class 1 hydrate accumulations, the hydratebearing layer is underlain by the mobile seawater zone and the two-phase zone involving mobile gas, respectively.
In this work, the gas production from an artificially fractured Class 3 hydrate reservoir through depressurization was investigated. The effects of fracture length and conductivity on production performance were investigated by numerical simulations.

Reservoir Specification
Based on the published geological information of hydrate reservoir in the Elbert area on the Alaska North Slope, a hypothetical fractured hydrate reservoir model was established, and numerical simulation experiments were conducted. In February 2007, four modular dynamic formation tester (MDT) tests were conducted at the methane hydrate-bearing  [28][29][30]. Hence, a conceptual 2D radial model was established based on the information of Mount Elbert C2 hydrate zone, as shown in Figure 1; the dimension of the model is 450 m in the radial direction and 152.5 m in the vertical direction.
Based on the published geological information of hydrate reservoir in th on the Alaska North Slope, a hypothetical fractured hydrate reservoir mod lished, and numerical simulation experiments were conducted. In Februa modular dynamic formation tester (MDT) tests were conducted at the meth bearing units at the Mount Elbert Prospect on the Alaska North Slope. Acc classification of natural gas hydrate deposits [26,27], the hydrate deposits in are examples of Class 3. With the highest likelihood of providing interpretab C2 MDT test was selected by the International Methane Hydrate Simulato parison Study for history matching simulation, and informative parameter mined [28][29][30]. Hence, a conceptual 2D radial model was established based mation of Mount Elbert C2 hydrate zone, as shown in Figure 1; the dimension is 450 m in the radial direction and 152.5 m in the vertical direction. In the radial direction, there are 82 grids in total, which include the inne 0.11 m to represent the wellbore. Loosely speaking, the size of grids except block increases logarithmically in the radial direction. In the vertical direction were considered in the numerical model, with 12.5 m of hydrate-bearing wiched by 70 m of impermeable hydrate-free layers divided into 10 nonunif capture the features of dissociation and reformation of hydrate accurately layer has a finer discretization, divided into 110 grids, and most of the grids of 0.1 m. Time-independent Dirichlet boundary conditions were applied a bottom blocks of the model, while closed boundary conditions with no mass were set in the radial direction.
The properties of the Mount Elbert hydrate reservoir in this work we from the related literature, which is listed in Table 1 [29][30][31][32]. Though the rese acterized by coarse-grained sand sediment with high intrinsic permeability mD (1 mD = 10 −15 m 2 ), it has a relatively cold temperature and low initial hyd permeability of 0.12 mD, which brings forth a few challenges for gas produ employing artificial fracture can be an effective method to enhance produ mance. In the radial direction, there are 82 grids in total, which include the innermost one of 0.11 m to represent the wellbore. Loosely speaking, the size of grids except the wellbore block increases logarithmically in the radial direction. In the vertical direction, three layers were considered in the numerical model, with 12.5 m of hydrate-bearing layer sandwiched by 70 m of impermeable hydrate-free layers divided into 10 nonuniform grids. To capture the features of dissociation and reformation of hydrate accurately, the hydrate layer has a finer discretization, divided into 110 grids, and most of the grids have a height of 0.1 m. Time-independent Dirichlet boundary conditions were applied at the top and bottom blocks of the model, while closed boundary conditions with no mass and heat flow were set in the radial direction.
The properties of the Mount Elbert hydrate reservoir in this work were referred to from the related literature, which is listed in Table 1 [29][30][31][32]. Though the reservoir is characterized by coarse-grained sand sediment with high intrinsic permeability of about 1000 mD (1 mD = 10 −15 m 2 ), it has a relatively cold temperature and low initial hydrate-bearing permeability of 0.12 mD, which brings forth a few challenges for gas production. Hence, employing artificial fracture can be an effective method to enhance production performance. Numerical simulations were conducted with TOUGH + HYDRATE (T + H) reservoir simulator that was developed by Lawrence Berkeley National Laboratory, aiming to model the non-isothermal methane release, phase behavior, and fluids and heat flow under conditions typical of methane hydrate accumulations by solving the coupled equations of mass and heat balance [33]. The reliability of T + H has been validated in the two International Methane Hydrate Simulator Code Comparison Studies [29,34]. T + H can be used to fit the experimental results of methane hydrate formation and dissociation in the lab [35], and it can also be applied in large-scale simulation to evaluate the accumulation [36,37] and exploitation [38,39] of hydrate reservoirs. There are two models of hydrate dissociation and formation in T + H, i.e., the equilibrium model and kinetic model; the equilibrium model was recommended for field-scale hydrate production for the justified results and less computational demands [40], which was why it was adopted in the current simulations. The detailed theoretical background of T + H can be found in the user's manual [33].

Artificial Fracture Model
As the depth of the natural gas hydrate reservoir of Mount Elbert is very shallow, approximately a little deeper than 600 m [28], the artificially induced fractures tend to propagate in the radial direction in the hydrate reservoir, originating in the vertical production well [41]. Ito et al. [19] also showed that fracture appeared at the interface between sand and mud layers in unconsolidated sands. In this paper, only the main fracture in the radial direction was considered, and the secondary fractures and fissures were ignored, as shown in Figure 1. Specifically, the artificial fracture plane was assumed located around the center of the hydrate-bearing layer and with a height of 0.1 m. Moreover, it was assumed that only an aqueous phase exists in the fracture because the hydrate will dissociate instantly during the fracturing operation by changing the temperature and pressure of the fracture. In addition, the flow in the fracture was considered to obey the Darcy's Law during the whole simulation process, which was computed [33].
where F is the mass flux; k is the intrinsic permeability; k r is the relative permeability; µ is the viscosity; p is the pressure; g is the gravitational acceleration vector. Defined as the product of fracture permeability and fracture width, fracture conductivity plays an important role in evaluating fracture efficiency. Under the equivalent conductivity hypothesis [42,43], the conductivity of pseudo fracture grids is the same as the actual fracture that can be given as where F cd is the fracture conductivity; k f is the fracture permeability; w f is the fracture width; k e is the equivalent grid permeability; w grid is the equivalent grid width.
In the base case, the conductivity of fracture was set to be 10 D · cm, and the width of fracture grids was 10 cm; as a result, the equivalent permeability was 1000 mD, which was equivalent to the intrinsic permeability of the reservoir.

Relative Permeability and Capillary Pressure
The relationship between the phase saturation and relative permeability of the fluid phase can be given as [33] where k rA and k rG are the relative permeability indices of water and gas; S A and S G are the saturation values of water and gas; S irA and S irG are the irreducible saturation values of water and gas; n A and n G are the relative permeability indices for water and gas. Capillary pressure is defined as [44] where P cap is the capillary pressure; P 0 is the gas entry pressure; λ is the capillary pressure index; S mxA is the maximum of water saturation. Due to the high porosity and permeability of the fracture, it shows different hydraulic properties, compared with the hydrate-bearing sediments. As a result, different hydraulic parameters were employed for sediments and fracture, which are shown in Table 2 [30,45].

Initial Conditions
The temperature and pressure in the middle of the hydrate-bearing layer were 2.95 • C and 6.78 MPa, respectively. The geothermal gradient of this area was 0.0355 K/m, and the hydrostatic pressure gradient was set to be 0.0098 MPa/m [30]. Moreover, it was assumed the pressure and temperature were distributed evenly along each layer in the radial direction. After initialization processes, thermal and hydrostatical conditions achieved a steady state.

Production Strategy
In this study, a single vertical production well was assumed to penetrate into the whole hydrate layer. The target wellbore pressure was 2.7 MPa (the equilibrium temperature is about 0.48 • C) in the production process, which is slightly above the pressure of the quadruple point of methane hydrate (about 2.6 MPa and 0 • C) to prevent ice formation. The simulations were conducted for 5-year gas production with the depressurization Energies 2021, 14, 7513 6 of 13 method. Two simulation groups were designed to investigate the effect of the fracture length and fracture conductivity on hydrate production performance, as summarized in Table 3. In group A, all the artificial fractures had the same conductivity, although they were different in length. On the other hand, in group B, all the fractures had an identical length but with different conductivities.

Effect of Fracture Length on Production
To investigate the effect of fracture length on production, three cases of fractures with different lengths were investigated, and the results were compared. In all the cases, the fracture conductivity was identical, with the value of 10 D · cm.

Gas Production
The gas production rate is a crucial parameter in gas hydrate production, and its variation with fracture length is shown in Figure 2. The gas production rate increases with time although with some fluctuations for all cases. Generally, the longer the fracture is, the higher is the gas production rate. However, a better return from a longer fracture is not presented at the early production stage, which is different from traditional oil and gas production in that the gas production rate increases with fracture length [46]. To be more specific, the gas production rate with a fracture length of 50 m is higher than that of 100 m and 150 m before the production reaches time point A. After that, the advantage of longer fracture of 100 m and 150 m over 50 m appears. Additionally, it is not until point B that the fracture with the maximum length yields the highest gas production rate. In other words, the artificial fracture can enhance gas production rate, but the advantage of longer fracture can only be observed after a sufficiently long production time. In conclusion, the effect of fracture length on production is closely related to the time frame of production, and the benefits from larger fracture length are not presented at the early production stage.
To find out the mechanism of fracture length on the gas production rate, the physical parameters of the grids that are located 50 m away from the production well in the radial direction and adjacent to the fracture plane were monitored (as depicted in Figure 1), and the evolution of pressure and hydrate saturation is shown in Figure 3. As depicted in Figure 3, the monitored grid in the R = 50 m case has the lowest pressure and the fastest hydrate saturation decreasing rate, and the hydrate dissociated completely before 1 year. During the early production stage, as fracture length increases, the influence range of pressure propagation extends due to the high conductivity of artificial fracture, while the pressure drop decreases, and the amount of gas released from hydrate dissociation also reduces. As a result, the gas production rate of R = 50 m case is higher than the other two cases in the early stage before point A. After that, the gas production rate in case of R = 100 m overpasses the fracture length of 50 m. The same mechanism also applies to the interpretation of the phenomenon that happens before point B. words, the artificial fracture can enhance gas production rate, but the advantage o fracture can only be observed after a sufficiently long production time. In conclus effect of fracture length on production is closely related to the time frame of prod and the benefits from larger fracture length are not presented at the early pro stage. To find out the mechanism of fracture length on the gas production rate, the p parameters of the grids that are located 50 m away from the production well in th direction and adjacent to the fracture plane were monitored (as depicted in Figure the evolution of pressure and hydrate saturation is shown in Figure 3. As depicted ure 3, the monitored grid in the R = 50 m case has the lowest pressure and the hydrate saturation decreasing rate, and the hydrate dissociated completely before During the early production stage, as fracture length increases, the influence r pressure propagation extends due to the high conductivity of artificial fracture, w pressure drop decreases, and the amount of gas released from hydrate dissociati reduces. As a result, the gas production rate of R = 50 m case is higher than th two cases in the early stage before point A. After that, the gas production rate in R = 100 m overpasses the fracture length of 50 m. The same mechanism also ap the interpretation of the phenomenon that happens before point B.  In the production process the evolution of (a) pressure (MPa) and (b) hydrate saturation (fraction) at the monitoring points that are 50 m away from the production well in the radial direction and adjacent to the fracture plane.

Evolution of Physical Property Distribution in Production Process
The evolution of the spatial distribution of hydrate saturation in the gas production process is shown in Figure 4, indicating that hydrate dissociation progresses over time. As illustrated in the first column of Figure 4, i.e., a1, b1, and c1, after gas production for 1 year, a similar dissociation front can be found for all cases due to limited hydrate dissociation. This is the case up to point A in Figure 2, gas production rates with different fracture lengths are similar at this stage. After production lasts 3 years, the hydrate dissociation front progresses beyond the fracture when the fracture length is 50 m, as depicted in a2 in Figure 4. In this case, flow channels for gas production are insufficient in the shortest length case, and the gas production rate is limited, compared with longer fracture cases. Therefore, after point A in Figure 2, the gas production rate in the case with a fracture length of 50 m lags behind the other two cases. Similarly, after 5-year gas production, 150 m long fracture can provide the most flow channels for gas migration and yield the highest gas production rate. Therefore, the effect of fracture length on production is closely related to the time frame of production. This finding implies that the fracture length In the production process the evolution of (a) pressure (MPa) and (b) hydrate saturation (fraction) at the monitoring points that are 50 m away from the production well in the radial direction and adjacent to the fracture plane.

Evolution of Physical Property Distribution in Production Process
The evolution of the spatial distribution of hydrate saturation in the gas production process is shown in Figure 4, indicating that hydrate dissociation progresses over time. As illustrated in the first column of Figure 4, i.e., a1, b1, and c1, after gas production for 1 year, a similar dissociation front can be found for all cases due to limited hydrate dissociation. This is the case up to point A in Figure 2, gas production rates with different fracture lengths are similar at this stage. After production lasts 3 years, the hydrate dissociation front progresses beyond the fracture when the fracture length is 50 m, as depicted in a2 in Figure 4. In this case, flow channels for gas production are insufficient in the shortest length case, and the gas production rate is limited, compared with longer fracture cases. Therefore, after point A in Figure 2, the gas production rate in the case with a fracture length of 50 m lags behind the other two cases. Similarly, after 5-year gas production, 150 m long fracture can provide the most flow channels for gas migration and yield the highest gas production rate. Therefore, the effect of fracture length on production is closely related to the time frame of production. This finding implies that the fracture length should be designed properly according to the expected production time to better balance the production performance and the cost of artificially fracturing operations.
Energies 2021, 14, 7513 9 of 14 fracture and boundary layers can be found. The results obtained from Figure 6 support the aforementioned findings that the highest gas production rate after 5 years of gas production appears when the fracture length is 150 m. As with the fracture, length increases, the scope of influence of pressure and temperature drop both increase, and maximum extra gas flow channels are generated by the artificial fracture in this case.
(1) 1 year (2) 3 years (3) 5 years  Detailed analyses of the spatial distributions of hydrate saturation evolution in the base case are depicted in Figure 5. As shown in the figure, hydrate dissociates around the fracture heterogeneously in the production process, and the hydrate dissociation area expands in both vertical and radial directions. This is because hydrate dissociation around the fracture will increase the permeability of hydrate-bearing layer; thus, new flow channels for gas and water migration are generated. As a result, the dissociation of hydrate is accelerated, and the gas production rate is enhanced. It is noteworthy that through the simulation, almost no secondary hydrate formation in the fracture is observed, and the fracture conductivity is kept high, good for long-term industrial production.
The distributions of physical properties of the hydrate-bearing sediments can offer detailed information to interpret the phenomenon of increase in gas production by artificial fracture. In our simulation, the evolutions of the spatial distributions of pressure, temperature, and gas saturation after 5 years of gas production were monitored, as shown in Figure 6. According to the pressure distributions, shown in all of the three cases in Figure 6a, the pressure decrease propagates laterally within the hydrate layer due to the impermeability of the hydrate-free layers. As the dissociation process of gas hydrate is endothermic, it can be inferred from Figure 6b that hydrate dissociation happens around the perforated interval with the largest temperature drop. Moreover, from the spatial distributions of gas saturation in Figure 6c, remarkable gas accumulation between the artificial fracture and boundary layers can be found. The results obtained from Figure 6 Energies 2021, 14, 7513 9 of 13 support the aforementioned findings that the highest gas production rate after 5 years of gas production appears when the fracture length is 150 m. As with the fracture, length increases, the scope of influence of pressure and temperature drop both increase, and maximum extra gas flow channels are generated by the artificial fracture in this case.

Effect of Fracture Conductivity on Production Performance
To investigate the effect of fracture conductivity on production behavior, the length of fracture was taken as the same value of 100 m. As indicated in group B of Table 3, the conditions of three fracture conductivities were investigated, and the gas production performances were compared accordingly.
The gas production rates and cumulative gas production volumes are shown in Figure 7 for different conductivity values at a fracture length of 100 m. As can be found, the gas production rate increases over time in all of the three cases, although the magnitude of enhancement is somewhat limited. After 5 years of production, the cumulative gas production volumes of three cases are 168 × 10 4 ST m 3 , 184 × 10 4 ST m 3 , and 192 × 10 4 ST m 3 , respectively. To make a quantitative comparison, the gas production volume increment from the conductivity of 5 D·cm to 10 D·cm is 16 × 10 4 ST m 3 , corresponding to 9.5% of enhancement. When the conductivity is doubled to 20 D·cm, the growth rate is only 4.3%, with 8 × 10 4 ST m 3 increased. In other words, the fracture with higher fracture conductivity can reduce the fluid flow resistance, then hydrate dissociation can be accelerated by lowering the pressure in fracture. However, this enhancement effect will be compromised because of the reservoir hydrate-bearing permeability improvement induced by hydrate dissociation around the fracture. Based on the simulation results, increasing the fracture conductivity may have limited influence on improving gas production performance. Therefore, in practical gas hydrate production, fracture conductivity should be designed to better balance the cost for increasing the fracture conductivity and the return from gas production increase.

Effect of Fracture Conductivity on Production Performance
To investigate the effect of fracture conductivity on production behavior, the length of fracture was taken as the same value of 100 m. As indicated in group B of Table 3, the conditions of three fracture conductivities were investigated, and the gas production performances were compared accordingly.
The gas production rates and cumulative gas production volumes are shown in Figure 7 for different conductivity values at a fracture length of 100 m. As can be found, the gas production rate increases over time in all of the three cases, although the magnitude of enhancement is somewhat limited. After 5 years of production, the cumulative gas production volumes of three cases are 168 × 10 4 ST m 3 , 184 × 10 4 ST m 3 , and 192 × 10 4 ST m 3 , respectively. To make a quantitative comparison, the gas production volume increment from the conductivity of 5 D cm ⋅ to 10 D cm ⋅ is 16 × 10 4 ST m 3 , corresponding to 9.5% of enhancement. When the conductivity is doubled to 20 D cm ⋅ , the growth rate is only 4.3%, with 8 × 10 4 ST m 3 increased. In other words, the fracture with higher fracture conductivity can reduce the fluid flow resistance, then hydrate dissociation can be accelerated by lowering the pressure in fracture. However, this enhancement effect will be compromised because of the reservoir hydrate-bearing permeability improvement induced by hydrate dissociation around the fracture. Based on the simulation results, increasing the fracture conductivity may have limited influence on improving gas production performance. Therefore, in practical gas hydrate production, fracture conductivity should be designed to better balance the cost for increasing the fracture conductivity and the return from gas production increase. Industrial exploitation of natural gas hydrate requires long-term and efficient production, some factors related to fracturing should be considered. For example, as the reservoir pressure drops with gas production, the effective pressure may increase to induce fracture closure [47]. Moreover, the deformation of sediment induced by hydrate dissociation [48] can result in a decrease in fracture conductivity. In such conditions, how to Industrial exploitation of natural gas hydrate requires long-term and efficient production, some factors related to fracturing should be considered. For example, as the reservoir pressure drops with gas production, the effective pressure may increase to induce fracture closure [47]. Moreover, the deformation of sediment induced by hydrate dissociation [48] can result in a decrease in fracture conductivity. In such conditions, how to maintain fracture conductivity is challenging in long-term industrial production. As a result, further understanding of the interaction between the artificial fracture and the gas production process should be investigated through numerical and practical approaches. Furthermore, horizontal wells [49] and multi-well [50] configurations can be employed to enhance gas production rate to achieve commercial exploitation of hydrate reservoir.

Conclusions
To study the effect of fracture length and conductivity on gas production performances from an artificially fractured hydrate reservoir, a hypothetical model of fractured Class 3 hydrate reservoir was established, and numerical simulations were conducted. The following conclusions can be made from the simulation results obtained: (1) An artificial fracture can enhance the gas production rate, and the enhancement effect will be more pronounced with increasing fracture length. However, the effect of fracture length on production is closely related to the time frame of production, and the benefits from larger fracture length are not presented at the early production stage. (2) The enhancement effect of artificial fracture results from the extra flow channels provided by the fractures. The fracture acts as a flow channel, and the hydrate dissociation is accelerated with increased permeability of hydrate-bearing sediments in the dissociation area. (3) New flow channels for gas and water migration are generated in the area adjacent to artificial fracture, which helps gas migrate into the wellbore. (4) It is found that the secondary hydrate formation is absent in the artificial fracture, and its high conductivity will not be compromised during the entire simulation period. (5) The effect of the artificial fracture conductivity is limited in enhancing the gas production rate. As a result, fracture length and conductivity should be designed properly in practical production projects to balance cost and benefit.

Conflicts of Interest:
The authors declare no conflict of interest.