地下巷道对近区爆破动荷载的响应特性研究*

2021-09-15阎要锋王公忠严敏嘉柴修伟彭亚利金胜利

爆破 2021年3期

阎要锋，王公忠，严敏嘉,3，柴修伟，徐亮，彭亚利，金胜利

(1.武汉理工大学土木与建筑工程学院，武汉 430070；2.河南工程学院安全工程学院，郑州 450000；3.广西民族大学建筑工程学院，南宁 530006；4.武汉工程大学资源与安全工程学院，武汉 430073；5.湖北兴发化工集团有限责任公司，兴山 443700)

随着世界经济的高速发展，深部矿产资源开采、交通隧道、深埋存储、深地军事等工程将向深部岩体开挖转变，“高地应力”是深部岩石力学的重要研究方面，如锦屏二级水电站引水隧洞最大地应力可达到72 MPa；锦屏一级和官地等水电工程地下洞室的最大主应力也达到20～40 MPa。目前开挖施工中钻爆法依然是深部工程开挖最常用的方式，相比于浅层岩体开挖，在“三高一扰动”的影响下，深部岩体钻爆开挖会呈现出更多的工程技术难题。国内外许多学者对深部岩体工程进行了大量的研究探索，岳万友[1]、蔡济勇等采用有限元模拟方法研究深部裂隙岩体在爆破动荷载作用下的力学特性[2]；宋凯文等从理论上分析和数值模拟角度讨论了开挖过程中岩体内应力波传播规律以及加卸载应力波之间作用特征[3,4]，以及在显式分析里面地应力施加和初始地应力对掏槽爆破有明显的抑制作用[5,6]。Mandal和Singh发现在钻爆开挖过程中[7]，初始地应力的瞬态卸荷会大幅度增加隧洞的超挖量。严鹏等采用ANSYS的动力计算模块研究了巷道围岩的动态卸荷过程[8]，并将数值解与理论解进行了对比验证。卢文波等基于等效弹性边界概念[9]，提出了一种群孔起爆条件下爆炸荷载与开挖荷载瞬态卸荷耦合作用的等效数值模拟方法。

对于深部开挖卸荷下围岩力学特征，董春亮等认为主要表现为最小主应力卸荷、最大主应力集中[10]，而主应力差的瞬间增大，诱发开挖面内裂隙扩展、贯通，形成由强及弱连续分布的破坏区。Carter等针对长隧道的开挖研究发现[11]，初始地应力的瞬态卸载会产生强烈的拉应力，其幅值与卸荷速率有关。同时深部岩体往往存在高地应力，Toks和Kehrer认为高地应力诱导引起应变能快速释放诱发的地表振动甚至超过爆炸地震波引起的地表振动[12]。Cao等研究不同地应力状态下地应力瞬态卸荷诱发围岩振动的特性[13]。Yang等采用幅值谱分析和时频分析方法对地应力瞬态卸荷诱发的围岩振动进行了识别[14]。卢文波等认为深埋高地应力岩体爆破开挖，伴随着爆破破岩裂纹扩展以及新开挖面的形成[15]，开挖面上的地应力在爆破瞬间也随之释放。

国内外学者对深部岩体开挖过程进行了全面系统的研究，而对于逐步向深部开拓采掘的矿山，对于相邻开拓采掘巷道的爆破开挖施工的受力特征和变形破坏情况研究较少。因此通过构建数值计算模型，研究临近巷道爆破施工对于既有倒梯形隧道的影响，合理确定其安全距离有助于地下工程结构安全稳定，同时指导优化钻爆施工技术参数。

1 矿山地质背景及计算模型

1.1 地质背景

某磷矿地质构造复杂，矿层呈缓倾斜产出，为薄-中厚层，矿体形态基本呈层状分布且两矿层相距较近，形成缓倾斜叠层矿体。矿层顶板为中厚层状夹薄层状结构，其中夹有相对软弱薄层状云质泥岩。断面形状采用倒梯形布置，断面顶板宽为4.2 m，底板宽为4.0 m，边墙高为4.0 m，断面面积为16.4 m2，矿石坚固性系数f=10，巷道开挖断面特征如图1所示。

