Full Text:   <812>

Summary:  <257>

CLC number: TV5

On-line Access: 2018-08-03

Received: 2017-08-30

Revision Accepted: 2018-03-22

Crosschecked: 2018-07-09

Cited: 0

Clicked: 1088

Citations:  Bibtex RefMan EndNote GB/T7714


Xiang Sun


-   Go to

Article info.
Open peer comments

Journal of Zhejiang University SCIENCE A 2018 Vol.19 No.8 P.600-623


A coupled thermal–hydraulic–mechanical–chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics

Author(s):  Xiang Sun, Hao Luo, Kenichi Soga

Affiliation(s):  Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China; more

Corresponding email(s):   shyaansun@outlook.com, hl423@cam.ac.uk, soga@berkeley.edu

Key Words:  Hydrate bearing sediments, Coupled thermal–, hydraulic–, mechanical–, chemical (THMC) model, COMSOL, Gas production

Xiang Sun, Hao Luo, Kenichi Soga. A coupled thermal–hydraulic–mechanical–chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics[J]. Journal of Zhejiang University Science A, 2018, 19(8): 600-623.

@article{title="A coupled thermal–hydraulic–mechanical–chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics",
author="Xiang Sun, Hao Luo, Kenichi Soga",
journal="Journal of Zhejiang University Science A",
publisher="Zhejiang University Press & Springer",

%0 Journal Article
%T A coupled thermal–hydraulic–mechanical–chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics
%A Xiang Sun
%A Hao Luo
%A Kenichi Soga
%J Journal of Zhejiang University SCIENCE A
%V 19
%N 8
%P 600-623
%@ 1673-565X
%D 2018
%I Zhejiang University Press & Springer
%DOI 10.1631/jzus.A1700464

T1 - A coupled thermal–hydraulic–mechanical–chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics
A1 - Xiang Sun
A1 - Hao Luo
A1 - Kenichi Soga
J0 - Journal of Zhejiang University Science A
VL - 19
IS - 8
SP - 600
EP - 623
%@ 1673-565X
Y1 - 2018
PB - Zhejiang University Press & Springer
ER -
DOI - 10.1631/jzus.A1700464

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.

This manuscript deals with the interaction of the dissociation of MH and deformation etc. The authors try to compare the results by several other computer codes.

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

目的:水合物沉积物开采过程是一个热-水-力-化多场耦合过程,该过程包含了不同土层间的热对流、压缩引起的局部变形以及胶结结构破坏引起的应力松弛. 不适当的开采会引起出砂、塌孔等破坏问题. 本文旨在建立天然气水合物沉积物多场耦合计算模型,以量化由开采引起的地质灾害风险.
创新点:1. 通过COMSOL Multiphysics实现水合物开采过程多场耦合有限元控制方程的计算; 2. 建立的模型考虑变形-渗流双向全耦合过程.
方法:1. 通过理论推导,给出开采天然气水合物过程模拟的控制方程;采用偏微分方程模块实现除力学之外其他物理场的耦合计算;采用结构力学模块实现变形计算. 2. 通过与试验数据进行比较验证模型的可靠性. 3. 通过对比全耦合模型与半耦合模型,分析双向耦合对水合物开采过程中沉积物物理力学行为的影响.
结论:1. 所建立模型能够精确模拟水合物开采过程中沉积物的物理力学行为. 2. 当考虑压缩对渗流的影响时,由于孔隙率的降低,计算得到的水合物分解速度要小于不考虑该影响时的速度. 3. 由于存在层间对流效应,非均质模型计算得到的水合物分解速度要快于均质模型.


Darkslateblue:Affiliate; Royal Blue:Author; Turquoise:Article


[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 CH4 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 CH4-CO2 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 CO2 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.

Open peer comments: Debate/Discuss/Question/Opinion


Please provide your name, email address and a comment

Journal of Zhejiang University-SCIENCE, 38 Zheda Road, Hangzhou 310027, China
Tel: +86-571-87952783; E-mail: cjzhang@zju.edu.cn
Copyright © 2000 - Journal of Zhejiang University-SCIENCE