Spatio-temporal characteristics of b value in the northern margin of Ordos block and its earthquake prediction efficiency
-
摘要: 使用中国地震台网1970—2017年地震目录,在评估目录完整性和删除余震的基础上,使用最大似然法对鄂尔多斯地块北缘进行b值时空扫描. 扫描结果显示,鄂尔多斯地块北缘的低b值区域主要分布于阿拉善和晋冀蒙地区,低b值区域的空间分布随时间呈现一定变化. 低b值区域与中强地震存在一定的对应关系,对14次震例的回溯性检验显示,其中10次地震发生在低b值区域及其边缘. 使用Molchan图表法和R值评分法对b值时空扫描结果进行预测效能检验,以低于阈值的区域为预测区域,计算不同阈值的漏报率、报准率和时空占有率。Molchan图表法检验结果显示,本文所使用的b值时空扫描方法的预测效果优于泊松分布,最优结果的概率增益达到2,显著性水平达到0.25. R值曲线显示,当阈值为0.8时,达到最高R值0.28,并达到97.5%置信水平.
-
关键词:
- b值 /
- 时空扫描 /
- 预测效能 /
- Molchan图表法 /
- 鄂尔多斯地块北缘
Abstract: Using earthquake catalogue of China Earthquake Networks Center from 1970 to 2017, the b values of the north edge of Ordos block were scanned by the maximum likelihood method after evaluating completeness and deleting aftershocks. The results show that the low b value region of the north margin of Ordos block is mainly distributed in the Alxa area and the Shanxi-Hebei-Inner Mongolia area, and the spatial distribution of low b value region changes with time to a certain extent. Meanwhile, there is a certain corresponding between the low b value region and the location of moderate-strong earthquakes. The retrospective test of 14 earthquake cases shows that ten earthquakes occurred in the low b value regions and their edges, and four earthquakes occurred beyond the low b value region. And then Molchan chart method and R value scoring method were used to test the prediction efficiency of b value spatio-temporal scanning. Taking the area below the threshold as the predicted area, the missing rate, the hit rate and the spatio-temporal occupation rate of different thresholds were calculated. The result of Molchan chart test shows that the prediction efficiency of spatio-temporal scanning is better than that of Poisson distribution, the probability gain of the optimal result reached 2, and the significance level reached 0.25. The R value curve shows that the highest R value is 0.28 at the threshold of 0.8 with 97.5% confidence level. -
-
表 1 鄂尔多斯地块北缘14次震例的回溯情况
Table 1 Retrospective test for the fourteen earthquake examples in the northern margin of Ordos block
编号 年-月-日 北纬/° 东经/° MS 地点 低b值 1 1976-04-06 40.20 112.10 6.2 和林格尔 × 2 1976-09-23 39.90 106.40 6.2 巴音木仁 √ 3 1979-08-25 41.20 108.10 6.0 五原 √ 4 1981-08-13 40.58 113.41 5.6 丰镇 √ 5 1989-10-19 39.92 113.91 5.9 山西大同 √ 6 1991-01-13 40.40 105.60 5.2 阿拉善左旗 × 7 1991-03-26 40.00 113.80 5.8 山西大同 √ 8 1991-09-14 40.20 105.50 5.1 阿拉善左旗 × 9 1996-05-03 40.80 109.60 6.4 包头 √ 10 1997-10-21 41.20 107.10 5.0 乌拉特后旗 √ 11 1998-01-10 41.10 114.30 6.2 河北张北 × 12 1999-03-11 41.20 114.60 5.6 河北张北 √ 13 1999-11-01 39.80 113.90 5.6 山西浑源 √ 14 2015-04-15 39.80 106.30 5.8 阿拉善左旗 √ -
曹刚. 2001. 内蒙古地震研究[M]. 北京: 地震出版社: 5–6. Cao G. 2001. Seismological Research of Inner Mongolia[M]. Beijing: Seismological Press: 5–6 (in Chinese).
邓起东, 程绍平, 闵伟, 杨桂枝, 任殿卫. 1999. 鄂尔多斯块体新生代构造活动和动力学的讨论[J]. 地质力学学报, 5(3): 13-21. Deng Q D, Cheng S P, Min W, Yang G Z, Ren D W.1999.Discussion on Cenozoic tectonics and dynamics of ordos block[J]. Journal of Geomechanics, 5(3): 13-21(in Chinese)
韩晓明, 张文韬, 王树波, 赵星, 包金哲, 李娟. 2016. 河套地震带的b值时空变化特征分析[J]. 中国地震, 32(3): 522-532. Han X M, Zhang W T, Wang S B, Zhao X, Bao J Z, Li J. 2016. Analysis of spatial and temporal variation characteristics of b-value in the Hetao seismic belt[J]. Earthquake Research in China, 32(3): 522-532 (in Chinese).
蒋长胜, 张浪平, 韩立波, 来贵娟. 2011. 中长期地震危险性概率预测中的统计检验方法Ⅰ: Molchan图表法[J]. 地震, 31(2): 106-113. Jiang C S, Zhang L P, Han L B, Lai G J. 2011. Probabilistic forecasting method of long-term and intermediate-term seismic hazard I: Molchan error diagram[J]. Earthquake, 31(2): 106-113 (in Chinese).
李全林, 陈锦标, 于渌, 郝柏林. 1978. b值时空扫描: 监视破坏性地震孕育过程的一种手段[J]. 地球物理学报, 21(2): 101-125. Li Q L, Chen J B, Yu L, Hao B L. 1978. Time and space scanning of the b-value: A method for monitoring the development of catastrophic earthquakes[J]. Chinese Journal of GeophysicsSinica, 21(2): 101-125 (in Chinese).
李正芳, 周本刚. 2014. 利用断裂带上的低b值识别凹凸体方法的探讨:以龙门山断裂带和鲜水河断裂带为例[J]. 2014. 震灾防御技术, 9(2): 213-225. Li Z F, Zhou B G. 2014. Asperity identification based on low b-value: Application to the Longmenshan and Xianshuihe fault zone[J]. Technology for Earthquake Disaster Prevention, 9(2): 213-225 (in Chinese).
马鸿庆. 1978. 华北地区几次大震前的b值异常变化[J]. 地球物理学报, 21(2): 126-141. Ma H Q. 1978. Variations of the b-values before several large earthquakes occurred in north china[J]. Chinese Journal of Geophysics, 21(2): 126-141 (in Chinese).
梅世蓉. 1993. 中国地震预报概论[M]. 北京: 地震出版社: 65–69. Mei S R. 1993. Introduction to Earthquake Prediction in China[M]. Beijing: Seismological Press: 65–69 (in Chinese).
孙加林. 1985. 鄂尔多斯块体北部边缘发生强震的危险性[J]. 西北地震学报, 7(增刊1): 13-23. Sun J L. 1985. Strong earthquake assessment in the northern edge of Ordos[J]. Northwestern Seismological Journal, 7(S1): 13-23 (in Chinese).
孙丽娜, 齐玉妍, 温超, 张合. 2012. 中长期地震预测中的PI算法改进研究及应用[J]. 地震, 32(4): 44-52. Sun L N, Qin Y Y, Chao Wen CW, Zhang H. 2012. Application of improved PI algorithm in medium to long-term earthquake prediction[J]. Earthquake, 32(4): 44-52 (in Chinese).
王熠熙, 张辉, 刘双庆, 曹井泉, 李赫, 陈嵩, 高也. 2015. 河北平原地震带b值时空变化特征[J]. 地震工程学报, 37(1): 188-195. Wang Y X, Zhang H, Liu S Q, Cao J Q, Li H, Chen S, Gao Y. 2015. Spatial and temporal variations of b-value in the Hebei plain seismic belt[J]. China Earthquake Engineering Journal, 37(1): 188-195 (in Chinese).
许绍燮. 1989. 地震预报能力评分[G]//地震预报方法实用化研究文集: 地震学专辑. 北京: 学术书刊出版社: 586–590. Xu S X. 1989. Earthquake prediction ability score[G]//Earthquake Prediction Method Practical Research Anthology: Seismology Album. Beijing: Academic Book Press: 586–590 (in Chinese).
徐锡伟, 程国良, 马杏垣, 孙宇航, 韩竹君. 1994. 华北及其邻区块体转动模式和动力来源[J]. 地球科学: 中国地质大学学报, 19(2): 129-138. Xu X W, Cheng G L, Ma X Y, Sun Y H, Han Z J. 1994. Rotation model and dynamics of blocks in North China and its adjacent areas[J]. Earth Science: Journal of China University of Geosciences, 19(2): 129-138 (in Chinese).
薛艳, 陈琳荣, 黎明晓, 刘双庆. 2017. 华北地区震群活动与强震[J]. 华北地震科学, 35(2): 1-9. Xue Y, Chen L R, Li M X, Liu S Q. 2017. Relationship between earthquake swarms and the strong earthquake in North China[J]. North China Earthquake Sciences, 35(2): 1-9 (in Chinese).
易桂喜, 闻学泽, 张致伟, 龙锋, 阮祥, 杜方. 2010. 川南马边地区强震危险性分析[J]. 地震地质, 32(2): 282-293. Yi G X, Wen X Z, Zhang Z W, Long F, Ruan X, Du F. 2010. Study on potential strong earthquake risks inAround the Mabian area, Southern Sichuan[J]. Seismology and Geology, 32(2): 282-293 (in Chinese).
易桂喜, 闻学泽, 辛华, 乔慧珍, 龙锋, 王思维. 2011. 2008年汶川MS8.0地震前龙门山—岷山构造带的地震活动性参数与地震视应力分布[J]. 地球物理学报, 54(6): 1490-1500. Yi G X, Wen X Z, Xin H, Qiao XH Z, Long F, Wang S W. 2011. Distributions of seismicity parameters and seismic apparent stresses on the Longmenshan-Minshan tectonic zone before the 2008 M(S)8.0 Wenchuan earthquake[J]. Chinese Journal of Geophysics, 54(6): 1490-1500 (in Chinese).
张琳琳, 唐兰兰. 2015. 新疆天山地区b值时空变化特征[J]. 地震, 35(3): 93-101. Zhang L L, Tang L L. 2015. Temporal and spatial variations of b-value in Tianshan, Xinjiang[J]. Earthquake, 35(3): 93-101 (in Chinese).
郑兆苾, 刘东旺, 沈小七,汪雪泉. 2001. 华北b值全时空扫描结果的可靠性及与地震的相关性[J]. 地震, 21(3):8-14. Zheng Z B, Liu D W, Shen X Q, Wang X Q. 2001. Reliability of time-space scanning results of b value and correlativity between b vales and earthquakes with M≥5.0 in North China[J]. Earthquake, 21(3):8-14 (in Chinese).
Aki K. 1984. Asperities, barriers, characteristic earthquakes and strong motion prediction[J]. J Geophys Res, 89(B7): 5867-5872. doi: 10.1029/JB089iB07p05867
Amitrano D. 2003. Brittle-ductile transition and associated seismicity: Experimental and numerical studies and relationship with the b value[J]. J Geophys Res, 108(B1): 2044.
Bak P, Tang C. 1989. Earthquakes as a self-organized critical phenomenon[J]. J Geophys Res, 94(B11): 15635-15637. doi: 10.1029/JB094iB11p15635
Chen S, Jiang C S, Zhuang J C. 2016. Statistical evaluation of efficiency and possibility of earthquake predictions with gravity field variation and its analytic signal in western China[J]. Pure Appl Geophys, 173(1): 305-319. doi: 10.1007/s00024-015-1114-x
Gutenberg B, Richter C F. 1954. Seismicity of the Earth and Associated Phenomena[M]. 2nd ed. Princeton: Princeton University Press: 1–310.
Ito K, Matsuzaki M. 1990. Earthquakes as self-organized critical phenomena[J]. J Geophys Res, 95(B5): 6853-6860. doi: 10.1029/JB095iB05p06853
Kamer Y, Hiemer S. 2015. Data-driven spatial b value estimation with applications to California seismicity: to b or not to b[J]. J Geophys Res, 120(7): 5191-5214. doi: 10.1002/2014JB011510
Molchan G M. 1997. Earthquake prediction as a decision-making problem[J]. Pure Appl Geophys, 149(1): 233-247. doi: 10.1007/BF00945169
Mori J, Abercrombie R E. 1997. Depth dependence of earthquake frequency-magnitude distributions in California: implications for rupture initiation[J]. J Geophys Res, 102(B7): 15081–15090. doi: 10.1029/97JB01356
Reasenberg P. 1985. Second‐order moment of central California seismicity, 1969–1982[J]. J Geophys Res, 90(B7): 5479-5495. doi: 10.1029/JB090iB07p05479
Scholz C H. 1968. The frequency-magnitude relation of microfracturing in rock and its relation to earthquakes[J]. Bull Seismol Soc Am, 58(1): 399-415.
Schorlemmer D, Wiemer S, Wyss M. 2005. Variations in earthquake-size distribution across different stress regimes[J]. Nature, 437(7058): 539-542. doi: 10.1038/nature04094
Urbancic T I, Trifu C I, Long J M, Young R P. 1992. Space-time correlations of b values with stress release[J]. Pure Appl Geophys, 139(3/4): 449-462.
Utsu T. 1978. Estimation of parameters in formulas for frequency-magnitude relation of earthquake occurrence: in cases involving a parameter c for the maximum magnitude[J]. Zisin, 31(4): 367-382. doi: 10.4294/zisin1948.31.4_367
Wiemer S, Wyss M. 1997. Mapping the frequency‐magnitude distribution in asperities: An improved technique to calculate recurrence times[J]. J Geophys Res, 102(B7): 15115-15128. doi: 10.1029/97JB00726
Wiemer S, Schorlemmer D. 2007. ALM: An asperity-based likelihood model for California[J]. Seismol Res Lett, 78(1): 134-140. doi: 10.1785/gssrl.78.1.134
Woessner J, Wiemer S. 2005. Assessing the quality of earthquake catalogues: estimating the magnitude of completeness and its uncertainty[J]. Bull Seismol Soc Am, 95(2): 684-698. doi: 10.1785/0120040007
Wyss M. 1973. Towards a physical understanding of the earthquake frequency distribution[J]. Geophys J Int, 31(4): 341-359. doi: 10.1111/j.1365-246X.1973.tb06506.x
Wyss M, McNutt S R. 1998. Temporal and three-dimensional spatial analyses of the frequency-magnitude distribution near Long Valley Caldera, California[J]. Geophys J Int, 134(2): 409-421. doi: 10.1046/j.1365-246x.1998.00561.x
Wyss M, Wiemer S. 2000. Change in the probability for earthquakes in southern California due to the Landers magnitude 7.3 earthquake[J]. Science, 290(5495): 1334-1338. doi: 10.1126/science.290.5495.1334
Wyss M. 2004. Fractal dimension and b-value on creeping and locked patches of the san Andreas fault near Parkfield, California[J]. Bull Seismol Soc Am, 94(2): 410-421. doi: 10.1785/0120030054
Zechar J D, Jordan T H. 2008. Testing alarm-based earthquake predictions[J]. Geophys J Int, 172(2): 715-724. doi: 10.1111/gji.2008.172.issue-2