一、关于“油藏数值模拟技术”
（一）基本概念及作用
（二）数据准备
（三）模型初始化
（四）生产史拟合
（五）生产动态预测
二、油藏数值模拟的主流软件系统简介
三、油藏数值模拟技术的进展及发展方向
（一）进展
（二）发展方向
一、关于“油藏数值模拟技术”
油藏数值模拟技术是一门将油田开发重大决策纳入严格科学轨道的关键技术。从油田投产开始，无论是单井动态，还是整个油田动态，都要进行监测与控制。油藏数值模拟是油田开发最优决策的有效工具。
油藏数值模拟技术从20世纪50年代开始研究至今，已发展成为一项较为成熟的技术，在油气藏特征研究、油气田开发方案的编制和确定、油气田开采中生产措施的调整和优化以及提高油气藏采收率方面，已逐渐成为一种不可欠缺的主要研究手段。油藏数值模拟技术经过几十年的研究有了大的改进，越来越接近油气田开发和生产的实际情况。油藏数值模拟技术随着在油气田开发和生产中的不断应用，并根据油藏工程研究和油藏工程师的需求，不断向高层次和多学科结合发展，将得到不断的发展和完善。
（一）基本概念及作用
（1）基本概念
油藏数值模拟：从地下流体渗流过程中的本质特征出发，建立描述渗流过程基本物理现象、并能描述油藏边界条件和原始状况的数学模型，借助计算机计算求解渗流数学模型，结合油藏地质学、油藏工程学重现油田开发的实际过程，用来解决实际问题。
油藏数学模型的分类，一般有四种方法：
1)按流体中相的数目，划分为：单相流模型、两相流模型、三相流模型。
2)按空间维数，划分为：零维模型、一维模型、二维模型、三维模型。
3)按油藏特性类型，划分为：气藏模型、黑油模型、组分模型。气藏模型按其组分的贫富，可以用黑油数值模型模拟，也可以用组分类型的数值模拟模型模拟。所以，气藏模型也可以划进黑油或组分模型。故数学模型一般分为黑油型和组分型两类模型。
4)按油藏结构特点、开采过程特征，分类为：裂缝模型、热采模型、化学驱
模型、混相驱模型、聚合物驱模型等。
其中：
数学模型：通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。
黑油模型：描述油、气、水三相同时存在的油藏数学模型，一般认为，只有天然气可以溶于油中或从油中分离出来，油和水及气和水之间不发生质量交换。
组分模型：描述油藏内碳氢化合物化学组分的数学模型。
（2）油藏数值模拟技术的主要内容
油藏数值模拟技术的主要内容包括以下四部分：
第一部分：建立数学模型，也就是要建立一套描述油藏渗流的偏微分方程组解此方程组，还要有相应的辅助方程、初始条件和边界条件。
第二部分：建立数值模型。需要三个过程：首先，通过离散化，将偏微分方程组转换成有限差分方程组；然后，将其非线性系数项线性化，从而得到线性代数方程组；再通过线性方程组解法，求得所需求的未知量(压力、饱和度、温度、组分等)的分布及变化。
第三部分：建立计算机模型，也就是将各种数学模型的计算方法编制成计算机程序，以便用计算机进行计算，得到所需要的各种结果。
第四部分：油藏数值模拟应用研究。有了计算机模型后，下一步就是对实际研究对象，根据不同的研究目的进行模拟计算。
简单地说，油藏模拟就是把油藏在三维的空间里分为许多离散的单元，并且模拟油藏及流体在空间及一系列离散的时间步里的发展变化。与物质平衡方法一样的是，系统遵循物质守横原理。
建立油藏模拟软件，一般需要基于流体渗流力学原理等，如：
1）质量守恒原理
2）能量守恒原理
3）运动方程（达西定律）
4）状态方程
5）辅助方程（如饱和度方程，毛管力方程等）
（3）油藏数值模拟的作用
油藏模型分黑油模型和组分模型。它在石油工业中有着广泛的应用领域，可解决油田开发中的许多问题。由于模型不同，则侧重面不相同。现以黑油模型为例，阐述油藏数值模拟的用途。
A、模拟初期开发方案。
1)实施方案的可行性评价。
2)选择井网、开发层系、井数和井位。
3)选择注水方式。
4)对比不同的产量效果。
5)对油藏和流体性质的敏感性进行研究。
B、对已开发油田的历史进行模拟。
1)确定产液量和生产周期。
2)确定油藏和流体特性，拟合全油田和单井的压力、含水率(气油比)等历史动态。
3)指出问题及潜力所在区域。
C、动态预测。
动态预测包括枯竭开采、注水开发、钻加密井等，首先要确定限定条件，比
如：生产井最大采出量和注水井最小注水量、最小井口压力、井底流压或生产井最大压降、注入井最大井底压力。根据预测情况来确定最佳方案，以达到最大采收率。
1)评价提高采收率的方法。
①一次采油；
②注水；
③注气。
④注聚合物；
⑤注表面活性剂；
⑥注C02和其他混相驱；
⑦注蒸汽；
⑧火烧；
⑨几种方法的混合。
2)研究较高的剩余油饱和度分布。
①研究剩余油饱和度的分布范围和类型；
②单井进行调正，改变液流方向，改变注采井别
③扩大水驱油波及系数；
改变注水层位；
④回答油田开发中所遇到的问题并提出解决问题的方法
3)评价潜力和提高采收率的方向。
①确定井位、加密井的位置；
②确定产量、开采方式；，
③确定地面和井的设备；
④各种调整开发方案指标对比及经济评价。
D、专题和机理问题的研究。．
1)对比注水、注气和天然枯竭开采动态；
2)研究各种注水方式的效果；
3)研究井距、井网对油藏动态的影响；
4)研究不同开发层系对油藏动态的影响；
5)研究不同开发方案的各种指标；，
6)研究单井产量对采收率的影响；
7)研究注水速度对产油量和采收率的影响；
8)研究油藏平面性质和层间非均质性对油藏动态的影响；
9)验证油藏的面积和地质储量；
10)检验油藏数据；
(二)数据准备
（1）数模所需的数据类型
1)网格数据。
网格数据包括网格步长、油藏顶面海拔、每个油层的总厚度、有效厚度、渗透率、孔隙度、初始饱和度、初始压力和岩石类型，除网格步长、岩石类型和初始压力外，其他数据都来自测井解释，不同的岩石类型提供不同的相对渗透率和毛管压力数据。
2)表格数据。
包括油气pVT数据表(来自高压物性分析)；水和岩石性质(也可借用标准参考值)；油水相对渗透率曲线也来自室内实验；毛管压力曲线来自室内的岩心实验；井筒流动数据来自井筒参数计算。
对于一个黑油油藏，pVT数据是非常重要的。pVT数据是由地层体积系数、溶解油气比和粘度作为压力的函数表所组成的。
各项表格数据可以不止一组，关键是变量与自变量之间的关系要光滑。
3)动态数据。
动态数据是一切与时间有关的数据，包括：
完井数据：射孔、补孔、压裂、堵水、解堵日期、层位、井指数等。
生产数据：平均日产油、日产水、日产气、平均气油比和含水比等。
压力数据：井底流压、井网格压力等。
还有累积时间和输出控制等参数。
4)其他参数。
这些参数包括算法选择和使用，输入输出控制，油水井约束界限，油井的定压、定产参数等。
模型中油井数据是由产量、井底压力和附加于每口井的约束条件所组成的，如果规定油井产量，则产水量和产气量将根据井节点的相对渗透率曲线算出，如果规定了井底压力，那么产量可以根据所维持的压力来计算。
（2）数据准备要点
1)建立网格参数场。
各大数值模拟软件(如：WorkBench、DESKTOP—VIP、CMG等)都有前后处理软件，可用它们自身的前处理模块建立各种油藏性质的网格参数场，也可以用油藏描述结果来建立网格参数场。如何用这两种方式来建立网格参数场，将在后面详细介绍。但无论用哪种方式，首先都要建立合理的网格系统，建立网格系统。建立网格必须考虑以下几个方面的问题：
①选择合理的网格系统。
在前面已叙述了网格的发展状况。如何选择网格系统，一是要根据实际油藏的几何形态和地质特征，如果油藏的几何形态和地质结构都很复杂，则选用垂直平分网格；如果油藏的几何形态和地质特征都很简单，则选用直角坐标。二是要根据研究的问题及目的，对于大的油藏模拟问题，仅仅在油藏中饱和度或压力变化剧烈的区域需要细网格，使用标准的均匀网格将导致油藏某些部位不希望出现的小网格，在这种情况下，使用局部网格加密；对于研究锥进问题或井周围要获得精确计算等，可用局部网格加密，也可用混合网格。
②网格定向。
网格的边界要与天然的非流动边界相符合；应包含有效的井位；考虑液体流动的主要方向和油藏内的天然势能梯度；考虑油藏性质的方向；尽量减少死网格节点数。
③网格尺寸。
选择网格尺寸，要根据实际模拟对象、目的及本单位计算机的情况而定。若只是进行一些机理性研究，可选择较大尺寸的网格；若要作精细油藏数值模拟，
提供加密井井位等研究时，就要用较小的网格尺寸，尤其是加密井周围，要用较小的网格；若要作非常大的油藏数值模拟研究，考虑到计算机的内存及计算速度，一般用较粗的网格。
2)建立网格属性参数场。通过内插得到属性数据
根据各小层在参数井点处的属性参数值，内插出各节点处的属性值(如孔隙度、渗透率、油藏顶部深度、有效厚度、砂岩厚度、初始饱和度、初始压力等)，插值方式有最近井插值、确定性距离反比加权平均插值、统计性加权平均插值、克利金插值等。
3)建立表格数据。
用油藏数值模拟软件的前处理模型建立各种表格数据(如平衡区、相对渗透率、毛管压力、油气水的高压物性、输出区域等)
4)建立动态模型。
可用数模软件的前处理模型或手工编辑等方式建立动态模型，如Eclipse软件可用Schedule软件建立动态模型；WorkBench软件可用“生产数据分析”软件建立动态模型；DESKTOP—VIP软件用Prexec建立动态模型。网格数据、表格数据和动态模型就构成了模拟模型。
（三）模型初始化
在这一阶段，主要是对输入数据进行检查，并对油藏储量进行拟合，其步骤是：
1.软件依次读入及处理各部分数据。在开始读入下部分内容前，必须对本部分各种数据进行一致性检查。只有最后一部分例外，因为该部分数据是时间相关的，因此，总体上初始化时并不读入与处理。
2.软件开始运行的第一个任务是为输入的数据分配内存。
3.处理模拟网格的几何情况及属性，得到更利于流动计算的三维空隙体积、传导率、单元中部深度及可产生流动的网格间的连接。
4.在初始化时，计算油藏每一层的静压力剃度，并给每一个网格赋每一相的饱和度值。
5.储量拟合。
（四）生产史拟合
用已知的地质、流体性质和特殊岩心分析资料和实测的生产历史(产量或井底压力随时间变化)，输入计算机程序中，将计算结果与实际观测和测定的开发指标(油层压力和综合含水率等)相比较。若发现两者间有相当大的差异，则说明所用的资料与实际油田资料差异很大，逐步修改输入数据，使计算结果与实测结果一致。
(1)基本概念。
所谓历史拟合，就是用已有的油藏参数(如渗透率、有效厚度、孔隙度、饱和度等)去计算油田的开发历史，并将计算的开发指标(如压力、产量、含水率等).
与油田开发实际动态相对比。若计算结果与实测不一致，则说明对油田的认识还不清楚，输人参数与地下情况不符，必须作适当调整，修改后再进行计算，直到计算结果与实际动态相吻合或在允许的误差范围内为止。通过历史拟合，可以比较客观地认识油田的过去和现状，为开发动态预测打下基础。
(2)历史拟合的步骤。
历史拟合是个复杂的、消耗人力和机时的工作，如不遵循一定步骤，可能陷入难以解脱的矛盾之中。一般认为，同时拟合全区和单井的压力、含水和气油比是没有希望的。必须将历史拟合过程分解为相对比较容易处理的步骤。历史拟合采取以下五个步骤：
1)确定模型参数的可调范围，
2)对模型参数全面检查，
3)全区和单井压力拟合，
4)全区和单井饱和度拟合，
5)单井生产指数拟合。
(3)确定模型参数的可调范围。
在历史拟合过程中，由于模型参数数量多，可调的自由度很大，而实际油藏动态数据的种类和数量有限，不能够惟一确定油藏模拟模型的参数。为了避免或减少修改参数的随意性，在历史拟合开始时，必须确定模型参数的可调范围，使模型参数的修改在合理的、可接受的范围内，这是历史拟合的原则。
确定参数的可调范围是一项重要而细致的工作，需要工程师和地质师共同努力，需要收集和分析一切可以利用的资料。首先分清哪些参数是确定的，即准确可靠的，哪些参数是不确定的，即不准确、不可靠的，然后根据情况，确定可调范围。确定的参数一般不允许修改，或只允许少量修改，不定的参数则允许修改，或在较大范围内修改。
1)孔隙度。如果油层大量岩心分析资料表明，油层部分孔隙度在19％到2％之间，平均为±20％，变化范围不大。则把孔隙度视为确定参数，不做修改，或允许改动范围在±3％。
2)渗透率。渗透度在任何油田都是不定参数。这不仅是由于测井解释的渗透率值和岩心分析值误差较大，而且根据渗透率的特点，井间的渗透率分布也是不确定的。因此对渗透率的修改，允许范围较大，可放大或缩小2～3倍或更多。
3)岩石与液体的压缩系数。液体的压缩系数是实验测定的，变化范围很小，认为是确定的。而岩石的压缩系数虽然也是实验室测定的，但受岩石内饱和液体和应力状态的影响，有一定变化范围，而且与有效厚度相连的非有效部分，也有一定孔隙和流体在内，在开发过程中也起一定弹性作用。考虑这部分影响，允许岩石的压缩系数可以扩大一倍。
4)初始流体饱和度和初始压力。和通常做法一样，认为是确定参数。必要时允许小范围内修改。
5)相对渗透率曲线。由于油藏模拟模型的网格粗，网格内部存在严重非均质，其影响不可忽视，这与均质岩心的情况不同。因此相对渗透率曲线应看作是不定参数。在拟曲线的研究中，给出了较好的初始值，但仍允许做适当修改。
6)油气的PVT性质。视为确定参数。
7)油水界面，在资料不多的情况下，允许在一定范围内修改。
(4)模型的数据检查。
油藏模拟模型的数据很多，也许会有几万至十几万个数据，出现错误的可能性大，甚至是不可避免的。在正式进行历史拟合之前，对模型数据进行全面检查是十分必要的。
数据检查包括模拟器自动检查和人工检查两个方面，缺一不可。
模拟器自动检查包括：
1)各项参数上下界面的检查，如：
①检查原始地质储量并跟容积计算进行比较；
②检查所有原始油藏的性质和输人数据。
发现某一参数超过界限，打出错误信息。
2)平衡检查。在全部模型井的产率(或注入率)都指定为零的情况下，进行一次模拟计算，模拟的时间应大于或等于油藏已经开发的时间(或历史拟合的时间)加上动态预测的时间。经过这么长时间的模拟，油藏状态参数(压力场和饱和度场)应该没有任何明显变化。这与油藏初始状态参数应是稳定状态的假定是一致的。如果发现状态变量发生了明显变化，则表明模拟参数有了问题。这是不允许的。
模拟器的自动检查是重要的，可以减少人工检查的工作量，但是不能代替人工检查。它不能发现所有的错误。人工检查可以把全部参数打印出来，进行肉眼检查，与原始数据核对。也可以用一些可视化软件将参数绘成网格块状图、等值线、曲线等进行检查。
(5)压力拟合。
压力拟合分两步，首先是拟合全区压力，然后拟合单井压力。
1)全区压力拟合方法。
①修改岩石压缩系数、孔隙度、厚度参数，增加或减少压力异常带的储量。改变孔隙度声值，可以改变地质储量。如果油藏压力水平过高，则往往表示油藏地质储量过高，此时若减少孔隙度，就可降低地质储量，而使压力水平降低；如果油藏压力水平过低，则相反，相应增加孔隙度值。
改变油层厚度、原油饱和度或压缩系数，也可改变地质储量。以上几个参数可以同时进行修改，若降低高压带的压力，则同时把上述几个参数减小，若增加低压带的压力，则同时把上述几个参数增加。
②修改渗透率值，改变流体流动方向，从而改变油层压力。
改变渗透率可以使流体从高压带向低压带流动，这样，可以改变压力异常区，改变水区渗透率，也可以改变水区附近的压力异常带的压力。
一般采用对低压带增加渗透率，对与低压带相连的水体增加水区渗透率或增加水体体积。
③检查原始地层压力梯度、原油体积系数、脱气油密度，校正地层压力水平。
原始地层压力与深度关系的准确性直接影响整个地层压力的水平。如果油藏压力水平普遍偏高或偏低，首先要检查输入的基准面深度和相应的油相压力，是否符合压力梯度关系如果不符，则需要修改。其次是从压力梯度数据中校核地下原油密度。压力梯度Dp／dD =(P2-P1)/(D2-D1),通过单位换算为地下原油密度，如果ρo有误差，则需要修改。第三,是检查原油的体积系数。检查地面脱气原油的密度，如果不正确，则按B0=ρ地面/ρ地下，求出原油的体积系数，进行修改。
总之，先拟合压力水平，后拟合压力形式。
2)单井压力拟合方法。
全区压力拟合后，注意力转移到单井压力拟合上。单井压力拟合主要靠修改井局部地区的渗透率或方向渗透率。如果邻近的井都做了类似的修改，则井间地区的渗透率也随着做相应的修改。
以上两个步骤不是截然分开的，在进行全区压力拟合时，也考虑单井情况，附带做局部修改。
压力拟合中有时也使用虚拟井。虚拟井只有在遇到特殊困难，其他办法无效时才使用，这时应尽量查明工程上或地质上的原因。例如作为模拟模型的边界——断层可能不密封，模型内的流体可能向断层的另一侧漏失或从另一侧得到补充。这时就在断层附近设虚拟井反映这种影响。又如固井质量不好，可能有流体来自外层系，则在井网格内增设虚拟井。一般情况下，应尽量避免使用。因为虚拟井的未来动态不好估计，不利于动态预测。
压力拟合时，同时要照顾到含水拟合的情况，有时可从含水拟合情况得到某些启示，助修改方向渗透率。
(6)饱和度拟合。
压力拟合达到满意的效果后，注意力转到饱和度拟合，即含水和气油比的拟合。油藏中流体饱和度的分布，影响井的注入量和产量，从而也影响油水比和气油比。
饱和度拟合也分全区和单井拟合两步。
1)全区饱和度拟合方法：
主要是修改拟相对渗透率(有时包括油水界面)。全区拟合基本满意后，再作单井拟合。
2)单井含水率拟合方法：
①调整相对渗透率数据。
如果使用实验室相对渗透率数据计算的结果提示含水比偏高，应从两个方面去检查：一个是输入实验室相对渗透率数据是否平滑，先平滑后，再运转模型。另外，用拟相对渗透率曲线代替实验室相对渗透率曲线，拟相对渗透率曲线可用Hearn方法计算或用J．R．KYTE等人的方法计算。如果拟合还有问题，可直接调整水的相对渗透率数据。
②调整毛管压力曲线。
拟合油田含水还需检查毛管压力曲线，改变束缚水饱和度，可改变初始含水饱和度分布。拟合时还应注意使用拟毛管压力曲线。
③在局部地区含水拟合相差较大，可调整渗透率
在调整局部地区含水时也有两种方法：一是改变含水区地质储量，调整这些地区的孔隙度、渗透率，即可改变含水饱和度；二是控制含水，对需要调整地区调整渗透率增加流动性，可不增加地质储量，或控制边水推进，减少与水区连通部位的渗透率或传导率值，也可达到控制含水上升的目的。
在这个阶段花费的时间可能要很长，在拟合含水饱和度(或拟合含气饱和度)的时候也可能出现油藏压力拟合的问题，所以还要兼顾油藏压力拟合。在上述压力和饱和度拟合基本完成后，拟合单井的流动压力。
(7)生产指数的拟合。
历史拟合的最后一步是单井生产指数拟合，即调整各个井的井指数，使模型
计算的井底流压与井的实际流压相拟合。流压偏小，可增大采油指数；反之，可减小采油指数。生产指因为只有最后一个时间步的生产指数对预测有意工作。单井拟合的好坏，关系到模拟质量的高低，预测方案的成败，所以搞好油藏模拟除建立地质模型(包括参数处理)比较重要之外，最重要的是单井的历史拟合。一个模拟方案不搞历史拟合就无法检验模型的可信程度，而一个历史拟合的模拟方案只搞全区拟合，不搞单井拟合，就降低了模拟方案的可靠程度，这样的模拟也只是油田指标的概算，不能寻找剩余油饱和度的分布地区，也不能进行单井调正。匡，即单井历史拟合的符合率占全油田油井总数的百分数。拟合率多少为合格?根据经验，在油田的主要开采地区的油井都应该做到拟合，或者从油田范围来讲，占全油田总油井75％以上的井拟合比较好，为合格。这个数字基本上能控制整个油、气藏的饱和度的变化。油井不可能做到100％井拟合上，如果是100％油井拟合，也可能是不真实的。个别井数据不准，地质情况特异，拟合有些相差是可以放弃的。
怎么叫拟合上了?在这方面没有一个严格的、科学的规定。一般根据现场资料的准确程度而定。比如：测地层压力的压力计误差0．5—1．0大气压，拟合标准就不能小于1．0大气压，气油比测定，有的矿场相差20％～30％，拟合标准要求过于严格，也就失去它的意义讳的是为了拟合而拟合，单井拟合的目的是为了修改油层参数，使模型参数通过动态数据修改，比较符合实际油、气藏的真实情况，就是比较准确描述地下油、气、水的分布。如果单纯调整井眼数据，比如搞井函数的特殊处理，其结果很快就会使油井全部拟合，曲线重叠远远好于上述所讲的拟合标准，误差极小。这种井函数的处理只是在井筒网格起作用，远离井筒的网格数据并没有什么改变，地下的饱和度仍然没有拟合上，所以这种拟合是不可取的。
（五）生产动态预测
在历史拟合的基础上，对未来的开发指标进行计算
应用油藏数值模拟技术决不是以再现油田开发历史为目的，油田的开发历史已成为过去。开发过程的不可逆性决定了再现开发历史只不过是为了加深对油藏现状的认识，总结以往的经验教训，确定剩余油分布规律等。而油藏数值模拟的主要任务是预测油田的未来，制定最佳开发方案或调整方案，获得最佳的效益。
新区和老区的方案预测有所不同。新区的方案预测主要是对开发方案的可行性评价、井网层系模式研究、开采方式、注水方式的选择，通过对比不同速度、不同井网模式、不同方式对开发方案的影响，提供合理的开发方案。老区的方案预测主要是在历史拟合的基础上，制定一些调整方案，通过预测计算，提供合理的调整方案。
(1)新区开发方案预测。
1)模拟初始压力和含水饱和度。
建立地质模型、核实储量，模拟地层压力分布，模拟原始含水饱和度分布。
2)开发指标预测。
利用天然能量开发预测，注水开发预测，不同井距、不同井网、不同层系、不同开采速度预测。
3)对比不同开发方案。
经济效果分析，不同方案敏感性分析，选择合理的开发方案。
(2)老区开发方案预测。
1)维持地层压力。
油藏工程师根据动态分析，制定产量预测的方案。在模拟预测阶段，注水井采用给定维持的地层压力，通过地层压力来控制和选择注入量。油井给定流压，调整产液水平，也可计算相应的地层压力，或者给定产液水平，调整流压，来维持地层压力。
2)加密井或局部新井投产。，
根据方案部署新井，可以对比一次加密和局部加密投产的开发效果，新井的工作制度按地层条件，给定产液量，再限定一个合理流压，这样，新井按产液量生产，如果达不到产液量，则会按流压限制，产液量自动调整，这样会得到比较满意的预测效果。
3)单井调整。
油田调整除了对带有气顶的油藏注意油气界面的调整外，大部分注水开发的油田都是调整油井含水，如何调整油井见水、高含水，一般分为两类：一类为不改变现有井(油田)开发方式，不钻新井的开发调整。
改善油、水井的水动力学特征(开采层段补孔、压裂、酸化、表面活性剂处理)，在模拟模型上主要是调整开采层的Kh值。
封堵油井高含水层，在模型上暂时关闭高含水层段。
改变工作制度，强采或周期性采油，在模型上也比较好解决，主要是增大生产压差或改间歇开采等。
改变注水井工作制度，调整液流方向，增加或限制注水，根据高含水层突进方向进行平面调整，这个过程在模拟中是比较困难的，要仔细分析油井的来水方向和层段(原因)，经多次调整、试算，反复改变注水井的工作制度，才能达到较好的效果。
对多油层实行细分，改变出油、出水剖面，如果不重新划分模型层的网格，在模拟计算中一般是比较容易实现的。
另一类是通过局部或全面的改变注水方式和钻新井(油、水井)调整油藏的开采过程。开发部署有重大改变，钻后备井。
通过钻新井和老井重新配产配注(老井转注等)，调整注水动态。
这个调整也含有前一类的调整，在老井单井调整的同时，要注意新井的动态调整新井，新井与老井调整相一致。
4)采收率预测。
也要随时模拟预测不同开发方案的最终采收率，并结合一般常用的油藏工程计算采收率的方法进行对比，数值模拟模型能比较方便地模拟计算各种开发方案的采收率。但有一个问题值得讨论，就是最终采收率的含义如何确定，一般油藏工程师习惯使用油田综合含水98％时的采收率为最终采收率，实际上这个数值是比较理想化的，毫无意义，因为当油田生产低于开采经济极限时，这个油田长期高含水、低产能已毫无价值，所以预测油田最终采收率的含水界限，应采用经济分析求得油田最大的含水率，来对比不同开发方案的最终采收率。
名词解释 油藏模拟油藏数值模拟数学模拟物理模型数值模型质量守恒定律适定问题初始条件黑油模型组分模型网格节点块中心网格点中心网格离散化有限差分法显示差分 隐式差分前差分后差分中心差分点交替排列格式交替对角排列格式标准排列格式 对角排列格式隐式差分格式差分方程稳定性截断误差松弛法IMPES方法历史拟合 动态预测灵敏度实验 选择题 由于油藏各点的渗透率不同，束缚水饱和度不同，因而需要对相对渗透率曲线进行归一化处理 以X方向为例，传导系数为 块中心网格是用（）来表示小块坐标的 A网格块中心B节点C网格块边缘D网格块夹角 下述表达式表示定产量内边界条件的是 认识油田的主要方法有直接观察法和模拟法 相对渗透率取值一般取上游权的处理方法 IMPES方法是（）的求解方法 A隐式压力B隐式饱和度C全隐式 历史拟合在含水拟合时主要是对（）的修改 A孔隙度B相对渗透率曲线C渗透率D地层厚度 在隐式差分格式中，有多个未知数，当已知第n时刻的值P i n时，为了求出第n+1时刻的P i n+1，需要（） A解n个方程B解一个线性代数方程组C直接求解D解一个方程 根据每一组分的质量守恒建立的渗流数学模型称为（）模型 A热采B化学驱C黑油D组分 一维径向模拟时r=10cm，r=40cm，那么可以推断r s的大小是 A120 B200 C400D 640 下列哪一种方法不属于迭代求解方法 A雅克比法B超松弛法CLU分解法D交替方向隐式法 对于二位6*4网络系统，如果按行标准排列，气半带宽W= A6 B4 C12 D8 克兰克?尼克森差分格式的截断误差为（） 块中心网格和点中心网格的差分方程相比较，结果（） A一样的B有半个网格的误差C相差流动项系数D维数不同 三．判断题2分*10 1.黑油模型中水相与其他两相不发生质量转移，气可以从油中出入，但不能汽化液相 2.离散化的核心是把整体分为若干单元来处理，它是油藏对象的空间离散 3.显式差分格式是有条件收敛的 4.差分方程组的直接解法的特点是计算工作量小，精确度较高，计算程序简单 5.差分方程组的迭代解法主要用于处理系数矩阵阶数较高的问题 6.相对渗透率取值一般取上游权的处理方法 7.油藏模拟的基础在于油藏描述和生产动态，若油层参数和生产数据不准确，通过数值模 拟的算法也可以消除 8.显示差分格式的稳定条件是△t/△x2≤0.5 9.有限差分法就是用差商来代替微商
油藏数值模拟 油藏数值模拟是随着电子计算机的出现和发展而成长的一 门新学科，在国内外都取得了迅速的发展和广泛的应用。 1953年美国G..H.BUCE等人发表了《孔隙介质不稳定气体渗流的计算》后，为用数值方法计算油气藏渗流问题开辟了道路。三十多年来，由于大型快速电子计算机的迅速发展，大大地促进了数值模拟方法的广泛应用。20世纪60年代初期研究了多维多相的黑油模型；20世纪70年代初期研究了组分模型、混相模型和热力采油模型；20世纪70年代末期研究各种化学驱油模型。目前，黑油、混相和热力采油模型已经投入工业性应用，并已经成为商业性软件，化学驱油模型也正日趋完善。 油藏数值模拟方法是迄今为止定量地描述在非均质地层中 多相流体流动规律的惟一方法。例如许多常规方法要假定油层为圆形的均匀介质，如油藏几何形状稍复杂一些，且为非均质介质，则求解非常困难，甚至无法求解。而对油气藏数值模拟而言，计算形态复杂的非均质油藏和计算简单形态的均质油藏工作量几 乎是一样的。因此油藏数值模拟可解决其它方法不能解决的问题。对于其它方法能解决的问题，用数值模拟方法可以更快、更省、更方便、更可靠地解决，并增加其它分析方法的可信度。 一个油气藏，在现实中只能开发一次。但应用油藏数值模拟，可以很容易地重复计算不同开发方式的开发过程，因此人们可以从中选出最好的开发方法。
因此，对油藏工程师而言，数值模拟给动态分析提供了一种快速、精确的综合性方法；对管理者而言，数值模拟提供了不同开采计划的比较结果；对尚无经验的工程师而言，数值模拟则是有效的培训工具。 数值模拟研究的主要工作程序对一个油气藏进行综合的数模研究，往往需要花较大的精力和较长时间（有时会达一年甚至更长的时间），同时还对计算机硬件和技术人员有很高的要求，然而尽管在不同的项目中，面对的问题会千差万别，但大多数油藏数值模拟的基本研究过程是一样的。为了使读者一开始就对数模研究工作有一个明确的整体概念，下面简要地介绍一下油藏数值模拟的主要工作程序。 问题的定义：开展油藏数模工作的第一步，是确定研究的目标和范围。即首先要给本次数模研究一个明确的定位，明确本次模拟要解决的主要问题是什么，需要研究哪些油藏动态特性，这些项目的完成对油藏的经营管理者会产生什么影响等等。从而根据项目的要求进行数值模拟研究程序设计，并收集有关的油藏基础地质、流体及生产动态数据。 数据的检查：一旦把数据收集起来以后，必须对这些来自不同渠道的数据进行鉴别，再组织和再检查，看收集到的数据是否足够，是否都合格。如果取得的数据，依靠经验和评价方法进行修正和补充后仍不合要求，那就需要修正或重新确定研究目标。
通过了几节课的“油藏数值模拟课”的学习，我知道了“油藏数值模拟”是应用计算机研究油气藏中多相流体渗流规律的数值计算方法,它能够解决油气藏开发过程中难以解析求解的极为复杂的渗流及工程问题,是评价和优化油气藏开发方案的有力工具。它主要是让我们石油石油工程专业的学生掌握一些基本的油藏数值模拟技术和技巧,学习基本的油藏渗流数学模型及其解法、计算方法和应用方法,培养我们用计算机解决油藏开发问题的能力。 “油藏数值模拟”涉及的学科较多,利用数学知识和计算机知识较多,我认为是非常难的。虽然教师教的很认真也很耐心，我仍然不能跟着老师的节奏。因为一开始就知道这个软件很有实际应用价值，所以我也就特别的想好好的学习它。可惜现在我面临着考研这座大山，我实在是没有充分的时间课下来好好的温习与研究老师上课所讲的东西。很遗憾，后来老师讲的东西我有些就不会了。好在前三四节课讲的内容还学会了，学会了模拟三层的油层概况。也许这点知识对我以后的再次学习会有不错的基础作用吧！总之还是很感谢老师的耐心教导。 在学习的过程中，我觉得油藏原始参数，如渗透率、孔隙度等的收集，以及油藏原始数据是否齐全准确非常重要，尤其是一开始填date时的单位的选择，这些都关系到数值模拟的效果。如果原始资料很少，数值模拟的效果就不可能好。数值模拟方法越复杂，所需的原始资料也越多。收集资料时，如发现必需的资料不够或不准确，应采取补救措施。通常要求准备的参数包括：①油藏地质参数。产层构造图，油、气、水分布图，油层厚度、孔隙度、渗透率、原始含油饱和度的等值图等。②流体物理性质参数。地面性质和地层状态下的物性数据，原始压力和地层温度数据，对凝析气田还需要相图和相平衡的资料。③专项岩心分析资料。油水相渗透率曲线，油气相渗透率曲线，油层润湿性，吸入和排驱毛细管压力曲线；对碳酸盐岩孔隙裂缝双重介质储层，还需渗吸曲线。④单井和分层分区的生产数据和有关测试资料。⑤油田建设和经济分析的有关数据。 将收集的油藏地质资料进行系统整理后，要将油藏的地质特征模式化，以充分反映油藏的构造特征和沉积特征，如油层物理性质参数的分布、油气水的分布、油气水在地面和地下的性质、驱油动力、压力系统和地温梯度等。油藏地质模型是否符合实际情况，直接影响数值模拟成果的准确性。 由于人们对油田实际地质条件的认识有一定的限度，计算时所用的参数也就有一定的局限性，因此，第一次模拟计算的结果，如压力、产量、气油比、含水率等与油田实际生产状况常有较大的出入。必须进行分析，修改相关的计算参数，重新进行计算。通常，经过多次修改可使计算结果与实际生产历史基本相符，误差在允许范围以内。从工程应用的角度看，可认为此时所应用的计算参数，反映了油田地下的实际状况，使用这些参数来计算和预测油田未来的动态，能够达到较高的精度。在油田开采过程中这类历史拟合要进行多次，使油田的模型逐步更接近实际而得到更适用的结果。
第一章油藏数值模拟方法分析 油藏数值模拟 油藏数值模拟简述 油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律 的数学模型，并利用计算机求得数值解来研究其运动变化规律。其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。其基础理论是基于达西渗流定律。 油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模 拟实际的油田开采的一个过程。基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值（生产动态）与实际吻合。其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。这组流动方程组由运动方程、状态方程和连续方程所组成。油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。 油藏数值模拟技术从50年代的提出到90年代间历经40年的发展，日益成熟。现在进入另外一个发展周期。近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用 非常广泛。 油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模 拟研究，且可重复、周期短、费用低。
图1油藏数值模拟流程图 油藏数值模拟的类型 油藏数值模拟类型的划分方法有多种， 划分时最常用的标准是油藏类型、 需要模拟的油 藏流体类型和目标油藏中发生的开采过程， 也可以根据油气藏特性及开发时需要处理的各种 各样的复杂问题而设定， 油气藏特性和油气性质不同， 选择的模型也不同， 还可以根据油藏 数值模拟模型所使用的坐标系、空间维数和相态数来划分。 以油藏和流体类型来划分，其模型有：气体模型、黑油模型和组分模型； 以开采过程来 划分，其模型包括：常规油藏、化学驱、热采和混合驱模型。 以油藏和流体描述为基础的油藏模型分为两类：黑油模型和组分模型。 （1） 黑油模型，是常规油田开发应用的油藏数值模型，用于开采过程中，对油藏 流体组分变化不敏感的情况， 是最完善、最成熟的。黑油模型假设质量转移完全取决于 压力变化，适应于油质比较重的油藏类型，在这些模型中，流体性质 E O 、B g 、R S 决定PVT 的 变化，如普通稠油及中质油的油气藏。 （2） 组分模型，应用于开采过程中对组分变化敏感的情况。这些情况包括：挥发性油 藏和凝析气藏的一次衰竭采油阶段， 用组分模型进行模拟。在组分模型中，适用于油质比较轻、气体组分比较高的油气藏， 使用 数据化 流体的PVT 数据、相 渗曲线、岩石数据 建立地质模型 建立网格 参数场 表格数据 油水井产量、井史 数据 T 动态模拟 含油边界拟合 非井点地质静态参数拟合 区块、单井压力拟合 生产指数拟合 以及压力保持阶段。同时，多次接触混相过程通常也采
[转]【推荐】油藏数值模拟入门指南 尝试写一写油藏数值模拟入门指南，希望对那些刚刚开始进入油藏数值模拟领域的工作者有所帮助。 第一：从掌握一套商业软件入手。 我给所有预从事油藏数值模拟领域工作的人员第一个建议是先从学一套商业数值模拟软件开始。起点越高越好，也就是说软件功能越强越庞大越好。现在在市场上流通的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的集成平台。 对于初学者，不但要学主模型，也需要学前后处理。对于ECLISPE的初学者，应该先从ECLISPE OFFICE学起，把ECLISPE OFFICE的安装练习做完。然后再去学Flogrid，Schedule 和SCAL。PVTi主要用于组分模型，做黑油模型可以不用。 第二：做油藏数值模拟都需要准备什么参数 在照着软件提供的安装例子做练习时经常遇到的问题是：虽然一步一步按照手册的说明做，但做的时候不明白每一步在做什么，为什么要这么做。这时候的重点在于你要知道你一开始做的工作都是为数值模拟计算提供满足软件格式要求的基础参数。有了这些基础参数你才能开始进行模拟计算。这些基础参数包括以下几个部分： 1。模拟工作的基本信息：设定是进行黑油模拟，还是热采或组分模拟；模拟采用的单位制（米制或英制）；模拟模型大小（你的模型在X,Y,Z三方向的网格数）；模拟模型网格类型（角点网格，矩形网格，径向网格或非结构性网格）；模拟油藏的流体信息（是油，气，水三相还是油水或气水两相，还可以是油或气或水单相，有没有溶解气和挥发油等）;模拟油田投入开发的时间；模拟有没有应用到一些特殊功能（局部网格加密，三次采油，端点标定，多段井等）；模拟计算的解法（全隐式，隐压显饱或自适应）。 2。油藏模型：模型在X,Y,Z三方向的网格尺寸大小，每个网格的顶面深度，厚度，孔隙度，渗透率，净厚度（或净毛比）。网格是死网格还是活网格。断层走向和断层传导率。
数值模拟的目的 (一)、为什么开展油藏数值模拟工作 研究和开发一个油田是一个复杂的综合性的科技问题，高精度的地震资料的处理解释提供研究区域的构造、断层、边界及其走向，但地震纵向分辨率受到限制，不能很好的反映一个同相轴(地震道) 中沉积砂体的物性变化特征；测井可较好的反映到小于1米以下沉积砂体的物性特征，提供可靠的地层对比结果。但作为新老油田开发方案的研究及剩余油分布的研究，是地震、地质、测井理论方法都无法做到的。地质上仅定性或半定量分析，测井用于生产监测不能以点带面。惟独油藏数值模拟工作可再现生产历史，定量分析剩余油潜力；并做到室内研究投入少、时间短，还可进行开发方案优选及经济评价工作。所以总公司强调开发方案的部署一定要开展数值模拟工作。值得强调的是油藏数值模拟工作提倡一体化，注重前期的地震解释和测井解释即油藏描述工作。 (二)、油藏数值模拟的目的 在进行油藏数值模拟工作前，首先应根据油田开发过程中存在难以解决的实际问题，提出开展此项工作的目的及意义，即最终所要达到解决问题的目标是什么？一般通过油藏数值模拟可进行以下研究工作： 1. 初期开发方案的模拟 1) .评价开发方式；如：枯竭开采、注水开发等。 2) .选择合理井网、开发层系、确定井位； 3) .选择合理的注采方式、注采比； 4) .对油藏和流体性质敏感性研究。 2. 对已开发油田历史模拟 1) . 核实地质储量，确定基本的驱替机理(如：是天然驱，还是注水开发。)； 2) .确定产液量和生产周期； 3) .确定油藏和流体特性； 4) .提出问题、潜力所在区域。 3. 动态预测 1) .开发指标预测及经济评价 2) .评价提高采收率的方法(如：一次采油、注水、注气、化学驱等) 3) . 剩余油饱和度分布规律的研究,再现生产历史动态诸如：研究剩余油饱和度分布范围和类型； ?单井调整：改变液流方向、注采井别、注水层位； ?扩大水驱油效率和波及系数； 4) .潜力评价和提高采收率的方向 诸如： ? 确定井位、加密井的位置；
前言： 油藏数值模拟是随着计算机的发展，而在石油行业中逐步成为一门成熟的技术。追溯油藏数值模拟的发展史，从30年代开始研究渗流力学到50年代在石油工业方面得以应用，到70年代进入商品化阶段，而80年代油藏数值模拟又向完善、配套、大型多功能一体化综合性软件飞跃发展。近十年油藏数值模拟已成为油田开发研究，解决油田开发决策问题的有力工具。在衡量油田开发好坏、预测投资、对比油田开发方案、评价提高采收率方法等方面应用都极为广泛。 油藏数值模拟就是应用数学模型再现实际油田生产动态。具体通过渗流力学方程借用大型计算机，结合地震、地质、测井、油藏工程学等方法在建立的三维地层属性参数场中，对数学方程进行求解，实现再现油田生产历史，解决油田实际问题。 油藏数值模拟是一门综合性很强的科学技术，涉及油田地质、油层物理、油藏工程、采油工程、测井、数学、计算机及系统等学科。而油藏数值模拟工作又以其繁重的前期准备和上机历史拟合运算工作让人望而生畏。 那么如何做好前期资料准备工作和尽快掌握模拟技巧？使得今后的油藏数值模拟工作在作业区顺利开展，便是出此书的目的所在。 本书结合以往工作中的实际经验教训，成功与失败，参考诸多资料从前期数据准备工作开始到模拟技巧做了较为的详细介绍，以舐读者。有不妥之处，请予指证。同时，今后不定期的将更新的模拟技术及方法推荐给大家。 目录 一、数值模拟发展概况 二、数值模拟的基本原理 二、选择适当的数值模型及相类 三、数据录取准备工作 （一）建立油藏地质模型 （二）网格选择 （三）数据录入准备 四、历史拟合方法及技巧 （一）确定模型参数的可调范围 （二）对模型参数全面检查 （四）历史拟合 附件1：关于实测压力的皮斯曼校正 附件2：关于烃类有效孔隙体积的计算 一、数值模拟发展概况 30年代人们开始研究地下流体渗流规律并将理论用于石油开发； 50年代在模似计算的方法方面，取得较大进展； 60年代起步，人们开始用计算机解决油田开发上的一些较为简单间题，由于当时计算机的速度只有每秒几万到几十万次，实际上只能做些简单的科学运算； 70 年后主要体现于计算机的快速升级带动了油藏数模的迅猛发展，大型标量机计算速度达到100--500万次，内存也高增主约16兆字节。在理论上黑油模型计算方法更趋成熟，D. W.
目录 1. 前言 (1) 上机实践的目的及要求 (1) 主要完成的实践内容 (2) 2. 油藏特征分析 (2) 储层物性特征 (2) 流体物性特征 (2) 储层岩石物性特征 (2) 气藏数值模型建立 (2) 模型网格的划分 (2) 模型物性 (3) 模型流体性质及相渗曲线 (3) XX气藏地质储量 (3) 4. XX气藏方案优选 (3) 开发方案的优选 (3) 采速与稳产时间的关系 (4) 5. 结论认识 (4) 结论 (4) 对本实践课程的建议 (4) 1. 前言 上机实践的目的及要求 1. 掌握油藏数值模拟的上机操作流程； 2. 掌握ECLIPSE软件的数据录入、编辑和修改方法； 3. 掌握ECLIPSE软件结果输出及三维可视化方法；
4. 掌握机理模型研究方案设计的思路及方法 主要完成的实践内容 1. 油藏数值模拟数值整理； 2. 依据现有数据，应用块中心网络系统建立一个三维油藏数值模拟模型； 3. 预测单口气藏天然能量开发的最终采收率（20年）（不考虑水体能量）； 4. 预测多口气井采收率（20年）； 5. 预测不同稳产年限下，气井的合理产量（稳产5年）； 6. 水平井开发和直井开发效果对比； 2. 油藏特征分析 储层物性特征 表2-1 储层物性特征 流体物性特征 气藏数值模型建立 模型网格的划分
模型流体性质及相渗曲线 XX气藏地质储量 4. XX气藏方案优选开发方案的优选 水平井方案
水平井方案 采速与稳产时间的关系 采油速度越快，稳产时间越短。采油速度越慢，稳产时间越长。由此可见采油速度与稳产时间成反比。 5. 结论认识 结论 通过这个实验，我们了解了eclipse软件的基本操作，并且建立了一个简单的均质油藏的模型，并且成功计算了产量。这个实验然我们获益匪浅。 对本实践课程的建议 建议增加实验课的课时，其余的方面都很好。老师讲的不错，需要学习的内容都学会了。
第一章油藏数值模拟方法分析 令狐采学 1.1油藏数值模拟 1.1.1油藏数值模拟简述 油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。其基础理论是基于达西渗流定律。 油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值(生产动态)与实际吻合。其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。这组流动方程组由运动方程、状态方程和连续方程所组成。油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层
模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。 油藏数值模拟技术从50 年代的提出到90 年代间历经40 年的发展，日益成熟。现在进入另外一个发展周期。近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用非常广泛。 油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。 图1 油藏数值模拟流程图 1.1.2油藏数值模拟的类型 油藏数值模拟类型的划分方法有多种，划分时最常用的标准是油藏类型、需要模拟的油藏流体类型和目标油藏中发生的开采过程，也可以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定，油气藏特性和油气性质不同，选择的模型也不同，还可以根据油藏数值模拟模型所使用的坐标系、空间维数和相态数来划分。 以油藏和流体类型来划分，其模型有：气体模型、黑油模型和组分模型；以开采过程来划分，其模型包括：常规油藏、化学驱、热采和混合驱模型。 以油藏和流体描述为基础的油藏模型分为两类：黑油模型
第一章油藏数值模拟方法分析 1.1油藏数值模拟 1.1.1油藏数值模拟简述 油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。其基础理论是基于达西渗流定律。 油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值(生产动态)与实际吻合。其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。这组流动方程组由运动方程、状态方程和连续方程所组成。油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。 油藏数值模拟技术从50 年代的提出到90 年代间历经40 年的发展，日益成熟。现在进入另外一个发展周期。近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用非常广泛。 油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。 图1 油藏数值模拟流程图 1.1.2油藏数值模拟的类型 油藏数值模拟类型的划分方法有多种，划分时最常用的标准是油藏类型、需要模拟的油藏流体类型和目标油藏中发生的开采过程，也可以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定，油气藏特性和油气性质不同，选择的模型也不同，还可以根据油藏数值模拟模型所使用的坐标系、空间维数和相态数来划分。 以油藏和流体类型来划分，其模型有：气体模型、黑油模型和组分模型；以开采过程来划分，其模型包括：常规油藏、化学驱、热采和混合驱模型。 以油藏和流体描述为基础的油藏模型分为两类：黑油模型和组分模型。 （1）黑油模型，是常规油田开发应用的油藏数值模型，用于开采过程中，对油藏 流体组分变化不敏感的情况，是最完善、最成熟的。黑油模型假设质量转移完全取决于压力变化，适应于油质比较重的油藏类型，在这些模型中，流体性质B o、B g、R s决定PVT 的变化，如普通稠油及中质油的油气藏。 （2）组分模型，应用于开采过程中对组分变化敏感的情况。这些情况包括：挥发性油藏和凝析气藏的一次衰竭采油阶段，以及压力保持阶段。同时，多次接触混相过程通常也采用组分模型进行模拟。在组分模型中，适用于油质比较轻、气体组分比较高的油气藏，使用三次状态方程表示PVT变化，如轻质油或凝析气藏。 （3）根据一些特殊开采方式的需要而形成的其他类型的数值模型，如热采模型、注聚
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 g o g rw w w w w w w 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 ρρρ=+，地面油的相对密度为地面油与溶解气相对密度之和。 3、,,og o gd o o gd gd g g γγγγργρ=+== 辅助方程： 饱和度（三相）1o g w S S S ++= 饱和度（两相）1o g S S += 毛管力（三相）：() ()o w cow w g o cog g p p p S p p p S -=???-=?? 毛管力（两相）：()o w cow w p p p S -=
油藏数值模拟技术现状与发展趋势 摘要：介绍了当前国内外油藏数值模拟的现状，简述了并行算法、网格技术、粗化技术、数值解法、动态油藏模型建立、动态跟踪模拟及三维显示等技术，指出了数值模拟的发展趋势。 关键词：并行算法；网格技术；网格粗化；分阶段模拟；动态跟踪模拟；数值解法 引言 近年来，随着计算机、应用数学和油藏工程学科的不断发展，油藏数值模拟方法得到不断的改进和广泛应用。通过数值模拟可以搞清油藏中流体的流动规律、驱油机理及剩余油的空间分布；研究合理的开发方案，选择最佳的开采参数,以最少的投资,最科学的开采方式而获得最高采收率及最大经济效益[1]。经过几十年的发展,该技术不断成熟和完善并呈现出一些新的特点。 1 国内外现状 1.1 并行算法 并行算法是一些可同时执行的诸进程的集合这些进程互相作用和协调动作从而达到给定问题的求解[2]。并行算法首先需合理地划分模块,其次要保证对各模块的正确计算,再次为各模块间通讯安排合理的结构,最后保证各模块计算的综合效果并行机及并行软件的开发和应用将极大地提高运算速度,以满足网格节点不断增多的油藏数值模型。在并行计算机上使用并行数值解法是提高求解偏微分方程的计算速度,缩短计算时间的一个重要途径[3,4]。在共享内存的并行机上把一个按向量处理的通用油藏模拟器改写成并行处理是容易的,但硬件扩充难；分布内存并行机编程较共享式并行机困难,但硬件扩充容易,关键是搞好超大型线形代数方程组求解的并行化。并行部分包括输入输出、节点物性、构造矩阵、节点流动及井筒等。 1.2 网格技术 为了模拟各种复杂的油藏、砂体边界或断层渗透率在垂向或水平方向的各向异性,以及近井地区的高速、高压力梯度的渗流状态,近年来在国外普遍发展了各种类型的局部网格加密及灵巧的网格技术。这种系统大体可以分为二类：一类称控制体积有限元网格(CVFE),这是将油藏按一定规则剖分为若干个三角形以后,把三角形的中心和各边的中点连接起来所形成的网格。另一类则称垂直等分线排比网格(PEBI),其剖分方法是将油藏分成若干三角形后,使三角形各边的垂直等分线相交而形成网格。这些方法在处理复杂几何形状油藏及进行局部网格加密时简单而一致。在多相流情况下,参照某一给定的几何准则时该方法是单调的,这保证了其稳定性和收敛性。这两种方法都能以直观的控制体积的概念出发并且采用一致的上游权而推导得出这些方法对网格的方向不敏感,在某些情况下比九点差分格式的效果好。 1.3 计算机辅助历史拟合技术
裂缝性油藏数值模拟方法 姚军 （中国石油大学山东东营 257061） 摘要：目前对天然裂缝性油藏的数值模拟可以大致分为连续性模型和离散性模型两大类；连续性模型又可以分为双重介质模型和单介质模型，双重介质模型主要是以Barrenblatt和Warren-Root在20世纪60年代提出的双重孔隙/双重渗透模型为基础，在这类模型中认为油藏中每一点都存在有基岩和裂缝两种介质，基岩被相互平行排列的裂缝分割称为单个的岩块，每种介质存在独立的水动力场，通过两种介质间的窜流的将其联系起来；而对于单介质模型，则是通过一定的方法将裂缝的渗透率和基岩的渗透率进行综合的考虑，得出整个油田的有效渗透率，该有效渗透率考虑了裂缝的密度、方位等的影响，然后将该有效渗透率输入到普通的单一介质模拟器中来对裂缝性油藏进行模拟； 由于双重介质模型不能够对不连续且控制着流体流动的大裂缝进行准确的模拟等原因，离散性模型在近段时间逐渐发展起来，而其又可以分为离散裂缝网络模型和离散管网模型；在离散裂缝网络模型中，对地质上描述出来的每个裂缝都进行了离散的显式的表示，同时根据局部裂缝的形状决定基岩的几何形状，由于地质上描述的裂缝数目一般较多，相应的在数值模拟中需要的离散点数目也就十分巨大，对模拟造成了一定的困难，所以目前很多的专家和学者又对该方法进行了进一步的改进，有许多简化的方法存在；离散管网模型则是先对所要模拟的区域进行了网格的划分，进而采用管子连接两个网格块，相应的两个网格块之间的传导率也采用管子的传导率来代替，这种方法的特点是数学上比较简单，灵活性较强，同时由于管子只对其连接的两个网格有影响，所以改变管子的传导率只会影响一个方向的传导性，而不会像常规的模拟器那样要同时影响两边的传导性，但是该方法目前研究较少。 0 前言 随着世界碳酸盐岩油气田的大规模开发，系统深入研究这类油气田的渗流模式及其在开发中的应用已成为重要课题。地质学家通过岩芯分析，确认碳酸盐岩(灰岩、白云岩)具有明显可见的裂缝、孔洞，含有密集的树枝状构造的粗裂缝以及连接的孔洞和孔隙。这类特殊的储集层结构不仅造成了井的高产、不稳定、跃变等开采特征，而且也造成各异的油气井压力降或压力恢复曲线特征。 碳酸盐岩油藏在孔隙结构和渗流机理上同砂岩油藏相比都存在很大的差别，由于天然裂缝的发育十分的不规则，裂缝的密度、长度、方位等参数都会因沉积
eclipse油藏数值模拟一些入门心得 记得上大学最早学围棋时总感觉无从入手，看身边的朋友下棋时学着聂卫平从容入定，潇洒自如的样子，很是羡慕。后来从书店买来围棋入门指南，夜深人静时照着指南慢慢学如何吃子，如何做眼，什么是打劫，怎么样布局。掌握了一点基本知识以后开始找水平最差的下，输了一定不能弃擂，脸皮要厚，缠着对方接着下。赢了水平最差的人后去找中等水平的人下。这样经过一年半载，再看以前那些学着聂卫平从容入定，潇洒自如下棋的同学，心想他们原来不过如此，赶老聂差十万八千里哪。在这里也有许多人把我叫大师，专家，如果哪一天你觉得其实我的水平也很一般，那你就到了专业段位了。 市场上有不少关于油藏数值模拟的书，但好像没有类似围棋入门指南那样从基础开始一步一步介绍的书。我收到不下二十个问油藏数值模拟如何入门的问题。我尝试写一写油藏数值模拟入门指南，希望对那些刚刚开始进入油藏数值模拟领域的工作者有所帮助。 第一：从掌握一套商业软件入手。 我给所有预从事油藏数值模拟领域工作的人员第一个建议是先从学一套商业数值模拟软件开始。起点越高越好，也就是说软件功能越强越庞大越好。现在在市场上流通的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同时也是EC LIPSE的集成平台。 对于初学者，不但要学主模型，也需要学前后处理。对于ECLISPE的初学者，应该先从ECLISPE OFFICE 学起，把ECLISPE OFFICE的安装练习做完。然后再去学Flogrid，Schedule和SCAL。PVTi主要用于组分模型，做黑油模型可以不用。 第二：做油藏数值模拟都需要准备什么参数 在照着软件提供的安装例子做练习时经常遇到的问题是：虽然一步一步按照手册的说明做，但做的时候不明白每一步在做什么，为什么要这么做。这时候的重点在于你要知道你一开始做的工作都是为数值模拟计算提供满足软件格式要求的基础参数。有了这些基础参数你才能开始进行模拟计算。这些基础参数包括以下几个部分： 1。模拟工作的基本信息：设定是进行黑油模拟，还是热采或组分模拟；模拟采用的单位制（米制或英制）；模拟模型大小（你的模型在X,Y,Z三方向的网格数）；模拟模型网格类型（角点网格，矩形网格，径向网格或非结构性网格）；模拟油藏的流体信息（是油，气，水三相还是油水或气水两相，还可以是油或气或水单相，有没有溶解气和挥发油等）;模拟油田投入开发的时间；模拟有没有应用到一些特殊功能（局部网格
精细油藏数值模拟研究现状及发展趋势 教师：王庆 学生：扎紫拉 学号：2012030012 中国石油大学（北京）
前言 油藏数值模拟是一门工程应用学科。它立足于流体在多孔介质中的渗硫理论，利用数值物理方法，通过编制计算机软件来求解油藏流体渗硫问题。油藏数值模拟技术自20世纪50年诞生至今，随计算机应用数学和油藏工程的发展而不断发展，目前已成为油田开发方案设计，动态分析和油藏开发中后期方案调整的有效工具，在各油田开发生产中得到了广泛应用。 油藏数值模拟是油藏管理环节中必不可少的部分。油藏管理的根本是在一个良好开发的油田中决定采用什么样的油气开采方式来获得最大的经济效益。要想达到油藏管理的基本目标，油藏数值模拟是最精密复杂的可行方法。开展数值模拟研究有许多原因。从商业的角度考虑，可能最重要的一点是油藏模拟可以进行大概的现金流动预测。数值模拟从产量的方面预测经济效益。将产量和价格结合起来就可以估计将来的现金流动。 表1.数值模拟的原因 油气藏是在单一圈闭中具有同一个压力系统的油气聚集单元。在原始条件下，油气藏处于平衡状态 ; 当受到干扰（如打井，生产）时，原来的平衡状态被打破，油气藏处于动态变化中。油气藏从投入开发到最后废弃就是一个不断变化的动态过程。
描述或实现油气藏动态变化的过程称为模拟（或仿真）。要描述或实现这一动态变化，可以有两种方法： (1)采用物理实体的方法，称为物理模拟; (2)采用数学描述的方法，称为数学模拟; 物理模拟是指根据同类现象或相似现象的一致性，利用某种模型来观察和研究其 原型或原现象的规律性。物理模型包括相似模型和单元模型两种。 数学模拟是指用数学模型来进行研究，即通过求解某一物理过程的数学方程组来 研究这个物理过程变化规律的方法。数学模型的核心问题是把地层流体在孔隙介质中 的渗流机制描述清楚，并施加一定的初始，边界条件，则可以相应的求出地下流体的 压力，饱和度等参数， 从而认识地下流体运动的规律。数学模型有3类：水电相似模型，解析模型，数值模型。 一.油藏数值模拟的主要内容 藏数值模拟的主要内容可以概括为三大部分，即建立数学模型（Mathematical Model ），数值模型（Numerical Model）和计算机模型（Computer Model ）。 1.建立书数学模型，就是要建立一套描述油藏中 流体渗流的偏微分方程组一起构成一格完整的数学模型。 2.建立数值模型，即对所建立的数学模型进行数值求解，一般要经过以下三个步骤： （1）离散化，将连续的偏微分方程组转化成离散的有限差分方程组，即将连续 函数变为离散函数。 （2）线性化，将有限差分方程组中的非线性系数项性化，从而得到线性代数方 程组。对线性代数方程组进行求解，包括直接求解法和迭代求解法。
-1- 一选择题 1．双模包括物理模拟和（B 数学模拟） 2．随着计算机的迅速发展，求解数学方程组常用（A 数值模拟） 3．油藏模拟的基础在于油藏描述和（A 生产动态） 4．模拟是用（C 油藏模型）来研究油藏的各种物理性质和流体在其中的流动规律 5．数学模型来研究某个物理过程变化规律是通过求解某一物理过程的（A 数学方程组） 6．数值模型用离散化方法将偏微分方程组转化为（D 有限差分方程组） 7．计算机模型是将各种数学模型的计算方法编制成（C 计算机程序） 8．油藏数值模拟用（B 数值）方法求解油藏数学方程组 9．离散化就是把整体分割为若干（D 单元） 10．有限差分法是对（C 网格范围内的各点）求解。 11．有限差分法使偏微方程被（B 代数方程组）所取代 12．块中心网格是用（A 网格块中心）来表示小块坐标的 13．点是心网格是用（B 结点）的位置来确定小块的中心 14．根据每一组份的质量守恒建立的渗流数学模型称为（A 组份）模型 15．动态预测是在历史拟合的基础上对（C 未来的开发指标）进行计算 16．IMPES 方法是（C 隐式求解压力方程，显式求解饱和度）的 17．半隐式方法是（D 同时求出压力和饱和度）的 18．适定问题是指一个问题的解（D 存在、唯一且稳定） 19．认识油田的主要方法有直接观察法和（B 模拟法） 20．块中心网格和点中心网格的差分方程是(A 一样的) 21．加权六点格式中,要使之成为显式差分格式,必须满足(A θ=0 22．加权六点格式中,要使之成为隐式差分格式,必须满足(C θ=1/2 ) 23．加权六点格式中,要使之成为克兰克.尼克森差分格式,必须满足(B θ=1/2 ) 24．显式差分格式是( A 有条件)收敛的. 25．隐式差分格式是( B 无条件)收敛的 26．显示差分格式的稳定条件是(A △t/△x 2<=1/2) 27．对于各向同性的正方形网格的等值供给半径为(C 0.208△x) 28．差分方程组的直接解法适用于处理(B 系数矩阵阶数不太高)的问题. 64．差分方程组的迭代解法主要用于处理(A 系数矩阵阶数较高)的问题 29．对渗透率取值一般取(A 上游权)的处理方法.) 30．将含水饱和度归一化的公式为(C Sor Sw c Sw c Sw S W ---=1) 31．对数学模型按空间维数分有无零维模型(A 有) 32．历史拟合在含水拟合时,主要是对(C 相对渗透率曲线)的修改. 二、解释概念 1.油藏模拟:是用油藏模型来研究油藏的各种物理性质和流体在其中的流动规律，以便更好地认识油层，作出正确的评价，确定合理的开发方案和提高采收率的措施。 2.数学模型：通过求解某一物理过程的数学方程组来研究这个物理过程变化规律的方法。 3.数值模型：用离散化方法将偏微分方程组转化为有限差分方程组，将其非线性系数线性化，得到线性方程组，然后求解。 4.油藏数值模拟：用数值方法求解油藏数学方程组，就是油藏数值模拟。 5.离散化：离散化就是把整体分割为若干单元来处理。 6.离散空间：离散空间就是把所研究的空间范围套上某种类型的网格，将其划分成一定数量的单元。 7.离散时间：离散时间就是在所研究的时间范围内离散成一定数量的时间段。 8.有限差分法：有限差分法是对网格范围内的各点求解。即原先表示连续的、足够光滑函数的偏微分方程，被一套对每个离散点的、与该点近似解有关的代数方程组所取代。 9.块中心网格：用网格分割成小块的中心来表示小块坐标。 10.点中心网格：由网格的交点，即结点的位置来确定小块的中心。 11.一阶向前差商：对于函数),(t x p ， x p p x p i i ?-=??+1为 一阶向前差商。 12.一阶向后差商：对于函数),(t x p ， x p p x p i i ?-=??-1 为一阶向后差商。 13.一阶中心差商：对于函数),(t x p ， x p p x p i i ?-=??-+211为一阶中心差商。 14.二阶差商：对于函数 ),(t x p ， 2 1 122 2x p p p x p i i i ?+-= ??-+为二阶中心差商。 15.显式差分格式：在这种差分格式中，只有一个未知 数。当已知第n 时刻的n i p 值时，由一个方程就可以求出第n+1时刻的值 1+n i p ，不需联立求解。 16.隐式差分格式：在这种差分格式中，有多个未知数。当已知第n 时刻的n i p 值时，为了求出第n+1时刻的 值 1+n i p ，必须解一个线性代数方程组。 17.网格等值供给半径：当网格中有一口井时，均质地层的网格等值供给半径为：22e 14.0y x r ?+?= 18.黑油模型：黑油模型是简化的组份模型。烃类系统只考虑两个组份：“油”组份是地层油经微分蒸发后在大气压下的残存液（即黑油），而“气”组份是剩余的流体。水相与其它两相不发生质量转移；气可以从油中出入，但油不能汽化为气相。 19.组份模型：油藏内的碳氢化合物是由多种化学成分组成的，在流动过程中由于温度压力的变化，各流动相的各组份之间可能会发生质量转换。根据每一组份的质量守恒建立的渗流数学模型称为组份模型。 20.历史拟合：利用已知的地质、流体性质和特殊岩心分析资料和实测生产历史，输入计算机中，将计算结果和测定的开发指标相比较，修改油层静态参数，直到计算结果和实际动态参数相当接近，达到允许的误差范围为止。 21.动态预测：在历史拟合的基础上对未来的开发指标进行计算。 22.灵敏度试验：将影响开发指标的油层静态资料输入到计算机中并人为地加以改变，观察它们对开发指标的影响，从中找出其影响比较大的参数。 23.IMPES 方法：是指隐式求解压力方程，显式求解饱和度方法。 24.半隐式方法： 联立求解 25.零维模型：描述均质岩石、统、而且系统内的饱和度分布和压力分布是连续的，油藏内任意处的压力发生变化时，整个油藏系统内的压力都随着同时发生变化的一类数学模型。 26.差分方程稳定性：解差分方程组时，如果计算开始时引入的误差在逐层计算过程中的影响逐渐消失或保持有界，则此差方程是稳定的。 27.显式处理：在n+1时刻求解方程组时，若其系数直接用n 时刻的值，为显式处理。 28.半隐式处理：在n+1时刻求解方程组时，若将其系数用泰勒级数展开，并忽略二阶小量，一阶导数用n 时刻的值，则称为半隐式处理。 29.隐式系数处理：在n+1时刻求解方程组时，若将其系数用泰勒级数展开，展开式中的一阶导数用n+1时刻的值，则称为隐式系数处理。 30.定产条件：油藏数值模拟中内边界条件的处理方法之一，即井以一定产量生产。 31.定压条件：油藏数值模拟中假设井以一定流动压力生产，需要把q 用网络压力P 和流动压力wf P 来表示。 32.直接法：直接法就是经有限次数的运算即可求得方程组精确解的方法。 33.迭代法：迭代法是将方程组的求解问题构造为一个无限序列，其极限就是方程组的解。 34.松弛法：松弛法是在赛德尔迭代法的基础上再进一步加快收敛速度，它一般采取在余项上乘一个松弛因子来加快迭代敛速的处理技巧。 35.点松弛法：点松弛法为对网络的每个单元逐次应用松弛法，它是从区域的左下角起由左到右，由下到上的逐点进行计算。 36.最优松弛因子：使逐次松弛法的渐进收敛速度最快的松弛因子通常称为最优松弛因子。 37.不均匀网格：为了模拟油藏的实际情况，划分网格时，在靠近井的附近网格取密一些，而沿径向相外逐渐稀疏，这种网格称为不均匀网格。 38.顺序求解：求解两相流差分方程组中的四个未知数2时，先求得压力项，然后再求饷度的解法。 39.联立求解：求解两相流差分方程组中的四个未知数时，同时求解压力项和饱和度项的解法。 三、填空题 1.认识油田的主要方法有直接观察法和模拟法两种。 2.物理模型可分为定性模型和定量模型两种；前者主要是为了了解油层中发生的各种现象，后者主要是为了得到油田开发过程中的有关定量指标。 3.数值模型就是电解模型和电网模型。 4.物理模拟和数值模拟简称为双模， 5.求微分方程数值解的方法有有限差分法、变分法、有限元法等。 6.油藏数值模拟包括四部分内容：地质模型、 数学模型、数值模型、计算机模型。 7.值模拟输入的数据包括油藏的静态数据和生产井/注入井的动态数据。 8.油藏数学模型的命名，一般依据流体的相态、空间维数、使用功能与特点三个方面。 9.如果一个问题的解存在、唯一且稳定时称该问题为适定问题。 10.建立数学模型常用的物理原理包括：质量守恒原理、能量守恒原理、达西定律。 11.二阶微分方程三种基本类型为：抛物型、椭圆型和双曲型。 12.黑油模型是简化的组份模型。烃类系统只考虑 两个组份。 13.黑油模型中水相与其它两相不发生质量转移；气可以从油中出入，但油不能汽化为气相。 14.定解条件一般包括边界条件和初始条件，前者包括内边界条件和外边界条件。 15.离散化的核心是把整体分为若干单元来处理，它包括空间离散和时间离散。 16.有限差分的网格系统分为块中心网格和点中心网格，二者的离散点数是不同的，但差分方程是一样的。 17.一阶差商包括：一阶向前、一阶向后差商和一阶中心差商。 18.以()t x p ,关于的差商为例，一阶向前差商的截断误差为0 (△x )，一阶向后差商的截断误差为0(△x)， 一阶中心差商的截断误差为0(△x 2 )。 19.以()t x p ,关于x 的差商为例，二阶中心差商的截断误差为0(△x 2 )。 20.微分方程的差分格式有显示差分格式、隐式差分格式、克兰克·尼克森格式和加权六点格式。 21.加权六点格式中，当0=θ时为显示差分格式，当2 1 = θ时为克兰克·尼克森差分格式，当1=θ时为隐式差分格式。 22.显示差分格式是有条件收敛的，隐式差分格式无条件收敛的。 23.显示差分格式是有条件稳定的，隐式差分格式无条件稳定的。 24.显示差分格式的稳定条件是 2/12≤??x t 。 25.对于各向同性的正方形网格的等值供给半径为 0.208△x 。 26.差分方程组的直接解法的特点是计算工作量小，精确度较高、计算程序复杂。 27.差分方程组的迭代解法主要用于处理系数矩阵阶数较高的问题。 28.差分方程组的迭代解法的特点是计算程序较为简单、工作量有时较大。 29.差分方程组的解法分为直接解法和迭代解法两类。 30.传导系数的三种处理方法是算术平均、几何平均和调和平均。 31.渗透率的四种平均取值方法是算术平均、加权平均、几何平均和调和平均。 32.相对渗透率取值一般取上游权的处理方法。 33.差分方程组的系数有显示、半隐式和隐式的处理方法。 34.IMPES 方法是隐式压力显示饱和度的处理方法。 35.IMPES 方法是隐式压力显示饱和度的处理方法。 36.对数学模型按相态分类可分为单相流模型、两相流模型和三相流模型。 37.对数学模型按空间维数分可分为零维模型、一维模型、二维模型、三维模型。 38.对数学模型按模型的使用功能及特点可分为气藏模型、黑油模型、组份模型。 39.常用研究差分方程的稳定性的方法有误差图解法和 Von Neumann 法。 40.由于油藏各点的渗透率不同，束缚水饱和度不同，因而需要对相对渗透率曲线进行归一化处理。 41.将含水饱和度归一化的公式为o r wc wc w w S S S S S ---= 1 42.历史拟合的对象油层平均压力和单井压力、见水时间和含水变化和油气比的变化。 43.历史拟合在压力拟合时，主要调整 孔隙度、饱和度、油层综合压缩系数、渗透率及相对渗透率等。