希尔伯特-黄变换在浅海沉积体系地震精细识别刻画中的应用

鲁银涛，曹晓初，冉伟民，栾锡武，许小勇，李东，杨涛涛，邵大力，魏新元

( 1. 中国石油杭州地质研究院，浙江 杭州 310023；2. 中国地质调查局青岛海洋地质研究所，山东 青岛 266071；3. 青岛海洋科学与技术试点国家实验室 海洋矿产资源评价与探测技术功能实验室，山东 青岛 266237；4. 中国海洋大学 海洋地球科学学院，山东 青岛 266100)

1 引言

浅海、海陆过渡相沉积体系一般发育复杂网状水下分支河道系统，这类水下分支河道往往规模小、数量大、期次多、切割关系复杂，常规的地震资料较难识别这类沉积体系的边界和空间展布特征。由于垂向、侧向迁移快，沉积结构单元的空间展布特征复杂，对有利储集体空间展布的刻画存在较大困难[1-4]。现今对这类沉积体系的研究主要基于地震资料，所以对沉积体系的空间展布特征的刻画很大程度地受限于地震资料的品质[5-7]。

希尔伯特-黄变换（HHT）是一种处理复杂的非线性、非平稳地震信号的工具[8-10]。针对复杂地震信号进行时频分析，HHT 能够有效提高地震资料分辨率、地震瞬时属性、消除噪声，对识别隐蔽型地质体、小型沉积体、小断裂等具有很好的效果[11-12]，特别是对检测地质体边缘具有较高的实用性。HHT 通过经验模态分解（EMD）方法将信号分解成有限个固有模态函数（IMF）分量之和，根据有效波和噪声分布在不同的分辨尺度范围内进行数据重构，提高地震资料品质；或者对IMF 分量进行希尔伯特变换求取其瞬时函数[12]，得到高品质地震属性，识别地质体。

本文针对西纳土纳盆地浅海相沉积中的小型分支河道开展识别研究，旨在刻画水下分支河道的平面展布特征。在常规三维地震资料的基础上，通过HHT对三维叠后地震资料进行重新处理，提高了地震分辨率，有助于在地震剖面上识别河道的边界；通过对HHT 变换后的三维地震提取瞬时地震属性，河道的平面展布特征在属性图上更加清晰。与常规地震属性对比，经过 HHT 的地震资料对沉积体系的精细刻画研究提供了更高分辨率的地震数据。

2 区域地质背景

研究区位于马来盆地与西纳土纳盆地交汇处（图1）。晚始新世至早渐新世，研究区由于受印澳联合板块与欧亚板块的相互挤压，形成SE 向走滑运动，导致该区发生陆内断陷作用[13-18]，发育一系列NW-SE 向地堑-半地堑，发育裂陷期沉积，沉积以湖相沉积为主，是烃源岩发育主要层系，研究区就位于该地堑-半地堑体系的东缘[19-21]；而后早渐新世-晚渐新世，构造作用很弱，发育后裂陷期沉积，以沉积河流相（砂岩为主）与湖相（泥页岩为主）交互沉积地层为主，是储层发育主要层系，期间发育的湖相泥岩是比较好的层间盖层[22-25]；晚渐新世-中中新世，南海扩张停止导致该区域由前期拉张应力场改变为挤压和压扭应力场环境[26-27]。NW-SE 向的挤压应力场重新激活了裂谷期的控盆断裂并使其发生反转，形成一系列NE-SW 方向的反转构造[28-30]，使沉积地层向上抬升，发育大型不整合界面，发育反转期沉积，主要由一系列海进海退旋回沉积组成，沉积环境从湖相到滨浅海相再到含煤层较多的海陆过渡相变化；中中新世后，构造相对稳定，发育后反转期沉积，主要发育滨浅海相-海陆过渡相沉积[25-28]（图1）。

图1 研究区构造位置及综合地层柱状图[13-20]Fig. 1 The location and comprehensive stratigraphic histogram of the study area [13-20]

