1、 题目：地震预测数学建模姓名：张志鹏 学号：12 291233 学院：电气工程学院姓名： 赵鑫 学号：10291033 学院：电气工程学院数学建模竞赛论文姓名：张书铭 学号：12291232 学院：电气工程学院目录摘要 3一、问题重述 4二、问题的分析 4三、建模过程 5问题 1：地震时间预测 51、问题假设 .52、参数定义 .63、求解 .6问题 2：地震地点预测 71、 .问题假设：72、 .参数定义83、 .求解过程：8四、模型的评价与改进 12参考文献 13摘要大地振动是地震最直观、最普遍的表现。在海底或滨海地区发生的强烈地震，能引起巨大的波浪，称为海啸。在大陆地区发生的强烈地震，会
2、引发滑坡、崩塌、地裂缝等次生灾害。对人们的生产生活成巨大影响，严重威胁人们的生命和财产安全，所以，对地震的预测是十分必要的。本文根据从 1900 年以来中国发生的八级以上地震的时间和地点分析，利用合理的数学建模方法，对下一次中国可能发生的八级以上地震的和时间和地点进行合理的预测。建模方法分为对于时间的预测和地点的预测两个方面。问题 1：对于时间的预测采用的方法为指数平滑法，它是通过计算指数平滑值，配合一定的时间序列预测模型对现象的未来进行预测。其原理是任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。问题 2：对于地点的预测根据长久的数据表明，八级以上地震主要发生在东经 701
3、10，北纬 2050这个范围内，据此将整个地震带划分为100 个区域，按顺序进行编号。建立时间与地震区域编号的数学模型，利用线性回归的方法对下次地震地点预测。关键词：地震，预测，数学建模，指数平滑法，线性回归一、问题重述地震预报问题，大地震的破坏性是众所周知的，为了减少大地震带来的灾难，人们提出了各种预报地震的方法，以求减少大地震产生的破坏。本赛题请大家用数学建模的方式预报下一次大地震发生的时间和地点。为了减少问题的复杂性，请根据 20 世纪以来我国发生 8 级及以上大地震的时间和地点预报下一次我国发生 8 级及以上大地震的时间和地点。二、问题的分析由题意可知，目的就是为了建立一种模型，预测下
4、一次八级医生地震发生的时间、地点。问题一中，首先建立两次地震相隔时间是按照一定的时序发生的，将此时序抽象为“1”，对下一次地震的时间间隔进行预测。在这个问题中，顺理成章的就会想到利用指数平滑法，指数平滑法是生产预测中常用的一种方法。也用于中短期经济发展趋势预测，所有预测方法中，指数平滑是用得最多的一种。简单的全期平均法是对时间数列的过去数据一个不漏地全部加以同等利用；移动平均法则不考虑较远期的数据，并在加权移动平均法中给予近期资料更大的权重；而指数平滑法则兼容了全期平均和移动平均所长，不舍弃过去的数据，但是仅给予逐渐减弱的影响程度，即随着数据的远离，赋予逐渐收敛为零的权数。也就是说指数平滑法是
5、在移动平均法基础上发展起来的一种时间序列分析预测法，它是通过计算指数平滑值，配合一定的时间序列预测模型对现象的未来进行预测。而问题二，是利用模块化分析以及线性回归分析结合对下一次地震发生的地点的预测。根据长久的数据表明，八级以上地震主要发生在东经 70110，北纬 2050这个范围内，利用模块化将东经 70110，北纬 2050这个地震范围划分为 100 个区域，按顺序进行编号。建立时间与地震区域编号的数学模型，利用线性回归的方法对下次地震地点预测。最后，根据网上提供的知识，再结合自己的亲身体验，提出了预测地震时间、地点的可行性方案。三、建模过程问题 1：地震时间预测1、问题假设（1） 首先根
6、据中国地震网的数据，算出两次地震的时间差时间 地点 震级 经纬度 相邻两次地震的时间 差（年）1902.8 新疆阿图什 8.3 E76.12 N39.731920.12 宁夏海原 8.5 E105.64 N36.56 181927.5 甘肃古浪 8 E102.86 N37.43 71931.8 新疆富蕴 8 E89.44 N47.05 41950.8 西藏察隅 8.5 E97.49 N28.62 191951.11 西藏当雄 8 E91.05 N30.51 12001.11 昆仑山口西 8.1 E90.90 N36.20 502008.5 四川汶川 8 E103.61 N31.46 7（2）
7、相邻时间差只是按照时间序列排列的，与其他因素无关。（3） 利用指数平滑法对下一次地震相邻时间差进行预测2、参数定义式中， -t+1 期的预测值，即本期（t 期）+1=+(1) +1的平滑值 ； -t 期的实际值； -t 期的预测值，即上期的S 平滑值 。S13、求解利用 excel 中的指数平滑数据分析方法得到下表 +118 =0.2 =0.4 =0.6 =0.87 18 18 18 184 9.2 11.4 13.6 15.819 5.04 6.96 9.76 13.441 16.208 14.184 13.456 14.55250 4.0416 6.2736 8.4736 11.84167
8、 40.80832 32.50944 25.08416 19.473281 2 3 4 5 6 70204060实 际 值预 测 值指 数 平 滑=0.2数 据 点值1 2 3 4 5 6 70204060实 际 值预 测 值指 数 平 滑=0.4数 据 点值1 2 3 4 5 6 70204060实 际 值预 测 值指 数 平 滑=0.6数 据 点值1 2 3 4 5 6 70204060实 际 值预 测 值指 数 平 滑=0.8数 据 点值根据分析与求解，在 =0.8 的时候预测值与实际值的拟合度最好，所以根据公式 +1=+(1)下一次地震时间差为：=0.8x7+(1-0.8)x19.47
9、3=9.49+1据此可知下次地震时间差为 9.49，因此下一次地震事件在20172018 年之间。问题 2：地震地点预测1、 问题假设（1） 、根据长久的数据表明，八级以上地震主要发生在东经 70110，北纬 2050这个范围内，据此将整个地震带划分为 100 个区域，下次地震一定发生在这个范围内。（2） 、地震发生的区域只与时间有关系（3） 、按顺序进行编号。建立时间与地震区域编号的数学模型，利用线性回归的方法对下次地震地点预测。2、 参数定义y 为地震区域，x 为地震发生时间3、 求解过程（1） 、根据假设将八级以上地震范围划分为 100 个区域，并按照顺序编号。（2） 、然后将八次地震的
10、区域记为 y，发生时间为 x，建立数学模型（3） 、利用线性回归的方法求得 x 与 y 的拟合方程（4） 、有问题 1 中求出的 x 算出 y，并根据 y 的值确定所在区域的经纬度，即为下次地震所发生的区域。5047.54542.54037.53532.53090 91 92 93 94 95 96 97 98 9980 81 82 83 84 85 86 87 88 8970 71 72 73 74 75 76 77 78 7960 61 62 63 64 65 66 67 68 6950 51 52 53 54 55 56 57 58 5940 41 42 43 44 45 46 47 4
11、8 4930 31 32 33 34 35 36 37 38 3920 21 22 23 24 25 26 27 28 2910 11 12 13 14 15 16 17 18 190 1 2 3 4 5 6 7 8 927.52570 74 78 82 86 90 94 98 102 106 110【经纬度分区图】根据上表列 x 与 y 关系，并用 excel 求线性回归方程：x y1902 511920 481927 481931 841950 161951 252001 452008 281900 1920 1940 1960 1980 2000 2020020406080100Seri
12、es 1预 测 51Linear (预 测 51)拟 合 曲 线190251SUMMARY OUTPUT回 归 统 计Multiple R0.376165R Square 0.1415Adjusted R Square-0.0302标 准 误 差 22.7185观 测 值 7方 差 分 析 df SS MS F Significance F回 归 分 析 1425.3491425.34910.8241120.405608残 差 52580.651516.1302总 计 6 3006Coefficients标 准 误 差 t Stat P-valueLower 95%Upper 95%下 限 9
13、5.0%Intercept506.1405511.34920.9898140.367716-808.3241820.605-808.3241902-0.237360.261465-0.907810.405608-0.909480.434758-0.90948RESIDUAL OUTPUT观 测 值 预 测 51 残 差150.40932-2.40932248.74781-0.74781347.7983736.20163443.28853-27.2885由上表，去除最大的坏值点即（1931，84 ）之后重新求拟合直线。x y1902 511920 481927 481950 161951 252
14、001 452008 281900 1920 1940 1960 1980 2000 202001020304050Series 1预 测 51Linear (预 测 51)修 正 后 的 拟 合 曲 线190251无明显坏值！由分析计算可知：y = -0.0868x + 205.12则下一次地震区域在 y= -0.0868x2017 + 205.12=30.044 区域 31，对应实际地理位置即：东经 7478，北纬 32.535的范围内。四、模型的评价与改进模型最大优点在于对原始数据拟合时, 采用多种方法进行, 使之愈来愈完善, 具有很高的拟合精度和适度性在此基础上, 对模型作进一步讨论便可得到一系列可靠而实用的信息并且, 所得结论与客观事实很好地吻合, 从而进一步说明模型是合理的。不足在于忽略了很多其他影响因素，如地壳变化等因素，使模型在合理性与可靠性上可能存在一定的缺陷。本文建模，我小组方法对于前段时间尼泊尔大地震具有一定的预见性。根据我小组的模型计算，与尼泊尔地震位置时间相距很近，因此我小组认为我们的地震预测数学模型具有一定的可行性！参考文献【1】沙海军、陈虹，时间 震级可预测模型在中国大陆的应用，地壳构造与地壳应力文集（15） ，2002【2】张忠平，指数平滑法，北京统计出版社，1996【3】王松桂，线性统计模型，高等教育出版社，1999
