答G.Stavrakakis和J.Drakopoulos的“VAN方法:自1981年以来矛盾和误导的结果”

把Stavrakakis和Drakopoulos(1996)(此后称作SD)中关于VAN表述与VAN真实内容作一简单比较,即可证明对方的批评是错误的。典型例子有: (1) 与Drakopoulos(1993)选择并加以讨论的一样,SD所讨论的VAN预报也是VAN预报中的一小部分。Varotsos等(1993b)当时已给予详细回答,并指出Drakopoulos等(1993)的7个例子中有4个没有考虑VAN电报提出的两个可供选择的震中位置,令人吃惊的是没有考虑被强调为最可能的那一个。在该文中SD对VAN预报的真实含意再次做了较大的改动。 (2) SD仅选择了那些VAN预报法具有较大偏差的情形来进行讨论。 (3) SD通过选择地区内的某些预报后在该地区出现的地震活动性而提出怀疑。本文将证明这些预报指的是远离有早期地震活动的另外地区。SD把VAN预报的两个例子说成是对应的余震,我们将证明这种观点是错误的。 (4) SD改变了VAN公布的单个SES和SES电活动性的超前时间。 (5) SD增加(或删除)对VAN观点及预报的批评措词,因此说明对方真正意思是急剧变化的。 (6) 虽然识别SES的标准已经由VAN发表,可SD仍称这种标准依据尚未形成。 (7) SD的示图是错误的,因为他们忽略了M_s=5.0地震(大地震)标准仅发生在预测地区,但也包括远处M_s=4.0～4.2的地震。 (8) SD称“没有一个主震被预报”,但VAN成功地预报了1987年至1989年间的3个处于调查中的最大主震,震中偏差小于50 km。 (9) SD认为如果已经发生了M_s≈4.0的“前震”,再成功地预报M_s=5.0的地震是毫无意义的(不管偏差Δr、ΔM达到多小)。这种观点是不对的。 (10) 虽然VAN已经阐明了只有M_s的预计值大于(或等于)5.0时才提出预报,可SD曲解了实验误差ΔM=±0.7的含意,并要求M_s≥4.3的地震都应被预报出来。这一要求是不对的。此外,我们将进一步证明SD关于“偶然预报成功的概率”计算是错误的,因为它违反了基本原理。     

对Y.Honkura和N.Tanaka的“专门用VAN方法预报希腊地震发生的概率”一文的答复

Honkura 和 Tanaka(1996)以及 Aceves 等(1996)的计算为阐明这一论点的主要问题,即是否 VAN 预报可以归因于偶然,提供了重要的工具。Honkura 和 Tanaka(1996)的计算表明,对半径为120km 的圆形区域和22天的时间窗,在希腊 M_s≥5.0地震(EQ)发生的概率 P 小于0.25,甚至对于11天的时间窗来说更小。对更大震级阈值,即 M_s≥5.5或 M_s≥5.8(及由古登堡-里克特关系式来看),它们的 P 值更小。简单比较这些 P 值与 Mulargia 和 Gasperini(1992;1996a)的表格会立刻显示出 VAN 预报不能归因于偶然。记住,检查后者的表格会得到:1)VAN 成功率为40%～45%(当仅考虑与 M_s≥5.0地震相关时);2)VAN 警报率随地震震级闶值的提高而增加,对 M_(EQ)≥5.5和5.8,分别达到50%和60%。Honkura 和 Tanaka(1996)计算中出现的另一要点是“余震必须仔细处理。”这一点加强了我们在 Varotsos 等(1996a)的原则4和原则5中对 Mulargia 和 Gasperini(1992)对余震的不适当处理(基于泊松分布)的评论:①大大改变了显著性水平值;②预报与地震间的“正向时间相关”转变为“逆向时间相关”。这后一观点分别也被Honkura 和 Tanaka(1996)检查过,他们得出结论:“……考虑这种逆向时间相关的观点……,我们不能在目标区域里找到 M_s≥5.0的地震发生的高概率事件区。”在此答复中,我们还进行一些必要的阐述即在预报方法中,正如所料想的那样,在震级确定中存在试验误差时计算“成功率”和“警报率”有关的内容。     

从巴塘现场地震分析谈持续时间定震级

