深部金属矿超前序次释压机理与调控方法研究进展.pdf

深部金属矿开采已经成为金属矿床开采的重要组成。在深部“三高一扰动”复杂工程环境下,采动地压灾害日趋凸显,严重制约矿山安全高效开采。结合国内外深部采矿研究进展和多年工程经验,提出深部金属矿超前序次释压机理与调控方法,将采动地压调控分为基本地压调控和二次地压调控。基本地压调控含采矿方法选择、回采顺序优化和隔离矿柱留设,二次地压调控含卸压爆破、锚喷支护和释能支护,并就地压调控技术分别进行详细阐释。基本地压调控在采前即已完成,二次地压调控在开采中或开采后进行,调控策略选择基于采动地压灾害风险评估等级。地压监测技术是确保地压调控技术成功实施的关键,诸多现场工程实践表明多点位移计、微震监测和钻孔电视是可行的监测手段。随着数值模拟技术发展,矿山尺度精细三维数值模拟融合微震监测信息具有良好的应用前景,可以为深部矿山采动地压调控与防灾减灾提供借鉴。关键词 深部开采 采动地压 序次调控 地压灾害 释能支护 中图分类号TD325 文献标志码A 文章编号1001-1250(2023)-08-153-09DOI 10.19614/ki.jsks.202308019Research Progress of Advanced Sequential Pressure Relief Mechanism and Regulation Method for Deep Metal MinesZHAO Xingdong1 ZHU Qiankun1 DAI Bibo1,2 LI Huaibin3(1.Laboratory of Safe Mining in Deep Metal Mine,Northeastern University,Shengyang 110819,China;2.State Key Laboratory of Safety and Health for Metal Mines,Maanshan 243000,China;3.School of Safety Science and Engineering,Anhui University of Science&Technology,Huainan 232001,China)Abstract Deep metal mining has become an important part of metal deposit mining.In the deep complex engineering en-vironment of three high and one disturbance,mining induced ground pressure disaster has becoming increasingly prominent,which seriously restricts the safe and efficient mining of deep mines.Combined with the research progress of deep mining at home and abroad and many years of engineering experiences,the mechanism and control method of advanced sequential pres-sure relief in deep metal mines is put forward.The regulation and control of mining ground pressure is divided into basic ground pressure regulation and secondary ground pressure control.The basic ground pressure control includes the selection of mining method,the optimization of stoping sequence and the setting of isolated pillars.The secondary ground pressure control includes destress blasting,bolt-shotcrete support and energy-absorbing support.The paper also gives a detailed explanation on the ground pressure control technology.The basic ground pressure control has been completed before mining,and the secondary ground pressure control is carried out during/after mining.The selection of control strategy is based on the risk assessment level of mining induced ground pressure hazards.Ground pressure monitoring technology is the key to ensure the successful imple-mentation of ground pressure control technology.Many field engineering practices show that multi-point displacement meters,microseismic monitoring and borehole cameras are feasible monitoring technology.With the development of numerical simula-tion technology,mine scale fine three-dimensional numerical simulation incorporate with microseismic monitoring information has a good application prospect,which can provide reference for deep mine ground pressure control and hazard prevention and mitigation.351Keywords deep mining,ground pressure,sequence regulation of mining ground pressure,ground pressure hazards,ener-gy-absorbing support 深部金属矿开采已经成为金属矿床开采的重要组成,并且持续向更深处推进。我国开采深度达到或超过千米的矿山已达 60 座,未来 10 年内,我国三分之一的地下金属矿山开采深度将达到或超过千米1-2。与浅部开采相比,深部开采处于“三高一扰动”工程背景3,致使在采掘活动中遇到诸多采动地压灾害(诸如层裂、岩爆、冒顶),严重制约深部金属矿安全高效开采。由于采矿采掘工程是动态过程,采动应力也处于动态变化之中,再加上深部金属矿所处复杂地质条件,使得深部开采采动地压灾害防控变得日趋复杂和困难。当前,采动地压调控大多局限于局部地压调控,取得不少成功实践。然而随着深部开采的深入和发展,传统地压控制策略变得捉襟见肘,因此需要创新地压调控理念。赵兴东4提出基于采动地压灾害风险评估的调控策略、机理及调控方法,建立了采动地压序次调控理论框架。吴姗等5从矿山开发投资价值出发,从战略与战术规划两个维度探讨了深井矿山从整体到局部设计的过程,并重点分析了岩石强度时效性对于开展深井矿山整体规划的重要意义。于世波等6认为深部区域采矿时序是深部地压管控的战略方法之一,提出基于数值模拟与微震监测相结合的分析技术来评价深部区域采矿时序地压调控卸荷效果。本文首先阐释了深部金属矿超前序次释压机理及技术体系,从开采前、开采中和开采后三个阶段进行采动地压灾害风险评估及序次地压调控策略,即基本地压调控与二次地压调控。然后进一步论述采动地压灾害风险评估方法、基本地压调控和二次地压调控相关方法和深部采动地压监测技术,着重强调矿山尺度数值模拟和微震监测技术及其交互显示在矿山尺度地压监测与调控、采动应力分析和采矿顺序优化领域的前景。1 深部金属矿超前序次释压机理与调控方法对于深部采动地压研究,需要从矿体赋存空间形态的整体出发,科学、系统地研究矿山开采潜在地压灾害风险,包括地压灾害发生概率及其危害程度,与矿山采动过程应力场动态演化的关联性,并对矿山地质灾害进行系统和定量的分析和评估;依据深部采动地压灾害风险等级,对深部采动地压进行合理调控,以最大程度降低地压灾害发生频率和危害程度。所谓采动地压超前序次调控,依据开采前、开采中和开采后采动地压灾害评估结果,分别采取不同地压调控策略,消除或降低采动地压灾害等级,超前序次释压流程如图 1 所示。基本地压调控是核心,关系到矿山开采全局地压调控成败,主要包括采矿方法选择、回采顺序优化和隔离矿柱留设,在开采前即已完成。二次地压调控是关键,在基本地压调控后,根据采掘工程揭露的矿岩实际情况,进行详细地压灾害评估,采取合适的地压调控方法,消除或降低局部高应力区或岩爆灾害区危险,主要包括卸压爆破、锚喷支护和释能支护,其是在开采中或开采后进行。图 1 超前序次释压流程Fig.1 Technical process of advanced sequential pressure relief 2 深部采动地压风险评估2.1 采动地压灾害风险评估基础资料对于深部采动地压灾害风险评估,基础资料搜集大致可分为 3 个阶段:开采前、开采中和开采后阶段,采动地压风险评估所需基础资料如图 2 所示。勘探与采矿设计阶段完成工程地质、采矿设计,同时还要完成矿岩稳定性评估、采矿顺序优化与采动应力分析及地压灾害风险初步评估,优化采矿设计方案、回采顺序和开拓工程布置。开采阶段完成工程地质资料补充、采空区处理、围岩变形破坏监测、支护设计与评估、开采期间微震监测等工作和局部地压灾害风险详细评估,完善局部地压控制措施。开采后主要进行微震数据、充填体稳定性监测数据及地表沉降数据等资料搜集和采空区地压风险及对地表和后续开采影响的评估。2.2 采动地压灾害风险评估方法深部开采围岩破坏主要分为结构面控制型和应力驱动型破坏形式。对于不同地压显现形式,需要采取相应的地压控制方法。地压评估采用地压灾害风险评估定量矩阵。以采场地压灾害风险评估为例,评价流程为:工程地质调查,确定岩体特性、应力边界条件、采场形状和支护;确定潜在的破坏模式及其诱发因素;确定损伤破坏范围及破坏的概率值;451总第 566 期 金 属 矿 山 2023 年第 8 期图 2 采动地压风险评估基础资料Fig.2 Basic data for mining induced ground pressure hazard assessment 确定采场围岩失稳造成的经济损失;绘制采场稳定性风险评估定量矩阵。在分析采场冒落型破坏时,需要统计冒落体积、分布特征,以确定采动地压灾害风险发生概率值。在分析采动应力导致围岩损伤破坏时,根据二维或三维弹性、弹塑性数值模拟或者经验方法确定采动应力损伤范围,根据数值模拟确定破坏概率(PoF)的方法有:点估计方法、响应面法和响应影响因子法。采动地压灾害风险=PoF 潜在破坏的后果,(1)式中,潜在破坏的后果以采动地压灾害造成的损失来计算。则可绘制出相应的地压灾害风险评估定量矩阵,如表 1 所示。表 1 采动地压灾害风险评估定量矩阵(改自文献7)Table 1 Quantitative matrix of mining induced ground pressure hazard assessment(modified from 7)发生的概率灾害风险不显著较小中等较大灾难性几乎肯定发生低中等高极高极高很可能发生低中等高高极高可能发生低低中等高高不太可能发生低低中等中等高罕见低低低中等中等 依据上述工作流程,完成矿山采动地压风险评估工作后,即可应用风险评估矩阵对矿山采动地压灾害进行评估,然后将采动地压灾害空间信息可视化,以便在生产实践中及时采取相应的地压调控措施。3 深部金属矿超前序次释压技术体系3.1 基本地压调控3.1.1 采矿方法地下采矿方法是指从地下矿床的矿块里采出矿石的过程,通过矿床开拓、矿块的采准、切割和回采 4个步骤实现8。地下采矿方法分类繁多,以地压管理方法为依据,分为矿柱支护采矿法、人工支护采矿法以及无支护采矿法(图 3)。不同采矿方法会在矿体近场和远场产生不同类型和程度的岩体力学响应9。从图3 可以看出,岩体采动位移与应变能储能特性呈负相关,无支护采矿法围岩会产生较大的位移,但围岩储存的应变能较低;矿柱支护围岩会产生较小的位移,但围岩储存的应变能较高,极易产生应力驱动型地压灾害,诸如层裂、屈曲、岩爆。对于深部开采金属矿山,充填采矿法是国内外主流选择,其原因为:改善围岩受力条件;填充结构面,胶结岩体形成完整结构;让压作用,减缓采场围岩地压释放和围岩积聚能量释放速率10-11。当然,采矿方法还受矿体形态、围岩性质、矿石品位、环保、产能、经济等因素影响。图 3 地下采矿方法与采动岩体响应9Fig.3 Underground mining methods and response of rock mass93.1.2 回采顺序地下矿床开采及采空区处理过程中,不同回采顺序致使围岩经受复杂的加卸载过程,进而会导致岩体产生不同类型的力学响应特征,选择合理的回采顺序是控制深部采矿地压灾害的有效途径之一。设计回采顺序核心思想是以压力拱理论和应力转移理论为指导,确保开采活动处于低应力区。图 4 所示为开挖顺序下采场周边高应力区演化规律,可以看出采场 1551 赵兴东等:深部金属矿超前序次释压机理与调控方法研究进展 2023 年第 8 期开挖后其周边存在宽度接近采场宽度的卸压区,采场周边同时存在一高应力区域;当回采采场 5 时,采空区两侧卸压区宽度基本未变,而顶底板出现大范围卸压区;从模拟结果看采空区周边卸压区近似呈圆形。目前,常用的深部矿体开采回采顺序有盘曲交错式回采、交错式回采、金字塔形阶梯式回采、倒阶梯式回采、不连续走向回采顺序和连续回采顺序4。在矿山回采顺序研究方面,工程上比较常用的研究方法主要有:工程类比法(即经验法)、数值分析法等。在数值分析方面,针对矿山尺度采矿顺序优化,大型边界元软件(如 MAP3D 软件)具有独特的优势,其具有单元个数少、计算速度快、数据准备简单、建模方便等特点,且相同离散精度下,边界元法解的精确度优于有限元法12。图 4 采动条件下高应力区演化规律Fig.4 Evolution law of high stress area under mining conditions 3.1.3 隔离矿柱在深部采矿过程中,每隔一定垂直距离留设一个或多个特定厚度的水平隔离矿层,支撑上下盘围岩,确保空区围岩整体稳定,防止空区上部发生大规模塌落时对下部采场产生的动力冲击及次生气浪对工作人员造成伤害,同时也起到降低采空区围岩高应力区分布范围和控制地表沉降的作用13。隔离矿柱的留设必须进行严格岩石力学计算以确定其合理位置及尺寸。判断隔离矿柱留设是否合理的依据是:当采动应力作用到该矿柱上时,隔离矿柱不会产生变形、屈服和破坏4。隔离矿柱可与采矿顺序、采场充填协同应用以增强地压调控效果。3.2 二次地压调控3.2.1 卸压爆破卸压爆破是金属矿山地下开采常用的调控高应力分布和防治岩爆灾害的措施,与水力压裂、注水、卸压孔和 CO2相变致裂技术相比,卸压爆破具有卸压范围大、效率高的优点,比较适合金属矿开采领域。其原理是在高应力岩体内利用爆破产生大范围爆生裂隙区,使岩体力学参数弱化,并将高应力转移至邻近岩体结构。图 5 所示为巷道掘进面卸压爆破前后垂直应力分布示意图,可以看出,卸压爆破后掘进工作面前方垂直应力会显著降低且峰值应力点向工作面前方深部转移14。Crouch15研究了卸压爆破前后应力和应变条件,并提出卸压爆破的亚临界、临界和超临界程度(图 6)。图 5 巷道掘进工作面爆破卸压应力调控机制14Fig.5 Regulation mechanism of destress blasting in the working face14 数值模拟方法是设计或优化卸压爆破方案的有效手段,其中卸压爆破后岩体力学参数取值方法是关键,为此诸多专家学者对此进行研究。Blake16给出卸压爆破后岩体弹性模量 Edb建议计算方法:Edb=E,(1)式中,为岩石破碎系数,取值范围为 01。Tang 等17给出卸压爆破后岩体的泊松比和应力值计算方法:db=(2-),(2)ijdb=(1-)ij,(3)式中,为应力耗散系数,取值范围为 01,反映爆破卸压区应力释放程度;ij为初始应力状态。Saharan 等18提出卸压爆破的本构方程:651总第 566 期 金 属 矿 山 2023 年第 8 期图 6 爆破卸压程度(改自文献15)Fig.6 Degree of destress blasting(Modified from 15)Eidb=iE,(4)idb=(1-i)i,(5)式中,i 表示主应力方向。该本构方程假设岩体为正交各向异性材料。此外,应用微震手段可以有效地评价卸压爆破效果19。3.2.2 锚喷支护国内外工程实践表明,锚喷支护可用于深部矿山岩体稳定性控制20-21。深部开采围岩控制应能够提高岩体强度和峰后强度,同时使支护结构形成整体,防止岩块冒落或滑落,共同抵抗围岩进一步变形破坏,同时能够承受一定程度的大变形。深部井巷、采场围岩控制设计方法多以经验或工程类比方法为主,或者采用悬吊理论、组合梁理论、压缩拱理论、组合拱理论、最大水平应力理论等进行设计。本研究提出以下支护设计思路(图 7):进行工程地质调查和基础岩石力学实验,构建精细数值模型;综合应用Q、RMR、GSI、MRMR 等岩体质量分级方法评估岩体质量等级,估算岩体力学参数;根据经验图表初选支护方案,校核支护能力,并应用数值模拟进行优化分析;现场支护施工,施工过程进行严格质量控制;综合地压监测,综合应用多种地压监测手段(多点位移计、钻孔应力计和钻孔电视等)监测评估支护效果,及时发现存在的问题和不足,及时反馈并再次优化支护设计方案。图 7 深部开采围岩支护设计思路Fig.7 Flowchart of surrounding rock support for deep mining3.2.3 释能支护对于深部开采金属矿山,受地质构造、高地应力和强采动应力影响,会出现不同程度的动力破坏形式,普通支护难以达到理想的支护效果(图 8),而释能支护是针对具有岩爆倾向性巷道支护的有效支护方式。释能支护是通过释能支护结构释放岩爆灾害产生的动能,因此释能结构要具备自身能产生较大的变形或能够产生足够的滑移变形的特点,同时自身还要具备高强度,在承受动载时能保持结构完整性。图 8 岩爆造成支护结构破坏22Fig.8 Support structure damage caused by rock burst22 目前,释能锚杆已经成为世界深井开采矿山的重要支护形式。通常释能锚杆自身具有高强度特点,具有较高的锚固力,能够抵抗动载荷多次冲击,同时在承受动载荷时能够产生杆体大变形或杆体滑移,这样才能有效释放岩爆灾害产生的动能,确保围岩稳定。当前世界上主流释能锚杆见表 2。与传统锚杆相比,释能锚杆通常具有较高的释能能力,具有高拉拔力和静载拉拔位移。岩爆等动力冲击巷道围岩释能支护设计过程为:评估潜在岩爆震级及释放能量值;选择合适的支护结构与支护参数,相应支护结构改为释能支护结构;计算支护系统支护能力与岩爆等动力冲击巷道围岩支护需求;逐一校验载荷、释能量和位移指标达到设定安全系数值;现场施工并观测支护效果,若存在问题应及时调整支护结构与支护参数,并对其进行载荷、能量与位移等指标的再校核,直至达到设定安全系数值。此外,也可以根据 Q 系统岩体质量分级结果,应用图 9 选择合理释能支护方式。4 深部采动地压监测现场工程应用表明,对于局部地压监测,钻孔应力计应用于深部硬岩应力监测效果不太理想,单(多)点位移计和钻孔电视可以更好地用于监测围岩变形和岩体内部破裂程度。图 10 所示为红透山铜矿751 赵兴东等:深部金属矿超前序次释压机理与调控方法研究进展 2023 年第 8 期表 2 释能锚杆类型及其典型技术特征Table 2 Types and typical technical characteristics of energy absorbing bolts释能锚杆类型主要技术特征Cone 锚杆23锚杆一端为圆锥形体,在锚杆杆体喷涂脱脂剂;锚固力为 160 kN 以上,静载拉拔位移 200 mm,总释能量 1025 kJ。Durabar 锚杆24对锥体锚杆进行了改进,在光滑杆体设计几个褶皱,锚杆的尾部设计成一个光滑的圆环,锚固力为 120 kN,静载拉拔位移600 mm,总释能量 48 kJ。Hybrid 锚杆25为管缝锚杆与螺纹钢锚杆组合而成,锚固力 160 kN,静载拉拔位移 120 mm 以上。Swellex 锚杆26是一种水力膨胀锚杆,注入高压水致使锚杆杆体变形,锚杆管壁与锚杆孔壁产生高摩擦力锚固岩体;锚固力为 200 kN,静载拉拔位移 80 mm,总释能量 1829 kJ。Garford 锚杆27由圆钢、锚头及粗牙螺纹钢套组成,采用树脂锚固,锚固力 190 kN,静载拉拔位移 390 mm,总释能量 2733 kJ。Roofex 锚杆28由锚固端和圆钢组成,采用树脂锚固,锚固力 80 kN,静载拉拔位移 200 mm;总释能量 1227 kJ。D 锚杆29-31由圆钢带一定数量的具有一定间隔的锚固点组成,锚固力 200 kN,静载拉拔位移 100120 mm,总释能量 3639 kJ。Yield-Lok 锚杆32杆体直径 17.2 mm,端部布置有若干锥形体,锥形体外部套有高分子涂层;静载拉拔位移 175230 mm;总释能量 16.4 kJ。VersaBolt33锚杆有2 种型号直径(20 mm 和22 mm),锚固器的位置和数量可以根据需求加工;静载拉拔位移86.55205.72 mm;总释能量37 kJ。He-Bolt34恒阻大变形锚杆由恒阻装置、连接套、杆体、托盘、球垫和螺母组成;恒阻体通过连接套与杆体连接成一体,恒阻套管内外表面均有螺纹;该锚杆通过杆体在连接套内滑动和杆体弹性变形释能(伸长量达 1 000 mm),锚杆具有恒阻特性;总释能量高达 200 kJ。MCB-Superbolt 和Versa-Superbolt35和 Cone 锚杆形式一样,锚杆端部有一屈服模块(滑动或杆体伸长),区别在于杆体外部有波纹状套管,端部与 MCB 锚杆相同的是 MCB-Superbolt,端部与 Versa 锚杆相同的是 Versa-Superbolt;静载拉拔位移约为278 mm;总释能量分别高于50 kJ和 70 kJ。J 锚杆36J 释能锚杆由锚固(阻尼)模块、变形模块、搅拌模块、锚固端四部分组成;静态拉拔位移 15.215.3 mm;总释能量 24.046.5 kJ。图 9 基于 Q 系统的释能支护设计经验图表37Fig.9 Support design chart for rock burst rock mass based on Q-system37 某斜坡道在周边采场分层回采期间围岩内部破裂规律,可以看出随着开采的进行,斜坡道岩体表面层裂深度逐渐增加,至第 20 分层时层裂剥落深度接近0.6 m,而此时深部岩体破裂位置距初始斜坡道表面深度达 1.01.1 m。而对于大范围地压监测,微震监测技术具有无可比拟的优势,当前,微震监测已成为分析采动岩体破裂失稳和岩爆预警评估的有效工具。典型的微震监测系统由若干传感器、接线盒、授时装置、网络通信、数据采集器和配套处理软件组成,图 11 所示为 ESG图 10 采动岩体破裂深度演化规律38Fig.10 Evolution law of fracturing depth of rock mass38 微震监测系统。Ge39提出矿山微震监测系统设计的851总第 566 期 金 属 矿 山 2023 年第 8 期5 个原则:监测范围应当为三维空间;传感器不仅要安装在目标研究区域,还要布置在研究区域周围岩体;传感器均衡布置;避免出现传感器二维布置,会导致垂直于传感器阵列面精度降低;仔细研究背景噪声和微震信号的频率分布特征以实现高精度滤波。在后处理过程中,微震事件的统计学规律、空间分布特征、b 值演化规律、能量和震源机制解均会提升对采动岩体响应的理解程度40-44。图 11 典型微震监测系统45Fig.11 Typical microseismic monitoring system45 在油气开采领域,地表微地震监测系统已被广泛应用于水力压裂特征监测和诊断。微地震系统记录特定的频率范围内的事件以覆盖具体工程应用相关的主频率,并据此将其分为实验室声发射系统、微地震监测系统、区域短周期地震台网、宽频地震台网和长周期全球台网46。微地震监测系统范围通常用于观测和分析数千公里范围内震级小于 0 的事件。微地震系统通常采用被动记录地表震动的地震检波器,成本较低,典型的频率为 5 Hz 或更高。从图 12 可以看出,注入(油气水力压裂、地热开采等)诱发地震在震级、断裂尺寸和角频率等关键参数与采矿诱发声发射和地震活动有大范围重叠,说明地面微地震系统同样适用于采矿诱发岩体破裂监测和分析领域,从而为采动诱发岩体破裂过程和采动地压灾害防控提供新思路。与常规微震监测系统相比,地表微地震监测具有成本低廉、抗(机械、电气等)干扰性强、覆盖范围大、维护成本低、迁移性好等优点;其缺点是垂直方向定位精度相对较低,但从地震统计学上来看,对分析采动灾害形成过程及其防控是有效的。地表微地震监测系统常用的定位算法是相似加权算法47,该方法具有时效性强、精度高等优点,可以实现微地震实时监测。图 12 不同尺度岩石破裂关键特性及参数48Fig.12 Key characteristics and parameters associated with rock fracture at different scales48 随着数值模拟技术发展,矿山尺度高效数值模拟技术已经出现。大型边界元软件 MAP3D 是目前分析矿山尺度采动岩体响应的高效分析工具,其具有建模简便、计算速度快、可进行并行计算的优点,而最具特色之处在于数值模拟微震交互功能49。能够在完成数值模拟分析时导入微震监测数据,实现应力、岩体破坏等与微震监测数据的高度融合交互显示(图 13),可以有效分析开采高应力分布区范围及评估潜在地压灾害,也可为下一步回采顺序设计提供精准指导。图 13 矿山尺度采动应力模拟(左)及与微震监测交互显示(右)49Fig.13 Simulation of mining-induced stress at mine scale(left)and interactive display with microseismic monitoring information(right)49951 赵兴东等:深部金属矿超前序次释压机理与调控方法研究进展 2023 年第 8 期5 结 论(1)提出深部金属矿超前序次释压机理与调控方法,依据开采前、开采中和开采后采动地压灾害评估结果,分别采取不同地压调控策略,消除或降低采动地压灾害等级。基本地压调控包括采矿方法选择、回采顺序优化和隔离矿柱留设。二次地压调控包括卸压爆破、锚喷支护和释能支护。(2)给出地压灾害风险评估所需资料及评估方法和流程。(3)详细阐释基本地压调控和二次地压调控核心技术及其应力调控机制。(4)根据工程实践经验,多点位移计、钻孔电视和微震监测均为有效的地压监测手段。随着数值模拟技术发展,矿山尺度采动应力模拟与微震监测信息融合显示,可有效辨识高应力区域及评估地压调控效果。参 考 文 献1 蔡美峰,薛鼎龙,任奋华.金属矿深部开采现状与发展战略J.工程科学学报,2019,41(4):417-426.CAI Meifeng,XUE Dinglong,REN Fenhua.Current status and devel-opment strategy of metal minesJ.Chinese Journal of Engineering,2019,41(4):417-426.2 赵兴东.超深竖井建设基础理论与发展趋势J.金属矿山,2018(4):1-10.ZHAO Xingdong.Basic theory and development trends of ultra-deep shaft constructionJ.Metal Mine,2018(4):1-10.3 何满潮,谢和平,彭苏萍,等.深部开采岩体力学研究J.岩石力学与工程学报,2005,24(16):2803-2813.HE Manchao,XIE Heping,PENG Suping,et al.Study on rock me-chanics in deep mining engineeringJ.Chinese Journal of Rock Mechanics&Engineering,2005,24(16):2803-2813.4 赵兴东.深部硬岩矿床采动地压与控制M.北京:冶金工业出版社,2019.ZHAO Xingdong.Deep Mining Induced Stress and Controlling Meth-od at Hard Rock SepositsM.Beijing:Metallurgical Industry Press,2019.5 吴姗,杨小聪,郭利杰.高应力环境下深部金属矿整体规划的思考与展望J.煤炭学报,2019,44(5):1432-1436.WU Shan,YANG Xiaocong,GUO Lijie.Thinking and prospect of o-verall planning for deep metal mine in high stress environmentJ.Journal of China Coal Society,2019,44(5):1432-1436.6 于世波,杨小聪,原野,等.深部区域采矿时序的地压调控卸荷效应研究J.黄金科学技术,2020,28(3):345-352.YU Shibo,YANG Xiaocong,YUAN Ye,et al.Research on destress effect of ground pressure control for the time-space mining sequence at depthsJ.Gold Science&Technology,2020,28(3):345-352.7Wael ABDELLAH,Hani S MITRI,et al.Geotechnical risk assess-ment of mine development intersections with respect to mining se-quence J.Geotechnical and Geological Engineering,2014,32:657-671.8 煤炭科学研究总院北京开采研究所.地下开采现代技术理论与实践M.北京:煤炭工业出版社,2002.Beijing Mining Research Institute,China Academy of Coal Sciences.Theory and Practice of Modern Underground Mining TechnologyM.Beijing:China Coal Industry Publishing House,2002.9BRADY B H G,BROWN E T.Rock Mechanics:for Underground MiningM.3rd Ed.Dordrecht:Springer,2004:355.10 王明旭.胶结充填体与围岩相互作用机理及变形演化规律研究D.武汉:武汉科技大学,2018.WANG Mingxu.Study on the Interaction Mechanism and Deforma-tion Evolution Law of Cemented Filling Body and Surrounding RockD.Wuhan:Wuhan University of Science and Technology,2018.11Zaid ALDHAFEERI,Mamadou FALL.Time and damage induced changes in the chemical reactivity of cemented paste backfillJ.Journal of Environmental Chemical Engineering,2016(4):4038-4049.12 孙炳楠.工程中边界元法及其应用M.杭州:浙江大学出版社,1991.SUN Bingnan.Boundary Element Method and Its Application in EngineeringM.Hangzhou:Zhejiang University Press,1991.13 赵兴东,朱乾坤,赵一凡.板庙子金矿深部开采留设隔离矿柱控制地压数值优化J.采矿技术,2019,19(1):61-66.ZHAO Xingdong,ZHU Qiankun,ZHAO Yifan.Numerical optimiza-tion of controlling ground pressure by setting isolation pillars in deep mining of Banmiaozi Gold Mine J.Mining Technology,2019,19(1):61-66.14ROUX A J A,LEEMAN E R,DENKHAUS H G.Destressing:a means of ameliorating rockburst conditions.Part 1:the concept of destressing and results obtained from its applicationJ.Journal of South Africa Institute of Mining and Meta