Fuzzy Decision Theory of Rock Mass Explodability and Prediction of Explosive Unit Consumption
-
摘要: 为推进岩体可爆性分级后的进一步工程实际应用,利用回归分析方法分析可爆性指标的相关性,最终确定以岩体的抗拉强度、密度、脆性指数、完整性系数4个基本不相关的参数作为评级指标;通过正交试验设计确定各评级指标的敏感性,并确定其权重;运用模糊决策理论对岩体进行可爆性评级,采用评级指标结合岩体基本质量指数,推出地下采矿浅孔爆破的炸药单耗预测公式,进一步推导出可爆性等级对应的炸药单耗预测区间。结果表明:模糊决策方法为岩体可爆性评级提供了一种新思路,同时也证明了权重分配与评级指标选取的正确性。现场爆破试验证明了炸药单耗预测区间的合理性,研究结果可为生产爆破和类似工程实践提供一定的指导。Abstract: In order to promote the further practical application of rock mass explosivity classification in engineering, regression analysis method was used to analyze the correlation of explosivity indicators, and finally determine the four basically unrelated parameters of rock mass tensile strength, density, brittleness index, and integrity coefficient as rating indicators. The sensitivity of each rating index and its weight were determined using orthogonal design experiments. By using the method of fuzzy decision-making, the explosivity of rock mass was rated, and based on the rating index and the basic quality index of rock mass, a formula for predicting the explosive unit consumption of shallow hole blasting in underground mining was derived. The corresponding explosive unit consumption prediction interval for the explosivity level was further derived. This research shows that the fuzzy decision-making method provides a new approach for rock mass explosivity rating, and also proves the correctness of weight allocation and rating index selection. The on-site blasting test has proven the rationality of the predicted range of explosive unit consumption, which can provide certain guiding significance for production blasting and similar engineering practices.
-
岩体可爆性是指岩石或岩体受到炸药爆炸作用时对炸药释放的能量进行吸收或阻挡的能力。岩体可爆性可以用来衡量岩体对爆破作业的适应性和阻碍自身破碎的难易程度。基于岩体可爆性评级结果,可以进行爆破参数设计和优化,以提高爆破效率、降低爆破风险,达到预期的工程目标。这包括确定合适的爆破药量、装药方式、装药布置以及起爆方式等。总之,岩体可爆性评级和基于评级结果的爆破参数设计是提高爆破效率和安全性的基础。
国内外学者对爆破作用下岩体力学性质和岩体可爆性进行了大量研究[1–4]。张紫晗等[5]将岩体脆性指数引入可爆性分级指标中,使用熵权理论对岩体可爆性进行定性和定量评级,该方法可以提供系统的评级结果,帮助爆破参数设计与优化;史涵虚等[6]基于加权聚类分析原理,对魏家峁煤矿区的岩体可爆性进行了研究,主要探讨了不同岩体的可爆性特征,并通过聚类分析的方法进行评级;王文军等[7]运用爆破指数法对镜铁山矿业的黑沟山矿区进行了可爆性分级,爆破指数法是一种常用的爆破参数设计方法,通过评估岩体的可爆性特征来确定适当的爆破参数;Wu等[8]提出了一种基于概念格和粗糙集理论的可爆性分级方法,该方法通过引入概念格和粗糙集理论,提高了分级效率和准确性;Zhou等[9]建立了一种改进的RES-多维云岩团爆破分类模型,该模型结合了不同的评估指标,现场实际应用效果良好;Alipour等[10]采用Mamdani模糊算法,用模糊集来表示爆破稳定性指数,并分析了岩体对爆破碎裂的抵抗能力。此外,还有其他研究方法应用于可爆性分级,如神经网络法[11]、集对分析理论[12]、CRITIC法与Vague集理论[13]、综合赋权云模型[14]等。这些方法通过不同的理论和模型,对岩体可爆性进行多方位、多角度的评估与分级,取得了良好的效果。
以上学者的研究为岩体可爆性分级提供了多种方法和工具,有助于爆破参数的设计与优化,从而提高爆破效率和安全性。然而,在可爆性分级研究中,大多数学者只是对某种岩体进行可爆性评级,评级结果在工程实际中的应用却鲜见报道。基于此,综合考虑指标获取的难易程度,本研究利用相关性分析提高评级指标的独立性,采用正交试验设计确定评级指标的权重,运用模糊决策理论进行岩体可爆性评级,并基于评级结果预测岩体爆破所需炸药单耗的范围。
1. 岩体可爆性分级方法及指标选取
1.1 模糊决策原理
模糊决策方法可以处理不确定和模糊性的问题,适用于岩体可爆性评级。在模糊决策中,可以建立模糊集和隶属函数来表示评级指标与可爆性等级之间的关系。通过模糊推理和模糊逻辑运算,可以实现对岩体可爆性的评级。具体分为5个步骤。
1.1.1 确定研究对象
对于模糊决策而言,首先确定其研究对象:(x1, x2, ··· , xn)为待分类对象的全体。每个样本xn=(x1n, x2n, ··· , xmn)T可用m个指标的特征值表示,则样本集可用m×n阶指标特征值矩阵X表示
X=[x11x12...x1nx21x22...x2n............xm1xm2...xmn]=(xij) (1) 式中:xij为第j个样本对应的第i个指标的特征值。
1.1.2 确定评级指标与评级标准
假设有m个评级指标,共分为c个等级,那么其评级标准表对应的特征值矩阵Y为
Y=[y11y12...y1cy21y22...y2c............ym1ym2...ymc]=(yih) (2) 式中:yih为第i个指标对应的第h个等级所对应的特征值。
1.1.3 评级标准的规格化处理
由于评级指标存在量纲上的差异,为消除其对结果的影响,需要对指标进行规格化处理,使其数据大小处于0~1之间,且量纲均为1。规格化处理可采用
sih=yih−yi1yic−yi1 (3) 式中:sih为第i个指标对应的第h个等级进行规格化处理后的值,且0≤sih≤1;yic、yi1分别为第i项指标对应的第c个等级和第1个等级的特征值。
利用式(3)对评级标准进行规格化处理,得到
S=[s11s12...s1cs21s22...s2c............sm1sm2...smc]=(sih) (4) 1.1.4 确定样本的相对隶属度矩阵
设模糊概念β分为c级,在进行模糊决策的时候应先消除指标量纲的影响,基本思路如下。对模糊概念β而言,特征值越大,级别越高,即极大型指标,可以认为:当xij大于或等于c级特征值yic时,xij对于β的相对隶属度为1;反之,当xij小于或等于1级特征值yi1时,其相对隶属度为零;若xij介于yic与yi1之间时,则用线性插值法确定其隶属度,隶属度介于0~1之间。相反,若对模糊概念β而言,特征值越小,级别越高,即极小型指标,则认为:当xij小于或等于c级特征值yic时,xij对于β的相对隶属度为1;当xij大于或等于1级特征值yi1时,其相对隶属度为零;xij介于两者之间时,则用线性插值确定。
综上,样本对模糊概念β的相对隶属度rij的公式为
rij={1xij⩾ (5) 式中:rij为样本第i个指标的第j个特征值对β的相对隶属度;左边公式对应极大型指标,右侧公式对应极小型指标。
利用式(5)对样本进行处理,得到其相对隶属度矩阵R
\boldsymbol R = \left[ {\begin{array}{*{20}{c}} {{r_{11}}}&{{r_{12}}}&{...}&{{r_{1c}}} \\ {{r_{21}}}&{{r_{22}}}&{...}&{{r_{2c}}} \\ {...}&{...}&{...}&{...} \\ {{r_{m1}}}&{{r_{m2}}}&{...}&{{r_{mc}}} \end{array}} \right] = ({r_{ij}}) (6) 1.1.5 确定模糊决策矩阵
样本与评级标准之间的贴近程度用加权距离来衡量,在已知模糊标准矩阵以及各指标的权重矩阵后,模糊决策矩阵中的元素可用下式得出
{\mu _{hj}} = \dfrac{1}{{\displaystyle\sum\limits_{k = 1}^c {{{\left[ {\dfrac{{{{\displaystyle\sum\limits_{i = 1}^m {\left[ {{\omega _i}({r_{ij}} - {s_{ih}})} \right]^p} }}}}{{{{\displaystyle\sum\limits_{i = 1}^m {\left[ {{\omega _i}({r_{ij}} - {s_{ik}})} \right]^p} }}}}} \right]}^{{2}/{p}}}} }} (7) 式中:μhj为样本j从属于类别h的隶属度;ωi为第i项指标的权重,
\sum\limits_{i = 1}^m {{\omega _i}} = 1 ;p为系数。1.2 岩体可爆性分级指标的选取
岩体可爆性的影响因素复杂多样,针对不同的研究内容,选择的评价指标也有所不同。岩石的单轴抗压强度、单轴抗拉强度、动载冲击强度、抗剪强度、密度、脆性指数、岩体完整性系数等是目前众多学者采用相对较多的评级指标。本研究拟通过相关性分析找出相关性较低的评级指标用于可爆性评级,以避免冗余信息,提高评级指标的独立性。
1.2.1 评级指标的相关性分析
根据岩体爆破理论,岩体破坏是爆破冲击波与爆生气体共同作用的结果。当爆破冲击波和爆生气体对岩体的作用力超过岩体本身的抗剪强度时,岩体就会发生破坏。因此,岩石的抗剪强度可以作为评估岩体可爆性的指标。
脆性[15]既是一种变形特性,也是一种材料特性。脆性可以用脆性指数度量,脆性指数等于抗压强度与抗拉强度的比值。爆破是一种对岩体施加动态冲击的过程,脆性程度代表岩体对外部冲击的抵抗能力,所以将脆性指数作为评估岩体可爆性的指标是合理的。当岩石的脆性较高时,岩石受到爆破冲击时更容易发生破碎。
璩世杰等[16]通过对14种岩石的多种物理力学参数进行相关性分析,得到了岩石的静载抗拉强度、密度、完整性系数可以较好地评价岩体的可爆性的结论。在选取岩石单轴抗拉强度、密度、完整性系数作为可爆性评级指标的基础上,为确定能否将岩石的抗剪强度和脆性指数也引入评级指标,选取文献[17–18]中岩石的抗拉强度和抗剪强度以及文献[5, 15, 19–20]中岩石的物理力学参数进行相关性分析。岩石的抗拉强度(σt)和抗剪强度(τ)见表1。岩石的物理力学参数见表2,其中:ρ为岩石密度,η为岩体的完整性系数,B为岩石脆性指数。共分析测试39块岩石样品的物理力学参数,表2仅列出12块岩石样品的数据。
表 1 岩石的抗拉强度和抗剪强度Table 1. Tensile and shear strength of rocksRock category σt/MPa τ/MPa Rock category σt/MPa τ/MPa Magnetite ore 12 43 Migmatite 8.5 38 Chlorite schist 2.6 12 Quartz schist 7 16 Magnetite rich ore 11.5 40 Mixed metamorphic rock 4.81 21.66 Amphibolite 7.07 22 Mixed granite 12.93 47.85 k shell 8 32 Fault breccia 8.53 33.02 表 2 岩石的物理力学参数Table 2. Physical and mechanical parameters of rockRock No. ρ/(g·cm−3) σt/MPa η B 1 2.618 9.30 0.767 12.74 2 2.674 14.25 0.760 5.46 3 3.810 5.85 0.510 16.20 4 2.710 7.05 0.410 22.70 5 2.646 14.62 0.780 4.40 6 2.588 7.61 0.830 5.27 7 2.919 15.33 0.793 3.94 8 2.631 9.71 0.763 16.55 9 3.418 16.30 0.230 11.87 10 3.029 21.04 0.197 7.27 11 2.725 14.98 0.480 8.01 12 3.479 18.26 0.779 10.55 分别对表1与表2中的数据进行相关性分析,结果见表3,其中:β0、β1为回归方程系数,R为相关度。从结果看:岩石的抗拉强度与抗剪强度的相关度为0.8138,相关度较高,不可将两者同时用作可爆性评级指标;岩石脆性指数与岩石密度、单轴抗拉强度、完整性系数的相关度分别为0.0084、0.2386、0.0273,基本线性无关,即可以在岩石密度、单轴抗拉强度、完整性系数作为评级指标的前提下,将脆性指数作为评级指标的一部分。
表 3 相关性分析结果Table 3. Correlation analysis resultsResults β0 β1 R τ=0.3658+3.5369σt 0.3658 3.5369 0.8138 B=5.2149+2.0844ρ 5.2149 2.0844 0.0084 B=19.9728−0.7139σt 19.9728 −0.7139 0.2386 B=14.3477−5.4339η 14.3477 −5.4339 0.0273 相对于抗拉强度,抗剪强度(内聚力与内摩擦角)的获取途径复杂,为此,本研究最终确定以岩石密度、单轴抗拉强度、完整性系数、岩石脆性指数作为可爆性的评级指标。
1.3 权重系数的确定
在岩体可爆性评级过程中,各分级指标权重的确定是至关重要的一步,它体现岩体的物理力学参数对最终评级结果的贡献度。本研究选择应用每个因子3个水平的正交试验设计确定4个评级指标的敏感性,进而确定其权重。通过正交设计,可以高效地确定权重,减少主观性和随机性对权重的影响。正交设计中各因素水平见表4,L9(34)型正交试验设计见表5,正交试验设计结果见表6。利用极差分析方法对结果进行分析,将各列水平数相同的结果相加,记为K,3个水平的计算结果分别为K1、K2、K3。表6中的等级基于文献[12]中提到的方法得出。
表 4 正交设计各因素水平Table 4. Levels of various factors in orthogonal designLevel ρ/(g·cm−3) σt/MPa B η 1 2.618 8.61 10.77 0.360 2 2.975 13.38 7.79 0.513 3 3.270 20.14 5.46 0.783 表 5 正交设计试验L9(34)Table 5. Orthogonal design test L9(34)Test No. Level ρ σt B η 1 1 1 1 1 2 1 2 2 2 3 1 3 3 3 4 2 1 2 3 5 2 2 3 1 6 2 3 1 2 7 3 1 3 2 8 3 2 1 3 9 3 3 2 1 表 6 正交试验设计结果Table 6. Orthogonal design test resultsTest No. ρ/(g·cm−3) σt/MPa B η Level[12] 1 2.618 8.61 10.77 0.360 2 2 2.618 13.38 7.79 0.513 4 3 2.618 20.14 5.46 0.783 5 4 2.975 8.61 7.79 0.783 7 5 2.975 13.38 5.46 0.360 3 6 2.975 20.14 10.77 0.513 4 7 3.270 8.61 5.46 0.513 6 8 3.270 13.38 10.77 0.783 6 9 3.270 20.14 7.79 0.360 5 K1 11 15 12 10 K2 14 13 16 14 K3 17 14 14 18 Range 2 0.667 1.333 2.667 Sensitiveness 2 4 3 1 通过表6可以看出,对于评价岩体可爆性的4个指标,按照敏感性由大到小依次为:完整性系数、岩体密度、岩石脆性指数、岩石单轴抗拉强度,因此在确定评级指标的权重时应予以考虑。根据极差分析结果,完整性系数、岩体密度、岩石脆性指数、岩石单轴抗拉强度的极差比值为4∶3∶2∶1,由此确定其权重为(0.4, 0.3, 0.2, 0.1)。文献[6,18]指出,岩石的单轴抗拉强度对可爆性评级的影响较大,最终将权重调整为(0.35, 0.29, 0.21, 0.15)。
2. 岩体可爆性等级评定
选取文献[5]中的12种岩体试样作为样本,样本集指标特征值矩阵X为
{ \boldsymbol {X}} ={\left[ {\begin{array}{*{20}{c}} {2.618}&{2.674}&{2.604}&{3.070}&{2.875}&{2.813}&{2.962}&{3.368}&{3.388}&{3.445}&{3.325}&{3.479} \\ {9.30}&{14.25}&{11.12}&{17.38}&{9.41}&{6.61}&{9.00}&{15.83}&{25.66}&{9.34}&{15.27}&{18.26} \\ {12.74}&{5.46}&{11.64}&{9.02}&{7.79}&{3.13}&{8.69}&{9.11}&{6.41}&{13.90}&{10.77}&{10.55} \\ {0.767}&{0.760}&{0.796}&{0.783}&{0.783}&{0.767}&{0.783}&{0.783}&{0.774}&{0.783}&{0.763}&{0.779} \end{array}} \right]} = {({x_{ij}})_{4 \times 12}} (8) 基于以上分析,以岩石单轴抗拉强度、岩石密度、岩石脆性指数、岩体完整性系数为评级指标,以文献[5, 21]为参考,建立可爆性分级标准,如表7所示,共分为7个等级,即最易爆Ⅰ、易爆Ⅱ、较易爆Ⅲ、中等Ⅳ、较难爆Ⅴ、难爆Ⅵ、极难爆Ⅶ。
表 7 可爆性分级标准Table 7. Explosivity classification standardsRock mass explosivity level ρ/(g·cm−3) σt/MPa B η Level Ⅰ 2.50 6.6 17.5 0.0494 Easiest Ⅱ 2.60 10.0 15.5 0.2555 Easy Ⅲ 2.75 13.0 13.3 0.3654 Easier Ⅳ 2.90 17.0 11.8 0.5122 Medium Ⅴ 3.16 20.0 9.5 0.6021 More difficult Ⅵ 3.30 23.0 6.0 0.7122 Difficult Ⅶ 3.45 26.0 3.5 0.8232 Most difficult 取表7中的数据建立可爆性评级标准对应的特征值矩阵Y
\boldsymbol Y = \left[ {\begin{array}{*{20}{c}} {2.50}&{2.60}&{2.75}&{2.90}&{3.16}&{3.30}&{3.45} \\ {6.6}&{10.0}&{13.0}&{17.0}&{20.0}&{23.0}&{26.0} \\ {17.5}&{15.5}&{13.3}&{11.8}&{9.5}&{6.0}&{3.5} \\ {0.049\,4}&{0.255\,5}&{0.365\,4}&{0.512\,2}&{0.602\,1}&{0.712\,2}&{0.823\,2} \end{array}} \right] = {({y_{ij}})_{4 \times 7}} (9) 采用式(3)对Y进行规格化处理,得到消除量纲影响后的标准特征值矩阵S
\boldsymbol S = \left[ {\begin{array}{*{20}{c}} 0&{0.105\,3}&{0.263\,2}&{0.421\,1}&{0.694\,7}&{0.842\,1}&1 \\ 0&{0.175\,3}&{0.329\,9}&{0.536\,1}&{0.690\,7}&{0.845\,4}&1 \\ 0&{0.142\,9}&{0.300\,0}&{0.407\,1}&{0.571\,4}&{0.821\,4}&1 \\ 0&{0.266\,3}&{0.408\,4}&{0.598\,1}&{0.714\,3}&{0.856\,6}&1 \end{array}} \right] = {({s_{ih}})_{4 \times 7}} (10) 利用式(5)对X进行规格化处理,确定样本的相对隶属度矩阵R
\boldsymbol R = {\left[ {\begin{array}{*{20}{c}} {0.124\,2}& {0.183\,2}& {0.109\,5}& {0.600\,0}& {0.394\,7}& {0.329\,5}& {0.486\,3}& {0.913\,7}& {0.934\,7}& {0.994\,7}& {0.868\,4}& 1.000\,0 \\ {0.139\,2}& {0.394\,3}& {0.233\,0}& {0.555\,7}& {0.144\,8}& {0.000\,5}& {0.123\,7}& {0.475\,8}& {0.9825}& {0.141\,2}& {0.446\,9}& {0.601\,0} \\ {0.340\,0}& {0.860\,0}& {0.418\,6}& {0.605\,7}& {0.693\,6}& {1.000\,0}& {0.629\,3}& {0.599\,3}& {0.792\,1}& {0.257\,1}& {0.480\,7}& {0.496\,4} \\ {0.927\,4}& {0.918\,3}& {0.964\,8}& {0.948\,0}& {0.948\,0}& {0.927\,4}& {0.948\,0}& {0.948\,0}& {0.936\,4}& {0.948\,0}& {0.922\,2}& {0.942\,9} \end{array}} \right] } = {({r_{ij}})_{4 \times 12}} (11) 通过式(7)进行计算,取p=2,计算得到样本的模糊决策矩阵U
\boldsymbol U = {\left[ {\begin{array}{*{20}{c}} {0.061\,5}& {0.047\,3}& {0.059\,2}& {0.017\,2}& {0.037\,0}& {0.047\,8}& {0.032\,8}& {0.013\,9}& {0.004\,4}& {0.034\,9}& {0.016\,9}& {0.016\,1} \\ {0.124\,7}& {0.088\,1}& {0.116\,2}& {0.031\,2}& {0.069\,2}& {0.083\,4}& {0.060\,9}& {0.022\,7}& {0.007\,0}& {0.055\,3}& {0.028\,5}& {0.025\,7} \\ {0.194\,2}& {0.145\,2}& {0.180\,2}& {0.055\,7}& {0.119\,7}& {0.132\,0}& {0.106\,3}& {0.037\,0}& {0.011\,0}& {0.084\,2}& {0.047\,6}& {0.040\,7} \\ {0.282\,2}& {0.250\,8}& {0.274\,9}& {0.144\,6}& {0.242\,1}& {0.212\,9}& {0.227\,0}& {0.075\,4}& {0.021\,6}& {0.145\,5}& {0.102\,6}& {0.079\,2} \\ {0.163\,6}& {0.205\,3}& {0.173\,9}& {0.355\,3}& {0.247\,7}& {0.210\,9}& {0.273\,5}& {0.237\,1}& {0.038\,2}& {0.270\,1}& {0.324\,1}& {0.233\,2} \\ {0.105\,5}& {0.161\,6}& {0.118\,1}& {0.284\,1}& {0.180\,8}& {0.188\,6}& {0.193\,6}& {0.417\,2}& {0.519\,2}& {0.243\,6}& {0.333\,0}& {0.379\,6} \\ {0.068\,3}& {0.101\,8}& {0.077\,3}& {0.112\,1}& {0.103\,5}& {0.124\,4}& {0.105\,7}& {0.196\,7}& {0.368\,6}& {0.166\,4}& {0.147\,3}& {0.225\,5} \end{array}} \right]} = {({u_{hj}})_{7 \times 12}} (12) 根据最大隶属度原则得到可爆性等级结果,见表8。参考文献[5, 21],采用熵权理论法和属性识别法也得到可爆性等级,见表8。
表 8 可爆性分级结果Table 8. Explosion classification resultsRock No. Proposed method Attribute recognition method Entropy weight theory method 1 Ⅳ Ⅲ Ⅲ 2 Ⅳ Ⅲ Ⅲ 3 Ⅳ Ⅴ Ⅴ 4 Ⅴ Ⅴ Ⅵ 5 Ⅴ Ⅴ Ⅴ 6 Ⅳ Ⅳ Ⅳ 7 Ⅴ Ⅳ Ⅴ 8 Ⅵ Ⅵ Ⅵ 9 Ⅵ Ⅵ Ⅵ 10 Ⅴ Ⅵ Ⅶ 11 Ⅵ Ⅵ Ⅵ 12 Ⅵ Ⅵ Ⅵ 从结果来看,本研究得到的评级结果与其他两种方法得到的结果具有较好的一致性,说明模糊决策方法具有一定的可行性,同时也证明采用正交试验设计确定各评级指标的权重具有合理性。然而,从实际应用的角度来看,模糊决策方法更适合处理数据,它基于模糊集理论,可以处理模糊的、有限的、不完全的和不精确的数据,使决策过程更加灵活,可以适应各种实际情况。此外,模糊决策方法还可以方便地植入数学软件中进行计算和分析,有许多基于模糊决策方法的软件工具和数据库可供使用,它们提供了专门的算法和函数,能够快速进行模糊推理和评级计算。这些软件工具使实施模糊决策方法变得简单,大大减少了计算工作量和人工操作的复杂性。
岩体的可爆性评级是多个因素共同作用的结果,任何一种因素发生变化都会对评级结果产生影响。在本研究的评级方法中,岩石1和岩石2的评级结果比其他两种方法高一级,而岩石3和岩石10的评级结果则比其他两种方法低,这是因为3种评级方法的权重分配方案不同。对于岩石10,其密度、抗拉强度、完整性系数、脆性指数分别为3.445 g/cm3、9.34 MPa、0.783、13.90,与其他岩石相比,其脆性指数较大,且脆性指数是极小型指标,其值越大代表岩体可爆性等级越低,岩体越容易爆破。相对其余两种方法,本研究赋予脆性指数的权重略大,因此岩石10的评级结果比其他两种方法低。
3. 可爆性等级对应炸药单耗区间的预测
基于可爆性评级结果预测岩体爆破所需炸药单耗的范围是评价岩体可爆性的一种补充应用,可为工程实际应用提供更详细的指导信息。本研究利用多元回归分析方法,对可爆性等级与炸药单耗之间的关系建立数学模型,从而实现对炸药单耗范围预测。
3.1 炸药单耗预测公式推导
对于非煤矿山地下采矿浅孔爆破,爆破每立方米岩体所需的炸药量q通常要根据不同岩石的坚固性系数f进行调整,具体如表9[22]所示。
参考《工程岩体分级标准》(GB/T 50218—2014)[23],基于岩体基本质量指标BQ,将岩体分为5级,将其与表9结合得到表10。
表 10BQ 与炸药单耗Table 10.BQ and explosive unit consumptionBQ q/(kg·m−3) BQ q/(kg·m−3) ≤250 0.25–0.65 451–550 1.60–2.60 251–350 0.65–1.00 >551 >2.80 351–450 1.00–1.60 利用MATLAB将表10中的数据进行拟合,得到炸药单耗q与BQ之间的关系,见表11,qll为炸药单耗的下限,qul为炸药单耗的上限。
表 11 炸药单耗与BQ 拟合结果Table 11. Fitting results of explosive unit consumption andBQ Result β0 β1 R qll=–0.6180+0.0055BQ –0.6180 0.0055 0.9629 qul=–0.7425+0.0075BQ −0.7425 0.0075 0.9425 岩体基本质量指标BQ计算公式为
BQ=100+3R_{\mathrm{c}}+250\eta (13) 式中:Rc为岩石的单轴饱和抗压强度。当Rc>90η+30时,将η和Rc=90η+30代入式(13)求BQ;当η>0.04Rc+0.4时,将η=0.04Rc+0.4和Rc代入式(13)求BQ的值。
基于此,得到39种岩石的BQ,如表12所示,与表2一致,这里仅列其中12种岩石的BQ。
表 12 岩石的BQ Table 12. RockBQ Rock No. BQ Rock No. BQ 1 588 7 479 2 523 8 588 3 455 9 309 4 403 10 597 5 489 11 440 6 428 12 595 根据工程实际,对岩体的BQ进行修正,利用表2中线性相关度较弱的数据,将岩体密度、抗拉强度、脆性指数、完整性系数与BQ进行拟合,得到的拟合公式为:BQ=57.4-5.5ρ+7.41σt+460.2η+5.22B,相关度为0.913。将该式与表11中炸药单耗拟合公式结合,得到基于评级指标的炸药单耗预测公式。
完整性系数可以由岩体与岩石波速计算得到,即知道岩体与岩石波速、岩体的密度、抗拉强度、抗压强度,就可通过式(14)和式(15)计算出炸药单耗的上限和下限。
q_{\mathrm{ul}}=-0.311\,9-0.041\,2\rho+0.055\,6\sigma\mathrm{_t}+3.451\,4\eta+0.039\,1B (14) q\mathrm{_{ll}}=-0.302\,2-0.030\,2\rho+0.040\,7\sigma\mathrm{_t}+2.531\eta+0.028\,7B (15) 3.2 可爆性等级对应的炸药单耗
以39种岩石为样本,采用本研究的方法对其进行评级,利用式(14)和式(15),计算炸药单耗的范围,结果见表13。
表 13 评级与炸药单耗计算结果Table 13. Rating and explosive unit consumption calculation resultsLevel Number q/(kg·m−3) Level Number q/(kg·m−3) Ⅰ 0 <0.77 Ⅴ 12 1.91–2.71 Ⅱ 3 0.77–1.22 Ⅵ 5 2.48–3.49 Ⅲ 8 1.17–1.69 Ⅶ 0 >3.49 Ⅳ 11 1.58–2.26 由表13可以看出,39种岩石中可爆性等级为Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ级的分别有3、8、11、12、5种。将评级结果与式(14)和式(15)计算得到的炸药单耗区间进行对应,即可得到各可爆性等级对应的炸药单耗预测范围。
3.3 炸药单耗预测的有效性验证
为验证炸药单耗预测范围的有效性,在云南某锡矿开展岩石的物理力学实验,并进行现场爆破试验。试验地点位于三坑1480中段,爆破试验的对象岩体为大理岩和氧化矿。
3.3.1 岩石的室内实验
采用微机电液伺服试验机HYE-2000对矿岩试样进行单轴抗压和单轴抗拉实验,采用ZT801岩体参数测定仪测量矿岩试样的波速,通过称量法测量矿岩试样的密度,得到的岩石物理力学数据见表14。
表 14 岩矿试样的物理参数Table 14. Physical parameters of rock specimenRock type ρ/(g·cm−3) σt/MPa B η Marble 2.72 4.66 12.9 0.566 Oxidized ore 2.47 4.86 15.1 0.156 根据表14所列数据,采用本研究的方法对岩矿进行可爆性评级:大理岩的可爆性等级为Ⅳ级,氧化矿的可爆性等级为Ⅰ级。由表13即可得到其对应的炸药单耗范围,大理岩的炸药单耗范围为1.58~2.26 kg/m3,氧化矿的炸药单耗小于0.77 kg/m3。
3.3.2 现场爆破试验
分别进行1次大断面大理岩巷道掘进的爆破试验以及1次氧化矿采场爆破试验。试验采用大直径空孔直眼掏槽,装药孔与空孔的距离为200 mm,为空孔直径D的2.0倍,起爆网路采用微差起爆技术。大理岩掏槽孔布置见图1,爆破效果见图2。氧化矿装药作业以及连线作业分别见图3和图4。
爆破后,对出矿、进尺、孔深等相关数据进行统计,发现无特殊大块产生,大理岩爆破的平均炸药单耗为2.04 kg/m3,氧化矿爆破的平均炸药单耗为0.71 kg/m3,处于预测的炸药单耗范围内,这在一定程度上证明了预测范围的合理性,为相似工况的炸药单耗设计提供了一定的参考。
4. 结 论
(1) 通过相关性分析得到了抗拉强度、岩石密度、脆性指数和岩体完整性系数基本不相关的结论,4个量能提供独立的信息,可同时作为可爆性评级指标;
(2) 通过设计一系列正交试验,确定了评级指标对结果的影响,在此基础上确定了各指标权重;
(3) 模糊决策理论是一套科学合理的决策理论,为岩体可爆性评级提供了一种新思路,通过评级结果的对比,证明了该方法的可靠性;
(4) 结合可爆性评级指标和岩体基本质量指标,推导出了炸药单耗的预测公式,进一步给出了各可爆性等级对应的炸药单耗预测区间,现场爆破试验证明了其合理性。
-
表 1 岩石的抗拉强度和抗剪强度
Table 1. Tensile and shear strength of rocks
Rock category σt/MPa τ/MPa Rock category σt/MPa τ/MPa Magnetite ore 12 43 Migmatite 8.5 38 Chlorite schist 2.6 12 Quartz schist 7 16 Magnetite rich ore 11.5 40 Mixed metamorphic rock 4.81 21.66 Amphibolite 7.07 22 Mixed granite 12.93 47.85 k shell 8 32 Fault breccia 8.53 33.02 表 2 岩石的物理力学参数
Table 2. Physical and mechanical parameters of rock
Rock No. ρ/(g·cm−3) σt/MPa η B 1 2.618 9.30 0.767 12.74 2 2.674 14.25 0.760 5.46 3 3.810 5.85 0.510 16.20 4 2.710 7.05 0.410 22.70 5 2.646 14.62 0.780 4.40 6 2.588 7.61 0.830 5.27 7 2.919 15.33 0.793 3.94 8 2.631 9.71 0.763 16.55 9 3.418 16.30 0.230 11.87 10 3.029 21.04 0.197 7.27 11 2.725 14.98 0.480 8.01 12 3.479 18.26 0.779 10.55 表 3 相关性分析结果
Table 3. Correlation analysis results
Results β0 β1 R τ=0.3658+3.5369σt 0.3658 3.5369 0.8138 B=5.2149+2.0844ρ 5.2149 2.0844 0.0084 B=19.9728−0.7139σt 19.9728 −0.7139 0.2386 B=14.3477−5.4339η 14.3477 −5.4339 0.0273 表 4 正交设计各因素水平
Table 4. Levels of various factors in orthogonal design
Level ρ/(g·cm−3) σt/MPa B η 1 2.618 8.61 10.77 0.360 2 2.975 13.38 7.79 0.513 3 3.270 20.14 5.46 0.783 表 5 正交设计试验L9(34)
Table 5. Orthogonal design test L9(34)
Test No. Level ρ σt B η 1 1 1 1 1 2 1 2 2 2 3 1 3 3 3 4 2 1 2 3 5 2 2 3 1 6 2 3 1 2 7 3 1 3 2 8 3 2 1 3 9 3 3 2 1 表 6 正交试验设计结果
Table 6. Orthogonal design test results
Test No. ρ/(g·cm−3) σt/MPa B η Level[12] 1 2.618 8.61 10.77 0.360 2 2 2.618 13.38 7.79 0.513 4 3 2.618 20.14 5.46 0.783 5 4 2.975 8.61 7.79 0.783 7 5 2.975 13.38 5.46 0.360 3 6 2.975 20.14 10.77 0.513 4 7 3.270 8.61 5.46 0.513 6 8 3.270 13.38 10.77 0.783 6 9 3.270 20.14 7.79 0.360 5 K1 11 15 12 10 K2 14 13 16 14 K3 17 14 14 18 Range 2 0.667 1.333 2.667 Sensitiveness 2 4 3 1 表 7 可爆性分级标准
Table 7. Explosivity classification standards
Rock mass explosivity level ρ/(g·cm−3) σt/MPa B η Level Ⅰ 2.50 6.6 17.5 0.0494 Easiest Ⅱ 2.60 10.0 15.5 0.2555 Easy Ⅲ 2.75 13.0 13.3 0.3654 Easier Ⅳ 2.90 17.0 11.8 0.5122 Medium Ⅴ 3.16 20.0 9.5 0.6021 More difficult Ⅵ 3.30 23.0 6.0 0.7122 Difficult Ⅶ 3.45 26.0 3.5 0.8232 Most difficult 表 8 可爆性分级结果
Table 8. Explosion classification results
Rock No. Proposed method Attribute recognition method Entropy weight theory method 1 Ⅳ Ⅲ Ⅲ 2 Ⅳ Ⅲ Ⅲ 3 Ⅳ Ⅴ Ⅴ 4 Ⅴ Ⅴ Ⅵ 5 Ⅴ Ⅴ Ⅴ 6 Ⅳ Ⅳ Ⅳ 7 Ⅴ Ⅳ Ⅴ 8 Ⅵ Ⅵ Ⅵ 9 Ⅵ Ⅵ Ⅵ 10 Ⅴ Ⅵ Ⅶ 11 Ⅵ Ⅵ Ⅵ 12 Ⅵ Ⅵ Ⅵ f q/(kg·m−3) f q/(kg·m−3) 0–4 0.25–0.65 10–15 1.60–2.60 5–8 0.65–1.00 >15 >2.80 8–10 1.00–1.60 表 10
BQ 与炸药单耗Table 10.
BQ and explosive unit consumptionBQ q/(kg·m−3) BQ q/(kg·m−3) ≤250 0.25–0.65 451–550 1.60–2.60 251–350 0.65–1.00 >551 >2.80 351–450 1.00–1.60 表 11 炸药单耗与
BQ 拟合结果Table 11. Fitting results of explosive unit consumption and
BQ Result β0 β1 R qll=–0.6180+0.0055BQ –0.6180 0.0055 0.9629 qul=–0.7425+0.0075BQ −0.7425 0.0075 0.9425 表 12 岩石的
BQ Table 12. Rock
BQ Rock No. BQ Rock No. BQ 1 588 7 479 2 523 8 588 3 455 9 309 4 403 10 597 5 489 11 440 6 428 12 595 表 13 评级与炸药单耗计算结果
Table 13. Rating and explosive unit consumption calculation results
Level Number q/(kg·m−3) Level Number q/(kg·m−3) Ⅰ 0 <0.77 Ⅴ 12 1.91–2.71 Ⅱ 3 0.77–1.22 Ⅵ 5 2.48–3.49 Ⅲ 8 1.17–1.69 Ⅶ 0 >3.49 Ⅳ 11 1.58–2.26 表 14 岩矿试样的物理参数
Table 14. Physical parameters of rock specimen
Rock type ρ/(g·cm−3) σt/MPa B η Marble 2.72 4.66 12.9 0.566 Oxidized ore 2.47 4.86 15.1 0.156 -
[1] 坦加耶夫 И А. 岩石的可钻性和可爆性 [M]. 王伟德, 周树良,译. 北京: 冶金工业出版社, 1987: 55–113.TANGAYEV И А. Drillability and explosive of rock [M]. Translated by WANG W D, ZHOU S L. Beijing: Metallurgical Industry Press, 1987: 55–113. [2] 张强. 岩体爆破性分级研究进展 [J]. 工程爆破, 1998, (3): 75–79,35.ZHANG Q. Advance in study on rock blastability classification [J]. Engineering Blasting, 1998, (3): 75–79,35. [3] XIAO S S, LI K M, DING X H, et al. Rock mass blastability classification using fuzzy pattern recognition and the combination weight method [J]. Mathematical Problems in Engineering, 2015, 12: 724619. [4] OZGUR Y. Rock factor prediction in the Kuz-Ram model and burden estimation by mean fragment size [J]. Geomechanics for Energy and the Environment, 2023, 33: 100415. doi: 10.1016/j.gete.2022.100415 [5] 张紫晗, 胡光球, 郑建礼, 等. 基于脆性指数和熵权理论的岩体可爆性分级方法 [J]. 工程爆破, 2021, 27(5): 65–71.ZHANG Z H, HU G Q, ZHENG J L, et al. Rock blastability classification method based on brittleness index and entropy weight theory [J]. Engineering Blasting, 2021, 27(5): 65–71. [6] 史涵虚, 周传波, 蒋楠, 等. 魏家峁煤矿区深孔台阶剥离爆破施工岩体可爆性分析 [J]. 爆破, 2023, 40(1): 37–44.SHI H X, ZHOU C B, JIANG N, et al. Analysis on rock mass explosibility during deep-hole bench stripping by blasting in Weijiamao coal mine area [J]. Blasting, 2023, 40(1): 37–44. [7] 王文军, 李润然, 纪旭波, 等. 黑沟矿区露天开采矿岩可爆性分级研究 [J]. 金属矿山, 2021(4): 64–69.WANG W J, LI R R, JI X B, et al. Blastability classification of open pit mining in heigou mining area [J]. Metal Mine, 2021(4): 64–69. [8] WU S L, YANG S, WANG Q Y. Classification of open pit iron mine rock mass blastability based on concept lattice and rough set [J]. Geotechnical and Geological Engineering: An International Journal, 2020, 38(1): 449–458. doi: 10.1007/s10706-019-01036-3 [9] ZHOU T, HU H J, KUANG Y, et al. Improved rock engineering system (RES)-multidimensional cloud evaluation model and its application to the rock mass blastability [J]. IEEE Access, 2019, 7: 100305. doi: 10.1109/ACCESS.2019.2930629 [10] ALIPOUR A, MOKHARIAN M, CHEHREGHANI S. An application of fuzzy sets to the blastability index (BI) used in rock engineering [J]. Periodica Polytechnica Civil Engineering, 2018, 62(3): 580–589. [11] 冯夏庭. 岩石可爆性神经网络研究 [J]. 爆炸与冲击, 1994, 14(4): 298–306.FENG X T. A study on neural network on rock blastability [J]. Explosion and Shock Waves, 1994, 14 (4): 298–306. [12] 马红贝, 赵国彦, 路凡. 岩体可爆性分级评价的集对分析模型 [J]. 爆破, 2016, 33(2): 28–31, 38.MA H B, ZHAO G Y, LU F. Set pair analysis model for rock blastability classification estimation [J]. Blasting, 2016, 33(2): 28–31, 38. [13] 白玉奇, 李振阳, 李传增, 等. 基于CRITIC-Vague模型的岩体可爆性评价 [J]. 金属矿山, 2020(12): 27–32.BAI Y Q, LI Z Y, LI C Z, et al. Assessment for rock mass blastability based on CRITIC-Vague model [J]. Metal Mine, 2020(12): 27–32. [14] 丁小华, 原文杰, 解祯, 等. 基于综合赋权云模型的露天矿岩体可爆性分级识别 [J]. 煤炭科学技术, 2019, 47(10): 96–101.DING X H, YUAN W J, XIE Z, et al. Classification and identification of rock blastability in open-pit mine based on comprehensive weighted cloud model [J]. Coal Science and Technology, 2019, 47(10): 96–101. [15] 刘恩龙, 沈珠江. 岩土材料的脆性研究 [J]. 岩石力学与工程学报, 2005, 24(19): 51–55.LIU E L, SHEN Z J. Study on brittleness of geomaterials [J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(19): 51–55. [16] 璩世杰, 辛明印, 毛市龙, 等. 岩体可爆性指标的相关性分析 [J]. 岩石力学与工程学报, 2005, 24(3): 468–473.QU S J, XIN M Y, MAO S L, et al. Correlation analyses of blastability indexes for rock mass [J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(3): 468–473. [17] 潘勇, 杨天祥. 岩体可爆性分级的神经网络判别方法 [J]. 现代矿业, 2013, 29(1): 21–23.PAN Y, YANG T X. Neural network discriminant method of rock blasting classification [J]. Modern Mining, 2013, 29(1): 21–23. [18] 璩世杰, 齐宝军, 许文耀, 等. 水厂铁矿岩体可爆性分级方法与应用 [J]. 金属矿山, 2009, 39(11): 25–27, 178.QU S J, QI B J, XU W Y, et al. Rock mass blastability classification and its application in shuichang open pit mine [J]. Metal Mine, 2009, 39(11): 25–27, 178. [19] 张德明, 王新民, 郑晶晶, 等. 基于模糊综合评判的矿岩体可爆性分级 [J]. 爆破, 2010, 27(4): 43–47.ZHANG D M, WANG X M, ZHENG J J, et al. Blastability classification of rock and mine based on fuzzy comprehensive evaluation [J]. Blasting, 2010, 27(4): 43–47. [20] 高成宇. 巴鲁巴铜矿基于加权聚类分析的岩体可爆性分级 [J]. 铜业工程, 2014(2): 28–31.GAO C Y. Rock-mass blastability classification of baluba copper mine based on weighted clustering analysis [J]. Copper Engineering, 2014(2): 28–31. [21] 李蓉, 宋娟, 何永延. 基于属性识别理论的岩体可爆性分级方法 [J]. 金属矿山, 2008, 38(5): 32–34,48.LI R, SONG J, HE Y Y. Classification for rockmass blastability based on attribute recognition theory [J]. Metal Mine, 2008, 38(5): 32–34,48. [22] 汪旭光. 爆炸设计与施工 [M]. 第2版. 北京: 冶金工业出版社, 2011: 303–305.WANG X G. Explosion design and construction [M]. 2nd ed. Beijing: Metallurgical Industry Press, 2011: 303–305. [23] 长江水利委员会长江科学院. 工程岩体分级标准: GB/T 50218—2014 [S]. 北京: 中国计划出版社, 2015. -