中国地震  2020, Vol. 36 Issue (2): 324-332
2017年四川九寨沟MS7.0地震及余震定位研究
徐旭1, 徐锦承2, 张伟2     
1. 中国科学技术大学, 地球和空间科学学院, 合肥 230026;
2. 南方科技大学, 地球与空间科学系(深圳市深远海油气勘探技术重点实验室), 广东深圳 518055
摘要:采用特定震源台站校正定位方法对九寨沟MS7.0地震及其余震进行了精定位研究,结果显示九寨沟地震主震深度为19km左右,余震主要分布于5~15km深度范围;与单事件定位方法相比,其余震展布和断层的线性特征更吻合,更具有平面特征;与双差定位法相比,精定位事件覆盖范围更广。精定位结果显示,断层整体为虎牙断裂西北向的延伸,根据余震分布和特定震源区的校正项分布,认为此发震断层应分南、北两段;余震区长度为40km左右,西北段宽度为10km左右,东南段宽度为6.5km左右。
关键词特定震源台站校正    九寨沟地震    余震分布    发震断层    
Relocation of Main Shock and Aftershocks of the MS7.0 Sichuan Jiuzhaigou Earthquake
Xu Xu1, Xu Jincheng2, Zhang Wei2     
1. School of Earth and Space Sciences, University of Science and Technology of China, Hefei 230026, China;
2. Department of Earth and Space Sciences, Southern University of Science and Technology(Shenzhen Key Laboratory of Deep Offshore Oil and Gas Exploration Technology), Shenzhen 518055, Guangdong, China
Abstract: Relocation of the Jiuzhaigou earthquake and its aftershocks by using the source-specific station terms (SSTT) method has been performed in this paper. The results show that the Jiuzhaigou main shock is located at depth of about 19km, and the aftershocks are mainly distributed at the depth range of 5~15km. Compared with the results of single-event location method, the distribution of the aftershocks in our study is more consistent with the linear character of faults strike, and more likely to have plane character. Compared with the results of double-difference method, the spatial coverage of events could be wider. The results of relocation indicate that the faulting plate is from the northwest extension of the Huya fault. Based on the distribution of aftershocks and correction terms in the study area, we believed that the seismogenic fault could be divided into two sections, north and south. The length of aftershock area is approximately 40km, the width of northwest section is approximately 10km, and the width of southeast section is approximately 6.5km.
Key words: SSST method     Jiuzhaigou earthquake     Aftershock sequence     Seismogenic fault    
0 引言

2017年8月8日21时19分,四川省九寨沟地区发生MS7.0地震,主震发生后,震中区域陆续发生大量余震,部分余震震级达MS4.0。此次地震成为继2013年芦山MS7.0地震后,巴颜喀拉块体东部地区的又一次7.0级地震。九寨沟地区地质结构复杂,又处于强震多发地带,故有必要对该区域的地震序列时空展布进行详细的分析与研究。

地震定位对震区活动构造、地球内部结构研究等问题有重要意义(田玥等,2002)。高精度的地震定位结果有利于分析地震分布与活动构造之间的关系,更完备的地震事件检测有助于断裂带三维几何结构和地震演化过程的精确描述。地震定位精度的主要影响因素有台站覆盖度、走时拾取精度和速度模型精度。

传统的绝对定位法是利用走时信息通过网格搜索、高斯牛顿法和共轭梯度法等反演震源信息,但受速度模型和走时拾取精度的影响较大。为减少速度模型对定位精度的影响,前人提出了一系列相对定位法来提高定位精度,目前主要的精定位方法有:双差定位法(Waldhauser et al,2000)、主事件定位法(Spencer et al,1980)、震源与速度结构联合定位法(Crosson,1976)、联合震源定位法(Dewey,1972)、特定震源台站校正定位法(Richards-Dinger et al,2000Lin et al,2005)和交叉双差法(Tian et al,2016)等。其中,双差定位法已被国内的学者广泛应用于不同地区(华卫等,2003张洪艳等,2013房立华等,2015王健等,2016刘建明等,2017曹颖等,2018)。双差定位法利用2个地震走时差的观测值与理论计算值的残差最小来提高地震的相对定位精度,但是,双差定位法仅对地震群进行定位,且限制条件较为苛刻,部分事件不满足精定位条件。如果为了提高精定位事件数量而降低双差定位法的限制条件,又会降低精定位效果。

Lin等(2005)研究显示,特定震源台站校正定位法(Source-specific Station Term,简称SSST定位法)与双差定位法的定位效果相当,精定位事件的覆盖范围广,校正量计算与定位计算过程相互独立,因而可以使用任何单事件的定位算法。SSST定位法可以处理大量的数据,采用并行计算(Lin et al,20052007Nooshiri et al,2017)。然而,震相到时的精度会显著影响精定位的结果。

