Abstract: Methane gas extraction by a deep well installed in methane hydrate bearing sediments (MHBS) found in deep subsea and permafrost regions is a coupled thermal–hydraulic–mechanical–chemical (THMC) process. The key processes include heat convection between layers, local deformation due to compaction, and stress relaxation caused by damage of the bonded structure. As improper production may induce formation compaction, sand production, and wellbore failures, a numerical code is needed to simulate the THMC processes during methane gas production so that geomechanics and production risks can be quantified. In this study, a nonlinear THMC model was implemented in the partial differential equations (PDE) and structural mechanics module of the COMSOL Multiphysics^® finite element code. This paper describes the non-linear coupled governing equations of the mechanical behavior during hydrate dissociation. In particular, it introduces a new thermodynamics-based constitutive model to simulate the mechanical behavior of hydrate bearing sediments. The performance of the newly developed code was examined by comparing the computed results with test data and other simulation results. The differences between fully coupled and semi-coupled models were analyzed. For example, heterogeneous turbidite layers observed in the Nankai Trough were modeled, and behaviors such as heat convection between different layers, shear stress and strain concentration were examined.

基于COMSOL Multiphysics天然气水合物沉积物热-水-力-化多场耦合模型研究

[1]Ahmadi G, Ji C, Smith DH, 2004. Numerical solution for natural gas production from methane hydrate dissociation. Journal of Petroleum Science and Engineering, 41(4):269-285.

[2]Bishop AW, 1959. The principle of effective stress. Teknik Ukeblad, 39:859-863.

[3]Borowski WS, Paul CK, 1997. The gas hydrate detection problem: recognition of shallow-subbottom gas hazards in deep-water areas. Offshore Technology Conference.

[4]COMSOL, 2017. Reference Manual Version 5.3. COMSOL Corporation, Stockholm, Sweden.

[5]Gens A, 1996. Constitutive modelling: application to compacted soils. Proceedings of the 1st International Conference on Unsaturated Soils, p.1179-1200.

[6]Hannegan D, Todd RJ, Pritchard DM, et al., 2004. MPD-uniquely applicable to methane hydrate drilling. SPE/IADC Underbalanced Technology Conference and Exhibition.

[7]Hyodo M, Kajiyama S, Yoshimoto N, et al., 2013a. Triaxial behaviour of methane hydrate bearing sand. Proceedings of the 10th ISOPE Ocean Mining and Gas Hydrates Symposium.

[8]Hyodo M, Li YH, Yoneda J, et al., 2013b. Mechanical behavior of gas-saturated methane hydrate-bearing sediments. Journal of Geophysical Research: Solid Earth, 118(10):5185-5194.

[9]Hyodo M, Li YH, Yoneda J, et al., 2014. Effects of dissociation on the shear strength and deformation behavior of methane hydrate-bearing sediments. Marine and Petroleum Geology, 51:52-62.

[10]Jommi C, di Prisco C, 1994. Un semplice approccio teorico per la modellazione del comportamento meccanico dei terreni granulari parzialmente saturi. Conference Il Ruolo dei Fluidi nei Problemi di Ingegneria Geotecnica, p.167-188 (in Italian).

[11]Kakumoto M, Tenma N, Sakamoto Y, et al., 2011a. Development of the geo-mechanical simulation code COTHMA. Proceedings of the 9th ISOPE Ocean Mining Symposium.

[12]Kakumoto M, Yoneda J, Miyazaki K, et al., 2011b. Preliminary study on well integrity during gas production from methane hydrate sediments. Proceedings of the International Symposium on Earth and Technology.

[13]Kamath VA, 1998. A perspective on gas production from hydrates. JNOC’s Methane Hydrate International Symposium, p.20-22.

[14]Kayen RE, Lee HJ, 1993. Slope stability in regions of sea-floor gas hydrate: Beaufort Sea continental slope. US Geological Survey Bulletin, USA, p.97-103.

[15]Khalili N, Khabbaz MH, 1998. A unique relationship for χ for the determination of the shear strength of unsaturated soils. Geotechnique, 48(5):681-687.

[16]Kim HC, Bishnoi PR, Heidemann RA, et al., 1987. Kinetics of methane hydrate decomposition. Chemical Engineering Science, 42(7):1645-1653.

[17]Kim J, Moridis G, Yang D, et al., 2012. Numerical studies on two-way coupled fluid flow and geomechanics in hydrate deposits. SPE Journal, 17(2):485-501.

[18]Kimoto S, Oka F, Fushita T, et al., 2007. A chemo-thermo-mechanically coupled numerical simulation of the subsurface ground deformations due to methane hydrate dissociation. Computers and Geotechnics, 34(4):216-228.

[19]Klar A, Soga K, Ng MYA, 2010. Coupled deformation–flow analysis for methane hydrate extraction. Geotechnique, 60(10):765-776.

[20]Kowalsky MB, Moridis GJ, 2007. Comparison of kinetic and equilibrium reaction models in simulating gas hydrate behavior in porous media. Energy Conversion and Management, 48(6):1850-1863.

[21]Li YH, Liu WG, Song YC, et al., 2016a. Creep behaviors of methane hydrate coexisting with ice. Journal of Natural Gas Science and Engineering, 33:347-354.

[22]Li YH, Liu WG, Zhu YM, et al., 2016b. Mechanical behaviors of permafrost-associated methane hydrate-bearing sediments under different mining methods. Applied Energy, 162:1627-1632.

