地震预测

·编者按·

地震预测

·编者按·

地震预测是指同时指明未来地震的位置、大小、时间和概率四个参数的过程。其中，位置、大小和时间被称为“地震三要素”。尽管研究有所进展，到目前为止，地震预测仍然是公认的世界性科学难题。

地震是地壳运动的一种表现，是构造变动急剧的地块所发生的破坏活动。其在孕育的过程中，带有应力、应变、重、磁、电、地下水等异常值的出现。地震预测专家通常需要捕捉这些异常信息来给出地震预测所需的四个参数，并根据异常信息来源不同，采用不同的地震预测方法。主要包括：地应力观测法、钻孔应变法、GPS观测法、地下流体异常观测、氢同位素法、电磁法、电离层法、潮汐力谐振共振波法等。除此以外，人们还试图以各种自然现象预测地震，例如天气状况、地震云、井水水位及动物行为等等。但是多数方法是经验式的，很多方法尚缺乏能证实因果联系的科学理论。

尽管困难重重，自20世纪60年代以来，中长期地震预测还是取得了一些有意义的进展，比如板块边界大地震空区的确认、“应力影区”、地震活动性影像、图像识别、由帕克菲尔德地震预报试验场的预报实践获得的正反两方面的经验等。

本专题得到吴忠良研究员（中国地震局地球物理研究所）、李德威教授（中国地质大学）的大力支持。

·热点数据排行·

截至2016年4月20日，中国知网（CNKI）和Web of Science（WOS）的数据报告显示，以“地震预测”为词条可以检索到的期刊文献分别为2277条与1213条，本专题将相关数据按照：研究机构发文数、作者发文数、期刊发文数、被引用频次进行排行，结果如下。

研究机构发文数量排名（WOS）

作者发文数量排名（CNKI）

作者发文数量排名（WOS）

期刊发文数量排名（CNKI）

期刊发文数量排名（WOS）

根据中国知网（CNKI）数据报告，以“地震预测”为词条可以检索到的高被引论文排行结果如下。

国内数据库高被引论文排行

根据Web of Science统计数据，以“地震预测”为词条可以检索到的高被引论文排行结果如下。

国外数据库高被引论文排行

·经典文献推荐·

基于Web of Science检索结果，利用Histcite软件选取LCS（Local Citation Score，本地引用次数）TOP 30文献作为节点进行分析，得到本领域推荐的经典文献如下。

本领域经典文献

来源出版物: Earth-Science Reviews, 1968, 4: 245-282

Prediction of haicheng earthquake

Raleigh, B; Molnar, P; Hanks, T; et al.

Abstract: At 1936 LT, the evening of February 4, 1975. A magnitude 7.3 earthquake occurred near the town of Haicheng. Liaoning Province, in the People’s Republic of China. This earthquake was widely felt through northeastern China and caused coniderable damage to existing structures and cultivated land in the epicentral region. Approximately 90% of the structures in Haicheng, a city of 90000 people, were destroyed or seriously damaged at the time of the earthquake. Given the time of day, subfreezing ambient temperatures, and the population density through the region in which significant damage occurred, casualties in excess of 100000 would have ordinarily been anticipated—Had this earthquake not been predicted accurately enough to effect massive evacuation of the local population from their dwellings on the afternoon and evening of February 4. The prediction of the Haicheng earthquake is an extraordinary achievement of the seismological and geophysical workers of the People's Republic of China, whose national program in earthquake research was less than 10 years old at the time. At the same time the Haicheng earthquake was the first major shock to have been accurately predicted anywhere in the world, and our delegation is one of several groups of scientists whose trips to China have been highlighted by the accounts of this remarkable success.

来源出版物: Transactions-American Geophysical Union,1977, 58: 236-272

The Parkfield, California, earthquake prediction experiment

Bakun, WH; Lindh, AG

Abstract: Five moderate(magnitude 6)earthquakes with similar features have occurred on the Parkfield section of the San Andreas fault in central California since 1857. The next moderate Parkfield earthquake is expected to occur before 1993. The Parkfield prediction experiment is designed to monitor the details of the final stages of the earthquake preparation process; observations and reports of seismicity and aseismic slip associated with the last moderate Parkfield earthquake in 1966 constitute much of the basis of the design of the experiment.

来源出版物：Science, 1985, 229(4714): 619-624

Statistical short-term earthquake prediction

Kagan, YY; Knopoff, L

Abstract: A statistical procedure, derived from a theoretical model of fracture growth, is used to identify a foreshock sequence while it is in progress. As a predictor, the procedure reduces the average uncertainty in the rate of occurrence for a future strong earthquake by a factor of more than 1000 when compared with the Poisson rate of occurrence. About one-third of all main shocks with local magnitude greater than or equal to 4.0 in central California can be predicted in this way, starting from a 7 year database that has a lower magnitude cutoff of 1.5. The time scale of such predictions is of the order of a few hours to a few days for foreshocks in the magnitude range from 2.0 to 5.0.

来源出版物: Science, 1987, 236(4808): 1563-1567

Strategies in strong earthquake prediction

Molchan, GM

Abstract: The prediction problem of a stochastic point process in terms of loss function γ is considered. Losses depend on prediction errors:the fraction of failures-to-predict; andthe fraction of alarm time. The structure of γ-optimal prediction strategies for some functions γ is described. This allows analysis of the announced prediction of strong earthquakes in the southern part of the San Andreas fault. We consideras a property of prediction strategy and set problems on the description of ‘prediction capability’ and comparison of prediction algorithms. In addition, the distribution of empirical prediction errors for the renewal process model and the type of attenuation of inter-event time distribution, are considered.

来源出版物: Physics of the Earth and Planetary Interiors, 1990, 61(1): 84-98

·推荐综述·

鲁甸、景谷、康定地震预测的原理、方法及其意义

李德威

近年来青藏高原东缘相继发生系列强震，包括2008 年5月12日汶川Ms8.0级地震、2010年4月14日青海玉树Ms7.1级地震、2013年4月20日四川芦山Ms7.0级地震、2014年8月3日云南鲁甸Ms6.5级地震，2014 年10月7日云南景谷Ms6.6级地震，2014年11月22日四川康定Ms6.3级地震。这些沿着川滇菱形块体东北部边缘构造带串珠状分布的地震发震频率高、震源深度浅、链式灾害多、灾害损失重、关联研究少、预测环节弱。

地震预测不仅是科学问题，也是社会问题。《礼记·中庸》中说：“凡事预则立，不预则废”。搞好预测、作好准备是成功的基础。地震是地球物质运动的一种构造表现形式，自然界任何物质运动都是有规律的，取决于能量的源、汇、释。因此，地震活动是在某种能量作用下物质发生规律性运动的结果，必然会出现与地震机理密切关联的前兆异常，据此可以科学地预测地震。

笔者曾在东昆仑、汶川、玉树、通海地震区分别进行了3次、5次、4次、3次野外调查，在芦山、鲁甸、康定地震区各进行了1次实地调查，并与50多位地震专家进行了长达3年的广泛、深入交流。在此基础上，本文以下地壳层流与盆山耦合大陆动力学（李德威，1993，1995， 2003； 李德威和纪云龙，2000；李德威等， 2009；Li，2008）、多级物质循环地球系统动力学（李德威，2005，2011，2012，2014a，2014b）、热灾害链（李德威，2011，2012，2014a，2014b）和陆内地震热流体撞击（李德威， 1995， 2008a，2010，2011，2012，2014a，2014b；李德威等，2014）等学术思想为指导，通过芦山、鲁甸、景谷、康定等地震的长期和中期预测（李德威，2008a，2010，2011，2012；李德威等，2013a），进一步总结陆内地震预测的思路和方法，阐明地震的形成机理，强调取出灾能、变害为宝的重要意义。

1地震能够预测

自从Gelleret al.（1997）在Science上发表“地震不能预测”一文以来，主流基本认同“地震是一种自组织临界现象”，支持地震不可预测。也有一些学者对地震自组织临界提出质疑（Knopoff，1999；Mainand Kindy，2002；吴忠良和蒋长胜，2007；陈运泰，2009；许绍燮，2011）。Cyranoski（2004）在Nature上发表了“A seismic shift in thinking”（地震预测观念的大转变）一文，令人鼓舞。

我国成功预测了海城、松潘—平武、龙陵—潞西、乌恰等二十多次6级以上的地震，还创造了唐山大地震中的“青龙奇迹”。这些成功案例充分说明地震是可以预测的。正因如此，我国一些专家克服许多困难，长期致力于地震预测事业，并取得丰硕的研究成果（陈立德和付虹，2014；陈运泰，2009；傅承义等，1979；邓志辉和马瑾，1993；耿庆国，1984；郭增建，1997；胡敦宽等，1995；刘德富和康春丽，2003；强祖基等， 1990；任振球等，2001；许绍燮，2011；曾小苹和林云芳，1995；张小涛等，2009；张元生等，2010；郑治真，1994；Hu et al.，2013；Zhao and Qian，1994）。笔者曾与马宗晋、许绍燮、李玶、郭增建、耿庆国、汪成民、李均之、赵玉林、钱复业、林云芳、曾小苹、杨巍然、强祖基、杜乐天、徐道一、任振球、邓志辉、石绍先、付虹、张元生、曾佐勋、汤懋苍、高晓清、陈维升、夏雅琴、李正心、胡辉、韩延本、沈宗丕、高建国、陈辉、岳中琦、池顺良、罗灼礼、刘德富、顾国华、查志远、马鸿钧、丁跃军、李志平、张建国、宋期等就地震成因与预测进行过不同程度的交流，结合笔者在多个地震区的调查和地震地质的研究，初步总结大陆地震不同演变阶段的前兆特征。

上述前兆异常是由多个震例（特别是大震）高度综合而成，而不是每个地震必然出现所有的异常，因此，前兆异常是复杂的，有待进一步深入研究。另外，还需要说明和探讨的是：

（1）异常范围大小与单个地震的震级和震群分布有关，也与中下地壳非均匀热结构有关。初步估计，对于里氏7级地震的地面异常而言，长期异常区域（如干旱）震中距可达400 km，中期异常区域震中距约300 km，短期异常区域震中距约200 km，临震异常区域震中距约100 km，宏观异常极强区可能在数十公里范围内。因此，根据中长期异常范围确定短临地震监测区，分析增量热能（灾能）对应的地震数和震级，然后分步骤开展立体监测和加密流动监测，根据前兆异常发展的时空强结构基本能够确定震中和发震时间，进行临震预测和预报。

（2）地震前兆异常具有时空和强度的有序结构，长中期预测阶段前兆异常从外围向震中定向迁移，进入短期和临震预测阶段前兆异常则从震中向外定向迁移；异常强度从热流体准备撞击到进行撞击之间出现物理量突变；流体相关异常早于固体相关异常；震级越大，前兆异常的时空尺度和强度变化越大，出现的前兆异常种类越多，前兆异常演变的规律性越强，地震预测的难度越小；显著的地应力相关突变异常，多与同震形变和位移有关。

（3）晚期小范围异常区出现在早期大范围异常区边缘，晚期异常种类多、频率高、分布集中的区域往往与震中有关。在地质构造上，表现为流动的“热河”与相对固定的“河岸”之间的接触部位或潜在的撞击点。