本文研究目标主要集中在上新世-更新世的后反转期沉积层，该时期构造活动稳定，发育滨浅海相沉积，水下分支河道发育，多期叠置特征复杂，是刻画河道沉积结构单元的有利层段。

3 数据和方法

3.1 地震资料特征

本次研究使用的三维地震数据为叠后时间域地震数据，主频约30 Hz，主频较低，频带较窄（图2）。研究目的层上新统-全新统地层对应地震剖面上0.5 s以上层段，构造相对稳定，发育上新统-更新统海陆过渡相沉积，以水下分支河道为主（图3）。该层段地震剖面主要呈平行-亚平行强弱相间的反射特征，夹弱反射地震反射体（图4）。

图2 原始地震资料频谱（蓝线）与HHT 变换后第一分量地震频谱（红线）对比（振幅归一化处理）Fig. 2 Comparison between the spectrum of original seismic data (blue line) and the spectrum of the first component after HHT transformation (red line) （amplitude normalization）

图3 研究区三维原始叠后地震资料0.3 s 等时切片（A 为振幅切片，B 为相干切片）Fig. 3 0.3 s isochronous slice of 3D original post stack seismic data in the study area (A is amplitude slice and B is coherent slice)

这些弱反射地震反射体为下分支河道，在地震剖面上表现为地震反射特征杂乱、不连续，且振幅变弱。水下分支河道底界在地震剖面上有明显的河道下切侵蚀特征（图4），呈“U”型、“V”型或“W”型的强振幅反射；而水道顶面为水平状强振幅反射特征；水道内部为弱振幅-空白反射，频率较低，这种反射特征可能与充填物有关，岩性较均一的沉积物导致波阻抗差异较小，不易形成波阻抗界面，所以整体呈较均一反射，呈弱振幅甚至空白反射。基于原始叠后地震资料，尽管识别出河道大致形态，但受限于地震资料分辨率影响，河道与围岩的边界仍刻画不清，且多期河道发育交错切割关系及其内部精细充填特征仍不明确，需要进一步提高地震资料分辨率进行精细特征研究。

3.2 希尔伯特-黄变换

希尔伯特-黄变换时频分析方法是黄变换与希尔伯特变换的结合算法，黄变换可将原始信号自适应地分解成一系列本征模态函数分量（图5），对这些本征模态函数分量进行希尔伯特变换就可以获得有意义的瞬时频率，从而给出频率随时间变化的精确表达[3-12,31-33]。由于相比较小波变换提取的分量，黄变换分解的本征模态函数分量在频率成分方面更加纯净，是一种较优越的信号分解算法，黄变换被广泛运用在信号处理中[3-12,31-33]。结合希尔伯特变换后，HHT 具有很好的时频分辨能力，在地震属性剖面提取方面的应用有助于对原始剖面中的细微构造特征或岩性界面进行较好的分辨和验证，在含噪地震数据的噪声剔除方面能够进一步得到信噪比品质较好的地震剖面，有利于对地下地层的正确识别[3-7,11-12]。所以，HHT 变换现在被广泛地运用在叠后地震资料处理上。

图5 理论信号黄分解本征模态分量图Fig. 5 Eigenmode component diagram of theoretical signal Huang transform

4 HHT 效果

经HHT 后，信号频率有了明显增高，主频提高至50 Hz 左右，且频带宽度明显增加（图2）。经HHT的IMF1 分量具有高频优势，分辨率明显增高，有助于表征薄层信息，特别能够反映岩性变化的细节特征。在砂泥岩互层的浅海、海陆过渡相沉积体系下，HHT 后的地震数据能够检测薄层砂泥岩岩性变化的特征，能有效检测沉积体的边缘，区分水下分支河道和围岩的边界和砂泥岩沉积界面，为识别多期河道、刻画河道横向迁移、垂向叠置等特征提供了地震资料基础。

