油藏数值模拟基础培训

油藏数值模拟基础培训油藏数值模拟基础培训求实创新超越REALITY，INNOVATION,TRANSCENDENCY目录第一讲油藏数值模拟总论第二讲概念模型数据文件解剖第四讲RUNSPEC/GRID要点分析油藏数值模拟基础培训求实创新超越REALITY，INNOVATION,TRANSCENDENCY油藏数值模拟总论油藏数值模拟总论求实创新超越REALITY，INNOVATION,TRANSCENDENCY油藏数值模拟的概念为什么要做油藏数值模拟油藏数值模拟的方法原理油藏数值模拟的作用油藏数值模拟的大体流程油藏数值模拟都需要什么参数油藏数值模拟的概念求实创新超越REALITY，INNOVATION,TRANSCENDENCY标准定义：应用已有规律，采用数值方法求解描述油藏内流体流动问题，并利用计算机研究油藏开发及动态规律的一门技术。■油藏数值模拟工程定义：从地下流体渗流过程中的特征出发，建立描述渗流过程的基本物理现象，并能描述油藏边界条件和原始状况的数学模型，借助计算机计算求解描述油气藏渗流数学模型，并结合油藏地质学、油藏工程学等学科知识重现油田开发的全过程，主要用于解决油田开发实际问题。简单地说就是在电子计算机上开发油藏。油藏数值模拟是一种运用较复杂的数学方法预测油藏动态的一门特殊的油藏工程研究方法，本身赋予较重比例的数学元素，但应用领域属于油藏工程。油藏数值模拟的概念求实创新超越REALITY，INNOVATION,TRANSCENDENCY■油藏工程定义一：是一门以油田地质学和渗流力学为基础，以油藏数值模拟为手段，研究油气田开发设计和工程分析方法的综合性学科。定义二：是一门高度综合的技术学科。它综合应用地球物理、油藏地质学、油层物理、渗流理论和采油工程等方面的成果以及所提供的信息资料，对油藏开发方案进行设计和评价，以及应用这种预测结果提出相应的技术措施，以获得最大的经济采收率。油藏工程含义宽泛，数值模拟是油藏工程研究技术方法之一，理论基础与其它的油藏工程动态预测数学计算方法一致。随着油藏数值模拟技术正规化与综合性的发展进步，其本身的内涵与功能逐步扩大，与传统油藏工程的含义逐步融合交叉，成为“现代油藏工程”的核心，但仍不能取代传统的油藏工程方法。为什么要做油藏数值模拟求实创新超越REALITY，INNOVATION,TRANSCENDENCY■油藏工程与数值模拟油藏工程应用观测：信息采集与分析假设：油藏特征判断计算：定量计算与评估决策：提出措施或方案类比法：资料相对缺乏时的有效研究方法，但无法判断不同开发策略。水动力学法：机理明确，但过于简化，只适合于早期开发指标概算。物质平恒法：忽略非均质性和动态参的分布，适合于宏观开发指标计算，不适合于油田局部动态预测。数值模拟法：能考虑油藏的复杂几何形状、非均质性、岩石和流体性质变化、井网方式和产量等因素，是迄今为止油藏动态研究中考虑因素最多的一种方法。采收率统计关系式、水驱特征曲线、产量递减方程及经验公式等。垂向非均质管流模型、等值渗流阻力法、Buckey-Leverett方程非活塞驱替模型等。为什么要做油藏数值模拟求实创新超越REALITY，INNOVATION,TRANSCENDENCY■油藏工程与数值模拟油藏数值模拟技术优势：可以实现预测风险最小化。油藏非均质与岩石各相异性流体性质及岩石流体渗流特征空间变化烃类油藏开采机理的复杂性各种油藏工程方法的适应性影响预测风险的因素输入油藏地质模型数据选择合适的模拟器可靠的工程实践正确的模拟应用地质建模软件功能应用技能数值模拟技术实现为什么要做油藏数值模拟求实创新超越REALITY，INNOVATION,TRANSCENDENCY■油藏工程与数值模拟传统油藏工程方法对油藏数值模拟的作用在于：弥补收集资料数据的不充足减小输入数据的误差深化对模型合理性的认识辅助动态历史拟合分析利于模拟结果的评价与解释

以上这些都是影响数值模拟应用研究质量的关键，也是体现从业者能力大小的重要方面。油藏数值模拟的方法原理求实创新超越REALITY，INNOVATION,TRANSCENDENCY开采过程非线性偏微分方程单/多相流公式非线性代数方程离散化线性代数方程线性化①建立数学模型②建立数值模型A、通过离散化将偏微分方程组转换为有限差分方程组。B、将非线性系数线形化，得到线形代数方程组。Ａ、通过质量/能量守恒方程、状态方程、运动方程、辅助方程建立基本方程组。B、根据所研究的具体问题建立相应的初始和边界条件。求实创新超越REALITY，INNOVATION,TRANSCENDENCY油藏数值模拟的方法原理压力、饱和度分布及井流量求解油藏模拟器程序化模拟应用商业化②建立数值模型C、应用线性代数方程组的解法求得未知量（压力、饱和度等）的数值分布。③建立计算模型A、将求解数值模型进行程序化。B、连同前后处理软件，形成油藏模拟器。④油藏模拟应用利用油藏模拟器研究和解决具体的油气田开发问题。求实创新超越REALITY，INNOVATION,TRANSCENDENCY油藏数值模拟的方法原理1.建立数学模型

2.建立数值模型

