落凼矿分段空场嗣后充填法采场结构参数研究

1、 落凼矿分段空场嗣后充填法采场结构参数研究贾学元1，乔登攀1，何名声2，王俊1，李广涛1，杨天雨1（1.昆明理工大学 国土资源工程学院，昆明 650093；2.凉山矿业股份有限公司，凉山 615141）摘要： 分段空场嗣后充填法是两步骤开采，采场结构参数主要是指一步骤矿柱与二步骤矿房的结构尺寸。针对拉拉铜矿落凼矿区深部矿段地下开采遇到高地压释放致采场结构不稳定问题，采用RHINO-FLAC3D进行数值模拟的方法对矿区的采场结构参数进行研究。对落凼矿矿体厚度在30m以上的倾斜及急倾斜矿体采用两种不同的采场结构参数方案进行对比，用Flac3D模拟开挖，对采区内的塑性区进行统计与分析。研究结果表明：

2、随着采深的加大，针对该矿矿体厚度在30m以上的倾斜及急倾斜矿体，一步骤矿柱采用50m11m，二步骤矿房采用50m25m的结构参数更加安全稳定。关键词： 分段空场嗣后充填法；塑性区；采场结构参数；Flac3D数值模拟Study on Structural Parameters of Stope with Sublevel Open Stope and Subsequent Filling in Luodang MineJIA Xue-yuan1 QIAO Deng-pan1 HE Ming-sheng2 WANG Jun1 LI Guang-tao1 YANG Tian-yu1(1.School

3、 of Land Resources Engineering,Kunming University of Science and Technology,Kunming 650093,China; 2. Liangshan Mining Co Ltd, Liangshan 615141, China.) Abstract: Sublevel open stope and subsequent filling is a two-step mining method,stope structural parameters mainly refer to the structural dimensio

4、ns of one-step pillar and two-step room. Aiming at the instability of stope structure caused by high ground pressure release during underground mining in Luodang mining area of Lala copper mine, the numerical simulation method of RHINO- 基金项目：云南省基础研究专项-青年项目（编号：202101AU070022）；昆明理工大学人培基金（编号：NO.KKZ3202

5、021040）；昆明理工大学人培基金（编号：KKSY201921017）。1.贾学元 (1996-)男，采矿工程，硕士研究生，主要研究矿山充填技术，昆明理工大学 HYPERLINK TelTelEmail：952757663通讯作者：乔登攀(1969-)男，采矿工程，教授，博士生导师，主要研究地下开采工艺，矿山充填技术，昆明理工大学，Email： HYPERLINK mailto:2724783849 2724783849FLAC3d was used to study the stope structure parameters of th

6、e mining area. Two different schemes of stope structure parameters are used to compare the inclined and steeply inclined ore bodies with ore body thickness of more than 30m in Luodang mine. FLAC3D is used to simulate the excavation, and the statistics and analysis of the plastic zone in the mining a

7、rea are carried out. The results show that: with the increase of mining depth, for the inclined and steeply inclined ore body with the thickness of more than 30m, 50m 11m pillar is adopted in one step, 50m 25m for two-step room is more safe and stable.Key words: sublevel open stope and subsequent fi

8、lling；plastic zone；the parameters of stope structure；FLAC3D numerical simulation0引言采矿方法是矿山安全生产的核心，采矿方法的选择直接决定了矿山生产运营过程中的效益及安全。现阶段随着国家政策的支持及矿山自我认识的提高，我国地下矿山采用充填法的比例逐渐提升 REF _Ref72850418 r h * MERGEFORMAT 1。分段空场嗣后充填采矿法采用两步骤进行开采，先采一步骤矿柱，隔一采一，对一步骤矿柱进行充填，充填体达到一定强度，进行二步骤矿房的开采，对二步骤矿房进行充填，二步骤矿房是在一步骤充填后留下的充填

9、体矿柱支撑下进行开采，一步骤矿柱的充填质量是二步骤矿房安全开采的保障，该方法经济高效，同时又能保障井下安全作业。采场结构参数是影响采场稳定性的重要因素之一 REF _Ref72912258 r h * MERGEFORMAT 2。在采矿方案设计时，若设计的矿柱尺寸过大、矿房跨度过小，将导致矿石损失大、回收率低；若设计的矿柱尺寸过小，矿房稳定性降低、生产成本提高，且易发生大规模地压活动 REF _Ref72912198 r h * MERGEFORMAT 3。国内学者对采场结构参数优化进行了大量研究，初步范围的确定有Mathew法 REF _Ref75441792 r h * MERGEFORM

