青藏高原东缘的四川地区位于南北地震带中部，巴颜喀拉地块、华南地块里的四川盆地和川滇地块在区域内交汇，区域内分布有多条地震活断层，如龙门山断裂、鲜水河断裂等(如图 1).2008年5月12日，在龙门山断裂带上发生了汶川M8.0地震.时隔五年之后，2013年4月20日，龙门山断裂带南段发生了芦山M7.0地震.2017年8月8日，在岷江断裂与塔藏断裂交汇区发生了九寨沟M7.0地震.中国地震台网中心地震目录显示，自2010年01月01日至2017年12月31日，四川地区共发生M≥5.0地震17次，5.0＞M≥4.0地震大约有100多次.汶川M8.0地震的发生，对区域内构造应力场有调整作用，研究汶川震后中强地震活动和震源机制及区域应力场分布特征，分析中强地震深度分布有重要意义.

图 1

Fig. 1

图 1 研究区构造背景及地震台站和震中分布 图中黑色三角形为地震台站位置，字母为台站代码，红色圆点为M≥5.0地震位置，红色小五角星为芦山地震和九寨沟地震震中位置，红色大五角星为汶川地震位置. F1龙门山前山断裂；F2龙门山中央断裂；F3龙门山后山断裂；F4岷江断裂；F5塔藏断裂；F6虎牙断裂；F7雪山断裂；F8鲜水河断裂；F9峨边断裂；F10西河—美姑断裂；F11甘洛—竹核断裂；F12德格—乡城断裂；F13理塘—义敦断裂；F14金沙江断裂；F15德钦—中甸断裂；F16米亚罗断裂. Fig. 1 The tectonic setting and the distributions of the seismic stations and epicenters The black triangles indicate the seismic stations and the characters are station coda. The red dots indicate the earthquake with M≥5.0. The small red stars indicate Lushan earthquake and Jiuzhaigou earthquake. The big red star indicates Wenchuan earthquake. F1 Front-range of Longmenshan fault; F2 Central-range of Longmengshan fault; F3 Back-rang of Longmenshan fault; F4 Mingjiang fault; F5 Tazang fault; F6 Huya fault; F7 Xueshan fault; F8 Xianshuihe fault; F9 Ebian fault; F10 Xihe-Meigu fault; F11 Ganluo-Zhuhe fault; F12 Dêgê-Xiangcheng fault; F13 Litang-Yidun fault; F14 Jinshajiang fault; F15 Dêqên-Zhongdian fault; F16 Miyaluo fault.

中强地震的震源机制解和准确的震源深度有益于对区域应力特征、地震活动特性、发震断层结构等的理解.研究显示，震源机制解和震源深度分别是认识断层产状结构，确定发震断层，了解震源区应力状态，判定发震机理和后续发震趋势等的重要依据(Waldhauser and Ellsworth, 2000; 易桂喜等, 2016, 2017a; Rubin, 2002; Hauksson et al., 2012; 高原等, 1997; 赵博等，2013).震源深度定位的准确性一直是地震学研究的难题.传统的走时定位受到台站分布和区域地震波速度模型等的影响，在深度方向上的定位精度较低.经典的盖格法在线性化、最小二乘求解以及观测系统的误差分布上有局限性(许力生等, 2013a, b).特别是深度定位，受到定位方法、台站分布、区域速度结构的影响较大，本研究拟根据面波能量谱特征进行深度定位，以避免目标函数线性化、台站分布、区域速度结构等对深度的影响.本研究首先利用波形拟合反演法对四川地区17次M≥5.0地震和100多次5.0＞M≥4.0震源机制解进行反演(Herrmann, 2013)，然后根据震源机制解计算理论地震图，根据面波振幅谱特征，尝试对震源深度进行准确定位(Tsai and Aki, 1970; Fox et al., 2012; Heyburn et al., 2013).Tsai和Aki提出了利用地震面波振幅谱精确确定震源深度的方法, 并论证了方法的可行性(Tsai and Aki, 1970).地震面波能量与震源机制、观测方位和震源深度有密切关系.当震源机制与观测方位确定后，可以利用面波振幅谱与震源深度唯一依赖关系，确定震源深度.

