Laboratory Study on Hydrate Production Using a Slow, Multistage Depressurization Strategy

Li, Yanlong; He, Chuqiao; Wu, Nengyou; Chen, Qiang; Liu, Changling; Sun, Zhixue; Jin, Yurong; Meng, Qingguo

doi:https://doi.org/10.1155/2021/4352910

Geofluids

On this page

Abstract Introduction Numerical Simulation Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Multiscale and Multiphysical Approaches to Fluids Flow in Unconventional Reservoirs

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 4352910 | https://doi.org/10.1155/2021/4352910

Laboratory Study on Hydrate Production Using a Slow, Multistage Depressurization Strategy

Yanlong Li,^1,2Chuqiao He,^1,3Nengyou Wu ,^1,2Qiang Chen,^1,2Changling Liu,^1,2Zhixue Sun ,³Yurong Jin,^1,4and Qingguo Meng^1,2

Academic Editor: Keni Zhang

Received31 Jan 2020

Revised04 Apr 2020

Accepted03 Feb 2021

Published24 Feb 2021

Abstract

Optimization of the depressurization pathways plays a crucial role in avoiding potential geohazards while increasing hydrate production efficiency. In this study, methane hydrate was formed in a flexible plastic vessel and then gas production processes were conducted at constant confining pressure and constant confining temperature. The CMG-STARS simulator was applied to match the experimental gas production behavior and to derive the hydrate intrinsic dissociation constant. Secondly, fluid production behavior, pressure-temperature () responses, and hydrate saturation evolution behaviors under different depressurization pathways were analyzed. The results show that integrated gas-water ratio (IGWR) decreases linearly with the increase in depressurizing magnitude in each step, while it rises logarithmically with the increase in the number of steps. Under the same initial average hydrate saturation and the same total pressure-drop magnitude, a slow and multistage depressurization strategy would help to increase the IGWR and avoid severe temperature drop. The pore pressure rebounds logarithmically once the gas production is suspended, and would decrease to the regular level instantaneously once the shut-in operation is ended. We speculate that the shut-in operation could barely affect the IGWR and formation response in the long-term level.

1. Introduction

Natural gas hydrate (NGH) is an ice-like compound formed from water and gas molecules under relatively low-temperature and high-pressure conditions [1] NGH is considered as the most promising alternative fossil fuel due to its great energy potential [2]. Hydrate-related field exploration and exploitation attract attentions from both the international ocean discovery program (IODP) and national plans such as the United States, Japan, China, South Korea, and India [3, 4].

To extract natural gas from hydrate-bearing sediments (HBS), four methods have been proposed based on thermodynamic conditions of NGH, i.e., depressurization [5, 6], thermal stimulation [7], CO₂-CH₄ exchange [8], and chemical injection. These methods aim to promote hydrate in-situ decomposition and then produce gas and water by using the methods that are widely applied in conventional oil-gas industry. Depressurization is considered to be the most efficient and technically feasible method [9, 10]. However, field production test in Mallik 2L-38project (2007-2008) [11], AT1-P project (2013) [10], AT1-P2~AT1-P3 projects (2017) [12], and SHSC-4 project (2017) [13] showed relatively low gas productivity, which is far to reach the commercial requirements.

Depressurization is also restricted by potential geohazards such as severe sand production [14–16], wellbore instability [17], and deformation of the sediment [18, 19]. These potential geohazards are mainly controlled by the mechanical properties of HBS [20, 21]. Both numerical simulation [22] and experiments [23] reveal that sand production rate tends to become severe with the increase in the magnitude of total pressure-drop and depressurizing rate. Hydrate production test can be forced to terminate once the sand production rate exceeds the tolerance of wellbore or platform [10]. Besides, the influence of geological parameters, such as compressibility [24, 25], absolute permeability [26, 27], relative permeability [25, 28], and clayey content [29], on the gas production behavior has been covered both in lab-scale [7, 9] and field-scale [10] studies.

The depressurization parameters, which mainly include pressure-drop magnitude and depressurizing rate, would also affect gas production behaviors [9, 30–32]. Aggressive depressurization would increase gas production rate in a certain extent [5, 33], but causes unpredictable occurrence of hydrate reformation [34]. There are many pathways to reach the same pressure-drop magnitude such as one-step quick depressurization and slow-stepwise depressurization [30]. The number of steps, the pressure-drop magnitude in each step, and the average depressurizing rate could be used to characterize the depressurization pathways under the same total pressure-drop magnitude.

Therefore, balancing the gas productivity and potential geohazards is crucial during depressurization [35, 36]. The influences of depressurization pathways on gas-water production behaviors and formation response need to be clarified. In this paper, we form methane hydrate in a flexible rubber sleeve. Two depressurizing gas production experiments were conducted at constant confining pressure and temperature. A numerical simulation model based on CMG-STARS module was developed to match the measured gas production, as well as pressure evolutionary behaviors. Then, the verified numerical model was employed to investigate the influence of depressurization pathways on gas-water production behaviors, as well as formation pressure-temperature (P‐T) and hydrate saturation responses.

