BOAST黑油用户手册范本.docx

《BOAST黑油用户手册范本.docx》由会员分享，可在线阅读，更多相关《BOAST黑油用户手册范本.docx（22页珍藏版）》请在三一文库上搜索。

1、1程序概述软件介绍BOAST程序模拟以三维情况的等温达西流，它假设油藏中的流体是具有恒定成分，且物性只由压力决定的三相流体油、气、水。这种对油藏流体的近似适用于世界上的大部分油、气藏。因此， BOAST能被广泛应用。例如： BOAST能模拟由流体的膨胀、驱替、重力驱和毛管力吸入机理进行的油和（或）气的开采。能用 BOAST处理的典型油田生产问题有：一次开采的研究，注水和（或）气保持压力的开采。注水、驱替等二次开采的评价。从技术上看： BOAST是一个隐式压力，显式饱和度（ IMPES）有限差分数值模拟软件，在井动态方面， BOAST是井模型允许规定井的产率或压力， 模拟过程中， 用户可以任意增

2、加或减少井数， 在BOAST里有压力初始化的算法，泡点轨迹线，时间步长自动控制和检测解的稳定性的平衡法，此外，许多由用户控制的输出选择项也是有效的。数据输入要求：数据输入段分为两部分： 初始化数据段和循环数据段。初始化数据包括： 油藏模型网格的大小及形状、孔隙度、渗透率在油藏中的分布，流体的 PVT数据，岩石的相对渗透率和毛管压力数据， 油藏的原始压力及饱和度分布， 解法说明和各种运行控制参数。循环数据包括： 模型中井的位置及其初始规定， 在时间上继续模拟下去所需要的时间步长控制信息，单井产率和（或）压力动态表，完井和井下作业随时间的变化以及对模拟程序输出资料和类型和输出次数的控制。本段将只讨

3、论模型所要求的输入数据的类型的格式。注意：下面有关数据输入格式 （特别是对输入数组的规定） 的注释适用于所有数据输入。1在本手册中的输入数据行称作“卡”，其含意是：在卡片输入计算机上表示计算机输入卡， 而在输入数据的显示器或行式打印机上表示单独的一行， 每张卡片的输入格式在紧随卡片描述的圆括号中说明。2标题卡先于每个主要的和许多次要的数据输入段读入。除非在文件中有别的说明，必须在每种情况下都读入卡片， 为了使输入数据文件易读易编辑， 将这些卡片设计成易看的图型。3在很多情况下，读入控制码来规定将要输入的数据的类型和数值个数，为便输入文件可使用选择项具有灵活性，必须读入控制码。4除标题卡外， 所

4、有数值都按固定格式输入， 考虑到要让输入文件易读易编辑，大多数的输入格式为 I5 ，F10.0 或 F8.0 。注意：在程序中，空格将按0 读入5若某一特定参数要求输入所有网格的值，则必须按如下顺序输入：首先读第一层（ K1）对每一层的数据输入是按行进行的，先读第一行（ J 1），对每一行以是以第一列 （ I 1）开始的，第一行读入 I 1 到 II 列的数值，读完第一行，再读第二行（ J2），第二行仍从 I 1 到 II ，以此类推，直至读完 JJ 行，再对第二层重复上述过程，直至读完 KK层。6记住 BOAST用的坐标原点在左边上角，因此每一层的网格分布如下：I1I 2J=1J=2.Z 方

5、向的值向下为增。对运行黑油模型的提示：在解决油藏工程和生产问题方面，BOAST是一个复杂的工程工具，为了更灵活地处理在油藏管理和生产过程中遇到的大量问题就更需要象 BOAST这样的一个通用的且有些复杂的（因为程序中包含很多选择项）油藏模拟软件。也许尽快熟悉如何运行 BOAST的最佳方法是先浏览一下输入数据和熟悉输入文件的一般格式。BOAST具有对时间步长自动控制的特点，建议在大多数运行中使用这种自动控制，当地层条件随时间变化不大时， 这种特点可通过增加步长使程序最有效地工作，而当地下条件变化较大时（如泡点压力从未饱和状态变到饱和状态），较小的时间步长将使程序具有较大的稳定性， 建议将时间步长自