四川地区中强地震大多分布在区域内主要构造块体的交汇区和活动断层附近.川西—青藏交汇区地震活动背景一直很强(Gao et al., 2000).如图 1，汶川M8.0地震、芦山M7.0地震和九寨沟M7.0地震均发生在巴颜喀拉块体和华南块体的交汇区.根据震源机制解资料和剪切波分裂参数研究，在巴颜喀拉地块东部、龙门山断裂带南段的主压应力方向为NWW—近EW向(Zhao et al., 2012; 赵博等，2013；高原等，2013).四川盆地—龙门山—巴颜喀拉一线的地震勘探剖面显示，该地区地壳厚度在龙门山两侧变化很大(滕吉文等，2008; 朱介寿，1984；Zhang et al., 2010).在2008—2017年期间，研究区内破坏性地震发生频繁.汶川地震和芦山地震的发震断层为龙门山断裂带，九寨沟地震发生在虎牙断裂、岷江断裂、塔藏断裂所围限的区域.四川地区南部为川滇块体和巴颜喀拉块体的交汇区.道孚M5.0地震、炉霍M5.4、康定M6.4和M5.9地震均发生在该区域的鲜水河断裂上，白玉M5.5地震、理塘M5.2地震、得荣M5.2地震则发生在德格—乡城断裂、理塘断裂、金沙江断裂和德钦—中甸断裂的交汇区.在华南块体和川滇块体的交汇区，发生过越西M5.2地震和金口河M5.0地震.

本研究收集了2010年1月1日—2017年12月31日，四川地区发生的中强地震的波形数据(图 2)，其中包括60多个四川区域地震台网和140多个国家地震台网波形数据.震源机制反演所用的波形数据为区域台网近场波形数据(震中距≤500 km)，震源深度确定所用的数据为国家台网远场波形数据(震中距≥500 km).

图 2

Fig. 2

图 2 研究使用的地震台站分布 蓝色三角形为国家地震台网台站，黑色三角形为区域地震台网台站. Fig. 2 The distribution of seismic networks The blue triangles indicate the national seismic networks and the black triangles indicate the regional seismic networks.

震源机制反演方法有多种，例如，P波初动法、振幅比法、波形拟合法等.本研究通过计算理论地震图，采用CPS程序进行波形拟合反演(Herrmann, 2013).对于M6.0以上的强震，可以利用远场波形记录和大尺度速度模型(如AK135模型)进行反演，但对于区域中等强度地震，如5.0≤M≤5.9的地震，远场波形记录信噪比不高，因此可以选择区域台网记录的近场波形进行反演.本研究利用四川省区域台网的波形记录，对2010年1月1日—2017年12月31日四川地区中强地震进行震源机制反演.M≥5.0地震参数见表 1(包括汶川地震主震震源参数)，所采用的速度模型为表 2(赵珠等，1997).

表 1 Table 1

表 1 地震参数 Table 1 The parameters of the earthquakes 序号 发震时刻 北京时间 经度 (°E) 纬度 (°N) 深度 /km 震级 /MS 地名 0* 2008-05-12 14:28:04 103.39 31.00 18 8.0 四川汶川 1 2010-04-28 04:22:27 101.45 30.60 8 5.0 四川道孚 2 2010-05-25 14:11:53 103.49 31.17 20 5.0 四川汶川 3 2011-04-10 17:02:42 100.80 31.28 10 5.4 四川炉霍 4 2011-11-01 05:58:15 105.30 32.60 6 5.2 四川青川 5 2013-01-18 20:42:50 99.40 30.95 15 5.5 四川白玉 6** 2013-04-20 08:02:47 102.97 30.29 17 7.0 四川芦山 7 2013-04-20 11:34:16 102.94 30.24 15 5.4 四川芦山 8 2013-04-21 04:53:44 103.05 30.35 27 5.4 四川芦山 9 2013-04-21 17:05:23 103.00 30.34 17 5.4 四川芦山 10 2013-08-28 04:44:52 99.33 28.20 9 5.2 四川得荣 11 2014-10-01 09:23:29 102.74 28.38 10 5.2 四川越西 12 2014-11-22 16:55:28 101.68 30.29 20 6.4 四川康定 13 2014-11-25 23:19:09 101.75 30.20 16 5.9 四川康定 14 2015-01-14 13:21:40 103.20 29.30 20 5.0 四川金口河 15 2016-09-23 00:47:13 99.60 30.08 19 5.2 四川理塘 16 2017-08-08 21:19:48 103.82 33.20 10 7.0 四川九寨沟 17 2017-09-30 14:14:37 105.05 32.25 10 5.4 四川青川 注：*汶川主震地震参数引自Zhao等(2011). **芦山主震地震参数赵博等(2013)，其他地震参数引自中国地震台网中心地震目录.

