当前位置:首页 >> 数学 >>

地震模型Word版


地震预测模型

摘要

地震预报是当代自然科学领域里一个直接关系到人民生命安全和社会经济 的发展,是科学难度很大的前沿课题。近几十年来,人民的努力虽然取得了丰硕 的成果,预报的实践有些有限的成功。 在地震发生前有很多前兆性指标,如磁场强度、氡值、电压、雨量、水温等, 这些指标都从不同的侧面反映了地震活动的各项特征。在正常情况下,也常常有 些指标有明显的异常, 而令一些指标并不出现异常。这些都给实际的前兆指标数 据特征分析以及地震预报工作带来了困难。 本文针对地震影响因素多,数据多,联系弱的特点,分别建立了针对各个任 务的数学模型, 首先,处理初始数据时选取日平均变化来消除一天中数据的随机 因素的影响, 然后更具这些数据建立了地震发生前后各个指标的不同程度度量模 型解决了地震对指标的敏感度的度量问题,进而找出了 EW,气氡,水位,电压, 雨量这些衡量地震的主要因素, 并分析出了这些指标的重要程度,在解决找统计 量时, 利用上次任务所得的主要因素的组合来衡量一次地震的地震前兆指标数量 特征,并得出了电压达到 26.8754,气氡达到 17.685 左右时就有地震发生的可 能性, 考虑全面周到, 效果较好, 在模型改进中所提出的判别分析法, 科学有效, 对数据利用率较高。 而除此之外的其他几个指标与震级有关,这样就这些指标的 分析在任务三中分别建立了地磁前兆异常的动态从属函数模型以及非线性综合 模型, 进而找出了地震级数与这些因素之间的非线性关系,对于上述数学方法在 地震预测预报中的应用没有达到预期效果,须进一步研究改进。综合这些工作, 我们应用易语言编写了分析地震数据的平台,通过分析地震数据生成了报告。

1

一、 问题重述
地震是地壳快速释放能量过程中造成的振动。虽然预测地震是世界性难题, 但迄今科学界普遍认为,有可能反映地震前兆特征的指标可能不少于 10 个。已 经有专业仪器在多个定点实时按秒记录这些指标的数据, 期望通过对记录数据的 分析研究找到地震的前兆特征。 现已采集到某地 2005 年 1 月 1 日至 2010 年 6 月 30 日按小时观测的 10 多个 指标的数据,和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的 数据。 这些数据中隐藏着地震发生的前兆特征。 科学地截取这些数据的有用片段, 对数据进行合理地预处理, 用数学方法揭示地震前兆的数据特征,是一项很有意 义的研究工作。 题给数据中的这 10 多个指标,究竟哪些与地震的发生有关,有何种关系, 是单一关系还是复合关系;除这 10 多个指标外还有哪些因素及含题给指标在内 的哪些指标的哪种数学模型更能反映地震的前兆特征等等,人们迄今仍不很清 楚,需要进行深入地研究。地震数据的观测是持续进行的,随着时间的推移数据 的规模会不断扩大。从中挖掘地震的前兆特征,必须有合理的数学模型,也必须 有科学高效的算法分析平台。因此,请参赛队结合附件中给出的实际记录数据, 尝试完成以下任务。 任务一:分析数据特征,建立数学模型以度量各指标对地震发生的敏感程度。 任务二:构造由某些或全部指标构成的综合指标,使其尽可能地集中反映地 震发生前的数据特征的统计规律。 任务三:结合题给数据,广泛查阅与地震相关的其它指标的数据和分析方法, 建立数学模型来研究地震发生前的数量特征。 任务四:将前述各项任务的计算程序集结成地震数据分析平台,使其能够完 成形如题给数据的其它地震数据的分析,并能自动输出前述任务的重 要的分析结果。 任务五:对于进一步的研究设想写一篇切实可行的报告。

