基于水压致裂法的深埋山岭隧道地应力特征分析及应用_王坛华.pdf

DOI:10.13379/j.issn.1003-8825.202210032开放科学（资源服务）标识码（OSID）基于水压致裂法的深埋山岭隧道地应力特征分析及应用王坛华（1.福建省交通规划设计院有限公司，福州350004；2.近海公路建设与养护新材料技术应用交通运输行业研发中心，福州350004）摘要：以莆田至永定高速公路白鸽岭隧道地应力测试为例，介绍地应力测试基本原理和过程。通过分析多年积累的福建地区 16 个隧道实测地应力测试成果，采用回归分析法建立最大和平均侧压力系数与埋深的函数关系，并用于预测隧道地应力水平和判定岩爆等级。研究表明：乘幂函数能较好地反映 600 m 以内埋深地应力变化规律。其成果可供福建地区隧道地应力评价参考。关键词：水压致裂法；地应力；回归分析；侧压力系数；岩爆中图分类号：U452.1文献标志码：A文章编号：1003 8825(2023)03 0165 06 0 引言根据我国现有工程岩体分级理论和工程实践 1，地应力是岩体基本质量指标 BQ 计算的 3 个主要修正系数之一。在地应力场与渗流场耦合作用下，地应力的大小、方向能影响裂隙张开度，进而影响裂隙渗透特性，基于裂隙发育特征的岩体渗透系数可能因地应力场发生弱化或强化，因此地应力场可间接影响岩体基本质量指标 BQ 值的地下水修正系数。因地应力场与断层节理等主要控制性结构面的相互关系，主要控制性结构面对围岩稳定性的影响可能加剧或减轻。地应力直接或间接影响岩体基本质量指标的 3 个修正系数，意味着地应力对围岩分级结果影响可达 1 至 2 级。孙广忠等2认为不仅高地应力状态，低地应力状态对隧道围岩稳定性和灾变性也有重要影响。众多学者在重大公路、铁路选线或工程实践中，对高地应力诱发工程风险进行了深入研究，潘智锋等3对新疆地区某深埋特长11 km 隧道钻孔，进行地应力实测并用于判定岩爆风险；陈仕阔等4研究了新川藏线穿越强构造区隧道工程面临的高地应力问题，建立了高应力背景下的选线模型；张朝强等5以某深埋特长 13 km 隧道为例，研究软岩地应力场与围岩破坏之间的关系，构建了以地应力为核心属性识别模型，成功用于判定岩爆风险等级；魏英华等6以地应力作为边界条件，对隧道岩爆机理进行数值模拟用以推测岩爆发生具体位置。在重大工程理论与实践研究中，地应力都是重要的控制性条件。因此，如何测定和获取地应力值对地下工程岩体分级分类、灾害预防、稳定性分析和预判具有重要意义。天然地应力状态通常是指未经人工扰动，以重力场和构造应力场为主，或有时叠加岩体物理化学变化及岩浆侵入作用等形成的应力，是一种初始应力状态。一般情况下天然应力处于相对稳定的状态，但在地壳升降、断层推挤及地震等活动下，天然应力场会发生大小和方向变化，这种改变极为复杂，后期应力场往往起主导作用。边坡开挖、洞室掘进、爆破等人工活动对天然应力场均造成不同程度扰动，形成二次应力场。福建省地处我国东南沿海，受华夏系、新华夏系构造活动影响强烈。受燕山期以来造山运动影响，侏罗系、白垩系地层广泛分布福建省北部、中部和东部沿海地区。福建地区大规模的火山运动和岩浆侵入活动，深刻改变福建地区地壳结构。地壳抬升剥蚀，以及太平洋板块、菲律宾板块向西挤压俯冲，诸多复杂的地质活动均在不同程度上改变天然地应力场，使天然地应力形成路径更加复杂。福建是多山的省份，山地密度居全国前列，在 收稿日期：2022 11 28基金项目：国家自然科学基金青年项目（42007261）；福建省自然科学基金项目（2021J05026）作者简介：王坛华（1981），男，福建宁德人。高级工程师，博士，主要从事公路工程勘察设计咨询工作。E-mail：。王坛华：基于水压致裂法的深埋山岭隧道地应力特征分析及应用 165 华南地区极具代表性，福建与浙赣粤山区地貌形态极为相似，但构造应力历史有所不同。隧道是福建地区公路常见的构造物类型，为满足长大深埋隧道勘察、设计和施工需要，必须获得隧道工程围岩地质条件，而地应力是最重要的隐性工程条件之一。地应力测试方法有钻孔孔径变形测量法、水压致裂法、钻孔孔壁应变测量法和岩体表面应力恢复测量法等7，但在勘察设计阶段最适用且实测效果较好的是水压致裂法，这也是国际岩石力学学会推荐的方法。本文以莆田至永定高速白鸽岭隧道 ZCS1A地应力测试为例，介绍水压致裂法地应力测试的基本原理与过程。由于勘测阶段地应力测试完全依赖深孔钻探，而深孔钻探风险高难度大，钻孔常因掉块、塌孔等导致试验失败，工程人员常面临无有效地应力数据的问题。通过对分布福建各地的 16 个山岭隧道实测地应力测试数据进行统计分析，寻找地应力发育特征和分布规律，为相似条件下福建地区山岭隧道工程应用提供技术支撑。1 工程概况 1.1 地形地貌与构造条件白鸽岭隧道全长 4.2 km，进出口高程约 295 m，最大埋深约 465 m，隧道区属低山地貌，地形起伏大，天然山坡坡度约 2040，局部达 50，隧址地表沟谷发育，切割较深，大多呈“V”型谷；主山脊走向为 NWW 向，其次以 NE 向为主，深沟四周山脊呈放射状。隧道所在地区是福建省中东部典型的火山地质地貌区，隧道北部、东部及南部分别为侏罗系白垩系火山地貌。隧道周边环状火山断裂与隧道呈较小角度相交，断裂带为其南部火山环状断裂带的一段，为具一定走滑性质的正断层，断层内岩石破碎，节理裂隙发育。隧道以西发育北西向断裂，可见花岗斑岩呈脉状大规模侵入。因此，白鸽岭隧址区地质构造条件较为复杂，地应力场受火山、断层、岩浆侵入、剥蚀等叠加影响。1.2 岩体工程地质性质隧址区地表上覆第四系残坡积黏性土层，层厚小于 4 m；下伏基岩为侏罗系南园组（J3n）凝灰熔岩、熔结凝灰岩，因火山喷发沉积作用，局部夹凝灰质粉砂泥岩。因岩浆侵入活动，局部发育花岗斑岩、辉绿岩脉等。经勘探和试验测试，隧道洞身围岩以微风化凝灰熔岩为主，抗压强度约 74150 MPa，饱和劈裂强度 8.7 MPa，黏聚力 5 MPa，内摩擦角 53。2 水压致裂法地应力测试 2.1 基本原理水压致裂法基于 3 个假定条件：岩体为均质、连续、各向同性的弹性介质；围岩为多孔介质并且岩体孔隙流服从达西定律；岩体中一个主地应力为铅垂方向并且大小等于上覆岩层压力8。设 H和 h分别为钻孔横截面上最大和最小水平主应力，Pw为钻孔承受注液压力，P0为孔隙水压力，t为围岩的抗拉强度，Pb为破裂压力。考虑岩体中的孔隙水压力等作用，钻孔承压段周围岩壁围岩产生破裂的应力条件为PbP0=3(hP0)(HP0)+tK（1）式中：K 为渗透系数，对于完整岩石可取 1。试验中加压至孔壁围岩破裂后，即关闭压裂泵，此时维持裂缝张开的瞬时关闭压力 Ps与裂纹面相垂直的最小水平主应力 h平衡，从而有h=Ps（2）H=3PsPb+tP0（3）当钻孔周边围岩初次破裂以后，对钻孔进行二次加压至破裂缝继续张开，此时的压力为重张压力 Pr。因围岩已破碎，抗拉强度为零，此时Pr=3hHP0（4）代入变换，得到H=3PsPrP0（5）结合试验中测量仪表读数和压力传感器各压裂参数特征值 Pb、Pr、Ps，可以确定钻孔最大主应力 H、最小水平主应力 h大小，同时根据印模器裂纹破裂方向记录可以确定最大水平主应力方向。典型压裂过程曲线，见图 1。15105测试压力/MPa4080120160200测试时间/sPbPrPs图1典型压裂过程曲线 2.2 适用条件根据岩体卸荷分带理论2，目标测试段在一定埋深以后受卸荷作用影响就逐渐减弱，构造应力居路基工程 166 Subgrade Engineering2023 年第 3 期（总第 228 期）于主导地位。因此，福建地区山岭隧道地应力测量段一般选择埋深 150 m 以下至隧道洞身范围。2.3 试验设备与步骤测试主要设备：钻孔承压段封隔系统，由两个串联封隔器组成，在两个封隔器间形成承压段承受液压；液压加压系统，实现对封隔器和钻孔承压段加压的管路系统，以及控制压力液体流向的推拉阀；测量和记录系统，包括压力传感器、流量传感器、压力表等。测试系统示意，见图 2。封隔器推拉阀承压试验段孔壁连接杆高压油泵测量和记录系统电脑传感器液压系统图2测试系统示意 地应力测试的主要步骤有：座封、注水加压、岩壁破裂、关泵、卸压、重张、解封、记录破裂缝方向等。2.4 试验成果针对 ZCS1A 钻孔，选取其中 11 段岩体较完整部位进行试验，剔除无效数据后获得了 8 段试验曲线，从测试曲线上读取各压裂参数特征值，并带入式（1）式（5），得到 ZCS1A 钻孔地应力测试成果，见表 1。各试验段的压力-时间关系曲线，见图 3，均较符合典型压裂过程曲线。表1ZCS1A 钻孔地应力测试成果深度/m抗拉强度/MPaH最大水平主应力/MPah最小水平主应力/MPav垂直应力/MPa最大水平主应力方向171.53.825.594.684.46216.83.567.235.355.64NE47253.12.207.736.106.58298.40.829.437.077.76334.61.3510.147.908.70370.90.8611.078.519.64407.14.8012.449.5910.58NE43425.24.5812.9110.0711.06 ZCS1A 钻孔实测最大、最小水平主应力与自重应力比值（侧压力系数）随岩层深度变化关系曲线，见图 4。测孔的应力值随深度的变化均呈增大趋势，但并不完全呈线性增加。侧压力系数随埋深增大整体呈减小趋势，但变化趋缓。10525 50 75 100 125150175 200测试压力/MPa测试时间/s（a）测深171.5 m105测试压力/MPa20025 50 75 100 125150175测试时间/s（c）测深253.1 m175 200105测试压力/MPa25 50 75 100 125150测试时间/s（d）测深298.4 m175 200105测试压力/MPa25 50 75 100 125 150测试时间/s（f）测深370.9 m15105测试压力/MPa175 20025 50 75 100 125150测试时间/s（g）测深407.1 m15105测试压力/MPa25 50 75 100 125 150175 200测试时间/s（h）测深425.2 m105测试压力/MPa2575125175225测试时间/s（b）测深216.8 m225105测试压力/MPa2575125175测试时间/s（e）测深334.6 m00000000图3ZCS1A 钻孔测试段压力时间曲线 最大水平应力最小水平应力自重应力450400350300250200150468101214应力/MPa埋深/m0.80.91.01.11.21.31.41.5最大侧压力系数侧压力系数最小侧压力系数图4钻孔地应力深度变化曲线 3 福建地区地应力变化规律分析 限于深部岩体地应力测试数据稀缺，以往研究难以获取具有统计意义的地区地应力信息，而将地应力变化和侧压力系数简单概括为线性关系。现将分布于福建省不同地区的 16 座山岭隧道深部地应力实测数据进行统计分析。其山岭隧道地应力测试点，见表 2。各个隧道实测地应力成果，最大主应力与自重应力的比值，即最大侧压力系数 并不是固定的，王坛华：基于水压致裂法的深埋山岭隧道地应力特征分析及应用 167 同一个钻孔在不同深度 值是变化的，最大约2.43，最小 0.66，平均 1.30。最小主应力与自重应力比值也有类似特点，最大约 1.79，最小约 0.61，平均 0.96。可见地应力分布是复杂多变的，大部分情况下第一、第三主应力是近水平的，自重应力近似为第二主应力。也有少部分情况下自重应力成为第一主应力，如宁武高速地坪隧道，分析其原因可能与造山运动、断层活动和沟谷深切侧向卸荷等因素有关。为研究水平地应力与自重应力间的变化规律，考虑到各测段实测点数量差异较大，现将各隧道 值按不同深度分段归类统计分析，即将各实测数据按 50 m 为一档，求其平均值和最大值，获得不同深度条件下的侧压力系数