[发明专利]一种基于多次挤压小波变换的地震时频分析方法及系统

权利要求书

1.一种基于多次挤压小波变换的地震时频分析方法，其特征在于，包括以下内容：

1)选取待测三维地震数据体中某一剖面的若干道地震信号；

2)对选取的地震信号分别进行傅里叶变换，得到对应的平均谱图；

3)根据得到的平均谱图，确定选取的各道地震信号所对应的主频区间，并在该主频区间内选取参考频率数值；

4)对选取的各道地震信号分别作多次挤压小波变换，得到各道地震信号多次挤压小波变换的时频谱，具体过程为：

4.1)对地震信号f(t)进行小波变换，得到地震信号的时频表示W(a,b)：

其中，ψ(t)为母小波函数；a和b分别为尺度因子和时间变量；

4.2)根据地震信号的时频表示W(a,b)，计算地震信号的频率估计值

…

其中，为取实部；为求关于变量b的偏导数；为瞬时频率估计算子；M为同步挤压操作的次数；

4.3)根据地震信号的频率估计值计算地震信号多次挤压小波变换的时频谱T^[M](ω,b)：

其中，δ(·)为单位冲激函数，ω为频率，K为实数，T为时间；

5)根据选取的参考频率数值和得到的时频谱，提取待测三维地震数据体中该剖面的三维常频率数据体；

6)对提取的三维常频率数据体进行沿层切片；

7)进入所述步骤1)重新选取待测三维地震数据体中某一剖面的若干道地震信号，直至完成待测三维地震数据体中所有剖面三维常频率数据体的沿层切片。

2.如权利要求1所述的一种基于多次挤压小波变换的地震时频分析方法，其特征在于，所述步骤5)的具体过程为：

5.1)根据选取的参考频率数值和得到的时频谱，提取出待测三维地震数据体中该剖面某道地震信号的常频率切线T^[M](ω₀,b)，其中，ω₀为选取的参考频率数值；

5.2)重复所述步骤5.1)，直至提取出待测三维地震数据体中该剖面各道地震信号的常频率切线，各道地震信号的常频率切线集合为待测三维地震数据体中该剖面的三维常频率数据体。

3.一种计算机设备，其特征在于，包括计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现权利要求1-2中任一项所述的基于多次挤压小波变换的地震时频分析方法对应的步骤。

4.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现权利要求1-2中任一项所述的基于多次挤压小波变换的地震时频分析方法对应的步骤。