目前,人工拾取的震相到时精度为0.1~0.5s,由于人工拾取产生的定位误差可达数千米。采用波形互相关技术可以在一定程度上提高走时拾取的精度(Poupinet et al,1984),进而提高精定位的精度。使用CrazySeismic软件可以实现交互式互相关的走时拾取(Yu et al,2017),提高多震相拾取的精度。

本文首先采用特定震源台站校正定位法对九寨沟MS7.0地震及其余震序列进行精定位,其次对比了双差定位法,讨论并分析了九寨沟地震序列的时空展布规律,为进一步分析潜在的断层结构及断层形态提供基础支撑。

1 特定震源台站校正定位方法

基于走时的地震定位方法通常是通过最小化计算走时与理论走时的残差来实现的。走时残差主要受到速度结构误差和震相拾取随机误差的影响。三维速度结构误差是主要的影响因素,且其和震源区域的空间位置相关。考虑到相近事件的射线路径受到速度结构的影响相似,特定震源台站校正定位方法通过将校正项加到走时上来弥补速度模型带来的走时误差,可以部分解决速度模型误差带来的残差(Richards-Dinger et al,2000),从而提高定位精度。其走时残差的目标函数为

$ \varphi \left({{t_i}, {x_i}, {s_{ij}}} \right) = {\left\| {T_{ij}^{{\rm{obs}}} - {t_i} - {\tau _j}\left({{x_i}} \right) - {s_{ij}}} \right\|_p} $ (1)

式中,Tijobs为第i个事件在第j个台站的拾取到时;ti为第i个事件的发震时刻;τj(xi)为第i个事件在第j个台站的计算走时;sij为第i个事件在第j个台站的走时校正项;‖ · ‖p代表p的范数。

该目标函数可以采用两步迭代法求解,第一步利用静态台站校正方法获取初始的台站校正量,第二步利用事件和台站对时差更新校正量,计算流程如图 1所示,迭代过程如下:

图 1 SSST法计算流程

(1) 利用该地区的速度模型计算选定震相的理论走时表,进行单事件定位;

(2) 采用静态台站校正方法进行几次迭代,直至得到稳定的各个台站校正值,作为初始的台站校正量;

(3) 计算走时校正项,即SSST值;

(4) 校正震源台站对的到时(Tijobs =Tijobs-Sij);

(5) 利用校正到时进行精定位,直至收敛到预设的残差。

在计算过程中,第一次迭代需要人为设定起始的事件选取范围,记为Rmax,一般要求要覆盖所有余震的展布范围。另外,最后一次迭代的Rmax以及最少的临近事件个数N也需要测试确定来保证校正量空间变化的连续性。

2 资料选取与数据处理 2.1 构造背景

九寨沟MS7.0地震的主要断裂发生在巴颜喀拉块体东缘的活动构造带,处于塔藏断裂南东段,附近还有龙门山断裂带、岷江断裂、东昆仑断裂带、龙日坝断裂、虎牙断裂、白龙江断裂、阿万仓断裂、迭部-武都断裂、临潭-宕昌断裂等一系列断裂和断裂带(闻学泽,2018李陈侠等,2011高明星等,2015徐锡伟等,2017Kirby et al,2000Ren et al,2013Li et al,2011)。

图 2 研究区域构造情况及强震的历史分布
2.2 速度模型

九寨沟震源区海拔4km,附近的部分台站海拔也在1km以上。采用“全球AK135速度模型”进行定位,若不考虑地表地形的影响,则走时残差较大,定位误差也会较大。为了提高定位精度,本文采用的一维速度模型参考了周民都等(2006)嘉世旭等(2008)杨宜海等(2015)朱介寿等(2017)基于接收函数和人工测深剖面的结果,速度模型如图 3所示。

图 3 用于定位的一维速度模型
2.3 数据处理

自发震时刻至2017年9月1日0时,共筛选地震事件2521个,筛选规则为:选取ML≥3.5的地震事件450km内的台站,余震选取事件350km内的台站,其中18个台站震相清晰、噪声小,此外还有6个流动台站的走时资料,共计24个台站;方位覆盖角较为均匀,要求每个地震事件至少有5个走时记录;特定震源区选取最小半径为1km左右,满足最少的临近事件个数N=30;要求初始特定源区半径Rmax包含所有事件,取为55km。

2.4 精定位结果分析

