[发明专利]基于对数时频域双曲平滑的动态反褶积方法

摘要

本发明涉及一种基于对数时频域双曲平滑的动态反褶积方法，时频谱平滑方法的选择直接决定了吸收衰减和震源子波的质量，对动态反褶积结果具有至关重要的影响。常规时频域的双曲平滑方法主要借助除法实现吸收衰减和震源函数的分离，需要在分母中增加白噪因子以确保数值稳定性。这将明显降低迭代效率，致使吸收衰减与震源子波估求不准确，必将对最终结果造成不利影响。鉴于此，本发明将双曲平滑方法推广到对数时频域，运用线性拟合(舍去截距)与减法实现了吸收衰减和震源函数的完美分离，避免了白噪因子的使用，提高了迭代效率。另外，本发明提出了一种分段线性化方法，克服了实际处理中吸收衰减曲线并非呈单一线性的瓶颈。模型测试和实际资料处理结果从诸多方面都证明了本发明相对于传统方法的优越性。

主权项

基于对数时频域双曲平滑的动态反褶积方法，其特征在于采用以下步骤实现：步骤1：根据POU条件，确定一套分析时窗与合成时窗；步骤2：借助分析时窗与快速傅里叶变换实现Gabor变换，获得输入地震道的时频谱，并对其值谱取对数获得对数频谱；步骤3：根据地震资料的信噪比确定有效振幅水平，确定双曲平滑的时频范围并执行双曲划分；步骤4：用本发明所采用的算法在对数时频域执行双曲平滑，运用迭代思想，通过分段线性拟合并舍去截距的方法，不断优化对吸收衰减的估计。随后移除吸收衰减的影响，根据纵向(时间方向)平均、横向(频率方向)平滑求震源子波；步骤5：采用最小相位假设构建传播子波谱，形成反褶积因子，在时频域执行反褶积；步骤6：借助合成时窗与逆傅里叶变换实现反Gabor变换，将反褶积结果转换到时间域，取其实部作为最终输出结果。