6、动控制的最小值设为 0.1 天。根据实际应用的要求还可选用更小的时间步长，例如，模拟用实验岩心注水时要求的最小时间步长为 0.01 天或更小，本手册提供的所有缺省都刻意于油田规模的实际应用， 这种小的时间步长将只用于用户指定的一个时间步长里发生的最大压力和饱和度变化。2初始化数据初始化数据卡描述的是油藏模型网格的大小和形状， 以及孔隙度、 渗透率、相对渗透率、毛管压力数据、流体的 PVT值。初始压力和饱和度在油藏中的分布，同时还描述了使用的解法及各种控制、诊断参数。这些数据卡只在模拟开始时读入一次，它们必须按如下的数据段中的顺序读入。网格尺寸和形状网格尺寸1标题卡（ 40A2）2重启动控制码2

7、I5IW 1 建立重启动文件0 不建立重启动文件IR 1 需重启动运行2 不需重启动3标题卡（ 40A2）4模型的网格块数（ 3I5 ）II X 方向的网格块数JJ Y 方向的网格块数KK Z 方向的网格块数5标题卡（ 40A2）6规定输入类型的控制码（3I5 ）KDX 输入 X 方向的网格尺寸的控制码KDY 输入 Y 方向的网格尺寸的控制码KDZ 输入 Z 方向的网格尺寸的控制码控制码含意KDX 1 所有网格块的 X 方向网格尺寸相同（只读一个值）KDX 0读入第一层（ K1）第一行（ J 1）的每个网格块的X 方向网格尺寸，并将同样的值赋给第一层其它行和其它层的所有行 X 方向的对应网格（

8、注：必须读II个值）KDX 1 读入第一层（ K1）的每个网格块的X 方向网格尺寸，并将同样的值赋给其它层的对应网格。（注：必须读 II JJ 个值）KDY 1 所有网格块的 Y 方向网格尺寸相同（只读一个值）KDY= 0 读入第一层（ K 1）第一列（ I 1）的每个网格 Y 方向尺寸，并将相同的值赋给第一层的其余列及其它所有层的各列的对应网格块。（注：必须读JJ 个值）KDY 1 读第一层（ K 1）的所有网格块的Y 方向尺寸，并将同样的值赋给其余各层对应网格。（注：必须读II JJ 个值）KDZ 1 所有网格块 Z 向尺寸（块厚）相同。（只读一个值）KDZ 0每一层读入一个不变的厚度值，

9、不同层可以有不同的厚度。（注：必须读KK个值）KDZ 1 读模型的每个网格Z 方向的尺寸（厚度）（注：必须读II JJKK个值）7网格 X 方向的尺寸（ DX）（ 10F8。0）如果 KDX 1，只读一个值如果 KDX 0，读 II个值（第列一个值）如果 KDX 1，读 II JJ 个值（第一层的每个网格块有一个值）8网格 Y 向的尺寸（ DY）（ 10F8。 0）如果 KDY 1，只读一个定值如果 KDY 0，读 JJ 个值（第行一个值）如果 KDY 1，读 II JJ 个值（第一层的每个网格有一个值）9网格 Z 向的尺寸（ DZ）（ 10F8。 0）如果 KDZ 1，只读一个值如果 KDZ

10、 0，读 KK个值（第一层一个值）如果 KDZ 1，读 II JJ KK个值（整个模型的每个网格有一个值）网格尺寸的修改1标题卡（ 40A2）2要修改尺寸的网格数和打印控制码（4I5 ）NUMDX 需要修改 X 方向尺寸（ DX）的网格数NUMDY 需要修改 Y 方向尺寸（ DX）的网格数NUMDZ 需要修改 Z 方向尺寸（ DX）的网格数IDCODE 打印控制码IDCODE 0 不打印修改的网格尺寸IDCODE 1 打印修改的网格尺寸3网格 X 向尺寸（ DX）的修改（ 3I5 ， F10。0）如果 NUMDX0 就不读这些数据卡I 被改网格的 X 坐标J 被改网格的 Y 坐标K 被改网格的

