[發明專利]基于深度波數譜的瞬變波包提取方法、存儲介質與終端

有效

• 申請號: 201910040925.1
• 申請日: 2019-01-16
• 公開（公告）號: CN109856672B
• 公開（公告）日: 2019-09-20
• 發明（設計）人: 張繁昌;胥勛勇;朱磊;桑凱恒;唐杰;張佳佳
• 申請（專利權）人: 中國石油大學（華東）
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 濟南圣達知識產權代理有限公司 37221
• 代理人: 李琳

本公開提供了一種基于深度波數譜的瞬變波包提取方法、存儲介質與終端，通過對深度域地震信號進行積分變換，得到深度波數譜；將每一深度點的波數譜變換到同態域，得到對應的波數復賽譜，經低通濾波后反變換為波數譜；確定波數譜中每一深度點的波數譜最大值，根據該最大值對應的波數，計算對應深度的瞬變波包；合成深度域地震記錄，并與觀測地震記錄對比，若差別超過設定值則返回積分變換，改變積分變換的調節因子，重新得到深度波數譜，直到最終計算的差值滿足設定值。

技術領域

本公開涉及一種基于深度波數譜的瞬變波包提取方法、存儲介質與終端。

背景技術

本部分的陳述僅僅是提供了與本公開相關的背景技術信息，不必然構成在先技術。

隨著疊前深度偏移技術的發展，深度域地震資料越來越多。與時間域資料相比，深度域地震資料具有構造真實直觀的特點，可以綜合利用測井資料得到更精確的反映地下地層情況的資料。在深度地震資料反演中，影響反演過程的至關重要一步是深度域合成地震記錄，進行深度域層位標定。深度域層位標定是聯系深度域地震資料和測井資料的橋梁，其準確與否取決于深度域合成地震記錄的好壞，而精確的深度域合成地震記錄依賴于準確的地震瞬變波包的提取。

目前，地震波包提取方法一般是先通過時間域地震資料提取地震子波，然后通過時深轉換及深度域重新采樣，將時間域的子波轉換為深度域的波包，然而由于時間域和深度域地震的非線性特征，用這種方法提取到的波包用于深度域合成記錄誤差較大。同時，由于地層的吸收衰減作用，地震波包的波數在每個深度點是不一樣的，隨深度變化而變化。在深度域用同一個波包進行合成地震記錄不符合實際情況，因此合成地震記錄與觀測地震記錄相似度不高，影響了深度域反演的準確性。

林伯香在《利用變換深度域速度函數制作深度域合成地震記錄》一文中指出，深度域波包是地層速度的函數，其傳播過程不滿足線性時不變條件，必須對深度域速度函數進行適當的變換，使其滿足褶積方法計算合成地震記錄的條件。此后，很多學者通過不同的方法變換速度函數，在深度域提取波包并合成地震記錄。例如張杰在《基于Gibbs抽樣方法的深度域地震子波提取》一文中提出基于褶積理論，首先對層速度和層厚度用一個常速度進行變換，在變換后的常速度下，利用抽樣方法進行波包提取；陳學華等人在《深度域地震子波的優化提取方法》中提出一種深度域速度變換與脊回歸算法相結合的深度域地震子波提取方法。此類深度域波包提取方法都依賴地層速度模型，速度變換方法不同則波包提取結果不同，在實際應用中還存在諸多問題。

發明內容

本公開為了解決上述問題，提出了一種基于深度波數譜的瞬變波包提取方法、存儲介質與終端，本公開利用深度域地震記錄的深度波數譜來提取瞬變波包，實現了深度域瞬變波包的精確提取。利用此瞬變波包通過瞬變褶積模型合成深度地震記錄，與實際地震記錄匹配良好。

根據一些實施例，本公開采用如下技術方案：

一種基于深度波數譜的瞬變波包提取方法，包括：

對深度域地震信號進行積分變換，得到深度波數譜；

將每一深度點的波數譜變換到同態域，得到對應的波數復賽譜，經低通濾波后反變換為波數譜；

確定波數譜中每一深度點的波數譜最大值，根據該最大值對應的波數，計算對應深度的瞬變波包；

合成深度域地震記錄，并與觀測地震記錄對比，若差別超過設定值則返回積分變換，改變積分變換的調節因子，重新得到深度波數譜，直到最終計算的差值滿足設定值。

作為進一步的限定，對深度域地震信號進行積分變換的具體過程為：

式中S(γ,ω)表示深度波數譜；γ表示目標深度點；ω表示圓波數；s(h)表示深度域地震信號；h表示對應的深度范圍；W(h,ω)是深度域高斯窗函數：