Coefficients of Cavett  correlation.
The characterization of petroleum fluids is fundamental for the calculation of their thermodynamic properties. Laboratory experiments are able to identify a limited number of pure components present in a sample. All remaining species, the so called “cut”, are characterized by its molecular weight and density. The thermodynamic calculations performed using cubic equations of state require the critical properties and the acentric factor, which are unknown for the petroleum “cut.” In this chapter, different correlations are used to calculate the critical properties and the acentric factor of the “cut” fraction. The performance of the correlations is evaluated through the comparison of a simulated pressure-volume-temperature (PVT) experiment using an equation of state and experimental data of two reservoir fluids.
- correlations of critical properties
- petroleum heavy fractions
- PVT experiments
- phase equilibrium
Petroleum is a complex mixture of several chemical components, mainly hydrocarbons. In addition to the hydrocarbons, it may also contain some inorganic contaminants such as carbon dioxide (CO2), nitrogen (N2) and hydrogen sulfide (H2S).
The physical properties of reservoir fluids are related to the concentration of their components. Some properties such as bubble point pressure, oil formation volume factor, solubility ratio, oil density, gas formation volume factor and gas specific gravity are of particular interest in black oil reservoir engineering studies. These properties are generally obtained in laboratory using reservoir fluids samples. These experiments seek to replicate the isothermal recovery path of an oil field.
Pressure-volume-temperature (PVT) experiments are carried out in liquid mixtures of hydrocarbons and during the pressure reduction steps, dissolved gas is released. There are two types of experiments that simulate the constant temperature depletion of a reservoir fluid, “flash” and “differential”. In a “flash” experiment, the overall composition of the system is kept constant, whereas in a “differential” experiment the gas phase is removed from the system at each pressure step.
These experiments can be simulated using equations of state in order to evaluate different recovery schemes without carrying out one experiment for each possible scenario, especially in the cases of enhanced oil recovery techniques. The correct identification of the species and their concentrations is fundamental for the success of the simulation. However, in any laboratory test, the heavy fraction of the oil, the so-called “cut,” is characterized by its molecular weight and density. In order to simulate the experiment, correlations are necessary to calculate the critical properties and the acentric factor of the “cut” from its molecular weight and density. In this chapter, different correlations are used to determine the critical properties and acentric factor of the heavy fraction of two reservoir fluids, the PVT experiments of these fluids are simulated using an equation of state and the results are compared with laboratory data.
Cubic equations of state are widely used to describe the volumetric properties of pure substances and mixtures in the petroleum industry. Furthermore, these models can be used for equilibrium calculations and to estimate PVT properties, such as bubble point pressure, oil formation volume factor, solubility ratio, oil specific weight, gas formation volume factor and gas density.
The input parameters of the most used equations of state are the acentric factor, ω, critical temperature, Tc, and critical pressure, Pc, of the mixture components. These parameters are tabulated for a large number of chemical compounds, but for the heavy fractions of petroleum fluids, they are determined from correlations. Most correlations are functions of density, γ, molecular mass, M, and/or normal boiling temperature, Tb, of the fractions .
Edmister  proposed a correlation to estimate the acentric factor of pure liquids and petroleum fractions. This correlation is given by:
Cavett  presented equations to estimate the critical temperature and pressure of hydrocarbon fractions. These correlations are a function of the normal boiling point and the API gravity:
Tb: normal boiling temperature [°F].
|1||1.7133693000||0.94120109 × 10−3|
|2||−0.0010834003||−0.30474749 × 10−5|
|3||−0.0089212579||−0.20876110 × 10−4|
|4||0.3889058400 × 10−6||0.15184103 × 10−8|
|5||0.5309492000 × 10−5||0.11047899 × 10−7|
|6||0.3271160000 × 10−7||−0.48271599 × 10−7|
|7||0.13949619 × 10−9|
Kesler and Lee  developed correlations for critical temperature and pressure, molecular weight and acentric factor of oil fractions. These expressions are given by:
for Tbr > 0.8, and
for Tbr ≤ 0.8.
In Eq. (6), Kw is the Watson characterization factor, given by:
Standing  developed the following correlations based on experimental data:
Riazi and Daubert  proposed a simple two-parameter equation to calculate the physical properties of pure components and mixtures of hydrocarbons:
where θ is the property (Tc, Pc, vc or M), vc is the critical volume (ft3/lb) and a to c are constants for the correlation of each property.
|θ||a||b||c||Average deviation (%)||Max deviation (%)|
|M||−4.56730 × 10−5||2.19620||−1.0164||2.6||11.8|
|Pc||−3.12281 × 109||−2.31250||2.3201||3.1||−9.3|
|Vc||−7.52140 × 10−3||0.28960||−0.7666||2.3||−9.1|
Sim and Daubert  determined expressions to calculate the critical pressure and critical temperature of petroleum “cuts,” given by:
Twu  proposed a set of correlations, based on the perturbation-expansion theory with normal paraffins as the reference state, to determine the critical properties of heavy fractions of hydrocarbons. The method is based on the selection of a normal paraffin with normal boiling temperature, TbPi, identical to the normal boiling temperature of the Cn+ fraction. Once the normal paraffin is chosen, the heavy fraction critical properties are calculated through the following steps:
1. Calculate the critical properties of the normal paraffin
|TcPi||0.5332272||0.191017 × 10–3||0.779681 × 10−7||−0.284376 × 10−10||0.959468 × 10−2||0.01|
2. Calculate the heavy petroleum fraction properties from the equations
where TcPi is the critical temperature of the normal paraffin (°R), TcC+ is the critical temperature of the heavy petroleum fraction (°R), PcPi is the critical pressure of the normal paraffin (psia), PcC+ is the critical pressure of the heavy oil fraction (psia), γPi is the density of the normal paraffin, γC+ is the density of the heavy petroleum fraction, vcPi is the critical volume of the normal paraffin (ft3/lbmol) and vcC+ is the critical volume of the heavy oil fraction (ft3/lbmol).
Riazi and Daubert  developed a general correlation given by the following expression:
In Eq. (24), θ1 and θ2 are two parameters accounting for the molecular forces and molecular size of a component. It was found that (Tb, γ) and (M, γ) are appropriate sets of input parameters for the correlation. Based on these results, two expressions were proposed:
The coefficients were adjusted to experimental data of 38 pure hydrocarbons with carbon number in the range of 1–20, including paraffins, olefins, naphthenes and aromatics with molecular weight between 70 and 300, and boiling point in the range of 80–650°F.
Magoulas and Tassios  correlated the critical properties of heavy fractions through the following expressions:
Riazi and Alsahhaf  developed a correlation for the acentric factor:
Sancet  proposed correlations for Pc and Tc as a function of the molecular weight of the heavy oil fractions:
From the aforementioned correlations, 22 sets were created and used to calculate the critical properties and acentric factor of the heavy fraction of two reservoir fluids (Table 7). These sets were evaluated by comparing the experimental data and simulated PVT analysis of these fluids. The PVT experiment was simulated using the Peng-Robinson equation of state :
where R is the universal gas constant, T is the temperature and v is the molar volume.
|Set||Pc and Tc||ω||Set||Pc and Tc||ω|
|1||Cavett ||Riazi and Alsahhaf ||12||Twu ||Riazi and Alsahhaf |
|2||Cavett ||Edmister ||13||Twu ||Edmister |
|3||Standing ||Riazi and Alsahhaf ||14||Riazi and Daubert ||Riazi and Alsahhaf |
|4||Standing ||Edmister ||15||Riazi and Daubert ||Edmister |
|5||Kesler and Lee ||Kesler and Lee ||16||Riazi and Daubert ||Riazi and Alsahhaf |
|6||Kesler and Lee ||Riazi and Alsahhaf ||17||Riazi and Daubert ||Edmister |
|7||Kesler and Lee ||Edmister ||18||Magoulas and Tassious ||Magoulas and Tassious |
|8||Riazi and Daubert ||Riazi and Alsahhaf ||19||Magoulas and Tassious ||Riazi and Alsahhaf |
|9||Riazi and Daubert ||Edmister ||20||Magoulas and Tassious ||Edmister |
|10||Sim and Daubert ||Riazi and Alsahhaf ||21||Sancet ||Riazi and Alsahhaf |
|11||Sim and Daubert ||Edmister ||22||Sancet ||Edmister |
In order to choose the most appropriate correlation for the calculation of the critical properties of heavy fractions of a reservoir fluid, 22 sets of correlations were used in the estimative of properties of the “cut” of fluids A and B . The compositions of the fluids A and B are presented in Tables 8 and 9, respectively.
|Component||Mole fraction (%)||Component||Mole fraction (%)||Component||Mole fraction (%)|
|Component||Mole fraction (%)||Component||Mole fraction (%)||Component||Mole fraction (%)|
Fluid A contains approximately 40% of methane and 6% of hydrocarbons with 20 or more carbons, while fluid B contains 13% of heavy fractions and 27% of methane.
|Set||Average Bo (%)||Average Rs (%)||Average ρo (%)||Average Bg (%)||Average Zg (%)||Average γg (%)||Average Pb (%)||Total deviation (%)|
|Set||Average Bo (%)||Average Rs (%)||Average ρo (%)||Average Zg (%)||Average Pb (%)||Total deviation (%)|
The best correlation for fluids A and B is chosen based on the sum of the average deviation for each property (Tables 10 and 11). For fluid A, the correlation of Twu  led to the best results in terms of the critical properties while for the fluid B it was the correlation of Kesler and Lee . For both fluids, the correlation of Riazi and Alsahhaf  was chosen to calculate the acentric factor (Table 12).
Figures 1–6 compare the experimental and calculated PVT data of fluid A using the Peng-Robinson equation of state with the critical pressure and temperature calculated from the correlation of Twu  and the acentric factor calculated from the Riazi and Alsahhaf  correlation. It is possible to note the excellent agreement between the experimental and simulated data. Figures 7–10 show the same results for fluid B, with similar performance.
Overall, 22 sets of correlations to calculate the critical properties and acentric factors of heavy oil fractions were tested in order to simulate a PVT experiment of two reservoir fluids. The results obtained from the simulation were compared with experimental laboratory PVT data. The Twu  correlation showed the best performance for one case, while the Kesler and Lee  correlation led to better results for the other fluid. In both cases, the Riazi and Alsahhaf  correlation was used to calculate the acentric factor. These results show the importance of evaluating the correlations for each reservoir fluid since no model for the critical properties and acentric factor calculations can be applied in any case.