（4）尽管震前下地壳因增量热能引起半固态流变物质发生异常运动，由于中地壳吸收了大量的热能发生固态相变，顺层调节的韧脆性转换带储存和积累了由热能转换的应变能，因此，上地壳表及其浅表层次变形不明显，位移速率变化不大，地下至空中热流体变化较明显，局部甚至出现由热隆作用产生的抛射式滑坡、地陷、矿难等自然灾害；同震块体位移和地表破裂突然加著，热流体异常活动更加显著，热流体撞击后的上冲作用常导致抛射式碎屑状远程高速滑坡；震后地下水汽释放和局部大气环境改变，常出现异常降雨，伴生泥石流。

（5）不仅要重视热灾害链异常区立体监测网络体系的建设，还要大力加强新型群测群防网络体系建设，宏观异常信息对短临地震预测具有极其重要的作用。

导致现今地震预测困局的核心问题是流行的地震成因理论与客观的前兆异常不相容：大量的前兆异常用流行的弹性回跳模式无法解释，而根据弹性回跳模式应有的地应力、地形变异常在震前不一定出现或者不明显。

地震前兆异常之间是否存在本质上的关联性？它们到底指向哪种地震机理？广泛分布在中地壳的震源与刚性板块和脆性断层有关吗？地震的时、空、物、能有规律和联系吗？本文试图以鲁甸、景谷、康定地震预测为实例，解开这些谜团。

2鲁甸、景谷、康定地震的预测方法

2.1地震预测方法简评

目前地震预测方法主要建立在测震、形变、电磁、地下流体、地震地质的基础上。由于执行专群结合、群测群防路线，我国曾取得了成功预测海城地震这一开天辟地的伟绩。我国震灾惨烈，损失惊人，一些地震专家在困境中奋争，自创或尝试多种地震预测方法。在长达3年多的《国家自然灾害空间信息基础设施专项》（简称623国家专项）论证过程中，我们了解、调查和评价了许多地震预测方法，包括地倾斜法、地应力法、地应变法、钻孔应变观测法、GPS法、地震活动性图像分析法、旱震关系及五项气象指标法、地电法（包括HRT波法、地电阻率法、自然电位法、土地电法等）、地磁法（包括地磁转换函数法、地磁加卸载响应比法、地磁空间相关法、地磁异常极值环交汇法、磁喷法等）、电磁辐射法、MDCB法、重力法、天文—重力法、地声法、地温法、地下流体法（包括水温、水位、水硬度、水氡等）、地气法、次声波法、动物异常法、氢氧碳同位素法、卫星热红外图像地震预测法、卫星红外潜热通量法、卫星热红外相对功率谱法、卫星热红外长波辐射法、三性（周期性、倍周期性和黄金分割性）法、磁暴（月相）二倍法、加卸载响应比法、电离层扰动法、可公度法或信息有序法、三星一线法、天文周期分析法、天文活动法、时纬残差法、震象云法、天象—干支序列图法、日食—地震法等等。经筛选后分析了一些预测方法的科学性、准确性、技术可行性、发展潜力和存在问题。总体而言，有成功的案例，也存在误报、漏报现象。常常是异常出现，但背景不清、机理不明、关联不紧，致使灾种难辨、定位不准，精度不高。

归纳起来，地震预测方法大致可分为经验预测、统计预测、物理预测和综合预测。实际操作中，这些方法各有侧重，互有重叠。

目前，有人倾向将海城地震预测归为经验预测，从观察到的小震群异动和宏观异常预测主震。因此“小震闹、大震到”，成为海城经验。而这种模式在接下来的唐山地震中没有复制。但是，唐山地震之前出现了宏观异常组合，还创造了“青龙奇迹”，可谓青龙经验。

既然是经验，就应当积极总结，不断提高认识水平。不是所有大震和强震都有小震群作为前震；宏观异常的构造物理属性、地震专属性、机理关联性还不清楚；宏观异常与微观异常的系统监测和关联分析仍不够。如果解决了这些问题，科学地建设立体监测网络体系和新型群测群防体系，地震预测水平一定能够大幅提升。

统计预测是将一定时空结构中地震活动性的统计规律用于地震预测。目前西方国家倾向于研究特定区域的地震序列，建立地震概率模型，预测未来地震发生的概率。主要问题是：地震系统的时空结构是动态变化的；地震概率模型参数选取具有主观性；地震概率预测结果具有不确定性。

物理预测是从地震形成机理出发，监测地震发生之前物理量异常并进行地震预测。这是地震预测的发展方向，也是地震学家努力奋斗的重要目标。但是，目前地震机理不明；难以阐明常见的前兆异常与流行机理之间的关联性。因此，建立在断层活动、弹性回跳等物理机理之上的地震预报试验场接连受挫，滋生和助长了地震不可预测论。

综合预测是整合多种地震监测、预测手段，对前兆信息和预测建议进行综合判断。将理论与实践相结合，优选多种地震监测方法和手段，获取地下、地面、天基、空基地震异常监测数据，进行多类异常的四维动态立体填图，提取关联、耦合的前兆信息，有序开展长、中、短、临相结合的综合预测（李德威，2010，2011，2012；李德威等，2013a），是地震预测的发展方向。需要解决前兆是否有效可靠、方法是否合理权重、灾种鉴别是否有科学依据、整体与个体如何协调、人机能否有机结合等问题。

笔者认为，物理预测是地震预测的根本，综合预测是地震预测的根基。如果地震预测的思路由破裂物理预测转向热流物理综合预测，许多问题可能迎刃而解。芦山、鲁甸、景谷、康定地震的中期预测在一定程度上支持了上述设想。要深入认识这一预测思想，离不开下地壳层流、盆山耦合大陆动力学和多级物质循环地球系统动力学这个大视野。

2.2陆内地震的预测思路

初步总结陆内地震预测的思想、方法和步骤。这是一个以热流为主导、循序渐进的地震监测预测思路，目前只经过了芦山、鲁甸、景谷、康定等地震中长期地震预测的检验，还需要大量的验证和不断地完善，特别是在短临地震监测和预测方面。

对于短临地震预测，应当采取围震思路。在调查和研究热灾害链和“热河”精细结构及地震空区的基础之上，圈定潜在震区，在已有监测系统基础上及时建设针对热流体活动有关的立体监测体系和流动加密监测体系，逐步缩小地震发生区的时空范围，精准预测预报地震。本文强调必须立体监测与热流体灾变异常有关的前兆，在此只作综合分析和理论探讨，有待实践检验。因此，应当深入研究不同类型前兆之间的关联性、前兆与灾种之间的关联性、前兆异常与灾害机理之间的关联性，针对热异常关联前兆，建立专群结合、群测群防的立体监测系统和预测体系，能够不断提高地震预测水平。

由于活动“热河”中固态、半固态流变物质在重力、天体引潮力等因素作用下异常运动，撞击到强弱块体的分界面，形成地震。根据热流体撞击形成地震的模式，热和流体是地震的直接前兆，机理相关的气象异常、物理异常、化学异常、水文异常、地质异常、生物异常、天文异常是地震的间接前兆。地震监测和预测应当在热流体地震撞击机理指导下科学部署，对通过热灾害链圈定的短临地震预测区，除了运用优选的地震监测手段获得综合前兆异常，并进行异常动态立体填图外，还要在短期地震预测区布设一系列横穿“热河”边界过渡带和几条平行“热河”并靠近边界过渡带的高精度地震探测剖面，抓住S波对流体敏感的特点，重点动态监测震源层及其下部的热软化物质的异常运动，就像医学CT或核磁共振持续监测人体一样，能够准确把握“血管瘤”的破裂过程（热流体撞击产生地震的过程）。在临震期，合理布置按上述思路研发的灾害监测和航拍无人机、综合异常多功能监测车，开展短临地震预测区的流动加密监测和实时摄像，并自动传输到监测预测预报中心，这样就能够灵活提供急需区段的前兆信息，而且能够在第一时间发布地震灾情，绘制烈度图，及时开展科学救灾。

因此，地震预测应当是在符合客观事实的地球科学系统理论指导下采用科学方法而开展的一项系统工程。地震预报还需相应的科学决策与应急管理系统。

值得指出的是，地下剩余热能可形成多种灾害，火山、跨季度干旱、抛射式滑坡、天然矿难、天然森林大火、突发性雾霾等许多灾害可释放剩余热能和热流体，产生各种相似异常，因此灾种鉴别十分重要。高温部分熔融或高度熔融形成的岩浆由于低密度产生的浮力，只能向上运动，喷出地表形成火山，会出现区域相对固定的点式异常。而地震的形成机理不同于火山，是下地壳固态、半固态流变物质非均匀层流导致顺层滑脱的中地壳不断积累应变能，下地壳热能通过中地壳应变能转换成上地壳机械能之后发生脆性破裂，形成地震。当局部受热的上地壳经过冷加工硬化变脆后更易于破裂，因此，跨年度的地源热干旱之后出现异常降雨或冰冻，当地下热能再度异常积累，易于发生地震。在灾种鉴别和灾源定位过程中，必须以致灾机理为指导，对各种前兆异常进行系统的组合分析、时空分析、关联分析和主次（权重）分析。还要深入研究直接前兆与间接前兆的关系，区分必然前兆与偶然前兆，挖掘机理关联的前兆异常数据，分层块提取有效前兆异常信息，进行大数据定量处理和关联分析，定能科学地预测地震。

地球如人，热流似血。地震是中下地壳热河（“血管”）局部热流体异常运动击破较坚硬的围岩（“血管瘤破裂”）造成的，因此，对由热灾害链圈定的异动热河高危区段进行高分辨率多波段多分量四维地震探测（“CT动态扫描”），查明热河潜在发震区的精细结构、流动轨迹和流速变化，能够准确地预测地震。

查明异常热活动构造区带特定时空结构中自然灾害连发、群发、多发的致灾因子，辨认有效、可靠、确定、关联的前兆，搞清前兆异常和前兆组合所指向的致灾机理，是经验预测上升到科学预测的关键环节。纵观各大地震，地气耦合现象明显，震前闷热、震后降雨（或下雪）似乎形成规律。震前热流体异常与物理、化学、气象、水文、生物、天文等异常之间可以建立成因联系，但是用断层活动、弹性回跳难以合理解释。因此，从地质、自然灾害、地球物理、气象等事实出发，认识地震和关联灾害的形成机理，并以此作为创立大陆动力学和地球系统动力学理论体系的组成部分，应当是当代地球科学的重要发展方向。

2.3鲁甸、景谷、康定地震的预测思路

笔者在参与汶川地震之后45天的科技赈灾过程中，几乎天天下雨，不仅看到了地表破裂、喷水冒沙、抛射式碎屑状远程高速滑坡等地震相关现象，而且在多地采访了幸存者。综合各种宏观前兆异常，主要有：震前天气闷热；地下闷响；河井池塘的水位、水温、颜色、气味发生变化；动物异常；地震当天北川中学初二（1）班上物理课时张家春老师演示的指南针不停地乱摆；地震发生时看到拖着黑烟的火球，砂土像火山爆发或大炮发射一样喷出，闻到硫磺等异味。

笔者耳闻目睹地震的惨状和现象，下定了探索地震机理、开展地震预测的决心。在地球系统动力学多级循环假说、热灾害链思想和地震热流体撞击模式的指导下，终于对西南地区2013年以来发生的芦山、鲁甸、景谷、康定等强震均作了准确的长期和中期预测。