表 1 地震参数 Table 1 The parameters of the earthquakes

表 2 Table 2

表 2 速度模型 Table 2 Velocity model 层厚度/km α/(km·s-1) β/(km·s-1) 3 4.88 2.86 5 5.80 3.40 14 6.04 3.55 21 6.82 3.98 上地幔顶部 7.61 4.45

表 2 速度模型 Table 2 Velocity model

以2014年11月22日四川康定M6.4地震为例，简单介绍震源机制反演过程.首先，计算理论波形并与观测波形进行波形拟合(图 3)，同时根据最佳拟合结果进行震源机制随深度变化的拟合度分布分析(图 4).最终反演得到本次地震震源机制解为：节面一：走向238°，倾角80°，滑动角-165°；节面二：走向145°，倾角75°，滑动角-10°；矩震级M_W=5.81，标量地震矩：M₀=6.53×10¹⁹N·cm.对于本次地震，前人也给出了反演结果：节面一：走向234°，倾角81°，滑动角-172°；节面二：走向143°，倾角82°，滑动角-9°(易桂喜等，2015)；康定地震发生在鲜水河断裂，震源机制解显示，本次地震为高倾角走滑型，结合震源区断层走向，推断该次地震为左旋走滑型地震.除本次康定M6.4地震外，鲜水河断裂上的大部分M≥4.0地震均为左旋走滑型.鲜水河断裂是川滇块体和巴颜喀拉块体的交界区，川滇菱形块体的顺时针旋转及鲜水河断裂的高倾角断层产状结构，是该地区左旋走滑型地震的主要成因.其他地震震源机制反演结果见图 5，表 2给出了本研究结果与GCMT(Global Centroid Moment Tensor)结果对比.

图 3

Fig. 3

图 3 2014年11月22日四川康定M6.4地震波形拟合 红色为观测波形，蓝色为理论波形.Z、R、T分别表示垂直、径向和切向分量，右侧字母是台站代码.波形上方数字为互相关系数，波形下方数字为时间偏移量(单位：s). Fig. 3 The synthetic waveforms of Kangding M6.4 earthquake on 22 Nov. 2014 The red waveforms are observation waveforms and the blue waveforms are synthetic waveforms. Z, R, T indicate vertical, radial, and tangential cannel. The characters are stations code. The numbers above the waves are the cross-correlation coefficients. The number under the wave are the time shift.

图 4

Fig. 4

图 4 拟合度-深度分布图 箭头位置为最佳拟合结果. Fig. 4 Fitness-depth distributions The location of the arrow is the best fitting result.

图 5

Fig. 5

图 5 震源机制解与应力场反演结果 黑色箭头代表主压应力场方向，红色沙滩球为M≥5地震震源机制，蓝色沙滩球为M＜5.0地震震源机制. Fig. 5 The results of focal mechanism and stress filed inversions The black arrows indicate the directions of the principal compressive stress. The red beach balls are the focal mechanism of earthquake with magnitude M≥5.0, and the blue beach balls are the focal mechanism of earthquake with magnitude M < 5.0.

通过大量震源机制解资料，可以得到区域的应力场分布情况(许忠淮等，1984)，如果对一个地区许多地震的P轴方向和T轴方向做统计平均，就有可能获得该地区构造应力的方向(陈运泰和顾浩鼎，2007).本研究采用力轴张量计算法(钟继茂和程万正，2006)对四川地区构造应力方向进行了反演.

