Vip Core
Vip Core
Vip Core
Reference Manual
© 2001, 2002 Landmark Graphics Corporation
Trademark Notice
Landmark, the Landmark logo, 3D Drill View, 3D Drill View KM, 3DVIEW, Active Field Surveillance, Active
Reservoir Surveillance, ARIES, Automate, BLITZ, BLITZPAK, CasingSeat, COMPASS, Contouring Assistant,
DataStar, DBPlot, Decision Suite, Decisionarium, DecisionDesktop, DecisionSpace, DepthTeam, DepthTeam
Explorer, DepthTeam Express, DepthTeam Extreme, DepthTeam Interpreter, DESKTOP-PVT, DESKTOP-VIP,
DEX, DFW, DIMS, Discovery, Drillability Suite, DrillModel, DrillVision, DSS, Dynamic Surveillance System,
EarthCube, EdgeCa$h, eLandmark, EPM, e-workspace, FastTrack, FZAP!, GeoDataLoad, GeoGraphix, GeoGraphix
Exploration System, GeoLink, GES, GESXplorer, GMAplus, GrandBasin, GRIDGENR, I2 Enterprise, iDims,
IsoMap, LandScape, LeaseMap, LMK Resources, LogEdit, LogM, LogPrep, Make Great Decisions, MathPack,
Model Builder, MyLandmark, MyWorkspace, OpenBooks, OpenExplorer, OpenJournal, OpenSGM, OpenTutor,
OpenVision, OpenWorks, OpenWorks Well File, PAL, Parallel-VIP, PetroBank, PetroWorks, PlotView, Point
Gridding Plus, Pointing Dispatcher, PostStack, PostStack ESP, PRIZM, PROFILE, ProMAX, ProMAX 2D, ProMAX
3D, ProMAX 3DPSDM, ProMAX MVA, ProMAX VSP, pStaX, QUICKDIF, RAVE, Real Freedom, Reservoir
Framework Builder, RESev, ResMap, RMS, SafeStart, SCAN, SeisCube, SeisMap, SeisModel, SeisSpace,
SeisVision, SeisWell, SeisWorks, SeisXchange, SigmaView, SpecDecomp, StrataMap, Stratamodel, StratAmp,
StrataSim, StratWorks, StressCheck, STRUCT, SynTool, SystemStart, T2B, TDQ, TERAS, Total Drilling
Performance, TOW/cs, TOW/cs The Oilfield Workstation, Trend Form Gridding, Turbo Synthetics, VIP, VIP-COMP,
VIP-CORE, VIP-DUAL, VIP-ENCORE, VIP-EXECUTIVE, VIP-Local Grid Refinement, VIP-POLYMER, VIP-
THERM, WavX, Web OpenWorks, Well Editor, Wellbase, Wellbore Planner, WELLCAT, WELLPLAN,
WellXchange, wOW, Xsection, ZAP!, Z-MAP Plus are trademarks, registered trademarks or service marks of
Landmark Graphics Corporation.
All other trademarks are the property of their respective owners.
Note
The information contained in this document is subject to change without notice and should not be construed as a
commitment by Landmark Graphics Corporation. Landmark Graphics Corporation assumes no responsibility for any
error that may appear in this manual. Some states or jurisdictions do not allow disclaimer of expressed or implied
warranties in certain transactions; therefore, this statement may not apply to you.
Chapter
Table of Contents
Table of Contents
Preface
About This Manual
Purpose . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xix
The Modules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xix
The Chapters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xx
Data Formatting Conventions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xx
Compatibility . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxi
Related Documentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxi
Chapter 1
Data Overview
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1
1.1.1 VIP-COMP Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1
1.1.2 VIP-ENCORE Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-2
1.1.3 VIP-DUAL Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-2
1.1.4 VIP-POLYMER Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-3
1.1.5 VIP-THERM Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-3
1.1.6 Shared Features . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-4
1.2 Typical Data Requirements to Initialize a Simulation Study . . . . . . . . . . . . 1-6
1.2.1 Geological Descriptions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-6
1.2.2 Reservoir Rock Characteristics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-6
1.2.3 Hydrocarbon Fluid Properties . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-6
1.3 Data Deck Layout . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-7
1.4 Input Data Template . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-9
1.5 Data Format . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-11
1.5.1 General Utility Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-13
1.5.1.1 Comment Lines (C) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-13
1.5.1.2 Read Data from an External File (INCLUDE) . . . . . . . . . 1-13
1.5.1.3 Stop Reading Data from the Current INCLUDE File
(ENDINC) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-14
Chapter 2
Initialization Data
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-47
2.2 Initialization Utility Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-49
2.2.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-49
2.2.1.1 Initialization (INIT) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-49
2.2.1.2 Change Default Dimensions (DIM) . . . . . . . . . . . . . . . . . . 2-49
iv Landmark - R2003.4
VIP-CORE REFERENCE MANUAL
R2003.4 - Landmark v
VIP-CORE REFERENCE MANUAL
2.2.7.5
Gas-Oil Capillary Pressure Hysteresis (PCHYSG) . . . . . 2-90
2.2.7.6
Oil Relative Permeability Hysteresis (RPHYSO) . . . . . . 2-91
2.2.7.7
Gas Relative Permeability Hysteresis (RPHYSG) . . . . . 2-91
2.2.7.8
Relative Permeability Hysteresis Tolerances
(RPHYST) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-92
2.2.7.9 Leverett J-Function (JFUNC) (Not available with
SDFUNC option) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-92
2.2.7.10 Rock Compaction (COMPACT) . . . . . . . . . . . . . . . . . . . 2-93
2.2.7.11 Freeze Pcwo at Initial Value (FRZPCW) . . . . . . . . . . . . . 2-93
2.2.7.12 Freeze Pcgo at Initial Value (FRZPCG) . . . . . . . . . . . . . . 2-94
2.2.8 Initialization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-94
2.2.8.1 Nonequilibrium Initialization (NONEQ) . . . . . . . . . . . . . 2-94
2.2.8.2 Gridblock Center Initialization (GBC) . . . . . . . . . . . . . . . 2-94
2.2.8.3 Integrated Saturation Initialization (INTSAT) . . . . . . . . . 2-94
2.2.8.4 Integrated Saturation Initialization (VAITS) . . . . . . . . . . 2-95
2.2.8.5 Thickness Center (THCNTR) . . . . . . . . . . . . . . . . . . . . . . 2-96
2.2.8.6 Do Not Initialize (NOINIT) . . . . . . . . . . . . . . . . . . . . . . . . 2-96
2.2.8.7 Totally Refined ROOT Grid (NOROOT) . . . . . . . . . . . . . 2-96
2.2.8.8 Grid Deactivation (DEACTIVATE) . . . . . . . . . . . . . . . . . . 2-97
2.2.9 Off-Band Connections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-97
2.2.9.1 Pinchout Gridblock Connections (PINCHOUT) . . . . . . . 2-97
2.2.9.2 Fault Modeling (FAULTS) . . . . . . . . . . . . . . . . . . . . . . . . . 2-98
2.2.9.3 Completing the Circle in Radial Grids (FLOW360) . . . . 2-98
2.2.10 Vertical Equilibrium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-99
2.2.10.1 Water-Oil Vertical Equilibrium (VEWO) . . . . . . . . . . . . 2-99
2.2.10.2 Gas-Oil Vertical Equilibrium (VEGO) . . . . . . . . . . . . . 2-100
2.2.10.3 Vertical Equilibrium Directional Relative
Permeability (DRELPM) . . . . . . . . . . . . . . . . . . . . . . . . 2-100
2.2.10.4 Capillary Gravity Equilibrium (VEITS) . . . . . . . . . . . . 2-101
2.2.11 Fluid Property Options . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-102
2.2.11.1 Energy Minimization Phase Equilibrium (GIBBS)
(VIP-COMP) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-102
2.2.11.2 Near Critical Fluid Property Adjustment (IFT)
(VIP-COMP) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-102
2.2.11.3 Suppresses Table Data Checking (NOCHK) . . . . . . . . 2-103
2.2.11.4 PVT Interpolation for VIP-ENCORE (BOTINT) . . . . . 2-103
2.2.11.5 Flash Calculation Method (FLASH) (VIP-COMP) . . . 2-104
2.2.11.6 Super-Critical Equilibration (CRINIT) (VIP-COMP) . 2-105
2.2.11.7 Dry Gas Simulation (DRYGAS) . . . . . . . . . . . . . . . . . . . 2-106
2.2.11.8 Limit on Rate of Increase of Solution GOR (DRSDT) . 2-106
2.2.12 Corner-Point Grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-107
2.2.12.1 Corner-Point Simulation Grid (CORNER) . . . . . . . . . 2-107
2.2.12.2 Fault Connections from Corner-Point Data
(CORTOL) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-108
2.2.12.3 Data Checking Corner-Point Grid Data (CORCHK) . 2-109
2.2.13 Dual Porosity with Optional Dual Permeability (VIP-DUAL) . 2-110
vi Landmark - R2003.4
VIP-CORE REFERENCE MANUAL
Chapter 3
Print Control
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-125
3.2 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-125
3.2.1 Print Everything (ALL) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-125
3.2.2 Print Nothing (NONE) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-125
3.3 Individual Group Print Controls . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-126
3.3.1 Input Array Printing (PRINT ARRAYS) . . . . . . . . . . . . . . . . . . . . 3-126
3.3.2 Coefficient Array Printing (PRINT COEFS) . . . . . . . . . . . . . . . . . 3-127
3.3.3 Compositional Data Printing (PRINT COMP) . . . . . . . . . . . . . . . 3-131
3.3.4 Corner-Point Data Printing (PRINT CORNER) . . . . . . . . . . . . . . 3-131
3.3.5 Equilibrium Data Printing (PRINT EQUIL) . . . . . . . . . . . . . . . . . 3-132
3.3.6 Fault Data Printing (PRINT FAULTS) . . . . . . . . . . . . . . . . . . . . . . 3-132
3.3.7 Influx Data Printing (PRINT INFLUX) . . . . . . . . . . . . . . . . . . . . . 3-133
3.3.8 Initialization Arrays . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-134
3.3.9 Region Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-139
3.3.10 Separation Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-139
3.3.11 Tabular Data - Saturation and PVT tables . . . . . . . . . . . . . . . . . . 3-139
3.4 Rescaled Saturation Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-141
Chapter 4
Tables
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-145
4.2 Equilibrium Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-145
4.2.1 Saturation Pressure is Constant by Regions . . . . . . . . . . . . . . . . . 4-147
4.2.1.1 IEQUIL for Three-Phase . . . . . . . . . . . . . . . . . . . . . . . . . . 4-147
4.2.1.2 IEQUIL for GASWATER Option . . . . . . . . . . . . . . . . . . . 4-148
4.2.1.3 IEQUIL for WATEROIL Option . . . . . . . . . . . . . . . . . . . 4-149
4.2.2 Saturation Pressure Varies with Depth . . . . . . . . . . . . . . . . . . . . . 4-151
4.2.3 Equilibrium for User-Specified Saturations . . . . . . . . . . . . . . . . . 4-151
4.2.4 Saturation Pressures for VIP-ENCORE (BPTAB) . . . . . . . . . . . . . 4-151
4.2.5 Saturation Pressure Variation with Depth for
Modified Black Oil . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-152
4.3 Saturation-Dependent Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-155
4.3.1 Saturation-Dependent Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-157
4.3.1.1 Water-Oil Saturation for the Matrix (SWT) . . . . . . . . . . 4-158
4.3.1.2 Gas-Oil Saturation for the Matrix (SGT) . . . . . . . . . . . . . 4-163
4.3.1.3 Gas-Dependent Water Relative Permeability for the
Matrix (SGWT) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-167
4.3.1.4 Water-Oil Saturation for the Fracture (SWTF) . . . . . . . . 4-168
4.3.1.5 Gas-Oil Saturation for the Fracture (SGTF) . . . . . . . . . . 4-169
4.3.1.6 Gas-Dependent Water Relative Permeability for the
Fracture (SGWTF) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-170
4.3.1.7 Oil Phase Hysteresis Option (SOTR) . . . . . . . . . . . . . . . 4-171
4.3.1.8 Gas Phase Hysteresis Option (SGTR) . . . . . . . . . . . . . . 4-173
4.3.1.9 Gas Remobilization Option (GASRM) . . . . . . . . . . . . . . 4-175
4.3.1.10 Normalized Saturation-Dependent Functions . . . . . . 4-179
4.3.1.11 Vertical Equilibrium Pseudo-Relative
Permeabilities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-179
4.3.2 Saturation-Dependent Functions (VIP-THERM) . . . . . . . . . . . . . 4-180
4.3.3 Temperature-Dependent Endpoints (VIP-THERM) . . . . . . . . . . 4-183
4.3.4 Water-Oil Hysteresis (Reference 32) (VIP-THERM) . . . . . . . . . . 4-185
4.4 Equation of State PVT Property Data (VIP-COMP or VIP-THERM) . . . . 4-188
4.4.1 Reservoir Equation of State (EOS) . . . . . . . . . . . . . . . . . . . . . . . . . 4-189
4.4.2 Component Names (COMPONENTS) . . . . . . . . . . . . . . . . . . . . . 4-191
4.4.3 Component Characteristics (PROPERTIES) . . . . . . . . . . . . . . . . . 4-192
4.4.4 Binary Interaction Coefficients (DJK) . . . . . . . . . . . . . . . . . . . . . . . 4-196
4.4.5 Lohrenz-Bray-Clark Viscosity Coefficients . . . . . . . . . . . . . . . . . . 4-198
4.4.6 HSTAR Card (VIP-THERM) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-198
4.4.7 Separator Equation of State (EOSSEP)
(VIP-COMP or VIP-THERM) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-199
R2003.4 - Landmark ix
VIP-CORE REFERENCE MANUAL
x Landmark - R2003.4
VIP-CORE REFERENCE MANUAL
Chapter 5
Grid Data Arrays
5.1 Array Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-355
5.2 Start of Array Data (ARRAYS) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-357
5.3 Grid Definition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-358
5.3.1 X (R) Non-Corner Point Grid Dimension (DX, DXB, DR, R) . . . 5-358
5.3.2 X Direction Corner Point Location (XCORN) . . . . . . . . . . . . . . . . 5-359
5.3.3 Y (THETA) Non-Corner Point Grid Dimension
(DY, DYB, DTHETA) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-361
5.3.4 Y Direction Corner Point Location (YCORN) . . . . . . . . . . . . . . . . 5-361
5.4 Gross Thickness - Z grid dimension . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-363
5.4.1 Gross Vertical Thickness, Non-Corner Point Grid (DZ) . . . . . . . 5-363
5.4.2 Gross Stratum Thickness, Non-Corner Point Grid (DZB) . . . . . 5-363
5.4.3 Corner Point Gross Vertical Thickness (DZCORN) . . . . . . . . . . . 5-364
5.4.4 Corner Point Gross Stratum Thickness (DZBCOR) . . . . . . . . . . . 5-364
5.4.5 Depth Corner Point Gross Stratum Thickness (DZVCOR) . . . . . 5-365
5.5 Net Thickness - Z Grid Dimension . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-365
5.5.1 Net Vertical Thickness, Non-Corner Point Grid (DZNET) . . . . . 5-365
5.5.2 Net Stratum Thickness, Non-Corner Point Grid (DZBNET) . . . 5-366
5.5.3 Ratio Net Vertical Thickness to Gross Thickness (NETGRS) . . . 5-366
5.5.4 Fracture Block Net to Gross Vertical Thickness Ratio (NETGF) 5-366
5.6 Depth - Non-Corner Point Grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-368
5.6.1 Depth to Top of Gridblock (DEPTH) . . . . . . . . . . . . . . . . . . . . . . . 5-368
5.6.2 Depth to Center of Gridblock (MDEPTH) . . . . . . . . . . . . . . . . . . . 5-368
R2003.4 - Landmark xi
VIP-CORE REFERENCE MANUAL
Chapter 6
Fault Data
6.1 Fault Option . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-445
6.2 Start of Fault Data (FAULTS) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-446
6.3 Standard Fault Data (FX, FR, FY, FTHETA, FXCORN, FYCORN) . . . . . . 6-448
6.4 Specification of Layers to Which a Shift Applies (LAYERS) . . . . . . . . . . . 6-451
6.5 Specification of a Conductive Fault (LEAKY) . . . . . . . . . . . . . . . . . . . . . . 6-452
R2003.4 - Landmark xv
VIP-CORE REFERENCE MANUAL
Chapter 7
Overread Options
7.1 Transmissibility / Pore Volume Modification Options . . . . . . . . . . . . . . . 7-459
7.2 Override Modification (OVER) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7-460
7.3 Override Modification for VIP-DUAL (OVER) . . . . . . . . . . . . . . . . . . . . . . 7-463
7.4 Value Override (VOVER) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7-464
7.5 Value Override for VIP-DUAL (VOVER) . . . . . . . . . . . . . . . . . . . . . . . . . . . 7-467
Chapter 8
Grid Coarsening
8.1 Grid Coarsening (COARSEN) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8-469
Chapter 9
Region Data
9.1 Assign Output Region Names (REGION) . . . . . . . . . . . . . . . . . . . . . . . . . . 9-473
9.2 Assign Output Regions Separator Batteries (REGSEP) . . . . . . . . . . . . . . . 9-474
9.3 Specify Datum Depth Each Output Region (REGDTM) . . . . . . . . . . . . . . 9-475
Chapter 10
Grid Boundary Flux
10.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-477
10.2 Analytical Aquifer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-477
10.2.1 Carter-Tracy Aquifer Influx (INFLUX) . . . . . . . . . . . . . . . . . . . 10-478
10.2.2 Fetkovich Aquifer Influx (INFLUX) . . . . . . . . . . . . . . . . . . . . . . 10-484
10.3 Coarse Grid, Fine Grid Boundary Flux . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-488
10.3.1 Flux Across a Grid Perimeter (FLUX) . . . . . . . . . . . . . . . . . . . . 10-488
Chapter 11
Local Grid Refinement
11.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-497
11.2 Grid Definition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-497
11.3 Grid Refinement (LGR) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-497
11.3.1 Grid Refinement Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-498
11.3.1.1 Cartesian Grid Refinement . . . . . . . . . . . . . . . . . . . . . . 11-498
11.3.1.2 Radial Grid Refinement . . . . . . . . . . . . . . . . . . . . . . . . 11-499
11.4 Array Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-509
11.5 Array Data Propagation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-509
11.6 Array Input Option . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-509
11.7 Saturation and Relative Permeability Endpoint Arrays . . . . . . . . . . . . . 11-510
11.8 Grid Geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-510
11.9 Corner Point Data (CORP) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-510
11.9.1 Modify by a Constant (MODX,MODY,MODZ) . . . . . . . . . . . . 11-511
11.10 Handedness of Coordinates (RIGHTHANDED) . . . . . . . . . . . . . . . . . 11-513
11.11 Transmissibility Calculations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-514
11.11.1 Harmonic Integration (HARTRAN) . . . . . . . . . . . . . . . . . . . . . 11-514
11.11.2 No Integration (NEWTRAN) . . . . . . . . . . . . . . . . . . . . . . . . . . 11-514
11.11.3 Rectangular or Radial ROOT Grid (BLOCKTR) . . . . . . . . . . . 11-514
11.12 Minimum Radius of Radial Refinements (RMIN) . . . . . . . . . . . . . . . . 11-515
11.13 Connection Transmissibility Modification (MULT) . . . . . . . . . . . . . . . 11-515
11.14 Function Input Option (FUNCTION) . . . . . . . . . . . . . . . . . . . . . . . . . . 11-516
11.15 Arbitrary Grid-block Connections (FTRANS) . . . . . . . . . . . . . . . . . . . . 11-516
11.16 Override Modification (OVER) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-516
11.17 Value Override (VOVER) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-517
11.18 Half-Transmissibility Override (TOVER) . . . . . . . . . . . . . . . . . . . . . . . . 11-517
11.19 Pinchout Gridblock Connections (PINCHGRID) . . . . . . . . . . . . . . . . . 11-518
11.20 Pore Volume Cuttoff (TOLPV) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-519
Chapter 12
Tracer Option
12.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-521
12.2 New Input Data for Initialization Module VIP-CORE . . . . . . . . . . . . . . 12-521
12.2.1 Activate Tracer Option (TRACER) . . . . . . . . . . . . . . . . . . . . . . . 12-521
12.2.2 Additional INFLUX Data (INFLUX) . . . . . . . . . . . . . . . . . . . . . 12-522
Chapter 13
Heat Loss Data (VIP-THERM)
13.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13-523
13.2 Gridding of Over/Underburden . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13-523
13.3 Method of Vinsome and Westerveld (Reference 10) . . . . . . . . . . . . . . . 13-523
13.4 Method of Coats, George, Chu, and Marcum (Reference 40) . . . . . . . . 13-524
13.5 Heat Loss Data Specification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13-525
13.5.1 Specified Index Option . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13-525
13.5.2 Automatic Index Option . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13-526
Chapter 14
Parallel Computing
14.1 Automatic Grid Decomposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14-529
14.1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14-529
14.1.2 Domain Decomposition of Cartesian Grids (DECOMP) . . . . . 14-529
Appendix A
References
Keyword Index
Subject Index
Purpose
VIP-CORE® is the initialization module of the VIP-EXECUTIVE® Family
of simulators. 00
The Modules
The VIP-CORE module calculates the initial reservoir conditions for the
following simulation modules: 00
The Chapters
The VIP-CORE input data stream consists of keywords and data values
which invoke the features of the simulator. 00
■ The subsequent chapters describe the initialization data, which for the
most part, are order dependent. (Any restrictions are described in the
appropriate section.)
Keyword Variables
Data entry formats are always shown between horizontal gray bars, as
illustrated above. Keywords used to label the data are denoted by upper-
case letters. Variable names are shown in lower-case. For example, an
actual line of data based on the above format might look like this:
DATE 31 12 1992
xx Landmark - R2003.4
VIP-CORE REFERENCE MANUAL About This Manual
The following examples show the four different ways in which the above
format can be used:
PRINT ARRAYS KX KY KZ
PRINT ARRAYS ALL
PRINT ARRAYS NONE
PRINT ARRAYS EXCEPT KZ
Compatibility
This documentation is compatible with Version 1998 of the VIP-CORE
software. 00
Internal calculations are carried out in customary oil field units, but input
and output can be in either customary units or the International System of
Units (SI) metric system. Throughout this Reference Manual, units are first
listed as the customary units followed by the SI units in parentheses. The
user may choose metric pressure units of kg/cm2 instead of the default
kPa. In this case, whenever the documentation reads kPa, kg/cm2 will be
expected. 00
Related Documentation
The following manuals provide more information related to the material
in this manual. For more information, please consult the appropriate
manual listed below.
00
1
00000Data Overview
1.1 Introduction
VIP-CORE® is the initialization module of the VIP-EXECUTIVE® family
of simulators. It is used to calculate the initial reservoir conditions to be
used by several simulation modules, including VIP-COMP®, VIP-
ENCORE®, VIP-DUAL®, VIP-POLYMER®, and VIP-THERM. For
example, VIP-ENCORE and VIP-DUAL could be combined to simulate a
dual-porosity, dual-permeability, "black-oil" reservoir. If VIP-COMP were
included in the same program, the user could convert to a fully
compositional version of the dual-porosity, dual-permeability model
simply by substituting the compositional specific data for the black-oil
specific data.
VIP-CORE will only accept data for those modules which have been
purchased. They are listed on the first page of the computer output in the
title box.
4. Replace VIP-ENCORE PVT data with either EOS data (Section 4.4) or
Dead Oil PVT data (Section 4.7). If PCHOR was specified in VIP-
COMP EOS data, that column must be removed before the data will be
accepted by VIP-THERM.
Wells are controlled by a variety of options which allow both rate and
pressure constraints. Wells can also be shut in, or recompleted,
automatically. Rates can be adjusted to meet production and injection
targets at any of the target levels. These targets levels include gathering
center, flow station, area, and field.
Separator conditions are also taken into consideration for all of the
simulator modules. In VIP-COMP (or in the VIP-THERM compositional
model), a multi-component model, separators are required to determine
surface production rates. In VIP-ENCORE, separators allow for the
additional flexibility of treating flash separation conditions at the surface
versus the differential calculations which take place in the reservoir.
Internal calculations are carried out in customary oil field units, but input
and output can be in either customary units or the International System of
Units (SI) metric system. Throughout this Reference Manual, units are first
listed as the customary units followed by the SI units in parentheses. The
user may choose metric pressure units of kg/cm2 instead of the default
kPa. In this case, whenever the documentation reads kPa, kg/cm2 will be
expected.
C
INIT
C ------------------------------------------------------------
C IDENTIFICATION DATA
C ------------------------------------------------------------
C
TITLE1
VIP-CORE BATCH DATA INPUT TEMPLATE
TITLE2
THIS LIST INCLUDES ONLY THE MOST FREQUENTLY USED OPTIONS
C
C
C
C
C ************************************************************
C INITIALIZATION DATA CHAPTER 2
C
C ***********************************************************
C
C ------------------------------------------------------------
C UTILITY DATA SECTION 2.2
C ------------------------------------------------------------
DATE DY MO YR
C
C ------------------------------------------------------------
C PRINT CONTROL CHAPTER 3
C ------------------------------------------------------------
PRINT (ALL)(EQUIL)(COMP)(TABLES)(ARRAYS)(COEFS)(INFLUX)
(INIT)(FAULTS)(CORNER)(NONE)
C NOTE: PRINT CARDS ARE NOT REQUIRED; THE DEFAULT IS PRINT ALL
C
C -----------------------------------------------------------
C GRID SYSTEM OPTIONS SECTION 2.2.3
C ------------------------------------------------------------
NX(NR) NY(NTHETA) NZ (RI) NCOMP (NCV)
C
NO.X-DIR NO.Y-DIR NO.Z-DIR (INNER RADIUS)NO. OF
INCREMENTSINCREMENTS INCREMENTS RI≠0 COMPONENTS
C THE PREVIOUS 2 CARDS ARE REPLACED BY 1 DATA CARD
C NOTE: NCV IS REQUIRED IN VIP-THERM ONLY
C
C -----------------------------------------------------------
C PHYSICAL PROPERTY CONSTANTS SECTION 2.2.4
C ------------------------------------------------------------
DWB BWI VW CW CR TRES TS PS
WATER WATER WATER WATER ROCK RES STAND. STAND.
DENSITY F.V.F. VISC COMPR COMPR TEMP TEMP PRESS
C THE PREVIOUS 2 CARDS ARE REPLACED BY 1 DATA CARD
C NOTE: THE ABOVE TABLE FORMAT DIFFERES IN VIP-THERM
C
C -----------------------------------------------------------
C TABLES CHAPTER 4
C ------------------------------------------------------------
TABLES
C
C -----------------------------------------------------------
C EQUILIBRIUM INITIALIZATION TABLE(S) SECTION 4.2
C ------------------------------------------------------------
IEQUIL PINIT DEPTH PCWOC WOC PCGOC GOC PSAT
C
EQUIL INIT PRES CAP WATER CAP GAS INIT
REGION PRES REF PRES OIL PRES OIL SATR’N
NUMBER DEPTH AT WOC CONTACT AT GOC CONTACT PRES
C THE PREVIOUS 3 CARDS ARE REPLACED BY 1 DATA CARD FOR EACH
C REGION
C
C NOTE: THERE ARE TWO INPUT OPTIONS FOR THIS TABLE: (A)
C SATURATION PRESSURE CONSTANT BY REGIONS, (B)
C SATURATION PRESSURE VARIES WITH DEPTH.
C
C ------------------------------------------------------------
C WATER SATURATION TABLE(S) SECTION 4.3
C ------------------------------------------------------------
SWT 1
SW KRW KROW PCWO
SWR 0.0 KROCW PCWOCW
. . . .
SWRO KRWRO 0.0 PCWORO
. . 0. .
SWMX KRWMX 0.0 PCWOMN
C NOTE: IF ONLY PCWOCW AND PCWOMN ARE INPUT, THE MISSING
C VALUES OF PCWO WILL BE GENERATED BY LINEAR
C INTERPOLATION
C ------------------------------------------------------------
C GAS SATURATION TABLE(S) SECTION 4.3
C ------------------------------------------------------------
SGT 1
SG KRG KROG PCGO
SGMN 0.0 KROCW PCGOMN
. 0. . .
SGC 0.0 KROCG PCGOCG
. . . .
SGRO KRGRO 0.0 PCGORO
. . 0. .
SGMX KRGMX 0.0 PCGOMX
C
C NOTE: IF ONLY PCGOMN AND PCGOCW ARE INPUT, THE MISSING
C VALUES OF PCGO WILL BE GENERATED BY LINEAR
C INTERPOLATION
C
C ---------------------------------------------------------------
C PVT PROPERTY DATA SECTION 4.4
C THE FOLLOWING TABLE APPLIES ONLY TO VIP-ENCORE
C ---------------------------------------------------------------
BOTAB 1
DOS WTOS PSAT
DENSITY OF MOLECULAR WEIGHT INITIAL
SATURATED OIL OF SATURATED OIL SATURATION PRESSURE
C THE PREVIOUS 2 CARDS ARE REPLACED BY 1 DATA CARD
PSAT RS BO BG(ZG) GR VO VG
SATURATIONSOLUTION OIL GAS GAS SATURATED GAS
PRESSURE GAS-OIL FORMATION FORMATION GRAVITY OIL VISCOSITY
RATIO VOLUME (GAS VISCOSITY
FACTOR COMPRESSI-
BILITY)
VOLUME
FACTOR
C THE PREVIOUS 7 CARDS ARE REPLACED BY AT LEAST 2 DATA CARDS
PSAT SATURATION PRESSURE VALUES (AT LEAST 1 VALUE)
DP BOFAC VOFAC (BOFAC VOFAC). . . (BOFAC VOFAC)
RELATIVE OIL OIL
PRESSURE FORMATION VISCOSITY
(P - PSATI) VOLUME FACTOR
C
C NOTES: THE PREVIOUS 3 CARDS ARE REPLACED BY AT LEAST 1 DATA
C CARD. A PAIR OF BOFAC, VOFAC VALUES MUST APPEAR FOR
C EACH VALUE OF PSAT ON THE PSAT CARD.
C
SEPTEST IBAT : SEPARATOR BATTERY NUMBER
PVTTABLE 1 : PVT TABLE NUMBER
PSATF BOF
SATURATION PRESSURE OIL FORMATION VOLUME FACTOR
P T GOR BOSTG GR
STAGE STAGE STAGE GAS STAGE OIL STAGE GAS
PRESSURES TEMPERATURES OIL RATIOSVOLUME FACTORS GRAVITIES
C REPLACE THE PREVIOUS 2 CARDS WITH 1 DATA CARD FOR EACH STAGE
C ------------------------------------------------------------
GRID DATA ARRAYS CHAPTER 5
C ------------------------------------------------------------
ARRAYS
DX(DR) INPUT OPTION
VALUES
DY(DTHETA)INPUT OPTION
VALUES
DZ INPUT OPTION
VALUES
DEPTH INPUT OPTION
VALUES
POR INPUT OPTION
VALUES
KX(KR) INPUT OPTION
VALUES
KY(KTHETA)INPUT OPTION
VALUES
KZ INPUT OPTION
VALUES
C
C NOTE: ALTERNATIVELY, THE TRANSMISSIBILITIES TX, TY, AND
C TZ; THE MDEPTH ARRAY; AND THE PV ARRAY MAY BE
C ENTERED INSTEAD OF THE PRECEDING ARRAYS
C
C NOTE: THE INPUT OPTIONS ARE DESCRIBED IN CHAPTER 1
C
C
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C THE FOLLOWING ARRAYS ARE REQUIRED FOR MULTIPLE TABLES
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C
ISAT INPUT OPTION
VALUES
IEQUIL INPUT OPTION
VALUES
C
C
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C FAULT OPTION DATA (NOT REQ’D) CHAPTER 6
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C
C TRANSMISSIBILITY/PORE VOLUME MODIFICATION OPTION DATA (NOT
C REQ’D) CHAPTER 7
C
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C REGION DATA (NOT REQ’D) CHAPTER 8
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C
C
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C INFLUX OPTION DATA (NOT REQ’D) SECTION 9.1
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
STOP
END
The data stream includes both numbers and alpha keywords; the latter are
used to identify subsequent numbers or select program options. Generally,
each new type of data is introduced by an alpha keyword. Secondary
keyword items in [ ] are "required mutually exclusive" while those in ( )
are "optional". The brackets and parenthesis should not be input with the
data, they simply help describe data input options. In the data
descriptions that follow, keywords are shown in upper case letters (they
can be input in upper or lower case). The names of the variables that are
entered as numbers are shown in lower case letters. For example, one of
the utility data cards indicates the date at which the simulation is to begin.
In the data description, this is written as:
This indicates that the data card must contain the alpha label DATE
followed by three numbers that indicate the day, month, and year. An
example of a valid data card follows:
DATE 15 2 1980
None of the various forms may contain an imbedded blank, since the
simulator interprets them as two words rather than one. There is no
distinction between integer and floating point representations of numbers.
(All numbers are decoded as if they were floating point, then are stored as
either integer or floating point variables, depending on use.)
The single "word" 3*12.5 is decoded as 12.5 12.5 12.5 and 2*15 becomes 15
15. On the other hand, 3* 12.5 could not be decoded properly because of
the blank between * and 12.5. In this case, the simulator issues an error
message. If an error occurs in the initialization data, the run stops prior to
initializing.
Any word beginning with a number (or a decimal point) must be a valid
numeric form, or it causes the run to terminate before the first timestep.
Any word beginning with a character other than a number (or a decimal
point) is treated as alpha data. A # sign before a number causes the
number to be interpreted as an alpha string.
Any word beginning with an exclamation point (!) indicates the beginning
of inline comments. All text after the exclamation point is ignored by the
simulator. The following is an example of the use of inline comments:
C comment
Makes a “comment” of the field which follows. The alpha label C must be
the first word on the card and must be followed by a blank.
See also the use of the inline comment character “!,” which is discussed in
the introduction to this section.
INCLUDE file-name
Definition: 00
When the INCLUDE card is encountered in the input file, the named file is
opened and it becomes the current input file. A fatal error occurs if the file
cannot be opened. Reading from the include file stops when either a
physical end-of-file is encountered or an ENDINC card is encountered.
Reading then continues from the previous input file. Include files can be
nested (i.e. contain other INCLUDE cards). However, the nesting level
cannot exceed nine.
1.5.1.3 Stop Reading Data from the Current INCLUDE File (ENDINC)
ENDINC
The ENDINC card indicates the end of data for the current include file.
When the ENDINC card is encountered on an INCLUDE file, the file is
closed and data continues to be read from the previous input file. If
INCLUDE files are nested then the nesting level is decreased by one. The
ENDINC card is optional in that a physical end-of-file also indicates end
of data. A warning message is given when an ENDINC card is
encountered on the primary input file.
Examples: ENDINC
Definitions:
file-name The pathname to the VDB file from which grid data
should be read.
case-name The case name of the data being read from the VDB
file.
Array data must be of class CALC in order to be read from the VDB file. If
the case name is entered, obviously it must exist in the VDB file. If the case
name is not entered, then only one case name may exist in the VDB file.
The file name may be contained in double quotes. The name may contain
blanks; in this case, it must be contained in double quotes.
LIST
LIST and NOLIST cards control printing of card images of the data read.
Until a NOLIST Card is entered, LIST is assumed.
Examples: LIST
NOLIST
Examples: NOLIST
SKIP
A skip card indicates that all subsequent data is ignored until a NOSKIP
card is encountered; i.e., it is as if all the cards were comment lines. The
card images are not printed.
Example: SKIP
.
.
.
NOSKIP
NOSKIP
NLINES nlines
Definition: 00
Examples: NLINES 60
NCOL ncol
Definition: 00
Only the columns 1 to ncol will be processed. Items beyond column ncol
will be ignored (e.g. comments).
Most arrays have one entry per gridblock. However, several of the corner-
point arrays have more than this. In general, the array input options
extend to these arrays in the obvious way. In a few cases, however, input
options behave differently for corner-point data arrays; the differences are
described in this section.
Definitions: 00
CON requires only one value for input. The entire array is constant for all
gridblocks.
array CON
value(s)
Definition: 00
array XVAR
value(s)
array YVAR
value(s)
ZVAR requires NZ data values. For corner-point arrays NZ+1 data values
are required. ZVAR varies the array property in the z direction while
holding it constant in the x and y directions. This is useful for reading data
that varies by layers.
array ZVAR
value(s)
array VALUE
value(s)
or 00
POR VALUE
C LAYER 1
.10 3*.12
.15 .20 .25 .32
.32 .35 2*.36
C LAYER 2
.12 .13 .14 .15
.15 .16 .16 .17
.25 .22 .20 .21 00
The LAYER option may be used with any input data arrays, but it is
particularly useful for the DEPTH, MDEPTH, or ZCORN arrays. With all
but these three arrays, the LAYER option replicates the values provided
for the first layer for each of the other layers of the grid. When used with
non-corner-point arrays or with the ZCORNE, ZCORNW, or ZCORSW
arrays, the LAYER option requires NX*NY values. When used with other
corner-point arrays, the LAYER option requires (NX+1)*(NY+1) data
values.
To use the LAYER option with the ZCORN array, the XCORN and
YCORN arrays and one of the corner-point thickness arrays DZCORN,
DZVCOR, or DZBCOR must have been previously read. In this case the
LAYER option not only adjusts the values of the ZCORN array, but also it
may adjust the values of the XCORN and YCORN arrays depending on
which of the corner-point thickness arrays are read. If the DZCORN array
is read, then no adjustment is made to the XCORN or YCORN arrays; the
LAYER option works in this case as it does for DEPTH or MDEPTH.
However, if either the DZBCOR or DZVCOR arrays are read, then only the
top layer values of the XCORN and YCORN arrays are kept; the values in
the other layers are adjusted so that the dip of the top layer is mimicked in
succeeding layers and the specified thickness data is honored. The
difference between the DZBCOR and DZVCOR is the interpretation of the
thickness: for DZBCOR, thickness is the thickness perpendicular to the
bedding plane, while for DZVCOR it is the depth difference between a
corner point in one layer and the corresponding corner point in the next
layer.
When MOD cards (Section 1.5.4.1) are used with this option only the
specified locations are changed; the depths of the locations not specified
on the MOD cards are not recalculated. For example, if the depths of layer
1 are modified, then the depths of the other layers are not recalculated.
Thus,
array LAYER
value(s)
The depths for layer 2 are calculated using the DZ values from the ZVAR
example.
The DIP option can be used only to enter the ZCORN, DEPTH or
MDEPTH arrays. When used for the ZCORN array, the XCORN and
YCORN arrays and one of the corner-point thickness arrays must have
been previously read. When used with the DEPTH or MDEPTH arrays,
either the DZ or the H array must have been previously read. Using the
DIP option, VIP-CORE calculates the data values for the first layer of
gridblocks based upon the DIP information, then calculates the depths for
the remaining gridblocks based upon the thicknesses using the LAYER
option.
When MOD cards (Section 1.5.4.1) are used with this option only the
specified locations are changed; the depths of the locations not specified
on the MOD cards are not recalculated. For example, if the depths of layer
1 are modified, then the depths of the other layers are not recalculated.
Thus,
Using the DIP option for radial coordinates does NOT result in a cone.
Instead, it results in a dipping cylinder in which elevations at 180 degrees
will decrease moving outward from the well at the same slope that
elevations will increase along the 0 degree plane.
DEPTH DIP
dref xdip ydip
Definitions: 00
DEPTH DIP
dref rdip theta
Definitions: 00
The MULT option allows the user to define a new array from a previously-
defined array. Each element of the new array is a constant multiple of the
previously defined array. If the new array is the ZCORNE, ZCORNW, or
ZCORSW array, and if the previously-defined array is the ZCORN array,
then the program automatically matches up corner points in ZCORN and
the new array. In all other cases, the new array and the previously-defined
array should have the same dimensions.
array MULT
mult array
Definitions: 00
Example: KZ MULT
.1 KX 00
array dir
value(s)
Definitions: 00
Data values are entered as necessary for the array option being used. 00
Several data conventions should be noted. If array data is entered for the
(+) direction and not for the (-) direction, the (-) direction array will default
to the (+) direction data. If neither (+) direction or (-) direction arrays are
entered, both arrays will default to the array of the same type without
direction specification. For example, if the ISAT X+ array is entered, the
values for ISAT X- default to the ISAT X+ values. However, if the ISAT X+
array is not entered, the default values for both the ISAT X+ and ISAT X-
arrays are the values in the ISAT array. The saturation table assignment for
well calculations defaults, as before, to the ISAT array. If a directional
array is entered for either SWR or SWRO, both arrays must be entered for
that direction. The same requirement holds for the SGR and SGRO arrays.
The areal (A) and vertical (V) direction options cannot be mixed with the
(+) and (-) direction options for the same type of array data.
In the data given above, well relative permeability calculations will use
saturation Table 1, interblock flow calculations in the areal direction will
use saturation Table 2, and interblock flow calculations in the vertical
direction will use saturation Table 3. Water-oil saturation function
endpoint values in the x direction are specified, while the equivalent
values for the y direction default from saturation Table 2, and the z
direction default from saturation Table 3.
array LNXVAR
value(s)
Definition: 00
in the NL+1, NL+2, ....,NZ layers are copied from the NL layer. This option
is usually used with the XCORN array for the LINE corner point position
option.
array LNYVAR
value(s)
Definition: 00
array LNZVAR
value(s)
Definition: 00
array LNVAL
value(s)
Definition: 00
array COPY
Definition: 00
This option can only be used with the ZBOT, ZBOTNW, ZBOTNE, or
ZBOTSW array. No data input is required. The depths of the corner points
in the bottom of the every block are set to be equal to the depths of the
corner points in the top of the next block below it. It is required that
ZCORN, or ZCORNW, ZCORNE, ZCORSW arrays must have been
previously read.
is equivalent to 00
ZCORN ZVAR
10 15.3 18
ZBOT ZVAR
15.3 18 00
MOD
i1 i2 j1 j2 k1 k2 #v
Definitions: 00
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
k1 ≤ K ≤ k2
+ add
- subtract
/ divide
* multiply
= equal
MODLYR
i1 i2 j1 j2 #v
Definitions:
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
K=1
+ add
- subtract
/ divide
* multiply
= equal
Example: DZ CON
10
DEPTH LAYER
3*2000 3*2020
Consider the following MOD and MODLYR cards:
MOD
1 3 1 2 1 1 +4
MODLYR
1 3 1 2 +4
VMOD modifies the array data immediately preceding the VMOD card
with an individual value for each changed gridblock. A minimum of two
cards must follow the VMOD card. The first card contains the locations
describing the gridblocks to be changed. The second card contains the
altered values for those gridblocks. A new VMOD card and its
corresponding data cards are read for each different portion of the grid
system being altered.
VMOD
i1 i2 j1 j2 k1 k2 (op)
Values as necessary
Definitions: 00
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
k1 ≤ K ≤ k2
ADD add
SUB subtract
DIV divide
MULT multiply
EQ equal. This is the default.
2 4 2 3 1 1 EQ
0.25 0.26 2*.27
0.29 0.30 00
Example:
POR VALUE
INCLUDE por.dat
CORP VALUE
INCLUDE corp.dat
The options on the MULT card can be specified in any order. One array
name must be included on each card. Other keywords are optional. If
some option is not included the corresponding value defined in the
previous MULT card is used. For the first MULT card the following default
values are used: ALL MINUS MULT.
The cards are order-dependent. The MULT card allows the user to modify
transmissibility multipliers for standard and non-standard connections.
Connections between Block (I, J, K) and Blocks (I-1, J, K), (I+1, J, K), (I, J+1,
K), (I, J-1, K), (I, J, K-1), (I, J, K+1) are defined as standard connections. All
other connections are non-standard. The transmissibility multipliers for
left and right block faces in three coordinate directions can be specified. If
a nonstandard transmissibility multiplier for a block face is set to zero,
then faulted connections are not automatically generated for this face.A
minimum of one card must follow the MULT card.
The MULT cards must appear after all array input and before any fault
data.
STD
MINUS
MULT array NONSTD (operator)
PLUS
ALL
(FNAME fname)
I1 I2 J1 J2 K1 K2 (val)
(Repeat as necessary)
Definitions: 00
array One of the following: TX, TY, TZ, TR, TTHETA, TXF,
TYF, TZF, TRF or TTHETF. In VIP-THERM, TXT0,
TYT0, TZT0, TRT0, TTT0 are also allowed.
ADD - add
SUB - subtract
DIV - divide
MULT - multiply
EQ - equal.
Example: 00
C
C TRANSMISSIBILITY MULTIPLIERS FOR NON-STANDARD
C CONNECTIONS
C
MULT TX PLUS NONSTD MULT
1 1 3 3 1 1 0.5
MULT TX MINUS DIV
2 2 3 3 2 2
2
2 2 4 5 2 2 3
MULT TY
2 2 4 4 1 1 2
C
CTRANSMISSIBILITY MULTIPLIERS FOR STANDARD
CCONNECTIONS
C
MULT TZ STD PLUS MULT
1 10 2 2 2 2
1 2 3
4 5 6
7 8 9 10
MULT TX
2 3 3 3 2 2
2 0.5
MULT TY MINUS ADD
2 2 3 3 2 2 3
MULT TZ DIV
2 2 3 3 2 3
2
4
NOTE: 1. The PLUS and MINUS options can be mixed only if LGR is on,
or if STD and the operators are MULT or DIV.
MULTIR
itr1 itr2 tmul (X) (Y) (Z) (STD) (NONSTD)
(Repeat as necessary)
Definition: 00
Note: All the transmissibility multipliers are cumulative. They are applied
on top of any previously defined using other options. 00
Definition: 00
2
00000Initialization Data
2.1 Introduction
Initialization data include all data defined at time zero in a simulation. These
include: 00
■ title cards
■ fluid properties.
The data are input after the INIT card and are terminated by the END card. The
data are checked for consistency as they are read. If no errors occur, the data are
then processed. If no initialization errors are generated, the program initializes
and writes an initialization restart file. When the program initializes, the average
pressure, fluids in place, saturation pressure, oil, water, and gas saturation, and oil
and gas composition arrays are calculated. 00
All of the data on each card are free field, following the description given in
Section 1.5. However, the order of data cards is important. The first card in the
initialization data must be the INIT card. 00
The data groups must be arranged as shown in Figure 2-1. Some restrictions on
the order of data cards within each data group may apply. These restrictions are
described in subsequent sections of this chapter. 00
In all of the data descriptions that follow, parentheses are used to indicate optional
items of data. Parentheses are never included in the actual data stream. Items of
data that are aligned vertically in the description of a single data card indicate a
choice; these items are mutually exclusive. 00
00
2.2.1 General
INIT
The first card in the initialization data contains only the keyword INIT.
Initialization data are terminated by the appearance of an END card. 00
Example: 00
INIT
TITLE1
HISTORY MATCH NO FAULTS ZEROED OUT
TITLE2
SAG12 WITH RELICT OIL
DATE 1 6 1977 00
The DIM card allows the user to change the default dimensions on any
initialization run. Multiple sets of DIM cards may be entered, one after the other. 00
Definitions: 00
NOTE: If data is read for oil only or gas only, the minimum required is 2 curves. If
data is read for both oil and gas, the minimum required is 4 curves.
◆ Dual Porosity/Permeability
NBMAX = NBMAX * 2
Example: 00
TITLE cards contain descriptive information about the run that will be printed in
the output title blocks. A maximum of three titles can be read, each on a separate
card following its corresponding alpha label card (TITLE1, TITLE2, or TITLE3).
Title cards defined during initialization are used throughout the simulation unless
redefined in a restart run (see Section 2.1 of the Simulation Modules Manual).
Title cards are not required. Any or all of the three title cards contained on a
restart record can be overwritten by entering new title cards. 00
TITLE1
title
TITLE2
title
TITLE3
title
Definition: 00
The initialization date is input using a DATE card. Time is initialized to zero at
this date. VIP-CORE accounts for leap years, and provides for proper operation
into the 21st century. 00
DATE day mo yr
Definitions: 00
This primary keyword appears in both the Initialization and Simulation Modules 00
END
The END card is required and must be the last card in the data stream; it acts as an
end-of-file marker. 00
The MAP card causes initialization arrays to be written to the vdb file or to the
map file. 00
One of the map cards MAP, MAPX, MAPY, or MAPZ must be entered if
recurrent arrays are to be mapped from the simulation module. 00
Definition: 00
POR Porosity.
KX X-direction permeability.
KY Y-direction permeability.
KZ Z-direction permeability.
P Initial pressure.
PV Pore volume.
TX X-direction transmissibility.
TY Y-direction transmissibility.
TZ Z-direction transmissibility.
KH Permeability thickness.
PD Dilation pressure.
PR Recompaction pressure.
T Temperature.
NOTE: 1. The pore volume array and the corner-point array are always mapped.
2. If no array names are entered on the MAP card (i.e, MAP or MAP
FORM is input), the three saturation arrays and the x-direction, y-
direction, and z-direction transmissibilty arrays will by default be
mapped.
5. Writing map arrays to the vdb file is the default. The NOVDB card
(Section 2.2.2.4) must be entered to write arrays to the map file.
The MAPX, MAPY, and MAPZ cards cause the appropriate initialization mole
fraction arrays to be written to the vdb file or to the map file. 00
One of the map cards MAP, MAPX, MAPY, or MAPZ must be entered if
recurrent arrays are to be mapped from the simulation module. 00
Definitions: 00
NOTE: 1. If the user wishes to map arrays of mole fractions, the components to
be mapped must be explicitly selected on one of the above cards.
3. Writing map arrays to the vdb file is the default. The NOVDB card
(Section 2.2.2.4) must be entered to write arrays to the map file.
The NOVDBPACK card is used to cause full size gridblock arrays to be written to
the vdb file in both VIP-CORE and the simulator module. The default is to
write array data for only the active gridblocks. 00
NOVDBPACK
2.2.2.4 Map File Instead of VDB File (NOVDB) (VIP-COMP and VIP-
ENCORE)
The NOVDB card is used to cause, in both VIP-CORE and the simulation
module, map arrays to be written to the map file rather than the vdb file. Writing
to the vdb file is the default. 00
NOVDB
Arrays 00
T TXT TYT TZT TRT TTT DENW VISW YW HLOS CHLS SOR HOIL HGAS
HWAT αFβγ (VIP-THERM only) 00
TDL TDR PVML TXML TYML TZML TRML TTML (VIP-THERM executable
only) 00
Definitions: 00
P Pressure.
SG Gas saturation.
SW Water saturation.
SO Oil saturation.
WCUT Water-cut.
API API gravity of the liquid phase, available only if the PVT
interpolation option is in use.
The following arrays are available only for the polymer injection option (VIP-
POLYMER): 00
The following arrays are available only in the thermal option (VIP-THERM): 00
T Temperature.
NOTE: The MAPOLD card is needed if records of grid data, in the SIMOUT
form, are to be written during the simulation for subsequent post-
processing. There are two formats. If none of the array names are
specified, all recurrent output arrays are written to the map file at
WMAPOLD frequency (see Simulation Module). If the array names are
specified, only those arrays are written to the map file.
For additional information, see the SIMOUT Map file format description
in the SIMULATION MODULE, OUTPUT CONTROL section.
Example: 00
MAPOLD FORM P SG SW SO 00
NX NY NZ NCOMP
nx ny nz nc
For VIP-THERM, 00
NX NY NZ NCOMP NCV
nx ny nz nc ncv
Definitions: 00
NOTE: For the PVT interpolation option (Section 2.2.11.4) the number of
components must be set equal to 3.
Example: 00
NX NY NZ NCOMP
100 1 25 2 00
2.2.3.2 Radial (Cylindrical) Grid System (NR, NTHETA, NZ, RI, NCOMP)
NR NTHETA NZ RI NCOMP
nr ntheta nz ri nc
For VIP-THERM, 00
Definitions: 00
NOTE: 1. For the PVT interpolation option (Section 2.2.11.4) the number of
components must be set equal to 3.
2. For the VIP-THERM dead oil option, nc must be set to 1 and ncv
must be set to 0.
Example: 00
NR NTHETA NZ RI NCOMP
37 4 3 .5 2 00
WBSIM
This option is used to invoke a special radial initialization for gridded wellbore
simulation in which the first ring of gridblocks is used to model the wellbore. The
equations for wellbore flow are transformed into a form similar to Darcy flow, and
then applied to the vertical flow within the first ring of the grid system. Numerous
correlations are available for vertical or inclined pipe flow, including Hagedorn
and Brown, Dunns and Ross, Beggs and Brill, Aziz and Govier, Orkiszewski, and
Griffith, Lau, Hon, and Pearson. In these correlations, flow conditions are divided
into patterns or flow regimes. Using the user-specified flow correlation and the
computed flow regime at each interval, the simulator transforms the wellbore flow
equations into Darcy-type vertical flow coefficients for each interval at the start of
each timestep.
When the key word WBSIM is specified, the first column of gridblocks is
initialized to a porosity of 1.0 and the connate water saturation is set to 0.0, and it
is sealed off from the reservoir. All other input data for this first column
(wellbore) except for depth and gross thickness will not be used, except if a zero
vertical permeability is specified for any gridblock in the first column, the
wellbore will be sealed at that point.
The RI input variable, which is normally the wellbore radius, should be set to 0.0
for the case of simulating the wellbore as the first column of gridblocks. If the first
column of gridblocks is to be used for the annulus, then RI should be the outside
diameter of the tubing, and the first radius should be the inside diameter of the
casing.
The first radius specified will be the radius of the wellbore, followed by the
normal progression of radii to define the areal extent of the drainage area for the
well.
A separate equilibrium region should be defined for the wellbore [first column of
gridblocks] using the PVT data from the initially producing horizon with the
highest common datum pressure, but lowering the water level to below the deepest
productive horizon. This ensures that the wellbore is in pressure equilibrium with
and contains the fluid from the initial producing horizon.
NOTE: This option is compatible only with non-corner point rectangular grids.
Data given in this section must immediately follow that given in Sections
2.2.3.1 or 2.2.3.2.
Automatic grid setup is available for 1/8 and 1/4 elements of symmetry of 5- and
9-spot patterns and 1/12 and 1/6 elements of symmetry of 7-spot patterns. Both
diagonal and parallel grids are available for the 5- and 9-spot patterns. Use of the
nine-point finite difference option (NINEPT card, Section 2.2.6.5) is
recommended to reduce grid orientation effects. Patterns are allowed to be non-
square (rectangular) only for 1/8 or 1/4 elements of 5- or 9-spot patterns with
parallel grids. 00
00
A diagonal grid for the 1/8 element of symmetry is shown in Figure 2-3. This grid
is called diagonal because the flow from the injector to the corner producer is in
the diagonal grid direction. The grid is drawn such that the standard well locations
are in the center of their respective gridblocks and such that the lines of symmetry
(lines, AB, AD and BD) cut the gridblocks on the boundary by exactly half. This
introduces restrictions upon the grid dimensions. In this case NX must equal NY,
and for all types of pattern elements DX and DY must be constant. 00
In order to properly set up the simulation grid for the pattern element shown in
Figure 2-3, the following must be provided for and are accomplished
automatically with this option: 00
In addition, all total field, well management level, and regional extensive output
data (production, injection, fluids in place, etc.) are scaled to (an optional multiple
of) the full pattern value. All extensive output data for individual wells are scaled
to the full well value. 00
The following sections describe the generated grid, required data, and restrictions
for each pattern element option. 00
A 5-spot pattern is the same as the 9-spot pattern shown in Figure 2-2 except that
no wells exist at points B. For the diagonal grid only, the well spacing is allowed
to differ in the xp and yp pattern coordinate directions. A well spacing ratio (wsr)
is defined as L/W in Figure 2-2. 00
The diagonal grid for 1/8 of a 5- or 9-spot pattern is shown in Figure 2-3. The
parallel grid is shown in Figure 2-4. The following data must immediately follow
the grid dimension data in Section 2.2.3.1 or 2.2.3.2. 00
5
PATTERN (np) EIGHTH (DIAGONAL) area (wsr)
9
Definitions: 00
Restrictions: 00
Diagonal Grid
NX = NY ≥ 2
Parallel Grid
NX = 2 * NY - 1 ≥ 3
wsr may not be specified
NOTE: 1. All field and regional extensive output data are scaled to values for np
full patterns.
2. All extensive well data, such as maximum rates, are input on a full-
well basis.
3. All extensive output data for individual wells are scaled to full well
values.
00
The diagonal grid for 1/4 of a 5- or 9-spot pattern is shown in Figure 2-5. The
parallel grid is shown in Figure 2-6. For the diagonal grid only, well spacing is
allowed to differ in the xp and yp pattern coordinate directions. The following data
must immediately follow the grid dimension data given in Section 2.2.3.1 or
2.2.3.2. 00
5
PATTERN (np) FOURTH (DIAGONAL) area (wsr)
9
Definitions: 00
Restrictions: 00
Diagonal Grid
NX = NY ≥ 2
Parallel Grid
NX = NY = odd number ≥ 3
wsr may not be specified
NOTE: 1. All field and regional extensive output data are scaled to values for np
full patterns.
2. All extensive well data, such as maximum rates, are input on a full-
well basis.
3. All extensive output data for individual wells are scaled to full well
values.
The 7-spot pattern is in the shape of a regular hexagon and is shown in Figure 2-7.
Only one type of grid is available for 1/12 of a 7-spot pattern. This grid is parallel
with respect to the pattern coordinate directions but diagonal with respect to inter-
well flow, and is shown in Figure 2-8. 00
00
The following data must immediately follow the grid dimension data given in
Section 2.2.3.1 or 2.2.3.2. 00
Definitions: 00
7 Designates 7-spot.
Restrictions:
NX = NY ≥ 2 00
NOTE: 1. All field and regional extensive output data are scaled to values for np
full patterns.
2. All extensive well data, such as maximum rates, are input on a full-
well basis.
3. All extensive output data for individual wells are scaled to full well
values.
For 1/6 of a 7-spot pattern, the parallel and diagonal grids are identical since the
wells at points B and C in Figure 2-7 are identical. The grid is shown in Figure 2-
9. 00
The following data must immediately follow the grid dimension data given in
Section 2.2.3.1 or 2.2.3.2. 00
Definitions: 00
7 Designates 7-spot.
Restrictions:
NX = 2 * NY -1 ≥ 3 00
NOTE: 1. All field, regional, and well management level extensive output data
are scaled to values for np full patterns.
2. All extensive well data, such as maximum rates, are input on a full-
well basis.
3. All extensive output data for individual wells are scaled to full well
values.
The constant properties of the reservoir and its fluids are input in one of the
following forms: 00
b. This form assumes that water properties will be specified with either
PVTW (Section 4.11.1) or PVTWSAL (Section 4.11.2).
Definitions: 00
2.2.4.2 VIP-THERM
3. Specified by gridblock using the TEMP array (Section 5.18), overriding options 1 and
2.
If the second or third option is used, some value for TRES must be input in the
physical property constant data, although it will not be used.
Definitions:
KT = KT0 (1 - DKTDSG*SG)
where KT0 is the reference liquid filled thermal conductivity. The derivative term,
DKTDSG, is therefore defined as
– δ KT 1
DKTDSG = -------------- -----------
δSG KT 0
The partial derivatives appearing in the above equations are assumed to be
constants.
Example:
where:
PVMUL = Optional pore volume multiplier given as a function of
pressure (Compaction Option, Section 4.12),
PVREF = input pore volume, or input porosity times total (constant)
block volume,
CR = input rock compressibility,
Porosity is given by the same expressions, with PHI (true porosity) and POR
(input, or reference, value of porosity) replacing PV and PVREF. Note that this
value of porosity is not the true rock porosity, but the ‘simulator’ porosity
corresponding to the assumption of constant total block volume, such that the pore
volume expression predicts the true pore volume behavior as a function of
pressure. 00
PVEXP
PVLINEAR
Definitions: 00
2.2.5 Output
The METRIC card indicates that all data being read and all printed output are in
metric units.
KPA
METRIC
KG ⁄ CM2
BAR
The pressure units for the metric option may be kPa, bar or kg/cm2; if neither KG/
CM2 nor BAR is specified, the default is kPa.
Example:
METRIC KG/CM2
The LAB card indicates that all data being read and all printed output are in
laboratory units. Please refer to VIP-EXECUTIVE Technical Reference, Chapter
31.
LAB
When laboratory units are used, the program automatically adjusts some of the
internal tolerances used by the well calculations for testing of shut-in wells.
CROSS
REGNZ
Initialization region reports usually provide totals for each region, but REGNZ
causes layer values to be output, as well. 00
Definitions: 00
2.2.6 Formulation
GASWATER
This option is used to invoke the two phase gas-water option. This option is useful
for cases with no oil or condensate in the reservoir (dry gas). The gas phase has
constant composition and hence only one hydrocarbon equation is solved in
IMPLICIT cases - coefficient generation is unaffected. This option provides
reduced data input of the saturation and PVT tables, where only the gas properties
are required. 00
NOTE: 1. Simplified equilibrium data for this option is possible (see Section
4.2.1.2).
2. The saturation tables are specified using a modified SGT card, and
SWT data is not required (as described in Section 4.3.1.2).
3. The PVT data is specified by the BGTAB card, rather than the
BOTAB card (as described in Section 4.5.3).
WATEROIL
This option is used to invoke the two phase water-oil option. This option is useful
for cases with no gas in the reservoir. The oil phase has constant composition and
hence only one hydrocarbon equation is solved in IMPLICIT cases - coefficient
generation is unaffected. This option provides reduced data input of the saturation
and PVT tables, where only the oil properties are required. 00
NOTE: 1. Simplified equilibrium data for this option is possible (see Section
4.2.1.3).
3. The PVT data is specified by the BOTAB card, but gas parameters
may be omitted (as described in Section 4.5.2.2).
BLACKOIL
The black-oil option allows the simulation of a true black-oil system in which the
K-value of the oil component is near zero (on the order of 10-10). The option will
prevent the residual oil (or the relict oil in a gas cap) from vaporizing during gas
injection. In the standard black-oil model, the K-value of the oil component is on
the order of 10-3 to 10-2 even if the oil is specified as non-volatile, i.e., input of a
constant GR (gas gravity) in the BOTAB table. In addition to this keyword, the
other data requirements are noted below. 00
NOTE: 1. All input GR values in the BOTAB table must be the same.
TWOPT
Two point upstream has a stability limit only two-thirds the size of the limit of the
default single-point upstream method. Care must be taken when timestep size is
determined. When using the IMPSTAB option to control timestep size,
parameters "ststar" and "stslim" must be explicitly defined (on IMPSTAB card)
as 0.6 and 0.66 respectively. 00
NINEPT
1. Cartesian grid system or LGR model with Cartesian root grid and Cartesian
refinements only.
2. NX > 1, NY > 1.
The NINEPT card causes the program to utilize nine-point finite difference
approximations on the X-Y plane, while the Z-direction (if present) is discretized
using the usual three-point approximation. Unless input, diagonal
transmissibilities are internally calculated using the method of Coats and Modine
(Reference 3). For models with local grid refinements, the diagonal
transmissibilities are computed only inside each grid but not across grid
boundaries. 00
When using the NINEPT option in cartesian or corner point grids, the areal outer
reservoir boundary may be specified as passing through the outer gridblock
centers rather than along the outer edges of the outer gridblocks. This option is
particularly useful when manually setting up pattern elements of symmetry. The
areal transmissibilities parallel to these boundaries are then automatically
computed using the corrrect boundary conditions (for non-homogenous and/or
non-isotropic grids, the correct transmissibilities cannot be obtained by manual
scaling when using the NINEPT option). The SEBOUND option affects only the
areal flow trans-missibility calculations - scaling of pore volumes, bulk volumes,
vertical transmissibilities, thermal transmissibilities, etc., to account for partial
gridblocks must be performed manually. 00
SEBOUND dir
Definitions: 00
In this option, the residual hydrocarbons will expand (or flash into two phases) or
contract with the pressure change. If the pressure drops below the bubblepoint, the
total molar density of the residual hydrocarbons will be calculated using internally
generated tables. Alternatively, the user may elect to approximate the total molar
density using direct extrapolation of the oil phase molar density at the
bubblepoint. When using internally generated tables, fine grids in LGR models
may not be deactivated and/or activated during the run. 00
Definitions: 00
Example: 00
END2P
In the two-point case, each curve is scaled over its entire length (from residual/
irreducible saturation to the saturation at which it attains a maximum). In the
three-point case, all curves in a table are scaled together (retaining the relative kr
and Pc characteristics of the curves). In this case, which is the default scaling
method, all endpoints serve to break the table up into partitions, with scaling
being done in each partition independently. As an example, the water-oil table
would be scaled in two sections; from water saturations of Swir to Swor and from
Swor to 1.
The first keyword defines the calculation method for 3-phase oil relative
permeability. The second keyword defines the calculation method for 3-phase
water relative permeability when the two-curve water relative permeability option
is selected. The default method for the two-curve water relative permeability
option is the same as the calculation method for oil relative permeability.
STONEm
STONE1 (st1exp (st1eps (st1nu)))
KROINT
STONEm
STONE2
KROINT
Definition:
If not one of the STONE1, STONE2, or KROINT cards is input, the default
method is STONE2.
STONE1
STONE2
Examples: STONE1
This option will cause three phase oil relative permeability to be calculated using
a saturation weighted interpolation method. (see Reference 25).
The calculation model for 3-phase water relative permeability when the two-curve
water relative permeability option is selected is defined in the second keyword.
The default is the same as the calculation method for oil relative permeability. 00
STONEm
KROINT
KROINT
Definitions: 00
The PCHYSW card is used to invoke the water-oil capillary pressure hysteresis
option and to define certain control parameters for use in the simulation modules. 00
Definitions: 00
Examples: 00
The PCHYSG card is used to invoke the gas-oil capillary pressure hysteresis
option and to define certain control parameters for use in the simulation modules. 00
Definitions: 00
minsg The minimum gas saturation value for which the program
will calculate hysteresis. Default is 0.
Examples: 00
The RPHYSO card is used to invoke the oil phase relative permeability hysteresis
option for use in the simulation modules. 00
RPHYSO (meth)
Definitions: 00
The RPHYSG card is used to invoke the gas phase relative permeability hysteresis
option for use in the simulation modules. 00
RPHYSG (meth)
Definitions:
The RPHYST card can be used to modify the tolerances to be used for the oil and
gas relative permeability hysteresis calculations. These tolerances must be
exceeded before the hysteresis calculations are performed; otherwise, only the
drainage relative permeability curves are used. 00
Definitions:
NOTE: For both oil and gas relative permeability hysteresis, the calculated trapped
saturation from Land’s equation is restricted to be no greater than the critical
saturation plus 70 percent of (adjusted historical maximum saturation minus
the critical saturation).
JFUNC
The JFUNC card invokes the option to use the Leverett J-function for the
calculation of capillary pressures. The J-function is input in place of capillary
pressure in the SWT table as a function of water saturation. For each gridblock,
the program calculates the square root of permeability divided by porosity, in
order to scale the J-function to a unique capillary pressure curve. The J-function is
defined as follows: 00
Pc K
J ( Sw ) = ------------------- ---
σ cos θ c φ
Pc wo ( Sw ) = J ( Sw ) σ cos θ c 00
Pc wo ( Sw )
Pc = ------------------------
- 00
Kx
------
φ
The COMPACT card invokes the compaction option. Compaction tables are
described in Section 4.12.
COMPACT (REVERSE)
Definition:
FRZPCW
This option freezes the water-oil capillary pressure at its initial value for the entire
simulation. Applies only to saturation tables, does not apply to saturation-
dependent functions option.
FRZPCG
This option freezes the gas-oil capillary pressure at its initial value for the entire
simulation. Applies only to saturation tables, does not apply to saturation-
dependent functions option.
2.2.8 Initialization
NONEQ
The GBC card invokes the gridblock center initialization algorithm. In this case
the saturation distribution of the gridblock is determined by the fluid located at the
gridblock center. This is the default. 00
GBC
Definitions: 00
The INTSAT card invokes the integrated saturation initialization algorithm. In this
case, initial fluid saturations are calculated based upon the actual fluid distribution
throughout the gridblock, by integrating the capillary pressure on the block
thickness (default is gridblock center initialization, see GBC card). The MOBILE
keyword turns on calculation of modified oil and gas residual values, to account
for initial contacts in some blocks. This is especially useful for grids with large
blocks downdip, near the water-oil contact, where a small oil thickness above the
contact will contain mobile oil, while the oil saturation is much less than Sorw. 00
Examples: 00
INTSAT MOBILE 00
Definitions: 00
The VAITS card invokes the volume averaged integrated saturation initialization
algorithm. In this case, the saturation values are calculated based on the actual
Examples: 00
THCNTR
The THCNTR card uses the thickness along the bedding plane to calculate block
properties. The default is to use true vertical thickness to calculate block
properties; this more accurately represents dipping reservoir blocks. 00
NOINIT
The NOINIT card will cause the program to stop after all the coefficents are
calculated. This card is useful when used in conjunction with the PRINT COEFS
FILE, Section 3.3.2. 00
NOROOT
The keyword NOROOT will reduce memory and disk space by ignoring the
ROOT grid, when LGR’s result in no active grid blocks remaining in the ROOT
grid. The keyword NOROOT should not be used if grids are planed for
deactivation in VIP-EXEC such that grid blocks in the ROOT grid will become
active. 00
DEACTIVATE
The keyword DEACTIVATE is required when LGR’s of the same level (i.e.,
siblings) touch, and are planed to be deactivated separately in VIP-EXEC.
Without this keyword, the transmissibility connections between the LGR’s and the
parent grid at the interface between the sibling grids are not created. 00
Definition: 00
The FAULTS card invokes the fault modeling option. Additional fault data are
described in Section 6.1. 00
LATERAL
FAULTS
NONE
Definition: 00
The FLOW360 card indicates that, for a radial model, flow is permitted between
the first and last gridblocks in the theta direction. There is no internal boundary in
the angular direction. This option requires that the total angular span be 360
degrees. In radial models sweeping less than 360 degrees it is assumed that no
flow occurs across the external boundaries of the first and last gridblocks in the
theta direction. Without a FLOW360 card, a full 360 degree model has a no-flow
boundary at θ = 0 degrees = 360 degrees. The FLOW360 option is subject to the
same formulation choices and solution technique choices as the FAULT option. 00
FLOW360
The vertical equilibrium with completely segregated fluids may be invoked for the
water and oil phases, the gas and oil phases, or both. Initial fluid distributions are
calculated which are consistent with total gravity segregation at the phase contacts
and with the gridblock geometry, if that option is used. Pseudo capillary pressures
are calculated to ensure that the segregated fluids remain in equilibrium. Each
gridblock is divided into several sublayers. Fractional volumes and fractional face
areas of the sublayers are calculated and stored in tables. Oil and gas saturations in
the sublayers are defined to be consistent with the segregated flow assumption and
the sublayer fractional volumes. Pseudo relative permeabilities, krw, krow, krg,
krog, for six faces of each gridblock are defined from saturation values and
fractional face areas of the sublayers. 00
Oil relative permeability, kro, is calculated from the usual three phase
correlations. 00
The degree of vertical equilibrium in each block may be varied by using the
FVEWO and FVEGO array keywords (Section 5.33.1 and Section 5.33.2). 00
VEWO (nvelev)
Definitions: 00
NOTE: If the grid dimension NZ is greater than one and VE corner point geometry
option is off then DRELPM (Section 2.2.10.3) should also be used.
VEGO (nvelev)
Definitions: 00
NOTE: If the grid dimension NZ is greater than one and VE corner point geometry
option is off then DRELPM (Section 2.2.10.3) should also be used.
DRELPM (drpfrc)
Definition: 00
The DRELPM card is used to select the directional relative permeability option
for a vertical equilibrium model. DRELPM turns off the VE corner point
geometry option (nvelev > 1, sections 2.2.10.1 and 2.2.10.2). If the DRELPM card
is omitted in a vertical equilibrium model, and nvelev = 1, the same effective
saturation will be used for vertical flow as for areal flow. If the model does not use
vertical equilibrium, the DRELPM card has no effect. 00
Example: 00
VEITS (nvelev)
Definitions: 00
NOTE: A capillary gravity equilibrium option can be applied only in corner point
geometry models if capillary pressure curves are not flat. This assumes
that phase hydrostatic potential is independent of depth within a cell. But
the assumption about complete fluid segregation is not used. The VAITS
option is applied for the calculation of initial saturation distributions.
Accuracy of the initial saturation calculations can be controlled by the
VAITS card. Each gridblock is divided into several sublayers. Fractional
volumes and fractional face areas of the sublayers are calculated and
stored in tables. Phase pressure and saturations in the sublayers are
accurately defined, taking into account the sublayer fractional volumes,
varying ratios of capillary, viscous and gravity forces at different times
during simulation. Pseudo relative permeabilities, krw, krow, krg, krog, for
the six faces of each gridblock are determined from saturation values and
fractional face areas of the sublayers. Oil relative permeability, kro, is
calculated from the usual three phase correlations.
Example: 00
VEITS 20 00
GIBBS
This card invokes a phase stability test, and calculation of phase equilibrium,
using a GIBBS energy minimization algorithm. The technique provides superior
near critical fluid property values, but at the expense of higher computing cost. It
should be tried in any compositional case where near critical behavior is
suspected, and program performance is unsatisfactory. 00
Definitions: 00
This card invokes the near critical relative permeability and capillary pressure
adjustment option. The adjustments are based on interfacial tension (ift)
calculations - capillary pressure is multiplied by 00
[ift / ift(threshold)]**xex. 00
The IFT option is compatible with both GIBBS and non-GIBBS phase
equilibrium calculation options. Additionally, PCHOR data must be included in
the fluid properties table. 00
Definitions: 00
The NOCHK card can selectively suppress PVT table and saturation table data
checks. By default, all table checks listed above are suppressed if a NOCHK card
is entered. This card should be used cautiously since PVT and saturation table
nonconformity can cause instability in the simulation. 00
Examples: 00
NOCHK PVTTAB 00
The BOTINT card is used to invoke the PVT interpolation option. In this option,
three hydrocarbon components are required, instead of the usual two components
for VIP-ENCORE with black-oil data. PVT data is still specified in terms of two
pseudo-components, just as for ordinary black-oil models, but the PVT data is
parameterized by the API gravity or density of the liquid oil phase at stock tank
conditions. In each PVT region, BOTABS data is specified for several different
values of API gravity (i.e., for several different oil types), and PVT properties are
calculated by interpolation between the specified tables depending on the actual
stock tank API gravity of the insitu oil. Internally, the simulator treats the oil
phase as a mixture of three pseudo-components, and the PVT properties vary
based on the liquid phase split between the heavier two components as well as on
the values of pressure and light component composition. Initially, the split
between the heavier components is chosen to match the API gravity of the stock
tank liquid as specified in a table of API gravity versus depth, so the overall
treatment is equivalent to varying PVT properties with stock tank density. 00
The use of the PVT interpolation option requires changes to some of the following
data cards: NX card (Section 2.2.3.1), NR card (Section 2.2.3.2), PRINT INIT
card (Section 2.2.1.1), MAP card (Section 2.2.2.5), IEQUIL card (Section 4.2.1),
BPTAB card (Section 4.2.4), BOTAB card (Section 4.5.2), SEPTEST data
(Section 4.8.3), and other associated cards. In the simulation modules, the option
may cause a change to the OUTPUT card (Section 2.5.1 of the Simulation
Modules Manual). 00
BOTINT ntypes
Definition: 00
ntypes The number of oil types for which BOTABS data must be
provided in each PVT region. The product of ntypes and
the actual number of PVT regions must be no larger than
the maximum number of PVT tables for which the
program is dimensioned. Default is 1.
Example: 00
BOTINT 2 00
NOTE: The BOTAB tables for each PVT region must be introduced in decreasing
order of API gravity.
The FLASH card allows the user to define the flash calculation method and
maximum allowable iterations for reservoir and surface flashes. The two available
methods are accelerated successive substitution and Newton-Raphson with
successive substitution preconditioning. The default method for reservoir flashes
is accelerated successive substitution. The default method for surface flashes is
Newton-Raphson. 00
RES
FLASH maxss maxnr (ACC)
SURF
Definitions:
NOTE: 1. Only one of the RES or SURF labels may be specified. Separate
FLASH cards can be input for RES and SURF. If neither the RES nor
SURF label is specified, the data are used for both reservoir and
surface flash calculations.
2. If a FLASH card with the RES option is not input, then ACC
(accelerated successive substitution) is the default. If a FLASH card
with the SURF option is not input, accelerated successive substitution
calculations are not performed. When a FLASH card is input, ACC
must be specified to invoke the option.
Example:
FLASH 20 10 SURF
Definition: 00
The CRINIT card is used to invoke the super-critical equilibration option. This
option is used to initialize fluid systems which are always one phase but are gas
condensate overlaying volatile oils. The compositions are such that they are
always above the equilibrium two phase envelope at the reservoir temperature and
pressure. The gas-oil contact is defined by looking at the heaviest component k-
value to insure it is less than one because there is no two phase region in the
reservoir at initial conditions. The value of PINIT at DEPTH from the IEQUIL
card is used as the starting point for the initial pressure profile calculations. 00
If the CRINIT card is used with no arguments, all equilibrium regions will be
initialized using the super-critical option. If the CRINIT keyword is followed by a
list of equilibrium region numbers, the super-critical option will only be applied to
the specified equilibrium regions. 00
DRYGAS
The dry gas option enables the compositional modeling of a hydrocarbon system
that is always in a single phase state. The water phase can co-exist with this single
phase. The hydrocarbon fluid properties are calculated by an equation of state. 00
The dry gas option should only be used if fluids always remain in a single phase
state without crossing any phase boundaries. For example, injection and
production from a gas storage reservoir would be an appropriate situation for this
option, depending on the compositions of the in-situ and injected gas. 00
The dry gas option is required for a fluid system existing at a reservoir
temperature above its cricondentherm. For such a system, no saturation pressure
can be found. This would result in termination of the VIP initialization procedure
without the dry gas option. The dry gas option also allows the simulator to run
considerably faster because phase transition calculations are avoided. 00
Definitions: 00
the vapor phase. The equivalent separate oil and gas data tables are BOOTAB and
BDGTAB. 00
The DRSDT option controls the rate at which the solution gas-oil ratio (Rs) within
a block can increase. The rate at which Rs can decline during depletion is
unaffected. Additionally, the rate at which Rs increases will still be limited by the
availability of free gas, and by the saturated values in the tables. 00
At the extreme where DRSDT has a limit of 0, Rs will stay constant within each
gridblock, and free gas cannot dissolve into the undersaturated oil. At the other
extreme, where DRSDT has a very large limit, Rs will rise to that allowed by the
table or until no free gas remains. This is equivalent to not using the option at all. 00
The use of the FREE keyword limits the application of the DRSDT constraint to
blocks that have free gas. This prevents the limiting of changes in Rs due to the
mixing of fluids from adjacent unsaturated blocks that have different starting
values of Rs. This will allow gas flowing from a two-phase cell into an
undersaturated cell to dissolve into solution in an unrestrained manner. 00
Example: 00
DRSDT 10
C Limits increase of solution GOR to 10 SCF/STB/day
The CORNER card is used to set optional parameters for corner point grids. The
user may set the variables iquads, jquads, and kquads to control the accuracy of
the integration process used for pore volume and transmissibility calculations.
Corner point grids allow the user to describe more general grid systems than those
describable through the specification of block-size arrays only. Some examples of
the use of the corner-point option include the description of nonorthogonal grids
and the modeling of sloping faults and pinchouts. The user should take great care
with the use of this option, since less data checking is performed on the grid
system description. The corner-point option is not compatible with radial or
cylindrical grids. 00
When the FAULTS card (see Section 2.2.9.2) is also specified, the corner-point
option requires somewhat different specification of faults than the standard fault
data description (see Section 6.3). If the user specifies corner-point depths (instead
of explicit faults), the program will automatically generate fault connections. 00
When corner point grids are used, the program will automatically recognize
pinchouts and generate the appropriate nonstandard interlayer gridblock
connections. The FTRANS fault option (see Section 6.6.2) allows the user to
overread the transmissibilities generated in this way. 00
Definitions: 00
Example: 00
CORNER NEWTRAN 3 3 3 00
Definitions: 00
The CORTOL card is used to control under what conditions block connections
will be generated from corner-point data. If a block connection value is less than
one or more of the tolerances, a connection is not calculated. 00
Examples: 00
NOTE: The CORTOL data must be in put in the order shown, but may truncate at
any point after tolt.
The program automatically checks to make sure that each face of each gridblock
projects as a convex quadrilateral on to an average plane "parallel" to the face.
Faces that fail the convexity condition are checked for the angle between each of
the pair of opposite sides. If any angle is larger than 90 degrees a warning message
is printed. The check is performed initially excluding blocks with zero porosity.
After the pore volume has been calculated the check is repeated for blocks with
non-zero pore volume. The CORCHK card allows the user to control the amount
of printout. 00
The program checks that x coordinates of corner points (XCORN array) do not
decrease with increasing i index. The optional keyword NOX can be used to
bypass this check. 00
The program checks that y coordinates of corner points (YCORN array) do not
increase with increasing j index. The optional keyword NOY can be used to
bypass this check. 00
Definition: 00
Example: 00
CORCHK 300 00
This card invokes the dual porosity/permeability option. Two distinct options are
available in VIP-DUAL. By default this option models the flow of fluids in two
continuous media which represent the matrix rock and fractures. Exchange of
fluids between the fractures and matrix rock is based on the Warren and Root
theory and includes the effects of imbibition and gravity drainage. 00
Definition: 00
Example: 00
C INITIAL DATE
DATE 23 6 1983
C DUAL POROSITY
DUAL 00
VIP-DUAL accounts accurately for the effects of imbibition and gravity drainage
from matrix gridblocks that contain many matrix blocks by invoking an automatic
calculation of matrix pseudo capillary pressure. 00
Gas-oil gravity drainage is modeled using either a method which is parallel to that
described for the water-oil imbibitions above, or Coats’ method (Reference 29).
Coats’ method requires additional data as described in Section 4.15. 00
The PSEUDO card allows the user to change the various parameters associated
with these calculations. 00
Definitions: 00
Example: 00
NOTE: The data on the PSEUDO card is order-dependent. Except for the COATS/
NOCOATS keyword, all previous fields must be entered to specify a later
one. The COATS/NOCOATS keyword may be omitted entirely if COATS
is to be used.
The DIFF card is used to enter values of oil phase and gas phase diffusivities.
These are used to compute mass transfer between fracture and matrix via the
mechanism of molecular diffusion. The diffusion calculation is not used for
matrix-matrix or fracture-fracture mass transfer. If the DIFF card is omitted, no
diffusive mass transfer will take place. 00
Definitions: 00
NOTE: The oil and gas phase diffusivities are converted into mass transfer
coefficients using the matrix block shape factors (Section 5.37.2 and
5.37.3) and gridblock dimensions:
Alternatively, TDg and TDo may be entered directly via the TDIFFG and
TDIFFO arrays (Section 5.37.2 and 5.37.3).
Example: 00
DIFF 1 0 00
TRACK nfl
Definition: 00
The TRACK card is used to invoke the fluid tracking option and to define the
number of tracked fluids. 00
Example: 00
TRACK 6 00
The NAMES card is used to assign alphanumeric names to the tracked fluids. 00
Definition: 00
iflnm Tracked fluid name. The first character in the name must
be alphabetic unless the name is immediately preceded
by the character #. Only the first six (6) characters of the
name are retained.
NOTE: The number of iflnm names must equal the number of tracked fluids.
Example: 00
Definitions: 00
nfldwn Number of the tracked fluid assigned to the oil below the
gas-oil contact in transition blocks.
The CONTACT card is used to invoke the special option of assigning oil in
transition blocks to two tracked fluids. 00
Example: 00
If the RELICT and BLACK Oil are tracked fluids 1 and 2 respectively,
and GOC is near the bottom of the grid block: 00
CONTACT 1 2 0.077 00
The TRACKW card is used to invoke the water tracking option and to define the
number of tracked water types. This feature is available with the IMPES
formulation option only. 00
Definition: 00
nsitu Index to the water type to be used for tracking the insitu
water.
The NAMESW card is used to assign alphanumeric names to the tracked water
types. 00
Definition: 00
NOTE: The number of tknamw names must equal the number of tracked water
types.
Definitions: 00
beta The water blocking factor. The value must be greater than
zero.
Input for the miscible option, in addition to the MIS and the ALPHA cards, is
discussed in the Tabular Data section and the Matrix Grid Data Array section. The
miscible option contains both a three-component (the original Todd and
Longstaff) option and a four-component option. To select the three component
(water, oil and solvent) option, two hydrocarbon components must be selected on
the grid dimension card (NX). To select the four component (water, oil, gas and
solvent) option, three hydrocarbon components must be selected. The miscible
option is compatible with both the IMPES formulation and the implicit
formulation. 00
or 00
ALPHA (gridname1)...(gridnamei)dpmis
Definitions: 00
The user can turn on the miscible transition zone option by entering an ALPHA
card. If this card is omitted, a first contact miscible process will be represented. In
a first contact miscible process, to avoid the solution gas and oil being treated as
miscible when the solvent is not present (i.e., during depletion or waterflood), an
ALPHA card with zero values for pmis1 and pmis2 and a non-zero value for
ssmin must be input. 00
There are no default values for pmis1, pmis2, ssmin, and dpmis. The employment
of Format (1) automatically invokes a miscible-immiscible transition where the
interpolation parameter increases linearly between the lower- and the upper-
bound miscibility pressures, provided that the solvent pseudo-saturation is larger
than the value of ssmin. Using the IMPES formulation, the value of pmis1 must
be less than or equal to pmis2, and the user may specify a step function for
interpolation by entering identical values for pmis1 and pmis2. However, if the
fully implicit formulation is to be used, the value of pmis2 must be greater than
pmis1, i.e., a step function for interpolation (with pmis1 = pmis2) is not allowed.
Todd and Longstaff suggest that a value of 0.01 for the parameter ssmin is
sufficient to account for the loss of miscibility. 00
CREEP REVERSE
Definition: 00
The grid description must be defined with the DX, DY, DZ, KX arrays, and the
simulator must be run in the fully implicit, single porosity mode. 00
The required data for this option are HYDFRAC and HYDBETA (Section 4.16),
along with appropriate compaction (CMT) tables to define the fracture closure
effects and permeability reduction as a function of pressure. 00
HYDFRAC nf
Definition: 00
NOTE: The DY’s for Y=1 should be adjusted to define the tapering half-fracture
widths.
The porosities for gridblocks nf+1 through NX (for Y=1) should be set to
zero.
Fracture closure stress can be simulated using the compaction option, with
a separate table(s) for the fracture gridblocks. The pore volume reduction
factor will also be used to reduce the cross-sectional area in the flow
velocity calculations. The transmissibility reduction factor will also be
used to reduce the original KX for the Beta (turbulence) factor
calculations.
POLYMER
The POLYMER card will initialize the model for polymer injection. By default
the model does not initialize for polymer injection. Initializing for polymer
injection will increase the size of the restart files. 00
The physical properties required to model the polymer transport can be introduced
at any point of the recurrent data in the simulation module, preceding polymer
injection. 00
THERMAL
■ Water mole fraction in the vapor phase yw is given by VPW/P where VPW is
the vapor pressure of water as a function of temperature.
In the default method, water is treated as a component in the vapor phase by the
equation of state. Water EOS parameters are internally fixed and were adjusted to
match vapor pressure, steam z-factor, and steam enthalpy along the saturation
curve. Liquid water fugacity coefficient is obtained from an internal table which
contains EOS-predicted values which were adjusted to match the vapor values
along the saturation curve. All other liquid water properties are obtained from
steam tables. Saturation pressure and flash calculations are performed using three
phase algorithms. If k-value options are used, input k-values are assumed to be
three-phase values and are evaluated at total system pressure.
WATIDEAL
The FLOWS card is required in order to map, print, or write to the spreadsheet
file (in VIP-EXECUTIVE) flow rate arrays or region boundary flow rates or
cumulatives. If all arrays are mapped, flow rate arrays will be mapped only if a
FLOWS card is input. Individual flow rate array names may be specified on the
MAP card only if a FLOWS card is input.
FLOWS
The formulation of the velocity dependent models have been developed in the
Department of Petroleum Engineering at Heriot-Watt University under the
directorship of Professor A. Danesh and Professor D. H. Tehrani in a research
project sponsored by the UK Department of Trade and Industry and 11 oil and gas
companies. 00
Velocity dependent effects may be applied to both the production well model and
for flow between gridblocks. 00
Parachor values must be input in the equation of state property tables (see section
4.4.3). 00
WELL
CN
VELCTY FIELD SAND
NDARCY CARBON
BOTH
ALL
N cbg
( m o ) ( no ) 00
Vg
N cbg
( m g ) ( ng ) 00
Vg
(NDARCY (B0) (B1) (B2) (B3) (B4) (B5) ) 00
Definitions: 00
RPCNO, RPCNG, and NDARCY are optional input fields that can be used to
override the default parameters. They are discussed below.
RPCNO
RPCNG
Examples 00
VELCTY
VELCTY CN ALL
3. Apply both capillary number effects and non-Darcy effects at the well
only, but using carbonate default parameters. Also, selectively modify
some of the non-Darcy exponents, removing the dependency on
interfacial tension.
VELCTY CARBON
NDARCY B1 B3 B5
-0.7 0.7 0.
CNDBNK TRAP 20
RSM
3
00000Print Control
3.1 Introduction
Because of the large volume of data that can be generated, PRINT cards
should be used judiciously.
The processed initialization input data are divided into the following
groups: equilibrium data, compositional data, tables, input arrays,
coefficient arrays, fault descriptions, influx data, initialization arrays,
corner-point arrays, separation data, and region data.
3.2 General
PRINT ALL
Definition:
PRINT NONE
Definition:
To continue the table or array names on a PRINT card, list the remaining
table or array names on successive cards. Do not repeat the PRINT label.
Vertically aligned data indicate a choice; these data are mutually exclusive.
The ARRAYS option controls the printing of the input arrays of the grid
data (Section 5.1) as specified, after any modifications have been applied.
If no PRINT ARRAYS card is read, none of the arrays that were read are
printed. Array names can be used with the E as described below. If
EXCEPT is omitted, only the arrays listed are printed. The array names
can appear in any order. The parentheses indicate that the labels within
them are optional.
ALL
PRINT ARRAYS
NONE
EXCEPT array 1 (array 2 ...)
Definitions:
NONE None of the input array data is printed. This has the
same result as not reading a PRINT ARRAYS card.
Example:
PRINT ARRAYS DX DY KX
ALL
PRINT COEFS NONE
EXCEPT array 1 (array 2 ...)
FILE
Arrays
MDEPF PVF TEX TXYLF TXYRF TDIFFG TDIFFO TXF TYF TRF TZF
TTHETF (Dual Porosity Only)
Definitions:
Examples:
PRINT COEFS PV TX TY TZ
ALL
PRINT COMP
NONE
Definitions:
Example:
The CORNER option controls the printing of corner-point data and the
writing of a formatted corner-point output file (FORTRAN Unit 12). The
parentheses indicate that the labels within them are optional.
( PRINT ) ( FILE )
PRINT CORNER ALL
NONE
Definitions:
Example:
ALL
PRINT EQUIL
NONE
Definitions:
Example:
The FAULTS option controls the printing of fault data (Section 6.1) in
tabular form, array form, or both. This card also controls the printing of
any PINCHOUT (Section 2.2.9.1) warning messages. If no PRINT FAULTS
card is read, none of the data are printed. If EXCEPT is omitted, the data
are printed only in the form(s) specified. The parentheses indicate that the
labels within them are optional.
( TABLES ) (ARRAYS)
Definitions:
ALL The fault data are printed in both tabular form and
array form.
NONE None of the fault data are printed. This has the
same result as not reading a PRINT FAULTS card.
Example:
The INFLUX option controls the printing of the tabular and/or array
influx data (Section 10.2) as specified. If no PRINT INFLUX card is read,
none of the influx data are printed. If EXCEPT is omitted, the data are
printed only in the form(s) specified. The parentheses indicate that the
labels within them are optional.
( TABLES ) (ARRAYS)
Definitions:
ALL Both the tabular and array influx data are printed.
NONE None of the influx data are printed. This has the
same result as not reading a PRINT INFLUX card.
Example:
ALL
PRINT INIT
NONE
EXCEPT table 1 (table 2 ...)
Arrays
PF PSATF SOF SWF SGF DENOF DENGF VISOF VISGF VFF XIF YIF ZIF
PVF KHF PDATF SOMF SWMF PCSWF PCSGF IFTF TXF TYF TRF
TTHETF TZF TEX TXYLF TXYRF SGMF PWSATF WCO2F (Fractured
Cases Only)
Definitions:
P Initial pressure.
TEMP Temperature.
Examples:
ALL
PRINT REGION
NONE
Example:
ALL
PRINT SEPARATOR
NONE
Example:
ALL
PRINT TABLES
NONE
EXCEPT table 1 (table 2 ...)
Table Names
Definitions:
NONE None of the tabular data are printed. This has the
same result as not reading a PRINT TABLES card.
Examples:
Endpoint Names:
Definitions:
Examples:
1 .1 .9 .8
Endpoint Names:
Definitions:
4
00000Tables
4.1 Introduction
00 VIP-CORE tabular data must be preceded by the alpha label TABLES.
TABLES
00 The labels on the left side of the figure correspond to specific entries in the
equilibrium initialization table (Section 4.2.1 and Section 4.2.2). The
reference values of capillary pressure are shown in the right hand figure.
VIP-CORE constructs pressure versus depth tables for each phase and
assigns to each gridblock an oil pressure equal to the oil phase pressure at
the depth of the gridblock center. The difference between phase pressures,
capillary pressure, is then used to compute average saturation values for
each gridblock.
Data for each region are input in one of two ways: (1) with a constant
saturation pressure, Section 4.2.1, or (2) with a variable saturation pressure
that is a function of depth, Section 4.2.2. If one equilibrium region has a
constant saturation pressure while another region has a varying saturation
pressure, use the variable saturation pressure format to enter all of the
equilibrium data. 00
IEQUIL PINIT (TINIT) DEPTH PCWOC WOC PCGOC GOC PSAT (API)(IPVT)(IPVTW)
iequil pinit (tinit) depth pcwoc woc pcgoc goc psat (api)(ipvt)(ipvtw)
(one data card for each equilibrium region)
Definitions: 00
The titles on this card must appear in the order shown. The title card is not
repeated for multiple equilibrium regions. 00
If the goc is within the reservoir, then the model calculated psat at goc will
be used for initializing the reservoir pressure, and PINIT will be ignored.
In the VIP-THERM dead oil model, the value of psat is arbitrary and is not
used. 00
Definitions: 00
The titles on this card must appear in the order shown. The title card is not
repeated for multiple equilibrium regions. 00
Definitions: 00
The titles on this card must appear in the order shown. The title card is not
repeated for multiple equilibrium regions. 00
Examples: 00
00 TABLES
C
IEQUIL PINIT DEPTH PCWOC WOC PCGOC GOC PSAT
1 4335 8575 6.0 9022 0.0 8575 4335
00 C
IEQUIL PINIT DEPTH PCGWC GWC
1 5000. 8000. 0.0 8000.
IEQUIL PINIT (TINIT) DEPTH PCWOC WOC PCGOC GOC (IPVT) (IPVTW)
iequil pinit(tinit) depth pcwoc woc pcgoc goc (ipvt) (ipvtw)
(one data card for each equilibrium region)
Data items are defined exactly as described in Section 4.2.1. Note that the
only difference is the absence of the PSAT, and possibly API, data. 00
Depth versus psat tables using the BPTAB card only apply when the PVT
property data are input using the BOTAB card (Section 4.5.2). 00
Saturation pressure tables are read immediately after the cards defining
the fluid contacts for VIP-ENCORE. The values on these cards must
appear in the order shown. One saturation pressure table must be read for
each equilibrium region. 00
Definitions: 00
00 Example:
00 BPTAB 1
DEPTH PSAT
6600 2850
7800 2647
8600 2540
4.2.5 Saturation Pressure Variation with Depth for Modified Black Oil
00 PSTAB or RSTAB and RVTAB cards
These cards may not be specified if PSAT was entered on the IEQUIL card. 00
Depth versus psat or solution oil/gas and solution gas/oil tables using the PSTAB,
RSTAB, or RVTAB cards only apply when the PVT property data are input using
the BOETAB, or some combination of BODTAB, BOOTAB, BOGTAB,
BDGTAB cards. 00
These tables are read immediately after the cards defining the fluid contacts for
VIP-ENCORE. The values on these cards must appear in the order shown. Either
a PSTAB or a RVTAB/RSTAB combination must be read for each equilibrium
region. For a dead oil, only a RVTAB card is required. For a dry gas, only a
RSTAB card is required. 00
00 Definitions:
00 Example:
00 PSTAB 1
DEPTH PSAT
6600 2850
7800 2647
8600 2540
00 Definitions:
(1) RSTAB Indicates the data being read is a table of solution gas-oil
ratio versus depth.
00 Example:
00 RSTAB 1
DEPTH RS
6600 700
7800 750
8600 800
00 Definitions:
(1) RVTAB Indicates the data being read is a table of solution oil-gas
ratio versus depth.
Example:
00 RVTAB 1
DEPTH RV
6600 0.04
7800 0.042
8600 0.045
Stone Model I: 00
00 where
00 S orm = S orw – γ S g ,
S orw – S org
00 γ = -------------------------------- .
1 – S wc – S org
Sg A Sg A
S orm = S org ---------- + S orw 1 – ---------
2 2 3
S gro - – ε ( S gro S g – S g ) + ν ( S gro S g – S g )
S gro
00
k row k rog
k ro = k rocw ------------ + k rw ------------ + k rg – k rg – k rw ,
k rocw k rocw
00
( S g * k rog + ( S w – S wc ) * k row )
K ro = ------------------------------------------------------------------------
( S g + S w – S wc )
00
where: 00
Often the shape of the relative permeability and capillary pressure curves
are identical for several core samples even though the endpoints are not
the same. VIP-CORE allows the user to model spatial differences in curve
endpoints (such as connate water and critical gas saturation) with a single
set of generic relative permeability and capillary pressure curves. 00
1. All water saturation tables must precede the gas saturation tables,
which must precede the gas-water tables, if present.
2. The table numbers for both the water and gas tables must be
consecutively increasing.
3. The number of water saturation tables must equal the number of gas
saturation tables.
concert with the ISAT array for drainage tables). RPHYSO and
RPHYSG cards can optionally be specified to control the method used
to calculate oil and gas relative permeability hysteresis, respectively.
Water saturation tables define the rock properties that depend on water
saturation: relative permeability of water, relative permeability of oil in
the presence of water, and water-oil capillary pressure. If capillary
pressure hysteresis is selected, secondary water-oil capillary pressure
must also be defined for a drainage table. These tables are not needed for
the GASWATER option, Section 2.2.6.1. 00
Definitions: 00
00 (3) The values on the data cards following this title card appear in the
order defined by this card. They must appear in the order shown.
00 KRW = 0.0
pressure at swc).
00 (b) SW = swr
00 KRW = 0.0
00 PCWO = krorw
00 PCWOS = pcworw
00 (c) SW = swd
00 KRW = krwd
00 KROW = krowd
00 PCWO = pcwod
00 PCWOS = pcwod
00 (d) SW = swro
00 KRW = krwro
00 KROW = 0.0
00 PCWO = pcworo
00 PCWOS = pcworos
00 (e) SW = swmx
00 KRW = krwmx
00 KROW = 0.0
00 PCWO = pcwomn
00 PCWOS = pcworos
For PCWO, if only the values pcwocw and pcwomn are read and the
capillary pressure entries between these two values are left blank, the
program linearly interpolates between these two values to complete the
table entries. 00
Example: 00
00 C
C SATURATION DEPENDENT TABLES
C
C
SWT 1
SW KRW KROW PCWO
.05000 0.00000 1.00000 56.30945
.32500 .00005 .64416 39.83473
.36250 .00071 .50332 32.91141
.40000 .00292 .38467 26.80040
.43750 .00758 .29631 21.45391
.47500 .01563 .20633 16.82418
.51250 .02798 .14283 12.86340
.55000 .04555 .09391 9.52381
.58750 .06928 .05768 6.75761
.62500 .10008 .03221 4.51702
.66250 .13888 .01563 2.75426
.70000 .18659 .00601 1.42156
.73750 .24414 .00147 .47111
.77500 .31245 .00009 .17500
.81250 .39245 0.00000 .13750
.85000 .48506 0.00000 .10000
.88750 .59119 0.00000 .06250
.92500 .71178 0.00000 0.02500
1.0000 1.00000 0.00000 0.00000
C
00 Example:
00 C
C SATURATION DEPENDENT TABLES WITH HYSTERESIS
C
00 C
SWT 1
SW KRW KROW PCWO PCWOS
.25000 0.00000 1.0000 56.30945
.32500 .00005 .64416 39.83473
.36250 .00071 .50332 32.91141
.40000 .00292 .38467 26.80040
.43750 .00758 .28631 21.45391
.47500 .01563 .20633 16.82418
.51250 .02798 .14283 12.86340
.55000 .04555 .09391 9.52381
.58750 .06928 .05768 6.75761
.62500 .10008 .03221 4.51702 4.51702
.66250 .13888 .01563 2.75426 2.0
.70000 .18659 .00601 1.42156 -1.0
.73750 .24414 .00147 .47111 -3.0
.77500 .31245 .00009 .17500 -8.5
.81250 .39245 0.00000 .13750 -15.0
.85000 .48506 0.00000 .10000 -15.0
.88750 .59119 0.00000 .06250 -15.0
.92500 .71178 0.00000 0.00000 -15.0
00 Do not use these tables with the WATEROIL option (Section 2.2.6.2).
Gas saturation tables define the rock properties that depend on gas
saturation: relative permeability of gas, relative permeability of oil in the
presence of gas, and gas-oil capillary pressure. The rock properties defined
with the GASWATER option (Section 2.2.6.1) are: relative permeability of
gas; relative permeability of water in the presence of gas; and gas-water
capillary pressure. 00
GASWATER Option 00
SGT nsgt
SG KRG KRWG PCGW
0.0 0.0 1.0 pcgwmn
. . . .
sgc 0.0 krwgc .
. . . .
. . . .
sgcw krgcw 0.0 pcgwcw
Definitions: 00
00 (2) The values on the data cards following this title card appear in the
order defined by this card. They must appear in the order shown.
00 (a) SG = sgmn
00 KRG = 0.0
00 PCGOS = pcgomns
00 (b) SG = sgc
00 KRG = 0.0
00 KROG = krocg
00 PCGO = pcgocg
00 PCGOS = pcgocgs
00 (c) SG = sgd
00 KRG = krgd
00 KROG = krogd
00 PCGO = pcgod
00 PCGOS = pcgod
00 KROG = 0.0
00 PCGO = pcgoro
00 (e) SG = sgmx
00 KRG = krgmx
00 KROG = 0.0
00 PCGO = pcgomx
For PCGO, if only the values pcgomn and pcgomx are read and the
capillary pressure entries between these two values are left blank, the
program linearly interpolates between these two values to complete the
table entries. If capillary pressure hysteresis is selected, secondary gas-oil
capillary pressure must also be defined. 00
00 Examples:
00 C
C GAS SATURATION FUNCTIONS - DRAINAGE
C NOTE: ENTRIES OMITTED IN THE PCGOS COLUMN
C WILL HAVE VALUES SET TO PCGO AT THAT
C GAS SATURATION.
C
SGT 1
SG KRG KROG PCGO PCGOS
0. 0. 1. 0. 0.
.03 .007 .91 .21 0.
.10 .04 .77 .30 0.
00 C
C GAS WATER OPTION
C
SGT 1
SG KRG KRWG PCGW
0.0 0.0 1.0 0.0
0.05 0.0 0.92
0.1 0.0 0.83
0.2 0.1 0.650
0.3 0.2 0.430
0.4 0.37 0.250
0.5 0.45 0.110
0.6 0.57 0.032
0.7 0.75 0.006
0.75 0.85 0.00 4.0
SGWT nsgwt
SG KRWG
0.0 1.0
. .
sgrw 0.0
. .
1.0 0.0
00 Definitions:
00 The values on the data cards following this title card must appear in the
order shown.
A minimum of three data cards, for gas saturations of 0.0, sgrw and 1.0 are
required for each table. SG must increase and the KRWG must be
monotonically decreasing. 00
When the DUAL option is invoked the user must specify fracture
saturation tables in addition to the matrix saturation tables. The water
saturation tables for the fracture must follow the normal (matrix) gas
saturation tables. These tables define the rock properties that depend on
water saturation in the fractures: relative permeability of water, relative
permeability of oil in the presence of water, and water-oil capillary
pressure. 00
SWTF nswtf
(SWMNI swmni)
SW KRW KROW PCWO
swc 0.0 krocw pcwocw
. . . .
swr 0.0 krorw .
. . . .
swro krwro 0.0 .
. . . .
swmx krwmx 0.0 pcwomn
00 Definition:
00 Example:
00 SWTF 1
SW KRW KROW PCWO
0. 0. 1. 0.
.5 .5 .5 0.
1. 1. 0. 0.
When the DUAL option is invoked the user must specify fracture gas
saturation tables. These must follow the fracture water saturation tables.
These tables define the rock properties that depend on gas saturation in
the fractures: relative permeability of gas, relative permeability of oil in
the presence of gas, and gas-oil capillary pressure. 00
SGTF nsgtf
SG KRG KROG PCGO
0.0 0.0 krocw pcgomn
. . . .
sgc 0.0 krocg .
. . . .
. . . .
. . . .
sgro krgro 0.0 .
. . . .
sgmx krgmx 0.0 pcgomx
00 Definitions:
All other definitions are as given in Section 4.3.1.2, Gas Saturation Tables
for the Matrix. 00
00 Examples:
00 SGTF 1
SG KRG KROG PCGO
0. 0. 1. 0.
.5 .5 .5 0.
.9999 .999 0. 0.
1. 1. 0. 0.
SGWTF nsgwtf
SG KRWG
0.0 1.0
. .
sgrw 0.0
. .
1.0 0.0
00 Definitions:
To invoke the oil hysteresis option, the SOTR card must be inserted in the
water-oil saturation table as shown, unless user defined tables are entered
for imbibition functions. All other water saturation data (table entries) are
unchanged. 00
SWT nswt
(SWMNI swmni)
SOTR sotr
SW KRW KROW PCWO PCWOS
swc 0.0 krocw pcwocw
. . . . .
swr 0.0 krorw pcworw
. . . . .
(swd) (krwd) (krowd)(pcwod)(pcwod)
. . . . .
swro krwro 0.0 pcworo pcworos
. . . . .
swmx krwmx 0.0 pcwomn pcworos
00 Definitions:
00
00 krowi(1-So) = krowd(1-Sodr),
To invoke the gas hysteresis option, the SGTR card must be inserted in the
gas saturation table as shown unless user defined tables are entered for
imbibition functions. All other gas saturation data (table entries) are
unchanged. 00
SGT nsgt
SGTR sgtr
SG KRG KROG PCGO (PCGOS)
sgmn 0.0 krocw pcgomn (pcgomns)
. . . . .
sgc 0.0 krocg pcgocg pcgocgs
. . . . .
(sgd) (krgd) (krogd) (pcgod) (pcgod)
. . . . .
sgro krgro 0.0 pcgoro .
. . . . .
sgmx krgmx 0.0 pcgomx .
00 GASWATER Option
SGT nsgt
SGTR sgtr
SG KRG KRWG PCGW
0.0 0.0 1.0 pcgwmn
. . . .
sgc 0.0 krwgc .
. . . .
. . . .
sgcw krgcw 0.0 pcgwcw
00 Definitions:
general, each gridblock may follow its own imbibition curve. Trapped gas
saturation is calculated from Sgmax by using Land’s formula (Reference
6). 00
S gmax * C VIP
S gtr = -------------------------------
S gmax + C VIP
SGTR * (1 – S wc )
C VIP = -------------------------------------------
( 1 – S wc ) – SGTR
krg
Sg
D
krg
I
krg
k rg ( S g ) = k rg ( S gdr ) ,
l D
1
--
1 2 2
S gdr = -- ( S g – S gtr ) + [ ( S g – S gtr ) + 4C VIP ( S g – S gtr ) ]
2
00 Example:
00 C
C GASWATER OPTION
C
SGT 1
SGTR 0.25
SG KRG KRWG PCGW
0.0 0.0 1.0 0.0
0.5 0.0 0.92
0.1 0.0 0.83
0.2 0.1 0.650
0.3 0.2 0.430
0.4 0.37 0.250
0.5 0.45 0.110
0.6 0.57 0.032
0.7 0.75 0.006
0.75 0.85 0.00 4.0
00 Tabular Form:
(4) GASRMT
(5) SG KRG
(6a) sgrmi 0.
. .
(6b) sgmx krgmx
00 Analytical Function:
(4) C1 C2 C3
(5) c1 c2 c3
00 Definitions:
00 (2) The values on the data cards following this title card appear in the
order defined by this card. They must appear in the order shown.
00 (3) Data cards A minimum of two data cards is required for each
saturation table.
00 Tabular Form:
00 (5) The values on the data cards following this title card appear in the
order defined by this card. They must appear in the order shown.
00 (6) Data cards A minimum of two data cards, including each of the
types (a) and (b), is required for each saturation
table. The first data card must be type (a), while the
last must be type (b).
(a) SG = sgrmi
KRG = 0.0
KRG = krgmx
00 Analytic Function:
00 (4) The values on the data card following this title card appear in the
order defined by this card. They must appear in the order shown.
00 (5) Data card Input c1, c2, and c3 are values of C1, C2, and C3.
00 Note:
K rgmax, m o
K rg ( S g ) = ------------------------K rg2d, r ( S g )
krgmx
00 where
o ( sgmx – sgrmi )
S g = sgrmi + ( S g – S grm, m ) ------------------------------------------------
( S gma, m – S grm, m )
00 Here Sgma,m is the gridblock maximum gas saturation in the SGT table.
Symbols Sgrm,m and Krgmax,m are gridblock sgrm and krgmax values
determined from the table lookup of cards (2) and (3) above. The input
sgtr for the table lookup is the gridblock trapped gas saturation
corresponding to the historic maximum gas saturation.
by
C
( 1 + C 2 )S g11
00 K rg ( S g ) = K rgmax, m ----------------------------
C3
-
1 + C 2 S g1
00 where
S g – S grm, m
S g1 ≡ -------------------------------------------
S gma, m – S grm, m
00 The output of the input tables and the reconstructed tables may be
requested through keywords SGT and RSGT in the PRINT card.
00 Example:
C
C Gas Remobilization Option Input
C
GASRM 1
SGTR SGRM KRGMAX
0.05 0.10 0.9
0.10 0.15 0.85
0.17 0.23 0.81
0.24 0.30 0.75
0.33 0.40 0.7
0.4 0.48 0.6
0.5 0.6 0.5
GASRMT
SG KRG
0.21 0.0
0.30 0.02
0.43 0.10
0.55 0.20
0.66 0.33
0.77 0.5
0.985 0.65
C
GASRM 2
SGTR SGRM KRGMAX
0.05 0.10 0.9
0.10 0.15 0.85
0.17 0.23 0.81
0.25 0.31 0.74
0.33 0.40 0.7
0.4 0.48 0.6
0.5 0.6 0.5
GASRMT
SG KRG
0.21 0.0
0.30 0.02
0.43 0.10
0.55 0.20
0.66 0.33
0.78 0.51
0.9 0.65
Often the shape of the relative permeability and capillary pressure curves
are similar for several core samples even though the endpoints are not the
same. VIP-CORE allows the user to model spatial differences in connate
water and critical gas saturations with a single set of generic relative
permeability and capillary pressure curves. 00
SWT nswt
SW KRW KROW PCWO
swc 0.0 krocw 0.0
1-sor krwro 0.0 .
1.0 1.0 0.0 0.0
00 Definitions:
The definitions in Section 4.3.1.1 apply, except only the three data
cards shown are entered.
SGT nsgt
SG KRG KROG PCGO
0.0 0.0 krocw 0.0
sgr 0.0 krogr .
sgro krgro 0.0 .
sgcw krgcw 0.0 0.0
00 Definitions:
The definitions in Section 4.3.1.1 apply, except only the four data cards
shown are entered. 00
ew
Sw – Swir
krw = krwro ----------------------------------------------
1 – Sorw – Swir
eow
1 – Sorw – Sw
krow = krocw --------------------------------------------
1 – Sorw – Swc
eg
Sg – Sgc
krg = krgro ----------------------------------------------------------
1 – Swc – Sorg – Sgc
eog
1 – Swc – Sorg – Sg
krog = krocw -------------------------------------------------------
1 – Swc – Sorg
3
Pcwo = Po – Pw = – ( Pcw1 + Pcw2 ( 1 – Sw ) + Pcw3 ( 1 – Sw ) ) ⋅ ( 1 – Pcw4 ( T – Tini ) )
3
Pcgo = Pg – Po = ( Pcg1 + Pcg2 Sg + Pcg3 Sg ) ⋅ ( 1 – Pcg4 ( T – Tini ) )
where
Sw = water saturation
Sg = gas saturation
T = temperature
Definitions:
TENDPT isat
TEMP t1 t2 t3 ...
SWIR swir1 swir2 swir3 ...
SWC swc1 swc2 swc3 ...
SWRO swro1 swro2 swro3 ...
SORW sorw1 sorw2 sorw3 ...
SGR sgr1 sgr2 sgr3 ...
SGRO sgro1 sgro2 sgro3 ...
SORG sorg1 sorg2 sorg3 ...
KRWRO krwro1 krwro2 krwro3 ...
KROCW krocw1 krocw2 krocw3 ...
KRGRO krgro1 krgro2 krgro3 ...
Definitions:
imbibition
krw
drainage
Swr Swro
drainage
krow
imbibition
Swr Swro
Figure 4-6: Hysteretic Oil Relative Permeability Bounding Curves for a Water-
wet System.
drainage
krw
imbibition
Swr Swro
Figure 4-7: Hysteretic Water Relative Permeability Bounding Curves for an Oil-
wet System.
imbibition
krow
drainage
Swr Swro
Figure 4-8: Hysteretic Oil Relative Permeability Boundary Curves for an Oil-
wet System.
For a water-wet system, the ISAT array (Section 5.13.1) is used to assign
the drainage curves, and the ISATI array (Section 5.13.2) is used to assign
the imbibition curves. The definitions of ISAT and ISATI are reversed for
an oil-wet system. In other words, the ISAT array defines the bounding
curves for decreasing water saturation and the ISATI array defines the
bounding curves for increasing water saturation. Table ISATI is used for
krw and krow for water saturations less than Swr or greater than Swro
and for krg, krog, and all capillary pressures.
For Sw increasing:
Sw – Swr ewi
1 – --------------------------------
Swro – Swr
krw = krwi – Kwrev ( krwi – krwd ) --------------------------------------------
Swrev – Swr
1 – -----------------------------------
Swro – Swr
Sw – Swr eoi
1 – --------------------------------
Swro – Swr
krow = krowi + Korev ( krowd – krowi ) --------------------------------------------
Swrev – Swr
1 – -----------------------------------
Swro – Swr
For Sw decreasing:
Sw – Swr ewd
--------------------------------
Swro – Swr
krw = krwd + ( 1 – Kwrev ) ( krwi – krwd ) -----------------------------------
Swrev – Swr
-----------------------------------
Swro – Swr
Sw – Swr eod
--------------------------------
Swro – Swr
krowd – ( 1 – Korev ) ( krowd – krowi ) -----------------------------------
Swrev – Swr
-----------------------------------
Swro – Swr
where
krwi – krw
kwrev = ---------------------------------- at Swrev
krwi – krwd
krow – krowi
korev = ------------------------------------------ at Swrev
krowd – krowi
By default, Kwrev and Korev are initialized to 1.0 (to drainage curves for
water wet, to imbibition curves for oil-wet). If water saturation is at or
below Swr initially or if water saturation decreases to Swr or below during
simulation, kwrev and korev are reset to 0.0 in order to follow the
imbibition curves (water-wet) or drainage curves (oil-wet) when water
saturation increases. Similar logic applies to the endpoint at Swro. Kwrev
and Korev can be specified with the KWHYS and KOHYS arrays (Sections
5.45.1 and 5.45.2) to initialize to either bounding curve or a scanning
curve.
Table 4-1 and 4-2 summarize the internally coded properties for all
components that can be automatically determined. Using one of the
component identification codes contained in these tables as data in the
"COMPONENTS" input data will cause the tabulated values to be loaded
as data. Unless these values are overridden by the user in the
"PROPERTIES" data section, they will be used by default. 00
PR
PRORIG
EOS RK (ipvt)
SRK
ZJRK
00 Definitions:
00 Example:
00 C*********************************************
EOS PR 1
COMPONENTS
CO2 N2 C1 C2 C3 IC4 NC4 IC5 NC5 C6 C7 C8 C9 C10
C11 C14P C20P C27P C36P
C
The alphanumeric labels that will be used to identify the components are
defined here. If any of the component names contained in either Table 4-1
or Table 4-2 are used, the properties data will be picked from the tables,
subject to user override. If unrecognized labels are used, the user must
directly specify all of the properties data. 00
For multiple equation of state tables, the component names need only be
entered for the first table. Component names may also be entered for
subsequent tables, but the names must be the same as those entered for
the first table. 00
COMPONENTS
cmp1 cmp2 . . . cmpnc
00 Definition:
00 Example:
00 C*********************************************
EOS PR
COMPONENTS
CO2 N2 C1 C2 C3 IC4 NC4 IC5 NC5 C6 C7 C8 C9 C10
C11P C14P C20P C27P C36P
The properties table (Card 2 and Cards 3) may be truncated at any point
after molecular weight. All unread data is defaulted. The column titles on
the table title card (Card 2) must appear in the order shown. The number
of entries on Card 2 and Cards 3 must be the same. 00
F PSIA
(1) PROPERTIES C KPA
K PSIG
R K ⁄ C2
00 Definitions:
00 (1) Alpha label indicating that the units of the critical temperature values
(TC) are:
Alpha label indicating that the units of the critical pressure values (PC)
are:
PSIG Psig.
K/C2 Kg/cm2.
00 (2) The titles on this card must appear in the order shown. PCHOR is not
allowed in VIP-THERM.
00 Example:
00 EOS PR
C************ COMPONENT PROPERTIES **************
C
PROPERTIES
COMP MW TC PC ZC ACENTRIC OMEGAA OMEGAB VSHFT
DJK cards are used to define binary interaction coefficients djk. These are
used in the mixing rules that determine for mixtures the "A" parameter of
the equation of state. 00
Table 4-3 summarizes values of djk that have been internally coded for
various combinations of components. These values will be automatically
loaded as data when the corresponding component names are entered in
the COMPONENTS data. All binary combinations that involve an
unrecognized component will be assigned default values of zero. 00
DJK cmpj
cmpk djk
00 Definitions:
NOTE: 1. Repeat the data card containing cmpk for each component k
that interacts with component j.
00 Example:
00 EOS PR
00 COMPONENTS
C1N2 NO2C3 C4C5 C5C8 C9C19 C2035 C36P
C
C****************** COMPONENTS PROPERTIES *********************
C
PROPERTIES
COMP MW TC PC ZC ACENTRIC OMEGAA OMEGAB
00 C9C19 -0.7258570E-01
C2035 -0.9795381E-01
ENDEOS
C
00 By default, the coefficients C1, C2, C3, C4, and C5 are equal to 0.1023,
0.023364, 0.058533,-0.040758, and 0.0093324, respectively. Use of the LBC
keyword allows the user to modify these default coefficients to obtain a
better match with experimental viscosities.
00 The LBC keyword is not required if the user wishes to retain the default
coefficients.
00 The LBC keyword can appear anywhere after the PROPERTIES keyword
and data and before the ENDEOS keyword. The modified coefficients will
only apply for the current EOS table being defined.
LBC
C1 C2 C3 C4 C5
H i* = h i0
* + h * T + h * T2 + h * T3 + h * T4 + h * T5
i1 i2 i3 i4 i5
For components not listed in Table 4-4, default ideal gas state enthalpy
coefficients are calculated using the Kesler and Lee correlation (References
35, 36, and 37). This correlation is not always applicable. It fails if the
Any of the default values for ideal gas state enthalpy coefficients may be
overridden by means of HSTAR data.
HSTAR * * * * * *
cmpi h io h i1 h i2 h i3 h i4 h i5
Definitions:
NOTE: Water ideal gas state enthalpy coefficients are internally hardcoded
and may not be input.
PR
RK
EOSSEP
SRK
ZJRK
00 Definitions:
The equation of state to be used for PVT properties for surface calculations
is specified on the EOSSEP data card. The default equation of state is the
reservoir equation of state. 00
00 The additional data required for three parameter EOS treatment is input
by supplying VSHFT data as a part of the PROPERTIES data.
00 Examples:
00 EOS PR
COMPONENTS
C02 N2 C1 C2 C3 IC4 NC4 IC5 NC5 C6 C7 C8 C9 C10
C11P C14P C20P C27P C36P
C
PROPERTIES
COMP MW TC PC ZC ACENTRIC OMEGAA OMEGAB
.
.
.
DJK
.
.
.
EOSSEP SRK
PROPERTIES
COMP MW TC PC ZC ACENTRIC OMEGAA OMEGAB
.
.
.
STKZDN
.
DJKSEP
.
.
.
ENDEOS
STKZDN
den1 . . . dennc
00 Definition:
NOTE: This data can be entered only after the EOSSEP card. If the equation
of state data for surface separators is identical to that for the
reservoir, the EOSSEP card can be followed by STKZDN card.
00 Examples:
00 C
EOSSEP SRK
STKZDN
0.632 0.641 0.658 0.70 0.80 0.777
ENDEOS
DJKSEP cmpj
cmpk djk
00 Definitions:
DJKSEP cards are used to define binary interaction coefficients djk for
surface separators. These are used in the mixing rules that determine for
mixtures the ’A’ parameter of the equation of state. 00
Table 4-3 summarizes values of djk that have been internally coded for
various combinations of components. These values will be automatically
loaded as data when the corresponding component names are entered in
the COMPONENTS data. All binary combinations that involve an
unrecognized component will be assigned default values of zero. 00
The data entered using DJKSEP card is used for all stages. The data can be
selectively overridden by using the DJKSEP card as a part of the separator
stage data. 00
Repeat the data card containing cmpk for each component k that interacts
with component j. 00
00 Examples:
00 EOS PR
COMPONENTS
C1 N2 C02C3 C4C5 C6C8 C9C19 C2035 C36P
C
C COMPONENTS PROPERTIES
C
PROPERTIES
COMP MW TC PC ZC ACENTRIC OMEGAA OMEGAB
.
.
.
C
DJK C02C3
C1 N2 -0.508225E-01
DJK C4C5
C1 N2 0.4172094E-01
C02C3 -0.2031559E+00
DJK C6C8
C1 N2 -0.3255523E-01
C02C3 -0.6801761E-01
C4C5 -0.1434877E+00
DJK C9C19
C1 N2 0.1238792E+00
C02C3 -0.8605607E-03
C4C5 0.2762086E+00
C6C8 -0.6246630E-01
DJK C2035
C1 N2 0.8639371E-01
2C3 0.1466657E-01
C4C5T -0.2710768E+00
C6C8 -0.1425948E+00
C9C19 0.8496494E-02
DJK C36P
C1 N2 0.8916950E-01
C02C3 0.1177855E+00
C4C5 -0.2947721E-01
C6C8 0.5526345E-01
C9C19 -0.7258570E-01
C2035 -0.9795381E-01
C
EOSSEP SRK
C
00 PROPERTIES
COMP MW TC PC ZC ACENTRIC OMEGAA OMEGAB
C
DJKSEP C02C3
C1 N2 -0.508225E-01
DJKSEP C4C5
Cl N2 0.4172094E-01
C02C3 -0.2031559E+00
DJKSEP C2035
.
. .
.
C1 N2 0.8639371E-01
C02C3 0.1466657E-01
C4C5 -0.2710768E+00
C6C8 -0.1425948E+00
C9C19 0.8496494E-02
ENDEOS
ENDEOS
00
MOL. ACENTRIC
COMPONENT TC, ˚F. PC, psia ZC
WT. FACTOR
* Adjusted to fit the vapor pressure curve using the P-R equation.
00
00
00
Described in this section are oil and gas phase viscosity data and data for
options in which EOS calculations (oil and gas phase fugacity coefficients,
densities, and enthalpies) may be selectively replaced by tables or
correlations.
Four methods are available for modeling oil phase viscosity - (A) the
Pederson Model28, (B) mixing rule option, (C) the molecular weight-
dependent intercept logarithmic correlation option, and (D) the Lohrenz,
Bray, and Clark model38 (not recommended for non-isothermal problems).
VISOIL PEDERSON
VISK
i ki
. . .
VISKJ
i j Xij
. . . . . .
Definitions:
ki Overwriting ki value.
i Component i.
j Component j.
k1 = 9.74602
k2 = 18.0834
k3 = 4126.66
k4 = 44.6055
k5 = 0.976544
k6 = 81.8134
k7 = 15649.9
Examples:
ENDEOS
VISOIL PEDERSON
VISK
1 42.5
2 0.0002279
3 11770
4 600.3
5 16.49
6 1614
7 0.05552.
VISKJ
35 -0.36409
36 -0.38517
37 -0.8
38 -0.62763
OILMF 1
0.0001 0.2917 0.0045 0.0024 0.1538 0.3003 0.1095 0.1377
GASMF 1
X X X X X X X X
µ 01
x1
* µ 02
x2
* µ 03
x3
. . . µ oxNC
NC
is used to obtain the oil phase viscosity from the component viscosities µoi
which are interpolated/extrapolated from the tabular data using the
equation
VOMIN vomin
Definitions:
tj Temperature, ˚F (˚C).
min
00 where µo is oil viscosity in centipoise, µ o is minimum oil viscosity, T is
temperature in degrees Rankin, B is a constant, and A is a function of oil
min
molecular weight. Note that the use of µ o = .200001 is equivalent to the
standard ASTM viscosity relationship in which viscosity approaches a
VISOIL
(VOMIN vomin)
B b
A MW
a1 mw1
a2 mw2
. .
. .
. .
Definitions:
VISOIL LBC
Three methods are available fo modeling gas phase viscosity - (A) the
Pederson Model28, (B) the component gas viscosity option, and (C) the
Lohrenz, Bray, and Clark model38.
In all three methods, steam tables are used to obtain the steam viscosity,
which is combined with the hydrocarbon gas viscosity (determined by
one of the methods) by mole fraction average. Steam viscosity is given by
µ s ( T, P ) if P < Pwsat
µs =
µ sat ( T ) if ≥ Pwsat
s
SAT
where: µ s ( T, P ) is the steam table value at T,P ; µ s ( T ) is the saturated
value at T, and Pwsat is the water saturation pressure as a function of
temperature.
VISGAS PEDERSON
VISK
i ki
. . .
VISKJ
i j Xij
. . . . . .
Definitions:
ki Overwriting ki value.
i Component i.
j Component j.
k1 = 9.74602
k2 = 18.0834
k3 = 4126.66
k4 = 44.6055
k5 = 0.976544
k6 = 81.8134
k7 = 15649.9
Examples:
ENDEOS
VISGAS PEDERSON
VISK
1 42.5
2 0.0002279
3 11770
4 600.3
5 16.49
6 1614
7 0.05552.
VISKJ
35 -0.36409
36 -0.38517
37 -0.8
38 -0.62763
OILMF 1
0.0001 0.2917 0.0045 0.0024 0.1538 0.3003 0.1095 0.1377
GASMF 1
X X X X X X X X
This option is identical to the component oil viscosity option except that
the component viscosities are linearly interpolated/extrapolated from the
tabular data and that a mole fraction average mixing rule is used:
NCV
∑
*
µg = µ gi y i
i=1
Definitions:
tj Temperature, ˚F(˚C)
VISGAS LBC
and are used to obtain the oil phase density in the volume-average mixing
rule
1
ρ o = ---------------------
NC
-
xi Mi
∑ -----------
ρ oi
-
i=1
DENO COMP
COIL coil1 coil2 . . . coilNC
CTEOIL cteoil1 cteoil2. . . cteoilNC
DOSTD dostd1 dostd2 . . . dostdNC
Definitions:
In this option, specific heat capacity coefficients for components in the oil
and gas phases and volatile component heats of vaporization are used to
obtain oil and gas phase enthalpies, thus eliminating the equation-of-state
enthalpy calculations.
Component specific heat capacities in the oil and gas phase are defined as
T
h oi = ∫ Cpoi dT
Ts
TB i T
h gi = ∫ C poi dT + ∫ Cpgi dT + ∆H vi
Ts TB i
The oil and gas phase enthalpies are then given by mole fraction averages
of the component values:
NC
ho = ∑ xi hoi
i=1
NCV
hg = ∑ y i h gi + y s h gs
i=1
where ys is the steam mole fraction in the gas phase and hgs is the steam
enthalpy, which is taken from the steam tables as:
H ( P, T ) + ∆H vw ( T )if P > Pwsat
h gs = WL
H WV ( P, T ) if P ≤ Pwsat
where: HWL (P,T) and HWV (P,T) are liquid water and steam enthalpies at
P,T; ∆HVW(T) is the saturated water heat of vaporization at T; and Pwsat is
the water saturation pressure as a function of temperature.
ENTHALPY COMP
CPO0 cpo01 cpo02 . ..cpo0NC
CPO1 cpo11 cpo12 . ..cpo1NC
CPG0 cpg01 cpg02 . ..cpg0NCV
CPG1 cpg11 cpg12 . ..cpg1NCV
CPG2 cpg21 cpg22 . ..cpg2NCV
CPG3 cpg31 cpg3s . . . cpg3ncv
HVAP hvap1 hvap2 . . .havapNCV
TB tb1 tb2 . . .tbNC
Definitions:
cpo0i, cpo1i Values of the 0th and 1st degree specific heat
capacity coefficients for oil component i, Btu/
(LB ˚Fdegree+1) (KJ/(KG ˚Cdegree+1)).
CPG0, CPG1, CPG2, Apha labels indicating that the next NCV
CPG3 entries on each card are the 0th or 1st or 2nd or
3rd degree specific heat capacity coefficients for
the gas phase components. CPG3 data is
optional.
cpg0i, cpg1i, cpg2i Values of the 0th and 1st and 2nd and 3rd
cpg3 degree specific heat capacity coefficients for gas
component i, Btu/(LB ˚Fdegree+1) (KJ/(KG
˚Cdegree+1)).
In the default treatment of water in the vapor phase (using the Equation of
State), the k-values are assumed to be three phase values which account
for the presence of water in the vapor phase.
B –Di
KV i = A i + -----i + C i *P *EXP --------------
P T – E i
KVCOR
COMP A B C D E
1 a1 b1 c1 d1 e1
2 a2 b2 c2 d2 e2
. . . . . .
. . . . . .
. . . . . .
NCV aNCV bNCV cNVC dNCV eNCV
KVTBS
P p1
COMP 1 2 3 ... NCV
TEMP
t11 kv111 kv211 kv311 ...
t21 kv121 kv221 kv321 ...
t31 kv131 kv231 kv331 ...
. . . .
. . . .
. . .
P p2
COMP 1 2 3 ...
TEMP
t21 kv112 kv212 kv312 ...
t22 kv122 kv222 kv322 ...
t32 kv132 kv232 kv332 ...
P p3
.
.
.
Definitions:
00
00 Definitions:
00 Definitions:
NOTE: 1. Both the OILMF and GASMF cards must be entered when
defining an equilibrium region composition that is constant.
However, the composition for the equilibrium phase may be
defaulted by using the alpha label X in place of the
composition values xi or yi. In this case, the primary phase
composition must be entered first, with the defaulting
equilibrium phase following.
00 Example:
00 In VIP-THERM, all mole fractions are input on a water-free basis. For gas
compositions, enter zero values for non-volatile components.
00 Definitions:
di Depth, ft (m).
00 Example:
00 EOS PR
COMPONENTS
CO1 N2 C1 C2 C3 IC4 NC4 IC5 NC5 C6 C7 C8 C9 C10
C11P C14P C20 C27P C36P
. . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . .
COMPOSITION 1
DEPTH PSAT Z
8393.0 4312.0 0.114694 0.00596171 0.0535574 0.0253688
0.00375082 0.00931554 0.00287086 0.00407180 0.0452978 0.00641069
0.00608870 0.00382681 0.00266587 0.00422379 0.00241988
0.000426779 0.0000480376 0.00000520974
8568.5 4312.0 0.092183 0.00286226 0.449443 0.0467961 0.0290827
0.00647076 0.0177852 0.00652344 0.00992164 0.0138857 0.0260174
0.0308275 0.026336 0.0222263 0.0454738 0.0655906 0.0668725
0.0214402 0.0219302
8784.0 4312.0 0.0725313 0.00335000 0.444713 0.0461563 0.0293813
0.00583125 0.0140063 0.00593125 0.00847500 0.0126438 0.0246996
0.0294541 0.0232041 0.0197861 0.0495547 0.0694539
COMPOSITION ieq
DEPTH PSAT X Y (T) Z
d1 ps1 x1 y1 (t1) z1,1 . . . z1,nc
. . . . .
. . . . .
dn psn xn yn (tn) zn,1 . . . zn,nc
00 Definitions:
di Depth, ft(m).
VISPE (ipvt)
VISK
i ki
. . .
VISKJ
i j Xij
. . . . . .
ENDVIS
00 Definitions:
ki Overwriting ki value.
i Component i.
j Component j.
00 k1 = 9.74602
k2 = 18.0834
k3 = 4126.66
k4 = 44.6055
k5 = 0.976544
k6 = 81.8134
k7 = 15649.9
00 Examples:
00 OILMF 1
0.0001 0.2917 0.0045 0.0024 0.1538 0.3003 0.1095 0.1377
GASMF 1
X X X X X X X X
VISPE
VISK
1 42.5
2 0.0002279
3 11770
4 600.3
5 16.49
6 1614
7 0.05552.
VISKJ
35 -0.36409
36 -0.38517
37 -0.8
38 -0.62763
ENDVIS
GASPLANT
NKEY ikey ibat
KEYCMP
vkcmp1 (vkcmp2... vkcmpj ... vkcmpNI)
(enter number of KEY component plus composition values for
interpolation, j = 1 to number of interpolation points, NI)
PLNTRY
pr1,1 (pr1,2 ... pr1,j ... pr1,NI)
pr2,1 (pr2,2 ... pr2,j ... pr2,NI)
. . . . . .
. . . . . .
pri,1 (pri,2 ... pri,j ... pri,NI)
. . . . . .
prNC,1 (prNC,2 ... prNC,j ... prNC,NI)
(enter the plant liquid molar recovery fractions for each
interpolation point, j = 1 to number of interpolation values, NI
and repeat for all components, i = 1 to the number of components,
NC)
00 Definitions:
00 Notes:
■ Input to a gas plant is the total well stream, while output is determined
by the molar liquid recovery fractions. There are no surface flash
calculations as are carried out with a normal surface separator.
■ The user may optionally enter the surface separator equation of state
parameters. These parameters will be used for the stock tank density
calculations to obtain the surface rates. The new equation of state
parameters must follow the last stage of the battery to which they
apply. The values of the separator equation of state parameters will
default to the reservoir if not given.
00 Example:
C=================================
C GAS PLANT SURFACE SEPARATOR
C=================================
GASPLANT
NKEY 6 1
KEYCMP
C DEFINE KEY COMPONENT PLUS FRACTIONS (NC = 6 TO 8)
C NUMBER OF INTERPOLATION POINTS (NI = 11)
.9999 .108 .104 .098 .075 .065 .047 .028 .018 .010 .000
C DEFINE COMPONENT LIQUID RECOVERIES (NI = 11, NC = 8)
.0240 .0240 .0240 .0220 .0170 .0140 .0100 .0050 .0030 .0020 .0020
.0070 .0070 .0070 .0060 .0050 .0040 .0030 .0010 .0010 .0000 .0000
.0610 .0610 .0590 .0560 .0430 .0370 .0270 .0140 .0090 .0060 .0060
.1790 .1790 .1750 .1790 .1370 .1220 .0920 .0550 .0400 .0290 .0290
.4680 .4680 .4640 .4530 .4000 .3710 .3050 .2200 .1770 .1480 .1480
.9960 .9960 .9960 .9960 .9940 .9930 .9890 .9790 .9690 .9590 .9590
1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
NOTE: CO2 solubility tables must follow all mole fraction and surface
separation data.
00 Use of this option allows carbon dioxide to dissolve in water and provides
for the dependence of water volume factor and viscosity upon carbon
dioxide content.
00 The input data describes water CO2 content and water properties for
water at equilibrium with 100% CO2 vapor.
00 Definitions:
sat
cw = cw (rsw)
sat
pw = psw (rsw)
sat sat
Bw = 1 ⁄ bw
00 Notes:
■ The last data card (pswn = 0, rswn = 0) defines pure water properties.
■ The number of oil PVT tables read must be greater than or equal to the
maximum value in the IPVT array (Section 5.14).
■ Only one PVT table can apply to any one equilibrium region (see
Section 5.28). However, one PVT table can apply to any number of
different equilibrium regions.
00
Bo Vo
00 Two distinct types of black oil tables may be entered. One type of input
uses the BOTAB table. This approach is best when the results of a
differential liberation experiment are available. The user can input the
experimental results directly, and the simulator can account for variations
in produced gas gravity. Separator calculations are treated independently,
so that changes in separator conditions are allowed. Internally, data is
converted to K-value data and internal calculations are performed using
K-values and compositions.
00 Alternatively, black oil table data may be input with a combination of the
BOETAB, BOOTAB, BODTAB, BOGTAB, and BDGTAB tables. The
internal representation of these tables is more traditional, because black oil
00 The user may choose this alternate form of black oil table data, which we
refer to as the modified black oil approach, for a number of reasons.
Among them are:
3. The modified black oil tables are more flexible because they allow for
separate input of oil phase and gas phase tables. This can be
considerably more convenient, especially for computer-generated
black oil data. The modified tables are also considerably more flexible
for the input of undersaturated data, both in terms of the form of the
data and the location of the data points. They are also more
convenient in that they allow data to be monotonically increasing or
decreasing.
00 The original BOTAB table does not allow for the direct specification of a
heavy component in the gas phase, so it cannot be used for gas condensate
systems. The modified tables allow the user to specify solution oil-gas
ratios, which makes them suitable for modeling gas condensate systems.
5. Gas gravity relative to air for the gas removed at each stage of the
differential expansion. In most "black-oil" models this is treated as a
constant, but gas gravity actually changes with pressure and is
reported at each stage by the PVT laboratory. If the correct pressure
variation is unknown, a constant value can be entered at each
pressure.
Since the molecular weight of the residual oil is not usually reported in the
differential expansion test, the preferred method is to use values of oil
density and molecular weight at the original bubblepoint of the reservoir
oil. The molecular weight of the bubblepoint oil is easily calculated from
the compositional analysis of the reservoir fluid usually given in the
reservoir fluid study. Residual oil molecular weight can be approximated
but will lead to inaccuracies in produced oil molecular weight and molar
production rates. Surface volumes will not be affected. 00
The data for undersaturated conditions is normally available only for the
original oil sample. (Oil that is saturated at the measured bubblepoint of
the reservoir.) This is the only data required by the program, but in the
event that undersaturated data has been measured for samples at various
stages of depletion, the additional data can be used to better define
volumetric behavior in the undersaturated region. 00
the gas z-factor will also be derived. During the simulation, K-values and
z-factors are used for phase behavior and density calculations, rather than
Rs, Bo, and Bg values.
The PVT data for saturated conditions (Type 5 Data Cards) is subjected to
a consistency check that determines whether the data could lead to a
nonpositive total hydrocarbon compressibility. Total hydrocarbon
compressibility is given by the following equation:
1 dV h
00 c h = ------ ---------- ,
V h dp
dR s dB o
00 B g --------- – --------- > 0 .
dp dp
This simply states that; for an increase in pressure, the reduction in total
volume due to gas going into solution must be greater than the increase
due to oil swelling. If Bg is in RB/MSCF, dRs/dp is in (MSCF/STB)/psi
and dBo/dp is in (RB/STB)/psi, the units will be correct. Bg should be
evaluated at the lower pressure.
DOB WTRO
DOR WTRO
(2) API WTRO
DOS WTOS PSAT
APIS WTOS PSAT
dob wtro
dor wtro
(3) api wtro
dos wtos psat
apis wtos psat
ZG
(4) PSAT RS BO GR VO VG
BG
zg
(5) psat rs bo gr vo vg
bg
(At least two data cards are required here.)
00 Definitions (1):
00 (1) BOTAB Alpha label indicating the data being read are
differential expansion data.
00 Definitions (2):
00 Definitions (3):
00 Definitions (4):
00 Definitions (5):
00 Definitions (6):
00 Definitions (7):
00 Definitions (8):
00 Notes:
b. The last psat value must equal the standard pressure PS (Section
2.2.4).
e. If applicable, gas formation volume factor must increase (bgj > bgj-
1).
c. Relative oil volume factor must be less than 1.0 (bofacij < 1.).
5. The BOFAC and VOFAC alpha labels on the Type 7 Data Card are
repeated as a pair for each value of psati entered on the Type 6 Data
Card.
6. The bofaci and vofaci values on the Type 8 Data Card are repeated as a
pair for each value of psati entered on the Type 6 Data Card.
7. A non-zero value for solution gas-oil ratio (Type 5 Data Card) in the
last entry of the table (at standard pressure) is allowed but is not
physically possible. Gas cannot be liberated by cooling a saturated
fluid from reservoir temperature to standard temperature at constant
(standard) pressure. A non-zero value would indicate a composition
change between the conditions for the last entry in the table and
standard conditions. This implied composition change is not
considered.
8. If the initial saturation pressure is not input (psat, Type 3 Data Card),
the default value is equal to the saturation pressure at the gas-oil
contact in the IEQUIL table (Section 4.2.1.2).
00 Example:
00 C
BOTAB 1
API WTRO
27.40 200
PSAT RS BO ZG GR VO VG
DOB WTRO GR
(2) DOR WTRO GR
API WTRO GR
dob wtro gr
(3) dor wtro gr
api wtro gr
(4) PSAT RS BO VO
(5) psat rs bo vo
(At least two data cards are required here.)
00 Definitions (1):
00 Definitions (2):
00 Definitions (3):
00 Definitions (4):
00 Definitions (5):
00 Definitions (6):
00 Definitions (7):
00 Definitions (8):
4.5.3 Gas PVT Data for the GASWATER Option (BGTAB) (VIP-ENCORE)
zg
(5) psat vg
bg
(At least two data cards are required here.)
00 Definitions (1):
BGTAB Alpha label indicating the data being read are gas
properties.
00 Definitions (2):
00 Definition (3):
00 Definitions (4):
00 Definitions (5):
4.5.4.1 BOETAB
(1)BOETAB ipvt
(4)PSAT RS RV BO BG VO VG
(5)psat rs rv bo bg vo vg
(At least two Type 5 data cards are required.)
psat
(6) PSAT
psat 1
2 . . . psat i
RS rs 1 rs rs
2 i
(Only one Type 6 data card is required. Cards 7 and 8 must be
repeated for each Card 6.)
VO . . VO
P BO VO . BO
(7)
BO
DP BOFAC VOFAC BOFAC VOFAC BOFAC VOFAC
p bo 1 vo 1 bo 2 vo 2 bo i vo i
(8) . . .
dp bofac vofac bofac vofac
bofac 1 vofac 1 2 2 i i
(Only one Type 8 data card is required for each Type 7 data
card.)
p1 p
(9) P 2 . . . pi
RVSAT rvsat 1 rvsat rvsat
2 i
BG VG BG VG BG VG
(10) RV . . .
BGFAC VGFAC BGFAC VGFAC BGFAC VGFAC
bg 2 vg 2 bg i vg i
. . .
bg 1 vg 1
(11) rv
bgfac 1 vgfac 1 bgfac vgfac bgfac vgfac
2 2 i i
(Only one Type 11 data card is required for each Type 10 data
card.)
00 Definitions (1):
(1) BOETAB Alpha label indicating the data being read are extended
black oil data, with both oil and gas data read in the same
table.
00 Definitions (2):
DOB Alpha label indicating that the data value that appears in
this location on the following card is the density of the
residual oil at standard conditions as defined in Section
2.2.4.
API Alpha label indicating that the data value that appears in
this location on the following card is the API gravity of
the residual oil at standard conditions as defined in
Section 2.2.4.
WTRO Alpha label indicating that the data value that appears in
this location on the following card is the molecular
weight of the residual oil. This alpha label must be used
if the previous alpha label on this data line is either DOB,
DOR, or API.
DGB Alpha label indicating that the data value that appears in
this location on the following card is the density of the
gas at standard conditions as defined in Section 2.2.4.
Definitions (3): 00
Definitions (4): 00
00 Definitions (5):
data table. For dead oils, data card types 6, 7, and 8 are
not allowed.
NOTE: If all values of zero in the RS column, the oil phase does not contain any
solution gas. Therefore, input of Data Card Types 6, 7, and 8 is not
permitted.
NOTE: Data Card Types 6, 7, and 8 may be repeated as a group as many times as
are necessary.
Definitions (6): 00
00 Definitions (7):
00 Definitions (8):
NOTE: If all values of zero in the RV column, a dry gas is defined. Therefore,
input of Data Card Types 9, 10, and 11 cannot be permitted.
NOTE: Data Card Types 9, 10, and 11 may be repeated as a group as many times
as are necessary. However, one set is required.
00 Definitions (9):
Definitions (10): 00
00 Definitions (12):
4.5.4.2 BOOTAB
(1)BOOTAB ipvt
(4)RS PSAT BO VO
(5)rs psat bo vo
(At least two Type 5 data cards are required.)
psat
psat 1 psat
(6) PSAT 2
. . . i
RS rs 1 rs rs
2 i
(Only one Type 6 data card is required. Cards 7 and 8 must be
repeated for each Card 6.)
P VO
BO VO . BO VO
(7) . . BO
DP BOFAC VOFAC BOFAC VOFAC BOFAC VOFAC
p bo 1
vo 1bo 2 vo 2 bo i vo i
(8) . . .
dp bofac 1 vofac 1 bofac vofac bofac vofac
2 2 i i
(Only one Type 8 data card is required is for each Type 7
data card.)
00 Definitions (1):
(1) BOOTAB Alpha label indicating that live black oil data is being
read.
00 Definitions (2):
DOB Alpha label indicating that the data value that appears in
this location on the following card is the density of the
residual oil at standard conditions as defined in Section
2.2.4.
API Alpha label indicating that the data value that appears in
this location on the following card is the API gravity of
the residual oil at standard conditions as defined in
Section 2.2.4.
WTRO Alpha label indicating that the data value that appears in
this location on the following card is the molecular
weight of the residual oil. This alpha label must be used
Definitions (3): 00
Definitions (4): 00
00 Definitions (5):
NOTE: If all values are zero in the RS column, the oil phase does not contain any
solution gas. Therefore, input of Data Card Types 6, 7, and 8 is not
permitted.
NOTE: Data Card Types 6, 7, and 8 may be repeated as a group as many times as
are necessary.
Definitions (6): 00
00 Definitions (7):
00 Definitions (8):
4.5.4.3 BODTAB
00 (4) PSAT BO VO
00 (5) psat bo vo
(At least two Type 5 data cards are required.)
00 Definitions (1):
(1) BODTAB Alpha label indicating that dead black oil data is being
read.
00 Definitions (2):
DOB Alpha label indicating that the data value that appears in
this location on the following card is the density of the
residual oil at standard conditions as defined in Section
2.2.4.
API Alpha label indicating that the data value that appears in
this location on the following card is the API gravity of
the residual oil at standard conditions as defined in
Section 2.2.4.
WTRO Alpha label indicating that the data value that appears in
this location on the following card is the molecular
weight of the residual oil. This alpha label must be used
if the previous alpha label on this data line is either DOB,
DOR, or API.
Definitions (3): 00
Definitions (4): 00
00 Definitions (5):
4.5.4.4 BOGTAB
00 (2) DGB
GR
00 (3)
dgb
gr
00 (4) PSAT RV BG VG
00 (5) psat rv bg vg
p1 p p
(6) P 2 i
00
BG VG BG VG BG VG
00 (7) RV . . .
BGFAC VGFAC BGFAC VGFAC BGFAC VGFAC
bg 2 vg 2 bg i vg i
. . .
bg 1 vg 1
(8) rv
bgfac vgfac bgfac vgfac
00
bgfac 1 vgfac 1 2 2 i i
(Only one Type 8 data card is required for each Type 7 data
card.)
00 Definitions (1):
(1) BOGTAB Alpha label indicating that wet gas data is being read.
00 Definitions (2):
DGB Alpha label indicating that the data value that appears in
this location on the following card is the density of the
gas at standard conditions as defined in Section 2.2.4.
00 Definitions (3):
00 Definitions (4):
00 Definitions (5):
NOTE: If all values of zero in the RV column, a dry gas is defined. Therefore,
input of Data Card Types 6, 7, and 8 cannot be permitted.
NOTE: Data Card Types 6, 7, and 8 may be repeated as a group as many times as
are necessary. However, one set is required.
00 Definitions (6):
Definitions (7): 00
00 Definitions (8):
4.5.4.5 BDGTAB
00 (2) DGB
GR
00
dgb
00 (3)
gr
00 (4) PSAT BG VG
00 (5) psat bg vg
00 Definitions (1):
(1) BDGTAB Alpha label indicating that dry gas data is being read.
00 Definitions (2):
DGB Alpha label indicating that the data value that appears in
this location on the following card is the density of the
gas at standard conditions as defined in Section 2.2.4.
00 Definitions (3):
00 Definitions (4):
00 Definitions (5):
4.5.4.6 Examples
00 RV BG VG BG VG BG VG
00 0. 12.5694 0.0143 4.4010 0.0147 2.6341 0.0154
00 P 1780. 2280.
00 RV BG VG BG VG
00 0. 1.8760 0.0164 1.4631 0.0175
00 P 2780. 3280. 3780.
00 RV BG VG BG VG BG VG
00 0. 1.2086 0.0189 1.0388 0.0203 0.9191 0.0218
00 P 4280 4780
00 RV BG VG BG VG
00 0. 0.8309 0.0234 0.7636 0.0249
00
00 BOETAB 1
00 DOB WTRO DGB
00 42.5 200. 0.495
00 C
00 C SATURATED DATA
00 C
00 PSAT BO VO RS RV BG VG
00 280. 1.1258 0.1962 0.0 0.02207 12.5134 0.0144
00 780. 1.1132 0.1966 0.0 0.01806 4.3571 0.0149
00 1280. 1.1024 0.2051 0.0 0.01797 2.5938 0.0157
00 1780. 1.0927 0.2100 0.0 0.01963 1.8368 0.0168
00 2280. 1.0841 0.2200 0.0 0.02590 1.4251 0.0184
00 2780. 1.0764 0.2300 0.0 0.02590 1.1728 0.0202
00 3280. 1.0693 0.2452 0.0 0.02590 1.1728 0.0202
00 3780. 1.0630 0.2815 0.0 0.03406 0.8917 0.024
00 4280. 1.0570 0.3378 0.0 0.03779 0.8094 0.0270
00 4780. 1.0516 0.4435 0.0 0.04004 0.7486 0.0292
00 C
00 C UNDERSATURATED DATA
00 C
00 P 280. 780. 1280.
00 RV BG VG BG VG BG VG
00 0. 12.5694 0.0143 4.4010 0.0147 2.6341 0.0154
00 P 1780. 2280.
00 RV BG VG BG VG
00 0. 1.8760 0.0164 1.4631 0.0175
00 P 2780. 3280. 3780.
00 RV BG VG BG VG BG VG
00 0. 1.2086 0.0189 1.0388 0.0203 0.9191 0.0218
00 P 4280 4780
00 RV BG VG BG VG
00 0. 0.8309 0.0234 0.7636 0.0249
00
00
00 Example 3:
00 C Separate black live oil and live gas tables, with different pressure
C entries. RSM keyword invoked, so solution gas-oil ratio divided by
C 1000.
00 C
00 BOOTAB
00 DOB WTRO
00 48.25 200.
00 C
00 C SATURATED DATA
00 C
00 RS PSAT BO VO
00 0.00001 14.7 1.103 1.126
00 0.040 264.7 1.126 0.960
00 0.084 514.7 1.152 0.809
00 0.187 1014.7 1.203 0.703
00 0.292 1514.7 1.252 0.615
00 0.399 2014.7 1.303 0.535
00 0.511 2514.7 1.357 0.450
00 0.627 3014.7 1.415 0.387
00 0.757 3514.7 1.483 0.316
00 0.884 3900.0 1.539 0.270
00 1.400 5000.0 1.639 0.270
00 C
00 C UNDERSATURATED OIL DATA SET 1
00 C
00 RS 0.00001 0.040 0.084
00 P BO VO BO VO BO VO
00 5514.7 0.9688 1.127 0.9979 0.961 1.0302 0.810
00 C
00 C UNDERSATURATED OIL DATA SET 2
00 RS 0.187 0.292 0.399
00 P BO VO BO VO BO VO
00 5514.7 1.0932 0.704 1.1544 0.616 1.2176 0.536
00 C
00 C UNDERSATURATED OIL DATA SET 3
00 C
00 RS 0.511 0.627 0.757
00 P BO VO BO VO BO VO
00 5514.7 1.2838 0.451 1.354 0.388 1.4342 0.317
00 RS 0.884 1.40
00 P BO VO BO VO
00 5514.7 1.503 0.271 1.624 0.280
5. Oil phase mole fraction of the first component and/or the gas phase
mole fraction of the last component at each pressure step.
00 The K-value tabular data must be preceded by the alpha label KVALUES.
KVALUES
00 The alphanumeric labels that will be used to identify the components are
defined here.
COMPONENTS
cmp1 cmp2 . . . cmpnc
00 Definitions:
00 Example:
00 C************************
KVALUES
COMPONENTS
CO2 N2 C1 C2 C3 IC4 NC4 IC5 NC5 C6 C6 C8 C9 C10
C11P C14P C20P C27P C36P
PROPERTIES
COMP MW
cmp1 mw1
. .
. .
. .
cmpnc mwnc
00 Definitions:
NOTE: Repeat the data card containing cmpi and mwi for each
component.
For any value when the default is acceptable, enter the alpha label
X.
00 Example:
00 C KVALUES CASE
COMPONENTS
ONE TWO THREE
PROPERTIES
COMP MW
ONE mw1
TWO mw2
THREE mw3
ENDKV
ENDKV
X1
(6) PSAT ZG VG ZO VO YNC KV
X1 YNC
x1
(7) psat zg vg zo vo ync kv1 . . . kvNC
x1 ync
(At least two data cards are required here.)
00 Definitions (1):
00 Definitions (2):
00 Definitions (3):
00 Definitions (4):
00 Definitions (5):
00 Definitions (6):
00 Definitions (7):
00 Definitions (8):
00 Definitions (9):
00 Definition (10):
00 Definitions (11):
00 Definitions (12):
00 Definition (13):
00 Notes:
3. The following data restrictions apply to the Type 8 (Type 11) Data
Card:
4. The following data restrictions apply to the Type 10 (Type 13) Data
Cards:
c. Relative oil (gas) volume factor must be less than 1.0 (bofacij <1.).
e. Relative oil (gas) viscosity must be no less than 1.0 (vofacij ≥ 1.).
5. The BOFAC and VOFAC alpha labels on the Type 9 (Type 12) Data
Card are repeated as a pair for each value of psati entered on the Type
8 (Type 11) Data Card.
6. The bofaci and vofaci values on the Type 10 (Type 13) Data Card are
repeated as a pair for each value of psati entered on the Type 8 (Type
11) Data Card.
7. The optional residual oil Z-factor (Type 5 Data Card) is used to define
the oil Z-factor for the default separator. If not input, the default
separator oil Z-factor is set equal to that in the last Type 7 Data Card
(which should be at standard pressure) corrected for temperature.
Since this approximation neglects oil density variation between
reservoir and standard temperature, errors in surface volumes will
result if the default separator is used for separation and a Type 5 Data
Card is not input.
00 KVALUES
COMPONENTS
ONE TWO THREE
PROPERTIES
COMP MW
ONE MW1
TWO MW2
THREE MW3
ENDKV
KVTAB 1
PSAT ZG VG ZO VO X1 KV
psat zg vg zo vo x1 kv1 kv2 kv3
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
PSAT psat1
DP BOFAC VOFAC
dp bofac vofac
. . . .
. . .
. . .
GASMF 1 1
y1 y2 y3
OILMF 1 1
x1 x2 x3
In the Dead Oil PVT Property Tables option, oil density, enthalpy, and
viscosity are input in tabular form as functions of temperature and
pressure. In the Dead Oil PVT Properties Correlations option, oil density is
represented by compressibility and coefficient of thermal expansion, oil
enthalpy is represented by a constant heat capacity, and oil viscosity is
input as a function of temperature.
OILTABLES
MWOIL mwoil
DOSTD dostd
(RSF GR )
(rsf gr )
DOTAB
NP np
T t1 t2 t3 . . .
P
p1 d11 d12 d13 . . .
p2 d21 d22 d23 . . .
p3 d31 d32 d33 . . .
. . . .
. . . .
. . . .
HOTAB
NP np
T t1 t2 t3 . . .
P
p1 h11 h12 h13 . . .
p2 h21 h22 h23 . . .
p3 h31 h32 h33 . . .
. . . .
. . . .
. . . .
VOTAB
(VOMIN vomin)
NP np
T t1 t2 t3 . . .
P
p1 v11 v12 v13 . . .
p2 v21 v22 v23 . . .
p3 v31 v32 v33 . . .
. . . .
. . . .
. . . .
(TTAB (ieqr))
(DEPTH T)
(depth1t1)
(depth2t2)
Definitions:
OILTABLES Alpha label indicating that the data being read are
Dead Oil PVT Property Tables.
NOTE: 1. The number and values of the temperature and pressure entries
in the DOTAB, HOTAB, and VOTAB tables must be the same.
Required data for the oil include molecular weight, density at standard
conditions, compressibility, coefficient of thermal expansion, heat capacity
and values of viscosity as a function of temperature.
δDO 1
CTEOIL = ------------- ----------------------
δP DOSTD
δDO 1
CTEOIL = – ------------- ----------------------
δT DOSTD
HO = CPOIL * (T - TS)
OILPROPS
MWOIL mwoil
DOSTD dostd
COIL coil
CTEOIL cteoil
CPOIL cpoil
VISOIL
(VOMIN vomin)
T VISCOSITY
t1 v1
t2 v2
. .
. .
. .
(RSF GR)
(rsf gr)
(TTAB ieqr)
(DEPTH T)
(depth1 t1)
(depth2 t2)
(. .)
(. .)
(. .)
Definitions:
OILPROPS Alpha label indicating that the data being read are
Dead Oil PVT Correlations data.
Both the default and the input separators can be accessed for surface
volume calculations by means of the REGSEP Card (Section 9.2). 00
SEPARATOR ibat
(PVTTABLE ipvt)
(1) STAGE TEMP PRES VFRAC VDEST LFRAC LDEST
(2) n tn pn vfn1 vdn1 lfn1 ldn1
(3) X X X vfn2 vdn2 lfn2 ldn2
(Type 2 and 3 Data Cards are repeated as
necessary to define all the stages of separation.)
00 Definitions:
m ( stage number )
vd n1 = GAS ( alpha label )
VENT ( alpha label )
lfn1 Fraction of the liquid stream leaving separator stage
n to be sent to the destination indicated by ldn1.
Values must lie in the range 0-1. If lfn1 is less than
1.0, then one Type 3 Data Card must be provided for
this stage so that lfn1 and lfn2 sum to 1.0.
m ( stage number )
ld n1 =
OIL ( alpha label )
00 Example:
00 C ========================================
C COMPOSITIONAL SEPARATOR BATTERY
C ========================================
SEPARATOR 1
STAGE TEMP PRES VFRAC VDEST LFRAC LDEST
1 100.0 665.0 1. GAS 1. 2
2 132.0 100.0 1. GAS 1. 3
3 126.0 40.0 1. GAS 1. 4
4 149.0 15.0 1. GAS 1. 4
5 60.0 14.7 1. GAS 1. OIL
OMEGAS istg
Ωa1 . . . Ωanc
OMEGBS istg
Ωb1 . . . Ωbnc
00 Definitions:
00 Definitions:
VSHFTS istg
vshfts1... vshftsnc
00 Definitions:
When “black-oil” type problems are run in VIP-ENCORE, the PVT data
are converted into a multicomponent format, including the use of K-
values and z-factors to calculate the phase behavior and volumetric
behavior of both oil and gas. This treatment makes it possible for
separator conditions to be exactly modeled, while using differential
expansion data to describe fluid behavior in the reservoir. This eliminates
the conflict between differential and flash volumetrics that creates
difficulty for conventional black-oil simulators. 00
The types and format of surface separation data depend upon the type of
PVT Property Data which is input (black-oil laboratory data or K-value
tabular data, see Section 4.4). 00
Both the default and the input separators can be accessed for surface
volume calculations by means of the REGSEP Card (Section 9.2). 00
For black-oil laboratory PVT data input, there are four options for the
treatment of separators. 00
00 These options are not mutually exclusive, i.e., any number of input
separator data options may be used in the same run. Order is not
important.
Use of the default separator implies that the differential liberation data
have been adjusted for flash conditions in the conventional manner (see
Section 2.2.2.2 (VIP-ENCORE)). If no adjustments are made and the
default separator is used, significant errors in surface volumes will result. 00
K-values are calculated from the differential liberation data using the
calculated composition of the total gas liberated from the initial saturation
pressure to standard pressure and the composition of the liquid at
standard pressure. A liquid compressibility factor at standard conditions
is calculated for each battery using the residual oil density and molecular
weight input in or derived from the differential liberation data. 00
If a non-zero value for gas-oil ratio at standard pressure was input in the
PVT data, use of the default separator will result in a match between the
separation gas-oil ratio at the initial saturation pressure and the “effective”
input gas-oil ratio at the initial saturation pressure. The “effective” gas-oil
ratio is equal to the difference between the input gas-oil ratio and the gas-
oil ratio at standard pressure. See Note 7 at end of Section 4.5.2.1. 00
Input data format for this option is similar to that given in the reference
manual for the simulation modules. In earlier documentation it has been
suggested that one use the K-values derived for the last step in the
differential liberation experiment. This may lead to large errors in surface
volumes. This option should be used only when accurate separator K-
values and liquid density are known. When using this option, one must
take great care to insure that the input K-values are consistent with the
internally defined component molecular weights, which are printed out in
the default separator properties table. 00
2. The density of the oil product (stock tank oil) leaving the last separator
stage.
3. The molecular weight of the stock tank oil or the saturation pressure of
the saturated oil feed.
4. The equilibrium K-values for each component in the system for each
stage of separation.
SEPARATOR ibat
PVTTABLE ipvt
(1) STAGE VFRAC VDEST LFRAC LDEST
(2) n vfn1 vdn1 lfn1 ldn1
(3) X vfn2 vdn2 lfn2 ldn2
(Data Cards are repeated as necessary to define all of the
stages of separation.)
PSATF
(4) DLIQ
MWL
psatf
(5) dliq
mwl
(6) KVALUES
(7) COMP STAGE 1 (STAGE 2 . . . STAGE n)
(8) comp kval1 (kval2 . . . kvaln)
(Enter one of these cards for each component in the fluid
system.)
00 Definitions:
m ( stage number )
vd n1 = GAS ( alpha label )
VENT ( alpha label )
lfn1 Fraction of the liquid stream leaving separator stage
n to be sent to the destination indicated by ldn1.
Values must lie in the range 0-1. If lfn1 is less than
1.0, then one Type 3 Data Card must be provided for
this stage, so that lfn1 and lfn2 sum to 1.0.
m ( stage number )
ld n1 =
OIL ( alpha label )
Each separator battery may contain any number of stages. Each stage
contains one feed stream and two output streams: vapor and liquid. Each
of the two output streams can itself be split into two streams, each of
which may be fed to any downstream separator stage, the gas sales line or
the oil sales line. 00
00 Example:
00 C==============================================
C BLACK-OIL AND K-VALUE SINGLE STAGE SEPARATOR
C==============================================
SEPARATOR 1
STAGE VFRAC VDEST LFRAC LDEST
1 1. GAS 1. OIL
DLIQ MWL
0.8966 200
00 C
KVALUES
COMP STAGE 1
1 89.4
2 0.0056
SEPTEST ibat
PVTTABLE ntab (itype)
PSATF BOF
psatf bof
P T GOR BOSTG GR
p1 t1 gor1 bostg1 gr1
. . . . .
. . . . .
. . . . .
pns tns gorns bostgns grns
(enter one data card for each stage)
00 Definitions:
ntab PVT table number for the fluid to which these data
apply.
This option is provided for cases in which laboratory separator test data is
not available or applicable, but a flash oil formation volume factor and
gas-oil ratio are known or estimated. These values must correspond to
feed of saturated oil at some pressure in one of the BOTAB tables. A stock
tank oil compressibility factor is calculated using the residual oil density
from PVT table ntab and a stock tank oil molecular weight computed
internally by material balance.
If the stock tank oil density or the stock tank oil molecular weight or the
surface gas gravity is known or can be estimated, then it is preferable to
use the Separator Test Data input option (with a single stage) as given in
the previous subsection. The required surface gas gravity can be
calculated from stock tank oil density or molecular weight by using the
equations:
S
5.6146 ⋅ 62.428 ⋅ RT ( BOF ⋅ DES – DENO )
GR = --------------------------------------------------------------------------------------------------------
-
GOR ⋅ 29 ⋅ P
S
00 and
00 where
BOSEP ibat
PVTTABLE ntab
PSATF BOF GOR
psatf bof gor
00 Definitions:
ntab PVT table number for the fluid to which these data
apply.
NOTE: This option may not be used in conjunction with the PVT
interpolation option.
A default separator is created for each input PVT property table. Each
default separator is single stage at standard conditions with a vapor
fraction of 1.0 sent to the destination GAS (gas sales line) and a liquid
fraction of 1.0 sent to the destination OIL (oil sales line) (see separator
configuration data below). 00
Default separator K-values are equal to those in the last entry of the
KVTAB table (at standard pressure). If a residual oil Z-factor is input in the
KVTAB table (Section 4.6), this is the default separator value. If not, the
default separator oil Z-factor is given by the oil z-factor in the last entry of
the table (at standard pressure) corrected for temperature. Since this
neglects oil density variation between reservoir and standard
temperature, some error in oil surface volume will result. 00
Separator properties may be input using the same format as given in the
reference manual for the simulation modules. The required data include: 00
2. The density of the oil product (stock tank oil) leaving the last separator
stage.
4. The equilibrium K-values for each component in the system for each
stage of separation.
SEPARATOR ibat
(1) STAGE VFRAC VDEST LFRAC LDEST
(2) n vfn1 vdn1 lfn1 ldn1
(3) X vfn2 vdn2 lfn2 ldn2
(Type 2 and 3 Data Cards are repeated as necessary to
define all of the stages of separation.)
(6) KVALUES
(7) COMP STAGE 1 (STAGE 2 . . . STAGE n)
(8) comp kval1 (kval2 . . . kvaln)
(Enter one of these cards for each component in the fluid
system.)
00 Definitions:
m ( stage number )
vd n1 = GAS ( alpha label )
VENT ( alpha label )
lfn1 Fraction of the liquid stream leaving separator stage
n to be sent to the destination indicated by ldn1.
Values must lie in the range 0-1. If lfn1 is less than
1.0, then one Type 3 Data Card must be provided for
this stage, so that lfn1 and lfn2 sum to 1.0.
m ( stage number )
ld n1 =
OIL ( alpha label )
00 Example:
00 C=============================================
C BLACK-OIL AND K-VALUE SINGLE STAGE SEPARATOR
C=============================================
SEPARATOR 1
STAGE VFRAC VDEST LFRAC LDEST
1 1. GAS 1. OIL
DLIQ MWL
0.8966 200
C
KVALUES
COMP STAGE 1
1 89.4
2 0.0056
00 Definitions:
00 Definitions:
The compaction tables provide for user input of pore volume and
permeability changes based on fluid pressure history. When the
compaction option is invoked with the COMPACT card in utility data, and
compaction tables are provided as shown in this section, then the pore
volumes and transmissibilities used by the simulator are multiplied by
reduction factors prior to use. The factors are looked up in the compaction
tables using either the current gridblock pressure (reversible option) or the
minimum historical gridblock pressure (irreversible option; default) in
each block. If more than one compaction table is provided, the ICMT array
(Section 5.19) is used to determine which one applies to each gridblock. 00
CMT ncmt
P PVMULT TAMULT (TVMULT)
p pvmult tamult (tvmult)
(At least two data cards are required.)
00 Definitions:
The values on the data cards following this title card appear in the order
defined by this title card. The alpha labels must appear in the order
shown. 00
00 Definitions:
The alpha labels on title Card 3 must appear in the order shown. 00
The values on Card 4 must appear in the order shown. There must be a
minimum of two Card 4 data cards 00
00 Example:
00 C
C
WIRCT 1
SWINIT 0.06
DSW PVMULT TAMULT TVMULT
.01 1.0 1.0 1.0
.05 0.95 0.95 0.95
.1 0.9 0.9 0.9
1. 0.8 0.8 0.8
00 SWINIT 0.12
DSW PVMULT TAMULT TVMULT
.01 1.0 1.0 1.0
.05 0.95 0.95 0.95
.1 0.9 0.9 0.9
00 SWINIT 0.2
DSW PVMULT TAMULT TVMULT
.01 1.0 1.0 1.0
.05 0.95 0.95 0.95
1. 0.8 0.8 0.8
00 SWINIT 0.3
DSW PVMULT TAMULT TVMULT
.01 1.0 1.0 1.0
1. 0.8 0.8 0.8
00 WIRCT 2
SWINIT 0.06
DSW PVMULT TAMULT TVMULT
.01 1.0 1.0 1.0
.05 0.95 0.95 0.95
1. 0.8 0.8 0.8
00 SWINIT 0.3
DSW PVMULT TAMULT TVMULT
.01 1.0 1.0 1.0
.1 0.9 0.9 0.9
1. 0.8 0.8 0.8
00
card. In this case, the light oil component is assigned as the first
component, the heavy oil component (dead oil) as the second component
and the injected solvent as the third component. In addition to the typical
black-oil table BOTAB that specifies the PVT properties and VLE for the
light and the heavy oil components, the user is required to input a solvent
property table, SLVTAB, immediately following the black oil table. 00
Note that in black-oil mode, with the multiple PVT table option, the user is
required to enter consecutively all BOTAB tables, followed by the same
number of SLVTAB and MISPTB tables (in pairs), i.e. for each PVT region,
a SLVTAB table must be immediately followed by an MISPTB table. After
all SLVTAB and MISPTB tables are entered, an MWS card must be input. 00
For the miscible model, there are two options for separators; default and
separator K-value input. Detailed descriptions of the options and the
input requirements are given in this section. The other two options
(separator test data input and black-oil separator data) for the
conventional black-oil model are not applicable. For the default option, a
large K-value (10E6) is used for the solvent component, if short solvent
property tables (P, BS and VS only in SLVTAB tables) are entered. The
component names for the three hydrocarbon components to be entered in
the separator K-value input option are “light,” “heavy,” and “solvent,”
which correspond to Components 1, 2, and 3, respectively. 00
The data cards that specify the injected solvent composition should take
into account the position of the solvent component. In the simulation
modules, pure solvent injection in the three and four-component cases
would be input as: 00
00 C 3 component
YINJ wl
1. 0.
and
C 4 component
YINJ wl
0.0 0.0 1.0
respectively, where wl is the solvent injection well list. Note that the
solvent is identified as the first hydrocarbon component in the three-
component case and the third hydrocarbon component in the four-
component case. 00
SLVTAB (itab)
P BS (or ZS) VS (KS ZSO VSO)
p bs (or zs) vs (ks zso vso)
(at least two data cards are required here)
00 Definitions:
If Format (2) in the ALPHA card is used, the user must enter a miscibility
pressure table immediately after each solvent table. The format for the
miscibility pressure table is described on the following card. 00
MISPTB (itab)
Z1
z11 z12 ... z1n
(at least two z1 values are required here)
Z3 z3
MISP
misp1 misp2 ... mispn
(the MISP lines are repeated for each new Z3 line)
00 Alternatively:
MISPTB (itab)
Y1
y11 y12 ..... y1n
(at least two y1 values are required here)
MISP
misp1 misp2 ........ mispn
MISPTB (itab)
Z1
z11 z12 ....... z1n
(at least two z1 values are required here)
MISP
misp1 misp2 ........mispn
00 Definitions:
At least two values of z1 (or y1) must be entered and the values of z1 (or
y1) must be monotonically increasing. In three-component cases, no gas is
allowed and misp is simply a function of solvent mole fraction (z1). In
four-component cases, when gas (z1) and solvent (z3) mole fractions are
both specified, misp values must be input for each combination (one set
for each z3 value). The values of z3 must also be monotonically increasing.
Alternatively, the gas pseudo-phase (y1) can be specified. The number of
misp values must be equal to the number of z1 values in three-component
cases. Values not specified in four-component cases will be filled with the
last misp value entered. 00
MWS mws
00 Definition:
SIGT itab
P SIGR (DGOG)
p sigr (dgog)
. . .
. . .
. . .
(At least two data cards are required.)
00 Definitions:
σ(p)
sigr ( p ) = ---------------
σ ( po )
σ = surface tension
po = pressure used to generate the matrix capillary
pressure curve entered in the SGT table, see Section
4.3.1.2. This is usually the initial reservoir pressure.
00 Example:
00 SIGT 1
P SIGR DGOG
14.7 100.0 .7795
1688.7 66.6667 .5934
2044.7 52.2222 .5666
2544.7 36.6607 .5299
3005.7 24.4444 .4972
3567.7 14.2222 .4597
4124.7 8.0 .4247
4558.7 4.9333 .3986
4949.7 2.8333 .3760
5269.7 1.7222 .3570
5559.7 1.0 .3411
7014.7 0.5556 .2676
00 Definitions:
NOTE: Rock types with no specified values of betag will default to betag
= 1.0.
HYDBETA
a1 a2 ... anf*NZ
b1 b2 ... bnf*NZ
00 Definitions:
NOTE:
b
Beta = ----------------------------------------------a-
( tamult ( P ) ⋅ k D )
KX
where KD = ------------ Darcies.
1000
■ separator batteries
CORE for the automatic generation of the composition entries. The fluid
compositions calculated in different stages of simulated PVT tests are
automatically included in the list of the composition entries. Also, these
compositions are output in a file (Fortran Unit 77) in the format of the
CMP card. Liquid and vapor densities, viscosities, and compressibility
factors in different stages of the differential expansion and multiple
contact tests are also included in this file.
00 The user can divide the composition entries into paths. Path numbers can
be defined in the CMP, SWELLTEST, DIFEXPTEST, CONVDPTEST, and
MULCONTEST cards.
nc
F(Z) = ∑ ( ci × zi ) , (Eq. 14-1)
i=1
nc
∑ ( zi – zi )
j 2
D ( Z j, Z ) = × ci .
i=1
00 The MDI keyword must be included on the EOSINT card to apply the
multidimensional interpolation procedure based on the distance. The path
definitions are ignored in this case and the interpolation between
compositional paths is not applied.
nc
Fp ( Z ) = ∑ ( cpi × zi ) .
i=1
00 The coefficients c p1, c p2, …, c pnc can be input by the user on the
00 All input cards of the EOSINT option must be input in the TABLES section
after the EOSINT card. All input cards required for the simulation of a
PVT test must be input after the corresponding DIFEXPTEST,
CONVDPTEST, SWELLTEST, or MULCONTEST card. Several PVT tests
can be simulated.
00 The CPU memory requirements for the EOS interpolation option are
determined internally. Therefore, input on the DIM card is not required.
Definitions:
00 Example:
00 C
C The EOS interpolation option is to be used for the phase-equilibrium
C calculations in separator batteries and well tubing strings
C (surface facility system).
C
EOSINT NORESERVOIR
TEMPERATURE t1 t2 …
Definition:
00 Example
C
C Temperature entries of the EOS interpolation tables.
C
TEMPERATURE 100 150 180
PMAX pmax
Definition:
NOTE: 12,000 psia is applied as the maximum pressure entry of the EOS
interpolation tables if the PMAX card is not included.
00 Example:
C
C Maximum pressure entry of the EOS interpolation tables.
C
PMAX 4600
PMIN pmin
Definition:
00 Example:
C
C Minimum pressure entry of the EOS interpolation tables.
C
PMIN 940
ALL
CMP (ipath)
iequil
psat1
1
z1
1
z2 … 1
z nc
psat2
2
z1
2
z2 … 2
z nc
.
.
.
psatnz
nz
z1
nz
z2 … nz
z nc
Definitions:
00 Example
C
C Composition entries of the EOS interpolation table.
C
CMP 1 1
C psat Z
2235.24 .4892781 .0303435 .0714230 .2044471 .1533799 .0511285
2163.28 .4783313 .0306918 .0728783 .2089891 .1568294 .0522801
2091.33 .4671521 .0310449 .0743667 .2136294 .1603513 .0534556
2019.37 .4557324 .0314028 .0758894 .2183714 .1639480 .0546560
1947.41 .4440639 .0317654 .0774476 .2232186 .1676224 .0558822
COEFFICIENTS (PATH)
c1 c2 … cn
d d2 … d n
1
Definitions:
∑ ci zi
F(z) = i=1
-------------------
n
∑ di zi
i=1
00 Example:
C
C Coefficients of the path interpolation function.
C
COEFFICIENT PATH
0. 10. 0. 0. 0. 0.
ITNMAX itnmax
Definition:
00 Example:
C
C Number of outer iterations in
C which recovery and compressibility factors are recalculated.
C
ITNMAX 4
DELTA delta
00 Definition:
00 Example
C
DELTA 0.00001
ROILCM
RGASCM
OUTPUT (iequil isepnum istnum)
SOILCM
SGASCM
Definitions:
00 Example:
C
C Output of oil compositions in the entries of the EOS interpolation
C table in Equilibrium Region 2.
C
OUTPUT ROILCM 2
C
C Output of oil compositions in the entries of the EOS interpolation
C table in oil sales line of Separator Battery 2.
C
OUTPUT SOILCM 1 2 0
ALL
SWELLTEST (ipath)
iequil
OILCM x 1 x 2 … x nc
GASCM y 1 y 2 … y nc
(GASFRAC f 1 f 2 … f nf )
(NGASFR ngasfr)
(PSAT psat)
Definitions:
NOTE: If GASFRAC and NGASFR data are both entered, the NGASFR
data is ignored.
00 Example:
C
C Swelling test simulation
C
SWELLTEST ALL 1
C Oil Composition
OILCM .50 .03 .07 .20 .15 .05
C Gas Composition
GASCM .77 .20 .03 0.0 .0 .0
OIL
ALL
DIFEXPTEST GAS (ipath)
iequil
BOTH
COMP z 1 z 2 … z nc
(PRINPS prinps)
(NPRNPS nprnps)
(PRES p 1 p 2 … p np
(NPRES npres)
(GASFRAC f 1 f 2 … f nf )
(NGASFR ngasfr)
(PSAT psat)
(TEMP t t … t )
1 2 nt
Definitions:
NOTE: If PRES and NPRES data are both entered, the NPRES data is
ignored. If GASFRAC and NGASFR data are both entered, the
NGASFR data is ignored.
00 Example:
C
C Differential expansion test simulation
C
DIFEXPTEST 1 BOTH 1
C Number pressure entries
NPRES 20
C Fluid composition in the differential expansion test
COMP .50 .03 .07 .20 .15 .05
Number of gas fractions in the swelling test.
NGASFR 10
OIL
ALL GAS
CONVDPTEST (ipath)
iequil BOTH
TOTAL
COMP z 1 z 2 … z nc
(PRINPS prinps)
(NPRNPS nprnps)
(PRES p 1 p 2 … p np )
(NPRES npres)
(GASFRAC f f … f )
1 2 nf
(NGASFR ngasfr)
(PSAT psat)
(TEMP t t … t )
1 2 nt
Definitions:
NOTE: If PRES and NPRES data are both entered, the NPRES data is
ignored. If GASFRAC and NGASFR data are both entered, the
NGASFR data is ignored.
00 Example:
00 C
C Constant volume depletion test simulation
C
CONVDPTEST 1 TOTAL 1
C Number pressure entries
NPRES 25
C Fluid compositon
COMP .50 .03 .07 .20 .15 .05
PRINPS 20.
NPRNPS 11
FORWARD
ALL
MULCONTEST BACKWARD (ipath)
iequil
BOTH
OILCM x x … x
1 2 nc
GASCM y y … y
1 2 nc
NFLASH nflash
(PRES p p … p )
1 2 np
(NPRES npres)
(GASFRAC f f … f )
1 2 nf
(NGASFR ngasfr)
(PSAT psat)
(TEMP t t … t )
1 2 nt
Definitions:
NOTE: If PRES and NPRES data are both entered, the NPRES data is
ignored. If GASFRAC and NGASFR data are both entered, the
NGASFR data is ignored.
00 Example
C
C Multiple contact test simulation.
C
MULCON ALL BACKWARD 1
C Oil Composition
OILCM .50 .03 .07 .20 .15 .05
C Injected Gas Composition
GASCM .77 .20 .03 0.0 .0 .0
C Number of the stages in the PVT test
NFLASH 10
C Pressure in the PVT test
PRES 2000
C Gas fraction in the PVT test
GASFRAC 0.525
Ce
φmax
FR = B/A
Ce 2
3
Ce
Cr
A 4 Cd
B Ce
φo 1
R I
Pbase P Pd
r
The default value of Pbase shown in Figure 4-10 is PINIT for equilibrium
region 1 in the EQUIL table (Section 4.2). If PINIT(1) is greater than Pd,
then Pbase must be reset to some value less than Pd as shown below.
Definitions:
KMUL ⋅ ( φ – φ 0 )
k = k 0 ⋅ EXP -------------------------------------------
( 1 – φ0 )
Values of KMUL are specified as array data (Section 5.46) for each flow
direction (including diagonal directions for NINEPT option) in such a
manner as to define the fracture plane.
Rock heat capacity table numbers are assigned to gridblocks with the
ICPRTB array which is described in Section 5.47.
CPRTAB icprtb
T CPR
t1 cpr1
t2 cpr2
. .
. .
. .
Definitions:
ti Temperature, ˚F (˚C).
00 nreactants
Rate = V A e(-Ea / (R T)) Π (Ci)ni
i=1
00 The phase volume option should be used for, and is allowed only for,
homogeneous reactions. The fluid volume option is recommended for
heterogeneous reactions. The total volume option is not recommended
except to match results from other simulators which may use this
representation. This option results in reactant concentrations which
depend on porosity, which is intuitively incorrect. Adding rock volume to
a system while maintaining fluid volume and saturations constant should
not change the reaction rate. This is not true in the total volume option
except for first-order reactions (reaction order is equal to the sum of
reactant orders).
00 The input data allow specification of model parameter units for ease of
use. The units of rate are moles of first reactant consumed per unit time
(MUNIT/TUNIT), concentrations are in moles of reactant per unit volume
(MUNIT/VUNIT), molar energies (activation energy and heat of reaction,
units may be independently specified) are in energy per mole (EUNIT).
00 Any of the above units given for any quantity may be specified for each
reaction. VUNIT may be additionally specified as barrels or cubic
centimeters. EUNIT may be additionally specified as calories/g mole. The
units of the rate constant A are (VUNIT / MUNIT)n-1/ TUNIT where n is
the order of reaction.
REACTION nr (card 1)
si cmpi PLUS sj cmpj PLUS ... = sk cmpk PLUS sl cmpl PLUS ...
(card 2)
ORDER ni nj ... (card 3)
PHASE W W ... (card 4)
O O ...
G G ...
VOLUME PHASE (VUNIT) (card 5)
FLUID
TOTAL
A a (card 6)
EACT eact (EUNIT) (card 7)
HR hr (EUNIT) (card 8)
( MOLES MUNIT ) (card 9)
( TIME TUNIT ) (card 10)
00 Definitions:
W water,
O oil, or
G gas.
BBL barrels,
M3 cubic meters,
L liters, or
CC cubic centimeters.
LBM lb moles,
KGM kg moles, or
GM g moles.
DAYS days, or
HRS hours.
NOTE: 1. Reaction numbers must be greater than zero and less than or
equal to the total number of reactions specified.
5
00000Grid Data Arrays
Arrays described in this section are input using the array input options
described in Section 1.5.2, and can be modified using the array
modification cards described in Sections 1.5.4 and 6.1. 00
In VIP-THERM, only combinations 1-4 are valid, and in each of those the
following arrays are additionally required: CPR0, DTX0 (or DTR0), DTY0
(or KTTH0) and KTZ0. 00
KTHETA KZ
R DZNET KR DEPTH
1) DTHETA DZ POR
DR NETGRS KRF KTF KZF MDEPTH
DZNET KX KYF KZ DEPTH
2) DX DY DZ POR
NETGRS KXF KYF KZF MDEPTH
4) DZCORN
XCORN YCORN ZCORN DZBCOR ( ZCORNE ) ( ZCORNW ) ( ZCORSW ) ( NETGRS ) KX KY KZ POR
DZVCOR
5) DZCORN
XCORN YCORN ZCORN DZBCOR ( ZCORNE ) ( ZCORNW ) ( ZCORSW ) ( NETGRS ) TX TY TZ PV
DZVCOR
Certain arrays must be read before other arrays. For example, the
thickness array must precede the depth array. Order restrictions are
specified in the array description. 00
When using the DUAL option (Section 2.2.13.1) these arrays represent the
matrix properties. Fracture properties are input using the following sets of
keywords: 00
TEX
1) R ... KRFEFF KTFEFF KZFEFF PORF DEPF
SIGMA
DR MDEPF
LX LY LZ
TEX
2) DX ... SIGMA KRFEFF KTFEFF KZFEFF PORF DEPF
MDEPF
LX LY LZ
TEX
3) DXB ... SIGMA KRFEFF KTFEFF KZFEFF PORF DEPF
MDEPF
LX LY LZ
TEX
4) XCORN ... POR SIGMA KRFEFF KTFEFF KZFEFF PORF DEPF
MDEPF
LX LY LZ
ARRAYS
Definitions: 00
Certain arrays must be read before other arrays. For example, the
thickness array must precede the depth array. Order restrictions are
specified in the individual array writeups. 00
Examples: 00
ARRAYS
C
C GRID BLOCK DIMENSIONS
C
DX CON
66.67007
C
DY CON
66.67007
C
DZ ZVAR
7.2 4.4 15.8
C
C GRID BLOCK TOPS
C 00
These grid dimensions are usually input using the XVAR (RVAR) option. 00
Definitions: 00
Rectangular Grid: 00
Radial Grid: 00
Example: 00
ARRAYS
C
DX XVAR
100 150 200 00
Definition: 00
Example: 00
ARRAYS
XCORN VALUE
C FIRST LAYER
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
C SECOND LAYER
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.2
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
0.00 956.52 1913.04 2869.57 3826.09 4782.61 5739.13 6695.65 7652.17
8608.70 9565.22
C THIRD LAYER
0.00 886.96 1773.91 2660.87 3547.83 4434.78 5321.74 6208.70 7095.65
7982.61 8869.57
0.00 886.96 1773.91 2660.87 3547.83 4434.78 5321.74 6208.70 7095.65
7982.61 8869.57
0.00 886.96 1773.91 2660.87 3547.83 4434.78 5321.74 6208.70 7095.65
7982.61 8869.57
0.00 886.96 1773.91 2660.87 3547.83 4434.78 5321.74 6208.70 7095.65
7982.61 8869.57
0.00 886.96 1773.91 2660.87 3547.83 4434.78 5321.74 6208.70 7095.65
7982.61 8869.57
0.00 886.96 1773.91 2660.87 3547.83 4434.78 5321.74 6208.70 7095.65
7982.61 8869.57 00
These grid dimensions are usually input using the YVAR (THVAR)
option. 00
Definitions: 00
Rectangular Grid: 00
Radial Grid: 00
Definition: 00
Example: 00
ARRAYS
YCORN VALUE
C FIRST LAYER
10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000
8000 8000 8000 8000 8000 8000 8000 8000 8000 8000 8000
7000 7000 7000 7000 7000 7000 7000 7000 7000 7000 7000
6000 6000 6000 6000 6000 6000 6000 6000 6000 6000 6000
5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000
4000 4000 4000 4000 4000 4000 4000 4000 4000 4000 4000
3000 3000 3000 3000 3000 3000 3000 3000 3000 3000 3000
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
00000000000
C SECOND LAYER
8021.74 8021.74 8021.74 8021.74 8021.74 8021.74 8021.74 8021.74 8021.74
8021.74 8021.74
7130.43 7130.43 7130.43 7130.43 7130.43 7130.43 7130.43 7130.43 7130.43
7130.43 7130.43
6239.13 6239.13 6239.13 6239.13 6239.13 6239.13 6239.13 6239.13 6239.13
6239.13 6239.13
5347.83 5347.83 5347.83 5347.83 5347.83 5347.83 5347.83 5347.83 5347.83
5347.83 5347.83
4456.52 4456.52 4456.52 4456.52 4456.52 4456.52 4456.52 4456.52 4456.52
4456.52 4456.52
3565.22 3565.22 3565.22 3565.22 3565.22 3565.22 3565.22 3565.22 3565.22
3565.22 3565.22
2673.91 2673.91 2673.91 2673.91 2673.91 2673.91 2673.91 2673.91 2673.91
2673.91 2673.91
1782.61 1782.61 1782.61 1782.61 1782.61 1782.61 1782.61 1782.61 1782.61
1782.61 1782.61
891.30 891.30 891.30 891.30 891.30 891.30 891.30 891.30 891.30 891.30
891.30 891.30
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0 0
C THIRD LAYER
6456.52 6456.52 6456.52 6456.52 6456.52 6456.52 6456.52 6456.52 6456.52
6456.52 6456.52
5739.13 5739.13 5739.13 5739.13 5739.13 5739.13 5739.13 5739.13 5739.13
5739.13 5739.13
5021.74 5021.74 5021.74 5021.74 5021.74 5021.74 5021.74 5021.74 5021.74
5021.74 5021.74 00
4500 4500 4500 4500 4500 4500 4500 4500 4500 4500 4500
4000 4000 4000 4000 4000 4000 4000 4000 4000 4000 4000
3500 3500 3500 3500 3500 3500 3500 3500 3500 3500 3500
3000 3000 3000 3000 3000 3000 3000 3000 3000 3000 3000
2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500
2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000
1500 1500 1500 1500 1500 1500 1500 1500 1500 1500 1500
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
500 500 500 500 500 500 500 500 500 500 500
0 0 0 0 0 0 0 0 0 0 0 00
Gross thickness usually is input using the ZVAR or VALUE option. It must
precede the DEPTH or MDEPTH array (non-corner-point grids). 00
Definition: 00
Gross thickness usually is input using the ZVAR or VALUE option. It must
precede the DEPTH or MDEPTH array (non-corner-point grids). 00
Definition: 00
Example: 00
ARRAYS
DZB ZVAR
1.0 2.0 3.0 00
Definition: 00
This gross thickness array must precede the ZCORN array. This array
requires (NX+1) * (NY+1) * NZ values and is used only to process the
LAYER or DIP options for the ZCORN array. 00
Definition: 00
This gross thickness array is usually input using the VALUE option. It
must precede the ZCORN array. This array requires (NX+1) * (NY+1) * NZ
values and is used only to process the LAYER or DIP options for the
ZCORN array. 00
Definition: 00
This gross thickness array is usually input using the VALUE option. It
must precede the ZCORN array. This array requires (NX+1) * (NY+1) * NZ
values and is used only to process the LAYER or DIP options for the
ZCORN array. 00
Net thickness usually is input using the ZVAR or VALUE option. The net
pay is assumed to be evenly distributed across the entire gross interval. If
both gross and net thicknesses are input, gross thickness is used only in
the calculations of depth and vertical transmissibility; net thickness is used
to determine pore volume and areal transmissibility. 00
Definition: 00
Net thickness usually is input using the ZVAR or VALUE option. The net
pay is assumed to be evenly distributed across the entire gross interval. If
both gross and net thicknesses are input, gross thickness is used only in
the calculations of depth and vertical transmissibility; net thickness is used
to determine pore volume and areal transmissibility. 00
Definition: 00
Definition: 00
Net/Gross ratio usually is input using the ZVAR or VALUE option. Net/
Gross ratio is used to calculate net from gross thickness. The net pay is
assumed to be evenly distributed across the entire gross interval. Gross
thickness is used in the calculations of depth and vertical transmissibility;
net thickness is used to determine pore volume and areal transmissibility. 00
Definition: 00
Net/Gross ratio usually is input using the ZVAR or VALUE option. Net/
Gross ratio is used to calculate net from gross thickness. The net pay is
assumed to be evenly distributed across the entire gross interval. Gross
thickness is used in the calculations of depth and vertical transmissibility;
net thickness is used to determine pore volume and areal transmissibility. 00
Definition: 00
Definition: 00
Example: 00
MDEPTH LAYER
20*120.0 00
Since the reservoir simulator deals with block properties and interblock
connections only, any system which can be described in these terms can be
modelled. Corner point gridding permits the user to define the eight
corners of each gridblock and generates the required simulation
properties. The user is cautioned that generation of non-orthogonal
connections between gridblocks can lead to errors in flow calculations,
because the resulting cross terms in the equations are ignored. 00
Definition: 00
Definition: 00
Example: 00
ZCORNE ZVAR
100 200 250 00
Definition: 00
Example: 00
ZCORNW LAYER
20*110.0 00
Definition: 00
Example: 00
ZCORSW ZVAR
50*200.0 00
Definition: 00
Example: 00
ZBOT ZVAR
15.3 18 00
Definition: 00
Example: 00
ZBOTNE ZVAR
15.3 18 00
Definition: 00
Example: 00
ZBOTNW ZVAR
15.3 18 00
Definition: 00
Example: 00
ZBOTSW ZVAR
15.3 18 00
Definition: 00
Example: 00
ZLNCOR ZVAR
15.3 18 00
Definitions: 00
Definitions: 00
Example: 00
POR CON
0.3 00
Definitions: 00
The PORF array must be used with the POR array. The PVF array must be
used with the PV array. 00
Definitions: 00
Rectangular Grid:
Radial Grid:
Example: 00
KXF ZVAR
10 20 25 00
Definitions: 00
Rectangular Grid:
Radial Grid:
Definitions: 00
The TXYL and TXYR arrays may only be entered if the nine-point finite
difference operator option is in use (Section 2.2.6.5). 00
No corner points
NX > 1, NY > 1
Definitions: 00
k = kf * φf
array option
Enter data values as required.
Definitions: 00
Rectangular grid:
Radial grid:
array option
Enter data values as required.
Definitions: 00
Rectangular grid:
Radial grid:
array option
Enter data values as required.
Definitions: 00
Rectangular grid:
array option
Enter data values as required.
The TXYLF and TXYRF arrays may only be entered if the nine-point finite
difference operator option is in use (Section 2.2.6.5). 00
No corner points
NX > 1, NY > 1
Definitions: 00
Definition: 00
Definition: 00
The ISATF array is only required when the DUAL option is in use and the
user requires different numbered rock types in the fracture than the
matrix. In this case, both the ISAT and the ISATF arrays are entered. If the
ISATF array is omitted the rock type for the fractures defaults to the rock
type for the matrix. This does not, however, mean that by default the same
saturation tables are used for the fracture and matrix but just that if a
gridblock uses matrix saturation table 1 (SWT and SGT) it also uses
Definition: 00
The ISATIF array is only required when the DUAL option is in use and the
user requires different numbered imbibition rock types in the fracture than
the matrix. In this case, both the ISATI and the ISATIF arrays are entered.
If the ISATIF array is omitted the rock type for the fractures defaults to the
rock type for the matrix. This does not, however, mean that by default the
same saturation tables are used for the fracture and matrix but just that if a
gridblock uses matrix saturation table 1 (SWT and SGT) it also uses
fracture saturation table 1 (SWTF and SGTF). 00
Definition: 00
The IPVT array is required if more than one set of PVT property tables is
input (see Section 4.5). At present the IPVT array is inactive in VIP-COMP
and VIP-THERM. It is recommended that the IPVT number should be
entered on the IEQUIL card instead of entering this array (Section 4.2.). 00
Definition: 00
The IPVTW array is required if more than one set of water property tables
(PVTW, Section 4.11.1 or PVTWSAL, Section 4.11.2) is input. It is
recommended that the IPVTW number should be entered on the IEQUIL
card (Section 4.2) instead of entering this array. 00
Definition: 00
00
Definition: 00
The extra region (XREG) option may be used to assign gridblocks to more
than one output region. See Section 5.16. 00
Example: 00
C--------------------------------------------------------
C ARRAY DATA
C--------------------------------------------------------
ARRAYS
...................
IREGION ZVAR
1 2 3
MOD
35 37 2 2 1 1 =4 00
34 37 3 3 1 1 =4
34 37 4 4 1 1 =4 00
The XREG cards must follow the IREGION array specification and any
MOD or VMOD cards that apply to it. These extra regions are identified
by an X next to the region number on the Initial Fluids in Place report in
VIP-CORE and the Region Summary reports in the simulation modules.
The totals printed on these reports do not include the calculated values of
the extra regions. 00
Only one title card containing the keyword XREG is required, but the data
cards may be repeated as necessary. 00
XREG
i1 i2 j1 j2 k1 k2 (= v)
Definitions: 00
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
k1 ≤ K ≤ k2 00
Examples: 00
IREGION ZVAR
1 2 3 4 5
XREG
1 4 2 3 1 2 =6
1 4 2 3 3 5 =7
1 2 1 2 1 5 =8
1 2 3 4 1 5 =9 00
XREG
1 4 2 3 1 2
6 6 6 6 6 6 6 6
7 7 7 7 7 7 7 7 00
The IREGF array is only required when the DUAL option is in use and the
user requires different numbered output regions in the fracture than the
matrix. In this case, both the IREGION and the IREGF arrays are entered.
If the IREGF array is omitted, the output region for the fractures defaults
to the output region for the matrix. 00
Definition: 00
The extra region (XREG) option may be used to assign fracture blocks to
more than one output region by using the XREG card after the IREGF
array. See Section 5.16. 00
Definition: 00
The ICMT array is required if more than one set of COMPACTION tables
is input (see Section 4.12). 00
Definition: 00
The ICMT array may be entered only if the COMPACT card is entered in
the utility data (Section 2.2.7.10). 00
The ICMTF array is only required when the DUAL option is in use and the
compaction properties for the fracture are different from those for the
matrix. In this case, both the ICMT and the ICMTF arrays are entered. If
the ICMTF array is omitted, the compaction table used for each fracture
block will be the same as the table used for the corresponding matrix
block. 00
Definition: 00
The IWIRC array is required if more than one set of WATER INDUCED
ROCK COMPACTION tables are input (Section 4.13). 00
Definition: 00
The IWIRCF array is only required when the DUAL option is in use and
the water induced rock compaction properties for the fracture are different
from those for the matrix. In this case, both the IWIRC and IWIRCF arrays
are entered. 00
Definition: 00
Definition: 00
The CRF array is only required when the DUAL option is in use and the
fracture compressibilities are different from the matrix compressibilities.
In this case, both the CR and the CRF arrays are entered. If the CRF array
is omitted, the fracture compressibilities will be set equal to the
corresponding matrix compressibilities. 00
Definition: 00
Definition: 00
Definition: 00
The turbidite reservoir option is compatible with both the IMPES and
IMPLICIT formulations in the simulation modules. It can be used in both
regular Cartesian and radial grid systems. The use of this option in a
corner point grid system is not recommended because of its irregular
geometry. The option should not be used in conjunction with the dual
porosity option and is not yet compatible with the water tracking option. 00
Definition: 00
Examples: 00
SCLFCT CON
0.
MOD
6 10 1 3 1 1 =6
1 5 1 2 2 2 = 10. 00
6 10 1 3 3 3 = 20. 00
Definitions: 00
where
The TCTBD array must be specified if the SCLFCT array is entered. Non-
zero TCTBD values must be entered for all gridblocks with non-zero
SCLFCT values. Gridblocks with zero SCLFCT values will be treated as
non-turbidite gridblocks without sand-shale fluid exchange, regardless of
their TCTBD values. 00
Examples: 00
TCTBD ZVAR
6.7645E-3 1.6911E-3 6.7645E-3 0. 00
Definitions: 00
where
The BTBD array must be specified if the SCLFCT array is entered. Non-
zero BTBD values must be entered for all gridblocks with non-zero
SCLFCT values. Gridblocks with zero SCLFCT values will be treated as
non-turbidite gridblocks without sand-shale fluid exchange, regardless of
their BTBD values. 00
Examples: 00
BTBD ZVAR
3.3840E-3 6.7680E-3 3.3840E-3 0. 00
The IEQUIL array is required if more than one set of fluid contacts is input
(see Section 4.2). 00
Definition: 00
Definition: 00
The SAL array may be entered only if the PVTWSAL data has been
entered (Section 4.11.2). 00
Definitions: 00
The mole fraction for each component in the gas phase for each gridblock
may be input. The format of the array input is very similar to other array
data input. 00
Definitions: 00
NOTE: If the mole fractions for any component are overread, then
the mole fractions must be overread for all components.
The mole fraction for each component in the oil phase for each gridblock
may be input. The format of the array input is very similar to other array
data input. 00
Definitions: 00
NOTE: If the mole fractions for any component are overread, then
the mole fractions must be overread for all components.
Definitions: 00
Example: 00
PF CON
20 00
The mole fraction for each component in the fracture gas phase for each
gridblock may be input. The format of the array input is very similar to
other array data input. 00
Definitions: 00
NOTE: If the mole fractions for any components are overread, then
the mole fractions must be overread for all components.
Example: 00
The mole fraction for each component in the fracture oil phase for each
gridblock may be input. The format of the array input is very similar to
other array data input. 00
Definitions: 00
NOTE: If the mole fractions for any component are overread, then
the mole fractions must be overread for all components.
All of the above arrays may be entered with the directional relative
permeability option. See Section 1.5.2.9 for details. 00
All of the above arrays are duplicated for DUAL cases, with an "F"
appended to the array name (e.g., SWL becomes SWLF). 00
The user can select the two-point scaling of relative permeability table
endpoints, or the default three-point scaling (four-point for capillary
pressures), by use of the END2P option. The two-point scaling approach
was the only method available in VIP-EXECUTIVE Version 1.6R, and
earlier versions. In the two-point case, each curve is scaled over its entire
length (from residual/irreducible saturation to the saturation at which it
attains a maximum). In the three-point case, all curves in a table are scaled
together (retaining the relative kr and Pc characteristics of the curves). In
this case, which is the default scaling method, all endpoints serve to break
the table up into partitions, with scaling being done in each partition
independently. As an example, water-oil table relative permeabilities
would be scaled in two sections; from water saturations of Swr to Swro
and from Swro to 1. 00
If any saturation endpoint is not specified, rock data will be used. The
following consistency checks are performed for each gridblock: 00
NOTE: If the SWR array is entered, the SWRO array must also be
entered.
NOTE: If the SWRO array is entered, the SWR array must also be
entered.
NOTE: 1. The SWLF array is required only when the DUAL option
is in use and the user requires different values than
would be derived from the appropriate fracture
saturation tables (based on the ISATF array).
NOTE: 1. The SWRF array is required only when the DUAL option
is in use and the user requires different values than
would be derived from the appropriate fracture
saturation tables (based on the ISATF array).
NOTE: 1. The SWUF array is required only when the DUAL option
is in use and the user requires different values than
would be derived from the appropriate fracture
saturation tables (based on the ISATF array).
NOTE: If the SGR array is entered, the SGRO array must also be
entered.
NOTE: If the SGRO array is entered, the SGR array must also be
entered.
NOTE: 1. The SGLF array is required only when the DUAL option
is in use and the user requires different values than
would be derived from the appropriate fracture
saturation tables (based on the ISATF array).
NOTE: 1. The SGRF array is required only when the DUAL option
is in use and the user requires different values than
would be derived from the appropriate fracture
saturation tables (based on the ISATF array).
NOTE: 1. The SGUF array is required only when the DUAL option
is in use and the user requires different values than
would be derived from the appropriate fracture
saturation tables (based on the ISATF array).
Definition: 00
NOTE: The FVEWO array may be entered only if the VEWO card is
entered in the utility data (Section 2.2.10.1).
Definition: 00
NOTE: The FVEGO array may be entered only if the VEGO card is
entered in the utility data (Section 2.2.10.2).
Definition: 00
Definition: 00
Definitions: 00
Example: 00
LX CON
100
Definition: 00
Definition: 00
Example: 00
TEX CON
0.1
Definition: 00
Definition: 00
Example: 00
SIGMAD CON
0.8 00
Definition: 00
Definition: 00
Definition: 00
Definition: 00
Definition: 00
Definition: 00
The OMGV array allows the user to input variable values for viscosity
mixing parameter in miscible option (invoked by the MIS card). This array
overrides the omegav value specified in the MIS card of the Initialization
Data section. The omegav value will be used if this card is omitted. 00
Example: 00
OMGV ZVAR
0.5 0.6 0.7 00
Definition: 00
The OMGD array allows the user to input variable values for density
mixing parameter in miscible option (invoked by the MIS card). This array
overrides the omegad value specified in the MIS card of the Initialization
Data section. The omegad value will be used if this card is omitted. 00
Example: 00
OMGD ZVAR
0.5 0.6 0.7 00
Definition: 00
Example: 00
CREEPB ZVAR
6.169 1.227 3.000
Definition: 00
Example: 00
CREEPC CON
1.E-5
Definition: 00
Example: 00
CREEPM ZVAR
17.1 15.65 16.3
Transmissibilities and pore volumes can also be modified using OVER and
VOVER cards. Separate options are available in PRINT ARRAYS output or
in the initialization MAP file to report transmissibility multiplier values
before and after the effects of OVER/VOVER cards. 00
Definition: 00
NOTE: Omit if NX = 1.
Definition: 00
NOTE: Omit if NY = 1.
Definition: 00
NOTE: Omit if NZ = 1.
Definition: 00
NOTE: Omit if NR = 1.
Definitions: 00
Definition: 00
Definition: 00
The TMXF array is only required when the DUAL option is in use and the
fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMX and the TMXF arrays are entered. If
the TMXF array is omitted, the multipliers will be set equal to the
corresponding matrix multipliers. 00
Definition: 00
The TMYF array is only required when the DUAL option is in use and the
fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMY and the TMYF arrays are entered. If
the TMYF array is omitted, the multipliers will be set equal to the
corresponding matrix multipliers. 00
Definition: 00
The TMZF array is only required when the DUAL option is in use and the
fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMZ and the TMZF arrays are entered. If
the TMZF array is omitted, the multipliers will be set equal to the
corresponding matrix multipliers. 00
Definition: 00
The TMRF array is only required when the DUAL option is in use and the
fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMR and the TMRF arrays are entered. If
the TMRF array is omitted, the multipliers will be set equal to the
corresponding matrix multipliers. 00
Definition: 00
The TMTHF array is only required when the DUAL option is in use and
the fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMTH and the TMTHF arrays are
entered. If the TMTHF array is omitted, the multipliers will be set equal to
the corresponding matrix multipliers. 00
Definition: 00
The TMXYLF array is only required when the DUAL option is in use and
the fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMXYL and the TMXYLF arrays are
entered. If the TMXYLF array is omitted, the multipliers will be set equal
to the corresponding matrix multipliers. 00
Definition: 00
The TMXYRF array is only required when the DUAL option is in use and
the fracture transmissibility multipliers are different from the matrix
multipliers. In this case both the TMXYR and the TMXYRF arrays are
entered. If the TMXYRF array is omitted, the multipliers will be set equal
to the corresponding matrix multipliers. 00
Definition: 00
Definition: 00
NOTE: Omit if NX = 1.
Definition: 00
NOTE: Omit if NY = 1.
output generated by PRINT ARRAYS and in the MAP file. The values of
this array will not be modified if thermal transmissibilities are further
modified by OVER or VOVER cards. 00
Definition: 00
NOTE: Omit if NZ = 1.
Definition: 00
NOTE: Omit if NR = 1.
Definitions: 00
00
The entries of the function table can be input in arbitrary order. However,
they are internally sorted for each grid block in increasing order of a
distance between values of the input variables in the grid block and in the
function table. The first m entries of the sorted function table which satisfy
the constraints defined in the DRANGE card are used for the
interpolation. 00
The function procedure can be executed only for selected blocks. These
blocks are defined by the FUNCTION, BLOCKS, and/or RANGE INPUT
cards. 00
Input or output variables in the function option can be integer arrays like
ISAT, ISATI, IREGION, IPVT, ICMT, and/or IEQUIL. They are internally
transformed in a real data type and, then, the standard function option is
applied. 00
Note that when IREGION is used as an input or output array, and extra
regions have been defined, the extra region values are not involved in the
FUNCTION calculations. For IREGION as input, only the original region
number is used in the function table lookup. For IREGION as output, only
the original region number is replaced; the extra region numbers are
retained. 00
Analytical functions can be applied using the ANALYT card. In this case,
the function table entries, DRANGE and VOLAVR cards are ignored. In
the ABS, EXP, EXP10, LOG, LOG10, SQRT, SIN, COS, GE, LE, POLYN
analytical functions, only the first input array is used. All other input
arrays are ignored. In the DIV, MULT, ADD, SUBT, MIN, MAX analytical
functions, the first two input arrays are used. All other input arrays are
ignored. If the number of the input arrays in the DIV, MULT, ADD, SUBT,
MIN, MAX analytical functions is equal to one, the second input array is
assumed to be equal to the first input array. 00
Definitions: 00
p1...pN Shifts of the first, second, ... , and N-th input arrays.
They are deducted from the values of the input
variables in the function table. The default values
are zero.
i1 ≤ I ≤ i 2,
j1 ≤ J ≤ j 2,
k1 ≤ K ≤ k2.
i1 = j1 = k1 = 1,
i2 = NX, j2 = NY, k2 = NZ.
Every function table entry should be input on one line. It should consist of
N values of the input variables and K correspondent values of the output
variable. The function table entries can be input in arbitrary order. 00
NOTE: 1. FUNCTION cards must appear after all array input and
before any MULT card or fault data. Several FUNCTION
cards can be included in the data set.
Examples: 00
C Problem 1.
C Define the volume-averaged values of the permeability KX in
C Blocks (5,1,1)and (6,1,1) using its values in seven reservoir locations
FUNCTION
BLOCKS 5 6 1 1 1 1
VOLAVR
XC YC MDEPTH OUTPUT KX
4000 9000 8335 300
4210 9200 8335 304.455
4300 9300 8335 303.3
4400 9400 8335 304.4
4500 9600 8335 305.8
4600 9800 8335 306.6
5000 10000 8335 311
C 00
C Problem 2.
C Define the permeability KX as a function of the porosity POR
C and the critical water saturation SWR in grid blocks from the
C second and fourth output regions, in which the critical water
C saturation values
C are in the range [0.18,0.2].
C
FUNCTION IREGION
24
RANGE INPUT 0.18 0.2
SWR POR OUTPUT KX
0.1909 0.2 10.4767
0.1909 0.4 83.8135
0.1909 0.6 282.871
0.1909 0.8 670.508
C
C Problem 3.
C Define permeability KX as the following analytical function
C of porosity POR :
C KX = 10(17.06715 POR - 2.70903).
C Set permeability to zero, if its calculated value is less than one.
C
FUNCTION
ANALYT POLYN 17.06715 - 2.70903
POR OUTPUT KX
C
FUNCTION
ANALYT EXP10
KX OUTPUT KX
C
FUNCTION
ANALYT GE 1 0
KX OUTPUT KX 00
CPR0, Btu/FT3 ˚F
(KJ/M3˚C) Reference rock heat capacity of each gridblock.
KT = KT0 (1 - DKDSG*SG)
Rectangular Grid:
Radial Grid:
Rectangular Grid:
Radial Grid:
5.44.3 Z Direction
Omit if NZ=1
Rectangular Grid:
Radial Grid:
6
00000Fault Data
1. The standard Fault option assumes that the faults occur at gridblock
faces and that the fault throw (depth difference along the fault) is
constant in the z direction. (That is, if a fault occurs at the interface
between blocks (I-1,J,1) and (I,J,1), then it also occurs, with the same
throw, at the interfaces between blocks (I-1,J,K) and (I,J,K) for all K,
NP1 ≤ K ≤ NP2. NP1 and NP2 default to 1 and NZ, respectively, and
may be altered using the LAYERS card.) By using the corner-point
option, it is possible to model a sloping fault by defining the blocks
along the fault to have sloping faces. In this case, the fault is still
"logically vertical" since blocks in any column of the grid remain
adjacent only to blocks in one of the four surrounding (North, South,
East, or West) columns.
FAULTS
(FNAME fname)
Definitions:
NOTE: 1. FNAME identifiers can also be assigned using the OVER, and
VOVER keywords. A gridblock will be assigned based on the
last identifier encountered.
FX
FR i j1 (-j2) fshift (*tm *tmt)
or
FY
FTHETA j i1 (-i2) fshift (*tm *tmt)
or
FXCORN i j1 (-j2)
fshift(j
1) fshift(j
1+1) ... fshift(j 2)
(tm(j1) tm(j1+1) ... tm(j2-1))
(tmt(j1) tmt(j1+1) ... tmt(j2-1))
or
FYCORN j i1 (-i2)
fshift(i1) fshift(i 1+1) ... fshift(i2)
(tm(i1) tm(i1+1) ... tm(i2-1))
(tmt(i1) tmt(i1+1) ... tmt(i2-1))
Definitions:
Example: X (I) 00
/
/
/
/
//////// /
/
/////
/
/
/
FAULT
FX 3 1 -2 -20.1
FY 3 3 -4 20.1
FX 5 3 -15.5
FY 4 5 10. *1.5
FX 6 4 -5 -5 *2 00
or 00
FAULT
FXCORN 3 1 -3
-20.1
FYCORN 3 3 -5
20.1
FXCORN 5 3 -4
-15.5
FYCORN 4 5 -6
10 00
1.5
FXCORN 6 4 -6
-5
2 00
Definitions: 00
Example 1: 00
Example 2: 00
LEAKY ileaky
Definitions:
Example:
C Conductive fault #1
FX 3 1 -2 -20
LEAKY 1
FY 3 3 -4 20.
LEAKY 1
FX 5 3 -15.
LEAKY 1
FY 4 5 10.
LEAKY 1
C Non-Conductive fault
FX 10 2 -5 -10.
C Conductive-fault #2
FY 8 3 -8 10.
LAYERS
LEAKY 2
FLTXC ( TRNS0 )
FLTRC ( TRNS0 )
i j k1 k2 t tt
(repeat as necessary)
or
FLTYC ( TRNS0 )
FLTTC ( TRNS0 )
i j k1 k2 t tt
(repeat as necessary)
Definitions:
NOTE: 1. If the same fault connection is described more than once, the last
transmissibility value specified will be used.
FTRANS
( FNAME fname )
i1 j1 k1 i2 j2 k2 t (tt)
(repeat as necessary)
Definitions:
FLTXCF ( TRNS0 )
FLTRCF ( TRNS0 )
i j k1 k2 t
(repeat as necessary)
or
FLTYCF ( TRNS0 )
FLTTCF ( TRNS0 )
i j k1 k2 t
(repeat as necessary)
Definitions: 00
NOTE: 1. If the same fault connection is described more than once, the
last transmissibility value specified will be used.
FTRANF
( FNAME fname )
i1 j1 k1 i2 j2 k2 t
(repeat as necessary)
Definitions: 00
7
00000Overread Options
Arrays
TX TY TZ TR TTHETA TXYL TXYR PV
Definitions:
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
k1 ≤ K ≤ k2
+ add
- subtract
/ divide
* multiply
= equal
LE values larger than v will be set equal to v
GE values smaller than v will be set equal to v
NOTE: 1. The TXYL and TXYR arrays will be used only if the nine-point
option has been invoked (Section 2.2.6.5).
Examples:
OVER TZ
1 37 1 29 1 3 *0.
OVER TX TY TZ
12 56 18 *.1 *.1 *.05
Arrays
Definitions:
NOTE: 1. The TXYLF and TXYRF arrays will be used only if the nine-
point option has been invoked (Section 2.2.6.5).
VOVER array
(FNAME fname)
i1 i2 j1 j2 k1 k2 (op)
values as necessary
Arrays
Definitions:
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
k1 ≤ K ≤ k2
ADD add
SUB subtract
DIV divide
MULT multiply
EQ equal (this is the default).
NOTE: 1. Only one array can be changed with each VOVER card.
2. The TXYL and TXYR arrays will be used only if the nine-
point option has been invoked (Section 2.2.6.5).
Example:
VOVER TX
57 67 1 7 2 2 EQ
77*0.
VOVER array
(FNAME fname)
i1 i2 j1 j2 k1 k2 (op)
values as necessary
Arrays
Definitions:
NOTE: 1. Only one array can be changed with each VOVER card.
2. The TXYLF and TXYRF arrays will be used only if the nine-
point option has been invoked (Section 2.2.6.5).
8
00000Grid Coarsening1
COARSEN (gridname)
i1 i2 j1 j2 k1 k2 (nxc nyc nzc)
(X nx1 nx2 nx3 .........nxnxc )
(Y ny1 ny2 ny3 .........nynyc )
(Z nz1 nz2 nz3 .........nznzc )
( i1 i2 j1 j2 k1 k2 (nxc nyc nzc) )
.
.
.
( COARSEN (gridname) )
.
.
.
Definitions:
NOTE: 1) When X,Y,Z data are not specified, then nxc,nyc,nzc need to be
integer fractions of (i2-i1+1),(j2-j1+1),(k2-k1+1) respectively.
nxc
∑ nxi = i2 – i1 + 1
i=1
nyc
∑ nyi = j2 – j1 + 1
i=1
nzc
∑ nzi = k2 – k1 + 1
i=1
Example:
COARSEN
2 3 1 1 1 3 1 1 3
4 5 1 2 1 3 1 1 3
1 7 6 10 1 3 2 2 3
X 3 4
Y 3 2
11 20 11 20 1 3 5 5 3
COARSEN REFINE1
2 3 1 1 1 3 1 1 3
4 5 1 2 1 3 1 1 3
1 7 6 10 1 3 2 2 3
X 3 4
Y 3 2
11 20 11 20 1 3 5 5 3
9
00000Region Data
Any output region data must follow all grid and fault data and any
OVER/VOVER cards.
Definitions:
Example:
REGION 1 SAG
REGION 2 SHUB
REGION 3 SAD
REGION 4 WESAG
REGION 5 WESHUB
REGION 6 WESAD
Definitions:
NOTE: 1. The REGSEP card must follow the REGION cards and must
precede any REGDTM data.
dtm 1 dt m nreg
REGDTM . . .
X X
Definitions:
Example:
A problem with four output regions will use the default datum for
regions 2 and 3.
10
00000Grid Boundary Flux 1
10.1 Introduction
In VIP, boundary flux is handled by including source/sink terms in the
interblock flow equations for edge blocks. The outer boundaries of the
grid are normally treated as sealing barriers to flow. Two specific
boundary flux cases are considered: 1) use of an aquifer influence function
to represent a surrounding body of water, and 2) inclusion of flux from a
course grid simulation run as the boundary conditions for a fine grid
model.
3. Influx data should not refer to zero pore volume blocks. Such a
reference is a fatal error for data entry when using the VALUE option.
If the WINDOW option is used in the influx description, the zero pore
volume blocks are ignored, depending on the sinf option used.
4. To use the XCALC, YCALC, and ZCALC options, the DX, DY, DZ,
XCORN...or CORP, KX, KY, and KZ arrays must be entered in the
array data (Section 4.2).
TD PD
td pd
(repeat as required)
Definitions:
INFLUX Indicates that the data being read are influx data.
2
2π φ c t hr e s
binf = ------------------------------
α1
where
α2 k
tc = --------------------2-
ct φ µ d
where
The titles on the fifth card must appear as shown. The optional fifth
keyword controls possible connections to non-aquifer gridblocks
The number of data cards following the fifth card must equal nbinf,
unless an "ENDAQ" card is given.
i1 i2 j1 j2 k1 k2 sinf 00
NOTE: If TD, PD data are not supplied, an infinite radial aquifer will be
assumed, and the corresponding tables for an infinite radial aquifer
will be used by default.
The first data card must have both td and pd equal to zero. Enough
values should be entered to extend to dimensionless times beyond the
time to be simulated or until a linear extrapolation of the last two
values will provide satisfactory results. 00
td Dimensionless time.
(td = tc * t)
Examples: 00
INFLUX 1
WTR CT
...
IINF JINF KINF SINF
WINDOW
1 10 1 15 6 6 ZCALC 00
9 10 14 15 6 6 *2
10 10 15 15 6 6 *2 00
1 2 1 2 6 6 =0
VALUE
1 3 6 =0
3 1 6 =0
ENDAQ
TD PD
... 00
INFLUX 2
WTR CT
... 00
INFLUX 1
WTR CT
NBINF BINF TC
(enter values) 00
1 10 01 15 6 6 ZCALC
9 10 14 15 6 6 *2
10 10 15 15 6 6 *2
1 02 01 02 6 6 =0
VALUE
1 3 6 =0
3 01 6 =0
ENDAQ
TD PD
(enter td pd values) 00
Definitions:
INFLUX Indicates that the data being read are influx data.
7.08 ⋅ k ⋅ h
pi = -------------------------------------, for radial flow
re
µ ln ---- – α 3
r o
α4 ⋅ k ⋅ b ⋅ h
pi = ----------------------------, for linear flow
µd
where
where
The titles on the fifth card must appear as shown. The optional fifth
keyword controls possible connections to non-aquifer gridblocks
The number of data cards following the fifth card must equal nbinf,
unless an "ENDAQ" card is given.
i1 i2 j1 j2 k1 k2 sinf 00
The combination of these two modes allows the user great flexibility in
developing a reservoir study. An initial program run in OUTPUT mode
for a large, coarsely gridded reservoir can provide boundary flux
information for a subset of the reservoir that is to be studied in greater
detail. The subsequent run of the finely gridded portion of the reservoir
uses the boundary flux information to include the effects of gridblocks in
the reservoir that are outside the area of detailed interest. The program
may also be used effectively in OUTPUT mode to determine the direction
and type of fluid flow across any gridblock boundary in the reservoir.
FLUX nflux
INPUT
(MOBWT)
ADJUST
FINE
FLXOIL
(VEOFF)
(CFXOFF)
COARSE
NX NY NZ
nxc nyc nzc
I1 I2 J1 J2 K1 K2
Output Mode: 00
FLUX nflux
OUTPUT
(CUMFLUX)
(VEDIST)
CMPFLX
ON
OFF
i1 i2 j1 j2 k1 j2
Definitions: 00
Alpha label indicating the phase to which each of these options applies
The next six cards define the correspondence of fine gridblocks to coarse
gridblocks. Each data card contains an alpha label followed by an
appropriate number of values, depending upon the direction indicated by
the label. If necessary, more than one card may be used to enter the values,
but the label should not be repeated on subsequent cards.
i1 < I < i2
j1 < J < j2
k1 < K < k2. 00
The data given above are used to define a rectangular subregion of the
reservoir for which boundary flux calculations will be performed in the
simulation module. This data must appear in the order shown.
1 2 3 4 5 6
2
1 2 3 4 5 6 7 8 9 10
1
3
2
3
4 4
5
6
7
5 8
9
Fine Grid
6
Coarse Grid
Figure 10-1: Schematic Representation of Boundary Flux Feature
Examples: 00
The data given below completely describes the coarse and fine grid
models, assuming that boundary flux will be allocated based upon
fluid mobility. 00
FLUX 1
OUTPUT
CUMFLUX
VEDIST
CMPFLX ON
2 5 3 5 1 1 00
FLUX 1
INPUT
MOBWT
ADJUST FINE
COARSE
NX NY NZ
6 8 1
l1 l2 J1 J2 K1 K2
2 5 3 5 1 1
NOVEAD O
2 5 3 5 1 1
VEAREA G
2 4 4 4 1 1
2 2 5 5 1 1
VEONLY W
3 5 5 5 1 1
SEGAREA O
3 5 5 5 1 1
FINE
I1 1 4 6 9
I2 3 5 8 10
J1 1 3 6
J2 2 5 9
K1 1
K2 1 00
11
00000Local Grid Refinement
11.1 Introduction
The LGR option is compatible with rectangular, radial and corner-point
coarse (ROOT) grids. Once a ROOT grid has been defined, it can be
refined and each refined grid can be further refined, with the exception
that if the refinement is a radial grid, then the radial grid cannot be further
refined. There is no limit to the number of grids or levels of refinements.
LGR (name)
ENDLGR
Definitions: 00
It is strongly recommended that all grid refinements are continued through all
communicating layers in the grid system, i.e., normally from layer 1 through NZ.
CARTREF name
i1 i2 j1 j2 k1 k2
nx1 nx2 nx3 .........nx(i2-i1+1)
ny1 ny2 ny3 .........ny(j2-j1+1)
nz1 nz2 nz3 .........nz(k2-k1+1)
(OMIT io1 io2 jo1 jo2 ko1 ko2)
(INGRID in1 in2 jn1 jn2 kn1 kn2)
(CARTREF name) (RADXREF name) (RADYREF name) (RADZREF name)
nested refinement
(ENDREF)
ENDREF
Definitions:
i1, i2, j1, j2, k1, k2 Indices defining the portion of the coarse grid to be
refined.
io1, io2, jo1, . . . Indices defining the portion of the coarse grid to be
removed from the refinement.
i1 ≤ io1, io2 ≤ i2
j1 ≤ jo1, jo2 ≤ j2
k1 ≤ ko1, ko2 ≤ k2
in1, in2, jn1, . . . Indices defining the portion of the coarse grid to be
included in the refinement.
i1 ≤ in1, in2 ≤ i2
j1 ≤ jn1, jn2 ≤ j2
k1 ≤ kn1, kn2 ≤ k2
RADZREF name
i j k1 k2
nr ntheta ri (RMIN rmin)
nz1 nz2 nz3 .........nz(k2-k1+1)
xw1 xw2 xw3 .........xw(k2-k1+2)
yw1 yw2 yw3 .........yw(k2-k1+2)
(OMIT i i j j ko1 ko2)
(INGRID i i j j kn1 kn2)
ENDREF
Definitions:
k1 ≤ ko1, ko2 ≤ k2
k1 ≤ kn1, kn2 ≤ k2
RADYREF name
i j1 j2 k
nr ntheta ri (RMIN rmin)
nz1 nz2 nz3 .........nz(j2-j1+1)
xw1 xw2 xw3 .........xw(j2-j1+2)
zw1 zw2 zw3 .........zw(j2-j1+2)
(OMIT i i jo1 jo2 k k)
(INGRID i i jn1 jn2 k k)
ENDREF
Definitions:
j1 ≤ jo1, jo2 ≤ j2
j1 ≤ jn1, jn2 ≤ j2
RADXREF name
i1 i2 j k
nr ntheta ri (RMIN rmin)
nz1 nz2 nz3 .........nz(i2-i1+1)
yw1 yw2 yw3 .........yw(i2-i1+2)
zw1 zw2 zw3 .........zw(i2-i1+2)
(OMIT io1 io2 j j k k)
(INGRID in1 in2 j j k k)
ENDREF
Definitions:
i1 ≤ io1, io2 ≤ i2
i1 ≤ in1, in2 ≤ i2
Example 1:
NX NY NZ NCOMP
5 5 1 2
LGR BASEGRID
CARTREF REF1
2 4 2 4 1 1
2 3 2
2 2 2
1
RADZREF RAD
3 3 1 1
3 4 .25
1
2*0.5
2*0.5
ENDREF
ENDREF
ENDLGR
Example 2:
NX NY NZ NCOMP OR NX NY NZ NCOMP
5 5 1 2 5 5 1 2
LGR BASEGRID LGR BASEGRID
CARTREF REF1 CARTREF REF1
2 4 2 4 1 1 2 4 2 4 1 1
2 3 2 2 3 2
2 2 2 2 2 2
1 1
OMIT 4 4 3 3 1 1 OMIT 3 4 3 4 1 1
OMIT 3 4 4 4 1 1 INGRID 3 3 3 3 1 1
ENDREF ENDREF
ENDLGR ENDLGR
Example 3:
NX NY NZ NCOMP
5 5 1 2
LGR BASEGRID
CARTREF REF1
2 4 2 4 1 1
2 3 2
2 2 2
1
OMIT 2 4 2 4 1 1
INGRID 3 4 2 2 1 1
INGRID 2 3 3 3 1 1
INGRID 2 2 4 4 1 1
ENDREF
ENDLGR
Example 4:
NX NY NZ NCOMP
5 5 1 2
LGR BASEGRID
CARTREF REF1
1 2 1 2 1 1
2*3
2*3
1
ENDREF
CARTREF REF2
4 5 4 5 1 1
2*3
2*3
1
ENDREF
ENDLGR
Example 5:
NX NY NZ NCOMP
5 5 1 2
LGR BASEGRID
CARTREF REF1
2 4 2 4 1 1
2 3 2
3*2
1
CARTREF REF2
2 3 3 4 1 1
3 2
4 2
1
ENDREF
ENDREF
ENDLGR
Arrays (name)
Definitions:
Example:
ARRAYS
POR CON
0.3
.
.
.
ARRAYS CHILD1
POR NONE
VMOD
2 4 2 3 1 1 EQ
0.25 0.26 2*0.27 0.29 0.3
CORP (EIGHT)
Enter data values as required.
Definitions:
5 6
8 7
MODX
MODY
MODZ
i1 i2 j1 j2 k1 k2 #v
Definitions: 00
i1 ≤ I ≤ i2
j1 ≤ J ≤ j2
k1 ≤ K ≤ k2
+ add
- subtract
/ divide
* multiply
= equal
Example: CORP
INCLUDE grid.inc
MODZ
1 10 1 10 1 3 *3.048
RIGHTHANDED
Definitions:
5 6
8 7
Definitions:
NEWTRAN
Definitions:
BLOCKTR
Definitions:
RMIN rmin
Definitions:
[STD ]
MULT array [NONSTD] [MINUS] [operator]
[ALL ] [PLUS ]
(GRID name)
i1 i2 j1 j2 k1 k2 val
Definitions: 00
FUNCTION ..................
(GRID name)
..
..
..
..
Definitions:
FTRANS
(GRID name1 (name2))
i1 j1 k1 i2 j2 k2 t (tt)
(repeat as necessary)
Definitions:
Definitions:
VOVER array
(GRID name)
i1 i2 j1 j2 k1 k2 (op)
values as necessary
Definitions:
The TOVER cards must follow the array data and precede any FTRANS,
OVER and VOVER cards.
TOVER array
(GRID name)
i1 i2 j1 j2 k1 k2 (op)
Definitions:
Definitions:
Definitions:
12
00000Tracer Option1
12.1 Introduction
This is a short description of a set of tools which have been developed to
improve and expand the simulation capabilities in the analysis and
interpretation of tracer tests and in the design and performance analysis of
waterflood projects. The tools are built around the particle tracking
method which allows accurate simulation of tracer flow associated with
convection and physical dispersion. The method is nearly numerical
dispersion free and allows accurate simulation of tracer flow in field scale
simulation. The algorithm is implemented in VIP-EXECUTIVE and allows
simulation of tracer flow within the framework of three-dimensional,
multi-phase, non-steady state reservoir simulation. In addition to accurate
simulation of tracer flow the software allows: (1) tracking of water fronts
in waterflood operations; (2) construction of three-dimensional flow
trajectories and streamlines of velocity field; (3) calculation of the areal
sweep; (4) visualization and animation of tracer flow.
TRACER
The card notifies the simulator that later in the run the tracer option may
be used. As a result, VIP-CORE generates some grid related information
required by the tracer option. This information is passed to the simulation
module through the restart record.
Definition:
13
Heat Loss Data (VIP-THERM)
13.1 Introduction
Heat loss data must appear at the end of the initialization data.
Definitions:
HTLOSS Alpha label indication that the data being read are
heat loss data.
Definitions:
HTLOSS Alpha label indicating that the data being read are
heat loss data.
Two options for definition of heat loss surfaces are available. In the
“Specified Index Option” (Section 13.5.1), each gridblock face through
which heat loss is to occur is specified by giving the gridblock indices
(i,j,k) and a direction flag (except for radial grids for which heat loss is
allowed only in the z-direction). In the “Automatic Index Option” (Section
13.5.2), a single heat capacity and thermal conductivity are given for the
bordering rock along each reservoir boundary (indicated by a direction
flag). The program automatically assigns heat loss to the outermost active
(non-zero total volume) gridblock face on the specified boundary.
Definitions:
Definitions:
14
00000Parallel Computing
14.1.1 Introduction
Any cartesian grid defined with the LGR option can be automatically
decomposed into subdomains for the purpose of running on a network of
CPU’s (parallel computing).
DECOMP
gridname1 (KEEP) ndx ndy ndz
(X npx1 npx2 npx3 ... npxndx)
(Y npy1 npy2 npy3 ... npyndy)
(Z npz1 npz2 npz3 ... npzndz)
(gridname2 (KEEP) ndx ndy ndz)
.
.
.
ENDDEC
Definitions:
NOTE: If the KEEP option is used then the decomposed grid itself will be
the parent of the new grids. Otherwise the parent of the
decomposed grid will also be the parent of the new grids.
ndx
NOTE: ∑ npxi = NXP
i=1
ndy
∑ npyi = NYP
i=1
ndz
∑ npzi = NZP
i=1
or
2
D1 D2
3
7
D3 D4
8
10
I
1 2 3 4 5 6 7 8 ROOT 10
J
1
2
D1 D2
3
7
D3 D4
8
10
7
D4 D5
8
10
A
References
1. Stone, H.L., ’Probability Model for Estimating Three-phase Relative
Permeability’, Trans. SPE of AIME, 249, pp. 214-218 (1970).
8. Frick, T.C. and Taylor, R.W., Petroleum Production Handbook, Vol. II,
Chap. 35, Society of Petroleum Engineering (AIME), Dallas, Texas
(1962).
9. Carter, R.D. and Tracy, G.W., ’An Improved Method for Calculating
Water Influx’, Trans. AIME, Vol. 219, pp. 415-417 (1960).
12. Jain, A.K., ’An Accurate Explicit Equation for Friction Factor’, J.
Hydraulics Div. ASCE, Vol 2, No. Hy5 (May, 1976).
14. Wallis, J.R., Kendall, R.P., and Little, T.E., ’Constrained Residual
Acceleration of Conjugate Residual Methods’, SPE paper 13536
presented at the Eighth SPE Symposium on Reservoir Simulation,
Dallas (1985).
17. van Everdingen, A. F. and Hurst, W., ‘The Application of the Laplace
Transformation to Flow Problems in Reservoirs’, Trans., AIME (1949)
186, 305-324.
18. Warren, J.E., and Root, P.J., ‘The Behariam of Naturally Fractured
Reservoirs’, Soc. Pet. Engr. J., pp. 245-255 (September, 1963).
19. Thomas, L.K., Dixon, T.N., and Pierson, R.G., ‘Fractured Reservoir
Simulation’, Soc. Pet. Engr. J., pp. 42-54 (February, 1983).
21. Wallis, J.R., Foster, J.A., and Kendall, R.P., ‘A New Parallel Iterative
Linear Solution Method for Large Scale Reservoir Simulation’, SPE
paper 21209 presented at the Eleventh SPE Symposium on Reservoir
Simulation, Anaheim (1991).
23. Redich, O., and Kwong, J.N.S., ‘On the Thermodynamics of Solutions.
V. - An Equation of State. Fugacities of Gaseous Solutions’, Chemical
Review (1949) Vol. 44, pp.52-63.
26. Penelous, A., Rauzy, E., and Freze, R., ‘A Consistent Correction for
Redlich-Kwong-Soave Volumes’, Fluid Phase Equilibrium (1982) Vol.
8, pp. 7-23.
28. Pedersen, K.S., et al, ‘Vicosity of Crude Oil’, Chem. Eng. Sci., 39, (1984)
pp. 1011-1016.
31. Joffe, J., Schroeder, G.M. and Zudkevich, D., “Vapor-Liquid Equilibria
with the Relich-Kwong Equation of State,” (3) AIChE Jour. (1970) 16,
pp. 496-498.
32. Beattie, C.I., Boberg, T.C., and McNab, G.S., “Reservoir Simulation of
Cyclic Steam Stimulation in the Cold Lake Oil Sands,” SPE 18752
presented at the SPE-AIME, 1989, California Regional Meetings.
33. Passut, C.A., and Danner, R.P., I & EC Process Design and
Development, 11, 543 (1972).
34. Hwang, P.K. and Daubert, T.E., I & EC Process Design and
Development, 13, 193 (1974).
36. Kesler, M.G. and Lee, B.I., Hydrocarbon Processing, 55 (3), 153, 1976.
38. Lohrenz, J., Bray, B.G., and Clark, C.R., “Calculating Viscosities of
Reservoir Fluids from their Compositions,” SPE Reprint Series No. 15,
Phase Behavior, 1981, pp. 224-229.
39. Settari, A. and Ito, Y., “Coupling of a Fracture Mechanics Model and a
Thermal Reservoir Simulator for Tar Sands,” Jounal of Canadian
Petroleum Technolgy, V.31, No.9, November 1992, pp. 20-27.
40. Coats, K.H., George, W.D., Chu, C., and Marcum, B.E., "Three
Dimensional Simulation of Steamflooding," SPE 4500 presented at the
41. Robinson, D.B., Peng D., and Chung, Y., “The Development of the
Peng-Robinson Equation of State and it’s Application to Phase
Equilibrium in a System Containing Methanol,” Fluid Phase Equilibra
(1985), Vol. 24, pp. 25-41.
000000Keyword Index
A
A 4-212
ACENTRIC 4-192
ADJUST FINE 10-488
ALL 3-126, 3-127, 3-131, 3-132, 3-133, 3-140
ALPHA 2-116, 2-117
AMULT 13-526
API 4-147, 4-246
APIRO 4-280
ARRAYS 5-357
AUTO 13-526
B
B 4-212
BDGTAB 4-270
BETAG 4-322
BG 4-239, 4-251, 4-259, 4-264, 4-266, 4-270
BGTAB 4-249
BINF 10-478, 10-484, 12-522
BLACKOIL 2-84
BLOCKTR 11-515
BO 4-239, 4-246, 4-251, 4-259, 4-264, 4-266, 4-270
BODTAB 4-264
BOETAB 4-251
BOF 4-303, 4-306
BOFAC 4-239, 4-246, 4-251, 4-259, 4-280
BOGTAB 4-266
BOOTAB 4-259
BOSEP 4-306
BOSTG 4-303
BOTAB 4-239, 4-246, 4-259, 4-264, 4-266, 4-270
BOTINT 2-104
BPTAB 4-151, 4-152
BTBD 5-395
BWI 2-78, 4-310, 4-311
C
C 1-13
CARTREF 11-498
CDPKRH 2-116
CMP 4-330
CMT 4-312, 4-313
CO2TAB 4-232
COARSEN 8-469
COATS 2-111
COEFFICIENTS 4-331
COIL 4-216, 4-291
COMP 4-192, 4-211, 4-215, 4-219, 4-220, 4-278, 4-300, 4-308, 4-337, 4-340
COMPACT 2-93
COMPONENTS 4-191, 4-278
COMPOSITION 4-224, 4-226
COMPSAT 4-280
CON 1-19
CONTACT 2-114
COPY 1-33
CORCHK 2-110
CORNER 2-108
CORP 11-510
CORTOL 2-108
COURSE 10-488
CPG0 4-218
CPG1 4-218
CPG2 4-218
CPG3 4-218
CPO0 4-218
CPO1 4-218
CPOIL 4-291
CPR 13-525, 13-526
CPR0 5-440
CPRTAB 4-348
CR 2-78, 2-79, 5-392
CREEP 2-118
CREEPB 5-418
CREEPC 5-418
CREEPM 5-419
CRF 5-392
CRINIT 2-105
CROSS 2-82
CSORM 2-86
CT 10-478, 10-484, 12-522
CTEOIL 4-216, 4-291
CUMFLUX 10-490
CW 2-78, 4-310, 4-311
D
DATE 2-53, 1-12
DCPRDT 2-79
DEACTIVATE 2-97
DECOMP 14-529
DELTA 4-333
DENO COMP 4-216
DEPF 5-374
DEPTH 1-28, 4-147, 4-148, 4-149, 4-151, 4-152, 4-153, 4-154, 4-224, 4-226, 4-288, 4-291, 5-368
DGOG 4-321
DIAGONAL 2-71, 2-73
DIFEXPTEST 4-337, 4-340
DIFF 2-112
DIM 2-49
DIP 1-28
DIRCN 13-525, 13-526
DJK 4-196
DJKSEP 4-202, 4-297
DKTDSG 2-79
DLIQ 4-300, 4-308
DOB 4-246
DOR 4-246
DOSTD 4-216, 4-291
DOTAB 4-288
DP 4-239, 4-246, 4-251, 4-259, 4-280
DR 5-358
DRELPM 2-100
DRO 4-280
DRSDT 2-106
DRYGAS 2-106
DTHETA 5-361
DUAL 2-111
DWB 2-78, 2-79, 4-310
DX 5-358
DXB 5-358
DY 5-361
DYB 5-361
DZ 5-363
DZB 5-363
DZBCOR 5-364
DZBNET 5-366
DZCORN 5-364
DZNET 5-365
DZOB 13-524
DZVCOR 5-365
E
EACT 4-351
EG 4-182
EIGHTH 2-71
END 2-53
END2P 2-87
ENDDEC 14-529
ENDEOS 4-205
ENDINC 1-14
ENDKV 4-279
ENDLGR 11-497
ENDREF 11-498, 11-499, 11-501, 11-502
ENDSEP 4-297
ENDVIS 4-228
ENTHALPY COMP 4-218
EOD 4-188
EOG 4-182
EOI 4-188
EOS 4-189
EOSINT 4-327
EOSSEP 4-199
EOW 4-182
EW 4-182
EWD 4-188
EWI 4-188
EXCEPT 3-126, 3-127, 3-131, 3-133, 3-140
F
FAULTS 2-98, 6-446
FINE 10-489
FLASH 2-104
FLOW360 2-98
FLOWS 2-120
FLTRC 6-453
FLTRCF 6-456
FLTTCF 6-456
FLTXC 6-453
FLTXCF 6-456
FLTYCF 6-456
FLUX 10-488, 10-490
FLXOIL 10-488
FORM 2-59, 2-60
FOURTH 2-73
FRZPCG 2-94
FRZPCW 2-93
FTRANF 6-457
FTRANS 6-454, 11-516
G
GASCM 4-335, 4-343
GASFRAC 4-337, 4-340, 4-343
GASMF 4-223
GASPLANT 4-230
GASRM 4-175
GASRMT 4-175
GASTRF 5-416
GASWATER 2-83
GBC 2-94
GIBBS 2-102
GOC 4-151
GOR 4-303, 4-306
GR 4-239, 4-246, 4-249, 4-288, 4-291, 4-303
GRID 11-515, 11-516, 11-517
GWC 4-148
H
HARTRAN 11-514
HCPVTAB 2-82
HOTAB 4-288
HR 4-351
HSTAR 4-199
HTLOSS 13-523
HVAP 4-218
HYDBETA 4-323
HYDFRAC 2-118
HYSWO 4-188
I
I 13-525
ICMT 1-44, 1-45, 5-390, 5-392
ICMTF 5-390, 5-391, 5-393
IEQUIL 4-147, 4-148, 4-149, 4-151, 5-396, 5-397
IFT 2-102
IINF 10-478, 10-484, 12-522
INCLUDE 1-13
INFLUX 10-478, 10-484, 12-522
J
J 13-525
JFUNC 2-92
JINF 10-478, 10-484, 12-522
K
K 13-525
KEEPSG 5-398
KEEPSW 5-398
KEYCMP 4-230
KINF 10-478, 10-484, 12-522
KOHYS 5-442
KR 5-376
KRF 5-376
KRFEFF 5-380
KRG 4-163, 4-169, 4-173, 4-175, 4-180
KRGMAX 4-175
KRGRO 4-182, 4-183, 5-409
KRGROF 5-410
KRGRW 5-409
KRGRWF 5-410
KROCW 4-182, 4-183
KROG 4-163, 4-169, 4-173, 4-180
KROINT 2-89
KROLW 5-408
KROLWF 5-409
KROW 4-158, 4-168, 4-171, 4-179
KRW 4-158, 4-168, 4-171, 4-179
KRWG 4-163, 4-167, 4-170, 4-173
KRWRO 4-182, 4-183, 5-408
KRWROF 5-409
KT 5-441
KTF 5-377
KTFEFF 5-381
KTHETA 5-377
KTR 13-525, 13-526
KTR0 5-441
KTTH0 5-441
KTX0 5-441
KTY0 5-441
KTZ0 5-442
KV 4-280
KVALUES 4-278, 4-300, 4-308
KVCOR 4-219
KVTAB 4-280
KVTBS 4-220
KWHYS 5-442
KX 5-376
KXF 5-376
KXFEFF 5-380
KY 5-377
KYF 5-377
KYFEFF 5-380
KZ 5-378
KZF 5-378
KZFEFF 5-381
L
LAB 2-82
LAYER 1-26
LAYERS 6-451
LBC 4-212, 4-215
LDEST 4-294, 4-300, 4-308
LEAKY 6-452
LFRAC 4-294, 4-300, 4-308
LGR 11-497
LINE 2-108
LINFAC 12-522
LIST 1-15
LNVAL 1-33
LNXVAR 1-31
LNYVAR 1-32
LNZVAR 1-32
LX 5-412
LY 5-412
LZ 5-412
M
MAP 2-53, 2-60
MAPX 2-59
MAPY 2-59
MAPZ 2-59
MDEPF 5-374
MDEPTH 5-368
METRIC 2-81
MIS 2-116
MISP 4-318, 4-319
MISPTB 4-318, 4-319
MOBILE 2-94
MOBWT 10-488
MOD 1-35, 11-511
MODLYR 1-36
MULCONTEST 4-343
MULT 1-29, 1-41, 11-515
MULTFL 1-45
MULTIR 1-44
MW 4-192, 4-212, 4-278
MWL 4-308
MWOIL 4-288, 4-291
MWS 4-320
N
NAMES 2-113
NAMESW 2-115
NBINF 10-478, 10-484, 12-522
NCOL 1-16
NCOMP 2-65, 2-66
NCV 2-65, 2-66, 4-219
NETGF 5-366
NETGRS 5-366
NEWTRAN 2-108, 11-514
NFLASH 4-343
NGASFR 4-337, 4-340, 4-343
NINEPT 2-85
NKEY 4-230
NLINES 1-16
NOCHK 2-103
NOCOATS 2-111
NOINIT 2-96
NOLIST 1-15
NONE 3-126, 3-127, 3-131, 3-132, 3-133, 3-140
NONEQ 2-94
NOPRINT 2-111
NOROOT 2-96
NOSKIP 1-15
NOVDB 2-60
NOVDBPACK 2-60
NOVEAD 10-489
NOX 2-110
NOY 2-110
NP 4-288
NPRES 4-337, 4-340, 4-343
NPVUPD 2-95
NR 2-66
NTHETA 2-66
NX 2-65, 10-488
NY 2-65, 10-488
NZ 2-65, 10-488
NZOB 13-524
O
OILCM 4-335, 4-343
OILMF 4-223
OILPROPS 4-291
OILTABLES 4-288
OILTRF 5-416
OMEGAA 4-192
OMEGAB 4-192
OMEGAS 4-296
OMGD 5-417
OMGV 5-417
OMIT 11-498, 11-499, 11-501, 11-502
ORDER 4-351
OUTPUT 4-334, 10-490
OVER 11-516
P
P 4-220, 4-288, 4-303, 4-310, 4-311, 4-312, 4-313, 4-321, 5-397
PATTERN 2-71, 2-73, 2-75, 2-76
PBASEW 4-310
PC 4-192
PCGO 4-163, 4-169, 4-173, 4-180
PCGO1 4-182
PCGO2 4-182
PCGO3 4-182
PCGO4 4-182
PCGOC 4-151
PCGOS 4-163, 4-173
PCGW 4-163, 4-173
PCGWC 4-148
PCHOR 4-192
PCHYSG 2-90
PCHYSW 2-89
PCWO 4-158, 4-168, 4-171, 4-179
PCWO1 4-182
PCWO2 4-182
PCWO3 4-182
PCWO4 4-182
PCWOC 4-147, 4-149, 4-151
PCWOS 4-158, 4-171
PD 10-478, 10-484, 12-522
PF 5-400
PHASE 4-351
PINCHGRID 11-518
PINCHOUT 2-97
PINIT 4-147, 4-148, 4-149, 4-151
PLNTRY 4-230
PMAX 4-329
PMIN 4-330
POLYMER 2-119
POR 2-111, 5-374
PORDEF 4-347
PORF 5-374
PRES 4-337, 4-340, 4-343
PRINT 2-111, 3-131
PRINT ALL 3-125
PRINT ARRAYS 3-126
PRINT COEFS 3-127
PRINT COMP 3-131
PRINT CORNER 3-131
PRINT EQUIL 3-132
PRINT FAULTS 3-132
PRINT INFLUX 3-133
PRINT INIT 3-134
PRINT NONE 3-125
PRINT SEPARATOR 3-139
PRINT TABLES 3-140, 3-141, 3-143
PROPERTIES 4-192, 4-278
PRSTAB 3-141
PS 2-78, 2-79
PSAT 4-147, 4-149, 4-151, 4-152, 4-224, 4-226, 4-239, 4-246, 4-249, 4-251, 4-259, 4-264, 4-266,
4-270, 4-280, 4-337, 4-340, 4-343
PSATF 4-303, 4-306
PSEUDO 2-111
PSTAB 4-153
PV 5-374
PVEXP 2-81
PVF 5-374
PVLINEAR 2-81
PVMULT 4-310, 4-311, 4-312, 4-313
PVTTAB 2-103
PVTTABLE 4-294, 4-300, 4-303, 4-306
PVTW 4-310
PVTWSAL 4-311
R
R 5-358
RADXREF 11-498, 11-502
RADYREF 11-498, 11-501
RADZREF 11-498, 11-499
REACTION 4-351
REGDTM 9-475
REGION 9-473
REGNZ 2-82
REGSEP 9-474
REVERSE 2-93
RI 2-66
RIGHTHANDED 11-513
RMIN 11-515
RPCNO 2-121
RPHYSG 2-91
RPHYSO 2-91
RPHYST 2-92
RS 4-239, 4-246, 4-251, 4-259, 4-264, 4-266, 4-270
RSF 4-288, 4-291
RSM 2-124
RSTAB 4-153
RVAR 1-20
RVTAB 4-154
S
SALINT 4-311
salint 4-311
SATTAB 2-103
SCLFCT 5-394
SDFUNC 4-182
SEBOUND 2-85
SEGAREA 10-489
SEPARATOR 4-294, 4-300, 4-308
SEPTEST 4-303
SG 4-163, 4-167, 4-169, 4-170, 4-173, 4-175, 4-180, 5-397
SGC 4-182
SGF 5-400
SGL 5-405
SGLF 5-406
SGR 4-183, 5-405
SGRF 5-407
SGRM 4-175
SGRO 4-182, 4-183, 5-405
SGROF 5-407
SGRW 5-406
SGRWF 5-408
SGT 4-163, 4-173, 4-180
SGTF 4-169
SGTR 4-173, 4-175
SGU 5-406
SGUF 5-408
SGWT 4-167
SGWTF 4-170
SIDE 12-522
SIGMA 5-413
SIGMAD 5-415
SIGR 4-321
SIGT 4-321
SINF 10-478, 10-484, 12-522
SIXTH 2-76
SKIP 1-15
SLVTAB 4-317
SORG 4-182, 4-183
SORW 4-182, 4-183
SOTR 4-171
STAGE 4-294, 4-300, 4-308
STKZDN 4-201
STONE1 2-87, 2-88
STONE2 2-87, 2-88
SW 4-158, 4-168, 4-171, 4-179, 5-397
SWAPMF 2-111
SWC 4-182, 4-183
SWELLTEST 4-335
SWF 5-400
SWINIT 4-313
SWIR 4-182, 4-183
SWL 5-403
SWLF 5-403
SWMNI 4-158, 4-168, 4-171
SWR 5-403
SWRF 5-404
SWRO 4-182, 4-183, 5-403
SWROF 5-404
SWT 4-158, 4-171, 4-179
SWT nswt 4-171
SWTF 4-168
SWU 5-403
SWUF 5-405
T
T 4-211, 4-215, 4-224, 4-288, 4-291, 4-303
TABLES 4-145
TAMULT 4-310, 4-311, 4-312, 4-313
TB 4-218
TC 4-192, 10-478, 12-522
TCTBD 5-394
TD 10-478, 10-484, 12-522
TDIFFG 5-415
TDIFFO 5-415
TEMP 4-183, 4-220, 4-337, 4-340, 4-343, 5-389
TEMPERATURE 4-329
TENDPT 4-183
TEX 5-413
THCNTR 2-96
THERMAL 2-119
THVAR 1-21
TITLE1 2-52
TITLE2 2-52
TITLE3 2-52
TKWEXP 5-416
TMR 5-421
TMRF 5-424
TMTH 5-422
TMTHF 5-424
TMX 5-420
TMXF 5-423
TMXYL 5-422
TMXYLF 5-425
TMXYR 5-422
TMXYRF 5-425
TMY 5-420
TMYF 5-423
TMZ 5-421
TMZF 5-424
TOLPV 11-519
TOVER 11-517
TR 5-376
TRACER 12-521
TRACK 2-113
TRACKW 2-115
TRES 2-78, 2-79
TRF 5-380
TS 2-78, 2-79
TTAB 4-288, 4-291
TTHETA 5-378
TTHETF 5-381
TTMR 5-427
TTMTH 5-428
TTMX 5-426
TTMY 5-426
TTMZ 5-427
TVMULT 4-310, 4-311, 4-312, 4-313
TWELFTH 2-75
TWOPT 2-84
TX 5-376
TXF 5-380
TXYL 5-379
TXYLF 5-382
TXYR 5-379
TXYRF 5-382
TY 5-377
TYF 5-380
TZF 5-381
V
VAITS 2-95
VALUE 1-23
VDEST 4-294, 4-300, 4-308
VEAREA 10-489
VEDIST 10-490
VEGO 2-100
VEITS 2-101
VELCTY 2-121
VEOFF 10-488
VEONLY 10-489
VEWO 2-99
VFRAC 4-294, 4-300, 4-308
VG 4-239, 4-249, 4-280
VISCOSITY 4-291
VISGAS 4-215
VISGAS PEDERSON 4-213
VISK 4-198, 4-209, 4-213, 4-228
VISKJ 4-209, 4-213, 4-228
VISOIL 4-211, 4-212, 4-291
VISOIL PEDERSON 4-198, 4-209
VISPE 4-228
VMOD 1-39
VO 4-239, 4-246, 4-251, 4-259, 4-264, 4-266, 4-270, 4-280
VOFAC 4-239, 4-246, 4-251, 4-259, 4-280
VOLUME 4-351
VOMIN 4-210, 4-212, 4-288, 4-291
VOTAB 4-288
VOVER 7-467, 11-517
VSHFT 4-192
VSHFTS 4-297
VW 2-78, 4-310, 4-311, 13-523
VWP 4-310, 4-311
W
WATEROIL 2-83
WATIDEAL 2-120
WBSIM 2-67
WIRCT 4-313
WOC 4-149, 4-151
WTR 10-478, 10-484, 12-522
WTRO 4-246
X
X1 4-280
XCORN 5-359
XI 5-399
XIF 5-401
XREG 5-387
XVAR 1-20
Y
YCORN 5-361
YI 5-398
YIF 5-400
YVAR 1-21
Z
Z 4-224
ZBOT 5-371
ZBOTNE 5-371
ZBOTNW 5-372
ZBOTSW 5-372
ZC 4-192
ZCORN 5-369
ZCORNE 5-369
ZCORNW 5-370
ZCORSW 5-370
ZG 4-239, 4-251, 4-259, 4-264, 4-280
ZLNCOR 5-373
ZO 4-280
ZRO 4-280
ZVAR 1-22
000000Subject Index
A
accelerated successive substitution 2-104
aquifers
fluid influx 10-477
number in each gridblock 5-394
areal sweep 12-521
arithmetic operations
applying to a grid 7-460
applying to a grid (VIP-DUAL) 7-463
array data
arithmetic operations on 1-35, 11-511
automatic generation 1-25
constant value 1-19
define new from previous 1-29
directional variation 1-20
directional variation (X,r) 1-20
directional variation (Y,theta) 1-21
directional variation (Z) 1-22
faults 6-448
fracture model 5-442
full array input 1-23
general format 1-17
gridblock depths 1-26
in grid refinements 11-509
modifying 1-34
options for input 1-16
printing by cross-section 2-82
printing coefficients 3-127
printing corner points 3-131
printing of 3-125
printing of input 3-126
replication of depths 1-33
see also grid data
selective replacement 1-39
suppress checking 2-103
transmissibility 1-42
B
Beattie fracture model 5-442
Beattie, Boberg, and McNab 4-346, 4-348
Beggs and Brill correlation 2-67
bibliography A-535
binary interaction coefficient 4-196
for surface separator 4-297
of separators 4-202
of various compounds 4-207
black-oil
laboratory data 4-236
surface separation data with PVT data 4-298
black-oil option 2-84
black-oil reservoirs
how simulated 1-1
boundary flux 10-477
defining regions 10-488
bulk volume tables 2-82
C
capillary gravity equilibrium 2-101
capillary pressure
endpoint scaling 2-87
hold steady 2-93
hysteresis 4-157
Leverett-J calculation 2-92
matrix pseudo 2-111
near-critical 2-102
capillary pressures
how described 4-180
Carter-Tracy method
for estimating acquifer influx 10-478, 10-484
Cartesian grids
how specified 2-65
refinement of 11-498
Chemical Reactions 4-349
CO2 solubility 4-232
Coats, George, Chu, and Marcum 13-524
coefficients
binary interaction 4-196
Coefficients of Interpolation Function 4-331
columns
to be scanned in input 1-16
comments
entering as full lines 1-13
entering inline 1-13
compaction option 2-93
compaction regions 1-44, 1-45, 5-390, 5-391, 5-392, 5-393
compaction tables 4-312
components
naming 4-191
properties data for 4-205
composition
overreads,fracture gas 5-400
overreads,fracture oil 5-401
overreads,gas 5-398
overreads,oils 5-399
variation with area 4-226
variation with depth 4-224
compositional model
definition of 1-3
compositional simulator
overview 1-1
compressibility
creep option (time-dependent) 5-418
for each gridblock 5-392
time dependent 2-118
conductive faults
how specified 6-452
connection transmissibility
modification of 5-419
connection transmissibility modification 1-41
constant equilibrium region
gas 4-223
oil 4-223
coordinate system
handedness 11-513
correlations
D
Darcy vertical flow coefficients 2-67
data
formatting requirements 1-12
how organized 1-7
initialization utility 2-49
repeated values 1-12
sample input 1-9
data checking
of corner-point data 2-109
date
how to specify 2-53
datum depth
for output regions 9-475
printing for regions 3-139
dead oil model
definition of 1-3
deck layout 1-7
defaults
how specified 2-49
density
correlation of 4-201
mixing parameter 5-417
dependency
functional 5-429
depth
arithmetic operations on 1-36
effect on composition 4-224
of gridblocks 1-27
replication of 1-33
descriptive information 2-52
diagonal grids 2-69
diffusivity
of oil and gas phase 2-112
dimensions
how specified for model 2-49
E
echo print
how to activate 1-15
endpoint arrays 11-510
endpoints
temperature variations in 4-183
energy balance equation
how solved 13-523
energy minimization 2-102
enthalpy
how calculated for components 4-217
ideal gas state 4-198
enthalpy coefficients
Passut-Danner 4-208
equation of state interpolation 4-324
equilibration
super-critical 2-105
equilibrium
capillary gravity 2-101
vertical 2-99
vertical equilib. fraction 5-410
vertical equilib. fraction (VIP-DUAL) 5-411
equilibrium data
printing 3-132
equilibrium initialization
of water/gas saturation 4-151
equilibrium regions 5-396
for gas-water option 4-149
for three-phase 4-147
how initialized 4-145
exchange saturation tables
for VIP-DUAL 5-414
external files
how to include 1-14
extra regions
F
fault data
arbitrary gridblock connections 6-453
arbitrary gridblock connections (VIP-DUAL) 6-455
how input 6-448
overview 6-445
printing of 3-132
fault modeling 2-98
fault shift
how applied 6-451
faults
automatic generation of 6-458
conductive 6-452
connection from corner-point data 2-108
five-spot pattern
how specified 2-71
flash calculation 2-104
flow
complete 360 flow 2-98
flow rate arrays
writing to spreadsheet 2-120
flow trajectories 12-521
fluid composition
initial 4-222
fluid properties
constants 2-78
general requirements for 1-6
fluid tracking 2-113, 5-416
flux
at boundaries 10-477
FORTRAN units 1-7
fracture
porosity/pore volume 5-374
fracture blocks
depth 5-373
fracture compaction regions 5-390
fracture compressibility 5-392
fracture modeling 4-345
fracture option 2-118
fracture pressure
overreads 5-399
fracture saturation
overreads 5-399
fracture transmissibility regions 5-393
fracture water induced rock compaction regions 5-391
fractured reservoirs
how modeled 1-2
fractures
hydraulic fracture option 4-323
permeability of 5-379
transfer to, from matrix 4-321
fracturing
effect on porosity 4-346
free field format
explained 1-11
function input option 5-429
functional dependencies 5-429
G
gas composition
constant equilibrium region 4-223
gas hysteresis 4-173
gas phase relative permeability hysteresis 2-91
gas phase viscosity
how modeled 4-213
gas plant
input data for 4-230
gas relative permeability hysteresis 4-173
gas remobilization 4-175
gas-oil capillary pressure hysteresis 2-90
gas-oil drainage 4-322
gas-oil normalized saturations 5-405
gas-oil saturation
for fracture 4-169
gas-oil saturation tables 4-163
gas-water option 2-83, 4-163
geological descriptions
data requirements for 1-6
Gibbs energy minimization 2-102
grid data
arrays of 5-355
see also array data
grid patterns
how setup 2-68
grids
360 flow 2-98
rectangular 2-65
refinement of 11-497
refinement, minimum radius 11-515
rock properties 5-382
root 11-514
saturation-dependent functions 5-402
temperature dependency 5-388
thermal conductivity 5-441
thickness 2-96
thickness ratio 5-366, 5-367
value override 7-464
value override (VIP-DUAL) 7-467
Griffith, Lau, Hon, and Pearson correlation 2-67
H
Hagedorn and Brown correlation 2-67
harmonic integration 2-107
heat capacity
of rock 4-348
variations 5-443
heat loss
data specification 13-525
how calculated 13-523
heterogeneous reservoirs
how modeled 1-2
homogeneous reservoirs
gridding 2-68
hydraulic fracture option 2-118, 4-323
hydrocarbon pore volume table 2-82
hydrocarbons
tracked 2-113
hysteresis
gas 4-173
gas phase relative permeability 2-91
gas relative permeability 4-173
gas-oil capillary pressure 2-90
how specified 4-185
oil 4-171
oil phase relative permeability 2-91
tolerances 2-92
water-oil capillary pressure 2-89
I
ideal gas state enthalpy 4-198
imbibition relative permeability 4-157
imbibition saturation 4-157
immiscible flow
how modeled 1-2
include files
how to specify 1-14
influx data
printing of 3-133
influx option 10-477
initial composition 4-222
initial reservoir temperature
how specified 2-79
initialization
nonequilibrium 2-94
initialization arrays
printing 3-134
initialization data
date 2-53
defaults 2-49
first card 2-49
last card 2-53
title 2-52
what it includes 2-47
initialization module
explanation of 1-1
initialization region reports 2-82
input data
array options 1-16
columns to be scanned 1-16
continued lines 1-16
formatting 1-12
how organized 1-7
repeated values 1-12
skipping lines 1-15
input/output
diagram of 1-8
interations
maximum allowable 2-104
interfacial tension 2-102
K
kPa 1-5
k-values
equilibrium, function of pressure 4-277
example of input 4-286
how calculated 4-219
L
laboratory data
for black-oil 4-236
laboratory units 2-82
Leverett J-function 2-92
line continuation
in input data 1-16
lines per page
how to specify 1-16
local grid refinement (LGR) 11-497
Lohrenz-Bray-Clark correlation 4-209, 4-215
M
map file
how written 2-53
matrix fracture exchange transmissibility 5-412
matrix pseudo capillary pressure 2-111
matrix-fracture diffusion
for VIP-DUAL 5-415
matrix-fracture transfer 4-321
metric units 2-81
miscibility
transition zones 2-117
miscibility pressure 2-117
miscible data
pressure table 4-318
solvent molecular weight 4-320
solvent PVT properties 4-317
miscible flood data 4-315
miscible option 2-116
mobility weighting 2-84
mole fraction
weighting 2-84
multiple reservoirs
in same model 4-145
multipliers
for transmissibility 5-419
N
named fault/region transmissibility multiplier 1-45
near-critical fluid properties 2-102
Newton-Raphson method 2-104
nine-point finite difference approximation 2-85
nine-spot pattern
how specified 2-71
non-Darcy flow behavior 4-323
nonequilibrium initialization 2-94
numbers
how entered 1-12
O
off-band connections 2-97
oil composition
constant equilibrium region 4-223
oil density
how computed 4-215
oil hysteresis 4-171
oil phase relative permeability hysteresis 2-91
oil phase viscosity
how modeled 4-209
oil-vapor phase behavior 2-119
OMEGA data
for surface separation 4-296
Orkiszewski correlation 2-67
output regions
assign separator batteries to 9-474
assigning names to 9-473
for multiple groups of output data 5-385
for VIP-DUAL 5-388
specify datum depth 9-475
overreads
fracture gas composition 5-400
fracture oil composition 5-401
fracture pressure and saturation 5-399
gas composition 5-398
P
parallel computing 14-529
Passut-Danner ideal gas state enthalpy coefficients 4-208
patterns
in grids 2-68
Pedersen viscosity correlation 4-228
Pederson correlation 4-209
Peng-Robinson equation of state 4-188, 4-199
permeability
normalized relative endpoints 5-408
of fractures 5-379
relative 1-29
X/R direction 5-376
Y/theta direction 5-377
Z direction 5-378
permeability behavior
how modelled 4-348
phase equilibrium
calculation of 2-102
phase stability test 2-102
pinchouts
automatic detection of 6-458
nonstandard connections for 2-97
warning messages 3-132
polymer flooded reservoirs
how modeled 1-3
polymer injection 2-119
pore volume
modification options 7-459
pore volume calculation 2-107
pore volume table
for hydrocarbons 2-82
porosity
dual 2-110
porosity deformation model 4-346
pressure
overreads 5-397
pressure table
for miscible data 4-318
printing
R
radial grids 2-65
flow in 2-98
how specified 2-66
refinement of 11-499
Redlich-Kwong equation of state 4-188, 4-199
refinement
of grids 11-497
region names
printing of 3-139
region reports
for initialization 2-82
regions
see output regions, extra regions
relative permeability
from two-phase saturation table 4-155
hysteresis 4-157
hysteresis tolerance 2-92
hysteresis, gas phase 2-91
hysteresis, oil phase 2-91
near-critical 2-102
of water, gas-dependent 4-167
of water, gas-dependent, for fracture 4-170
saturation weighted interpolation 2-89
scaling of endpoints 2-87
two-phase 4-180
using Stone method 2-87
repeated values
how entered 1-12
restart record
defined/how written 1-7
rock characteristics
data requirements for 1-6
rock heat capacity 4-348
rock properties
assignment of 5-382
root grid 11-514
run date 2-53
run title 2-52
S
sand-shale fluid exchange 5-393
saturation
equilibrium initialization 4-151
for fracture (gas-oil) 4-169
for fracture (water-oil) 4-168
gas-oil 4-163
gas-oil, normalized 5-405
initialization 2-95
overreads 5-397
printing of tables 3-139, 3-141, 3-143
water-oil 4-158
skip data
how to specify 1-15
Soave-Redlich-Kwong equation of state 4-188, 4-199
solubility
of CO2 4-232
solvent molecular weight 4-320
specific heat capacity
how calculated 5-440
spreadsheet file
writing to 2-120
Standing-Katz density correlation 4-201
Stone method 2-87
Stone Model 4-155
streamlines of velocity 12-521
summary records
how written 1-7
super-critical equilibration 2-105
surface tension
ratio tables for 4-321
T
temperature
dependency 5-388
initial
how specified 2-79
temperature-dependent endpoints 4-183
thermal conductivity
how calculated 5-441
how specified 5-441
thickness
use in calculating block properties 2-96
Thomas, Dixon, and Pierson 1-2
three-phase reservoir simulator
overview 1-2
time
conversion constant for turbidite models 5-394
title
of run 2-52
Todd and Longstaff option 4-315
tracer tests
analysis of 12-521
tracking
of fluids 2-113, 5-416
of hydrocarbons 2-113
of water 2-115
water types 2-115
transition block oil 2-114
transmissibility
diagonal 5-378
for local grid refinements 11-514
matrix fracture exchange 5-412
modification for connection 5-419
modification options 7-459
multipliers 1-42
override in grid refinement 11-517
use of permeability to calculate 5-376
transmissibility calculation 2-107
Transmissibility Regions 5-392
transmissibility regions 5-392
trapezoidal grids 2-65
trapped hydrocarbon saturation 4-157
turbidite reservoir option 5-393
two phase gas-water option 2-83
two phase water-oil option 2-83
two-phase relative permeability
how described 4-180
two-point scaling 2-87
two-point upstream component mole fraction weighting 2-84
two-point upstream mobility weighting 2-84
U
units
laboratory 2-82
metric 2-81
type allowed/conventions 1-5
utility data 2-49
V
Van Everdingen and Hurst analytical solution
for aquifer influx 10-478
vapor phase
behavior of 2-120
vdb 2-53
vdb file
how written 2-53
velocity dependent relative permeability 2-121
W
Warren and Root theory 1-2, 2-110
water induced rock compaction 4-313
water induced rock compaction regions 5-390
water tracking 2-115
waterflood projects
design and analysis 12-521
water-oil capillary pressure hysteresis 2-89
water-oil hysteresis
how specified 4-185
water-oil normalized saturations 5-403
water-oil option 2-83
water-oil saturation 4-158
for fracture 4-168
wellbore simulation 2-65, 2-67
wells
general modeling requirements 1-5
Z
Zudkevitch- Joffe-Redlich-Kwong equation of state 4-188
Zudkevitch-Joffe equation of state 4-199