WinProp相态程序 Windows版本用户
WinProp相态程序 Windows
WinProp .............................................................................................................................................. 1
第一章 数据集结构和控制 ......................................................................................................... 1
第二章 组分 ................................................................................................................................ 4
第三章 常用选择 ........................................................................................................................ 6
第四章 两相饱和度和相边界计算 ............................................................................................. 7
第五章 闪蒸计算 ...................................................................................................................... 10
第六章 三相边界计算 .............................................................................................................. 12
第七章 组分的劈分和归并 ....................................................................................................... 13
第八章 实验室计算 .................................................................................................................. 14
第九章 回归 .............................................................................................................................. 17
什么是WinProp
WinProp是CMGProp的Windows
WinProp是CMG状态方程多相平衡特性软件包,它包括流体特征化、组分归并、实验室数据回归拟合、相态图计算、沥青沉淀等等。在WinProp中考虑的实验包括分离器油和气的合
可用WinProp分析油藏油气系统相态,产生CMG组分模拟器GEM TM WinProp具有图形接口,图形由Excel
WinProp也可产生CMGprop的关键字数据文件,除CMGProp需要的常规关键字外,还包括图形接口的特殊控制字符。因此,WinProp不能解释用编辑程序(editor)产生或修改的数据。
第一章
第一节
数据输入表可以由双击某一行(或一方框)激活。在每一张表格中可以输入和编辑数
若插入一计算选择,则(1)选择将要插入新的计算的某一行,单击该行。(2)由菜单
在Stat(状态)列中的“U”表示相应的数据未定义。一旦表格数据输入,“U”则由状
第二节
1
除第一节描述的插入计算选择以外,也可由菜单选择Edit
undo
Cut 剪切选择,与它们有关的数据进入Clipboard,但对Title和Control及Component
Copy 复制选择。与它们有关的数据拷入CLipboard
Paste 该选择从Clipboard
可使用Copy和Paste选择在一WinProp内传递信息和计算选择或从一个数据集传入另一
在很多表格上,可将数据输入表(或网格),大多数可使用菜单项Table插入一新行或
用左鼠标键点黑不同
中的数据,然后可删除,剪切或粘贴相应单元的数据:删除DELETE,剪切CTRL+x,粘贴CTRL+V,注意:DELELE删除单元中的数据,但并不删除相应的行。如要删除行,使用Table,Delete rows
在编辑一网格单元之后,必须按一Enter键,使得修改有效。对与网格无关的文本箱而言,Enter
第三节 创建、打开和存贮数据文件
选择File,New创建一新的数据文件。WinProp在数据集中预设了三个未定义的表:Title and controls,Component和Composition。现存的文件可由选择File,open„来打开。
存放数据文件用File,save„。存放不同名称时用File,Save as
按常规,所有数据集名称具有?dat后缀。在WinProp 文件后缀
?out
?gem 组成性质输出,
适用于组分模拟器GEM
?rls
?Srf TM ?xls Excel文件,包含由?srf
结果文件在?out文件中,当进行回归、归并和劈分过程时,WinProp产生一?rls文件。可打开该文件产生一新的数据集(参见组分劈分、归并及回归)。如果在表Component properties print option表中,选择用于CMG组分模拟器GEM的组分打印,组分性质则输出到?Gem文件中,格式直接适用于GEM
选File,Exit退出WinProp
第四节 运行、查看输出和产生图形
2
运行数据集,选File,Run(或使用F2
硬拷贝输出在?out文件中。查看硬拷贝输出,选File,View output(或使用F3键)。缺省编辑是DOS编辑器edit。通过在文件CMGLCH.INI 中包含“user_editor=”的重定义行,可以再指定一编辑器。请参见Default user_defined editor. TM 某些计算选择可产生绘图数据。绘图信息在?srf文件中。为用Excel产生图形,选TMFile,Create Excel plots(或键入F4)。WinProp具有它自己的Excel模板和宏,调用TMTMExcel,并产生相应的图形。Excel图形存在.xls文件中,如果没有绘图数据,则输出一信
TM 查看已经用WinProp产生的Excel图形,选File,View Excel plots
第五节 获得帮助信息
选help
Contents
Search for help on
help on current form
当前表格帮助信息也可由F1
第六节
Comments
输入关于该数据集的注释。这些注释将在数据集 Title1,Title2,title3
可输入最多3
Equation of state
选择油气相的状态方程。缺省为PR(1978
PR(1978) Peng-Ronbinson状态方程,
常数“a”的1978
PR(1976) Peng-Ronbinson状态方程,常数“a”的1976表达式,这是原始状态方程。
SRK(G&D) Soave-Redlich-Kwong状态方程,带常数“a”,由Grabosk和Daubert(1978
SRK 原始Soave-Redlich-Kwong
Units
psia & deg F 全部压力用psia,温度用degF
KPa & deg C 全部压力用kPa,温度用deg C 22kg/cm & deg C 全部压力用kg/cm,温度用deg C
Feed
3
Mole 表格Composition
Mass 表格Composition用质量(例:kg
当选Mass时,WinProp使用组分分子量,将所有质量分数转换成摩尔分数。全部输出为相应摩尔分数,而不是输入的质量分数。
第二章 组分
第一节
状态方程需要每个组分的下列性质:临界压力(P)、临界温度(T)、偏心因子(ω)cc
和不同组分之间的相互影响系数(δ)。还需要分子量计算质量密度。对每个组分还要定义ij
附加因子,例如,体积偏移量τ、参数Ωa和Ωb,以增加状态方程预测精度。τ的缺省值是零。Peng-Ronbinson状态方程的Ωa和Ωb的缺省值分别为0.45723553和0.077796074,对Soave-Redlich-Kwong状态方程而言,Ωa和Ωb的缺省值分别为0.4274802和0.08664035。
选择或编辑组分时,在数据集结构表的“Component”行双击,随后出现组分定义(component definition
可从WinProp
提取库组分时选Option,Insert Library component
一、Library component
同时按下CTRL键和鼠标
按下View selection order
如果有一已知性质的组分,可选option,insert own component,Criticalproperties
三、已知SG,Tb和MW
对很多重烃分数而言,已知性质为比重(SG)、常压沸点(Tb)和分子量(MW)。这些组分可选择Option,insert own component ,SG,TB,MW
至少要求其中两种性质。为估计状态方程要求的临界性质和偏心因子,WinProp有三种选择方式:Lee-Kesler相关式(Lee和Kesler,1975),Riazi和Daubert相关式(Riazi和Daubert,1980),或Twu相关式(Twu,1984)
第二节
在组分选择和定义之后,其参数和性质在表格Component definition中表示出来。
表 头 参数或性质
Component 组分名
HC 烃标志(,1,表示烃)
4
P(atm) 临界压力(大气压) c
T(K) 临界温度(K) c
Acentric fact. 偏心因子
Mol.Weight 分子量
Vol.shift 体积偏移(无量纲)
Vc(1/mol) 临界体积(1,mol)
Vc(Viscosity) 用于粘度计算的临界体积(1/mol)
Omega A Ωa因子
Omega B Ωb因子
SG 比重(水,1)
Tb(deg F)
上述表中所列参数在表Component
SG和Tb之值在其定义的组分中表示。一旦SG和T曾用于计算临界性质和偏心因子,则b
改变表中的SG和T值不会影响其它性质。如果希望用修改过的SG和Tb值重新计算给定组分b
的性质,则需使用option,Delete Component从表中删除那个组分。再使用option,Insert own component,SG,TB,and MW,用修改过的SG和T b
第三节
引入相互影响系数(δij
测的饱和压力与实验数据拟合而得。组分之间的相互影响系数在表Int.coef.
一、烃类之间相互影响系
烃组分(HC)由表component中HC列的值为1来标识。烃类之间的相互影响系数由下式计
θ1/61/61/31/3 δ=[(2VV)/(V+V)] ijcicjcicj
其中,Vci是组分i的临界体积,θ
已经证明,θ,1.2能很好地拟合链烷烃之间的相互影响系数(oellrich等,1981)。然而,θ
为避免相互影响系数表的混乱,当装载该表时,未显示烃类的相互影响系数。显示按show HC lnt.Coef.隐含按Hide HC lnt. Coef
WinProp库中的非烃类之间、烃类和非烃类之间的相互影响系数在表中显示出来。也可按需要进行编辑,注意:由于δ=δ ijji
油气之间的粘度由Jossi,stiel和Thodos(JST)相关式(Reid等,1977)计算而得:
*-40.252345 〔(μ,μ)?+10〕,a+aρ+aρ+aρ+aρ+aρ 01r2r3r4r5r
μ ,油气粘度,厘泊或mPa?s *μ,低压粘度,厘泊或mPa?s ,1,21,6,2,3?,T MP, cc
5
T单位为K,P cc
M
ρ,摩尔密度,ρ,ρ,V,V rcc
混合临界体积V c
1/4nc,,4Vc,XV,ici,,i,,,1
其中α是混合指数参数,X是组分组成,Vcomponentici
中的Vc(viscosity
α,a,a,a,a和a,在表Viscosity中输入。当表第一次激活时,缺省值即显示出o1234
第四节
组分性质通常不在输出文件中列出。当需要输出组分性质时,选择Option,print Component properties.在数据集中出现Component properties print Options.可选择详细的组分性质输出或以适用于CMG的GEM格式输出。其输出的组分性质形成(gem文件。
第三章
这一章描述大多数计算使用的常用选择。包括Composition specificaton,Initial K-
values,Output level和Stability test level。大多数实验室计算常用的饱和压力计算在下一章讨论。
第一节
组成以摩尔或重量单位输入,指定为分数或百分数。值在内部归一化。如果输入重量分数或重量百分数,则在内部将它们转换为摩尔分数。重量单位的选择在Titles and controls。表Composition有两列组分组成输入。主组成一般指油藏中油或气组成。主组成必须输入。
在一数据集中,可定义几个组成部分。该组成之后的计算使用这一组成,直到运行结束
?前次计算的n
组成在Feed
6
第二节 初始K
要求初始K
?由Wilson方程(internal
lnK,5(37(1+W)(1-T,T)+ln(P,P iicici
?由前次两相计算(previous
?由前次多相计算的n相(phase
第三节
通常输出级为1。如果要求每一步计算的更多信息,输出级定为2
第四节
在相态计算中,相的数目一般预先是未知的。WinProp假定系统初始为单相,并进行稳定性测试。稳定性测试指确定一个系统是否需要劈分另外的相而获得稳定性。稳定性测试在多维Gibbs自由能表面寻找平稳点。对一稳定相而言,Gibbs自由能必须小于所有平稳点的值。稳定性测试确定平稳点搜寻的彻底性。其值从0到4。通常1值对大多数两相油气系统而言已经足够。对多于两相的系统而言,可能需要其值等于4。当一个数据表首次激活时,稳定性测试级置为缺省值。如果怀疑系统可能具有更多的相,增加其值,重新运行数据集。
第四章
?相态图计算(压力—温度、压力—组成、温度—
相边界计算还生成相的等摩尔分数线或等体积分数线(特征线)
第一节
7
选Calculation,2P Saturation Pressure调用该选择。有一数据集例子为pressat.dat.
Feed,K-values,Output和Stability test level
如果估计值误差很大,可以请求WinProp内部产生一饱和压力计算的更好的初始假设,只需点一下Improve saturation pressure estimate
与饱和压力计算有关的实验数据(它们可以放入一回归中拟合),在表Experimental中表示。包括饱和压力、液体和蒸汽质量密度、压缩系数和粘度,也可指定对每个实验数据的
第二节
选:Calculation,2P Saturation temperature.
数据例:tempsat.dat.
Feed,K-values,Output leveL和stability test level,
饱和温度的估计在合适的文本箱输入,如果估计很差,则可请求WinProp产生一较好的
与温度压力计算有关的实验数据在表Experimental中,包括饱和温度、液体和蒸汽质量
第三节
选Calculation ,2P Envelopeenvel 2P1.dat(PT图)和envel2P2.dat (PX
两相包络图计算产生单相和两相区的边界。泡点包络图指单相液区和两相气液区的边
?压力—温度(PT
?压力—组成(PX
?温度—组成(TX
在表Type中可指定计算的包络图的类型。在这张表中,可指定相应X轴和Y轴的变量。还可输入一X变量(当前为温度)和Y变量的估计值,从而指定要计算的第一个点。X变量和Y变量分别称为独立或依赖型变量,Y变量的输入选择为“unknown”,“previous”,或一输入值,也可输入X和Y的最大最小值。当达到最大或最小值时或计算到最大点数时,计算停止(参看Vapor,Control
在计算相态图过程中,WinProp还通过插值计算临界点的位置。如果临界点位于相包络图
二、Vapor/control
1
8
蒸汽摩尔分数为零,计算泡点包络图,为1计算露点包络图。为方便起见,WinProp一般
2
蒸汽体积分数值为0,1,计算与F和1—F相对应的曲线,两条曲线在临界点混合。 vv
当输入一蒸汽体积分数值时,WinProp计算与V和1—V相对应的曲线,等体积分数曲线vv
3
4
初始步长控制包络图计算点之间的间隔。可使用正值和负值。正值表示相态图初始轨迹为X值增加方向,负值则为X值减少方向。WinProp
5
指要计算Y值的X值。因为包络图计算中的步长是自动的,因此,输入这些值使得在这些点上计算。
第四节
临界凝析压力指PT相包络图上的最大压力,而临界凝析温度指最大温度。这两个值在两相PT包络图计算中估计(参看相边界和特征线计算)或直接选Calculation |Criconden-
bar/therm计算。例子为cricon.dat
对Feed,K-values,Output level和Stability test level见常用选择。压力和温度的初值,可选unknown或Previous
第五节
相边界和特征线计算通过插值计算临界点。该方法是有效的,产生相边界和临界点。然而,如果想直接计算或想在回归计算中拟合临界点,应当使用Critical points计算选择。临界点计算由选择Calculation,Critical points调用。例子为Critical.dat,该选择使用Heidemann和Khelil(1980)计算方法。需要的数值数据包括无量纲体积下限和无量纲体积上限。无量纲体积等于摩尔体积V与状态方程参数b之比。计算两个体积限之间的全部临界点。缺省值分别为1.0和5.0。可以输入实验临界压力和温度,用于回归目的。
第六节
气和溶解剂注入过程通常包括多次接触计算,研究蒸发或萃取过程。在结果解释时产生拟三角相图。大多数气或溶解剂注入是汽化,凝析过程(Zick,1986),结果三角相图呈砂漏
9
多次接触计算选Calculation,Multiple contacts来调用,例子为multi .dat.计算采用下述步骤,压力和温度在表conditions
1
解剂。主气和补气的摩尔组成,补气的分数在表composition中输入,油的组成在表Feed中
2
量比在表condition的solvent increment ratio输入。缺省值为0.05。闪蒸计算最多100次混合。如果未检测到两相区,则认为是一次接触混相,计算停止。在两相区形成后,计算继续第3步和第4
32步检测的两相区中的第一个点,取掉所有液体。剩余气与原始油以气油比Ratio:(1—Ratio)结合,其中Ratio的值在Conditions中的Equilibrium gas/original
oil mixing ratio 输入(缺省值为0.01),执行闪蒸计算,取掉液体,过程重复进行。该过程模拟蒸发或萃取过程。最多进行50
42步检测到的两相区的第一个点,取掉所有气体。剩余的液体与原始溶解剂以溶解剂—液体比RATIO?(1—RATIO)结合。进行闪蒸计算,取掉蒸汽。过程继续,直到油不能进一步变富,或执行50
表composition的最后一列确定作三角相图三个拟组分的分组。
第五章
闪蒸计算决定系统在给定压力、温度和组成下的劈分。计算相的数目和每个相的性质。WinProp
1—
2—液体1—液体2
3—液—
4
5
在上述计算中,流体相用状态方程模拟,除第3个以外。在第3个计算中,水相中的组分溶解度用亨利定律模拟。计算5
第一节
调用方法为calculation,2P flash,实例为flash 2P dat.
Feed,K-values,Output level和stability test level
闪蒸计算在给定的压力和温度下进行。通过指定压力或温度步及步数执行一系列计算。
回归计算中的实验数据在表Experimental和Exp.K-values中输入。实验数据包括质量密度、摩尔分数、体积分数、压缩系数、气相和液相的粘度。实验K—值在表Exp.K-values中
10
第二节
调用方法:Calculation,Multiphase flash
Feed,K-values,output level和Stability test level,
在flash type中选择多相计算的类型。三相和四相计算使用的计算技术描述见Nghiem和Li(1984)。这是一个有戏剧效果的过程,其相的数目逐步增加。全部相用状态方程模拟。选择的相的数目指最大相数。因此,对两相系统选四相计算,其结果与两相闪蒸计算相同。
OGW计算包括三相计算,其中气液相用状态方程模拟,水相用Henry定律模拟。由于状态方程是为类气烃系统而推导的,因而不可能精确模拟水相。Li和Nghiem(1986)推荐,在水相中对组分溶解度使用Henry定律常数。水相中组分i的逸度系数Φ iw
lnΦ=ln(Hi/p iw
其中Hi是组分i的Henry定律常数,每个组分的Hi在表Henry′s law中输入。如果不指定Hi,WinProp
三相和OGW计算的实验数据在表Experimental中输入,包括质量密度、摩尔分数、体积
WinProp
第三节 具有固相的多相闪蒸计算
1
模拟沥青沉淀的模型详见Nghiem等(1993)。表示沉淀沥青为单组分固相。固相中沥青
** lnfs=lnfs+Vs(P,P),(RT * 其中fs和fs是纯固体沥青在压力P和P*下的逸度,Vs为纯沥青的摩尔体积,R是通用**气体常数,T为温度,fs,P和Vs
2
模拟沥青沉淀的关键步骤是沥青组分的特征化,包括溶解状态和固相。已经发现,将该
调用方法Calculation,SLV asphaltene flash,实例为asphalt.dat.在表Calculation中的数据输入与其它闪蒸计算类似。在表Ref.state
1
三相闪蒸算法交替进行闪蒸和稳定性计算,以对单相流体做稳定性测试作为开始。如果出现不稳定,则进行两相闪蒸计算,紧接着对两相系统稳定性进行测试等等。三种计算方法均可用于计算,其差别在于稳定性测试执行的顺序不同。在方法1中,其稳定性测试首先对固相。在方法2(缺省)中,稳定性测试首先对全部流体相,然后对固相。方法3是方法1的特例;方法3对收敛的两相流体—固体系统不做稳定性测试。因此,方法3更有效,但不如方法1严密。对大多数情况而言,方法1和方法2收敛结果相同。也有例外情况,一种方法
11
2ln(固相逸度(atm)
Calculate:在给定压力下的参考逸度等于系统收敛于单液相或气液两相后沥青组分的逸
previous
input value
3
计算参考逸度的参考压力。仅当输入参考逸度时才需输入该值。当选Calculate时,参
41,mol
?参考逸度选择为Calculate,
?参考逸度选择为Previous
第四节
单相区的闪蒸计算生成单相系统。如果预先知道流体为单相,则其性质可用单相计算直接计算而得。调用方法为Calculation,Single-phase fluidWinProp假定全部计算为单相。实例为:Single.dat。
第六章 三相边界计算
对于显示三相相态的系统而言,存在条件使三相之一的摩尔分数趋于零。在此条件下,存在两相与一无限小量的第三相之间的平衡。所有这些条件的轨迹为三相边界。Nghiem和Li(1984)描述了构造三相边界的计算技术;这些是两相边界计算的扩展。可计算下述包络图。
?压力—温度(PT
?压力—组成(PX
?温度—组成(TX
第一节
调用方法Calculation,3P Envelope,例子为Envel,3P1.dat和enuel3P2.dat要求给予恰当的初值以使其收敛。因此,压力、温度和K—值必须由边界附近的前次三相闪蒸
通过选择X轴和Y轴的变量,在这张表上可指定要计算的PT、PX或TX包络图的类型。还可输入一X变量值作为要计算的第一个点并估计一Y变量值。X变量和Y变量还分别称为独立和非独立变量。Y变量的选择为“previous”,或输入一值。X变量和Y变量的最大和最小12
值在适当的文本格输入。当最大最小值达到或最大计算点数达到,计算停止(参看相,控
在该表中,可输入气摩尔分数或第二液相摩尔分数的初值。输入的气摩尔分数值作为与第二液相零摩尔分数相对应的边界计算的初值。相反,输入的第二液相摩尔分数作为相应于零气相边界计算的初值。选择use values from previous calculations,WinProp计算相应
1
2
初始步长控制相包络图上计算点之间的间隔。可使用正值和负值。正值的相图初始轨迹沿X值增加方向。负值的相图初始轨迹沿X值减小方向。WinProp内部估计随后点的步长。
3
该值调整相图上两相邻点之间的距离。如果实际迭代次数小于输入值,距离增加;反之
这些表值指想要计算Y值的X值。因为在相图计算中步长自动,这些输入值迫使计算在
四、K—
边界图上第一个点的K值估计值可以取自前次计算,或以上述表格输入。
第七章
经验证明,模拟油藏流体相态通常需要两个阶段。第一阶段,流体系统由大量组分表示。(例如:C1,C2,C3,„C29C30+
简单计算以验证状态方程的正确性。我们发现,在大多数情况下,状态方程能够准确预测饱和压力,而只需稍微调整一下烃类相互影响系数指数(HICE)。然而,多组分模型对组分模拟并不实用,因为它将耗费额外的运行时间和存贮空间。第二阶段,包括将多组分系统归并为几个组分(例如:10个)。我们建议在油藏主要条件下,基于多组分系统的K值归并组
1CSCN 7+
2
当需要详细分析C1。油藏流体典型地由纯净的、意义明确的组分,7+
比如,CO2,N2,C1,C2,等等和上百种庚烷和重组分(C)组成。油藏流体实验室分析一般7+
包括将C7+馏分分离为单烃组分的气色谱分析,例如,直到C30+
为油藏流体的多组分表示。如果实验室没有分析,则需调用劈分选择将重组分劈分。劈分后,将单烃组分基于Wilson相关式产生的K
在第二步,调用归并计算减少组分数目。正如前面所讨论的,建议依据在油藏条件下状态方程的K
第一节 劈分“plus
13
调用方法Lump,split,splitting,实例为splitdat.
在表properties输入plus
在表controls
一、plus
指WinProp将plus
有下面几
No lumping 单烃组分不做归并
Determined internally WinProp内部估计plus馏分的拟组分数
input value
1Lee-Kesler(1975
2Riazi-Daubert(1980)
3Twu(1984)
第二节
调用方法为Lump/split,lumping.实例为lump.dat.
第一列是多组分系统的组分名,主组成和次组成列在第二列和第三列。在第四列按鼠标左键指定归并的设置。在第i行第一次按左鼠标键将多组分系统的组分1到组分i归并为拟组分1。第二次在第j行将多组分系统的组分i+1到组分j归为拟组分2,依此类推。undo取消前一次鼠标操作。Reset
第三节 将结果文件转到其它数据集
劈分和归并计算结果是随后计算中将要使用的组分性质和组成。当做这些计算时,比如对test.dat,除产生
输出test.out外,WinProp还产生一输出文件test.rls?rls文件含有组分性质和组成,可传输到其它数据集。这可以通过调用File,open„„,选用test.rls。该文件由三部分组成:Title and controls,components和composition。计算选择可由菜单插入,在定义了计算选择之后,使用File,Save as„„将该文件存为另一数据组,例如:test1.dat。类似的过程适用于Regression选择。
第八章 实验室计算
WinProp
1
2
3
14
4
5
6
7
实验数据以表格形式输入。一旦在某一列输入一些数据,要求在那一列的所有行均输入数据。对实验值没有的行,在相应的位置输入“,1.0”,在回归中程序忽略全部负值。WinProp对某些行有数据的列的空行自动赋“,1.0
也可以在任何一个数据前赋一“,”号,而将该数据在
第一节
当无法获取井底流体样品时,就将分离器油气混合用来代表油藏流体。该过程按分离器
选Lab,Recombination
分离器压力、温度和汽油比在Recombination表中指定。分离器油气组成在表Composition
该选择依据油气各自组成,分离器压力和温度,计算油密度和气密度。从状态方程角度来看,气油并不处于平衡状态。因此,对混合流体的分离器试验计算并不给出测量的汽油比
混合的另一方法是使用Separator test计算,通过尝试法,找出气油混合物中气的部分。气油组成在composition
第二节
选择Lab,compressibility calculation来调用。例子为compressdat。压力和温度数据在表calculation
可以输入原始压力、压力增量、混合物中气的摩尔分数和压力级数,实验压缩系数在表Experimental
第三节
油藏流体样品放在实验室容器中。调整压力等于或大于原始油藏压力。温度设置为油藏
在每一步,测量压力和油藏流体(油和气)的总体积。该过程也称为闪蒸、闪蒸释放、
调用方法Lab,constant Composition expansion.例子为cce.dat.在表Sat.pres和Pres levels
等组分膨胀的压力级和相对油体积的实验数据在表Pres.Levels中输入。相对油体积(ROV
15
第四节
一般对黑油流体做该实验,模拟其油藏条件。实验室中油藏流体样品加压到泡点压力,
增加容器体积,压力减少。去掉全部气体,同时减少容器体积,保持压力不变。收集排
随后,温度减少为15?或60?F,测量剩余液体体积,这就是微分脱气的剩余油。
每次容器中液体体积除以剩余油体积就是相对油体积或地层体积系数B。气的压缩因子o(Z因子)、气密度和地层体积系数Bg也可测量而得。整个过程中去掉的气的总体积是在泡点压力下溶解在油中的气量。溶解气油比(R)等于总的气体积除以剩余油体积。任意一个较s
调用方法Lab,Differential liberation.例子为:diflib.dat。输入数据在表Sat.Pres,Std.Conditions,和Pres.Levels中.
饱和压力数据输入说明见Saturation pressure o标准条件在表Std.condition中输入。缺省值为14.696psia和60F或101.325kPa和15
在表PresLevels,输入微分脱气压力数据,实验数据,油比重(SG),油的地层体积系数(FVF),溶解气油比,气的Z因子,气的地层体积系数,气的比重,油和气的粘度以及
第五节
该实验通常模拟凝析气的油藏条件。实验室容器中的油藏流体样品加压到露点压力,温
增加容器体积,减小压力。去掉一部分气体,直到容器体积等于露点时体积。收集气送
调用方法Lab,constant Volume depletion,实例为cvd.dat.在表sat, pres.Separator,
和Pres.Levels
多级分离器条件和标准条件在表separator中输入。多级分离器并不包括Standard conditions中输入的储罐条件。也可指定较低碳数目和较高碳数目组分。对每一压力级计算
在表Pres.Levels中输入压力级,也可输入实验数据:产出气的摩尔百分数,析出液体量和气体偏差因子,用于回归计算。
第六节
16
油藏流体送入满足一定压力和温度条件的多级分离器,使得液体产量最大,分离器的最
测量总的气油比,储罐比重或API
二、输入
调用方法Lab,Separator test,例子为sep.dat,在表Sat.Pres.,Separators和Experimental
分离器条件在表Separators
不同级的压力和温度在表的左边输入,不包括储罐条件,储罐条件在右边输入。如果除
第七节
该实验可提供注气过程中关于流体相态的信息。当气注入油藏中时,气体溶解于油中,
油藏油放入容器中,温度设定为油藏温度。测量油的泡点和相应体积。注入少量气体。确定新的饱和压力,记录新的饱和体积。该过程重复数步,直到达到注入气浓度的上界(例如:80
调用方法Lab,Swelling Test,例子为swell.dat,在表Sat.Pres, Pres.Levels 和CCE
在表Mixtures中,第一列指注入气的摩尔分数,第二列指注入气和油混合物的饱和压力估计值。第三列和第四列也可输入实验饱和压力和膨胀
如果必要,在表CCE中输入等组分膨胀过程的压力级数据。
第九章
可使用WinProp的回归特征拟合实验数据。在回归中使用的参数是组分性质和相互影响系数。WinProp
WinProp使用Agarwal等(1990)的回归过程。对用户输入的参数,模型将参数排序,使得最敏感参数位居前列。每次回归少量参数,一旦达到参数允许的最大值或最小值,或者改变参数值拟合无任何改进,则按参数排序拟合下一个参数。因此,用户可以指定大量参数,WinProp
第一节
注意Form中的Regression Parameters和End regression参数。两种表之间的计算为回归实验数据,表Regression Parameters用Regression,Start插入,而表End Regression用Regression,End来调用。在每个计算选择的适当表格输入实验数据。权与每一个数据或数据组有关。权W给每个数据分配一重要度。缺省值是1.0。大值表示该数据重要i
性较大,而小值则较小。回归使得下列目标
为最小:
17
2F,Wi(X,X),icalcimeas,,,i1
其中XXi,calci,meas
在Regression Parameters之前的全部计算使用回归之前的参数。在End Regression之后的全部计算使用回归后的参数。对回归中包含的计算,输出回归前和回归后的结果和图
第二节
在表Comp.Prop
在网格上按鼠标键选择或取消组分性质。第一次按键选择该性质,第二次按键取消该性质。
在表Int.Coef
烃组分之间的相互影响系数由第二章第三节描述的方程计算。HC int.Coef.Exponent
在表Num.Params
Convergence tolerance
两次迭代之间的目标函数小于该值,回归
Maximum no of iterations
Noof Simultaneous reg.Parameters
每次回归的参数个数。缺省值为5。如果回归参数个数小于5
第三节
回归结果是随后计算中使用的组分性质。比如,对regreess.dat进行回归,除标准输出regress.out外,WinProp还输出一文件regress.rls。?rls文件包括回归后的组分性质,其格式可直接用于另一数据集。调用方法File,open„,打开文件regress.rls 。该文件由三部分组成:Title and controls,Component 和Composition,然后可插入计算。使用File,Saveas„将该文件存入另一数据集,例如test1.dat。类似的做法适用于Component splitting and lumping。
18