[發明專利]一種疊前地震反演方法及裝置

有效

• 申請號: 201710106720.X
• 申請日: 2017-02-27
• 公開（公告）號: CN106772589B
• 公開（公告）日: 2018-11-02
• 發明（設計）人: 陳雙全;李向陽;支麗霞;宋寶山
• 申請（專利權）人: 中國石油大學(北京)
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 李輝

本申請實施例公開了一種疊前地震反演方法及裝置。所述方法包括：獲取目的工區的疊前地震數據和測井數據；基于所述疊前地震數據和測井數據，確定所述目的工區的初始正演地震數據；基于所述疊前地震數據和所述初始正演地震數據，確定誤差函數項；基于所述誤差函數項，以及與所述目的工區的彈性參數的改變量相關聯的第一特征項、第二特征項和第三特征項，構建所述目的工區的反演目標函數；基于所述反演目標函數，確定所述目的工區的彈性參數的目標值。可以提高疊前地震反演所得到的反演結果的準確度。

技術領域

本申請涉及地震勘探技術領域，特別涉及一種疊前地震反演方法及裝置。

背景技術

將地表接收到的來自地下同一反射點的地震波數據經過動校正后進行疊加，疊加后的地震波數據可以稱為疊后地震數據。疊加前的地震波數據可以稱為疊前地震數據。相對于疊后地震數據，疊前地震數據保留了更多的地下地質信息。例如，與儲層中某一地層位置處的縱波速度、橫波速度或密度等相關的信息。正因如此，基于疊前地震數據的疊前AVO(Amplitude Versus Offset，振幅隨偏移距的變化)反演可以提供儲層中目的層位的縱波速度、橫波速度或密度，其反演精度對儲層預測和油氣藏描述有至關重要的作用。

佐普利茲(Zoeppritz)方程是疊前AVO反演的理論基礎。但由于其復雜性，通常將Zoeppritz方程簡化后得到用于表示地震波反射系數與彈性參數之間的映射關系的Zoeppritz近似公式，基于Zoeppritz近似公式進行疊前AVO反演。彈性參數可以包括：縱波速度、橫波速度和密度。

現有技術中通常采用的疊前AVO反演方法的主要過程是：獲取目的工區的疊前地震數據；在疊前地震數據的縱波入射角為30度以內的假設條件下，基于Zoeppritz近似公式，可以構建目的工區中目的層位的地震波反射系數與彈性參數之間的映射關系；基于疊前地震數據和映射關系，進行AVO反演處理，從而確定目的工區中目的層位的彈性參數的目標值，即縱波速度、橫波速度和密度的目標值。

發明人發現現有技術中至少存在如下問題：針對縱波入射角大于30度的疊前地震數據，基于Zoeppritz近似公式得到的地震波反射系數的準確度較低，可能導致最終得到的反演結果存在較大偏差，即得到目的工區中目的層位的縱波速度、橫波速度和密度的準確度較低。

發明內容

本申請實施例的目的是提供一種疊前地震反演方法及裝置，以提高疊前AVO反演所得到的反演結果的準確度。

為解決上述技術問題，本申請實施例提供一種疊前地震反演方法及裝置是這樣實現的：

一種疊前地震反演方法，包括：

獲取目的工區的疊前地震數據和測井數據；基于所述疊前地震數據和測井數據，確定所述目的工區的初始正演地震數據；

基于所述疊前地震數據和所述初始正演地震數據，確定誤差函數項；

基于所述誤差函數項，以及與所述目的工區的彈性參數的改變量相關聯的第一特征項、第二特征項和第三特征項，構建所述目的工區的反演目標函數；第一特征項用于調整所述反演目標函數的矩陣條件數；第二特征項用于控制所述彈性參數的改變量的范圍；第三特征項用于對所述彈性參數的值進行平滑處理；

基于所述反演目標函數，確定所述目的工區的彈性參數的目標值。

優選方案中，所述基于疊前地震數據和測井數據，確定所述目的工區的初始正演地震數據，包括：

基于所述疊前地震數據，確定所述目的工區的地震子波數據；

基于所述測井數據，確定所述目的工區中各個層位的彈性參數的初始值；

根據預設縱波入射角和所述各個層位的彈性參數的初始值，確定所述目的工區中各個層位的地震波反射系數的初始值；