通过离散化，将连续的偏微分方程组转换成离散的有限差分方程组，再用多种方法将非线性系数线性化，成为线性代数方程组，然后求解线性代数方程组。

即建立一套描述油藏中流体渗流的偏微分方程组，包括初、边值问题。偏微分方程组线性代数方程组得到压力、饱和度等有限差分方程组离散化线性化解方程组求实创新超越REALITY，INNOVATION,TRANSCENDENCY油藏数值模拟的方法原理3.建立计算机模型

将资料（静、动态）的输入，系数矩阵和常数项的形成，多种解法和结果的输出等，编制成计算机程序。

数值模拟的关键是计算的精度和速度。计算的精度取决于离散的程度、数值计算误差、离散方程的稳定性，计算的速度取决于计算机速度、解法速度和模型规模。从离散的程度看，速度和精度相矛盾，要根据解决问题的需要选择离散化程度和计算速度。油藏数值模拟的方法原理求实创新超越REALITY，INNOVATION,TRANSCENDENCY质量守恒方程/连续性方程：根据物质平衡原理，流入单元体中的流体流量减去流出单元体流体流量等于单元体流体质量变化。运动方程：达西定律。流动方程：状态方程：流体和岩石体积与压力的关系。辅助方程：饱和度归一，毛管压力、密度、粘度及相渗曲线的相关函数。

求实创新超越REALITY，INNOVATION,TRANSCENDENCY——几点说明●描述油藏流体渗流这一具体物理过程的完整的数学模型是非线性的偏微分方程，不宜直接求解，需要通过离散转化成比较容易求解的代数方程组。离散方法一般为有限差分法。●离散后形成的代数方程组是非线性的差分方程组，还要采用某种线性化方法将其线性化，然后求解。常用的线性化方法有显示方法、半隐式方法或全隐式方法等。●对数学模型离散化并线性化，得到每个网格节点上的一个（单相）或多个（多相）线性代数方程。每个方程除含有本点上的未知变量外，一般还含有相邻节点上的未知变量。因此，为了求得线性方程的解，需要将各点上的方程联立，形成联立代数方程组。该方程组一般为大型稀疏方程组。●求解线性代数方程组所用的方法有直接法和迭代法两大类，直接法常用的有高斯消去法、主元素消去法、D4方法等；迭代法常用的有交替方向隐式方法、超松弛达代法、强隐式方法等。油藏数值模拟的方法原理ECLIPSE软件介绍求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE100:黑油模拟器

ECLIPSE300:组分模拟器

Flogrid:建模型

PVTi:EOS分析

VFPi:垂直管流

Schedule:动态数据

SCAL:岩芯数据

OFFICE:项目管理

Frontsim:流线法

Floviz:3维显示

Simopt:历史拟合

R2SL:地面管网与模拟模型偶合■ECLIPSE主要模块油藏数值模拟的作用求实创新超越REALITY，INNOVATION,TRANSCENDENCY在理论上：探索多孔介质中各种复杂渗流问题的规律；在工程上：作为开发方案设计、动态监测、开发调整、反求参数、提高采收率的有效手段，能为油气田开发中的各种技术措施的制定提供理论依据。——有效的油田开发科学决策工具！渗流机理研究开发可行性评价参数敏感性分析开发方案优化剩余油分布研究提高采收率研究动态跟踪研究反演油藏地质模型油藏数值模拟的作用求实创新超越REALITY，INNOVATION,TRANSCENDENCY1、模拟初期开发方案1）实施方案的可行性评价；2）选择井网、开发层系、井数和井位；3）选择注水方式；4）对比不同产量效果；5）对油藏和流体性质的敏感性研究。油藏数值模拟的作用求实创新超越REALITY，INNOVATION,TRANSCENDENCY2、对已开发油田历史模拟1）证明地质储量，确定基本的驱替机理及驱替类型（是溶解气驱、注水驱、蒸汽驱或是重力驱？）2）确定产液量和生产周期；3）确定油藏和流体特性，拟合全油田和单井的压力、含水（气油比）动态历史；4）指出问题、潜力所在区域。油藏数值模拟的作用求实创新超越REALITY，INNOVATION,TRANSCENDENCY3、动态预测1）评价提高采收率的方法（一次采油、注水、注气、注聚等）2）研究剩余油饱和度分布规律

●研究剩余油饱和度分布的范围和类型

●单井进行调整，改变液流方向、改变注采井别、改变注水层位的效果；

●扩大水驱效率和波及系数的方法；

●回答油田开发中所遇到的问题并致力解决问题的方法。3）评价潜力和提高采收率方向

●确定井位和加密井的位置；

●确定产量、开采方式；

●确定地面和井的设备；

