(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202111486586.3 (22)申请日 2021.12.07 (65)同一申请的已公布的文献号 申请公布号 CN 114239431 A (43)申请公布日 2022.03.25 (73)专利权人 中国石油大 学(北京) 地址 102249 北京市昌平区府学路18号 (72)发明人 贾品　郭红鑫　曹冲　程林松　 王鹏　 (74)专利代理 机构 北京三友知识产权代理有限 公司 11127 专利代理师 谢层层　徐焕 (51)Int.Cl. G06F 30/28(2020.01) G06F 113/08(2020.01)G06F 119/14(2020.01) (56)对比文件 CN 108536982 A,2018.09.14 CN 111062129 A,2020.04.24 CN 111104766 A,2020.0 5.05 姜永等.裂缝性油藏注水井 压降试井解释方 法及应用. 《断块油气田》 .2017,(第04期), 审查员 慈丽雁 (54)发明名称 裂缝发育储层水驱油的模拟方法、 装置及设 备 (57)摘要 本说明书实施例公开了一种裂缝发育储层 水驱油的模拟方法、 装置及设备， 基于裂缝性储 层中不同产状裂缝岩心注水驱替流线特征， 将基 质划分若干线性子区域， 将水驱油二维问题转化 为一维线性子区域耦合求解问题， 将研究区域划 分为注入区、 裂缝区以及采出区， 分别创建对应 的线性流模型， 依次对各个线性流模型进行求 解， 求解注入端基质区线性驱替解析解， 裂缝线 性导流数值解， 出口端基质驱替解析解， 并耦合 基质与裂缝系统线性流， 将水驱油二维问题转化 为一维线性子区域耦合求解问题， 获得带缝岩心 水驱油三区线性流动半解析解， 进而获得裂缝 发 育储层压力分布和饱和度分布， 实现了对裂缝 发 育储层水驱油过程的准确模拟。 权利要求书3页 说明书16页 附图4页 CN 114239431 B 2022.08.26 CN 114239431 B 1.一种裂缝发育储层水驱油的模拟方法， 其特 征在于， 所述方法包括： 根据待处理裂缝发育储层中的裂缝产状分布将所述待处理裂缝发育储层对应的双重 介质模型中的基质离散划分为多个一维线性子区域， 其中， 每个一维线性子区域中包括有 裂缝微元； 将每个一维线性子区域划分为注入区、 裂缝区、 采出区， 根据 所述注入区、 所述裂缝 区、 所述采出区对应的流动特征， 分别建立各个一维线性子区域对应的三线性流模型， 所述三 线性流模型中包括： 注入区线性 流模型、 裂缝区线性 流模型、 采出区线性 流模型； 依次对各个一维线性子区域对应的三线性流模型中的注入区线性流模型以及采出区 线性流模型进行求解， 获得各个一维线性子区域对应的注入区的油水流量、 注入区压力 分 布、 注入区饱和度分布和采出区的油水流 量、 采出区压力分布、 采出区饱和度分布； 基于各个一维线性子区域对应的注入区的油水流量和采出区的油水流量， 依次对各个 一维线性子区域对应的三线性流模型中的裂缝区线性流模型进行求解， 获得各个一 维线性 子区域对应的裂缝区的裂缝压力和裂缝含水饱和度； 根据获得的裂缝压力、 所述裂缝含水饱和度、 所述注入区压力分布、 所述注入区饱和度 分布、 所述采出区压力分布、 所述采出区饱和度分布， 对所述待处理裂缝发育储层的饱和度 场和压力场进行模拟。 2.如权利要求1所述的方法， 其特征在于， 所述根据 所述注入区、 所述裂缝区、 所述采出 区对应的流动特 征， 分别建立各个一维线性子区域对应的三线性 流模型， 包括： 按照下述公式建立各个一维线性子区域对应的三线性 流模型中的注入区线性 流模型： 初始条件： pm|t＝0＝pi,So|t＝0＝Soi 内边界条件： 外边界条件： 其中， kro为油相相对渗透率， krw为水相相对渗透率， μo为油相粘度， μw为水相粘度， x为 流体线性流动方向， Bo为原油体积系数， Bw为水体积系数， pm为基质压力， pi表示初始压力， qo,m为单位体积下一 维线性子区域间窜流的油量， qw,m为单位体积下一 维线性子区域间窜流 的水量， φ为孔隙度， km为基质绝对渗透率， So为含油饱和度， Sw为含水饱和度， t为时间， Soi 为原始含油饱和度， qw为单个一维线性子区域的注水量， LL为裂缝位置， pf为裂缝的压力。 3.如权利要求1所述的方法， 其特征在于， 所述根据 所述注入区、 所述裂缝区、 所述采出 区对应的流动特 征， 分别建立各个一维线性子区域对应的三线性 流模型， 包括： 按照下述公式建立各个一维线性子区域对应的三线性 流模型中的裂缝区线性 流模型：权　利　要　求　书 1/3 页 2 CN 114239431 B 2初始条件： 内边界条件： 外边界条件： 其中， kf为裂缝绝对渗透率， kro为油相相对渗透率， krw为水相相对渗透率， μo为油相粘 度， μw为水相粘度， ε为裂缝延伸方 向， Bo为原油体积系数， Bw为水体积系数， qo,L为单位体积 下注入区供给的油量， qo,R为单位体积下采 出区供给的油量， 取负值， qw,L为单位体积下注入 区供给的水量， qw,R为单位体积下采出区供给的水量， 取负值， So为含油饱和度， Sw为含水饱 和度， t为时间， Soi为原始含油饱和度， pf为裂缝的压力， pi表示初始压力， Lf表示裂缝的总长 度， 裂缝的内外边界为封闭边界。 4.如权利要求1所述的方法， 其特征在于， 所述根据 所述注入区、 所述裂缝区、 所述采出 区对应的流动特 征， 分别建立各个一维线性子区域对应的三线性 流模型， 包括： 按照下述公式建立各个一维线性子区域对应的三线性 流模型中的采出区线性 流模型： 初始条件： 内边界条件： 外边界条件： 其中， kro为油相相对渗透率， krw为水相相对渗透率， μo为油相粘度， μw为水相粘度， x为 流体线性流动方向， Bo为原油体积系数， Bw为水体积系数， pm为基质压力， pi表示初始压力， qo,m为单位体积下一 维线性子区域间窜流的油量， qw,m为单位体积下一 维线性子区域间窜流 的水量， φ为孔隙度， km为基质绝对渗透率， So为含油饱和度， Sw为含水饱和度， t为时间， Soi 为原始含油饱和度， qw为单个一维线性子区域的注水量， LL为裂缝位置， pw为出口端压力。权　利　要　求　书 2/3 页 3 CN 114239431 B 3