油藏数值模拟技术课件.ppt

1、油藏数值模拟技术与应用油藏数值模拟技术与应用常晓平河南油田勘探开发研究院 油藏数值模拟技术是一门将油田开发重大油藏数值模拟技术是一门将油田开发重大决策纳入严格科学轨道的关键技术。决策纳入严格科学轨道的关键技术。从油田投产开始，无论是单井动态，还是从油田投产开始，无论是单井动态，还是整个油田动态，都要进行监测与控制。油藏整个油田动态，都要进行监测与控制。油藏数值模拟是油田开发最优决策的有效工具，数值模拟是油田开发最优决策的有效工具，成为油田开发中不可缺少的一项重要工具。成为油田开发中不可缺少的一项重要工具。 模拟 仿真： 油藏数值模拟就是把油藏在油藏数值模拟就是把油藏在三维的空间里分为许多离散的

2、三维的空间里分为许多离散的单元，并且模拟油藏及流体在单元，并且模拟油藏及流体在空间及一系列离散的时间步里空间及一系列离散的时间步里的发展变化。的发展变化。建立油藏模拟软件，一般包括：（1）质量守恒原理（2）能量守恒原理（3）运动方程（达西定律）（4）状态方程（5）辅助方程（如饱和度方程，毛管力方程等） 1) 1)剩余油分布研究。剩余油分布研究。 2)2)优化井网、开发层系、井数和井位。优化井网、开发层系、井数和井位。 3)3)选择注水方式。选择注水方式。 4)4)对油藏和流体性质的敏感性进行研究对油藏和流体性质的敏感性进行研究 5)5)实施方案的可行性评价实施方案的可行性评价数据文件准备初始化

3、计算生产史拟合开发预测结果整理输出每一步的目的？一、所需数据基本类型一、所需数据基本类型：2 2、PVTPVT及岩石性质数据及岩石性质数据 提供流体的地层体积系数、黏度、密度、气油比、岩石及油水提供流体的地层体积系数、黏度、密度、气油比、岩石及油水 的压缩系数、相对渗透率及界面毛管力曲线。的压缩系数、相对渗透率及界面毛管力曲线。3 3、初始化数据、初始化数据提供平衡区，如流体界面、参考深度、压力及毛管力。提供平衡区，如流体界面、参考深度、压力及毛管力。 1 1、静态油藏描述数据、静态油藏描述数据 用离散的网格单元构造油藏的一个几何模型。用离散的网格单元构造油藏的一个几何模型。 对每一个网格单元

4、提供尺寸、海拔高程、孔隙度及渗透率。对每一个网格单元提供尺寸、海拔高程、孔隙度及渗透率。4 4、井数据、井数据 提供井、完井位置、生产提供井、完井位置、生产/ /注入速率、井组及其它数据如表皮注入速率、井组及其它数据如表皮 系数、井半径及井控制等系数、井半径及井控制等5 5、其它：、其它： 聚合物等。聚合物等。1 1、网格数据、网格数据地质模型等地质模型等2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区数据、平衡区数据5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例( (一一) ) 网格的概念网格的概念( (二二) ) 网格设计的发展阶段网格设计的发展阶

5、段( (三三) ) 网格的类型及描述方法网格的类型及描述方法( (四四) ) 网格数据的读取顺序网格数据的读取顺序( (五五) ) 各类网格的优缺点各类网格的优缺点 1、网格的概念：网格的概念： 2、网格描述的目的网格描述的目的: 建立能较准确描述油藏建立能较准确描述油藏特征的地质模型。特征的地质模型。建立地质模型的网格设计方法：建立地质模型的网格设计方法：1 1、选择模型的几何描述、选择模型的几何描述（1 1）研究区域的大小及形态）研究区域的大小及形态（2 2）需要的研究精度）需要的研究精度（3 3）所得数据的详细程度）所得数据的详细程度（4 4）断层结构的复杂性）断层结构的复杂性（5 5）

6、断层两边地层的接触关系及属性）断层两边地层的接触关系及属性（6 6）存在倾斜及下降断层）存在倾斜及下降断层（7 7）建模的时间限制）建模的时间限制3 3、设计区域网格、设计区域网格1 1、基于顶面构造。、基于顶面构造。2 2、只包含一层。、只包含一层。3 3、井尽可能位于网格中心、井尽可能位于网格中心4 4、断层尽可能位于网格边界、断层尽可能位于网格边界5 5、边界为不流动边界、边界为不流动边界6 6、网格单元的大小及形态的变化随着：、网格单元的大小及形态的变化随着：（1 1）井附近的网格一般较小。）井附近的网格一般较小。（2 2）水体位置的网格一般较大）水体位置的网格一般较大在垂向上网格分层

