据中国地震台网中心测定，北京时间2015年4月25日14时11分26秒，尼泊尔发生M_S8.1大地震，震中位置28.2°N、84.7°E，震源深度20 km.截至5月6日，尼泊尔境内已经造成至少7566人丧生，超过14500人受伤(中国地震局，2015).强震还引发珠穆朗玛峰雪崩，部分基地大本营被掩埋.此次地震震中距离尼泊尔首都加德满都约80 km，其余震由近西北向近东南方向迁移，震源深度多集中在20 km以上(图 1)，属于浅源型地震，因此造成地表建筑破坏严重.

图 1

Fig. 1

图 1 本研究震源机制反演所用台站(蓝色三角)及尼泊尔MS8.1地震序列分布.地震序列数据来自USGS. 星号代表大于7级地震震中，小圆圈代表7级以下余震震中，填充的颜色代表震源深度，色标位于图底 Fig. 1 Distribution of the seismic stations (blue triangles) used in focal mechanism inversion in the present study and the Nepal MS8.1 earthquake sequences. The earthquake sequence data are from USGS. Stars denote the epicenters of the earthquakes with magnitude greater than 7.0, whereas small circles denote the epicenters of events with magnitude smaller than 7.0, the color filled in which represents the focal depth. The scale for focal depth is shown at the bottom

此次尼泊尔M_S8.1地震引发多次中强余震，精确的震源参数，如发震断层面、震源质心深度、矩震级等是地震学研究中重要的基础资料，因此快速求取可靠的震源机制解对后续研究非常重要.对于M_S≥7.0强震，其破裂过程具有一定的尺度，点源模型不适应于描述震源破裂的真实情况，故通常采用非点源模型，比如有限断层破裂模型(Chen et al.，2005; Zhang et al.，2009; 王卫民等，2013; 张勇等，2015).而对于中等地震，其破裂尺度小，将震源简化为点源可以很好地反映震源破裂特征(郑勇等，2009；陈伟文等，2012；Chen et al.，2015).本研究就针对尼泊尔M_S8.1主震后发生的中等地震，求取震源机制解，探讨震源破裂的变化特征，讨论尼泊尔M_S8.1地震对青藏高原及中国西南地区的影响作用.

本研究采用国家数据备份中心提供的西藏台网波形资料(郑秀芬等，2009)，选择4月25日—5月16日尼泊尔M_S8.1地震5次余震(5.0≤M_S≤6.5)及西藏定日M_S5.9地震，总共6次中等地震事件.通过波形筛选，挑选震中距6°以内共11个台站的三分量记录(图 1).速度模型采用基于人工源地震剖面结果给出的平均速度结构分层模型(滕吉文等，2015).

震源机制解采用gCAP(generalized Cut and Paste)方法求取(Zhu and Ben-Zion，2013).该方法是在CAP法(Zhu and Helmberger，1996)的基础上增加各向同性源(ISO)和补偿线性偶极子(CLVD)两个分量，求解地震全矩张量解(李圣强等，2013; 张广伟和雷建设，2015).本研究将非双力偶分量ISO和CLVD约束为0，求解地震的最佳双力偶节面解.与CAP方法相同，gCAP法在反演过程中将三分量全波形分为Pnl和面波两部分，对Pnl和面波两部分赋予不同权重，采用不同频段滤波后参与反演，通过计算理论和实际波形的拟合误差函数，利用网格搜索得到最小误差的最优解.该方法允许每个时间窗理论波形和实际观测波形相对时移拟合，大大减少了速度模型不精确及地壳速度横向不均匀性的影响(Zhu and Helmberger，1996；郑勇等，2009)，并且前人研究表明，仅使用二个台站记录的区域波形资料也可以得到较好的震源机制结果(Tan et al.，2006).本研究反演过程中，Pnl和面波滤波范围分别为0.05~0.2 Hz和0.02~0.1 Hz，走向、倾角和滑动角的搜索间隔均为5°，深度为2 km.格林函数采用频率-波数法(FK)计算(Zhu and Rivera，2002)，采样间隔设为0.1 s，采样点2048个.