二、问题分析
地震是地下岩层受应力作用错动破裂造成的地面震动, 是一种破坏性极强的 自然灾害,是自然灾害之首恶。地震可以在很短的时间内使一座城市夷为平地, 使无数的家庭支离破碎。 为了有效地减小地震带来的损失,对地震进行预测则很 重要。 本题中给出了地震发生前 2005 年 1 月 1 日至 2010 年 6 月 30 日按小时观测 的 10 多个指标的数据和该地区该时期内已发生地震的时刻、经纬度、震级及震 源深度的数据。 这些数据中隐藏着地震发生的前兆特征。在地震发生前有很多前 兆性指标,如磁场强度、氡值、电压、雨量、水温等,这些指标都从不同的侧面 反映了地震活动的各项特征。在正常情况下,也常常有些指标有明显的异常,而 令一些指标并不出现异常。 这些都给实际的前兆指标数据特征分析以及地震预报 工作带来了困难。 2.1、问题的性质 本文主要是对所给数据进行处理,用数学方法揭示地震前兆的数据特征。 2.2、解决问题的难点和关键
2

1、如何确定前兆个指标对地震发生的敏感程度。 2、如何构造地震发生的前兆综合指标,尽可能的揭示地震发生前的数据特 征的统计规律。 3、构建地震数据分析平台。

三、模型假设

1、假设地震发生前 2005 年 1 月 1 日至 2010 年 6 月 30 日按小时观测的 10 多指 标的数据是准确无误的。 2、 假设由于其他因素而引起 10 多个指标数据的变化以及非正常波动可以忽略不 计。 3、假设第二问剔除由第一问求出的敏感度较低的指标对地震的敏感度可以忽略 不计。 4、假设地震的前兆指标的数据特征符合一定的概率统计分布。

四、符号说明
?1i (t )
i ?2 (t )
? ( j) 2

第 i 次地震的地震前的检测指标随时间的变化规律 第 i 次地震的地震后的检测指标随时间的变化规律 第 i 次地震的地震后的预测指标随时间的变化规律 地震级别组成的集合 地震前兆测量指标集合 各因素所占的权重 地震前兆预测异常特征量 观测值的自相关系数 观测曲线 Y (t ) 对应于时间轴的斜率

?

(t )

D E

gi
? r K

五、模型的建立与求解 5.1 数据的处理与分析
5.1.1、各指标数据的分析 在地震活动性分析中有许多前兆性的数据指标,例如题中数据给出的氡值、 水位、磁场强度等因素。这些指标在不同的侧面反映了地震活动的各种特征。但 在实际的预报中, 常常有些参数在一些中强以上的地震前出现比较明显。在正常 情况下,也常常有些参数出现较明显的异常,而另一些参数并不出现异常。这些 都给实际预报带来了困难。 5.1.2 数据的预处理
3

由于题目中给出地震每小时测数据受其他因素受天气、 气候等其他外在因素 的影响造成的数据波动较大, 为了剔除偶然因素的影响,本文采取了求日平均的 方法来对数据进行预处理。如下图分别列举了 2005/2006 的 EW 人均值图像。由 于篇幅有限,其他年份指标的图像都附在附录里面。

图 1: 2005 年 EW 日平均变化趋势图

4

图 2: 2006 年 EW 日平均变化趋势图

本数据预处理使用了均值结构模型消除了其他偶然因素的影响,使地震前兆 指标数据更加准确的反应地震发生的前兆信息。

5.2 模型的建立与求解
5.2.1 指标敏感度模型 设变量 x 的变化规律如下所示:

??1 (t ) t ? ? t 1 t , ? , 0 ? x ? ? (t ) ? , 1 ?? 2 (t ), t ? ? t 0 t ? ?
(1) 求 t 0 前后 ?1 (t ) 与 ?2 (t ) 的不同敏感程度。 现已知一函数序列 ? (i ) 2 (t ) , t ? ? t0 , t1 ? 表示第 i 年 (i ? j ) 的 ?2 (t ) 的变化规律, 现在结合 ? ( j )1 (t ) 来求 ? ( j )2 (t ) ,现在若 ?
? ( j) 2 (t ) 已知,如何求出 ? ? ( j) 2

(t ) 与 ? ( j )2 (t ) 的

关联程度。若关联程度较大,则 ? ( j )2 (t ) 与 ? (i )1 (t ) 不同程度表示: 即定义 ? ( j )1 (t ) 与 ? ( j )2 (t ) 的不同程度为:

5

p?

1

??

? ( j) 2

(2)
( j) 2

:?

分析变量 X 1 与 X 2 的关联程度 Cov( X1 , X 2 ) ? E( X1 ? X 2 ) ? E( X1 ) E( X 2 )
1 Cov(?
?
? ( j) 2

P?

(3)
,?
( j) 2

)

则由分析可得该问题便为 ? ( j ) 2 如何求的问题。 现已知 ? ( j ) 2 (i ? j ) 与 ? ( j )1 (t ) ,如何根据已知的以上信息求出 ? ( j ) 2 , 已知 ? (i )1 (t ),(i ? j ) 的变化规律,可知令
?

? ( j ) 2 (t ) ? E (?

(i ) 2

)

E ( E?1( i ) (t )) , E (?1( j ) (t ))

(i ? j )

(4)

模型求解: 由上述结论我们分析数据得出各个因素在地震发生前后的敏感程度分别如 下: 电压 0.2 EW 1.8 NS 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:主要地震前兆指标敏感度的表

标准化后有: 气氮 0.45

EW 0.29

NS 0.13

雨量 0.08

电压 0.03

地温 0.01

表 3:主要地震前兆指标敏感度归一化表

如下图分别画了敏感度较高的 EW 与敏感度较低的低温年平均曲线,在途中 可以看出在前兆指标因素发生剧烈变化时,会有地震的发生。

6

图 3: 2005,2006 年 EW 日平均变化趋势图

图 4: 2005、2006 年地温日平均变化趋势图 5.2.2 综合指标模型 方法及步骤 用两个集合,一是地震级别组成的集合,用 D 表示,二是地震前兆车辆指标 集合,用 E 来表示,且有 D ? [ D1, D2 ??? Dm ] , E ? [ E1 , E2 ??? Em ] ,每个指标因素都 有 m 个状态级,如此需要对个指标影响程度分敏感、一般、不敏感、很不敏感五
7

个状态。 有 n 个指标因素分别用 U1 ,U 2 ???U n 来表示,其中 Ui ? [U i1 ,Ui 2 ???Uim ]T ,于是 可得到敏感指标分析数学模式
, ?U1 1 U 2??? U ?1 1 n ?U , U ??? U ? D ? ? 1 2 2 2 n ?2 ?? ? ? ? ??? ?U1m , U1 2 Un m ?

(5)

有了这个模型,分别给予地震级别的隶属度函数值 P (i ? 1, 2, ???m) 再根据指 i 标因素间的关联度及其重要性,分配权数

gi (i ? 1, 2, ???n)
这样就得到了集合 D 上的模糊关系:
, 2 1 ? g1 P1 g P ??? g n P 1? ? g P , g P ??? g P ? 2 U ?? 1 2 2 2 n ? (6) ?.... ? ? ? ? g1 Pm , g 2 Pm ??? g n Pm ? 基于地震的综合指标步骤有下列步骤: (1) 详细分析影响地震发生的因素, 、 并筛选出若干重要因子作为分析指标。 (2) 、在分析地震各指标以及波形的基础上,综合得到地震的等级指标。

(3) 、给出适当的隶属度 Pi 和权数 gi 。 (4) 、修正矩阵 U。 确定分析指标 地震发生前的检测指标较多,而各个指标的敏感程度不同,经过第一问的各 因素的关联度分析,筛选出气氮、EW、NS 等六个指标。 隶属度 Pi 和权数 gi 的确定 应用综合指标进行判别,隶属度函数是一个关键,用它来反映地震发生。设
N 是评定集合上的一个模糊子集,则隶属度函数为 uN 。 uN ( Di ) ? P 为隶属度, i

其值可以有专家评判方法来确定。假设等级划分与指标间基本呈线性关系,则根 据 第 一 位 求 出 的 关 联 度 来 确 定 隶 属 度 , 因 此 选 用
P ? [ 0 . 4 5 , 0 . 2 9 , 0 ,为了方便起见将隶属度扩大一百倍后可得:] . 1 3 , 0 . 0 8 , 0 . 0 3 , 0 . 0 2

