镇沅金矿中深孔爆破参数数值模拟研究.pdf

第3 6 卷第1 期 爆破 V 0 1 ．3 6 N o 1 2 0 1 9 年3 月B L A S T I N G M a r ．2 0 1 9 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 9 ．0 1 ．0 1 2 镇沅金矿中深孑L 爆破参数数值模拟研究 姜永恒1 ，雷恒永2 ，杨⋯、2 ，唐学义1 ，宋士生1 ，叶光祥3 1 ．长春黄金研究院有限公司，长春1 3 0 0 1 2 ；2 ．云南黄金有限责任公司镇沅分公司，普洱6 6 5 0 0 0 ； 3 ．赣州有色冶金研究所，赣州3 4 1 0 0 0 摘要镇沅金矿一直采用无底柱分段崩落法进行回采，但爆破中存在大块率高、缩口比较严重等问题。 为了确定合理的爆破参数，运用A N S Y S ／L S ．D Y N A 数值模拟软件对扇形炮孔爆破过程进行数值模拟，再现扇 形炮孔爆破应力场变化过程，分析了爆炸仿真过程中的V o nM i s e s 有效应力峰值，确定了中深孔爆破最优的 最小抵抗线和最大孔底距，为中深孔爆破提供理论依据。并对此参数进行了现场试验，现场中深孔爆破后， 爆破效果良好，进一步证明了所确定的最小抵抗线为1 ．5m ，最大孔底距为2 ．0m 为最优爆破参数。 关键词中深孔爆破；抵抗线；孔底距；爆破参数；数值模拟 中图分类号T D 2 3 5 ．3文献标识码A 文章编号 1 0 0 1 4 8 7 X 2 0 1 9 0 1 0 0 7 7 0 7 N u m e r i c a lS i m u l a t i o no fZ h e n y u a nG o l dM i n eB l a s t i n gP a r a m e t e r s r A N GY o n g ．h e n 9 1 ，L E It t e n g y o n 9 2 ，Y A N GJ i e 2 ，T A N GX u e y i l ，S O N GS h i s h e n 9 1 ，Y EG u a n g x i a n 9 3 1 ．C h a n g c h u nG o l dR e s e a r c hI n s t i t u t e ，C h a n g c h u n13 0 0 12 ，C h i n a ； 2 ．Z h e n y u a nB r a n c hC o m p a n yY u n n a nC oL t d ，P u e r6 6 5 0 0 0 ，C h i n a ； 3 ．G a n z h o uN o n f e r r o u sM e t a l l u r g yR e s e a r c hI n s t i t u t e ，G a n z h o u3 4 1 0 0 0 ，C h i n a A b s t r a c t T h en o n p i l l a rs u b l e v e lc a v i n gh a sb e e nu s e di nZ h e n y u a nG o l dM i n e ．H o w e v e r ，t h e r ea r es o m ep r o b l e m ss u c ha sh i g hb u l kr a t ea n ds e r i o u ss h r i n k a g e ．I no r d e rt od e t e r m i n er e a s o n a b l eb l a s t i n gp a r a m e t e r s ，t h en u m e f i c a ls i m u l a t i o no fb l a s t i n gp r o c e s sw i t hf a n s h a p e db o r e h o l el a y o u tw a ss t u d i e db yA N S Y S ／L S D Y N A ．T h ec h a n g eo f t h es t r e s sf i e l df o rt h eb l a s t i n gw a so b t a i n e d ．T h ep e a kv a l u eo fV o nM i s e se f f e c t i v es t r e s sd u r i n gt h ee x p l o s i o ns i m u l a t i o nw a sa n a l y z e d ．T h em i n i m u mb u r d e na n dm a x i m u mh o l eb o t t o mp i t c ho ft h eb l a s t i n gw e r ed e t e r m i n e d ．T h er e s e a r c hr e s u l t sc a np r o v i d et h e o r e t i c a lb a s i sf o ri n d u s t r i a lt e s t ss u b s e q u e n t l y ，a n dt h eb l a s t i n ge f f e c ti sv e r yw e l l ．I ti s f u r t h e rp r o v e dt h a tt h em i n i m u mb u r d e no f1 ．5m e t e r sa n dt h em a i m u mh o l eb o t t o mp i t c ho f2 ．0m e t e r sa r et h eo p t i m a lb l a s t i n gp a r a m e t e r s ． K e yw o r d s d e e p h o l eb l a s t i n g ；b u r d e n ；h o l eb o t t o mp i t c h ；b l a s t i n gp a r a m e t e r ；n u m e r i c a ls i m u l a t i o n 扇形中深孔爆破在国内外金属矿山中已经广泛 应用，其主要特点有生产效率高、工作能力大、劳动 作业安全、生产成本低等’1J 。实践证明，影响爆破效 果的因素很多，除了爆破工艺以外还包括有地质构 造、爆破参数、装药方式、炸药性能等心，3J 。在爆破参 数设计中孔底距和最小抵抗线是两个关键参数HJ ，由 收稿日期2 0 1 8 1 l 0 5 作者简介姜永恒 1 9 8 9 一 ，男，辽宁省建平县人，工程师、硕士研究 生， E m a i l j i a n g y o n g h e n g 1 2 6 ．c o r n 。 于爆破过程的复杂性，在现有的实验条件和测试手段 下，难以对这两个重要参数准确确定。仅凭借于经验 类比法确定爆破参数，在实际生产中针对性不强，具 有一定的局限性。因此在实际生产中需要根据岩石 条件进行数值模拟试验，对比分析确定合理的爆破参 数以降低矿石损失、贫化率，提高生产效率。 镇沅金矿矿岩主要为凝灰质粉砂岩及岩屑石英 杂砂岩，凝灰质杂砂岩岩石力学强度低，为泥质粉砂 状结构，中厚层状构造，岩屑石英杂砂岩，岩石胶结 万方数据 7 8爆破2 0 1 9 年3 月 紧密坚硬，造成岩体结构不均，给爆破参数的确定带 来较大的困难。镇沅金矿一直采用无底柱分段崩落 法进行回采，但爆破中存在大块率高、缩口比较严重 等问题。为了提高炸药的利用率，改善爆破效果，采 用A N S Y S ／L S ．D Y N A 数值模拟软件对扇形炮孔爆破 过程进行数值模拟，对爆破参数进行优化设计研究， 确定出中深孔爆破的最优最小抵抗线和最大孔底 距，为矿山生产爆破提供依据。 1 爆破参数数值模拟计算 1 ．1 模型的建立 在爆破过程中，岩石内部的受力状态非常复杂， 难以用M i s e s 屈服条件判断岩石内部某点是否发生 破坏【5o ；然而，自由面附近的岩石可认为处于接近 简单拉伸的状态。所以，在模型建立中选取M i s e s 屈服准则和拉伸准则共同作用，从而可以采用M i s e s 屈服准则近似来判断岩石是否发生破坏。根据M i s e s 屈服准则M 7 J ，对于爆破模拟结果，可以提取自 由面上关键位置的有效应力曲线，读取曲线的峰值， 用该峰值与岩石的动态抗拉强度做对比，确定岩石 是否已经爆开哺J ，通过分析镇沅金矿的矿岩特性， 确定其矿岩动抗拉强度约为9 0M P a 。扇形孔爆破 示意图见图1 。 扇 图1 扇形孔爆破示意图 F i g ．1 T h eb l a s t i n gd i a g r a mo ft h ef a n - s h a p e dh o l e 注释岩土材料理论上不能单独采用M i s e s 屈 服准则来判断岩石是否发生破坏，在计算中采用了 M i s e s 屈服准则和拉伸准则同时对岩体进行控制，在 A N S Y A 计算过程中两个准则同时起作用，在计算中 可以采用M i s e s 屈服准则近似来判断岩石是否发生 破坏 爆破最小抵抗线是指第一排扇形炮孔到前排面 自由面的距离，图1 中爆破抛掷方向即为前排面自 由面，为了简化模型和计算过程，用一个炮孔来模拟 爆破过程。 确定最小抵抗线的爆破模型如图2 所示，由于 爆破模拟的模型只包含凿岩硐室顶板以上部分，模 型不包含凿岩硐室，对结果几乎不产生影响，却大大 节省计算耗时，故在计算过程中，除了前排面自由面 及下表面之外，其他表面均设置无反射边界条 件旧J ，以表示模型在该方向上无限延伸。 凿岩巷道顶板 图2 确定最小抵抗线的爆破模型示意图 F i g ．2 T h eb l a s t i n gm o d et od e t e r m i n e t h em i n i m u mr e s i s t a n c el i n e 确定最大孔底距的爆破模型如图3 所示，其中 a 为整个模型， b 为切开一个剖面看到的内部情 况。模型的上表面、左右表面及里表面均设置无反 射边界条件，表示这些方向模型无限延伸。模型的 下表面即为凿岩硐室顶板，为节省计算耗时，未将凿 岩硐室建在计算模型中，且仅用两个炮孔来进行模 拟，未将所有扇形炮孔都放在计算模型中。 岩石和炸药选择s o l i d l 6 4 单元归J ，岩石选木m a t _ p l a s t i ck i n e m a t i c 材料，炸药选用木m a t h i g h e x p l o - s i v e b u r n ，有状态方程控制其起爆过程0 | ，岩石与 炸药的参数见表1 和表2 。 表1 岩石参数 T a b l e1R o c kp a r a m e t e r 。， 密度／ 弹性模量泊松抗拉强度／动态抗拉 石臼 k g ．m 一3 G P a 比M P a 强度／M P a 镇沅金矿 2 7 0 07 ．50 ．3 01 3 ．99 0 万方数据 第3 6 卷第1 期 姜永恒，雷恒永，杨杰，等镇沅金矿中深孔爆破参数数值模拟研究 7 9 表2 乳化炸药参数 T a b l e2E x p l o s i v ep a r a m e t e r s a 模型整体图 a M o d e ld i a g r a m b 模型剖面图 b M o d e lp r o f i l e 风扇炮孔 图3 确定最大孔底距的爆破模型示意图 F i g ．3 T h eb l a s t i n gm o d et od e t e r m i n e t h em a x i m u mh o l es p a c i n g 1 ．2 最小抵抗线的确定 根据生产经验，首先取最小抵抗线为1 。3m ，若 能爆开岩体，则将最小抵抗线每次增加0 ．2m ，直至 岩体不能爆开为止，本次模拟选取最小抵抗线分别 为1 ．3m 、1 ．5m 、1 ．7m 。 爆破过程中的有效应力分布情况见如图4 所 示，从图4 可看出，扇形炮孔从孔口起爆，爆炸应力 场从孔口向炮孔内部传播1 ’1 2 1 ，之后应力场到达前 排自由面，最后爆炸应力场逐渐衰减。 提取前排自由面的单元的有效应力曲线，读取 其峰值，采用V o nM i s e s 有效应力屈服准则判断单 元的破坏情况。图5 显示了最小抵抗线分别为 1 ．3m 、1 ．5m 、1 ．7m 时自由面单元及有效应力曲 线。 F r i n g eL e v e l s 95 2 0 e 0 8 ss e s s l 7 ．6 1 6 1 8 一l 66 “e 0 8 雪 图4 爆破过程有效应力分布图 F i g ．4 T h ee f f e c t i v es t r e s sd i s t r i b u t i o no ft h eb l a s t i n g 通过图5 可以看出，最小抵抗线为1 ．3m 时，单 元的有效应力峰值为9 8 ．2M P a ，大于岩石动态抗拉 强度9 0M P a ，岩石被破坏；最小抵抗线为1 ．5m 时， 单元的有效应力峰值为9 5 ．5M P a ，大于岩石动态抗 拉强度9 0M P a ，岩石被破坏；最小抵抗线为1 ．7m 时，单元的有效应力峰值为8 8 ．4M P a ，岩石未破坏。 因此，应选择爆破最小抵抗线为1 ．5m 。 1 ．3 最大孔底距的确定 在最小抵抗线为1 ．5m 的情况下，根据矿山经 验，初定孔底距为1 ．8m ，最大孔底距每次增加0 ．2 ， 直至岩体不能爆开为止。本次模拟所选取的最大孔 底距分别为1 ．8m 、2 ．0m 、2 ．2m 。 爆破过程中的有效应力分布情况见图6 ，从图 中可看出，扇形炮孔从孔口起爆，爆炸应力场首先在 炮孔口产生，并向炮孔深部传播，然后两炮孔的应力 场开始叠加，共同作用破碎岩石2 | 。 F r i n g eL e v e l s 。，。z e o s l 36 8 Ie o , 1 32 2 l e 0 8 鞠 图7 显示了关键位置单元的有效应力瞳线，从 图中可以看出当孔底距为1 ．8m 时，单元的有效应 力峰值为9 9 ．8M P a ，大于岩石动态抗拉强度 9 0M P a ，岩石被破坏；当孔底距为2 ．0m 时，单元的 有效应力峰值为9 0 ．8M P a ，大于岩石动态抗拉强度 9 0M P a ，岩石被破坏；当孔底距为2 ．2m 时，单元的 有效应力峰值为8 6 ．3M P a ，岩石未破坏，此处出现 大块的可能性较大。因此，应选择最大孔底距为 2 ．0m 。 2 工程实践 根据数值模拟的采场爆破参数在现场进行试 验，试验采用的爆破参数为数值模拟所确定的最优 爆破参数，即最小抵抗线为1 ．5m ，最大孔底距为 2 ．0m ，并在试验中根据实际的爆破效果，对爆破参 数进行适当的调整。 ●●．■ o；甾疆■露叠曩■■■■ s 一 一 一 一 一 一 一 一 一 一 一 e i i i 7 0 v O O 0 0 0 0 0 0 0 0 O m 十 ； e e e e e e e ；8 6 4 2 9 7 6 3 O g 2 7 ●6 ●6 O 5 4 2 0 }8 4 ，7 4 O O 5 O ■■- 器誊叠■量叠誓I■一 s 一 一 一 一 一 一 一 一 一 一 一 e i i i 7 0 j 0 0 0 O O O O 0 O m ； B e C e e e e 1 7 8 0 2 3 5 7 3 0 目8 l 4 7 1 4 7 O 3 8 O 5 j 6 0 3 6 0 3 6 O 6 5 4 4 3 2 2 l 6 0 万方数据 爆破 2 0 1 9 年3 月 镇沅金矿在冬瓜林3 0 8 1 采场1 6 3 3 分层1 2 一1 6 线沿脉第2 1 排～3 0 排进行现场工业试验，试验地 点如图8 所示。采场采用Y G Z - 9 0 型凿岩机进行钻 凿上向扇形中深孔，孔径为6 5m m ，从孔口起爆。此 采场矿石的岩性主要为构造角砾岩以及砂岩成矿， i 自由面单元 H 17 8 3 3 1 l ／ i ／ 卦”1 3 3 平均地质品位约为2 ．5g ／t 。采场内每次爆破两排 炮孔，采用逐排起爆的形式，每排炮孔之间采用毫秒 微差起爆，雷管采用跳段使用，当前面两排炮孔爆破 并出矿结束后再近些后面炮孔的钻凿和爆破工作。 其中第2 1 、2 2 排炮孑L 的布置形式示意图如图9 所示。 T i m e ／ I i 一0 3 1 a 1 ．3 m 抵抗线时自由面单元及效应力曲线 a T h ef r e es u r f a c ee l e m e n ta n de f f e c tf o r c er e s i s t a n c el i n ew h e nt h er e s i s t a n c el i n ei s1 ．3m e t e r s 自由面单元 H 1 8 8 5 3 j H 1 1 1 1 5 3 b 1 ．5 m 抵抗线时自由面单元及有效应力曲线 b T h ef r e es u r f a c ee l e m e n ta n de f f e c tf o r c er e s i s t a n c el i n ew h e nt h er e s i s t a n c el i n ei s1 ．5m e t e r s 自由面单元 H 1 9 5 7 3 l c 1 ．7 m 抵抗线时自由面单元及有效应力曲线 c T h ef r e es u r f a c ee l e m e n ta n de f f e c tf o r c er e s i s t a n c el i n ew h e nt h er e s i s t a n c el i n ei s1 ．7m e t e r s F r i n g eL e w i s 63 5 0 e 0 8 1 5 “5 e 0 8 J 4 9 3 9 e 0 8 一薹 4 2 3 4 e “n 蓄 35 2 8 “”| 2 8 2 2 e “8 | 2 1 1 7 ””1 1 。4 ““”| 70 2 ， o ；t 0O O O e O O 一- 图5自由面单元的有效应力曲线图 F i g ．5 T h ee f f e c t i v es t r e s sc u r v eo ft h eu n i ti nf r e es u r f a c e F r i n g eL e v e l s 6 5 ‰ 0 8 】 57 8 l Ⅲ8 J ’”8 “” 4 3 3 6 。“一鋈 3 6 1 3 8 7 0 8 2 ”“”0 1 t 2 1 6 8 ⋯8 | 1 。4 4 5 “8 1 7 2 2 6 ””● 00 0 0 e 0 0 J F r i n g eL e v e l s 53 2 l c 0 8 1 47 3 0 e 0 8 1 35 4 7 。 0 8 一蓄 2 ’9 5 6 e 0 8 一I 2 “5 “8 一I 。7 7 ”8 { “”⋯8 | 5 ”2 ””1 o0 0 0 e O O J 图6 爆破过程有效应力分布图 F i g ．6T h ee f f e c t i v es t r e s sd i s t r i b u t i o no ft h eb l a s t i n g F r i n g eL e v e l s 35 S l e 0 8 3 1 9 0 8 】 2 ．8 4 0 e 0 8 】 2 ”“ 0 8 一蔷 7 7 “”8 | 4 2 仉 0 8 - ● 10 6 5 ”0 8 | 7 1 0 1 ”0 7 | 35 5 1 “0 7 | oO O O e 0 0 一● 万方数据 第3 6 卷第1 期姜永恒，雷恒永，杨杰，等镇沅金矿中深孔爆破参数数值模拟研究 4 6 2 0 0 0 0234 T i m e ／ E 一0 3 J a 1 ．8m 孔底距时关键位置单元及有效应力曲线 a T h ek e yp o s i t i o nu n i ta n de f f e c t i v es t r e s sc urvew h e nt h eh o l es p a c i n gi s1 ．8m e t e r s 10 0 关键位置宗 单元 益o ．0 8 H 2 5 3 8 2 i 三O ．0 6 高0 ．0 4 20 ．0 2 0l234 T i m e ／ E 一0 3 b 2 ．0m 孔底距时关键位置单元及有效应力曲线 b T h ek e yp o s i t i o nu n i ta n de f f e c t i v es t r e s sc u r v ew h e nt h eh o l es p a c i n g i s2 ．0m e t e r s 1 ．0 0 曲0 ．0 8 关键位置立o0 6 单元 2 H 2 5 3 8 2 量0 0 4 童0 ．0 2 01234 乃m e ／ E 0 3 c 2 ．2m 孔底距时关键位置单元及有效应力曲线 C T h ek e yp o s i t i o nu n i ta n de f f e c t i v es t r e s sc u r v ew h e nt h eh o l es p a c i n g i s2 ．2m e t e r s 图7 关键位置单元的有效应力曲线图 F i g ．7 T h ee f f e c t i v es t r e s sc u r v eo ft h eu n i ti nf r e es u r f a c e 4 6 4 0 ＼锄 1 2 ‘＼ 图8 试验地点炮排设计 F i g ．8 T h eh o l ed e s i g ni nt e s ts i t e 4 0 4 0 1 6 4 6 1 6 4 3 1 6 3 3 ④⑤⑥⑦ ③ t | ② 巡＼ ① 夕 厂 a 第2 1 排炮孔布置图 a T h eh o l el a y o u to f t h el i n e2 1 1 6 4 4 1 6 4 3 1 6 3 3 ①②③ ，钰 ④ ⑤ b 第3 0 排炮孔布置图 b T h eh o l el a y o u to f t h e l i n e3 0 图9 炮孔布置图 F i g ．9 T h eb l a s th o l el a y o u tp l a n ∞ 嘶 舛 眈 m 仉 m m 一讲蛐一一目l A鲁∞皇∞。∞％时 置 ” 型元％键单拍 关 哪 瓣蘸骥瓣鬻溅瓣瓣鬻熏 万方数据 8 2 爆破 2 0 1 9 年3 月 冬瓜林3 0 8 1 采场1 6 3 3 分层1 2 ～1 6 线沿脉第 2 1 3 0 排炮孔，通过现场工业试验并完成爆破后， 爆破效果较好，矿石的爆破块度大小适中，大块率较 低，各排炮孔的爆破参数及爆破效果统计如表3 所 示，爆破后采场爆破效果如图1 0 所示。 表3 爆破参数及爆破效果统计表 T a b l e3D e m o l i t i o np a r a m e t e r sa n db l a s t i n ge f f e c t s 3 结论 图1 0 爆破效果图 F i g ．1 0 T h eb l a s t i n ge f f e c td i a g r a m 镇沅金矿矿岩条件变化较大，主要采用无底柱 分段崩落法进行采场回采，运用A N S Y S ／L S D Y N A 数 值模拟软件对扇形中深孔炮孔爆破过程进行数值模 拟，通过对最小抵抗线和孔底距进行模拟分析，根据 关键单元的应力曲线分析爆炸仿真过程中的V o nM i s e s 有效应力峰值是否达到岩石的动态抗拉强度，从 而确定爆破效果，通过对不同抵抗线和孔底距模拟效 果的对比分析，确定出了适合镇沅金矿中深孔爆破的 最优参数最小抵抗线为1 ．5m ，最大孔底距2 ．0m 。 根据数值模拟分析所确定的爆破参数在冬瓜林 3 0 8 1 采场1 6 3 3 分层1 2 1 6 线进行现场试验，其爆 破效果较好，大块率较低。通过运用A N S Y S ／L S D Y N A 数值模拟进行优化后的爆破参数，有效的降 低了采场回采的综合成本，提高了劳动生产率。 参考文献 R e f e r e n c e s [ 1 ] 柴修伟，张电吉，周文勇，等．厚大磷矿体井下中深孔 爆破参数优化研究[ J ] ．爆破，2 0 1 1 ，2 8 2 3 5 - 3 8 ． [ 1 ] C H A IX i u - w e i ，Z H A N GD i a n _ j i ，Z H O UW e n - y o n g ，e ta 1 ． R e s e a r c ho nb l a s tp a r a m e t e r so fm i d d e e ph o l eu n d e r - g r o u n db l a s t i n gi nt h i c ka n dl a r g ep h o s p h o f i t e [ J ] ．B i a s t i n g ，2 0 1 1 ，2 8 2 3 5 - 3 8 ． i nC h i n e s e [ 2 ] 杨才亮．单自由面无限岩体中深孔微差爆破数值模拟 [ J ] ．矿业研究与开发，2 0 1 0 ，3 0 4 9 3 - 9 6 ． [ 2 ] Y A N GC a i ．1 i a n g ．N u m e r i c a ls i m u l a t i o no fd e e ph o l ed i f - f e r e n t i a lb l a s t i n gi nas i n g l ef r e es u r f a c ei n f i n i t er o c km a s s [ J ] ．M i n i n gR e s e a r c ha n dD e v e l o p m e n t ，2 0 1 0 ，3 0 4 9 3 - 9 6 ． i nC h i n e s e [ 3 ] 任高峰，王威，冯海昀，等．融冠1 号矿中深孔爆破 参数优化研究[ J ] ．爆破，2 0 1 1 ，2 8 3 3 4 - 3 5 ，3 8 ． [ 3 ] R E NG a o - f e n g ，W A N GW e i ，F E N GH a i - y u n ，e ta 1 ．O p t i m i z a t i o ns t u a yo nm i d - l e n g t hh o l eb l a s tp a r a m e t e r si n R o n g g u a nN o ．1M i n e [ J ] B l a s t i n g ，2 0 1 1 ，2 8 3 3 4 ．3 5 ， 3 8 ． i nC h i n e s e [ 4 ]阎南．扇形中深孔爆破大块率与爆破参数的关系 [ J ] ．中国矿山工程，2 0 1 6 ，4 5 5 7 4 刀6 ，8 4 ． [ 4 ]Y A NN a n ．R e l a t i o n s h i pb e t w e e nb o u l d e ry i e l da n dp a r a m - e t e r si nt h eb l a s t i n gw i t hr a d i a lm e d i u m - l e n g t hh o l e [ J ] ． C h i n aM i n eE n g i n e e r i n g ，2 0 1 6 ，4 5 5 7 4 - 7 6 ，8 4 ． i n C h i n e s e [ 5 ] 潘淼昌，武永猛，史秀志，等．凡口铅锌矿扇形中深孔 爆破参数数值模拟研究[ J ] ．采矿技术，2 0 1 3 ，1 3 5 8 7 - 8 9 ，1 4 2 ． [ 5 ] P A NM i a o c h a n g ，W UY o n g m e n g ，S H IX i u - z h i ，e ta 1 ． N n u m e r i c a ls i m u l a t i o no fd e e ph o l eb l a s t i n gp a r a m e t e r s i n F a n K o ul e a da n dz i n cm i n e [ J ] ．M i m n gT e c h n o l o g y ， 2 0 1 3 ，1 3 5 8 7 - 8 9 ，1 4 2 ． i nC h i n e s e [ 6 ] 钱七虎，戚承志．岩石、岩体的动力强度与动力破坏准 则[ J ] ．同济大学学报 自然科学版 ，2 0 0 8 ，3 6 1 2 万方数据 第3 6 卷第1 期姜永恒，雷恒永，杨杰，等镇沅金矿中深孔爆破参数数值模拟研究 8 3 上接第6 2 页 [ 1 0 ] [ 1 0 ] [ 1 1 ] [ 1 2 ] [ 1 2 ] 李鹏，周 佳，李振．爆炸应力波在层状节理岩 体中传播规律及数值模拟[ J ] ．长江科学院院报， 2 0 1 8 ，3 5 5 9 7 - 1 0 2 ． L IP e n g ，Z H O UJ i a ，L IZ h e n ．B l a s t i n gs t r e s sw a v ei nl a y e r e dj o i n t e dr o c km a s s p r o p a g a t i o nl a wa n dn u m e r i c a l s i m u l a t i o n [ J ] ．J o u r n a lo fY a n g t z eR i v e rS c i e n t i f i cR e ． s e a r c hI n s t i t u t e ，2 0 1 8 ，3 5 5 9 7 - 1 0 2 ． i nC h i n e s e 周剑，张路青，胡瑞林，等．大型结构面产状影响下 应力波传播规律研究[ J ] ．岩石力学与工程学报， 2 0 1 1 ，3 0 4 7 6 9 - 7 8 0 ． Z H O UJ i a n ，Z H A N GL u - q i n g ，H UR u i l i n ，e ta 1 ．S t u d yo f r u l e so fs t r e s sw a v e sp r o p a g a t i o nu n d e rv a r i o u sa t t i t u d e s o fl a r g e - s c a l ef r a c t u r e s [ 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 ＆E n g i n e e r i n g ，2 0 1 l ，3 0 4 7 6 9 - 7 8 0 ． i nC h i n e s e 胡建华，高宏伟，薛小蒙，等．岩体结构面对爆炸冲击 工程响应的影响[ J ] ．中南大学学报 自然科学版 ， 2 0 1 7 6 1 5 6 0 ．1 5 6 7 ． H UJ i a n - h u a ，G A OH o n g w e i ，X U EX i a o - m e n g ，e ta 1 ．I n f l u e n c eo fs t r u c t u r a l p l a n e so nr e s p o n s e o f e x p l o s i o n s h o c ke n g i n e e r i n g [ J ] ．J o u r n a lo fC e n t r a lS o u t hU n i v e r s i - [ 1 3 ] [ 1 3 ] [ 1 4 ] [ 1 4 ] [ 1 5 ] [ 1 5 ] [ 1 6 ] t y ，2 0 1 7 ，4 8 6 1 5 6 0 - 1 5 6 7 ． i nC h i n e s e 马冲，詹红兵，姚文敏，等．爆破振动作用下含软弱 夹层边坡稳定性及安全判据[ J ] ．爆炸与冲击，2 0 1 8 ， 3 8 3 5 6 3 - 5 7 1 ． M AC h o n g ，Z H A NH o n g - b i n g ，Y A OW e n m i n ，e ta 1 ．S t a - b i l i t ya n ds a f e t yc r i t e r i o no fas l o p ew i t hw e a ki n t e r l a y e r u n d e rb l a s t i n gv i b r a t i o n [ J ] ．E x p l o s i o n ＆S h o c kW a v e s ， 2 0 1 8 ，3 8 3 5 6 3 - 5 7 1 ． i nC h i n e s e 张建华，黄刚．爆破震动对H M L 露天矿边坡影响 的研究[ J ] ．爆破，2 0 1 2 ，2 9 2 1 1 4 - 1 1 8 ． Z H A N GJ i a n - h u a ，H U A N GG a n g ．S t u d yo fb l a s t i n gv i - b r a t i o ne f f e c to nH M Lo p e n p i tm i n es l o p e [ J ] ．B l a s t i n g ， 2 9 2 1 1 4 一1 1 8 ． i nC h i n e s e 杨金华，夏元友，祝文化．风化岩体边坡开挖爆破试 验研究[ J ] ．爆破，2 0 1 2 ，2 9 1 5 1 - 5 3 ． Y A N GJ i n h u a 。X I AY u a n - y o u ．Z H UW e n - h u a ．E x p e r i m e n ts t u d yo f b l a s t i n ge x c a v a t i o no fw e a t h e r i n gr o c k s l o p e [ J ] ．B l a s t i n g ，2 0