论文-地震预测模型.docVIP

地震预测模型 摘 要： 本文通过研究图表中的10多个指标，来调查其与地震发生的关系，通过对各指标数据量分析建立相应的模型。 对于问题一，通过数据分析，度量各指标对地震发生的敏感程度。我们利用逻辑回归分析法，设定地震是否发生为0-1变量，建立模型描述因变量（地震是否发生）与自变量（12个指标）之间的线性关系，以此来来得出各个指标的敏感程度。经计算各个指标的敏感程度为：电压倾斜仪EW电磁波辐射EW气氡地温 水温电磁波辐射NS 气温气压倾斜仪NS雨量水位。 对于问题二，通过构建综合指标，反映地震发生前的数据特征的统计量。 我们建立了层次分析法模型来进行求解，由可以求出各指标对地震的权重，再通过可得出综合 对于问题三，通过建立数学模型来研究地震发生前的数量特征。我们应用平均值﹑方差﹑中位数等数量特征来作为模型研究的依据，应用软件excel绘制每次地震发生前后各数量特征的变化情况，并对其进行数据分析。 对于问题四，通过分析各项指标集结地震数据分析平台。我们使用Windows2000server、MS-SQLServer2000数据库地震是地壳快速释放能量过程中造成振动……12 Y 0-1变量，代表地震是否发生 表示十一个相关指标的判断矩阵的元素() ：矩阵B中每行元素 的乘积() 的次方根() 由组成的矩阵的转置 表示各评价指标的权重() 矩阵B的最大特征值 由组成的矩阵的转置 综合指标 三﹑问题分析 题中给出了地震发生前2005年1月1日至2010年6月30日按小时观测的12个指标的数据和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。这些数据中隐藏着地震发生的前兆特征。这12个指标都从不同的侧面反映了地震活动的各项特征。 本题主要是对所给数据进行分析处理，建立适当的数学模型来揭示地震前兆的数据特征。 四﹑模型的建立与求解 问题一： 设地震是否发生为0-1变量，即Y=0（地震不发生），Y=1（地震发生），我们线性回归模型： E(Y)=b0+b1X1+b2X2+…+b12X12 …………(1) 因为Y只能取0或1，而E(Y)的取值是连续的，又E(Y)=P（Y=1）=p,因而我们可用E(Y)来预测Y=1的概率，即： p=a0+a1X1+a2X2+…+a12X12 …………(2) 为了得到各指标的敏感系数，我们通过对P进行Logit变换，即 Logit（P）=ln= 利用SPSS软件得如下两表： 表一 表二 由上述结论我们分析数据得出各个指标在地震发生前的敏感程度分别如下： 电压 电磁波辐度EW 电磁波辐度NS 地温 水位 气温 气压 水温 气氡 雨量 倾斜仪NS 倾斜仪EW 3.922 1.127 0.984 1.000 28499.278 0.926 0.926 1.00 1.006 0.001 0.890 1.647 表三 在有效的敏感系数范围之内，由于考虑到水位的敏感系数太大，我们将不予考虑，也就是说水位对地震没什么影响. 我们将上述指标-与地震关联度从大到小排列： 电压 倾斜仪EW 电磁波辐度EW 气氡 地温 水温 电磁波辐度NS 气温 气压 倾斜仪NS 雨量 3.922 1.647 1.127 1.006 1.000 1.00 0.984 0.926 0.926 0.890 0.001 表四 标准化后有： 电压 倾斜仪EW 电磁波辐度EW 气氡 地温 水温 电磁波辐度NS 气温 气压 倾斜仪NS 雨量 0.292055 0.122645 0.083923 0.074913 0.074466 0.074466 0.073274 0.068955 0.068955 0.066274 7.44657E-05 表五 由以上结果可知各指标敏感程度的顺序如下： 电压倾斜仪EW电磁波辐射EW气氡地温 水温电磁波辐射NS 气温气压倾斜仪NS雨量水位。 问题二： 在反映地震发生前的数据特征的统计规律时，指标的权重对研究统计规律具有重要规律。我们利用AHP来确定权重来对地震指标进行评估。 ( 1) 建立层次结构。根据对问题的初步分析将所包含的因素按总目标、准则层和措施层进行分组, 每一组作为一个层次, 然后以连线表示各层次元素之间的关系, 构成一个从上至下的递阶层次结构( 如下图) , 包括目标层( A层) 、准则层( B层) 、措施层( C层) 共三层。 （注：在B层中，由上问可知，水位对地震的影响很小，可忽略不计） ( 2) 建立判断矩阵。通过两两比较的方法确定各层次中各因子的相对重要性,建立判断矩阵。按九标度法对同一层元素的重要程度进行赋值。 结合已有的资料、专家的意见及自己的认识程度对赋值,