[發明專利]一種基于空間約束的地震反演方法

有效

• 申請號: 201910935313.9
• 申請日: 2019-09-29
• 公開（公告）號: CN112578439B
• 公開（公告）日: 2023-06-02
• 發明（設計）人: 劉嬋娟
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于空間約束的地震反演方法，其特征在于，包括：

基于褶積模型和井旁地震道，提取地震子波；

構建反演目標函數；

基于所述反演目標函數，構建最終反演目標函數；

對所述最終反演目標函數進行求解，獲得反射系數；

其中，所述對所述最終反演目標函數進行求解，獲得反射系數包括：

步驟401：設置對角加權矩陣；

步驟402：基于所述對角加權矩陣求取反射系數；

步驟403：對所述反射系數進行FX空間濾波，計算最終反演目標函數值；

步驟404：更新所述對角加權矩陣，重復執行步驟402-步驟403，直到所述最終反演目標函數值滿足收斂條件，所對應的反射系數即為最終反演的反射系數；

其中，通過以下公式更新對角加權矩陣：

Q_kk＝1/(1+(L_FXm)_k²/σ²)

其中，Q為對角加權矩陣，k為數據列的第k個元素，L_FX為FX濾波算子，m為反射系數，σ為尺度因子。

2.根據權利要求1所述的基于空間約束的地震反演方法，其特征在于，所述基于褶積模型和井旁地震道，提取地震子波包括：

基于所述褶積模型，利用測井資料中的速度信息和密度信息計算反射系數序列，結合所述反射系數序列和井旁地震道提取所述地震子波。

3.根據權利要求1所述的基于空間約束的地震反演方法，其特征在于，所述構建反演目標函數包括：

基于貝葉斯定理，結合似然函數與先驗分布，獲得后驗概率分布；

使所述后驗概率分布極大，獲得所述反演目標函數。

4.根據權利要求3所述的基于空間約束的地震反演方法，其特征在于，所述后驗概率分布為：

其中，p(md)是后驗概率分布，d為觀測數據，G為地震子波褶積矩陣，m為反射系數，C_d為協方差，μ為正則化參數，F(m)是由先驗分布得到的正則化項，其為柯西分布，p為待反演數據的長度，σ為尺度因子。

5.根據權利要求4所述的基于空間約束的地震反演方法，其特征在于，所述反演目標函數為：

其中，J(m)表示反演目標函數。

6.根據權利要求5所述的基于空間約束的地震反演方法，其特征在于，在所述反演目標函數中引入空間平滑約束，獲得最終反演目標函數。

7.根據權利要求6所述的基于空間約束的地震反演方法，其特征在于，所述最終反演目標函數為：

其中，L_FX為FX濾波算子。

8.根據權利要求1所述的基于空間約束的地震反演方法，其特征在于，在所述步驟402中，通過以下公式求取反射系數：

m^j＝(G^TG+μQ^j-1)^-1G^Td^j，

其中，m為反射系數，d為觀測數據，G為地震子波褶積矩陣，j為迭代次數，μ為正則化參數，Q為對角加權矩陣。