煤层底板分层装药爆破数值模拟研究与应用.pdf

第3 2 卷第2 期 爆破 V 0 1 ．3 2N o - 2 2 0 1 5 年6 月B L A S T I N G J u n ．2 0 1 5 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 5 ．0 2 ．0 1 4 煤层底板分层装药爆破数值模拟研究与应用 朱宽，周桂松，吴欣欣，曹迸军 葛洲坝易普力股份有限公司，重庆4 0 1 1 2 1 摘要通过有限元分析软件A N S Y S ／L S D Y N A 对含3m 厚煤层底板的岩层分别建立5 0 0c m 5 0 0c m 1 4 0 0e m 大小的连续装药模型和分层装药模型，采用拉格朗日算法计算得到两种模型在不同位置的岩石应 力分布规律。对比两种模型的计算结果可知分层装药结构对于改善底部装药段岩石破碎效果不明显，但对 于增加中间充填段以上岩石应力峰值次数、延长有效应力作用时间具有明显作用。分层装药技术在汝箕沟 无烟煤矿羊齿采区煤层底板爆破应用实践表明分层装药结构对于提高炸药能量利用率、改善爆破效果具有 重要意义 关键词A N S Y S ／L S D Y N A ；煤层底板；分层装药；拉格朗日；爆破效果 中图分类号T D 8 2 4文献标识码A文章编号1 0 0 1 4 8 7 X 2 0 1 5 0 2 0 0 7 2 0 6 N u m e r i c a lS i m u l a t i o nR e s e a r c ha n dA p p l i c a t i o n o fL a y e r e dC h a r g ei nS e a mF l o o rB l a s t i n g Z H UK u a n ，Z H O UG u i s o n g ，W UX i n x i n ，C A OJ i n - j u n G e z h o u b aE x p l o s i v eC oL t d ，C h o n g q i n g4 0 1121 ，C h i n a A b s t r a c t B yt h ef i n i t ee l e m e n ta n a l y s i ss o f t w a r eA N S Y S ／L S D Y N A ，t h ec o n t i n u o u sc h a r g i n gm o d e la n dl a y e r e d c h a r g i n gm o d e lw i t h5 0 0c m 5 0 0e m 1 4 0 0c mw e r es e tu pi ns e a mf l o o rw i t h3 mt h i c k n e s s ．T h e s t r e s sd i s t r i b u - t i o nr e g u l a t i o n si nr o c kw e r ef o u n db yu s i n gL a g r a n g ec a l c u l a t i n gm e t h o d ．R e s u l t ss h o w e dt h a tl a y e r e dc h a r g i n gs t r u c t u r ep r o d u c e dl i t t l ee f f e c to ni m p r o v i n gf r a g m e n ts i z ea tt h eb o t t o mo fb l a s th o l e ，h o w e v e r ，w h i c hp l a y e dr o l e si ni n c r e a s i n gf r e q u e n c yo fr o c ks t r e s sp e a k sa n de x t e n d i n ga c t i n gt i m eo fe f f e c t i v er o c ks t r e s s ．W h e na p p l i e di nY a n g c h i a n t h r a c i t em i n i n ga r e aa tR uQ i g o ur o a lm i n e ，l a y e r e dc h a r g i n gm e t h o da l s op l a y e dv i t a lr o l ei ne n h a n c i n gt h eu t i l i z a t i o nr a t i oo fe x p l o s i v ea n di m p r o v i n gb l a s t i n ge f f e c t ． K e yw o r d s A N S Y S ／L S D Y N A ；s e a mf l o o r ；l a y e r e dc h a r g e ；l a g r a n g e ；b l a s t i n ge f f e c t 随着煤矿开采技术的发展以及大型采装挖运设 备的出现，井下开采的方式逐渐被露天开采的方式 所取代。露天开采通常采用中深孔台阶进行爆破， 孔深多为1 0m 以上，有的甚至高达2 0 ～3 0m 。因 台阶跨度较大，对岩性较复杂的矿区，很有可能会跨 越具有复合岩性的岩体。煤层底板主要为沉积岩， 厚度通常为2m 左右，因煤层的差异局部厚度可达 收稿日期2 0 1 5 0 3 一0 9 作者简介朱宽 1 9 8 8 一 ，男，硕士，2 0 1 4 年毕业于武汉科技大学 工程力学专业，现任葛洲坝易普力股份有限公司爆破技术 员， E - m a i l 4 1 7 4 9 8 8 8 3 q q ．e o m 。 3 ～4m ，其硬度和力学性能与其毗邻的岩层相差巨 大，采用连续装药的方式进行爆破时因底板层正好 位于炮孔堵塞段，爆破后底板部分易形成整体抬起。 增加单耗后底板以下部分产生过破碎，而底板层只 是局部产生微小裂纹，破碎效果改善不明显，增加了 挖运难度和二次破碎成本。研究表明，调整装药结 构能够调整应力场的分布从而有效改善爆破效果 朱红兵、卢文波等利用爆轰波理论分析了空气间隔 装药孔内一维不定长激波的相互作用及其在孔底与 堵头的反射过程o ，同时对孔内各点的压力随时间 万方数据 第3 2 卷第2 期 朱宽，周桂松，吴欣欣，等煤层底板分层装药爆破数值模拟研究与应用 7 3 变化过程进行分析，计算得到空气间隔装药在梯段 爆破工程中合理的空气层长度比例约为3 0 ％～ 4 2 ％。辜大志对采用孔底间隔装药结构的爆破效果 进行了试验研究“ 。，结果表明使用合理的孔底空气 间隔装药结构可以使大块降低率达到3 0 ％以上。 神宁集团汝箕沟无烟煤分公司羊齿采区位于宁 夏回族自治区贺兰山北部腹地，汝箕沟矿区中部，采 区盛产的“太西无烟煤”以其优良煤质深受国内外 用户的青睐，在国际市场上享有盛名。采区目前主 要爆破区域为 1 9 6 0 ～ 2 1 4 0 平台，使用1 0m 台 阶梯段爆破，2 0 1 4 年爆破工程量约9 0 0 万m 3 。矿区 岩石以半坚硬、坚硬砂岩为主， 1 9 8 0 平台出煤完 成后，遇到平均厚度约3 m 的煤层顶板，主要成分为 坚硬页岩，顶板以下为泥质砂岩。基于以上研究基 础，采用动力分析软件L S D Y N A 对该岩层使用分层 装药结构爆破的合理性进行了理论研究，并进行了 数次工程爆破实践。 1 ．2 爆轰波作用机理 分层药包和底部连续装药柱起爆后，爆轰波沿 药柱开始传播，与此同时产生的爆生气体开始膨胀， 当爆轰波达到堵塞段表面时将产生强烈的冲击压 缩，形成的冲击波将沿炮孔轴向继续向前传播。在 轴向压缩作用下，堵塞段将向径向膨胀，造成对孔壁 岩石的压缩作用。同时，两段药柱在爆炸过程中产 生的冲击波一部分将穿过填塞层而在孑L 底发生反射 和透射，在爆轰波相互叠加作用下促使岩石破碎。 另一方面，由于上层药包的起爆也为下层药柱的爆 破创造了新的自由面，使最小抵抗线减小，爆破的作 用范围增大。此外，两段药柱爆破后破碎的岩块会 在运动过程中发生碰撞，借助自身的动能进行再次 破碎。碰撞能量的消耗有利于减弱岩块的运动速度 和抛掷距离，使破碎下来的岩块更加均匀，大块率降 低，岩块分布更加集中。 1 分层装药作用机理 2 数值模拟分析 1 ．1 分层装药能量分布 采用连续耦合装药结构实施爆破时，炸药爆炸 后产生的压缩应力波和准静态压力直接作用于岩石 表面，瞬问压力高达5 0 0 0 ～1 00 0 0M P a ，在孑L 壁周围 将形成粉碎区。粉碎区不仅消耗掉大量能量，而且 温度会骤然升高，造成炸药能量的进一步损失。使 用分层装药技术其本质是以岩屑为间隔介质，改变 药柱与炮孔壁直接接触的关系，降低压缩应力波和 爆轰气体产物作用于孔壁的初始压力，延长压力的 作用时间，使炮孔周围不产生粉碎区或使粉碎区范 围减少，从而在减少炮孔装药量的前提下，使能量分 布更均匀，提高能量利用率，连续装药与分层装药结 构孔壁压力与作用时间关系如图1 所示。同时，采 用分层装药还能提高药柱上部的作用高度，使药柱 以上部分岩石破碎更加均匀，从而达到降低大块的 目的‘3 引。 钆 R 出 到 J h 时间t 图1孔壁压力与作用时间关系 F i g ．1 V a r i a t i o n so ft h ee x p l o s i o n p r e s s u r e P w i t ht i m e t 2 ．1 理论模型 在L S D Y N A 中采用c m - g - I x s 制分别建立连续 装药和分层装药爆破结构模型，模型大小均为 5 0 0c m 5 0 0c m 1 4 0 0c m ，上层建立厚度为3 0 0e m 煤层顶板，炮孔直径为1 4c m ，炮孔深度1 2 0 0c m ，堵 塞总长度均为3 ．5m ，分层装药部分从煤层顶板与泥 质砂岩交接面开始建立长度为1 ．0m 药柱，如图2 所 示。模型竖直方向对称，为节省资源，取1 ／4 模型进 行简化计算，计算中使用共节点方式计算，即划分网 格时炸药与岩石具有相同的节点，模型使用拉格朗日 算法可以追踪到炮孔膨胀过程。其中分层装药的方 式设两处起爆点，两种方式均采用孑L 底起爆方式，计 算时长取5 0 0 0I x s 。炸药和岩石选用实体单元S O L I D l 6 4 ，模型在X Z 面和泫面均设置为对称边界，上 表面为自由边界，其它3 个边界均使用无反射边界定 义约束，防止应力波在边界上发生反射”墙o 。 2 ．2 材料模型及参数 1 炸药及其状态方程 炸药采用高能炸药模型 M A T H I G H E X P L O S I V E B U R N ，用J W L 状态方程来描述爆炸过程爆 轰产物的压力、体积及能量特性一，1 0 。 P e 0 。 A 一燕 e 也～B ，一蒜 e 啦～等 1 万方数据 7 4爆破 2 0 1 5 年6 月 式中P 一为由J W L 状态方程决定的压力；V 为相对 体积；玩为初始比内能；A 、B 、R 。、R 和t O 为描述 J w L 方程的5 个独立物理常数，2 岩石乳化炸药的 参数取值见表1 。 o薹 煤层底板 N ， 糕 瑟 C n 一 雾 e 软弱砂岩 o ’寸 一 e n 卜 7 ● 霎 煤层底板 一 高 一 t o 雾 o软弱砂岩 o 甘 I n一 。c 一 7 图2 计算模型示意图 单位c m F i g ．2 T h em o d e lo fc a l c u l a t i o n u n i t e m 2 岩石材料模型 爆炸过程岩石达到屈服应力极限以致破碎，选 择随动硬化模型 M A T P L A S T I C K I N E M A T I C 来模 拟煤层顶板及软弱砂岩破坏过程，可通过屈服应力 极限和失效应变来控制材料失效。在屈服应力中引 进应变率因子1 e ／C 9 ，岩体屈服应力o r ，与应 变率善的关系如下⋯o 盯， [ 1 - 占C 9 ] o - 。 脾占≥ 2 其中 E p E o E 。／ E o E 。。。 3 式中d r 。为岩体的初始屈服应力；E 。为杨氏模量；0 为加载应变率；C 为C o w p e r ．S y m o n d s 应变率参数； E 。为岩石塑性硬化模量；E 。为切线模量；卢为各向 同性硬化和随动硬化贡献的硬化参数，0 ≤口≤1 ；， 为岩石塑性应变分量。堵塞段材料选用软弱砂岩材 料进行近似，岩石力学特性参数见表2 。 表1 炸药参数 T a b l e1 E x p l o s i v ep a r a m e t e r s 3 模拟结果及分析 分别选取连续装药和分层装药模型上的单元， 研究其最大主应力分布情况及随时间变化规律。以 分层装药模型为对象依次沿竖直方向在Z 6 0 0c m 平面、Z 1 0 5 0c m 平面、Z 1 1 5 0e m 平面取三组单 元，三组单元分别代表了在连续装药段、分层装药段 及两段药柱间充填段的岩石受力情况。每一组沿炮 孑L 径向分别在距离炮孔中心1 0 0e m 、2 0 0c m 、 3 0 0c m 、4 0 0c m 、5 0 0c m 位置取5 个单元，相应的在 连续装药模型的同样位置选取单元进行对比分析， 计算单元示意图如图3 所示。 Z 6 0 0e m 平面沿径向单元主应力时程曲线如 图4 、图5 所示 图4 为连续装药模型，图5 为分层 装药模型 ，图中A 、口、C 、D 、E 均表示距离炮孔中心 1m 、2 m 、3m 、4 m 、5m 位置的单元。对比两图可以 得到在应力峰值及曲线走势上两者基本一致在离 炮孑L 中心1m 处应力峰值为2 7 0M P a ，2m 处应力峰 值为1 3 0M P a ，3 ～5m 处应力峰值相差不大，均为 7 5M P a 左右，软弱砂岩主要受压应力破坏。结果表 明使用分层装药结构后在连续装药段主应力分布 规律变化不大。 图3 计算单元示意图 F i g ．3 T h ee l e m e n to fc a l c u l a t i o n 万方数据 第3 2 卷第2 期 朱宽，周桂松，吴欣欣，等煤层底板分层装药爆破数值模拟研究与应用 7 5 Z 1 0 5 0c m 平面沿径向单元主应力时程曲线 如图6 、图7 所示 图6 为连续装药模型，图7 为分 层装药模型 。对于连续装药模型在0 ～1m s 间由 于爆轰波还未传播到药柱顶端，故受力基本为零；对 于分层装药模型由于上、下两段药柱同时起爆，上部 药柱冲击波已率先对此部位发生作用。从应力峰值 分析，对于炮孔近区lI T I 处，连续装药模型应力峰值 为2 7 5M P a ，分层装药模型在1m 处的应力峰值虽 然只有1 3 0M P a ，但由于应力叠加作用在整段曲线 中出现了多个压应力峰值。从有效应力作用时间分 析，在1i n 处连续装药模型应力峰值超过1 0 0M P a 作用时间为0 ．6m s ，分层装药模型的作用时间为 1 ．2m s ，有效应力作用时间延长l 倍。结果表明分 层装药结构在两段药柱充填段对于增加应力峰值次 数、延长有效应力作用时间具有明显效果。 图4z 6 0 0c m 平面连续装药单元主应力时程曲线 F i g ．4E l e m e n tp r i n c i p a ls t r e s st i m eh i s t o r yc u r v eo f c o n t i n u o u sc h a r g i n gm o d e lo nZ 6 0 0c mp l a n a r 图5z 6 0 0e m 平面分层装药单元主应力时程曲线 F i g ．5 E l e m e n tp r i n c i p a ls t r e s st i m eh i s t o r yc u r v eo f l a y e r e dc h a r g i n gm o d e lo nZ 6 0 0c mp l a n a r 图6Z 1 0 5 0c m 平面连续装药单元主应力时程曲线 F i g ．6 E l e m e n tp r i n c i p a ls t r e s st i m eh i s t o r yc u r v eo f c o n t i n u o u sc h a r g i n gm o d e lo nZ 1 0 5 0c mp l a n a r 图7Z 1 0 5 0c m 平面分层装药单元主应力时程曲线 F i g ．7 E l e m e n tp r i n c i p a ls t r e s st i m eh i s t o r yc u r v eo f l a y e r e dc h a r g i n gm o d e lo nZ 1 0 5 0c mp l a n a r Z 11 5 0 3 1 1 1 平面沿径向单元主应力时程曲线 如图8 、图9 所示 图8 为连续装药模型，图9 为分 层装药模型 。以距离炮孑L11 1 1 处的单元为例分别 对两种模型的主应力时程曲线进行分析，在0 ～ 1 ．5m s ，由于应力波尚未传递到该部位，连续装药模 型基本不受力的作用，对于分层装药模型，在上部药 柱的作用下，首先受压应力作用，其后在自由面和两 种岩石交界面发生反射的过程中形成多个应力峰 值，最高拉峰值应力达到1 0 5M P a ，远远超过岩石的 抗拉屈服强度，故岩石将受拉破坏。在1 ．5 ～2m s 范围内，连续装药模型受到压应力作用，产生峰值压 应力5 0M P a ，但尚未达到岩石的压应力屈服极限， 不会对岩石造成破坏；而对于分层装药模型，在该区 间受底部药柱作用也产生一段压应力作用曲线，其 峰值压应力为2 5M P a 。在2 5m s 范围内，连续装 药模型由于在自由面发生反射产生反射稀疏波，其 峰值为5 5M P a ，岩石将达到拉应力屈服极限而发生 破坏，对于分层装药模型，由于受两段药柱的应力波 叠加作用，虽然其峰值拉应力有所下降，但是其有效 作用时间明显增加，有利于加强对上部岩石的破坏 和粉碎。 时间／I 上一I E 3 l 图8 z l1 5 0c m 平面连续装药单元主应力时程曲线 F i g ．8 E l e m e n tp r i n c i p a ls t r e s st i m eh i s t o r yc u r v eo f c o n t i n u o u sc h a r g i n gm o d e lo nZ 115 0c mp l a n a r 4 工程应用实例 汝箕沟无烟煤分公司羊齿采区2 0 1 4 年8 月2 6 日首次对 1 9 8 0 水平煤层底板采取连续装药结构 5 0 5 0 5 0 S 0 o 加一。之j．o coI邑芒錾粤ff型r丰嚼 万方数据 7 6爆破2 0 1 5 年6 月 实施爆破，平均孔深1 3 ．0m ，顶板厚2 3m ，爆破后 大块严重，最大块度直径达5m 以上。为降低大块 改善爆破效果，根据现场实际情况采取了分层装药 结构来实施爆破。设计具体如下第一步，根据岩石 性质、炸药种类等确定炸药单耗q ，k g ／m 3 ；第二步， 根据孔排距与台阶高度确定单孑L 装药量Q ，k g ，Q 口b H q ，式中n 为孑L 距，m ，b 为排距，m ，H 为台 阶高度，m ；第三步，根据延米装药量确定总充填长 度L ，m ，其中粉乳西1 4 0m m 炮孔延米装药量为 11 ．5k g ／m ，L Q ／11 ．5 ；第四步，确定各充填段长 度，为防止飞石并控制大块，上部充填长度￡，至少 为1 ．5 1 ．6m ，实际工程中取1 ．6m ，下部充填长度 L L L 。，m ；第五步，确定各装药段长度和药量， 令分层段的炸药单耗为q ，，k g ／m 3 ，要达到满意的破 碎效果且不产生飞石，则必须满足0 1 2 ] 1 ／3 q ≤q ．≤ 1 ／2 q ，按照q 取0 ．4 5k g ／m 3 计算，则分层段装药量 Q 。为6 ～8k g ，实际工程中采取3 支咖9 0m m 膏状 乳化炸药绑扎在一起装药 每支长3 0c m ，重2k ， 共小k g ，需要指出的是分层装药段起爆药包的位 置需要根据岩层走向和煤层底板厚度酌情处理，避 免将起爆药包置于煤层底板与下面软岩交界处。底 部使用粉状乳化炸药连续装药，装药量Q 为Q Q Q ．，k g 。分层装药结构图如图l O 所示。 2 0 1 4 年9 月～1 0 月葛洲坝易普力宁夏分公司 采用分层装药结构对 1 9 8 0 水平煤层底板共实施 了1 3 次爆破，共计1 9 ．8 万m 3 ，除编号为0 9 0 5 7 爆 破单耗选择过大，堵塞质量不佳导致冲孔现象较严 重，产生大块较多外，其它试验大块率均得到有效控 制，如表3 所示。针对前几次爆破过程中出现的后 排拉裂较严重，调整方案后将后排孔装药量适当减 小，后几次爆破对保留区的破坏得到有效控制，取得 了较好的爆破效果。 营1 5 士 暮l 0 景 呈0 ．5 R 蠹 o 堪一O5 图9Z 1 1 5 0e m 平面分层装药单元主应力时程曲线 F i g ．9 E l e m e n tp r i n c i p a ls t r e s st i m eh i s t o r yc u r v eo f l a y e r e dc h a r g i n gm o d e lo nZ 1 15 0e mp l a n a r 导爆管雷管 图1 0 分层装药结构图 F i g ．10T h ec h a r to fl a y e r e db l a s t i n g 表3 分层装药爆破参数 T a b l e3 L a y e r e db l a s t i n gp a r a m e t e r s 万方数据 第3 2 卷第2 期朱宽，周桂松，吴欣欣，等煤层底板分层装药爆破数值模拟研究与应用 7 7 5 结论 相比连续装药结构，采用分层装药爆破对于改 善底部岩石破碎效果不明显，但能有效提高上部岩 石的破碎效果。 采用分层装药后，相比连续装药结构炮孔近区 压力峰值略有降低，但应力峰值出现次数和有效应 力作用时间明显增加，有利于增加爆生气体作用时 间，增强对岩体的充分破碎。 将分层装药结构的理论分析与数值模拟结果运 用到爆破实践取得了较好的技术和经济效果，提高 了炸药利用率，降低了挖运和二次破碎成本。 [ 2 ] [ 2 ] [ 3 ] [ 3 ] [ 4 ] [ 4 ] 参考文献 R e f e r e n c e s 朱红兵，卢文波，吴亮．空气间隔装药爆破机理研究 [ J ] ．岩土力学，2 0 0 7 ，5 5 9 8 6 9 9 0 ． Z H UH o n g - b i n g 。L UW e n b o ．W UL i a n g ．R e s e a r c ho nm e c h - a n i s mo fa i r - d e c k i n gt e c h n i q u ei nb e n c hb l a s t i n g [ J ] ．R o c k a n dS o i lM e c h a n i c s ，2 0 0 7 ，5 5 9 8 6 - 9 9 0 ． i nC h i n e s e 辜大志，谢圣权，陈寿如．孑L 底空气间隔装药改善爆破 震动和效果的研究[ J ] ．采矿技术，2 0 0 4 ，4 4 6 4 ．6 6 ． G UD a z h i ，X I ES h e n g - q u a n ，C H E NS h o u 一1 3 1 ．A i ri n t e r v a l c h a r g ea tt h eb o t t o mo ft h eb o r et oi m p r o v eb l a s t i n gv i b r a t i o na n de f f e c tr e s e a r c h [ J ] ．M i n i n gT e c h n o l o g y ，2 0 0 4 ， 4 4 6 4 6 6 ． i nC h i n e s e 李晓杰，曲艳东，闫鸿浩，等．中深孔爆破分层装药分 层填塞研究[ J ] ．岩石力学与工程学报，2 0 0 6 ，2 2 5 3 2 6 9 3 2 7 5 ． L IX i a o - j i e ，Q UY a n d o n g ，Y A NH o n g h a o ，e ta 1 ．S t u d yo n t h el a y e r e dc h a r g i n ga n dl a y e r e dt a m p i n gf o rm e d i u ma n d d e e ph o l eb l a s t i n g [ 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 sa n dE n g i n e e r i n g ，2 0 0 6 ，2 2 5 3 2 6 9 3 2 7 5 ． i nC h i n e s e 璩世杰，陈广平，尚峰华，等．台阶爆破垂直中深孔间 隔装药技术的理论分析与应用[ J ] ．工程爆破，1 9 9 9 ， 5 4 5 7 6 1 ． Q US h i - j i e ，C H E NG u a n g p i n g ，S H A N GF e n g h u a ．e ta 1 ． T h e o r e t i c a la n a l y s i sa n df i e l dt e s t so ne f f e c to fd e c k e d c h a r g e si nv e r t i c a lh o l e so fb e n c hb l a s t i n g [ J ] ．E n g i n e e r - i n gB l a s t i n g ，1 9 9 9 ，5 4 5 7 6 1 ． i nC h i n e s e [ 5 ]时党勇，李裕春，张建民．基于A N S Y S ／I S D Y N A8 ．1 进行显式动力学分析[ M ] ．北京清华大学出版社， 2 0 0 4 ． [ 6 ]白金泽．L S D Y N A 3 D 理论基础与实例分析[ M ] ．北京 科学出版社，2 0 0 5 ． [ 7 ] 杨军，杨国梁，张光雄．建筑结构爆破拆除数值模拟 [ M ] ．北京科学出版社，2 0 1 2 ． [ 8 ] L S T C ．L S D Y N Ak e y w o r du s e r ’Sm a n u a l [ z ] ．C a l i f o r n i a L i v e r m o r eS o f t w a r eT e c h n o l o g yC o r p o r a t i o n ，2 0 0 3 ． [ 9 ]刘红岩，王新生．地表垂直爆破震动速度的数值计算 [ J ] ．爆炸与冲击，2 0 0 8 ，2 8 3 2 6 5 ．2 6 7 ． [ 9 ] L I UH o n g y a n ，W A N GX i n s h e n g ．N u m e r i c a lc a l c u l a t i o n o fg r o u n dv e r t i c a lb l a s t i n gv i b r a t i o nv e l o c i t y [ J ] ．E x p l o s i o na n dS h o c kW a v e s ，2 0 0 8 ，2 8 3 2 6 5 2 6 7 ． i nC h i n e s e [ 1 0 ] 夏祥，石永强，李海波，等．岩体单孔及群孔齐发爆 破爆炸荷载数值分析[ J ] ．岩石力学与工程学报， 2 0 0 7 ，2 6 S 1 3 3 9 0 3 3 9 6 ． [ 1 0 ] X I AX i a n g ，S H IY o n g - q i a n g ，L IH a i b o ，e ta 1 ．N u m e r i c a l a n a l y s i so fe x p l o s i v el o a do fs i n g l e h o l eb l a s t i n ga n d m u l t i p l e h o l es i m u l t a n e o u sb l a s t i n g i nr o e k m a s 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 sa n dE n g i n e e r i n g ， 2 0 0 7 ，2 6 S 1 3 3 9 0 3 3 9 6 ． i nC h i n e s e [ 1 1 ] 殷秀红，钟冬望，黄小武，等．水下钻孔爆破岩石各区 应力的数值模拟[ J ] ．爆破，2 0 1 4 ，1 2 4 1 1 - 1 5 ． [ 1 1 ] Y I NX i u h o n g ，Z H O N GD o n g - w a n g ，H U A N GX i a o W U ，e t a 1 ．N u m e r i c a ls i m u l a t i o no fr o c ks t r e s so fu n d e r w a t e r d r i l l i n gb l a s t i n g [ J ] ．B l a s t i n g ，2 0 1 4 ，1 2 4 1 1 1 5 ． i n C h i n e s e [ 1 2 ]邓端正．新型间隔装药技术降低露天矿深孔爆破炸 药单耗的应用[ J ] ．生产技术，2 0 0 9 2 3 9 4 0 ． [ 1 2 ] D E N GD u a n z h e n g ．T h eA p p l i c a t i o no f ”n e w ”i n t e r v a l c h a r g et e c h n o l o g yt or e d u c et h ec o n s u m p t i o no fe x p l o s i v e sa b o u td e e p h o l eb l a s t i n gi no p e n ．p i tm i n e [ J ] ．c e m e n tE n g i n e e r i n g ，2 0 0 9 2 3 9 4 0 ． i nC h i n e s e 万方数据