利用双差层析成像方法反演青藏高原东南缘地壳速度结构

*东南缘地壳速度结构刘伟吴庆举’张风雪（中国北京100081中国地震局地球物理研究所）摘要本文利用云南及周边区域地震台网2010-2016年记录到的近震资料.采用双差层析成像 方法进行地震重定位并获得了青藏高原东南缘的维地壳速度结构•结果显示：巫定位后的震源

位置精度得到明显提高，震源主要分布于20 km深度以上的中上地壳；地戾分布与速度结构存在 一定的相关性，大多数地震发生在中上地壳的低速异常区内以及高、低速异常区域之间；研究区 上地壳速度结构存在明显的横向不均匀性.其速度异常与地表地形及地质特征密切相关；中下地

壳分布着两条主要的低速带，一条沿着安宁河断裂、小江断裂分布在川滇菱形地块的东侧；另一

条主要分布在川西北次级地块内，并穿过痂江断裂向南延伸，推测这两条低速带町能是青藏高原

中下地壳物质向南逃逸的两条通道.关键词双差层析成像地壳速度结构青藏高原东南缘重定位地壳流doi： 10.11939/jass.20180083 中图分类号:P315.3* 1 文献标识码：ACrustal structure of southeastern Tibetan Plateau inferred from double-difference tomographyLiu Wei Wu Qingju* Zhang Fengxue（Institute of Geophysics, China Earthquake Administration, Beijing 100081, China）Abstract： We applied the double-difference tomography to relocate seismic events and deter­

mined the crustal structure of the southeastern Tibetan Plateau from the inversion of the local

seismic data recorded by regional networks of Yunnan and its surrounding regions from 2010 to 2016・ Our results demonstrate that the accuracy of events relocation is significantly improved, and the events are mainly distributed in upper-mid crust at the depths shallower than 20 km. The locations of events are closely related to the velocity structure of upper-mid crust. Most of the earthquakes occurred in some regions with low-velocity anomaly or in the regions between low- velocity and high-velocity anomalies. Tomography results also show that the velocity structure of upper crust exhibits apparent horizontal heterogeneities, which is largely consistent with sur­face geological and topographic features. In the mid-lower crust, there are two main NS-trend- ing low-velocity zones, one is located in the east of Sichuan-Yunnan diamond block along the*基金项目国家I'l然科学基金（41674094）资助.收稿日期 2018-06-13收到初稿.2018-09-30决定采用修改稿.A 通信作者 e-mail： wuqj@cea-igp.ac.cn156地 震 学 报41卷trace of Anninghe and Xiaojiang faults, and the other is mainly located in the northwestern Sichuan sub-block and extends to the south across the Lijiang fault. These low-velocity zones are deduced to be the two crustal channels for the southward escape of crustal matter in the middle and lower crust of Tibetan Plateau.Key words： double difference tomography； crustal velocity structure； southeastern Tibetan

Plateau； relocation； crustal flow引言印度板块与欧亚板块的碰撞导致青藏高原的快速隆升，而强烈的挤圧作用对青藏高原 周边地块的地质地貌及构造演化也造成了巨大影响.青藏高原东南缘由于地处高原的边界 地区，地壳高度破碎，深大断裂发育，地震活动频繁且剧烈，在过去的30年里发生M6.0以

上地震百余次，是现今地壳形变与地震活动最为强烈的地区之一.前人为了解释青藏高原的 抬升过程与形变机制提出了刚性地块挤出模型(Tapponier, Molnar, 1976； Tapponnier et al,

1982； Avouac, Tapponnier, 1993)和地壳流模型(Royden et al, 1997; Clark, Royden, 2000),这

两种模型均认为青藏高原东南缘是高原地壳物质转移的重要场所与通道.因此，研究该地区 的地壳结构，对进一步厘清地震孕育环境以及青藏高原构造动力学过程等均具有重要意义.青藏高原东南缘地区因其独特的构造背景以及频繁的地震活动一直是研究的热点，近

年来众多研究人员利用多种地球物理手段开展了大量研究，揭示了该区域地壳结构的基本 特征，例如：体波层析成像的结果显示，研究区的中下地壳存在大范围的低速异常，多数强 震发生在中上地壳内的高速异常区以及高、低速异常转换带上(王椿铺等，2002；韦伟等，

2010)；面波与背景噪声层析成像的结果显示，研究区内壳内低速层的分布受断裂与构造边

界的限制(Yao etal, 2008),强震大多发生在中上地壳内的低速异常区以及高、低速异常分界

面上(王琼，高原，2014；郑定昌等，2014；潘佳铁等，2015);接收函数的研究结果则显示研 究区内S波速度结构具有很强的横向不均匀性，低速层分布的深度和范围差别很大(吴建平

等，2001；胡家富等，2003;李永华等，2009).综合以上研究成果，在青藏高原东南缘的中下 地壳存在低速层已得到公认，但是低速层的位置与分布形态仍然存在较大争议，地震分布与 速度结构之间的关系也尚未明确，因此仍有必要继续开展该区域深部结构的研究工作.为此，本文收集整理了云南及周边地区81个固定台站的观测数据，拟采用双差层析成 像方法(Zhang, Thurber, 2003, 2006)进行地震重定位，获取青藏高原东南缘的三维地壳速度

