2024 年 4 期目次

2024 年 4 期目次

计量
- 文章访问数: 77
- HTML全文浏览量: 85
- PDF下载量: 42
- 被引次数: 22
出版历程
- 刊出日期: 2024-04-07

20世纪以来，工程爆破已深入到国民经济建设的各个领域，工程爆破是完成人力和机械力所不能胜任的一种非同寻常的施工方法。但是爆破引起的振动是最突出的爆破公害之一，所以进行爆破振动的安全评价是实现对它的准确预报、有效控制和安全实施爆破的迫切需要。爆破引起的质点振动速度峰值，常用萨道夫斯基公式进行计算，但是爆破引起的岩土质点振速实际峰值v往往随机性比较大，它与Q^1/3/R的关系并不是一般意义上的“函数关系”，v与Q^1/3/R的关系是“随机变量”的“相关关系”。其原因是，岩土是经过漫长地质年代形成的地质体，其内部包含大量的裂隙，这些裂隙相当复杂，造成岩石内部的不连续和不均质性，从爆破测振用回归分析计算K、α，也说明了这一点。

饶运章等^[1]利用SPSS(statistical product and service solutions)软件进行了爆破振速衰减规律计算；胡建华等^[2]利用多元线性回归方法，分析了单孔爆破条件下的振动衰减规律；卢文波等^[3]基于柱面波理论、长柱状装药中的子波理论以及短柱药包的应力波场Heelan解的分析，推导了岩石爆破中质点峰值振动速度的衰减公式，对现有的公式进行了改进；言志信等^[4]探讨了爆破振动峰值速度预报的公式法和Fourmap法，特别地还尝试了利用人工神经网络预报了爆破振动峰值速度；高振儒等^[5]对测试数据用不同的置信度进行拟合，确定系数K、α，用安全系数对回归式预报爆破振动强度的安全性进行评估；吕涛等^[6]通过线性回归法和非线性回归法，得到萨道夫斯基公式和其修正公式，研究爆破振动衰减规律。上述研究基本解决了各类爆破振动衰减公式的求解问题，即对萨道夫斯基公式中K、α的确定，但对质点振速实际峰值v与Q^1/3/R的关系是“随机变量”的“相关关系”的特性并没有研究。

本文中，基于概率论中的正态分布对“随机变量”爆破振速峰值v和Q^1/3/R的“相关关系”特性进行分析，利用随机变量v的分布函数提出爆破振动安全评价的方法，给出安全炸药量的计算方法，并通过案例进行验证^[7-8]。

1. 随机变量v的分布函数

由于萨道夫斯基公式所表达的是v与Q^1/3/R的函数关系，即对于给定的Q^1/3/R，通过该公式计算的v₀只是质点振速实际峰值v的期望值或估计值，实际值v_a落在v₀的附近，具体在v₀的多远处与概率有关，可以看作近似服从正态分布。通过分析可以看出，只计算出K、α，求出的v，可靠度是不高的，要获得更高的可靠度，必须得出随机变量v的分布函数，当然求解过程仍然离不开萨道夫斯基公式。

1.1 线性回归法确定K、α

为了确定萨道夫斯基公式中的K、α，需要根据最小二乘法原理^[9]，对实测数据使用线性回归法进行拟合。萨道夫斯基公式为：

$v = K{\left( {\frac{{{Q^{1/3}}}}{R}} \right)^\alpha } = k{\rho ^\alpha }$

(1)

式中:v为质点的最大振速(cm/s)；Q为单段最大炸药量(kg)；R为测点至爆源的距离(m)；K、α为与爆破点地形、地质等条件相关的系数和衰减指数；ρ为比例药量。

从式(1)可以看出，v与Q^1/3/R不是线性关系，为了便于回归分析，需处理成线性关系。对式(1)的两边取对数，得到：

$\ln v = \ln K + \alpha \ln \left( {\frac{{{Q^{1/3}}}}{R}} \right)$

(2)

设：

