光学天文时纬残差异常在强震预测中的实践
胡辉1, 李语强1, 苏有锦2, 尹志强3, 韩延本3, 王锐1, 王博3     
1. 中国科学院云南天文台, 云南 昆明 650011;
2. 云南省地震局, 云南 昆明 650224;
3. 中国科学院国家天文台, 北京 100012
摘要: 为了进一步研究光学天文时纬残差异常在地震预测中的作用,介绍了光学天文时纬残差震前异常的发现和在地震预测中的研究实践,以及云南天文台光电等高仪的时纬残差异常变化与其周邻强震的对应关系,最后以讨论的形式给出了这种关联的可能的地球物理机制,目前在地震预测中的局限性和可能的解决途径。重要的是,2010年以来的预测实践进一步表明,利用光学天文时纬残差的同步异常提供地震预测信息,既没有虚报,也没有漏报。这就说明,它是完全可以作为地震预测的一种手段,投入地震预测实践中,值得更多重视和更加深入的研究。
关键词: 光学天文时纬残差     铅垂线变化     地下物质运动     地震    
Practice of Anomalies of Optical Astronomic Time-latitude Residuals in Strong Earthquake Prediction
Hu Hui1, Li Yuqiang1, Su Youjin2, Yin Zhiqiang3, Han Yanben3, Wang Rui1, Wang Bo3     
1. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650011, China;
2. Seismological Bureau of Yunnan Province, Kunming 650224, China;
3. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China
Abstract: To investigate the practical application of the anomalies of optical astronomic time-latitude residuals in earthquake prediction, the article introduces the discovery of the anomalies of optical astronomic time-latitude residuals before earthquake, its research practice in earthquake prediction, and the actual relation between the anomalies of the time-latitude residuals of photoelectric astrolabe in Yunnan Observatory and the strong earthquakes around, and finally discusses its possible geophysical mechanism, current limitations in earthquake prediction and possible solution methods. The forecasting practice since 2010 especially shows that neither does it generate false forecasts, nor misses forecasts of major earthquakes. This indicates that it is can be used as a method of earthquake prediction for practical usage, such work requires more attention and further study.
Key words: Optical astronomic time-latitude residuals     Plumb line variation     Motion of groundmass     Earthquake    

1966年以来,一个问题一直困扰着北京天文台(即现在的国家天文台)授时组的科研人员,观测结果一直很好的目视等高仪,为什么在1966年2~5月期间如此反常。所谓反常,就是这台仪器测定的地球自转参数(主要是地球自转速率)与全国其它仪器测定的综合值的偏差超出了正常的误差范围。到了1976年,这个问题又一次出现。正在大家感到困惑不解的时候,1976年7月28日唐山地震爆发了。这次大地震不仅震憾了中国,而且震动了全世界。正当全国上下支援唐山人民抗震救灾的同时,全国的科研人员也抓紧对地震的研究。正是在这种背景下,北京天文台的科研人员突然意识到,1966年与1976年北京天文台目视等高仪的观测结果反常,可能与1966年的邢台地震和这次唐山地震有关。因为这两次地震发生时,北京都有明显的震感,并且邢台地震震中位于北京天文台的西南,相距约300 km,而唐山地震震中位于北京天文台的东偏南,相距只有160 km。在经过一番认真思考讨论之后,他们潜身心对观测数据进行仔细的分析,终于发现了北京天文台光学天文测定的时间和纬度在强震前的异常现象。为了进一步证实这一现象的客观性和真实性,他们收集了全球10个天文台站的11台观测稳定、质量较好的仪器自1966年以来的观测结果,并进行了系统的分析。结果发现,在无震的年份和无震的台站,观测结果基本上没有出现这种异常现象,从而肯定了异常现象与地震有关,并称其为强烈地震前天文时纬残差的短期异常[1]。这一发现说明,精确的光学天体测量结果不仅可用于基本天体测量和地球自转参数的测定,而且可能为台站附近的地震预测提供警戒信息。

研究结果得到了国内外同行的重视,地震学家、中国科学院院士、原国家地震局副局长陈顒认为,该研究属于探索地震前兆复杂性方面涌现的“新思路、新方法、新理论”,自1985年以来,国际天文学联合会大会在有关光学天体测量仪器未来工作的决议中多次指出,可以开展垂线变化和地震预测研究。

在1990年11月举行的第5届全国地球动力学学术讨论会上,王仁、马宗晋等30多位地震、地球动力学专家紧急呼吁国家地震局和有关部门从速采取措施,使我国各天文台的光学天体测量手段能够继续为地震预测作贡献。他们还指出,“地震预测实践表明,光学天体测量观测结果是一项有希望的预测手段”。

