改进的F—X域EMD去噪技术及分布式并行实现.pdfVIP

石油天然气学报 2012年6月 第34卷 第6期 JournalofOilandGasTechnology Jun．2012 Vol-34 No-6 改进的F—x域EMD去噪技术及分布式并行实现 邴萍萍 ，曹思远 ，中国石油大学 (北京)油气资源与探测国家重点实验室、 胡言防，路交通 kCNPC物探重点实验室，北京102249 [摘要]在地震数据处理领域 ，较高的计算效率及较好 的处理效果一直是地球物理工作者追求 的 目标 。介 绍 了改进的 F—X域 (频率一空 间域)经验模态分解 (EmpiricalModeDecomposition，EMD)方法 ，首先将 地震信号转换到F—X域后分解成一系列 固有模态 函数 (IntrinsicModeFunctions，IMF)，然后通过小波闽 值滤波对 固有模态函数除第一个分量之外的其他分量进行滤波处理来达到去噪 目的。此外，为解决改进 后的FX域 EMD去噪方法计算效率问题 ，提 出了分布式并行算法，兼顾 了运算效率与计算精度 。数值模 拟结果表 明，与F-X域预测滤波技术相比，改进后的F—X域 EMD去噪方法可 以更加有效地衰减随机噪声 和压制线性干扰 ；最后将该方法应用于实 际地震数据处理 ，进一步验证 了该方法的有效性和优越性 。 [关键词]F-X域 EMD；IMF；小波 阈值去噪；分布式并行；随机噪声；线性干扰 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752 (2012)06—0061—04 经验模态分解_1](EmpiricalModeDecomposition，简称 EMD)是 由Huang等人在 1998年提出 的一种全新的自适应信号时频分析方法 ，它不需要设定基函数，可 自适应地将信号分解成若干个本征模 态函数分量 (IntrinsicModeFunction，简称 IMF)，是一种适用于分析非线性非平稳信号的方法 。EMD 滤波可以作为一个非常有效的二进制滤波器口]，利用 EMD进行噪声滤除的同时，不会损害有效信号。 在此基础上 ，F—X域 EMD去噪l4 技术应运而生 ，然而 ，F—X域 EMD去噪法实际上是对 EMD分解 出 的 IMF1分量进行了去除，即认为噪声都集中在 IMF1分量 ，而忽略了IMF2分量和 IMF3分量中也可 能存在振荡较快噪声的事实。其次 ，由于对地震信号做 EMD分解是一个不断筛选的过程 ，并针对单个 地震道展开 ，传统情况下，计算机 中的程序运行是基于串行 的，只有一小部分资源被利用 ，因此针对实 际地震资料的F—X域 EMD去噪将会非常耗时 。 ， 为解决上述问题，笔者对F-X域 EMD去噪技术进行如下改进 ：首先将小波阈值去噪算法应用于沿 X 方向EMD分解后的IMF2~IMFn分量中，尽可能消除混杂在信号中的噪声成分 ；其次，为了提高计算效 率 ，提出分布式并行算法，实现了改进后的F-X域 EMD分布式并行去噪。理论模型和实际地震数据处理 结果表明，笔者方法不仅可以有效地衰减随机噪声和压制线性干扰，而且能够明显地提高计算效率。 1 改进的F—x域 EMD去噪与分布式并行实现 1．1 改进的F-X域 EMD去噪方法 利用改进的F—X域 EMD方法对地震信号去噪的基本思想是 ：首先将地震信号由时间域变换到频率 域 ，然后沿空间 X方向对每一个频率值做 EMD分解 ，信号经 EMD分解后 ，得到的是 r／个本征模态函 数 IMF1～IMFn和一个残余项。对于 IMF1分量 ，直接去除；对于剩余 的IMF分量 ，利用小波变换把 含噪分量再次分解到尺度域 ，然后在尺度域对含噪声的小波系数进行 阈值处理 ，最后将处理后的信号进 行重构即可得去噪后的结果。改进 的F—X域 EMD去噪技术与 F—K滤波 、F—X预测滤波相比，其优点是 [收稿 日期]2012—03—20 [基金项 目]国家科技重大专项 (2011Z1)；国家自然科学基金项 目 。 [作者简介]邴萍萍 (1981一)，女，2004年大学毕业，博士生，现主要从事地质资源与地质工程研究工作。 石油天然气学报 石油物探与测井 2012年 6月 信号的非线性和非平稳性在滤波后得到了最大限度的保留。利用 F—X域