CMG中IMEX中文操作手册(二).doc

目 录 imex 的介绍1指导段5关键字输入系统中的数据段5如何建立数据文件的文档方式6如何执行重新启动运行7控制打印输出文件的内容8图形文件(sr2)的内容控制9网格系统描述10局部加密网格描述11双孔/双渗模型的使用12死结点的输入17水区选项的使用18拟混相选项的使用19注聚合物模型的使用20程序执行时间长或时间步长太小问题21单相油藏的模拟23水平井25垂直平衡计算25定义多个 pvt 区29井的定义30井的类型定义31如何关井以及重新开井32亏空填充(voidage replacement)34井筒流动模型的使用36操作及监测限制37井指数的输入39中止模拟运行41在井的列表中使用通配符41关于井数据段设置的指导42在循环数据段内可用的其他段关键字44关键字数据输入系统46关键字系统介绍46注释行 (可选择)52空行52包含文件53控制数据文件列表54基岩网格性质的输入57裂缝网格性质的输入57加密网格性质的输入58由 i 方向确定 j 和 k 方向的数据59常数值数组59以 ijk 方式输入数组60数组输入值沿 i 方向变化61数组输入值沿 j 方向变化62数组输入值沿 k 方向变化62大多数或所有网格的值都是不同的63输入/输出控制段67基本网格定义99k 坐标方向105i 方向的网格尺寸107j 方向的网格尺寸108k 方向的网格尺寸110网格的顶部深度112加密网格的位置123双孔介质minc方法126死结点标识符128孔隙度130岩石压缩系数131渗透率134压力影响函数147流体组分性质数据段153流体模型153油藏温度154油和气的 pvt 表155作为压力函数的油压缩系数168密度173水相粘度179原始气油比180聚合物粘度混合（条件）183岩石流体数据段186相对渗透率表186油水相对渗透率表187气液相对渗透率表194评价三相共存时油相相对渗透率的方法（任选）203初始条件数据段209初始条件标识209油藏初始油相压力212初始泡点压力212初始泡点压力与深度关系213初始含油饱和度215初始含水饱和度216初始聚合物浓度216参考深度和参考压力217油水界面深度218数值计算方法数据控制段225井数据段241井的改变日期241井的改变时间242设定井底流压初始化的频率248附录320151 输入/输出控制段通常情况下，用户不必修改关键字的任何数据就能满足使用要求，所以*dim通常被省略掉。但对于下列两种情况，要注意改变*dim的参数。一种情况是存在尖灭层。当输出文件计算的最小维数与实际情况相差很大时，要在数据文件中设置*dim *active_blocks 。另一种情况是有效网格数占总网格数的比例很小（小于）使用*dim *active_blocks关键字增加有效网格节点的维数。输入时该数据段位于输入输出控制段的前面。*filenames输入/输出文件名 (可选择)_目的： 确定输入和输出文件名。格式 *filename(s) file_types name_options 这里 file_types 是下面的一种或多种: *output *index-out *main-results-out *rewind-out *index-in *main-results-in *rewind-in 对于 file_types 是 *index-out, *main-results-out, *rewind-out, *main-results-in 和 *rewind-in， 必须确定下面的 name_options 之一。 filename *prompt 对于 file_type 是 *output ，必须确定下面的 name_options 之一。 filename *prompt *screen 对于 file_type 是 *index-in ，必须确定下面的 name_options 之一。 filename *prompt 定义： *filename(s) 文件名关键字，s为可选择字符。*output 输出文件，将写入模拟结果。*index-out 结果索引文件(irf)，将写入 ascii 码模拟结果数据。*main-results-out 主结果文件(mrf)，将写入二进制码模拟结果数据。*rewind-out 可回绕结果文件(rrf)，当使用 *rewind 选项时，将写入重新启动数据。*index-in 结果索引文件，将从中读出模拟结果和重新启动记录，只有在重新启动运行时才需要这个文件。*main-results-in 主结果文件，将从中读出模拟结果和重新启动记录(二进制码)，只有在重新启动运行时才需要这个文件。*rewind-in 可回绕结果文件，将从中读出回绕的重新启动记录(二进制码)，只有在重新启动运行时才需要这个文件。 空字符串，表示使用缺省文件名，如果输入数据文件中没有 *filenames 关键字，此为缺省情况。filename 文件名字符串，最多80个字符，文件名是否被接收取决于使用的操作系统。*prompt 表示如果需要这个文件，将通过屏幕/键盘对用户提示文件名。*screen 表示文件将输出到屏幕(标准输出设备)上。缺省： 如果任何输入/输出文件名未通过*filename关键字确定，则缺省为*prompt。例如：*filename 隐含 *filename *output *prompt。如果使用了*filename file_type 确定一个文件名，只是未给出 name_options ，则缺省为 ，即使用缺省文件名。例如：*filename *output 隐含 为 *filename *output 。 缺省文件名：仅当确定了 name_option 时，才使用缺省文件名，对于任何文件都可以独立地确定缺省，然而某些文件的缺省名依赖于其它文件名的选取，例如输入文件。 一个缺省文件名的形式为xxxxxxxx.yyy，具有个根xxxxxxxx和后缀yyy，这个后缀取决于与其相关的文件选项，例如，out 用于主要输出，irf 用于 sr2 索引，而 mrf 用于 sr2 的主结果文件，根是由另一个文件名导出的，如下： (a) 去掉后缀和 . 分隔符(b) 除了最右边的字母数字组和 - 字符，去掉其它字符(c) 去掉 8 个字符后的任何字符 例如：unix 路径名/cmg/data/test-3.dat具有根 test-3。 下面的表为缺省文件名的根和后缀的汇总： 文件 根的来源 后缀 - - - *output 数据文件名 .out *index-out *output .irf *main-results-out *output .mrf *rewind-out *output .rrf *index-in (无缺省) .irf *main-results-in *index-in .mrf *rewind-in *index-in .rrf 使用这个缺省系统，用户只需要在每次运行中改变数据文件名和 *index-in 文件名就可执行重新启动运行。 例 1 ：数据文件为 cycle.dat，因此 *output 的缺省为 cycle.out。 如果使用 *output 的缺省，那么 *index-out 的缺省是 cycle.irf，*main-results-out 的缺省是 cycle.mrf。 例 2 ：数据文件为 cycle.dat。 在第一次运行时，数据段内含有 *filename *output triall。 那么 *index-out 的缺省是 trial1.irf，*main-results-out 的缺省是 trial1.mrf。 对于第二次运行数据段，修改cycle.dat，并加入 *restart 关键字，然后使用 *filename *output trial2.out，和 *filename *index-in trial1.irf，不改文件名或删除文件。 *index-out 的缺省为trial2.irf，*main-results-out 的缺省为trial2.mrf，*main-results-in 的缺省为trial1.mrf，条件： *filename 在输入数据文件中必须是第一个关键字，否则将对用户提示文件名。 *filename 在输入/输出控制数据段必须是第一个关键字。 所有输出文件皆以unknown方式打开，所以并不对覆盖进行保护。 所有输入文件皆以old方式打开，所以在模拟运行时必须存在。解释： cmg 新的模拟结果文件系统(第二代)也称作 sr2，sr2 文件系统由在一起工作的三个文件组成，它们是结果索引文件(irf)，主结果文件(mrf)和可回绕结果文件(rrf)。对于图形后处理软件“results”，其中两个文件是必须有的，它们是 irf 文件和 mrf 文件。对于重新启动运行也同样需要这些文件，假如使用了 *rewind 选项写重新启动记录，那么对于重新启动运行也需要可回绕结果文件 rrf。请注意： sr2 文件系统与先前的 srf 文件系统不兼容，为了从 sr2 文件系统中生成一个 srf 文件，请使用这个发行版所提供的 sr2 到 srf 的转换程序。 例： * 使用 index-out-file 的缺省文件名 root.irf * 采用不写入输入数据文件的方法 * 让程序提示输入 mrf 文件 *main-results-in *filenames *output root.out *index-out *main-results-out name.mrf *index-in input.irf*dim运行时间定义（任选）_-目的：给维数变量赋值。格式：*dim dim_variable 变量值其中 dim_variable可以是： *memory_model *active_blocks *null_blocks *max_faults *max_aquifers *max_sectors *max_leaselines *maxblks_in_leaselines *max_pvt_regions *max_rock_types *maxentries_in_kr_table *maxentries_in_pvt_table *max_wells *max_layers *max_well_constraints *max_gathering_centers *max_specialvars_in_srfoutput *maxpercent_of_fullyimplicitblocks *solver_dimensioning *solver_dimensioning 必须在下列选项读入之前读入： *max_orthogonalizations *max_interblock_connections *max_offdiagonal_entries *max_blockentries_in_lu_factors *maxsize_solver_vector *max_solver_diag_entries *max_nonzero_l_or_u_entries定义：*dim 关键字*memory_model设置维数的参数，使imex以一定的内存空间运行，可设为8，16，26，32，64或128数字大致反映出内存的大小（以mb为单位）如运行一个128mb的模型,可输入*dim *memory_model 128。若出现维数方面的问题，将memory_model设为128，或许问题就会解决。*active_blocks 有效网格数。有效网格包括双孔网格，混合网格和加密网格，不包括死节点（孔隙度为零的网格）。例如 *dim *active blocks 2000*null_blocks 死节点数。缺省情况是没有死结点，即所有的网格都当作有效网格。例如 *dim null_blocks 500*max_faults 最多断层数。例如 *dim max_faults缺省值为15。*max_aquifers最多侵入区数 例如*dim max_aquifers 15。缺省值为15。*max_sectors最多分区数。例如*dim max_sectors 50。缺省值为50。*max_leaselines最多租赁线数。例如*dim max_leaselines 12。缺省值为15。*maxblks_in_leaselines 一条租赁线上最多网格数。例如 *dim maxblks_in_leaselines 5000。缺省值为150*max_pvt_regions 允许的最多pvt区域 例如*dim max_pvt_regions 50 缺省值为15。*max_rock_types 允许的最多岩石类型 例如*dim max_rock_types 100 缺省值为30。*maxentries_in_kr_table 岩石-流体数据表中最多输入项 例如*dim maxentrie_in_kr_table 51。缺省值为51。*maxentries_in_pvt_table 每个pvt表中最多输入项 例如*maxentries_in_pvt_table 51。缺省值为51。*max_wells 最多井数（包括关闭的井） 例如*dim max_wells 600 。 缺省值取决于数据文件中定义的井的信息。*max_layers 完井段包含的最多层数 例如 *dim max_layers 5。 缺省值为5。*max_well_constraints 最多井的操作或监测限制条件。 例如 *dim max_well_constraints 5。 缺省值为5。*max_gathering_centers 最多集输中心数。 例如 *dim max_gathering_centers 4。 缺省值为15。*masx_specialvars_in_srfoutput srf 输出中允许的最多*special 变量数。例如，*dim *masx_specialvars_in_srfoutput 500。缺省值为100。*maxpercent_of_fullyimlitblocks 全隐式网格的最大百分数。取值在10到100之间。若全隐式选100。例如,*dim maxpercent_of_fullyimlitblocks 100。对于油水模型和黑油模型缺省值为100，对拟混相驱模型和聚合物驱模型缺省值为42。*solver_dimensioning 例如，*dim solver_dimensioning on解法数组按常规定义维数；*dim solver_dimensioning off解法数组不按常规定义维数。缺省为off。*max_orthogonalizations线性解法中正交化次数。减少正交化次数能节约存储空间，但效果不好。只有当solver_dimensioning设为on时，选项才有效。例如，*dim max_orthogonalizations 5.缺省值为11。*max_interblock_connections 最多的网格间连接。网格间的连接包括网格间两两相连，以及由于井而产生的连接。只有solver_dimensioning置为on，该选项才有效。例如，*dim maxlnter_block_connection 5。缺省值是根据基础网格数ni、nj和nk，死节点数以及模拟的问题类型来估算的。*max_offdiagonal_entries 非对角jacobi矩阵中最多输入项数只有当solver_dimensidoning置为on,该选项才有效。例如 *dim maxdffdiagonalenteries 5。缺省值是根据基础网格数，ni、nj和nk，死节点数以及模拟的问题类型来估算的。*max_blockentries_in_lu_factors lu分解输入的最多网格数。只有当solver_dimensioniong置为on，该项才有效。例如，*dim max_blockentries_in_lu_faotors 5。缺省值是根据基础网格数，ni、nj和nk，死节点数以及模拟的问题类型来估算的。*max_solver_vector 解向量的最大维数。只有当solver_dimensioning置为dn，该选项才是有效的。例如*dim max_solver_vector 5。缺省值是根据基础网格数，ni、nj和nk，死节点数以及模拟的问题类型来估算的。*max_solver_diag_entries 对角线上的最多输入项数。只有当solver_dimensioning置为on，该选项才有效。例如，*dim max_solver_diag_entries 5缺省值是根据基础网格数，ni、nj和nk，死节点数以及模拟的问题类型来估算的。*max_nonzero_l_or_u_entries lu矩阵中非零项的最多项数。只有当solver_dimensioning置为on，该选项才有效。例如 *dim max_nonzero_l_or_u_entries 5。缺省值是根据基础网格数，ni、nj和nk，死节点数以及模拟的问题类型来估算的。*title1项目主标题 (可选择)_目的： *title1 用于项目标识。格式： *title1 string 定义：string 包含在单引号之间的任何字母数字字符串。(最多40个字符)缺省： 可选择关键字，缺省为空格。条件： 此项关键字必须出现在输入文件开始的输入/输出控制数据段中、说明： 此项关键字说明一个用于标识项目的字母数字字符串，它将出现在输出文件和结果索引文件中。 任何含有空格和逗号的字符串必须使用单引号括住。 例如： *title1 dual porosity/dual permeability run no.1 *title2项目第二标题 (可选择)_目的： *titel2 提供方案标识的第二行，作为 *title1 的补充。格式： *title2 string 定义：string 任何用单引号括起的字母数字字符串，最多80个字符。 缺省： 此项为可选择关键字，缺省为空格。条件： 此项关键字必须出现在输入文件开始的输入/输出控制数据段中、说明： 此项关键字说明一个用于标识项目的字母数字字符串，它将出现在输出文件和结果索引文件中。 任何含有空格和逗号的字符串必须使用单引号括住。 例如： *title1 dual porosity/dual permeability run no.1 *title2 run by f.r. sangiovanni, december 16, 1988. ibm *title3项目第三标题 (可选择)_目的： *titel3 提供方案标识的第三行，作为*title1 和 *title2 的补充。格式： *title2 string 定义：string 任何用单引号括起的字母数字字符串，最多80个字符。 缺省： 此项为可选择关键字，缺省为空格。条件： 此项关键字必须出现在输入文件开始的输入/输出控制数据段中、说明： 此项关键字说明一个用于标识项目的字母数字字符串，它将出现在输出文件和结果索引文件中。 任何含有空格和逗号的字符串必须使用单引号括住。 例如： *title1 dual porosity/dual permeability run no.1 *title2 run by f.r. sangiovanni, december 16, 1988. ibm *title3 imex 4.0; 4200 grid blocks; variable thickness*caseid方案标识 (可选择)_目的： *caseid 用于标识方案运行。格式： *caseid string定义：string 任何用单引号括起的字母数字字符串，最多8个字符。 缺省： 此项为可选择关键字，缺省为空格。条件： 此项关键字必须出现在输入文件开始的输入/输出控制数据段中、说明： 此项关键字说明一个用于标识方案运行的字母数字字符串，它将出现在输出文件和结果索引文件中。 任何含有空格和逗号的字符串必须使用单引号括住。 例如： *caseid case 21c*checkonly只进行检查 (可选择)_目的： *checkonly 表示只检查井数据的格式错误，不做模拟运行。格式： *checkonly缺省： 此项为可选择关键字，如不存在则不检查井数据。条件： 这一关键字必须位于输入/输出控制段内。*interupt 陷阱unix #2中断（任选）_目的：表示当用户同时按下ctr_c键或键入”kill_2进程号（pid）”来终止一个运行时，可以转而采取下列行动：1. 继续进行模拟。2. 完成当前时间步，写重启记录并结束模拟。3. 刷新缓冲区，关闭文件并立即结束模拟。4. 写一个临时的重启记录，并继续运行。5. 仅对pc机而言，立即无条件结束模拟（control_break键）。上述1-3条，可以防止.irf和.mrf文件在退出模拟时被破坏，同时可以有选择地在模拟结束前写重启记录。第4条可以写一个临时的重启记录。格式：*interupt (*restart_stop) (default) (*wrst_cont) (*stop) (*stop_cold) (*interactive)定义：*restart_stop 表示完成当前时间步，所有定义的输出数据写到输出文件和sr2文件，同时重启记录也写到sr2文件。*wrst_cont 确定写一个临时的重启记录，也就时说，所有定义的输出数据写到输出文件和sr2文件，重启记录也写到sr2文件，模拟继续进行。*stop 确定模拟立即终止。当前时间步不必完成，但在模拟终止前，所有的缓冲区被刷新，文件关闭以防破坏。*stop_cold 模拟立即无条件结束运行。contyol_c和交互模式下的control_break都能够产生这样的结果。*interactive 提示用户输入要采取的措施。缺省：该关键字为任选关键字。缺省时表示要提示用户输入上面三种措施中的一种。若只输入*interupt而没有输入措施，缺省措施为*restart_stop。条件：该关键字必须位于输入/输出控制数据段内。*sr2ascii sr2文件ascii码形式文件（任选）_目的：产生一个主要结果文件（.mrf）的一个ascii码拷贝。该备份文件扩展名为.asc。格式：*sr2ascii缺省：该关键字是任选关键字。无缺省值条件：该关键字是必须出现在输入/输出控制数据段内。*inunit输入数据单位 (可选择)_ 目的： *inunit 用于确定输入数据单位。格式： *inunit (*si) (*field) (*lab) (*modsi) 定义：*si 这个选项确定输入数据使用国际标准单位。 *fifld 这个选项确定输入数据使用矿场单位(英制)。 *lab 这个选项确定输入数据使用实验室单位。 *modsi 这个选项确定输入数据使用修改的国际标准单位。缺省： 这一关键字为可选择项，缺省为 *si。条件： 这一关键字必须位于输入数据文件开始的输入/输出控制段内。说明： 数据可使用下面单位表中四种选项中的一种输入。 例如： *inunit *field 单位 表 _ | | | | | | | | *si | *field | *lab | *modsi | |_|_|_|_|_| | | | | | | | 时间 | days | days | mins | days | |_|_|_|_|_| | | | | | | | 液体体积es | m3 | bbl | cm3 | m3 | |_|_|_|_|_| | | | | | | | 气体体积 | m3 | ft3 | cm3 | m3 | |_|_|_|_|_| | | | | | | | 压力 | kpa | psi | kpa | kg/cm2 | |_|_|_|_|_| | | | | | | | 绝对渗透率 | md | md | md | md | |_|_|_|_|_| | | | | | | | 粘度 | mpa-s | cp | mpa-s | mpa-s | |_|_|_|_|_| | | | | | | | 质量密度 | kg/m3 | lbm/ft3 | g/cm3 | kg/m3 | |_|_|_|_|_| | | | | | | | 长度，距离 | m | ft | cm | m | |_|_|_|_|_| | | | | | | | 溶解气油比 |