爆破开挖扰动下锚固节理岩质边坡位移突变特征与能量机理*

杨建华，代金豪，姚 池，胡英国，张小波，周创兵

（1. 南昌大学建筑工程学院，江西 南昌 330031；2. 长江水利委员会长江科学院，湖北 武汉 430010）

中国西南地区大型水电工程高边坡工程规模巨大、地质环境复杂、工程作用强烈，高边坡的变形与稳定控制是高坝工程建设与运行安全的重要保障。钻孔爆破是水电工程岩石高边坡开挖的主要手段，节理岩质边坡在爆破开挖过程中常发生位移突变，影响边坡的稳定。预应力锚索作为一种主动支护技术，因其施工方便、见效快和造价低等优点在边坡支护中得到了广泛的应用，有效地控制了边坡岩体的位移突变。

对于节理岩质边坡的位移突变问题，国内外学者从卸荷松动和爆炸荷载扰动的角度进行了广泛的研究。Yoshida 等、Sheng 等通过数值模拟与现场实测数据的对比分析，认为开挖引起边坡岩体内应力释放及重分布是造成边坡表面岩体位移增加、节理张开及滑移的主要原因。卢文波等的研究表明，爆破开挖过程中岩体初始地应力的释放具有瞬态特性，地应力的瞬态释放会导致岩体节理张开；罗忆等采用室内试验验证了地应力瞬态卸载条件下的节理岩体位移突变问题，并认为突变位移与初始地应力的平方成正比。Gran 等通过对节理岩体模型进行爆破试验，发现爆炸应力波传播至节理处发生反射会将节理拉开；金李等利用一维应力波理论对节理岩体的爆破松动机理进行了分析；代金豪等则从能量的角度研究了节理岩体的爆破松动问题，认为爆炸荷载作用期间岩体积累的应变能及其释放速率决定了岩体的爆破松动位移。以往的研究多将地应力瞬态卸荷和爆炸荷载扰动效应分开讨论，较少涉及二者耦合作用下节理岩质边坡的位移突变特征。

在工程现场，边坡爆破开挖施工与锚固支护常穿插进行。众多学者在预应力锚索对节理岩质边坡变形控制的机理方面开展了研究，如丁秀丽等提出预应力锚索锚固效应使得岩体力学性能劣化在边坡开挖过程中得到缓解，进而抑制了边坡变形；Yang 等通过数值模拟研究发现，锚索预应力的施加可有效改善岩体的应力状态并提高软弱面的抗滑力，从而达到控制边坡变形的目的；李剑等则认为对于高陡岩质边坡，预应力锚索主要在于限制潜在可动滑体在开挖过程中的位移发展。上述有关预应力锚索控制边坡变形的机理研究多针对静态或拟静态力学条件进行。对于深切河谷底部地应力集中区的边坡岩体爆破开挖，存在爆炸荷载与地应力瞬态卸荷2 种动力扰动。有关爆炸荷载与地应力瞬态卸荷耦合作用下的锚固节理岩质边坡位移突变机理，目前研究还较为少见。

本文中，首先，对白鹤滩水电站左岸坝基节理岩质边坡爆破开挖过程中的岩体位移及锚索轴力进行现场监测资料分析；然后，建立三维节理岩质边坡数值模型，分析爆炸荷载与地应力瞬态卸荷耦合作用下节理岩质边坡的位移突变特性及其能量机理，并从能量的角度研究预应力锚索对节理岩质边坡位移突变的控制机理。

1 节理岩质边坡位移突变与预应力锚索轴力监测分析

1.1 工程概况

白鹤滩水电站位于金沙江下游四川省宁南县和云南省巧家县境内，左岸坝基边坡岩体由峨眉山玄武岩构成，其中834～820 m 高程的岩体主要为弱风化下段、弱卸荷的Ⅲ类岩体，820～600 m 高程的岩体主要为弱风化下段～微新无卸荷的Ⅲ类岩体，600 m 高程以下的岩体主要是微新无卸荷的Ⅱ类岩体，如图1 所示。680 m 高程出露的NE 向断层F斜切坝基边坡，在F上盘岩体中发育有NNW 向张扭性断层f和f、层间错动带C和C以及层内错动带LS、LS和LS等。在F下盘岩体中存在有顺坡向的层内错动带LS、LS和LS等，还发育有纵横交错且易松弛的柱状节理和微裂隙。

