荐言献策

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于GA-GRNN的地震震级预测模型

刘新 尹康达 寇海川

刘新,尹康达,寇海川. 基于GA-GRNN的地震震级预测模型[J]. 华北地震科学,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
引用本文: 刘新,尹康达,寇海川. 基于GA-GRNN的地震震级预测模型[J]. 华北地震科学,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
LIU Xin,YIN Kangda,KOU Haichuan. The Prediction Model of Earthquake Magnitude Based on GA-GRNN[J]. North China Earthquake Sciences,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
Citation: LIU Xin,YIN Kangda,KOU Haichuan. The Prediction Model of Earthquake Magnitude Based on GA-GRNN[J]. North China Earthquake Sciences,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006

基于GA-GRNN的地震震级预测模型

doi: 10.3969/j.issn.1003-1375.2023.02.006
基金项目: 河北省地震科技星火计划项目红山野外站科研专项(DZ202111040001)
详细信息
    作者简介:

    刘新(1978—),男,河北景县人,高级工程师,主要从事地震观测研究和台站运行维护工作. E-mail:hbdsn_lx@163.com

  • 中图分类号: P315.6

The Prediction Model of Earthquake Magnitude Based on GA-GRNN

  • 摘要: 为科学描述地震震级与其敏感因子之间复杂的非线性关系,将遗传算法(Genetic Algorithm,GA)与广义回归神经网络(General Regression Neural Network,GRNN)相结合,利用主成分分析法(Principal Component Analysis,PCA)对地震震级敏感因子进行降维处理,然后对提取出的主成分进行归一化,将归一化的主成分数据作为预测模型的输入向量,地震震级作为预测模型的输出向量;以20个地震数据作为学习样本进行训练,运用GA寻优获得最优光滑因子,建立基于PCA-GA-GRNN的地震震级预测模型,并对8个测试样本进行预测。结果表明:PCA-GA-GRNN模型的最小误差、最大误差和平均误差分别为1.563 0%、4.878 0%和2.647 0%,其平均误差相比于GA-GRNN模型和GRNN模型分别降低5.666 7%和5.026 4%,具有较高的预测精度。
  • 图  1  遗传算法流程图

    Figure  1.  The flow chart of GA

    图  2  PCA-GA-GRNN模型流程图

    Figure  2.  The flow chart of PCA-GA-GRNN

    图  3  GA寻优

    Figure  3.  GA optimization

    图  4  学习样本回判结果

    Figure  4.  Learning sample feedback results

    表  1  归一化地震数据

    Table  1.   Normalized seismic data

    序号地震累积频度累积释放能量平均震级ηb相关区震级实际震级
    10.139 5000.098 60.390 800
    20.085 60.000 10.166 70.251 60.630 50.250.138 6
    30.090 20.000 60.50.203 10.501 40.350.310 8
    40.230 80.000 600.528 60.947 50.100.345 2
    500.000 80.333 30.490 30.240 10.450.344 8
    60.227 60.002 30.50.118 90.100 30.300.447 6
    70.140 4000.068 70.652 00.050.072 0
    80000.323 00.791 80.050
    90.174 20.005 80.50.498 60.727 50.300.310 3
    100.086 20.000 60.50.205 20.510 30.350.310 8
    110.968 30.001 70.666 700.985 710.344 8
    1200.001 9100.400 50.350.448 3
    130.160 30.110 50.666 70.650 10.084 210.793 1
    140000.319 50.810 60.050
    150.079 9000.321 80.809 20.050.140 2
    160.219 30.005 00.50.208 20.438 50.550.380 2
    170.908 50.981 40.833 30.401 20.620 30.801
    180.053 10.074 20.666 70.783 100.950.758 6
    190.168 70.000 50.50.051 20.740 10.300.310 3
    200.040 20.000 600.987 10.240 20.300.344 8
    210.099 60.000 70.50.013 10.837 60.300.345 2
    220.140 40.019 50.333 30.648 50.330 10.500.655 2
    230.227 60.000 700.532 10.953 10.100.345 2
    240.158 30.000 90.50.013 20.850 10.300.379 3
    250.230 50.001 70.50.121 50.962 10.300.448 3
    260.460 30.003 20.333 30.330 50.621 30.550.482 8
    270.140 4000.075 90.641 70.050.069 0
    280.140 4000.101 30.389 000
    下载: 导出CSV

    表  2  相关系数阵

    Table  2.   Correlation coefficient matrix

    累积释放能量地震累积频度平均震级ηb相关区震级
    累积释放能量10.5940.3670.140.0430.384
    地震累积频度0.59410.355−0.1430.2860.526
    平均震级0.3670.3551−0.14−0.2320.743
    η0.14−0.143−0.141−0.3960.283
    b−0.0430.286−0.232−0.3961−0.352
    相关区震级0.3840.5260.7430.283−0.3521
    下载: 导出CSV

    表  3  归一化特征值、贡献率和累计贡献率

    Table  3.   Normalized eigenvalue, contribution rate and cumulative contribution rate

    成分归一化特征值贡献率/%累计贡献率/%
    Y12.52842.1442.14
    Y21.61426.969.04
    Y30.95715.94484.985
    Y40.5268.75993.744
    Y50.3145.23798.981
    Y60.0611.019100
    下载: 导出CSV

    表  4  三种模型预测结果对比

    Table  4.   Comparison of prediction results of three models

    编号实际值GRNN模型GA-GRNN模型PCA-GA-GRNN模型
    预测值相对误差/%预测值相对误差/%预测值相对误差/%
    215.45.87.407 45.15.556 05.31.852 0
    225.24.95.769 05.67.692 35.03.846 0
    237.07.45.714 36.57.143 07.22.857 1
    245.14.77.843 05.45.882 45.01.961 0
    256.46.06.250 06.86.250 06.31.563 0
    265.05.48.000 05.48.000 05.12.000 0
    274.14.612.195 14.59.756 14.34.878 0
    284.55.113.333 35.011.111 14.62.222 2
    下载: 导出CSV
  • [1] 李龙师, 金波, 冯志仁. 改进的层次分析法与模糊综合评价在评估水库诱发地震震级上限中的应用[J]. 世界地震工程, 2020, 36(3): 223-232.
    [2] 吴清, 高孟潭. 利用烈度数据点估算中强历史地震震级与震中的方法研究[J]. 震灾防御技术, 2014, 9(1): 12-28. doi:  10.3969/j.issn.1673-5722.2014.01.002
    [3] 刘家豪, 朱国维, 周俊杰. 舟坝水库诱发地震震级预测[J]. 高原地震, 2020, 32(3): 13-20.
    [4] 苏义鑫, 沈俊, 张丹红, 等. 神经网络和改进粒子群算法在地震预测中的应用[J]. 计算机应用, 2011, 31(7): 1793-1796,1807.
    [5] 陈一超, 曾三友, 张好春, 等. 基于遗传神经网络的地震预测研究[J]. 计算机应用与软件, 2008, 25(4): 135-137.
    [6] 项月文, 饶泓, 汤兰荣, 等. 基于GA-BP神经网络的地震预报技术[J]. 国际地震动态, 2016(8): 4-10.
    [7] 韩晓飞, 潘存英, 罗词建. 基于遗传算法的广义回归神经网络在地震预测中的应用[J]. 华北地震科学, 2012, 30(1): 48-53. doi:  10.3969/j.issn.1003-1375.2012.01.010
    [8] 王威, 马东辉, 苏经宇, 等. 利用GRNN方法分析地下管道抗震性能影响因素的敏感性[J]. 地震工程与工程振动, 2009, 29(2): 138-142.
    [9] 周德红, 冯豪, 程乐棋, 等. 遗传算法优化的BP神经网络在地震死亡人数评估中的应用[J]. 安全与环境学报, 2017, 17(6): 2267-2272.
    [10] 郭小东, 田杰, 王威, 等. 基于GA-SVR的建筑物液化震陷预测方法[J]. 北京工业大学学报, 2011, 37(6): 829-835. doi:  10.11936/bjutxb2011060829
    [11] 胡瑞卿, 王彦春, 尹志恒, 等. 结合CEEMDAN和主成分分析的低信噪比微地震初至信号检测[J]. 石油地球物理勘探, 2019, 54(1): 45-53.
    [12] 邹冠贵, 任珂, 吉寅, 等. 基于主成分分析和最近邻算法的断层识别研究[J]. 煤田地质与勘探, 2021, 49(4): 15-23.
    [13] 宫凤强, 李嘉维. 基于PCA-DDA原理的砂土液化预测模型及应用[J]. 岩土力学, 2016, 37(S1): 448-454.
    [14] 杨波, 刘泽民, 隆爱军, 等. 遗传算法在提高非包围地震定位精度中的应用[J]. 华北地震科学, 2017, 35(2): 43-47.
    [15] 丁风和, 高立新. 用遗传算法反演内蒙古中西部地区的Q值和场地响应[J]. 华北地震科学, 2005, 23(1): 22-26. doi:  10.3969/j.issn.1003-1375.2005.01.004
    [16] 王立琦, 姚静, 王睿莹, 等. 基于PLS-GRNN的豆粕品质近红外光谱检测研究[J]. 光谱学与光谱分析, 2022, 42(5): 1433-1438. doi:  10.3964/j.issn.1000-0593(2022)05-1433-06
    [17] 李浩, 王晓原, 韩俊彦, 等. 基于高德导航数据与FOA-GRNN模型的驾驶倾向性辨识方法[J]. 交通信息与安全, 2022, 40(2): 63-72. doi:  10.3963/j.jssn.1674-4861.2022.02.008
    [18] 李新春, 谷永延, 黄朝晖, 等. 基于高阶累积量和改进GRNN的CSI手臂行为识别[J]. 重庆邮电大学学报:自然科学版, 2022, 34(2): 331-340.
    [19] 张研, 邝贺伟. 地震震级预测的相关向量机模型[J]. 世界地震工程, 2020, 36(1): 212-221.
    [20] 何冰琛, 杨薛明, 王劲松, 等. 基于PCA-GPR的锂离子电池剩余使用寿命预测[J]. 太阳能学报, 2022, 43(5): 484-491. doi:  10.19912/j.0254-0096.tynxb.2022-0422
    [21] 张晓阳, 侯博超. 基于PCA-SVM模型的饲料企业财务精细化管理研究[J]. 饲料研究, 2022, 45(10): 122-127. doi:  10.13557/j.cnki.issn1002-2813.2022.10.026
    [22] 乔守旭, 钟文义, 谭思超, 等. 基于PCA-GA-SVM的竖直下降两相流流型预测[J]. 核动力工程, 2022, 43(3): 85-93. doi:  10.13832/j.jnpe.2022.03.0085
  • [1] 王嘉琦, 周剑青, 王曰风, 袁国旭, 常玉柱.  同址观测洞体倾斜仪器背景噪声对比 . 华北地震科学, 2023, 41(2): 90-97. doi: 10.3969/j.issn.1003-1375.2023.02.013
    [2] 黄惠宁, 阎春恒, 向巍.  公众对地震预报认识和理解程度调查分析 . 华北地震科学, 2022, 40(1): 106-111, 121. doi: 10.3969/j.issn.1003-1375.2022.01.017
    [3] 王晨晖, 袁颖, 刘立申, 陈凯男, 吴鹤帅.  基于PCA-GSM-SVM的地震伤亡人数预测 . 华北地震科学, 2019, 37(3): 25-30. doi: 10.3969/j.issn.1003-1375.2019.03.004
    [4] 张亮娥, 闫计明, 陈常俊, 刘素珍, 贾海英, 田文清, 曹文强.  主成分分析在山西地磁台阵数据分析中的应用 . 华北地震科学, 2018, 36(1): 64-68. doi: 10.3969/j.issn.1003-1375.2018.01.011
    [5] 杨波, 刘泽民, 隆爱军, 谢石文, 韩成成.  遗传算法在提高非包围地震定位精度中的应用 . 华北地震科学, 2017, 35(2): 43-47. doi: 10.3969/j.issn.1003-1375.2017.02.007
    [6] 韩晓飞, 潘存英, 罗词建.  基于遗传算法的广义回归神经网络在地震预测中的应用 . 华北地震科学, 2012, 30(1): 48-53.
    [7] 李强, 徐戈, 杨彩霞.  人工神经网络在潜在地震危险区估计中的应用 . 华北地震科学, 2007, 25(2): 36-41.
    [8] 丁风和, 高立新.  用遗传算法反演内蒙古中西部地区的Q值和场地响应 . 华北地震科学, 2005, 23(1): 22-26.
    [9] 俞言祥, 许忠淮, 黄雨蕊.  用斜井钻孔崩落资料反演上地壳应力状态的遗传算法 . 华北地震科学, 2000, 18(1): 1-8.
    [10] 冯德益, 郑熙铭, 汪德馨, 蒋淳, 田山, 罗兰格.  神经网络在短期地震跟踪预报中的应用 . 华北地震科学, 1995, 13(2): 16-25.
  • 加载中
