1問題的提出.doc

地震的分类和预测 中南大学 孙跃、吴俊儒、徐苗苗 摘要 本文旨在研究地震的分类和预测问题。根据国家地震科学共享中心提供的世界1970—2011年的数据按照多种指标进行了分类。在此基础上通过数学模型和统计方法相结合的方式对地震的预测问题进行了重点研究。分地震带对地震活动总体趋势、7级以上强震活动的趋势建立了数学模型，并对预测结果进行了分析，依据分析结果对模型进行了评价。 首先依照地震本身特征对地震进行了多种分类。根据地震系列(主震、余震、前震等)和震源所处地理位置的不同，都采取了动态聚类方法－K均值法对其进行了分类，并得到了我国基本的大地震带；根据震源深度、震级大小的不同，也对原始数据进行简单的分类。 然后用探索性数据分析(EDA)初步研究了地震活动的规律和趋势，对原始的海量数据进行了数据可视化处理。引入地震活动度的概念后，在R2.7.1软件上通过编程绘出全中国的地震度时序图(图3)和四个大地震区的时序图(图4)，对地震规律和趋势进行了定性的分析。用Tukey Smoothing方法对上面不平稳的活动度时序图进行了平滑，得到了四个地震带大致的活动度的趋势图(图5)，得出了一些粗略的结论。 考虑到地震活动度时间序列数据的自相关性，提出了ARIMA(p,d,q)模型对地震活动度进行了20个月的趋势预测。再用R2.7.1软件求解，得到中国2009年9月至2011年4月的地震活动度的预测值(表2)。且预测值通过了白噪声检验。并根据预测值画出了四个大地震带20个月的活动度趋势图(图7)。 考虑到强震对人类造成的巨大危害，基于地震活动系统的灰色性，建立了灰色预测模型预测7级以上的地震趋势。通过MATLAB7.0编程求解，预测了1990年——2011年及未来20年内，西部、华北、台海、南北四大地震带7级以上大地震发生时间(表6)。 最后本文根据地震数据不稳定的特点，创新性地提出了基于均生函数的时间序列外推模型，对1990年——2011年7级以上的地震数据进行了拟合优度检验，并将预测结果与灰色预测模型的结果进行了对比，得出了该预测模型优于灰色预测模型的结论。 在模型的求解之后，对每个模型都进行了评价。3个预测模型都很好的将影响地震的因素模糊化，通过地震系统的内在规律对未来进行预测。 关键字：Tukey Smoothing ARIMA 模型灰色预测 均生函数 时间序列外推 地震活动度 数据可视化 目录 1 问题的提出 1 1.1基本情况 1 1.2相关信息 1 1.3需要解决的问题 1 2 模型假设与符号说明 1 2.1 模型假设 1 2.2 符号说明 2 3 问题分析 2 3.1 问题一： 2 3.2 问题二 2 4 地震的分类模型 3 4.1 按地震系列划分 3 4.2 按震源位置划分 3 4.3 其它分类方法 5 5 地震的预测模型 5 5.1探索性数据分析 5 5.2 ARIMA（p,d,q）模型 8 5.2.1 模型的提出 8 5.2.2 模型原理 8 5.2.3 建模步骤 8 5.2.4 模型求解 9 5.2.5 模型的诊断和预测 10 5.2.6 模型评价 11 5.3 灰色预测模型 11 5.3.1 模型的提出 11 5.3.2 模型原理 11 5.3.3 模型求解 13 5.3.4 预测值精度的检验 15 5.3.5 灰色预测模型的评价 15 5.4 基于均生函数的时间序列外推模型 15 5.4.1 模型的建立 15 5.4.2 模型求解 17 5.4.3 模型评价 17 6 参考文献 18 7 附录 18 附录1：地图以及地震发生散点汇集图 18 附录2 地震带的划分 19 附录3 活动度及相关程序 20 附录4 ARIMA程序 22 附录5 灰色预测程序 23 附录6 地震系列聚类结果 24 1 问题的提出 1.1基本情况 地震是地下岩层受应力作用错动破裂造成的地面震动，是一种破坏性极强的自然灾害，是自然灾害之首恶。地震可以在很短的时间内使一座城市夷为平地，使无数的家庭支离破碎。为了有效的减小地震带来的损失，对地震进行预测则很重要。我们从以往已经发生过的地震记录中寻求规律，以达到对未来时间地震发生情况的预测。全球强震带主要集中分布在环太平洋地震带和欧亚地震带，而且身处其中的不同的地区间地震发生情况是互相关联的。我国亦处在其中，因此我国的地震活动和全球的地震活动是密切相关的，尤其是和板块边界的强震活动更是如此。 1.2相关信息 国家地震数据科学共享中心的中国地震台网数据目录中给出的全球从1970年1月1日到2011年4月25日的地震数据。 1.3需要解决的问题 1)根据文中给出的参考数据，以及地震的性质，对地震划分种类进行划分。 2)根据文中给出的参考数据，以