一种高精度反q滤波地震资料处理方法
【专利摘要】本发明提供了一种高精度反Q滤波地震资料处理方法，属于地震信号处理领域。本发明方法包括：(1)从原始地震道信号u0(τ)的频率域中提取品质因子Q值；(2)通过Gabor正变换将地震信号u0(τ)转换成时频谱<img file="DDA00002263127300011.TIF" wi="177" he="67" />其中τ为旅行时，ω为角频率；(3)对所述时频谱<img file="DDA00002263127300012.TIF" wi="150" he="67" />进行反Q滤波，得到新的时频谱U(τ，ω)；(4)通过Gabor反变换，将U(τ，ω)转换成反Q滤波补偿后的时间域地震道信号u(τ)。利用本发明不仅能够实现全频Q值稳定补偿，而且能够使地震振幅信息尽可能地恢复，在振幅补偿方面，是一种有效的地震信号处理方法。
【专利说明】一种高精度反Q滤波地震资料处理方法
【技术领域】
[0001]本发明属于地震信号处理领域，具体涉及一种高精度反Q滤波地震资料处理方法。
【背景技术】
[0002]当今石油勘探的主要地质目标已经从简单构造油气藏变成了隐蔽性油气藏和复杂构造油气藏，这对勘探的精度提出了更高的要求。然而，地震波在地下介质中传播会产生吸收衰减现象，其振幅随着传播时间路径增大而快速减弱，从而降低了地震资料的信噪比和分辨率。为了解决这种现象，提出了反Q滤波，它能够恢复高频信息，提高弱反射波的能量和地震资料的信噪比、分辨率。
[0003]针对反Q滤波算法，国内外专家、学者进行了大量研究，如Hale提出的基于Futterman数学模型反Q滤波、McCarley提出的基于常数Q值滤波的自适应模型算法、Hargreaves与Calvert提出的基于常数Q值模型的相位反Q滤波、裴江云提出的基于Kjartansson模型的反Q滤波、王仰华提出的一个全局稳定有效反Q滤波等等。其算法大致可以分成三大类:用级数展开作近似高频补偿的反Q滤波方法、基于波场延拓的反Q滤波方法和其他反Q滤波方法。 相比较而言，基于波场延拓的反Q滤波方法由于使用了快速傅里叶变换，计算速度较快，生产效率高，能够有效地校正由频散引起的相位畸变，而且无条件稳定，但是在振幅补偿方面却存在一定的不稳定性，还有一些尚待解决的问题，需要进一步研究探讨，以满足复杂地震资料高分辨率处理的要求。本发明属于基于波场延拓的反Q滤波方法范畴，下面介绍该类型反Q滤波基本原理。
[0004]基于波动理论，平面波U(X，ω)沿波前方向传播距离Λχ后，当前能量为:
[0005]U (χ+ Δ χ, ω) = U (χ, ω) exp [-ik (ω) Δ χ](1)
[0006]其中，i为单位虚数，ω为角频率，k(co)为复数波数。
[0007]当引入品质因子Q，k(co)变形为:
【权利要求】
1.一种高精度反Q滤波地震资料处理方法，其特征在于:所述方法包括以下步骤: (1)，从原始地震道信号1^(τ )的频率域中提取品质因子Q值； (2)，通过Gabor正变换将地震道信号u0(x)转换成时频谱，其中τ为旅行时，ω为角频率； ⑶，利用下式对所述时频谱?}(^_进行反Q滤波，得到新的时频谱υ(τ，ω):
2.根据权利要求1所述的高精度反Q滤波地震资料处理方法，其特征在于:所述步骤(3)中的指数参数kO的取值范围为6~18。
3.根据权利要求2所述的高精度反Q滤波地震资料处理方法，其特征在于:所述步骤(3)中的稳定系数I的取值范围为10_2~10,。
4.根据权利要求1至3任一所述的高精度反Q滤波地震资料处理方法，其特征在于:所述步骤(1)中的提取品质因子Q值采用的是谱比法。
5.根据权利要求1至3任一所述的高精度反Q滤波地震资料处理方法，其特征在于:所述步骤(2)中的Gabor变换采用的是短时窗Gabor变换。
【文档编号】G01V1/28GK103728661SQ201210392936
【公开日】2014年4月16日 申请日期:2012年10月16日 优先权日:2012年10月16日
【发明者】余波 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院