1989年4月16日北京时间04时34分在四川省巴塘县境内(北纬29°54′,东经99°11′)发生 M_s6.7级地震,继后再度于4月25日和5月3日发生 M_s6.7级和 M_s6.4、M_s6.2级地震,在这次混合型强震群型地震中,M≥1.0级地震超过2万次,截至9月30日止,M_L≥2.0级地震就有3428次,面对如此众多的地震,台站地震分析仍按常规方法进行,不仅人手有限,而且时间上也不允许。大震后震区人民要求地震工作者回答的第一个问题就是还有没有更大的地     

1993年6月份首都圈地区M_L≥2.0地震

1993年6月份首都圈(38.0°—41.0°N,113.0°—120.0°E)范围内,共发生M_L≥2.0地震23次(见表),其中2—2.9级地震19次,有4次3级以上地震,最大震级为M_L3.4。     

1993年8月份首都圈地区M_L≥2.0地震

1993年8月份在首都圈(38.0°-41.0°N,113.0°-120.0°E)范围内,共发生M_L≥2.0地震30次(见下表)。其中2—2.9级的23次,3级以上的7次,无5级以上地震发生。震级最大为8月30日山西浑源东北M_L4.5地震。大同余震区(山西浑源东北)自8月29日起至9月5日出现一组震群活动,最大震级为M_L4.5,其中4级地震4次,3级5次,总计发生1级以上地震131次。     

1989年四川巴塘地震序列特征

1.概况1989年4月16日04时34分在四川省甘孜藏族自治州巴塘县境内发生M_s6.7地震,随后于4月25日10时13分、5月3日13时53分和23时41分分别发生M_s6.7、M_s6.4和M_s6.2地震,到9月30日止,共发生M_L≥1.0地震21000余次,M_L≥2.0地震3428     

不同类型短周期地震仪器对测定M_L的影响

本文用广东台网26个测震台站,10种类型短周期地震仪记录的291次地震,1700多组震级资料,研究了不同类型地震仪测定的近震震级 M_L(K)值与台网值 M_L 之间的系统差。还对 M_L4.0级以上的地震做了单相关回归分析。结果发现①在震中距500km 范围内,用不同类型短周期地震仪测定5级以下中小地震的 M_L 值,一般不出现系统性偏差,一但超越上述检测范围,例如△>600km,M_L>5以后,各类仪器的检测效果将呈现较大差异。②对于 M_L4级以上地震,任何一类地震仪器测定的 M_L(K)值与台网值 M_L之间的单相关性很差,相关系数都比较小。而且台网 M_L 值与 M_s之间的相关系数也很低,因此,不宜相互换算。     

山西省1988年地震简况

据山西省地震台网测定,1988年山西省共发生M_L≥10地震214次~＊.其中1.0至1.9级112次;2.0至2.9级88次;3.0至3.9级14次(见表3).最大震级3.9(M_L),即:1月26日发生在太原盆地文水县开栅镇和4月18日发生在大同盆地朔县城关镇的两次地震.另外9月26日至11月21日和8月23表1 1988年山西境内地震活动月统计表     

全球大地震及中国4级以上地震目录

应广大读者要求,我们定期公布全球7级以上大地震和中国4级以上地震目录。地震参数取自中国地震局分析预报中心的“中国地震速报台网观测报告”和“中国地震月报目录”。其中:发震时间为北京时(BJT);震中位置除给出经纬度外还给出了参考地区名,以供查阅时参考;震源深度33N表示该地震为通常意义上的浅源地震;震级MS为面波震级、ML为近震震级、MW为矩震级(取自美国PDE报告)。表1　世界地震目录(2003年4～6月,M≥7.0)编号发震时间　月 日 时:分:秒　震中位置　φ/(°)　　　　λ/(°)　深度/km震　级　MS　　　MW　地　　区105 26 1…     

滇西地区震级与地震矩标度

本文使用云南省滇西地区,红河断裂带北端1982—1984年由PDR—2数字事件记录仪记到的77个地震的记录图,用美国加州的技术计算了这些地震的地震矩M_0,伍—安地方震震级M_L。对于1(1/2)≤M_L≤3(1/2),得到logM_0与M_L的关系为logM_0=16.97+(1.17±0.05)M_L。与美国中加州的logM_0与K_L的关系一致。将得到的M_L与云南省使用的震级M(持续时间震级),M(振幅震级)比较,发现对这些地震,云南省的震级一般都比M_L大,对于1.7≤M≤5.3,M_L=0.49+(0.67±0.04M),logM_0=17.30+(0.86±0.04),对于2.0≤M≤5.4,M_L=0.24+(0.76±0.05)M,logM_0=17.06+(0.95±0.05)M。     

