深井冻结黏土爆破振动传播规律研究

于建新，郭 敏，高帅杰，王 晓，董保立

（1.河南理工大学 土木工程学院，河南 焦作 454000；2.河南迅达爆破有限公司，河南 焦作 454000）

0 引 言

借鉴立井冻结基岩爆破的成功经验，深部立井冻结黏土爆破施工在矿井建设中逐步使用。 冻结黏土爆破主要面临2 个问题：①冻结黏土的密实坚硬程度不如岩石，导致土体爆破效率低。 ②冻土爆破装药量要合理的控制，防止造成冻结管和立井井壁的损伤。 关于冻土爆破的研究，早期主要研究青藏高寒区的冻土爆破，张俊兵等［1］依托青藏高原多年冻土区进行了冻土爆破漏斗试验，确定了高含冰量冻土的爆破参数。 随着人工冻结法的应用，宗琦等［2］最早提出一种宽孔距崩落技术，该技术主要控制炮眼密集系数来提高爆破掘进效率。 文献［3－4］根据理论研究与现场试验，总结了多年冻土和季节冻土的爆破方法。 李志敏等［5］探究了冻结砂土的爆破机理，将爆破作用区域划分为空腔、挤压区、破碎区、裂纹区和弹性振动区，细化了冻结砂土爆破破坏形式。 李涛等［6］利用TDW－200 冻土三轴试验机对饱和冻结灰砂岩进行了试验，获得了不同围岩条件下冻结灰砂岩弹性模量的演化规律。 张慧梅等［7］研究了冻融循环下岩石的物理力学性质，研究发现：随着冻融循环次数增多和围压增加，岩石破坏形式逐渐由脆性转化为延性。 李栋伟等［8］通过构建人工冻土BP 神经网络本构模型，较好地模拟了冻土中复杂应力路径变形特征。 但对冻土中爆破应力波的振动传播规律及能量特征，没有更深一步地进行研究。 因此开展冻结黏土爆破振动监测，研究冻结黏土中爆破应力波的传播规律和爆破能量分布，进行爆破方案的合理优化，对于立井冻结爆破施工具有重要意义。 随着小波（包）分析在工程中的应用，林大超等［9］研究发现小波分析不仅能给出爆破振动能量在不同频率范围内的相对分布规律，还可以给出不同频率带上振动分量的分布和实际衰减规律，在工程中的应用越发广泛。 针对不同的工程进行小波（包）能量分析，发现不同的爆破介质能量的衰减规律是不同的。 例如在露天采矿场边坡［10］、地 下 工 程 爆 破［11－12］、采 石 场［13］等 进 行 小波（包）分析，得出的能量衰减规律呈现出差异性。而目前对冻结黏土单一介质的爆破能量分析还鲜见研究，大多是针对砂土、岩石等类型的研究。 其中仅有单仁亮［14］等结合立井模型试验，对冻结岩壁的能量衰减进行了预测，与冻结黏土爆破较为接近。

针对立井冻土单一介质爆破振动传播规律的研究较少，尤其是应力波在冻结黏土中的传播。 以赵固二矿西风井冻结爆破掘进为背景，在深690 m 处开展冻结黏土爆破试验，分析爆破振动在冻土中的传播规律，并结合小波（包）分析得到振动时能量分布特征，以期指导深井冻结黏土爆破的工程实践，保证周围冻结管和已浇筑井壁的安全。

1 试验方案

根据爆破漏斗理论，在立井冻结黏土段设计单孔爆破漏斗试验。 该土层冻结温度－15 ℃，极限强度为8.16 MPa，黏聚力16.6～214.1 kPa。

试验炮孔孔径设计为50 mm，采用水胶炸药，药卷长度320 mm，每孔半卷炸药。 相邻炮孔设计间距不小于1.8 m，炮孔深度按0.05 m 逐渐增大，最浅炮孔深度为0.35 m，最深炮孔深度为0.50 m，共布置了4 个炮孔（Ⅰ、Ⅱ、Ⅲ、Ⅳ）、5 台振动传感器。 炮孔编号1～4，传感仪依次为A、B、C、D、E 具体布置如图1所示。

图1 炮孔及传感仪布置示意Fig.1 Schematic of holes and sensor arrangement

2 振动测试及分析

2.1 测试参量

根据炸药的爆破作用理论和大量的现场爆破表明，炸药爆破产生冲击波和地震波，二者引起的振动，是对周围岩土体产生破坏作用的重要因素。 而振动速度的大小与其破坏程度是成正相关的，常把爆破振动速度作为当前国内外爆破振动效应和爆破安全控制监测的标准。

结合目前国内爆破振动标准，本次试验将垂直、水平径向和水平切向3 个速度分量作为试验参量。

2.2 现场测点布置

测试前，将爆破振动传感器安装在选定的测点位置，第1 个测点布置在距离爆破孔2.5 m 的位置，共布置了5 台测振仪，每个测点相距125 cm，最后一个测点距离井壁内缘1.8 m，仪器测点与爆点径直方向为x，切向方向为y，垂直方向为z，现场布置如图2 所示。

