厚大矿体采场结构参数优化研究.pdf

矿 业 工 程黄金 年第 期 第 卷厚大矿体采场结构参数优化研究收稿日期：；修回日期：基金项目：福建省自然科学基金（）；紫金矿业集团股份有限公司 年第二批科技计划项目（，）作者简介：严鹏（），男，正高级工程师，硕士，从事金属矿山安全高效开采及技术标准化方向的研究工作；：严鹏，陈绍民，黄敏，黄志国，张宝（紫金矿业集团股份有限公司；贵州新恒基矿业有限公司；富蕴金山矿冶有限公司；紫金（长沙）工程技术有限公司；东华理工大学地球科学学院）摘要：合理的结构参数对维护采场稳定、提高单个采场产能、降低采切工程量具有重要意义。以蒙库铁矿为研究对象，利用岩体遥测与结构分析系统 对采场顶板、侧帮及上盘围岩节理面进行调查和数字化识别，获取采场节理组的空间方位统计信息及优势节理产状；基于 稳定图解法对沿走向布置与垂直走向布置采场的顶板、侧帮及上盘开展稳定性分析，得到采场顶板、侧帮及上盘的实际水力半径和允许水力半径的相关关系；采用数值分析方法，通过分析矿柱的位移、最大剪切应变率及采场的主应力等因素，确定了 种典型采场结构矿柱的最小安全厚度。研究结果表明：采场沿矿体走向布置时，当矿体厚度小于 时，长度控制在 之内的采场稳定性较好；采场垂直矿体走向布置时，矿体厚度不超过 ，宽度小于 的采场稳定性较好；对于沿走向布置的 典型采场，矿柱最小安全厚度为 ；对于垂直走向布置的 典型采场，矿柱最小安全厚度为 。研究结果可为蒙库铁矿采场结构参数优化提供依据，并能为国内外同类矿山采场结构设计提供参考。关键词：厚大矿体；稳定图解法；采场结构参数；数值分析；矿柱厚度中图分类号：文章编号：（）文献标志码：引言中国经济从高速增长阶段转变为高质量发展阶段，经济结构的改变导致其对矿产资源的依赖程度增加，需要消耗大量的矿产资源来支撑 。然而，中国矿产资源整体开发利用率较低，存在一定程度的资源浪费现象，制约矿产资源开发向规模化、集约化方向发展 。蒙库铁矿东段矿体总体走向 ，倾向 ，倾角 ，整体呈不连续分布，勘探线 勘探线矿体厚 ，是典型的急倾斜厚大矿体。因此，结合矿山工程地质条件、工艺技术水平及现场工程布置等情况，从保障矿山生产能力、降低采矿成本、提高采矿效率及确保采场的稳定性和安全性等方面出发，对厚大矿体合理采场结构参数的确定及优化开展研究工作具有重要意义。目前，国内外相关领域学者利用工程类比法、理论分析法和数值模拟法等方法，对采场结构参数的确定及优化进行了丰富的研究 。其中，稳定图解法和数值模拟法因其实用性及准确性得到了广泛应用 。赵兴东等 利用 种常见的岩体质量分级方法进行岩体质量分级综合评价，采用 稳定图解法对采场结构参数进行初步计算，在此基础上借助数值模拟软件对采场稳定性进行分析，最终确定了最优的采场结构参数。刘建东等 针对矿山开采过程中出现的地压灾害现象，通过理论计算、稳定图解法及数值模拟分析等方法对采场稳定性进行综合评价，得到合理的采场结构参数，确保高应力条件下矿体高效开采的稳定性。郭进平等 为了解决矿体开采过程中采场顶板的稳定性问题，应用 稳定图解法对采场结构参数进行初步选择，借助 软件对采场结构参数进行针对性的分析，为矿山提供合理的采场结构参数。本研究通过现场工程地质调查获取采场不同区域内岩体分布特征，利用 系统分级评价方法对矿（岩）体进行综合评价，在此基础上，结合 稳定图解法对不同赋存条件下的矿体采场结构参数进行初步确定，利用 数值模拟软件对采场稳定性进行分析，最终确定合理的采场结构参数，并为类似矿山采场结构参数设计提供参考依据。采场矿岩工程地质调查采用三维数字摄影测量和岩体结构分析系统 年第 期 第 卷矿 业 工 程（），对深度为 左右的矿体顶板、侧帮及下盘围岩典型区域进行节理裂隙扫描测量，统计节理裂隙分布规律并进行优势节理划分，如图 所示。根据调查结果，采场顶板矿岩岩组中控制性的节理产状主要发育有 组：，；采场顶板为近水平状（倾角取 ），产状为 。采场侧帮的矿岩岩组中控制性的节理产状主要发育有 组：，；采场侧帮为竖直状，产状为 。采场上盘围岩均为变粒岩，变粒岩岩组中主要发育 组控制性节理：，；采场上盘产状与矿体产状类似，取 （取矿体平均倾角 ）。图 典型区域节理赤平投影图根据 系统分级评价方法 ，对矿岩和变粒岩进行分级评价，值计算见式（）：（）式中：为岩石质量指标；为岩体的节理数；为岩体的节理粗糙度系数；为岩体的节理蚀变系数；为节理渗水折减系数；为应力折减系数。根据 系统参数评分表，并结合矿岩节理裂隙的实际调查结果，矿岩和变粒岩的评价结果如表 所示。表 矿岩和变粒岩的评价结果岩性 等级矿岩 级变粒岩 级 计算采场稳定性系数本文基于 稳定图解法设计采场结构参数，稳定图解法是预测采场稳定性的实用方法，在工程中得到了广泛应用。稳定图解法设计过程以稳定系数和水力半径的计算为基础，根据这 个因子将采场划分为稳定区、无支护过渡区、支护稳定区、支护过渡区和崩落区。稳定性系数反映了在一定应力条件下岩体的自稳能力，其计算见式（）：（）式中：为稳定性系数；为修正的 值；为岩石应力系数；为节理方位系数；为重力调整系数。水力半径考虑了单独采场暴露表面的尺寸和形状，水力半径可用表面积与暴露面周长的比值来表示，计算见式（）：（）（）式中：为水力半径（）；为采场的长度（）；为采场的宽度或高度（）。稳定性系数与水力半径之间的关系如图 所示，通过该图来确定容许的水力半径。图 稳定系数与水力半径相关关系）值的计算。稳定图解法使用修正后的 隧道质量指标值 ，值的计算见式（）：（）根据表 ，可计算出矿岩和变粒岩的 值分别为 和 。）系数 、的计算。岩石应力系数 与完整岩体的单轴抗压强度和平行开挖面的最大诱导应力的比值呈线性关系。值为 ，其值可根据图 ）确定。结合室内岩石力学试验及现场应力测试得到矿岩及变粒岩的值分别为 ，结合图 ），可以确定矿岩及变粒岩的岩石应力系数 的取值分别为 ，。节理方位系数 受不连续面的影响，其值可根 矿 业 工 程黄金图 系数 、的计算图据控制性节理与采场表面的相对方位确定（如图 ）所示）。根据前述的地质调查结果，可以得到采场节理方位系数 的计算结果如表 所示。表 节理方位系数 计算结果节理方位系数分析区域优势节理产状与节理方位对照产状夹角值顶板：产状夹角值侧帮：产状夹角值上盘：结合最不利原则，进行采场稳定性分析时，顶板取 （矿岩）；侧帮取 （矿岩）；上盘取 （变粒岩）。重力调整系数 考虑了重力对采场暴露表面崩落、滑落等稳定性的影响（如图 ）所示），重力调整系数和采场表面倾角的关系式为：（）式中：为采场工作面倾角（）。根据采场顶板、侧帮及上盘的产状，通过式（）可计算出重力调整系数 的值分别为 ，。根据式（）及上述系数的取值，得到相关的参数及计算结果见表 。由表 可知，采场顶板、侧帮及上盘的稳定性系数 分别为 ，。表 容许水力半径计算区域岩性 容许水力半径 稳定区无支护过渡区支护稳定区支护过渡区崩落区顶板矿体 侧帮矿体 上盘变粒岩 采场结构参数设计根据蒙库铁矿东段 以上采场实际布置情况，并结合矿体的赋存条件，采用 稳定图解法分析采场稳定性。现场工程地质调查结果表明，蒙库铁矿东段 以上矿体直接围岩以变粒岩为主，夹少量石榴子石，矿体主要为磁铁矿或块状石榴子石磁铁矿。矿体 中段采用阶段空场采矿法开采，阶段高度为 ，当矿体厚度小于 时，采场沿矿体走向布置，当矿体厚度大于 时，采场垂直矿体走向布置。沿矿体走向布置当采场沿矿体走向布置时，根据水力半径的计算方法，得出不同采场结构参数下采场顶板的水力半径，以及采场侧帮及上盘水力半径，计算结果分别如图 、表 及表 所示。图 沿矿体走向布置采场顶板暴露尺寸与水力半径关系曲线表 采场高度为 时沿矿体走向布置侧帮水力半径采场宽度 水力半径 年第 期 第 卷矿 业 工 程 表 采场高度为 时沿矿体走向布置上盘水力半径采场长度 水力半径 由图 可知，当矿体厚度小于 ，采场长度小于 时，采场顶板最大水力半径仅为 ，小于顶板稳定区的允许水力半径 ，表明采场顶板的稳定性较好；同理，根据表 、表 可以看出，采场侧帮的最大水力半径为 及上盘的最大水力半径为 ，也都小于它们稳定区的允许水力半径 和 。因此，当矿体厚度小于 时，采场沿矿体走向布置，当采场长度控制在 之内，稳定性较好，能确保开采时采场安全。垂直矿体走向布置当采场垂直矿体走向布置时，根据水力半径的计算方法，得出不同采场结构参数下采场顶板水力半径，以及采场侧帮及上盘水力半径，计算结果如图 、表 及表 所示。图 垂直矿体走向布置采场顶板暴露尺寸与水力半径关系曲线表 采场高度为 时垂直走向布置侧帮水力半径采场长度 水力半径 表 采场高度为 时垂直走向布置上盘水力半径采场宽度 水力半径 根据 稳定图解法采场顶板容许水力半径的计算，稳定区的水力半径小于 ，由图 可知，当矿体厚度小于 时，宽度为 的采场对应的水力半径为 ，表明采场顶板稳定性较好；高度为 ，矿体厚度小于 的采场侧帮最大水力半径为 ，小于侧帮稳定区允许的水力半径 ；而宽度小于 的采场上盘最大水力半径为 ，小于上盘稳定区允许的水力半径 。因此，采场垂直矿体走向布置，长度小于 采场，其顶板、侧帮及上盘的稳定性均较好，能确保开采时采场安全。矿柱宽度设计在矿岩质量分级基础上采用 稳定图解法，通过考虑采场顶板、侧帮及上盘的稳定性，得到了采场安全回采的结构参数。稳定图解法是一种经验的参数选择方法，尽管没有考虑应力环境对采场结构参数的影响，但为采场结构参数选取提供了可供参考的范围。本节基于 数值分析软件，利用 稳定图解法确定采场结构参数，确定合理的矿柱宽度。数值模拟方案根据采场布置形式选择典型的采场参数，分析采场间矿柱的宽度。对于沿矿体走向布置采场，选择的典型采场参数为 ，如图 ）所示；垂直走向布置的采场，选择的典型采场参数为 ，如图 ）所示。图 采场布置示意图 矿 业 工 程黄金运用岩体工程地质力学方法进行工程地质岩组划分、岩体结构分类，结合岩体分级结果及原有资料综合选取岩体力学参数，得到蒙库铁矿东段 以上矿体及围岩的物理力学参数，如表 所示。表 岩体的物理力学参数岩性容重（）弹性模量 泊松比抗压强度 抗拉强度 内聚力 内摩擦角（）矿体 变粒岩 矿柱是决定采场稳定状态的重要结构单元，矿柱结构破坏必然引发采场力学状态发生变化，影响整个采空区稳定性。本次数值模拟方案在固定跨度条件下，通过研究不同矿柱厚度对应的采场及矿柱的稳定性，进而确定矿柱最小安全厚度，分析方案如表 所示。表 矿柱安全厚度分析方案沿走向布置垂直走向布置长 宽（）矿柱 长 宽（）矿柱 沿走向布置采场矿柱最小安全厚度的确定位移是反应采场结构稳定性的重要参数，采场某些区域的位移过大预示着这些区域很可能发生局部失稳。不同厚度矿柱的位移云图如图 所示。当矿柱厚度为 时，矿柱的最大位移为 ；矿柱厚度为 、对应的最大位移分别为 、。这表明，随着矿柱厚度的减小，矿柱的最大位移逐渐增加，并且矿柱厚度越小，矿柱最大位移的增量越大。最大剪切应变率则可反映采场中发生剪切破坏的区域，剪切应变率大的区域则形成剪切破坏带。不图 矿柱的最大位移云图同厚度矿柱的最大剪切应变率分布云图如图 所示。矿柱厚度为时，矿柱表面绝大部分的剪切应变率处于随机分布的状态，矿柱整体剪切变形比较均匀，只有矿柱中间部位出现了剪切应变率较大的区域；矿柱厚度变为时，中间部位的剪切变形区域向矿柱顶底部扩展，但是没有形成连续的剪切变形条带；随着矿柱厚度的进一步降低，矿柱表面已经形成了连续的剪切变形条带，整个矿柱的剪切变形呈现明显的不均匀状态，矿柱容易沿着形成的剪切变形条带发生剪切破坏。图 沿走向布置采场矿柱的最大剪切应变率分布云图最大主应力和最小主应力的分布对于判断采场稳定性是至关重要的。采场的最大主应力分布云图如图 所示。由图 可知：不同矿柱采场的最大主应力大小几乎相同，分布十分相似，尽管矿柱侧帮及采场侧壁处都出现了拉应力，拉应力的最大值约为 ，但都小于围岩及矿石的抗拉强度，因而采 年第 期 第 卷矿 业 工 程 场不会出现拉伸破坏。随着矿柱厚度的不断减小，矿柱的最小主应力呈逐渐增大的趋势，如图 所示。但 矿柱采场相较于 矿柱采场的最小主应力变化不大，而矿柱采场的最小主应力却有了较为明显的增加，结合图 ）可知，矿柱采场在较高的最小主应力下发生压剪破坏的可能性增大。图 沿走向布置采场矿柱的最大主应力分布云图图 沿走向布置采场矿柱的最小主应力分布云图综上分析，矿柱采场与 矿柱采场的最大位移、最大主应力及最小主应力的大小基本相同、分布大致相似，没有表现出明显的差异。尽管 矿柱采场的最大剪切应变率相较于 矿柱采场而言有了进一步的增加，但增加的幅度不大，没有形成贯穿的剪切带，可以认为 矿柱采场仍然处于稳定状态。而 矿柱采场的最大位移和最小主应力相较于 矿柱采场有了较为明显的不同，并且 矿柱采场形成贯穿的剪切带，据此可以得出，对于 的采场而言，矿柱的最小安全厚度为 。垂直走向布置采场矿柱最小安全厚度确定垂直走向布置时采场的最大位移分布云图如图 所示，最大位移总体也呈现出随矿柱厚度变小最大位移增大的规律，并且矿柱厚度越小，矿柱最大位移的增量越大，这与沿走向布置采场的位移变化规律一致。最大剪切应变率的分布也表现出了由无序到规则的变化趋势（如图 所示），但 矿柱采场的剪切带在矿柱顶底已经大致形成，虽未贯通，但剪切破坏的可能性很大。图 垂直走向布置采场矿柱的最大位移云图图 垂直走向布置采场矿柱的最大剪切应变率分布云图 矿 业 工 程黄金采场最小主应力、最大主应力均随矿柱厚度的减小呈现出增大的趋势，并且增加的幅度较为明显，如图 、图 所示。矿柱采场的最大主应力为 ，但未达到矿柱或者围岩的最大抗拉强度，而 矿柱采场的最大主应力达到了 ，非常接近矿柱或者围岩的最大抗拉强度。结合图 ）分析可知，矿柱采场很可能发生拉伸与剪切联合破坏，处于不稳定的临界状态。图 垂直走向布置采场矿柱的最小主应力分布云图图 垂直走向布置采场矿柱的最大主应力分布云图综上分析，垂直走向布置时，矿柱位移、最大剪切应变率及采场的最大、最小主应力都随着矿柱厚度的变小而增大，但是矿柱采场的顶部出现了较大范围的剪切带，并且采场周边的最大拉应力接近于围岩或矿体的抗拉强度，在回采过程中受到其他外力扰动，很有可能发生剪切或拉伸破坏，据此可以得出，对于 的采场而言，矿柱的最小安全厚度为。结论）采用岩体结构分析系统对矿体顶板、侧帮及下盘围岩典型区域进行节理裂隙扫描测量，得到采场顶板的 组控制性节理产状：，；采场 侧 帮 的 组 控 制 性 节 理 产 状：，；采场上盘围岩的 组控制性节理产状：，。）基于 稳定图解法针对沿走向布置与垂直走向布置采场的顶板、侧帮及上盘开展稳定性分析，采场沿矿体走向布置时，矿体厚度小于 ，采场长度控制在 之内，稳定性较好；采场垂直矿体走向布置时，矿体厚度不超过 ，宽度小于 的采场，稳定性较好。）通过分析矿柱的位移、最大剪切应变率及采场的主应力等因素，确定了 种典型采场结构的矿柱最小安全厚度。对于沿走向布置的 的典型采场，矿柱最小安全厚度为 ；对于垂直走向布置的 的典型采场，矿柱最小安全厚度为。参 考 文 献 王安建，高芯蕊 中国能源与重要矿产资源需求展望 中国科学院院刊，（）：薛小蒙，莫东旭，罗佳，等 大直径深孔落矿工艺在低品位厚大矿床开采中的应用 矿业研究与开发，（）：李启月，刘恺，李夕兵 基于协同回采的深部厚大矿体分段充填采矿法 工程科学学报，（）：陈霖，黄明清，唐绍辉，等 大直径深孔空场嗣后充填法采场结构参数优化及稳定性分析 金属矿山，（）：邹平，贺超，李爱兵，等 某缓倾斜铁矿采场结构参数与开采顺序优化研究 矿冶工程，（）：罗先伟，韦军，高峰，等 深部急倾斜厚大矿体采场宽度优化 金属矿山，（）：徐帅，安龙，李元辉，等 基于 的深埋厚大矿体采场结构参数优化研究 采矿与安全工程学报，（）：刘嘉伟，黄明清，谭伟 基于拓展的 稳定图法的采场结构参数优化 有色金属工程，（）：罗来和，夏国进，王沉，等 基于 的采矿方法优选及采场 年第 期 第 卷矿 业 工 程 结构参数优化 矿冶工程，（）：赵兴东，周鑫，魏慧 三山岛金矿分段空场嗣后充填采场结构参数设计 金属矿山，（）：刘建东，解联库，曹辉 大规模充填采矿采场稳定性研究与结构参数优化 金属矿山，（）：郭进平，晏承园，张雯，等 分段矿房嗣后充填采矿采场结构参数优化研究 矿业研究与开发，（）：张小辉，姜春雨，程小辉 基于 系统的岩体质量快速评价优化创新方法及应用 有色金属（矿山部分），（）：牞 牞 牞 牞 牞 牞 牞 牗 牞 牷 牞 牷 牞 牷 牗 牘 牞 牷 牞 牘 牶 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牞 牶 牷 牷 牷 牷