11、 Z 坐标DX 网格块（ I ，J，K）X 方向网格尺寸（ DX）的新值。注：必须读 NUMDX卡4网格 Y 向尺寸（ DX）的修改（ 3I5 ， F10。0）如果 NUMDY0 就不读这些数据卡I 被改网格的 X 坐标J 被改网格的 Y 坐标K 被改网格的 Z 坐标DY 网格块（ I ，J，K）Y 方向网格尺寸（ DY）的新值。注：必须读 NUMDY卡5网格 Z 向尺寸（ DZ）的修改（ 3I5 ， F10。0）如果 NUMDZ0 就不读这些数据卡I 被改网格的 X 坐标J 被改网格的 Y 坐标K 被改网格的 Z 坐标DZ 网格块（ I ，J，K）Z 方向网格尺寸（ DZ）的新值。注：必须读

12、 NUMDZ卡第一层网格块顶部的深度记住，在 BOAST所使用的坐标系下， Z 方向的值是向下递增的，因此，读入的深度应为低于用户选定的基准面的值（高于此基准面为负）1标题卡（ 40A2）2输入深度值的控制码（I5 ）KEL 输入控制码KEL 0 意指对第一层的所有网格块的顶部深度只读入一个定值（即第一层为水平面）KEL 1 意指对第一层的每个网格块各读入一个深度值，共计读II JJ 个值。3深度值（ 8F10。0）ELEV 网格块顶部注意：如果 KEL0，则只读一个值。如果KEL1，则读入 II JJ个值，对第一行（ J 1读 II 个值，然后对第二行（ J2）又读 II 个值，如此读至第

13、JJ 行。除第一层外，各层网格块顶部深度用其前一层网格顶部深度加上前一层网格厚度得到：即： TOP（I ，J，K1）TOP（I ，J，K） DZ（I ，J，K）孔隙度和渗透率的分布孔隙度和渗透率1 标题卡（ 40A2）2 输入类型控制码（ 4I5 ）KPH 控制孔隙度数据输入的控制码KKX 控制 X 方向渗透率输入的控制码KKY 控制 Y 方向渗透率输入的控制码KKZ 控制 Z 方向渗透率输入的控制码控制码含意 1 读入一个定值并赋给模型的所有网格块（只读一个值）0 对模型 KK层的每一层读入一个定值，不同层可以有不同的定值。（必须读 KK个值） 1 对模型的每一个网格块都读入一个值（必须读

14、II JJ KK个值）3孔隙度值（ 10F8。 0）孔隙度不是按百分数，而是按小数形式读入的。如果 KPH 1，则只读一个定值如果 KPH 0 ，则读 KK个值（一层一个值）如果 KPH 1，则读 II JJ KK个值（模型之第一个网格块都有一个值）4X 方向渗透率值（ KX）（ 10F8。0）渗透率值按毫达西（ md）读入如果 KKX 1，则读一个定值如果 KKX 0 ，则读 KK个值（一层一个值）如果 KKX 1，则读 II JJ KK个值（模型的每个网格块都有一个值）5Y 方向渗透率（ KY）（ 10F8。0）Y 方向渗透率值（ KY）（ 10F8。0）渗透率值按毫达西（ md）读入如果

15、 KKY 1，则读一个定值如果 KKY 0 ，则读 KK个值（一层一个值）如果 KKY 1，则读 II JJ KK个值（模型的每个网格块都有一个值）6Z 方向渗透率值（ KZ）（ 10F8。0）渗透率值按毫达西（ md）读入如果 KKZ 1，则读一个定值如果 KKZ 0 ，则读 KK个值（一层一个值）如果 KKZ 1，则读 II JJ KK个值（模型的每个网格块都有一个值）对孔隙度、渗透率分布的修改1 标题卡（ 40A2）2 要进行孔隙度和（或）渗透率值修改的网格数和打印控制码（5I5 ）NUMP要修改孔隙度的网格块数NUMKX要修改 X 方向渗透率值（ KX）的网格块数NUMKY 要修改 Y

16、 方向渗透率值（ KY）的网格块数NUMKZ 要修改 Z 方向渗透率值（ KZ）的网格块数IPCODE 打印控制码IPCODE 不打印修改了的孔隙度和渗透率分布IPCODE 打印修改了的孔隙度和渗透率分布3孔隙度的修改（ 3I5 ， F10。0）如果 NUMP0，则不读这张数据卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K 被改网格块的 Z 坐标PHI 网格块（ I ， J， K）的新的孔隙度值注：孔隙度必须以小数形式输入，必须读NUMP卡4X 方向渗透率（ KX）的修改（ 3I5 ，F10。 0）如果 NUMKX0，则不读这张卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K

