CN118981973B 一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法 （中國(guó)地質(zhì)大學(xué)（北京））

HeterogeneousPore-ThroatNetwork一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快本發(fā)明公開了一種多孔介質(zhì)微觀油水運(yùn)動(dòng)基于最短算法計(jì)算連通入口和出口的路徑及帶和位置預(yù)測(cè)模型得到模型中各類型剩余油分布2制作不同物性特征的孔喉網(wǎng)絡(luò)物理模型，在將所對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行圖像二值化處理，將包含孔喉結(jié)構(gòu)的孔喉網(wǎng)絡(luò)物理采用機(jī)器學(xué)習(xí)方法確定調(diào)節(jié)系數(shù)入、剩余油類型飽和度預(yù)測(cè)模型Sroi喉位置預(yù)測(cè)模型Z，基于調(diào)節(jié)系數(shù)入建立喉道權(quán)重模型ρ,通過(guò)喉道權(quán)重模型ρ建立孔喉網(wǎng)絡(luò)所述調(diào)節(jié)系數(shù)入的確定方法為：對(duì)帶權(quán)有向圖應(yīng)用狄克斯特拉最短路徑算法，將帶權(quán)路徑孔喉順序和根據(jù)路徑集合中孔喉編號(hào)和對(duì)應(yīng)的孔喉屬性參數(shù)計(jì)算各條路徑上孔喉的覆蓋面積與設(shè)定一個(gè)帶權(quán)路徑截止值Lthreshold，得到[L0,Lthreshold]區(qū)在所述水流路徑上重構(gòu)剩余油類型，并且在所述水流路徑范2.根據(jù)權(quán)利要求1所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特3.根據(jù)權(quán)利要求2所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特34.根據(jù)權(quán)利要求3所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特;5.根據(jù)權(quán)利要求1所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特;計(jì)算不同實(shí)驗(yàn)的穩(wěn)定狀態(tài)下不同采收率時(shí)的各類型剩余油飽和度S6.根據(jù)權(quán)利要求1所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特;7.根據(jù)權(quán)利要求6所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特根據(jù)不同實(shí)驗(yàn)分析各類型剩余油賦存在孔隙和根據(jù)剩余油類型孔喉位置預(yù)測(cè)模型Z計(jì)算模型驅(qū)替波及范圍內(nèi)所有喉道和孔隙的判定參數(shù)，分別記為Rcoi和pupo?,并從大到小排序;以Rcoi和pupo的前列大值與Rcoro和puporo誤差最小為目標(biāo)函數(shù)進(jìn)行機(jī)器學(xué)習(xí)方法48.根據(jù)權(quán)利要求7所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特將水流路徑之內(nèi)的剩余油判定為柱狀剩余油和多孔狀9.根據(jù)權(quán)利要求8所述的一種多孔介質(zhì)微觀油水運(yùn)動(dòng)及水流路徑快速模擬方法，其特在所述水流路徑范圍圖上添加柱狀剩余油和多孔狀剩余所述柱狀剩余油和多孔狀剩余油的添加數(shù)量依據(jù)向量Rway和向量Pway，56余油類型飽和度預(yù)測(cè)模型Sroi以及剩余油類型孔喉位置[0017]在進(jìn)行孔喉連通性計(jì)算時(shí)建立鄰接矩陣，若編號(hào)為i的孔體和編號(hào)為j的孔體相[0020]水流路徑的權(quán)重與水流路徑的流動(dòng)阻力成正比，則所述喉道權(quán)重模型ρ計(jì)算方式長(zhǎng)度L0[0025]根據(jù)路徑集合中孔喉編號(hào)和對(duì)應(yīng)的孔喉屬性參數(shù)計(jì)算各條路徑上孔喉的覆蓋面[0026]設(shè)定一個(gè)帶權(quán)路徑截止值Lthreshold，得到[L0,Lthreshold]區(qū)間內(nèi)所有路徑孔隙和喉74+b3+c+e；[0040]進(jìn)一步地，建立水流路徑波及范圍內(nèi)剩余油類型孔喉位置預(yù)測(cè)模型Z的具體方法[0044]θ為喉道中軸線和驅(qū)替壓力方向之間的夾角，τ為喉道的迂曲度，r為相對(duì)喉道半[0051]根據(jù)剩余油類型孔喉位置預(yù)測(cè)模型Z計(jì)算模型驅(qū)替波及范圍內(nèi)所有孔隙和喉道的[0052]以PCO_L和PMPO_L的前列大值與PCO_ro和PMPO_ro誤差最小為目標(biāo)8片狀剩余油，并進(jìn)一步通過(guò)R_CO和P_MPO在水流路徑范圍內(nèi)添加柱狀剩余油和多孔狀剩余[0072]圖4為本發(fā)明實(shí)施例提供微觀孔喉網(wǎng)絡(luò)模型孔喉屬性參數(shù)一孔隙半徑和喉道半徑[0074]圖6為本發(fā)明實(shí)施例提供微觀孔喉網(wǎng)絡(luò)數(shù)字模型孔隙喉道識(shí)別、分割與編號(hào)后圖9[0076]圖8為本發(fā)明實(shí)施例基于狄克斯特拉算法計(jì)算的最短水流路徑及覆蓋孔喉示意[0078]圖10為本發(fā)明實(shí)施例孔喉網(wǎng)絡(luò)內(nèi)水流路徑及剩余油類型快速預(yù)測(cè)重構(gòu)結(jié)果示意[0094]第五步，基于狄克斯特拉(Dijkstra)等最短路徑算法計(jì)算入[0109]所有剩余油類型飽和度與非油相飽和度加和為1，剩余油類型可在此原則下進(jìn)行[0126]喉道的權(quán)重模型定義為ρ,權(quán)重ρ由影響流動(dòng)阻力的孔喉屬性參數(shù)計(jì)算獲[0132]根據(jù)路徑上孔喉編號(hào)和對(duì)應(yīng)孔喉屬性參數(shù)計(jì)算各條路徑上孔喉的覆蓋面積(體孔隙和喉道編號(hào)、路徑覆蓋面積(體積)與孔喉總面積(體積)的百分比Sthreshold；此時(shí)Lthreshold與Sthreshold存在一ρ;[0143]假設(shè)各類型剩余油飽和度與采收率的函數(shù)關(guān)系為多項(xiàng)式關(guān)系，為表述和計(jì)算方3+c[0148]進(jìn)一步地，建立水流路徑波及范圍內(nèi)剩余油類型孔喉位置預(yù)測(cè)模型Z的具體方法[0152]θ為喉道中軸線和驅(qū)替壓力方向之間的夾角，τ為喉道的迂曲度，r為等效喉道半[0157]根據(jù)預(yù)測(cè)模型Z計(jì)算模型驅(qū)替相波及范圍內(nèi)所有孔喉的PCO_L和PMPO_L，從大到小排在目標(biāo)采收率α?xí)r的各類剩余油飽和度，根據(jù)水流路徑之外的剩余油類型飽和度確定水流[0160]確定帶權(quán)路徑長(zhǎng)度閥值Lthreshold，使得最短帶權(quán)路徑長(zhǎng)度到帶權(quán)路徑長(zhǎng)度閥值區(qū)間[L0,Lthreshold]內(nèi)路徑所通過(guò)的孔隙和喉道面積(體積)百分比等于水流路徑波及體積]內(nèi)路徑所通過(guò)的所有孔隙和喉道即為設(shè)置目標(biāo)采收率α下的[0170]其中，在RCO最大的喉道內(nèi)添加柱狀剩余油，在PMPO最大的孔[0171]依據(jù)向量Rway和向量Pway，從大到小向?qū)?yīng)孔喉中填充柱狀剩余油和多孔狀剩余實(shí)驗(yàn)證明，傳統(tǒng)方法適用于孔喉數(shù)量較少的小尺寸模型，而且計(jì)算速度慢，如Navier_方案具有更快的計(jì)算速度，可進(jìn)一步應(yīng)用于孔喉數(shù)量多的大尺寸模型，如在AMDRyzen種修改或等