程潮铁矿联合开采隔离矿柱合理厚度研究.pdf
程潮铁矿联合开采隔离矿柱合理厚度研究 ① 裴明松, 许梦国, 程爱平, 王 平 (武汉科技大学 资源与环境工程学院,湖北 武汉 430081) 摘 要 程潮铁矿矿体随开采水平下降而逐渐西移,导致选厂下保安矿柱的矿量逐渐增加。 为了保护地表选厂同时尽可能回收矿 石资源,程潮铁矿采用充填法回收保安矿柱和无底柱分段崩落法开采塌陷坑下矿体联合开采回收-430~-500 m 阶段矿体。 为了防 止两种采矿方法交界处的采矿活动相互影响,采用数值模拟方法,选取 15 m、20 m 和 25 m 厚度的隔离矿柱分别建立了联合开采数 值模型模拟开采,得到联合开采后地表沉降以及充填采场的安全系数,并对其进行比较。 结果表明,隔离矿柱厚度为 20 m 时,选厂 地表监测点的沉降值较小,充填采场安全系数较高,能够满足选厂和充填采场安全性要求,同时矿柱压矿量少,经济效益较优。 研 究结果能够为矿山的实际生产提供一定的指导。 关键词 隔离矿柱; 合理厚度; 地表沉降; 安全系数 中图分类号 TD853文献标识码 Adoi10.3969/ j.issn.0253-6099.2016.05.002 文章编号 0253-6099(2016)05-0005-05 Reasonable Thickness of Partition Pillar for A Combined Mining Method in Chengchao Iron Mine PEI Ming⁃song, XU Meng⁃guo, CHENG Ai⁃ping, WANG Ping (School of Resource and Environmental Engineering, Wuhan University of Science and Technology, Wuhan 430081, Hubei, China) Abstract The location of ore body moves toward west gradually with the mining going deeper, leading to safety pillar with ore reserves under the processing plant gradually increased. In order to recover the mineral resources as much as possible while protecting the safety of plant surface, a combined mining method including back filling to recover safety pillar and pillarless sublevel caving to stope orebody under collapsed pit, was adopted to recover orebody at the level -430~-500 m in Chengchao Iron Mine.Considering the interaction in the cross⁃operation area, numerical simulation model for the combined mining method with partition pillar in thickness of 15 m, 20 m and 25 m was constructed, resulting in safety factors for surface subsidence after combined mining and for the backfill stope. Based on the comparison of those results, it is found that with the partition pillar in the thickness of 20 m, the dressing plant surface is monitored to have smaller subsidence and the backfill stope has higher safety factor, meeting the safety requirement. In addition, ore quantity of partition pillar can be reduced to be small, bringing in better economic benefit. This research result can provide some guidance for the actual production on the mine. Key words partition pillar; reasonable thickness; subsidence; safety factor 在国际铁矿石价格出现断崖式下跌后,程潮铁矿 在完成集团公司的生产任务之外,也要考虑如何保证 矿山经营指标的实现。 在开采进入中深部以后,矿体 逐渐西移,随之而来的是选厂预留保安矿柱所占的资 源量越来越大。 为了保护地表选厂同时尽可能回收国 家矿产资源,在经过一系列优选之后,程潮铁矿采用充 填法回收保安矿柱、无底柱分段崩落法回收塌陷坑范 围内矿体联合开采的方法。 为了防止两种采矿方法交 界处的采矿活动相互影响,有必要在两种采矿方法衔 接区域设置隔离矿柱降低两种方法联合开采对采场的 安全性和选厂区域地表的影响。 确定隔离矿柱的合理 厚度不仅能够保证采场和地表的稳定性,同时能够尽 ①收稿日期 2016-03-25 基金项目 国家自然科学基金项目(51604195) 作者简介 裴明松(1991-),男,河南南阳人,硕士研究生,主要研究方向为采矿工艺优化及采矿方案评价。 通讯作者 程爱平(1986-),男,湖北仙桃人,博士,讲师, 主要研究方向为采矿工艺及矿压监测。 第 36 卷第 5 期 2016 年 10 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.36 №5 October 2016 万方数据 可能多地回收矿石,增加经济效益。 1 工程概况 隔离矿柱合理厚度的确定是一个复杂的安全、经 济多目标系统问题,受采场围岩强度、联合开采回采顺 序、采场安全系数、地表稳定性以及出矿量制约。 在同 一开采水平采用联合开采生产组织管理难度较大,对 设备和人员的管理要求较高,因此只有部分矿体产状 变化较大的矿山在实际生产中采用,且基本不涉及两 种或多种采矿方法的衔接处隔离矿柱的布置问 题[1-2]。 目前国内隔离矿柱的研究多存在于露天转地 下开采以及崩落转充填开采,隔离矿柱的布置多为水 平布置方式[3-4]。 因此,程潮铁矿采用两种采矿方法 联合开采隔离矿柱的合理厚度研究没有可参考的实际 矿山案例。 本文以隔离矿柱的合理厚度作为研究目 标,以充填采场及隔离矿柱的安全系数以及地表沉降 作为决策要素,研究在不同隔离矿柱厚度下充填采场 的安全以及地表沉降的规律,考虑出矿量以及经济效 益,寻找能够满足充填采场的安全系数以及地表稳定 性的隔离矿柱最小厚度,为矿山生产提供一定的参考。 由于-430~-500 m 阶段矿体是程潮铁矿下一阶 段的重点开采对象,-482~-500 m 分段为该阶段的最 后一个分段,保安矿柱的矿量以及崩落法采场覆盖层 的厚度也很大,故将该阶段作为联合开采隔离矿柱厚 度的研究对象,具有代表性。 目前,矿山实际生产中保 安矿柱岩层移动角为 62,根据矿山实际情况,选取典 型剖面,建立开采示意图如图 1 所示。 图 1 矿体开采示意 根据程潮铁矿充填法回收保安矿柱的试验段的相 关经验,将矿块划分为矿房和矿柱,其宽度分别为 13 m 和 7 m,采用充填法采矿时考虑到地表变形不能过大, 仅回收矿房,留矿柱与充填体支撑上部岩体荷载。 由于 隔离矿柱同时作为附近矿块的矿柱使用,根据专家建 议,选择隔离矿柱厚度从 15 m 开始,按照 5 m 的梯度增 加,进行研究的隔离矿柱厚度分别为 15,20 和 25 m。 2 联合开采数值模拟 2.1 模型的建立 数值分析是分析岩土工程复杂问题的常用有效研 究手段,通过数值模拟矿体开采过程能够直观分析采 矿过程中的问题,从而达到研究目的[5]。 根据矿区实 际情况,由图 1 所示矿体开采示意图建立数值模型,如 图 2 所示。 利用 Ansys workbench 并导入 FLAC3D则可 以避免 FLAC3D建立不规则模型的缺点[6]。 利用 Ansys workbench 建立三维模型并划分网格,模型的长宽 高为 2 000 m 10 m 1 200 m, 网格为四面体网格,总 计单元 84 199 个,节点 27 223 个。 将上述模型导入 FLAC3D中进行数值计算。 为了研究在不利条件下采 场与隔离矿柱的稳定性情况,模拟中充填法和崩落法 均一次开采完毕。 图 2 矿体开采数值模型 2.2 模型初始条件及边界条件 自重应力以及构造应力相互叠加形成模型边界的 应力约束,构造应力中,最大主应力 σ1为水平方向, σ1= 1.27γH ,最小主应力 σ3基本与矿体走向垂直, σ3=0.44γH,次主应力 σ2= 12.7 MPa。 模型的自重应 力及其水平分量可以通过下列算式[7]获得 σv= γH(1) σH= μ 1 - μγH (2) 式中 σv为自重应力;γ 为岩体容重;H 为岩体赋存深 度;σH为自重应力水平分量;μ 为覆岩的泊松比。 在模型的左右侧面施加水平位移约束,在其前后 侧面施加固定位移约束,底部施加竖直约束。 由于模 型的上表面为实际地表,故不施加任何约束。 2.3 物理力学参数 根据矿区地质资料显示,在矿体上盘围岩以大理 岩和矽卡岩为主,下盘以花岗岩、闪长岩以及矽卡岩为 主。 围岩及矿体的物理力学参数通过现场采样经实验 室中岩石力学实验获得。 数值计算中岩体的物理力学 参数经过实验室岩样的参数折减获得[8]。 如表 1 所示。 6矿 冶 工 程第 36 卷 万方数据 表 1 数值模型中的岩体力学参数 岩性 容重 / (kNm -3 ) 弹性模量 / GPa 泊松比 内摩擦角 / () 粘聚力 / MPa 大理岩27.34.50.2131.36.48 花岗岩26.34.80.2336.23.50 磁铁矿36.55.20.2833.65.32 崩落散体18.60.60.3320.00.60 充填体22.10.60.2928.00.60 矽卡岩27.43.40.3220.03.50 闪长岩25.94.50.2834.05.12 2.4 选择本构模型 由于上述围岩及矿体均属于弹塑性材料,故选取 摩尔⁃库伦屈服准则进行计算[9] fs = σ 1 - σ 3N - 2c N (3) 式中 σ1为最大主应力;σ3为最小主应力;N 的取值 为N= (1+sinφ) / (1-sinφ);c 为粘聚力;φ 为内摩擦 角。 一般认为,在 fs<0 时材料处于剪切屈服阶段,在 fs>0 时材料处于拉伸屈服阶段。 3 数值模拟结果分析 通过对保安矿柱及崩落开采矿体运用两种采矿方 法进行联合开采过程进行数值模拟,在选厂区域地表 自西向东每隔 30 m 依次布设 18 个监测点来观测地表 的变形情况。 在模拟矿体开采过程中记录各个监测点 的沉降情况,分析在不同隔离矿柱情况下下部采矿生产 活动对选厂区域地表沉降的影响。 同时在记录开采结 束后充填采场的应力分布情况并得到其安全系数,从而 分析在隔离矿柱厚度的不同对采场稳定性的影响。 3.1 隔离矿柱厚度对地表沉降影响分析 利用 FLAC3D软件模拟-482~ -500 m 分段开采, 得到两种方法同时联合开采后模拟断面垂直位移情况 如图 3 所示。 图 3 不同厚度隔离矿柱下开采模型垂直位移情况 从图 3 中可以看到,在采用两种采矿方法同时开 采结束后,选厂区域的地表均出现了不同程度的沉降, 且影响的区域面积较开采区域面积大,自塌陷坑西侧 边缘往西地表沉降值逐渐减小。 塌陷坑区域的沉降值 远大于选厂区域的沉降,这与塌陷坑下部矿体采用无 底柱分段崩落法开采有关。 选厂区域地表监测点沉降 结果如表 2 所示。 根据表 2 得到监测点沉降值与距离 塌陷坑的距离关系曲线如图 4 所示。 由表 2 及图 4 可 以得到,选厂区域各监测点的垂直位移值随其与塌陷 坑距离增大而逐渐减小。 监测点的垂直位移值整体与 两种采矿方法之间隔离矿柱的厚度负相关。 对比 15 m 厚度隔离矿柱下与其他两种厚度下监测点的垂直位移 情况可以看到,在距离塌陷坑 0~110 m 之间的监测点 垂直位移值明显增大,增大量为 1.01~5.23 mm。 表明 在该区域隔离矿柱为 15 m 厚度时对地表的影响明显 较大。 对比隔离矿柱厚度为 20 m 和 25 m 时监测点的 垂直位移值可以看到,随着隔离矿柱厚度由 20 m 增加 到 25 m,地表垂直位移值减小量较小,近似减小 0.3~ 0.6 mm。 此外,监测点垂直位移结果还表明,隔离矿柱 厚度为 20 m 时,选厂区域与塌陷坑西侧边缘交界位 表 2 不同厚度隔离矿柱下选厂区域地表监测点沉降值 点号 距塌陷坑距离 / m 垂直位移/ mm 15 m20 m25 m 1#0 -35.71 -31.01 -30.48 2#20 -30.05 -26.86 -26.31 3#50 -24.56 -21.93 -21.41 4#80 -19.98 -17.92 -17.43 5#110 -16.26 -15.24 -14.78 6#140 -13.75 -13.02 -12.60 7#170 -10.79 -10.27 -9.88 8#200 -8.82 -8.43 -8.06 9#230 -7.49 -7.15 -6.78 10#260 -6.30 -6.00 -5.64 11#290 -5.28 -5.01 -4.66 12#320 -4.56 -4.30 -3.97 13#350 -4.10 -3.85 -3.55 14#380 -3.77 -3.53 -3.24 15#410 -3.44 -3.22 -2.93 16#440 -3.09 -2.88 -2.59 17#470 -2.77 -2.57 -2.28 18#500 -2.56 -2.37 -2.08 置的监测点垂直位移值为 31.01 mm;隔离矿柱厚度为 25 m 时该值为 30.48 mm,但是当隔离矿柱厚度为 15 m 时该值为 35.71 mm,较前两者分别增大 4.7 mm 和 5.23 mm,在深部开采仍然沿用两种方法同时开采 7第 5 期裴明松等 程潮铁矿联合开采隔离矿柱合理厚度研究 万方数据 02223m 0 -5 -10 -15 -20 -25 -30 -35 -40 0200100300400500 0*DAmm 15 m 20 m 25 m 图 4 监测点沉降值与其距塌陷坑距离的关系 时,已经可能危及到塌陷坑附近矿区公路以及选厂围 墙附近区域的安全。 此外,在研究剖面上,选厂主要构 (建)筑物距离塌陷坑边缘的距离均在 170 m 以上,地 表的垂直位移值为 2 ~ 10 mm,基本不会对选厂的构 (建)筑物产生影响。 综上所述,从不同厚度隔离矿柱 引起地表垂直位移值的角度而言,选取 20 m 及以上的 隔离矿柱厚度比较合适,能够保证选厂主要构(建)筑 物、矿区公路的安全。 但是,在开采进入深部以后,主要构(建)筑物所 处范围内的地表垂直位移值可能会逐渐增大从而影响 地表构(建)筑物的安全,因此应该对该区域加强监测 并根据国家规范建筑物、水体、铁路及主要井巷煤柱 留设与压煤开采规程对砖混结构建筑物破坏等级的 损坏分类,对该砖混结构建筑损坏等级进行评级,及时 采取相应的措施进行处理。 3.2 隔离矿柱厚度对采场稳定性的影响 根据数值计算结果,将不同厚度隔离矿柱下充填 采场附近区域围岩的弹塑性区域标示出来,如图 5 所 示。 由图 5 可以看到,厚度为 15 m 隔离矿柱中塑性区 域面积最大,且在进行-482~-500 m 分段开采时塑性 区域已经到达上一阶段。 厚度为 20 m 隔离矿柱塑性 区域面积次之,厚度为 25 m 隔离矿柱最小,且在本阶 段。 但是厚度为 25 m 隔离矿柱下采场周围的隔离矿 柱围岩的塑性区域与厚度为 20 m 时减小面积不大。 由于模拟中未考虑充填采场的支护以及顶底柱,因此, 在不同厚度隔离矿柱下充填采场均为 shear⁃p,即在开 采模拟过程中发生剪切破坏。 同时,从充填采场及隔离矿柱周围的安全系 数[10-12]来看,厚度为 15 m 的隔离矿柱中安全系数在 1.2~2.2 之间,但是在与无底柱分段崩落法交界面之间 存在多处区域安全系数小于1;采场安全系数在1.8~2.2 之间,采场上部已充填完成的区域安全系数在 1.0~1.8 之间。 厚度为 20m 的隔离矿柱中安全系数在 1.8~2.2 之间,但是在与无底柱分段崩落法交界面之间存在个别 区域安全系数小于 1;采场安全系数在 1.8~2.2 之间,采 场上部已充填完成的区域安全系数在 1.2~1.8 之间。 厚度为 25 m 的隔离矿柱中安全系数为 1.8~2.6,与无底 柱分段崩落法交界面之间存在极少数区域安全系数小 于 1;采场安全系数在 1.8~2.4 之间,采场上部已充填 完成的区域安全系数在 1.2~1.8 之间。 因此,就隔离 矿柱及充填采场的安全系数而言,隔离厚度为 20 m 和 25 m 均可满足要求,并有一定的余量。 如图 6 所示。 图 5 不同厚度隔离矿柱下充填采场及隔离矿柱塑性区域 图 6 不同厚度隔离矿柱下充填采场及隔离矿柱安全系数 8矿 冶 工 程第 36 卷 万方数据 根据模拟开采结果及上述分析,就隔离矿柱厚度对 采场稳定性而言,20 m 厚度的隔离矿柱即可满足要求。 3.3 分 析 综上数值模拟分析结果可以看到,当隔离矿柱厚 度为 15 m 时,选厂区域地表监测点的垂直位移值较 大,隔离矿柱中塑性区域面积较大,充填采场和隔离矿 柱的安全系数也较其他两种厚度方案低,一旦发生较 大扰动,可能危及采矿生产的安全。 当隔离矿柱厚度 为 20 m 时,选厂区域地表监测点的垂直位移值不大, 隔离矿柱中塑性区域面积较小,充填采场及隔离矿柱 的安全系数较高,能够保证采矿生产的安全。 当隔离 矿柱厚度为 25 m 时,选产区域地表监测点的垂直位移 值很小,隔离矿柱中塑性区域面积较小,充填采场和隔 离矿柱的安全系数高,采矿生产活动的安全能够保证。 一般而言,在能够满足安全要求的前提下,选取经 济上最优的方案。 在隔离矿柱厚度为 20 m 和 25 m 时,就回收矿产资源角度而言,隔离矿柱为 20 m 可多 回收铁矿石 26.24 万吨,可为企业多创造经济效益。 因此程潮铁矿回收保安矿柱的充填法和回收塌陷坑下 矿体的无底柱分段崩落法之间衔接区域预留 20 m 厚 度的隔离矿柱为最优。 4 结 论 程潮铁矿采用无底柱分段崩落法回收塌陷区下矿 体和采用充填法回收选厂下保安矿柱时会面临两种采 矿方法衔接区域相互影响的问题,为了减弱该影响需 要在衔接区域设置隔离矿柱,隔离矿柱的合理厚度关 系着地表选厂以及充填采场的安全,经过对不同厚度 隔离矿柱下地表垂直位移以及充填采场安全性的数值 模拟分析,可以得到以下结论 1) 隔离矿柱厚度为 20 m 选厂区域地表的垂直位 移值较小,且充填采场的安全性能够得到保证,与厚度 为 25 m 的隔离矿柱相比可提高出矿量,增加经济效 益,因此应将隔离矿柱的合理厚度设置为 20 m。 2) 在隔离矿柱厚度为 20 m 时充填采场及隔离矿 柱的塑性区域较小,隔离矿柱中安全系数在 1.8~2.2 之间,但是在与无底柱分段崩落法交界面之间存在个 别区域安全系数小于 1;采场安全系数在 1.8~2.2 之 间,采场上部已充填完成的区域安全系数在 1.2~1.8 之间。 能够满足安全要求,但是还应该合理安排出矿 计划,防止在衔接区域单次爆破量过大产生强地震冲 击波对隔离矿柱以及充填采场的安全性造成威胁。 3) 本文采用有限元的方法进行数值模拟,得到的 结果与实际可能会有一定的差异。 由于程潮铁矿两种 采矿方法联合开采衔接区域预留的隔离矿柱类型可借 鉴的经验模型较少,因此在实际运用时可将本研究作 为借鉴,但还应该结合采矿实际情况对隔离矿柱的厚度 进行更合理地取值,以期达到更经济、安全的结果。 参考文献 [1] 赵忠权,陈红刚. 同一采场多种采矿方法的联合应用[J]. 中国锰 业, 2011,29(2)43-46. [2] 朱先艳,杨庆雨,李兴平. 浅眼留矿与向上水平分层充填法联合开 采试验与应用[J]. 金属矿山, 2009,11(5)199-203. [3] 刘洪磊,杨天鸿,许宏亮,等. 桓仁铅锌矿隔离矿柱留设方案模拟 研究[J]. 地下空间与工程学报, 2012,8(4)785-790. [4] 曾令义,潘 东. 崩落法转充填法隔离层安全厚度的确定[J]. 矿 冶工程, 2015,35(4)17-20. [5] 陈育民,徐鼎平. FLAC/ FLAC3D基础与工程实例[M]. 北京中国 水利水电出版社, 2009. [6] 廖秋林,曾钱帮,刘 彤,等. 基于 ANSYS 平台复杂地质体 FLAC3D模型的自动生成[J]. 岩石力学与工程学报, 2005,24 (6)1010-1013. [7] 姚高辉. 金属矿山深部开采岩爆预测及工程应用研究[D]. 武汉 武汉科技大学资源与环境工程学院, 2008. [8] 陈顺满,许梦国,王 平,等. 某铁矿挂帮矿开采对地表拦洪坝的 影响[J]. 金属矿山, 2014(2)150-153. [9] 徐文彬,宋卫东,杜建华,等. 崩落法转阶段嗣后充填法采场稳定 性分析[J]. 北京科技大学学报, 2013,35(4)415-421. [10] 杨宇江,庄文广,王照亚,等. 基于强度折减法的地下采场稳定性 分析[J]. 东北大学学报(自然科学版), 2011,32(6)864-867. [11] 赵兴东. 谦比希深部开采隔离矿柱稳定性分析[J]. 岩石力学与 工程学报, 2010,29(s1)2616-2622. [12] 江 权,冯夏庭,向天兵. 基于强度折减原理的地下洞室群整体 安全系数计算方法探讨[J]. 岩土力学, 2009,30(8)2483- 2488. (上接第 4 页) [2] 黄 牧. 太平洋深海沉积物稀土元素地球化学特征及资源潜力初 步研究[D]. 广州国家海洋局第一海洋研究所, 2013. [3] 张 鑫,栾振东,阎 军,等. 深海沉积物超长取样系统研究进展 [J]. 海洋地质前沿, 2012,28(12)40-45. [4] 任江波,姚会强,朱克超,等. 稀土元素及钇在东太平洋 CC 区深 海泥中的富集特征与机制[J]. 地学前缘, 2015,22(4)200-211. [5] 朱坤杰,王金莲,邓希光,等. 中⁃西太平洋海盆浅层沉积物的物理 力学性质初探[J]. 土工基础, 2015,29(3)149-152. [6] 魏 巍. 南海中沙天然气水合物资源远景区海底沉积物的物理力 学性质研究[J]. 海岸工程, 2006,25(3)33-38. [7] 王 凯,吴建政,安永宁,等. 渤海湾北部表层沉积物的物理力学 性质[J]. 海洋地质前沿, 2011,27(1)14-18. [8] 宋连清. 大洋多金属结核矿区沉积物土工性质[J]. 海洋学报中 文版, 1999,21(6)47-54. 9第 5 期裴明松等 程潮铁矿联合开采隔离矿柱合理厚度研究 万方数据