基于偏振分析方法的微动勘察技术应用与研究.pdf

收稿日期:2023-02-01;修回日期:2023-08-28第一作者简介:李鹏(1983),男,湖北襄阳人,2014 年毕业于中国地质大学(武汉),地球探测与信息技术专业,博士,高级工程师,现从事工程物探工作。E-mail:187196605 。引用格式:李鹏,张建清,陈爽爽.基于偏振分析方法的微动勘察技术应用与研究J.隧道建设(中英文),2023,43(9):1565.LI Peng,ZHANG Jianqing,CHEN Shuangshuang.Application and research of polarization analysis method-based microtremor survey methodJ.Tunnel Construction,2023,43(9):1565.基于偏振分析方法的微动勘察技术应用与研究李 鹏1,2,张建清1,2,陈爽爽1,2(1.长江地球物理探测(武汉)有限公司,湖北 武汉 430010;2.水利工程健康诊断技术创新中心,湖北 武汉 430010)摘要:为解决线形微动方法在地下工程勘察数据采集和处理中所面临的噪声信号源偏离测线方向的问题,采用偏振分析方法,结合理论研究和现场实践,提出一种城市强噪声环境下的微动勘察技术,可实现数据的有效采集和准确处理。该技术已在深圳市罗田水库铁岗水库输水隧洞工程中得到应用,应用结果表明:1)偏振分析方法可用于城市强噪声环境下线形微动数据方位角的计算,得出不同优势来源噪声信号的方位;2)当噪声源偏线角度较大时,可对数据进行分频偏振分析,求取方位角信息并用于微动数据处理,进而得到测区更加准确的视横波速度剖面。关键词:输水隧洞;微动勘察技术;数据校正;偏振分析DOI:10.3973/j.issn.2096-4498.2023.09.014文章编号:2096-4498(2023)09-1565-08中图分类号:U 45 文献标志码:A开放科学(资源服务)标识码(OSID):A Ap pp pl li ic ca at ti io on n a an nd d R Re es se ea ar rc ch h o of f P Po ol la ar ri iz za at ti io on n A An na al ly ys si is s MMe et th ho od d-B Ba as se ed d MMi ic cr ro ot tr re em mo or r S Su ur rv ve ey y MMe et th ho od dLI Peng1,2,ZHANG Jianqing1,2,CHEN Shuangshuang1,2(1.Changjiang Geophysical Exploration&Testing Co.,Ltd.(WUHAN),Wuhan 430010,Hubei,China;2.Technological Innovation Center for Health Diagnosis of Water Conservancy Project,Wuhan 430010,Hubei,China)A Ab bs st tr ra ac ct t:The noise signal source deviates from the survey line direction during data collection and processing of underground engineering survey using the linear microtremor method.Hence,the polarization analysis method is adopted,and a microtremor survey technology in urban strong noise environment is proposed based on theoretical research and on-site practice,allowing for effective data collection and accurate processing.The proposed method is applied to the water conveyance tunnel project from Luotian reservoir to Tiegang reservoir in Shenzhen,China.The outcomes indicate the following:(1)The polarization analysis method can calculate the azimuth angle of linear microtremor data in strong urban noise environments,obtaining the azimuth of noise signals from different sources.(2)When the deviation angle of the noise source is large,frequency division polarization analysis can be performed on the data to obtain azimuth information for micro motion data processing,resulting in a more accurate apparent shear wave velocity profile in the measurement area.This technology improves the quality of survey results by accurately correcting the orientation of survey data.K Ke ey yw wo or rd ds s:water conveyance tunnel;microtremor survey technology;data correction;polarization analysis0 引言随着城市化进程的加速,城市内的工程物探勘察工作变得越来越重要和频繁。城市中存在大量的地下设施和基础设施,如建筑物、道路、管线、桥梁等,因此需要进行详细的地下勘探和调查,确保工程施工安全。为了获得准确的地质信息,在城市工程物探勘察中通常会使用多种地球物理探测技术和方法,但是在城市电磁、噪声等环境的干扰下,电磁类方法、电法、地震反隧道建设(中英文)第 43 卷射法等常规地球物理方法使用受限1。近几年,微动勘察技术应用越来越广泛,该技术可以利用城市中广泛存在的振动噪声进行数据处理和成像,常用于水利2-5、交通、市政等行业的地质勘探6,其原理是利用自然界存在的震动,通过专业的仪器设备采集震动信号,然后采用专业的数据软件处理得到地层的面波速度结构,从而为规划、设计和地质等专业提供数据支撑。微动法的台阵布设方式主要有嵌套三角形、圆形、L 形和线形等,其中,线形的排列方式采集效率最高、布设方式最简单7。但目前微动法所采用的线形采集方式存在不足之处。微动采集的是天然源的噪声,在自然条件下,天然源的噪声来自四面八方,但在某些情况下,由于人为的干扰(汽车、施工等)或突发的自然灾害(地震、火山爆发等),天然源的噪声具有优势来源方向。因此,理论上来说采用线形排列时排列的方向应指向噪声的优势来源方向,以保证采集的资料信噪比最高、反演的速度最真实。但是在实际施工过程中,很难将排列指向噪声的优势来源方向,原因是:1)天然源的噪声很微弱,施工现场很难有条件进行噪声源调查;2)线形排列的方向一般与设计好的工作测线方向重叠,不允许偏离设计好的工作测线。当线形排列的方向不指向天然源的优势来源方向时,采集的数据经过处理以后,得到的横波速度值与真实的横波速度值不符,从而影响资料解释的精度8-9。针对微动所存在的噪声源偏离排列的问题,已有学者进行了相关研究,研究显示:空间自相关法可以消 除 噪 声 源 方 向 的 影 响;F-K(frequency-wavenumber)法可以推断波的到来方向,但一般要求各传感器在研究区内尽量呈平面展布10-11。Louie12在2001 年提出了微动折射法用于线形排列方式,但是需要假设面波的传播方向沿测线方向入射。Park 等13在 2008 年提出了路边噪声法,可用于偏线噪声数据的分析,通过在 180范围内扫描,得到给定频率不同方位角的相速度信息。Cheng 等14针对定向噪声源数据提 出 了 MAPS(multi-channel analysis of passive surface waves)处理方法。Liu 等15提出在原有的线形排列中加入几个垂直测线的检波器,形成一个伪线形排列,增加方位角的覆盖范围,求取噪声源方位角。以上研究均是针对线形排列的沿线或偏线的微动,采用单分量数据进行处理,而本文采用三分量数据进行处理,使用三分量数据可以提供更多信息来解决噪声源偏离排列的问题;同时,通过考虑 3 个方向上的振动信号,可以增加噪声源方位角的准确度和稳定性。相比于单分量数据,三分量数据能够提供更全面的观测特征,从而更有效地分析和定位噪声源。本文依托深圳市罗田水库铁岗水库输水隧洞工程,结合理论研究和现场实践,研究采用偏振分析方法解决线形微动技术存在的噪声源偏线问题,对线形微动排列数据进行方位角校正,直接通过线形排列的三分量数据进行数据的分析和处理,无需其他限定条件,可以提高数据采集和数据处理的效率,同时提高勘察成果质量。1 问题提出1.1 工程背景深圳市罗田水库铁岗水库输水隧洞工程(简称罗铁工程)线路全长约 21.6 km,隧洞埋深 50190 m,输水线路从罗田水库进水口起,近南北向接至铁岗水库出水口16。物探工作采用微动法,测线覆盖整条输水线路(含局部比选线路),微动测线总长约 53 km。微动测线布置示意如图 1 所示。图 1 罗铁工程微动测线布置示意图Fig.1 Layout of microtremor survey line of Luotie project1.2 线路特点隧洞沿线地形地质条件复杂,地形地貌为低丘陵、河流冲积平原、丘陵山地。丘陵山地地面高程一般为 100350 m,地表被林地覆盖;冲积平原段平原地面高程为 430 m,主要为城区,部分为农田。隧洞下穿深圳市宝安区繁华街区、京广深客运专线铁路、南光龙大高速公路、民生路、洲石公路以及茅洲河、溪沟和小水库等地表水体17。工程区第四系地层分布广泛,基岩岩性复杂,变质岩、岩浆岩、沉积岩均有分布18。测区地质构造发育,地层岩层褶皱剧烈、接触关系复杂,断层多、规模较大,破碎带普遍存在岩体硅化的现象。隧洞穿越 F1111、F4291、F1121、F3411、F3341 等断裂,断裂影响带宽数米至数十米,工程特性差。断层 F3341 具有活动性19。复杂的地形地质条件和构造、多样化的地层岩性分布、强烈的城市环境干扰为勘察方法的选择和勘察工作的布置提出了较大的挑战20。2 方法与原理2.1 微动基本原理从微动数据中提取频散曲线的方法主要有频率波数法(F-K 法)和空间自相关法等。空间自相关法包括 SPAC(spatial autocorrelation)法和 ESPAC(extended spatial autocorrelation)法。SPAC 法仅适用于圆形台阵6651第 9 期李 鹏,等:基于偏振分析方法的微动勘察技术应用与研究观测,如图 2 所示,位于圆心的接收点为中心点,其余接收点等角度分布于圆周上;而 ESPAC 法结合了SPAC 法和 F-K 法的优点,适用于任意形状的台阵(如图 3 所示)。图 2 圆形台阵Fig.2 Circular array图 3 非圆形台阵Fig.3 Non-circular arraySPAC 方法由 Aki21在 1957 年提出,根据空间自相关系数和零阶第一类贝塞尔函数的关系计算面波相速度。其基本原理是:对于圆形台阵,其中心点 C(0,0)与圆周上的任一接收点 X(r,)接收的角频率为 w的面波信号分别为 u(0,0,w,t)和 u(r,w,t),则其空间自相关函数(r,w)=u(0,0,w,t)u(r,w,t)。(1)式中:r 为接收点与中心点之间的距离;为接收点的方位角;t 为时间。空间自相关系数 p(r,w)为空间自相关函数在所有方向上的平均,即(r,w)=1220(r,w)d。(2)式(2)的积分结果可表示为(r,w)=J0wrv(w)=J02frv(f)。(3)式中:J0(x)为第 1 类零阶贝塞尔函数;v(w)为面波相速度;f 为频率。从式(3)中可见,空间自相关系数是面波相速度和频率的函数,通过拟合可以导出面波相速度。由于 SPAC 法要求采用圆形台阵,在实际工作中受野外地形地表条件的限制有时很不方便。Ling 和Okada 在 1993 年提出了 ESPAC 方法,成功用于非规则台阵。根据式(3),对于常量 f,有 f=f0,v=v0,则(f0,r)=J0(2f0v0)。(4)式(4)等号左边为随 r 变化的单值函数。在(f0,r)随 r 变化的情况下,可以通过最小二乘法求得最优的 Bessel 函数。ESPAC 法就是基于这一原理,可以通过在较大范围内增加r,(f0,r)数据以提高解的精度22。2.2 偏振分析设 Xj=xi(j=1,2,3;i=1,N)为三分量地震数据。其中,N 为采样点数,1、2、3 分别代表 Z、N 和 E 3 个方向的分量23。构建协方差矩阵 S,即S=Var X1Cov X1,X2Cov X1,X3Cov X1,X2Var X2Cov X2,X3Cov X1,X3Cov X2,X3Var X3。(5)式中 Var 和 Cov 分别代表方差和协方差。求取矩阵 S 的特征值 j(j=1,2,3)和对应的特征向量 uj(j=1,2,3),其中,最大特征值对应的特征向量为 u1,表征了波的传播方向,同时也表征了数据在各个方向的能量分布。微动主要是针对瑞雷面波进行分析,瑞雷面波沿自由表面传播,介质质点位移由水平位移和垂直位移叠加做逆时针的椭圆运动,水平方向的质点振动与波传播方向在一条直线上,而方位角的值由水平方向的分量计算得出。因此,优势震源的方位角 P 可由式(6)计算得出。P=tan-1u12sign(u11)u13sign(u11)。(6)式中 u1j(j=1,2,3)为特征向量 u13 个方向的余弦。由式(6)可知,方位角的计算使用了特征向量 u13个方向的余弦,其中,垂直方向的余弦 u11决定了方位角的象限(相当于在线形测线的小桩号或大桩号产生的振动),而方位角的值由水平方向的余弦 u12和 u13计算得出。2.3 强噪数据校正在城市环境中进行长距离线路勘探时,为了提高数据采集的效率,常常采用线形的数据采集方式。城市中的噪声来源丰富,汽车行驶、工地打桩、行人行走等均可产生振动,特别是沿道路布置的测线,会受到汽7651隧道建设(中英文)第 43 卷车等强噪声的干扰。采用线形排列时测线不可能一直指向优势噪声源方向,如果采用常规的微动处理方法对线形排列采集的数据进行处理,则得到的视横波速度剖面存在较大误差;如果能够了解每条测线优势噪声的来源方向,对强噪数据进行校正,则可保证数据处理的准确性。通过带通滤波将原始数据提取为一系列窄频带,对各频带数据进行偏振分析,可求得每个频带数据对应的优势方位角 Pf,得到各频率成分的方位校正因子 Af。Af=cos Pfdx。(7)式中 dx 为道间距。采用方位校正因子对各台阵之间的距离进行校正,再采用相应的频散谱分析方法提取频散曲线。2.4 理论效果分析为了验证本文数据校正方法的有效性,设计了如图 4 所示的城市环境下的路边微动观测系统。噪声源分为远距离低频噪声、近距离高频噪声和马路噪声,检波器沿 y 轴线形布置,远距离低频噪声来源于 360方向,近距离高频噪声来源于 4560方向,马路噪声来源于近 90方向。图 4 微动震源分布图Fig.4 Distribution of microtremor sources选取圆心的检波器采集的三分量地震数据进行分析,理论数据的 Z 分量、X 分量和 Y 分量如图 5 所示。其中,Z 分量指向垂直方向,X 分量指向 E 方向,Y 分量指向 N 方向。采用窄带滤波器将理论数据分解为低频信号、高频信号,并根据路边噪声的信号特点分解出路边噪声信号,然后采用本文所述的偏振分析方法分别对低频、高频和路边噪声信号进行偏振分析。(a)Z 分量(b)X 分量(c)Y 分量图 5 理论数据原始记录Fig.5 Raw record of theoretical data理论数据的偏振分析结果如图 6 所示。图中颜色代表归一化后的方位角分布密度,其中,暖色代表分布密集,冷色代表分布稀疏,半径代表微动数据的采集时间。低频信号来源于自然界的随机噪声,这种信号分布广泛、传播距离远,偏振分析的结果也验证了其分布的随机性;高频信号主要来源于城市人类活动造成的噪声,偏振分析显示其主要分布于 45偏向 60方向的扇形区域,与正演模拟数据一致;马路噪声主要来源于穿行的车辆,由于噪声源离检波器距离较近,因此其频率较高、能量较强,偏振分析显示其主要分布于 0偏向 10的扇形区域,与正演模拟数据一致。以上理论数据分析证明了采用偏振分析的方法可以获取不同类型噪声源的来源方向,从而可以准确计算出不同类型噪声的方位角。在实际应用中,偏振分析方法可以作为微动数据处理流程的一个环节,为微动勘探提供准确的结果,为研究城市线形微动方法提供参考和依据。3 应用效果分析3.1 沿线噪声数据分析应用地点位于罗铁工程西线比选线路,数据采集采用 0.2150 Hz 频带宽度的三分量检波器。为了保证数据采集的质量,在检波器上覆盖沙袋进行压重耦合,采用路边噪声的观测系统,测线平行于马路布置于路边,工作布置如图 7 所示(红线为测线)。检波器间距为 10 m,检波器的 Y 分量指向正北方向,马路走向为北微东约 10。实际数据原始记录如图 8 所示。偏振分析结果如图 9 所示。测线附近的机械作业产生高频的噪声(如图 9(a)所示),其主要分布于 30偏向 40方向的扇形区域;马路噪声主要分布于北微偏东方向的扇形区域内(如图 9(b)所示)。马路噪声方位角统计如图 10 所示。由于车辆在马路上双向行驶,因此马路噪声的方位角主要集中于 0、180和 360 3 个区域,马路噪声相对于测线的偏线角度为 510,8651第 9 期李 鹏,等:基于偏振分析方法的微动勘察技术应用与研究校正前后的频散曲线如图 11 所示。图中的彩色代表频散能量。由图可知,校正前后的相速度差较小,大部分在 10 m/s 以内。这充分印证了在城市环境中进行地下工程勘察时线形的排列方式可以满足勘察的需要。(a)低频噪声(b)高频噪声(c)马路噪声图 6 理论数据的偏振分析结果Fig.6 Polarization analysis results of theoretical data上为正北方向。图 7 工作布置图Fig.7 Work layout(a)Z 分量(b)X 分量(c)Y 分量图 8 实际数据原始记录Fig.8 Raw record of actual data(a)高频噪声(b)马路噪声图 9 实际数据偏振分析结果Fig.9 Polarization analysis results of actual data9651隧道建设(中英文)第 43 卷图 10 马路噪声方位角统计图Fig.10 Polarization analysis statistical chart of road noise图 11 校正前后的频散曲线Fig.11 Dispersion curve before and after correction3.2 偏线噪声数据分析部分测线横穿马路或测线周边有明显噪声源时,噪声相对于测线的偏线角度较大(如图 12 所示),对不同频率的噪声数据进行偏振分析。图中颜色表示归一化的方位角分布密度,冷色代表分布稀疏,暖色代表分布密集。结果显示:19 Hz 的噪声源相对于测线的偏线角度较大(1558),1020 Hz 的噪声源相对于测线的偏线角度较小(07)。将不同频率的方位信息应用于噪声数据并进行频散能量成像,结果如图 13 所示。图中分别示出了校正前后的频散曲线,其中,19 Hz 频段的数据由于偏线角度较大,校正后的相速度明显小于校正前的相速度;大于 9 Hz 频段的数据校正前后相速度基本无差异。校正前后的反演成像如图 14 所示。波速值是工程中用于刻画地下介质特性的重要参数,可以提供关于地质结构、岩性的重要信息。校正前的反演成像图显示整体剖面具有较高的波速,而校正后的反演成像图显示从浅区域到深区域的波速值有所降低,这意味着在校正过程中对波速进行了修正。为了验证校正方法的效果,在桩号 320 m 处布置了钻孔 Z1 进行验证。钻孔 Z1 揭露全风化层、强风化层和弱风化层的深度分别为 24、27.8、41.5 m,同时在深度约 70 m 发育一倾角约为 65的断层。校正后的反演成果与钻孔成果基本一致,而校正前的反演成果在风化层的划分上与钻孔成果存在较大差异,如果不对数据进行方位角校正,将无法准确刻画出层位信息。图 12 不同频率偏振分析结果Fig.12 Polarization analysis results at different frequencies图 13 不同频率校正前后的频散曲线Fig.13 Dispersion curve before and after correction at different frequencies0751第 9 期李 鹏,等:基于偏振分析方法的微动勘察技术应用与研究(a)校正前(b)校正后图 14 校正前后的反演成像图Fig.14 Inversion profile before and after corrction4 结论与展望本文针对城市环境下地下工程勘察所面临的噪声干扰大、方法应用受限等问题,提出了基于偏振分析的微动勘察技术,实现了对城市环境下偏线微动数据的方位角校正,得出如下结论:1)偏振分析方法可用于城市强噪声环境下线形微动数据方位角的计算,准确得出不同来源噪声信号的方位。2)在城市环境中使用线形微动方法是可行的,特别是沿马路布置的微动观测系统,当主要噪声来源于测线沿线时,频散能量拾取误差较小。3)当噪声源偏线角度较大时,可对数据进行分频偏振分析,求取方位角信息并用于微动数据处理,从而得到测区更加准确的视横波速度剖面。基于偏振分析的城市强噪环境下的线形微动数据方位角分析方法的成功应用,可为城市地下工程的勘察提供科学和有效的解决方案,下一步可通过进一步的数据分析,不断完善数据处理流程,提高数据信噪比和成果质量。参考文献(R Re ef fe er re en nc ce es s):1 刘宏岳,黄佳坤,孙智勇,等.微动探测方法在城市地铁盾构施工“孤石”探测中的应用:以福州地铁 1 号线为例J.隧道建设,2016,36(12):1500.LIU 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.1J.Tunnel Construction,2016,36(12):1500.2 钮新强,张传健.复杂地质条件下跨流域调水超长深埋隧洞建设需研究的关键技术问题J.隧道建设(中英文),2019,39(4):523.NIU Xinqiang,ZHANG Chuanjian.Some key technical issues on construction of ultra-long deep-buried water conveyance tunnel under complex geological conditionsJ.Tunnel Construction,2019,39(4):523.3 邓娜,张正军.深圳市地下交通发展现状与规划思考J.隧道建设(中英文),2022,42(12):1985.DENG Na,ZHANG Zhengjun.Current status and planning of underground transportation in Shenzhen,China J.Tunnel Construction,2022,42(12):1985.4 洪开荣.我国隧道及地下工程近两年的发展与展望J.隧道建设,2017,37(2):123.HONG Kairong.Development and prospects of tunnels and underground works in China in recent two yearsJ.Tunnel Construction,2017,37(2):123.5 杜峰.长距离盾构法输水隧道设计要点及改进建议J.隧道建设,2011,31(4):484.DU Feng.Key design points and optimization recommendations for long-distance water-supplying tunnels bored by shield machinesJ.Tunnel Construction,2011,31(4):484.6 张建清,熊永红,李鹏.微动勘探技术在城市轨道交通勘察中的应用J.人民长江,2016,47(1):39.ZHANG Jianqing,XIONG Yonghong,LI Peng.Application of micro tremor technology in urban rail transit surveyJ.Yangtze River,2016,47(1):39.7 程建设,李鹏.微动勘探技术在水库大坝隐患探测中的应用J.人民长江,2017,48(3):57.CHENG Jianshe,LI Peng.Application of micro tremor exploration technology in detection of reservoir and dam latent hazardsJ.Yangtze River,2017,48(3):57.8 李洋,党如姣,朱培民,等.花岗岩风化残留体地震法1751隧道建设(中英文)第 43 卷探测应用研究J.隧道建设,2015,35(5):435.LI Yang,DANG Rujiao,ZHU Peimin,et al.Application of seismic method in detection of weathered-granite bouldersJ.Tunnel Construction,2015,35(5):435.9 陈爽爽,李鹏,张智,等.线性台阵被动源面波法在地下空间探测中的应用J.人民长江,2022,53(2):77.CHEN Shuangshuang,LI Peng,ZHANG Zhi,et al.Application of linear-array passive source surface wave method in underground space exploration J.Yangtze River,2022,53(2):77.10 庞景尹.浅地表被动源面波成像中数据筛选方法研究D.北京:中国地质大学,2019.PANG Jingyin.Data selection for imaging near-surface passive surface wavesD.Beijing:China University of Geosciences,2019.11 王振东.面波勘探技术要点与最新进展J.物探与化探,2006,30(1):1.WANG Zhendong.Essentials and recent advances of the surface wave exploration technique J.Geophysical&Geochemical Exploration,2006,30(1):1.12LOUIE J N.Faster,better:Shear-wave velocity to 100 meters depth from refraction microtremor arrays J.Bulletin of Seismological Society of America,2001,91(2):347.13 PARK C B,MILLER R D.Roadside passive multichannel analysis of surface waves(MASW)J.Journal of Environment and Engineering Geophysics,2008,13(1):1.14CHENG Feng,XIA Jianghai,LUO Yinhe,et al.Multichannel analysis of passive surface waves based on cross correlationsJ.Geophysics,2016,81(5):57.15LIU Ya,XIA Jianghai,CHENG Feng.Pseudo-linear-array analysis of passive surface waves based on beam formingJ.Geophysical Journal International,2020,221(1):640.16 陈令明,贾宁霄,冯敏,等.罗田水库铁岗水库输水隧洞工程 BIM 正向设计及应用研究J.土木建筑工程信息技术,2022,14(6):121.CHEN Lingming,JIA Ningxiao,FENG Min,et al.Research on BIM forward design and application of Luotian reservoir-Tiegang reservoir water conveyance tunnel engineeringJ.Journal of Information Technology in Civil Engineering and Architecture,2022,14(6):121.17 李粮纲,徐玉胜,江辉煌,等.深圳地区地质环境特征与地质灾害防治J.安全与环境工程,2007(4):28.LI Lianggang,XU Yusheng,JIANG Huihuang,et al.Characters of geological environment and the precaution of geological hazards in Shenzhen area J.Safety and Environmental Engineering,2007(4):28.18 王锐,胡坤生,张延仓,等.罗田-铁岗输水隧洞围岩蚀变花岗岩特性J.人民长江,2021,52(增刊 2):79.WANG Rui,HU Kunsheng,ZHANG Yancang,et al.Characteristics of erosion granite in surrounding rock of Luotian-Tiegang water conveyance tunnel J.Yangtze River,2021,52(S2):79.19 孙杰,贾建业,詹文欢,等.深圳断裂带构造活动性分析J.地球科学进展,2007(3):234.SUN Jie,JIA Jianye,ZHAN Wenhuan,et al.A study of the tectonic activity of Shenzhen fracture zoneJ.Advances in Earth Science,2007(3):234.20 王栩,王志辉,陈昌昕,等.城市地下空间地球物理探测技 术 与 应 用 J.地 球 物 理 学 进 展,2021,36(5):2204.WANG Xu,WANG Zhihui,CHEN Changxin,et al.Geophysical exploration and application for urban underground spaceJ.Progress in Geophysics,2021,36(5):2204.21 AKI Keiiti.Space and time spectra of stationary stochastic waves,with special reference to microtremorsJ.Bulletin of the Earthquake Research Institute,1957,35:415.22 赵东.被动源面波勘探方法与应用J.物探与化探,2010,34(6):759.ZHAO Dong.Passive surface waves:Methods and applicationsJ.Geophysical&Geochemical Exploration,2010,34(6):759.23 JURKEVICS A,贺冬梅,钟南才.三分向台阵数据的偏振分析J.世界地震译丛,2002(3):59.JURKEVICS A,HE Dongmei,ZHONG Nancai.Polarization analysis of three-component array data J.Translated World Seismology,2002(3):59.2751