一种基于小波变换的高频延拓方法
【技术领域】
[0001] 本发明属于地震勘探资料处理等数字信号处理领域，具体涉及一种基于小波变换 的高频延拓方法。
【背景技术】
[0002] 吸收衰减严重、地震资料分辨率低是我国西部沙丘和西南碎屑岩探区油气勘探的 一个突出问题。常规基于權积理论的高分辨率处理手段，通常无法保证处理后的保真性和 信噪比。因此，开展能够在保护振幅特征和信噪比前提下，有效提高叠前叠后资料的分辨率 非常重要，可W为高分辨率地震成像和精细储层解释提供高质量基础资料。
[0003] 目前高分辨率处理方法的发展趋势集中在W下几个方面；地层吸收补偿和相位校 正；反權积方法研究；高精度速度分析和动校正方法研究；高精度静校正方法研究；高精度 同相叠加方法研究；叠前、叠后信噪比增强方法研究；叠前、叠后频谱展宽方法研究；高精 度成像方法研究；地震岩性处理技术等，从物探的发展角度来看，高分辨率技术正从中浅层 向深层、从简单构造区块向复杂区块延伸，从时间域成像到深度域成像，而且要求高分辨率 处理与岩性处理相结合。另一方面，Η维VSP、高分辨率井间地震勘探、多波多分量地震勘探 等新技术也在蓬勃发展。
[0004] 高分辨率处理技术已经达到较高水平，各种改进的反權积算法层出不穷，如确定 性子波反權积方法、时变谱白化方法、特征值提高分辨率法、独立分量法（ICA)和盲反權 积等。众多的地球物理工作者多年来致力于高分辨处理的研究，总结出大量的实用处理方 法，并将很多现代科学成就如混沛理论、分形科学、小波理论、神经网络方法、模式识别方法 等应用于地震资料处理领域，有不少方法已经应用于生产实际并取得了较为显著的实际效 果。
[0005] 在地球物理信号处理中，常规基于權积理论的高分辨率处理手段，通常无法保证 处理后的保真性和信噪比，对偏移成像及反演、储层描述带来严重影响。
【发明内容】
[0006] 本发明的目的在于解决上述现有技术中存在的难题，提供一种基于小波变换的高 频延拓方法。利用连续小波变换将输入的地震信号分解为不同尺度的信号，根据优势区域 信号的主频频率对不同的谐波进行频谱外推，得到主谐波和延拓后的谐波，重建地震信号， 得到高分辨率处理后的地震数据，达到了好于商业软件的处理效果。
[0007] 本发明是通过W下技术方案实现的：
[0008] 一种基于小波变换的高频延拓方法，包括：
[0009] (1)逐道输入地震数据，对地震数据进行连续小波变换，从而将地震数据分解在不 同小波尺度上，并获得主轴频率；
[0010] (2)按倍频程逐个频段延拓；在各个小波尺度上，对主轴频率按倍频程在小波域 进行频谱延拓，得到主谐波和一系列延拓谐波；
[0011] (3)重构地震数据：将主谐波和延拓谐波利用小波变换反变换进行重构，得到高 频延拓后的地震数据。
[0012] 所述步骤（1)是送样实现的：
[0013] 首先，利用Morlet小波对地震数据进行连续小波变换：对每一道地震数据逐道进 行小波变换，将地震数据分解在不同小波尺度上；
[0014] 然后，选取优势区域数据进行分析，获得小波变换的频率范围和期望延拓的频谱 范围，从而确定主轴频率值。
[0015] 所述Morlet小波的时间域和频率率显示如下：
[0016]
(4)
[0017] 其中，W是小波变换后的地震信号，f(t)是输入的地震信号，t是时间，τ > 0是 尺度因子，S反映位移，其值为正或负，Ψ (t)是基本小波。
[0018] 所述步骤（2)是送样实现的：
[0019] 外推W主轴频率的小波为基准，按倍频程向高频端延拓，得到延拓谐波，具体如 下：
[0020] 振幅用每个谐波振幅替换超谐波振幅，相位用Φ = Φι 来更新，送样便得 到主谐波和一系列延拓谐波；
[0021 ] 其中妍是主轴频率相位，從是噪声信号相位。
[002引所述步骤（3)中的小波变换反变换如下：
[0023]
(6)
[0024] 其中，f(t)是重构的地震信号，W是小波变换后得到主谐波和延拓谐波信号，t是 时间，τ > 0是尺度因子，S反映位移，其值为正或负，Ψ ω是基本小波。
[0025] 与现有技术相比，本发明的有益效果是：从模拟地震数据拓频效果来看，本方法具 有明显的提高分辨率的优势，能在提高分辨率的同时不产生任何假同相轴，可信度高。
【附图说明】
[0026] 图IMorlet小波波形图
[0027] 图2本发明方法的步骤框图 [002引图3a原始模拟道集；
[0029] 图3b原始模拟道集去掉高频后的道集；
[0030] 图3c高频延拓后的（FSE)道集
[0031] 图4a某区缝洞型储层原始剖面 [003引 图4b F沈后的剖面
[0033] 图4c原始剖面频谱
[0034] 图4d FSE后的剖面频谱
[0035] 图5a实际地震数据实际地震数据FSE效果图
[003引 图化F沈后的剖面
[0037] 图5c原始剖面频谱与FSE后的剖面频谱对比图
[00測图6a原始剖面 [003引 图化F沈后的剖面
[0040] 图6c反權积后的剖面
[0041] 图6d原始剖面频谱与FSE后的剖面频谱对比图
[0042] 图6e原始剖面频谱与反權积后的剖面频谱对比图 [004引图7a原始剖面
[0044] 图化F沈后的剖面
[0045] 图7c原始剖面频谱与FSE后的剖面频谱对比图
【具体实施方式】
[0046] 下面结合附图对本发明作进一步详细描述：
[0047] 本发明能够在保护振幅特征和信噪比前提下，有效提高叠前叠后资料的分辨率， 为高分辨率地震成像和精细储层解释提供高质量基础资料。
[0048] 本发明为信号处理提供了一种高频延拓方法，具体实现步骤是：
[0049] (1)经典小波变换
[0050] 设x(t)是平方可积函数记作（x(t) e L2(R))，抑0是被基本小波或母小波（mother wavelet)的函数。则
[0051]
(1)
[0052] 称为x(t)的小波变换。式中a > 0是尺度因子，τ是反映位移，其值可正可负。 符号〈X，y>代表内积，它的含义是（上标*代表共辆）
[005引 <x(t)，y(t)〉= / x(t)，y*(1:)dt 似
[0054]
量基本小波的位移与尺度伸缩。式（1)中不但t是连续变量， 而且a和τ也是连续变量，因此称为连续的小波变换（CWT)。关于式（1)，有W下几点补充 说明：
[0055] (1)基本小波口的可能是复数信号，特别是解析信号。
[0056] 似尺度因子a的作用是将基本小波料的作伸缩。
[0057] 做AtW前加的因子的目的是使不同a值下的能量保持相等。
[005引 （4)式中的内积往往被不严格的解释成卷积。
[0059] 小波变换在频域上的特点：
[0060] 公式（1)的等效频域表示是：
[0061] 斌
[00的]关于式（3)，可见：
[0063] (1)如果Ψ (ω)是幅频特性比较集中的带通函数，则小波变换便具有表征待分析 信号Χ(ω)频域上局部性质的能力。
[0064] 似采用不同a值作处理时，各Ψ (aw)的中必频率和带宽都不一样，但品质因数 (即（中必频率）/(带宽））却不变。
[0065] 总之，通过连续小波变换，将信号分解在不同的小波尺度上，在各个尺度上进行频 谱延拓，得到一系列延拓谐波，将主谐波和延拓谐波进行重构，得到高频延拓后的地震数 据。
[0066] (2)技术思路及技术实现
[0067] 本发明为为信号处理提供了 一种高频延拓方法。
[0068]