7、基于：在垂向上网格分层基于：1 1、可得的层面数据。、可得的层面数据。2 2、渗透性随深度的变化。、渗透性随深度的变化。3 3、影响垂向上流动的隔挡：、影响垂向上流动的隔挡： 如泥岩夹层的影响程度如泥岩夹层的影响程度 及范围等。及范围等。优点：夹层描述精确； 弱点：网格多，占用空间大；容易与流量相关的导致不收敛问题优点：网格少，省空间；与流量相关的不收敛问题较少。 弱点：垂向上描述较粗，需花较大的功夫调整垂向传导率4 4、设计分层结构、设计分层结构(二二) 网格设计的发展阶段：网格设计的发展阶段：1、块中心、块中心2、角点、角点3、PEBI1、块中心网格、块中心网格2、角点网格、角点网格1 8

8、1 61 41 21 08642642i n j 1p r o d 1p r o d 2( (三三) ) 网格的类型及描述方法网格的类型及描述方法到目前主要有三种网格类型： 块中心、角点及PEBI网格。块中心网格：给出DX、DY、DZ及深度角点网格：需要指定组成每一个网格的四条坐标线坐标（COORD）及八个角点的深度（ZCORN）在块中心网格系统，虽然相邻网格的深度错开，但ECLIPSE还是认为其间是相连的，有流动发生。PEBI网格与角点网格对比直角网格直角网格PEBI网格与角点网格对比Prod1Prod2Inj1( (一一) ) 网格的概念网格的概念( (二二) ) 网格设计的发展阶段网格设

9、计的发展阶段( (三三) ) 网格的类型及描述方法网格的类型及描述方法( (四四) ) 网格数据的读取顺序网格数据的读取顺序( (五五) ) 各类网格的优缺点各类网格的优缺点( (六六) ) 各类网格的优缺点各类网格的优缺点网格类型优点缺点块中心简单、正交性好描述断层等误差较大角点描述地层、特别是断层的能力较强较复杂，正交性不好，需借助专用软件生成PEBI描述各类地层的能力强，正交性好。较复杂，需借助专用软件生成1 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例 来源于来源于

10、PVT分析，数据以表的方式分析，数据以表的方式描述有关参数与压力、温度的关系。描述有关参数与压力、温度的关系。数据项：数据项： 1、各相的地层体积系数（、各相的地层体积系数（Bo、Bg、Bw） 2、各相的黏度、各相的黏度 3、气油比、气油比Rs 和（或）油气和（或）油气比比Rv）区域A为黑油，该区域或许经过泡点线，并且包含溶解气，但是距离临界点较远DD线：表示接近临界点的流体。这时很难确切的知道流体是液体还是气体，因为流体沿着DD经过临界点时，变成了两相混合物并且不易确定过度区。HH线：非等温过程，如在分离器中时。l严格地说，严格地说，黑油黑油在地面条件下不包含在地面条件下不包含溶解气。溶解气

11、。l黑油模拟器最好用在远离临界点的单黑油模拟器最好用在远离临界点的单相区域。相区域。l如果流体的组成变化比较大，就不应如果流体的组成变化比较大，就不应该使用黑油模拟器。该使用黑油模拟器。黑油模型：黑油模型： 黑油模型是描述含有非挥发组分的黑油和挥发性组分的黑油模型是描述含有非挥发组分的黑油和挥发性组分的原油溶解气在油藏中运动规律的数值模型。黑油模型也称低挥原油溶解气在油藏中运动规律的数值模型。黑油模型也称低挥发油双组分模型。发油双组分模型。组分模型：组分模型： 组分模型是用来研究高挥发性烃类的系统在油气藏中运动组分模型是用来研究高挥发性烃类的系统在油气藏中运动规律的数值模型，也成多组分模型。规

12、律的数值模型，也成多组分模型。1、因为凝析及压力降低而释放出来的气体、因为凝析及压力降低而释放出来的气体的总量只占总含烃量的很少一部分。的总量只占总含烃量的很少一部分。2、由于凝析及压力降低而释放出来气体后，、由于凝析及压力降低而释放出来气体后，流体的烃组成变化不明显。流体的烃组成变化不明显。3、路径远离临界点。、路径远离临界点。4、整个过程是等温的。、整个过程是等温的。如果上述条件不能完全满足，则需要使如果上述条件不能完全满足，则需要使用全组分模拟器。用全组分模拟器。为了保证使用黑油模拟器可以达到足够的精度，为了保证使用黑油模拟器可以达到足够的精度，油藏流体必须具备下列条件：油藏流体必须具备