$\ln v = y, \;\;\;\;\;\;\;\ln k = {a_0}, \;\;\;\;\;\;\;\ln \left( {\frac{{{Q^{1/3}}}}{R}} \right) = x, \;\;\;\;\;\;\alpha = {a_1}$

(3)

则式(2)变成：

$y = {a_0} + {a_1}x$

(4)

式(4)是线性关系，根据最小二乘法原理以及实测数据，求常数a₀和a₁。那么其方程组为：

$\left\{ \begin{gathered} {a_0}\sum\limits_{i = 1}^m {{w_i}} + {a_1}\sum\limits_{i = 1}^m {{w_i}{x_i}} = \sum\limits_{i = 1}^m {{w_i}{y_i}} \hfill \\ {a_0}\sum\limits_{i = 1}^m {{w_i}{x_i}} + {a_1}\sum\limits_{i = 1}^m {{w_i}x_i^2} = \sum\limits_{i = 1}^m {{w_i}{x_i}{y_i}} \hfill \\ \end{gathered} \right.$

(5)

式中：m为现场爆破实验次数，w_i=1。通过方程组得到a₀、a₁，再根据式(3)的变量代换关系，可得：

$K = {{\text{e}}^{{a_0}}}, \;\;\;\;\;\;\;\;\;\;\alpha = {a_1}$

(6)

1.2 正态分布函数

随机变量ξ与任一实数x的关系式“ξ≤x”表示一个事件，即为{w|ξ(w)≤x}。其概率P(ξ≤x)与实数x有关，应为实数变量x的函数，称为随机变量的分布函数^[10]，记作F_ξ(x)：

${F_\xi }\left( x \right) = P\left( {\xi \leqslant x} \right)$

(7)

根据v与Q^1/3/R关系分析可以看出，对于任意给定的Q^1/3/R，实际值v_a会落在萨道夫斯基公式计算的v₀附近，偏离v₀越远，可能性越小，越靠近v₀，可能性越大。根据这个性质，可以近似视随机变量v服从正态分布N(μ，δ²)，μ指按最小二乘法曲线拟合公式的计算值，δ²指实际爆破振速峰值偏离拟合公式计算值的程度：

$\mu = v = K{\left( {\frac{{{Q^{1/3}}}}{R}} \right)^\alpha }$

(8)

${\sigma ^2} = \frac{{\sum\limits_{i = 1}^m {{{\left( {{v_{{\text{mea, }}i}}-{v_{{\text{cal}}.i}}} \right)}^2}} }}{m}$

(9)

式中：v_mea为实际爆破振速峰值测量值；v_cal为通过最小二乘法拟合的萨道夫斯基公式计算的爆破振速峰值；m为现场爆破实验次数。所以随机变量v服从一个数学期望为μ、方差为σ²的正态分布，记作v~(μ, σ²)；相应的分布函数是：

$F\left( v \right) = \frac{1}{{\sqrt {2\pi } \sigma }}\int_{-\infty }^v {{{\text{e}}^{-\frac{{{{\left( {t-\mu } \right)}^2}}}{{2{\sigma ^2}}}}}{\text{d}}t}$

(10)

2. 爆破振动评价和安全炸药量计算

根据随机变量v的分布函数，可以计算振速小于目标设施安全振速的概率，也就是对目标设施在爆破过程中得到保护的情况进行评价，同时根据分布函数，可以推导目标设施满足一定可靠性条件的安全炸药量计算的概率公式。

2.1 振动评价

在概率论与数理统计中，称N(0, 1)为标准正态分布，其分布函数常记为：

$\mathit{\Phi} \left( x \right) = \frac{1}{{\sqrt {2\pi } }}\int_{-\infty }^x {{{\text{e}}^{-\frac{{{t^2}}}{2}}}{\text{d}}t}$

(11)

由于正态分布在概率理论与应用中特殊的重要地位，一般的概率统计著作往往都附有Φ(x)的函数表。而有关任何正态分布N(μ，σ²)的概率计算问题，常常需要借助这些数表来解决。事实上，根据式(11)，设ζ~N(μ，σ²)，在计算P(ξ < b)时，可作如下变换：

$\begin{gathered} P\left( {\xi < b} \right) = F\left( b \right) = \frac{1}{{\sqrt {2\pi } \sigma }}\int_{-\infty }^b {{{\text{e}}^{-\frac{{{{\left( {x-\mu } \right)}^2}}}{{2{\sigma ^2}}}}}{\text{d}}x} = \frac{1}{{\sqrt {2\pi } }}\int_{ - \infty }^b {{{\text{e}}^{ - \frac{{{{\left( {x - \mu } \right)}^2}}}{{2{\sigma ^2}}}}}{\text{d}}\left( {\frac{{x - \mu }}{\sigma }} \right)} \hfill \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;=\frac{1}{{\sqrt {2\pi } }}\int_{ - \infty }^{\frac{{b - \mu }}{\sigma }} {{{\text{e}}^{ - \frac{{{t^2}}}{2}}}{\text{d}}t} = \mathit{\Phi} \left( {\frac{{b - \mu }}{\sigma }} \right)\;\;\;\;\;\;\;\;\;\;\;\;t = \frac{{x - \mu }}{\sigma } \hfill \\ \end{gathered}$

(12)

所以为了便于描述，常将正态变量作数据转换。将一般正态分布转化成标准正态分布。若

$X \sim N\left( {\mu, {\sigma ^2}} \right), \;\;\;\;\;\;\;\;\;\;\;\;Y = \frac{{X-\mu }}{\sigma } \sim N\left( {0, 1} \right)$

服从标准正态分布，通过查标准正态分布表就可以直接计算出原正态分布的概率。设保护目标设施的安全振速为v_saf，单段最大炸药量为Q₀，爆心距R₀，根据式(12)得：

$P\left( {v < {v_{{\text{saf}}}}} \right) = F\left( {{v_{{\text{saf}}}}} \right) = \mathit{\Phi} \left( {\frac{{{v_{{\text{saf}}}}-\mu }}{\sigma }} \right)$

(13)

$\mu = v = K{\left( {\frac{{Q_0^{1/3}}}{{{R_0}}}} \right)^\alpha }, \;\;\;\;\;\;\;\;\;\;\;\sigma = \sqrt {{\sigma ^2}}$

(14)

令：

$\frac{{{v_{{\text{saf}}}}-\mu }}{\sigma } = {m_0}$

(15)

通过查表得：

$\mathit{\Phi} \left( {{m_0}} \right) = {P_0}$

(16)

所以在对爆破设计的振动安全评价过程中，认为单段最大炸药量为Q₀、爆心距R₀处的目标设施得到保护的概率为P₀。

2.2 安全炸药量计算

爆破测振确定振动公式常数的目的是为了用它来控制振速，使它小于目标设施的安全振速v_saf。在已知爆心距R时，只有通过控制炸药量Q来控制振速。使质点振速峰值v小于或等于安全振速v_saf的炸药量称为安全炸药量。设使目标设施得到保护的概率达到P₁，通过查表得到Φ(m₁)=P₁；根据式(13)，得到：

$\mathit{\Phi} \left( {{m_1}} \right) = \Phi \left( {\frac{{{v_{{\text{saf}}}}-\mu }}{\sigma }} \right)$

(17)

即：

${m_1} = \frac{{{v_{{\text{saf}}}}-\mu }}{\sigma }, \;\;\;\;\;\;\;\;\;\;\;\;\mu = {v_{{\text{saf}}}}-{m_1}\sigma$

(18)

根据式(14)，得到：

$\mu = K{\left( {\frac{{{Q^{1/3}}}}{{{R_0}}}} \right)^\alpha }$

(19)

所以：

$Q = {\left( {\sqrt[\alpha]{{\frac{\mu }{K}}}{R_0}} \right)^3}$

(20)

将式(15)代入, 得到：

$Q = {\left( {\sqrt[\alpha]{{\frac{{{v_{{\text{saf}}}} -{m_1}\sigma }}{K}}}{R_0}} \right)^3}$

(21)

通过式(21)，可以求出指定设施安全概率条件下的单段最大炸药量，即使目标设施得到保护的概率^[11]达到P₁的安全炸药量为Q。该式称为计算爆破安全炸药量的概率公式。

3. 应用案例

3.1 工程概况

西部矿业股份公司锡铁山铅锌矿位于青海省柴达木盆地北缘。在铅锌矿井下2 702 m水平1025采场爆破过程中，距爆破点50 m处是运输大巷，里边有一些重要的设备设施。

为了防止由于爆破振动而使巷道损坏，从而造成巷道内设备设施的破坏，需要对此次爆破设计进行振动评价和单段最大炸药量计算。

3.2 振速峰值的分布函数

根据锡铁山铅锌矿之前相同水平其他采场的爆破振动监测数据来进行对1025采场爆破设计，其数据见表 1。表中，Q为最大一段药量，R为爆心距，v为水平径向爆破振动速度峰值。

表 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

| 显示表格

根据表 1，采用线性回归法^[9]拟合萨道夫斯基公式，得到数据拟合曲线图(见图 1)：

图 1 数据拟合曲线

Figure 1. Data fitting curve

下载: 全尺寸图片幻灯片

$v = 210.2{\left( {{Q^{1/3}}/R} \right)^{1.59}}$

(22)

根据式(14)~(15)，对表 1数据进行处理，得到数据处理结果, 见表 2。其中，方差为0.122 cm²/s²。

表 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

| 显示表格

所以随机变量v服从正态分布^[10]，记作v~(v_cal, 0.122)；相应的分布函数是：

$F\left( v \right) = \frac{1}{{\sqrt {2\pi } 0.349}}\int_{-\infty }^v {{{\text{e}}^{-\frac{{{{\left( {t-{v_{{\text{cal}}}}} \right)}^2}}}{{2 \times 0.122}}}}{\text{d}}t}$

(23)

3.3 振动评价和安全药量

根据GB 6722-2014《爆破安全规程》，矿山巷道的安全允许振速15~30 cm/s。为了安全起见，在这里选取其下限作为安全判据，即v_saf=15 cm/s。根据保护目标的重要性，设此次爆破中巷道得到保护的概率应大于95%。通过查表得到：

$\mathit{\Phi} \left( {1.65} \right) = 0.950\;5$

(24)

根据式(21)，取m₁=1.65，σ=0.349，R₀=50 m，v_saf=15 cm/s，求解得出：

$Q = 797.89\;{\text{kg}}$

(25)

那么在爆破设计的过程中，单段最大药量控制在797.89 kg以内，巷道得到保护的概率在95%以上。如果按照萨道夫斯基公式即式(22)计算，求解得出：

$Q = 858.73\;{\text{kg}}$

(26)

即单段最大炸药量控制在858.73 kg以内。此时巷道在爆破过程中得到保护的概率根据式(13)、式(15)~(16)计算：

$P\left( {v < {v_{{\text{saf}}}}} \right) = F\left( {15} \right) = \mathit{\Phi} \left( 0 \right)$

(27)

通过查表得：

$\mathit{\Phi} \left( 0 \right) = 0.50$

(28)

即在爆破过程中，巷道得到保护的概率只有50%。

在萨道夫斯基公式和正态分布函数公式计算下，单段最大药量和爆破后设施得到保护概率见表 3。可以看出，用萨道夫斯基公式计算出的单段最大炸药量比概率公式多60 kg左右，但是根据概率公式计算出的单段最大炸药量，其可靠性在95%以上，比应用萨道夫斯基公式计算值的可靠性(仅为50%)要高的多。因此，根据之前爆破震动监测数据及分析结果，建议采用式(21)的概率公式来计算安全药量，以保证施爆期间1025采场运输大巷的安全。

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

Table 3. Single biggest dosage andfacilities security probability

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

下载: 导出CSV

| 显示表格

4. 结论

通过对正态分布函数分析、爆破振动评价与安全药量计算，现得出如下结论：

(1) 岩土内部含有大量复杂的裂隙，岩石内部是不连续和不均质性的，爆破引起的质点振动的实际振速峰值v是随机变量，v=K(Q^1/3/R)^α是实际测点质点振动速度峰值的期望值或估计值。

(2) 对任一确定的Q^1/3/R，爆破实际振速峰值v小于目标设施安全振速的概率可以计算出，即对保护对象的振动安全性可以作出评价。

(3) 为了保证使目标设施得到保护，有较高的可靠性，安全炸药量的计算要应用概率公式。

(4) 该计算方法用于爆破的振动安全评价和安全药量计算是可行的、合理的，为爆破振动安全评价和安全药量计算提供了新的方法。

参考文献(0)

施引文献

期刊类型引用(15)

1.	梁富，张鹏飞，何运华，孙雪. 基于综合赋权的矿岩可爆性分级识别及炸药单耗预测研究. 煤. 2025(02): 48-53+92 . 百度学术
2.	孙鹏昌，杨招伟，卢文波，孟海利，薛里. 白鹤滩水电站坝肩边坡开挖爆破振动概率分布研究. 工程科学与技术. 2024(03): 179-188 . 百度学术
3.	王剑辉. 宏岩煤矿10105工作面大巷煤柱切顶卸压技术. 能源与节能. 2024(07): 201-204 . 百度学术
4.	张文涛，汪海波，高朋飞，王梦想，程兵，宗琦. 大型露天矿山深孔台阶爆破振动特征与安全控制. 金属矿山. 2024(09): 151-160 . 百度学术
5.	李启月，李丽，黄海仙，肖宇航，魏新傲. 基于概率论的高边坡爆破振动安全预测及控制研究. 中国安全生产科学技术. 2024(12): 12-18 . 百度学术
6.	李辉，王董东，杜文康，弓煜. 安家岭煤矿爆破振速预测与安全允许距离研究. 煤炭工程. 2023(10): 156-161 . 百度学术
7.	黄聪，雷振，韦善阳，高正华，雷兴海. 基于AHP和正态分布的高程爆破振动评价与预测. 工程爆破. 2022(02): 121-127+134 . 百度学术
8.	贾海鹏，刘殿书，方真刚，田帅康. 地铁隧道钻爆法施工中敏感区间及安全药量确定. 北京理工大学学报. 2021(01): 23-29 . 百度学术
9.	钱安康，陈志敏，朱煊，文勇. 短天窗点间隔下邻近既有隧道爆破方案的优化. 工程爆破. 2020(01): 36-42 . 百度学术
10.	吴波，兰扬斌，杨建新，韩亚龙. 基于AHP-GRA新建隧道爆破对邻近铁路隧道振动影响研究. 铁道科学与工程学报. 2020(03): 668-675 . 百度学术
11.	孙赛赛，池恩安，牛国廷，张修玉，明悦，雷振. 标准正态分布在高程爆破振动评价与预测中的应用. 矿业研究与开发. 2020(04): 35-40 . 百度学术
12.	王林峰，邓冰杰，莫诎，赵精富，肖弘光. 基于概率论的爆破振动安全评估与控制. 振动与冲击. 2020(14): 122-129 . 百度学术
13.	魏雯. 隧道爆破施工的关键技术浅析. 江西建材. 2020(10): 159-160+162 . 百度学术
14.	管晓明，聂庆科，李华伟，安建永. 隧道爆破振动下既有建筑结构动力响应及损伤研究综述. 土木工程学报. 2019(S1): 151-158 . 百度学术
15.	张伟，艾广建，陈仁山. 浅埋地下爆炸的地面运动研究. 山东交通学院学报. 2019(04): 54-59+77 . 百度学术

其他类型引用(7)

资源附件(0)

访问统计

计量

文章访问数: 77
HTML全文浏览量: 85
PDF下载量: 42
被引次数: 22

1. 随机变量v的分布函数
1.1 线性回归法确定K、α
1.2 正态分布函数
2. 爆破振动评价和安全炸药量计算
2.1 振动评价
2.2 安全炸药量计算
3. 应用案例
3.1 工程概况
3.2 振速峰值的分布函数
3.3 振动评价和安全药量
4. 结论