图 5给出了汶川M8.0主震和四川地区2010年1月1日至2017年12月31日发生的M≥4.0中强地震震源机制解的空间分布.在龙门山地区，震源机制呈现明显的空间分布差异，这与前人的研究一致(易桂喜等，2012；郑勇等，2009).龙门山断裂带中部，震源机制以逆冲型为主，与青藏高原向东推覆挤压运动有关.断裂带西北段青川—平武地区，震源机制以逆冲为主，兼有走滑分量.青川—平武段的剪切波分裂研究显示，该地区快剪切波偏振方向呈现与断裂带走向一致的特点(石玉涛等，2009)，这与该地区震源机制解兼有走滑分量相符.在理县地区，震源机制解以走滑型为主.汶川地震余震序列重新定位结果显示，该地区余震与龙门山断裂带走向共轭分布(即呈NW向)，即沿米亚罗断裂走向(Zhao et al., 2011; 易桂喜等，2012).米亚罗断裂为高倾角断层，这与左旋走滑型震源机制解反映结果一致.利用龙门山断裂带地震震源机制解(除芦山序列)，反演得到该区域主压应力场方向为105.9°.前人剪切波分裂研究结果显示，该区域快剪切波偏振方向也为近WNW向(石玉涛等，2009).芦山地震发生在龙门山前山断裂南端(易桂喜等，2016；赵博等，2013)，为逆冲型地震.利用芦山地震主震及其余震震源机制解反演得到了本地区的主压应力场方向为112.0°.剪切波分裂结果显示，MDS(雅安台)快剪切波偏振方向为106°(石玉涛等，2009).剪切波分裂快波偏振方向与区域主压应力场一致，通过震源机制解应力场反演和剪切波分裂分析结果推断，龙门山地区不同位置应力场方向有差异，但总体方向为WNW—近EW向.

九寨沟M7.0地震也发生在巴颜喀拉地块，但其发震断层不在龙门山断裂带区域内，其震源机制显示为左旋走滑型.九寨沟地震发生在虎牙断裂、岷江断裂、塔藏断裂和雪山梁子断裂四条断裂所围限区域的树正断裂(易桂喜等，2017b)，震中距龙门山后山断裂的平武—青川段的直线距离约为100 km.本研究反演得了主震震源机制解，节面一：走向330°，倾角70°，滑动角-15°；节面二：走向65°，倾角76°，滑动角-159°.根据区域断层走向推断，节面一为发震断层面，本次地震为左旋走滑型震源机制.利用主震及余震序列震源机制解反演的到该地区构造主压应力场方向为106.7°，与龙门山地区压应力轴方向相近.

川滇地块与巴颜喀拉地块交汇区的鲜水河断裂也是四川地区主要的中强地震发断层.鲜水河断裂是分割川滇块体和巴颜喀拉块体的一条呈弧形的大型活动断裂带(Allen et al., 1991).该断裂带北西起与甘孜北，向南东经炉霍、道孚、乾宁、康定等地(唐荣昌和韩渭滨, 1993; 易桂喜等, 2015).康定M6.4和M5.9地震、炉霍M5.4地震、道孚M5.0地震均发生在鲜水河断裂上，震源机制都为左旋走滑型，这与川滇菱形地块顺时针旋转有密切关系.剪切波分裂结果显示，DFU(道孚台)快剪切波偏振方向为108°，SMI(石棉台)快剪切波偏振方向为95°(石玉涛等，2013).鲜水河断裂带自晚更新世以来，构造应力场为近东西向挤压和近南北向拉张，在此应力作用下，鲜水河断裂开始大规模水平错动，并以左旋走滑为主(谢富仁等，1995).2010—2017年，鲜水河断裂上共发生了多次中强地震(M≥4.0)，由北向南分别有炉霍、道孚、康定、泸定、九龙等地震，这些地震均为左旋走滑型，利用震源机制解反演的主压应力场方向为108.5°，近EW向.

越西M5.2地震发生在四川东南部，是川滇块体的东部边界.该地区除了越西M5.2地震，还发生多次M≥4.0地震，例如昭觉M4.6，木里M4.2、M4.4、M4.6，普格M4.0等.这些地震震源机制均为走滑型，反演得到该区域主压应力轴方向为147.5°，近NW向.该地区快剪切波偏振方向也呈现为近NW向：MLI(木里台)为161.4°，PGE(普格台)为144.3°，YYC(园艺场台)为126.6°(石玉涛等，2013).剪切波分裂结果和本研究结果均显示该地区的主压应力轴方向为近NW向.

川滇块体西部存在着金沙江断裂、德格—乡城断裂、理塘—义敦断裂、德钦—中甸断裂等一系列地震断层.GPS资料显示，巴塘、理塘一带呈明显的张性运动(李煜航等, 2014)，本研究显示，该地区发生的理塘M5.2和得荣M5.2地震以及多次5.0≥M≥4.0地震，其震源机制解都为正断层.川滇地块西部位于扩展裂隙区，既有红河断裂尾端的拉伸区(虢顺民等，1996)又有横切川滇菱形块体的丽江—小金河断裂在区域内交汇，青藏高原东扩及中下地壳物质流动的拖曳作用与构造部位的特殊性叠加形成了显著的拉张应力(王晓山等，2015).剪切波分裂快波偏振方向显示：LTA(理塘台)为48.3°，LGH(泸湖台)为5.5°(石玉涛等，2013).本研究结果反演的到的主压应力轴方向为12.2°，其仰角近40°，该区域的压应力轴方向为NEN向.

