[發(fā)明專利]薄層識別方法、裝置、設備及存儲介質(zhì)

在審

• 申請?zhí)枺?/td>: 202011324343.5
• 申請日: 2020-11-23
• 公開（公告）號: CN114527505A
• 公開（公告）日: 2022-05-24
• 發(fā)明（設計）人: 張本健;王宇峰;胡欣;裴森奇;尹宏;楊華;鄧波;楊迅;陳驍;鄭超;孫志韻;王旭麗;李榮容
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三高永信知識產(chǎn)權代理有限責任公司 11138
• 代理人: 寧立存

1.一種薄層識別方法，其特征在于，所述方法包括：

獲取三維地震數(shù)據(jù)中m個地震道的輸入數(shù)組，m為正整數(shù)；

對于所述m個地震道中的每個地震道的輸入數(shù)組，計算所述輸入數(shù)組在頻率域的廣義希爾伯特變換，得到所述地震道的廣義希爾伯特數(shù)組；

對所述地震道的廣義希爾伯特數(shù)組在時間域進行相位移計算，得到所述地震道的分數(shù)階希爾伯特變換數(shù)組；

根據(jù)所述m個地震道的分數(shù)階希爾伯特變換數(shù)組，得到所述三維地震數(shù)據(jù)的分數(shù)階希爾伯特變換三維數(shù)組；

提取所述三維地震數(shù)據(jù)的分數(shù)階希爾伯特變換三維數(shù)組中的二維剖面進行薄層識別。

2.根據(jù)權利要求1所述的方法，其特征在于，所述對于所述m個地震道中的每個地震道的輸入數(shù)組，計算所述輸入數(shù)組在頻率域的廣義希爾伯特變換，得到所述地震道的廣義希爾伯特數(shù)組，包括：

對于所述m個地震道中的每個地震道的輸入數(shù)組x(n)，進行加窗傅里葉變換，得到二維數(shù)組x(n，ω)，ω代表頻率；

確定所述廣義希爾伯特變換的階次c；

基于所述階次c計算所述二維數(shù)組x(n，ω)對應的廣義希爾伯特變換數(shù)組X_ght(n)；

其中，ω＝0,1,…,N-1，且ω＝N/2,…,N-1對應負頻率。

3.根據(jù)權利要求2所述的方法，其特征在于，所述基于所述階次c計算所述二維數(shù)組x(n，ω)對應的廣義希爾伯特變換數(shù)組X_ght(n)，包括：

基于如下公式計算所述二維數(shù)組X(n,ω)對應的廣義希爾伯特變換數(shù)組X_ght(n)：

其中，Re[X(n,ω)]代表取復數(shù)X(n,ω)的實部，Re[X(n,0)]代表取復數(shù)X(n,0)的實部，c為階次。

4.根據(jù)權利要求1至3任一所述的方法，其特征在于，所述對所述地震道的廣義希爾伯特數(shù)組在時間域進行相位移計算，得到所述地震道的分數(shù)階希爾伯特變換數(shù)組，包括：

按照如下公式對所述地震道的廣義希爾伯特變換數(shù)組X_ght(n)在時間域進行相位移計算，得到所述地震道的分數(shù)階希爾伯特變換數(shù)組X_fht(n)：

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

其中，v是所述分數(shù)階希爾伯特變換的階數(shù)。

5.根據(jù)權利要求1至3任一所述的方法，其特征在于，所述提取所述三維地震數(shù)據(jù)的分數(shù)階希爾伯特變換三維數(shù)組中的二維剖面進行薄層識別，包括：

基于所述薄層的地震振幅響應特征，提取所述三維地震數(shù)據(jù)的分數(shù)階希爾伯特變換三維數(shù)組中的二維剖面進行薄層識別。

6.一種薄層識別裝置，其特征在于，所述裝置包括：

獲取模塊，用于獲取三維地震數(shù)據(jù)中m個地震道的輸入數(shù)組，m為正整數(shù)；

計算模塊，用于對于所述m個地震道中的每個地震道的輸入數(shù)組，計算所述輸入數(shù)組在頻率域的廣義希爾伯特變換，得到所述地震道的廣義希爾伯特數(shù)組；對所述地震道的廣義希爾伯特數(shù)組在時間域進行相位移計算，得到所述地震道的分數(shù)階希爾伯特變換數(shù)組；根據(jù)所述m個地震道的分數(shù)階希爾伯特變換數(shù)組，得到所述三維地震數(shù)據(jù)的分數(shù)階希爾伯特變換三維數(shù)組；

識別模塊，用于提取所述三維地震數(shù)據(jù)的分數(shù)階希爾伯特變換三維數(shù)組中的二維剖面進行薄層識別。

7.根據(jù)權利要求6所述的裝置，其特征在于，

所述計算模塊，用于對于所述m個地震道中的每個地震道的輸入數(shù)組x(n)，進行加窗傅里葉變換，得到二維數(shù)組x(n，ω)，ω代表頻率；確定所述廣義希爾伯特變換的階次c；基于所述階次c計算所述二維數(shù)組x(n，ω)對應的廣義希爾伯特變換數(shù)組X_ght(n)；

其中，ω＝0,1,…,N-1，且ω＝N/2,…,N-1對應負頻率。