结构，分析地震分布与速度结构之间的关系，并试图结合地质构造背景来讨论速度结构特征 及其包括的动力学意义，以期为该区的动力学模型提供更多地震学证据的支持.1 数据与方法1.1数据本文采用的数据源于中国地震台网中心提供的2010-2016年期间云南及周边地区(图1)

81个固定台站记录到的初至Pg波震相报告资料.为了保证反演结果的可靠性，我们对原始 地震数据进行严格的筛选，仅选取观测报告中Ml》2. 0且走时残差W0.5s的震相数据，并根

据时距曲线拟合的方法剔除初始数据中存在明显异常的走时信息，即图2中两条绿线所括区 间以外的走时.双差层析成像方法需要利用事件对之间的相对走时数据，而满足一定条件的2期刘伟等：利用双差层析成像方法反演青藏高原东南缘地壳速度结构97° E 98°

99°

157100° 101 102° 103° 104° 105°

10630°N图1研究区域构造背景及台站分布戸：怒江断裂；△：澜沧汀断裂；Fy.南汀河断裂；F4：无量山断裂；Fs：金沙江断裂；F6：红河断裂；

F=：丽江断裂；F8；程海断裂；F9：绿汁江断裂；F10：安宁河断裂；F11：则木河断裂；F12：小江断裂Fig. 1 Regional tectonic settings and distribution of seismic stations in the studied areaF[： Nujiang fault； F2： Lancangjiang fhult； Nantinghe fault； F4： Wuliangshan fault； F5： Jinshajiang

fault； F6： Honghe fault； F7： Lijiang fault； F8： Chenghai fault； F9： Liizhijiang fault；F[o： Anninghe fault； Zemuhe fault； Fi2： Xiaojiang fault两个事件才能够组成事件对.本文将事件对 的参数选取标准设定为：事件对与台站的最

大间距为800 km.事件对之间的最大间距为 50 km,每个事件最多可以与50个事件组成 事件对，每个事件对所需震相数的最小值与

最大值分別为8和30.经过筛选，最终获得

1万3 135个地震事件的9万7 908条绝对到 时数据和209万1 189条相对到时数据参与反

演，射线路径分布如图3所示.可见射线均 匀地覆盖了整个研究区域.可以得到较好的 反演结果.1.2双差层析成像方法双差层析成像方法是Zhang和Thurber (2003, 2006)在双差定位方法(Waldhauser,

图2 时距曲线拟合图Fig. 2 The diagram for fitting travel timeEllsworth, 2000)的基础上发展而来，该方法with epicentral distance同时使用绝对走时数据和相对走时数据进行三维速度结构与震源位置联合反演，在反演过程15897° E 98°

99°

地 震 学 报41卷100° 101° 102° 103° 104° 105°

106°30°N中，先赋予绝对走时数据较高的权重以 建立整个区域内的三维速度结构；几次

29°28°27°迭代后赋予相对走时数据较高的权重来

更好地约束速度结构与震源位置（肖卓, 高原，2017）.由于考虑了介质速度结构

的空间变化，该方法克服了双差定位对 台站到事件对之间速度恒定的假设，所

26°25°24°得定位结果精度更高（王长在等，2013）; 而相对走时数据的加入，与传统的走时

层析成像方法相比可以提高速度结构反 演的精度，因此能够揭示更多的精细结 构信息（Zhang, Thurber, 2003）.1.3参数选择23°22°反演所选用的初始模型参考了研究

图3射线路径分布图区域内人工地震测深的结果（Zhang, Wang, 2009）,并根据数据的残差分布适 当调整而来（表1）.水平方向上采用0.5°X

Fig. 3 Distribution of ray paths表1 P波初始速度模型Table 1 Initial model of P wave velocity深度/km0.5。的网格间隔.竖直方向上的网格节点 分别位于0, 5, 10, 20, 30. 40 km.针对 阻尼最小二乘问题，双差层析成像方法

vp/( kms* 1)5.505.896.02深度/kmvP/(km-s 1)020306.236.506.905采用带阻尼的最小二乘正交分解（least

1040squares QR factorization,简写为 LSQR）算

法，以总走时残差的2范数为目标函数进行迭代求解.由于平滑因子和阻尼参数的大小对反演结果的稳定性有较大影响.因此在反演之前需要对不同平滑因子和阻尼参数的组合进行权衡分析（Eberhart-Phillips. 1986；王小娜 等，2015）,选取模型方差变化较小而数据方差显著降低时所对应的参数组合为最优值参与1.00.90.80.70.60.50.40.30.20.1解的归一化范数解的归-化范数图4利用L曲线法所选的最优平滑因子(a)和阻尼参数(b)Fig. 4 The optimum smoothing parameter (a) and damping parameter (b) selected by L curve method2期刘伟等：利用双差层析成像方法反演青藏高原东南缘地壳速度结构159最终反演.本文利用L曲线法搜索最优参数值，设定平滑因子的搜索范围为1-600,阻尼参 数的搜索范围为10-1 000,最终选取的最优平滑因子为40,阻尼参数为400 （图4）.1.4分辨率测试分辨率测试选用棋盘格测试方法（Spakmanercf/, 1993）.该方法首先在初始模型的基础

上添加正负相间的扰动合成理论棋盘模型，然后利用棋盘模型计算理论走时数据充当实际 观测数据并反演模型参数，最后观察反演结果对理论棋盘模型的恢复程度.若在…些区域内 恢复程度良好，则说明这些区域的结果是可信的.本文以0・5。乂0・5。网格和±5%高低速异常 的棋盘格模型作为理论模型正演计算理论走时数据，然后利用该数据和实际反演所用的初

始速度模型进行反演.图5为棋盘格测试的恢复结果，图像显示在地壳5, 10, 20, 30 km深

度上，研究区内大部分区域均得到良好的恢复•证明反演结果相对可靠.98 E 100

102°

104’100°2')N27°252321°l°E 10()102° 104°00°29° N27°252'21-5%图5不同深度上剖面的棋盘格测试结果Fig. 5 The checkboard resolution test at different depths2 主要结果双差层析成像方法利用地震事件对的相对走时数据来提高震源间相对位置的精度，得

到的结果更加准确.重定位前后的地震分布如图6所示.重定位前地震集中分布在5—10 km

范围内.重定位后浅层地震数目增多，地震主要分布在20 km深度以上的中上地壳（图7）； 而且重定位后所有地震事件的走时残差均方根显著减小，由0.380 s降至0.163 s.

160地 震 学 报41卷图6 重定位前(a)、后(b)的地震空间分布变化Fig. 6 Spatial variation of seismic events before (a) and after (b) relocation图7重定位前(a)、后(b)的震源深度变化Fig. 7 Focal depth variation before (a) and after (b) relocation图8给出了 5, 10, 20, 30 km深度处水平方向的P波速度扰动分布，可见在上地壳5—10 km 的深度范围内，速度结构呈现明显的横向不均匀性，其分布特征与地表的地质构造相关.5km

深度上(图8a),西昌东部、大理地区、楚雄地区、景谷一思茅地区均表现为低速异常，这些 地区处于不同的盆地内部，表明这些盆地内具有较厚的中生代沉积.在川西北次级地块内显 示的低速特征可能与该区域的复理石杂岩沉积有关(Yin. Harrison. 2000；潘佳铁等，2015).

熊绍柏等(1993)得到的丽江一者海深地震剖面显示丽江至永胜之间为低速异常，攀枝花以东 至安宁河断裂为高速异常.安宁河断裂至小江断裂之间为低速异常，本文的成像结果与之基

本吻合.小江断裂的东川地区所表现的低速异常可能是受到了深部热作用的影响.高速异常

多集中在攀枝花地区、昭通地区以及南汀河断裂与澜沧江断裂的交会处.2期98° E 99°

刘伟等：利川双差层析成像方法反演青藏高原东南缘地壳速度结构100°

161101° 102° 103° 104° 105'98° E 99° 100° 101° 102° 103° 104° 105'21----------1___I------1------1-------5% -4% 一3% -2% -1% 0 1% 2% 3% 4% 5%dvp/vp地震图8不同深度力处P波速度扰动和地震分布图(粗实线AA', BB\\ CC'为地震剖面的位置)Fig. 8 P wave velocity disturbance and earthquake distribution at different depth h

The thick lines AA ； BB', and CC' represent the position of the seismic sections(a) h = 5 km； (b) h = 10 km； (c) /» = 20 km； (d)方=30 km；10 km深度上(图8b),高低速异常分布主要受到断裂带的控制，例如：红河断裂、南汀

河断裂、澜沧江断裂等整体处于低速异常；小江断裂两侧差异明显，其东侧表现为高速异

常，西侧的昆明一玉溪表现为明显的低速异常，这与吴建平等(2013)的层析成像结果相近. 腾冲地区的低速异常反映了新生代火山的活动特征，而攀枝花附近的区域仍然表现为高速

异常.在中下地壳20-30 km的速度扰动图上出现大范围的低速异常.20km深度上(图8c), 红河断裂的东西两侧速度结构差异明显.低速异常主要分布在红河断裂东侧的两条带状区

域内.其中：一条低速带沿着安宁河断裂、小江断裂分布在川滇菱形地块东侧；另一条主要

162地 震 学 报41卷分布在川西北次级地块内.并穿过丽江断裂向南延伸.两条低速带之间的攀枝花地区表现为 高速异常.红河断裂的西侧，除了临沧和腾冲地区表现为低速异常外，其余大部分地区均表

现为高速异常.30 km深度上(图8d),高低速异常分布与20 km深度相近，低速带的分布范 围较20 km深度的更广.前人的近震层析成像结果同样显示，在研究区的中下地壳内存在大 范围的低速带(王椿铺等，2002； Huang et al, 2002；韦伟等，2010)，但是受观测资料的限制，

不同研究显示的低速带位置及形态相差较大.3 讨论与结论3.1地震分布与速度结构为了讨论研究区内的地震分布与速度结构之间的关系，我们将重定位后的震源位置分 别投影至不同深度的速度扰动剖面上，即将震源深度处于0-7.5 km的地震投影至5 km深度 的剖面上，7.5—15 km的地震投影至10 km深度的剖面上，15—25 km的地震投影至20 km深

度的剖面上，大于25 km的地震投影至30 km深度的剖面上.从图8可以看出，区域内的地震

活动与速度结构的变化明显相关，大多数地震发生在中上地壳的低速异常区内以及高、低速 异常区域之间.钱晓东等(2011)统计分析了研究区域内的震源机制解参数认为，该地区主要

受到3个方向的应力作用：印缅地块北东向的作用力、川滇菱形地块东南向的作用力以及华 南地块北北西向的作用力.由于低速异常区域的强度相对较弱，在这种横向挤压的应力场作 用下，这些低速异常区域更容易发生破裂导致地震.一些研究认为地震活动集中分布在中上

地壳的低速区域内，可能与断层或微裂隙中存在流体有关(Zhao可刃，1996； 2000；齐诚等， 2006；李永华等，2014；潘佳铁等，2015).在研究区域内.地壳高度破碎，断层广泛发育，流

体的存在会降低断层的强度从而造成孕震区的弱化，使地震更容易发生.另外，高、低速异

常交界的区域往往是介质物性发生变化的区域，由于结构和应力场的变化，这些相对脆弱的 区域更有利于能量的积累与释放进而引发地震.研究区的中下地壳存在大范围的低速异常

区域，黄金莉等(2001)认为下部的低速区域是上部发生强烈地震的重要构造背景，王椿籀等 (2002)认为中下地壳低速区域的存在有利于应力集中于上部地壳从而导致地震的发生.将距离剖面两侧25 km范围内重定位后的地震震中投影至图8c所示的纵向速度剖面图

中，以进一步分析地壳速度结构与地震分布的关系，结果如图9所示，可以看出，重定位后 的地震沿着断裂带呈现明显的带状或簇状分布特征，地震大多发生于低速异常区域内以及

高、低速异常区域之间.青藏高原东南缘地区的地壳结构复杂，地震活动频繁且剧烈.研究 分析地震分布与速度结构的对应关系，对于进一步研究地震孕育环境等具有重要意义.3.2低速异常带与地壳流Royden等( 1997)提出的地壳流模型认为，在青藏高原东缘中下地壳存在一个黏性流动

层，青藏高原的物质沿此向外流出.这一模型很好地解释了青藏高原东缘的地形变化以及上 地壳无明显缩短的现象(Royden et a/, 1997； Clark. Royden, 2000),获得了较多研究人员的支

持，而寻找地壳流存在的证据则成为深部地球物理学必须面对的一个科学问题.判断地壳流的主要依据是偏低的地震波速度，Yao等(2008)利用背景噪声成像的结果显

示.研究区壳内低速层分布十分复杂，认为中下地壳物质不会发生大规模的流动，而是被断 裂及构造边界限制在一定的区域内；韦伟等(2010)的近震P波层析成像结果显示，研究区中 下地壳深度处的川滇菱形地块内低速异常广泛发育.认为该低速异常很可能就是青藏高原2期刘伟等：利用双差层析成像方法反演青藏高原东南缘地壳速度结构163丽江断裂 安丫河断裂4.53.L.0S5W.50

O

E<

102030BB'99° E100°澜沧江断裂 无駐山断裂4.5 r3.0 [红河断裂01o

E<«嶷

2o

30 I99° E100° 1015.4

5.6

102°6.2 (km»s ')103°6.6

104'5.8 6.0 6.4 6.8 7.0图9 P波速度纵向剖面与地震分布(剖面位置见图8)Fig. 9 P wave velocity vertical profiles (their locations are shown in Fig. 8) and earthquake distribution下地壳流的通道；Liu等(2014)利用背景噪声与接收函数联合反演的结果显示，松潘一甘孜 地块内存在大范围的低波速、高泊松比区域，认为是高原内部物质流出所致.近年来一些研究表明，青藏高原东南缘可能存在两条地壳流通道(Bai etal, 2010； Bao et al. 2015). Bai等(2010)的大地电磁观测结果显示，青藏高原东部存在两条壳内高导带，-条环绕喜马拉雅东构造结向南转折通过滇西南地区，另一条沿鲜水河断裂、安宁河断裂、 小江断裂向南延伸至印支地块，且这两条高导带被认为是两条独立的地壳流通道；Bao等

(2015)关于面波频散与接收函数联合反演的结果表明，研究区壳内20-40 km深度存在两条 S波低速带，且这两条低速带的分布较好地对应了上述高导带的位置.然而，另外一些研究

164地 震 学 报41卷人员采用面波频散与接收函数联合反演的方法得到的结果所显示的低速带位置和形态分布 均存在差异.例如：郑晨等（2016）的结果显示两条低速带仅向南延伸至24呗左右，并未穿过 红河断裂；而Li等（2016）的结果表明两条低速带穿过红河断裂，且最终交会于云南南部.这

些结果的不一致一定程度上说明了研究区下方低速体的形态分布和位置较为复杂.本文结果显示，在研究区的中下地壳同样存在着两条低速带，其中一条沿安宁河断裂、

小江断裂分布在川滇菱形地块的东侧，另一条主要分布在川西北次级地块内.并穿过丽江断

裂向南延伸，这两条低速带在前人的研究中均有所体现（Bao et al, 2015;郑晨等，2016； Lie/ al. 2016）.在低速带分布的位置同样存在高衰减（周龙泉等.2009； Zhao et al, 2013）、高热流 （Hu etal. 2000）、高*/（ Sun et a/, 2014）等异常，这表明壳内物质存在部分熔融，其在重力

的作用下M发生塑性流动.此外，王琼等（2015）基于背景噪声的瑞雷波方位各向异性的研究

显示，在中下地壳深度范围内，川西北次级地块和小江断裂附近表现为强各向异性，本文中 低速带的延展方向与其快波优势方向基本一致.进一步证明了处于低速带分布位置的壳内

物质可以发生塑性流动，造成矿物的定向排列从而导致各向异性现象的产生.此外，这两条 低速带并未完全连通，而是被攀枝花及周边区域的高速异常体分隔.地质研究表明，中晚二 叠纪时期发生了一次巨大的玄武岩喷发事件，形成了著名的峨眉山大火成岩省（徐义刚，钟 孙霖，2001；张招崇等，2006）,这次巨大的玄武岩喷发事件很可能是由古地幔柱的活动作用

所引起，攀枝花地区位于玄武岩喷发前地幔柱活动的核心部位（H e et al, 2003）.地幔柱活动 导致的基性和超基性岩体残留在地壳内部，表现为高速异常（吴建平等，2013；徐涛等，

2015）,这些髙速物质可能对青藏高原物质向南转移起到了一定的阻挡作用.在图9所示的纵向速度剖面中均能观测到壳内低速带的存在，且低速带主要沿着研究区

内绿汁江断裂、小江断裂等走滑断裂分布.Leloup等（1999）认为走滑断裂相对运动产生的热 量能够降低中下地壳的黏度和地震波速度，从而形成低速层；而低黏度的低速层也易与上地 壳发生相对运动.本文地震重定位的结果显示少量地震发生在低速带内部，绝大多数地震发 生在低速带的边界区域.表明壳内低速带的存在可能会促使断层发生运动而引发地震.另 外，沿小江断裂分布的低速带的延伸方向在26呗发生偏转.由与断裂平行的近南北向逐渐

转变为西南向.与全球定位系统显示的速度场方向（相对于华南地块）一致，这表明低速带的

分布可能对I ：地壳运动产生影响.张培丧（2008）认为造成该地区地壳发生顺时针旋转的深部 动力学机制可能是中下地壳软弱层物质的流动.结合前人的研究结果，本文认为青藏高原东 南缘存在两条低速带，它们可能是青藏高原物质向南逃逸的两条通道.3.3结论

本文采用双差层析成像方法，利用2010-2016年中国地震台网中心提供的云南及周边 地区Pg波走时资料，进行地震重定位并获得了青藏高原东南缘的三维地壳速度结构，主要

结果如下：1） 经重定位后，走时残差明显减小，地震主要分布在20 km深度以上的中上地壳；地震 分布与速度结构存在一定关系，大多数地震发生在中上地壳的低速异常区内以及高、低速异 常区域之间；2） 层析成像结果显示，研究区上地壳速度结构存在明显的横向不均匀性.速度分布特 征与地衣的地质构造相关；3） 研究区的中下地壳存在两条低速带，一条沿安宁河断裂、小江断裂分布在川滇菱形

2期刘伟等：利川双差层析成像方法反演青藏高原东南缘地壳速度结构165地块的东侧；另一条主要分布在川西北次级地块内，并穿过丽江断裂向南延伸.两条低速带

被攀枝花及周边区域的高速异常体分隔，它们可能是青藏高原中下地壳物质向南逃逸的两 条通道.中国地震台网中心提供了震相数据，中国科学技术大学张海江教授提供了 tomoDD-SE 程序，两位评审专家在稿件修改过程中提出了宝贵意见，本文图件绘制采用GMT软件和

Matlab软件，作者在此一并表示衷心的感谢.参考文献胡家富.苏沧協.朱雄关.陈诜.2003.云南的地壳S波速度泊松比结构及其蠢义[J].中国科学：Dtt. 33(8)： 714-722.

Hu J F, Su Y J, Zhu X G, Chen Y. 2005. S-wave velocity and Poisson's ratio structure of crust in Yunnan and its implication[J].Science in China： Series D, 48(2) ： 210 218.黄金莉.赵大鹏.郑斯华.2001.川滇活动构造区地震层析成像[J].地球物理学报.44(增刊1)： 127-135.Huang J L, Zhao D P. Zheng S H. 2001. Seismic tomography of the Sichuan-Yunnan active tectonic region[J]. Chinese Journal

of Geophysics, 44( SI)： 127 - 135 (in Chinese).李永华.吴庆举，田小波，张瑞青，潘佳铁，曾融生.2009.用接收函数方法研究云南及其邻区地売上地幔结构[J].地球物

理学报,52( 1)： 67-80.Li Y H, Wu Q J. Tian X B, Zhang R Q. Pan J T, Zeng R S. 2009. Crustal structure in the Yunnan region determined by model­

ing receiver functions[J], Chinese Journal of Geophysics, 52( 1 ): 67-80 (in Chinese).李永华，徐小明，张恩会，高家乙.2014.青藏高原东南缘地壳结构及厶南鲁甸、最谷地農深部孕丧环境[J].地虎地质，

36(4) : 1204 -1216.Li Y H, Xu X M, Zhang E H. Gao J Y. 2014. Three-dimensional crust structure beneath SE Tibetan Plateau and its seismotcctonic

implications for the Ludian and Jinggu earthquakes[J]. Seismology and Geology, 36( 4) : 1 204 1216 (in Chinese).潘佳铁，李永华，吴庆举，丁志峰.2015.青藏高原东南部地区瑞雷波相速度层析成像[J].地球物理学报.58( 1 1 )：

3993-4006.Pan J T, Li Y H, Wu Q J. Ding Z F. 2015. Phase velocity maps of Rayleigh waves in the southeast Tibetan Plateau卩].Chinese

Journal of Geophysics, 58( 11)： 3993 4006 (in Chinese).齐诚.赵大鹏.陈顒.陈棋福.E宝眸.2006.首都圈地区地说P波和S波〔维速度结构及其与大地隠的关系[J].地球物理

学报,49(3)： 805-815.Qi C, Zhao D P, Chen Y, Chen Q F, Wang B S. 2006. 3-D P and S wave velocity structures and their relationship to strong

earthquakes in the Chinese capital rcgion[J]. Chinese Journal of Geophysics, 49( 3): 805-815 (in Chinese).钱晓东，秦嘉政.刘丽芳.2011.云南地区现代构造应力场研究[J].地震地质.33( 1)： 91-106.Qian X D. Qin J Z, Liu L F. 201 1. Study on recent tectonic stress field in Yunnan region[J]. Seismology and Geology, 33( 1 ):

91-106 (in Chinese).王椿铺.Mooney W D.王溪莉.吴建平，楼海.王E. 2002.川滇地|x•地壳地幔二维速度结构研究[J].地虎学报.

24( 1)： 1-16.Wang C Y. Mooney W D, Wang X L, Wu J P. Lou H. Wang F. 2002. Study on 3-D velocity structure of crust and upper mantle

in Sichuan-Yunnan region, China[J]. Acta Seismologica Sinica. 24( 1) : 116 (in Chinese).王长在.吴建平.房立华，王未来.2013.玉树地虎艰源区速度结构9余丧分布的关系[J].地球物理学报.56( 12)：

4072-4083.Wang C Z. Wu J P, Fang L H . Wang W L. 2013. The relationship between wave velocity structure around Yushu earthquake

source region and the distribution of aftershocks[J]. Chinese Journal of Geophysics. 56( 12)： 4072 4083 (in Chinese).£琼，高原.2014.青藏东南缘背景噪尙的瑞利波相速度层析成像及强震活动[J].中国科学：地球科学，44(11)：

2440-2450.Wang Q. Gao Y. 2014. Rayleigh wave phase velocity tomography and strong earthquake activity on the southeastern front of the166地 震 学 报41卷Tibetan Plateau[J]. Science China Earth Sciences, 57( 10) : 2532-2542.F琼，高原，石E涛.2015.青藏高原东南缘基于背景噪声的Rayleigh面波方位各向异性研究[J].地球物理学报.58( I I)：

4068-4078.Wang Q, Gao Y, Shi Y T. 2015. Rayleigh wave azimuthal anisotropy on the southeastern front of the Tibetan Plateau from seis­

mic ambient noise[J]. Chinese Journal of Geophysics. 58( 11)： 4068-4078 ( in Chinese).匸小娜，于湘伟，章文波.2015.芦山震区地壳：.维P波速度精细结构及地震殖定位研究[J].地球物理学报，58(4)：

1179-1193.Wang X N, Yu X W, Zhang W B. 2015. 3D P-wave velocity structure of the crust and relocation of earthquakes in the Lushan

source area[J]. Chinese Journal of Geophysics. 58( 4) ： 1 1 79- 1 1 93 (in Chinese).韦伟.孙若昧.石耀霖.2010.青藏高原东南缘地震层析成像及汶川地震成因探讨[J].中国科学：地球科学，40( 7)：

831-839.Wei W, Sun R M. Shi Y L. 2010. P-wave tomographic images beneath southeastern Tibet： Investigating the mechanism of the

2008 Wenchuan earthquake[JJ. Science China Earth Sciences, 53( 9) : 1252-1259.吴建平，明跃红，匸椿铺.2001.云南数字地震台站下方的S波速度结构研究[J].地球物理学报，44(2)： 228-237.Wu J P, Ming Y P, Wang C Y. 2001. The S wave velocity structure beneath digital seismic stations of Yunnan Province inferred

from teleseismic receiver function modeling[J]. Chinese Journal of Geophysics, 44( 2): 228-237 (in Chinese).吴建平，杨婷，王未来，明跃红，张天中.2013.小江断裂带周边地区〔维P波速度结构及具构造意义[J].地球物理学报，

56(7): 2257-2267.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).肖卓，高原.2017.利用双差成像方法反演青藏高原东北缘及其邻区地壳速度结构[J].地球物理学报.60(6)： 2213-2225.Xiao Z, Gao Y. 2017. Crustal velocity structure beneath the northeastern Tibetan Plateau and adjacent regions derived from

double difference tomography[J]. Chinese Journal of Geophysics, 60( 6) : 2213-2225 (in Chinese).熊绍柏，郑晔.尹周勋，曾晓献，全幼黎.孙克忠.1993.丽江一攀枝花一者海地带二维地壳结构及其构造意义[J].地球物

理学报.36(4)： 434-444.Xiong S B, Zheng Y. Yin Z X, Zeng X X, Quan Y L. Sun K 乙 1993. The 2-D structure and it's implications of the crust in

the Lijiang-Panzhihua-Zhehai region[J]. Chinese Journal of Geophysics, 36( 4) : 434-444 (in Chinese).徐涛，张忠杰，刘宝峰，陈费，张明辉.田小波，徐义刚，滕吉文.2015.峨眉山大火成岩省地壳速度结构与書地幔柱活动

遗迹：来自丽江一清镇宽角地震资料的约束[J].中国科学：地球科学.45(5)： 561-576.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) : 1 133-1 147.徐义刚，仲孙霖.2001.峨眉山大火成岩省：地幔柱活动的证据及其熔融条件[J].地球化学.30(1)： 1-9.Xu Y G, Zhong S L. 2001. The Emeishan large igneous province： Evidence for mantle plume activity and melting conditions[J].

Geochimica. 30( 1) : 1 9 (in Chinese).张培震.2008.青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程[J].中国科学：D辑，38(9)：

1041-1056.Zhang P Z. 2008. The present day tectonic deformation, strain partitioning and deep dynamic process of the western Sichuan

region along eastern margin of Qinghai-Tibet PlateaufJ]. Science in China'. Series D, 38( 9) : 1041-1056 (in Chinese).弓长招崇，Mahoney J J.王福生，赵莉，艾羽.杨铁铮.2006.峨眉山大火成岩省西部苦橄岩及其共生玄武岩的地球化学：地