我们曾经阐述了芦山地震的预测过程（李德威等，2013a）。鲁甸、康定地震与芦山、玉树等地震属于同一个地震构造系统和发震系列（李德威等，2014）。因此，本文不再详细叙述鲁甸、景谷、康定地震的预测过程。概括其要点是：在地质构造上，它们沿着中新世以来从青藏高原下地壳非均匀流出产生的川滇菱形块体的东北部边缘构造带，与玉树地震一起，呈串珠状分布。在长期和中期异常上，它们经历了共同的跨年度干旱和后续的气候突变和波动期，与地源热能、水汽等作用有关，“热河”异常流动导致地气耦合，形成热灾害链。而且，在川滇菱形块体的东北部活动的边缘构造带“热河”流域出现了强震空区。由于我们没能掌握“热河”的精细结构，目前还不能十分精准地确定异动热流物质的可能撞击点，这种预测尚处于理论分析阶段。

地震长期和中期预测的基础是热灾害链。初步认为，地震是开放的复杂地球系统中热能驱动不同相态的物质发生多级循环运动的表现形式之一，与其他自然灾害一起，构成有序的时—空—物—能结构。地球四维非均匀热结构导致热流体非均匀流动，不仅成矿成藏，而且致灾致难，并决定了自然灾害的区域性、群发性、连发性、关联性、有序性、迁移性、突发性等特点。因此，海洋之中、大陆内部、洋陆之间、地气之间、天地之间的各种自然灾害常存在着关联性。地球系统在非均匀热动力作用下产生的致灾因子构成时空有序、机理关联的灾害结构，称为热灾害链（李德威，2011，2012，2014a，2014b）。

2004年以来，西南地区自然灾害频发（陈桂凡等，2013）。震前3～6年出现跨年度干旱，伴生森林大火（长期前兆异常）；震前1～2年常表现为异常降雨和气候大幅波动，甚至出现大面积冰冻（中期前兆异常）。因此，地源热流体异常活动及其地气耦合是地震中长期预测的基础，也是预测极端气候事件的一个新思路。

自然灾害连发、群发区域与地热异常区域吻合。在特定地热构造单元热能的源—汇—释量级和过程决定了热灾害链的强度和时空结构。大陆盆山体系中下地壳具有四维非均匀流动特征。地貌、地质、地球物理、GPS、地震、温泉、矿床的综合特征表明（孙洁等，1989；曾融生和孙为国，1992；Royden et al.，1997；Clark and Boyden，2000；Wang et al.，2001；周伏洪等，2002；蔡学林等，2003；李德威，2008b；罗文行等，2008；马宏生等，2008；Burchfiel et al.，2008；Hubbard and Shaw，2009；Wang et al.，2009；Bai et al.，2010；Zhang et al.，2010；郑勇等，2013），青藏高原地壳分层流变显著，深部地壳物质向北流动受阻后转向东流，盆山边界是地貌、地壳厚度、岩性强度的突变带，也是地震多发带。近年来西南地区灾害频发与源于恒河盆地流经亚东、尼木、羊八井、错那、雁石坪、沱沱河、玉树、鲜水河、安宁河、小江的下地壳“热河”中下游异常流动有关，有关内容下面将作较详细的说明。该区段上一次地震活跃期从1970年1月5日号通海7.8级地震开始，其后相继发生了1973年2月6日炉霍7.6级地震、1974年5月11日云南永善7.1级地震、1976年11月7日和12 月13日盐源6.7级和6.4级地震、1979年3月15日普洱6.8级地震和1981年1月24日道孚6.9级地震，地震烈度等震线长轴为NW—SE向或近南北向，与“热河”流动方向和热液体撞击形式有关。

上述地震大爆发事件之后，经过了23年的灾害平静期，于2004年开始出现灾害群发和连发事件，根据其时空结构，可分解成两个次级的热灾害链。早期的热灾害链偏北，晚期的热灾害链偏南。

总之，根据热灾害链的时空结构，结合异常活动的“热河”流域灾能（剩余热能）聚散规律及其地震空区分析，对芦山地震和鲁甸、景谷、康定地震进行了较准确的长期和中期预测。

3鲁甸、景谷、康定地震的预测原理

3.1流行的地震机理及其存在的问题

自从1901年Reid提出弹性回跳模式以来，人们相信断层运动产生地震。一般认为，大陆构造地震是在板块运动构造背景下，水平挤压力远程作用于活断层，导致活断层闭锁段应力达到强度极限时发生破裂，断层两盘滑动位移，出现反向的弹性弯曲，断面破裂处产生冲击卸载，激发的弹性波向外传播，断层两盘变形岩石弹性回跳并发生位移。

这个模式的核心是弹性体内发生位错和位移。因此，人们努力监测震前断层两盘弹性断块“应有”的显著变形和断裂带“应有”地应力集中。按照此思路一些多震发达国家精心布置了高密度、高精度的监测网络，至今没能成功预测地震。特别是受到美国Parkfield地震预测试验场的严重打击之后，地震预测的风向突变，地震不可预测成为主流观点。

地震真的是弹性体位错而成吗？若如此，震前断层两盘应有十分显著的弹性变形，为何用高精度的GPS等手段测量不到？人们是否忽视了真正的震源物理和震源之下的能量—物质状态而过分强调表浅层次直接观察到的同震脆性破裂现象？震源应当是认识地震成因的关键因素，地震体的物理状态、构造性质、能量来源、物质运动是绝对不能被忽视的。因此，震源孕育过程和灾能聚散规律是解析地震机理的关键。

3.2热流体地震机理简介

震源体的发育背景、基本结构、物质组成、物理性质、构造样式，是认识地震成因、开展地震监测和预测的基础。震源体周围震前常区域性面状分布与热有关的异常，傅承义等（1979）将这种现象称为“红肿说”，认为断层理论阻碍了地震机理研究和地震预测。正如滕吉文（2010）指出：“一个强烈地震的孕育发生和发展，必须具有一个特异的深部介质和构造环境，及其在力源作用下物质的重新分异调整和运移的深层动力过程。”陈胜早（2006）认为震源一般发育在300～400℃的地壳脆韧性转换带。

震源体物质没有直接喷出地表，肉眼无法观察，深钻遥不可及，难以确定其精细结构。然而，根据不同构造单元、构造层次、构造阶段展现的地质现象总结的地质规律，可建立包括地震甚至灾害链的构造模式，阐明地球不同层块之间能量转换和物质运动的基本规律。地球物理探测、高温高压实验、数字模拟也是认识震源物理、震源机理的重要途径。

笔者自从1990年以来对青藏高原及邻区持续开展地质调查和研究，先后创立了下地壳层流驱动盆山耦合的大陆动力学假说、软流圈层流驱动洋陆耦合并制约大陆盆山耦合的地球内部系统动力学假说、多级物质循环关联运动的地球系统动力学假说、热流体撞击地震成因模式和热灾害链成因模式。这些地学基础理论和应用基础理论的新尝试是芦山、鲁甸、景谷、康定等地震预测的基础，可为大陆地震预测提供理论指导。

概括起来，简述如下：地球大地构造单元与地貌单元吻合，分为一级的大洋与大陆和二级的造山带与盆地，同级构造单元结构互补、构造转换、建造互换、机理统一、演化同步。两级构造单元之间具有结构、构造、建造和机理的自相似性，并具有关联成因：地球内核在高温熔融的外核中极不稳定而偏离地心（黄定华等，2001），由此引起外核巨量流体的顺层流动，外核的非均匀极高温热流体在核幔边界汇流形成巨量地幔上升热流。不同形态的低密度深地幔热流物质在浮力作用下呈墙状、柱状上涌。墙状深地幔热流物质上涌（地幔墙）构成热线，形成洋中脊；柱状深地幔热流物质上涌（地幔柱）构成热点，形成洋岛。地幔墙热流物质上涌造成上部地幔部分熔融形成软流圈，软流圈溢出低密度的玄武岩形成洋中脊，密度相对较高的地幔半固态热流物质随着不断倾斜的软流圈底面从洋中脊顺层流向大陆，带动仍具有陆壳成分的洋盆水平扩张。流入大陆的地幔软流圈低密度物质造成大陆垂向生长，陆洋差异升降，大陆剥蚀物主要通过河流搬运到海洋。大陆软流圈加厚部位发生线性或点状地幔底辟作用，一方面改变了莫霍面的产状，另一方面引起大陆下地壳部分熔融和固态流变，下地壳固态流变物质侧向流失的区域形成大陆裂谷和沉积盆地，下地壳热流物质侧向流入的区域形成造山带或高原，山脉剥蚀物充填到盆地中，盆山之间地壳物质发生循环运动。

地球内部时空分布不均匀的热能经过一定时间积累之后达到临界状态，在天文事件的触发下，过剩热能快速释放产生不同尺度的构造活动及其灾害—环境效应，一系列关联的突发性灾害连发和群发，引起不同级别的气候环境变化。对于导致恐龙灭绝这样的巨型热灾害链（李德威，2014a，2014b）与洋陆系统热构造强烈活动有关，地下热能经过长时间积累到一定程度，岩浆侵位形成大陆超级干旱，海洋出现热喷泉，洋陆同步发生超级火山喷发和岩体侵位，巨量温室气体溢出，海水大幅度升温和缺氧，风暴潮频发，海洋生物灭绝。陆地火山发生碎屑流和岩熔流，引发森林草原大火，火山和植物燃烧释放巨量二氧化碳，导致地球内部释放出的热量不能向外层空间扩散，地球气温升高，氧气含量急剧减少；火山喷发释放大量的二氧化硫，与大气层中水汽混合形成酸雨。随着火山岩浆作用持续进行，地壳释放大量热能后逐渐变冷，同时火山灰、二氧化硫和燃烧产生的灰烬遍布地球上空，反射太阳辐射热能，进一步造成地球浅表层变冷，进入冰期。同时大批陆生植物死亡，在低温环境下易于保存，形成煤层。随着地球内部热能再积累，在地幔和下地壳热流作用下冷而脆的大陆上地壳发生公里级的断裂，洋陆相互作用区发生更大规模的变形，形成无数群发和连发的里氏12级以上地震，部分地震引起海啸，热喷泉、滑坡、岩崩和泥石流。通过构造活动、地震、热喷泉等释放的热能、水汽改造冷凝的大气混合物，形成大范围巨量降雨，超级洪水强烈冲刷和剥蚀山脉，盆地同步形成巨厚（km 级）的砾石层不整合覆盖在下伏地层之上，陆生动物最终灭绝。这种巨型热灾害链能够彻底改变地球环境，可能是海洋生物、陆生植物和陆生动物锐减甚至分批灭绝的根本原因。

西南地区正在发生微型热灾害链，其时空结构、演变过程和形成机理基本类似于巨型、大型、中型和小型热灾害链（李德威，2014b），但是其孕育背景、热能量级、致灾作用显著不同。陆内微型热灾害链与地幔墙、地幔柱底辟和软流圈层流没有直接的关系，而且热能积累时间短，到不了洋陆耦合、盆山耦合峰期大规模熔融岩浆垂向上流的程度，而是地壳局部非均匀固态流变物质侧向异常流动产生的灾害效应。