13、下列条件：1 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例l来源于来源于SCAL（special core analysis）分析，）分析，数据主要是以表的方式描述属性与饱和度的数据主要是以表的方式描述属性与饱和度的关系。关系。l对于油藏的不同部分，可分用不同的相渗曲对于油藏的不同部分，可分用不同的相渗曲线。线。l在每个区域，需要设置最大、最小及临界饱在每个区域，需要设置最大、最小及临界饱和度值。和度值。l用于定义过渡带的饱和度。用于定义过渡带的饱和度。 岩心分析岩心分析

14、岩石数据是特定的岩心分析试验的结果，岩石数据是特定的岩心分析试验的结果，该数据用于：该数据用于：l设置每一流体相的最大、最小饱和度，该设置每一流体相的最大、最小饱和度，该值用于定义平衡区的相饱和度。值用于定义平衡区的相饱和度。l定义过渡带的范围及属性。定义过渡带的范围及属性。l描述各相在网格块间流动时的动态表现。描述各相在网格块间流动时的动态表现。1 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例 在原始状况下，认为在同一个在原始状况下，认为在同一个平衡区中的流体的静力学是平

15、衡平衡区中的流体的静力学是平衡的。的。 用于定义油藏初始条件下的用于定义油藏初始条件下的每一个网格单元的压力、饱和度每一个网格单元的压力、饱和度值。值。初始压力计算：l也可以在PROPS 部分中使用SWATINIT 定义每个单元的原始含水饱和度值。SW变化相对较为剧烈在。1 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法l井流量方程：l单井的完井传导率系数：l井的压力等值系数：如果有效的井眼半径超过了re的一半，将会出现警告；如果有效的井眼半径超过了re，Twj变成了负值从而导致了错误的结果l井的生产指数：初

16、始化阶段应完成的主要工作初始化阶段应完成的主要工作l数据查错： 1、文本方式 - *.PRT 2、OFFICE 环境中。l储量拟合： 1、油水界面深度检查。2、油藏每一层边界检查。3、地质模型检查： Q=A*H*POR*SO1、 软件依次读入及处理各部分数据。在开始读入下部软件依次读入及处理各部分数据。在开始读入下部分内容前，必须对本部分各种数据进行一致性检查。分内容前，必须对本部分各种数据进行一致性检查。只有最后一部分例外，因为该部分数据是时间相关只有最后一部分例外，因为该部分数据是时间相关的，因此，总体上初始化时并不读入与处理。的，因此，总体上初始化时并不读入与处理。2、软件开始运行的第一

17、个任务是为输入的数据分配内存。、软件开始运行的第一个任务是为输入的数据分配内存。3、处理模拟网格的几何情况及属性，得到更利于流动计算、处理模拟网格的几何情况及属性，得到更利于流动计算的三维空隙体积、传导率、单元中部深度及可产生流动的三维空隙体积、传导率、单元中部深度及可产生流动的网格间的连接。的网格间的连接。4、 在初始化（在初始化（initialization）时，计算油藏每一层的静压力）时，计算油藏每一层的静压力剃度，并给每一个网格赋每一相的饱和度值剃度，并给每一个网格赋每一相的饱和度值l给定平衡区数据： 平衡区数据被用来初始化平衡区数据被用来初始化饱和度分布、基准深度的压饱和度分布、基准

18、深度的压力及压力梯度。力及压力梯度。l用枚举的方式给出每一个单元的初始饱和度及压力。 修改现有的模拟数据（如渗透率、相渗曲线、修改现有的模拟数据（如渗透率、相渗曲线、井的生产指数等）直到能与观察数据（如压力、井的生产指数等）直到能与观察数据（如压力、产量、含水等）进行合理的对比，这种过程叫做产量、含水等）进行合理的对比，这种过程叫做生产史拟合。生产史拟合。 生产史拟合过程也就是一个反复修改参数、反复试算生产史拟合过程也就是一个反复修改参数、反复试算的过程，需要耗费大量的机时和人力。的过程，需要耗费大量的机时和人力。拟合结果的好坏，不拟合结果的好坏，不但取决于人们对油藏的认识，还取决于油藏工程师