表 1给出6次地震事件最佳双力偶节面解及震源质心深度，图 2详细展示事件1、2、3和6理论与实际波形拟合图以及震源机制解(下半球投影)，整体上波形拟合效果较好.每个事件至少有7个台站参与反演(表 1)，我们依据波形平均拟合相关系数对震源机制结果质量进行分类，如果平均拟合相关系数＞65%则为A，在60%~65%之间为B，小于60%则为C.由表 1可以看出，6次地震事件震源机制解的质量均在B及以上，说明我们所获得的震源参数较为可靠.

表 1 (Table 1)

表 1 中等地震震源机制解 Table 1 Focal mechanism solutions of the moderate-sized earthquakes 代号 发震时刻(UTC)年-月-日 时∶分∶秒 纬度/(°N) 经度/(°E) 质心深度/km MW 节面I 节面II 使用台站个数 震源类型 质量分类 走向/(°) 倾角/(°) 滑动角/(°) 走向/(°) 倾角/(°) 滑动角/(°) 1 2015-04-25 09∶17∶05.0 28.40 87.30 6 5.7 172 52 -98 6 39 -79 7 正断 A 2 2015-04-25 17∶42∶53.0 28.20 85.90 16 5.0 138 81 -111 25 22 -25 9 正断 B 3 2015-04-25 23∶16∶14.0 27.80 85.00 10 4.8 312 69 -141 206 54 -26 7 走滑 B 4 2015-04-26 16∶26∶05.0 27.76 85.75 10 4.9 145 85 133 241 43 8 9 走滑 A 5 2015-05-12 07∶36∶53.0 27.50 86.10 8 5.9 126 70 104 269 24 56 9 逆冲 B 6 2015-05-16 11∶34∶10.0 27.60 86.20 6 5.3 132 73 106 268 23 48 9 逆冲 A

表 1 中等地震震源机制解 Table 1 Focal mechanism solutions of the moderate-sized earthquakes

图 2

Fig. 2

图 2 表1中事件1、2、3和6地震震源机制解及理论(红色)和实际(黑色或蓝色)波形对比图，拟合较好观测波形用黑色表示，拟合较差用蓝色表示；波形下方数字表示理论相对实际波形的相对移动时间和二者的相关系数；左侧大写字母表示台站名，台站下方数字为震中距(km)和相对偏移时间(s) Fig. 2 The focal mechanisms and the comparison between synthetic (red) and observed (black for better fits or green for bad fits) waveforms of the events 1、2、3 and 6 in Table 1. The numbers below each traces are relative time shifts and cross-correlation coefficients. The station names are given on the left and the numbers below each station are epicentral distances (km) and relative time shifts (s)

事件1西藏定日地震表现为正断型，质心深度6 km，与韩立波反演得到结果(中国地震局地球物理研究所，2015)较为一致.事件2也为正断型地震，但与事件1有所区别，其不仅兼有少量走滑分量，且质心深度也较深，为16 km.事件1和2位于发震断层的上盘区域(图 3)，表明该区域受到主震同震位移的影响，表现出应力拉张的变化特征.事件3和4在M_S8.1主震后2天内发生，均表现为走滑型地震(表 1)，震中位于主震的东南向，其他余震也多位于主震的东南方向，说明主震破裂可能是东南向为主的单侧破裂.事件5和6发生于5月12日M_S7.5强余震之后，表现为逆冲类型(图 3)，与M_S7.5强余震的逆冲破裂机制相同；另外，事件5和6在全球矩张量解中也给出结果(www.globalcmt.com，简写为GCMT)(图 3)，我们将这两次事件与GCMT结果进行对比，均为逆冲性质，也证实了本研究结果的可靠性.

图 3

Fig. 3

图 3 本研究获得的中等地震震源机制解(黑色震源球)和事件5、6的全球矩张量解(红色震源球，由GCMT下载).星号代表大于7级地震，小圆圈代表7级以下余震，填充的颜色代表震源深度，色标位于图底 Fig. 3 Focal mechanism solutions of the moderate-sized earthquakes obtained by present study (black focal sphere) and downloaded from global central moment tensor (red focal sphere from GCMT).Stars denote the earthquakes with magnitude greater than 7.0, whereas small circles denote the events with magnitude smaller than 7.0, the color filled in which represents the focal depth. The scale for focal depth is shown at the bottom

