增强地热系统停止运行后温度恢复过程的数值模拟

doi:10.7523/j.ucas.2022.067

摘要/Abstract

摘要： 在增强地热系统(EGS)运行过程中，干热岩(HDR)人工储层的温度会迅速下降，外界的热量补充速度远远跟不上系统热量提取的速度，同时也会对人工储层的流体压力产生巨大的影响，EGS停止运行以后人工储层温度和压力场会逐渐的恢复到初始状态。前人对干热岩开采过程中的温度场和流体压力场变化做了大量研究，但很少关注停止开采以后人工储层的温度恢复过程和流体压力场变化，而停止开采以后的温度场恢复过程对于新的地热井选址和废弃地热井的重新利用具有重要意义。利用数值模拟软件COMSOL Multiphysics建立双井EGS运行过程和停止运行之后两个阶段的二维热-孔隙流体耦合数值模型，模型考虑了放射性生热的贡献和不同热导率的影响。研究表明：在EGS运行期间，热对流是热量交换的主要方式，人工储层温度下降过程很快，19 a后EGS已经到达运行寿命，68 a后HDR人工储层已经失去开采价值。而地热开采停止以后，人工储层温度的自然恢复过程长达万年，热传导是主要的热量传递和温度恢复方式，热导率越高，人工储层温度恢复越快。

关键词: 干热岩, 温度恢复, 热传导, 热导率

Abstract: During the operation of the enhanced geothermal system (EGS), the temperature of the artificial reservoir in the hot dry rock (HDR) will rapidly decline, and the external heat supply will be unable to keep up with the system’s rate of heat extraction. At the same time, it will also have a significant impact on the artificial reservoir’s fluid pressure. The predecessors conducted extensive research on the temperature field and fluid pressure field changes during HDR exploitation process, but little attention has been paid to the temperature recovery process and fluid pressure field changes of artificial reservoirs after production has ceased, this process is critical for the site selection of new geothermal wells and the reuse of abandoned geothermal wells. The numerical simulation software COMSOL Multiphysics is used in this paper to establish a two-dimensional thermal-hydraulic coupling numerical model of the two-well EGS during operation and after shutdown, considering the contribution of radioactive heat generation and the influence of different thermal conductivities. According to the findings, heat convection is the primary mode of heat exchange during EGS operation, and the temperature of the artificial reservoir declines rapidly. The EGS has reached the end of its operational life after 19 years, while the artificial reservoir has lost its exploitation value after 68 years. The natural recovery process of artificial reservoir temperature lasts 10 000 years after geothermal exploitation is terminated. The primary technique of heat transfer and temperature recovery is conduction. The higher the thermal conductivity, the sooner the artificial reservoir temperature recovers.

Key words: hot dry rock (HDR), temperature recovery, heat conduction, thermal conductivity

中图分类号:

P315

刘汉青, 胡才博, 赵桂萍. 增强地热系统停止运行后温度恢复过程的数值模拟[J]. 中国科学院大学学报, 2024, 41(2): 222-230.

LIU Hanqing, HU Caibo, ZHAO Guiping. Numerical simulation of temperature recovery process after enhanced geothermal system shutdown[J]. Journal of University of Chinese Academy of Sciences, 2024, 41(2): 222-230.

参考文献

