二次小波变换的阈值去噪算法及其在地震信号处理中的应用.pdfVIP

石油天然气学报 (江汉石油学院学报) 2010年6月 第32卷 第3期 JournalofOilandGasTechnoLogy(J．JPI) Jun．2010 Vo1．32 No．3 二次小波变换的阈值去噪算法 及其在地震信号处理中的应用 张 涛 (长江大学信息与数学学院，湖北荆州434023) 周秀君 (广东纺织职业技术学院，广东佛山528041) 杨世 明 (长江大学地球物理与石油资源学院，湖北荆州434023) [摘要]设计 了一种二次小波变换 的阈值去噪算法 ，利用该算法可 以合理去掉地震信号中面波分布 区域 的 低频成分，保留中、高频有效信息，从而既压制 了面波 ，又减少低频有效波的损失。该方法较之传统的 域阈值去噪方法有更好 的去噪效果，且对信号保真度好 ，基本不会使有效信号发生畸变，不会损失有效 信号 。 [关键词]地震信号；二次小波变换 ；阚值 ；去噪 [中图分类号]P631．4；TP391．41 [文献标识码]A [文章编号]1000—9752(2010)03—0224—04 提高地震信号的分辨率是地震勘探所要解决的一项主要任务，而提高地震信号的信噪比是提高地震 信号分辨率的先决条件。要提高地震资料的信噪比，需要去除地震资料中的随机噪声 ，而随机噪声是没 有统一的规律 ，在整张记录上随机出现 ，比较难 以去除[1q]。因此如何有效消除随机噪声是人们一直追 求的 目标 。为此，笔者提出了一种基于二次小波变换的小波域阈值去噪方法，该方法不仅利用信号和噪 声在小波变换下奇异性截然不同的表现特征来去除噪声，还利用地震信号道间相关性进行去噪处理 。该 方法在进行相关分析时对不同频率分量选择不同宽度的时窗进行相关分析，并根据地震剖面情况灵活选 取进行相关分析 的小波系数尺度个数 、时窗的宽度、倾角的个数等参数 ，因此该方法较之传统 的域 阈值 去噪方法有更好的去噪效果，且对信号保真度好 ，基本不会使有效信号发生畸变 ，基本不损失有效 信号 。 1 算 法 算法具体步骤如下： (说明：z为含噪信号；HO、F0分别表示分解及重建低通滤波器 ；HI、F1分别表示分解及重建高通滤 波器 ；Y为去噪后的信号 。) Stepl 读入含噪信号 NN]一Ix ， ，…，XNN]。 输入滤波器组 ： HOVK]一 [hol，h02，…，厅0K] HI[KK]一 [hl，hi2，…，1KK] FOEK]一 [fO，fo：，…，_厂O] FI[KK]一 [fl，fl，…，f1KK] 计算 L—max{K，KK}。 Step2 分解含噪信号 x[NN]。 Step2．1 x[NN]周期延拓至： zoF]一 [zNN—f+2，…，zNN，．；171，57／2，…，zNN，z】，X2，…，z[一】] [收稿日期]2009—05—2O [基金项目]国家 自然科学基金项目 D0206)；石油科技中青年创新基金项目 (2002{70104)。 [作者简介]张涛 (1978一)，男，2002年江汉石油学院毕业 ，硕士，讲师，现主要从事最优化理论与算法研究工作。 第32卷第3期 张涛等：二次小波变换的阈值去噪算法及其在地震信号处理中的应用 Step2．2 xOENN+2L一2]延拓至： 0oL —LLO二 ，XNNL+2，…，XNN，z1，X2，…，XNN，1，z2，…，Ll，0，…，0，] K一’】个 K—V1个 Step2．3 低通滤波： O口与H0口的卷积 ∞[]：一0 Step2．3．1 i—J，h一0，ca[力 一 0 Step2．3．2 caEJ]一xOOEi]HOEK—h一 1]4-∞[] h — h+ 1 i— i+ 1 Step2．3．3 如果 hK ，转 Step2．3．2；否则 ，J—J+ 1。如果J≤ N一2L+K一4转 ，Step2．3．1继 续 ；否则 ，终止。 St