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

有效

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

1.一種盲反演方法，其特征在于，包括：

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

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

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

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

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

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

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

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

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

2.如權利要求1所述的方法，其特征在于，將構建的時間維度、線維度及道維度上的一階差分矩陣分別轉換為時間維度、線維度及道維度上的總變分正則化矩陣，包括：

構建時間維度、線維度及道維度上的一階差分矩陣；

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

3.如權利要求1所述的方法，其特征在于，將構建的時間維度、線維度及道維度上的二階差分矩陣分別轉換為時間維度、線維度及道維度上的吉洪諾夫正則化矩陣，包括：

構建時間維度、線維度及道維度上的二階差分矩陣；

根據時間維度、線維度及道維度上的二階差分矩陣和時間維度、線維度及道維度上的單位矩陣，將時間維度、線維度及道維度上的二階差分矩陣分別轉換為時間維度、線維度及道維度上的吉洪諾夫正則化矩陣。

4.如權利要求1所述的方法，其特征在于，根據時間維度、線維度及道維度上的二階差分矩陣和時間維度、線維度及道維度上的單位矩陣，將時間維度、線維度及道維度上的二階差分矩陣分別轉換為時間維度、線維度及道維度上的吉洪諾夫正則化矩陣，包括：

根據道維度上的單位矩陣和線維度上的單位矩陣，及時間維度上的二階差分矩陣，確定時間維度上的吉洪諾夫正則化矩陣；

根據道維度上的單位矩陣和線維度上的二階差分矩陣，及時間維度上的單位矩陣，確定線維度上的吉洪諾夫正則化矩陣；

根據道維度上的二階差分矩陣和線維度上的單位矩陣，及時間維度上的單位矩陣，確定道維度上的吉洪諾夫正則化矩陣。

5.如權利要求1所述的方法，其特征在于，根據反射系數向量和總變分正則化矩陣，及吉洪諾夫正則化矩陣構建波阻抗的目標函數，包括：

根據反射系數向量構建針對橫向突變波阻抗和橫向漸變波阻抗的共同約束項；

根據總變分正則化矩陣構建針對橫向突變波阻抗的總變分正則化約束項；

根據吉洪諾夫正則化矩陣構建針對橫向漸變波阻抗的吉洪諾夫正則化約束項；

根據共同約束項和總變分正則化約束項，及吉洪諾夫正則化約束項構建波阻抗的目標函數。

6.如權利要求1所述的方法，其特征在于，根據波阻抗的目標函數取得最小值時的橫向突變波阻抗和橫向漸變波阻抗確定所述地震記錄的波阻抗，包括：

利用交替方向乘子確定波阻抗的目標函數取得最小值時的橫向突變波阻抗和橫向漸變波阻抗；

將波阻抗的目標函數取得最小值時的橫向突變波阻抗和橫向漸變波阻抗相加，得到所述地震記錄的波阻抗。