天津滨海爆炸事件检测与相对定位研究

摘要/Abstract

摘要： 2015年8月12日发生在天津滨海新区的化学品爆炸事故造成了严重的人员伤亡和财产损失。基于区域地震台网记录的数字化连续地震波形资料, 对爆炸发生过程的精细分析是深入调查爆炸事故发生原因、快速估计事故所造成灾害与损失的重要手段之一。本文首先利用震相到时和质点运动轨迹信息对记录波形中能量较强的震相性质进行分析, 再通过模板匹配算法检测爆炸过程中是否存在两次大爆炸以外的小爆炸事件, 并对检测方法的效能进行检验, 最后依据波形互相关得到的震相到时差估计两次较大爆炸的相对位置。计算结果显示区域台网记录波形中能量最大的体波震相可能为沿沉积层顶部传播的Pg震相, 垂直向能量最强的面波震相为Rg震相, 而能量较弱的初至震相可能为结晶基底首波; 爆炸过程中除了两次能量较大的爆炸以外还能检测到两次较小的爆炸事件, 其当量相当于约M_L0.5的地震; 基于两次大爆炸Pg震相相对到时信息得到第二次大爆炸发生在第一次大爆炸的北西侧, 两次爆炸位置距离约50~55 m。本文研究结果为爆炸事故的调查分析提供了地震学依据。

关键词: 相对定位, 震相识别, 天津滨海爆炸事件

Abstract: A massive chemical explosion accident hit Tianjin Binhai on August 12, 2015, which caused serious casualties and property losses, and also caused a larger social repercussion. Tianjin seismic network center and China Earthquake Networks Center first reported in an official post on social media that there are two massive explosions, the first of which occurred at 11:34:04 (Beijing time) and only about 30s later another more powerful one occurred nearby, and the preliminary magnitude of these accidental explosions were estimated to be about M_L2.3 and M_L2.9 respectively. In this paper, we first indentified the properties of seismic phase with stronger energy based on the arrival time information and the trajectories of particles. And then detected whether there is other smaller explosion event except the two massive explosions using matched filter technique, and inspected the effects of detection method by 1000 times Bootstrap test. Further we estimated the relative position of two massive explosions using the differences of seismic phase arrival time information obtained from waveform cross-correlation phase detection.
Seismic phase identification results show that the body wave seismic phase has the strongest energy recorded by regional seismic network may be the Pg phase which spread along the top of sediments. And the surface wave phase has the strongest energy recorded in vertical component may be the Rg phase. The first arrival seismic phase might be head wave phase spread along the crystalline basement. We detected two small suspected explosions, which occurred about 2 s before and 84 s after the first massive explosion, and their equivalent approximate to M_L0.5 earthquakes. The relocation result based on Pg phase arrival time differences shows that the second massive explosion is in the north west side of the first one, and the two explosions have about 50~55 m apart. According detailed seismology analysis of the explosion accident, we note that seismic network can play an important role in the process of accident emergency disposal. In the short time after the accident happened, when investigators cannot enter the explosion core region, analysis based on seismic waveform can help the government to understand the detailed process of explosion quickly and take timely and reasonable disposal method.

Key words: Matched filter technique, Relative relocation, Seismic phase identification, Massive explosions in Binhai, Tianjin

中图分类号:

P315.7

邓莉, 谭毅培, 刘双庆, 马婷, 卞真付, 曹井泉. 天津滨海爆炸事件检测与相对定位研究[J]. 地震, 2018, 38(3): 158-169.

DENG Li, TAN Yi-pei, LIU Shuang-qing, MA Ting, BIAN Zhen-fu, CAO Jing-quan. Event Detection and Relative Relocation of the Massive Explosions in Binhai, Tianjin[J]. EARTHQUAKE, 2018, 38(3): 158-169.

参考文献