19、的拟合经但取决于人们对油藏的认识，还取决于油藏工程师的拟合经验和处理技巧。验和处理技巧。 由于油藏地质情况的复杂性及反求参数的多解性，不同由于油藏地质情况的复杂性及反求参数的多解性，不同的人拟合同一油田，其试算次数和拟合结果可能均不同。的人拟合同一油田，其试算次数和拟合结果可能均不同。 为了减少拟合工作的盲目性，有必要了解各项油藏参数为了减少拟合工作的盲目性，有必要了解各项油藏参数与拟合对象的相互关系，同时还应一些拟合经验和处理技巧与拟合对象的相互关系，同时还应一些拟合经验和处理技巧。 拟合的对象有：拟合的对象有：压力、流量、油压力、流量、油气比、油水比等气比、油水比等。l在油藏的历史拟合中，

20、一般把产量作为已知条件来拟合其它动态参数。拟合步骤为：开始油藏原始平衡状态检查（零流量模拟）拟合井底压力拟合油田平均压力拟合单井压力拟合油田综合含水率拟合单井含水率油藏压力已经拟合好了？结束生产指数控制着单井流量开始下降的时间，对预测单井的动态意义较大。 可以单个的或集体的改变几个参数。通常按下述范可以单个的或集体的改变几个参数。通常按下述范围进行修改：围进行修改： 1、岩石数据的修改：、岩石数据的修改： a.渗透率渗透率 b.孔隙度孔隙度 c.厚度厚度 d.饱和度饱和度 2、流体数据的修改：、流体数据的修改： a.压缩性压缩性 b.PVT数据数据 3、相对渗透率数据：、相对渗透率数据： a.

21、相对渗透率曲线的移位相对渗透率曲线的移位 b.临界饱和度数据的移位临界饱和度数据的移位 4、单井完井数据：、单井完井数据： a.表皮效应表皮效应 b. 井底流动压力。井底流动压力。 在历史拟合过程中，由于模型参数数量多，在历史拟合过程中，由于模型参数数量多，可调的自由度很大。为了避免或减少修改参数的可调的自由度很大。为了避免或减少修改参数的随意性，在历史拟合开始时，必须确定模型参数随意性，在历史拟合开始时，必须确定模型参数的可调范围，使模型参数的修改在合理的、可接的可调范围，使模型参数的修改在合理的、可接受的范围内，这是历史拟合的原则。受的范围内，这是历史拟合的原则。 确定参数的可调范围是一项

22、重要而细致的确定参数的可调范围是一项重要而细致的工作，需要油藏工程师和地质师共同努力，需要工作，需要油藏工程师和地质师共同努力，需要收集和分析一切可以利用的资料。首先分清哪些收集和分析一切可以利用的资料。首先分清哪些参数是确定的，即准确可靠的，哪些参数是不确参数是确定的，即准确可靠的，哪些参数是不确定的，即不准确、不可靠的，然后根据情况，确定的，即不准确、不可靠的，然后根据情况，确定可调范围。通常只修改不确定的参数，而确定定可调范围。通常只修改不确定的参数，而确定的参数一般不允许修改。的参数一般不允许修改。 一般情况下，认为：一般情况下，认为： （1 1）孔隙度：）孔隙度： 如果油层大量岩心分

23、析资料表明，油层部分孔隙度如果油层大量岩心分析资料表明，油层部分孔隙度在在1919到到2121之间，平均为之间，平均为2020，变化范围不大。则，变化范围不大。则把孔隙度视为确定参数，不做修改，或允许改动范围在把孔隙度视为确定参数，不做修改，或允许改动范围在3 3。 （2 2）渗透率：渗透率： 渗透度在任何油田都是不定参数。这不仅是由于测渗透度在任何油田都是不定参数。这不仅是由于测井解释的渗透率值和岩心分析值误差较大，而且根据渗井解释的渗透率值和岩心分析值误差较大，而且根据渗透率的特点，井间的渗透率分布也是不确定的。因此对透率的特点，井间的渗透率分布也是不确定的。因此对渗透率的修改，允许范围较

24、大，可放大或缩小渗透率的修改，允许范围较大，可放大或缩小2 23 3倍或倍或更多。更多。 （3）岩石与液体的压缩系数：）岩石与液体的压缩系数： 液体的压缩系数是实验测定的，变化范围很小，认为液体的压缩系数是实验测定的，变化范围很小，认为是确定的。而岩石的压缩系数虽然也是实验室测定的，是确定的。而岩石的压缩系数虽然也是实验室测定的，但受岩石内饱和液体和应力状态的影响，有一定变化范但受岩石内饱和液体和应力状态的影响，有一定变化范围，考虑这部分影响，允许岩石的压缩系数可以扩大一围，考虑这部分影响，允许岩石的压缩系数可以扩大一倍。倍。 （4）初始流体饱和度和初始压力：）初始流体饱和度和初始压力： 和通