17、 被改网格块的 Z 坐标KX 网格块（ I ，J，K）的新的 X 向渗透率值。注：必须读 NUMKX卡5Y 方向渗透率（ KY）的修改（ 3I5 ，F10。 0）如果 NUMKY0，则不读这张卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K 被改网格块的 Z 坐标KY 网格块（ I ，J，K）的新的 Y 向渗透率值。注：必须读 NUMKX卡4Z 方向渗透率（ KZ）的修改（ 3I5 ，F10。 0）如果 NUMKZ0，则不读这张卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K 被改网格块的 Z 坐标KZ 网格块（ I ，J，K）的新的 Z 向渗透率值。注：必须读 NUMKX卡传

18、导系数的修改注：记住在修改传导系数中对方向的规定是十分重要的。例如，对网格块（ I ， J， K）： TX（ I ， J， K）指的是穿越网格块 I 1 与 I 的边界的流动。 TY（I ，J，K）指的是穿越网格块 J1 与 J 的边界的流动以及 TZ（I ，J，K）指的是穿越网格块 K1 与 K 的边界的流动。1标题卡（ 40A2）2要修改传导的网格块数，加之打印控制码（4I5 ）NUMTX 要修改 X 方向传导系数（ TX）的网格块数NUMTY 要修改 Y 方向传导系数（ TY）的网格块数NUMTZ 要修改 Z 方向传导系数（ TZ）的网格块数ITCODE 打印控制码ITCODE 0，不打

19、印修改了的传导系数分布ITCODE 1，印修改了的传导系数分布3X 方向传导系数（ TX）的修改（ 3I5 ， F10。0）如果 NUMTX0，则不读这张卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K 被改网格块的 Z 坐标TX 网格块（ I ，J，K）X 方向传导系数（ TX）的新值。注：必须读 NUMTX卡4Y 方向传导系数（ TY）的修改（ 3I5 ， F10。0）如果 NUMTY0，则不读这张卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K 被改网格块的 Z 坐标TY 网格块（ I ，J，K）X 方向传导系数（ TY）的新值。注：必须读 NUMTY卡5Z 方向传导系数

20、（ TZ）的修改（ 3I5 ， F10。0）如果 NUMTZ0，则不读这张卡I 被改网格块的 X 坐标J 被改网格块的 Y 坐标K 被改网格块的 Z 坐标TZ 网格块（ I ，J，K）X 方向传导系数（ TZ）的新值。注：必须读 NUMTZ卡相对渗透率和毛管力表1标题卡（ 40A2）2读相对渗透率和毛管力表（6F10。0）SAT1 KRO1 KRW1 KRG1 PCOW1 PCGO1 SATN KRON KRWN KRGN PCOWN PCGONSAT 相饱和度值。注意： SATI 应为 0.01 ，而 SATN必须是1.10 。每个饱和度均按小数读入KRO 油相相对渗透率（小数）KRW 水相

21、相对渗透率（小数）KRG 气相相对渗透率（小数）PCOW 油 / 水毛管压力（ at ）PCGO 气 / 油毛管压力（ at ）SAT指的是每一特定相的饱和度，即：在SAT0.2 的数据行上：KRO指的是当油的饱和度为20时，油的相对渗透率。KRW指的是当水的饱和度为20时，水的相对渗透率KRG指的是当气的饱和度为20时，气的相对渗透率PCOW指的是当水的饱和度为20时的油 / 水毛管压力PCGO指的是当气的饱和度为20时的气 / 油毛管压力流体的 PVT数据表1 标题卡（ 40A2）2 泡点压力，未饱和油的性质，PVT表中的最大压力（ 5F10。 0）PBO 油的原始泡点压力VSLOPE 未