根据单事件定位的初始位置,计算每个台站的一组P波台站校正量,取其平均值作为SSST计算的初始校正量。最终精定位的结果相较于单事件定位法,其事件分布更加线性,特别是沿着SE—NW方向展布更为清晰(图 4),与易桂喜等(2017)推测的树正断裂方向基本一致;与房立华等(2018)的双差精定位法相比,其定位结果相近,但是在震源深度上更集中于5~15km。余震近垂直分布特征明显,可以推测发震断层陡,整体余震分布沿着虎牙断裂的西北向展布。余震区长度为40km左右,西北段宽度为10km左右,东南段宽度为6.5km左右。

图 4 余震的分布和剖面图

进一步通过研究地震分布的规律,对比定位方法的可靠性。图 5为单事件定位法、SSST定位法和双差定位法定位的P波震源深度(房立华等,2018)的统计直方图,经过精定位后,SSST定位法和双差定位法的深度分布于5~15km,且SSST定位法的优势深度为7.5~10km,与梁建宏等(2018)的优势深度50%集中在10km附近的结论基本一致,说明SSST定位法对于震源深度的定位有一定优势。需要指出的是,相比双差定位法,SSST定位法的震源深度似乎整体更浅一些,主要原因可能是速度模型的差异,本文使用的速度模型近地表速度相对要大一些。

图 5 不同定位方法的深度分布直方图

图 6为选取的4个典型台站(JZGSC、SPASC、ZHQGS、WXTGS)P波的特定源区台站SSST校正项,相应的校正值绘制在该事件的位置上,由图可见,相对于真实的速度结构,一维速度模型存在强烈的速度异常结构。台站校正项为正值,说明该震源区到台站的射线路径上有低速异常;校正项为负值,说明震源区到台站的射线路径上有高速异常。为了更好地查看分布特征及特定源区相应台站的台站校正量,沿着AB线做剖面(A点:104.00°E、32.96°N,B点:103.70°E、33.40°N),选取测线左右1.5km内的事件,如图 7所示。由图 7的剖面可见,本次九寨沟MS7.0地震的余震主要集中于5~15km,多数余震分布于主震的西北方向,释放了大量能量,东南方向余震分布较少。从震源区校正项来看,余震具有明显的分区特征,有可能佐证了梁姗姗等(2018)给出的震源区域介质结构具有强烈的不均匀性的结论。

图 6 特定台站SSST校正项

图 7 AB线剖面的SSST校正量分布
3 结论

本文应用特定台站校正定位方法(SSST定位法)对九寨沟MS7.0地震及其余震进行了精定位。由于区域内地表地形复杂、速度变化较大,采用了基于接收函数和人工测深剖面的高分辨率速度模型进行定位,并与单事件定位法及双差定位法的结果进行了对比。通过对震后约3周的数据进行精定位,认为主震震源深度为19km左右,余震分布更集中,带状分布明显,发震断层较陡。西北段余震分布较多,深度浅;东南段余震分布较少,深度深。东南段与西北段的台站校正项差异较大,我们倾向于将发震断层以主震为中心分为南、北两段进行研究。尽管SSST定位法减少了邻近事件之间的相对误差,但绝对位置仍然存在偏差,而SSST定位法仅能反应震源区到台站射线路径上的速度异常。此外,三维速度模型的精度也会影响事件定位的绝对位置的精度。