对于陆内地震而言，大陆非均匀的物性结构和热结构决定下地壳流动具有非均匀时空结构。在下地壳流层中出现热流体线性富集带（称为“热河”），当“热河”中热流体在天体触发或因自身流量、流速、流向突变而发生异常流动，以俯冲、上冲、仰冲、侧冲和斜冲的方式撞击周边坚硬围岩，受到爆轰的中地壳形成点状震源，震源附近热能转化为机械能，通过上地壳脆性断层活动和地表破裂释放出来，形成地震。上地壳破裂系统加剧了包括水汽在内的地下气体的释放，改变了局部的大气格局，形成异常天气。

根据上述假说或模式，震前不同阶段会产生一系列与热流体异常有关的前兆，以此为基础建立地震及其关联灾害的动态立体监测与预测系统，不仅可作为地震及其关联灾害预测的一个新思路和新途径，而且是检验地学基础理论与基础应用理论新体系的重要标准。然而，要创立系统的地学理论体系，有赖于大量事实的归纳和新事实、新预测的反复检验，科学假说才有可能上升到科学理论，更好地指导实践。

地震预测和地震机理的复杂性在于震源通常出现在地下10 km之下的中地壳，人们只看到地震发生后的地表破裂，易于片面强调断层的作用。虽然我们无法直接观察到现今的震源物质，直接测量震源物理量，但是大量的地质和地球物理现象、高温高压实验、数值模拟和地震前兆表明：垂向上，地震常发生在脆性上地壳与韧性下地壳之间的韧脆性转换带；横向上，地震常发生在地壳厚度和岩石强度突变带；物质上，沿着下地壳四维不均匀流变带（“热河”）的半固态、固态物质点式作用于具有显著强度差异的分层分块物性界面；构造上，下地壳韧性流变与中地壳韧—脆性转换和上地壳脆性破裂密切关联；能量上，震源之下高能域的地热能通过震源层应变能积累与转换制约震源之上低能域的机械能。高分辨率地震层析成像可获得震源区精细的速度结构，一些大地震的震源位于低速—高导层与高速—高阻层之间，与韧脆性转换及其热流体非均匀作用有关。

需要强调的是，地震机理必须遵循热力学第二定律，发育在大陆地壳分层流变构造系统韧—脆性转换带的地震能量，只能来自更高能的下部韧性区域，不可能来自更低能的上部脆性区域。若果真如此，必然出现重大的观念转变：断裂是表，热流是因，地震是果。仅就地震与断裂而言，前者是因，后者是果。如今因果关系倒置，因此得出地震无法预测的结论。

在大陆地壳非均匀分层流变构造系统中，局部热能集中区带的下地壳韧性流层中，半固态、固态流变物质持续缓慢地流动产生活动的顺层韧性剪切带，通过中地壳韧—脆性滑脱层的变形调节，实现热能—应变能—机械能转换，导致上地壳断裂和地表层破裂。震前地壳表浅层次构造变形及其增量应变并不明显，可能出现微量的地倾斜、地面升降，有时产生地表塌陷和抛射式滑坡。上地壳及表浅层突变性的破裂、位移只能是同震构造。地壳深部热流体的异常流动，会在相应的时空结构中出现热、甲烷、二氧化碳、氡、地下水、电性、磁性、重力、波速、地声、气象、某些元素的同位素、水文、天文观测、生物等异常，它们是地震预测的基础。从热流体撞击地震模式出发，应当加强短临地震监测区多层次地源异常热能及其相关的尾波Qc值、S波分裂、泊松比等的动态监测和综合分析。例如，Aki（2004）研究圣安德烈斯带地震与尾波Qc值之间存在时间关联性，认识到地壳韧性层与上地壳脆性层之间发生相互作用，地震的动力来源可能是脆—韧性转化带。

3.3鲁甸、景谷、康定地震的机理分析

上述从地球系统到地震机理的系列创新思想是认识汶川、玉树、芦山、鲁甸、景谷、康定等地震机理的基础。限于篇幅，除去细节，阐述如下要点：

（1）青藏高原东南部地震孕育的构造背景。大陆一般经历基底演化、洋陆演化和陆内盆山演化三阶段地质构造过程。青藏高原及邻区经历了太古宙—古元古代结晶基底和中—新元古代褶皱基底形成阶段、原特提斯—古特提斯—中特提斯—新特提斯及其相关古陆向南有序迁移阶段和青藏高原北东部燕山期、中部喜马拉雅早期和南部喜马拉雅晚期的盆山耦合阶段和3.6 Ma以来青藏高原整体均衡隆升阶段（李德威，2008b；Li，2010，2013）。青藏高原东南部地震活动与燕山期以来的陆内构造过程有关：燕山早期带动太平洋生成和扩张的地幔软流圈物质层流至四川盆地聚集加厚，底辟上升，四川盆地下地壳热软化物质从中心向外层流，周边山脉下地壳加厚并同步隆升，龙门山造山带同期形成。燕山晚期以来，四川盆地逐渐冷却，基底固结硬化。喜马拉雅晚期（中新世）印度洋地幔软流圈物质层流至恒河盆地加厚并底辟上升，恒河盆地下地壳热流物质向北不均匀流向青藏高原，喜马拉雅造山带同步形成。由于源于恒河盆地的下地壳层流仍在进行，下地壳不均匀流动在青藏高原及邻区地壳加厚过程中形成三角形发震构造域（李德威等，2014）。从地质、地震、地球物理综合分析，三角形发震构造域内部存在7条源于恒河盆地向北流动的下地壳“热河”（李德威等，2009，2013a；李德威，2010），驱动多震域向北漂移，从东到西为：（1）错那—桑日—墨竹工卡—嘉黎—波密—察隅—保山—耿马“热河”；（2）亚东—尼木—羊八井—错那—雁石坪—沱沱河—玉树—鲜水河—安宁河—小江“热河”；（3）定结—谢通门—申扎—双湖—西金乌兰湖—库赛湖—冬给错拉湖—迭部—舟曲及平武“热河”；（4）老定日—锁作—许如错—当惹雍错—依布茶卡—朝阳湖“热河”；（5）浪强错—打加错—扎日南木错—扎西错“热河”；（6）亚热—森里错—仓木错—美马错—清澈湖构“热河”；（7）普兰—扎达—扎西岗—印度河“热河”。这些“热河”上游的宽度、流量、能量较大，下游宽度、流量和能量较小，但是支流和撞击点多，因而在中新世以来形成的青藏高原仍具活性的中下地壳与燕山期形成现已固结的塔里木、柴达木、四川、鄂尔多斯等盆地基底之间构造边界带易于发震。李德威等（2009）首次提出下地壳存在“热河”，主要依据强震沿着高热流异常带和磁异常带的地堑（裂谷）分布，也是中新世以来金属矿床成矿带的集中部位，地壳尺度的线性热隆伸展构造显著。大地电磁测深获得的青藏高原东南部电性结构（Bai et al.，2010）证实两条与推断的“热河”吻合的下地壳非均匀流动通道。亚东—尼木—羊八井—错那—雁石坪—沱沱河—玉树—鲜水河—安宁河—小江“热河”中下游呈串珠状分布的玉树、汶川、芦山、鲁甸、景谷、康定地震正是孕育在这种地壳热活动构造背景中。应当一体化开展“热河体系”的地质—地球物理—地热—地震及其相关灾害—矿床—地壳（工程）稳定性—环境变化的理论研究、系统勘查、效益评估和综合利用。

（2）垂向上的强弱边界与多震层。近20年来，确认大陆地壳分层流变是超越刚性岩石圈板块的突破口之一。陆内地震不是脆性断层造成的，而是流固耦合的结果。大陆下地壳因异常热活动而富含流体的韧性层抗剪能力弱，由于盆山体系莫霍面倾斜，在重力、剪切力、惯性等作用下半固态—固态流变物质易于发生顺层流动。非均匀流动的下地壳与其上的固体壳层之间通过顺层韧—脆性剪切带与上地壳脆性层调节，因此，具有顺层滑脱性质的中地壳韧—脆性剪切带才具有吸热储能和热能—应变能转换功能，当积蓄的地震能量达到边界条件或外界触发而发生地震。过量热能导致热流体异常运动，在非均匀活动地壳内部和边缘多点撞击构成震源层。脆性上地壳先期产生的断层带是震前热流体良好的外溢通道，也是发震时优先但不是唯一的灾能突然释放区。

统计分析表明，青藏高原地震震源集中分布在15～40 km的深度范围，其中30～33 km深度是一个优势层（罗文行等，2008，2012）。南北构造带南段多震层主要分布在地下10～15 km，33 km附近地震也较多。壳内多震层与其下的低速层、低阻层密切相关，青藏高原东部下地壳存在向东非均匀运动的高温富含流体的部分熔融层，与刚性的四川盆地形成鲜明反差，中下地壳波速、电性的各向异性与地表GPS的动向耦合特征反映了下地壳非均匀流动的方向。

（3）横向上的强弱边界与撞击点。活动的大陆下地壳存在四维非均匀流变特性，不仅垂向分层流变，而且横向分层流变，还存在时间上的动态演变。下地壳“热河”与两侧的固体岩块之间存在流固作用，“热河”的产状、形态及其内部固态、半固态流变物质的温度、流体类型与含量、黏滞性等决定了异常流动状态和撞击时的入射角和撞击点位。从热流物质异常流动到爆轰围岩，由连续变形运动到加速突变运动，在地壳、土壤圈、生物圈、水圈、大气圈都有不同程度的动态响应。地声、火球、大气甲烷含量的突变、碎屑状抛射滑坡等说明中上地壳在富含甲烷、二氧化碳等热流体的撞击下可能伴随甲烷爆炸。瞬时高压冲击波改变了震源体及周边岩石的力学状态，甚至发生构造和岩石重组，如假玄武玻璃的形成。主流线方向、强弱边界产状、热流体入射角与撞击点的关系决定了震源机理及其固体上地壳弹性波（纵波、横波、表面波）的传播方向、路径和范围，影响地震烈度分区。烈度图等震线长轴方向一般平行于强弱边界。主流线与强弱边界平行形成走滑型（或侧冲式）地震，玉树、鲁甸、康定地震位于上述“热河”的东北侧，应属于此类。主流线与强弱边界直交形成逆冲型（包括仰冲式和俯冲式）地震，芦山地震属于逆冲型（或仰冲式）地震。主流线与强弱边界斜交形成逆冲兼走滑型（或斜冲式）地震，汶川地震属于此类，而且沿着龙门山从映秀主震区向北东方向，热流体撞击方向与强弱边界之间的交角变小，走滑分量增大。

