[发明专利]一种考虑聚合物沉淀及吸附的注聚井伤害计算方法在审
申请号: | 202310515315.9 | 申请日: | 2023-05-08 |
公开(公告)号: | CN116629105A | 公开(公告)日: | 2023-08-22 |
发明(设计)人: | 刘平礼;黄城熹;杜娟;陈祥;李骏;刘金明;王铭;王冠;李策;娄凤成 | 申请(专利权)人: | 西南石油大学 |
主分类号: | G06F30/27 | 分类号: | G06F30/27;G06F30/28;E21B43/16;G06F111/10 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 610500 四*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 考虑 聚合物 沉淀 吸附 注聚井 伤害 计算方法 | ||
1.一种考虑聚合物沉淀及吸附的注聚井伤害计算方法,其特征在于,包括如下步骤:
S1、取待模拟的目标储层岩心,通过对岩心进行室内聚合物驱伤害评价实验,得到岩心孔渗及尺寸、聚合物性质及驱替参数:初始渗透率K0和孔隙度伤害后的渗透率K1和孔隙度室内岩心伤害时间t1、稳态下的注入流量Q1、聚合物粘度μ1、岩心长度L和半径r;
S2、实验岩心r方向等比例划分、θ和z方向等间距划分网格,实验时间等步长划分,采用高斯分布模型、Garman-Kozeny半经验公式分别构建初始孔隙度和渗透率场:
式中:f(x)为概率密度函数,无量纲;
x为孔隙度随机变量,f;
s为孔隙度标准差,无量纲;
为平均孔隙度,f;
K为渗透率,m2;
为孔隙度,f;
β为参数,由实验室测得数据,无量纲;
S3、依次求解压力场和速度场;
初始条件:
Pr,θ,z=0 (4)
边界条件:
式中:P1为岩心出口端压力,Pa;
式中:u为r方向流体流速,m/s;
V为θ方向流体流速,m/s;
ξ为z方向流体流速,m/s;
S4、求解聚合物浓度场;
S5、考虑聚合物沉积和吸附伤害,求解聚合物驱替过程中岩心孔隙度和渗透率场;
S6、若计算时间等于室内岩心伤害时间t1,则计算结束,得出岩心的平均渗透率和平均孔隙度否则,继续执行步骤S3~S5;计算目标函数
S7、采用NSGA-Ⅱ多目标遗传算法优化目标函数,优化目标函数A和B;
S8、选取目标函数A和B所对应的α和β,参照步骤S2~S6对注聚井伤害数值模拟计算。
2.如权利要求1所述的考虑聚合物沉淀及吸附的注聚井伤害计算方法,其特征在于,步骤S4求解聚合物浓度场,具体过程如下:
式中:为聚合物吸附量,ppm;
σ为聚合物沉积量,ppm;
C为聚合物浓度,ppm;
U为聚合物流速,m/s;
D为聚合物扩散系数,m2/s;
为梯度算子;
聚合物沉积量σ:
聚合物吸附量
式中:λ0为聚合物过滤系数,m-1;
b为模型参数,无量纲;
为最大无量纲吸附聚合物浓度,ppm;
初始条件:
Cr,θ,z=0 (10)
边界条件:
式中:C1为岩心端面入口处的聚合物浓度,ppm;
C2为岩心端面出口处的聚合物浓度,ppm。
3.如权利要求1所述的考虑聚合物沉淀及吸附的注聚井伤害计算方法,其特征在于,步骤S5考虑聚合物沉积和吸附伤害,求解聚合物驱替过程中岩心孔隙度和渗透率场,具体过程如下:
式中:α为聚合物沉淀参数,ppm-1;
β为聚合物吸附参数,ppm-1。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西南石油大学,未经西南石油大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202310515315.9/1.html,转载请声明来源钻瓜专利网。