2019年6月17日四川长宁MS6.0地震早期序列参数分析及强余震概率预测

doi:10.12196/j.issn.1000-3274.2020.02.011

摘要/Abstract

摘要： 2019年6月17日四川长宁发生M_S6.0地震, 之后发生了一系列的强余震, 为更好地分析此次地震的序列特征以及强余震的可预测属性, 采用国际上对复杂序列拟合相对较好的“传染型余震序列”(ETAS)模型以及基于Reseanberg-Jones(R-J)模型发展的Omi-R-J模型, 通过连续滑动、拟合和余震发生率预测, 对地震序列的模型参数稳定性、预测结果进行了比较研究, 并利用N-test、 T-test检验方法对预测结果进行了效能评估。结果表明, 相比于其他中强震序列参数, 此次长宁M_S6.0地震序列参数中反映激发能力的αETAS较其他序列明显偏小, 而反映衰减能力的pORJ值和应力累积水平的bORJ值相对较小, 与此次余震序列丰富、持续时间相对较长相吻合; ETAS和Omi-R-J模型对于复杂序列在[3.0, 3.5, 4.0]三个震级档的强余震仍具有一定的预测能力; 总体的“每个地震的信息增益”(IGPE)计算结果显示, ETAS模型略优于Omi-R-J模型, 前者或更适合复杂地震序列的余震预测。

关键词: ETAS模型, Omi-R-J模型, 序列参数, 余震预测, 效能评估

Abstract: An M_S6.0 earthquake occurred on June 17, 2019 in Changning, Sichuan Province, which was followed by a series of strong aftershocks. In order to better analyze the sequence characteristics and predictability of the strong aftershocks, this paper focuses on the comparative study of the applicability of the “Epidemic Type Aftershock Sequence” (ETAS) model which have been widely adopted internationally for complex sequence, the Omi-R-J model which based on the Reseanberg-Jones (R-J) model. We compare the stability of model parameters, forecasting results by continuous sliding, fitting and aftershock occurrence rate forecasting, and evaluate the effectiveness of forecasting by using the N-test and T-test with multiple time windows. The results show that the sequence parameters αETAS of response excitation ability of Changning M_S6.0 earthquake sequence is obviously smaller than other parameters of moderate-strong earthquake sequence, while the pORJ value of response attenuation ability and the bORJ value of response stress accumulation level are relatively small, which are consistent with the rich and relatively long duration of this aftershock sequence. ETAS and Omi-R-J models still have certain forecasting ability for strong aftershocks of complex sequence at three magnitude ranges [3.0, 3.5, 4.0]. The overall “Information Gain Per Earthquake” (IGPE) calculation results show that the ETAS model is slightly superior to the Omi-R-J model, thus the former is more suitable for aftershock forecasting of complex earthquake sequences.

Key words: Epidemic Type Aftershock Sequence (ETAS) model, Omi-Reseanberg-Jones (Omi-R-J) model, Sequence parameters, Aftershock forecasting, Effectiveness evaluation

中图分类号:

P315.7

毕金孟, 蒋长胜, 马永. 2019年6月17日四川长宁M_S6.0地震早期序列参数分析及强余震概率预测[J]. 地震, 2020, 40(2): 140-154.

BI Jin-meng, JIANG Chang-sheng, MA Yong. The Study on Early Sequence Parameters and Probability Forecasting of Strong Aftershocks of Changning M_S6.0 Earthquake on June 17, 2019, Sichuan Province[J]. EARTHQUAKE, 2020, 40(2): 140-154.

参考文献

