[發明專利]多物理模型測井橫波速度曲線重構方法、系統及電子設備

在審

• 申請號: 202111518405.0
• 申請日: 2021-12-13
• 公開（公告）號: CN114280670A
• 公開（公告）日: 2022-04-05
• 發明（設計）人: 李偉東;彭蘇萍;殷裁云;鄒冠貴;柳寶平;莫仕林;王海軍;顧雷雨;趙清全;付康國;曹運飛;李金鑫
• 申請（專利權）人: 華能煤炭技術研究有限公司;中國礦業大學（北京）;云南滇東雨汪能源有限公司
• 主分類號: G01V1/48
• 分類號: G01V1/48
• 代理公司: 北京薈英捷創知識產權代理事務所(普通合伙) 11726
• 代理人: 張陽

1.一種多物理模型測井橫波速度曲線重構方法，其特征在于，所述方法包括：

獲取目標地層的測井數據及巖性物性參數；所述測井數據包括以下至少一項：自然伽馬測井數據、自然電位測井數據、電阻率測井數據、密度測井數據；

根據所述測井數據計算所述目標地層的泥質含量、地層密度和地層孔隙度；

根據所述泥質含量、所述地層密度及所述巖性物性參數，利用預設巖石物理模型，計算所述目標地層的巖層基質模量；所述預設巖石物理模型包括以下任一項：自洽模型、微分等效介質模型、xu-white模型；

將所述地層孔隙度、所述巖層基質模量代入Chapman模型中，得到所述目標地層的巖層等效模量；

根據所述巖層等效模量、巖層等效模量與地層縱橫波速度的預設關系式，構建所述目標地層的橫波速度曲線。

2.根據權利要求1所述的方法，其特征在于，若所述測井數據包括所述目標地層的局部橫波速度測井數據，則所述根據所述測井數據計算所述目標地層的地層孔隙度，包括：

基于所述局部橫波速度測井數據，利用所述預設巖石物理模型對所述目標地層的橫波速度進行反演計算，得到最優收斂結果；

根據所述最優收斂結果，確定所述目標地層的地層孔隙度。

3.根據權利要求2所述的方法，其特征在于，對所述目標地層的橫波速度進行反演計算的公式如下：

其中，縫隙孔隙度表示裂縫孔隙度與裂隙孔隙度之總和，ɑ表示孔隙縱橫比，ε_f表示裂縫密度，v_p—real表示真實縱波速度，v_S—real表示真實橫波速度，v_p-model、v_S-model分別表示基于所述預設巖石物理模型獲得的初步的縱波速度、橫波速度重構結果，a、b、c為各項參數平滑項的權值，分別表示縫隙孔隙度、孔隙縱橫比、裂縫密度在某一深度點的反演值。

4.根據權利要求1所述的方法，其特征在于，所述巖層等效模量包括以下至少一項：巖層等效體積模量、巖層等效剪切模量；所述預設關系式為：

V_p為縱波速度，V_s為橫波速度，K_dry為巖層等效體積模量，μ_dry為巖層等效剪切模量，ρ為電子密度。

5.根據權利要求1所述的方法，其特征在于，所述根據所述測井數據計算所述目標地層的泥質含量，包括：

根據所述自然伽馬測井數據計算所述目標地層的泥質含量；其中，

V_sh為所述泥質含量，GR為自然伽馬測井值，GR_min為純巖性地層自然伽馬最小測井值，GR_max為純巖性地層自然伽馬最大測井值，SH為自然伽馬相對值，GCUR為關于地層年代的經驗系數。

6.根據權利要求1所述的方法，其特征在于，所述根據所述測井數據計算所述目標地層的泥質含量，包括：

根據所述自然電位測井數據計算所述目標地層的泥質含量；其中，

V_sh為所述泥質含量，SP為自然電位測井值，SP_min為純巖性地層自然電位最小測井值，SP_max為純巖性地層自然電位最大測井值。

7.根據權利要求1所述的方法，其特征在于，所述根據所述測井數據計算所述目標地層的泥質含量，包括：

根據所述電阻率測井數據計算所述目標地層的泥質含量；其中，

V_sh為所述泥質含量，R_im為所述目標地層的最大電阻率值，R_sh為泥巖電阻率，R_τ為所述目標地層的電阻率，b為預設系數。