[发明专利]考虑变启动压力梯度的水驱稠油数值模拟方法在审
申请号: | 201910962175.3 | 申请日: | 2019-10-11 |
公开(公告)号: | CN110685653A | 公开(公告)日: | 2020-01-14 |
发明(设计)人: | 周文胜;耿艳宏;刘凡;姜瑞忠;张春光;王凯;刘晨;李珂 | 申请(专利权)人: | 中海石油(中国)有限公司;中海石油(中国)有限公司北京研究中心 |
主分类号: | E21B43/20 | 分类号: | E21B43/20;E21B47/00;G06Q50/02 |
代理公司: | 11245 北京纪凯知识产权代理有限公司 | 代理人: | 李晓红 |
地址: | 100010 北*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 稠油区块 启动压力 数值模拟模型 饱和度场 稠油 饱和度 关系模型 模拟模型 收敛条件 数学模型 数值模拟 物性参数 物质平衡 初始化 梯度场 压力场 迭代 流度 求解 渗流 隐式 油藏 检查 输出 重复 | ||
1.一种考虑变启动压力梯度的水驱稠油数值模拟方法,其特征在于,其包括以下步骤:
S1、确定水驱稠油区块的油藏物性参数,形成水驱稠油区块的模拟模型;
S2、建立水驱稠油变启动压力梯度数值模拟模型;
1)获得油相启动压力梯度和流度的关系模型
2)将步骤1)中的油相启动压力梯度和流度的关系式代入水驱稠油区块的渗流运动方程中,并结合连续方程和状态方程建立水驱稠油区块的渗流数学模型;
其中,渗流运动方程中的油相运动方程为
式中,K为渗透率张量;Kro为油相的相对渗透率,μo为油相的粘度;Φo为油相的流动势;G为启动压力梯度;为Hamilton算子;
3)对步骤2)中的水驱稠油区块的渗流数学模型进行有限差分离散化,并采用全隐式方法最终获得变启动压力梯度数值模拟模型;
S3、初始化压力和饱和度场;
S4、采用全隐式迭代方法求解水驱稠油区块的压力和饱和度;
S5、进行收敛条件与物质平衡检查;
若满足收敛性要求和物质平衡则输出当前时间步的启动压力梯度、压力以及饱和度等参数,进入下一个时间步的循环计算,若不满足上述条件,则需要调整迭代步长,重复步骤S4;
S6、检查是否达到模拟结束时刻,是则输出水驱稠油区块的启动压力梯度场、压力场以及饱和度场,否则重复步骤S1~S5。
2.如权利要求1所述的考虑变启动压力梯度的水驱稠油数值模拟方法,其特征在于:上述步骤S1中,油藏物性参数包括渗透率、孔隙度、流体的高压物性以及相渗曲线。
3.如权利要求1所述的考虑变启动压力梯度的水驱稠油数值模拟方法,其特征在于,在上述步骤S2的1)中,获得油相启动压力梯度和流度的关系模型的具体内容包括:
根据稠油区块的岩心驱替试验获得油相流度与启动压力梯度数据,对油相流度和启动压力梯度数据进行拟合处理,获得油相启动压力梯度随流度变化的关系曲线图及拟合方程;
拟合方程为:
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中海石油(中国)有限公司;中海石油(中国)有限公司北京研究中心,未经中海石油(中国)有限公司;中海石油(中国)有限公司北京研究中心许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910962175.3/1.html,转载请声明来源钻瓜专利网。