薄层识别方法、装置、设备及存储介质制造方法及图纸

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

【技术实现步骤摘要】
薄层识别方法、装置、设备及存储介质


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

技术介绍

[0002]薄层一般定义是小于地震分辨率极限的地层，如威德斯把薄层定义为小于1/8入射波主波长的地层。随着勘探开发的不断深入，薄层越来越受到重视，作为储集层，其蕴藏着较大规模的油气资源量，作为泥岩盖层，在一定的地质条件下其控制着油气的运移及富集程度。
[0003]薄层识别技术是基于对三维地震数据进行高分辨率处理，再提取振幅类属性或采用波阻抗反演方法进行薄层识别。在对三维地震数据进行高分辨率处理的过程中，常使用频率域的分数阶希尔伯特变换方法来进行处理。该分数阶希尔伯特变换进行任意角度的相位移变换，其纵向分辨率受制于常规希尔伯特变换的纵向分辨率精度，存在精度不佳的问题。

技术实现思路

[0004]本申请提供了一种薄层识别方法、装置、设备及存储介质。所述技术方案如下：
[0005]根据本申请的一方面，提供了一种薄层识别方法，所述方法包括：
[0006]获取三维地震数据中m个地震道本文档来自技高网...

【技术保护点】

【技术特征摘要】
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

【专利技术属性】
技术研发人员：张本健，王宇峰，胡欣，裴森奇，尹宏，杨华，邓波，杨迅，陈骁，郑超，孙志韵，王旭丽，李荣容，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：