幔柱头部熔融的证据[J].岩石学报,22(6)： 1538-1552.Zhang Z C, Mahoney J J, Wang F S, Zhao L, Ai Y. Yang T 乙 2006. Geochemistry of picritic and associated basalt flows of the

western Emeishan flood basalt province, China： Evidence for a plume-head origin[J]. Acta Petrologica Sinica, 22( 6 ): 1538-1552 (in Chinese).郑晨，丁志峰.宋晓东.2016.利用面波频散与接收函数联合反演青藏高原东南缘地壳上地幔速度结构[J].地球物理学

报,59(9)： 3223-3236.Zheng C, Ding Z F, Song X D. 2016. Joint inversion of surface wave dispersion and receiver functions for crustal and uppermost

2期刘伟等：利用双差层析成像方法反演青藏高原东南缘地壳速度结构167mantle structure in southeast Tibetan Plateau[J]. Chinese Journal of Geophysics, 59( 9): 3223-3236 (in Chinese).郑定昌，盖增喜，杨润海.闵照旭.唐有彩，姜明明.庞卫东.2014.云南地区背景噪声层析成像[J].地震学报.36(4)：

602-614.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-6 1 4 (in Chinese).周龙泉，刘杰，苏有锦，马宏生，周俊杰.2009.利丿US波高频衰减参数对厶•南地区地壳0值成像[J].地球物理学报.

