基于诱导冒落开采的强制崩落区范围确定方法.pdf

Series No.568October 2023 金 属 矿 山METAL MINE 总 第568 期2023 年第 10 期收稿日期 2023-04-27基金项目 河北省自然科学基金项目(编号:E2020209009);华北理工大学博士科研启动基金项目(编号:28414099);华北理工大学教育教学改革研究与实践项目(碳中和专项)(编号:T-L2246)。作者简介 马姣阳(1985),男,讲师,博士。基于诱导冒落开采的强制崩落区范围确定方法马姣阳1,2 陈生鑫3 魏殿恩4 张庆嵩5(1.华北理工大学矿业工程学院,河北 唐山 063210;2.河北省矿业开发与安全技术重点实验室,河北 唐山 063210;3.北京科技大学土木与资源工程学院,北京 100083;4.北京交通大学土木建筑工程学院,北京 100044;5.东北大学资源与土木工程学院,辽宁 沈阳 110819)摘 要强制崩落区范围是否合理对诱导冒落法的冒落范围、过程及矿石回收影响重大。为此,以急倾斜破碎中厚矿体诱导冒落法开采条件为工程背景,采用数值模拟的方法,在确定崩落区采幅、高度及开采方式的前提下,构建动态开采模型,分析不同崩落区开采范围、埋藏深度条件下的诱导冒落矿体的致冒机理,形成强制崩落区确定方法。研究结果表明:结合崩落形状、应力分布状态和岩体破坏情况可有效确定强制崩落区,以及特定条件下强制崩落区范围最优形式;增加崩落隅角改变强制崩落区范围,能够有效地提高诱导矿层冒落能力;近崩落区位置存在复杂应力分布状态,崩落区顶板的垂直应力及水平应力的常值受深度影响呈现规律性阶梯型增长。研究结果可为诱导冒落采矿工程提供理论支撑。关键词 诱导冒落采矿法 急倾斜矿体 强制崩落区 采深 中图分类号TD853.36 文献标志码A 文章编号1001-1250(2023)-10-010-08DOI 10.19614/ki.jsks.202310002Method for Determining the Extent of Forced Collapse Zone Based on Induced Fall MiningMA Jiaoyang1,2 CHEN Shengxin3 WEI Dianen4 ZHANG Qingsong5(1.School of Mining Engineering,North China University of Science and Technology,Tangshan 063210,China;2.Hebei Mining Development and Safety Technology Laboratory,Tangshan 063210,China;3.School of Civil and Resource Engineering,University of Science and Technology Beijing,Beijing 100083,China;4.School of Civil Engineering,Beijing Jiaotong University,Beijing 100044,China;5.School of Resources and Civil Engineering,Northeastern University,Shenyang 110819,China)Abstract Whether or not the extent of the forced collapse zone is reasonable has a significant impact on the extent of the induced fracture process and ore recovery.To this end,the mining conditions of the induced fall method for sharply inclined and broken medium-thick ore bodies are used as the engineering background,and a dynamic mining model is constructed by numerical simulation to analyze the mechanism of induced fall ore bodies under the conditions of mining range and burial depth of different collapse zones,and to develop a method to determine the forced collapse zone by determining the width,height and mining method of the collapse zone.The research results show that the forced collapse zone can be effectively determined by combining the collapse shape,stress distribution state and rock damage,and the optimal shape of the forced collapse zone under specific conditions;increasing the collapse angle to change the scope of the forced collapse zone can effectively improve the a-bility to induce the ore layer to fall;the location of the near collapse zone has a complex stress distribution state,and the con-stant values of vertical stress and horizontal stress of the top plate of the collapse zone are affected by depth.The research re-sults can provide theoretical support for the induced fall mining project.Keywords induced fall mining method,sharply inclined ore body,forced collapse zone,mining depth 诱导冒落采矿法作为一种技术要求高、大规模开采方法,其生产能力甚至可以与露天采矿法媲美,是目前乃至未来几十年中开采低品位、深埋中厚以上矿体的最理想方法之一。该法实质是通过人为的开掘01措施工程1,促使诱导冒落区(简称诱导区或冒落区)矿体形成冒落条件,即矿体底部达到一定的暴露面积,在重力和内部应力共同作用下,冒落区矿体内原生节理扩展、贯通,最终以矿块的形式发生冒落,达到矿体回收和形成覆盖层的目的。该采矿技术历经西石门铁矿南区2、桃冲矿业公司盲矿体3-4、书记沟铁矿5、北洺河铁矿6等多个矿山的应用和完善已日臻成熟。20022005 年,北洺河铁矿利用第一分段进路回采提供的空间,诱导上部顶板矿石、围岩先后自然冒落,形成诱导冒落与强制崩落相结合的无底柱高效采矿方法,开创了投产 1 a 零 7 个月即达到设计生产能力的国内最好纪录。此后和睦山铁矿7、小汪沟铁矿8在生产中进一步完善并扩展了该采矿方法的应用范围,均实现了安全高效的开采目标,取得了良好的经济效益。伴随其应用取得成绩,该法的相关理论研究也有了进一步的进展,李海英等9-10利用诱导冒落技术控制挂帮矿地采岩移,根据采空区冒落过程,建立了采空区冒落诱发地表塌陷范围的计算方法。任凤玉等11-12针对复杂条件下不规则空区围岩冒落时空演化特征,提出了一种空区围岩冒落范围预测方法;通过非连续变形分析方法研究了不同边坡岩体结构下,诱导冒落法回采挂帮矿时引发的边坡失稳过程。胡颖鹏等13结合极限平衡分析方法,揭示出崩落法开采扰动下高陡边坡渐进破坏的规律,利用诱导冒落控制技术实现了塌陷坑有序承接边坡滑移体。伴随对冒落控制的研究逐步深入,该法聚焦于急倾斜破碎中厚矿体开采方面的试验研究,马姣阳等14在分析此类矿体可冒性基础上,提出了进路诱导冒落法开采方案,在处理近矿体中小空区基础上,取得了较好的工业指标。但在现场试验过程中,崩落区范围若选定不合理,极易引起周边岩体在回采中先于矿体冒落,增加混岩率。为此,以某铁矿急倾斜中厚矿体诱导冒落法开采条件为工程背景,在确定崩落区主要参数及开采方式基础上,采用数值模拟的方法,构建急倾斜破碎矿体的动态开采模型,通过改变崩落区剖面形式,研究不同埋深条件下崩落区开采过程中诱导冒落矿体的致冒机理,依此形成强制崩落区确定方法。1 动态连续开挖模型的构建1.1 工程背景该矿矿区内共有 8 条铁矿体,、号矿体规模较大,号矿体是北区的主矿体,是本次研究矿体。该矿体厚度 810 m,呈层状分布,倾角多为 7085,整体稳定性较差,矿石密度约为 3.56 t/m3,普氏硬度系数为 814,松散系数 1.75。回采过程中,由于矿体破碎,切割难以形成,在已施工工程中,切割工程不到位,从而影响后排炮孔爆破效果,形成了多个形状各异的小空区。这些小空区极易受岩层扰动、地应力影响,不断地发展变大并冒落垮塌积压更多的矿石,严重影响了矿山的正常回采工作,对矿山地下工作人员安全造成了一定的风险隐患。因此,采矿方法也由阶段矿房法改为诱导冒落法。但介于如何更有利地诱导冒落区矿体冒落回收,提高回采率,依旧需要进一步地确定崩落范围。1.2 崩落区采幅、崩落高度确定通过现场调研及临界冒落跨度理论计算,该矿体的可冒性良好,且临界冒落跨度小于 6.5 m,满足诱导冒落条件。此外,考虑矿体冒落后的采场矿石流动问题,结合回采过程中冒落矿石块度、湿度,通过结拱实验分析,不结拱的跨度约 9 m。因此,崩落区采幅取两者的最大值,不应小于 9 m15,故结合矿体厚度,模拟的采幅选择为 10 m。空区高跨比即为采空区截面高度与宽度之比,诱导冒落采矿过程中,拉底崩落空间需要为上方冒落矿体提供足够缓冲层及冒落空间,综合考量国内代表性的崩落法矿山选取数值模型的拉底高垮比为 2.5115,崩落高度定为 25 m。1.3 崩落区剖面主要形式为了探究诱导冒落矿石的最佳冒落效果,能够使矿体在最短时间内发生持续、稳定的冒落,从而可实现破碎矿体的安全高效回采。因此,在不改崩落高度的条件下,通过爆破技术对崩落剖面形状进行控制,即崩落下盘侧部分围岩,通过边孔角控制爆破范围,尽可能使诱导冒落区矿石较多处在预测冒落线内,同时减少崩落岩石基础上,增加隅角的尖锐程度,降低诱导冒落区矿体的支撑条件,改变其受到的应力状态致其冒落。此外,结合急倾斜空区冒落及顶板岩体应力分布规律,少崩落部分矿体,即通过改变上盘侧炮孔长度,使崩落区域上盘侧矿体,通过自身拉应力破坏达到自行冒落。由此本研究共设计 3 种崩落剖面见图 1。1.4 崩落区开采方式崩落区可采用从矿块中部切割拉槽或由矿体端部切割拉槽方式进行开采。为了更好地实施诱导冒落,应将切割工程尽可能布置在最易破坏矿体下方,并以其为中心向四周不断扩大崩落面积,利用爆破震动互相扰动因素,增加上层诱导冒落矿体内部裂隙程度,降 低 其 抗 拉 强 度,实 现 最 佳 的 矿 石 冒 落 效果16-17。鉴于该矿段整体相对规则,矿体破碎易崩,因此,选择采用从矿块中部切割拉槽的方式,从矿块中部开切做切割巷道,然后向两端同时开采方式,炮11 马姣阳等:基于诱导冒落开采的强制崩落区范围确定方法 2023 年第 10 期图 1 75矿体倾角的崩落剖面Fig.1 Collapse profile of orebody with 75 dip angle孔排距 22.5 m,见图 2。图 2 由矿块中部向两端开采方式Fig.2 Mining method from the middle to both ends of the block1诱导冒落矿体;2回采方向;3切割巷道;4切割槽1.5 模型的构建及不同采深的应力取值根据现有阶段的地质资料及动态模拟目标,建立矿体三维动态开挖模型。矿体呈层状分布,模拟倾角为 75。模型南北方向长 120 m,东西两端间距厚度为 60 m,整个模型高度为 85 m,矿块高度为 50 m,崩落区高度为 25 m,宽度为 10 m,模型由矿体中部向两端爆破开采,开采 4 排炮孔作为 1 步模拟开挖,即 10 m,分 5 步开采完成。为了开挖层运算结果较高的准确性,加强主要分析区域(空区围岩及顶板)的网格密度,其余部分网格适中,模型网格共计 205 250 个节点,338 483 个单元,如图 3(a)所示。结合已测地质资料及力学实验对模型进行赋值,模型参数见表 1。诱导冒落区矿体整体遵循拱形冒落机理,为了更好探寻3种崩落剖面下开采过程中顶板岩体冒落机图 3 动态开挖模型网络及主要监测点位Fig.3 Dynamic excavation model network and main monitoring points1-11-5、2-1、2-2、3-13-3、4-1监测点编号理,沿截面以拱形状在矿岩体内布置 10 个监测点,沿截面上以拱形状设置监测点 10 个,2个截面的主要监测点位见图 3(b)、(c),由此监测 3种方案下顶板岩体应力分布及破坏情况。为探究不同采深对崩落区开采过程中的影响,改变地应力,分析其对崩落区的影响性,深地平均水平应力与深地垂直应力均呈增加趋势。E.T.Brown 和E.Hoke18两位学者归纳总结了世界各地区地应力的测量结果,即深地垂直应力 v随深度 H 变化的规律:v=0.027H.(1)表 1 动态开挖模型材料参数Table 1 Dynamic excavation model material parameters矿岩体体积模量/GPa剪切模量/GPa质量密度/(kg/m3)内聚力/MPa内摩擦角/()抗拉强度/MPa底板2.961.83 1003.3942.51.3诱导岩体1.891.233 5600.88360.88矿体层2.181.83 5601.75371.25上盘围岩4.512.462 7003.82401.34下盘围岩4.732.82 7003.49391.28顶板 15.483.12 8003.45421.69顶板 25.873.22 6003.24411.5821总第 568 期 金 属 矿 山 2023 年第 10 期 依据大量理论分析及工程实际可知,随着采深的增加平均水平地应力与垂直地应力的比值 k 呈逐渐减小趋势,见式(2)19。k 值一般为 0.55.0,大多数为 0.81.5。不同岩性的硬岩、软岩等两者比值判别式不尽相同,但都满足该规律。国内常以 k=1 作为地下深部与浅部的分界线。k=h,avv,(2)h,av=h,max+h,min2.(3)赵德安等20基于大量国内实测数据,分别拟合了岩浆岩、变质岩、沉积岩等多种常见岩层的深地水平平均应力。E.T.Brown 等研究的深地垂直应力变化曲线和赵德安等研究的变质岩深地拟合曲线能够很好地解释初始应力的变化情况。为此,结合式(2)与式(3),作为增加模型初始边界应力条件的计算曲线,参考国内地质情况相似岩层的矿山资料21对数值模型进行不同深度的赋值,边界赋值情况见表 2。表 2 初始边界应力条件模拟赋值Table 2 Initial boundary stress condition simulation assignment深度/m垂直应力 v/MPa水平平均应力 h,av/MPak2005.410.7461.9940010.813.0681.2160016.219.441.2080021.623.761.101 00027271.02 基于连续开挖模型崩冒范围确定2.1 应力分布情况随开挖过程的变化不同类型崩落剖面模型在开挖过程中垂直应力的变化结果差别较小,且水平应力及最大剪切力随动态开挖规律性变化不明显。故选取常规崩落模型的截面监测的垂直应力结果进行分析,如图 4所示。图 4 常规崩落模型动态开挖过程截面垂直方向应力Fig.4 Conventional collapse model dynamic excavation process section vertical direction stress 由图 4 中监测点数据可知,随着崩落区的开采,所形成的空区上盘围岩及顶板的垂直方向应力减小,而两侧崩落区矿体和上下盘承压岩体的垂直应力急剧增加,承担开采部位所应承担的垂向应力。该现象在监测点 1-4(27.5,30,40)点较为明显。在第 13步开挖时,监测点 1-4(27.5,30,40)处作为承压矿区垂直应力逐渐增大。在第 4 步开挖后,该点对应的下方崩落区矿体被开采,不再作为承压矿体,垂向应力急剧下降。由常规崩落模型截面垂直方向应力云图及监测信息可知:在整个动态开挖的过程中,垂向应力重新分布,主要会形成应力降低区、应力增高区和原岩应力区 3 个应力区域22。开挖的卸载作用致使崩落区矿体及围岩压应力增大。开采形成的空区上方应力减少,底板甚至会出现应力集中现象23-24。由图 5 可知,空区上方竖向应力变化在水平分布上均出现 3 个区域:空区中心上端为相对受拉区,空区中心两端及下端为相对受压区。随着开采的进行,空区形状由“窄高”向”宽矮”过渡,应力集中系数逐渐增加,但趋势增加幅度逐渐减小。诱导冒落区矿体应力释放系数逐渐增大,崩落区未崩落矿体的承压能力随着空区扩大而逐渐达到极限25。空区临界失稳的破坏因素逐渐发生改变,由原来的崩落区侧壁切落向崩落区拉剪破坏转化26。图 5 常规崩落模型动态开挖过程截面垂直方向应力云图Fig.5 Stress clouds in the vertical direction of section during dynamic excavation of conventional collapse model在整个开挖过程中。承压岩体的垂直应力随着空区的开挖而逐渐增大。空区上盘在应力降低达到最小值时,甚至会反向出现拉应力。此刻未崩落区极易受拉致裂破坏。2.2 应力分布及扰动破坏情况受崩落剖面形状的影响 选取常规崩落、增加崩落隅角和减少部分崩落 3种不同的崩落剖面的模型实验方案的第五步动态开挖后,对截面的垂直应力、水平应力、剪切应力及扰动破坏云图进行比较分析。31 马姣阳等:基于诱导冒落开采的强制崩落区范围确定方法 2023 年第 10 期表 3 动态开挖中截面垂直应力及最大剪切应力变化Table 3 Variation of vertical stress and maximum shear stress in section in dynamic excavation剖面编号75动态开挖模型第 5 步截面应力云图及破坏情况垂直应力水平应力剪切应力破坏情况(a)(b)(c)注:同种性质应力的图例已做规范化处理。均为(a)模拟图例 对比表 3 中 75动态开挖模型第 5 步截面(a)、(b)、(c)3 种崩落剖面形式结合破坏分布分析可发现,在上盘侧多开隅角会提高其空区上方矿体的扰动破坏程度,有助于上部矿体诱导冒落,而在垂直及水平拉应力较大区域少崩落矿体,则会减缓其部分矿体受扰动的破坏情况。上盘围岩的破坏程度降低,不利于诱导冒落空区上方矿体。对比每次开挖后的新增破坏情况,可视作此次开挖步骤中矿岩的扰动破坏情况。由此可知不同形状崩落剖面的围岩和需诱导冒落矿体的扰动破坏情况有较大差别。增加崩落隅角导致其上部矿体失去支撑力,即在区域内增加了空区的水平跨度,提高了其冒落概率。减少崩落部分矿体更多的依靠水平及垂直拉应力的集中致使矿体破坏冒落。对于提高诱导矿层冒落能力来说,增加崩落隅角的剖面形式更优。2.3 应力分布受采深的影响200 m 开采深度下截面剪切应力的变化如图 6 所示。随着采矿深度增加,监测点最大剪切力监测值如图 7 所示。200 m 开采深度下 截面剪切应力显著增加,41总第 568 期 金 属 矿 山 2023 年第 10 期图 6 200 m 开采深度下截面剪切应力变化图Fig.6 Shear stress variation diagram of section at a mining depth of 200 m图 7 监测点受剪切应力随着采深变化Fig.7 Variation of shear stress with mining depth at the monitoring points随开挖的进行,剪切应力以“椭球形”向外部岩体扩张,剪切应力峰值集中于开挖区底部。由图 7 可见,监测点 1-1 在 5 步开挖过程中均在剪切滑移线范围内部,故没有出现剪应力峰值突增现象;监测点 1-2、1-3 分别在第 1 步、第 3 步开挖时位于剪切滑移线上,此时监测点的剪应力峰值会大于其他开挖阶段的峰值;监测点 2-1、3-1、4-1 在不同开挖阶段均靠近剪切滑移线,剪应力峰值会处于较高水平。剪切滑移线的范围是影响剪应力大小的主要原因。监测点最大剪切应力受与开挖区域间距的影响,处于滑移圈的剪应力明显大于其他范围,深度的增加相应地也增加了剪应力峰值的大小。监测点 4-1 在各个开挖阶段没有出现在剪切滑移线上的情况,最大剪切应力随采深的增加而增加,且增加幅度也随之增加。就垂直应力及水平应力而言,其常值受深度影响呈现规律性阶梯型增长,赵德安等20对实测地应力场分布规律统计分析也得出了相同的成果。评判诱导冒落采矿法的效果有冒落过程、冒落形式、控制冒落边界、边界残留矿量等指标。考虑最佳的开采崩落截面,主要有以下几个方面,以本参考矿山工程为例,需要在前期形成持续稳定的矿体冒落,并且尽量使得冒落矿体落于空区中部,便于放矿回收,形成一个良好的放矿椭球体形状。崩落剖面图 1(b)空区在前期便可实现稳定大量冒矿。而其后期由于隅角下方受剪切力较大更容易产生岩帮崩落。崩落剖面图 1(c)由于垂直应力等线值逐渐向下方向扩展,且弧形中心也逐渐向下方方向迁移,导致开采前期上盘岩体及上方需冒落矿体受扰破坏较小,整体不如崩落剖面图 1(b)出矿稳定。3 结 论(1)在诱导冒落区,矿体整体遵循拱形冒落机理。通过结合采幅、崩落高度和开采方式,分析不同崩落剖面形状条件下的垂直、水平和剪切应力分布及扰动破坏情况,可以确定特定条件下强制崩落区范围51 马姣阳等:基于诱导冒落开采的强制崩落区范围确定方法 2023 年第 10 期的最优形式。(2)增加隅角可改变强制崩落区范围,促使降低上部矿体所受支撑力,且增加了空区的水平跨度,提高了其冒落概率。水平及垂直拉应力的集中致使矿体破坏冒落,进而减少崩落部分矿体。在诱导冒落采矿法开采急倾斜破碎中厚矿体过程中,增加崩落隅角能够有效地提高诱导矿层冒落能力。(3)剪切滑移线的范围是影响剪应力大小的主要原因,处于滑移圈的剪应力明显大于其他范围,深度的增加相应地增加了剪应力峰值的大小。崩落区顶板的垂直应力及水平应力的常值受深度影响呈现规律性阶梯型增长。参 考 文 献1 何荣兴,韩智勇,任凤玉,等.诱导冒落处理采空区技术适用条件分析J.矿业研究与开发,2020,40(8):36-40.HE Rongxing,HAN Zhiyong,REN Fengyu,et al.Analysis on appli-cable conditions for the treatment of goaf with induced cavingJ.Mining Research and Development,2020,40(8):36-40.2 何鑫,任凤玉.西石门铁矿南区高效开采工艺研究J.金属矿山,2013(1):53-55.HE Xin,REN Fengyu.Study on the high efficient mining techniques in South Section of Xishimen Iron MineJ.Metal Mine,2013(1):53-55.3 周宗红,任凤玉,袁国强.桃冲铁矿采空区处理方法研究J.中国矿业,2005,13(2):17-18.ZHOU Zonghong,REN Fengyu,YUAN Guoqiang.Study on treatment method of the mined area in Taochong Iron MineJ.China Mining Magazine,2005,13(2):17-18.4 周宗红,任凤玉,袁国强,等.诱导冒落技术在空区处理中的应用J.金属矿山,2005(12):73-74.ZHOU Zonghong,REN Fengyu,YUAN Guoqiang.,et al.Application of induced caving technique in void area treatmentJ.Metal Mine,2005(12):73-74.5 任凤玉,李海英,任美霖,等.书记沟铁矿相邻空区诱导冒落技术研究J.中国矿业,2012,21(S1):378-380.REN Fengyu,LI Haiying,REN Meilin,et al.Technique of induced caving on adjacent mined-out areas in Shujigou iron mineJ.China Mining Magazine,2012,21(S1):378-380.6 任凤玉,韩智勇,赵恩平,等.诱导冒落技术及其在北洺河铁矿的应用J.矿业研究与开发,2007,27(1):17-19.REN Fengyu,HAN Zhiyong,ZHAO Enping,et al.Induced caving technique and its application in Beiminghe Iron MineJ.Mining Research and Development,2007,27(1):17-19.7 周宗红,任凤玉,王文潇,等.后和睦山铁矿倾斜破碎矿体高效开采方案研究J.中国矿业,2006,14(3):47-50.ZHOU Zonghong,REN Fengyu,WANG Wenxiao,et al.Study on the high efficient mining scheme of inclined fractured orebody in Houhe-mushan Iron MineJ.China Mining Magazine,2006,14(3):47-50.8 李楠,常帅,任凤玉.小汪沟铁矿露天转地下分区高效开采技术J.金属矿山,2014(5):21-23.LI Nan,CHANG Shuai,REN Fengyu.High efficient divisional min-ing technology during transition from open-pit to underground mining in Xiaowanggou Iron MineJ.Metal Mine,2014(5):21-23.9 李海英,任凤玉,严国富,等.露天转地下过渡期岩移危害控制方法J.东北大学学报(自然科学版),2015,36(3):419-422,447.LI Haiying,REN Fengyu,YAN Guofu,et al.Control method of rock caving hazard during transition from open pit to underground miningJ.Journal of Northeastern University(Natural Science),2015,36(3):419-422,447.10 李海英,任凤玉,赵云峰,等.金属矿空区冒落诱发地表塌陷范围的计算模型探讨J.矿冶工程,2020,40(5):5-8.LI Haiying,REN Fengyu,ZHAO Yunfeng,et al.Calculation model of range of surface subsidence caused by the caving in mined-out area of metal mine J.Mining and Metallurgical Engineering,2020,40(5):5-8.11 任凤玉,张晶,何荣兴,等.复杂条件下不规则空区围岩冒落时空演化特征研究J.金属矿山,2022(3):28-34.REN Fengyu,ZHANG Jing,HE Rongxing,et al.Study on the tem-poral and spatial evolution characteristics of surrounding rock ca-ving in irregular goaf under complex conditionsJ.Metal Mine,2022(3):28-34.12 任凤玉,谭宝会,付煜,等.诱导冒落法回采挂帮矿引发的边坡失稳J.东北大学学报(自然科学版),2019,40(2):273-277,283.REN Fengyu,TAN Baohui,FU Yu,et al.Slope failure caused by hanging-wall ores mining by induced caving methodJ.Journal of Northeastern University(Natural Science),2019,40(2):273-277,283.13 胡颖鹏,谭宝会,丁航行,等.露天地下楔形转接过渡诱导冒落控制边坡岩移危害技术研究J/OL.金属矿山:1-122023-03-08.HU Yingpeng,TAN Baohui,DING Hangxing,et al.Study on the technology of induced caving to control the hazard of slope rock movement in the open-pit underground wedge transition methodJ/OL.Metal Mine:1-122023-03-08.14 马姣阳,任凤玉,张庆嵩,等.基于“崩、冒”耦合效应的近矿体小空区处理方法J.金属矿山,2020(3):42-49.MA Jiaoyang,REN Fengyu,ZHANG Qingsong,et al.Method of near orebody and small cavity treatment for coupling effects based on breaking and cavingJ.Metal Mine,2020(3):42-49.15 马姣阳.急倾斜破碎中厚矿体进路诱导冒落法及其应用研究D.沈阳:东北大学,2017.MA Jiaoyang.Research and Application of Drift Induced Caving Method for Steeply Inclined Fractured Orebody with Medium Thick-nessD.Shenyang:Northeastern University,2017.16 吴季洪.我国充填开采技术发展现状与展望J.山西焦煤科技,2022,46(1):7-11.WU Jihong.Development status and prospect of filling mining tech-nology in ChinaJ.Shanxi Coking Coal Science&Technology,2022,46(1):7-11.17 张文静.采矿方法环境影响的经济比较评价J.科技创新与61总第 568 期 金 属 矿 山 2023 年第 10 期生产力,2015(9):116-118.ZHANG Wenjing.Economic evaluation and comparison on environ-mental impact of mining methodJ.Sci-tech Innovation and Pro-ductivity,2015(9):116-118.18 BROWN E T,HOEK E.Trends in relationships between measured in-situ stresses and depthJ.International Journal of Rock Me-chanics&Mining Sciences&Geomechanics Abstracts,1978,15(4):211-215.19 谢和平,周宏伟,薛东杰,等.煤炭深部开采与极限开采深度的研究与思考J.煤炭学报,2012,37(4):535-542.XIE Heping,ZHOU Hongwei,XUE Dongjie,et al.Research and consideration on deep coal mining and critical mining depthJ.Journal of China Coal Society,2012,37(4):535-542.20 赵德安,陈志敏,蔡小林,等.中国地应力场分布规律统计分析J.岩石力学与工程学报,2007,26(6):1265-1271.ZHAO Dean,CHEN Zhimin,CAI Xiaolin,et al.Analysis of distri-bution rule of geostress in China J.Chinese Journal of Rock Me-chanics and Engineering,2007,26(6):1265-1271.21 蔡美峰,乔兰,于波,等.金川二矿区深部地应力测量及其分布规律研究J.岩石力学与工程学报,1999,18(4):414-418.CAI Meifeng,QIAO Lan,YU Bo,et al.Results and analysis of in-situ stress measurement at deep position of No.2 mining area of Jin-chuan nickel mineJ.Chinese Journal of Rock Mechanics and Engineering,1999,18(4):414-418.22 陈勇.沿空留巷围岩结构运动稳定机理与控制研究D.徐州:中国矿业大学,2012.CHEN Yong.Study on Stability Mechanism of Rockmass Structure Movement and Its Control in Gob-side Entry Retaining D.Xuzhou:China University of Mining and Technology,2012.23 贺昌友,王进,王洪武.昭通铅锌矿采空区稳定性 FLAC(3D)数值模拟研究J.采矿技术,2008,8(5):41-44.HE Changyou,WANG Jin,WANG Hongwu.Zhaotong city lead-zinc mine goafs FLAC3D numerical stabilityJ.Mining Technology,2008,8(5):41-44.24 刘树新,孙会龙.重复