震源深度的准确定位是地震学研究的一个基本但重要的内容，准确的震源深度对理解发震断层结构、板块运动和构造，识别非天然地震等有重要意义(Waldhauser and Ellsworth, 2000; Stein and Wiens, 1986; Heyburn et al., 2013; 高原等, 1997).本研究采用面波振幅谱能量衰减特征定位方法，对研究区域内中强地震进行深度定位，该方法能够准确测定震源深度.Rayleigh波在0.01~0.08 Hz频段内振幅谱特征与震源深度有密切关系，其频谱形状随深度变化比较敏感.在某一频段内，Rayleigh波振幅谱能量会出现陷波特征，即能量突然衰减.能量衰减对应的频段与震源深度有密切关系，可以根据这种关系进行深度定位.

Aki和Richards给出了Rayleigh波辐射花样(Aki and Richards, 2002)，公式(1)：

(1)

其中，ϕ为震中-台站方位角，ω_n为角频率，k_n为波数，M为地震矩张量，α, β分别为P波和S波速度，ρ为密度，r₁为水平位移向量特征函数，r₃、r₄为剪切应力和正应力特征函数，h为震源深度.

Rayleigh波的振幅公式显示，振幅A(ϕ, ω_n)是观测方位ϕ与频率ω_n的函数，受震源机制和震源深度的影响.矩张量M可以通过其他反演方法获得，或者利用GCMT已知等矩张量目录数据.P波和S波速度，本研究选用全球大尺度速度结构，例如，AK135.r₁、r₃、r₄可以通过解运动应力向量的微分方程获得，该微分方程包含密度、波数以及拉梅常数(λ和μ).

图 6为理论逆冲型地震(倾角为45°、滑动角为90°；当滑动就为-90°时为正断层)和走滑型地震(倾角90°、滑动角0°)在震源深度为15 km和5 km时，Rayleigh波振幅谱在360°方位角上分布特征.图中虚线范围内为面波能量锐减区，该区域出现的频率范围与震源深度有关.

图 6

Fig. 6

图 6 理论面波振幅谱 (a)和(b)中震源深度为15 km和5 km，震源机制解走向、倾角、滑动角为200°，45°和90°的理论面波振幅谱；图 6(c, d)中震源深度为15 km和5 km，震源机制解走向、倾角、滑动角为200°，90°和0°的理论面波振幅谱(A为振幅值). Fig. 6 The amplitude spectrum of the synthetic surface wave (a) and (b) are the synthetic suface wave amplitude spectrum with the depths are 15 km and 5 km which focal mechanisms are 200°, 45°, 90°; (c, d) (c) and (d) are the synthetic surface wave amplitude spectrum with the depths are 15 km and 5 km which focal mechanisms are 200°, 90°, 0°.

以芦山M7.0地震为例，该地震震源机制解为：走向200°、倾角45°、滑动角95°(赵博等, 2013).为了选择振幅能量衰减明显的观测台站，利用震源机制解计算5 km、10 km、20 km深度全方位角理论振幅谱(图 7).图 7显示不同深度振幅谱特征有明显变化，在震中-台站方位角为40~170°和220~350°范围内，振幅能量衰减特征较为明显.因此选取该方位角范围内的国家地震台网台站用于深度定位.

图 7

Fig. 7

图 7 芦山地震面波振幅谱 (a), (b), (c)分别为以芦山地震震源机制解计算5 km，10 km和20 km时理论振幅谱(A为振幅值). Fig. 7 The amplitude spectrum of Lushan earthquake (a), (b), (c) are the synthetic amplitude spectrum which are calculated by the Lushan earthquake focal mechanism.