52(6) : 1500-1507.Zhou L Q. Liu J, Su Y J. Ma H S, Zhou J J. 2009. Tomography for Q of Yunnan region from high-frequency attenuation of S

wave[J]. Chinese Journal of Geophysics, 52( 6): 1 500- 1507 (in Chinese).Avouac J P, Tapponnicr P. 1 993. Kinematic model of active deformation in central Asia[J]. Geophys Res Lett, 20( 10 ):

895-898.Bai D H , Unsworth M J, Meju M A, Ma X B, Teng J W, Kong X R, Sun Y. Sun J, Wang L F, Jiang C S, Zhao C P, Xiao P

F, Liu M. 2010. Crustal deformation of the eastern Tibetan Plateau revealed by magnetotelluric imagingfJJ. Nat Geosci,

3(5) ： 358-362.Bao X W, Sun X X, Xu M J, Eaton D W, Song X D, Wang L S, Ding Z F, Mi N, Li H, Yu D Y, Huang Z C, Wang P. 2015.

Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver func-

tions[J]. Earth Planet Sci Lett、415： 1 6-24.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.Eberhart-Phillips D. 1986. Three-dimensional velocity structure in northern California Coast Ranges from inversion of local earth­

quake arrival times[J]. Bull Seismol Soc Am, 76( 4) : 1025-1052.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 basaits[J], Earth Planet Sci Lett、213( 3/4) : 391—405.Hu S B. He L J. Wang J Y. 2000. Heat flow in the continental area of China： A new data set[J]. Earth Planet Sci Left、179( 2):