[1] 沈伟森, 罗艳, 倪四道, 等. 天然地震频率范围内首都圈地区地表S波速度结构[J]. 地震学报, 2010, 32(2): 137-146.
[2] 刘渊源, 崇加军, 倪四道. 基于井下摆天然地震数据测量首都圈近地表波速结构[J]. 地震学报, 2011, 33(3): 342-350.
[3] 武岩, 丁志峰, 朱露培. 利用接收函数研究渤海盆地沉积层结构[J]. 地震学报, 2014, 36(5): 837-849.
[4] Lay T, Wallace T C. Modern Global Seismology[M]. San Diego: Academic Press, 1995: 121-127.
[5] Kafka A L. Rg as depth discriminant for earthquakes and explosion: A case study in New England[J]. Bull Seism Soc Amer, 1990, 80(2): 373-394.
[6] 罗艳, 倪四道, 曾祥方, 等. 一个发生在沉积盖层里的破坏性地震: 2010年1月31日四川遂宁—重庆潼南地震[J]. 科学通报, 2011, 56(2): 147-152.
[7] Peng Z, Zhao P. Migration of early aftershocks following the 2004 Parkfield earthquake[J]. Nature Geoscience, 2009, 2(12): 877-881.
[8] Lengline O, Enescu B, Peng Z, et al. Decay and expansion of the early aftershock activity following the 2011, M_W9.0 Tohoku earthquake[J]. Seismol Res Lett, 2012, 39(18): 18309.
[9] Tang C C, Lin C H, Peng Z. Spatial-temporal evolution of early aftershocks following the 2010 M_L6.4 Jiashian earthquake in southern Taiwan[J]. Geophys J Int, 2014, 199(3): 1772-1783.
[10] Kato A, Obara K. Step-like migration of early aftershocks following the 2007 M_W6.7 Noto-Hanto earthquake, Japan[J]. Geophy Res Lett, 2014, 41(11): 3864-3869.
[11] Wu Q, Chapman M C, Beale J N. The aftershock sequence of the 2011 Mineral, Virginia, earthquake: temporal and spatial distribution, focal mechanisms, regional stress, and the role of Coulomb stress transfer[J]. Bull Seismol Soc Amer, 2015, 105(5): 2521-2537.
[12] Kato A, Obara K, Igarashi T, et al. Propagation of Slow Slip Leading Up to the 2011 M_W 9.0 Tohoku-Oki Earthquake[J]. Science, 2012, 335(6069): 705-708.
[13] Kato A, Nakagawa S. Multiple slow-slip events during a foreshock sequence of the 2014 Iquique, Chile M_W8.1 earthquake[J]. Geophys Res Lett, 2014, 41(15): 5420-5427.
[14] Sugan, M, Kato A, Miyake H, et al. The preparatory phase of the 2009 M_W6.3 L′Aquila earthquake by improving the detection capability of low-magnitude foreshocks[J]. Geophys Res Lett, 2014, 41(17): 6137-6144.
[15] Zhang M, Wen L X. An effective method for small event detection: match and locate (M&L)[J]. Geophys J Int, 2015, 200(3): 1523-1537.
[16] Kato A, Igarashi T, Obara K. Detection of a hidden Boso slow slip event immediately after the 2011 M_W 9.0 Tohoku-Oki earthquake, Japan[J]. Geophys Res Lett, 2014, 41(16): 5868-5874.
[17] 谭毅培, 曹井泉, 刘文兵, 等. 2013年3月涿鹿微震群遗漏地震事件检测和发震构造分析[J]. 地球物理学报, 2014, 57(6): 1847-1856.
[18] Zhao X, Feng W, Tan Y P, et al. Seismological Investigations of Two Massive Explosions in Tianjin, China[J]. Seismol Res Lett, 2016, 87(4): 826-836.
[19] Waldhauser F, Ellsworth W L. A double-difference earthquake location algorithm: method and application to the northern Hayward Fault California[J]. Bull Seism Soc Am, 2000, 90(6): 1353-1368.
[20] Schaff D P, Waldhauser F. Waveform cross-correlation based differential travel-time measurements at the Northern California Seismic Network[J]. Bull Seismol Soc Am, 2005, 95(6): 2446-2461.
[21] Akaike H. A new look at the statistical model identification[J]. IEEE Trans Autom Control, 1974, AC-19(6): 716-723.
[22] 徐锡伟, 吴卫民, 张先康, 等. 首都圈地区地壳最新构造变动与地震[M]. 北京: 科学出版社, 2002, 214-217.
[23] Bormann P. IASPEI New Manual of Seismological Observatory Practice[M]. Potsdam: GeoForschungZetntrum, 2002, 43-46.
[24] 刘瑞丰, 陈翔, 沈道康, 等. 宽频带数字地震记录震相分析[M]. 北京: 地震出版社, 2014.
[25] Li J, Tian B. Refined locations of major explosions in Tianjin harbor[J]. Chinese Sci Bull, 2015, 60(21): 1868-1870.