吉林大学学报(地球科学版) ›› 2015, Vol. 45 ›› Issue (2): 552-563.doi: 10.13278/j.cnki.jjuese.201502204

• 地质工程与环境工程 • 上一篇    下一篇

鄂尔多斯盆地深部咸水层二氧化碳地质储存热-水动力-力学(THM)耦合过程数值模拟

雷宏武, 李佳琦, 许天福, 王福刚   

  1. 吉林大学地下水资源与环境教育部重点实验室, 长春 130021
  • 收稿日期:2014-05-06 发布日期:2015-03-26
  • 作者简介:雷宏武(1985-),男,博士,主要从事多相流数值模拟程序的开发和应用研究,E-mail:hongwulei2008@aliyun.com
  • 基金资助:

    中国地质调查局工作项目(12120113006300);国土资源部公益性行业科研专项项目(201211063-06);吉林大学博士交叉学科研究项目(2012JC014)

Numerical Simulation of Coupled Thermal-Hydrodynamic-Mechanical (THM) Processes for CO2 Geological Sequestration in Deep Saline Aquifers at Ordos Basin, China

Lei Hongwu, Li Jiaqi, Xu Tianfu, Wang Fugang   

  1. Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China
  • Received:2014-05-06 Published:2015-03-26

摘要:

注入CO2到深部咸水层(CO2地质储存)被认为是一种直接有效地减少CO2向大气排放的途径。CO2地质储存涉及到热、水动力和力学耦合过程,该耦合过程是预测CO2在储层中的迁移转化、评价储层储存能力和分析潜在风险的关键。基于Terzaghi固结理论,在热-水动力(TH)耦合软件TOUGH2框架中加入了力学模块,形成了新的热-水动力-力学(THM)模拟器。结合鄂尔多斯盆地CO2捕获和储存(CCS)示范工程场地的地质、水文地质条件,采用新的THM模拟器数值分析了CO2注入后地层中的温度、压力、CO2饱和度、位移和有效应力的时空变化特征。结果显示:在井口保持8 MPa和35℃情况下,能够实现10万 t/a的CO2注入量;压力上升的范围远远大于CO2运移和温度降低的范围,注入20 a后,其最大距离分别达到接近边界10 km、620 m和100 m;位移和应力变化主要与压力变化相关,注入引起最大抬升为0.14 m,在注入井附近位置储层中有效应力变化水平方向要大于垂直方向,而在远井位置相反;注入引起井附近有效应力明显减小,从而导致了孔隙度和渗透率的增大,增强了CO2注入能力。

关键词: CO2地质储存, 热-水动力-力学耦合过程, 数值模拟, 鄂尔多斯盆地

Abstract:

Injecting CO2 into deep saline aquifers, referred to as CO2 geological sequestration (CGS), is considered to be a promising method to reduce the emission of anthropogenic CO2 to atmosphere. CGS involves coupled thermal-hydrodynamic-mechanical (THM) processes, which are important for predicting migration and transformation of injected CO2, evaluating the reservoir performance, and assessing the risk associated. Based on the Terzaghi consolidation theory, a coupled mechanical module is developed that is incorporated into the existing simulator TOUGH2, which is a well-established code for TH processes in subsurface flow systems. Based on the geological and hydrogeological conditions of the Ordos CCS Demonstration Project, the new THM simulator is used to numerically analyze the spatial and temporal distribution of the temperature, pressure, CO2 saturation, vertical displacement, and effective stress. The results show that 105 metric tons of CO2 per year can be finished under the injection condition of wellhead pressure 8 MPa and 35℃. The lateral distance with pressure buildup is more than that of CO2 plume and that of temperature decrease. They are 10 km, 620 m and 100 m, respectively after 20 years' injection. The vertical displacement and change of the effective stress are mainly related to the change of the pressure. The maximum surface uplift is about 0.14 m. The effective stress change is more significant in the horizontal direction than that in the vertical direction near the injection well, while it's in the opposite away from the injection well. Injection induces an obvious decrease in the effective stress but enhances the porosity and permeability, this, in turn increases the CO2 injectivity ultimately.