P ? 45 , P ? 24 , P3 ? 13 , P4 ? 8 , P ? 3 , P ? 0.02 , 1 6 2 5
权数 gi 的确定是人给定的,选取的方法有很多,本文选用特费尔法,并确定 气氮、EW、NS、雨量、电压地温的权数分配为:
g ? [1,0.95,0.85,0.80,0.70,0.60]
8

这里 g 为 E ? ( E1, E2 , ???, E6 ) 上的模糊子集,按模糊数学的常规技法

g? 1

E1

?0 . 9 5 ? 0 . 8 5 ? E2 E3

0.80 ? E4

0.70 ? E5

0 . 6 0 (7) E6

模糊关系矩阵 U 的确定 根据 P ? [ P , P , ???, P ] , g ? [ g1, g2 , ???g6 ] 就可以得到在集合 D 上所需的模糊 1 2 6 关系矩阵:
? g1 P , g 2 P ??? g n P ? 1 1 1 ? g P , g P ??? g P ? n 2 ? U ?? 1 2 2 2 ?.... ? ? ? ? g1 Pm , g 2 Pm ??? g n Pm ? 由第一问可得各因素的敏感度为: 电压 EW NS 地温 水位 气温 气压 0.2 1.8 0.8 0.1 2595.3 2.3 4.3
表 4:地震前兆指标敏感度的表

(8)

水温 40.2

气氮 2.8

雨量 0.5

每次地震的等级: D=[2.3 3.6 3.1 3.1 3.2 3.0 4.4]; 稳定性分析数学模式: 0.0761 5.0850 0.0589 1.3422 0.4320 0.9064 0.3576 0.8787 0.0456 1.8595 0.0715 0.8196 0.0719 0.6806 地震指标权重 -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 26.9050 -2.0332 4.7336 17.5866 8.8819 27.0236 -3.0234 3.3379 17.5680 12.0243 27.0298 -2.7623 3.6132 17.5712 12.4378
9

0.3290 1.2721 0.5601 0.6075 1.6741 0.4019 0.1246

0.0920 0.0548 0.1643 0.1105 0.1224 0.0172 0.0048

3.2349 2.2143 3.7511 3.2558 1.8006 2.7385 3.0955

0.1262 0.5503 0.9895 0.2348 0.2730 0.8949 0.0307

表 5:每次地震前兆稳定性系数

0 0.0277 0 0

26.7648 26.7152 26.7021

0.1313 -1.3634 1.8757

6.9721 6.4789 8.5844

17.6188 17.6202 17.6550

33.2737 17.9615 27.2101

0 0.6602 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

?1 ? ? (ti ) ? 1 ?
式中:

?
ki ri

?1

(9)

Ki ?

?t j ? y j ?n?t j y j
j

(

?t j )
j

j

2

?n

?t j2
j

j

(10)

ri ?

?t ?t
j 2 j

j

yj ?

1 (? t j ? y j ) n j j
1 2

(11)
1 2

?

1 (? t j ) 2 n j

?y
j

2 j

?

1 (? y j ) 2 n j

y j 是前兆观测随时间变化的数列, t j 是相应的时间序列
公式(1)中的 k 物理意义是观测曲线 Y(t)对于时间轴的斜率,反映观测速 率的变化;r 是观测值的自相关系数,反应观测值内在质量的好坏。从(1)式 可知,只有 k 和 r 都大时 ? ,才可能大。在 [0,1] 内取值,表示在所研究的资料在
t j 时段属于地震异常的程度。 ?i ? 0.5 作为模糊界点,只有当 ?i ? 0.5 时,才视为

