特种油气藏 ›› 2024, Vol. 31 ›› Issue (2): 112-119.DOI: 10.3969/j.issn.1006-6535.2024.02.013

• 油藏工程 • 上一篇    下一篇

天然气水合物注热分解诱发储层变形破坏的正交数值模拟研究

翟诚1,2, 吴迪3, 秦冬冬1,2   

  1. 1.山西工程技术学院,山西 阳泉 045000;
    2.山西工程技术学院矿区生态修复与固废资源化省市共建山西省重点实验室培育基地,山西 阳泉 045000;
    3.辽宁工程技术大学,辽宁 阜新 123000
  • 收稿日期:2022-10-24 修回日期:2024-01-20 出版日期:2024-04-25 发布日期:2024-07-26
  • 通讯作者: 吴迪(1981—),男,教授,2003年毕业于沈阳农业大学水利水电工程专业,2013年毕业于辽宁工程技术大学工程力学专业,获博士学位,现主要从事渗流力学及应用方面的教学与科研工作。
  • 作者简介:翟诚(1981—),男,副教授,2005年毕业于辽宁工程技术大学理论与应用力学专业,2018年毕业于该校工程力学专业,获博士学位,现主要从事天然气水合物开采方面的研究工作。
  • 基金资助:
    国家自然科学基金“压裂页岩吸附超临界CO2驱替CH4渗流与损伤耦合机理研究”(51974147)、“新疆缓斜巨厚煤层季节性地表水保护开采研究”(52104140);山西工程技术学院教师科研启动基金“THM耦合作用下天然气水合物分解诱发储层变形破坏的机理研究”(2021QD-02);2021年第二批博士毕业生、博士后研究人员来晋工作奖励经费项目“THM耦合作用下天然气水合物分解诱发储层变形破坏的机理研究”(2021PT-13)

Qrthogonal Numerical Simulation Study of the Reservoir Deformation Damage Induced by the Thermal Injection Decomposition of Natural Gas Hydrate

Zhai Cheng1,2, Wu Di3, Qin Dongdong1,2   

  1. 1. Shanxi Institute of Technology, Yangquan, Shanxi 045000, China;
    2. The Provincial and Municipal Key Laboratory Cultivation Base of Shanxi Province for Ecological Rehabilitation of Mining Areas and Solid Waste Resource Utilization, Shanxi Institute of Technology, Yangquan, Shanxi 045000, China;
    3. Liaoning Technical University, Fuxin, Liaoning 123000, China
  • Received:2022-10-24 Revised:2024-01-20 Online:2024-04-25 Published:2024-07-26

摘要: 针对天然气水合物热分解而引起的储层失稳破坏机理及影响因素不明确的问题,考虑因水合物热分解引起的地层传热、孔隙流体渗流和固体骨架的变形破坏,建立了水合物沉积层的热-流-固耦合模型及其作用下储层固体骨架的弹塑性本构模型。基于ABAQUS软件的USDFLD子程序进行二次开发,开展正交实验数值模拟,分析注热温度差、水合物沉积层的绝对渗透率和有效主应力差对水合物沉积层近井储层变形破坏影响的敏感性及影响机理。结果表明:对近井储层变形破坏影响的敏感程度由大到小依次为有效主应力差、绝对渗透率、注热温度差;水合物热分解引起的近井储层力学性质劣化和有效应力下降是导致其发生变形破坏的主要原因;有效主应力差越大,等效塑性应变最大值越大,并且受水平地应力非均匀性的影响,近井储层变形破坏最严重的位置始终位于井眼最小水平地应力方向上。该研究成果可为热激法开采天然气水合物或含水合物地层的钻井作业提供借鉴。

关键词: 天然气水合物, 注热分解, 热-流-固耦合模型, 弹塑性本构模型, 正交数值模拟, 等效塑性应变最大值

Abstract: In response to the problem that the damage mechanism and influencing factors of the reservoir destabilization caused by the thermal decomposition of natural gas hydrate are not clear, the coupled heat-fluid-solid model of the hydrate deposition layer and the elastic-plasticity intrinsic model of the solid skeleton of the reservoir under the action of such coupled model have been established in consideration of the formation heat transfer, pore fluid seepage and deformation damage of solid skeleton due to the thermal decomposition of hydrate. Based on the secondary development of the USDFLD subroutine of ABAQUS software, the orthogonal numerical simulation experiments were carried out to analyze the sensitivity of the thermal injection temperature difference, the absolute permeability of the hydrate deposition layer, and the effective principal stress difference to the influence of the deformation damage of near-well reservoir in the the hydrate deposition layer and the influence mechanism. The results show that the sensitivity to the influence of the deformation damage of near-well reservoir is in descending order of the effective principal stress difference, absolute permeability, and thermal injection temperature difference; the deterioration of the mechanical properties of the near-well reservoir and the decrease of the effective stress caused by the thermal decomposition of hydrate are the main reasons for the deformation damage; the larger the effective principal stress difference is, the larger the maximum value of the equivalent plastic strain is, and the location of the most severe deformation damage in the near-well reservoir is always located in the direction of the minimum horizontal geostress in the borehole under the influence of the the non-uniformity of the horizontal geostress. This study provides a reference for the natural gas hydrate extraction with the thermal stimulation method and drilling operations of hydrate-bearing formations.

Key words: natural gas hydrate, thermal injection decomposition, thermo-hydro-mechanical coupled model, elasto-plastic constitutive model, orthogonal numerical simulation, maximum equivalent plastic strain

中图分类号: