基于正态分布的爆破振动评价与安全药量计算

张小军; 汪旭光; 王尹军; 于亚伦; 吴春平; 杨德强

doi:10.11883/bzycj-2017-0091

基于正态分布的爆破振动评价与安全药量计算

doi: 10.11883/bzycj-2017-0091

1.
北京科技大学土木与资源工程学院, 北京 100083
2.
北京矿冶研究总院, 北京 100160

基金项目:

国家自然科学基金项目 50704005

中国工程院咨询研究项目 2016-XZ-09

详细信息

作者简介:
张小军(1991-), 男, 博士研究生, 1024770807@qq.com

中图分类号: O383.1
计量
- 文章访问数: 5355
- HTML全文浏览量: 2001
- PDF下载量: 114
- 被引次数: 26
出版历程
- 收稿日期: 2017-03-24
- 修回日期: 2017-04-17
- 刊出日期: 2018-09-25

Blasting vibration evaluation and safety dose calculation based on normal distribution

1.
School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China
2.
Beijing General Research Institute of Ming and Metallurgy, Beijing 100160, China

摘要

摘要: 工程爆破是非常重要的一种施工方法，但同时爆破引起的岩石振动也是爆破公害之一。由于爆破场地的复杂性，爆破地震引起的某位置质点振速峰值v与单段最大药量Q、爆心距R没有严格的函数关系，只能将振速视为随机变量，因此振动只能从概率的角度来描述。为了控制振动强度，达到较高的可靠度，必须计算振速小于目标设施安全振速的概率。本文中基于概率，将振速峰值近似视为服从正态分布，对目标设施进行安全分析以及安全炸药量计算，最后通过案例应用，解释概率公式计算炸药量的合理性。
- 爆破 /
- 振动 /
- 振速峰值 /
- 正态分布 /
- 安全药量
Abstract: Blasting is an indispensable construction method to human beings. However, in the meanwhile, the rock vibration caused by blasting is one of blasting hazards. Due to the complicated geology, the peak particle velocity is a random variable approximately. It does not have a general function relationship with the maximum dose Q and the distance R from blasting center. Therefore, the vibration can only be described by the theory of probability. In order to control the vibration intensity and get high reliability, it is necessary to calculate the probability that the peak particle velocity is less than the safety vibration velocity of the target facility. Based on the probability, the distribution of peak particle velocity is approximately regarded as the normal distribution. Then safety analysis of target facilities and calculation of safe explosives are carried out. Finally, through the application of the case, it is reasonable to calculate the dose by probability formula.
- blasting /
- vibration /
- peak particle velocity /
- normal distribution /
- safety dose

HTML全文

炸药在密闭空间内爆炸，冲击波在空间内往返运动，并在该过程中逐渐减弱并消失，之后空间内的压力分布较均匀，随时间变化缓慢，这种压力状态被称为准静态压力状态。准静态气压状态下空间内的压力因气体与外界热交换而缓慢变化，而一般受关注的是最初的准静态压力值或者准静态压力峰值，下文将准静态压力峰值简称为准静态压力。

准静态压力受到的关注远不及爆炸冲击波形成的脉冲压力，但其仍是非常重要的力学参数，比如：准静态压力是爆炸容器安全设计的重要参考数据^[1-2]；另外，对于地下核爆炸，爆室内的准静态压力是核爆炸形成的剧毒物质在岩土孔隙和裂缝中扩散的动力^[3]。国外实验研究工作较丰富^[4-8]，Anderson等^[7]基于文献中的实验数据，利用相似理论拟合得到容器内爆炸准静态气压无量纲峰值的经验公式，Anthony等^[8]通过实验得到了炸药LX-13在容器内爆炸形成的准静态压力经验公式，Bultman^[1]在容器设计时也引用了基于实验数据的经验公式，但这些经验公式差异较大。国内的研究工作不多，特别是实验工作比较缺乏。王等旺等^[9]、林俊德^[10]开展了柱形容器内准静态压力的实验研究，钟方平^[11]通过数值计算研究了柱形容器内的准静态压力，钟巍等^[12]研究了考虑化学反应动力学影响的约束爆炸准静态压力的计算方法。

本文中，在球形爆炸容器内开展了4种不同当量的爆炸加载实验，实验中压力传感器采用了2种安装方式，均获得了理想数据；还理论推导了准静态压力的表达式，并通过拟合实测数据得到具体的经验公式。

1. 实验

实验在内径523 mm的球形爆炸容器内进行，容器如图 1所示。炸药采用球形装药，装药包括主装药和起爆微型药球两部分，如图 2所示。主装药由2个半球组成，为60%RDX和40%TNT浇铸而成。起爆药球的当量为1 g TNT，直径为10 mm，由微米级的PETN粉压制而成。起爆药球连接一根柔爆索，其中心为柔爆索的一端。柔爆索直径为1 mm，线装药密度为0.5 g RDX/m，为铅皮柔爆索。

图 1 实验用球形爆炸容器

Figure 1. The spherical explosion containment vessel used in experiment

下载: 全尺寸图片幻灯片

图 2 实验装药

Figure 2. Explosive installation

下载: 全尺寸图片幻灯片

实验时采用细绳吊挂装药，装药置于容器中心。容器内除了炸药、导爆索和细绳，再无其他物体，这样的做法可以尽量减少容器内除炸药以外的多余物体。雷管处于雷管罩内，与容器内部体积隔开。由雷管引爆较细的柔爆索，柔爆索在中心起爆微型药球，微型药球再起爆球形主装药，该起爆方式能保证形成较理想的球面冲击波。采用27、60、185、370 g TNT 4种当量，对应的比距离范围为0.36~0.87 m/kg^1/3。

容器赤道面上均布4个用于传感器安装的接口。压力传感器安装包括齐平安装和导孔安装2种方式，如图 3所示。齐平安装方式中传感器敏感面与容器内壁齐平；导孔安装方式设计了一个空腔，传感器敏感面处于空腔内，该空腔通过6个直径1 mm的导孔连接容器。导孔安装方式可以避免传感器遭受强冲击波或高速破片的破坏，但也导致整个测试的频响降低^[13]，因此该方式仅适合准静态压力这种变化较缓慢的参数测试。

图 3 压力传感器安装结构示意图

Figure 3. Diagrams of pressure sensor installations

下载: 全尺寸图片幻灯片

实验中采用了某型压电式压力传感器和某型压阻式压力传感器，参数见表 1。由于实验中压阻式传感器在采用齐平安装方式下容易损坏或实测信号质量不高，最终压阻式传感器采用导孔安装方式，压电式传感器采用齐平安装方式。2种压力测量系统通过同步机实现时间同步。

表 1 压力传感器参数

Table 1. Performance parameters of pressure sensors

类型	量程/MPa	频响/kHz	上升时间/μs	非线性度	精度
压电式传感器	34.48, 68.95	500	≤1.0	≤1.0%Fs	0.14, 0.28 kPa
压阻式传感器	15, 40	600	≤0.3	≤1.0%Fs	0.5%Fs

下载: 导出CSV

| 显示表格

2. 实验结果及分析

2.1 实验结果

图 4为齐平安装压力传感器测得的典型波形。可见压力波形的前期(爆炸后1 ms之前)存在3次较明显的脉冲，表明炸药在容器中心爆炸后，形成的冲击波在容器壁和容器中心之间往返运动，在容器壁上反射了3次；在3个脉冲之后，压力波形变化较缓，表明容器内的压力趋于平静，进入准静态压力状态。当然，在比距离比本文的更小的工况下，也可能出现冲击波往返次数多于3次，容器内压力进入准静态状态的现象。准静态压力的取值是对压力波形的准静态压力部分进行平均处理，图 4中冲击波反射压力达到约16 MPa，准静态压力仅为1.5 MPa。

图 4 齐平安装传感器获得典型压力波形(当量60 gTNT)

Figure 4. Pressure-time curve obtained by the pressure sensor of flush installation (when explosion equivalent is 60 gTNT)

下载: 全尺寸图片幻灯片

图 5为导孔安装压力传感器获得的典型压力波形。可见波形早期有一个压力最大值(约3.5 MPa)，之后压力缓慢变化，且压力在初期时间内(爆炸后0.2 s之前)近似地随时间线性下降。

图 5 导孔安装传感器获得的典型波形(当量60 gTNT)

Figure 5. Pressure-time curve obtained by the pressure sensor of guide-hole installation (when explosion equivalent is 60 gTNT)

下载: 全尺寸图片幻灯片

虽然导孔安装传感器未与容器内部空间直接接触，而是通过6个直径1 mm的小孔连接，但爆炸冲击波通过小孔仍会对传感器加载一个较弱的冲击载荷。因此，简单地读取图 5中压力波形的最大值作为准静态压力值是不合适的，该峰值应该是冲击波压力在低频测试系统中的表现。准静态气压是由冲击波压力连续变化而来的，如何确定准静态气压值显得比较困难，这可能也是不同文献中实测准静态气压差别较大的原因之一。本文的取值方法如图 6所示，按照线性下降斜率进行反向回推，与波形上升沿交于一个点，该点的压力值即为准静态气压。这种取值方法也是Anderson等^[7]推荐的方法。

图 6 准静态压力值的选取方法

Figure 6. The method for determining quasi-static pressure

下载: 全尺寸图片幻灯片

表 2为两种安装方式的传感器测量到的准静态压力值，可见两种测量方式得到的准静态压力值相近。

表 2 压力传感器参数

Table 2. Quasi-static pressures obtained by pressure sensors

当量/(gTNT)	准静态压力/MPa
当量/(gTNT)	压电式传感器	压阻式传感器
27	0.7, 0.7, 0.6
60	1.4, 1.5	1.5, 1.2
185	2.4	3.7
370	5.5, 4.9, 5.3	5.4

下载: 导出CSV

| 显示表格

值得注意的是，由于冲击波的作用，压力波形的早期存在一个远高于准静态压力值的压力峰，如图 4~5所示。传感器的量程必须根据该压力峰值来进行选取，但传感器量程过高必然导致准静态压力值的测量精度下降。

2.2 数据分析

密闭空间内的准静态气压主要是爆炸气体产物在爆炸释放热量的作用下温度升高而引起的气体热运动压力。影响准静态气压的主要因素包括两部分：一是由爆炸形成的大量气体产物使得密闭空间内气体量增大；二是爆炸释放的热量使气体的温度升高。下面基于理想气体的假设，简单地推导密闭空间内准静态压力的表达式。

理想气体的状态方程为:

$p = \rho RT$

(1)

式中：ρ为气体密度；R为理想气体参数，为普适气体参数R₀除以气体摩尔质量M，即R=R₀/M；T为绝对温度。

忽略气体与外界热交换带来的能量损失，且忽略空间内介质气体的质量(因为其质量远小于炸药的质量)，同时假设爆炸释放的能量全部用来加热球形容器内的气体物质(因为大部分爆炸能量转换为了气体内能)，则气体的温升为:

$\Delta T = \frac{E}{{Q{c_V}}}$

(2)

式中：E为爆炸释放的总能量，Q为爆炸当量，c_V为容器内气体的平均定容比热。

实验初始状态为常温、常压，因此爆炸后密闭空间内的准静态压力为:

$\begin{array}{l} {p_{\rm{q}}} = \rho RT = \frac{Q}{V}\frac{{{R_0}}}{M}({T_0} + \Delta T) = \\ \frac{Q}{V}\frac{{{R_0}}}{M}\left( {{T_0} + \frac{E}{{Q{c_V}}}} \right) = \frac{Q}{V}\left( {\frac{{{R_0}{T_0}}}{M} + \frac{{e{R_0}}}{{M{c_V}}}} \right) \end{array}$

(3)

式中：T₀为初始温度，V为容器体积，e为炸药释放的质量能量，即比内能。对于特定的炸药来说，气体摩尔质量M、定容比热c_V、比内能e皆为常数，从式(3)可以看出，空间内准静态压力与当量容积比Q/V呈正比例关系。该关系式与Pastrnak等^[5]采用的表达式形式相同。