由HHT 第一分量即高频分量（IMF1）结果可见（图6），虽然由于原始叠后地震资料弱振幅和空白反射范围较大，导致了经HHT 后的第一分量也存在局部的空白反射和弱振幅反射带，但是整体来看，经HHT后的IMF1 分量整体频率有了明显的提高，直观上表现为地震反射同相轴的数量有了明显的增加。比较原始叠后地震资料可见，在原始叠后地震资料上表现为空白反射的区域，有了强弱交替的地震反射波组出现，反映了砂泥岩岩性变化的特征。特别是在水道内部，由原始地震剖面上呈弱振幅-空白反射区增强为内部“V”型和“U”型充填反射特征，反映了河道内部多期充填变化过程。

通过对比原始叠后地震剖面和HHT 变换后的地震剖面，HHT 变换后的地震剖面无论是在地震反射连续性和频率上都有了明显的提高。对比过水下分支河道的地震剖面（图6，图7），HHT 后的地震剖面上，河道顶底面的强反射仍然清晰，同时，能够追溯河道在横向上的延伸范围，可较清楚地识别河道内部层状充填与河道边缘楔状反射特征，确定河道砂体尖灭位置。

图6 过水下分支河道地震剖面HHT 第一分量结果Fig. 6 HHT first component results of underwater branch channel seismic profile

图7 HHT 变换后地震剖面（b）与原始叠后地震剖面（a）河道内部充填细节对比（位置见图4 和图6）Fig. 7 Comparison of filling details between HHT transformed seismic profile (b) and original post stack seismic profile (a) (see Fig. 4 and Fig. 6 for location)

除了河道内部反射有了较大提高外，围岩的反射特征也有了较大程度地提高，原始地震剖面上不连续的反射体能够连续追踪，同时，原本振幅较弱或者空白反射区域，代表岩性变化不大或者岩性较均一的层段，出现了高频的地震反射体，说明这些层段和区域的岩性有细微的变化，符合浅海相、海陆过渡相沉积的特征。

5 水下分支河道地震剖面特征

利用HHT 后的地震剖面仔细追踪上述地震剖面所横切的河道，有明显的底面侵蚀特征，通过主河道内部波组特征与相互间的切割交错关系，能够识别5 个不同期次河道发育特征反射体（图8）。

图8 河道地震剖面解释成果图与水道发育剖面模式Fig. 8 Interpretation results of channel seismic profile and the development model of channel

一般来讲，水下分流河道与辫状河道特征相似，具有垂向加积、侧向迁移、透镜体叠置等典型地震反射特征[2,17]。在横切分流河道的地震剖面上，除有典型的侵蚀下切特征外，一般表现为双向上超调整，或明显的“底凹顶平”特征，期内部多为强振幅反射，为双向上超充填地震相模式，顶部被上覆泥岩地层所覆盖。在河道拐角部位，由于侧向侵蚀，凹岸一侧产生侧向侵蚀，把泥砂带到凸岸，由凸岸向凹岸倾斜向上发生堆积，砂岩超覆于泥岩之上，形成泥夹砂沉积，后期被上覆的泥岩所覆盖，形成岩性圈闭，在地震剖面上表现为典型单项上超充填地震相。正是由于侧向侵蚀的存在，使得水下分流河道侧向迁移，并且频繁改道。由于频繁迁移，在水下分流河道分支交汇处，往往在若干个以泥岩为主的围岩背景中发育多个透镜体且相互叠置。透镜体在地震剖面上表现为弱振幅-空白反射背景下的一段孤立的、能量向两侧迅速减弱的强振幅弧状地震相。强振幅反射体代表河道砂体，弱振幅-空白反射体代表河道间的富泥质沉积[2]。

