Cyclic methane hydrate production stimulated with CO 2 and N 2

. The cyclic methane hydrate production method was proposed with CO 2 and N 2 mixture stimulation. The cyclic production model was established based on actual hydrate reservoir parameters, accordingly, the production characteristics were analyzed, and a sensitivity analysis was conducted. The results show the following: (1) The depressurization mechanism is dominant in the cyclic production. CH 4 production and CH 4 hydrate dissociation can be greatly enhanced because the cyclic process can effectively reduce the partial pressure of CH 4 (gas phase). However, there is a limited effect for CO 2 storage. (2) Heat supply is essential for continuous hydrate dissociation. The CH 4 hydrate dissociation degree is the highest in the near-wellbore area; in addition, the ﬂ uid porosity and effective permeability are signi ﬁ cantly improved, and the reservoir temperature is obviously decreased. (3) The initial CH 4 hydrate saturation, absolute permeability, intrinsic CO 2 hydrate formation kinetic constant, injection time and production time can signi ﬁ cantly in ﬂ uence the production performance of the natural gas hydrate reservoir.


Introduction
Natural gas hydrates (mainly methane hydrates) are crystalline solids stabilized under low-temperature and highpressure environments, which have been considered a potential energy supply since the mid-1960s due to their wide distribution and enormous reserves [1]. Therefore, the effective development of hydrate reservoirs (i.e., the methane recovery) is of great significance.
Currently, depressurization, heat stimulation, chemical inhibitor stimulation and gas exchange (CO 2 ) are the most investigated methods for hydrate reservoir development [2]. Among these methods, the depressurization method is acknowledged as the most potential one due to its simplicity and efficiency. Many studies in the depressurization method have been conducted, such as sensitivity analysis of laboratory-scale parameters [3], heat transfer behavior [4], hydraulic fracturing enhanced depressurization [5], decreasing bottom-hole pressure [6], hydrate production trials in the Mallik area [7], Eastern Nankai Trough [8] and South China Sea [9], and relevant simulation studies [10][11][12].
The gas exchange method was initially proposed by Ohgaki et al. [13]. When the environment temperature was lower than 10°C (283 K), the CO 2 hydrate was more stable than the CH 4 hydrate (Fig. 1). The CH 4 hydrate will dissociate, and the CO 2 hydrate will form if the reservoir condition falls into the shadow area in Figure 1. Therefore, CO 2 injection into hydrate reservoirs can enhance CH 4 recovery and sequester CO 2 (the well-known greenhouse gas) under a suitable P-T range. This is the main reason why the gas exchange method has received special attention in recent years. Lee et al. [14] investigated the influence of the CO 2 addition on CH 4 + C 3 H 8 hydrates from thermodynamic, microscopic and kinetic aspects, and found high replacement degree at the high CO 2 pressure. Gharasoo et al. [15] developed a model of CO 2 injection into a gashydrate-filled pressure vessel, and evaluated the effective gas hydrate dissociation rate into gas and water by fitting the experimental results. Shagapov et al. [16] developed the mathematical model of the warm CO 2 injection into a CH 4 -saturated hydrate reservoir and concluded that the CO 2 hydrate formation heat had a great impact on the CH 4 hydrate dissociation rate. Khasanov et al. [17] also developed a similar mathematical model of liquid CO 2 injection into the hydrate reservoir and analyzed four different regimes during the process. All these investigations showed the feasibility of the methane recovery from hydrate reservoirs through CO 2 injection.
However, the effect of pure CO 2 injection into the hydrate reservoir may be quite limited, and experimental studies have shown that the existence of N 2 in the CO 2 stream can improve the sweeping region and recovery efficiency. Koh et al. [18] investigated the fluid flow characteristics and exchange mechanisms of continuous CO 2 /N 2 injection into a one-dimensional high-pressure reactor containing the hydrate, and they found that the exchange efficiency was inversely proportional to the gas injection rate. Yang et al. [19] investigated the dissociation of the CH 4 hydrate stimulated with flue gas flooding through a sand pack model, and they concluded that the flue gas injection was potential for hydrate production and CO 2 sequestration. Li et al. [20] investigated the fracture-filled CH 4 hydrate dissociation by CO 2 and CO 2 /N 2 injection, respectively, and they found that the presence of N 2 could enhance the final accumulative methane recovery ratio and exchange rate. In addition to N 2 , Sun et al. [21] also studied the hydrate dissociation using the CO 2 /H 2 mixture injection, and they observed good CH 4 hydrate recovery ratio and CO 2 storage capability.
Although the gas exchange method is promising, the low production rate or efficiency is the main disadvantage of this method. Merey et al. [22] reviewed the CH 4 -CO 2 exchange studies and pointed out that the combination of depressurization and gas exchange was promising for hydrate reservoir development. Okwananke et al. [23] experimentally studied the enhanced-depressurization for CH 4 hydrate dissociation by injection of compressed air and nitrogen, and they found that the CH 4 hydrate was quickly dissociated and concluded that the gas injection was a potential approach to improve conventional depressurization method for the hydrate reservoir development. Based on extensive experimental and numerical investigations, the Ig_nik Sikumi Field production trial was carried out to assess the CH 4 production potential with CO 2 /N 2 stimulation in 2012. In this production trial, a huff and puff method was used which was a combination of depressurization and gas exchange, and the "bulk exchange" may have occurred [24].
To investigate the possibility of simultaneous CH 4 production and CO 2 sequestration in a field scenario, this paper proposed a cyclic CH 4 hydrate production method stimulated by CO 2 and N 2 mixture based on the production trial mentioned above. An actual-parameter based simulation model was established, the CH 4 hydrate production characteristics were evaluated, and a sensitivity analysis was carried out.