根据表 2中数据，绘制准静态压力与当量容积比Q/V的变化关系，如图 7所示。可见，2种安装方式传感器获得的数据变化规律一致。另外，也可以看出，球形爆炸容器内准静态压力与当量容积比Q/V呈近似的线性关系。

图 7 球形容器内准静态压力与当量容积比的关系

Figure 7. Relation between quasi-static pressure and explosion equivalent-to-vessel volume ratio in spherical vessel

下载: 全尺寸图片幻灯片

基于式(3)和表 2中的实测数据，可以拟合得出准静态压力经验公式：

${p_{\rm{q}}} = 1.10Q/V\;\;\;\;\;\;0.36\;{\rm{kg}}\;{\rm{TNT}}/{{\rm{m}}^3} \le Q/V \le 4.94\;{\rm{kg}}\;{\rm{TNT}}/{{\rm{m}}^3}$

(4)

式中p_q单位为MPa，当量单位为kg TNT，容器体积单位为m³。

在理论推导准静态压力表达式—式(3)的过程中，忽略了容器内介质气体的质量，因为大部分情况下介质气体的质量远小于炸药质量，但当炸药质量较小，与容器内介质气体质量相当时，这种处理方式是不合适的，此时经验公式(4)不再适用。另外，文献[10]表明，容器内的氧气环境对准静态压力有影响，在爆炸后的高温环境中，爆轰产物与氧气化合燃烧，增大容器内气压，如果爆炸当量比本文实验的爆炸当量更大，可能导致氧气相对不足，爆轰产物与氧气的化合燃烧不充分，此时经验公式(4)也不再适用。因此，基于现有的实验结果，当量容积比范围0.36~4.94 kg TNT/m³是式(4)成立的重要前提条件，如果当量容积比不在该范围，式(4)可能不再适用。

再者，式(3)显示，炸药种类不变的情况下准静态压力和当量容积比之间关系的线性系数才是一个常数。因此，式(4)可能不适用于性质与本文实验所用炸药相差较大的炸药。

3. 结束语

在球形容器内开展了4种当量的爆炸加载实验，实验中压力传感器采用了齐平和导孔两种安装方式，均获得了准静态压力数据，且两组数据一致。研究结果表明：

(1) 在比距离0.36~0.87 m/kg^1/3、当量容积比0.36~4.94 kg TNT/m³范围内，爆炸冲击波在球形容器内往返3次后，容器内气体进入准静态压力状态。

(2) 球形容器内准静态压力与当量容积比呈正比例关系，系数为1.10 MPa·m³/ kg TNT。

图 1 数据拟合曲线

Figure 1. Data fitting curve

下载: 全尺寸图片幻灯片

表 1 锡铁山铅锌矿爆破振动监测数据

Table 1. Blasting vibration experiment data

测点	Q/kg	R/m	(Q^1/3·R^-1)/(kg^1/3·m^-1)	v/(cm·s^-1)
1	30.2	20	0.156	11.82
2	30.2	40	0.078	3.47
3	30.2	60	0.052	1.77
4	30.2	80	0.039	1.28
5	30.2	100	0.031	0.84
6	74.8	20	0.211	17.24
7	74.8	40	0.105	5.62
8	74.8	60	0.070	2.72
9	74.8	80	0.053	1.84
10	74.8	100	0.042	1.22

下载: 导出CSV

表 2 数据处理结果

Table 2. Data processing results

测点	v_exp/(cm·s^-1)	v_cal/(cm·s^-1)	(v_exp-v_cal)/(cm·s^-1)
1	11.82	10.92	0.90
2	3.47	3.63	-0.16
3	1.77	1.90	-0.13
4	1.28	1.21	0.07
5	0.84	0.85	-0.01
6	17.24	17.67	-0.43
7	5.62	5.87	-0.25
8	2.72	3.08	-0.36
9	1.84	1.95	-0.11
10	1.24	1.37	-0.13

