大规模充填采矿采场稳定性研究与结构参数优化_刘建东.pdf

大规模充填采矿采场稳定性研究与结构参数优化 刘建东 1 解联库 2 曹辉 2 （1. 中国矿业大学矿业工程学院， 江苏 徐州221116； 2. 北京矿冶科技集团有限公司， 北京102628） 摘要甲玛铜多金属矿二期地下开采采用大直径深孔阶段空场嗣后充填采矿法， 为有效控制地压灾害， 实现 矿体安全高效开采， 采用理论计算、 Mathews稳定图和FLAC3D数值模拟多种方法对采场顶板及矿柱的稳定性进行综 合研究， 在此基础上确定合理的采场结构参数。结果表明， 12 m、 15 m和18 m 的采场宽度均能满足矿柱稳定的要 求； 当采场宽度为18 m时， 最大允许采场长度为77.3 m； 当采场长度为60 m时， 最大允许采场宽度为22.4 m。综合考 虑采场稳定性以及深部高应力环境等因素， 甲玛铜多金属矿采场宽度采用15 m较为合适， 推荐的采场结构参数为 15 m60 m。 关键词Mathews法数值模拟采场稳定性结构参数 中图分类号TD853文献标志码A文章编号1001-1250 （2018） -12-010-04 DOI10.19614/ki.jsks.201812002 Study on Structural Parameters Optimization and Stability of Stope for Large-scale Backfill Mining Liu Jiandong1Xie Lianku2Cao Hui2 （1. School of Mines， China University of Mining and Technology， Xuzhou 221116， China； 2. BGRIMM Technology Group， Beijing 102628， China） AbstractThe two-stage underground mining of Jiama Copper Polymetallic Mine adopts the large diameter deep hole stage mining with subsequent backfilling. In order to effectively control the ground pressure disaster and realize the safe and efficient mining of the ore body， the stability of the roof and pillar of the stope is synthetically studied by a variety of s including the theoretical calculation，Mathews stability map and FLAC3Dnumerical simulation. On this basis，the reasonable structure parameters of the stope are determined. The results show that the width of 12 m， 15 m and 18 m can meet the requirements of pillar stability. When the width of the stope is 18 m， the maximum allowable stope length is 77.3 m. When the length of the stope is 60 m， the maximum allowable width of the stope is 22.4 m. Considering the stability of the stope and the deep high stress environment，15 m is more suitable for the width of the Jiama Copper Polymetallic Mine，and the stope structure parameters are recommended as 15 m60 m. KeywordsMathews ， Numerical simulation， Stope stability， Structural parameters 收稿日期2018-11-15 作者简介刘建东 （1980-） ， 男， 高级工程师， 博士研究生。 总第 510 期 2018 年第 12 期 金属矿山 METAL MINE Series No. 510 December 2018 大直径深孔采矿技术以其高效率、 高强度、 集中 作业和改善作业环境等特点直接推动了采矿大型 化、 连续化的发展， 在世界范围内的地下矿山得到了 广泛的应用和发展 ［1］。以大直径深孔阶段空场嗣后 充填采矿法为代表的大规模采矿技术在大量开采和 集中强化开采方面发挥了重要作用 ［2］。甲玛铜多金 属矿矿体品位高、 赋存量大、 厚度适中、 矿体连续性 好， 适合采用大直径深孔阶段空场嗣后充填采矿法 回采。实际应用中， 选择合理的采场结构参数是控 制地压灾害、 实现矿体安全高效开采的重要保证 ［3］， 必须结合具体的工程地质条件和地应力环境对采场 稳定性和结构参数进行分析和优化。本研究采用理 论计算、 Mathews稳定图和数值模拟多种方法对甲玛 铜多金属矿的采场稳定性进行综合分析， 确定合理 的采场结构参数， 为大直径深孔阶段空场嗣后充填 采矿法的实施应用提供参考和指导。 1工程概况 甲玛铜多金属矿位于海拔4 350 ～ 5 407.5 m， 属 于典型的高原矿山， 矿体主要为矽卡岩型， 赋存标高 3 700 ～ 4 550 m， 矿体厚度10 ～50 m， 最厚达199 m。 受推覆构造控制， 矿体呈上陡下缓特征， 较陡部分矿 体靠近地表， 主要为铅锌矿体， 倾角60 ～ 70 ； 较缓部 10 ChaoXing 分为地下深部隐伏铜钼矿体， 倾角小于20 。矿体、 顶底板属坚硬、 中硬岩石， 浅部节理裂隙较发育， 总 体稳固程度为中等以上。 矿山二期地下开采设计采用大直径深孔高分段 空场嗣后充填采矿法， 包括北区4 400～4 500 m矿段 范围内的4 420～4 465 m、 4 460～4 520 m两区段。先 采4 460～4 520 m区段， 后采4 420～4 465 m区段。盘区 内采场回采顺序采用 “隔三采一” 的方式， 一、 二步骤回 采矿房采场， 采用灰砂比为1 ∶ 8的尾砂胶结充填， 三、 四步骤回采矿柱采场， 采用尾砂非胶结或废石充填。 2工程岩体质量RMR分级 工程岩体质量的好坏直接关系到岩体工程稳定 性， 工程岩体质量评价是进行工程设计与现场施工 的基础依据。根据对甲玛铜多金属矿大规模开采区 域的工程地质调查， 结合该区域的工程地质特征， 采 用国内外常用的RMR （CSIR） 岩体分类系统对角岩和 矽卡岩进行分类， 并使用经验公式 ［4］ Qe RMR-44/9分别 计算对应的Q值， 结果如表1所示。 3基于矿柱稳定的合理采场宽度计算 矿柱稳定性是大规模开采安全的关键问题之 一。采用国内外常用的Bieniawski矿柱强度公式 ［5］进 行大规模开采的采场及盘区矿柱合理参数的论证优 化。由于矿柱沿走向长度较长， 可简化为条带式矿 柱计算矿柱所受到的垂直应力， 按下式计算 σPγH1wo/wp，（1） 式中，γ为覆岩容重， kN/m3； H为顶板埋深， m；wp为矿 柱宽度， m；wo为矿房宽度， m。 Bieniawski与Vanheerden对南非Witbank煤矿宽 高比为0.5～34的多个煤柱试件进行了大量的原位测 试， 得出的矿柱强度计算公式为 SPSl[0.640.36 wP h ]α，（2） 式中，Sl为矿柱抗压强度， MPa； h为矿柱高度， m；α取 值根据矿柱高宽比h/wP确定， 当h/wP≥5时， 取α1.4， 当h/wP＜5时， 取α1.0。 甲玛矿大规模开采区域矿柱高宽比均小于5， 故 取α1.0。因此， 矿柱的安全系数计算公式如下 K SP σP Sl[0.640.36 wP h ] γH1wo/wp .（3） 根据岩石力学试验结果，Sl取值为75 MPa； 将采 高50 m， 采场宽分别为12 m、 15 m、 18 m的参数代入 式 （3） 中， 可以得出不同采场宽度、 不同埋深条件下 矿柱安全系数值， 如图1所示。根据甲玛铜多金属矿 矿体条件， 采场顶板最大埋深约600 m， 在此情况下 矿柱安全系数仍能满足大于临界安全系数1.5的要 求， 说明12 m、 15 m、 18 m宽度均能够满足矿柱稳定 的要求。上述计算是在矿柱两侧临空的极限状态下 得出的， 若在采场充填状态下， 矿柱将更加稳定。 4基于Mathews稳定图的采场稳定性分析 4. 1Mathews稳定图法 Mathews稳定图法是由Mathews等人于1981年最 早提出 ［6］， 它是依据大量工程实例提出的经验公式 法。1988年Potvin在分析242个案例和重新定义某 些调整系数的基础上， 对Mathews的稳定性图表进行 修正 ［7］。1992年Potvin， Nickson等人通过收集更多深 部采矿现场数据资料， 对该方法的合理性进行验证 并提出修正 ［7-9］； 2000年Trueman等人根据大量新增 的实例资料， 采用对数回归的方法重新定义稳定区 和严重破坏区 ［10］； 2004年Mawdeskey等人给出了该方 法稳定区、 破坏与严重破坏区的等概率图 ［11］。 Mathews稳定图法实质是利用NGI岩体分级指 标Q计算岩体稳定性指数N， 依据矿山开拓和采准工 程确定采场结构参数并计算采场暴露面形状系数S， 将N和S值投影到Potvin修改后的稳定性图表上， 即 可初步判断采场的总体稳定性； 或者是根据岩体稳 定性指数N在稳定性图表上求出总体稳定的采场形 状系数S， 在初步选定采场某一结构参数后即可确定 其他结构参数。 刘建东等 大规模充填采矿采场稳定性研究与结构参数优化2018年第12期 11 ChaoXing 4. 2Mathews法计算参数选取 （1）Q′值 以节理裂隙水折减系数Jw与应力折减 系数SRF之比等于1时的Q值即为Q′值。根据前述 计算矽卡岩Q4.74， 依据地应力计算结果， SRF 取 1.556，Jw取0.66， 计算得Q′为11.17。 （2） 岩石应力系数A 通过岩石单轴抗压强度σc和 次生压应力σ1之比确定。根据模拟分析， 采场顶板中 央位置的压应力平均为12.85 MPa， 顶板角岩单轴抗压 强度为113.53 MPa， 计算得σc/σ18.84， 故A值取0.847。 （3）节理产状调整系数 B 角岩平均倾角取 74.5， 故B值取0.897。 （4） 重力影响系数C 采场顶板均为水平状态， 故 C值取1.0。 4. 3采场结构参数Mathews法分析 根据以上分析， 按公式NQ ′ABC计算得， N8.49。根据图1计算得出采场顶板保持稳定的最 大水力半径S为7.65 m， 在此允许水力半径条件下， 采场顶板暴露面积和长度随采场宽度的变化趋势见 图 2。设计推荐的 45 m18 m、 52.5 m18 m 和 60 m18 m 三种结构尺寸， 对应的水力半径分别为 6.43、 6.7、 6.92， 均小于允许水力半径， 故可以满足采 场稳定的要求。从图2中可以看出， 当采场宽度为 18 m时， 最大允许采场长度为77.3 m； 当采场长度为 60 m时， 最大允许采场宽度为22.4 m。 5采场结构参数的数值模拟 5. 1计算模型 综合考虑矿体及围岩分布特点和岩层移动角等 因素， 建立的三维模型尺寸为 200 m420 m200 m， 模型共划分110 880单元和118 728节点。根据甲 玛铜多金属矿的实测地应力结果， 赋予模型实测地 应力。采用理想弹塑性Mohr-Coulumb模型。计算的 范围位于4 420～4 465 m中段， 计算模型的岩体全部 采用矽卡岩， 矿房采用胶结充填， 矿柱采用非胶结充 填， 具体物理力学参数如表2所示。本次计算假定充 填体是接顶的， 在计算中仅考虑地应力的作用， 忽略 地震波、 爆炸冲击波及地下水等因素对岩体稳定性 的影响。计算过程采用隔 “三采一” 的开挖顺序， 每 个采场采完后立即充填。 5. 2应力变化分析 在矿房回采期间， 矿柱中发生了明显的应力集 中现象， 且首采、 次采矿房开采时， 矿柱中所产生的 应力集中程度不同， 主要原因是首采矿房开采后， 应 力转移到邻近的围岩体中， 围岩体完整性相对较好， 能够赋存住高应力； 待次采矿房开采后， 矿房之间的 岩体受力进一步恶化， 高应力无法全部得到赋存， 出 现部分应力损失和能量耗散， 这部分应力和能量转 移到矿房的胶结充填体中， 使充填体发展为三向受 力状态。对于高应力区域的盘区采场， 采场宽度较 小时矿柱容易发生屈服破坏， 故采用15 m、 18 m宽度 的矿柱尺寸相对较好。 5. 3顶板变形分析 如图3所示， 在矿块内各个采场回采过程中， 位 于矿块中央位置采场顶板变形量呈线性增长趋势。 矿房、 矿柱的回采对于邻近区域的覆岩都具有一定 程度的影响。在矿块充填后一段时间， 矿块中央位 置采场顶板变形量趋于稳定。采场宽度采用12 m、 15 m时， 顶板变形量较小； 采场宽度采用18 m时， 顶 板变形相对较大， 且对比明显。因此， 从覆岩移动的 剧烈程度来看， 采场宽度采用12 m或15 m有利于控 制覆岩变形。 金属矿山2018年第12期总第510期 12 ChaoXing ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ ［8］ ［9］ ［10］ ［11］ 5. 4塑性区分布 在矿房开采胶结充填后， 3种采场宽度下的顶板 覆岩的塑性区发展深度较小， 大多数塑性单元恢复 到弹性工作状态； 采场结构参数为12 m时， 矿柱的剪 切破坏塑性区分布明显， 以低围压下发生的压剪、 矿 柱表面的拉剪破坏为主， 大量塑性单元未恢复到弹 性工作状态， 而采场宽度15 m、 18 m时矿柱的受力状 态相对较好， 尤其是18 m矿柱部分核区仍以弹性工 作状态为主， 这与前述的应力分析结果一致。当矿 块被完全采空非胶结充填后， 覆岩变形相对剧烈， 但 最终恢复到弹性工作状态； 胶结充填体发挥一定的 支承作用， 充填体出现大量的塑性单元， 但最终会趋 于稳定。 5. 5采场结构参数数值模拟结果分析 通过对4 420～4 465 m中段矿块开采过程的采场 应力、 变形和塑性区分析可以看出， 对于高应力条件 下的采场， 采用12 m的采场宽度， 在矿房采空条件下 发生矿柱屈服失稳的可能性相对较大； 而采用18 m 的采场宽度， 在开采过程中及采后覆岩移动相对剧 烈， 尤其是在采矿过程中采场顶板容易发生破断。 综合分析， 在深部高应力的条件下， 甲玛铜多金属矿 采场宽度采用15 m较为合适。结合Mathews法分析 结果， 推荐的采场结构参数为15 m60 m。 6结论 从多方法、 多角度对甲玛铜多金属矿大规模开 采采场的稳定性进行分析， 对采场结构参数进行优 化， 研究结果对大直径深孔阶段空场嗣后充填采矿 法在甲玛矿的实施应用具有一定指导意义和参考价 值。主要结论如下 （1） 角岩RMR值为56， 矽卡岩RMR值为58， 二者 均属Ⅲ类岩体， 质量较好。 （2） 通过理论计算， 在确保矿柱稳定的前提下， 18 m宽的采场结构参数是可行的。通过Mathews法 分析可知， 设计推荐45 m18 m、 52.5 m18 m、 60 m18 m三种结构参数均可行。在确保采场稳定的 前提下， 当采场宽度为18 m时， 最大允许采场长度为 77.3 m； 当采场采场长度为60 m时， 最大允许采场宽 度为22.4 m。 （3） FLAC3D数值模拟分析结果表明， 采用12 m的 采场宽度， 在矿房采空条件下发生矿柱屈服失稳的 可能性相对较大； 而采用18 m的采场宽度， 在开采过 程中及采后覆岩移动相对剧烈。综合分析， 在深部 高应力的条件下， 甲玛铜多金属矿采场宽度采用15 m较为合适， 推荐的采场结构参数为15 m60 m。 （4） 采用理论计算、 Mathews稳定图和数值模拟 多方法、 多角度对甲玛铜多金属矿采场稳定性和结 构参数进行分析与优化， 不同方法的分析结果互相 验证、 互为补充， 克服了单一方法分析的不足。 参 考 文 献 孙忠铭， 刘庆林， 余斌， 等. 地下金属矿山大直径深孔采矿技 术 ［M］ . 北京冶金工业出版社， 20142-3． Sun Zhongming， Liu Qinglin， Yu Bin，et al. Large Diameter Long Hole Mining Technology in Underground Metal Mines ［M］ .Beijing Metallurgical Industry Press， 20142-3. 陈何， 孙忠铭.束状孔大量高效采矿技术的开发与应用 ［J］ .金 属矿山， 2010 （11） 1-4. Chen He， Sun Zhongming.Development and application of the high efficiency bunch hole mass mining technique ［J］ . Metal Mine， 2010 （11） 1-4. 王新民， 李洁慧， 张钦礼， 等. 基于FAHP 的采场结构参数优化 研究 ［J］ . 中国矿业大学学报， 2010， 39 （2） 163-168． Wang Xinmin， Li Jiehui， Zhang Qinli， et al. Optimizing mining stope structural parameters using a FAHP ［J］ . Journal of China Uni- versity of Mining Technology， 2010， 39 （2） 163-168． Bieniawski Z T. Engineering Rock Mass Classifications A Com- plete Manual for Engineers and Geologists in Mining， Civil and Pe- troleum Engineering ［M］ . New York Wiley， 1989． 贾喜荣， 王丽. 回采巷道煤柱临界宽度理论计算方法 ［J］ .太原 理工大学学报， 2011， 42 （1） 102-103. Jia Xirong， Wang Li. The theoretical calculation of critical width for the coal pillar of the roadway ［J］ . Journal of Taiyuan Uni- versity of Technology， 2011， 42 （1） 102-103. Mathews K E， Hoek E， Wyllie D C， et al. Prediction of Stable Exca- vations for Mining at Depths below 1 000 m in Hard Rock［R］ . CANMET Report DSS Serial No OSQ80-00081， DSS File No 17SQ 23440-0， 1981， 9020． Potvin Y. Empirical Open Stope Design in Canada ［D］ . Vancouver University of British Columbia， 1988． Stewarts B V， Forsyth W W. The Mathews for open stope design ［J］ . CIM Bulletin， 1995， 88 （1） 45-53． Brown E T. Block Caving Geomechanics ［M］ . Queensland， Austra- lia Julius Kruttschnitt Mineral Research Centre， 2003 3-12． Trueman R， Mikula P， Mawdesley C， et al．Experience in Australia with the Mathews for open stope design ［J］ . CIM Bulletin， 2000， 93 （2） 162-167． Mawdesley C. Using logistic regression to investigate and improve an empirical design ［J］ . International Journal of Rock Me- chanics and Mining Sciences， 2004， 41 （S1） 756-761． （责任编辑石海林） 刘建东等 大规模充填采矿采场稳定性研究与结构参数优化2018年第12期 13 ChaoXing