图 1 某巷道轮廓Fig. 1 The outline of a roadway

1.2 计算模型

本数值模型采用ANSYS/LSDYNA显式算法，为节省计算机时和提升计算效率，对该模型进行相应简化如下：(1)将三维模型抽象为平面应变模型来模拟在距离既有运输巷道不同距离处采准巷道爆破施工对既有运输巷道的影响。模型的尺寸为：X、Y、Z分别为5000 cm、2400 cm、6 cm；(2)考虑到建立群孔炸药的复杂性以及两个巷道距离相对较远，为便于数值计算，用集中装药取代施工现场的群孔装药；(3)在模型边界施加无反射边界条件实现以有限的模型模拟无限域；(4)针对不同工况下特征单元精确位置选取问题，考虑模型尺寸较大且相邻巷道距离较远特点，将岩石区域中的炸药位置不进行开挖设置炸药而是另行建立炸药单元附着在岩石上。模型材料采用2号岩石乳化炸药参数如表1所示，磷矿石的力学参数如表2所示，爆破数值模型如图2(a)所示，从图2(b)中单元时程压力曲线的特点可以看出，模型的简化处理是合理的。

表 1 炸药参数

表 2 岩石力学参数

图 2 模型概况及单元应力曲线图Fig. 2 Model overview and element stress curve

2 爆破数值计算模型结果分析

地下矿山、交通隧道、水利工程等在施工过程中，同一水平或者不同水平相邻巷道在不同时间甚至同一时间进行钻爆开挖施工时有发生，高能量的爆破荷载的输入严重威胁到巷道结构物的安全稳定，巷道结构物在施工过程中及其后续运行过程中安全稳定性成为亟待研究和解决的问题，因此以某个矿区某开采水平施工现场为背景进行针对性分析来阐明不同相对水平位置处巷道爆破施工对既有运输巷道的影响。

2.1 爆源在巷道右侧不同位置典型时刻应力云图分析

如图3所示，不同时刻下爆源位置在既有运输巷道右侧不同相对位置处的有效应力云图(依次记录500 μs、3000 μs、4200 μs、9000 μs处的四个典型时刻)，根据在应力波对岩石的做功过程中，根据炸药对岩石做功的理论和原理，岩石在炸药附近最先出现环向裂纹，然后出现径向裂纹，裂纹交叉贯通实现了岩石破碎，在爆生气体巨大的准静态气体压力作用下实现碎石的抛掷飞溅。图3中的4个模型显示爆源位置距离既有运输巷道依次是10 m、14 m、18 m、22 m，四个计算模型的有效应力云图形态的变化过程大致经历以下三个阶段：第一阶段在均质各向同性的物理环境下应力波呈现以炸药为中心以柱面波的形式向四周传播特点，其具体形态见500 μs时刻应力云图且随着时间继续应力范围不断增大。第二阶段：随着应力波继续传播，呈显现出如图3(a)中第2幅图片所示的应力集中现象，当应力波在行进中遇到结构物，应力云图产生剧烈的“结构效应”，见图3(a)第3幅图，应力在巷道右侧壁及其上隅角出现剧烈的应力集中，同时应力波在巷道侧壁处产生强烈的反射，行进的应力波和反射的应力波在削弱应力的同时改变应力波的而前进方向。第三阶段：当应力波波阵面通过巷道结构后，应力波的形态发生显著变化，波阵面由最初凸形状变为凹形状且在最凹处应力强度最低，在巷道左侧壁处显示为应力降低区，仅在巷道左侧上下隅角出现微弱的应力集中，随着爆源和既有运输巷道距离加大，到达既有运输巷道轮廓周围的应力波成分变得异常复杂，如图3(a)的第2、3幅图所示，出现这样情况这可能与巷道结构形状、应力波的叠加相关也可能与应力波传播一定距离后因强度衰减频率分散有关，限于篇幅，仅对爆源位置在巷道右侧10 m进行具体研究和讨论，后续的14 m、18 m、22 m的工况产生的有效应力波云图均出现上述特征，这里不再赘述。

