露采平台安全厚度与隐伏采空区跨度关系的数值模拟研究

康伟权， 袁海平，2， 吴贤振， 刘建伟， 尹丽冰

（1.江西理工大学资源与环境工程学院，江西 赣州 341000；2.合肥工业大学土木与水利工程学院，合肥 230009）

0 引 言

对于地下转露天开采的矿山，地采遗留下的隐伏采空区是露天开采工作的重大隐患[1]，为了保证矿山持续发展和安全[2]开采，找出露采平台安全厚度与隐伏采空区跨度之间的规律[3]愈发重要，尤其对于危险采空区，难以实行地下实时安全监测[4]，数值模拟[5]可作为分析露采平台安全厚度的一种有效手段.薛涛等[6]从应力和位移多角度，考虑空区与相邻空区相互之间的影响作用，分析开采境界下空区的稳定性状态；李地元等[7]综合采用结构力学梁法、鲁佩涅依特理论估算法等多种方法计算，并对理论计算结果和三维数值模拟分析结果进行比较，得出露采平台安全厚度与隐伏采空区跨度的对应关系，等等[8].国内外对于露采平台安全厚度与隐伏采空区[9]跨度关系的相关研究并不多，了解隐伏空区跨度与露采平台安全厚度的关系对实际工程具有一定的指导意义.

以珠江某钨钼矿为例，运用FLAC3D软件[10]进行数值模拟稳定性分析，得出露采平台临界安全厚度隐伏空区跨度之间的对应关系，提出了采空区跨度与露采平台安全厚度对应线性拟合关系.

1 岩石物理力学参数测定

矿区内的岩体主要是在燕山旋回早期第2阶段第3次岩浆侵入形成的花岗岩，岩性多为中细粒白云母花岗岩及少量黑云母花岗岩，成为钨、钼、铋、铍、钽、铌等有色金属和稀有金属矿床的成矿母岩.根据室内岩石力学参数试验测试与分析结果，可获得典型岩样物理参数及力学参数，拟取用5种岩体不同的抗拉强度参数进行研究，如表1所示.

表1 矿岩物理力学参数

2 数值模拟方案

2.1 模型建立

图1 模型计算域选取

根据计算精度与计算机计算能力的要求，选择数值计算模型[11]区域范围，如图1所示.网格划分采用四面体单元[12].为便于找出隐伏空区几何参数与露采平台厚度关系规律，计算模型略作简化.如图2所示，模型四面体单元数为472 635个，节点数80 551个.模型底面和四周均受链杆约束，顶面受荷载约束.

图2 整体模型切割后的计算模型

2.2 方案设计

一般情况下完整顶板其破坏方式有拉伸破坏和剪切破坏2种.由于岩体抗拉强度很低，一般会先发生拉伸破坏，然后拉裂隙扩展并与地表贯通而后发生冲剪破坏[13-15].实际空区形状很不规整，为了简化计算，这里将采空区视为矩形洞室进行考虑，设计计算方案如表2所示.

表2 计算方案参数表

由于钻孔设备、运输汽车、人员及辅助设备等外部荷载的存在，采空区中央上部平台载入合计0.5 MPa荷载，以摩尔-库仑（Mohr-Coulomb）强度理论准则[16]为数值模拟计算的本构模型，依据岩体的强度指标对顶板安全状况进行判别.观测和采集模型中破坏区关键点数据.

3 数值模拟结果与分析

3.1 空区顶板应力和位移分析

（1）主应力分布规律.限于篇幅原因，仅将方案6各步骤结果示列如下.如图3所示给出了方案6各开挖进度完成后的最大与最小主应力σ1，σ3分布情况.从模拟结果看出，主应力场具有以下特征：①空区形成后，与周围环境应力场比较，附近围岩应力释放明显.在空区顶板出现明显的拉应力集中区，且随着采空区高度的增加，顶板的拉应力逐渐增大，顶板拉应力大部分位于顶板中间附近.各开采步所形成空区都具有这种规律性；②空区顶板中央形成圆弧状的应力等值线拱，靠近顶板中央为拉应力，往上等值线拱径逐渐变大，拉应力减小，最终变为压应力.

图3 方案6各开采步完成后的最大与最小主应力分布情况

（2）剪应力分布规律.在空区周边产生剪应力集中，所不同的是最大应力集中区域不是分布在顶部或侧壁的中央，而是在空区的4个隅角处.在空区的右上角和左下角为正剪应力集中，左上角和右下角为负剪应力集中.

（3）围岩变形与位移特征.空区变形随空区空间的增大而不断加剧，甚至出现顶板冒落.分布特征为：空区周边围岩的位移最大，往外距离空区边界越远，围岩位移就越小，且顶部围岩位移比底部的要大得多，围岩移动方向均指向空区.

3.2 露采平台安全厚度与隐伏空区跨度关系分析

