深部松散矿体开采技术数值模拟研究

徐培良，王建国，石 飞，周山山，杨德源，何 进

（1.云南锡业股份有限公司大屯锡矿，云南 个旧 661000；2.云南农业大学建筑工程学院，云南 昆明 650201；3.昆明勘晟矿业有限公司，云南 昆明 650100；4.昆明钢铁集团公司大红山铁矿，云南 新平 653100）

深部松散矿体开采技术数值模拟研究

徐培良1，王建国2，石 飞3，周山山1，杨德源4，何 进2

（1.云南锡业股份有限公司大屯锡矿，云南 个旧 661000；2.云南农业大学建筑工程学院，云南 昆明 650201；3.昆明勘晟矿业有限公司，云南 昆明 650100；4.昆明钢铁集团公司大红山铁矿，云南 新平 653100）

云锡松矿1360中段以下30-14矿体松散破碎，开采难度较大。为满足深部安全及经济合理的开采要求，通过对房柱法和下向进路式分层胶结充填法进行技术经济比较，以及两种方法的数值模拟分析，确定采用下向进路式分层胶结充填法作为30-14矿体的采矿方法。该方法虽然成本高于房柱法，但回采的矿石资源多，创造的经济价值大，比房柱法多获得利润总额5 408万元，且开采的安全性更能得到保障。研究成果为类似矿体开采积累了经验，对矿山生产建设和深部资源开发利用都有实际的指导意义。

深部开采；松散矿体；下向进路式分层胶结充填法；经济指标；安全率

金属矿山深部开采面临一系列新的技术难题，如采准巷道维护更加困难、底板突水事故增加、冲击地压频发等[1]。而深部破碎矿体形态比较复杂，矿体内穿插有围岩夹层，顶底板的起伏变化大，易发生局部或区域性的地压灾害，巷道和采场垮塌、冒落严重，掘进和支护十分困难，所以迫切需要对深部复杂条件下的矿体开采进行研究[2]。杨清林等[3]根据某金矿极破碎矿床的地质条件，对采场参数进行数值模拟优化，确定了可行的回采方案。张宝等[4-5]对矿体在大体积充填体下的回采顺序和胶结充填采矿法盘区回采工艺进行数值模拟研究，确定了最优回采顺序和充填工艺。黄文等[6-8]针对地下钨矿山出现的巷道开裂、垮塌、冒落及采空区大面积塌方和地表沉陷等问题，对采空区进行了废石或尾砂胶结充填处理，不仅有效控制了矿区开采过程中地压活动对井下工程的破坏，还降低了生产成本。但本研究矿体开采条件极为复杂，需对开采方案进行设计比较和数值模拟分析[9]。

1 采矿方法的初步确定

1.1 工程概况

云南锡业集团云锡松矿1360中段以下30-14矿体产出于花岗岩与大理岩接触带，控制走向长约200 m，宽120 m，铅垂厚1～18 m，矿体赋于1 360 m中段潜水面以下，上部有1360m平台主运输巷道和排水巷道穿过，如图1所示。矿体围岩顶板为大理岩，节理裂隙发育，底板为夕卡岩、长英岩、花岗岩。花岗岩的顶部呈云英岩化、绢云母化、绿泥石化，矿石呈强风化～半风化，稳固性极差，黄铁矿、磁黄铁矿、黄铜矿等硫化物呈稀疏－稠密状、细脉状浸染分布，矿体多呈透镜状、似层状，节理裂隙发育，稳固性差。

图1 矿体与巷道投射关系Fig.1 Projection relation between orebody and roadway

1.2 采矿方法初步分析

根据矿体开采技术条件，初步可选采矿方法有：房柱法；下向进路式分层胶结充填法；上向水平分层充填法。采用上向水平分层充填法开采时，回采作业面顶板安全性得不到保障，施工安全投入较大，需要从矿体最低端开始回采，矿体的生产准备时间较长，难以满足当期的生产需求，可以排除。房柱法和下向进路式分层胶结充填法主要技术经济指标对比如表1所示。