●各种调整开发方案和开发指标对比及经济评价。油藏数值模拟的作用求实创新超越REALITY，INNOVATION,TRANSCENDENCY4、专题和机理问题的研究1）对比注水、注气和天然枯竭开采动态；2）研究各钟注水方式的效果；3）研究井距、井网对油藏动态的影响；4）研究不同开发层系对油藏动态的影响；5）研究不同开发方案的各种指标；6）研究单井产量对采收率的影响；7）研究注水速度对产油量和采收率的影响；8）研究油藏平面和层间非均质性对油藏动态的影响；9）验证油藏的面积和地质储量；10）检验油藏数据资料；11）为谈判和开发提供必要的数据资料。油藏数值模拟的作用求实创新超越REALITY，INNOVATION,TRANSCENDENCY地质工程师油藏工程师数值模拟工程师现场工程师一体化工作模式和协作式解决方案是油藏模拟成功的保证！！油藏数值模拟的流程求实创新超越REALITY，INNOVATION,TRANSCENDENCY制定研究目标选择模拟的目标区域收集数据检验数据的准确性和数据间的相关性选择重要的信息简化数据量，以便于管理用平均压力和地面条件的地质储量设计并建立一个粗化的模型建立平面或剖面模型来估算驱替效率和开采机理，并设计拟函数建立单井模型如果井太多，难以处理，则把井合并为“重叠井”或“虚拟井”拟合储量对一个3D的全油藏模型进行历史拟合对全油藏模型分区进行历史拟合把所有拟合好的部分组合起来进行一个基本的生产预测计算为优化生产进行多重预测计算记录每一个阶段的运算结果如果可能，删去模型中多余的部分油藏数值模拟的流程求实创新超越REALITY，INNOVATION,TRANSCENDENCY确定研究目标（5%）获取检查数据（20%）建立油藏模型（15%）开展历史拟合（40%）进行动态预测（20%）明确目的、确定基本策略、划分可用资源、决定研究所需收集数据、分析数量和品质校正错误、合理转换选择模拟器、测算工作量制定辅助措施、设计网格落实拟合目标、确定拟合质量提出油藏认识、分析剩余潜力制定预测方案、选择优化方法确定控制条件、分析预测指标

合理的时间分配是确保模拟研究质量的必要条件。油藏数值模拟的流程求实创新超越REALITY，INNOVATION,TRANSCENDENCYEclipseOffice网格建立FloGrid(GRID)PVT数据PVTi井数据准备Schedule岩石数据

SCAL井筒模拟VFPi地质模型GocadEV

GEOFRAMEF3D,P3DEclipse

Result3Dview历史拟合

SimOpt生产数据模拟运算分析数据处理ECL工作流程油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY（一）静态资料

1.小层数据表或等值线图（包括砂层厚度、有效厚度（或净毛比）、顶部深度、孔隙度、渗透率等）；

2.地质储量及地层、油藏特点的总结报告；

3.油、气、水高压物性PVT数据；

4.油水、油气相渗曲线数据和毛管压力曲线数据；

5.原始地层压力、温度、压力系数数据；

6.油、气、水分布（原始饱和度）或压力分布或油水界面和油气界面；

7.井位分布图；

8.流体和岩石化验分析报告；油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY（二）动态资料

1.射孔完井报告；

2.井史报告、压裂等措施；

3.系统测压资料；

4.试油、试井和试采资料（压力恢复曲线）；

6.油水井别，调整井位示意图；

7.油井生产（水井注水）数据报表：日产油、日产液、日产气、综合含水、压力累积产油（气、水、液）日注水、累积注水

8.区块综合生产数据统计报表：日产油（水、气、液）、采出程度、综合含水累积产油（气、水、液）日注水、累积注水油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备1。模拟工作的基本信息2。油藏模型3。流体PVT属性4。岩石属性5。油藏分区参数6。初始化计算参数7。输出控制参数8。生产参数油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY设定是进行黑油模拟，还是热采或组分模拟；模拟采用的单位制（米制或英制）；模拟模型大小（你的模型在X,Y,Z三方向的网格数）；模拟模型网格类型（角点网格，矩形网格，径向网格或非结构性网格）；模拟油藏的流体信息（是油，气，水三相还是油水或气水两相，还可以是油或气或水单相，有没有溶解气和挥发油等）;模拟油田投入开发的时间；模拟有没有应用到一些特殊功能（局部网格加密，三次采油，端点标定，多段井等）；模拟计算的解法（全隐式，隐压显饱或自适应）。数据准备—模拟工作的基本信息：油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—油藏模型：模型在X,Y,Z三方向的网格尺寸大小；每个网格的顶面深度，厚度，孔隙度，渗透率，净厚度（或净毛比）；网格是死网格还是活网格；断层走向和断层传导率。油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—流体PVT属性：油，气，水的地面密度或重度；油，气的地层体积系数，粘度随压力变化表；溶解油气比随压力的变化表；水的粘度，体积系数，压缩系数；岩石压缩系数。如果是组分模型，需要提供状态方程。油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—岩石属性：相对渗透率曲线和毛管压力曲线。如果是油，气，水三相，需要提供油水，油气相对渗透率曲线和毛管压力曲线（软件会自动计算三相流动时的相对渗透率曲线）；如果是油，水两相或气，水两相，只需要提供油水或气水两相相对渗透率曲线和毛管压力曲线。油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—油藏分区参数：如果所模拟的油田横向或纵向流体属性，岩性变化比较大，或者存在不同的油水界面，这时需要对模型进行PVT分区（不同区域用不同的PVT流体参数表），岩石分区（不同区域用不同的相对渗透率曲线和毛管压力曲线）或者平衡分区（不同平衡区用不同的油水界面）。另外如果想掌握油藏不同断块的储量或采收率，可以对模型进行储量分区（不同储量区可以输出不同的储量，产量，采收率，剩余储量等）。油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—初始化计算参数：油藏模型初始化即计算油藏模型初始饱和度，压力和油气比的分布，从而得到油藏模型的初始储量。这部分需要输入模型参考深度，参考深度处对应的初始压力，油水界面以及气水界面；油气比或饱和压力随深度的变化；如果是组分模型，需要输入组分随深度的变化。油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—输出控制参数：即要求软件在计算时输出哪些结果参数。比如要求输出模型计算油田的油，气，水产量变化曲线；油田压力变化曲线；单井油，气，水产量变化曲线；单井井底压力变化曲线；单井含水，油气比变化曲线等。油藏数值模拟的数据要求求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据准备—生产参数：

