Application of Rayleigh Wave Elliptic Polarization in Shield Tunnel Stratum Disturbance Exploration
-
摘要:
为解决传统物探方法在城市盾构隧道工程中难以对地层扰动达到较好探测和监测目的，提出了主动源瑞利面波椭圆极化方法. 采用二维有限元数值模拟对实际盾构隧道施工背景下的地层情况进行全波场正演模拟，并对其波场传播、速度频散和椭圆极化频散等特征进行分析；结合工程实例，进一步将椭圆极化方法和传统速度频散方法进行对比研究. 研究结果表明：主动源瑞利面波椭圆极化方法在探测地下介质结构分布及监测地下空洞等异常结构的变化方面具有和传统速度频散方法同样的探测效果，但由于其不需要长大检波器排列，对周围环境震动等噪声抗干扰能力相对较强，因此更适合城市等场地受限和具有一定震动和噪声干扰的环境，从而具有更广阔的应用前景.
Abstract:In order to solve the problem that traditional geophysical methods are difficult to achieve better detection and monitoring of stratum disturbance in urban shield tunnel engineering, the application of active source Rayleigh surface wave elliptical polarization method is proposed. Firstly, the 2D finite element numerical simulation is used to simulate the full-wave field forward modeling of the stratum under the actual shield tunnel construction background, and the characteristics of wave field propagation, velocity dispersion and elliptical polarization dispersion are analyzed. Then, the elliptical polarization method and the traditional velocity dispersion method are further analyzed and compared on the basis of a combination of examples. The results indicate that the active source Rayleigh surface wave elliptic polarization method has the same detection effect as the traditional velocity dispersion method in detecting the distribution of underground medium structure and monitoring the changes of underground cavity and other abnormal structures. However, because it does not need the arrangement of long and long geophones, it has relatively strong anti-interference ability to the vibration and other noise of the surrounding environment. Therefore, it is more suitable for the environment with limited sites such as cities and certain vibration and noise interference, so it has a broader application prospect.
-
随着国内城市盾构隧道的发展，盾构施工引发的地面塌陷事故对人民生命财产安全造成严重威胁[1-2]. 为避免此类事故的发生，在盾构施工过程中对其产生的变形或沉降[3]等进行了监控测量[4]，但这些测量都是在沉降或变形等发生后才能进行的，因此很难实现提前预警. 近年来，人们开始采用物探方法，通过对介质物性参数变化的观测，实现对介质变形、沉降以及可能形成的空洞和塌陷等的探测和监测，如探地雷达、瑞利面波法等[5-9]. 但探地雷达精度虽高，勘探深度却有限，传统瑞利面波速度频散方法由于需要长大检波器排列，因此在数据采集上易受场地限制，尤其是人口密集、交通繁忙的城市等地. 为解决传统速度频散瑞利面波法在应用中存在的问题，研究者们开始关注瑞利面波另一重要特征，即椭圆极化特征，该方法可通过单个三分量检波器采集信息[10]，不易受场地限制，能更好适应城市背景下的检测或勘探任务.
早在1969年，Boore等[11]首先提出利用瑞利面波椭圆极化率来研究地壳结构的可行性；Taner等[12]计算了单个三分量检波器采集的地脉动信号水平与垂直分量的傅里叶变换频谱比，并提出了单点H/V谱比法（H、V分别为地震记录面波时间域信号水平分量、垂直分量）；Tokimatsu和Miyadera[13]发现，地脉动H/V谱比对场地横波速度Vs变化的响应特征与基阶瑞利面波椭圆极化率对Vs的响应特征是一致的；蒋通等[14]讨论了地脉动实测H/V计算方法的可靠性，并据此评价场地的卓越周期；Arai等[15]指出地脉动H/V谱比主要反映的是瑞利面波或勒夫波的椭圆极化率，其对地层速度变化的灵敏度比速度频散高；张立[16]推导了水平层状介质中瑞利面波质点位移解析解，结合地脉动单点谱比法通过数值模拟得到层状介质中主动源瑞利面波的椭圆极化频散，展示了利用单点瑞利面波的多分量评价地层泊松比结构的可行性；于文福[17]利用有限元方法开展了横向非均匀介质中瑞利面波椭圆极化率频散数值模拟研究，验证了瑞利面波椭圆极化率与介质泊松比之间的线性反比关系.
目前对瑞利面波椭圆极化的相关研究和应用尚处于探索阶段，在此前提下，依据城市盾构隧道工程实例，通过结合传统瑞利面波速度频散方法与主动源椭圆极化频散方法对其进行数值模拟研究和实际案例应用对比，探讨主动源瑞利面波椭圆极化频散方法在城市盾构隧道上覆地层扰动探测和地面塌陷监测中应用的可行性，为今后类似探测和监测工作提供参考.
1. 基本原理
1.1 有限元数值模拟原理
基于地震波原理，应用有限元数值模拟软件，对城市地下盾构隧道开挖引起地层扰动的地震全波场特征进行数值模拟研究[18-19].
本文研究采用的是零相位、主频24 Hz的宽带雷克子波作为震源函数，如式（1）.
u(t)=[1−2π 2f2(t−t1)2]e−π 2f2(t−t1)2， (1) 式中：u(t)为载荷大小，t为时间；t1为雷克子波主峰时间；f为雷克子波主频.
1.2 瑞利面波频散曲线提取方法
通常情况下地震波并非单一频率的简谐波，而是由多个不同频率的波叠加而成，此频散现象存在于非均匀介质中[20-21]. 由于瑞利面波产生在自由表面附近，其振幅随深度衰减，穿透深度约为一个波长，因此不同频率穿透深度不同. 高频即较短波长的波采样浅层，相反，低频即较长波长的波采样深层. 因此，不同频率波的传播速度各异，地层速度也随深度变化，这种传播速度随频率变化的现象，即为速度频散Vr(f).
频率-波数（f-k）变换是目前提取面波频散曲线应用较为广泛的方法之一. 通过对原始x-t域地震记录作二维傅里叶正变换使其变换到f-k域，其中二维傅里叶正变换式为
G(ω,k)=∫+∞−∞∫+∞−∞u(x,t)e−i(ωt−kx)dtdx, (2) 式中：ω为圆频率；u(x，t)为原始地震记录，G(ω，k)为经过二维傅里叶变换后得到的地震数据.
f-v （频率-速度）变换则是在f-k变换基础上，利用波数和速度、频率之间的关系变换到f-v域再进行频散曲线提取，其关系式为
k=fv. (3) 1.3 瑞利面波椭圆极化提取方法
根据前期研究成果，与速度频散特性相似，在非均匀介质中瑞利面波椭圆极化率Er也会随着频率f的变化而变化[22]，即为椭圆极化频散Er（f），然而对于均匀介质，椭圆极化率是一个固定常数. 目前计算瑞利面波椭圆极化频散的方法主要是H/V谱比法[23-24]，其计算式为
Er(f)=H(f)V(f), (4) 式中：H(f)为水平位移频谱；V(f)为垂直位移频谱.
对于瑞利面波在自由表面传播，其椭圆极化率E0的计算式为
E0=(n−m)2+1−m2m(1−n2), (5) 式中： m = 1−VR2/Vp2；n = 1−VR2/Vs2 ；Vp、Vs、VR分别为纵波速度、横波速度和面波速度.
2. 有限元数值模拟
2.1 数值模拟模型
盾构隧道施工时开挖面支护压力较小，掘进后容易在盾构上方形成空洞，随着土拱作用，空洞顶部土体逐渐向下脱落，空洞逐渐向上方地表发展，最后导致地表塌陷[25-26]. 本文根据上述盾构隧道施工引发的滞后地面塌陷机理，结合长沙万家丽路220 V电力隧道里程300—370段实例施工方案及实际地质等情况，首先从均匀模型（图1（a））瑞利面波数值模拟研究出发，并根据实际上覆地层埋深8~15 m和盾构隧道外径4.1 m情况，建立盾构隧道模型（图1（b）和图1（c））、盾构开挖后形成的空洞初期模型（图1（d））以及空洞向地表发展模型（图1（e））.
图 1 盾构隧道数值模型观测排列示意Figure 1. Observation arrangement for numerical model of shield tunnel本文数值模拟研究采用的是多道瞬态瑞利面波观测排列[27]，最小偏移距为40 m，共51个检波器沿测线依次排列，道间距为1 m. 采样时间1.024 s，时间步长为0.0005 s. 模型中原状粉质黏土和疏松状粉质黏土分别用编号1、2代表，介质参数的取值以实际地质情况为参考，因无特殊岩性与地质构造，故采用一般均质岩土弹性参数，其模型主要参数如表1.
表 1 数值模型介质参数Table 1. Medium parameters of numerical model编号 弹性
模量/
MPa密度/
（kg•m−3）纵波
波速/
（m•s−1）横波
波速/
（m•s−1）面波
波速/
（m•s−1）泊松比 1 280 1800 630 234 221 0.42 2 11 1500 524 160 151 0.45 2.2 结果与分析
基于有限元软件对图1中各模型进行二维模拟仿真，并对模拟结果进行波场传播、速度频散以及椭圆极化频散等特征的相关研究，特征分析如表2所示.
表 2 数值模拟结果特征分析Table 2. Characteristic analysis of numerical simulation results模型 波场传播特征 速度频散特征 椭圆极化特征 a 直达纵波和面波同相轴均为直线，且随着偏移距增大，纵波衰减较快，面波衰减较慢 Vr 基本为200~220 m/s，与均匀介质模型Vr 理论值 221 m/s 基本吻合 椭圆极化率值 Er 趋于均一，整体 Er 值分布在 0.8 附近 b 直达纵波和面波之间以及面波下方都出现了弯曲弧形同相轴（实线圈处），且弧形同相轴顶点位于第25 个接收点附近 黑色实线区内（偏移距65 m、深度10 m 附近）介质 Vr 为130~160 m/s，呈相对低速特征，其周围介质整体速度相对较为均匀，约190~220 m/s，为相对高速特征，该高速值与设置的围岩 Vr 值221 m/s 相符合 隧道断面（黑色实线区）周围出现分别向两侧斜下方倾斜的 Er 异常，其左侧（左虚线圈）为向左下方微倾斜的散射状Er（1.0~1.4）区域，右侧（右虚线圈）为向右下方微倾斜的散射状低 Er（0.1~0.5）区域，黑色实线区 Er 为 0.6~0.8，介于两侧 Er 之间 c 直达纵波 DP 下方，每隔一定 x-t 间隔出现一组振幅较强的同相轴波形（实线圈范围内） 以10 m 处白色实线为界有明显的速度差异，隧道（白色实线区域）及以下呈明显低速特征，其 Vr 约为110~160 m/s，隧道上层介质速度为200~235 m/s，与模型设置围岩 Vr 221 m/s 接近 深度10 m 即隧道界面附近，椭圆极化特征大致呈上下两层分布，其中上层椭圆极化率值高低较不均匀，呈区块状分布，其中较高（1.3~2.5）或较低（0.1~0.7）的 Er 值区块，整体上，上层 Er 大于下层 d 空洞位置即第25个检波器偏移距65 m附近，体波下方及面波分布区中叠加有弯曲的弧形同相轴（实线圈处）波形 空洞区域及隧道（上下两个白色实线框）下方介质均为低速特征，面波波速 Vr 为120~170 m/s，上方介质Vr 较高为195~235 m/s，与设置的围岩面波速度 221 m/s 较吻合 偏移距 65 m 处即空洞（白色实线小框）下方 Er 值分别向两侧斜下方呈放射状高低相间排列的分布规律，其中左侧为 Er 高值区（左虚线圈所示， 1.0~2.0），右侧为 Er 低值区（右虚线所示，0.3~0.9），空洞异常在两侧放射状高低 Er 之间的顶部位置；且 Er 分布规律整体上呈上下两层分布，上侧 Er 明显高于下侧 e 空洞位置即第25 个检波器偏移距65 m 附近，体波下方及面波分布区中叠加有弧形同相轴（实线圈处） 空洞运动轨迹范围（白色实线区域和白色虚线区域）内堆积土体介质的低速特征，其速度约为 120~170 m/s，周围介质速度仍为原来设定的190~230 m/s 左侧 Er 高值区（左虚线圈所示）为1.6~2.5，右侧 Er 低值区（右虚线圈所示）0.3~1.2，空洞异常仍在两侧高低 Er 之间的顶部位置，且 Er 分布规律整体上呈上下两层分布，上侧 Er 值明显高于下侧 2.2.1 均匀介质模型
在图2（a）中：通过偏移距与时间的关系可以算出波速，与模型设定（表1）的理论速度630 m/s （Vp）和221 m/s （Vr）基本吻合；根据研究发现，对于瑞利面波无论选用水平分量还是垂直分量，提取的频散曲线基本是一致的[28]. 图2（b）中：由于体波等干扰，图中近源场区可看到面波速度有一定的低速区，而在偏移距较大、面波能量占优的远场区，速度频散剖面显示则较为均匀，即没有频散现象；由已有研究[17]可知，均匀介质中震源近场区由于体波等干扰，椭圆极化率值Er会呈现放射状高低值交替现象，随着偏移距增大，即在远场区，这种现象逐渐减弱. 从图2（c）中可看出：随着偏移距的增大，Er整体趋于均一，说明在均匀介质震源远场区， Er与速度频散特征类似，也是一个常数，即没有频散现象.
2.2.2 盾构隧道横剖面模型
图 3 盾构隧道横剖面模型数值模拟成果Figure 3. Numerical simulation results of transverse section model of shield tunnel由图3可知：通过能量大小关系可以清楚辨别出最上方直线同相轴为直达纵波DP，其下方能量最强的直线同相轴是面波R，其弧形同相轴顶点位于第25个接收点附近，即对应偏移距65 m处、隧道横剖面位置（虚线处），则可知弧形曲线是直达波和面波在遇到隧道横断面时发生的反射波或转换波等；速度频散分析能够较准确地反映出隧道横剖面及周围介质的相对位置；此模型模拟结果（图3（c））与邓瑞[29]研究中所述均匀介质中局部空洞模型的Er异常呈“左高右低”现象[29]相一致. 因此，由数值模拟结果可知，盾构隧道横剖面上Er具有特定异常分布现象，且Er与介质泊松比之间存在一定对应关系[17,22]，即椭圆极化频散方法能够有效识别隧道在横剖面上的分布结构.
2.2.3 盾构隧道纵剖面模型
图 4 盾构隧道（纵剖面）模型数值模拟成果Figure 4. Numerical simulation results of longitudinal section model of shield tunnel图4（a）中：由于沿隧道纵剖面方向，震源激发产生的地震波主要在上层区域内即隧道上壁界面和地表之间来回反射叠加，并在界面和地表处发生波型转换和生成新的面波，从而上层区域内各种波的叠加使得局部范围内波振幅增大或减小，因此在地震记录图上呈现为隔一定时间和距离即x-t间隔出现一组同相轴的波形[30]. 图4（b）中：由于隧道挖空，沿隧道纵剖面方向地震波能量难以穿透空气向下传播，可看出隧道及以下区域都呈明显低速特征；可发现传统速度频散方法能够较好地反映出隧道所处位置及介质分层情况. 图4（c）中：椭圆极化值Er呈区块状不均匀分布，与上述地震记录中出现的时强时弱的同相轴地震波密切相关，均为地震波在上层区域内来回反射叠加使局部范围内波振幅增大或减小，从而造成局部区域内按H/V谱比方法计算所得的Er值偏大或偏小的区块现象；隧道上方介质Er虽呈较不均匀的区块分布，但整体明显比下方介质高. 因此， Er在隧道纵剖面结构分布上具有明显差异特征，能够对隧道纵剖面分层情况进行有效识别.
2.2.4 盾构隧道上方空洞初期模型
对空洞形成初期模型（图1（d））进行数值模拟，结果如图5所示.
图 5 盾构隧道上方空洞初期模型数值模拟成果Figure 5. Numerical simulation results of initial cavity model above shield tunnel由图5可知：图5（a）中地弧形曲线是地震波遇到隧道上方空洞时发生的反射波所形成；速度频散方法对隧道纵剖面分层及开挖形成的初期空洞都能有效识别；与图3（c）和图4（c）对比可发现，图5（c）上新形成的空洞位置附近出现了明显类似图3（c）的椭圆极化率异常现象，且隧道上层介质中其他高Er值区块的形成与图4（c）类似，是地震波在地表和隧道壁之间相互反射叠加的结果. 因此，椭圆极化频散方法能够较好地分辨出盾构隧道在纵剖面上的结构情况，并且能有效识别出隧道上方形成的局部空洞异常.
2.2.5 盾构隧道上方空洞向地表发展模型
盾构隧道开挖引起的上方空洞随着时间推移逐渐向地表发展，对此时空洞发展模型（图1（e））的数值模拟结果如图6所示.
图 6 盾构隧道上方空洞向地表发展模型数值模拟成果Figure 6. Numerical simulation results of development model of shield tunnel cavity to ground surface图6（a）通过与图5（a）的对比可发现：如表2模型e所述的弧形同相轴特征出现的时间更早，弧形特征更为明显. 造成该现象的原因一方面是因为空洞上移后埋深变浅，另一方面是因为空洞发展过程中掉落的土体相对周围介质较为疏松，从而使空洞发展轨迹范围内都属于低速区域，因此随着低速区域范围的变浅和增大，弧形反射波特征出现的时间也更早、更明显. 通过对比图5（b）可见：图6（b）的速度频散
方法对隧道纵剖面结构分层及隧道上方局部空洞和空洞发展状态均有非常好的识别效果. 如图6（c）所显示，通过与图4（c）及图5（c）对比可看到：随着空洞位置上移，空洞和低速异常区域范围扩大，空洞下方Er值分别向两侧斜下方呈高低相间放射状排列的现象出现的时间更早、特征更加明显. 通过Er分布规律能够较好分辨盾构隧道纵剖面结构的基础上，能够进一步有效识别出隧道上方局部空洞分布及规模，即该方法能够对介质中空洞的发展状态进行有效监测. 通过研究发现：Er频散和FV方法都能够对隧道开挖引起的上方空洞及空洞发展引起的地表塌陷进行有效探测和监测，且Er频散方法具有不需要长大的检波器排列的优点，因此更适合城市等场地受限的情况. 对图1（a）、（b）、（c）模型的数值模拟正是参照本文实际工程案例的施工方案和地质情况，分别对盾构隧道开挖前和开挖后纵剖面与横剖面实际情况进行仿真研究，而下文则是对实际案例采集数据的相同处理分析.
3. 应用实例
3.1 工程概况
以长沙万家丽路220V电力隧道里程300—370段为例. 本工程北起国防科技大学三号院北侧马栏山变电站，向东铺设至万家丽路、特立路交叉路口，本段地貌单元属于丘陵地貌，在特立路口向南拐弯至开福区福元路，沿既有万家丽路向南，纵贯开福区、止于芙蓉区之火炬路，主要为浏阳河冲积阶地，地貌单元为典型河流侵蚀堆积地貌. 该段采用盾构法施工，隧道内径3.6 m，外径4.1 m，电力隧道顶面覆土厚度为8~15 m.
3.2 采集参数
由于现场公路交通繁忙，测线布置在隧道顶部地表沿轴线走向. 采用DAQlink-Ⅲ分布式地震仪，主频4 Hz的三分量检波器，选定主动源方式采集数据，采样时间0.512 s，采样率0.0625 ms，测量时间分别在开挖前原状态时和开挖7 d后地表相对位移沉降稳定时进行，具体测线布置如图7所示.
3.3 数据处理及解释
首先对开挖前和开挖后第7 d两次采集的数据进行前期编辑和预处理后：对垂直分量数据经过FK滤波（率波数域滤波）等处理后提取速度频散曲线，绘制速度频散剖面图；对水平分量和垂直分量经过频谱计算等处理后，对其进行频谱比即Er计算，并绘制Er频散剖面图，最后结合地质资料等对隧道开挖前后速度频散和Er频散成果图进行分析和解释.
3.3.1 盾构隧道开挖前
隧道开挖前垂直分量速度频散剖面及HV分量椭圆极化Er频散剖面分别如图8（a）、（b）所示，从图8可知：在隧道深度范围以内，整个测区在隧道开挖前速度和Er值分别为190~250 m/s和4.4~5.0，面波速度Vr和Er整体分布都较为均匀，没有明显大范围的空洞等低速或低泊松比的异常现象，该现象与前期勘察地质资料结果相符.
图 8 实例（盾构隧道开挖前）解译Figure 8. Projecte interpretation diagram （before shield tunnel excavation）3.3.2 盾构隧道开挖后
盾构隧道开挖后垂直分量速度频散剖面和HV分量Er频散剖面如图9所示，其中白色水平线之间为隧道位置. 由图9可知：在隧道上壁界面附近出现V和Er分界面，其中速度频散剖面中隧道下层介质呈明显低速特征，Er则为低泊松比特征.
图 9 实例解译（盾构隧道开挖后）Figure 9. Projecte interpretation diagram （after shield tunnel excavation）由两种方法对比可看到，瑞利面波椭圆极化频散方法同样可以达到速度频散方法对地下岩土介质结构有效判识的效果. 此外，根据前面对隧道开挖形成的空洞初期及发展中状态的数值模拟特征分析，此次工程实例中未见有明显类似空洞的特征出现，结合现场地质勘察等资料证实确实没有空洞存在.
图9（a）中低速区块可能是因为在交通繁忙的公路边，现场采集数据质量较差所造成的速度频散曲线拾取产生偏差，相对而言，通过对HV分量频谱比计算所得的Er频散方法对此类干扰具有较强的抵抗性，因此在对应的Er频散图上没有明显的空洞Er放射状现象.
根据以上工程实例研究发现：椭圆极化频散方法不仅具有传统速度频散方法类似的探测效果，能够探测和监测地下介质的结构分布及比变化状态，并且不需要长大检波器排列，对现场震动噪声也具有较强的抗干扰性，由此可见，椭圆极化频散方法具有更广的应用范围和应用前景.
4. 结 论
本文在以城市盾构隧道上覆地层扰动探测和地面塌陷监测为例，对其进行了正演数值模拟研究，并结合工程实例进行了应用研究,对比分析得到以下结论：
1） 数值模拟结果中，对于地下介质的结构分布以及横向非均匀体如空洞等异常体的探测，相比传统成熟的速度频散方法，主动源利面波椭圆极化方法，具有类似的探测效果，且该方法可实现单点采集，不受场地限制，并且处理步骤更为快捷，受周围噪声干扰影响相对较小，因此在采集和处理上具有更高的效率.
2） 相比传统方法不能对城市盾构隧道施工中出现的地面塌陷等情况进行及时预警的问题，瑞利面波椭圆极化方法可以在不受场地限制的情况下，对上覆地层受到扰动后结构发生变化的情况进行无损探测和监测，因此可以对可能发生的塌陷等灾害进行及时预警，从而可以避免重大的人民生命财产损失.
-
图 1 盾构隧道数值模型观测排列示意
Figure 1. Observation arrangement for numerical model of shield tunnel
图 3 盾构隧道横剖面模型数值模拟成果
Figure 3. Numerical simulation results of transverse section model of shield tunnel
图 4 盾构隧道（纵剖面）模型数值模拟成果
Figure 4. Numerical simulation results of longitudinal section model of shield tunnel
图 5 盾构隧道上方空洞初期模型数值模拟成果
Figure 5. Numerical simulation results of initial cavity model above shield tunnel
图 6 盾构隧道上方空洞向地表发展模型数值模拟成果
Figure 6. Numerical simulation results of development model of shield tunnel cavity to ground surface
图 8 实例（盾构隧道开挖前）解译
Figure 8. Projecte interpretation diagram （before shield tunnel excavation）
图 9 实例解译（盾构隧道开挖后）
Figure 9. Projecte interpretation diagram （after shield tunnel excavation）
表 1 数值模型介质参数
Table 1. Medium parameters of numerical model
编号 弹性
模量/
MPa密度/
（kg•m−3）纵波
波速/
（m•s−1）横波
波速/
（m•s−1）面波
波速/
（m•s−1）泊松比 1 280 1800 630 234 221 0.42 2 11 1500 524 160 151 0.45 
下载: 导出CSV
表 2 数值模拟结果特征分析
Table 2. Characteristic analysis of numerical simulation results
模型 波场传播特征 速度频散特征 椭圆极化特征 a 直达纵波和面波同相轴均为直线，且随着偏移距增大，纵波衰减较快，面波衰减较慢 Vr 基本为200~220 m/s，与均匀介质模型Vr 理论值 221 m/s 基本吻合 椭圆极化率值 Er 趋于均一，整体 Er 值分布在 0.8 附近 b 直达纵波和面波之间以及面波下方都出现了弯曲弧形同相轴（实线圈处），且弧形同相轴顶点位于第25 个接收点附近 黑色实线区内（偏移距65 m、深度10 m 附近）介质 Vr 为130~160 m/s，呈相对低速特征，其周围介质整体速度相对较为均匀，约190~220 m/s，为相对高速特征，该高速值与设置的围岩 Vr 值221 m/s 相符合 隧道断面（黑色实线区）周围出现分别向两侧斜下方倾斜的 Er 异常，其左侧（左虚线圈）为向左下方微倾斜的散射状Er（1.0~1.4）区域，右侧（右虚线圈）为向右下方微倾斜的散射状低 Er（0.1~0.5）区域，黑色实线区 Er 为 0.6~0.8，介于两侧 Er 之间 c 直达纵波 DP 下方，每隔一定 x-t 间隔出现一组振幅较强的同相轴波形（实线圈范围内） 以10 m 处白色实线为界有明显的速度差异，隧道（白色实线区域）及以下呈明显低速特征，其 Vr 约为110~160 m/s，隧道上层介质速度为200~235 m/s，与模型设置围岩 Vr 221 m/s 接近 深度10 m 即隧道界面附近，椭圆极化特征大致呈上下两层分布，其中上层椭圆极化率值高低较不均匀，呈区块状分布，其中较高（1.3~2.5）或较低（0.1~0.7）的 Er 值区块，整体上，上层 Er 大于下层 d 空洞位置即第25个检波器偏移距65 m附近，体波下方及面波分布区中叠加有弯曲的弧形同相轴（实线圈处）波形 空洞区域及隧道（上下两个白色实线框）下方介质均为低速特征，面波波速 Vr 为120~170 m/s，上方介质Vr 较高为195~235 m/s，与设置的围岩面波速度 221 m/s 较吻合 偏移距 65 m 处即空洞（白色实线小框）下方 Er 值分别向两侧斜下方呈放射状高低相间排列的分布规律，其中左侧为 Er 高值区（左虚线圈所示， 1.0~2.0），右侧为 Er 低值区（右虚线所示，0.3~0.9），空洞异常在两侧放射状高低 Er 之间的顶部位置；且 Er 分布规律整体上呈上下两层分布，上侧 Er 明显高于下侧 e 空洞位置即第25 个检波器偏移距65 m 附近，体波下方及面波分布区中叠加有弧形同相轴（实线圈处） 空洞运动轨迹范围（白色实线区域和白色虚线区域）内堆积土体介质的低速特征，其速度约为 120~170 m/s，周围介质速度仍为原来设定的190~230 m/s 左侧 Er 高值区（左虚线圈所示）为1.6~2.5，右侧 Er 低值区（右虚线圈所示）0.3~1.2，空洞异常仍在两侧高低 Er 之间的顶部位置，且 Er 分布规律整体上呈上下两层分布，上侧 Er 值明显高于下侧
下载: 导出CSV
-
[1] 张成平,张顶立,王梦恕,等. 城市隧道施工诱发的地面塌陷灾变机制及其控制[J]. 岩土力学,2010,31(增刊1): 303-309.ZHANG Chengping, ZHANG Dingli, WANG Mengshu, et al. Catastrophe mechanism and control technology of ground collapse induced by urban tunneling[J]. Rock and Soil Mechanics, 2010, 31(S1): 303-309. [2] 侯艳娟,张顶立,李鹏飞. 北京地铁施工安全事故分析及防治对策[J]. 北京交通大学学报,2009,33(3): 52-59. doi: 10.3969/j.issn.1673-0291.2009.03.011HOU Yanjuan, ZHANG Dingli, LI Pengfei. Analysis and control measures of safety accidents in Beijing subway construction[J]. Journal of Beijing Jiaotong University, 2009, 33(3): 52-59. doi: 10.3969/j.issn.1673-0291.2009.03.011 [3] 郭汉章. 武汉地铁七号线盾构隧道施工对地面沉降及建筑物变形影响分析[D]. 武汉: 武汉工程大学, 2017. [4] 周松,荣建,陈立生,等. 单管双线大直径地铁盾构隧道施工与监测[J]. 地下空间与工程学报,2013,9(2): 365-372.ZHOU Song, RONG Jian, CHEN Lisheng, et al. Construction and monitoring of single underground railway tunnel with dual carriageway using large-diameter shield[J]. Chinese Journal of Underground Space and Engineering, 2013, 9(2): 365-372. [5] 王红庭. 盾构隧道注浆病害探地雷达检测实验研究[D]. 武汉: 华中科技大学, 2019. [6] 刘宏岳,黄佳坤,孙智勇,等. 微动探测方法在城市地铁盾构施工“孤石”探测中的应用: 以福州地铁1号线为例[J]. 隧道建设,2016,36(12): 1500-1506. doi: 10.3973/j.issn.1672-741X.2016.12.014LIU Hongyue, HUANG Jiakun, SUN Zhiyong, et al. Application of microtremor method to boulders detection in urban metro shield construction: case study of Fuzhou metro line No. 1[J]. Tunnel Construction, 2016, 36(12): 1500-1506. doi: 10.3973/j.issn.1672-741X.2016.12.014 [7] SHEN L Y, CHEN X W, ZHANG X, et al. Seismic characterization for underground construction projects using microtremor survey method: a case study in Shenzhen Subway Line 10[C]//SEG Technical Program Expanded Abstracts 2017. Houston: Society of Exploration Geophysicists, 2017: 5453-5457. [8] SHEN L Y, CHEN X W, ZHENG Y, et al. Seismic characterization for underground construction projects using microtremor survey method: a case study in Chengdu Line 18[C]//SEG Technical Program Expanded Abstracts 2018. Anaheim: Society of Exploration Geophysicists, 2018: 1-5. [9] 张新. 基于随机源模型的新型微动探测技术在上软下硬地层盾构隧道中的应用研究[J]. 现代隧道技术,2020,57(5): 43-50.ZHANG Xin. Research on application of random source model based new microtremor survey technology to shield tunnelling in upper soft and lower hard strata[J]. Modern Tunnelling Technology, 2020, 57(5): 43-50. [10] 张立,刘争平,熊自英. 利用地脉动H/V谱比反演地层Vs结构剖面[J]. 工程地球物理学报,2009,6(增1): 4-5.ZHANG Li, LIU Zhengping, XIONG Ziying. S-wave velocity profiling by inversion of microtremor H/V spectra[J]. Chinese Journal of Engineering Geophysics, 2009, 6(S1): 4-5. [11] BOORE D M, NAFI TOKSÖZ M. Rayleigh wave particle motion and crustal structure[J]. Bulletin of the Seismological Society of America, 1969, 59(1): 331-346. doi: 10.1785/BSSA0590010331 [12] TANER M, SHERIFF R. Application of amplitude, frequency, and other attributes to stratigraphic and hydrocarbon determination: section 2, application of seismic reflection configuration to stratigraphic interpretation[J]. Railway Technical Reasearch Institute, Quarterly Reports, 1989, 30(1): 20-24. [13] TOKIMATSU K, MIYADERA Y. Characteristics of Rayleigh waves in microtremors and their relation to underground structures[J]. Journal of Structural and Construction Engineering (Transactions of AIJ), 1992, 439: 81-87. doi: 10.3130/aijsx.439.0_81 [14] 蒋通,崔哲,李文艺. 层状介质中的瑞雷面波理论及其在场地脉动观测中的应用[J]. 上海力学,1999,20(3): 252-258.JIANG Tong, CUI Zhe, LI Wenyi. Theory of Rayleigh wave in multi-layered media and its application in microseism observation[J]. Shanghai Journal of Mechanical, 1999, 20(3): 252-258. [15] ARAI H. S-wave velocity profiling by joint inversion of microtremor dispersion curve and horizontal-to-vertical (H/V) spectrum[J]. Bulletin of the Seismological Society of America, 2005, 95(5): 1766-1778. doi: 10.1785/0120040243 [16] 张立. 层状介质中瑞利面波波场特征分析和反演方法研究[D]. 成都: 西南交通大学, 2009. [17] 于文福. 点源横向非均匀介质瑞利面波波场数值模拟研究[D]. 成都: 西南交通大学, 2017. [18] 陈伟,徐义贤,张煜. 瑞雷面波在弹性楔状体中的传播特性研究[J]. 工程地球物理学报,2009,6(2): 143-149. doi: 10.3969/j.issn.1672-7940.2009.02.003CHEN Wei, XU Yixian, ZHANG Yu. Study of transmission characteristics of Rayleigh waves on an elastic wedge[J]. Chinese Journal of Engineering Geophysics, 2009, 6(2): 143-149. doi: 10.3969/j.issn.1672-7940.2009.02.003 [19] 苏传行,刘争平,谭屹,等. 二维起伏地形瑞利面波速度频散特性数值模拟研究[J]. 地球物理学进展,2019,34(2): 810-816. doi: 10.6038/pg2019BB0384SU Chuanhang, LIU Zhengping, TAN Yi, et al. Study to the dispersion characteristics of Rayleigh waves trvalling along topographies with numerical modeling[J]. Progress in Geophysics, 2019, 34(2): 810-816. doi: 10.6038/pg2019BB0384 [20] 张碧星,鲁来玉,鲍光淑. 瑞利波勘探中“之”字形频散曲线研究[J]. 地球物理学报,2002,45(2): 263-274. doi: 10.3321/j.issn:0001-5733.2002.02.013ZHANG Bixing, LU Laiyu, BAO Guangshu. A study on zigzag dispersion curves in Rayleigh wave exploration[J]. Chinese Journal of Geophysics, 2002, 45(2): 263-274. doi: 10.3321/j.issn:0001-5733.2002.02.013 [21] 张凯,刘争平,熊自英,等. 声电效应在瑞利面波频散特性物理模拟中的应用[J]. 地球物理学进展,2015,30(6): 2955-2962. doi: 10.6038/pg20150665ZHANG Kai, LIU Zhengping, XIONG Ziying, et al. Application of seismoelectric effect to the experiment measurememt of Rayleigh surface waves[J]. Progress in Geophysics, 2015, 30(6): 2955-2962. doi: 10.6038/pg20150665 [22] 张立,刘争平. 水平层状介质中基阶瑞利面波椭圆极化特征数值分析与研究[J]. 地球物理学报,2013,56(5): 1686-1695. doi: 10.6038/cjg20130526ZHANG Li, LIU Zhengping. A study of the elliptic polarization characteristics of fundamental mode Rayleigh wave based on numerical simulation[J]. Chinese Journal of Geophysics, 2013, 56(5): 1686-1695. doi: 10.6038/cjg20130526 [23] 刘争平,张立,阎胜. 局部非均质体的瑞利面波动力学响应数值模拟研究[J]. 工程地球物理学报,2009,6(增刊1): 1-3.LIU Zhengping, ZHANG Li, YAN Sheng. The dynamic feature analysis of Rayleigh surface wave on a local inhomogeneity using the numerical method[J]. Chinese Journal of Engineering Geophysics, 2009, 6(S1): 1-3. [24] 郭明珠,谢礼立,凌贤长. 弹性介质面波地脉动单点谱比法研究[J]. 岩土工程学报,2004,26(4): 450-453. doi: 10.3321/j.issn:1000-4548.2004.04.004GUO Mingzhu, XIE Lili, LING Xianzhang. Research on spectral ratio of horizontal to vertical component for elastic model and surface microtremors[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(4): 450-453. doi: 10.3321/j.issn:1000-4548.2004.04.004 [25] 白永学,漆泰岳,吴占瑞. 砂卵石地层盾构施工引发的滞后地面塌陷机理[J]. 土木建筑与环境工程,2013,35(1): 12-19.BAI Yongxue, QI Taiyue, WU Zhanrui. Lagged surface collapse mechanism of shield construction in sandy pebble stratum[J]. Journal of Civil, Architectural and Environmental Engineering, 2013, 35(1): 12-19. [26] 罗松,张浩然. 成都富水砂卵石地层盾构施工滞后沉降防控措施探讨[J]. 隧道建设,2010,30(3): 317-319,335.LUO Song, ZHANG Haoran. Discussion on prevention and control of delayed settlement induced by shield tunneling in water-rich sandy cobble stratum in Chengdu[J]. Tunnel Construction, 2010, 30(3): 317-319,335. [27] 李哲生. 瞬态多道瑞利波勘探技术在岩土工程勘察中的应用[J]. 工程勘察,1996,24(3): 33-34.LI Zhesheng. Application of transient multichannel Rayleigh wave exploration technique in geotechnical engineering investigation[J]. Geotichnical Investigation and Surveying, 1996, 24(3): 33-34. [28] FÄH D, KIND F, GIARDINI D. Inversion of local S-wave velocity structures from average H/V ratios, and their use for the estimation of site-effects[J]. Journal of Seismology, 2003, 7(4): 449-467. doi: 10.1023/B:JOSE.0000005712.86058.42 [29] 邓瑞. 横向非均匀介质多模态瑞利面波波场数值模拟研究[D]. 成都: 西南交通大学, 2020. [30] 李梦. 水平层状介质中多模式瑞利面波波场有限元数值模拟研究[D]. 成都: 西南交通大学, 2020. 期刊类型引用(1)
1. 刘乃飞，周浩，宋战平，刘廉柏超，陶磊. 盾构隧道围岩特性实时获取技术研究进展. 长沙理工大学学报(自然科学版). 2024(05): 86-103+135 . 
百度学术其他类型引用(1)
-
下载:
下载:
百度学术
点击查看大图
计量
- 文章访问数: 349
- HTML全文浏览量: 203
- PDF下载量: 27
- 被引次数: 2