Model description 2.1 Production method
The injected gas in this paper is a mixture of CO 2 and N 2 , and the existence of a N 2 hydrate is ignored because it is difficult to form under typical hydrate reservoir conditions (Fig. 1). The cyclic production scheme [25] is adopted to conduct this research, which is widely used in the heavy oil production by high-temperature steam stimulation in the petroleum industry. The whole production process is composed of several production cycles and each production cycle includes the injection stage, the soaking stage and the production stage. The well is opened in the injection and production stage, and during the soaking stage, the well is shut in in order that the injected gas and pressure disturbance can originate in the depth of the reservoir.

Simulator selection
The methane recovery from the hydrate reservoir is a complicated process involving phase transition, heat and mass transfer, multi-phase flow, etc. Numerical simulation is an effective way to investigate the possible hydrate production behaviors and several codes/simulators have been used to carry out relevant analyses. In this paper, the CMG-STARS [26] was selected to conduct the simulation study. This simulator was used to simulate hydrate reservoir production in the Mount Elbert area [27], and it obtained long-term production characteristics similar to other tested simulators. In addition, it has also been used for sensitivity studies of CH 4 hydrate production [28], economic studies [29], gas production behavior [30], and CO 2 hydrate formation in depleted gas reservoirs [31]. The base model was established using CMG-STARS by referring to the reservoir and operation parameters of the Ig_nik Sikumi Field production trial.

