[發明專利]彈性參數反演方法及裝置

有效

• 申請號: 202011462534.8
• 申請日: 2020-12-11
• 公開（公告）號: CN112596107B
• 公開（公告）日: 2021-09-28
• 發明（設計）人: 鄒安祺;王彥飛
• 申請（專利權）人: 中國科學院地質與地球物理研究所
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/28
• 代理公司: 北京超凡宏宇專利代理事務所(特殊普通合伙) 11463
• 代理人: 安衛靜

1.一種彈性參數反演方法，其特征在于，包括：

獲取多道非平穩地震記錄；

根據建立完成的第一關聯關系模型確定與所述多道非平穩地震記錄對應的彈性波阻抗；其中，所述第一關聯關系模型為：

S_θ＝W_aD_tZ_θ+N_θ

其中，S_θ表示已知的地震剖面，θ表示入射角，N_θ表示隨機噪聲，W_a表示入射角下非平穩地震子波的矩陣形式，D_t是一個差分矩陣，表示在多道非平穩地震記錄縱向時間維度上的差分算子，其中，EI_θ為待求解的彈性波阻抗；

根據建立完成的第二關聯關系模型確定與所述彈性波阻抗對應的彈性參數；其中，所述第一關聯關系模型和所述第二關聯關系模型構成多道非平穩地震記錄模型。

2.根據權利要求1所述的方法，其特征在于，所述根據建立完成的第一關聯關系模型確定與所述多道非平穩地震記錄對應的彈性波阻抗的步驟，包括：

根據所述第一關聯關系模型，將求解所述彈性波阻抗的不適定問題轉化為帶約束的最優化問題；

根據所述多道非平穩地震記錄和目標優化方法求解所述最優化問題，得到所述彈性波阻抗。

3.根據權利要求2所述的方法，其特征在于，所述目標優化方法包括Split Bregman迭代方法，所述根據所述多道非平穩地震記錄和目標優化方法求解所述最優化問題，得到所述彈性波阻抗的步驟，包括：

確定初始化參數；

基于所述初始化參數和所述多道非平穩地震記錄，利用所述Split Bregman迭代方法求解所述最優化問題，得到所述彈性波阻抗。

4.根據權利要求1所述的方法，其特征在于，所述根據建立完成的第二關聯關系模型確定與所述彈性波阻抗對應的彈性參數的步驟，包括：

確定所述第二關聯關系模型的表達形式；其中，所述表達形式包括離散矩陣形式；

根據所述彈性波阻抗和離散矩陣形式的第二關聯關系模型，利用最小二乘法求解，得到所述彈性參數；其中，所述第二關聯關系模型表示彈性波阻抗和彈性參數以離散矩陣的形式對應關聯。

5.根據權利要求1所述的方法，其特征在于，所述多道非平穩地震記錄為不同入射角下非平穩地震子波與地層反射系數褶積的集合。

6.根據權利要求1所述的方法，其特征在于，所述彈性參數包括以下參數中的一種或幾種：縱波速度、橫波速度、密度。

7.一種彈性參數反演裝置，其特征在于，包括：

獲取單元，用于獲取多道非平穩地震記錄；

第一確定單元，用于根據建立完成的第一關聯關系模型確定與所述多道非平穩地震記錄對應的彈性波阻抗；其中，所述第一關聯關系模型為：

S_θ＝W_aD_tZ_θ+N_θ

其中，S_θ表示已知的地震剖面，θ表示入射角，N_θ表示隨機噪聲，W_a表示入射角下非平穩地震子波的矩陣形式，D_t是一個差分矩陣，表示在多道非平穩地震記錄縱向時間維度上的差分算子，其中，EI_θ為待求解的彈性波阻抗；

第二確定單元，用于根據建立完成的第二關聯關系模型確定與所述彈性波阻抗對應的彈性參數；其中，所述第一關聯關系模型和所述第二關聯關系模型構成多道非平穩地震模型。

8.根據權利要求7所述的裝置，其特征在于，第一確定單元，包括：

轉化模塊，根據所述第一關聯關系模型，將求解所述彈性波阻抗的不適定問題轉化為帶約束的最優化問題；

第一求解模塊，根據所述多道非平穩地震記錄和目標優化方法求解所述最優化問題，得到所述彈性波阻抗。