图2 传感器的布置Fig.2 Arrangement of measuring sensors

2.3 冻土爆破振动结果的计算与分析

冻结黏土的松动爆破效果如图3 所示。

图3 爆破效果Fig.3 Blasting effect

4 个炮孔中，Ⅰ号炮孔爆破引起的振动速度最大，且振动具有规律性；Ⅱ号和Ⅳ号炮孔振动速度小，且振动速度无规律；Ⅲ号炮孔的垂直速度曲线跳跃较大，无规律性。 因此，把Ⅰ号炮孔作为典型爆源，对各测点的监测结果进行分析，其中测距R利用余弦公式计算得出，三向振动速度取振动时程中的最大值，监测结果见表1。

表1 监测数据Table 1 Monitoring data

从表1 中看，各向振动速度呈现出差异，垂向方向振动速度＞径向方向振动速度＞切向方向振动速度，说明了冻土爆破后应力波的传播主要是沿某个特定方向。 单仁亮［15］在对冻结岩壁的振速研究中，也得出岩体振动具有明显方向差异性，说明了冻土和岩石爆破的振速传播均具有方向的差异性。 对于爆破振动速度v，我国普遍使用萨道夫斯基经验公式［16］：

式中：K为场地系数；α为衰减系数；Q为单段最大装药量，kg；R为测点与爆破位置距离，m。

表2 K 和α 拟合结果Table 2 Fitting results of K and α

图4 各分量及合速度拟合曲线Fig.4 Fitting curves of each component and resultant velocity

曲线拟合程度：径向方向、切向方向、垂向方向以及合速度方向的决定系数R2分别为0.97、0.97、0.98 和0.98，相关性较好，其中径方向相关性最好。 将拟合得到的K、α值代入萨道夫斯基公式，得到各分量及合速度的拟合公式如下：

式中：vx为水平径向的振动速度，cm／s；vy为水平切向的振动速度，cm／s；vz为垂直方向的振动速度，cm／s；vu为合速度，cm／s。

由于传感仪距离爆破点较近，且传感仪是插入土体中的，爆破后垂直方向上受到的冲击大，振动信号的采集出现了失真，导致拟合后的K、α过大。 但其他2 个方向上的拟合较为正常，可为冻土爆破振动传播的衰减规律提供一定的借鉴。 《爆破安全规程》中对爆破频率超过50 Hz 的巷道振动速度规定不能超过30 cm／s，将30 cm／s 作为振动速度控制值代入拟合公式，在一定药量下可以计算出安全距离R，根据R的取值，合理布置冻结管和立井井壁的浇筑位置。

3 冻土黏土振动信号的小波（包）分析

3.1 小波（包）的分解尺度

小波（包）分析是基于小波分析发展而来的，它能根据信号特性和分析要求自适应的选择相应频带与信号频谱相匹配，是比小波分解更为精细的一种分解方法［9］。 爆破振动信号进行小波（包）分解时，最优分解层数的确定依据信号分析的特征和振动记录仪的最小工作频率来确定［17］。 监测采用仪器为TC4850 爆破测振仪，其最小工作频率为5 Hz，由于爆破振动信号的频率一般集中在200 Hz 范围以内，根据采样定律［18］，设置信号的采样频率为16 kHz，其奈奎斯特频率为8 kHz。 根据小波（包）分析的原理将信号分解到8 层，共有28个小波（包），则其对应于的最低频段为0～31.25 Hz。

3.2 小波基的选择

小波基的选择在爆破振动信号分析中是至关重要的，不同的小波基会导致分析结果千差万别。Daubechies 小波系列具有良好的紧支撑性、光滑性以及对称性［19］。 根据算法进行matlab 编程，用db5～db10 的小波基对A测点的振动信号进行8 层小波（包）分解与重构，得出重构误差见表3。

表3 小波重构误差Table 3 Reconstruction errors of wavelet packet analysis

从表3 可以看出，db6 小波基对应的重构误差较小，因此选用db6 对振动信号进行小波（包）分解。

3.3 各频带能量的表征

依据小波分析原理将振动信号分解至第8 层，设S8，j所对应的能量为E′8，j，则有：

式中：xj，k（j＝0，1，2，…，2i－1；k＝1，2，…，m） 为重构信号S8，j的离散幅值。S8，j（t）为8 层第j个频带上的重构信号；k为离散点的个数；t为时间变量。

设被分析振动信号的总能量为E′，则有：

各频带能量的百分比为：

其中：j＝0，1，2，…，28－1。

通过计算爆破振动信号经小波分解后信号总能量和不同频带能量所占百分比，分析冻结黏土爆破振动信号在传播过程中的能量分布。 并通过小波（包）的分解，获得不同频带的能量及百分比。

3.4 冻结黏土爆破应力波能量的传播

3.4.1 冻结黏土爆破振动能量的衰减

根据试验得到的振动信号总能量，作出不同测点处的振动能量分布如图5 所示，从图5 中可以看出爆破作用下各点振动能量随测点距离的变化。