对于已开发油田，这部分的数据量非常大。包括油田每口井的井位，井轨迹，井的射孔位置，井的生产或注入历史（油，气，水产量，注入量，井底压力，井口压力等），井的作业历史等。作业练习求实创新超越REALITY，INNOVATION,TRANSCENDENCY

根据数值模拟数据要求，按照以下八个方面落实准备数据，建立一个项目文件夹：1。模拟工作的基本信息2。油藏模型3。流体PVT属性4。岩石属性5。油藏分区参数6。初始化计算参数7。输出控制参数8。生产参数THANKS谢谢!油藏数值模拟基础培训（第三讲）张世明，2009年4月，职工大学油藏数值模拟基础培训求实创新超越REALITY，INNOVATION,TRANSCENDENCY概念模型数据文件解剖概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCYRUNSPEC（6个必要关键字）GRID（8个必要关键字）EDITPROPS（7个必要关键字）REGIONSSOLUTION（2个必要关键字）SUMMARYSCHEDULE（9个必要关键字）组建一个概念模型概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCYRUNSPEC（6个必要关键字）----RunspecSection-------------RUNSPECTITLEEclipsetestexampleDIMENS--NXNYNZ1011216/--PhasespresentOILWATERGASDISGAS--UnitsMETRIC--WelldimensionWELLDIMS--MaxNoMaxPerfMaxGroupMaxWell/Group20020010100/START1'JAN'1980/1234561—标题，2—模型网格维数，3—流体相，4—单位制，5—井及井组维数，6—起始时间。概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCY1234561—构造顶深，2—网格X方向步长，3—网格Y方向步长，4—模型厚度，5—X方向渗透率，6—Y方向渗透率，7—垂向渗透率，8—孔隙度。2、GRID（8个必要关键字）----------GridSection-----------------------GRIDEQUALS'TOPS'20001101112111/'DX'101101112116/'DY'101101112116/'DZ'21101112116/'PERMX'10001101112116/'PERMY'10001101112116/'PERMZ'1001101112116/'PORO'0.301101112116//78概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCY1234561—流体密度，2—气物性，3—油物性，4—水物性，5—岩石物性，6—油水相渗，7—油气相渗。74、PROPS（7个必要关键字）---------------PropertiesSection------------------PROPSDENSITY934.001000.000000.90/--2PVDG--PRESBGVISC--BARSARM3/SM3CPOISE/PVTO--RSPRESBOVISC--SM3/SM3BARSARM3/SM3CPOISE/PVTW--PrefBwCwVwCvw--BARSARM3/SM31/BARSCPOISE1/BARS155.001.017544.3E-050.450.0//ROCK--RefPressureCompressibility--BARSA1/BARSA155.01.8E-05/SWOF--DatafromrecordKrPc(OW)(ID=5)--SwKrwKroPc/SGOF--DatafromrecordKrPc(OG)(ID=6)--SgKrgKroPc/概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCY121—溶解气油比随深度变化表，2—平衡初始化设定。6、SOLUTION（2个必要关键字）------------------------------SolutionSection--------------------------------SOLUTIONRSVD100087.41300087.41/--EquilibriumdataEQUIL--DATUMDATUMOWCOWCGOCGOCRSVDRVVDACCURACYOPT--DEPTHPRESSDEPTHPCOWDEPTHPCOGTABLETABLENoofLAYERS2000.000200.002020.00.0001410.000.000100/概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCY121—井为置定义，2—井射孔定义，3—井组（集输）定义，4——生产井控制，5—生产井组控制，6—井经济限制条件，7—注水井控制，8—注水井组控制，9—时间步长设定。8、SCHEDULE（9个必要关键字）--------SCHEDULE--------SCHEDULEWELSPECS

/COMPDAT

/GRUPTREE

/WCONPROD

/GCONPROD

/WECON

/WCONINJE

/GCONINJE

/TSTEP50*60/3456789概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE数据文件由以下几个部分组成：必需的RUNSPEC标题，问题的维数，开关，存在相，组成等。必需的GRID描述要计算的网格的几何形状（网格块角的位置），各网格的岩石物性（孔隙度，绝对渗透率等）。可选的EDIT

修改计算的孔隙体积，网格块中心深度和传导能力。必需的

PROPS

油层岩石和流体性质（密度、粘度、相对渗透率、毛管压力等）数据表，是流体压力、饱和度和组分的函数。包括组分计算中的状态方程描述。可选的REGIONS把计算网格分成几个区，用于计算：

§PVT属性（流体密度和粘度）

§饱和度属性（相对渗透率和毛管压力）

§初始条件（静压和饱和度）

§该处流体（该处流体和区中心的流动）

§EOS区域（用于组分计算）如果忽略这部分，所有的网格都放在1区中。概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE数据文件由以下几个部分组成：必需的SOLUTION油层初始条件描述—可以是：

§通过给定流体接触深度得到(potentialequilibrium)平衡来计算

§从前面运行建立的重启动文件读取

§由用户指定每个网格块（一般不推荐）可选的SUMMARY指定每一时间步后要写入摘要文件的数据。当运算结束后要生成某些图形输出（如时间有关的含水率）时才是必需的。如果忽略这一部分就不生成摘要文件。必需的

SCHEDULE

