Parameter Optimization of Presplitting Blasting Based on Model Test
-
摘要: 为了有效地降低预裂爆破振动对边坡岩体和周边建构筑物的扰动,除了关注预裂缝的减振效应外,还需对预裂孔本身的爆破参数和起爆方式进行优化。为此,采用正交试验法,选用C30混凝土为爆破对象的相似模拟材料,开展预裂孔爆破成缝过程的相似模型试验,选取不耦合系数、延期时间、最大单响药量3个因素,每个因素设置4个水平,以有效半孔率、预裂缝宽度和原岩损伤率为评价指标,通过极差与方差计算分析各因素对评价指标的敏感度,确定了本次模型试验达到最佳爆破效果的预裂孔爆破参数:不耦合系数为1.33,延期时间为12 ms,最大单响药量为1.8 g。研究结果可为现场预裂孔爆破参数选择和精确延期时间设计提供指导。Abstract: In order to effectively reduce the disturbance derived from the vibration of the presplitting blasting in the slope rock mass and surrounding structures, it is important not only to pay attention to the vibration damping effect of the presplitting, but also to optimize the blasting parameters and initiation modes of the presplitting hole. The study adopts the orthogonal experimental method, which selects the C30 concrete as the similar simulation material of the blasting object to carry out the similar model test of the blasting process with presplitting hole. In this model test, three factors are considered as the non-coupling coefficient, which are delay time, and maximum single-shot charge, and four levels are set for each factor, taking effective half porosity, pre-fracture width and original rock damage rate as evaluation indicators. After the analysis of the sensitivity of each factor to evaluation indicators through range and variance calculations, the best blasting effect with the blasting parameters of the presplitting hole are finally achieved in this model test: non-coupling coefficient is 1.33, delay time is 12 ms, and maximum single-shot charge is 1.8 g. The research results provide guidance for the design of the parameters of on-site blasting with presplitting hole and precise delay time.
-
近年来,随着社会对爆破安全的要求越来越严苛,受地质条件、工程条件和技术条件等限制,为了降低爆破对周边岩体的振动破坏,保障周边岩体和建筑的安全性,在施工中广泛采用预裂爆破技术。研究预裂孔参数对预裂爆破技术应用具有重要意义[1-3]。邓珂等[4]利用数值模拟软件LS-DYNA,对初始地应力对预裂缝成缝的影响进行了研究。朱强等[5]通过数值模方法研究了空气间隔装药预裂爆破对岩体的损伤影响,提出采用径向不耦合药卷和反向起爆可以优化空气间隔装药预裂爆破效果。He等[6]应用数值模拟软件LS-DYNA优化了预裂爆破模型,并将其应用于现场,使有效半孔率超过20%,减振率达28.2%。杨仁树等[7]基于数字激光动态焦散线,研究了预裂缝对爆生裂纹和原生裂纹动态断裂特性的影响。璩世杰等[8]进行了不同角度的节理对预裂爆破效果的试验,得到了不同角度下预裂缝成缝效果。陈俊桦等[9]结合岩石爆破理论、爆生气体和应力波共同作用等理论,提出了基于岩石损伤的预裂爆破参数计算公式。Ma等[10]对预裂爆破的成缝机理进行了探讨,得到了预裂孔孔距的计算公式,并验证了预裂缝的减振效果。饶宇等[11]通过分析预裂孔爆破和主爆破的功率谱和能量谱,得出预裂孔爆破形成的预裂缝能很好地阻挡主爆破能量的传播,揭示了预裂缝对爆破振动传播及其频谱特征变化的影响。严鹏等[12]通过对爆破开挖振动监测和爆破损伤的声波监测,得到了预裂爆破的振动衰减规律。Xiao等[13]结合爆炸力学,揭示了预裂爆破的开裂机理,提出了一套预裂爆破参数设计公式,最终确定了黑岱沟的预裂爆破参数。
以往学者主要通过数值模拟、室内实验及理论推导等手段对预裂爆破的影响因素开展了研究,然而针对预裂孔参数的模型试验研究却较少。本研究拟采用相似材料模拟预裂孔爆破过程,通过正交分析确定爆破成缝效果最佳的预裂孔爆破参数。
1. 试验准备与设计
(1) 模型的制备。本试验的相似模拟对象为云南某露天矿山,现场岩体主要为片岩、大理岩体、矽卡岩。根据相似理论,选用C30混凝土作为相似材料,其水、水泥、砂、石子的配合比为0.60∶1.00∶2.51∶3.93,平均密度为2.382 g/cm3,抗压强度为27.6 MPa,弹性模量为27.6 GPa。一共浇筑16组模型,通过插入钢筋预制炮孔,如图1所示。为了增加约束,爆破时不脱模。
(2) 不耦合系数(因素A)水平。徐颖等[14]的模型试验研究表明,当不耦合系数为1.67时,爆炸裂纹数和裂纹扩展长度最佳,为此选用1.33、1.67、2.00和3.00这4个水平,分别对应A1、A2、A3、A4。
(3) 延期时间(因素B)水平。何理等[15]通过浇筑混凝土模型模拟预裂爆破实验,确定孔间最佳延期时间为11~15 ms,司剑峰等[16]通过模型研究得出最佳延期时间约为8 ms,因此本研究中延期时间选用8、10、12、15 ms 4个水平,分别对应B1、B2、B3、B4。
(4) 最大单响药量(因素C)水平。本次模型试验选用数码电子雷管,导爆索为主爆药,模型试验的线装药密度为2 g/m,孔深为80 cm,因此单孔装药量为1.6 g,选用1.6、3.6、5.4 g,分别对应C1、C2、C3。考虑到因素C只有3个水平,为此采用拟水平设计将因素C补充至4个水平,用3.6*表示为因素C的拟最大单响药量3.6 g水平,即C4。孔距统一为90 mm,填塞长度为1.2倍孔径。考虑误差而增设误差项(D),对试验没有影响。正交试验因素水平见表1。
表 1 正交试验因素水平Table 1. Factors and levels for orthogonal testsLevel Factors A B/ms C/g D 1 1.33 8 1.8 1 2 1.67 10 3.6 2 3 2.00 12 5.4 3 4 3.00 15 3.6* 4 本次模型试验的爆后效果用预裂缝宽度(d)、有效半孔率(δ)和原岩损伤率(η)来评价。预裂缝宽度采用卷尺直接测量;有效半孔率为爆破后较完整半孔长度与炮孔总深度的比值;原岩损伤率则用相同部位的爆破后波速(
Cp2 )与爆破前波速(Cp1 )表示,即η=1−Cp2Cp1 (1) 2. 正交试验分析
本试验设有3个因素,每个因素均设计了4种水平。需采用4因素4水平的正交表
L16(44) ,共计有16组试验,不同因素组合试验的正交表见表2。表 2 相似材料模拟试验的正交设计Table 2. Orthogonal test design of the simulation tests with similar materialTest No. Factors A B/ms C/g D 1 1.33 8 1.8 1 2 1.33 10 3.6 2 3 1.33 12 5.4 3 4 1.33 15 3.6* 4 5 1.67 8 3.6 4 6 1.67 10 1.8 3 7 1.67 12 3.6* 2 8 1.67 15 5.4 1 9 2.00 8 5.4 2 10 2.00 10 3.6* 1 11 2.00 12 1.8 4 12 2.00 15 3.6 3 13 3.00 8 3.6* 3 14 3.00 10 5.4 4 15 3.00 12 3.6 1 16 3.00 15 1.8 2 根据设计的16组方案进行相似模型试验,有效半孔率、预裂缝宽度和原岩损伤率的测试过程如图2、图3和图4所示,具体的试验结果列于表3。
表 3 爆破效果统计Table 3. Statistics of blasting effectTest No. Evaluating indicator Test No. Evaluating indicator δ/% d/mm η/% δ/% d/mm η/% 1 45 51 14.0 9 65 16 24.8 2 48 63 14.5 10 76 15 21.4 3 51 80 19.3 11 70 10 10.0 4 41 54 13.1 12 65 12 10.6 5 56 41 19.4 13 76 12 12.3 6 30 45 14.0 14 90 21 13.5 7 78 62 16.2 15 87 17 12.8 8 62 51 19.8 16 76 15 12.0 2.1 极差分析
Ki 为该列因素的第i个水平所对应分析的评价指标数值之和。ki为该列因素Ki具有的水平数的均值。极差R为该列k1、k2、k3、k4中的最大值减去最小值。极差越大,则该因素对应的评价指标的影响越显著。
根据16组试验的试验结果,对有效半孔率、预裂缝宽度和原岩损伤率下的各因素进行极差分析,见表4。由表4可知,当只考虑有效半孔率时,因素A的ki指标有k1 < k2 < k3 < k4,可以判断A4为只考虑有效半孔率评价指标下的最优水平。同理,可得因素B和因素C在只考虑有效半孔率时的最优水平,分别为B3、C3;只考虑预裂缝宽度时,因素A、B、C的最优水平分别为A1、B2、C2;只考虑原岩损伤时,因为原岩损伤程度越低,对周围岩体的保护越好,所以因素A、B、C的最优水平分别为A4、B4、C1。
表 4 极差分析Table 4. Range analysisRange analysis
indexFactors Range analysis of δ Range analysis of d Range analysis of η A B C D A B C D A B C D K1 185.0 242.0 231.0 270.0 248.0 120.0 121.0 134.0 60.9 70.5 50.0 68.0 K2 226.0 244.0 266.0 277.0 199.0 144.0 133.0 156.0 69.4 63.4 57.3 67.5 K3 286.0 286.0 268.0 232.0 53.0 169.0 168.0 149.0 66.8 58.3 77.4 56.2 K4 339.0 264.0 271.0 257.0 65.0 132.0 143.0 126.0 50.6 55.5 63.0 56.0 k1 46.3 60.5 57.8 67.5 62.0 30.0 30.3 33.5 15.2 17.6 12.5 17.0 k2 56.5 61.0 66.5 69.3 49.8 36.0 33.3 39.0 17.4 15.8 14.3 16.9 k3 71.5 71.5 67.0 58.0 13.3 42.3 42.0 37.3 16.7 14.6 19.4 14.1 k4 84.8 66.0 67.8 64.3 16.3 33.0 35.8 31.5 12.7 13.9 15.8 14.0 R 38.5 11.0 9.3 11.3 33.5 12.3 11.8 7.5 4.7 3.8 6.9 3.0 对比分析极差可知影响预裂爆破效果的各个因素的主次关系。对于有效半孔率,因素A的极差最大,说明不耦合系数对预裂爆破效果中有效半孔率的影响程度最大,因素B、因素C依次减小,各个因素的影响程度为:不耦合系数的影响最大,延期时间次之,最大单响药量的影响最小。对于预裂缝宽度,因素A的极差最大,说明不耦合系数对预裂爆破效果中预裂缝宽度的影响程度最大,因素B、因素C依次减小,各个因素的影响程度为:不耦合系数的影响最大,延期时间次之,最大单响药量的影响最小。对于原岩损伤程度,因素C的极差最大,说明最大单响药量对预裂爆破效果中原岩损伤程度的影响程度最大,因素A、因素B依次减小,各个因素的影响程度为:最大单响药量的影响最大,不耦合系数 次之,延期时间的影响最小。
2.2 方差分析
方差分析是考察各因子在各水平下的效益有无显著差异的一种分析方法,与极差分析相互验证补充。在方差分析中偏差平方和QT的表达式为
QT = 1nn∑i=1K2i−T2mn (2) 式中:n为正交试验水平数,m为水平重复次数,T为所分析评价指标数值的总和。设μ为自由度,Si为均方差,ST为误差均方差,则
μ=n−1 (3) Si=QTμ (4) F=SiST (5) 通过计算得出F后,与选定置信度
α 的F分布中的定值进行比较,当F > Fα 时,该因素对该评价指标影响显著。 根据分析得到统计结果,如表5所示。最大单响药量仅有3个水平,故计算时自由度为2。
表 5 方差分析结果Table 5. Results of variance analysisFactors Variance analysis of δ Variance analysis of d Variance analysis of η QT μ Si F QT μ Si F QT μ Si F A 3423.5 3 1141.2 11.60 7083.2 3 2361.1 50.35 52.2 3 17.4 1.54 B 317.0 3 105.7 1.08 328.7 3 109.6 2.34 32.5 3 10.8 0.96 C 261.3 2 130.7 1.33 286.7 2 143.3 3.06 97.0 2 48.5 4.28 Error 294.5 3 98.2 140.7 3 46.9 34.0 3 11.3 查询F分布表可知,F0.10(2, 3) = 5.46,F0.05(2, 3) = 9.55,F0.01(2, 3) = 30.82,F0.10(3, 3) = 5.39,F0.05(3, 3) = 9.28,F0.01(3, 3) = 29.46。
根据表5,对于预裂缝宽度,因素A的F > F0.01(3, 3),对于有效半孔率,因素A的F > F0.05(3, 3),与极差分析结果基本一致。其他两个因素虽然没有达到显著水平,但是总体上仍反映并证明了极差分析结果。
对于有效半孔率和原岩损伤程度,3个因素不显著的原因是本试验的误差相对较大,而且误差的自由度相对较小,使得检验的灵敏度较低。
2.3 最优水平确定
(1) 因素A在有效半孔率指标和预裂缝宽度中的影响程度排第1,为最重要因素;在对原岩的损伤程度指标中排第2,在方差分析中,因素A 在预裂缝宽度情况下F > F0.01(3, 3),因素A在预裂宽度情况下达到显著水平,故可取A1水平,即1.33。
(2) 因素B在有效半孔率指标和预裂缝宽度中的影响程度排第2,为次要因素;在对原岩的损伤程度指标中排第3,为最不重要因素;故可取水平B2和水平B3,因为在因素B取B2时,3种指标总值分别为244、144、63.4,比B3时的286、169、58.6指标差,故取B3水平,即12 ms。
(3) 因素C在有效半孔率指标和预裂缝宽度中影响程度排第3,为最不重要因素;在对原岩的损伤程度指标中排第1,为最重要因素,故应取C1水平,即1.8 g。
3. 结 论
通过C30混凝土相似材料的爆破模型试验,对影响预裂孔爆破成缝效果的多因素进行了正交试验分析,得到以下结论。
(1) 方差分析表明,预裂缝宽度对不耦合系数最敏感,且不耦合系数越大,预裂缝宽度越小;有效半孔率与不耦合系数成正比,说明在药量一定的情况下,增大不耦合系数可提高对炮孔壁的保护。
(2) 延期时间对有效半孔率和预裂缝宽度的影响程度均排第2,说明延期时间对爆破成缝和孔壁的保护均有影响,因此选择合理的延期时间对爆破效果至关重要。
(3) 模型试验爆破效果评价表明,不耦合系数取1.33、延期时间取12 ms、最大单响药量取1.8 g的因素组合条件下,预裂爆破成缝和减振效果最优。
-
表 1 正交试验因素水平
Table 1. Factors and levels for orthogonal tests
Level Factors A B/ms C/g D 1 1.33 8 1.8 1 2 1.67 10 3.6 2 3 2.00 12 5.4 3 4 3.00 15 3.6* 4 表 2 相似材料模拟试验的正交设计
Table 2. Orthogonal test design of the simulation tests with similar material
Test No. Factors A B/ms C/g D 1 1.33 8 1.8 1 2 1.33 10 3.6 2 3 1.33 12 5.4 3 4 1.33 15 3.6* 4 5 1.67 8 3.6 4 6 1.67 10 1.8 3 7 1.67 12 3.6* 2 8 1.67 15 5.4 1 9 2.00 8 5.4 2 10 2.00 10 3.6* 1 11 2.00 12 1.8 4 12 2.00 15 3.6 3 13 3.00 8 3.6* 3 14 3.00 10 5.4 4 15 3.00 12 3.6 1 16 3.00 15 1.8 2 表 3 爆破效果统计
Table 3. Statistics of blasting effect
Test No. Evaluating indicator Test No. Evaluating indicator δ/% d/mm η/% δ/% d/mm η/% 1 45 51 14.0 9 65 16 24.8 2 48 63 14.5 10 76 15 21.4 3 51 80 19.3 11 70 10 10.0 4 41 54 13.1 12 65 12 10.6 5 56 41 19.4 13 76 12 12.3 6 30 45 14.0 14 90 21 13.5 7 78 62 16.2 15 87 17 12.8 8 62 51 19.8 16 76 15 12.0 表 4 极差分析
Table 4. Range analysis
Range analysis
indexFactors Range analysis of δ Range analysis of d Range analysis of η A B C D A B C D A B C D K1 185.0 242.0 231.0 270.0 248.0 120.0 121.0 134.0 60.9 70.5 50.0 68.0 K2 226.0 244.0 266.0 277.0 199.0 144.0 133.0 156.0 69.4 63.4 57.3 67.5 K3 286.0 286.0 268.0 232.0 53.0 169.0 168.0 149.0 66.8 58.3 77.4 56.2 K4 339.0 264.0 271.0 257.0 65.0 132.0 143.0 126.0 50.6 55.5 63.0 56.0 k1 46.3 60.5 57.8 67.5 62.0 30.0 30.3 33.5 15.2 17.6 12.5 17.0 k2 56.5 61.0 66.5 69.3 49.8 36.0 33.3 39.0 17.4 15.8 14.3 16.9 k3 71.5 71.5 67.0 58.0 13.3 42.3 42.0 37.3 16.7 14.6 19.4 14.1 k4 84.8 66.0 67.8 64.3 16.3 33.0 35.8 31.5 12.7 13.9 15.8 14.0 R 38.5 11.0 9.3 11.3 33.5 12.3 11.8 7.5 4.7 3.8 6.9 3.0 表 5 方差分析结果
Table 5. Results of variance analysis
Factors Variance analysis of δ Variance analysis of d Variance analysis of η QT μ Si F QT μ Si F QT μ Si F A 3423.5 3 1141.2 11.60 7083.2 3 2361.1 50.35 52.2 3 17.4 1.54 B 317.0 3 105.7 1.08 328.7 3 109.6 2.34 32.5 3 10.8 0.96 C 261.3 2 130.7 1.33 286.7 2 143.3 3.06 97.0 2 48.5 4.28 Error 294.5 3 98.2 140.7 3 46.9 34.0 3 11.3 -
[1] MELNIKOV N V. Charge construction influence on explosion operations efficiency [C]//Reports of the Science Symposium on Drilling, Explosives, Explosion Operations and Study of Physical and Mechanical Properties of Rocks. Rolla, USA: Gosgortekhizdat, 1962: 35−38. [2] FOURNEY W L, BARKER D B, HOLLOWAY D C. Model studies ofexplosive well stimulation techniques [J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 1981, 18(2): 113–127. doi: 10.1016/0148-9062(81)90737-3 [3] PALROY P, SINGH R B. Air-decking principle and its applications inproduction and presplit blasting [J]. CIM Bulletin, 2001, 94(1): 60–64. [4] 邓珂, 陈明, 卢文波, 等. 地应力对坝肩槽预裂爆破成缝的影响研究 [J]. 岩土力学, 2019, 40(3): 1121–1128.DENG K, CHEN M, LU W B, et al. Investigation of influence of in-situ stress on presplitting induced fracture in abutment slot [J]. Rock and Soil Mechanics, 2019, 40(3): 1121–1128. [5] 朱强, 陈明, 郑炳旭, 等. 空气间隔装药预裂爆破岩体损伤分布特征及控制技术 [J]. 岩石力学与工程学报, 2016, 35(Suppl 1): 2758–2765.ZHU Q, CHEN M, ZHENG B X, et al. Distribution and control technology of rock damage induced by air-deck charge presplitting blasting [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(Suppl 1): 2758–2765. [6] HE L H, WANG J G, XIAO J Q, et al. Pre-splitting blasting vibration reduction effect research on weak rock mass [J]. Disaster Advances, 2013, 6: 338–343. [7] 杨仁树, 苏洪. 爆炸荷载下含预裂缝的裂纹扩展实验研究 [J]. 煤炭学报, 2019, 44(2): 482–489.YANG R S, SU H. Experimental study on crack propagation with pre-crack under explosion load [J]. Journal of China Coal Society, 2019, 44(2): 482–489. [8] 璩世杰, 刘际飞. 节理角度对预裂爆破成缝效果的影响研究 [J]. 岩土力学, 2015, 36(1): 189–194, 204.QU S J, LIU J F. Numerical analysis of joint angle effect on cracking with presplit blasting [J]. Rock and Soil Mechanics, 2015, 36(1): 189–194, 204. [9] 陈俊桦, 张家生, 李新平. 基于岩石爆破损伤理论的预裂爆破参数研究及应用 [J]. 岩土力学, 2016, 37(5): 1441–1450.CHEN J H, ZHANG J S, LI X P. Study of presplitting blasting parameters and its application based on rock blasting-induced damage theory [J]. Rock and Soil Mechanics, 2016, 37(5): 1441–1450. [10] MA L, LI K M, XIAO S S, et al. Research on effects of blast casting vibration and vibration absorption of presplitting blasting in open cast mine [J]. Shock and Vibration, 2016: 4091732. [11] 饶宇, 夏元友, 胡英国, 等. 预裂缝对爆破振动频谱分布特征的影响 [J]. 振动与冲击, 2017, 36(7): 191–198.RAO Y, XIA Y Y, HU Y G, et al. Influence of pre-splitting crack on spectrum distribution characteristics of blasting vibration [J]. Journal of Vibration and Shock, 2017, 36(7): 191–198. [12] 严鹏, 邹玉君, 卢文波, 等. 基于爆破振动监测的岩石边坡开挖损伤区预测 [J]. 岩石力学与工程学报, 2016, 35(3): 538–548.YAN P, ZOU Y J, LU W B, et al. Predicting the damage zone of rock slopes under blasting excavation based on vibration monitoring [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(3): 538–548. [13] XIAO S S, WANG H S, DONG G W, et al. A preliminary study on the design method for large-diameter deep-hole presplit blasting and its vibration-isolation effect [J]. Shock and Vibration, 2019: 2038578. [14] 徐颖, 孟益平, 程玉生. 装药不耦合系数对爆破裂纹控制的试验研究 [J]. 岩石力学与工程学报, 2002(12): 1843–1847. doi: 10.3321/j.issn:1000-6915.2002.12.020XU Y, MENG Y P, CHEN Y S. Study on control of blast crack by decoupling charge index [J]. Chinese Journal of Rock Mechanics and Engineering, 2002(12): 1843–1847. doi: 10.3321/j.issn:1000-6915.2002.12.020 [15] 何理, 钟冬望. 微差爆破地震波沿高程传播特性的试验研究 [J]. 化工矿物与加工, 2015, 44(5): 36–40.HE L, ZHONG D W. Experimental study on propagation characteristics of millisecond blasting seismic waves along elevation [J]. Industrial Minerals & Processing, 2015, 44(5): 36–40. [16] 司剑峰, 钟冬望, 黄小武. 钻孔爆破孔间最佳延时时间模型试验研究 [J]. 金属矿山, 2015, 468(6): 19–23. doi: 10.3969/j.issn.1001-1250.2015.06.005SI J F, ZHONG D W, HUANG X W. Experimental model of the optimal delay time in drilling blasting [J]. Metal Mine, 2015, 468(6): 19–23. doi: 10.3969/j.issn.1001-1250.2015.06.005 -