[1] Reasenberg P A, Jones L M. Earthquake hazard after a mainshock in California[J]. Science, 1989, 243: 1173-1176.
[2] Gerstenberger M C, Wiemer S, Jones L M, et al. Real-time forecasts of tomorrow’s earthquakes in California[J]. Nature, 2005, 435: 328-331.
[3] Nanjo K Z, Tsuruoka H, Yokoi S, et al. Predictability study on the aftershock sequence following the 2011 Tohoku-Oki, Japan, earthquake: First results[J]. Geophysical Journal International, 2012, 191: 653-658.
[4] Ogata Y, Katsura K, Falcone G, et al. Comprehensive and topical evaluations of earthquake forecasts in terms of number, time, space, and magnitude[J]. Bulletin of the Seismological Society of America, 2013, 103(3): 1692-1708.
[5] Schorlemmer D, Werner M J, Marzocchi W, et al. The Collaboratory for the Study of Earthquake Predictability: Achievements and Priorities[J]. Seismological Research Letters, 2018, 89: 1305-1313.
[6] Japan Meteorological Agency (JMA). The Iwate-Miyagi Nairiku Earthquake in 2008[Z]. Reports Coordinate Community Earthquake Prediction, 2009, 81: 101-131.
[7] Ogata Y. Estimation of the parameters in the modified Omori Formula for aftershock frequencies by the maximum likelihood procedure[J]. Journal of Physics of the Earth, 1983, 31: 115-124.
[8] Enescu B, Mori J, Miyazawa M. Quantifying early aftershock activity of the 2004 mid-Niigata Prefecture earthquake (M_W6.6)[J]. Journal of Geophysical Research, 2007, 112: B04310.
[9] Ogata Y. Statistical models for earthquake occurrences and residual analysis for point processes[J]. Journal of the American Statistical Association, 1988, 83(401): 9-27.
[10] Omi T, Ogata Y, Hirata Y, et al. Forecasting large aftershocks within one day after the main shock[J]. Scientific Reports, 2013, 3: 2218.
[11] Ogata Y, Katsura K. Immediate and updated forecasting of aftershock hazard[J]. Geophysical Research Letters, 2006, 33: L10305.
[12] Omi T, Ogata Y, Shiomi K, et al. Automatic aftershock forecasting: A test using real-time seismicity data in Japan[J]. Bulletin of the Seismological Society of America, 2016, 106(6): 2450-2458.
[13] Omi T, Ogata Y, Shiomi K. Implementation of a real-time system for automatic aftershock forecasting in Japan[J]. Seismological Research Letters, 2019, 90(1): 242-250.
[14] 徐锡伟, 韩竹军, 杨晓平, 等. 中国及邻近地区地震构造图[M]. 北京: 地震出版社, 2016.
XU Xi-wei, HAN Zhu-jun, YANG Xiao-ping, et al. Seismic tectonic maps in China and its adjacent areas[M]. Beijing: Seismological Press, 2016 (in Chinese).
[15] 何登发, 鲁人齐, 黄涵宇, 等. 长宁页岩气开发区地震的构造地质背景[J]. 石油勘探与开发, 2019, 46(3): 1-14.
HE Deng-fa, LU Ren-qi, HUANG Han-yu, et al. Seismic tectonic geological background of Changning shale gas development zone[J]. Petroleum Exploration and Development, 2019, 46(3): 1-14 (in Chinese).
[16] 顾功叙. 中国地震目录: 公元前1831年—公元1969年[M]. 北京: 地震出版社, 1983.
GU Gong-xu. Catalogue of Chinese Earthquakes: 1831BC—1969AD[M]. Beijing: Seismological Press, 1983 (in Chinese).
[17] 蒋长胜, 吴忠良. 玉树M_S7.1地震前的中长期加速矩释放(AMR)[J]. 地球物理学报, 2011, 54(6): 1501-1510.
JIANG Chang-sheng, Wu Zhong-liang. Intermediate-term medium-range Accelerating Moment Release (AMR) priori to the 2010 Yushu M_S7.1 earthquake[J]. Chinese Journal of Geophysics, 2011, 54(6): 1501-1510 (in Chinese).
[18] Ogata Y, Katsura K. Analysis of temporal and spatial heterogeneity of magnitude frequency distribution inferred from earthquake catalogues[J]. Geophysical Journal International, 1993, 113(3): 727-738.
[19] Ogata Y. Detection of precursory relative quiescence before great earthquakes through a statistical model[J]. Journal of Geophysical Research, 1992, 97: 19845-19871.
[20] Ogata Y. Increased probability of large earthquakes near aftershock regions with relative quiescence[J]. Journal of Geophysical Research, 2001, 106: 8729-8744.
[21] Omori F. On aftershocks of earthquakes[J]. Journal of the College of Science. Imperial University of Tokyo, 1894, 7: 11-200.
[22] Utsu T. A statistical study of on the occurrence of aftershocks[J]. Geophysical Magazine, 1961, 30: 521-605.
[23] Gutenberg R, Richter C F. Frequency of earthquakes in California[J]. Bulletin of the Seismological Society of America, 1944, 34: 185-188.
[24] Wiemer S, Katsumata K. Spatial variability of seismicity parameters in aftershocks zones[J]. Journal of Geophysical Research, 1999, 104 (13): 13135-13151.
[25] Enescu B, Enescu D, Ito K. Values of b and p: their variations and relation to physical processes for earthquakes in Japan and Romania[J]. Romanian Journal of Physics, 2011, 56 (3-4): 590-608.
[26] Bi J M, Jiang C S. Comparison of early aftershock forecasting for the 2008 Wenchuan M_S8.0 earthquake[J]. Pure and Applied Geophysics, 2020, 177: 9-25.
[27] 蒋长胜, 庄建仓, 龙锋, 等. 2013年芦山M_S7.0地震序列参数的早期特征: 传染型余震序列模型计算结果[J]. 地震学报, 2013, 35(5): 661-669.
JIANG Chang-sheng, ZHUANG Jian-cang, LONG Feng, et al. Statistical analysis of the early period of the 2013 Lushan M_S7.0 Earthquake Sequence by using the ETAS model[J]. Acta Seismologica Sinica, 2013, 35(5): 661-669 (in Chinese).
[28] 蒋长胜, 吴忠良, 韩立波, 等. 地震序列早期参数估计和余震概率预测中截止震级M_c的影响: 以2013年甘肃岷县漳县6.6级地震为例[J]. 地球物理学报, 2013b, 56(12): 4048-4057.
JIANG Chang-sheng, WU Zhong-liang, HAN Li-bo, et al. Effect of cutoff magnitude M_c of earthquake catalogues on the early estimation of earthquake sequence parameters with implication for the probabilistic forecast of aftershocks: The 2013 Minxian-Zhangxian, Gansu, M_S6.6 earthquake sequence[J]. Chinese Journal of Geophysics, 2013b, 56(12): 4048-4057 (in Chinese).
[29] 蒋长胜, 吴忠良, 尹凤玲, 等. 余震的序列参数稳定性和余震短期发生率预测效能的连续评估以2014年云南鲁甸6.5级地震为例[J]. 地球物理学报, 2015, 58(11): 4163-4173.
JIANG Chang-sheng, WU Zhong-liang, Yin Feng-ling, et al. Stability of early-estimation sequence parameters for continuous forecast of the aftershock rate: A case study of the 2014 Ludian, Yunnan M_S6.5 earthquake[J]. Chinese Journal of Geophysics, 2015, 58(11): 4163-4173 (in Chinese).
[30] Cheng J, Wu Z, Liu J, et al. Preliminary report on the 3 August 2014, M_W 6.2/M_S 6.5 Ludian, Yunnan-Sichuan border, southwest China, earthquake[J]. Seismological Research Letter, 2015, 86(3): 750-763.
[31] 蒋长胜, 韩立波, 郭路杰. 新疆于田地区2008年以来3个地震序列的参数早期特征[J]. 地震学报, 2014, 36(2): 165-174.
JIANG Chang-sheng, HAN Li-bo, GUO Lu-jie. Statistical analysis of the early period of three earthquake sequences in the Xinjiang Yutian region by using the ETAS model[J]. Acta Seismologica Sinica, 2014, 36(2): 165-174 (in Chinese).
[32] Fang L H, Wu J P, Liu J, et al. Preliminary report on the 22 November 2014 M_W6.1/M_S6.3 Kangding earthquake, western Sichuan, China[J]. Seismological Research Letter, 2015, 86(6): 1603-1613.
[33] 蒋长胜, 庄建仓, 吴忠良, 等. 两种短期概率预测模型在2017年九寨沟7.0级地震中的应用和比较研究[J]. 地球物理学报, 2017, 60(10): 4132-4144.
JIANG Chang-sheng, ZHUANG Jian-cang, WU Zhong-liang, et al. Application and comparison of two short-term probabilistic forecasting models for the 2017 Jiuzhaigou, Sichuan, M_S7.0 earthquake[J]. Chinese Journal of Geophysics, 2017, 60(10): 4132-4144 (in Chinese).
[34] 毕金孟, 蒋长胜. 华北地区地震序列参数的分布特征[J]. 地球物理学报, 2019, 62(11): 4300-4312.
BI Jin-meng, JIANG Chang-sheng. Distribution characteristics of earthquake sequence parameters in North China[J]. Chinese Journal of Geophysics, 2019, 62(11): 4300-4312 (in Chinese).
[35] 毕金孟, 蒋长胜. 全球部分强震序列参数的分布特征[J]. 地震, 2019.
BI Jin-meng, JIANG Chang-sheng. Distribution characteristics of partial strong earthquake sequence parameters in the world[J]. Earthquake, 2019 (in Chinese).
[36] 蒋长胜, 毕金孟, 王福昌, 等. 利用早期余震预测的Omi-R-J方法对2017年四川九寨沟M_S7.0地震的应用研究[J]. 地球物理学报, 2018, 61(5): 2099-2110.
JIANG Chang-sheng, BI Jin-meng, WANG Fu-chang, et al. Application of the Omi-R-J method for forecast of early aftershocks to the 2017 Jiuzhaigou, Sichuan, M_S7.0 earthquake[J]. Chinese Journal of Geophysics, 2018, 61(5): 2099-2110 (in Chinese).
[37] Kagan Y Y, Jackson D D. New seismic gap hypothesis: Five years after[J]. Journal of Geophysical Research, 1995, 100(B3): 3934-3959.
[38] Schorlemmer D, Gerstenberger M C, Wiemer S, et al. Earthquake likelihood model testing[J]. Seismological Research Letters, 2007, 78: 17-29.
[39] Zechar J D. Evaluating earthquake predictions and earthquake forecasts: a guide for students and new researchers, in: Community Online Resource for Statistical Seismicity Analysis[Z]. 2010, doi:10.5078/corssa-77337879.
[40] Imoto M. Information gain of a model based on multidisciplinary observations with correlations[J]. Journal of Geophysical Research, 2007, 112: B05306.