图1 白鹤滩水电站左岸坝基边坡地质剖面图Fig. 1 Geological section of the left-bank dam foundation slope at the Baihetan hydropower station

1.2 监测点布置

随着左岸坝基边坡爆破开挖施工向河谷底部推进，LS、LS、LS等顺坡向层内错动带和柱状节理岩体在620～580 m 高程逐渐出露，这些软弱结构面交错切割岩体形成的潜在可滑动块体与临空面构成了不利组合，严重威胁边坡的稳定性。为此，现场垂直于建基面布置长短不等的预应力锚索对坝基边坡潜在不稳定块体进行加固，锚索预应力等级为3000 kN。

为掌握开挖施工过程中坝基边坡的位移变化情况及预应力锚索的工作性态，现场开展了边坡变形与预应力锚索轴力的联合监测，如图2 所示。在615 m 高程处布置了一个表面位移计，编号为D1；考虑到层内错动带LS在开挖过程中可能会产生过大的滑动而引起边坡失稳，在1 号帷幕洞与1 号排水洞的侧壁上分别布置了一套裂缝计监测裂缝的开合度和错动位移，编号分别为F1 和F2；在607 和600 m 高程处3000 kN级的预应力锚索上分别安装了一套锚索测力计，编号分别为C1 和C2。

图2 边坡位移与预应力锚索轴力监测点布置Fig. 2 Arrangement of the observation points for monitoring the slope displacement and anchor cable axial force

1.3 边坡位移及锚索轴力监测数据分析

2016 年6 月1 日至2016 年8 月23 日，白鹤滩水电站左岸坝基边坡从605 m 高程开挖至580 m 高程，607 和600 m 高程处3000 kN 级预应力锚索分别在2016 年6 月1 日和7 月3 日前安装完成。该段时间内边坡表观测点D1 处的水平位移、内部裂缝监测点F1 处的错动位移和F2 处的张开位移随时间变化如图3 所示。从图3 可以看到，爆破开挖施工期间，各监测点的位移增长较快，特别是在595～590 m高程和590～585 m 高程梯段爆破时，由于爆区距F1 和F2 监测点较近，F1 监测点的错动位移和F2 监测点的张开位移均有在爆破后大幅突然增加的现象，即位移突变。与此同时，边坡表面监测点D1 也发生了多次较明显的位移突变。尽管位移监测点所在高程的预应力锚索已经安装完成，但位移突变现象仍有发生。2016 年8 月23 日梯段爆破施工停止后，各监测点的位移增长较为缓慢，且逐渐趋于稳定。截至2016 年10 月29 日，边坡表面测点D1 处的水平位移为29.9 mm，内部监测点F1 处的错动位移为7.0 mm、F2 处的张开位移为6.8 mm。

图3 各监测点位移随时间变化曲线Fig. 3 Slope displacement histories at the different observation points

锚索监测点C1 和C2 处的轴力变化如图4 所示，在595～590 m 高程梯段爆破开挖期间，607 m 高程处C1 监测点因距爆区较近且爆破规模较大，该监测点的锚索轴力出现了明显的突增现象，增量为54.6 kN，在随后的2 个台阶爆破开挖过程中，C1 监测点的锚索轴力也发生了不同程度的突增。位于600 m 高程的C2 监测点在590～580 m 高程岩体爆破开挖过程中发生了较大幅度的锚索轴力突增，最大增量为53 kN。待开挖至580 m 高程，爆破施工停止后，预应力锚索轴力增大变缓，轴力逐渐趋于稳定。C1 和C2 监测点锚索的轴力最终分别稳定在3383.9 和3321.2 kN，较初始值分别增大16.2%和8.9%。

图4 各监测点预应力锚索轴力随时间变化曲线Fig. 4 Axial force histories of the prestressed anchor cables at the different observation points