图 3 爆源位置在既有运输巷道右侧典型时刻应力云图Fig. 3 Stress nephogram of typical moments at the right side of the existing transportation roadway where the detonation source is located

2.2 爆源位置在巷道右侧不同位置处典型时刻单元应力分析

在2.1节中采取定性的分析方法通过用应力云图分析来衡量和评价距离不同距离处的临近巷道爆破产生的应力波对既有运输巷道结构安全和稳定性的分析，为了能够从不同方位和不同损伤破坏深度研究临近巷道爆破对既有运输巷道的影响，从单元拉、压应力曲线角度对模型在顶板、底板、巷道帮部三个特殊部位选取特征单元进行研究，图4所示的A、B、C、D表示临近巷道距离既有运输巷道依次是10 m、14 m、18 m、22 m，图中横轴为单元编号简写形式。

图4显示巷帮附近距离洞壁一定距离处选取的13个特征单元的应力状态，每个单元的尺寸为6cm，由表2参数可知岩石单元的应力状态均小于岩石的抗压强度，且爆源距离洞壁距越远其压应力越小，并且都在距离洞壁附近处其压力值变化幅度较大，在极小范围内单元由最小压应力变化为最大压应力，其最大压应力低于岩石的抗压强度；特征单元拉应力曲线亦在距离隧道洞壁42 cm处变化明显，尤其爆源在距离既有运输巷道10 m，14 m处拉应力变化最为明显，且拉应力数值较大超过岩石的抗拉强度。由上可知，在巷道巷帮附近处岩石的破坏以拉应力破坏为主，同时相邻巷道爆破施工巷帮的安全范围为14～18 m。

图 4 运输巷道巷帮附近单元应力Fig. 4 Unit stress near roadway side of transportation roadway

图(5)表明，巷道底板附近特征单元应力变化特征，对于单元压应力：随着爆源距离的增加，爆源在距离既有运输巷道10～18 m阶段对应单元的应力衰减梯度较为明显，爆源在距离既有运输巷道18～22 m阶段单元应力衰减幅度降低。各个模型中距离巷道轮廓越远压应力数值越小，这反映出与一般模型规律一致的特点，同时由图5和表2可知单元的压应力远远小于岩石的抗压强度不能引起岩石的直接破坏。对于单元拉应力：其变化特征基本遵循拉应力的变化规律，但爆源在距离既有运输巷道10～14 m阶段时，单元拉应力变化幅度较大，这可能与爆炸应力波作用特点有关，也与巷道结构特点有关，最大拉应力超过岩石的抗拉强度，爆源在距离既有运输巷道18～22 m阶段单元拉应力变化较为平缓，且拉应力数值较小，其数量级不足以造成岩石的破坏。因此，对于相邻巷道爆破施工的巷道底板来说巷道变形破坏仍然是以拉应力为主。

图 5 运输巷道底板附近单元应力Fig. 5 Unit stress near the floor of transport roadway

图(6)显示巷道顶板附近的单元应力特征，单元编号从1到13表示特征单元距离隧道轮廓面越来越近，对于单元压应力：在爆源距离既有运输巷道为10 m时，无论压应力还是拉应力特征单元曲线均出现震荡的特点，这与应力波的作用特点以及和巷道开挖后巷道结构特点有关；同样是在爆源距离既有运输巷道为18、22 m时，特征单元的应力曲线变化趋于平稳，巷道周边的岩体仍是以拉应力破坏为主，且既有运输巷道和临近巷道的安全范围在14～18 m之间。