以QZN台(琼中台)的观测波形为例.计算QZN台的理论地震图振幅谱并与观测波形振幅谱进行对比，在1~50 km深度搜索振幅谱能量衰减区位置的最佳拟合深度.图 8a显示了7~16 km深度理论波形振幅谱与观测波形振幅谱对比，在11.4 km深度时，频谱特征最为相似，能量衰减频率出现的位置最接近.图 8b为国家地震台网AXX、FZCM、GTA和JYG台的测量结果.M7.0芦山地震震源深度定位选用了48个国家地震台网台站，最终确定本次地震的震源深度为13.1 km，均方根误差为1.8 km.表 3为研究区域2010年1月1日至2017年12月31日M≥5.0地震震源深度定位结果和其他机构深度定位结果.

图 8

Fig. 8

图 8 芦山地震震源深度单台定位结果 (a) QZN台结果; (b) AXX, FZCN, GTA, JYG台定位结果；红色曲线为观测波形振幅谱，蓝色曲线为最佳拟合深度振幅谱.曲线右边的数字为深度. Fig. 8 The results of single station location of Lushan earthquake (a) The location result of QZN station; (b) The location results of AXX, FZCN, GTA, JYG stations. The red curves indicate the amplitude spectrum of observation waveforms. The blue curves indicate the amplitude spectrum of the synthetic waveforms. The numbers on the right are the depth.

表 3 Table 3

表 3 本研究震源机制反演结果与GCMT对比 Table 3 The results of focal mechanism inversions compared with GCMT 序号 本研究结果 GCMT结果 矩震级 MW 面波震级 MS 走向 倾角 滑动角(°) 走向 倾角 滑动角(°) 0 227 40 116 231 35 138 7.9 8.0 15 55 70 357 68 63 1 237 81 160 236 79 176 4.8 5.0 330 70 10 327 86 11 2 236 47 105 233 53 116 4.9 5.0 35 45 75 14 44 59 3 43 81 -160 45 83 -165 5.2 5.4 310 70 -10 313 75 -7 4 70 50 90 79 50 108 4.9 5.2 250 40 90 232 44 70 5 206 85 -170 208 86 179 5.3 5.5 115 80 -5 298 89 4 6 200 45 95 212 42 100 6.5 7.0 15 45 85 19 49 81 7 233 50 112 215 45 100 5.1 5.4 20 45 65 21 46 80 8 175 50 60 177 43 74 4.7 5.4 37 48 121 17 50 103 9 215 50 70 221 45 94 5.1 5.4 65 44 112 35 45 86 10 290 60 -75 294 51 -55 4.9 5.2 82 33 -114 66 50 125 11 355 75 35 349 80 25 5.0 5.2 255 56 162 254 66 169 12 238 80 -165 233 89 -175 6.0 6.4 145 75 -10 143 85 -1 13 242 80 -170 239 84 -175 5.7 5.9 150 80 -10 148 85 -6 14 335 50 75 340 46 91 4.7 5.0 178 42 107 158 44 89 15 39 64 -134 56 53 -122 4.9 5.2 285 50 -35 281 47 -55 16 65 76 -159 243 82 -168 6.5 7.0 330 70 -15 151 79 -8 17 20 70 45 21 54 42 5.0 5.4 271 48 153 263 57 135

表 3 本研究震源机制反演结果与GCMT对比 Table 3 The results of focal mechanism inversions compared with GCMT

汶川地震后，本地区发生了两次M7.0地震，即芦山地震和九寨沟地震.本研究给出了芦山M7.0震源深度为13.1 km，与USGS给出的震源深度14 km较为一致，CENC地震目录给出的深度为17 km，三者的差别不到4 km.对于九寨沟M7.0地震的震源深度，本研究给出的深度为8.2 km，USGS为9 km，CENC的定位深度为10 km，赵博等(2018b)利用干涉地震震源成像法定位的到本次地震震源深度为9 km，几种结果基本一致.这两次强震都发生在巴颜喀拉地块与华南地块交界附近，综合本研究结果和国内外其他研究机构的定位结果，在龙门山断裂带及其西部区域，地震震源深度的优势分布应该在10~20 km之间的中上地壳.这与龙门山断裂带的叠瓦状结构有关，其基地面大约在20 km以上倾斜展布.

康定M6.4地震和康定M5.9余震的震源深度分别确定为13.2 km和13.1 km，USGS结果皆为9 km，CENC结果分别为20 km和16 km；而ISC地震目录给出的结果较浅，分别为2 km和6.3 km.位于同一断层西北方向约50 km的道孚M5.0地震震源深度为11.5 km，USGS结果为10 km，CENC结果为8 km，ISC结果为12.5 km.鲜水河断裂西北端的炉霍M5.4地震，定位的震源深度为9.5 km，CENC结果为10 km.综合考虑整个鲜水河断裂地震，本研究认为，鲜水河断裂带上的地震震源深度优势分布应该为10 km左右(表 4).

