干热花岗岩压裂裂缝剪切滑移数值模拟研究.pdfVIP

干热花岗岩压裂裂缝剪切滑移数值模拟研究 1，2 1，2* 徐志鹏 ,盛茂 1.石油工程学院，中国石油大学（北京），北京 2.油气资源与探测国家重点实验室，中国石油大学（北京），北京 * 通讯作者，shengmao@cup.edu.cn 前言 本文 以干热花岗岩为研究对象，基于 COMSOL 多物理场计算平台，建立考虑热-流- 固 干热岩（简称 HDR ）是一种不含水或蒸汽的 耦合作用的干热岩天然裂缝剪切滑移模型，得到 热岩体，温度 150℃以上，其中蕴含着巨大的地 了裂缝剪切滑移过程中加速度、速度、位移等因 热资源，可用于发电和热水利用[1] 。据中国地质 素随滑移距离的变化规律，揭示了法向加载路径、 调查局数据[2] ，我国埋深在3-10 km 的干热岩资源 切向载荷、缝面摩擦系数等因素对极限滑移距离 量折合 856 万亿吨标准煤，是2016 年全国能源消 的影响规律。研究结果有望深入认识干热岩压裂 费总量（43.6 亿吨标准煤）的20 万倍，资源量巨 增产机理。 大。干热岩岩性致密且天然裂缝发育，但裂缝连 通性差，自然状态下很难获取具有经济价值的热 1. 考虑温度应力的裂缝剪切滑移模型 能[3] 。要达到干热岩地热能的经济高效开采，压 裂改造是目前首选的增产措施。 1.1 物理模型 干热岩压裂增产机理不同于常规油气压裂， 干热岩压裂通常保持注入压力低于最小水平 能否充分激活天然裂缝、剪切滑移成网、提高储 主应力，裂缝在缝内水压力作用下，缝面有效正 层连通性是获得良好改造效果的关键。干热岩压 压力先减后增，裂缝滑移速度先增后减。将花岗 裂时监测到的微地震事件，被解释为天然裂缝的 岩剪切滑移简化为平面应力问题，且不考虑体积 剪切滑移，验证了天然裂缝对干热岩裂缝系统起 决定性作用[4] 。Jung (2013) [5]提出翼状裂缝概念， 力，岩石为线弹性材料。 即天然裂缝错动一定程度后会在裂缝的两端沿一 定倾角扩展，然后沟通其他天然裂缝。Mcclure and Horne (2012) [6]对张开型和剪切型两种裂缝扩 展模式的缝间流体压力临界值进行了研究。 Mcclure and Horne (2013) [7] 总结前人成果提出四 种裂缝扩展模式：纯剪切型，纯张开型，混合型 和主裂缝尖端滤失型。Ghassemi 等 (2016) [8]]用理 图1 花岗岩剪切滑移过程示意图 论分析和数值模拟，研究储层的热-孔隙弹性效应 Fig. 1 Schematic diagram of hydro-shearing process 对渗透率和地震活动（岩石破裂）的影响。结果 表明：注入流体的水压作用和热应力作用使裂缝 1.2 材料属性 面正应力产生显著变化。随后，Bauer 与 [9] 花岗岩杨氏模量 10~85GPa （模拟不同硬度岩 Ghassemi 等(2016) 搭建了干热岩注采换热室内 3 石），泊松比 0.3，密度 2800kg/m 。上盘岩石的 模拟实验装置，结果表明，裂缝剪切滑移量受到 体积为 50