Hilbert-Huang变换在油气检测和工程物探中的应用.pdfVIP

石油天然气学报 (江汉石油学院学报) 2010年 6月 第32卷 第3期 JournalofOilandGasTechnology (J．JPI) Jun．2010 Vo1．32 No．3 Hilbert—Huang变换在油气检测 和工程物探 中的应用 何 峰 (中原油田分公司勘探开发科学研究院，河南NB8457001) 代金芝 (油气资源与勘探技术教育部重点实验室 (长江大学)，湖北NJJ1＼l434023) 关B 旭 (中原油田勘探开发科学研究院，河南濮阳457001) [摘要]希尔伯特 黄变换一种是针对非线性、非平稳信号的时频分析方法，该方法将信号分解为有限数 量的固有模态 函数 ，能清晰刻画信号的时频能量分布。研 究将此方法应用于地震信号和浅层工程物探信 号，分别做希尔伯特一黄燹换 ，从而得到振幅一频率一时间的三维谱分布 ，进而提取信号的频率一能量切片和 频率一相位切片，验证 了希尔伯特一黄变换 的时频域刻画能力，并得到了相应 的结论 。 [关键词]油气检测 ；希尔伯特一黄变换 ；固有模态函数；经验模态分解 ；希尔伯特谱 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752(2010)03—0243—05 美籍华人 N．E．Huang等l1于 1998年提出的 Hilbert—Huang变换 ，该方法主要包含经验模态分解 法 (EmpiricalModeDecomposition，简称 EMD)和 Hilbert谱分析 (HilbertSpectralAnalysis，简称 HSA)2个步骤 。EMD把数据分解成有限个 固有模态 函数 (IntrinsicModeFunction，简称 IMF)l2]； 对每个 IMF分量作 Hilbert变换 ，从而得到瞬时频率和能量 。该方法是一种新的具有 自适应性 的时频分 析方法[，可根据信号的局部时变特征进行 自适应的分解 ，消除了人为因素的干扰 ，克服了传统方法中 的谐波分量来表示非平稳 、非线性信号的缺陷 ，具有 良好 的时频聚集性 ，并得到 了极高 的时频分辨率 ， 非常适合对非平稳 、非线性信号进行分析 。 1 Hilbert—Huang变换原理 Hilbert—Huang变换假设信号是 由若干固有模态信号或固有模态函数组成，如果 固有模态信号之 间 相互重叠 ，便形成复合信号 。 1．1 经验模态分解 经验模态分解步骤如下 ：① 找出信号 ()中的局部极大值和局部极小值 ，利用 3次样条插值 ，分别得 到极大值包络和极小值包络 ；② 对每个时刻的局部极大值和极小值取平均 ，得到瞬时平均值 ；③ 用原始数 列 z(f)减去瞬时平均值 ，得到新数列 (￡)。 若 m(t)满足：① 局部极大值 以及局部极小值的数 目之和必须与零交点的数 目相等或是最多相差 1 个 ；② 在任何时间点，由局部极大值点和局部极小值点构成的 2条包络线平均值为0[。～，则 ”z(￡)为同有模 态 函数 ；若不满足，将 ()作为原始数列重复上述 3个步骤 ，直到满足2个条件为止 ；经过EMD处理 ，可从 原始数列中分离出 ”个固有模态函数和一个趋势项 ，如果把分离出来 固有模态函数 和趋势项加起来 ，则 得到原始数列，即： (￡)一∑IMF(￡)+ () (1) [收稿日期]2010 01 16 [作者简介]何峰 (1964一)，女，l996年大学毕业 ，工程师，现从事石油地质勘探综合研究工作 。 石油天然气学报 (江汉石油学院学报) 2O1o年 6月 1．2 Hilbert变换与 Hilbert谱 固有模态函数 (IMF)通常为窄带信号，满足 Hilbert变换 的要求。对每个 IMF做 Hilbert变换 ： H[IMF，()『]一IMF )一 一 dr (2) 7lJ … ￡ -／- 由此可构造…解析信号： Z (f)一 IMF ()+jIMF (f)一 A ()ejo,“