全球大地震及中国4级以上地震目录(2003年1～3月)

应广大读者要求 ,我们定期公布全球 7级以上大地震和中国 4级以上地震目录。地震参数取自中国地震局分析预报中心的“中国地震速报台网观测报告”和“中国地震月报目录”。其中 :发震时间为北京时(BJT) ;震中位置除给出经纬度外还给出了参考地区名 ,以供查阅时参考 ;震源深度 33N表示该地震为通常意义上的浅源地震 ;震级MS 为面波震级、ML 为近震震级、MW 为矩震级 (取自美国PDE报告 )。表 1　世界地震目录 (2 0 0 3年 1～ 3月 ,M ≥ 7.0 )编号 发震时间　月 日 时 :分 :秒　震中位置　φ/ (°)　　　　λ/ (°)　深度/km震　级　M…     

震源断裂规模与震级的定量关系

研究震源参数(震源断裂规模)与震级的定量关系必须考虑震级的物理含义。本文利用我国海城地震、唐山地震前后的151个小震(1.5≤M_L≤3.9)和国内外58个浅源的中、强震(4.4≤M_s≤8.2)资料,得到了一组稳定的、互相一致的震源断裂参数与震级M_L和M_S的函数关系。由此讨论了常用震级M_L、M_S与m_b之间的关系。 新的震源参数与震级关系表明,震源断裂规模如断裂长度、面积或错距等与面波震级M_S和近震震级M_L不存在一一对应的比例关系     

最大熵原理在地震重复关系中的应用

