[發明專利]基于層速度約束的初始模型構建方法及系統

在審

• 申請號: 201910959464.8
• 申請日: 2019-10-10
• 公開（公告）號: CN112649843A
• 公開（公告）日: 2021-04-13
• 發明（設計）人: 周單;呂慧
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/40
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于層速度約束的初始模型構建方法，其特征在于，包括：

針對工區內的每口井的每一地層，針對合成地震記錄進行層位標定，獲得井旁地震道；

計算所述井的所述地層處的地震記錄與井旁地震道之間的地層變換因子；

針對工區內的每口井，針對所述井的每一地層，計算地震記錄與每一個井旁地震道的地震記錄之間的相關系數，進而計算所述相關系數對應的權重；

計算合成縱波阻抗、合成橫波阻抗與合成密度；

獲得歸一化合成縱波阻抗、歸一化合成橫波阻抗、歸一化合成密度與歸一化層速度；

獲得約束合成縱波阻抗、約束合成橫波阻抗與約束合成密度；

獲得最終合成縱波阻抗、最終合成橫波阻抗與最終合成密度，進而獲得初始模型。

2.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(1)計算所述井旁地震道：

s(t)_il＝w(t)_il*r(t)_il (1)

其中，s(t)_il為第i口井、第l地層的井旁地震道，w(t)_il為第i口井、第l地層的地震子波，r(t)_il為第i口井、第l地層的反射系數序列。

3.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(2)計算所述地層變換因子：

其中，為第i口井、第l地層的地層變換因子，seis(t)_il為第i口井、第l地層的地震記錄，s(t)_il為第i口井、第l地層的井旁地震道。

4.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(3)計算所述相關系數：

其中，x_ijl為第i口井、第l地層的第j個地震記錄的樣點值，y_ijl為第i口井、第l地層的第j個井旁地震道的地震記錄的樣點值，γ_il為相關系數，n為為地震記錄的采樣點數。

5.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(4)計算所述相關系數對應的權重：

其中，ε_il為相關系數的權重，γ_il為相關系數。

6.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(5)計算所述合成縱波阻抗：

通過公式(6)計算所述合成橫波阻抗：

通過公式(7)計算所述合成密度：

其中，AI(t)^*為合成縱波阻抗，SI(t)^*為合成橫波阻抗，DEN(t)^*為合成密度。

7.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(8)計算所述歸一化合成縱波阻抗：

通過公式(9)計算所述歸一化合成橫波阻抗：

通過公式(10)計算所述歸一化合成密度：

通過公式(11)計算所述歸一化層速度：

其中，為第k個歸一化合成縱波阻抗，為第k個歸一化合成橫波阻抗，為第k個歸一化合成密度，為第k個歸一化層速度。

8.根據權利要求1所述的基于層速度約束的初始模型構建方法，其中，通過公式(12)計算所述約束合成縱波阻抗：

通過公式(13)計算所述約束合成橫波阻抗：

通過公式(14)計算所述約束合成密度：

其中，為第k個約束合成縱波阻抗，為第k個約束合成橫波阻抗，為第k個約束合成密度，weight為調節因子。