数学模型的-应用实例课件

1、数据驱动的地震余震序列分析模型,数据挖掘方法建模（ Data-Mining,数学建模工作室 2021/1/26,第1页 / 共19页 mecca_,数学模型的-应用实例,完成于2008年5月23日,科学研究的交叉思维方式,数学建模工作室 2021/1/26,第2页 / 共19页 mecca_,数学模型的-应用实例,解决问题的途径,数学建模工作室 2021/1/26,第3页 / 共19页 mecca_,现实问题,内在机理,外在表现,模型驱动（Model-Drived,数据驱动 （Data-Drived,地震余震发生趋势预测,Han J, Kamber M. Data Mining: Concep

2、ts and TechniquesM. Second Edition ed. Singapore: Elsevier, 2006,万有引力定律行星椭圆型运行轨迹,Brahe 的观测数据行星椭圆型运行轨迹,数学模型的-应用实例,1 问题重述（研究意义,2008年5月12日，四川汶川发生了8.0级的特大地震，给中国带来了特大灾难和无法估量的生命和财产损失。同时，从5月12日开始到22日为止，已经发生4级以上的余震达170余次，并继续不断的发生余震。能否快速、有效的根据现有已经发生的余震数据，分析和判断余震的发展趋势对安定人心和指导救灾工作的具有极其重要的意义。本文就此问题进行数学建模，并给出对余震

3、发展趋势的分析和判断,数学建模工作室 2021/1/26,第4页 / 共19页 mecca_,Introduction,数学模型的-应用实例,数学建模工作室 2021/1/26,第5页 / 共19页 mecca_,数学模型的-应用实例,数学建模工作室 2021/1/26,第6页 / 共19页 mecca_,数学模型的-应用实例,由于地震的物理机理尚不明确，大多数现有地震理论都具有一定的适用范围，对某地区发生的地震有效并不能保证对其它地区地震的有效1-11。我国几次准确的地震预报案例很大程度上是依赖“群测群防”这样密集“探测”的经验预测12。对地震余震序列的研究，还停留在“统计分析”这个水平，并

4、没有完全理解其真实物理含义1-4, 7, 13。同时，一个地区的地震发生频率具有一定的“混沌”特征14，使得对余震序列的分析和研究十分困难,数学建模工作室 2021/1/26,第7页 / 共19页 mecca_,2 问题分析（研究方法：前人,Introduction,数学模型的-应用实例,在这种背景下，将地震余震序列视为“黑箱”的基于数据自身特性的分析方法具有更好的效果。考虑到余震的发生具有很强的随机性，准确预测其发生时间尚不具备条件，故退而预测其发展趋势。借用灰色系统理论15的数据“累加”处理方式，在消除其随机性的基础上，以数据本身的特性驱动建立数学模型，是本模型的最大特色,数学建模工作室

5、2021/1/26,第8页 / 共19页 mecca_,2 问题分析（研究方法：本文,Introduction,数学模型的-应用实例,数学建模工作室 2021/1/26,第9页 / 共19页 mecca_,3 模型假设（研究基础,1. 假设地震的发生具有自身特定的规律，并具有一定的稳定性,2.假设主-余型地震的所有余震的发生与主震特性相关，并具有统一的特性,3.所有地震数据具有一定的精度和可靠性,Introduction,确保数学模型为“常”系数,确保反演计算时基础可靠,数学模型的-应用实例,数学建模工作室 2021/1/26,第10页 / 共19页 mecca_,4 模型建立,Method,

6、已知余震序列数据 ，其中 表示在 时刻发生的余震震级,数学模型的-应用实例,数学建模工作室 2021/1/26,第11页 / 共19页 mecca_,4 模型建立,Method,处理后的数据已经具有明显的规律性,旧方法，新问题,邓聚龙. 灰色控制系统（修订版）M. 华中理工大学出版社, 1993,数学模型的-应用实例,数学建模工作室 2021/1/26,第12页 / 共19页 mecca_,4 模型建立,Method,其中，(a,b,c)为模型参数，可通过具体数据反演计算出来,根据处理后的数据所反映出的特性和模型假设，我们建立如下常微分方程模型描述累加后数据曲线,1,姜启源. 数学模型（第二版

7、）M.高等教育出版社，1993 （传染病模型、人口Logistic模型,旧方法，新问题,假设1，2 （a,b,c是常数,数学模型的-应用实例,数学建模工作室 2021/1/26,第13页 / 共19页 mecca_,5 模型求解,Method,代入（1）式有,2,关于(a,b,c)的线性方程组,假设3原始数据可靠,数学模型的-应用实例,数学建模工作室 2021/1/26,第14页 / 共19页 mecca_,5 模型求解,Method,根据最小二乘原理，我们可反演计算出方程（1）中的参数,设,则方程（2）可写为,可利用Matlab软件中的ode45()常微分方程求解器来计算方程的数值解,数学模

