陈玉鑫 唐明帅

摘要：利用地震台记录的背景噪声数据，进行台站对之间数据互相关以获得反映地下介质信息的格林函数，是目前地震学领域理论研究和实际应用方面的前沿和热点。首先介绍了噪声互相关应用的发展历史，之后对其在面波成像、体波信号提取与体波成像、监测地下介质波速变化与衰减、地震定位、噪声源分析与时移校正等多个方面的应用进行了阐述，最后针对目前背景噪声互相关应用的发展进行了分析与讨论。

关键词：背景噪声;格林函数;频散曲线;层析成像;介质变化

中图分类号：P315.3文献标识码：A文章编号：1000-0666（2021）04-0594-13

0引言

长久以来，地震学家们利用地震波走时或振幅信息对地球内部结构及其变化进行研究，这种方法通常采用台站记录中的天然地震信息或人工地震产生的地震波。一般天然地震多集中在构造运动活跃的板块边界或断裂带，这使得地震活动性低的区域无法得到充足可靠的地震资料。人工地震往往采用爆破等方式，不仅费用昂贵，还面临着破坏地层结构、污染环境等风险。以上两种方法对于研究地下介质信息有很大的局限性，而背景噪声提供了新的途径。背景噪声是指那些通过各种拾震器采集的、按照常规的地震数据处理方法难以识别出有效信号的、常常作为干扰剔除或压制的地震数据（徐义贤，罗银河，2015）。过去背景噪声被认为是影响震相提取的无用随机信号。随着认识的加深，人们发现台站对记录的噪声互相关结果，与反映台站间地下介质信息的格林函数存在联系。由于地震噪声是一种连续、稳定的自然源，不存在传统研究方法的缺陷，因而可以很好地应用于地球内部结构研究。只要有足够的台站以及合理的布设，就可以在限制条件很多的地区展开应用（孔祥艳，2019），因而受到越来越多科学家的关注，为揭示世界各地的构造变形和地下动力学机制做出了重要贡献。

本文回顾了背景噪声互相关方法在地震学领域的发展历史，并对其应用进行了介绍，包括噪声面波成像、提取体波信号及体波成像、波速变化和衰减结构监测、地震定位、噪声源分析、时移校正等，最后分析了现有背景噪声互相关应用存在的问题，对其发展前景进行了展望。

1背景噪声互相关应用发展历史

背景噪声互相关方法利用时刻都在产生的噪声信号获取地下介质的信息，它的发展历史可以追溯到20世纪50年代。Aki（1957）提出利用空间自相关方法（Spatial AutoCorrelation method ，SPAC）从“微震”信号中获得研究区域相速度分布。Claerbout（1968）证明了在水平层状介质中，自由地表接收的来自地下的地震记录自相关等价于该点自激自收的模拟记录。随后，这种思想就被推广到地表任意两个接收器之间计算互相关信号，并在太阳地震学领域得到成功实践。Duvall 等（1993）率先完成了利用太阳表面的噪声信息互相关提取时距曲线的突破，这种方法被称为声波日光成像法。后来，Weaver 和Lobkis（2001，2002，2003）通过实验发现接收电路中热噪声自相关的时间导数与直接获得的脉冲信号相同，因此又对比了铝块上两点记录的热噪声互相关结果与格林函数的差异，并基于模式均分假设对此进行了解释。此后，互相关函数和格林函数之间的关系得到了众多研究者关注，并基于不同假设，提出了时间反对称理论（Derode et al ，2003）、稳相近似理论（Snieder，2004）等进行了论证。

噪声互相关方法先后在太阳地震学、声学等领域的成功实践，促使地震学家开始寻找与“热噪声”类似的信号。首先應用于这种方法的是一种能够形成漫射波场的多次散射信号，类似于背景噪声，被称为地震尾波。Campillo和Paul（2003）对墨西哥两个台站记录到的101次远震的尾波进行互相关，发现提取到的波形具有Ray-leigh 波和Love波的极化特性和群速度。紧接着，Shapiro 等（2004，2005）证明了应用背景噪声提取地球结构的相干信息是可行的，并利用美国加州62个USArray观测台所记录的30d背景噪声数据提取了短周期Rayleigh 波频散进行层析成像，成像结果显示速度异常区域与地表地质构造有着良好的对应关系。

Yao等（2006）利用青藏高原东南部台网记录到的噪声数据首次计算了10～30s的经验格林函数（Empirical Green Functions ，EGFs），进而提取Rayleigh 波相速度频散进行成像。Bensen等（2007）对背景噪声提取面波频散的流程进行了系统性总结，并对常见的五种时域归一化方法进行比较，肯定了滑动绝对平均法去除地震信号的有效性和灵活性（图1）。Lin等（2008）基于美国西部250多个台站一年所记录的三分量噪声数据首次提取了Love波的频散信息，并指出利用背景噪声可以研究地壳内部各向异性。

地震背景噪声互相关研究在近些年得到了不断地完善，已经成为研究地球内部结构和动力学演变的重要手段，加深了人们对地球内部的认识。

2背景噪声互相关方法的应用

2.1噪声面波成像

地震背景噪声信号来源比较复杂，如海洋与陆地的相互作用、大气扰动以及人类活动等。这些噪声源主要分布在地表，因此，沿自由表面传播且能量衰减较慢的面波是噪声互相关函数（Noise Cross -correlation Function ，NCF）中最主要的信号，从浅层地表到上地幔深度，面波成像被广泛应用于地球内部结构研究。短周期密集台阵仪器布设间隔近（几米到几千米）、观测时间短（1～2个月），因而成本较低、时效性强、数据质量好，其高频数据可以用于对地表浅层地下结构进行面波成像（王仁涛等，2019;范兴利等，2020;杨润海等，2020;Li et al ;2019;Zhang et al ，2020）。Lin等（2013b）利用2011年布设在美国加利福尼亚Long Beach 的密集台阵三周噪声记录进行程函函数层析成像，获得了0.5～4Hz相速度分布。Li等（2016）基于布设在五大连池火山区的台阵数据观测了火山锥下的岩浆房，并在尾山火山下方7～13km处发现了显著的低速异常。不同周期面波速度随深度的变化对地下介质密度、P波速度及S波速度具有不同的敏感性，其中受S波速度影响最大，所以利用面波频散反演地壳和上地幔中群速度或相速度分布（Sabra et al ，2005;Bensen et al ，2008;付媛媛，高原，2016）进而研究三维S波速度变化是一种常见的研究方法。Li等（2009）通过四川台网记录的噪声数据获得了8～44s川西及藏东地区Rayleigh 波群速度分布，并基于Herrmann 的线性反演程序计算了研究区10、25、45km深度S波速度分布图像。Tang等（2013）采用环境噪声和远震双重面波层析成像法得到了华北克拉通6～143s Rayleigh 波相速度图像，并构建出该区域三维地壳和上地幔S波速度结构。速度图像清晰地显示了华北盆地与鄂尔多斯地块不同形式的岩石圈、地壳、上地幔结构，表明它们具有不同的热性质和（或）化学性质，自新生代以来经历了不同的地幔过程和演化历史。

