1·程序概述

软件介绍

BOAST程序模拟以三维情况的等温达西流，它假设油藏中的流体是具有恒定成分，且物性只由压力决定的三相流体──油、气、水。这种对油藏流体的近似适用于世界上的大部分油、气藏。因此，BOAST能被广泛应用。例如：BOAST能模拟由流体的膨胀、驱替、重力驱和毛管力吸入机理进行的油和（或）气的开采。

能用BOAST处理的典型油田生产问题有：一次开采的研究，注水和（或）气保持压力的开采。注水、驱替等二次开采的评价。从技术上看：BOAST是一个隐式压力，显式饱和度（IMPES）有限差分数值模拟软件，在井动态方面，BOAST是井模型允许规定井的产率或压力，模拟过程中，用户可以任意增加或减少井数，在BOAST里有压力初始化的算法，泡点轨迹线，时间步长自动控制和检测解的稳定性的平衡法，此外，许多由用户控制的输出选择项也是有效的。

数据输入要求：

数据输入段分为两部分：初始化数据段和循环数据段。初始化数据包括：油藏模型网格的大小及形状、孔隙度、渗透率在油藏中的分布，流体的PVT数据，岩石的相对渗透率和毛管压力数据，油藏的原始压力及饱和度分布，解法说明和各种运行控制参数。

循环数据包括：模型中井的位置及其初始规定，在时间上继续模拟下去所需要的时间步长控制信息，单井产率和（或）压力动态表，完井和井下作业随时间的变化以及对模拟程序输出资料和类型和输出次数的控制。

本段将只讨论模型所要求的输入数据的类型的格式。

注意：下面有关数据输入格式（特别是对输入数组的规定）的注释适用于所有数据输入。

1·在本手册中的输入数据行称作“卡”，其含意是：在卡片输入计算机上表示计算机输入卡，而在输入数据的显示器或行式打印机上表示单独的一行，每张卡片的输入格式在紧随卡片描述的圆括号中说明。

2·标题卡先于每个主要的和许多次要的数据输入段读入。除非在文件中有别的说明，必须在每种情况下都读入卡片，为了使输入数据文件易读易编辑，将这些卡片设计成易看的图型。

3·在很多情况下，读入控制码来规定将要输入的数据的类型和数值个数，为便输入文件可使用选择项具有灵活性，必须读入控制码。

4·除标题卡外，所有数值都按固定格式输入，考虑到要让输入文件易读易编辑，大多数的输入格式为I5，F10.0或F8.0。

注意：在程序中，空格将按0读入

5·若某一特定参数要求输入所有网格的值，则必须按如下顺序输入：首先读第一层（K＝1）对每一层的数据输入是按行进行的，先读第一行（J＝1），对每一行以是以第一列（I＝1）开始的，第一行读入I＝1到II列的数值，读完第一行，再读第二行（J＝2），第二行仍从I＝1到II，以此类推，直至读完JJ行，再对第二层重复上述过程，直至读完KK层。

6·记住BOAST用的坐标原点在左边上角，因此每一层的网格分布如下：

I＝1 I＝2 ……

J=1

J=2

.

.

.

Z方向的值向下为增。

对运行黑油模型的提示：

在解决油藏工程和生产问题方面，BOAST是一个复杂的工程工具，为了更灵活地处理在油藏管理和生产过程中遇到的大量问题就更需要象BOAST这样的一个通用的且有些复杂的（因为程序中包含很多选择项）油藏模拟软件。

也许尽快熟悉如何运行BOAST的最佳方法是先浏览一下输入数据和熟悉输入文件的一般格式。

BOAST具有对时间步长自动控制的特点，建议在大多数运行中使用这种自动控制，当地层条件随时间变化不大时，这种特点可通过增加步长使程序最有效地工作，而当地下条件变化较大时（如泡点压力从未饱和状态变到饱和状态），较小的时间步长将使程序具有较大的稳定性，建议将时间步长自动控制的最小值设为0.1天。根据实际应用的要求还可选用更小的时间步长，例如，模拟用实验岩心注水时要求的最小时间步长为0.01天或更小，本手册提供的所有缺省都刻意于油田规模的实际应用，这种小的时间步长将只用于用户指定的一个时间步长里发生的最大压力和饱和度变化。

2·初始化数据

初始化数据卡描述的是油藏模型网格的大小和形状，以及孔隙度、渗透率、相对渗透率、毛管压力数据、流体的PVT值。初始压力和饱和度在油藏中的分布，同时还描述了使用的解法及各种控制、诊断参数。

这些数据卡只在模拟开始时读入一次，它们必须按如下的数据段中的顺序读入。

网格尺寸和形状

网格尺寸

1·标题卡（40A2）

2·重启动控制码 2I5

IW──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 读入第一层（K＝1）第一行（J＝1）的每个网格块的X方

向网格尺寸，并将同样的值赋给第一层其它行和其它层的

所有行X方向的对应网格（注：必须读II个值）

KDX＝＋1 读入第一层（K＝1）的每个网格块的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 每一层读入一个不变的厚度值，不同层可以有不同的厚

度。（注：必须读KK个值）

KDZ＝－1 读模型的每个网格Z方向的尺寸（厚度）（注：必须读II

