[發明專利]碳酸鹽巖儲層孔隙結構預測方法、裝置、設備及存儲介質

有效

• 申請號: 202010104055.2
• 申請日: 2020-02-20
• 公開（公告）號: CN111239821B
• 公開（公告）日: 2021-02-12
• 發明（設計）人: 趙建國;潘建國;李闖;孫朗秋;劉欣澤;歐陽芳;李智;肖增佳
• 申請（專利權）人: 中國石油大學(北京);中國石油天然氣股份有限公司勘探開發研究院西北分院
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V99/00
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 劉飛;周達

1.一種碳酸鹽巖儲層孔隙結構預測方法，其特征在于，包括：

獲取目標區塊基于地震疊前反演的縱波速度數據體和密度數據體；

根據所述目標區塊的測井數據確定孔隙度與縱波速度、密度的第一關系模型，并根據所述目標區塊的數字巖芯數據確定孔隙度與縱波速度、密度的第二關系模型；

將所述第一關系模型和所述第二關系模型的加權和作為第三關系模型；

基于所述測井數據獲取所述目標區塊的第一孔隙結構因子，并基于所述數字巖芯數據獲取所述目標區塊的第二孔隙結構因子；

將所述第一孔隙結構因子和所述第二孔隙結構因子的加權和作為第三孔隙結構因子；

將所述縱波速度數據體和所述密度數據體輸入所述第三關系模型，獲得所述目標區塊的孔隙度數據體；

將所述縱波速度數據體與所述孔隙度數據體進行交會，獲得所述目標區塊的孔隙結構數據體；

其中，所述基于所述測井數據獲取所述目標區塊的第一孔隙結構因子，包括：

根據公式確定所述目標區塊的第一孔隙結構因子；

其中，ρ₁為目標區塊基于測井數據的密度，ρ_m1為目標區塊基于測井數據的孔隙相中固體相礦物骨架的密度，ρ_f1為目標區塊基于測井數據的孔隙相中充填的流體或氣體的密度，φ₁為目標區塊基于測井數據的孔隙度，K₁為目標區塊基于測井數據的體積模量，V_P1為目標區塊基于測井數據在指定頻段下隨圍壓變化的縱波速度，V_S1為目標區塊基于測井數據在指定頻段下隨圍壓變化的橫波速度，μ₁為目標區塊基于測井數據的剪切模量，F_k1為第一中間變量，K_m1為目標區塊基于測井數據的孔隙相中固體相礦物骨架的體積模量，K_f1為目標區塊基于測井數據的孔隙相中充填的流體或氣體的體積模量，f₁為第二中間變量，γ₁為目標區塊基于測井數據的孔隙結構因子；

所述基于所述數字巖芯數據獲取所述目標區塊的第二孔隙結構因子，包括：

根據公式確定所述目標區塊的第二孔隙結構因子；

其中，ρ₂為目標區塊基于數字巖芯的密度，ρ_m2為目標區塊基于數字巖芯的孔隙相中固體相礦物骨架的密度，ρ_f2為目標區塊基于數字巖芯的孔隙相中充填的流體或氣體的密度，φ₂為目標區塊基于數字巖芯的孔隙度，K₂為目標區塊基于數字巖芯的體積模量，V_P2為目標區塊基于數字巖芯在指定頻段下隨圍壓變化的縱波速度，V_S2為目標區塊基于數字巖芯在指定頻段下隨圍壓變化的橫波速度，μ₂為目標區塊基于數字巖芯的剪切模量，F_k2為第三中間變量，K_m2為目標區塊基于數字巖芯的孔隙相中固體相礦物骨架的體積模量，K_f2為目標區塊基于數字巖芯的孔隙相中充填的流體或氣體的體積模量，f₂為第四中間變量，γ₂為目標區塊基于數字巖芯的孔隙結構因子。

2.如權利要求1所述的碳酸鹽巖儲層孔隙結構預測方法，其特征在于，所述根據所述目標區塊的測井數據確定孔隙度與縱波速度、密度的第一關系模型，包括：

獲取所述目標區塊的測井數據；所述測井數據包括孔隙度、縱波速度和密度；

基于預設的機器學習模型訓練所述測井數據，以獲取孔隙度與縱波速度、密度的第一關系模型。

3.如權利要求1所述的碳酸鹽巖儲層孔隙結構預測方法，其特征在于，所述根據所述目標區塊的數字巖芯數據確定孔隙度與縱波速度、密度的第二關系模型，包括：

獲取所述目標區塊的指定巖樣的巖石物理參數；

基于所述巖石物理參數及所述指定巖樣的CT掃描圖像構建數字巖芯；

讀取所述數字巖芯的指定數據；所述指定數據包括孔隙度、縱波速度和密度；

基于預設的機器學習模型訓練所述指定數據，以獲取孔隙度與縱波速度、密度的第二關系模型。