噪声面波成像是探究地下不同深度结构的重要方法，因此被广泛应用于小尺度、区域尺度乃至全球尺度的研究（Nishida et al ，2009）。李玲利等（2020）对合肥市区城市下方地壳浅层速度结构进行了成像，并对S波速度变化与多组断裂的关系进行分析，揭示了城市下方隐伏断裂的展布特征以及地壳浅部的横向不均匀性和纵向成层性。曾求等（2020）计算得到了四川威远地区5km以内的三维S波速度模型，通过比较2016年威远地震震源位置和矿物开采区域与当地地质构造的关系探讨了地震发生机制。地震波在传播过程中速度及质点偏振方向等特性都会随波传播的方向而产生变化，即产生所谓的地震各向异性（张忠杰，许忠淮，2013），具体表现为两种（彭艳菊等，2007）：径向各向异性（反映了同一路径上SV和SH波速度结构的差异）和方位各向异性（代表了同一类型的波沿不同方向传播速度不同）。各向异性反映了地下矿物的特定结构和排列分布，对探究地表地质构造演化具有重要意义，一直以来是研究地下结构的热点（欧阳龙斌等，2015;杨志高等，2019;顾勤平等，2020;王怀富等，2020）。Xie等（2013）利用青藏高原东部及周边地区Rayleigh 波（T=8～65s）和Love波（T=8～44s）相速度分布及贝叶斯蒙特卡罗反演方法获得了地壳和上地幔SV和SH分布，讨论了该区域径向各向异性速度场。Yao等（2010）使用中国西南地区宽频台站记录的数据，结合噪声互相关方法（T=10～40s）和双台法（T=20～150s）提取Rayleigh 波相速度频散曲线，得到了青藏高原东南部岩石圈至150km深处的方位各向异性。Shen等（2016）基于中国以及周边国家2000多个台站的数据制作了覆盖全中国的各向同性和方位各向异性Rayleigh 波速度分布图像，为进一步探究地下结构提供了有益的S波参考模型。

传统噪声面波成像反演过程通常采用两步法，在获得研究区不同周期二维面波速度分布的基础上利用线性或非线性反演方法得到S波速度随深度的变化（房立华，2009）。这种反演过程在计算速度和成像分辨率等方面有待提升，因此，不断有新方法和新思路被提出以弥补传统方法的不足。Fang等（2015）提出了基于射线追踪的面波直接反演方法，这种方法不需要构建二维面波速度分布，考虑了复杂介质中面波非大圆路径传播的情况，减少了路径假设误差，在国内很多地方得到了广泛应用（郭瑛霞等，2017;冯红武等，2019;李玲利等，2020）。另一些方法，全波形环境噪声成像（Gao，Shen，2014）、背景噪声伴随成像（Chen et al ，2014）则集中于提高成像分辨率，对计算机的性能有更高要求，不具有高效的计算速率，应用受到了限制。除此之外，随着研究深度的增加，面波成像分辨率会逐渐降低，无法有效反演深部结构细节。为降低面波反演的多解性，压制背景噪声成像方法中的反演假象，利用多种数据（郭震等，2015;吴萍萍等，2020）和方法联合反演获取地下介质速度结构已被应用到很多研究中。Fang等（2016）结合区域尺度体波双差成像和面波直接反演法提出了一种联合体波数据和面波数据反演的新方法（图2），并对南加州地区板块边界P波和S波速度结构进行了反演，这种方法结合了面波数据和体波数据的优势，能够更好地恢复P波和S波速度结构。Guo等（2015）利用由环境噪声提取的面波数据和远震P波接收函数联合反演方法获得了东北地区地壳和上地幔高分辨率S波速度模型，对兴安—蒙古造山带和长白山山脉中下地壳S波速度变化有了进一步认识。以往噪声面波成像的另一个不足在于大部分研究仅利用基阶频散进行结构反演，没有使用高阶频散信息。高阶振型面波虽然受制于能量、振幅变化不易提取（李雪燕等，2020），但在某些频率可能能量占优（罗银河等，2008），因此，很多研究尝试加入高阶面波频散，提高反演稳定性和精度（Yaoetal，2011;Gaurav et al ，2018）。Wang等（2019）提出了频率-贝塞尔变换法（F-J法）来提取噪声数据中的高阶频散信息。该方法是一种阵列分析方法，用F-J方法提取的频散曲线的质量取决于该阵列的范围和相邻台站之间的间隔，范围越大、结果分辨率越高，间隔越小、结果质量越好。李雪燕等（2020）采用此方法提取了上海市苏州河地区城市微动信号高阶频散信息，并对浅地表S波速度结构进行反演。将F-J方法与SPAC方法进行对比，可以发现F-J方法提取的高阶频散更清晰，高频部分（>13Hz）優势更明显。

2.2体波信号提取与体波成像

