2013年4月20日8时02分(北京时间),在龙门山断裂带的大川-双石断裂上发生了四川芦山7.0级强震(图 1).此次地震给芦山县及其周边地区造成重大的人员伤亡和财产损失,截至2013年4月22日18时,已经造成192人死亡,23人失踪,11470人受伤,其中重伤968人(中国地震信息网).至2013年4月22日10时,中国地震台网中心发布已有2089次余震,其中3级以上91次、4级以上20次、5级以上4次,最大余震5.4级.
芦山7.0级强震发生在龙门山构造带上,是2008年汶川8.0级地震以来震级最大的一次地震.龙门山构造带是松潘-甘孜地块的东边界,主要由3条大断裂构成,分别是龙门山前山断裂、中央主断裂和后山断裂,此次芦山地震震中位置位于龙门山前山断裂西南段的大川-双石断裂(DSF)上(图 1).大川-双石断裂为龙门山推覆构造带南段一条区域性大断裂,野外地质调查结果表明,其走向N43°E、倾向NW、倾角45°~65°不等,其东南侧为开阔的中新生代陆相盆地,西北侧为古生代地层组成的中高山区[2].地震发生后,不同的研究机构均给出主震的震源参数.哈佛大学矩张量结果显示,芦山7.0级强震断层面倾角约38°(www.globalcmt.org);中国科学院地质与地球物理研究所的震源破裂过程反演结果显示,断层面倾角约38.5°[3].为进一步深入认识本次芦山地震断层面倾角及余震扩展空间范围,本研究收集了主震后48小时余震到时资料,通过双差定位,给出较为精确余震空间分布,为今后判定余震活动趋势提供重要的参考.
本研究采用中国地震台网中心提供的2013年4月20日8时至4月22日8时,即芦山7.0级主震后48小时内的观测报告.截至该时间,主震区尽管发生了近2000次余震,但从台网中心能够下载到的地震事件共504个,其中M≥2级435个,占总事件的86%.我们选择震中距小于400 km的到时资料,其中P波5418条、S波4817条,数据集包括891条Pn到时.Pn波能够较好地约束震源深度,结合近台直达波数据,可以较好地改善震源深度,提高重定位精度[4-5].本研究用于重定位的地震台站总共44个,均来自四川台网,可以看出这些台站围绕余震有着较好的空间分布,特别是在震源区还有3个近台(图 1),能够进一步提高震源深度的精度.
本研究采用双差定位法[6]对余震序列进行精确定位.该方法使用相对走时残差来修定地震位置,其主要是基于两个震源之间的距离远小于事件到台站的距离,认为两个事件传播到台站的射线路径几乎相同,从而可以有效地降低由于速度模型的不确定性对定位结果的影响.本研究采用适合本研究区的一维速度模型[7](图 2)来开展双差定位,其中P波和S波速度比为1.71.由于S波到时资料拾取精度相对P波要低些,因此在定位过程中将P波权重设为1.0,而S波权重设为0.5.为了获得高精度的定位结果,设定事件对的距离阀值为3 km.该方法已在地震重定位工作中得到广泛应用,且获得较好的定位结果[8-11].
通过双差定位,我们最终获得328次地震事件的精定位结果,其均方根残差由重定位前的0.54 s下降到重定位后的0.17 s.由误差棒分布可以清楚看出,东西、南北及深度上平均标准误差分别为0.35 km、0.39 km和0.48 km,80%地震事件的三方向误差均在0.25 km以内,其中误差大于1.0 km的地震事件只有15个(图 3),说明重定位结果有意义地改进了初始震源位置精度.本次芦山7.0级主震重定位震中位置为102.96°E、30.29°N、20.1 km深,其东西、南北及深度上误差分别为0.23 km、0.28 km和0.48 km.
为进一步清楚了解重定位的效果,我们比较一下定位前后结果(图 4)的差别.重定位前,地震震中在平面图上零散地分布在断裂带两侧,在深度方向上震源分布也很零散,约在2~30 km较大深度范围内,看不出优势发震层(图 4a);重定位后,地震震中在平面图上明显地集中在断层两侧,在深度方向上也显示出明显的优势发震层(18~22 km)(图 4b).另外,余震主要分布在主震西南向、且较为集中,表明破裂主要沿西南向扩展,这一结果得到矩张量解及断层面滑动角结果显示的破裂面往西南向发展[3](www.globalcmt.org)结果的支持.从重定位后主震震源深度变化来看,7.0级主震震源深度从初始的17 km(中国地震台网中心提供的观测报告)改进到20.1 km.这一结果尽管与房立华等(www.cea-igp.ac.cn)重定位结果显示主震深度为12.3 km存在一定差别,但与张勇等[12]和哈佛大学矩张量结果反演震源破裂深度(约21.8 km)及韩立波等利用台阵数据采用CAP方法反演得到主震矩心深度(约24 km)(www.cea-igp.ac.cn)较为接近.
由于余震能够直观地展示强震发生的深部结构及断层的破裂延展范围,为此本研究详细展示余震沿断层破裂方向及垂直向深度剖面(图 5).沿着断裂带方向,余震分布在时间尺度上尽管看不出明显扩展过程,与2011年盈江5.8级地震[11]有些不同,这可能与本研究所用资料的时间较短有关,需进一步增加资料来证实.然而,在空间分布上主震西南方向余震分布较广、且较为集中,可能表明余震主要向西南方向扩展(AA′剖面).穿过主震垂直于断层的纵剖面显示出明显的铲状形态(BB′剖面),尽管主震西南的纵剖面依然保持了其铲状结构形态,但倾角变大(CC′剖面).这种随着余震的扩展断层面逆冲角度发生变化的现象,非常类似于汶川8.0级地震错断中央断层的倾角变化[13].为进一步定量了解断层面倾角的变化,本研究采用重定位后的余震来拟合发震断层倾角的变化,其中假定了发震断层可以用一个平面来模拟[14-15],因此其倾角为所用到的地震空间范围内的平均倾角.结果表明,大川-双石断裂主震处断层面倾角约44°(图 5BB′剖面),与曾祥方等[16]采用P波初动反演得到的断层面结果44°相一致,而且与哈佛大学矩张量反演倾角38°也较为接近;在主震西南CC′剖面处,大川-双石断层面倾角约73°,可能表明断层破裂往西南向扩展过程中走滑分量逐渐增强.
为试图了解余震是否存在明显的时空变化特征,我们以主震后每5个小时为时间段来展示余震空间分布(图 6).第一阶段为主震发生后0~5小时,余震主要往西南向扩展约22 km、东北向仅8 km左右,难以看出明显的余震随时间扩展规律,可能与本次主震造成断层破裂长度较短有关.震源深度主要集中在16~22 km范围内.该阶段余震震级相对较大,主要为M≥3级地震,最大余震为5.4级.第二阶段为主震发生后5~10小时,该阶段余震随时间扩展变化的规律也不明显,但范围在增大,整个扩展长度达约33 km,主要表现为东北向扩展约3 km.在该阶段,部分余震向浅部运移,最浅达10 km左右,同时也有地震向深部发展,最深可达25 km左右,与此相应的余震震级也有所降低.第三阶段为主震发生后10~15小时,该阶段主震西南向余震基本在前两个时间段内发生过的空间范围内发生,但东北向余震空间范围继续扩展约2 km.该阶段余震数量明显减少,震级也明显降低,主要为M≤2级地可达35 km左右,但这还需要收集更长时间的余震资料进一步确认.
4 结论与讨论采用双差定位法,结合近台及Pn波资料,获得芦山7.0级强震及48小时内余震序列精定位结果.结果显示,余震主要沿大川-双石断裂向西南方向扩展兼东北向扩展,扩展长度达35 km,震源优势分布深度在18~22 km之间.余震清楚地展示出主震的铲状断层破裂面,且断层面在不同位置其倾斜角度有所不同,主震处大川-双石断裂的断层面倾角约44°,往西南向逐渐增大.
芦山7.0级强震发生在龙门山构造带,是2008年汶川8.0级地震以来震级最大的一次地震.汶川8.0级地震发生在龙门山断裂带的中段,而此次芦山地震发生在西南段,两次地震的震中相距约90 km(图 1).此次芦山地震是汶川地震的强余震还是一次独立的地震?目前还存在着争议.有些研究者认为此次芦山地震可视为汶川地震余震[3],而另一些研究者则认为该次地震不是汶川地震余震,但两者关系密切[17].从破裂范围及方向上看,汶川地震造成近350 km长的断层破裂[8, 18],余震主要往北东向延展;而芦山地震破裂方向相反,往西南向,破裂长度也相对较短.本文认为此次芦山7.0级地震可能意味着补充了汶川地震的西南方向破裂,这与2004年12月26日的印尼Mw9.0级地震的余震主要向北北西向延展、但在2005年3月29日发生在相同断裂带上的MW8.6级地震的余震向南南东向的扩展补充了南东向的破裂[19-20]相类似.从前人的地震波速度异常来看,此次芦山地震和2008年汶川地震恰好位于明显的低波速和高泊松比异常体南端和北端[21-22],因此认为此次芦山地震的发生也与印度板块和欧亚板块相互作用形成的高温高压深部流体作用于大川-双石断裂带密切有关,这可能和1605年7月13日海南琼山7.5级地震及2011年3月10日云南盈江5.8级地震的发生与深部流体作用于断裂带相类似[23-24].关于此次芦山地震和2008年汶川地震之间的关系还有待进一步深入研究.
[1] | 邓起东, 张培震, 冉勇康, 等. 中国活动构造基本特征. 中国科学(D辑) , 2002, 32(12): 1020–1030. Deng Q D, Zhang P Z, Ran Y K, et al. Basic characteristics of active tectonics of China. Science in China (Series D) (in Chinese) , 2002, 32(12): 1020-1030. |
[2] | 杨晓平, 蒋溥, 宋方敏, 等. 龙门山断裂带南段错断晚更新世以来地层的证据. 地震地质 , 1999, 4: 341–345. Yang X P, Jiang P, Song F M, et al. The evidence of the south Longmengshan fault zones cutting late Quaternary stratum. Seismol. Geol. (in Chinese) , 1999, 4: 341-345. |
[3] | 王卫民, 郝金来, 姚振兴. 2013年4月20日四川芦山地震震源破裂过程反演初步结果. 地球物理学报 , 2013, 56(4): 1412–1417. Wang W M, Hao J L, Yao Z X. Preliminary result for rupture process of Apr. 20, 2013, Lushan Earthquake, Sichuan, China. Chinese J. Geophys. (in Chinese) , 2013, 56(4): 1412-1417. |
[4] | 张天中, 武巴特尔, 黄媛, 等. 近台资料对近震相对定位算法的影响. 地球物理学报 , 2007, 50(4): 1123–1130. Zhang T Z, Wu Bater, Huang Y, et al. Effect of the data recorded at nearby stations on earthquake relative location. Chinese J. Geophys. (in Chinese) , 2007, 50(4): 1123-1130. |
[5] | 黄嫒, 吴建平, 张天中, 等. 汶川8.0级大地震及其余震序列重定位研究. 中国科学(D辑) , 2008, 38(10): 1242–1249. Huang Y, Wu J P, Zhang T Z, et al. Relocation of the 8.0 Wenchuan earthquake and its aftershock sequence. Science in China (Series D) (in Chinese) , 2008, 38(10): 1242-1249. |
[6] | Waldhauser F, Ellsworth W. A double-difference earthquake location algorithm:method and application to the Northern Hayward Fault, California. Bull. Seismol. Soc. Am. , 2000, 90(6): 1353-1368. DOI:10.1785/0120000006 |
[7] | 赵珠, 范军, 郑斯华, 等. 龙门山推覆构造带地壳速度结构和震源位置的精确修订. 地震学报 , 1997, 19(6): 615–622. Zhao Z, Fan J, Zheng S H, et al. Precision determination of the crustal structure and hypocentral locations in the Longmenshan thrust nappe belt. Acta Seismologica Sinica (in Chinese) , 1997, 19(6): 615-622. |
[8] | 陈九辉, 刘启元, 李顺成, 等. 汶川Ms8.0地震余震序列重新定位及其地震构造研究. 地球物理学报 , 2009, 52(2): 390–397. Chen J H, Liu Q Y, Li S C, et al. Seismotectonic study by relocation of the Wenchuan Ms8.0 earthquake sequence. Chinese J.Geophys. (in Chinese) , 2009, 52(2): 390-397. DOI:10.1002/cjg2.v52.2 |
[9] | 张广伟, 雷建设, 谢富仁, 等. 华北地区小震精定位及构造意义. 地震学报 , 2011, 33(6): 699–714. Zhang G W, Lei J S, Xie F R, et al. Precise relocation of small earthquakes occurred in North China and its tectonic implication. Acta Seismologica Sinica9 (in Chinese) , 2011, 33(6): 699-714. |
[10] | 房立华, 吴建平, 张天中, 等. 2011年云南盈江Ms5.8地震及其余震序列重定位. 地震学报 , 2011, 33(2): 262–267. Fang L H, Wu J P, Zhang T Z, et al. Relocation of mainshock and aftershocks of the 2011 Yingjiang Ms5.8 earthquake in Yunnan. Acta Seismologica Sinica (in Chinese) , 2011, 33(2): 262-267. |
[11] | Lei J S, Zhang G W, Xie F R, et al. Relocation of the 10 March 2011 Yingjiang, China, earthquake sequence and its tectonic implications. Earthq. Sci. , 2012, 25: 103-110. DOI:10.1007/s11589-012-0836-4 |
[12] | 张勇, 许力生, 陈运泰. 芦山4.20地震破裂过程及其致灾特征初步分析. 地球物理学报 , 2013, 56(4): 1408–1411. Zhang Y, Xu L S, Chen Y T. Rupture process of the Lushan 4.20 earthquake and preliminary analysis on the disaster-causing mechanism. Chinese J. Geophys. (in Chinese) , 2013, 56(4): 1408-1411. |
[13] | 张勇, 许力生, 陈运泰. 2008年汶川大地震震源机制的时空变化. 地球物理学报 , 2009, 52(2): 379–389. Zhang Y, Xu L S, Chen Y T. Spatio-temporal variation of the source mechanism of the 2008 great Wenchuan earthquake. Chinese J. Geophys. (in Chinese) , 2009, 52(2): 379-389. DOI:10.1002/cjg2.v52.2 |
[14] | 万永革, 沈正康, 刁桂苓, 等. 利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用. 地球物理学报 , 2008, 51(3): 569–583. Wan Y G, Shen Z K, Diao G L, et al. An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence. Chinese J. Geophys. (in Chinese) , 2008, 51(3): 569-583. DOI:10.1002/cjg2.v51.3 |
[15] | 王福昌, 万永革, 胡顺田. 粒子群算法在主震断层面参数估计中的应用. 地震研究 , 2008, 31(2): 149–154. Wang F C, Wan Y G, Hu S T. Application of particle swarm optimization to the estimation of mainshock fault plane parameters. Journal of Seismological Research (in Chinese) , 2008, 31(2): 149-154. |
[16] | 曾祥方, 罗艳, 韩立波, 等. 2013年4月20日四川芦山Ms7.0地震:一个高角度逆冲地震. 地球物理学报 , 2013, 56(4): 1418–1424. Zeng X F, Luo Y, Han L B, et al. The Lushan Ms7.0 earthquake on 20 April 2013:A high-angle thrust event. Chinese J. Geophys. (in Chinese) , 2013, 56(4): 1418-1424. |
[17] | 刘杰, 易桂喜, 张致伟, 等. 2013年4月20日四川芦山M7.0级地震介绍. 地球物理学报 , 2013, 56(4): 1404–1407. Liu J, Yi G X, Zhang Z W, et al. Introduction to the Lushan, Sichuan M7.0 earthquake on 20 April 2013. Chinese J.Geophys. (in Chinese) , 2013, 56(4): 1404-1407. |
[18] | 徐锡伟, 陈桂华, 于贵华, 等. 5.12汶川地震地表破裂基本参数的再论证及其构造内涵分析. 地球物理学报 , 2010, 53(10): 2321–2336. Xu X W, Chen G H, Yu G H, et al. Reevaluation of surface rupture parameters of the 5.12 Wenchuan earthquake and its tectonic implication for Tibetan uplift. Chinese J.Geophys. (in Chinese) , 2010, 53(10): 2321-2336. |
[19] | Engdahl E, Villasen A, DeShon H, et al. Teleseismic relocation and assessment of seismicity (1918-2005) in the region of the 2004Mw9.0 Sumatra-Andaman and 2005Mw8.6 Nias Island Great Earthquakes. Bull. Seismol. Soc. Am. , 2007, 97: 43-61. DOI:10.1785/0120050614 |
[20] | Pesicek J, Thurber C, Zhang H, et al. Teleseismic double-difference relocation of earthquakes along the Sumatra-Andaman subduction zone using a 3-D model. J. Geophys. Res. , 2010, 115. |
[21] | 雷建设, 赵大鹏, 苏金蓉, 等. 龙门山断裂带地壳精细结构与汶川地震发震机理. 地球物理学报 , 2009, 52(2): 339–345. Lei J S, Zhao D P, Su J R, et al. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake. Chinese J.Geophys. (in Chinese) , 2009, 52(2): 339-345. |
[22] | Lei J S, Zhao D P. Structure heterogeneities of the Longmenshan fault zone and the mechanism of the 2008 Wenchuan earthquake (Ms8.0). Geochemistry Geophysics Geosystems , 2009, 10. |
[23] | Lei J S, Zhao D P, Steinberger B, et al. New seismic constraints on the upper mantle structure of the Hainan plume. Phys. Earth Planet. Int. , 2009, 173: 33-50. DOI:10.1016/j.pepi.2008.10.013 |
[24] | Lei J S, Xie F R, Mishar O P, et al. The 2011 Yingjiang, China, earthquake:A volcano-related fluid-driven earthquake?. Bull Seismol. Soc. Am. , 2012, 102: 417-425. DOI:10.1785/0120110143 |