基于OpenGeoSys的天然气水合物降压开采热-水-力耦合分析模型

doi:10.16285/j.rsm.2023.1012

岩土力学 ›› 2023, Vol. 44 ›› Issue (11): 3191-3202.doi: 10.16285/j.rsm.2023.1012

基于OpenGeoSys的天然气水合物降压开采热-水-力耦合分析模型

叶智刚^{1, 2, 3}，王路君^{1, 2, 3}，朱斌^{1, 2, 3}，陈云敏^{1, 2, 3}

1. 浙江大学软弱土与环境土工教育部重点实验室，浙江杭州 310058；2. 浙江大学岩土工程研究所，浙江杭州 310058； 3. 浙江大学超重力研究中心，浙江杭州 310058

收稿日期:2023-07-12 接受日期:2023-09-11 出版日期:2023-11-28 发布日期:2023-11-28
通讯作者: 王路君，男，1985年生，博士，副教授，主要从事水合物沉积物力学特性和土体多相多场耦合方面研究。E-mail: lujunwang@zju.edu.cn E-mail:zhigangye@zju.edu.cn
作者简介:叶智刚，男，1992年生，博士，博士后，主要从事海洋岩土工程多相多场相互作用研究。
基金资助:
中央高校基本科研业务费专项资金（No. 226-2023-00083）；国家自然科学基金（No. 51988101, No. 52078458）。

Thermo-hydro-mechanical coupling model for natural gas hydrate-bearing sediments with depressurization based on OpenGeoSys

YE Zhi-gang^{1, 2, 3}, WANG Lu-jun^{1, 2, 3}, ZHU Bin^{1, 2, 3}, CHEN Yun-min^{1, 2, 3}

1. Key Laboratory of Soft Soils and Geoenvironmental Engineering of the Ministry of Education, Zhejiang University, Hangzhou, Zhejiang 310058, China; 2. Institute of Geotechnical Engineering, Zhejiang University, Hangzhou, Zhejiang 310058, China; 3. Center for Hypergravity Experimental and Interdisciplinary Research, Zhejiang University, Hangzhou, Zhejiang 310058, China

Received:2023-07-12 Accepted:2023-09-11 Online:2023-11-28 Published:2023-11-28
Supported by:
This work was supported by the Fundamental Research Funds for the Central Universities (226-2023-00083) and the National Natural Science Foundation (51988101, 52078458).

摘要/Abstract

摘要： 天然气水合物是突破世界能源瓶颈的重要潜在途径，降压法是目前公认最经济有效的开采方法。降压开采时会引起固相水合物分解、液气产生与运移、热量传递及沉积物骨架变形，表现出温度场、液气渗流场、应力应变场间强烈的热-水-力（T-H-M）耦合响应，易诱发开采井倾斜、储层塌陷等灾变。基于开源有限元平台OpenGeoSys及其已有的饱和土T-H-M耦合模型，通过嵌入水合物动力分解方程表征固相水合物分解为气液相过程，引入水气质量守恒描述液气相变特性；在此基础上采用非线性互补方法（nonlinear complementary problem，简称NCP）结合特殊主变量选取严格框定液气饱和度范围，并将水合物动力分解速率与控制方程源汇项予以关联；构建了开源的水合物开采T-H-M耦合分析模型。通过与室内降压分解试验和大尺度数值模型对比，验证了该模型的适用性和正确性，并探讨了水合物开采过程中孔隙水压缩性对沉积物T-H-M耦合过程的影响。结果表明：本模型求解水合物分解引起的固液气相变与液气相生成或消失等强非线性问题时稳定性好；降压开采过程中水合物由近及远逐渐分解，产生气液逐渐向边界处运移并逐渐达到稳定运移状态，气、液相饱和度渐趋稳定；受气相溶解的影响，气相饱和度发展滞后于液相饱和度；高孔隙压力下孔隙水压缩性对远端孔压与骨架变形影响显著。

关键词: 水合物沉积物, 降压法, 分析模型, 热水力耦合, OpenGeoSys