2. Experiments

The hydrate production simulation device is shown in Figure 1. The device consists of (a) a self-developed cylindrical pressure vessel system, (b) a reactant supply system for injecting gas and water into the pressure vessel, (c) a thermostatic chamber to monitor the temperature, and (d) a system controlling and data acquisition system.

The main body of the pressure vessel is made from titanium alloy with maximum experimental pressure of 15 MPa. A flexible rubber sleeve () is installed inside the titanium alloy vessel. The annulus between the rubber sleeve and inner of the titanium alloy vessel is filled with confining oil. Sediment for HBS formation is packed inside the rubber sleeve. Five pressure sensors (~), with measurement accuracy of 0.01 MPa and measuring ranges of ~10 MPa, are installed linearly on the rubber sleeve to test the real-time pore pressure. The distance from to the outlet of the vessel is 50 mm, and the distance between each of the two sensors is 200 mm.

The maximum effective confining pressure (difference between confining pressure and pore pressure) for the rubber sleeve is 8 MPa when it was packed with sediment. During hydrate formation and dissociation, the temperature and pressure of the confining oil remain constant, indicating that a constant temperature and constant pressure boundary condition could be reached during the experiments. Depressurization process is achieved by adjusting the back pressure valves at the outlet.

Quartz sand with diameter of 0.12 mm~0.18 mm was used in the experiments. The average porosity of the sediment is 39.9%. Methane gas with purity of 99.99% and distilled water were used to form hydrate. Detailed sand packing and fluid injection procedures could be found in Jin et al. [37].

We tested the intrinsic permeability of the sediment before hydrate formation. The confining pressure was adjusted to 2 MPa after the sediment was packed into the rubber sleeve. Then, distilled water was injected into the sediment at a rate of 35 ml/min. During the process of water injection, the pressure in the sediment was recorded by a pressure sensor. When the pressure values remained constant, we considered that the water flow in the sediment is a steady-state flow. The pressure difference between and was about 0.11 MPa, providing that the viscosity of water at 300 K is 0.855 mPa·s. The intrinsic permeability of the sediment could be obtained according to Darcy’s law, which is about 606 mD in this manuscript.

The average hydrate saturation in the HBS, which is determined by the volume of gas consumption, was controlled around 30%. The confining pressure was set as 8 MPa, and the temperature was set as 3.5°C while forming HBS. Once the HBS formation process was completed, the confining pressure was adjusted to 6 MPa, while the effective confining pressure was about 2 MPa.

The one-step quick depressurization and two-step multistage depressurization schedules were applied to simulate hydrate production process. In the one-step quick depressurization experiment, the pressure at the outlet was controlled to decrease to the atmospheric pressure within 5 minutes. The pressure-drop magnitude was 4 MPa during gas production. In the two-step depressurization experiment, the pressure at the outlet was firstly controlled at 2 MPa for 4 hours. Then, it was decreased to the atmospheric pressure within 5 minutes.

3. Numerical Modeling and Verification

The TOUGH+HYDRATE (Lawrence Berkeley National Laboratory, USA) [38], MH21-Hydrates (AIST Group, Japan) [39], STOMP-HYD (Pacific Northwest National Laboratory, USA) [39], and CMG-STARS (Computer Modeling Group Ltd., Canada) [39] are widely adopted to simulate hydrate production processes. In these simulators, the hydrate intrinsic dissociation constant is normally taken from bulk hydrate dissociation. However, recent study documented that the hydrate intrinsic dissociation constant varies with the characteristics of porous media [40]. Therefore, numerical simulation should firstly verify the hydrate intrinsic dissociation constant based on experimental results.

The CMG-STARS is primarily developed for thermal recovery. The simulator considers the coupling effect of fluid flow, heat transfer, and fines transport in porous media and wellbore. In simulating hydrate production process, the hydrate is considered as “unmovable heavy oil” phase filled in the porous media [40, 41], while the other part of the pore space is saturated by either water or gas phase. Based on Kim-Bishnoi model, Uddin et al. [42] modified the hydrate dissociation kinetic equation into the following equation: where is the hydrate dissociation rate in gmol/[m³·d]. represents the intrinsic hydrate decay rate in gmol/[m²·kPa·d]. is the specific area of hydrate particles, . and are water saturation and hydrate saturation, respectively. is the porosity of the sediment, . is the activation energy, 76.516 kJ/[mol·K]. is the universal gas constant, 8.3144 J/[mol·K]. is the temperature in . is the water-hydrate-vapour equilibrium pressure in kPa. .

In Equation (1), the mass concentration of water and hydrate during dissociation can be written as , . Equation (1) can be written as follows: where is the modified hydrate intrinsic dissociation constant for CMG-STARS.

