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

说明书

本申请公开了一种薄层识别方法、装置、设备及存储介质，所述方法包括：获取三维地震数据中m个地震道的输入数组，m为正整数；对于所述m个地震道中的每个地震道的输入数组，计算所述输入数组在频率域的广义希尔伯特变换，得到所述地震道的广义希尔伯特数组；对所述地震道的广义希尔伯特数组在时间域进行相位移计算，得到所述地震道的分数阶希尔伯特变换数组；根据所述m个地震道的分数阶希尔伯特变换数组，得到所述三维地震数据的分数阶希尔伯特变换三维数组；提取所述三维地震数据的分数阶希尔伯特变换三维数组中的二维剖面进行薄层识别。

技术领域

本申请涉及油气田勘探领域，特别涉及一种薄层识别方法、装置、设备及存储介质。

背景技术

薄层一般定义是小于地震分辨率极限的地层，如威德斯把薄层定义为小于1/8入射波主波长的地层。随着勘探开发的不断深入，薄层越来越受到重视，作为储集层，其蕴藏着较大规模的油气资源量，作为泥岩盖层，在一定的地质条件下其控制着油气的运移及富集程度。

薄层识别技术是基于对三维地震数据进行高分辨率处理，再提取振幅类属性或采用波阻抗反演方法进行薄层识别。在对三维地震数据进行高分辨率处理的过程中，常使用频率域的分数阶希尔伯特变换方法来进行处理。该分数阶希尔伯特变换进行任意角度的相位移变换，其纵向分辨率受制于常规希尔伯特变换的纵向分辨率精度，存在精度不佳的问题。

发明内容

本申请提供了一种薄层识别方法、装置、设备及存储介质。所述技术方案如下：

根据本申请的一方面，提供了一种薄层识别方法，所述方法包括：

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

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

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

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

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

在本申请的一个可选实现中，所述对于所述m个地震道中的每个地震道的输入数组，计算所述输入数组在频率域的广义希尔伯特变换，得到所述地震道的广义希尔伯特数组，包括：

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

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

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

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

在本申请的一个可选实现中，所述基于所述阶次c计算所述二维数组x(n，ω)对应的广义希尔伯特变换数组X_ght(n)，包括：

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

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

在本申请的一个可选实现中，所述对所述地震道的广义希尔伯特数组在时间域进行相位移计算，得到所述地震道的分数阶希尔伯特变换数组，包括：