在隧道开挖的过程中，利用地震弹性波反射波法(Tunnel Seismic Prediction, TSP)对掌子面前方围岩的状况进行预报，可以有效地防止事故的发生，加快工程进度，降低成本，因此合理运用TSP超前预报系统对于隧道的安全施工具有重要意义。本文在介绍了TSP的基本原理的基础上，开展了隧道前方含低速异常体的TSP超前预报数值模拟，分析了异常体的地震波场响应特征，随后总结归纳了波速的影响因素。随后以山西某隧道TSP超前探测为例，推断了隧道前方的围岩情况，并通过现场开挖情况分析了该方法的应用效果，研究成果为该隧道以及类似隧道的安全施工提供有力保障。 In the process of tunnel excavation, Tunnel Seismic Prediction (TSP) is used to forecast the surrounding rock in front of the face, which can effectively prevent the occurrence of accidents, accelerate the progress of the project, and reduce the cost. Therefore, the reasonable application of TSP advanced prediction system is of great significance to the safe construction of tunnels. On the basis of introducing the basic principle of TSP, this paper carries out numerical simulation of TSP advanced prediction with low velocity abnormal body in front of tunnel, analyzes the characteristics of seismic wave field response of abnormal body, and then summarizes the influencing factors of wave velocity. Then, taking the TSP advanced detection of a tunnel in Shanxi Province as an example, the surrounding rock situation in front of the tunnel is inferred, and the application effect of the method is analyzed through the field excavation. The research results provide a strong guarantee for the safe construction of the tunnel and similar tunnels.
在隧道开挖的过程中，利用地震弹性波反射波法(Tunnel Seismic Prediction, TSP)对掌子面前方围岩的状况进行预报，可以有效地防止事故的发生，加快工程进度，降低成本，因此合理运用TSP超前预报系统对于隧道的安全施工具有重要意义。本文在介绍了TSP的基本原理的基础上，开展了隧道前方含低速异常体的TSP超前预报数值模拟，分析了异常体的地震波场响应特征，随后总结归纳了波速的影响因素。随后以山西某隧道TSP超前探测为例，推断了隧道前方的围岩情况，并通过现场开挖情况分析了该方法的应用效果，研究成果为该隧道以及类似隧道的安全施工提供有力保障。
TSP，超前地质预报，解释参数，数据处理板
Ruitao Ma
China Railway Langfang Geophysical Prospecting Co., Ltd., Langfang Hebei
Received: Nov. 25th, 2021; accepted: Dec. 9th, 2021; published: Dec. 27th, 2021
In the process of tunnel excavation, Tunnel Seismic Prediction (TSP) is used to forecast the surrounding rock in front of the face, which can effectively prevent the occurrence of accidents, accelerate the progress of the project, and reduce the cost. Therefore, the reasonable application of TSP advanced prediction system is of great significance to the safe construction of tunnels. On the basis of introducing the basic principle of TSP, this paper carries out numerical simulation of TSP advanced prediction with low velocity abnormal body in front of tunnel, analyzes the characteristics of seismic wave field response of abnormal body, and then summarizes the influencing factors of wave velocity. Then, taking the TSP advanced detection of a tunnel in Shanxi Province as an example, the surrounding rock situation in front of the tunnel is inferred, and the application effect of the method is analyzed through the field excavation. The research results provide a strong guarantee for the safe construction of the tunnel and similar tunnels.
Keywords:TSP, Advanced Geological Prediction, Interpretation Parameters, The Data Processing
Copyright © 2021 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着我国基础建设的不断推进，在公路以及铁路修建过程中，隧道工程建设项目越来越多。而在隧道施工建设过程中，由于地形地质条件复杂，导致前期勘察准确度较低，施工过程中的塌方、大变形事故频发。所以采用准确、快速的方法进行超前预报，获得隧道开挖前方的围岩情况，对隧道的安全施工具有重要意义。现阶段，常用的隧道超前预报方法有弹性波法(包括TSP、TRT、TST以及TGP)以及电磁法(GPR和TEM)两大类 [
尽管该方法已经广泛应用到全国范围内的隧道超前预报中，但实际上在隧道地质预报过程中，仍然存在一部分技术人员在为如何提高TSP法的准确性而大费周折，有时会出现数据处理解释准确性较差的问题，这就导致行业内认为TSP法的准确性距离要求还有差距，究其原因是多方面的 [
TSP法实质上就是地震反射波法，其与普通的地震方法主要区别在于施工空间 [
图1. TSP系统工作原理
地震波场常用的数值模拟方法有射线追踪法、积分方程法、有限单元法以及有限差分法。其中有限差分法具有计算速度快、精度高、适合复杂模型等优点，被广泛应用于求解波动方程 [
{ ρ ∂ v x ∂ t = ∂ τ x x ∂ x + ∂ τ x z ∂ z ρ ∂ v z ∂ t = ∂ τ x z ∂ x + ∂ τ z z ∂ z ∂ τ x x ∂ t = ( λ + 2 μ ) ∂ v x ∂ x + λ ∂ v z ∂ z ∂ τ z z ∂ t = λ ∂ v x ∂ x + ( λ + 2 μ ) ∂ v z ∂ z ∂ τ x z ∂ t = μ ( ∂ v x ∂ z + ∂ v z ∂ x ) (1)
其中， v x 和 v z 分别为x方向和z方向上的速度分量， ρ 为岩体的密度， λ 和 μ 为拉梅系数，t为采集时间， τ x x 、 τ z z 和 τ x z 分别为应力分量。
采用交错网格有限差分对上式进行离散，计算时分别将速度和应力放在交错网格的节点上，在分析了频散和稳定性条件之后，在边界出采用吸收边界条件，最终得到有限差分迭代公式 [
建立隧道前方存在低速异常体模型，二维模型尺寸为200 m × 200 m，隧道底板位于105 m位置处，高度为10 m，即顶板在95 m位置处；隧道长度为80 m，以模型边界为起点，掌子面位于横向80 m处。垂直低速异常体与掌子面的距离为40 m，该异常在二维剖面的横向位置为120~140 m、纵向位置为60~140 m，即低速异常体的宽为20 m，高80 m，如图2所示。围岩的具体物性参数如下：纵波速度为4000 m/s，横波速度为2000 m/s，密度为1800 kg/m3；隧道的物性参数如下：纵波速度为340 m/s，横波速度为0 m/s，密度为1 kg/m3；垂直低速异常体的物性参数如下：纵波速度为1500 m/s，横波速度为600 m/s，密度为540 kg/m3。观测系统为：震源采用300 Hz雷克子波，采样长度为0.02 s，采样点数为4000，采样间隔为0.05 ms。检波器和震源埋深在底板2 m深度处，偏移距为15 m，道间距1.5 m，检波器个数24个。
图2. 超前预报模型示意图
通过数值模拟可以得到不同时刻波场快照(t = 0.01 s、0.026 s、0.05 s和0.08 s)，如图3所示。从图中可以看出，当地震波没有传播到异常体时(图3(a))，不会产生反射波，而随着地震波的传播到异常体时(图3(b))，在异常体处会出现反射波、透射波以及绕射波，且随着波的传播，会出现多次反射波(图3(c)、图3(d))。图4所示为地震单炮记录，与波场快照显示一致，从地震记录中可以看到在0.026 s之前没有发现发射波，即地震波没有传播到异常体位置，随着采样时间增加，出现反射波以及多次波。对于隧道前方存在异常时的地震波场传播特征，可以为实际超前预报工作提供理论依据。
图3. 波场快照。(a) 0.01 s时刻波场快照；(b) 0.026 s时刻波场快照；(c) 0.05 s时刻波场快照；(d) 0.08 s时刻波场快照
图4. 地震记录
从上文数值模拟中可以发现，TSP法最直观的参数就是波速，包括纵波速度和横波速度，如果定义纵横波波速比为 A = v p / v s ，根据TSP法之间参数之间的关系，可以得到 [
{ E = ρ ⋅ v p 2 ⋅ 3 ⋅ A 2 − 4 A 2 ⋅ ( A 2 − 1 ) σ = A 2 − 2 2 ⋅ ( A 2 − 1 ) k = ρ ⋅ v p 2 ⋅ 3 ⋅ A 2 − 4 3 ⋅ A 2 μ = ρ ⋅ v p 2 ⋅ 1 A 2 λ = ρ ⋅ v p 2 ⋅ A 2 − 2 A 2 v s (2)
式中，E为杨氏模量， ρ 为岩体的密度， v p 和 v s 分别是纵波速度和横波速度， σ 为泊松比，k为弹性模量， μ 为剪切模量， λ 为拉梅系数。
从上式中可以发现其它参数包括泊松比、密度、动态杨氏模量等都是通过纵、横波速换算过来的，所以波速准确与否非常重要，如果波速不准确，其他的参数也不准确了，因此解释的基础参数应该还是纵波和横波的波速。而影响波速的因素较多，本文从岩性、孔隙度、岩石的密度、岩石的地质年代、地层的埋深以及初至波波速六个方面对地震波波速的影响进行总结归纳，为后续TSP法的数据解释提供理论依据，具体见表1。
| 序号 | 影响因素 | 影响效果分析 |
|---|---|---|
| 1 | 岩性 | 对于不同的岩石来讲，由于物性不同，导致地震波波速也不同。 ① 不同成因的岩石波速不同：火成岩和变质岩的波速要大于沉积岩的波速； ② 相同成因的岩石波速也不同：几种常见的沉积岩波速从大到小依次为石灰岩、页岩、砂岩。 ③ 同一种岩石的波速并不固定：受岩石的节理发育程度、风化程度、岩体的完整性以及强度等因素的影响，会导致同一种岩石的波速发生改变。 |
| 2 | 孔隙度 | 由于岩石的波速大于水中和空气中的速度，所以当岩石中的孔隙充填水或者是空气时，其波速会减小。岩石的孔隙度越大，地震波波速越小，即地震波波速与岩石的孔隙度成反比。除此之外，由于地震波在水中的速度要大于在空气中的速度，所以当孔隙度一定时，充填的空气越多时，波速越小。 |
| 3 | 岩石的密度 | 岩石的密度与孔隙度有一定的联系，当孔隙度越小，即密度越大，那么岩石的波速也就越大。通常情况，采用如下的经验公式来描述岩石密度与地震波波速之间的关系： ρ = 0.31 ⋅ v p 0.25 。 |
| 4 | 岩石的地质年代 | 通常情况下，岩石的地质年代对地震波波速的影响如下：地质年代越老，波速越大。除此之外，地震波波速还与构造运动有关：在挤压力作用下形成的褶皱区具有较大的速度，而在张力作用下形成的张性断裂则具有较小的速度。 |
| 5 | 地层的埋深 | 地层埋深实质是与岩石的密度呈正比，因为随着地层埋深的增加，受到上覆地层的重力作用下，岩石的密度越大，地震波波速也就越大。 |
| 6 | 初至波波速 | 由于在波速解释分析中，往往参照初至波波速(炮孔段波速)对该段数据进行分析解释，所以初至波波速的大小往往对后方探测的波速大小影响很大。因此，在探测后方数据的处理解释过程中，要注意不因波速的绝对值大小去判断，重点参照波速的相对变化趋势去判断前方地层变化。 |
表1. 影响地震波波速的因素总结
山西某隧道2斜11+94.0~2斜10+58.0段主要地层为黑云角闪斜长片麻岩，局部花岗岩脉，弱风化，节理裂隙很发育，岩体破碎，呈碎石状压碎结构，设计围岩级别为IV级。在隧道施工过程中，发现掌子面局部呈雨滴状出水，围岩稳定性差，易掉块、坍塌，所以采用TSP法对其进行超前预报工作，查明该段范围内的工程地质及水文地质情况，提供隧道施工建议，保证本段工程合理顺利的施工。
在布设激发孔时，只需要根据隧道内岩层的走向选择掌子面后方的左侧帮或者右侧帮一侧即可，激发孔的布设角度为斜向底板方向10˚~20˚，孔位高度为底板以上1.5 m，布置孔深为1.5 m，激发孔孔距为1.5 m，共24个。在距离掌子面最外侧的一个激发孔后方15~20 m位置处，分别在左、右两侧帮对称布置一个接收孔，孔位高度与激发孔一致，布置孔深为2 m，观测系统布置如图5和图6所示。
图5. 接收孔和激发孔平面布置图
图6. 接收孔和激发孔的截面图
在隧道2斜12+53.76的左、右边墙位置分别布置一个地震波信息接收孔，孔径约为50 mm，孔深均为1.92 m，孔高分别为1.05 m、0.95 m。在2斜12+04.40~2斜12+39.56段的左边墙，按约1.53 m的间距布置24个激发孔分别激发地震波，孔径约40 mm，孔深约为0.75~1.90 m；孔高约为1.05 m，每个激发孔装填的药量约150 g。
图7所示为Z分量地震原始记录，图8所示为TSP超前预报结果，从图中可以看到：1) 2斜11+94~2斜11+52段纵波速度和泊松比呈上升趋势，横波速度和动态杨氏模量呈下降趋势。推测该段围岩岩体破碎，局部较破碎，节理裂隙发育，围岩稳定性差。2) 2斜11+52~2斜11+07段纵波速度和泊松比呈下降–上升–下降趋势，横波速度和动态杨氏模量呈上升–下降–上升趋势。推测该段围岩岩体破碎，节理裂隙发育，围岩稳定性差。3) 2斜11+07~2斜10+58段纵横波速度和泊松比及动态杨氏模量均呈下降趋势。推测该段围岩岩体破碎，节理裂隙发育，局部节理裂隙很发育，围岩稳定性差。
图7. 实测地震记录
图8. TSP法超前预报成果图
对于2斜11+94~2斜11+52段，开挖揭示围岩为黑云角闪斜长片麻岩，节理裂隙发育~很发育，岩体破碎；2斜11+52~2斜11+07段，开挖揭示围岩为黑云角闪斜长片麻岩，节理裂隙很发育，岩体破碎；2斜11+07~2斜10+58段，当前开挖段落为黑云角闪斜长片麻岩，节理裂隙很发育，岩体破碎，其中2斜11+02出现塌方。现场验证情况与解释成果基本一致，表明本次超前预报结果准确。
通过介绍TSP法的基本原理和数值模拟，总结了影响波速的因素，结合工程实例可以得到以下结论。
1) 对比不同时刻的波场快照可以看到地震波在空间的传播特征，以及由异常体引起的反射波和绕射波，为TSP法在超前预报实际工作中的异常识别提供了理论依据。
2) 由于其他弹性参数包括泊松比、密度、动态杨氏模量等都是通过纵、横波速换算过来的，所以解释的最基础的参数是纵波和横波的波速。而由于地域及岩性的差异，同样的信号可能导致不同的解释结果。所以需要结合现场实际地质情况进行分析推断，只有不断地将地震波的岩石力学参数与地质知识有机地结合起来，才能更好地提高TSP探测的准确性。
3) 利用TSP法预报了山西某隧道2斜11+94.0~2斜10+58.0段的围岩情况，并通过了现场开挖证明，说明该方法可以准确预报出掌子面前方围岩破碎以及含水异常体，对于类似隧道超前预报工作具有重要的借鉴意义。
马瑞涛. TSP超前预报数值模拟及其在山西某隧道中的应用研究Numerical Simulation of TSP Advanced Forecast and Its Application in a Tunnel in Shanxi Province[J]. 土木工程, 2021, 10(12): 1247-1255. https://doi.org/10.12677/HJCE.2021.1012137