地球物理学进展  2016, Vol. 31 Issue (5): 2020-2026   PDF    
究竟有无可靠的大地震短临物理前兆?
杨百存, 秦四清, 薛雷, 吴晓娲, 张珂     
中国科学院地质与地球物理研究所 中国科学院页岩气与地质工程重点实验室, 北京 100029
摘要: 探寻可靠的地震前兆对大地震预测具有重要意义.地震前兆复杂多样,根据感知方式不同,可分为微观前兆和宏观前兆两大类.本文从大地震孕育的物理机制角度出发,结合典型震例,讨论前震事件与标志性地震事件之间已发现的可能地震前兆的可靠性.研究表明,已发现的物理异常现象可能反映了大地震孕育过程中应力的局部涨落效应,难以作为可靠的大地震短临前兆.我们强调,判断某种物理异常是否为可靠的大地震前兆,需进行长期观测才能得出明确结论.
关键词大地震     物理机制     标志性地震事件     前震     微观前兆     宏观前兆    
Do reliable short-imminent physical precursors exist before major earthquakes?
YANG Bai-cun , QIN Si-qing , XUE Lei , WU Xiao-wa , ZHANG Ke     
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
Abstract: It is very important to seek reliable earthquake precursors for earthquake prediction. The precursors can be divided into microscopic and macroscopic ones by human perception. Combined with typical earthquake cases, the reliability of possible earthquake precursors suggested by scholars between foreshock and characteristic earthquake event is discussed in the paper from the angle of the physical mechanism. The study indicates that the physical anomalies discovered previously may reflect the local fluctuation effect of stress during the seismogenic process of major earthquakes, but can hardly be regarded as reliable short-imminent earthquake precursors. We emphasize that the long-term observation about a certain physical anomaly is required for determining its reliability as earthquake precursor.
Key words: major earthquake     physical mechanism     characteristic earthquake event     foreshock     microscopic precursor     macroscopic precursor    
0 引言

《Science》杂志在创刊125周年之际,把“是否存在有助于预报的地震先兆”列为125个最具挑战性的科学问题之一,由此可见探寻可靠地震前兆的重要性.地震前兆指地震发生前出现的异常现象,本质上是岩石介质受力作用的产物(Main and Meredith, 1989严尊国和陈俊华,1999).大地震预测的科学探索思路是建立在地震有前兆,并可利用其预测未来地震发生时间、地点与强度的假定之上(张国民等,2001Cicerone et al., 2009).为寻找可靠的地震前兆,诸多学者(Dieterich,1979唐春安和徐小荷,1990尹祥础等,1994马瑾等,1996)基于岩石破裂实验研究,探寻应力(应变)等物理量与地震前兆信息之间的联系.图 1示出了岩样变形破坏过程中物理量响应特征,可看出当岩石损伤累积至体积膨胀点时,其物理和力学性质均发生显著变化(陈颙和黄庭芳,2001),这说明岩石宏观破裂前有前兆(耿乃光等,1986Rudajev et al., 2000),相应的大地震发生前亦应有前兆.断层运动导致岩石破裂发生地震,所以探寻地震前兆应直接从岩石破裂过程入手.秦四清等(2016abc)的研究表明,锁固段被加载至体积膨胀点时,震群事件开始发生,震群事件是监测技术手段可识别的唯一地震活动性前兆;体积膨胀点和峰值强度点对应的标志性地震事件时间间隔,一般为数十年或数百年,因此该地震活动性前兆可视为长期或超长期前兆.进一步的研究(李培等,2016)指出,前震事件可作为临近大地震发生的物理前兆.遗憾的是,前震发生后,即使某地震区已接近或到达临界状态,标志性地震事件并不会短期内发生,而是滞后数年乃至数十年,因此前震可视为标志性地震事件发生前的中长期或长期前兆,仍不是短临前兆.

那么,究竟有无可靠的大地震短临物理前兆呢?本文将从大地震孕育的物理机制角度出发(吴晓娲等,2016),讨论前震事件与标志性地震事件之间已发现的可能地震前兆的可靠性.

图 1 岩样变形破坏过程中监测物理量响应特征(据陈颙和黄庭芳(2001)修改) Figure 1 The response characteristics of monitored physical quantities in the deformation and failure process of rock specimen (modified after Chen and Huang (2001))
1 大地震前兆评述

