《经验模态分解法在地震资料处理中的应用2016年论文》.pdfVIP

第 l7卷 第 3期 重庆科技学院学报(自然科学版) 2015年 6月 经验模态分解法在地震资料处理中的应用 李世涛 张 勇 范存辉 梁 亮 (1．中海油能源发展股份有限公司上海采油技术服务分公司，上海 200030； 2．西南石油大学油气藏地质及开发工程 国家重点实验室，成都 610500) 摘 要：在连续小波基的平稳序列时频分析方法基础上，对非平稳时间序列的时频分析方法深入研究，引入经验模 态分解法，利用信号的局部特征 —— 时间尺度，从原信号中提取若干个固有模态函数(IMF)和一个残余量，用多个 固有模态函数表征信号的局部特征 ，用残余分量表征信号缓慢变化的趋势，并形成非平稳信号经验模态分解算法。 通过实例应用研究验证了EMD分解后所作的瞬时Hilbert幅值谱，比连续小波变换CWT分解后的谱的分辨率高，并 与地震道波形匹配 自然，可用于地震非平稳数据的时频分析。 关键词：非平稳序列；时频分析；经验模态分解；固有模态函数；地震资料 中图分类号：P631 文献标识码：A 文章编号：1673—1980(2015)03一O030—04 小波变换在石油勘探 中已被广泛应用，它本质 和最小值形成的包络均值为0。下面介绍非平稳信 上是一种 自适应加窗富氏变换，要求小波窗内的信 号经验模态分解算法。 号短时平稳，没有摆脱富氏变换的局限性；其次，小 (1)EMD分解步骤，分解过程采用以下 “筛分” 波基的多样性及零相位不适应非零相位子波生成的 算法 。 信号，即小波基的有限长度还会造成信号能量的泄 第 1步，对信号 (t)求取所有的局部极值点， 漏 ：因此，基于小波变换的时频分析很难适应非平稳 用样条函数分别拟合所有的局部极大值点和极小值 时间序列的时频分析。为此，本次研究应用 Huang 点，形成上下包络线b(t)和b(t)。 等人(1996)提出的经验模态分解法(EmpiricalMode 第2步，计算上下包络线的均值： Decomposition，简称 EMD)中的固有模态函数 (In． mI．1(t)=0．5(b ()+b(t)) (1) trinsicModeFunction，简称 IMF)作为基函数，对非平 并将原始序列 (t)减去包络均值 ： 稳时间序列的时频分析方法进行深入研究，把该研 hlll(t)= (t)一ml-l(t) (2) 究成果应用于地震和测井资料联合处理解释之中。 此时判别 h．．，()是否满足 IMF的2个条件 。如果 不满足，将 h¨(t)作为原始数据重复上述处理过 1 经验模态分解法时频分析 程 ，直至序列 经验模态分解是对一种非线性、非平稳信号的 hm (t)=h 一l()一m (t) (3) 新的时频分析方法 。它能够对数据进行线性化、平 满足 IMF的条件 。这样 ，就得到 IMF的第 一个 稳化处理，并在处理过程 中保留信号 自身特性；它利 分量 ： 用非平稳信号的统计特征尺度随时间信号变化的局 c。(t)=h．(t) (4) 部时变特征，从原信号中提取若干个固有模态函数 第3步，由信号 (t)分离出第一个分量C(t)， (IMF)和一个残余量，用多个固有模态函数表征信 得剩余量 ： 号的局部特征，用残余分量表征信号缓慢变化的趋 ，()= ()一c() (5) 势。固有模态函数需满足 2个条件 ：(1)整个数据