25、常做法一样，认为是确定参数。必要时允许小范和通常做法一样，认为是确定参数。必要时允许小范围内修改。围内修改。 （5）相对渗透率曲线：）相对渗透率曲线： 由于油藏模拟模型的网格粗，网格内部存在严重非均由于油藏模拟模型的网格粗，网格内部存在严重非均质，其影响不可忽视，这与均质岩心的情况不同。因此质，其影响不可忽视，这与均质岩心的情况不同。因此相对渗透率曲线应看作是不定参数。在拟曲线的研究中，相对渗透率曲线应看作是不定参数。在拟曲线的研究中，给出了较好的初始值，但仍允许做适当修改。给出了较好的初始值，但仍允许做适当修改。（6 6）油气的）油气的PVTPVT性质：性质： 视为确定参数。视为确定参数。

26、（7 7）油水界面：）油水界面： 在资料不多的情况下，允许在一定范围内修改。在资料不多的情况下，允许在一定范围内修改。 当油藏未投产时，油藏当油藏未投产时，油藏压力系统处于静止平衡状态，压力系统处于静止平衡状态，相邻点不应有压力势差；否相邻点不应有压力势差；否则会出现流动，打破静态平则会出现流动，打破静态平衡。衡。 当油藏局部区域出现高压异常与低压异常时：（1）沿高压异常区向低压异常区方向增大渗透率值，可使两类异常区同时得到改善。 但修改K值应在渗透率图上成片修改，以保持岩性的连续性。 （2）若修改K值效果不佳，则可减少高压异常区储量，增加低压异常区储量。 当油藏采出大量液体后压力无当油藏采出

27、大量液体后压力无明显下降，可考虑是否因岩石弹性明显下降，可考虑是否因岩石弹性能量过大所致。可减少能量过大所致。可减少岩石压缩系岩石压缩系数数；当油藏压力下降过快时，也可；当油藏压力下降过快时，也可增加岩石压缩系数加以改善。但此增加岩石压缩系数加以改善。但此时还应考率到是否是油藏储量偏小时还应考率到是否是油藏储量偏小所致。若是，则应增加所致。若是，则应增加油藏储量油藏储量。 单元含水率主要受相对渗透率曲线的影响。含水率主要受相对渗透率曲线的影响。油井的见水时间和水驱前缘的拟合比较难以实现。它油井的见水时间和水驱前缘的拟合比较难以实现。它还受到模型的时间步长和网格密度的限制还受到模型的时间步长和网

28、格密度的限制。当见水区域或油井含水上升过快时，可以通过调整油当见水区域或油井含水上升过快时，可以通过调整油相、水相相对渗透率曲线来控制。增大水相临界饱和度值相、水相相对渗透率曲线来控制。增大水相临界饱和度值可以减缓见水，而增大可以减缓见水，而增大Kro/KrwKro/Krw的值可降低含水上升速度。的值可降低含水上升速度。 在拟合上油藏综合含水率之后，还有必要对某在拟合上油藏综合含水率之后，还有必要对某些单井含水率进行拟合。些单井含水率进行拟合。 单井含水率的拟合主要通过调整单井含水率的拟合主要通过调整平面上平面上的渗透的渗透率场及率场及纵向上纵向上的生产指数两个方面来做。的生产指数两个方面来做

29、。 主要是通过修改生产井的采油指数（包括表皮系数）主要是通过修改生产井的采油指数（包括表皮系数）来实现。流动压力偏小，可增大采油指数；反之，可减小来实现。流动压力偏小，可增大采油指数；反之，可减小采油指数采油指数。 岩石数据通常是二维或三维网格矩阵输入值。修改网格数据从不按单个单元来进行，而是在需作必要变化的面积上进行的。 多孔岩石内膨胀能量的数量的大小直接与三种流体相的多孔岩石内膨胀能量的数量的大小直接与三种流体相的饱和度有关，油层能量的饱和度有关，油层能量的影响影响表现在几方面：表现在几方面： 1.油层压力水平油层压力水平 2.产量产量 3. 油气比油气比 当计算中无法保持一定地区给定的产