异常,反之,即为正常状态。? 是一个经过学习得到的常数。它的选择原则和方 法如下: (1) 使得“异常”频次降低; (2) 使“异常”和“地震”之间由最好的相关性; (3) 在满足上述两个原则下要尽可能的小; (4) 选择的“地震”样本要视监视区的活动水平和预报目标而定; (5) 选择后的 ? 值,应保持稳定,只有出现预报目标与实际情况有明显 差异时,才按照新的样本要求进行修改,以保持第二条原则的实现。 模型的检验 选择样本 从属函数中的 ? 的选择由样本学习而决定,不同学科不同地域选择的参数可 能会不同, 需要在实践中不断的完善,一些前兆场可直接根据原始观测数据分析 其变化,而地磁场本身是一个变化场,且磁层、电离层强烈扰动时引起的地磁场 变化非常剧烈, 远大于地震异常引起的变化。从属函数法以曲线斜率作为判断异 常的组要指标,如果直接分析某一点的从属函数值,会频繁出现“异常” ,判断 食物。 同样是前兆手段, 地磁学在应用从属函数方法时首先对原始地磁场资料预 处理,才由利于寻找震磁关系。 电磁波由地磁场产生,而地磁场值包含核远源场,外源场及地下局部磁异常 源场,提取电磁波 EW 的问题是如何找出因地磁场的地下局部磁异常发生变化的 问题, 并以此判断该变化量是否是由地震引起的局部异常。核源场是一个相对稳 定场, 外源场的变化受诸多因素影响,而且其变化幅度远大于因地震引起的异常 量,因此应用在应用资料预报之前,首先应剔除外源场引起的地磁场变化。任一
11

点的地磁场垂直强度 Z 由以下三部分组成:

Z ? ZC ? Ze ? Z N
C-地核源场,e-外源场,N-地下局部磁场。 由此可推得下面公式:

(12)

EW ? EW C ? EW e ? EW N

(13)

下面采用均值法消除偶然因素如太阳黑子等造成的突变, 并用对应一年中的 均和几年的均值作差以消除外源场和核源场的变化而引起地磁场的突变

?Z ? Z ? Z ? Z N ? Z N
由此可推得:

(14)
N

?EW ? EW ? EW ? EW N ? EW

(15)

正常情况下地点确定后局部磁场在很小的范围内波动;由于地震孕育而导致 局部磁场发生变化时,就在孕震区产生局部磁异常,这个范围比较小,即使两地 都处于孕震区, 异常将由于距离震源距离不同而有所差别, 最终导致 EW 的变化, ?EW 代表局部磁异常的相对变化。 从属函数的计算 对应年份观测数据和各年份数据的平均作差消除基本场与外空场而引起的 EW 变化,用于消除外源场和核源场的变化而引起地磁场的突变,取各年份的均值 是为偶然因素而引起的电磁波突变 经过多年样本学习,确定 ? 值,研究电磁波随时间的变化情况。 每次地震每个指标相对应 ? 的值
第一次地震 第二次地震 第三次地震 第四次地震 第五次地震 第六次地震 第七次地震

电压 0.5000 0.9883 0.9984 0.9786 0.9994 0.9927 0.9630

EW 0.9867 0.5000 0.9613 0.9916 0.9837 0.9950 0.9989

NS 0.9666 0.9669 0.9255 0.9882 0.5000 0.9943 0.9949

地温 0.9994 0.5000 0.7123 0.9889 0.9923 0.9613 0.9982

气氡 0.9988 0.9997 0.9937 0.9999 1.0000 0.9999 0.5000

雨量 0.9916 0.8144 0.9236 0.8648 0.9836 1.0000 0.5000

表 7:每次地震前兆预测异常特征量

? 值的分析: ? 的取值范围为 [0,1] ,0.5 为判别地震发生的界点。
当 ? ? 0.5 时,则地震没有发生的可能。 当 ? ? 0.5 时,地震有发生的可能, ? 值越大,则地震发生的可能性越大。 每个指标的 ? 值
12

1.0e-003 * 0.0117 0.1812 0.3049 0.0021 0.0142 0.0736

图 5: EW 两年相对变化趋势图

