接前文。刚到秦四清研究员博客看了一下,秦研究员不愧是科学网开博客的,比原来打过交道的民科好玩多了,预报了那么多的地震,时空强三要素全给出,还曾经到地震局去传经送宝,他的预报结果咱要是不再多分析两个、那就太不重视他的“学术研究结果”了。
秦研究员再上一个预报是5月11日甘肃西部4.9级地震之后做出的。该地震震中是北纬37.7度、东经102.1度。他的预报是:
“该震是一个强震发生前的标志性foreshock事件,快则10天内,慢则3个月内,该地震区的强震将发生。此外,把预测震中位置调整为该Ms4.9级的位置附近。”
秦研究员把时间和震级都给了(震级,我们把6.0以上就算强震);空间我们算经度两度(东经101到103度)、纬度两度(北纬37到39度),覆盖区域超过170公里×200公里。现在时空强都有了,检验一下,秦研究员测准的概率有多大?
log10(每年次数) = 3.4503 - 0.8994 × 震级
这个斜率0.8994就是这个地区的古登堡—里克特公式特征值b。
如果没有前震因素,将震级=6级带入,得到平均每年发生的次数是0.01131次;考虑到秦四清研究员预报的3个月幅度,所以预报成功的概率是0.01131/12×3 = 0.28%。
但是在5月11日有一个4.9级前震,按照美国地震研究人员统计,“在5天之内、10公里范围内发生更大地震的可能是6%”,其后地震概率衰减如常。根据上述特征斜率0.8994,我们可以得到,4.9级地震之后,“在5天之内、10公里范围内”发生6级以上地震的概率是0.61%。
所以在秦研究员发布他的地震预报的当天,他预报成功的概率是0.28% + 0.61% = 0.89%,也就是说,他的预报有大于99%的可能会失败。
而现在那个前震的前兆效应已过,而且时间窗口也从三个月缩小到只有70天了,所以他那个预报成功的概率,现在只剩下0.22%——他当初的预报有大于99.7%的可能会失败。
注意:本文思路看似简单,但这就是地震概率预测的基本方法(见Sachs et al., 2012, Evaluating the RELM Test Results, International Journal of Geophysics,
http://downloads.hindawi.com/journals/ijgp/2012/543482.pdf)。
https://blog.sciencenet.cn/blog-70036-577389.html
上一篇:
秦四清研究员命中巨震的概率有多大?下一篇:
地震预报的基本定律