30、量、进而油层压力当计算中无法保持一定地区给定的产量、进而油层压力降低较严重时，通常这是由于油层中储油量不充足和流入给降低较严重时，通常这是由于油层中储油量不充足和流入给定位置的流量不够所引起的。定位置的流量不够所引起的。 这时除了改变这时除了改变孔隙体积孔隙体积可以校正这种情况以外，也可通可以校正这种情况以外，也可通过过提高油饱和度提高油饱和度来进行校正。来进行校正。 在模拟研究中，流体数据一般是一知的。这些数据包括：在模拟研究中，流体数据一般是一知的。这些数据包括： 1. 地层体积系数地层体积系数 2. 粘度粘度 3. 压缩性压缩性 4. 溶解气数据溶解气数据 双河双河IV下数模研究下数模研

31、究 动态预测可用来研究预测包括枯竭开采、动态预测可用来研究预测包括枯竭开采、注水开发、钻加密井等各种方案的效果。注水开发、钻加密井等各种方案的效果。1.设置单井的生产目标2.设置井组的生产目标井底流压井口压力产油率产水率产气率产液率井组产率 压力控制一般用于预测阶段，一般不适用于生产史拟合阶段。 在压力控制下，模拟器计算相应的流量；如果计算的流量超过了指定的速率，则会切换回速率控制模式下。示例：一注一采两口井的剖面模型：BHP产油量产水量产油/水量压力控制方式控制方式定生产压差求产： 经济限制是通过修改生产策略来满足所定的最低经济目标（如产油率、油气比等）的一种方法。当抵达经济限制条件时，根据

32、事先的定义，一系列的措施就会自动实施。以经济极限制为约束求产：采用经济限制可设定：经济限制可应用于：油田/井组 产量GECON井的产量WECON层的产量（ratios only）CECON下面的措施在计算满足经济限制条件时，可自动实施：关井。关井。关层。关层。打水泥塞。打水泥塞。采用经济限制用法说明1用于在需要时自动进行作业及关井。当抵达含水率、GOR、OGR以及流量限制条件时，自动进行作业一旦井、层关闭后，缺省情况下，要一直保持关闭状态可以用WTEST通过周期性测试将井重新打开。采用经济限制用法说明2如果指定了第二含水率限制条件，井可以被再次打开。如果指定了第二含水率限制条件，井可以被再次打

33、开。如果井达到经济限制条件，模拟可以被合法地终止。如果井达到经济限制条件，模拟可以被合法地终止。在井名中可以使用通配符。在井名中可以使用通配符。井组的经济限制设置方法作业方法与作业方法与WECON一样一样层段的经济条件限制方法：自动测试开井：自动测试开井： 通常，井通常，井/ /层关闭后就不再被打开。自动层关闭后就不再被打开。自动测试用来重新打开原来因无力或经济原因关闭测试用来重新打开原来因无力或经济原因关闭的井的井自动测试开井：1、井名2、测试间隔时间3、测试当时关闭井的原因及目前的有效性检查 P：物理原因（如BHP 、THP） E： 经济原因 G： 井组或油田目标限制 D： 超过最大的TH

34、P C： 层的关闭原因4、测试的最大时间步5、重新开井时间，天1.设置单井的生产目标2.设置井组的生产目标井组控制：用WELSPECS将井放在井组中生产井的产量用GCONPROD关键字注入井的流量用GCONINJE关键字生产井与注入经应该分别放在不同的井组中FIELD 永远存在在井组中的井应分别在WCONPROD及WCONINJE中指定控制模式为GRUP。井组控制：井组井组/ /油田油田 生产控制，井的生产是以相对于指导速率生产控制，井的生产是以相对于指导速率（guide rateguide rate）的比率，该比率服从于他自己的速率及压）的比率，该比率服从于他自己的速率及压力限制，指导速率缺

35、省为该井的生产能力。力限制，指导速率缺省为该井的生产能力。井组控制：亏空补偿（亏空补偿（voidage replacemnetvoidage replacemnet）: :如确定注采比如确定注采比（voidagevoidage replacement ratio replacement ratio）优先级控制（优先级控制（prioritypriority）: :定义一个打开井次序的优先级定义一个打开井次序的优先级序列，在完成井组或油田目标的前提下，根据需要、按顺序列，在完成井组或油田目标的前提下，根据需要、按顺序打开不同的井。序打开不同的井。井组控制：本井组的产量在父井组中占的百分比井组中的生

36、产控制：该井组产量是由子井组及井完成的，其中各井组及井的产量该井组产量是由子井组及井完成的，其中各井组及井的产量服从各自定义的生产量百分数（服从各自定义的生产量百分数（GUIDE RATESGUIDE RATES）。）。在井组目标控制下不在井组目标控制下，以单井控制为准指定井及井组是否处于井组目标控制下：定气油比生产： 通过采取定气油比的方式，可以控制日产油量。定气油比生产：备注：如果base GOR 为零或负数，则自动关闭该井的气油比控制。WCONPROD O1 OPEN LRAT 3* 60 5* / O2 OPEN LRAT 3* 75 5* / O3 OPEN LRAT 3* 75 5

