刘 晓
(山东金洲矿业集团有限公司, 山东 威海 264500)
随着矿产资源的开发不断推进,浅部资源逐渐枯竭,深部开采已经成为矿山的主战场[1-2]。巷道担负着矿井运输、通风、行人等任务,其稳定性对矿井安全生产起到至关重要的作用[3-4]。随着矿井深部延伸,巷道受采动、矿山压力增大等影响开始变形或遭到破坏。爆破开挖是矿山生产的必要工序之一,其危害主要包括爆破振动、爆炸应力冲击及飞石等,诱导采场围岩及巷道失稳,降低矿山生产的安全稳定性。
一般主要进行巷道支护的静力学分析,但实际上,采空区围岩和巷道在承受静压的同时,还受到来自开采爆破的动力扰动[5-7]。姜明伟等[8]利用FLAC3D软件分析了动力扰动对顶板稳定性和断裂规律的影响,得到不同扰动强度下顶板破坏规律。唐礼忠等[9]利用Abaqus软件对爆破作用下巷道围岩的力学响应进行了数值分析,探讨了软弱夹层对巷道围岩稳定性的影响,提出了支护优化方案。朱万成等[10]采用数值分析方法,模拟了不同侧压条件下深部巷道失稳的动力触发过程,探索了动力扰动触发巷道岩爆的力学机制。邓红卫等[11]采用Midas/GTS与FLAC3D数值模拟技术,研究了不同扰动幅值下充填体的破坏规律及动力响应规律,发现动力扰动对充填体的振动特性有显著影响。研究动力扰动作用下巷道的失稳破坏过程及其振动特性,可以为矿山安全高效生产提供技术支持,但动力扰动作用下巷道的动态响应、失稳机制及其稳定性与时间相关的分析研究较少。
本文利用3DEC离散元软件对爆破动载荷作用下充填体内预留巷道的动态响应规律进行数值模拟分析,揭示动力扰动下巷道周围的变形特征以及应力场分布特征,为理论分析提供依据,提高矿山生产的安全性保障能力。
山东金洲矿业集团有限公司金青顶矿区位于山东省乳山市下初镇境内,坐落在黄垒河南岸南东庄南部的金青顶上,行政区划隶属乳山市下初镇管辖。矿区地理坐标:东经121°38′00″~121°41′15″,北纬37°04′00″~37°07′30″。矿区南距乳山市约25 km,距威海至桃村铁路下初站15 km,北距烟台港75 km,均有公路相接,交通方便,交通位置如图1所示。
图1 交通位置图
金青顶矿区地貌属低缓丘陵区,地面标高55.7~133.80 m,最高海拔133.8 m,地势较缓,地形坡度0°~15°,沟谷发育,植被广布。地表水主要为黄垒河,发源于昆嵛山南麓,自北向南从矿区西 2 km 处流过,向南注入黄海,当地历史最高洪水水位为43.8 mm。
成矿带内地层简单,构造及岩浆活动强烈,成矿地质条件优越。区内成矿前北北东、北东向断裂最为发育,成矿期北北东向断裂继续活动,成矿后以北西向断裂为主。主要控矿构造为北北东向断裂,纵贯全区,区域上有5条规模较大的断裂。这些断裂对金矿的形成和分布起控制作用。随着矿山开采深度的不断增加,矿体内外环境发生一系列变化,矿山压力逐渐增加,矿体赋存条件变得更加复杂,严重影响了采场的稳定性。
本文采用钢丝绳+钢丝网的形式进行充填体预留巷道,预留巷道的尺寸为3 m×3 m的三心拱,数值模拟采用的充填体高60 m,宽40 m。使用数值分析软件3DEC对爆破动荷载下充填体内预留巷道的稳定性进行模拟,3DEC软件是在二维软件UDEC基础上发展而来,内置几种常用的材料模型以及12种分析模式,可以对接触面的非连续力学行为进行模拟,除可模拟连续介质外,在处理非连续介质在外力作用下静、动态响应等问题上相对有限元软件有本质优势。在3DEC软件中可以将实体离散为黏结块体(Bonded Block Model,BBM)模型,具体如图2所示,其中,图2a为一柱状BBM模型,2b为其剖面示意图。
图2 BBM模型及块体之间的接触
采用灰砂比1∶8的充填体材料进行空区的充填,充填体及钢筋的力学参数见表1,其中1∶8充填体接触面的内摩擦角为36.94,黏聚力为0.171,法向刚度和切向刚度均为0.3 GPa/m。
表1 充填体及钢筋网力学参数
为了分析动荷载作用下充填体预留巷道的稳定性,分别在距离预留巷道中心轴线10 m、15 m、20 m以及25 m分别建立半径为10 cm、长度为1 m的炮孔模型。由于3DEC软件中需要对炮孔附近块体进行离散化处理,全模型采用BBM将极大地增加计算量,严重降低计算效率。因此,为了突出充填体损伤特征,仅在充填体预留巷道的重点研究区域,采用BBM模型进行充填体离散化处理,以获得充填体受动载作用下的裂隙扩展规律。在模型的BBM区域,块体遵守线性关系,接触面采用摩尔-库伦准则,模型其余部分采用连续体,其本构关系采用摩尔-库伦准则。
本文从巷道周围应力场分布特征以及破坏特征(塑性区)两方面对充填体巷道的动荷载稳定性做出全方位的定性定量分析评价。根据爆破动载理论,图3所示为简易爆破动载时程曲线,分别由上升阶段以及下降阶段组成,一般情况下,爆破动载的上升时间明显小于下降时间。根据工程实际,一般炸药爆破岩石产生的峰值应力取1.0 GPa,上升时间与总时间的比值为1∶5,爆破作用时间为600 us。
图3 数值仿真采用的爆破波理论曲线
本文仅考虑爆破距离为10 m的爆破动载荷对支护条件下预留巷道的稳定性影响,重点分析巷道周围的应力场分布特征以及巷道周围的破坏特征等。
图4所示为预留巷道周围最大主应力的分布云图,分析范围为0~1 MPa(-号表示压应力)。从图中可以看出,炮孔周围由于爆破作用,为主要的冲击受压区,在未受到动载作用之前,充填体预留巷道周围最大主应力表现为压应力,在应力波作用于巷道周围时,巷道上方植入的钢筋网以及悬吊钢丝绳等支护材料,能有效减小巷道顶板上方的受压范围,使钢筋网为主要受压承重体,在应力波影响后,巷道右侧上方及右侧底部位置的压应力区有所增加,但局限在1 MPa之内,表明爆破作用对充填体预留巷道产生的影响较小。
图4 距离10 m巷道周围最大主应力随时间变化
图5所示为预留巷道周围最小主应力的分布云图,分析范围为0~1 MPa(-号表示压应力,白色部位表示最小主应力为拉应力)。在未受到爆破作用前,BBM巷道模型中,除顶板钢筋网支护上端部位最小主应力表现为压应力,其余部位均变现为拉应力,随着爆破应力波的推进,爆破时间为30 ms时,应力波正好作用于巷道中心,上下顶板部位受到P波影响而表现为压应力,但也出现了靠近巷道表面的充填体表现为拉应力的现象,当应力波作用过后,时间为60 ms时,巷道顶板上方的压应力区相较爆破前有所减小,而巷道右侧下方压应力区有所增加。
图5 距离10 m巷道周围最小主应力随时间变化
以上说明,爆破作用在一定程度上改变了巷道周围的压应力、拉应力区的范围,对于巷道顶部而言,拉应力区有所增大,这对充填体的稳定性是不利的,另外,钢筋网+钢丝绳的支护方式对于减弱巷道周边由压应力区转变为拉应力区的范围起到了一定的作用。
1) 未支护(不加网)条件下巷道稳定性分析
图6、图7所示是在未支护条件下,爆破距离为10 m,爆破时间为30 ms、40 ms和50 ms时巷道周围裂隙的法向位移和切向位移。从图6和图7中可以看出,在爆破应力波传播过程中,巷道右侧帮首先出现较为明显的拉伸状态,而顶板右侧下部易出现较为明显的压缩状态,这种状态一直持续到爆破应力波传播之后,可以判定的是,巷道右侧帮容易出现拉伸破坏。底板部位由于应力波的作用,裂隙的剪切位移量也相对来说有所增大。
图6 未支护条件下不同爆破时间巷道周围裂隙法向位移图
图7 未支护条件下不同爆破时间巷道周围裂隙切向位移图
2) 支护条件下巷道稳定性分析
爆破距离为10 m时,预留巷道周围的充填体的裂隙张拉错动位移如图8、图9所示。从图8和图9中可以看出,在钢筋网及钢丝绳支护的上端,裂隙的剪切位移量增加较少,而在钢筋网及钢丝绳支护下端的巷道右上侧部位,裂隙的剪切位移量则明显增多。
图8 支护条件下不同爆破时间巷道周围裂隙法向位移图
图9 支护条件下不同爆破时间巷道周围裂隙切向位移图
对图6和图8支护前后预留巷道周围裂隙法向位移量进行对比分析。暂且将图6以及图8中30 ms 时方框作为钢筋网+钢丝绳柔性支护的影响范围区。从图中可以看出,在应力波的传播过程中,支护影响范围区的裂隙法向位移量要明显小于未支护条件下的裂隙法向位移量,且在未支护条件下,红色方框内的裂隙一般呈现出压缩变形,对应于负的裂隙法向位移量。但是,随着时间的推进,应力波的作用减弱,除了支护影响区域以外,巷道周围其余部位的裂隙法向位移量未出现较为明显的差异。
对图7和图9支护前后预留巷道周围裂隙切向位移量的进行对比分析。巷道周围裂隙切向位移量的空间分布特点与法向位移量的空间分布特点较为一致。在应力波传播的30 ms左右,巷道右侧上墙周边裂隙切向位移量的分布尤为明显,在未支护条件下,巷道直接受到应力波的强烈冲击影响。
这说明了采用钢筋网+钢丝绳的支护方式,对于削弱充填体预留巷道在爆破扰动下的剪切破坏,具有较强的效果。
采用3DEC软件对爆破动载荷作用下充填体内预留巷道的动态响应规律进行数值模拟分析,研究结果表明:
(1) 爆破作用在一定程度上改变了巷道周围的压应力、拉应力区的范围,对于巷道顶部而言,拉应力区有所增大,这对充填体的稳定性是不利的。
(2) 钢筋网+钢丝绳的支护方式对于减弱巷道周边由压应力区转变为拉应力区的范围起到了一定的作用。对于削弱充填体预留巷道在爆破扰动下的剪切破坏,具有较强的效果。
因此,为了降低爆破动载荷对充填体的影响,保持充填体稳定性,保证矿山生产的安全性,可以采用钢筋网+钢丝绳的预留巷道开采方式,合理有效地对巷道外爆破动力扰动进行很大程度上的削弱,这对于巷道稳定性的保持是非常有利的。