综上，节理岩质边坡位移的变化与锚索轴力的变化具有一致性，二者在附近有大规模岩体爆破时均出现了不同程度的突变，预应力锚索通过自身轴力的增大抑制边坡出现过大的位移突变。毫无疑问，岩体中结构面交错分布形成的具有后缘拉裂面、底滑面和侧滑面的潜在滑动块体是边坡发生位移突变的内在条件，而爆破开挖扰动是诱发节理岩质边坡位移突变的外在条件。众所周知，深切河谷具有应力降低、升高和原始三带的“驼峰形”及谷底存在“高应力包”的地应力分布基本规律，“高应力包”深度可达谷底以下150～200 m。刘国锋等的研究表明，白鹤滩水电站河床谷底应力在20 MPa 以上。对河床谷底高地应力岩体的爆破开挖，伴随着炸药爆轰、岩体破碎及新开挖面的形成，开挖边界上的岩体地应力在岩体破碎瞬间也随之突然释放；该过程为一瞬态卸荷力学过程，在边坡保留岩体中激发瞬态卸载应力波，引起近邻开挖面的岩体应力动态调整。因此，对于谷底高地应力岩体爆破开挖，诱发节理岩质边坡位移突变的动力扰动除爆炸荷载外，还包括地应力瞬态卸荷过程，二者几乎同步发生。然而，由于监测设备的局限性，上述实测数据仅能反映每日的边坡位移及锚索轴力变化值，无法反映爆破开挖动力扰动下边坡位移与锚索轴力的瞬态调整过程。下文将采用数值模拟方法，研究爆炸荷载与地应力瞬态卸荷联合作用下边坡位移与锚索轴力的动态变化过程及相关能量机理。

2 数值模拟

2.1 数值模型与材料参数

采用有限差分软件FLAC，研究爆破开挖动力扰动下锚固节理岩质边坡的位移突变特性及能量机制。白鹤滩左岸坝基边坡地质条件复杂，岩体结构面众多，影响边坡位移突变的因素多而复杂，若将众多因素均考虑在内，很难揭示爆破开挖扰动下锚固节理岩质边坡位移突变的力学本质。为此，采用概化的节理岩质边坡模型进行数值模拟研究，以揭示爆破开挖扰动下锚固节理岩质边坡位移突变的能量机理。为便于进行机理分析，考虑爆炸应力波和瞬态卸载应力波垂直入射节理面的情况，建立含竖直和水平节理的概化边坡模型，如图5 所示。模型整体尺寸为40 m×6 m×30 m，包括3 个高5 m 的直立台阶，上2 层台阶已开挖完成，第3 层台阶正在进行爆破开挖施工。第2 层和第3 层台阶岩体中存在一个由一组后缘竖直节理面、一组底部水平节理面和两组侧部竖直节理面切割而成的潜在滑动体，其尺寸为8 m×2 m×10 m。在第2 层台阶上施加2 根垂直于边坡坡面的预应力锚索，锚索总长设置为20 m，自由段与锚固段的长度分别为12 和8 m，2 根锚索到第3 层台阶顶部的距离分别为2 和4 m。卢文波等、罗忆等和金李等也采用类似的概化模型进行了相关的研究，并取得了有益的结论；此外，卢文波等的研究还表明，如底滑面为顺层的倾斜滑面时，节理岩体的位移突变特征与底滑面为水平时的情况类似，只是最终的位移值会更大一些。篇幅所限，本文中仅考虑底滑面为水平节理面的情况。

图5 锚固节理岩质边坡数值模型（单位：m）Fig. 5 The anchored jointed rock slope model (unit: m)

白鹤滩左岸坝肩及坝基边坡主要出露玄武岩中部Pβ及Pβ层，总厚度约350 m，岩性主要有斜斑玄武岩、隐晶质玄武岩、柱状节理玄武岩、杏仁状玄武岩、角砾熔岩、凝灰岩6 种。高程660～555 m 坝基部位出露 Pβ层第1 类柱状节理玄武岩，该层厚度约为55 m。本文中研究610～580 m 高程坝基边坡岩体开挖，该高程范围内坝基边坡岩体只含有第1 类柱状节理玄武岩这一种岩层，根据现场地质勘察和室内试验结果，第1 类柱状节理玄武岩岩石及节理物理力学参数如表1 所示。数值模拟中将岩石假定为线弹性材料；对于岩体中的节理面，采用FLAC提供的基于库伦剪切本构模型的无厚度interface 单元进行模拟，interface 单元可以模拟荷载作用下2 个接触面之间的错动滑移、张开、闭合等变形，模拟中不考虑节理面的抗拉强度和黏聚力；预应力锚索则采用cable 单元进行模拟，锚索的物理力学参数如表2 所示。根据Blair 的研究，为满足动力分析精度，数值模型的网格尺寸不应超过应力波波长的1/12～1/6，本模型的最小波长(Rayleigh 波的波长)约为2.2 m，因此动力加载边界附近的单元尺寸划分为0.25 m，模型总共划分为1 793 136 个单元。