To simulate the constant temperature and constant pressure boundary experimental conditions, a numerical model is built with three-dimensional cylindrical shape (Figure 2). The cross-sectional area and length of the model are the same with the rubber sleeve (). The model is divided into 20 grid blocks in the -direction and 25 grid blocks in the -direction. The cross section of the model is divided into 20 grid blocks in the -direction. Each grid block has the same lengths in the -direction (50 mm) and -direction (3 mm), and the same coverage angle of 18 degree in the -direction. In order to investigate physical response inside the formation, a longitudinal section ( in Figure 3) was taken along the -direction. The longitudinal section crosscuts the model in the -direction. There is no mass exchange between the confining oil and the HBS.

(a)

(b)

The initial average hydrate saturation is 30% according to the experiments. The input parameters for numerical simulation are listed in Table 1. Detailed parameters of adjusting and regressing methods can be found in Ajayi et al. [41]. The one-step depressurization experiment is used for matching the modified hydrate intrinsic dissociation constant. The two-step depressurization experiment is used for verifying the adaptability of the numerical model.

Relative permeability of the HBS was calculated from Stone’s equation [43]. where and are the relative permeabilities of water and gas, respectively. is the immobile water saturation, whereas is the immobile gas saturation. , , , and [44].

Comparisons between the experimental and numerical results for one-step quick depressurization and two-step multistage depressurization are shown in Figures 3 and 4, respectively. Solid lines in Figures 3 and 4 represent the numerical results, while the dash lines represent the experimental results.

(a)

(b)

Figures 3 and 4 indicate that the gas production behaviors between numerical and experimental results match favorable with each other. The one-step quick depressurization witnessed only one local peak in transient gas production rate during the whole production process. However, the two-step multistage pressurization procedure saw another local peak in transient gas production rate immediately after the second depressurization step.

Under the same initial average hydrate saturation, the same boundary temperature, the same pore pressure, and the same host sediment, the cumulative gas volume from two-step multistage depressurization would be a little more than that from one-step quick depressurization. However, the duration for gas production process from two-stepwise depressurization is shorter than that from one-step quick depressurization.

The pressure evolutionary behaviors from the experiments showed obvious hysteresis, though the trend is similar with that from numerical simulation. This could be attributed to the reformation of hydrate or icing at the inlet of pressure sensor in the experiment [45]. Once the hydrate particles accumulate at the inlet of the pressure sensor, the pressure sensor would not be able to record the real-time pore pressure. Further depressurization would dissociate hydrate accumulated at the inlet of pressure sensor; hence, a plummet in pressure would follow. The modified hydrate intrinsic dissociation constant is according to the regression.

4. Numerical Simulation and Result Discussion

4.1. Depressurization Pathways

Based on the verified numerical model, field hydrate reservoir temperature and pressure conditions in 2017 offshore methane hydrate production test in the Nankai Trough of Japan are taken as the initial conditions. The average pore pressure and temperature at the middle of the hydrate reservoir are 13 MPa and 13.5°C, respectively [33]. The confining pressure of the hydrate reservoir is set as 15 MPa according to hydrostatic pressure, and the confining temperature is 17°C. Three hydrate saturation conditions (~75%, ~50%, and~10%) would be considered. The size and boundary conditions of the model are the same with the experiments (Table 1).

During depressurization, the driving force for hydrate dissociation comes from the difference between water-hydrate-vapour equilibrium pressure and real-time bottom-hole flow pressure. Therefore, the water-hydrate-vapour equilibrium pressure at the reservoir temperature (10.0 MPa at 13.5°C) is taken as the reference to choose the depressurization pathways. The total pressure drop is fixed to 4 MPa, which is the same with the experiments. Hence, the final pore pressure for all simulation cases is 6 MPa. We defined six depressurization pathways to reach the total pressure-drop magnitude of 4MPa (Figure 5).

In Figure 5, the 1^st pathway represents one-step quick depressurization, while paths 2~6 are multistage depressurization pathways. In the paths 2^nd~4^th, the bottom-hole pressure is decreased to 6 MPa within 3 hours, indicating the same average depressurizing rate of 4 MPa/3 h. The depressurization processes are divided into 2 steps, 4 steps, and 16 steps equidistantly in the paths 2^nd~4^th, respectively. In the 5^th path, the bottom-hole pressure is decreased to 6 MPa within 4 hours, indicating an average depressurizing rate of 4 MPa/4 h. The 6^th path is taken as an example to investigate the influence of shut-in operation during production. The depressurization procedure for the 6^th path is the same with the 2^nd pathway except the shut-in operation (from 1 h to 2 h).

4.2. Cumulative Gas-Water Production

