超高温条件下花岗岩力学性质演化规律模拟研究_赫文豪.pdf

华南师范大学学报(自然科学版)Journal of South China Normal University(Natural Science Edition)2022，54(6):1622doi:106054/jjscnun2022080收稿日期:20220307华南师范大学学报(自然科学版)网址:http:journalnscnueducn基金项目:国家重点研发计划项目(2019YFA0708302);国家自然科学基金项目(U19B600305，52004296，52274016);地热资源开发技术与装备教育部工程研究中心开放课题(21018);中国石油大学(北京)科研基金项目(2462022YXZZ007，2462022BJC012)*通信作者:史怀忠，Email:shz cupeducn超高温条件下花岗岩力学性质演化规律模拟研究赫文豪1，2，臧雨3，刘利1，2，史怀忠3*，黄中伟3，宋先知3，王天宇3，4，刘子龙1，2(1 中国石油大学(北京)油气光学探测技术北京市重点实验室，北京 102249;2 中国石油大学(北京)理学院能源交叉学科基础研究中心，北京 102249;3 中国石油大学(北京)油气资源与探测国家重点实验室，北京 102249;4 地热资源开发技术与装备教育部工程研究中心，长春 130026)摘要:为提高深层干热岩破碎效率，以干热岩地层常见岩性花岗岩为研究对象，提出超高温条件下花岗岩力学性质演化规律研究。鉴于室内实验方法的局限性，主要采用数值模拟和理论计算等研究方法，建立了高温花岗岩巴西圆盘劈裂数值模型，探究了硬质花岗岩在温度压力联合作用下岩石应力场分布特征与扰动机制。研究发现:在单轴径向压缩载荷下，高温花岗岩各向应力分布受到显著扰动，岩石沿横轴拉伸应力、两端载荷施加位置附近压缩应力和剪切应力均显著减小，导致裂纹越容易扩展。温度越高，岩石损伤越明显，其岩石抗拉强度越低。研究结果揭示了温度作用下花岗岩力学性质演化规律，可为深层干热岩资源的高效开发提供理论依据。关键词:干热岩;花岗岩;温度;巴西劈裂;拉伸;应力场中图分类号:P314文献标志码:A文章编号:10005463(2022)06001607Numerical Study on Evolution Law of the Mechanical Properties ofGranite under UltraHigh Temperature ConditionsHE Wenhao1，2，ZANG Yu3，LIU Li1，2，SHI Huaizhong3*，HUANG Zhongwei3，SONG Xianzhi3，WANG Tianyu3，4，LIU Zilong1，2(1 Beijing Key Laboratory of Oil and Gas Optical Detection Technology，China University of PetroleumBeijing，Beijing 102249，China;2 Basic esearch center for Energy Interdisciplinary，China University of PetroleumBeijing，Beijing 102249，China;3 State Key Laboratory of Petroleum esources and Prospecting，China University of PetroleumBeijing，Beijing 102249，China;4 Engineeringesearch Center of Geothermal esources Development Technology and Equipment(Ministry of Education)，Jilin University，Changchun 130026，China)Abstract:To improve the rockbreaking efficiency of deep dry hot rock，taking granite，a common lithology of dryhot rock formation，as the research object，the changes in the mechanical properties of granite under ultrahightemperature was proposed In view of the limitations of laboratory experimental methods，the Brazilian disc splittingnumerical model of hightemperature granite was established by using numerical simulation and theoretical calcula-tion methods，and the distribution characteristics and disturbance mechanism of stress field of hard granite undertemperaturepressure conditions were investigated It is found that under uniaxial and radial compression load，thestress distributions of hightemperature granite get significantly disturbed，leading to lower tensile stress along thehorizontal direction，smaller compression and shear effect near the loading points，and easier fracture conductionThe higher the temperature，the more obvious the rock damage and the lower the tensile strength of the rock Theresults reveal the evolution law of mechanical properties of granite under the effect of temperature，which can pro-vide theoretical basis for efficient development of deep dry hot rock resourcesKeywords:hot dry rock;granite;temperature;Brazilian splitting;tension;stress field当前，我国对能源需求越来越大，在大力发展传统化石能源的基础上，也需要积极拓展其他能源形式。在“十三五”能源规划中，我国计划今后要大力开发清洁可再生能源，以满足能源需求，并实现节能减排的战略目标，“十四五”期间，习近平总书记提出并不断深化低碳经济发展战略，使得新能源相关技术及产业获得了进一步发展。为保障我国能源安全，地热资源因其清洁、高效、储量巨大且可持续发展等特点愈发凸显其优势，有望成为我国能源续接的重要路径。根据国土资源部 2015 年统计数据，我国中深层(2004 000 m)中低温资源量合 13 700 亿吨标准煤当量、高温资源发电潜力为 8 466 MW，深层(310 km)干热岩资源总计为 2091025J，合860 万亿吨标准煤，约为 2014 年全国能源消耗总量的26 万倍，但截至 2015 年底我国仅实现年替代标煤290 万吨13。因此，加快干热岩等深部地热资源的高效开发和利用，对改善我国能源战略布局、优化能源结构及促进转型发展具有重大政治、经济意义。但是干热岩地层，尤其是最典型的火成岩花岗岩地层，所处地层环境复杂，具有温度高、硬度大、抗压强度高、研磨性强、钻头难以吃入地层等特征，普通钻井技术不能实现其高效钻进。为有效推进干热岩地热资源的开发利用，明确超高温作用下的岩石力学性质演化规律尤为重要。围绕花岗岩这一干热岩常见岩性，前人研究结果表明，当花岗岩所处温度升高时，除了外观(如颜色等)的显著差异，岩石体积膨胀，骨架密度逐渐下降，内部水分不断蒸发，在岩石加热过程中岩石微裂纹开始产生并逐渐沟通形成宏观裂纹，岩石矿物骨架遭到破坏，且随着岩石损伤的不断增强，其宏观刚度和强度等力学指标显著降低，蠕变特征和塑性变形趋于明显，岩石的脆性降低，且存在阈值温度使得岩石力学、物理及化学指标(如单轴抗压强度、弹性模量等)发生明显改变49。同时，岩石颗粒尺寸差异越大，其阈值温度也就越低，岩石强度变化愈明显10。另外，加热源类型、升温速率、外界温差、循环次数等作业参数对花岗岩损伤效果、阈值温度和峰值强度具有显著影响9，1114。但目前开展的高温作用下干热岩力学性质研究多集中于花岗岩室内特性测试，对高温作用下岩石受载荷的应力场分布特征和裂纹延展机制尚缺乏理论研究，且干热岩室内力学实验对加热极限温度和实验设备要求较高，导致实验研究难度较大，干热岩破碎机理亟需进一步研究。为揭示干热岩破碎机理，以干热岩地层常见的花岗岩为研究对象，本研究提出超高温条件下花岗岩力学性质演化规律研究，建立高温花岗岩巴西圆盘劈裂数值模型，通过预置硬质花岗岩的岩石力学参数和温度参数，探究硬质花岗岩在温度载荷联合作用下岩石应力场分布特征与扰动机制，从而揭示温度对岩石单轴抗拉强度影响规律和岩石在超高温作用下裂缝延展机制，为温度作用下花岗岩力学性质演化规律和深层干热岩资源的高效开发提供理论依据。1有限元破岩模型构建为揭示巴西劈裂过程中干热岩破碎机制，基于岩石力学、断裂力学等学科基础理论，参照干热岩与温度相关力学参数，利用 Abaqus 有限元模拟平台，建立温度载荷联合作用下的硬质花岗岩巴西劈裂圆盘数值模型(图 1)。图 1硬质花岗岩巴西劈裂圆盘数值模型Figure 1The numerical modelling of Brazilian test of granites鉴于花岗岩脆性破坏模式，定义岩石破碎本构关系为 DruckerPrager 准则且设定其硬化特征，并以黏聚力单元定义岩石破坏模式，巴西劈裂过程中，将上下承压板视为刚体部件，为模拟实际岩石破碎过程，设定载荷过程中上下承压板均以 1 mm/s 压缩岩石，持续压缩时间为 1 s(包含 40 个等时间隔，间隔 0025 s)。由于加热过程中岩石密度与体积变化相对较小，模拟过程中统一岩石密度和泊松比分别为27 g/mm3和0200，且以室温(25)条件下岩石强度作为参照对象，模拟岩石所处温度范围为150650。为进行超高温条件下花岗岩力学损伤性质的模71第 6 期赫文豪等:超高温条件下花岗岩力学性质演化规律模拟研究拟，表 1 中数据包括文章调研数据和必要插值，具体数据详见文献 15 18。岩石与承压板三棱形切削齿的接触形式采用通用显式接触，切向摩擦系数设置为 03，并定义岩石与承压板接触面的接触关系为“硬接触”。通过观测岩石破碎过程中承压板受到的反作用力可以获得岩石在巴西劈裂过程中的受压载荷，而应力场动态分布特征则可帮助揭示干热岩裂纹延展机制。表 1巴西劈裂圆盘材料主要物性参数Table 1The main material parameters of Brazilian test modelling材料温度/密度/(g mm3)弹性模量/GPa泊松比单轴抗压强度/MPa承压板78210000077花岗岩2527393002009042150273226020087732002729220200863530027274302007936400271897020083775002720990200682060027198702006956650271616020059832有限元模拟结果分析以 25 花岗岩模拟岩样作为参照，本模拟研究重点阐述超高温作用下花岗岩巴西劈裂圆盘的岩石应力场分布特征、裂纹延展特征和巴西劈裂拉伸强度特征。21超高温作用下岩石应力场分布特征为明确干热岩应力场分布特征，本节模拟未考虑裂缝等非均质性影响因素，通过定义岩石在不同温度下的强度参数，可以获得花岗岩巴西劈裂圆盘在横向和纵向上的应力分布(图 2、图 3)，其中暖色(如红色)表示拉伸效应，冷色(如蓝色)表示压缩效应。通过观察不同加热处理样品的应力图谱分布，可以得到岩石的应力场分布特征。为对比不同温度处理下的花岗岩巴西劈裂圆盘应力场分布特征，需要对各模拟结果作统一图例处理，以便更加直观、高效地观察花岗岩圆盘各处实时应力大小与受力状态。图 2花岗岩巴西劈裂圆盘横向应力分布Figure 2The horizontal stress distribution of granite Brazilian disk从横向(x 轴方向)应力分布图(图 2)来看，若统一横向应力图谱显示图例为100100 MPa，岩石最大横向拉伸应力发生在圆盘中心处，沿 x 轴方向向外逐渐递减，但在靠近载荷的 y 轴纵向方向呈81华 南 师 范 大 学 学 报(自 然 科 学 版)第 54 卷不断减小趋势并逐渐减小至 0，在载荷施加处及其附近区域横向拉伸应力转变为压缩应力，载荷位置附近压缩效应显著且呈对称状扇形分布，且随着压缩时间的增大而越发明显。随着岩石处理温度从25 升高至 650，花岗岩巴西劈裂圆盘各处横向应力显著降低，圆盘中心横向拉伸效应减弱，两端载荷位置附近区域压缩效应亦显著减弱，表明超高温作用可以削弱岩石在横向的应力分布。从纵向(y 轴方向)应力分布图(图 3)来看，若统一纵向应力图谱显示图例为500500 MPa，岩石最大纵向压缩应力发生在圆盘载荷施加处，沿 y 轴方向向圆心逐渐递减，越靠近施加载荷初，载荷纵向压缩应力越大，相较于 x 轴方向压缩应力分布，沿载荷方向纵向压缩效应显著，且随着压缩时间的增大而越发明显。模拟发现，圆盘纵向拉伸效应微弱，拉伸应力可以忽略。随着岩石处理温度从 25 升高至 650，花岗岩巴西劈裂圆盘各处压缩应力显著降低，圆盘两端载荷位置附近区域压缩效应显著减弱，圆盘中心纵向压缩效应呈现不断减小趋势，表明超高温作用可以削弱岩石纵向的应力分布。图 3花岗岩巴西劈裂圆盘纵向应力分布Figure 3The vertical stress distribution of granite Brazilian disk在花岗岩圆盘巴西劈裂过程中，岩心各处沿 x 轴(横向)和 y 轴(纵向)应力分布的差异明显，故作 xy平面内剪切应力图谱分析如图 4 所示。若统一剪切应力图谱显示图例为300300 MPa，巴西劈裂圆盘载荷施加处剪切效应显著，其他位置(如圆心附近)剪切力可以忽略。从局部看，左部圆盘和右部圆盘靠近同一载荷施加位置分别呈现拉伸效应与压缩效应，其数值大小相等且较 y 轴呈对称分布;上部圆盘和下部圆盘靠近两端载荷施加位置分别呈现拉伸效应与压缩效应，其数值大小相等且较 x 轴呈对称分布;从整体上看，圆盘剪切应力图谱呈中心对称分布。随着岩石处理温度从 25 升高至 650，花岗岩巴西劈裂圆盘各处剪切应力显著降低，表明超高温作用可以削弱岩石在载荷施加位置处的剪切应力分布。图 4t=0025 s 时刻花岗岩巴西劈裂圆盘剪切应力分布Figure 4The shear stress distribution of granite Brazilian disk at t=0025 s22温度载荷联合作用下岩石裂纹延展机制基于超高温作用下岩石应力场分布特征，向花岗岩巴西劈裂圆盘全局网格添加 cohesive 粘结单元并设置粘结拉伸破坏属性。以 25 下巴西劈裂圆盘模91第 6 期赫文豪等:超高温条件下花岗岩力学性质演化规律模拟研究拟数据为例，取第 1(t=0025 s)、第 5(t=0125 s)、第15(t=0375 s)、第 25(t=0625 s)、第 40(t=1000 s)分析间隔数据作为参照，实现巴西劈裂过程中花岗岩模拟裂缝的动态捕捉(图 5)。若统一米塞斯(Mises)应力图谱显示图例为 001000 MPa，由图 5 可知，25处理的花岗岩在 t=0025 s 时刻载荷施加位置出现明显的应力集中现象。由于剪切效应(图 3)的存在，使得载荷施加位置逐渐产生微裂缝，并逐渐向岩石圆盘中部延展，在 0125 s 时在两端载荷位置形成裂纹簇。裂纹簇中沿载荷施加方向径向裂纹具有优势传播方向，在 0375 s 时延展形成裂纹网状结构，进一步延展沟通，在 0625 s 时刻显示为径向宏观裂纹。在载荷继续作用 1000 s 时形成最终网状纹络结构，圆盘中心损伤效果显著。相较于 25 的花岗岩，在同一劈裂时刻下，600 花岗岩受压损伤区域相对较小，裂缝沟通过程中岩石应力释放较为明显。图 525 和 600 下花岗岩巴西劈裂圆盘裂缝延展机制Figure 5The fracture dynamics comparison of granite Brazilian disks between 25 and 600 为进一步揭示温度载荷联合作用下岩石裂纹的延展机制，将不同温度处理下花岗岩巴西劈裂圆盘在 t=1000 s 时形成的最终裂纹汇总于图 6。可以发现，随着温度的升高，花岗岩圆盘沿载荷方向中央裂纹密度显著减小，裂纹沟通难度降低，表明温度升高后岩石横轴方向的应力较小，沿径向张开较为容易，超高温作用可以削弱岩石在横轴方向上的应力分布，与未添加粘结单元时模拟数据支撑结论一致(图 2);同时，花岗岩圆盘横向裂纹延展性显著增强，裂纹沟通较好，表明温度升高后岩石纵轴方向的应力较小，沿横向张开较为容易，超高温作用可以削弱岩石在纵轴方向上的应力分布，与未添加粘结单元时的模拟数据支撑结论一致(图 3)。另外，当温度从 25 升高至 650，花岗岩两端载荷位置的裂纹簇规模显著减小，破碎区域较小且裂纹数目减少，表明载荷位置的剪切效应被显著削弱，超高温作用可以削弱岩石在载荷施加位置的剪切应力分布，与未添加粘结单元时的模拟数据支撑结论一致(图 4)。图 6不同温度作用下花岗岩巴西劈裂圆盘裂缝延展机制Figure 6The fracture dynamics of granite Brazilian disks under various temperature conditions02华 南 师 范 大 学 学 报(自 然 科 学 版)第 54 卷在巴西劈裂过程中，花岗岩圆盘应力分布特征与裂纹延展机制互为印证，表明数值模型的有效性与准确性，该模型可被用于模拟花岗岩圆盘在超高温作用下的损伤破坏效果，因此可计算其在不同温度条件下的巴西劈裂强度。根据巴西劈裂数值模型模拟数据，由于超高温作用对花岗岩圆盘横向应力、纵向应力和剪切应力的削弱作用，在温度从 25 升高至 650 的过程中，温度的升高使得岩石沿圆盘中央方向的裂纹沟通更容易，载荷施加位置的裂纹簇规模显著减小，圆盘裂纹密度显著降低，岩石损伤效果更为明显，表明相较于低温度范围处理的岩石，超高温作用下岩石在巴西劈裂过程中受到单轴径向压缩载荷越容易产生拉伸破坏。23温度对岩石巴西劈裂强度的影响规律根据岩石模拟破裂过程中的载荷曲线，选取最大载荷计算硬质花岗岩在超高温作用下的巴西劈裂抗拉强度，具体计算方法依据国际通用测量标准(ASTM)执行:t=2PDL，(1)式中，P 为施载过程中最大载荷(N)，D 为岩石的直径(mm)，L 为岩石的厚度(mm)。花岗岩巴西劈裂圆盘的破坏过程中受力载荷可以通过收集承压板受到的反作用力获得。如表 2 所示，在温度从 25 升高至 650 的过程中，施载过程中最大载荷不断减小，25 花岗岩在压碎过程中的最大载荷为 62177 N，而 650 高温作用下花岗岩的最大载荷下降至 12271 N，降幅达 80%。该结果表明:在超高温作用下，花岗岩的巴西劈裂抗拉强度显著降低。表 2花岗岩的巴西劈裂圆盘的抗拉强度模拟结果Table 2The numerical results of Brazilian tensile strength of granite discs温度/弹性模量/GPa单轴抗压强度/MPa最大载荷/N巴西劈裂强度/MPa25393090426217713191503226877353213112920029228635494111049300274379364749310084001897837737406794500209968203823581160019876956140662996501616598312271260将巴西圆盘劈裂过程中最大载荷带入式(1)求解可知，在 25、150、200、300、400、500、600、650 温度作用下，花岗岩圆盘破坏过程中受到的巴西劈裂抗拉强度分别为 1319、1129、1049、1008、794、811、298、260 MPa(图 7)。图 7不同温度下花岗岩巴西劈裂圆盘的抗拉强度Figure 7The tensile strength of granite Brazilian disk undervarious temperature conditions随着岩石温度的升高，岩石的劈裂强度不断降低，且存在温度阈值约 500，可使岩石的拉伸强度发生显著变化，岩样强度显著下降。相较于 25 处理岩样的抗拉强度(1319 MPa)，600 处理所得岩样的劈裂强度下降至 260 MPa，降幅达 80%，仅为原有强度的 20%，表明超高温作用下花岗岩巴西劈裂抗拉强度显著降低。3结论为揭示超高温作用下花岗岩破碎特性，利用Abaqus 有限元分析建立了高温花岗岩巴西圆盘劈裂数值模型，探究了硬质花岗岩在温度载荷联合作用下岩石的应力场分布特征、裂纹延展机制和劈裂强度。数值模拟数据结果表明:当花岗岩温度从25 升高至 650，超高温作用可以显著削弱岩石圆盘的应力分布，扰动原有应力场分布特征，尤其是削弱圆盘中心沿横轴方向的拉伸效应、两端载荷施12第 6 期赫文豪等:超高温条件下花岗岩力学性质演化规律模拟研究加位置附近的压缩效应和剪切效应;伴随着岩石应力的削弱，岩石沿圆盘中央方向裂纹沟通更为容易，载荷施加位置的裂纹簇规模显著减小，岩石损伤效果更为明显;根据模拟数据分析，进一步处理载荷数据发现，超高温作用下岩石的劈裂拉伸强度降幅达80%，且存在阈值温度约 500 使得岩石拉伸强度发生显著变化。参考文献:1 林伯强 中国能源发展报告 2018 M 北京:北京大学出版社，2019 2 蔺文静，刘志明，马峰，等 我国陆区干热岩资源潜力估算 J 地球学报，2012，33(5):807811LIN W J，LIU Z M，MA F，et al An estimation of HDresources in Chinas mainland J Acta Geoscientica Sin-ica，2012，33(5):807811 3王贵玲，张薇，梁继运，等 中国地热资源潜力评价 J 地球学报，2017，38(4):449450WANG G L，ZHANG W，LIANG J Y，et al Evaluation ofgeothermal resources potential in China J Acta Geosci-entica Sinica，2017，38(4):449450 4 董晋鹏 高温后不同粒径花岗岩三轴力学特性及颗粒流模拟研究 D 北京:中国矿业大学，2020 5 李天斌，高美奔，陈国庆，等 硬脆性岩石热力损伤本构模型及其初步运用J 岩土工程学报，2017，39(8):14771484LI T B，GAO M B，CHEN G Q，et al A thermalmechani-caldamage constitutive model for hard brittle rocks andits preliminaryapplication J ChineseJournalofGeotechnical Engineering，2017，39(8):14771484 6 杜守继，刘华，职洪涛，等 高温后花岗岩力学性能的试验研究J 岩石力学与工程学报，2004，23(14):23592364DU S J，LIU H，ZHI H T，et al Testing study on mechani-cal properties of posthightemperature granite J Chi-nese Journal of ock Mechanics and Engineering，2004，23(14):23592364 7 ZHANG W，SUN Q，ZHU S，et al Experimental study onmechanical and porous characteristics of limestone affect-ed by high temperatureJ Applied Thermal Engineer-ing，2017，110:356362 8YIN T，LI X，CAO W，et al Effects of thermal treatmenton tensile strength of Laurentian Granite using BraziliantestJ ock Mechanics ock Engineering，2015，48(6):22132223 9CHEN Y L，WANG S，NI J，et al An experimentalstudy of the mechanical properties of granite after hightemperature exposure based on mineral characteristics J Engineering Geology，2017，220:234242 10 KANG F，LI Y，TANG C Grain size heterogeneity controlsstrengthening to weakening of granite over hightempera-ture treatment J International Journal of ock Mechan-ics and Mining Sciences，2021，145:104848/117 11 张宁 高温三轴应力下花岗岩蠕变渗透热破裂规律与地热开采研究 D 太原:太原理工大学，2013 12 程绍桐 加热方式对花岗岩热损伤影响的研究D 北京:北京工业大学，2020 13 李春，胡耀青，张纯旺，等 不同温度循环冷却作用后花岗岩巴西劈裂特征及其物理力学特性演化规律研究J 岩石力学与工程学报，2020，39(9):17971807LI CHUN，HU Y Q，ZHANG C W，et al Brazilian splitcharacteristics and mechanical property evolution of gran-ite after cyclic cooling at different temperatures J Chi-nese Journal of ock Mechanics and Engineering，2020，39(9):17971807 14 余莉，彭海旺，李国伟，等 花岗岩高温水冷循环作用下的试验研究 J 岩土力学，2021，42(4):10251035YU L，PENG H W，LI G W，et al Experimental study ongranite under high temperaturewater cooling cycleJ ock and Soil Mechanics，2021，42(4):10251035 15 LIU S，XU J An experimental study on the physicome-chanical properties of two posthightemperature rocks J Engineering Geology，2015，185:6370 16 MIAO S，PAN P Z，ZHAO X，et al Experimental study ondamage and fracture characteristics of Beishan Granitesubjected to hightemperature treatment with DIC and AEtechniques J ock Mechanics and ock Engineering，2021，54:721743 17 QIN Y，TIAN H，XU N X，et al Physical and mechanicalproperties of Granite after high temperature treatment J ock Mechanics and ock Engineering，2020，53(1):305322 18 郤保平，吴阳春，王帅，等 青海共和盆地花岗岩高温热损伤力学特性试验研究 J 岩石力学与工程学报，2020，39(1):6983XI B P，WU Y C，WANG S，et al Experimental study onmechanical properties of granite taken from Gonghe ba-sin，Qinghai province after high temperature thermal dam-ageJ Chinese Journal of ock Mechanics and Engi-neering，2020，39(1):6983【责任编辑:谭春林责任校对:谭春林英文审校:曾姝倩，邓乾霞】22华 南 师 范 大 学 学 报(自 然 科 学 版)第 54 卷