地震前兆复杂多样,根据感知方式不同,可分为微观前兆和宏观前兆两大类(吴中海和赵根模,2013).

1.1 微观前兆

微观前兆指人的感官无法察觉、只有用专门仪器才能测量到的地震异常现象(张国民等,2001).孙强等(2014)类比岩石破坏的临界现象,将微观前兆分为两类:(1)与物理力学和结构特征变化有关的异常,包括地应力异常、地形变异常、地下水异常、波速异常和地电异常等;(2)与能量释放及物理场参数变化有关的异常,包括电磁辐射异常和水氡异常等.

1.1.1 与物理力学和结构特征变化有关的异常

在岩石破裂过程中,随着内部损伤累积,其承载能力降低,变形速率增大,渗透率增大,波速参数与电阻率值出现突变现象.这些异常现象能否作为大地震的短临前兆呢?以下将结合典型震例,从物理机制角度进行分析.

目前关于地震地电异常的成因机制主要有两种观点,即压电效应机制和渗滤电场机制(刘心恒,1978郝建国和张云福,2001).压电效应机制是指:随应力增加,矿物晶格发生扭曲形变而极化,致使矿物表面荷电产生电流引起区域电场异常;渗流电场机制是指:当应力增加到某种水平,岩石产生扩容变形,导致裂隙中的流体渗滤,因流体中离子电荷移动而形成渗滤电场引起区域电场异常(申俊峰等,2010).这两种机制本质上均由岩石所处应力环境改变引起(Morrow and Brace, 1981).2008年5月12日汶川大地震前,四川地区地电和地磁类观测结果表明(程万正等,2010),出现的异常现象比例极低,且未观测到电磁类异常现象,引起地震学界的广泛争议.秦四清等(2016b)基于孕震断层多锁固段脆性破裂理论,分析了汶川大地震的孕育过程,图 2示出了汶川地震区标志性地震事件之间的力学联系,可看出1976年8月四川松潘-平武间MS 7.2级双震震群是汶川大地震前的显著前震事件,前震发生后直至汶川大地震前,该区未发生ML>5.9级地震.换言之,汶川地震区自1976年前震事件后,第2锁固段已处于临界状态,但直至汶川大地震前,其物理状态并无显著改变,因此震前未观测到可识别的地电异常实属必然.

图 2 汶川地震区638.2.14 -2015.11.21之间CBS值与时间关系(数据分析时选取MS≥5.5级地震事件;误差修正已被考虑) Figure 2 Temporal distribution of CBS in the period from 14 February 638 to 21 November 2015 for the Wenchuan seismic zone (The earthquake events with MS≥5.5 are selected for data analysis. The error correction is also considered.)

再以地形变异常为例,分析其观测结果能否反映短临前兆.图 3示出了汶川大地震前5个跨断层观测场地的原始观测值时序,可看出除耿达水准外,其他时序曲线波动较大,难以找到形变异常与大地震的对应关系.若假设近震中场地(台站)观测到的显著震前断层预滑动为大地震前兆(周硕愚等,2009),则距震中最近的耿达场地确实观测到了可识别的异常现象.我们认为,这种异常现象是否为可靠地震前兆有待商榷,其原因之一是断层运动导致岩石破裂才能发生地震,即断层必须是活动的,但每次断层滑动并不一定能引发大地震;原因之二是前震事件后,标志性地震事件通常滞后数年乃至数十年才能发生,汶川大地震约滞后32年.换句话说,判断某种物理异常是否为可靠大地震前兆,需进行长期观测才能得出明确结论.对汶川地震区,若在1976年后开始观测水准变化,监测曲线可能出现多次显著“波动”,但多数情况下短期内并未有大地震发生.显然,若某种异常现象与大地震的对应关系不具有唯一性,则不能视为可靠的地震前兆.

图 3 龙门山断裂带5个跨断裂面短水准测站的观测时序(周硕愚等,2009) Figure 3 The time sequences of 5 fault-crossing short leveling sites at Longmenshan fault zone (Zhou et al., 2009)