×JJ×KK个值）

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＝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）

如果NUMDX＝0就不读这些数据卡

I── 被改网格的X坐标

J── 被改网格的Y坐标

K── 被改网格的Z坐标

DX── 网格块（I，J，K）X方向网格尺寸（DX）的新值。

注：必须读NUMDX卡

4·网格Y向尺寸（DX）的修改（3I5，F10。0）

如果NUMDY＝0就不读这些数据卡

I── 被改网格的X坐标

J── 被改网格的Y坐标

K── 被改网格的Z坐标

DY── 网格块（I，J，K）Y方向网格尺寸（DY）的新值。

注：必须读NUMDY卡

5·网格Z向尺寸（DZ）的修改（3I5，F10。0）

如果NUMDZ＝0就不读这些数据卡

I── 被改网格的X坐标

J── 被改网格的Y坐标

K── 被改网格的Z坐标

DZ── 网格块（I，J，K）Z方向网格尺寸（DZ）的新值。

注：必须读NUMDZ卡

第一层网格块顶部的深度

记住，在BOAST所使用的坐标系下，Z方向的值是向下递增的，因此，读入的深度应为低于用户选定的基准面的值（高于此基准面为负）

1·标题卡（40A2）

2·输入深度值的控制码（I5）

KEL── 输入控制码

KEL＝0 意指对第一层的所有网格块的顶部深度只读入一个定值

（即第一层为水平面）

KEL＝1 意指对第一层的每个网格块各读入一个深度值，共计读

II×JJ个值。

3·深度值（8F10。0）

ELEV── 网格块顶部

注意：如果KEL＝0，则只读一个值。如果KEL＝1，则读入II×JJ

个值，对第一行（J＝1〕读II个值，然后对第二行（J＝2）

又读II个值，如此读至第JJ行。除第一层外，各层网格块

顶部深度用其前一层网格顶部深度加上前一层网格厚度得

到：

即：TOP（I，J，K＋1）＝TOP（I，J，K）＋DZ（I，J，K）

孔隙度和渗透率的分布

孔隙度和渗透率

1· 标题卡（40A2）

2· 输入类型控制码（4I5）

KPH── 控制孔隙度数据输入的控制码

KKX── 控制X方向渗透率输入的控制码

KKY── 控制Y方向渗透率输入的控制码

KKZ── 控制Z方向渗透率输入的控制码

控制码 含意

－1 读入一个定值并赋给模型的所有网格块（只读一个值）

0 对模型KK层的每一层读入一个定值，不同层可以有不同的定

值。（必须读KK个值）

＋1 对模型的每一个网格块都读入一个值（必须读II×JJ×KK个

值）

3·孔隙度值（10F8。0）

孔隙度不是按百分数，而是按小数形式读入的。

如果KPH＝－1，则只读一个定值

如果KPH＝ 0，则读KK个值（一层一个值）

如果KPH＝＋1，则读II×JJ×KK个值（模型之第一个网格块都有

一个值）

4·X方向渗透率值（KX）（10F8。0）

渗透率值按毫达西（md）读入

如果KKX＝－1，则读一个定值

如果KKX＝ 0，则读KK个值（一层一个值）

如果KKX＝－1，则读II×JJ×KK个值（模型的每个网格块都有一

个值）

5·Y方向渗透率（KY）（10F8。0）

Y方向渗透率值（KY）（10F8。0）

渗透率值按毫达西（md）读入

如果KKY＝－1，则读一个定值

如果KKY＝ 0，则读KK个值（一层一个值）

如果KKY＝－1，则读II×JJ×KK个值（模型的每个网格块都有一

个值）

6·Z方向渗透率值（KZ）（10F8。0）

渗透率值按毫达西（md）读入

如果KKZ＝－1，则读一个定值

如果KKZ＝ 0，则读KK个值（一层一个值）

如果KKZ＝－1，则读II×JJ×KK个值（模型的每个网格块都有一

个值）

对孔隙度、渗透率分布的修改

1· 标题卡（40A2）

2· 要进行孔隙度和（或）渗透率值修改的网格数和打印控制码（5I5）

NUMP──要修改孔隙度的网格块数

NUMKX──要修改X方向渗透率值（KX）的网格块数

NUMKY── 要修改Y方向渗透率值（KY）的网格块数

NUMKZ── 要修改Z方向渗透率值（KZ）的网格块数

IPCODE── 打印控制码

IPCODE── 不打印修改了的孔隙度和渗透率分布

IPCODE── 打印修改了的孔隙度和渗透率分布

3·孔隙度的修改（3I5，F10。0）

如果NUMP＝0，则不读这张数据卡

I── 被改网格块的X坐标

J── 被改网格块的Y坐标

K── 被改网格块的Z坐标

PHI── 网格块（I，J，K）的新的孔隙度值

注：孔隙度必须以小数形式输入，必须读NUMP卡

4·X方向渗透率（KX）的修改（3I5，F10。0）

如果NUMKX＝0，则不读这张卡

I── 被改网格块的X坐标

J── 被改网格块的Y坐标

K── 被改网格块的Z坐标

KX── 网格块（I，J，K）的新的X向渗透率值。

注：必须读NUMKX卡

5·Y方向渗透率（KY）的修改（3I5，F10。0）

