流體識別方法和裝置的制造方法
【技術領域】
[0001] 本發明涉及地質勘探技術領域，特別涉及一種流體識別方法和裝置。
【背景技術】
[0002] 隨著勘探開發程度的不斷提高，勘探目標逐漸由最初的構造油氣藏轉變為巖性油 氣藏和隱蔽性油氣藏，勘探策略也由最開始的尋找構造圈閉轉變為尋找巖性異常體或者直 接進行流體檢測。近年來，基于地震資料反演儲層彈性參數的方法已被廣泛應用于儲層預 測當中，并且儲層預測的目的也不再局限于構造、巖性圈閉識別，而是要對儲層含流體性 (包括流體類型以及流體飽和度)進行檢測，這就對流體識別技術提出了更高的要求。
[0003] 上世紀70年代，基于地震資料振幅信息的流體識別技術作為首個直接進行地下含 油氣性檢測技術被廣泛應用在實際的勘探開發中。隨后，依據地震波在地下介質中傳播時 會在儲層油氣水分界面處產生沿水平方向振幅強度穩定的反射波同相軸的平點勘探方式 也開始被使用。80年代提出的AV0技術，主要是利用地震波反射振幅隨偏移距或者入射角的 變化規律判斷含氣砂巖儲層，使得利用AV0分析技術檢測流體成為可能。流體因子是利用縱 橫波速度的加權組合來預測儲層流體，隨后將流體因子引入到反射系數域，計算縱橫波反 射系數的加權組合。1999年之后，對流體因子的研究逐漸從反射系數域轉變到波阻抗域，其 中λρ、μρ屬性是阻抗域最主要且應用效果最顯著的流體因子，λρ_μρ屬性交會圖要比常規的 縱橫波阻抗交會圖對儲層流體有更好的識別能力。
[0004] 隨后，Russell等對流體因子進行了進一步的改進得到了適用性更強的流體因子。 近年來，對于流體因子的研究不斷深入，應用領域逐漸從疊后發展到疊前、從各向同性介質 發展到各向異性介質、從時間域發展到頻率域。
[0005] 然而，如何準確有效地對流體進行識別，目前尚未提出有效的解決方案。
【發明內容】
[0006] 本發明實施例提供了一種流體識別方法，以達到準確識別流體的目的，該方法包 括：
[0007] 獲取研究區的疊前角道集地震數據體；
[0008] 對所述疊前角道集地震數據進行分角度疊加，得到分角度疊加數據體；
[0009] 采用所述分角度疊加數據體進行疊前AV0反演，得到縱波阻抗體和密度體；
[0010] 對所述縱波阻抗體進行時頻分解處理，得到分頻波阻抗體；
[0011] 結合所述密度體，對所述分頻波阻抗體進行反演運算，得到所述研究區的儲層地 震頻散屬性；
[0012] 根據所述儲層地震頻散屬性進行流體識別。
[0013] 在一個實施方式中，按照以下公式結合所述密度體，對所述分頻波阻抗體進行反 演運算，得到所述研究區的儲層地震頻散屬性：
[0014]
[0015] 其中，z (t，f)表示特定頻率時的波阻抗體，z (t，fo)表示參考頻率時的波阻抗體，P 表示所述密度體，D表示反演得到的速度頻散值。
[0016] 在一個實施方式中，所述參考頻率為所述疊前角道集地震數據體的主頻率。
[0017] 在一個實施方式中，對所述疊前角道集地震數據進行分角度疊加，得到分角度疊 加數據體，包括：
[0018] 確定分角度疊加的多個角度區間；
[0019] 根據所述多個角度區間對所述疊前角道集地震數據進行分角度疊加，得到所述多 個角度區間中各個角度區間對應的部分疊加數據體。
[0020] 在一個實施方式中，對所述縱波阻抗體進行時頻分解處理，得到分頻波阻抗體，包 括：
[0021 ]通過S變換對所述縱波阻抗體進行時頻分解處理，得到分頻波阻抗體。
[0022] 本發明實施例還提供了一種流體識別裝置，以達到準確識別流體的目的，該裝置 包括：
[0023] 獲取模塊，用于獲取研究區的疊前角道集地震數據體；
[0024] 疊加模塊，用于對所述疊前角道集地震數據進行分角度疊加，得到分角度疊加數 據體；
[0025] 反演模塊，用于采用所述分角度疊加數據體進行疊前AV0反演，得到縱波阻抗體和 密度體；
[0026] 時頻分解模塊，用于對所述縱波阻抗體進行時頻分解處理，得到分頻波阻抗體；
[0027] 運算模塊，用于結合所述密度體，對所述分頻波阻抗體進行反演運算，得到所述研 究區的儲層地震頻散屬性；
[0028] 識別模塊，用于根據所述儲層地震頻散屬性進行流體識別。
[0029] 在一個實施方式中，所述運算模塊具體用于按照以下公式結合所述密度體，對所 述分頻波阻抗體進行反演運算，得到所述研究區的儲層地震頻散屬性：
[0030]
[0031 ]其中，z (t，f)表示特定頻率時的波阻抗體，z (t，f 〇)表示參考頻率時的波阻抗體，P 表示所述密度體，D表示反演得到的速度頻散值。
[0032] 在一個實施方式中，所述參考頻率為所述疊前角道集地震數據體的主頻率。
[0033] 在一個實施方式中，所述疊加模塊包括：
[0034]角度區間確定單元，用于確定分角度疊加的多個角度區間；
[0035] 疊加單元，用于根據所述多個角度區間對所述疊前角道集地震數據進行分角度疊 加，得到所述多個角度區間中各個角度區間對應的部分疊加數據體。
[0036] 在一個實施方式中，所述時頻分解模塊具體用于通過S變換對所述縱波阻抗體進 行時頻分解處理，得到分頻波阻抗體。
[0037] 在本發明實施例中，基于研究區疊前AV0反演得到的縱波阻抗數據體與密度體，利 用時頻分析技術對波阻抗數據體進行時頻分解得到不同瞬時頻率的阻抗體，然后結合地震 波阻抗依賴分頻反演公式提取地震頻散數據體，最后根據頻散屬性實現流體識別，從而實 現對儲層流體分布規律的有效預測，提高了流體識別的準確性。
【附圖說明】
[0038] 此處所說明的附圖用來提供對本發明的進一步理解，構成本申請的一部分，并不 構成對本發明的限定。在附圖中：
[0039] 圖1是根據本發明實施例的流體識別方法流程圖；
[0040] 圖2是根據本發明實施例的儲層飽和不同流體時縱波速度值隨頻率變化關系圖；
[0041] 圖3是根據本發明實施例的儲層飽和不同流體時縱波頻散值隨頻率變化關系圖；
[0042] 圖4是根據本發明實施例的利用地震波阻抗依賴分頻反演得到的地震頻散屬性剖 面示意圖；
[0043] 圖5是根據本發明實施例的利用地震波阻抗依賴分頻反演得到的地震頻散屬性沿 目的層的切片示意圖；
[0044]圖6是根據本發明實施例的流體識別裝置的結構框圖。
【具體實施方式】
[0045] 為使本發明的目的、技術方案和優點更加清楚明白，下面結合實施方式和附圖，對 本發明做進一步詳細說明。在此，本發明的示意性實施方式及其說明用于解釋本發明，但并 不作為對本發明的限定。
[0046] 鑒于地震波在地下含流體介質中傳播時會在不同頻帶范圍和不同尺度上發生反 射波能量衰減與速度頻散，而這種現象主要是由儲層孔隙流體的存在導致的。頻散屬性作 為常規儲層彈性信息的一個延伸和補充，從另一個角度展示了地震數據中被主頻分量所掩 蓋的高頻和低頻信息，對于巖性及隱蔽油氣藏的識別尤其是對流體檢測具有獨特的優勢。 與頻率有關的速度變化主要受到儲層孔隙流體的影響，并且流體的流動性（即巖石滲透率 與流體粘滯性的比值)是影響地震波速度頻變特征的主要因素。
[0047] 發明人考慮到可以基于波阻抗依賴頻散屬性反演進行流體識別，具體的，可以綜 合利用波