厚大破碎矿体采场矿岩冒落规律及影响因素研究_李晓斌.pdf

收稿日期2019-12-24 作者简介李晓斌 （1980） ， 男， 高级工程师， 硕士。通讯作者勒治华 （1992） ， 男， 博士研究生。 厚大破碎矿体采场矿岩冒落规律及影响因素研究 李晓斌 1 勒治华 2 樊硕 2 李新成 11 （1. 中冶沈勘秦皇岛工程设计研究总院有限公司， 河北 秦皇岛066004； 2. 东北大学资源与土木工程学院， 辽宁 沈阳 110819） 摘要自然崩落法矿岩冒落是拉底卸荷、 矿岩冒落、 直至形成平衡拱的过程， 其本质则是不断重复该过程。矿 岩冒落规律受其地应力、 采场结构参数、 矿岩物理力学性质等地质因素和工程因素的影响， 目前的研究多是假定为二 维模型， 侧重于采用理论方法或数值模拟分析冒落过程和冒落拱形态， 较少涉及三维空间的冒落过程研究。根据矿岩 冒落过程的主要特征参数， 建立了描述该过程力学特征的本构模型， 并采用FLAC3D实现数值模拟， 模拟研究了矿体赋 存深度、 侧压力系数、 拉底形态及采场尺寸参数对矿岩冒落规律的影响， 研究结论有助于指导现场实践。 关键词地应力采矿结构参数冒落高度冒落体积 中图分类号TD853.36文献标志码A文章编号1001-1250 （2020） -06-172-05 DOI10.19614/ki.jsks.202006027 The Caving Process and Influential Factors in the Stope of Thick and Large Fractured Ore Body Li Xiaobin1Jin Zhihua2Fan Shuo2Li Xincheng12 （1. Shenkan Qinhuangdao General Engineering DesignResearch Institute Corporation， MCC， Qinhuangdao 066004， Chi- na； 2. School of Resources and Civil Engineering， Northeastern University， Shenyang 110819， China） AbstractOre rock caving in the natural caving is a process of undercutting and unloading，ore-rock caving and even ing a balance arch． Its essence is to repeat the process continuously． The regularity of ore rock caving is affected by geological and engineering factors such as ground stress，stope structural parameters，physical and mechanical properties of ore rock． At present，most of the research is assumed to be two-dimensional model，focusing on the theoretical or numerical simulation analysis of the caving process and the shape of the caving arch， less involving the research of the caving process in three-dimensional space． According to the main characteristic parameters of the ore rock caving process，a consti⁃ tutive model is established to describe the mechanical characteristics of the process， and FLAC3Dis used to realize the numer⁃ ical simulation on the effects of the ore body occurrence depth，side pressure coefficient，bottom and stope size parame⁃ ters on the regularity of ore rock caving． The research conclusions are helpful to guide the field practice． KeywordsGeostress， Structural parameters of stope， Falling height， Falling volume 我国金属矿产资源的特点是规模较小、 开采成 本较高、 市场核心竞争力弱。研究结果表明， 大型金 属矿山采矿成本的30～60与岩体破碎有关。通过 研究矿岩冒落的规律， 调控采动应力诱导矿岩自然 冒落， 实现利用地压破碎岩体， 可大幅度节约采准工 程量与落矿费用， 由此大幅度降低生产成本并提高 生产效率， 从而为国内金属矿山高效开采提供一条 新的思路 ［1］。 近年来， 人们逐渐认识到构造应力对采场上覆 岩层移动起重要作用。王军等 ［2］利用FLAC3D数值模 拟软件、 理论分析等方法对不同构造应力场中采空区 上覆岩层的破坏规律进行研究， 结果表明随侧向应力 的增加， 上覆岩层塑性区高度和体积先减小后增大； 范广勤等 ［3］研究了在地应力较高的情况下矿岩冒落 区域的轴比 （或纵横比） 与原始岩体的应力场的侧压 力系数成正相关关系， 研究发现， 拉底的大小、 形状和 布置， 以及切帮和预裂工程的存在， 都直接影响岩体 的崩落。徐腊明 ［4］利用三维有限元数值模拟程序 ADINA模拟了程潮铁矿拉底和割帮工艺作用下矿体 的应力分布， 在此基础上， 优选出合理的拉底方案。 总第 528 期 2020 年第 6 期 金属矿山 METAL MINE Series No． 528 June 2020 172 ChaoXing 影响矿岩可冒性的因素有很多， 主要包括矿岩 的强度、 岩体结构面特征、 地应力的相对大小、 采场 结构参数等， 本项目重点研究地应力与采场结构参 数对矿岩冒落的影响。为实现借助拉底等工程因素 诱导矿岩应力发生改变， 利用地压破碎岩体的目的， 就需要深入研究 “破碎厚大矿体采场矿岩冒落规律” 这一科学问题， 这对于合理确定诱导冒落工程， 充分 利用地压， 实现低成本、 高效率采矿， 具有重要的理 论价值和实际意义。 1矿岩冒落数值模型 1. 1矿岩冒落的FLAC3D系统 FLAC3D是基于连续介质力学理论的数值分析方 法， 虽然岩体中会产生塑性区， 但岩体不会发生冒 落， 不能模拟冒落的整个过程。为了模拟矿岩的冒 落过程， 引入了人为的假定。在FLAC3D数值模型中， 如果数值网格的单元位于拉底空间的上部塑性区 中， 且该单元中存在拉应力， 则认为该单元在重力的 作用下可以冒落， 利用FLAC3D的建模命令， 该单元被 杀死 （ “NULL” 掉） ， 被NULL掉单元的模型又处于非 平衡状态， 再次进行迭代计算， 计算平衡后如果塑性 区内依然有拉应力单元， 则重复该过程， 直至模型中 不再出现拉应力区、 平衡为止， 计算停止。以此来近 似模拟矿岩冒落过程。在计算过程中很难做到塑性 区内完全不存在拉应力， 所以在塑性区内基本无拉 应力即可将计算停止。具体计算步骤如图1。 1. 2矿岩冒落本构关系 研究表明， 利用自然崩落法开采矿石时， 矿岩冒 落以拉伸破坏为主， 因此采用拉伸准则描述冒落矿 岩。实际上， 利用自然崩落法开采矿岩体质量均稍 差， 在冒落过程中采场边帮围岩也会发生剪切破坏， 但由于堆积在采场内的冒落矿岩散体支撑作用， 往 往导致剪切片帮破坏并不明显。但采用FLAC3D模拟 时并不能考虑该作用， 因此选用基于摩尔库伦准 则的理想弹塑性模型描述采场围岩， 围岩可以发生 屈服， 但并不降低承载力 ［7］。在两个准则中， 拉伸准 则优先。 实际上， 到目前为止， 还没有一种数值模拟方 法， 可以完整地模拟矿岩冒落、 堆积、 并对采场边帮 起到支撑作用的全过程。本研究利用FLAC3D模拟矿 岩冒落也是在一些假定条件下实现的。 （1） 拉伸破坏准则 σ 3≥ σt, （1） 式中，σ3、σt分别为岩体受到的拉应力与岩体的单轴 抗拉强度。 （2） 摩尔库伦破坏准则 fs σ1-σ3 1 sinφ 1 -sinφ -2c 1 sinφ 1 -sinφ ,（2） 式中，σ1、σ3分别是最大和最小主应力；c、φ分别为黏 聚力和内摩擦角。当fs＞ 0时， 材料将发生剪切破坏。 1. 3数值计算模型与参数 由冒落拱理论可知， 矿体埋藏深度、 侧应力以及 冒落拱跨度等因素是影响矿岩冒落高度的重要因素。 基于上述FLAC3D的数值模拟方法和矿岩本构关系， 对不同埋藏深度、 不同侧压力系数、 拉底形态及采场 尺寸参数条件下的矿岩冒落规律进行了数值模拟研 究。各种工况条件的数值计算模型略有不同， 现以矿 体埋藏深度对矿岩冒落的影响为例， 简要介绍其数值 计算模型以及矿、 岩参数。建立三维数值模型如图2 所示， 数值模型长300 m， 宽150 m， 高350 m。模型共 分3个部分， 分组1为围岩部分、 分组2为矿体部分， 分组 3为拉底部分。矿体部分长 100 m， 宽 50 m， 高 100 m。矿体部分的底部距离模型底部50 m， 保证其 不受边界效应的影响。岩体力学参数如表1。 李晓斌等 厚大破碎矿体采场矿岩冒落规律及影响因素研究2020年第6期 173 ChaoXing 2数值模拟方案及边界条件设计 数值模拟分析了矿体不同埋藏深度、 不同侧压 力系数、 拉底形态及采场尺寸参数4种工况条件下矿 岩冒落规律。各工况条件下的模拟方案、 边界条件 和拉底尺寸等如表2。 数值计算过程中矿岩冒落达到稳定状态判定 以矿体埋深700 m计算分析为例， 如图3所示， 过拉底 区域中心点， 沿矿体长度方向的截面， 可以看出， 在 将塑性区开挖3次之后形成的冒落拱， 拱线的上方基 本无拉应力， 可以判定冒落拱达到稳定状态。其他 工况条件下矿岩冒落是否达到稳定状态判定完全相 同。 3矿岩冒落规律分析 3. 1 矿体埋藏深度对采场矿岩冒落的影响 各埋藏深度矿体拉底后矿体冒落部分如图4所 示， 其中下部区域为拉底部分， 上部区域为塑性区部 分。 图5为矿岩冒落高度和冒落体积随矿体埋藏深 度变化曲线。可知， 矿岩的冒落高度和冒落体积随 着埋藏深度的增加而减小， 但减小幅度均较小。因 此， 在实际生产中， 对于厚大破碎矿体， 小埋深反而 有利于矿岩的冒落； 矿体埋深较大时， 可通过割帮等 诱导工程使矿体水平应力减小， 促进冒落拱的发育。 同时矿体埋深较大时应注意冒落过程的片帮现象， 制定合理的拉底范围， 以免矿体开采过程中上下盘 围岩的片帮， 造成矿石贫化率增加。 3. 2侧压力系数对采场矿岩冒落的影响 不同侧压力系数下的矿体冒落部分如图6所示， 图中下部为拉底区域， 上部为冒落部分。发现侧压 力系数为0.4时冒落拱发育最好， 矿体冒落高度与冒 落体积在5个方案中最大。 矿体的冒落高度和冒落体积随侧压力系数变化 如图7所示， 矿岩的冒落高度和冒落体积的变化趋势 基本一致。当侧应力系数从0.4增加到1时， 矿体的 冒落高度仅减小了15.7， 但是矿体的冒落体积却减 少了43， 表现出在这一过程中冒落拱的发育变差。 金属矿山2020年第6期总第528期 174 ChaoXing 侧压力系数超过1后， 冒落高度保持稳定， 冒落体积 出现缓慢增长。 在实际生产中， 如果矿体的侧应力较大， 拉底完 成后片帮问题会较为严重， 应合理规划拉底的范围。 因此， 若使上赋矿岩冒落得更加充分， 可以通过切帮 等辅助工序来尽可能降低矿块的侧应力， 从而达到 使冒落高度和冒落体积增加的目的。 3. 3拉底长宽比对采场矿岩冒落的影响 图8为不同拉底长宽比条件下矿体冒落部分， 下 部为拉底区域， 上部为上部塑性区。由塑性区分布 图可知， 拉底形状不同， 各冒落拱的形状差异较大， 冒落高度随着拉底长宽比的增大而逐渐减小。 图9为矿体冒落高度和冒落体积随拉底长宽比 变化曲线， 矿体的冒落高度和冒落体积随拉底长宽 比增大均呈现逐渐减小趋势。拉底形状为正方形时 （长宽比为1） ， 冒落高度和冒落体积达到最大， 然而， 在拉底长宽比为2时， 冒落体积有小幅度的增加。因 此， 在拉底设计时应考虑拉底形状为正方形 （长宽比 为1） ， 在地质条件或开采条件的约束下， 拉底的长宽 比最好不大于2。 3. 4采场尺寸参数对矿岩冒落的影响 如图10所示， 随着拉底跨度的增加， 采场冒落高 度在前期增长较快， 拉底的跨度接近拉底宽度时， 冒 落高度的增长趋势开始减缓， 随着拉底跨度的持续 增加， 矿体的冒落高度保持稳定。拉底跨度较小时， 不同宽度的矿体冒落高度较为接近， 随着拉底跨度 的增加， 30 m宽度矿体的冒落高度增长速度首先降 低下来， 50 m和70 m宽度的矿体的冒落高度增长速 率在拉底跨度不断增加的过程中也相继降低， 最后 都趋于稳定。图11为不同矿体宽度的矿岩冒落高度 与采场拉底跨度比值 （高跨比） 随拉底跨度变化曲 线， 从图中可以看出， 拉底跨度相同时， 高跨比随矿 体宽度的增加而增加。相同宽度的采场， 高跨比随 拉底跨度的增加而减小。 4结论 以FLAC3D数值模拟软件为工具， 研究矿体在不 同地质条件和结构参数条件下的冒落规律， 利用所 得到的冒落规律指导生产。主要结论如下 李晓斌等 厚大破碎矿体采场矿岩冒落规律及影响因素研究2020年第6期 175 ChaoXing ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ （1） 矿体在不同地质赋存条件下呈现不同的冒 落规律。拉底条件相同， 埋藏深度增加矿岩冒落高 度， 冒落体积略有减小， 但总体变化不大； 冒落高度 随侧压力系数的增加呈现先减小后趋于平衡， 冒落 体积先减小后略有增加， 侧应力系数为1时， 冒落高 度和冒落体积最小。 （2） 拉底面积相同， 拉底长宽比越大， 矿体冒落 高度和冒落体积越小。因此在设计拉底形状时， 矿 体拉底的长宽比应最好小于2， 保证拉底后有足够的 冒落高度和冒落体积， 以减小冒落辅助工程的工作 量。 （3） 采场尺寸确定后， 冒落高度随拉底长度的增 加而增大， 增长速度随拉底长度的增加而逐渐减小， 最后冒落高度趋于平衡。如果是厚大矿体， 矿块尺 寸设计时适当增大矿块的宽度以增加冒落高度和冒 落体积。 参 考 文 献 何荣兴， 任凤玉， 谭宝会， 等．论诱导冒落与自然崩落 ［J］ ．金属矿 山， 2017 （3） 9-14． He Rongxing， Ren Fengyu， Tan Baohui， et al． Discussion on in⁃ duced caving and block caving ［J］ ．Metal Mine， 2017 （3） 9-14． 王军， 宁建国， 刘学生， 等． 构造应力下采空区上覆岩层破坏 规律数值模拟 ［J］ ． 煤矿安全， 2017， 48 （1） 183-186． Wang Jun， Ning Jianguo， Liu Xuesheng， et al． Numerical simula⁃ tion of failure law of overlying strata in goaf under tectonic stress ［J］ ． Safety in Coal Mines， 2017， 48 （1） 183-186． 范广勤， 张千里． 高地应力地下工程仿围岩自然崩落外廓特征 的洞形结构模拟试验研究 ［C］ //第15届全国结构工程学术会议 论文集 （第Ⅱ册） ．焦作 ［s．n．］ ， 2006 377-384． Fan Guangqin， Zhang Qianli． Study on simulation test of cave- shaped structure of underground engineering with high in-situ stress and imitation of surrounding rock′s natural collapse outline ［C］ //Proceedings of the 15th National Conference on Structural En⁃ gineering（The 2nd book） ． Jiaozuo［s．n．］ ， 2006 377-384． 徐腊明． 自然崩落法拉底上部矿体应力分布有限元分析 ［J］ ． 金 属矿山， 1996 （2） 11-15． Xu Laming． Analysls of orebody stress distribution over undercut dreft for natural block caving by finite element ［J］ ． Metal Mine， 1996 （2） 11-15． 王利， 陈亚娟， 高谦． 拉底诱导的岩体崩落过程区数值模拟 ［J］ ． 有色金属矿山部分， 2009， 61 （3） 44-48． Wang Li， Chen Yajuan， Gao Qian． Numerical on caving process zone induced by undeercut ［J］ ．Nonferrous MetalsMining Section， 2009， 61 （3） 44-48． 金川集团有限公司， 长沙矿山研究院． 自然崩落法矿岩崩落规 律研究 （研究报告之四） ［R］ ． 长沙长沙矿山研究院， 2006． Jinchuan Co． Ltd．， Changsha Institute of Mining Research． Study on the Law of Natural Caving Ore and Rock Caving （Research Report 4） ［R］ ． Changsha Changsha Institute of Mining Research， 2006． Hoek E，Brown E T． Practical estimates of rock mass strength［J］ ． International Journal of Rock Mechanics Mining Sciences，1997， 34 （8） 1165-1186． 金属矿山2020年第6期总第528期 176 ChaoXing