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

权利要求书

1.一种基于加权叠加的时频域地震资料处理方法，其特征在于，包括以下步骤：

步骤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)e^iφ(t)，则S变换的结果为：

ST{x}(τ,f)=ei(φ(τ)-2πfτ)A(τ)·W(2πf-φ′(τ)2πf)]]>

其中，A(t)和φ(t)分别是待分析单道信号的振幅和相位；W表示S变换中的窗函数的傅里叶变换结果，取实数值；A(τ)和φ(τ)分别是待分析单道信号的瞬时振幅和瞬时相位；

提取时频变换域内的瞬时相位方法如下：

(1)在S变换的结果上乘上一个时频域因子e^i2πfτ，使得S变换结果解析，得到下式：

ST{x}(τ,f)ei2πfτ=eiφ(τ)A(τ)·W(2πf-φ′(τ)2πf);]]>

(2)对上式左边求取瞬时相位可以得到：

Φ[ST{x}(τ,f)e^i2πfτ]＝φ(τ)；

步骤2)通过将瞬时相位叠加而成的时频平面上减去其最小值来构造权值函数；

步骤3)将权值函数进行加权叠加，应用叠加权值进行S反变换，得到叠加后地震数据，输出叠加后的地震波形。