随着油气勘探程度的不断提高，勘探目标越来越复杂，对地震资料分辨率的要求也越来越高。近年来，随着大数据和人工智能技术的发展，机器学习也与地球物理领域有机结合，提高了传统地震信号处理的运算效率，但在分辨率方面的研究和应用还很少。

匹配追踪(Matching Pursuit,MP)算法[1-3]是一种用于信号处理和机器学习的技术，特别是在稀疏表示和信号重构方面，能够把地震信号分解为一系列小波集合，实现信号的时频分解。但是，当前使用匹配追踪算法分解和重建多子波地震信号尚存在两方面问题[4]:(1)地震子波旁瓣振荡严重导致时域追踪精度不足，地震数据分辨率降低；(2)传统匹配追踪算法分解地震道后，频谱成分众多、数据冗余，不能针对性地凸显目标地层的特征，缺乏智能化重建方法。

1901年，Pearson[5]首次提出了主成分分析(Principal Component Analysis,PCA)方法，描述了一种通过维度变换把多元问题转化为更少元素的数学统计方法，用于提取数据特征和建立数理模型。但是PCA是基于线性变换的降维方法，仅适用于线性数据。为了更好地处理非线性数据，Schölkopf等[6]将PCA方法扩展到非线性领域，并引入核函数进行映射，形成了核主成分分析(Kernel Principal Component Analysis,KPCA)方法。KPCA方法更适合处理非线性数据[7]，通过非线性映射把原始空间的数据投影到高维空间，然后在高维空间提取数据特征，再用PCA方法进行数据处理。

本文将匹配追踪算法与KPCA方法应用于地震资料处理，提出了地震信号高分辨率处理方法。首先利用匹配追踪算法通过稀疏分解不断迭代，得到地震信号的最有效信息；然后，将子波替换为宽带Ricker子波进行整形处理，有效压制子波旁瓣，提高地震剖面分辨率；最后，采用KPCA方法将原始地震信号非线性映射到高维空间，在高维空间重建地震信号，消除冗余信息，可更准确地识别薄互层。最终实现地震信号的高精度、高效率、智能化重建，在保持原有信号准确性的同时，大大提高资料分辨率，为断层识别、层位追踪、薄层刻画等后续应用提供可靠的技术支持。

1、原理及实现方法

1.1 匹配追踪算法

在多子波地震道模型中，原始地震信号能够转化为各种相关联的子波原子集合，可以用以下形式刻画

式中：S(t)为地震信号；t为时间；Wi(t)表示第i个反射层的地震子波；Ri(t)表示第i个反射层的反射系数；N (t)表示噪声项；I表示反射层总数。

匹配追踪算法从大量的Morlet小波字典中，自适应地筛选一系列与原始信号匹配的原子[8]。Morlet小波能够精确地刻画实际地震道中存在的高频能量衰减现象[9]，其时间域表达式为

式中：ωm为平均角频率；u是时移；σ是尺度参数，决定波形形状；ψ是相移。

小波的频谱可以表示为

基于Morlet小波的匹配追踪算法步骤包括：(1)对输入地震信号进行Hilbert变换，构建复地震道；(2)刻画Morlet小波，参数u为复地震道最大振幅包络处对应的时间，中心频率ω以及相移ψ都可以用对应时刻的瞬时频率和瞬时相移代替[10];(3)在小波原子库里追踪最佳小波，计算小波振幅值；(4)从当前的残差信号中减去上一步输出的最佳小波，得到新的残差；(5)重复上述步骤，直至满足设定的迭代终止条件。

匹配追踪算法通过稀疏分解，利用小波原子库最大程度地逼近信号的有效分量，从而达到忽略噪声的效果。经过匹配追踪后，地震数据虽然保持了原始地震剖面的反射特征，但并不会改善原本地震子波旁瓣存在的振荡问题，导致时频追踪精度不足[11]、数据的分辨率降低。因此，本文利用子波替换的方法对分解后的子波进行整形。

1.2 子波整形