22、饱和油（ PPBO时的油）的粘度与压力关系曲线的斜率。 BSOLOPE 未饱和油的体积系数与压力关系曲线的斜率。注：与未饱和没的压缩系数不同，BSLOPE应为负值RSLOPE 未饱和油的溶解油气比与压力曲线的斜率（通常情况下该值为 0）PMAX 所有 PVT表的最大压力记录IREPRS3标题卡（ 40A2）4油的 PVT表（ 8F10。0）P1 MUO1 BO1 RSO1 PMAX MUOPMAX BOPMAX RSOPMAXP 压力（注：正如上表所示，最后一个压力记录必须是PMAX）MUO 饱和油粘度BO 饱和油地层体积系数RSO 饱和油溶解油气比重要问题：在整个压力范围，油的物性必须作为饱

23、和油的数据输入，一般来说，实验室的饱和油数据必须进行外推，直到高于所测得的泡点压力， 以覆盖预计将在模拟运行中出现的最大压力范围。要求输入饱和油的数据。是因为在 BOAST 中使用了泡点轨迹线。 请注意，只有当当地油藏压力增至高于原始泡点压力并且自由气进入该油藏时，才使用高于原始泡点压力的饱和油数据。 （例如：向油藏注气保持压力）5 标题卡（ 40A2）6 水的 PVT表（ 8F10。 0）P1 MUW1 BW1 RSW1 PMAX MUWPMAX BWPMAX RSWPMAXP 压力（注：正如上表所示，最后一个压力记录必须是PMAX）MUW 水的粘度BW 水的地层体积系数RSW 水的溶解油气

24、比在黑油模拟中， 通常假定气在水中的溶解度可以忽略不计， 在这种情况下， 对所有压力都设置 RSW 0，在 BOAST中有水的这个 PVT表是为处理如下情况用的 从地层含水层采气的情况以及很多被认为气在水中的溶解度对问题的解有较大影响的情况。7 标题卡（ 40A2）8 气的 PVT表及岩石的压缩系数（ 8F10。0）P1 MUG1 BG1 CR1 PMAX MUGPMAX BGPMAX CRPMAXP 压力（注：正如上表所示，最后一个压力记录必须是PMAX）MUG 气的粘度BG 气的地层体积系数CR 取决于压力岩石压缩系数9 标题卡（ 40A2）10 储罐流体密度（ 8F10。 0）RHOSC

25、O 储罐油的密度RHOSCW 储罐水的密度RHOSCG 标准状况下的气体的密度压力和饱和度的初始化BOAST含有两个有关压力和饱和度初始化的选择项， 初始压力分布既或由考虑了油、气和油水接触位置及接触位置处的压力的平衡条件的程序计算出（选择项1），也可象在非平衡初始化情况下那样一个一个网格块地读入（选择项 2），饱和度值（ SO，SW，SG）既可对整个网格块作为定值读入（选择项 1），也可在一个一个网格块的基础上读入整个 SO和 SW的分布，并由程序计算出 SG的分布。SG1SO SW（选择项 2）1 标题卡（ 40A2）2 压力及饱和度分布的控制码（2I5 ）KPI 压力初始化控制码KSI

26、饱和度初始化控制控制码含义KPI 0 使用平衡压力初始化，要求输入油水和油气接触处的压力以及各接触处的深度KPI 1 使用非平衡压力初始化，必须在一个一个的网格块基础上读入每个网格块的压力（注：必须读 II JJ KK个值）KSI=0 整个模型网格的原始油、气、水饱和度为定值注：必须读三个值： SOI，SWI， SGI）KSI=1 必须在一个一个的网格的基础上对每个网格块读入油、水饱和度，每个网格块的气体饱和度由程序计算出（注：必须读 2*II*JJ*KK个值）3压力初始化数据（ 10F8.0 ）如果 KPI=0（平衡初始化）PWOC 油水接触处的压力PGOC 油气接触处的压力WOC 油水接触

27、处的深度（低于基准面的英尺数）GOC 油气接触处的深度（低于基准面的英尺数）如果 KPI=1（非平衡初始化），则对模型的每个网格块都读入一个压力值（注：必须读 II*JJ*KK个压力）4饱和度初始化数据（ 10F8.0 ）如果 KSI=0（恒定饱和度情况），则SOI 把原始含油饱和度（小数）赋给模型的所有网格块SWI 把原始含水饱和度（小数）赋给模型的所有网格块SGI 把原始含气饱和度（小数）赋给模型的所有网格块如果 KSI=1（非平衡初始化）SO数组：对模型的每一网格块都读入原始含油饱和度值（注：必须读 II*JJ*KK个含油饱和度值）SW数组：对模型的每一个网格块都读入原始含水饱和度值（注