如果NUMKY＝0，则不读这张卡

I── 被改网格块的X坐标

J── 被改网格块的Y坐标

K── 被改网格块的Z坐标

KY── 网格块（I，J，K）的新的Y向渗透率值。

注：必须读NUMKX卡

4·Z方向渗透率（KZ）的修改（3I5，F10。0）

如果NUMKZ＝0，则不读这张卡

I── 被改网格块的X坐标

J── 被改网格块的Y坐标

K── 被改网格块的Z坐标

KZ── 网格块（I，J，K）的新的Z向渗透率值。

注：必须读NUMKX卡

传导系数的修改

注：记住在修改传导系数中对方向的规定是十分重要的。

例如，对网格块（I，J，K）：TX（I，J，K）指的是穿越网格块I－1与I的边界的流动。TY（I，J，K）指的是穿越网格块J－1与J的边界的流动以及TZ（I，J，K）指的是穿越网格块K－1与K的边界的流动。

1·标题卡（40A2）

2·要修改传导的网格块数，加之打印控制码（4I5）

NUMTX── 要修改X方向传导系数（TX）的网格块数

NUMTY── 要修改Y方向传导系数（TY）的网格块数

NUMTZ── 要修改Z方向传导系数（TZ）的网格块数

ITCODE── 打印控制码

ITCODE＝0，不打印修改了的传导系数分布

ITCODE＝1，印修改了的传导系数分布

3·X方向传导系数（TX）的修改（3I5，F10。0）

如果NUMTX＝0，则不读这张卡

I── 被改网格块的X坐标

J── 被改网格块的Y坐标

K── 被改网格块的Z坐标

TX── 网格块（I，J，K）X方向传导系数（TX）的新值。

注：必须读NUMTX卡

4·Y方向传导系数（TY）的修改（3I5，F10。0）

如果NUMTY＝0，则不读这张卡

I── 被改网格块的X坐标

J── 被改网格块的Y坐标

K── 被改网格块的Z坐标

TY── 网格块（I，J，K）X方向传导系数（TY）的新值。

注：必须读NUMTY卡

5·Z方向传导系数（TZ）的修改（3I5，F10。0）

如果NUMTZ＝0，则不读这张卡

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 PCGON

SAT── 相饱和度值。注意：SATI应为－0.01，而SATN必须是

1.10。每个饱和度均按小数读入

KRO── 油相相对渗透率（小数）

KRW── 水相相对渗透率（小数）

KRG── 气相相对渗透率（小数）

PCOW── 油/水毛管压力（at）

PCGO── 气/油毛管压力（at）

SAT指的是每一特定相的饱和度，即：在SAT＝0.2的数据行上：

KRO指的是当油的饱和度为20％时，油的相对渗透率。

KRW指的是当水的饱和度为20％时，水的相对渗透率

KRG指的是当气的饱和度为20％时，气的相对渗透率

PCOW指的是当水的饱和度为20％时的油/水毛管压力

PCGO指的是当气的饱和度为20％时的气/油毛管压力

流体的PVT数据表

1· 标题卡（40A2）

2· 泡点压力，未饱和油的性质，PVT表中的最大压力（5F10。0）

PBO── 油的原始泡点压力

VSLOPE── 未饱和油（P>PBO时的油）的粘度与压力关系曲线的斜

率。BSOLOPE── 未饱和油的体积系数与压力关系曲线的斜率。

注：与未饱和没的压缩系数不同，BSLOPE应为负值

RSLOPE── 未饱和油的溶解油气比与压力曲线的斜率（通常情况

下该值为0）

PMAX── 所有PVT表的最大压力记录

IREPRS──

3·标题卡（40A2）

4·油的PVT表（8F10。0）

P1 MUO1 BO1 RSO1

· ·

· ·

PMAX MUO@PMAX BO@PMAX RSO@PMAX

P── 压力（注：正如上表所示，最后一个压力记录必须是PMAX）

MUO── 饱和油粘度

BO── 饱和油地层体积系数

RSO── 饱和油溶解油气比

重要问题：在整个压力范围，油的物性必须作为饱和油的数据输入，一般来说，实验室的饱和油数据必须进行外推，直到高于所测得的泡点压力，以覆盖预计将在模拟运行中出现的最大压力范围。要求输入饱和油的数据。是因为在BOAST中使用了泡点轨迹线。请注意，只有当当地油藏压力增至高于原始泡点压力并且自由气进入该油藏时，才使用高于原始泡点压力的饱和油数据。（例如：向油藏注气保持压力）

5· 标题卡（40A2）

6· 水的PVT表（8F10。0）

P1 MUW1 BW1 RSW1

· ·

· ·

PMAX MUW@PMAX BW@PMAX RSW@PMAX

P── 压力（注：正如上表所示，最后一个压力记录必须是PMAX）

MUW── 水的粘度

BW── 水的地层体积系数

RSW── 水的溶解油气比

在黑油模拟中，通常假定气在水中的溶解度可以忽略不计，在这种情况下，对所有压力都设置RSW＝0，在BOAST中有水的这个PVT表是为处理如下情况用的── 从地层含水层采气的情况以及很多被认为气在水中的溶解度对问题的解有较大影响的情况。