37、* / O4 OPEN LRAT 3* 90 5* / O5 OPEN LRAT 3* 60 5* / O6 OPEN LRAT 3* 75 5* / / GCONINJE G1 WATER VREP 3* 1.2 / 定注采比为1.2/ WCONINJE -I1 WATER OPEN GRUP 150 / -I2 WATER OPEN GRUP 150 / I1 WATER OPEN GRUP / I2 WATER OPEN GRUP / / 定注采比生产：井组控制示例：井组控制示例：井组注入控制：亏空补偿系数：井组注入量=补偿系数*采液量-注入的其它相。亏空补偿系数被通常用来设定为注采比值

38、。井打开顺序优先级设置：PRIORITY设置确定每口井优先级的计算方法，格式为：如果A=D=E=F=H=0，B=C=G=1 则：P=(Q0+Qw)/Qw ,就是说含水率的倒数越高，优先级越高，也即是说，含水率越小的井优先级越高。示例井打开顺序优先级设置：当不用PRIORITY计算优先级时，可用WELPRI逐井定义各自的优先级。WELPRI-1 2井名 优先级号设置在优先级控制下（GCONPRI）的井组目标： 油田产油率为10000，以PLAT开头的井组的产液率限制为6000，井按优先级顺序被打开直至达到组产量限制值。井组产气的限制值为20000，如果达到限制值，则对最差的井进行作业措施。示例

39、在动态预测阶段，根据所处的开发阶段、计划采取的开发方式以及研究目的的不同所采取的动态预测方式也有所不同。主要有：定流量求产及定压求产两种预测方式。动态预测工作流程：1. 在sch.inc模块中使用WELSPECS、WCOMPDAT以及WCONPROD/WCONINJE 分别定义井位、生产层位及流量（如果需要定义井组的流量，则需使用GCONPROD及GCONINJE来定义）。2. 检查sum.inc文件中是否定义了足够的输出项3. 数模计算4. 数据整理 在预测阶段，首先要确定预测控制方式，在预测阶段，首先要确定预测控制方式，通常情况下，都是采用通常情况下，都是采用定流量定流量求产的方式求产的方

40、式-即：设定所有井的注采量预测开发效果。即：设定所有井的注采量预测开发效果。单井定义示例： 井组定义示例： 注采比1.01、ECLIPSE2、VIP3、WorkBench4、CMG5、SURE现象现象：局部压力异常。：局部压力异常。目标目标：降低高压区，提高低压区。：降低高压区，提高低压区。方法方法：要实现上述目标，可通过下述一个或几个方法来完成：要实现上述目标，可通过下述一个或几个方法来完成： 1 1、改变岩石渗透率，使流体从高压带流向低压带。、改变岩石渗透率，使流体从高压带流向低压带。 2 2、减少高压带得储油量：、减少高压带得储油量： 3 3、通过下述的任意方法增加低压带的储油量：、通过

41、下述的任意方法增加低压带的储油量： 方法：方法：最可行的方法是使流体从高压区穿过油层流到低压区。渗透率修改前后对比图压力场修改前后对比图：通过给孔隙度减小一定的比值，比如0.96：孔隙度新=孔隙度原始*比值,由图可见，整体压力水平降低了。整体压力水平降低了。现象现象：。现象现象：甚至在大量采出流体以后，模型中压力降不明：甚至在大量采出流体以后，模型中压力降不明显。显。原因原因：岩石压缩性较大，因而使饱和度变化的影响可：岩石压缩性较大，因而使饱和度变化的影响可以忽略不记，从有效压缩性方程就可看出：以忽略不记，从有效压缩性方程就可看出： Ce=(SoCo+SwCw+SgCg+Cj)/So处理处理：

42、调整岩石的压缩系数：调整岩石的压缩系数现象现象：没有任何注入的和流入的水，模型内含水饱和：没有任何注入的和流入的水，模型内含水饱和度出现增加。度出现增加。原因原因：岩石压缩性太低，引起孔隙空间自由体积的扩：岩石压缩性太低，引起孔隙空间自由体积的扩大。在模型中，这种自由体积充填着非流动相，通常大。在模型中，这种自由体积充填着非流动相，通常为水。因此，在水的物质平衡中出现了额外的数量。为水。因此，在水的物质平衡中出现了额外的数量。处理处理：校正岩石的压缩系数：校正岩石的压缩系数现象现象：一个油田一直在生产。并：一个油田一直在生产。并且总油气比与累计产量的关系曲且总油气比与累计产量的关系曲线如图线如