28、：必须读 II*JJ*KK个含水饱和度值）如果 IR=1 时， 3 与 4 二项（ P，S）不赋值调试和诊断控制码在程序调试方面， 本程序提供了几个控制诊断输出的控制码。 这些控制码通常被设置为 0，这些控制码将不给出有关输入数据调试的信息，活化这些控制码的任何一个都将产生极大量的输出。1标题卡（ 40A2）2控制诊断输出的控制码（5I5 ）KSN1 LSOR参数调试输出控制KSM1 解矩阵调试输出控制KCO1 压缩系数和体积系数调试输出控制KTR 传导系数调试输出控制KCOF 密度和饱和度调试输出控制控制码：0 不打印诊断输出1 打印诊断输出运行控制参数1标题卡（ 40A2）2运行控制参数（

29、 I5 ，710.0 ）NMAX 运行结束前允许的最大时间步数FACT1 在自动控制时间步长的情况下，增大时间步长的因子（对恒定时间步长，设置FACT1=1.0， FACT1的一个常用值为1.25 ）FACT2 在对时间步长自动控制的情况下，减小时间步长的因子（对固定时间步长设置FACT2=1.0）TMAX 在程序运行中模拟的最大实际时间，天（当模拟时间超过了 TMAX时，运行将结束）WORMAX 油田限定的最大油水比（当模拟的整个生产WOR超过了WORMAX时，运行结束）GORMAX 油田的最大油气比限定值（当模拟的总的GOR超过 GORMAX时，运行结束）PAMIN 油田限定的最低平均压力

30、值（当模拟的平均油藏压力低于PAMIN时，运行结束）PAMAX 油田限定的最高平均压力值（当模拟的平均压力超过了PAMAX时，运行结束）解法说明1标题卡（ 40A2）2解法控制参数（ 2I5 ， 6F10.0 ）KSOL解法控制码KSOL 解法1 直接解法条带算法（应用于一维问题）2 LSOR（应用于大的二维和三维问题）3 直接解法算法（应用于二维问题）MITR 每一时间步长 LSOR的最大迭代次数。MITR的一个典型值是250。OMEGA LSOR的初始加速参数。OMEGA的初值范围是：。 OMEGA的一个典型初始值为 1.7 ，随着解题过程的进行，程序将优化OMEGA值。TOL KSOR迭

31、代收敛性允许的最大压力变化。TOL的一个典型值为0.1。TOL1 决定何时改变 OMEGA的参数。 TOL1的一个典型值是0.001 。（注：如果 TOL1=0，则上面赋给 OMEGA的初值在整个程序运行中保持不变）DSMAX 一个时间步长所允许的最大饱和度变化（小数），如果在一个时间步长下，任何网格块的任一相饱和度变化超过了DSMAX，则时间步长将按FACT2缩小， DSMAX的一个典型值是 0.05 。DPMAX 一个时间步长所允许的最大压力变化。如果在一个时间步长里，任何一个网格块的压力变化超过了DPMAX，则时间步长将按 FACT2缩小， DPMAX的一个典型值是5at时间步长和输出控

32、制注：在模拟开始时（即当 TIME=0时），和在按（ ICHANG*DT）定义的每个时间间隔的开始必须读以下的两个时间和输出控制卡。1时间步长和输出控制码（9I5 ）IWLCNG告诉程序在这个时间步长里是否要读井资料的控制码控制码含意IWLCNG=0在这个时间步长里不读井资料IWLCNG=1在这个时间步长里读井资料注：如果 IWLCNG=1，则必须读井资料卡，在下一时间步长时将使用新的井资料ICHANG 将使用输出控制资料和时间步长控制资料的时间步长数IWLREP 控制打印井报告的输出控制码ISUMRY 控制打印时间步长综合报告的输出控制码IPMAP 控制打印网格块压力图的输出控制码ISOMA

