[發(fā)明專利]一種基于地震數(shù)據(jù)的時頻分解與氣藏檢測方法及系統(tǒng)

有效

• 申請?zhí)枺?/td>: 201810801169.5
• 申請日: 2018-07-20
• 公開（公告）號: CN109001800B
• 公開（公告）日: 2020-03-10
• 發(fā)明（設(shè)計）人: 黃昱丞;鄭曉東;段艷廷;潘豪杰;胡蓮蓮
• 申請（專利權(quán)）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 王濤

1.一種基于地震數(shù)據(jù)的時頻分解與氣藏檢測方法，其特征在于，包括：

根據(jù)讀入的地震數(shù)據(jù)，利用短時傅里葉變換及雷尼熵定義生成第一時頻譜；

根據(jù)讀入的地震數(shù)據(jù)，利用不同階數(shù)的連續(xù)小波變換生成各階數(shù)對應(yīng)的第二時頻譜，具體包括：

對每道地震數(shù)據(jù)進行不同階數(shù)的連續(xù)小波變換，生成每道地震數(shù)據(jù)對應(yīng)的不同階數(shù)的小波函數(shù)；

根據(jù)所述不同階數(shù)的小波函數(shù)及預(yù)設(shè)的采樣率與記錄長度，生成每道地震數(shù)據(jù)的各階數(shù)對應(yīng)的第二時頻譜；

將所述第一時頻譜分別與各所述第二時頻譜進行互相關(guān)運算，生成最優(yōu)時頻譜，具體包括：

將每道地震數(shù)據(jù)的對應(yīng)的所述第一時頻譜分別與每道地震數(shù)據(jù)的對應(yīng)的各所述第二時頻譜進行互相關(guān)運算，輸出每道地震數(shù)據(jù)的各所述第二時頻譜對應(yīng)的相關(guān)值；

將每道地震數(shù)據(jù)的最大的相關(guān)值對應(yīng)的第二時頻譜作為最優(yōu)時頻譜，輸出每道地震數(shù)據(jù)對應(yīng)的最優(yōu)時頻譜；

根據(jù)所述最優(yōu)時頻譜，生成衰減梯度數(shù)據(jù)及分頻數(shù)據(jù)，具體包括：

根據(jù)每道地震數(shù)據(jù)對應(yīng)的最優(yōu)時頻譜，利用吸能分析法生成衰減梯度數(shù)據(jù)；

根據(jù)每道地震數(shù)據(jù)對應(yīng)的最優(yōu)時頻譜，生成分頻數(shù)據(jù)；

根據(jù)讀入的層位數(shù)據(jù)、所述衰減梯度數(shù)據(jù)及所述分頻數(shù)據(jù)，生成沿層切片并進行氣藏檢測；

其中，根據(jù)讀入的地震數(shù)據(jù)，利用短時傅里葉變換及雷尼熵定義生成第一時頻譜，具體包括：

根據(jù)讀入的地震數(shù)據(jù)，利用短時傅里葉變換生成每道地震數(shù)據(jù)對應(yīng)的第三時頻譜；根據(jù)各所述第三時頻譜，利用雷尼熵定義生成各所述第三時頻譜對應(yīng)的雷尼熵，將最小的雷尼熵對應(yīng)的第三時頻譜作為所述第一時頻譜。

2.根據(jù)權(quán)利要求1所述的基于地震數(shù)據(jù)的時頻分解與氣藏檢測方法，其特征在于，由若干道所述地震數(shù)據(jù)組成疊后地震數(shù)據(jù)。

3.根據(jù)權(quán)利要求2所述的基于地震數(shù)據(jù)的時頻分解與氣藏檢測方法，其特征在于，所述利用短時傅里葉變換生成每道地震數(shù)據(jù)對應(yīng)的第三時頻譜，包括：根據(jù)選取的窗函數(shù)及時窗長度，對每道地震數(shù)據(jù)進行短時傅里葉變換生成各所述第三時頻譜。

4.一種基于地震數(shù)據(jù)的時頻分解與氣藏檢測系統(tǒng)，其特征在于，包括：

第一生成單元，用于根據(jù)讀入的地震數(shù)據(jù)，利用短時傅里葉變換及雷尼熵定義生成第一時頻譜；

第二生成單元，用于根據(jù)讀入的地震數(shù)據(jù)，利用不同階數(shù)的連續(xù)小波變換生成各階數(shù)對應(yīng)的第二時頻譜；

所述第二生成單元包括：

連續(xù)小波變換模塊，用于對每道地震數(shù)據(jù)進行不同階數(shù)的連續(xù)小波變換，生成每道地震數(shù)據(jù)對應(yīng)的不同階數(shù)的小波函數(shù)；

第二時頻譜生成模塊，用于根據(jù)所述不同階數(shù)的小波函數(shù)及預(yù)設(shè)的采樣率與記錄長度，生成每道地震數(shù)據(jù)的各階數(shù)對應(yīng)的第二時頻譜；

互相關(guān)單元，用于將所述第一時頻譜分別與各所述第二時頻譜進行互相關(guān)運算，生成最優(yōu)時頻譜；

所述互相關(guān)單元包括：

相關(guān)值輸出模塊，用于將每道地震數(shù)據(jù)的對應(yīng)的所述第一時頻譜分別與每道地震數(shù)據(jù)的對應(yīng)的各所述第二時頻譜進行互相關(guān)運算，輸出每道地震數(shù)據(jù)的各所述第二時頻譜對應(yīng)的相關(guān)值；

最優(yōu)時頻譜輸出模塊，用于將每道地震數(shù)據(jù)的最大的相關(guān)值對應(yīng)的第二時頻譜作為最優(yōu)時頻譜，輸出每道地震數(shù)據(jù)對應(yīng)的最優(yōu)時頻譜；

第三生成單元，用于根據(jù)所述最優(yōu)時頻譜，生成衰減梯度數(shù)據(jù)及分頻數(shù)據(jù)；

所述第三生成單元包括：

衰減梯度數(shù)據(jù)生成模塊，用于根據(jù)每道地震數(shù)據(jù)對應(yīng)的最優(yōu)時頻譜，利用吸能分析法生成衰減梯度數(shù)據(jù)；

分頻數(shù)據(jù)生成模塊，用于根據(jù)每道地震數(shù)據(jù)對應(yīng)的最優(yōu)時頻譜，生成分頻數(shù)據(jù)；

氣藏檢測單元，用于根據(jù)讀入的層位數(shù)據(jù)、所述衰減梯度數(shù)據(jù)及所述分頻數(shù)據(jù)，生成沿層切片并進行氣藏檢測；

其中，所述第一生成單元具體包括：

傅里葉變換模塊，用于根據(jù)讀入的地震數(shù)據(jù)，利用短時傅里葉變換生成每道地震數(shù)據(jù)對應(yīng)的第三時頻譜；

第一時頻譜輸出模塊，用于根據(jù)各所述第三時頻譜，利用雷尼熵定義生成各所述第三時頻譜對應(yīng)的雷尼熵，將最小的雷尼熵對應(yīng)的第三時頻譜作為所述第一時頻譜。

5.根據(jù)權(quán)利要求4所述的基于地震數(shù)據(jù)的時頻分解與氣藏檢測系統(tǒng)，其特征在于，由若干道所述地震數(shù)據(jù)組成疊后地震數(shù)據(jù)。