图(4) / 表 (4)
计量
  • 文章访问数:  1976
  • HTML全文浏览量:  970
  • PDF下载量:  15
出版历程
  • 收稿日期:  2022-06-15
  • 网络出版日期:  2023-04-20
  • 刊出日期:  2023-04-30

基于GA-GRNN的地震震级预测模型

doi: 10.3969/j.issn.1003-1375.2023.02.006
    基金项目:  河北省地震科技星火计划项目红山野外站科研专项(DZ202111040001)
    作者简介:

    刘新(1978—),男,河北景县人,高级工程师,主要从事地震观测研究和台站运行维护工作. E-mail:hbdsn_lx@163.com

  • 中图分类号: P315.6

摘要: 为科学描述地震震级与其敏感因子之间复杂的非线性关系,将遗传算法(Genetic Algorithm,GA)与广义回归神经网络(General Regression Neural Network,GRNN)相结合,利用主成分分析法(Principal Component Analysis,PCA)对地震震级敏感因子进行降维处理,然后对提取出的主成分进行归一化,将归一化的主成分数据作为预测模型的输入向量,地震震级作为预测模型的输出向量;以20个地震数据作为学习样本进行训练,运用GA寻优获得最优光滑因子,建立基于PCA-GA-GRNN的地震震级预测模型,并对8个测试样本进行预测。结果表明:PCA-GA-GRNN模型的最小误差、最大误差和平均误差分别为1.563 0%、4.878 0%和2.647 0%,其平均误差相比于GA-GRNN模型和GRNN模型分别降低5.666 7%和5.026 4%,具有较高的预测精度。

