薛 冰，凌 静，陈华东，张永明，赵 静，廖学燕，林谋金

(1.西南科技大学 环境与资源学院，绵阳 621010；2.泸州北方化学工业有限公司，泸州 646605；3.成都市城市安全与应急管理研究院，成都 610031)

爆破振动危害一直是爆破作业中研究和控制的重点，爆破振动达到一定强度会导致周围建构筑物和设施设备发生破坏[1,2]。现有的爆破振动控制研究主要关注减震沟降震[3,4]、预裂爆破[5,6]、装药结构和爆破网络优化等措施[7-10]。除此之外，通过采用新破岩技术降低振动效应也是有效的手段，比如采用静态膨胀破碎[11]、二氧化碳爆破以及金属燃烧剂弱爆破等技术取得了较好的减震效果[12,13]，但上述技术存在作用时间长、操作复杂或成本高昂等问题，限制其广泛应用。因此有必要寻找满足易操作、低成本、高效率工程应用需求的新型弱振动爆破器材。发射药作为一种常见的爆破器材，常被用于油气井致裂造缝。已有研究表明，其具有较好的破岩能力，做功形式多为燃烧或爆燃，且燃烧峰值压力小于传统工业炸药，在爆破振动控制中具有较好的应用前景[14,15]。当前，有关发射药爆破的研究主要关注其破岩和致裂性能，针对其爆破振动特性的研究较少。为拓展发射药爆破技术的应用范围，有必要进一步研究其爆破振动特性，为实际工程应用提供参考。

已有研究表明，爆破振动的峰值振动速度大小和主频范围是判断其破坏性的关键参数，因此在爆破器材的爆破振动特性研究中，主要关注其质点峰值速度及其随比例距离的变化规律、爆破振动的能谱分布等[16]。本研究通过系统开展不同药量的发射药爆破振动实验，采集不同爆心距的爆破振动信号，基于萨道夫斯基公式分析其振动速度衰减规律，利用希尔伯特-黄变换(Hilbert Huang Transform,HHT)方法获得其能谱分布规律，并与相同条件下的乳化炸药爆破振动信号进行对比研究，为发射药爆破振动预测与控制提供参考。

1 爆破振动实验

本次实验使用的3/1单基发射药由泸州北方化学工业有限公司提供，其采用电点火头引燃；乳化炸药由绵竹兴远特种化工有限公司提供，其采用8#电雷管起爆。为保证爆破材料燃烧和爆轰充分，使用内径30 mm、壁厚2 mm 的45#钢外壳加强约束，实验装药照片如图1所示。为获得不同比例距离条件，针对单基发射药和乳化炸药分别进行了3种药量的爆破振动测试，实验药量为20 g、30 g和50 g，并根据爆破材料类型进行相应编号：单基发射药(G20、G30、G50)，乳化炸药(E20、E30、E50)。

实验在混凝土地面进行，药包中心距混凝土自由面0.5 m，采用砂石填塞。爆破振动信号采集使用Blast-UM型测振仪，采样速率为10 000 Sps，采样时长为2 s，传感器X轴方向指向爆源，Y轴方向与爆源方向垂直，Z轴方向垂直于地面向上。测点爆心距为2～10 m，相邻测点间隔2 m，共计5个测点，布置如图2所示。

图 2 测点布置示意图(单位：m)Fig. 2 Arrangement of measuring points(unit:m)

2 爆破振动衰减规律及等效系数

2.1 爆破振动衰减规律

实验获得了不同药量时各测点水平径向、水平切向、垂直方向上的质点峰值振动速度以及峰值振动速度合速度，整理如表1所示。

表 1 不同测点爆破振动测试结果

由表1可知，两种爆破器材的垂直方向峰值速度均大于水平方向的峰值振动速度，且三个方向的峰值振动速度随爆心距的增加而减小。通过对比两种爆破器材振动数据，发现在三个速度方向上单基发射药峰值振动速度均小于同药量的乳化炸药，单基发射药各方向峰值振动速度为乳化炸药的36.05%～71.62%。

相关文献研究表明[17]，利用爆破振动合速度进行爆破振动衰减规律研究具有一定代表性。萨道夫斯基公式是研究爆破振动衰减规律的常用方法，基于萨道夫斯基公式(式1)，利用最小二乘法对表1中单基发射药(vG)和乳化炸药(vE)峰值振动合速度进行回归分析，得到相应的拟合结果(式2、3)，如图3所示。

(1)

(2)

(3)

式中:v为测点振动速度，cm·s-1；Q为药量，kg；R为测点到爆源的距离，m；K、α为与质地条件有关的相关系数和衰减系数。

对比两种爆破器材拟合结果，发现单基发射药和乳化炸药爆破振动衰减曲线形状相似，在衰减公式中表现为两者的α值接近。两者爆破振动衰减规律的差异主要体现于K值大小不同。该差异是由于两种爆破器材反应类型和对外做功特性不同而导致。乳化炸药主要反应形式为爆轰，其在起爆后迅速对周围介质产生高压破坏，形成强烈的冲击波，引起剧烈的振动效应。而单基发射药的主要反应形式为燃烧或爆燃，其被引燃后发生快速燃烧导致压力上升，从而进一步加速燃烧，直至达到峰值压力。其压力形成速度和峰值均小于乳化炸药，但压力持续时间长于乳化炸药，使得相同药量情况下其振动速度小于乳化炸药。

