[发明专利]地层旋回自动识别与转换方法、装置及存储介质

权利要求书

1.一种地层旋回自动识别与转换方法，其特征在于，包括：

对地震波形进行时频分析，形成时频谱；

从所述时频谱中心道处的时频波形中提取波形特征点序列；

以所述波形特征点序列中的波形特征点作为种子点，进行时频谱同相轴追踪，形成同相轴趋势线；

确定所述同相轴趋势线的斜率；

根据所述同相轴趋势线的斜率，将所述波形特征点序列转换为地层旋回特征点序列；

将所述地层旋回特征点序列转换为地层旋回曲线；

将地层旋回曲线转换为指定地层旋回单元序列；

其中，所述根据所述同相轴趋势线的斜率，将所述波形特征点序列转换为地层旋回特征点序列，包括：

取出所述波形特征点序列中一个波形特征点的同相轴趋势线；

确认该同相轴趋势线的斜率是否为大于或等于零；

如果该同相轴趋势线的斜率大于零，则将对应波形特征点的旋回类型确定为正旋回；否则，将对应波形特征点的旋回类型确定为反旋回；

按照以上方式，获得所述波形特征点序列中其余波形特征点的旋回类型，从而形成地层旋回特征点序列。

2.如权利要求1所述的地层旋回自动识别与转换方法，其特征在于，在所述对地震波形进行时频分析之前，还包括：

确定时频分析的频率范围和频段数；

根据所述频率范围和所述频段数将所述地震波形划分为多个频段，形成频段集合。

3.如权利要求2所述的地层旋回自动识别与转换方法，其特征在于，所述确定时频分析的频率范围和频段数，包括：

对所述地震波形进行频谱特征分析，以确定所述地震波形的主频带范围；

根据所述主频带范围的低截频和高截频，对应确定时频分析的频率下限和频率上限；

根据公式确定时频分析的频段数；

其中，p为频段数，f_H为频率上限，f_L为频率下限，a为预设的频段宽度。

4.如权利要求2所述的地层旋回自动识别与转换方法，其特征在于，所述对地震波形进行时频分析，形成时频谱，包括：

从所述频段集合中取一个频段；

根据该频段的频率范围对所述地震波形进行时频分析，获得该频段的窄带时频波形；

按照上述方式，获得所述频段集合中其余频段的窄带时频波形，从而形成所述地震波形的时频谱。

5.如权利要求1所述的地层旋回自动识别与转换方法，其特征在于，所述从所述时频谱中心道处的时频波形中提取波形特征点序列，包括：

确定所述时频谱中心道处的时频波形的一阶导数和二阶导数；所述一阶导数对应于所述时频谱中心道处的时频波形的波峰点，所述二阶导数对应于所述时频谱中心道处的时频波形的波谷点；所述波峰点及波谷点作为波形特征点形成波形特征点序列。

6.如权利要求1所述的地层旋回自动识别与转换方法，其特征在于，所述以所述波形特征点序列中的波形特征点作为种子点，进行时频谱同相轴追踪，形成同相轴趋势线，包括：

从所述波形特征点序列中取一个波形特征点作为种子点；

对该种子点进行时频谱同相轴追踪，获得该种子点处的同相轴趋势线；

按照上述方式，获得所述波形特征点序列中其余波形特征点的同相轴趋势线。

7.如权利要求1所述的地层旋回自动识别与转换方法，其特征在于，所述确定所述同相轴趋势线的斜率，包括：

取出所述波形特征点序列中一个波形特征点的同相轴趋势线；

对该同相轴趋势线进行最小二乘法线性拟合，确定该同相轴趋势线的梯度，并将其作为该同相轴趋势线的斜率；

按照以上方式，获得所述波形特征点序列中其余波形特征点的同相轴趋势线的斜率。

8.如权利要求1所述的地层旋回自动识别与转换方法，其特征在于，所述将所述地层旋回特征点序列转换为地层旋回曲线，包括：

对所述地层旋回特征点序列进行线性插值，形成与所述地震波形采样间隔一致的地层旋回曲线。