分数阶希尔伯特变换在提高纵向分辨率中的应用.pdfVIP

第31卷第4期 华 北 地 震 科 学 Vot．31，No．4 2013年12月 NORTHCHINA SCIENCES Dec．，2013 EARTHQUAKE DOI：10．3969／j．issn．1003—1375．2013．04．004 分数阶希尔伯特变换在提高纵向分辨率中的应用 宋洪勇1，别立珍2 (1．山东省煤田地质局第四勘探队，山东潍坊261206；2．山东省煤田地质规划勘察研究院，山东泰安271021) 摘要：提出一种广义希尔伯特变换——分数阶希尔伯特变换(FHT)，其可以进行任意角度的相位转换，并且可以提 高地震勘探的纵向分辨率。地震勘探中，受纵向分辨率的影响，同相轴往往是复合叠加的结果，在地震剖面上，层 位和测井曲线无法对应。分数阶希尔伯特变换通过任意角度移相，克服了常规移探相90。的缺陷——相位为零，解 决了相位任意角度移相问题，使同相轴的极性赋有岩性信息，同相轴的振幅具有了岩性地层的概念。从而能够在 地震剖面上使地层和测井曲线相对应。 关键词：分数阶希尔伯特变换；相位；地层；正演模拟 中图分类号：P631．443文献标志码：A 文章编号：1003—1375(2013)04一0013—09 得到是阶希尔伯特算子H。一H‘。对志阶希尔伯特变 0 引言 换有： 近年来，伴随着科学的发展和现代信号处理的 (3) 向前发展，人们从不同的研究领域和应用角度出发， 如果将上式中的是换成实数y∈(一。。，。。)得到的 提出了拓展经典希尔伯特变换(以下简称HT)，也算子H。=H”对应的传递(乘子)函数为： 即将其由整数阶推广到分数阶，这样得到的分数阶 左。(U)一 希尔伯特变换，即任意角度希尔伯特移相变换，其在 f exp(--西丌／2) U＆01 地震数据处理及解释中有着重要的应用口-3]。 U＆=0一 [-exp(--ivlr／2)+exp(iv，r／2)]／2 l J exp(iv，r／2) U＆0 1 方法原理 cos(埘／2)一isin(vx／2)·sgn(U6)，yE(一C,O，CO) 如果假设实信号s(￡)EL2(R)，其傅里叶变换 因此算子H。可分解成如下形式： v H。=cos(婀／2)·J+sin(圻／2)·HE(～ j(u＆)一F{s(z))一IS(z)P一扩dt，其中U表示信号 高 oo，oo) (4) 的模拟角频率，F为傅里叶变换算子，则s(￡)和 式中：J为恒等算子；实数v取正整数时，就得到通常 ；(U＆)是信号时间域和频率域的一对等价表 意义下的整数阶HT算子，并且算子H，关于v是连 示[4-6]。 续的，是整数阶希尔伯特变换算子的连续内插，称之 实信号s(z)的HT，即一阶HT定义为： 为广义HT算子，相应的