一种基于试爆数据预测水下爆破装药量的方法与流程

1.本发明涉及一种基于试爆数据预测水下爆破装药量的方法，属于工程爆破领域。


背景技术：

2.炸药在介质中爆炸时必然有一部分能量转换为地震波，并从爆源向周围介质迅速传播，导致地表和建(构)筑物产生震动，可能给附近非爆破目标的建(构)筑物造成不同程度上的破坏，这种现象称之为爆破震动效应。随着国民经济建设的迅猛发展，爆破作业得到广泛应用，相应的爆破震动效应带来的危害引起了人们的高度重视。
3.然而，工程实践中，爆破工程技术人员对爆破施工方案可能带来的震动危害评估主要通过以下两个措施来判定：(1)引用学者萨道夫斯基根据牛顿力学导出的经验公式进行核算，因核算公式中爆破点至被保护对象间的地形、地质条件有关的系数k和衰减指数α处于较大的范围，核算时参数值的选择完全依赖技术人员个人经验，缺乏明确的依据，使得核算结果合理性难以保障；(2)通过试爆进行实验测试以观察爆破效果及爆破危害，但通过试爆结果推算施工方案的爆破震动危害完全依赖于个人经验，缺乏试爆效果观察数据与爆破施工方案预期效果之间的科学关联，使得依据试爆测试数据反馈于爆破施工方案设计科学、合理性难以保障。
4.另外,以上两种措施因简便、易操作而被工程爆破行业广泛应用，且在爆破拆除领域应用效果良好，但在水下爆破领域带来众多困惑。特别的，爆破震动危害带来的损失与纠纷使得爆破作业常常面临难以实施的困境，极大的影响了正常的工程建设。


技术实现要素：