一般而言，噪声源分布与地表活动密切相关而缺乏一定深度，这使得互相关信号中面波占主导而体波振幅较弱，并且体波信号在台站间距较大时会更快衰减，因此互相关提取体波信号存在困难。尽管如此，体波信号包含很多关于地球深部核幔结构的信息，因此，不断有研究尝试利用噪声互相关方法提取体波信号。Roux等（2005）利用加州Parkfield 密集台阵记录的16个月的数据进行互相关提取，得到了P波信号。Zhan等（2010）在短周期频带（1～5s）成功恢复了莫霍面反射体波SmS和SmS2，证实了噪声源分布不均可以导致较弱的体波信号被其他波形掩盖。Poil等（2012）通过芬兰北部POLEnEt/LAPNEt地震阵列记录的数据进行互相关提取出了0.5～2Hz的高频体波，并与实际地震记录和合成地震图进行了对比，结果证明互相关可以提取出清晰的SmS和PmP波。Lin等（2013a）利用美国1181个台站垂直分量数据互相关成功提取了地核震相ScS和PKIKP 2，还识别出其它几个较弱的体波震相，包括P、S、PeP以及ScP/PeS。提取体波信号为噪声成像提供了一种新的可能，利用体波成像对地球内部结构进行研究发展起来。Nakata 等（2015）利用两种不同的滤波器从加州Long Beach 密集台阵10天记录数据计算的日互相关函数中提取体波，这两个滤波器能增强体波信号，同时压制非相干噪声，因此可以重建清晰的体波，在此基础上，又进行了三维体波成像，成像结果很好地显示了不同深度P波高低速异常。Wang等（2015）分析了1992—2012年57个全球宽带地震台阵测得的地震尾波的自相关信号，发现两种穿过地核的震相PKIKP 2、PKIIKP 2的走时残差在低纬度地区变化高达10s，认为内内核和外内核的各向异性分布差异导致了这种变化，这种不同取向和形式可能代表了地核演化过程中的一种转变。Feng等（2017）采用背景噪声互相关成功恢复了410km和660km不连续间断面反射信号P410P和P60P，揭示出华北克拉通两条剖面下方不连续间断面的分布变化，这种变化可能与热物质上涌和华北克拉通下方太平洋板块停滞的过渡带的影响有关。

2.3监测地下介质波速变化

通过监测地下介质波速变化推断地球内部成分、状态以及应力场的响应，进而预测火山爆发时间及其强度、探究断层深部应力变化、认识地震活动性，一直以来是地球物理学家们的研究重点。噪声相关方法不受地震活动性和震源参数变化影响，具有良好的时间连续性和成本优势，适宜于监测某时间段内地下介质速度变化。利用噪声相关研究波速变化，可以通过层析成像得到火山、地震爆发前后的速度结构图像，计算前后时间段的速度差值获得其相对扰动（严珊等，2016），还可以基于由噪声互相关或自相关提取的面波（Xu，Song，2009）或尾波EGFs信息，利用其计算介质受到扰动前后的走时偏移进而得到速度变化，常见方法有：尾波干涉法（CodaWave Interferometry， CWI）（Snieder， Science 2002;Gret et al ，2006）、拉伸法（Wegler，Sens-Schoenfelder， 2007： Hadziioannou et al， 2009）及移动窗口互谱法（Moving Window Cross -Spec-tral，MWCS）（Brenguier et al ，2008a，b;刘志坤，黄金莉，2010）。Gret等（2006）在实验室中利用CWI方法对不同岩石样品中单轴应力、水饱和度以及温度变化引起的波速变化进行了研究，证实了这种方法对研究板块边界及油气藏应力变化、监测地下水污染以及火山温度变化具有重要参考价值;Wegler和Sens-Schoenfelder（2007）通过Stretching 方法发现了日本新泻中部6.6级地震前后地壳内部速度的快速下降;Brenguier等（2008a，b）利用MWCS方法对富尔奈斯火山和圣安德烈斯断层地下介质波速变化进行了研究，发现火山喷发及大地震都会使面波速度明显降低;刘志坤和黄金莉（2010）与赵盼盼等（2012）分别采用互相关和自相关函数对汶川地震前后地壳速度变化进行了研究，与互相关函数不同的是自相關结果中显现了鲜水河断裂和龙门山断裂交汇区存在的同震速度增加区。

测量波速变化的过程中，不同的噪声相关方法、走时偏移测量的准确性等都会对结果产生影响。相较于互相关方法，自相关函数的尾波敏感区更靠近台站附近（赵盼盼等，2012），从而避免了仪器时移和噪声源分布对测量结果的影响，能够得到较浅层介质性质变化的响应。走时偏移的精准测量则与很多因素有关。首先，台站数量、间距会影响对地下介质的响应范围。台站数量多，反映的地下介质性质变化区域越大。台站间距越大，台站对之间EGFs的信噪比越低，测量结果越不稳定;其次，相关函数路径的数量直接影响测量结果的稳定性和精确性;最后，互相关函数叠加时长和滤波方法也对结果有影响，随着叠加时长增大，EGFs信噪比会逐渐提高，处理背景噪声时对NCF使用滤波方法，也可提高信噪比，如Hadziioannou（2011）、Stehly等（2015）以及温扬茂等（2019）分别利用S变换、曲波变换以及小波变换对NCF进行了处理，滤波的使用使得结果精度明显提高（图3）。

2.4地下衰减结构监测

地震波在传播的过程中受到地下岩石、矿物以及介质空隙、流体含量等的影响，能量会产生损耗而振幅发生衰减（周连庆，2016）。利用噪声相关可以恢复衰减结构，提升对地下介质的认识。Prieto 等（2009）通过测量作为距离函数的频率域相干性来恢复相速度和衰减系数a，由此计算了横穿南加州主要盆地的路径和其它路径的相速度与介质衰减的变化，与其它路径相比，穿过盆地的路径相速度更低、α更大、衰减更快。Cupillard和Capdeville（2010）采用数值实验研究了噪声源分布方式与one-bit和whitened 预处理对噪声相关振幅的影响，发现在噪声源均匀分布情况下，介质的几何扩展和固有衰减可以从相关信号中得到很好地恢复而不受预处理方式影响，但在非均匀分布情况下，几何扩展会受到影响。Weaver （2011）利用数值模拟的方法研究了one-bit归一化对EG-Fs振幅的影响，并提出了“temporal flattening ”的新方法，这种方法相比于one-bit能够更好地恢复衰减。周连庆（2016）在此基础上提出了一种“异步temporal flattening ”方法，该方法避免了对连续波形一致性和完整性的要求，且不会改变EG-Fs的相对振幅，因此更适宜于实际数据的处理。提取了相对振幅就可以用于研究地壳和上地幔衰减成像，这对了解地球内部温度和挥发性含量变化以及预测地震动幅度等现象至关重要。Lawrence 和Prieto （2011）采用扩展一维衰减方法对美国西部地壳和上地幔衰减的横向变化进行了成像。周连庆（2016）提出了背景噪声提取瑞利波振幅的流程，并对中国大陆进行了衰减结构二维成像，10s和20s的成像结果均表现出显著的横向不均匀性，历史强震活动性特征也与衰减结构有很好的对应性。

