[发明专利]一种对地震波信号进行反Q滤波的方法

权利要求书

1.一种对地震波信号进行反Q滤波的方法，该方法包括：

从一地震道采集的地震波信号中提取品质因子Q值；

对所采集的地震波信号进行时频变换得到地震波频谱其中τ为旅行时，ω为角频率；

利用下式对地震波频谱进行反Q滤波以得到反Q滤波后的地震波频谱U(τ，ω)：

上式中γ＝1/(πQ)，并且其中ω_h为地震带宽内的最高频率，σ²是稳定因子，i为虚部单位，指数项中τ′为积分变量，其取值范围为积分区间0～τ，品质因子Q与积分变量τ′有关。

2.根据权利要求1所述的对地震波信号进行反Q滤波的方法，其中从一地震道采集的地震波信号中提取品质因子Q值是基于振幅补偿方法的，包括以下步骤：

利用下式对地震波信号的频域振幅函数A(χ)进行归一化求得归一化后的振幅函数其中χ＝ωτ，A_a为最大振幅值，且当χ＝χ_a时，A(χ)＝A_a；

利用下式基于归一化后的振幅函数求取增益函数Λ_d(χ)：

Λd(χ)=1,χ≤χaA‾(χ)+σ2A‾2(χ)+σ2,χ>χa;]]>

利用下式求取理论补偿函数Λ(χ，Q)：其中β(χ，Q)＝exp[-χ/2Q]；

求取使得时的Q值，该Q值即为所提取的品质因子Q值。

3.根据权利要求2所述的对地震波信号进行反Q滤波的方法，该方法还包括在对频域振幅函数A(χ)进行归一化之前，先将频域振幅函数A(χ)进行中值滤波。

4.根据权利要求1所述的对地震波信号进行反Q滤波的方法，该方法还包括对反Q滤波后的地震波频谱U(τ，ω)进行反时频变换得出反Q滤波后的地震波时域信号。

5.根据权利要求1-4中任一权利要求所述的对地震波信号进行反Q滤波的方法，其中所述时频变换为加伯变换。

6.根据权利要求1所述的对地震波信号进行反Q滤波的方法，其中稳定因子σ²的取值范围为10^-2～10^-10之间。