表1 岩石与节理的物理力学参数Table 1 Physical and mechanical parameters of the rock and joints

表2 锚索的物理力学参数Table 2 Physical and mechanical parameters of the anchor cable

2.2 动力荷载

由于预裂孔距边坡开挖轮廓最近，且预裂孔爆破后形成的预裂缝阻碍了主爆孔爆破产生的应力波向边坡岩体内传播，预裂孔爆破对边坡保留岩体的动力扰动往往是最强烈的。因此，本文中仅研究预裂孔爆破动力扰动作用下锚固节理岩质边坡的位移突变。预裂爆破时每段4 个炮孔同时起爆，由于炮孔直径远小于数值模型的尺寸，若将所有的炮孔均建立在模型中，为保证计算精度，数值模型的单元数将达到上千万，从而导致计算效率低甚至无法计算。为提高计算效率，本文中将炮孔壁受到的爆炸荷载压力等效施加在图5(a)所示的爆破开挖面上。根据圣维南原理，这种等效处理方法在计算炮孔附近的岩体动力响应时会产生一定的偏差，而在计算炮孔中远区的岩体动力响应时则与爆炸荷载压力施加在炮孔壁上的计算结果基本一致，该方法对于模拟中远区岩体动力响应问题的有效性已被证实。本文中研究的竖直节理面距炮孔较远，属于中远区问题，采用该等效方法是可行的。对于预裂孔爆破，炮孔中心连线与炮孔轴线所组成的平面即为爆破开挖面。根据圣维南原理，爆破开挖面上的等效爆炸荷载压力()与炮孔壁上的爆炸荷载压力()满足如下关系：

式中：ρ为炸药密度；为炸药爆轰速度；γ 为炸药的等熵指数，一般取3；υ 为爆生气体的绝热膨胀常数，近似取1.5；为装药直径。

根据白鹤滩水电站左岸坝基边坡现场爆破采用的钻孔及炸药设计，取ρ=1300 kg/m,=4500 m/s,=76 mm,=25 mm,=0.6 m，炮孔装药长度=4.5 m。由式(1)和(2)计算可得爆破开挖面上的等效爆炸荷载峰值压力=15 MPa。目前获取爆炸荷载压力时程曲线的方法主要有3 种：(1)半经验半理论的计算公式；(2)通过简化模型分析爆生气体的运动过程；(3)采用数值模拟软件（如LS-DYNA）模拟炸药爆轰过程。由于炸药爆炸产生的高温高压环境，现有的测试手段还很难从现场准确获取炮孔壁上或爆破开挖面上的爆炸荷载压力时程曲线。在实际工程爆破的数值模拟中，目前大多采用半经验半理论的爆炸荷载压力时程曲线，其中以双指数函数和进一步简化的三角形函数压力时程曲线应用最为广泛。本文中采用三角形荷载曲线，如图6 所示，取爆炸荷载压力上升时间=1 ms、持续时间=9 ms，其计算过程详见文献[16] 和[23]，爆炸荷载压力时程曲线()的表达式为：

图6 爆破开挖面上等效爆炸荷载与地应力瞬态卸荷历程曲线Fig. 6 Time-histories of the equivalent blasting pressure and the transient unloading of in-situ stress on the blasting excavation boundary