10、AT 4，单一的数值模拟手段的有FLAC3D REF _Ref75441852 r h * MERGEFORMAT 5，ANSYS REF _Ref75441162 r h * MERGEFORMAT 6，秩和比法 REF _Ref75441239 r h * MERGEFORMAT 7，采用两种数值模拟结合的有ANSYS-FLAC3D REF _Ref75441308 r h * MERGEFORMAT 8，SUBPAC-FLAC3D REF _Ref75441315 r h * MERGEFORMAT 9，考虑多种因素优化采场结构参数的有响应面法 REF _Ref75441430 r h

11、* MERGEFORMAT 10，CRITIC法 REF _Ref75441374 r h * MERGEFORMAT 11，AHP法 REF _Ref75441344 r h * MERGEFORMAT 12和TOPSIS法 REF _Ref75441344 r h * MERGEFORMAT 12。 针对拉拉铜矿落凼矿区深部矿段地下开采遇到的高低压导致采场结构稳定性变差问题，研究基于Flac3D数值模拟对采场结构参数进行优化。Flac3D在数值计算过程中功能十分强大，但对于复杂赋存形态、不规则地质体的前处理建模上有不足之处 REF _Ref73001036 r h * MERGEFORMA

12、T 8,研究用RHINO软件建模，导入FLAC3D数值模拟软件进行计算，克服了单一采用FLAC3D存在建模，网格划分困难等问题。1工程概况拉拉铜矿兴建于1958年，为凉山矿业股份有限公司下属子公司之一，矿山经过多次改扩建，目前实际生产为能力5000t/d，其落凼矿区具有多年露天开采历史，已形成了大、小两个的东、西露采坑。矿山近年来开展了落凼矿区深部的探矿工作，扩大了矿权范围，为延长矿山开采服务年限提供了资源条件。地下开采矿体形态较为复杂，矿体倾角多在1540之间，矿体厚度多在2m32m之间。现阶段，落凼矿处于露天转地下开采的关键时期，矿体厚度在78m以上的水平至倾斜矿体和矿体厚度在30m以上的

13、倾斜及急倾斜矿体采用分段空场嗣后充填法进行开采，研究选取矿体厚度在30m以上的倾斜及急倾斜矿体用Flac3D数值模拟软件进行分析，优化采场结构参数，为进一步开采提供科学的指导。2 建模方案考虑到倾角厚度的变化及建模工作量，本次模拟选用平均倾角45、厚度30m作为矿体产状，为了克服FLAC3D建模、网格划分困难，选用RHINO建模后导入FLAC3D进行计算，如图1所示。为优化分段空场嗣后充填法采场结构参数，本次数值模拟设计两个方案进行对比，如表1所示。表 1 数值模拟方案对比Table 1 Comparison of numerical simulation schemes方案一方案二一步骤矿柱

14、（m）50115011二步骤矿房（m）50255020方案一一步骤矿柱沿走向、倾向均有布置，矿房宽度11m，二步骤矿房沿走向布置，矿房宽度25m，先回采一步骤矿柱，胶结充填后回采二步骤矿房并充填。本次建模区域中包含二次暴露的一步骤矿柱12个，其中沿走向6个，沿倾向6个。建模区域如下图2所示，总单元数（zone）为901050，节点数（gridpoint）为898571。方案二与方案一的区别在于二步骤矿房宽度为20m。图 1 RHINO建模区域示意图Fig. 1 Schematic diagram of rhino modeling area图 2 由RHINO导入FLAC3D建模区域示意图Fi

15、g. 2 Schematic diagram of modeling area imported from rhino to FLAC3D 2.1 计算模型的选择本次数值模拟选用的本构模型为mohr-coulomb，各种力学性质岩体视为各向同性弹塑性连续介质，并在数值模拟计算过程中，应用摩尔-库伦塑性本构模型表示这种各向同性的塑性连续介质，其中摩尔-库伦屈服准则为：fs=1-31+sin1-sin-2C1+sin1-sin式（2-1） ft=t-3式（2-2）式中，1，3，t分别为最大主应力、最小主应力和抗拉强度，c，分别为粘结力和内摩擦角，当fs0时，材料将发生剪切破坏；ft0时，材料将发生

16、拉伸破坏。2.2 力学参数的取值本次数值模拟选取的力学参数如下表2所示。表 2 计算岩石力学参数选取结果表Table 2 Table of selection results of calculated rock mechanics parameters名称弹性模量(Gpa)泊松比密度(g/cm3)内摩擦角()抗拉强度(MPa)内聚力(MPa)围岩5.5430.382.857380.961.15矿体7.7990.242.818380.921.2充填体（胶结）0.2610.312.2400.130.32.3 边界条件落凼矿区矿体埋深较浅（300m），故竖直应力取自重应力，水平应力随深度按梯度变化

17、。v=H式（2-3）h=1-v式（2-4）式中为矿岩容重，=g；g取9.81m/s2，N/m3；H为埋深，m；为泊松比。模型的四个侧面和底面均采用固定速度约束。2.4 方案步骤模拟矿房矿柱布置如图3所示。方案模拟开采过程共分为48步，见下表3。图 3 模拟区域矿房矿柱布置示意图Fig. 3 Layout of room and pillar in simulated area表 3 分段空场嗣后三维模型模拟计算步骤Table 3 Simulation steps of three-dimensional model for segmented open field步骤模拟计算初始地应力计算模拟

18、三维模型原岩应力状态第一步（开挖1）开挖沿走向步骤胶结矿柱A1第二步（充填1）充填沿走向步骤胶结矿柱A1第三步（开挖2）开挖沿倾向步骤胶结矿柱B1第四步（充填2）充填沿倾向步骤胶结矿柱B1第五步（开挖3）开挖沿走向步骤胶结矿柱A7第六步（充填3）充填沿走向步骤胶结矿柱A7第七步（开挖4）开挖沿走向步骤矿房A4第八步（充填4）充填沿走向步骤矿房A4第四十六步（充填23）充填沿走向步骤矿房A17第四十七步（开挖24）开挖沿走向步骤矿房A18第四十八步（充填24）充填沿走向步骤矿房A18备注胶结矿柱胶结充填、矿房非胶结尾砂充填2.5数值模拟计算结果分析分段空场嗣后充填法按表3所示的回采顺序对矿体进行

19、模拟开挖、充填，步骤胶结矿柱和步骤矿房开挖结束后按表2中的强度设计、尾砂充填体参数进行赋值。分析胶结矿柱在区段回采过程中最危险状态即一侧临空另一侧受尾砂主动压力作用时的塑性区分布，并统计塑性区体积，从而判断胶结充填体的稳定性。2.5.1数值模拟结果以下为两个方案步骤矿柱的塑性区分布情况，由于篇幅限制，无法列出所有开挖方案对应的所有矿柱塑性区分布图，每个方案只列举沿走向和倾向各两幅图展示。所有开挖完成后，对沿走向步骤矿柱A8进行塑性区统计，开挖17沿走向步骤矿柱A8塑性区分布如图4所示，方案一塑性区总体积6848.8m3，矿柱总体积23499.7 m3，塑性区占比29.14%，方案二塑性区总体积

20、6580.6 m3，矿柱总体积21308.7 m3，塑性区占比30.88%，方案二比方案一塑性区占比大，稳定性较差。开挖17沿走向步骤矿柱A8塑性区分布如图4所示。方案一 （二步骤矿房50m25m）（b）方案二（二步骤矿房50m20m）图 4 开挖17沿走向步骤矿柱A8塑性区分布Fig.4 Distribution of plastic zone of I step pillar A8 along strike in excavation17所有开挖完成后，对沿走向步骤矿柱A14进行塑性区统计，方案一塑性区总体积5380.8m3，矿柱总体积23499.7 m3，塑性区占比22.90%，方案二塑

21、性区总体积5334.1 m3，矿柱总体积21308.7m3，塑性区占比25.03%，方案二比方案一塑性区占比大，稳定性较差。开挖23沿走向步骤矿柱A14塑性区分布如图5所示。方案一 （二步骤矿房50m25m） （b）方案二 （二步骤矿房50m20m）图 5 开挖23沿走向步骤矿柱A14塑性区分布Fig. 5 Distribution of plastic zone of I step pillar A14 along strike in excavation 23所有开挖完成后，对沿倾向步骤矿柱B1进行塑性区统计，方案一塑性区总体积2300.1m3，矿柱总体积15549.6m3，塑性区占比14

22、.79%，方案二塑性区总体积2386.3 m3，矿柱总体积12931.0 m3，塑性区占比18.45%，方案二比方案一塑性区占比大，稳定性较差。开挖8沿倾向步骤矿柱B1塑性区分布如图6所示。（a）方案一 （二步骤矿房50m25m） （b）方案二 （二步骤矿房50m20m）图 6 开挖8沿倾向步骤矿柱B1塑性区分布Fig. 6 Distribution of plastic zone of I step pillar B1 along dip in excavation 8所有开挖完成后，对沿倾向步骤矿柱B6进行塑性区统计，方案一塑性区总体积878.1m3，矿柱总体积15549.6 m3，塑性区

23、占比5.65%，方案二塑性区总体积1620.0 m3，矿柱总体积12931.0m3，塑性区占比12.53%，方案二比方案一塑性区占比大，稳定性较差。开挖24沿倾向步骤矿柱B6塑性区分布如图7所示。 方案一 （二步骤矿房50m25m）（b）方案二 （二步骤矿房50m20m）图 7 开挖24沿倾向步骤矿柱B6塑性区分布Fig. 7 Distribution of plastic zone of I step pillar B6 along dip in excavation 24部分沿走向矿柱塑性区统计结果如表4所示，为了方案对比便利，绘制成柱状图，如图8所示。在进行开挖时，两个方案A15矿柱塑性

24、区占比近乎相同，其它矿柱方案二比方案一塑性区占比要大，说明方案一比方案二更加稳定。表 4 沿走向矿柱塑性区统计结果Table 4 Statistical results of plastic zone along strike pillar A7A8A9A13A14A15塑性区总体积（m3）方案一5528.36848.85418.22719.25380.84762.9方案二5799.26580.65621.54878.45334.14329.2矿柱总体积（m3）方案一23499.723499.723499.723499.723499.723499.7方案二21308.721308.721308

25、.721308.721308.721308.7塑性区占比方案一23.52%29.14%23.06%11.57%22.90%20.27%方案二27.22%30.88%26.38%22.89%25.03%20.32%图8沿走向矿柱塑性区占比Fig. 8 Proportion of plastic zone along strike pillar部分沿倾向矿柱塑性区统计结果如表5所示，为了方案对比便利，绘制成柱状图，如图9所示。在进行开挖时，方案二矿柱整体要比方案一矿柱塑性区占比大，说明方案一比方案二更加稳定。表 5沿倾向矿柱塑性区统计结果Table 5 Statistical results of

26、 plastic zone along dip pillar B1B2B3B4B5B6塑性区总体积（m3）方案一2300.11854.4911.12774.52159.8878.1方案二2386.32070.51508.12677.72269.11620.0 矿柱总体积（m3）方案一15549.615549.615549.615549.615549.615549.6方案二12931.0 12931.0 12931.0 12931.0 12931.0 12931.0 塑性区占比方案一14.79%11.93%5.86%17.84%13.89%5.65%方案二18.45%16.01%11.66%20

27、.71%17.55%12.53%图9沿倾向矿柱塑性区占比Fig. 9 Proportion of plastic zone along dip pillar2.6方案综合分析矿体厚度在30m以上的倾斜及急倾斜矿体选用分段空场嗣后充填法数值模拟，该类型矿体采用方案一和方案二数值模拟计算。对比两个方案的塑性区占比，对两个方案的稳定性情况进行综合分析，两个方案的数值模拟如上所示。由数值模拟分析结果来看，沿走向矿柱比沿倾向矿柱塑性区占比更大，破坏更严重，塑性区集中在采场顶部和临空面，矿柱整体性较好。随开采区域的增大，顶板暴露面积增大，作用于胶结矿柱的松脱地压增大，超过胶结矿柱顶部设计强度时，顶部发生破

28、坏。在模拟区域回采过程中，胶结矿柱受四周矿房回采的影响多次暴露，围岩应力重新分布过程中，胶结矿柱应力状态随之变化，由矿岩包裹的三维力学状态转变为单侧暴露的三维力学状态，胶结矿柱变形向空区内发展，临空侧单元体应力状态沿加载路径向强度曲线靠近，部分单元体应力状态达到极限应力状态发生破坏。由方案一和方案二对比可知，矿柱在50m11m参数不变条件下，在整体上方案二（矿房50m20m）沿走向和倾向上矿柱塑性区占比明显比方案一（矿房50m25m）要大，稳定性更差，故选择方案一。结论（1）随着采深的加大，拉拉铜矿落凼矿深部开采按原矿柱、矿房的采场结构参数设计势必导致地压显现频繁,需要优化采场结构参数指导进一

29、步的开采。（2）落凼矿深部采用分段空场嗣后充填法开采，塑性区并未贯穿胶结矿柱，胶结矿柱整体性较好，采用方案一（一步骤矿柱50m11m，二步骤矿房50m25m）和方案二（一步骤矿柱50m11m，二步骤矿房50m20m）的采场结构参数均不会出现失稳破坏。胶结矿柱塑性区主要分布于胶结矿柱顶部和临空面，而顶部出现塑性区的主要原因是充填体在分层充填过程中顶部选用的充填配比强度较低导致。（3）通过FLAC3D数值模拟软件对采区进行参数优化，综合安全稳定性方面考虑，最终选择了一步骤矿柱50m11m，二步骤矿房50m25m为最佳采场结构参数，为拉拉铜矿落凼深部矿区的进一步开采提供了科学的指导，也对采用分段空场

30、嗣后充填法的相似矿山有一定的借鉴意义。参考文献任卫东，范玉乾.分段空场嗣后充填采矿法集约化开采顺序优化研究J.金属矿山.2020，55(6):167-171.REN Weidong，Fan Yuqian. Study on Optimization of Mining Sequence in Intensive Mining with Sublevel Open Stope and Subsequent FillingJ.Metal Mine. 2020，55(6):167-171.周科平，苟国刚，罗佳，等.基于 Mathew与数值模拟法的采场结构参数优化及稳定性研究J.科技导报，2013，3

31、1（30）:34-38.ZHOU Keping，GOU Guogang，LUO Jia，et al. Optimization of Structural Parameters and Stability of Stope Based on Mathew Method and Numerical Simulation MethodJ. Science and Technology Review，2013，31(30)：3438.甯瑜琳.深部缓倾斜薄矿体采场结构参数优化及采动沉陷特性研究D.长沙:中南大学，2011.NING Yulin. Research on Optimization of

32、Stope Parameters of Deep Gently Inclined Thin Orebody and Mining Subsidence Characteristics D.Changsha: Central South University, 2011.刘泽州，吴伟伟，雷鸣.缓倾斜厚矿体采场结构参数优化研究J.矿业研究与开发，2020，40(9):6-10.LIU Zezhou,WU Weiwei,LEI Ming.Study on Optimization of Stope Structure Parameters and Support System for Gently-

33、inclined Thick Ore BodyJ. Mining Research and Development, 2020，40(9):6-10.王炳文，黎林，张书金，等.基于FLAC3D模型的上向进路采场结构参数优化J.中国矿业，2020,29(9):93-96.WANG Bingwen,LI Lin，ZHANG Shujin，et al. Optimization of Structural Parameters of Upward Drift Stope Based on FLAC3D ModelJ. China Mining Magazine, 2020,29(9):93-96.罗

34、朋，王振昌，张元庆等.紫金山深部铜矿采场结构参数研究J.矿业工程，2017(4)：21-25.LUO Peng,WANG Zhenchang,ZHANG Yuanqing,et al.Study on Stope Structural Parameters of Zijinshan Deep Copper MineJ.Mining Engineering, 2017(4)：21-25.周科平，王星星，高峰.秩和比法耦合数值模拟优化采场结构参数研究J.矿冶工程2013，33(1):1-4.ZHOU Keping，WANG Xingxing，GAO Feng. Optimization of Stope Structural Parameters Using Coupled Numerical Simulation with Rank-sum Ratio MethodJ. Mining and Metallurgical Engineering，2013，33(1): 1-4.马春德，徐家庆，刘焕新等.基于ANSYS-FLAC3