郑勇等（2013）利用地震台网数据采用层析成像方法和远震接收函数分析方法，获得了南北构造带中段精细的S波速度结构、地壳厚度和泊松比分布状况，显示汶川、芦山地震均位于地壳厚度和波速结构剧烈变化梯度带，地震深度处于从均匀波速结构向非均匀波速结构的过渡带。鲁甸、景谷、康定地震产于与汶川、玉树、芦山地震类似并关联的地壳分层流变构造环境，但是横向流变结构有所不同，青藏高原活动弱地壳与四川盆地热耗散固结硬化的强地壳之间的构造边界在雅安至石棉一带发生转折，燕山期四川盆地与周边造山带同步盆山耦合，盆地NW侧的盆山边界呈NE—SW走向，SW侧盆山边界为NW—SE走向。中新世以来，从恒河盆地向北不均匀流入青藏高原的下地壳热流物质，受到塔里木、柴达木盆地固结地壳的阻挡后转向东流；沿着“热河”支流大角度流向四川盆地的热流物质除了造成龙门山造山带地壳活化、伴生仰冲型地震外，大量的热流物质受阻后转向分流，一部分汇入东昆仑—西秦岭构造结外泄，沿着固结的鄂尔多斯盆地南部、东部和西部边界流动，形成线性热隆伸展地震活动带；另一部分汇入带动川滇菱形块体向SSE方向运动的鲜水河—安宁河—小江“热河”，“热河”流向与两侧强度差异显著的构造边界近于平行或呈小角度相交，在“热河”转弯或“河岸”突出部位易于发生侧冲型或斜冲型地震。鲁甸、康定地震主要受NNW—SSE走向的鲜水河—小江下地壳“热河”非均匀流动过程中半固态流变物质侧向撞击强弱构造边界，物质流变方向与等震线长轴平行，类似于玉树地震，为走滑型（或侧冲式）地震。景谷地震位于近南北向小江下地壳“热河”的南缘，处于NW—SE走向保山—耿马“热河”的结合部位，可能与小江下地壳“热河”西侧的侧冲作用及其对保山—耿马“热河”北侧硬块的撞击作用有关。

（4）灾能（剩余热能）的源、汇、释（泄）与地震平静期—活跃期转换。从本质上说，大陆下地壳四维非均匀流动是热能源、汇、释动态平衡过程。对于汶川、玉树、芦山、鲁甸、景谷、康定等地震和青藏高原的地震而言，下地壳是灾能源区，增量热能从恒河盆地通过非均匀下地壳流动实现远程输送；中地壳是灾能汇区，传导至中地壳的下地壳热能导致固态相变，转变成应变能；上地壳是灾能释区，当中地壳应变能积累到临界状态，在下地壳“热河”加速运动作用下，所积累的应变能通过上地壳断裂活动和地表破裂释放，发生系列地震；对于青藏高原“热河”而言，地幔软流圈层流和底辟是源区，恒河盆地下地壳是汇区，青藏高原及邻区的“热河”是泄区。还应当从时间上认识热能收支与灾害演变，对鲁甸、康定等地震有关“热河”中下游而言，1981 —2004年是灾害能量孕育期，不断积累增量热能，源、汇难分；2004年至今，通过跨年度干旱、潜热及水气蒸发、排气、地震、抛射式滑坡、天然矿难等灾害活动不断释放剩余热能，自从汶川地震以来总体上震级呈递减之势，剩余能量释放强度减弱。

4讨论与结论

本文归纳如下5点仍有待深入探讨的重要认识：

（1）地震能够预测。理论上，地震是一种构造现象，构造是一种地质现象，地质是一种地球现象。地球如人，热流似血。具有生命力的地球是在热能作用下由热流体非均匀运动产生各种关联现象，在开放的复杂地球系统中展现出能量—物质动态演变规律。因此，从地球系统动力学角度认识洋陆及其盆山构造与建造的物质运动规律和能量转换过程，解析地震及其相关自然灾害的机理，能够认清地球多尺度动态演变规律，也是科学预测地震的指导思想和理论保证。实践上，我国在20世纪70年代中期就成功地预报了海城地震和松潘地震。近年来笔者对芦山、鲁甸、景谷、康定等地震作了长期和中期预测，说明地震是有规律可循的。我们不能只盯着震后都很难确定、更不知道是否连接到温度大于300℃的震源区的发震断层，也不必无休止地争论地震能否预测。针对我国目前仍然十分严峻的震情，当务之急是联合跨部门、多学科的专家切实开展地震科学预测，建立新型专群结合的群测群防体系，再创地震预测预报的辉煌。今后工作重心应当前移，在热异常多灾区实施取热减灾减排宏伟工程。

（2）在创新地学思想指导下中长期预测了芦山、鲁甸、景谷、康定等地震。1992年以来，笔者提出了下地壳层流驱动盆山耦合的大陆动力学假说、软流圈层流驱动洋陆耦合并制约盆山耦合的地球内部系统动力学假说、地核和太阳热能共同驱动地球多级物质循环的地球系统动力学假说、热流体撞击地震成因模式和热灾害链成因模式，并在此基础上总结了地震预测的思路、方法和步骤，分析了川滇菱形块体东北部边界构造带2008年以来发生的串珠状地震及其相关自然灾害与亚东—羊八井—玉树—鲜水河—小江下地壳“热河”异常运动的关系。在热构造孕震环境下，根据热灾害链的时空结构，结合灾能（剩余热能）的源—汇—释（泄）及其关联的“热河”地震空区分析，进行地震的长期和中期预测。芦山、鲁甸、景谷、康定等地震的长期和中期预测就是按这一思路进行的。

（3）初步阐述了陆内地震热流物理综合预测的原理、方法和步骤。按照多级物质循环地球系统动力学假说，洋陆之间、盆山之间、流固之间、地气之间、天地之间都存在物质和能量的有序交换。不同温度的热作用导致地球固—液—气相物质转换，熔融的低密度岩浆在浮力作用下向上流动不仅形成侵入岩，而且导致火山活动；固态、半固态流变物质在重力作用下侧向流动，形成顺层韧性剪切带，大陆下地壳顺层韧性剪切带与上地壳脆性破裂带之间的韧—脆性转换带是地震能量积累和转换区，也是震源层。因此，不同级别的热灾害链会出现热流体异常及其各种关联异常有序的时空结构，根据热灾害链时空结构及其能量聚散规律，可开展长中期地震预测；根据热流体前兆异常及其关联前兆异常时空强结构，可进行短临地震预测。如果在中期预测圈定的短临地震预测区开展高精度动态监测获得精细连续的地震波速结构，并在临震期实施地空联合流动加密监测并进行灾情速报，将大幅提高地震预测水平，取得巨大的减灾实效。

（4）根据热灾害链的时空结构，进一步分析我国西南和华北的地震形势。笔者初步认为，西南地区正在发震的热流通道中30多年积累的灾能大部分已经释放，从汶川地震到康定地震，震级有减小的趋势，未来3年内向南漂动的川滇菱形块体东北部边缘带“热河”流域还可能发生2～3次里氏6.5±0.5级地震。此外，还要重视川滇菱形块体东南侧错那—桑日—墨竹工卡—嘉黎—波密—察隅—保山—耿马“热河”的异动，未来的强震群发带将转移到东构造结及其川滇菱形块体西南部边缘带“热河”流域。近期需要引起极度重视的是华北的震情，特别是东北的震情危急，热灾害链时空结构明显（李德威等，2013b），近期气象异常波动，已经演变到短期地震预测阶段。因此，在中长期预测区（李德威等，2013b）建立新型专群结合、群测群防地震预测体系，开展以热流体前兆异常为主线的立体监测和短临综合物理预测，已经刻不容缓！

（5）取热减灾减排，是防灾减灾的根本途径。地震不是弹性体内的位错，而是分层分块分时非均匀流变系统中的热流体撞击。正如人体血液流动一样，地球热流体也发生非均匀流动，大陆下地壳“热河”是具有剩余热能的高温系统，热以传导、对流、辐射的方式向低能的中上地壳和大气层传播，经过一定时间积累的增量热能以灾害链的方式释放。因此，热能转变为应变能、机械能的自然过程是成灾事件，热能转变为机械能及其电能的人为过程则是造福工程。我们应当充分地认识到主动取出剩余热能变宝、被动释放剩余热能成灾的重要性，从“热河”中系统开发地热能，特别是从震源层之上的干热岩中抽取产生地震的剩余热能，在“热河”关键部位建设大型梯级地热发电站，必定能够从根本上减灾减排，快速走上可持续发展之路，推动人类与自然的协调发展。♦

