[发明专利]一种基于拟流线法的底水油藏水平井水驱数值模拟方法有效
申请号: | 201910134624.5 | 申请日: | 2019-02-23 |
公开(公告)号: | CN109858177B | 公开(公告)日: | 2022-06-07 |
发明(设计)人: | 苏彦春;侯亚伟;李彦来;黄凯;罗宪波;蔡晖;张东;赵秀娟;牟松茹;谭捷 | 申请(专利权)人: | 中国海洋石油集团有限公司;中海石油(中国)有限公司天津分公司 |
主分类号: | G06F30/20 | 分类号: | G06F30/20 |
代理公司: | 天津市北洋有限责任专利代理事务所 12201 | 代理人: | 曹玉平 |
地址: | 100010 北京市东城区朝*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 流线 底水 油藏 水平 井水 数值 模拟 方法 | ||
1.一种基于拟流线法的底水油藏水平井水驱数值模拟方法,其特征在于,包括如下步骤:
步骤1、收集油藏水平井基本静态和动态信息获取数值模拟计算参数;
步骤2、对数值模拟计算参数中的原始相渗数值进行归一化计算获得模拟特征参数;
步骤3、根据模拟特征参数按照拟流线方法对拟流区域进行等体积网格划分建立拟流区域网格模型;根据水脊剖面描述公式,沿油水界面等距离确定拟流线,在两流线之间的区域即为拟流管,沿流管方向进行等体积剖分;
步骤4、对拟流区域模型中每个单独拟流区域进行驱替计算;
步骤5、计算整个拟流区域模型二维流场含水率叠加生成井口含水率;
步骤6、判断井口含水率是否大于或等于截止条件;如果满足,则进入下一步;否则返回步骤4;
步骤7、输出整个拟流区域模型数值;
所述步骤4对拟流区域模型中每个拟流区域进行驱替计算过程:
4.1、根据拟流区域模型初始化,按照如下公式计算拟流区域模型初始状态下的渗流阻力,计算公式:
式中,Rl为流管l的渗流阻力,MPa﹒d﹒m-3;Ln为流管l的第n个网格的长度,m;An为流管l的第n个网格的截面积,m2;kron为流管l的第n个网格的油相相对渗透率,mD;krwn为流管l的第n个网格的水相相对渗透率,mD;μo为原油地层粘度,mPa﹒s;μw为水的地层粘度,mPa﹒s;
4.2、设定生产压差为控制条件,计算第一个时间步的流量,计算公式:
式中,为流管l在t1时刻的流量,m3/d;pe为流管入口端压力,MPa;pf为流管出口端压力,MPa;
4.3、根据水驱前缘追踪理论,代入当前时间步下的流量,计算水驱前缘沿流管推进的距离,计算公式:
式中,ξ为表示沿流管的网格标号;fw表示出口端含水率,小数;Sw表示出口端含水饱和度,小数;
4.4、根据等饱和度面移动方程,计算拟流区域起点与水驱前缘之间的饱和度分布,建立拟流区域的饱和度剖面;
4.5、根据拟流区域出口处网格的饱和度,根据相渗曲线分流量方程计算出口端含水率;
4.6、按照上面的计算步骤对每个拟流区域进行计算,得到各个拟流区域出口处的流量、含水率及流管的饱和度分布,按照如下公式并对拟流区域出口处指标进行求和,即可得到整管个二维流场的含水率;
式中,f为二维流场的总含水率,小数;fwl为流管l的出口端含水率,小数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国海洋石油集团有限公司;中海石油(中国)有限公司天津分公司,未经中国海洋石油集团有限公司;中海石油(中国)有限公司天津分公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910134624.5/1.html,转载请声明来源钻瓜专利网。