Mass and energy conservation
Assuming that the hydrates in the reservoir are pure CH 4 hydrates, the system is composed of the immobile phase (hydrate phase) and the mobile phases (aqueous and gas phases), and the flows of the gas and aqueous phases follow Darcy's law. The relationships among the components and phases in this study are shown in Table 1. For each single cell, the spatially discretized conservation equations of each component and the energy are shown in equations (1), (2) and (3), respectively.
For the conservation of the flowing component i (i.e., CH 4 , CO 2 , N 2 or H 2 O): where V f is the volume occupied by the mobile phases; V b is the apparent volume of the cell; q A and q G are the density of the aqueous and gas phases, respectively; w i and y i are the mass percent of component i in the aqueous and gas phases, respectively; n f is the number of neighboring cell faces; A l À Á e is the ratio of the effective area and distance between interfaces; k e is the effective permeability at the interface; p A and p G are the pressure of the aqueous and gas phases, respectively; n t is the number of chemical reactions; s 0 ni and s ni are the product and reactant stoichiometric coefficients of component i, respectively; r n is the volumetric reaction rate; and q i is the mass source/sink from the injection/production wells.
For the conservation of the solid component i (i.e., where V v is the volume occupied by the mobile and immobile phases and c i is the volumetric concentration of the solid component i.

For energy conservation:
where V r is the rock volume (solid inert matrix, rock grains); c S is the total solid concentration; U r is the internal energy per rock volume; U A , U G and U S are the internal energy of the aqueous, gas and solid phases, respectively; H A , H G and H rn are the enthalpy of the aqueous phase, gas phases and reaction n, respectively; k e is the effective thermal conductivity at the interface; T is the temperature; and q e is the heat source/sink from the injection/production wells.

Permeability model
When the hydrate exists, the fluid phase permeability in the porous medium will decrease [32]. Currently, there are two methods to treat the permeability variation in the presence of the hydrate [33], namely, the variable absolute permeability method and the constant absolute permeability method. For the first method, the reservoir absolute permeability is modeled as a function of the hydrate saturation; in addition, the mobile phase relative permeability changes with the effective phase saturation. In the second method, the reservoir absolute permeability is constant, while the mobile phase relative permeability changes with the actual phase saturation. We adopted the first method to conduct our research. The definitions of the fluid phase permeability, effective phase saturation and actual phase saturation are shown in equations (4), (5) and (6), respectively: where k b is the effective permeability of phase b; k a is the absolute permeability of the hydrate reservoir; and k rb is the relative permeability of phase b: where S e b is the effective saturation of phase b.
where S b is the actual saturation of phase b.
(1) Absolute permeability model (Carman-Kozeny-type formula) where k a0 is the reservoir absolute permeability without the existence of the gas hydrate; u and u f are the reservoir porosity and fluid porosity, respectively; and m is the model parameter, which is set to 4.3413 by transforming the Civan's permeability-porosity relationship from Anderson et al. [34].
(2) Relative permeability model The flows of the mobile phases follow Darcy's law, and the relative permeability models are shown in equation (8): where k rA and k rG are the relative permeability of the aqueous and gas phases, respectively; S irA and S irG are the irreducible aqueous and gas saturation, respectively; n A is the model parameter, which is set to 5.04 [35]; and n G is the model parameter, which is set to 3.16 [35].

Capillary pressure model
The capillary pressure between the gas phase and the aqueous phase is shown in equation (9): where p c is the capillary pressure; p c0 is the model parameter, which is set to 10 4 Pa [35]; and k is the model parameter, which is set to 0.77437 [35].

CH 4 /CO 2 hydrate dissociation/formation model
Two chemical reactions shown in equations (10) and (11) are considered:  [37,38]. In their model, the dissociation rate was proportional to the hydrate particle surface area or dissociation area and to the difference in the methane fugacity at the equilibrium and dissociation pressures. The fugacity can be approximated with an equivalent pressure when setting the fugacity coefficient equal to 1.0 [39].
We assume that both the dissociation and formation of the CH 4 hydrate and CO 2 hydrate follow the Kim-Bishnoi model. The expression of the CH 4 /CO 2 hydrate dissociation rate can be expressed as follows: [40] dc Hyd dt where c Hyd is the mole quantity of the CH 4 /CO 2 hydrate per unit volume; k 0 d is the intrinsic dissociation rate constant of CH 4 /CO 2 hydrate (Tab. 2); A HS is the reaction specific area, which is 750,000 m 2 /m 3 (i.e., hydrate particles are assumed to be regular spheres with a diameter of 8 lm [41]); q w is the density of the aqueous phase, 1000 kg/m 3 ; q h is the density of the CH 4 hydrate or CO 2 hydrate, 919.7 kg/m 3 or 1100 kg/m 3 [40]; DE d is the activation energy of the dissociation reaction, which is 81 kJ/mol [41] and 102.88 kJ/mol [42] for the CH 4 hydrate and CO 2 hydrate, respectively; R is the gas universal constant; p e is the equilibrium pressure; y is the mole fraction of CH 4 /CO 2 in gas phase; and K is the equilibrium ratio, as follows: where p g is the gas phase pressure; a 1 , a 2 and a 3 are model parameters, and they are calculated based on the experimental results of Adisasmito et al. [43] (Tab. 3), where the average fitting error is 0.46% and 0.39% for the CH 4 hydrate and CO 2 hydrate, respectively.
(2) Hydrate formation rate model Similar to equation (12), the expression of the CH 4 /CO 2 hydrate formation rate (including nucleation and growth) can be expressed as equations (14) and (15).
For the hydrate nucleation process: For the hydrate growth process: where k 0 f is the intrinsic formation rate constant of the CH 4 /CO 2 hydrate (Tab. 2), and DE f is the activation energy of the formation reaction, which is assumed to be equal to that of the dissociation reaction for the CH 4 hydrate or CO 2 hydrate.

Reservoir parameters
As shown in Table 4, the model is established based on the reservoir parameters of "C-1 sand" [24,44]

Operation parameters
The hydrate reservoir is produced by cyclic gas injection with a composition of 22.5 mol% CO 2 and 77.5 mol% N 2 .
The maximum Bottom-Hole Pressure (BHP) was maintained lower than the formation breakdown pressure (10 MPa) for operation safety during the gas injection stage. A total of 5 production cycles (250 d) were conducted. The main operation parameters are shown in Table 5.     4 show the production volume and the production rate of each gas component, respectively. The gas production rate changes periodically, where the highest production rate occurs at the beginning in each cycle, because the diving force (i.e. the production pressure difference) is the largest at this time. A total of 63 250 m 3 (SCM) of gas is produced at the end of the simulation, where CH 4 accounts for 48.78%. The peak gas production rate of each production cycle increases with the increase of the cycle number, where the highest value reaches 8090 m 3 /d. From Table 6, the total production volume of CH 4 remains almost the same except the first cycle, and the percent of the CH 4 in the total gas production accounts for about 50% in each cycle. At a constant BHP of 4.5 MPa, if no gas was injected into the hydrate reservoir, there would be no CH 4 production (the curve labelled with "depressurization (4.5 MPa)" in Fig. 3) when the BHP maintains at 4.5 MPa. This is because it is higher than the corresponding equilibrium pressure of the CH 4 hydrate (4.34 MPa) at the reservoir condition. Besides, the total CH 4 produced by the cyclic production method is 50% more than that by the depressurization method at the BHP of 3 MPa (the curve labelled with "depressurization (3 MPa)" in Fig. 3). This shows the production potential of the cyclic production. Depressurization method obtains the CH 4 production by lowering the reservoir pressure which is economical and simple. However, for low-pressure hydrate reservoirs, the minimum BHP usually cannot be smaller than 2.75 MPa to avoid potential ice formation, and thus the available pressure-reduction space is quite limited, which affects the production. Under this situation, using the cyclic production can enhance CH 4 production. The injected CO 2 and N 2 can cause the partial pressure to change in the gas component and thus affect the stability of the gas hydrate. For the CH 4 hydrate, using the 1st cycle as an example, the partial pressure of CH 4 in the gas phase and the area below/above the corresponding equilibrium pressure at the end of the injection stage are shown in Figures 5 and 6, respectively. From Figure 5, the average radius of pressure propagation into the hydrate reservoir   is about 25 m. The partial pressure of CH 4 in the gas phase changes greatly due to the introduction of substitute gas molecules, and the partial pressure of CH 4 in the nearwellbore area is lower than the equilibrium pressure, leading to the dissociation of the CH 4 hydrate. At the end of the simulation, the CH 4 hydrate saturation in the near-wellbore area decreased significantly (Fig. 7), and the CH 4 hydrate was highly dissociated in the upper part of the HBL, where the injected gas preferentially sweeps due to its smaller gravity compared with water. At the end of production, the dissociation percent [45] of the CH 4 hydrate is approximately 0.759%. The maximum radius of the significantly hydrate-dissociated area is about 20 m in the plane with an average value of 5 m. Therefore, the cyclic gas stimulation method can enhance the dissociation of the CH 4 hydrate, and the depressurization mechanism (partial pressure reduction) is dominant because the CO 2 hydrate formation is quite limited as discussed below.
The partial pressure change in CO 2 can also cause the dissociation and formation of the CO 2 hydrate. At the end of the simulation, a total of 33 556 m 3 (SCM) of gas mixture is injected, where 7550 m 3 is CO 2 . Approximately 3.5% of CO 2 is stored in the reservoir (Tab. 7), and only a relatively small quantity of CO 2 hydrate is formed. The high-concentration CO 2 is mainly in the near-wellbore area, but cyclic production can cause near-wellbore pressure to rise and drop periodically each cycle which can affect CO 2 storage. On the other hand, the partial pressure of CO 2 in the gas phase is mostly lower than the equilibrium value during the cyclic production, and this is not beneficial for  CO 2 hydrate formation. Therefore, the cyclic production method may not effectively store CO 2 , especially as a CO 2 hydrate.

Change of temperature, permeability and fluid porosity
During cyclic production stimulation, the CH 4 hydrate dissociation is dominant compared with the CO 2 hydrate formation, which can be inferred from the temperature distribution (Fig. 8), because the hydrate dissociation is an endothermic process. At the end of the simulation, the reservoir temperature decreases significantly with the maximum radius of 32 m, especially in the near-wellbore area where the lowest temperature approaches 1.5°C. The similarity between the temperature distribution (Fig. 8) and the CH 4 hydrate concentration distribution (Fig. 7) shows that heat supply is quite essential for continuous hydrate dissociation. Besides the latent heat in the HBL, the latent heat contained in the overlying and bottom layers also plays an important part in the CH 4 hydrate dissociation (Fig. 8).
The effective permeability is the key factor affecting the pressure propagation, thus controlling the gas or water flow in the reservoir. It is correlated with the fluid porosity through equation (7), which is determined by the solid phase saturation, i.e. hydrate saturation. From  Figures 9 and 10, the effective permeability and the fluid porosity of the hydrate reservoir increase significantly due to the CH 4 hydrate dissociation, with the highest permeability close to the reservoir absolute permeability (1D) and the highest fluid porosity close to the reservoir porosity (40%). This shows that the seepage condition improves as the production continues, and it is favorable for the water or gas flow in the hydrate reservoir during the injection or production stage.

Sensitivity analysis
The values of reservoir or operation parameters can have significant influence on the production performance of the hydrate reservoir. Using the base model, the sensitivity analyses of several parameters were conducted, including the initial CH 4 hydrate saturation, absolute permeability, intrinsic CO 2 hydrate formation rate constant, injected gas composition and production time.

Initial CH 4 hydrate saturation
The initial CH 4 hydrate saturation value was set to 0.2, 0.4, 0.6 and 0.8. From Figure 11, the CH 4 production volume increases when the initial CH 4 hydrate saturation increases from 0.2 to 0.6, but the CH 4 hydrate dissociation percent and the CO 2 storage percent decrease. The production is dominated by the reserve (CH 4 hydrate saturation), but a high CH 4 hydrate saturation is not helpful for CO 2 storage due to low gas injectivity (Tab. 8). When the saturation value changes to 0.8, the influence of the poor seepage condition is more significant, and the CH 4 production volume and the total amount of injected gas decrease.

Absolute permeability
The absolute permeability was set to 0.01, 0.1, 1 and 10 D, and the results are shown in Figure 12. With the increase in the absolute permeability, the CH 4 production volume and the CH 4 hydrate dissociation percent increase due to the improvement of the seepage condition, while the CO 2 storage percent decreases significantly.
However, when the permeability increases to 10 D, the CH 4 production will decrease. This is because when the absolute permeability value is too high (10 D for example), the injected gas will mainly accumulate in the upper part of Initial CH 4 hydrate saturation CH production CH hydrate dissociation percent CO storage percent Fig. 11. Results of the sensitivity study of the initial CH 4 hydrate saturation.   the hydrate reservoir because of its lower gravity, which limits its effect in enhancing the CH 4 production. Taking the 1st cycle for an example (Fig. 13), the total sweeping volume of the injected gas is more uniform and larger at the end of gas injection when the absolute permeability is 1 D. Thus, the dissociation region of the CH 4 hydrate is larger in the plane but smaller in the vertical area ( Fig. 14) compared with the situation of 1 D (Fig. 7). The overall dissociation region of the CH 4 hydrate is limited when the absolute permeability is 10 D, resulting in a decline in CH 4 production.

Heterogeneity of absolute permeability
To investigate the influence of the vertical heterogeneity of the absolute permeability. Three situations are considered, i.e. scheme 1 (negative rhythm), scheme 2 (homogeneous), and scheme 3 (positive rhythm), as shown in Table 9. For scheme 1, the value of the absolute permeability increases from the top layer to the bottom layer of the HBL, and scheme 3 is just the opposite. The results of the sensitivity studies are shown in Figure 15. From Figure 15, for the CH 4 production, scheme 1 is the best, followed by scheme 2 and scheme 3 in sequence. This is because the injected gas prefers to accumulate in the upper layer of the hydrate reservoir, and high permeability can assist the injected gas to take effect to increase the CH 4 production, as shown in Figure 16. From Figure 16, however, the CH 4 hydrate dissociation area is more uniform in scheme 3 than in scheme 1. That is, the negative rhythm is the most preferred scheme for CH 4 production, but the least preferred one for CO 2 storage.

Intrinsic CO 2 hydrate formation kinetic constant
The parameter values of hydrate kinetics are taken from bulk phase systems in this paper because these values are less accessible in porous media [46][47][48]. Besides, the measurement of the hydrate formation kinetic constant is much more difficult, which also exhibits a weak correlation to the temperature [49,50]; thus, there is great uncertainty. For the reasons presented above, we carried out the sensitivity study on the intrinsic CO 2 hydrate formation kinetic constant.
The intrinsic CO 2 hydrate formation kinetic constant was set to 5, 10, 100 and 1000 k f 0 , where k f 0 is the value in the base model. The high CO 2 hydrate formation kinetic constant is unfavorable for hydrate production but favorable for CO 2 storage (Fig. 17), especially when it is higher than 10 k f 0 . This trend is observed because as the CO 2 hydrate formation kinetic constant increases, the CO 2 hydrate quick formation in the near-wellbore area weakens the CH 4 hydrate dissociation. The CO 2 hydrate formation degree is more significant than the CH 4 hydrate dissociation, and this can be deduced from the temperature distribution of Figure 18.

Injected gas composition
The CO 2 mole fraction in the injected gas was set to 0.3, 0.4, 0.5, 0.6 and 0.7, and the results are shown in Figure 19. The injected gas composition has little effect on the CH 4 hydrate dissociation percent and CO 2 storage percent. This result is observed because in the established model, there is no essential difference in the role of the pressure reduction by CO 2 and N 2 , as there is only a small amount of CO 2 hydrate formed. The CH 4 production volume slightly decreases with the increase of the CO 2 mole fraction in the injected gas, because these small-amount formed CO 2 hydrate may influence the seepage condition. From this aspect, smaller CO 2 mole fraction in the mixture is more beneficial to the CH 4 production, and pure N 2 can have the best production performance.

Injected gas temperature
The injected gas temperature was set to 50, 100, 150 and 200°C, and the results are shown in Figure 20. The CH 4 hydrate dissociation percent and CO 2 storage percent vary little with the injected gas temperature. The CH 4 production volume slightly increases with the increase of the gas temperature. This is because the injected gas conveys a small amount of heat due to its low heat capacity, and the latent heat contained in the reservoir provides the most consumed heat for the CH 4 hydrate dissociation.

Injection time
The injection time of each cycle was set to 5, 10, 15, 20 and 25 days. From Figure 21, a longer injection time can enhance the CH 4 production volume and the CH 4 hydrate dissociation percent. More gas will be injected into the reservoir as the injection time increases, and this can induce more hydrate to dissociate. The CH 4 production volume increases about 164% as the injection time increases from 5 days to 25 days, but the increase of the CH 4 production slows down when the injection time is beyond 15 days. However, longer injection time is unfavorable for CO 2 storage and the CO 2 storage percent decreases from 6.03% to 2.39% as the injection time increases from 5 days to 25 days.

Production time
The production time was set to 1, 2, 3 and 4 years. From Figure 22, a longer production time can enhance the CH 4 production volume and the CH 4 hydrate dissociation percent. After a four-year production, about 2.27% of the CH 4 hydrate has been dissociated compared to the dissociation percent of 1.03% after a one-year production. The CO 2 storage percent will decrease to a stable low level as the production time increases.

Summary and conclusion
This paper proposed the cyclic methane hydrate production method with CO 2 and N 2 mixture stimulation, and investigated the field-scale production using this method based on a numerical model according to actual hydrate reservoir parameters. The following conclusions were obtained: 1. The depressurization mechanism is dominant in cyclic stimulation with the CO 2 and N 2 mixture; in addition, the CH 4 hydrate can be significantly produced due to the effective reduction of the partial pressure of CH 4 in the gas phase, and the CH 4 hydrate in the upper part of the HBL is highly dissociated. However, the effect of CO 2 storage is limited, and a continuous CO 2 injection (i.e., CO 2 flooding) may be better for CO 2 storage. The waste gas produced from the plant (mainly composed of CO 2 and N 2 ) may be directly injected into the hydrate reservoir for production. 2. In the near-wellbore area, the hydrate dissociation degree is relatively high, and the fluid porosity and effective permeability increase significantly during the cyclic stimulation process. Heat supply plays an important role in the cyclic production. The reservoir temperature decreases significantly, and the overlying and bottom layers can supply some heat for CH 4 hydrate dissociation in addition to the HBL. External heat can be added to improve the hydrate reservoir production during the process. 3. Relatively high initial CH 4 hydrate saturation, high reservoir absolute permeability, and long injection or production time are helpful in improving the CH 4 hydrate production; however, they are adverse to CO 2 storage. Too high of an initial CH 4 hydrate saturation or reservoir absolute permeability will cause a decrease in the CH 4 production due to poor seepage conditions. Relatively speaking, the injected gas temperature and injected gas composition have limited effect on the CH 4 hydrate production or CO 2 storage.

Comments
The findings in this paper are based on the model established from publicly available data; however, some parameter values may not be quite accurate, especially in the porous media environment, which may have a significant influence on the results. For example, the hydrate formation kinetics are not easy to determine, and do not follow an Arrhenius temperature dependence. Moreover, the dissolution of CO 2 into the aqueous phase probably has some influence on the simulation results. However, this study can provide some insights into the cyclic CH 4 hydrate production process.