在北京天文台对天文观测异常与强震的可能关系进行深入研究的同时[2-4],上海天文台、云南天文台也开展了此项研究工作[5-9];上海天文台发现在1990年6月5日江苏太仓5级地震前该台中星仪的测时残差与同期光电等高仪的测时测纬残差均有明显异常[6]

云南天文台由于地处地震多发地区,自1985年禄劝地震以后,与云南省地震局合作,开展了云南天文台光电等高仪的时纬残差异常与其周邻强震的相关研究,并为地震预测提供过几次重要参考信息[7],2010年以来,更是为周邻地区的几次强震作了成功的预测。

1 云南天文台光电等高仪的时纬残差与周邻的强震

1985年4月18日云南省禄劝县发生了6.3级地震,震中距云南天文台只有93 km,昆明震感强烈。地震发生以后,北京天文台的科研人员和我们分析了震前云南天文台光电等高仪的观测结果,发现在1985年3月的观测结果中有明显的异常[5]

为了深入探讨云南天文台光电等高仪的时纬残差异常与周邻强震的关系,与云南省地震局合作,按照光学天文时纬残差的定义与数据处理方法,全面、系统地分析了云南天文台光电等高仪自有观测记录以来的全部观测结果,发现二者之间有良好的对应关系,而且发现时间残差异常与纬度残差异常同步与否,是有震与无震的重要标志。这些结果已列于文[7]的表 1中。为了使读者对二者之间的关系有更完整的了解,本文再次引用文[7]表 1的结果,并且随着时间的推移增加了更多实际证例。

表 1 云南天文台的时纬残差异常变化与周邻强震的对应关系 Table 1 A The actual correspondence of the anomalies of the time-latitude residuals and the strong earthquakes around YO
发震日期 地点 φ λ M D A S(t) S(φ) Y/N T(d)
1 1976-05-29 龙陵 24°33′ 98°45′ 7.4 411 262 >3σ >2σ Y 120
2 1976-11-07 宁蒗 27°30′ 101°05′ 6.7 322 328 无观测
3 1978 >2σ >3σ 无≥6.0 地震
4 1979-03-15 普洱 23°07′ 101°15′ 6.8 265 216 >2σ >3σ Y 100
5 1983-06-24 莱州 21°14′ 103°10′ 7.0 425 173 >3σ >3σ Y 160
6 1984-04-24 孟连 22°04′ 99°14′ 6.0 489 227 >2σ >2σ Y 70
7 1985-04-18 禄劝 25°52′ 102°51′ 6.3 93 4 >2σ >3σ Y 60
8 1986 >2σ >2σ 无≥6.0 地震
9 1988-11-06 澜沧 22°50′ 99°43′ 7.6 372 231 >3σ >3σ Y 170
10 1989-05-07 耿马 23°29′ 99°29′ 6.3 377 242 >2σ >3σ Y 30
11 1992-04-23 孟连 22°25′ 98°56′ 6.8 464 233 >3σ >3σ Y 90
12 1993-01-27 普洱 22°56′ 101°05′ 6.3 291 216 >3σ >3σ Y 50
13 1995-07-12 孟连外 22°00′ 99°18′ 7.3 476 225 >2σ >3σ Y 100
14 1995-10-24 武定 25°46′ 102°16′ 6.5 110 330 无观测
15 1996-02-03 丽江 27°18′ 100°06′ 7.0 360 310 >2σ >3σ Y 80
16 2003-07-21 大姚 25°57′ 101°14′ 6.2 210 315 >3σ >2σ Y 80
17 2008-05-12 汶川 31°00′ 103°24′ 8.0 670 24 >2σ >2σ Y 50
18 2013-04-20 芦山 30°18′ 103°00′ 7.0 580 12 >2σ >2σ Y 45
19 2014-05-30 盈江 25°00′ 97°48′ 6.1 480 280 >2σ >2σ Y 80
20 2014-08-03 鲁甸 27°06′ 103°18′ 6.5 240 20 >2σ >2σ Y 65
21 2014-10-07 景谷 23°24′ 100°30′ 6.6 280 230 >2σ >2σ Y 35
表中包括发震日期、震区地点、地理纬度φ和地理经度λ、震级M、相对于观测仪器的距离D和方位A、时纬残差异常变化量S(t)和S(φ),以及震前预兆天数T(d)。Y/N表示异常的同步与否

正如在文[7]中指出的,光学天文时纬残差的异常与其它和地震相关的地球物理现象有相似之处,相当多的异常出现在地震前两三个月。统计结果表明,当地震的震级大于6.0时,光学天文时纬残差警戒信号可以涉及到震中距达300或400 km的地震,但当震中距更大时,例如700 km,一般无反映。

