[发明专利]一种基于时间同步抽取广义Chirplet变换的地震时频分析方法

权利要求书

1.一种基于时间同步抽取广义Chirplet变换的地震时频分析方法，其特征在于，包括以下步骤：

1)假设三维地震资料为S(T,M,N)，其中，T为时间持续长度，M为线号总数，N为道号总数；

2)计算地震数据的若干道地震信号平均傅里叶谱，选取要提取的常频率值η₀；

3)对三维地震资料S(T,M,N)的每一道地震信号f(t)＝S(T,m,n)作时间同步抽取广义Chirplet变换，得到时频谱TGC_f(η,t)，其中，1≤m≤M；1≤n≤N，令S_TF(T,m,n)＝TGC_f(η₀,T)；

得到每一道地震信号的时间同步抽取广义Chirplet变换时频谱包括三个操作：对信号作广义Chirplet变换、计算群延时算子和时间同步抽取，其中：

对信号f(t)进行广义Chirplet变换，得到时频表示GC_f(η,t)，其表达式如下：

式中，g(t)为窗函数，为优化后的Chirp率，η为频率变量，t是时间变量；

计算群延时算子如下式：

其中，表示取虚部，表示求关于频率变量η的偏导数；

时间同步抽取计算如下式：

式中，的定义为：

4)得到常频率数据体为S_TF(T,M,N)，根据需求作岩层切片后续的处理及解释。

2.根据权利要求1所述的基于时间同步抽取广义Chirplet变换的地震时频分析方法，其特征在于，步骤2)中，画出时频谱，观察信号的主频区间，并在此主频区间内选取所要提取的频率值。