Love wave tomography in Sichuan-Yunnan area from ambient noise
-
摘要: 本文选取了川滇地区98个固定台站记录到的三分量地震背景噪声数据，采用互相关方法提取了勒夫波互相关函数，并利用自适应时频分析方法获取了勒夫波群速度频散曲线，经反演得到周期为8—30 s的勒夫波群速度分布图像.层析成像结果显示：短周期的勒夫波群速度分布图像呈明显的横向不均匀性，且与地表地质和构造特征基本一致，其中四川盆地呈不均匀的低速异常，盆地内成都平原地区的群速度相对低于盆地中部的丘陵地区，速度分界线为遂宁与峨眉山之间的连线，四川盆地内的群速度变化反映出沉积层厚度的变化情况；攀枝花地区呈高速异常，可能与古地幔活动有关，幔源物质以侵入岩和底侵岩浆的形式停留在地壳的不同深度，从而形成高速异常的特征.本文结果为了解川滇地区的构造运动提供了地震学线索，并为下一步研究地壳径向各向异性奠定了基础.Abstract: We collected seismic ambient noise data recorded at 98 permanent seismic stations in Sichuan-Yunnan area to image crustal structure. The Love wave cross-correlation functions are obtained by cross-correlation of ambient noise. The group velocity dispersion curves for Love wave are measured by automated time-frequency analysis. Love wave group velocity is mapped in the periods between 8 s and 30 s. The tomographies show seismic images at short periods exhibit apparent horizontal heterogeneities, which is largely consistent with geological features and geophysical studies. The Sichuan basin appears as low group velocity anomaly, and the Chengdu plain shows lower group velocity relative to the hilly area in Sichuan basin, with the line connecting Suining with Emeishan being the boundary between higher velocity and lower one. The variation of group velocity anomalies in Sichuan basin reflects variable sedimentary thicknesses. The Panzhihua area appears as high velocity anomaly, which is associated with an ancient mantle plume. The materials from mantle is remained at the different depth in crust in the form of intrusive rocks and magmatic underplating, resulting in the characteristics of high velocity anomaly. These seismic images provide necessary information for the research on crustal radial anisotropy and understanding of tectonic process.
-
Keywords:
- ambient noise /
- tomography /
- Love wave /
- crustal structure /
- Sichuan-Yunnan area
-
引言
被誉为“世界屋脊”的青藏高原是地球科学界最为关注的活动碰撞造山带，其形成与印度板块和欧亚板块的碰撞有关.在过去几十年里，为了解释青藏高原的隆升与演化机制，不同研究人员根据所获得的观测数据提出了多个板块碰撞模式, 其中下地壳流模型(Royden et al, 1997, 2008; Clark, Royden, 2000)和刚性块体挤出模型(Tapponnier et al, 2001)最具代表性, 这两个模型均认为青藏高原东南缘(川滇地区)是青藏高原地壳物质运移的重要场所和通道.青藏高原的隆升效应使川滇地区形成极其复杂且特殊的地质构造，而四川盆地作为扬子地台的一部分，是发育在青藏高原东缘的大型构造沉积盆地，对青藏高原东缘的构造模式和高原内部物质的运移起着重要的作用.四川盆地西南侧为峨眉山大火成岩省，其形成与古地幔柱有关.大火成岩省记录了地球某一特定时期巨量物质和能量由地球内部向外迁移的过程，其形成与地球表面的形态变化及矿产资源的形成有密切关系(徐义刚等, 2013).峨眉山大火成岩省所涉及的成矿作用和古地幔柱作用遗迹等的研究(如位置规模、组分来源等)仍需地球物理探测结果的进一步约束.此外，川滇地区内有大量的活动构造带(如龙门山断裂带、红河断裂带和小江断裂带等)活动强烈、地震频发，同样需要来自地球物理方面的约束和研究，而高分辨率地壳层析成像研究恰巧能够为探讨川滇地区的构造演化及动力学研究提供重要信息.
背景噪声层析成像已成为目前探测地球深部结构的有效手段之一(Shapiro et al, 2005; Yao et al, 2006; Yang et al, 2007; Fang et al, 2010; 唐有彩等, 2011; Saygin，Kennett, 2012; 房立华等, 2013; Tang et al, 2013; 付媛媛，高原, 2016).利用互相关方法可以从地震背景噪声数据中提取出瑞雷波和勒夫波互相关函数(Lin et al, 2008; Moschetti et al, 2010).一般认为, 由于噪声源形成机制的复杂性和水平向地震记录干扰较大，从噪声中提取高信噪比勒夫波较瑞雷波提取更为困难.因此，在背景噪声层析成像中，利用垂直向背景噪声记录提取瑞雷波的研究相对较多，而利用水平向背景噪声记录提取勒夫波的研究相对较少.相较于瑞雷波，同周期勒夫波的敏感深度更浅，有利于改善浅部地层的分辨率(付媛媛，高原, 2016).此外，利用勒夫波可以反演出SH波速度结构，结合由瑞雷波反演得到的SV波速度结构，可对壳幔介质的径向各向异性(Moschetti et al, 2010)进行深入研究.
Yao等(2006, 2008)首次将背景噪声层析成像方法引入到川滇地区的地壳结构和动力学模式研究中，利用垂直向的噪声记录提取了瑞雷波互相关函数，并结合天然地震资料构建了周期为10—150 s的相速度分布图，进而反演得到了SV波速度结构，其结果显示川滇地区低速层主要位于中地壳或下地壳，而且低速层的分布形态十分复杂. Yang等(2010, 2012)利用垂直向背景噪声数据构建了青藏高原及其周边地区的瑞雷波相速度分布图和SV波速度结构，分析了青藏地区低速层的分布范围及可能成因. Xie等(2013)利用位于青藏高原东侧326个台站的背景记录资料，构建了三维SH波和SV波速度结构. Zheng等(2015)收集了青藏高原东南缘170个台站的背景噪声数据，提取了瑞雷波互相关函数，最终构建了三维SV波速度结构模型，其结果表明，青藏高原东南缘地区的低速层埋深在20—35 km中下地壳内，厚度为5—20 km, 速度降幅为5.4%—10.8%.范莉苹等(2015)利用青藏高原东南缘地区557个台站的背景噪声资料，反演获得了瑞雷波群速度图像，并进一步分析探讨其分布特征和构造意义. Zheng等(2017)利用川滇地区121个台站的背景噪声数据, 采用贝叶斯方法反演获取了川滇地区8—40 s周期的相速度分布图，并对该地区的地壳流空间分布和动力学过程进行了深入分析.这些成果为川滇地区地震活动背景研究奠定了基础，促进和加深了对川滇地区壳幔结构及其动力学的认识.
鉴于目前利用背景噪声层析成像构建瑞雷波和SV波速度结构对川滇地区地壳结构和动力学模型开展的研究较多，而利用勒夫波相速度和群速度分布图进行详细分析解释的研究鲜有报道.鉴于此，本文拟基于川滇地区地震台站记录的三分量连续记录数据，构建勒夫波群速度分布图，并结合已有的地质与地球物理研究成果，对其速度分布及其可能存在的地球动力学及地质含义进行探讨，以期对川滇地区速度结构研究提供新的约束，更好地获取该地区地壳内部的深浅构造形态.
1. 数据和方法
本文使用2008年1月—2009年12月云南省和四川省区域地震台网的98个台站所记录的背景噪声连续波形，台站分布如图 1所示.由于固定台站存在时钟误差、极性相反及定位误差等问题，本文参考前人研究成果，剔除其中有问题台站的数据(Niu，Li, 2011; Zhou et al, 2012; Xie et al, 2013)，然后参照Bensen等(2007)和Lin等(2008)的数据处理流程，对连续波形三分量数据进行处理，从而提取瑞雷波和勒夫波互相关函数.数据处理流程主要分为5步：单台数据预处理、互相关函数计算和叠加、频散曲线测量、质量控制和误差分析、面波层析成像.在单台预处理阶段，首先将数据去除直流分量及线性趋势，重采样到10 Hz，截取成长度为6小时的数据单元，然后进行1—50 s带通滤波、时域归一化和谱白化等.
![]()
图 1 研究区区域地质构造简图(修改自Wang et al，2003)及台站分布Figure 1. Location of the seismic stations (denoted by triangles) and major faults (modified from Wang et al, 2003) in the studied areaPurple lines are major faults. F1: Longmenshan fault; F2: Xianshuihe fault; F3: Anninghe-Zemuhe fault; F4: Xiaojiang fault; F5: Honghe fault; F6: Langcangjiang fault; F7: Nujiang fault; F8: Lijiang-Xiaojinhe fault; F9: Mile-Shizong fault. The main geological units are: SG: Songpan-Garze fold system; YZ: Yangtze platform; SJ: Sangjiang fold system; SB: Sichuan basin; SC: South China fold system; ZG: Zuogong-Gengma fold system; BT: Bomi-Tengchong fold system为了提高计算效率，对于每一对台站对的三分量记录，我们首先计算出多个方向的互相关函数(Z-Z，E-E，E-N，N-N，N-E)，再根据台站间的方位角和反方位角计算出径向(R-R)和切向分量(T-T)(Lin et al, 2008)，其中垂向分量(Z-Z)和径向分量(R-R)代表瑞雷波互相关函数，切向分量(T-T)代表勒夫波互相关函数.图 2为台站JIH与其余各台站路径上的垂向分量和切向分量的互相关函数.
![]()
由于噪声源分布不均匀，互相关函数正负时间经常出现不对称.为了提高信噪比，研究中将负信号反序后叠加于正信号.在数据质量控制方面，我们只保留信噪比高于20且台间距大于3个波长的互相关函数(Bensen et al, 2007; 郑定昌等, 2014)，以得到可靠的频散曲线.本文采用自适应时频分析技术计算勒夫波群速度频散曲线(Levshin, Ritzwoller, 2001; Bensen et al, 2007; Lin et al, 2008)，图 3给出了基于台站对JIH-PWU背景噪声数据得到的勒夫波互相关函数和由自适应时频分析技术求取的勒夫波群速度频散曲线示意图.
![]()
本文将研究区划分为0.5°×0.5°的网格，采用面波层析成像反演方法(Ditmar，Yanovskaya, 1987; Yanovskaya，Ditmar, 1990)分别对不同周期的频散曲线进行反演，得到相应周期的群速度分布图像.该面波层析成像反演方法是传统Backus-Gilbert一维方法在面波二维反演情况下的推广.该面波层析成像程序在反演群速度的同时，还根据射线密度和方向分布来计算每个格点的空间平均分辨率核函数，以确定层析成像的空间分辨率.目前，该方法已成为面波层析成像研究中广泛应用的方法之一(Guo et al, 2009, 2013; 房立华等, 2009, 2013; Fang et al, 2010; 潘佳铁等, 2011, 2015; 郑定昌等, 2014).反演过程中，正则化参数选取为0.2，该值能够确保数据误差较小且结果比较光滑.每次迭代完成，将走时残差大于3倍均方根的频散曲线剔除后，再进行下一次迭代反演.
图 4和图 5分别给出了川滇地区不同周期的射线分布和射线路径数量统计.由图 5可以看出，在10—22 s的范围内，各周期的射线均超过了3500条；其余周期的射线数量也达到2000条以上.密集均匀的射线分布(图 4)为层析成像的可靠性提供了保证.面波层析成像的分辨率主要受到射线数量和射线方位覆盖的影响，地震射线的数量越多、交叉越好，分辨率则越高.图 6给出了不同周期勒夫波层析成像结果的横向分辨率，可以看出，中心区域的横向分辨率可达到40—50 km，总体上可以保证川滇地区地壳结构较高的横向分辨率.自20 s周期起，随着周期的增大，符合条件的射线数量逐渐减少，研究区边缘的分辨率逐渐变低；在高于30 s的频段内，勒夫波互相关函数的信噪比显著降低，符合质量控制条件的射线数量较小，故未对其进行后续研究.
![]()
图 4 不同周期(T=8, 12, 16, 20, 25, 30 s)的射线分布Figure 4. Ray path coverage at different periods (T=8, 12, 16, 20, 25, 30 s)![]()
图 5 不同周期的射线路径数量统计Figure 5. Statistics on the number of the ray paths at different periods![]()
图 6 不同周期(T=8, 12, 16, 20, 25, 30 s)勒夫波层析成像结果的横向分辨率Figure 6. Estimated resolution of Love wave tomography results at different periods (T=8, 12, 16, 20, 25, 30 s)2. 勒夫波群速度结构
勒夫波具有频散特性，不同周期勒夫波具有不同的穿透深度.一般情况下，群速度敏感核函数的峰值为相速度的两倍.因此，在探索介质的各向异性方面，面波群速度比相速度更有效(Suresh et al, 2015).为了更好地研究川滇地区地壳介质各向异性，便于定性地解释勒夫波群速度分布图，本文计算不同周期基阶勒夫波的敏感核函数(Rodi et al, 1975).计算敏感核函数时, 参考CRUST1.0模型(Laske et al, 2013)和Yang等(2012)的研究结果，上地幔模型则采用AK135模型(Kennett et al, 1995).图 7给出了不同周期勒夫波和瑞雷波群速度频散所对应的剪切波深度响应范围.可以看到，在相同周期上，勒夫波采样于更浅的深度.总体来看，勒夫波的敏感核主要受浅部地层速度的影响，有利于更准确地揭示浅部地层构造信息.
![]()
图 7 不同周期T勒夫波(a)和瑞雷波(b)的群速度敏感核函数Figure 7. Sensitivity kernels of group velocity for Love wave (a) and Rayleigh wave (b) at different periods T图 8a和图 8b分别为周期8 s和12 s的勒夫波群速度分布图，主要反映了地壳浅部速度结构特征，与地表区域地质构造(如沉积层厚度、基底埋深等)有较大相关性.由图 8a可以看出：最为显著的低速异常位于四川盆地西部，且四川盆地整体显示为低速异常；松潘—甘孜造山带为高速异常, 与大地电磁测深结果的“浅部高阻”特征对应，可能与该地区三叠纪复理石杂岩沉积有关(张乐天等, 2012)；龙门山断裂带为明显的高低速分界线，攀枝花附近区域和峨眉山—昭通一带，显示为高速异常，其成因推测与古地幔柱活动有关；红河断裂和小江断裂带显示为低速异常；弥勒—师宗断裂南侧的滇东南区域是稳定的抬升区，显示为高速异常.此外，低速异常区分布在楚雄盆地、兰坪—思茅弧后盆地、西昌盆地和滇西北裂陷区(大致位置对应于图 8a中盐源、剑川与大理围成的区域)，保山以南的地块为低速异常.弥勒—师宗断裂显示为高低速分界线.
![]()
图 8 不同周期的勒夫波群速度分布图(黑线为四川盆地的轮廓线)Figure 8. The Love wave group velocity maps at different periods where the black line indicates the boundary of Sichuan basin(a) T=8 s; (b) T=12 s; (c) T=16 s; (d) T=20 s (e) T=25 s; (f) T=30 s正如图 8b所示，随着周期的逐渐增大，低速区分布面积逐步缩小.群速度分布图中显示出两个最明显的极低值低速异常区：一个为四川盆地内的低速异常区，其面积相对于图 8a中的结果缩小，群速度值有所增大，速度最低区仍然位于四川盆地的西北侧(龙门山断裂带附近)；另一个为兰坪—思茅盆地内的低速异常区，群速度最低值的速度异常区位于思茅与元江之间.根据遂宁—阿坝和孟连—思茅—元江—马龙这两条人工测深剖面观测，这两个极低值速度异常区均分布有极厚的沉积层(张中杰等, 2005；嘉世旭等, 2014). Burchfiel等(2008)的研究表明，四川盆地西部(龙门山断裂附近)的沉积层厚度可达10 km以上)，思茅至元江西侧的沉积厚度达到12 km (张中杰等, 2005).从松潘—甘孜褶皱系经龙门山断裂带至四川盆地的群速度特征与遂宁—阿坝人工测线(嘉世旭等, 2014)以及大地电磁测深方法(张乐天等, 2012；Zhao et al, 2012)所得的地球物理特征相吻合.此外，小江断裂带西侧存在南北向条带状低速异常体，丽江—小金河断裂呈低速异常.
由勒夫波群速度敏感核函数(图 7a)可知，虽然相对长周期群速度图中不能完全剥离浅层物质对其的影响，但是随着周期的增大，群速度分布图可以反映出S波速度在深度上的渐变过程.图 8c和图 8d分别为周期16 s和20 s的群速度分布图，二者特征分布较为相似.四川盆地仍然为低速异常区域，但范围相对缩小；低速异常还分布在思茅—元江、剑川—盐源和保山地区；楚雄—攀枝花—康定北一线呈现不规则的南北向条带状高速异常分布，昭通地区的高速异常逐渐减弱；保山南部总体呈现出低速异常. 20 s周期的群速度分布(图 8d)虽然与Li等(2010)关于四川地区同周期勒夫波群速度结果所显示的形态相似，但仍存在一定差异，特别是在遂宁—重庆一带和Li等(2010)研究区域的西南边缘地区，这可能是由于不同分辨率所致. Li等(2010)的研究区域在本研究中处于相对中心的位置，射线密度更高，具有更高的分辨率.
图 8e和图 8f分别为周期25 s和30 s的群速度分布图，可见四川盆地的东南部地区逐渐呈现高速异常特征.简单来说，四川盆地内可以以遂宁—峨眉山的连线为速度差异的分界线，且群速度分布形态差异也可以用大地电磁测深的高阻异常体分布来解释(张乐天等, 2012).丽江—小金河断裂以北的川滇菱形地块北部逐渐变化为低速异常，攀枝花附近的区域仍然为高速异常.小江断裂—昭通—西昌一带呈低速异常，盐源—剑川一带和小江断裂附近持续为低速异常区(周期8—30 s).这两个区域的大地热流平均值均大于85 mW/m2(汪缉安等, 1990; 汪集旸，黄少鹏, 1990)，因此推测低速异常可能与深部热作用有关.虽然腾冲地区也为大地热流高值区，但由于其已处于本文研究区的边缘地区，射线密度较小，分辨率较差，故未作深入分析.
3. 讨论与结论
3.1 攀枝花地区
峨眉山大火成岩省分布于扬子克拉通西缘，主要由亚碱性及偏碱性的基性火山熔岩和火山碎屑组成，普遍认为其形成与古地幔柱活动有关(徐义刚等, 2013; 徐涛等, 2015).在距今256—259万年的晚古生代, 云南、四川和贵州等3省境内曾发生过大规模的玄武岩喷发，攀枝花地区位于峨眉山玄武岩喷发前由地幔柱活动引起的地壳穹隆构造的核心部位(He et al, 2003).古地幔柱的火山作用所引起的地壳热异常及岩浆入侵和喷发等，会造成地壳内物质成分的改变.这种古地幔柱活动的遗迹在地球物理性质上主要表现为地震波速度、vP/vS波速比、密度、各向异性和磁性等与周围地区的差异.本文中8—30 s的群速度分布图显示，在攀枝花附近，一直存在高速异常，且在不同周期上的形态各异：在8—20 s范围内，攀枝花地区的群速度主要在3.3—3.5 km/s之间变化；在25 s和30 s群速度分布图上，群速度在3.5—3.7 km/s之间变化，攀枝花地区下方连续的高速特征可能暗示了其深部构造背景.丽江—清镇宽角地震资料研究结果显示：攀枝花地区上地壳P波速度高至6 km/s以上，下地壳为6.9—7.2 km/s，在整个地壳内，P波速度明显高于剖面的其它地区；莫霍面局部上隆(徐涛等，2015).程远志等(2015)基于大地电磁测深剖面数据的研究显示了攀枝花地区从地表至30 km深度存在高阻体，相对于周围地区的高速、高阻可能反映了该地区的密度异常；Zhang等(2013)的研究认为其中的玄武岩自下而上由低钛玄武岩变为高钛玄武岩；熊绍柏等(1993)根据人工地震探测剖面的结果，认为攀枝花地区上地壳高速与成矿盐体密切相关，而且可能有深部构造背景；吴建平等(2013)关于P波速度层析成像的研究结果显示，攀枝花附近的高速异常体可能与地幔柱活动所导致的基性和超基性幔源物质侵入地壳有关；徐涛等(2015)推测高速异常可能是二叠纪地幔柱活动引起的底侵作用和岩浆上侵的结果.结合多种观测结果及已有认识，本文推测：由地幔柱在地幔深度下熔融产生的高温岩浆除一部分喷出地表外，其余部分以侵入岩和底侵岩浆的形式残留在地壳的不同深度；幔源物质的侵入形成了密度和力学强度大于周边的情形，而且地壳内的不均匀分布导致了不同周期上高速异常的形态差别.
3.2 四川盆地
四川盆地属扬子地块，其演化始于三叠纪晚期.该盆地的下陷以及沉积物的充填是扬子西缘与秦岭和松潘—甘孜褶皱带发生陆内汇聚作用的表现形式(王二七，尹纪云, 2009).四川盆地内部地形平缓，广布的沉积岩层近水平分布.以龙泉山断裂为分界线，可将四川盆地分为西部的成都平原和盆中丘陵地区(张乐天等, 2012)，其中成都平原地区为白垩系、第三系地层分布和第四系大片覆盖区，也是中新生代主要坳陷区，而盆中丘陵地区，由于其稳定的刚性基底，盖层褶皱受到限制，是盆地内褶皱最弱的地区，且构造极其平缓，均属低平构造.四川盆地由于褶皱弱、断层少、构造定向性差、多旋卷构造体系(江为为等, 2001)，其内部群速度值在群速度分布图上(8—20 s)呈现两个速度异常体，相对表现为“成都平原低，盆中丘陵地区高”，这种特性即反映了四川盆地内的沉积层变化特点，其原因可能与青藏高原的逆冲推覆作用有关.
3.3 结论
本文收集了川滇地区98个固定台站2008—2009年的三分量连续背景噪声数据，使用互相关技术，构建了川滇地区不同周期(8—30 s)勒夫波群速度的分布图像以揭示川滇地区上中地壳的结构特征.在研究区内的大部分地区，横向分辨率约为40—80 km.研究结果显示：川滇地区存在强烈的横向不均匀性，高低速异常交错分布；沉积层厚度变化是短周期群速度分布图中速度分布形态变化的主要原因；在浅层，四川盆地内部整体表现为“成都平原低，盆中丘陵地区高”的特征；随着周期的增大，四川盆地逐渐转变为高速的克拉通特征；攀枝花地区在不同周期均为高速异常，该异常可能为以侵入岩和底侵岩浆的形式停留在地壳不同深度的幔源物质所致.
对比已有的基于背景噪声的瑞雷波研究(Li et al, 2009; Yang et al, 2010, 2012; 范莉苹等, 2015；Zheng et al, 2015)，在相同短周期(如8 s，10 s)情况下，勒夫波群速度分布图揭示出了更多浅层地壳构造特征，如小尺度的西昌盆地、四川盆地的内部结构形态、思茅—元江低速区等；本研究中的结果与其它地球物理观测结果具有很好的一致性，在一定程度上提高了研究区浅层结构的分辨率.川滇地区不同周期群速度分布图所呈现的横向不均匀性和随周期变化的纵向特征对认识地块结构构造具有重要的参考价值.
当周期大于30 s时，勒夫波互相关函数的信噪比迅速降低.本文现阶段得到的勒夫波互相关函数仅适于对中上地壳的结构调查.今后的工作中将注重数据方法的改进，以获取高信噪比的长周期勒夫波格林函数.下一阶段将结合瑞雷波相速度和群速度分布图，进一步反演出SH波和SV波速度结构及方位各向异性，由此判明该地区深浅构造的地壳形变特征，了解该地区不同时期的地质作用，为该地区的地震孕震模式提供基础数据.
审稿专家为本文提供了诸多宝贵的修改意见，闵照旭、杨润海和闻学泽为本研究收集并提供了云南省和四川省固定台网的连续波形资料，研究中使用了SEIZMO程序包、杨英杰博士的部分程序和Yanovskaya教授的面波层析成像计算程序，作者在此一并表示感谢. -
![]()
图 1 研究区区域地质构造简图(修改自Wang et al，2003)及台站分布
Figure 1. Location of the seismic stations (denoted by triangles) and major faults (modified from Wang et al, 2003) in the studied area
Purple lines are major faults. F1: Longmenshan fault; F2: Xianshuihe fault; F3: Anninghe-Zemuhe fault; F4: Xiaojiang fault; F5: Honghe fault; F6: Langcangjiang fault; F7: Nujiang fault; F8: Lijiang-Xiaojinhe fault; F9: Mile-Shizong fault. The main geological units are: SG: Songpan-Garze fold system; YZ: Yangtze platform; SJ: Sangjiang fold system; SB: Sichuan basin; SC: South China fold system; ZG: Zuogong-Gengma fold system; BT: Bomi-Tengchong fold system
![]()
![]()
![]()
图 4 不同周期(T=8, 12, 16, 20, 25, 30 s)的射线分布
Figure 4. Ray path coverage at different periods (T=8, 12, 16, 20, 25, 30 s)
![]()
图 5 不同周期的射线路径数量统计
Figure 5. Statistics on the number of the ray paths at different periods
![]()
图 6 不同周期(T=8, 12, 16, 20, 25, 30 s)勒夫波层析成像结果的横向分辨率
Figure 6. Estimated resolution of Love wave tomography results at different periods (T=8, 12, 16, 20, 25, 30 s)
![]()
图 7 不同周期T勒夫波(a)和瑞雷波(b)的群速度敏感核函数
Figure 7. Sensitivity kernels of group velocity for Love wave (a) and Rayleigh wave (b) at different periods T
-
程远志, 汤吉, 陈小斌, 董泽义, 肖骑彬, 汪利波. 2015.南北地震带南段川滇黔接壤区电性结构特征和孕震环境[J].地球物理学报, 58(11): 3965-3981. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201511008&dbname=CJFD&dbcode=CJFQ Cheng Y Z, Tang J, Chen X B, Dong Z Y, Xiao Q B, Wang L B. 2015. Electrical structure and seismogenic environment along the border region of Yunnan, Sichuan and Guizhou in the south of the North-South Seismic Belt[J]. Chinese Journal of Geophysics, 58(11): 3965-3981 (in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201511008&dbname=CJFD&dbcode=CJFQ
范莉苹, 吴建平, 房立华, 王未来. 2015.青藏高原东南缘瑞利波群速度分布特征及其构造意义探讨[J].地球物理学报, 58(5): 1555-1567. doi: 10.6038/cjg20150509 Fan L P, Wu J P, Fang L H, Wang W L. 2015. The characteristic of Rayleigh wave group velocities in the southeastern margin of the Tibetan Plateau and its tectonic implications[J]. Chinese Journal of Geophysics, 58(5): 1555-1567 (in Chinese). doi: 10.6038/cjg20150509
房立华, 吴建平, 吕作勇. 2009.华北地区基于噪声的瑞利面波群速度层析成像[J].地球物理学报, 52(3): 663-671. doi: 10.1002/cjg2.v52.3 Fang L H, Wu J P, Lü Z Y. 2009. Rayleigh wave group velocity tomography from ambient seismic noise in North China[J]. Chinese Journal of Geophysics, 52(3): 663-671 (in Chinese). doi: 10.1002/cjg2.v52.3
房立华, 吴建平, 王未来, 王长在, 杨婷. 2013.华北地区勒夫波噪声层析成像研究[J].地球物理学报, 56(7): 2268-2279. doi: 10.6038/cjg20130714 Fang L H, Wu J P, Wang W L, Wang C Z, Yang T. 2013. Love wave tomography from ambient seismic noise in North-China[J]. Chinese Journal of Geophysics, 56(7): 2268-2279 (in Chinese). doi: 10.6038/cjg20130714
付媛媛, 高原. 2016.东北地区背景噪声的Rayleigh和Love波相速度层析成像[J].地球物理学报, 59(2): 494-503. doi: 10.6038/cjg20160209 Fu Y Y, Gao Y. 2016. Phase velocity tomography of Rayleigh and Love waves using ambient noise in Northeast China[J]. Chinese Journal of Geophysics, 59(2): 494-503 (in Chinese). doi: 10.6038/cjg20160209
嘉世旭, 刘保金, 徐朝繁, 刘志, 酆少英, 张建狮, 林吉焱, 田晓峰, 刘巧霞, 郭文斌. 2014.龙门山中段及两侧地壳结构与汶川地震构造[J].中国科学:地球科学, 44(3): 497-509. doi: 10.1007/s11430-013-4744-9 Jia S X, Liu B J, Xu Z F, Liu Z, Feng S Y, Zhang J S, Lin J Y, Tian X F, Liu Q X, Guo W B. 2014. The crustal structures of the central Longmenshan along and its margins as related to the seismotectonics of the 2008 Wenchuan earthquake[J]. Science China Earth Sciences, 57(4): 777-790. doi: 10.1007/s11430-013-4744-9
江为为, 刘伊克, 郝天珧, 宋海斌. 2001.四川盆地综合地质、地球物理研究[J].地球物理学进展, 16(1): 11-23. http://d.wanfangdata.com.cn/Periodical/dqwlxjz200101002 Jiang W W, Liu Y K, Hao T Y, Song H B. 2001. Comprehensive study of geology and geophysics of Sichuan basin[J]. Progress in Geophysics, 16(1): 11-23 (in Chinese). http://d.wanfangdata.com.cn/Periodical/dqwlxjz200101002
潘佳铁, 吴庆举, 李永华, 张风雪, 张广成. 2011.华北地区瑞雷面波相速度层析成像[J].地球物理学报, 54(1): 67-76. http://d.wanfangdata.com.cn/Periodical/dqwlxb201101008 Pan J T, Wu Q J, Li Y H, Zhang F X, Zhang G C. 2011. Rayleigh wave tomography of the phase velocity in North China[J]. Chinese Journal of Geophysics, 54(1): 67-76 (in Chinese). http://d.wanfangdata.com.cn/Periodical/dqwlxb201101008
潘佳铁, 李永华, 吴庆举, 丁志峰. 2015.青藏高原东南部地区瑞雷波相速度层析成像[J].地球物理学报, 58(11): 3993-4006. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201511010&dbname=CJFD&dbcode=CJFQ Pan J T, Li Y H, Wu Q J, Ding Z F. 2015. Phase velocity maps of Rayleigh waves in the southeast Tibetan Plateau[J]. Chinese Journal of Geophysics, 58(11): 3993-4006 (in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201511010&dbname=CJFD&dbcode=CJFQ
唐有彩, 陈永顺, 杨英杰, 丁志峰, 刘瑞丰, 冯永革, 李鹏, 俞春泉, 魏松峤, 范文渊, 王海洋, 周仕勇, 宁杰远. 2011.华北克拉通中部地区背景噪声成像[J].地球物理学报, 54(8): 2011-2022. http://d.wanfangdata.com.cn/Periodical/ctta201401007 Tang Y C, Chen Y S, Yang Y J, Ding Z F, Liu R F, Feng Y G, Li P, YU C Q, Wei S Q, Fan W Y, Wang H Y, Zhou S Y, Ning J Y. 2011. Ambient noise tomography in North China craton[J]. Chinese Journal of Geophysics, 54(8): 2011-2022 (in Chinese). http://d.wanfangdata.com.cn/Periodical/ctta201401007
汪缉安, 徐青, 张文仁. 1990.云南大地热流及地热地质问题[J].地震地质, 12(4): 367-377. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzdz199004012&dbname=CJFD&dbcode=CJFQ Wang J A, Xu Q, Zhang W R. 1990. Heat flow data and some geologic-geothermal problems in Yunnan Province[J]. Seismology and Geology, 12(4): 367-377 (in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzdz199004012&dbname=CJFD&dbcode=CJFQ
汪集旸, 黄少鹏. 1990.中国大陆地区大地热流数据汇编(第二版)[J].地震地质, 12(4): 351-363, 366. http://d.wanfangdata.com.cn/Periodical/dqwlxb200105005 Wang J Y, Huang S P. 1990. Compilation of heat flow data in the China continental area (2nd edition)[J]. Seismology and Geology, 12(4): 351-363, 366 (in Chinese).
王二七, 尹纪云. 2009.川西南新生代构造作用以及四川原型盆地的破坏[J].西北大学学报:自然科学版, 39(3): 359-367. http://d.wanfangdata.com.cn/Periodical/xbdxxb200903003 Wang E Q, Yin J Y. 2009. Cenozoic multi-stage deformation occurred in southwest Sichuan: Cause for the dismemberment of the proto-Sichuan basin[J]. Journal of Northwest University: Natural Science Edition, 39(3): 359-367 (in Chinese). http://d.wanfangdata.com.cn/Periodical/xbdxxb200903003
吴建平, 杨婷, 王未来, 明跃红, 张天中. 2013.小江断裂带周边地区三维P波速度结构及其构造意义[J].地球物理学报, 56(7): 2257-2267. doi: 10.6038/cjg20130713 Wu J P, Yang T, Wang W L, Ming Y H, Zhang T Z. 2013. Three dimensional P-wave velocity structure around Xiaojiang fault system and its tectonic implications[J]. Chinese Journal of Geophysics, 56(7): 2257-2267 (in Chinese). doi: 10.6038/cjg20130713
熊绍柏, 郑晔, 尹周勋, 曾晓献, 全幼黎, 孙克忠. 1993.丽江—攀枝花—者海地带二维地壳结构及其构造意义[J].地球物理学报, 36(4): 434-444. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx199304003&dbname=CJFD&dbcode=CJFQ Xiong S B, Zheng Y, Yin Z X, Zeng X X, Quan Y L, Sun K Z. 1993. The 2-D structure and it's tectonic implications of the crust in the Lijiang-Panzhihua-Zhehai region[J]. Chinese Journal of Geophysics, 36(4): 434-444 (in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx199304003&dbname=CJFD&dbcode=CJFQ
徐涛, 张忠杰, 刘宝峰, 陈赟, 张明辉, 田小波, 徐义刚, 滕吉文. 2015.峨眉山大火成岩省地壳速度结构与古地幔柱活动遗迹:来自丽江—清镇宽角地震资料的约束[J].中国科学:地球科学, 45(5): 561-576. doi: 10.1007/s11430-015-5094-6 Xu T, Zhang Z J, Liu B F, Chen Y, Zhang M H, Tian X B, Xu Y G, Teng J W. 2015. Crustal velocity structure in the Emeishan large igneous province and evidence of the Permian mantle plume activity[J]. Science China Earth Sciences, 58(7): 1133-1147. doi: 10.1007/s11430-015-5094-6
徐义刚, 何斌, 罗震宇, 刘海泉. 2013.我国大火成岩省和地幔柱研究进展与展望[J].矿物岩石地球化学通报, 32(1): 25-39. http://d.wanfangdata.com.cn/Periodical/kwysdqhxtb201301002 Xu Y G, He B, Luo Z Y, Liu H Q. 2013. Study on mantle plume and large igneous provinces in China: An overview and perspectives[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 32(1): 25-39 (in Chinese). http://d.wanfangdata.com.cn/Periodical/kwysdqhxtb201301002
张乐天, 金胜, 魏文博, 叶高峰, 段书新, 董浩, 张帆, 谢成良. 2012.青藏高原东缘及四川盆地的壳幔导电性结构研究[J].地球物理学报, 55(12): 4126-4137. doi: 10.6038/j.issn.0001-5733.2012.12.025 Zhang L T, Jin S, Wei W B, Ye G F, Duan S X, Dong H, Zhang F, Xie C L. 2012. Electrical structure of crust and upper mantle beneath the eastern margin of the Tibetan Plateau and the Sichuan basin[J]. Chinese Journal of Geophysics, 55(12): 4126-4137 (in Chinese). doi: 10.6038/j.issn.0001-5733.2012.12.025
张中杰, 白志明, 王椿镛, 吕庆田, 滕吉文, 李继亮, 孙善学, 王新征. 2005.冈瓦纳型和扬子型地块地壳结构:以滇西孟连—马龙宽角反射剖面为例[J].中国科学: D辑, 35(5): 387-392. doi: 10.1360/03yd0547 Zhang Z J, Bai Z M, Wang C Y, Lü Q T, Teng J W, Li J L, Sun S X, Wang X Z. 2005. Crustal structure of Gondwana-and Yangtze-typed blocks: An example by wide-angle seismic profile from Menglian to Malong in western Yunnan[J]. Science in China: Series D, 48(11): 1828-1836. doi: 10.1360/03yd0547
郑定昌, 盖增喜, 杨润海, 闵照旭, 唐有彩, 姜明明, 庞卫东. 2014.云南地区背景噪声层析成像[J].地震学报, 36(4): 602-614. http://www.dzxb.org/Magazine/Show?id=28988 Zheng D C, Ge Z X, Yang R H, Min Z X, Tang Y C, Jiang M M, Pang W D. 2014. Broadband ambient noise tomography in Yunnan Province[J]. Acta Seismologica Sinica, 36(4): 602-614 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28988
Bensen G D, Ritzwoller M H, Barmin M P, Levshin A L, Lin F, Moschetti M P, Shapiro N M, Yang Y. 2007. Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements[J]. Geophys J Int, 169(3): 1239-1260. doi: 10.1111/gji.2007.169.issue-3
Burchfiel B C, Royden L H, van der Hilst R D, Hager B H, Chen Z, King R W, Li C, Lü J, Yao H J, Kirby E. 2008. A geological and geophysical context for the Wenchuan earthquake of 12 May 2008, Sichuan, People's Republic of China[J]. GSA Today, 18(7): 4-11. doi: 10.1130/GSATG18A.1
Clark M K, Royden L H. 2000. Topographic ooze: Building the eastern margin of Tibet by lower crustal flow[J]. Geology, 28(8): 703-706. doi: 10.1130/0091-7613(2000)28<703:TOBTEM>2.0.CO;2
Ditmar P G, Yanovskaya T B. 1987. A generalization of Backus-Gilbert method for estimation of lateral variations of surface wave velocity[J]. Phys Solid Earth, 23(6): 470-477.
Fang L H, Wu J P, Ding Z F, Panza G F. 2010. High resolution Rayleigh wave group velocity tomography in North China from ambient seismic noise[J]. Geophys J Int, 181(2): 1171-1182.
Guo Z, Gao X, Yao H J, Li J, Wang W M. 2009. Midcrustal low-velocity layer beneath the central Himalaya and southern Tibet revealed by ambient noise array tomography[J]. Geochem Geophys Geosys, 10(5): Q05007. doi: 10.1029/2009GC002458.
Guo Z, Gao X, Shi H, Wang W M. 2013. Crustal and uppermost mantle S-wave velocity structure beneath the Japanese islands from seismic ambient noise tomography[J]. Geophys J Int, 193(1): 394-406. doi: 10.1093/gji/ggs121
He B, Xu Y G, Chung S L, Xiao L, Wang Y M. 2003. Sedimentary evidence for a rapid, kilometer-scale crustal doming prior to the eruption of the Emeishan flood basalts[J]. Earth Planet Sci Lett, 213(3/4): 389-405.
Kennett B L N, Engdahl E R, Buland R. 1995. Constraints on seismic velocities in the Earth from traveltimes[J]. Geophys J Int, 122(1): 108-124. doi: 10.1111/gji.1995.122.issue-1
Laske G, Masters G, Ma Z T, Pasyanos M. 2013. Update on CRUST1.0: A 1-degree global model of Earth's crust[J]. Geophys Res Abstr, 15: EGU2013-2658.
Levshin A L, Ritzwoller M H. 2001. Automated detection, extraction, and measurement of regional surface waves[J]. Pure Appl Geophys, 158(8): 1531-1545. doi: 10.1007/PL00001233
Li H Y, Su W, Wang C Y, Huang Z X. 2009. Ambient noise Rayleigh wave tomography in western Sichuan and eastern Tibet[J]. Earth Planet Sci Lett, 282(1/4): 201-211.
Li H Y, Su W, Wang C Y, Huang Z X, Lü Z Y. 2010. Ambient noise Love wave tomography in the eastern margin of the Tibetan Plateau[J]. Tectonophysics, 491(1/4): 194-204.
Lin F C, Moschetti M P, Ritzwoller M H. 2008. Surface wave tomography of the western United States from ambient seismic noise: Rayleigh and Love wave phase velocity maps[J]. Geophys J Int, 173(1): 281-298. doi: 10.1111/gji.2008.173.issue-1
Moschetti M P, Ritzwoller M H, Lin F, Yang Y. 2010. Seismic evidence for widespread western-US deep-crustal deformation caused by extension[J]. Nature, 464(7290): 885-889. doi: 10.1038/nature08951
Niu F L, Li J. 2011. Component azimuths of the CEArray stations estimated from P-wave particle motion[J]. Earthquake Science, 24(1): 3-13. doi: 10.1007/s11589-011-0764-8
Rodi W L, Glover P, Li T M C, Alexander S S. 1975. A fast, accurate method for computing group-velocity partial derivatives for Rayleigh and Love modes[J]. Bull Seismol Soc Am, 65(5): 1105-1114.
Royden L H, Burchfiel B C, King R W, Wang E C, Chen Z L, Shen F, Liu Y P. 1997. Surface deformation and lower crustal flow in eastern Tibet[J]. Science, 276(5313): 788-790. doi: 10.1126/science.276.5313.788
Royden L H, Burchfiel B C, van der Hilst R D. 2008. The geological evolution of the Tibetan Plateau[J]. Science, 321(5892): 1054-1058. doi: 10.1126/science.1155371
Saygin E, Kennett B L N. 2012. Crustal structure of Australia from ambient seismic noise tomography[J]. J Geophys Res, 117(B1): B01304.
Shapiro N M, Campillo M, Stehly L, Ritzwoller M H. 2005. High-resolution surface-wave tomography from ambient seismic noise[J]. Science, 307(5715): 1615-1618. doi: 10.1126/science.1108339
Suresh G, Bhattacharya S N, Teotia S S. 2015. Crust and upper mantle velocity structure of the northwestern Indian Peninsular shield from inter-station phase velocities of Rayleigh and Love waves[J]. Ann Geophys, 58(2): S0215.
Tang Y C, Chen Y J, Zhou S Y, Ning J Y, Ding Z F. 2013. Lithosphere structure and thickness beneath the North China Craton from joint inversion of ambient noise and surface wave tomography[J]. J Geophys Res, 118(5): 2333-2346.
Tapponnier P, Xu Z Q, Roger F, Meyer B, Arnaud N, Wittlinger G, Yang J S. 2001. Oblique stepwise rise and growth of the Tibet Plateau[J]. Science, 294(5547): 1671-1677. doi: 10.1126/science.105978
Wang C Y, Chan W W, Mooney W D. 2003. Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications[J]. J Geophys Res, 108(B9): ESE13-11.
Xie J Y, Ritzwoller M H, Shen W S, Yang Y J, Zheng Y, Zhou L Q. 2013. Crustal radial anisotropy across eastern Tibet and the western Yangtze craton[J]. J Geophys Res, 118(8): 4226-4252. doi: 10.1002/jgrb.50296
Yang Y J, Ritzwoller M H, Levshin A L, Shapiro N M. 2007. Ambient noise Rayleigh wave tomography across Europe[J]. Geophys J Int, 168(1): 259-274. doi: 10.1111/gji.2007.168.issue-1
Yang Y J, Zheng Y, Chen J, Zhou S Y, Celyan S, Sandvol E, Tilmann F, Priestley K, Hearn T M, Ni J F, Brown L D, Ritzwoller M H. 2010. Rayleigh wave phase velocity maps of Tibet and the surrounding regions from ambient seismic noise tomography[J]. Geochem Geophys Geosys, 11(8): Q08010. doi: 10.1029-2010GC003119/
Yang Y J, Ritzwoller M H, Zheng Y, Shen W S, Levshin A L, Xie Z J. 2012. A synoptic view of the distribution and connectivity of the mid-crustal low velocity zone beneath Tibet[J]. J Geophys Res, 117(B4): B04303.
Yanovskaya T B, Ditmar P G. 1990. Smoothness criteria in surface wave tomography[J]. Geophys J Int, 102(1): 63-72. doi: 10.1111/gji.1990.102.issue-1
Yao H J, van der Hilst R D, De Hoop M V. 2006. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis: Ⅰ. Phase velocity maps[J]. Geophys J Int, 166(2): 732-744. doi: 10.1111/gji.2006.166.issue-2
Yao H J, Beghein C, van der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis: Ⅱ. Crustal and upper-mantle structure[J]. Geophys J Int, 173(1): 205-219. doi: 10.1111/gji.2008.173.issue-1
Zhang Y, Ren Z Y, Xu Y G. 2013. Sulfur in olivine-hosted melt inclusions from the Emeishan picrites: Implications for S degassing and its impact on environment[J]. J Geophys Res, 118(8): 4063-4070. doi: 10.1002/jgrb.50324
Zhao G, Unsworth M J, Zhan Y, Wang L F, Chen X B, Jones A G, Tang J, Xiao Q B, Wang J J, Cai J T, Li T, Wang Y Z, Zhang J H. 2012. Crustal structure and rheology of the Longmenshan and Wenchuan MW7.9 earthquake epicentral area from magnetotelluric data[J]. Geology, 40(12): 1139-1142. doi: 10.1130/G33703.1
Zheng D C, Saygin E, Cummins P, Ge Z X, Min Z X, Cipta A, Yang R H. 2017. Transdimensional Bayesian seismic ambient noise tomography across SE Tibet[J]. J Asian Earth Sci, 134: 86-93. doi: 10.1016/j.jseaes.2016.11.011
Zheng X, Zhao C P, Zhou L Q, Zheng S H. 2015. 3D shear-wave velocity structure beneath the southeastern Tibetan Plateau from ambient noise[J]. Bull Seismol Soc Am, 105(3): 1371-1382. doi: 10.1785/0120140211
Zhou L Q, Xie J Y, Shen W S, Zheng Y, Yang Y J, Shi H X, Ritzwoller M H. 2012. The structure of the crust and uppermost mantle beneath South China from ambient noise and earthquake tomography[J]. Geophys J Int, 189(3): 1565-1583. doi: 10.1111/gji.2012.189.issue-3
-
期刊类型引用(10)
1. 马永，张海江，高磊，宋程，毕金孟，高也. 滇西地区地壳三维精细结构成像与构造特征研究. 地球物理学报. 2023(09): 3674-3691 . 
百度学术
2. 常英娜，梁春涛，曹飞煌，周鲁，廖江涛，陆威帆，王朝亮. 基于波场梯度法研究安宁河—则木河断裂带速度结构. 地球物理学报. 2022(08): 2886-2903 . 百度学术
3. 张丽晓，闫俊岗，张双凤. 基于背景噪声的晋冀鲁豫交界地区瑞利波群速度层析成像. 地震. 2022(01): 155-168 . 百度学术
4. 杨建文，叶泵，高琼，张华英，陈佳. 2021年云南漾濞M_S6.4地震前后主动源观测走时数据的变化. 地震工程学报. 2021(04): 767-776 . 百度学术
5. 马永，张海江，高磊，毕金孟. 漾濞Ms6.4地震三维精细速度结构与地震活动研究（英文）. Applied Geophysics. 2021(04): 579-591+595 . 百度学术
6. 张演，郑定昌，林国良，吕帅. 基于三分量数据的云南地区背景噪声源分析. 地震研究. 2020(04): 689-700+768 . 百度学术
7. ZhiGao Yang，XiaoDong Song. Ambient noise Love wave tomography of China. Earth and Planetary Physics. 2019(03): 218-231 . 必应学术
8. 杨志高，陈运泰，张雪梅，宋晓东. 青藏高原东缘及东北缘S波速度结构和径向各向异性. 地球物理学报. 2019(12): 4554-4570 . 百度学术
9. 周青云，陈俊磊. 气枪震源激发条件对走时变化观测结果的影响. 地震研究. 2018(02): 264-272+344 . 百度学术
10. 谭夏露，房立华，王未来，吴建平. 安宁河-则木河断裂带及周边地区Rayleigh波群速度背景噪声成像研究. 中国地震. 2018(03): 400-413 . 百度学术
其他类型引用(10)
下载:
计量
- 文章访问数: 1000
- HTML全文浏览量: 393
- PDF下载量: 81
- 被引次数: 20