同一排炮孔同时起爆时，爆生裂纹优先在炮孔连线方向上扩展，当裂纹完全贯通、开挖面上的爆炸荷载压力衰减至与开挖面上的地应力大小相等时，开挖面岩体才出现宏观上的卸荷效应。根据应力连续条件，爆破开挖面上的地应力瞬态卸荷历程在岩体出现宏观卸荷后与爆炸荷载作用曲线重合，如图6所示，地应力瞬态卸荷时程曲线σ()可表示为：

式中：σ为爆破开挖面上的地应力；为卸荷开始时间。

第3 层最后一个梯段预裂爆破时，取爆破开挖面上的地应力σ=5 MPa，根据图6，地应力瞬态卸荷持续时间=-=2.7 ms。需要说明的是，虽然白鹤滩水电站河床谷底的初始地应力较高（20 MPa 以上），但随着开挖工作面逐步向坝基边坡轮廓面推进，边坡轮廓面上的初始地应力也逐步释放，因此在最后一个梯段爆破时，坝基边坡轮廓面上的地应力取一个较小的值。

数值模拟过程中，首先进行地应力和预应力锚索作用下的静力求解，模型四周设置为法向位移约束，待计算达到平衡后打开动力求解模块，将上述爆炸荷载和地应力瞬态卸荷作用历程用FISH 语言编写后施加在爆破开挖面上进行动力求解。为防止动力荷载激发的应力波传播至模型四周边界发生反射而影响计算结果，动力计算中将模型四周设置为可吸收入射波的黏性边界。岩体阻尼采用FLAC提供的局部阻尼，局部阻尼系数设置为0.157。

3 节理岩质边坡位移突变特征与能量机理

为对比分析预应力锚索对节理岩质边坡位移突变的影响，本文数值模拟分2 类工况进行研究，一类未施加预应力锚索，另一类为施加预应力锚索工况。由于岩体出现朝向河谷的位移时对边坡稳定最为不利，因此本文在潜在滑动块体的前缘和后缘各选取一个监测点和（见图5）分析边坡的水平位移及应变能密度特性。

3.1 未施加预应力锚索工况

未施加预应力锚索时，爆破开挖动力扰动下两监测点的水平位移和应变能密度时程曲线如图7 所示。在初始地应力σ的作用下，节理岩体本身储存有一定量的应变能；爆炸荷载作用后，节理岩体被压缩，岩体中进一步积聚应变能。随着爆炸荷载压力降低以及地应力释放，岩体开始回弹，岩体中的应变能得以释放。在=9 ms 左右时，爆炸荷载压力降为零，被压缩的岩体回弹恢复到爆炸荷载作用前的状态，此时监测点的位移为零。但由于爆炸荷载压力快速降低和地应力瞬态释放所产生的惯性力作用，潜在滑动块体继续向河谷方向运动，导致后缘节理张开，产生了节理张开位移。由于底部及两侧滑面的摩擦作用，在=30 ms 左右时，岩块停止了运动，和两监测点的最终位移分别为6.5 和5.6 mm，该位移即为爆破开挖动力扰动过程中节理岩质边坡的突变位移。爆破开挖完成后，和两监测点的最终水平位移之差为0.9 mm；而爆破前，在初始地应力σ的作用下，、两点的水平位移也正好相差0.9 mm（即该岩块的压缩变形量为0.9 mm），这表明和两点最终位移之差是由于地应力释放、岩体松弛所产生的回弹位移。对于河谷底部高地应力边坡岩体爆破开挖，边坡表面岩体的突变位移包括节理张开位移和回弹位移两部分。对于本算例中边坡表面的监测点，总突变位移为6.5 mm，节理张开位移为5.6 mm，回弹位移为0.9 mm，节理张开位移占到了总突变位移的86.2%，是节理岩质边坡位移突变的主要组成部分。

图7 无预应力锚索时各监测点的岩体位移和应变能密度时程曲线Fig. 7 Rock mass displacement and strain energy density histories at the observation points without prestressed anchor cables

从能量的观点来看，岩体由于地应力作用本身所存储的应变能以及爆炸荷载挤压作用所集聚的应变能为节理边坡岩体位移突变提供了能量来源，上述2 种应变能的快速释放是节理岩体位移突的直接诱因。代金豪等的相关研究表明，岩体存储和集聚的应变能越高、应变能释放持续时间越短，产生的位移突变就越大，即位移突变与应变能释放率正相关。应变能释放率为：