2.5地震定位

传统的地震定位方法大多依赖体波信息，受地下速度结构、台站分布、到时拾取精度等因素的影响，较少利用同样包含震源信息且受事件深度影响振幅的面波。噪声互相关方法包含面波反映的震源参数信息，能够对体波定位方法进行弥补，从而为地震预警、地质构造研究、地震层析成像、核试验检测等多方面应用提供精确的定位信息。Barmin等（2011）提出了一种基于密集台阵EGFs的震中定位方法，计算了布设在局部区域的临时台阵和远距离固定台阵一定周期的EGFs，并生成了区域格网，通过比较格林函数与固定台阵记录的该区域发生的地震的相关性进行定位。他们利用这种方法对发生在美国西部的一些地震和矿井坍塌事件进行了测试，结果显示定位误差平均不到1km，但可能会随着震源机制和深度而变化。Zhan等（2011）分别计算了震中附近的台站以及合成地震记录与远台的互相关函数，利用两者的时间差来校正2008年南加州Chino Hills 地震定位结果，这种方法虽然需要事先确定震源机制和深度，但得出的定位结果与P波定位具有很好的一致性。Zeng等（2014）利用噪声互相关和地震记录中Rayleigh 波群速度走时时差进行地震定位。他们利用此方法对西澳大利亚Kalannie地震进行测试，精度达到2km以内，并且实验发现群速度走时受震源机制和震源深度的影响很小，因此利用该方法时未知震源参数不会给地震定位结果带来实质性偏差。Xie等（2020）利用基于群速度走时的方法，对2017年格林兰岛Nuugaatsiag滑坡进行了重新定位，定位结果距离卫星图像给出的地点2.5km，优于传统定位方法，证明该算法可以用于滑坡精定位，为次生灾害预警提供帮助。

2.6噪声源分析

噪声互相关应用的前提是噪声源近似均匀分布，而实际噪声源分布是不均的，甚至存在一些局部信号使得格林函数的可靠性和频散曲线的提取质量明显降低，对实际应用产生干扰，因此需要对噪声源实际分布进行探究。一些学者研究了全球不同区域的噪声源分布，Stehly等（2006）采用归一化噪声背景能量流（Normalized Background Energy Flux ，NBEF）方法，研究了北美、西欧和坦桑尼亚的几个台站噪声能量的季节性平均方位分布，结果发现5～10s信号是由海洋涌浪与海岸的非线性相互作用产生的，10～20s信號的产生则与深海中的海浪活动有关;Yang和ritzwoller 2008）通过分析欧洲、南非和西藏以及北美洲地区互相关函数信噪比，对第一类地脉动（10～20s）和第二类地脉动（5～10s）的噪声源分布进行了探究，认为观察到的噪声方向性差异是传播和衰减的结果。另一些学者则针对某研究区范围的噪声源进行了探究，Gerstoft和Tanimoto（2007）采用聚束分析方法观测了南加州台站记录的噪声源空间分布季节性变化，结果显示冬季和夏季第一类地脉动噪声源分别来自于西北方向北太平洋和210°方位角处的南加州海岸;鲁来玉等（2009）采用同样方法对华北科学探测台阵地震噪声方位变化和季节性变化进行了研究（图4），聚束结果显示一般情况下夏季噪声源强度在每个周期大于冬季;王伟涛等（2011）采用NBEF方法对云南地区地震噪声的性质进行了分析，认为第二类地脉动噪声在夏季和冬季分别来源于印度洋和北太平洋;田原等（2020）通过计算台站记录的功率谱密度度量了四川盐源盆地内布设的短周期密集台阵的噪声水平，结果发现，人类活动和浅部松散沉积层对台站记录的高频噪声强弱有很大的影响。Wang等（2018）、王芳等（2020）利用基于NCF的慢度和位置聚束分析方法对ChinArray一、二期台阵接收到的P波类型噪声源进行了分析，发现这些信号由P、PP、PKPbe波干涉产生，主要来自于北大西洋、北太平洋和南大洋凯尔盖朗深海高原。不同区域噪声源的成分和时空分布是不同的，因而需要研究人员在实际应用中对本区域的噪声分布有充分的认识以保证数据处理的准确性。

此外，局部噪声源也引起了研究人员的兴趣。Shapiro 等（2006）对全球不同观测台站噪声互相关记录中观测到的非常强的、窄带26s信号进行了定位研究，在假设面波沿海洋大弧路径传播的条件下，这种微震信号起源于几内亚湾，并在斐济盆地存在镜像点。与他们不同的是，Zeng和Ni（2014）认为面波沿大陆小弧路径传播更有效，在此条件下，西太平洋中26s信号源不再为镜像点而为独立源，该源靠近Vanuatu 岛火山，可能由岩浆过程所激发。Zeng和Ni（2010）还对东亚地区台站对NCF中0.07～0.12Hz频段内存在的视速度高于Rayleigh 波的持续局部强信号进行了研究，发现其振幅表现出年际变化而非季节变化，定位结果表明该信号可能发源于日本九州阿苏山。Ma等（2013）同样对秘鲁南部互相关记录中快于面波的连续信号进行了定位研究，确定了位于火山区的5～10s强地震能量散射体。

2.7时移校正

