[发明专利]一种基于自适应增益限反Q滤波的地震资料处理方法

权利要求书

1.一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：包括以下步骤：

S1.输入地震资料包含的所有地震数据；

S2.选择一道地震数据作为输入信号，通过时频分析求取输入信号不同时间点有效频带的截止频率；

S3.通过不同时间有效频带的截止频率算出对应的自适应增益限，得到的振幅补偿函数；

S4.若输入信号信噪比较高，选用频率域的实现方式实现吸收衰减补偿，若输入信号的信噪比较低，选用时频域的实现方式实现吸收衰减补偿；

S5.对吸收衰减补偿后的结果沿着频率方向求和，得到补偿后的时间域信号；

S6.重复步骤S2～S5，依次对地震资料中包含的每一道地震数据进行补偿。

2.根据权利要求1所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：所述步骤S2包括以下子步骤：

S201.选择一道地震数据作为输入信号，进行时频分析，得到输入信号的时频谱S(t,f)；

S202.计算t时刻时频谱S(t,f)的模G_t(f)，并求出t时刻的时频谱模G_t(f)的最大值G_t(f_p)：

其中，Max[]表示取最大值；|·|表示取模运算；

S203.计算地震资料有效频带的截止频率f_d(t)：

其中，d是G_t(f)的最大动态范围，单位分贝；L_t(f)表示所有满足的频率成分。

3.根据权利要求1所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：所述步骤S3包括以下子步骤：

S301.计算自适应增益限C(t)和稳定因子δ₃：

式中，p，a，b是用户预设的参数，p为实数且p＞1；a、b为实数且b≥a≥0；Q为地层品质因子，将步骤S2得到的f_d(t)作为频率f带入中进行计算，得到γ[t,f_d(t)]，由于f_d(t)为有效频带的截止频率，即希望进行吸收衰减补偿的有效频带的最大频率，故令γ[t,f_d(t)]表征任意时间t的振幅补偿函数的最大值，即自适应增益限C(t)；

S302.获取得到的振幅补偿函数：

其中，B₃(t,f)为自适应增益限反Q滤波振幅补偿函数。

4.根据权利要求1所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：

所述步骤S4中，判断输入信号信噪比高低的方式如下：

S201.测量输入信号中包含的信号能量P_signal和噪音能量P_noise；

S202.计算输入信号的信噪比SNR：

S203.判断信噪比SNR是否大于预设阈值；

若是，则输入信号的信噪比较高；

若否，则输入信号的信噪比较低。

5.根据权利要求1所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：所述步骤S4中，选用频率域的实现方法进行吸收衰减补偿的公式如下：

A(0,f)为输入信号的频谱，h₁(t)为频率域处理后的时域信号，Re[]表示取复信号的实部；

选用时频域的实现方法行吸收衰减补偿的公式如下：

式中，N_t是Re[H(t,f)]≠0的总数，h₂(t)为时频域处理后的时域信号。