(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111577672.5 (22)申请日 2021.12.2 2 (71)申请人 长江大学 地址 434023 湖北省荆州市南环路1号 (72)发明人 赵辉　周玉辉　盛广龙　孟凡坤　 饶翔　钟珣　刘伟　雷昇　 (74)专利代理 机构 武汉谦源知识产权代理事务 所(普通合伙) 42251 专利代理师 王力 (51)Int.Cl. E21B 43/20(2006.01) E21B 33/13(2006.01) G06F 30/28(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 基于微凝胶驱的油藏动态预测方法、 装置、 介质和终端 (57)摘要 本发明公开了一种基于微凝胶驱的油藏动 态预测方法、 装置、 介质和终端， 方法包括以下步 骤： 将油藏注采系统简化为井与井 之间的连通单 元； 以连通单元为对象建立物质 平衡方程， 并生 成井间流量和井点含水饱和度； 在 注水井中注入 预设浓度的微凝胶， 并建立微凝胶驱的浓度平衡 方程； 计算目标时刻注水井在各个连通单元方向 上的微凝胶颗粒相浓度； 调整水相渗透率下降系 数， 并计算目标时刻的井点含水饱和度， 以预测 加入微凝胶后目标区块的生产动态变化。 本发明 将微凝胶的封堵特性和井间连通性思想相结合， 对微凝胶颗粒相进入孔道后的浓度分布、 封堵情 况及驱油效果进行动态预测， 从而对现场开发方 案进行指导， 提高非均质油藏的采油量， 增加经 济效益。 权利要求书3页 说明书11页 附图5页 CN 115538997 A 2022.12.30 CN 115538997 A 1.一种基于微凝胶驱的油藏动态预测方法， 其特 征在于， 包括以下步骤： 步骤1， 建立用于模拟井间油水动态的初始连通性模型， 所述初始连通性模型将油藏注 采系统简化为井与井之 间的连通单元， 所述连通单元由传导率和连通体积两个特征参数表 征； 步骤2， 以所述连通单元为对象建立物质平衡方程， 并对所述物质平衡方程进行定压求 解生成井间流 量和井点含水饱和度； 步骤3， 在注水井中注入预设浓度的微凝胶， 并根据所述井间流量和所述井点含水饱和 度建立微凝胶驱的浓度平衡方程； 步骤4， 对所述浓度平衡方程进行求解， 计算注入微凝胶后目标时刻所述注水井在各个 连通单元方向上的微凝胶颗粒相浓度； 步骤5， 根据目标时刻连通单元的微凝胶颗粒相浓度调 整对应的水相渗透率下降系数， 并计算目标时刻的井点含水饱和度， 以预测 加入微凝胶后目标区块的生产动态变化。 2.根据权利要求1所述基于微凝胶驱的油藏动态预测方法， 其特征在于， 所述以连通单 元为对象建立物质平衡方程， 并对物质平衡方程进行定压求解生成井间流量和井点含水饱 和度， 具体为： S201， 以所述连通单 元为对象建立物质平衡方程 为： 其中， i、 j表示井序号； N为总井数量； n为当前时刻， n ‑1为前一时刻； 和 分别表示第 i井和第j井在n时刻泄油区内的平均油藏压力； Δt为时间间隔； 为第n时刻第i井的流速； Ct为油藏综合压缩系数； Tij表示第i井和第j井的平均传导率； Vpi为第i口井泄油区的连通体 积； S202， 采用定压方法对所述物质平衡方程进行隐式差分求解， 生成连通单元之间的井 间流量为： 其中， n为当前时刻， 表示第n时刻的第i井和第j井的井间连通单 元的流速； S203， 基于关停井及转注造成的井底液流反转， 利用插值方法 反算井点含水饱和度为： 其中， n为当前时刻， n ′为上一时刻， fw为含水率， f ′w(Swijk)为第k层从第j井追踪到第i 井处的含水率导数， f ′w(Swjk)为第j井在第k层的含水率导数， Swik为第i井在第k层的含水饱 和度， Swjk为第j井在第k层的含水饱和 度， Swijk为第k层i井与j井间连通单元的含水饱和 度； 代表从n′时刻到n时刻从j井 流向i井的反向无因次累计流量， 表示从0时刻直到 n′时 刻i井流向j井的无因次累计流量， 代表从0时刻直到n时刻i井流向j井的无因次累计流 量。 3.根据权利要求1或2所述基于微凝胶驱的油藏动态预测方法， 其特征在于， 所述在注 水井中注入预设浓度的微凝胶， 并根据井间流量和井点含水饱和度建立微凝胶驱的浓度 平权　利　要　求　书 1/3 页 2 CN 115538997 A 2衡方程为： 其中， Nw为注采井数， n为当前时刻， ypijk表示i井与j井在第k层连通单元 的微凝胶颗粒 相浓度； Swijk是i井与j井在第k层连通单元的含水饱和度； 和 分别为上一时刻i井与j 井在第k层的连通单 元流入与流出的流 量； Vijk为第k层i井与j井间的连通体积。 4.根据权利要求3所述基于微凝胶驱的油藏动态预测方法， 其特征在于， 所述对浓度平 衡方程进行求 解， 计算注入微凝胶后注水井在各个连通单 元方向上微凝胶的浓度， 具体为： S401， 采用油 藏的历史油水产出数据对所述井点含水饱和度进行拟合反演， 调整所述 初始连通 性模型的特 征参数， 并生成优化后的连通 性模型； S402， 根据所述优化后的连通性模型计算所述注水井向各方向上油井的注水劈分系 数； S403， 调用预设微凝胶浓度分布表达式， 并基于所述注水劈分系数计算初始时刻所述 注水井在各个连通单 元方向上的微凝胶颗粒相浓度； S404， 根据初始时刻各个连通单元方向上的微凝胶颗粒相浓度对所述浓度平衡方程进 行求解， 生成目标时刻所述注水井在各个连通单 元方向上的微凝胶颗粒相浓度。 5.根据权利要求4所述基于微凝胶驱的油藏动态预测方法， 其特征在于， 所述预设微凝 胶浓度分布 表达式为： 其中， yi＝0时对应的xi被定义为临界分流量， 记作 xi表示连通单元i的分流量， yi表 示连通单 元i的微凝胶颗粒相浓度， A为与注水井连通的连通单 元总数,a为颗粒 粒径参数。 6.根据权利要求5所述基于微凝胶驱的油藏动态预测方法， 其特征在于， 所述根据目标 时刻连通单元的微凝胶颗粒相浓度调整对应的水相渗透率下降系数， 并计算目标时刻的井 点含水饱和度的预设公式为： 其中， yi为连通单元i在目标时刻的微凝胶颗 粒相浓度， r、 s、 t为预设值， Rmax为最大渗透 率下降系数， kro、 krw分别为地层原油、 地层水的相对渗透率值， μo、 μw分别为第k层地层原油、 地层水的黏度值,Ck是渗透率相关系数,Rk是渗透率下降系数, μl是黏浓相关系数,fw是含水权　利　要　求　书 2/3 页 3 CN 115538997 A 3