|
陈一文顾问按:沈宗丕老师的这篇文章发表在 2002年10 月在中国地球物理学会天灾预测专业委员会《“磁暴二倍法”学术研讨会论文集》中,较详细地说明了中国人原始创新的“磁暴二倍法”到“磁偏角异常二倍法”以及“磁暴月相二倍法”预测地震的发展由来及其预测效果。
文摘∶利用许多地震前往往发生的两次磁暴之间的间隔时间,延长一倍,可以计算许多地震的发震日期。本文介绍张铁铮
1969年底开发出利用这种自然规律的“磁暴二倍法”,在中国首次应用于发震日期的预测。本文作者在此基础上开发出“磁偏角异常二倍法”,从
1970年9
月开始预测发震日期。 90年代,作者在预测环太平洋地震带上8
级左右大地震时,发现二个异常日期之间的天数,如果符合
29.6天的倍数,与8
级左右大地震发震日期有更好的对应关系。
29.6天近似月球的望、朔周期。作者开发的“磁暴月相二倍法”预测发震日期由此产生。
1991年12
月1日至1994年11
月30日内,全球共发生
Ms>=7.5的大震(主震),共12
次。作者预测了14次(其中1
次重复,应按13次统计)。对应MS>=7.5
大震有8次,虚报5
次(其中对应上7.4>=Ms>=7.0
的大震3次)
,漏4次。
1998
年5月
1日~2001年1月31日期间全球
Ms³7
.5大地震共发生
16次,发生于“磁暴月相二倍法
”得出的计算发震日期〔在±5天的范围内〕有13次
,无对应的有3次。在15次计算发震日期中有
11次对应全球Ms³7.5的大地震,2次对应7.4
£
Ms³7.0的大震。K指数大的“起倍磁暴”日与Ms
³
7.8
的巨震的发生有较好的相关性。这表明发生在月相的磁暴与
Ms³7.5
大地震的相关关系较好。
1.
张铁铮及其“磁暴二倍法”预测地震 1966年3
月邢台强烈地震发生后,原先从事石油物探工作的张铁铮工程师在参加地震预测的研究工作中,发现在地震之前地磁场往往出现大幅度的异常变化,一般认为是磁暴引起的。这样提出了一个问题∶磁暴同地震有没有联系? 通过以北京白家疃地磁台的磁暴资料与同期的地震资料对照分析,张铁铮发现∶“在168个磁暴中,与地震同时发生的有90
个,占总数的54%;震前一天出现的21
个,占12%;震后一天出现的16
个,占10%;也就是说有127
个磁暴,占总数76%,是在地震发生的同时或前后一天内出现的。没有发生地震的磁暴41
个,占24%。……上面的统计分析表明,磁暴的出现与地震发生有一定的联系。”〔耿庆国,1991〕 对此进行总结时,张铁铮指出〔张铁铮,1975〕∶“激发地震发生的外因条件并不只磁暴一种。因此,磁暴与地震之间的对应关系是有条件的,不是绝对的。” |
![]() |
对已经发生过的有对应关系的磁暴和地震,进行对比分析,张铁铮进一步发现∶“两个先后发生磁暴的时间
T1和T2
,其间隔的天数延长一倍,加上T1的日期,与它们所对应的地震发生日期MT
相近。用公式表示即∶MT = T1 + 2(T2 – T1)
”〔图1和表1
〕
图1 “起倍磁暴日”、“被倍磁暴日”和预测的地震发震日之间的关系
表1
张铁铮对60年代中国大地震与磁暴对应关系进行分析的数据
|
序号 |
起倍磁暴 日期 |
被倍磁暴 日期 |
推算发震 日期 |
实际发震 日期 |
震级 |
天数差 |
地震地点 |
|
1 2 3 4 5 6 7 8 |
1961.04.15. 1965.07.19. 1964.02.06. 1964.05.14. 1967.01.09. 1967.02.16. 1968.09.07. 1967.05.03. |
1961.10.01. 1965.09.16. 1965.02.06. 1965.04.18. 1967.02.16. 1967.05.25. 1969.02.11. 1968.06.13. |
1962.03.19. 1965.11.14. 1966.02.07. 1966.03.23. 1967.03.27. 1969.08.31. 1969.07.18. 1969.07.25. |
1962.03.19. 1965.11.13. 1966.02.05. 1966.03.22. 1966.03.27. 1967.08.30. 1969.07.18. 1969.07.26. |
6.1 6.6 6.5 7.2 6.3 6.8 7.4 6.4 |
0 1 2 1 0 1 0 -1 |
广东河源 乌鲁木齐 云南东川 河北宁晋 河北河间四川甘孜 渤海 广东阳江 |
接着,张铁铮“用这种方法进行未来地震的预报试验,选择了
1968年6
月11日和1969年3
月24日两个较大磁暴,推算1970年1
月4日将有大地震发生。结果,一月五日云南省通海地区发生了Ms7.7
强烈地震,推算的日期比实际发震日期早一天。这次没有报发震地区,不能起到预防以后。吸取了经验和教训以后,他们增加了多台对比同一天的同一个磁暴幅度的变化来确定发震地区,选择1968年6
月11日和1969年4
月13日两个磁暴,推算1970年2
月13日,华北地区要有
Ms5.0级左右的地震。
2月12
日渤海湾发生了Ms4.8
级地震,比推算的早了一天,预报的地区范围大了。尽管方法还不完善,但是初步验证的结果说明磁暴是可以用来预报地震的。”参看表
1。张铁铮将这种地震预测方法称为“磁暴二倍法”〔张铁铮,1975〕。
张铁铮是中国第一位总结出利用前兆资料预测大震方法取得成功的人。
表2 张铁铮采用“磁暴二倍法”在中国首次预测的两次地震及其对应的磁暴
|
序号 |
起倍磁暴 日期 |
被倍磁暴 日期 |
推算发震 日期 |
实际发震 日期 |
震级 |
天数差 |
地震地点 |
|
1 2 |
1968.06.11. 1968.06.11. |
1969.03.24. 1969.04.13. |
1970.01.04. 1970.02.13. |
1970.01.05. 1970.02.12. |
7.7 4.8 |
-1 -1 |
云南通海 渤海湾 |
在震前预测过
1970年1
月5日云南通海
Ms7.7级大地震的发震时间后,张铁铮受到过国务院周恩来总理的接见。
2.
“磁偏角异常二倍法”预测地震的由来
1970年,张铁铮根据“磁暴二倍法”较准确地预测了
1970年8
月10日山东曲阜地区发生的
5.0级地震。他用
1970年7
月9日和1970年7
月25日两个磁暴日的地磁垂直强度异常相隔时间的二倍,预测这次地震的发震时间。这个经验启发了我们。
1970年8
月,趁张铁铮到北京白家疃地磁台抄录磁暴目录的机会,作者向他了解有关“磁暴二倍法”的方法。
张铁铮采用地磁垂直强度这个要素搞地震预测,作者则采用磁偏角要素来搞地震预测。作者考虑到磁偏角有几个优点∶一、不受温度影响,磁针始终指南指北;二、只要仪器水平,室内干燥,仪器不会有零漂现象;二、容易被广大群众所使用,而且能够推广。
为了突出震磁信息,作者使用二台〔北京台与河北省红山台〕磁偏角幅度相减的办法来消除外空磁场的影响。与张铁铮的“磁暴二倍法”相同,计算方法也是用二倍的关系来预测地震。1970年9
月6日我在内部进行第一次试报∶1970年9
月9日±1
天,在邢台
150公里范围内,可能发生一次
3级以上地震。果然于
1970年9
月8日凌晨三时三十一分在牛家桥发生了一次
ML=4.0级的地震,离邢台
50公里。作者觉得这个方法不错,称这个方法为“磁偏角异常二倍法”。
第一次内部试报获得成功后,作者于
1970年9
月16日正式向中央地震办公室〔国家地震局前身〕预测∶1970年9
月29日±1
天,在邢台
150公里范围内可能发生一次
4级左右地震。
结果于
1970年9
月29日凌晨五时三十二分在河北省磁县发生了一次
ML=4.5级的地震,离邢台
70公里左右。后来就参加了中央地办组织的多次地震会商会,而且在会商会每次都作预测。有时能报准,有时也报虚。
2.1
“磁偏角异常二倍法”较好预报1972年1
月25日台湾省
8级大地震
用
1970年3
月9日北京台〔
D
D〕减去上海佘山台〔DD
〕的DRD=13.8¢,与
1971年2
月16日北京台〔
D
D〕减去上海佘山台〔DD
〕的DRD=4.4¢,相隔344
天二倍后得
1972年1
月26日。据此,作者第一次于1971年9
月13日向国家地震有关部门填好地震短临预报卡片预报∶1972年1
月26日±
1天,在地球上可能发生一次8
级左右大地震。
第二次于
1971年10
月28日向国家地震局有关部门填好地震短临预报卡片预报∶1972年1
月26日±
1天,在我们台湾省或日本,可能发生一次8级左右大震。
结果于
1972年1
月25日,在我国台湾省东部海域发生了台湾地区
50多年来最大的一次8
级大地震。地震发生后得到国家地震局领导的肯定。
1972年第三期《地震战线》刊登了一篇报导,题为“利用磁偏角二倍法较好地预报了台湾8级地震”。
2.2
“磁偏角异常二倍法”较好预报1973年9
月29日日本海
8级深震
用
1970年3
月9日北京台〔
D
D〕减去上海佘山台〔DD
〕的DRD=13.8¢,与
1971年12
月18日北京台〔
D
D〕减去上海佘山台〔DD
〕的DRD=6.5¢,相隔649
天二倍后得
1932年9
月27日。据此,作者第一次于
1972年3
月在北京召开的第二次全国地震工作会议上填好地震短临预报卡片预报∶1973年9
月27日±
1天,在日本到阿留申群岛之间可能发生一次8级左右大地震。
第二次于1972
年10月在山西省临汾召开的地震中期预报趋势会商上重申上述预报意见,并提交一篇题为“用磁偏角二倍法作中期预报的试探”一文。结果于1973年9
月29日,在日本海发生了一次
8级深震。
国家地震局《地震战线》编辑部编印的1973
年《地震》技术汇编2上发生了这篇文章。
“我们在预报环太平洋地震带上7.5
级左右的地震时,是用北京台溢佘山台的幅度,经过年变更正后,若
D
RD=3.5¢,则认为是异常。二个异常相隔的天数必须要符合一定的周期方可预报地震。其原则是起倍日期
D
RD=10.0¢上被倍日期DRD=3.5
¢
进行二倍〔表3〕。
表
3 起倍日期
1972年8
月5日北京台
-佘山台DRD=10.0
¢
|
被倍日期 |
D
RD |
相隔天数 |
T=29.6
天 |
计算发震日期 |
阴历 |
实际发生地震情况 |
|
1972.09.14 |
4¢.0 |
40
天 |
X
X |
1972.10.24 |
|
X
X |
|
1972.
11.01 |
4.0 |
88
天 |
(3) 88.8 |
1973.01.28 |
廿五 |
1973.01.31.墨西哥7.9级 |
|
1972.
12.16 |
4.0 |
133
天 |
X
X |
1973.
04.28 |
|
X
X |
|
1973.07.26 |
3.6 |
355
天 |
(12) 355.2 |
1974.07.16 |
廿七 |
1974.11.09.巴拿马7.5级
|
|
1973.
09.23 |
7.0 |
414
天 |
(14) 414.4 |
1974.
11.11 |
廿八 |
1974.07.13.秘鲁7.5级
* |
|
1973.
10.18 |
3.5 |
439
天 |
X
X |
1974.
12.31 |
|
X
X |
|
1973.
10.29 |
5.9 |
450
天 |
X
X |
1975.01.22 |
|
X
X |
|
1973.
11.22 |
3.6 |
474
天 |
(16) 473.6 |
1975.
03.11 |
廿九 |
1975.03.13
智利中部7
级* |
|
1973.
12.21 |
3.9 |
503
天 |
(17) 503.2 |
1975.
05.08 |
廿六 |
1975.05.10. |
|
1974.01.27 |
3.5 |
540
天 |
X
X |
1975.
07.21 |
|
X
X |
|
1974.
02.12 |
4.2 |
556
天 |
X
X |
1975.
08.22 |
|
X
X |
|
1974.
03.21 |
3.8 |
593
天 |
(20) 592.0 |
1975.
11.04 |
初一 |
1975.10.31菲律宾7.1级 |
|
1974.
04.19 |
3.5 |
622
天 |
(21) 621.6 |
1976.01.01 |
初一 |
1976.01.01
克
马德克7.1级 |
|
1974.
07.06 |
4.5 |
700
天 |
X
X |
1976.
06.05 |
|
X
X |
|
1974.
08.29 |
4.5 |
754
天 |
X
X |
1976.
09.21 |
|
X
X |
|
1974.
09.16 |
4.1 |
772
天 |
(26) 769.6 |
1976.
10.27 |
初五 |
1976.10.29伊里安7.2级 |
|
1974.
09.20 |
4.3 |
776
天 |
X
X |
1976.
11.04 |
|
X
X |
|
1974.
10.15 |
3.5 |
799
天 |
(27) 799.2 |
1976.
12.20 |
卅 |
1976.12.21加拿大7.0级 |
*
震前曾向国家地震局作过短临预报。
3.
“磁暴月相二倍法”预测地震的由来
“磁暴月相二倍法”预测地震是在上述“磁偏角异常二倍法”的基础上发展起来的〔沈宗丕,
1977〕。“磁偏角异常二倍法”是通过南北相距较远的两个地磁台同一天的磁偏角的幅度值相减,并经过适当的纬度“校正”,选出突出的地区性异常,来估计地震发生的大致地点;根据两次磁偏角异常出现的日期中间所包括的天数,从第二次异常日期算起,往后推同样的天数,这就是预测发震的时间;震级的大小是根据异常的大小来估计的。
在使用北京台与佘山台的磁偏角数据作“磁偏角异常二倍法”预测地震的过程中,发现可以预测环太平洋地震带上
8级左右的大地震,但是还存在着一定的虚报。通过不断的总结发现,二个异常日期之间的天数,必须要符合
29.6天的倍数,方可进行预测,才能对应上8
级左右的大地震,否则会带来虚报,(指小于
7级地震)。29.6
天正好近似月球的望、朔周期,从中又发现所选用的异常日期大多是太阳上发生的大耀斑和质子事件,所引起的磁暴日。因此,“磁暴月相二倍法”预测地震就从这里产生了〔沈宗丕,
1977〕。
4.
“磁暴月相二倍法”预测地震的方法介绍
4.1
“磁暴月相二倍法”要区分两种性质的磁暴∶
“起倍磁暴”(MS1)和“被倍磁暴”(MS2
)。预测地震时间的计算是求出“起倍磁暴日”与“被倍磁暴日”的时间间隔(
D),即D=MS2-MS1
,以天为单位。在“被倍磁暴日”的日期上加上
D值,即为预测“发震日期”(TC
),即TC=MS2+D,
误差一般为±3天或±
6天。
4.2
“起倍磁暴日”与“被倍磁暴日”的选取∶
大的磁暴日,大多是太阳上发生的大耀斑或质子事件所引起的,在选取“起倍磁暴日”的时候必须要选
K指数大的磁暴。我们国家的地磁台,以采用三个小时时段内(国际时)水平强度(H)最大幅度与K
指数之间的关系如下∶
|
R
(H幅度)= |
0 |
3 |
6 |
12 |
24 |
40 |
70 |
120 |
200 |
300
以上(NT) |
|
K
指数= |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
是以三小时的时段来量算的,分别为00h~03
h 为第
1时段,03 h
~06 h为第
2时段,……21 h
~24 h 为第
8时段。当K=5
时为中常磁暴(m);K=6
或7时为中烈磁暴(ms
);K=8或9
时为强烈磁暴(s)。在选取“起倍磁暴日”时,必须K
³
7。在选取“被倍磁暴日”时,必须K
³
6,而且都应该在月相的日期中选取(上弦日为初七~初九;望日为十四
~十六;下弦日为二十~
廿三;朔日为廿九(三十)~初一)。但“被倍磁暴日”的K
指数一般不能超过“起倍磁暴日”的
K指数〔个别情况例外〕。
5.
“磁暴月相二倍法”对全球
Ms³7.5大地震的预测效果(1991
年12月~1994
年11月)
5.1
资料的选取;
(1)地震目录∶表4
中列出1991年12
月~1994年11
月全球Ms³7.5
大地震目录。地震的参数取自中国地震局地球物理研究所陈培善在《地震学报》上发表的“全球
M³6
的地震目录”。
表4
、1991年12
月1日~1994
年11
月30
日
全球.Ms≥7.5
大震(主震)
目录
|
序号 |
发震日期 |
发震地点 |
震级 |
序号 |
发震日期 |
发震地点 |
震级 |
|
1 |
1991.12.27 |
南桑德韦奇群岛 |
7.5 |
7 |
1992.12.12 |
印度尼西亚 |
7.7 |
|
2 |
1992.05.26 |
牙买加地区 |
7.7 |
8 |
1993.07.12 |
日本北海道 |
7.6 |
|
3 |
1992.06.28 |
美国加州南部 |
8.0 |
9 |
1993.08.08 |
马里亚纳群岛 |
7.6 |
|
4 |
1992.08.19 |
吉尔吉斯斯坦 |
7.7 |
10 |
1993.09.11 |
墨西哥哈帕斯 |
7.6 |
|
5 |
1992.09.02 |
尼加拉瓜近海 |
7.5 |
11 |
1994.06.09 |
秘鲁 |
7.6 |
|
6 |
1992.10.18 |
哥伦比亚西岸 |
7.7 |
12 |
1994.10.04 |
日本北海道 |
7.9 |
(2) 磁暴日目录∶表5
中列出
1991年8
月20日~1993
年4
月30
日
在月相上弦日和下弦日K³6的磁暴日〔取中国地震局地球物理研究所编印的《磁暴报告》中的目录〕。
表
5、1991
年8月20
日~1993
年4
月30
日, K
³
6,上弦与下弦磁暴日目录
|
序号 |
磁暴日期 |
农历 |
K |
月相 |
序号 |
磁暴日期 |
农历 |
K |
月相 |
|
1 |
1991.08.30 |
七月廿一 |
7 |
下弦 |
7 |
1992.08.20 |
七月廿二 |
6 |
下弦 |
|
2 |
1991.10.28 |
九月廿一 |
9 |
下弦 |
8 |
1992.09.03 |
八月初七 |
6 |
上弦 |
|
3 |
1992.02.25 |
一月廿二 |
6 |
下弦 |
9 |
1992.09.17 |
八月廿一 |
6 |
下弦 |
|
4 |
1991.05.10 |
四月初八 |
8 |
上弦 |
10 |
1993.02.28 |
二月初八 |
6 |
上弦 |
|
5 |
1992.05.22 |
四月二十 |
6 |
下弦 |
11 |
1993.04.13 |
三月廿二 |
6 |
下弦 |
|
6 |
1992.08.05 |
七月初七 |
6 |
上弦 |
|
|
|
|
|
5.2
计算发震时间与大震对应情况∶
表6、表7
和表8表明作者利用
1991年8
月30日和1991年10
月28日两个“起倍磁暴”日与
1991年~1993
年一系列“被倍磁暴”日预测的地震发震日期和大震实际发震日的对应情况。
表
6、起倍磁暴日∶
1991年8
月30日,农历七月廿一,K=7
,与下弦磁暴日二倍
|
被倍磁暴日 |
相隔 天数 |
计算 发震日期 |
实际发生地震 |
误差 (天) |
||||
|
公历 |
农历 |
K |
公历 |
地点 |
震级 (Ms) |
|||
|
1991.10.28 |
九月廿一 |
9 |
59 |
1991.12.26 |
1991.12.27 |
南桑德韦奇 群岛 |
7.5 |
+1 |
|
1992.02.25 |
一月廿二 |
6 |
179 |
1992.08.22 |
1992.08.19 |
吉尔吉斯斯坦 |
7.7 |
-3 |
|
1992.05.22 |
六月二十 |
6 |
266 |
1993.02.12 |
1993.02.07 |
日本本州 西海岸 |
6.9 |
-5 |
|
1992.08.20 |
七月廿二 |
6 |
356 |
1993.08.11 |
1993.08.08 |
马里亚纳 群岛 |
7.6 |
-3 |
|
1992.09.17 |
八月廿一 |
6 |
384 |
1993.10.06 |
1993.10.02 |
中国 |
6.6 |
-4 |
预测
5次,
对应Ms³7.5
的大震有3次。
表
7、起倍磁暴日∶
1991年8
月30日,农历七月廿一,K=7
,与上弦磁暴日二倍
|
被倍磁暴日 |
相隔天数 |
计算 发震日期 |
实际发生地震 |
误差 (天) |
||||
|
公历 |
农历 |
K |
公历 |
地点 |
震级(Ms) |
|||
|
1992.05.10 |
六月初八 |
8 |
254 |
1993.01.19. |
1993.01.15 |
北海道 |
7.2* |
-4 |
|
1992.08.05 |
七月初七 |
6 |
341 |
1993.07.12. |
1993.07.12 |
北海道 |
7.8 |
0 |
|
1992.09.03 |
八月初七 |
6 |
370 |
1993.09.08. |
1993.09.11 |
墨西哥哈帕斯 |
7.6 |
+3 |
|
1993.02.28 |
二月初八 |
6 |
548 |
1994.08.30. |
1994.09.01 |
美国西部海中 |
7.0 |
+2 |
*
日本国内定此震级为7.8级。
预测
4次,对应Ms
³
7.5的大震有2
次
表
8、起倍磁暴日∶
1991年10
月28日,农历九月廿一,K=9,与下弦磁暴日二倍
|
被倍磁暴日 |
相隔天数 |
计算 发震日期 |
实际发生地震 |
误差 (天) |
||||
|
公历 |
农历 |
K |
公历 |
地点 |
震级 (Ms) |
|||
|
1992.02.25 |
一月廿二 |
6 |
120 |
1992.06.24 |
1992.06.28 |
美国加州南部 |
8.0 |
+4 |
|
1992.05.22 |
四月二十 |
6 |
207 |
1992.12.15 |
1992.12.12 |
印度尼西亚 |
7.7 |
-3 |
|
1992.08.20 |
七月廿二 |
6 |
297 |
1993.06.13 |
1993.06.08 |
堪察加半岛 |
7.3 |
-5 |
|
1992.09.17 |
八月廿一 |
6 |
325 |
1993.08.08 |
1993.08.08 |
马里亚纳群岛 |
7.6 |
0 |
|
1993.04.13 |
三月廿二 |
6 |
533 |
1994.09.28 |
1994.10.04 |
日本北海道 |
7.9 |
+6 |
预测
5次,对应Ms
³
7.5的大震有4
次,其中7.7级以上强震有3
次。
5.3
预测效果∶
在
1991年12
月1日~1994
年11月30
日内全球共发生Ms³7.5
的大震〔主震〕共12次。预测了14
次〔其中有1次重复)
,按13次统计,对应Ms
³
7.5的大震有8
次,虚报5次〔其中对应上7.4
£
Ms³7.0
的大震有3次)
,漏4次〔沈宗丕和徐道一,1996
;徐道一等,1994a;徐道一等,1994b
〕。
6.
“磁暴月相二倍法”对全球
Ms³7.5大地震的预测效果
(1998年5
月~2001年1
月)
6.1
资料的选取∶
(1)地震目录∶表9
中列出1998年5
月~2001年1
月,全球MS³7.5
大地震目录。地震的参数,取自中国地震局地球物理研究所陈培善在《地震学报》上发表的“全球
M³6
的地震目录”。
表
9、1998
年5月1
日~2001
年1
月31
日
全球MS³7.5大震〔主震〕目录
|
序号 |
发震日期 |
发震地点 |
震级 |
序号 |
发震日期 |
发震地点 |
震级 |
|
1 |
1998.05.04 |
台湾东南以远地区 |
7.6 |
9 |
1999.11.13 |
土耳其境内 |
7.6 |
|
2 |
1998.08.05 |
厄瓜多尔海岸 |
7.5 |
10 |
2000.06.05 |
苏门答腊 西南 |
7.9 |
|
3 |
1998.11.29 |
斯兰海 |
7.7 |
11 |
2000.06.18 |
南印度洋 |
8.0 |
|
4 |
1999.08.17 |
土耳其西部 |
8.0 |
12 |
2000.11.16 |
新不列颠 地区 |
7.7 |
|
5 |
1999.08.20 |
哥斯达利加海岸 |
7.6 |
13 |
2000.11.16 |
新爱尔兰 地区 |
7.8 |
|
6 |
1999.09.21 |
中国台湾南投 |
7.7* |
14 |
2001.01.11 |
科迪亚克岛地区 |
7.5 |
|
7 |
1999.10.01 |
墨西哥瓦哈卡 |
7.9 |
15 |
2001.01.14 |
中美洲海岸 |
8.4 |
|
8 |
1999.10.16 |
加利福尼亚以南 |
7.6 |
16 |
2001.01.26 |
印度 |
7.9 |
*
震级为Msz
(2) 磁暴日目录∶:
表10中列出
1997年5
月15日~2000
年9
月18
日
在月相中K³6的磁暴日,取自中国地震局地球物理研究所编辑的《磁暴报告》中的目录。
表11. “起倍磁暴日”、“被倍磁暴日”和全球Ms
³
7.5大地震的对应情况
|
起倍磁暴日 序号 |
被倍磁暴日期 |
相隔 时间天数 |
计算发震日期 |
实际发生地震对应情况 |
误差 天数 |
||||
|
公历 |
农历 |
K |
|
公历 |
公历 |
地点 |
Ms |
|
|
|
1 |
1997.11.07 |
十月初八 |
6 |
176 |
1998.05.02 |
1998.05.04 |
台湾东南以远地区 |
7.6 |
+2
Δ▼ |
|
|
1998.07.31 |
六月初九 |
6 |
442 |
1999.10.16 |
1999.10.16 |
加利福尼亚南部 |
7.6 |
0
▼ |
|
|
1998.12.25 |
十一月初七 |
6 |
589 |
2000.08.05 |
2000.08.05 |
萨哈林岛 |
7.4 |
0
▼ |
|
2 |
1998.07.31 |
六月初九 |
6 |
88 |
1998.10.27 |
1998.10.29 |
马鲁古海峡 |
6.3 |
+2 |
|
|
1998.12.25 |
十一月初七 |
6 |
235 |
1999.08.17 |
1999.08.17 |
土耳其西部 |
8.0 |
0
▼ |
|
|
|
|
|
|
|
1999.08.20 |
哥斯达黎加海岸 |
7.6 |
+3
▼ |
|
3 |
1999.03.01* |
一月十四 |
6 |
207 |
1999.09.24 |
1999.09.21 |
中国台湾南投 |
7.7 |
-3 |
|
|
1999.10.22 |
九月十四 |
7 |
442 |
2000.01.06 |
2000.01.11 |
科迪亚克岛 |
7.5 |
-5
Δ |
|
4 |
1998.11.09 |
九月廿一 |
6 |
21 |
1998.11.30 |
1998.11.29 |
斯兰海 |
7.7 |
-1 |
|
5 |
2000.05.24 |
四月廿一 |
7 |
244 |
2001.01.23 |
2001.01.26 |
印度 |
7.9 |
+3
Δ▼ |
|
6 |
2000.02.12 |
一月初八 |
7 |
113 |
2000.06.04 |
2000.06.05 |
苏门答腊西南 |
7.9 |
+1
Δ▼ |
|
|
2000.06.08 |
五月初七 |
6 |
230 |
2001.01.24 |
2001.01.26 |
印度 |
7.9 |
+2
Δ▼ |
|
7 |
2000.06.08 |
五月初七 |
6 |
119 |
2000.10.05 |
2000.10.06 |
日本本州南部 |
7.2 |
+1
Δ▼ |
|
8 |
2000.06.08 |
五月初七 |
6 |
15 |
2000.06.23 |
2000.06.18 |
南印度洋 |
8.0 |
-5 |
|
|
2000.06.23 |
五月廿二 |
6 |
30 |
2000.07.23 |
2000.07.21 |
哥斯达黎加海岸 |
6.6 |
-2 |
|
|
2000.09.18 |
八月廿一 |
8 |
117 |
2001.01.13 |
2001.01.14 |
中美洲海岸 |
8.4 |
+1 |
|
9 |
2000.09.18 |
八月廿一 |
8 |
64 |
2000.11.21 |
2000.11.16 |
新不列颠地区 |
7.7 |
-5 |
|
|
|
|
|
|
|
2000.11.16 |
新尔兰地区 |
7.8 |
-5 |
*未选入《磁暴报告》内。
Δ在震前填表“地震短临测卡片”寄给中国地震局有关部门。
▼在震前曾向中国地球物理学会“天灾预测专业委员会”提交发震时间与震级的书面预测报告。
表
10 在月相中K
≥6的磁暴日与9
个起倍磁暴日
|
磁暴日期 |
月相 |
K |
序号 |
磁暴日期 |
月相 |
K |
序号 |
||
|
公历 |
农历 |
公历 |
农历 |
||||||
|
1997.05.15 |
四月初九 |
上弦 |
7 |
1 |
1998.12.25 |
十一月初七 |
上弦 |
6 |
|
|
1997.10.02 |
九月初一 |
朔 |
6 |
|
1999.03.01* |
一月十四 |
望 |
6 |
|
|
1997.11.07 |
十月初八 |
上弦 |
6 |
|
1999.09.23 |
八月十四 |
望 |
8 |
5 |
|
1997.11.22 |
十月廿三 |
下弦 |
6 |
|
1999.10.22 |
九月十四 |
望 |
7 |
6 |
|
1998.05.04 |
四月初九 |
上弦 |
8 |
2 |
2000.02.12 |
一月初八 |
上弦 |
7 |
7 |
|
1998.07.16 |
五月廿三 |
下弦 |
6 |
|
2000.05.24 |
四月廿一 |
下弦 |
8 |
8 |
|
1998.07.31 |
六月初九 |
上弦 |
6 |
|
2000.06.08 |
五月初七 |
上弦 |
6 |
|
|
1998.08.06 |
六月十五 |
望 |
7 |
3 |
2000.06.23 |
五月廿二 |
下弦 |
6 |
|
|
1998.10.19 |
八月廿九 |
朔 |
7 |
4 |
2000.07.16 |
六月十五 |
望 |
9 |
9 |
|
1998.11.09 |
九月廿一 |
下弦 |
6 |
|
2000.09.18 |
八月廿一 |
下弦 |
8 |
|
“序号”∶起倍磁暴序号
*
未被选入《磁暴报告》内。
6.2
计算发震时间与大地震对应情况:
表11
〔前边一页〕表明作者利用1991年5
月4日和1991年11
月13日之间九个“起倍磁暴”日与1997年11
月7日~2000
年9
月18
日
之间一系列“被倍磁暴”日预测的地震发震日期和大震实际发震日的对应情况。
6.3预测效果
对
1998年5月1日~2001年1月31日期间全球
Ms³7
.5大地震共发生
16次,发生于“磁暴月相二倍法
”得出的计算发震日期(在±5天的范围内
)有13次,无对应的有
3次。在15次计算发震日期中有11次对应全球
Ms³7.5
的大地震,2次对应
7.4£Ms³7.0
的大震。K指数大的起倍磁暴日与Ms
³
7.8
的巨震的发生有较好的相关性。这表明发生在月相的磁暴与
Ms³7.5
大地震的相关关系较好〔沈宗丕等,2001,待发表〕。
参考文献
1.
〔报道〕“利用磁偏角二倍法较好地预报了台湾
8级地震”,《地震战线》1972
年第3期。
2.
沈宗丕,“用磁偏角二倍法作中期预报的试探”,《地震》技术汇编
2,《地震战线》编辑部,1973
。
3.
张铁铮,“磁暴二倍法预报地震”,《自然科学争鸣》,
1975年第2
期P35~40
。
4.
耿庆国,“地磁二倍法研究的历史回顾”
—“太阳活动—
磁暴—强震活动关系分析”,《天文与自然灾害》,1991
年P61~67,地震出版社
5.
沈宗丕“
谈谈磁偏角二倍法”,《地震战线》1977
年第3期P30
~32。
6.
沈宗丕、徐道一,“应用磁暴月相二倍法对全球
Ms³7.5
大地震的预测效果分析“,《西北地震学报》
1996年9
月本18卷第3
期P84~86
。
7.
徐道一、王湘南、沈宗丕,“
1994年9
月底10月初Ms
³
7.5大地震的预测依据”,《地震危险性预测研究》(1995年度),北京地震出版社,
1994,P187
~191。
8.
徐道一、王湘南,沈宗丕,“磁暴与地震跨越式关系探讨”,《地震地质》,
1994年3
月第16卷第1
期P21~25
。
9.
沈宗丕、徐道一、张晓东、汪成民,“磁暴月相二倍法的计算发震日期与全球
Ms³7.5
大地震的对应关系”,2001年,待发表。