›› 2014, Vol. 35 ›› Issue (8): 2415-2425.

• 数值分析 • 上一篇    下一篇

地下热-水动力-力学耦合过程数值模拟:以CO2地质储存为例

雷宏武,金光荣,石 岩,李佳琦,王福刚,许天福   

  1. 吉林大学 地下水资源与环境教育部重点实验室,长春 130021
  • 收稿日期:2013-08-27 出版日期:2014-08-12 发布日期:2014-08-14
  • 作者简介:雷宏武,男,1985年生,博士研究生,主要从事多相流数值模拟程序的开发和应用方面的研究工作。
  • 基金资助:

    中国地质调查局工作项目(No. 12120113006300);吉林大学博士交叉学科研究项目(No. 2012JC014);研究生创新基金(No. 20121069)。

Numerical simulation of subsurface coupled thermo-hydro-mechanical (THM) processes: Application to CO2 geological sequestration

LEI Hong-wu, JIN Guang-rong, SHI Yan, LI Jia-qi, WANG Fu-gang, XU Tian-fu   

  1. Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China
  • Received:2013-08-27 Online:2014-08-12 Published:2014-08-14

摘要: 在地下流动系统问题的研究中,热-水动力-力学(THM)耦合过程是研究的热点问题。在地下多相非等温数值模拟软件TOUGH2的框架内,基于Biot固结理论和摩尔-库仑破坏判定准则,建立了THM耦合模型;采用积分有限差和有限元联合的空间离散方法,开发了THM模拟器TOUGH2Biot。该模拟器中热和水动力过程是全耦合,力学过程是部分耦合。通过与解析解的对比,验证了其正确性。基于鄂尔多斯盆地CCS示范工程,采用TOUGH2Biot研究了CO2注入地层后的THM响应。结果显示CO2的注入引起流体压力急剧增加,地层有效应力减小,地表隆起,隆起大小在几十个厘米,同时孔渗增加,利于CO2注入引起的压力上升向外消散。CO2注入最有可能导致剪切破坏的位置位于最大速率注入点上部盖层,其次为靠近地表的位置。

关键词: 热-水动力-力学耦合过程, Biot固结理论, CO2地质储存, 数值模拟

Abstract: Coupled thermo-hydro-mechanical (THM) process is an important and hot-spot issue. Based on Biot consolidation theory and Morh-Coulomb failure criteria, coupled THM model is built in the framework of TOUGH2 which is a multiphase fluid flow, nonisothermal simulation software. At the same time, a new simulator-TOUGH2Biot is developed with fully coupling between thermal and hydrodynamic processes and partly coupling for mechanics, using mixed integral finite difference and finite element method. Analytical solution is employed to verify the simulator. Finally, TOUGH2Biot is used to numerically analyze the THM response of formations after CO2 injection. The results show that the CO2 injection induces to a quick increase in fluid pressure and corresponding decrease in effective stress. Consequently, there is an uplift of tens of centimeters at the surface. Also, the porosity and permeability are enhanced, which are favorable for dissipation of pressure buildup. The most possibility of shear slip failure is at the caprock above the maximum injection rate location. The second potential failure location is near the surface.

Key words: coupled thermo-hydro-mechanical (THM) process, Biot consolidation theory, CO2 geological sequestration, numerical simulation

中图分类号: 

  • O 351