从地震剖面上可识别，一期河道规模较大，河道宽度较大，下切深度也较深，呈一“V”字形态，说明改期是河道发育鼎盛期，建造了整个多期河道体系的格架。之后，二期河道在该处继承性发育，对一期河道充填沉积具有侵蚀作用，发育较清晰的侵蚀底界。二期河道在宽度和深度上均较一期河道小很多，呈“U”型叠置在一期河道内。三期河道发生了明显的侧向迁移，向南改道，并且对一期河道和二期河道均有明显的侵蚀，呈“V”型发育在一期河道充填内。三期河道由于向南迁移，侵蚀了二期河道的南部边界，但是规模发生了较大程度的萎缩。四期河道向北迁移，呈“U”型，发育位置与二期河道相当，宽度也与二期河道相近，但深度较二期河道浅，对二期河道也有侵蚀作用，二者之间存在清晰界面；该期河道对三期河道北侧有侵蚀作用，导致三期河道北侧边缘缺失。五期河道继续向北迁移，是该多期河道体系的最末期，也是最北端的河道单元。改期河道宽度较小，呈“U”字型，对四期、二期、一期河道均有侵蚀作用，导致了二期、四期河道的北部边缘缺失。

6 水下分支河道剖面展布特征

地震剖面上，水下分支河道初始发育在双程走时（TWT）0.3 s 层段，通过等时切片可以展示上新统-全新统水下分支河道平面展布特征。0.3 s TWT 地震振幅和相干属性等时切片上可观察到，该时期水下分支河道呈复杂网状：研究区中部发育一条相对大型的主水道，呈NE-SW 向展布，其他水道宽度较小，主要呈NE-SW、近E-W、NW-SE 向展布。基于HHT 后的地震数据体提取的地震属性刻画了水下分支河道沉积体系的细节特征，如水道的延续、水道的切割关系等（图9）；另外，一些在原始地震数据上忽略的小型河道，在HHT 后的地震数据提取的地震属性上被展现了出来。而且相干属性体较振幅属性更加清晰地刻画了水下分支河道平面展布特征。

通过HHT 后的地震等时切片和相干切片可观察到（图9），上新统-全新统水下分支河道体系呈“多向、多支、多态”分布特征。“多向”是指河道具有多个延伸方向，如研究区中部河道以NE-SW 向为主；研究区北部以NW-SE 向为主；研究区南部观察到一条近N-S 向河道。“多支”是指研究区内发育多条河道外，每条河道还具有多个分支，构成了研究区内的河道网格体系；“多态”是指河道的形态多样，既有比较平直的河道，也发育曲率比较大的弯曲河道。这种大量水道发育特征说明该时期河道沉积体系非常发育，且可能具有多个源头，反映该时期研究区可能发育于浅水、近源的沉积环境。

图9 研究区三维地震HHT 变换第一分量0.3 s 等时切片（a 为振幅切片，b 为相干切片）Fig. 9 Isochronous slice of the first component 0.3 s of HHT transformation of 3D seismic data in the study area(a is amplitude slice and b is coherent slice)

7 结论

通过对研究区地震资料进行HHT，提高了地震分辨率，有助于识别海陆过渡相水下分支河道的边界。HHT 后的地震剖面展示了沉积体内部反射特征；通过提取HHT 后的三维地震瞬时地震属性，河道的平面展布特征更加清晰。与常规地震属性对比，经过HHT 的地震资料对沉积体系的精细刻画研究提供了更高分辨率的地震数据。

（1）HHT 第一分量整体频率有了明显的提高，提高了地震分辨率，解释了沉积体系的内部细节特征和边界特征，与原始地震资料相比，更能展示小型水下分支河道的平面展布特征和河道之间的相互切割关系。

（2）HHT 后的地震剖面上，研究区内的典型水下分支河道可识别出5 期发育过程，各期河道底界面具有明显侵蚀特征，晚期河道侵蚀早期河道，具有多期叠置、侧向迁移的特征。

（3）研究区上新统-全新统目的层发育大量水下分支河道，呈复杂网状结构，具有“多向、多支、多态”的沉积体系特征。与海陆过渡相沉积环境相符，反映了沉积时期研究区可能发育于浅水、近源的沉积环境。

致谢：本文中所使用原始叠后地震数据来自于中油国际（新加坡）公司，在地震资料处理过程中，受到了中国石油杭州地质研究院各级领导的悉心帮助，在此一并感谢！