《武汉工程大学学报》 2016年3期
277-282
出版日期：2016-06-22
ISSN:1674-2869
CN:42-1779/TQ
基于FLAC3D大冶铁矿矿柱回采过程的应力分析
1 引 言为了增加资源回收率，金属非金属矿山在回采后期需要对残留矿柱进行二次回采. 矿柱结构破坏必然引发采场力学状态发生变化，甚至影响整个采空区的稳定性. 矿柱的回收对于采空区处理、资源节约等方面具有重要的经济和安全意义，正确合理确定矿柱回采方式是矿柱回收过程中必须解决的关键问题. 大量工程实践案例表明，矿山采空区任一矿柱的失稳易形成多米诺连锁效应的坍塌. 了解矿山采空区群的整体稳定性，研究承受高应力下矿柱对动力扰动的响应有非常重要的实践指导意义［1-3］. 基于上述考虑，在对现场详尽调查的基础上，笔者选用动力有限元法研究湖北大冶铁矿采空区动力扰动下的响应动力特征，通过使用FLAC3D（Fast Lagrangian Analysis of Continua）软件计算，综合分析采空区矿柱负载在20 MPa、40 MPa的情况下的应力分布规律，为矿柱安全回采、充填设计提供理论依据. 2 模型的建立2.1 工程概况大冶铁矿位于湖北省黄石市铁山区. 大冶铁矿为低山-丘陵组成的山丘-盆地地形，山脉走向NWW，与构造线走向一致. 地势北高南低，低山与丘陵之间走向为NWW-SEE. 铁门坎采区浅部矿体于1984年7月结束露天开采，坑底标高-36 m. -36 m以下及挂帮矿转为地下开采. 该矿区地下开采工程由长沙黑色冶金矿山设计院设计，阶段高60 m，分段高12 m，采矿方法为无底柱分段崩落法开采. 该采区基建工程于1995年竣工投产. 到2005年底，-62 m分段及以上分段已由大冶铁矿井下车间采用无底柱分段崩落法开采完毕. 各区段分层矿柱统计：其中-62 m到-50 m区段由于充填有大量黄泥，为防止-62 m分段的黄泥下泄影响下分段的采矿，对-62 m分段底柱与点柱不予开采，故不作统计. 各分段统计矿柱如表1所示. 由表1可知，铁门坎采区-50 m到-110 m区段，顶底柱占比84.5%，点柱占比15.5%，也就是说残留矿柱主要是顶底柱；残留矿柱约77 677 m3，矿石体质量按4.12 t/m3计算，资源储量32.0万吨. 从安全回采来看，-110 m分段底柱列入下区段回收较妥，故本区段可回收矿柱约63 217 m3，资源储量为26.0万吨. 2.2 数值模型网格的建立围岩和矿体物理力学参数根据前人研究的大冶铁矿岩体物理力学参数选取［4-9］，结合工程实践经验，各计算参数取值见表2. 选取矿柱直径3 m、高12 m的竖直圆形矿柱进行计算模拟分析. 为了模拟竖直方向的地应力，模型上边界施加一竖直方向的静载（负荷），下边界施加位移约束. 为考察动力扰动对矿柱的影响，在模型上边界施加动力荷载. 矿柱的静力模型、动力模型、计算网格见图1、图2和图3. 3 模拟结果及分析3.1 设计计算方案首先，为了考察矿柱的承载能力，计算当高径比为4∶1时，矿柱在不同静压力下的应力和变形情况. 基于此，为了进行动力扰动分析，在矿柱模型顶端施加应力波荷载. 选用正弦脉冲来分布荷载（见图4）. 当扰动应力波的峰值分别取pmax=10 MPa、30 MPa时，分析计算静载大小对矿柱的影响. 3.2 矿柱应力荷载模拟分析根据工程地质，利用FLAC3D程序强大的后处理功能［10-15］，对大冶铁矿矿柱回采过程进行安全性分析，限于篇幅，这里只给出矿柱动力模拟分布图. 分别取矿柱模型高0 m、6 m和12 m处横截面中心点作为动态信息记录对象，以研究应力波作用对矿柱的应力变化. 在时间变化时，负载为pmax=10 MPa、30 MPa时最大主应力的变化情况见图5～图16. 模拟过程显示，模型单元最大主应力在负载p=20 MPa、应力波峰值pmax=10 MPa时曲线如图5、图6和图7所示；主应力时程曲线在应力波峰值pmax=30 MPa时如图8、图9和图10所示. 对比表明，模型上部单元的应力在应力波作用后变化较为稳定，模型中部和下部单元的应力变化有少许波动；随着应力波峰值的增大，矿柱模型单元的最大主应力变化不明显. 单元最大主应力时程曲线在负载p=40 MPa、应力波峰值pmax=10 MPa时如图11、图12和图13所示；单元最大主应力的时程曲线在负载p=40 MPa、应力波峰值pmax=30 MPa时如图14、图15和图16所示. 根据模型上中下部单元的应力特征，上、中、下部单元的应力在应力波作用后有较大变化；随着应力波峰值的增大，矿柱模型单元的最大主应力变化不明显. 4 结 语以大冶铁矿采空区内矿柱作为研究对象，对数值分析方案进行了设计，对采空区矿柱在动力扰动下的响应动力特征选用动力有限元法进行了综合分析. 通过对回采矿柱的动力模拟分析，可以得出如下结论：1）当矿柱负载p=20 MPa时，矿柱模型单元的应力变化总体较稳定，模型单元应力值在应力波作用前后基本保持不变，说明负载20 MPa的该矿柱在选取的动力扰动下稳定性较差. 2）在负载p=40 MPa时，在应力波作用过程中及作用之后，矿柱模型单元应力变化的波动现象十分明显，而且随着应力波峰值的增大，单元应力波动幅度也有所增大. 3）以上2种荷载情况表明，随着矿柱荷载的增大，动力扰动对矿柱模型的影响越来越明显，模型也表现出越来越不稳定. 表明矿柱随着所受应力的增大，也越易引起失稳而破坏. 由此，当为矿柱的回采进行充填设计时，要考虑到矿柱的负载作用，矿柱的负载越小，整个采场的稳定性越好.