《神经网络在福建及其沿海地区地震预报中的应用》由会员分享，可在线阅读，更多相关《神经网络在福建及其沿海地区地震预报中的应用（5页珍藏版）》请在金锄头文库上搜索。
1、神经网络在福建及其沿海地区地震预报中的应用于波 叶振民地震的发生会给社会带来不同程度的破坏，造成生命、财产的损失。目前地震学只是一门观测科学，我们只能加强观测，通过各种经验方法来进行预测，确定重点监测区。由于地震学异常与地震之间存在较强的不确定性，通常一种地震学指标在一次地震中出现，而在另外的地震中不出现或不明显，因此，单一使用某一个地震学指标的预报方法就表现出它的局限性。而如何解决预测结论的权值以及各指标组合时相互间的非线性联系对最终预测结论产生影响就成为难题，神经网络技术的发展给我们提供了一个解决这个问题较好的方法。它通过分析各个指标，建立其内在联系，并根据这种联系对目标进行预测，能有效地
2、避免了采用单一指标进行预测的片面性，大大提高了预测的精度。本文就 BP 神经网络算法及其在地震预报中的具体运用进行初步的探讨。一、神经网络知识介绍人工神经网络，是一种旨在模仿人脑结构及其功能的信息处理系统。通常一个人工神经网络的功能主要是由两个方面决定的：一是网络的拓扑结构，也就是各人工神经元间相互连接的方式；二是网络的学习和运行规则，即网络中连接权值的调整规则。1、人工神经元的结构：（图 1 略）图 1 中 x1，x2xixn 分别代表来自其它神经元轴突的输入。wj1，wj2wjiwjn 则分别表示神经元 1，2，in 与第 j 个神经元的突触联结强度，即权值。这个权值可正可负，正权值表示兴
3、奋型突触，负权值表示抑制型突触。f()是转移函数,也称激活函数，其作用是模拟生物神经元所具有的非线性转移特性。2、人工神经网络的拓扑结构人的大脑皮层的横断面上有 36 层神经细胞，即具有多层结构。因此，人工神经网络也按层排列。多层网络由单层网络级联而成，即网络中各神经元接受前一层各神经元的输出。界于输入和输出层之间者称为隐层，它们不直接与外部打交道。隐层的加入大大地提高了人工神经网络对信息的处理能力。经过训练的多层网络，具有在任意精度下，实现矩阵 X 转成矩阵 Y 的任意非线性映射的能力，其根本原因就在于引入了隐层以及非线性转移函数。3、人工神经网络的工作过程人工神经网络最使人们感兴趣的特性之
4、一便是它的自学习的能力。如果要使人工神经网络具有某种智能，必须先对它进行训练，对网络来说就是要学会要它做的事情。学习就是相继给网络输入一些样本模式，并按照一定的规则（即学习算法）调整网络各层的权矩阵，待网络的各权值都收敛到一定值时，学习过程便告结束。二、误差反传播(Error Back Propagation)算法，BP 算法在了解了神经网络的基本知识后，我们开始介绍本文用到的 BP 棗误差反传播算法。BP 网络是应用最广泛的网络，现已成功地应用于图像识别、预测预估、语声变换、数据压缩、模式辨识和自动控制等领域。据统计，应用 BP 算法的神经网络高达 80%。因此本文选用它作为进行地震预测的方
5、法。下面介绍其工作原理，并给出其在福建及其沿海这一特定区域地震预报中的应用。BP 算法的基本思想是，学习过程由信号的正向传播与误差的逆向传播两个过程组成。正向传播时，模式作用于输入层，经隐层处理后，传向输出层。若输出层未能得到期望的输出，则转入误差的逆向传播阶段，将输出误差按某种子形式，通过隐层向输入层逐层返回，并“分摊”给各层的所有单元，从而获得各层单元的参考误差或称误差信号，以作为修改各单元权值的依据。这种信号正向传播与误差逆向传播的各层权矩阵的修改过程，是周而复始进行的。权值不断修改的过程，也就是网络的学习（或称训练）过程。此过程一直进行到网络输出的误差逐渐减少到可接受的程度或达到设定的
6、学习次数为止。每组数据有 6 个项目作为输入，目的是预测震级，因此选择输入层为 6 个节点，输出层为1 个节点，隐层节点采用几何金字塔规则确定为 2 个。因为单极性 Sigmoid 函数 g(s)=1/(1+e-s )s 可使同一网络既能处理小信号，也能处理大信号，所以作为转换函数。累计误差 BP 算法的具体步骤：先给出 P 个训练对（X 1，T 1）（X i，T i）（X p，T p）（1）预置较小的随机权矩阵；（2）施输入模式 Xp于网络，计算各层的输入；Y(s)j=fW(s) jiY(s-1)iY(s)为 s 层上第 j 个神经元的输出值，W(s) ji是(s-1)层第 i 个神经元至第
7、 s 层第 j 个神经元的连接权；（3）修改权值；W(s+1)ji = W(s)ji（t）+e(s) jXtiXti 是该权值的输入信号， 为学习率，t 为学习次数，e(s) j是 s 层第 j 个神经元的误差，对于中间层的神经元： e(s) j=Xi(1-Xi)e(s) kW(s)kj ；对于实际输出:e(s) j=Yi(1-Yi)(tj-yj), tj为第 j 个神经元的期望输出；（4）计算全局误差；E=1/2（t j,p - yj,p） 2= Ep（5）返回第 2 步，向网络加下一个模式对，直到 P 个模式对均循环一遍，再进行第 6 步；（6）若 EEmax(预先设定值)，则停止；否则，
8、令 E=0，返回第 2 步。三、神经元的选取、内符检验与预测效能评估1、神经元的选取神经网络要求选择那些影响输出的主要因素作为输入，各输入变量间应不存在密切的相关性，选定的输入变量数必须足够且具有代表性。输入变量数远远超过网络隐层的节点数，网络就能揭示其蕴藏其间的任意复杂的规律，即具有较好的外推能力。基于这一点我们选取了福建及其沿海东经 117120 度，北纬 2226 度范围内的地震在时间、空间、强度三方面具有代表性的六个指标：频度、蠕变、能量、b 值、缺震、 值指标的 29 组数据作为网络的输入，输出为这一地区第二年可能发生的地震的最大震级。2、迭代次数的确定训练次数直接影响网络的外推能力
9、，神经网络在隐层结点确定的情况下，存在一个最佳的训练次数问题，训练次数过大会出现“训练过头”现象，即训练误差减小，而测试误差反而增大（如图 2，略），外推能力下降。因此我们在训练过程中每经过一定训练次数后，就停止训练并测试其测试误差，当发现测试误差开始上升时，网络便达到最佳训练次数。此时网络便具有最佳的外推能力。3、内符检验与预测效能分析神经元、训练次数已经确定，那么这个网络对预报震级的效能如何，是否能达到要求的精度还必须进行内符检验。内符情况如表 1，可以看出目前我们构造的用于福建及其沿海地震活动预测的神经网络模型基本能识别训练样本，如取|M|0.5，内符检验正确率为 97%。外推能力检验采
10、用交叉法进行，把同样的训练样本分成 10，10，9 三个子集，每次抽出一个子集不参加训练，用它来测试外推能力，分析结果如表 2。可以看出最大预测震级与实际震级之差为0.6，如果限定|M|0.5 为预报正确，则预报震级的准确率为 86%。通过上述内符检验及预测效能分析，我们确定了适用于福建及其沿海这一特定区域地震预报的精度较高的神经网络模型。表 1 学习样本集及内符结果频度 蠕变 能量 B 值 缺震 值 期望输 出 实际输 出 震级 M0.144 0.5879 0.2937 0.71 0.7778 0.598 4.0 4.4 0.40.150 0.3990 0.2314 0.44 1 0.809
11、0 4.1 4.0 -0.10.190 0.7596 0.3680 0.62 0.8283 0.8542 4.1 4.3 0.20.271 0.1030 0.1452 0.59 0.8754 0.9950 3.9 4.0 0.10.230 0.5343 0.2183 0.70 0.7946 1 4.0 4.1 0.10.188 0.6384 0.8646 0.66 0.8114 0.9648 4.8 5.0 0.20.204 0.6697 0.9105 0.62 0.8283 0.9196 4.3 4.6 0.30.157 0.8152 0.5797 0.47 0.9360 0.8693 3.
12、8 4.0 0.20.211 0.7364 0.3865 0.69 0.8182 0.8291 4.4 4.7 0.30.213 0.7808 0.1648 0.63 0.8081 0.8342 4.0 4.2 0.20.126 0.3737 0.1823 0.75 0.7576 0.6935 5.0 4.9 -0.10.125 0.1990 0.1320 0.69 0.6936 0.7437 3.9 4.0 0.10.126 0.2939 0.5895 0.75 0.7374 0.7136 3.7 4.3 0.60.073 0.9414 0.4258 0.84 0.633 0.7940 4.
13、0 4.3 0.30.103 0.1778 0.2074 0.75 0.7104 0.8291 3.8 4.0 0.20.156 0.5768 0.5600 0.70 0.7879 0.7839 5.5 5.1 -0.40.159 0.8515 0.130 0.54 0.8586 0.8995 3.9 4.1 0.20.177 0.2798 0.3111 0.76 0.7240 0.8995 3.9 4.2 0.30.135 0.2859 0.5087 0.71 0.7441 0.8744 4.1 4.0 -0.10.094 0.5687 0.1714 0.55 0.8165 0.9096 3
14、.9 4.0 0.10.211 0.6899 0.4279 0.74 0.8081 0.7337 4.9 4.9 00.120 0.4899 0.3504 0.56 0.8215 0.9899 3.9 3.7 -0.20.154 0.6828 0.0156 0.67 0.7946 0.7889 4.4 4.0 -0.41 1 0.3940 0.70 0.9830 0.6784 5.5 5.5 00.330 0.1222 0.8603 0.68 0.8788 0.7186 4.7 4.6 -0.10.163 0.3162 0.2402 0.87 0.7003 0.6834 4.1 4.5 0.4
15、0.771 0.9121 1 1.03 0.7576 0.7437 5.1 5.1 00.238 0.1091 0.3734 0.69 0.8384 0.7136 4.3 4.1 -0.20.785 0.3512 0.0292 0.85 0.9339 0.7867 4.6 4.6 0表 2 预测效能分析编号 实际震级 检验震级 震级差 M1 4.0 4.3 0.32 4.1 3.6 -0.53 4.1 3.9 -0.24 3.9 3.3 -0.65 4.0 4.3 0.36 4.8 4.4 -0.47 4.3 4.6 0.38 3.8 4.0 0.29 4.4 4.1 0.310 4.0 3.
16、7 -0.311 5.0 4.6 -0.412 3.9 4.5 0.613 3.7 1.3 0.614 4.0 3.6 -0.415 3.8 3.6 -0.216 5.5 5.4 -0.117 3.9 4.1 -0.218 3.9 3.5 -0.419 4.1 4.6 0.520 3.9 4.0 0.421 4.9 4.8 -0.122 3.9 4.3 0.423 4.4 4.0 -0.424 5.5 5.6 0.125 4.7 4.7 026 4.1 3.9 -0.227 5.1 5.4 0.328 4.3 4.0 -0.329 4.6 5.2 0.6四、应用实例：应用该神经网络模型，以 1999 年福建及其沿海地区的地震学指标：频度、蠕变、能量、b值、缺震， 值来作为网络的输入，预测 2000 年的最大震级为 4.70.5，这一结果基本与实际情况相对应。实际上，漳州外海于 2000 年 5 月 27 日发生一次 4.5 级地震；运用 2000 年的数据对 2001 年福建省及