因此，在保证开挖岩体正常破碎的前提下，通过优化爆破方案，降低岩体应变能的释放率有助于减小节理岩质边坡的位移突变。如采用低密度、低爆轰波速的炸药降低爆炸荷载峰值压力，从而降低岩体集聚的应变能；采用小抵抗线降低开挖面上的地应力，从而降低岩体存储的应变能；加强炮孔堵塞减缓爆生气体从孔口逸出，从而延长爆炸荷载压力衰减时间，即延长应变能释放持续时间。

3.2 施加预应力等级锚索工况

除爆破方案优化外，工程中还主要采用锚杆、锚索等支护方式来控制节理岩质边坡的位移突变。为探讨预应力锚索对节理岩质边坡位移突变的控制机理，设置锚索预应力等级分别为1 000、1 500、2 000、2 500 和3 000 kN。图8 给出了不同等级预应力锚索作用下边坡表面点的位移和应变能密度变化时程曲线。可以看到，不同等级预应力锚索作用下，节理岩体的位移总体特征基本相同，均具有压缩→回弹→突变→稳定这4 个阶段。不同预应力等级下节理岩体的压缩与回弹路径完全重合，这表明预应力锚索并不影响节理张开前岩体的位移与能量转化，不同等级预应力锚索作用下点的应变能密度时程曲线完全重合也证明了这一点。在=9 ms 节理张开发生位移突变后，锚索预应力等级越高，节理岩体的突变位移越小。当预应力超过1 000 kN 时，节理岩体在达到最大突变位移后会出现明显的“被拉回”现象，最终使突变位移稳定在较小的值。当预应力等级超过1 500 kN 时，随着锚索预应力等级的提高，节理岩体的最终突变位移虽有所减小，但减小幅度不大。图9 给出了不同等级预应力锚索作用下监测点的节理最终张开位移和岩体最终回弹位移，随着预应力等级的提高，节理张开位移同样呈现先显著后缓慢降低的趋势，而岩体回弹位移基本保持不变，这表明预应力锚索主要是通过抑制节理张开来降低边坡岩体的位移突变。结合前述实测资料分析结果，预应力锚索的施加虽不能完全避免节理岩质边坡在爆破开挖过程中产生位移突变，但采用适当预应力等级的锚索进行加固后可有效减少岩体内节理的张开位移，从而防止边坡岩体位移突变过大。

图8 不同等级预应力锚索作用下监测点A 的位移和应变能密度的时程曲线Fig. 8 Rock mass displacement and strain energy density histories at the monitoring point A under different levels of the prestressed anchor cables

图9 不同预应力等级下的节理最终张开位移与岩体最终回弹位移Fig. 9 Joint opening displacement and rebound displacement under different levels of the prestressed anchor cables

针对预应力锚索在控制节理岩质边坡位移突变所发挥的作用，从能量的角度对其作用机理开展讨论。对于预应力锚索，其总的应变能由锚索索体的拉压弹性应变能和锚固剂的剪切弹性应变能组成，分别由下式计算：

式中：F为第个索体单元的轴力；L为第个索体单元的长度；为锚索索体的弹性模量；为锚索索体横截面面积；τ为第个锚固剂单元的剪切应力；V第个锚固剂单元的体积；为锚固剂剪切模量；和分别为索体单元和锚固剂单元的数量。

通过FLAC中的FISH 语言对锚索总应变能涉及的变量进行实时监测和计算，得到不同等级预应力锚索总应变能的时程变化曲线，如图10所示（篇幅所限，仅以1000 和3000 kN 等级为例）。锚索总应变能在开始突增前经历了降低和升高至初始值这2 个阶段，这与节理岩体在位移突变前发生压缩和回弹相对应；在锚索应变能突增开始后存在吸能和释能2 个阶段，定义锚索的平均吸能速率或释能速率为：

式中：Δt为第个时间段(=1, 2, 3,···,)；ΔW为第个时间段开始与结束时锚索的总应变能之差。锚索平均吸能和释能速率分别表示为和。

