ChemSim软件操作说明(中文).doc

油藏化学驱和调驱调剖模拟软件系统
（ChemSim）操作说明

2010年8月

目 录
一．软件概述
二、软件安装
三、模型建立与三维数值模拟的简捷操作说明
附录：ChemSim软件参数
1、标题和油藏描述
2、油藏性质
3、物理性质
4、历史拟合数据
5、输出选项

一、软件概述
ChemSim是美国IUT集团联合美国一个化学驱精英团队开发的最新一代三维油田化学驱和调驱调剖软件，在国内外油田大量需求中最近隆重推出。它是目前国际上功能最强，物理化学机理性计算最完备，实验室数据基础最充分的化学驱及调驱调剖专业性很强的软件系统。
ChemSim是集油藏建模、数值模拟于一体的大型综合油藏数值模拟软件。油藏注水调驱调剖模拟技术软件数值模拟引擎采用了目前油藏数值模拟研究的最新成果和多重网格高速算法，数值模拟运算速度较快。ChemSim除了用于数值模拟计算，还可以对模拟结果进行处理，因为它本身具有模型后处理显示功能及曲线可视化。利用OpenGL技术开发的功能强大的后处理三维图像显示模块可以三维显示油藏数值模拟器计算的结果以及油藏地质结构体、油藏厚度分布、有效厚度分布、孔隙度分布、渗透率分布、压力场分布、饱和度场分布等。

主要有以下特点：
1、三维非等温的油藏化学驱数值模拟器
2、含有直角坐标和柱坐标两种网格坐标
3、适用于直井、水平井、斜井
4、模型中有关交联调剖运算部分的数学描述主要基于质量作用定律
5、三阶有限差分网格
6、岩石表面存在离子交换
7、示踪剂的双重孔隙度模型以及示踪剂的分配，反应，吸附和放射性衰退
8、聚合物的粘弹性质模拟
9、颗粒胶体运移和膨胀
10、模拟中井的类型可以改变
11、自动时间步选择及精度、收敛、稳定性控制
12、考虑了全张量离散系数和分子弥散、粘土/活性剂阳离子交换

二、软件安装
1、启动ESSS-ChemSim Setup安装程序
2、按要求完成安装步骤后，打开ChemSim主界面，选择Help目录下的License Info，如图1所示。

图1 License Info
3、生成机器码（物理地址），如图2所示，并将其反馈给本公司相关人员

图2生成机器码
4、经公司总部审查，得到许可证，并加载许可证信息（文件类型为.LIC），如图3所示。

图3 加载许可证信息

5、许可证加载完成后，软件即可应用

三、模型建立与三维数值模拟的简捷操作说明
1、打开ChemSim软件到主界面
2、新建并输入软件中已存在的ChemSim模型文件夹，如图4所示。

图4新建并输入软件中已存在的ChemSim模型文件夹
3、一个新项目文件已经建立，可以对其中每一项进行描述，如参数的设置、选择和加载文件（可加载文件类型包括：.dat和.txt），如图5所示。
加载文件

图5 参数的设置、选择和加载文件
4、模型参数设置完成后，将其保存在已输入模型文件夹目录下，如图6所示。

图6 保存
5、进行数值模拟计算，如图7所示。

图7 数值模拟计算

6、计算完成之后，可以观看三维图像显示及曲线可视化，如图8、图9所示。

图8 三维图像显示

图9 曲线可视化