表 1可以看出,在2008年四川汶川8.0地震[7]]以前,有5次向云南省地震局提供了震前警戒信息,即表 1中的第11、12、13、15、16号震例,只因单台不成测网,还不能提供地震发生的三要素(即具体的时间、地点和震级)[7]

云南天文台所在地区的地质构造和相关的地震标示见图 1。图中只标出地震发生的日期和震级,没有标出地震发生的地点,因为在研究期间有多次地震发生在同一地点(如孟连)。

图 1 川滇菱形块体地质构造图 Figure 1 Figure of the geological structure of Xikang-Yunnan rhombic fault block (1)金沙江断裂;(2)红河断裂;(3)小江断裂;(4)安宁河断裂;(5)鲜水河断裂;(6)龙门山断裂
2 2010年以来的预测实践

自2010年2月,每个月的月初向云南省地震预报研究中心提供至上月底的光电等高仪时纬残差五日异常点日期,供他们分析研究和预测。

从2010年至今的几年时间里,在云南天文台周邻700 km以内,先后发生了4次大于6.0的地震,2013年4月20日四川芦山7.0地震,2014年云南省内3次大于6.0的地震,即5月30日盈江6.1地震、8月30日鲁甸6.5级地震和10月7日景东6.6级地震,这4次地震的震中距云南天文台分别是580 km、480 km、240 km和280 km,其方位分别是北偏东12°、西偏北10°、北偏东20°和南偏西48°(图 1)。这4次地震不同程度地造成了人员伤亡和财产损失。

同一期间,每月一报的云南天文台光电等高仪的时纬资料中,分别在这4次地震前的2013年3月和2014年3月、5月、9月有大于2σ的同步异常,而且仅有这4次同步异常,再无其它同步异常(见图 2)。所以从2010年以来报告的预警情况看,既没有虚报,也没有漏报。

图 2 2010.0~2016.0年云南天文台光电等高仪的时纬残差曲线 Figure 2 Residuals curves of the astronomical time (a) and latitude (b) of photoelectric astrolabe in YO from 2010.0 to 2016.0 ① 2013-04-20,芦山7.0;② 2014-05-30,盈江6.1;③ 2014-08-03,鲁甸6.5;④ 2014-10-07,景谷6.6
3 讨论

光学天体测量仪器是以地方铅垂线作为基准进行观测。因此,在观测结果中扣除地球整体运动(即地球自转参数影响)以后的部分,再经过适当的数据处理,消除偶然误差、星位误差、大气季节性等因素的影响以后,主要就是本地铅垂线的变化。而引起本地铅垂线变化的原因,主要可能来自震前地震孕育过程中的地下物质运动,如地下水活动和板内块体移动等[10-11],因此,它也可用于在较大范围内监测重力的变化。

光学天文时纬残差的同步异常,确实可以为台站周围即将发生的地震提前给出警戒信息,但是这种异常与地震发生之间的时间关系,目前也只是定性的。

表 1可以看到,异常信号的出现与大地震发生之间的时间延迟一般是2至3个月,但有时更长,特别是1988年孟连7.6地震与异常信号之间的时间延迟大约是170天。一些学者的研究结果表明,有些异常地球物理现象在地震前几个月、半年甚至更长的时间以前就出现了,例如1995年7月12日孟连7.3级地震,有些地球物理异常在1994年末就出现了。

在研究中还没有得到异常信号与地震震级、震中位置之间定量的关系,并且震级较小的地震和震中距较大的地震,在光学天文时纬残差中都还不能反映。如果用几架高精度的仪器,例如国家天文台研制的具有高精度、高自动化特征的CCD天顶望远镜[12],在地震活动区建立一个观测网,仪器之间相距约200~300 km,获得丰富的资料对于研究和估计地震的震级和地点将有积极的作用[11]

另外,光学天体测量仪器的观测只能在晴朗和少云的夜晚进行。在昆明,每年的5~10月是雨季,阴雨天严重影响光学天文观测资料的连续性(如图 2)。

总之,2010年以来的预测实践进一步表明,利用光学天文时纬残差的同步异常提供地震预测信息,既没有虚报,也没有漏报,光学天文观测的时纬残差异常可能是一种有效的地震前兆,值得进一步重视和更加深入的研究。

致谢:

感谢中国地球物理学会天灾预测专业委员会顾问、英籍华人陈一文博士和云南天文台李春晓博士的有益讨论和帮助。