[23]Masuda Y, Fujinaga Y, Naganawa S, et al., 1999. Modeling and experimental studies on dissociation of methane gas hydrates in Berea sandstone cores. Proceedings of the 3rd International Gas Hydrate Conference.

[24]Masui A, Haneda H, Ogata Y, et al., 2006. Triaxial compression test on submarine sediment containing methane hydrate in deep sea off the coast off Japan. Proceedings of the 41st Annual Conference, Japanese Geotechnical Society.

[25]Masui A, Haneda H, Ogata Y, et al., 2007. Compaction behavior of Toyoura sand during methane hydrate dissociation. Proceedings of the 7th ISOPE Ocean Mining Symposium.

[26]McIver RD, 1982. Role of naturally occurring gas hydrates in sediment transport. AAPG Bulletin, 66(6):789-792.

[27]Minagawa H, Nishikawa Y, Ikeda I, et al., 2008. Characterization of sand sediment by pore size distribution and permeability using proton nuclear magnetic resonance measurement. Journal of Geophysical Research: Solid Earth, 113(B7):B07210.

[28]Moridis GJ, Collett TS, Dallimore SR, et al., 2004. Numerical studies of gas production from several CH₄ hydrate zones at the Mallik site, Mackenzie Delta, Canada. Journal of Petroleum Science and Engineering, 43(3-4):219-238.

[29]Paull CK, Ussler W, Dallimore SR, et al., 2007. Origin of pingo-like features on the Beaufort Sea shelf and their possible relationship to decomposing methane gas hydrates. Geophysical Research Letters, 34(1):L01603.

[30]Ruppel C, 2011. Methane Hydrates and the Future of Natural Gas. MITEI Natural Gas Report, Supplementary Paper on Methane Hydrates, p.25.

[31]Sakamoto Y, Komai T, Miyazaki K, et al., 2010. Laboratory-scale experiments of the methane hydrate dissociation process in a porous media and numerical study for the estimation of permeability in methane hydrate reservoir. Journal of Thermodynamics, 2010:452326.

[32]Schrefler BA, 1984. The Finite Element Method in Soil Consolidation. MS Thesis, University College of Swansea, Swansea, UK.

[33]Song YC, Yu F, Li YH, et al., 2010. Mechanical property of artificial methane hydrate under triaxial compression. Journal of Natural Gas Chemistry, 19(3):246-250.

[34]Song YC, Wang XJ, Yang MJ, et al., 2013. Study of selected factors affecting hydrate-based carbon dioxide separation from simulated fuel gas in porous media. Energy & Fuels, 27(6):3341-3348.

[35]Song YC, Cheng CX, Zhao JF, et al., 2015. Evaluation of gas production from methane hydrates using depressurization, thermal stimulation and combined methods. Applied Energy, 145:265-277.

[36]Sun X, Guo XX, Shao LT, et al., 2015. A thermodynamics-based critical state constitutive model for methane hydrate bearing sediment. Journal of Natural Gas Science and Engineering, 27:1024-1034.

[37]Taylor DW, 1948. Fundamentals of Soil Mechanics. Chapman and Hall, London, UK.

[38]Uchida S, 2013. Numerical Investigation of Geomechanical Behaviour of Hydrate-bearing Sediments. MS Thesis, University of Cambridge, Cambridge, UK.

[39]Uchida S, Soga K, Yamamoto K, 2012. Critical state soil constitutive model for methane hydrate soil. Journal of Geophysical Research: Solid Earth, 117(B3):B03209.

[40]Uchida S, Xie XG, Leung YF, 2016a. Role of critical state framework in understanding geomechanical behavior of methane hydrate-bearing sediments. Journal of Geophysical Research: Solid Earth, 121(8):5580-5595.

[41]Uchida S, Deusner C, Klar A, et al., 2016b. Thermo-hydro-chemo-mechanical formulation for CH₄-CO₂ hydrate conversion based on hydrate formation and dissociation in hydrate-bearing sediments. Geo-Chicago 2016, p.235-244.

[42]van Genuchten MT, 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44(5):892-898.

[43]Wilder JW, Moridis GJ, Wilson SJ, et al., 2008. An international effort to compare gas hydrate reservoir simulators. Proceedings of the 6th International Conference on Gas Hydrates (ICGH 2008).

[44]Yang MJ, Song YC, Jiang LL, et al., 2014. Hydrate-based technology for CO₂ capture from fossil fuel power plants. Applied Energy, 116:26-40.

[45]Zhao JF, Yao L, Song YC, et al., 2011. In situ observations by magnetic resonance imaging for formation and dissociation of tetrahydrofuran hydrate in porous media. Magnetic Resonance Imaging, 29(2):281-288.

[46]Zhao JF, Yu T, Song YC, et al., 2013. Numerical simulation of gas production from hydrate deposits using a single vertical well by depressurization in the Qilian Mountain permafrost, Qinghai-Tibet Plateau, China. Energy, 52:308-319.

[47]Zhao JF, Liu D, Yang MJ, et al., 2014. Analysis of heat transfer effects on gas production from methane hydrate by depressurization. International Journal of Heat and Mass Transfer, 77:529-541.

[48]Zhao JF, Zhu ZH, Song YC, et al., 2015. Analyzing the process of gas production for natural gas hydrate using depressurization. Applied Energy, 142:125-134.