图 4 唐山地震前水氡异常观测值曲线(据冯德益和郑熙铭(1987)修改) Figure 4 The observation value of radon anomaly before Tangshan earthquake (modified after Feng and Zheng (1987))
1.1.2 与能量释放及物理场参数变化有关的异常

在岩石变形破坏过程中,其部分能量以弹性能、声能、电磁能等形式耗散,伴随有氡辐射脉冲增加、电磁辐射强度增加与红外辐射增强等,这些异常现象在室内岩石力学试验中能被识别(孙强等,2014),但前震事件后至标志性地震事件前,震源体附近是否会产生此类异常,需结合震例进行分析.

以水氡异常为例,氡射气脉冲异常变化与加载过程中岩石的裂隙和表面积大量增加密切相关(王俊峰,2010).秦四清等(2016abc)对全球62个主要地震区的震情分析表明,前震事件后至标志性地震事件前,大多数地震区无较大地震事件发生,也即震源体无较大的新裂缝产生,故难以把水氡异常现象视为大地震短临前兆.

图 4示出了1976年唐山大地震前棉4、田疃和务1井孔水氡含量变化曲线,可看出“波动”具有随机性特点,其与大地震的关系具有不确定性(冯德益和郑熙铭,1987).秦四清等(2016b)基于孕震断层多锁固段脆性破裂理论,分析了1976年唐山大地震孕育过程,图 5示出了唐山地震区标志性地震事件之间的力学联系.可看出1969年7月18日渤海MS 7.4级地震是1次显著前震事件,该震发生后直至1976年唐山大地震前,该区未发生MS >5.5级地震,意味着震源体没有新的较大破裂面产生.如此,即使出现显著水氡异常现象,也不能视为大地震的短临前兆.

图 5 唐山地震区公元前1767-2015.11.21之间CBS值与时间关系(数据分析时选取MS≥5.0级地震事件;误差修正已被考虑) Figure 5 Temporal distribution of CBS in the period from B.C. 1767 to 21 November 2015 for the Tangshan seismic zone (The earthquake events with MS≥5.0 are selected for data analysis. The error correction is also considered.)
1.2 宏观前兆

宏观前兆指人的感觉器官能直接感知的地震异常现象(张国民等,2001),如中小地震、地下水异常、气象异常、动物异常、地声和地光等.

1.2.1 中小地震能作为大地震的短临前兆吗?

大地震发生前,每个地震区都会发生数量众多的中小地震,根据中小地震活动增强特征在短期内预测大地震发生,属于小概率事件.

根据“小震闹,大地震到”的经验,我国“成功”预报了1975年海城大地震.秦四清等(2016abc)的研究表明,只有在锁固段峰值强度点发生的标志性地震事件和某些预震(preshock)事件能被预测,而大多数预震为随机事件,只能判断其震级上限.海城大地震发生在珲春地震区(秦四清等,2016b),图 6示出了该地震区标志性地震事件之间的力学联系.可以看出,海城大地震是未来MS 8.0~8.3级地震前的一次显著随机预震事件,理论上其不能被预测.因此,所谓海城大地震的“成功”预报,纯属偶然,不能认为是一次依据科学方法成功预报的案例.

图 6 珲春地震区19.2-2016.1.3之间CBS值与时间关系(数据分析时选取MS≥6.5级地震事件;误差修正已被考虑) Figure 6 Temporal distribution of CBS in the period from February 19 to 3 January 2016 for the Hunchun seismic zone (The earthquake events with MS≥6.5 are selected for data analysis. The error correction is also considered.)
1.2.2 地下水位异常

大地震前地下水水位的异常变化,是震源区及其附近岩石的突然变形所致(郭增建等,1974Roeloffs et al., 1989).2008年汶川大地震发生前后,对地下水位的观测结果加深了以往对大地震宏观异常的认识(陈立德和付虹,2010).图 7示出了汶川地震前泸州13号井水位日均值,未观测到显著的异常变化.我们认为其原因在于:(1)1976年前震事件后至汶川大地震前,汶川地震区未发生ML>5.9级地震,即震源体无较大的破裂变形产生,故不会出现明显的异常现象.(2)地下水位的动态变化,可能是区域构造活动的作用结果,也可能是水文因素的干扰所致(孙小龙,2016),如何排除干扰并提取有效前兆异常信息,目前仍是地震监测预报的一大难题(车用太等,2011).鉴于此,地下水异常仍不能作为可靠的大地震短临前兆.

