CMG软件STARS模块操作手册

关键字：Outprn，wprn，outsrf，gridcart,wsrf,wprn，vamod，vatypeDualpor，difrac,con，mod，djfrac，equalsi，por，matrix，permicpor，prpor，rockcp，thconr，thconw，thcono，thcong，hlosspropcmm，pcrit，tcrit，avg，bvg，Cartesian：直角坐标；defaulted：缺省；文件后缀：dat,lst一、前言概述STARS是三维、四相、多组分、热采、蒸汽添加剂的模拟器。网格系统：直角坐标、变深度/变厚度。分散组分一种组分稳定地分散在各相内，如聚合物、凝胶、微粒、乳化、泡沫。这些概念与灵活的组分物性软件包配合使用，如吸附、堵塞、非线性粘度、分散、非平衡的变量传递。特别是，泡沫有两种模拟方法：A）机理的，可模拟泡沫的产生、传播、破裂。B）经验的，但十分实用，如历史拟合。天然裂缝油藏，有四个模拟天然裂缝油内的流动：双孔隙的DP，双渗的DK，三维多层（从外向内）的MINC，垂向细分多层的VR。过程/机理不同，使用不同的模型。基本近似：裂缝与岩块二种介质。裂缝PV和储量小，是主要的流动通道，岩块恰好相反。自适应隐式STARS有全隐和自适应式两种。对大部分问题，仅有少数网格块要求全隐，绝大多数网格块要求显式。自适应隐式适合于锥进问题和非均质严重的层状油藏及高渗薄层。使用的△与全隐式的接近。因而节省CPU时间1/3~1/2，自适应的选取随时间、空间变化，因而是动态。全隐式的井用全隐式的方法求解PW和井块的变量。对于高度的层状油藏的多层完井，这种方法消除了收敛性的困难。井的控制、约束功能十分强，定速率/流压、经济限制、修井等。矩阵的解法解法称为AIMSOL，不完全的高斯分解作预处理，GMRES加速。AIMSOL是针对自适雅可比矩阵。通常使用解法的缺省参数即可。局部的CYZM网格加密有两种用途：近井加密，如锥进和水平井，静态裂缝的模拟。用户指明要细分网格的区域。块间连接、传导使用解法的缺省参数即可。灵活的网格系统有直角坐标、圆柱坐标、变深度/变厚度坐标，无论二维和三维均可以使用这三种坐标系。9、水层模型两种：小水层，且形状边界清楚，用水区块模拟水层。大的，甚至无限大的水层用半解析的水层模拟。单位三种：公制、油田、实验单位。支持用FORTRAN写成，可在以下机器上成功运行：IBM、CDC、CRAY、HONEYWELL、DEC、PRIME、OULD、UNISYS、APILLO、SUN、HP、FRS。绘图软件名为RESULTS，使用SF2文件系统绘图输出。它也可以作部分前处理，如形成网格以及φKSHD等值图的数字化输入。使用一个转换程序可将SR2文件转换成SRF文件。离散井筒当其水平井筒较长时，由于瞬变性、粘性压力降、不同的多相流方式，因而产生了沿井筒入、产出的非均质性。STARS采用离散井筒，且井筒块与油藏联立求解可解决这些问题。井筒流动为显式，可选取多种流动关系式。地质力学模型包括：塑性变形、剪切膨胀、吞吐的压实驱动、高压注入产生裂缝、近井地层的破裂与出砂。由三个子模型构造的地质力学模型可以处理以上的问题。地质力学模型与油藏模型的连接为显式的。可控体积有限元网格（CVFE）将求解区域离散成许多三角形的单元。细节看：“RESERVOIRSIMULATIONWITHACONTROL—VOLUMEFINITE—ELEMENTMETHOD”，SPERE（AUG1922）349—35F。包含CVEF网格的面可以是IJ、IK或JK面。还可以命名用混合网格，近进为锥进式网格，远井为粗大的CVRF网格。二、指导前言1、初始运算：输入数据文件为*.dat，生成的输出文件有：打印输出文件*.out,SR2的指标文件*irf和SR·的主要结果文件*mrf。重启运行：输入文件有*.dat、*.*irf、*mrf，输出文件有*.out、*irf和*mrf。输入数据的分组有7组输入数据。按可下顺序输入：I/O控制，油藏描述，组分性质，初始条件，数值方法模型，井数据与循环数据。个别的有例外，如I/O的关键字也可以出现在循环数据内。个别组内的关键字出现也有顺序限制。圆括号圆括号内的项是可选的。斜线表示任选其一，如*RANGEHCK（*ON/*OFF）。花括号内装一个数表，如{SWKRWKROW（PCOW）}饱和度函数表。网格数组如ARRAA：*POR或*POR{GRID}组分物性如*KEYWORD{NCOMP}*KEYWORD{NUMY}*KEYWORD{NUMX}冒号表示数据的范围，如il（：i2），1≤il＜i≤NX。组分的指定‘comp-des’组分号，或组分名。标题卡*TITLE1（至多40字符）’SimulationRun=1-1989-01-23’*TITLE2（至多80字符）’DualporosityproblemusingtheMINCoption’*TITLE3（至多80字符）‘ThisisA12X12X10Cartesiangridsystem’标题出现在输出的打印或绘图文件内。2、重启运行重启文件二进制文件，包括初始数据与基本变量，按一定的频率写入，写重启文件是可选的。重启的必要性敏感性计算，或历史拟合；改变井的说明；在作大运算之前，称作一个短的运算看正确性；从一个运算的尾部分出多个不同的运算；重启的含义不需每次从头开始计算，仅需从中间某点开始就可以了。如何作重启用初始运算重启文件，使用*WRST，*WRST的后面跟写重启的时间步数，如*WRST10。重启运算以*RESTARTn开始，如想从上次运算的最后一个时步开始重启，那么置n为空，即*RESTAT，否则用*RESTARTn，n表示时间步数。在重启运算时允许修初始文件的数据！把控制输出打印文件的内容输入数据文件内，不在I/O节内使用*WPRN和*OUTPRN，如果想改变二者的参数，它们还可以再次出现循环数据节内。*WPRN表示写网格数据，井数据，解法数据的频率，如果不写井与网格的数据，那么表为：*WPRN*WELL0*WPRN*GRID0这样可以大大减小文件的长度，节省存储单元。如果以上两行不出现，那么每个时步均要输出井、网格块的数据，文件迅速增大。*OUTPUT后跟输出的项目，如哪些井、网格、油藏的数据要输出。如要输出全部井的分井层指标应使用*OUTPUT*WELL*ALL如只输出全部井的一行摘要，应使用*OUTPUT*WELL*BRIEF。控制模拟结果文件的内容使用*OUTSRF控制模拟结果文件SR2的内容。*OUTSRF出现在I/O节，因参数的改变也可再次出现在循环节内。如果网格与井的数据不出现在输出打印文件内，那么频率数为0。*OUTSRF规定了哪些井数、网格数据、油藏数据要打印。可以对井数据、网格数据分别使用变量表。3、描述网格系统必须的：*GRID、*DI、*DJ、*DK。可选的：*DEPTH、*DTOP、*DIP。这些关键字位于油藏描述节内，在*POR之前出现。三种类型的网格系统，如*GRID*CART10106*GRID*CVRF10106*GRID*RADIAL101015网格尺寸用数组形式输入，如：*GRID*CART101012*DI*CON100.0*DJ*CON100.0*DK*CON100.0DKKVAR25.0250.0340.075.03*402*50如果*KDIR*UP，那DK的数据是从底部往上部排列。4、给定空网格块表示空块有两种方法：*NULL和*VAMOD，二者应出现在油藏描述节。在*NULL后跟一个数组，如果该块定空块则数为0，活块的数为1。如*NULL*NULL1：101：101：311：41：31：30如果*NULL不出现，那么全部块为活块。5、细分网格用*REFINE定义细分的区域，它出现在油藏描述节，且在*NULL和*POR之前。它的后面还要跟细分的要求。如网格10×10×3的基本块（1、1、3）细分为X方向2块，Y方向3块，Z方向2块，表为：*REFINE113into232**每个方向细分最大块数为4对双孔问题，不允许网格的细分。变深/变厚的网格可以细分，但细分的块与基本块的厚度相同，如果要求细分块的属性（如φKSHD）与基本块的物性不同，那么应使用*RG。6、时间步长太小用*OUTPRN*ITER*NEWTON输出内迭代（*ITER）、外迭代（*NEWTON）的信息，对分析太小的△T+分益收敛失败原因有：内收敛迭代失败；牛顿迭代失败，导致切小△T；物质守恒误差过大。如果在输出文件内多处出现“iterationroutinefailstoconverge”RGE”，那么试着采取如下的补救措施：缩小△T，办法是缩小*DTMAX或减小每步要求的变化*NORM*PRESS*NORM*SATUR。增加迭代次数，即放大*ITERMAX用*SDEGREE增大因式分解的级，但会增大存储量。如果对第一或第二次牛顿迭代的内迭代收敛失败，但至少在最后一次又收敛，那么问题并不严重。牛顿迭代收敛失败导致△T致的切小，失败的原因是基本变量的最大变化超过了要求量的2倍。可能会发生非物理的值，如负的压力或负的饱和度，或者超过了给定的最大的牛顿迭代次数。如果问题是由最大变化造成的，如不频繁发生问题就不大。如果出现频繁切小△T，那么应采取如下的补救办法：检查KR和PVT曲线的非线性，曲线应光滑。检查网格步长，以及网格块的物性数据的合理性。检查井的约束，生产井给MINPw，注入井给MAXPw。牛顿迭代次数大造成不收敛，那么增加*NEWTONCYC。如果解振荡，如Sg一会儿〈0，一会〉0，好么采取1和4的措施。置全油藏或油藏的某个部分全隐式解法。尤其变化大的区域。这种情况的例子有：1）没有使用重向平衡的初使化。这可能导致初期有大的P、S变化，即使全部井是关闭的。这种情况发生应用全隐式。2）有气顶。在气顶的底层设置隐式。3）有高K的网格或网格处。造成P变化小，S变化大。置这此块为全隐式。压力收敛客限置为0.1Kpa。如果收剑限比每步变量的实际变化大的多，那么会造成物质守恒误差。容限比变化小一个数量级。使用*MATBLTOL改变模型的灵敏度。在大部分情况下，要求变化（*NORM*PRESS和*SATUR）容限（*CONVERGE*PRESS和*COMVERGE*SATUR）的缺省值是合适的。但对水平井与裂缝问题，应用小值代替缺省值，锥进问题同样如此。7、定义井输入顺序为：*WELL（必须）PRODCER/INJECJOR/*SHUTIN/*OPEN（必须）*INCOMP（注油/报时必须）OPERATE（至少需要一个OPERATE）MONITOR（可选，不能与非权的注入井连用，后接*PERFV*GEO）PERF/PERFV/PERRG（必须）[全部处于井数据节内。]定义井类型有四类：*PRODUCER/*INJECTOR/*SAUTIN/*OPEN如WELL1’MEDRIVERP1’VERT11；WELL2’MEDRIVERP2’VERT1515；WELL3’MEDRIVERP1’VERT55；WELL4’MEDRIVERP21’VERT1010*PRODRCER1：2*INJECOR3：4*MOBWEIGHT*PERFV1：21：31.65*PERFV3：42：31.878、怎样关井与开井关井先用*OPERATE、*MONITOR、*PERFV定义一口井；而后在任何时间*TIME或*DATE就可在井，甚至刚一定义完就可以关井。开井;关闭的井可以在任何时候再一次打开。臂如，将一口吞吐井定义成一注、采两口井，井别的开关转换如下：wellI’Ingjector1’injectormobweight1’operatebhp1000tinjw450qual7perf1**Ijkwi11188well2’probucer1’probucer2operateLiquid1000perf2**Ijkwi11188shutin2time10shutin1time17open2time40stop9、操作和监控约束OPERATE和MONITOR表示对井的约束，至少需要一个*OPERATE，*MONITOR可选。每口井的第一个*OPERATE，称为基本的操作约束。如果违返了监控的操作约束，且带有*CONT，那么监控约束就变为操作约束。如果有多个操作约违反，那么采取的行动是：最严历的是*STOP或*SHUTIN，轻微的是*CONT，（1）生产井生产井的操作：定速率约束（基本操作约束），或定最小的井底流压。如果是油井就定油率或涂率，气井定气率。往后的约束可以是最小产率。如井1有两个操作约束：*OPERATE*MAX*OIL12000.0*CONT*OPERATE*MIN*NHP1500.0*CONT（2）注入井注入井的操作：定最大的注入速率（基本的操作约束），或定最大的井底流压约束。注气井、注水井分别选气率、水率作为约束。如*INJECTOR2*OPERATE*MAX*WATER1000.0*STOP*OPERATE*MAX*BHP2250.0*STOP（3）监控约束我们极力主张监控生产井的GOR或含水，如：*OPERATE*MAX*OIL1200.0*CONT*OPERATE*MIN*BHP2500.0*CONT*MONIOR*GOR15000.0*STOP10、给定井的指数为了给定井的指数，要使用*GEOMTRY/*PERF/*PERFV/*PERFRG。*PERF/*PERFV/*PERFRG必要，仅对*QEOMETRY可选，仅流度加权的生产井或注入井才用它。为了内部计算井的指数，*GEMOETRY必须后跟多个参数。完井关健字*PERF/*PERFV/*ERFRG必须后跟完井的位置和已外部算好的井指数。如果*GEOMETRY使用，那么必须后跟一个完井的关健字。*GEO与完井关键字一同使用。*PERF最适合用于水平井或斜井，当然也可用于垂直井。它的格式为*WELL1’12—09—18—56’*PERF1**ifjfkfwisetn112：41.241*WELL1’12—09—18—56’**radgeofacwfracskin*GEOMETRY*K37524881.00.0*PERF*GEO1**ifjfkfffsetn1·：411如果*VERT位于*WELL的后面，那么就应使用*PERFV。仅在*PERFV内给定完井段K1：K2。如：WELL2*VERT22*PERFV1**kfwisetn2：41.561GEOMETRY应与*PERFGEO/PERFV*GEO/*PERFRG*GEO联合使用！基本块（1，1，3）细分为（3，3，2）个小块，井处于基本块（1，1，3）和细分块（2，2，1：2）内，那么*REFINE332*RANGE113*WELL1*PERFRG1**ifjfkfirjrkrwisetn11、水平井模拟水平井有两种方法：方法1：注入井作为线源，生产井作为线汇。略去了井筒的摩阻压力降和流体的举升效应，（holdepeffects）。对有回流（backflow）的问题常给错误的结果，少量轻微的回流沿且不严重。K的平面差异可引起水平井筒的回流。*WELL*INJECTORor*PRODUCER*OPERATE*GEOMETRY*PERFGEOGEOMETRY和*PERF相当于输入了井的生产能力。运行以上的数据组，观察计算的产率，如果结果不对，那么就去掉*GEO，用*PERF直接输入用户计算的井指数。方法2：使用井筒的离散模型。这个模型既可用于水平井，又可用于垂井。它最适合摩阻压力降、举升效应显著的那些问题。要引用*WELLBORE，井筒作为块的次生孔隙度，因而必须指定相应的压缩系数、岩石类型、Kr表。停止一个模拟运算*STOP出现在数据文件的最后，只能出现一次，停放的时间或日期与它连用，如*DATE19980908*STOP12、建井数据的概要数据顺序如下：*TIME/*DATE是必要的用*DTWELL定义阶段的初始△t井别*PROUCER/*INJECTOR*OPERATE（MONITOR）*GEOMETRPERF/PERFV/*PERFRG*SHUTIN*OPEN可以使用不同的关键字以下是井约束随时间的改变定义新井：在调整现有井的参数之前，使用步（1）（3）（4）（5）（6）。可以用*ATER改变任何基本操作约束，且*ALTER与*ALTER与*TIME/*DATE同时使用。可以调动I/O控制的参数，以及传导率乘子。在以后井的改变中，也可使用*DTWELL和*DTMAX关键字数据输入系统13、关键字数据输入系统关键字系统前言格式：关键字、数据什。如*MAXERROR10，如果关键字不出现，就用缺省值。字符集注释的字符、引号内的字符不检查正误，只原样输出打印。CGM的字符集由大小写的字母、数字0——9、关键字指示符，算术运算符=、+、—、/构成。关键字指示符以下四种情况均有效：*POR/POR/*por/por。**字符串，表示注释行。可以用COMMENT重新定义新定义新的注释符。关键字顺序关键字分为若干组，不同组的关键字按要求顺序输入。某组关键字内的关键字可以按任意顺序输入，除非有特别的说明。LIST、MOLIST、*INCLUDE、*COMENT、*TRADNSLATE、*RANGECHCK可以出殃在数据文件的任何位置。某些如*WRST可以出现我有I/O节，也可出现在井数据节内。关键字分为为必须、可选的、条件必须的。可选的具有缺省值，如关键字不出现就用缺值。串连关键字如必要，一行内可以包含多个关键字。大小写关键字、字母数字串可以大写、小写或任意的混合，文件名按操作系统的要求书写，如IBM要求大写。行的长度最大长度130字符，以后的字符被略去。如需变更，可以在源程序文件KWCOML。INC内设变参数MDLINE，并重新编译源程序。分隔符关键字、数据、字符串用空格、逗号或空行分隔开。字串字串放入引号内，如‘5-35-48-W5’。以下情况是允许的‘THISISTHE“RIGHT”WAY。自己定义关键字如需要，用户可以用*TRANSLATE定义新的关键字的对应，等价关系。数值以自由格式输入数据，实数可以不带小数点、指数用E、e、D、d表示，一个数值内不能有空格。如果要求输入整数，而实际输入带有小数部分，那么打印出错误信息，程序停止运行。以下的实数均为有效：25.040，-3，1.23E+02，0.02D-4，34.e02，+2.3，+.3，-..3以下为无效：34.E0234.E0234.E02有空格，指数带小数。重复计数。23，。23，。23，。41，。27与3*23，。41，。27等价有效整数的范围如若干整数是相邻连接变化的，那么可以用，代替TO，如123456101112等价1：4610：12表格**Swkrwkow（pcow）0.20．01．045．00．28990．0220．676919．030．37780．0180．415310．070．46670．0610．21784．090．55560．1430．08351．800．67820．2890．0123．500．75610．4500．0．101．00001．0000．0．0使用表格的输入选择*INT，可将来源不同的表格自动合并。错误与警告信息数据文件回应打印输出时，带有警告或错误的信息。如果没有使用*NOLIST，那么即使有错误也不打印错误信息。因此，仅当数据文件完全正确以后，才建议使用*NOLIST。*COMMENT注释（可选）格式：*COMMENTTEXT如：*MAXERROR14**Changemaximumnuberoferrors空行（可选）为使输入的数据文件便于阅读，可以加入空行。*RANGECHECK数据的范围检查（可选）格式*RANGECHECK（*ON/*OFF）ON要检查，OFF不检查；无RANGECHECK时要检查。*RANGCHCK*OFF表示非常重要的数据不作范围检查，但重要的数据仍要检查。如要跳过压力变化限定的检查；*RANGECHECK*OFF*NORM*PRESS500*RANGHECK*ON每个数据的变化范围在CGM的内部自动设置。插入文件（可选）格式：*INCLUDE‘FILENAME’*LIST数据文件的列表（可选）LIST/NOLESTLIST与NOLIST可在一个输入数据文件内多次交替使用，缺省要列表。*COMMENT改变注释的定义符（可选）格式：*SOMMENT‘ab’缺省时‘ab’注释可以单占一行，或放在数据行的后面。*TRANSLATE定义等价的关键字（可选）格式：*TRANSLATE‘YOUR—KEYWORD’‘CMG—KEYWORD’目的是定义用户便于记忆新关键字，如：*TRANSLATE‘KX’‘PERMI’*TRANSLATE‘NX’‘NDIVIX’*GRID—ARRAY输入网格性质数组格式：*GRID—ARRAY（ARRY—QUALIFIER）READ—OTIONDATE（ARRAY—MODIFIER）ARRAY—QUALIFIER称为数组的限定符，包括：*MAERIX*FRACTURE*RGBLOCK—ADDRESS—RANGE*WELLBOREBLOCK—ADDREESS—RANGE*ANNULUSBLOCK—ADDDRESS—RANGE*TUBINGBLOCK—ADDRESS—RANGE如果ARRAY—QUALIFIER不出现，那么就是*ALLELEM。网格地址范围（BLOCK—ADDRESS—RANGE）指粗网格块的，形为：il（：i2）ji（：j2）kl（：k2）读入选择为*CON/*IVAR/*JVAR/*ALL/*IJK/*EQUALSIARRAY—MODIFIFR称为数组的修正，指*MOD。离散井筒处理成基本块的局部细分网格，基本块IJK，井筒块地址为ijk/11，如果井筒循环流动，那么油管地址是Ijk/111，环空地址是Ijk/21。如果给块IJK赋值，缺省的这个值也赋给它的全部细分块。*GRID—ARRAY输入细分网格的物性格式：*GRID—ARRAY*RGil（：i2）jl（：j2）kl（：k2）如像：*POR*RG111*ALL（细分为ni=3，nj=2nk=1）08、079、078、081、080760这种方式可为每个细分块输入不同的值，是与基本块的值不同。应与油藏描述节的*REINE联合使用（细分的定义）il（：i2）jl（：j2）kl（：k2）是被细分的基本块的区域。任何的读入选择对*RG有效。数据个数应与细分块的总数相同。*WELLBORE输入井筒的网格数据格式：*WELLBORE（：i2）jl（：j2）kl（：k2）*ANNULUS（：i2）jl（：j2）kl（：k2）*TUBIBNG（：i2）jl（：j2）kl（：k2）*WELLBORE表示将数据分给网块的井筒单元。对循环井就分给环空与油管。*ANNULUS（：i2）jl（：j2）kl（：k2）的环空单元。如果无油管也分给井筒。*TUBING（：i2）jl（：j2）kl（：k2）的渍管单元。仅当*WELLBORE/*CIRCWELL出现，以上三种数组限定符才有效。例如：*SO*MATRIX*CON0。7SOANNULUS1：625*CON0.05SOTUBING1：625*CON0例如：*SO*MATRIX*0.70SOWELLBORE1：625*CON0表示环空和油管内SO=0.0。有两种用法用法1：仅有*WELLBORE，无*ANNULUS与*TUBING。用法2：有*ANNULUS和*TUBING，无*WELLBORE。*CON常数值的数组格式：*CONValue例如：*POR*CON0.16*PERMI*CON100*IJK数组的分区域输入法格式：*JIK{il（：i2）jl（：j2）kl（：k2）value}{}表示一个网格区域给一个常数，可出现多个{}（表），*IJK读入法在循环数据节内为有用，或复写缺省的数据，或复写前给定的数据，如*KRTYPE*IJK1：51：51：14例如，网格为10×10×3，有三种输入法：方法1：*POR*IJK1：101：101：30.2461：51：510.17正确方法2：*POR*CON0.246*MOD1：51：51=0.17正确方法3*POR*CON0.246*POR*IJK1：51：510.17错误*IVARI方向变化的数组输入法格式：*IVARvalue（1）。。。value（ni）例如，nx=3DIIVAR3000.04000.05000.0*JVARJ方向变化的数组输入法格式：JVARvalue（1）。。。value（nj）*KVARK方向变化的数组输入法]格式：KVARvalue（1）。。。value（nk）例如：nz=3*PERMIKVAR2000.050.0500.0*PERMJ*KVAR2000.050.0500.0PERMKKVAR20.040.060.0*ALL全部网格块的数组输入法格式：*ALLvalue（1）。。。value（ni*nj*nk）数的输入顺序按先x后y最后z。例如：ni=10nj=3nk=2*POR*ALL.08.080.81.09.12.1509.097.087.011.15.13408.087.157.145.12.12.135.18.029.047.12.12.154.167.187.121.122.08.08.095.13.12.157.17.18.184.122.084.09.11.12.134.157.157.157.18.18.098.09.0908.09.144.143.123.16.165.102.10.10*Y\EPUAKSI数组的全等输入法（方赂性）格式：*EQUALSI（[*常用于传导率、渗透率、弥散系数。例如：*PERMIMATRIX*CON100.0*PERMJMATRIX*EQUALSI*RERMKMATRIX*EQUALST*2*PERMI*FRACTURE*CON10000*PERMJ*FEACTURE*EQUALSI*PERMK*FRACTURE*EQUALSI*2*MOD数组的修正（条件）格式：*MODil：I2jl：j2kl：k2（+）value（-）（*）（/）（=）或*MOD（*）value（-）（+）（/）（=）紧接修正的数组后面。如果网格区域不出现，则是整个网格。例如，网格*POR*ALL.08.08.081.09.12.15.09.097.087.011.15.134.08.087.157.145.12.135.18.092..047.12.12.154.167.187.187.121.122.08.08.095.13.12.157.17.18.184.122.084.09.11.12.134.157.157.18.18.098.09.09.08.09.144.143.123.16.165.102.10.10*MOD1:31:41+.01521=.13或*MOD*0.95*INT内插表格数据（可选0如果元素所在位置为*INT，那么用相邻元素的线性插值生成。例如：*SWT**SwKrwKrow0.20.01.00.30.05*INT0.4*INT0.70.50.40*INT0.6*INT0.50.70.8*INT0.81.00.01.01.00.0*INT只能出现在函数列，崦不能出现在自变量*SW列内。如果整个列仅有一个非0的数，那么其它元素也具有相同的值。输入输出控制摘要确定输入、输出行动的参数，如文件名、单位、标题、选择、写的频率等。选择表各种输入和输出的文件名只需要输入数据文件名重启运算仅需要一个输入的重启文件名。Irg全部相容文件名缺省。各文件名的给定或缺省是独立的。缺省依赖于输入和输出文件名的根据名。容易作重启运算的引导。输入、输出单位的选取有SI、油田和实验三种单位。个别量的单位可以改变。输入、输出的单位可以不同。基于mol的量，可用质量代替mol。输出文件有如下的写选择：可以用井、网格和数值解的数据。每个量的频率、；输出个数可变。对整个网格提供了长长的变量表。可以控制数组输出的方向排列。某些量可使用特殊的单位，如ppm，pH等。SR2文件具有如下的写选择：可以用井、网格、历史数据。每个量写出的频率及数量（个数）可变。对整个网格供了变量表。某些量可使用特殊的单位，如ppm，PH等。可以用二进制（无格式）或XDR格式写SR2文件。重启有如下的选择：写的频率。重绕频率。以上一个时步或特定时步，读重启。中断有如下的选择：在清除和关闭文件后立即结束运算。在完成当前时步和写重启记录后结束运行。用指令提示用户，要求回答（交互式）。必要的数据本节无必要的关键字，均为可选的。每个关键字均有缺省值可供使用。关键字顺序如果有*FILENAES，应该是第一个关键字，如果有*MASSBASIS，*PARTCLSIEE，*PARTCLWT，必须在*OUTPUT和*OUTSRF之前出现。其它节的使用*MAXTEPS可以出现本节和数值控制节内。本节的关键字可以或不可以出现在循环数据节，有下表：可出殃于数值节不能出现在数值节*MAXERROR*TIELE1*RESTART*SRFASCII*TITLE2*MASSBASIS*WRST*TITLE3*PARTCLSIXE*REWIND*CHECKONLY*MAXSTEPS*OUTSOLVR*INUNIT*OUTSRF*WEL*OUTPRN*OUTUNIT*OUTSRF*SPECIAL*WPRN*PRNTORIEN*XDR*OTSRFGRIS*DIM*WSRF*SR2PREC模拟结果文件系统SR·文件系统SR2包含三个文件：指标结果文件。Irf，主要结果文件。Mrf，重绕结果文件。Rrf。后处理的绘图软件RESULTS和编写报告的软件要使用。Irf和mrf文件。作重启运算也需要三个文件，如果用*rewind写重启记录，那么对重启运算rrf也是必须的。只有上次重绕烈军属入rrf的信息才可以得出（冲掉了老的重启信息）。SR2文件系统与以前的SRF文件系统不兼容，为了从SR2文件系统产生一个SRF文件系统，应该使用工具转换程序SR2RF。内部生成缺省的文件名（文件名包括路径名）有三种类型的文件：输入数据文件。Dat，输出文件。Out，输入重启文件。Irf。通过提示或命令行输入数据文件名，输出文件名为缺省。通过提示，*FIENAME或命令行输入‘输入重启文件名’。使用缺省功能，仅需改变每个运算的*INDEX—IN文件名，就可以完成一系列自展式的重启运算。例1、输入数据文件‘cycle。Dat’，*output的缺省是‘cycle。Out’*INDEX—OUT的缺省是‘cycle。Irf’，*‘runl。Mrf’。例2每一个运算部分用命令“stars。Ese—frunldat”完成，缺省的输出文件是‘runl。Out’，‘runl。Irf’。如果没有命令行的选择—frunl。Dat，程序会提无法用户回答。对第二个运算部分，将‘runl，dat’拷入‘runl2，dat’。在‘runl2，dat’。在‘runl2，dat’中加入*RESTART。执行命令“stars，ese—frun2，dat-frun2，dat—rrunl”。输出文件是run2，out，run2，irf和run2，mrf。事实上，runl，irf为run2，ir的你亲，当你在用RESULTS和报告编写软件（REPORTWRITER）查看run2时，井与特殊历史的第一部分数据就自动在SR2系统的r—unl内找到。对较小的数据文件，适合于如上所述的每个重启运算的单独COPY，但是对较大的数据文件，这就变得不实际了，可以注释掉无用的数据行，还可以将不改动的数据部分*INCLDE的形成插入这样文件就可变得较小了。运行模型STARS可以交互式或批处理式的运行。如果你采取了命令行无自变量的交互方式，那么你就键入提示要输入的文件名。每一次提示是对输入数据文件名，每二次提示是对输入重启文凭名（如果是重启运行的话）。仅当你同时用了*FILENAME和*PROMPT，才提示用户输入其它文件名。如果使用了命令一行自变量，就不会产生提示。执行批处理方式有两种方法，特别适合于多个作业的后台远行：方法1、把全部的键盘提示改为作业运行原本，原本包含了正确的响应。通过命令行一自变量和/或*FILENAME给定全部需要的文件名，这样就不需要提示了。方法2、放*SCREEN在输入数据文件的头部，这样数据就可以直接读入标准的输入装置，如象“stars，ese，〈run，dat”。这种情况下，必须给定*OUTPT的文件名，因为程序不知道输入数据文件名，且无缺省值。Stars。ese命令行一自变量（可选）格式：Stars。Ese（-finput-data）(-rinput-restart)(-checkonly)(-dimsum)(-onestep)Stars。Ese是STARS的请求命令，一般是一个可执行文件的名称。在UNIX内可以是一个局部的文件，一个同文件的联系，或仅通过你的寻找路线的一个进入。Input_data数据的文件名。Input_restart用前一次运行生成的输入重启文件名irf，缺省：如果Input_data不出现，那么将提示用户输入；如果这是一个重启运行，输入重启文件名INPUT—RETART不出现，将提示输入。*FILENAME输入/输出文件名（可选）仅用于覆盖缺省的文件名，或当其未使用命令行一自变量时给定输入的重启文件名。输入数据文件名必须通过标准的输入装置给定（键盘/作业运行Script）或命令行一自变量“_f”。格式：*FILENAME（S）{file_type（name_option）}文件类型有：*OUTPUT*INDEX—OUT*MAIN—RESULTS—OUT*REWIND—OUT*INDEX—IN*MAIN—RESULTS—IN*MAIN-RESULTS—IN*REWIND—IN*GEOMECHOUTname_option是：‘filename’*PROMPT如果文件类型是*OUTPUT，那么name

_option还可允许*SCREEN*OUTPUT是有格式的模拟结果输出文件（out）*INDEX—OUT是无格式的指标结果文件（irf）当其开始重启运行时，将将SR2输入文件的其些信息拷入SR2的输出文件。井的历史数据不拷贝。不抹去输入的SR2文件，直到它包括的信息不再有用为止。*MAIN—RESULTS—OUT是无格式（二进制）的主要结果文件（mrf）*REWIND—OUT以重绕的方式写入的重启数据，称为可重绕的结果文件（irf）。*INDEX—IN是指标结果文件，从它读入模拟结果和重启记录，此文件仅对重启运行是必要的。*MAIN—RESULTS—IN是主要结果文件，从它读入模拟结果和重启录（二进制），仅对重启运算它才是必要的。**REWIND—IN是可重绕的结果文凭，从它读入重绕、重启记录（二进制），仅对重启运算它才是必要的。GEOMECHOUT地质力学模型生成的有格式输出文件。如果GEOMECHOUT不存在，或者是PROMPT使用了，那么就把这个输出写到主要输出文件（用*OUTPUT）。‘’空串，将使用一个内部生成的缺省文件名。‘filenam’文件名，可接受的文件名依赖于操作系统。*PROMPT提示用户通过键盘打入文件名。除去INDEX—IN外，全部文件类型都有一个内部生成的缺省文件名，为了使用它，在提示下打入一个空白的反应（大概是立即回车）。*PROMPT提示用户通过键盘打入文件名。除去INDEX—IN外，全部文件类型都有一个内部生成的缺省文件名，为了使用它，在提示下打入一个空白的反应（大概立即回车）。*SCREEN对这种文件类型，数据被送到标准的输出装置（屏幕/作业日记文件）。缺省：输入数据文件无缺省的文件名。用户必须提供，或是通过标准的输入装置（键盘/作业运行Script），或命令行一自变量“_f”提供。如果任何其它的必要的文件名来通过*FILENAME（包括*PROMPT）给定，那么就命名用内部生成的缺省文件名。条件：如果*FILENAMES存在，那么它必须是数据文件的第一个关键字，如不是这样，就提示用户输入所输入所需的文件名，而任何以后的*FILENAME将不识别。打开全部带有‘UNKNOWN’状态的输出文件，因而不能避免重复的写入（OVERWITING）。依赖使用的操作系统，一个老的文件将重复写入或在后面附加。打开全部带有‘OLD’状态的输入文件，因此它们必须存在。*CHECKONLY检查错误的扫描方式（可选）格式：*CHECKONLY如果*CHECKONLY存在，作业只检查数据文件的错误，包括语法的和超范围的错误。如不存在，则要作时步的计算。*TITLE1*TITLE2*TITLE3方案和情况的标识（可选）格式：*TITLE1‘字串’max40*TITLE2‘字串’max80*TITLE3‘字串’max80*CASEID‘字串’max8*CSDATE‘字串’max10TITLE1最多40字符，TITLE2和*TITLE3最多80字符，它将出现在打印输出文件和SR2文件内。*CASEID最多8字符用于SR2文件标识曲线。*CSDATE最多10字符，模拟运算的日期，如缺省操作系统提供当天的日期。它产出现在I/O控制节内。例如：*TITLE1’DUALPOSITY/DUALPERMEABILTITYRUNNO。1’*TITLE2’RUNBYABSTARFF，DEC16，1988，C、DCO*TITLE3’4200gridblocks；var，thickeness’*CASEID’NOGAS’*INUNIT输入/输出数据的单位（可选）格式：INUNT（SI|*FIELD|*LAB）{*EXCEPTqnty_nounit_no}OUTUNIT（SI|*FIELD|*LB）{*EXCEPTqnty_nounit_no}定义：*EXCEPT、例外，它允许些量不命名用这种单位。qnty_no下表内的量编号。Unit_no下表内的单位骗号。缺省：如果*INUNIT不存在，则为*SI。如果OUTUNIT不存在，则与输入单位相同。本手册各关键字的缺省值而无单位，这个单位是*SI。解释：SR2文件内的数据使用的是STARS的内部单位，与输入单位无关。但是，指标结果文件。Irf数据使用了输出单位。例：INUNITSI*EXCEPT22**用0F，not单位表量名称*SI*FIELD*LAB1、时间daysdaysminutes2、温度degCDegFdegC3、压力kpaPsiKpa4、长度mFtcm5、体积m3Ft3Cm36、渗透率mdMdMd7、质量kgLbKg8、mol质量（质量基础）Gmole(kg)Lbmole(lb)Gmole(kg)9、粘度cpCpCp10、能量JoulesBtuJoules1、流量M3BblCm312、表面张力Dyne/cmByne/cmDyne/cm改变单位的选取名称量号0123时间1daysHrMinYr温度2degKDegCDegFDegR压力3kpaPsiAtmBar长度4mFtCm体积5M3Ft3BblCm3渗透率6M2DarcyMicro-m2Md质量7kgLbMol质量8GmoleLbmol粘度9Kpa-dayKpa-hrCp能量10jBtu流量11M3Ft3BblCm3表面张力12Kpa-mN/MDyne/cm*MASSBASIS质量基础指示符（可选）格式：*MASSBASIS解释：某些化学驱过程，要求用质量分数去代替mol分数，缺省单位为mol。例如，组分骤合物的分子量极大，对应的mol分数小，基于mol分数加权的混合规则可能不再合适。*MASSBAIS的作用：输入的mol数据解释成质量数据。组分的性质是单位质量，k值定义为相的质量分数，而不是mol分数比。相的组成为质量分数。本手册输入每个组分的数据进，包含了它的单位的一般定义。例如，密度单是molarmass/volume，molarmass解释为mole/massbasis。对SI单位，密度单位为gmol/m#，WJGJVKET*MASSBASIS，那么它为kg/m#。唯一的例外是分子量的定义，对SI单位是mass/mol/mole，即KG/GMOL。有组分汽化的地方不要使用*MASSBASIS选择，这是因为k=汽/液，必须用mole分数定义。对蒸汽驱过程特别是这样。*MAXERROR错误信息的最大个数（可选）格式：*MAXERRORnnmnum=1～100，缺省为20。数据检查错误，如果错误个数超过这儿的num，那么模拟停止运行。*RESTART开始和停算的时间步骤格式：*RESTART（nstart）*MAXSTEPSnstop定义：nstart重启运算开始计算的时间步数。Nstop是停止计算的时间步数。缺省：如果*RESTART不存在，不需要读重启良悄，那么第一个时间步数为1。如果它存在，但无nstart，就命名用*INDEX—IN内的最后一次重启录。如果*MAXSTEPS不存在，那么nstop=9999。解释：为了调试输入的数据文件，可以采用*MAXSTEPS1。*WRST写重启记录（可选）格式：*WRST（*TIME/freg）*REWIND（num）目的：在输出重启文件时，二者控制着写记录和重绕记录的频率。定义：fweg为整数，fweg是频率，表示时间步数。*TIME或为*DATE。台果freg=0不产生重启记录。TIME表示在每个TIME或*DATE时刻，要把网格的数组长写入重启文件。Num重绕重启文件的时步数。每隔num=0不重绕。如果num=1，那么仅保留最后一次重启记录，记一次重启又重绕一次。缺省：如果*WRST不存在，就不产生重启文件。如果后面无*TIME或freg，就认为有*EINE。如果*REWIND不存在，不后果绕重启文件。如果*REWIND存在但无num，就认为num=1。条件：以上二关键字可以出现I/O节与循环数据节内。*WPRN打印输出的频率（可选）格式：*WPRN（*GRID/*ITER）（*TIME/freg）目的：控制写入打印输出文件的频率，打例子文件名用*OUTPUT设置。定义:*GRID输出井和网格的数据.*ITER输出井墒要、迭代、守恒、数据。Freg每freg步写一次，为0不写。*TIME在*TIME或*TATE写也结果。条件：*WPRN可出现在I/O节或循环数据节。*OUTPRN打印输出文件的项目（可选）格式：*OUTPRN*GRID{*ALL|*NONE|*EXCEPT|*REMOVE|（special_untit）grid_var}*OUTPRN*WELL{*ALL|*BRIEF|NONE|well_var}*OUTPRN*ITER{BRIEF|NEWTON|UNCONV|TSS|NONE}*PARTCLSIXEvol*PARTCLMLWTcmm目的：给定写入输出的打印文件的项目。定义：*GRID指定风格变量写出。网格变量表为：PRES油相压力SW水饱和度（也可以WATSAT）SO油饱和度（也可以OILSAT）SG气饱和度（也可以GASSAT）TEMP温度Y：气相的组分组成X：油相的组分组成W：水相的组分组成Z、所有相的组分组成SOLCONC固体组分的浓度HEATLOS（已不使用）OBHLOSS热损失率VPOROSVOID孔隙度FPOROD流体孔隙度POREVOL孔隙体积VISO水粘度VISG油粘度VISG气粘度KRW水相相对渗透率KRO油相相对渗透率KRG气相相对渗透率COW水油抟管力PCOG油气毛管力MOLDENW水相mol密度MOLDENO油相mol密度MOLDENG气相mol密度MASDENW水相质量密度MASDENO油相质量量密度MASDENG气相质量密度RFW水相阻力因子RFD油相阻力因子RFG气相阻力因子FRCFLOW相的分相流ADSORP组分吸附KRINTERkr的插值值IFT局部表面张力CAPN局部PC数LOGIFT表面张力的自然对数LOGCAPNCAPN的自然对数FLUIDH流体的自然对数FLUIDH流体的热焓WATERHEAD到等价水柱项的深度（看*DTOP）AQWATCUM水层的净水浸量AQWATCUE水层的水浸速率AQWATCUM水层的净热浸量AQWATCUE水层的热浸速率以下的量受Special_unit的控制：VLKVCMP关键组分的成组，用于计算k=汽/液。LLKVCMP关键组分的成组，用于计算k=液/液。VISCCMP关键组分的成组，用于油水非线性混合粘度。ADSPCMP关键组分的成组，用于kr的计算。PLPMCMP关键组分的成组，用于kr的计算。通常y、x、w、x的单位为mol分数，如用*MASSBASIS则为质量分数。Special_unit物殊组分性质，如VLKVCMP，LLKVCMP，VIDCOMP，ADSPCMP，RLPMCMP在报告中用特殊单位。MOLFRmoL分数MASFR质量分数ppm百分之几VOLFR体积分数MOLAR体积克分子的浓度PHPH=14+LOG（体积克分子浓度）NUM单位相体积内的质点数用*SOLCONC或ADSORRP标识的固体吸附浓度单位：MOLEmoles/孔隙体积MASS质量。孔隙体积VOL固体体积/孔隙体积NUM质点/孔隙体积每种物性可以使用一种特殊的单位，例如：*OUTPUT*GRID*PPM*ADSPCMP*RLPMCMP造成*ADSPCMP组分为PPM，*RLPMCMP组分为number密度*WELL指定的井变量写入输出文件，井变量表：LAYPWF层号和BHPLAYPHASE层相率、累积LAYGAS层气率、累积LAYGAQ层油/水率、累积LAYGAS层气率、累积WELLCOMP井组分/相摘要TYPECONP井组分/相摘要TYPECOMP注、采组分/相的摘要*ITER打印以下迭代结果BRIEF收敛统计NEWTON牛顿迭代摘要TSSNEWTON+时步选择UNCONVTSS+不收敛变量的细节（仅用于调试）*ALL选取此表的全部项目*NONE无项目。*EXCEPT除去以后的例外变量表，其余项目全选。*REMOVE从表中去掉以后的项目VOL固相，吸附或重要的拖带组分的一个质点体积，（m3/ft3/cm3），用于计算*NUM的数值密度。缺省为1。E-11cm3，对应为半径=1。33E-4cm的地球。

CMM固相、吸附或捕集物质一个mol的质量（kg/gmole）此量仅对（MASS命名用，无缺省值。缺省：可选的关键字，如果它不存在，那么缺省为：*OUTPRN（*GRID/*WELL）*NONE*OUTPN*ITER*BRIEF如果*PARTCLSIXE不存在，vol=1，E-11cm3如果*MASSBASIS不存在，浓度的缺省单位为mol。如果*MASSBASIS存在，浓度的缺省单位为mass。条件：此关键字出殃在I/O、循环数据节。仅当*MASS出现，*PARTCLMLWT才是必须的。例如：*OUTPRN*GRID*SW*SO*SG*PRES*TEMP*PARTCLMLWT13**COKEmw=131b/lbmole*OUTPRN*GRIDMASSSOLCONE*WSRFSR2的输入频率（可选）格式：*WSRF（*GRID/*WELL）（*TIME/freg）目的：控制*OUTSRF内项目写入SR2的频率。定义：*GRID控制井数据网格数据的写出频率*WELL控制井数数据的写出频率。Freg频率，时间步数*WSRF*GRID0不写出*WSRF*WELL0不允许*TIME在每一个TIME或DATE时刻写出缺省：如果无*WSRF*GRID，就缺省为*WSRF*GRID*TIME如果无*WSRF*WELL，就缺省为*WSRF*WELL1。条件：*WSRF可出现在I/O，循环数据节内。*XDR模拟结果文件的项目（可选）格式：*OUTSRF*WELL{comp-untit|*DOWNHOLE|*BLOCKP|*COMPONENT(*ONE|*ALL|comp_list)|*LAYER(*NONE|*ALL|WELL_LIST)}OUTSRF*GRID{(*ALL*NONE|*EXCEPT|*REMOVE|(special_his)*SR2PREC(*SINGLE|*DOUBLE)*SRFASCII目的:标记要写入模拟结果文件的项目.定义:*XDR可以用外数据的表示,或以用户平台的二进制格式写入二进制数据文件。使用XDR后，就允许在一个平台上写成SR2二进文件，又可在另一个平台上读出SR2二进文件。例如，在UNIX工作站的服务器上生成SR2文件，然后在一台PC机上用RESULTS或PEPORTWRITER读出。如果SR2是XDR格式，那么*XDR将出现在指标文件，irf的顶部附近。XDR是一个非标准的选择*WELL全部井的数据写入SR2文件。Com_unit除去体积之外，在质量和/或mol单位存储井的动态，comp_unit可以是以下之一或之二：*MASS用质量单位存储*MOLE用mole单位存储*DOWNHOLE除去地面条件外，还以井底条件写入井的动态，使SR2文件变大。*BLOCKP如与*LAYER同用，输出每个井块层的P。*COMPONENT以组分的形成写出井的动态，缺省为*NONE，*ALL表示全部组分，增大SR2文件。*LAYER写出井组层的动态，缺省为*NONE，*ALL表示全部井，也可输入一个井名表。*GRID网格数据写入SR2文件。Srf—varsrf—var表为：由*OUTPRN提供的grid—var表全部（但要去掉FRCFLOW，POREVOL，POREVOL，增加WATFRFL，OILFRFL，GASFRFL），还有以下的变量：KVALXW组分的k=气/水（y/w）KVALYX组分的k=气/油（y/x）KVALXW组分的k=油/水（x/w）KVALWX组分的k=水/油（w/x）CMPDENW水相的组分质量密度CMPDENO油相的级玢质量密度CMPVISO水相的组分粘度CMPVISO油水相的组分粘度CCHLOSS热损失PERMIX方向的kPERMJY方向的KPERMKX方向的KWATMOBI方向的水相流度OILMOBI方向的油相流度GASMOBI方向相流度TOTMOBI方向的总流度WATFRFL水相的分相流OILFRFL油相的分相流GASFRFL气相的分相流VISOVMP油粘度的非线性混合中，关键组分的组成。VISWCMP水粘度的非线性混合中，关键组分的组成。IFTCMP不使用KRSETNKr数据的组分STRESIX方向的有效应力STRESJY方向的有效应力STRESKZ方向的有效应用力STRESSHy—Z或r—Z内的剪切应力STRESMYP最大主应力STRESMNP最小主应力STRNEPL有效塑性变形Special—unit特殊单位与*OUTPRN描述的相同。*ALL选取表内的全部项目。*NONE不选择取表内的项目。*EXCEPT除去后面列出的项目，表内其余项目均选定。*REMOVE从以前的表内去掉后继项目。*Special—his。。。。缺省：如无*OUTSRF*GRID等价于*OUTSRF*GRID*NONE如无*OUTSRF*WELL，结果为：仅写出体积仅写出地面条件仅写出相的仅写出井的如有*OUTSRF*WELL，结果为：如无*MOLE和*MASS，仅写体积。如无*DOWNHOLE，仅写地面条件。如无*COMPONENT，不写组会的。如无*；AYER，不写层的。如无*OUTSRF*SPECIAL，不写特殊历史。条件：*OUTSRF*GRID可出现在I/O和循环数据节内，对*GRID的缺省频率是*TIME。*WELL和*SPECIAL只能出现在I/O节。解释：例如*OUTSRF*GRID*OILSAT*GASSAT*WATSAT*PRES*PARTVCMLWT13**COKEMWIS13IB/IBMOLE*OUTSRF*GRID*MASS*SLXONC*OUTSRF*SPECIAL*BLKVAR*MASS*SOLCONC620*OUTSRF*SPECIAL*BLKVAR*MKLE*SKLCONC620*PRNTORIEN网格数组打印的方向（可选）格式：*PRNTORIENirotatijkord*PRINT_REF（*ON|*OFF）目的：提供打印的方向，提供细分网格的打印。定义：irotat打印的轴旋转标志旋转标志行号列号平面0（最紧凑的打印）1IJK2IKJ3KIJ4KJI5JIK6JKIIjkord轴的反向标志反向标志行号列号平面0底层处于页的底部1正常正常正常2正常正常反向3正常反向正常4正常反向反向5反向正常正常6反向正常反向7反向反向正常8反向反向反向*PRINT—REF用户可以或不可以用*OUTPRN*GRID打印细网格和基本网格的值。*ON可以打印，*OFF不能打印。缺省：如无*PRNTORIRN则为*PRNTORIEN00,给出最紧溱的打印方式,油藏底层在页底.如无*PRINT—REF，则为*PRINT—REF*ON。如果有*PRINT—REF，但后面页*ON或*OFF，则为*PRINT—ON。解释：缺省的打印方向：在打印纸上从左向右为I，从上往下为j，但后面页*ON或*OFF，则油层，从上一个xoy平面到下一个xoy平面上为k，油层的底层处于页底层处于页底。从左至右I增大，从上往下j增大，从上平面往下平面k增大，某些情况下可以用irotat改变打印的方向。通常，在一行内从左至右指数增大，一列内以上往下指增大，平面与列的增大方向同，可以用ijkrd改变指数增大的方向。提示：*PRNTORIEN造成网格数组以标的顺序写向输出文件，与网格数组输入选择*ALL对应。为使一个运算的输出文件作为另一个运算的输入文件，应使用*PRNTORIEN，11，拷贝和粘巾要求的数据到新文件去，而且自动抹去注释“k=、“j=”、“I=”。*OUTSOLVER矩阵解法的打印出（可选）目的：控制解法软件包AIMSOL结果的打印输出。格式：*OUTSOLVER（*ON/OFF）缺省：*OUTSOLVER*OFF解释：许多有关的量（职块间的连接）的个数与维值可以打印。这些个数用于建立公共存储，公共存储优化了可用的存储空间。还找印出了内迭代残量减小的细节。*INTERRUPT收集UNIX信号2中断（可选）目的：当其检测出中断信号后，给定应采取的行动。格式：*INTERRUPT（*INTRRACTIVE/*STOP/*RESTART—STOP）定义：*INTERACTIVE提示用户以交互方式打入指令，用户可选取*STOP或*RESTART—STOP。*STOP立即结束运算，当前时步未算完，但要关闭输出文件，为的是防止文件的破坏。*RESTART—STOP算完当前时步，写出用*OUTPRN和*OUTSRF给定的输出，定一个重启记录，停止运算。缺省：如无*INTERRRUP，或有但二级子关键字无，则为*INTERRUPT*INTERACTIVE。解释：可以用两种方式传送一个UNIX中娄（信号2）：手动同时按下Ctrl—c，可中断当前的交互作业用带有ID“pid”的命令“kill—2pid”可中断后台作业。中断用于失败一个作时业，文件可以正常的关闭，基些平台不清洗输出文件缓冲区，无中数处理，因些会丢失输出文件。*DATAECHO详细的数据回应控制（可选）目的：控制输入数据文件详细回应格式：*DATAECHO（*ON/*OFF）缺省：如无*DATAECHO，或有但无二级子关键字，则为*DATAECHO*ON解释：回应的数据构成主要输出打印文件的一部分分，包括网格定义，组分性质，岩石—流体性质，初始条件、数值控制等。回应打印增大了文件尺寸，待文件正确后可以去掉回应。油藏描述摘要本节数据的分类：模拟网格与网格的细分天然裂缝油藏井筒离散岩石性质盖层热损失井筒热损失热力水层关键字顺序*GRID是第一个关键字，上述（1）——（3）的关键字在*POR之前出现。网格的选择STARS支持如下两类网格：有限差分网格（FD）笛卡尔的、径向、变深/厚度。可控体积有限元网格（CVFE）对FD网格的选择，*GRID、*DI、*DJ、*DK是必须的，*NINEPOINT、*DIP、*REFINE、*VAMOD、*NULL是可选的。对变深网格，*DTOP可以与初始化的*DWOC、*DGOC和*REFDEPTH同时使用。对CVFE网格选择，*GRID、*ELEMTONODE、*COORDSI、*COORDSJ、*COORDSK是必须的，*DIP、*CVRFHYBRID是可选的。在使用CVFE选择时，以下的选择目前不能使用：。裂缝油藏。离散井筒。热力水层建议用CMG的交互式的网格生成器建立CVFE网格数据。裂缝油藏选择裂缝油藏选择为：*DUALPOR/*DUALPERM/*MINC/*SUBDOMAIC/*VERTNOSEG每个方向的裂缝间距为*DIFRAC、*DJFRAC、*DKFRAC当前不能使用CVFE网格选择。离散井筒的选择将井筒离散成若干网格块，这些网格块的流动方程与油藏流动方程联之。井筒离散用*WELLBORE实现，可后加子关键字*CIRCWELL表示环空与油管和循环流动，当然环空与油管也对直井部队发，井微离散可从油层通过上盖直到地面，非油层部分仅附加很少的风格块。从井筒到盖层热失可用*AQUIFER*WELLBORE模拟。当前井筒离散还不能与CVFE网格同时使。混合网格的离散井筒为了有效地模拟单井的SAGD（SteamAssititedGravityDrsinage），井筒必须直接与它上、下的网格块相连，不的是允许蒸汽上浮和液体的下沉。将离散的井筒加在混合的细分网格内即可实现。岩石性质孔隙度、渗透率和传导率的修正用*POR、*PERMI、*PERMJ、*PERMW、TRANSI、*TRACSJ、*TRANSK、*TRANSIJ+、*TRANSIJ-、*TRANSIK+、*TRANSIK-。要使用限空间*MATRIX和*FRATURE区分岩块和裂缝。当其细分网格用了*REFINE、*MINC、*SUBDOMAIN、*VERTNOSEG、*WELLBORE选择时，要加限空间*RG。岩石类型*ROCKTYRE、*THTYPE是岩石类型数组。岩石类型可按物性分组：压缩笥：*CPOR、*CTPOR、*CPORPD、*PORMAX热物性：*ROCKCP、*THCONP、*THCONW、*THCONO、*THCONG、*CPD盖层热损：*HLOSSPROP、*HLOSST、*HLOSSTDIFF盖层热选择对半解析的无限厚盖层热损模型来说，热损的方向与盖层的热物性用*HLOSSPROP、*HLOSST、*HLOSSTDIFF给定。对CVFR网格，要用*BOUNDELEM、*BOUNDNODE。井筒热损失选择指注蒸汽时井筒的热损失，须要最少的关键字为*RJI、*PTO、*RCI、*EH、*DEPYH，且满足rti〈rto〈rin热力水层选择水层或热力水层选择为*AQUIFER。水层的位置、方向、性质给定如下：AAQUIFER（*BOTTOM/*BOUNDARY/*REGION/*WELLBORE）*AQGEOM（*RADIAL/*RECTANG）*AQH*AQPOR*AQVISC*AQRCAP*AQRCND*AQCOMP*AQPERM如是裂缝油藏，仅裂缝与水层连结。当前水层选择不能与CVFE选择时使用。井至地表的选择在模拟从完井段到地表的井筒流动时，具有如下的油藏与井的组合选择：用*GRID和*DK定义地表到定义地表完井段之间的盖层网格快。臂如1600m厚的盖层，用16个100块模拟，只有两列井筒块参加时步的数值计算。对于*KDIR*UP、盖层块的编号在油层之后：对于*KDIR*DOWN在油层之前。除去包含井筒的两纵列网格块之外，盖层的其它块通通作为死块用*NULL标识。用*POR为其它盖层（非1、2列）列指定0孔隙度除去油藏的完井段外，在网格的定义部分用*WELLBORE定义穿过盖层块的离散井筒。典型地，它是一个斜井筒：油怪部分为水平，且垂直到地表。命名用*AQUIFFER选择，为盖层网格块列内的每块外面的区域指定解析热损。每个层的导体处理成径向块，这些径向块的中心即是井筒，按等体积算出各块的半径。热量从井筒向四周各块传递，使用仅热损失的水层模型，热量成径向地流向上盖层。在地表井的源汇项与离散井筒联系在一起（源汇项只在井口处）*GRID网格类型（必须）目的：*GRID标志油藏条件数据输出的开始。格式：*GRID*CARTninjnk（（*KDIR）*UP|*DOWN）GRID*RADIALninjnk（*RWrw）（（*KDIR）*UP|*DOWN）*GRID*CVFEneletnnodtncvlyr（*ONIJ/*ONIK/*ONJK）（*KDIR）*UP/*DOWN）定义：*CART长方形的笛卡网格*RADIAL径向圆柱网格*CVFE可控积有限元网格ninjnk三个方向的网格块数neletCVFEW全部层的三角形（单）元的总数mmodtCVFE全部层可控体积节点的总数ncvlyrCVFE网格的层数*RWRW最里面一个网格块摧半径（m），缺块省值不能为0。为了精确地复制