基于MATLAB经验模态分解在地球物理资料中的应用.
摘 要
经验模态分解(EMD)是由Huang等人提出的一种新的分析非线性、非平稳信号的方法。本文研究经验模态分解原理及其在地球物理资料中的应用。首先研究经验模态分解的基本原理和算法,对地球物理资料(地震资料,重磁资料)进行EMD分解试验分析,然后研究基于EMD的Hilbert变换原理及其在提取地震属性信息中的应用,对实际地震时间剖面和时间切片进行EMD时频分析试验。
本文的方法研究和数据试验分析表明:经EMD分解变换得到的IMF序列是直接从原始时序数据中分离出来的,事先无需确定分解阶次,能更好反映原始数据固有的物理特性,每阶IMF序列都代表了某种特定意义的频带信息;EMD分解获得的IMF序列具有稳态性,对IMF进行Hilbert变换,就可以得到单个固有模态函数的瞬时振幅、瞬时相位和瞬时频率,这些信息可以清楚的显示信号的时频特征;EMD分析方法用于分解地球物理资料和作时频分析是有效的。
关键词:经验模态分解;地球物理;Hilbert变换;固有模态函数;时频分析
目 录
第一章 绪 论 1
第二章 经验模态分解(EMD) 2
第一节 EMD基本原理 2
第二节 EMD算法流程 3
第三节 EMD示例 5
第三章 EMD分解地球物理信号 8
第一节 EMD分解位场资料 8
一、EMD分解一维重力异常数据 8
二、EMD分解二维重力异常数据 9
第二节 EMD分解地震资料 13
一、EMD分解一维地震数据 13
二、EMD分解地震时间剖面 14
三、EMD分解地震时间切片 17
第三章 基于EMD的Hilbert变换与地震信号属性信息 19
第一节 基于EMD的Hilbert变换原理 19
第二节 提取地震信号属性信息 20
一、提取地震时间剖面属性信息 20
二、提取地震时间切片属性信息 22
结 论 23
致 谢 24
参考文献 25
附录:经验模态分解MATLAB程序 26
参考文献
[1] 程佩青. 数字信号处理. 北京:清华大学出版社,1994.
[2] Huang, N.E., Shen, Z., Long, S.R., et. al.. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Royal Society London, 1998: 903-995.
[3] 熊学军, 郭炳火. EMD方法和Hilbert 谱分析法的应用与探讨. 黄渤海海洋, 2002, 20(2). 12-21.
[4] 武安绪, 吴培稚, 等. Hilbert-Huang 变换与地震信号的时频分析. 中国地震, 2005, 21(2): 207-215.
[5] Hassan H. Hassan, GEDCO. Empirical Mode Decomposition (EMD) of potential field data. SEG 72th Annual International Meeting Expanded Abstracts, 2005: 704-706.
[6] 胡劲松. 面向旋转械故障诊断的经验模态分解时频分析方法及实验研究. 博士学位论文. 杭州: 浙江大学, 2003.
[7] 朱金龙. 经验模态分解方法及在癫痫脑电波中的应用研究. 硕士学位论文. 南京: 南京邮电大学. 2006.
[8] 秦勃, 时鹏. 经验模式多尺度图象分解. 工程图学学报, 2002, 4: 73-78.
[9] 张志涌等编著. 精通Matlab6.5版. 北京: 北京航空航天大学出版社, 1998.
[10] 管志宁. 地磁场与磁力勘探. 北京:地质出版社,2005.
[11] 曾华霖. 重力场与重力勘探. 北京:地质出版社,2005.
[12] 罗孝宽, 郭绍雍. 应用地球物理教程—重力、磁法. 北京: 地质出版社, 1991.