(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111323681.1 (22)申请日 2021.11.10 (71)申请人 长江大学 地址 434023 湖北省荆州市南环路1号 (72)发明人 饶翔　赵辉　徐云峰　刘怡娜　 (74)专利代理 机构 武汉谦源知识产权代理事务 所(普通合伙) 42251 代理人 王力 (51)Int.Cl. G06F 30/28(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 基于迎风GFDM的多孔介质油水两相流计算 方法 (57)摘要 本发明公开了一种基于迎风GFDM的多孔介 质油水两相流计算方法， 包括以下步骤： S1、 建立 多孔介质油水两相流模型； S2、 基于迎风GFDM对 所述多孔介质油水两相流模型进行离散， 得到所 述多孔介质油水两相流模型的全隐式离散格式； S3、 利用基于牛顿迭代和自动微分的非线性求解 器对所述多孔介质油水两相流模型的全 隐式离 散格式进行求解， 实现压力和含 水饱和度的高精 度无网格计算。 本发明提供一种基于迎风 GFDM的 多孔介质油水两相流计算方法， 为形成不依赖网 格的多孔介质流动通用型无网格数值模拟器提 供了参考基础。 权利要求书2页 说明书9页 附图3页 CN 114021502 A 2022.02.08 CN 114021502 A 1.一种基于 迎风GFDM的多孔介质油水两相流计算方法， 其特 征在于， 包括以下步骤： S1、 建立多孔介质油水两相流模型； S2、 基于迎风GFDM对所述多孔介质油水两相流模型进行离散， 得到所述多孔介质油水 两相流模型的全隐式离 散格式； S3、 利用基于牛顿迭代和自动微分的非线性求解器对所述多孔介质油水两相流模型的 全隐式离 散格式进行求 解， 实现压力和含水饱和度的高精度无网格 计算。 2.如权利要求1所述的基于迎风GFDM的多孔介质油水两相流计算方法， 其特征在于， 所 述多孔介质油水两相流模型包括油相物质守恒方程、 水相物质守恒方程和辅助方程， 其中 多孔介质中油水两相流互不相溶； 其中， 所述油相物质守恒方程 为： 上式中： k是渗透率， mD； kro＝kro(Sw)是与含水饱和度Sw有关的油相相对渗透率； μo是油 相粘度， mPa ·s； po是油相压力， MPa； qo是油相源汇项， 1/天； t是时间， 天； So是油相饱和度； φ＝φ(p)是与压力有关的储层孔隙度； 所述水相物质守恒方程 为： 上式中： krw＝krw(Sw)是与含水饱和度有关的水相相对 渗透率； μw是水相粘度， mPa ·s； pw 是水相压力， MPa； qo是水相源汇项， 1/天； 所述辅助方程 为： So＝1‑Sw； pw＝po‑Pc(Sw)； φ＝φ(po)＝φ0+Cr(po‑po0)； krw＝krw(Sw)； kro＝kro(Sw)；                                                         (3) 上式中： po为油相压力； Sw为含水饱和度； Ct＝Ct(Sw)是综合压缩系数， 1/MPa； po0是储层 原始油相压力， φ0是po0对应的孔隙度。 3.如权利 要求1所述的基于迎风GFDM的多孔介质油水两相流计算方法， 其特征在于， S2 中得到多孔介质油水两相流计算方法如下： 2.1、 基于GFDM的油相物质守恒方程的离 散格式为： 权　利　要　求　书 1/2 页 2 CN 114021502 A 2上式中： Pc,(i,j)t+Δt＝Pc(Sw,(i,j)t+Δt)； Pc,it+Δt＝Pc(Sw,it+Δt)； φit+Δt＝φ(po,it+Δt)； φit＝ φ(po,it)； 2.2、 基于GFDM的水相物质守恒方程的离 散格式为： 上式中： Pc,(i,j)t+Δt＝Pc(Sw,(i,j)t+Δt)； Pc,it+Δt＝Pc(Sw,it+Δt)； 2.3、 油相的相对渗透率的迎风格式如下： 水相的相对渗透率的迎风格式如下： 权　利　要　求　书 2/2 页 3 CN 114021502 A 3