7· 标题卡（40A2）

8· 气的PVT表及岩石的压缩系数（8F10。0）

P1 MUG1 BG1 CR1

· ·

· ·

PMAX MUG@PMAX BG@PMAX CR@PMAX

P── 压力（注：正如上表所示，最后一个压力记录必须是PMAX）

MUG── 气的粘度

BG── 气的地层体积系数

CR── 取决于压力岩石压缩系数

9· 标题卡（40A2）

10· 储罐流体密度（8F10。0）

RHOSCO── 储罐油的密度

RHOSCW── 储罐水的密度

RHOSCG── 标准状况下的气体的密度

压力和饱和度的初始化

BOAST含有两个有关压力和饱和度初始化的选择项，初始压力分布既或由考虑了油、气和油水接触位置及接触位置处的压力的平衡条件的程序计算出（选择项1），也可象在非平衡初始化情况下那样一个一个网格块地读入（选择项2），饱和度值（SO，SW，SG）既可对整个网格块作为定值读入（选择项1），也可在一个一个网格块的基础上读入整个SO和SW的分布，并由程序计算出SG的分布。SG＝1－SO－SW（选择项2）

1· 标题卡（40A2）

2· 压力及饱和度分布的控制码（2I5）

KPI── 压力初始化控制码

KSI── 饱和度初始化控制

控制码 含义

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— 油水接触处的深度（低于基准面的英尺数）

GOC— 油气接触处的深度（低于基准面的英尺数）

如果KPI=1（非平衡初始化），则对模型的每个网格块都读入一个压力值（注：必须读II*JJ*KK个压力）

4．饱和度初始化数据（10F8.0）

如果KSI=0（恒定饱和度情况），则

SOI— 把原始含油饱和度（小数）赋给模型的所有网格块

SWI— 把原始含水饱和度（小数）赋给模型的所有网格块

SGI— 把原始含气饱和度（小数）赋给模型的所有网格块

如果KSI=1（非平衡初始化）

SO数组：对模型的每一网格块都读入原始含油饱和度值

（注：必须读II*JJ*KK个含油饱和度值）

SW数组：对模型的每一个网格块都读入原始含水饱和度值

（注：必须读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．运行控制参数（I5，710.0）

NMAX— 运行结束前允许的最大时间步数

FACT1— 在自动控制时间步长的情况下，增大时间步长的因子（对

恒定时间步长，设置FACT1=1.0，FACT1的一个常用值为

1.25）

FACT2— 在对时间步长自动控制的情况下，减小时间步长的因子（对

固定时间步长设置FACT2=1.0）

TMAX— 在程序运行中模拟的最大实际时间，天（当模拟时间超过

了TMAX时，运行将结束）

WORMAX— 油田限定的最大油水比（当模拟的整个生产WOR超过了

WORMAX时，运行结束）

GORMAX— 油田的最大油气比限定值（当模拟的总的GOR超过GORMAX

时，运行结束）

PAMIN— 油田限定的最低平均压力值（当模拟的平均油藏压力低于

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迭代收敛性允许的最大压力变化。

TOL的一个典型值为0.1 。

TOL1— 决定何时改变OMEGA的参数。TOL1的一个典型值是0.001。

（注：如果TOL1=0，则上面赋给OMEGA的初值在整个程序运行中保持不变）

DSMAX— 一个时间步长所允许的最大饱和度变化（小数），如果在

一个时间步长下，任何网格块的任一相饱和度变化超过了

DSMAX，则时间步长将按FACT2缩小，DSMAX的一个典型值

是0.05。

DPMAX— 一个时间步长所允许的最大压力变化。如果在一个时间步

长里，任何一个网格块的压力变化超过了DPMAX，则时间

步长将按FACT2缩小，DPMAX的一个典型值是5at

时间步长和输出控制

注：在模拟开始时（即当TIME=0时），和在按（ICHANG*DT）定义的每个时间间隔的开始必须读以下的两个时间和输出控制卡。

1．时间步长和输出控制码（9I5）

IWLCNG—告诉程序在这个时间步长里是否要读井资料的控制码

控制码 含意

IWLCNG=0 在这个时间步长里不读井资料

IWLCNG=1 在这个时间步长里读井资料

注：如果IWLCNG=1，则必须读井资料卡，在下一时间步长时将使用新的井资料

ICHANG— 将使用输出控制资料和时间步长控制资料的时间步长数

IWLREP— 控制打印井报告的输出控制码

ISUMRY— 控制打印时间步长综合报告的输出控制码

IPMAP— 控制打印网格块压力图的输出控制码

ISOMAP— 控制打印网格块含油饱和度的输出控制码

ISWMAP— 控制打印网格块含水饱和度的输出控制码

ISGMAP— 控制打印网格块含气饱和度的输出控制码

IPBMAP— 控制打印网格块饱和压力的输出控制码（通常设置IPBMAP=0）（IPBMAP=0不打印报告/数组，IPBMAP=1在这个时间间隔的每一时间步长都打印报告/数组）

2．时间步长控制资料（3F10.0）

DT— 初始时间步长大小，单位：天。（如果饱和度或压力超过

限定值，则自动时间步长控制将减小初始时间步长）

注：使用卡1的控制码和卡2的资料的总的时间是（ICHANG*DT），