[1] 汪集旸, 胡圣标, 庞忠和, 等. 中国大陆干热岩地热资源潜力评估[J]. 科技导报, 2012, 30(32): 25-31. DOI:10.3981/j.issn.1000-7857.2012.32.003.
[2] 许天福, 胡子旭, 李胜涛, 等. 增强型地热系统: 国际研究进展与我国研究现状[J]. 地质学报, 2018, 92(9): 1936-1947. DOI:10.3969/j.issn.0001-5717.2018.09.012.
[3] 甘浩男, 王贵玲, 蔺文静, 等. 增强型地热系统环境地质影响现状研究与对策建议[J]. 地质力学学报, 2020, 26(2): 211-220. DOI:10.12090/j.issn.1006-6616.2020.26.02.020.
[4] 关健, 方石. 地热系统的概念与传热机制综述[J]. 地质与资源, 2021, 30(2): 207-213, 206. DOI:10.13686/j.cnki.dzyzy.2021.02.013.
[5] 徐胜强, 张旭东, 张保平, 等. 测斜仪监测技术在共和盆地干热岩井压裂中的应用研究[J]. 探矿工程(岩土钻掘工程), 2021, 48(2): 42-48. DOI:10.12143/j.ztgc.2021.02.006.
[6] 周健, 曾义金, 陈作, 等. 青海共和盆地干热岩压裂裂缝测斜仪监测研究[J]. 石油钻探技术, 2021, 49(1): 88-92. DOI:10.11911/syztjs.2020139.
[7] Pandey S N. A three-dimensional coupled thermal-hydro model for enhanced geothermal systems[J]. International Journal of Renewable Energy Research, 2016, 6: 1516-1523. DOI:10.20508/ijrer.v6i4.4773.g6941.
[8] Asai P, Panja P, McLennan J, et al. Efficient workflow for simulation of multifractured enhanced geothermal systems (EGS)[J]. Renewable Energy, 2019, 131: 763-777. DOI:10.1016/j.renene.2018.07.074.
[9] Slatlem Vik H, Salimzadeh S, Nick H M. Heat recovery from multiple-fracture enhanced geothermal systems: The effect of thermoelastic fracture interactions[J]. Renewable Energy, 2018, 121: 606-622. DOI:10.1016/j.renene.2018.01.039.
[10] 肖鹏, 闫飞飞, 窦斌, 等. 增强型地热系统水平井平行多裂隙换热过程数值模拟[J]. 可再生能源, 2019, 37(7): 1091-1099. DOI:10.13941/j.cnki.21-1469/tk.2019.07.023.
[11] 孙致学, 姜传胤, 张凯, 等. 基于离散裂缝模型的CO₂增强型地热系统THM耦合数值模拟[J]. 中国石油大学学报(自然科学版), 2020, 44(6): 79-87. DOI:10.3969/j.issn.1673-5005.2020.06.010.
[12] Ma Y Y, Li S B, Zhang L G, et al. Analysis on the heat extraction performance of multi-well injection enhanced geothermal system based on leaf-like bifurcated fracture networks[J]. Energy, 2020, 213: 118990. DOI:10.1016/j.energy.2020.118990.
[13] Figueiredo B, Tsang C F, Niemi A. The influence of coupled thermomechanical processes on the pressure and temperature due to cold water injection into multiple fracture zones in deep rock formation[J]. Geofluids, 2020, 2020: 8947258. DOI:10.1155/2020/8947258.
[14] 段云星, 杨浩. 增强型地热系统采热性能影响因素分析[J]. 吉林大学学报(地球科学版), 2020, 50(4): 1161-1172. DOI:10.13278/j.cnki.jjuese.20190041.
[15] 陈炫沂, 姜振蛟, 徐含英, 等. 共和盆地干热岩体人工裂隙带结构的控热机理与产能优化[J]. 水文地质工程地质, 2022, 49(1): 191-199. DOI:10.16030/j.cnki.issn.1000-3665.202104058.
[16] Sun Z X, Zhang X, Xu Y, et al. Numerical simulation of the heat extraction in EGS with thermal-hydraulic-mechanical coupling method based on discrete fractures model[J]. Energy, 2017, 120: 20-33. DOI:10.1016/j.energy.2016.10.046.
[17] Lei Z H, Zhang Y J, Zhang S Q, et al. Electricity generation from a three-horizontal-well enhanced geothermal system in the Qiabuqia geothermal field, China: slickwater fracturing treatments for different reservoir scenarios[J]. Renewable Energy, 2020, 145: 65-83. DOI:10.1016/j.renene.2019.06.024.
[18] Liu H J, Wang H W, Lei H W, et al. Numerical modeling of thermal breakthrough induced by geothermal production in fractured granite[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2020, 12(4): 900-916. DOI:10.1016/j.jrmge.2020.01.002.
[19] 张盛生, 张磊, 田成成, 等. 青海共和盆地干热岩赋存地质特征及开发潜力[J]. 地质力学学报, 2019, 25(4): 501-508. DOI:10.12090/j.issn.1006-6616.2019.25.04.048.
[20] 张超, 张盛生, 李胜涛, 等. 共和盆地恰卜恰地热区现今地热特征[J]. 地球物理学报, 2018, 61(11): 4545-4557. DOI:10.6038/cjg2018L0747.
[21] 雷治红. 青海共和盆地干热岩储层特征及压裂试验模型研究[D]. 长春: 吉林大学, 2020.
[22] 王世永. 青海共和盆地干热岩岩体力学特征及渗透特性研究[D]. 北京: 中国石油大学(北京), 2019.
[23] 陈必光, 宋二祥, 程晓辉. 二维裂隙岩体渗流传热的离散裂隙网络模型数值计算方法[J]. 岩石力学与工程学报, 2014, 33(1): 43-51. DOI:10.13722/j.cnki.jrme.2014.01.005.
[24] 马哲民, 谭现锋, 郝俊杰, 等. 多种测温方法在青海共和干热岩GR1井中的应用[J]. 探矿工程(岩土钻掘工程), 2020, 47(12): 42-48. DOI:10.12143/j.tkgc.2020.12.007.
[25] 岳高凡, 邓晓飞, 邢林啸, 等. 共和盆地增强型地热系统开采过程数值模拟[J]. 科技导报, 2015, 33(19): 62-67. DOI:10.3981/j.issn.1000-7857.2015.19.010.
[26] 丁军锋, 王世民. 增强型地热系统的多区域多物理场耦合三维数值模拟[J]. 中国科学院大学学报, 2019, 36(5): 694-701. DOI:10.7523/j.issn.2095-6134.2019.05.015.
[27] 孙致学, 徐轶, 吕抒桓, 等. 增强型地热系统热流固耦合模型及数值模拟[J]. 中国石油大学学报(自然科学版), 2016, 40(6): 109-117. DOI:10.3969/j.issn.1673-5005.2016.06.014.
[28] 崔翰博, 唐巨鹏, 姜昕彤. 渗透率对干热岩开采过程储层变化规律的影响[J]. 水文地质工程地质, 2020, 47(1): 171-180. DOI:10.16030/j.cnki.issn.1000-3665.201811014.
[29] 张超, 胡圣标, 宋荣彩, 等. 共和盆地干热岩地热资源的成因机制:来自岩石放射性生热率的约束[J].地球物理学报, 2020, 63(7): 2697-2709. DOI:10.6038/cjg2020N0381.
[30] 张伟, 孙江, 曲占庆, 等. 高温地热开采热流固耦合模型及综合评价方法[J]. 地球物理学进展, 2019, 34(2): 668-675. DOI:10.6038/pg2019CC0445.
[31] John W L，Bloomquist R G，Tonya L B，et al.The united states of America country update[C]//Proceedings World Geothermal Congress 2005. Antalya, Turkey: Transactions-geothermal Resources Council, 2005: 1-19.