附录： ChemSim软件参数
Model Options:
RUNNO: Run number运行编号
TITLE: Title and run description标题和运行描述
IMES: Time-Step Mode时步长类型 - Flag indicating if a constant or automatic time-step is to be used.
1-Constant 常数
2-Automatic 自动
IMODE: Setup Mode模拟类型 - Flag indicating if the problem to be run is a first run or a restart simulation.
1-First Run 第一次运算
2-Restart 重新运算
IDISPC: Numerical Dispersion Control数值弥散控制条件-Flag indicating which type of numerical dispersion control is used.
0 - Single point upstream method is used
3 - Improved total variation diminishing third order method is used
IREACT: Gel Reactions有无凝胶反应- Flag indicating if gel reactions or alkaline options are used or not.
0 - Gel reactions are not used
1 - Gel reactions are used
ICOORD: Coordinate System坐标，直角坐标或柱坐标- Flag indicating which coordinate system is used.
1 - Cartesian
2 - Radial
ITREAC: Tracer Reactions有无示踪剂反应- Flag indicating if a tracer reaction is used or not.
0 - Tracer reactions are not used
1 - Tracer reactions are used
ITC: Second Order Time Approximation是否有二阶近似- Flag indicating if second-order time approximation is used or not.
0 - Second-order time approximation is not used
1 - Second-order time approximation is used
IENG: Thermal Variation是否考虑温度变化- Flag indicating if temperature variation is considered or not.
0 - Isothermal simulation
1 - Temperature equation is solved
IUNIT: Unit System单位，公制或英制- Flag indicating English or Metric units.
0 - English
1 – Metric
Components:
COMP_NAMES: Component Names组分名称-Names of the components in the model.
1.Water 2.Oil 3.Surfactant 4.Polymer
5.Chloride 6.Calcium 7.Alcohol1 8.Alcohol2
TRACER_NAMES: Names of the Water/Oil Tracers水/油示踪剂名称- Names of the Water/Oil Tracers in the model.
GEL_NAMES: Names of the Gel Components凝胶组分名称- Names of the Gel Components in the model.
ITRU_VALUES: Water Tracer Units水示踪剂单元-Flag indicating the units of the Ith water tracer.
Calculations:
ICF_VALUES: Component Included参与计算的组分选择- Flag indicating if given component is included in the calculations or not.
Reservoir Properties:
COMPR: Rock compressibility岩石可压缩性
PSTAND: Reference pressure基准压力- Reference pressure at which pore volume and fluid compressibilities are specified.
Grid:
NX,NY,NZ: Number Of Cells X/Y/Z - X/Y/Z方向的网格数
IDXYZ: Grid Geometry Mode网格几何参数类型- Flag indicating constant or variable grid size.
DX1,DY1,DZ1: Spacing X/Y/Z - X/Y/Z方向的网格大小