5.发明目的：为了克服现有技术中存在的不足，本发明提供一种基于试爆数据预测水下爆破装药量的方法，通过试爆数据的收集归纳整理，反推出预测公式模型的关键系数，为科学爆破提供依据。
6.技术方案：为解决上述技术问题，本发明的一种基于试爆数据预测水下爆破装药量的方法,先进行试爆，得到试爆的震动数据；所述震动数据包括爆破震动强度和爆破震动频率；然后结合试爆方案，基于爆破震动传播行为分析方法反推出爆破震动传播行为和衰减规律，从而得到爆破点与监测点之间的爆破震动强度和爆破震动主频衰减公式，具体的包括如下步骤:
7.步骤1、先进行试爆，用爆破测振仪对水下爆破近邻域堤岸建(构)筑物保护对象进行爆破检测从而得到震动强度v，以及爆破齐发炸药量q、监测点距爆破点空间距离r、爆破震动水域传播距离r
d
；
8.步骤2、将步骤1中得到的数据v、q、r，和r
d
代入以下预测模型公式(1)，通过采用线性最小二乘法回归拟合得到所述预测模型公式(1)中的待定参数k、β1和β2的数值：
[0009][0010]
式中，k表示的是场地系数，β1表示的是衰减系数，β2表示的是水土介质变化导致的衰减系数；
[0011]
步骤3、根据水下爆破近邻域堤岸建/构筑物保护对象特征，查阅《爆破安全规程》(gb6722
‑
2014)中该保护类别中v的最大值v
max
，将v
max
代入公式(2)中获得最大爆破齐发炸药量q
max
，为工程爆破进行指导：
[0012][0013]
进一步地，步骤一、由于最大许可爆破震动速度v与爆破震动主频相关联，将步骤3中得到的q
max
代入以下公式(3)获得f
max
:
[0014][0015]
步骤二、由于已知f
max
，根据查表可知在该保护类别中v
max’，将v
max’代入公式(2)中，从而得到q
max’。
[0016]
比较q
max
和q
max’，选取最小值作为该爆破作业的最大爆破齐发炸药量。
[0017]
步骤1～3已经求得了最大爆破齐发炸药量，通过步骤一和步骤二的进一步推算，可进一步计算出更为精确的最大爆破齐发炸药量。
[0018]
由许可爆破震动速度推算的爆破齐发药量q导出水下爆破震动传播至近邻域堤岸建(构)筑物保护对象的主频f。通过查阅《爆破安全规程》 (gb6722
‑
2014)，若此主频值f与许可爆破震动速度v于《爆破安全规程》 (gb6722
‑
2014)中对应相同的频率范围，则说明许可爆破震动速度v值选取合理；相应的，推算出的爆破齐发药量q可信，完成一个爆破参数方案设计逻辑。
[0019]
综上可知,(1)、(2)、(3)式形成一个闭环的计算逻辑。首先通过试爆数据代入(1)式拟合出待定参数k、β1和β2的数值；再根据水下爆破近邻域受保护建(构)筑物特征，参考《爆破安全规程》(gb6722
‑
2014) 查出各频率范围许可爆破震动速度v；再由(2)式推算出爆破齐发炸药量 q；进一步，将爆破齐发炸药量q代入(3)式，导出对应的爆破主频f；若此主频值f与许可爆破震动速度v于《爆破安全规程》(gb6722
‑
2014) 中对应相同的频率范围，则完成一个爆破参数方案方法。若要进一步精确爆破齐发炸药量，则再通过主频f
max
推算得到q
max’，比较q
max
和q
max’，选取最小值作为该爆破作业的最大爆破齐发炸药量。可见步骤一和步骤二是对步骤1～步骤3的进一步验算，为了更精确。
[0020]
水下爆破对堤岸近邻域结构的震动激励，较之于岩土爆破对其近邻域结构的震动激励而言，存在两个显著的差异：一是因流体状态的水介质不能承受剪切荷载而无法传递横波(又称为剪切波)，使得水下爆破对堤岸近邻域结构震动激励天然地过滤掉了横波成分，而主要呈现为纵波激励；二是因水介质与岩土介质波阻抗的差异，使得来自于水下爆破的纵波在水土界面上发生入射、反射和透射，具体的入射、反射、透射波强度与水、岩土介质波阻抗相关。总体而言，水中爆破对堤岸近邻域结构震动激励行为较之于岩土爆破而言，波在水介质及其与岩土界面上传播行为呈现显著差异。
[0021]
有益效果：(1)由于介质的差异，本发明针对水下爆破对堤岸近邻域结构震动激
励，认识到由于水、土介质受载响应行为及波阻抗差异，预测水下爆破对堤岸近邻域结构震动激励行为不能再直接引用当前工程实践中广泛采用的岩土爆破震动校核公式，而需计及水介质不能传递横波，且纵波在水土界面会发生入射、反射、透射的强度重新分配。因此，本发明通过量纲分析与最小二乘法回归拟合，在岩土爆破震动预测公式基础上补充了一个算子，见公式(1)右边第三项，以表征水土介质差异，并将其融入于预测公式中，使得水下爆破震动对堤岸近邻域结构震动激励预测公式得到相应改进；(2)本发明通过进一步验算，获得更为精确的爆破齐发炸药量，从而为水下工程爆破提供更为安全的指导；(3)本发明预测方法简单，通过实验测试以观察爆破效果及爆破危害，并通过试爆结果推算施工方案的爆破震动危害，让试爆效果观察数据与爆破施工方案预期效果之间具有的科学关联，有据可依；(4)本发明专门针对水下爆破做出改进，通过理论分析得到水中爆破对堤岸近邻域结构震动的激励行为，为后续研究提供基础。
附图说明
[0022]
图1为本发明基于试爆数据预测爆破施工方案的逻辑框图。
具体实施方式
[0023]
下面结合附图对本发明作更进一步的说明。
[0024]
实施例
[0025]
本发明的一种基于试爆数据预测水下爆破装药量的方法,先进行试爆，得到试爆的震动数据；震动数据包括爆破震动强度和爆破震动频率；然后结合试爆方案，基于爆破震动传播行为分析方法反推出爆破震动传播行为和衰减规律，从而得到爆破点与监测点之间的爆破震动强度和爆破震动主频衰减公式，具体的包括如下步骤:
[0026]
步骤1、先进行试爆，用爆破测振仪对水下爆破近邻域堤岸建(构)筑物保护对象进行爆破检测从而得到震动强度v，以及爆破齐发炸药量q、监测点距爆破点空间距离r、爆破震动水域传播距离r
d
；
[0027]
步骤2、将步骤1中得到的数据v、q、r，和r
d
代入以下预测模型公式(1)，通过采用线性最小二乘法回归拟合得到所述预测模型公式(1)中的待定参数k、β1和β2的数值：
[0028][0029]
式中，k表示的是场地系数，β1表示的是衰减系数，β2表示的是水土介质变化导致的衰减系数；
[0030]
步骤3、根据水下爆破近邻域堤岸建/构筑物保护对象特征，查阅《爆破安全规程》(gb6722
‑
2014)中该保护类别中v的最大值v
max
,将v
max
代入公式(2)中获得最大爆破齐发炸药量q
max
，为工程爆破进行指导：
[0031][0032]
步骤一、由于最大许可爆破震动速度v与爆破震动主频相关联，将步骤3中得到的q
max
代入以下公式(3)获得f
max
:
[0033][0034]
步骤二、由于已知f
max
，根据查表可知在该保护类别中v
max’，将v
max’代入公式(2)中，从而得到q
max’。
[0035]
比较q
max
和q
max’，选取最小值作为该爆破作业的最大爆破齐发炸药量。
[0036]
步骤1～3已经求得了最大爆破齐发炸药量，通过步骤一和步骤二的进一步推算，可进一步计算出更为精确的最大爆破齐发炸药量。
[0037]
通过查阅《爆破安全规程》(gb6722
‑
2014)，如表2所示，若此主频值 f与许可爆破震动速度v在表2中对应于相同的频率范围，则说明许可爆破震动速度v值选取合理。
[0038]
表2爆破振动安全允许标准
[0039][0040]
预测模型公式(1)是在考虑水下爆破震动在堤岸及其近邻域传播所伴随的水、土介质变化特征，通过筛选与爆破震动主频有关的物理量并运用量纲分析推导出水下爆破爆破震动强度泛函关系式，具体推导过程如下：
[0041]
研究水下爆破震动在堤岸及近邻域传播强度衰减规律时，由于震动传播介质涉及到水、土中，故得考虑水、土介质变化特征对爆破震动传播强度的影响。本发明采用量纲分析的方法推导震动强度衰减公式。现将影响爆破震动强度的自变量给出，分别是：炸药量q，爆心距r，测点到堤岸的距离r
d
，介质密度ρ，介质中震动波传播速度c以及爆轰时间t。爆破震动位移μ、震动速度v、震动加速度a以及震动主频f作为因变量，与其影响的变量可以根据布金汉定理(π定理)得到以下函数式：
[0042]
v＝f(q,r,c,μ,a,f,ρ,r
d
,t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0043]
由式(5)可知，分析问题的物理量总数目为n＝10。根据π定理，在自 变量中取基本量为q,r,c，所以基本量的量纲数目为m＝3，导出量与因变量的 量纲数目为n
‑
m＝7，π、π1～π7代表的是无量纲变量，则有：
[0044][0045]
将式(6)代入式(5)有：
[0046][0047]
根据对因变量π不起作用的无量纲变量可记为常数1，则有：
[0048][0049]
对于爆破对象场地介质属性，介质密度ρ与震动波传播速度c均可记为常数，由式(8)可进一步化简得：
[0050][0051]
考虑爆破震动与炸药量q，爆源距r，测点与堤岸的距离r
d
的关系，可将式(9)写成：
[0052][0053]
令则有：
[0054]
lnv0＝α1 β
′1lnq
‑
3β
′1lnr
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0055]
其中，α1 β
′1lnq表示的是爆破对象场地的介质条件与炸药量对震动质点速度的综合影响；
‑
β
′1lnr表示的是震动质点速度与爆源距r的衰减关系；β
′1表示的是衰减系数，主要反映场地介质条件对震动质点速度的影响。
[0056]
令lnk1＝α1，β1＝
‑
3β
′1，则有：
[0057][0058]
其中，表示的是露天爆破条件下的萨道夫斯基公式。
[0059]
将式(12)代入式(10)有：
[0060][0061]
令lnk2＝
‑
ε1，β2＝
‑
η1，则在考虑水域因素的条件下，通过引入距离比值系数的爆破震动质点速度的计算公式为：
[0062][0063]
其中，k1表示的是爆破对象场地系数；k2表示的是地形地貌影响系数；β1表示的是衰减系数；β2表示的与距离比值系数；
[0064]
在式(14)中，可以用系数k来代替k1,k2，有：
[0065]
[0066]
其中，在步骤2中，根据实测爆破震动数据，采用线性最小二乘法对预测模型公式中的待定参数进行回归拟合，具体拟合过程为：
[0067]
对式(2)等号两边取对数，有：
[0068][0069]
令y＝lnv，β0＝lnk，x2＝ln(r
d
/r)，则式(16)可化为：
[0070]
y＝β0 β1x1 β2x2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0071]
为了求出待定参数β0,β1,β2的最优估计值，需要使得震动实测数据y
i
与式 (17)拟合值的残差平方和最小，即求解目标函数：
[0072][0073]
由多元函数的极值定理知，f取最小值的必要条件：
[0074][0075]
那么，可求得非齐次线性方程组解β0，β1，β2。再通过变量代换可求得那么具体的水下爆破装药量表达式为：
[0076][0077]
其中，在步骤3中计算模型公式为：
[0078][0079]
式(21)中，根据预测点最大允许的震动速度值v，可以计算出爆破齐发炸药量q，从而指导工程爆破。
[0080]
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。