[发明专利]中高渗油藏注水开发过程中物性时变数值模拟方法及装置在审
| 申请号: | 202111579518.1 | 申请日: | 2021-12-22 |
| 公开(公告)号: | CN116335599A | 公开(公告)日: | 2023-06-27 |
| 发明(设计)人: | 钟小刚;罗波;朱红云;王博;郭丽丽;孙海涛;李淑静;李春雷;宋辉;何丕祥;杨衍云 | 申请(专利权)人: | 中国石油天然气股份有限公司 |
| 主分类号: | E21B43/20 | 分类号: | E21B43/20;E21B49/00;G06F30/20 |
| 代理公司: | 天津才智专利商标代理有限公司 12108 | 代理人: | 庞学欣 |
| 地址: | 100007 北京市*** | 国省代码: | 北京;11 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | |||
| 搜索关键词: | 中高 油藏 注水 开发 过程 物性 数值 模拟 方法 装置 | ||
1.一种中高渗油藏注水开发过程中物性时变数值模拟方法,其特征在于,具体包括如下步骤:
将油藏或储层进行划分以获取网格单元;
获取网格单元的地层油粘度、不同注水冲刷倍数下油水相对渗透率曲线,同时间步下相邻网格单元之间的传导率建立油藏模拟模型;
基于油藏模拟模型,获取网格单元的注水倍数,并基于所述网格单元的注水倍数获取所述网格单元渗透率、地层油粘度和油水相对渗透率曲线。
2.根据权利要求1所述的方法,其特征在于,采用如下方法获取同时间步下相邻网格单元之间的传导率:
获取网格单元的初始渗透率、注水冲刷倍数、注水冲刷后网格单元的最大渗透率以获取网格单元的渗透率,并获取不同时间步下的单元格渗透率;
基于不同时间步下的单元格渗透率获取不同时间步下相邻网格单元之间的传导率。
3.根据权利要求2所述的方法,其特征在于,采用如下公式构建渗透率时变模型以获取不同时间步下的单元格渗透率:
k(x)=kmax-(kmax-kint)·exp(-bkx)
kint—网格单元的初始渗透率;
kmax—注水冲刷后网格单元的最大渗透率;
k(x)—网格单元渗透率(mD),
x—注水冲刷倍数(注水冲刷孔隙体积倍数),无因次;
bk—渗透率时变因子,无因次;
采用如下公式获取不同时间步下相邻网格单元之间的传导率:
Ai—网格单元i与相邻的网格单元j之间的渗流面积;
Aj—网格单元j与相邻的网格单元i之间的渗流面积;
di—网格单元i到相邻网格单元j间的渗流距离;
dj—网格单元j与相邻的网格单元i之间的渗流距离;
Ti—网格单元i到相邻的网格单元j的传导率;
且
Tj—网格单元j到相邻的网格单元i的传导率,
且
Tij—相邻网格单元i和j之间的时变传导率;
u—数值模拟计算的迭代步,n+1为时间步;
ki(xi)—网格单元i在注水冲刷倍数xi时的渗透率;
kj(xj)为网格单元j在注水冲刷孔隙体积倍数为xj条件下的渗透率。
4.根据权利要求1所述的方法,其特征在于,所述获取网格单元的地层油粘度的步骤包括:
获取网格单元油藏开发初期地层原油粘度、注水冲刷后地层油的最大粘度,注水冲刷倍数以获取网格单元的地层油粘度。
5.根据权利要求4所述的方法,其特征在于,构建如下公式以获取网格单元的地层油粘度:
μ(x)=μmax-(μmax-μint)·exp(-bμx);
μint—网格单元油藏开发初期地层原油粘度;
μmax—注水冲刷后得到的最大地层油粘度;
μ(x)—网格单元地层油粘度;
x—注水冲刷倍数,无因次;
bμ—地层油粘度时变因子,无因次。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气股份有限公司,未经中国石油天然气股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202111579518.1/1.html,转载请声明来源钻瓜专利网。