Rock Properties:
IPOR1: Reservoir Porosity 油藏孔隙度类型-Flag indicating constant or variable porosity for reservoir.
0 - Constant porosity for whole reservoir
2 - Variable porosity over reservoir
4 - Porosity values for each gridblock are read from POR input file
PORC1: Reservoir Porosity油藏孔隙度
PORI_VALUES: 加载孔隙度场
IPERMX: Flag indicating constant or variable X direction permeability-X方向的渗透率类型
0 - Constant permeability for whole reservoir
2 - Variable permeability over reservoir
4 - X direction permeability values for each gridblock are read from PERMX input file
IPERMY: Flag indicating constant or variable Y direction permeability for reservoir- Y方向的渗透率类型
0 - Constant permeability for whole reservoir
2 - Variable permeability over reservoir
3 - Y direction permeability is dependent on X direction permeability
4 - Y direction permeability values for each gridblock are read from PERMY input file
IPERMZ - Flag indicating constant or variable Z direction permeability for reservoir. - Z方向的渗透率类型
0 - Constant permeability for whole reservoir
2 - Variable permeability over reservoir
3 - Z direction permeability is dependent on X direction permeability
4 - Z direction permeability values for each gridblock are read from PERMZ input file
PERMX_VALUES、PERMY_VALUES、PERMZ_VALUES: Permeability of the each gridblock in the X/Y/Z direction.分别加载X/Y/Z方向渗透率场
FACTY: Constant Permeability multiplier for Y direction -Y方向的渗透率乘数
FACTZ: Constant Permeability multiplier for Z direction Z方向的渗透率乘数
ITRANZ: Transmissibilities Modified传导率是否被修正- Flag indicating whether the transmissibilities are modified or not.
0 - No transmissibility modifier
1 - Use transmissibility modifier
INTG: Net-to-gross Values Modify Gridblock Pore Volumes是否用净毛比值修正网格孔隙体积- Flag indicating whether net-to-gross values are input to modify gridblock pore volumes.
0 - Net-to-gross values are not input
1 - Net-to-gross values are read from NTG input file
Initial Properties:
IDEPTH: Top Layer Depth Measurement Type顶层深度类型- Flag indicating type of depth measurement of the top layer.
0 - Single value for depth of the top layer is specified
1 - Depth of top gridblock (1,1,1) and the reservoir dip angles are specified
2 - Depth of each gridblock in the top layer is specified
4 - Depth values for each gridblock in the top layer are read from TOP input file
D111: Measured Reservoir Top Layer Depth顶层深度-Depth of the top layer of the reservoir measured from the surface (reference plane), positive downward.
D111_DIP: First Gridblock Depth网格坐标为（1,1,1）的深度- Depth of the first gridblock (1,1,1).
THETAX: Reservoir Dip Angle in X Direction X方向的油藏倾角- Reservoir dip angle in X direction, positive downward.
THETAY: Reservoir Dip Angle in Y Direction Y方向的油藏倾角- Reservoir dip angle in Y direction, positive downward.
ELI_VALUES: Depth of each gridblock in the top layer加载顶层深度场
ELI_VALUES should have 10000 values defined.
IPRESS: Reservoir Initial Pressure Measurement初始油藏压力类型- Flag indicating type of reservoir initial pressure measurement.
0 - Single value for reservoir initial pressure is used for all gridblocks
1 - Initial pressure for a point at a specified depth is specified by user
2 - Initial pressure for each gridblock is specified by user
4 - Initial pressure for each gridblock is read from PRES input file
PRESS1: Initial Reservoir Pressure初始油藏压力值
PINIT: Initial Reservoir Pressure at HINIT Depth指定深度点所对应的初始油藏压力
HINIT: PINIT Initial Pressure Depth指定的深度点- Depth of the point where PINIT initial pressure is specified.
PI_VALUES: Initial pressure of each gridblock加载初始油藏压力场-Initial pressure of each gridblock in the reservoir.
PI_VALUES should have 50000 values defined.
ISWI: Initial Water Saturation Measurement初始含水饱和度类型- Flag indicating type of initial water saturation measurement.
0 - Single value for initial water saturation is used for all gridblocks
2 - Initial water saturation for each gridblock is specified by user
3- Initial water saturation according to oil water contact (OWC)
4 - Initial water saturation for each gridblock is read from WSAT input file
SWI: Initial Water Saturation for all Gridblocks 初始含水饱和度-Initial water saturation for all gridblocks of the reservoir.
SI_VALUES: Initial water saturation of each gridblock加载初始含水饱和度场
SI_VALUES should have 50000 values defined.
OWC: Oil-Water Depth Contact油水界面深度-Depth to oil-water contact
SWC: Residual water Saturation for oil-water zone油水区束缚水饱和度
Brine Properties:
C50：Constant initial brine salinity初始盐度
C60：Initial divalent cation concentration of brine初始二价阳离子浓度
Property Modification:
IMOD: Reservoir Properties Modified油藏性质修正-Flag indicating whether the reservoir properties are modified or not.
0 - No property is modified
1 - Allow for property modification
IMPOR: Porosity Modified是否修正孔隙度值- Flag indicating whether the porosity is modified or not.
0 - No modification is considered in porosity values
1 - Allow modification in porosity
REGIONS_PORO: Porosity indexes孔隙度修正区域
IMIN_PORO: The first index in X direction. 孔隙度修正区域的X方向首坐标
IMAX_PORO: The last index in X direction.孔隙度修正区域的X方向末坐标
JMIN_PORO: The first index in Y direction. 孔隙度修正区域的Y方向首坐标
JMAX_PORO: The last index in Y direction.孔隙度修正区域的Y方向末坐标
KMIN_PORO: The first index in Z direction. 孔隙度修正区域的Z方向首坐标
KMAX_PORO: The last index in Z direction.孔隙度修正区域的Z方向末坐标
IFACT_PORO: How porosity is modified.孔隙度修正方式
1 - Replace porosity with FACTX
2 - Multiply porosity by FACTX
3 - Add FACTX to porosity
FACTX_PORO: The constant used to modify the porosity value.用于修正孔隙度值的常数
IMKX: Permeability in the X Direction Modified. 是否修正X方向渗透率值-Flag indicating whether the permeability in the X direction is modified or not.
IMKY: Permeability in the Y Direction Modified. 是否修正Y方向渗透率值-Flag indicating whether the permeability in the Y direction is modified or not.
0 - No modification is considered in Y permeability
1 - Allow modification in Y permeability
IMKZ: Permeability in the Z Direction Modified. 是否修正Z方向渗透率值-Flag indicating whether the permeability in the Z direction is modified or not.
0 - No modification is considered in Z permeability
1 - Allow modification in Z permeability
IMSW: Initial water saturation modified. 是否修正初始含水饱和度值- Flag indicating whether the initial water saturation is modified or not.
0 - No modification is considered in initial water saturation values
1 - Allow modification in initial water saturation
General Physical Properties:
Phase Behavior:
Surfactant:
EPSME: Critical micelle concentration (CMC)临界胶束浓度-minimum surfactant concentration for the formation of micelles.
HBNC70: Intercept of maximum height of binodal curve at zero salinity. 在零盐度时双结点曲线最高点的截距
HBNC71: Intercept of maximum height of binodal curve at optimal salinity.在最佳盐度时双结点曲线最高点的截距
HBNC72: Intercept of maximum height of binodal curve at at twice optimal salinity.在2倍最佳盐度时双结点曲线最高点的截距
CSEL7: Lower effective salinity limit for Type III phase region determined when Alcohol 1 and calcium approach zero.当醇1与Ca浓度接近0时，模型Ⅲ的有效盐度下限值
CSEU7: Upper effective salinity limit for Type III phase region determined when Alcohol 1 and calcium approach zero.当醇1与Ca浓度接近0时，模型Ⅲ的有效盐度上限值
BETA6: Ca有效盐度的斜率- The effective salinity slope parameter for calcium.
Cosolvent1:
HBNS70: Slope for maximum height of binodal curve vs.fraction of Alcohol 1 associated with surfactant at zero salinity. 双结点曲线最高点的斜率与在零盐度时添加表活剂的醇1分数之比
HBNS71: Slope for maximum height of binodal curve vs.fraction of Alcohol 1 associated with surfactant at optimal salinity. 双结点曲线最高点的斜率与在最佳盐度时添加表活剂的醇1分数之比
HBNS72: Slope for maximum height of binodal curve vs.fraction of Alcohol 1 associated with surfactant at twice optimal salinity. 双结点曲线最高点的斜率与在2倍最佳盐度时添加表活剂的醇1分数之比
BETA7: The effective salinity slope parameter for Alcohol 1. 醇1有效盐度的斜率
Cosolvent2:
HBNS80: Slope for maximum height of binodal curve vs. fraction of Alcohol 2 associated with surfactant at zero salinity.
HBNC80: Intercept of maximum height of binodal curve at zero fraction of Alcohol 2 associated with surfactant at zero salinity.
HBNS81: Slope of maximum height of binodal curve vs. fraction of Alcohol 2 associated with surfactant at optimal salinity.
HBNC81: Intercept of maximum height of binodal curve at zero fraction of Alcohol 2 associated with surfactant at optimal salinity.
HBNS82: Slope for maximum height of binodal curve vs. fraction of Alcohol 2 associated with surfactant at twice optimal salinity.
HBNC82: Intercept of maximum height of binodal curve at zero fraction of Alcohol 2 associated with surfactant at twice optimal salinity.
CSEL8: 当醇2与Ca浓度接近0时，模型Ⅲ的有效盐度下限值- Lower effective salinity limit for Type III phase region determined when Alcohol 2 and calcium approach zero.
CSEU8: 当醇2与Ca浓度接近0时，模型Ⅲ的有效盐度上限值- Upper effective salinity limit for Type III phase region determined when Alcohol 2 and calcium approach zero.
BETA8: 醇2有效盐度的斜率- The effective salinity slope parameter for Alcohol 2.
Alcohol Partitioning:
IALC: 乙醇分离模型选择- Flag indicating choice of alcohol partition model to use.
0 - Hirasaki's model will be used
1 - Prouvost's model will be used
OPSK70: 对于醇1的乙醇分离系数（油/水）- Alcohol partition coefficient (oil/water) for Alcohol 1.
OPSK7S: 对于醇1的乙醇分离系数（表活剂/水）- Alcohol partition coefficient (surfactant/water) for Alcohol 1.
OPSK80: 对于醇2的乙醇分离系数（油/水）- Alcohol partition coefficient (oil/water) for Alcohol 2.
OPSK8S: 对于醇2的乙醇分离系数（表活剂/水）- Alcohol partition coefficient (surfactant/water) for Alcohol 2.
NALMAX: 乙醇分离的最大迭代次数- Maximum number of iterations for alcohol partitioning for two alcohol system.
EPSALC: 迭代收敛公差- Tolerance for convergence of iterations for two alcohol system.
AKWC7、AKWS7: 决定水、油两相之间单体醇1的分配系数的参数- Parameters used to determine partition coefficient of monomeric Alcohol 1 between aqueous and oleic pseudophases.
AKM7: 表面活性剂和油相之间单体醇1的分配系数- Partition coefficient of monomeric Alcohol 1 between surfactant and oleic pseudophases.
AK7: 油相中醇1的自缔合常数- Self-association constant of Alcohol 1 in oleic pseudophase.
PT7: 醇1的摩尔体积与等量的表面活性剂摩尔体积之比- Ratio of molar volume of Alcohol 1 to equivalent molar volume of surfactant.
AKWC8、AKWS8: 决定水、油两相之间单体醇2的分配系数的参数- Parameters used to determine partition coefficient of monomeric Alcohol 2 between aqueous and oleic pseudophases.
AK8: 油相中醇2的自缔合常数- Self-association constant of Alcohol 2 in oleic pseudophase.
PT8: 醇2的摩尔体积与等量的表面活性剂摩尔体积之比- Ratio of molar volume of Alcohol 2 to equivalent molar volume of surfactant.
Temperature Effect:
HBNT0: Slope of height of binodal curve versus temperature at zero salinity
HBNT1: Slope of height of binodal curve versus temperature at optimal salinity
HBNT2: Slope of height of binodal curve versus temperature at twice optimal salinity
CSET: The Slope parameter, , for temperature dependency of the three-phase window
Interfacial Tension:
IFT: 界面张力关系类型- This flag indicates which interfacial tension correlation is used.
0 - Healy and Reed's correlat