露天采场高陡边坡下深部采空区的爆破处理研究_解治宇.pdf

收稿日期2019-12-17 基金项目 “十三五” 国家重点研发计划项目 （编号 2016YFC0801603） 。 作者简介解治宇 （1970） ， 男， 高级工程师。 总第 523 期 2020 年第 1 期 金属矿山 METAL MINE 露天采场高陡边坡下深部采空区的爆破处理研究 解治宇 1 李翰林 2 房洪亮 1 韩连生 1 张德辉 1 沙成满 31 （1. 鞍钢集团矿业有限公司， 辽宁 鞍山 114001； 2. 辽宁科技大学矿业工程学院， 辽宁 鞍山 114051； 3. 东北大学资源与土木工程学院， 辽宁 沈阳110819） 摘要鞍本地区是我国最大的BIF型铁矿集中分布区， 其铁矿资源量/储量约占全国总量的1/4， 但由于日伪 时期的掠夺式开采和矿业整合开发前的无序开发利用， 在一定深度空间内遗留了大量隐伏采空区， 为后期实现规 模化露天开采带来极大的安全隐患， 成为制约露天铁矿资源安全高效开采的难题。以弓长岭露天铁矿大砬子采区 高陡边坡下的深部采空区为研究对象， 利用LS_DYNA软件进行采空区的安全爆破模拟研究， 结合现场深部采空区 爆破参数试验结果， 提出了采用孔中起爆处理采空区的方法， 具体参数设置为空区预留顶板厚度3 m， 孔间延期时 间42 ms， 排间延期时间75 ms， 平均炸药单耗0.237 kg/t。在上述分析的基础上， 针对大砬子露天采场高陡边坡内的 深部不规则多层采空区提出了分层、 分区安全爆破处理的技术方案， 且空区边界用切割爆破方法处理， 爆破处理效 果良好， 采空区得到了安全有效治理。 关键词隐伏采空区露天铁矿爆破处理爆破模拟 中图分类号TD85， TD235文献标志码A文章编号1001-1250 （2020） -01-1049-07 DOI10.19614/ki.jsks.202001006 Study on Blasting Treatment of Deep Goaf under High and Steep Slope in Open-pit Mining Xie Zhiyu1Li Hanlin2Fang Hongliang1Han Liansheng1Zhang Dehui1Sha Chengman22 （1. Ansteel Mining Co.， Ltd.， Anshan 114001， China； 2.School of Mining Engineering， University of Science and Technology Liaoning， Anshan 114051， China； 3.School of Resources and Civil Engineering， Northeastern University， Shenyang 110819， China） AbstractThe largest concentrated distribution area of BIF-type iron deposits is Anben area in China， and its iron ore resources reserves account for about 1/4 of the countrys total.The development and utilization left a large number of concealed goafs in a certain depth space because of the predatory mining during the Japanese and puppet periods and the disordered ex- ploitation and utilization before the integrated exploitation of the mining industry，which brought great safety hazards to the large-scale open-pit mining in the later period， and became a problem restricting the safe and efficient mining of open-pit iron ore resources.In order to solve the above problems， taking the deep goaf under the high and steep slope in the Dalazi mining ar- ea of the Gongchangling Open-pit Iron Mine as an example， and using LS_DYNA software to simulate the safe blasting of the mined-out area.Based on the simulaiton results and the test results of blasting parameters in deep goaf， the safety treatment for explosion in the hole is proposed.The specific setting parameters are the reserved roof thickness of 3 m in the goaf， the delay time between the holes is 42 ms， the delay time between the rows is 75 ms， and the average unit consumption of ex- plosives is 0.237 kg/t.Based on the above discussion results， a layered， zoned and safe blasting technical solution is proposed for the deep irregular multi-layer goafs in the high and steep slopes of the Dalazi open-pit stope.The boundary of the goaf is treated with cutting blasting， and the blasting treatment effect is ideal， therefore， the goaf has been effectively and safety man- aged. KeywordsConcealed goaf， Open-pit iron mine， Blasting treatment， Blasting simulation Series No. 523 January2020 根据初步统计， 截至2015年底， 中国有色金属非 金属地下矿山共有采空区12.8 亿m3， 且地下矿山采 空区塌陷引发的安全事故接连发生， 造成了重大人 员伤亡和财产损失 ［1-3］。鞍本地区作为我国最大的 金属矿采空区致灾机理与安全高效处理技术 49 ChaoXing 金属矿山总第523期 BIF型铁矿集中区， 其铁矿资源量/储量占全国总量的 1/4， 由于日伪时期的掠夺式开采和矿业整合开发前 的无序开发和偷采滥挖， 在一定深度空间遗留了大 量隐伏采空区。上述采空区的存在给铁矿山规模化 的露天开采带来极大的生产安全巨大隐患， 制约了 露天采场资源安全高效开采。 上述现象在富铁矿发育的弓长岭地区特别典 型， 尤以其中的大砬子采区最为突出。露天铁矿大 砬子采区位于弓长岭二矿区东南端， 是一个开采多 年的露天老采区， 处于生产的中老年期。目前， 在大 砬子采区扩帮开采的南帮280 m标高采矿平台高陡 边坡下探测发现了大量隐伏采空区， 呈不规则多层 状分布， 且已探明的采空区范围大、 跨度广和形态复 杂， 并已经严重影响了矿山安全生产。本研究以大 砬子采区多层采空区为研究对象， 利用LS_DYNA软 件建立爆破处理的数值计算模型， 模拟在爆破荷载 作用下采空区的塌陷过程 ［4-8］， 并对模拟结果进行分 析， 在此基础上提出合理有效的采空区治理方案， 开 展多层采空区现场爆破处理研究 ［9-13］， 为类似矿山采 空区安全处理提供理论依据和实践参考。 1研究区地质概况 研究区处于弓长岭铁矿二矿区端部， 构造上位 于区域构造弓长岭背斜北翼， 矿区西北端以寒岭断 裂为界， 东南到地质勘探线30号剖面的大砬子， 全长 4 850 m， 面积4.01 km2， 其中大砬子采区 （24～30号勘 探线之间） 即为研究区 （图1） 。弓长岭铁矿带构造发 育， 褶皱、 断层齐全， 它们控制着矿床的形成与分 布 ［14］。铁矿带整体呈NW走向， 长约12 km， 但由于受 寒岭断裂、 偏岭断裂等一系列近平行的NE向断层的 影响， 铁矿带被切割分为一矿区、 二矿区、 三矿区、 老 岭八盘岭矿区 （包括老弓长岭、 独木、 哑叭岭和八盘 岭） 。弓长岭铁矿带整体为一NW向反S型褶皱， 但 NE向断层将NW向铁矿带错断成为断陷区和断隆区。 断陷区由于下降而有利于含铁岩系保存， 断隆区由于 上升遭受剥蚀强烈而不利于含铁岩系保存 ［14］。 弓长岭二矿区地层属于茨沟组， 位于寒岭断裂 与老岭断裂之间， 处于断陷区。二矿区含铁岩系呈 单层状残留体产于大片混合花岗岩中， 岩层走向 N20～60W， 倾向 NE， 倾角 60～85。含铁岩系两 端受断层影响产状有些变化， 其北西端受寒岭断裂 影响， 地层走向变为N80W， 倾向SE； 在其东南端因 受老岭断裂影响， 产状变为SN走向， 向东陡倾斜 ［14］。 2露天采场采空区爆破处理数值模拟研究 2. 1采空区概况 本研究爆破处理采空区位于大砬子采场南 帮280 m标高采矿平台下方40～70 m， 采空区地面投 影面积约5 000 m2， 按采空区平均厚度30 m计算得最 小体积为15万m3， 影响安全作业面积达上万平方米 （图1） 。本研究通过钻孔式三维激光扫描探测获取 采空区的三维点云数据 ［15-17］， 而后采用地学建模技术 形成采空区的三维虚拟实体。探测成果显示， 该平 台的采空区主要呈不规则的多层状分布特点 （图2） 。 2020年第1期 50 ChaoXing 解治宇等 露天采场高陡边坡下深部采空区的爆破处理研究 2. 2采空区爆破数值模拟及分析 2. 2. 1采空区的数值计算模型 根据采空区平面分布图 （图1） ， 作一条经过采空 区的剖切线A - A（图3） ， 由其切出的剖面可通过软 件LS_DYNA开展爆破数值模拟计算 ［18-21］。由剖面图 所切的空区剖面建立的二维理论计算模型如图4所 示。在模型左、 右、 下3个面上施加无反射边界条件 以模拟无限岩体， 并在底面边界施加X、Y、Z方向的 位移约束， 模型上表面均为自由边界， 对整体模型施 加-Y方向的重力， 重力加速度为9.8 m/s2。岩体和填 塞材料选用HJC模型， 采用高能炸药材料和状态方 程*EOS_JWL模拟乳化炸药， 选用NULL模型和线性 多项式模拟空气材料。岩石和填塞材料采用 La- grange算法， 炸药和空气采用ALE算法。起爆时， 采 用孔中起爆的方式， 预留空区顶板厚度为3 m， 采用 排距7 m， 为节约计算成本， 排间延期时间设置为17 ms。 2. 2. 2数值模拟结果及有效应力分析 由不同时刻采空区在爆破荷载作用下的有效应 力云图 （图5） 可知， 炸药起爆后爆破应力波对炮孔右 侧和下侧岩体的作用较强， 这两部分岩体所受到的 有效应力较大， 原因是炮孔右侧和下侧方向的抵抗 线最短， 爆轰波优先从这两个自由面溢出， 有利于达 到理想的工程爆破效果， 即促使采空区塌陷和台阶 发生破坏。 应力波遇到自由面发生反射， 形成反射拉伸波， 产生反向拉应力， 自由面处的岩体发生片落。前排 炮孔能为后排炮孔创造出新的自由面， 应力波在3个 自由面之间来回反射、 叠加， 岩体在复杂的多重因素 作用下发生破坏。右侧台阶坡面的破坏主要是第一 排炮孔炸药爆炸作用的结果， 后排炮孔也会对坡面 产生影响， 但主要是使破碎后的岩体产生侧向位移。 随着后排炮孔对采空区顶板的持续作用， 空区顶板 应力集中的区域越来越多， 顶板受到的最大应力为 124.6 MPa， 超过了岩体的屈服强度。 起爆后 50 ms 时， 也就是第 3 个炮孔起爆后 16 ms， 爆轰波对上部边坡的作用范围最广， 这一时刻应 作为判断边坡是否稳定的最不利时刻。爆炸进行到 70 ms时， 第一排炮孔下方的顶板出现了贯通裂隙。 随着时间的推移， 采空区左上角出现了较严重的应 力集中， 上部边坡的尖角部位也出现了小范围的应 力集中现象， 在进行现场爆破时对其应引起足够重 视。由图5 （h） 可知 爆炸进行到0.1 s时， 应力集中的 部位在炮孔中心连线处、 炮孔底部和空区的左上角， 应力集中部位的最大应力达到125.7 MPa， 远大于岩 体的抗拉强度， 因此这部分岩体发生屈服。在后续 采矿作业过程中， 应特别注意空区左上角出现的应 力集中现象， 做好监测和防护措施， 必要时可进行二 次爆破， 消除隐患。 特别是起爆后50 ms是爆破应力波对边坡底部 作用范围最大的时刻， 故本研究据此来分析爆破荷 载作用对上部边坡稳定性的影响。由图6可知 浅色 区域岩体的有效应力值为47.41 MPa， 超过了岩体的 动态抗拉强度， 岩体处于屈服状态， 深色区域的有效 应力值小于41.92 MPa， 岩体不会发生屈服。边坡底 部岩体的屈服范围仅限于比空区边界大2 m的范围 内， 该范围内部分岩体达到了屈服强度。因此边坡 在爆破荷载作用下整体处于稳定状态， 不会发生边 2020年第1期 51 ChaoXing 金属矿山总第523期 坡失稳的次生灾害。为了更准确地了解各个部位岩 体的受力情况， 取3个部位的代表单元作为监测点， 这3个部位分别为 空区预留顶板中、 炮孔之间的部 分， 代表空区顶板应力薄弱部位； 炮孔填塞部分之间 的部位， 代表易出现大块岩石的部位； 空区左边界以 上部位， 代表伸入坡面的采空区顶板。监测点位置 和单元编号如图7所示。 图7所示的10个测点监测结果分别如图8所示。 由图8 （a） 可知 预留顶板处的岩体所受到的有 效应力远大于岩体本身的抗拉强度， 403962点的有 效应力在起爆后 60 ms时已经达到 71.31 MPa， 岩体 处于屈服状态。而其他两点 （400024 点和 433361 点） ， 所受到的有效应力达到80 MPa， 并且还在不断 增大， 这说明在爆破荷载作用下， 空区顶板的破坏过 程既不是一蹴而就的达到破坏状态， 也不是单因素 作用产生的破坏， 顶板岩体首先受到爆轰波和反射 拉伸波的巨大破坏， 顶板不会完全塌陷， 随后岩体在 叠加应力波、 爆生气体、 碰撞以及重力等综合作用下 发生破碎和塌陷， 岩体的塌落过程是以秒来衡量的。 由图8 （b） 可知 对比空区顶板预留部分的应力， 位于 炮孔顶部之间的岩体所受到的有效应力较小， 421793点的应力峰值为66.51 MPa， 小于预留部分应 力最小的点349250点 （71.31 MPa） ， 是因为对于炮孔 炸药来说， 炮孔底部与采空区自由面的距离小于炮 孔填塞， 因此炸药对空区顶板预留部分的作用更强， 2020年第1期 52 ChaoXing 解治宇等 露天采场高陡边坡下深部采空区的爆破处理研究 有利于空区得到有效处理， 在采空区实际工程处理 中， 炮孔填塞可以适当加长。第一排和第二排炮孔 填塞之间的岩体， 在起爆后极短时间内的有效应力 随即达到了峰值37.50 MPa， 之后在70 ms之前， 都维 持在较高的水平， 然后下降， 这部分岩体达到了屈服 强度。第二排和第三排炮孔填塞之间的岩体， 在起 爆后37.6 ms时达到了有效应力峰值67 MPa， 处于屈 服状态， 但之后有效应力迅速下降， 100 ms时已经下 降到10 MPa， 低于岩体的屈服强度， 这部分岩体在爆 破后容易出现大块岩石， 塌落后随即填充采空区。 由图8 （c） 可知 当采空区的形状存在尖角或直角时， 在角的部位会产生应力集中现象， 有效应力峰值为 66.81 MPa， 范围以交点为中心， 半径3 m以内。其他 各点距离最后一排炮孔11.6 m且位于坡底， 有效应力 变化趋势一致， 大小相近， 起爆后40 ms时达到了峰 值54.98 MPa， 岩体发生屈服， 说明采空区得到了有效 处理。 基于上述理论数值模拟研究， 可以得出采空区 处理的最佳爆破参数， 即空区预留顶板厚度3 m， 孔 间延期时间42 ms， 排间延期时间75 ms。 3露天采场采空区现场爆破处理实例 3. 1采空区爆破崩落处理原则 针对大砬子采场采空区空间分布特点和顶板厚 度变化规律 （图1、 图2） ， 并依据上述采空区爆破处理 数值模拟研究成果， 提出了分区、 分层处理多层不规 则采空区的原则 （图9） 。首先， 针对顶板厚度≤41 m 的区域， 如图9 （a） 中I分区， 可采用一次爆破诱导空 区冒落， 爆破后产生的松散岩石将顺着空区底板的 坡面充填部分剩余空区， 并可通过钻探验证充填程 度， 而后研究进一步处理方案。其次， 顶板厚度为 41～72 m的区域， 如图9 （a） 中阴影区域II和III分区， 实施分层崩塌处理。第一分层高度15 m， 待该阶段 处理形成260 m标高水平后， 剩余部分则与IV分区 同水平进行穿孔爆破处理， 可满足一次性崩塌处理 条件。同时， 为确保处理效果， 空区边界用切割爆破 处理。具体实施时， 切割孔穿孔深度按照比空区顶 板厚度少3 m进行设计， 孔位在空区内比空区边界收 缩 5 m， 并采用耦合装药， 装药高度按照 20 m设计。 切割孔内侧按照三角形布孔方式布置辅助孔， 孔网 参数设计为7 m7 m， 穿孔深度比切割孔浅18 m， 最 大限度的改善空区顶板破碎程度。特别值得一提的 是， 每次安全爆破处理结束后， 都应对未爆破空区区 域进行补充勘探， 以确定未爆破空区顶板变化在可 控范围内， 后续再对下次爆破的采空区区域实施穿 孔爆破作业。 3. 2采空区安全爆破处理 本研究以IV分区为例进行采空区安全爆破处理 流程分析。 （1） 爆破处理孔网参数。采空区爆破采用孔径 为250 mm的牙轮钻机和孔径为200 mm的潜孔钻机 联合钻进， 并根据岩石的可爆性以及数值模拟结果， 爆破孔网参数取7 m7 m， 并按照矩形孔网或者三 角形孔网布置。炮孔深度按照设计的预留空区顶板 厚度为3 m， 最大孔深达35 m。对于已经穿透顶板的 炮孔， 装药前将气体间隔器放置在顶板底部上方2 m 处或吊袋固定在顶板底部。炮孔填塞长度一般不低 于孔径的25倍， 按照5～7 m设计， 对于炮孔深度在30 m 以上的超大深度炮孔， 可取填塞长度为 7～10 m。 针对IV分区， 涉及到的爆破处理采空区面积约2 500 m2， 设计炮孔64个。 2020年第1期 53 ChaoXing 金属矿山总第523期 （2） 爆破处理炸药参数。采空区爆破处理主要 采用乳化铵油炸药， 炸药的爆轰性能经过具体现场 试验确定。在现场试验时， 选取合理的装药长度和 填塞长度， 并在孔底放置两发 400 g 起爆具。起爆 后如监测到振动数据， 则说明炸药可完全爆轰； 如 果只观察到水柱喷出， 而无振感， 则说明炸药处于 爆燃状态， 并没有完全爆轰。通过现场试验得出每 个炮孔平均需要装填乳化铵油炸药 400 kg， 总药量 为 25 600 kg， 预计爆破量可达108 000 t， 延米爆破量 为 58.8 t/m， 总米道为 960 m， 平均炸药单耗为 0.237 kg/t。 （3） 爆破处理的起爆方式。采空区爆破处理均 采用高精度毫秒导爆管起爆方式， 具体为孔中起爆， 每个炮孔1发500 g起爆弹和2发澳瑞凯高精度导爆 管， 爆区外连接50 m导爆管。同时， 根据上述数值模 拟结果和现场处理条件， 孔间延期时间取42 ms， 排 间延期时间取 75 ms。在炮孔没有穿透空区顶板的 前提下， 如果孔内无水， 则优先选装起爆能量较大的 乳化铵油炸药。具体起爆顺序按照正常台阶爆破方 式进行， 临近台阶自由面的炮孔先起爆， 后逐孔起 爆， 以确保空区顶板彻底坍塌， 大块率控制在1以 内， 同时不留根底， 控制爆破震动。 通过上述采空区爆破处理流程和参数设计， 对 大砬子采区高陡边坡内的不规则多层采空区进行了 安全处理， 实现了采空区处理的数值模拟计算与现 场安全处理有机结合， 保障了矿山安全生产。同时， 在安全处理采空区的前提下释放了部分铁矿资源， 取得了经济效益与安全效益的双重收益。 4结论 （1） 采空区爆破模拟成果显示起爆后 50 ms， 即 第3个炮孔起爆后16 ms， 爆轰波对上部边坡的作用 范围最广， 并随着时间的推移在采空区左上角出现 了较严重的应力集中， 爆炸进行到0.1 s时， 应力集中 在炮孔中心连线处、 炮孔底部和空区的左上角， 其最 大应力可达到125.7 MPa， 远大于岩体的抗拉强度， 可 使部分岩体发生屈服。 （2） 提出了分区、 分层处理深部多层采空区的方 案， 对于顶板厚度小于41 m的区域可一次爆破诱导 空区冒落， 且爆破后产生的松散岩石将充填剩余部 分空区。对于顶板厚度为41～72 m的采空区区域， 可 实施分层崩塌处理， 并需在空区边界进行切割爆破 处理。 （3） 提出了合理可行的大砬子多层采空区爆破 处理参数， 即爆破处理孔径250 mm， 孔网7 m7 m， 孔中起爆方式， 空区预留顶板厚度3 m， 孔间延期时 间 42 ms， 排间延期时间 75 ms， 平均炸药单耗 0.237 kg/t， 并在采空区进行了爆破实施， 使得采空区得到 了安全有效治理。 参 考 文 献 马海涛， 刘宁武， 王云海， 等.金属矿山采空区灾害防治技术研 究综述 ［J］ .中国安全生产科学技术， 2014， 10 （10） 75-80. Ma Haitao， Liu Ningwu， Wang Yunhai， et al.Review on research status of controlling techniques for goaf disaster in metal mine ［J］ . Journal of Safty Science and Technology， 2014， 10 （10） 75-80. 李夕兵， 李地元， 赵国彦， 等.金属矿地下采空区探测、 处理与安 全评判 ［J］ .采矿与安全工程学报， 2006， 23 （1） 24-29. Li Xibing， Li Diyuan， Zhao Guoyan， et al.Detecting， disposal and safty ution of the underground goaf in metal mines ［J］ .Journal of Mining Safety Engineering， 2006， 23 （1） 24-29. 王海君， 李克民， 陈树召， 等.地下转露天开采的采空区危害与 ［1］ ［2］ ［3］ 2020年第1期 54 ChaoXing 探测技术探讨 ［J］ .煤炭工程， 2010 （10） 69-71. Wang Haijun， Li Kemin， Chen Shuzhao， et al.Disscussion on goaf danger and detection technology for mining converted from under- ground to surface ［J］ .Coal Engineering， 2010 （10） 69-71. 秦国震.中深孔爆破诱导冒落处理采空区技术研究 ［D］ .唐山 华北理工大学， 2017. Qin Guozhen.Study on Induction Caving Technique of Disposing of Goaf by Medium Deep Hole Blasting［D］ .TangshanNorth China University of Science and Technology， 2017. Zang L， Jia C， Zhang S， et al.Application on the open-air deep hole blasting caving in the underground mined-out area ［J］ .Mod- ern Mining， 2014 （4） 10-14. 贾宝珊 .露天开采中深孔爆破与地下采空区处理一体化实践 ［J］ .金属矿山， 2013 （2） 47-49 Jia Baoshan.The integrated practice of medium-depth hole blasting in open-pit mining and the treatment of underground goaf ［J］ .Metal Mine， 2013 （2） 47-49. Wang F， Tu S， Yuan Y， et al.Deep-hole pre-split blasting mecha- nism and its application for controlled roof caving in shallow depth seams ［J］ .International Journal of Rock Mechanics Mining Sci- ences， 2013， 64 （6） 112-121. 白俊， 王彦飞.大型叠层分布采空区爆破处理技术 ［J］ .金属矿 山， 2015 （9） 1-4. Bai Jun， Wang Yanfei.Blasting disposal technology in large layer- distributed mined-out area ［J］ .Metal Mine， 2015 （9） 1-4. 林卫星.井下竖向深孔逐孔等值微差爆破技术研究及应用 ［J］ . 矿业研究与开发， 2016， 36 （10） 35-38. Lin Weixing.Research and application on sequenced and equiva- lent millisecond blasting technology of vertical long hole in under- ground mine ［J］ .Mining Research and Development， 2016， 36 （10） 35-38. 胡从倩， 王德胜， 崔帅立， 等.多层采空区稳定性分析及处理措 施 ［J］ .金属矿山， 2013 （11） 29-33. Hu Congqian， Wang Desheng， Cui Shuaili， et al.Stability analysis and treatment of poly-laminate goaf ［J］ .Metal Mine， 2013 （11） 29- 33. 贾宝珊， 闫伟峰.露天正常台阶深孔爆破处理地下采空区的实 践 ［J］ .爆破， 2012， 29 （4） 65-69. Jia Baoshan， Yan Weifeng.Application of bench deep hole blasting in disposing of underground goaf ［J］ .Blasting， 2012， 29 （4） 65-69. Xia H， Yan B.Application of deep-hole pre-cracking blasting in the preserved roadway beside the gob for working Gaolin in coal seams ［J］ .Energy Procedia， 2012， 16828-835. Liu X L， Luo K B， Li X B， et al.Cap rock blast caving of cavity un- der open pit bench ［J］ .Transactions of Nonferrous Metals Society of China， 2017， 27 （3） 648-655. 王恩德， 夏建明， 赵纯福， 等.弓长岭铁矿床磁铁富矿形成机制 探讨 ［J］ .地质学报， 2012， 86 （11） 1761-1772. Wang Ende， Xia Jianming， Zhao Chunfu， et al.ing mechanism of high-grade magnetite bodies in Gongchangling， Liaoning Prov- ince ［J］ .Acta Geologica Sinica， 2012， 86 （11） 1761-1772. 贾三石， 付建飞， 门业凯， 等.深凹露天铁矿隐伏空区三维激光 探测技术方法应用研究 ［J］ . 安全与环境学报， 2019， 19 （5） 1581-1586. Jia Sanshi， Fu Jianfei， Men Yekai， et al. 3-D laser scanning ap- proach to the concealed goaf under an open-pit iron mine ［J］ .Jour- nal of Safety and Environment， 2019， 19 （5） 1581-1586. Miller F， Potvin Y， Jacob D.Laser measurement of open stope dilu- tion ［J］ .CIM Bulletin， 1992， 8596-102. Scheding S， Dissanayake G， Nebot E M， et al.An experiment in au- tonomous navigation of an underground mining vehicle［J］ . IEEE Transactions on Robotics and Automation， 2001， 15 （1） 85-95. Yan Y F， Zhang T， Cheng L F， et al.The simulation on millisecond blasting of submarine hard rock using LS-DYNA dynamic analysis ［J］ .Journal of Computational and Theoretical Nanoscience， 2013， 19 （6） 1567-1571. Zhang Y B， Chen C， Zhang Y P.Numerical simulation study on hy- draulic blasting based on LS-DYNA［J］ . Applied Mechanics and Materials， 20123245-3249. 李志忠 . 基于 LS-DYNA 的隧道爆破模拟分析 ［J］ . 公路， 2017 （8） 315-320. Li Zhizhong.Simulation analysis of tunnel blasting based on LS-DY- NA ［J］ .Highway， 2017 （8） 315-320. 叶海旺， 薛江波， 房泽法.基于LS-DYNA的砖烟囱爆破拆除模拟 研究 ［J］ .爆破， 2008， 25 （2） 39-42. Ye Haiwang， Xue Jiangbo， Fang Zefa.Numerical simulation of brick chimney blasting demolition based on LS-DYNA［J］ . Blasting， 2008， 25 （2） 39-42. ［4］ ［5］ ［6］ ［7］ ［8］ ［9］ ［10］ ［11］ ［12］ ［13］ ［14］ ［15］ ［16］ ［17］ ［18］ ［19］ ［20］ ［21］ 解治宇等 露天采场高陡边坡下深部采空区的爆破处理研究2020年第1期 55 ChaoXing