Figure 6 shows the evolutionary behaviors of cumulative gas and water when the initial hydrate saturation is 10% (the initial water saturation is 90%). It could be concluded from Figure 6 that depressurization pathways could dramatically affect the final cumulative gas and water volume, even though the total pressure-drop magnitude is fixed. Compared with the one-step quick depressurization (the 1^st path), stepwise depressurization could increase the final cumulative gas volume and decrease the final cumulative water volume to some extent. This implies that depressurization mode would have some significant influences on formation multiphase seepage field, as well as the residual fluid saturation.

Therefore, we define the ratio of final cumulative gas volume and water volume as the integrated gas-to-water ratio (IGWR, Equation (4)). The IGWR could be used to evaluate the influence of depressurization pathways on fluid production behaviors.

When the total pressure-drop magnitude is fixed, we use the depressurizing magnitude in each step and the number of steps to identify and characterize different depressurization pathways. Under the same total pressure drop, the influence of depressurizing magnitude in each step and the number of steps on IGWR are shown in Figure 7. Generally, the IGWR increases greatly with the increase in initial hydrate saturation. Under the same initial hydrate saturation, the IGWR decreases linearly with the increase in depressurizing magnitude in each step, while increases logarithmically with the increase in the number of steps.

In Figure 7(a), if the linear fitting curves for the IGWR and depressurizing magnitude in each step are extended leftwards to the -axis, we obtain the theoretical maximum IGWRs for stepwise depressurization. The maximum IGWR is mainly affected by the initial hydrate saturation (Figure 8). The theoretical maximum IGWR could be reached only if the depressurizing magnitude in each step becomes minimum, indicating that the depressurization procedure is conducted linearly at a constant depressurizing rate. However, this is almost impossible to be realized in a field operation.

On the other hand, when the number of steps is low, the IGWR increases rapidly with the increase in the number of steps. However, the IGWR keeps almost constant once the number of steps exceeds a certain value. The logarithmic relationships between IGWR and the number of steps indicate that it is unnecessary to maximize the number of steps in fieldwork. In short, stepwise depressurization is important for increasing the IGWR at the same total magnitude of pressure drop. However, it is unnecessary to maximize the number of steps, which would affect the difficulties of field operation.

A comparison between the 5^th path and the 4^th path indicates that a slower average depressurizing rate would improve the IGWR within a certain scope, although the depressurizing magnitude in each step and the number of steps are the same. However, the IGWR for the 2^nd path almost equals to that for the 6^th path, indicating that the shut-in operation could barely affect the IGWR. When the average depressurizing rate is low, the transient hydrate dissociation rate would be relatively low. Hence, gas phase would be separated from the water phase. Under the condition of relatively low seepage rate, the water phase cannot be carried by gas phase. As a result, the IGWR increases with the decrease in the average depressurizing rate.

4.3. Temperature Evolutionary Behavior

The temperature of the HBS during hydrate production is mainly determined by the combination of endothermic hydrate dissociation process, Joule-Thomson effect, and heat supply from the confining formation. The 3^rd path, with an initial hydrate saturation of 75%, is taken as an example to explore temperature field changing behaviors. Figure 9 shows the temperature distribution on the section (Figure 3). Point A in Figure 9 represents the geometric core of the section . The right side of the model (marked by rightwards arrow) represents the outlet of the model.

The low-temperature area initiates at the outlet (production well) and propagates in a nonpiston style towards the far side of the model in the -direction. However, the temperature field shows obvious axisymmetrical evolutionary behaviors in the -direction once the hydrate dissociation process is initiated. We attribute the axisymmetrical characteristics to the cylindrical heat supply from the confining oil. These evolutionary phenomena in temperature field were also observed in other simulation cases.

To further investigate the influence of depressurization mode on temperature, the depressurization pathways in Figure 5 are divided into two categories. The average depressurizing rate for the 2^nd~4^th paths is 4 MPa/3 h, while the numbers of steps and pressure drop in each step are different. Therefore, the 2^nd~4^th paths can be used to determine the influence of depressurization mode on temperature under the condition of the same total pressure-drop magnitude and the same average depressurizing rate. On the other hand, the average depressurizing rates for the 1^st, 4^th, and 5^th paths are 4 MPa/5 min, 4 MPa/3 h, and 4 MPa/4 h, respectively. Hence, they could be set into the same category to explore the influence of overall depressurizing rate.

Temperature of the confining oil is constant and equivalent for all simulation cases. Therefore, temperature differences among different simulation cases are probably caused by the combination of hydrate dissociation and Joule-Thomson effect, which are both endothermic processes. Figure 10 shows that the minimal temperature value decreases with the increase in hydrate saturation. It is easy to understand because higher hydrate saturation implies more heat would be consumed by hydrate dissociation.

(a)

(b)

(c)

(d)

(e)

(f)

Most interestingly, it can be seen from Figures 10(a), 10(c), and 10(e) that a decrease in the overall depressurizing rate could significantly alleviate the temperature drop. This is in consistency with that obtained by Zhao et al. [46]. Furthermore, a slower depressurization mode would postpone the appearance of minimal temperature value to a certain extent.