然后再对下向进路式分层胶结充填法和房柱法进行数值模拟计算，以确定最优的开采技术方案。

表1 房柱法和下向进路式分层充填法主要技术经济指标对比表Tab.1 Comparison table of main technical economic indicators of two schemes

2 数值模拟准备

2.1 岩石室内试验

针对30-14矿体顶板（大理岩）、氧化矿、底板（花岗岩）结构及强度的特点，选择一定数量的岩块进行加工并进行强度测定，得到数值计算的基本参数如表2所示。

表2 数值模拟基本参数Tab.2 Basic parameters of numerical simulation

2.2 数值计算指标

安全率是由摩尔-库仑强度准则所决定的极限应力状态与实际应力状态的比值。安全率为1时处于临界状态，且安全率越大，安全性越好（对地下矿山，一般认为安全率大于1.15是偏于安全的）。对于地下矿山而言，顶板、间柱、围岩等是否破坏，除了分析应力应变、塑性区分布及位移外，安全率的分布，在某种程度上，可以非常直观明了地对开挖所造成的损害程度进行说明。

3 充填效果数值模拟分析

采用FLAC3D进行数值模拟，模拟进路规格为4m×4m，上下分层的进路方向垂直交错布置。由于胶结充填进路式采矿开采[9]及充填步骤较多，数据量较大，在此对中间回采、充填等不再进行分析，仅对整个分层充填后充填体的作用及顶板稳定性进行分析。

3.1 应力分布模拟结果分析

充填体及充填完成后顶板应力分布模拟结果见图2、图3。图2模拟结果显示，充填体所承受的最大主应力为48.548 MPa；图3模拟结果显示，整个分层充填结束后，顶板最大主应力达到45.87 MPa。充分说明充填体的强支护作用；矿体顶板没有拉应力出现，此时顶板是趋于稳定的。

图2 充填体最大主应力分布Fig.2 The maximum principal stress distribution of the filling body

图3 充填后顶板最大主应力分布Fig.3 Main stress distribution of the back top plate after filling

3.2 安全率分布模拟结果分析

从充填体及充填完成后顶板安全率分布模拟结果（图4、图5）看，整个分层充填结束后，充填体的最小安全率为1.197，充填后顶板最小安全率为1.237，均大于临界状态，满足地下矿山安全的经验值，这也充分说明了充填体的存在大大提高了顶板的稳定性。

图4 充填体安全率分布Fig.4 Distribution of backfill safety rate

图5 充填后顶板安全率分布Fig.5 Safety rate distribution of top plate

3.3 塑性区分布模拟结果分析

从充填体及充填完成后顶板塑性区分布模拟结果（图6）看，整个分层充填结束后，充填体大部分进入了塑性状态，即进入承载状态，说明充填体虽然不能完全阻止围岩移动，但可以最大限度地限制围岩的变形和移动。充填完成后顶板几乎没有塑性区出现，仅在局部有零星的塑性区出现，此时顶板是稳定的，说明充填体在抑制顶板变形中起到了至关重要的作用。

图6 充填体塑性区分布Fig.6 Distribution of pack plastic zone

3.4 房柱采矿法数值模拟分析

宽度10 m，强度C20的混凝土房柱采矿法的数值模拟计算结果显示：开采后顶板拉应力最大值为2.16 MPa，大于大理岩折减后的抗拉强度值（约为1.8 MPa）；顶板位移指向采空区；最小安全率为1.12，不满足地下矿山安全率要求；顶板大理岩塑性区比较多，且呈连片现象。以上结果均说明此时顶板暴露面积超出了大理岩自身的极限暴露面积，会造成顶板变形、垮塌。相比之下，采用进路式分层胶结充填采矿能最大限度地抑制顶板和围岩的变形，满足矿体开采的安全要求。

4 下向进路式分层胶结充填采矿法

4.1 采场结构参数确定

矿块沿走向和垂直走向布置，长15～74 m，自上而下分为八层，各分层宽度如下：第一层46 m、第二层80m、第三层56m、第四层62m、第五层56m、第六层82 m、第七层35 m、第八层46 m，每层高4 m，进路宽4 m，不留间柱和顶底柱。采场搬运方式为铲运机及矿用汽车。在各分层内设置分层联道和回采进路，分层联道用于分层内运输、充填和通风，回采进路用于采矿。

