露天开采是国内矿产资源开采的重要方式（王忠鑫等，2024），爆破作为露天开采的高效开挖手段得到了广泛应用（吴世然等，2021）。但是，爆破产生的振动具有幅值大和频率高等特点（董英健等，2019），其危害控制是矿山安全领域的关键难点（杜军等，2024）。随着开采深度的增加，爆破振动严重影响高陡边坡的稳定性，可能引发滑坡灾害。例如：贵州某铝土矿因爆破振动与雨水侵蚀耦合作用发生崩塌（陈阳等，2024），加里曼丹矿因爆破引发多起边坡破坏事件（Deb et al，2011），四川某矿因地质条件与爆破振动共同作用导致顺层滑坡（胡卸文等，2004）。因此，有必要对爆破振动进行控制及优化（何理等，2024），以保证矿山边坡的安全。
以往研究表明，爆破振动对顺倾边坡稳定性的影响显著且持久，监测显示其影响周期可达10 d，累积变形呈阶梯状加速增长（姜旭桐等，2023）。通过Sadovsky回归方程和数值模拟分析，明确了振动波的传播范围（Su et al，2022）。Li et al（2022）揭示了爆破振动波三维传播规律，提出了小直径缓冲减振爆破技术可有效控制边坡动力响应。基于爆破荷载的尖点突变模型显示，失稳风险与荷载幅值成正相关，与频率成负相关（周子涵等，2020）。随着爆破工程的不断增多，爆破安全预测也日益受到重视（Yin et al，2018），其中伪静态法适用于快速评估，而动力学模型能够更精确地表征动力响应特征（Singh et al，2024）。值得注意的是，机器学习预测模型通过整合多源参数显著提升了预测精度（黄晶柱等，2023）。此外，现场监测与数值模拟协调研究发现，水平径向振动速度对稳定性影响最为显著（王子一等，2023），数值模拟进一步揭示黏聚力是影响稳定性的关键因素，提高内摩擦角，可有效降低失稳风险（林令鑫等，2023）。这些成果为优化爆破设计和边坡防护提供了重要的技术支撑。
综上所述，目前关于爆破振动对矿山边坡稳定性影响及安全预测的研究较为丰富，但存在明显的三维分析不足：首先，现有研究大多基于二维简化模型，难以准确地反映三维边坡的空间动力响应特征；其次，缺乏对单段药量变化引起的边坡三维振动场分布规律的研究。本文基于Midas-GTS有限元软件，对某露天矿山边坡在爆破荷载作用下的三维边坡动力响应特征进行了分析，并结合安全系数及质点峰值振动速度对边坡稳定性进行评价。通过分析不同单段药量对边坡质点振速和稳定性的影响规律，确定了单段最大药量，并探讨了其对边坡稳定性的影响。
1 工程背景
某矿山地形属于中低山构造侵蚀山地，山系属武夷山脉南段东列山的隆升地段，强烈的侵蚀作用造成深沟峡谷纵横，悬崖常高达百米以上，山脉走向NE，区内中部为最高标高（+1 063 m），最大高差为865 m。矿区范围内因海拔高，地形落差大，气候条件较复杂，风力和雨量偏大，夏秋多有灾害性风雨袭击，主峰顶周围为雷击区，夏季常有雷电袭击。矿区岩层较为复杂，包含花岗斑岩、花岗岩、英安玢岩和隐爆角砾岩等，岩石物理力学参数见表1。
表1 岩石物理力学参数Table 1 Rock physical and mechanical parameters |
| 岩石名称 | 弹性模量/GPa | 泊松比 | 容重/（kN·m-3） | 黏聚力/MPa | 内摩擦角/（°） | 抗压强度/MPa |
|---|---|---|---|---|---|---|
| 中细粒花岗岩 | 66.37 | 0.23 | 25.19 | 6.56 | 43.0 | 38.63 |
| 隐爆角砾岩 | 89.33 | 0.20 | 25.88 | 7.82 | 42.1 | 34.08 |
| 花岗斑岩 | 20.17 | 0.33 | 22.74 | 2.01 | 40.4 | 26.83 |
| 地开石化细粒花岗岩 | 15.34 | 0.36 | 22.64 | 2.10 | 39.5 | 23.14 |
2 爆破荷载及模型
2.1 模型建立
本文采用Midas-GTS软件对爆破荷载作用下边坡的动力响应特征及稳定性进行分析，具体岩石物理力学参数见表1。根据剖面图建立三维模型，模型尺寸严格遵守实际工程比例，在网格划分过程中，依据王海兵等（2016）的研究成果，网格尺寸应小于波长的1/16。经计算，本模型爆破主频对应的波长为53.33 m，最大网格尺寸应控制在3.33 m以内，因此本模型单元网格长度设置为2 m，以提高计算精度。最终本模型共建立了59 430个节点，48 075个单元，重力为Z方向（边坡竖直方向）。在岩层接触面处，设置了摩擦接触命令，以准确模拟不同岩层间的力学相互作用。爆破荷载施加于坡脚处的岩体表面，具体位置为高程48.5~75.5 m范围内的坡面区域。该区域对应于实际工程中的钻孔装药段，荷载方向垂直于坡面指向岩体内部。在数值模型中，通过曲面动力模块在该区域施加三角波形的爆破荷载，随后进行弹性边界的施加。边界条件设置考虑了相邻边坡的影响，模型左右两侧约束Y方向（横向）位移，背面约束X方向（水平方向）位移（Chugh，2003）。监测系统沿边坡倾向布置，从坡顶开始每隔一个台阶设置一个监测点，最后一个监测点位于距离爆源30 m处，以确保全面覆盖边坡的关键区域，并准确捕捉爆破振动的空间分布特征，三维边坡有限元模型如图2所示。
2.2 爆破荷载确定
3 初始爆破条件下边坡稳定性分析
3.1 动力响应特征分析
边坡在爆破荷载作用下监测点振速即可表征动力响应特征，图3所示为初始爆破条件下边坡的振速特征。从边坡的X方向振速图［图3（a）］中可以看出曲线波动强烈，表现为快速的振动响应。爆破振动速度峰值最大值约为16.40 cm/s，并在0.1 s左右迅速衰减至0，与施加的三角形爆破荷载的规律吻合。从边坡的Y方向振速图［图3（b）］中可以看出，虽然存在明显的初始波动，但总体振速较低，振动较为均匀且衰减速度较快。从边坡的Z方向振速图［图3（c）］中可以看出，爆破振动速度峰值小于X方向、大于Y方向爆破振动速度峰值，振动响应也较为均匀，衰减速度与Y方向一致。合成振速显示的是各个方向振速的综合响应，最初在极短时间内振速达到峰值，振速峰值为16.89 cm/s，与X方向峰值差值较小，且衰减极为迅速，随着时间的推移振速快速减小，在0.2 s内降低到接近0［图3（d）］。总体表明，在爆破荷载的作用下，边坡动力响应具有明显的方向异性，X方向振动最为显著，尤其在爆破初期和爆源附近区域响应最为强烈。
对监测点振动特征进行分析，边坡动力响应呈现明显的空间衰减特性。距离爆源最近的坡脚区域（测点7~测点9）振动最为剧烈，X方向峰值振速达16.40 cm/s，直接承受爆破冲击作用；随着传播距离的增加，坡中区域（测点4~测点6）X方向振速衰减至1.44~2.88 cm/s，但仍可能引发应力重分布和裂隙扩展；最远的坡顶区域（测点1~测点3）受爆破影响最小，X方向振速维持在1 cm/s以下。这种“近强远弱”的振动传播规律表明，坡脚是爆破荷载的主要作用区，若不加以控制，可能会导致局部崩塌或不稳定，而坡中区域作为振动传播的关键过渡带，其潜在稳定性问题同样不容忽视。
3.2 峰值振速特征分析
峰值振速是爆破工程中的一个关键参数，可作为评价结构安全的指标（李振阳等，2024），初始爆破条件下边坡各监测点的峰值振速如图4所示。由图4可知，各监测点的振速峰值随爆源距离的增加呈衰减趋势，其中1号监测点振速最低，9号监测点振速最大，表明振速峰值与爆源距离之间存在负相关关系，即距离爆源越近，振动影响越显著。通过对比不同方向的振速峰值发现，X方向振速峰值最大，数值为16.40 cm/s；Z方向振速峰值次之，为6.37 cm/s，Y方向振速峰值最小，为2.24 cm/s，表明爆破荷载在X方向的影响最突出，Y方向影响最小，X方向对合成振速起主导作用。上述结果表明，爆破振动能量在传播过程中存在明显的方向性差异，X方向是振动控制的关键方向。
在初始单段药量作用下，振动波频率如图5所示，主振频率为15.56 Hz，参考《爆破安全规程》（GB 6722-2014）可知：当主频处于10 Hz≤f≤50 Hz范围时，永久性岩石高边坡安全允许质点振动速度范围为8~12 cm/s。然而，质点峰值振动速度通常指3个方向中的最大分量值，所在方向被称为主方向，X方向峰值振速最大，因此主方向为X方向，主方向峰值振速为16.40 cm/s，超出安全允许范围，为确保边坡安全，需降低振速峰值。
3.3 边坡稳定性分析
边坡位移和塑性区是评估边坡稳定性的关键指标，前者反映爆破荷载下的变形情况，后者揭示潜在滑坡位置。图6所示为初始爆破条件下边坡位移云图，从图中可以看出，在爆破荷载作用下X方向位移最大值为2.94 cm，Y方向位移最大值为4.43×10-2 cm，Z方向位移最大值为4.61×10-1 cm，总位移最大值为3.43 cm。其中，X方向位移在3个方向中最大，表明边坡变形主要集中在X方向上，即爆破荷载对水平方向的影响最大。这是由于爆破荷载施加方向为X反方向，直接激发岩体的水平振动，根据应力波传播理论，P波在自由面反射后产生拉伸波，导致水平方向质点位移幅值增大。
4 单段药量优化及降振效果分析
4.1 单段药量的优化
基于上述分析，为确保边坡安全，需降低振速峰值并提高稳定性。单段药量是影响爆破振速的主要因素之一，通过优化单段药量，可以有效控制振速峰值，从而提升边坡稳定性。
根据萨道夫斯基公式（李启月等，2024）反推出最小单段药量公式，表示为
式中：Q min为最小单段药量；V为允许振速；k为场地系数；α为衰减系数。
根据《爆破安全规程》（GB6722-2014）规定，振速V取最小值8 cm/s；k为场地系数，硬岩取50~150，根据场地实际情况取150；α取1.5，此处爆心距R为测点9的爆心距，取30 m。经过计算得出最小单段药量Q=74.84 kg，因此，为优化单段药量，分别选取75，100，125，150，175，200，225 kg进行三维数值模拟计算，在保持其他爆破参数不变的情况下，分析单段药量对边坡动态响应特性及稳定性的影响。由于测点9的振速在所有测点中最大，因此重点分析不同单段药量下测点9的振速变化规律。计算结果如图8所示，从图中可以看到，不同单段药量引起的振速曲线变化趋势相似，但峰值振速大小有所差异。随着单段药量的增大，各方向振速峰值逐渐增大，其中，75 kg振速峰值最低，225 kg振速峰值最高，表明单段药量的增加导致爆破荷载强度增大，进而使边坡振动响应幅度显著增强。
表3 不同单段药量下的接收波主频Table 3 Dominant frequency of received wave under different single-segment dosages |
| 药量/kg | 主频 | 药量/kg | 主频 |
|---|---|---|---|
| 75 | 17.6 | 175 | 13.4 |
| 100 | 14.5 | 200 | 16.4 |
| 125 | 15.7 | 225 | 14.3 |
| 150 | 20.0 |
基于振速分析结果，单段药量为200 kg时振动速度满足《爆破安全规程》（GB6722-2014）的允许阈值。为进一步确定单段最大允许药量，本研究对不同单段药量作用下边坡安全系数变化规律进行了定量分析，如图10所示。由图10可知，单段药量从75 kg增大至100 kg，安全系数降低了3.4%；单段药量从100 kg增大至125 kg，安全系数降低了2.0%；单段药量从125 kg增大至150 kg，安全系数降低了3.5%；单段药量从150 kg增大至175 kg，安全系数降低了2.7%；单段药量从175 kg增大至200 kg，安全系数降低了1.9%；单段药量从200 kg增大至255 kg，安全系数降低了2.7%。其中，单段药量在175~200 kg范围内，安全系数降幅最为平缓，表明该区间具有最佳的安全—效益平衡特性。基于前文所述，在爆破荷载作用下边坡安全系数应处于1.18~1.23，计算结果表明，单段药量在75~200 kg范围内，边坡安全系数均满足规范要求，而单段药量为225 kg时，边坡安全系数为1.17，处于欠稳状态。从边坡安全角度考虑，选择200 kg作为单段最大药量。
4.2 优化后边坡降振效果分析
基于200 kg单段药量作用下边坡的振速特征（图11）和200 kg单段药量作用下边坡的各监测点峰值振速（图12），通过分析不同方向的振速变化，得出以下规律特征：X方向峰值振速为10.60 cm/s，较250 kg单段药量爆破荷载作用下边坡呈现显著降低，且振速衰减曲线与初始单段药量工况保持高度相似，这是由于爆破荷载波形虽幅值增大，但仍保持三角波特征，未改变荷载作用时程特性；Y方向振速峰值仅降低至1.61 cm/s，降低比例明显低于X方向，表明该方向对药量变化的敏感性较弱；Z方向振速则表现出中等敏感性，峰值振速降低至4.46 cm/s，其空间分布呈现典型的三维衰减特征。从总振速矢量合成结果来看，最大合成振速与X方向振速峰值差值较小，这一现象从能量分布角度印证了边坡动力响应的各向异性特征：爆破荷载产生的振动能量主要沿X方向传播，这与边坡岩体结构面的优势方位及爆破源—坡体空间几何关系密切相关，进一步验证了X方向作为爆破振动主效应方向的控制地位。
图11 200 kg单段药量作用下边坡的振速特征Fig. 11 Vibration velocity characteristics of slope under the action of 200 kg single-stage charge |
对200 kg单段药量作用下边坡稳定性进行了分析，图13为200 kg单段药量作用下边坡的位移云图。由图13可知，相较于初始爆破条件下边坡的位移（图8），200 kg单段药量引起的位移场范围显著减小。X方向位移影响范围从坡脚延伸至380 m处，最大位移值为1.04 cm；Y方向位移影响范围延伸至330 m处，最大位移值为0.37 cm；Z方向位移影响范围变化较小，但最大位移值减小为0.18 cm。总位移影响范围扩展至420 m，最大位移值为1.12 cm。通过对比3个方向的位移分布特征可以发现，X方向位移值最大，表明边坡在X方向的振动效应最为显著，这一现象与振速分布规律具有一致性，进一步验证了X方向作为爆破振动主效应方向的控制地位。图14展示了200 kg单段药量作用下边坡的塑性区变形，与初始爆破条件下边坡的塑性区云图（图8）相比，此时边坡塑性区范围显著减小，塑性区仍然贯通坡脚处。同时，边坡安全系数提升至1.20，符合规范要求，整体处于稳定状态。
5 结论
对初始爆破条件下三维边坡的动力响应特征及稳定性进行了初步评价，开展了爆破单段药量的优化研究，最后对单段最大药量条件下三维边坡进行了稳定性分析。得出如下主要结论：
（1）在初始爆破条件下，三维边坡水平方向（X方向）的振动响应最为显著，主方向的质点峰值振速为16.40 cm/s，超出安全允许范围（8~12 cm/s），边坡的安全系数为1.13，低于规范要求的1.18，处于局部失稳状态。
（2）随着单段药量的增加，振速峰值不断增加，边坡安全系数逐渐降低。单段药量从75 kg增加至225 kg，边坡的主方向质点峰值振速从2.64 cm/s升高至12.51 cm/s，安全系数从1.38降低至1.17，最终确定单段最大药量为200 kg。
（3）当单段药量为200 kg时，边坡主方向峰值振速为10.60 cm/s，安全系数为1.20，处于规范要求范围，边坡处于稳定状态。结合边坡安全系数和质点峰值振动速度，对边坡稳定性进行了更准确的评价。
甘公网安备 62010202000672号