33、P 控制打印网格块含油饱和度的输出控制码ISWMAP 控制打印网格块含水饱和度的输出控制码ISGMAP 控制打印网格块含气饱和度的输出控制码IPBMAP 控制打印网格块饱和压力的输出控制码（通常设置IPBMAP=0）（ IPBMAP=0不打印报告 / 数组，IPBMAP=1在这个时间间隔的每一时间步长都打印报告 / 数组）2时间步长控制资料（ 3F10.0 ）DT 初始时间步长大小，单位：天。（如果饱和度或压力超过限定值，则自动时间步长控制将减小初始时间步长）注：使用卡 1 的控制码和卡 2 的资料的总的时间是（ ICHANG*DT），单位：天DTMIN 在此期间使用的最小时间步长，单位：天（

34、该值通常设为0.1 天）DTMAX 在此期间使用的最大时间步长，单位：天（ DTMAX不要超过 30 天）井资料注：只有当前面一对时间步长和输出控制卡中的IWLCNG=1时，才读下面的套卡1标题卡（ 40A2）2模型中的总井数（ 3I5 ）NVQN 要读井资料的总井数NQO 油井数NQW 水井数注：在模拟过程的任何时间都可增加井或再完井，但是，一旦一口井被给定了，则在每次读这张卡时都必须描述它，即使这中井目前关闭了也是这样。3井资料（ A5，5I5 ，4F10.0 ）WELLID 含 5 个字符的井名I 井所在网格块的 X 坐标J 井所在网格块的 Y 坐标PERF1 这中井最上面的完井层的层位

35、NLAYER 连续完井的总层数，从PERF1开始并包含 PERF1，因此，在一个五层模型网格里，如果一口井是在第2、3和5层完井的，则设置PERF1=2和 NLAYER=4，注意在这种情况下，即使第 4 层没有射孔也必须包括在NLAYER中。在下面的卡 7 中，第 4 层的 PID 值将被设置为 0。KIP 是既说明井的类型，又说明井的生产（注入）动态是取决于给定的产率还是取决于井底流压的控制码。该控制码还说明程序将进行的是显示压力计算还是隐式压力计算，对大多数情况，建议使用显示压力计算。控制码含意KIP=3 气井 - 规定的产率KIP=2 水井 - 规定的注水率KIP=1 生产井（油井） -

36、 规定的产率KIP=-1 生产井（油井） -PI及 FBHP控制码（显式压力计算）KIP=-2 水井 -PI及 FBHP控制码（显式压力计算）KIP=-3 气井 -PI及 FBHP控制码（显式压力计算）KIP=-11 生产井（油井） -PI及 FBHP控制码（隐式压力计算）KIP=-12 水井 -PI及 FBHP控制码（隐式压力计算）KIP=-13 气井 -PI及 FBHP控制码（隐式压力计算）注：气井（ KIP=3，-3 ，或 -13 ）将只生产或注入单一气相；水井（ KIP=2，-2 ，或 -12 ）将只生产或注入单一水相；生产井（油井）（ KIP=1，-1 ，-11 ）将生产三相流体，它

37、们的比例由油藏条件和井的约束条件决定。注：这张卡的下面 4 个值是这口井给定的油、水、气的产率及总的流体产率 .4 个值中的任一个最大值可能是非零的数，如果 KIP（上面的）的值为负，则这四个值都按 0 读入。在读产率时的符号规定如下；负产率表示注入流体，正产率表示产出流体。 QO 产油率（只有当 KIP=1 且 QT为 0 时为非零值）QW 产水率（只有当 KIP=2 时为非零值）QG 产气率（只有当 KIP=3 时为非零值）QT 总的流体产率（只有当KIP=1 且 QO为 0 时为非零值）由 QT给出的总的流体产率是井的油、气、水产量之和，或者是油藏总的亏空体积。4井资料卡（ 2F10.0 ）注：对每一个 WELLID都必须读这些卡的 NLAYER（即使是用放率控制的），这些卡的每一张都给定了一个完井层的流动指数（ PID）和井底流压（ FBHP），因此，必须读这些卡的 NLAYER，所读的第一张卡应用于最上面的完井层（ PERF1），另外的卡应用于下面的各层，若给定了井的产率（ KIP=+1，+2 或+3），则将不用