下载: 导出CSV

表 3 单段最大药量和设施安全概率

Table 3. Single biggest dosage andfacilities security probability

方法	Q/kg	P/%
萨道夫斯基公式	858.73	50.00
正态分布函数	797.89	95.05

下载: 导出CSV

参考文献(11)

[1]	饶运章, 汪弘.爆破振动速度衰减规律的多元线性回归分析[J].金属矿山, 2013(12):46-51. http://d.old.wanfangdata.com.cn/Periodical/jsks201312014 RAO Yunzhang, WANG Hong. Multiple regression linear analysis on attenuation formula of blasting vibration velocity[J]. Metal Mine, 2013(12):46-51. http://d.old.wanfangdata.com.cn/Periodical/jsks201312014
[2]	胡建华, 尚俊龙, 罗先伟, 等.单孔爆破振动监测与衰减规律多元线性化回归[J].振动与冲击, 2013, 32(16):49-53. doi: 10.3969/j.issn.1000-3835.2013.16.009 HU Jianhua, SHANG Junlong, LUO Xianwei, et al. Monitoring of single-hole blasting vibration and detection of its attenuation law by using multiple linear regression[J]. Journal of Vibration and Shock, 2013, 32(16):49-53. doi: 10.3969/j.issn.1000-3835.2013.16.009
[3]	卢文波, HUSTRULID W.质点峰值振动速度衰减公式的改进[J].工程爆破, 2002, 8(3):1-4. doi: 10.3969/j.issn.1006-7051.2002.03.001 LU Wenbo, HUSTRULID W. An improvement to the equation for the attenuation of the peak particle velocity[J]. Engineering Blasting, 2002, 8(3):1-4. doi: 10.3969/j.issn.1006-7051.2002.03.001
[4]	言志信, 严浬, 江平, 等.爆破振动峰值预报方法探讨[J].振动与冲击, 2010, 29(5):179-182. doi: 10.3969/j.issn.1000-3835.2010.05.039 YAN Zhixin, YAN Li, JIANG Ping, et al. Prediction methods for blasting induced ground vibration velocity[J]. Journal of Vibration and Shock, 2010, 29(5):179-182. doi: 10.3969/j.issn.1000-3835.2010.05.039
[5]	高振儒, 雷俊丽, 范磊, 等.爆破振动强度预报的置信度和安全系数[J].工程爆破, 2005, 11(3):69-71. doi: 10.3969/j.issn.1006-7051.2005.03.021 GAO Zhenru, LEI Junli, FAN Lei, et al. Faith degree and safety coefficient of predicting blasting vibration intensity[J]. Engineering Blasting, 2005, 11(3):69-71. doi: 10.3969/j.issn.1006-7051.2005.03.021
[6]	吕涛, 石永强, 黄诚, 等.非线性回归法求解爆破振动速度衰减公式参数[J].岩土力学, 2007, 28(9):1871-1878. doi: 10.3969/j.issn.1000-7598.2007.09.019 LV Tao, SHI Yongqiang, HUANG Cheng, et al. Study on attenuation parameters of blasting vibration by nonlinear regression analysis[J]. Rock and Soil Mechanics, 2007, 28(9):1871-1878. doi: 10.3969/j.issn.1000-7598.2007.09.019
[7]	贾晓强, 方向, 张卫平, 等.基于置信度理论的岩石爆破振动安全控制研究[J].爆破器材, 2011, 40(3):31-34. doi: 10.3969/j.issn.1001-8352.2011.03.009 JIA Xiaoqiang, FANG Xiang, ZHANG Weiping, et al. Analysis on security control of rock blasting vibration based on faith degree theory[J]. Explosive Materials, 2011, 40(3):31-34. doi: 10.3969/j.issn.1001-8352.2011.03.009
[8]	李庆扬, 王能超, 易大义.数值分析[M].北京:清华大学出版社, 2001.
[9]	王国政, 刘洋, 肖继红, 等.概率论与数理统计[M].重庆:重庆大学出版社, 2015.
[10]	NATEGHI R. Prediction of ground vibration level induced by blasting at different rock units[J]. International Journal of Rock Mechanics and Mining Science, 2011, 48(6):899-908. doi: 10.1016/j.ijrmms.2011.04.014
[11]	陈钧璠.爆破振动公式中K值上下限的探讨[J].爆破, 1996, 13(2):13-15. doi: 10.3969/j.issn.1001-8352.2011.03.009 CHEN Junfan. Determination of lower and upper limits of K in equation of blasting vibration[J]. Blasting, 1996, 13(2):13-15. doi: 10.3969/j.issn.1001-8352.2011.03.009