致谢: 中国地震局参加流动台站布设的所有观测人员付出了辛勤劳动,中国地震台网中心提供了固定台站数据,中国地震局地球物理研究所房立华研究员提供了流动台站走时数据和相应的指导,在此一并表示衷心的感谢。
参考文献
曹颖、黄江培、付虹, 2018, 云南小湾水库蓄水后P波速度结构的双差地震层析成像研究, 中国地震, 34(4): 652-666.
房立华、吴建平、王未来等, 2015, 2014年新疆于田MS7.3地震序列重定位, 地球物理学报, 57(3): 802-808.
房立华、吴建平、苏金蓉等, 2018, 四川九寨沟MS7.0地震主震及其余震序列精定位, 科学通报, 63(7): 649-662.
高明星、陈桂华、徐锡伟, 2015, 地貌参数指示的临潭-宕昌断裂带最新构造隆升差异与地震活动, 地震地质, 37(3): 709-718.
华卫、刘杰、郑斯华等, 2006, 2003年云南大姚6.2、6.1级地震序列特征分析及地震触发研究, 中国地震, 22(1): 10-23.
嘉世旭、张先康, 2008, 青藏高原东北缘深地震测深震相研究与地壳细结构, 地球物理学报, 51(5): 1431-1443.
李陈侠、徐锡伟、闻学泽等, 2011, 东昆仑断裂带中东部地震破裂分段性与走滑运动分解作用, 中国科学:地球科学, 41(9): 1295-1310.
梁建宏、孙丽、刘杰, 2018, 2017年四川九寨沟MS7.0地震及余震精定位研究, 地球物理学报, 61(5): 2152-2162.
梁姗姗、雷建设、徐志国等, 2018, 2017年四川九寨沟MS7.0强震的余震重定位及主震震源机制反演, 地球物理学报, 61(5): 2163-2175.
刘建明、高荣、王琼等, 2017, 2017年8月9日精河6.6级地震序列重定位与发震构造初步研究, 中国地震, 33(4): 663-670.
田玥、陈晓非, 2002, 地震定位研究综述, 地球物理学进展, 17(1): 147-155.
王健、李春峰、雷建设等, 2016, 华北地区地震活动与地壳热结构关系研究, 地震学报, 38(4): 618-631.
闻学泽, 2018, 巴颜喀拉块体东边界千年破裂历史与2008年汶川、2013年芦山和2017年九寨沟地震, 地震学报, 40(3): 255-267.
徐锡伟、陈桂华、王启欣等, 2017, 九寨沟地震发震断层属性及青藏高原东南缘现今应变状态讨论, 地球物理学报, 60(10): 4018-4026.
杨宜海、梁春涛、苏金蓉, 2015, 用接收函数建立区域模型的震源机制反演及其在芦山地震序列研究中的应用, 地球物理学报, 58(10): 3583-3600.
易桂喜、龙锋、梁明剑等, 2017, 2017年8月8日九寨沟M7.0地震及余震震源机制解与发震构造分析, 地球物理学报, 60(10): 4083-4097.
张洪艳、张广伟、王晓山等, 2015, 2013年吉林前郭5.8级震群精定位及发震构造分析, 中国地震, 31(3): 518-528.
周民都、张元生、石雅鏐等, 2006, 青藏高原东北缘地壳三维速度结构, 地球物理学进展, 21(1): 127-134.
朱介寿、王绪本、杨宜海等, 2017, 青藏高原东缘的地壳流及动力过程, 地球物理学报, 60(6): 2038-2057.
Crosson R S, 1976, Crustal structure modeling of earthquake data:1.Simultaneous least squares estimation of hypocenter and velocity parameters, J Geophys Res, 81(17): 3036-3046.
Dewey J W, 1972, Seismicity and tectonics of western Venezuela, Bull Seismol Soc Am, 62(6): 1711-1751.
Kirby E, Whipple K X, Burchfiel B C, et al, 2000, Neotectonics of the Min Shan, China:Implications for mechanisms driving Quaternary deformation along the eastern margin of the Tibetan Plateau, Geol Soc Am Bull, 112(3): 375-393.
Li C X, Xu X W, Wen X Z, et al, 2011, Rupture segmentation and slip partitioning of the mid-eastern part of the Kunlun Fault, north Tibetan Plateau, Sci China:Earth Sci, 54(11): 1730-1745.
Lin G Q, Shearer P, 2005, Tests of relative earthquake location techniques using synthetic data, J Geophys Res:Solid Earth, 110(B4): B04304.
Lin G Q, Shearer P M, Hauksson E, 2007, Applying a three-dimensional velocity model, waveform cross correlation, and cluster analysis to locate southern California seismicity from 1981 to 2005, J Geophys Res:Solid Earth, 112(B12): B12309.
Nooshiri N, Saul J, Heimann S, et al, 2017, Revision of earthquake hypocentre locations in global bulletin data sets using source-specific station terms, Geophys J Int, 208(2): 589-602.
Poupinet G, Ellsworth W L, Frechet J, 1984, Monitoring velocity variations in the crust using earthquake doublets:An application to the Calaveras Fault, California, J Geophys Res:Solid Earth, 89(B7): 5719-5731.
Ren J J, Xu X W, Yeats R S, et al, 2013, Millennial slip rates of the Tazang fault, the eastern termination of Kunlun Fault:Implications for strain partitioning in eastern Tibet, Tectonophysics, 608(8): 1180-1200.
Richards-Dinger K B, Shearer P M, 2000, Earthquake locations in southern California obtained using source-specific station terms, J Geophys Res:Solid Earth, 105(B5): 10939-10960.
Spencer C, Gubbins D, 1980, Travel-time inversion for simultaneous earthquake location and velocity structure determination in laterally varying media, Geophys J R Astron Soc, 63(1): 95-116.
Tian X, Zhang W, Zhang J, 2016, Cross double-difference inversion for microseismic event location using data from a single monitoring well, Geophysics, 81(5): KS183-KS194.
Waldhauser F, Ellsworth W L, 2000, A double-difference earthquake location algorithm:method and application to the Northern Hayward Fault, California, Bull Seismol Soc Am, 90(6): 1353-1368.
Yu C Q, Zheng Y C, Shang X F, 2017, Crazyseismic:A MATLAB GUI-based software package for passive seismic data preprocessing, Seismol Res Lett, 88(2A): 410-415.