王卫民1 何建坤1 郝金来2 姚振兴2
1大陆碰撞与高原隆升实验室,中国科学院青藏高原研究所,北京,100101
2中国科学院地质与地球物理研究所,北京,100029
北京时间2017年11月18日6点34分 (UTC 2017-11-17 22:34:20) 在我国西藏林芝地区(USGS定位结果为北纬29.8347°N,东经94.9092°E,震源深度10km)发生了Mw6.4级地震。地震发生后,我们从IRIS数据中心下载了地震数据资料用于研究地震震源机制和震源破裂过程。选取其中信噪比较高并且沿方位角分布比较均匀的20个远场P波波形(震中距位于 范围之内)和22个SH波波形资料进行点源模型的震源机制解反演;根据反演结果构建有限断层模型进行震源过程反演,获得了这次地震破裂滑动分布的初步结果。计算得到的地震矩为5.0×1018Nm,Mw=6.4。最大滑动约60cm。
初步结果表明:这次地震为震级Mw6.4,震源深度约10km的高倾角逆冲型地震,破裂持续时间约13秒,破裂没有明显的方向性,断层面上的滑动分布比较集中。计算震中区的理论最大地震烈度(中国地震烈度表,2008)约为VIII-IX度。
图 1 西藏林芝M6.5级地震震源机制解
采用下半球投影,同时给出了点源模型的P波垂向位移理论图(红线)与资料(黑线)的拟合情况。图形下方给出了两组节面解(左下,λ,δ,θ,h分别表示错动倾伏角、断层倾角、断层走向、震源深度)和点源模型的震源时间函数(右下)。
Figure 1 Focal mechanism of the M6.5 earthquake, Nyingchi, Tibet, China.
Lower-hemisphere projection is used here. The observed P wave records (black line) and the synthetic seismogram (red line) based on the simple point model are compared. The parameters of two possible fault planes are listed also, with l, d, q ,h indicating the rake angle, dip angle, strike direction and source depth respectively. The obtained source time function are plotted.
图2 有限断层模型的远场P和SH波形拟合
左列为P波垂向位移理论图(红线)与记录(黑线)的对比;右列为SH波切向位移理论图(红线)与记录(黑线)的对比;并给出每个记录的方位角(左上)、震中距(左下)、振幅比(上)和台站名(右)。
Figure 2 Comparison of the observed teleseismic P and SH records with the synthetic waveforms.
Comparison of observed (black) and synthetic waveforms (red) for P vertical components, and SH components are shown in the panel, the azimuth, epicentral distance, amplitude ratio and station name are also indicated
图 3 地震断层面上滑动分布
图中分别给出了有限断层模型的地表投影(左上)和三维示意图(右上)
Figure 3 Inverted slip distribution on the fault.
Surface projection of fault model ( up-left ) and 3D view of the fault model ( up-right ) are shown
图 4 地震断层滑动分布投影
Fig. 4 The map shows the projection of slip distribution on the ground
图 5 根据破裂过程模型计算的理论烈度分布
Fig. 5 Predicted Shakemap calculated from slip model