一般来说，隐伏采空区上覆岩体的破坏方式主要有2种：一种为在上部荷载作用下顶板两端围岩产生大变形导致上部覆岩发生剪切破坏，顶板整体滑落；一种为顶板中央产生拉应力超出岩体抗拉强度，出现冒顶、累积性破坏最后发展到地表形成大面积塌陷.根据模拟结果，采空区顶板破坏方式起决定性作用的是拉伸破坏．方案6各开采步围岩的拉伸不安全度等值云图如图4所示．

图4 方案6各开采步拉伸不安全度等值云图

拉伸不安全度是指岩体的拉应力与抗拉强度的比值乘以储备安全度，用Ft表示．当Ft＞1时，露采平台将发生破坏；当Ft＜1时，露采平台处于稳定状态，当Ft=1时，露采平台处于临界状态．各方案结果如表3所示．

利用FLAC3D数值模拟软件计算，地下空区储备安全度取值2．0时，分别对6种方案用线性插值法求出Ft=1时，对应露采平台的临界安全厚度值和空区跨度值，5组数据6种方案的线性回归函数关系曲线如图5所示．

表3 各计算方案结果

图5 露采平台临界安全厚度与隐伏空区跨度的关系

露采平台的临界安全厚度值与隐伏空区跨度及抗拉强度的线性回归函数如式（1）所示．

式（1）中：x为隐伏采空区跨度，m；y为露采平台临界安全安度，m；σc为岩体抗拉强度，MPa．

由此，可以得出露天平台安全厚度与隐伏采空区跨度及抗拉强度的关系，如式（2）所示．

式（2）中，D表示露采平台安全厚度，m．

4 结 论

（1）针对矿山地下转露天开采生产中存在隐伏采空区的安全隐患问题，应用FLAC3D数值计算软件，提出了露采平台安全厚度与空区跨度关系的线性拟合公式，为类似岩性和生产条件矿山的采矿工作提供借鉴意义.

（2）由于矿山现场存在地下水、雨水、自然风化等不确定因素，会使岩石的各项力学强度有所减弱，对研究结果造成一定的影响；同时论文仅做了静荷载分析，对于动态荷载作用（如爆破震动）未进行分析，需要在今后的工作中深入研究.

[1]Cooley T.Geological and geotechnical context of cover collapseand subsidence in mid-continent US clay-mantled karst[J].Environmental Geology，2002，42：469-475.

[2]Stump D.Minimizing the adverse impacts of surface coal mining on the environment[J].World Water&Environmental Resources Congress，2003（6）： 1-9.

[3]Swift G M，Reddish D J.Stability problems associated with an abandoned ironstone mine[J].Bulletin of Engineering Geology and the Environment，2002， 61（3）： 227-239.

[4]刘希灵，李夕兵，宫凤强，等．露天开采台阶面下伏空区安全隔离层厚度及声发射监测[J]．岩石力学与工程学报，2012，31（1）：3357-3362．

[5]Singh G S P，Singh U K，Murthy V M S R.Application of numerical modeling for strata control in mines[J].Geotech Geol Eng，2010， 28： 513-524．

[6]薛涛，李夕兵，周士霖，等．露天矿开采境界内下部复杂隐伏采空区上部稳定性研究[J]．矿业研究与开发，2012，32（2）：29-31．

[7]李地元，李夕兵，赵国彦．露天开采下地下采空区顶板安全厚度的确定[J]．露天采矿技术，2005（5）：17-20．

[8]张敏思，朱万成，侯召松，等．空区顶板安全厚度和临界跨度确定的数值模拟[J]．采矿与安全工程学报，2012，29（4）：543-548．

[9]Niehoff J.The use of geophysical methods to detect abandoned mine workings[C]//James W，Niehoff J.The progress of geological and geotechnical engineering in colorado at the cusp of a new decade.Lakewood：ASCE，2010：119-128.

[10]Zarnani S， Sriskandakumar S， Arenson L，et al.The use of FLAC software for assessing deformation rates during ice excavation for open pit mining in glaciers[J].Cold Regions Engineering，2012（8）： 705-714.

[11]廖秋林，曾钱帮．基于ANSYS平台复杂地质体FLAC3D模型的自动生成[J]．岩石力学与工程学报，2005，24（6）：1010-1013.

[12]GUANZhen-qun，SONGCao，GUYuan-xian.Theboundaryrecovery and sliver elimination algorithmsof three-dimensional constrained delaunaytriangulation[J].International Journal for Numerical Methodsin Engineering，2006，68（2）： 192-209.

[13]张耀平．矿山空区诱发的岩移特征及覆盖层冒落效应研究[D]．长沙：中南大学，2010．

[14]黄晖．缓倾斜硬岩矿床采空区顶板冒落规律研究[D]．赣州：江西理工大学，2011．

[15]张耀平，曹平，袁海平，等．复杂采空区稳定性数值模拟分析[J]．采矿与安全工程学报，2010，27（2）：233-238．

[16]蔡美峰，何满潮，刘东燕．岩石力学与工程[M]．北京：科学出版社，2002．