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

摘要

本发明涉及一种提高地震数据分辨率的方法，其包括以下步骤：根据输入的地震数据，建立实际衰减地震道的时间域信号u(t)；对建立的实际衰减地震道的时间域信号u(t)，构建一系列Gabor片u(τ,t)；对所有时间点τ处的Gabor片u(τ,t)做FFT变换，生成Gabor谱在Gabor谱上应用稳定反Q滤波方法的同时联合时变带通滤波器F(τ,ω)，生成的稳定反Q滤波后的Gabor谱U(τ,ω)；对经稳定反Q滤波后的Gabor谱U(τ,ω)沿深度位置时间τ方向进行傅立叶逆变换，重构得到的每一Gabor片u(τ,t)；输出采用时变带通滤波器与稳定反Q滤波相结合的方法处理后的高分辨率地震数据。本发明可以广泛应用于对地震数据的处理中。

主权项

1.一种提高地震数据分辨率的方法，其包括以下步骤： 1）根据输入的地震数据，建立实际衰减地震道的时间域信号u(t)； 2）对建立的实际衰减地震道的时间域信号u(t)，构建的一系列Gabor片u(τ,t)为： u(τ,t)＝u(t)h(t-τ)， 其中，h(t-τ)是中心在时间τ处的Gabor分析窗，选用高斯窗函数，即 式中，T为高斯窗的半宽度，t为时间变量，τ为任一时间点，π为圆周率，exp为以无限不循环小数e为底的指数函数； 3）对所有时间点τ处的Gabor片u(τ,t)做FFT变换，生成Gabor谱其中，代表不同深度水平的地震波延拓波场的Gabor谱为：式中，ω为角频率； 4）在Gabor谱上应用稳定反Q滤波方法的同时联合时变带通滤波器F(τ,ω)，生成的稳定反Q滤波后的Gabor谱U(τ,ω)为：式中，为无条件稳定的相位校正函数，其中，ω_h为有效信号带宽的最高频率；γ(τ)=π^-1Q^-1(τ)，Q(τ)为相应地震道的Q模型，τ′为积分过程的中间变量，j为虚数单位；为稳定的振幅补偿函数，其中，是振幅衰减函数，σ²是和信噪比相关的稳定常数，σ²与给定振幅增益限G_lim的经验关系式为：σ²=exp[-(0.23G_lim+1.63)]；是时变带通滤波器的函数表达式；5）对经稳定反Q滤波后的Gabor谱U(τ,ω)沿深度位置时间τ方向进行傅立叶逆变 换，重构得到的每一Gabor片u(τ,t)为： 6）对重构得到的每一Gabor片u(τ,t)沿时间τ积分，得到采用时变带通滤波器与稳定反Q滤波相结合的方法处理后的时间域地震道为：式中，为以Gabor分析窗h(t)的形式表示的Gabor合成窗；7）输出采用时变带通滤波器与稳定反Q滤波相结合的方法处理后的高分辨率地震数据。