1. 引言
地震破坏性极强,它的破坏力足以使大片城市夷为平地,给社会和人民带来巨大损失 [1] 。为了减小地震带来的破坏,对地震进行准确预测十分必要。我们既可以从地壳的运动规律,来研究未来时间内的地震发生情况 [2] ;也可以从以往发生过的地震记录中寻求规律,以达到对未来时间地震发生情况的预测。本文利用后者,通过分析往年的全球统计数据,建立地震预测模型,从而对地震的发生进行科学的预测。之后,通过对近四十年来的全球地震数据进行研究,运用神经网络方法预测未来中国的地震发生情况。
全球强震(此处指大于或等于里氏6.0级的地震)带主要集中分布在环太平洋地震带和欧亚地震带,而且身处其中的不同的地区间地震发生情况是互相关联的,我国亦处在其中 [3] 。因此我国的地震活动和全球的地震活动是密切相关的,尤其是板块边界的强震活动更是如此。大地震往往发生在断裂转折处或断裂交叉处,因为这些区域受特殊力学原因的影响,往往不易断裂或容易造成断裂的阻碍,只有当致断能量达到一定程度时,才可能造成这些区域最后剧烈断裂 [4] 。这一点说明,发生一次大的地震需要一定的时间积累,同一地震带内一般不会连续爆发两次或以上的强震 [5] 。考虑到地震发生是受自然因素影响的,是一个能量释放的过程,所以可以假设在同一个区域相对较短的时期内发生一次强震后,不会连续发生第二次强震。而据全球的强震活动和板块边界分布,通常可把全球分为16个强震活动区(图1)。本文分析使用国家地震数据科学共享中心和国家台网大震速报目录数据。数据从1999年1月1日起,实时更新,发震时刻为北京时。通过分析数据,找出16个强震活动区的强震活动次数与其后一年我国强震发生与否的关系,用它预测我国2018年是否会有强震(6.0级以上)发生,并预测震级范围。
2. 地震长期预测神经网络模型
2.1. 模型的假设
在建立预测地震模型的过程中,我们假设:
![](//html.hanspub.org/file/3-2620604x9_hanspub.png)
Figure 1. Global pleistoseismic zone distribution map
图1. 全球强震带分布图
1) 不考虑远离16个强震活动区的区域所发生的强震对其后一年我国强震发生的影响;
2) 在短期时间(15天)内,各地震带不发生两次或以上强震(发生强震需长时间积蓄大量能量);
3) 将两个月定为统计强地震的单位时间;
4) 在考虑强震影响时,忽略6.0级以下地震的影响;
5) 假设一段时期内,弱震和中强震的产生有一定的规律,而强震的产生在短期时间内无规律。
2.2. 符号说明
Ms:表示地震等级;
:表示第i年j区的最大震级强度;
:归一化处理后第i年j区的震级指数。
2.3. 神经网络模型的建立与求解
BP神经网络(Back Propagation Network)是一种迭代梯度算法,用于求解前馈网络的实际输出与期望输出间的最小均方值 [6] 。BP神经网络的分类功能来自神经元节点的非线性,所以,它可以实现由N维空间向M维空间的高度非线性映射,它具有较强的自适应性和自学习功能,而且具有高度的容错性。当个别输入信号值出现畸变或干扰时一般都不会影响识别结果。因此从理论上讲,利用神经网络模型,可以更好地进行地震预测 [7] 。
首先,我们进行数据的归一化处理 [8] 。本文所用地震目录为中国地震局分析预报中心提供的全国地震目录。从地震目录中挑选出Ms ≥ 6的地震。以年为单位进行地震筛选,选取每年的最大震级,构成地震时间序列。将选出地震的震级按照下面公式归一化:
BP神经网络的输入节点,输出节点为下一个年的最大地震震级。根据各年中各强震活动区以及中国的最大震级统计数据,可以通过选取近40年的34个样本,以15年以及其后一年的数据作为训练集,对神经网络模型进行训练,通过选择不同的函数进行拟合,得出神经网络模型。将剩下的样本作为测试集合,通过测试验证,由模型得出的预测值与真实值之间的误差平均为0.5级。由模型得出的预测值与真实值之间的误差统计如表1所示。
2.4. 模型的修正
利用剩余的样本进行预测,发现预测值与真实值之间有一定的误差,但是总体上误差相对较小。误差产生的原因是由于之前筛选数据时,只把6.0级以上的地震数据考虑在内,6.0以下的数据全部被剔除处理造成的 [9] 。在神经网络地震预测模型中加入6.0级以下地震数据,经过多次调整修正,得到的部分结果如表2所示。
以误差阈值为0.5,在阈值内的预测认为是准确预测。
通过计算,数据误差均值在0.5级以内的占总体测试样本的百分之八十,震级误差较小。在真实情况中,这样的误差对于预测和判断地震情况的影响小,拟合效果好。通过对误差值的计算,我们认为神经网络模型所得结果比较令人满意,此神经网络模型基本可用于预测地震。
根据2017年及之前的地震数据,用已经训练好的神经网络模型对中国2018年的强震发生情况进行预测,预测中国2018年发生强震的可能性较大,预测最高震级达到6.53级或以上。
![](Images/Table_Tmp.jpg)
Table 1. Error statistical table of the predicted value of neural network model and the real value
表1. 神经网络模型的预测值与真实值的之间的误差统计表
![](Images/Table_Tmp.jpg)
Table 2. The predicted value of neural network model and the real value
表2. 神经网络模型的预测结果与真实结果
3. 地震短期预测灰度模型
3.1. 短期预测的理论依据
原则上说,地震发生的时间主要由应变能和致应变能增加点的机械承受力这一对矛盾来共同决定。从这个层次上分析,地震的爆发是不定时的,随时都可能发生。但是,地球还受到月球、太阳及其它行星的吸引。这个引力作用于地球,也作用于地壳 [10] 。各天体运转都有相对固定的周期,所以从某种程度上来说,同一地点地震的发生也是有据可依的 [11] 。因此,我们可以根据一个地区前五年的地震发生情况,观察时间与震级之间的规律,然后利用灰度预测方法对数据进行拟合,建立拟合模型。然后利用后五年的数据对模型进行检验。
3.2. 灰度预测模型的建立与求解
通过数据收集,做出中国近年来地震发生时间与震级强度散点图(图2),应用所得地震模型对数据进行分析。通过模型分析,选择用以下函数式拟合地震级数:
其中,y表示发生地震的时间,x表示地震的级数,
分别为拟合模型中的三个参数。
![](//html.hanspub.org/file/3-2620604x16_hanspub.png)
Figure 2. The time and magnitude of earthquake in China from 2006 to 2015
图2. 2006~2015年中国地震发生时间与震级强度散点图
使用灰度预测方法,将
分别赋值以[37,45]、[
]、[5,6]区间中的任意值,通过数次拟合后的对比分析,发现当
时得到的预测结果与真实值最能达到吻合。
得到拟合模型为:
使用地震短期预测灰度模型将近年来地震情况进行拟合,分析拟合结果得到:拟合数据平均误差在0.5级以内的数据占总体数据的近80%,均值误差在0.8级以内的占总体数据的90%,拟合效果好,我们认为建立的地震模型较为理想。
2008年是中国地震多发年,为验证模型效果,以2008至2009上半年为例,模型结果与真实值对比统计如表3所示。
4. 地震预测模型评价与讨论
4.1. 模型的优点
本文在建立模型时,首先通过常规的方法对往年的数据进行统计分析,根据地震的不可预测性以及随机性,分别选取神经网络模型和灰度预测模型进行长期和短期地震预测。通过分析以往的数据即可对未来一段时间内的地震情况进行预测,模型适用性强,预测效果好。另外,本文所建立的神经网络模型和灰度预测模型对于统计无规律事件较为通用,可以扩展应用于预测其他无规律自然灾害的发生,例如火山爆发等。
4.2. 模型的缺点
本文虽然用了大量地震数据进行统计分析,但是在做预处理的过程中,忽略了爆发地震点不在16个强震活动区的地震数据,可能会因此造成一定误差;同时,为了便于建立模型,寻求主要因素,我们忽略了许多其他的相关因素。地震的发生与地质构造以及地壳运动紧密相关,所以在进行预测时,加入对地质等自然因素的考虑,才能更加有效地对地震的发生以及震级的强度进行预测。
![](Images/Table_Tmp.jpg)
Table 3. The real value of earthquake magnitude and fitted value of gray prediction model
表3. 实际震级与函数拟合的理论震级比对表
5. 结语
地震作为最严重的自然灾害之一,一直威胁着人类和平稳定的生活,给人类社会带来重大的灾难,威胁人类的文明发展进程。本文分别使用神经网络模型建立了长期地震预测模型,使用灰度预测方法建立了短期地震预模型,从而对地震进行科学预测,以减少地震带来的损失。
本文首先建立了地震长期预测模型。将每15年划分为一个预测时间单位,并将16个地震活动区分别编号,以此对数据进行初步筛选。考虑到地震发生是一个能量释放的过程,因此假设在同一个区域相对较短的时期内发生了一次强震后,不会连续发生第二次强震。考虑到地震发生的随机性与不可预测性,我们选用BP神经网络模型进行拟合预测 [12] 。通过多次调整参数拟合得到一个结果相对较好的模型,通过模型可知,16个强震活动区的地震活动次数以及震级强度对中国下一年的地震发生情况有较明显影响。之后我们利用该模型对2018年中国地震发生情况进行预测,预测结果为2018年中国极有可能发生强震,最高震级达到6.53级或以上。
之后,本文通过灰度预测方法,建立了短期地震预测模型,利用前5年的地震数据进行预测。首先通过筛选得到10年间发生地震的情况,然后对发生地震的时间以及震级强度两组数据对比分析,发现地震发生的趋势并进行拟合。利用灰度预测方法所得的理论震级与真实值进行比对,对模型进行验证,通过修正参数,得出在误差允许范围内的最佳拟合模型。
本文所建立的模型基本完成了对未来长期及短期内地震的发生时间以及震级强度的预测。而且根据不同的情况分别采取了神经网络模型和灰度预测方法,并对参数进行了修正,使得模型的精度提高。但是,我们的拟合结果仍存在一定误差,这是由于最初对数据的处理方式和简化模型所做假设造成的。在将来的模型改进中,将加入全球范围内所有的地震数据进行统计,并建立修正模型。同时,应更多的将地理因素,如月亮赤纬角极值、厄尔尼诺和拉尼娜对地震的激发作用考虑进模型中 [13] ,进一步从地球构造、地壳运动等理论角度完善预测模型。
基金项目
本文得到国家自然科学基金(51709022),水力学与山区河流开发保护国家重点实验室(四川大学)开放合作基金(SKHL1623),成都信息工程大学科研基金(KYTZ201710)资助项目。