指定要模拟的操作（生产和注水的控制和约束）以及需要输出结果的时间点。垂向流动特性曲线和模拟程序中的参数调整也可以在这部分指定。概念模型数据文件解剖求实创新超越REALITY，INNOVATION,TRANSCENDENCY各部分之间的关联模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY关键字数据输入文件的各个关键字（包括各部分的头关键字）长度最多8个字符并且必须从列1开始。这些所有最多8个字符都是有意义的。与关键字同一行的从第9列开始的任何字符都将被当作注释。关键字数据关键字的数据都应该位于关键字后的新行。数据可以用空格或新行分开，而且在行上的位置也不重要。模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCYSATNUM11111111112222222222/SATNUM10*110*2/=星号两边不能有空格。重复次数关键字之后的数据中，星号可以用来表示重复次数。可以在数据前面加上需要重复的次数和一个星号来重复数据。例如，下面两个例子是相等的：模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY默认值

有些数据项可以默认为内建的默认值。如果可以使用默认值，关键字描述中会有说明。有两种方式可以设定默认值的数值。第一种，通过用一个斜杠（/）提前结束数据，剩下没指定的数值将使用其默认值。第二种，位于斜杠之前的选定数值可以通过输入n*来使用默认值，这里n是要连续使用默认值的数值的数量。例如，3*使关键字数据中下3个数值使用其默认值。数字和星号之间不能有空格。如果一次只有一个数值要使用默认值，必须输入1*。单独一个星号是不够的。

下面是一个例子，这口井中PROD1被给定最小的经济产油量2000，最小的经济产气量默认为0，最大含水率极限0.8，最大气油比和水气比极限默认为无穷大，油井维修程序是关闭最差的连接，所有其他项使用默认值：

WECONPROD120001*0.82*CON//模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY字符串当需要输入字符信息，例如井名（如上面的例子）或记号，可以在引号中随意输入。因此下面两个RPTSCHED关键字是相同的：

RPTSCHEDPRESSURESOILSGAS/RPTSCHED

’PRESSURE’

’SOIL’

’SGAS’/

这些引号通常只在一个名称包含内部空格、以数字开头或包含非字母符号时才是必需的。引号在使用通配符*表示井和组的根名称或井名列表时也是必须的。模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY注释任何以两个字符“--”开始的行都作为注释，将被ECLIPSE忽略。注释行（以及空行）可以插入在数据文件中的任何位置。注释也可以加在数据行的结尾，以“—”开始，但这种情况下注释不能包含任何引号。注释也可以在用斜杠（/）结束数据录入的同一行后面而不用“—”这两个字符。模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY全局关键字有些关键字可以出现在输入文件的任何部分。这些关键字及其功能的简要描述如下。读取和重复输入文件

ECHO和NOECHO关键字打开和关闭向打印文件重复输入文件。初始默认为打开重复。INCLUDE关键字使数据可以从另一个文件中读取。它后面紧跟要获取输入的文件的名称。读取之后，文件关闭，然后从INCLUDE之后的关键字开始继续从主文件输入。下面是一个例子：

INCLUDECASE6G.DATA/模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY激活特殊选项

OPTIONS和OPTIONS3关键字用于分别激活ECLIPSE100和ECLIPSE300的特殊功能中的特殊选项。这些选项主要是临时或试验性的，或者用于向后兼容早期版本的代码。设置消息输出和终止限制ECLIPSE有一个消息处理系统，根据严重程度分六个级别：

1消息2注释

3警告4问题

5错误6程序缺陷（臭虫）任一个消息发生后，它就被记录到简易格式输出，打印文件和调试文件。当特定类型的消息发生一定的次数（即“打印限制”）后，程序就不再打印这种类型的消息。同样，当特定类型的消息发生一定的次数（即“终止限制”）后，程序就自动终止运行。MESSAGES关键字允许用户修改默认的打印和终止限制。模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY结束输入文件

END关键字将先于输入文件的实际结束终止数据的读取。程序将不会重复或处理其后的数据。END可以在包含文件中使用。可以使用ENDINC关键字先于包含文件的实际文件尾结束包含文件并返回到主文件的控制。END和ENDINC都将在实际文件尾自动产生，因此不必每次都使用。模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY全局关键字列表COLUMNS重新设置读取数据输入文件时的左右边界。DEBUG设置ECLIPSE100中的调试输出控制开关。ECHO打开每次运行开始时重复打印数据的开关。END终止读取数据。ENDINC终止读取包含文件并返回主数据文件。EXTRAPMS

要求PVT和VFP表外插法的警告信息。FORMFEED在打印文件中设置走纸字符。GETDATA

从前一步的初始文件或重启动文件导入命名数组数据（GRID,EDIT,PROPS,REGIONS和SOLUTION部分。）。INCLUDE插入指定文件的内容。MESSAGES重新设置各个严重程度级别消息的打印和终止限制。NOECHO关闭每次运行开始时重复打印数据的开关。NOWARN禁止警告信息。OPTIONS激活ECLIPSE100中特殊用法的特殊选项。WARN在被NOWARN禁止警告信息之后重新激活警告信息。模型数据文件说明求实创新超越REALITY，INNOVATION,TRANSCENDENCY123456789012345678901234567890132--此注释表示数据文件的开始--关键字必须从第一列开始写起RUNSPEC可以从第8列开始写注释--这是一个注释--另一个注释EDIT该部分可选PROPS该部分必选DENSITY--OILWaterGas95010000.87/之后可以写注释--这是另外一个注释SOLUTIONColumnsECLIPSE关键字对大小写不敏感--FirstLast33/最后一列以后的任何语句都是注释SUMMARYSCHEDULEEND作业练习求实创新超越REALITY，INNOVATION,TRANSCENDENCY