图5 爆破总能量随测点变化的曲线Fig.5 Total energy of blasting vibration signal changing with measuring point

由图5 可得，立井冻结黏土爆破振动信号与测点环境、爆破距离有关。 由于测点距离爆点仅3.0 ～7.43 m，炸药爆破后能量首先沿冻土传递到A测点，然后随着距离的增加，能量在冻土中逐渐耗散，传递到测点E时，能量已经大幅降低。 总体曲线呈现出由A测点向E测点逐渐衰减的过程。 此外，从图5 中3个方向上的总能量变化曲线中可以看出，距离爆点0～6.243 m 以内，3 个分量上的总能量呈现出E′z＞E′x＞E′y的规律，从另一个方面说明了立井冻结黏土爆破能量的传播是从内到外的径向传播。

3.4.2 不同频率的能量占比

为分析不同方向上的冻结黏土能量占比分布，了解冻结黏土能量的扩散规律，下面以A测点为例，选取0～500 Hz 频率内的能量，具体分析冻结黏土爆破3 个方向上的能量占比分布规律（图6）。

根据图6 可以看出，能量处于0 ～187 Hz 时，能量占比急剧下降，当能量在218 Hz 左右时，其占比又有所回升。 在218～250 Hz，冻土爆破振动能量发生了剧烈衰减，之后一直保持平稳的衰减状态。 可以得出：此次爆破振动信号能量的变化是在高频段衰减快，但是波动很大，低频段衰减慢，但是稳定。在高频段，能量先急剧下降，然后会频繁波动，会有所回升，达到峰值，之后能量一直衰减。 从能量占比来看，冻结黏土爆破振动能量主要集中在0 ～250 Hz。 在低频段切向分量能量占比较大，而在高频段径向分量能量占比较大。 张广辉等［20］对循环加载下冲击倾向性煤能量耗散进行了研究，从能量的角度发现：耗散能同样呈现先迅速降低，后缓慢增长的趋势。 依据式（8）将测点的能量占比进行求解，做出各点在不同频率的能量占比如图7 所示。

从图7 中可以得出，径向和切向的两个分量的能量占比主要集中在250 Hz 以内。 而各个测点的能量占比趋势也是符合先逐渐减小，然后在125 Hz左右增加，在218 Hz 处达到顶点，然后呈波动逐渐衰减，向低频段移动。 虽然个别测点存在能量波动较大的现象，与爆破现场不稳定有关。 值得关注的是轴向的能量变化趋势较于径向和切向更为明显，变化更大。 即下降比径向和切向更快，突增也比径向和切向更为明显，在高频段的频率占比也较大，易引起共振现象。

3.4.3 爆破振动能量的回归分析

爆破振动能量的衰减过程极其复杂，李洪涛等［21］将单段、多段爆破最大药量与总药量视为定值，推导出爆破振动能量衰减公式，该公式在应用研究中，取得了积极的成果。 但由于药量定值的限制，还存在一定误差。 单仁亮等［14］利用量纲分析法推导出的公式，较好弥补了这个缺点。 参考单仁亮等［14］研究，将岩石爆破振动参数转换为冻结黏土爆破参数进行研究。 该公式将能量E′与等效药量Q1、弹性模量E、密度ρ联系起来，函数关系如下：

式中：K1为与冻结黏土的密度、弹性模量有关的参数；φ为相关系数。

结合萨道夫斯基公式，将爆破振动能量衰减公式简化为：

根据式（8）编写拟合函数，利用Origin 软件求解的总能量值，进行数据的拟合。 以径向为例，如图8 所示，拟合回归系数见表4。

图8 冻土爆破能量衰减回归曲线Fig.8 Regression of blasting energy attenuation in frozen soil

表4 中径向、垂直方向上的拟合决定系数达到0.8 以上，切向的相关性较差，冻结黏土爆破引起冻土各方向上的能量随距离的增加呈衰减趋势。 垂向的K值较大，是因为垂向为炮孔自由面方向，爆破近区能量做功大，引起的振动也较大，由于能量快速释放，爆破远区急剧下降。

表4 能量衰减拟合系数Table 4 Fitting coefficient of energy attenuation

4 结 论

1）冻结黏土爆破各向振动速度存在一定差异，垂向振动速度＞径向振动速度＞切向振动速度，冻土爆破后应力波的传播主要是沿某个特定方向。

2）运用萨道夫斯基公式对单一冻结黏土介质爆破振动速度进行了拟合，获得了冻结黏土爆破的萨氏拟合公式。

3）深井冻结黏土爆破时，振动信号能量在高频段波动较大但衰减快，低频段相对平稳且衰减缓慢。在高频段，能量首先急剧下降，然后出现频繁波动，达到峰值后一直衰减。 垂向的能量变化趋势较于径向和切向更为明显，变化更大，下降速率径向和切向更快，突增也比径向和切向更为明显。

4）运用量纲分析法对冻结黏土爆破能量信号进行回归分析，得出了冻结黏土爆破各方向上的能量衰减规律。