Under the same total pressure-drop value (4 MPa) and average depressurizing rate (4 MPa/3 h), increase in the number of steps (or decrease in the pressure-drop magnitude in each step) would benefit for the smoothness of the temperature behaviors. The minimal temperature value decreases in a very narrow range with the increase in the number of steps. However, the depressurization modes have little influences on temperature rebounding process once the production duration exceeds the minimal temperature.

In field operation, temperature drop is an unwanted factor during hydrate production because of potential hydrate reformation risk [34]. However, we infer from the results that we could hardly avoid temperature drop by only increasing the number of steps or decreasing the magnitude of pressure drop in each step. The best choice for avoiding temperature drop would be decreasing the average depressurizing rate. A slow, multistage depressurization strategy would be the best choice.

4.4. Hydrate Saturation Evolutionary Behavior

The hydrate saturation field evolutionary behaviors for the 3^rd depressurization pathway when initial hydrate saturation is 75% are shown in Figure 11. Figure 11 suggests that the hydrate dissociation initiates at the edge of the model near the outlet () and then propagates uniformly and randomly (). This would be probably caused by the abundant heat supply from the confining oil.

Hereafter, the hydrate dissociation front extends towards the far side of the model in the -direction. The hydrate dissociation front showed obvious nonpiston characteristics, which corresponds to the temperature field changing behaviors in Figure 9. In the -direction, the hydrate saturation field shows obvious axisymmetric evolution behaviors, indicating that the hydrate dissociation process is mainly controlled by the temperature field.

To determine the influence of depressurization modes on hydrate dissociation processes, we define the hydrate dissociation duration as the time when the hydrate saturation at the core of the far side of model (point B in Figure 11) is decreased to less than 0.1%. We collected the hydrate dissociation duration for all simulation cases under different hydrate saturations. The results are shown in Figure 12.

(a)

(b)

(c)

It could be concluded from Figure 12 that hydrate dissociation duration varies greatly with depressurization mode. The hydrate dissociation duration would increase remarkably with the decrease of average depressurizing rate (paths 1, 4, and 5). However, when the average depressurizing rate is fixed, the depressurization mode could hardly affect the hydrate dissociation duration (paths 2, 3, and 4). However, the shut-in operation (path 6) would postpone the hydrate dissociation duration greatly (compared with path 2).

4.5. Influence of Well Shut-In Operation

During marine hydrate production process, the gas production process might have to be suspended due to failure in down-hole equipment or irresistible factors such as typhoon [13]. Therefore, shut-in operation is unavoidable in fieldwork, especially for a long-term production. The temperature and pressure recovery behaviors at point A (Figure 9) are taken to explore the influence of shut-in operation (Figure 13).

It could be concluded from Figure 13 that the pore pressure would rebound logarithmically once the gas production is suspended. The higher the hydrate saturation, the faster the pressure recovery rate would be. However, at the late stage of shut-in period, the pressure for the lower hydrate saturation (10%) case exceeds the higher hydrate saturation (50% or 75%) cases. This is because the temperature of the sediment for lower hydrate saturation (10%) case is much higher than that when hydrate saturation is 50% and 75%.

When the shut-in operation is ended, the pore pressure would decrease to the regular level (path 2) instantaneously, whereas the temperature decreases logarithmically. This is quite important for field operation because it can be inferred from Figure 13 that the shut-in operation could only affect the pressure-temperature () changing behaviors in a short range. Therefore, we infer from the combination results from Figures 7 and 13 that the shut-in operation could hardly affect the fluid production behaviors (represented by IGWR) and formation response in the long run.

All the simulation cases in this paper were conducted under the constant confining pressure and constant confining temperature conditions at lab scale. The hydrate production processes are compared qualitatively between lab scale and field scale, although the lab-scale results might be influenced by the boundary effect. The conclusions obtained in this paper would have significant importance for field hydrate production operation.

5. Summary and Conclusions