407-419.Huang J L, Zhao D P, Zheng S H. 2002. Lithospheric structure and its relationship to seismic and volcanic activity in southwest

China[J]. J Geophys Res, 107(B10)： ESE 13-1-ESE 13-14.Leloup P H, Ricard Y, Battaglia J, Lacassin R. 1999. Shear heating in continental strike-slip shear zones： Model and field

examples[J]. Geophys J Int, 136( 1): 19-40.Li M K. Zhang S X, Wang F, Wu T F, Qin W B. 2016. Crustal and upper-mantle structure of the southeastern Tibetan Plateau

from joint ana】ysis of surface wave dispersion and receiver functions[J]. J Asian Earth Sci, 117( 9) : 52-63.Liu Q Y. van tier Hilst R D, Li Y, Yao H J, Chen J H, Guo B, Qi S H, Wang J, Huang H, Li S C. 2014. Eastward expansion

of the Tibetan Plateau by crustal flow and strain partitioning across faults[J]. Nat Geosci, 7( 5) ： 361-365.Royden L H, Burchfiel B C, King R W, Wang E, Chen Z L, Shen F, Liu Y P. 1997. Surface deformation and lower crustal

flow in eastern Tibet[J]. Science, 276( 5313)： 788-790.Spakman W, van der Lee S, van der Hilst R. 1993. Travel-time tomography of the European-Mediterranean mantle down to 1 400 km