43、图原因原因：油气比太高，说明：油气比太高，说明Kg/Ko曲线太乐观。曲线太乐观。处理处理：要降低这过高的采气量，：要降低这过高的采气量，就要将就要将Kg/Ko曲线从左向右移位，曲线从左向右移位，以减小给定流体饱和度下气体流以减小给定流体饱和度下气体流动数量，动数量，现象现象：情况见图。模型自由气过：情况见图。模型自由气过早开始流动；但曲线具有通常观早开始流动；但曲线具有通常观察到的曲线形状。察到的曲线形状。原因原因：临界气饱和度偏低。：临界气饱和度偏低。处理处理：临界气饱和度移位：临界气饱和度移位,推迟气推迟气的生产，直到气饱和度达到较高的生产，直到气饱和度达到较高值为止。值为止。现象现象：情

44、况见图。：情况见图。原因原因：这说明由两个或两个以上：这说明由两个或两个以上参数对计算结果有较大影响，必参数对计算结果有较大影响，必须至少修改两个参数才能满足要须至少修改两个参数才能满足要求。求。处理处理：这种情况的拟合难度较大。：这种情况的拟合难度较大。具体需改那些参数，要视实际情具体需改那些参数，要视实际情况而定。况而定。现象现象：在某一油层中，即使在可靠范围：在某一油层中，即使在可靠范围内改变油层数据后，下述异常情况仍内改变油层数据后，下述异常情况仍未得到改变：未得到改变：1.矿场观察到的压力仍然太高而不能匹矿场观察到的压力仍然太高而不能匹配。配。2.矿场观察到的产量，与模型数字相比，矿

45、场观察到的产量，与模型数字相比，仍然太高。仍然太高。原因原因：油层能量未得到认识。油层从目前所确定的范围外得到某种形式的能：油层能量未得到认识。油层从目前所确定的范围外得到某种形式的能量。这种油层能量的来源可以是：量。这种油层能量的来源可以是：1.1. 其它为开发的产层，或其它为开发的产层，或2.2. 影响油层的未认识的水区。影响油层的未认识的水区。处理处理：1. 1. 在存在砂层的地区扩大其模型面积，再运算；在存在砂层的地区扩大其模型面积，再运算； 2. 2. 在给定区域模拟影响油层的含水域。在给定区域模拟影响油层的含水域。注意注意：连通地段的存在，应作为最后的拟合手段来加以考察，如果证据不

46、是：连通地段的存在，应作为最后的拟合手段来加以考察，如果证据不是很充分的话，工程师不应该明显地增加生产面积。在这些情况下，建议很充分的话，工程师不应该明显地增加生产面积。在这些情况下，建议补充钻井已探明该面积，通常是个好办法。补充钻井已探明该面积，通常是个好办法。油藏物性受沉积相带的影响比较大。不同的相带，由于沉积环境不同，其岩石颗粒大小、分选等均不同，因此油藏流体流动时所适用的相对渗透率曲线也不同 。1、从建模软件中输出沉积相建模三维数据体。2、在.DATA文件中，根据沉积相带类型的数目定义关键字TABDIMS中所用的相渗曲线个数。3、在scal文件中给出每一个沉积相带类型所使用的相渗曲线。

47、4、在reg文件中，定义每个网格所使用的相渗曲线。1、精度提高。2、工作效率更高调整更加便利，可以分区域调整相应的相渗曲线。1、地质建模前期需作深入的沉积相研究、研究周期较长。2、进行相控地质建模的技术难度较大。 油藏数值模拟是一门技术密集、工作难度较大的应用学科，它涉及地质、油藏、采油、计算机等多学科知识，只有经过多次的应用实践研究才能真正地掌握该项技术。主运行文件: *.DATA*.DATAEXAMPLE_GGO.INC地质框架描述EXAMPLE_GHDR.INC 头文件,不用改。EXAMPLE_GOPP.INC网格参数打印输出定义EXAMPLE_GOTH.INC杂项定义，不用改。EXAMPLE_GPRO .INC 属性描述（dznet/poro/perm）EXAMPLE_REG .INC区域定义，如储量区、PVT区等EXAMPLE_SUM .INC动态数据输出定义EXAMPLE_INIT .INC初始化定义，如平衡区等EXAMPLE_SCAL .INC岩心试验数据定义，如相渗曲线等EXAMPLE_PVT .INCPVT试验EXAMPLE_SCH .INC动态数据定义