[發明專利]一種基于疊前瞬時頻率屬性分析的儲層烴類預測方法

在審

• 申請號: 201510088511.8
• 申請日: 2015-02-26
• 公開（公告）號: CN104597497A
• 公開（公告）日: 2015-05-06
• 發明（設計）人: 高亞力
• 申請（專利權）人: 浪潮電子信息產業股份有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 濟南信達專利事務所有限公司 37100
• 代理人: 姜明

技術領域

本發明涉及一種石油地球物理勘探技術領域，具體地說是一種基于疊前瞬時頻率屬性分析的儲層烴類預測方法。

背景技術

地震資料解釋是地震數據處理的重要環節，它是將地震資料轉化成我們對勘探區地下地質情況的認識，并從中找到與油氣有關的信息。在相當長的一段時間內，振幅屬性以及AVO特征都是地震資料解釋的重要信息，然而在實際地層中，由于上覆下浮地層介質參數的不同，含氣砂巖所表現出來的AVO異常也不盡相同，除了亮點異常外，還有暗點、相位反轉等特征，從而導致了基于振幅屬性的資料信息在實際應用中的多解性問題。對于某個特定地層，需要有針對性的選擇能夠反映地質現象或油氣信息的屬性加以分析，才能更準確的預測儲層、揭示地下異常。

在疊前地震信息中，瞬時頻率獨立于振幅屬性，有著明確的物理意義，地層的吸收作用使得反射波能量發生衰減，而高頻成分衰減速度大于低頻成分，使得主頻逐漸向低頻移動，特別是當儲層空隙中含有烴類物質時，品質因素Q會變大，從而導致地震反射波主頻有更明顯的降低，這種頻率在空間的變化也已成為指示油氣藏存在的重要指標。

Hilbert變換，也稱希爾伯特變換，是以著名數學家大衛·希爾伯特(David?Hilbert)來命名。希爾伯特變換(hilbert?transform)一個連續時間信號x(t)的希爾伯特變換等于該信號通過具有沖激響應h(t)＝1/πt的線性系統以后的輸出響應x_h(t)。所以信號經希爾伯特變換后，在頻域各頻率分量的幅度保持不變，但相位將出現90°相移。即對正頻率滯后π/2，對負頻率導前π/2，因此希爾伯特變換器又稱為90°移相器。現在的復地震道分析技術主要是建立在Hilbert變換基礎上，將地震波看作是一種解析信號，通過Hilbert變換，獲得解析信號的虛部，進而求得復地震道屬性，根據屬性特點判斷儲層巖性變化，預測油氣。該方法在相當長的一段時間內，被廣泛應用于疊后地震資料中。如何基于疊前瞬時頻率屬性進行分析，以得到儲層烴類預測，是目前還空白的技術。

發明內容

本發明的技術任務是提供一種減小了解釋誤差，提高了烴類預測精度的一種基于疊前瞬時頻率屬性分析的儲層烴類預測方法。

本發明的技術任務是按以下方式實現的：一種基于疊前瞬時頻率屬性分析的儲層烴類預測方法，先對疊前地震信號做加窗處理，消除干擾波對有效信號的影響；再利用復地震道技術，通過Hilbert變換獲取地震信號的瞬時頻率屬性及其隨偏移距的變化信息，通過地震信號的瞬時頻率屬性及其隨偏移距的變化信息得到FVO曲線(FVO--Frequency?Versus?Offset頻率隨偏移距)；將每條FVO曲線順時針旋轉90°，然后將旋轉后多個道集的FVO曲線按順序排列在一起，掛在一個目的儲藏層位上，得到基于瞬時頻率屬性的疊前偏移距體；結合模型和測井資料，找到與油氣相關的變化特征；在基于瞬時頻率的疊前偏移距體中，烴類異常區的瞬時頻率值明顯小于同一偏移距非異常區的瞬時頻率，而且隨著偏移距的增加而減小；從而進行儲層烴類預測。

實現步驟為：

(1)、采用人工激發并由檢波器記錄地震波，取得原始疊前地震數據，去噪保幅處理，抽取疊前地震道集并對其做同相軸拉平處理；

(2)、通過Hilbert變換獲取復地震道信號X(t)＝f(t)+ih(t)，繼而求得瞬時頻率拾取每個CMP道集的FVO曲線；

(3)、將每條FVO曲線順時針旋轉90°，此時數值方向是每個CMP道集的偏移距，水平方向是偏移距對應的瞬時頻率信息；

(4)、將旋轉后的多個CMP道集的FVO曲線順序排列在一起，掛在目的儲層上，由此得到基于瞬時頻率屬性的疊前偏移距體，縱坐標表示偏移距信息，橫坐標表示CMP道號；

(5)、瞬時頻率的突變位置對應了含烴類異常儲層的位置，當儲層中含有某種烴類物質時，地層衰減因子變大，從而使反射波主頻有更明顯的降低，利用這種頻率在空間的變化可指示烴類異常。

在基于高密觀測系統的三層地層模型疊前資料解釋上的應用，具體步驟如下：

(1)、在基于高密觀測系統的三層地層模型中有四個區域，區域1縱波速度Vp＝4600m/s，橫波速度Vs＝Vp/2，密度ρ＝2.2g/cm³，品質因數Q＝100；區域2中Vp＝4500m/s，s＝Vp/2，密度ρ＝2.2g/cm³，品質因數Q＝80；區域3中Vp＝4200m/s，s＝Vp/2，密度ρ＝2.0g/cm³，品質因數Q＝20，區域4中介質參數與區域1相同；