2022年地震模型 .pdf

1 地震预测模型摘要地震预报是当代自然科学领域里一个直接关系到人民生命安全和社会经济的发展，是科学难度很大的前沿课题。近几十年来， 人民的努力虽然取得了丰硕的成果，预报的实践有些有限的成功。在地震发生前有很多前兆性指标，如磁场强度、 氡值、电压、雨量、水温等，这些指标都从不同的侧面反映了地震活动的各项特征。在正常情况下， 也常常有些指标有明显的异常， 而令一些指标并不出现异常。 这些都给实际的前兆指标数据特征分析以及地震预报工作带来了困难。本文针对地震影响因素多， 数据多， 联系弱的特点， 分别建立了针对各个任务的数学模型， 首先，处理初始数据时选取日平均变化来消除一天中数据的随机因素的影响，然后更具这些数据建立了地震发生前后各个指标的不同程度度量模型解决了地震对指标的敏感度的度量问题，进而找出了EW ，气氡，水位，电压，雨量这些衡量地震的主要因素， 并分析出了这些指标的重要程度，在解决找统计量时，利用上次任务所得的主要因素的组合来衡量一次地震的地震前兆指标数量特征，并得出了电压达到26.8754，气氡达到 17.685 左右时就有地震发生的可能性，考虑全面周到， 效果较好，在模型改进中所提出的判别分析法，科学有效，对数据利用率较高。 而除此之外的其他几个指标与震级有关，这样就这些指标的分析在任务三中分别建立了地磁前兆异常的动态从属函数模型以及非线性综合模型，进而找出了地震级数与这些因素之间的非线性关系，对于上述数学方法在地震预测预报中的应用没有达到预期效果，须进一步研究改进。综合这些工作，我们应用易语言编写了分析地震数据的平台，通过分析地震数据生成了报告。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 20 页 - - - - - - - - - 2 一、 问题重述地震是地壳快速释放能量过程中造成的振动。虽然预测地震是世界性难题，但迄今科学界普遍认为，有可能反映地震前兆特征的指标可能不少于10 个。已经有专业仪器在多个定点实时按秒记录这些指标的数据，期望通过对记录数据的分析研究找到地震的前兆特征。现已采集到某地 2005年 1 月 1 日至 2010年 6 月 30 日按小时观测的 10 多个指标的数据， 和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。 这些数据中隐藏着地震发生的前兆特征。科学地截取这些数据的有用片段，对数据进行合理地预处理， 用数学方法揭示地震前兆的数据特征，是一项很有意义的研究工作。题给数据中的这 10 多个指标，究竟哪些与地震的发生有关，有何种关系，是单一关系还是复合关系；除这10 多个指标外还有哪些因素及含题给指标在内的哪些指标的哪种数学模型更能反映地震的前兆特征等等，人们迄今仍不很清楚，需要进行深入地研究。 地震数据的观测是持续进行的，随着时间的推移数据的规模会不断扩大。 从中挖掘地震的前兆特征， 必须有合理的数学模型， 也必须有科学高效的算法分析平台。因此，请参赛队结合附件中给出的实际记录数据，尝试完成以下任务。任务一：分析数据特征， 建立数学模型以度量各指标对地震发生的敏感程度。任务二：构造由某些或全部指标构成的综合指标，使其尽可能地集中反映地震发生前的数据特征的统计规律。任务三：结合题给数据， 广泛查阅与地震相关的其它指标的数据和分析方法，建立数学模型来研究地震发生前的数量特征。任务四：将前述各项任务的计算程序集结成地震数据分析平台，使其能够完成形如题给数据的其它地震数据的分析，并能自动输出前述任务的重要的分析结果。任务五：对于进一步的研究设想写一篇切实可行的报告。二、问题分析地震是地下岩层受应力作用错动破裂造成的地面震动，是一种破坏性极强的自然灾害，是自然灾害之首恶。地震可以在很短的时间内使一座城市夷为平地，使无数的家庭支离破碎。 为了有效地减小地震带来的损失，对地震进行预测则很重要。本题中给出了地震发生前2005 年 1 月 1 日至 2010 年 6月 30 日按小时观测的 10 多个指标的数据和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。 这些数据中隐藏着地震发生的前兆特征。在地震发生前有很多前兆性指标，如磁场强度、氡值、电压、雨量、水温等，这些指标都从不同的侧面反映了地震活动的各项特征。 在正常情况下， 也常常有些指标有明显的异常，而令一些指标并不出现异常。 这些都给实际的前兆指标数据特征分析以及地震预报工作带来了困难。21、问题的性质本文主要是对所给数据进行处理，用数学方法揭示地震前兆的数据特征。2.2 、解决问题的难点和关键名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 20 页 - - - - - - - - - 3 1、如何确定前兆个指标对地震发生的敏感程度。2、如何构造地震发生的前兆综合指标，尽可能的揭示地震发生前的数据特征的统计规律。3、构建地震数据分析平台。三、模型假设1、假设地震发生前2005年 1 月 1 日至 2010 年 6 月 30 日按小时观测的 10 多指标的数据是准确无误的。2、 假设由于其他因素而引起10 多个指标数据的变化以及非正常波动可以忽略不计。3、假设第二问剔除由第一问求出的敏感度较低的指标对地震的敏感度可以忽略不计。4、假设地震的前兆指标的数据特征符合一定的概率统计分布。四、符号说明1( )it第i次地震的地震前的检测指标随时间的变化规律2( )it第i次地震的地震后的检测指标随时间的变化规律()2( )jt第i次地震的地震后的预测指标随时间的变化规律D地震级别组成的集合E地震前兆测量指标集合ig各因素所占的权重地震前兆预测异常特征量r观测值的自相关系数K观测曲线( )Y t对应于时间轴的斜率五、模型的建立与求解5.1 数据的处理与分析5.1.1 、各指标数据的分析在地震活动性分析中有许多前兆性的数据指标，例如题中数据给出的氡值、水位、磁场强度等因素。 这些指标在不同的侧面反映了地震活动的各种特征。但在实际的预报中， 常常有些参数在一些中强以上的地震前出现比较明显。在正常情况下，也常常有些参数出现较明显的异常，而另一些参数并不出现异常。这些都给实际预报带来了困难。5.1.2 数据的预处理名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 20 页 - - - - - - - - - 4 由于题目中给出地震每小时测数据受其他因素受天气、气候等其他外在因素的影响造成的数据波动较大， 为了剔除偶然因素的影响， 本文采取了求日平均的方法来对数据进行预处理。如下图分别列举了2005/2006的 EW 人均值图像。由于篇幅有限，其他年份指标的图像都附在附录里面。图 1： 2005 年 EW 日平均变化趋势图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 20 页 - - - - - - - - - 5 图 2： 2006 年 EW 日平均变化趋势图本数据预处理使用了均值结构模型消除了其他偶然因素的影响，使地震前兆指标数据更加准确的反应地震发生的前兆信息。5.2 模型的建立与求解5.2.1 指标敏感度模型设变量 x的变化规律如下所示：110201() ,()(),ttttxttttt（1）求0t前后1( ) t与2( ) t的不同敏感程度。现已知一函数序列( )2( )it ，01,ttt表示第i年()ij的2( ) t 的变化规律，现在结合( )1( )jt 来求( )2( )jt ，现在若( )2( )jt 已知，如何求出()2( )jt 与()2( )jt 的关联程度。若关联程度较大，则( )2( )jt 与( )1( )it 不同程度表示：即定义( )1( )jt 与( )2( )jt 的不同程度为：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 20 页 - - - - - - - - - 6 ( )( )221:jjp（2）分析变量1X 与2X 的关联程度121212(,)()()()Cov XXE XXE XE X( )( )221(,)jjPCov（3）则由分析可得该问题便为( )2j如何求的问题。现已知( )2j()ij与( )1( )jt ，如何根据已知的以上信息求出( )2j，已知( )1( ),()itij 的变化规律，可知令( )( )2( )21()1()( )( ),( )ijijEtE EtEt()ij（4）模型求解：由上述结论我们分析数据得出各个因素在地震发生前后的敏感程度分别如下：电压EW NS 地温水位气温气压水温气氮雨量0.2 1.8 0.8 0.1 2595.3 2.3 4.3 40.2 2.8 0.5 表 1：各因素与地震敏感度的表但是考虑地震因素影响上述指标外，可能还有其它的因素影响上述指标。其中查阅资料可得出：水位，水温，气压，气温受到外界影响最大，这样我们除上述几个指标剩下的指标基于与地震关联度从大到小排列：气氮EW NS 雨量电压地温2.8 1.8 0.8 0.5 0.2 0.1 表 2：主要地震前兆指标敏感度的表标准化后有：气氮EW NS 雨量电压地温0.45 0.29 0.13 0.08 0.03 0.01 表 3：主要地震前兆指标敏感度归一化表如下图分别画了敏感度较高的EW 与敏感度较低的低温年平均曲线，在途中可以看出在前兆指标因素发生剧烈变化时，会有地震的发生。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 20 页 - - - - - - - - - 7 图 3： 2005，2006 年 EW 日平均变化趋势图图 4： 2005、2006 年地温日平均变化趋势图5.2.2 综合指标模型方法及步骤用两个集合， 一是地震级别组成的集合， 用D表示，二是地震前兆车辆指标集合，用E来表示，且有12,mDD DD，12,mEE EE，每个指标因素都有 m个状态级，如此需要对个指标影响程度分敏感、一般、不敏感、很不敏感五名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 20 页 - - - - - - - - - 8 个状态。有n个指标因素分别用12,nU UU 来表示，其中12,TiiiimUUUU，于是可得到敏感指标分析数学模式1 12 111 22 2211 2,nnmn mUUUUUUDUUU（5）有了这个模型，分别给予地震级别的隶属度函数值(1,2,)iP im 再根据指标因素间的关联度及其重要性，分配权数(1,2,)ig in这样就得到了集合D 上的模糊关系：112111222212,.,nnmmnmg P g Pg Pg P g Pg PUg Pg Pg P（6）基于地震的综合指标步骤有下列步骤：（1） 、 详细分析影响地震发生的因素， 并筛选出若干重要因子作为分析指标。（2） 、在分析地震各指标以及波形的基础上，综合得到地震的等级指标。（3） 、给出适当的隶属度iP 和权数ig 。（4） 、修正矩阵 U。确定分析指标地震发生前的检测指标较多， 而各个指标的敏感程度不同， 经过第一问的各因素的关联度分析，筛选出气氮、EW、NS 等六个指标。隶属度iP 和权数ig 的确定应用综合指标进行判别， 隶属度函数是一个关键， 用它来反映地震发生。 设N是评定集合上的一个模糊子集，则隶属度函数为uN。()iiuN DP 为隶属度，其值可以有专家评判方法来确定。假设等级划分与指标间基本呈线性关系，则根据第 一 位 求出的 关 联度 来 确定隶 属 度 ，因 此 选 用0 .4 5 ,0. 2 9 ,0 .P，为了方便起见将隶属度扩大一百倍后可得：145P，224P，313P，48P，53P，60.02P，权数ig 的确定是人给定的， 选取的方法有很多， 本文选用 特费尔法， 并确定气氮、 EW、NS、雨量、电压地温的权数分配为：1,0.95,0.85,0.80,0.70,0.60g名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 20 页 - - - - - - - - - 9 这里g为126(,)EE EE上的模糊子集，按模糊数学的常规技法1234560. 9 50. 8 50. 8 00. 7 00. 6 01gEEEEEE（7）模糊关系矩阵U的确定根据126,PP PP ，126,gg gg就可以得到在集合D上所需的模糊关系矩阵：112111222212,.,nnmmnmg P g Pg Pg Pg Pg PUg Pg Pg P（8）由第一问可得各因素的敏感度为：电压EW NS 地温水位气温气压水温气氮雨量0.2 1.8 0.8 0.1 2595.3 2.3 4.3 40.2 2.8 0.5 表 4：地震前兆指标敏感度的表每次地震的等级：D=2.3 3.6 3.1 3.1 3.2 3.0 4.4; 稳定性分析数学模式：0.0761 5.0850 0.3290 0.0920 3.2349 0.1262 0.0589 1.3422 1.2721 0.0548 2.2143 0.5503 0.4320 0.9064 0.5601 0.1643 3.7511 0.9895 0.3576 0.8787 0.6075 0.1105 3.2558 0.2348 0.0456 1.8595 1.6741 0.1224 1.8006 0.2730 0.0715 0.8196 0.4019 0.0172 2.7385 0.8949 0.0719 0.6806 0.1246 0.0048 3.0955 0.0307 表 5：每次地震前兆稳定性系数地震指标权重-6.1024 电压- 0.5740 东西磁场强度1.0214 南北磁场强度2.0882 水位1.6259 气氡-0.8900 雨量假设地震发生等级与地震发生前数据关系为D=Ub+g; RI=7337;4948;1932;2259;3265;4473;7401; 在 7 次地震中，选中的每个指标在地震发生前的平均值：26.8611 3.7737 4.5667 17.7488 12.0109 0 26.9050 -2.0332 4.7336 17.5866 8.8819 0.0277 27.0236 -3.0234 3.3379 17.5680 12.0243 0 27.0298 -2.7623 3.6132 17.5712 12.4378 0 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 20 页 - - - - - - - - - 10 26.7648 0.1313 6.9721 17.6188 33.2737 0 26.7152 -1.3634 6.4789 17.6202 17.9615 0.6602 26.7021 1.8757 8.5844 17.6550 27.2101 0 表 6：每次地震前每个指标在地震发生前的平均值在这 7 次中上述每个指标的均值为：26.8574 -0.4859 5.4695 17.6241 17.6857 0.0983 在这 7 次中上述每个指标的方差为：0.1369 2.5458 1.9290 0.0631 9.1568 0.2480 由此得出第 2,9 个指标与地震的级数有相当大的关系，而其他的指标为稳定指标。即当电压达到26.8754，气氡达到 17.685左右时就有地震发生的可能性。而 1,3,4,10这几个指标与震级有关 ;假设是线性关系。即：D=a1*E(1)+a2*E(3)+a3*E(4)+a4*E(10), 有上述地震发生前的平均值可得系数为：a1 a2 a3 a4=3.1119 0.4654 -4.6999 -0.4878; 5.2.3 地震前兆指标数量特征模型先查阅与地震相关的其他指标如震级，频率等，而后先分析单个指标的数量特征，并以电磁波EW为例，利用动态从属函数建立模型，对于气温，水位等其他指标，类推可得。而后利用震级，频率等指标的数学研究方法，将已知条件的10 个指标综合起来，利用非线性组组合建立模型。电磁波 EW 的数量特征研究。地震活动是一种复杂现象，各种观测数据之间缺乏对比性，难以利用统一的物理场进行描述。 为了将各种观测量的异常信息进行表达，可把各前兆异常量转化为无量纲量以便进行对比。 各种地震前兆观测量虽然是不同的物理量，但其共同的特点都是随时间变化的数值，即都是时间函数， 各种地震前兆异常也都表现为各种观测量随时间的突出变化，异常形态虽多种多样， 但究其本质共同点都是观测曲线随时间的斜率变化。 因此，观测曲线随时间的斜率变化将是判断异常的重要指标。由以上分析从各种前兆数据中提取地震信息量的方法应具以下特点：（1） 、地震信息是一个无量纲量；（2） 、用前兆观测量Y（t ）对时间轴的斜率变化K，反映观测值的速率变化，一般情况下， K值越大，异常量越大；（3） 、用前兆观测值的自相关系数r 来显示观测值质量的好坏，r 值越大，观测数据质量越高。原理和方法从地震单项前兆观测值中识别地震异常，并计算它属于地震异常的程度， 需要遵循的基本条件是：（1） 、观测本身内在质量要好；（2） 、 “异常”出现的频次要低，且出现“异常”与“地震”之间对应程度要高；（3） 、 “异常”与“正常”相比，有明显的差异性，且界限要清楚。按照地震前兆异常的普适性和上述基本条件，提出如下地震前兆异常特征量：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 20 页 - - - - - - - - - 11 11( )1iiitk r（9）式中：22()jjjjjjjjjjjtyntyitntK（10）112222221()11()()jjjjjjijjjjjjjjt ytynrttyynn（11）jy 是前兆观测随时间变化的数列，jt 是相应的时间序列公式（ 1）中的 k 物理意义是观测曲线Y(t) 对于时间轴的斜率，反映观测速率的变化； r 是观测值的自相关系数，反应观测值内在质量的好坏。从（1）式可知，只有 k 和 r 都大时，才可能大。在0,1内取值，表示在所研究的资料在jt时段属于地震异常的程度。0.5i作为模糊界点，只有当0.5i时，才视为异常，反之，即为正常状态。是一个经过学习得到的常数。它的选择原则和方法如下：（1） 使得“异常”频次降低；（2） 使“异常”和“地震”之间由最好的相关性；（3） 在满足上述两个原则下要尽可能的小；（4） 选择的“地震”样本要视监视区的活动水平和预报目标而定；（5） 选择后的值，应保持稳定，只有出现预报目标与实际情况有明显差异时，才按照新的样本要求进行修改，以保持第二条原则的实现。模型的检验选择样本从属函数中的的选择由样本学习而决定，不同学科不同地域选择的参数可能会不同，需要在实践中不断的完善， 一些前兆场可直接根据原始观测数据分析其变化，而地磁场本身是一个变化场，且磁层、电离层强烈扰动时引起的地磁场变化非常剧烈， 远大于地震异常引起的变化。 从属函数法以曲线斜率作为判断异常的组要指标，如果直接分析某一点的从属函数值，会频繁出现“异常”，判断食物。同样是前兆手段， 地磁学在应用从属函数方法时首先对原始地磁场资料预处理，才由利于寻找震磁关系。电磁波由地磁场产生，而地磁场值包含核远源场，外源场及地下局部磁异常源场，提取电磁波EW的问题是如何找出因地磁场的地下局部磁异常发生变化的问题，并以此判断该变化量是否是由地震引起的局部异常。核源场是一个相对稳定场，外源场的变化受诸多因素影响， 而且其变化幅度远大于因地震引起的异常量，因此应用在应用资料预报之前，首先应剔除外源场引起的地磁场变化。任一名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 20 页 - - - - - - - - - 12 点的地磁场垂直强度Z 由以下三部分组成：CeNZZZZ（12）C-地核源场， e-外源场， N-地下局部磁场。由此可推得下面公式：CeNEWEWEWEW（13）下面采用均值法消除偶然因素如太阳黑子等造成的突变，并用对应一年中的均和几年的均值作差以消除外源场和核源场的变化而引起地磁场的突变NNZZZZZ（14）由此可推得：NNEWEWEWEWEW（15）正常情况下地点确定后局部磁场在很小的范围内波动；由于地震孕育而导致局部磁场发生变化时， 就在孕震区产生局部磁异常，这个范围比较小， 即使两地都处于孕震区，异常将由于距离震源距离不同而有所差别，最终导致EW的变化，EW代表局部磁异常的相对变化。从属函数的计算对应年份观测数据和各年份数据的平均作差消除基本场与外空场而引起的EW 变化, 用于消除外源场和核源场的变化而引起地磁场的突变，取各年份的均值是为偶然因素而引起的电磁波突变经过多年样本学习，确定值，研究电磁波随时间的变化情况。每次地震每个指标相对应的值电压EW NS 地温气氡雨量第一次地震0.5000 0.9867 0.9666 0.9994 0.9988 0.9916 第二次地震0.9883 0.5000 0.9669 0.5000 0.9997 0.8144 第三次地震0.9984 0.9613 0.9255 0.7123 0.9937 0.9236 第四次地震0.9786 0.9916 0.9882 0.9889 0.9999 0.8648 第五次地震0.9994 0.9837 0.5000 0.9923 1.0000 0.9836 第六次地震0.9927 0.9950 0.9943 0.9613 0.9999 1.0000 第七次地震0.9630 0.9989 0.9949 0.9982 0.5000 0.5000 表 7：每次地震前兆预测异常特征量值的分析：的取值范围为0,1，0.5 为判别地震发生的界点。当0.5时，则地震没有发生的可能。当0.5时，地震有发生的可能，值越大，则地震发生的可能性越大。每个指标的值名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 20 页 - - - - - - - - - 13 1.0e-003 * 0.0117 0.1812 0.3049 0.0021 0.0142 0.0736 图 5： EW 两年相对变化趋势图5.2.3 地震前兆指标数量特征模型针对地震预测研究中， 变量的指标较多， 数据携带的信息较弱且预测变量是非线性关系，因此造成预测结果稳定性不强的问题， 构造了非线性组合预测模型，解决地震预测的上述问题，应用于本题中的地震前兆指标数量特征。即使数据携带的信号很弱， 但是我认为这些数据足以通过统计方法来探求深层次的东西， 因为数据不服从标准的平稳时间序列和线性模型，所以我们要小心谨慎地选择统计方法。 下面构造的多变量非线性模型在一定程度上解决了上述问题且在地震预报应用中取得了较满意的效果。非线性组合模型非线性组合模型是将多个看用于预测的一元线性模型进行线性组合，组成非线性组合模型，其优点是可以增加预测稳定性。非线性组合模型的应用背景设y是预测变量，12,px xx 是影响y的指标。（1） 、 指标数据携带信号较少、 较弱，需要选择携带信号较强的指标用于预测。（2） 、对于选定的指标需要建立较好的非线性预测模型。（3） 、为增强预测的稳定性，最后构建非线性组合模型用于预测。最佳指标和最佳曲线选择对于每一个指标ix(1,2,)ip，对应预测变量y，预先设计若干非线性模型，并进行曲线拟合，根据统计检验和拟合优度选择最优指标和最佳曲线。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 20 页 - - - - - - - - - 14 设所选的最优指标为ix 相应的最佳曲线函数是()(1,2, )if xik ，k是最优指标和最佳拟合曲线个数。非线性组合模型1()kiiiiyw fx（16）其中y是预测变量，()iif x是选择出的最佳非线性预测曲线，ix 是自变量指标，iw 是相应的权重，11kiiw。非线性组合模型参数iw 的估计一般参数估计的方法可选择：（1）等权平均法；（2）加权最小二乘法；（3）方差 -协方差方法；（4）最小二乘法；非线性组合模型的应用在本题中，选择最优指标和最佳曲线， 通过最小二乘法估计最佳曲线的参数，并构建非线性组合模型。最佳指标和最佳曲线的选择初步选择影响最大的若干地震参数指标ix 等，见表一为选择最优指标和最佳拟合曲线设计了七种曲线模型，见表二。选择拟合方式以及统计检验， 选出 6 个最优指标以及他们相应的最佳曲线。见表三。2.2 、非线性组合模型6111()6kiiiiiiyw fxy考虑到增加预测的稳定性，参数估计采取等权平均法16iw，(1,2,6)i经过 Matlab 编程得出了地震各因素的权重如下：电压EW NS 地温气氡雨量-0.0000 -0.0029 0.0001 0.0057 -0.0014 -2.3620 画出部分地震前兆指标与地震发生的关系图如下：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 20 页 - - - - - - - - - 15 图 6：相对震级图5.2.4 地震数据分析平台的设计把前三问的程序集成，基于易语言编程，系统开发使用Windows2000server、MS-SQLServer2000数据库，对地震目录、前兆观测数据、分析报告等地震信息进行了系统管理， 组成了地震数据分析平台， 本平台在实现了地震数据的数字化和自动化管理功能的基础上，系统运用数据挖掘的方法对地震预报进行了尝试，使其能够完成形如题给数据的其它地震数据的分析，并能自动输出前述任务的重要的分析结果。其主界面如下图. 图 7：地震数据分析平台主界面基于此地震数据分析平台，分析了2009 的地震指标数据，得出的结果是地名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 20 页 - - - - - - - - - 16 震发生时间为 10 月 31 日，震级为 4.1 级，在误差范围内， 实现了较好的地震前兆数据分析与得出重要的分析结果。六、模型的推广1、在地震分析中经常要分析各个变量间相互依存关系，本文模型只能分析每一个变量对因变量的直接影响，且要求其余自变量之间相互独立。事实上，变量之间的关系错综复杂，一个自变量可能对一个因变量有直接影响，也可能通过其他的自变量对因变量有间接地影响，显然本模型无能为力。2、可选用通径分析来求解各个变量之间的关系，在通径分析模型中，把影响其他变量的变量叫做原因变量，被影响的变量叫做结果变量。七、关于进一步研究报告地震预报预测之艰难不言而喻， 然而通过此次研究分析发现， 地震发生前并非不显示任何蛛丝马迹， 笔者基于促进其预测预报研究的目的提出以下建议和设想：确保各指标观测数据的真实性、准确性。数据是供人们研究的平台， 其真实性、准确性直接关乎研究方向的正确性及研究结果的实效性。优化数据挖掘方法。从各指标对事件的函数图像可知，各指标值振荡频繁，然而不是所有的振荡都是因为地震的发生而引起的，如电磁波的振荡容易受太阳黑子等因素影响而发生明显振荡。如何采取科学、 合适的数据挖掘方法关系到所得模型预测预报的准确率， 我们所采取的是均值结构模型以达到消除偶然误差的目的，其效果显著，然而由于原理上的限制，还是有所偏差，所以笔者设想能否利用数理统计方法进行估计，以使偏差尽量降低。研究地震发生时各指标异常机理。研究过程发现， 各指标异常会出现一定的提前量或滞后性， 如何估计这些提前量关乎预测预报的准确性，只有通过研究各指标因地震变化的机理才能用数学方法较为准确的估计这些提前量。优化数学研究方法。 本模型在求解任务一时采用均值结构模型达到较良的消除误差的效果；采用综合指标法求解任务二，考虑全面周到，效果较好，在模型改进中所提出的判别分析法， 科学有效， 对数据利用率较高。 在任务三中采用各年份指标值和历年各指标均值作差以消除地核、外空对电磁波的影响， 巧妙科学，适宜推广使用。 任务三中采用非线性组合模型， 较好的解决数据携带信号较弱的缺点。限于知识、 能力，对于上述数学方法在地震预测预报中的应用没有达到预期效果， 须进一步研究改进。由于地震的发生及受其影响的各指标间的复杂关系，很难采用现有的技术探测，笔者试想，能否应用Bp 神经网络解决此难题，设定科学的学习规则，利用神经网络的自我学习特性，达到预测预报的目的。参考文献名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 20 页 - - - - - - - - - 17 1 徐文耀，我们地震观测研究发展，地球物理学报，997,40:217 。2 冯德益，等模糊数学方法应用，北京，地震出版社。3 郭德科， 强震前兆异常的模糊识别函数与中弱地震的关系，地壳形变与地震。4 张岩波，潜变量分析，北京，高等教育出版社。5 王清河，随机数据处理方法（第三版） ，石油大学出版社。6 张继红、乔慧珍，地震前兆异常的动态从属函数，地球物理学进展。7 沈良峰，地震稳定性分析的综合性指标模型，建筑科学。附录1：日平均值的程序% 每年数据的平均值RI=7344 4968 1944 2280 3288 4488 7416; year=2005 2006 2007 2007 2008 2008 2009; S = load(matlab.mat); er=zeros(40*7,6); for i=1:7 u=strcat(S.t,num2str(year(i),(,num2str(RI(i)-40*24+1),:,num2str(RI(i),1 2 3 4 9 10); data=eval(u); temp1=zeros(40,6); for j=1:40 temp1(j,:)=mean(data(24*j-23):(24*j),:); end er(40*i-39):(40*i),:)=temp1; end K1=zeros(7,6);% 每次地震每个指标的K值r1=zeros(7,6);%每次地震每个指标的r 值for i=1:7 K1(i,:)=sum(1:40)*sum(er(40*i-39):(40*i),:)-40*sum(diag(1:40)*er(40*i-39):(40*i),:); K1(i,:)=K1(i,:)/(sum(1:40)2-40*sum(1:40).2); r1(i,:)=sum(diag(1:40)*er(40*i-39):(40*i),:)-(1/40)*sum(1:40)*sum(er(40*i-39):(40*i),:); r1(i,:)=r1(i,:)./(sqrt(sum(1:40).2)-(1/40)*sum(1:40)2)*sqrt(sum(名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 20 页 - - - - - - - - - 18 er(40*i-39):(40*i),:).2)-(1/40)*sum(er(40*i-39):(40*i),:).2); end gh=zeros(7,6); for j=1:6 gh(:,j)=K1(:,j).*r1(:,j); end miu=zeros(7,6);% 每次地震每个指标的miu 值for j=1:7 miu(j,:)=1./abs(1+min(gh)./gh(j,:); end 2 年平均曲线的 Matlab 程序function s=years_avg() yinsu=data VI EW NS ET WH AT AP WTAID RAIN PNS PEW; year=2005 2006 2007 2008 2009 2010; qs=input(请输入你所需要的第几年的编号：); qw=input(请输入您要的因素编号： ); coona=database(sql server,sa,xuenhappy); setdbprefs(DataReturnFormat,numeric); s=zeros(1,372); for i=1:12 for j=1:31 if i10&j10 t=strcat(yearqs,-0,num2str(i),-0,num2str(j),%); t=strcat(char(39), t,char(39); t= ,t; else if i9 t=strcat(yearqs,-0,num2str(i),-,num2str(j),%); t=strcat(char(39), t,char(39); t= ,t; else if i9&j10 t=strcat(yearqs,-,num2str(i),-0,num2str(j),%); t=strcat(char(39), t,char(39); t= ,t; else t=strcat(yearqs,-,num2str(i),-,num2str(j),%); t=strcat(char(39), t,char(39); t= ,t; end end end 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 20 页 - - - - - - - - - 19 s(j+31*(i-1)=fetch(coona, strcat(use earthquck select avg(,yinsuqw,) from years,yearqs,$ where Convert(Varchar(10),data,120) like,t) ); if isempty(s(j+31*(i-1) s(j+31*(i-1)=0; end end end close(coona); plot(1:max(size(s),s) 3 预测分析图的 Matlab 程序RI=7337;4948;1932;2259;3265;4473;7401; year=2005 2006 2007 2007 2008 2008 200