形变数据中天然扰动在时频域的响应研究.pdfVIP

第29卷　 第1期 内 陆 地 震 Vol.29　 No.1 2015年　 　 3月　 　 　 　 　 　 　 　 　 　 　 INLAND　 EARTHQUAKE　 　 　 　 　 　 　 　 　 　 　 Mar. 　 2015 文章编号:1001-8956(2015)01-0065-06 中图分类号:P315.7　 　 文献标识码:A　 　 形变数据中天然扰动在时频域的响应研究① 1 2 2 1 王 宁 ,吴 云 ,张 燕 , 吕健 (1.中国地震局第一监测中心,天津300180;2.中国地震局地震研究所,湖北 武汉430071) 摘要:对形变观测中可能影响数据质量的天然扰动现象进行了分析。 不同的自然现象对形变数据 在时频域上的响应不同。 气压对分钟值数据的影响明显,其在时频域上的表现为优势频率低于 0.05,即周期大于20 min,且在气压变化明显区域,信号能量较大;台风对秒值数据的影响较大,其优 势频率集中在0.15~0.3 Hz;降水对整时值数据影响明显,其影响集中在低频部分。 关键词:形变数据;天然扰动;时频响应 doi:10.16256/ j.issn.1001-8956.2015.01.008 　 　 地应变、地倾斜等定点形变测量指标能直接反应地壳介质的微动态变化,捕捉到地壳介质 -6 -10 破裂前的力学变化信息,观测优势在 10 ~10 量级,适用于监测地震的中短和短临阶段的 [1] [2-4] 前兆 。 汶川地震发生后,郝晓光等在文献 中最先指出汶川地震存在约48小时的“震前 扰动”现象。 但是,在形变观测过程中,还存在着降雨降雪、台风、气压等自然观测条件以及人 为因素等干扰现象,这些干扰因子的存在影响了地震前兆信息的获取。 为此,本文中针对形变 观测过程中可能存在的降雨降雪、台风以及气压现象,对形变资料进行了分析,通过其在时频 域的形态研究揭示自然现象作为干扰因子对形变数据的影响。 1　 数据分析方法 自然界中普遍存在着变化的非平稳信号,这种变化主要是指信号振幅随时间的变化和信 [5] 号频率随时间的变化 ,地震信号是典型的非平稳信号,因而单一的运用传统的傅里叶变换 不能显示出地震资料中干扰信号以及地震信号随时间的变化,时频分析方法很好地解决了这 个问题。 时频分析方法得到的能量谱图揭示了信号频率以及能量随时间的变化情况,具有很 好的应用价值。 本文中采用了两种时频分析方法:小波分析以及S变换。 1.1　 小波分析 小波变换最早由法国的地球物理学家于上世纪80年代提出。 信号h(t)的小波变换定义为 1 ∞ t -b C(a,b) = a∫-∞ h(t)Ψ × ( a )dt　 . ① 收稿日期:2014-01-07;修回日期:2014-03-24. 课题项目:中国地震局第一监测中心科技创新主任基金(2014008). 作者简介:王宁(1988~),女,助理工程师(实习),硕士,2014年6月毕业于中国地震局地震研究所大地测量学与测量 工程专业,从事地形变数据处理与分析研究. E-mail:kingning3344@126.com 66　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 内　 陆　 地　 震　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 29卷 式中,a表示尺度因子,b表示平移因子,Ψ(t)为基本小波。 在这采用的是连续小波变换。 连 续小波变换是在傅立叶变换和短时傅立叶变换的基础上提出来的,时频窗口可自行调节,同时 刻变换,在高频部分时间分辨率较好,在低频部分频率分辨率也较好,因而适宜分析地形变数 据这类非平稳信号。 在连续小波变换中,同时进行简谐波和高斯函数的伸缩和平移。 基本小 波必须满足容许性条件,即 ∞ | Ψ(ω)