Both experimental and numerical simulations were conducted in the lab-scale to determine the gas-water production behaviors and formation physical responses to depressurization modes. Two experiments were used to verify the numerical model and to get the accurate value of hydrate intrinsic dissociation constant under constant confining temperature and pressure conditions. Based on the verified lab-scale numerical model, we carried out a series of numerical simulations to determine the influences of depressurization modes on cumulative gas-water production behaviors, formation temperature, and hydrate saturation. All simulation cases were conducted under the same total pressure drop (4 MPa). The influences of shut-in operation were also discussed in this paper. This study yields the following results and conclusions. (1)The gas production behaviors match well with those of numerical simulation. However, the pressure evolutionary behaviors from experiments showed obvious hysteresis, although the trend is similar with that from numerical results. We inferred that the inlet of pressure sensor in the experiments could be blocked by secondary hydrate. The modified hydrate intrinsic dissociation constant is (mol/m²·Pa) according to the matching(2)The influences of depressurization mode on multiphase seepage field and the residual fluid (gas or water) saturation field could be characterized by the integrated gas-water ration (IGWR). The IGWR increases greatly with the increase in initial hydrate saturation. Under the same initial hydrate saturation and the same total pressure-drop magnitude, the IGWR decreases linearly with the increase in depressurizing magnitude in each step, while increase in a logarithmic style with the increase in the number of steps. Therefore, a multistep depressurization is vital for increasing the IGWR(3)Both the temperature decrease and hydrate dissociation initiate at the outlet and propagate towards the far side of the model in a nonpiston style in the -direction. In the -direction, the hydrate saturation and temperature evolve axisymmetrically under the combination of endothermic hydrate dissociation process, Joule-Thomson effect, and heat supply from the confining oil. Slow and multistage depressurization strategy would help to alleviate temperature drop; hence, it probably benefits for avoiding potential hydrate reformation(4)The pore pressure rebounds logarithmically if the gas production is suspended. Once the shut-in operation is ended, the pore pressure would decrease to the regular level instantaneously. Furthermore, the shut-in operation could only affect the IGWR and formation response in a short range, rather than in the long-term level. However, the shut-in operation would prolong the hydrate dissociation duration greatly

Data Availability

Data are available on request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the Taishan Scholar Special Experts Project (No. ts201712079), the National Natural Science Foundation of China (No. 41976074), and the National Key Research and Development Plan (No. 2017YFC030 7600). These supports are gratefully acknowledged.

References

