pdfcoffee.com_eclipse-pvti-refernce-manual-pdf-free
pdfcoffee.com_eclipse-pvti-refernce-manual-pdf-free
Reference Manual
2008.1
Proprietary Notice
Copyright © Schlumberger. All rights reserved.
No part of this document may be reproduced, stored in a retrieval system, or translated in any form or by any means, electronic or
mechanical, including photocopying and recording, without the prior written permission of Schlumberger.
Use of this product is governed by the License Agreement. Use of this product is governed by the License Agreement. Schlumberger.
Schlumberger makes no warranties, express, implied, or statutory, with respect to the product described herein and disclaims without
limitation any warranties of merchantability or fitness for a particular purpose.
Patent information
Schlumberger ECLIPSE reservoir simulation software is protected by US Patents 6,018,497, 6,078,869 and 6,106,561, and UK Patents
GB 2,326,747 B and GB 2,336,008 B. Patents pending. Schlumberger FrontSim reservoir simulation software is protected by US Patent
2004/0015295A1.
The Calculator, Charisma, ConPac, ECLIPSE 100, ECLIPSE 200, ECLIPSE 300, ECLIPSE 500, ECLIPSE Office, EDIT, Extract, Fill,
Finder, FloGeo, FloGrid, FloViz, FrontSim, GeoFrame, GRAF, GRID, GridSim, Nodal, NWM, Open-ECLIPSE, PetraGrid, PIPESIM,
PIPESIM FPT, PIPESIM GOAL, PlanOpt, Prodman, Pseudo, PVTi, RTView, SCAL, Schedule, SimOpt, VFPi, Weltest 200.
Trademark information
Silicon Graphics and IRIX are registered trademarks of Silicon Graphics, Inc. OpenGL® and the oval logo are trademarks or registered
trademarks of Silicon Graphics, Inc. in the United States and/or other countries worldwide. OpenInventor and WebSpace are trademarks
of Silicon Graphics, Inc. IBM, AIX and LoadLeveler are registered trademarks of International Business Machines Corporation. Sun,
SPARC, Solaris, Ultra and UltraSPARC are trademarks or registered trademarks of Sun Microsystems, Inc. Macintosh is a registered
trademark of Apple Computer, Inc. UNIX is a registered trademark of UNIX System Laboratories. Motif is a registered trademark of the
Open Software Foundation, Inc. The X Window System and X11 are registered trademarks of the Massachusetts Institute of Technology.
PostScript and Encapsulated PostScript are registered trademarks of Adobe Systems, Inc. OpenWorks and VIP are registered trademarks
of Landmark Graphics Corporation. Lotus, 1-2-3 and Symphony are registered trademarks of Lotus Development Corporation. Microsoft,
Windows, Windows NT, Windows 95, Windows 98, Windows 2000, Windows XP, Internet Explorer, Intellimouse and PowerPoint are
either registered trademarks or trademarks of Microsoft Corporation in the United States and/or other countries. Netscape is a registered
trademark of Netscape Communications Corporation. AVS is a registered trademark of AVS Inc. ZEH is a registered trademark of ZEH
Graphics Systems. Ghostscript and GSview are Copyright of Aladdin Enterprises, CA. GNU Ghostscript is Copyright of the Free Software
Foundation, Inc. Linux is Copyright of the Free Software Foundation, Inc. IRAP is Copyright of Roxar Technologies. LSF is a registered
trademark of Platform Computing Corporation, Canada. VISAGE is a registered trademark of VIPS Ltd. Cosmo is a trademark and
PLATINUM technology is a registered trademark of PLATINUM technology, inc. PEBI is a trademark of Veritas DGC Inc./HOT Engineering
GmbH. Stratamodel is a trademark of Landmark Graphics Corporation. GLOBEtrotter, FLEXlm and SAMreport are registered trademarks
of GLOBEtrotter Software, Inc. CrystalEyes is a trademark of StereoGraphics Corporation. Tektronix is a registered trade mark of
Tektronix, Inc. GOCAD and JACTA are trademarks of T-Surf. Myrinet is a trade name of Myricom, Inc. This product may include software
developed by the Apache Software Foundation (http://www.apache.org). Copyright (c) 1999-2001 The Apache Software Foundation. All
rights reserved. MPI/Pro is a registered trademark of MPI Software Technology, Inc. The TGS logo is a trademark of TGS, Inc. POSC, the
POSC logo and Epicentre are registered trademarks of Petrotechnical Open Standards Consortium, Inc. Red Hat is a registered
trademeak of Red Hat, Inc. This product may include software developed using LAPACK (http://www.netlib.org/lapack/), which is copyright
of its authors. Scali is a trademark of Scali Inc.
Table of Contents
List of Figures ..... ...................................................................................................................................................................7
List of Tables ...... ...................................................................................................................................................................8
Batch Mode
For the 2004A version of PVTi the batch mode has undergone a significant revamp. Over the
last few years the user interface of PVTi has evolved rapidly and the existing batch mode facility
no longer adequately supports more recent functionality. There have been 3 significant
modifications to the PVTi batch mode:
1 The way a batch mode is executed has changed. The new way to launch a batch mode run
on a PC is to use the command $pvti -batch filename where filename is the
name of your PVTi project. See "General information" on page 152 for more details on
running batch mode with other platforms.
Panel Addition/Modification
1 There is a new LBC Viscosity Coefficients panel available under the menu option Edit |
Fluid Model | LBC Viscosity Coefficients.... This panel shows, if using the Lorentz-
Bray-Clark (LBC) viscosity model, the current values of the five coefficients.
2 The Thermal Properties panel has been modified to include the new properties introduced
as a result of the new ECLIPSE Thermal export functionality. Specifically, the properties
Heat Cap. A and Heat Cap. B have been changed to Heat Cap. A/G and Heat Cap. B/H.
There is also a new property called Heat of Vaporization which corresponds to the
HEATVAPS keyword exported for ECLIPSE Thermal.
Component Constraints
1 It is now possible to read in projects where fluid samples can have up to 100 components.
However, no functionality involving the EoS flash can be used until a group operation has
been performed in order to reduce the number of components in all samples to <=50.
2 Two new components are available in the PVTi library - Benzene (C6H6) and Toluene
(C7H8). The short names for use in the fundamentals or components panel are BEN for
Benzene and TOL for Toluene. See "Library" on page 95 to view the components in the
PVTi library. To view the library within PVTi itself use the View | Library option.
Sections
There is a new section called APITRACK. This is used when preparing batch files if the API
Tracking export facility is to be made use of in batch mode. It is essentially the same as the
BLACKOIL section but has an extra keyword called SAMPLES, which records the samples for
which PVT tables will be exported.
Manual
1 The section "The Most Asked Questions About PVTi" on page 15 has been updated with
questions regarding the Batch Mode, ECLIPSE Thermal Export and API Tracking
functionality.
2 A number of the tutorials have been amended - in particular the tutorial entitled “Using the
ECLIPSE Thermal Support Module” has been replaced by a tutorial called "Exporting an
ECLIPSE Thermal model" on page 73.
Introduction
This section has been designed as a reference section so that you can quickly access information
about common problems encountered with PVTi without having to spend time looking through
the manual for the relevant section. The questions in this section have been constructed using
the most common support issues and also the InTouch database. Cross-references are provided
where necessary so that readers can access the appropriate parts of the manual for more detailed
information on a particular topic if required.
The questions are:
• "What is PVTi used for? Why do we need it?" on page 16
• "Where do I start? How do I set up a project within PVTi?" on page 16
• "How do I create an experiment along with a series of observations?" on page 17
• "What are the data limitations in PVTi?" on page 18
• "What is the Fluid Properties Estimation facility in PVTi?" on page 19
• "How do I perform regression on multiple fluid samples?" on page 19
• "What regression parameters should I choose?" on page 21
• "What is the difference between normal regression, special regression and automatic (PVTi
selects) regression?" on page 20
• "How does PVTi support gas condensate simulation?" on page 22
• "Can Black Oil tables be extended above the liquid bubble point in PVTi?" on page 22
• "Can PVTi be used if you know the composition of a fluid but do not have any
observations? And vice-versa?" on page 22
• "What black oil correlations are available in PVTi?" on page 23
• "How do I generate the asphaltene phase envelope using PVTi?" on page 23
• "How does PVTi support ECLIPSE Thermal?" on page 24
• "How do I Use PVTi’s Batch Mode?" on page 25
Warning If you try to perform operations on a fluid with incorrectly normalized mole/
weight fractions then the operation may fail.
Note Just because a component is defined in the ZI sample it does not mean there has to be
any of it there. It only has to be defined in the ZI sample to be used in other fluid
samples. If the mole fraction of a C8 component in the ZI sample is set to be 0.0 then
the C8 component can then be use in the OIL sample and the mole fractions set as
required.
Once at least one fluid sample (the ZI sample) has been defined then any experiment supported
within PVTi can be simulated as well as operations such as phase plots, fingerprint plots and
splitting.
For more information on creating fluid samples see "Defining Samples" on page 107.
Creating Observations
If there are no observations at all for a particular experiment then to create one you need to select
the Edit | Observations... option and the Observations panel appears.
On the Experiments column on the left-hand side there is a list of all the experiments that are
available within PVTi and * symbols are next to the ones that you currently have defined within
your project.
Note Currently defined observations for an experiment can be edited in the Observations
folder on the experiment Entry panel.
For more information on creating and editing experiments/observations see "Simulation using
PVTi" on page 117 and/or the tutorial "Simulating experiments" on page 43.
Note *When a splitting operation was performed it was possible to have more than 50
components (up to 100 in fact) but the components had to be grouped back so that there
were less than 50 before any experiment simulation could take place.
2003A
These pre-2003A data constraints have been present in PVTi for 4 to 5 years and, in-line with
the huge increase in computing power in the last few years, we have decided to enhance the data
constraint capability of PVTi so that the following is now available:
• 100 fluid samples
• 100 components per fluid sample* (see below)
• 100 experiments per fluid sample
• 300 observations per experiment
Note It is now possible to read in, save, split and group with fluids containing up to 100
components. However, the limit is still 50 components for any functionality involving
the EoS flash.
Note If two of your fluid samples are not relevant to each other for example they come from
different wells/reservoirs then a separate project should be created for each one.
Disabling Experiments/Observations
You can prevent PVTi from including an experiment in the regression by right-clicking on the
experiment and selecting Don’t use in Regression. A cross appears on the experiment
indicating it is not currently available within the regression facility.
You can disable an observation so that it is not used within the regression by again right-clicking
and selecting Don’t use in Regression. Alternatively, by right-clicking and selecting Set
Weight and then entering zero the observation is also not included in the regression.
Regression Weights
In general there will be a set of values in an observation. For example, if we have a differential
liberation (DL) experiment defined then a viscosity observation would have a value for each
pressure. We have two types of weight: there are single weights for each value of an observation
and global weights that apply to every value in an observation.
By right clicking on an experiment observation the global weight can be set. As mentioned
above, by setting this to zero none of the values in the observation would be used. Alternatively,
you may want to set a global weight for an experimental observation particularly high, for
example, matching the bubble point of a fluid is normally very important if one wants to ensure
that it is a single-phase liquid at the temperature and pressure of the reservoir.
Or maybe you do not trust the accuracy of a particular observation value, for example an oil
formation volume factor (FVF) value in a DL experiment.You may then not want to use a global
weight as all the other observation values look ok. In this case setting a single-value weight to
a very low value helps you match all the other values in the observation during regression as the
rogue, inaccurate value no longer inhibits convergence.
Both the single-value and global weights for an experimental observation can be set in the
Observations panel by selecting the Edit | Observations... option, highlighting the
appropriate observation and then simply typing in your chosen weights.
For a good example of how to use the regression facility, see the tutorial "Fitting an equation of
state to experimental results" on page 50.
• Pc, Tc and ω of any component with mole weight of C7 or heavier (as these are effectively
mixtures of different molecule types and so may differ from library values).
• Ω A and Ω B of any component with mole weight C7 or heavier. Again because these are
mixtures.
• No binary inter-active coefficients because of the risk of over-fitting.
• No viscosity-specific parameters, again because of the risk of over-fitting.
The variables mentioned above are all normal regression variables. The following set of special
regression variables can also often prove useful to get a match between samples:
• Do a multi-feed split to split the plus fraction into 2 or 3 pseudo-components. The α ,
SCTMW and Kw parameters, which control the splitting are then very good choices.
• The Cheuh-Prausnitz A binary parameter if using Cheuch-Prausnitz binaries.
• The mole weight of a plus fraction (if no split on the plus fraction has been performed).
For a detailed explanation of why some of the above are good and bad choices for regression
parameters, see "Regression in PVT analysis" on page 386. There are also more specific
guidelines for choosing regression parameters depending on whether one is dealing with an oil
reservoir, see "Recommended PVT analysis for oil reservoirs" on page 372 - or a gas-
condensate reservoir, see "Recommended PVT analysis for gas condensate reservoirs" on
page 377.
Finally, for a tutorial illustrating the use of the normal and special regression facilities in a
typical workflow see the new "Workflow Tutorial" on page 61.
Note If you have black oil tables exported from before and including 2001A_2 PVTi, or you
have no black oil tables at all then, unfortunately, PVTi will not be able to construct the
compositional model for you. !
Caution You cannot plot theoretical predictions from PVTi of an experimental unless at
least one observation has been defined.
Make sure there is at least one observation defined for this experiment and then right click on
the experiment and select plot. You can super-impose your curve onto a standard phase plot by
using the super-impose button on the toolbar.
Hint There is a shortcut for doing standard phase plots. This is accessed by clicking and
holding the left mouse button down on the appropriate fluid sample icon and then
dragging the cursor into the plot viewer.
There is a more detailed section on the theory behind how PVTi models the presence of wax and
asphaltene phases in "Wax and asphaltene precipitation in PVTi" on page 394.
General information
The PVTi program is an Equation of State based package for generating PVT data from the
laboratory analysis of oil and gas samples. The program may be used through an interactive
menu system or run in a batch mode. An interactive session can be saved as a batch input file,
which contains commands to reproduce the interactive operations. Alternatively, a batch input
file can be run from an interactive session.
Fluid definition
Multiple fluid samples can be defined by specifying components as one of three types. Library
components require only that the appropriate component mnemonic be entered. Characterized
components define properties of plus fractions from a limited set of information. Finally all the
properties of a component can be defined, a facility which can be used selectively to edit the
properties of existing components.
Simulation of experiments
Experiments may be performed on the fluid systems defined using the equation of state model.
Possibilities are:
• saturation pressures
• flash calculations
• constant composition expansions
• constant volume depletions
• differential liberations
• swelling tests
• multi-stage separator simulations.
Other experiments available are:
• composition versus depth
• vaporization test
• multiphase flash
• critical point
• saturation temperature
• first contact miscibility
• multiple contact miscibility (condensing and vaporizing).
• wax appearance temperature
• asphaltene appearance pressure
Hint Almost any result from the experiments mentioned can be used as an observation
against which to regress.
Starting PVTi
Windows platforms
ECLIPSE Program Launcher
1 Start the ECLIPSE Program Launcher.
2 Click on the PVTi button.
3 Select the version and working directories as required.
Command line
1 Type the command $PVTI in a command prompt window.
UNIX platforms
1 Type the command @pvti at the command prompt.
Overview
The tutorials provide a step-by-step introduction to the functionality of PVTi.
Note These tutorials are not intended to teach PVT analysis, but instead concentrate on
illustrating typical work-flows for PVTi.
Each tutorial is divided into a number of distinct sections intended to highlight a specific aspect
of the analysis process. To avoid repetition, later tutorials assume familiarity with some used in
the first tutorials, so it is strongly recommended that you work through them in the order they
are presented.
Available tutorials
The following tutorials are available:
• "Fluid Properties Estimation" on page 34
• "Creating a fluid system" on page 37.
• "Simulating experiments" on page 43.
• "Fitting an equation of state to experimental results" on page 50.
• "Exporting ECLIPSE Black Oil PVT tables" on page 54.
• "Converting a black oil run to compositional" on page 58.
• "Workflow Tutorial" on page 61.
• "Multiphase Flash" on page 69.
• "Exporting an ECLIPSE Thermal model" on page 73.
• "Data analysis and quality control" on page 77.
• "Removing contamination from samples" on page 84.
Introduction
Fluid properties estimation can provide quick-look PVT tables at the well site.
A saturation pressure (bubble or dew-point) together with a reservoir composition are sufficient
inputs to provide a quick-look simulation, giving an initial estimation of fluid properties in
advance of a full fluid analysis in the lab.
After completing this tutorial you should be able to use PVTi as a simulation tool for fluid
properties estimation.
Hint When a new, empty project is created in PVTi, the Fundamentals panel opens
automatically. To re-enter this panel at any later time, select PVTi: Edit |
Fundamentals...
The Fundamentals panel allows you to enter the minimum information required to create
a complete equation of state model.
3 Click on the Enter Weight Fractions check box.
Hint The mole fractions that you see in lab reports are derived from weight fractions and the
mole weights of the components. It is weight fractions that are actually measured.
PVTi allows a choice of either form in the Fundamentals panel.
Hint Only enter mole weights for components whose properties will be characterized, the
other components come from the library. Also, specific gravity is an optional
additional parameter that can be used in the characterization, if it is not specified it is
calculated using a correlation.
5 Click Apply
PVTi adds the mole fractions and the specific gravity of the plus fraction.
6 Click on the Fluid Properties Estimation (FPE) check box
The temperature and pressure fields are now active.
a Enter a Temperature of 220 F.
This is the temperature of the saturation pressure (dew or bubble point) and the
temperature that is used in the generated depletion experiments.
b Enter a Saturation Pressure of 2800 psia.
Hint PVTi uses this saturation pressure to fit the fluid model. The weight of the plus-fraction
is varied, the weight fraction being constant, until the saturation pressure predicted by
the equation of state matches the entered pressure.
7 Set the Project Units to Field, this sets the units that are used on the plots.
8 Click OK.
This is all the information required to fit the equation of state and to generate the Constant
Composition Expansion (CCE), depletion experiment (either differential liberation or constant
volume depletion) and the optimized separator. After fitting the equation-of-state and creating
the experiments, default plots from the depletion experiments are drawn along with the phase
curve for the fitted fluid.
The methods used in Fluid Properties Estimation are explained in "Fluid Properties Estimation"
on page 384.
9 PVTi: Run | Simulate
This opens a complete report for the project including the results from all the created
experiments.
Hint By clicking on one of the experiments in the sample tree with the right mouse button,
and selecting Report..., you can view the reports for individual experiments
separately.
Flash calculations
1 Right-click on ZI in the project tree-view and select Properties Estimation (FPE)...
a Enter a temperature of 60 F
b Enter a pressure of 15 psia
c Click OK.
2 Right-click on the newly created Flash simulation (FLASH1) and select Report to see the
results of flashing the reservoir fluid to standard conditions.
This allows you to attempt any Flash calculation on the reservoir fluid.
Hint The Properties Estimation panel can also be used to create additional separators,
saturation pressure or depletion experiments, for example at other temperatures.
Discussion
Fluid Properties Estimation is a useful tool, particularly in situations where full lab analysis of
the fluid is not available for a complete equation-of-state matching project.
For full details of the operations performed during fluid property estimation see "Fluid
Properties Estimation" on page 384.
During fluid properties estimation, the project created is a complete PVTi project. This means
that a more experienced user has access to the rich range of features within the product. At the
same time, the less experienced user can use PVTi for Fluid Properties Estimation without
requiring in-depth knowledge of equation-of-state methods or PVT analysis.
Introduction
The PVT report for this fluid contains details of three experiments: a Constant Composition
Expansion experiment, a Differential Liberation experiment, and a Bubble Point experiment.
The later tutorials describe how the experimental results may be used to fit an equation of state
to the experimental behavior, and how this fitted equation of state can be used to generate PVT
tables for use in reservoir simulations.
This tutorial shows how to set up basic fluid properties in PVTi and how to display the phase
envelope for the defined fluid.
Defining a fluid
PVT analysis involves fitting an Equation of State to experimental data and then using this
Equation of State to produce PVT tables for use in reservoir simulations. The first step is to start
up PVTi, and import the component and fluid definitions.
Table 5.2 shows the component and fluid definitions that are used in this tutorial.
1 Start PVTi (if you are unsure about this see "Starting PVTi" on page 31).
2 Select PVTi: File | New...
3 Enter BLACK.PVI as the project name in the file selection window.
4 Click on Open on PC or OK on UNIX platforms.
Note This has defined BLACK as the prefix for any files that are written by PVTi.
The Fundamentals panel opens so that basic project information can be entered.
1 Enter CO2, N2, C1 and C6 into the Components column.
2 Click Apply.
3 Click Yes so that PVTi fills in the library component names.
4 Enter the mole fractions from Table 5.2 and the details for the C7+ component into the
Fundamentals panel and click OK.
Note The components for which no mole weight or specific gravity has been specified are
automatically set to use the PVTi component properties library (see "Component
types" on page 102.)
Hint The component properties can be examined by selecting PVTi: Edit | Fluid Model |
Components.... This panel can also be used to add additional components, the select
alternative characterization methods and to manually defined component properties.
Hint Only one sample, ZI, is defined in the Fundamentals panel. Additional samples can
be created using PVTi: Edit | Samples | Names... and mole fractions can be entered
using PVTi: Edit | Samples | Compositions...
Program options
1 PVTi: Utilities | Program | Options...
This opens the Program Options panel.
2 Set the Separator GOR calculation to use Liquid at Stock Tank Conditions.
3 Set the Temperature-dependence for volume shifts to be calculated by Polynomial
correlations. (See "Shift parameters" on page 321.)
The Program Options panel should now contain the following data:
Field Data
Definition of Liquid Saturation in CCE Sliq=Vliq/Vsat
Treatment of volume shifts Dependent
Separator GOR Calculation Liquid at Stock Tank Conditions
Temperature dependence for volume shifts Polynomial correlations
Specify/Calculate density and molar volume units User units
Specific Heat Capacity Coefficients and Calorific Values
Calculated compositions No Save to samples
Component Library Katz-Firoozabadi
Experimental Compositions Output to Screen/PVP
Experimental Results Always Output to PVP
Plot Vectors No Output to file
Print File Output A4 format
Definition of GOR in Diff Lib Normal
Definition of Oil relative volume in Diff Lib Oil FVF = Voil(p)/Voil(stc)
Black Oil Table Output All Data
Flash Calculations E300 Flash
Sample mole fractions when regressing Keep Fixed
Note A default phase plot with a single quality line can be generated by dragging ZI from
the tree-view of the project (in the left pane of the main window), and dropping it on
to the main plot workspace.
Hint The complete project can be saved using PVTi: File | Save... This, effectively, saves
a history of the project. The original fluid description is saved, along with SPLIT or
GROUP sections for split and group operations you perform. By choosing to save
current modifications, the system is saved in its current state, after all splits, groups,
etc., have been performed. For work in progress it is usually better to use Save so that
past steps can be recovered. For a final fluid model, the Save (Concise) option allows
a complete description of the final model to be saved, without the steps taken to get
there.
This file can now be read in using the PVTi: File | Import | SYSTEM... option.
Introduction
This tutorial describes how experimental observations can be entered into PVTi and how the
experiments can then be simulated from an existing fluid definition.
1 Start PVTi (if you are unsure about this see "Starting PVTi" on page 31).
2 File | Open...
3 Open the file FLUID_DEF.PVI that was created in the last tutorial. (Alternatively open
the supplied tutorial file FLUID_CORRECT.PVI).
Setting units
1 Utilities | Units...
2 Set the Unit Type to Field
3 Set the Temperature Unit Type to Fahrenheit
4 Set Mole Fraction or Percentage to Percentage
5 Set Absolute or Gauge Pressure to Gauge.
6 Click on OK.
The multi-stage separator or SEPS experiments can also be called a separation test in
PVT Reports.
The Experiment Entry window now shows three folders: General, Observations and
Components. These folders are used to define the experiment entry form.
3 Select the Observations folder.
4 Click in the top left cell of the table and select Pressure from the drop-down list in that
cell.
5 In the second column select Relative Vol. from the drop-down list.
Hint By making the column headings the same as those in Table 5.4, the task of data entry
is simplified. The ability to tailor the table means that data entry can then be further
accelerated by importing observations from a text file or the clipboard.
6 Click on Next.
The table now shows two folders. The Components folder has disappeared as there were
no component observations selected; the General folder now shows an entry field to select
fluid type and another to enter the temperature of the experiment.
7 In the General folder, enter the temperature from Table 5.4 (220 F).
8 Select the Observations folder.
The Observations folder now shows a two-column table with the columns selected
previously. The table resembles Table 5.4.
Table 5.4 is provided in the file CCE_TABLE.TXT
9 Right-click in the table and select Table Import | From file...
10 Select CCE_TABLE.TXT and click on Open.
11 In the Text Import Wizard turn on Ignore Records and set the number of records to ignore
to 1 (since we want to ignore the column headings).
The view of the table should no longer contain the first row.
12 Click on OK.
Note The error message “Cannot delete rows from this table” appears This
is because the table has a fixed length and the file we are importing from has fewer
rows than the table. This message can be safely ignored.
13 Click on OK to remove the message “Cannot delete rows from this table”.
The table now contains the same information as Table 5.4. As the experiment information
is complete, the experiment can be created.
14 Click on Next to create the experiment.
Hint The information about which experiments have been defined, and for which
observations have been entered for those experiments, is contained in the Data Tree.
Hint Double-clicking on one of the small plots swaps it with the large plot.
4 Examine each of the plots and note how well (or badly) the simulation has matched the
data.
Hint You can right-click on an axis and select Show Edit Box from the drop-down menu.
This opens the Axis Property Editing panel. In this panel you can customize the axes,
for example by changing the units used in plotting.
Discussion
In this tutorial an existing project was extended to include experiments. Constant Composition
Expansion, Bubble Point and Differential Liberation experiments were imported and simulated
for the defined fluid. The match between the experimental observations and the simulated
results was examined using the plotting facilities in PVTi. The fluid model can then be adjusted
so that it provides the best fit (in a least-squares sense) to the experimental observations (see
"Fitting an equation of state to experimental results" on page 50). The fitted fluid definition is
finally used to generate PVT tables for ECLIPSE (see "Exporting ECLIPSE Black Oil PVT
tables" on page 54).
Introduction
This tutorial illustrates the fitting of the fluid definition to the experimental observations.
The fluid definition and experiments are read in from an existing PVI file and the regression
facilities in PVTi are used to generate an improved fit between the two.
Hint The Data Tree shows the contents of the opened project.
Hint The properties of C1 and C2 are well known and generally do not differ significantly
from the library properties.
Grouping the C3, IC4, NC4, IC5, NC5 and C6 components into a single regression variable
preserves monotonicity between the components, in addition to creating a variable that accounts
for 19.5% of the total composition.
Hint The plus fraction (C7+) contains a mixture of components C7+ and higher, so its
properties may not be so well-determined. This makes the plus fraction an ideal
candidate for regression to fit the equation of state to the experimental results.
Sensitivity analysis
Sensitivity analysis is used to establish which fluid properties most affect the difference
between the observed and simulated values. The sensitivities are calculated for critical
temperature and pressure for each experiment, for both regression variables. Finally the most
sensitive properties will be selected for use in the regression.
Hint In any regression, having a few very sensitive parameters is preferable to having
hundreds of insensitive ones. Always look for parameters that can be discarded (this is
called conditioning the problem - an ill-conditioned problem is difficult to solve).
Hint The second part of the Select EOS parameters for regression panel relates to
binary interaction coefficients.
12 Click on OK.
Hint Right-clicking on an experiment allows you to choose to turn off that experiment
during the regression process.
Discussion
In this tutorial a fluid definition, an Equation of State and some experiments along with their
observations were imported from an existing PVI file (the file created in "Simulating
experiments" on page 43). The definition was adjusted so that it provided the best fit (in a least-
squares sense) to the experimental observations ("Fitting an equation of state to experimental
results" on page 50). This fitted fluid definition can now be used to generate PVT tables for
ECLIPSE (see "Exporting ECLIPSE Black Oil PVT tables" on page 54).
Introduction
Once the fluid description has been fitted to the experimental observations, it may be used in a
reservoir simulation. PVTi facilitates the transition between a fluid description and the PVT
keyword description required by the ECLIPSE family of simulators.
In this tutorial PVT tables are created for the fluid definition developed in the tutorials 2, 3 and
4. The output tables are then used in an ECLIPSE simulation.
1 Start PVTi (if you are unsure about this see "Starting PVTi" on page 31).
2 Open REGRESS_SECTION.PVI created in the last tutorial (alternatively, open
REGRESS_CORRECT.PVI).
Hint PVTi: File | Export Keywords | Oil reservoir... produces an export panel for all
available Differential Liberation experiments.
2 Select PVTO and PVDG (Live oil and dry gas) on the radio button menu.
The Separators drop-down menu becomes active. This is because the produced fluid from
the Differential Liberation experiment must be passed through a surface separator to
calculate, for example, surface gas-oil ratios. The default is a separator at Standard
Conditions. If any separator experiments were defined for this sample, they would also
appear here.
3 Click OK
4 In the File Selection box, enter ECLIPSE100 as the name of the export file.
The keywords are generated and the Display Output module shows the generated keywords.
Note The comments prefixed with --PVTi that appear before each keyword are the concise
version of the current PVTi project. This is the minimum information PVTi requires to
create the tables. This information can be used to rapidly convert an ECLIPSE
BlackOil data-set to an ECLIPSE Compositional data-set.
Caution Avoid editing the --PVTi prefixed comments. Any changes may invalidate the
file, preventing PVTi from reading it.
Note The sample SWELLSAM has been added to the sample tree. This sample is the swelled
sample that was obtained by swelling the original sample with vapor that was split off
just below the bubble-point of the fluid. PVTi automatically swells the sample with the
vapor from the bubble-point so that the table can be extended to values above the
original bubble point.
The information in the keywords is also shown in the main plot space, keyword PVTO is shown
in Figure 5.4.
Note In this case, only RSVD is generated. This is because the reservoir is all initially in the
liquid phase. If there were a gas-oil contact, both RSVD and RVVD would have been
generated. If the reservoir were all in the gas phase, only RVVD would be generated.
Note At this stage you may want to view the keywords or plot them. For details on how to
do this, refer the "ECLIPSE Office User Guide".
7 Close the PVT Section, saving the file with the new keywords.
8 Open the Initialisation Section
9 Initialisation Section: File | Import | Append..., and import the equilibration keyword file
(RSVD.PVO). Click OK to remove the warning message.
10 In the Equilibration Data Specification keyword, set the Rs/Pb v Depth Table to 1, so
that the imported RSVD keyword is used.
11 Close the Initialisation Section, saving the file with the new keywords.
12 Run the simulation.
Discussion
In this tutorial an existing fluid definition was imported into PVTi and exported as PVT tables
that were used in an ECLIPSE Black Oil reservoir simulation. The basic requirements are that
PVTi must simulate a Constant Volume Depletion or Differential Liberation experiment and a
Separator experiment to allow the generation of black oil tables from the fluid definition.
Note Many conversion projects require conversion of SCHEDULE section keywords, which
is outside the scope of this tutorial. This tutorial covers conversion of fluid properties
and equilibration. Specifically, the PROPS and SOLUTION sections.
Note Some familiarity with ECLIPSE Office is assumed for this tutorial.
The data for this project are provided with the standard installation of PVTi under the directory:
$ECL/2007.1/pvti/tutorials
and should be copied to the local directory before starting the tutorial.
Note You should choose a short name for your directory and the name must not contain
spaces. You can use underscore characters. ECLIPSE does not recognize directory
names that are long or that contain spaces.
Introduction
In this tutorial, the black oil PVT tables (PVTO and PVDG) and the Equilibration table (RSVD)
are converted into a full compositional model and composition versus depth table (ZMFVD).
This allows the ECLIPSE data-set from the previous tutorials to be run as a full compositional
case.
Caution The --PVTi comments written out with the keywords are used by PVTi to
reconstruct the original fluid model. Without these there is not enough information
to convert blackoil projects to compositional models. It is important that the lines
prefixed by --PVTi in the ECLIPSE data-set are not edited or moved around
in the file.
1 Start ECLIPSE Office with the project created in "Exporting ECLIPSE Black Oil PVT
tables" on page 54 or create a new project and load in the data-set
ECLIPSE100_FULL.DATA.
2 Select the case.
Note The launch button has a default selection of launching PVTi with the PVT section of
the current case. PVTi reads this PVT section, creating a PVI file from the --PVTi
comments.
Hint The reservoir temperature is the temperature in the Differential Liberation experiment
definition. You can right-click on the DL1 experiment and select Edit... from the drop-
down menu to view the definition of the Differential Liberation experiment.
4 Click OK.
5 Export the fluid model to FLUID.PVO
6 PVTi: File | Exit (There is no need to save the PVI file as it can be created from the
ECLIPSE Office case).
Hint Again, PVTi searches for the --PVTi comments and uses them to construct a PVI
project file.
Discussion
In this tutorial, an ECLIPSE BlackOil simulation data-set was converted to ECLIPSE
Compositional using the integration of PVTi and ECLIPSE Office.
The insertion of the --PVTi comments into the keyword export from PVTi is a powerful tool,
not just for converting data-sets, but also for developing projects in either black oil or
compositional models.
Introduction
This tutorial illustrates a typical workflow for an oil or gas condensate. It involves splitting the
C7+ fraction into 2 pseudocomponents, special regression, normal regressing, grouping
components, and matching viscosity data.
We have an oil PVT case, with two fluid samples ZI and W2 (Well 2). There is a C7+
characterization with CO2 present. Well 2 has C7+ has a different MW and Specific Gravity,
but the C7+ has been characterized with the ZI fraction only at this point. They are going to
inject CO2 into this field, so there is a Swelling Test with CO2.
The files for this tutorial are provided in the default PVTi installation in the following directory:
$ECL/2007.1/pvti/tutorials
and should be copied to the local directory before starting the tutorial.
This tutorial contains the following sections:
• "Comparing the Default EOS Calculations to the Observations" on page 61.
• "Splitting the C7+ Component" on page 62.
• "Special regression to adjust the tail in the splitting calculation" on page 63.
• "Normal regression to fine tune 11 component match" on page 64.
• "Grouping Like Components to Reduce Nc" on page 64.
• "Regressing to match viscosities" on page 67.
• "Discussion" on page 68.
4 Click OK.
The C7+ has been split.
5 Select Edit | Samples | Compositions.
Check the mole fractions of the 2 pseudo heavy components. The split creates FRC2 with
small mole fraction (0.0477). We would rather have a splitting that has more of the mole
fraction in the heaviest component so we will perform another split.
6 Close this project and do not save the project.
7 Open WORKFLOW.PVI again.
8 Plot the Phase Diagram.
9 Select the Split panel and repeat step 6.
10 This time change the Mole Weight of Heaviest Pseudo Component to 300.
This gives us a larger mole fraction for FRC2.
11 Click OK.
12 Check the Samples.
The FRC2 mole fractions are 8.1% and 13.9%
Hint Regressing on the viscosity coefficients after the phase behavior of the fluid has been
matched is always strongly recommended.
We will also increase the weights of the Bubble Point Pressure (Sat. pressure) Observation as
this is a very important criterion to match.
1 Right-click on Vapor Visc and enter 0 as the Set Weight.
2 Right-click on Liquid Visc and enter 0 as the Set Weight.
Note Setting the weight of an observation to zero turns off that observation.
Note You may also view the Sensitivities and the Hessian, Covariance and Correlation
matrices in the Report panel. These matrices are explained in the "Reference section"
on page 89.
Note Do not attempt to get a perfect match since the grouping process described below
changes the match. The key here is to get the Swelling Test; Saturation Pressures to
match better before grouping.
CO2 1
C1 2
C2 3
C3 3
IC4 3
NC4 3
IC5 4
NC5 4
C6 4
FRC1 5
FRC2 6
4 Click Update.
5 Check the component (Group) names. You change them if you do not like a particular
name.
6 Click Update again, and then click OK.
7 Click the Edit sample composition’s button to view the mole fractions of the new group.
The usual rule is if one of the group mole fractions is significantly smaller than the others
group it into one of its neighbor groups.
8 Plot a Phase diagram on top of the un-grouped diagram.
If they are close to each other then the grouping is good. Your phase match should still be
close to the observed data.
9 Click GO and view the bubble point pressures
10 Plot the experiments to see how much the grouping has changed the calculated results.
11 If the match looks reasonable save this characterization with a new file name.
12 If the match has changed substantially, close this project without saving and read in the
.PVI file from the save at the end of the previous section.
13 Re-group with a different selection of groups or number of pseudocomponents.
Second iteration
1 In the Regression panel press Reject to return the characterization back to the pre-
regression values.
2 In the Select EOS parameters for regression panel enter 1 in the first box, 2 in the
second box, 3 through 11 in the remainder of boxes under the heading VcritV.
3 Click OK.
4 In the Regression panel click Run.
5 Plot the Liquid and Vapor viscosity and compare the new match with the observations.
Note this new match is better than the first match.
6 Press Accept to accept these results.
7 Save this characterization with a new file name.
The phase matching process is now complete. You are ready to export the PVT properties or
characterization for ECLIPSE simulations.
Discussion
This tutorial illustrated a typical workflow for an oil or gas condensate. It involved splitting the
C7+ fraction into 2 pseudocomponents, special regression, normal regressing, grouping
components, and matching viscosity data.
Introduction
The multiphase flash experiment tends to find more than two phases in systems with
Asphaltene/Waxes and/or with CO2 rich fluids at low temperatures. This tutorial demonstrates
multiphase flashes with both systems. The files for this tutorial are provided in the default PVTi
installation in the following directory:
$ECL/2007.1/pvti/tutorials
and should be copied to the local directory before starting the tutorial.
The tutorial is split into the following sections:
• "Asphaltene and wax system" on page 69.
• "CO2 Rich Fluids" on page 70.
• "Summary" on page 72.
We will attempt to verify the phases present with a 0.8 fraction of MI with PVTi multiphase
flash.
1 Load MULTIPHASE-CO2.PVI into PVTi.
This contains the fluids and characterization from the SPE paper,[Ref. 1].
2 View the compositions of the fluid sample by selecting Edit |Samples - Compositions.
3 Click OK to close the panel.
Note The MFLASH5 experiment can sometimes label both the phases as liquid. However,
one of them is clearly a vapor as the density value is 5.64323 lb./ft.3.
Summary
This tutorial demonstrated the multiphase flash experiment. It tends to find more than two
phases in systems with Asphaltene/Waxes and/or with CO2 rich fluids at low temperatures.
This tutorial demonstrates multiphase flashes with both systems.
References
Guler B. et al, "Three- and Four-Phase Flow Compositional Simulations of CO2/NGL EOR" [Ref. 1]
SPE 71485,
Introduction
This tutorial demonstrates using the new ECLIPSE Thermal export facility where a file can be
exported containing a fluid model suitable for use in ECLIPSE Thermal. For technical
information on the ELCLIPSE Thermal export facility see "ECLIPSE Thermal Export Module"
on page 401 and for more general workflow guidelines see "Compositional Data for ECLIPSE
Thermal" on page 367. The files for this tutorial are provided in the default PVTi installation in
the following directory:
$ECL/2007.1/pvti/tutorials
and should be copied to the local directory before starting the tutorial.
The tutorial is split into the following sections:
• "Verifying the Validity of the Fluid Model" on page 73.
• "Fitting the Component K-values" on page 74.
• "Viewing the K-value Fits" on page 75.
• "Exporting the Model" on page 76.]
Note Simulations using ECLIPSE Thermal tend to use fluids consisting of two or three
components.
Note The default values here are Pmax=1000psia, Pmin=50psia, Tmax=400F and Tmin=50F
and are considered reasonable max/min parameters within a reservoir. However, every
reservoir is different and any knowledge of these parameters for your particular
reservoir should be entered.
Note The Fit Crookston Coefficients panel enables you to calculate the optimum values of
the coefficients A-E in Crookston’s formula, so that the best fit is found to PVTi’s EoS
predicted K-values for each component over the temperature and pressure range
defined by the user.
Hint You can also access the Superimpose option using the toolbar.
10 Now click Apply on the Fit Crookston Coefficients panel to start the run.
Caution Make sure you turn off the Superimpose option before moving to the next
section.
Hint Plots can either be performed of K-values versus pressure (at constant temperature) or
K-values versus temperature (at constant pressure).
First we will look at the K-value versus temperature fits, which are dictated by the D and E
coefficients (just D in our case).
2 Enter 1250psia in the Enter Constant Pressure box and 400F and 200F as the
Maximum and Minimum Temperatures respectively. Now click Apply.
The PVTi EoS-based K-values are shown by the points and the K-values calculated using
Crookston’s equation are shown by the curves.
3 Experiment by changing the value of the constant pressure in the range 1000<P<1500 (our
chosen pressure range) to see how well Crookston’s formula models the temperature
dependence of the K-values at a given pressure.
In general, for the C20+ component, the D coefficient, models the observations pretty well
over the whole pressure-temperature range, although using the E coefficient as well may
well help slightly
4 Click the Plot at Constant Temperature box and enter 300F.
5 Now enter the appropriate pressure range, which is 1500 psia and 1000 psia for the
maximum and minimum values. Click Apply.
6 Again, experiment by changing the value of the constant temperature in the range
200<T<400 to see how well Crookston’s formula models the pressure dependence of K-
values.
Note If we had decided to export the KVWI keyword (that is use Wilson’s formula to model
K-values) then we would not have needed to fit the coefficients of Crookston’s
equation and could have exported straight away from the Export for ECLIPSE
Thermal panel. Although this may be quicker, Crookston’s formula models K-values
much better and spending this extra time is worthwhile.
Introduction
In addition to allowing an equation of state to be fitted to laboratory results and facilitating the
generation of ECLIPSE BlackOil/ Compositional PVT data, PVTi also provides material
balance checks to assess data quality. For information on the calculations involved in material
balance checking see "Compositional material balance" on page 307.
Note Problems with the observations in a PVT report equate to problems with the fitted fluid
model. It is therefore recommended that material balance checks are carried out on all
PVT data.
In this tutorial an existing project file (GAS.PVI) is read into PVTi and the data are checked
and modified for material balance errors.
1 Start PVTi (if you are unsure about this see "Starting PVTi" on page 31).
2 PVTi: File | Open...
3 Open GAS.PVI
Hint The Data Tree should show that there are two experiments, CVD1, defined with 10
different types of observations, SEPS1 with observations of fluid mole fractions
(liquid and vapor) and CCE1 with observation of Vapor Z-Factor.
4 Click and drag the ZI node from the Data Tree and drop it into the Main Plot Window.
The phase envelope should look like Figure 5.6.
The messages show that some mole fractions were calculated as negative, so there are
clearly problems with the data. PVTi supplies various options for plotting the data to
ascertain the source of the errors. The first type of data check to perform is to view the
pressure variation of the gas compositions.
3 Output Display: File | Close
4 Click Plot in the Material Balance panel.
a Select Vapor Compositions v Pressure in the Select Plot Type panel and click on
Plot.
b Click Close
Many of the components have non-monotonically varying gas compositions. In general, there
are several fluids or analyses available, and bad data can be discarded. However, if no other data
is available PVTi offers tools to make modifications to the bad data.
Note The other components are modified in proportion to their existing mole fractions.
Placing a letter in the thin column to the left of each column of modifiers allows the
proportion of that component to be fixed and thus not modified in proportion to its
existing mole fraction.
The K-values should plot monotonically in that N2 should be the largest, followed by C1, etc.
This is clearly not the case, so although there are now no calculated negative compositions, the
modified fluid definition is not fully consistent.
In this plot, one line is generated for each pressure stage. The Hoffman F coefficients
correspond to C1, C2 etc. and the lowest to C11, C12+. In general, these lines should be
monotonic with pressure, with the highest pressure at the top. This plot shows most of the error
to be in the first stage.
Recovery calculations
PVTi can allow recovery calculations to be performed if a valid Constant Composition
Expansion, Constant Volume Depletion and Separator test exist.
1 Right-click on CCE1 in the project tree-view and select Recovery... on the pop-up menu.
2 Click on Report to perform the recovery calculation.
Hint This assumes that there is no direct production of reservoir liquid. If you want to
include direct production of reservoir liquid, you need a relative permeability table,
which you can enter be clicking on Rel. Perm.
Note PVTi runs the material balance check on the Constant Volume Depletion experiment
selected, and performs recombination on the Separator selected, before performing the
recovery calculation.
Discussion
This tutorial has illustrated how fluids may be examined for consistency and, if necessary,
modified, within a PVTi project.
Introduction
Oil-based muds are in widespread use and often contaminate PVT samples taken at the well-
site.
This tutorial involves the cleaning of a sample that is contaminated by an oil-based mud.
The PVI file CLEAN.PVI is used for this tutorial and is provided in the default PVTi
installation in the following directory:
$ECL/2007.1/pvti/tutorials
and should be copied to the local directory before starting the tutorial.
The tutorial is split into several sections:
• "Introduction" on page 84
• "Removing oil-based mud contamination by skimming" on page 84
• "Removing oil-based mud contamination by subtraction" on page 85
• "Discussion" on page 86
Introduction
The VERSION keyword, introduced in 2000A, allows a systematic method for updating old
PVI files to be compatible with the latest version of PVTi. This tutorial describes how the
keyword can be used to convert an old PVI file into the current version.
Caution Files in FIELD units containing Differential Liberation (DL) experiments that
have Gas formation volume factor (GFVF) observations must be updated to the
current version.
Hint If you are unsure of the form of the VERSION keyword see "VERSION" on page 280.
Note You could now use the file CONVERT.PVI in a normal session and PVTi interprets it
according to the version specified by the VERSION keyword.
Hint You can compare the file NEW.PVI to OLD.PVI to see the differences (the DL
observation GFVF is converted from rb/stb to rb/Mscf and the heat capacity keywords
are added in the SYSTEM section).
Discussion
In this tutorial an old PVI file was converted to the latest version. This is important for files
using FIELD units, containing Differential Liberation (DL) experiments that have Gas
formation volume factor (GFVF) measurements as the units, for this type of observation was
changed in 99B from rb/stb to rb/Mscf to make PVTi’s units systems consistent with those of
the ECLIPSE simulators.
General information
• "Main PVTi window" on page 90
• "File" on page 92.
• "View" on page 95
• "The fluid model" on page 98
• "COMB - Compositional Material Balance" on page 112
• "Simulation using PVTi" on page 117
• "Regression in PVTi" on page 126
• "Exporting keywords" on page 133.
• "Utilities" on page 144.
• "Batch system and keywords" on page 152.
• "Error handling" on page 165.
General information
PVT analysis involves fitting an Equation of State to experimental data and then using the
Equation of State to produce PVT tables for use in reservoir simulations. PVTi contains
facilities to allow you to import experimental data, fit the data to an Equation of State, and
finally produce the PVT tables for reservoir simulation studies.
The menu bar of the main PVTi window has the following options:
• "File" on page 92.
• "Edit" on page 94
• "View" on page 95
• "Run" on page 96
• "Utilities" on page 97
• "Graph" on page 97
• "Window" on page 150.
• "Help" on page 150.
The main window contains all the tools necessary for Equation of State model fitting.
Basic features
The Data Tree provides a view of the current project’s contents. Each fluid sample is identified
with its experiments as sub-nodes in the tree. Likewise, each experiment has its observations as
sub-nodes. The Log Window is updated with pertinent information relating to actions taken in
PVTi. The Equation of State, upon which the current fluid model is based, is indicated in the
status bar. The Main Plot Window and the Sub-plots provide an area for viewing project
information graphically.
Hint The .NEW file contains all the changes made during the last session. If you change the
file extension to .PVI you can use it to recover the session.
Note Although only one project may be in use by the program, different sections of different
.PVI files may be read in.
Library
The properties of library components are preset by the program.
To display the current list of library components select View | Library...
C1 Methane C2 Ethane
C3 Propane C4 Butane
C5 Pentane iC 5 Iso-Pentane
C6 H6 Benzene C7 H8 Toluene
C7 Heptanes C8 Octanes
C9 Nonanes C 10 Decanes
C 11 Undecanes C 12 Dodecanes
C 13 Tridecanes C 14 Tetradecanes
C 15 Pentadecanes C 16 Hexadecanes
C 17 Heptadecanes C 18 Octadecanes
C 19 Nonadecanes C 20 Eicosanes
C 23 C23’s C 24 C24’s
C 25 C25’s C 26 C26’s
C 27 C27’s C 28 C28’s
C 29 C29’s C 30 C30’s
C 31 C31’s C 32 C32’s
C 33 C33’s C 34 C34’s
C 35 C35’s C 36 C36’s
C 37 C37’s C 38 C38’s
C 39 C39’s C 40 C40’s
C 41 C41’s C 42 C42’s
C 43 C43’s C 44 C44’s
C 45 C45’s
Note For components C 6 to C 45 , the properties stored in the internal library correspond to
the “grouped” properties of Single Carbon Number Groups (SCN), [Ref. 5].
Obvious candidates for the pseudoisation of components for use in large regressions or
compositional simulation are iso-butane and normal butane, and iso-pentane and normal
pentane, into single butane and pentane components. A study of many PVT reports [Ref. 19]
has shown that the typical ratios of iC4 : nC 4 , iC5 : nC 5 are 0.67:0.33 and 0.60:0.40 respectively.
The library also contains two other components, with the mnemonics C4 and C 5 , whose
properties are mole-weighted averages of the respective iso and normal component properties.
Run
The Run menu provides simulation and regression facilities.
The following options are available:
• Check Fluid System.
This provides a consistency check of the current fluid, the results of which are posted to the
log window. If there are a lot of fluid errors, the results are also displayed in a text window.
• Save As Samples .
If this option is turned on, any samples created by an experiment can be saved as additional
project samples.
Utilities
The Utilities menu option provides access to miscellaneous information relating to the project
and program set-up.
• Units...
See "Units..." on page 144.
• Standard Conditions...
See "Standard conditions..." on page 145.
• Program
This opens the sub-menu of program configuration options.
• Options... This opens the Options panel which mimics the OPTIONS keyword in the
PVI file. See "Program options" on page 145.
• Debug... See "Debug..." on page 150.
• Text Editor
This opens the text editor used for displaying simulation results, etc. It can be used to view
any ASCII file.
Graph
The Graph menu provides options related to the plotting of graphs.
• Add New Graph...
Adds a new graph to the existing plot windows.
• Superimpose
When the superimpose option is switched on, indicated by a tick next to the menu option,
subsequent graphs are superimposed on the current main graph.
• Tabulate...
This option creates a table showing the values plotted in the current main graph.
• Remove All
This option deletes all graphs from the window.
Hint You can load the first two sections of a PVI file as a system specification, rather than
using menu options. Additionally, you can choose to echo the contents of the whole
PVI file to the current print file, PVP.
Note EoS regions each have an EoS model defined within them that is an EoS plus a list of
critical properties defined for each component. An Equilibration region is a group of
cells where the initial pressure and saturation is defined. PVTi needs to know which
Equilibration region to read in if there are any composition versus depth (specified by
the ZMFVD or COMPVD keywords) tables in the ECLIPSE file. There is one table for
each of the Equilibration regions. By specifying which Equilibration region to use this
tells PVTi which table to read in.
Hint If you are uncertain as to the contents of the selected section, use File | View PVI
Section to display the section to the screen.
The syntax of the external file is similar to that of ECLIPSE. The data file is free format, except
for keywords, which must start in column 1. For further information on the keywords see "PVTi
keywords" on page 167. An example of such a file for a trivial two-component CO2 -isoButane
system is as follows:
---
---- RUNSPEC section: specific number of components and the EoS
----
RUNSPEC
NCOMPS
2 /
EOS
PR /
----
---- SYSTEM section: define hydrocarbon properties and composition
----
SYSTEM
---- Unit conventions
UNITS
METRIC ABSOL FRACTION /
DEGREES
KELVIN /
---- Component names (library defaults)
LNAMES
CO2 IC4 /
---- Overwrite default omega values by component
OMEGAA
0.4572 0.4572 /
OMEGAB
0.0778 0.0778 /
---- Initial sample composition
ZI
0.6 0.4 /
---- Binary Interaction Coefficients
BIC
2 1 1 0.13 /
/
----
---- This defines the fluid, EoS etc.
---- COMB, SIMULATE, REGRESS, BLACKOIL sections may now follow
---- See examples in Appendix C
----
END
Equation of State
Hint The default Equation of State is the Peng-Robinson three-Parameter equation. This is
suitable for most requirements, so generally you do not need to set the equation of
state.
This panel allows you to choose one of five main equations of state, to specify the required
viscosity correlation, and to decide whether or not to activate editing of specific heat capacities.
The equations of state are described in "Equation of state" on page 316; the viscosity options
are described in "Viscosity correlations" on page 329 and in [Ref. 5], [Ref. 7] & [Ref. 41].
The available equations of state are:
• Peng-Robinson
• Soave-Redlich-Kwong
• Redlich-Kwong
• Zudkevitch-Joffe
• Schmidt-Wenzel
Hint It has been our experience that the Peneloux et al. three-parameter equations of state,
PR3 and SRK3, generally give much better predictions of liquid properties and
saturations.
They also allow you an additional set of regression parameters, namely the component
volume shifts, making for an easier match to measured data.
Viscosity correlations
The Lohrenz-Bray-Clark, Pedersen and Aasberg-Petersen viscosity models are available. Select
the appropriate viscosity model.
Note You can re-select the equation of state or viscosity correlation at any stage. However,
the default EoS parameters for each component are dependent upon the EoS, and the
program re-initializes these if you change the EoS. Alternatively you can manually
reset the parameters to the default values at any time.
Components
PVTi: Edit | Fluid Model | Components...
User
This option allows you to define components. Enter the required properties into the panel:
critical pressure and temperature, acentric factors, etc.
You should enter the components in order of increasing molecular weight, and non-
hydrocarbons before hydrocarbons:
Non-Hydrocarbons
• H2
• H2 O
• CO
• N2
• H2 S
• CO 2
Hydrocarbons
• C1
• C2
• CN+
Hint By selecting PVTi: Run | Check fluid system the fluid is re-ordered into increasing
mole weights.
PVTi allows you to input a user component even if you know only the critical temperature and
pressure. It calculates the other properties as follows:
• Tb from T c and P c - Riazi-Daubert. For further information see [Ref. 30] and [EQ 8.11],
[EQ 8.12].
• Sg from T c and T b - Riazi-Daubert. For further information see [Ref. 30] and [EQ 8.11],
[EQ 8.12].
• Mw from T b , T b and Sg - Riazi-Daubert. For further information see [Ref. 30] and [EQ
8.33], [EQ 8.34].
• ω from P c , T b and T c - Edmister. For further information see [Ref. 30] & [EQ 8.10].
• P from Macleod and Sugden. For further information see [Ref. 14].
Use the Update button to calculate the other properties of the component.
Characterization
If you give a characterization, you must generally specify at least two out of the following (these
are specified in the Characterization folder):
• molecular weight M w ,
• specific gravity Sg ,
• reference temperature K
Hint If you have more than two of the set M w , S g , T b and K w , we recommend that you enter
the best two first, as the order of entry decides which pair the program selects. For
example, if you enter M w , T b and K w then the program uses M w and T b .
You can choose from the following correlations for estimating the physical properties and
acentric factors:
Critical properties
• Kesler-Lee. See [Ref. 10].
• Cavett. See [Ref. 11].
• Riazi-Daubert. See [Ref. 12].
• Winn. See [Ref. 43].
• Pedersen. See [Ref. 43], [Ref. 41] and [Ref. 45].
Acentric factors
• Kesler-Lee. See [Ref. 10].
• Edmister. See [Ref. 14].
• Thomassen. See [Ref. 30].
• Pedersen. See [Ref. 43], [Ref. 41] and [Ref. 45].
Note When reading in a file the critical volumes (Vc) and critical Z factors (Zc) for each
component must satisfy the relation PcVc=ZcRTc (where Tc, Pc, R are the critical
temperatures, critical pressures and universal gas constant respectively). If this is not
the case then PVTi will alter the values of the relevant critical Z-factors in order that
this relation is satisfied.
Volume shifts
Note Volume shifts are only available if you use a three-parameter Equation of State.
Use this option to enter the dimensionless volume shifts. The actual volume shifts in the
equation of state are displayed beside them.
1 PVTi: Edit | Fluid Model | Volume Shifts... .
This displays the Volume Shifts and Thermal Expansion Coefficient panel.
2 Enter the volume shifts for the required components.
3 Click on OK.
Note If the "Temperature dependence for volume shifts" on page 147 option is set then you
can enter a value for THERMX, the thermal expansion coefficient.
Thermal properties
Note You can only use this option if the program option "Specify/Calculate density and
molar volume units" on page 147 is switched on. It can be switched on in the Equation
of State panel (see "Equation of State" on page 100).
Specific heat capacity coefficients and calorific values for each component are the thermal
properties used in PVTi.
1 PVTi: Edit | Fluid Model | Thermal properties....
This opens the Thermal Properties panel.
2 Amend the thermal properties for the components, as required.
3 Click on OK.
Note LBC Viscosity coefficients are only available if you are using the LBC Viscosity
model.
Splitting components
This menu allows for the automatic splitting of the plus fraction into a required number of sub-
fractions for subsequent use in a large regression or for output to a compositional simulator such
as the one in ECLIPSE.
Splitting is also used to accommodate different plus-fraction properties for different fluid
samples. This process is often known as a multi-feed split.
This option allows you to input data for splitting the plus fraction. There are three methods
available from this option for splitting the plus fraction, which must be the last component:
• Constant Mole Fraction splitting (CMF)
• Whitson
• Multi-feed split or Semi-Continuous Thermodynamic (SCT) splitting
1 To choose the splitting method, select PVTi: Edit | Fluid Model | Split and select the
splitting option.
Note Note that splitting is not necessarily the opposite of grouping. Splitting the plus
fraction into two or more pseudo-components, followed by a re-grouping of those
pseudo-components back into a single plus fraction, generally results in a different set
of critical properties, etc., from those originally possessed by the plus fraction.
PNA Distribution
This splits all components heavier than the library C6 component into paraffinic (P),
naphthalenic (N), and aromatic (A) components. This is done according to the method outlined
in "The PNA distribution of heavy components" on page 394. The critical properties assigned
to the PNA components are those described in "Critical properties of PNA species" on
page 395.
Group
This menu allows for the automatic grouping of sub-fractions for subsequent use in a large
regression or for output to a compositional simulator such as the one in ECLIPSE.
This option allows you to choose components to group and perform the grouping operation. The
default scheme for grouping is to group to the default sample ZI using the mole fraction
weighting to group components. Other schemes of grouping include grouping by molecular
weight and by mixing rule, see [Ref. 44]. Also the sample to group to can be changed to any in
the defined set, or to an average of all samples.
To group components, select PVTi: Edit | Fluid Model | Group...
This displays the current component system, each component having an associated index. The
first time you enter this option, all these indices are set to 0, indicating that they do not belong
to any group.
1 To create a new pseudo component, give a new index of greater than zero to two or more
components.
2 Select the Grouping Technique.
3 Give the group or pseudo-component a new component mnemonic, if required.
Hint You can perform several groupings from the same original component description by
specifying the new components with ascending indices, 1, 2, etc.
Note Note that splitting is not necessarily the opposite of grouping. That is, splitting the plus
fraction into two or more pseudo-components, followed by a re-grouping of those
pseudo-components back into a single plus fraction, generally results in a different set
of critical properties, etc., from those originally possessed by the plus fraction.
Defining Samples
Sample names
1 PVTi: Edit | Samples | Names...
Defines sample names. Use this option to enter mnemonics for each component. You can enter
more than one sample for later use; to do this, reference each sample by its mnemonic, of up to
8 characters.
Note Note that the mnemonic for the default sample is “ZI”, for “z initial”.
For alternative samples, you may specify a line of text to give additional information. For
example: from different depths in the hydrocarbon column, a “saved” calculated composition
from a simulation, etc.
Sample compositions
1 PVTi: Edit | Samples | Compositions....
Enter the compositions for each defined sample. PVTi ensures that they all add up to unity. If a
sample does not add up to unity, a message appears asking whether or not the program should
redistribute the difference across the components.
Sample salinities
1 PVTi: Edit | Samples | Salinities...
This option allows you to enter sample salinities. If you have entered H 2 O as a component then
use this option to add the salinity of each sample.
Note This information is used in the MFLASH experiment in the "Simulation using PVTi"
on page 117.
Viewing samples
Sample fingerprint plot
This option allows you to generate fingerprint plots. This consists of plotting the logarithm of
the component mole fractions against the component molecular weights.
Hint Fingerprint plots give an idea of the nature, that is condensate or volatile oil, of a given
fluid sample. Providing a reasonable split of the Heptanes plus is available, then a
condensate typically has straight line or down-turning slope proceeding towards the
heavier fractions, whilst a volatile oil has an up-turning slope as it usually contains
more heavy fractions.
Hint A default phase plot (with one quality line) can be generated by simply dragging a
sample name from the Data Tree and dropping it into the Main Plot Window.
Note If depletion experiments or separators exist, they are plotted onto the phase plot too.
Hint The default component groupings are: C1 and the non-hydrocarbons, C2-C6 and C7
and heavier hydrocarbons. This is typically the best choice, so you should only need to
change the groupings from the default in special cases.
Introduction
Material balance checks provide an important tool in analyzing the quality of the data found in
a PVT report. In PVTi, material balance checking is provided for Constant Volume Depletion
(CVD) experiments, Separators by recombination, and the calculation of liquid and vapor
recovery.
Caution It is important to always examine PVT laboratory data for material balance errors.
If the reported observations contain serious errors, these will be reflected in the
fitted equation of state model.
It is generally the case that for gas condensate and volatile oil samples, a constant volume
depletion experiment (CVD) is performed as part of a laboratory analysis. Using material
balance considerations, it is possible to calculate liquid compositions and hence K -values,
molar masses, densities, etc. This can be very useful for two reasons:
1 With appropriate separator data, estimates of oil and gas recovery can be performed without
recourse to an equation of state model.
2 The consistency checks provide a measure of the quality of laboratory data and
consequently its value or otherwise in any subsequent regression analysis.
Plot
Use this option to perform plots.
The set of plots available from this option are:
• Vapor composition (input)
• Liquid composition
• K -values (1): log( K ) versus pressure
• Produced moles
• Vapor moles left in cell
• Liquid moles left in cell
Modify
Use this option to modify data. Errors in the input data may become evident after you performed
the CVD material balance. You can use this option to rectify these errors. You may consider
changing the saturation pressure liquid saturations, Z -factors, etc.
Hint You can avoid calculated liquid mole fractions remaining in the cell by changing the
composition of the feed stream (well stream) or one or more of the removed gas
streams.
Separator recombination
Use this panel to directly test the quality of recombinations. To open this panel,
right-click on the separator experiment in the project-tree and select Material Balance... from
the pop-up menu.
Report
Use this option to test recombinations. It generates the calculated and Standing's K -values and
provides a report on the results.
Testing combinations
If sufficient volumetrics data has been entered in the previous item the recombined sample is
calculated and output to the report and PVP file.
The output details the calculated feed to either separator using the given values of GOR, oil
density and liquid and vapor compositions from each separator stage and is printed together
with the given feed composition for comparison.
Plot
Use this option to plot the recombination results. It produces a Hoffmann-Crump-Hocott plot,
which gives a measure of the quality of the separator data.
Note This option is only available if there are valid depletion and separators experiments in
the project, in addition to the CCE experiment.
Report
Use this option to perform recovery calculations. It allows you to estimate recoveries of vapor
and liquid.
Plot
Use this option to produce a plot of the gas and oil recovered as a function of pressure.
Introduction
PVTi allows you to perform experiments on your defined hydrocarbon system. The experiments
available include:
• Flash calculation
• Bubble and dew point determination
• Simulation of laboratory experiments such as constant composition expansion, constant
volume depletion, and differential liberation
• Swelling tests
• Separators
• Variation of composition and pressure with depth
• Vaporization tests
• multiphase flash.
• Wax appearance temperature
• Asphaltene appearance pressure
Note You can enter any experimentally determined observations available to compare with
the output produced by the Simulate module, apart from the multiphase flash
experiment.
Defining Experiments
1 PVTi | Edit | Experiments...
This opens the Define Simulate Experiments panel.
2 To add a new experiment select it on the Add menu. To edit an existing experiment, select
its name from the drop-down list, which appears in the panel, and click OK.
3 Once an experiment has been created for definition or selected for editing, you will see a
customizable form that can be tailored to suit your data. The information is entered in a
number of folders.
4 In the first folder, General, you can select various single-valued pieces of information for
entry. The required information is automatically selected and cannot be deselected, so an
information line informs that this data is required.
5 The second folder, Observations, shows a table where you can customize column
headings to match your observation data. If you are editing an existing experiment, the
observations in the currently defined set are already shown as column headings. In this way
you can match the data entry panel to your own data-set. If you forget to enter a required
quantity (for example Pressure in a Differential Liberation experiment - see "Differential
liberation" on page 119) PVTi warns you and does not create the customized entry form.
Note The gas-oil ratio reported by the calculation is defined as gas volume at standard
conditions divided by liquid volume at flash conditions. The gas volume is obtained
using a Z -factor of unity.
Note If the temperature is such that no bubble point can be found (above the critical
temperature) the program returns a warning message.
Note If the temperature is such that no dew point can be found (above the critical
temperature) the program returns a warning message.
Note When obtaining relative volumes the program uses saturation volume as a
normalisation volume, if one exists, or the volume at the highest pressure, if not.
Hint You can apply this experiment to a liquid (bubble point) or vapor (dew point) system.
The program tests for both possibilities. It is not, however, possible to apply this
experiment to samples that are above the cricondentherm.
Note The relative volume reported by the program is the fraction of the cell filled with liquid
at the end of the constant volume step, that is after the original volume has been
restored by removing vapor.
Differential liberation
For this experiment you must specify a temperature and a series of pressures. You do not need
to give a value at the bubble point. PVTi provides this pressure point.
The program also provides automatically the last step in the differential liberation process, the
reduction to standard conditions. However, the program does not provide the pressure point at
standard pressure (usually 14.7psia) and at reservoir temperature and the user must enter this for
the final stage.
Note You may only apply this experiment to a liquid (bubble point) system.
Hint The relative volume reported by the program is the ratio of the oil volume at each step
to the oil volume at the final (standard conditions) step.
Note There are alternative definitions of the GOR and the relative oil volume available using
the program options "Definition of GOR in Diff. Lib." on page 148 and "Definition of
Oil relative volume in Diff. Lib." on page 149.
Swelling test
A swelling test consists of adding increasing amounts of a lean gas to a reservoir fluid and
determining the swelling of the mixture relative to the original fluid composition.
For this experiment you must specify:
Separators
Separators consist of a set of connected equilibrium flashes at user prescribed pressures and
temperatures.
For this experiment you must specify:
1 The composition of the feed stream from the defined sample mnemonics.
2 A number of stages (up to seven) for which you must give a pressure and temperature.
Additionally, you must connect the vapor and liquid outputs of each stream to some subsequent
stream. The default routing is to connect the liquid output of stage (j) to stage (j+1), and to take
the vapor output to stock tank conditions (as defined by the STCOND keyword or by the
Standard Conditions menu option under OPTIONS).
Note A stage output can be fed back to any previous stage though not back to the current
stage. No stock tank stage is defined automatically, and whereas it is customary to
quote vapor properties at stock tank conditions, liquid properties will be quoted at the
“final stage” conditions. Therefore, if liquid properties at stock tank conditions are
required, this should be the final (additional) stage which must be defined by the user.
For example if we have a separator with 3 stages with the last stage being stock tank
conditions, then a liquid FVF at stage 1 of the separator will be the volume of liquid
divided by the final liquid volume (stock tank conditions in this case) after flashing the
liquid feed of stage 1 through the remaining 2 stages of the separator chain.
Hint The "Definition of Oil relative volume in Diff. Lib." on page 149 program option
allows you to quote GORs as volume of gas at standard conditions per volume of stock
tank oil as opposed to the default calculation of volume of gas at standard conditions
per volume of separator liquid at separator conditions.
The program can calculate oil formation volume factors, that is the volume of reservoir fluid at
initial or bubble point conditions per stock tank volume (SRELV) and, by separator stage,
volume of separator liquid at separator conditions per stock tank oil volume (ORELV).
To use this option tick the box in the panel or use the FVFREF keyword in batch mode.
Note The assumptions made in the performance of this experiment, that there are no
asphaltenes and that the reservoir is in thermal, gravitational and diffusive equilibrium,
are probably not achieved in any real reservoir. However, despite these reservations,
this is a useful test of the depth-variation of a particular fluid.
Vaporization test
This is a somewhat specialised test performed for gas-injection on reservoir fluids, but in fact it
is rather similar to a swelling test.
For this experiment you must specify:
1 The composition of the reservoir fluid and injection fluid from the currently defined sample
mnemonics.
2 The reservoir pressure and temperature
3 The number of moles of injection gas to be added to the reservoir fluid.
Note The gas-oil ratio reported by the calculation is defined as gas volume at standard
conditions divided by liquid volume at flash conditions. The gas volume is obtained
using a Z -factor of unity. Note also no facility currently exists for comparing these
against observed values.
Note If the sample you select for the experiment contains water, you should enter the salinity
in the PROPS section.
Note In the case of the dew point, program returns the retrograde (highest pressure) dew
point.
Hint You can use this experiment in the REGRESS section. For further information see
"Regression in PVTi" on page 126.
Defining Observations
Observations can be defined at the same time as the experiment; see "Defining Experiments" on
page 117.
1 PVTi: Edit | Observations...
This opens the Define Simulate Observations panel.
2 Choose the index of the experiment
3 Select abbreviation for the observation to be entered
Abbreviation Observation
Liquid Z-Factor Liquid Z -factor
Vapor Z-Factor vapor Z -factor
Liquid Density Liquid density
Vapor Density Vapor density
Liquid Mol. Wght. Liquid molecular weight
Vapor Mol. Wght Vapor molecular weight
Liquid Visc. Liquid viscosity
Vapor Visc. vapor viscosity
Liquid Sat. Liquid saturation
Vapor Sat. Vapor saturation
Vapor Mol. Frac. Vapor mole fraction
Sat. Pressure Saturation pressure: gas - p dew , oil - pbub
Sat. Temperature Saturation temperature: not currently available
Gas-Oil Ratio GOR: SEPS - gas (STC)/oil (stage/STC); DL gas (STC)/oil (STC/ p sat )
Relative Vol. Relative volume (SWELL = swelling factor)
Total Gas-oil ratio Cumulative separator GOR: (Gas at STC/final stage Oil)
Ternary Plot Ternary Plot (Multi-Contact Test and plotting only)
Mole. wght. plus Mole weight of plus fraction (COMB Mat. Bal.)
Rel. oil sat. Relative oil saturated volume ( B o ( p bub ) in DL)
K Values K -values
Abbreviation Observation
Gas FVF Gas formation volume factor (DL)
Gas Vol. Ext. Gas volume extracted (at STC) (DL)
2-phase Z Two phase Z -factor (CVD)
Oil rel. FVF Oil FVF from p init ⁄ p bub to p stock (SEPS)
Note Note that not all the observed data types are available for all experiments.
Running a simulation
Note Simulations are automatically run on creation, so the results are immediately available.
Hint PVTi recognizes which experiment simulations are up to date and then only performs
necessary calculations. This means that to view the simulation results you should
always use PVTi: Run | Simulate.
Introduction
Performing a regression
To perform a regression you must specify:
• The experiments to be used in the regression. You can choose from the experiments
mentioned in "Defining Experiments" on page 117.
• The weighting for the observations associated with those experiments. You can use most of
the observations given in a laboratory experiment as observations to match against
predicted data.
• The Equation of State parameters you wish to vary to match predicted to observed
quantities. Most of the Equation of State parameters are available as regression variables.
Note The time taken for the regression operation rises rapidly with the number of variables
chosen, and the use of the minimum possible set is suggested. That said, any
combination of critical point data, Ωa and Ωb values, acentric factors, binary
interaction coefficients, δ ij and volume shift parameters (if the PR3 or SRK3 forms for
the Equation of State are being used) may be chosen to be modified. There is a
maximum total of 50 regression variables.
Regression function
The regression function to be minimized is the normalized root mean square (RMS) error of
predicted experiment results to the given (weighted) observed experiment results. See "Details
Folder" on page 131 for a description of the RMS value used.
Note In order to run the regression, there must be at least as many observations as chosen
regression variables.
Regression Panel
1 To open the Regression panel select PVTi: Run | Regression...
2 Use this menu to set up and perform regression. To open this menu, select PVTi: Edit |
Regression
Variables
The Variables section has the following options:
• Normal (component properties and BICs)
Normal variables are individual fluid component properties and binary interaction
coefficients. For further information see "Defining regression variables..." on page 127.
Report
The Report section has the following buttons:
• Regression: Opens the Regression Report panel. See"Regression Report" on page 130.
• Simulation: Opens the simulation report of all experiments.
Regress
The Regress section is for running regressions. The buttons in this section are:
• Run: Perform regression. For further information see "Performing a regression" on
page 130.
• Accept / Reject: Accept or reject the last regression. For further information see
"Accepting or reject regression results" on page 130.
Additional Information
Defining regression variables...
Use this option to define/re-define the set of variables for use in the regression. There may be
two sets of variables available for use in regression, depending on the state of the "Program
options" on page 145 and whether modified Whitson splitting (SCT) has been used on the plus
fraction. The two sets are denoted normal and special.
• The normal variables are the component dependent ones, that is variables such as critical
properties, acentric factors, etc.
• The special variables are system-wide or multi-component variables such as the thermal
expansion coefficient or the Cheuh-Prausnitz A -factor for binary interaction coefficients.
b Tc
c ω
d Ωa
Hint If you wish to vary a given EoS parameter, say T c , of two or more components as
separate independent quantities, you should give them different values. For example
1,2,..., etc.
b Vc
c Vc
Hint If you wish to vary the parameters as one or more groups of variables, you should give
the required group members the same integer. This may be particularly useful when
trying to vary V c values to match to viscosity data using the LBC correlation, for
example.
Note All of these data fields can take an integer value, 0,1,2,..., and so on. The default of
zero (or blank/null field) implies that the particular component’s EoS parameter is not
to be used in any subsequent regression.
Use the second panel to define the status of the lower half of the (symmetric, zero-diagonal)
matrix of binary interaction coefficients.
Note Note that the rules regarding choice of groups for binary interaction coefficients are
slightly different in that groups may be specified down columns or along rows of the
lower half matrix but not both.
Note Note that the mole weight and skewness variables apply on a sample-by-sample basis,
therefore there must be the appropriate number of measurements defined to allow this
option to be used.
Regression target
The regression target is the size of the objective function at which the regression terminates as
having achieved a match of calculated to observed values. Enter the regression target. The
default value is 0.000001.
Hint Normally you do not need to change any of these limits, except for the maximum
number of iterations (which might be reduced, especially for larger problems) and the
limits on some variables (for example to prevent violation of monotonicity
relationships).
Note For all variables specified in the regression the program displays the lower and upper
limits.
Note For all the equation of state parameters, except the acentric factor and binaries, these
lower and upper limits are scaled (to unity) variables, with default settings of 0.5 and
1.5. That is, the program allows the variable to decrease/increase by up to 50% before
it terminates the regression.
Running a regression
Use this option to perform regression.
Note This option is only available if you first define a set of experiments, observations and
variables.
Performing a regression
1 PVTi: Run | Regression
The program performs the regression after first checking that there are at least as many
regression points as variables.
Hint The results of regression can be viewed in the Sensitivity Analysis panel. See
"Regression Report" on page 130.
Regression Report
The regression report contains detail of the current fit between the model and the observations.
Hint In any regression, having a few very sensitive parameters is preferable to having
hundreds of insensitive ones. Always look for parameters that can be discarded. This
is called conditioning the problem - an ill-conditioned problem is difficult to solve.
Details Folder
The first folder, Details, shows the current fit. The numbers at the top of the folder show the
Total (normalized) RMS fit and the Weighted (normalized) RMS fit.
Hint The RMS values in PVTi are normalized (the difference between the observed and
calculated values is divided by the observed value) so that, for example, pressures
(which could be thousands of psi) are treated similarly to saturations (which are
between zero and one). To calculate the weighted RMS, each normalized difference is
multiplied by the assigned weighting. The default weight is 1 and so, initially, the two
RMS values are the same.
The remainder of this first folder shows the observations, their weighting in the regression, and
the percentage difference between the observed and calculated values.
Modifiers folder
The second folder, Modifiers, shows the selected regression parameters, the minimum and
maximum allowed values for each modifier, and the percentage change made to the modifier
during regression (initially this is zero as no regression has been run).
Sensitivities folder
The Sensitivities folder shows the sensitivity of each observation to changes in each regression
variable.
Hint A large sensitivity indicates that changing that regression variable has a large effect on
the fit to that observation. Likewise, consider discarding regression variables to which
the observations are insensitive, since a large modifier is needed to obtain a fit and, in
general, large modifications lead to unrealistic fluid models.
Hessian folder
The Hessian matrix is a good indication of the conditioning of the inversion problem
(regression).
Hint In a well-conditioned problem the leading diagonal of the Hessian matrix is dominant.
Specifically, look at an element on the leading diagonal. If it is larger than the other
values in that row (or column) of the Hessian, then it indicates that the regression is
likely to succeed.
Correlation folder
The Correlation matrix is very important as it can indicate links between regression variables
that are not obvious on first inspection. If two variables are strongly correlated (correlation close
to 1), they both move the fit in the same direction; and so changing one is similar, in effect, to
changing the other. If two regression variables are strongly anti-correlated, changing one has the
opposite effect to changing the other. This latter case can cause a difficulty in regression as the
two variables could be changed an unlimited amount in opposite directions without having a
noticeable effect on the fit.
Hint Look at off-diagonal elements in the correlation matrix. If any are close to 1 or -1,
consider removing one of the two regression variables that are correlated. This
improves the likelihood of a good final fluid model being created.
(If the variables are of the same type, for example if they are both Tcrit, you could
consider combining them into a single regression variable.)
General information
PVTi can be used to generated output for ECLIPSE BlackOil, ECLIPSE Gi option, ECLIPSE
Compositional, ECLIPSE Thermal, VFPi and the API Tracking option in ECLIPSE BlackOil.
The Export modules are used to produce models for the ECLIPSE simulators (see "Output for
ECLIPSE simulators" on page 353 for background information). You can generate data files for
exporting to ECLIPSE BlackOil , ECLIPSE Compositional, ECLIPSE Thermal, API Tracking
and VFPi from PVTi.
To open this menu, select PVTi: File | Export.
The Export menu has the following options:
• ECLIPSE Compositional Fluid Model...
See "Export for ECLIPSE Compositional" on page 133.
• ECLIPSE Thermal Fluid Model...
See "Export for ECLIPSE Thermal" on page 134.
• API Tracking option in ECLIPSE BlackOil...
See "Export for API Tracking option in ECLIPSE BlackOil" on page 134.
• Oil reservoir...
See "Export Oil Reservoir" on page 135.
• Gas reservoir...
See "Export Gas Reservoir" on page 136.
• Equilibration...
See "Export Equilibration" on page 136
• Water...
See "Export Water" on page 137.
• VFPi
Use this option to generate VFPi tables. For further information see "VFP module" on
page 138.
Units
Allows you to export in the ECLIPSE unit set of choice.
Sample
If a sample is selected, it will be exported in the ZI keyword.
Reservoir temperature
This is exported in the RTEMP keyword.
Sample
The sample for which the keywords will be exported for.
Max/Min Pressure/Temperature
You are recommended to enter the maximum and minimum values of pressure and temperature
for your reservoir.
Units
Allows you to export in the ECLIPSE unit set of choice.
A detailed explanation of the workflow required for this export option is described in the
"Compositional Data for ECLIPSE Thermal" on page 367. For a technical review of PVTi’s
export facility for ECLIPSE Thermal see "ECLIPSE Thermal Export Module" on page 401.
Properties Keyword
This sets the keywords that are exported, for example Live Oil (PVTO) and Dry Gas (PVDG).
Hint All the keywords are described in the "ECLIPSE Reference Manual".
Plot Results?
You can tell PVTi not to plot the tables if you wish. PVTi only has room for 8 plots so if there
are many samples in the list then it may be useful not to plot the tables.
Separator Experiment
You can select a separator in your project if you wish instead of the default standard conditions
separator (usually 14.7psia and 60F).
Units
Allows you to export in the ECLIPSE unit set of choice.
Units
Allows you to export in the ECLIPSE unit set of choice.
Properties Keyword
This sets the keywords that are exported.
Hint All the keywords are described in the "ECLIPSE Reference Manual".
Separator experiment
If creating tables for ECLIPSE BlackOil, the output from the Differential Liberation can be
passed through any separator. The default is to use a single stage at standard conditions (that is,
a stock-tank only - normally 14.7 psia, 60F).
Units
Allows you to export in the ECLIPSE unit set of choice.
Properties Keyword
This sets the keyword that is exported.
Hint All the keywords are described in the "ECLIPSE Reference Manual".
Separator experiment
If creating tables for ECLIPSE BlackOil, the output from the Constant Volume Depletion can
be passed through any separator. The default is to use a single stage at standard conditions (that
is, a stock-tank only).
Injection fluid
For the Gi option only, this provides the injection sample to be used in creating the pseudo-
compositional keywords. See "Pseudo-compositional tables for ECLIPSE GI option" on
page 361.
Export Equilibration
This panel exports keywords from a Composition versus Depth experiment.
The keywords are suitable for the ECLIPSE BlackOil INIT section and ECLIPSE
Compositional PROPS section.
Units
Allows you to export in the ECLIPSE unit set of choice.
Properties Keyword
This sets the keyword that are exported.
Hint All the keywords are described in the "ECLIPSE Reference Manual".
Export Water
Description
This panel exports water properties.
The keyword PVTW or PVTWSALT can be exported for the ECLIPSE simulators. The keyword
Hint All the keywords are described in the "ECLIPSE Reference Manual". WATPVT is
described in the "VFPi User Guide".
Units
Allows you to export in the ECLIPSE Unit set of choice.
Reservoir temperature
This sets the reservoir temperature to be used in generating the keywords.
Reservoir pressure
See above.
Brine
If this option is checked, salt concentration(s) can be accounted for.
VFPi
If this option is checked, the exported keyword will be WATPVT. The low temperature and low
pressure fields become active, these values correspond to the lower values to be used in the flow
table (for example, the top of the pipe).
For more information see "Water properties" on page 362.
Introduction
The VFP module generates blackoil tables for VFPi. To open this menu, select PVTi: File |
Export | VFPi from the main PVTi window.
Figure 6.7 The VFP module
VFP toolbar
The module toolbar contains the following buttons:
Display section from PVI file. Use this option to display VFP sections from the PVI file.
For further information see "Display VFP section from PVI..." on page 140.
Experiments
Use this option to define VFP experiments. For further information see "Experiments..." on
page 141.
Generate.
Use this option to select the method of generation for the VFPi blackoil tables. For further
information see "Method" on page 142.
Simulate
Use this option to generate the tables. For further information see "Perform" on page 142.
View results
View plots
Use this option to view plots of the tables. For further information see "Plot" on page 143.
File
Use this menu to read and display VFP sections from a PVI file.
The File menu has the following options:
Note If the program does not find a VFP section in the PVI file it produces an error message
and stops processing.
Print preview
Previews the printed plot in the main module area.
Print layout
Allows you to set the text styles and sizes for printing.
Print setup
Sets up the printer specifications. This option is specific to Windows.
Print
Prints the whole module print area, the main window and allows you to set the print type, for
example color, postscript.
Close
Use this option to close the VFP module. To close the VFP module, select VFP | File | Close.
Note The data is not lost. If you reselect the module all data that was previously set is
available.
Experiments...
Use this option to define experiments.
Defining experiments
1 Select VFP | Define | Experiments... .
This opens the Constant Composition Experiment panel.
2 Define a Constant Composition Expansion (CCE) experiment.
This is a depletion experiment where the moles entering the wellbore must leave the
wellbore. For further information on defining experiments see "Data requirements for the
experiments" on page 118.
Hint The first temperature you give for the CCE should be the highest temperature in the
production string, which can safely be taken to be the reservoir temperature. The
second temperature you give should be the lowest temperature in the production string,
for example a sub-sea temperature of 4 °C .
Note You do not need to give a saturation pressure as the program gives this data. It also
calculates undersaturated properties for all pressures. See "Output for ECLIPSE
simulators" on page 353 for further details
3 Define the separator conditions that take the reservoir liquid and/or vapor to surface
conditions to define the various ratios in the blackoil tables.
Note Note that to comply with the blackoil definition strictly, the last stage in the separator
train should be at standard conditions.
4 Specify the composition of the injection gas using one of the currently defined sample
mnemonics.
5 Give the gas-oil ratios in which the lean gas is to be added to the reservoir fluid.
6 Press OK.
Generate
This menu allows you to generate blackoil tables for VFPi. To open this menu, select VFP |
Generate.
The Generate menu has the following options:
• Perform
Generates the blackoil tables.
Perform
Use this option to generate blackoil tables.
Note You must define a depletion experiment and the separator configuration (which only
need be a one stage default system to standard conditions) to generate the blackoil
tables.
Note If the process is successful the program also writes the data to the PVO file which you
can add to an VFPi input file.
Note Note each set of tables is repeated twice, for the set of pressure nodes defined on the
CCE experiment, at the specified high (reservoir) and low temperatures.
Water properties
1 Click Yes when asked if you wish to generate water properties, using the VFPi keyword
WATPVT.
Plot
Use this option to produce plots.
Producing plots
1 Select VFP | View | Plot.
You can generate any of the following plots:
• oil formation volume factor, Bo
• oil viscosity μ o
• gas viscosity μg
• oil Rs
Help menu
Gives access to help with PVTi.
Help
Opens the context-sensitive on-line help panel.
Introduction
This menu allows you to change or redefine various program settings. It has the following
options:
• Units...
This allows you to set unit types. For further information see "Units..." on page 144.
• Standard Conditions...
This allows you to set the standard temperature and pressure. For further information see
"Standard conditions..." on page 145.
• Program Option
This allows you to change the program options For further information see "Program
options" on page 145.
• Debug...
This allows you to set the debug flags. For further information see "Debug..." on page 150.
Units...
This option allows you to select various unit types. You can:
• Choose the unit type within PVTi.
• Choose the temperature unit type; this can be different from the one selected by the above.
• Set mole fractions or percentages.
• Select absolute or gauge pressures.
For further information on the units see "Units" on page 409.
Program options
This option allows you to set various program options. Various options which cannot be set
elsewhere in PVTi have been collected together under this option.
1 To set the program options, select Utilities | Program Option.
This opens the Set PVTi Program Options panel.
You can set the following:
Component Library
This option allows you to specify the component library to use. The choice is between the Katz-
Firoozabadi and the Oil PVTi Library. The Katz-Firoozibadi is the default and recommended
choice.
Experiment compositions
This option allows you to turn off lengthy output to the print file of liquid and vapor
compositions calculated in experiments.
1 Switch off the output of the liquid and vapor compositions if they are not required.
PVTi prints the liquid and vapor compositions to the screen and PVP file, by default.
Experiment results
On definition of fluid properties, completion of experiments, etc., you can optionally choose to
write all the data to the PVP print file.
1 Select the appropriate output option.
By default, PVTi always writes all the data to the PVP print file; however, you can choose
to make this write optional or never done.
Plot vectors
For more advanced graphical manipulation of PVTi plots, you may choose to output data
vectors to a file in GRAF user vector format. The vectors are written when the plots are
performed.
1 Choose whether you wish to produce the vectors.
PVTi does not produce these vectors by default.
Flash calculations
This option allows the Flash and Psat algorithms to be changed from the default ECLIPSE
Compositional algorithms to the old PVTi (pre-99B) algorithms.
1 Choose to use the old PVTi algorithms for Flash and Psat.
Note This option must be set in the ON state for you to use these regression variables. You
must have characterized and/or user components, or be mixing samples, for the
variables to become of use.
Debug...
This option sets the debug flags.
Note This is a programmer test facility to request additional information from the program
at a debug level. You do not need to set one or more of the debug flags unless help is
required in tracing an apparent anomaly.
Monitor option
This is a programmer test facility to trace the root of a problem in the program at a subroutine
level.
Note You do not need to need to set this flag ON unless asked to do so.
Window
This menu allows you to control the size and appearance of the program windows. To select this
menu, select Window from the main PVTi window.
Tile
This puts all the visible windows into a “tiled” formation.
Cascade
This puts all the visible windows into a “cascade” formation.
Minimize children
This option minimizes all windows except the main PVTi window.
Restore children
This opens all minimized windows so that they are visible.
Help
This menu gives you access to the help for PVTi. To select this menu, select Help from the main
PVTi window.
ToolTips enabled
The ToolTips provide a single line of text about each toolbar button when the mouse pointer is
stationary over the icon. This option turns the ToolTips feature on and off.
General information
For the 2004A version of PVTi the batch mode has undergone a significant revamp. Over the
last few years the user interface of PVTi has evolved rapidly and many of the PVTi sections are
now no longer written out to the PVI file as they are no longer required when PVTi reads in a
file. For example a BLACKOIL section is no longer written out by PVTi when black oil tables
are exported using the user interface. Although they can be inserted into the file manually it was
felt that a more user friendly way of constructing the .PVI files containing all the appropriate
sections required for the batch mode was needed.
The next section outlines how it is now possible to perform a workflow interactively in order
for PVTi to be able to reproduce this workflow in batch mode at a later time.
The way in which you run a file in batch mode has also changed. Pre-2004A you had to enter
the keyword TESTCASE anywhere in the RUNSPEC section of the PVI file. You then launched
batch mode the command line, for example using $pvti filename (if using a PC).
This has changed for the 2004A release. It was felt that opening files and adding keywords and
then remembering to remove them at a later data was cumbersome. The new way to run in batch
mode is to launch PVTi from the command line and specify the word ‘-batch’ before the
filename in the line command instruction.
Hint To launch the file TEST.PVI to run in PVTi’s batch mode on a PC use the command
‘$pvti -batch TEST.PVI’. On a UNIX machine ‘@pvti -batch
TEST.PVI’ runs the file TEST.PVI in batch mode.
Sometimes users have more than 1 version of PVTi installed. It is also possible to specify which
version of PVTi to use on the command line using the -ver command.
Hint On a PC, to specify the 2004A version of PVTi to run the file TEST.PVI in batch
mode, use the command ‘$pvti -ver 2004a -batch TEST.PVI’.
Note The workflow must be pre-defined. The batch mode is not designed to be able to
reproduce workflows where you have been experimenting with a particular project.
1 Start PVTi in interactive mode with the file that you eventually wish to run in batch mode.
2 Save the file concisely by doing File | Save (concise)....
This erases the history of any workflows stored in the file.
3 Open the Options panel using the File | Utilities | Program | Options....
4 In PVTi 2004A a new option has been added at the bottom called Write Keywords for
Batch Mode. Select Yes and then close the panel.
5 Perform the required workflow and, when finished, save the file (but not concisely) using
the File | Save... option.
Note The file must not be saved concisely as this would erase all the history in the .PVI file
that PVTi uses to reproduce your workflow in batch mode.
6 The final task is to actually run the file in batch mode. To do this on a PC launch PVTi from
the command line using the statement ‘$pvti -batch filename’ where filename is
the name of your PVTi project, for example TEST.PVI.
Note The word ‘-batch’ can be put after or before the filename but it has to be somewhere
on the command line in order to tell PVTi to run in batch mode.
When you are running in batch mode the program automatically sends printed output such as
experiment simulation results to a print file with the same root name as the input .PVI file. For
example, if the input file is CRUDE.PVI, the print file is CRUDE.PVP. See "PVI file" on
page 155 for further information. In batch mode all the experiments are automatically simulated
in the project by default. If any regression is to be performed in the batch mode then the
experiment simulation is performed after this has been done.
If keywords are exported for ECLIPSE during the batch run then .PVO files are created as
normal, but they are named using a convention that is outlined in the next section.
A file is also created called BATCH_OUT.PVI. This is a saved version of your project after all
the steps in the workflow have been performed. It can be useful to have this file after the batch
run has finished if the fluid model has changed (for example during regression) in your project
during the course of the workflow
Regression
You are allowed a maximum of two REGRESS sections in the batch mode file. This is so that
regression on both special and normal variables is possible. Any further regression sections are
ignored by the batch mode.
Splitting/Grouping
You are allowed two GROUP sections and two SPLIT sections in the file. If there are REGRESS
sections then one SPLIT and one GROUP section is allowed before the first REGRESS section
and one SPLIT and one GROUP section is allowed after this REGRESS section.
Export
It is assumed that all exports would be performed at the end of the workflow. An unlimited
amount of export sections (such as BLACKOIL, OUTECL3 sections) are allowed as long as they
are after the last REGRESS section. Because multiple exporting is allowed a naming convention
has been invented to stop PVTi just writing each exported .PVO file over the top off the last one
that was written out.
The naming convention depends on what kind of export is being performed:
BLACKOIL
If a BLACKOIL section has been read in the .PVI file then the naming convention of the .PVO
file is “filename_samplename_experimenttype_keyword1keyword2.PVO”
where:
filename is the rootname of the project,
samplename is the sample name used for the export,
experimenttype is the type of depletion experiment used in the export,
keyword1 is the name of the first keyword exported
and keyword2 is the name of the second keyword exported.
For example if the file CRUDE.PVI was used to export sample ZI based on experiment DL1
using Live Oil and Dry Gas keywords then the name of the produced file would be
CRUDE_ZI_DLLIVEDRY.PVO.
OUTECL3
If an OUTECL3 section has been read in the .PVI file then the naming convention of the
outputted .PVO file is filename_samplename_FLUIDMODEL.PVO.
APITRACK
If an APITRACK section has been read in the .PVI file then the naming convention of the
outputted .PVO file is filename_experimenttype_APITRACK.PVO.
Note The COMB, PSEUCOMP and VFP sections-type keywords cannot currently be
prepared in a .PVI file for batch mode using interactive mode.
PVI file
The PVI file consists of a number of sections, each introduced by a section keyword.
• The first section must be RUNSPEC, which specifies the number of components, the
equation of state option and the run title.
• The SYSTEM section must follow the RUNSPEC section.
• Other sections may be in arbitrary order, and may occur more than once.
Keywords by section
Each of the main sections, for example the RUNSPEC section, has its own specific keywords.
Note the following:
• The keywords UNITS, DEGREES and STCOND are normally be specified in the SYSTEM
section, but can occur elsewhere, and can occur more than once.
• The keywords DEBUG and OPTIONS can be specified in the RUNSPEC section or in the
SYSTEM, COMB, SIMULATE or REGRESS sections.
• The keyword MESSAGE can be specified anywhere, its function being merely to echo the
argument to the print file at the time of a batch run.
Keyword errors
If an error occurs during keyword input, the program displays the offending line with ?
characters under the field that is causing difficulties.
The error numbers given by the program may be of two types:
• Up to and including 100.
These are errors detected within the PVTi data parser. Descriptions of these errors are given
in Table 6.17.
• Above 100
System data errors, usually caused by internal read operations used to convert a character
to its value.
PVTi keywords
This chapter contains details of all the keywords in PVTi. The keywords are listed in
alphabetical order.
• "Keywords A-D" on page 168
• "Keywords E-K" on page 191
• "Keywords L- O" on page 213
• "Keywords P- S" on page 241
• "Keywords T - Z" on page 267
First format
GROUP
COMB
SIMULATE
REGRESS
BLACKOIL • The keyword should be followed by up to ( N c -1) lines of data followed by a forward slash
PSEUCOMP
OUTECL3
character (/).
VFP • Each line should consist of three integers followed by the set of binary interaction
coefficients between two specified components, terminated with a forward slash character
(/). The first integer specifies the index of the first component and the second and third
integers the lower and upper indices of the other component.
Second format
The program uses this format when saving.
This format in effect displays the lower triangle of the BIC matrix, and there is no need for the
first three columns informing the program of the component and the lower and upper indices.
PVTi assumes the coefficients are read in as second, third component, etc., by row with the
columns as the first, etc., component.
Note For both formats the interaction coefficient between any given component and itself
must be zero. PVTi always ensures that this is the case.
Examples
First format
Binary interaction coefficients for a six component system.
BIC
2 1 1 -0.02000 /
3 1 2 0.10000 0.03600 /
4 1 2 0.13000 0.05000 /
5 1 2 0.13500 0.08000 /
6 1 5 0.12700 2* 2*0.0600 /
/
Here there are zero binary interaction coefficients between components 4 to 3, 5 to 3,4 and 6 to
2,3.
BIC
-0.02000
0.10000 0.03600
0.13000 0.05000
0.13500 0.08000
0.12700 2* 2*0.0600
/
Example
Characterizes two plus fractions, the first of which has known specific gravity and normal
boiling point temperature, the second of which has only molecular weight. Riazi-Daubert
critical properties and Edmister acentric factor correlations to be used:
CHARACT
C7+1 1* 0.7500 390.0 2* R E /
C7+2 200.0 /
/
If more than two of the four items, M w , S g , T b and K w are available, then the program uses the
first two, that is given data for M w , T b and K w , then PVTi uses M w and T b .
Example
For a nine-component condensate system, the first five being from library, the last four to be
user defined:
CNAMES
5* C4-6 C7+1 C7+2 C7+3 /
LNAMES
CO2 N2 C1 C2 C3 4* /
COMBINE
icl1 icu1 ifr1 zfr1 /
icl2 icu2 ifr2 zfr2 /
/
where:
• icl1 is a lower existing component index;
• icu1 is an upper existing component index;
• ifr1 is a group index;
• zfr1 is the mnemonic for a fraction.
Example
In a 19-component system, group the components 10 to 14 inclusive, and the components 15,
17 and 18. Accept default names for the fractions:
COMBINE
10 14 1 /
15 15 2 /
17 18 2 /
/
Example
OUTECL3
VFP
CORRACF
E /
Example
OUTECL3
VFP
CORRCP
C /
Example
Request debug from Newton and p sat routines:
CORRCP
C /
8 Observations
9 Saturation pressure/temperature
10 Flash routine
11 Newton routines
12 Fugacity evaluation (set as -1 to obtain Fugacities of components - batch
mode only)
13 Stability Check/Michelsen
14 Opening/closing of files
15 Regression output
16 LBC coefficients with regression
17 COMB routines
Example
Request debug from Flash and Michelsen routines only:
DEBUG
9* 10 2* 13 /
Example
Set to degrees Fahrenheit:
DEGREES
F /
Note
Example
For a two-component CO2/I-C4 system:
DREF
777.00 557.00 /
Example
EOS
PR3 /
Experiment Data
FLASH T P
BUBBLE T
DEW IRT T
CCE FLUID T P P2... ...PN
CVD FLUID T P1 P2... ...PN
DL T P P2... ...PN
SWELL FLUID YMF T [PRESS P] STYPE
X1... ...XN
SEPS T1 P1 L1 V1
(option 1) T2 P2 L2 V2
. . . .
TN PN LN VN
SEPS T1 P1 FL11 FV11 FV1N
(option 2) T2 P2 FL21 FV21...FL2N FV2N
. . . . .
TN PN LN FVN1...FLNN FVNN
COMPG FLUID T P D D1
TERN T P ZG1 ZG2 ZG3 IG1
IG2... ...IGN
VAPOR FLUID YMF T P M1... ...MN
MFLASH T P
CRIT
Experiment Data
TSAT P TSAT
FCMP YMF T FCMP
MCMP YMF ICORV T MCMP
PSAT T PSAT
MCT YMF ICORV T P FL2...FLN-1N MCT
(option 1) (option
1)
MCT YMF ICORV T P FV2...FVN-1N MCT
(option 2) (option
2)
TCLOUD T P
PCLOUD T P
where:
Keyword
Data requirements
arguments
T Single temperature entry
P Single pressure entry
IRT Retrograde (IRT=1) or normal (IRT=0) dew point
FLUID GAS or OIL system
PI I’th pressure
YMF Mnemonic of gas sample composition defined in SAMPLES or ZI.
STYPE Specifies subsequent data is MOLE fraction of gas in the mixture or GOR
gas volume at STP in the oil at original pressure or other pressure.
ICORV Condensing (ICORV=0) or vaporizing (ICORV=1) multiple contact
miscibility pressure required
XI Mole fraction or GOR (see STYPE) for step I.
MI Number of moles added cumulatively at step I.
TI Temperature of I’th separator stage.
LI Destination of liquid output from stage. Default is next stage, I+1.
VI Destination of vapor output from stage. Default is the cumulative vapor
output, 0.
Note, for separators, cumulative output stage is labelled 0, with vapor at
STC, liquid at final stage conditions.
FLIJ Fraction of liquid output from stage I to stage J. Default is 1.0 to the next
stage.
FVIJ Fraction of vapor output from stage I to stage J. Default is 1.0 to the
cumulative vapor output, 0.
D Reference depth.
Keyword
Data requirements
arguments
DI Depth I.
ZGI Names of the I’th ternary group.
IGI Ternary group index (1-3) for component I.
Example
Consider the following experiments:
1 Dew point pressure at 200 °F , no pressure estimate given;
2 Swelling test at 200 °F , injection gas composition specified in SAMPLES with mnemonic
ZINJ, four stage addition with given mole percentage of gas in reservoir fluid;
3 Constant composition expansion of condensate (GAS) at 200 °F , pressures in psi;
4 Constant volume depletion of condensate at 200 °F , pressures in psi;
5 2-stage separator at 500 psi, 125 °F and 50 psi, 90 °F with second stage fed from liquid
output of stage one;
6 Composition with depth: a volatile oil (OIL) at reference conditions of 220 °F , 4000 psi and
7500 ft. Require compositions and pressures at 7000, 7100,..., 8000 ft;
7 Ternary diagram calculation for an eight component system, explicit grouping;
8 Ternary diagram calculation for an eight component system, default grouping;
9 First contact miscibility pressure at 160 °F on ZI, injecting ZINJ;
10 Multiple contact miscibility pressure at 160 °F on ZI, injecting ZINJ, vaporizing drive;
11 Multiple contact miscibility pressure at 160 °F on ZI, injecting ZINJ, condensing drive;
12 Multiple contact test at 121.1 °F and 2738 psi on ZI, injecting ZINJ, condensing drive. At
each stage all of the remaining oil from the flash is contacted with one mole of ZINJ.
Note At pressures above the saturation pressure where no compositional changes take place,
the Constant Volume Depletion (CVD) and the Differential Liberation (DL) are
effectively a Constant Composition Expansion (CCE) experiment.
Restrictions
The EXP keyword arguments are subject to restrictions in the BLACKOIL, COMB, OUTECL3,
REGRESS and VFP sections:
Section Experiment
BLACKOIL CVD, DL & SEPS only.
PSEUCOMP CVD & SEPS only.
COMB CVD, SEPS & CCE only.
OUTECL3 CVD & COMPG only.
REGRESS Not COMPG, TERN, MFLASH, CRIT, TSAT, FCMP, MCMP.
VFP CCE & SEPS only.
BLACKOIL section
In the BLACKOIL section, to generate blackoil tables you must define a depletion experiment
and a separator network.
Number these experiments; 1 (CVD or DL) and 2 (SEPS).
OUTECL3 section
In the OUTECL3 section, you are restricted to a CVD for the purposes of generating either
KVTABLE or XMFVP and YMFVP tabular data.
Note Only two pressures should be entered, the maximum and minimum pressures required
in the table. The other pressures are calculated by interpolation (up to 50 pressures in
all).
You can also specify a COMPG experiment for the purpose of generating ZMFVD tables for
ECLIPSE Compositional.
Note Only two depths should be entered, the maximum and minimum depths required in the
table. The other depths are calculated by interpolation (up to 50 depths in all).
VFP section
In the VFP section, to generate blackoil tables you must define a depletion experiment and a
separator network.
Number these experiments 1 (CCE) and 2 (SEPS).
Example
X BLACKOIL
X PSEUCOMP
X OUTECL3
X VFP
We have 20 experiments in total in our project and wish that they all be ‘turned on’ during
regression except the 5th, 10th, 14th, 16th and 18th, which will be ‘turned off’:
OPTIONS
1 1 1 1 0 1 1 1 1 0 1 1 1 0 1 0 1 0 1 1/
FRAC
3 0.72 0.23 0.05 C7+1 C7+2 C7+3 /
Example
FVFREF
2 220.0 4000.0 /
4 200.0 0.0 /
7 210.0 /
/
Example
BLACKOIL
X PSEUCOMP
OUTECL3
VFP
Inject sample ZINJ at GI’s of 0.1, 0.2, 0.3 mscf/rb:
GI
ZINJ 0.1 0.2 0.3 /
By default, all components are assumed to be hydrocarbons, except for specific non-
hydrocarbons in the “LIB” set, that is, CO2 , H 2 S , H2 O , N2 , H2 and CO .
Example
An 11 component system consisting of CO2, N2, C1, C2, C3, IC4, NC4, IC5, NC5, C6, C7+:
HYDRO
N N 9*H /
Example
OUTECL3
VFP
LBCCOEF
1* 0.025 1* -0.04 0.01 /
Example
For a nine-component condensate system, with the first five from library and the last four user-
defined:
CNAMES
5* C4-6 C7+1 C7+2 C7+3 /
LNAMES
CO2 N2 C1 C2 C3 4* /
Example
MAXSTEP
0.2 /
MDP
1* 92.0 /
Example
MINDELP
14.7 /
Note This is not the minimum step that may be taken; it is the minimum upper limit imposed
on a step.
Example
MINSTEP
0.000001 /
• For the GOR option, the volume of the initial sample is usually the volume at its P sat at
the specified temperature. However, an alternative pressure can be specified; this is shown
in Example 2, below.
Each record must be terminated with a slash (/).
The set of records must end with a blank record, containing only a slash (/).
Examples
Example 1
Mix sample Z1 with Z2 to form Z3. The amount to mix is such that sample Z2 is 0.45 of the
mole fraction in the new sample Z3.
MIX
Z1 Z2 Z3 80.000 MOL 0.45 /
/
Example 2
Mix sample T1 with T2 to form T3. The amount to mix is represented as a GOR of 5.0 MSCF/
stb, and the volume of the first sample T1 is measured at a pressure of 5000 PSI at the mix
temperature of 670 °F . If the saturated volume of T1 was required then the PRES could be set
as 0.0, or those two items removed from the keyword.
MIX
T1 T2 T3 670.0 PRES 5000.0 GOR 5.0 /
/
MW
44.010 28.013 16.043 30.070 44.097 58.124
58.124 72.151 72.151 86.178 106.09 152.68
209.26 281.29 462.30 /
MWS
140.0 /
Example
BLACKOIL
PSEUCOMP
OUTECL3
VFP
Write out new PVI file with name SAVE:
NEWPVI
SAVE /
NEWPVO
SAVE /
OBS
no mnem weight data /
mnem weight data /
:
/
no mnem weight data /
mnem weight data /
:
/
/
where:
no specifies the experiment number given in the EXP keyword
mnem specifies the observation mnemonic (see below)
weight global weight to be given to subsequent data.
Several observations can be associated with a given experiment, each of which must be
terminated with a forward slash (/) character. The entry of observations for a particular
experiment is then terminated with a forward slash (/) character. A forward slash character (/)
also terminates the entry of all observations.
OBSIND
no mnem data /
mnem data /
:
/
no mnem data /
mnem data /
:
/
/
where:
no specifies the experiment number given in the EXP keyword ("EXP" on page 195);
mnem specifies the observation mnemonic (see OBS keyword, "OBS" on page 234);
data individual weights to be given to each point described in the previous OBS keyword.
Example
Observation weights for the experiments defined for the EXP keyword
(compare with example in OBS keyword description):
1 Saturation pressure and vapor Z -factor of dew point
2 Saturation pressure and swelling factor (relative volume) at each stage of swelling test
3 Relative volumes at each stage of CCE, vapor Z -factors at pressures above saturation
pressure only
4 Liquid saturation, vapor Z -factor and gas mole fractions at each stage of CVD (nine
component system)
5 GOR at each stage and stock tank GOR (in Mscf/stb) from two stage separators.
1 PSAT 50.0 /
ZV 1.0 /
/
2 PSAT 10.0 10.0 10.0 10.0 /
RELV 1.0 1.0 1.0 1.0 /
/
3 RELV 1.0 1.0 1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0 1.0 1.0
2.0 2.0 3.0 4.0 /
ZV 1.0 1.0 1.0 1.0 1.0 1.0 10* /
/
4 SL 1.0 1.0 1.0 1.0 /
ZV 1.0 1.0 1.0 1.0 /
YMF 1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0
0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
1.0 1.0 1.0 1.0 1.0 /
/
5 GOR 100.0 10.0 /
TGOR 10.0 /
/
/
Example
Ωa and Ωb values for two component system:
OMEGAA
0.457 0.456 /
OMEGAB
0.0780 0.0781 /
Option 11
Optional printing of results to the PVP file also takes specific arguments:
• 0 - ALways print;
• 1 - OPTionally print (that is user prompted);
• 2 - NEVer print.
Option 14
Alternative definitions of GOR in differential liberation experiments also takes specific
arguments:
• 0 - Default definition;
• 1 - Last stage removed and oil volume at its bubble point pressure;
• 2 - Incremental GOR;
• 3 - As default but volume of oil at its bubble point pressure rather than stock tank
conditions.
OPTIONS
0 0 1 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 /
Example
Parachors from [Ref. 18] for N2, C1, C2 and C3:
PARACHOR
41.0 77.0 108.0 150.3 /
PCRIT
72.90 36.00 /
Example
RETARG
0.1D-8 /
RTEMP
200.0 /
Example
Lean gas composition for swelling test:
SAMPLES
ZINJ 2* 0.9468 0.0527 0.0005 4* /
/
Example
Second sample from deeper in the column. Lean gas composition for swelling test:
SAMTITLE
Z2 ‘Second sample from deeper in the column’ /
ZINJ ‘Lean gas composition for swelling test’ /
/
Example
Save the total composition for experiment 4, a COMPG experiment, at stage (here depth) 4, and
label it Z80. Save the liquid output from stage 3 of experiment 2, a SEPS experiment, and label
it S23, and the stock tank vapor output of the same experiment, labelling it ST2:
SAVCOMP
4 ZMF 4 Z80 /
2 XMF 3 S23 /
2 YMF ST ST2 /
/
Example
The default and one alternative sample from different depths in the hydrocarbon column, and
therefore different mole weights and specific gravities for the C7+ plus fraction. Different
compositions are on the ZI and SAMPLES keyword, sample ZALT, which preceded this
keyword and data. Three pseudo-component split requested with minimum and maximum mole
weights of 90.0 and 300.0. Different skewness parameters reflects compositional grading in
fluid from top (ZI - low Alpha) to bottom (ZALT - high Alpha).
SCT
C7+ 3 90.0 300.0 /
ZI 1.10 130.0 0.75 /
ZALT 2.30 150.0 0.80 /
/
Example
For a two component system, in PVT-Metric units:
Example
For a two component CO2/I-C4 system, using degrees K :
TBOIL
194.70 261.30 /
TCRIT
305.6 408.1 /
THERMX
0.0005 /
Example
OUTECL3
X VFP
TLOW
40.0 /
Example
For a two component CO2/I-C4 system, using degrees K :
TLOW
40.0 /
Example
Set to field units:
UNITS
F /
1 The strings S1,S2... identify variables for regression, the indices I1, J1 (and K1, L1 for
binaries) identify component ranges, lower and upper.
Possible strings identifying variables (that is, S1, S2 etc.) are:
• TCRIT: Critical temperatures
• PCRIT : Critical pressures
• VCRIT: Critical volumes in LBC viscosity correlation
• ZCRIT: Critical Z -factors in LBC viscosity correlation
• ACF: Acentric factors
• OMEGAA: Ωa values
• OMEGAB: Ω b values
The quantities LL and UL specify limits imposed on the variables during the regression. All
variable limits are applied as a percentage of the current value, scaled to unity. Thus a maximum
value of 1.5 corresponds to allowing a 50% increase in any variable; a negative value implies
that the variable is allowed to change sign (restricted set!). Defaults are as follows:
Tc 0.5 1.5
Vc 0.5 1.5
Zc 0.5 1.5
Ωa 0.5 1.5
Ωb 0.5 1.5
ω 0.5 2.0
BIC -5.0 +5.0
si -5.0 5.0
Grouped regression variables can be defined by giving them the same variable index.
VAR
TCRIT 1 1 1 /
OMEGAA 2 3 1 /
BIC 1 1 10 10 1 0.30 2.00 /
SSHIFT 1 3 1 /
SSHIFT 4 4 2 /
SSHIFT 5 5 1 /
SSHIFT 6 10 2 /
/
Example 2
The default settings for the first five special regression variables with a system of three samples
is as follows:
VAR
’SCTALF’ 1 1 1 0.500000 5.000000 /
’SCTALF’ 2 2 2 0.500000 5.000000 /
’SCTALF’ 3 3 3 0.500000 5.000000 /
’SCTMW’ 1 1 1 0.900000 1.100000 /
’SCTMW’ 2 2 2 0.900000 1.100000 /
’SCTMW’ 3 3 3 0.900000 1.100000 /
’SCTKW’ 0.750000 1.250000 /
’CHEUHA’ 0.500000 1.500000 /
’THERMX’ -5.000000 5.000000 /
/
Note The special regression variables SCTALF and SCTMW cannot be grouped into one
variable and must be entered in the manner shown above, that is one variable for each
sample.
VCRIT values may be entered as an alternative to ZCRIT values - that is one of VCRIT or
ZCRIT should be entered. If VCRIT values are entered, critical Z -factors are obtained
internally using the usual relationship:
p crit V crit
-,
Z crit = ------------------ [EQ 7.1]
RT crit
Example
For a 20-component system, in field units:
VCRIT
1.5057 1.4417 1.5698 2.3707 3.2037
4.0847 4.9337 4.9817 5.8948 6.0710
7.8058 8.8374 9.8465 10.830 11.799
13.779 4.759 15.684 16.018 16.018
/
VCRITVIS values may be entered as an alternative to ZCRITVIS values; that is, one of
VCRITVIS or ZCRITVIS should be entered. If VCRITVIS values are entered, critical Z -
factors are obtained internally using the usual relationship:
p crit V critv
Z critv = --------------------- , [EQ 7.2]
RT crit
Note If VCRITVIS or ZCRITVIS is not entered then values entered with VCRIT or
ZCRIT are used in the viscosity correlation. The form of VCRITVIS is the same as
that of VCRIT.
Example
For a 20-component system, in Field units:
VCRITVIS
1.5057 1.4417 1.5698 2.3707 3.2037
4.0847 4.9337 4.9817 5.8948 6.0710
7.8058 8.8374 9.8465 10.830 11.799
13.779 14.759 15.684 16.018 16.018
/
Examples
Example 1
For PVTW output, salt present of salinity 10000 ppm, with gas dissolved. Temperature 200 °F ,
pressure 5514.7 psi.
Example 2
For PVTWSALT output, salt present with salinities 10, 20 and 30 gpl, with no gas dissolved.
Temperature 200 °F , pressure 5514.7 psi.
Examples
Example 1
For PVTW output, salt present of salinity 10000 ppm, with gas dissolved. Temperature 200 °F ,
pressure 5514.7 psi.
Example 2
For PVTWSALT output, salt present with salinities 10, 20 and 30 gpl, with no gas dissolved.
Temperature 200 °F , pressure 5514.7 psi.
Example
For PVTW output, salt present of salinity 10000 ppm, with gas dissolved. Temperature 200 °F ,
pressure 5514.7 psi.
Example
For PVTW output, salt present of salinity 10000 poem, with gas dissolved. Calculated at
temperatures 200 °F and 80 °F , pressure 6000.0 psi.
Example
Split C 7+ with z7+ = 0.1 , M 7+ = 140.0 , γ 7+ = 0.85 into three pseudo-components:
WHIT
7 0.10 140.0 0.85 3 /
Example
For a CO2/Methane/Ethane/I-C4 system:
ZCRIT
0.2709 0.2809 0.2808 0.279 /
Example
For a CO2/Methane/Ethane/I-C4 system:
ZCRIT
0.2709 0.2809 0.2808 0.279 /
ZI
0.6 0.4 /
Overview
This section of the manual contains information on the following:
• "Theoretical background of PVT" on page 295.
• "Analysis techniques" on page 371.
• "Regression in PVT analysis" on page 386
• "Wax and asphaltene precipitation in PVTi" on page 394
• "Cleaning samples contaminated with oil-based mud" on page 398
• "Mixing and recombination of samples" on page 400
Introduction
This section of the manual contains information on the following:
• "General background" on page 295.
• "Properties of pure components and mixtures" on page 296.
• "Characterization, splitting and grouping" on page 300.
• "Material balance" on page 316.
• "Flash calculations" on page 316.
• "Equation of state formulation" on page 317.
• "Multiphase flash" on page 323.
• "Viscosity correlations" on page 329.
• "Binary interaction coefficients" on page 336.
• "Blackoil systems" on page 338.
• "Gas condensate systems" on page 340.
• "Process simulation" on page 342.
• "Gas injection processes" on page 343.
• "Variation of composition with depth" on page 345.
• "Practical considerations" on page 347.
• "Theoretical model" on page 350.
• "Blackoil model" on page 353.
• "Pseudo-compositional tables for ECLIPSE GI option" on page 361.
• "Compositional data for ECLIPSE Compositional" on page 362.
• "Water properties" on page 362.
• "Model for API Tracking option in ECLIPSE BlackOil" on page 364.
• "Compositional Data for ECLIPSE Thermal" on page 367.
General background
Both in the fully compositional and extended black oil formulations of reservoir simulation, an
accurate description of the hydrocarbon system and its properties is important.
In the compositional case the phase behavior of the hydrocarbon system is an integral part of
the simulation. Equal phase fugacity conditions or flash calculations are used to determine the
phase split and composition.
In the black-oil and extended black-oil formulations averaged properties of the oil and gas
phases such as Bo (oil volume factor), Rs (solution gas-oil ratio), B g (gas volume factor) and R v
(condensate gas ratio) are obtained through laboratory experiments or by using an equation of
state description of the hydrocarbon system.
Pure components
The natural gas and crude oil found in underground natural deposits are mixtures of a large
number of different hydrocarbon components and some additional non-hydrocarbons such as
CO2 , H 2 S , N 2 , H 2 and CO .
Alkanes or paraffins
Alkanes or paraffins are very stable chemical compounds. The chemical formula for compounds
in the group is C n H 2n + 2 . The group contains hydrocarbon components such as:
Formula Name
C1 H4 Methane
C2 H6 Ethane
C3 H8 Propane
C 4 H 10 Butane
C 5 H 12 Pentane
C 6 H 14 Hexane
C 7 H 16 Heptane
C 10 H 22 Decane
Each of these molecules can exist as a single straight chain of carbon atoms attached to three
(end carbons), or two (non-end carbons) hydrogen atoms (except Methane, which has four
attached hydrogens).
Isomers
For molecules containing four carbons or more, additional structures other than the straight
chains are possible. The number of these branched molecules or isomers increases with
increasing carbon number.
• Butane has one isomer in addition to the straight-chained normal-Butane, the branched
molecule iso-Butane.
• Pentane has two isomers in addition to the normal-Pentane, known as iso- and neo-Pentane.
The neo- isomer (cross-shaped) is rarely found in naturally-occurring petroleum. Most
fluid analyses include the iso- and normal- components of Butane and Pentane.
• Hexane (six carbon atoms) has four isomers in addition to the normal-Hexane.
Note The number of isomers increases rapidly with carbon number making the
identification of the increasingly small concentration of these isomers impossible.
Formula Name
C3 H6 Cyclopropane
C 5 H 10 Cyclopentane
They share the same formula as the alkenes, which are the same as the alkanes except that one
or more pairs of carbon atoms are linked by double bonds with the consequent loss of one or
more hydrogen atoms.
Note Alkenes and their triple bond equivalents, the alkynes, are rarely found in naturally
occurring petroleum.
Any one of the hydrogens can be replaced by a methyl-( CH3 ) or higher group, for example
C 7 H 14 , Methylcyclohexane.
Aromatics
The third family of hydrocarbons are aromatics. These are ring-type structures. Some of the best
known members of this family are:
Formula Name
C6 H6 Benzene
C7 H8 Toluene
SCN groups
All of napthenes, aromatics, more complicated hydrocarbons (polynuclear molecules),
hydrocarbons compounds containing other species, namely Nitrogen, Sulphur, Oxygen and
certain trace metals, are all capable of existing as isomers. Identification of even a relatively
small number of the possible isomers within a petroleum mixture is a complicated and therefore
expensive task.
It is standard practice within the petroleum industry to lump all isomers together on the basis of
the boiling point of the molecule.
Single Carbon Number (SCN) groups are defined for components, for example hexanes,
heptanes, etc. Components are pure hydrocarbon components with normal boiling point
temperatures, that is at one atmosphere pressure, between consecutive normal-paraffin boiling
points.
For example, the SCN hexanes group consists of those hydrocarbons that boil between the
normal boiling points of n-pentane and n-hexane.
Term Nomenclature
Tc Critical temperature
pc Critical pressure
Vc Critical volume
Zc Critical Z -factor
Mw Molecular weight
ω Acentric factor
Tb Normal boiling point
The component library in PVTi contains properties for some of the more common pure
components.
Multi-component mixtures
As mentioned previously, natural gas and crude oil contains literally thousands of different
components. No attempt is usually made to identify all the hydrocarbons beyond C 5 : rather,
SCN groups denoted hexanes, heptanes, etc., are used. Even so, you may be left with ten or more
components in any given petroleum mixture.
To describe a multi-component system thermodynamically using an equation of state model, the
system must be defined in terms of the properties of the components and their mole fractions in
the mixture.
The terms and nomenclature used to characterize multi-component mixtures may be
summarized as follows, where the subscript i denotes the i th component:
Term Nomenclature
n Number of components in the mixture
T c, i Critical temperature
p c, i Critical pressure
M w, i Molecular weight
ωi Acentric factor
Term Nomenclature
xi Mole fraction in the liquid phase
Phase diagrams
In terms of a phase diagram plotted against pressure and temperature, the upper part of the phase
envelope, up to the critical point, represents the bubble-point curve.
From the critical point the phase envelope continues as dew point curve. The phase envelope
encircles the two phase region in the phase diagram.
A phase diagram is often characterized by the maximum pressure, the cricondenbar, and the
maximum temperature, cricondentherm.
The area within the phase envelope bounded by a vertical line through the critical point and the
line represented by the cricondentherm is the retrograde area. A reduction of pressure into this
area results in retrograde condensation.
For the RK and ZJ equations of state the critical point calculation implements the theory
developed in [Ref. 50].
For the PR and SRK equations of state the calculations use the theory developed in the follow-
up paper in [Ref. 47].
Components
All components in PVTi are labeled as being one of three possible types:
• LIB: Library
LIB components
Lib components are the commonly occurring hydrocarbons C 1 , C 2 ,..., C 45 and the specific non-
hydrocarbons H2 O , H 2 S , CO2 , N 2 , H 2 and CO .
CHAR components
Components designated as CHAR are typically the last component in a PVT report and as such
comprise the plus fraction, that is, contains all remaining hydrocarbons. For such fractions, a
laboratory would typically measure only the mole weight and specific gravity, denoted here by
M N+ , γ N+ , where N is the number of the first carbon group in the fraction.
User components
User components are those for which you must supply all the data necessary to define the
equation of state parameters, that is T c , pc , Vc , M w and ω .
Characterization
Clearly, correlations are required that are capable of generating the critical properties and
acentric factors from just the mole weight and specific gravity.
In general, any p - V - T report gives at least the mole weight of the plus fraction. If this is the
only information available, PVTi estimates a value of the specific gravity from known (SCN)
distributions. Otherwise, any two out of the four of M w , γ , T b , K w is sufficient to characterize
a component.
PVTi employs five sets of correlations for generating critical properties by:
• Kesler and Lee [Ref. 10]
• Cavett [Ref. 11]
• Riazi and Daubert [Ref. 12]
• Winn [Ref. 42] and [Ref. 43]
• Pedersen [Ref. 41].
It also uses four sets of correlations for acentric factors by:
• Lee and Kesler [Ref. 13]
• Edmister [Ref. 14]
• Thomassen [Ref. 30]
• Pedersen [Ref. 41].
Splitting
The plus fraction often has an importance that appears to far outweigh its relatively small mole
fraction of a fluid sample. In particular, saturation pressure calculations can be extremely
sensitive to the mole fraction and properties of the plus fraction. More accurate predictions
requiring less regression of equation of state parameters can be achieved if a thorough
description of the plus fraction can be made.
Ideally, a complete true-boiling-point distillation (TBP) should be made that yields not only the
detailed composition of the plus fraction but also the boiling points, specific gravities and
molecular weights of the constituent components. However, this analysis is rarely performed
and so a general procedure to describe the distribution of components and properties within a
plus fraction is required.
Three techniques are available within PVTi for the splitting of the plus fraction into sub-
fractions:
• Constant Mole Fraction (CMF);
• Whitson.
• Semi-Continuous Thermodynamics (Modified Whitson)
All of these techniques rely on a probability density function (PDF) to relate mole fraction to
mole weight. This approach is due to Whitson [Ref. 4].
where:
• α, β, η are parameters defining the distribution
• Γ is the gamma-function
• α gives a measure of the shape of the distribution
• η is the lowest mole weight in the plus fraction
• β is a normalization condition that can be determined from the condition:
M N + – η = αβ [EQ 8.3]
where M N + is the average mole weight of the plus fraction.
The cumulative probability function, P ( X ≤ x ) is the integral of p ( x ) from η to x :
x
P(X ≤ x) = ∫η p ( x ) dx [EQ 8.4]
where y = ( x – η ) ⁄ β
The frequency, f i of a component i having mole weight boundaries M i – 1 , M i is given by the
integral:
Mi
fi = ∫M i–1
p ( M ) dM [EQ 8.6]
= P ( M ≤ M i ) – P ( M ≤ Mi – 1 ) [EQ 8.7]
and the mole fraction zi is related to its frequency by:
zi = zN + fi [EQ 8.8]
The average mole weight in the same interval is given by:
P ( M ≤ M i ,α + 1 ) – P ( M ≤ M i – 1 ,α + 1 )
M i = η + αβ -----------------------------------------------------------------------------------------------
P ( M ≤ M i ,α ) – P ( M ≤ M i – 1 ,α )
- [EQ 8.9]
where the P ( X ≤ x ) functions all use the same value of β regardless of the value of α , (that is α
or α + 1 ) and where z N + is the total mole fraction of the plus fraction.
In the absence of any other data, you should assume that α is unity. The parameter η is the
minimum mole weight that occurs in the plus fraction, therefore if the plus fraction were C 7+ ,
then a good estimate to η is 92.0.
Generally, Whitson recommends:
η = 14n – 6 [EQ 8.10]
where n is the first SCN group in the plus fraction, for example 7.
Note Both splitting techniques available in PVTi use the Whitson PDF to define the mole
fraction/mole weight distribution of the plus fraction, however, they differ slightly in
the way they sub-divide into the required number of pseudo-components.
Generally, you may consider skewing the distribution thus: highest mole fraction for the lightest
pseudo-component to smallest mole fraction for the heaviest pseudo-component.
Hint If you have evidence on the shape of the plus fraction distribution, it may be prudent
to manually change the skewness parameter α , that is less than unity for a condensate
and greater than unity for a crude oil.
Whitson
The Whitson technique consists of splitting the plus fraction into SCN groups N (the first in the
plus fraction) to some high number, for example 45.
The mole weight boundaries are estimated from the user specified plus fraction mole weight and
specific gravity from which one can calculate a plus fraction Watson factor:
0.15178 – 0.84573
K w, N+ = 4.5579M w, N+ γ N+ [EQ 8.11]
Assuming that an average Watson factor K w, N+ can be applied to each individual SCN groups,
the SCN groups’ specific gravities can be estimated from equation (1), namely:
1⁄3
T b, i
γ i = --------------
- [EQ 8.12]
K w, N+
where the T b, i are the normal boiling point temperatures of the SCN groups, which in fact define
the SCN groups. That is SCN C6 is all hydrocarbons that have normal boiling point
temperatures between those of normal- C 5 and normal- C 6 , etc.
Then, with known boiling points and specific gravities, mole weights can be estimated using the
Riazi-Daubert correlation, which is of the general form:
b c
θ = aT b γ [EQ 8.13]
where θ is the property to be evaluated, say M w and ( a ,b ,c ) are tabulated coefficients, see [Ref.
4].
With average mole weights for the SCN groups, the boundary mole weights are estimated from:
1
M i = --- ( M i – 1 + M i )
2
[EQ 8.14]
These can then be used to integrate the PDF to give mole fractions for each of the SCN groups.
Hint Generally, you will not want to work with 20-30 additional components. Therefore the
SCN groups are pseudoised down into Multi-Carbon Number (MCN) groups. A
method for estimating the required number of MCN groups is that due to Sturge’s,
discussed in [Ref. 4]. However, this often predicts the use of four or five pseudo-
components. It has been our experience that two or three pseudo-components is
sufficient for most purposes.
where the summation is performed over the SCN groups i in the MCN group l and z l is the
mole fraction of the l th MCN group.
N
∞
∫0 ∑ wi f ( xi )
–x
f ( x )e dx = , [EQ 8.16]
i=1
where the weighting factors w i and quadrature points x i are determined from a class of Laguerre
polynomials for a given order N , see [Ref. 31].
For our mole weight/mole composition model, we can associate the quadrature points xi with
pseudo-component mole weights M i = η + βx i and the mole fractions z i = w i f ( xi ) .
To account for compositional variations with depth, Whitson extended the standard model to
account for these variations, or plus fraction mole weight and skewness parameter. Since
different plus fraction mole weights and skewness parameters lead to different β parameters,
the following modification was introduced:
β0
β = --------------------------
- [EQ 8.17]
[ 1 + ln ( δ ) ]
where β 0 and δ are parameters in the modified distribution function p0 ( M ) given by:
(α – 1) α
(M – η) – ( M – η ) ( 1 + ln ( δ ) ) -
p 0 ( M ) = --------------------------------
α
exp ---------------------- -----------------------------
β0 [ ( M – η ) ⁄ β0 ]
= p(M) [EQ 8.18]
β0 Γ ( α ) δ
x0 = M – η-
------------- [EQ 8.19]
β0
with dM = β 0 dx to give:
∞ ( α – 1 ) –x0
x0 e
∫0 ---------------------------------------------------
x
Γ ( α )δ 0 ( 1 + ln ( δ ) )
–α 0
dx = 1 [EQ 8.20]
and thus the mole fraction and mole weight are given by:
Procedure
The procedure for using this method is particularly attractive.
Firstly we choose the number of pseudo-components required, N , which is normally in the
range 2 ≤ N ≤ 5 : this then fixes the values of weights and points, w i and xi .
The minimum mole weight η is chosen as before.
Next we specify β 0 which is done by selecting the maximum mole weight component to be
used: Whitson recommends M N = 500.0 but we prefer 2 × M N+ . Either way we now have:
MN – η
β 0 = ----------------- [EQ 8.24]
x 0N
Next we estimate a value for the skewness parameter, which as before is taken as unity in the
absence of any other information. Then to satisfy the total number of moles and mass, we
calculate:
β
δ = exp ⎛⎝ ----0- – 1⎞⎠ [EQ 8.25]
β
and calculate the mole fractions and mole weights. Finally, the average plus fraction mole
weight is checked:
∑ zi Mi
M N+ = i-------------------
=1
N
- [EQ 8.26]
∑ zi
i=1
and the δ is adjusted if [EQ 8.26]does not equate to the measured value.
The appeal of this model is it allows a variable plus fractions mole weight and skewness
parameter by sample composition. That is different plus fraction mole weight and mole fraction,
but allows us to use a common set of pseudo-components with fixed mole weights and hence
critical properties.
Hint Having used this procedure to characterize the plus fraction(s), you can then use the
special regression facility, which allows the sample plus fraction mole weights and
skewness parameters to be variables.
In particular, where there is a known variation of composition with depth and you have
at least two samples from different depths in the hydrocarbon column, this technique
can be used with some success, see [Ref. 32].
where Fcal
c is the value calculated by the Whitson model, that is assuming SCN-cuts, and F c
reg
is
the potentially regressable function which defaults to unity. This allows you to make your fluid
more or less aromatic with respect to the standard distribution (equivalent to a K w ∼ 12.0 ), which
in turn yields different critical properties, etc., and hence different fluid behavior.
Grouping
Grouping of components is performed by one of three techniques:
• Molar averaging
• Weight averaging
• Mixing rule.
All of these are explained in the comprehensive text of Joergensen and Stenby [Ref. 44].
Consistency tests
Most fluids that show some form of compositional behavior, namely gas condensates and
volatile oils, are subject to a constant volume depletion experiment (CVD) as part of their
analysis. For further information see "Compositional material balance" on page 307.
A material balance calculation can be performed using the data that is generated from such an
experiment from which quantities such as liquid compositions, K -values, molecular weights
and densities of vapor and liquid, etc., can be evaluated. This data can then be examined to look
for any inconsistencies, such as lack of monotonicity.
Additional data often given in a laboratory PVT report are separator compositions for the liquid
and vapor streams. If such data is available, then K -values can be constructed and plotted
against theoretical fits, again as a test of consistency.
Finally, estimates of recovery of gas and oil from the reservoir can be made as well as generation
of blackoil tables, without recourse to the use of the equation of state model.
Note The model employed in PVTi was adapted from a program developed by Pearce [Ref.
20].
Equil
Equilibrium vapor Pressures (psi)
Liquid
Component 6764.7 5514.7 4314.7 3114.7 2114.7 1214.7 714.7 714.7
Carbon Dioxide 2.37 2.40 2.45 2.50 2.53 2.57 2.60 .59
Nitrogen 0.31 0.32 0.33 0.34 0.34 0.34 0.33 .01
Methane 73.19 75.56 77.89 79.33 79.62 78.9 77.8 12.42
Ethane 7.80 7.83 7.87 7.92 8.04 8.40 8.70 3.36
Propane 3.55 3.47 3.40 3.41 3.53 3.74 3.91 2.92
isoButane .71 .67 .65 .64 .66 .72 .78 .91
n-Butane 1.45 1.37 1.31 1.30 1.33 1.44 1.56 2.09
isoPentane .64 .59 .55 .53 .54 .59 .64 1.4
n-Pentane .68 .62 .58 .56 .57 .61 .66 1.6
Hexanes 1.09 .97 .88 .83 .82 .85 .9 3.68
Heptanes plus 8.21 6.2 4.09 2.64 2.02 1.84 2.12 71.01
Total 100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00
Mw(C7+) 184.0 160.0 142.0 127.0 119.0 115.0 114.0 213.0
Sg(C7+) .816 .799 .783 .770 .762 .758 .757 .833
Zvap 1.238 1.089 .972 .913 .914 .937 .960
Nprod 0.000 9.024 21.744 38.674 55.686 72.146 81.301
Sliq 0.00 14.1 19.70 21.60 21.30 20.20 19.3
This gas condensate with an assumed reservoir composition equal to composition at the first
pressure, which corresponds to the saturation pressure of the fluid, psat = 6764.7 psia is put
through the stages of the CVD.
The CVD consists of reducing the pressure to the required stage pressure, and then removing
gas so that the volume remaining in the cell is equal to volume at the saturation pressure.
where n pk is the number of moles of gas produced at stage, denoted Nprod in the table above,
so that Δn pk is the incremental gas recovered during the kth stage, that is, reducing the pressure
from p k – 1 to p k .
Similarly, the mole fraction of the i th component remaining in the cell after the j th stage is:
Cell volumes
To convert moles and mole fractions to volumes, one must have some reference volume, such
as the volume of the cell.
For a gas condensate, the cell volume, for one mole, can be calculated from the real gas law:
Z 1 RT
V cell = ------------
- [EQ 8.32]
p1
where:
• Z1 and p1 are the vapor Z -factor and pressure at the first (saturation) pressure
• R is the universal gas constant
• T is the temperature at which the CVD is performed, assumed constant throughout the
experiment.
Knowing the cell volume and the liquid saturation at the j th stage, S lj (denoted Sliq in the table
above), the volume of the cell occupied by liquid is:
V lj = S lj V cell [EQ 8.34]
and thus the volume of the cell occupied by gas is:
V vj = ( 1 – S lj )V cell [EQ 8.35]
Thus, the composition of the liquid remaining in the cell at the j th stage can be determined from:
n tj z ij – n vj y ij
x ij = ------------------------------
- [EQ 8.37]
n tj – n vj
Equilibrium K-values
The K -values can now be estimated from the ratio of the vapor and liquid mole fractions:
y
K ij = ----ij- [EQ 8.38]
x ij
The best test of the consistency of the CVD data is to plot the logarithm of the p j K ij product
against the characterization factor F , [Ref. 20], where:
F = b i ⎛⎝ ------
1 - – --1-⎞ [EQ 8.39]
T bi T⎠
Here:
• T ci and T bi are the critical and normal boiling point temperatures
• p ci is the critical pressure (in psia)
• p ref is the reference pressure, usually 14.7 psia.
Plots
As pointed out by Bashbush, [Ref. 19], the curves should plot in a parallel-like trend, that is,
there should be no humps in any individual curve or crossing of any pair or pairs of curves.
The highest curve should belong to nitrogen, followed by methane and carbon dioxide.
Physical properties
As well as performing material balances on moles and mole fractions, material balance can be
done on the total, liquid and vapor masses:
m tj = m lj + m vj [EQ 8.41]
The total mass remaining the cell at the j th stage is given by:
where M s and M vk are the initial or saturation pressure molecular weight and the vapor
molecular weight at the kth stage, respectively, both of which can be evaluated from Kay’s rule:
Nc
Correlations
Alani-Kennedy liquid densities
The liquid molar volume of a pure hydrocarbon component is given by the smallest real positive
root of the Van der Waal’s equation:
RT = ⎛⎝ p + -----
a-⎞ ( V – b )
2⎠
[EQ 8.46]
V
C N+ C N+
Note The non-hydrocarbons N 2 , CO2 and H 2 S are added into the C 1 fraction.
Vapor Z-factor
The vapor Z -factor, Z vap , can be estimated from the correlation due to Dranchuk, Purvis and
Robinson [Ref. 21]:
0.27p
Z vap = ---------------r [EQ 8.50]
ρr Tr
where T r and p r are the reduced temperature ( = T ⁄ T c ) and reduced pressure ( = p ⁄ pc ) , and ρ r
is the pseudoreduced density, which is found by iteratively solving:
f ( ρr ) = 0 [EQ 8.51]
where:
6 3 2 3 2 2
f ( ρ r ) = Aρ r + Bρ r + Cρ r + Dρ r + Eρ r ( 1 + Fρ r ) exp ( – F ρ r ) – G [EQ 8.52]
and:
A = 0.06423 [EQ 8.53]
B = 0.5353T r – 0.6123 [EQ 8.54]
D = Tr [EQ 8.56]
0.6816-
E = --------------- [EQ 8.57]
2
Tr
where γ HC
g is the hydrocarbon gas gravity that is related to the mixture gas gravity, γ mix
g by:
mix
HC γ g – 0.9672y N – 1.5195y CO – 1.1765y H S
γg 2 2 2
= ------------------------------------------------------------------------------------------------------------ [EQ 8.63]
1 – y N – y CO – y H S
2 2 2
with:
mix
γg
∑ yi Mi
= ------------------- [EQ 8.64]
M air
and
HC
p pc = ( 1 – y N – y CO – y H S )p c + 493.0y N + 1071.0y CO + 1306.0y H
2 2 2 2 2 2S
[EQ 8.66]
The correction to T pc and ppc in the presence of sour gas ( CO2 and H2 S ), due to Wichert and
Aziz, is:
to give:
T∗ pc = T pc – ε [EQ 8.68]
and
p pc ( T pc – ε )
p∗ pc = ----------------------------------------------------
- [EQ 8.69]
T pc + y H S ( 1 – y H S )ε
2 2
B = 3.5 + 986.0
------------- + 0.01M w [EQ 8.72]
T
with:
–4 –8 2
a = 1.20 + 4.5 ×10 p + 15.0 ×10 p [EQ 8.76]
–4 –8 2
c = 0.89 – 1.7 ×10 p – 3.5 ×10 p [EQ 8.77]
where F is the Hoffmann et al. characterization factor, refer to [EQ 8.39].
As most separators are operated at pressures less than 1000 psia, the so-called Standing K -
values provide a reasonable test of the quality of the calculated separator data and can be used
as an alternative for blackoil table generation etc., when no such actual data is available.
• where:
• G is the GOR (gas at STC, oil at stage conditions)
Recovery calculations
For reservoirs whose initial pressure, p init , is greater than its saturation pressure, psat , estimates
of fluid recovery can be made if the Z -factors in the range p init ≥ p ≥ psat are known.
Note Above the saturation pressure, no compositional changes in the fluid take place so for
a unit volume of the fluid at isothermal conditions:
n1 Z1 n2 Z2
----------
- = ----------
- [EQ 8.79]
p1 p2
Starting with 100 moles of fluid at the initial reservoir pressure, the number of moles at some
pressure p j above the saturation pressure is:
left rem
nj = 100 – n j [EQ 8.80]
where:
p Z
= 100 ⎛⎝ ----j⎞⎠ ⎛⎝ -----1⎞⎠
rem
nj
Zj p1
[EQ 8.81]
At pressures below the saturation pressure, the produced moles are known from the CVD
experiment.
The fluid at each stage can then be flashed to surface conditions, using separator or Standing’s
K -values to get the volumes of liquid and vapor that could be expected at the surface, assuming
100% efficiency in the separators.
An option exists to produce liquid along with gas, in a ratio defined by input relative
permeabilities, using the technique described in Reudelhuber and Hinds, refer to [Ref. 37].
General information
PVTi allows you to fit data to an Equation of State. See "The fluid model" on page 98 and
"Batch system and keywords" on page 152 for further information on how to define Equations
of State in PVTi.
This section of the manual contains information on:
• "Material balance" on page 316.
• "Flash calculations" on page 316.
• "Equation of state formulation" on page 317.
• "Surface tensions" on page 320.
• "Three-parameter equation of state" on page 320.
• "Binary interaction coefficients" on page 336.
Material balance
We consider a hydrocarbon system consisting of one mole of mixture or feed of composition z i .
In general, this consists of:
• L moles of liquid of composition xi
• V moles of vapor of composition y i
• at some pressure and temperature, p and T .
The conditions exist that:
L+V = 1 [EQ 8.82]
∑ xi = ∑ yi = ∑ zi = 1 [EQ 8.83]
i i i
Flash calculations
For a thermodynamic system to be in equilibrium it is further required that the fugacities in the
liquid and vapor phases must be equal for each component:
f iL = f iV [EQ 8.85]
The equations [EQ 8.82] to [EQ 8.85] in the unknowns L , V , x i and yi may be solved to find
the equilibrium state of a system of known composition z i .
The fugacities are functions of temperature, pressure and composition (see "Equation of state
formulation" on page 317),
f i = f i ( T , p, x i)
the mole fractions of each component in the liquid and vapor phases are defined as:
zi
x i = ------------------------------------
- [EQ 8.86]
[ 1 + ( K i – 1 )V ]
and
Ki zi
y i = ------------------------------------
- [EQ 8.87]
[ 1 + ( K i – 1 )V ]
The flash calculation determining the equilibrium conditions in the two-phase region is
performed in two stages:
Use Michelsen’s stability criterion, [Ref. 16], to establish how many phases are present. The
stationary condition used to determine the Gibb’s energy minimum is obtained using successive
substitution accelerated with the general dominant eigenvalue method, [Ref. 8], switched to
BFGS minimization near the minimum, [Ref. 9].
If two phases are present, solve the phase split calculation to obtain equal fugacity liquid and
vapor states. The solution variables used are Ji = log ( K i ) and the vapor mole fraction. The
stability check provides an initial set of values, which are refined using accelerated successive
substitutions, and finally solved using a full Newton-Raphson method.
In the case of finding a saturation pressure, that is a bubble point for an oil or a dew point for a
condensate, the same N -fugacity equations are used, [EQ 8.85], and a constraint based upon the
difference in the sum of vapor and liquid mole fractions, [EQ 8.83]. Rather than iterating to find
the vapor fraction, V , or liquid fraction L = 1 – V , from [Ref. 1], V is set to 0 (zero - bubble
point) or 1 (one - dew point) and we vary pressure until the set of conditions is satisfied.
The cubic equation for the Z -factors may be solved to obtain Z -factors for liquid and vapor
phases. Generally three solutions are obtained. The distinction between the liquid and vapor
phase is then made by choosing the smallest root as the Z -factor for the liquid phase and the
largest root as the Z -factor for the vapor phase.
Fugacity coefficients are calculated using:
f A 2S B ( Z + m2 B ) B
ln ⎛⎝ ------i -⎞⎠ = – ln ( Z – B ) + ---------------------------- -------i – -----i ln ------------------------
- + -----i ( Z – 1 ) [EQ 8.93]
px i ( m 1 – m 2 )B A B ( Z + m1 B ) B
where
n n
B = ∑ ( xj Bj ) [EQ 8.96]
j=1
1---
2
A jk = ( 1 – k jk ) ( A j A k ) [EQ 8.97]
and kjk are binary interaction coefficients, normally between hydrocarbons and non-
hydrocarbons.
These four equations, [EQ 8.94] to [EQ 8.97], express the mixing laws used in all the equations
of state.
The variables A j and B j are defined by the following equations:
p rj
A j = Ω A ( T, j ) ------
2
[EQ 8.98]
T rj
and ΩB ( T, j ) used in equations [EQ 8.98] and [EQ 8.99] are functions of the acentric
Ω A ( T, j )
factor ω and the reduced temperature T rj .
For Redlich-Kwong
1-
– --
2
Ω A ( T, j ) = Ω A T rj
0
[EQ 8.100]
Ω B ( T, j ) = Ω B
0
[EQ 8.101]
For Soave-Redlich-Kwong
1 2
2 ⎛
--- ⎞
2
Ω A ( T, j ) = Ω A 1 + ( 0.48 + 1.57ω j – 0.17ω j ) ⎜ 1 – T rj⎟ [EQ 8.102]
0
⎝ ⎠
Ω B ( T, j ) = Ω B
0
. [EQ 8.103]
For Zudkevitch-Joffe
– 1---
2
Ω A ( T, j ) = Ω A F aj ( T )T rj
0
[EQ 8.104]
Ω B ( T, j ) = Ω B F bj ( T )
0
[EQ 8.105]
For Peng-Robinson
1 2
2 ⎛
--- ⎞
2
Ω A ( T, j ) = Ω A 1 + ( 0.37464 + 1.54226ω j – 0.2669ω j ) ⎜ 1 – T rj⎟ [EQ 8.106]
0
⎝ ⎠
Ω B ( T, j ) = Ω B
0
[EQ 8.107]
The normal PR form can be optionally modified for large acentric factor, using the factor
2 3
rather than ( 0.37464 + 1.54226ωj – 0.2699ω 2j ) for
( 0.379642 + 1.48503ω j – 0.164423ω j + 0.016666ω j )
ω j > 0.49 . This correction is invoked by use of the PRCORR keyword or interactively from the
Equation of State panel. For further information see "The fluid model" on page 98.
ΩA
0
and Ω B0 are constants depending upon the equation of state:
Equation ΩA ΩB
0 0
Zudkevitch-Joffe equation
The Zudkevitch-Joffe equation contains additional temperature factors, denoted F aj ( T ) and
F bj ( T ) , multiplying the usual RK temperature dependence. These are adjusted to match the pure-
component fugacity values along the vapor pressure line, and to observe the correct component
liquid density.
In the program the required variation of component saturation pressure and liquid density with
temperature are obtained using correlations of Reidel, and Gunn and Yamada. These
correlations require the input of the normal boiling point and the reference liquid density at a
specified reference temperature.
Surface tensions
The surface tension between the liquid and vapor phase of a multi-component mixture can be
estimated by the Macleod-Sugden relationship, [Ref. 18]:
Nc
1---
∑ [ Pi ] ( ρm
4 Liq Vap
σ mix = xi – ρm yi ) [EQ 8.108]
i=1
where [ P i ] is the parachor of the i th component, see [Ref. 18], which has a liquid and vapor mole
fraction of xi and yi respectively, and the liquid and vapor molar densities are ρ Liq Vap
m and ρ m
respectively. If the parachors are assigned values consistent with [Ref. 18] and the molar
volumes are expressed in gmoles/ cm3 then the surface tension σ have units of dynes/cm.
Note In PVTi, this correction is available for the Peng-Robinson and the Soave-Redlich-
Kwong Equation of State.
∑ zi ci
EoS
V mol,p = V mol,p – [EQ 8.109]
i=1
where:
RT ci
-.
b i = Ω b, i ---------- [EQ 8.111]
p ci
Shift parameters
where, for each component, s i is the dimensionless Volume Shift, ω i the acentric factor and,
as mentioned above, T r =0.7.
For heavier components ( M w > 90 ) the shifts are estimated on a component by component basis
as the difference in the calculated single component molar volume using the two-parameter
Equation of State at the components reference pressure and temperature (often 60 °F and 14.7
psia), and the reference molar volume (mole weight divided by reference density).
The shift parameters of the heavier components are potential regression variables and are added
to the allowed set if you request either the PR3 or the SRK3 Equation of State. Experience has
shown that simultaneous matches to saturation pressure and density can be achieved more easily
as a result, and therefore the use of these Equations of State in PVTi is recommended.
κis a special regression function in PVTi, and takes the default value of 0.0005 °C –1 ,
determined from a fit to a crude oil. Physically, it represents a thermal expansion factor, and has
keyword THERMX within PVTi.
Polynomial Correlations
If the temperature dependence of the Volume Shifts is set to Polynomial Correlations then the
equations of Soreide are used which were introduced in the case of Tr<1 in equations [EQ 8.109]
and [EQ 8.109]. When the None option was set we always set Tr=0.7 and so these equations
were fine. However, if a temperature above Tc is supplied to the algorithm with the Polynomial
Correlations option set so that Tr>1, then PVTi chooses from a set of equations depending on
the mole weight of the component. These equations can be found in [Ref. 34].
For heavier components ( M w > 90 ) the thermal expansion scheme outlined in the "Linear
Expansion Only" on page 322 is applied. However, since the scheme is now applied to regions
where the fitted functional κ i is no longer valid, κ i is set to unity, and the thermal expansion
factor, κ , applies directly. For this reason it is often a good choice as a regression variable for
heavy hydrocarbon components if one is finding difficulty matching liquid densities.
The Linear Expansion Only model has some advantages in that since the modification to
volume shifts is linear, the volume shifts may be regressed for all components.
For the Polynomial Correlations model, since the correlations of Søreide are not linear, you
cannot regress the shift parameters of lighter components to which the correlations have been
applied.
Note If the volume shifts are Dependent, then they are calculated so that the calculated
reference density (the density at the reference temperature and standard conditions
pressure) is guaranteed to match the observed value. This means that the volume shifts
are dependent on the critical properties and acentric factor of the component. If the
volume shifts are Independent this guarantee is broken and the volume shifts can be
varied during regression. In this situation they are not affected by changes in the
critical properties or acentric factors.
Multiphase flash
Standard algorithm
The conventional two phase flash calculation starts with a stability test for the feed composition,
based on Michelsen’s (1982) stability criterion, . If unstable, a flash phase split calculation is
performed to determine the two phase compositions. The multiphase flash continues in a step-
wise fashion: if any one of the current N phases exhibits instability, an N+1 phase flash split
calculation is performed, the resulting tested for stability. The maximum allowed number of
phases is five, after that is reach no more phases are split off.
Stability analysis
The stability of each phase present in a multiphase regime can be tested with the Michelsen
analysis, [Ref. 16], precisely as for the two phase case, using the Gibbs energy tangent plane
distance criterion. However, the initial estimates for phases that will be split off are not limited
to the vapor/liquid equilibrium phases characterised by the Wilson K-values: additional
estimates are chosen from the following phases:
• Each inorganic component
• Lightest hydrocarbon
• Heaviest hydrocarbon
Systems of Equations
∑ βm = 1 [EQ 8.116]
m
These consistency rules can be checked at the end of the multiphase flash calculation.
Defining the multiphase equilibrium ‘K’-values as:
z ij
K ij = ----------
-
z ir ( i )
Where r ( i ) is the reference phase index chosen for component i, usually so that the largest
amount of component i is in the reference phase.
G- =
g° = ------
RT ∑ ∑ βmzim ln ( fim ) [EQ 8.119]
i m
= z j ( ln z jk – ln z jr ( j ) + ln Φ jk – ln Φ jr ( j ) ) [EQ 8.121]
= z j ( ln K jk + ln Φ jk – ln Φ jr ( j ) ) [EQ 8.122]
where Φ jk is the fugacity coefficient of component j in the k’th phase. In equilibrium, the Gibbs
energy is at a minimum, and the first order stationary condition thus implies fugacity balance
between like components in all phases.
The second derivative of the Gibbs free energy is:
δ mk – δ mK⎞ ⎛ δ ij ∂ ln Φ im ∂ ln Φ im⎞
g·· ° = z i z j ⎛⎝ -----------------------
- ⎜ ------- – 1 + -----------------
-– ∑ zpm -----------------
-⎟
βm ⎠ ⎝ z im ∂z im ∂z pm ⎠
p
δ Mk – δ MK⎞ ⎛ δ ij ∂ ln Φ iM ∂ ln Φ iM⎞
– ⎛⎝ ------------------------
- ⎜ -------- – 1 + ------------------
⎠ ⎝z -– ∑ zpM ------------------
-⎟
βM iM ∂z jM ∂z pM ⎠
p
Trace Elements
If a phase possesses a near-zero component composition, then the variable for that component
and phase can be dropped from the variable list (that is not changed in the SS or NR) and its
value calculated after each iteration by invoking fugacity balance between it and the fugacity of
a variable that was update in the usual way. In doing so, the trace element fugacity coefficient
is calculate by setting its composition to zero. If the component transcends the numerical zero
boundary, its variable is reintroduced into the variable set ready for the next iteration.
o
ln x 2i = ln x i + ln Φ i, L ( x i ) – ln Φ i, L [EQ 8.124]
o
ln x 3i = ln x i + ln Φ i, W ( x i ) – ln Φ i, L [EQ 8.125]
o
where Φ denotes pure fugacity, L,V and W denote liquid, vapor and wax. y i component in the
vapor phase. x i , x 2i and x 3i are the components in the lightest liquid, the heavier liquid and
the wax phase.
where the w denotes the component index of the water. The difference between the options is
then implemented in the ‘bring-back’ treatment for trace components used in the phase split
calculations, either none, CO 2 only, or all.
Water- CO 2 bic-fitting
A particular modification to the use of PWR Equation of State is made when only CO2 is
dissolved in water, in allowing for salinity effect CO 2 solubility decrease with salinity, but the
PRW Equation of State makes no allowance for this. It has been possible, however, to establish
a polynomial form for the H 2 O – CO 2 τ ij by iterating adjusting (with a NR scheme) the
H 2 O – CO 2 τ ij until the solubility predicted by the flash calculation matches experimentally
determined solubilities. With a set of such τ ij s one can fit a 3 degree surface polynomial in p,
T, S (salinity):
τ H2 O – CO2 = f1 P, T ) + Sf 2 P, T )
2
fi = a 1i + a 2i T + a 3i P + a4i TP + a 5i T
2 2 3 3
+a 6i P + a 7i T P + a 9i T + a 10i P
1⁄2
a = ∑ ∑ xi xj ( 1 – δij )( ai aj )
i j
b = ∑ xi bi
i
a i = a ci α i
2 2
R T ci
a ci = 0.45724 ---------------
-
P ci
2
κ = 0.37464 + 1.54226ω i – 0.26992ω i
RTci
b i = 0.07780 ----------
-
P ci
Two modification are made for water systems. Firstly, for the water component in any phase, a
correction is made to improve predicted water vapor pressures:
1⁄2 1⁄2
α = 1.0085677 + 0.82514 ( 1 – Tr )
1⁄2
when T r < 0.85
1⁄2
a= ∑ ∑ xi xj ( 1 – τij ( T ) )( ai aj )
i j
when Tr ≥ 0.85 , or the phase is not the aqueous liquid phase, the original PR Equation of State
is used.
The temperature τ ij s required may be found from the published curves given by Peng &
Robinson. The following forms have been derived:
τ C1 ↔ H2 O = 1.659Tr P rw – 0.761
τ C2 ↔ H2 O = 2.109Tr P rw – 0.607
2 2
τ C3 ↔ H2 O = – 18.032T r P rw + 0.9441Tr P rw – 1.208
τ C4 ↔ H2 O = 2.800Tr P rw – 0.488
T
τ Cn ↔ H2 O = ⎛⎝ 0.4 -----c⎞⎠ T r P rw – 0.8 , n > 4
Pc
τ N2 ↔ H2 O = 0.402Tr – 1.586
τ H2 S ↔ H2 O = 0.22T r – 0.19
2
τ CO2 ↔ H2 O = – 0.074T r + 0.478T r – 0.503
Viscosity correlations
Three correlations are available in PVTi for the estimation of viscosities, namely that due to
"Lohrenz, Bray and Clark" on page 330, [Ref. 7], the method due to "Pedersen et al." on
page 331 (PED), [Ref. 5] and [Ref. 36], and also that of "Aasberg-Petersen et al" on page 333,
[Ref. 63].
a 2 = 0.0233640
a 3 = 0.0585330
a 4 = – 0.0407580
and η∗ is the low-pressure gas mixture viscosity. ξ is the viscosity-reducing parameter, which
for a fluid mixture is given by:
N 1⁄6 N –1 ⁄ 2 N –2 ⁄ 3
ξ = ∑ z i T ci ∑ z i M wi ∑ z i p ci [EQ 8.128]
i=1 i=1 i=1
The dilute gas mixture viscosity is as given by Herning and Zippener, [Ref. 28]:
N N –1
1⁄2 1⁄2
η∗ = ∑ z i η i∗ M wi ∑ z i M wi [EQ 8.131]
i=1 i=1
where the dilute gas viscosities of the individual components, ηi∗ are derived from expressions
due to Stiel and Thodos, [Ref. 29]:
– 5 1 0.94
η i∗ = 34 ×10 ---- T ri T ri < 1.5 [EQ 8.132]
ξi
–5 1 0.625
η i∗ = 17.78 ×10 ---- ( 4.58T ri – 1.67 ) T ri > 1.5 [EQ 8.133]
ξi
where
1⁄6 –1 ⁄ 2 – 2 ⁄ 3
ξ i = T ci M wi p ci [EQ 8.134]
Pedersen et al.
Viscosities can be calculated from a modified form of the corresponding states method. A group
of substances obey the corresponding states principle if the functional dependence of the
reduced viscosity, η r , say, on reduced density and temperature, ρ r and T r , say, is the same for
all components within the group, namely:
η r ( ρ, T ) = f ( η r ,T r ) [EQ 8.135]
in which case comprehensive viscosity data is only needed for one component of the group,
which is denoted as the reference substance, to be given the subscript ( o ) all other components
are identified with the subscript ( x ).
Generally, the critical viscosity, ηc , is not known but it can be estimated from the inverse of [EQ
8.135]),
–1 ⁄ 2 2 ⁄ 3 –1 ⁄ 2
η r ( ρ, T ) = η ( ρ, T )T c Vc Mw [EQ 8.136]
Thus, the viscosity of component x at temperature T and a pressure with density ρ , is given by:
1⁄2 –2 ⁄ 3 1⁄2
T cx V cx M wx
η r ( ρ ,T ) = ----------
- ------------- ------------ η ( ρ ,T )
T co V co M mo o o o
[EQ 8.137]
Oil mixtures contain a range of components with molecular weights ranging from 16 (Methane)
to about 1100 ( C 80 ). It cannot be expected for C 1 and C80 to belong to a group where a simple
corresponding states principle applies. Generally, a third parameter must be added to account
for the shape of the molecules, such as the acentric factor.
Pedersen et al., [Ref. 5], presented a corresponding states principle depending on reduced
temperature and reduced pressure, η r = f ( p r, T r ) where:
η
η r = ηξ = ------------------------------------- [EQ 8.138]
–1 ⁄ 6 2 ⁄ 3 1 ⁄ 2
Tc pc Mw
The deviation from the simple corresponding states principle is expressed in terms of a
rotational coupling coefficient, denoted α , to give:
–1 ⁄ 6 2⁄3 1⁄2
T c, mix p c, mix M w, mix α mix
η mix ( p, T ) = ------------------------
T co
---------------------- ------------------------- ----------- η o ( p o, T o )
p co M wo αo
[EQ 8.139]
where
p co α o
p o = p -------------------------
p c, mix α mix
[EQ 8.140]
T co α o
T o = T -------------------------
T c, mix α mix
- [EQ 8.141]
The critical temperature and volume for unlike pairs of molecules are given by:
where the critical volume of a component can be expressed in terms of the critical temperature
and pressure using the real gas law and the critical Z -factor.
Assuming a constant Z c for all components,[EQ 8.143] becomes:
1 T c, i⎞ 1 ⁄ 3 ⎛ T c, j⎞ 1 ⁄ 3⎞ 3
V c, ij = --- constant ⎛⎝ ⎛⎝ --------
- + ⎝ ---------⎠ ⎠ [EQ 8.144]
8 p c, i⎠ p c, j
N N
∑ ∑ zi zj Tc, ij Vc, ij
i = 1j = 1
T c, mix = -------------------------------------------------
N N
[EQ 8.145]
∑ ∑ zi zj Vc, ij
i = 1j = 1
N N
1⁄3 1⁄3 3
⎛T c, i⎞ ⎛T c, j⎞ 1⁄2
∑ ∑ zi zj --------
⎝ p c, i⎠
- --------
⎝ p c, j⎠
- [ T c, i T c, j ]
T c, mix = i----------------------------------------------------------------------------------------------------------------
= 1j = 1
N N
[EQ 8.146]
T c, i⎞ 1 ⁄ 3 ⎛ T c, j⎞ 1 ⁄ 3 3
∑∑ z i z j ⎛⎝ -------- - ---------
p c, i⎠ ⎝ p c, j⎠
i = 1j = 1
N N
1⁄3 1⁄3 3
⎛T c, i⎞ ⎛T c, j⎞ 1⁄2
8 ∑ ∑ zi zj --------
-
⎝ p c, i ⎠
--------
-
⎝ p c, j ⎠ [ T c, i T c, j ]
i = 1j = 1
p c, mix = -------------------------------------------------------------------------------------------------------------------
2
- [EQ 8.147]
⎛ N N ⎞
⎜ ⎛ T c, i⎞ 1 ⁄ 3 ⎛ T c, j⎞ 1 ⁄ 3 3⎟
⎜
⎝i = 1 j = 1
∑∑ z i z j ⎝ ---------⎠ ⎝ ---------⎠
p c, i p c, j ⎟
⎠
where M w and M n are the weight average and number average mole weights, respectively.
This mixing rule is derived empirically on the basis of available viscosity data and assigns a
larger influence to heavier components.
The α -parameter for the mixture is found from:
– 3 1.847 0.5173
α mix = 1.0 + 7.378 ×10 ρ r M w, mix [EQ 8.149]
where the α of Methane, the reference substance, is given by:
1.847
α 0 = 1.0 + 0.031ρ r [EQ 8.150]
It has been our experience that the method of Pedersen gives much better prediction of
viscosities than does the method of Lohrenz, Bray and Clark. Neither method is particularly
good but typical errors based on un-regressed Equation of State data are, LBC ∼ 50 % and
PED ∼ 90 % of the reported values. However, you should treat experimentally measured
viscosities with some caution, as viscosities are often calculated.
Aasberg-Petersen et al
The Aasberg-Petersen model is also based on the principle of corresponding states, but uses two
references fluids, methane and decane, instead of just the one (methane) in the case of the
Pedersen et al. model. The usual problem with the Pedersen model is that inaccurate predictions
occur for fluids with components that are significantly different in molecular weight to the
reference component (methane).
The idea of this model is to use a heavier second reference component and effectively create an
optimum reference component by using the molecular weight of the fluid as an interpolation
parameter between the two reference components. Decane is chosen as the second reference
component because it is the heaviest alkane for which a significant amount of experimental
viscosity data is known. Methane is a natural choice as the first reference component because
of its presence in large mole fractions in most reservoir fluid mixtures
The interpolation law is used to calculate the reduced viscosity of the optimum reference
component (denoted with the subscript x) using the reduced viscosities of methane and decane
and is obtained using the following expression:
MW x – MW 1 ⎛ η r2⎞
ln η rx = ln η r1 + ------------------------------
- ln -------- [EQ 8.151]
MW2 – MW1 ⎝ η r1⎠
MW is the molecular weight and subscripts 1 and 2 refer to the reference components. The
functional form of equation [EQ 8.151] was originally suggested by Teja and Rice (1981), [Ref.
64], using the acentric factor instead of MW. This is not possible in the present work, since the
acentric factor decreases with increasing molecular weight for heavy oil fractions.
The reduced properties are determined from:
E
Er = -----, E = T, P, η [EQ 8.152]
Ec
Subscripts r and c indicated reduced and critical properties respectively. The following
expression is used to evaluate the critical viscosity (Pedersen et al. 1989):
1
--- 2
--- – 1
---
2 3 6
ηc = C ⋅ MW P c Tc [EQ 8.153]
C is a constant. From equations [EQ 8.151] and [EQ 8.142][EQ 8.142] the following equations
may be derived for determination of the viscosity:
η cx η 1 ( T 1, P 1 ) η 2 ( T 2, P2 )ηc1 K
ηx = ---------------------------------
- ---------------------------------- [EQ 8.154]
ηc1 η 1 ( T 1, P1 )ηc2
MW x – MW 1
K = ------------------------------
- [EQ 8.155]
MW2 – MW1
P ci P
- , i = 1, 2
Pi = ---------- [EQ 8.157]
P cx
The model outlined in the above is extended to mixtures using the same mixing rule as the
Pedersen et al. model, that is equations [EQ 8.146] and [EQ 8.147]. The mixture molecular
weight is calculated using the formula:
1.56079 1.56079
MWmix = MW n + 0.00867358 ( MW w – MW n ) [EQ 8.158]
N
2
∑ zi MWi
MWw = i-------------------------
=1 -
N
∑ zi MWi
i=1
and
The constants in equation [EQ 8.142] are determined by regression using experimental viscosity
data for binary mixtures and oils.
The following equations are used for determination of the viscosity of the reference
components:
η0 = ηk ( T ) + ρη 1 ( T ) + η2 ( ρ, T ) [EQ 8.161]
9 (-------------
i – 4 -)
3
ηk ( T ) = ∑ GVi T [EQ 8.162]
i=1
T 2
η1 ( T ) = A + B C – ln ⎛⎝ ---⎞⎠ [EQ 8.163]
F
j
η2 ( T, ρ ) = H 2 exp ⎛⎝ j1 + ---4-⎞⎠ [EQ 8.164]
T
0.1 ⎛ j3 ⎞ j j7 ⎞
- + θρ ⎛⎝ j5 + ---6- + -----
0.5
H 2 = – 1 + exp ρ ⎝ j2 + -------- [EQ 8.165]
1.5⎠ T T 2⎠
T
These equations were originally developed by Hanley et al. (1975), [Ref. 65], to correlate the
methane viscosity. For methane the values of the GVi parameters in equation [EQ 8.152] given
by Hanley et al. were maintained. For decane GV4-GV9 are equal to zero and GV1-GV3 were
taken from the DIPPR tables (1985). All the parameters in equations [EQ 8.59]-[EQ 8.165] were
estimated. In Table 8.9, data of the parameter estimation are given and Table 8.10 lists all
parameters for the pure component viscosity correlation. Of the parameters in Table 8.10 are
used the viscosity is obtained in μP if the density is given in g/cm3.
Methane Decane
N 881 252
T-range (K) 91-523 244-477
P-range (atm.) 0-680 0-1000
Dev. (%) 3.1 3.8
Methane Decane
GV1 -209097 0.2640
j2 30.574 16.092
j3 -14929 -18464
Methane Decane
j4 1061.5 -811.3
j5 -1.4748 1.9745
j6 290.62 898.45
j7 30396 119620
The density of the reference components as a function of the pressure and temperature is also
required in equations [EQ 8.48] and [EQ 8.165]. For methane, as with the Pedersen model, the
33-parameter MBWR-equation given by McCarthy (1974), [Ref. 66], was used to calculate the
density. For decane we decided to used the flash available within PVTi to obtain the density at
a given pressure and temperature.
Finally, the following physical properties for methane and decane were used in Table 8.11 below
when required.
Methane Decane
Tc 190.55 617.40
Pc 45.39 20.18
MW 16.043 142.284
ω 0.008 0.484
ρc 0.1649 0.2269
The predictions of the Aasberg-Petersen model have been shown to agree well with
experimental data over large pressure and temperature ranges. In particular this model is able to
handle mixtures with CO2, paraffinic and aromatic components with better accuracy than the
Pedersen model. A better match is also obtained than the Pedersen model for heavier oils, where
the size and shape of the molecules differ substantially from the single reference component
(methane).
Note It should be noted that the Aasberg-Petersen model is not recommended for fluids with
large concentrations of napthalenic components.
p = RT a- + …
------- + ----- [EQ 8.167]
V V2
Caution Many authors have suggested that binaries are the obvious Equation of State
parameter to adjust to match Equation of State to laboratory results, especially the
Methane to plus-fraction binary. However, Pedersen et al., [Ref. 30], have shown
that this is problematic.
Given that our preferred Equation of State is the Peng-Robinson, we derive binaries from one
of two sources.
The current default is to assume they are zero for all hydrocarbon-hydrocarbon interactions
except for Methane to heavier hydrocarbons that vary like:
KC
1, j = 0.14γ j – 0.06 [EQ 8.168]
where γ j is the (liquid) specific gravity of the heavier component. In addition, the hydrocarbon
to non-hydrocarbon interactions take certain fixed values, see [Ref. 4].
The alternative set, activated by the appropriate option switch is a set due to Cheuh and
Prausnitz, see [Ref. 32].
The hydrocarbon to non-hydrocarbon interactions are as above but the hydrocarbon to
hydrocarbon binaries are given by:
⎛ 2 ( V c, i V c, j ) 1 ⁄ 6⎞
6
K i, j = A 1 – ⎜ ----------------------------------
-⎟ [EQ 8.169]
⎝ V 1c, ⁄i3 + V 1c, ⁄ j3 ⎠
where V c, i is the molar volume of the i th component and A is a special regression variable which
is generally in the range 0.15 ≤ A ≤ 0.25 .
There is some appeal in using the pre-multiplying coefficient to regress all binaries together.
This preserves symmetry, which might be lost using the Methane to plus fraction binary
approach of [EQ 8.169], which (see [EQ 8.169] and [Ref. 28]) can lead to strange results.
These particular binaries are also given a temperature dependence by the use of an additional
multiplying factor that has the form:
1.0 + TC ( T – T std ) [EQ 8.170]
where TC has the fixed value of 0.0025.
Introduction
This section describes the basic laboratory experiments and how PVTi attempts to simulate
them. More detail on experimental methods can be found in the excellent text by Pedersen et al.
[Ref. 30].
Schematic diagrams of the apparatus used for the standard laboratory techniques to be described
here can be found in the above text.
See "The fluid model" on page 98 and "Batch system and keywords" on page 152 for further
information on defining experiments in PVTi.
This section contains information on the following:
• "Blackoil systems" on page 338.
• "Gas condensate systems" on page 340.
• "Gas injection processes" on page 343.
• "Process simulation" on page 342.
Blackoil systems
Essentially two experiments are performed on black or crude oil systems:
• "Bubble point evaluation" on page 338
• "Differential liberation" on page 339 (sometimes referred to differential expansion).
PVTi
1 PVTi finds the bubble point by taking one mole of the reservoir fluid and using estimates
for the K -values, constructs vapor compositions at a set of pressures, testing for the stability
of the second (vapor) phase using the method of Michelsen [Ref. 16].
2 Once the highest-pressure, two-phase state has been found, the N equal fugacity conditions
and the one mole fraction constraint (see "Material balance" on page 316 and "Flash
calculations" on page 316) are used to solve for the N + 1 unknowns of N K -values and the
bubble point pressure by iteration.
Laboratory experiment
1 Starting from the bubble point pressure, the pressure is dropped by several hundred psia,
whereupon a volume of gas is evolved that is removed from the cell at the (new) constant
pressure. This process is repeated several times noting the volume of gas evolved (at
reservoir and surface conditions) and the volume of liquid remaining.
2 The volume of liquid remaining at the last stage, which should be at standard or
atmospheric pressure, is then reduced to standard or atmospheric temperature and re-
measured. This liquid is discharged and its density is obtained. The density of the liquid at
the other pressure (and temperature) stages can be calculated from mass balance knowing
the volumes and molecular weight of the removed gas streams.
The data usually quoted is:
• Bo Oil formation volume factor
• Bg Gas formation volume factor
• Rs Gas-oil ratio
• ρo Oil density
• γg Gas gravity (density)
• Zg Gas deviation factor
• μo Oil viscosity
• μg Gas viscosity
where:
Vo ( p )
B o ( p ) = -------------
std
- [EQ 8.171]
Vo
Vg ( p )
B g ( p ) = ----------------
std
- [EQ 8.172]
Vg ( p )
i = std
1-
∑ Vg
std
R s ( p ) = --------
std
(i) [EQ 8.173]
Vo
i=p
and
Mg ( p )
γ g = --------------- [EQ 8.174]
M air
Here, Vo , Vg are volumes of oil and gas, M g , M air are mole weights of hydrocarbon gas and
air, the superscript std implies standard conditions, and the summation for the gas-oil ratio
is taken over the volumes of evolved gas from the current pressure, p to the final pressure
at std .
Note The "Definition of GOR in Diff. Lib." on page 148 program option allows three further
definitions of the GOR.
The first removes the stage to standard conditions and normalizes gas volumes to the
volume of oil at reservoir conditions, that is:
i = plast
1
∑
std
R s ( p ) = ---------------------------
V o ( T, p bub )
Vg ( i ) [EQ 8.175]
i=p
where T is the temperature of the depletion experiment, and p last is the last pressure stage
specified.
2 The second definition specifies the GOR as an incremental one, that is:
std
Vg ( p )
R s ( p ) = ----------------
std
- [EQ 8.176]
Vo
3 The third definition is the same as the default but the volume of oil is at its bubble point
pressure rather than at stock tank conditions, that is:
i = std
1
∑ Vg
std
R s ( p ) = ---------------------------
V o ( T, p bub )
(i) [EQ 8.177]
i=p
The program option "Definition of Oil relative volume in Diff. Lib." on page 149 allows an
alternative definition of the oil relative volume where the volume of oil is normalized to the
initial volume of oil at reservoir rather than standard conditions:
Vo ( p )
B o ( p ) = --------------------------- [EQ 8.178]
V o ( T, p bub )
Hint It is not uncommon for CCEs to be performed on all fluids, and it is recommended to
perform CVDs on more volatile oils.
Note This can be affected by imperfections or grease in the cell, poor experimental
procedures, etc., and it is not uncommon for errors of 100 psia to be associated with
this measurement.
PVTi
In PVTi, this experiment is simulated in much the same way as the bubble point evaluation
except that the trial second (liquid) phase is heavier than the original fluid.
1 The same stability test is done to find the highest two phase state and the same set of
equations and variables iterated to solution.
2 As an alternative to this high-pressure dew point, usually referred to as the retrograde dew
point, the low-pressure or normal dew point can be found by starting low in pressure and
iterating higher.
PVTi
In PVTi, the CCE is simulated by:
1 Finding the saturation pressure which for one mole of feedstream defines the cell or control
volume.
2 Then at pressures above the saturation pressure it is sufficient to find the single phase liquid
or vapor Z -factor to calculate the volumetric behavior.
3 Below the saturation pressure, the feedstream is flashed at each required pressure stage to
determine the phase split and other volumetric properties.
and occasionally:
x i, N : Composition of liquid stream left in cell at last pressure stage;
vap
M N+ : Mole weight of removed vapor plus fraction;
vap
γ N+ : Specific gravity of removed vapor plus fraction.
PVTi
In PVTi this experiment is simulated by:
1 First finding the saturation pressure. The volume occupied by one mole of fluid at its
saturation pressure then defines the cell or control volume.
2 At some pressure less than the saturation pressure, the fluid is flashed into two phases and
any excess gas is removed to return the volume to the control volume.
3 The number of moles of vapor and its properties are noted.
4 The composition of the remaining fluid is calculated by volumetric balance and this forms
the feedstream for the next pressure stage, and so on.
Process simulation
For the purposes of process-type simulation and/or for the definition of blackoil tables for
reservoir modeling, separator tests are performed to see what phase splits are achieved when a
fluid is flashed at a series of pressures and temperatures in some prescribed sequence.
PVTi
1 In PVTi this process is simulated by a set of flash calculations at the required set of
pressures and temperatures, taking the reservoir fluid as the feedstream and routing the
liquid and vapor outputs to other stages.
2 The default output streams are liquid to the next stage and vapor to the stock-tank
accumulation, but PVTi is capable of solving feedback loops, where, for example, vapor
output is routed back to an earlier stage.
3 In addition, a user option allows the output streams to be split and directed to more than one
other stage.
Optimized separators
In the Fluid Properties (FPE) workflow, PVTi offers the automatic generation of the optimum
separator configuration for a fluid. The optimum separator configuration is defined as the two-
stage separator for which the stock-tank formation volume factor (Bo) is minimized.
Given the reservoir temperature and the maximum separator pressure, the following separator
chain is constructed:
1 A separator stage at some temperature below the reservoir temperature and some pressure
below the highest separator pressure. This stage is optimized.
2 A second stage at standard conditions, the liquid output from Stage 1 passes into this second
stage.
3 The stock tank. The vapor output from Stages 1 and 2 both pass into the vapor stock tank
and the liquid output from Stage 2 passes into the liquid stock tank.
A special regression procedure is used to determine the optimum temperature and pressure for
Stage 1, such that Bo is minimized in the stock tank. This point coincides with minimum Total
Gas-Oil Ratio and minimum stock tank vapor Gravity.
For more information on fluid properties estimation see "Fluid Properties Estimation" on
page 384.
Swelling test
Laboratory experiment
1 The swelling test consists of finding the saturation pressure and hence volume of a reservoir
fluid.
2 Followed by adding, in a series of steps, prescribed volumes (or moles) of lean injection
gas, re-pressuring the resulting mixture to return to a single phase system and measuring
the new saturation pressure and volume.
The data quoted is the set of saturation pressures for the original fluid and the mixtures and the
ratio of the saturation volume of the mixtures to the saturation volume of the original fluid,
usually referred to as the swelling factor.
PVTi
In PVTi this is simulated by:
1 Firstly finding the saturation pressure of one mole of the required reservoir fluid and hence
the saturation volume.
2 Then prescribed volumes of a lean injection gas are added as a GOR (volume of injection
gas at standard conditions per volume of original reservoir fluid at its Psat or other
prescribed pressure) or a mole% (moles of lean gas per moles of mixture) to give a new
fluid composition.
3 The saturation pressure and volume of the new mixture are found, and hence the swelling
factor.
Vaporization test
This is similar to the swelling test, except that it is performed at constant volume and pressure.
Laboratory experiment
1 A volume of reservoir fluid at some pressure and temperature, usually below the saturation
pressure and hence two phase, is contacted by a series of lean gas injections.
2 After each contact, a volume of (enriched) gas and/or liquid is removed to return the system
to the original volume. The composition of the removed gas is measured.
The experiment seeks to measure the extent of vaporization of intermediate and heavy
components from the reservoir liquid phase by stripping into the injected gas stream.
In PVTi, the procedure is as follows.
1 One mole of reservoir fluid at a prescribed pressure and temperature is flashed to calculate
the phase split and volumes of liquid and vapor.
Note However, the compositional gradient experiment assumes many conditions: namely
thermal, gravitational and diffusive equilibrium. Any or all of these conditions can be
violated in a given reservoir, thus invalidating the use of [EQ 8.179]).
With the above caveats, the experiment may still be of use in the prediction of the existence and
possible location of either gas-oil contacts or so-called critical transitions. This latter transition
consists of the grading of a fluid from an oil to a gas, or vice-versa, without passing through a
gas-oil contact. This occurs when an over-pressured reservoir has a fluid composition at some
depth whose critical temperature is equal to the reservoir temperature at that depth. There is
evidence that such behavior exists in some North Sea fields.
FCMP experiment
Refer to Jensen and Michelsen [Ref. 38] and Pedersen et al. [Ref. 30] for detail.
MCMP experiment
Refer to Jensen and Michelsen [Ref. 38] and Pedersen et al. [Ref. 30] for detail.
Tsat experiment
Refer to Michelsen [Ref. 15] for further details.
Introduction
This section of the manual contains information on:
• "Practical considerations" on page 347.
• "Theoretical model" on page 350.
Practical considerations
There are no set rules for how to do regression of an equation of state model to match to
laboratory measurements.
The paper by Coats and Smart, [Ref. 27], contains an appendix on the choice, selection and
range limits of regression variables.
However, the Coats and Smart model is limited in its choice of regression variables to the Ωa ’s,
Ω b ’s and the binary interaction coefficients.
Another limitation of their model is the use of the standard two-parameter equation of state. It
is well known that the two-parameter Equation of State is in error in prediction of liquid
properties by as much as 10%. Therefore, to get matches to saturation pressures and densities it
becomes necessary to change the properties of “well-defined” components, say Ωa ( C 1 ) and
Ωb ( C1 ) .
Such a problem is avoided in PVTi by use of the PR3 or SRK3 Equation of State, which allows
the volume shift parameter to be a possible regression variable.
Additionally, rather than varying just the Ωa ’s and Ωb ’s, PVTi allows you to change the critical
pressures and temperatures. This has the advantage in that monotonicity tests can readily be
applied to the set of T c ’s, p c ’s,..., etc., for the hydrocarbon components to ensure that critical
temperatures increase with increasing mole weight, critical pressures decrease (except C 1 , C 2 ,
C 3 ) etc. For further information on performing a regression in PVTi see "Regression in PVTi"
on page 126 and "REGRESS section keywords" on page 160.
The present section also contains information on:
• "Consistency and quality of measured data" on page 347
• "Plus fraction" on page 348
• "Performing regressions in PVTi" on page 349.
Plus fraction
Having determined that the data is reliable, or having rejected either poor data or adjusted it to
be consistent, one will usually find that the Equation of State predictions differ from the
measured data.
Most of this error can probably be associated with the incomplete fluid description, namely the
failure to fully characterize the plus fraction, for example C 7+ , although some error has to be
associated with the inadequacies of a cubic Equation of State.
Clearly, the plus fraction, consisting as it does of many hundreds or possibly thousands of
components, cannot be represented by just one component without some modification. Even
then, just to represent a C 7+ of a condensate or volatile oil, which may be 5-20 mole percent of
the fluid, by a single component with regressed properties may be insufficient for describing a
multi-pressure process.
On some CVD analyses the mole weight and specific gravity of the produced gas plus fraction
are measured and reported.
In such a case you may note that the plus fraction gets lighter as one would expect. Then it is
clearly not adequate to represent the plus fraction by a single component and one should
consider splitting it into two or three pseudo-components, say using the techniques available in
PVTi.
The properties of the plus fraction, or its pseudo-components if splitting has been performed are
the obvious candidates for regression. However, one should avoid regressing the property of a
component with a small mole fraction, say 1%.
Note The LBC correlation can be as much as 150% in error and is regularly 50% in error.
Theoretical model
Generally the results of an equation of state model must be tuned by regression of one or more
variables, x = ( x1, …, x N )T to a set of laboratory experiments, y = ( y1, …, y M )T [EQ 8.180]
˜ ˜
where M ≥ N and:
yi = yi ( x )
˜
or in residual form:
ri ( x ) = wi [ ( y i – yi ( x ) ) ⁄ yi ] [EQ 8.181]
˜ ˜
Regression algorithm
The regression algorithm in PVTi seeks to minimize the least squares residual ( l 2 norm) given
by the objective function:
M
1 1
∑ ri ( x˜ )
2 T
f ( x ) = ---
2
= --- R ( x ) R ( x )
2 ˜
[EQ 8.182]
˜ ˜
i=1
Termination conditions
A variety of termination conditions from the regression algorithm are possible.
Note The required condition is that the residual goes to some small value, such as 10 –6 , but
this is rarely achieved for larger problems.
An alternative is that the gradient of the objective function goes to zero, indicating that the
regression is approaching a minimum and no further improvement can be expected.
As a result of one or more steps made by the algorithm, one or more of the selected regression
parameters can be pushed out of physical range.
What constitutes a physical range is debatable, but it is suggested that an increase or decrease
in excess of 50% of the original value causes termination. You can change these limits but care
should be taken.
You can request a premature termination by allowing only a small number of iterations, say five
or ten.
Hint This is a good practice to follow for any regressions which involve several variables,
experiments and observations.
Finally, if any errors are detected in any of the experiments during the regression, then
depending on their nature, the program may terminate the regression.
General information
Our suite of reservoir simulation software supports five different models for fluid behavior.
There is the ECLIPSE extended blackoil model, the pseudo-compositional model (GI option),
the ECLIPSE compositional model., the ECLIPSE Thermal model and the API Tracking option
used in ECLIPSE BlackOil.
p -V -Tdata suitable for use in any of these models can be generated from the PVTi package.
Each of the models and typical output is discussed in the appropriate section.
This section of the manual contains information on:
• "Blackoil model" on page 353.
• "Differential and composite from differential tables" on page 360.
• "Pseudo-compositional tables for ECLIPSE GI option" on page 361.
• "Compositional data for ECLIPSE Compositional" on page 362.
• "Water properties" on page 362.
• "Model for API Tracking option in ECLIPSE BlackOil" on page 364.
• "Compositional Data for ECLIPSE Thermal" on page 367.
Blackoil model
ECLIPSE has a so-called extended blackoil model. That is, in addition to the standard blackoil
parameters, Rs , B o and B g , it contains the parameter R v for modeling oil vaporization in the gas
stream.
Methods of obtaining black oil tables from compositional data are described by Whitson and
Torp, [Ref. 6] and Coats, [Ref. 3].
In both cases the basis of the method is a constant volume depletion experiment, used to supply
reservoir liquid and vapor compositions at a series of pressures. (In the case of crude oil
samples, a differential liberation process is used in the same way.)
The blackoil model, used by ECLIPSE, can be viewed as a two-component compositional
model. The “components” are stock tank oil and stock tank gas that are assumed to be invariant
and are assigned constant densities.
PVT model
The actual form of the PVT model for the hydrocarbons depends on whether:
• There is or could be liquid and/or vapor in the reservoir during its production.
• The reservoir phases produce stock tank oil and/or gas on flashing to surface conditions.
M gr = M og + M gg [EQ 8.185]
M or = M go + M oo [EQ 8.186]
where ( M gr, M or ) are the reservoir molar masses and ( M og, M gg ) , ( M go, M oo ) are the corresponding
actual surface masses after flashing.
Since ρ = M ⁄ V , we can write equations [EQ 8.185] and [EQ 8.186] as:
ρ gr V gr = ρ og V og + ρ gg V gg [EQ 8.187]
ρ or V or = ρ go V go + ρ oo V oo [EQ 8.188]
or:
1
ρ gr = ------ ( ρ gg + R v ρ og )
Bg
[EQ 8.189]
1
ρ or = ------ ( ρ oo + R s ρ go )
Bo
[EQ 8.190]
where:
V go V og
R s = -------
- R v = -------
- [EQ 8.191]
V oo V gg
Note Strictly, the blackoil model requires that the stock tank “components” are constant and
invariant with time.
CVD process
In the Coats method, only the reservoir vapor is taken through the separator, the oil B o and R s
values being obtained by solving the mass conservation equations:
ρ sto V 2 ( b g S g + b o R s S o ) 2 = ρ sto V 1 ( b g S g + b o R s S o ) 1 [EQ 8.193]
for each expansion step of the CVD process from V 1 to V2 , where ρ sto and ρ stg are fixed surface
densities, and b o = 1 ⁄ Bo , b g = 1 ⁄ B g .
The stock tank densities are obtained from the output of the separators at the saturation pressure.
Solving these equations yields values at all pressures except the saturation pressure, at which
the Whitson method can be used.
Alternatively, a small initial pressure step from the saturation pressure may be specified.
Note The constant volume expansion method is usually applied to condensates, but can also
be used for volatile oils.
Differential liberation
The differential liberation approach is only suitable for oils.
In this case, the analogue of the Whitson method runs both reservoir oil and vapor through the
separators as before.
The analogue of the Coats method uses the conservation equations, [EQ 8.193] & [EQ 8.194]
(although in this case S g1 is zero for each step), ρ sto and ρ stg are obtained from the reservoir oil
at the saturation pressure.
Hint The "Choosing the unit type for PVTi" on page 144 program option gives you a choice
of outputting and plotting all the data (saturated and undersaturated) or just the
saturated data.
Both the PVTO and PVTG tables are extended to include properties of the undersaturated oil and
gas, respectively.
In the case of the PVTG table the model computes, at all values of Rv between the saturated
value and the dry gas value ( Rv = 0.0 ), a value for Bg by adding sufficient stock tank oil to the
stock tank gas (both from the flash of the saturated gas to stock tank) to give the required R v and
then taking the ratio of this fluid's volume at reservoir and stock tank volumes.
In the case of the PVTO table, the model compresses the saturated fluid from a given Rs node at
all pressures from the saturated pressure up to the highest pressure. The ratio of these volumes
to the resulting stock tank oil volume (from flashing the saturated oil) gives the undersaturated
B o values.
1 – dB dR ( B o – R s B g )
C g = ------ -----------g- + --------v ----------------------------
- [EQ 8.198]
B g dp dp ( 1 – R s R v )
Alternatively if you do not want these generated points, PVTi allows you to increase the
saturated pressure of the sample to above the maximum pressure used in the depletion
experiment by mixing the sample with a lighter gas sample.
The gas sample used is the vapor split-off obtained from performing a p sat calculation on the
original sample.
The model mixes this vapor with the sample and performs a new p sat calculation on the new
sample.
It then mixes the resulting vapor split off with the new sample and repeats the process until a
p sat forms that is higher than the maximum pressure stage. If this process is successful you can
save the final sample and use it in the calculation of the black oil tables.
Note Oil and gas properties below the saturation point generated as described above are also
subjected to the ECLIPSE total compressibility test.
Note If negative compressibilities are detected, PVTi (with confirmation) solves for
intermediate points, which provide the necessary resolution to pass the compressibility
test.
Note These definitions are detailed in the manual (refer to "Basic laboratory experiments"
on page 338 and "Gas condensate systems" on page 340, and to equations [EQ 8.171]
and [EQ 8.173]). The quantities are known here as B Do and RDs .
From these differential quantities, it is possible to estimate composite values of Bo and R s using
a number of formulae.
Two of the more popular methods are due to Moses, , and Fragor, [Ref. 40] (both now termed
CM and CF respectively).
Both methods start by estimating a value for Bo and R s at the saturated pressure ( p sat ) of the
initial oil. This oil is passed through a separator chain and the values of B CF CM
o ( p sat ) and B o ( p sat )
are given by the oil formation volume factor from the separator. Likewise the values of R CF
s ( p sat )
and RCM
s ( p sat ) are given the value of the gas-oil ratio from the separator.
Once these two values are computed, the rest of the values from the lowest pressure up to psat
are found from the following formulae.
Firstly the Fragor definitions:
CF D
CF ( B o ( p sat ) – 1 ) ( B o ( p ) – 1 )
B o ( p ) = 1 + -----------------------------------------------------------------
D
- [EQ 8.199]
( B o ( p sat ) – 1 )
and
CF D
CF R s ( p sat )R s ( p )
R s ( p ) = --------------------------------------
D
[EQ 8.200]
R s ( p sat )
and
CM CM D D ⎛ B CMo ( p sat )⎞
R s ( p ) = R s ( p sat ) – ( R s ( p sat ) – R s ( p ) ) ⎜ -----------------------
-⎟ [EQ 8.202]
⎝ BD o ( p sat ) ⎠
PVTi can provide all four types of black oil output, that is, the usual composite, the differential,
and the two composite from differential tables.
Note The output form for the differential and composite from differential tables is by way
of the PVTO/PVDO and PVDG keywords only.
ECLIPSE GI option
One possible solution to this problem is to extend the three-parameter blackoil model as used in
ECLIPSE BlackOil, that is p , S w and So , by adding a fourth parameter/equation.
In the ECLIPSE pseudo-compositional model, this fourth parameter, denoted GI, is the
cumulative gas volume, which has passed over the volume of reservoir liquid in a gridblock. As
such, GI is equivalent to a gas-oil ratio of volume of injection gas at standard conditions to
volume of oil at reservoir conditions:
gas
V stc
GI = ---------
oil
[EQ 8.203]
V res
This model determines the oil volume at reservoir conditions, V oilres , for the first stage of gas
addition, from the GI=0.0, that is the standard blackoil model. It then calculates the gas volume
at standard conditions, Vgas
stc , from the gas law:
add
gas n ZRT stc
V stc = -------------------------
- [EQ 8.204]
p stc
where n add is the number of moles of gas added. By convention, at standard pressure and
temperature, p stc = 14.7 psia, T stc = 60.0 °F , and the gas-compressibility factor, Z = 1 .
Having defined a Constant Volume Depletion (CVD) experiment and a separator network with
the last stage corresponding to stock tank conditions, the model then performs a calculation of
the conventional extended blackoil tables. This defines the GI=0.0 data.
Compositional mode
The Equation of State mode is essentially the same as that employed by PVTi for its flash,
saturation pressure, etc., calculations.
You must specify the number of components, the Equation of State required (the same set as
available in PVTi), and the critical properties, acentric factors, binary interaction coefficients,
compositions, etc.
Clearly, the more components used, the more accurate the simulation (assuming the component
set has been optimised to measured data); however, such simulations require more CPU time.
Hint If running ECLIPSE Compositional in fully implicit mode, more than six or seven
components may become prohibitive for all but the largest of supercomputers or
workstations, in which case you should consider pseudoisation of this fluid system to
fewer components.
Water properties
ECLIPSE and VFPi consider water to be non-volatile (only exists as liquid water) and
immiscible with the hydrocarbon phases (water cannot dissolve in the hydrocarbons or vice-
versa). This simplifies the PVT treatment of water.
Vw ( p )
B w ( p ) = --------------
- [EQ 8.205]
V stw
which has units of rm3 /sm 3 (metric), rb/stb (field) or rcm3 /scm3 (lab).
Note Note that VFP only supports metric and field units.
Within the models, B w is defined in terms of the Bw at some reference pressure, p ref :
B w ( p ref )
B w ( p ) = -----------------------------
2
- [EQ 8.206]
(1 + X + X )
where
X = C w ( p – p ref ) [EQ 8.207]
and C w is the water compressibility which is defined as:
1 ⎛d
C w = – ------
B
---------w⎞ [EQ 8.208]
B w ⎝ dp ⎠ p = p ref
where
Y = C μw ( p – p ref ) [EQ 8.210]
and C μw is the water viscosibility which is defined as:
dμ w⎞
1- ⎛ --------
C μw = – ----- - [EQ 8.211]
μ w ⎝ dp ⎠ p = p ref
that also has units of barsa–1 (metric), psia –1 (field) or atmos–1 (lab).
An additional complication can arise if the effect of salt (brine) concentration is to be modeled
(ECLIPSE BlackOil and GI option only). See the ECLIPSE keyword PVTWSALT for further
details in this case.
Note Correlations for generating default values for all the above correlations can be found
in [Ref. 33].
Hint Since the API Tracking functionality in ECLIPSE models the mixing of live oil
properties the usual procedure in PVTi is to use the API Tracking export facility to
write out a series of tables containing the PVTO (live oil) and PVDG (dead gas).
You can then specify whether they want a gas table for each sample. Normally the API Tracking
option in ECLIPSE only uses a different oil table for each sample and just a single gas table is
used, as explained in the next section, and so this is the default. If this default option is used then
PVTi exports a single gas table for the sample with the median gas density at surface conditions.
The separator and units to use for the export can be specified in the normal way. You can also
specify a few other straightforward options such as whether to write to full double precision and
whether to plot the tables.
When the export is performed PVTi orders the tables so that the oil surface densities increase
monotonically with table number, which is a requirement of ECLIPSE. Each table also has
comments associated with it specifying which sample the table is linked to and what the oil
surface density is for this particular table.
Note In general the quantities TREFT, PREFT, DREFT are not the same as the related
quantities of reference temperature, TREF and references density, DREF used in
PVTi. Please see
Workflow
1 To access the Thermal Support module select Edit | Samples | ECLIPSE Thermal
Support | Export for ECLIPSE Thermal.... Alternatively right- click on a fluid sample
of interest and select Export ECLIPSE Thermal Model....
2 In the panel type in the maximum and minimum values of the pressure in your reservoir
and also the maximum and minimum temperature.
3 Select the sample you wish to export the model for by using the drop-down menu.
4 You need to decide if you want to estimate K-values in ECLIPSE Thermal using
Crookston’s equation or using Wilson’s formula. Wilson’s formula gives much less
accurate results than Crookston’s formula but if you have time constraints then much less
work is involved in exporting this information from PVTi. If you want to use Wilson’s
formula then your work is finished and you do not need to go though steps 5-12. Simply
select which units you want to export and click OK.
Hint If you have time and you know the rough limits on temperature and pressure within
your reservoir then Crookston’s equation will, in general, give you much more
accurate estimates of K-values than Wilson’s formula.
Note The default file extension when exporting for ECLIPSE Thermal is .PVOT. The save
panel assumes you want to call your exported file filename.PVOT where
filename.PVI is the name of your PVTi project file.
See the next section for some hints on performing the optimization of the Crookston
coefficients.
–D
K ( P, T ) = ⎛⎝ A + B
--- + CP⎞⎠ exp ⎛⎝ ------------⎞⎠ [EQ 8.212]
P T–E
Introduction
• "Recommended PVT analysis for oil reservoirs" on page 372.
• "Recommended PVT analysis for gas condensate reservoirs" on page 377.
• "Consistency tests and correlations" on page 381.
Oil reservoirs
Oil reservoirs have been divided into three classes in the following discussion:
Under-saturated reservoirs
Refers to those reservoirs that are under-saturated not only at the time of their discovery but also
throughout their development.
Note This class does not include reservoirs whose pressures fall below the bubble point at
some time during their development, or reservoirs that are subject to gas injection.
Under-saturated reservoirs can be represented by two-phase (oil and water) simulation models,
usually using the black oil formulation. For further information see "Under-saturated
reservoirs" on page 374.
Note All reservoirs with initial gas-caps are regarded as saturated (but may not be
adequately represented by the black oil formulation).
Note Volatile oil usually refers to a reservoir fluid whose critical temperature is only slightly
above the reservoir temperature under initial conditions.
Oil reservoirs
Almost all under-saturated reservoirs, and most saturated oil reservoirs where the reservoir fluid
has a critical temperature far above the reservoir temperature, do not normally require the
complexity of compositional simulation, and a black oil formulation is normally of sufficient
accuracy.
Note Such a formulation does not allow the stock tank oil density or gas density to vary with
time.
1 γ g ⋅ GOR ⋅ 0.00122
141.5 - + ---------------------------------------------
B o ( p b ) = -------- ---------------------------------
ρ res ( 131.5 + API )
- [EQ 8.213]
5.6146
where ρ res is the reservoir density in gm/cc, and the GOR in scf/rb.
4 You can find the oil formation volume factor at pressures higher than the saturation pressure
using the compressibility factor in the PVT report.
5 If the field has only been tested and there is no data on producing GOR and gas gravity then
often a separator test for appropriate separator conditions are reported. This supplies a
formation volume factor for bubble point oil at these conditions of separation. You can use
the reported compressibilities to derive formation volume factors at higher pressures.
If a separation test is reported and the field is under production, then you should rationalize
the formation volume factors derived from point 2 with the separator test.
6 If no separation test was performed by the laboratory and there is no production data, then
you can either use a suitable set of correlations, or, if the sample is a re-combination, it may
be possible to use the re-combination data to derive data at separation conditions using
PVTi, providing it is possible to characterize the fluid.
Note It is important to avoid using what is often referred to as the formation volume factor
(it is actually a relative volume) during a differential liberation experiment. The oil
remaining at the end of this experiment has a completely different composition from
stock tank oil and this formation volume factor is usually higher (by as much as 10%)
than the correct value.
For the simulation of undersaturated reservoirs this differential liberation experiment
is usually superfluous, although laboratories often perform it anyway.
Note Again it must be stressed that you should not use differential liberation data alone.
What is often referred to as the formation volume factor during the liberation makes
no reference to separator conditions.
To be completely accurate the laboratory should perform a separation test on the fluid resulting
from each pressure stage of the differential liberation experiment.
In practice this is rarely done, and the volumetric data is determined from the single separation
test (at the bubble point pressure), and the relative volumes and gas volumes of the differential
liberation experiment.
The correct oil formation volume factor is this relative volume, divided by the relative volume
of the bubble point fluid (during the same differential liberation), multiplied by the correct oil
formation volume factor of the bubble point fluid for the separator conditions operating. This
derivation, of course, uses the black oil approximation that the stock tank oil density is
unchanging.
Solution gas-oil ratio is derived in a similar way, the value at the bubble point pressure being
exactly that derived from the separator test ( R s ( p b ) ).
Simple material balance of gas and oil shows that the correct gas-oil ratio ( Rs ) at a pressure
below the bubble point is given by:
R s = R s ( p b ) – ( R s' ( p b ) – R s' )B o' ( p b ) ⁄ B o ( p b ) [EQ 8.214]
where Rs' and Bo' refer respectively to the reported solution gas-oil ratio and formation volume
factor from the differential liberation experiment.
Introduction
• "Blackoil or compositional simulation" on page 377.
• "Data analysis" on page 378.
• "Compositional grading with depth" on page 378.
• "Data analysis" on page 378.
• "Equation of state model" on page 379.
Note The results of such modeling are to over-estimate the condensate recovery from the
field by a few percent for full and partial pressure maintenance. Low pressure
recycling will probably be very poorly modeled.
Data analysis
A common failing when analyzing gas condensate reservoirs is to attempt to establish an
equation of state representation without thorough analysis of the data on which it is to be based.
A thorough comparison of samples, analysis of inter- and intra- sample consistency, and
clarification of the aims of establishing an equation of state model, allows you to determine
which data is worth trying to match, what components to use, and how best to establish a
matched model.
You should closely scrutinize all the available samples. Some of the questions that should be
answered are as follows:
1 Are the samples properly constituted (examine recombination data and compare with
correlations)?
2 Have experiments been performed consistently and accurately (examine compositional
material balances)?
3 Are other experimental data such as Z -factors, implied K -values and liquid densities
reasonable (correlations and tests of consistency)?
4 What are the characteristics of the heavy fractions (fingerprint plots)?
During this analysis it may seem reasonable to establish alternative (modified) data by adjusting
sample data to achieve material balance consistency and a reasonable fluid description. This can
be a useful exercise but a good appreciation of the main sources of inaccuracy is required for it
to be successful.
If the sample is a recombination, the report usually contains recombination data that can be used
for calibrating correlations, to obtain an analysis of separation using the COMB section of PVTi.
For further information see "COMB - Compositional Material Balance" on page 112 and
"COMB section keywords" on page 159.
It is possible to obtain gas and condensate recovery estimates for different abandonment
pressures by combining this with constant composition data above the dew point, and CVD data
below the dew point.
After you have analyzed each sample individually you should examine all the samples together
by, for example, plotting dew point pressures and condensate-gas ratios against sampling depth,
comparing compositions, fingerprint plots etc.
Before rejecting any particular sample that appears to be anomalous, you should investigate the
source of anomalies. Such an investigation may be useful in making recommendations
concerning sampling techniques or laboratory analysis.
Note A cubic equation of state can only supply approximations to the behavior of fluids.
This must be borne in mind when matching laboratory data, but especially when
extrapolating outside the area of data control. By the same token, laboratory analyses
of samples only normally supply an approximation to the sample behavior and the
sample itself is only an approximation to the reservoir fluid.
Introduction
During initial data analysis a number of tests are available in the COMB (COmpositional
Material Balance) section of the program for examining the consistency and quality of a
particular laboratory sample.
For further information see "COMB - Compositional Material Balance" on page 112 and
"COMB section keywords" on page 159.
This section contains information on the following:
• "Compositional material balances during the CVD experiment" on page 381.
• "Fluid density correlations" on page 381.
• "Viscosity calculations" on page 382.
• "Separator calculations" on page 382.
• "K-value checks" on page 382.
• "Recovery calculations" on page 383.
Viscosity calculations
Laboratories generally do not have the equipment for measuring gas viscosities and usually
calculate viscosities based on the Carr, Kobaysahi and Burrows correlation. Other correlations
exist, such as the "Lohrenz, Bray and Clark" on page 330 method and the "Pedersen et al." on
page 331 corresponding states technique. These different correlations often yield markedly
different predicted viscosities.
Ideally, a consistent correlation should be used throughout the reservoir engineering analysis,
from the well test analysis used to derive permeabilities to the reservoir simulation. ECLIPSE
accepts a set of critical volumes (entered with VCRITVIS or ZCRITVIS), which are only used
for viscosity calculations, in order that significant changes in Vc values used to match viscosity
correlations do not effect the remainder of the simulation.
Separator calculations
Standing evolved a method for calculating K -values for separators, based on the HCH method.
For each component, a b -factor is calculated at the separator temperature.
HCH suggested that this b -factor is based on pure component normal boiling points and critical
temperatures.
Standing suggested a slightly different set but this probably only applied to the particular oil he
was considering. b -factors obtained by the HCH method may be used in PVTi to generate
separator K -values.
An alternative is to use those derived from recombination data, though these should first be
checked against those obtained using the HCH method.
From the K -values derived in this way, PVTi can calculate the condensate and oil recoveries for
a given separator system. Separator liquid densities are based on the Alani-Kennedy correlation
(see "Alani-Kennedy liquid densities" on page 311)and stock tank densities are calculated using
Amagat's law (addition of specific volumes).
Recovery calculations
From input CCE, CVD data and the separation system, PVTi can calculate gas and condensate
recovery factors from any supplied initial reservoir pressure to any supplied abandonment
pressure, for a depletion scenario.
Introduction
Fluid properties estimation is used in situations where a full equation of state matching across
numerous samples is not available.
Typically, there is only one sample and limited information about the fluid behavior. From this,
fluid properties estimation simulates typical experiments and provide analysis of the fluid ahead
of full laboratory experiments.
Minimum information
In fluid property estimation the following data provide sufficient information to construct a
complete set of example experiments:
• Weight percentage of each component (or mole fractions)
• The mole weight of the plus fraction
• A salutatorian pressure (bubble or dew point)
• The maximum pressure to use when constructing a depletion experiment (the default is to
use the saturation pressure.
Technical Information
The default settings for the PVTi fluid properties estimation are in Table 8.12.
Property Value
Equation of state Peng-Robinson 3-parameter (corrected)
Library for component properties Katz-Firoozabadi
Correlations for plus-fraction properties Kesler-Lee
Viscosities Lohrenz-Bray-Clark
The equation-of-state and the viscosity model can be changed using PVTi: Edit | Fluid Model
| Equation of State...
The component library in use can be selected in PVTi: Utilities | Program | Options...
The correlation used for the plus-fraction can be changed in PVTi: Edit | Fluid Model |
Components...
Introduction
The goal of PVT analysis is to provide a tuned Equation of State that can model the reservoir
fluid in simulations. The selection of regression parameters is crucial in determining the quality
of the tuned fluid model.
This chapter discusses the selection of regression variables in PVT analysis and consolidates
traditional variable selection, based on physical principles ("Physical selection of regression
parameters" on page 386), with variable selection based on mathematical analysis of the
problem ("Mathematical analysis of the regression problem" on page 388). "A consistent
methodology that can be applied automatically" on page 389, describes how information in the
preceding two sections are used by the PVTi Quick Fit option.
Fundamentally, the fitting process twists
the standard cubic-equation of state to fit the supplied data. Across the temperature and pressure
ranges of the observations the fitted model may be useful, at very different temperatures and
pressures the model may not be representative of the fluid behavior at all.
Note Always use data at reservoir conditions, such as depletion experiments, and at surface
conditions, such as separator tests, when fitting a fluid model that is to be used in
reservoir simulation.
In "Weighting observations for regression" on page 390, there is a discussion on weighting the
observation data to get the desired fit.
This chapter does not contain a mathematical description of the Levenberg-Marquardt
algorithm, that can be found in "Weighting observations for regression" on page 390.
Critical properties
For pure library components, particularly the non-hydrocarbons and lighter hydrocarbons, the
acentric factor (ω) and critical temperature (Tc) and pressure (Pc) are well known. It is generally
not justifiable to tune these properties.
Hint It is very unlikely that the Pc, Tc and ω of H2O, H2S, N2, CO2, CO, C1, C2, C3, IC4,
NC4, IC5, NC5, and C6 differ from the library values. So do not use them in
regression.
Hint The Pc, Tc, and ω of the plus fraction usually come from characterization. This means
that they are not well determined, making them ideal candidates for tuning.
Pedersen et al. have suggested that the molecular weights of the characterized components
might be a good choice of regression variable as this varies Pc, Tc and ω consistently, in line
with the chosen characterization method, see [Ref. 59].
Caution Careless tuning of binary interaction coefficients often leads to a fluid model that
has convergence problems in compositional simulations.
The basic cubic equation-of-state uses the acentric factor to consider slight deviations in
molecule shape; the assumption is that all molecules are close to being spherical. The binary
interaction coefficient between light components and the heaviest hydrocarbon can be used to
compensate for the non-sphericity of the heavy hydrocarbon molecules.
ΩA and ΩB coefficients
The default values of ΩA and ΩB are based on the following assumptions:
1 The critical isotherm of a component has a slope of zero and an inflection point at the
critical point. In layman’s terms this means that the component is pure.
2 That the determination carried out, for the limited number of pure components available,
can be extrapolated to all heavier and lighter components.
3 The values were valid for pure-component density and vapor pressure below critical
temperature.
In developing his equation of state, Tareq Ahmed explained why these basic assumptions are
violated in many fluid models, see [Ref. 61]:
1 Pseudo-components and the plus fraction are not pure components, so assumption 1 is
violated by these.
2 The heavier hydrocarbons and light components may not have the same ΩA and ΩB as the
pure components tested. Specifically, assumption 2 might not be true.
Hint The assumptions upon which the default ΩA and ΩB constants are based do not hold
for light components, pseudo-components, and the plus fraction. This makes them
candidates for regression.
Hint Varying the critical volume and Z-factor are equivalent in this case, so it is futile to
vary both at once.
This critical volume or Z-factor is only used to calculate the viscosity, so remember that in a
regression only the viscosity observations themselves are used to tune these parameters.
Caution As only a few observations apply to the tuning of the critical volumes for
viscosity, make sure that there are no more than one or two regression variables of
this type in a regression run. It is much better to put all critical volumes for
viscosity into one regression variable than to vary all of them at once.
There is another option in PVTi that allows the constants used in the fourth order polynomial in
density to vary.
Caution This option is available so that the original Lohrenz-Bray-Clark analysis can be
repeated with new data-sets. It should only be used if you have a large number of
samples and lots of viscosity observations. Fitting too many parameters to too few
observations leads to convergence problems for the fluid in a compositional
simulator.
High sensitivities
The first rule means that interest should focus on the most sensitive parameters, those for which
a small change has a large effect on the fit.
The Hessian tab of the Sensitivity Analysis panel gives an indication of the sensitivity of the
different parameters. The values along the leading diagonal of the Hessian matrix are most
useful, a high value means high sensitivity.
Hint By removing insensitive regression variables from the set, the problem becomes less
non-linear and thus more solvable.
Kai Liu has reported success with an automatic regression technique based solely on selecting
the eight most sensitive regression parameters, see [Ref. 61].
Hint Amalgamating strongly correlated regression variables of the same type (for example
Tc) in a single regression variable or removing one of the two strongly correlated
variables makes the problem less non-linear and thus more solvable.
Note This method takes no advantage of amalgamating regression variables of the same type
into a single variable.
Hint Use the saturation pressure measurement from the depletion experiment to create a
dew-point or bubble-point experiment. Weight the observation of saturation pressure
very high (1000+).
Hint Weight the separator densities or GOR highly (500+) to guarantee that the right amount
of oil and gas are produced.
Hint Weight the worst-fit observation the most and successive observations less and less
through to the best fit observations that retain their original weighting of 1. This guides
the regression to take more account of the worst parts of the fit.
Caution Be aware of the errors in the observed data and do not be tempted to over-fit.
Models that are over-fitted cause convergence problems in compositional
simulations.
Introduction
The solid precipitation model in PVTi handles waxes as a solid solution and asphaltenes as an
additional liquid phase.
The wax model assumes that paraffins and naphthalenes can precipitate as a solid solution. In
the asphaltene model it is assumed that the heaviest aromatic component will form an
asphaltenic liquid phase.
To use the solid precipitation model, the paraffin, naphthalene and aromatic fractions for each
carbon number are needed. These can be entered by hand, if known, or generated from the
correlations described in "The PNA distribution of heavy components" on page 394 and
"Critical properties of PNA species" on page 395.
The theory behind the wax model is explained in "Wax precipitation" on page 395.
3660
A% = 430v + ------------ ( for v < 0 )
MW
3660
A% = 670v + ------------ ( for v > 0 )
MW
10000
R = 820w + --------------- ( for w > 0 )
MW
10600
R = 1440w + --------------- ( for w < 0 )
MW
N% = R – A%
n = ⎛1 + 2I⎞
-------------- [EQ 8.216]
⎝ 1–I⎠
Hint If a full PNA specification of the sample is available, the mole fractions can be used in
the normal way, over-riding the default estimates from this splitting procedure. First
enter the single carbon number (SCN) fluid; perform the PNA split, then type in the
measured mole fractions.
Note The critical properties of the PNA species can be tuned just as for any other
components.
In addition to the usual set of critical properties, the melting points of the PNA species are also
determined, for use in the wax precipitation.
For paraffins, the melting point is given by the correlation of Won (1986), see [Ref. 51]:
For naphthalenic and aromatic species, the melting point is given by (Pan, Firoozabadi and
Fotland, 1997 - see [Ref. 52]):
f – 0.00855MW
T = 333.45 – 419e [EQ 8.219]
Wax precipitation
In PVTi it is considered that the paraffinic and naphthalenic species of components heavier than
C15 can form waxes and that the heaviest aromatic component forms the asphaltenic liquid
phase.
Chemically, the waxes that drop out of hydrocarbon fluids at lower temperatures, are known to
contain paraffins and, to a lesser extent, naphthalenes.
In PVTi the method of Pedersen et al. (see [Ref. 55]) has been improved to use the paraffin and
naphthalene components, rather than a “wax forming component”. This means that the critical
properties correlations for the paraffins and naphthalenes, as given by the correlations in
"Critical properties of PNA species" on page 395, are used in determining the solid
precipitation.
f f f f
Δh T Δc Pi ⎛ T Δc Pi ⎛ T i⎞
f pure i ( P,T ) = f pure i ( P,T ) × exp ---------if ⎛⎝ 1 – -----i⎞⎠ – ----------
- 1 – -----i⎞ – ----------
S L
R ⎝ ⎠ R
- ln -----
⎝ T⎠ [EQ 8.220]
RT i T T
where Δh fi is the enthalpy of fusion for component i; ΔcPi is its heat capacity of fusion; and Tfi
is its melting point temperature.
Enthalpy of fusion
The enthalpy of fusion for paraffins comes from the correlation by Won [Ref. 51]:
f f
Δh i = 0.1426MW i T i [EQ 8.221]
For napthalenes the correlation of Lira-Galeana et. al. (1996) [Ref. 56] was used:
f f
Δh i = 0.0527MW i T i [EQ 8.222]
For aromatics the correlation of Pan et al. [Ref. 52] was used:
f f
Δh i = 11.2T i [EQ 8.223]
θ
⎛ 2V 1ci⁄ 6 V 1cj⁄ 6 ⎞
k ij = A 1.0 – ⎜ -----------------------------
-⎟ [EQ 8.225]
⎝ V 1ci⁄ 3 + V 1cj⁄ 3⎠
Discussion
PVTi uses a consistent, single fluid model for all calculations, so the critical properties used in
matching, say, the differential liberation experiment, are used in calculating the wax appearance
temperature. Consequently all the observations, including those of the solid phase, can be
regressed at once, leading to a more consistent and complete single description for use in both
reservoir and process simulations.
Introduction
Oil-based muds are in widespread use and often contaminate PVT samples taken at the well-
site.
PVTi offers two methods for cleaning oil-based muds, "Removing oil-based mud contamination
by skimming" on page 398 describes a method that can be used when the composition of the
contaminant is not known. If the contaminant composition is known, "Removing oil-based mud
contamination by subtraction" on page 399 provides an accurate method for removing the
contaminant.
The methods used in PVTi are based on the work of Gozalpour et al, see [Ref. 53].
Caution It is important that the sample being cleaned has heavier components than the oil-
based mud, so that there is an uncontaminated point to use in constructing the
straight line.
If the composition of the oil-based mud is known, it is better to use the subtraction method
described in "Removing oil-based mud contamination by subtraction" on page 399.
Note The restriction that the oil-based mud contain components C10 - C23 does not apply
here. The mud can be of any composition.
Note The assumption that there are uncontaminated components, heavier than the oil-based
mud components, is not necessary for this method.
Hint If the composition is known, the subtraction method is a better choice than the
skimming method.
Introduction
PVTi offers both sample mixing and sample recombination. There is a subtle difference between
mixing and recombining. In mixing, you provide a mole fraction or gas-oil ratio for the mix;
these are used and the mixed sample is created directly. In recombining, you enter a target gas-
oil ratio for the mixture. In that case the proportion of mixing is determined iteratively using the
scheme outlined in "Recombination" on page 400.
Mixing
Mixing is an addition of the two samples. If the mole fraction of the second sample is specified,
the mixture is a weighted sum. If gas-oil ratio (GOR) is specified, this value is first translated
to a mole fraction at the temperature and pressure requested.
The conversion from GOR to mole fraction (MF) is as follows:
n -
Assuming n moles of Sample 2 are mixed with 1 mole of Sample 1, then F = -----------
1+n
RT
F -⎞ ------------------
The GOR is then given by GOR = ⎛⎝ -----------
1 – F⎠ P std V oil
std
- where R is the universal gas constant, Tstd and
Pstd are standard temperature and pressure respectively and Voil is the oil volume at the
temperature and pressure you specified. If you do not specify the pressure, it is taken as the
saturation pressure of the sample
Recombination
In recombination, the aim is to create a mixture with the stock tank gas-oil ratio specified.
In this case, the temperature and pressure you enter are used as the first stage in a two-stage
separator. The second stage is always at standard conditions. The GOR you enter is the target
GOR for the mixture. This is used as an initial guess and the fluids are mixed in the usual way
(see "Mixing" on page 400). When this mixture is passed through the separators, a stock-tank
GOR is calculated. A regression is performed, adjusting the mixing combination until the stock-
tank GOR is equal to the target value you requested.
The mix applied at each iteration is recorded in the log window.
Introduction
In PVTi and ECLIPSE Compositional we deal with an isothermal flash. This means that, for a
particular cell in a simulation, we know the composition of the fluid summed over all the phases
and the pressure and temperature. We try to minimize the Gibbs Free Energy in order to
determine how each component splits across the different phases present. In effect we try to find
the K-values, which are the unknown variables.
ECLIPSE Compositional assumes that the temperature of each cell stays constant over time,
even if there is a distribution of temperatures across the cells to start with. It is assumed that if
a particular drop of fluid moves from one cell with temperature, T1 to another cell with
temperature T2 that the fluid takes on the temperature T2. When the Thermal option is used in
ECLIPSE Compositional this is not the case. Thermal diffusion is allowed to take place over
time, which means that the temperature in each cell is a free parameter that needs to be
determined using a different type of flash; one at constant energy. We therefore know the
pressure, functional form of the K-Values (K=K(P,T)), and the total energy of the system but the
temperature is unknown.
The long-term goal of this ECLIPSE Thermal module is to use PVTi’s powerful Equation of
State (EoS) functionality to provide extensive support for the Thermal option in the ECLIPSE
Compositional simulator. The first step was to introduce functionality into PVTi 2003A to
calculate an optimal K-value functional relation. This can then be used within ECLIPSE
Thermal to perform these flashes at constant energy.
For PVTi 2004A the module has been extended to write out a whole file containing a series of
keywords that can be used as the PVT input for an ECLIPSE Thermal simulation. In addition
to the standard keywords written out for an ECLIPSE Compositional run, for example acentric
factors (ACF), critical temperatures (TCRIT) etc., PVTi writes out additional keywords that
enables ECLIPSE Thermal to calculate accurate values for K-values oil/gas densities and oil/
gas phase enthalpies.
• "K-Values" on page 401, "Oil Density" on page 403,
• "Gas phase density" on page 405,
• and "Enthalpy" on page 405 explain how the relevant keywords are used in ECLIPSE
Thermal to calculate the appropriate quantities.
For a brief non-technical summary of these keywords see "Outline of keywords for ECLIPSE
Thermal" on page 367.
K-Values
For PVTi 2004A you have the option of either using Crookston’s equation (the KVCR keyword)
or Wilson’s formula (the KVWI keyword).
–D
--- + CP⎞⎠ exp ⎛⎝ ------------⎞⎠
K ( P, T ) = ⎛⎝ A + B [EQ 8.226]
P T–E
If a simplified version of equation [EQ 8.226] is used, given by setting A, C, E=0 that is:
B
K ( P, T ) = ⎛⎝ ---⎞⎠ exp ⎛⎝ – D
----⎞⎠ [EQ 8.227]
P T
then B and D can be calculated by using at least two experimentally determined observations.
Alternatively, Wilson’s formula can be used to estimate K-Values, which is given by the
formula:
P T
K ( P, T ) = ⎛⎝ -----c⎞⎠ exp 5.372697 ( 1 + ω ) ⎛⎝ 1 – ----c-⎞⎠ [EQ 8.228]
P T
where P, T are the pressure and temperature and K , P c , T c and ω are the K-value, critical
pressure, critical temperature and acentric factor respectively for a particular component. In fact
Crookston’s equation is a generalized form of Wilson’s formula, which can be obtained by
setting:
A = 0
B = Pc
C= 0
T
D = 5.372697T ( 1 + ω ) ⎛⎝ ----c- – 1⎞⎠
T
E = 0
in equation [EQ 8.226].
There are also tabulated values for B and D for certain components. The central problem we face
though is that K-values for a particular component are not functions of just temperature and
pressure, as Wilson’s formula would suggest, but also of the types and quantities of other
substances present. Consider the following. If we have a fluid of say 50% methane and 50%
decane then each component has its own K-value at a particular temperature and pressure. If
we now add, say toluene, to the mixture then the K-values of methane and decane changes and
we also have the K-value for the new component. This illustrates that changing the composition
of a fluid clearly effects the K-values of respective components.
Also, for heavier hydrocarbon components (anything above about C7), the K-values are an
increasing monotonic function of pressure. It is clear that Wilson's formula is completely
inadequate here as K=Pc/P is clearly a decreasing monotonic function of pressure. The full
Crookston expression can cope with such components by using the A and C terms.
PVTi can provide valuable assistance here because it provides its own experimental data using
the flash; as we know that the flash accurately reproduces experimental observations. The
ECLIPSE Thermal export facility then calculates the values of a chosen set of Crookston
coefficients so that the correlation best approximates a set of flashes performed by PVTi in the
pressure and temperature range you specified.
Oil Density
PVTi writes out the keywords PREFT, TREFT, DREFT, CREF and THERMEX1, which
ECLIPSE Thermal uses internally to calculate the oil density of the fluid at a specified pressure
and temperature.
Algorithm
ECLIPSE Thermal calculates the molar oil density b o of the fluid using the following set of
equations:
bo is given by:
where the molar specific volume Voil of the oil phase is calculated using Amagat’s law of partial
volumes:
Nc
c c
V oil = ∑ x Voil [EQ 8.230]
c=1
and xc is the mole fraction of component c in the oil phase. The component oil phase volume
c
V oil is given by:
c
Voil = MW
c
------------ [EQ 8.231]
c
ρ
where MW c is the molecular weight of component c given by the MW keyword, and the
component oil phase density ρ c is given by:
c
ρ
c ref
ρ = ----------------------------------------------------------------------------------------- [EQ 8.232]
( 1 + C T1 ( T – T ref ) ) ( 1 – C p ( P – P ref ) )
The reference density ρ c is defined by ECLIPSE Thermal’s DREF keyword (or PVTi’s
ref
DREFT keyword); the standard temperature and pressure, T ref and Pref are defined by keywords
TREF (TREFT in PVTi) and PREF (PREFT in PVTi); the thermal expansion coefficient C T1 is
defined by the THERMEX1 keyword; and the component isothermal compressibility C P is
defined by the CREF keyword.
The quantities C T1 and C p are defined by the usual thermodynamic relations:
1 ∂V
C T1 = --- ------ [EQ 8.233]
V ∂T P = PREF
Remarks
It is important to note that PVTi’s DREF, TREF keywords are, in general, not the same as
ECLIPSE Thermal’s DREF, TREF keywords.
Algorithm
The molar density of the gas phase bg is given by:
Nc
c c
V gas = ∑ y Vgas [EQ 8.236]
c=1
assuming no water is present. yc is the mole fraction of component c in the gas phase. Each
hydrocarbon component gas phase molar volume V gas c is obtained from the gas law,
c c
PV gas = Z RT [EQ 8.237]
Remarks
The algorithm outlined above assumes that Z c does not change as a function of temperature.
c
Z is calculated by PVTi for each component by performing an EoS flash on a fluid consisting
purely of the relevant component. If a single gas phase is found then Z c is set to the value of the
Z-factor returned by PVTi’s EoS code. If a liquid phase is returned then Z c is set to ECLIPSE
Thermal’s default value of 0.96.
Enthalpy
• "Algorithm" on page 405 outlines the three different approaches ECLIPSE Thermal can
take in calculating oil and gas phase enthalpies.
• "Remarks" on page 405 then explains which of these approaches PVTi supports.
Algorithm
ECLIPSE Thermal calculates the enthalpy of the oil phase using a mole fraction weighted
average of the component enthalpies:
Nc
c c c
H oil = ∑x ⋅ MW ⋅ H
oil
[EQ 8.238]
c=1
Nc
c c c
H gas = ∑y ⋅ MW ⋅ H gas [EQ 8.239]
c=1
Specific heats
The specific heats C coil , C cgas of a component c in the oil and gas phases are given by
c c c
C oil = C oil1 + C oil2 ⋅ ( T – T st )
c c c
C gas = C gas1 + C gas2 ⋅ ( T – T st )
where the oil phase coefficients are defined by keywords SPECHA and SPECHB and the gas
phase coefficients are defined by keywords SPECHG and SPECHH.
T st is the standard temperature defined by STCOND. The specific heats must always be positive.
If oil phase specific heats are defined, the molar enthalpy of component c in the oil phase is:
c c c 1 c 2
H oil = ∫ Coil dT = Coil ⋅ ( T – Tst ) + --2- ⋅ Coil ⋅ ( T – Tst )
1
2
[EQ 8.240]
T st
For temperatures above the critical temperature, T > T ccrit , the oil phase molar enthalpy is set
equal to the gas phase molar enthalpy.
If gas phase specific heats are defined, the molar enthalpy of component c in the gas phase is
where Hcvaps is the heat of vaporization at the standard temperature, defined by HEATVAPS.
If a solid phase is simulated, the molar enthalpy of component c in the solid phase is
c c c 1 c 2
H sol = ∫ Csol dT = Csol ⋅ ( T – Tst ) + --2- ⋅ Csol ⋅ ( T – Tst )
1
2
[EQ 8.242]
T st
Heat of vaporization
A temperature dependent heat of vaporization ΔHc ( T ) can be specified instead of either the oil
or gas specific heat.
If the specific heat of a component c in the oil phase is not specified, then the oil phase enthalpy
is calculated from the gas component enthalpy by
c c c
H oil = H gas – ΔH ( T ) [EQ 8.243]
Alternatively, if the specific heat of a component c in the gas phase is not specified, then the
gas phase enthalpy is calculated from the oil component enthalpy
c c c
H gas = H oil + ΔH (T) [EQ 8.244]
c c B T < T crit
c
ΔH ( T ) = A ⋅ ( 1 – T ⁄ Tcrit ) [EQ 8.245]
c c
ΔH ( T ) = 0 T ≥ T crit
[EQ 8.246]
where the constant A is defined by the HEATVAP keyword, the exponent B is defined by the
HEATVAPE keyword, and T crit c is the component critical temperature defined by the TCRIT
keyword.
Heats of vaporization are usually obtained at the normal boiling point T nb . The constant A ,
defined by the HEATVAP keyword can be found by inverting :
c
ΔH ( T nb )
A = ---------------------------------------
B
- [EQ 8.247]
( 1 – T nb ⁄ T crit )
where ΔHc ( T nb ) is the heat of vaporization at the normal boiling point. The exponent B for each
component, entered with keyword HEATVAPE, is usually set to a value in the range 0<B<=1.
Units
General information
This option allows the current unit and temperature conventions to be changed. The options
available are:
• Metric unit system
• Field units
• Laboratory units
• PVT-metric units
The temperature conventions are:
• Degrees Kelvin
• Degrees Celsius
• Degrees Rankine
• Degrees Fahrenheit
Any of the unit conventions are compatible with any of the temperature options. The unit
conventions may be changed at any point in a session, existing quantities being converted to the
new units. PVT-metric units are the same as metric, except in that pressure is measured in
atmospheres. In addition, it is possible to enter mole compositions as fractions (of unity) or
percentages (up to 100%), and pressures can be specified in absolute or gauge units.
Reservoir volume rm
3 rb rcc rm
3
Formation volume rm ⁄ sm
3 3 rb ⁄ mscf rcc ⁄ scc rm ⁄ sm
3 3
factor (liquid)
Gas-oil ratio sm ⁄ sm
3 3 mscf ⁄ stb scc ⁄ scc sm ⁄ sm
3 3
Oil-gas ratio sm ⁄ sm
3 3 stb ⁄ mscf scc ⁄ scc sm ⁄ sm
3 3
Volume m
3
ft
3 cc m
3
Specific volume 3
m ⁄ kg
-ml 3
ft ⁄ lb
-ml
cc ⁄ gm
-ml 3
m ⁄ kg
-ml
Energy kJ btu J kJ
Enthalpy kJ ⁄ kg
-ml
btu ⁄ lb
-ml
J ⁄ gm
-ml
kJ ⁄ kg
-ml
Specific heat kJ ⁄ kg
– ml
⁄K btu ⁄ lb
– ml
⁄R J ⁄ gm
– ml
⁄K J ⁄ kg
– ml
⁄K
Transmissibility 3
cPm ⁄ d ⁄ bar cPrb ⁄ d ⁄ psi cPcc ⁄ h ⁄ a cPm ⁄ d ⁄ a
3
1ft
3
= 0.02831685 m 3 = 0.1781076 bbl
1lb ⁄ ft
3
= 16.0184600 kg ⁄ m 3 = 0.01601846 gm ⁄ cc
Symbols
p Pressure of mixture.
p ri p
Reduced pressure of component i : p ri = ------i .
p ci
T Temperature of mixture.
Tc Critical temperature of component i .
T ri T
Reduced temperature of component i : T ri = ------i- .
T ci
Paper No. 15C Ohio, presented at the AIChE 65th National Meeting, Cleveland. May 4-7, (1969).
SPE paper no. 10512, presented at the Sixth SPE Symposium on Reservoir Simulation, New
Orleans , January 31st - February 3rd, (1982).
Whitson, C.H. “Topics on: Phase Behaviour and Flow of Petroleum Reservoir Fluids”, [Ref. 4]
Lohrenz, J., “Calculating Viscosity of Reservoir Fluids from their Composition”, [Ref. 7]
Bray, B.G., and J. Pet. Tech. Page 1171 , (1964); Trans., AIME,. 231.
Clark, C.R.
Crowe, C.M., and “Convergence Promotion in the Simulation of Chemical Processes - The General Dominant
Nishio, M. Eigenvalue Method”, [Ref. 8]
Cavett, R.H. “Physical Data for Distillation Calculations - Vapor-Liquid Calculations”, [Ref. 11]
Part 4: Compressibility Factors and Equation of State, Pet. Refiner , Page 173-179 , April, (1958)
Michelsen, M.L. “The Isothermal Flash Problem. Part I. Stability”, [Ref. 15]
Bashbush, J.L. “A Method to determine K-values from Laboratory Data and its Applications”, [Ref. 18]
Alani, G.H., and “Volumes of Liquid Hydrocarbons at High Temperatures and Pressures”, [Ref. 20]
Kennedy, H.T. Page Trans. AIME, 219, Page 288-292 , (1960)
Dranchuk, P.M., “Computer Calculations of Natural Gas Compressibility Factors using the Standing and Katz
Purvis, R.A., and Correlation”, [Ref. 21]
Robinson, D.B. Institute of Petroleum Technical Series, No. IP 74-008, (1974).
Søreide, I., “Procedures for Reservoir Fluid Characterisation using an Equation of State Model”, [Ref. 26]
Reffstrup, J., and , Report RE 88-8, Laboratory for Energetics, Danish Tech. Univ., Lyngby.
Whitson, C.H.
Coats, K.H., and “Application of a Regression-Based EOS PVT Program to Laboratory Data”, [Ref. 27]
Smart, G.T. SPE Res. Eng. Page 277-299 , May (1986)
Herning, F., and “Calculation of the Viscosity of Technical Gas Mixtures from the Viscosity of the Individual
Zippener, L. Gases”, [Ref. 28]
Stiel, L.I., and Thodos, “The Viscosity of Nonpolar Gases at Normal Pressures”, [Ref. 29]
G. AIChEJ , Page 611-615 , (1961)
Whitson, C.H., “ C 7+ Characterisation of Related Equilibrium Fluid using the Gamma Distribution”, [Ref. 32]
Anderson, T.F., and Paper in C 7+ Characterisation, Edited by Mansoori G.A., and Chorn L.G., Taylor and Francis, New
Søreide, I. York, (1989).
Søreide, I. “Improved Phase Behaviour Predictions of Petroleum Reservoir Fluids from a Cubic Equation
of State”, [Ref. 34]
Søreide, I., and “Peng-Robinson Predictions for Hydrocarbons, CO2 , N 2 , and H2 S With Pure Water and NaCl -
Whitson, C.H. Brines”, [Ref. 35]
Jensen, F., and “Calculation of First Contact and Multiple Contact Minimum Miscibility Pressures”, [Ref. 38]
Michelsen, M.L. In Situ, 14(1) , Page 1-17, (1990)
Gmehling, J., “Vapor-Liquid Equilibrium Data Collection DECHEMA Chemistry Data Series”, [Ref. 39]
Onken, U., and DECHEMA, Frankfurt/Main, (1981)
Arlt, W.
Sim W.J. and “Prediction of vapor-Liquid Equilibrium of Undefined Mixtures”, [Ref. 43]
Daubert T.E. Ind. Eng. Chem. Process. Des. Dev. 19 , Page 386-393, (1980)
1997
Michelsen, M.L. and “Calculation of critical points from cubic 2-constant equation of state”, [Ref. 46]
Heidemann, R.A. AIChE Journal (Vol 27 No3) , Page 769, (1980)
Nes, K. and Westerns, “Aspects of the Constitution of Mineral Oils”, [Ref. 48]
H. A. van Elsevier, New York, (1951)
Riazi, M. R. and Al- “Physical properties of n-alkanes and n-alkylhydrocarbons: Application to petroleum
Sahhaf, T. A. mixtures”, [Ref. 50]
Ind. Eng. Chem. Res., Vol. 34, Page 4145 -4148, (1995)
Gozalpour, F., Densh, “Predicting reservoir fluid phase and volumetric behaviour from samples contaminated with
A., Tehranio, D. H., oil-based mud”, [Ref. 53]
Todd, A.C. and Tohidi, SPE 56747, (1999)
B.,
Katz, D. “Overview of Phase Behavior in Oil and Gas Production”, [Ref. 54]
Pedersen, K. S. “Prediction of cloud point temperatures and amount of wax precipitation”, [Ref. 55]
Pedersen, K. S., “Wax precipitation from North Sea crude oils. 4. Thermodynamic Modeling”, [Ref. 57]
Skovborg, P. and Energy and Fuels, Vol. 5 (6), Page 924, (1991)
Ronningsen, H. P.
Trebble, M.A. “A preliminary evaluation of two and three phase flash initiation procedures”, [Ref. 62]
Teja, A.S. & Rice, P. “Generalised Corresponding States Method for the Viscosities of Liquid Mixtures”, [Ref. 64]
Hanley, H.J.M., “Equation of the Viscosity and Thermal Conductivity Coefficients of Methane”, [Ref. 65]
McCarty, R.D. & Cryogenics, 15, Page 413,(1975),
Haynes, W.M.
McCarthy, R.D., “A Modified Benedict-Webb-Rubin Equation of State for Methane Using Recent Experimental
Data” [Ref. 66]
GI
definition . . . . . . . . .148
LNAMES keyword . . . . . 177 O
LOG files . . . . . . . . . . . . . 93
nodes . . . . . . . . . . . .205 OBS keyword . . . . . . . . . 234
Lohrenz-Bray-Clark214 to 215, 245
GOR Observations . . . . . . . . . . 234
definition . . . . . . . . .148 OBSIND keyword . . . . . . 235
GRAF . . . . . . . . . . . . . . .148 Oil density. . . . . . . . . . . . 339
GROUP section . . . . . . . .159 M Oil formation volume factor339
keywords . . . . . . . . .159
Martin . . . . . . . . . . . . . . 317 Oil viscosity . . . . . . . . . . 339
Grouping . . . . . . . . . 106, 180
Material balance checks . . . 28 OMEGAA, OMEGAB keywords237
components . . . . . . .206
Maximum step size . . . . . 219 On-line help . . . . . . . . . . 151
MES file. . . . . . . . . . . . . . 93 Opening a project. . . . . . . . 95
OPTIONS keyword . . . . . 238
H Metric units . . . . . . . . . . 409
Michelsen’s stability criterion317 OUTECL3 keyword . . . . . 240
Hoffman-Crump plot . . . . .81 Mixing . . . . . . . . . . . . . . 224 OUTECL3 section . . . . . . 162
Hydrocarbon . . 102, 211, 297 keywords . . . . . . . . . 162
Mnemonics. . . . . . . . . . . 177
Modified form of Peng-Robinson246
Modified Whitson . . . . . . 302
I Modified Whitson Splitting305 P
MODSPEC keyword . . . . 225 PARACHOR keyword . . . 242
Individual observation weights235
MODSYS keyword . . . . . 226 Paraffins . . . . . . . . . . . . . 297
Internal library . . . . . . . . .217
Molar volume correction . 320 PCRIT keyword. . . . . . . . 243
Isomers . . . . . . . . . . . . . .297
Mole fractions. . . . . . . . . 292 PEARCE keyword . . . . . . 244