准确的到时信息对于地震走时成像、波速变化等应用研究以及地震定位、地震预警等现实需求具有重要意义。但是，受制于GPS、场地和仪器软硬件状态等多种因素的影响，地震仪内部时间可能无法及时与外界进行同步，特别是海底地震仪（Ocean Bottom Seismometer ，OBS），受电磁波传播能力的影响，记录的数据往往存在线性和非线性的时间漂移。很多学者试图通过NCF校正时移，Stehly等（2007）提出了基于NCF的时间对称性分析（Time Symmetry Analysis ，TSA）方法，利用美国加利福尼亚州3个地震台数据计算的格林函数与参考格林函数之间的相位差来测定5～10s和10～20s周期的面波走时变化，以此为基础对仪器时移进行校正，这种方法随后得到了广泛应用，特别是对OBS记录包括水听器分量在内的不同分量线性漂移的校正（Hannemann et al ，2014;Liu et al ，2018;郑宏等，2020），提高了数据的准确性和可靠性。TSA方法要求NCF能够代表真实地下介质的格林函数，即需要各向同性的噪声源产生具有稳定峰值且波形相似的互相关函数正负因果分支。针对这种局限性，Gouedard等（2014）提出了Virtual Doublet Analysis 方法，通过寻找不同数据下计算的当前NCF和参考NCF之间的相对拉伸来监测介质中的时间速度变化，这种方法不需要NCF收敛到格林函数，只要其在研究期内是稳定的就是有效的，但是，当满足各向同性的条件下（即长时间尺度下），TSA方法仍具有更好的分辨率。洲际地震台间的短周期（<20s）背景噪声互相关结果是十分微弱的，对地震仪走时同步不太有效，因此，Xia等（2015）尝试利用洲际台站能够清晰记录的几内亚湾26s局部持续微震信号同步非洲、北美和欧洲台站的到时。通过与远震P波到时残差和TSA方法进行比较，证实了在满足一定信噪比和假设（26s信号空间位置固定）的情况下，该信号用于检测洲际和地区台站对的时钟漂移的可能性。时移校正对于地震学应用意义非凡，利用噪声互相关来校正时移为我们提供了一条崭新的路径，有待进一步深入研究。

3结论与讨论

背景噪声互相关方法被认为是21世纪地震学领域的重要发展之一，为研究人员提供了一种全新的视角来研究地球内部结构。随着固定台站以及流动台站的布设和加密，背景噪声互相关方法将推动人们对浅层地表和地壳的认识，为进一步了解地幔及更深的内部结构奠定基础。这种方法已经被运用于面波、体波速度结构成像，地下介质波速变化和衰减结构检测、地震定位、噪声源分析以及时移校正等多个方面。尽管国内外学者利用背景噪声进行了一系列有益的研究，但是仍存在一些问题需要进一步解决和完善：

（1）如何更好地提取高阶频散和体波信息，进行噪声成像。目前的研究工作中，提取高阶频散和体波信息受到台站间距、数据处理方法等多种因素的限制，应用较少，不易推广，有待更深入的研究。

（2）随着密集台阵的布设，庞大的数据量已经对数据处理方法提出了新的要求，如何在提高结果精度的同时自动化处理，并加快计算速度也是未来噪声应用必须要考虑的问题之一。

（3）联合反演能够有效地弥补面波成像方法的不足，但是目前面波成像局限于和地震学方法或数据联合反演，缺少同电磁、重力等多学科数据交叉反演，这是联合反演进一步发展的方向之一

（4）互相關函数振幅信息对于了解地下结构性质至关重要，但是现有的数据预处理方式，特别是涉及振幅信息的one-bit等归一化操作能多大程度地保留振幅信息、如何更高效保留振幅信息仍需要进一步研究。

（5）噪声互相关方法用于研究地下介质波速变化、衰减结构以及地震定位，能够为监测地下结构、流体、温度变化等提供约束，为灾害预警提供参考。但这需要实时监测和足够高的精度，目前仍无法达到实际应用程度，需要进一步努力。

（6）背景噪声的来源及分布是保证提取出的格林函数是否可靠的基础，对研究成果影响很大，目前对于背景噪声来源的研究还很粗浅，限制于局部，对于某个地方具体噪声信号来源的分析还未有充分的研究策略。

背景噪声互相关使得人们可以从原本杂乱无章的噪声中获取有关地下介质的重要信息，为研究地球内部结构和动力学演变提供约束。未来随着研究思路、方法的进步，在更密集的台站布设、更高的数据质量和更优的处理方法的支持下，背景噪声互相关将拥有更广阔的应用空间，在现有应用研究深化的同时，扩展出新的研究领域，为人们更加深刻地认识地球内部结构、构造演化提供更多值得关注的研究成果。

参考文献：

范兴利，陈棋福，郭震，2020。长白山火山区高精度Rayleigh 面波相速度结构与岩浆系统[J].岩石学报，36（7）：2081-2091。

房立华。2009。华北地区瑞利面波噪声层析成像研究[D].北京：中国地震局地球物理研究所

冯红武，颜文华，严珊，等。2019。背景噪声和地震面波联合反演渭河盆地及邻区壳幔S波速度结构[J].地震地质，41（5）：1185-1205.

付媛媛，高原。2016。东北地区背景噪声的Rayleigh 和Love波相速度层析成像[J].地球物理学报，59（2）：494-503。

顾勤平，康清清，张鹏，等。2020。郯庐断裂带中南段及邻区Rayleigh 波相速度与方位各向异性[J].地震地质，42（5）：1129-1152。

郭瑛霞，张元生，颜文华，等。2017。甘东南地区基于射线追踪面波频散三维成像[J].地震工程学报，39（2）：268-277。

郭震，陈永顺，殷伟伟。2015。背景噪声面波与布格重力异常联合反演：山西断陷带三维地壳结构[J].地球物理学报，58（3）：821-831.

孔祥艳。2019。天山地区地壳上地幔速度结构研究[D].北京：中国地震局地球物理研究所，

李玲利，黄显良，姚华建，等。2020。合肥市地壳浅部三维速度结构及城市沉积环境初探[J].地球物理学报，63（9）：3307-3323。

李雪燕，陈晓非，杨振涛，等。2020。城市微动高阶面波在浅层勘探中的应用：以苏州河地区为例[J].地球物理学报，63（1）：247-255.

刘志坤，黄金莉。2010。利用背景噪声互相关研究汶川地震震源区地震波速度变化[J].地球物理学报，53（4）：853-863。

鲁来玉，何正勤，丁志峰，等。2009。华北科学探测台阵背景噪声特征分析[J].地球物理学报，52（10）：2566-2572。

罗银河，夏江海，刘江平，等。2008。基阶与高阶瑞利波联合反演研究[J].地球物理学报，51（1）：242-249。

欧阳龙斌，李红谊，吕庆田，等。2015。长江中下游成矿带及邻区地壳剪切波速度结构和径向各向异性[J].地球物理学报，58（12）：4388-4402.