[1] 李翻翻, 陈卫忠, 雷江, 于洪丹, 马永尚, . 基于塑性损伤的黏土岩力学特性研究[J]. 岩土力学, 2020, 41(1): 132-140.
[2] 夏 坤, 董林, 蒲小武, 李璐, . 黄土塬地震动响应特征分析[J]. 岩土力学, 2020, 41(1): 295-304.
[3] 郭院成, 李明宇, 张艳伟, . 预应力锚杆复合土钉墙支护体系增量解析方法[J]. 岩土力学, 2019, 40(S1): 253-258.
[4] 闫国强, 殷跃平, 黄波林, 张枝华, 代贞伟, . 三峡库区巫山金鸡岭滑坡成因机制与变形特征[J]. 岩土力学, 2019, 40(S1): 329-340.
[5] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[6] 金爱兵, 刘佳伟, 赵怡晴, 王本鑫, 孙浩, 魏余栋, . 卸荷条件下花岗岩力学特性分析[J]. 岩土力学, 2019, 40(S1): 459-467.
[7] 韩征, 粟滨, 李艳鸽, 王伟, 王卫东, 黄健陵, 陈光齐, . 基于HBP本构模型的泥石流动力过程SPH数值模拟[J]. 岩土力学, 2019, 40(S1): 477-485.
[8] 吴锦亮, 何吉, . 岩质边坡动态开挖模拟的复合单元模型[J]. 岩土力学, 2019, 40(S1): 535-540.
[9] 吴凤元, 樊赟赟, 陈剑平, 李军, . 基于不同侵蚀模型的高速崩滑碎屑 流动力过程模拟分析[J]. 岩土力学, 2019, 40(8): 3236-3246.
[10] 孙峰, 薛世峰, 逄铭玉, 唐梅荣, 张翔, 李川, . 基于连续损伤的水平井射孔-近井筒三维破裂模拟[J]. 岩土力学, 2019, 40(8): 3255-3261.
[11] 穆锐, 浦少云, 黄质宏, 李永辉, 郑培鑫, 刘 旸, 刘 泽, 郑红超, . 土岩组合岩体中抗拔桩极限承载力的确定[J]. 岩土力学, 2019, 40(7): 2825-2837.
[12] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[13] 张 聪, 梁经纬, 阳军生, 曹 磊, 谢亦朋, 张贵金, . 堤坝脉动注浆浆液扩散机制及应用研究[J]. 岩土力学, 2019, 40(4): 1507-1514.
[14] 严 健, 何 川, 汪 波, 蒙 伟, . 高地温对隧道岩爆发生的影响性研究[J]. 岩土力学, 2019, 40(4): 1543-1550.
[15] 李世俊, 马昌慧, 刘应明, 韩玉珍, 张 彬, 张 嘎, . 离心模型试验与数值模拟相结合研究 采空边坡渐进破坏特性[J]. 岩土力学, 2019, 40(4): 1577-1583.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 谭贤君,陈卫忠,杨建平,杨春和. 盐岩储气库温度-渗流-应力-损伤耦合模型研究[J]. , 2009, 30(12): 3633 -3641 .
[2] 魏 星,王 刚,余志灵. 交通荷载下软土地基长期沉降的有限元法[J]. , 2010, 31(6): 2011 -2015 .
[3] 温世亿,李静,苏霞,姚雄. 复杂应力条件下围岩破坏的细观特征研究[J]. , 2010, 31(8): 2399 -2406 .
[4] 毛 宁,张尧亮. 经验公式简便求法典型实例[J]. , 2010, 31(9): 2978 -2982 .
[5] 刘 杰,李建林,屈建军,陈 星,李剑武,骆世威. 基于卸荷岩体力学的大岗山坝肩边坡水平位移发育的多因素影响分析[J]. , 2010, 31(11): 3619 -3626 .
[6] 李伟华,赵成刚,杜楠馨. 软弱饱和土夹层对地铁车站地震响应的影响分析[J]. , 2010, 31(12): 3958 -3963 .
[7] 韩现民. 西格二线关角隧道浅埋砂层段施工技术及力学效应研究[J]. , 2010, 31(S2): 297 -302 .
[8] 蒋臻蔚,彭建兵,王启耀. 西安市地铁3号线不良地质问题及对策研究[J]. , 2010, 31(S2): 317 -321 .
[9] 刘用海,朱向荣,常林越. 基于Casagrande法数学分析确定先期固结压力[J]. , 2009, 30(1): 211 -214 .
[10] 李兴高,刘维宁. 挡土结构上水-土压力分算的进一步探讨[J]. , 2009, 30(2): 419 -424 .