露天转地下开采不同坡高影响下覆岩采动响应特征研究_罗浪.pdf

收稿日期2020-01-14 基金项目国家自然科学基金青年基金项目 （编号 41702327） ， 江西省研究生创新专项资金项目 （编号 YC2018-S316） 。 作者简介罗浪 （1994） ， 男， 硕士研究生。 露天转地下开采不同坡高影响下覆岩采动 响应特征研究 罗浪 1， 2 李小双 1， 2 王运敏 31 （1. 江西理工大学资源与环境工程学院， 江西 赣州 341000； 2. 江西省矿业工程重点实验室， 江西 赣州 341000； 3. 中钢集团马鞍山矿山研究总院股份有限公司， 安徽 马鞍山 243000） 摘要以云南滇池周边典型的缓倾斜含软弱夹层中厚磷矿床露天转房柱法地下开采为研究背景， 通过现场 勘查取样和室内岩石力学试验， 基于AutoCAD、 ANSYS和FLAC3D等软件进行耦合建模， 对倾角为45时边坡高度为 50 m和80 m的上覆岩层位移与应力响应特征进行动态模拟和对比分析。研究结果表明 在露天终了边坡为50 m 时， 地下开采对覆岩影响相较80 m坡高小， 此时围岩相对稳定， 塑性区面积较小， 整个岩体处于相对稳定状态。露 天终了边坡为80 m时， 开采对覆岩的扰动较大， 回采完毕后矿柱处于临界破坏状态， 塑性区较50 m坡高进一步发 育并延伸至坡面， 不利于矿山的安全稳定。研究成果确定了该矿山露天终了边坡的合理高度范围， 对矿山的安全、 经济生产具有重要的理论指导意义。 关键词露天转地下房柱法覆岩采动响应数值模拟 中图分类号TD325文献标志码A文章编号1001-1250 （2020） -08-032-06 DOI10.19614/ki.jsks.202008006 Research on Mining Response Character of Underground Overlying Strata under Different Slope Heights during Transition from Open Pit to Underground Mining Luo Lang1， 2Li Xiaoshuang1， 2Wang Yunmin32 （1. Faculty of Resources and Environment Engineering， Jiangxi University of Science and Technology， Ganzhou 341000， China； 2. Jiangxi Key Laboratory of Mining Engineering， Ganzhou 341000， China； 3. Sinosteel Maanshan General Institute of Mining Research Co.， Ltd.， Maanshan 243000， China） AbstractBased on the engineering background of open pit mining transiting to the underground mining of room and pillar for the typical gently-inclined soft interlayer medium-thick phosphate rock bed surrounding Yunnan Dianchi lake， through the field survey sampling and indoor rock mechanics tests，multi-software coupled modeling was made，includ⁃ ing CAD，ANSYS，FLAC3Detc．，to realize dynamic simulation and comparative analysis on the displacement and stress re⁃ sponse characteristics of overlying strata with a slope height of 50 m and 80 m when the inclination angle is at 45． The re⁃ sults show that when the final slope is 50 m at the open pit，the impact of underground mining on the overburden is smaller than that of the 80 m final slope． In this case，the surrounding rock is relatively stable，the plastic zone area is small，and the whole rock mass is in a relatively stable state． When the final slope at the open pit is 80 m，the mining has a great distur⁃ bance to the overburden rock． After the mining，the ore pillar is in a critical failure state，and the plastic zone further devel⁃ ops and extends to the slope compared with that of 50 m， which is not conducive to the safety and stability of the mine． Based on this study， the reasonable height range of the slope was determined， which has important theoretical guiding significance for the safety and economic production of the mine． KeywordsTransition from open pit to underground mining，Room and pillar ，Overlying strata mining re⁃ sponse，Numerical simulation 总第 530 期 2020 年第 8 期 金属矿山 METAL MINE Series No． 530 August2020 32 ChaoXing 磷矿是我国重要的战略资源， 自上世纪50年代 开始， 经过数十年持续高强度的开采， 绝大多数露天 矿山已经进入深凹露天开采阶段， 大部分矿山正在 或者已经转入深凹地下开采。露天矿进入深凹开采 阶段后， 一方面采场水平不断延深， 采场作业尺寸逐 渐缩小， 工作宽度变窄， 运输效率降低； 另一方面露 天边坡越来越高， 边坡控制和管理的难度愈加困难， 对安全生产的威胁日益加剧。因此， 随着露天坑的 延伸扩大， 剥采比增幅显著， 并形成了深凹底坑和陡 峭的露天终了边坡， 使得矿山开采不仅技术难度加 大， 而且无经济效益可言。矿山露天转地下开采后， 露天终了边坡及坑底周围岩体与地下开采组成统一 的复合开采系统。露天终了边坡及坑底周围岩体在 露天开采时期受到开挖卸荷的应力扰动， 在局部地 段产生应力集中， 并在坡体内形成潜在滑坡体。在 此基础上再转入地下开采阶段， 地下矿体开挖对其 造成二次扰动， 并随着开挖空间的推进， 二次扰动效 果的动态叠加， 不断地影响并改变其应力与变形状 态。因此， 露天转地下开采是一项复杂的系统工程， 它不是一个孤立的过程， 需要把露天和地下作为一 个整体， 通过科学规划和完美衔接， 实现总体安全高 效开采 ［1-2］。 本项目针对云南某缓倾斜含软弱夹层中厚磷矿 床露天转地下开采存在的问题， 通过现场勘查取样， 室内岩石力学试验以及 FLAC3D数值模拟 ［3］， 对在边 坡影响下缓倾斜中厚磷矿体覆岩的采动响应特征进 行研究分析， 以确定露天终了边坡的合理高度范围。 1工程概况 云南某缓倾斜含软弱夹层中厚磷矿体的埋藏深 度100～240 m， 矿体厚度12～14 m， 平均倾角16， 属于 典型的近浅埋缓倾斜中厚磷矿床。边坡倾向与矿体 和岩质倾向反向， 目前露天开采高度约为40 m， 坡度 为45， 坡体内无大的构造结构和软弱面， 组成边坡 的岩层均属于中硬坚硬岩体， 属于典型的中高稳 定性反倾向岩质边坡。 2模型与开采方案建立 2. 1参数确定 按一定密度在磷矿体、 采场围岩及上覆岩体中 取块状岩样， 依据国土资源部2015年颁布的 岩石物 理力学性质试验规程 ， 将获取的岩石样品在实验室 中加工成径高尺寸为 50 mm25 mm和50 mm100 mm 的圆柱形标准试件 ［4］， 通过巴西劈裂试验 （岛津 AGI-250 伺服材料实验机， 每种岩石取4 组， 每组4 个标样） 计算得岩体抗拉强度， 其破坏后的试件如图 1所示； 通过三轴压缩试验 （TAW-2000微机控制电液 伺服岩石三轴试验机， 取岩石4组， 每组4 个标样） 测 得在不同围压下的轴向破坏强度， 破坏后试件如图2 所示。计算得出边坡、 磷矿体及地下采场围岩及上 覆岩体的主要岩石力学性质与变形特征参数 （弹性 模量、 内聚力、 内摩擦角、 抗拉强度、 抗压强度、 泊松 比） 。然后根据Hoek-Brown经验公式 ［5-6］， 对计算出 的参数进行适当折减， 得到最终所使用的岩石力学 参数见表1。 2. 2三维建模 运用 AutoCAD制图软件画好平面模型， 然后导 入到ANSYS进行网格划分， 最后将生成的网格模型 导入到FLAC3D软件进行数值计算。根据研究内容， 对模型进行非均匀的网格划分， 对需要重点研究的 岩层如磷矿层及其直接和间接顶底板进行网格加 密， 其他岩层网格单元尺寸可按照规则适当加大， 坡 高50 m最终模型尺寸为350 m60 m150 m， 模型 共划分123 058个节点、 117 125个六面体单元； 坡高 80 m为350 m60 m180 m， 模型共划分128 882个 节点、 122 350个六面体单元 ［7］。 2. 3边界条件 为了尽可能模拟工程实际情况， 约束模型前后 左右四面的水平向位移， 底部约束三向位移， 顶部自 由， 在计算的时候不考虑其他构造应力的作用， 只考 虑自重应力， 因此在模型垂直方向施加反向应力加 速度。在自重应力下的初始地应力平衡结果见图3， 2020年第8期罗浪等 露天转地下开采不同坡高影响下覆岩采动响应特征研究 33 ChaoXing 其应力值整体上随着坡面的起伏而起伏， 基本符合 实际情况。 2. 4开采方案 根据控制单一变量思想， 除坡高这一变量不同 之外， 其他条件相同。开采方式采取房柱法开采， 境 界矿柱厚度为30 m， 矿房长度10 m， 矿柱宽5 m， 为连 续型矿柱。整个矿体开挖长度为145 m， 矿房分10步 从左至右依次开挖， 矿柱分5步从右至左隔一回采。 同时对边坡和开采区域设置11个变形监测点进行位 移分析， 监测点设置如图4所示。 3FLAC3D数值模拟与结果分析 3. 1应力分析 为了研究不同边坡高度影响下的采场覆岩应力 分布特征， 需要对模型在每个开挖步下的应力变化 进行计算、 监测和记录。2种模型矿房5步开挖完后 垂直应力分布如图5所示； 矿房10步开挖完后垂直应 力分布如图6所示； 矿柱回采完毕后垂直应力分布如 图7所示。 通过以上应力云图可以看出， 随着地下开挖的 进行， 开挖过后的矿房形成采空区， 其上方顶板形成 应力卸压区， 顶底板位置应力值降低， 这是由于磷矿 层开挖后形成临空面， 围岩应力得到释放导致。同 金属矿山2020年第8期总第530期 34 ChaoXing 时， 采场内部应力向前方岩体转移， 在前方形成应力 增压区。5步开挖后采场顶板垂直最大应力值50 m 坡高为 43 kPa， 80 m 坡高为 47 kPa， 增幅达到 9.5； 10步开挖完毕后50 m坡高垂直应力增大到55 kPa， 80 m坡高增大到64 kPa， 80 m模型相对50 m模型增 幅达16.3。随着开挖的逐步推进， 50 m坡高和80 m 坡高最大垂直应力值也逐步增大， 但总体坡高80 m 较坡高50 m顶板应力值更大， 围岩稳定性更低， 安全 性下降。 3. 2塑性区分析 图 8～图 10为各阶段开挖下采场塑性区分布情 况。由图可知， 矿房开挖过后形成采空区， 围岩应力 得到释放， 顶板受到拉应力作用， 出现拉伸破坏屈服 区， 而矿柱受到剪切力作用出现剪切破坏屈服区。 随着开挖的逐步进行， 塑性区面积也在逐步扩大。 相比较于坡高50 m， 坡高80 m回采完毕后塑性区进 一步发育， 已波及坡面， 这将不利于采场稳定， 在外 部扰动影响下， 加大了岩体失稳风险。同时回采完 毕后的矿柱处于临界破坏状态， 主要为剪切破坏， 此 时在外部因素的扰动下， 可能导致矿柱发生剪切破 断， 从而引发多米诺骨牌效应导致其他矿柱的损坏， 进而使得采场整体垮塌破坏。 3. 3位移分析 模型在未进行开挖之前， 位移分布基本跟岩层 的分布一致， 并出现一定程度下沉， 基本符合客观规 律。图11～图12分别为2种坡高模型矿房开挖完毕 及矿柱回采完毕后的位移沉降云图； 图13、 图14分别 为2种模型的监测点在各开挖步下的顶板沉降曲线 图。 由图11～图14可以看出， 随着开挖的逐步进行， 岩体内的应力也在不断进行转移， 在前10步开挖中， 顶板的最大沉降位置随着开采的向前推进而向前移 动， 无论是坡高 50 m 还是坡高 80 m， 规律都比较明 显， 回采完毕后位移沉降最大值都发生在6号监测点 的位置， 且呈现出 “拱” 型。 坡高 50 m 开挖完毕后最大位移沉降值为 172 2020年第8期罗浪等 露天转地下开采不同坡高影响下覆岩采动响应特征研究 35 ChaoXing mm， 回采完毕后最大位移沉降值为260.8mm， 采场覆 岩相对稳定； 坡高80 m开挖完毕后最大位移沉降值 为275 mm （相比坡高50 m增幅达59.8） ， 回采完毕 后最大位移沉降值为397.0 mm （相比坡高50 m增幅 达52.7） ， 此时采场覆岩达到比较充分的采动阶段， 部分采场可能出现顶板冒落、 片帮等现象 ［8］。 4讨论 目前国内外学者对矿山露天转地下开采的研究 主要集中在境界矿柱厚度的选取、 开采后坑底覆盖 层与隔离顶柱的稳定性与适宜厚度、 采场结构参数 优化等， 对于地下采场在露天终了边坡坡高及坡角 影响下的覆岩采动变形破坏特征研究较少 ［9］。由 于矿山露天转地下开采是一项复杂的系统工程， 受 到矿山地质赋存条件、 采矿工艺、 采矿方法等多方 面影响， 当前露天转地下开采后矿山边坡与地下采 场围岩及其上覆岩体的应力分布特点、 变形破坏特 征及失稳机理的研究成果尚未形成系统的理论体 系， 露天终了边坡与地下采场围岩在开采扰动下的 相互作用关系尚不明确， 这也是下一步需要研究的 内容。 5结论 以云南滇池周边典型的缓倾斜含软夹层中厚磷 矿床露天转房柱法地下开采为研究背景， 通过现场 勘查取样和室内岩石力学试验， 利用CAD、 ANSYS和 金属矿山2020年第8期总第530期 36 ChaoXing FLAC3D等软件构建在倾角为45条件下露天终了边 坡为50 m和80 m的2种不同坡高的仿真分析模型， 在房柱法开采扰动下对上覆岩层位移与应力响应特 征进行分析研究。研究表明 在露天终了边坡为50 m时， 地下开采对覆岩影响相较80 m坡高小， 此时围 岩相对稳定， 塑性区面积较小， 整个岩体处于相对稳 定状态。露天终了边坡为80 m时， 开采对覆岩的扰 动较大， 回采完毕后矿柱处于临界破坏状态， 塑性区 较50 m坡高进一步发育并延伸至坡面， 不利于矿山 的安全稳定。因此， 综合经济和安全性考虑， 建议控 制露天终了边坡高度在50 m至80 m之间。 参 考 文 献 胡建军．基于FLAC3D的露天转地下磷矿山矿块结构参数优化 ［J］ ．矿业研究与开发， 2013， 33 （4） 1-3． Hu Jianjun． Optimization of structural parameters of ore block for the transition from open-pit to underground mining in a phosphate mine based on FLAC3D［J］ ． Mining Research and Development， 2013， 33 （4） 1-3． 张东明， 郑彬彬， 李耀基， 等．磷矿露天转地下开采围岩稳定性 数值模拟研究 ［J］ ．矿业研究与开发， 2014， 34 （4） 24-26． Zhang Dongming，Zheng Binbin，Li Yaoji，et al． Numerical simu⁃ lation for stability of phosphate mine under transition from open-pit to underground mining［J］ ． Mining Research and Development， 2014， 34 （4） 24-26． 王云飞， 钟福平．露天转地下开采边坡失稳数值模拟与实验研 究 ［J］ ．煤炭学报， 2013， 38 （S1） 64-69． Wang Yunfei，Zhong Fuping． Study on slope instability in transi⁃ tion from open-pit to underground mining by similar experiment and numerical simulation ［J］ ． Journal of China Coal Society， 2013， 38 （S1） 64-69． 刘艳章， 蔡原田， 张奎， 等．露天转地下开采充填体变形参数 对边坡稳定性影响研究 ［J］ ．矿冶工程， 2019， 39 （1） 1-4． Liu Yanzhang， Cai Yuantian， Zhang Kui， et al．Study on influence of deation parameters of backfill on slope stability in open-pit mining ［J］ ．Mining and Metallurgical Engineering， 2019， 39 （1） 1-4． 吴恩桥， 刘艳章， 潘杨月， 等．隔离层厚度对露天最终边帮破坏 形式及稳定性影响研究 ［J］ ．矿冶工程， 2018， 38 （6） 12-16． Wu Enqiao， Liu Yanzhang， Pan Yangyue， et al． Study on the influ⁃ ence of the thickness of isolation layer on the failure and sta⁃ bility of the open end fringe ［J］ ． Mining and Metallurgical Engineer⁃ ing， 2018， 38 （6） 12-16． 杨振增， 张纪伟， 亓中华．露天转地下开采顶板安全性分析 ［J］ ． 中国矿山工程， 2018， 47 （5） 44-46． Yang Zhenzeng，Zhang Jiwei，Qi Zhonghua． Safety analysis of roof from open to underground mining［J］ ． China Mine Engineering， 2018， 47 （5） 44-46． 张鹏， 于兆成， 宋晓丽， 等．复杂应力条件下露天矿高边坡动 态开采模拟与响应分析 ［J］ ．煤矿安全， 2018， 49 （11） 233-236． Zhang Peng，Yu Zhaocheng，Song Xiaoli，et al． Dynamic mining simulation and response analysis of high slope in open-pit mine un⁃ der complex stress condition［J］ ．Safety in Coal Mines， 2018， 49 （11） 233-236． 王宇， 胡天寿．基于数值模拟的露天-地下联合开采隔离层稳 定性分析 ［J］ ．矿业研究与开发， 2019， 39 （5） 24-29． Wang Yu， Hu Tianshou．Stability analysis on isolation layer in com⁃ bined mining of open-pit and underground based on numerical sim⁃ ulation［J］ ． Mining Research and Development， 2019， 39 （5） 24- 29． 尹光志， 李小双， 魏作安， 等．边坡和采场围岩变形破裂响应特 征的相似模拟试验研究 ［J］ ． 岩石力学与工程学报， 2011， 30 （S1） 2913-2923． Yin Guangzhi，Li Xiaoshuang，Wei Zuoan，et al． Similar simula⁃ tion study of deation and failure response features of slope and stope rocks ［J］ ． Chinese Journal of Rock Mechanics and Engineer⁃ ing， 2011， 30 （S1） 2913-2923． ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ ［8］ ［9］ 2020年第8期罗浪等 露天转地下开采不同坡高影响下覆岩采动响应特征研究 37 ChaoXing