4.2 矿体回采顺序

将矿体分为南北两个盘区，两个盘区可同时回采，互不影响。南盘区从第一层往下回采，一共回采五个层；北盘区从第五层的西坡开始往下回采，一共回采四个层。南北盘区回采时，严格遵守由上而下的顺序进行开采，即上一分层采充结束之后方可进行下一分层的开采。

4.3 采切工程布置

在矿体主斜坡道上按各分层的高度，施工分层斜坡道达到矿体，沿矿体边界掘进，探清矿体形态，首先施工进路与回风井贯通，每条进路回采的进路规格为4 m×4 m，即每一分层的高度为4 m，为了保证充填体稳定性，相邻的上下两分层回采进路呈近似垂直交错布置（第一、三、五、七层进路为25°～205°方向，第二、四、六、八层为 295°～115°方向）。

4.4 回采工艺

采用自上而下回采，进路高度为4.0 m，宽度为4.0 m。采用阿特拉斯S1D掘进凿岩台车进行凿岩爆破，采用大直径空眼桶形掏槽爆破，空眼直径100mm，孔深3.4 m，其他装药炮孔直径45 mm，孔深3.2 m。偏进路底板布置大孔径掏槽，非电导爆管微差爆破，这样的爆破设计可减小爆破震动对顶板充填体的破坏作用，确保充填体的安全稳定；同时，在顶板充填体下部的碎石垫层能有效降低爆破冲击波对充填体的影响。

4.5 采场充填

采场清理完毕后，从里向外按排距1.2 m，间距1.2 m打吊挂锚杆，锚杆内焊接吊筋，吊筋的直径为16 mm，巷道底板钢筋以焊接方式与吊筋相连，由主筋和副筋构成，主筋直径16 mm，网格为1.2 m×1.2 m，副筋直径10～12 mm，网格为0.3 m×0.3 m。沿邻侧回采进路的巷道帮预留0.5 m长的主筋，并用直径不小于为12 mm的钢筋将留出的主筋连成整体，回采结束铺设底筋时，与预留主筋相焊接，使两个相邻的进路成为一个整体，以保证充填体的稳定性。

5 结语

通过对房柱法和下向进路式分层胶结充填法的技术经济比较和数值模拟分析可以看出：

（1）房柱法在生产能力上要略占优势，但回采可能会诱发运输、排水坑道的变形，导致透水事故发生，不适宜在开采中使用。

（2）下向进路式分层胶结充填法虽然成本要高，但是开采中的安全性能得到充分保障，回采的矿石资源多，创造的经济价值大（比房柱法多获得利润总额约5 408万元），符合“安全、高效、低成本、低贫损”的资源开发理念。因此，本次设计选择下向进路式分层胶结充填法作为30-14矿体的采矿方法。

[1] 徐培良.破碎、含水复杂条件矿体开采技术研究[D].昆明:昆明理工大学，2015.XU Peiliang.Research on mining technology of mining body with brokenandcomplicatedconditions[D].Kunming：KunmingUniversity of Science and Technology，2015.

[2] 张卫中，胡泽图，倪小山，等.复杂条件下厚大磷矿体安全高效采矿方法研究[J].矿业研究与开发，2017，37（7）：1-4.ZHANG Weizhong，HU Zetu，NI Xiaoshan，et al.Study on the safe and efficient mining method for thick and large phosphate rock in complex conditions [J].Mining Research and Development，2017，37（7）：1-4.

[3] 杨清林，王 星，孙国权，等.极破碎矿床开采的数值模拟与实践[J].现代矿业，2014（5）：1-3，27.YANG Qinglin，WANG Xing，SUN Guoquan，et al.Application and numerical simulation of extremely fractured ore deposit mining[J].Modern Mining，2014（5）：1-3，27.