Abstract: Natural gas hydrates are considered as one of the most important potential alternatives for energy shortage, and depressurization is believed to be the most economic effectiveness for hydrate exploitation. The exploitation induces solid hydrate dissociation, liquid / gas generation and migration, heat transfer and skeleton deformation, showing strong thermo-hydro-mechanical (THM) coupling behavior of sediments. These processes will trigger geotechnical hazards, such as borehole inclination and reservoir collapse. A THM coupling numerical model for hydrate exploitation was developed based on the open-source FEM platform OpenGeoSys. In detail, the kinetic reaction equation was implemented to represent the liquid / gas generation during gas hydrate dissociation. Air and water mass balance equations were introduced to illustrate the phase change between liquid and gas phases. Besides, the nonlinear complementary problem (NCP) combined with the choice of special primary variables was adopted to strictly constrain the liquid / gas saturation, and kinetic reaction rate was connected to the source / sink terms of governing equations. The model was validated and verified through a test and a large-scale numerical model, and the effects of pore water compressibility on the THM coupling response of hydrate-bearing sediments were discussed. Results show that the model was numerically stable in dealing with the nonlinear problems of phase change among solid, liquid and gas phases, and of phase appearance / disappearance of pore fluids caused by hydrate dissociation; the hydrates dissociate from near to far regions, and the produced gas / liquid gradually migrated to the driving well and reached a steady state, where the gas / liquid saturation tended to be stable; due to the influence of air dissolution, the evolution of gas phase saturation lagged behind that of the liquid saturation; the pore water compressibility had a non-negligible effect on the dissipation of pore pressure and skeleton deformation under large pore pressure.

Key words: hydrate-bearing sediments, depressurization, numerical model, thermo-hydro-mechanical coupling, OpenGeoSys

中图分类号:

TE53

叶智刚, 王路君, 朱斌, 陈云敏, . 基于OpenGeoSys的天然气水合物降压开采热-水-力耦合分析模型[J]. 岩土力学, 2023, 44(11): 3191-3202.

YE Zhi-gang, WANG Lu-jun, ZHU Bin, CHEN Yun-min, . Thermo-hydro-mechanical coupling model for natural gas hydrate-bearing sediments with depressurization based on OpenGeoSys[J]. Rock and Soil Mechanics, 2023, 44(11): 3191-3202.

参考文献

相关文章 15

[1]	刘燕晶, 王路君, 朱斌, 陈云敏, . 考虑填充和黏结作用的含水合物沉积物弹塑性本构模型[J]. 岩土力学, 2022, 43(9): 2471-2482.
[2]	叶智刚, 王路君, 朱斌, 黄家晟, 徐文杰, 陈云敏, . 考虑热渗效应的高温管道−饱和地基相互作用研究[J]. 岩土力学, 2021, 42(3): 691-699.
[3]	王辉, 周世琛, 周博, 薛世峰, 林英松, 吴海明, . 不同破坏准则下水合物沉积物的统计损伤模型[J]. 岩土力学, 2020, 41(12): 4015-4026.
[4]	张小玲, 夏飞, 杜修力, 许成顺, . 考虑含水合物沉积物损伤的多场耦合模型研究[J]. 岩土力学, 2019, 40(11): 4229-4239.
[5]	颜荣涛, 张炳晖, 杨德欢, 李杨, 陈星欣, 韦昌富, . 不同温-压条件下含水合物沉积物的损伤本构关系[J]. 岩土力学, 2018, 39(12): 4421-4428.
[6]	吴起, 卢静生, 李栋梁, 梁德青, . 降压开采过程中含水合物沉积物的力学特性研究[J]. 岩土力学, 2018, 39(12): 4508-4516.
[7]	颜荣涛，梁维云，韦昌富，吴二林, . 考虑赋存模式影响的含水合物沉积物的本构模型研究[J]. , 2017, 38(1): 10-18.
[8]	张强勇，王保群，向文,. 盐岩地下储气库风险评价层次分析模型及应用[J]. , 2014, 35(8): 2299-2306.
[9]	杨期君，赵春风 , . 含气水合物沉积物弹塑性损伤本构模型探讨[J]. , 2014, 35(4): 991-997.
[10]	杨期君，赵春风 , . 水合物沉积物力学性质的三维离散元分析[J]. , 2014, 35(1): 255-262.
[11]	关进安，梁德青，吴能友，万丽华，李栋梁，苏正. 高压环境下流体非饱和渗流对含水合物沉积物的影响[J]. , 2013, 34(3): 639-644.
[12]	吴二林，韦昌富，魏厚振，颜荣涛. 含天然气水合物沉积物损伤统计本构模型[J]. , 2013, 34(1): 60-65.
[13]	李忠友，刘元雪，刘树林，谭仪忠，葛增超 . 火灾作用下隧道衬砌结构变形理论分析模型[J]. , 2012, 33(S2): 307-310.
[14]	吴二林，魏厚振，颜荣涛，韦昌富. 含天然气水合物沉积物分解过程的有限元模拟[J]. , 2012, 33(9): 2811-2821.
[15]	张旭辉，鲁晓兵，王淑云，李清平. 四氢呋喃水合物沉积物静动力学性质试验研究[J]. , 2011, 32(S1): 303-0308.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

基于OpenGeoSys的天然气水合物降压开采热-水-力耦合分析模型

Thermo-hydro-mechanical coupling model for natural gas hydrate-bearing sediments with depressurization based on OpenGeoSys

PDF (Chinese)

PDF (English)

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

Metrics

本文评价

推荐阅读 10