结合图8(a)和图10 可以看到，节理岩体发生位移突变后，节理岩体与锚索相互作用，节理岩体的动能转化为锚索的应变能，锚索应变能急剧增加，表现为吸能，此时锚索被拉伸且轴力增大，这与实测资料反映的规律相对应，即节理岩体位移突变与锚索轴力突增具有同步性；锚索的预应力等级越高，吸能速率越高，节理岩体以更高的速率达到较小的突变位移峰值；在节理岩体的位移达到最大值后，锚索的应变能开始释放，释放的应变能作用于节理岩体，部分转化为节理岩体的动能，将滑动的岩块逐渐“拉回”，随着底部及两侧滑面的摩擦耗能，岩块最终停止运动，位移趋于稳定；锚索的预应力等级越高，释能速率也越高，节理岩体“被拉回”的效应也就越明显，即节理岩体以更高的速率稳定至较小的最终位移值。综上，预应力锚索通过自身具有的吸能与释能特性来控制边坡节理岩体的位移突变。

图10 不同等级预应力锚索的总应变能时程变化曲线Fig. 10 Time-histories of the total strain energy of the anchor cables under different prestress levels

受构造应力、边坡结构、开挖方案等多方面影响，爆破开挖面上的初始地应力大小多有不同，此外，不同工程的岩体变形参数也各不相同。不同地应力水平和岩体弹性模量条件下监测点的最终突变位移如图11 所示。从图11(a)可以看到，随着地应力水平的提高，节理岩质边坡的突变位移增大，这是由于较高的地应力会使岩体内储存更多的应变能，相同释放时间情况下具有更高的应变能释放率；当爆破开挖面上的初始地应力由2.5 MPa 升高12.5 MPa 时，与未施加预应力锚索相比，施加3 000 kN 等级的预应力锚索后边坡突变位移分别减小了64.2%和58.7%，控制效果显著。从图11(b)可以看出，随着岩体弹性模量的提高，节理岩质边坡的突变位移逐渐减小，这是由于在相同的应力作用下，较高弹性模量的岩体应变较小，岩体所存储的应变能较低，从而使应变能释放率降低。对于弹性模量为20 GPa 的岩体，未施加预应力锚索和施加3 000 kN 等级的预应力锚索情况下，边坡突变位移由11.2 mm 减小为5.9 mm，减小了47.3%；而对于弹性模量为60 GPa 的岩体，施加3 000 kN 等级的预应力锚索后位移突变降低了59.2%，预应力锚索对岩体弹性模量不同的岩质边坡均体现出了较好的控制效果。由此可见，对于地应力水平较高或岩体弹性模量较低的节理岩质边坡爆破开挖，尤其要注意节理岩体的位移突变问题，可考虑采用适当等级的预应力锚索进行加固。

图11 岩体地应力和弹性模量对不同等级预应力锚索作用下监测点A 的最终突变位移的影响Fig. 11 Effects of in-situ stress level and elastic modulus of the rock mass on final mutation displacement of the monitoring point A under different levels of the prestressed anchor cables

4 结 论

（1）白鹤滩水电站左岸坝基锚固节理岩质边坡在爆破开挖过程中出现了明显的岩体位移突变和锚索轴力突增的现象，二者具有同步性。

（2）对于深切河谷底部高地应力边坡岩体爆破开挖，爆炸荷载对岩体挤压作用所集聚的应变能以及地应力作用下岩体所存储的应变能为节理岩质边坡位移突变提供了能量来源，爆破过程中这2 种应变能的快速释放是节理岩体位移突变的直接诱因。

（3）深切河谷底部高地应力边坡岩体在爆破开挖过程中产生的突变位移包括节理张开位移和岩体回弹位移，爆破开挖面上初始地应力水平越高、岩体弹性模量越低，产生的突变位移越大；预应力锚索主要控制节理张开位移，锚索预应力等级越高，其吸能和释能速率越高，位移突变控制效果越明显，但当锚索的预应力等级高到一定程度后，节理岩体的突变位移不再明显减小。

为揭示爆破开挖扰动下锚固节理岩质边坡位移突变的能量机理，本文的数值模拟研究中做了一些简化和假定，岩体结构面实际形态、爆破开挖扰动下岩体及锚固结构的损伤对边坡位移突变的影响还需要开展进一步的研究。