[4] 张 宝，刘 涛，李向东，等.上向水平分层充填法试验盘区回采顺序数值模拟研究[J].矿业研究与开发，2016，36（10）：11-15.ZHANG Bao，LIU Tao，LI Xiangdong，et al.Numerical simulation on tested panel mining sequenceby upward horizontal slice stoping and backfilling[J].Mining Research and Development，2016，36（10）：11-15.

[5] 焦文宇，徐必根，唐绍辉，等.充填体下矿体回采顺序数值模拟优化研究[J].矿业研究与开发，2012，32（5）：9-11.JIAO Wenyu，XU Bigen，TANG Shaohui，et al.Simulation and optimization on stoping sequenceof orebody under backfill[J].Mining Research and Development，2012，32（5）：9-11.

[6] 黄 文.废石充填采空区在章源公司矿山的应用[J].中国钨业，2008，23（6）：5-7.HUANG Wen.Application of mullock stowing in a tungsten mine goaf[J].China Tungsten Industry，2008，23（6）：5-7.

[7] 万林海.全尾砂胶结充填在某白钨矿山的应用 [J].中国钨业，2014，29（2）：1-4.WAN Linhai.Application of whole tailings backfill in a scheelite mine[J].China Tungsten Industry，2014，29（2）：1-4.

[8] 曾凡珍，石求志，郑建新.某黑钨矿采空区废石充填研究[J].中国钨业，2015，30（4）：32-36，52.ZENG Fanzhen，SHI Qiuzhi，ZHENG Jianxin.Waste rock filling for themined-outareainawolframitemine[J].ChinaTungstenIndustry，2015，30（4）：32-36，52.

[9] 任少峰，杨 伟，张卅卅，等.不同采矿方法下采场稳定性数值分析[J].武汉理工大学学报，2015，37（3）：94-97.REN Shaofeng，YANG Wei，ZHANG Sasa，et al.Numerical analysis of stability of stope with different mining methods[J].Journal of Wuhan University of Technology，2015，37（3）：94-97.

[10] 李 力，李占金，杨立刚，等.不稳定矿体高效进路充填采矿方法优化研究与应用[J].矿业研究与开发，2014，34（4）：1-3.LI Li，LI Zhanjin，YANG Ligang，et al.Research and application on optimization of high-efficiency drift backfill mining method in unstable orebody[J].Mining Research and Development，2014，34（4）：1-3.

Numerical Simulation on the Mining Technology for Deep Loose Ore Body

XU Peiliang1,WANG Jianguo2,SHI Fei3,ZHOU Shanshan1,YANG Deyuan4,HE Jin2
(1.Datun tin ore,Yunnan Tin Co.,Ltd.,Gejiu 661000,Yunnan,China;2.College of Civil and Architectural Engineering,Yunnan Agricultural University,Kunming 650201,Yunnan,China;3.Kunming Kansheng Mining Co.,Ltd.,Kunming 650100,Yunnan,China;4.Dahongshan Iron Mine Bureau,Kunming Ironamp;Steel Group Corporation,Xinping 653100,Yunnan,China)

The 30-14 ore body under 1 360 middle section of Songshujiao mine is characteristics of loose,broken and large water content,which makes it difficult to be exploited.The method of downward direction layered cementing and filling mining technology is put forward to meet the requirements of deep safety and economic reasonable mining,after analyzing the stability of the mining field through on-site exploration and technical economic comparison.The new technology is more profitable and secure than the room pillar method.The results of this study have accumulated experience for similar ore mining,and have practical guiding significance for mine production and construction and deepresource developmentandutilization.

deep mining;loose ore body;method of downward direction layered cementing and filling mining technology;economic indicators;safety factor

TD862

A

10.3969/j.issn.1009－0622.2017.05.007

2017-08-18

云南省科技厅基础研究青年项目（2016FD029）；云南省教育厅科学研究基金项目（2016ZZX108）

徐培良（1985-），男，云南双柏人，硕士，工程师，主要从事地下工程和采矿技术方面的研究工作。

王建国（1987-），男，河南信阳人，博士，副教授，主要从事工程爆破和岩石动力学方面的研究工作。

（编辑：游航英）