5.2.3 地震前兆指标数量特征模型 针对地震预测研究中, 变量的指标较多,数据携带的信息较弱且预测变量是 非线性关系, 因此造成预测结果稳定性不强的问题, 构造了非线性组合预测模型, 解决地震预测的上述问题,应用于本题中的地震前兆指标数量特征。 即使数据携带的信号很弱, 但是我认为这些数据足以通过统计方法来探求深 层次的东西, 因为数据不服从标准的平稳时间序列和线性模型,所以我们要小心 谨慎地选择统计方法。 下面构造的多变量非线性模型在一定程度上解决了上述问 题且在地震预报应用中取得了较满意的效果。 非线性组合模型 非线性组合模型是将多个看用于预测的一元线性模型进行线性组合,组成非 线性组合模型,其优点是可以增加预测稳定性。 非线性组合模型的应用背景 设 y 是预测变量, x1 , x2 ??????x p 是影响 y 的指标。 (1) 指标数据携带信号较少、 、 较弱, 需要选择携带信号较强的指标用于预测。 (2) 、对于选定的指标需要建立较好的非线性预测模型。 (3) 、为增强预测的稳定性,最后构建非线性组合模型用于预测。 最佳指标和最佳曲线选择 对于每一个指标 xi (i ? 1, 2, ???, p ) ,对应预测变量 y ,预先设计若干非线性模 型,并进行曲线拟合,根据统计检验和拟合优度选择最优指标和最佳曲线。
13

设所选的最优指标为 xi 相应的最佳曲线函数是 f ( xi )(i ? 1, 2, ??????, k ) , k 是最 优指标和最佳拟合曲线个数。 非线性组合模型
y ? ? wi fi ( xi )
i ?1 k

(16)

其中 y 是预测变量, fi ( xi ) 是选择出的最佳非线性预测曲线, xi 是自变量指 标, w i 是相应的权重, ? wi ? 1 。
i ?1 k

非线性组合模型参数 wi 的估计 一般参数估计的方法可选择: (1)等权平均法; (2)加权最小二乘法; (3)方差-协方差方法; (4)最小二乘法; 非线性组合模型的应用 在本题中, 选择最优指标和最佳曲线, 通过最小二乘法估计最佳曲线的参数, 并构建非线性组合模型。 最佳指标和最佳曲线的选择 初步选择影响最大的若干地震参数指标 xi 等,见表一 为选择最优指标和最佳拟合曲线设计了七种曲线模型,见表二。 选择拟合方式以及统计检验, 选出 6 个最优指标以及他们相应的最佳曲线。见表 三。 2.2、非线性组合模型
y ? ? wi fi ( xi ) ?
i ?1 ? k

1 6 ? ? yi 6 i ?1
1 , (i ? 1, 2, ???,6) 6

考虑到增加预测的稳定性,参数估计采取等权平均法 wi ? 经过 Matlab 编程得出了地震各因素的权重如下: 电压 EW NS 地温 气氡 -0.0000 -0.0029 0.0001 0.0057 -0.0014 画出部分地震前兆指标与地震发生的关系图如下:

雨量 -2.3620

14

图 6:相对震级图

5.2.4 地震数据分析平台的设计 把前三问的程序集成,基于易语言编程,系统开发使用 Windows2000server、 MS-SQLServer2000 数据库,对地震目录、前兆观测数据、分析报告等地震信息 进行了系统管理, 组成了地震数据分析平台,本平台在实现了地震数据的数字化 和自动化管理功能的基础上,系统运用数据挖掘的方法对地震预报进行了尝试, 使其能够完成形如题给数据的其它地震数据的分析, 并能自动输出前述任务的重 要的分析结果。其主界面如下图.

图 7:地震数据分析平台主界面 基于此地震数据分析平台,分析了 2009 的地震指标数据,得出的结果是地
15

震发生时间为 10 月 31 日,震级为 4.1 级,在误差范围内,实现了较好的地震前 兆数据分析与得出重要的分析结果。

六、模型的推广
1、在地震分析中经常要分析各个变量间相互依存关系,本文模型只能分析每一 个变量对因变量的直接影响,且要求其余自变量之间相互独立。事实上,变 量之间的关系错综复杂,一个自变量可能对一个因变量有直接影响,也可能 通过其他的自变量对因变量有间接地影响,显然本模型无能为力。 2、可选用通径分析来求解各个变量之间的关系,在通径分析模型中,把影响其 他变量的变量叫做原因变量,被影响的变量叫做结果变量。

