岩体受到开挖卸荷等扰动时，伴随着裂缝的产生及扩展，储存的弹性能以应力波的形式释放，从而产生微震事件[1]。微震信号包含了有关岩体变形破坏的大量信息，微震监测通过微震系统布置、微震信号采集、到时拾取以及微震事件定位、微震事件时空分布规律研究，推测出岩体损伤程度、微裂隙位置，从而可以对岩体稳定性作出评估和预测，进而降低施工风险。

随着各类地下工程向着深部发展，微震监测技术得到了广泛的应用。与传统监测技术相比，微震监测技术具有远距离、动态、三维、整体、实时监测，能够分析不稳定体位置等优点[2]。目前，世界各国普遍将微震监测技术用于岩爆[3]、冲击地压[4-5]、突涌水[6]、滑坡[7]等地质灾害的预报预警，页岩气等非常规油气资源的开采[8]及二氧化碳储存库中[9]。Young等[10]基于加拿大原子能有限公司(Atomic Energy of Canada Limited, AECL)地下实验室(Underground Research Laboratory, URL)微震监测数据，提出了定量描述隧洞开挖过程中岩体弹性模量和破裂密度变化情况的方法。Storka等[11]论证了微震监测技术在二氧化碳捕获和储存库中的重要作用，它能够优化CO2注入参数，并能对破裂区活化、CO2泄漏等预报预警。陈炳瑞等[12]在锦屏Ⅱ级水电站3#引水隧洞隧道掘进机施工洞段开展了微震监测，获得了岩爆发生前微震活动的演化特征与规律，表明了微震技术预测预报岩爆发生的可行性。张文东等[13]对锦屏Ⅱ级水电站引水隧洞岩爆特征及微震监测规律进行了研究。李庶林等[14]将多通道微震监测系统应用于凡口铅锌矿，对深部采区的大爆破余震进行了监测。徐奴文等[7]将微震监测用于锦屏Ⅰ级水电站左岸边坡，提出了基于微震信息反演的岩体损伤模型，实现了边坡稳定性的定量评价。张伯虎等[15]将微震监测系统引入水电站地下厂房，建立了基于微震数据的地下厂房安全性评价方法。这些成果进一步发展了微震监测技术，扩展了微震监测技术的应用范围，积累了较多有关微震监测技术工程应用的经验。

尽管国内外学者对微震监测技术的工程应用进行了较多的研究，目前该技术的应用及相关理论的研究仍处于探索和发展阶段，主要侧重于对微震源定位等方面的研究，而基于隧道掘进机(tunnel boring machine, TBM)施工隧道开展微震监测规律的研究较少，如何充分挖掘监测数据以获得更多有用信息需要进一步研究。

本文以西藏在建的某公路隧道为研究对象，构建了微震监测系统，首先运用数值模拟方法验证了监测方案的可行性；进而，总结了深埋隧道TBM开挖过程中微震活动的时空演化规律及矩震级分布情况，并结合TBM掘进参数分析了微震事件数目及其能量释放与TBM掘进进尺并未同时达到最大值的原因；最后，揭示了微震事件b值与围岩失稳风险水平之间的关系，同时表明了b值表征隧道围岩稳定性的可行性。

西藏某公路隧道总长约为4.7 km，开挖洞径9.13 m，最大埋深832 m，采用双护盾TBM单头掘进施工，隧道地质剖面图如图 1所示。隧道沿线置于南迦巴瓦岩群多雄拉组的混合岩中，主要岩性为片岩和片麻岩。根据规划的试验成果，片岩的饱和抗压强度60~120 MPa，为坚硬岩；片麻岩的饱和抗压强度30~80 MPa，为中硬~坚硬岩。隧道最大垂直应力为22.4 MPa、最大水平主应力为28.89 MPa、最小水平主应力为18.74 MPa，二者叠加之后最大主应力超过30 MPa，为高地应力洞段。洞身段岩体新鲜、完整，且由于位于印度大陆楔入和俯冲到欧亚大陆的前缘，构造地应力高，围岩失稳问题突出，局部区域可能发生岩爆。因此，在现场搭建微震监测系统，对TBM施工过程中围岩稳定性进行监测，确保双护盾TBM施工安全十分必要。

数据采集仪和传感器是微震监测系统的重要硬件，如图 2所示。传感器采用单分量和三分量2种方法，其灵敏度较高(200 V/(m·s-1))，频带宽度较大(1.5~1 500 Hz)，固有频率为4.5 Hz。数据采集仪采用Daqlink3-24，它具有高分辨率、高采样率、高触发精度等特点。根据现场围岩情况，采集仪采样率选为2 000 Hz。采集仪内置有GPS接收器(时间精度为0.001 μs)，为仪器提供准确的时间和位置信息。