Key words: CO2 geological sequestration, coupled thermal-hydrodynamic-mechanical processes, numerical simulation, Ordos basin

中图分类号: 

  • P641.69

[1] IPCC (Intergovernmental Panel on Climate Change). Carbon Dioxide Capture and Storage[R]. New York: Cambridge University Press, 2005.

[2] IPCC (Intergovernmental Panel on Climate Change). Climate Change 2007: The Physical Science Basis: Contribution of Working Group I to the Fourth Assessment Report of the IPCC[R]. New York: Cambridge University Press, 2007.

[3] Bachu S, Adams J J.Sequestration of CO2 in Geological Media in Response to Climate Change: Capacity of Deep Saline Aquifers to Sequester CO2 in Solution[J]. Energy Conversion and Management, 2003, 44(20): 3151-3175.

[4] Rutqvist J, Wu Y S, Tsang C F, et al. A Modeling Approach for Analysis of Coupled Multiphase Fluid Flow, Heat Transfer, and Deformation in Fractured Porous Rock[J]. International Journal of Rock Mechanics & Mining Sciences, 2002, 39: 429-442.

[5] Rutqvist J, Tsang C F. A Study of Caprock Hydromechanical Changes Associated with CO2 Injection into a Brine Formation[J]. Environmental Geology, 2002, 42: 296-305.

[6] Rutqvist J, Birkholzer J, Cappa F, et al. Estimating Maximum Sustainable Injection Pressure During Geological Sequestration of CO2 Using Coupled Fluid Flow and Geomechanical Fault-Slip Analysis[J]. Energy Conversion and Management, 2007, 48:1798-1807.

[7] Rinaldi A P, Rutqvist J. Modeling of Deep Fracture Zone Opening and Transient Ground Surface Uplift at KB-502 CO2 Injection Well, In Salah, Algeria[J]. International Jounal of Greenhouse Gas Control, 2013, 12:155-167.

[8] Taron J, Elsworth D, Min K B. Numerical Simulation of Thermal-Hydrologic-Mechanical-Chemical Processes in Deformable, Fractured Media[J]. International Journal of Rock Mechanics & Mining Sciences, 2009, 46: 842-854.

[9] 于子望, 张延军, 张庆, 等. TOUGHREACT搭接FLAC3D算法[J]. 吉林大学学报:地球科学版, 2013, 43(1): 199-206. Yu Ziwang, Zhang Yanjun, Zhang Qing, et al. Algorithm of TOUGHREACT Links to FLAC3D[J]. Journal of Jilin University: Earth Science Edition, 2013, 43(1): 199-206.

[10] Winterfeld P H, Wu Y S. Parallel Simulation of CO2 Sequestration with Rock Deformation in Saline Aquifers[C].[S.l.]: SPE Reservoir Simulation Symposium, Society of Petroleum Engineers, 2011.

[11] Hu L T, Winterfeld P H, Fakcharoenphol P, et al. A Novel Fully-Coupled Flow and Geomechanics Model in Enhanced Geothermal Reservoirs[J]. Journal of Petroleum Science and Engineering, 2013, 107:1-11.

[12] Pruess K, Oldenburg C, George M. TOUGH2 User's Guide, Version 2.0[R]. Berkeley: Earth Science Division, Lawrence Berkeley National Laboratory, University of California, 1999.

[13] Pruess K, Garcia J, Kovscek T, et al. Code Intercomparison Builds Confidence in Numerical Simulation Models for Geological Disposal of CO2[J]. Energy, 2004, 29: 1431-1444.

[14] Pruess K, Garcia J. Multiphase Flow Dynamics During CO2 Disposal into Saline Aquifers[J]. Environmental Geology, 2002, 42: 282-295.

[15] Xu T, Sonnenthal E, Spycher N, et al. TOUGHREACT:A Simulation Program for Non-Isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geological Media: Application to Geothermal Injectivity and CO2 Geological Sequestration[J]. Computers & Geosciences, 2006, 32: 145-165.

