深部破碎矿体分段空场嗣后充填法采场结构参数优化.pdf

Series No.566August 2023 金 属 矿 山METAL MINE 总 第566 期2023 年第 8 期收稿日期 2023-02-20基金项目 中国有色集团“揭榜挂帅”项目(编号:2022JBGS01)。作者简介 刘鹏鹏(1998),男,硕士研究生。深部破碎矿体分段空场嗣后充填法采场结构参数优化刘鹏鹏1 王贻明1 周科礼2 熊国雄2 石大庆1(1.北京科技大学土木与资源工程学院,北京 100083;2.大冶有色金属集团控股有限公司,湖北 黄石 435005)摘 要 为确定深部破碎矿岩条件下采场合理结构参数,通过理论分析、现场测试以及数值模拟研究相结合的方法,对深部矿体采用分段空场嗣后充填法的采场结构参数进行分析与模拟。首先,通过对井下不同岩体进行取样实验与工程地质调查,获得不同岩体的岩石力学参数与节理裂隙参数,采用 RMR 法与 BQ 法进行岩体分级,开展岩体质量评价,评估矿岩稳定性。其次,采用 Mathews 图法,确定了铜绿山铜铁矿现有赋存条件下的采场结构参数合理范围。最后,利用 Flac3D 三维数值模拟软件,对开采过程进行仿真模拟,结合采矿过程中的应力、应变等的变化,分析采场的稳定性,验证采场结构参数的合理性。关键词 破碎矿体 采场结构参数 岩体分级 Mathews 图法 数值模拟 中图分类号TD853 文献标志码A 文章编号1001-1250(2023)-08-181-08DOI 10.19614/ki.jsks.202308022Optimization of Stope Structure Parameters of Sublevel Open Stope with Delayed Filling Method for Deep Broken OrebodyLIU Pengpeng1 WANG Yiming1 ZHOU Keli2 XIONG Guoxiong2 SHI Daqing1(1.School of Civil and Resource Engineering,University of Science and Technology Beijing,Beijing 100083,China;2.Daye Nonferrous Metals Group Holdings Co.,Ltd.,Huangshi 435005,China)Abstract In order to determine the reasonable structural parameters of the stope under the condition of deep broken rock,the stope structural parameters of deep orebody with sublevel open stope with subsequent filling method are analyzed and simulated by combining theoretical analysis,field test and numerical simulation.Firstly,the rock mechanics parameters and joint fracture parameters of different rock masses were obtained through sampling experiments and engineering geological survey of different underground rock masses.RMR method and BQ method were used for rock mass classification,rock mass quality e-valuation and rock stability evaluation.Secondly,the reasonable range of stope structural parameters under existing occurrence conditions of Tonglushan copper iron mine is determined by using Mathews diagram method.Finally,Flac3D numerical simula-tion software is used to simulate the mining process.Combined with the changes of stress and strain in the mining process,the stability of stope is analyzed and the rationality of stope structural parameters is verified,which can provide reference for similar mines.Keywords broken orebody,stope structural parameters,rock mass classification,Mathews stability diagram approach,numerical simulation 随着金属矿山的机械化水平逐步提高以及开采深度逐渐加深,矿山开采面临着许多全新挑战。比如:原来设计的采场结构参数无法匹配当下的机械化水平;深部矿体开采时所受到的地应力增大、赋存的工程地质条件复杂;深部矿体形状与早期勘探的结果相比变化较大等。这都会大大制约矿山进入深部的开采效率,矿山往往需要根据自身条件,推动深部开采的技改工程。其中最核心的问题之一是深部采场的结构参数1。针对深部矿岩特征对采场结构参数进行合理优化是深部技改工程的关键技术问题。针对该问题,姜伟明、朱训国、苏龙等2-4通过现场岩体条件对岩体质量等级采用单一评价方式或模型进行了讨论,评价岩体的质量等级;郭进平、胡高建、李胜辉等5-7利用 Mathews 稳定图表法,针对浅部矿体采场结构参数进行了优化设计;李闯、徐永文、周亚博、路燕泽等8-11采用 Flac3D 数值模拟软件对薄181矿体采场结构参数进行了优化研究,确定了采场结构参数的合理性。本文主要研究背景为铜绿山铜铁矿深部矿体,其平均厚度为 30 m,为厚矿体,在评价岩体质量时采用了两种评价方法。铜绿山铜铁矿地表标高+80 m,已开采至-605 m中段水平,开拓水平达到-800 m 左右,整体回采深度已达到 900 m 左右。随着开采深度的增加,岩性与工程地质条件变化大,采场应力赋存条件复杂,破碎带多,矿山迫切需要提高采矿效率与生产能力,原有采场结构参数设计难以满足深部后期工程的需求,亟需进行优化。基于此,本研究采用理论分析、现场测试、数值模拟研究相结合的方法,基于现场取样开展岩石力学实验与工程地质调查,对深部岩体质量进行分级评价,在此基础上,采用 Mathews 稳定图法与 Flac3D数值模拟,从解析解和数值解相结合验证优化采场结构参数,研究可为深部矿体高效开采提供理论依据和工程借鉴。1 岩体质量分级评价岩体质量分级主要针对铜绿山铜铁矿深部矿岩特征开展定性和定量评价,深部矿岩类型包括铜铁矿、矽卡岩和大理岩,这 3 种岩体涵盖了深部绝大多数矿岩类型。1.1 岩石力学参数研究根据对深部矿体进行钻芯取样,开展岩石力学基础实验,具体工作如下:在井下-545 m 中段与-605 m 中段不同地点取得铜铁矿石、上盘大理岩与下盘矽卡岩岩样,加工成 50 mm100 mm 的标准岩石试件,利用 MTS 岩石力学试验机等开展单轴压缩、三轴压缩、巴西劈裂等实验,测试的力学参数包括:岩体的密度、单轴抗压强度、弹性模量、泊松比、单轴抗拉强度、单轴抗剪强度(内聚力和内摩擦角)等。岩石力学实验的具体过程依据标准 DZ/T 0276.42015 进行,如图 1。表 1 为测试后得到的不同岩性的力学参数结果。图 1 单轴抗压强度实验Fig.1 Uniaxial compressive strength test1.2 节理裂隙参数研究采用详细线法对深部岩体的节理裂隙展开调查。表 1 不同岩性岩石力学参数测试结果Table 1 Test results of mechanical parameters of rocks with different lithologies矿岩类型密度/(g/cm3)抗拉强度/MPa抗压强度/MPa弹性模量/GPa泊松比内聚力/MPa内摩擦角/()铜铁矿4.0359.1858.83 48.833 0.27015.3843.50矽卡岩3.2097.9242.00 37.772 0.20911.9044.99大理岩2.6479.3977.50 59.795 0.2458.8843.65具体在-545 m 中段、-605 m 中段、-665 m 中段、-725 m 中段选取结构面发育状况较好的地点对不同岩性岩体进行工程地质调查,内容包括:结构面的倾向、倾角、节理密度、节理间距等,表 2 为调查得到的节理裂隙参数结果。表 2 节理裂隙参数调查统计结果Table 2 Statistical results of investigation of nodal fracture parameters矿岩类型节理密度/(条/m)结构面平均间距/m优势节理组数优势节理产状铜铁矿0.333.044547517571361420523矽卡岩0.432.303210315525615大理岩0.224.494219217822141825941.3 岩体质量分级评价结果根据实验获得的不同岩石力学参数与节理裂隙参数,采用 RMR 法和 BQ 法对铜绿山矿矿岩进行岩体质量分级,最终铜铁矿为 级岩体,岩体质量为差中,属于坚硬岩体且岩体较为完整;矽卡岩为级岩体,岩体质量为差 中,属于坚硬岩体且岩体较为破碎;大理岩为 级岩体,岩体质量为中良,属于坚硬岩体且岩体较为破碎。表 3 为不同岩体质量评价等级结果。表 3 不同岩体质量评价等级Table 3 Different rock quality assessment levels矿岩类型RMR 法BQ 法岩体等级质量评价铜铁矿差中矽卡岩差中大理岩中良 RMR 法和 BQ 法在对岩体质量进行评价时,充分考虑了岩体结构面的发育状况,采用这 2 种方法对281总第 566 期 金 属 矿 山 2023 年第 8 期铜绿山铜铁矿深部岩体进行质量评价,结果更加接近工程实际。2 采场结构参数优化Mathews 稳定图方法用于空场法采矿设计,是一种关于岩体质量、开采深度、采场尺寸和稳定性间的一种经验关系。铜绿山铜铁矿深部采用的采矿方法70%以上均为分段空场嗣后充填法,因此结合深部矿岩特性12-14,可以采用 Mathews 稳定图方法对采场结构参数进行合理验算。2.1 Mathews 稳定图法参数介绍应用 Mathews 稳定图方法需要计算 2 个参数:形状因子和稳定性系数。(1)稳定性系数 N。描述采场条件的稳定性系数 N 定义为N=A B C Q,(1)式中,Q为修正的 Q 系统分级值;A 为岩石强度因子;B 为节理方向调整参数;C 为重力调整因子。具体计算如下:按照 Mathews 稳定图方法的要求对 Q 值进行修正,将地应力影响因素 SRF 设为 1,地下水的影响因素 Jw=1,最终得到 Q(修正后的 Q 值)。岩石强度因子:A=0.1ci 20.112 5ci-0.1252 ci 10,(2)式中,c为完整岩石单轴抗压强度,MPa;i为平行于暴露面方向的采矿诱导应力,MPa。c的值可通过室内岩石力学实验测得,i值可利用数值分析技术求出。当数值分析得出的 i为拉应力时,取 A 值为 1。节理方向调整参数 B=0.21。为了保证采场的安全,取节理对采场最不利的条件下进行考虑,因此,取 B=0.2 来进行计算。重力调整因子 C=8-6cos,为暴露面与水平面的夹角。(2)形状因子 S。形状因子 S 定义为S=L L12(L+L1),(3)式中,L 为采场宽度,m;L1为采场长度,m。当暴露面长宽之比41 时,此时的形状因子基本保持不变,即此时暴露面的稳定性主要受暴露面宽度的影响。2.2 采场结构参数计算根据现场调查结果并结合铜绿山矿相关地质资料与分段空场嗣后充填法设计资料,确定铜绿山矿深部矿体平均倾角为 45,矿房平均厚度为 30 m,矿块垂直矿体走向布置。按照 Q 系统分级修正后,得到 Q值见表 4。其中,RQD 为岩体质量指标,Jn为节理组数;Jr为节理粗糙度;Ja为节理蚀变、充填机胶结程度;Jw为节理裂隙水折减系数;SRF 为应力折减系数。表 4 修正后的 Q 值(Q值)Table 4 Corrected Q value(Qvalue)矿岩类型RQD/%JnJrJaJwSRFQ铜铁矿806111113.3矽卡岩601211115大理岩764111119 按照式(1)计算岩体稳定性系数 N,结果见表 5。表 5 稳定性系数 N 值Table 5 Stability coefficient N value矿岩类别QABCN铜铁矿13.310.23.7610矽卡岩510.23.763.76大理岩1910.23.7614.27 Mathews 稳定图方法的设计过程以 2 个因子稳定性系数 N 和形状因子(或水力半径)S 的计算为基础,然后将这 2 个因子绘制在划分为预测稳定区、潜在不稳定区和崩落区的图上。稳定性系数代表岩体在给定应力条件下维持稳定的能力,形状因子或水力半径 S 则反映了采空区尺寸和形状。图 2 为稳定性系数 N 与水力半径 S 相关关系的 Mathews 图。图 2 稳定性系数 N 与水力半径 S 相关关系的 Mathews 图Fig.2 Mathews plot of the correlation between stability coefficient N and hydraulic radius S根据 Mathews 稳定图结合不同矿岩类型修正的Q值与稳定性系数 N 值确定了形状因子 S 值,见表 6。任何井下的暴露面均可认为是由两个方向的跨度组成,即认为是一个长方形。由于矿块沿走向布381 刘鹏鹏等:深部破碎矿体分段空场嗣后充填法采场结构参数优化 2023 年第 8 期表 6 形状因子 S 值Table 6 Shape factor S value矿岩类型QN稳定 S崩落 S铜铁矿13.3105.788.14矽卡岩53.764.236.26大理岩1914.275.858.22置,故矿体平均厚度即为矿块长度。根据式(3)结合形状因子 S 值,确定采场结构参数。矿体平均厚度为 30 m,根据铜绿山矿分段空场嗣后充填法可知,采场垂直走向布置。若仅考虑采场顶板跨度,在矿体中布置采场,当其岩体无支护时,采场顶板跨度需要限制在 18.81 m 之内,此时允许暴露面积为 564.21 m2;当其岩体被支护时,采场顶板跨度需要限制在 35.6 m 之内,此时允许暴露面积为1 067.93 m2。但考虑到矿体平均倾角为 45时,还要考虑采场上盘围岩的稳定性。根据表 6 可知,无论是无支护情况或者支护情况下,3 种岩性的形状因子 S 大小关系均为:大理岩铜铁矿矽卡岩。因此只需考虑矽卡岩为采场上盘时采场的稳定性是否满足要求。而当矽卡岩为采场上盘围岩,当其岩体无支护时,采场上盘围岩跨度需要限制在 12.01 m 之内,此时允许暴露面积为 360.3 m2;当其岩体被支护时,采场上盘围岩跨度需要限制在 21.49 m 之内,此时允许暴露面积为644.62 m2。3 基于数值模拟的采场结构参数优化研究利用 Flac3D 对分段空场嗣后充填采矿法进行模拟优化,主要针对无支护采场,分析采场稳定性,优化采场结构参数。铜绿山铜铁矿的分段空场嗣后充填采矿法,采准切割布置如图 3 所示。中段开采划分为 4 个分段,分段高度 12 m,回采顺序为从下至上,矿体回采一般划分矿块(即矿房和矿柱),矿块一般采用“隔一采一”,每分段最多同时回采 2 个矿房。矿房和矿柱长度为矿体平均厚度 30 m,宽度为 10 m。图 3 分段空场嗣后充填采矿法示意(单位:m)Fig.3 Schematic of sublevel open stope with subsequent filling method1阶段运输平巷;2分段平巷;3溜井;4分段凿岩巷道;5充填回风天井;6出矿联络道;7出矿进路;8受矿堑沟;9炮孔;10矿堆;11阶段斜坡道;12充填体3.1 数值模拟方案介绍目前铜绿山铜铁矿主要开采中段为-545 m 中段以及-605 m 中段,因此本次模拟对象选择为-545 m中段矿块回采。本次模拟主要从模拟各中段盘区出矿完成及充填后位移、应力及塑性区的变化。每个开挖步骤客观地反映了矿体逐步开采过程中,围岩应力叠加、岩体变形延续与破坏逐渐发展的力学进程和效果。本次模拟主要研究采场宽度对于采场稳定性的影响,根据 Mathews 图法确定的采场宽度范围,综合考虑了顶板与上盘围岩稳定性,设计模拟采场宽度分别为 15、12、10 m。其中采场宽度 10 m 为目前矿山生产正在使用的采场结构参数。采场长度为矿体平均厚度 30 m,分段高度 12 m。深部岩体都属于弹塑性材料,所以这里采用 Mo-hr-Coulomb 模型。3.2 数值模拟参数确定根据岩石力学物理实验所得的力学参数,采用Hoek-Brown 法15-16对实验参数进行参数折减,使实验所得力学参数更加接近实际的矿岩力学参数,而后对数值模拟模型进行参数赋值。充填体力学参数根据现场试验得出。所有数值模拟力学参数见表 7。3.3 模型边界条件和基本假设(1)对矿岩性质的假设。获得折减后的岩体参481总第 566 期 金 属 矿 山 2023 年第 8 期表 7 折减后的矿岩和充填体物理力学参数值Table 7 Discounted values of physical and mechanical parameters of ore rocks and fillings名称密度/(g/cm3)抗拉强度/MPa内摩擦角/()内聚力/MPa弹性模量/GPa泊松比充填体2.3600.24022.71.3203.5200.320铜铁矿4.0351.01536.524.62230.1890.270矽卡岩3.2090.72434.213.30125.5110.209大理岩2.6471.33737.544.99034.6540.245数后就可以假设矿体为各向同性、均质且符合摩尔-库仑弹塑性模型的介质,其自重应力场如下:v=H,(4)h=v1-vH,(5)式中,v为自重应力,MPa;h为自重应力的水平分量,MPa;v 为上覆岩层泊松比;为岩体的容重,g/cm3;H 为距离地表的深度,m。(2)对计算模型的假设。地下开采过程中各阶段巷道、斜井、联络道等工程对矿体开采的稳定性有一定的影响,但对矿山整体来说影响不大。从宏观上看,对采空区的稳定性影响很小,故在建模及数值模拟分析过程中,可以相应忽略其对工程的影响。(3)矿 房结构的简化。为模拟方便,已较大程度模拟矿体形态,同时对尖锐、零散区域进行了简化。3.4 模型介绍及模拟分析本次模拟重点为矿体开采及充填,矿体沿走向为模型 Y 方向。如图 4 为建立的铜绿山矿部分中段矿体的三维数值计算模型,模型尺寸为 240 m120 m180 m(长宽高),网格共计 146 520 格,节点数量为 155 250。模型顶部施加的荷载 22.8 MPa,为模拟上部 565 m 深覆岩体的自重荷载。图 4 数值模拟 Flac3D 模型Fig.4 Flac3D model numerical simulation3.4.1 矿房开挖应力分析如图 5图 7 为-545 m 中段四分段 2 个矿房采场取不同采场宽度的开采及充填后应力变化云图,空区开采后相同高度相邻区域矿岩应力变化相似,分别通过平行于矿体倾角方向进行切片分析及对矿体剖面切片,分析讨论了空区周边应力变化。图 5 矿房宽度 10 m 开采与充填最小主应力云图Fig.5 Minimal principal stress cloud map for excavation and filling of 10 m width of mine room图 6 矿房宽度 12 m 开采与充填最小主应力云图Fig.6 Minimal principal stress cloud map for excavation and filling of 12 m width of mine room 从一分段至四分段,在每分段 2 个矿房同时回采过程中,矿房宽度分别为 10、12、15 m 时,四周应力集中主要在顶板和侧帮,顶板所受到的应力值均小于5 MPa,侧帮所受到的应力值均小于 2.5 MPa,均在岩体承受范围内。并且在经过充填后,采场四周的应力集中均得到了改善。581 刘鹏鹏等:深部破碎矿体分段空场嗣后充填法采场结构参数优化 2023 年第 8 期图 7 矿房宽度 15 m 开采与充填最小主应力云图Fig.7 Minimal principal stress cloud map for excavation and filling of 15 m width of mine room 但由于回采方式采用“隔一采一”,必须考虑矿柱受力情况,为后期矿柱回采打好基础。根据不同宽度矿房回采和充填时,两矿房之间矿柱所受最大应力,获得关于矿柱承受最大应力与矿房宽度曲线图,如图 8 所示。根据图 8 可知,矿房宽度越大,两矿房之间矿柱所受应力也逐渐增大,但应力值均在 15 MPa 以内。并且通过充填后,矿柱所承受最大应力均有 1 3 MPa 的降低。3.4.2 矿房开挖位移分析由于当前模拟矿体开采,采空区暴露面积较大,在不同分段中,皆出现了较大的位移,为准确反映空区周边矿岩位移变化,通过垂直走向方向剖面图综合图 8 矿柱承受最大应力与矿房宽度曲线Fig.8 Curves of the maximum stress on the pillar versus the width of the mine room分析,-545 m 中段四分段矿体开采及空区充填位移变化如图 9图 11 所示。图 9 矿房宽度 10 m 开采与充填位移云图Fig.9 Displacement cloud map of mining and filling with 10 m width of mine room图 10 矿房宽度 12 m 开采与充填位移云图Fig.10 Displacement cloud map of mining and filling with 12 m width of mine room 在模拟 3 种矿房宽度回采时,从一分段至四分段,空区顶板区域出现了较大的位移,靠近空区顶板中心位置矿岩下沉量较小,靠近下盘位置的底板出现底鼓以及巷道壁出现片帮较大,受空区暴露面较大的681总第 566 期 金 属 矿 山 2023 年第 8 期图 11 矿房宽度 15 m 开采与充填位移云图Fig.11 Displacement cloud map of mining and filling with 15 m width of mine room影响。不同矿房宽度模拟回采时,矿房宽度与顶板、下盘底板的形变位移关系如图 12 所示。矿房宽度为12 m 时,顶板、下盘底板的形变位移最大,但 3 种矿房宽度下的形变位移量相差不到 10 cm。图 12 矿房最大位移与矿房宽度曲线Fig.12 Curves of the maximum displacement of the mine room versus the width of the mine room 根据位移云图可知,对空区采取充填措施后,充填体及周边均出现了较小的下沉,下沉量普遍小于1 cm,说明对该矿房采场采空区采取充填措施效果显著。4 结 论(1)根据在铜绿山铜铁矿井下现场取样实验与现场工程地质调查,通过 RMR 法与 BQ 法进行岩体分级,确定铜铁矿为 级岩体,属于坚硬岩体且岩体较为破碎;矽卡岩为 级岩体,属于坚硬岩体且岩体较为破碎;大理岩为 级岩体,属于坚硬岩体且岩体较为完整。(2)根据 Mathews 稳定图结合不同矿岩类型修正的 Q值与稳定性系数 N 值确定了形状因子 S 值,确定了铜绿山铜铁矿分段空场嗣后充填法深部采场宽度。为保证顶板稳定性,在无支护条件下采场宽度需小于 18.81 m;在支护条件下采场宽度需小于 35.6 m。为保证上盘稳定性,在无支护条件下采场宽度需小于 12.01 m;在 支 护 条 件 下 采 场 宽 度 需 小 于21.49 m。(3)采用 Flac3D 对铜绿山矿-545 m 中段的 2 个矿房采场,进行回采模拟。综合考虑顶板与上盘围岩稳定性,设计模拟采场宽度分别为 10、12、15 m。在采用 3 种采场宽度回采时,矿柱所承受最大应力均在岩体承受范围之内;顶板、下盘底板的形变位移相差均在 10 cm 左右。而且矿房充填后,矿柱所承受的应力与顶板、下盘底板的形变位移均有良好的改善。在这种条件下,若采用 15 m 采场宽度可以避免更多采切工程,提高矿山回采效率,并且采场稳定性相比采用 10 m 或 12 m 采场宽度也仅有少量降低。所以最终推荐采场宽度为 15 m。参 考 文 献1 董金奎,冯夏庭,张希巍,等.地下采场破碎岩体稳定性评价与参数优化J.东北大学学报(自然科学版),2013,34(9):1322-1326.HUANG Jinkui,FENG Xiating,ZHANG Xiwei,et al.Stability evalu-ation and parameter optimization on the fractured rock mass around underground stopeJ.Journal of Northeastern University(Natural Science),2013,34(9):1322-1326.2 姜明伟,刘焕新,王玺,等.基于 Sirovision 节理岩体调查及岩体质量分级与支护研究J.有色金属工程,2022,12(1):115-121.JIANG Mingwei,LIU Huanxin,WANG Xi,et al.Joint rock mass sur-vey and rock mass quality classification and support research based on SirovisionJ.Nonferrous Metals Engineering,2022,12(1):115-121.3 朱训国,夏洪春,王忠昶.工程岩体分级标准在深部巷道围岩分级中的应用及分析J.煤田地质与勘探,2017,45(2):118-125.ZHU Xunguo,XIA Hongchun,WANG Zhongchang.Application and analysis of standard of classification of engineering rocks in the clas-sification of surrounding rock of deep roadwayJ.Coal Geology and Exploration,2017,45(2):118-125.4 苏龙,刘爱华.考虑岩体质量分级的地下采场稳定性研究J.中国安全科学学报,2014,24(12):9-15.SU Long,LIU Aihua.Study on underground mining stope stability based on quality classification of rock massesJ.China Safety Sci-ence Journal,2014,24(12):9-15.5 郭进平,晏承园,张雯,等.分段矿房嗣后充填采矿采场结构参781 刘鹏鹏等:深部破碎矿体分段空场嗣后充填法采场结构参数优化 2023 年第 8 期数优化研究J.矿业研究与开发,2021,41(10):1-7.GUO Yanping,YAN Chengyuan,ZHANG Wen,et al.Optimization study on stope structure parameters of sublevel open stoping with backfillJ.Mining Research and Development,2021,41(10):1-7.6 胡高建,杨天鸿,胡忠强,等.基于 Mathews 稳定图等方法的多角度采空区群稳定性分析评价J.采矿与安全工程学报,2017,34(2):348-354.HU Gaojian,YANG Tianhong,HU Zhongqiang,et al.Mathews-sta-bility-method-based multi-angle analysis and evaluation of mined-out zones stabilityJ.Journal of Mining and Safety Engineering,2017,34(2):348-354.7 李胜辉,王立杰,刘志义,等.复杂破碎矿体试验采场稳定性分析及结构参数优化J.金属矿山,2021(8):41-45.LI Shenghui,WANG Lijie,LIU Zhiyi,et al.Stability analysis and structural parameter optimization of test stope for complex broken ore depositJ.Metal Mine,2021(8):41-45.8 李闯,李江,唐文平,等.铧厂沟金矿试验采场数值模拟稳定性分析J.黄金,2016,37(6):39-42.LI Chuang,LI Jiang,TANG Wenping,et al.Analysis of numerically simulated stability of the test stope in Huachanggou Gold MineJ.Gold,2016,37(6):39-42.9 徐永文,吴炳书,李官兵,等.缓倾斜极薄矿体采场结构参数优化研究J.金属矿山,2022(11):52-57.XU Yongwen,WU Bingshu,LI Guanbing,et al.Study on the optimi-zation of the stope structure parameters of gently inclined extremely thin orebodyJ.Metal Mine,2022(11):52-57.10周亚博.基于 Flac3D 数值模拟的深部采场结构参数优化及应用J.黄金,2022,43(5):51-55.ZHOU Yabo.Optimization and application of deep stope structural parameters based on Flac3D numerical simulationJ.Gold,2022,43(5):51-55.11 路燕泽,郭斌,杨志强,等.复杂采场结构参数数值模拟及采矿方法优化J.金属矿山,2022(8):69-74.LU Yanze,GUO Bin,YANG Zhiqiang,et al.Numerical simulation on complex stope structure parameters and optimization on mining methodJ.Metal Mine,2022(8):69-74.12 周科礼,尹东升,熊国雄.高分段嗣后充填采矿法在铜绿山矿的实践和应用J.中国有色金属,2018(S1):352-354.ZHOU Keli,YIN Dongsheng,XIONG Guoxiong.Practice and appli-cation of high-segment subsequent backfill mining method in Tonglushan Copper Mine J.China Nonferrous Metals,2018(S1):352-354.13 成涌,戴宏辉,周科礼,等.盘区机械化上向中深孔高分层充填法在铜绿山矿的应用J.采矿技术,2015,15(6):3-4,64.CHENG Yong,DAI Honghui,ZHOU Keli,et al.Application of me-chanized upward-to-medium-depth hole and high-layer filling meth-od in Tonglushan MineJ.Mining Technology,2015,15(6):3-4,64.14 戴宏辉.铜绿山矿采矿方法优化研究及应用J.采矿技术,2021,21(1):12-14.DAI Honghui.Research and application of mining method optimiza-tion in Tonglushan MineJ.Mining Technology,2021,21(1):12-14.15 王永增,孟磊磊,王润,等.齐大山铁矿露天边坡岩体力学参数确定J.有色金属(矿山部分),2021,73(1):98-102.WANG Yongzeng,MENG Leilei,WANG Run,et al.Determination of mechanical parameters of rock mass from Qidashan Iron Mine open-pit slopeJ.Nonferrous Metals(Mining Section),2021,73(1):98-102.16 尹幸乐,苏华友,雷国荣,等.基于 RMR 与 Hoek-Brown 准则的矿岩力学参数估算J.化工矿物与加工,2020,49(6):15-18.YIN Xingle,SU Huayou,LEI Guorong,et al.Estimation of mechan-ical parameters of ore and rock based on RMR and Hoek-Brown cri-terionJ.Industrial Minerals&Processing,2020,49(6):15-18.881总第 566 期 金 属 矿 山 2023 年第 8 期