8、型的-应用实例,数学建模工作室 2021/1/26,第15页 / 共19页 mecca_,6 计算结果及分析,Results,实际计算时， 是由Internet资源得到汶川地震余震序列,数学模型的-应用实例,数学建模工作室 2021/1/26,第16页 / 共19页 mecca_,6 计算结果及分析,Results,建立模型数据,检验模型数据,模型预测数据,数学模型的-应用实例,上图给出本数学模型的计算结果，其中蓝色加号表示参与反演计算的原地震余震序列的发生时间和震级，绿色叉号表示未参与反演计算的原地震余震序列的发生时间和震级；蓝色叉号表示累加后的参与反演计算的余震序列的发生时间和震级，红色圆

9、圈表示未参与反演计算余震序列的发生时间和震级；红色虚线表示模型预测的 地震震级累加变化率曲线，蓝色虚线表示实际地震余震 序列震级累加变化率曲线,显然，本数学模型的预测结果与实际发生的地震余震序列（红色圆圈）十分吻合，验证了模型的有效性。根据数学模型的计算结果，汶川地震余震还将持续很长一段时间，但总趋势是明显的在有规律减弱,数学建模工作室 2021/1/26,第17页 / 共19页 mecca_,6 计算结果及分析,Discussion,数学模型的-应用实例,最后，由于建立预测模型对时间轴也进行了累加处理，本文模型无法预测余震的发生时间和震级，只能预测余震的发展趋势，这是将来需要改进的研究方向。

10、 （缺点,本文以数据自身特性为驱动力建立了地震发生余震的发展趋势预测模型，并以实际数据验证了模型的有效性，能对余震发展趋势的分析和预测提供有效的数据分析工具。所建模型理论基础非常简单，可使大多数具有高等数学知识的工程技术人员快速分析地震余震序列数据，有效避免谣言的传播和产生。（优点,数学建模工作室 2021/1/26,第18页 / 共 19页 mecca_,7 结论及模型改进,Conclusion,同时，如果后续发生的地震数据严重偏离模型预测曲线，可迅速判定余震发生的性质改变趋势：出现新的余震诱发因素（数据向上偏离曲线）或余震能量由于特殊原因停止释放，可能诱发更大的地震（数据向下偏离曲线）。

11、（优点,数学模型的-应用实例,参考文献,1 蒋海昆, 郑建常, 吴琼, et al. 中国大陆中强以上地震余震分布尺度的统计特征J. 地震学报. 2007, 29(2): 151164. 2 蒋海昆, 郑建常, 代磊, et al. 中国大陆余震序列类型的综合判定J. 地震. 2007, 27(1): 1725. 3 毛燕. 姚安地震余震序列研究J. 地震研究. 2007, 30(2): 142145. 4 吕晓健, 高孟潭, 高战武, et al. 强余震和主震地面运动分布比较研究J. 地震学报. 2007, 29(3): 295302. 5 李琼. 复杂储层地震预测理论及方法研究D. 博士

12、, 成都理工大学, 2007. 6 余怀忠. 地震前兆物理模型与地震预测初步研究R. Technical 中国 北京: 中国地震局地质研究, 2004. 7 刘正荣. 极限时间法与临震预报J. 地震. 2005, 25(1): 4150. 8 李树华, 陈燕. 云南地区地磁加卸载响应比短临异常特征J. 地震研究. 2003, 26(增刊): 133139. 9 钱晓东, 秦嘉政, 叶建庆. 云南地区的环境应力参数及其在地震短临预报中的应用J. 地震研究. 2006, 29(1): 713. 10 焦青, 刘耀炜, 杨选辉, et al. 文安5 . 1级地震前兆变化特征分析J. 大地测量与地球动力学. 2007, 27(6): 9699. 11 夏雅琴, 陈维升, 李均之, et al. 地震预测的新进展J. 北京工业大学学报. 2006, 32(1): 1114. 12 刘成龙, 车用太, 王广才. 大规模宏观异常的双重性及其在地震预报中的意义J. 地震地质. 2004, 26(2): 340346. 13 蒋海昆, 郑建常, 吴琼, et al. 传染型余震序列模型震后早期参数特征及其地震学意义J. 地球物理学报. 2007, 50(6):