图 3 单基发射药与乳化炸药振动衰减公式拟合曲线图Fig. 3 Fitting curve of vibration attenuation formula of single-base gun propellant and emulsion explosive

2.2 爆破振动等效系数

(4)

等效系数能直观反映出相同振动速度时单基发射药与乳化炸药的药量关系。对所有比例距离时的等效系数进行统计分析，如图4所示。

图 4 单基发射药振动的乳化炸药等效系数计算结果Fig. 4 Calculation of emulsion-explosive equivalent ratio of single-base gun propellant

从图4可以看出，在本次实验条件下等效系数分布较为集中，其分布范围为0.31～0.52，平均值为0.43。而在相同振动速度时，等效乳化炸药药量与实际单基发射药药量比值为等效系数的立方η3，说明单位质量单基发射药与乳化炸药爆破振动效应差异巨大。单基发射药爆破振动的乳化炸药等效系数可为其实际应用提供参考。

3 爆破振动主频及能量分布

已有研究表明，除峰值振动速度外振动频率也是影响爆破振动破坏效应的主要因素，这是因为当爆破振动频率与建构筑物自振频率接近时易出现共振导致破坏加剧，而建构筑物的自振频率多为几赫兹到二三十赫兹的低频，所以开展爆破振动研究需要关注其频率分布[18]。为了更好地研究单基发射药和乳化炸药爆破振动特性，利用希尔伯特-黄变换对比分析了两种爆破器材的能频谱。通过数据预处理发现实验中乳化炸药和单基发射药的振动信号频率均主要集中于0～50 Hz以内。因此，本研究主要讨论单基发射药和乳化炸药爆破振动在0～10 Hz、10～50 Hz两个频率区间的能量分布。

3.1 爆破振动能谱分布

希尔伯特-黄变换(式5)是一种非线性非平稳信号的数据处理方法[19]。该方法主要由经验模态分解(Empirical Mode Decomposition,EMD)和Hilbert谱分析(Hilbert Spectrum Analyze,HSA)两部分组成。其中，Hilbert谱具有较好的局部化能力，而Hilbert能量谱显示了能量随时频变化的具体分布。为分析单基发射药与乳化炸药爆破振动能量随时频变化的规律，采用EMD算法进行计算，得到原始信号的多个本征模态函数(IMF)，对这些IMF进行Hilbert变换，得到每个IMF分量的瞬时频率，即Hilbert谱。

(5)

对IMF分量的瞬时频率进行进一步处理(式6)，即可得到对应边际能量谱

(6)

上式表达了每个频率在全局上的能量，为统计意义上的全部累加能量，即描述能量随频率的分布情况。根据爆破振动实验垂直方向振动信号HHT变换结果，绘制出边际能量谱(图5)如下。

由图5可知，单基发射药振动信号能量主要集中于0～5 Hz和15～25 Hz频率范围，且随着爆心距增加0～5 Hz频段内的能量逐渐减小，而主要集中于15～25 Hz频率范围。当爆心距为6 m、8 m、10 m处的高频振动能量趋近于0，这符合地震波在粘弹性介质中的传播规律，即高频波被不均匀且具有粘弹性特性的混凝土介质吸收[20]。而乳化炸药的能量分布则较分散，在0～50 Hz范围内呈现多峰分布。对比相同药量时单基发射药和乳化炸药的能谱可知，在20 g药量时两者能谱分布及幅值接近。而随着药量的增加，两种爆破器材的能谱分布差别增大，且各频段的乳化炸药能量幅值均要高于单基发射药。

3.2 爆破振动能量占比分析

为了更直观分析两种爆破器材能量累积和分布特性并方便工程应用，参考爆破安全规程分别讨论两种爆破器材在0～10 Hz与10～50 Hz两个频段内的能量占比，整理得到两频段能量相对总能量占比图(图6)。

由图6可知，爆心距和实验药量对爆破振动总能量和不同频段能量占比有显著影响。单基发射药与乳化炸药振动信号总能量与10～50 Hz频段能量均随着爆心距的增加而迅速减小，且单基发射药爆破振动信号总能量水平远小于乳化炸药，这与振动速度分析结果一致。但随着实验药量的增加，单基发射药爆破振动总能量与乳化炸药的差距在减小。此外，随着乳化炸药药量的增加，其0～10 Hz频段的能量占比在逐渐减小。而单基发射药0～10 Hz频段能量随爆心距的变化则存在一定的波动，有时甚至超过乳化炸药，且与药量的变化存在耦合作用。因此，在实际应用中应关注单基发射药火药爆破振动在0～10 Hz频段的能量和速度。

图 5 实验振动信号的HHT边际能量谱Fig. 5 HHT marginal energy spectrum of experimental vibration signals

图 6 实验振动信号的能量占比图Fig.6 Energy ratio diagram of experimental vibration signals

4 结论

2)通过HHT获得了爆破振动信号的边际能量谱，发现单基发射药振动信号能量水平要低于乳化炸药，且单基发射药振动信号主要集中于0～5 Hz和15～25 Hz频率范围，而乳化炸药在0～50 Hz范围内能量分布较分散；且随着爆心距增加，单基发射药0～5 Hz频段内的能量逐渐减小，而主要集中于15～25 Hz频率范围。

3)参考爆破安全规程发现，两种爆破器材振动总能量及10～50 Hz频段能量均随爆心距的增加而减小；但随着药量的增加，单基发射药与乳化炸药振动信号总能量的比值在减小，且在实际应用中需关注单基发射药在0～10 Hz频段能量波动的影响。