[发明专利]一种基于加权叠加的时频域地震资料处理方法

摘要

本发明提供了一种基于加权叠加的时频域地震资料处理方法，包括以下步骤1)将原始地震叠加前信号转化到S变换构成的时频变换域内，提取时频变换域内的瞬时相位；2)构造叠加权值；3)应用叠加权值，得到最终的叠加结果，完成对地震信号的叠加。本发明通过将瞬时相位构成的权值引入到加权叠加过程中，衰减噪声干扰以及增强弱能量有效信号，在弱振幅有效信号能量保持、随机干扰压制和剩余时差校正方面进一步地提高地震资料的叠加效果。

主权项

一种基于加权叠加的时频域地震资料处理方法，其特征在于，包括以下步骤：步骤1)将原始地震叠加前信号转化到S变换构成的时频变换域内，提取时频变换域内的瞬时相位；S变换的表达式为：ST(τ,f)=∫-∞∞x(t)g(τ-t,f)e-i2πftdt=∫-∞∞x(t)|f|2πe-(τ-t)2f2/2e-i2πftdt]]>其中，x(t)为待分析单道信号；ST(τ,f)为时频变换结果；g(τ‑t,f)为窗函数；t为待分析单道信号的时间位置；τ为时频变换结果的时间位置；f为频率；i为虚数单位；假设输入待分析单道信号x(t)＝A(t)eiφ(t)，则S变换的结果为：其中，A(t)和φ(t)分别是待分析单道信号的振幅和相位；W表示S变换中的窗函数的傅里叶变换结果，取实数值；A(τ)和φ(τ)分别是待分析单道信号的瞬时振幅和瞬时相位；提取时频变换域内的瞬时相位方法如下：(1)在S变换的结果上乘上一个时频域因子ei2πfτ，使得S变换结果解析，得到下式：(2)对上式左边求取瞬时相位可以得到：Φ[ST{x}(τ,f)ei2πfτ]＝φ(τ)；步骤2)通过将瞬时相位叠加而成的时频平面上减去其最小值来构造权值函数；步骤3)将权值函数进行加权叠加，应用叠加权值进行S反变换，得到叠加后地震数据，输出叠加后的地震波形。