本发明公开了一种基于自适应增益限反Q滤波的地震资料处理方法，包括以下步骤：S1.输入地震资料包含的所有地震数据；S2.选择一道地震数据作为输入信号，通过时频分析求取输入信号不同时间点有效频带的截止频率；S3.通过不同时间有效频带的截止频率算出对应的自适应增益限，得到的振幅补偿函数；S4.根据输入信号的品质，选择吸收衰减补偿的实现方式；S5.对吸收衰减补偿后的结果沿着频率方向求和，得到补偿后的时间域信号；S6.重复步骤S2～S5，依次对地震资料中包含的每一道地震数据进行补偿。本发明可以根据需要灵活的调节振幅补偿函数和合理的选择吸收衰减补偿的实现方式，且自适应增益限可以自适应于地震资料有效频带的截止频率，从而提高地震资料的分辨率。
1.一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：包括以下步骤：S1.输入地震资料包含的所有地震数据；
S2.选择一道地震数据作为输入信号，通过时频分析求取输入信号不同时间点有效频带的截止频率；
所述步骤S2包括以下子步骤：S201.选择一道地震数据作为输入信号，进行时频分析，得到输入信号的时频谱S(t,f)；其中，t表示t时刻，f表示频率；
S202.计算t时刻时频谱S(t,f)的模Gt(f)，并求出t时刻的时频谱模Gt(f)的最大值Gt(fp)：
其中，Max[]表示取最大值；|·|表示取模运算；
S203.计算地震资料有效频带的截止频率fd(t)：其中，d是Gt(f)的最大动态范围，单位分贝；Lt(f)表示所有满足 的频率成分；
S3.通过不同时间有效频带的截止频率算出对应的自适应增益限，得到振幅补偿函数；
S4.若输入信号信噪比较高，选用频率域的实现方式实现吸收衰减补偿，若输入信号的信噪比较低，选用时频域的实现方式实现吸收衰减补偿；
S5.对吸收衰减补偿后的结果沿着频率方向求和，得到补偿后的时间域信号；
S6.重复步骤S2～S5，依次对地震资料中包含的每一道地震数据进行补偿。
2.根据权利要求1所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：所述步骤S3包括以下子步骤：S301.计算自适应增益限C(t)和稳定因子δ3：式中，p，a，b是用户预设的参数，p为实数且p＞1；a、b为实数且b≥a≥0；
Q为地层品质因子，将步骤S2得到的fd(t)作为频率f带入中进行计算，得到γ[t,fd(t)]，由于fd(t)为有效频带的截止频率，即希望进行吸收衰减补偿的有效频带的最大频率，故令γ[t,fd(t)]表征任意时间t的振幅补偿函数的最大值，即自适应增益限C(t)；
S302.获取得到的振幅补偿函数：其中，B3(t,f)为自适应增益限反Q滤波振幅补偿函数。
3.根据权利要求1所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：
所述步骤S4中，判断输入信号信噪比高低的方式如下：S401.测量输入信号中包含的信号能量Psignal和噪音能量Pnoise；
S402.计算输入信号的信噪比SNR：S403.判断信噪比SNR是否大于预设阈值；
4.根据权利要求2所述的一种基于自适应增益限反Q滤波的地震资料处理方法，其特征在于：所述步骤S4中，选用频率域的实现方法进行吸收衰减补偿的公式如下：A(0,f)为输入信号的频谱，h1(t)为频率域处理后的时域信号，Re[]表示取复信号的实部；
选用时频域的实现方法行吸收衰减补偿的公式如下：式中，Nt是Re[H(t,f)]≠0的总数，h2(t)为时频域处理后的时域信号。
一种基于自适应增益限反Q滤波的地震资料处理方法
技术领域
[0001] 本发明属于地震资料数字处理领域，特别是涉及一种基于自适应增益限反Q滤波的地震资料处理方法。
背景技术
[0002] 地震波在地层中传播时，要经受地层Q滤波对地震波造成的与频率和传播时间有关的能量衰减和速度频散；因此随传播时间增大，地震波分辨率逐渐降低且相位畸变越来
越严重。因此，我们可以利用反Q滤波来提高地震资料的分辨率。
[0003] 针对反Q滤波，这种方法能够有效的校正畸变的相位，而且是无条件稳定的，但是振幅补偿函数是一与频率和旅行时相关的e指数函数，在资料处理过程中经常导致数值非
稳定性问题，从而难以获得高分辨率的结果。
[0004] 为了控制反Q滤波的数值非稳定性问题，James and Knight(2003)和王仰华(2006)分别提出稳定因子法反Q滤波。稳定因子法反Q滤波能够控制数值非稳定性问题，但
是以上方法的增益限是时不变的，与地震资料无关；若增益限太小，深层分辨率将降低；若
增益限太大，提高分辨率的同时也会放大高频噪音。
发明内容
[0005] 本发明在现有技术的基础上，提出一种基于自适应增益限反Q滤波的地震资料处理方法，可以根据需要灵活的调节振幅补偿函数和合理的选择吸收衰减补偿的实现方式，
且自适应增益限可以自适应于地震资料有效频带的截止频率，从而提高地震资料的分辨
率。
[0006] 本发明的目的是通过以下技术方案来实现的：一种基于自适应增益限反Q滤波的地震资料处理方法，包括以下步骤：
[0007] S1.输入地震资料包含的所有地震数据；
[0008] S2.选择一道地震数据作为输入信号，通过时频分析求取输入信号不同时间点有效频带的截止频率；
[0009] S3.通过不同时间有效频带的截止频率算出对应的自适应增益限，得到振幅补偿函数；
[0010] S4.若输入信号信噪比较高，选用频率域的实现方法实现吸收衰减补偿，若输入信号的信噪比较低，选用时频域的实现方法实现吸收衰减补偿；
[0011] S5.对吸收衰减补偿后的结果沿着频率方向求和，得到补偿后的时间域信号；
[0012] S6.重复步骤S2～S5，依次对地震资料中包含的每一道地震数据进行补偿。
[0013] 进一步地，所述步骤S2包括以下子步骤：
[0014] S201.选择一道地震数据作为输入信号，进行时频分析，得到输入信号的时频谱S(t,f)；其中，t表示t时刻，f表示频率；
[0015] S202.计算t时刻时频谱S(t,f)的模Gt(f)，并求出t时刻的时频谱模Gt(f)的最大值Gt(fp)：
[0016]
[0017] 其中，Max[]表示取最大值；|·|表示取模运算；
[0018] S203.计算地震资料有效频带的截止频率fd(t)：
[0019]
[0020] 其中，d是Gt(f)的最大动态范围，单位分贝，通常由用户设定；Lt(f)表示所有满足的频率成分。
[0021] 优选地，一般情况下，用户设定的最大动态范围 其中AMAX表示设备可准确记录到的数据的最大值，AMIN表示设备可准确记录到的数据的最小值。
[0022] 进一步地，所述步骤S3包括以下子步骤：
[0023] S301.计算自适应增益限C(t)和稳定因子δ3：
[0024]
[0025] 式中，p，a，b是用户预设的参数，p为实数且p＞1；a、b为实数且b≥a≥0；Q为地层品质因子，将步骤S2得到的fd(t)作为频率f带入
中进行计算，得到γ[t,fd(t)]，由于fd(t)为有效频带的截止频率，即希望
进行吸收衰减补偿的有效频带的最大频率，故令γ[t,fd(t)]表征任意时间t的振幅补偿函
数的最大值，即自适应增益限C(t)；
[0026] S302.获取得到的振幅补偿函数：
[0027]
[0028] 其中，B3(t,f)为自适应增益限反Q滤波振幅补偿函数。
[0029] 进一步地，所述步骤S4中，判断输入信号信噪比高低的方式如下：
[0030] S401.测量输入信号中包含的信号能量Psignal和噪音能量Pnoise；
[0031] S402.计算输入信号的信噪比SNR：
[0032]
[0033] S403.判断信噪比SNR是否大于预设阈值；
[0034] 若是，则输入信号的信噪比较高；
[0035] 若否，则输入信号的信噪比较低。
[0036] 进一步地，所述步骤S4中，选用频率域的实现方法进行吸收衰减补偿的公式如下：
[0037]
[0038] A(0,f)为输入信号的频谱，h1(t)为频率域处理后的时域信号，Re[]表示取复信号的实部；
[0039] 选用时频域的实现方法行吸收衰减补偿的公式如下：
[0040]
[0041] 式中，Nt是Re[H(t,f)]≠0的总数，h2(t)为时频域处理后的时域信号。
[0042] 本发明的有益效果是：
[0043] (1)本发明通过引入三个参数p，a，b，可以根据需要灵活的调节振幅补偿函数，保证了本发明在吸收衰减补偿过程中的灵活性；
[0044] (2)本发明可以根据地震资料品质的差异，合理的选择在计算效率更高的频率域或抗噪性更好的时频域进行吸收衰减补偿，保证了本发明在吸收衰减补偿过程中的高效性
和稳定性。
[0045] (3)本发明所用的自适应增益限能够自适应于地震数据有效频带的截止频率，能够充分补偿有效频带范围内损失的能量，压制有效频带范围外的高频背景噪音，保证了本
发明的吸收衰减补偿的准确性；
附图说明
[0046] 图1为本发明的方法流程图；
[0047] 图2为无噪的合成信号与不同处理方法的处理结果示意图；
[0048] 图3为添加的随机噪音的条件下不同处理方法的处理结果示意图；
[0049] 图4为实施例中输入剖面一及本发明处理结果示意图；
[0050] 图5为实施例中输入剖面二及本发明处理结果示意图。
具体实施方式
[0051] 下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
[0052] 本发明对地震资料进行时频分析求取有效频带的截止频率，根据有效频带的截止频率计算自适应增益限，然后确定该方法的振幅补偿函数，选择频率域或者时频域的实现
方式后，对输入信号进行自适应增益限的反Q滤波，即可得到补偿后的时域信号。
[0053] 一种基于自适应增益限反Q滤波的地震资料处理方法的基本思想：若没有地层Q效应的吸收衰减，深、浅层反射波具有波形相似性；地层Q效应使得在不同时刻的地震波的频
率成分按一定的规律进行衰减，结合地震资料有效频带的变化规律对其乘以时频变的振幅
补偿项，使得深浅层反射波的波形具有相似性即可完成补偿，具体地：
[0054] 如图1所示，一种基于自适应增益限反Q滤波的地震资料处理方法，包括以下步骤：
[0055] S1.输入地震资料包含的所有地震数据；
[0056] S2.选择一道地震数据作为输入信号，通过时频分析求取输入信号不同时间点有效频带的截止频率；在本申请的实施例中，采用改进的短时傅里叶变换进行时频分析，时频
谱S(t,f)表示为：
[0057]
[0058] 其中，τ为时移因子，z(t)为输入信号经过Hilbert变换后的解析信号，w(t)为窗函数；
[0059] S3.通过不同时间有效频带的截止频率算出对应的自适应增益限，得到振幅补偿函数；
[0060] S4.若输入信号信噪比较高，选用频率域的实现方式实现吸收衰减补偿，若输入信号的信噪比较低，选用时频域的实现方式实现吸收衰减补偿；
[0061] S5.对吸收衰减补偿后的结果沿着频率方向求和，得到补偿后的时间域信号；
[0062] S6.重复步骤S2～S5，依次对地震资料中包含的每一道地震数据进行补偿。
[0063] 经过以上处理过程以后，得到了分辨率和信噪比均比较理想的地震剖面，与原始的输入相面相比，处理后的剖面减弱了地层Q效应的影响，精细的刻画了该区域的构造特
征；
[0064] 在本申请的实施例中，图2示出了无噪的合成信号与不同处理方法的处理结果示意图，图2中的(a)图为无噪的合成信号；图2中的(b)图为地层Q效应的结果；图2中的(c)图
为稳定因子法反Q滤波的结果；图2中的(d)图为本发明反Q滤波的结果；可见本发明在频率
域的实现效果优于稳定因子法；图3为添加的随机噪音的条件下不同处理方法的处理结果
示意图，图3中的(a)图为添加的随机噪音；图3中的(b)图为将图3中的(a)图添加到图2中的
(b)图的结果；图3中的(c)图为稳定因子法反Q滤波的结果；图3中的(d)图为本发明反Q滤波
的结果可见本发明在时频域实现的效果优于稳定因子法。图4为实施例中输入剖面一及本
发明处理结果示意图；在本申请的实施例中，图4中的(a)图为输入剖面一；图4中的(b)图为
本发明(频率域)的处理结果；图5为输入剖面二及本发明处理结果示意图，图5中的(a)图为
输入剖面二；图5中的(b)图为本发明(时频域)处理结果；可见，经过本发明处理过程以后，
得到了分辨率和信噪比均比较理想的地震剖面。
[0065] 以上所述是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应该看作是对其他实施例的排除，而可用于其他组合、修改和环境，并能够在本文所
述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的
改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。
