[发明专利]基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法

说明书

本发明涉及一种基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法，属于地震资料处理与分析领域。本发明提出的基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法，具体的实施过程主要分为以下四步：1)高分辨率倾斜叠加谱的创建；2)基于高分辨率倾斜叠加谱进行线性同相轴的追踪；3)利用FK扇形滤波进行线性干扰同相轴衰减；4)进行多次迭代的线性干扰同相轴追踪压制。相较于传统FK视速度滤波方法，本方法仅在线性干扰同相轴所属的短时窗内进行压制处理，降低了损伤有效信号的可能性；通过记录重排将倾斜线性干扰同相轴校正为水平，再利用FK扇形滤波法予以压制，有效避免了滤除大倾角线性干扰时所产生的空间假频。

技术领域

本发明涉及地震资料处理与分析领域，是一种应用于地震数据的倾斜干扰剔除方法。

背景技术

应用于油气勘探的地震勘探方法主要为反射波法，其通常认为一次反射是有效的，而直达波、折射波和面波等均为干扰波。这些同相轴呈线性的强振幅相干噪音，不但会干扰随机噪音衰减、鬼波与多次波压制等步骤的处理效果，还能影响地震成像的真实性与可靠性，进而误导后续的地震地质解释。

对于直达波及浅层折射波等干扰，一般可直接在记录上予以切除，但会损失远炮检距道中的有效反射信号；若采用FK扇形滤波方法对其进行消除，相应压制结果中会产生振幅较强的空间假频。因此，研发一种保幅、高效的强线性干扰的剔除方法，具有重要的现实意义。

发明内容

本发明要解决的技术问题是提供一种基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法。其首先在输入原始炮集记录的基础上，利用同相加权技术创建高分辨率的倾斜叠加谱；然后通过同相轴追踪技术确定时空域中的线性干扰同相轴，最后应用短时窗的FK扇形滤波予以压制。由于仅在线性干扰同相轴所属的短时窗内进行压制处理，不但降低了损伤有效信号的可能性，而且避免了空间假频的产生。

本发明是通过如下技术方案来实现的：

基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法，所述方法的具体步骤如下：

1)高分辨率的倾斜叠加谱创建在野外地震勘探中，通过人工震源激发地震波，获得包含直达波、折射波或面波在内的线性干扰波的地震记录；针对渐变的零偏移距时τ，将输入的原始炮集记录按照一系列的射线速度进行线性时差校正和叠加处理，这个过程即为倾斜叠加变换的正过程

式中：x、t、p和τ分别为偏移距、旅行时、射线速度和零偏移距时，x₁和x₂分别表示最小偏移距和最大偏移距，d和u分别代表偏移距域记录和倾斜叠加记录；

计算各射线速度值的同相加权因子s(p,τ)，即当地震记录中存在该射线速度的同相轴时此加权因子的值较大，否则较小或者为零

式中，x、t、p和τ分别为偏移距、旅行时、射线速度和零偏移距时，d表示偏移距域记录，l表示时窗长度，N_x为偏移距域记录的道数，C为阻尼因子，取地震记录平均振幅的0.01～0.001；

利用式(2)计算的因子s(p,τ)对倾斜叠加道集进行同相加权，并取绝对值即可获得高分辨率的倾斜叠加谱E(p,τ)

E(p,τ)＝|u(p,τ)s(p,τ)| (3)

为了避免对有效波同相轴的误追踪，需要确定其射线速度阈值p_min，然后将倾斜叠加谱E(p,τ)中pp_min的范围充零；