[发明专利]一种提高地震数据分辨率的方法

说明书

技术领域

本发明涉及一种地震数据处理方法，特别是关于一种时变带通滤波器与稳定反Q滤波相结合的提高地震数据分辨率的方法。

背景技术

高分辨率地震数据是油藏描述的先决条件，精细油气藏描述、空间非均质性的判断以及监测油气田开发过程中储层流体变化等都需要高分辨率地震数据作为基础。在不改变原有采集方法的情况下，对地震波的衰减进行能量补偿是一种经济有效的获取高品质地震数据的方法。

地下介质的内摩擦力和非均质性引起高频能量耗散和速度频散效应，表现为振幅衰减和波形畸变。经过振幅衰减补偿后的地震数据可以为振幅反演及随后的油藏描述提供可靠的相对振幅信息；经过相位校正后的高分辨率地震数据可以为岩性识别提供正确的时间层位信息。反Q滤波是一种补偿地层的黏弹性衰减提高地震资料分辨率的有效方法。

稳定性和计算效率是任何反Q滤波方法都必须考虑的问题。Gabor域稳定反Q滤波方法，使用Gabor域内的波场代替直接向下延拓的波场，避免了波场延拓的不稳定性，而Gabor变换是一系列Gabor片的傅立叶变换，计算效率很高。对于不含噪声的地震数据，稳定反Q滤波方法能够恢复原则上可以恢复的频率成分；但是，对于含噪地震数据，应用该方法时为了不放大环境噪声，需要使用较大的稳定常数，而这种做法将降低对地震数据的补偿效果。

发明内容

针对上述问题，本发明的目的是提供一种时变带通滤波器与稳定反Q滤波相结合的提高地震数据分辨率的方法。

为实现上述目的，本发明采取以下技术方案：一种提高地震数据分辨率的方法，其包括以下步骤：1）根据输入的地震数据，建立实际衰减地震道的时间域信号u(t)；2）对建立的实际衰减地震道的时间域信号u(t)，构建一系列Gabor片u(τ,t)，即

u(τ,t)＝u(t)h(t-τ)，

其中，h(t-τ)是中心在时间τ处的Gabor分析窗，选用高斯窗函数，即

h(t-τ)=1Tπexp[-(t-τ)2T2],]]>

式中，T为高斯窗的半宽度，t为时间变量，τ为任一时间点，π为圆周率，exp为以无限不循环小数e为底的指数函数；3）对所有时间点τ处的Gabor片u(τ,t)做FFT变换，生成Gabor谱其中，代表不同深度水平的地震波延拓波场的Gabor谱U~(τ,ω)]]>为：

式中，ω为角频率；4）在Gabor谱上应用稳定反Q滤波方法的同时联合时变带通滤波器F(τ,ω)，生成的稳定反Q滤波后的Gabor谱U(τ,ω)为：