下载文档原格式
1 黑油模型理论基础1.1 基本假设（1）油藏中渗流是等温渗流；（2）油藏中有油、气、水三相，各相流体的渗流均符合达西定律； （3）模型考虑油组分、气组分、水组分三组分； （4）气组分在油气相、水气相之间发生质量交换； （5）相平衡瞬间完成；（6）水组分只存在于水相中，与油气相之间没有质量交换； （7）油藏岩石微可压缩，各向异性；（8）油藏流体可压缩，且考虑渗流过程中重力、毛管力的影响。
1.2 数学模型()()()()()rw w w w vw w w w ro o o o vo o o o rg so ro g g o o g g o o g sw rw so o sw w w w vg w w go w kk s p gD q B t B kk s p gD q B t B kk R kk p gD p gD B B s R kk R s R s p gD q B t B B B φρμφρμρρμμρφμ⎡⎤⎛⎫∂∇∇-+= ⎪⎢⎥∂⎣⎦⎝⎭⎡⎤⎛⎫∂∇∇-+= ⎪⎢⎥∂⎣⎦⎝⎭⎡⎤⎡⎤∇∇-+∇∇-+⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎛⎡⎤∂∇∇-+=++ ⎢⎥∂⎣⎦⎝⎧⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎛⎫⎫⎪ ⎪⎪ ⎪ ⎪⎪⎭⎝⎭⎩(1)辅助方程：1o w g cow o w cgo g o s s s p p p p p p ++=⎫⎪=-⎬⎪=-⎭(2)初始条件：()()()()()()000000,,,,,,,,,,,,,,,t w t w o t o p x y z t p x y z s x y z t s x y z s x y z t s x y z ===⎧=⎪⎪=⎨⎪=⎪⎩(3)边界条件：()()()()0,,,()(,,),,,,,L v v wf wf pn Q x y z t Q t x y z p x y z t p t x y z δδ∂⎧=⎪∂⎪=⎨⎪=⎪⎩(4)2 黑油模型程序整体结构图3 组员及分工4 主程序4.1 主程序主要功能（1）定义运算所需数组；（2）需要调用和生成的文件的打开和关闭；（3）通过调用子程序给模型赋基础数据和初始数据；（4）通过调用子程序给模型的运行确定各种控制；（5）在运行过程中反复读入井点数据（包括产量和注入量、井底流压、流动指数等）以及打印输出控制码；（6）分层计算油气水地质储量；（7）进行井点产量项处理和形成压力矩阵；（8）通过调用子程序求解压力方程；（9）显式计算饱和度；（10）进行过泡点处理；（11）根据压力和饱和度增量控制，自动调节计算时间步长；（12）未饱和网格块饱和度计算；（13）变泡点处理；（14）在每一运算时间步末进行物质平衡检验；（15）打印油井、水井、气井的分层报告和分井报告；（16）通过调用子程序进行计算结果打印输出和形成文件；（17）重启动运行方式的选择和重启动文件的生成和调用；（18）运算终止的控制及错误信息的提示。
CMG 油气藏数值模拟软件简介该软件能完成油藏开发方案设计即开发方案概念设计、详细开发方案设计和开发方案调整等工作。
主要包括：储量计算、油气井产能评价、开发层系划分、井网形式、井距、生产史拟合及地质模型修正、开发指标预测、开发方案指标的优选。
软件包括以下功能模块：（1）地质建模：主要通过测井解释成果、地震解释成果、储层地质及实验分析构造气藏三维地质模型，建立油藏网格属性。
（2）组分模型（用于凝析气藏和稠油热采）（3）黑油模型（用于油藏和一般气藏模拟计算）（4）相态计算模块（5）后处理模块。
是一个考虑重力及毛细管力的三相黑油模拟软件,网络系统可采用直角坐标,径向坐标,变深度/变厚度坐标,在任何网络系统中.都可建立两维或三维模型.在处理气相的出现及消失情况时,程序采用了变量替换方法。
一些主要特征和功能为：（1）、自适应隐式方法& nbspIMEX可以在显示,全隐式以及自适应隐式三种方式下运行。
在大多数情况下,只有很少一部分网格需要采用全隐式求解,而大部分网格都可采用显式方法求解.自适应隐式方法正是适合于这种情况的解法,并且在井附近以及层状油藏的薄层中,开采时会产生高速流动的锥进问题,采用自适应隐式处理这类问题是很有效的。
采用自适应隐式选项可节省三分之一到一半的运行时间。
计算时可采用和全隐式方法同样大的时间步长.用户可以指定采用全隐式方法计算的网格，可根据用户确定的界限或矩阵转换临界值，动态地选择采用全隐式计算的网络网格。
（2）、双孔/双渗双孔隙度选项允许采用两种方法对基岩模型进行离散化处理，其中一种为嵌套格式，成为“多重内部作用连续域”(MINC)方法，另一种为层状格式，称作“子区域”方法。
双孔隙模型对裂缝油藏进行了理想化的近似处理，认为裂隙油藏由两部分组成：主要孔隙度和次要孔隙度，主要孔隙度(基岩)代表岩块中的微小粒间孔隙，次要孔隙度(裂缝)由裂缝，通道和溶洞组成。
双孔隙模型将油藏分为两个连续域，裂隙是流体流动的主要通道，只具有很小的储集性能；而基岩具有较低的流体传导能力，但具有较大的储存能力。
名词解释:1油藏模拟：是用油藏模型来研究油藏的各种物理性质和流体在其中的流动规律，以便更好地认识油层，作出正确的评价，确定合理的开发方案和提高采收率的措施。
2 数值模型：用离散化方法将偏微分方程组转化为有限查分方程组，将其非线性系数线性化，得到线性方程组，然后求解。
3 油藏数值模拟：用数值方法求解油藏数学方程组，就是油藏数值模拟。
4 动态预测:在历史拟合的基础上对未来的开发指标进行计算。
5 黑油模型：黑油模型是简化的组份模型。
烃类系统只考虑两个组份:“油”组份是地层油经微分蒸发后在大气压的残存液（即黑油），而“气”组份是剩余的流体。
水相与其他两相不发生质量转移；气可以从油中出入，但油不能汽化为气相。
6 适定问题：一个问题的解存在，唯一且稳定时就称问题为适定问题。
7 三相流模型：描述有三相流体同时流动的数学模型。
8 三维模型：描述油藏流体沿三个方向上同时发生流动的数学模型。
9 气藏模型：描述天然气气藏的数学模型，有的气藏只有天然气的存在，而有的气藏不仅有天然气存在还有水存在。
10 离散化：离散化就是把整体分割为若干单元来处理。
11 有限差分法：有限差分法是对网格范围内的各点求解。
即原先表示连续的、足够光滑函数的偏微分方程，被一套对每个离散点的、与该点近似解有关的代数方程组所取代。
12 块中心网格：用网格分割成小块的中心来表示小块坐标。
13 一阶向前查商:对于函数p(x,t) ,x p p x p ii ∆-=∂∂+1 为一阶向前查商。
14 截断误差：当微商用查商表示时，把泰勒级数的余项截断，由于截断了泰勒级数的余项所产生的误差称为截断误差。
15 网格节点：网格的交点称为节点。
16 显式处理：在n+1时刻求解方程组时，若其系数直接用n 时刻的值，为显式处理。
17不均匀网格：为了模拟油藏的实际情况，划分网格时，在靠近井的附近网格取密一些，而沿径相外逐渐稀疏，这种网格称为不均匀网格。
18 IMPES 方法：是指隐式求解压力方程，显式求解饱和度方法。
1第六章黑油模型及其应用2第一节黑油模型及求解思路一、假设条件1、考虑油、气、水三相2、考虑油组分、气组分、水组分三个组分3、气组分在油、气相中要发生质量交换压力增加时，气组分可溶解在油相中（溶解气）压力降低时，气组分可从油相中分离出来（自由气）4、水相与气、油两相间无质量交换5、考虑毛管力、重力；油、气、水、岩石均可压缩6、油藏温度不变3二、数学模型1. 组分质量守恒方程()()⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛+∂∂=⎥⎦⎤⎢⎣⎡∇∇⋅∇⎥⎦⎤⎢⎣⎡∇−∇⋅∇−B S R B S t q D g p B kk R D g p B kk o o ggo so so gv o oo rog g g g rg φρμρμ＋＋()⎟⎟⎠⎞⎜⎜⎝⎛∂∂=+⎥⎦⎤⎢⎣⎡∇−∇⋅∇ooov o o o o ro B S t D g p B kk q φρμ式中R so —气油比（1）（2）（3）油组分水组分气组分()⎟⎟⎠⎞⎜⎜⎝⎛∂∂=+⎥⎦⎤⎢⎣⎡∇−∇⋅∇wwwv w w w w rw B S t D g p B kk q φρμ4辅助方程：1=++s s s w o g pp pp p p ogcgow o cow −=−=(4)(5)(6)未知量：g w o g w o s s s p p p 、、、、、2.未知量和辅助方程分析53. 初始条件和边界条件假设边界不规则的油藏中有若干口井生产或注入，求油藏中的压力和饱和度分布。
I.C 0xxL yL y()()(),,0,,00,,0O Oi w wi o oi P x y P S x y S t S x y S =⎫⎪=>⎬⎪=⎭B.C 1) 外边界封闭2)内边界>=∂∂Γt xP •定产•定流压P iwf P wft>0()ηζδ−−⋅=y x Q Q v v ,0>t61.隐式求压力1利用毛管压力消去则未知量减少为2利用饱和度归一化方程将油、气、水方程进行适当的组合和化简，最后得到一个只含有油相压力Po 的方程，称为压力方程。
黑油模型
黑油模型（Black Oil Model）是描述含有非挥发组分的黑油和挥发性组分
的原油溶解气两个系统在油藏中运动规律的数学模型。
也称低挥发油双组分模型。
黑油只是为了阐述油、气、水三相流体在油藏中的渗流规律而假设的一系列情况下的油藏模型。
这种模型不是指“石油”本身在物理性质和化学性质上有什么不同，而是指“油”在渗流规律和与油藏其它性质上的不同。
它只是通过物理、数学的手段人造的一个模型，便于进行油藏的模拟开发。
模型中烃类系统可用两组分描述：⑴非挥发组分(黑油)；⑵挥发组分，即溶于油中的气。
黑油模型也称低挥发油双组分模型。
IMPES方法原理方法：处理过程：方程左边系数（毛管力和传导系数）以显示形式取值；求解过程：压力和饱和度求解交替进行（先求压力再求饱和度）。
求解过程（差分方程线性化过程）：对于油气水三相的方程组：两两（3）压力和饱和度联立求解。
求解步骤：（1）方程的系数以隐式的形式展开，对方程的求解过程中，进行若干次迭代：第L+1次迭代的系数是以L次为基础进行泰勒展开，取一阶小项，忽略掉二阶小项，获得线性方程组，联立起来求压力和饱和度（相当于一个半隐式的步骤）。
（2）重复这个迭代过程，满足精度的要求，进入下一时间步。
（3）重复上述步骤，直到模拟时间结束。
IMPIMS方法基本原理：（1）毛管压力和传导系数的取值确定：以显式确定压力，以隐式确定饱和度，显式和隐式的取值过程交替进行；（2）求解压力和饱和度均以隐式迭代求得；（3）压力和饱和度的求解交替进行。
求解步骤：一、求解压力，与IMPES求压力的过程一样：对于油气水三相的方程组：（1）方程左边系数（毛管压力和传导系数等）以显式的形式取值（全都写成n方），将P n+1在P n处做泰勒展开，取一阶小项，即P n+1=P n+δp。
（2）方程右边以两相或三相的形式对求解变量进行分步微分展开（求解变量：三相：δp、δSw、δSg ||两相：δp、δSw、δPb（饱和压力））。
（3）将左右两边的表达式合并在一起，把与求解变量有关的内容移到等号左边，把常数项移动到等号右边，形成一个与三个求解变量有关的线性方程组。
（4）求解这个线性方程组。
考虑到δSw和δX只和本节点有关，消去这两项，得到一个与压力有关的方程，对压力用隐式求解。
得到δp。
二、求解饱和度，执行一遍半隐式方法，得到油气水三相方程。
对于油气水三相的方程组：（1）对方程左边中的相毛管压力和传导系数以隐式形势取值，以两相或三相的形式（三相以δp、δSw、δSg ||两相以δp、δSw、δPb）对毛管压力和传导系数进行处理：以n时刻为基准对n+1时刻的变量进行泰勒级数展开，取一阶小项，忽略掉两个微小项的乘积项。
1 黑油模型理论基础1.1 基本假设（1）油藏中渗流是等温渗流；（2）油藏中有油、气、水三相，各相流体的渗流均符合达西定律； （3）模型考虑油组分、气组分、水组分三组分； （4）气组分在油气相、水气相之间发生质量交换； （5）相平衡瞬间完成；（6）水组分只存在于水相中，与油气相之间没有质量交换； （7）油藏岩石微可压缩，各向异性；（8）油藏流体可压缩，且考虑渗流过程中重力、毛管力的影响。
1.2 数学模型()()()()()rw w w w vw w w w ro o o o vo o o o rg so ro g g o o g g o o g sw rw so o sw w w w vg w w go w kk s p gD q B t B kk s p gD q B t B kk R kk p gD p gD B B s R kk R s R s p gD q B t B B B φρμφρμρρμμρφμ⎡⎤⎛⎫∂∇∇-+= ⎪⎢⎥∂⎣⎦⎝⎭⎡⎤⎛⎫∂∇∇-+= ⎪⎢⎥∂⎣⎦⎝⎭⎡⎤⎡⎤∇∇-+∇∇-+⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎛⎡⎤∂∇∇-+=++ ⎢⎥∂⎣⎦⎝⎧⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎛⎫⎫⎪ ⎪⎪ ⎪ ⎪⎪⎭⎝⎭⎩(1)辅助方程：1o w g cow o w cgo g o s s s p p p p p p ++=⎫⎪=-⎬⎪=-⎭(2)初始条件：()()()()()()000000,,,,,,,,,,,,,,,t w t w o t o p x y z t p x y z s x y z t s x y z s x y z t s x y z ===⎧=⎪⎪=⎨⎪=⎪⎩(3)边界条件：()()()()0,,,()(,,),,,,,L v v wf wf pn Q x y z t Q t x y z p x y z t p t x y z δδ∂⎧=⎪∂⎪=⎨⎪=⎪⎩(4)2 黑油模型程序整体结构图3 组员及分工4 主程序4.1 主程序主要功能（1）定义运算所需数组；（2）需要调用和生成的文件的打开和关闭；（3）通过调用子程序给模型赋基础数据和初始数据；（4）通过调用子程序给模型的运行确定各种控制；（5）在运行过程中反复读入井点数据（包括产量和注入量、井底流压、流动指数等）以及打印输出控制码；（6）分层计算油气水地质储量；（7）进行井点产量项处理和形成压力矩阵；（8）通过调用子程序求解压力方程；（9）显式计算饱和度；（10）进行过泡点处理；（11）根据压力和饱和度增量控制，自动调节计算时间步长；（12）未饱和网格块饱和度计算；（13）变泡点处理；（14）在每一运算时间步末进行物质平衡检验；（15）打印油井、水井、气井的分层报告和分井报告；（16）通过调用子程序进行计算结果打印输出和形成文件；（17）重启动运行方式的选择和重启动文件的生成和调用；（18）运算终止的控制及错误信息的提示。
油藏数值模拟基本过程一、数值模拟发展概况30年代人们开始研究地下流体渗流规律并将理论用于石油开发；50年代在模似计算的方法方面，取得较大进展；60年代起步，人们开始用计算机解决油田开发上的一些较为简单间题，由于当时计算机的速度只有每秒几万到几十万次，实际上只能做些简单的科学运算;70 年后主要体现于计算机的快速升级带动了油藏数模的迅猛发展，大型标量机计算速度达到100——500万次，内存也高增主约16兆字节.在理论上黑油模型计算方法更趋成熟，D. W。
Peaceman的<油藏数值模似基础〉以及K。
Aziz和A。
Settari的〈油藏模似〉等主要著作都是在这个阶段出版的，但仍受到计算机速度和内存的限制,使用的方法一般仅限于IMPES及半隐式等，只能解决中小型油藏的模拟应用问题；80年代则是油藏数值模似技术飞跃发展的年代，解决不同类型油藏的数模计算方法及软件相应问世,同时超级向量机的诞生，使计算机速度达到亿次，甚至几十亿次,内存高达10—20亿字节。
90年代特别是后期，油藏模似软件各模块功能也有了惊人的发展，主要体现为向一体化方面发展;即集地震、测井、油藏工程（数模）、工艺及地面集输、经济评价等为一体的大型软件方面发展。
目前油藏数值模似软件基本上形成了一套能处理各种类型油气藏和各种不同开采方式的软件系列。
?黑油模型已被广泛用于各种常规油气藏的模拟；?裂缝模型可用来解决除砂岩以外的灰岩、花岗岩、凝灰岩和变质岩的裂缝性油气藏开发问题；?组分模型用于凝析气藏、轻质油、挥发油藏的开发设计和混相驱的研究；？热采模型用于稠(重)油油藏蒸气吞吐、蒸汽驱和就地燃烧的设计；？化学驱模型用于在注入水中添加聚合物、表面活性剂、碱等各种化学剂进行三次采油提高采收率的计算和设计。
油藏数值模拟方法的新突破随着计算机运算速度的提高,向量算法的出现和应用是软件设计上一个划时代的发展。
预处理共轭梯度法更快速、有效地解各种更为复杂和困难的大型稀疏线性方程组.网格化方面不局限于静态和动态的局部网格加密技术，不规则网格、PEBI网格的出现更好的解决了在边界、断层插值计算以及面与面垂直正交的新型数模计算方法，更快速收敛。
油藏数值模拟中⼏种主要的数学模型1、⿊油模型（Black Oil ）：⿊油模型是指⾮挥发性原油的数学模型，是相对于油质极轻的挥发性油⽽⾔，因油质重⽽⾊泽较深，故称之为⿊油其基本假设为：<1> 油藏中的渗流为等温渗流； <2> 油藏中最多只有油⽓⽔三相，每⼀相的渗流均遵守达西定律；<3> 油藏烃类只含有油⽓两个组分，油组分是指将地层原油在地⾯标准状况下经历分离后所残存的液体，⽽其组分是指全部分离出来的天然⽓。
油藏状况下油⽓两种组分可能形成油⽓两相，油组分完全存在于油相中，⽽⽓组分则可以以⾃由⽓的形式存在于⽓相内，也可以以溶解⽓的⽅式存在于油相中，所以地层中油相应为油组分和⽓组分的某种组合。
常规⿊油模型⼀般不考虑油组分向⽓组分的挥发过程；<4> 油藏中⽓体的溶解和逸出是瞬间完成的，即认为油藏中油⽓两相瞬时地达到相平衡状态；<5> 油⽔之间不互溶；<6> 由于天然⽓在⽔中溶解度很⼩，可以认为它不溶于⽔。
油⽓⽔三相渗流基本微分⽅程：g ()()()()[()]()()ro o o o o o o o ro gd rg g gd o g g o og g g s o go g rw w w w w ww w kk S P D q t kk kk S S P D P D R q q t kk S P D q t ρφργµρρφρφργγµµρφργµ-?+=??+?-?+??-?++=?-?+=油相：⽓相：⽔相：油⽔两相渗流基本微分⽅程：g()()()()ro og og o o o o o rw w w w w w w w kk S P D q t kk S P D q t ρφργµρφργµ-?+=-?+=?油相：⽔相：注意：1、式中的产量项是以质量计的单位时间内单位地层体积的产出（注⼊）量；2、og o gd ρρρ=+，地⾯油的相对密度为地⾯油与溶解⽓相对密度之和。
ANS黑油模型油藏数值模拟软件说明书(仅供内部使用)刘月田中国石油大学(北京)复杂科学研究中心2001—2011前言本软件可用来模拟油、气、水三相同时渗流（三相），地层和流体性质及压力、饱和度在X、Y、Z三个方向上变化（三维）条件下的油田注水、注气开发过程。
本软件数据输入采用固定格式方法。
本软件输出结果除综合数据文件外，还可单独输出压力、饱和度分布等单项数据文件。
本软件具有多次连续模拟功能。
本软件共分三部分：输入文件，模拟软件，输出文件。
一、数据输入（文件）系统采用固定输入格式，数据文件名：ANS1.DAT； ANS2.DAT二、模拟运算系统：ANS.EXE三、数据输出（文件）系统1、综合结果文件：ANS.RES2、典型时刻单一生产项指标和压力、饱和度等的分布3、重新连续运行文件：ANSTRT.DAT三维三相黑油油藏数值模拟软件数据输入说明Ⅰ《网格系统》1、标题（文字说明行—在使用者编辑数据时起提示作用，对程序运行不提供任何信息）2、网格块数（3个整数）II=X方向网格块数JJ=Y方向网格块数KK=Z方向网格块数3、标题（文字说明行—在使用者编辑数据时起提示作用，对程序运行不提供任何信息）4、输入用到的代码(3个整数)KDX—控制X方向网格大小输入方式的代码KDY—控制Y方向网格大小输入方式的代码KDZ—控制Z方向网格大小输入方式的代码代码含义：KDX=-1 X方向网格大小相同（只需输入一个值）。
KDX=0 第一层（K=1）第一行的每一网格大小须读入。
X方向相同列的网格相同（须读入II个值）KDY=-1 y方向上网格大小相同（只须读入一个值）。
KDY=0 y方向第一层（K=1）第一行的每一网格大小都需读入。
Y方向相同列的网格相同（须读入JJ个值）。
KDZ=-1 Z方向网格大小（网格厚度）相同（只须读入一个值）。
KDZ=0 网格中各层读入不同厚度值（须读入KK个值）。
KDZ=+1 每个网格的Z方向步长都须一一读入（须读入II×JJ×KK个值）。
一.做油藏数值模拟都需要准备什么参数1。
模拟工作的基本信息：设定是进行黑油模拟，还是热采或组分模拟；模拟采用的单位制（米制或英制）；模拟模型大小（你的模型在X,Y,Z三方向的网格数）；模拟模型网格类型（角点网格，矩形网格，径向网格或非结构性网格）；模拟油藏的流体信息（是油，气，水三相还是油水或气水两相，还可以是油或气或水单相，有没有溶解气和挥发油等）;模拟油田投入开发的时间；模拟有没有应用到一些特殊功能（局部网格加密，三次采油，端点标定，多段井等）；模拟计算的解法（全隐式，隐压显饱或自适应）。
2。
油藏模型：模型在X,Y,Z三方向的网格尺寸大小，每个网格的顶面深度，厚度，孔隙度，渗透率，净厚度（或净毛比）。
网格是死网格还是活网格。
断层走向和断层传导率。
3。
流体PVT属性：油，气，水的地面密度或重度；油，气的地层体积系数，粘度随压力变化表；溶解油气比随压力的变化表；水的粘度，体积系数，压缩系数；岩石压缩系数。
如果是组分模型，需要提供状态方程。
4。
岩石属性：相对渗透率曲线和毛管压力曲线。
如果是油，气，水三相，需要提供油水，油气相对渗透率曲线和毛管压力曲线（软件会自动计算三相流动时的相对渗透率曲线）；如果是油，水两相或气，水两相，只需要提供油水或气水两相相对渗透率曲线和毛管压力曲线。
5。
油藏分区参数：如果所模拟的油田横向或纵向流体属性，岩性变化比较大，或者存在不同的油水界面，这时需要对模型进行PVT分区（不同区域用不同的PVT流体参数表），岩石分区（不同区域用不同的相对渗透率曲线和毛管压力曲线）或者平衡分区（不同平衡区用不同的油水界面）。
另外如果想掌握油藏不同断块的储量或采收率，可以对模型进行储量分区（不同储量区可以输出不同的储量，产量，采收率，剩余储量等）。
6。
初始化计算参数：油藏模型初始化即计算油藏模型初始饱和度，压力和油气比的分布，从而得到油藏模型的初始储量。
这部分需要输入模型参考深度，参考深度处对应的初始压力，油水界面以及气水界面；油气比或饱和压力随深度的变化；如果是组分模型，需要输入组分随深度的变化。
第一章油藏数值模拟方法分析1.1油藏数值模拟1.1.1油藏数值模拟简述油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。
其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。
其基础理论是基于达西渗流定律。
油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。
基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值(生产动态)与实际吻合。
其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。
充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。
这组流动方程组由运动方程、状态方程和连续方程所组成。
油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。
具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。
油藏数值模拟技术从50 年代的提出到90 年代间历经40 年的发展，日益成熟。
现在进入另外一个发展周期。
近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。
在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用非常广泛。
油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。
图1 油藏数值模拟流程图1.1.2油藏数值模拟的类型油藏数值模拟类型的划分方法有多种，划分时最常用的标准是油藏类型、需要模拟的油藏流体类型和目标油藏中发生的开采过程，也可以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定，油气藏特性和油气性质不同，选择的模型也不同，还可以根据油藏数值模拟模型所使用的坐标系、空间维数和相态数来划分。
Eclipse黑油模型数值模拟入门指南Eclipse黑油模型数值模拟入门指南记得上大学最早学围棋时总感觉无从入手，看身边的朋友下棋时学着聂卫平从容入定，潇洒自如的样子，很是羡慕。
后来从书店买来围棋入门指南，夜深人静时照着指南慢慢学如何吃子，如何做眼，什么是打劫，怎么样布局。
掌握了一点基本知识以后开始找水平最差的下，输了一定不能弃擂，脸皮要厚，缠着对方接着下。
赢了水平最差的人后去找中等水平的人下。
这样经过一年半载，再看以前那些学着聂卫平从容入定，潇洒自如下棋的同学，心想他们原来不过如此，赶老聂差十万八千里哪。
在这里也有许多人把我叫大师，专家，如果哪一天你觉得其实我的水平也很一般，那你就到了专业段位了。
市场上有不少关于油藏数值模拟的书，但好像没有类似围棋入门指南那样从基础开始一步一步介绍的书。
我收到不下二十个问油藏数值模拟如何入门的问题。
我尝试写一写油藏数值模拟入门指南，希望对那些刚刚开始进入油藏数值模拟领域的工作者有所帮助。
第一：从掌握一套商业软件入手。
我给所有预从事油藏数值模拟领域工作的人员第一个建议是先从学一套商业数值模拟软件开始。
起点越高越好，也就是说软件功能越强越庞大越好。
现在在市场上流通的ECLIPSE,VIP和CMG都可以。
如果先学小软件容易走弯路。
有时候掌握一套小软件后再学商业软件会有心里障碍。
对于软件的学习，当然如果能参加软件培训最好。
如果没有机会参加培训，这时候你就需要从软件安装时附带的练习做起。
油藏数值模拟软件通常分为主模型，数模前处理和数模后处理。
主模型是数模的模拟器，即计算部分。
这部分是最重要的部分也是最难掌握的部分。
它可以细分为黑油模拟器，组分模拟气，热采模拟器，流线法模拟器等。
数模前处理是一些为主模拟器做数据准备的模块。
比如准备油田的构造模型，属性模型，流体的PVT参数，岩石的相渗曲线和毛管压力参数，油田的生产数据等。
数模后处理是显示模拟计算结果以及进行结果分析。
以ECLIPSE软件为例，ECLIPSE100,ECLIPSE300和FrontSim是主模拟器。
油藏数值模拟求解中的IMPIMS方法
肖阳; 张茂林; 梅海燕
【期刊名称】《《石油工业计算机应用》》
【年(卷),期】2006(14)1
【摘要】在油藏数值模拟中对数值模型进行求解,目前基本有两类方法一类是顺序求解法,如IMPES方法和IMPIMS方法;一类是联立求解法,如半隐式法和全隐式法。
IMPIMS方法是一种顺序求解方法,是IMPES和半隐式方法的混合变种。
它既具有半隐式求解饱和度的特点,又保留了IMPES方法省内存,省计算工作量的优势。
【总页数】3页(P26-28)
【作者】肖阳; 张茂林; 梅海燕
【作者单位】西南石油大学石油工程学院四川成都 610500
【正文语种】中文
【中图分类】TE3
【相关文献】
1.油藏数值模拟方法在复杂小断块油藏开发调整中的应用 [J], 邓江洪
2.油藏水平井数值模拟研究中一种新的求解方法 [J], 孙洪林;陈钦雷;刘(忄尉)宁
3.油藏数值模拟求解中的IMPIMS方法 [J], 肖阳;张茂林;梅海燕
4.岩石分类方法在油藏数值模拟中的应用 [J], 姜彬;杨希濮;李珂;陈翰;李雪
5.油藏数值模拟求解中的IMPIMS方法 [J], 肖阳; 张茂林; 梅海燕
因版权原因，仅展示原文概要，查看原文内容请购买。
