地震背景噪声成像研究综述.pdf

SEISMOLOGICAL AND GEOMAGNETICOBSERVATION AND RESEARCH第 44 卷 第 2 期2023 年 4 月Vol.44 No.2Apr.2023地 震 地 磁 观 测 与 研 究doi:10.3969/j.issn.1003-3246.2023.02.003第一作者简介：郭瑛霞（1990），女，助理工程师，毕业于中国地震局兰州地震研究所，获硕士学位，主要从事地震活动性和数字地震学研究工作。E-mail：基金项目：青海格尔木青藏高原内部地球动力学野外科学观测研究站专项；2023 年度震情跟踪定向工作任务（项目编号：2023010114）；江苏省地震局青年科学基金（项目编号：202109）本文收到日期：2022-03-08地震背景噪声成像研究综述郭瑛霞1），2）张丽峰1），2）胡维云1），2）黄 浩1），2）李振凯2）1）中国西宁 810001 青海格尔木青藏高原内部地球动力学野外科学观测研究站 2）中国西宁 810001 青海省地震局摘要 地震背景噪声成像方法已成为 21 世纪地震学的伟大突破之一，其原理是，对 2 个台站记录的连续背景噪声信号进行互相关计算，得到台站间的格林函数，利用经验格林函数得到面波速度函数，获取面波频散特性，用地震层析成像方法对面波速度进行反演，得到地球内部的速度结构。其应用近年来日益广泛，涉及各向同性和各向异性的速度结构成像、大地震前后速度结构变化监测、体波联合成像、衰减结构成像、地震定位精度提高、噪声源分布和物理起源探究、强地面运动评估等。其优点是不需要等待天然地震或使用对环境构成威胁的人工爆破，所有台阵都可以当作源，拓宽了频带范围，可以获得较多短周期频散数据，提高反演分辨率。关键词 背景噪声；互相关函数；层析成像0 引言野外布设台阵越来越多，不仅可以记录到地震信号，也能记录到丰富的连续噪声信息，这些信息来自地球震动、固体潮汐、大气气压变化、海浪拍击海岸等。起初这些连续的背景噪声记录被视为无用的干扰信号，在进行天然地震研究时，通过各种方法手段去除噪声而凸显地震信号，后来经过科学家地不断研究，发现背景噪声记录也能用于地下速度结构研究（Duvall et al，1993），且具有较高的分辨率。早在 1957 年，Aki（1957）就通过布设合理的台站，使用背景噪声记录获取了浅层地壳中面波频散信息，提出可以采用空间自相关方法研究地表浅层介质；1968 年，Claerbout（1968）提出台站间的反射响应，通过在三维各向异性介质中 2 个台站的投射波场计算得到互相关函数，该思路在太阳地震学中得到成功应用。发展到 1993 年，Duvall 等（1993）通过互相关计算太阳表面的噪声场，提出声波时距曲线，并得到太阳外层的三维流速度结构，这一研究成果奠定了背景噪声技术发展的里程碑，越来越多的研究和学术领域，包括超声波、水底声波场，均发现从两点间的噪声信号提取的互相关函数与格林函数的相似度较高。2001 年，Lobkis 等（2001）在声学研究领域，对测量得到的热起伏噪声进行互相关计算，得到两点间的格林函数，这一开创性的研究成果被运用到地震学中。基于之前的理论指导，一种新的地震学方法诞生并逐渐发展起来，这种方法已证实可用于研究壳幔速度结构，19第 2 期即背景噪声层析成像。2003 年，Campillo 等（2003）发现，在地震尾波记录中也能提取格林面波函数；2004 年，Shapiro 等（2004）通过台站对地震噪声垂向分量信息进行互相关计算，得到地球结构信息。2005 年，Shapiro 等（2005）用此方法得到美国南加州周期为 7.5 s和 15 s 的瑞利面波群速度分布图。之后地震背景噪声成像方法得到飞速发展，应用范围小到地方、大到全球尺度，研究内容包括了瑞利面波相速度和群速度、勒夫面波群速度和相速度，研究周期可以是几秒或者上百秒，对应深度从地球近地表到地壳浅层到整个地壳以及上地幔（Shapiro et al，2005；Yao et al，2006；房立华等，2009；Fang et al，2010，2015；郑定昌等，2014；Li et al，2016）。研究人员证实，层析成像结果与已知地质构造吻合较好，表明背景噪声成像方法是有效的（姚华建等，2004）。本文在此基础上，对其在地震学中采用的方法理论原理、应用效果和未来发展前景进行综述和展望。1 背景噪声成像方法理论基础1.1 互相关函数与格林函数的关系利用两点间记录的连续背景噪声信号进行互相关函数计算，可以得到两点间的格林函数（Weaver et al，2005）。格林函数在地震学中的物理意义是：把源当做一个接收点产生的记录，在另一接收点记录到的脉冲响应就是 2 个接收点之间的格林函数，此脉冲响应可以认为是由无限多个具有不同周期的连续正弦波列组成，始于震源事件发生，终止于尾波结束，包含 2 个接收点间地震波传播的全部信息，反映出地下介质结构的信息，互相关函数的导数就是经验格林函数（图 1），经验格林函数不是真实的格林函数，是真实格林函数中的面波部分。同时，在理论方面，促使不同领域的研究者来证实其基于不同物理模型的假设，相关假设有：模式均分假设（Weaver，2005）、时间反对称理论（Derode et al，2001，2003）、稳相近似理论（Snieder，2004）、互相关类型的表示定理（Wapenaar et al，2004），对互相关函数和经验格林函数的关系给出不同理论解释，以证实二者具有内在一致性。1.2 面波频散曲线测量方法面波具有频散特性，可从经验格林函数中提取相速度或群速度频散曲线，地震学家通常采用单台法或双台法来提取面波频散特性（Brune et al，1960）。20 世纪 60 年代初，Alexander 等（1963）把数值滤波技术用于面波频散测量；Pliant 等（1964）首次提出时变滤波法（TVF）来测量相速度，并得到双台相速度频散，此项技术应用奠定了面波频散测量发展的基础。到1969年，Landisman等（1969）提出移动窗分析法；紧接着，Dziewonski等（1969）提出多重滤波法；Herrin 等（1977）为了使提取的相速度或群速度频散曲线更加精确和稳定，提出相位匹配滤波法。图 1 用背景噪声计算互相关函数提取格林 函数示意图（引自 Yao et al，2006）黑色箭头代表散射波场，是震源产生的记录经过多次散射 形成的，其中噪声信号先后被 2 个台站接收，震相不发生变化，通过对 2 个台站进行互相关计算，得到两接收点之间的格林函数Fig.1 Schematic diagram of extracting Greens function by calculating cross-correlation function of background noise(From Yao et al,2006)背景噪声记录背景噪声记录地震检波器地震检波器互相关函数 U1(t)U2(t+)dt=C()郭瑛霞等：地震背景噪声成像研究综述地震地磁观测与研究44 卷 20基于时变滤波思想（徐果明等，2000），姚华建等（2004）提出基于图像分析的相速度频散曲线快速提取方法，是一种基于双台法提取面波频散特性的有效方法，其原理是：对 2 个台站的连续波形记录做零漂校正、去仪器响应，归一化处理后采用多重滤波法，得到双台间不同周期 T 的波对应的不同到时，对归一化双台记录做加可变宽移动时窗的窄带滤波处理，得到频率上基本同频的简谐波（姚华建等，2004）。与传统双台法等值线手工描点勾画方法相比，该方法只需在频散曲线附近点动鼠标，便可快速追踪到整条频散曲线，提高了面波频散特性提取精度和效率，不仅可以清晰展现频散特征，而且可以有效识别和分析频散曲线（图 2）。人机交互为处理大批量背景噪声数据资料提供了快速可靠的软件基础。1.3 三维 S 波速度反演理论以往面波成像一般都假设面波沿大圆路径传播，路径是固定的，不随频率变化，也不随结构变化。但是做体波层析成像时，由于结构差异比较大，研究人员必须考虑体波传播路径的射线追踪。在全球地震学中，面波反演时由于群速度或相速度变化都在 7%以内，这时面波沿大圆路径传播与所考虑面波走的实际路径接近，对反演结果几乎没什么影响，便不考虑面波射线路径的弯曲。但据前人研究结果（Huang et al，2010；李想等，2015）表明，当地下介质速度结构变化较大，且某一周期的面波相速度相对于平均相速度的异常大于 10%时，面波会偏离大圆路径传播，导致反演结果与假设面波沿大圆路径反演的结果出现偏差。使用面波频散测量数据进行层析成像或背景噪声层析成像通常采取两步法：基于射线追踪理论或二维有限频敏感核使用走时层析成像方法反演得到二维相速度或群速度分布图（Yao et al，2006）；在每个网格节点逐点反演面波频散数据，作为剪切波速度随深度的函数的一维剖面，然后用一维模型反演出三维模型（Yao et al，2008）。Fang 等（2015）提出基于射线追踪的面波频散直接反演方法，把所有双台路径频散曲线各个不同周期的走时信息合在一起，直接反演出三维剪切波速度结构，这种方法省去了传统两步法（Ritzwoller et al，2002；Yao et al，2006）中通过由频散曲线做走时成像得到 2D 相速度/群速度分布图这一步，并考虑了复杂介质中面波会偏离大圆路径传播，受到地下结构的控制，不同频率的面波走得路径不一样，用快速行进法（Rawlinson et al，2004）进行面波传播路径的射线追踪。图 2 相速度频散曲线测量示意（引自 Yao et al，2006）(a)对每一周期的双台窄带滤波结果做互相关，求取面波经验格林函数的振幅矩阵图像，黑色波形为归一化处理的经验格林函数，中心周期 20 s；(b)反映周期和相速度关系的振幅矩阵图像，P 处黑色实线为相速度频散曲线测量图Fig.2 Schematic of phase velocity dispersion curve measurement(From Yao et al,2006)1201401601802002205.04.54.03.53.02.510152025303540455055601015202530354045505560周期/s速度/（kms-1）走时/sPP*(b)(a)21第 2 期1.4 背景噪声源特征原理分析理论研究表明，背景噪声源分布的均匀性对互相关函数和经验格林函数的近似程度十分重要。由于地球内部结构的不均匀性，实际观测到的噪声分布并不均匀，所以通过合理布设台阵、改进数据处理方法、分析背景噪声特征等方法，可以进一步完善噪声源的均匀分布，进而提高反演结果的分辨率。任意 2 个台站间计算的互相关函数，其振幅取决于 2 个台站连线方向的能量流（Paul et al，2005），如果台站两侧的噪声源分布均匀，说明 2 个方向的能量流是相同的，则互相关函数的振幅和到时大小是完全相同的；如果台站两侧的噪声源分布不均匀，则会发现互相关函数的到时相同、振幅不对称，说明互相关函数正负分支能量流的来源方位不同（图3）。同时，互相关函数 2 个方向的能量流在分布时间上表现为一定的季节性变化特征，会影响频散曲线提取的稳定性。因此，有必要分析背景噪声源的分布特征和季节性变化特征。Stehly 等（2006）选取了美国东部、欧洲西部、加利福尼亚地震台网记录的垂直分量背景噪声数据，通过分析互相关函数因果信号和非因果信号的不对称性，发现不同频带的噪声起源是不同的：环境背景噪声在 5 10 s 和 10 20 s 周期之间存在 2 个明显的峰值，10 20 s 周期之间存在的峰值称为第一类地脉动，会随季节的变化而变化；5 10 s 周期之间存在的峰值称为第二类地脉动，不随季节变化。图 3 噪声源分布对所得互相关函数的影响（引自 Stehly et al，2006）(a)当 1、2 两侧的噪声源分布