CMG数模软件STARS模块使用入门教程.pptx

CMGCMG数模软件培训数模软件培训庞占喜庞占喜2007.3.17 中国石油大学（北京）中国石油大学（北京）中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 目目 录录*CMGCMG软件简介软件简介*STARSSTARS模块主要关键字模块主要关键字*STARSSTARS模块泡沫的模拟模块泡沫的模拟*STARSSTARS模块所需数据的准备及处理模块所需数据的准备及处理*STARSSTARS模块油藏热采模型的建立模块油藏热采模型的建立*油藏热采模型的运行及结果后处理油藏热采模型的运行及结果后处理*氮气及氮气泡沫压水锥数值模拟氮气及氮气泡沫压水锥数值模拟中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 一、一、CMGCMG软件简介软件简介STARSSTARSIMEXIMEXGEMGEMMODEL MODEL BUILDERBUILDERGRID GRID BUILDERBUILDER3D3D2D2D CMG可可以以进进行行常常规规黑黑油油模模拟拟、稠稠油油热热采采模模拟拟、组组分分模模拟拟以以及及泡泡沫沫模拟。模拟。STARS模模块块是是三三维维、四四相相、多多组组分分、热热采采、蒸蒸汽汽添添加加剂剂模拟器。模拟器。其其数数据据体体文文件件 为为*.dat.dat,计计算算产产生生的的文文件件包包括括:输输出出文文件件*.out.out(用用 户户查查看看),SR2二二进进制制索索引引文文件件*.irf.irf(数数据据 后后 处处 理理),SR2二二进进制制结结果果文文件件*.mrf.mrf(二二 进进制制计计算算结结果果存储存储)。中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 油藏模型数据体包含内容油藏模型数据体包含内容STARS数据体数据体INPUT/OUTPUT CONTROLINPUT/OUTPUT CONTROL：输入输入输入输入/输入控制输入控制输入控制输入控制，定义控制模拟器输入和输，定义控制模拟器输入和输出行为的各个参数，例如，文件名、单位、出行为的各个参数，例如，文件名、单位、outout文件和文件和SR2SR2文件写入频率，文件写入频率，重启文件的定义等。重启文件的定义等。GRID AND RESERVOIR DEFINITIONGRID AND RESERVOIR DEFINITION：网格和油藏定义网格和油藏定义网格和油藏定义网格和油藏定义，这部分包括：，这部分包括：模拟网格的定义、天然裂缝油藏选项、离散化井筒定义、基本油层岩石特模拟网格的定义、天然裂缝油藏选项、离散化井筒定义、基本油层岩石特性、区块选项，其他油藏特性描述（岩石压缩系数、岩石热物性参数、顶性、区块选项，其他油藏特性描述（岩石压缩系数、岩石热物性参数、顶底盖层热损失系数、井筒热损失系数、水体）。底盖层热损失系数、井筒热损失系数、水体）。FLUID AND COMPONENT DEFINITIONSFLUID AND COMPONENT DEFINITIONS：流体和组分定义流体和组分定义流体和组分定义流体和组分定义，定义组分，定义组分名称、个数，相应的名称、个数，相应的K K值，各组分的基本参数（摩尔质量、密度、粘度、临值，各组分的基本参数（摩尔质量、密度、粘度、临界温度、临界压力，化学反应式等）。界温度、临界压力，化学反应式等）。一、一、CMGCMG软件简介软件简介中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 油藏模型数据体包含内容油藏模型数据体包含内容STARS数据体数据体ROCK-FLUID PROPERTIESROCK-FLUID PROPERTIES：岩石岩石岩石岩石-流体特性流体特性流体特性流体特性，定义相渗曲线，毛管压力、，定义相渗曲线，毛管压力、组分的吸附和扩散特性；（组分的吸附和扩散特性；（*泡沫的定义以及相渗插值的定义）。泡沫的定义以及相渗插值的定义）。INITIAL CONDITIONSINITIAL CONDITIONS：初始条件初始条件初始条件初始条件，这部分包括：初始压力（或者参考压，这部分包括初始压力（或者参考压力及参考深度），初始温度，初始的饱和度场（或者油水界面及油气界面）力及参考深度），初始温度，初始的饱和度场（或者油水界面及油气界面）WELL AND RECURRENT DATAWELL AND RECURRENT DATA：井定义和生产动态数据井定义和生产动态数据井定义和生产动态数据井定义和生产动态数据，定义井名、井，定义井名、井位和完井层位，设置相对应的生产动态数据。位和完井层位，设置相对应的生产动态数据。NUMERICAL METHODS CONTROLNUMERICAL METHODS CONTROL：数值方法控制数值方法控制数值方法控制数值方法控制，这部分定义模拟，这部分定义模拟器数值方法参数：时间步数、非线性迭代解法、误差控制；（器数值方法参数：时间步数、非线性迭代解法、误差控制；（*等温、非等等温、非等温控制项）温控制项）一、一、CMGCMG软件简介软件简介中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 INPUT/OUTPUT CONTROLINPUT/OUTPUT CONTROL输入输入/输出文件名输出文件名(可选可选)*FILENAMES(*FILENAMES(建立重启文件用建立重启文件用建立重启文件用建立重启文件用)错误检查的模式错误检查的模式(可选可选)*CHECKONLY项目的主要名称项目的主要名称(可选可选)*TITLE1,*TITLE2,*TITLE3,*CASEID*TITLE1,*TITLE2,*TITLE3,*CASEID输入输入/输出单位控制输出单位控制(可选可选)*INUNIT,*OUTUNIT*INUNIT,*OUTUNIT质量守恒选项质量守恒选项(可选可选)*MASSBASIS错误信息的最大数量错误信息的最大数量(可选可选)*MAXERROR重启文件选项和最大计算步数重启文件选项和最大计算步数*RESTARTRESTART,*MAXSTEPS写入重启文件选项写入重启文件选项(可选可选)*WRST,*REWIND*WRST,*REWINDOUT文件写入频率文件写入频率(可选可选)*WPRN写入写入OUT文件中的信息选项文件中的信息选项(可选可选)*OUTPRN*OUTPRN,*PARTCLSIZESR2文件写入频率文件写入频率(可选可选)*WSRF写入写入SR2文件中的信息选项文件中的信息选项(可选可选)*OUTSRF*OUTSRF,*SR2PREC,*SRFASCII,*XDR网格数组输出方向选项网格数组输出方向选项(可选可选)*PRNTORIEN,*PRINT_REF方程组求解信息输出选项方程组求解信息输出选项(可选可选)*OUTSOLVR模拟计算终止方式选项模拟计算终止方式选项(可选可选)*INTERRUPT二、二、STARSSTARS模块主要关键字模块主要关键字中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 GRID AND RESERVOIR DEFINITIONGRID AND RESERVOIR DEFINITION网格类型选项网格类型选项 *GRID,*KDIRGRID,*KDIRI方向网格步长方向网格步长(必需必需)*DI*DI J方向网格步长方向网格步长(必需必需)*DJ*DJK方向网格步长方向网格步长(必需必需)*DK*DK网格中心深度选项网格中心深度选项(条件条件)*DEPTH网格顶部深度数据网格顶部深度数据(条件条件)*DTOP*DTOP油层中部深度数据油层中部深度数据(条件条件)*PAYDEPTH网格倾角数据网格倾角数据(条件条件)*DIP局部网格加密选项局部网格加密选项(条件条件)*REFINE,*RANGE网网格格几几何何特特征征修修正正选选项项(可可选选)*VAMOD,*VATYPE（修修正正部部分分网网格格的的面面积积及体积）及体积）无效网格选项无效网格选项(可选可选)*NULL*NULL离离散散井井筒筒选选项项(条条件件)*WELLBORE,*RELROUGH,*LAMINAR,*TRANSIENT,*CIRCWELL,*WELLINFO,*REGIME,*WELLWALL,*TUBINSUL,*ANNULUSWAL,*CASING,*FILM_COND,*RANGE,*WBZ,*WBZADJ二、二、STARSSTARS模块主要关键字模块主要关键字中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 二、二、STARSSTARS模块主要关键字模块主要关键字GRID AND RESERVOIR DEFINITIONGRID AND RESERVOIR DEFINITION孔隙度数据输入标志孔隙度数据输入标志(必需必需)*POR*POR渗透率数据输入标志渗透率数据输入标志(必需必需)*PERMI,*PERMJ,*PERMK*PERMI,*PERMJ,*PERMK孔隙体积修改因子孔隙体积修改因子(可选可选)*VOLMOD有效厚度有效厚度(可选可选)*NETPAY*NETPAY净毛比净毛比(可选可选)*NETGROSS*NETGROSS传导率因子传导率因子(可选可选)*TRANSI,*TRANSJ,*TRANSK,*TRANSIJ+,*TRANSIJ-,*TRANSIK+,*TRANSIK-尖灭网格输入尖灭网格输入(可选可选)*PINCHOUT,*PINCHOUTARRAY断层选项断层选项(可选可选)*FAULT,*FAULTARRAY门门限限孔孔隙隙体体积积(可可选选)*PVCUTOFF（如如如如果果果果网网网网格格格格的的的的孔孔孔孔隙隙隙隙体体体体积积积积小小小小于于于于此此此此值值值值，此此此此 网格将被设定为无效网格）网格将被设定为无效网格）网格将被设定为无效网格）网格将被设定为无效网格）油层分段选项油层分段选项(可选可选)*SECTOR*SECTOR,*SECTORARRAY网格定义结束标志网格定义结束标志(必需必需)*END-GRID*END-GRID 油藏描述所需资料：油藏描述所需资料：油藏描述所需资料：油藏描述所需资料：小层数据表或等值线图小层数据表或等值线图小层数据表或等值线图小层数据表或等值线图 砂层厚度砂层厚度砂层厚度砂层厚度 有效厚度（或净毛比）有效厚度（或净毛比）有效厚度（或净毛比）有效厚度（或净毛比）顶部深度顶部深度顶部深度顶部深度 孔隙度孔隙度孔隙度孔隙度 渗透率（水平及垂向）渗透率（水平及垂向）渗透率（水平及垂向）渗透率（水平及垂向）断层参数断层参数断层参数断层参数中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 二、二、STARSSTARS模块主要关键字模块主要关键字OTHER RESERVOIR DEFINITIONOTHER RESERVOIR DEFINITION岩石类型岩石类型 *ROCKTYPE,*THTYPE岩石压缩系数岩石压缩系数(必需必需)*PRPOR,*CPOR,*CTPOR*PRPOR,*CPOR,*CTPOR,*CPORPD,*PORMAX岩岩石石热热物物性性参参数数(可可选选)*ROCKCP,*ROCKCP,*THCONR,*THCONR,*THCONW,*THCONW,*THCONO,*THCONO,*THCONG*THCONG,*THCONMIX盖底层热损失参数盖底层热损失参数(可选可选)*HLOSSPROP,*HLOSST*HLOSSPROP,*HLOSST,*HLOSSTDIFF井井筒筒热热损损失失参参数数(可可选选)*RTI,*RTO,*RIN,*RCI,*RCO,*RH,*ETO,*ECI,*EIN,*EE,*XKE,*XKIN,*XKCM,*XAE,*DEPTH,*CDEPTH,*ANG,*AGD水水体体模模型型 *AQUIFER,AQUIFER,*AQGEOM,*AQGEOM,*AQH,*AQH,*AQPERM,*AQPERM,*AQCOMP,*AQPOR,*AQVISC*AQCOMP,*AQPOR,*AQVISC,*AQRCND,*AQRCAP 其他油藏描述参数所需资料：其他油藏描述参数所需资料：其他油藏描述参数所需资料：其他油藏描述参数所需资料：1.1.1.1.油层热物性及流体热物性数据油层热物性及流体热物性数据油层热物性及流体热物性数据油层热物性及流体热物性数据 岩相热容岩相热容岩相热容岩相热容 岩石导热系数岩石导热系数岩石导热系数岩石导热系数 油油油油/气气气气/水导热系数水导热系数水导热系数水导热系数 顶、底层的热损失参数顶、底层的热损失参数顶、底层的热损失参数顶、底层的热损失参数 2.2.2.2.流体和岩石化验分析报告流体和岩石化验分析报告流体和岩石化验分析报告流体和岩石化验分析报告 3.3.3.3.边、底水数据报告边、底水数据报告边、底水数据报告边、底水数据报告中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 FLUID AND COMPONENT DEFINITIONSFLUID AND COMPONENT DEFINITIONS组分类型及名称组分类型及名称(必需必需)*MODEL,*COMPNAME*MODEL,*COMPNAMEK值相关系数值相关系数 *KV1,*KV2,*KV3,*KV4,*KV5K值表值表 *GASLIQKV,*LIQLIQKV,*KVTABLIM,*KVTABLE,*KVKEYCOMP摩尔质量摩尔质量(必需必需)*CMM*CMM临界特性临界特性(必需必需)*TCRIT,*PCRIT*TCRIT,*PCRIT,*IDEALGAS参考条件参考条件 *PRSR,*TEMR,*PSURF,*TSURFPRSR,*TEMR,*PSURF,*TSURF,*SURFLASH流体焓系数流体焓系数 *CPL1,*CPL2,*CPL3,*CPL4,*CPG1,*CPG2,*CPG3,*CPG4,*HVR,*EV,*HVAPR固相特性固相特性(必需必需)*SOLID_DEN,*SOLID_CP*SOLID_DEN,*SOLID_CP二、二、STARSSTARS模块主要关键字模块主要关键字中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 FLUID AND COMPONENT DEFINITIONSFLUID AND COMPONENT DEFINITIONS液液相相标标识识 *LIQPHASE,*WATPHASE,*OILPHASE（随随随随后后后后的的的的特特特特性性性性参参参参数数数数赋赋赋赋值值值值对象）对象）对象）对象）液相密度液相密度(必需必需)*MOLDEN,*MASSDEN,*MOLVOL,*CP,*CT1,*CT2*MOLDEN,*MASSDEN,*MOLVOL,*CP,*CT1,*CT2,*GASSYLIQ液相密度非线性计算参数液相密度非线性计算参数 *DNMIXCOMP,*DNMIXENDP,*DNMIXFUNC气相粘度气相粘度 *AVG,*BVG,*GVISCOR液相粘度液相粘度(必需必需)*AVISC,*BVISC,*VISCTABLE*AVISC,*BVISC,*VISCTABLE,*XNACL液相粘度非线性计算参数液相粘度非线性计算参数 *VSMIXCOMP,*VSMIXENDP,*VSMIXFUNC指定网格的粘度类型指定网格的粘度类型(可选可选)*VISCTYPE,*VSTYPE关键化学反应方程式参数关键化学反应方程式参数 *STOREAC,*STOPROD,*FREQFAC非关键化学反应方程式参数非关键化学反应方程式参数 *RENTH,*RPHASE,*RORDER,*EACT,*O2PP,*O2CONC,*RTEMLOWR,*RTEMUPR二、二、STARSSTARS模块主要关键字模块主要关键字中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 ROCK-FLUID PROPERTIESROCK-FLUID PROPERTIESRock-Fluid数据起始标志数据起始标志(必需必需)*ROCKFLUID*ROCKFLUID岩石类型数据岩石类型数据 *RPT,*KRTYPE,*RTYPERPT,*KRTYPE,*RTYPE插值参数插值参数 *INTCOMPINTCOMP界面张力数据界面张力数据 *INTLIN,*INTLOG,*IFTTABLEINTLIN,*INTLOG,*IFTTABLE泡沫插值参数泡沫插值参数泡沫插值参数泡沫插值参数 *FMSURF,*FMCAP,*FMOIL,*FMGCP,*FMOMF,FMSURF,*FMCAP,*FMOIL,*FMGCP,*FMOMF,*FMMOB,EPSURF,*EPCAP,*EPOIL,*FMMOB,*EPSURF,*EPCAP,*EPOIL,*EPGCP,*EPOMF *EPGCP,*EPOMF二、二、STARSSTARS模块主要关键字模块主要关键字有关泡沫参数的设置及泡沫模拟方法后面介绍。有关泡沫参数的设置及泡沫模拟方法后面介绍。有关泡沫参数的设置及泡沫模拟方法后面介绍。有关泡沫参数的设置及泡沫模拟方法后面介绍。泡沫定义参数所需资料：泡沫定义参数所需资料：泡沫定义参数所需资料：泡沫定义参数所需资料：表面活性剂类型表面活性剂类型表面活性剂类型表面活性剂类型,摩尔质量摩尔质量摩尔质量摩尔质量,密度密度密度密度,临界压临界压临界压临界压 力力力力,临界温度临界温度临界温度临界温度,压缩系数压缩系数压缩系数压缩系数 油藏注表活剂后表面张力随浓度变化数据油藏注表活剂后表面张力随浓度变化数据油藏注表活剂后表面张力随浓度变化数据油藏注表活剂后表面张力随浓度变化数据 表活剂在油层中的吸附数据表活剂在油层中的吸附数据表活剂在油层中的吸附数据表活剂在油层中的吸附数据 表活剂的最大吸附量表活剂的最大吸附量表活剂的最大吸附量表活剂的最大吸附量 表面活性剂在油水两相间的分配数据表面活性剂在油水两相间的分配数据表面活性剂在油水两相间的分配数据表面活性剂在油水两相间的分配数据 泡沫阻力因子和残余阻力因子泡沫阻力因子和残余阻力因子泡沫阻力因子和残余阻力因子泡沫阻力因子和残余阻力因子 中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 ROCK-FLUID PROPERTIESROCK-FLUID PROPERTIES插值序号及参数插值序号及参数 *KRINTERPKRINTERP,*DTRAPW,*DTRAPN*DTRAPW,*DTRAPN,*WCRV,*OCRV,*GCRV,*SCRV水油相渗数据表水油相渗数据表 *SWTSWT气液相渗数据表气液相渗数据表 *SLTSLT相渗临界点相渗临界点 *SWR,*SORW,*SGR,*SORG,*SWRGSWR,*SORW,*SGR,*SORG,*SWRG,*KRWRO,*KROCW,*KRGCW,*KRWRO,*KROCW,*KRGCW,*PCWEND,*PCGEND相渗临界点随时间的变化数据相渗临界点随时间的变化数据 *KRTEMTAB组分弥散系数组分弥散系数 *DISPI_WAT,*DISPJ_WAT,*DISPK_WAT,*DISPI_OIL,*DISPJ_OIL,*DISPK_OIL,*DISPI_GAS,*DISPJ_GAS,*DISPK_GAS吸附组分函数吸附组分函数 *ADSCOMPADSCOMP,*ADSLANG,*ADSTABLE依赖于岩石性质的吸附数据依赖于岩石性质的吸附数据 *ADSROCK,*ADMAXT*ADMAXT,*ADRT*ADRT,*PORFT,*RRFT*RRFT,*ADSTYPE 二、二、STARSSTARS模块主要关键字模块主要关键字流体定义参数所需资料：流体定义参数所需资料：流体定义参数所需资料：流体定义参数所需资料：粘度温度关系曲线粘度温度关系曲线粘度温度关系曲线粘度温度关系曲线 不同温度下的油水相渗曲线和气液相渗曲线不同温度下的油水相渗曲线和气液相渗曲线不同温度下的油水相渗曲线和气液相渗曲线不同温度下的油水相渗曲线和气液相渗曲线 高压物性高压物性高压物性高压物性PVTPVTPVTPVT数据数据数据数据 毛管压力曲线数据（压汞曲线）毛管压力曲线数据（压汞曲线）毛管压力曲线数据（压汞曲线）毛管压力曲线数据（压汞曲线）组分摩尔质量组分摩尔质量,密度密度,临界压力临界压力,临界温度临界温度,压缩系数压缩系数 中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 INITIAL CONDITIONSINITIAL CONDITIONS初始条件输入标识初始条件输入标识(必需必需)*INITIAL*INITIAL垂向平衡输入标识垂向平衡输入标识(可选可选)*VERTICAL*VERTICAL,*REFPRES,*REFDEPTH,*REFBLOCK原始地层压力及温度原始地层压力及温度 *PRES,*TEMPPRES,*TEMP原始饱和度原始饱和度 *SW,*SO,*SG,*DWOC,*DGOCSW,*SO,*SG,*DWOC,*DGOC原始各相的摩尔分数原始各相的摩尔分数 *MFRAC_WAT,*MFRAC_OIL,*MFRAC_GAS,*PBC原始固相浓度原始固相浓度 *CONC_SLD二、二、STARSSTARS模块主要关键字模块主要关键字初始条件定义所需资料：初始条件定义所需资料：初始条件定义所需资料：初始条件定义所需资料：原始地层压力、温度、压力系数数据原始地层压力、温度、压力系数数据原始地层压力、温度、压力系数数据原始地层压力、温度、压力系数数据 原始油、气、水分布原始油、气、水分布原始油、气、水分布原始油、气、水分布 原始油水界面和油气界面原始油水界面和油气界面原始油水界面和油气界面原始油水界面和油气界面 地质储量报告地质储量报告地质储量报告地质储量报告中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 NUMERICAL METHODS CONTROLNUMERICAL METHODS CONTROL数值控制输入标志数值控制输入标志(可选可选)*NUMERICAL最大时间步长最大时间步长(可选可选)*DTMAX数值模型算法选项数值模型算法选项(可选可选)*TFORM,*ISOTHERMAL*TFORM,*ISOTHERMAL每时步变量的正常变化量每时步变量的正常变化量(可选可选)*NORM松弛选项松弛选项(可选可选)*CONVERGE,*MATBALTOL最大牛顿迭代次数最大牛顿迭代次数(可选可选)*NEWTONCYC低松弛选项低松弛选项(可选可选)*UNRELAX上游权计算选项上游权计算选项(可选可选)*UPSTREAM线性解法的收敛精度线性解法的收敛精度(可选可选)*PRECC最大正交化次数最大正交化次数(可选可选)*NORTH网格块的求解排序网格块的求解排序(可选可选)*SORDER最大迭代次数最大迭代次数(可选可选)*ITERMAX自适应隐式求解方法自适应隐式求解方法(可选可选)*AIM回流关闭开关回流关闭开关可选可选)*BAKFLOSW压压力力和和温温度度的的限限制制值值(可可选选)*MINPRES,*MAXPRES,*MINTEMP,*MAXTEMP每时步最大相转变次数每时步最大相转变次数(可选可选)*PVTOSCMAX二、二、STARSSTARS模块主要关键字模块主要关键字中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 WELL AND RECURRENT DATAWELL AND RECURRENT DATA井和动态数据输入标识井和动态数据输入标识(必需必需)*RUN*RUN模拟参考时间模拟参考时间 *TIME,*DATE,*DTWELLTIME,*DATE,*DTWELL井定义标识井定义标识(必需必需)*WELL*WELL井组定义井组定义(可选可选)*GROUP,*GROUPWT井类型定义井类型定义(必需必需)*PRODUCER,*INJECTOR,*SHUTIN,*OPEN*PRODUCER,*INJECTOR,*SHUTIN,*OPEN注入蒸汽特性参数注入蒸汽特性参数 *TINJW,*QUALTINJW,*QUAL,*PINJW,*TINJOV*TINJOV,*HEATLOSS注入相的含量注入相的含量 *INCOMPINCOMP井控制选项井控制选项(必需必需)*OPERATE*OPERATE控制条件转换开关控制条件转换开关(可选可选)*CONSTRNCHK井监测选项井监测选项(可选可选)*MONITOR井筒几何特征井筒几何特征(条件条件)*GEOMETRY*GEOMETRY井的完井层位井的完井层位(条件条件)*PERF*PERF垂直井的完井层位垂直井的完井层位(条件条件)*PERFV加密网格完井层位加密网格完井层位(条件条件)*PERFRG改变井动态数据选项改变井动态数据选项 *ALTERALTER模拟终止标志模拟终止标志(必需必需)*STOP*STOP 二、二、STARSSTARS模块主要关键字模块主要关键字井定义即生产数据所需资料：井定义即生产数据所需资料：井定义即生产数据所需资料：井定义即生产数据所需资料：射孔完井及井史报告；射孔完井及井史报告；射孔完井及井史报告；射孔完井及井史报告；试油、试井（压力恢复曲线）和试采资料试油、试井（压力恢复曲线）和试采资料试油、试井（压力恢复曲线）和试采资料试油、试井（压力恢复曲线）和试采资料 系统测压资料（压力恢复曲线）系统测压资料（压力恢复曲线）系统测压资料（压力恢复曲线）系统测压资料（压力恢复曲线）吸水剖面和产液剖面资料吸水剖面和产液剖面资料吸水剖面和产液剖面资料吸水剖面和产液剖面资料 井位示意图；井位示意图；井位示意图；井位示意图；区块油井生产数据报表：区块油井生产数据报表：区块油井生产数据报表：区块油井生产数据报表：日产油、日产液、日产气日产油、日产液、日产气日产油、日产液、日产气日产油、日产液、日产气 综合含水、压力综合含水、压力综合含水、压力综合含水、压力 累积产油（气、水、液）累积产油（气、水、液）累积产油（气、水、液）累积产油（气、水、液）区块综合生产数据统计报表：区块综合生产数据统计报表：区块综合生产数据统计报表：区块综合生产数据统计报表：日产油（水、气、液）日产油（水、气、液）日产油（水、气、液）日产油（水、气、液）采出程度、综合含水采出程度、综合含水采出程度、综合含水采出程度、综合含水 累积产油（气、水、液）；累积产油（气、水、液）；累积产油（气、水、液）；累积产油（气、水、液）；热采模型热采模型热采模型热采模型中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置FLUID AND COMPONENT DEFINITIONSFLUID AND COMPONENT DEFINITIONSROCK-FLUID PROPERTIESROCK-FLUID PROPERTIESWELL AND RECURRENT DATAWELL AND RECURRENT DATA组分特性的选择及计算组分特性的选择及计算组分特性的选择及计算组分特性的选择及计算两相相渗内插函数的参数选择两相相渗内插函数的参数选择两相相渗内插函数的参数选择两相相渗内插函数的参数选择泡沫特性参数的选择泡沫特性参数的选择泡沫特性参数的选择泡沫特性参数的选择吸附特性参数的选择吸附特性参数的选择吸附特性参数的选择吸附特性参数的选择注入井注入参数的选择及计算注入井注入参数的选择及计算注入井注入参数的选择及计算注入井注入参数的选择及计算中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置组分特性的选择及计算组分特性的选择及计算组分特性的选择及计算组分特性的选择及计算*MODEL 6 6 4 2 *组分总个数组分总个数 油气水中组分总个数油气水中组分总个数 油水中组分总个数油水中组分总个数 (水相中组分总个数水相中组分总个数)*COMPNAME WATER SURFACT OIL SOL_GAS N2 LAMELLA*-*CMM 0.018 0.308 0.4 0.018 0.028 0.018 *各组分的摩尔质量，各组分的摩尔质量，Kg/mol*MASSDEN 1000.0 1137.0 965.0 14.46 *各组分的质量密度各组分的质量密度mol/m3*PCRIT 22048 992.8 1475.5 4544 3394 3394 *临界压力，注意单位临界压力，注意单位KPa*TCRIT 374.0 527.0 434.0 -82.72 -147.0 -147.0 *临界温度，单位，摄氏度临界温度，单位，摄氏度 在在STARSSTARS中中，把把气气泡泡和和液液膜膜看看作作一一种种稳稳定定的的分分散散组组分分，存存在在于于某某相相或或某某几几相相中中的的。把把泡泡沫沫的的各各种种特特征征作作为为该该分分散散组组分分的的性性质质来来处处理理，包包括括吸吸附附特特性性、封封堵堵孔孔道道的的特特性性、非非线线性性粘粘度度特特性性、扩扩散散传传播播等等。通通过过合合理理选选择择输输入入数数据据，适适当当的的选选取取该该分分散散组组分分的的性性质质，可可以以描描述述泡泡沫沫流流现象。现象。中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置组分特性的选择及计算组分特性的选择及计算组分特性的选择及计算组分特性的选择及计算 已已知知：发发泡泡剂剂原原液液密密度度为为1.044g/cm3，有有效效浓浓度度为为35，发发泡泡剂剂的的摩摩尔尔质质量量为为0.308kg/mol，水水的的摩摩尔尔质质量量为为0.018kg/mol，现现在在求求纯发泡剂的质量密度。纯发泡剂的质量密度。假假设设单单位位体体积积发发泡泡剂剂溶溶液液中中，发发泡泡剂剂的的体体积积为为x，则则水水的的体体积积为为（1-x），若若纯纯发发泡泡剂剂的的质质量量密密度度为为DENs，水水的的密密度度为为DENw，则则可可以以求求出出单单位位体体积积发发泡泡剂剂溶溶液液中中纯纯表表活活剂剂所所占占的的体体积积和和质质量量密密度度（MASSDENMASSDEN），同同样样可可以以求求出出对对应应的的摩摩尔尔密密度度（MOLDENMOLDEN）和摩尔体积和摩尔体积（MOLVOLMOLVOL）。各组分质量密度、摩尔密度、摩尔体积的计算方法各组分质量密度、摩尔密度、摩尔体积的计算方法各组分质量密度、摩尔密度、摩尔体积的计算方法各组分质量密度、摩尔密度、摩尔体积的计算方法中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置两相相渗内插函数的参数选择两相相渗内插函数的参数选择两相相渗内插函数的参数选择两相相渗内插函数的参数选择 高界面张力和低界面张力间的插值计算：高界面张力和低界面张力间的插值计算：考考虑虑一一种种岩岩石石类类型型，常常规规的的油油水水相相渗渗与与高高界界面面张张力力（无无表表活活剂剂存存在在）所所对对应应；当当表表活活剂剂加加入入该该系系统统时时，残残余余饱饱和和度度Sr降降低低，润润湿湿性性发发生生改改变变。最最后后，由由于于高高的的表表活活剂剂浓浓度度，极极低低的的界界面面张张力力，导致导致Sr降低，降低，Kr成直线。成直线。利利用用内内插插函函数数DTRAPWDTRAPW和和DTRAPNDTRAPN在在高高界界面面张张力力时时Kr曲曲线线和和极极低低界界面面张张力力时时Kr曲曲线线间间插插值值，从从而而得得到到不不同同界界面面张张力力时时的的相相渗渗曲曲线线。所所需需的的参参数数包包括括：插插值值组组分分（INCOMPINCOMP）、界界面面张张力力数数据据（IFTABLEIFTABLE）、泡泡沫沫插插值值函函数数（FMSURFFMSURF）、相相渗渗内内插插函函数数（DTRAPWDTRAPW和和和和DTRAPN DTRAPN）、高高低低界界面面张张力力时时的的相相渗渗曲曲线线（SWTSWT和和和和SLT SLT）。中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置两相相渗内插函数的参数选择两相相渗内插函数的参数选择两相相渗内插函数的参数选择两相相渗内插函数的参数选择 高界面张力和低界面张力间的插值计算高界面张力和低界面张力间的插值计算有无表活剂时油水相渗曲线图有无表活剂时油水相渗曲线图有无表活剂时气液相渗曲线图有无表活剂时气液相渗曲线图中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置泡沫插值函数泡沫插值函数泡沫插值函数泡沫插值函数*FMSURF fmsurf 临界表活剂浓度临界表活剂浓度*FMCAP fmcap 临界毛管数临界毛管数*FMOIL fmoil 临界油饱和度临界油饱和度*FMGCP fmgcp 临界总毛管数临界总毛管数*FMOMF fmomf 临界油摩尔分数临界油摩尔分数*FMMOB fmmob 流度降低因子流度降低因子*EPSURF epsurf 临界摩尔分数指数临界摩尔分数指数*EPCAP epcap 临界毛管数指数临界毛管数指数*EPOIL epoil 临界含油饱和度指数临界含油饱和度指数*EPGCP epgcp 临界总毛管数指数临界总毛管数指数*EPOMF epomf 临界油摩尔分数指数临界油摩尔分数指数 只只有有当当*INCOMP和和*IFTABLE都都出出现现时时，泡泡沫沫插插值值函函数数有有效效。FMMOB的的取取值值：在在临临界界表表活活剂剂浓浓度度FMSURF时时，So=0.0，油摩尔分数为，油摩尔分数为0时测出。时测出。典型值：典型值：FMSURF=0.00001,FMCAP=0.0001,FMOIL=0.2,FMGCP=1.0*10-6,FMOMF=0.2,EPSURF=1.0,EPCAP=0.5,EPOIL=1.0,EPGCP=1.0,EPOMF=1.0。当当考考虑虑毛毛管管数数的的变变化化时时，收收敛敛性性较较差差，可取可取EPCAP=0中国石油大学（北京）石油天然气工程学院油藏数值模拟组中国石油大学（北京）石油天然气工程学院油藏数值模拟组 三、三、STARSSTARS模块泡沫的模拟模块泡沫的模拟与泡沫有关的设置与泡沫有关的设置吸附特性参数吸附特性参数吸附特性参数吸附特性参数*ADSCOMP SURFACT *WATER 吸附组分名吸附组分名 组分所在相组分所在相*ADMAXT 4.0E-