表 4 Table 4

表 4 本研究深度定位结果和其他研究结果对比 Table 4 The depth location results compared with others results 序号 本研究结果 km 均方根误差 CENC结果 km USGS结果 km ISC结果 km 1 11.5 1.7 8 10 12.5 2 13.0 2.0 20 10 24.6 3 9.5 2.0 10 42.7 18.2 4 11.7 1.7 6 39.5 16.3 5 9.4 1.2 15 8.0 14.9 6 13.1 1.8 17 14 18 7 8.5 0.8 15 10 15.2 8 21.9 2.6 27 10 26.8 9 7.6 1.2 17 15.7 16 10 7.6 1.3 9 10 10 11 10.4 2.4 10 9.3 17.6 12 13.2 2.1 20 9 2 13 13.1 2.3 16 9 6.3 14 12.0 1.9 20 10 17.8 15 8.9 1.6 19 22 - 16 8.2 1.1 10 9 - 17 9.1 1.2 10 10 -

表 4 本研究深度定位结果和其他研究结果对比 Table 4 The depth location results compared with others results

位于川滇块体西南的理塘断裂、巴塘断裂、金沙江断裂、德格—乡城断裂等的地震：白玉M5.5、得荣M5.2、理塘M5.2，震源深度分别为9.4 km、7.6 km、8.9 km.川滇块体西南的这几次地震的震源深度均小于10 km.川滇块体西部，由于下地壳物质东移，对上地壳的拖曳作用，使该地区上地壳呈张性应力状态，而震源机制解反演与应力场反演也证明了这一点.几次中强地震震源深度小于10 km，也表明该地区上地壳受张性拉张应力，其中强地震发震孕震层位于10 km以上的上地壳.

青藏高原东缘的四川地区是中国大陆的地震多发区，自2010—2017年，该地区共发生18次M5.0以上地震和100多次5.0＞M≥4.0地震，其中两次M7.0和一次M6.4强地震.这些地震主要分布在四川西部和南部的巴颜喀拉地块和川滇地块内部以及与四川盆地(华南地块)的交汇区.青藏高原物质向东运移在龙门山地区受到华南地块的阻挡，产生挤压推覆作用.本研究震源机制解反演结果表明，龙门山地区的地震主要以逆冲型震源机制为主，兼有少量走滑分量.在理县附近，震源机制显示为走滑型，这与区域内NW向的高倾角米亚罗断裂有关.而川滇块体西南部的金沙江断裂、巴塘断裂、理塘断裂、德格—乡城断裂等地区，地壳运动显示为张性运动，该地区的震源机制表现为正断层性质.位于川滇块体和巴颜喀拉块体交界的鲜水河断裂，是两大块体的分界，该断层上发生的地震多为走滑型地震.

震源机制解反演的得到的龙门山地区构造应力场主压应力轴方向为105.9°，芦山地震震源区应力场反演结果为112°，九寨沟地震震源区主压应力轴为106.7°.因此，龙门山及其周边地区构造应力场方向应为WNE—近EW向.鲜水河断裂被认为是川滇块体和巴颜喀拉块体的分界.2010年以来，鲜水河断裂上共发生了4次地震：道孚M5.0地震、炉霍M5.4地震和两次康定(M6.4、M5.9)地震以及在九龙、泸定等发生的M4.0左右的地震.这些地震震源机制非常相似，均为走滑型，这与该地区近东西向挤压和近南北向拉张的构造应力有关.应力场反演得到的主压应力轴为108.5°WNW向.川滇块体东部边界地区，其主压应力轴反演结果为147.5°NW向.川滇地块的西部地区，地壳构造应力场呈现近EN向拉张特性，反演得到主压应力轴为12.2°.

中强地震深度定位显示，在龙门山断裂带及其西部区域，地震震源深度的优势分布应该在10~20 km之间的中上地壳.这与龙门山断裂带的叠瓦状结构有关，其基地面大约在20 km以上倾斜展布.鲜水河断裂带上的地震震源深度优势分布在为10 km左右.川滇块体西部，上地壳收到中下地壳物质东向流动的拖曳拉张，使得上地壳受到近EW向的拉张应力控制，其中强地震震源深度均在10 km以上.