基于Hilbert-Huang变换的地震瞬时属性提取方法及应用.pdfVIP

2011 年6 月 APPLIED GEOPHYSICS （应用地球物理） 第8 卷，第2 期 基于Hilbert-Huang 变换的地震瞬时属性提取方法及应用//Seismic attribute extraction based on HHT and its application in a marine carbonate area, Huang Ya-Ping, Geng Jian-Hua, Zhong Guang-Fa, Guo Tong-Lou, Pu Yong, Ding Kong-Yun, and Ma Ji-Qiang. APPLIED GEOPHYSICS, 2011, 8(2), P. 125 − 133. DOI:10.1007/s11770-010-0279-z 基于Hilbert-Huang 变换的地震瞬时属性提取方法及应用 1 1 1 2 2 1 1 黄亚平 ，耿建华 ，钟广法 ，郭彤楼 ，蒲勇 ，丁孔芸 ，麻纪强 （1. 同济大学 海洋地质国家重点实验室 海洋与地球科学学院，上海 200092 ；2. 中国石化勘探南方 分公司，成都 610041 ） 摘要：希尔伯特—黄变换（Hilbert-Huang transform ，HHT ）是一种新的适合非平稳和非线性信号的分析方法， 由于地震信号一般呈现出非平稳与非线性特性，因此HHT 非常适合地震信号的分析。本文首先介绍了HHT 中 关于经验模态分解(Empirical Mode Decomposition, EMD)的实现过程，在此基础上分析了几种基于EMD 获得本 征模态函数（Intrinsic Mode Functions ，IMF ）来计算瞬时频率的算法，其中利用了两个采样间隔瞬时频率的平 均来计算瞬时频率，较好地反映了地震信号频率成分随时间变化的特征。将该方法应用于四川东北部某地区海 相碳酸盐岩地层三维地震叠后偏移数据处理，提取“三瞬”地震属性，与传统的希尔伯特变换提取的“三瞬”地 震属性进行对比，结果表明基于HHT 的“三瞬”地震属性结果具有更高的分辨率，IMF2 的瞬时相位能够较好 地刻画台地边缘生物礁相，IMF2 的瞬时频率亦具有较好的分带性。将IMF2 的“三瞬”地震属性与钻井等资料 结合分析，能够更好地识别沉积相的分布。 关键词：希尔伯特—黄变换，经验模态分解，瞬时频率，地震属性 的分析非线性、非平稳信号的方法，称为希尔伯 引言 特—黄变换（Huang 等, 1998; Magrin-Chagnolleau 和Baranuik, 1999; Hassan, 2005 ）。它将信号分解 地震信号是非平稳信号，随时间而变化，因此 为有限个具有固有模态的函数之和，并认为这些 时频分析是分析这类信号的有效手段 （Huang 等, 固有模态函数均为窄带信号。国内外一些学者都 1998; Taner 等, 1979; 高静怀等, 1997; 李春峰和 对 HHT 在地震信号处理中的应用进行了研究， Liner, 2005 ），但大部分方法都是依赖于傅里叶变 Battista 等（2007 ）将EMD 作为一种滤波器，用 换，因此同样受到傅里叶变换局限性的影响。短 于消除地震数据中拖揽的影响；周挚等（2008 ） 时傅立叶变换通常假设信号在窗函数的有效持续 基于 HHT 提取了昆明、下关重力固体潮的地震 时间内是平稳的，但此条件通常无法满足或近似 前兆信息。Wen 等（2009 ）将经验模态分解与关 满足，因而它的时间和频率分辨率都很低，不能 联维相结合进行储层预测，取得了较好的效果； 满足分析精度上的要求；小波变换是近十几年才