参考文献
[1] 李致森, 张国栋, 张焕志, 等. 天文测时测纬的一种短时间异常与台站附近强震的对应关系[J]. 地球物理学报, 1978, 21(4): 278–291
Li Zhisen, Zhang Guodong, Zhang Huanzhi, et al. Relationship between the short-term anomaly of astronomical time-latitude observation and the strong earthquakes around the observatory[J]. Acta Geophysical Sinica, 1978, 21(4): 278–291.
[2] Li Z S, Han Y B, Tian J. A possible warning sign of major earthquake through observations of astronomical time and latitude[J]. Astrophysics and Space Science, 1991, 177(1): 383–384.
[3] 韩延本. 北京天文台Ⅱ型光电等高仪时纬残差异常与地震的相关性[J]. 科学通报, 1992, 37(12): 1105–1108
Han Yanben. Relationship between residual anomalies of time-latitude and earthquakes[J]. Chinese Science Bulletin, 1992, 37(12): 1105–1108.
[4] Li Z S, Han Y B, Tian J, et al. Monitoring of anomalies of residuals of astronomical time-latitude before major earthquakes[J]. Publications of Beijing Astronomical Observatory, 1992, 19: 58–63.
[5] 韩延本, 胡辉, 杜红荣. 禄劝6.3级地震前云南天文台时纬残差的初步分析[J]. 科学通报, 1986, 31(16): 1244–1246
Han Yanben, Hu Hui, Du Hongrong. Du Hongrong. Occurrence of short-period anomaly of residuals of astronomical time-latitude at Yunnan observatory preceding the Luquan earthquake (ML=6.3)[J]. Chinese Science Bulletin, 1986, 31(16): 1244–1246.
[6] 廖德春, 郑大伟. 江苏太仓五级地震与上海天文台测时异常[J]. 科学通报, 1990, 35(23): 80–81
Liao Dechun, Zheng Dawei. Earthquake (MS=5) in Taicang and anomalous residuals of time determination at Shanghai Observatory[J]. Chinese Science Bulletin, 1990, 35(23): 80–81.
[7] 阚荣举, 胡辉. 云南及其周邻大地震与云南天文台时纬残差变化[J]. 中国地震, 1995, 11(4): 299–309
Kan Rongju, Hu Hui. The major earthquakes in and around Yunnan and variation of astronomical time-latitude residuals determined in Yunnan Observatory[J]. Earthquake Research in China, 1995, 11(4): 299–309.
[8] 胡辉, 苏有锦, 付虹, 等. 大姚、普洱地震前云南天文台时纬残差异常[J]. 自然灾害学报, 2007, 16(5): 106–110
Hu Hui, Su Youjin, Fu Hong, et al. Anomalies in time-latitude residuals at Yunnan Observatory before Dayao and Puer earthquakes[J]. Journal of Natural Disasters, 2007, 16(5): 106–110.
[9] 胡辉, 韩延本, 苏有锦, 等. 2008~2009年云台光电等高仪时纬残差的异常波动[J]. 天文研究与技术——国家天文台台刊, 2011, 8(1): 91–94
Hu Hui, Han Yanben, Su Youjin, et al. Anomalous variations of daily residual fluctuations of astronomical time of the astrolabe in the Yunnan Observatory during 2008-2009[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2011, 8(1): 91–94.
[10] 张国栋. 强震前地下水活动引起的垂线变化[J]. 地震学报, 1981(2): 50–57
Zhang Guodong. The variation in the vertical line caused by the activity of the underground water before a strong earthquake[J]. Acta Seismological Sinica, 1981(2): 50–57.
[11] 胡辉, 林辉, 尹志强, 等. 汶川地震的天文背景[J]. 天文研究与技术——国家天文台台刊, 2009, 6(2): 161–166
Hu Hui, Lin Hui, Yin Zhiqiang, et al. Astronomical background of Wenchuan earthquake[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2009, 6(2): 161–166.
[12] 王博, 田立丽, 王政, 等. 数字化天顶望远镜观测图像及数据处理[J]. 科学通报, 2014, 59(12): 1100–1107
Wang Bo, Tian Lili, Wang Zheng, et al. Image and data processing of Digital Zenith Telescope (DZT-1) of China[J]. Chinese Science Bulletin, 2014, 59(12): 1100–1107.
由中国科学院国家天文台主办。
0

文章信息

胡辉, 李语强, 苏有锦, 尹志强, 韩延本, 王锐, 王博
Hu Hui, Li Yuqiang, Su Youjin, Yin Zhiqiang, Han Yanben, Wang Rui, Wang Bo
光学天文时纬残差异常在强震预测中的实践
Practice of Anomalies of Optical Astronomic Time-latitude Residuals in Strong Earthquake Prediction
天文研究与技术, 2017, 14(3): 401-406.
Astronomical Research and Technology, 2017, 14(3): 401-406.
收稿日期: 2016-08-31
修订日期: 2016-09-12

工作空间