[發明專利]縱波反射系數確定方法、裝置、電子設備及存儲介質

有效

• 申請號: 202011513807.7
• 申請日: 2020-12-18
• 公開（公告）號: CN112649871B
• 公開（公告）日: 2021-08-24
• 發明（設計）人: 王康;彭蘇萍;盧勇旭;崔曉芹
• 申請（專利權）人: 中國礦業大學（北京）
• 主分類號: G01V1/30
• 分類號: G01V1/30;G06F17/11
• 代理公司: 北京超凡宏宇專利代理事務所(特殊普通合伙) 11463
• 代理人: 張萌

1.一種縱波反射系數確定方法，其特征在于，所述方法包括：

獲取表征巖石模型對地震波傳播影響的模量表達式方程組；其中，所述巖石模型中的孔隙之間，以及裂縫與孔隙之間存在液體流動；所述模量表達式方程組中包括初始縱波模量表達式；

基于預先確定的所述巖石模型的地質參數的值和所述模量表達式方程組，得到所述地震波的頻率趨近于0時的多個第一模量值；其中，所述地質參數包括：所述巖石模型的體積模量、所述巖石模型中的巖石顆粒的體積模量、所述巖石模型中的孔隙的孔隙率、孔隙的滲透率、孔隙內的第一流體的體積模量、所述第一流體中的第一液體的粘度、所述第一流體中的氣體的粘度、所述第一流體的飽和度、所述第一流體中的氣泡的平均直徑，縫隙內的第二流體的體積模量、所述第二流體的粘度、裂縫弱度、過量法向柔度、過量切向柔度和巖石模型的拉梅系數；所述第一流體包括：液體；所述第二流體為液體或氣體；

基于所述地質參數的值和所述模量表達式方程組，得到所述地震波的頻率趨近于無窮大時的多個第二模量值；

基于所述地震波的真實頻率、所述地質參數的值和所述初始縱波模量表達式，得到第三模量值；

基于所述多個第一模量值、所述多個第二模量值、所述第三模量值、預先確定的所述地震波的入射角度和預先確定的縱波反射系數的表達式，得到表征所述地震波在所述巖石模型中的衰減程度的縱波反射系數的值。

2.根據權利要求1所述的方法，其特征在于，所述基于所述多個第一模量值、所述多個第二模量值、所述第三模量值、預先確定的所述地震波的入射角度和預先確定的縱波反射系數的表達式，得到表征所述地震波在所述巖石模型中的衰減程度的縱波反射系數的值，包括：

基于松弛函數、所述多個第一模量值、所述多個第二模量值和所述第三模量值，得到多個剛度張量的值；

基于所述多個剛度張量的值、所述第三模量值、所述地震波的入射角度和所述縱波反射系數的表達式，得到所述縱波反射系數的值。

3.根據權利要求2所述的方法，其特征在于，所述基于所述多個剛度張量的值、所述第三模量值、所述地震波的入射角度和所述縱波反射系數的表達式，得到所述縱波反射系數的值，包括：

基于所述多個剛度張量的值、所述第三模量值和垂直橫向各向同性VTI介質的各向異性系數確定表達式，得到多個第一各向異性系數的值；

根據所述多個第一各向異性系數的值、所述地震波的入射角度和所述縱波反射系數的表達式，得到所述縱波反射系數的值。

4.根據權利要求3所述的方法，其特征在于，所述巖石模型包括：第一層模型和第二層模型，所述第二層模型位于所述第一層模型的正下方；所述第一層模型為未包括孔隙或裂縫的巖石層，所述第二層模型為包括孔隙和裂縫的巖石層；所述根據所述多個第一各向異性系數的值、所述地震波的入射角度和所述縱波反射系數的表達式，得到所述縱波反射系數的值，包括：

基于預先確定的拉梅系數矩陣和所述各向異性系數確定表達式，得到多個第二各向異性系數的值；其中，所述拉梅系數矩陣C為λ和μ為所述第一層模型的拉梅系數；

根據所述多個第一各向異性系數的值、所述多個第二各向異性系數的值、所述地震波的入射角度和所述縱波反射系數的表達式，得到所述縱波反射系數的值。