English Abstract

刘新,尹康达,寇海川. 基于GA-GRNN的地震震级预测模型[J]. 华北地震科学,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
引用本文: 刘新,尹康达,寇海川. 基于GA-GRNN的地震震级预测模型[J]. 华北地震科学,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
LIU Xin,YIN Kangda,KOU Haichuan. The Prediction Model of Earthquake Magnitude Based on GA-GRNN[J]. North China Earthquake Sciences,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
Citation: LIU Xin,YIN Kangda,KOU Haichuan. The Prediction Model of Earthquake Magnitude Based on GA-GRNN[J]. North China Earthquake Sciences,2023, 41(2):37-42. doi:10.3969/j.issn.1003−1375.2023.02.006
    • 地震震级敏感因子众多,单一的数学公式很难表达这种复杂的非线性关系,近年来,许多专家学者进行了大量的探索。李龙师等[1]利用改进的层次分析法和模糊评判法评估了水库诱发地震震级上限;吴清等[2]采用统计学原理对烈度估计线进行统计回归并建立烈度分布模型,利用该模型对历史震级进行了估算;刘家豪等[3]利用类比法、概率统计检验法、综合影响参数预测法等对水库诱发地震震级进行了整体和分区预测。以上方法均取得了良好的效果,但模型预测依赖主观经验,导致预测精度较低。针对上述问题,神经网络等机器学习算法应运而生,并迅速应用于地震震级预测。苏义鑫等[4]提出了一种基于神经网络与改进粒子群算法的地震预测方法;陈一超等[5]将遗传算法和具有深度搜索能力的BP算法相结合实现地震震级预测;项月文等[6]构建GA-BP神经网络模型,并以东南沿海地区的地震目录作为研究样本进行了预测。而BP神经网络易陷入局部最优,还需探索更为有效的科学方法。

      广义回归神经网络(General Regression Neural Network,GRNN)是典型的前馈式局部逼近神经网络,运行效率高,且不易陷入局部最优。韩晓飞等[7]利用GRNN选取了100组样本数据进行仿真实验;王威等[8]基于广义回归神经网络建立了地下管道地震反应预测模型,并对地下管道抗震性能影响因素进行敏感性分析,均取得了良好效果,但GRNN模型参数不易收敛,还需探索更加科学的寻优方法。遗传算法具有较强的全局搜索能力,周德红等[9]利用遗传算法对BP神经网络的权值与阀值进行优化,建立了地震死亡人数预测模型。郭小东等[10]利用遗传算法自动确定ν-SVR的最优模型参数,建立了建筑物液化震陷量与其各种影响因素之间的非线性关系,表明遗传算法能够改善模型的寻优效果。因此,本文选取遗传算法全局搜索GRNN最优参数,同时为了进一步降低模型整体复杂度,运用主成分分析法对地震数据进行降维处理,该方法已广泛应用于地震信号检测[11]、断层识别[12]和地震砂土液化预测[13]等。基于此,本文建立了基于PCA-GA-GRNN的地震震级预测模型,首先利用地震震级敏感因子进行降维处理,然后运用GA寻优获得GRNN最优光滑因子,并对8个测试样本进行预测,取得了一定的效果。

    • 遗传算法[14-15]是一种汲取自然进化基本思想寻求最优解的方法,有较强的全局优化特性和全局搜索能力。主要步骤如图1所示:

      图  1  遗传算法流程图

      Figure 1.  The flow chart of GA

    • GRNN理论分析如下[16-18],设随机变量xy的观测值分别为XY,其联合概率密度函数为fx,y),则Y的预测输出$ \hat {\boldsymbol Y} $为:

      $${\boldsymbol Y}_{id}^{k + 1} = {\boldsymbol Y}_{id}^k + {\boldsymbol V}_{id}^{k + 1} $$ (1)
      $$ {\hat {\boldsymbol Y}} = E(y/{\boldsymbol X}) = \frac{{\displaystyle\int_{ - \infty }^{ + \infty } {yf({\boldsymbol X},y){\rm d}y} }}{{\displaystyle\int_{ - \infty }^{ + \infty } {f({\boldsymbol X},y){\rm d}y} }} $$ (2)

      假设fx,y)服从正态分布,则有:

      $$ \begin{aligned} {\hat f} ({\boldsymbol X},y) =& \frac{1}{{n{{(2{\text{π )}}}^{\frac{{p + 1}}{2}}}{\sigma ^{p + 1}}}}\sum\limits_{i = 1}^n {\exp } \\& \left[ { - \frac{{{{({\boldsymbol X} - {X_i})}^{\text{T}}}({\boldsymbol X} - {X_i})}}{{2{\sigma ^2}}}} \right]\exp \left[ { - \frac{{{{({\boldsymbol X} - {Y_i})}^2}}}{{2{\sigma ^2}}}} \right] \\ \end{aligned} $$ (3)

      式中:$ {X_i} $$ {Y_i} $分别为xy的样本观测值;X为输入变量;n为样本数量;p为随机变量x的维度;$ \sigma $为光滑因子。将$ {\hat f} ({\boldsymbol X},y) $代替fx,y)代入式(3)得下式:

      $$ {\hat {\boldsymbol Y}} ({\boldsymbol X}) = \frac{{\displaystyle\sum\nolimits_{i = 1}^n {{{\boldsymbol Y}_i}\exp \left[ { - \frac{{{{(X - {X_i})}^{\text{T}}}(X - {X_i})}}{{2{\sigma ^2}}}} \right]} }}{{\displaystyle\sum\nolimits_{i = 1}^n {\exp \left[ { - \frac{{{{(X - {X_i})}^{\text{T}}}(X - {X_i})}}{{2{\sigma ^2}}}} \right]} }} $$ (4)

      GRNN的网络结构由以下4个步骤实现。

      1)输入层。该层神经元个数对应输入变量xi即主成分个数,本文将PCA降维得到的4个主成分作为GRNN模型输入变量,因此输入层神经元个数为4。

      2)模式层。该层神经元个数等于训练样本的数量。每个神经元传递函数为:

      $$ {P_i} = \exp \left[ { - \frac{{{{(X - {X_i})}^{\text{T}}}(X - {X_i})}}{{2{\sigma ^2}}}} \right]{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} i = 1,2, \cdots ,n $$ (5)

      3)求和层。该层由2种求和方式,一种是算数求和$ {S_D} $,另一种是加权求和$ {S_{Nj}} $,传递函数分别为:

      $$ {S_D} = \sum\nolimits_{i = 1}^n {{P_i}} {\kern 1pt} $$ (6)
      $$ {\kern 1pt} {\kern 1pt} {\kern 1pt} {S_{Nj}} = \sum\nolimits_{i = 1}^n {{y_{ij}}{P_i}{\kern 1pt} } $$ (7)

      式中:$ {P_i} $为模式层第i个神经元的传递函数;$ {y_{ij}} $为模式层第i个神经元与求和层第j个分子求和神经元连接权值。

      4)输出层。该层神经元个数只有1个,即地震震级。将$ {S_{Nj}} $$ {S_D} $相除,可得输出结果。

      $$ {y_i} = \frac{{{S_{Nj}}}}{{{S_D}}} $$ (8)
    • 基于PCA-GA-GRNN的地震震级预测模型构建流程图如图2所示。

      图  2  PCA-GA-GRNN模型流程图

      Figure 2.  The flow chart of PCA-GA-GRNN

    • 通过查阅文献并结合前人经验选取文献[19]的地震数据作为研究样本,随机选取20个地震样本作为学习样本,剩余8个样本作为测试样本,地震样本归一化数据如表1所示,其中敏感因子为:半年内M≥3的地震累计频度、半年内能量释放积累值、b值(震级和频度关系式中的比例系数)、平均震级(某一次地震中主震和余震的平均值)、η值(震级-频度关系偏离G-R关系式的程度)及相关地震区地震震级(某地震活动区特定时间段内中所发生的除主震之外的最大地震震级)。这些因子反映了地震活动在时空和强度2个方面的特征,具有较好的中期或中短期预测效果。

      表 1  归一化地震数据

      Table 1.  Normalized seismic data

      序号地震累积频度累积释放能量平均震级ηb相关区震级实际震级
      10.139 5000.098 60.390 800
      20.085 60.000 10.166 70.251 60.630 50.250.138 6
      30.090 20.000 60.50.203 10.501 40.350.310 8
      40.230 80.000 600.528 60.947 50.100.345 2
      500.000 80.333 30.490 30.240 10.450.344 8
      60.227 60.002 30.50.118 90.100 30.300.447 6
      70.140 4000.068 70.652 00.050.072 0
      80000.323 00.791 80.050
      90.174 20.005 80.50.498 60.727 50.300.310 3
      100.086 20.000 60.50.205 20.510 30.350.310 8
      110.968 30.001 70.666 700.985 710.344 8
      1200.001 9100.400 50.350.448 3
      130.160 30.110 50.666 70.650 10.084 210.793 1
      140000.319 50.810 60.050
      150.079 9000.321 80.809 20.050.140 2
      160.219 30.005 00.50.208 20.438 50.550.380 2
      170.908 50.981 40.833 30.401 20.620 30.801
      180.053 10.074 20.666 70.783 100.950.758 6
      190.168 70.000 50.50.051 20.740 10.300.310 3
      200.040 20.000 600.987 10.240 20.300.344 8
      210.099 60.000 70.50.013 10.837 60.300.345 2
      220.140 40.019 50.333 30.648 50.330 10.500.655 2
      230.227 60.000 700.532 10.953 10.100.345 2
      240.158 30.000 90.50.013 20.850 10.300.379 3
      250.230 50.001 70.50.121 50.962 10.300.448 3
      260.460 30.003 20.333 30.330 50.621 30.550.482 8
      270.140 4000.075 90.641 70.050.069 0
      280.140 4000.101 30.389 000
    • 利用PCA[20-22]对6个敏感因子进行相关性分析和降维处理,相关系数矩阵见表2,归一化特征值,各成分贡献率和累计贡献率见表3。由表3可以发现前4个主成分的累计贡献率为93.744%,涵盖了原始数据的主要信息,其中4个主成分的表达式见式(9)。

      表 2  相关系数阵

      Table 2.  Correlation coefficient matrix

      累积释放能量地震累积频度平均震级ηb相关区震级
      累积释放能量10.5940.3670.140.0430.384
      地震累积频度0.59410.355−0.1430.2860.526
      平均震级0.3670.3551−0.14−0.2320.743
      η0.14−0.143−0.141−0.3960.283
      b−0.0430.286−0.232−0.3961−0.352
      相关区震级0.3840.5260.7430.283−0.3521

      表 3  归一化特征值、贡献率和累计贡献率

      Table 3.  Normalized eigenvalue, contribution rate and cumulative contribution rate

      成分归一化特征值贡献率/%累计贡献率/%
      Y12.52842.1442.14
      Y21.61426.969.04
      Y30.95715.94484.985
      Y40.5268.75993.744
      Y50.3145.23798.981
      Y60.0611.019100
      $$ \begin{aligned} {Y_1} = &0.281S{{1}} + 0.279S{{2}} + 0.314S{{3}} + 0.063S{\text{4}} +\\& 0.155S{{5 - }}0.101S{\text{6}} \\ {Y_2} = &0.11S{\text{1}} + 0.341S{\text{2 + }}0.007S{{3}} - 0.452S{{4 + }}\\&0.518S{{5}} - 0.135S{\text{6}} \\ {Y_3} =& 0.47S{\text{1 + }}0.268S{\text{2 - }}0.55S{\text{3}} + 0.615S{\text{4 + }}\\&0.23S{{5 - }}0.14S{\text{6}} \\ {Y_4} =& - 0.927S{\text{1 + }}0.285S{\text{2 - }}0.015S{\text{3}} + 0.488S{\text{4 + }}\\&0.595S{\text{5}} + 0.609S{\text{6}} \\ \end{aligned} $$ (9)

      其中:S1为累积释放能量;S2为地震累积频度;S3为平均震级;S4为η值;S5为b值;S6为相关区震级。

    • 采用遗传算法确定GRNN的光滑因子,设置种群数量n=20,终止代数为100,GA寻优过程见图3,确定GRNN的最优光滑因子δ=3.457,并对学习样本进行回判检验(图4)。

      图  3  GA寻优

      Figure 3.  GA optimization

      图  4  学习样本回判结果

      Figure 4.  Learning sample feedback results

      图4可以看出,预测准确率平均水平达89.2%,因此该模型可以用于测试样本的预测。

    • 为验证模型的适用性,在保证相同学习样本的前提下,PCA-GA-GRNN采用降维后的主成分数据作为输入向量;而GSM-GRNN模型和GRNN模型的学习样本未经PCA降维处理。3种模型同时对测试样本进行预测,利用预测值与实际值的绝对误差与真实值之比乘以100%求出相对误差,所得的预测结果见表4

      表 4  三种模型预测结果对比

      Table 4.  Comparison of prediction results of three models

      编号实际值GRNN模型GA-GRNN模型PCA-GA-GRNN模型
      预测值相对误差/%预测值相对误差/%预测值相对误差/%
      215.45.87.407 45.15.556 05.31.852 0
      225.24.95.769 05.67.692 35.03.846 0
      237.07.45.714 36.57.143 07.22.857 1
      245.14.77.843 05.45.882 45.01.961 0
      256.46.06.250 06.86.250 06.31.563 0
      265.05.48.000 05.48.000 05.12.000 0
      274.14.612.195 14.59.756 14.34.878 0
      284.55.113.333 35.011.111 14.62.222 2

      表4可以看出,GRNN模型的最小误差、最大误差和平均误差分别为5.714 3%、13.333 3%和8.314 0%;GA-GRNN模型的最小误差、最大误差和平均误差分别为5.882 4%、11.111 1%和7.674 0%;PCA-GA-GRNN模型的最小误差、最大误差和平均误差分别为1.563 0%、4.878 0%和2.647 0%。因此PCA-GA-GRNN模型预测结果平均误差相比于GA-GRNN模型和GRNN模型分别降低5.666 7%和5.026 4%,同时PCA-GA-GRNN模型由于进行了主成分降维处理和遗传算法的优化,能够明显提高运行效率,降低模型复杂度。

    • 1) 采用PCA对影响地震震级的6个敏感因子进行降维处理,模型输入维度由6维降至4维,使模型运行效率明显提高,模型结构得以简化。

      2) 利用GA寻优得到最优光滑因子δ=3.457,建立了基于PCA-GA-GRNN的地震震级预测模型,并对测试样本进行预测,同时与GRNN模型和GA-GRNN模型的预测结果进行对比,结果表明PCA-GA-GRNN模型预测结果平均误差相比于GA-GRNN模型和GRNN模型分别降低5.666 7%和5.026 4%,可以为地震震级预测提供参考。

参考文献 (22)

目录

    /

    返回文章
    返回