[发明专利]基于同相轴分频追踪平滑的浅剖数据机械干扰压制方法

权利要求书

1.基于同相轴分频追踪平滑的浅剖数据机械干扰压制方法，其特征在于所述方法的具体步骤如下：

1)为了获得分频的高信噪比浅剖剖面，设计N组频带由窄到宽的带通滤波参数

式中，f_low(n)与f_high(n)分别为分频滤波的低截止频率与高截止频率；F_low为起始的低截止频率，F_high为起始的高截止频率，要求F_lowF_high；Δf为频带间隔；n为分频序号，且1≤n≤N，随着n的增大频带逐渐变宽；

根据f_low(n)与f_high(n)确定带通滤波窗，对输入的原始剖面s(x_i,t)进行带通滤波，其中i为道号，x_i为第i道记录的坐标，t为旅行时，得到滤波后的剖面s_n(x_i,t)，具体过程可表示为：

公式(2)中的f表示频率；上式是对原始剖面s(x_i,t)做傅里叶变换，公式(2)中的下式是根据给定的频带范围[f_low(n),f_high(n)]进行反傅里叶变换，即在反傅里叶变换过程中实现了基于带通滤波的分频处理；

2)对于每一个分频剖面s_n(x_i,t)中具有一定延续长度的强反射同相轴，通过追踪其最大振幅值来确定该同相轴的旅行时，从而获得一条随坐标x_i变化的旅行时t_n(x_i)曲线，即

t_n(x_i)＝max[s_n(x_i,t)] (3)

3)不存在机械干扰的强反射同相轴及其旅行时曲线是平缓变化的，因此通过样条函数对旅行时t_n(x_i)曲线进行平滑，采用三次样条插值公式进行处理

其中，x表示插值道的坐标，x_i为第i道记录的坐标，x_i+1为第i+1道记录的坐标；t为旅行时，t_n(x_i)是样点x＝x_i的旅行时，t_n(x_i+1)是样点x＝x_i+1的旅行时；m_n(x_i)是t_n(x_i)在样点x＝x_i时的导数，m_n(x_i+1)是t_n(x_i+1)在样点x＝x_i+1时的导数；

4)计算出各分频剖面其旅行时曲线的时间校正量后，所述的时间校正量为平滑前、后的旅行时差，并通过取平均来消除旅行时曲线追踪过程中的随机误差，总的计算公式为

其中Δt(x_i)为时间校正量，t_n(x_i)是随坐标x_i变化的旅行时，t_n′(x_i)是平滑后的旅行时；

5)基于剖面中各地震道的时间校正量Δt(x_i)进行旅行时校正，以消除浅剖数据中因机械干扰导致的同相轴错动现象。