[發明專利]高分辨率時變譜模擬方法、裝置、電子設備及介質

在審

• 申請號: 202210143307.1
• 申請日: 2022-02-16
• 公開（公告）號: CN116643315A
• 公開（公告）日: 2023-08-25
• 發明（設計）人: 劉炯;劉喜武;霍志周;張豐麒;劉宇巍;張金強;劉倩
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種高分辨率時變譜模擬方法，其特征在于，包括：

針對地震數據進行Gabor變換，獲得時頻振幅譜與時頻對數譜；

針對所述時頻對數譜進行非線性平滑處理，獲得均值濾波后的時頻對數譜；

根據所述均值濾波后的時頻對數譜，獲得時頻對數譜時間方向上的衰減趨勢，進而計算時變子波時頻對數譜；

根據所述時變子波時頻對數譜，計算時變子波振幅譜；

計算時變的低頻端與高頻端的頻帶拓寬量，確定時變子波期望輸出；

根據所述時變子波振幅譜與所述時變子波期望輸出，計算高分辨率算子；

根據所述時頻振幅譜與所述高分辨率算子，計算高分辨率時變譜結果。

2.根據權利要求1所述的高分辨率時變譜模擬方法，其中，針對所述時頻對數譜進行非線性平滑處理，獲得均值濾波后的時頻對數譜包括：

針對所述時頻對數譜的每一個頻率進行如下步驟，獲得所述均值濾波后的時頻對數譜：

針對所述時頻對數譜開時窗，將所述時窗內數據由小到大排序；

刪除數據總數中設定比例的極小值和極大值部分，計算剩余部分的均值作為時窗中心位置處的濾波輸出；

滑動時窗得到當前頻率的所述時頻對數譜的alpha-trim均值濾波輸出。

3.根據權利要求1所述的高分辨率時變譜模擬方法，其中，針對所述時頻對數譜時間方向上的衰減趨勢，沿頻率方向通過公式(1)進行多項式擬合，獲得所述時變子波時頻對數譜：

其中，Q_tf(t_i,f)為時變子波時頻對數譜，i＝1,2,…,M，M為每道數據的采樣點個數，p_f為多項式階數，f為頻率，c_k為擬合系數。

4.根據權利要求1所述的高分辨率時變譜模擬方法，其中，通過公式(2)計算所述時變子波振幅譜：

w(t,f)＝exp(Q_tf(t,f))???????????????????(2)

其中，w(t,f)為時變子波振幅譜，Q_tf(t_i,f)為時變子波時頻對數譜。

5.根據權利要求1所述的高分辨率時變譜模擬方法，其中，計算時變的低頻端與高頻端的頻帶拓寬量，確定時變子波期望輸出包括：

計算時變的低頻端與高頻端的頻帶拓寬量，將瞬時子波振幅譜的主頻兩側振幅譜分別向低頻端與高頻端移動相應的頻帶拓寬量，移動后中間缺失的振幅譜幅值置為主頻所對應的幅值，以此作為所述時變子波期望輸出。

6.根據權利要求1所述的高分辨率時變譜模擬方法，其中，通過公式(3)計算所述高分辨率算子：

其中，o(t,f)為高分辨率算子，w(t,f)為時變子波振幅譜，w′(t,f)為時變子波期望輸出，σ為穩定性系數。

7.根據權利要求6所述的高分辨率時變譜模擬方法，其中，通過公式(4)計算所述高分辨率時變譜結果：

x(t)＝G^-1[d(t,f)o(t,f)]?????????????????(4)

其中，x(t)為高分辨率時變譜結果，d(t,f)為時頻振幅譜，G^-1表示Gabor反變換。