【作者单位：(中国地质大学（武汉）地球科学学院，中国地质大学（武汉）重大地质灾害研究中心】

（摘自《大地构造与成矿学》2015年1期）

·高被引论文摘要·

被引频次：117

地震预测新途径的探索

尹祥础

为了寻求解决地震预报问题的新途径，本文分析了各种现有主要地震前兆。虽然它们在地震预报工作中曾经起过，并且今后还将起着重要的作用，但是，毋庸讳言，现有的各种前兆与地震的发生都不存在唯一性的对应关系。这表明，到目前为止，还没有找出这些地震前兆和地震孕育过程之间真正内在的本质联系。本文从岩石的固有特征及稳定性问题的本质出发，借鉴经典力学中的量纲分析与现代信息论的概念，提出了一个表征地壳力学稳定性的新的无量纲参数Fy，Fy=1-Mp/Mo，其中Mp，为地壳材料加载时的广义切线模量，而Mo是其弹性模量。所以Fy表征着地壳材料在同一载荷增量条件下，在加载与卸载过程中力学响应的比值。这一参数正好反映了稳定问题的本质。最后讨论了参数Fy的测量。地壳介质的加载与卸载可以通过固体潮来实现。由此引起的力学响应（广义位移）可以借助于多种现代技术进行测量，例如地面激光测距，SLR，VLBI和GPS等。

地震预测；地震前兆；地壳介质；地震孕育过程；力学响应；地震学家；无量纲参数；力学稳定性；切线模量；量纲分析

来源出版物：中国地震, 1987, 3(1): 1-7

被引频次：85

利用P波和S波的初动和振幅比计算中小地震的震源机制解

刘杰，郑斯华，康英，等

摘要：应用Snoke最新发展的利用P波、SV波、SH波的初动和振幅比联合计算地震震源机制解的程序，使用一些省地震局数字台网得到的数字波形资料，尝试计算了中小地震的震源机制解，并与原来仅用P波初动得到的结果进行对比研究。结果表明，这种方法是可行的。该方法为数字地震资料在地震预测中的应用提供了一种新的可能。

关键词：震源机制；初动；振幅比

来源出版物：地震, 2004, 24(1): 19-26

被引频次：73

孕震断层的多锁固段脆性破裂机制与地震预测新方法的探索

秦四清，徐锡伟，胡平，等

摘要：通过建立孕震断层中锁固段的临界破坏—断裂力学与其蠕变过程联系，发现每个锁固段断裂点与加速应变能释放起点的累积Benioff应变之比，有着依赖于锁固段个数的指数律关系，指数函数的底数为1.48。在最后一个锁固段破坏以前，每一个锁固段的临界破坏断裂过程表示一个中等强度或更大震级预震（preshock）的产生过程；在最后一个锁固段破坏之后，主震将发生。通过相当数量地震实例的回溯性预测分析，表明这一指数律能可靠地应用于中等强度预震和主震的预测。该方法能够在地震的中期、短期以及临震预报方面发挥作用。常数1.48很可能是表征中等强度或更大震级的预震和主震产生过程的普适常数。我们还对地震预报研究中某些有争议的问题进行了讨论，从理论和技术层面论证了本方法的适用性。

关键词：地震；断层锁固段；Benioff应变；常数

来源出版物：地球物理学报, 2010, 53(4): 1001-1014

被引频次：71

有关活断层分段的一些问题

丁国瑜

摘要：在地震预测、潜在震源的确定和重大工程的地震危险性评估中，活断层的的分段已成为日益引起人们重视的一个问题。但目前有许多涉及分段问题的概念和认识还有待深入的讨论。本文就活断层分段问题中有关不同类型的分段、段及分段的定义和理解、段的稳定性、分段尺度与段的层次结构、分段标志以及表部与深部破裂段的关系等六个方面的问题提出了看法与讨论。

关键词：段；分段；段的稳定性段的层次结构；段的标志

来源出版物：中国地震, 1992, 8(2): 1-10

被引频次：67

地震地层压力预测

云美厚

摘要：目前，地震地层压力预测方法归纳起来可以分为图解法和公式计算法两大类10余种。本文对各种地震地层压力预测方法进行了系统地归纳和总结，并对各种方法的特点、适用性以及存在的问题进行分析和讨论。在此基础上，就如何提高压力预测的精度，提出了一种简单适用的改进措施，经JLK地区的实测资料的验证，效果良好。

关键词：地层压力；地震预测；正常压实；异常压实

来源出版物：石油地球物理勘探, 1996, 31(4): 575-586

被引频次：66

地震预测——进展、困难与前景

陈运泰

摘要：地震预测是一个既紧迫要求予以回答、又需要通过长期探索方能解决的地球科学难题。20世纪60年代以来，中期与长期地震预测研究取得了一些有意义的进展，如板块边界大“地震空区”的确认、“应力影区”、地震活动性图像、图像识别以及由美国帕克菲尔德地震预报实践获得的正反两方面的经验等。但是地震预测尚处于初期的科学探索阶段，目前总体水平仍然不高，特别是短期与临震预测的水平与社会需求相距甚远。地震预测的进展主要受到地球内部的“不可人性”、大地震的“非频发性”以及地震物理过程复杂性等困难的制约。依靠科技进步，强化对地震及其前兆的观测，选准地点、开展并坚持以地震预测试验场为重要方式的地震预测科学试验，系统地进行基础性的对地球内部及对地震的观测、探测与研究，坚持不懈，对实现地震预测的前景是可以审慎地乐观的。

关键词：地震预测；地震预报；地震前兆；地球内部

来源出版物：地震地磁观测与研究, 2007, 28(2): 1-24

被引频次：57

地震预测：回顾与展望

陈运泰

摘要：本文概要回顾自20世纪60年代以来国际地震预测研究与地震预报实践的进展情况，指出地震预测这一既紧迫要求予以回答、又需要通过长期探索方能解决的地球科学难题目前尚处于初期的科学探索阶段，虽然总体水平仍然不高，特别是短期与临震预测的水平与社会需求相距甚远，但是近半个世纪以来并非毫无进展.文中以板块边界大“地震空区”的确认、“应力影区”、地震活动性图像、图像识别等方法以及美国帕克菲尔德（Parkfield）的地震预报实践为例，说明在中期与长期地震预测方面，地震预测研究均取得了一些有意义的进展。文中分析了地震预测在科学上面临的困难，阐述了为解决这些困难所应当采取的科学途径，展望了地震预测的前景，指出地震预测的进展主要受到地球内部的“不可入性”、大地震的“非频发性”以及地震物理过程的复杂性等困难的制约；地震预测虽然困难，但并不是不可能的;依靠科技进步，强化对地震及其前兆的观测，选准地点、开展并坚持以地震预测试验场为重要方式的地震预测科学试验，坚持不懈地、系统地进行基础性的对地球内部及对地震震源区的观测、探测与研究，对实现地震预测的前景是可以审慎地乐观的。

关键词：地震预测；地震预报；地震前兆；地球内部；地震震源

来源出版物：中国科学：D辑, 2009(12): 1633-1658

被引频次：53

电磁卫星和地震预测

赵国泽，陈小斌，蔡军涛

摘要：地震前几周、几天、几小时甚至几分钟，电磁场的强度、相位或者谱密度等将发生异常变化，电离层的电子、离子浓度和温度等会发生异常扰动，这些现象的出现向人们发出地震即将来临的信号。这些现象不仅在地面可观测到，在卫星上也能观测到，它们在确定异常现象发生的地点和时间等方面具有各自的优势，把两种观测相结合，对异常信号进行立体监测，是实现地震预测预报的发展趋势。至今，对地面观测已经进行了大量的研究和实践，利用卫星进行观测近几年刚刚开始，但已经显示了它监测地震异常的独特能力，引起各方面的高度重视，成为各国建立立体预警系统不可缺少的组成部分。本文重点介绍利用卫星等观测空间电磁异常现象的研究结果。

关键词：地震；电磁；卫星；监测；预警

来源出版物：地球物理学进展, 2007, 22(3): 667-673

被引频次：53

地球磁场对太阳风的加卸载响应与地震

曾小苹，续春荣，赵明，等

摘要：将磁暴过程作为地球磁场对太阳风的加卸载响应，计算分析了北京地磁台1965.1—1979.12和1989.1 —1991.12共18年垂直分量Z的暴时场Dst（Z）加卸载响应比值fD（Z）的变化｡□□fD（Z）异常与其周围550 km范围内发生的Ms≥5.5的中强主震有较好的对应，对结果和在地震预报中应用的可能性进行了初步讨论。

关键词：地震预报；太阳风；暴时场；加卸载响应比

来源出版物：地震地磁观测与研究, 1996, 17(1): 49-53

被引频次：48

寂静的前震与地震预测

赵根模，杨港生，陈化然

摘要：寂静地震的研究具有重要的理论与实际意义。统计表明，与寂静地震有关的前驱波比前震有更大的普遍性，出现的时间主要集中在主震前7 d之内，从一部分中等地震到特大地震之前都已观测到前驱波，这就说明大的脆性破裂之前，较小的缓慢破裂在自然界是存在的而且比较普遍。根据1969年渤海地震和1976年唐山地震的前驱波记录进行定位试验，证明在观测条件具备时，概略的定位是可行的，而且表明前驱波是来源于大震震源附近的信号。可为地震预报特别是短临预报提供重要的前兆信息。

关键词：前驱波；寂静的前震；缓慢破裂；普遍性；地震预测

来源出版物：地震, 2001, 21(1): 69-77

被引频次：646

Earthquake prediction: Physical basis

Scholz, CH; Sykes, LR; aggarwal, yp

Abstract: Earthquake prediction, an old and astrologers alike, appears to be on the verge of practical reality as a result of recent advances in the earth and materials sciences.

来源出版物：Science, 1973, 181(4102): 803-810

被引频次：276

Complex critical exponents from renormalization group theory of earthquakes: Implications for earthquake predictions

Sornette, D; Sammis, CG

Abstract: Several authors have proposed discrete renormalization group models of earthquakes, viewing them as a kind of dynamical critical phenomena. Here, we propose that the assumed discrete scale invariance stems from the irreversible and intermittent nature of rupture which ensures a breakdown of translational invariance. As a consequence, we show that the renormalization group entails complex critical exponents, describing log-periodic corrections to the leading scaling behavior. We use the mathematical form of this solution to fit the time to failure dependence of the Benioff strain on the approach of large earthquakes. This might provide a new technique for earthquake prediction for which we present preliminary tests on the 1989 Loma Prieta earthquake in northern California and on a recent build-up of seismic activity on a segment of the Aleutian-Island seismic zone. The earthquake phenomenology of precursory phenomena such as the causal sequence of quiescence and foreshocks is captured by the general structure of the mathematical solution of the renormalization group.

来源出版物：Journal de Physique I, 1995, 5(5): 607-619

被引频次：274

Aftershock zones of great earthquakes, seismicity gaps, and earthquake prediction for Alaska and the Aleutians

Sykes, LR

Abstract: Aftershocks of shallow earthquakes larger than magnitude 7 in the Aleutians, southern Alaska, southeast Alaska, and offshore British Columbia from 1920 to 1970 were relocated by computer in an attempt to delineate the rupture zones of large earthquakes. Plate tectonic theory indicates that gaps in activity for large earthquakes for the past 10’s to 100’s of years are likely sites of future large earthquakes. Three prominent gaps of this type are delineated: one in southeast Alaska; another in southern Alaska near the epicenters of the great earthquakes of 1899 and 1900; and one in the far western Aleutians. These gaps deserve high priority for study and instrumentation. Large earthquakes appear to be much more regular than smaller shocks in their distributions with respect to space, time, and size. Aftershock zones of events since 1930 that are larger than magnitude 7.8 are longer than 250 km and those less than 7.5 are shorter than 125 km. The rupture zones of events that occurred before 1930 could not be delineated from aftershock locations. Aftershock zones of large earthquakes tend to abut without significant overlap even for rupture zones as long as 1200 km. Nearly the entireAlaska-Aleutian zone from 145°W to 171°E has broken since 1938 in a series of large earthquakes. The rupture zones of five large events appear to form a space-time sequence that progressed from 155°W in 1938 to 171°E in 1965. This sequence is much like the well-known westward progression of activity since 1939 along the North Anatolian fault. Shocks with long rupture zones tend to occur along those parts of the Alaska-Aleutian zone that are relatively simple tectonically. The ends of many aftershock zones of large earthquakes are located at the intersection of major transverse features with the Aleutian arc. Large earthquakes rarely, if ever, reoccur along the same part of a fault zone in less than several tens of years, i.e. within a time less than that for substantial strain accumulation. Events of comparable magnitude that occur soon after some great earthquakes usually involve rupture in a region adjacent to but different from that of the main shock. The March 30, 1965, earthquake of magnitude 7.5, which involved normal faulting in the Aleutian trench, appears to have been triggered by thrust faulting along the adjacent inner margin of the trench in the magnitude 7.9 earthquake of February 4, 1965. Large events of the thrust type are commonly followed within ten years by events involving normal faulting in the adjacent part of the trench. Estimates of average displacements and of the repeat times of great earthquakes from measurements of 20-sec surface waves are systematically too small and do not agree with the meager historic record of great shocks. Other estimates of repeat times vary from 30 to 850 years, but neither of these extremes appears to be typical. The aftershock zone of the April 1, 1946, Aleutian earthquake, which generated one of the largest and most widespread seismic sea waves in the Pacific during this century, was very small. A large displacement of the ocean floor may be responsible for the generation of the large sea wave. An average displacement of 2.4 to 4.1 meters was calculated from amplitudes of 100-sec waves.

来源出版物：Journal of Geophysical Research, 1971, 76(32): 8021-8041

被引频次：235

Rupture zones of large South American earthquakes and some predictions

Kelleher, JA

Abstract: This study attempts to forecast likely locations for large shallow South American earthquakes in the near future by examining the past space-time pattern of occurrence of large(M≥7.7)earthquakes, the lateral extent of their rupture zones, and, where possible, the direction of rupture propagation. Rupture zones of large shallow earthquakes generally abut and do not overlap. Patterns of rupture propagation appear to follow certain trends. These facts, plus the nonrandom behavior of the space-time history of seismic activity, present consistencies that may permit prediction, in a gross sense, of future events. By mapping the rupture zones of large earthquakes(in contrast with plotting only epicenters), it is possible to identify segments of the shallow seismic zone that have not ruptured in many decades. Limited experience elsewhere indicates that these gaps between rupture zones tend to be filled by large-magnitude earthquakes. In certain places it is possible to make approximate estimates of the time of occurrence of the next large earthquake. For at least 300 or 400 years, the entire fault segment near the Central Valley province of central and southern Chile(about 32°–46°S)has fractured about once each century from a generally N-S progression of several large(M≥8)earthquakes. Large earthquakes in this region have almost always occurred to the south of a previous large earthquake. In addition, it is possible to infer a direction of rupturing for two large earthquakes in this century(1928 and 1960). Both these earthquakes fractured southward away from the rupture zone of an earlier earthquake. It would be consistent with these observations if a new series started about the end of this century near Valparaiso(33°S)and progressed southward. In other sections of South America there are several extensive segments of the active seismic belt that have not ruptured during this century. Northern Chile and southernmost Peru(about 17°–25°S)have been relatively aseismic for about 100 years. South of Lima(about 12.5°–14°S), between the rupture zones of the 1940 and 1942 Peruvian earthquakes, there is another significant gap in recent activity. Both these regions are probably areas of relatively high earthquake risk. The northern Peru and southern Ecuador region(about 9°–1°S)has also been relatively aseismic during this century. However, this region differs from the two previously mentioned gaps in that this coastal zone was a region of moderate seismicity during historic times. Perhaps aseismic creep is an unusually important factor in relieving tectonic strain alongthis particular segment of the shallow seismic zone. Another possibility is that large shallow earthquakes in this region have an extremely long recurrence time. Much of the shallow seismic zone of northern Ecuador and southwestern Colombia has ruptured twice during this century. During large earthquakes in this region, the rupturing tends to be directed toward the north or NE. The data for this region suggest that the area to the NE of the 1958 Colombian earthquake may be a region of relatively high earthquake risk.

来源出版物：Journal of Geophysical Research, 1972, 77(11): 2087-2103

被引频次：233

The Parkfield, California, earthquake prediction experiment

Bakun, WH; Lindh, AG

Abstract: Five moderate(magnitude 6)earthquakes with similar features have occurred on the Parkfield section of the San Andreas fault in central California since 1857. The next moderate Parkfield earthquake is expected to occur before 1993. The Parkfield prediction experiment is designed to monitor the details of the final stages of the earthquake preparation process; observations and reports of seismicity and aseismic slip associated with the last moderate Parkfield earthquake in 1966 constitute much of the basis of the design of the experiment.

来源出版物：Science, 1985, 229(4714): 619-624

被引频次：225

Latest aspects of earthquake prediction in Greece based on seismic electric signals

Varotsos, P; Lazaridou, M

Abstract: Since 1983, continuous monitoring of the electrotelluric field has been carried out using an array of measuring stations located at various sites in Greece. The basic physical properties of the transient changes-seismic electric signals(SES)-in the electrotelluric field that are forerunners of earthquakes were first described six years ago. Since then a large body of data has been collected resulting in new insight into various aspects of the method. The present paper reviews the latest developments in SES-based earthquake prediction and describes the current procedures used to predict the epicenter and magnitude of an impending earthquake. A detailed list of the predictions officially issued in Greece during the past 3 years(January 1, 1987-November 30, 1989)is also given. Public warnings were issued well before the most destructive seismic activity.

来源出版物：Tectonophysics, 1991, 188(3): 321-347

被引频次：222

Latest aspects of earthquake prediction in Greece based on seismic electric signals, II

Varotsos, P; Alexopoulos, K ; Lazaridou, M

Abstract: The latest aspects of the physical properties of seismic electric signals(SES)that are used for earthquake prediction in Greece are described. The procedure currently followed for the selection of a site appropriate for the collection of SES and for the electric dipole configuration of such a station are reviewed. The process of constructing a map indicating the seismic areas for which such a station is sensitive is also described. A review of other electrical precursors observed in Greece is given, along with a summary of the newest theoretical models of the generation of SES. The compatibility of these models with the existing data on SES is discussed. Furthermore, a number of unsolved problems are listed, along with suggestions for future experimentation.

来源出版物：Tectonophysics, 1993, 224(1): 1-37

被引频次：211

Earthquake prediction: Variation of seismic velocities before San-Francisco earthquake

Whitcomb, JH; Garmany, JD; Anderson, DL

Abstract: A large precursory change in seismic body-wave velocities occurred before the earthquake in San Fernando, California. The discovery that this change is mainly in the P-wave velocity clearly relates the effect to the phenomenon of dilatancy in fluid-filled rocks. This interpretation is supported by the time-volume relation obtained by combining the present data with the data from previous studies. The duration of the precursor period is proportional to the square of an effective fault dimension, which indicates that a diffusive or fluid-flow phenomenon controls the time interval between the initiation of dilatancy and the return to a fully saturated condition which is required for rupture.

来源出版物：Science, 1973, 180(4086): 632-635

被引频次：153

Premonitory changes in seismic velocities and prediction of earthquakes

Aggarwal, YP; Sykes, LR; Armbrust, J; et, al.

Abstract: Large premonitory changes in ratio of seismic P and S velocities correlate with time of occurrence and size of earthquakes at Blue Mountain Lake, New York. These changes may be related to the opening of cracks or to changes in pore pressure.

来源出版物：Nature, 1973, 241: 101-104

被引频次：149

Some characteristics of foreshocks and their possible relationship to earthquake prediction and premonitory slip on faults

Jones, LM; Molnar, P

Abstract: Foreshocks occur before a large fraction of the world’s major(M≥7.0)earthquakes. Teleseismically located events before major earthquakes from 1914 to 1973 were considered together to examine possible average temporal and spatial patterns of foreshock occurrence. Several days before the main shocks and apparently near the epicenters of them（Δ ≲ 30 km)the activity begins to increase, culminating in a final rapid acceleration of activity in the last day. The acceleration continues up to the time of the main shocks, except for a possible temporary decrease about 6 hours before them. The seismicity increases approximately as the inverse of time before main shock. This relationship is essentially unrelated to the magnitude of the main shock. The magnitude of the largest foreshock is also unrelated to the magnitude of the main shock. In addition, pairs of major events are common. Ten percent of the world’s major events are preceded by other major events within 100 km and 3 months. For foreshocks within each of three sequences studied, the ratio of the amplitudes of the P and S waves were approximately the same, suggesting that the faulting mechanisms are the same for events in each sequence. By assuming an inhomogeneous fault plane on which asperities fail by static fatigue, we derived an equation for accelerating premonitory slip as a function of time, which agrees with the observed time dependence of foreshocks.

来源出版物：Journal of Geophysical Research: Solid Earth, 1979, 84(B7): 3596-3608

·推荐论文摘要·

鲁甸、景谷、康定地震预测的原理、方法及其意义

李德威

摘要：继中长期预测了芦山地震之后，笔者中期预测了鲁甸、景谷和康定地震。例如，鲁甸地震的预测震级为7级左右或6.5级以上（实为里氏6.5级），地点为北纬26°～29°、东经101.5°～105°（实际震中北纬27.1°，东经103.3°），发震时间可能是2014年5月至2015年5月（实际为2014年8月3日）。本文总结了鲁甸、景谷、康定地震和陆内地震热流体物理综合预测的原理、方法和步骤;阐明了在开放复杂地球系统多级物质循环热构造背景下，大陆地壳非均匀流动（“热河”）过程中热能的源、汇、释过程与热灾害链及其地震之间的关系；提出了根据热灾害链时空结构和活动“热河”地震空区相结合进行长期和中期地震预测，与根据热流体直接和间接前兆异常开展立体监测和短临地震预测有机结合的新思路。根据当前热灾害链的演变规律和异动“热河”地震空区分布，进一步分析了西南和华北地震的发展形势，强调华北（特别是东北）的震情极为严峻，短临地震监测和预测已刻不容缓。

关键词：下地壳非均匀流动；热流体撞击；地震预测；鲁甸、景谷、康定地震

来源出版物：大地构造与成矿学, 2015, 39(1): 1-14

联系邮箱：李德威，dewei89@sina.com

由震源到动力学系统──地震模式百年演化

周硕愚，吴云

摘要：地震模式受制于学术思想和科技水平，并随时间演化。研究了从“弹性回跳说”以来的演化轨玶：仅考虑震源区（断层段）→震源区和周围构造相互作用→地壳上部脆性层和深部韧性层与上地幔间的耦合→将地震视为现今大陆变形动力系统自组织演化过程的一种行为。百年演化总趋势是从震源到“动力系统”。特征为：模式时空域扩大、维数上升；从还原论、机械论到整体论、演化论；从经典力学到复杂系统理论及非线性动力学；新模式包容并超越已有模式，向大自然的本性逐渐接近。受益于地震大地测量学观测成就，关联性强的跨学科新地震模式正在探索与成长中。

关键词：地震模式；震源；动力系统；现今大陆变形动力系统模式；地震大地测量学；地震预测

来源出版物：大地测量与地球动力学, 2015, 35(6): 911-918

联系邮箱：周硕愚，cwhezsyw2007@163.com

芦山Ms7.0地震前水平位移和同震水平位移研究

顾国华，王武星，占伟，等

摘要：本文利用“中国地壳运动观测网络（二期）”多个GPS连续观测站观测数据处理结果，将2013年4 月20日四川芦山Ms7.0地震区域参考框架同震水平位移与全球参考框架同震水平位移进行比较，结果表明两组框架解一致，说明两种参考框架均可当作位错参考框架，也即全球参考框架同震水平位移也可视为区域参考框架同震水平位移。区域参考框架下GPS连续观测站地震前的水平位移和同震水平位移结果表明，震前数年，SCTQ站西侧的GPS站构造运动十分显著，而该站水平位移却很小，即出现反常的闭锁。但该站的同震水平位移使其弹性回跳至正常构造水平位移水平，因此SCTQ站震前的位移闭锁是水平位移空间分布中的异常，是芦山Ms7.0地震的前兆。水平位移时空变化表明，该站震前和震时位移完全符合里德的弹性回跳理论。区域参考框架中位移时间系列和同震水平位移的综合研究有助于对芦山地震地壳运动前兆的认识和解释。尽管本文未能直接采用其它GPS连续观测站的资料，但结合本文和其它研究结果可以证实，震中附近其它站地震前后的变化与SCTQ站类似。基于芦山地震前水平位移和同震水平位移及其与前兆关系的研究，本文进一步讨论了GPS监测网的布设、数据处理和分析等问题。

关键词：GPS；GNSS；同震水平位移；弹性位错模型；芦山地震；地震前兆；地震预测

来源出版物：地震学报, 2015, 37(1): 53-64

联系邮箱：顾国华，ggh@seis.ac.cn

地震预测方法Ⅲ：固体潮与地震

赵永红，张丹，王航

摘要：本文从分析研究固体潮与地震发生的相关性的物理机制出发，探讨了用观测固体潮的方法来预测地震发生的可能性.首先，根据地震发生的物理力学过程，给出了地震产生所需要的应力应变条件。实验室岩石力学实验结果和地震应力的观测以及计算结果给出，完整地壳岩石的强度界于10～103MPa之间，断层介质的强度界于1～102MPa之间。其次，将固体潮所引起的全球地壳变形在某个活动断裂带附近产生的形变值与地震触发的应变值进行对比，讨论了固体潮触发地震的可能性。计算分析表明：固体潮在地表对应的应变在10-8的量级，对岩石介质所引起的应力最高可以达到1.6×10-2MPa量级，这样的载荷水平不足以引发地震。最后，本文结合固体潮观测资料表明，地震活动性与固体潮的周期性变化无明显相关性，固体潮对于临震前的含有活动断裂的岩体在某个方向发生地震可能有触发作用，但是由于这种触发作用与地震发震的时间、地点和震级关系的不确定性，从而无法用来作为地震预测的方法。

