地震活动性的时间序列与
陕西省地震的790年周期
陕西省气象局 张时钊
试用自己的质周期分析法,我曾对浙江省地震局逐月发的三年(1976─78年)
地震目录作过第一次周期分析和预报,看来还可能有些用处,自然想用更长期、更
详细的历史资料再作一次试验研究。陕西省地震队分析组于1973年已汇编出版《陕西
省地震目录(公元1177─1972年)》,共编录4级以上地震68个,4级
以下地震416个。78年,西安地震办又续编了1973─78年部分,记录地
震110个(全部为4级以下)。本文对这两份历史地震资料进行了初步分析,得到一些
结果。首先在理论和实践上建立了地震周期性分析基础,即论证了在历史资料不完整的条
件下(长期的地震历史资料总是非等精度的记录,总是不完整的),如何构成可供分析的
地震活动性时间序列;其次,发现了陕西省地震活动性的一个确凿无疑的790年长周期
。不过由于缺乏近期的地震资料,没有完成中、短周期的分析,因而也未作出长期预报。
㈠、陕西省地震区的划分
把陕西省的各个历史地震点绘到地图上,立即可以发现震中分布的点聚性质。陕西省
地震队76年已经出版过包括所有历史地震的震中分布图,后面附的图1则仅点有
1957─77年全部由仪器测定的中小震,它们表现出来的分布性质几乎没有两样。从
图1很容易划分出五个区:1、2、3、4、5;每个区内的地震分别团聚在一起,因而
可以彼此分开。只有区2大而杂,但不易再细分了。区6很少有地震,历史上记录也不多
,姑且另作一区。我们也可以看出,所有地震好象总是分布在河流两侧。实际上,地震的
发生是一种地质运动,所以它的分布必然在一定程度上反映了地质构造。正由于地质构造
、首先是断层或皱摺的走向,既左右河流的发育,也规定着地震的分布,三者是非常一致
的。
我们这样划分地震区,是为了把地震序列分解为属于各别地区的子序列,对这些子序
列分别分析得出预报,使预报落实到较小的区。为了使序列预报的落区更正确些,当然把
这些区分得愈细愈好。但是,当地震区划得太小时,每区的地震资料也就少了,将使统计
分析变得不可靠或不可能。另外,如图1所示,我们是简单地用经、纬线作为区界,这完
全是为了各别地震的归属更容易些。由于按地震点聚性质划的各个区,其边缘地带地震少
而小,使这样不自然的人为做法,对分析的正确性不会有很大的影响。
现把这5个区简述如下:
①.汉水上游区:北纬34°以南,东经107°41以西。断层走向为东东北─
─西西南。
②.汉水下游区:北纬34°以南,东经107°41以东。地质构造复杂,在许
多西北──东南向的断层间交错着一些东西向的断层。这应该是震中分布比较散的原因。
③.渭河上游区:北纬34°以北,东经108°以西。断层走向为西北──西南。
④.泾河流域区:北纬34──35.5°,东经108──109.5°之间。地
质构造的特点是,在几条偏东西向的大断层间,还贯穿着许多西北──东南向的短断层。
⑤.洛河下游区:北纬34──36°之间,东经109.5°以东。断层走向,北
部都为东北──西南向,向南则逐渐转为东西向。
②.陕北丘陵区:东经109.5°以西为北纬35.5°以北,东经109.5°
以东则为北纬36°以北。本区属于伊陕台地,仅北部有少数东北──西南向的短断层,
近期几乎没有地震。
这6个区域中,无论从发震频数或最大震级来看,都以第④、⑤为最重要,见表一。
该表只统计了公元600──1977年为止的大于3级的地震。(对于解放前的4级以
下的中、小地震,既无仪器测量记录,又没有详细研究,资料中也没有列出具体的
经、纬度和震级;但既见之于前人记录,必为有感地震,姑且认为大于3级。另外,原资
`&100`料中一个编号,常包含数个地震,凡容易分辨者,则已区分为各别地震)。
震 级 ┃ 3 -3.9 ┃ 4 -4.9 ┃ 5 -5.9 ┃ 6 -6.9 ┃ 7 -7.9 ┃ ≥ 8 ┃ 合 计
──┰────╂────╂────╂────╂────╂────╂───╂────
┃ 1 ┃ 14 ┃ 1 ┃ 6 ┃ 0 ┃ 0 ┃ 0 ┃ 21
┠────╂────╂────╂────╂────╂────╂───╂────
区 ┃ 2 ┃ 43 ┃ 4 ┃ 4 ┃ 1 ┃ 0 ┃ 0 ┃ 52
┠────╂────╂────╂────╂────╂────╂───╂────
┃ 3 ┃ 28 ┃ 5 ┃ 2 ┃ 1 ┃ 0 ┃ 0 ┃ 36
┠────╂────╂────╂────╂────╂────╂───╂────
┃ 4 ┃ 117 ┃ 8 ┃ 5 ┃ 2 ┃ 0 ┃ 0 ┃ 132
┠────╂────╂────╂────╂────╂────╂───╂────
号 ┃ 5 ┃ 57 ┃ 2 ┃ 8 ┃ 1 ┃ 1 ┃ 1 ┃ 70
┠────╂────╂────╂────╂────╂────╂───╂────
┃ 6 ┃ 39 ┃ 1 ┃ 6 ┃ 0 ┃ 0 ┃ 0 ┃ 46
──┸────╂────╂────╂────╂────╂────╂───╂────
合计: N ┃ 298 ┃ 21 ┃ 31 ┃ 5 ┃ 1 ┃ 1 ┃ 357
───────╂────╂────╂────╂────╂────╂───╂────
N ┃ 2.474 ┃ 1.322 ┃ 1.491 ┃ 0.699 ┃ 0.0 ┃ 0.0 ┃
㈡、表征地震活动性的度量
我开始时完全从预报的角度出发,选用单位时间T(旬、月、年或十年)内的震数和
最大震级,构成时间序列,再作周期分析和预报。用这两个数来表征某个地区在某个时段
内的相对地震活动性,似乎都不如用单位时间内所有地震释放的总能量ΣE来得合理,
而每个地震的能量E则可由其震级M算得(古登堡──李希特公式):
E = 11.8 + 1.5 M ┅┅┅┅┅┅ ⑴
如果定义地震的功率密度函数为:
ΣE
P = ────── ┅┅┅┅┅┅┅┅┅┅┅┅┅ ⑵
S·T
娛街校女槲玫厍婊雍偷ノ皇奔洌阅诘母鞲龅卣鹉芰俊5ゾ凸娇矗耆梢园眩印
T都取得很小,从而把P值作为表征时、空上某个点的绝对活动性。不过,地震在时空上
的分布是不连续的,一般又很稀少,为使时间序列有意义,S、T都必须取比较大的数值。
在S、T取定之后,P值的时间序列就等于总能量ΣE的的序列。但这样的序列
,各元素数值相差太悬殊了,难以分析。在作降水量序列的分析预报时,我曾把各元素事
先取对数而获得成功,所以自然也把这些总能量之和ΣE,再按公式⑴取对数:
(ΣE)-11.8
M = ────────────────
1.5
并特别称之为能级M,以区别各别地震的震级M。但这次采用这个新数量后,使用的
结果竟与最大震级序列无大差别。原来,一般说来,单位时间内的各个地震中,总有一个
震级较大,其它小震的能量相对地就显得微不足道,以致在计算能级只保留小数一位时这
些小震几乎不影响计算结果,仅当几个较大的地震震级相差很小(≤0.2)时,能级才
比最大地震的震级大 0.1 - 0.2 左右。 例如一个7级地震的能量为1.995×10
尔格,两个7级地震的总能量为3.990×10尔格,其相应的能级为7.2,只增
加0.2。
当S、T充分大时,在确定的S、T内的地震在震级上的分布服从(古登堡):
N = a - b·M ┅┅┅┅┅┅┅ ⑶
其中N为震级属于M级左右某个区间(如M±0.5)内的地震数目。或者,按公式⑴
用K= E 代替震级M,则有:
N = α - r·K
按文献所说,这里r=0.43(相当于b=1.5×0.43=0.65)是相
当稳定的。将表一中近1400年的全省中──大震资料,点绘成图2甲。表二是陕西省
`&100`近11年的中、小震资料:
表二. 67──77年陕西省地震的震级分布
震级 M ┃ 0.0-0.9┃ 1.0-1.9┃ 2.0-2.9┃ 3.0-3.9┃ 4.0-4.9┃ 5.0-5.9
──────╂────╂────╂────╂────╂────╂────
震数 n ┃ 29 ┃ 149 ┃ 43 ┃ 7 ┃ 0 ┃ 1
──────╂────╂────╂────╂────╂────╂────
n ┃1.48┃2.18┃1.64┃0.90┃ / ┃0.00
由表二的数据点绘成图2乙。由甲、乙两图可见,它们的线性关系都比较好,只有一点偏
离直线较远,这是因为表一中4─4.9级和表二中0─0.9级地震数太少了,看来都
是由于资料不够精确造成的, 文献中对解放前的 4─4.9 级地震 只录有一种
4.75,当然不够确切,而近期资料则遗漏了一些<1级的小震,也是显然的。其次,
甲乙两直线的斜率是相等的,由图求得陕西省的b=0.6,对于同一地区来说,b值果
然很稳定。
地震震级分布的这种统计性质,使得任一相当大的地震区面积S和相当长的时间单位
T内,最大地震的震级只有一个(最多两个)。按上所述,这时完全可以用最大震级代替
总能量的能级来表征相对的地震活动性。显然,这最大能级由图2中的直线与横坐标的交
点所决定。又由于这些直线的斜率,至少对同一个地区来说是稳定不变的,所以直线与纵
坐标的交点(即式⑶中的截距a,相当于0级地震的频数之对数),以及横坐标一定时的崐纵坐标值(相当于某个规定的标准震级的频数,黎兹尼琴科正是以它作为地震活动度的定
义的),甚或某级以上的所有地震的总频数(累计地震法),都可以作为相对地震活动性
相互等价(仅差一个因子)的度量(参见《地震学基础》220──224页)。
应该指出,上述文献中所说的地震活动度,是一个地震区在很长的、常常是全部有记
录时期内的统计特性,势必被认为是永远不变的,至少也被认为在今后很长时期内是不变
的。而我们作地震序列分析预报时,需要考察的却是地震活动性在时间上的变化。当序列
的时间单位取得很小,以更高的精度描述这种变化时,公式⑶的统计特性将不复存在,这
时只能用式⑵定义的地震功率密度或总能量的对数(即能级M)作为地震活动性的唯一
度量。
㈢、陕西省地震的790年周期
《陕西省地震目录》提供了有史以来长达3000年的地震资料,这允
许我们去分析、寻找年以上的长周期。但是这两份目录完全不是等精度的观察记录,古代
中小震被遗漏的必然很多。而要作周期分析,记录必须是长期连续且在同一条件下取得的
。好在我们只研究其中的长周期,时间单位可以取得很大,例如十年。我们已看到,十年
时间内的震级分布已较稳定地满足古登堡公式⑶,如上所述,这允许我们用震数或最大震
级作为地震活动性的度量,以构成有意义的时间序列。只是在公元600年前的1700
年中,该资料只录有8个4级以上和23个4级以下的地震,资料太少,作分析时只好割
爱了。
历史文献中被遗漏的必然是比较小的地震,对最大震级这个数值似乎没有什么影响,
本来应当是较为实用而方便的活动性度量,可惜的是该资料中对4级以下的历史中小震,
均未列出具体的震级,所以只能用震数这个数值。若用4级以上的地震累计数,则从公元
600年至今的1400年中又只有60个地震,资料也太少了,最后不得不假定该资料
中所有震级不明的地震(有感地震)都是3级以上的,而近年来的地震也只挑取该级以上
者,这就是说采用3级以上的地震累计数作为相对活动性的度量。3级地震是有感的,在
史料中应有记载,只因年代久远,资料失散,遗漏和错误依然免不了。尤其是因为记录不
详,同时期的地震更常常混而为一,使震数减少,我仅揣度原文含义,能分辨者尽量区别
为各别地震。虽然如此,解放前的地震还是会比实际少了许多的。
这样,从原始资料中按第一节所分的6个区分别摘录出6个十年震数(
每十年内3级以上地震总数)序列,然后连续作3次平滑化时程为10时间单位(100
年)的滑动平均,绘得1──6区的6条曲线如图3所示,其中最下面的第7条是全省震
数的滑动平均曲线。从图⑶可看出:
①.全省的地震活动性在8世纪和16世纪有两个十分显著的高峰期,这个近800
年的长周期在第4条曲线上(泾河流域区)看得更清楚,两个高峰的高度和形状都相近。
这显然是由于:该区是历代政治文化中心──长安(西安)的所在地,历史资料特别丰富
,早期记录很少遗漏之故。曲线4的第一个高峰约在785年,这时有西安788年2月
12日的5.75级地震,第二个高峰在1575年,这时有西安1568年5月15日
6.75级地震。姑且把这两个高峰的时间差1575-785=790年作为精确的周
期长度,那么从785年上溯790年,即在公元前5年,应是另一个高峰。查资料,公
元前7年11月11日确有长安6级地震。这三个地震也是本区(4区)的最大地震,另
一次1487年8月10日临潼6.25级地震当属于第二个高峰,其震中靠近第5区。
实际上,这800年左右的大周期,在资料后面所附的关中地区地震年度分布图中
,也可以明显地看出来。
②.其它5个区,由于相对来说属于边远地区,古代记录势必残缺,所以8世纪的高
峰期均不明显。但是,大震记录还是有的。如第2区的788年3月8日安康6.5级地
震,第5区的753年5月27日华山6级地震。只因4级左右的地震大部被遗漏,高峰
期未能很好地显示出来。不过16世纪的高峰期,对所有6个地区来说,都有,而且很一
致,可以认为这790年的长周期是全省性的,只是位相略有不同。从图可见,波峰到来
的先后次序大概是:5区、2区、4区、1区、3区,最后才是6区。如2区的1501
年1月19日朝邑7级地震,1556年1月23日华山8级大震(1457年临潼地震
虽列在4区,但位置也靠近2区),都比较早。这与文献上说的汾渭地震带内的地震,有
从东北向西南迁移的现象是一致的。
当然也可能有另有稍短的长周期,其周期长度对6个地区来说也可能不完全相同,但
从图3看来,它们至少是不明显而不重要的。所有曲线的后部,都有上抬现象,尤以第2
区曲线为最,这大概是由于近年来用仪器记录到的3级以上的地震较过去无仪器观察时为
多所造成的假象。另外,2区地大而构造复杂,2、3、5区都由于省界的限定,不足以
代表一个完整的地震单元,对曲线形状也可能有所歪曲。
`&120` 最后,把本省全部6个区6级以上的9次地震摘录于下表:
表三 陕西省≥6级地震与所对应的高峰期
─────┰──────┰──────┰───────┰───────────
高峰期序号┃ 1 ┃ 2 ┃ 3 ┃ 4
─────╂──────╂──────╂───────╂───────────
高峰年份 ┃ 前 790 年 ┃ 前 5 年 ┃ 785 年 ┃ 1575 年
─────╂──────╂──────╂───────╂───────────
┃⑴前 780 年 ┃⑴前 7 年 ┃⑴ 788 年 ┃⑴1487年 临潼 6.25级
┃ ┃ ┃ ┃
≥6级的 ┃ ┃ ┃ 安康 6.5 级 ┃⑵1501年 朝邑 7 级
┃ ┃ ┃ ┃
┃ 岐山 6-7级 ┃ 长安 6 级 ┃⑵ 793 年 ┃⑶1556年 华山 8 级
┃ ┃ ┃ ┃
地 震 ┃ ┃ ┃ 华县 6 级 ┃⑷1568年 西安 6.75 级
┃ ┃ ┃ ┃
┃ ┃ ┃ ┃⑸1704年 陇县 6 级
─────┸──────┸──────┸───────┸───────────
由表3可以看到,从全省来说,这790年长周期已有4次高峰记录。第一个高峰是公元前
780年的岐山6-7级地震。这次地震,在周幽王三年,记载称泾、渭、洛三川皆震,
当遍及3、4、5区。再往前在公元前1600年应是高峰期,可惜已无历史记载。最早
的地震记录是公元前1177年的岐山地震,时当周初,而岐山是周的发源地,才被记录
下来,可是它的震级仅达4-5级,所以不是高峰期。在有史以来的3000年中,竟能
完整地找出这4个高峰期,相邻的两高峰期间的间隔都是790年左右,且已知的9次大
震都落在这些高峰期中而无一例外,您能不相信这790年的长周期是确凿无疑的吗?
从第4个高峰期(1575)再外推790年,下一次高峰期将出现在2365年前
后,即要到24世纪中叶,陕西省才有可能发生6级以上的大震,在此之前是不可能发生
象华山大震那样的有严重破坏性的地震的。正因为眼下处在790年大周期的低潮期,所
以象1976年全国地震中心由康藏高原和台湾(这无疑是地震最严重的省区,但每当地
震活动中心北移时,它的地震活动性总是大大下降了)迁向华北时,中国北部许多地方确
也发生了破坏性大震,陕西省的地震活动性也有所提高,但还是没有发生破坏性地震。尤
其是今后数十年内,另一个约40年的周期也将转入低潮期,5-6级的地震更不会发生
,我认为这个看法是有根据而可信的。
1979年7月底写 10月4日改 于彬县气象站
参考文献:
1.质周期分析和地震预报 张时钊 1979年
见《浙江地震》 1979年第一期
2.陕西省地震目录(公元前1177年至公元1972年)
陕西省地震分析组汇编 1973年
3.陕西省地震目录(1973年至1977年)
西安市革命委员会地震办公室 1978年
4.陕西省地震震中分布图
陕西省地震队编 1976年
5.地震学基础 中国科学院地球物理研究所编著
科学出版社 1977年