单位：天

DTMIN— 在此期间使用的最小时间步长，单位：天

（该值通常设为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— 这中井最上面的完井层的层位

NLAYER— 连续完井的总层数，从PERF1开始并包含PERF1，因此，

在一个五层模型网格里，如果一口井是在第2、3和5层

完井的，则设置PERF1=2和NLAYER=4，注意在这种情况

下，即使第4层没有射孔也必须包括在NLAYER中。在下

面的卡7中，第4层的PID值将被设置为0。

KIP— 是既说明井的类型，又说明井的生产（注入）动态是取决

于给定的产率还是取决于井底流压的控制码。该控制码还

说明程序将进行的是显示压力计算还是隐式压力计算，对

大多数情况，建议使用显示压力计算。

控制码 含意

KIP=3 气井­­­­---规定的产率

KIP=2 水井---规定的注水率

KIP=1 生产井（油井）---规定的产率

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）将生产三相流体，它们的

比例由油藏条件和井的约束条件决定。

注：这张卡的下面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），则将不用PWF，并给PWF读入零值。但是，要使用PI来计算该井的FBHP，并且在井报告中打印出FBHP，并非在任何情况下都用FBHP来控制井的动态。

PID—层的流动指数，其计算如下：

式中：k—层的渗透率，md

h— 层厚度，m

dx— 网格块X方向的尺寸，m

dy— 网格块Y方向的尺寸，m

— 井眼半径，m

S— 层的表皮系数

PWF— 层的井底流压（FBHP），at

只有当该井的KIP为负数时才用这个值，如果给定了该井的产率（KIP=+1，+2或+3），则设置PWF=0。

注：在一层一层的基础上改变PI值能及时地处理在任何地方发生的地层伤害或刺激。

重要问题：一旦一口井在一定层上完井了，即使该层或该井关闭了，在后面所有的井资料卡上都必须描述这一层，要关闭一层，就让该层的PID=0，要关一口井，就让该井的所有层的PID值为0。

重启动操作：当IW=1时，产生重启动文件F8.DAT。进行重启动运行时，将IR赋1，将F8.DAT改为F18.DAT，同时不赋压力和饱和度初始化值。

白天，或某个偶然的时间，就想写一点，然而等到到回家后坐在电脑前时，一点写的欲望都没有，一颗字也想敲，以前随手在手机里记录一点，现在手机又不好用了，也就写的更少了。以前坚持写了很多年的日记，成习惯了，总想到处又涂又写几下，后来工作了，也不方便写了，在日记本上写字的日子也就成了历史。虽然网上有blog，但有些东西还是属于自己内心深处不想让人知道的。现在想找个能让自己安静的畅开心扉的地方真的很难。

很喜欢的一首歌，感觉其中某些东西和自己能产生共鸣.Beyond：海阔天空 今天我寒夜里看雪飘过 怀著冷却了的心窝飘远方 风雨里追赶 雾里分不清影踪 天空海阔你与我 可会变(谁没在变) 多少次迎著冷眼与嘲笑 从没有放弃过心中的理想 一刹那恍惚 若有所失的感觉 不知不觉已变淡 心里爱..........

   当中石油和中石化受制于成品油价格和原油价格倒挂导致炼油亏损的时候，中海油得益于没有炼油企业，利润大幅飘红。公司上半年净利润同比增长８９．３％，达到创纪录的２７５．４亿元。其他两家公司利润分别下降３４．５％和７７．３％。中石化主要专注上游业务，业绩真是惨不忍睹。不过这一切都是原油价格或者说国家政策的错，如非如此，中石化必定笑得开心的很。

    ２００８年上半年，中海油油气产量达９２４０万桶油当量，同比增长８．３％。其中，原油净产量为７２９０万桶，天然气净产量为１１２５亿立方英尺，分别同比增长７．１％和１２．８％。

    公司平均实现油价达１０２．４９美元／桶，同比涨幅为７４．３％，平均实现气价为３．６６美元／千立方英尺，同比上涨１４．０％。

    由于国际油价走高，上半年中海油共缴纳石油特别收益金９７．５亿元，比去年同期增加近７５亿元。

    上半年中海油共获得６个油气新发现，还成功地评价了３个含油气构造。公司共有三个新项目，即涠洲１１－４北、西江２３－１和文昌油田群在上半年成功投产。傅成玉表示，下半年几个重大项目的投产将继续构成中海油中长期增长的重要动力。

    中海油上半年基本每股盈利达人民币０．６２元。董事会批准派发２００８年中期股息每股０．２０港元。