图 7 泸州13号井水位日均值(据陈立德和付虹(2010)修改) Figure 7 The daily value curves of No. 13 well water level in Luzhou (modified after Chen and Fu (2010))
1.2.3 其他宏观异常

对于气象、动物、地声和地光等异常现象,目前尚缺乏合理的物理机制解释(Rikitake,1994),需借助统计检验分析判断这些异常现象作为地震前兆的可靠性.应该指出,统计检验无法对观测现象本身是否为地震前兆给出肯定或否定的结论,但观测现象若要被视为可信地震前兆,则必须首先通过统计检验(吴忠良,1999).

卢双苓等(2015)以《中国震例》为主要参考资料,统计了中国大陆(范围: 17. 0°~55. 0°N,70°~150°E)在1966-2002年间,发生的1576次M≥5.0级地震的宏观异常现象,其中动物习性异常58次、地光异常18次、地声异常17次和气象异常9次,分别占地震总数的3.68%、1.14%、1.08%和0.57%,表明这些宏观异常现象与大地震的相关性很差,不能作为可靠的大地震短临前兆.

2 讨论与结论

多年来,前人基于诸多地震前观测到的异常现象,总结出多种地震孕育模式和地震预测方法.但目前尚未发现可重复观测到的地震短临前兆,也没有可信的大地震短临预测成功案例,其原因在于:

(1)异常现象能否作为大地震前兆仍待考证.某物理量监测值在正常背景值范围内出现异常变化后,应首先鉴别其是否为地震异常(张国民等,2001). 秦四清等(2016abc)的研究表明,每个地震区主震事件的孕育历时均很长,一般为数千年乃至上万年;标志性地震事件的时间间隔一般为数十年或数百年.目前各种异常现象的观测时间最长仅为几十年,通常为震前几年或几个月,因此我们对事后分析认为的异常是否为可信大地震前兆持严重怀疑态度,估计某些异常可能反映了大地震孕育过程中应力的局部涨落效应,而不是大地震前兆.再次强调,判断某种物理异常是否为可靠的大地震前兆,需进行长期观测才能得出明确结论.否则,就会陷入“事后总能找到异常,而事前则难以判明某种异常是否为地震前兆”的怪圈,且难以自拔.

(2)产生异常现象的物理机制不明确.若认为某种异常现象是地震前兆,必须严格论证异常现象的真实性,且阐明其产生机制,并接受可重复性检验,如此才能找到可靠的地震前兆.秦四清等(2016abc)基于大量震例分析发现,前震事件后至标志性地震事件前,处于临界状态的大多数地震区再无较大地震发生,即震源体物理状态基本不变,相应地监测物理量也基本不变,故不存在可识别的短临物理前兆.

(3)可能混淆了异常现象与大地震发生前、震时与震后的对应关系.以电磁或电离层扰动异常为例,若在某次大地震震时或震后短期内观测到显著异常,容易从物理机制上给予解释.然而,根据上述分析知,标志性地震事件发生前,大多数地震区长期处于物理稳定状态,若大地震前能观测到显著异常,不仅不能证明其是短临前兆,反而证明是太阳高能辐射扰动所致.

综上分析,目前已发现的物理异常均不能作为可靠的大地震短临前兆.至于是否存在有助于地震预报的化学前兆或其他尚未可知前兆,本文并未涉及,仍需进一步探索.我们对前震与标志性地震事件发生时间关系的研究表明,从岩石剪切破裂的粘滞性入手,可能是突破大地震短临预测的关键所在.