彭艳菊，黄忠贤，苏伟，等。2007。中国大陆及邻区海域地壳上地幔各向异性研究[J].地球物理学报，50（3）：752-759。

田原，瞿辰，王伟涛，等。2020。四川盐源盆地短周期密集台阵背景噪声分布特征分析[J]。地球物理学报，63（6）：2248-2261。

王芳，王伟涛，袁松湧。2020。大孔径地震台阵噪声互相关函数中体波信号的研究——以ChinArray二期数据为例[J].地球物理学报，63（9）：3370-3386。

王怀富，吴建平，周仕勇，等。2020。川西地区瑞雷波方位各向异性[J].地震学报，42（3）：293-305，377。

王仁涛，李志伟，包丰，等。2019。松辽盆地沉积层结构的短周期地震背景噪声成像研究[J].地球物理学报，62（9）：3385-3399。

王伟涛，倪四道，王宝善。2011。云南地区地脉动噪声特征分析研究[J].地震，31（4）：58-67。

温扬茂，高松，许才军，等。2019。利用双台站背景噪声分析2017年墨西哥Mw7.1地震震源区的地震波速变化[J].地球物理学报，62（8）：3024-3033.

吴萍萍，谭捍东，陶涛，等。2020。基于交叉梯度约束的电阻率法和背景噪声法三维联合反演研究[J].地球物理学报，63（10）：3912-3930.

徐义贤，罗银河。2015。噪声地震学方法及其应用[J].地球物理学报，58（8）：2618-2636。

嚴珊，杨立明，周龙泉，等。2016.2013年岷县漳县M6.6地震前后南北地震带北段面波速度变化特征[J].地震工程学报，38（3）：423-430.

杨润海，谭俊卿，向涯，等。2020。基于改进的线性调频Z变换的高精度地震波速干涉测量[J].地震研究，43（1）：1-9。

杨志高，陈运泰，张雪梅，等。2019。青藏高原东缘及东北缘S波速度

结构和径向各向异性[J].地球物理学报，62（12）：4554-4570曾求，储日升，盛敏汉，等。2020。基于地震背景噪声的四川威远地区浅层速度结构成像研究[J].地球物理学报，63（3）：944-955。

张忠杰，许忠淮。2013。地震学百科知识（五）——地震各向异性[J].国际地震动态，（6）：34-41。

赵盼盼，陈九辉，Campillo M ，等。2012。汶川地震区地壳速度相对变化的环境噪声自相关研究[J].地球物理学报，55（1）：137-145.

郑宏，范建柯，LeBM，等。2020。基于背景噪声互相关技术的雅浦俯冲带海底地震仪时间校正[J].地球物理学进展，35（2）：799-806.

周连庆。2016。地球介质衰减特性层析成像[D].北京：中国地震局地球物理研究所

Aki K.1957. Space and time spectra of stationary stochastie waves， with special reference to microtremors[ J]. Bulletin Earthq. res.inst.tokyo Uni，35：415-416.

Barmin M P， Levshin AL， Yang Y， et al.2011. Epicentral location based on Rayleigh wave Empirical Green's Funetions from ambient seismic noise[ J. Geophysical Journal International， 184（2）： 869-884

Bensen G D Ritzwoller MH.Barmin M p et a.2007.Processing seismic ambient noise data to obtain reliable broad-band surface wave dis-persion measurements [ J]. Geophysical Journal International， 169（3）：1239-1260.

Bensen G D， Ritzwoller M H， Shapiro N M.2008. Broadband ambient noise surface wave tomography across the United States[ J].Joumal of Geophysical Research， 113（ B5）： B05306

Brenguier F， Campillo M， Hadziioannou C， et al.2008a.Postseismic relax-ation along the sanandreas fault at parkfield from continuous seis-mologicalobservations[ J].Science， 321（5895）： 1478-1481

Brenguier F， Shapiro NM， Campillo M， el al.2008b.Towards forecasting voleanic eruptions using seismie noise J]. Nature Geoscience，I （2）：126-130.

Campillo M， Paul A.2003. Long-range correlations in the diffuse seismic coda[J].Science ，299（5606）：547-549.

Chen M， Huang H， Yao H， et al.2014. Low wave speed zones in the crust beneath SE Tibet revealed by ambient noise adjoint tomography[]Geophysical Research Letters， 41（2）： 334-340

Claerbout， Jon F.1968. Synthesis of a layered medium from its acoustic transmission responseLJ. Geophysics 33（2）： 264-269

Cupillard P， Capdeville Y.2010. On the amplitude of surface waves oltained by noise correlation and the capability to recover the attenua-tion： a numerieal approach [ J]. Geophysical Journal Intemational 181（3）：1687-1700.

Derode a.larose e tanter m. et a.2003. Recovering the Green's fune-tion from field-field correlations in an open scattering medium（L）[ J]. Journal of the AcoustiealSoeiety of America， 113（6）： 2973-2976.

DuvallT， Jefferies SM， Harvey Jw， et al.1993. Time-distance helio-seismology [ J]. Nature， 362（6419）： 430-432.

Fang h. Yao h. Zhang h. et a 2015. Direct inversion of surface wave dispersion for three - dimensional shallow crustal structure based on ray tracing： methodology and applieation[]. Geophysical Journal International ，201（3）：1251-1263.

Fang h. Zhang h Yao h. et a 2016. a new algorithm for three - dimen-sional joint inversion of hody wave and surface wave data and its ap-plication to the Southern California plate boundary region [ J]. Jour-nal of Geophysical Research： Solid Earth， 121（5）： 3557-3569

Feng J， Yao H， Poli P， et al.2017. Depth variations of 410 km and 660 km discontinuities in eastern North China Craton revealed by ambi-ent noise interferometry [ J]. Geophysical Research Letters， 44（16）8328-8335.

Gao H， Shen Y.2014. Upper mantle structure of the Cascades from full wave ambient noise tomography： Evidence for 3D mantle upwelling in the back-are [ J]. Earth and Planetary Science Letters， 390：222-233.

Gaurav T Eleonore S， Aurelien M， et al.2018. Joint inversion of first over-tone and fundamental mode for deep imaging at the Valhall oil field using ambient noise[ J]. Geophysical Journal International：214（1）：122-132.

