中imex中文操作手册（三）

目 录 流体组分性质数据段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附录320207流体组分性质数据段*MODEL流体模型_目的： *MODEL 标志使用的流体组份模型.格式： *MODEL (*BLACKOIL) (*LIGHTOIL) (*OILWATER) (*MISCG) (*MISNCG) (*POLY) (*POLYOW)定义：*BLACKOLL 使用黑油模型,模拟油,气,水的流动.*LIGHTOIL 使用黑油模型,模拟油,气,水的流动。分离 Bo, Rs 和两个参数 alpha 和 beta用于将储罐油产量转换为油藏条件产量， alpha 是气体收缩因子，而 beta 是在储罐条件下气中析出的液量，这个选项对某些种流体是有用的，它们在油藏条件下是油，但是很轻，在地面条件时有液体从气中析出，这个选项使用的方法大致来源于“用于挥发油藏的改进黑油模型”。*OILWATER 使用油,水两相模型，不模拟自由气或溶解气的变化.*MISCG 使用带有驱动气的拟混相模型，假设溶解气始终保持溶解状态.*MISNCG 使用没有驱动气的拟混相模型,假设注入气(如果存在)组份与溶解气相同. *POLY 使用聚合物模型，模拟油,气,水和聚合物的流动过程. *POLYOW 使用聚合物模型,不考虑气相流动或溶解气的变化过程.缺省： 此项为要求关键字,不得缺省。条件： 这一关键字必须位于组份性质数据段的开头,紧随于油藏描述数据段之后.说明： 在IMEX中,有 6 种不同的流体模型,但每次只能使用其中之一,当使用 *OILWATER 关键字时,对于每个网格,模型只求解两个方程,在这种情况下,不存在气相,或不模拟气相. 在使用带有驱动气的拟混相模型 *MISCG 时,溶解气将始终保持溶解状态.用户使用 *GORINT 关键字输入初始溶解油气比,这些数值在计算中将保持不变,而 PVT 表中的泡点压力与溶解油气比曲线用于描述驱动气,这样溶解气总量就是由两部份组成的,由 *GORINT 给出的溶解气和 PVT 关系给出的驱动气溶解量,即 总油气比=初始油气比驱动汽油气比PVT 数据： 可直接使用 *PVT 关键字输入 PVT 表，也可选择使用 *DIFLIB 关键字输入差异分离表，如果使用了 *DIFLIB 选项，那么在这个数据段中也必须存在 *BUBBLE，*BOBF 和 *RSIF 关键字，在数据文件中，*PVT 或 *DIFLIB 必须有一个存在。 如果使用了 *LIGHTOIL 选项，那么必须使用 *PVTLO 选项确定 PVT 性质。*TRES 油藏温度 _目的： *TRES 表示输入油藏温度. 格式： *TRES value定义：value 油藏温度值。(C|F|C)缺省： 要求关键字，没有缺省值。条件： 这些关键字必须位于组份性质段，如果在 PVT 表中,使用 *ZG 选项输入气体压缩系数,则此项为要求关键字.说明： 这里确定的油藏温度用于将气的 Z 因子 Zg 转换为气体膨胀系数 Eg: 假设储罐条件为 14.7psi,(101.325 kPa) 和 60 F,(15.56C). 那么以矿场单位表示: Eg = 198.588 * P/ZT (SCF/Res.BBL) P 的单位是 psia T 的单位是 R 以 SI 单位表示: Eg = 2.84935 * P/ZT (Std.m3/Res.m3) P 的单位是 kPa T 的单位是 K 例： *TRES 180.0 对于油藏温度，可接受的数值范围是： SI Field Lab C F C min -100 -100 -100 max 1000 1000 1000*PVT油和气的 PVT 表_目的： *PVT 说明开始输入油和气的 PVT 表.表格式： *PVT *EG set_number p rs bo eg viso visg (co) (srftn) . . . . . . . . . . . . . . . . . . . . . . . . 或者 *PVT *BG set_number p rs bo bg viso visg (co) (srftn) . . . . . . . . . . . . . . . . . . . . . . . . 或者 *PVT *ZG set_number p rs bo zg viso visg (co) (srftn) . . . . . . . . . . . . . . . . . . . . . . . . 定义：set_number 用于这个特定 PVT 区的标号,可通过在 *PTYPE 关键字输入时使用,对不同的区域网格设置计算时使用的 PVT 表.*EG 关键字说明将使用气的膨胀系数,如果 *EG, *BG ,*ZG 都不存在于数据文件之中，则此项为缺省关键字.*BG 关键字说明使用气的体积系数而不使用气的膨胀系数。*ZG 关键字说明使用气的压缩因子而不使用气的膨胀系数，使用这一选项之前必须输入油藏温度 *TRES. p 压力(kPa | psi | kPa | kg/cm2)。rs 在压力 P 时,饱和油的溶解油气比(m3/m3|scf/STB|cm3/cm3)，如使用驱动气选项,则此项为驱动气的溶解油气比.bo 饱和油在压力 P 时的体积系数(m3/m3|RB/STB|cm3/cm3)。eg 压力 P 时的气膨胀系数(m3/m3|scf/RB|cm3/cm3).当使用驱动气选项时,此项为驱动气的膨胀系数. bg 压力 P 时的气体积系数(m3/m3|RB/scf|cm3/cm3).当使用驱动气选项时,此项为驱动气的体积系数.zg 压力 P 时的气压缩因子,无因次量.如使用驱动气选项,则此项为驱动气的压缩因子. viso 压力 P 时的饱和油粘度(mPa.s|cp|cp).visg 压力 P 时的气粘度(mPa.s|cp|cp).如使用驱动气选项,则此项为驱动气粘度. srfth 油气界面张力(dyne/cm|lbf/ft|dyne/cm).必须单调递减并且大于零. co 油的压缩系数,必须保持常数或单调递增.也可选择使用 *CO 关键字而不使用此项输入(1/kPa).缺省： 此项为要求关键字,没有缺省值。*PVT，*PVTLO 或 *DIFLIB 三者之一必须存在于数据文件之中。条件： 此项关键字必须位于组份性质数据段。 此表必须按压力增加的顺序输入,并且应当扩展到初始泡点压力之上,表面张力为可选择项.说明： 一般来说,这张表必须至少输入两行,有时则要求有更多的行输入。 当使用油水两相模型时,有两种方法可用于确定初始油气比.第一种方法是通过 PVT 表,第一行必须是饱和压力行并且只需输入一行,使用第一行输入的油气比作为初始油气比；第二种方法是利用 *GORINT 关键字确定初始油气比.在输入时如未发现 *GORINT 关键字,则使用第一种方法. 模拟中使用的典型 PVT 数据曲线由附录 D 的图7(a),(b),(c),(d)和(e)中的实线给出，注意曲线包括无限气存在的情况(饱和情况)。 图7(a),(b),(d)的虚线表示当压力超过泡点压力 PB 时得到的曲线，假设它为线性的。 在泡点压力以上 Rs 曲线的斜率对于零，IMEX 自动使用一种变量替换技术模拟未饱和情况，见附录 A。 PVT 表中的最低压力应低于任何生产井的压力限制,正常情况下应选择零值.而最高压力应在任何注入井所确定的压力之上. 程序将利用函数关系计算水的粘度和体积系数. STC 表示标准的储罐条件。(14.7 psia 和 60 F| 101.3 kPa 和 288.15 K | 1.03 kg/cm2 和 288.15 K) rs = STC 溶解气体积 / STC 油的体积 bo = 油藏条件下油的体积 / STC 油的体积 eg = STC 气的体积 / 油藏条件下气的体积 bg = 油藏条件下气的体积 / STC 气的体积 缺省为 *EG，使用 *BG 关键字为 bg 输入数值。 set_number 用于对油藏的网格设置不同的 PVT 表,如只使用一个 PVT 表,则可忽略表号，因为缺省是 1。 如使用多个 PVT 表,则第一个 PV T表号必须是 1,而其它表号顺序增加。 例如： *PVT 1 * 矿场单位 * p rs bo eg viso visg 14.7 1.0 1.062 6.0 1.04 0.0080 264.7 90.5 1.15 82.692 0.975 0.0096 514.7 180.0 1.207 159.388 0.91 0.0112 1014.7 371.0 1.295 312.793 0.83 0.0140 2014.7 636.0 1.435 619.579 0.69 0.0189 2514.7 775.0 1.5 772.798 0.64 0.0208 3014.7 930.0 1.565 925.926 0.594 0.0228 4014.7 1270.0 1.695 1233.046 0.51 0.0268 5014.7 1618.0 1.827 1540.832 0.449 0.0309 9014.7 2984.0 2.36 2590.674 0.203 0.0470 对于压力可接受的数值范围是： SI Field Lab Mod. SI kPa psi kPa kg/cm2 min 0.0 0.0 0.0 0.0 max 1.0E+10 1.45E+9 1.0E+10 1.01E+8 对于溶解气油比可接受的数值范围是： SI Field Lab m3/m3 scf/STB cm3/cm3 min 0.0 0.0 0.0 max 1.0E+3 5,617 1.0E+3 对于地层体积系数可接受的数值范围是： SI Field Lab m3/m3 Rb/STB cm3/cm3 min 0.5 0.5 0.5 max 5.0 5.0 5.0 对于气的膨胀系数可接受的数值范围是： SI Field Lab m3/m3 scf/RB cm3/cm3 min 1.0E-20 5.617E-20 1.0E-20 max 1.0E+3 5617 1.0E+3 对于油的粘度可接受的数值范围是： SI Field Lab mPa*sec cp mPa*sec min .0095 .0095 .0095 max 86,400. 86,400. 86,400. 对于气的粘度可接受的数值范围是： SI Field Lab mPa*sec cp mPa*sec min 9.5E-5 9.5E-5 9.5E-5 max 86,400. 86,400. 86,400. 对于界面张力可接受的数值范围是： SI Field Lab dyne/cm lbf/ft dyne/cm min -1.1E-20 -6.852E-25 -1.1E-20 对于油的压缩系数可接受的数值范围是： SI Field Lab Mod. SI 1/kPa 1/psia 1/Pa 1/(kg/cm2) min 0.0 0.0 0.0 0.0 max 1.0E-3 6.89E-3 1.0E-3 1.0E-2 在 IMEX 中有一个总得烃类压缩性检查，这个检查确保了当气从溶解状态析出时结果具有物理意义，油气混和引起的体积变化，即： Bg*dRS/dP - dBo/dP 0.0 必须成立， 或者以 i 和 i+1 的方式输入： Bg(i+1)*(Rs(i+1)-Rs(i)/(P(i+1)-P(i) 大于或等于 (Bo(i+1)-Bo(i)/(P(i+1)-P(i) 如果这个检查出错将打印警告信息，然而将继续运行。 *PVTLO轻质油和气的PVT表（任选）_目的：表示开始输入轻质油和气的PVT数据表数据表：*PVTLO*EGset_number prsboegvisovisgboocalphabeta或者*PVTLO*BGset_numberprsbobgviso visg booc rsoc alpha beta或者*PVTLO*ZGset_numberprsbozgviso visg booc rsoc alpha beta定义：set_numberPVT表号，可通过在*PTYPE关键字输入时使用，对不同的区域节点设置计算时使用的PVT表。*EG说明使用的是气体的膨胀系数。如果EG，*BG和*ZG都不存在，则此项为缺省关键字。*BG说明使用的是气体体积系数，而不是膨胀系数。*ZG说明使用的是气体的压缩因子，使用这一选项之前必须输入油藏温度*TRES。p压力（kpapsikpakg/cm2）rs由差异分离试验得到的饱和油在压力p时的溶解气油比（m3/m3scf/STB）。该rs用于流动方程的计算。bo由差异分离试验得到的饱和油在压力p时的体积系数（m3/m3RB/STB），该项用于建立流动方程。eg压力p时的气体膨胀系数（m3/m3scf/RBcm3/cm3）。当使用驱动气选项时，此项为驱动气的膨胀系数。bg压力p时的气体体积系数（m3/m3 RB /scf cm3/cm3）。当使用驱动气选项时，此项为驱动气的体积系数。zg压力p时的气体压缩因子，无因次量。若使用驱动气选项，此项为驱动气的压缩因子。viso压力p时的饱和油粘度（mPascpcp）。visg压力p时的气体粘度（mPascpcp）。当使用驱动气选项时，此项为驱动气的体积系数。booc压力为p时的分离器饱和油的体积系数（m3/m3 RB /STB cm3/cm3）。该选项用于地面流量与油藏条件流量之间的转换。rsoc压力为p时的分离器饱和油的溶解气油比（m3/m3 RB /STB cm3/cm3）。该选项用于地面流量与油藏条件流量之间的转换。alpha油藏条件下的气体体积与分离器条件下气体体积之比，即气体收缩率（m3/m3 scf /STB cm3/cm3）。beta每单位自由气产生的分离器油（m3/m3 STB / scf cm3/cm3）。缺省：任选关键字。无缺省值。条件：该关键字必须位于组分性质数据段。此表必须按压力增加的顺序输入，并且应当扩展到初始泡点压力之上，表面张力为任选项。解释：一般来说，这张表至少要输入两行，也可以按需要输入更多的行。模拟软件中所用的典型PVT数据曲线列在附录D、图7(a)、(b)、(c)、 (d)和(e)中。这些曲线是在饱和条件下得出的。图7(a)、(b)、(d)中的虚线是通常状况下的压力大于泡点压力时的曲线，这些曲线假设为直线。Rs曲线中泡点压力以上部分的梯度为0。IMEX自动使用一个变量代换来模拟未饱和油藏（见附录A）。PVT表中的最低压力应低于任何生产井的压力限制，通常情况选择零值。而最高压力应大于任何注入井的注入压力。程序将利用函数关系计算水的粘度和体积系数。STC表示储罐条件（14.7 psi 600F 101.3 kPa 288.15K 1.03kg/cm2 288.15K）。rs=bo=eg=bg=缺省为*EG，可以使用*BG来代替*EG。将输入的产量转换为模型使用的产量：若 GORi*Qgialpha/betaQgi=Qoi*alpha/betaQom=(Qoi-Qgi*beta/alpha)/(bo/booc)*(1-rsoc*beta/alpha)若 GORialpha/betaQom=0.0Qgm=Qom*(rs-(rsoc/aloha)+Qgi/alpha若GORi0.0 成立以表格形式表示时，满足：Bg(i+1)*(Rs(i+1)-Rs(i)/(P(i+1)-P(i)(Bo(i+1)-Bo(i)/(P(i+1)-P(I)若检查出不符合上述条件,则显示出警告信息。*DIFLIB油气的差异分离数据表（任选）_目的：表示开始输入油和气的差异分离数据表。表：*DIFLIB *EG set-number*BUBBLEpressure*BOBFbobf*RSIFrsifprsdbodegvisovisg(co)(srftn)或者：*DIFLIB *BGset_number *BUBBLFbobf *RSIF rsifprsdbodegvisovisg(co)(srftn)或者：*DIFLIB *ZGset_number*BUBBLFpressure*BUBFbobf*RSIF rsifprsdbodegvisovisg(co)(srftn)定义：set_numberPVT表号，可通过在*PTYPE关键字输入时使用，对不同区域节点设置计算时使用的PVT表。*EG，*BG，*ZG的定义和*PVTLO中所定义的相同。*BUBBLE*DIFLIB表中的泡点压力，用于将*DIFLIB数据转换成*PVT数据。pressure泡点压力值（kPapsikPakg/cm2）。BOBF泡点处油的体积系数，通过在适当的地面分离器条件闪蒸来确定。bobf泡点处的油体积系数值，（m3/m3RB/STB cm3/cm3）。*RSIF泡点处溶解气油比，通过在适当的地面分离器条件闪蒸来确定。rsif泡点处溶解气油比，（m3/m3 scf /STB cm3/cm3）。P压力（kPapsikPakg/cm2）。rsd压力为p时饱和油的差异溶解气油比（m3/m3scf/STB cm3/cm3）。当使用驱动气选项时，该值是驱动气溶解气油比Rs。bod压力为P时饱和油的差异地层体积系数（m3/m3RB/STB cm3/cm3）。eg压力为P时的气体膨胀系数（m3/m3 scf/STB cm3/cm3）。当使用驱动气选项时，该项是驱动气的膨胀系数。bg压力为P时的气体体积系数（m3/m3 RB/scf cm3/cm3）。当使用驱动气选项时，该项是驱动气的体积系数。zg压力为P时的气体压缩因子（无量纲）。当使用驱动气选项时，该项是驱动气的压缩因子。viso 压力为P时的饱和油的粘度(mPascpcp)。visg压力为P时的气体粘度(mPascpcp)。当使用驱动气选项时，该项是驱动气的粘度。srftn油-气界面张力。该项必须单调递减且大于0(dyne/cm1bf/ftdyne/cm)。co油相压缩系数。该项必须保持不变或单调递增。也可以使用*co关键字来输入油相压缩系数（1/kPa）。缺省：任选关键字。无缺省值。*PVT，*PVTLO或*DIFLIB至少存在一个。条件：该关键字必须位于组分性质数据段。此表必须按压力增加的顺序输入，并且应当扩展到初始泡点压力之上。表面张力为任选项。解释：rsd=bod= eg= bg=缺省为*EG。也可以用*BG代替EG。*DIFLIB数据转换成*PVT数据：bo=bod*（bobf/bobd）rs=rsif-（rsid-rsd）*（bobf/bobd）其中：bo - 校正的PVT Bo；bod - *DIFLIB表中的bod项。bodf - *BOBF关键字的bodf项。dobd - *DIFLIB表中*BUBBLE关键字确定的压力时bobd项。rs - 校正的PVT Rs。rsif - *RSIF关键字rsif项。rsid - *DIFLIB表中*BUBBLE所确定的压力时的rsid项。rsd - *DIFLIB表中rsd项。表号用于给油藏中的网格节点设置不同的*DIFLIB表。若只有一个PVT区，那么表号可以忽略，缺省为1。若使用多个PVT区域，第一个DIFLIB表的表号必须为1，其它的表号顺序增加。PVT表号在*PTYPE关键字输入时使用。举例： *DIFLIB 1 * 矿场单位 *BUBBLE 2014.7 *RSIF 600.0 *BOBF 1.35 * p rsd bod eg viso visg 14.7 1.0 1.062 6.0 1.04 0.0080 264.7 90.5 1.15 82.692 0.975 0.0096 514.7 180.0 1.207 159.388 0.91 0.0112 1014.7 371.0 1.295 312.793 0.83 0.0140 2014.7 636.0 1.435 619.579 0.69 0.0189 2514.7 775.0 1.5 772.798 0.64 0.0208 3014.7 930.0 1.565 925.926 0.594 0.0228 4014.7 1270.0 1.695 1233.046 0.51 0.0268 5014.7 1618.0 1.827 1540.832 0.449 0.0309 9014.7 2984.0 2.36 2590.674 0.203 0.0470*COT作为压力函数的油压缩系数_目的： *COT 用以说明，当油的压缩系数即是泡点压力 Pb 的函数，又是压力 P 的函数时，输入的油压缩系数表。 格式： *COT set_number p co . . . . . .定义：set_number 对应于 co 表的 PVT 区。P 对应于 PVT 表中泡点压力以上的压力(kPa|psi|kPa|kg/cm2).这个表的第一行是泡点压力时的油压缩系数值，后面输入的压力必须大于泡点压力并且是单调递增的。 co 对应于压力的油压缩系数。(1/kPa|1/psi|1/kPa|1/(kg/cm2)缺省： 此项为可选择关键字，无缺省。条件： 此项关键字必须位于组份性质段，当油的压缩系数是压力和泡点压力两者的函数是，则需要使用 *COT 关键字。它对应于油的体积系数，其在泡点压力以上是依赖于压力的非线性函数。在使用 *COT 关键字时，不能同时使用 *CO 关键字，也不能在 PVT 表中输入 co 值。说明： 油的压缩系数 co，给出形式为：在泡点压力以上时，油的体积系数 bo 曲线的斜率与 bo(Pb) 的比值。对于一个典型的黑油系统，当压力 P 大于 Pb 时假设 co 为常数。当 co 相对于 Pb 变化时，显示只具有很小的组份依赖性；而在 P 大于 Pb 时，co 对压力 P 的变化则具有非线性特征。假如将 co 作为常数，那么: bo(P) = bo(Pb) * 1 - co *(P-Pb) 如果 co 只是泡点压力的函数，则对于每一个泡点压力必须使用不同的 co 值。如果同时 co 又是压力的函数，则 co 汇集于 Pb 和 P 之间的压力范围内，在 IMEX 中采用数值方法处理这一问题。 在组份性质数据段中，允许出现多个 *COT 表。如果定义了多个 PVT 表，则 *COT 关键字后必须有一个整数以确定它所代表的 PVT 表号。如果只有一个 PVT 表，则没有必要设置表号。 每个 *COT 表至少要有两行数据，如果计算的网格压力在表数据范围之外，则使用标的最后一行计算网格压力。 要求最少的 *COT 关键字是一个，如果有多个 *COT 表输入，在计算网格的泡点压力时，则采用线性插值方法确定 co 相对于 P 的改变量。如果计算网格的泡点压力落在表的范围之外，则使用最接近于计算值的表。换句话说，就是不在确定的