Hilbert--Huang变换在地震信号瞬时参数提取中应用.pdfVIP

Hilbert-Huang 变换在地震信号瞬时参数提取 1 中的应用 刘财，皮红梅，王典 吉林大学地球探测科学与技术学院，吉林长春（130026 ） E-mail ：liucai@ 摘 要：Hilbert-Huang 变换是一种新的分析非平稳非线性信号的方法，它先将信号分解为 有限个固有模态函数之和，再对固有模态函数作 Hilbert 变换求取时频谱，求取的时频谱在 时域和频域都具有很高的分辨率。现有的地震信号瞬时参数提取方法未考虑瞬时参数的物理 意义，本文将 Hilbert-Huang 变换应用到地震瞬时参数的提取上，使得计算得到的瞬时参数 有物理意义。实例表明对地震剖面做EMD 可以得到不同时间尺度上的特征，Hilbert 谱比传 统的时频谱在时间和频率域上的分辨率都要高，强反射层在第一阶 IMF 瞬时频率剖面上比 原瞬时频率剖面上表现得更为明显。 关键词：Hilbert-Huang 变换，经验模态分解，地震信号，瞬时频率 1 引言 在地震勘探中，地震信号瞬时参数的提取具有十分重要的意义。地震瞬时参数（主要指 瞬时振幅、瞬时相位和瞬时频率三参数）不仅可以直接用来研究地球介质的构造、岩性等方 面，还可以反演介质的品质因子等参数。计算信号瞬时参数有多种方法，常用的是 1946 年 Gabor 提出的解析信号法，即通过Hilbert 变换构造复信号求取瞬时参数的方法。1979 年 Taner 等人首次应用解析信号法构造复地震道求得地震瞬时参数，并探讨了瞬时参数剖面在地震解 释中的应用。Taner 等人的工作引发了人们研究地震属性的热情。随后，Robertson 等人将复 地震道方法应用于薄层解释，Barnes 的工作进一步推动和扩展了地震信号瞬时参数的理论基 础、计算方法以及应用范围。国内方面高静怀(1997)提出了利用解析小波求地震信号瞬时参 数的方法，虽然该方法对克服 Hilbert 变换对噪声干扰比较敏感等问题有一定效果，但应用 中仍有其局限性。事实上，通过 Hilbert 变换求取的信号瞬时参数并非对任何信号都有物理 意义。这种方法通常要求被分析信号是窄带或平稳的[Huang,1998]，而且实际信号通常含有 噪声，Hilbert 变换对噪声很敏感，如实际地震信号既非平稳又含有噪声，因此有必要寻找 新的途径来提取地震信号瞬时参数。Huang 于 1998 年提出了一种新的分析非平稳非线性信 号的方法－Hilbert-Huang 变换，该方法首先将信号分解为有限个固有模态函数之和，并认 为这些固有模态函数均为窄带信号，利用 Hilbert 变换对其求得的瞬时参数具有明确的物理 意义。本文将探讨 Hilbert-Huang 变换在提取地震信号瞬时参数中的应用。 2 提取信号瞬时参数的方法 提取信号瞬时参数有多种方法，如 Teager-Kaiser 算法、Shekel 算法、解析信号法等。 Vakman （1996）比较了各种提取信号瞬时参数的方法，认为解析信号法是最为合理的。李 世雄(2000,2003)讨论了解析信号法的适用范围。事实上，只有对窄带信号且在噪声干扰比较 小的情况下解析信号方法提取的信号瞬时参数才有物理意义。解析信号法是通过 Hilbert 变 换求取信号对应的解析信号，构成复信号从而提取出信号的瞬时参数。 1本课题得到高等学校博士学科点专项科研基金(20030183056)的资助。 -1- 设实信号为 x(t) ，用Hilbert 变换求出与其对应的虚部 y(t) [Taner,1979]: 1 ∞ ( ) x τ y(t) [ ( )] HT x t ∫ dτ