[16] 许天福, 金光荣, 岳高凡, 等. 地下多组分反应溶质运移数值模拟:地质资源和环境研究的新方法[J]. 吉林大学学报:地球科学版, 2012, 42(5): 1410-1425. Xu Tianfu, Jin Guangrong, Yue Gaofan, et al. Subsurface Reactive Transport Modeling: A New Research Approach for Geo-Resources and Environments[J]. Journal of Jilin University: Earth Science Edition, 2012, 42(5): 1410-1425.

[17] Pruess K. ECO2N: A TOUGH2 Fluid Property Module for Mixtures of Water, NaCl, and CO2[R]. Berkeley: Earth Science Division, Lawrence Berkeley National Laboratory, University of California, 2004.

[18] 方云, 林彤, 谭松林. 土力学[M]. 武汉: 中国地质大学出版社, 2003. Fang Yun, Lin Tong, Tan Songlin. Soil Mechanics[M]. Wuhan: China University of Geosciences Press, 2003.

[19] Jaeger J C, Cook, N G W, Zimmerman R W.Fundamentals of Rock Mechanics[M].Malden: Blackwell Publishing, 2007.

[20] 任相坤, 崔永君, 步学朋, 等. 鄂尔多斯盆地CO2地质封存潜力分析[J]. 中国能源, 2010, 32(1): 29-32. Ren Xiangkun, Cui Yongjun, Bu Xuepeng, et al. The Potential Analysis of CO2 Geological Storage in Ordos Basin[J]. Energy of China, 2010, 32(1): 29-32.

[21] 吴秀章, 崔永君. 神华10万t/a CO2盐水层封存研究[J].石油学报: 石油加工, 2010(增刊1): 236-239. Wu Xiuzhang, Cui Yongjun. Study on 0.1 Mt/a CO2 Sequestration in Saline Formation[J]. Acta Petrolei Sinica: Petroleum Processing Section, 2010(Sup.1): 236-239.

[1] 李冠, 何文祥, 文志刚, 胡勇, 高小洋, 肖娜. 鄂尔多斯盆地陇东地区长71-2亚段不同源储组合的页岩油富集机理及富集模式[J]. 吉林大学学报(地球科学版), 2026, 56(3): 768-787.
[2] 夏英杰, 陈华斌, 陈健, 姚明宇, 杨海. 不同施工参数下深层页岩水力压裂数值模拟[J]. 吉林大学学报(地球科学版), 2026, 56(3): 924-936.
[3] 涂君汕, 陈辉, 邓居智, 陈康, 李焱, 余辉. 水库堤坝对高密度电法的影响及校正技术[J]. 吉林大学学报(地球科学版), 2026, 56(3): 1013-1025.
[4] 王健, 张媛瑗, 吴楠, 徐清海, 崔子岳, 刘显凤, 付清萌. 鄂尔多斯盆地志靖—安塞地区长7段夹层型页岩油储层特征及分类评价[J]. 吉林大学学报(地球科学版), 2026, 56(2): 453-468.
[5] 束龙仓, 位书静, 澈丽木格, 温中琦, 刘波. 西辽河平原生态输水的地下水响应及不确定性量化[J]. 吉林大学学报(地球科学版), 2026, 56(2): 647-660.
[6] 刘涛, 刘宗堡, 张可佳, 张岩, 张瑞雪, 刘晓文, 徐翠云.

基于深度学习的致密砂岩储层薄片图像生成与识别方法 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 724-738.

[7] 黄兴国, 翁央央, 韩丽. 基于广义递归卷积的孔隙黏弹地震波正演模拟[J]. 吉林大学学报(地球科学版), 2026, 56(1): 377-385.
[8] 刘静怡, 陈守民, 申战勇, 王筱晔, 董宪鹏, 张雯, 朱玉双, .

