密山地电场数据小波包去噪研究.pdfVIP

地 震 地 磁 观 测 与 研 究 第 36 卷 第 1 期 SEISMOLOGICAL AND GEOMAGNETIC Vol .36 N o . 1 20 15 年 2 月 OBSERVATION AND RESEARCH Feb . 20 15 doi: 10. 3969/j. issn. 1003-3246. 2015. 01. 014 密山地电场数据小波包去噪研究 1), 2) 1) 1) 3) 梁 跃 董桂菊 崔天时 张东海 1）中国哈尔滨 150030 东北农业大学 2 ）中国黑龙江 157009 黑龙江省牡丹江地震台 3 ）中国黑龙江 158300 黑龙江省密山地震台 摘要 选取 2009 年 11 月 2 1 日密山市知一镇地震发生当天密山地 电场数据 ，采用小波 包阈值去噪算法 ，对地 电场数据进行去噪研究 ，还原地震发生后地 电场数据变化信息 ， 并结合去噪效果及误差分析结果 ，选出适合密山地 电场数据去噪的最优小波函数。 关键词 地 电场 ；小波包 ；小波 ；阈值去噪；小波函数 0 引言 小波变换 自出现以来因其优越的时频分析特性，已经广泛应用在信号分析、图像处理、 地震勘探数据处理等诸多领域。小波包变换是由小波变换进一步发展而来的，同小波变换 相比，小波包变换具有更高的时频分辨率及更加优越的去噪性能，尤其对于非平稳信号， 小波包去噪已经成为信号处理的主流方向之一。魏红梅等（2008 ）将小波包去噪方法应用于 地震信号的预处理中，削弱了地震信号中的高、低频干扰信号。邱颖等（2009 ）应用小波域 阈值滤波方法对叠加在大地电场信号上的地电阻率人工供电干扰信号进行滤波处理。王培 茂等（20 11）提出基于匹配算法的小波包去噪方法，很好地压制了地震信号中的白噪声。 地 电观测资料信息量大，各种高频信息丰富、干扰因素复杂（张彩艳等，20 11）。密 山地 电场观测数据一直存在类似尖脉冲干扰信号 ，可能是由多种不同因素叠加到一起造 成的（高研等，20 12 ）。2009 年 11 月 2 1 日密山市知一镇发生 M L 3.7 地震 ，密山地震台测 定震 中距为 27.8 km ，地震发生前后密山地 电场数据均未体现任何变化信息。本文以此 次地震为例 ，采用小波包阈值去噪方法 ，还原地震发生后地 电场数据中的变化信息，并 选出适合密山地 电场数据去噪的小波函数。 1 原理及方法 自 1822 年傅里叶变换发表以来，傅里叶变换一直是信号处理领域中应用广泛的分析 手段。但是，由于傅里叶变换是一种全局变换，无法表述信号的时频局部特性。尽管窗口 傅里叶变换的提出实现了信号在时频域内局部化的联合分析，但因其时频窗的宽度是固定 不变的，无法根据信号的频率变化进行调整。小波变换继承和发展了窗口傅里叶变换的局 部化思想，克服了窗口大小不随频率变化的缺点，在低频部分具有较低的时间分辨率和较 作者简介：梁跃 （1988—），男，黑龙江五常人，硕士，主要从事地震监测工作，E-mail: 381665839@ 基金项 目：中国地震局 20 14 年度 “地震监测、预测、科研三结合”课题 （14080 1）资助 本文收到 日期：20 14-09-29 第 1 期 梁 跃等 ：密山地 电场数据小波包去噪研究 85 高的频率分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率，是一种比较理 想的信号处理工具。小波包变换是基于小波变换的进一步发展，其优一点在于，小波包分 解对小波分解中没有细分的高频部分做进一步分解，具有更高的时频分辨率，图 1 和图2 分别为小波及小波包的 3 层分解结构图。 s s a 1 d 1