关键词：地震；固体潮；相关性；物理机制

来源出版物：地球物理学进展, 2015(1): 77-84

联系邮箱：赵永红，zhaoyh@ pku.edu.cn

地震预测方法Ⅰ：综述

赵永红，杨家英，惠红军

摘要：为了解各种地震预测方法的预测情况，本文参考大量国内外研究成果，以地应力法、钻孔应变观测法、GPS法、地下水异常观测法、氢同位素法、电磁异常观测法、电离层异常观测法、潮汐力谐振共振波法等多种方法为载体，对这些地震预测方法进行总结和讨论。探究这些预测方法的原理、异常信息提取、与地震三要素的关系、预测实例及相关的进展。

关键词：地震预测方法；综述；总结

来源出版物：地球物理学进展, 2014, 29(1): 129-140

联系邮箱：赵永红，zhaoyh@pku.edu.cn

地震预测方法Ⅱ：评述

赵永红，王航，杨家英，等

摘要：本文分析了各种地震预测方法对地震的预测情况。首先通过岩石力学实验和有限元模拟计算结果分析了含裂纹岩石材料的变形破坏过程；其次，将实验室岩石破坏过程的观测结果和有限元模拟计算结果外推到地壳岩石的变形破坏过程，从地震发生的物理力学机理出发，对多种地震预测方法进行了评述，研究这些地震预测方法的有效性和尚待解决的问题，同时对利用物理化学方法预测地震的前景进行了讨论。

关键词：地震预测方法；评述；岩石力学实验；有限元模拟

来源出版物：地球物理学进展, 2014(2): 630-637

联系邮箱：赵永红，zhaoyh@pku.edu.cn

云南地区水温异常与地震关系

张彬，方震，刘耀炜，等

摘要：云南地区地震频发，水温观测点分布也最集中，为研究水温异常与地震的关系提供了便利条件。通过收集已经公开发表关于该区域“十五”之前的水温异常与地震的对应关系的文章，较系统地研究了水温异常持续时间、异常空间位置与地震之间的关系。结果表明：水温异常主要是地震短临异常，强震前也存在水温中期趋势异常；一般情况下，地震震级越大，异常范围越广，发震地点通常出现在水温异常集中的区域。

关键词：水温异常；地震预测；震级；正常动态特征；天然地震

来源出版物：地球科学—中国地质大学学报, 2014, 39(12): 1880-1886

联系邮箱：方震，zhen215@126.com

地震预测：从芦山地震到大陆地震

李德威，陈桂凡，陈继乐

摘要：自从1990年以来，通过对青藏高原的调查和研究，认识到下地壳流动同步形成盆地和造山带，并受控于相关洋盆地幔软流圈向大陆的顺层流动和底辟作用。下地壳不均匀流动通过韧脆性中地壳热能—应变能转换孕育地震，部分发震能量通过上地壳脆性断层释放。在地震孕育过程中通常会伴生跨年度干旱和异常降雨，构成热灾害链。近5年内青藏高原东部连续发生汶川、玉树、芦山大地震，形成于从亚东流经羊八井、安多、玉树并分支流向汶川和芦山—康定的下地壳“热河”的仰冲式和侧冲式撞击作用。从2008年9月以来连续发表5篇论文，根据地壳热构造和热灾害链的时空结构对芦山地震的三要素进行了长期和中期预测。2008年9月预测从2013年开始可能发生大地震，2012年9月将鲜水河—安宁河—小江异常热流构造带5年内将发生多个7级地震的首个大震锁定在芦山或西昌。芦山地震只释放了亚东—羊八井—安多—玉树—鲜水河—安宁河—小江“热河”剩余热能中的一小部分，在西昌—会理—昭通地区、道孚—康定地区、通海—石屏地区近5年内很可能发生4个7级左右的地震。此外，华北典型的热灾害链结构表明震情严峻，环渤海地区近3年内很可能发生大地震。从地震热流体撞击机理与地震异常之间的关联性出发，提出了动态立体监测及短临预测地震的思路和方法。

关键词：地震机理；地震预测；下地壳流动；芦山地震；青藏高原

来源出版物：地学前缘, 2013, 20(3): 1-10

联系邮箱：李德威，dewei89@sina.com

加卸载响应比——地震预测与力学的交叉

尹祥础，刘月

摘要：地震预测是世界性科学难题。地震现象虽然复杂其物理实质倒是明确：地震就是地壳块体的快速剪切脆断。相应地，地震的孕育过程就是震源区介质的损伤、演化，并最终导致破坏的过程，这一过程主要是力学过程。抓住这一点就抓住了问题的物理本质。但是，在研究地震预测时遇到的力学问题和通常的（工程）力学问题有所不同。根据地震问题的特点，紧扣地震孕育过程的物理本质，提出了加卸载响应比这一地震预测新思路。文中介绍了加卸载响应比理论的基本科学问题，包括如何对地壳加载/卸载，如何选择适当的地球物理参数作为响应量，以及怎样定义加卸载响应比。用实验研究、数值模拟和理论分析3种基础研究手段，揭示了地震孕育过程中加卸载响应比共同的演化规律:孕震初期加卸载响应比在1附近涨落，之后上升至峰值点，地震不在峰值点发生，而是在下降过程中发生。从峰值点到地震发生这段滞后时间称为T2，T2和震级有关。为了预测地震，必须在全国范围内作加卸载响应比的时空扫描。在时空扫描基础上，结合量纲分析，更多地考虑当地的地球物理条件，拟定了全面预测未来地震的时、空、强的思路。回顾了用这种思路进行地震预测实践的历程。

关键词：加卸载响应比；地震预测；峰值点；量纲分析

来源出版物：力学进展, 2013, 43(6): 555-580

联系邮箱：刘月，liuyue@lnm.imech.ac.cn

Cannot earthquakes be predicted

Wyss, M

Abstract: Robert J. Geller et al. are on shaky ground when they state, in the title of their Perspective, that earthquakes cannot be predicted. In spite of their advice, we should not stop studying the physics of preparation for catastrophic rupture, in the field, the laboratory, and theoretically; neither should we stop measuring crustal parameters that might furnish constraints for physical models; and we should continue researching statistical methods to evaluate prediction claims and to test hypotheses quantitatively.

来源出版物：Science, 1997, 278(5337): 487-490

Earthquakes cannot be predicted

Geller, RJ; Jackson, DD; Kagan, YY; et, al.

Abstract: Research suggests that earthquakes cannot be predicted. Earthquakes result from sudden slips along geological faults, and they are notoriously intractable.

来源出版物：Science, 1997, 275(5306): 1616

Statistical analysis of an ionospheric parameter as a base for earthquake prediction

Li, M; Parrot, M

Abstract: This paper is related to the use of ionospheric density variations to tentatively predict earthquakes. The results of this statistical analysis are presented as a function of various parameters. The ion density was recorded by the low-altitude satellite DEMETER during more than 6 years, and a search for anomalies was automatically conducted with the complete data set. In a second time, a software checked if each anomaly could correspond to an earthquake. The search was conducted at less than 1500 km from the anomaly positions, and until 15 days after the anomaly time. The earthquakes have been classified depending on their magnitude, depth, and position(below the sea or inland). This attempt to predict earthquakes of course generates a lot of false alarms and wrong detections. Nevertheless, it is shown that the number of good detections increases with the magnitude of the earthquakes. In average the number of perturbations is higher the day of the earthquake, and then smoothly decreases for the days before. Earthquakes below the sea are better detected. There are seismic areas close to the South Atlantic Magnetic Anomaly and at high latitudes where the number of natural perturbations is too important to expect a high number of good detections. Finally, when there are several perturbations corresponding to a single earthquake, it ispossible to combine their positions to have a better estimation of the location of the future epicenter. However, uncertainties about the time and the magnitude are large.

来源出版物：Journal of Geophysical Research: Space Physics, 2013, 118(6): 3731-3739

联系邮箱：Parrot, M; mparrot@cnrs-orleans.fr

Earthquake prediction: 20 years of global experiment

Kossobokov, VG

Abstract: Earthquake professionals have for many decades recognized the benefits to society from reliable earthquake predictions, but uncertainties regarding source initiation, rupture phenomena, and accuracy of both the timing and magnitude of the earthquake occurrence have oftentimes seemed either very difficult or impossible to overcome. The problem is that most of these methods cannot be adequately tested and evaluated either because of(a)lack of a precise definition of “prediction” and/or(b)shortage of data for meaningful statistical verification. This is not the case for the pattern recognition algorithm M8 designed in 1984 for prediction of great, Magnitude 8, earthquakes, hence its name. By 1986, the algorithm was rescaled for applications aimed at smaller magnitude earthquakes, down to M5+ range, and since then it has become a useful tool for systematic monitoring of seismic activity in a number of test seismic regions worldwide. After confirmed predictions of both the 1988 Spitak(Armenia)and the 1989 Loma Prieta(California)earthquakes, a “rigid test” to evaluate the efficiency of the intermediate-term middle-range earthquake prediction technique has been designed. Since 1991, each half-year, the algorithm M8 alone and in combination with its refinement MSc has been applied in a real-time prediction mode to seismicity of the entire Earth, and this test outlines, where possible, the areas in the two approximations where magnitude 8.0+ and 7.5+ earthquakes are most likely to occur before the next update. The results of this truly global 20-year-old experiment are indirect confirmations of the existing common features of both the predictability and the diverse behavior of the Earth’s naturally fractal lithosphere. The statistics achieved to date prove(with confidence above 99%)rather high efficiency of the M8 and M8-MSc predictions limited to intermediate-term middle- and narrow-range accuracy. These statistics support the following general conclusions:(1)precursory seismic patterns do exist;(2)the size of an area where precursory seismic patterns show up is much larger than that of the source zone of the incipient target earthquake;(3)many precursory seismic patterns appear to be similar, even in regions of fundamentally different tectonic environments; and(4)some precursory seismic patterns are analogous to those in advance of extreme catastrophic events in other complex nonlinear systems(e.g., magnetic storms, solar flares, “starquakes”, etc.)—that are of high importance for further searches of the improved earthquake forecast/prediction algorithms and methods.

Keywords:extreme events; statistics; forecast; prediction; earthquakes; precursory patterns

来源出版物：Natural Hazards, 2013, 69(2): 1155-1177

联系邮箱：Kossobokov, VG; volodya@ipgp.fr

编辑：王微

Earthquake prediction

Rikitake, T

Earthquake prediction research programmes in a number of countries are reviewed together with achievements in various disciplines involved in earthquake prediction research, i.e., geodetic work, tide gauge observation, continuous observation of crustal movement, seismic activity and seismological method, seismic wave velocity, geotectonic work, geomagnetic and geoelectric work and laboratory work and its application in the field. Present-day development of earthquake prediction research suggests that actual prediction of some class of earthquakes, if not all, may be possible within a period of a few tens of years provided that basic data could steadily be accumulated. Mention is also made about the difficulties in issuing an actual forewarning of earthquake on the basis of the experience of the 1965–1966 Matsushiro earthquakes during which a few long-term predictions of moderately destructive earthquakes were officially issued to the public for the first time in history.

ionosphere; earthquake