[發明專利]地球物理勘探密度信息提取的方法和裝置

有效

• 申請號: 201711455745.7
• 申請日: 2017-12-27
• 公開（公告）號: CN108345035B
• 公開（公告）日: 2020-06-02
• 發明（設計）人: 馮清源
• 申請（專利權）人: 中國神華能源股份有限公司;神華地質勘查有限責任公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京康信知識產權代理有限責任公司 11240
• 代理人: 趙囡囡

本發明公開了一種地球物理勘探密度信息提取的方法和裝置。其中，該方法包括：獲取目標區的目標地震剖面，其中，目標地震剖面用于表征目標區的地層界面密度變化規律；通過反演模型對目標地震剖面進行疊后反演，得到目標區的地層密度剖面，其中，反演模型基于預設的密度測井數據來建立。本發明解決了現有技術中難以獲取地層密度信息的技術問題。

技術領域

本發明涉及地球物理勘探數據處理領域，具體而言，涉及一種地球物理勘探密度信息提取。

背景技術

縱波速度(Vp)、橫波速度(Vs)、密度(ρ)是表征儲層的三個基本地震彈性參數，其中，密度與儲層孔隙度，飽和度、泥質含量相關系數最高，能夠最直觀地反映儲層和流體的變化規律。

目前，疊前反演是獲取上述參數的最有效手段，其理論基礎為Knott-Zoeppritz方程，由于該方程形式復雜，物理意義不明確，學者們進行了大量的近似處理，其中Aki-Richards近似(1980年)運用最為廣泛，確切地表達了PP波、PS波反射振幅與彈性參數之間的關系。對于縱波(PP波)而言，縱波AVO特征主要取決于地層Vp、Vs、ρ三者的變化率，當地震偏移距有限(近道、中道)時，密度項對反射系數的貢獻權重較小，在實際反演過程中，通常采用速度——密度的擬合關系(如Gardner公式、各工區擬合公式)進行密度項替換或直接省略密度項，將縱波三參數反演降維成兩參數反演，這種近似處理明顯地改善了中小偏移距道集的AVO反演穩定性，不僅提高了反演效率，而且使反演方程的病態性大大削弱，但直接損失就是無法實現精確的密度參數求解。

實現高精度的密度反演將對疊前地震數據用于勘探階段儲層精細描述、開發階段油氣藏動態監測產生深遠的影響，但由于密度反演存在嚴重的病態性，密度參數一直難以獲取。

針對現有技術中難以獲取地層密度信息的問題，目前尚未提出有效的解決方案。

發明內容

本發明實施例提供了一種地球物理勘探密度信息提取的方法和裝置，以至少解決現有技術中難以獲取地層密度信息的技術問題。

根據本發明實施例的一個方面，提供了一種地球物理勘探密度信息提取的方法，包括：獲取目標區的目標地震剖面，其中，目標地震剖面用于表征目標區的地層界面密度變化規律；通過反演模型對目標地震剖面進行疊后反演，得到目標區的地層密度剖面，其中，反演模型基于預設的密度測井數據來建立。

進一步地，獲取目標區中角道集的最優疊加區間；通過對最優疊加區間的角道集進行疊加，得到目標區的目標地震剖面。

進一步地，獲取目標區中角道集的最大疊加區間和最大疊加區間對應的地震信息；根據目標區中角道集的最大疊加區間和最大疊加區間對應的地震信息確定最優疊加區間。

進一步地，獲取目標區地震波入射臨界角和入射角角度改變量；根據入射臨界角和入射角角度改變量確定最大疊加區間。

進一步地，對目標區進行層位標定，得到目標層段；根據目標層段的橫波速度、縱波速度以及第一關系，得到地震波的入射臨界角，其中，第一關系用于表征在剪切模量項系數為零的情況下，入射角與縱橫波速度比的關系。

進一步地，根據目標層段的橫波速度、縱波速度以及第二關系，得到入射角角度改變量，其中，入射角角度改變量包括第一入射角角度改變量和第二入射角角度改變量，第二關系為M((θ₀-Δθ₁),γ)+M((θ₀+Δθ₂),γ)＝0，其中，M用于表征剪切模量項系數，θ₁用于表征第一入射角，θ₂用于表征第二入射角，γ用于表征縱波速度與橫波速度之比。

進一步地，獲取目標區的地震波數據，其中，地震波數據包括：角道集數據、縱波速度、橫波速度和測井解釋結果；對地震波數據進行層位標定。