Gerstoft p Tanimoto t.2007. a vear of microseisms in southern Califomia [ J]. Geophysical Research Letters， 34（20）： 120304

Gouedard P， Seher T.Meguire JJ. et al.2014. Correction of ocean-bot-tom seismometer instrumental clock errors using ambient seismic noise[ J]. Bulletin of the Seismological Society of America， 104（3）1276-1288.

GretA， Snieder R， Seales J.2006. Time-lapse monitoring of rock proper-ties with coda wave interferometry []Joumal of Geophysical Re-search-Solid Earth， 111（ B3）： B03305.

Guo Z， Chen YJ， Ning J， et al.2015. High resolution 3-D crustal struc-ture beneath NE China from joint inversion of ambient noise and re-ceiver functions using NECESS Array data[ J. Earth and Planetary Science letters 416.-1

Hadziioannou C L. Larose e， Coutant o. et al.2009. Stability of monito-ring weak changes in multiply scattering media with ambient noise correlation： laboratory experiments [ J]. Journal Of The Acoustical Society Of Ameriea.125（6）： 3688-3695.

Hadziioannou C， Larose E， Baig A， et al.2011. Improving temporal resolu-tion in ambient noise monitoring of seismie wave speed[ J].Joumal of Geophysical Research， 116（ B7）： B07304.

Hannemann K， Kriiger F， Dahm T.2014. Measuring of clock drift rates and statie time offsets of ocean bottom stations by means of ambient noise[ J]. Geophysical Joumnal International， 196（2）： 1034-1042

Lawrence JF， Prieto G A.2011. Attenuation tomography of the westem U-nited States from ambient seismienoise[ J].Joumal of Geophysical Research， 116（ B6）： B06302.

Li D， Wang G， Lin R， et al.2019. The shallow erustal s-velocity strue-ture of the Longmenshan fault zone using ambient noise tomography of a seismic dense array []Earthquake Science， 32（5-6）： 197-206.

Li H， Su w. Wang C-Y， et al.2009.Ambient noise Rayleigh wave tomo-graphy in western Sichuan and eastern Tibet[ J. Earth and Planetary Science Letters， 282（1-4）： 201-211

Li Z， Ni S， Zhang B， et al.2016. Shallow magma chamber under the Wu-dalianchi Volcanic Field unveiled by seismic imaging with dense array[ J]. Geophysical Research Letters， 43（10）： 4954-4961

Lin FC， Li D， Clayton R W， et al.2013b.High-resolution 3D shallow crustal structure in Long Beach， Califomnia： Application of ambient noise tomography on a dense seismiearray[ J]. Geophysics， 78（4）：Q45-Q56.

Lin F-C， Moschetti MP， Ritzwoller M H.2008. Surface wave tomography of the western United States from ambient seismic noise： Rayleigh and Love wave phase velocity maps[ J]. Geophysical Journal Intema-tional ，173（1）：281-298.

Lin F-c. Tsai v c.schmandt b. et al 2013a. Extracting seismic core phases with array interferometry [ J]. Geophysical Research Letters 40（6）：1049-1053.

Liu Y， Liu C， Tao C， et al.2018. Time correction of the ocean bottom seis-mometers deployed at the southwest Indian ridge using ambient noise cross-correlation J]. Acta OceanologieaSinica， 37（5）： 39-46.

May.ClavtonRw. Tsai vc. et al.2013. Locating a seatterer in the ac-tivevolcanie area of southern peru from ambient noise cross-corre-lation[ J]. Geophysical Journal International， 192（3）： 1332-1341

NakataN， Chang JP， Lawrence JF， et al.2015. Body wave extraction and tomography at long Beach， Califomnia， with ambient-noise interferom-etry [ J].Joumal of Geophysical Research： Solid Earth， 120（2）：1159-1173.

Nishida K， Montagner JP， Kawakatsu H.2009. Global surface wave tomo-graphy using seismic hum[ J]. Science， 326（5949）： 112-112

Poli P， Pedersen HA， Campillo M.2012. Emergence of body waves from cross-correlation of short period seismic noise [ J]. Geophysical lournal Intemational.188（2）： 549-558.

PrietoGA. Lawrence JF， Beroza G C.2009. Anelastic Earth structure from the coherency of the ambient seismiefield[ J]. Journal of Geo-physical Research， 114（ B7）： B07303

Roux P， Sabra KG， Gerstoft P， et al.2005. P-waves from cross-correla-tion of seismie noise [ J].Geophysieal Research Letters， 32（19）：L19303.

Sabra KG， Gerstoft P， Roux P， et al.2005. Surface wave tomography from microseisms in Southern Califomia[ J]. Geophysical Research Let-ters，32（14）：L14311.

Shapiro NM， Campillo M， Stehly L， et al.2005. High-resolution surface - wave tomography from ambient seismic noise [ J]. Science， 307（5715）：1615-1618.

ShapiroN M， Campillo M.2004. Emergence of broadband Rayleigh waves from correlations of the ambient seismienoise[ J]. Geophysical Re-search Letlers， 31（7）： L07614.

ShapiroNM， Ritzwoller MH， Bensen G D.2006. Source location of the 26 see microseism from cross - correlations of ambient seismie noise [ J.Geophysieal Research Letters， 33（18）： L18310

Shen w.Ritzwoller MH. d Kang. et a.2016. a seismic reference model for the crust and uppermost mantle beneath China from surface wave dispersion [ J]. Geophysical Journal International， 206（2）： 954-979.

Snieder R 2004. Extracting the Green's function from the correlation of co-da waves： A derivation based on stationary phase [ J]. Physical Re-viewF ，69（4）：046610.

Snieder， Seience R J 2002. Coda wave interferometry for estimating non-linear behavior in seismic velocity[ J]. Science， 295（5563）： 22532255.

Stehly L， Campillo M， Shapiro N M.2006. A study of the seismie noise from its long - range correlation properties []. Journal of Geophysical Research， 111（ B10）： B10306.

Stehly L， Campillo M， Shapiro N M.2007.Traveltime measurements from noise correlation stability and detection of instrumental time - shifts ]. Geophysical Journal International， 171（1）： 223-230

Stehly， Laurent， Froment， et al.2015. Monitoring seismic wave velocity changes associated with the Mw7.9 Wenchuan earthquake： increas-ing the temporal resolution using curvelet filters [ J]. Geophysical Journal Intemational， 201（3）： 1939-1949

Tang Y， Chen Y J， Zhou S， et al.2013. Lithosphere structure and thick-ness beneath the North China Craton from joint inversion of ambient noise and surface wave tomography [ J]. Journal of Geophysical Re-search： Solid Earth， 118（5）： 2333-2346

Wang J， Wu G， Chen X.2019. frequency-bessel transform method for ef-fective imaging of higher-mode rayleigh dispersion curves from ambi-entseismie noise data[ J]. Journal of Geophysical Research：Solid Earth ，124（4）：3708-3723.

Wang T， Song x， Xia HH .2015. Equatorial anisotropy in the inner part of Earth's inner core from autocorrelation of earthquake coda[J].Na-tureGeoseience， 8（3）： 224-227.

Wang w， Gerstoft P， Wang B.2018. Seasonality of p wave microseisms from NCF-based beamforming using ChinArray [ J]. Geophysical Journal International， 213（3）： 1832-1848.

Weaver r L Lobkis o 2003.Flastie wave thermal nuctuations ultra-sonie waveforms by correlation of thermal phonons J.Joumnal of the Acoustical Society of America， 113（5）： 2611-2621.

Weaver r l lobkis o 2001. ultrasonics without a source： thermal luc-tuation correlations at MHz frequencies [ J]. Physical Review Letters 87（13）：134301.

Weaver r l lobkis o L 2002. On the emergence of the green's function in the correlations of a diffuse field： pulse - echo using thermal pho-nons[ J]. Ultrasonics， 40（1-8）： 435-439

Weaver R L.2011. On the amplitudes of correlations and the inference of attenuations， specific intensities and site factors from ambient noise [ J].ComptesRendus Geoscience， 343（8-9）： 615-622.

Wegler U， Sens - Schoenfelder C.2007. Fault zone monitoring with pas-sive image interferometry[ J]. Geophysical Journal International， 168（3）：1029-1033.

Xia Y， Ni S， Zeng x， et al.2015. Synchronizing intercontinental seismic networks using the 26 s persistent localized microseismiesource[]Bulletin of the Seismological Society of America， 105（4）： 2101-2108.

Xie J. Chu R. Ni S.2020.Relocation of the 17 June 2017 Nuugaatsiag Greenland） landslide based on Green s functions from ambient seis-mic noises [. Journal of Geophysical Research： Solid Earth， 125（5）：e2019JB018947.

Xie J， Ritzwoller MH. Shen W， et al.2013. Crustal radial anisotropy across Eastem Tibet and the Westem Yangtze Craton[ J]. Journal of Geophysical Research： Solid Earth， 118（8）： 4226-4252

Xu ZJ， Song X.2009. Temporal changes of surface wave velocity associat-ed with major Sumatra earthquakes from ambient noise correlation [ J. Proceedings of the National Aeademy of Seiences of the United States of America， 106（34）： 14207-14212

Yang Y.Ritzwoller M h.2008. Characteristics of ambient seismic noise as a source for surface wave tomography [ J]. Geochemistry Geophysics Geosystems， 9： Q02008

YaoH.Gouedard P. Collins JA. et al.2011. Structure of young East Pa-cific Rise lithosphere from ambient noise correlation analysis of fun-damental-and higher - mode Scholte - Rayleigh waves [ ComptesRendus Geoscience， 343（8-9）： 571-583

Yao H， Van Der Hilst RD， De Hoop MV.2006. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-I. Phase velocity maps [ J]. Geophysical JoumnalIntema-tional ，166（2）：732-744.

Yao H， Van Der Hilst RD， Montagner J-P.2010. Heterogeneity and ani-sotropy of the lithosphere of sE Tibet from surface wave array tomo-graphy[ J. Journal of Geophysical Research， 115（ B12）： B12307

Zeng x， Ni S.2010. A persistent localized microseismie source near the Kyushu Island， Japan [ J ] Geophysical Research Letters 37：L24307.

Zeng X， Ni S.2014. Evidence for an independent 26-s microseismic source near the Vanuatu IslandsJ]. Pure and Applied Geophysics 171（9）：2155-2163

Zeng x， Xie J， Ni S.2014. Ground truth location of earthquakes by use of ambient seismic noise from a sparse seismic network： A case study in Western Australia[ J]. Pure and Applied Geophysics， 172（6）： 1397-1407.

Zhan Z， Ni S， Helmberger D V， et al.2010. Retrieval of Moho-reflected shear wave arrivals from ambient seismie noise [ J].Geophysica Journal International， 182（1）： 408-420

Zhan Z， Wei S， Ni S， et al.2011. Earthquake centroid locations using cali-bration from ambient seismienoise[]. Bulletin of the Seismologica Society of Ameriea， 101（3）： 1438-1445

Zhang Y， Li H， Huang Y， et al.2020. Shallow structure of the Longmen Shan Fault Zone from a high - density， short-period seismie array [ J. Bulletin of the Seismological Society of America， 110（1）38-48.

Review of the Application of Seismic Ambient NoiseCross-correlation Method

CHEN Yuxin'，TANG Mingshuai3

（1.Lanzhou Institute of Seismology，China Earthquake Administration，Lanzhou 730000，Gansu，China）

（2.Xinjiang Pamir Intracontinental Subduction National Field Observation and Research Station，Urumqi 830011，Xinjiang，China）

（3.Urumqi Institute of Central Asia Earthquake，China Earthquake Administration，Urumqi 830011，Xinjiang，China）

Abstract

Obtaining Green's function，which can reflect the information about the underground media by cross correla- ting the ambient noise data recorded by seismic stations， is the frontier and hot spot of theoretical research and practical application in the field of seismology. This paper firstly introduces the evolution of the noise cross-correla- tion application， and then expounds its application to surface-wave tomography，body-wave signal extraction and body-wave imaging，monitoring of variation and attenuation of the wave-velocity of underground media， earthquake location，analysis of noise source， and time-shifting correction. Finally，theprospeet of the applica- tion of ambient noise cross-correlation is discussed.

Keywords： ambient noise;Green's function; dispersion curves; tomography imaging;medium change