超埋深高应力复杂地层巷道围岩损伤分区与承载结构特征研究.pdf

Series No.566August 2023 金 属 矿 山METAL MINE 总 第566 期2023 年第 8 期收稿日期 2023-02-28基金项目 中铝集团科技重大专项项目(编号:ZLZX20210102)。作者简介 雷显权(1985),男,正高级工程师,博士后。超埋深高应力复杂地层巷道围岩损伤分区与承载结构特征研究雷显权1,2,3 刘福春1,2 严庆文4 张帅帅1,2 郭奇峰5 王勇兵4 汪小东1,2(1.长沙有色冶金设计研究院有限公司,湖南 长沙 410019;2.深井矿山安全高效开采技术湖南省工程研究中心,湖南 长沙 410019;3.中南大学资源与安全工程学院,湖南 长沙 410083;4.云南驰宏锌锗股份有限公司会泽矿业分公司,云南 曲靖 654200;5.北京科技大学土木与资源工程学院,北京 100083)摘 要 巷道围岩损伤分区与承载结构是巷道岩层控制与支护设计的重要依据。以云南某深井开采矿山 1 500 m 埋深巷道为例,采用弹塑性理论分析、松动圈测试和数值模拟相结合的方法,研究了不同地层巷道围岩的损伤分区和承载结构特征,分析了不同围岩承载结构的差异及其与巷道稳定性的关系。结果显示:巷道围岩松动圈深度范围为 0.91.9 m,平均为 1.38 m;灯影组、宰格组、大塘组和摆佐组松动圈平均深度分别为 1.7 m、1.15 m、1.17 m 和1.5 m。数值模拟揭示巷道围岩破碎区深度范围为 0.921.63 m,平均为 1.20 m;塑性区深度范围为 1.482.50 m,平均为 1.87 m;不同地层围岩损伤区深度由大到小排序为灯影组栖霞茅口组宰格组大塘组摆佐组。不同地层围岩承载区厚度差别不大,但承载区深度存在明显差异;关键承载区厚度范围为 1.091.34 m,平均为 1.28 m;关键承载区的内边界深度范围为 1.222.18 m,平均为 1.58 m;不同地层围岩关键承载区深度由大到小排序为灯影组栖霞茅口组宰格组大塘组摆佐组。研究表明:在巷道规格及所处环境相同条件下,岩体强度是影响围岩承载结构和巷道稳定性的最主要因素;岩体强度越大,承载区深度越小,损伤区深度越浅,围岩变形量越小,巷道越稳定;对于锚杆支护,锚杆有效长度等于围岩关键承载区的内边界深度是比较合理的取值。关键词 深井开采 巷道支护 围岩控制 承载结构 松动圈 高应力 中图分类号TD353 文献标志码A 文章编号1001-1250(2023)-08-171-10DOI 10.19614/ki.jsks.202308021Research on Damage Zone and Load-bearing Structure Characteristics of Surrounding Rock Mass of Roadway in Complex Strata with Overburied Depth and High in Situ StressLEI Xianquan1,2,3 LIU Fuchun1,2 YAN Qingwen4 ZHANG Shuaishuai1,2 GUO Qifeng5WANG Yongbing4 WANG Xiaodong1,2(1.CINF Engineering Co.,Ltd.,Changsha 410019,China;2.Hunan Provincial Engineering Research Center for Safe and Efficient Mining Technology of Deep Mines,Changsha 410019,China;3.School of Resources and Safety Engineering,Central South University,Changsha 410083,China;4.Huize Mining Branch,Yunnan Chihong Zn&Ge Co.,Ltd.,Qujing 654200,China;5.School of Civil and Resources Engineering,University of Science and Technology Beijing,Beijing 100083,China)Abstract The damage zone and load-bearing structure of roadway surrounding rock are the important basis for the road-way rock control and support design.Taking 1 500 m buried roadway of a deep mining mine in Yunnan province as an exam-ple,through elastoplastic theory analysis,rock loose circle test and numerical simulation,the damage zone and load-bearing structure characteristics of surrounding rock of different strata have been studied,and the difference of the load-bearing struc-ture of different surrounding rock and its relationship with the stability of roadway have been analyzed.The results show that:The depth of the rock loose circle ranges from 0.9 to 1.9 m,with an average of 1.38 m;and the average depth of the rock loose circle in Dengying Formation,Zaige Formation,Datang Formation and Baizuo Formation is 1.7 m,1.15 m,1.17 m and 1.5 m,respectively.The numerical simulation reveals that the depth of the broken zone of surrounding rock ranges from 0.92 to 1.63 m,with an average of 1.20 m;and the depth of the plastic zone of surrounding rock ranges from 1.48 to 2.50 m,with an average of 1.87 m;the depth of damage zone in different strata is ranked as Dengying Formation Qixia-Maokou For-171mationZaige FormationDatang FormationBaizuo Formation.The difference in the thickness of the load-bearing zone of different strata is little,but the difference in the depth of the load-bearing zone is obvious;the thickness of the critical load-bearing zone ranges from 1.09 to 1.34 m,with an average of 1.28 m;the depth of the inner boundary of the critical load-bear-ing zone ranges from 1.22 m to 2.18 m,with an average of 1.58 m;the depth of the critical load-bearing zone in different stra-ta is ranked as Dengying FormationQixia-Maokou FormationZaige FormationDatang FormationBaizuo Formation.The re-search shows that:Under the same conditions of roadway size and environment,the strength of rock mass is the most important factor affecting the load-bearing structure of surrounding rock and the stability of roadway.The greater the rock strength,the smaller the load-bearing zone depth,the shallower the damage zone depth,the smaller the surrounding rock deformation,the more stable the roadway.For rockbolting support,the reasonable effective length of rockbolt is equal to the inner boundary depth of the critical bearing zone of surrounding rock.Keywords deep mining,roadway support,surrounding rock control,load-bearing structure,rock loose circle,high in situ stress 随着浅部矿产资源的日益枯竭,我国矿产资源开发即将全面进入第二深度空间(1 0002 000 m),深部资源开采逐渐成为矿产资源开发新常态1-3。进入深部开采后,地应力显著增大,地质环境复杂多变,矿岩非线性力学特征显著,开采扰动强烈,巷道变形、片帮、塌方、冒顶、岩爆等灾害问题突出,给深部巷道围岩稳定性控制与矿山生产安全防控造成巨大挑战4-6。为了应对深部高应力复杂环境条件下巷道失稳致灾问题,相关学者开展了大量研究,在深部开采岩体力学特性、围岩变形破坏机理、岩层稳定性控制方法、巷道支护技术和地压灾害监测预警等方面取得了长足进展7-12。巷道围岩的损伤分区和承载结构特征是巷道稳定性评价的重要内容,也是巷道岩层控制与支护设计的重要依据。康红普13-14于1997 年首次提出巷道围岩关键承载圈概念,分析了围岩承载圈的分布特征及其影响因素,阐述了关键承载圈对巷道支护的工程意义。李树清等15-16采用弹塑性理论分析和数值模拟方法研究了巷道围岩的承载结构特征,分析了深部巷道与浅部巷道围岩承载结构的差别,探讨了支护阻力对深部巷道围岩承载结构的影响。赵光明等17-18通过理论分析和数值模拟,将软弱破碎巷道围岩的承载结构划分为浅支撑层、深支撑层和关键支撑层,分析了各支撑层对围岩变形的控制作用以及岩体残余强度和支护作用对围岩承载结构演化和巷道变形的影响。Qin 等19采用理论分析和数值模拟相结合的方法研究了动压作用下深部软岩巷道围岩的承载结构,分析了不同支护方式对围岩承载结构的控制效果。Zhang 等20采用理论分析、数值模拟和现场试验相结合的方法研究了巷道开挖空间效应产生的应力调整和重分布对围岩力学特性及承载结构的影响,揭示了岩石破裂发展状态对围岩承载结构的调节机制。王玺等21通过理论计算和数值模拟实验,研究了不同岩体质量等级巷道围岩的承载结构特征,分析了岩体质量等级对各承载区主应力分布和塑性破坏的影响。这些研究深化了对巷道围岩承载结构特征及其作用机理的认识,对巷道围岩稳定性控制与支护具有重要指导作用。充分发挥围岩的自身强度和自承能力以达到采矿工程的稳定是矿山岩石力学最基本的原理22-23。为了实现这一目的,首先就要了解和掌握围岩的损伤分区和承载结构特征。虽然前人对巷道围岩承载结构相关问题进行了大量研究并取得了丰富成果,但由于深井矿山开采环境的复杂性,还存在许多问题需要作进一步研究和探讨,比如超埋深高应力复杂地层巷道围岩承载结构的差异性及稳定性问题。针对该问题,本研究以云南某深井开采矿山 1 500 m 埋深中段石门巷道为例,采用弹塑性理论分析、松动圈现场测试和数值模拟计算相结合的方法,研究了超埋深高应力条件下不同地层巷道围岩的损伤分区和承载结构特征,分析了不同地层围岩损伤区和承载结构的差异及其原因,探讨了不同地层围岩承载结构与巷道稳定性的关系,提出了一种确定深部巷道锚杆支护长度的新方法,为该矿山深部巷道围岩稳定性控制与支护提供理论指导。1 工程概况云南某矿山是目前国内已建成投产的开采深度最深的矿山,最低开采中段埋深超过 1 500 m。矿山地层、构造条件复杂,深部巷道埋深大、地应力高,开采扰动强烈,具有典型的“三高一扰动”特征。现场调查结果显示,在深部复杂环境条件下,巷道围岩变形、喷层开裂和支护破坏现象多见,发生片帮、冒顶和大规模失稳的安全风险突出,给巷道稳定性控制造成极大挑战。以该矿山深部某中段石门巷道为研究对象(图1)。该中段埋深约 1 500 m,巷道断面形状为三心拱271总第 566 期 金 属 矿 山 2023 年第 8 期形,断面规格为 3.40 m3.20 m。巷道从竖井马头门往 SE 方向开拓掘进,依次穿越震旦系陡山沱组、灯影组、寒武系筇竹寺组,泥盆系宰格组,石炭系大塘组、摆佐组、威宁组、马坪组,二叠系梁山组、栖霞茅口组等多组地层以及多条断层。图 1 研究区巷道布置Fig.1 Roadway layout in research area 根据巷道埋深和上覆岩层重度估算垂向应力约为 42 MPa。矿山 1 500 m 深处的地应力测量结果显示,最大主应力为 42.34 45.95 MPa,平均 43.78 MPa,最大主应力方向为 NNW-SSE24。现场节理裂隙调查和岩体质量评价结果显示,巷道围岩的岩体质量等级以 级为主。通过室内基础岩石力学试验,获得了不同地层岩块试样的力学参数。以此为基础,结合岩体节理裂隙调查和岩体质量评价结果,采用等效 Hoek-Brown 强度准则与 Mohr-Coulomb 强度准则的岩体力学参数估算方法25-26,获得了深部巷道围岩的岩体力学参数,见表 1。2 深部巷道围岩损伤分区与承载结构模型假设巷道围岩为均质各向同性连续介质,本构关系为弹塑性软化模型。基于圆形开挖断面和小孔扩张理论,建立如图 2 所示的巷道围岩受力分析模型。巷道开挖后,围岩应力发生二次重分布,并产生应力集中,导致围岩发生不同程度的损伤破裂。根据围岩应力分布特征和损伤破裂程度,围岩变形损伤区由表及里依次划分为破碎区(或松动区)、塑性软化区、弹性区和原岩应力区8,17-18。其中,破碎区和塑性软化区为塑性变形区,即广义上的塑性区;弹性区和原岩应力区为弹性变形区。表 1 深部岩体力学参数Table 1 Mechanical parameters of deep rock mass地层岩性岩体质量等级密度/(kg/m3)弹性模量 E/GPa泊松比 黏聚力 c/MPa内摩擦角/()抗压强度 c/MPa抗拉强度 t/MPa灯影组白云岩2 8164.380.285.4920.358.720.88宰格组白云岩2 81010.150.266.5823.5212.261.92大塘组灰岩夹泥层2 74811.430.266.8124.0511.602.34摆佐组白云岩2 72516.110.267.1725.0213.232.65栖霞茅口组灰岩2 8326.610.266.0722.2512.311.78图 2 巷道围岩受力分析模型Fig.2 Mechanical model of roadway surrounding rockP0原岩应力;Pi支护阻力;R0巷道半径;Rb破碎区半径;Rp塑性(软化)区半径根据巷道开挖后围岩二次应力分布规律,并结合围岩变形损伤分区特征,建立如图 3 所示的深部巷道围岩承载结构模型。将 P0的区域定义为承载区。如图 3(a)所示,破碎区(A-B)为应力降低区(kP0),是承受围岩压力的关键区域,将其划为关键承载区(E-F)。研究表明,k 值范围一般为 1.21.517-18,21。综合考虑计算得到的巷道围岩二次应力分布特征及其与最大峰值应力和原岩应力的关系,本研究 k 值取1.3。图 3(b)显示了深部巷道围岩承载区与变形损伤区的对应关系。3 深部巷道围岩损伤分区与承载结构特征3.1 弹塑性理论计算塑性区深度基于圆形巷道力学模型和 Mohr-Coulomb 屈服准则,通过弹塑性力学理论分析,可求得巷道围岩塑性371 雷显权等:超埋深高应力复杂地层巷道围岩损伤分区与承载结构特征研究 2023 年第 8 期图 3 深部巷道围岩承载结构模型Fig.3 Load-bearing structure model of deep roadway surrounding rockr径向应力;环向应力;max弹塑性边界上的最大峰值应力区半径 Rp的计算公式,即卡斯特奈(H.Kastner)公式或修正芬纳(R.Fenner)公式27:Rp=R0P0+ccotPi+ccot1-sin()1-sin2sin.(1)原岩应力可通过地应力实测获得,或采用下式进行估算:P0=H,(2)式中,为岩体容重;H 为巷道埋深。当支护阻力 Pi为 0,即巷道开挖未支护时,围岩塑性区半径 Rp计算公式为Rp=R0P0+ccotccot1-sin()1-sin2sin.(3)为了便于计算,基于最小外接圆原则,采用几何作图法得到矿山深部三心拱形巷道的等价圆形断面,以获得巷道等效半径,即 2.08 m。然后根据研究矿山深部巷道围岩的岩体力学参数(表 1)及其所处的地应力环境(原岩应力 P0取 42 MPa),计算得到不同地层巷道围岩的塑性区深度,见表 2。计算结果显示,不同地层围岩塑性区深度范围为 1.623.03 m,平均 2.16 m。3.2 声波法现场测试松动圈深度根据巷道围岩松动圈支护理论28-29,巷道开挖后,原岩应力在围岩中发生重分布并产生应力集中,当集中应力超过围岩强度时,巷道周边围岩发生破坏,形成一条环向破裂带,该破裂带被称为围岩松动表 2 不同地层围岩塑性区理论计算结果Table 2 Theoretical calculation results of plastic zone of surrounding rock in different strata地层巷道等效半径/m 塑性区半径/m塑性区深度/m灯影组2.084.923.03宰格组2.083.861.97大塘组2.083.721.83摆佐组2.083.511.62栖霞茅口组2.084.242.35平均2.084.052.16圈。松动圈的范围通常采用现场实测方法获得,作为巷道支护分类和参数设计的依据。本研究采用声波法进行松动圈测试。声波在岩体中的传播速度与岩体的结构及其物理力学性质密切相关。在裂隙发育、低密度、高声阻抗岩体中,声波传播速度较慢;在高应力、高密度岩体中,声波传播速度较快。因此,通过测试超声波在巷道围岩钻孔内的传播速度,绘制波速孔深曲线,分析波速沿孔深的变化规律,就可以判断围岩松动圈的范围。现场测试采用武汉中科智创生产的 RSM-SY5 型智能声波检测仪以及与之配套的一发双收探头。该检测仪具备自动判读功能,使得首波声时、声幅读取更为准确。测点布置在巷道侧壁,针对不同地层巷道围岩,共完成了 11 个点的松动圈测试,结果见图 4 和表 3。实测结果显示,不同地层围岩松动圈深度范围为 0.91.9 m,平均 1.38 m。图 4 松动圈测试波速孔深曲线Fig.4 Wave velocity-hole depth curve of rock loose circle test3.3 数值模拟分析塑性区特征根据研究矿山深部巷道几何特征,采用 FLAC3D软471总第 566 期 金 属 矿 山 2023 年第 8 期表 3 不同地层围岩松动圈测试结果Table 3 Test results of rock loose circle in different strata地层测点/个松动圈深度范围/m平均/m灯影组31.51.91.7宰格组41.11.31.15大塘组30.91.51.17摆佐组11.51.5栖霞茅口组0合计/平均110.91.91.38件,建立巷道开挖数值计算模型。采用 8 节点六面体单元划分网格,为了提高计算精度,对巷道周边网格进行了加密,划分后的网格模型包括 46 820 个单元和 52 030 个节点。模型几何尺寸和网格划分情况如图 5 所示。图 5 巷道模型网格Fig.5 Mesh grid of roadway model 巷道围岩的变形破坏行为采用耦合张拉屈服准则的 Mohr-Coulomb 强度准则来描述,屈服函数如下:Fs=1-1+sin1-sin3+2c1+sin1-sinFt=3-t,(4)式中,1和 3分别为最大和最小主应力;Fs和 Ft分别为剪切和张拉屈服函数。模型物理力学参数见表 1。根据地应力实测和理论估算结果,并考虑巷道走向与最大主应力方向近似平行的关系,在模型 x、y 方向上分别施加 42 MPa和 43.78 MPa 的水平应力,在模型 z 方向上施加 42 MPa 的垂向应力和 9.8 m/s2的重力加速度。采用 FLAC3D软件计算模拟巷道开挖过程。根据数值模拟结果,绘制了不同地层巷道开挖后的围岩塑性区分布结果,见图 6。在此基础上,根据上述深部巷道围岩损伤分区理论模型,获得了不同地层巷道围岩的损伤区(破碎区和塑性区)深度,见表 4。数值模拟结果显示,不同地层围岩破碎区深度范围为0.921.63 m,平均 1.20 m;不同地层围岩塑性区深度范围为 1.482.50 m,平均 1.87 m。3.4 深部巷道围岩承载结构划分基于数值模拟得到的巷道围岩应力分布规律和损伤区分布结果,根据上述深部巷道围岩承载结构理论模型,划分了不同地层巷道围岩的承载区,见图 7,图 6 不同地层围岩塑性区分布特征Fig.6 Distribution characteristics of plastic zone of surrounding rock in different strata并获得了不同承载区的结构参数,见表 5。研究结果显示,巷道围岩塑性、弹性和关键承载区平均深度范围分别为 1.201.88 m、1.888.09 m 和 1.582.86 m,平均厚度分别为 0.68 m、6.21 m 和 1.28 m。571 雷显权等:超埋深高应力复杂地层巷道围岩损伤分区与承载结构特征研究 2023 年第 8 期表 4 不同地层围岩损伤区数值模拟结果Table 4 Numerical simulation results of damage zone of surrounding rock in different strata地层破碎区深度/m塑性区半径/m塑性区深度/m灯影组1.634.202.50宰格组1.103.421.72大塘组1.033.341.64摆佐组0.923.181.48栖霞茅口组1.303.732.03平均1.203.571.874 分析与讨论4.1 不同围岩损伤区与承载结构的差异4.1.1 不同地层围岩损伤区深度对比图 8 所示为采用上述 3 种方法得到的不同地层巷道围岩损伤区深度对比。从图 8 可以看出,不同地层围岩的损伤区深度存在明显差异,塑性区及破碎区深度由大到小排序为:灯影组栖霞茅口组宰格组大塘组摆佐组。不同方法得到的损伤区深度虽有差异但变化趋势基本一致。由图 8(a)可以看出,图 7 不同地层巷道围岩最大主应力分布特征Fig.7 Distribution characteristics of maximum principal stress in surrounding rock of different strata表 5 不同地层围岩承载区划分结果Table 5 Results of load-bearing zone division of surrounding rock mass in different strata地层塑性承载区深度范围/m厚度/m弹性承载区深度范围/m厚度/m关键承载区深度范围/m厚度/m灯影组1.632.50 0.87 2.508.87 6.37 2.183.27 1.09宰格组1.101.72 0.62 1.727.92 6.20 1.452.79 1.34大塘组1.031.64 0.61 1.647.81 6.17 1.362.69 1.33摆佐组0.921.48 0.56 1.487.61 6.13 1.222.54 1.32栖霞茅口组1.302.03 0.73 2.038.22 6.19 1.713.05 1.34平均1.201.88 0.68 1.888.09 6.21 1.582.86 1.28对于同一地层,围岩塑性区深度理论计算结果比数值模拟结果偏大,主要原因是理论计算没有考虑巷道围岩的三维结构。由图 8(b)可以看出,对于同一地层围岩,灯影组、宰格组和大塘组的破碎区深度数值模拟结果略小于松动圈实测平均结果,但都在松动圈实测结果(表 3)范围内;摆佐组破碎区深度数值模拟结果与松动圈实测结果相差较大,主要是因为摆佐组只进行了 1 次松动圈实测,没有反映松动圈深度的统计特征。上述 3 种方法研究结果的对比分析也表明,数值模拟方法能够较为客观准确地揭示深部巷道围岩的损伤破坏特征。图 9 所示为巷道围岩损伤区深度与岩体强度的关系。由图 9 可以看出,围岩塑性区、破碎区深度均与岩体抗压、抗拉强度呈负线性相关关系。总体而言,岩体抗压、抗拉强度越大,围岩损伤区深度越小。这表明,在相同的地应力条件下,巷道围岩的损伤破坏程度主要取决于岩体强度。虽然这一结论早已从理论上得到预测,但本研究进一步基于深井矿山高应力巷道实例并通过围岩松动圈现场实测和三维数值模拟实验验证了这个结论。4.1.2 不同地层围岩承载结构对比图 10 所示为不同地层巷道围岩承载结构参数对比。由图 10(a)可以看出,不同地层围岩的承载区厚度没有明显差异,除灯影组外,宰格组、大塘组、摆佐组和栖霞茅口组围岩的塑性、弹性和关键承载区厚671总第 566 期 金 属 矿 山 2023 年第 8 期图 8 不同地层围岩损伤区对比Fig.8 Comparision of damage zone of surrounding rock in different strata图 9 围岩损伤区与岩体强度关系Fig.9 Relationship between damage zone and strength of surrounding rock度均非常接近。由图 10(b)可以看出,不同地层围岩的承载区深度存在显著差异,但塑性、弹性和关键承载区深度变化趋势一致,由大到小排序为:灯影组栖霞茅口组宰格组大塘组摆佐组。这也与上述不同地层围岩损伤区深度排序相一致。图 11 所示为巷道围岩承载结构参数与岩体强度的关系。从图 11 可以看出,围岩塑性、弹性和关键承图 10 不同地层围岩承载结构对比Fig.10 Comparison of load-bearing structures of surrounding rock in different strata载区深度均与岩体抗压、抗拉强度呈负线性相关关系。总体而言,岩体抗压、抗拉强度越大,围岩承载区深度越小。图 11 围岩承载结构与岩体强度关系Fig.11 Relationship between load-bearing structures and strength of surrounding rock mass4.2 围岩承载结构与巷道稳定性的关系图 12 所示为不同地层巷道开挖后围岩变形特征对比。从图 12 可以看出,不同地层围岩的总体变形特征相似,但位移量存在明显差异,顶板、巷帮最大位771 雷显权等:超埋深高应力复杂地层巷道围岩损伤分区与承载结构特征研究 2023 年第 8 期移由大到小排序为灯影组栖霞茅口组宰格组大塘组摆佐组。图 12 不同地层围岩变形特征对比Fig.12 Comparison of deformation characteristics of surrounding rock in different strata 图 13 所示为巷道开挖后围岩位移、损伤区深度与关键承载区深度的关系。从图 13 可以看出,顶板、巷帮位移量以及围岩塑性区、破碎区深度均与关键承载区深度呈正线性相关关系。总体而言,关键承载区深度越大,围岩位移量和损伤区深度越大。图 13 围岩变形损伤与承载结构关系Fig.13 Relationship between deformation damage and load-bearing structure of surrounding rock4.3 对深部巷道围岩稳定性控制的启示了解和掌握巷道围岩的损伤分区和承载结构特征对巷道稳定性控制具有重要意义。研究结果显示,在深部高应力条件下,不同地层围岩的损伤区和承载区结构具有显著差异,而且关键承载区深度对围岩损伤变形具有明显控制作用。这就要求在制定深井矿山高应力巷道支护方案时,需要根据各类围岩的实际损伤破坏特征,有针对性地选择恰当的支护方法和确定合理的支护参数。以目前矿山常用的锚杆支护为例,锚杆长度是锚杆支护最重要的技术参数,但如何确定锚杆长度是锚杆支护设计的难点。通常做法是将锚杆端头锚固到围岩深部稳定岩层,或取锚杆有效长度大于塑性区深度,但研究和实践表明这 2 种方法确定的锚杆长度都偏安全,容易造成支护浪费。另一种工程实践性很强的方法是根据围岩松动圈深度来确定锚杆长度,但松动圈测试现场实施难度大,而且很多测量结果并没有明显的松动圈界线。根据承压拱(或压力拱)支护理论30-31,锚杆支护的主要作用是在围岩损伤区内重建人造承压拱,提高损伤围岩的整体强度和自承能力,以此来抵抗围岩变形压力和维护巷道稳定。这也符合深部巷道围岩稳定性控制的基本原理23。结合承压拱支护理论和本文研究结果,综合分析认为较为合理的做法将锚杆锚固到围岩关键承载区内,即取锚杆有效长度等于围岩关键承载区的内边界深度(图871总第 566 期 金 属 矿 山 2023 年第 8 期14)。这样既避免了锚杆过长浪费,也避免了因锚杆支护长度不足所导致的围岩控制不稳。图 14 深部巷道锚杆支护示意Fig.14 Schematic of rockbolt support in deep roadway5 结 论针对深井矿山高应力巷道围岩稳定性控制问题,以云南某深井开采矿山 1 500 m 埋深中段石门巷道为例,采用弹塑性理论分析、松动圈现场测试和数值模拟计算 3 种方法,研究了深井高应力复杂地层巷道围岩的损伤分区和承载结构特征及其对巷道稳定性控制的意义。(1)根据巷道开挖后围岩损伤分区特征和应力分布规律,建立了深部巷道围岩承载结构理论模型,将深部巷道围岩承载区划分为塑性承载区、弹性承载区和关键承载区 3 个部分,提出了确定各承载区范围的方法。(2)围岩松动圈现场实测结果显示,该矿山深部巷道围岩松动圈深度总体变化范围 0.91.9 m,平均1.38 m。其中:灯影组松动圈深度范围 1.51.9 m,平均 1.7 m;宰格组松动圈深度范围 1.11.3 m,平均 1.15 m;大塘组松动圈深度范围 0.91.5 m,平均1.17 m;摆佐组松动圈深度 1.5 m。(3)通过 3 种不同方法的综合研究和对比分析,揭示了该矿山深部巷道围岩的损伤区发育特征。数值模拟结果显示不同地层围岩的破碎区深度范围为0.921.63 m,平均 1.20 m;塑性区深度范围为 1.482.50 m,平均 1.87 m。虽然所处的地应力环境相同,但不同地层围岩的损伤区深度存在明显差异,由大到小排序为:灯影组栖霞茅口组宰格组大塘组摆佐组。(4)基于深部巷道围岩承载结构理论模型和数值模拟计算结果,划分了该矿山深部巷道围岩的承载结构。不同地层围岩的承载区厚度差别不大,但承载区深度存在显著差异。关键承载区厚度范围为 1.091.34 m,平均 1.28 m;关键承载区内边界深度范围为 1.222.18 m,平均 1.58 m;不同地层关键承载区深度由大到小排序为:灯影组栖霞茅口组宰格组大塘组摆佐组。(5)在巷道规格和地应力环境相同条件下,岩体强度是影响围岩承载结构和巷道稳定性的最主要因素。总体而言,岩体强度越大,关键承载区深度越小,围岩损伤区深度越浅,围岩变形量越小,巷道越稳定。(6)巷道围岩损伤分区和承载结构研究结果为矿山深部高应力巷道稳定性控制与支护参数选取提供了理论依据。针对锚杆支护,提出锚杆有效长度等于巷道围岩关键承载区的内边界深度是比较合理的取值。参 考 文 献1 李夕兵,周健,王少锋,等.深部固体资源开采评述与探索J.中国有色金属学报,2017,27(6):1236-1262.LI Xibing,ZHOU Jian,WANG Shaofeng,et al.Review and practice of deep mining for solid mineral resourcesJ.The Chinese Journal of Nonferrous Metals,2017,27(6):1236-1262.2 谢和平.深部岩体力学与开采理论研究进展J.煤炭学报,2019,44(5):1283-1305.XIE Heping.Research review of the state key research development program of china:deep rock mechanics and mining theoryJ.Jour-nal of China Coal Society,2019,44(5):1283-1305.3 蔡美峰,谭文辉,吴星辉,等.金属矿山深部智能开采现状及其发展策略J.中国有色金属学报,2021,31(11):3409-3421.CAI Meifeng,TAN Wenhui,WU Xinghui,et al.Current situation and development strategy of deep intelligent mining in metal minesJ.The Chinese Journal of Nonferrous Metals,2021,31(11):3409-3421.4 RANJITH P G,ZHAO J,JU M H,et al.Opportunities and challenges in deep mining:a brief reviewJ.Engineering,2017,3(4):527-537.5 FAIRHURST C.Some challenges of deep miningJ.Engineering,2017,3(4):546-551.6 赵兴东,周鑫,赵一凡,等.深部金属矿采动灾害防控研究现状与进展J.中南大学学报(自然科学版),2021,52(8):2522-2538.ZHAO Xingdong,ZHOU Xin,ZHAO Yifan,et al.Resea