受限自由面爆破装药量计算公式的优化研究.pdf

第3 4 卷第1 期 2 0 1 7 年3 月 爆破 B L A S T I N G V 0 1 ．3 4N o ．1 M a r ．2 0 1 7 d o i 1 0 ．3 9 6 3 ／j ．i s s n ．1 0 0 1 - 4 8 7 X ．2 0 1 7 ．0 1 ．0 0 7 受限自由面爆破装药量计算公式的优化研究术 李启月，张成君，吴正宇，陈英，韦佳瑞 中南大学资源与安全工程学院，长沙4 1 0 0 8 3 摘要为提高巷道掘进爆破速度及精确控制巷道轮廓。针对现有装药量公式计算结果与实际装药量存 在较大误差的现状，分别运用M A T L A B 及1 s t O p t 数值分析软件对装药量公式的参数进行拟合，采用P o l y n o m i a l 拟合方法对装药量公式进行拟合，并从残差平方和 S S E 、均方根误差 R M S E 以及相对误差均值三个 方面进行拟合结果的有效性和优越性分析。结果表明根据1 s t O p t 优化公式计算出的残差平方和、均方根误 差、相对误差均值分别为0 ．5 9 7 、0 ．0 8 1 4 、0 ．0 2 2 5 ，得出的装药量与实际更吻合，其结果明显优于其他方法。 关键词受限自由面；巷道爆破；装药量公式；拟合 中图分类号T D 2 3 5 ．3文献标识码A文章编号1 0 0 1 4 8 7 X 2 0 1 7 0 1 0 0 3 7 0 5 O p t i m i z a t i o nS t u d yo nC a l c u l a t i o nF o r m u l ao fE x p l o s i v e C h a r g ei nL i m i t e dF r e eS u r f a c eB l a s t i n g L lQ i 一r u e ，Z H A N GC h e n g - j u n ，W UZ h e n g y u ，C H E NY i n g ，W E IJ i a r u i S c h o o lo fR e s o u r c e sa n dS a f e t yE n g i n e e r i n g ，C e n t r a lS o u t hU n i v e r s i t y ，C h a n g s h a4 1 0 0 8 3 ，C h i n a A b s t r a c t T oi m p r o v et h es p e e do ft u n n e le x c a v a t i o nb l a s t i n ga n dt oc o n t r o lt h et u n n e lc o n t o u rp r e c i s e l y ，a c c o r d - i n gt ot h ec a l c u l a t i o nr e s u l t so ft h ee x i s t i n gf o r m u l a s ，t h e r ei sab i ge r r o rb e t w e e nt h et h e o r e t i c a le x p l o s i v ec h a r g ea n d t h ea c t u a le x p l o s i v ec h a r g e ．T h e n ，t h en u m e r i c a la n a l y s i ss o f t w a r eM A T L A Ba n d1s t O p tw e r eu s e dt of i tt h ep a r a m e t e r so ft h ec a l c u l a t i o nf o r m u l ao fe x p l o s i v ec h a r g e ，a n dan e wf o r m u l aw a sf i t t e db yP o l y n o m i a l ．T h ev a l i d i t ya n ds u p e f i o r i t yo ft h ef i t t e dr e s u l t sw e r ea n a l y z e df r o mt h r e ea s p e c t s ，n a m e l yr e s i d u a ls q u a r es u m S S E ，r o o tm e a ns q u a r eo r r o r R M S E a n dm e a nr e l a t i v ee r r o r ．T h er e s u l t ss h o wt h a tt h e1s t O p to p t i m i z a t i o nf o r m u l aw a sm o r ep r e c i s et h a n o t h e rm e t h o d s ，b yw h i c ht h eS S E ，t h eR M S Ea n dt h em e a nr e l a t i v ee r r o rw e r eo b t a i n e da s0 ．5 9 7 ，0 ．0 8 1 4a n d 0 ．0 2 2 5 ，r e s p e c t i v e l y ．F i n a l l y ，t h ec a l c u l a t e dr e s u l t sf i t t e dt h ea c t u a lc h a r g ea m o u n tw e l l ． K e yw o r d s l i m i t e df r e es u r f a c e ；t u n n e lb l a s t i n g ；c a l c u l a t i o nf o r m u l ao fc h a r g ea m o u n t ；f i t t i n g 在岩石爆破技术中，隧道和平巷的掘进速度和 爆破效果受到特别的重视，掘进速度的快慢对施工 进度起着决定性的作用。为提高掘进速度，改善爆 破效果，许多学者针对掏槽方式做了一些研究，余永 强等通过实验提出了分阶楔形加中心孔的掏槽方 式，得出了增加中心孔能提高掏槽孔利用率和掘进 收稿日期2 0 1 6 1 1 2 7 作者简介李启月 1 9 6 8 一 ，男，博士，湖南衡阳人，1 9 9 0 年毕业于沈 阳黄金学院矿山系采矿专业，现任教授、博士生导师，主要 从事非煤矿山采矿方法与岩土工程爆破方面的教学与研 究工作， E m a i l h e a d j u n 2 0 1 0 1 6 3 ．c o m 。 基金项目国家自然科学基金 5 1 3 7 4 2 4 3 、4 1 3 7 2 2 7 8 速度o 。杨国梁等研究出复式楔形深孔掏槽爆破 的方式嵋J 。陈士海等通过在坚硬岩石巷道进行中 深孔掏槽爆破实验研究得出复式掏槽完全适合于中 深孔爆破∞J 。单仁亮等提出了一种新的掏槽方 法准直眼掏槽方式H 1 。汪海波等认为增大中 心空孔直径能提高竖井井筒坚硬岩石段施工的爆破 效率H j 。廖涛等提出切缝药包装药结构能有效改 善掏槽效果MJ 。另外一些学者对爆破参数进行优 化来改善爆破效果“ ’10 f 。以上研究是对掏槽形式及 爆破参数优化方面的研究，然而在受限自由面爆破 过程中，精确的装药量是控制巷道轮廓并提高巷道 万方数据 3 8爆破 2 0 1 7 年3 月 掘进速度的保障。很多学者在装药量方面做了一些 研究1 ’1 2 j ，其中UL a n g e f o r s 和BK i h l s t o r m 根据实验 数据推算出受限自由面爆破装药量的公式。13 I ，对巷 道爆破有一定的指导意义，但存在一定的误差。为 提高受限自由面爆破时装药量的精确性，运用 M A T L A B 及1s t O p t 数值分析软件对UL a n g e f o r s 和B K i h l s t o r m 提出的装药量公式从参数和公式两方面 进行优化。目前参数的预测方式主要有最小二乘 法4 | 、神经网络法[ 1 5 1 6 1 等算法，但这些算法比较复 杂，而且需要提供经验初值，从而影响求取参数的准 确性。M A T L A B 及1 s t O p t 在参数求解过程中使用 通用全局优化算法，克服了优化迭代时必须给出的 初值难题，使参数求解得到简化。 1 参数优化 为减小传统装药量在计算装药量过程中的误 差，运用M A T L A B 以及1s t O p t 两种方法对传统装药 量中的参数进行拟合得到优化参数，并构建了拟合 参数和实验矩阵。 1 ．1 构建参数以及实验矩阵 如果用￡表示每米装药量，“ 表示空心孔角的 一半，V m 表示炮眼到自由面的距离，B 1 1 1 表示矩 形槽口的宽度，矩形槽口的宽度是受限的为窄自由 面，向受限自由面爆破如图1 所示。UL a n g e f o r s 和 BK i h l s t o r m 根据实验数据得出装药量公式n 3 1 L 0 ．3 5V ／ s i nt 、 。 1 炮 图1同党限臼由向爆破l 刳 F i g ．1B l a s t i n gt o w a r d sn a r r o wo p e n i n g 因为口可以用y 和B 表示，因此装药量可以表 示成y 和B 的函数，即 L - o ．3 5 ∥ 南 埘 ㈦ 因此可以构建参数a 、b ，则装药量公式 ￡⋯∥ 蠡 6 3 根据UL a n g e f o r s 和BK i h l s t o r m 的进行试验得 到了使岩石完全破碎的装药量数据3 | ，绘制成装药 量试验数据表见表1 。 表1 装药量试验数据 T a b l e1 E x p e r i m e n t a ld a t ao fc h a r g e 0 ．1 0 0 ．1 5 0 ．2 0 0 ．2 5 0 ．3 0 O ．3 5 0 ．4 0 0 ．5 0 O ．6 0 0 ．7 0 0 ．8 0 1 ．0 0 1 ．2 0 1 ．4 0 1 ．6 0 2 ．0 0 为使得不同的抵抗线长度和自由面长度和装药 量对应起来，将装药量建成3 9 0 的矩阵如表2 所 示，为了更方便的数据处理和过程操作用x 、l ，、z 替 代B 、y 、￡。 ∞粥加∞∞∞ c ；m L L 互t 嬲弘配加∞∞∞∞舳 c ；m m m L L 互王乱 墙巧∞∞加∞∞∞舳∞ m m m m m m L 互王i 毖∞∞∞∞如跚∞∞ m m m n L L L 王乱 屹坞撕巧∞加∞踯加∞ m m c ；c ；m m L L Z 钆 H尥孔∞∞∞∞加加 m m m n m L L Z 孓 ∞怕拍弘撕加如如∞ m m m m m m L L 互 加如如酪帅∞加 m m m c ；m m L Z 嘶B抖”∞如加∞ m m m c ；c ；m L 互 吣博药∞∞加∞ m m m m m L L 他∞∞∞如 眦眦眦L L 万方数据 第3 4 卷第1 期李启月，张成君，吴正宇，等受限自由面爆破装药量计算公式的优化研究 3 9 表2 装药量实验数据矩阵 T a b l e2E x p e r i m e n t a ld a t am a t r i xo fc h a r g e 变量 X y Z 数据 0 ．81 ．41 ．41 ．4 1 ．60 ．81 ．01 ．2 4 ．80 ．60 ．90 ．2 1 ．4 2 ．0 4 ．0 1 ．2 M A T L A B 参数拟合 根据推算的装药量公式可以通过M A T L A B 中 自定义方程进行参数拟合。在拟合方程式中自定义 方程“Z a Y [ X 2 ／ X 2 4 护 ] 6 ’’，式中x 、y 、z 分别 为B 、y 、L 。在拟合选项 F i to p t i o n s 中选用非线性 最小二乘法进行拟合，参数口、b 的上下限不限定，在 d a t a 选项中Xd a t a 、Yd a t a 、Zd a t a 分别对应x 、l ，、z 。 通过设定的拟合方式得到拟合参数n 、b 的值， 并得到拟合后装药量公式的残差平方和、均方根 误差以及相关系数并计算了相对误差均值如表3 所示。 表3 参数及拟合优度 T a b l e3P a r a m e t e r sa n dg o o d n e s so ff i t 1 ．3 l s t O p t 参数拟合 传统参数优化方法在预测过程中提供经验初 值，从而影响求解参数的准确性。1s t O p t 中的通用 全局优化算法，克服了优化迭代时必须给出初值的 难题，使参数优化得到简化。 在1s t O p t 参数优化过程给定了相应的参数、变 量、和公式，并按照软件自带的编程规则编写相应的 参数求解程序代码，格式如下 T i t l e 受限自由面装药量公式的参数优化 P a r a m e t e r sa ，b ； V a r i a b l ex ，Y ，Z ； F u n c t i o nZ a 术y ／ x 2 ／ x 2 4 术y “ 2 6 b ； d a t a ； X l ，X 2 ，X 3 ，⋯，X 9 0 ； Y I ，Y 2 ，Y 3 ，⋯，Y 9 0 ； Z l ，Z 2 ，Z 3 ，⋯，Z 9 0 ； 根据代码，将实测的装药量数据输入，选用通用 全局优化法模式进行自动拟合计算。得到拟合参数 值，以及残差平方和、均方根误差以及相关系数并计 算了相对误差均值如表4 所示，得出了装药量实测 值及计算值如图2 所示。 表4 参数及拟合优度 T a b l e4P a r a m e t e r sa n dg o o d n e s so ff i t 图2 装药量参数优化图 F i g ．2C h a r g ep a r a m e t e r so p t i m i z a t i o n 从图2 中可以看出拟合效果良好，证明优化后 的装药量公式得出的药量与实际药量基本吻合。 2 装药量公式拟合 根据装药量实验数据，通过M A T L A B 曲线拟合 工具箱，运用多项式 P o l y n o m i a l 拟合出新的装药量 公式来进行比较，从而得到更精确的装药量公式。 2 ．1 P o l y n o m i a l 拟合 将抵抗线以及自由面长度拟合成装药量的公式 即￡ 八B ，V ，拟合时对拟合参数进行设置，设定自 变量X 。、y ，的最高拟合次数，因为拟合次数越高误 差就越小，把D e g r e e s 选项中戈、Y 的最高拟合次数设 为5 。R o b u s t 稳定性控制 选项用于设定求解算法 所采用的回归模式，“O F F ”选项是常规的最小二乘 法进行拟合，“L A R ”选项是最小绝对残差法， 4 6 4 1 l 2 4 4 6} 1 3 3 0 0 1 O 5 0 1 2 0 O 0 l l 2 6 O O 0 O 5 O 1 1 3 0 0 O O O 2} O O O 万方数据 爆破2 0 1 7 年3 月 “B i s q u a r e ”为加权最小二乘法。加权最小二乘法是 在拟合过程中含有异常值时最健壮的拟合方法，与 常规最小二乘法相比它对异常值不敏感，少数的异 常值对整体的影响很小甚至不会影响整体的拟合效 果。因为本次拟合的装药量数据是进行多次实验得 出的最合理的数据因此基本上没有异常值，所以选 择常规最小二乘法进行拟合。在拟合过程中为了避 免计算过程中出现病态矩阵选择“C e n t e ra n ds c a l e ” 居中和缩放 。 在拟合工具箱d a t a 选项中X d a t a 、Y d a t a 、Z d a t a 分别对应前面的装药量数据x 、y 、z 后进行拟合。 通过拟合得到装药量曲面图如图3 图3 装药量曲面图 F i g ．3 C h a r g ec u r v e ds u r f a c ef i g u r e 2 ．2 拟合结果 通过实验数据基于M A T L A B 拟合出最高拟合 次数为5 的装药量公式为 L 厂 B ，V P 0 0 P l o B P 0 1 V P 2 0 8 2 P 1 1 B V P 0 2 v 2 ⋯ P 0 5 矿 4 式中，p i 伪多项式系数。并给出了多项式系数 表如表5 所示，以及处理后的残差平方和、均方根误 差以及相关系数并计算了相对误差均值见表6 。 3 拟合结果的有效性及优越性分析 为了得到更精确的装药量公式以及对参数拟合 和装药量公式拟合结果的可靠性作出评价，对拟合 结果从有效性和优越性两方面进行分析。 3 ．1 拟合结果的有效性分析 对拟合结果的有效性评判常采用相关系数尺2 进行检验。在统计中常用相关系数R 来衡量变量 之间线性关系的强弱，对变量戈。、Y 。不全为零时相关 系数为叫 尺。蒜2 一 √厶i 1 ” 7 “‘“ 5 式中戈为变量戈的平均值；Y 为变量Y 均值。 拟合效果的判定系数为7 。 R 2 监L x x L y y 潞 6 式中S S R 表示回归值关于它们的平均值的离 差平方和；S S T ∑ y 。一y 2 反应的是样本关于它 们的均值Y 的离差平方和；判定系数R 2 反应了回归 平方和在总离差平方和中的比例。因此，I 别愈大， 相关程度愈高，线性回归愈好，相关系数的相关性判 定表如7 所示。 表5 拟合公式系数 T a b l e5C o e f f i c i e n to ff i t t i n gf o r m u l a 表6 拟合优度 T a b l e6G o o d n e s so ff i t 表7 相关性判定表 T a b l e7C o r r e l a t i o nd e c i s i o n RJ 的取值范围 J R l 的意义 0 ．0 0 ～0 ．1 9 0 ．2 0 ～0 ．3 9 0 ．4 0 ～0 ．6 9 0 ．7 0 ～0 ．8 9 0 ．9 0 ～1 ．0 0 极低相关 低度相关 中度相关 高度相关 极高相关 此次拟合三种方法的拟合过程中相关系数分别 为0 ．9 8 6 、0 ．9 9 7 6 、0 ．9 9 4 9 ，均接近1 属于极高相关 性，所以拟合的效果好相关性强，因此这三种方式得 出的装药量公式均是有效地。 3 ．2 拟合结果的优越性分析 为了验证拟合结果的优越性，通常从残差平方 和 S S E 、相对误差均值、均方根误差三个指标进行 比较，残差平方和计算公式为 S S E ∑ 厶一￡。 2 7 J _ 。 i 1 式中L 。为实测的实际装药量；厶为装药量公 式得到的药量。将通过三种方式得出的装药公式与 万方数据 第3 4 卷第1 期李启月，张成君，吴正宇，等受限自由面爆破装药量计算公式的优化研究 4 l 传统装药量公式进行比较，结果见表8 。 表8 不同装药量公式优度比较表 T a b l e8G o o d n e s sc o m p a r i s o no fd i f f e r e n tc h a r g ef o r m u l a 从表8 可以看出，三种方式得到的装药量公式 计算的装药量比传统装药量计算的装药量更接近实 际值，其中通过1s t O p t 参数拟合出的新装药量公式 以及P o l y n o m i a l 拟合出的新装药量公式计算出的结 果更精确，由于P o l y n o m i a l 拟合出的装药量公式比 较复杂，实用性需要进一步探讨。因此通过优越性 分析得出的更精确的装药量公式为 L 0 ．4 1 0 4 V ／ s i n 秽 7 巧 8 4 结论 1 UL a n g e f o r s 和BK i h l s t o r m 提出的单孔受限 自由面下爆破装药量公式的计算装药量与工程实际 装药量存在着一定误差。因此针对上述问题，通过 三种不同拟合方法进行了优化分析，并证明了优化 后的装药量计算公式是可行的。 2 从装药量公式参数优化及装药量公式拟合 两个方面进行了优化对比。证明拟合过程中三种方 法的拟合的装药量公式的线性相关系数均接近1 ， 属于高度相关，拟合的效果好。 3 通过优越性分析最终得出的装药量公式为 L 0 ．4 1 0 4 W s i n 秽 7 力，利用该公式计算出装药量的 残差平方和、均方根误差以及误差均值分别为 0 ．5 9 7 、0 ．0 8 1 4 、0 ．0 2 2 5 ，这三项评优指标均为最低， 计算出的装药量更为精确，对工程实际具有良好的 决策依据及广泛的应用价值。 4 仅对单孔受限自由面下爆破装药量公式进 行优化，对于多孔条件下的受限自由面爆破，爆破的 装药量受炮孔个数、炮孔间距、炮孔排距以及起爆时 间等因素的影响较大，例如，对于四个或更多炮眼爆 破，当炮眼间距等于最小抵抗线时，每个炮眼的装药 量可减少为计算所得装药量的8 0 ％左右3 。。因此 对于多孑L 条件下的受限自由面爆破仍需进一步研 究。此外此次拟合得出的参数是基于有限的实验数 据进行的优化，其结果可能存在一定的偏差，仍需要 更多的实验数据来验算其准确性。 参考文献 R e f e r e n c e s 余永强，吴帅锋，褚怀宝，等．非均质硬岩巷道爆破研 究方法[ J ] ．工程爆破，2 0 1 3 ，1 9 6 2 5 - 2 8 ． Y UY o n g q i a n g ，W US h u m f e n g ，C H UH u a i - b a o ，e ta 1 ． S t u d yo nb l a s t i n gm e t h o do fr o a d w a ye x c a v a t i o ni nh e t e r o - g e n e o u sh a r dr o c k [ J ] ．E n g i n e e r i n gB l a s t i n g ，2 0 1 3 ， 1 9 6 2 5 2 8 ． i nC h i n e s e 杨国梁，姜琳琳，杨仁树．复式楔形深孔掏槽爆破研 究．[ J ] ．中国矿业大学学报，2 0 1 3 ，4 2 5 7 5 4 - 7 6 0 ． Y A N GG u o l i a n g ，J I A N GL i n l i n ，Y A N GR e n s h u ．I n v e s t i g a t i o no fc u tb l a s t i n gw i t hd u p l e xw e d g ed e e ph o l e s [ J ] ． J o u r n a lo fC h i n aU n i v e r s i t yo fM i n i n g T e c h n o l o g y ． 2 0 1 3 ，4 2 5 7 5 4 _ 7 6 0 ． i nC h i n e s e 陈士海，魏海霞，薛爱芝．坚硬岩石巷道中深孔掏槽爆 破实验研究[ J ] ．岩石力学与工程学报，2 0 0 7 ，2 6 1 3 4 9 8 ．3 5 0 2 ． C H E NS h i h a i ，W E IH a i x i a ，X U EA i - z h i ．T e s t i n gs t u d y o nm i d d l ed e e pc u t h o l eb l a s t i n gi nh a r dr o c kt u n n e l [ J ] ． C h i n e s eJ o u r n a lo fR o c kM e c h a n i c s a n dE n g i n e e r i n g ， 2 0 0 7 ，2 6 1 3 4 9 8 - 3 5 0 2 ． i nC h i n e s e 单仁亮，黄宝龙，高文蛟，等．岩巷掘进准直眼掏槽爆 破新技术应用实例分析[ J ] ．岩石力学与工程学报， 2 0 1 1 ，3 0 2 2 2 4 - 2 3 2 ． S H A NR e n l i a n g ，H U A N GB a o l o n g ，G A OW e n - j i a o ，e t a 1 ．C a s es t u d yo fn e wt e c h n o l o g ya p p l i c a t i o no fq u a s i p a r - a l l e lc u tb l a s t i n gi nr o c kr o a d w a yd r i v a g e [ J ] ．C h i n e s e J o u r n a lo fR o c kM e c h a n i c sa n dE n g i n e e r i n g ，2 0 11 ， 3 0 2 2 2 4 - 2 3 2 ． i nC h i n e s e 汪海波，宗琦，赵要才．立井大直径中空孔直眼掏槽 爆炸应力场数值模拟[ J ] ．岩石力学与工程学报， 2 0 1 5 ，3 4 S 1 3 8 1 3 8 6 ． W A N GH a i b o ，Z O N GQ i ，Z H A OY a o c a i ．N u m e r i c a la n a l y s i sa n da p p l i c a t i o no fl a r g ed i a m e t e rc a v i t yp a r a l l e l c u ts t r e s sf i e l di nv e r t i c a ls h a f t 『J ] ．C h i n e s eJ o u r n a lo f R o c kM e c h a n i c sa n dE n g i n e e r i n g ，2 0 1 5 ，3 4 S 1 3 8 1 3 8 6 ． i nC h i n e s e 下转第6 6 页 ] J 1J]J]J]j]j]J]J]J ] 1 1 2 2 3 3 4 4 5 5 口 阻 心 心 b b M M 瞪 瞪 万方数据 6 6爆破 2 0 1 7 年3 月 上接第4 1 页 [ 6 ] [ 6 ] [ 7 ] [ 7 ] [ 8 ] [ 8 ] [ 9 ] [ 9 ] [ 1 0 ] [ 1 0 ] 廖涛，蒲传金，王俊青，等．岩巷切缝药包楔形掏槽 爆破模型试验与分析[ J ] ．金属矿山，2 0 1 5 ，4 6 4 4 0 4 4 ． L I A OT a o ，P UC h u a n - j i n ，W A N GJ u n q i n g ，e ta 1 ．M o d e l e x p e r i m e n ta n da n a l y s i so fw e d g e s h a p e dc u t t i n gb l a s r i n gw i t hc u t t i n gs e a mc a r t r i d g ei nr o c kt u n n e l [ J ] ．M e t a l M i n e ，2 0 1 5 ，4 6 4 4 0 4 4 ． i nC h i n e s e 曹胜祥，邱贤阳，史秀志．大直径深孔爆破参数优化研 究[ J ] ．有色金属工程，2 0 1 5 ，5 S 1 1 6 0 1 6 4 ． C A OS h e n g - x i a n g ，Q I UX i a n y a n g ，S H IX i u z h i ．R e s e a r c h o nt h ep a r a m e t e r so p t i m i z a t i o no fl a r g e d i a m e t e rl o n g h o l e b l a s t i n g [ J ] ．N o n f e r r o u sM e t a l sE n g i n e e r i n g ，2 0 1 5 ， 5 S 1 1 6 0 1 6 4 ． i nC h i n e s e 张建华，王琪，黄其冲．基于A H P 和灰色关联分析法 的爆破参数优化[ J ] ．爆破，2 0 1 4 ，3 1 4 7 2 - 7 6 ． Z H A N GJ i a n h u a ，W A N GQ i ，H U A N GQ i - c h o n g ．O p t i m i z a t i o no fb l a s t i n gp a r a m e t e r sb a s e do na n a l y s i sh i e r a r c h y p r o c e s sa n dg r e yc o r r e l a t i o na n a l y s i sm e t h o d [ J ] ．B i a s t i n g ，2 0 1 4 ，3 1 4 7 2 - 7 6 ． i nC h i n e s e 徐帅，彭建宇，李元辉，等．急倾斜薄矿脉中深孔落 矿爆破参数优化[ J ] ．爆炸与冲击，2 0 1 5 ，3 5 5 6 8 2 6 8 8 ． X US h u a i ，P E N GJ i a n y u ，L IY u a n h u i ．e ta 1 ．B l a s t i n gp a r a m e t e ro p t i m i z a t i o no fm e d i u m d e p t hh o l ec a v i n gf o r s t e e p l yi n c l i n e d t h i nv e i n s [ J ] ．E x p l o s i o na n d S h o c k W a v e s ，2 0 1 5 ，3 5 5 6 8 2 _ 6 8 8 ． i nC h i n e s e 张钦礼，郑晶晶，张德明，等．新桥硫铁矿凿岩爆破参 数优化实验研究[ J ] ．爆破，2 0 0 9 ，2 6 3 1 5 ． Z H A N GQ i n l i ，Z H E N GJ i n g - j i n g ，Z H A N GD e m i n g ，e t a 1 ．T e s tr e s e a r c ho nd r i l l i n ga n db l a s t i n gp a r a m e t e r so p t i m i z a t i o ni nX i nQ i a op y r i t em i n e [ J ] ．B l a s t i n g ，2 0 0 9 ， [ 1 1 ] [ 1 2 ] [ 1 2 ] [ 1 3 ] [ 1 4 ] [ 1 5 ] [ 1 5 ] [ 1 6 ] [ 1 6 ] [ 1 7 ] 2 6 3 1 - 5 ． i nC h i n e s e 孟国良，金勇．基于爆破漏斗实验装药量的计算公 式[ J ] ．煤矿爆破，2 0 1 0 ，4 5 2 1 5 1 7 ． M E N GG u o - l i a n g ，J I NY o n g ．S t u d yo nc a l c u l a t i o nf o r m u l ao fe x p l o s i v ec h a r g eo nb l a s t i n gt e s t [ J ] ．C o a lM i n e B l a s t i n g ，2 0 1 0 ，4 5 2 1 5 1 7 ． i nC h i n e s e 郑明焦．中深孔爆破炸药量计算方法的应用[ J ] ．中 国矿山工程，2 0 0 5 ，4 3 4 1 5 1 7 ． Z H E N GM i n g j i a o ．A p p l i c a t i o no fe x p l o s i v ea m o u n tc a l c u l a t i o nm e t h o di nm e d i u m l e n g t hh o l eb l a s t i n g [ J ] ．C h i n aM i n eE n g i n e e r i n g ，2 0 0 5 ，4 3 4 1 5 1 7 ． i nC h i n e s e UL a n g e f o r s ，BK i h l s t r o m ．T h em o d e mt e c h n i q u eo fr o c k b l a s t i n g [ M ] ．W i l e y ，1 9 6 3 ． G U P T AR a j e n d e r ．A d v a n c e dc o a lc h a r a c t e r i z a t i o n ar e ． v i e w [ J ] ．E n e r g y ＆F u e l s ，2 0 0 7 ，2 1 2 4 5 1 4 6 0 ． 孙文彬，刘希亮，谭正龙，等．基于抛掷爆破预测的B P 神经网络参数优化方法[ J ] ．煤炭学报，2 0 1 2 ， 3 7 S 1 5 9 6 4 ． S U NW e n b i n ，L I UX i l i a n g ，T A NZ h e n g l o n g ，e ta 1 ．P a r a m e t e ro p t i m i z a t i o no fB P n e u r a ln e t w o r kb a s e do nt h e f o r e c a s to fc a s tb l a s t i n g [ J ] ．J o u r n a lo fC h i n aC o a lS o c i e t y ，2 0 1 2 ，3 7 S 1 5 9 6 4 ． i nC h i n e s e 任芳，杨兆建，熊诗波．国内外煤岩界面识别技术 研究动态综述[ J ] ．煤，2 0 0 1 4 4 7 - 5 3 ． R E NF a n g ，Y A N GZ h a o - j i a n ，X I O N GS h i b o ．C o m m e n t o nt h er e s e a r c hd y n a m i co fr o c ki n t e r f a c ei d e n t i f i c a t i o n t e c h n o l o g ya th o m ea n da b r