尼泊尔位于印度和欧亚板块交界处，北部为喜马拉雅山，地形起伏明显，在其境内分布有多条逆冲断裂，分别有喜马拉雅主前锋断裂(MFT)、主边冲断裂(MBT)和主中冲断裂(MCT)(许志琴等，2013).尼泊尔近代地震活动频繁，比如1934年1月15日发生的M_S8.4尼泊尔—比哈尔邦地震，造成约10600人死亡，而近期发生的地震主要集中在尼泊尔东部和西部，中部地区不太活跃，这次M_S8.1地震就发生在该地区，说明此次地震是印度板块北向俯冲挤压欧亚板块(Huang and Zhao，2006; Li et al.，2008; Lei et al.，2014)应力长期积累的结果.

我们获取的6个中等地震震源机制解，分别有2个正断类型、2个走滑类型和2个逆冲类型(表 1).其中，2个正断类型地震(事件1和2)位于主震的东北方向，即发震断层的上盘(图 3)，这与2011年日本9级大地震上盘主要分布正断型地震非常类似(Asano et al.，2011)，主要受到主震同震位移的影响，上盘表现出应力拉张的特征.但事件1与事件2有所不同，事件1与尼泊尔M_S8.1主震相距较远，约250 km，位于中国境内的西藏定日县，震源机制解显示发震断层面近南北向(表 1)；静态应力触发研究表明，尼泊尔地震产生的应力加载在中国大陆主要集中在本次地震邻近的西藏及新疆地区部分断层上(盛书中等，2015)；同震位移结果也显示，主震对雅鲁藏布江地区和拉萨块体南北走向的正断层影响较大(张贝等，2015)；同时，从余震序列分布特征可以看出(图 3)，余震主要沿东南向扩展，与事件1所处的东北方向不在一个方位.因此本研究认为事件1不是尼泊尔地震的余震，而是由于主震发生后引发邻近区域构造应力调整所造成的触发型地震.事件2虽然也为正断类型，但其包含一定的走滑分量，且其节面I走向(138°)与余震序列的走向也较为一致(图 3和表 1)，因此我们认为该地震为尼泊尔地震的余震，其正断型破裂机制是由于主震同震位移影响导致发震断层上盘应力拉张而引起.

事件3和4为走滑型地震，这2个地震位于主震破裂延展的东南方向上，震中处在主震与5月12日M_S7.5强余震之间(图 3)，可能表明随着破裂往东南方向延伸，断层的走滑特征愈加明显.这样的余震分布特征与2008年汶川M_S8.0地震余震序列相类似(张勇等，2009)，在主震处为逆冲型，沿破裂方向分布的余震其走滑分量逐渐增强.但与汶川M_S8.0地震扩展约300 km不同(Zhang et al.，2009)，此次尼泊尔M_S8.1地震破裂长度相对较短，大约150 km(张勇等，2015)，我们推测主震的东南向破裂并未扩展完全，可能受到了阻挡.5月12日在主震破裂的东南向边缘又发生了一个M_S7.5强余震，我们反演的事件5和6也位于该区域，与M_S7.5强余震破裂机制相同，均表现为逆冲类型地震(图 3).上述不同时间、不同空间分布的中等余震震源破裂机制表明，在尼泊尔M_S8.1主震发生后，受到主震同震位移的影响，不同区域处于不同的应力状态，从而使中等余震表现出不同的震源类型.

此次尼泊尔M_S8.1地震位于青藏高原的南边界，主震的破裂过程(张勇等，2015)以及余震序列分布特征(图 1)显示，破裂沿东南方向迁移，表明此次M_S8.1地震对青藏高原及我国西南地区具有应力加载作用，建议应加强这些地区的地震监测预测工作.自2000年以来，在青藏高原及周边地区发生了一系列强震活动，如2001年昆仑山M_S8.1地震、2010年玉树M_S7.0地震、2008年汶川M_S8.0地震和2013年芦山M_S7.0、于田M_S7.3地震(张广伟等，2014)，说明印度板块与欧亚板块的持续碰撞对青藏高原及周边构造运动和地震活动的强烈影响作用(邓起东等，2010；徐锡伟等，2008，2011；闻学泽等，2011；Jia et al.，2012).此次尼泊尔地震，更印证了印度板块与欧亚板块的持续碰撞作用，说明青藏高原仍处于强震活跃期(邓起东等，2014).

致谢 感谢中国地震局地球物理研究所国家数字测震台网数据备份中心(doi:10.7914/SN/CB)为本研究提供地震波形数据.