七、关于进一步研究报告
地震预报预测之艰难不言而喻,然而通过此次研究分析发现,地震发生前并 非不显示任何蛛丝马迹, 笔者基于促进其预测预报研究的目的提出以下建议和设 想: 确保各指标观测数据的真实性、准确性。数据是供人们研究的平台,其真实 性、准确性直接关乎研究方向的正确性及研究结果的实效性。 优化数据挖掘方法。从各指标对事件的函数图像可知,各指标值振荡频繁, 然而不是所有的振荡都是因为地震的发生而引起的, 如电磁波的振荡容易受太阳 黑子等因素影响而发生明显振荡。如何采取科学、合适的数据挖掘方法关系到所 得模型预测预报的准确率, 我们所采取的是均值结构模型以达到消除偶然误差的 目的,其效果显著,然而由于原理上的限制,还是有所偏差,所以笔者设想能否 利用数理统计方法进行估计,以使偏差尽量降低。 研究地震发生时各指标异常机理。研究过程发现,各指标异常会出现一定的 提前量或滞后性, 如何估计这些提前量关乎预测预报的准确性,只有通过研究各 指标因地震变化的机理才能用数学方法较为准确的估计这些提前量。 优化数学研究方法。 本模型在求解任务一时采用均值结构模型达到较良的消 除误差的效果;采用综合指标法求解任务二,考虑全面周到,效果较好,在模型 改进中所提出的判别分析法,科学有效,对数据利用率较高。在任务三中采用各 年份指标值和历年各指标均值作差以消除地核、 外空对电磁波的影响, 巧妙科学, 适宜推广使用。 任务三中采用非线性组合模型,较好的解决数据携带信号较弱的 缺点。限于知识、能力,对于上述数学方法在地震预测预报中的应用没有达到预 期效果, 须进一步研究改进。 由于地震的发生及受其影响的各指标间的复杂关系, 很难采用现有的技术探测,笔者试想,能否应用 Bp 神经网络解决此难题,设定 科学的学习规则,利用神经网络的自我学习特性,达到预测预报的目的。

参考文献
16

[1] [2] [3] [4] [5] [6] [7]

徐文耀,我们地震观测研究发展,地球物理学报,997,40:217。 冯德益,等模糊数学方法应用,北京,地震出版社。 郭德科, 强震前兆异常的模糊识别函数与中弱地震的关系, 地壳形变与地震。 张岩波,潜变量分析,北京,高等教育出版社。 王清河,随机数据处理方法(第三版) ,石油大学出版社。 张继红、乔慧珍,地震前兆异常的动态从属函数,地球物理学进展。 沈良峰,地震稳定性分析的综合性指标模型,建筑科学。

