多种面波层析成像方法及其在青藏高原的应用与对比  

Multiple surface wave tomography methods and their applications to the Tibetan Plateau

作　　者：李伦 蔡晨 付媛媛[4] 方洪健 Li Lun;Cai Chen;Fu Yuanyuan;Fang Hongjian(Department of Earth Sciences and Engineering,Sun Yat-Sen University,Guangzhou 510275,China;Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai),Zhuhai 519000,China;Guangdong Provincial Key Lab of Geodynamics and Geohazards,Guangzhou 510275,China;Key Laboratory of Earthquake Prediction,Institute of Earthquake Forecasting,China Earthquake Administration,Beijing 100036,China)

机构地区：[1]中山大学地球科学与工程学院,广州510275 [2]南方海洋科学与工程广东省实验室(珠海),珠海519000 [3]广东省地球动力作用与地质灾害重点实验室,广州510275 [4]中国地震局地震预测研究所,北京100036

出　　处：《地球与行星物理论评（中英文）》2023年第2期174-196,共23页Reviews of Geophysics and Planetary Physics

基　　金：第二次青藏高原综合科学考察研究资助项目(2019QZKK0701);国家自然科学基金青年科学基金资助项目(41804043,42106067);广东省自然科学基金面上资助项目(2019A1515011729);广州市科技计划项目(202102020456);广东省引进人才创新创业团队资助项目(2017ZT07Z066)。

摘　　要：面波层析成像是一种广泛应用的获取地壳与上地幔地震波速度与各向异性结构的地球物理方法.本论文综述了面波层析成像方法的简要历史,阐述了多种常用的面波层析成像方法(双台法、双平面波法、程函方程法、Helmholtz方程法、背景噪声成像法与直接面波层析成像法)的基本原理及其优缺点.双台法的理论与实际使用简单,但该方法要求震源与两个台站需近似位于同一大圆弧路径,在台站较少且记录时间较短的区域,获取的相速度水平分辨率偏低.双平面波法能一定程度克服地震波多路径传播与散射对频散的影响,但其对面波波形数据要求较为严格,且通常适用于区域地震台网.程函方程法和Helmholtz方程法可直接从地震记录同时获取相速度与方位各向异性,计算方便快速,无需经过正演与反演过程,但这两种方法要求台站分布密度要高,不适用于台站间距大且分布稀疏的区域.与程函方程法相比,Helmholtz方程法不仅考虑了波形的相位,还利用了其振幅,能获取更准确的相速度与方位各向异性信息.背景噪声成像法的优点是无需震源就可获取高分辨率地壳尺度的成像结果,但通常缺乏长周期面波的信息,难以约束岩石圈深部与软流圈结构.直接面波层析成像法能直接从台站间的面波频散数据通过反演获取三维剪切波(S波)速度结构与方位各向异性信息,省去了反演相速度图的中间步骤.此外,我们对比了多种方法在青藏高原获取的相速度结果.结果表明,多种面波层析成像方法获取的同一周期相速度结果高度吻合,主要特征表现在:在中长周期的相速度图中,青藏高原内部主要以低速为主,而周缘区域(如,柴达木盆地、四川盆地等)以相对高速为主,这表明青藏高原中下地壳与上地幔的流变强度均比其周缘区域要弱,青藏高原的岩石圈变形受控于周缘块体的阻挡.在青藏高原东�Surface wave tomography is a widely used geophysical method to measuring seismic velocity and anisotropy in the crust and upper mantle.This paper briefly reviewed the history of surface wave tomography,and summarized the principles and advantages of multiple surface wave tomography methods(i.e.,two-station method,two-plane-wave method,Eikonal and Helmholtz surface wave tomography,ambient noise tomography,and direct surface wave tomography).The theory and application of the two-station method are straightforward,but this method requires the earthquake and seismic stations to be located in the same great circle.This restriction results in a low lateral resolution phase velocity map for regions where stations are sparse and deployment times short.The twoplane-wave method can overcome the effect of scattering and multipathing of seismic wave propagation on phase velocity dispersions,but this method requires high-quality surface wave waveforms and is usually suitable for regional seismic networks.The Eikonal and Helmholtz surface wave tomography can directly compute the phase velocities and azimuthal anisotropy straightforwardly without any processes of forward modeling and inversion,but this method is limitedly suitable for high-density and orderly seismic arrays.In comparison to the Eikonal surface wave tomography method,the Helmholtz surface wave tomography approach analyzes both the phase and amplitude of the waveform and can produce better results.Ambient noise tomography can obtain high-resolution crustal structure without seismic events.However,this method is difficult to obtain long-period phase velocities,resulting in a lack of constraints on the mantle and lithosphere structure.Direct surface wave tomography can obtain the shear-wave velocity and azimuthal anisotropy from dispersive curves without the inversion process of phase velocity maps.We conducted detailed comparisons of the phase velocity maps at the short-intermediate periods(20~40 s)previously obtained from these surface wave tomography methods in

关 键 词：面波层析成像 瑞利波 勒夫波 双台法 双平面法 程函方程法 背景噪声 青藏高原 

分 类 号：P315[天文地球—地震学]