致谢 感谢国家自然科学基金资助项目(41572311;41302233)对研究工作的资金支持.
参考文献
[] Che Y T, Yu J Z, Liu C L, et al .2011. Principles on distinguishing interference from seismic precursor of underground water variation and its application[J]. Acta Seismologica Sinica (in Chinese), 33 (6) : 800–808. DOI:10.3969/j.issn.0253-3782.2011.06.010
[] Chen L D, Fu H .2010. The contrast study of Wenchuan M 8.0 earthquake precursor anomalies[J]. Inland Earthquake (in Chinese), 24 (4) : 289–297. DOI:10.3969/j.issn.1001-8956.2010.04.001
[] Chen Y, Huang T F .2001. Rock Physies[M]. Beijing: Peking University Press .
[] Cheng W Z, Ren Y X, Wu X P .2010. Applying the module intelligent sensor sentinel to the power supply in seismic stations[J]. Earthquake Research in Sichuan (in Chinese) (2) : 1–9. DOI:10.3969/j.issn.1001-8115.2010.02.001
[] Cicerone R D, Ebel J E, Britton J .2009. A systematic compilation of earthquake precursors[J]. Tectonophysics, 476 (3-4) : 371–396. DOI:10.1016/j.tecto.2009.06.008
[] Dieterich J H .1979. Modeling of rock friction:1. Experimental results and constitutive equations[J]. Journal of Geophysical Research:Solid Earth, 84 (B5) : 2161–2168. DOI:10.1029/JB084iB05p02161
[] Feng D Y, Zheng X M .1987. Analysis on the characteristic of precursory anomaly and the fuzzy recognition of the variations of radon contents in ground water before the great Tangshan eathquake in 1976[J]. Earthquake Research in China (in Chinese), 3 (1) : 40–48.
[] Geng N G, Li J H, Hao J S, et al .1986. Precursory deformation and stress drop in rupture experiments of brittle rocks[J]. Earthquake (in Chinese) (1) : 9–12.
[] Guo Z J, Qin B Y, Feng X C .1974. Discussion on the change of ground-water level preceding a large eathquake from an earthquake source model[J]. Acta Geophysica Sinica (in Chinese), 17 (2) : 99–105.
[] Hao J G, Zhang Y F .2001. Earthquake Static Electricity for the Prediction[M]. Beijing: Petroleum University Press .
[] Li P, Qin S Q, Xue L, et al .2016. Identification method of foreshock events[J]. Progress in Geophys (in Chinese), 31 (4) : 1450–1455. DOI:10.6038/pg20160406
[] Liu X H .1978. Discussion on earthquake prediction by using underground current anomaly variation[J]. Journal of Seismological Research (in Chinese) : 1–15.
[] Lu S L, Qu B A, Cai Y, et al .2015. Statistical analysis of the relationship between macroscopic anomaly and earthquake[J]. Earthquake Research in China (in Chinese), 31 (1) : 141–151.
[] Ma J, Ma S L, Liu L Q, et al .1996. Geometrical texture of faults evolution of physical field and instability characteristics[J]. Acta Seismological Smica (in Chinese), 18 (2) : 200–207.
[] Main I G, Meredith P G .1989. Classification of earthquake precursors from a fracture mechanics model[J]. Tectonophysics, 167 (2-4) : 273–283. DOI:10.1016/0040-1951(89)90078-4
[] Morrow C, Brace W F .1981. Electrical resistivity changes in tuffs due to stress[J]. Journal of Geophysical Research:Solid Earth, 86 (B4) : 2929–2934. DOI:10.1029/JB086iB04p02929
[] Qin S Q, Li G L, Xue L, et al .2014. Analysis of the future earthquake situation for some seismic zones in the northeast China, north China and Taiwan[J]. Progress in Geophys (in Chinese), 29 (2) : 535–554. DOI:10.6038/pg20140210
[] Qin S Q, Li P, Yang B C, et al .2016a. The identification of mainshock events for main seismic zones in seismic belts of the Circum-Pacific, ocean ridge and continental rift[J]. Progress in Geophysics (in Chinese), 31 (2) : 574–588. DOI:10.6038/pg20160209
[] Qin S Q, Yang B C, Wu X W, et al .2016b. The identification of mainshock events for some seismic zones in mainland China (II)[J]. Progress in Geophysics (in Chinese), 31 (1) : 115–142. DOI:10.6038/pg20160114
[] Qin S Q, Yang B C, Xue L, et al .2016c. The identification of mainshock events for main seismic zones in the Eurasian seismic belt[J]. Progress in Geophysics (in Chinese), 31 (2) : 559–573. DOI:10.6038/pg20160208
[] Rikitake T .1994. Nature of Macro-Anomaly Precursory to an Earthquake[J]. J. Phys. Earth, 42 : 149–163. DOI:10.4294/jpe1952.42.149
[] Roeloffs E A, Burford S S, Riley F S, et al .1989. Hydrologic effects on water level changes associated with episodic fault creep near Parkfield, California[J]. Journal of Geophysical Research:Solid Earth, 94 (B9) : 12387–12402. DOI:10.1029/JB094iB09p12387
[] Rudajev V, Vilhelm J, Ek T L .2000. Laboratory studies of acoustic emission prior to uniaxial compressive rock failure[J]. International Journal of Rock Mechanics and Mining Science, 37 (4) : 699–704. DOI:10.1016/S1365-1609(99)00126-4
[] Shen J F, Shen X H, Liu Q .2010. The thermoelectric effect of magnetite:a new model for abnormal geo-electricity from the formation and occurrence of earthquake[J]. Journal of Mineralogy and Petrology (in Chinese), 30 (4) : 21–27.
[] Sun Q, Qin S Q, Xue L, et al .2014. Critical phenomena and information identification of rock failure[M]. Xuzhou: China University of Mining and Technology press .
[i] Sun X L. 2016. Study on the relationship between dynamic changes of groundwater and seismic zctivity[D](in Chinese). Beijing:China University of Geosciences.
[] Tang C A, Xu X H .1990. Precursory deformation rate and stress drop in fracture experiments of brittle rocks and earthquake precursory[J]. Crustal Deformation and Earthquake (in Chinese), 10 (3) : 12–18.
[] Wang J F. 2010. Random migration in overlying strata during spontaneous combustion of coal underground and its application[D](in Chinese). Taiyuan:Taiyuan University of Technology.
[] Wu X W, Qin S Q, Xue L, et al .2016. Physical mechanism of major earthquake by earthquake cases[J]. Chinese J. Geophys. (in Chinese), 59 (10) : 3696–3710. DOI:10.6038/cjg20161016
[] Wu Z H, Zhao G M .2013. The earthquake prediction status and related problems:A review[J]. Geological Bulletin of China (in Chinese), 32 (10) : 1493–1512.
[] Wu Z L .1999. Seismological Problems of the Statistical Test of Earthquake Precursors-A Review on the Current Discussions about Earthquake Prediction (II)[J]. Earthquake Research in China (in Chinese), 15 (1) : 14–22.
[] Yan Z G, Chen J H .1999. The difficulty and limitation of the anomaly signal recognition of earthquake precursor observation[J]. South China Journal of Seismology (in Chinese), 19 (1) : 8–14. DOI:10.13512/j.hndz.1999.01.002
[] Yin X C, Chen X Z, Song Z P, et al .1994. The load-unload response ratio theory:a new approach to earthquake prediction[J]. Acta Geophysica Sinica (in Chinese), 37 (6) : 767–775.
[] Zhang G M, Fu Z X, Gui X T, et al .2001. Introduction to earthquake prediction[M]. Beijing: Science Press .
[] Zhou S Y, Wu Y, Shi S Y, et al .2009. Anomalous change of fault deformation before Wenchuan MS 80 earthquake and its comparison with response of other large earthquakes[J]. Acta Seismologica Sinica (in Chinese), 31 (2) : 140–151. DOI:10.3321/j.issn:0253-3782.2009.02.003
[] 车用太, 鱼金子, 刘成龙, 等.2011. 判别地下水异常的干扰性与前兆性的原则及其应用实例[J]. 地震学报, 33 (6) : 800–808. DOI:10.3969/j.issn.0253-3782.2011.06.010
[] 陈立德, 付虹.2010. 汶川8.0级地震前兆异常对比研究[J]. 内陆地震, 24 (4) : 289–297. DOI:10.3969/j.issn.1001-8956.2010.04.001
[] 陈颙, 黄庭芳.2001. 岩石物理学[M]. 北京: 北京大学出版社 .
[] 程万正, 任越霞, 吴小平.2010. 对汶川8.0级地震前四川地区地电地磁观测异常的分析与思考[J]. 四川地震 (2) : 1–9.
[] 冯德益, 郑熙铭.1987. 1976年唐山地震前水氡含量变化的模糊识别和前兆异常的特性分析[J]. 中国地震, 3 (1) : 40–48.
[] 耿乃光, 李纪汉, 郝晋升, 等.1986. 脆性岩石破裂的前兆变形和前兆应力降[J]. 地震 (1) : 9–12.
[] 郭增建, 秦保燕, 冯学才.1974. 从震源孕育模式讨论大地震前地下水的变化[J]. 地球物理学报, 17 (2) : 99–105.
[] 郝建国, 张云福.2001. 地震静电预测学[M]. 北京: 石油大学出版社 .
[] 李培, 秦四清, 薛雷, 等.2016. 前震事件判识方法[J]. 地球物理学进展, 31 (4) : 1450–1455. DOI:10.6038/pg20160406
[] 刘心恒.1978. 利用地下电流异常变化预报地震的探讨[J]. 地震研究 : 1–15.
[] 卢双苓, 曲保安, 蔡寅, 等.2015. 宏观异常与地震关系的统计分析[J]. 中国地震, 31 (1) : 141–151.
[] 马瑾, 马胜利, 刘力强, 等.1996. 断层几何结构与物理场的演化及失稳特征[J]. 地震学报, 18 (2) : 200–207.
[] 秦四清, 李国梁, 薛雷, 等.2014. 东北、华北与台湾地区某些地震区地震趋势分析[J]. 地球物理学进展, 29 (2) : 535–554. DOI:10.6038/pg20140210
[] 秦四清, 李培, 杨百存, 等.2016a. 环太平洋、大洋海岭与大陆裂谷地震带主要地震区主震事件判识[J]. 地球物理学进展, 30 (2) : 574–588. DOI:10.6038/pg20160209
[] 秦四清, 杨百存, 吴晓娲, 等.2016b. 中国大陆某些地震区主震事件判识(II)[J]. 地球物理学进展, 31 (1) : 115–142. DOI:10.6038/pg20160114
[] 秦四清, 杨百存, 薛雷, 等.2016c. 欧亚地震带主震事件判识[J]. 地球物理学进展, 31 (2) : 559–573. DOI:10.6038/pg20160208
[] 申俊峰, 申旭辉, 刘倩.2010. 磁铁矿热电效应:地震地电异常的新模式[J]. 矿物岩石, 30 (4) : 21–27.
[] 孙强, 秦四清, 薛雷, 等.2014. 岩石破坏临界现象与信息识别[M]. 徐州: 中国矿业大学出版社 : 145 .
[] 孙小龙. 2016.地下水动态变化与地震活动的关系研究[D].北京:中国地质大学. http://cdmd.cnki.com.cn/Article/CDMD-11415-1016068641.htm
[] 唐春安, 徐小荷.1990. 岩石破裂的前兆变形率、前兆应力降及临震前兆[J]. 地壳形变与地震, 10 (3) : 12–18.
[] 王俊峰. 2010.煤地下自燃时覆岩中氡气运移规律及应用研究[D].太原:太原理工大学. http://cdmd.cnki.com.cn/Article/CDMD-10112-2010143547.htm
[] 吴晓娲, 秦四清, 薛雷, 等.2016. 基于震例讨论大地震的物理机制[J]. 地球物理学报, 59 (10) : 3696–3710. DOI:10.6038/cjg20161016
[] 吴中海, 赵根模.2013. 地震预报现状及相关问题综述[J]. 地质通报, 32 (10) : 1493–1512.
[] 吴忠良.1999. 地震前兆统计检验的地震学问题--对目前地震预测问题争论的评述(之二)[J]. 中国地震, 15 (1) : 14–22.
[] 严尊国, 陈俊华.1999. 地震前兆观测异常信号识别的困难与局限性[J]. 华南地震, 19 (1) : 8–14. DOI:10.13512/j.hndz.1999.01.002
[] 尹祥础, 陈学忠, 宋治平, 等.1994. 加卸载响应比--一种新的地震预报方法[J]. 地球物理学报, 37 (6) : 767–775.
[] 张国民, 傅征祥, 桂燮泰, 等.2001. 地震预报引论[M]. 北京: 科学出版社 .
[] 周硕愚, 吴云, 施顺英, 等.2009. 汶川8.0级地震前断层形变异常及与其它大地震的比较[J]. 地震学报, 31 (2) : 140–151. DOI:10.3321/j.issn:0253-3782.2009.02.003