[J]. Phys Earth Planet Inter, 79( 1/2) ： 3-74.Sun X X. Bao X W. Xu M J, Eaton D W, Song X D, Wang L S, Ding Z F. Mi N, Yu D Y, Li H. 2014. Crustal structure be­

neath SE Tibet from joint analysis of receiver functions and Rayleigh wave dispersion[J]. Geophys Res Lett. 41 ( 5 ):

1479-1484.Tapponnier P. Molnar P. 1976. Slip-line field theory and large-scale continental tectonics[JJ. Nature. 264( 5584) : 3 1 9-324.Tapponnier P. Peltzer G. Le Dain A Y, Armijo R, Cobbold P. 1982. Propagating extrusion tectonics in Asia： New insights from

simple experiments with plasticinc[J]. Geology, 10(12)： 61 1-616.Waldhauser F. Ellsworth W L. 2000. A double-difference earthquake location algorithm： Method and application to the northern

Hayward fault, California[J]. Bull Seismol Soc Am, 90( 6): 1353-1368.Yao H J. Bcghein C. van der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-

station analysis： II. Crustal and upper-mantle structure[J]. Geophys J Int, 173( 1 ): 205-219.168地 震 学 报41卷Yin A, Harrison T M. 2000. Geologic evolution of the Himalayan-Tibetan orogen[J]. Ann Rev Earth Planet Sci, 28： 2 1 1-280.Zhang H J, Thurber C H. 2003. Double-difference tomography： The method and its application to the Hayward fault. California

[J]. Bull Seismol Soc Am, 93(5) ： 1875-1889.Zhang H J, Thurber C H. 2006. Development and applications of double-difference seismic tomography[J], Pure Appl Geophys,

163(2/3) ： 373-403.Zhang X. Wang Y H. 2009. Crustal and upper mantle velocity structure in Yunnan, Southwest China[J]. Tectonophysics,

471(3/4) : 1 71-185.Zhao D P. Kanamori H, Negishi H, Wiens D. 1996. Tomography of the source area of the 1995 Kobe earthquake： Evidence for

fluids at the hypocenter?[J]. Science, 274( 5294): 1891-1 894.Zhao D P, Ochi F, Hasegawa A, Yamamoto A. 2000. Evidence for the location and cause of large crustal earthquakes in

Japan [J]. J Geophys Res, 105( B6) ： 1 3579-13594.Zhao L F, Xie X B. He J K, Tian X B, Yao Z X. 2013. Crustal flow pattern beneath the Tibetan Plateau constrained by regional

Lg-wave Q tomography[J]. Earth Planet Sci Lett, 383( 4): 1 13-122.作者简仆王耀中国科学技术大 惠少兴陕西省地震局工 学地球和空间科学学院地

程师.2009年云南大学资

球物理学专业在读硕士研 源环境与地球科学学院地 究生.2015年吉林大学地 球物理学匕业毕业，获理 球探测科学与技术学院地

学学士学位；2012年中国 球物理学专业毕业，获理

地震局兰州地震研究所固 学学士学位.主要从事地体地球物理学专业毕业.壳速度结构模型评价等方面的研究.获理学硕士学位.现主耍从爭:地震监测方法.非天然地震口动识别及卫星数据等的应用研

究.陕西省地震学会会员.刘伟中国地震局第二 李君中国地震局第 监测中心助理工程师.2015

二监测中心助理工程师. 年河北大学物理科学与技 2010年中国地质大学（武

术学院应用物理学专•业毕 汉）地球物理与空间信息学 业,获工学学士学位;2018

院地球信息科学与技术专 年中国地震局地球物理研 业毕业.获工学学士学位； 究所地球探测与信息技术2016年中国地震局地震预&业毕业.获丁.学硕七学位.主要从审地震定位测研究所地球物理学专业毕业.获理学硕上学 与层析成像方面的研究.位.曾从事地球物理勘探作，现从书数■:地

震学相关研究工作.韩佳东中国地震局地球 崔子健中国地震局地

物理研究所地球探测与信 震预测研究所副研究员. 息技术专业在读硕士研究 2011年中国地震局地震预

生.2015年吉林大学地球探

测研究所固体地球物理学

测科学与技术学院勘查技 专业毕业.获理学硕七学 术与工程（应用地球物理） 位；2018年中国地震局地

Q业毕业.获工学学士学球物理研究所固体地球物位.主要从事地震监测1-J地震活动性方面的研理学专业毕业.获理学博七学位.现主要从书数究.中国地震学会会员.字地震学和地震活动性等方面的研究.I