2008年8月26日17：40左右，旅行去都匀途中，在贵新高速距贵阳100公里处发生的一起车祸。天雨路滑，两辆出事的长安面包车似乎是一家人，车上有五六名未成年人，没看到驾驶员，估计伤势严重已送医，但白色面包车上还有一妇女满脸血不省人事。
生命宝贵，行车注意安全。“宁慢三分，不抢一秒”，这可不仅仅只是挂的标语。

    奥运会结束了，留给我们很多宝贵的财富和经验，也让民众心态开始变得平和，有了真正的大国心态，当然也留下了很多不甚好的记忆，不过随着时间，都会过去。奥运会同样留给我们精彩的奥运语录，这些语录反映着国家和人民的变化。有的很值得咀嚼和回味。

   “我去香港时，你们不用为我订酒店房间，给我一个马厩住就行了。”——国际奥委会主席罗格对北京奥运会香港马术赛场的“六星级”马厩大为赞赏。不能不说老外很会夸人，但你要真给他个马厩估计很难让他满意。

   “请大家相信我，我会回来的。”——退赛后刘翔在博客上写下这句话。其实是否能够回来并不是什么道德宣判的尊则，毕竟运动员的身体不是自己能左右的，不过可以欣慰的看到，国人的心态变得宽容多了。

   “对不起我的师弟们，因为我，他们牺牲了参加奥运会的资格。”——丢失男子10米气步枪金牌，朱启南三次落泪。比较人性化，不再对不起党对不起国家对不起人民了。

   “我脸上在笑，但我的心在流泪。”——银牌是中国女曲的历史最好成绩，但韩籍教练金昶伯还是很难过。虽然银牌是历史性的突破，但是，这个韩国来的教练并不满意。

   “尽200%努力拼巴西，要让对手害怕我们。”——小组赛最后一战，谢亚龙动员讲话。很好奇，谢亚龙有没有去问问巴西人怕了吗，可能会比较怕中国功夫踢碎睾丸吧。

    “虽然我们无权制定规则，只是让上场比赛就上场，但我还是觉得在正午踢球不是个好主意。”——虽然获得男足比赛银牌，但尼日利亚队主教练西亚西亚对在中午12时、在42℃的高温下比赛不解。想不到非洲朋友也怕热，不过我想可能对手更怕热吧。

    “如果换作是我，我也许会去找个心爱的人结婚，刘翔也应该这样。人生并不只有跑道和栏架。”——罗伯斯安慰因伤退赛的刘翔.确实如此，生活不仅仅是跑到，让自己的生活丰富点，多彩点，不是坏处。为自己或者不要为那些期盼你成为一个苦行僧一样的运动员而活着。

• 编年史【石油地质工作者博客】
• 李传亮老师 || 勃勃娱乐投资|| 非常秀
• 六度空间|| 爱博客|| 子说博客
• 涌泉博客 ||Mymike ||许许如生
• 尖笔头 ||纯粹灵魂 ||小天博客
• 帝王难民 ||花果山寨||LeoChung
• 不是俺||虫儿的童话|| 流浪者的心情
• 星空NET|| Nipao博客 ||安拉博客
• {博客派} ||中国时尚街拍|| Copterfly
• 和意博客|| 好心情Blog ||博X's Blog
• 破草屋的大杂院||水自流|| 左右单行道
• Marshal的博客|| HR'Blog||起步者'Blog
• 黔首报||博客久久|| 非礼勿视
• 流金岁月|| Oh-Cool ||积水成渊|
• 迷走部落格||网海天涯|| 暗潮天空
• 岩苏博客|| 北街博客||飞信博客
• 魔术空间||朝闻博客|| 老Q的博客
• 37计|| 听风看海||爱地质爱生活|
• 林网博客 ||新开博客|| 主角博客
• 傻猫网络日志|| 梅菜扣肉||完美生活
• 怪鸟博客 很抱歉，更换主题丢失了6月份后的部分链接，如果你发现自己不在上面请告诉我。

　　刚开始建立博客的时候，没有友情链接，很不容易有人肯给做一个，几乎每天点过去看几遍。现在链接多了，也是要坚持每个星期都要去看一趟，不然这个友情链接就变成无情连接了。链接有朋友来申请的，也有我出去看见好博客申请的。随着链接的增长，也遇到这样的问题，就是链接太多，放在首页超级难看，博客名字有长有短，有中文有英文，还有混合体的。真希望国家能立法，博客名字只能四个字，^_^。

　　虽然不是全部，但肯定不少朋友是非常在意链接的ｐｒ值，所以，尽管放在首页很难看，我也一直没有把它单独开个页面。不过这终究不是个办法，我望着下面这些链接，不禁长叹，怎么摆弄你们好呢？

o(∩_∩)o...哈哈，各位朋友，在这个页面，一人送个外链。

• 编年史【石油地质工作者博客】
• 李传亮老师 || 勃勃娱乐投资|| 非常秀
• 六度空间|| 爱博客|| 子说博客
• 涌泉博客 ||Mymike ||许许如生
• 尖笔头 ||纯粹灵魂 ||小天博客
• 帝王难民 ||花果山寨||LeoChung
• 不是俺||虫儿的童话|| 流浪者的心情
• 星空NET|| Nipao博客 ||安拉博客
• {博客派} ||中国时尚街拍|| Copterfly
• 和意博客|| 好心情Blog ||博X's Blog
• 破草屋的大杂院||水自流|| 左右单行道
• Marshal的博客|| HR'Blog||起步者'Blog
• 黔首报||博客久久|| 非礼勿视
• 流金岁月|| Oh-Cool ||积水成渊|
• 迷走部落格||网海天涯|| 暗潮天空
• 岩苏博客|| 北街博客||飞信博客
• 魔术空间||朝闻博客|| 老Q的博客
• 37计|| 听风看海||爱地质爱生活|
• 林网博客 ||新开博客|| 主角博客
• 傻猫网络日志|| 梅菜扣肉||完美生活
• 怪鸟博客 

