黑油模型数值模拟入门指南
作作作
作
者者者
者
:::
:
gulfmoon79
整整整
整
理理理
理
:::
:
Henry Li
二 零 零 八 年 八 月
精准石油论坛
Precision Petroleum Forum
本文字内容由精准石油论坛 gulfmoon79提供...
.
版权所有
The Following Content is Presented by gulfmoon79 @ Precision Petroleum Forum.
All Right Reserved
http://forum.petro-china.com
黑油模型数值模拟入门指南黑油模型数值模拟入门指南黑油模型数值模拟入门指南
黑油模型数值模拟入门指南
作者:ggg
g
ulfmoon79ulfmoon79ulfmoon79
ulfmoon79
第一:从掌握一套商业软件入手...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
111
1
第二:做油藏数值模拟都需要准备什么参数...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
222
2
第三:如何准备各部分参数...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
333
3
1.1.1.
1.
油藏模型...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
333
3
2.2.2.
2.
流体 PVTPVTPVT
PVT
参数...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
555
5
3.3.3.
3.
相渗和毛管压力曲线...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
777
7
4.4.4.
4.
分区设置...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
888
8
5.5.5.
5.
模型初始化...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
999
9
6.6.6.
6.
生产历史拟合...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
101010
10
7.7.7.
7.
如何最快完成历史拟合...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
111111
11
8.8.8.
8.
模型重启计算...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
131313
13
9.9.9.
9.
产量预测...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
131313
13
10.10.10.
10.
如何加快数模计算以及如何解决数模计算的收敛性问题?...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
151515
15
11.11.11.
11.
结束语...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
...
.
181818
18
HenryHenryHenry
Henry
LiLiLi
Li
整理于 2008-08-222008-08-222008-08-22
2008-08-22
精准石油论坛 数值模拟软件应用
- 1 -
黑油模型数值模拟入门指南黑油模型数值模拟入门指南黑油模型数值模拟入门指南
黑油模型数值模拟入门指南
Gulfmoon79
记得上大学最早学围棋时总感觉无从入手,看身边的朋友下棋时学着聂卫平从容入定,潇
洒自如的样子,很是羡慕。后来从书店买来围棋入门指南,夜深人静时照着 指南慢慢学如何吃
子,如何做眼,什么是打劫,怎么样布局。掌握了一点基本知识以后开始找水平最差的下,输了
一定不能弃擂,脸皮要厚,缠着对方接着下。赢了水平最差的人后去找中等水平的人下。这样经
过一年半载,再看以前那些学着聂卫平从容入定,潇洒自如下棋的同学,心想他们原来不过如
此,赶老聂差十万八千里哪。在这里也有许多人把我叫大师,专家,如果哪一天你觉得其实我的
水平也很一般,那你就到了专业段位了。
市场上有不少关于油藏数值模拟的书,但好像没有类似围棋入门指南那样从基础开始一步
一步介绍的书。我收到不下二十个问油藏数值模拟如何入门的问题。我尝试写一写油藏数值模拟
入门指南,希望对那些刚刚开始进入油藏数值模拟领域的工作者有所帮助。
第一:从掌握一套商业软件入手
我给所有预从事油藏数值模拟领域工作的人员第一个建议是先从学一套商业数值模拟软件
开始。起点越高越好,也就是说软件功能越强越庞大越好。现在在市场上流通的 ECLIPSE,VIP和
CMG都可以。如果先学小软件容易走弯路。有时候掌握一套小软件后再学商业软件会有心里障
碍。
对于软件的学习,当然如果能参加软件培训最好。如果没有机会参加培训,这时候你就需
要从软件安装时附带的
做起。油藏数值模拟软件通常分为主模型,数模前处理和数模后处
理。主模型是数模的模拟器,即计算部分。这部分是最重要的部分也是最难掌握的部分。它可以
细分为黑油模拟器,组分模拟气,热采模拟器,流线法模拟器等。数模前处理是一些为主模拟器
做数据准备的模块。比如准备油田的构造模型,属性模型,流体的 PVT参数,岩石的相渗曲线和
毛管压力参数,油田的生产数据等。数模后处理是显示模拟计算结果以及进行结果分析。
以 ECLIPSE 软件为例,ECLIPSE100,ECLIPSE300 和 FrontSim 是主模拟器。ECLISPE100
是对黑油模型进行计 算,ECLISPE300是对组分模型和热采模拟进行计算,FrontSim是流线法模
拟器。前处理模块有 Flogrid,PVTi,SCAL,Schedule,VFPi 等。Flogrid 用于为数值模拟建立模拟模
型,包括油田构造模型和属性模型;PVTi用于为模拟准备流体的 PVT参数,对于黑油模型,主要
是流体的属性随地层压力的变化关系表,对于组分模型是状态方程;SCAL为模型准备岩石的相渗
曲线和毛管压力输入参数;Schedule处理油田的生产数据,输出 ECLIPSE需要的数据格式(关键
字);VFPi是生成井的垂直管流曲线表,用于模拟井筒管流。ECLIPSE OFFICE和 FLOVIZ是后
处理模块,进行计算曲线和三维场数据显示和分析,ECLIPSE OFFICE同时也是 ECLIPSE的集成
平台。
精准石油论坛 数值模拟软件应用
- 2 -
对于初学者,不但要学主模型,也需要学前后处理。对于 ECLISPE的初学者,应该先从
ECLISPE OFFICE 学起,把 ECLISPE OFFICE 的安装练习做完。然后再去学 Flogrid,Schedule和
SCAL。PVTi主要用于组分模型,做黑油模型可以不用。
第二:做油藏数值模拟都需要准备什么参数
在照着软件提供的安装例子做练习时经常遇到的问题是:虽然一步一步按照手册的说明
做,但做的时候不明白每一步在做什么,为什么要这么做。这时候的重点在于你要知道你一开始
做的工作都是为数值模拟计算提供满足软件格式要求的基础参数。有了这些基础参数你才能开始
进行模拟计算。这些基础参数包括以下几个部分:
1. 模拟工作的基本信息:设定是进行黑油模拟,还是热采或组分模拟;模拟采用的单位制
(米制或英制);模拟模型大小(你的模型在 X,Y,Z三方向的网格数); 模拟模型网格
类型(角点网格,矩形网格,径向网格或非结构性网格);模拟油藏的流体信息(是油,
气,水三相还是油水或气水两相,还可以是油或气或水单相,有没有溶解气和挥发油等);
模拟油田投入开发的时间;模拟有没有应用到一些特殊功能(局部网格加密,三次采油,
端点标定,多段井等);模拟计算的解法(全隐式,隐压显饱或自适应)。
2. 油藏模型:模型在 X,Y,Z三方向的网格尺寸大小,每个网格的顶面深度,厚度,孔隙度,
渗透率,净厚度(或净毛比)。网格是死网格还是活网格。断层走向和断层传导率。
3. 流体 PVTPVTPVT
PVT
属性:油,气,水的地面密度或重度;油,气的地层体积系数,粘度随压力变
化表;溶解油气比随压力的变化表;水的粘度,体积系数,压缩系数;岩石压缩系数。如
果是组分模型,需要提供状态方程。
4. 岩石属性:相对渗透率曲线和毛管压力曲线。如果是油,气,水三相,需要提供油水,油
气相对渗透率曲线和毛管压力曲线(软件会自动计算三相流动时的相对渗透率曲线);如
果是油,水两相或气,水两相,只需要提供油水或气水两相相对渗透率曲线和毛管压力曲
线。
5. 油藏分区参数:如果所模拟的油田横向或纵向流体属性,岩性变化比较大,或者存在不同
的油水界面,这时需要对模型进行 PVT分区(不同区域用不同的 PVT流体参数表),岩
石分区(不同区域用不同的相对渗透率曲线和毛管压力曲线)或者平衡分区(不同平衡区
用不同的油水界面)。另外如果想掌握油藏不同断块的储量或采收率,可以对模型进行储
量分区(不同储量区可以输出不同的储量,产量,采收率,剩余储量等)。
6. 初始化计算参数:油藏模型初始化即计算油藏模型初始饱和度,压力和油气比的分布,从
而得到油藏模型的初始储量。这部分需要输入模型参考深度,参考深度处对应的初始压
力,油水界面以及气水界面;油气比或饱和压力随深度的变化;如果是组分模型,需要输
入组分随深度的变化。
7. 输出控制参数:即要求软件在计算时输出哪些结果参数。比如要求输出模型计算油田的
油,气,水产量变化曲线;油田压力变化曲线;单井油,气,水产量变化曲线;单井井底
压力变化曲线;单井含水,油气比变化曲线等。
精准石油论坛 数值模拟软件应用
- 3 -
8. 生产参数:对于已开发油田,这部分的数据量非常大。包括油田每口井的井位,井轨迹,
井的射孔位置,井的生产或注入历史(油,气,水产量,注入量,井底压力,井口压力
等),井的作业历史等。
第三:如何准备各部分参数
上面介绍了做油藏数值模拟所需要的参数,那么这些参数是如何得来的?又应该如何输入
到数模模型中哪?下面具体介绍各部分数据的来源即处理办法。
1.1.1.
1.
油藏模型
大多数油公司现在都在采用一体化工作流程,数模工程师不再需要自己去建立油藏模型。
地质学家对油田进行详细的油藏描述工作,基于地震解释数据,测井解释数据,岩芯数据以及结
合地质家对油田的认识建立三维地质模型。根据数模工程师的要求,地质家对三维地质模型进行
粗化处理,可以直接为数模工程师输出符合数模工程师需要的油藏模拟模型。
如果你不幸工作在一个没有很好地采用一体化工作流程的油公司或学校,那你就需要自己
建立油藏模拟模型了。这时你需要用到数模软件的前处理建模模块(比如Flogrid).如果你从地质家
处拿到的是三维地质模型,你的工作量还不是很大,主要是采用前处理建模模块对地质模型进行
粗化,生成数值模拟模型。如果你拿到的是两维数据体,那你的工作量就要大的多。首先你要确
定你必须要得到以下几方面的两维数据体:
每层的顶面深度
每层的厚度
每层的孔隙度分布
每层的渗透率分布
每层的净厚度或净毛比分布
断层数据
有了这些数据体,你就可以在前处理建模模块中建立三维数值模拟模型。
还有更不幸的情况,你有可能工作在软硬件环境都很落后的油田(比如国内的四川油
田),你根本连两维数据体也拿不到,你能拿到的只是一张一张的等值线图,或者只是单井井点
数据,这时候你首先需要做的是将这些等值线图件用数值化仪数值化为两位数据体,或者根据井
点数据生成等值线,然后再建三维模型。
你在开始建三维数模模型时首先需要确定你的模型应该采用什么样的网格类型。目前数值
模拟软件常用的网格有正交网格,角点网格,径向网格,非结构化网格(比如 PEBI)和动态网格等。你
需要知道这些网格的不同之处及其优缺点:
A:A:A:
A:
正交网格
精准石油论坛 数值模拟软件应用
- 4 -
正交网格是最常见网格,也是最早用来描述油藏的网格类型,目前仍然被广泛应用。由于其
计算速度快的特点,一些大型油气田经常采用此网格类型.有研究认为正交网格计算结果比其他网
格精确.
正交网格的数学描述也比较简单。以 ECLISPE为例,TOPS描述油藏顶部深度,DZ描述
油藏每层厚度,DX描述每个网格 X方向长度,DY描述每个网格 Y方向长度。
B:B:B:
B:
角点网格
角点网格的特点是网格的走向可以延着断层线,边界线或尖灭线,也就是说网格可以是扭
曲的。这样角点网格克服了正交网格的不灵活性,可以用来方便地模拟断层,边界,尖灭。但由于
角点网格网格之间不正交,这种不正交一方面给传导率计算带来难度,增加模拟计算时间,另一方面
也会对结果的精度有影响。
角点网格的数学模型很复杂,必须由前处理软件来生成。以 ECLISPE为例,COORD用来
描述模型网格的顶底坐标线(X,Y,Z),ZCORN描述每个网格八个角点的深度。
C:C:C:
C:
径向网格
径向网格比较简单,主要用于单井模拟。径向网格可以更合理的描述井附近流体的径向流
动。
D:D:D:
D:
非结构网格(PEBIPEBIPEBI
PEBI
网格)
PEBI网格源于 1908年就产生的 Voronoi网格。起主要特点是灵活而且正交.PEBI网格体系
提供了方便的方法来建立混合网格,比如模型整体采用正交网格,而对断层,井,边界等采用径
向,六边型或其他网格。网格间的传达率可以自动计算.PEBI网格的灵活性对模拟直井或水平井的
锥进问题非常有用。另外 PEBI网格可以用来精确模拟试井问题。还有 PEBI网格降低了网格走向
对结果的影响。PEBI网格的缺点是矩阵比其他网格要复杂的多,需要更加有效的解法。
E:E:E:
E:
动态网格
动态网格是指网格可以随时间而改变.通常用于动态网格加密或动态粗化。比如说在井生产
时采用局部加密而当井关闭时则采用正常网格。
有时候建立全油田整体模型后,对于压力及饱和度变化快的区域,常常需要进行局部网格加
密。采用局部加密可以准确的描述井附近流体的细微变化。网格局部网格可以是正交网格,或是
径向网格。Aziz认为(JPT 1993年)在正交网格中进行正交网格局部加密,有时并不会对结果有改
善。他建议采用混合网格,及在正交网格内采用径向网格加密,这样可以精确地模拟含水和气油比的
变化规律。
知道了这几种网格类型,那么在你建立模型时应该选择什么网格类型哪?你在建立网格时
又应该注意些什么哪?
A:A:A:
A:
在条件许可情况下尽量采用正交网格,而且尽量使网格保持均匀.尽量避免大网格直接连接小网
格,这样会带来严重的收敛问题。如果你的模型很大,最好采用正交网格。
B:B:B:
B:
角点网格已经非常成熟,但在建立角点网格时不要过分扭曲网格,要尽量保持模型的正交性。
如果你的边界与你的主断层相对平行,那么建立的网格系统正交性会比较好。你在建立网格后可
以用前处理软件计算模型正交性。
C:C:C:
C:
目前 PEBI网格在解法上还不成熟,应避免使用。但相信在五年内它会成为主导网格。
精准石油论坛 数值模拟软件应用
- 5 -
D:D:D:
D:
使用局部网格加密要小心,最小的网格不能小于井半径。而且局部网格加密部分要覆盖饱和度
变化大的网格。如果是水平井,局部网格加密要覆盖水平段。
E:E:E:
E:
网格越多模拟结果就越精确的概念是不对的。可以建立单井模型研究多大网格尺寸足够描述地
质上的非均质性。
F:F:F:
F:
网格走向会影响计算结果。在天然裂缝油气藏,要使网格走向与主裂缝方向一致。
G:G:G:
G:
DX/DY应接近于 1,不要大于 3。
H:H:H:
H:
井之间应有不少于三个网格。
I:I:I:
I:
模型的属性分布也很重要。尤其是数模人员自己插值,属性分布经常不合理。我见过许多模型,
井只是穿过单个的网格,井穿过网格的附近却全是死网格,这时模型计算根本不能收敛。数模人
员如果需要自己建模,最好要找地质人员参谋,或者自己有很好的地质背景。
初学者还有几个基本概念需要掌握:
AAA
A
:网格的 I,J,KI,J,KI,J,K
I,J,K
在做模拟时你需要知道每个网格的位置。模型每个网格的位置用网格的 I,J,K表示。在 ECLISPE软
件中,左上角的第一个网格为(1,1,1),在 X方向的第二个网格为(2,1,1),在 Y方向的
第二个网格为(1,2,1),在 Z方向的第二个网格为(1,1,2)。其他类推。
BBB
B
:死网格
死网格是模型中不参予流动计算的网格。通常把模型中的泥岩设为死网格。模型自动会把孔隙度
和净厚度为 0的网格设为死网格。如果你的油藏水体很大,你也不需要把水体都建在模型中。你
可以把水体部分设为死网格,然后用解析水体模拟油藏水体的贡献。
C:C:C:
C:
在数值模拟模型中如何正确应用有效孔隙度,净毛比?要明确区分定义.
总孔隙度:孔隙体积占总体积的百分数。
有效孔隙度:连通孔隙体积占总体积的百分数。
毛有效孔隙度:平均有效孔隙度。(泥质含量截至值)
净有效孔隙度:产层平均有效孔隙度(渗透率截至值)
净毛比:净厚度(渗透率截至值)与总厚度之比。
数值模拟模型中要用净有效孔隙度和净毛比,这样模型的体积计算是正确的。
2.2.2.
2.
流体 PVTPVTPVT
PVT
参数
我们前面讲过,数值模拟模型包括黑油模型,组分模型和热采模型。对不同的模型类型需
要用不同的模拟器。这些模型类型的区别主要就是在流体的 PVTPVTPVT
PVT
描述方面采用的方法不同。
对黑油模型,流体的 PVT属性描述方法是你直接给模型提供油,气,水 PVTPVTPVT
PVT
属性表(油
气体积系数,粘度,压缩系数随压力的变化;溶解油气比随压力的变化;水在参考压力下的体积
系数,压缩系数,粘度;岩石在参考压力下的压缩系数)。模型在计算过程中对每一个网格,根
据当前时间步的网格压力来查你提供的表得到每个网格当前的PVT属性。
组分模型要复杂的多。你不是提供 PVT表,而是提供状态方程(EOS)参数(每个组分的
临界温度,临界压力,临界 Z因子,分子量,偏心因子等),模型通过闪蒸计算来得到每个网格
精准石油论坛 数值模拟软件应用
- 6 -
当前的油气 PVT属性(粘度,体积系数,压缩系数,溶解油气比)。另外水和岩石的属性还需要
单独提供。
黑油模型和组分模型都假定油藏的温度在开采过程中是常数,不随时间变化。如果油藏温
度发生了变化,就需要用热采模型。对热采模型,最重要的流体属性变化是流体粘度随温度的变
化,即大家常说的粘温曲线。
那么你的模型应该采用哪种模型类型?这一方面取决于你拿到的 PVT实验
,如果你的
油藏是凝析气藏,挥发油藏,即你的油藏流体组分在生产过程中会发生明显变化,那你应该用组
分模型。另外如果你的油藏将来会考虑注气混相驱,那你也应该用组分模型。
重油热采用热采模型,其他大多数油气藏都可以用黑油模型。本部分只介绍黑油模型的
PVT属性,对组分和热采模型,希望以后有机会再介绍。
对于黑油模型,在 PVTPVTPVT
PVT
部分需要提供的参数包括:
油,气,水的地面密度或重度。
油的体积系数,粘度,溶解油气比随压力的变化。
气的体积系数,粘度随压力的变化。
参考压力下水的粘度,压缩系数,体积系数。
参考压力下岩石的压缩系数。
如果你的油藏在不同部位流体属性不同,那你需要提供多个 PVT表,每个表对应于模型的
某一部分。在后面的分区部分我还要做详细介绍。
对于黑油模型,根据不同的流体属性,可以分为以下四种类型:
死油油藏:油藏在整个开发阶段压力一直保持在泡点压力以上,在油藏中不会发生脱气。
生产油气比是常数(脱气发生在井筒或地面)。
活油油藏:油藏在开发阶段压力会降到泡点压力以下,在油藏中发生脱气,溶解油气比降
低,生产油气比增加(自由气会生产出来)。
干气藏:气藏压力在生产过程中不会穿过露点线,气藏中没有挥发油产生。
湿气藏:气藏压力在生产过程中会穿过露点线,气藏中产生挥发油。
对不同的黑油流体,PVT描述方法也不同,比如对