致密储层微观孔隙结构特征及其对渗流特征的影响——以鄂尔多斯盆地姬塬地区长6储层为例 [J]. 吉林大学学报(地球科学版), 2025, 55(5): 1434-1444.

[9] 周柳湘, 余思琴, 陈俊华, 刘城, 陈义, 张鑫鑫, . 潜孔锤与冲击钻机联合钻进工艺成孔碎岩过程数值模拟[J]. 吉林大学学报(地球科学版), 2025, 55(5): 1608-1618.
[10] 苑成旺, 张敏, 张少龙, 秦磊, 郭海洋, 墨海滢, 信旸. 长春地铁基坑支护结构变形与地表沉降特征[J]. 吉林大学学报(地球科学版), 2025, 55(3): 879-892.
[11] 于子望, 卢帅屹, 白林, 郑天琪. CO2地质封存岩石力学问题研究进展[J]. 吉林大学学报(地球科学版), 2025, 55(3): 930-942.
[12] 程瑶, 陆丹丹, 赵龙飞.

微波加热油页岩储层的热响应 [J]. 吉林大学学报(地球科学版), 2025, 55(2): 387-400.

[13] 刁国君, 何昕.

地铁明挖车站下穿越既有建筑物的局部暗扩挖施工关键技术 [J]. 吉林大学学报(地球科学版), 2025, 55(2): 536-549.

[14] 蔺学旻, 肖红琳. 缝洞型储层声波远探测测井响应模拟及应用——以塔河油田碳酸盐岩缝洞型储层为例[J]. 吉林大学学报(地球科学版), 2025, 55(1): 312-327.
[15] 魏虹羽, 李世超, 王伟安. 地球动力学数值模拟算法的应用现状与展望[J]. 吉林大学学报(地球科学版), 2025, 55(1): 98-124.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 吴远坤, 刘成林, 于春勇. 松辽盆地双城断陷深层原油成藏模式[J]. 吉林大学学报(地球科学版), 2024, 54(5): 1443 -1456 .
[2] 马荣,石建省,刘继朝. 人工内分泌网络模型在水文地质参数研究中的应用[J]. 吉林大学学报(地球科学版), 2013, 43(3): 914 -921 .
[3] 李宁, 王成文. 东北及邻区晚古生代地层接触关系与佳-蒙地块的形成和演化[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1331 -1340 .
[4] 胡大千,初凤友,姚 杰. 中太平洋YJA海山富钴结壳矿物组成与元素地球化学[J]. J4, 2006, 36(01): 32 -0037 .
[5] 姜 雪, 程日辉,于民凤. 裂谷地层的气候和构造控制:Zscape模型分析与在松辽盆地北安断陷的应用[J]. J4, 2006, 36(01): 54 -0059 .
[6] 李春柏,张新涛,刘 立,任延广,孟 鹏. 布达特群热流体活动及其对火山碎屑岩的改造作用--以海拉尔盆地贝尔凹陷为例[J]. J4, 2006, 36(02): 221 -0226 .
[7] 孟宪纲,薄万举,刘志广,刘勇,畅柳,李朝柱,王子平. 芦山7.0级地震与巴颜喀拉块体中东段的活动性[J]. 吉林大学学报(地球科学版), 2014, 44(5): 1705 -1711 .
[8] 陈欢庆, 梁淑贤, 舒治睿, 邓晓娟, 彭寿昌. 冲积扇砾岩储层构型特征及其对储层开发的控制作用——以准噶尔盆地西北缘某区克下组冲积扇储层为例[J]. 吉林大学学报(地球科学版), 2015, 45(1): 13 -24 .
[9] 贾大成,邢立新, 潘 军, M. J. van Bergen, H. van Roermund. 伊通上地幔剪切带捕虏体中富铝尖晶石的地球化学特征[J]. J4, 2006, 36(04): 497 -502 .
[10] 谢忠雷,杨佰玲,包国章,董德明. 茶园土壤不同形态镍的含量及其影响因素[J]. J4, 2006, 36(04): 599 -604 .