地震活动性的时间序列与

   陕西省地震的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 ┃ 合 

──┰────╂────╂────╂────╂────╂────╂───╂────

            14                                21

    ┠────╂────╂────╂────╂────╂────╂───╂────

          43                                52

    ┠────╂────╂────╂────╂────╂────╂───╂────

            28                                36

    ┠────╂────╂────╂────╂────╂────╂───╂────

          ┃ 117 ┃                            ┃ 132

    ┠────╂────╂────╂────╂────╂────╂───╂────

         57                                70

    ┠────╂────╂────╂────╂────╂────╂───╂────

            39                                46

──┸────╂────╂────╂────╂────╂────╂───╂────

合计:      ┃ 298 ┃  21    31                  ┃ 357

───────╂────╂────╂────╂────╂────╂───╂────

        2.474 ┃  1.322 ┃  1.491 ┃  0.699 ┃   0.0    0.0 ┃

 

 

                ㈡、表征地震活动性的度量

 

    我开始时完全从预报的角度出发,选用单位时间T(旬、月、年或十年)内的震数和

最大震级,构成时间序列,再作周期分析和预报。用这两个数来表征某个地区在某个时段

内的相对地震活动性,似乎都不如用单位时间内所有地震释放的总能量ΣE来得合理,

而每个地震的能量E则可由其震级M算得(古登堡──李希特公式):

                    11.8    1.5      ┅┅┅┅┅┅ 

如果定义地震的功率密度函数为:

                             ΣE

                      ──────      ┅┅┅┅┅┅┅┅┅┅┅┅┅ 

                             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        b·M    ┅┅┅┅┅┅┅   

其中N为震级属于M级左右某个区间(如M±0.5)内的地震数目。或者,按公式⑴

用K=    代替震级M,则有:

                  N    α    r·K

按文献所说,这里r=0.43(相当于b=1.5×0.43=0.65)是相

当稳定的。将表一中近1400年的全省中──大震资料,点绘成图2甲。表二是陕西省

`&100`近11年的中、小震资料:

 

表二.    67──77年陕西省地震的震级分布

 

    震级    ┃ 0.0-0.9┃ 1.0-1.9┃ 2.0-2.9┃ 3.0-3.9┃ 4.0-4.9┃ 5.0-5.9

  ──────╂────╂────╂────╂────╂────╂────

    震数      29  ┃ 149 ┃  43                

  ──────╂────╂────╂────╂────╂────╂────

     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级地震与所对应的高峰期

 

─────┰──────┰──────┰───────┰───────────

高峰期序号┃                                      

─────╂──────╂──────╂───────╂───────────

 高峰年份 ┃  前 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年