熟悉一个基础的概念模型文件组成，掌握必要的关键字。有条件的可以运行一下。BASE.DATA油藏数值模拟基础培训（第四讲）张世明，2009年4月，职工大学油藏数值模拟基础培训求实创新超越REALITY，INNOVATION,TRANSCENDENCYRUNSPEC/GRID要点分析RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY主要作用（分配存储空间并启动应用程序选项）在主存储空间中顺序为模拟的各个组成部分分配内存；模拟的组成部分包括井、表格数据、模拟网格及解算器存储栈；一些关键字是必选的，大部分关键字是可选的；忽略的关键字等价于将这些关键字设为默认值。

模型大小：DIMENS

模型流体：OIL,GAS,WATER,DISGAS,VAPOIL,根据油藏实际情况选用不同组合。

模型单位：METRIC,FIELD,LAB

模型维数：WELLDIMS(井，射孔和井组维数），VFPIDIMS(注水井VFP表），VFPPDIMES(生产井VFP表），EQLDIMS(平衡区），TABDIMS(流体及岩石维数）

模型输出格式控制：UNIFOUT,FMTOUT,UNIFIN,FMTIN

特殊功能：API,BRINE,DUALPERM,DUALPORO,MISCIBLE,TRACES

标题：TITLES

模拟开始时间：STARTRUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY一般性的两相PVT相图A区：代表我们通常所说的黑油。线AA：代表泡点以上的油（死油）。线BB：代表初始时高于泡点的黑油（活油）。线CC：代表油气两相混合状态（活油）。线DD：代表一种初始状态接近于临界点的流体线EE：代表初始状态是单相的气体（干气）。线FF：也是初始是油气两相混合状态（湿气）。线GG：代表初始状态是高于露点的挥发油（湿气）。RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY•死油油藏：AA油藏在整个开发阶段压力一直保持在泡点压力以上，在油藏中不会发生脱气。生产油气比是常数（脱气发生在井筒或地面）。

•活油油藏：BB油藏在开发阶段压力会降到泡点压力以下，在油藏中发生脱气，溶解油气比降低，生产油气比增加（自由气会生产出来）。

•干气藏：EE气藏压力在生产过程中不会穿过露点线，气藏中没有挥发油产生。

•湿气藏：FF气藏压力在生产过程中会穿过露点线，气藏中产生挥发油。一般性的两相PVT相图黑油模型可以模拟的流体类型RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY流体类型的关键字选择

黑油模拟器，如ECLIPSE，是不能显式模拟组分的变化的。如果需要模拟气体的分离或油滴的析出过程，则要间接的通过改变溶解气的GOR（Rs）和挥发油的OGR（Rv）来实现。具有这种特征的流体分别称为liveoil和湿气。RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY文件输出格式汇总RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE输出文件格式主要是以文件管理形式，并可定义为有格式和无格式，也可定义为单个文件或合成文件（特别是与时间有关的结果文件）。⑴无格式：在缺省状态条件下，文件存储和输出为不可读，便于程序读取，但转换有困难。⑵有格式：可读，便于编辑、交换，但文件较大。使用定义关键字：FMTOUT、FMTIN。⑶单个文件：在缺省状态条件下，用户可根据需要自行定义某些文件，对中途终止运行的作业，不会影响已做完的计算步，单个文件文件总数可达999个。⑷合成文件：为与时间有关的文件。如：重启文件、计算总结文件。使用定义的关键字：UNIFOUT、UNIFIN。已合成一个文件，对时间步无限定，但对中途终止运算的作业没有用于GRAF的文件。RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE输出文件RUNSPEC部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE输出文件

.GRID或.FGRID:网格文件（NOSIM控制输出）

.EGRID:网格文件,与GRID格式不同，文件要小的多。（用关键字GRIDFILE来控制输出类型）

.INIT或.FINIT:属性文件。（用关键字INIT来控制输出）

.PRT:报告输出。文件很大，模型处理及计算结果详细报告。（RPTGRID,RPTPROP,RPTSOL,RPTSCHED控制输出）

.LOG:后台作业时的输出报告，文件比PRT要小很多。可用于错误检查。

.DBG:Debug文件，一般不用。可用于检查ECLIPSE如何处理输入参数。

.SAVE:用于快速重启。（用关键字SAVE来控制输出）

.RFT:RFT计算结果。（用关键字WRFTPLT来控制输出)

.FLUX:流动边界。(用关键字DUMPFLUX来控制输出)

.Snnnn或.UNSMRY:图形文件输出（在SUMMARY部分定义）

.Xnnnn或.UNRST:重启文件输出（用RPTRST，RPTSOL或RPTSCHED来控制输出）RUNSPEC部分作业练习求实创新超越REALITY，INNOVATION,TRANSCENDENCY1、根据自己实际区块的流体性质，选择合适的关键字；2、改变RUSPEC中文件输出类型关键字FMTOUT和UNIFOUT，看输出文件格式的改变；3、运行概念模型后分析产生的文件类型，找出对应的控制关键字。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY主要作用油藏的几何模型和基本的岩石属性；网格块的大小尺寸和深度，DX、DY、DZ、TOPS或COORD和ZCORN；建立网格间的流体流动连接关系；每个网格的孔隙度、净毛比、渗透率（PERMX/Y/Z或PERMR、PERMTHT、PERMZ）；传导系数修改（MULT）；创建各网格孔隙体积、深度和传导系数数组；数值水体和网格水体的定义；有效网格定义，尖灭控制和最小有效网格孔隙体积定义；GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY最简化数据需求网格大小和深度DX或DXV，DY或DYV，DZ，TOPS

或COORD，ZCORN孔隙度PORO渗透率PERMX，PERMY，PERMZ

或PERMR，PERMTHT，PERMZ净毛比或净厚度NTG或DZNET（默认为1）GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY关于净毛比NTGh1h2WOCOilwater储层总厚:H=h1+h2油层总厚：h1水层总厚：h2假设油层净毛比：ntg1假设水层净毛比：ntg2则储层净毛比：NTG=(h1Xntg1+h2Xntg2)/(h1+h2)地质概念净毛比：NTG’=ntg1模拟概念净毛比：NTG≠NTGGRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY主要作用达西方程：适于地面条件下水平单相流动，其中KA/L项视为表征流体在模拟网格间流动能力的量，即传导系数。考虑净毛比后的传导系数定义：GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY主要作用在离散网格中，X方向传导系数：传导系数不能在GRID部分显示设定，但可以通过其它方式修改。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据排列读取网格数据时是X方向优先，然后是Y方向，再是Z方向网格数据读取规则GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据排列笛卡儿网格：起始点在左上角。在页面上，数值从左至右沿着X方向排列，从上到下沿着Y方向排列。对于3D网格，网格块和节点是按照从左到右，从后到前，从上到下的顺序排列的。网格的起始点并不是一定要从网格块（1，1，1）开始。X，Y，Z轴并不是一定要与I，J，K方向平行。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据排列径向网格数据读取规则GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据排列径向网格：起始点是网格模型的中心点。在R方向上，最内部的一圈是第1圈，然后圈数从内到外依次增加。角是沿着顺时针方向度量的；Z方向保持不变。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY数据排列规则对于每一个网格都有的数据，应用下列读写规则：第一个读写的网格编号为（1，1，1）。对于笛卡儿网格的显示，总是把网格块（1，1，1）放在上、后、左的位置。径向网格的网格块（1，1，1）在最接近网格模型中心的位置。在读取笛卡儿网格数据时，总是X方向循环最快，其次是Y方向，最后是Z方向。在读取径向网格数据时，总是R方向循环最快，然后是THEATA方向，最后是Z方向。所有的GeoQuest模拟软件都使用这种数据读取规则。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格几何形状块中心和角点网格GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCYECLIPSE支持多种网格类型。块中心（BC）网格。角点（CP）网格。每一种网格类型都各有优劣。BC或CP的选择取决于模型的类型。在同一个数据文件中不能够同时使用BC和CP网格。网格几何形状GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格几何形状块中心网格：块中心网格需要知道各网格的顶深及X、Y、Z方向上的网格尺寸。各网格的上、下底面都是水平的，各侧面都是平的且垂直的。所有的网格块都是矩形的。要创建笛卡儿BC网格需要指定关键字TOPS，DX（或DXV），DY（或DYV）和DZ。要创建径向BC网格需要指定关键字DR（或DRV），DTHETA（或DTHETAV）和DZ。以V结尾的关键字指的是向量格式的关键字，是其他关键字的替代。由于每个网格块都只用四个实数来定义，所以BC网格比它响应的CP网格的数据量要少一些。用BC网格建立的简单模型甚至可以不用GRID之类的前处理软件就可以创建。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格几何形状块中心网格表示的含有断层的倾斜构造GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格几何形状考虑相邻的BC网格块。对于倾斜结构，各网格的顶深TOPS是不同的。那么断层两边的网格的顶深TOPS也是不同的。BC网格有效地忽略了倾斜。面和断层的差别。例如，在没有构造图指导的情况下，图中的横断面可以看作是倾斜面，也可以看作是断层。BC网格没有足够的信息来计算相邻网格间的交迭面积，这是因为不知道网格块拐点处的深度。流体通常在相邻网格间流动，所以需要建立相邻网格之间的连接关系。ECLIPSE假设有相邻标号的网格（如图中箭头所连接的网格）之间是有连接的，即使它们之间实际可能并没有，如图中所示的有交叠面的网格间就没有连接。这就是因为在BC网格中不区分断层和倾斜面。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格几何形状角点网格表示的含有断层的倾斜构造GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格几何形状角点网格：角点（CP）网格是建立在坐标线和拐点深度的概念基础上的。一个坐标线定义一个网格柱的一个边。坐标线总是直线，但不一定是垂直的。给出直线上分别在网格上下的两个点的X、Y、Z坐标就可以定义一条坐标线。然后，只要确定其各拐点在各坐标线上的位置就可以确定网格块了。这种定义网格的方式使得网格块可以是任何可能的形状：倾斜表面、断层面、尖灭和剥蚀面都可以正确的表示出来。因为每一个网格块都是用4条坐标线和8各拐点坐标定义的，所以CP网格比BC网格有更大的数据量，而且几乎只能借助于如FloGrid/GRID之类的前处理软件才可以创建。角点网格包含足够的信息来计算相邻网格块间的交迭区域，因为知道各网格块的拐点深度。这意味着，在图中，只有可以看到真正有公共交接面的网格之间才有流体的交换。因此，箭头所指的网格间的流体不可能穿过断层流动。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCYBC（块中心网格）CP（角点网格）网格描述简单网格描述复杂前处理软件不是必须的必须使用前处理软件与其它许多模拟器兼容与很少的其它模拟器兼容很难为非常规构造建立模型能够准确的为非常规构造建立模型网格数据量小网格数据量大无法区分倾斜面和断层区分倾斜面和断层很难如实的建立尖灭和剥蚀面的模型如实的建立尖灭和剥蚀面的模型沿着断层面建立的网格连接不正确，需要手工进行修改。能够正确的模拟沿着断层面层之间的接触关系容易建立径向网格模型若没有前处理软件，则很难建立径向网格模型。角点网格与块中心网格GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格模型格块属性定义网格模型格块属性GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY描述各网格格块的大小尺寸和深度的参数称为几何参数。用来指定几何参数的关键字是：对于笛卡儿网格是TOPS，DX（或DXV），DY（或DYV）；对于角点网格是TOPS，DR（或DRV），DTHETA（或DTHETAV）和DZ；对于角点网格是COORD和ZCORN。描述孔隙度、渗透率的参数称为属性参数。描述属性的参数是：PORO（），PERMX（Kx），PERMY（Ky），PERMZ（Kz）。净厚比可以用关键字NTG（净厚比率）或DZNET（净厚度）来定义。用来显式定义一个网格是否是有效网格的关键字是ACTNUM。对于无效网格，此值为0；对于有效网格则有一个统一的值。网格模型格块属性GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY网格模型格块属性属性关键字是一个网格一个值。有许多种方法，可以用不同的ECLIPSE关键字的组合来来给这些关键字赋值。各属性值认为是网格块中心处的值，而且是网格块内各点处此属性的平均值。这种求平均值的运算是在数据输入ECLIPSE之前做的，而且求平均值的方法是由客户自己决定的。同时注意，虽然有些网格会被排除在模拟运算之外（无效网格），但是仍然要提供计算其孔隙体积和传导系数的数据，这是因为ECLIPSE的内部工具是根据其最小孔隙体积来设定一个网格是无效网格的。ECLIPSE有一个规则，那就是所有提供的数据都必须是显式的，因为ECLIPSE没有提供输入函数形式数据的工具。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性典型格式为：KEYWORDValue1value2value3……..value(NX*NY*NZ)/如何为每个网格格块设定一个属性值GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性如何用Boxes设定网格格块属性值

输入的box是用I、J、和K的范围来定义的，是一个由多个网格格块组成的矩形块。Box可以是网格的一个区域、一层、一列或一行。可以为Box内的网格块设置所有属性值，应该为Box内的每一个网格块（包括区域端点的网格块）都指定一个参数值。一个Box会一直保持有效，直到又读入了一个新的Box或读入了ENDBOX关键字。这个新的box会关闭前一个box，并打开另一个box。默认的box是模型中所有网格格块的集合。BOX--I1 I2 J1 J2 K1 K21 20 1 1 1 10 /PERMX200*100 /ENDBOXGRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性如何用EQUALS设定网格格块属性值EQUALS关键字是用数组的形式对网格的属性值进行操作，可以用作关键字BOX的替代品。EQUILS同样可以用来对第3层的网格格块设置其PERMX值：BOX--I1 I2 J1 J2 K1 K21 20 1 1 3 3/EQUALS‘PERMX’100/ENDBOX不用BOX关键字，EQUILS也可以显式定义box。例如：EQUALS--Array value I1 I2 J1 J2 K1 K2‘PERMX’ 2000 / 默认是当前打开的box--i.e.entirereservoir‘PERMX’ 10 1 20 1 5 1 1 //GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性Box作为EQUALS关键字的一部分被打开，然后被退出EQUALS关键字的第二个正斜杠关闭。为了保证为每一个网格格块都指定了参数值，建议现给全体网格格块赋值，然后再用BOX和EQUILS选择性的覆盖其中某些值。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性如何拷贝网格格块属性值COPY关键字用来在当前的box内从一个数组向另一个数组拷贝数据。如果油田的渗透率是统一的，并且当前的box是整个油田时，可以这样设置其渗透率值：COPY'PERMX' 'PERMY'/'PERMX' 'PERMZ'//如果只有一个特定区域的关键字是各向同性的，则可以把BOX关键字和COPY关键字结合起来使用。例如：BOX--I1 I2 J1 J2 K1 K21 20 1 5 1 1/COPY‘PERMX’ ‘PERMY’/‘PERMX’ ‘PERMZ’//ENDBOX在COPY关键字内部，对数组拷贝的个数没有限制。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性如何加、减、乘、除网格格块属性值<OPERATION>‘KEYWORD’ Value I1 I2 J1 J2 K1 K2/‘KEYWORD’ Value I1 I2 J1 J2 K1 K2/‘KEYWORD’ Value I1 I2 J1 J2 K1 K2//或BOXI1 I2 J1 J2 K1 K2/<OPERATION>‘Keyword’ Value/‘Keyword’ Value/‘Keyword’ Value/ENDBOXGRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENCY如何指定网格格块属性如何用MULTPV对网格格块孔隙体积做乘法对于网格块的孔隙体积可以用MULTPV做乘法操作，例如：MULTPV200*1.01/能够增加前200个网格格块1％的孔隙体积。这种乘法操作是很必要的，例如当要确保拟合的地下流体（FIP）与来自于其它来源的估计相一致时。然而，要注意的是，孔隙体积越大，则这个网格所提供的压力支持度就越高。因此，孔隙体积的修改不应该过大，而且对油藏全区压力的拟合应该在调整孔隙体积之后进行。在历史拟合过程中，有时需要调整孔隙体积来达到拟合一个或多个目标的目的。GRID部分要点分析求实创新超越REALITY，INNOVATION,TRANSCENDENC