地震资料的垂向分辨率可以通过子波主旁瓣比评价。当旁瓣能量被压制，地震资料的变化体现在同相轴减少，分辨率提高[12]。

地震子波和反射系数序列在复赛谱上具有可分离性[13]，首先将地震记录在频率域表示出来，可以用与式（1）类似的褶积模型描述[14]

式中：W (t)为地震子波；R (t)为地震反射系数序列。

根据对分辨率和保真度的需要构建均衡Ricker子波，再对地震信号进行反褶积处理，提取反射系数序列，把子波替换为提前构建好的宽带Ricker子波。宽带Ricker子波表示为

式中q和h为Ricker子波频率参数的积分限。

相较于Morlet子波，宽带Ricker子波在提高地震信号的垂向分辨率方面具有显著的优势，在相同条件下，具有主瓣窄、旁瓣弱的特性[15]，有效提高了地震信号的垂向分辨率，解决了子波旁瓣振荡严重导致的精度不足问题。两种子波的波形如图1。

图1 宽带Ricker子波与Morlet子波对比

1.3 KPCA方法

KPCA和PCA有着相似的理论基础，都可以通过维度变换把多元问题转化为更少元素。KPCA是PCA在非线性领域的延伸。

1.3.1 PCA方法原理

PCA的基本思想是提取原始线性数据主要的特征代替原有数据[16-17]，这些特征必须包含原始数据中绝大部分信息，并且彼此之间互不相关。假设有包含了n个样本的样本集{x1,x2,⋯,xn}，其中每个样本点本身还具备p个变量，则

式中

PCA的过程是将原始数据的p个变量重新进行线性组合，从而得到一组p个变量的新组合[18]，

1.3.2 KPCA方法原理

KPCA的核心思想就是引入一个非线性函数ϕ(x)，将原始样本集合X={x1,x2,⋯,xn}从样本空间映射到更高维的特征空间H中[19]。图2给出了KPCA方法的基本思路，不同形状代表不同类别的样本。其中，映射函数为

在特征空间H内再使用线性的PCA对样本集合ϕ(x1),ϕ(x2),⋯,ϕ(xn)进行数据降维，这样问题就变成在特征空间内求解样本协方差矩阵的特征值与特征函数。

特征空间内，样本的协方差矩阵可以表示为

用λ表示协方差矩阵CH的特征值，v表示对应的特征向量，则有

定义核函数

特征向量v在ϕ(x1),ϕ(x2),⋯,ϕ(xn)构成的空间内，因此必定存在系数集as(s=1,⋯,k)，使得

式中k表示组成系数集的子集的数量。

将式(13)和式(14)代入式(12)中，得到

式中a=[a1,⋯,ak]T。因此式(15)等价于

需要指出的是，虽然假设原始数据已经进行过中心化处理，但经过映射后的数据并非零均值化的结果，因此用代替K，对其进行中心化，得到

式中ln为n×n的矩阵且其中元素值都为1/n。

1.4 模型试算

使用模型对子波整形方法进行测试，模型选取主频为50 Hz的Morlet子波，输入后进行子波整形，与同极性的反射系数序列合成为地震记录。处理前、后的合成记录如图3所示，显然，整形后的子波具有更窄的主瓣，调谐现象变弱，能更好地与真实反射系数序列对应，具有较高的分辨率。

为进一步验证方法对薄层的识别能力和降噪效果，利用表1参数构建楔状地质模型。其中，地层(2)为楔状层。使用25 Hz的Ricker子波作为震源，得到原始地震记录，添加45 d B的高斯噪声后，得到如图4a的模型记录。

图2 KPCA方法基本思路示意图

图3 合成记录和反射系数

表1 楔状地质模型参数

对模型数据进行匹配追踪处理，选出20、25、30 Hz的剖面数据，进行子波整形，整形后的数据经KPCA方法处理，得到楔状模型记录(图4b)。

模型中实际尖灭点位于第22道，合成记录中调谐作用出现在第49道(图4a中红色虚线位置)，并将其作为识别出的尖灭点位置；由图4b可见，尖灭点位置出现在第42道。同时，处理后数据增强了弱反射，并剔除了大部分高斯噪声，有利于识别薄层。