“将骨瓷拿在手里的时候，你会觉得自己的双眼仿佛具有穿透力，对着光观察瓷的通透，用手感觉瓷的细腻和坚硬，再用手指轻轻一弹，听到“叮”的一声脆响，犹如久叩的璃扉终于得到回应，灵犀刹那间相通，这便是今后放在唇边，朝夕相伴的美瓷了。”

一种美丽的瓷器，一种中国传统陶瓷艺术，在一个特定的历史年代，那便是照片中本周“景德镇陶瓷新大陆礼品瓷厂专卖店”推出的，骨瓷手绘水点桃花球形壶茶具。这套作品瓷质细腻通透，器型美观典雅，彩面清高别致，釉色白瑕如玉、晶莹润澈，不论是其白度、透光度还是手绘画面及色彩均达到了较高水准，充分展现出当年“水点桃花”独具的“薄如纸、白如玉、明如镜、声如磬”的瓷艺境界。而骨瓷的材质，又使这套作品更洁白、细腻、通透、轻巧和清新脱俗。“水点桃花”瓷诞生于特定的历史年代，继承了釉下彩的烧制传统和釉上彩的装饰工艺，创烧出了“洁白如玉、晶莹剔透”的釉下五彩瓷器，是景德镇20世纪制瓷工艺上的一个里程碑。因而“水点桃花”专用瓷被被誉为“天下第一瓷”，收藏界称作“红色官窑”。

【什么是水点桃花】
说到“水点桃花”这个词，可能年轻朋友不太了解，但上了些年纪的都知道。就算不知道“水点桃花”也一定有听说过“7501瓷”、“毛瓷”、“毛主席专用瓷”。

1949年后，中南海里用的瓷器，主要都是江西景德镇烧制的，中央机构和国务院机构用的瓷器，分别以在碗底和杯底烧上的“甲”和“政”（政务院）两个字来区别，这种用瓷延续到20世纪70年代。1975年，中央办公厅考虑毛泽东年事已高，吃饭咀嚼很慢，又由中南海丰泽园搬到游泳池住，厨房离用餐的地方有一段距离，常常出现饭菜送来了，不能当时用餐，把饭菜放凉了的情况。所以，为毛泽东制作既保温、防尘，又实用美观的优质餐具的工作就提上了议事日程。1975年4月，中央办公厅向江西省委下达任务，要求原轻工业部陶瓷工业科学研究所在3个月内完成毛泽东用瓷的研制任务。这个任务是该所1975年的第一项任务，故称7501任务。这次研制的产品称7501瓷，又称“毛主席专用瓷”，或简称“毛瓷”。

7501瓷的外形装饰是高级工艺美术师李雨苍设计完成的。7501瓷器的外饰在风格上力求简洁、典雅、端庄、稳重，有中国特色，线条流畅优雅，器形饱满，古朴大方，达到了艺术与实用的完美统一。器具画面有水点桃花或水点梅花，正式绘制是把水点程序按步骤划分为：拍图、点玻白、点洋红、画枝干、画叶、点托、结果（点蕊）、填色、贴底款、检查下位、烤花等；从1975年5月25日——8月31日，共烧制22窑，入窑14103件，实际成品4200余。从中精选出千余件“7501”专用瓷于当年9月送达北京中南海。中南海收到这批瓷器后，时任中央办公厅主任的汪东兴确曾下令将未送中南海的“7501瓷”制品全部销毁。但主持这批瓷器研制的原轻工部陶瓷研究所并没有按中央办公厅的意见办，他们认为，送到中南海的瓷器在使用过程中，肯定会有损耗，必须留下一批作为将来的补充储备。另外，这批瓷器不但凝聚了研究所200多人半年多时间的心血，它还是瓷中精品，代表了建国以来瓷器制作的最高水平，大家实在舍不得将其销毁。最后，研究所决定将其全部封存，使这批在中南海以外的“毛瓷”幸运地保存了下来。

1980年，在毛泽东逝世了数年以后，研究所按照以往的惯例，在当年的春节，从库房里拿出了一部分封存的“毛瓷”作为职工福利，每人分发了10件左右，近3000件的“毛瓷” 就这样流入了社会。在流入社会的仅3000件“毛瓷”中，有600余件已被一位新加坡收藏家收入囊中，现在理论上能在世上流传的“毛瓷”最多不过2000件左右。1996年，首次亮相拍卖会的87件“7501”毛主席专用瓷总成交额达到882.86万元，每件平均价约为10万多元；1997年广州嘉德拍卖会上，一件湖南醴陵釉下彩红月季花碗虽有釉裂，但因收藏价值高，最终以8.8万元成交；1998年，香港一家拍卖行以178万元的天价拍出一对“毛瓷”茶具对杯。

“7501毛泽东主席专用瓷” 诞生于特定的历史年代，继承了釉下彩的烧制传统和釉上彩的装饰工艺，创烧出了“洁白如玉、晶莹剔透”的釉下五彩瓷器，是景德镇20世纪制瓷工艺上的一个里程碑。这批专用瓷被称为“主席用瓷”，被誉为“天下第一瓷”，收藏界称作“红色官窑”。