在隧道内分3个断面布置传感器，每个断面不同位置处布设3个传感器，将传感器置于深为1 m的孔底，并用腻子粉或锚固剂进行封堵，确保传感器与岩体充分耦合，具体布设如图 3所示。当TBM掘进20 m左右时，回收最后一个断面的传感器，经检测可靠后将其移至最前方，形成新的监测断面，从而确保掌子面始终处于传感器的监测范围之内，实现连续监测。

微震监测系统的选择关系到监测项目的成败，因此，对微震监测系统方案合理性进行评价十分必要。微震监测系统灵敏度及定位精度是评估微震监测方案可行性的重要指标，系统灵敏度是指系统能够精确定位的微震事件的最小矩震级，定位精度取决于系统灵敏度和监测目的，如表 1所示。

为了评估上述方案的合理性，对系统灵敏度及定位精度进行了数值模拟，结果分别如图 4和5所示。结果表明：距离传感器越远，系统灵敏度越小，定位误差越大；在隧道区域一定范围内，传感器断面之间系统灵敏度最大，能够精确定位最小矩震级为-2.22；隧道区域一定范围内，定位误差能控制在10 m以内。因此，上述微震系统能够满足深埋隧道中围岩失稳预警的需要。

微震监测能否在地下洞室中发挥岩爆预报预警的作用，主要取决于定位精度能否满足工程需求。选择合理的微震源定位方法对于提高微震源定位精度十分重要。

本次监测采用网格搜索法进行微震源定位，在速度模型(采用定位爆破震源进行反演)已知情况下，首先给出解(微震源位置坐标及发震时刻)的存在空间，然后将解空间进行网格划分，最后以台站实际走时与计算走时差为目标函数，通过遍历所有网格点来寻找解空间中使目标函数值取得最小的位置(全局最优解)，相应的网格点即为微震源的位置。监测期间(2016年9月10日—2017年1月7日)228个微震事件的定位误差分布情况如图 6所示，定位误差小于10 m的事件数占了70%，根据表 1中的定位精度推荐值，此结果能够满足隧道围岩稳定性评价的需要。

围岩失稳过程实质上是围岩内部微裂隙不断扩展演化的过程，而微震事件是岩体微破裂的反映，因此，通过分析微震活动的时空演化规律有助于了解岩体内部微裂隙的扩展情况，从而为评估地下洞室围岩稳定性提供前兆信息。

隧道开挖过程中，径向应力释放，岩体中积累的能量也得以释放，因此，微震事件的数量及其释放的能量，一定程度上能够反映隧道实际开挖情况。图 7和8分别为微震事件数量、微震事件释放能量与TBM日掘进进尺的关系曲线，其中，2016年10月28日至11月17日由于现场施工等原因，未开展微震监测。

图 7和8中反映了微震事件数量及其释放的能量随TBM掘进进尺的增大而增多的趋势，然而二者之间并非严格的正相关关系，且变化趋势在有些地方稍有滞后或提前，这可能是由于受到围岩地质条件的变异性的影响。微震事件是岩体在开挖卸荷过程中，伴随着裂缝的产生及扩展，储存的弹性能以应力波的形式释放而产生的。当岩体储能条件较好时，相同掘进进尺下因能量释放较多而产生较多微震活动，而掘进进尺会有所减少。隧道开挖过程中的微震活动性是围岩地质条件和掘进进尺以及其他因素综合作用的结果。归一化(min-max标准化)后的TBM掘进参数(见图 9)是围岩地质条件的反映，波动的掘进参数曲线反映了围岩地质条件的变异性。根据图 9可知，TBM掘进进尺最大时(图中阴影区域)，掘进参数表现为推力下降，刀盘扭矩增加，这反映出围岩储能条件变差的情况。这种情况下单位进尺微震事件较少，而有利于TBM掘进，表现出TBM掘进进尺最大，而微震事件未达最大的情况。因此，在围岩地质条件变异性的影响下，掘进进尺与微震事件数表现出不同步的现象。

2016年10月2—8日微震监测结果显示每天都有大量微震事件聚集于掌子面前方的桩号K10+253附近区域(图 10中黑色虚线圆圈)，并沿着图 10d中2条黑色实线聚集，微震事件能量释放突增，并进一步增加，如图 10a—10d所示。图 10中的圆表示微震事件，圆的直径表示微震事件的矩震级，直径越大，矩震级越大。微震事件的聚集说明桩号K10+253附近区域隧道围岩应力调整强烈，岩体内部微破裂较多，围岩失稳风险较高，该区域可能存在岩石剥落或局部岩体垮塌等现象。基于微震事件时空演化情况，10月8日16：00对桩号K10+253处隧道上方及右前方区域进行了围岩失稳预警，并建议密切关注掘进出渣量及渣体形态，选择合理的掘进参数和支护措施，及时掘进通过。根据预警在掘进过程中控制了出渣量、TBM的姿态，选择了重型管片进行及时支护，并在渣体连续长时间出现大块时通过超前钻孔及时加固了围岩。