图4 楔状模型记录

2、实际应用

现有M区实际模型数据，该区储层整体较薄，平均厚度为17.2 m，主要为颗粒泥晶灰岩和颗粒灰岩，具有多层叠置、横向变化快的特点。同时，该区断层非常发育，有12条主要断层，断层产状较陡。常规的地震勘探技术不能满足油藏描述要求，增加了勘探开发的风险。

因此，本文以该区地震资料为基础，采用基于匹配追踪与KPCA的高分辨率处理方法，从断层识别和薄层刻画的角度，对地震资料开展智能重建研究，以期提高资料分辨率。

2.1 断层识别

在地震勘探过程中，细微断层会导致注水效果变差，影响采油效果[20]。但受地下复杂构造以及地震资料分辨率的影响，很难分辨地层中的细微断层及微小的错动[21-22]。在图5a原始地震剖面中，红色椭圆部分解释为同相轴变形，但细微的错动很难判别是否为断层；经过重建后(图5b)，可以看出相对于原始数据的连续结构，地层存在明显的位移，同相轴发生错断，可以证实为细微断层，蓝色线条为断层解释结果。

图5 原始地震剖面

2.2 薄层刻画

薄互层储集层是油气藏勘探的重要目标，当层的厚度较小、间隔较近时，同相轴受到相邻地层的影响而被掩盖[23]。在资料解释过程中，很难分辨薄互层的底部和顶部界面反射，两个薄层的同相轴相互干扰，从而呈现低频、强振幅反射特征，很容易导致错误的解释。

在原始地震剖面（图6a）中，T10层位的“粗壮”同相轴之间十分模糊，聚焦性较差，薄层分辨率很低；经过本文方法处理后，刻画的层位特征更准确。在重建剖面（图6b）的紫色区域，原本难以判别的薄层变得清晰可辨，同相轴连续完整，伽马测井曲线呈低值响应特征，有效揭示了厚度约为16 m的砂岩薄层顶底界面。因此，高分辨处理后的剖面具有较强的井震相关性且数据可靠。

2.3 频谱分析

在地震记录中，地震信号的能量会聚焦在一定的时频范围内，而噪声信号往往分布在整个时频空间中[24]。在重建信号的过程中，选择部分优势频段，提取原始地震信号的有效频段，使信号能量集中在优势频段内，更能凸显目标储层的构造特征。图7b为高分辨率处理后的井震标定剖面频谱分析结果，相较于原始信号(图7a)，信号能量集中在优势频带附近，分辨地层的能力得到明显提升，有利于特殊地质构造的识别和储层的精准预测。

图6 井震标定剖面

图7 井震标定剖面高分辨率处理前

3、结论

本文将匹配追踪算法和核主成分分析方法应用于地震信号重建，采用子波整形技术对匹配追踪后的子波进行旁瓣抑制，在保证信号保真度的情况下压制了噪声。通过核主成分分析方法重建地震信号，最大程度凸显了信号的主要特征，信号叠加的同相轴更清晰，分辨率得到明显提升。

实际资料应用表明，该方法适用于地震信号的高分辨率处理。经过重建信号，对目标储层的分辨能力得到提升，对地层的构造特征刻画更精细，能有效识别细微断层、薄互层等较难判别的地质现象，为后续的地震资料解释、储层预测提供了可靠的依据。

在实际应用过程中发现方法仍存在一定的局限性，主要体现在以下几个方面：(1)本文的研究目的主要是针对细微断层和薄层的细节刻画，并未验证对其他地质构造的处理效果；(2)匹配追踪及子波整形对较大断层响应较好，经过KPCA重建后对小断层的响应更敏感；(3)该方法虽然提高了薄层的分辨率，但重构信号时选取部分优势频段会在一定程度上丢失一些细节，针对资料保真、保幅仍需进一步攻关研究。

参考文献:

[1]徐天吉,沈忠民,文雪康.多子波分解与重构技术应用研究[J].成都理工大学学报(自然科学版), 2010,37(6):660-665.

[2]朱博华,向雪梅,张卫华.匹配追踪强反射层分离方法及应用[J].石油物探, 2016, 55(2):280-287.

[3]印兴耀,许璐,宗兆云,等.基于局部频率约束的动态快速匹配追踪方法[J].中国石油大学学报(自然科学版), 2018, 42(6):59-66.

[4]刘汉卿,罗明,孙辉,等.应用匹配追踪算法的高分辨率处理方法[J].石油地球物理勘探, 2022, 57(6):1325-1331.

[5]邱娜.地震子波分解与重构技术研究[D].山东青岛:中国海洋大学, 2012.

[7]印兴耀,孔国英,张广智.基于核主成分分析的地震属性优化方法及应用[J].石油地球物理勘探, 2008,43(2):179-183.

[8]张繁昌,兰南英,李传辉,等.地震匹配追踪技术与应用研究进展[J].石油物探, 2020, 59(4):491-504.

[11]段中钰,李婷婷,肖勇,等.基于压缩感知的SRADMM地震数据重建[J].石油地球物理勘探, 2021,56(6):1220-1228.

[12]曹思远,袁殿.高分辨率地震资料处理技术综述[J].新疆石油地质, 2016, 37(1):112-119.

[14]罗有春.提高地震勘探分辨率方法理论与应用研究—反褶积[D].四川成都:成都理工大学, 2009.

[15]赵宝银,陈思远,陶钰,等.应用宽带Ricker子波的期望目标频谱整形[J].石油地球物理勘探, 2020,55(3):541-547.

[16]赵秀红.基于主成分分析的特征提取的研究[D].陕西西安:西安电子科技大学, 2016.

[17]张家乐,程冰洁,徐天吉,等.应用地震属性主成分分析识别古河道[J].石油地球物理勘探, 2023, 58(1):190-195.

[18]郑和忠,魏长江,王树华.基于主成分分析和核主成分分析的地震属性优化的研究[J].青岛大学学报(自然科学版), 2017, 30(3):76-80.

[19]陈国超.基于KPCA和LSSVM的锅炉主汽温异常检测[D].河北保定:华北电力大学, 2021.

[20]房宝才,王长生,刘卿,等.微小断层识别及其对窄薄砂体油田开发的影响[J].大庆石油地质与开发,2003(6):24-26, 72.

[21]杨平,宋强功,詹仕凡,等.基于深度学习的高效构造解释技术研发及工业化应用[J].石油地球物理勘探, 2022, 57(6):1265-1275.

[22]陈俊安,陈海东,龚伟,等.深度学习与边缘增强相结合的断裂综合检测技术——顺北地区超深走滑断裂检测应用实例[J].石油地球物理勘探, 2022, 57(6):1304-1316.

[23]陈晨.基于匹配追踪算法的地震多子波分解与重构技术研究[D].黑龙江大庆:东北石油大学, 2015.

[24]张超.微弱地震信号时频消噪算法研究[D].吉林长春:吉林大学, 2016.

[25]徐龙秀,辛超山,牛东晓,等.基于自适应粒子群参数优化的最小二乘支持向量机用电量预测模型[J].科学技术与工程, 2019, 19(6):136-141.

[26]廖庆陵,窦震海,孙锴,等.基于自适应粒子群算法优化支持向量机的负荷预测[J].现代电子技术, 2022,45(3):125-129.

基金资助:中国石油天然气集团有限公司重大科技专项“多物理场高精度油气地球物理探测技术与装备研究”(2023ZZ05)子课题“多波场地震成像与弹性参数同步反演理论与技术”; 四川省自然科学基金项目“深层—超深层页岩储层“甜点”智能预测及融合表征方法”(2023NSFSC0255)联合资助;

文章来源:党腾雲,徐天吉,钱忠平,等.基于匹配追踪与核主成分分析的地震信号高分辨率处理方法[J].石油地球物理勘探,2024,59(04):782-789.