CN118981973B 一种多孔介质微观油水运动及水流路径快速模拟方法 （中国地质大学（北京））

HeterogeneousPore-ThroatNetwork一种多孔介质微观油水运动及水流路径快本发明公开了一种多孔介质微观油水运动基于最短算法计算连通入口和出口的路径及带和位置预测模型得到模型中各类型剩余油分布2制作不同物性特征的孔喉网络物理模型，在将所对实验数据进行图像二值化处理，将包含孔喉结构的孔喉网络物理采用机器学习方法确定调节系数入、剩余油类型饱和度预测模型Sroi喉位置预测模型Z，基于调节系数入建立喉道权重模型ρ,通过喉道权重模型ρ建立孔喉网络所述调节系数入的确定方法为：对带权有向图应用狄克斯特拉最短路径算法，将带权路径孔喉顺序和根据路径集合中孔喉编号和对应的孔喉属性参数计算各条路径上孔喉的覆盖面积与设定一个带权路径截止值Lthreshold，得到[L0,Lthreshold]区在所述水流路径上重构剩余油类型，并且在所述水流路径范2.根据权利要求1所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特3.根据权利要求2所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特34.根据权利要求3所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特;5.根据权利要求1所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特;计算不同实验的稳定状态下不同采收率时的各类型剩余油饱和度S6.根据权利要求1所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特;7.根据权利要求6所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特根据不同实验分析各类型剩余油赋存在孔隙和根据剩余油类型孔喉位置预测模型Z计算模型驱替波及范围内所有喉道和孔隙的判定参数，分别记为Rcoi和pupo⃞,并从大到小排序;以Rcoi和pupo的前列大值与Rcoro和puporo误差最小为目标函数进行机器学习方法48.根据权利要求7所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特将水流路径之内的剩余油判定为柱状剩余油和多孔状9.根据权利要求8所述的一种多孔介质微观油水运动及水流路径快速模拟方法，其特在所述水流路径范围图上添加柱状剩余油和多孔状剩余所述柱状剩余油和多孔状剩余油的添加数量依据向量Rway和向量Pway，56余油类型饱和度预测模型Sroi以及剩余油类型孔喉位置[0017]在进行孔喉连通性计算时建立邻接矩阵，若编号为i的孔体和编号为j的孔体相[0020]水流路径的权重与水流路径的流动阻力成正比，则所述喉道权重模型ρ计算方式长度L0[0025]根据路径集合中孔喉编号和对应的孔喉属性参数计算各条路径上孔喉的覆盖面[0026]设定一个带权路径截止值Lthreshold，得到[L0,Lthreshold]区间内所有路径孔隙和喉74+b3+c+e；[0040]进一步地，建立水流路径波及范围内剩余油类型孔喉位置预测模型Z的具体方法[0044]θ为喉道中轴线和驱替压力方向之间的夹角，τ为喉道的迂曲度，r为相对喉道半[0051]根据剩余油类型孔喉位置预测模型Z计算模型驱替波及范围内所有孔隙和喉道的[0052]以PCO_L和PMPO_L的前列大值与PCO_ro和PMPO_ro误差最小为目标8片状剩余油，并进一步通过R_CO和P_MPO在水流路径范围内添加柱状剩余油和多孔状剩余[0072]图4为本发明实施例提供微观孔喉网络模型孔喉属性参数一孔隙半径和喉道半径[0074]图6为本发明实施例提供微观孔喉网络数字模型孔隙喉道识别、分割与编号后图9[0076]图8为本发明实施例基于狄克斯特拉算法计算的最短水流路径及覆盖孔喉示意[0078]图10为本发明实施例孔喉网络内水流路径及剩余油类型快速预测重构结果示意[0094]第五步，基于狄克斯特拉(Dijkstra)等最短路径算法计算入[0109]所有剩余油类型饱和度与非油相饱和度加和为1，剩余油类型可在此原则下进行[0126]喉道的权重模型定义为ρ,权重ρ由影响流动阻力的孔喉属性参数计算获[0132]根据路径上孔喉编号和对应孔喉属性参数计算各条路径上孔喉的覆盖面积(体孔隙和喉道编号、路径覆盖面积(体积)与孔喉总面积(体积)的百分比Sthreshold；此时Lthreshold与Sthreshold存在一ρ;[0143]假设各类型剩余油饱和度与采收率的函数关系为多项式关系，为表述和计算方3+c[0148]进一步地，建立水流路径波及范围内剩余油类型孔喉位置预测模型Z的具体方法[0152]θ为喉道中轴线和驱替压力方向之间的夹角，τ为喉道的迂曲度，r为等效喉道半[0157]根据预测模型Z计算模型驱替相波及范围内所有孔喉的PCO_L和PMPO_L，从大到小排在目标采收率α时的各类剩余油饱和度，根据水流路径之外的剩余油类型饱和度确定水流[0160]确定带权路径长度阀值Lthreshold，使得最短带权路径长度到带权路径长度阀值区间[L0,Lthreshold]内路径所通过的孔隙和喉道面积(体积)百分比等于水流路径波及体积]内路径所通过的所有孔隙和喉道即为设置目标采收率α下的[0170]其中，在RCO最大的喉道内添加柱状剩余油，在PMPO最大的孔[0171]依据向量Rway和向量Pway，从大到小向对应孔喉中填充柱状剩余油和多孔状剩余实验证明，传统方法适用于孔喉数量较少的小尺寸模型，而且计算速度慢，如Navier_方案具有更快的计算速度，可进一步应用于孔喉数量多的大尺寸模型，如在AMDRyzen种修改或等