高分辨率地震勘探综述.docx
《高分辨率地震勘探综述.docx》由会员分享，可在线阅读，更多相关《高分辨率地震勘探综述.docx（24页珍藏版）》请在冰豆网上搜索。
高分辨率地震勘探综述
高分辨率地震勘探综述
摘要高分辨率是地震勘探的一个重要研究方向，涉及地震数据采集、处理与解释等各个方面。
在回顾高分辨率地震勘探发展历程及存在问题的基础上，重点阐述了高分辨率的评价机制，并对近年来发展的高分辨率方法原理及应用实例进行了详细介绍。
高分辨率是一个系统工程，实际生产中的各个环节都有可能对分辨率造成影响，因此，高分辨率不仅仅局限于某个单独的技术，需要同时发展采集、处理与解释各方面的技术，尤其是借鉴交叉学科的新方法。
关键词：
采集；处理；解释；高分辨率；评价机制
1概述
1.1高分辨率勘探的目的及技术发展历程
地震勘探是一种应用地震波在地下介质中的传播来对地下地质构造与岩性进行测量的技术，经过近一个世纪的发展，该方法已经成为最有成效的油气勘探物探方法。
纵观地震勘探的发展历程，高分辨率一直是科研、生产的重点与难点。
诚然，高分辨率地震勘探是一个系统工程，从地震资料采集、处理到解释，每一个环节都对分辨率有着重要的影响。
虽然采集、处理与解释分属不同的环节，考量高分辨率的角度也有所不同，但三者是有机联系的。
首先，野外地震数据的采集质量直接关系着地震勘探的成败，只有在采集质量得到保证的前提下，处理技术（诸如静校正、拓频与压噪技术等）才有发挥的空间，而地震处理得到的剖面又是解释的基础，解释成果则是高分辨率地震勘探的最终目标，三者环环相扣，紧密联系；其次，采集、处理与解释的方法也是相互影响与促进的，例如，采集观测方式的改变有可能对处理方法或参数提出新的要求（如可控震源采集对处理提出了谐波压制的要求等），解释方法的突破也有可能对处理提出新的标准（如AVO解释技术要求处理方法具有高保真度等）。
在阐述高分辨率地震勘探之前，有必要先介绍一下分辨率的概念及主要影响因素。
地震勘探分辨率是基于地震测量技术对地下构造进行空间测量的精度描述，在反射波地震勘探中可以概括如下：
可分辨的最小地质体的厚度或最窄地质体的宽度，前者称为垂（纵）向分辨率，后者称为横向分辨率[1-2]。
马在田[3]（2004年）又提出了偏移成像广义空间分辨率的概念，扩展了分辨率的范畴。
对于垂向分辨率，不少学者从不同的角度进行了研究，形成了几个具有代表性的分辨率标准，如Rayleigh准则[4]（分辨率极限为
）、Ricker准则（分辨率极限为
）与Widess准则[5]（小于
），其中，
是目前生产中最常用的分辨率判断标准。
上述准则都是以地震子波波长
作为定量估算垂向分辨率的参量，这表明，在地震波速度一定的情况下，地震子波的优势频带越高，波长（主周期）越短，则垂向分辨率越高。
对于横向分辨率，人们一般用菲涅尔带的半径进行描述[6]。
横向分辨率主要取决于勘探深度、地层速度与地震波频率，速度越小、深度越浅、频率越高，则横向分辨率越高，反之亦然。
广义空间分辨率除了与频率、速度有关，还受深度、距成像点地面投影的水平距离及炮检距等因素的影响[7]。
高分辨率地震数据采集技术主要体现在激发与接收两个方面。
首先，激发效果的好坏，直接决定地震数据品质的高低。
激发震源所产生的地震子波需要满足四个条件：
1.具有较宽的频带，包含足够多的高频与低频成分；2.具有较高的主频，能够满足纵、横向的分辨率要求；3.具有足够强的能量，保证深层反射波能被地面检波器接收；4.具有较好的信噪比水平，尤其是高频成分的信噪比，为处理阶段改善信噪比提供必要条件。
在科研与生产中，人们对激发因素的研究重点集中在埋藏井深与药量（以炸药震源为例）等方面[8-13]。
其次，接收是另一个影响地震数据采集质量的关键点，包括检波器类型及组合方式、面元大小、炮检距、道距、覆盖次数等因素[14-16]。
在具体生产中，需要对工区进行反复试验，获得最佳激发接收参数，尽可能提高采集数据的信噪比、拓宽有效波频带范围以及提高地震子波的主频等。
高分辨率地震资料处理技术是在数据有效采集的基础上拓宽频带、提高主频，其本质是对弱有效信号（一般指高低频成分）进行真振幅恢复。
常用的方法主要分三大类：
反褶积、吸收补偿及基于时频谱的拓频技术。
反褶积的基础是维纳滤波与褶积模型，实质是从输出信号中提取或恢复原始信号，通过压缩地震子波，达到展开频谱、提高资料分辨率的目的。
Robinson等[17]（1967）建立了预测反褶积的理论基础，有效地压制了海上鸣震与多次波。
UlrychT.J.[18]（1971）提出了同态反褶积，规避了地震子波最小相位及反射系数白噪的假设，可以同时提取地震子波与反射系数。
Wiggins[19]（1978）在反褶积中引入了熵的概念，提出了最小熵反褶积（MED），可以有效加强尖脉冲的特性，增强剖面“亮点”的识别能力。
20世纪80年代出现了最大似然反褶积[20]，该方法有效地解决了传统反褶积无法进行信噪分离的缺点，实际资料测试表明，该方法可以求出高分辨与高信噪比的反射系数。
随着盲系统理论的成熟，HaykinS[21]（1994）详细探讨了盲反褶积。
地震反褶积基本是一个盲过程，因为地震子波通常是未知的，而且检波器接收到的地震子波是重叠的，常规反褶积假设很多理想条件，因此在无任何假设条件下完成地震反褶积（即盲反褶积），具有较好的实用价值。
20世纪90年代后，随着小波分析的深入研究与应用，章珂等[22]（1999）基于二进小波变换提出了一种新的多分辨率地震信号反褶积，在地震信号二进小波变换域中的各尺度上分别进行其分辨率随小波尺度变化的反褶积，利用不同分辨率反褶积结果之间的相关性，以及测量噪声随尺度的衰减特性，从低分辨率反褶积结果逼近高分辨率反褶积结果。
GaryMargrave[23-25]（2002）提出了非稳态地震道褶积模型，并基于Gabor时频谱进行反褶积处理，突破了传统的时间域或频率域反褶积模式。
李芳等[26]（2013）认为径向道域数据比时间域数据更符合非稳态模型关于地震波垂直入射的假设，将非稳态反褶积拓展到径向道域实现，该技术处理的同相轴在横向上具有更好的连续性，对远炮检距与深层的地震资料恢复效果更好。
李国发等（2012）[27]考虑了实际反射系数序列的非白噪特点，提出了有色反褶积，将信号纯度谱作为反褶积输出的期望振幅谱，并进行有色补偿，进一步提高了地震资料的分辨率与反映薄层真实细节的能力。
吸收补偿是第二类高分辨率处理技术，以吸收与衰减为基础，具有物理意义。
由于地下介质的黏弹性，地震波在传播过程中能量不可避免地减弱，且高频成分衰减比低频严重，导致中深层地震资料的分辨率减小、成像精度降低。
补偿这种地层衰减常用的方法是反Q滤波，即在合理估算地下Q场分布的前提下，按照理论衰减模型（如Kjartansson衰减模型[28]），对各频率成分进行振幅补偿与相位校正。
在反Q滤波的具体实现方面，不同学者提出了不同的算法。
Hale[29]（1981）根据Futterman物理模型提出了用泰勒级数展开近似高频补偿的反Q滤波。
Hargreeaves[30]（1991）以一维方程为基础提出了与Stolt偏移相仿的针对常Q模型的相位反Q滤波。
Bano[31]（1996）将常Q值相位反Q滤波拓展到层Q值模型。
裴江云与何樵登[32]（1994）利用Fourier快速算法实现了基于波场延拓的反Q滤波处理，计算速度快、生产效率高。
王仰华[33]（2002）基于Gabor变换提出了稳定高效的反Q滤波，其后[34]（2006）对理论进行了完善改进，该算法目前在生产上颇受青睐。
姚振兴[35]（2003）根据非弹性介质中波的传播规律提出了在深度域剖面进行反Q滤波的方法，不仅考虑了介质吸收对地震波振幅的影响，还保证了所造成的波形畸变满足因果规律，即地震波具有某种频散性质。
陈增保等[36]（2014）将非稳态地震道褶积模型引入到反Q滤波中，用于消除地震记录的非稳态特性，取得了较好的处理效果。
经过几十年的发展，反Q滤波从理论到实践、从不稳定到稳定、从只补偿相位到振幅相位同时补偿、从低效到高效，逐步得到了完善。
另外，还有一类拓频技术，在高分辨率处理中起到了较为重要的作用。
该类拓频技术主要以时频分解为基本手段，在时间-频率（尺度）域进行高低频成分的恢复处理，达到压缩子波、拓宽频带的效果。
例如，汪小将等[37]（2009）将希尔伯特-黄方法（HHT）引入到地震数据的处理中，通过统计不同时间、不同频率的能量分布趋势，求取时频域的补偿因子，在保持地震数据相对振幅关系的同时，较好地提高了资料的分辨率。
尚新民等[38]（2014）利用谱拟合方法在地震数据的广义S变换谱上估算具有时变特征的地震子波谱，并通过消除地震子波谱的影响实现高分辨率处理。
周怀来等[39-40]（2014，2015）将非稳态褶积模型（Gabor域）推广到广义S变换时频域，通过估算并消除吸收衰减项、初始地震子波谱与相位校正量，实现地震资料的拓频处理。
张军华等[41]（2015）根据压缩感知理论用有限带宽地震数据恢复低频信息，实现低频成分的有效拓宽。
高分辨率地震资料解释技术是指在解释过程中为了获得更高分辨率而采取的各种技术手段，涵盖面较宽泛，这里主要介绍谱分解与谱反演技术。
谱分解（Spectraldecomposition）技术是Partyka[42]（1999）提出来的一项基于时频分析的地震属性分析技术，即对目的层段地震数据进行时频分解，并在频率域对数据进行重新成像，排除了时间域中不同频率成分的相互干扰，可以得到高于传统分辨率的解释成果。
进一步地，Partyka（2004）将谱分解技术与地震反演方法结合，提出谱反演的理论。
Castagna等[43-46]（2005-2008）对谱反演进行了系统的研究，丰富并完善了谱反演理论。
该理论认为，任何一个反射系数序列都可以分解成奇、偶分量，奇分量不利于检测薄层，而少量的偶分量可以提高薄层的分辨能力，其实质就是利用偶分量在厚度趋于零时的有效干涉提高地震资料的分辨率，该技术在薄层识别应用上取得了较好的效果。
1.2地震分辨率的判别准则及评价机制
1.2.1判别准则
这里简单介绍垂向与横向分辨率的判别准则。
（1）垂向分辨率
对于垂向分辨率，一般有三个准则（假设地震子波为零相位）：
Rayleigh准则、Ricker准则与Widess准则。
Rayleigh准则：
地震反射波的分辨率极限是
。
当两个子波的到达时间差大于或等于子波的半个视周期时，两个子波可以分辨，否则，不可分辨。
半个视周期指子波主极值与相邻异号次极值的时间间隔。
Rayleigh准则适用于主极值幅度显著大于次极值幅度的零相位子波。
Ricker准则：
地震反射波的分辨率极限是
。
当两个子波的到达时间差大于或等于子波主极值两侧的两个最大陡度点之间的时间间隔时，两个子波可分辨，否则，不可分辨。
对比发现，Ricker极限为子波导数两个异号极值点的间隔，Rayleigh极限为子波导数两相邻过零点的间隔。
Widess准则：
当两个极性相反的子波到达时间差小于1/4视周期时，复合波形接近子波的时间导数，不随时间差改变，两个异号极值间隔等于子波视周期的一半；复合波形的幅度与到达时间差近似呈正比关系，利用振幅可解释地层厚度。
显然，Widess准则适用于厚度小于
的地层。
（2）横向分辨率
横向分辨率判别准则按照未偏移剖面与偏移剖面可以分为2大类。
对于未偏移剖面，人们普遍用菲涅尔带（FresnelZone）来衡量横向分辨率。
菲涅尔带的定义：
反射波实际上是由反射面上相当大的一个面积内返回的能量叠加而成的，产生相干干涉的反射区域称为菲涅尔带。
菲涅尔带半径：
（1）
其中，
为平均速度，
为反射时间，
为地震子波主频。
需要注意的是：
上式是一个近似的公式，其假设前提条件为
，当深度
较小（与波长
接近）时，该公式是不适用的。
对于偏移后的地震剖面就不能再照搬上面的菲涅尔带准则了。
很多人误认为偏移后的剖面好像是将震源放在了反射面上，并且接收点与震源重合，此时反射界面深度为零，菲涅尔带的大小也为零，这种理解是错误的，错误地将该前提条件忽视了，当
时，根据精确的公式，可推导出
，在理想情况下这是地震剖面经过偏移处理后地震横向分辨率可达到的极限值。
对于偏移后的剖面，其剖面的横向分辨率可以根据Chen与Schuster[47]（1999）提出的公式进行定量的计算。
已知
为地层深度、
为波数、
、
分别为
、
方向的偏移孔径的一半。
（2）
另外，还可以根据Lindsey提出的理论估算偏移后的横向分辨率：
菲涅尔带半径大致变为偏移前的十分之一。
1.2.2评价机制
地震资料垂向分辨率高低的评价可以从两个角度进行：
子波频谱与波形形态。
（1）子波频谱（绝对频宽与相对频宽）
绝对频宽是最高频率与最低频率之差，相对频宽是两者之比。
绝对频宽决定子波的包络，；相对频宽决定子波的相位数。
通过零相位带通子波的频移与缩放可以得到两个重要认识：
当绝对宽度不变时，子波包络基本不变，主频越高，振荡越多；
当相对宽度不变时，子波形态基本不变，主频越高，波形越窄。
频宽与分辨率的关系，已有相当经典的研究，读者可以参阅相关的文献，这里不作赘述。
需要说明的是，在构造解释中，强调高频成分的恢复，主要看绝对频宽；在反演过程中，强调低频信息的恢复，主要看相对频宽。
因此，从频谱带宽的角度，高分辨率可概括如下：
拓宽高频（主要影响绝对频宽），提高构造解释的分辨率；拓宽低频（主要影响相对频宽），提高反演结果的分辨率。
（2）波形形态（波长与旁瓣）
传统的分辨率判别准则主要以子波波长为标准，这容易让人产生误解，即高分辨率地震子波的要素就是波长。
图1展示了不同带宽的理论子波波形（上）与振幅谱（下），频带宽度从左往右逐渐增加，分别为40Hz，50Hz，60Hz与80Hz。
随着频宽的增加，子波长度明显变短，但同时，还有一个特征发生了变化，即旁瓣大小（见红色基准线）。
因此，我们把零相位子波形态简化为两个要素来定量描述，即波长与主旁瓣极值比。
下面重点考察子波形态与地震垂向分辨率的关系。
图2a是两个理论子波波形，两者波长不同、主旁瓣极值比相同，图2b是仿真反射系数序列，图2c是子波频谱，图2d是合成地震记录（蓝色记录由蓝色子波合成，红色记录由红色子波合成，下同）。
对比反射系数序列与两个合成记录发现，蓝色记录在反射系数较为密集的地方（见红色框部分）分辨率明显低于红色记录，在不存在反射界面的地方（见黑色框部分）却产生了较强的波峰（谷），这是由子波复合调谐引起的假象，与波长有关。
图3a是两个理论子波波形，两者波长相同、主旁瓣极值比不同，图3b是仿真反射系数序列，图3c是子波频谱，图3d是合成地震记录。
对比发现，蓝色记录在不存在反射界面的地方（见黑色框部分）产生了较强的波峰（谷），这是由子波旁瓣过强引起的假象，与主旁瓣极值比有关。
图1不同频宽的理论子波波形（上）及频谱（下）
（a）（b）
（c）（d）
图2a：
子波波形（第一组）；b：
反射系数序列；c：
子波频谱；d：
合成地震记录
（a）（b）
（c）（d）
图3a：
子波波形（第二组）；b：
反射系数序列；c：
子波频谱；d：
合成地震记录
图4是一组模型数据的对比，其中，图4a是根据波动方程正演的宽频带高分辨率剖面（10-160Hz），图4b是窄频带低分辨率剖面（10-80Hz）。
对比发现，高分辨率剖面中存在的同相轴（见红色圈部分），在低分辨率剖面中不明显；高分辨率剖面中不存在的同相轴（见蓝色圈部分），在低分辨率剖面中反而出现了。
这个现象恰好可以用图2与图3来解释，同时道出了高分辨率处理的真相：
压缩子波长度，同相轴数增加；压制旁瓣能量，同相轴数减少。
因此，高分辨率地震子波具有两个判别要素：
波长与旁瓣。
（a）（b）
图4a：
宽频带高分辨率剖面（10-160Hz）；b：
窄频带低分辨率剖面（10-80Hz）
1.3当前高分辨率勘探存在的问题
高分辨率勘探是一个系统工程，牵一发而动全身。
首先，在地震资料处理过程中，某些环节有可能对地震资料的分辨率造成影响。
例如，叠加处理，如果时差校正不准确，横向非一致性叠加，将损失地震数据的高频成分，降低分辨率。
图5是仿真数据测试，a为同相一致的叠前CRP道集，b为同相非一致的叠前CRP道集，各道波峰随机错开N毫秒（N≤2），c为叠加波形，d为叠加频谱，红线为一致叠加结果，蓝线为非一致叠加结果。
数值计算发现，一致叠加波形相对长度为4.86ms，非一致叠加波形为5.79ms，拉伸了19.1%；一致叠加谱频带宽度为68.77Hz，非一致叠加谱带宽为58.02Hz，降低了15.6%。
对比表明，非一致叠加改变了叠加记录的波形，损失了相当程度的分辨率（尤其是高频成分）。
因此，需要采用各种技术改善叠前道集，如高精度静校正、精细速度分析、子波整形与剩余时差校正等技术，尽可能实现同相一致叠加，保护高低频成分。
同理，在其他处理过程中，应注意有可能造成高低频损失的因素。
其次，信噪比水平是制约高分辨率的一大因素。
在实际生产中，噪音是不可避免的，其对资料的品质与可信度造成非常大的影响。
地震记录的高频成分能量相对较弱，信噪比较低，如果没有合理有效的噪音滤除手段，即使在采集或处理过程中将高频信息恢复，也将因为信噪比过低而不得不放弃这些信息，高分辨率就不能很好地实现。
因此，噪音压制技术是实现高分辨率的一个很重要的前提。
第三，保真度是影响高分辨率的又一大因素。
随着岩性等非常规圈闭逐渐替代构造圈闭，地震资料的保真度越来越受到重视。
因此，地震资料处理过程中的保幅性是需要特别注意的一个方面，如去噪、偏移等技术，甚至连高分辨率处理本身也要考虑保幅性的问题。
（a）（b）
（c）（d）
图5模型分析（a：
同相一致的CRP道集；
b：
同相非一致的CRP道集；c：
叠加记录；d：
叠加振幅谱）
最后，就高分辨率处理技术而言，目前还存在一定的问题。
反褶积技术对地震模型存在较多的假设（如地震子波最小相位且非时空变、反射系数白噪、噪音干扰较小等），后来发展的反褶积技术都是尽可能避免或弱化某些假设，但始终不能精确地提取地震子波，处理效果受到不同程度的压制，不可避免地存在误差。
反Q滤波存在三个关键问题：
吸收衰减理论、稳定高效补偿与地层Q值估算。
反Q滤波的基础是吸收衰减理论，目前使用的模型是基于一定的假设推导的，实际地震资料与理想模型存在一定的差异，当资料偏离模型较大时，反Q滤波存在适用性的问题。
在补偿过程中，如何进行稳定高效的补偿，人们一直在尝试与探索新的算法。
地层Q值的估算一直是个难题，而Q值又是反Q滤波的基础，其准确性将直接影响补偿效果的好坏。
基于时频分解的拓频技术存在的问题与反褶积类似，即拓频因子及相位校正量的准确提取。
2高分辨率勘探技术的应用效果
这里，摘录一些（不包括全部）高分辨率勘探技术的应用效果，分别作展示与说明。
2.1宽频地震勘探技术
宽频地震勘探技术是实现高精度勘探的重要方法之一[48]，涵盖了采集、处理与反演等各个方面，低频与高频端的拓宽显著提高了地震资料品质，改进了对深部地质体的穿透力与照明度，能够实现深部目标体的清晰成像，获得更多的地层结构及细节等信息，提供更稳定的反演结果，提高地震资料的解释精度。
目前，国外先进的宽频技术采用单点激发、单点接收、室内组合处理的方式，形成了采集—处理—解释一体化的宽频地震勘探技术方案。
要获得宽频的地震信息，在震源激发时尽可能产生较宽的频谱，这需要对可控震源进行设计，定制扫描激发信号。
在设计信号时，并不是盲目增加信号频宽、任意选取扫描信号，宽频信号受到激发条件、接收条件、大地介质条件与处理方法等因素的影响，只有针对工区地形地质条件与目的层要求，综合分析各要素的影响才能制定出最优的扫描信号，达到预期结果。
同时，选择与宽频信号相匹配的采集设计与处理方法才能最大限度地发挥宽频信号的优势，改善地震资料的品质，提高勘探精度[49]。
2.2基于信号有色纯度谱的反褶积
线性反褶积技术的核心是通过展开子波频谱达到压缩子波的目的。
其中涉及两个主要方面：
1.如何从地震记录中估算子波；2.如何确定期望输出子波。
众所周知，地震资料的分辨率与信噪比是两个相互矛盾的指标。
由于反褶积本身并不改变资料的信噪比谱，高频成分的信噪比较低，高频能量的抬升将降低资料整体的信噪比水平。
为了兼顾分辨率与信噪比，李国发等（2012）提出了基于信号纯度谱的有色反褶积技术，重点考虑了期望输出的问题。
该技术以信号纯度谱（由信噪比谱定义）作为反褶积的期望输出，在输出的子波频谱中，高信噪比频段得到增强、低信噪比频段得到削弱，具有自适应的特点。
同时，对输出频谱进行有色补偿，进一步改善了地震数据反映薄层真实细节的能力。
图7是东部某油田地震数据的处理实例，a为原始叠加剖面，b为基于信号纯度谱的反褶积处理剖面，c为在b的基础上进行有色补偿的剖面，红线为根据测井信息合成的地震记录。
对比发现，b剖面整体的分辨率比原始剖面a高，且信噪比得到较好的兼顾；在2200ms附近（蓝色线圈内），a剖面为一根强轴，b剖面分离为一强一弱两个反射，c剖面则分离为两个能量相近的同相轴，在视觉分辨率上b与c数剖面与合成记录一致，不同的是，c在反射特征上与合成记录更匹配，更适用于薄储层的预测。
（a）（b）
（c）
图7东部某油田实际资料处理（a：
原始剖面：
b：
基于信号纯度谱的反褶积处理剖面；c：
基于信号纯度谱的有色反褶积处理剖面）
（摘自《基于信号有色纯度谱的反褶积》，李国发等[27]，2012）
2.3Gabor反褶积
褶积模型是反褶积的基础，模型可靠与否，直接影响反褶积的效果好坏。
传统的褶积模型认为地震子波是时不变的，即不考虑大地滤波对子波的影响。
为了适应子波时变的情况，Margrave（1998）提出非稳态滤波方法，随后（2004，2011）将传统的稳态地震道模型推广为非稳态地震道模型，建立了基于Gabor变换的非稳态褶积模型。
非稳态反褶积具有常规反褶积的优点，并且有效补偿地层吸收引起的能量衰减。
近些年，不少学者对Gabor反褶积进行了研究，取得较好的应用实效。
为了更好地满足自激自收的假设条件，李芳等（2013）将地震数据从XT域（时空域）变换到RT域（径向道域），基于广义自激自收模型提出了径向道域变步长采样非稳态反褶积技术，在模型测试与资料处理中取得了优于常规方法的效果。
图8是某工区单炮道集的处理实例，a为原始单炮道集，b为常规XT域非稳态反褶积处理道集，c为变步长XT域非稳态反褶积处理道集，d为变步长RT域非稳态反褶积处理道集。
从处理效果中可以发现，变步长采样可以更有效地恢复中深层地震数据，从XT域变换到RT域能有效恢复大炮检距数据，增强横向连续性。
（a）（b）
（c）（d）
图8某工区单炮处理结果（a：
原始单炮道集；b：
常规XT域非稳态反褶积处理道集；c：
XT域变步长采样非稳态反褶积处理道集；d：
RT域变步长采样非稳态反褶积处理道集）
（摘自《径向道域变步长采样叠前非稳态反褶积处理方法研究》，李芳等[26]，2013）
2.4基于HHT的拓频技术
在基于时频谱的拓频技术中，时频分解的精度对处理效果有着重要的影响。
希尔伯特-黄（HHT）是由Huang等[50]（1998）提出的，具有非常高的时频分解能力，其核心思想是将时间序列的信号通过经验模态分解算法（EMD）分解为多个固有模态函数（IMF），利用Hilbert变换构造各分量的解析信号，计算出瞬时频率与振幅，即得到Hilbert谱。
由于HHT没有固定的先验基底，IMF是基于信号的时间特征得到，其分解过程是自适应的，特别适用于非稳态信号（如地震记录）的处理。
汪小将等（2009）将HHT引入高分辨率处理，基于HHT谱统计各频率能量的展布，通过估算合理的补偿因子对各频率的振幅进行有效补偿。
实际资料处理结果表明，该技术在提高地震分辨率的同时能较好地保持地震信号
升级会员 