施引文献

期刊类型引用(16)

1.	陈青华，陶彦光，梁振刚，杨佳岐. 双舱室模型内爆炸临舱压力载荷特性分析. 舰船科学技术. 2024(11): 75-79 . 百度学术
2.	李营，杜志鹏，陈赶超，王诗平，侯海量，李晓彬，张攀，张伦平，孔祥韶，李海涛，郭君，姚术健，王志凯，殷彩玉. 舰艇爆炸毁伤与防护若干关键问题研究进展. 中国舰船研究. 2024(03): 3-60 . 百度学术
3.	胡俊华，董奇，胡八一，任逸飞，黄广炎. 抗爆容器的内部爆炸效应和动态力学行为研究进展. 含能材料. 2024(09): 986-1008 . 百度学术
4.	王昭，吴祖堂，杨军，李焰，刘文祥. 新型薄膜式压力传感器的参数设计. 爆炸与冲击. 2023(07): 174-184 . 本站查看
5.	冯蕴雯，林心怡，薛小锋，杨祥，刘佳奇. 高可靠单向爆破的民机防爆结构设计. 航空学报. 2023(18): 191-205 . 百度学术
6.	刘文祥，张德志，程帅，马艳军. 爆炸容器研究进展. 现代应用物理. 2023(03): 61-69 . 百度学术
7.	刘欣，顾文彬，蔡星会，王涛，刘建青，王振雄，沈慧铭. 圆柱形爆炸容器的内壁爆炸载荷. 爆炸与冲击. 2022(02): 19-30 . 本站查看
8.	刘正，聂建新，徐星，朱英中，刘攀，郭学永，闫石，张韬. 密闭空间内六硝基六氮杂异伍兹烷基复合炸药能量释放特性. 兵工学报. 2022(03): 503-512 . 百度学术
9.	张鹏宙，董奇，杨沙. 爆炸载荷特征参数对无限长圆柱壳弹性动态响应的影响. 爆炸与冲击. 2021(06): 48-57 . 本站查看
10.	孙琦，董奇，杨沙，张刘成. 内爆炸准静态压力对球形容器弹塑性动态响应的影响. 含能材料. 2020(01): 25-31 . 百度学术
11.	张翔，崔春生，刘双峰，刘嘉颖，崔光强. 子弹在爆炸容器内的爆燃压力测试技术研究. 中国测试. 2020(01): 39-43 . 百度学术
12.	张龙，邹虹，张宝国，张继军，张东亮，孔德骞. 有限空间爆炸静态压力的温度补偿方法. 爆炸与冲击. 2020(03): 100-109 . 本站查看
13.	徐景林，顾文彬，刘建青，王振雄，陆鸣，徐博奥. 圆柱形爆炸容器内爆炸载荷的分布规律. 振动与冲击. 2020(18): 276-282 . 百度学术
14.	夏彬伟，高玉刚，刘承伟，欧昌楠，彭子烨，刘浪. 缝槽水压爆破破岩载荷实验研究. 工程科学学报. 2020(09): 1130-1138 . 百度学术
15.	张明明，张连生，王鑫. TNT内爆准静态压力实验和数值模拟研究. 兵器装备工程学报. 2019(05): 195-199 . 百度学术
16.	孙琦，董奇，杨沙，张刘成. 内爆炸准静态压力对球形容器弹性动态响应的影响. 含能材料. 2019(08): 698-707 . 百度学术