附录
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),':',num2st r(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((4 0*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(e r((40*i-39):(40*i),:)); r1(i,:)=r1(i,:)./(sqrt(sum((1:40).^2)-(1/40)*sum((1:40))^2)*sqrt(sum(
17

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' 'WT' 'AID' '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 i<10&&j<10 t=strcat(year{qs},'-0',num2str(i),'-0',num2str(j),'%'); t=strcat(char(39), t,char(39)); t=[' ',t]; else if i<10&&j>9 t=strcat(year{qs},'-0',num2str(i),'-',num2str(j),'%'); t=strcat(char(39), t,char(39)); t=[' ',t]; else if i>9&&j<10 t=strcat(year{qs},'-',num2str(i),'-0',num2str(j),'%'); t=strcat(char(39), t,char(39)); t=[' ',t]; else t=strcat(year{qs},'-',num2str(i),'-',num2str(j),'%'); t=strcat(char(39), t,char(39)); t=[' ',t]; end end end
18

s(j+31*(i-1))=fetch(coona, strcat('use earthquck avg(',yinsu{qw},') from years',year{qs},'$ 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 2009]; S = load('matlab.mat'); er=zeros(7,6); for i=1:7

select where

u=strcat('S.t',num2str(year(i)),'(',num2str(RI(i)-30),':',num2str(RI( i)),',[1 2 3 4 9 10])'); data=eval(u); data=mean(data); er(i,:)=data; end 4、相关系数的 Matlab 程序 %计算每次地震中指标相关程度; RI=[6555,8018,7337;4346,5834,4948;1418,2882,1932 ;2162,2882,2259;2186 ,3650,3265 ;3650,5114,4473;6554,8028,7401]; year=[2000511 2000611 2000711 2000711 2000811 2000811 2000911]; vv=zeros(7,10); for sd=1:7 S = load('matlab.mat'); t2000511=S.t2005(RI(sd,1):RI(sd,2),1:10); t2000611=S.t2006(RI(sd,1):RI(sd,2),1:10); t2000711=S.t2007(RI(sd,1):RI(sd,2),1:10); t2000811=S.t2008(RI(sd,1):RI(sd,2),1:10); t2000911=S.t2009(RI(sd,1):RI(sd,2),1:10); tem1=strcat('t',num2str(year(sd))); s=(t2000611+t2000711+t2000811+t2000911+t2000511-eval(tem1))/4; yu=RI(sd,3)-RI(sd,1); tem=strcat('t',num2str(year(sd)),'(1:',num2str(yu),',:)'); k=mean((t2000511(1:yu,:)+t2000611(1:yu,:)+t2000711(1:yu,:)+t2000811(1 :yu,:)+t2000911(1:yu,:))/5)./mean(eval(tem)); y=s*0; for p=1:10
19

y(:,p)=s(:,p)/k(p); end y=y(yu:end,:); rt=t2000511(yu:end,:); vv(sd,:)=std(y-rt,0,1); end b=mean(vv);

20


相关文章:
地震科普知识手抄报 word电子小报模板
地震科普知识手抄报 word电子小报模板_公务员考试_资格考试/认证_教育专区。电子小报模板手抄报模板 中国板报网 www.cnbanbao.cn 震多发地带: 东南沿海地区,就是亚欧...
地震工程学概论复习(word版)_图文
地震工程学概论复习(word版)_工学_高等教育_教育专区。老师课件,精心整理。...功率谱密度函数 为: 2)有限带宽白噪声过程/模型: 理想的白噪声具有无限带宽,...
地震受灾报告Microsoft Word 文档
地震受灾报告Microsoft Word 文档 隐藏>> xx 乡遭受地震灾害袭击的情况报告 2012 年 9 月 7 日 11 时 19 分至 13 时,我乡遭受地震灾 害袭击,经受灾村上报...
汶川地震Microsoft Word 文档
汶川地震Microsoft Word 文档 2008年汶川地震的一次亲身经历2008年汶川地震的一次亲身经历隐藏>> 我所经历的汶川地震三年前的今天,一个平常的不能再平常的午时 2:28...
地震总结oft Word 文档
地震总结oft Word 文档 总结总结隐藏>> 西沟小学防震演练总结为了预防突如其来的地震灾害,保证师生安全脱险,把 损失降低到最低限度,我校在组织了一次防震演练,其效...
云南地震Microsoft Word 文档
云南地震Microsoft Word 文档 earthquakeearthquake隐藏>> 1. Where did the earthquake heppen? 2. How many people are killed and injured? 3. How strong is...
地震演练 Microsoft Word 文档
江南实验学校地震逃生演练方案一、演练目的 通过地震应急演练,使全校师生掌握应急避震的正确方法,熟悉震后我校紧急疏散的 程序和线路,确保在地震来临时,我校地震应急...
新建Microsoft Word日本地震最新数据.doc
新建Microsoft Word日本地震最新数据.doc 日本地震最新数据日本地震最新数据隐藏>> 日本大地震满百日 灾民生活困苦经济短时难复原 2011 年 06 月 18 日 03:07 来...
地震后救援物资供给模型
地震模型Word版 20页 5财富值 地震模型波形模拟 5页 5财富值 地震人员伤亡快速...地震后救援物资供给问题摘要地震作为对人类最大的自然危害之一,震后救援是一个非...
张靓颖 地震主题 Word 文档 (3)
张靓颖 地震主题 Word 文档 (3)_计算机软件及应用_IT/计算机_专业资料 暂无评价|0人阅读|0次下载|举报文档 张靓颖 地震主题 Word 文档 (3)_计算机软件及应用_...
更多相关标签: