缓倾斜厚大矿体采场结构参数优化研究与应用.pdf

Series No.569November 2023 金 属 矿 山METAL MINE 总 第569 期2023 年第 11 期收稿日期 2022-08-01基金项目 河南省高等学校重点科研项目计划(编号:21A560009)。作者简介 柴 红(1977),女,副教授,硕士。缓倾斜厚大矿体采场结构参数优化研究与应用柴 红1 郭俊超1 杨淑慧2 刘成博3(1.黄河水利职业技术学院土木与交通工程学院,河南 开封 475004;2.郑州大学水利与土木工程学院,河南 郑州 450001;3.山西省第三地质工程勘察院有限公司,山西 晋中 030620)摘 要缓倾斜厚大矿体采场结构参数的确定和优化对采场安全、经济、可靠回采等工作意义重大。以山东半岛某地下矿井为例,基于 Mathew 稳定图表法对采场地质参数和开采参数进行计算,确定了采场尺寸范围。利用FLAC3D数值模拟软件,对采场参数取值进行了进一步分析,综合经济、技术因素确定了采场参数合理取值。研究表明:在综合考虑利用围岩自身支撑能力的基础上,给出了合理的采场宽度范围为 28.3642.17 m。利用 FLAC3D数值模拟软件对采场进行数值模拟,通过对不同宽度方案条件下应力场与位移场的分析,明确了该采场 34 m 宽度为最优方案。关键词 缓倾斜厚大矿体 结构参数 Mathew 稳定图表 FLAC3D 矿柱 中图分类号TD853 文献标志码A 文章编号1001-1250(2023)-11-142-06DOI 10.19614/ki.jsks.202311016Study and Application on the Optimization of Stope Structural Parameters of Gently Inclined Thick OrebodyCHAI Hong1 GUO Junchao1 YANG Shuhui2 LIU Chengbo3(1.Department of Civil Engineering and Transportation Engineering,Yellow River Conservancy Technical Institute,Kaifeng 475004,China;2.Department of Water Conservancy and Civil Engineering,Zhengzhou University,Zhengzhou 450001,China;3.Shanxi The Third Institute of Geological Investigation Co.,Ltd.,Jinzhong 030620,China)Abstract The determination and optimization of stope structural parameters of gently inclined thick large orebody is of great significance to stope safety,economy and reliable mining.Taking an underground mine in Shandong Peninsula as the ex-ample,this paper calculates the stope geological parameters and mining parameters based on Mathew stability chart method,and determines the range of stope size parameters.Using FLAC3D numerical simulation software,the stope parameters are further an-alyzed,and the reasonable values of stope parameters are determined by integrating economic and technical factors.The study results show taht:Based on the comprehensive consideration of the support capacity of surrounding rock,a reasonable stope width range of 28.3642.17 m is given.Using FLAC3D numerical simulation software to simulate the stope,through the a-nalysis of stress field and displacement field under the conditions of different width schemes,it is indicated that the 34 m width of the stope is the optimal scheme.Keywords gently inclined thick orebody,structural parameters,Mathew stability chart,FLAC3D,ore pillar 在矿山实际生产过程中,采场结构布置参数优化必须从经济效益和安全条件两个方面考虑。良好的结构布置方案不仅能够提升矿石的采出率、提高采出矿石的品位,而且能够有效地保证矿柱、矿房内施工安全。矿山的安全与经济效益是矿山设计的核心要义,因此,矿山设计过程中,根据实际情况确定与优化采场结构参数就显得尤为重要。近年来,国内外诸多学者对采场结构参数优化进行了深入研究1-4。目前,常用的确定采场结构参数的方法主要有经验类比法5、数值模拟分析法6-9、理论分析法9-10等。刘泽洲等11基于查表法确定了厚大矿体的采场宽度,并利用经验方法与数值模拟方法对采场参数进行了优化。焦国芮等12、罗来和等13设计了采场结构参数方案,应用 FLAC3D软件进行数值模拟,探究了采场结构参数合理取值。周宝坤等14为了提高矿石品位,利用响应面法对矿房参数241进行了优化。陈顺满等15为了取获最优的采场结构参数,采用中心复合试验研究了采场高度、矿柱直径、矿柱排距和矿柱间距对破碎围岩条件下的矿柱最大拉应力、最大压应力和矿石回收率的影响规律。徐帅等16运用 SOM 模型,设计了多水平参数试验,并通过数值模拟软件,实现了扩展方案的正向预测、结构参数反向预测以及影响采场稳定性的主因素分析。张亚伟等17根据矿山实际情况,利用三维有限元动态仿真模拟方法对采场结构及推进方向进行了优化分析。兰明等18基于弹性力学理论,建立了板状结构力学模型,分析了厚板结构发生变化时所对应的拉应力变化特征,确定了阶段嗣后充填采场最优开采参数,并通过多因素分析法对理论计算结果进行了验证。本研 究 以 山 东 半 岛 地 区 某 地 下 矿 为 例,将Mathew 稳定图表法与数值模拟计算方法相结合,首先利用经验分析法综合考量采场的安全性与经济性,对采场尺寸参数进行预估;然后基于 FLAC3D数值模拟软件,以矿柱垂直应力、顶板拉应力与顶板位移为研究对象,对采场参数取值进行进一步分析,在此基础上,综合经济、技术因素确定最优采场结构参数。1 工程概况山东半岛某地下矿所开采的某矿段中富含了大量的金、钨、铜等多种金属矿体,开采段长约 143 m,倾向长度为 128186 m,厚度为 18.736.8 m,平均24.8 m。矿体上盘为花岗岩,普氏系数 f 为 1318;下盘为大理岩,普氏系数 f 为 1014;矿体较为致密且相对坚硬,普氏系数 f 为 9。矿块沿矿体走向方向布置,由于矿块在中间位置存在一定的夹矸,降低了矿石品位,因此在矿段中间位置、北偏东 12方向设置了宽 19 m 的永久性保护矿柱(图 1)。开采顺序分为两个阶段进行,首先是回采矿柱,待回采完毕后进行充填;其次是矿房回采,回采后再次进行充填。根据设计,矿柱宽度为 16 m,左右两个盘区的走向长度为 62 m,矿房宽度(即矿房的倾斜长度)需要进一步确定。采场矿房宽度越大将会导致整个采场的悬顶面积增大,矿柱的稳定性与顶部围岩的稳定性将同时减弱,因此本研究通过进一步分析给出合理的矿房宽度取值。采矿方案如图 1所示。2 基于 Mathew 稳定图表法的矿房宽度确定Mathew 稳定图表法是由英国 Golder 公司工程人员 MATHEW 提出的,在判断采场稳定性方面,是目前应用较为广泛的方法之一。其基本思想为综合考虑采场的岩石力学性质和开采技术因素经验性确定图 1 矿块开采尺寸Fig.1 Mining size of ore blocks采场的稳定性情况。随后 POTVIN 基于工程实例对Mathew 稳定性图表进行了发展,得到新的 Mathew 稳定性图表,如图 2 所示。图 2 中横坐标为采场暴露面形状系数 S,纵坐标为岩体稳定性指数 N,通过现场确定的采场尺寸和稳定性系数,从而可以确定坐标点位置,据此判断采场的稳定性。其中,稳定性区域分为 3 类,即稳定区、不稳定区和崩落区。因此,在确定采场区域时,根据采场稳定性系数 N,可以反算出稳定情况下的采场尺寸。图 2 基于 Mathew 稳定图表法的采场稳定性区域划分Fig.2 Regional division of stope stability based on Mathew stability chart method2.1 基于 Mathew 方法确定矿房宽度2.1.1 稳定性指数 N 确定稳定性指数 N 可进行如下计算:N=QABC,(1)式中,Q 为岩体质量指数,综合考虑了岩石的力学特性、所处的应力场以及岩体的空间关系;A 为岩石应力系数,与采场所受的应力有关;B 为岩体缺陷方位修正系数,综合考虑了岩体中的缺陷,例如弱面、节理、孔洞等影响;C 为岩石暴露的修正系数,与采场暴露情况有关。341 柴 红等:缓倾斜厚大矿体采场结构参数优化研究与应用 2023 年第 11 期岩体质量指数 Q 计算公式为Q=RJrJwJnJaSf,(2)式中,R 为岩体质量指标,即 RQD(Rock Quality Des-ignation),表示岩体良好度的指标;Jr为节理粗糙度系数,表征节理面形貌;Jw为节理裂隙水折减系数,表征节理裂隙水对节理的影响;Ja为节理蚀变、充填及胶结程度系数,表征节理初始状态的胶结程度;Jn为节理组数;Sf为应力折减系数。岩石应力系数 A 综合表征了两个方面因素:外在因素,由于矿石开采导致的应力场重新分布形成应力集中;内在因素,岩石自身的单轴抗压强度对应力系数 A 的影响。岩体缺陷方位修正系数 B 综合反映了顶板方位对稳定性的影响,其表达式为B=1.57 10-53-0.001 42+0.042 6-0.105,(3)式中,为顶板倾角,()。对于岩石暴露的修正系数 C,在水平情况下,C=1;其他情况下,取值为C=8-6cos,(4)式中,为暴露面与水平方向的夹角,()。2.1.2 采场形状系数 S形状系数 S 是主要表征暴露面形貌特征的参数,一般情况下采场为矩形,主要参数为矩形的长和宽。S 值可进行如下计算:S=L1L22 L1+L2(),(5)式中,L1为采场长度,m;L2为采场宽度,m。2.2 采场尺寸确定根据上述关于稳定性指数 N 和采场形状系数 S计算过程的描述,此时根据采场的空间力学特性能够计算得到稳定性指数 N,在此基础上,判断出稳定区、不稳定区以及崩落区的范围。由于设计过程中采场长度固定,因此通过 S 值能够确定采场宽度。基于采场的围岩力学性质与应力特征,根据式(1)至式(4)得到稳定性指数 N 与计算参数取值,见表 1。表 1 计算参数取值Table 1 Values of calculation parametersQABCN20.110.52.83828.38 根据确定的稳定性系数 N 取值,并根据图 2 可知,若采场处于稳定区,此时形状系数 S 为 09.73;若采场处于不稳定区,S 为 9.7312.55;若采场处于崩落区,S 为 12.5525。为防止采场崩落,在局部支护的情况下,S 为9.7312.55,既保证了安全,同时确保达到较高的产量,因此由式(6)计算得到的采场宽度 L2为 28.3642.17 m。L2=2SL1L1-2S(),(6)式中,根据采场设计参数,L1取 62 m。3 数值模拟计算3.1 模拟方案基于 Mathew 稳定图表法确定的采场宽度为28.3642.17 m,在保证矿柱宽度为 19 m、采场长度为 62 m 的前提下,进一步对采场宽度参数取值进行优化,根据采场宽度范围,设计的数值模拟方案见表 2。表 2 数值模拟方案Table 2 Numerical simulation schemes方案采场长度/m采场宽度/m悬顶面积/m2形状系数 SA62281 7369.64B62311 99210.71C62342 10810.98D62372 29411.59E62402 48012.163.2 FLAC3D数值模型构建FLAC3D软件是美国 ITASCA 公司开发的三维有限差分程序,能够进行土质、岩石和其他材料的受力分析和塑性流动分析。该款软件采用显式拉格朗日算法和混合离散分区技术,在塑性破坏和塑性流动方面有着显著优势19-20。同时由于其无须整合刚度矩阵,在解算三维问题时所需内存较小。故本研究采用该款软件进行数值模拟分析。根据采场尺寸,模拟所用到的网格模型通常为 35 倍的研究范围,因此构建的网格模型尺寸为 600 m600 m500 m(长宽高),如图 3 所示。模型边界条件为:底部约束垂直位移,前后左右 4 个边界采用固定界面法向位移,由于模型尺寸延伸至地表,因此不再对顶部施加载荷,采用自然边界。岩体采用Mohr-Coulomb 模型,模型参数取值见表 3。3.3 计算结果分析3.3.1 应力分析由于矿房中间位置为最容易发生崩落的位置,故本研究取中间位置的应力云图进行分析,结果如图 4所示。对矿柱上的垂直应力分别进行了统计,结果如图 5 所示。由图 4、图 5 可知:随着采场宽度增加,矿柱上垂直方向的应力增大;最大值出现在方案 E 中,矿柱中垂直应力为27.4 MPa,显然最大垂直应力小441总第 569 期 金 属 矿 山 2023 年第 11 期图 3 数值模拟网格模型Fig.3 Grid model for numerical simulation于矿柱岩石的抗压强度(31 MPa)。矿房上部的拉应力是预测岩体崩落的重要指标,一旦拉应力超过抗拉强度,往往将会从矿房中部位置产生岩块的崩落,造成安全隐患。对不同方案中的最大拉应力进行了统计,结果如图 6 所示。由图 6 可知:随着矿房宽度增加,拉应力逐渐增大,5 个方案对应的拉应力分别为 0.23、0.28、0.32、0.53、0.62 MPa,最大拉应力为 0.62 MPa(方案 E),显然 5 个方案中的最大拉应力小于矿石的抗拉强度。综合最大垂直应力和最大拉应力可以看出,尽管5 个方案的应力值均小于岩石强度,但对于方案 D、E中的应力值接近岩石强度。同时在方案 A 至方案 C中垂直应力与拉应力增长相对缓慢,从方案 C 到方案 D 垂直应力与拉应力增长幅度较大。3.3.2 位移分析矿房中的矿体一旦开采后,矿房上方的顶板将失去原始支撑,出现顶底板向开采空间挤压、顶底相互移近现象。因此,顶底板的移近量是采场空间安全性表 3 数值模拟力学参数Table 3 Mechanical parameters of numerical simulation类型岩石密度/(kg/m3)弹性模量/GPa泊松比黏聚力/MPa摩擦角/()抗拉强度/MPa单轴抗压强度/MPa矿石2 5208.20.271.82330.8731变质大理石2 70059.00.363.48421.4042花岗岩2 66074.00.334.23481.6058图 4 各方案的垂直应力云图Fig.4 Vertical stress nephogram of each scheme的重要指标。5 个方案的位移云图如图 7 所示。分析可知:随着矿房宽度增加,顶底板移近量逐渐增大。移近量最大值发生在矿房中部位置,对 5 个方案开采后的顶底板移近量进行了统计,结果如图 8 所示。由图 8 可知:随着矿房宽度增大,最大移近量分别为11.83、13.43、14.9、17.53、19.25 cm。除了方案 C、D541 柴 红等:缓倾斜厚大矿体采场结构参数优化研究与应用 2023 年第 11 期图 5 不同矿房宽度方案中的最大垂直应力Fig.5 Maximum vertical stress in different schemes图 6 不同矿房宽度方案中的最大拉应力Fig.6 Maximum tensile stress in different schemes存在微小波动外,顶底板移近量与矿房宽度近似呈线性关系。由上述分析可知:5 个方案的围岩基本能够保持稳定,都具备了一定的自支撑能力。方案中 D、E 的应力值趋近于岩石强度,故不建议采纳两者。综合考虑经济效益与安全性,认为方案 C(宽度为 34 m)为最佳方案。该方案充分利用了围岩自身承载能力,同时在一定支护的前提下,围岩稳定性相对可控。4 结 论(1)基于山东半岛某矿现场地质条件,利用Mathew 稳定图表法,计算得到在相对稳定情况下的采场结构参数,在综合考虑利用围岩自身支撑能力的基础上,给出了采场宽度合理取值范围为 28.36 42.17 m。(2)利用 FLAC3D软件对采场进行了数值模拟分析,不同宽度方案下应力场与位移场的数值模拟结果表明,34 m 采场宽度为最优参数。图 7 不同矿房宽度方案中顶底位移云图Fig.7 Nephogram of roof and floor displacement in different room widths schemes图 8 不同矿房宽度方案中顶底板移近量Fig.8 Displacement of roof and floor under different room width schemes参 考 文 献1 王新民,李洁慧,张钦礼,等.基于 FAHP 的采场结构参数优化研究J.中国矿业大学学报,2010,39(2):163-168.WANG Xinmin,LI Jiehui,ZHANG Qinli,et al.Optimizing mining stope structural parameters using a FAHPJ.Journal of China Uni-versity of Mining and Technology,2010,39(2):163-168.2 汪伟,罗周全,秦亚光,等.无底柱深孔后退式崩矿法采场结构参数优化J.东北大学学报(自然科学版),2016,37(4):578-582.WANG Wei,LUO Zhouquan,QIN Yaguang,et al.Stope parameters 641总第 569 期 金 属 矿 山 2023 年第 11 期optimization of non-pillar longhole retreat caving J.Journal of Northeastern University(Natural Science),2016,37(4):578-582.3 李启月,刘恺,李夕兵.基于协同回采的深部厚大矿体分段充填采矿法J.工程科学学报,2016,38(11):1515-1521.LI Qiyue,LIU Kai,LI Xibing.Sublevel filling method for a heavy orebody in deep mining based on collaborative stoping J.Chinese Journal of Engineering,2016,38(11):1515-1521.4 阮运斌.瓮福磷矿阶段嗣后采场结构参数优化J.中国矿业,2017,26(S1):385-389,398.RUAN Yunbin.The structural parameters optimization of stage back-filling of Wengfu Phosphateore MineJ.China Mining Magazine,2017,26(S1):385-389,398.5 吕明伟,王忠盛,汤永平,等.急倾斜厚大矿体采矿方法优化选择J.采矿技术,2021,21(4):5-9.L Mingwei,WANG Zhongsheng,TANG Yongping,et al.Optimiza-tion of mining methods for steep thick orebodiesJ.Mining Tech-nology,2021,21(4):5-9.6 赵增山,张舟,金爱兵,等.张家洼铁矿崩落转充填采矿工艺及结构参数研究J.金属矿山,2021(8):1-9.ZHAO Zengshan,ZHANG Zhou,JIN Aibing,et al.Study on the min-ing technology and structure parameters of Zhangjiawa Iron Mine from caving to fillingJ.Metal Mine,2021(8):1-9.7 罗周全,管佳林,冯富康,等.盘区隔离矿柱采场结构参数数值优化J.采矿与安全工程学报,2012,29(2):261-264.LUO Zhouquan,GUAN Jialin,FENG Fukang,et al.Stope structural parameters of panel isolation pillar numerical optimizationJ.Jour-nal of Mining&Safety Engineering,2012,29(2):261-264.8 周科平,王星星,高峰.秩和比法耦合数值模拟优化采场结构参数研究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.9 徐海.崩落转充填采矿方法粗糙模糊集理论优选J.中国矿业,2020,29(3):110-114.XU Hai.Optimization from caving to filling method based on rough fuzzy set theoryJ.China Mining Magazine,2020,29(3):110-114.10 陶干强,刘振东,任凤玉,等.无底柱分段崩落法采场结构参数优化研究J.煤炭学报,2010,35(8):1269-1272.TAO Ganqiang,LIU Zhendong,REN Fengyu,et al.Optimization re-search of stope structural parameters in sublevel caving with non-pillarJ.Journal of China Coal Society,2010,35(8):1269-1272.11 刘泽洲,吴伟伟,雷明.缓倾斜厚矿体采场结构参数优化研究J.矿业研究与开发,2020,40(9):6-10.LIU Zezhou,WU Weiwei,LEI Ming.Study on optimization of stope structural parameters of gently inclined thick orebodyJ.Mining Research and Development,2020,40(9):6-10.12 焦国芮,陈才贤,杨福斗,等.缓倾斜极薄矿体采场结构参数优化研究J.矿业研究与开发,2020,40(6):1-4.JIAO Guorui,CHEN Caixian,YANG Fudou,et al.Study on optimi-zation of stope structural parameters of gently inclined extremely thin orebody J.Mining Research and Development,2020,40(6):1-4.13 罗来和,夏国进,王沉,等.基于 FLAC3D的采矿方法优选及采场结构参数优化J.矿冶工程,2021,41(6):129-133.LUO Laihe,XIA Guojin,WANG Chen,et al.Optimization of mining method and stope structure parameters based on FLAC3DJ.Min-ing and Metallurgical Engineering,2021,41(6):129-133.14 周宝坤,陈晓青,田迎春,等.基于响应曲面法的崩落法采场结构参数优化J.金属矿山,2021(3):67-73.ZHOU Baokun,CHEN Xiaoqing,TIAN Yingchun,et al.Optimiza-tion of stope structure parameters by caving method based on re-sponse surface methodJ.Metal Mine,2021(3):67-73.15 陈顺满,吴爱祥,王贻明,等.基于响应面法的破碎围岩条件下采场结构参数优化研究J.岩石力学与工程学报,2017,36(S1):3499-3508.CHEN Shunman,WU Aixiang,WANG Yiming,et al.Optimization research on stope structure parameters in broken rock conditions based on the response surface methodJ.Chinese Journal of Rock Mechanics and Engineering,2017,36(S1):3499-3508.16 徐帅,安龙,李元辉,等.基于 SOM 的深埋厚大矿体采场结构参数优化研究J.采矿与安全工程学报,2015,32(6):883-888.XU Shuai,AN Long,LI Yuanhui,et al.SOM-based optimization of stope structural parameters of deep&large-sized orebodyJ.Jour-nal of Mining&Safety Engineering,2015,32(6):883-888.17 张亚伟,郭进平.缓倾斜极薄矿体空场法开采方案优化分析J.化工矿物与加工,2016,45(6):37-40.ZHANG Yawei,GUO Jinping.Optimization analysis of open stoping scheme for gently inclined and extremely thin ore bodies J.Chemical Minerals and Processing,2016,45(6):37-40.18 兰明,刘志祥,李夕兵.阶段嗣后充填采场结构参数的多目标多属性优化J.中南大学学报(自然科学版),2019,50(2):375-383.LAN Ming,LIU Zhixiang,LI Xibing.Multi-objective optimization and multi-attribute decision making on structural parameters of stage backfilling stope J.Journal of Central South University(Science and Technology),2019,50(2):375-383.19 陈霖,黄明清,唐绍辉,等.大直径深孔空场嗣后充填法采场结构参数优化及稳定性分析J.金属矿山,2022(11):44-51.CHEN Lin,HUANG Mingqing,TANG Shaohui,et al.Optimization and stability analysis of stope structure parameters of large-diameter deep hole open stoping with subsequent fillingJ.Metal Mine,2022(11):44-51.20 苗磊刚,牛园园,潘泱波.某矿尾砂胶结充填体强度演化特征及力学效应研究J.金属矿山,2023(3):52-58.MIAO Leigang,NIU Yuanyuan,PAN Yangbo.Study on strength e-volution characteristics and mechanical effect of cemented tailings backfill in a mineJ.Metal Mine,2023(3):52-58.741 柴 红等:缓倾斜厚大矿体采场结构参数优化研究与应用 2023 年第 11 期