E. D. Sloan, “Fundamental principles and applications of natural gas hydrates,” Nature, vol. 426, no. 6964, pp. 353–359, 2003.
View at: Publisher Site | Google Scholar
Z. R. Chong, S. H. B. Yang, P. Babu, P. Linga, and X.-S. Li, “Review of natural gas hydrates as an energy resource: prospects and challenges,” Applied Energy, vol. 162, pp. 1633–1652, 2016.
View at: Publisher Site | Google Scholar
W. Zhang, J. Liang, J. Wei, P. Su, L. Lin, and W. Huang, “Origin of natural gases and associated gas hydrates in the Shenhu area, northern South China Sea: results from the China gas hydrate drilling expeditions,” Journal of Asian Earth Sciences, vol. 183, p. 103953, 2019.
View at: Publisher Site | Google Scholar
W. Zhang, M. Shao, C. Jiang, and Q. Tian, “World progress of drilling and production test of natural gas hydrate,” Marine Geology & Quaternary Geology, vol. 35, no. 8, pp. 1–13, 2018.
View at: Publisher Site | Google Scholar
M. Terzariol, G. Goldsztein, and J. C. Santamarina, “Maximum recoverable gas from hydrate bearing sediments by depressurization,” Energy, vol. 141, pp. 1622–1628, 2017.
View at: Publisher Site | Google Scholar
Z. R. Chong, J. W. R. Moh, Z. Yin, J. Zhao, and P. Linga, “Effect of vertical wellbore incorporation on energy recovery from aqueous rich hydrate sediments,” Applied Energy, vol. 229, pp. 637–647, 2018.
View at: Publisher Site | Google Scholar
B. Wang, H. Dong, Z. Fan et al., “Numerical analysis of microwave stimulation for enhancing energy recovery from depressurized methane hydrate sediments,” Applied Energy, vol. 262, p. 114559, 2020.
View at: Publisher Site | Google Scholar
D. Lim and H. Ro, “Electrical resistivity measurements of methane hydrate during N₂/CO₂ gas exchange,” Energy & Fuels, vol. 31, no. 1, pp. 708–713, 2016.
View at: Publisher Site | Google Scholar
C. Liu, Y. Li, L. Liu, G. Hu, Q. Chen, and N. Wu, “An integrated experimental system for gas hydrate drilling and production and a preliminary experiment of the depressurization method,” Natural Gas Industry B, vol. 7, no. 1, pp. 56–63, 2020.
View at: Publisher Site | Google Scholar
K. Yamamoto, Y. Terao, and T. Fujii, “Operational overview of the first offshore production test of methane hydrates in the Eastern Nankai Trough,” in 2014 Offshore Technology Conference, Houston, Texas, USA, 2014.
View at: Google Scholar
K. Heeschen, E. Spangenberg, J. M. Schicks, M. Priegnitz, R. Giese, and M. LH, Eds.“Simulating the gas hydrate production test at Mallik using the pilot scale pressure reservoir LARS,” in EGU General Assembly Conference Abstracts, K. Heeschen, E. Spangenberg, J. M. Schicks, M. Priegnitz, R. Giese, and M. LH, Eds., Vienna, Austria, 2014.
View at: Google Scholar
K. Yamamoto, X. Wang, M. Tamaki, and K. Suzuki, “Hydrate reformation characteristics in natural gas hydrate dissociation process: a review,” RSC Advances, vol. 256, pp. 113878–116013, 2019.
View at: Publisher Site | Google Scholar
J. F. Li, J. L. Ye, X. W. Qin et al., “The first offshore natural gas hydrate production test in South China Sea,” China Geology, vol. 1, no. 1, pp. 5–16, 2018.
View at: Publisher Site | Google Scholar
L. I. Yanlong, H. U. Gaowei, L. I. Changling et al., “Gravel sizing method for sand control packing in hydrate production test wells,” Petroleum Exploration and Development, vol. 44, no. 6, pp. 1016–1021, 2017.
View at: Publisher Site | Google Scholar
Y. L. Li, N. Y. Wu, F. L. Ning et al., “A sand-production control system for gas production from clayey silt hydrate reservoirs,” China Geology, vol. 2, no. 2, pp. 1–13, 2019.
View at: Publisher Site | Google Scholar
Y. Li, F. Ning, N. Wu et al., “Protocol for sand control screen design of production wells for clayey silt hydrate reservoirs: a case study,” Energy Science & Engineering, vol. 8, no. 5, pp. 1438–1449, 2020.
View at: Publisher Site | Google Scholar
Ş. Merey, “Evaluation of drilling parameters in gas hydrate exploration wells,” Journal of Petroleum Science and Engineering, vol. 172, pp. 855–877, 2019.
View at: Publisher Site | Google Scholar
J. Sun, F. Ning, H. Lei et al., “Wellbore stability analysis during drilling through marine gas hydrate-bearing sediments in Shenhu area: a case study,” Journal of Petroleum Science and Engineering, vol. 170, pp. 345–367, 2018.
View at: Publisher Site | Google Scholar
L. Liu, X. Lu, X. Zhang, C. Liu, and B. Du, “Numerical simulations for analyzing deformation characteristics of hydrate-bearing sediments during depressurization,” Advances in Geo-Energy Research, vol. 1, no. 3, pp. 135–147, 2017.
View at: Google Scholar
X. Sun, L. Wang, H. Luo, Y. Song, and Y. Li, “Numerical modeling for the mechanical behavior of marine gas hydrate-bearing sediments during hydrate production by depressurization,” Journal of Petroleum Science and Engineering, vol. 177, pp. 971–982, 2019.
View at: Publisher Site | Google Scholar
L. Dong, Y. Li, H. Liao et al., “Strength estimation for hydrate-bearing sediments based on triaxial shearing tests,” Journal of Petroleum Science and Engineering, vol. 184, p. 106478, 2020.
View at: Publisher Site | Google Scholar
S. Uchida, A. Klar, and K. Yamamoto, “Sand production model in gas hydrate-bearing sediments,” International Journal of Rock Mechanics and Mining Sciences, vol. 86, pp. 303–316, 2016.
View at: Publisher Site | Google Scholar
L. Haojia, L. Yanlong, L. Changling, D. Changyin, W. Nengyou, and S. Jianye, “Calculation model for critical veolocity of sand movement in decomposed hydrate cemented sediment,” Marine Geology $ Quanternary Geology, vol. 37, no. 5, pp. 166–173, 2017.
View at: Google Scholar
X. Sun, Y. Li, Y. Liu, and Y. Song, “The effects of compressibility of natural gas hydrate-bearing sediments on gas production using depressurization,” Energy, vol. 185, pp. 837–846, 2019.
View at: Publisher Site | Google Scholar
J. Cai, Y. Xia, S. Xu, and H. Tian, “Advances in multiphase seepage characteristics of natural gas hydrate sediments,” Chinese Journal of Theoretical and Applied Mechanics, vol. 52, no. 1, pp. 208–223, 2019.
View at: Google Scholar
X. Sun, T. Luo, L. Wang, H. Wang, Y. Song, and Y. Li, “Numerical simulation of gas recovery from a low-permeability hydrate reservoir by depressurization,” Applied Energy, vol. 250, pp. 7–18, 2019.
View at: Publisher Site | Google Scholar
J. Cai, Y. Xia, C. Lu, H. Bian, and S. Zou, “Creeping microstructure and fractal permeability model of natural gas hydrate reservoir,” Marine and Petroleum Geology, vol. 115, p. 104282, 2020.
View at: Publisher Site | Google Scholar
L. Dong, Y. Li, C. Liu et al., “Mechanical properties of methane hydrate-bearing interlayered sediments,” Journal of Ocean University of China, vol. 18, no. 6, pp. 1344–1350, 2019.
View at: Publisher Site | Google Scholar
Y. Sun, X. Ma, W. Guo, R. Jia, and B. Li, “Numerical simulation of the short- and long-term production behavior of the first offshore gas hydrate production test in the South China Sea,” Journal of Petroleum Science and Engineering, vol. 181, p. 106196, 2019.
View at: Publisher Site | Google Scholar
S. C. Phillips, P. B. Flemings, K. You, D. W. Meyer, and T. Dong, “Investigation of in situ salinity and methane hydrate dissociation in coarse- grained sediments by slow, stepwise depressurization,” Marine and Petroleum Geology, vol. 109, pp. 128–144, 2019.
View at: Publisher Site | Google Scholar
Z.-F. Sun, S. Jia, Q. Yuan, C.-Y. Sun, and G.-J. Chen, “One-dimensional study on gas production characteristics of methane hydrate in clayey sediments using depressurization method,” Fuel, vol. 262, p. 116561, 2020.
View at: Publisher Site | Google Scholar
Y. Li, H. Sun, Q. Meng, C. Liu, Q. Chen, and L. Xing, “2-D electrical resistivity tomography assessment of hydrate formation in sandy sediments,” Natural Gas Industry, vol. 39, no. 10, pp. 133–139, 2017.
View at: Publisher Site | Google Scholar
T. Yu, G. Guan, A. Abudula, A. Yoshida, D. Wang, and Y. Song, “Gas recovery enhancement from methane hydrate reservoir in the Nankai Trough using vertical wells,” Energy, vol. 166, pp. 834–844, 2019.
View at: Publisher Site | Google Scholar
M. Yang, J. Zhao, J. Zheng, and Y. Song, “Hydrate reformation characteristics in natural gas hydrate dissociation process: a review,” Applied Energy, vol. 256, article 113878, 2019.
View at: Publisher Site | Google Scholar
Y. Li, C. Liu, L. Liu, J. Sun, H. Liu, and Q. Meng, “Experimental study on evolution behaviors of triaxial-shearing parameters for hydrate-bearing intermediate fine sediment,” Advances in Geo-Energy Research, vol. 2, no. 1, pp. 43–52, 2018.
View at: Publisher Site | Google Scholar
M. T. Reagan, A. F. Queiruga, and G. J. Moridis, “Simulation of gas production from multilayered hydrate-bearing media with fully coupled flow, thermal, chemical and geomechanical processes using TOUGH+Millstone. Part 3: production simulation results,” Transport in Porous Media, vol. 129, no. 1, pp. 179–202, 2019.
View at: Publisher Site | Google Scholar
Y. Jin, D. Yang, and S. X. Li, “Hydrate dissociation conditioned to depressurization below the quadruple point and salinity addition,” Fuel, vol. 255, p. 115758, 2019.
View at: Publisher Site | Google Scholar
G. J. Moridis, M. B. Kowalsky, and K. Pruess, TOUGH+HYDRATE v1.0 User’s Manual: A Code for the Simulation of System Behavior in Hydrate-Bearing Geologic Media. Report LBNL-00149E, Lawrence Berkeley National Laboratory, Berkeley, CA, 2008.
Computer Modeling Group (CMG) Ltd, STARS User Guide–Advanced Processes and Thermal Reservoir Simulator, Computer Modeling Group (CMG) Ltd, Calgary, Alberta, Canada, 2016.
J. Sun, F. Ning, L. Zhang et al., “Numerical simulation on gas production from hydrate reservoir at the 1st offshore test site in the eastern Nankai Trough,” Journal of Natural Gas Science and Engineering, vol. 30, pp. 64–76, 2016.
View at: Publisher Site | Google Scholar
T. Ajayi, B. J. Anderson, Y. Seol, R. Boswell, and E. M. Myshakin, “Key aspects of numerical analysis of gas hydrate reservoir performance: Alaska North Slope Prudhoe Bay Unit “L-Pad” hydrate accumulation,” Journal of Natural Gas Science and Engineering, vol. 51, pp. 37–43, 2018.
View at: Publisher Site | Google Scholar
M. Uddin, D. Coombe, D. Law, and B. Gunter, “Numerical studies of gas hydrate formation and decomposition in a geological reservoir,” Journal of Energy Resources Technology, vol. 130, no. 3, pp. 0195–0738, 2008.
View at: Publisher Site | Google Scholar
H. L. Stone, “Probability model for estimating three-phase relative permeability,” Journal of Petroleum Technology, vol. 22, no. 2, pp. 214–228, 1970.
View at: Publisher Site | Google Scholar
T. Yu, G. Guan, A. Abudula, A. Yoshida, D. Wang, and Y. Song, “Heat-assisted production strategy for oceanic methane hydrate development in the Nankai Trough, Japan,” Journal of Petroleum Science and Engineering, vol. 174, pp. 649–662, 2019.
View at: Publisher Site | Google Scholar
K. A. Birkedal, C. M. Freeman, G. J. Moridis, and A. Graue, “Numerical predictions of experimentally observed methane hydrate dissociation and reformation in sandstone,” Energy & Fuels, vol. 28, no. 9, pp. 5573–5586, 2014.
View at: Publisher Site | Google Scholar
J. Zhao, Y. Liu, X. Guo et al., “Gas production behavior from hydrate-bearing fine natural sediments through optimized step-wise depressurization,” Applied Energy, vol. 260, p. 114275, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Yanlong Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

607

Downloads

980

Citations