2016年10月9日掘进至桩号K10+253区域时，TBM刀盘顶部及右侧出现塌方，塌腔沿洞线方向长8 m，顶部位于刀盘上方5 m左右，如图 11a和11b所示；伸缩护盾左侧岩体锤击清脆，强度较高，未见明显结构面，如图 11c所示；渣体以块状、粉状为主，粒径不均匀，如图 11d所示。

10月2—8日监测到的微震事件空间分布特征与10月9日发生的隧道岩体垮塌情况基本吻合。该微震监测系统投入使用前，发生的3次垮塌事故由于未提前采取措施，平均卡机时间长达30多d，本次垮塌事故因预警而提前采取了相应措施，对TBM损坏较小，卡机处理难度较小，恢复施工处理时间大大缩短(15 d)。

监测期间，微震事件的矩震级主要位于-2.8~0之间，如图 12所示。大量基于现场观测的研究结果表明[17-18]，微震事件的矩震级与频度也服从地震领域中的矩震级-频度(G-R)关系式，该公式是1944年Gutenberg与Richter[19]提出的有关地震震级与频度的统计关系，表达式为

其中：M为地震震级，N为震级大于等于M的地震次数，a、b为常数。b值表示大小震级事件数的比例关系，b值越小，大震级事件相对越多，b值的变化可以用来反映应力场的变化。b值可以反映微震活动震源特征，与断层滑动相关的微震活动b值较低，一般小于0.8；应力调整导致的微震活动b值通常在1.2~1.5之间[17]。

监测期间，2 128个微震事件的矩震级-频度曲线如图 12所示。由图 12可以看出，在最大震级、最小震级处存在两个取样偏差区域(虚线框处)。最小震级处存在取样偏差是由于监测系统灵敏度及小震级事件本身特点所致。由图 13可知，监测到的矩震级≤-1.7的小震级微震源多聚集在传感器附近，这是由于距离传感器阵列越远，系统灵敏度越小，能够精确定位微震事件的最小矩震级也越大。离传感器相对较远的小震级事件由于传播过程的衰减与耗散而未被监测到，从而导致最小震级处出现取样偏差。最大震级事件处存在取样偏差是由于大震级微震事件发生频度本身较低，导致大震级微震事件的不完备。因此，求取b值时应设置最大、最小震级阈值。最小震级阈值通常取震级-非累积频度曲线上频度最高处对应的震级(图 12中数值虚线所示的Mc)，最大震级阈值通常取震级-累积频度曲线上偏离线性趋势所对应的震级。由图 12可知，最小震级阈值取值偏小时，会导致b值变小；而最大震级阈值取值偏大时，会导致b值变大。

监测期间全部微震事件b值为1.24，如图 12所示，这表明监测期间微震活动主要由开挖卸荷引起的。这和图 7—8显示的9、10月份掘进进尺较少的情况并不矛盾，这是由于9、10月份TBM处于卡机状态，为了清除作用于TBM护盾上的岩体，掏渣清理、钻孔爆破、开挖导洞、围岩支护等TBM脱困措施实施频繁，这是另一种形式的开挖卸荷作用，在此过程中围岩应力重分布，产生较多微裂隙，由此产生较多微震事件。

不同月份微震事件b值、微震事件数及震级阈值如表 2所示。9月份监测开始时，TBM处于被卡状态，监测区域隧道顺右侧边墙大范围分布有长英质脉体，遇水易崩解，岩石强度较低，该区域围岩处于再平衡状态，发生大震级事件的可能性较小，b值最大。10月份TBM脱困后因掌子面附近围岩发生垮塌而再次被卡，微震事件矩震级相对较大且数目较多，围岩失稳风险较高，b值最小。随着TBM正常掘进，11—12月份微震事件均主要由因开挖卸荷、应力调整导致的，b值有所增加并保持相对稳定，反映出围岩失稳风险的降低。

综上，微震活动b值可以反映不同震级微震事件的活动情况，b值的变化可以反映围岩失稳风险水平，可以将其作为反映区域围岩整体稳定性的重要指标。

将微震监测系统引入高海拔深埋TBM施工隧道中，对隧道围岩稳定性进行实时监测与分析。基于对微震事件活动规律的研究，得出以下结论：建立了微震事件数及其释放能量与TBM日掘进进尺的关系，微震事件数目及其释放的能量随TBM掘进进尺的增大而增多，微震事件活动性能够反映TBM施工动态，然而当掘进进尺最大时，微震事件数并未达到最大值，这是由于围岩地质条件变异性的影响。岩体垮塌发生前，有大量微震事件聚集，并且随着时间的推移微震事件进一步增多，微震活动特征与现场实际情况符合较好，微震事件聚集区域与岩体垮塌位置吻合。监测期间，2 128个微震事件b值为1.24，这与监测期间微震活动基本由开挖卸荷引起的相吻合。b值可以反映不同震级微震事件的活动情况，b值的变化可以反映围岩失稳风险水平，可将其作为反映区域围岩稳定性的重要指标。