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

权利要求书

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

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的范围充零；

2)基于高分辨率倾斜叠加谱进行线性同相轴的追踪对E(p,τ)进行适度平滑后，d(x,t)中射线速度值p_i、零偏移距时τ_i的线性同相轴，将在谱E(p,τ)中形成以(p_i,τ_i)为中心极值的团状结构能量；为了保证同相轴追踪过程的稳定性，通过对地震记录观察分析给定同相轴密度N_m参量，其为单位长度时窗内混叠同相轴数目的平均值；应用等值线追踪方法求出各能量团的分布范围，并搜索出其极值点位置，根据该点的坐标(p_i,τ_i)拟合出时空域中的相应同相轴，其所经各道的旅行时为

式中，i表示追踪同相轴的序号，其中i≥1，n表示炮集记录中各道数据的道号，x_n为该道数据的偏移距，p_i和τ_i分别代表所追踪同相轴的射线速度值和零偏移距时；

3)利用FK扇形滤波进行线性干扰同相轴衰减基于同相轴追踪过程确定出多个线性同相轴，将其准确参数(p_i,τ_i)代入公式(4)计算所经地震道的旅行时通过FK扇形滤波法予以消除；针对每个同相轴的滤波过程为：①给定进行线性干扰压制的起始道号n₀；②对于道号n≥n₀的地震道，以t_n为中心截取给定的一个短时窗长度的记录段，使各记录段沿起点位置对齐，从而将目标同相轴校正为水平；③以截取的多道记录段作为输入，通过FK扇形滤波法消除已被校正为水平的同相轴；④将滤波后的记录反重排，并放回各地震道的原时窗位置，其中道号n≥n₀；

4)进行多次迭代的线性干扰同相轴追踪压制采用迭代的追踪压制过程：①通过倾斜叠加谱的观察分析给定谱能量阈值E₀参量，用以界定所追踪同相轴的叠加能量范围；②进行多次迭代的同相轴追踪与压制处理，对于第k次迭代，其中k≥2，基于剩余记录d^k(x,t)创建倾斜叠加谱E^k(p,τ)，若谱中振幅极值E_max不小于E₀表示记录中仍存在较强的线性干扰同相轴，则进行再次追踪压制；③重复步骤②，直至剩余叠加谱中的极值E_max小于阈值E₀为止。