本文应用最大熵原理确定震级概率密度分布函数P(M)=βe~(-βM)/(e~(-βM_o)-e~(-βM_u)) M_o≤M_(?)(1)这里 M_o 是非负震级下限,M_u 是震级上限,参数β可由1/β+(M_oe~(-βM_o)-M_ue~(-βM_u)/(e~(-βM_o)-e~(-βM_u)=(?)(2)确定。用数值方法求解方程(2)并且把β的数值代入方程(1),得到震级 M 的最小无偏概率密度分布。这种分布是与 M_o、M_u 以及平均震级(?)的可利用信息相一致的。震级大于或等于 M 的地震数N(M)=T(e~(-βM)-e~(-βM_u))/(e~(-βM_o)-e~(-βM_u))(3)这里 T 是 M≥M_o 的地震总数。当 M_o→0并且 M_u→∞时,(3)式简化为N(M)=Te~(-βM)(4)方程(4)精确地等价于古登堡—李克特关系log N=α-bM(5)只需令a=log T,并且b=β/ln 10。本文把基于最大熵原理得到的重复关系与全球 M_s≥7、全国 M_s≥6以及鄂尔多斯地块周边 M_s≥5级的地震观测资料作了对比。结果表明,在较大的震级跨度范围内,应用最大熵原理得到的重复关系与实际资料更为吻合。最后求得了各级地震的复发周期。上述结果显然可应用于地震危险性分析中。     

对P.W.Burton的“地震赌博”一文的答复

Burton(1996)的文章除了它的奇怪的标题外,还有很多不精确或错误的地方。Burton 明显地滥用了统计学,因为他处理一个参数(例如,震级),并不注意基本点,例如:涉及到一个区域中的一个主震的预报震级值,这样的事件是非常少的,不能被等同处理为有同样震级的预报,但涉及的却是另一区域的余震。利用统计的简单规律,我们证明了 Burton 所用的统计处理显然是错误的,甚至可以“否定”一个理想的预报方法。此外,我们证明了 Burton 对 VAN 震级预报的批评是基于对实验误差真实含义以及古登堡-里克特关系明显错误解释之上的;Burton 的论点与物理预期的结果相矛盾。对地震预报合适的空间解也进行了讨论,对于5.5～6.0之间的震级,VAN 所得到的空间解大约是50km。如果我们考虑到源的尺度,这样一个解从实际的观点看是合理的。对于较大的震级,例如:约7.5级,约150km 的空间解应当是可以接受的。因此 Burton 对 VAN 增加空间解的论点是不合理的,因为他们仅对较小事件(～5.0个单位)有一些根据,不过这些没有多少实际的意义。最后但不是最不重要的,Burton 利用了不寻常的措辞,如:“有效预报的事后论证”、“不清楚的回顾性滤波”、“实际的东西太少而宣传的太大”,但没有提到事实,例如:在讨论的3年中(即1987～1989),VAN 仅仅发出了一个公开的警报,这一警报后面有一个这一时期最强烈的破坏性活动。     

1992年5月份首都圈地区M_L≥2.0地震

1992年5月份在首都圈(38.5°—41.0°N,113.0°—120.0°E)范围内,共发生 M_L≥2.0地震17次,其中2至2.9级地震15次,3至3.9级地震2次,震级最大的为5月15日滦县北 M_L3.7地震,详见下表。     

华北地震活动对辽西的影响

1.1966年至1987年华北M_s≥5.5级地震取自中国地震目录和中国地震台网观测报告以及1981—1987年辽宁省地震局公布的国内5.0级以上地震简录;辽西M_L≥2.0级地震取自辽宁省地震局1970—1987年地震目录和中国东部地震目录。 2.华北和辽西的划分及华北M_s≥5.5级地震分布 文中所说的华北范围大致为北纬31°至43°,东经105°至123°(不包括海洋地区);辽西为北纬40°至42°30′,东经118°至120°并用斜线表示,详见图1。     

全球大地震及中国4级以上地震目录

应广大读者要求,我们定期公布全球7级以上大地震和中国4级以上地震目录。地震参数取自中国地震局分析预报中心的“中国地震速报台网观测报告”和“中国地震月报目录”。其中:发震时间为北京时(BJT);震中位置除给出经纬度外还给出了参考地区名,以供查阅时参考;震源深度33N表示该地震为通常意义上的浅源地震;震级MS为面波震级、ML为近震震级、MW为矩震级(取自美国PDE报告)。表1　世界地震目录(2003年7～9月,M≥7.0)编号发震时间　月 日 时:分:秒　震中位置　φ/(°)　　　　λ/(°)　深度/km震　级　MS　　MW　地　　区107 16 04:…     

全球大地震及中国4级以上地震目录

应广大读者要求,我们定期公布全球7级以上大地震和中国4级以上地震目录。地震参数取自中国地震局分析预报中心的“中国地震速报台网观测报告”和“中国地震月报目录”。其中:发震时间为北京时(BJT);震中位置除给出经纬度外还给出了参考地区名,以供查阅时参考;震源深度33N表示该地震为通常意义上的浅源地震;震级MS为面波震级、ML为近震震级、MW为矩震级(取自美国PDE报告)。表1　世界地震目录(2004年1～3月,M≥7.0)编号发震时间　月 日 时:分:秒　震中位置　φ/(°)　　　　λ/(°)　深度/km震　级　MS　　MW　地　　区102 07 10:…     

2020年河北唐山5.1级地震前震源参数变化

利用谱分析方法对河北省测震台网2016年1月—2020年7月唐山老震区(117.5°～119.5°E,39°～41°N)85次M_L≥2.5地震进行震源参数分析,发现该时段唐山老震区应力降范围在0.1～10.3MPa,主要集中在0.1～1.0MPa,其中应力降小于2MPa的地震占81%。当2.5≤M_L≤3.5时,应力降、拐角频率与震级之间无明显的对应关系;ML≥4.0地震较少,但应力降有随震级增大而增大的趋势,拐角频率有随震级增大而减小的趋势。2019年12月5日河北丰南M4.5地震和2020年7月12日唐山M5.1地震之间发生了 7次M_L≥2.5地震,震级范围为2.5≤M_L≤2.8,提取该震级范围的震源参数进行时空分析,发现唐山M5.1地震前有2次明显的高应力降事件,2次高应力降事件呈现出应力降数值高、拐角频率高、震源尺度小的特点,表明单位体积所释放能量较多。2次高应力降事件均发生在唐山M5.1地震震中17km范围内,可能是唐山M5.1地震前震源区构造应力较高的一个标志。     

1992年7月份首都圈地区M_L≥2.0地震

1992年7月份在首都圈地区(38.5°—41.0°N,113.0°—120.0°E)范围内,共发生 M_L≥2.0地震9次,其中小于等于2.5级地震7次,震级最大为7月22日宁河 M_L4.9地震,活动水平相对偏低(见下表)。