【什么是骨瓷】
陶瓷起源于中国，但骨瓷始创于英国，曾长期是英国皇室的专用瓷器，尊贵而华美。就在如今，骨瓷也是身份与地位的象征。而一件骨瓷作品，应该用触觉、视觉和听觉与之交流。骨瓷手感好、轻而薄、强度高、韧性好，敲打起来声音清脆，并且晶莹剔透，色泽亮丽，有着说不出的典雅和皇家气派。

所谓骨瓷，是在其黏土中加入牛、羊等食草动物骨灰（以牛骨粉为佳）而得名，是环保的绿色消费品。骨质瓷的形成主要依靠氧化硅、氧化铝和氧化钙，其中氧化钙的含量越高，色泽越好。在自然界中，氧化钙的来源不多，所以选择动物的骨粉作为氧化钙的来源。由于骨瓷工艺复杂，因此只有英国、中国、日本、德国、俄罗斯和泰国有生产。

关于骨瓷的发明，大多的说法是1794年英国人在制作瓷器的时候，也许为了让瓷器更轻巧，便尝试加入动物的骨粉和骨灰。当第一件骨瓷出品的时候，它在灯光下的光泽不是散射，而是透射的。于是，全世界赞叹之后都开始将牛、羊、猪骨磨成粉，将骨粉加入到瓷器里。这其中，牛骨粉因为能使骨瓷最坚韧，因而最有价值。

骨瓷色泽呈天然骨粉独有的自然奶白色，白度高、透明度高、瓷质细腻。骨质瓷在烧制过程中，对它的规整度、洁白度、透明度、热稳定性等诸项理化指标均要求极高，由此废品率很高。由于用料考究、制作精细、标准严格，所以价值高于其它瓷种。独特的烧制过程和骨碳的加入，使瓷土中的杂质被消除，骨瓷显得更洁白、细腻、通透、轻巧，极少瑕疵，并且比一般瓷器薄，在视觉上有一种特殊的清洁感，强度高于一般瓷器，是日用瓷器的两倍。骨粉的含量越高，黏土的成分就相对降低，在制作过程中就越易烧裂，在成形上需要更高的技术，增加了烧制难度，所以更加珍贵。有人做过这样的实验，把4只骨瓷杯子垫在一辆劳斯莱斯车轮下，这4只晶莹的杯子竟能胜任支撑房车的重量！

【骨瓷挑选与养护】
1.一定用手清洗，不可用洗碗机。若真不愿意用手洗，则要选择有“瓷器及水晶”类洗涤功能的洗碗机。
2.有金边的餐具，勿放入微波炉，以免腐蚀。
3.洗涤PH值必须在11-11.5之间。
4.用清水洗净时，水温不超过80摄氏度。
5.勿将热杯直接浸入冷水中，以免温度迅速改变损伤瓷质。
6.如有刮花，可以用牙膏略打磨。
7.如有茶渍，可用柠檬汁或食醋清洗。
8.最好不要骤泠骤热使用，以免炸裂。
9.不要用明火直接烧用。

【陶瓷小知识】
如何分辨手绘图案和机器印刷图案
只要把两个杯子放在一起比较就会发现，机器印刷的图案完全是一模一样的，因为都是一个模子里印出来的了；而手工绘制的图案无论画师如何绘制，都不可能画的一模一样的。

    引用一个网友的话说，伦敦八分钟之前我以为英国实在欧洲，之后我怀疑是不是在欧洲。如果说开幕式过于奢华浪费，哪么这个200万英镑预算的八分钟是不是英镑变卢布了，钱都花在那里了。当然也有很多精英们号称看懂了，并且被这个表演感动了，我深表佩服。但我很好奇，如果中国号称花了2000W人民币搞8分钟。出现的是：
　　一辆公交车开过来，一群骑自行车的，一堆人车来了不排队，报纸乱扔，把残疾人挤一边去。汉族小女孩踩着少数民族的背上车，李宇春和崔健献唱，姚明把一个篮球扔出去，最后车开走了，精英们会怎么解读？坦率的说，精英们总是双重标准。

    说说几个印象深刻的情节吧，如同生化危机般向巴士涌上去的僵尸，彻底把我搞败，如果真的一定要表现什么，那也不用非得摆出这样的造型吧，跟僵尸一般。等巴士的时候，手中报纸乱扔...这个。贝克汉姆的角色太让我吃惊，搞的他象某种珍贵稀有动物一样，放在橱窗里亮个相。似乎英国人想表达这个意思，“你们看啊，我们有这个，有这个啊，贝克汉姆啊！！”

转自BBC论坛

   The closing ceremony was just fine apart from the British bit which was simply atrocious!
　　Why did we have some prehistoric rock guitarist and a talent show winner miming to some abysmal sound track? What has David Beckham got to do with the Olympics? Hardly portrays a young, vibrant and exciting vision of what the 2012 games might offer.
　　If this is an example of how good London's Olympics is going to be, then perhaps Boris should ask the Chinese to keep the Olympic flame in Beijing.

　　这哥们最后说，‘如果伦敦奥运会将来就是这样，还不如叫中国人继续留着旗子呢！’

另外这个会徽花费２０万英镑，来欣赏一下