[发明专利]薄层识别方法、装置、设备及存储介质

权利要求书

1.一种薄层识别方法，其特征在于，所述方法包括：

获取三维地震数据中m个地震道的输入数组，m为正整数；

对于所述m个地震道中的每个地震道的输入数组，计算所述输入数组在频率域的广义希尔伯特变换，得到所述地震道的广义希尔伯特数组；

对所述地震道的广义希尔伯特数组在时间域进行相位移计算，得到所述地震道的分数阶希尔伯特变换数组；

根据所述m个地震道的分数阶希尔伯特变换数组，得到所述三维地震数据的分数阶希尔伯特变换三维数组；

提取所述三维地震数据的分数阶希尔伯特变换三维数组中的二维剖面进行薄层识别。

2.根据权利要求1所述的方法，其特征在于，所述对于所述m个地震道中的每个地震道的输入数组，计算所述输入数组在频率域的广义希尔伯特变换，得到所述地震道的广义希尔伯特数组，包括：

对于所述m个地震道中的每个地震道的输入数组x(n)，进行加窗傅里叶变换，得到二维数组x(n，ω)，ω代表频率；

确定所述广义希尔伯特变换的阶次c；

基于所述阶次c计算所述二维数组x(n，ω)对应的广义希尔伯特变换数组X_ght(n)；

其中，ω＝0,1,…,N-1，且ω＝N/2,…,N-1对应负频率。

3.根据权利要求2所述的方法，其特征在于，所述基于所述阶次c计算所述二维数组x(n，ω)对应的广义希尔伯特变换数组X_ght(n)，包括：

基于如下公式计算所述二维数组X(n,ω)对应的广义希尔伯特变换数组X_ght(n)：

其中，Re[X(n,ω)]代表取复数X(n,ω)的实部，Re[X(n,0)]代表取复数X(n,0)的实部，c为阶次。

4.根据权利要求1至3任一所述的方法，其特征在于，所述对所述地震道的广义希尔伯特数组在时间域进行相位移计算，得到所述地震道的分数阶希尔伯特变换数组，包括：

按照如下公式对所述地震道的广义希尔伯特变换数组X_ght(n)在时间域进行相位移计算，得到所述地震道的分数阶希尔伯特变换数组X_fht(n)：

X_fht(n)＝cos(vπ/2)·x(n)+sin(vπ/2)·X_ght(n)；

其中，v是所述分数阶希尔伯特变换的阶数。

5.根据权利要求1至3任一所述的方法，其特征在于，所述提取所述三维地震数据的分数阶希尔伯特变换三维数组中的二维剖面进行薄层识别，包括：

基于所述薄层的地震振幅响应特征，提取所述三维地震数据的分数阶希尔伯特变换三维数组中的二维剖面进行薄层识别。

6.一种薄层识别装置，其特征在于，所述装置包括：

获取模块，用于获取三维地震数据中m个地震道的输入数组，m为正整数；

计算模块，用于对于所述m个地震道中的每个地震道的输入数组，计算所述输入数组在频率域的广义希尔伯特变换，得到所述地震道的广义希尔伯特数组；对所述地震道的广义希尔伯特数组在时间域进行相位移计算，得到所述地震道的分数阶希尔伯特变换数组；根据所述m个地震道的分数阶希尔伯特变换数组，得到所述三维地震数据的分数阶希尔伯特变换三维数组；

识别模块，用于提取所述三维地震数据的分数阶希尔伯特变换三维数组中的二维剖面进行薄层识别。

7.根据权利要求6所述的装置，其特征在于，

所述计算模块，用于对于所述m个地震道中的每个地震道的输入数组x(n)，进行加窗傅里叶变换，得到二维数组x(n，ω)，ω代表频率；确定所述广义希尔伯特变换的阶次c；基于所述阶次c计算所述二维数组x(n，ω)对应的广义希尔伯特变换数组X_ght(n)；

其中，ω＝0,1,…,N-1，且ω＝N/2,…,N-1对应负频率。