[發明專利]盲反演方法及裝置

有效

• 申請號: 201910192024.4
• 申請日: 2019-03-14
• 公開（公告）號: CN111694052B
• 公開（公告）日: 2023-04-25
• 發明（設計）人: 郭欣;高建虎;王洪求;雍學善;李勝軍;劉炳楊
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 王天堯;任默聞

本發明公開了一種盲反演方法及裝置，該方法包括：利用稀疏脈沖反演確定地震記錄的反射系數矩陣，將反射系數矩陣轉換為反射系數向量，將構建的時間維度、線維度及道維度上的一階差分矩陣和二階差分矩陣分別轉換為總變分正則化矩陣和吉洪諾夫正則化矩陣，根據反射系數向量和總變分正則化矩陣，及吉洪諾夫正則化矩陣構建波阻抗的目標函數，波阻抗包括橫向突變波阻抗和橫向漸變波阻抗；根據目標函數取得最小值時的橫向突變波阻抗和橫向漸變波阻抗確定地震記錄的波阻抗。本發明利用總變分正則化和吉洪諾夫正則化分別對橫向突變波阻抗和橫向漸變波阻抗進行正則化約束，利用構建的目標函數求解波阻抗，可以提高波阻抗橫向能量的連續性。

技術領域

本發明涉及油氣田勘探技術領域，尤其涉及盲反演方法及裝置。

背景技術

本部分旨在為權利要求書中陳述的本發明實施例提供背景或上下文。此處的描述不因為包括在本部分中就承認是現有技術。

波阻抗反演實際上是從地震剖面上消除子波的影響，留下反射系數，再由反射系數計算出能夠反映地層物性變化的物理參數波阻抗。波阻抗反演在地震勘探領域非常重要，廣泛用于地震解釋及儲層預測，同時在巖石物理方面也得到了很好的應用。

波阻抗反演問題通常會涉及到反卷積的問題，通過從地震記錄中提取地震子波求解出反射系數，然后根據遞推原理得到波阻抗。但是，由于子波的不確定性，導致反射系數在橫向上連續性較差，同時遞推的過程對于噪聲干擾往往非常敏感，使得反演后的波阻抗在橫向上存在能量變化大、連續性差的問題，在地震波阻抗剖面上就會出現“掛面條”現象，不利于儲層預測分析。

因此，現有的波阻抗反演存在波阻抗橫向能量連續性差的問題。

發明內容

本發明實施例提供一種盲反演方法，用以提高波阻抗橫向能量的連續性，該方法包括：

利用稀疏脈沖反演確定地震記錄的反射系數矩陣，將反射系數矩陣按照時間維度、線維度及道維度的順序轉換為反射系數向量；

將構建的時間維度、線維度及道維度上的一階差分矩陣分別轉換為時間維度、線維度及道維度上的總變分正則化矩陣；

將構建的時間維度、線維度及道維度上的二階差分矩陣分別轉換為時間維度、線維度及道維度上的吉洪諾夫正則化矩陣；

根據反射系數向量和總變分正則化矩陣，及吉洪諾夫正則化矩陣構建波阻抗的目標函數；所述波阻抗包括橫向突變波阻抗和橫向漸變波阻抗；

根據波阻抗的目標函數取得最小值時的橫向突變波阻抗和橫向漸變波阻抗確定所述地震記錄的波阻抗；

根據時間維度、線維度及道維度上的一階差分矩陣和時間維度、線維度及道維度上的單位矩陣，將時間維度、線維度及道維度上的一階差分矩陣分別轉換為時間維度、線維度及道維度上的總變分正則化矩陣，包括：

根據道維度上的單位矩陣和線維度上的單位矩陣，及時間維度上的一階差分矩陣，確定時間維度上的總變分正則化矩陣；

根據道維度上的單位矩陣和線維度上的一階差分矩陣，及時間維度上的單位矩陣，確定線維度上的總變分正則化矩陣；

根據道維度上的一階差分矩陣和線維度上的單位矩陣，及時間維度上的單位矩陣，確定道維度上的總變分正則化矩陣。

本發明實施例還提供一種盲反演裝置，用以提高波阻抗橫向能量的連續性，該裝置包括：

反射系數確定模塊，用于利用稀疏脈沖反演確定地震記錄的反射系數矩陣，將反射系數矩陣按照時間維度、線維度及道維度的順序轉換為反射系數向量；

總變分正則化轉換模塊，用于將構建的時間維度、線維度及道維度上的一階差分矩陣分別轉換為時間維度、線維度及道維度上的總變分正則化矩陣；