[发明专利]一种基于高精度时频瞬时相位的地震数据时频谱处理方法

权利要求书

1.一种基于高精度时频瞬时相位的地震数据时频谱处理方法，其特征在于：包括以下步骤：

S1.采集记录待处理的地震原始数据；

S2.对地震原始数据进行时频分析，得到地震原始数据的时频谱G(t,f)，其中t表示时间；f表示频率；

S3.计算时频瞬时相位θ(t,f)和高精度时频瞬时相位θ_h(t,f)；

所述步骤S3包括以下子步骤：

S301.计算时频瞬时相位θ(t,f)；

其中，Im[]和Re[]分别表示取时频谱G(t,f)的虚部和实部；

S302.计算高精度时频瞬时相位θ_h(t,f)；

其中，Δt表示时间采样间隔；||表示取绝对值；Δθ(t,f)为时频瞬时相位增量，表示为：

Δθ(t,f)＝θ(t,f)-θ(t-Δt,f)；

S4.通过高精度时频瞬时相位θ_h(t,f)计算高精度瞬时频率V_h(t,f)；

所述步骤S4中：通过高精度时频瞬时相位计算高精度瞬时频率V_h(t,f)的方式如下；

S5.计算高精度时频瞬时相位的同步挤压变换结果S_h(t,v)，v表示同步挤压变换结果中的频率。

2.根据权利要求1所述的一种基于高精度时频瞬时相位的地震数据时频谱处理方法，其特征在于：所述步骤S2包括以下子步骤：

S201.对地震原始数据进行解析，得到地震原始数据的解析信号z(t)，其中对地震原始数据的解析方法为Hilbert变换；

S202.对地震原始数据的解析信号z(t)进行改进的短时傅里叶变换，得到地震原始数据的时频谱G(t,f)，改进的短时傅里叶变换表示为：

其中，t表示时间；f表示频率；τ是时移量；p(t)表示时窗函数；

3.根据权利要求1所述的一种基于高精度时频瞬时相位的地震数据时频谱处理方法，其特征在于：所述步骤S5中计算高精度时频瞬时相位的同步挤压变换结果S_h(t,v)方式如下；

S_h(t,v)＝∫G(t,f)df,if V_h(t,f)∈[v-Δv,v+Δv]

其中，Δv表示频率增量。