分段空场嗣后充填采矿法集约化开采顺序优化研究_任卫东(1).pdf

收稿日期2020-02-25 基金项目国家重点研发计划项目 （编号 2016YFC0801606） 。 作者简介任卫东 （1980） ， 男， 高级工程师。 总第 528 期 2020 年第 6 期 金属矿山 METAL MINE 分段空场嗣后充填采矿法集约化开采顺序 优化研究 任卫东范玉乾 1 （长沙有色冶金设计研究院有限公司， 湖南 长沙 410019） 摘要某矿山是一个水文地质条件复杂， 生产多年的老矿山， 矿山采用空场法开采， 多中段同时生产， 采场 多， 生产效率低， 时有采场顶板冒落等问题发生， 严重影响企业的安全生产。为了解决这些问题， 企业采用了分段 空场嗣后充填采矿法， 虽然大大改善了采场回采的安全， 但矿山采用多中段生产， 造成中段间相互影响， 生产效率 低， 且存在潜在安全隐患， 本研究在分段空场嗣后充填采矿法的前提下， 采用数值模拟软件， 依据摩尔库伦破坏 准则， 通过对拟定的不同开采顺序采场回采进行数值模拟， 分析不同开采顺序采场的应力及位移的变化情况， 从而 能选择适合矿山的开采顺序， 为矿山实现集约化、 规模化开采提供了依据， 对矿山的安全、 高效、 经济的回采有至关 重要的意义， 也为类似矿山的开采提供了借鉴依据。 关键词分段空场法嗣后充填采矿法集约化开采开采顺序优化 中图分类号TD853文献标志码A文章编号1001-1250 （2020） -06-167-05 DOI10.19614/ki.jsks.202006026 Study on Optimization of Mining Sequence in Intensive Mining with Sublevel Open Stope and Subsequent Filling Ren WeidongFan Yuqian2 （Changsha Nonferrous Metallurgy Design and Research Institute Co.， Ltd.， Changsha 410019， China） AbstractA mine is an old mine with complex hydrogeological conditions and many years of production． It is mined by open stope ． Many middle sections are produced at the same time． There are many stopes and low production efficien⁃ cy． Sometimes there are problems such as stope caving， which seriously affect the safety production of the enterprise． In order to solve these problems，the enterprise adopts the sublevel open stope followed by filling mining ． Although the safety of stope mining is greatly improved，the mine adopts the multi sublevel production，which causes the interaction between the middle sections， low production efficiency and potential safety risks． In this paper， on the premise of sublevel open stope fol⁃ lowed by filling mining ，and according to the Moore Coulomb failure，the numerical simulation software is used to make numerical simulation on different mining sequences of the stopes． The change of stress and displacement of stopes in different mining sequences is analyzed，so that the mining sequence suitable for the mine can be selected，which provides a basis for the mine to realize intensive and large-scale mining． It is of vital significance for the safe，efficient and economic mining of the mine， and also for the mining of similar mines． The reference basis is provided for the similar mines． KeywordsSublevel open-stoping ， Subsequent filling， Mining ， Intensive mining， Optimization of min⁃ ing sequence Series No． 528 June 2020 采矿方法是矿山安全生产的核心， 采矿方法的 选择直接决定了矿山生产运营过程中的效益及安 全。现阶段随着国家政策的支持及矿山自我认识的 提高， 我国地下矿山采用充填法的比例逐渐提升。 地下采矿对矿体的开挖过程其实就是地应力的释放 过程 ［1-2］， 开挖岩体应力的释放使岩体应力进行了重 新调整和分配， 而充填方法的使用又使应力重新发 生了变化， 不同的采矿、 充填过程即形成了不同的开 采顺序， 开采范围内的矿岩介质也就有了不同的变 化过程， 对未开采区域以及已充填的区域， 就有不同 167 ChaoXing 金属矿山2020年第6期总第528期 的加载和卸载路径 ［3-4］， 即存在不同的应力场和位移 场的变化。本项目通过采用数值模拟软件， 依据摩 尔库伦破坏准则， 通过采用分段空场嗣后充填采 矿法的前提下对拟定的不同开采顺序进行数值模 拟， 分析不同开采顺序采场的应力及位移的变化情 况， 从而能选择适合矿山的开采顺序， 为矿山实现集 约化、 规模化开采提供依据， 对矿山的安全、 高效、 经 济的回采有至关重要的意义， 也为类似矿山的开采 提供了借鉴。 1工程概况 某矿床成因类型属接触交代高中温热液型矿 床， 矿体以铜铁矿石 （矽卡岩） 为主， 矿体内夹石主要 为大理岩， 次为矽卡岩。顶、 底盘围岩随矿体赋存部 位不同而不同。矿体与大理岩的接触带部位的岩体 强度相对较低， 但接触影响宽度较小， 而且结合比较 紧密。矿体与花岗闪长斑岩的接触带条件较恶劣， 岩体破碎， 容易水解， 稳定性最差。 矿山采用竖井开拓， 有轨运输。矿山现阶段主 要采用分段空场嗣后充填的采矿方法， 其中 采场垂 直矿体走向布置， 长度为矿体水平厚度， 中段厚度50 m， 分段高度12.5 m。留顶柱和底柱， 顶柱高5 m， 底 柱高8 m。在底柱内布置出矿进路， 采用铲运机平底 出矿底部结构。 2数值模拟分析 2. 1数值模型的建立 本次数值模拟采用数值软件进行建模和网格化 划分及后处理计算 ［5］。为方便建模和数值模拟计算， 建模前做以下的假设 ［6］ ①矿体和围岩为各向同性的 连续介质， 不考虑裂隙、 节理和断层的影响； ②忽略井 巷工程对采场稳定性的影响； ③只考虑重力对模型的 影响， 不考虑爆破振动、 地震波及地下水的影响。 矿体主要为大理岩、 灰岩等岩石组成， 因此适用 摩尔库伦破坏准则。 在分析时， 首先进行矿床开采前的初始应力 场 ［7-10］分析， 最终开挖引起的位移和应力为各阶段步 骤 （各分步） 引起的位移ui（i1， ， n） 与应力σi（i1， ， n） 之和， 即 { ｝u { ｝u1{ ｝u2 { ｝un，（1） { ｝σ { ｝σ0{ ｝σ1 { ｝σn，（2） 式中，{ ｝σ0为初始地应力。显然位移结果中不包括原 岩应力状态下的位移值。 根据现场工程地质调查、 室内岩石力学参数试 验结果， 数值计算采用的材料参数见表1。 为研究不同中段回采顺序对采场稳定性的影 响， 对-270 m、 -320 m、 -370 m、 -420 m4个中段的回 采顺序进行数值模拟， 共设计了2种回采方案 方案1 （M-1） 从上向下回采和方案 2 （M-2） 从下向上回 采。几何模型见图1所示。 方案1、 方案2模拟步骤如表2所示。 2. 2模拟结果分析 在采场回采顺序一致的条件下进行中段开采顺 序的比较， 计算结果从最大主应力、 垂直位移2个方 面进行分析， 以此对采空区的稳定性提供定量依据。 168 ChaoXing 2020年第6期任卫东等 分段空场嗣后充填采矿法集约化开采顺序优化研究 2. 2. 1应力分布状态的分析与比较 应力分布中， 上盘岩体最大主应力和顶板主应 力最为重要， 因其对采场稳定性的影响最大。 由表 3、 表 4 和图 2～图 5 可知， 对于采场顶板而 言， 在-270 m中段回采时， 两方案采场顶板的主应力 和集中应力差别不大； 在-320 m、 -370 m、 -420 m中 段回采时， 采用从上至下的中段回采顺序 （方案1） ， 顶板最大主应力、 中间主应力和最小主应力均较小， 顶板最大主应力分别为 4.60 MPa、 5.17 MPa、 5.96 MPa， 均未超过矿体最大抗拉强度6 MPa， 且拉应力区 面积也小； 而采用从下至上的中段回采顺序 （方案2） 时， 在-320 m、 -370 m、 -420 m中段回采时， 采场顶柱 最大主应力分别为5.49 MPa、 6.19 MPa （超过矿体最 大抗拉强度6 MPa） 、 8.02 MPa （超过矿体最大抗拉强 度 6 MPa） ， 且方案 1 较方案 2 分别小 19.4、 19.7、 34.6。需要注意的是 在采场两侧的端角处出现最 大应力集中， 而且最大主应力都是拉应力， 过大的拉 应力是造成顶柱垮塌的重要原因， 方案2 应力集中现 象比方案1明显， 且方案2在回采-370 m、 -420 m中 段时， 顶板最大拉应力已超过矿体的最大抗拉强度， 相比之下， 方案1较优。对于采场上盘而言， 在-270 m、 -320 m、 -370 m中段回采时， 采用从上至下的中段 回采顺序 （方案1） ， 顶板最大主应力、 中间主应力和 最小主应力相对较小， 尤其是顶板最大主应力， 其大 小分别为0.184 MPa、 0.188 MPa、 0.192 MPa； 而采用从 下至上的中段回采顺序 （方案2） 时， 采场顶柱最大主 应力分别为0.628 MPa、 0.615 MPa、 0.552 MPa， 且方案 1 较方案 2 分别小 241、 227、 188。而在-420 m 中段回采时， 由于-420 m中段是方案2的首采中段， 上盘最大主应力相对方案1小， 而中间主应力和最小 主应力仍然比方案1大， 总的来说， 方案1较优。 注表中压应力为 “-” ， 拉应力为 “” ； 垂直位移为 “-” ， 表示沿Y轴 负方向， 反之为 “” 。 2. 2. 2位移分布状态的比较 由表 3、 表 4和图 6～图 9可知， 对于采场顶板而 言， 在-270 m、 -320 m、 -370 m、 -420 m4个中段回采 169 ChaoXing 金属矿山2020年第6期总第528期 时， 采用从上至下的中段回采顺序 （方案1） ， 顶板最 大 垂 直 位 移 均 较 小 ， 分 别 为 -44.81 mm、 -36.26 mm、 -28.76 mm、 -21.75 mm， 而采用从下至上的中 段回采顺序（方案 2）时顶板最大垂直位移分别 为-48.81 mm、 -41.18 mm、 -34.47 mm、 -26.98 mm， 方 案 1较方案 2分别小 8.9、 13.6、 19.9、 24； 对于 采场上盘而言， 同样， 方案1的上盘最大垂直位移也 均比方案2小。很明显， 由顶板和上盘的最大垂直位 移的比较， 方案1优于方案2。 3结论 本项目结合某矿山所采用的分段空场嗣后充填 法所采用的采场结构参数， 拟定了2种回采顺序， 通 过数值模拟方法， 对2个方案模型的应力、 位移等分 析可知， 采取从上向下的中段回采顺序， 各中段均处 于较低的应力状态， 对安全高效开采非常有利； 相 反， 若采取从下向上的中段回采顺序， 在开采初期， 采场受力状态较差， 由于处于高应力状态， 顶板较易 产生冒落， 甚至于发生岩崩， 从而危及采场回采安 全。故宜采取从上向下的中段回采顺序。由于不同 矿山矿岩类型和采用的采矿方法不一样， 因此本项 目得出的结果也并不一定适合其他矿山， 但通过建 立矿体几何模型， 通过数值模拟手段比较直观地模 拟回采过程及采场应力及位移的变化， 可为类似矿 山提供借鉴。 170 ChaoXing ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ ［8］ ［9］ ［10］ 参 考 文 献 陈敏， 史秀志．某石灰岩矿采场结构参数数值模拟分析 ［J］ ．采 矿技术， 2015， 15 （1） 16-21． Chen Min，Shi Xiuzhi． Numerical simulation analysis of stope structural parameters of a limestone mine［J］ ． Mining Technology， 2015， 15（1） 16-21． 杜永亮， 包东程．红岭铅锌矿采场结构参数数值模拟与优化 ［J］ ． 采矿技术， 2019， 19 （3） 7-9． Du Yongliang，Bao Dongcheng． Numerical simulation and optimi⁃ zation of stope structural parameters in Hongling lead zinc mine ［J］ ． Mining Technology， 2019， 19（3） 7-9． 谭军， 史秀志， 陈佳耀．基于数值分析的中深孔爆破采场结构 优化研究 ［J］ ．采矿技术， 2016， 16 （6） 3-6． Tan Jun，Shi Xiuzhi，Chen Jiayao． Study on optimization of stope structure based on numerical analysis［J］ ． Mining Technology， 2016， 16（6） 3-6． 苏先锋， 陈顺满．某铜铁矿二期工程采场结构参数优化研究 ［J］ ． 采矿技术， 2016， 16 （2） 18-21． Su Xianfeng，Chen Shunman． Study on optimization of stope struc⁃ tural parameters of a copper iron ore phase II project［J］ ． Mining Technology， 2016， 16（2） 18-21． 孟庆彬， 王从凯， 韩立军， 等．无底柱分段崩落法采场围岩稳定 性分析与采矿进路修复技术 ［J］ ．采矿与安全工程学报， 2016， 36 （2） 364-372． Meng Qingbin，Wang Congkai，Han Lijun，et al． Stability analysis of surrounding rock in stope with sublevel caving without sill pillar and repair technology of mining access［J］ ． Journal of Mining and Safety Engineering， 2016， 36（2） 364-372． 叶义成， 施耀斌， 王其虎， 等．缓倾斜多层矿床充填法开采围岩 变形及回采顺序试验研究 ［J］ 采矿与安全工程学报， 2015， 32 （3） 407-413． Ye Yicheng， Shi Yaobin， Wang Qihu， et al． Experimental study of deation of wall rock and stoping sequence in mining gently in⁃ clined and multilayer deposits by backfill mining［J］Journal of Mining and Safety Engineering， 2015， 32（3） 407-413． 陶干强， 任青云， 罗辉， 等．无底柱分段崩落法采场稳定性分 析 ［J］ ．岩土力学， 2011， 32 （12） 3768-3772． Tao Ganqiang ， Ren Qingyun， Luo Hui， et al． Stability analysis of stope in pillarless sublevel caving［J］ ． Rock and Soil Mechanics， 2011， 32 （12） 3768-3772． 吴少华．无底柱分段崩落法采矿进路受力过程分析 ［J］ ．化工矿 山技术， 1993 （4） 16-18． Wu Shaohua． Analysis of the stress process of sublevel caving with⁃ out sill pillar［J］ ． Chemical Mining Technology， 1993（4） 16-18． 杨继胜， 余伟健．复杂地质条件下矿体开采顺序的数值分析 ［J］ ． 金属矿山， 2011 （7） 13-18． Yang Jisheng，Yu Weijian． Numerical analysis of orebody mining sequence under complex geological conditions［J］ ． Metal Mine， 2011 （7） 13-18． 蔡晓盛， 李胜， 卢宏建， 等．石人沟铁矿地下三期矿体开采顺 序优化 ［J］ ． 矿业研究与开发， 2019， 39 （7） 62-65． Cai Xiaosheng，Li Sheng，Lu Hongjian， et al． Mining sequence op⁃ timization of the third underground ore body in Shirengou Iron Mine［J］ ． Mining Research and Development，2019， 39（7） 62- 65． 2020年第6期任卫东等 分段空场嗣后充填采矿法集约化开采顺序优化研究 171 ChaoXing