人们普遍认为地震是自然界的突发事件，地震预报是全球的科学难题，而且有的地震学家断言地震不能预报，认为“处于自组织临界状态的大地，任何一次小地震都有可能灾变为一次大地震”，而“小地震发展成为大地震将决定于不仅仅是其断层附近，而且是其整个（震源体）空间的物理状态的无数细微结构”，而人们根本无法掌握深部无数细微结构的临界状态，因之地震根本不能预报[1]。
近半个世纪以来，曾有大量各种各样关于地震前出现异常现象的报道，如地电、地磁、地声、重力、应力（应变）、水位、倾斜、惰性气体散发、气象、高空离子浓度、太阳活动、大范围动物异常等等，有些现象不是专一的，有些现象是在震后报道的，每次地震异常的形态相差很大，有的仅在某一地区，并非整个震中区域都能观察到[1]。由于这些各具特点的、孤立的、离散的、间断的异常现象均不具备确定性，只能视为前兆现象，不能称其为地震前兆。那么地震究竟有没有前兆？有！笔者将地震前兆用孕震物理模型表示。
《1 孕震物理模型》
1 孕震物理模型
通过采用自行研制的 SW 系列地震前兆监测仪，在长期的观测中发现，5级以上的中、强地震从孕育到发生是一个完整的、具有阶段性的连续过程，每一阶段都有其固有的特征信息，地震是这个过程的结果。只要这个过程连续、完整，各阶段的特征信息充分、明确，就构成具有确定性的地震前兆。地震前兆不是前兆现象，但前兆现象可能是地震前兆某个阶段的特征信息。
《1.1 单一震源作用下的孕震物理模型》
1.1 单一震源作用下的孕震物理模型
图1是通过对长期观测中得到的多个震例（图 2）进行分析总结，得到了单一震源作用下的孕震物理模型。其中 ξ 为用 mV 记录的微应变， t 为以天为单位的时间，异常面积 S = ξ × t ，孕震过程的异常面积 S 与该地震发生时，波及监测台站的地震烈度 I 相关，即与震级 M 和震中距相关。孕震物理模型分为加载、相持、卸载、短临和临震五个阶段。加载和相持阶段的持续时间较长，一般在1个月到数月之间，趋势也比较一致；卸载和短临阶段的持续时间较短，一般在数天到数十天，且短临阶段的趋势各异（
，b，c，d )。导致短临阶段趋势各异的因素很多，如监测台站与震中的相对位置、震中距及地质条件等，然而，无论这些因素多么复杂，相同阶段的特征信息基本相同。室内岩石破坏试验的结果与实际观测总结得到的模型类似[2]，岩石破坏试验的结果反映的正是弹性物体在力的作用下，首先发生弹性变形，继而发生弹塑性变形，然后产生失稳，最终发生断裂即岩体破裂发生地震的力学机制。也是图１所示的孕震模型的物理基础。图2所示的在监测实践中记录到的典型孕震过程就是很好的见证。
《图 1 》
图 1 单一震源作用下的孕震物理模型
Fig.1 The physical model of earthquack gestation under single quack
《图 2》
图 2 实测的孕震过程
Fig.2 The recorded earthquacke gestation processes
所以，地震前兆是一个完整的、具有阶段性的连续过程，每一阶段都有其固有的特征信息。只要孕震物理模型的阶段性连续和完整，各阶段的特征信息充分和明确，就构成具有确定性的地震前兆。
《1.2 多震源作用下孕震物理模型》
1.2 多震源作用下孕震物理模型
由于发震构造宏伟，地震孕育过程可能有多处地应力集中，同时或先后发生地震；或者在不同发震构造上，也可能同时或先后发生地震，单一震源作用下的孕震物理模型在自然界中比较少见，所见到的是许多单一震源作用下的孕震物理模型的叠置，即多震源作用下孕震物理模型。任何一台仪器在同一时间里记录到的是在不同地区、不同地质条件下、孕育着不同发展阶段的多个大小不同地震前兆信息的总合（图 3 和图 4）。
《图 3 》
图 3 1975 — 1976 年包头市记录的多震源作用下孕震物理模型
Fig.3 The physical model of earthquake gestation under multiple quakes，recorded during 1975—1976，Baotou，China
《图 4 》
图 4 2003 年美国旧金山台记录的多震源作用下孕震物理模型
Fig.４ The physical model of earthquank gestation under multiple quakes，recorded in 2003，San Francisco，USA
从图3和图4中每分解出1个单一震源作用下的孕震物理模型就对应1个中、强地震，分解的原则是力学平衡和能量守恒（分解方法见另文）。
图3包括内蒙古和林格尔（6.3）、云南龙陵（7.5+7.6）、河北唐山（7.8）、四川松潘（7.2）、宁夏巴音木仁（6.2）、内蒙古阿巴嘎（5.3）、四川盐源（6.9）以及天津宁河（6.9）等地震的孕育信息，此外还有许多其他影响较小的远震或小震的孕育信息。
图4包括：2003 年 9 月和 11 月的 2 个 4 级小震群，2003 年 12 月 23 日 6.5 级地震和 2004 年强震加载阶段的信息。
比较图 3 和图 4，尽管形状不完全相同，但形态却极为相似。图3是在中国大陆块内记录到的，图4是在美国西海岸圣安德烈斯（San Andreas）大断裂带附近记录到的，两地相隔万余公里，地理地质条件差异很大，时隔 27 年,用笔者研发的仪器能够重复和再现；两处的地震类型不同，分别为块内破裂型地震和断裂粘滑型地震，也能够重复和再现。足以说明地震前兆在不同时间和不同地理地质条件下，对于不同类型的地震都能够重复和再现。用孕震物理模型表示的地震前兆，一切都是相对的，由于地壳处于自组织临界状态及其对力作用的线性和非线性响应，使得绝对尺度失去意义。
孕震物理模型（地震前兆）的连续性和阶段性及其各阶段的特征信息，是地震前兆具有确定性的判据。
《2 地震前兆的五个阶段及其特征信息》
2 地震前兆的五个阶段及其特征信息
《2.1 正常日变规律》
2.1 正常日变规律
识别地震前兆的首要问题在于要知道什么是正常和什么是异常。地震前兆各阶段的特征信息是在正常日变规律背景参照和比较下识别的。
用 24 h 为周期的极座标圆图显示日变规律，始点0为当地时间每天的0点，时间用角度表示，每 15° 为 1 h ，从 0 开始逆时针计时 360° 为 1 d ；用 mV 记录应力应变的变化幅度用半径表示。没有地震孕育时，正常的日变规律是一个基本闭合的、体现地壳脉动的、显现固体潮的偏心圆，偏心长轴与监测台站所处地理经度相关（图 5 )。
《图 5》
图 5 正常日变规律
Fig.5 The normal daily variation
《2.2 地震前兆的各个阶段及其特征信息的识别》
2.2 地震前兆的各个阶段及其特征信息的识别
通过对图6中3个震例的分析，来说明识别地震前兆各个阶段及其特征信息的方法。
2.2.1加载阶段加载阶段的特征是应变量随着时间的增长而增长，即使发生小震也不影响总的增长趋势。日变圆图开口向外扩张，半径逐渐增大显示应力应变量在逐渐增加（图 7 )，在此阶段，连续数天乃至数月的长图可以更好地表现加载趋势。
2.2.2相持阶段相持阶段的典型特征是出现大幅度的应变扰动，在日变圆图上表现为梅花图状（图 8），它是可能发生波及到监测台站有感地震的重要信息。这种梅花图扰动有时也会发生在加载阶段和卸载阶段。
《图 6》
图 6 地震前兆各阶段的特征信息
Fig.6 The characteristics of earthquake precursor at each stage
《图 7 》
图 7 加载阶段的特征信息-扩张
Fig.7 The characteristics of the loading stage-expanding
《图 8 》
图 8 相持阶段的特征信息-扰动“梅花图”
Fig.8 The characteristics of the stalemate stage-perturbating （The daily variation at this stage looks like plum blossom）
已经发现的异常扰动有4种形式：
1）超低频扰动。即扰动梅花图，是有感地震孕育的重要特征，不见扰动梅花图就不会发生波及观测台站的有感地震。
2）断层滑移。断层滑移信息的特征是成组的、周期性的、幅度变化不大的、呈张或压应力锯齿波形，是扰动梅花图的一种特例。在 2004 年 3 月 24 日内蒙古 5.9 级地震前；2005 年 1 月 19 日日本 7.0 级地震前；2005 年 3 月 20 日日本 7.0 级地震前；2005 年 3 月 29 日印度洋 8.5 级大震前，2005 年 3 月 31 日印尼 7.0 级地震前，都曾记录到断层滑移信息，显示断裂粘滑型地震的特征。
3）小震群。即“小震闹”，小震群也是应变量扰动的一种形式。
4）地动（也称“慢地震”）。地震仪记录不到，人能感觉到，时而出现在长周期地震仪的记录图上，频率很低，这种超低频波动在震前断断续续，规律性尚待探索。
一般来说，扰动幅度的大小与震级和震中距相关。震级越大，距离越近，扰动幅度亦越大（见表 1）。如果在台网的支持下同步观察，判断震中位置和震级并不困难。
《表 1 》
表 1 梅花图扰动程度与震级、震中距和波及观测台站的烈度关系
Table 1 The relationship between the size of perturbation and magnitude ， epicenter distance ， and earthquake intensity at the site of observation
2.2.3 卸载阶段卸载阶段的特征是应变量随着时间的增长而降低，恰与加载特征相反，日变圆图开口向内收缩，即半径逐渐减小显示应力应变量在逐渐减小（图 9），该阶段用趋势性长图可以得到更明确的表示。持续出现卸载信息，表示将要进入短临阶段。
《图 9 》
图 9 卸载阶段的特征信息收缩
Fig.9 The characteristics of unloading stage-shranking
2.2.4 短临阶段从卸载阶段到短临阶段的过渡是连续的，这时看不清地壳脉动，固体潮发生异变，日变圆图呈紧绷状态，是“震前平静”的一种表现，故称短临“绷圆”（图 10），这种现象是有感地震即将发生的重要标志，在此阶段，正常日变规律中的地壳脉动及周围环境的干扰已基本不可见。短临绷圆持续时间不等，长则几天或十几天，短则只有几个小时。
《图 10》
图 10 短临阶段的特征信息绷圆
Fig.10 The characteristics of short-term and impending stage-circle tightening
2.2.5 临震阶段临震阶段的典型特征是出现几个（或几组）大于脉动幅度３倍的脉冲，称为临震信号。它出现的时间很短，转瞬即逝，及时排除干扰，才能准确无误地识别和判断。地震多发生在临震信号出现后的72 h 之内。临震信号是不可或缺的、确定性的，在前述各阶段连续完整、特征信息明确的情况下，临震信号是地震前兆确定性的充分和必要条件（图 11 和图 12）。
《图 11 》
图 11 1976 年在北京用模拟仪表记录的临震信号
Fig.11 The Signal of Impending Quake Recorded with Analog Instrument in 1976，Beijing，China
《图 12 》
图 12 用数字式仪器记录的临震信号
Fig.12 The signal of impending quake recorded with digital instrument
《2.3 临震信号的识别》
2.3 临震信号的识别
临震信号发生在扰动梅花图和短临绷圆之后；地震仪记录不到的超低频脉冲，幅值大于3倍脉动幅度，要排除小震和人为干扰；临震信号多为 2 个以上，值得注意的是第2个临震信号具有确定性。
图13为临震信号的重复和再现。
图14是一个完整的实际震例（详见另文），时间顺序为自上而下由左至右，特征信息明确，易于识别，只要熟悉特征信息，及时看图分析，预报地震并不困难。
《图 13 》
图 13 4 个临震信号所对应的 4 次强震（2004.12.01 — 2005.03.23 由新民台记录）
Fig.13 The four signals of impending quake corresponding four strong earthquakes （Recorded at Xinming from Dec.1，2004 to Mar. 23，2005）
《图 14》
图 14 1999 年 11 月 29 日辽宁岫岩 5.6 级地震，营口台模拟仪表记录，Δ43 km ，台站烈度 6 度
Fig.14 The earthquake of Ms 5.6，Xiuyan，Liaoning，Nov.11，1999（recorded at Yingkou，epicenter distance is 43 km，the earthquake intensity at the site is degree of 6）
综上所述，这种按地震孕育阶段追踪特征信息的方法，就是在实践中总结出来的“追踪地震前兆信息，循序渐进预报地震”的分析方法。形象思维，看图识字，一目了然。
《3 监测仪器及有效监测范围》
3 监测仪器及有效监测范围
监测地震前兆的 SW 系列仪器都是自行研制的，包括应力仪、谐振仪、倾斜仪、地磁仪和大地电位仪。以应力仪为主，其他仪器辅助记录各种物理前兆现象，补充和印证特征信息。
应力传感器（图 15）应用了笔者研制的高灵敏度微位移传感元件，其特殊性能是：
具有对外力微小作用的敏感性－灵敏度高；具有自组织能力－能自主适应地震发生前后的地壳活动特性；
对力的作用具有线性和非线性响应－与地壳的固有特性一致；
具有能与地壳和地震过程相耦合的可能性（慎重选择台址和使传感器与土层密切耦合）。
《图 15 》
图 15 应力传感器
Fig.15 The stress sensor
仪器的有效监测范围是 5
Ms 9。即 5 级为 100 km ，6 级为 200 km ，7 级为 500 km ，8 级为 1000 km ，9 级为 3500 km 。
在有效监测范围内，科学配置仪器，合理分布台网，慎重选择台址，能给出地震三要素：时间（几天）、震级（±1.0 级)和震中位置（±2/3 台距）。
超出有效监测范围记录到的前兆不够完整，只能视为前兆现象，没有确定性。小于 5 级的地震，尤其是 4 级以下的地震，很难观测到完整的地震前兆。
以监测6级地震为主的台网，平均台距为 200 km ，5 级地震可能会丢失；若台距 100 km 就不会丢失 5 级地震；地震重点监视区台网的台距以 70~80 km 为宜。
《4 结语与讨论》
4 结语与讨论
笔者用 SW 系列地震前兆监测专用仪器对地震前兆进行了 30 年的观测，发现了地震孕育的全过程。震级和震中位置的判断是根据仪器有效监测范围、孕震物理模型异常面积和地震前兆扰动程度确定的，科学配置的台网有助于确定震中位置和震级。地震预报研究有可能从经验性预测、统计性预测提高到以物理机制为基础的确定性预测。地震前兆在不同时间、不同地区、不同地质条件下能够重复和再现，破坏性地震是可以预测的。
临震信号是不可或缺的、确定性的，是地壳失稳的标志，还是地震已经发生了，目前还不能确定，有待深入研究。
是什么机制能使仪器接收到数千公里之外的地震孕育过程？如果地震孕育过程的信息是从震源区发出的，是怎么传递到数千公里之外的？
不能因为这些现象目前还得不到确切的解释而放弃扩大试验和实施地震预测，深入研究有助于探索地震成因。
致谢：包头钢铁设计院对地震观测点的支持和投入，国家地震局领导的支持和鼓励，辽宁省地震局的合作试验，中科院物理研究所和北京电业中学的支持与合作，特别是美国科学家为在圣安德烈斯断层附近的试验提供了难得的条件和支持，是本报告获得成果的根本保证。笔者借此机会向有关方面致以敬意。特别要感谢李院士和许绍燮院士为本文进行了严肃认真的指导，在此一并致谢。
京公网安备 11010502051620号
