2022年地震模型 .pdf
《2022年地震模型 .pdf》由会员分享,可在线阅读,更多相关《2022年地震模型 .pdf(20页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、1 地震预测模型摘要地震预报是当代自然科学领域里一个直接关系到人民生命安全和社会经济的发展,是科学难度很大的前沿课题。近几十年来, 人民的努力虽然取得了丰硕的成果,预报的实践有些有限的成功。在地震发生前有很多前兆性指标,如磁场强度、 氡值、电压、雨量、水温等,这些指标都从不同的侧面反映了地震活动的各项特征。在正常情况下, 也常常有些指标有明显的异常, 而令一些指标并不出现异常。 这些都给实际的前兆指标数据特征分析以及地震预报工作带来了困难。本文针对地震影响因素多, 数据多, 联系弱的特点, 分别建立了针对各个任务的数学模型, 首先,处理初始数据时选取日平均变化来消除一天中数据的随机因素的影响,
2、然后更具这些数据建立了地震发生前后各个指标的不同程度度量模型解决了地震对指标的敏感度的度量问题,进而找出了EW ,气氡,水位,电压,雨量这些衡量地震的主要因素, 并分析出了这些指标的重要程度,在解决找统计量时,利用上次任务所得的主要因素的组合来衡量一次地震的地震前兆指标数量特征,并得出了电压达到26.8754,气氡达到 17.685 左右时就有地震发生的可能性,考虑全面周到, 效果较好,在模型改进中所提出的判别分析法,科学有效,对数据利用率较高。 而除此之外的其他几个指标与震级有关,这样就这些指标的分析在任务三中分别建立了地磁前兆异常的动态从属函数模型以及非线性综合模型,进而找出了地震级数与这
3、些因素之间的非线性关系,对于上述数学方法在地震预测预报中的应用没有达到预期效果,须进一步研究改进。综合这些工作,我们应用易语言编写了分析地震数据的平台,通过分析地震数据生成了报告。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页,共 20 页 - - - - - - - - - 2 一、 问题重述地震是地壳快速释放能量过程中造成的振动。虽然预测地震是世界性难题,但迄今科学界普遍认为,有可能反映地震前兆特征的指标可能不少于10 个。已经有专业仪器在多个定点实时按秒记录这些指标的数
4、据,期望通过对记录数据的分析研究找到地震的前兆特征。现已采集到某地 2005年 1 月 1 日至 2010年 6 月 30 日按小时观测的 10 多个指标的数据, 和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。 这些数据中隐藏着地震发生的前兆特征。科学地截取这些数据的有用片段,对数据进行合理地预处理, 用数学方法揭示地震前兆的数据特征,是一项很有意义的研究工作。题给数据中的这 10 多个指标,究竟哪些与地震的发生有关,有何种关系,是单一关系还是复合关系;除这10 多个指标外还有哪些因素及含题给指标在内的哪些指标的哪种数学模型更能反映地震的前兆特征等等,人们迄今仍不很清楚,需要
5、进行深入地研究。 地震数据的观测是持续进行的,随着时间的推移数据的规模会不断扩大。 从中挖掘地震的前兆特征, 必须有合理的数学模型, 也必须有科学高效的算法分析平台。因此,请参赛队结合附件中给出的实际记录数据,尝试完成以下任务。任务一:分析数据特征, 建立数学模型以度量各指标对地震发生的敏感程度。任务二:构造由某些或全部指标构成的综合指标,使其尽可能地集中反映地震发生前的数据特征的统计规律。任务三:结合题给数据, 广泛查阅与地震相关的其它指标的数据和分析方法,建立数学模型来研究地震发生前的数量特征。任务四:将前述各项任务的计算程序集结成地震数据分析平台,使其能够完成形如题给数据的其它地震数据的
6、分析,并能自动输出前述任务的重要的分析结果。任务五:对于进一步的研究设想写一篇切实可行的报告。二、问题分析地震是地下岩层受应力作用错动破裂造成的地面震动,是一种破坏性极强的自然灾害,是自然灾害之首恶。地震可以在很短的时间内使一座城市夷为平地,使无数的家庭支离破碎。 为了有效地减小地震带来的损失,对地震进行预测则很重要。本题中给出了地震发生前2005 年 1 月 1 日至 2010 年 6月 30 日按小时观测的 10 多个指标的数据和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。 这些数据中隐藏着地震发生的前兆特征。在地震发生前有很多前兆性指标,如磁场强度、氡值、电压、雨量、水
7、温等,这些指标都从不同的侧面反映了地震活动的各项特征。 在正常情况下, 也常常有些指标有明显的异常,而令一些指标并不出现异常。 这些都给实际的前兆指标数据特征分析以及地震预报工作带来了困难。21、问题的性质本文主要是对所给数据进行处理,用数学方法揭示地震前兆的数据特征。2.2 、解决问题的难点和关键名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页,共 20 页 - - - - - - - - - 3 1、如何确定前兆个指标对地震发生的敏感程度。2、如何构造地震发生的前兆综合指标
8、,尽可能的揭示地震发生前的数据特征的统计规律。3、构建地震数据分析平台。三、模型假设1、假设地震发生前2005年 1 月 1 日至 2010 年 6 月 30 日按小时观测的 10 多指标的数据是准确无误的。2、 假设由于其他因素而引起10 多个指标数据的变化以及非正常波动可以忽略不计。3、假设第二问剔除由第一问求出的敏感度较低的指标对地震的敏感度可以忽略不计。4、假设地震的前兆指标的数据特征符合一定的概率统计分布。四、符号说明1( )it第i次地震的地震前的检测指标随时间的变化规律2( )it第i次地震的地震后的检测指标随时间的变化规律()2( )jt第i次地震的地震后的预测指标随时间的变化
9、规律D地震级别组成的集合E地震前兆测量指标集合ig各因素所占的权重地震前兆预测异常特征量r观测值的自相关系数K观测曲线( )Y t对应于时间轴的斜率五、模型的建立与求解5.1 数据的处理与分析5.1.1 、各指标数据的分析在地震活动性分析中有许多前兆性的数据指标,例如题中数据给出的氡值、水位、磁场强度等因素。 这些指标在不同的侧面反映了地震活动的各种特征。但在实际的预报中, 常常有些参数在一些中强以上的地震前出现比较明显。在正常情况下,也常常有些参数出现较明显的异常,而另一些参数并不出现异常。这些都给实际预报带来了困难。5.1.2 数据的预处理名师资料总结 - - -精品资料欢迎下载 - -
10、- - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页,共 20 页 - - - - - - - - - 4 由于题目中给出地震每小时测数据受其他因素受天气、气候等其他外在因素的影响造成的数据波动较大, 为了剔除偶然因素的影响, 本文采取了求日平均的方法来对数据进行预处理。如下图分别列举了2005/2006的 EW 人均值图像。由于篇幅有限,其他年份指标的图像都附在附录里面。图 1: 2005 年 EW 日平均变化趋势图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - -
11、名师精心整理 - - - - - - - 第 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 ,
12、现在若( )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(,)j
13、jPCov(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:各因素与地震敏感度的表但是考虑地震因素影响上述指标外,可能还有其它的因素影响上述指标。其中
14、查阅资料可得出:水位,水温,气压,气温受到外界影响最大,这样我们除上述几个指标剩下的指标基于与地震关联度从大到小排列:气氮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 与敏感度较低的低温年平均曲线,在途中可以看出在前兆指标因素发生剧烈变化时,会有地震的发生。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - -
15、 名师精心整理 - - - - - - - 第 6 页,共 20 页 - - - - - - - - - 7 图 3: 2005,2006 年 EW 日平均变化趋势图图 4: 2005、2006 年地温日平均变化趋势图5.2.2 综合指标模型方法及步骤用两个集合, 一是地震级别组成的集合, 用D表示,二是地震前兆车辆指标集合,用E来表示,且有12,mDD DD,12,mEE EE,每个指标因素都有 m个状态级,如此需要对个指标影响程度分敏感、一般、不敏感、很不敏感五名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 -
16、 - - - - - - 第 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)
17、 、 详细分析影响地震发生的因素, 并筛选出若干重要因子作为分析指标。(2) 、在分析地震各指标以及波形的基础上,综合得到地震的等级指标。(3) 、给出适当的隶属度iP 和权数ig 。(4) 、修正矩阵 U。确定分析指标地震发生前的检测指标较多, 而各个指标的敏感程度不同, 经过第一问的各因素的关联度分析,筛选出气氮、EW、NS 等六个指标。隶属度iP 和权数ig 的确定应用综合指标进行判别, 隶属度函数是一个关键, 用它来反映地震发生。 设N是评定集合上的一个模糊子集,则隶属度函数为uN。()iiuN DP 为隶属度,其值可以有专家评判方法来确定。假设等级划分与指标间基本呈线性关系,则根据第
18、 一 位 求出的 关 联度 来 确定隶 属 度 ,因 此 选 用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 页 - - - - - - - - -
19、 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
20、 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.
21、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 -
22、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
23、 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(1
24、0), 有上述地震发生前的平均值可得系数为:a1 a2 a3 a4=3.1119 0.4654 -4.6999 -0.4878; 5.2.3 地震前兆指标数量特征模型先查阅与地震相关的其他指标如震级,频率等,而后先分析单个指标的数量特征,并以电磁波EW为例,利用动态从属函数建立模型,对于气温,水位等其他指标,类推可得。而后利用震级,频率等指标的数学研究方法,将已知条件的10 个指标综合起来,利用非线性组组合建立模型。电磁波 EW 的数量特征研究。地震活动是一种复杂现象,各种观测数据之间缺乏对比性,难以利用统一的物理场进行描述。 为了将各种观测量的异常信息进行表达,可把各前兆异常量转化为无量纲量
25、以便进行对比。 各种地震前兆观测量虽然是不同的物理量,但其共同的特点都是随时间变化的数值,即都是时间函数, 各种地震前兆异常也都表现为各种观测量随时间的突出变化,异常形态虽多种多样, 但究其本质共同点都是观测曲线随时间的斜率变化。 因此,观测曲线随时间的斜率变化将是判断异常的重要指标。由以上分析从各种前兆数据中提取地震信息量的方法应具以下特点:(1) 、地震信息是一个无量纲量;(2) 、用前兆观测量Y(t )对时间轴的斜率变化K,反映观测值的速率变化,一般情况下, K值越大,异常量越大;(3) 、用前兆观测值的自相关系数r 来显示观测值质量的好坏,r 值越大,观测数据质量越高。原理和方法从地震
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年地震模型 2022 地震 模型
限制150内