青海共和盆地干热岩热储层人工水力致裂裂缝扩展规律
2023-07-20郭茂生姬长发刘宗鑫张鹏举
郭茂生 姬长发 刘宗鑫 张鹏举
摘 要：干熱岩热储层裂隙的渗透性、连通性直接影响着其换热效果，水力压裂是增强裂隙特性的重要手段。为了掌握水力压裂时裂缝的扩展规律，以青海共和盆地GR1地热井3 500～3 705 m深度段为研究对象，通过对ABAQUS水力压裂模拟软件进行二次开发，模拟研究了基质的弹性模量、抗拉强度、水平主应力差及压裂液排量等主要因素对水力裂缝长度和宽度的影响规律。结果表明：4个因素对水力裂缝长度和宽度的影响范围由大到小依次为压裂液排量、基质的弹性模量、基质的抗拉强度和水平主应力差；在压裂液排量一定时，水力裂缝长度随基质弹性模量、水平主应力差的增大而增大，随基质抗拉强度的增大而减小，宽度随基质弹性模量、水平主应力差的增大而减小，随基质抗拉强度的增大而增大；在弹性模量、抗拉强度和水平地应力差一定时，水力裂缝的长度和宽度均随着压裂液排量的增大而增大。研究结果可为干热岩热能的开发利用提供一定依据。
关键词：干热岩；水力压裂；水力裂缝；扩展规律；数值模拟
中图分类号：TK 521文献标志码：A
文章编号：1672-9315（2023）03-0514-09
DOI：10.13800/j.cnki.xakjdxxb.2023.0308开放科学（资源服务）标识码（OSID）：
Artificial hydraulic fracture propagation law of hot dry rock reservoir in Gonghe Basin，Qinghai province
GUO Maosheng，JI Changfa，LIU Zongxin，ZHANG Pengju
（College of Energy Science and Engineering，Xian University of Science and Technology，Xian 710054，China）
Abstract：The permeability and connectivity of cracks in hot dry rock reservoir directly affect the heat transfer effect.Hydraulic fracturing is an important means to enhance the fracture characteristics.In order to master the fracture propagation law during hydraulic fracturing，this paper takes GR1 geothermal well in Gonghe Basin，Qinghai province at a depth of 3 500～3 705 meters as the research object.Through the secondary development of ABAQUS hydraulic fracturing simulation software，the effects of matrix elastic modulus，tensile strength，horizontal principal stress difference and fracturing fluid displacement on the length and width of hydraulic fractures were simulated and examined.The results show that the influence range of the four factors on the hydraulic fracture length and width in descending order is the fracturing fluid displacement，elastic modulus of the matrix，tensile strength of the matrix and horizontal principal stress difference.The length of hydraulic fracture increases with the increase of matrix elastic modulus and horizontal principal stress difference，and decreases with the increase of matrix tensile strength; the width decreases with the increase of matrix elastic modulus and horizontal principal stress difference，and increases with the increase of matrix tensile strength.And when the elastic modulus，tensile strength and horizontal ground stress difference are certain，the length and width of hydraulic fracture increase with the increase of fracturing fluid discharge.The research results can provide some basis for the exploitation and utilization of thermal energy of dry hot rock.
Key words：hot dry rock；hydraulic fracturing；hydraulic fracture；propagation rules；numerical simulation
0 引 言
开发干热岩地热资源是“国家地热能开发利用‘十三五规划”的重点任务，也是建立中国能源战略体系的重要一环[1]。通过水力压裂等人工手段在干热岩中建立增强型地热系统（EGS）是提取热量的关键技术，也是难点技术之一。主要原因是压裂液排量、注水压力、水平主应力差、基质的弹性模量及抗拉强度等因素都对水力裂缝的扩展具有一定的影响[2-4]，而裂缝扩展形态的变化表现为裂缝长度、宽度及扩展方向的不同。
目前，国内外学者关于干热岩水力压裂裂缝形态的研究主要集中在数值模拟和物理试验两方面。王素玲等学者基于断裂力学和损伤力学理论，描述了裂缝表面上岩体的力学行为，分析岩石力学性质、压裂液粘度及排量对水力裂缝形态的影响[5-6]。李连崇等对水力裂缝的扩展过程进行真三维模拟研究，结果表明岩体介质的非均匀性是造成裂缝扩展的主要因素[7]。物理试验方面主要是通过真三轴水力压裂模拟试验的方式进行研究，国内外学者研究了不同储层中水力裂缝的形态变化[8-9]，以及水力裂缝遇到随机天然裂缝等弱面[10-12]时的扩展规律。程远方[13]利用真三轴水力压裂模拟试验系统，研究了主应力差对裂缝形态的影响[14]。郝家兴也通过试验研究，证明了水力裂缝扩展方向是层理方向及应力差公共作用的结果[15]。然而该方式选用的样本尺寸较小（300 mm×300 mm×300 mm），雖然通过试验可以较为明确地得到试样的起裂压力，但是由于试样较小，并不能完整的获取大尺寸干热岩破裂时裂缝长度及宽度等数据。并且，目前的研究主要是通过CT层析扫描、剪切波衍射技术、显微观测、添加染色剂[16]等水力压裂监测手段来获取裂缝的长度、面积、宽度等参数信息，因设备的局限性和试验的不确定性，在研究过程中，以这种方式获取裂缝面参数存在一定的困难和缺陷[17]。
限于试验研究的局限性以及水力压裂裂缝破裂面信息的研究缺乏定量分析[18]，以青海共和盆地GR1地热井3 500～3 705 m深度段干热岩热储层为地质背景，通过对ABAQUS水力压裂模拟软件进行二次开发，用以提取水力裂缝的长度和宽度，分析基质的弹性模量、抗拉强度、水平主应力差及压裂排量对裂缝长度及宽度的影响，为青海共和盆地干热岩的开发提供一定的参数依据。
1 物理模型
1.1 GR1地热井概况及岩体力学参数
高品位干热岩资源主要分布在青海共和盆地、海南北部及云南腾冲等地，其中青海共和盆地的地热资源占全国总地热资源的20.5%[19]。2017年5月中国科学家在青海共和盆地GR1地热井3 705 m处钻获的236 ℃高温岩体是国内首次发现的埋藏最浅且温度最高的干热岩体[20]。GR1地热井作为中国在青海共和盆地恰卜恰场地钻获的五口达到干热岩标准的地热井之一，就温度（236 ℃）和深度（3 705 m）而言，是国内目前最成功的地热钻井。目前，恰卜恰地热田已经圈定可开采干热岩面积达246.9 km2，资源总量13.66 EJ，折合标准煤约4.66亿t[21]。
通过对GR1地热井的钻井岩芯分析可得，在3 500～3 705 m深度段为印支期花岗岩，岩性以黑云母花岗岩、花岗闪长岩等深层侵入岩为主[22]。储层中存在较多的天然裂缝和断层，裂缝间距为0.3～50 m，裂缝宽度为0.1～1 mm[20]，如图1所示。储层的岩石力学及压裂过程中的参数[21]设置见表1。
1.2 建立物理模型
在ABAQUS中建立50 m×50 m的二维压裂模型，模型中设置有随机生成的天然裂缝。在模型中心位置设置有注入点和2条0.5 m长的预置裂缝，预置裂缝沿y轴方向分布，如图2所示。为了模拟真实的压裂过程，重现裂纹的扩展路径，使用PYTHON编程对ABAQUS模拟软件进行二次开发，在模型中全局嵌入0厚度的Cohesive单元，并将模拟时间设置为100 s。
在实际的工程压裂过程中，影响裂缝扩展状态的因素多种多样，如钻井的扰动、压裂液的注入、部分岩体的非均匀性和各向异性等因素都会对裂缝的扩展形态造成较大影响[23]，为了方便计算，做出如下假设。
1）假设储层岩体为各向同性的均匀材料。
2）假设压裂介质为不可压缩流体，且不与岩石发生化学反应。
3）忽略压裂过程对储层地应力的影响。
4）忽略高温储层对流体和岩石物性参数及物理性质的影响。
1.3 连续型方程
由质量守恒定律可知，一定时间内流入岩石内部流体质量等于流进与流出该单元的流体质量差。该过程可由达西渗流定律表述，推导可得流体渗流连续型方程见下式[24]。
1.4 边界条件
1）流量边界条件
式中 n为流量边界的单位法线方向；k为渗透率系数张量。
2）孔压边界条件。孔隙压力Pp=P0，即孔隙压力为定值P0。模拟均采用有效应力原理，即Pp=P0=0。
水平主应力差见式（3）。
Δσ=σH-σh（3）
式中 σH为最大水平主应力，MPa;σh为最小水平主应力，MPa。
3）位移边界条件。约束X和Y方向的节点位移为0，即Ux=0；Uy=0。
1.5 裂缝起裂准则
目前在国内外学者的研究中，判断裂缝起裂的常用准则一般有最大主应力准则、最大主应变准则、最大正应力准则、最大正应变准则、二次牵引准则等[25]，由于最大主应力准则具有计算容易收敛，稳定性较高，兼容性较强的优点，因此文中的模拟均使用最大主应力准则，只有当最大应力比达到某一临界值时，模型才会起裂。最大主应力准则见式（4）。
f={〈σmax〉/σamax}（4）
式中 σamax为临界最大主应力；〈〉表示在纯压缩应力作用下模型不会产生初始损失。
1.6 模拟验证
为验证模拟方法的准确性，以其他学者[17]的试验参数为依据，模拟相同情况下水力裂缝扩展情况，模拟时间设置120 s。岩石力学参数见表2。由于试验结果中裂缝宽度为0.028 mm，为提高模拟的精确度，通过修改ABAQUS软件中的inp文件，将裂缝的初始宽度设置为2×10-6。
图3为模拟过程中水力裂缝的最大宽度随时间的变化曲线。由图3可以看出，水力裂缝在114 s左右达到最大宽度，最大宽度为0.029 5 mm，之后水力裂缝穿透模型，裂缝宽度呈逐渐下降趋势。
试验中获得的裂缝宽度为0.028 mm。模拟结果与试验误差为5.35%，模拟结果具有较高的准确性。
2 模拟工况
为了具体分析基质弹性模量、抗拉强度、水平主应力差和压裂液排量这4个因素对裂缝长度和宽度的影响情况，通过控制变量法设置对照模拟，模拟工况及参数见表3。
3 模拟结果
压裂时间设置为100 s，压裂完成后，将模型中的破裂单元放大120倍获得裂缝的扩展形态，通过PYTHON编程对ABAQUS软件的后处理进行二次开发，以便于实时提取水力裂缝的长度和宽度数据，并绘制不同变量下水力裂缝长度和宽度随时间的变化曲线。水力裂缝形态如图4所示。
当基质弹性模量不同时，压裂完成后裂缝形态如图5所示，水力裂缝的长度和宽度随时间的变化分别如图6，图7所示。从图6，图7可以看出，水力裂缝长度随基质弹性模量的增大而增大，宽度随基质弹性模量的增大而减小。当基质弹性模量为28 GPa时，水力裂缝长度为25.31 m，宽度为6.20 mm；当基质弹性模量为44 GPa时，水力裂缝长度为32.57 m，宽度为5.57 mm；当基质弹性模量为60 GPa时，水力裂缝长度为37.79 m，宽度为4.98 mm。
基质抗拉强度不同时，压裂完成后裂缝形态如图8所示，水力裂缝的长度和宽度随时间的变化分别如图9，图10所示。
从图9，图10可以看出，水力裂缝长度随基质抗拉强度的增大而减小，宽度随基质抗拉强度的增大而增大。当基质抗拉强度为11.5 MPa时，水力裂缝长度为35.90 m，宽度为5.24 mm；当基质抗拉强度为15.7 MPa时，水力裂缝长度为32.57 m，宽度为5.57 mm；当基质抗拉强度为19.9 MPa时，水力裂缝长度为27.10 m，宽度为5.90 mm。
压裂液排量不同时，压裂完成后裂缝形态如图11所示，水力裂缝的长度和宽度随时间的变化分别如图12，图13所示。
从图12、图13可以看出，水力裂缝长度和宽度均随压裂液排量的增大而增大，当压裂液排量为3×10-4m3/s时，水力裂缝长度为19.11 m，宽度为4.74 mm；当压裂液排量为6×10-4m3/s时，水力裂缝长度为32.57 m，宽度为5.57 mm；当压裂液排量为9×10-4m3/s，水力裂缝长度为44.50 m，宽度为6.15 mm。
不同的水平主应力差下，压裂完成后裂缝形态如图14所示，水力裂缝的长度和宽度随时间的变化曲线分别如图15，图16所示。
从图15，图16可以看出，水力裂缝长度随水平主应力差的增大而增大，宽度随水平主应力差的增大而减小。当水平主应力差为0 MPa时，裂缝的长度为30.90 m，宽度为5.74 mm；当水平主应力差为5 MPa时，裂缝的长度为32.57 m，宽度为5.57 mm；当水平主应力差为10 MPa时，裂缝的长度为34.60 m，宽度为5.18 mm。
为对比各个变量对GR1地热井3 500 m处水力裂缝长度和宽度的影响范围，以基本方案为基准，分别计算同一变量的最大值和最小值对应的裂缝长度和宽度的变化量，如图17，图18所示。
从图17，图18可得，当基质弹性模量从28 GPa增大到60 GPa时，水力裂缝长度变化量从-22.29%升至16.03%，宽度变化量从11.31%降至-10.59%；当基质抗拉强度从11.59 MPa增大到19.9 MPa时，水力裂缝长度变化量从10.22%降至16.77%，宽度变化量从-5.56%升至5.92%；当压裂液排量从3×10-4m3/s增大到9×10-4m3/s时，水力裂缝長度变化量从-41.33%升至36.63%，宽度变化量从-14.90%升至10.41%；当水平主应力差从0 MPa升至10 MPa时，水力裂缝长度变化量从-5.12%降至6.25%，宽度变化量从3.05%升至-7.00%。
4个因素对水力裂缝长度和宽度的影响范围由大到小依次为压裂液排量、基质的弹性模量、基质的抗拉强度和水平主应力差。
4 结 论
1）4个因素中对水力裂缝长度和宽度的影响范围由大到小依次为压裂液排量、基质的弹性模量、基质的抗拉强度和水平主应力差。
2）压裂液排量一定时，水力裂缝长度随基质弹性模量、水平主应力差的增大而增大，随基质抗拉强度的增大而减小；宽度随基质弹性模量、水平主应力差的增大而减小，随基质抗拉强度的增大而增大。
3）弹性模量、抗拉强度和水平地应力差一定时，水力裂缝的长度和宽度均随着压裂液排量的增大而增大。
4）當基质弹性模量为60 GPa，压裂液排量为9×10-4m3/s，基质抗拉强度为11.5 MPa，水平主应力差为10 MPa时，水力裂缝的长度最大，为48.85 m。当基质弹性模量为28 GPa，压裂液排量为9×10-4m3/s，基质抗拉强度为19.9 MPa，水平主应力差为0 MPa时，水力裂缝的宽度最大为6.35 mm。
参考文献（References）：
[1] 周舟，金衍，曾义金，等.青海共和盆地干热岩地热储层水力压裂物理模拟和裂缝起裂与扩展形态研究[J].吉林大学学报（地球科学版），2019，49（5）：1425-1430.
ZHOU Zhou，JIN Yan，ZENG Yijin，et al.Experimental study on hydraulic fracturing physics simulation，crack initiation and propagation in hot dry rock geothermal reservoir in Gonghe Basin，Qinghai[J].Journal of Jilin University（Earth Science Edition），2019，49（5）：1425-1430.
[2]康向涛.煤层水力压裂裂缝扩展规律及瓦斯抽采钻孔优化研究[D].重庆：重庆大学，2014.
KANG Xiangtao.Research on hydraulic fracturing pro-pagation rules in coal seam and gas extraction drilling optimization[D].Chongqing：Chongqing University，2014.
[3]谢晓深，侯恩科，龙天文，等.浅埋缓倾斜煤层开采覆岩及地表裂缝发育规律与形成机理[J].西安科技大学学报，2022，42（2）：200-209.
XIE Xiaoshen，HOU Enke，LONG Tianwen，et al.Development law and formation mechanism of overburden and surfacecracks induced by inclined coal seam mining in shallow buried area[J].Journal of Xian University of Science and Technology，2022，42（2）：200-209.
[4]刘超，程成，崔娜.弹性模量对岩石注浆裂缝形成影响分析[J].西安科技大学学报，2017，37（5）：760-764.
LIU Chao，CHENG Cheng，CUI Na.Influence analysis of the elastic modulus on fractures formation law in rock grouting[J].Journal of Xian University of Science and Technology，2017，37（5）：760-764.
[5]王素玲，姜民政，刘合.基于损伤力学分析的水力压裂三维裂缝形态研究[J].岩土力学，2011，32（7）：2205-2210.
WANG Suling，JIANG Minzheng，LIU He.Study of hydraulic fracturing morphology based on damage mecha-nics analysis[J].Rock and Soil Mechanics，2011，32（7）：2205-2210.
[6]周治东，程万，魏子俊，等.基于BEM的水力裂缝起裂与扩展数值模拟[J].地球物理学进展，2020，35（2）：807-814.
ZHOU Zhidong，CHENG Wan，WEI Zijun，et al.Nume-rical simulation of hydraulic fracture initiation and propagation based on BEM[J].Progress in Geophysics，2020，35（2）：807-814.
[7]李连崇，梁正召，李根，等.水力压裂裂缝穿层及扭转扩展的三维模拟分析[J].岩石力学与工程学报，2010，29（S1）：3208-3215.
LI Lianchong，LIANG Zhengzhao，LI Gen，et al.Three-dimensional numerical analysis of traversing and twisted fractures in hydraulic fracturing[J].Chinese Journal of Rock Mechanics and Engineering，2010，29（S1）：3208-3215.
[8]金衍，陈勉，张旭东.天然裂缝地层斜井水力裂缝起裂压力模型研究[J].石油学报，2006，27（5）：124-126.
JIN Yan，CHEN Mian，ZHANG Xudong.Initiation pressure models for hydraulic fracturing of inclined wells in naturally fractured formation[J].Acta Petrolei Sinica，2006，（5）：124-126.
[9]金衍，张旭东，陈勉.天然裂缝地层中垂直井水力裂缝起裂压力模型研究[J].石油学报，2005，26（6）：113-114，118.
JIN Yan，ZHANG Xudong，CHEN Mian.Initiation pressure models for hydraulic fracturing of vertical wells in naturally fractured formation[J].Acta Petrolei Sinica，2005，26（6）：113-114，118.
[10]丁立钦.层理地层各向异性对井壁稳定性的影响规律研究[D].北京：中国地质大学，2020.
DING Liqin.Research on the influences of anisotropy on borehole stability in bedding formations[D].Beijing：China University of Geosciences，2020.
[11]REINICKE A，RYBACKI E，STANCHITS S.Hydraulic fracturing stimulation techniques and formation damage Mechanisms-implications from laboratory testing of tight sandstone-proppant systems[J].Geochemistry，2010，70（3）：107-117.
[12]ZHOU J，CHEN M，JIN Y.Analysis of fracture propagation behavior and fracture geometry using a tri-axial fracturing system in naturally fractured reservoirs[J].International Journal of Rock Mechanics Mining Science and Geomechanics，2008，45（7）：1143-1152.
[13]程远方，徐太双，吴百烈，等.煤岩水力压裂裂缝形态试验研究[J].天然气地球科学，2013，24（1）：134-137.
CHENG Yuanfang，XU Taishuang，WU Bailie，et al.Experimental study on the hydraulic fracturemorphology of coal bed[J].Natural Gas Geoscience，2013，24（1）：134-137.
[14]潘俊锋，马文涛，刘少虹，等.坚硬顶板水射流预制缝槽定向预裂防冲技术试验[J].岩石力学与工程学报，2021，40（8）：1591-1602.
PAN Junfeng，MA Wentao，LIU Shaohong，et al.A prevention technology of rock burst based on directional presplitting of water jet prefabricated slot in hard roof[J].Chinese Journal of Rock Mechanics and Engineering，2021，40（8）：1591-1602.
[15]郝家兴，刘垒.煤岩水力压裂后裂缝形态及增透效果试验研究[J].能源技术与管理，2019，44（6）：19-21.
HAO Jiaxing，LIU Lei.Experimental study on crack form and anti-reflection effect after hydraulic fracture of coal rock mass[J].Energy Technology and Management，2019，44（6）：19-21.
[16]竇锦爱，林业青，邵丰，等.页岩气储层孔隙结构表征技术及试验方法研究进展[J].西安科技大学学报，2020，40（6）：1019-1030.
DOU Jinai，LIN Yeqing，SHAO Feng，et al.Advances in characterization techniques and experimental methods of shale gas reservoir pore structure[J].Journal of Xian University of Science and Technology，2020，40（6）：1019-1030.
[17]马耕，张帆，刘晓，等.地应力对破裂压力和水力裂缝影响的试验研究[J].岩土力学，2016，37（S2）：216-222.
MA Geng，ZHANG Fan，LIU Xiao，et al.Experimental study of impact of crustal stress on fracture pressure and hydraulic fracture[J].Rock and Soil Mechanics，2016，37（S2）：216-222.
[18]冯彦军，康红普.定向水力压裂控制煤矿坚硬难垮顶板试验[J].岩石力学与工程学报，2012，31（6）：1148-1155.
FENG Yanjun，KANG Hongpu.Test on hand and stable roof control by means of directional hydraulic frcturing in coal mine[J].Chinese Journal of Rock Mechanics and Engineering，2012，31（6）：1148-1155.
[19]王世永.青海共和盆地干热岩岩体力学特征及渗透特性研究[D].北京：中国石油大学，2019.
WANG Shiyong.Study on mechanical properties and permeability characteristics of hot dry rock in Gonghe Basin，Qinghai Province[D].Beijing：China University of Petroleum，2019.
[20]雷治红.青海共和盆地干热岩储层特征及压裂试验模型研究[D].长春：吉林大学，2020.
LI Zhihong.Study on the characteristics of hot dry rock reservoir and fracturing test model in the Gonghe Basin，Qinghai Province[D].Changchun：Jilin University，2020.
[21]张森琦，严维德，黎敦朋，等.青海省共和县恰卜恰干热岩体地热地质特征[J].中国地质，2018，45（6）：1087-1102.
ZHANG Senqi，YAN Weide，LI Dunpeng，et al.Characteristics of geothermal geology of the Qiabuqia HDR in Gonghe Basin，Qinghai Province[J].Geology of China，2018，56（6）：1087-1102.
[22]张盛生，张磊，田成成，等.青海共和盆地干热岩赋存地质特征及开发潜力[J].地质力学学报，2019，25（4）：501-508.
ZHANG Shengsheng，ZHANG Lei，TIAN Chengcheng，et al.Occurrence gelogical charaeteristics and development potential of got dry rocks in Qinghai Gonghe Basin[J].Journal of Geomechanics，2019，25（4）：501-508.
[23]李玮，闫铁，毕雪亮.基于分形方法的水力压裂裂缝扩展机理[J].中国石油大学学报（自然科学版），2008，32（5）：87-91.
LI Wei，YAN Tie，BI Xueliang.Mechanism of hydraulically created fracture breakdown and propagation based on fractal method[J].Journal of China University of Petroleum（Natural Science Edition），2008，32（5）：87-91.
[24]龚迪光，曲占庆，李建雄，等.基于ABAQUS平台的水力裂缝扩展有限元模拟研究[J].岩土力学，2016，37（5）：1512-1520.
GONG Diguang，QU Zhanqing，LI Jianxiong，et al.Extended finite element simulation of hydraulic fracture based on ABAQUS platform[J].Rock and Soil Mechanics.2016，37（5）：1512-1520.
[25]王孝慧，姚卫星.复合材料胶接结构有限元分析方法研究进展[J].力学进展，2012，42（5）：562-571.
WANG Xiaohui，YAO Weixing.Research progress of finite element analysis methods for bonded structures of composite materials[J].Advances in Mechanics，2012，42（5）：562-571.
（責任编辑：刘洁）
