Open access peer-reviewed chapter

PVT Properties of Black Crude Oil

Written By

Abdelaziz El-Hoshoudy and Saad Desouky

Submitted: 28 June 2018 Reviewed: 26 October 2018 Published: 18 December 2019

DOI: 10.5772/intechopen.82278

From the Edited Volume

Processing of Heavy Crude Oils - Challenges and Opportunities

Edited by Ramasamy Marappa Gounder

Chapter metrics overview

2,274 Chapter Downloads

View Full Metrics

Abstract

Precise PVT studies and behavior of phase-equilibrium of petroleum reservoir fluids are essential for describing these fluids and appraising their volumetric behavior at several pressure stages. There are numerous laboratory studies that can be performed on a reservoir sample. The amount of data desired determines the number of tests to be performed in the laboratory. Generally, there are three laboratory tests which characterize hydrocarbon fluids, namely primary study, constant mass depletion, and differential vaporization test. Generally, PVT properties are determined either experimentally or calculated theoretically through published correlations. This chapter presents different PVT laboratory tests that are required to understand the phase behavior of black oils.

Keywords

  • PVT experiments
  • black oil properties
  • petroleum reservoirs

1. Introduction

Petroleum (an equivalent term is a crude oil) is a complex mixture consisting predominantly of hydrocarbons and containing sulfur, nitrogen, oxygen, and helium as minor constituents. The physical and chemical properties of crude oils vary considerably and depend on the concentration of the various types of hydrocarbons and minor constituents present. Crude oil reservoirs are classified according to initial reservoir pressure into the following categories: (1) under-saturated oil reservoir; in which initial reservoir pressure is greater than the bubble point pressure of reservoir fluid; (2) saturated oil reservoir; in which initial reservoir pressure is equal to the bubble point pressure of reservoir fluid; (3) gas-cap reservoir; in which initial reservoir pressure is below the bubble point pressure of reservoir fluid, so reservoir is termed as gas-cap or two-phase reservoir, since the gas or vapor phase is underlain by an oil phase. Petroleum hydrocarbons exist as gaseous or liquid phase depending on reservoir temperature. If it is higher than the critical temperature of the fluid, the reservoir fluid is gas. Otherwise, the reservoir fluid is oil [1]. Black oil reservoir is considered as one of the most precious reservoir fluids. It constitutes the majority of oil reservoirs and exists in every basin. An accurate description of physical properties of crude oils is of considerable importance in the fields of both applied and theoretical science and, especially, in the solution of petroleum reservoir engineering problems. Some of these physical properties are of primary interest in petroleum engineering studies and detected through PVT tests which aim to determine reservoir fluid behavior at simulated reservoir conditions. Data on these fluid properties is usually determined by laboratory experiments performed on samples of actual reservoir fluids. In the absence of experimentally measured properties of crude oils, it is necessary for the petroleum engineer to determine the properties from empirically derived correlations. However; correlations are approximations and may be useful only in regional geological provinces.

Crude oils cover a wide range of physical properties and chemical compositions and are classified into black oils and near-critical or volatile-oils [2] depending on their phase behavior. The phase diagrams of two reservoir types are illustrated in Figures 1 and 2, while their physical properties are indicated in Table 1 according to classifications established by Moses [1], McCain [3], Whitson and Brule [4]. From the phase diagram, following the pressure reduction path as indicated by the vertical line EF, the iso-lines are scattered uniformly around the reservoir temperature which can be located at any region between points C&F. It is apparent that the reservoir temperature (Tres) in black oil reservoirs is less than the critical temperature (Tc). Moreover, the liberated gas below the bubble point pressure is considered as lean gas for industrial applications [5]. While in volatile oils, Tres is closest to the critical temperature. Another important differentiation lies in the volume of evolved gas below bubble point which is much greater in volatile oil than that observed for black oils [6]. Another important factor comprising gas oil ratio (GOR), and mole % of C7+, where GOR is obtained from field data and the mole % of C7+ is obtained from composition analysis. PVT properties for black oils in almost all reservoirs are required for production and surface facilities calculations to manage, develop, and forecast oil field behavior [7], as well as maximizing economic profit.

Figure 1.

Typical phase diagram for black oil [8].

Figure 2.

Typical phase diagram for volatile oil [8].

Property Black oil [7, 9] Volatile oil [1]
Fluid color Deep black Light colored turn to slightly reddish
Mole % of heptane plus (C7+) >17.5% up to 26.5% 12.5–17.5%
(GOR) ∼1500–1750 & up to 2000 scf/STB 2000–3000 scf/STB
API <45° ∼40° or higher
Co 3–150 × 10−6 psi−1 (under-saturated to highly saturated GOR oils) ……
βo <2 rbbl/STB >2 rbbl/STB

Table 1.

Criteria of black and volatile oil reservoirs.

PVT analysis is utilized by reservoir engineers to distinguish physical properties of reservoir fluids and variations in the volume and phase state that occurs during oil production [10, 11]. The phase envelope is characterized by some physical criteria such as oil formation volume factor, saturation pressure, gas oil ratio, oil density, oil viscosity, and oil isothermal compressibility. This phase envelope may be drawn through experimental data points or via software packages like PVTp, PVTsim, Eclipse and so on, based on the differential equation of states. The importance of PVT properties for reservoir performance analysis makes a dire need to develop a convenient way to predict these properties mathematically either by empirical correlations or equations of state. These empirical relations comprise two types. The first one belongs to the black oil type which forecast PVT behavior from the available PVT data, including saturation pressure, reservoir temperature, oil API index, gas-oil ratio, and reservoir pressure. The second one belongs to compositional models which rely on the equation of states and its hybrids. In these models besides PVT parameters, other measurements including fluid composition, critical temperature, molar masses and components acentric factor are utilized [12, 13]. Although accuracy of the empirical PVT correlations is often limited due to variations and complexity of multicomponent reservoir fluid systems, these correlations still used on field scale [14, 15, 16, 17].

Advertisement

2. PVT physical parameters and black oil modeling

In a black oil reservoir, the oil and gas formation volume factors, gas densities, the solution gas-oil ratio, and the viscosities of oil and gas are measured at reservoir temperature as a function of pressure. Once these measurements are, they can be applied in empirical correlations to obtain the relative in situ amounts of oil and gas during the production life of the reservoir.

2.1 Crude oil API gravity

Crude oil density is defined as the mass of a unit volume of the crude at a specified pressure and temperature and expressed in g/cc. In typical oil reservoirs, the oil density decreases with pressure depletion until it reaches the lowest value at the bubble point. The higher oil density at higher pressures resorts to the high oil molecules compressing per the unit volume. By reducing pressure beneath the bubble-point pressure, the gas is ejected from the solution which forms a free gas. The released gas is saturated with some of oil intermediate components, while the heavier one remains in the reservoir. These intermediates release results in the density reversal trend versus pressure beneath the saturation pressure. It is calculated either by compositional models or by correlating parameters. The API gravity is the preferred gravity scale, and related to the specific gravity which is defined as the ratio of the oil density to water density at 60°F and atmospheric pressure by the following relation:

API = 141.5 γ o 131.5 E1

2.2 Solution gas specific gravity

It is a dimensionless property that represents the ratio between hydrocarbon gas density and the air density at standard conditions. It is ascribed practically by the weighted average of the specific gravities of the separated gas from each separator, and expressed by the following relation:

γ g = Σ i = 1 n R sep i γ sep i + R st γ st Σ i = 1 n R sep i + R st E2

2.3 Gas solubility

It is defined as the number of standard cubic feet of gas that dissolve in one stock-tank barrel of crude oil at certain pressure and temperature and defined in SCF/STB [18].

R s = V g P , T V o sc E3

At standard conditions, dissolved gas is completely released from the oil, and therefore the oil contains no gas. Gas solubility increase with pressure increasing until it reaches the maximum value at saturation pressure. Below the bubble point pressure, gas is evolved with reduction in pressure and the gas solubility also decreases. Generally, the lighter the oil the higher the gas solubility, therefore volatile oils often has a higher solution GOR than black oils.

2.4 Bubble point pressure

It is defined as the highest pressure at which a bubble of gas is first liberated from the oil, and is determined experimentally by conducting a constant mass depletion test or estimated from empirical correlations. By reservoir pressure depletion beneath the saturation pressure, a released gas is formed in the reservoir. Since gas mobility is greater than the oil mobility, the produced GOR increase with continual pressure depletion beneath the saturation pressure. Other PVT properties are also greatly affected by reservoir pressure depletion.

2.5 Isothermal oil compressibility

It is the change in the fluid volume with respect to the change in pressure at isothermal conditions, expressed in psi−1 and calculated above the bubble point pressure from Eqs. (4 and 5).

C o = 1 d v d v d p E4
C o = 1 B o B o P T B g R s P T E5

Oil compressibility is determined during constant mass depletion test. In that procedure, the oil in the PVT cell is allowed to expand by reducing the pressure, and both the increase in oil volume and the decrease in oil pressure are measured. Oil compressibility above the saturation pressure greatly affects the material balance calculations for estimating volumetric reserve in oil reservoirs, as well as pressure transient analysis for under-saturated oils. Moreover, it is used in the estimation of vertical lift performance.

2.6 Gas compressibility factor

It is a dimensionless quantity used for gas correction at higher pressures and temperature, calculated from the basic equation of state as

Z = PV nRT E6

2.7 Oil formation volume factor, Bo

It is defined as the ratio of oil volume either at the reservoir or separator condition to the volume of oil at standard conditions.

B o res = V o p T res V o sc E7
B o sep = V o p T sep V o sc E8

The oil formation volume factor is expressed in units of reservoir volume over standard volume (bbl/STB). At the surface, the oil loses the dissolved gas it contained under high pressure and temperature in the reservoir. By pressure decline from reservoir to surface conditions, the following scenarios are generated: (1) The oil weight decrease owing to the loss of dissolved gas by pressure decrease; (2) The oil volume reduced slightly owing to the temperature reduction, however, the oil slightly expands owing to the pressure decrease. The volume reduction and expansion due to temperature and pressure successively usually cancel each other and no significant changes occur. Practically, (Bo) is determined from residual oil volume at 60°F at the end of differential liberation. The increase in Bo above bubble point pressure represents oil expansion by pressure decline so the oil volume at P & T (the numerator in Eq. 8) increase and consequently Bo increase. Below bubble point, the oil shrinkage due to gas evolving tends to overcome its expansion by pressure reduction, consequently Bo decline steeply.

2.8 Gas formation volume factor, Bg

Defined as the ratio of the gas volume at the temperature and pressure of any stage below the bubble point and the volume of the same gas at standard conditions through differential liberation (Eq. 9). Bg can also be calculated from composition analysis as a function of Z-factor as stated in Eq. 10.

B g = V g P , T V g sc E9
B g = 0.028269 Z T R P psia E10

2.9 Crude oil viscosity

It is defined as the resistance of the fluid to flow, expressed in centipoise (cp). It is affected by the oil composition, GOR, and reservoir temperature. Oil viscosity represents internal dynamic fluid friction relevant to fluid layers. Higher and lower viscosity depends on the gas/oil ratio. The behavior of viscosity is similar to that of density since oil viscosity shrinks with pressure decrease, owing to reduced friction between fluid layers, and reaches its minimal at the saturation pressure. With further pressure depletion more and more gas escaped from the oil, so the oil becomes denser and more viscous owing to loss of lighter components. Practically oil viscosity is measured by falling ball viscometer.

Table 2 summarizes some of the input physical parameters used to characterize black oil reservoirs.

# Property Input parameters
1 Bubble point pressure (Pb) Reservoir temperature (Tres); API; gas solubility (Rs); the specific gravity of the solution gas (γg); oil specific gravity (γo); separator gas gravity (γg sep); reservoir pressure (Pres); separator pressure (Psep); separator temperature (Tsep)
2 Crude oil gravity (API) The specific gravity of the solution gas (γg); oil specific gravity (γo)
3 The specific gravity of the solution gas (γg) Separator GOR (Rsp); separator gas gravity (γsep); GOR from the stock tank (Rst); gas gravity from the stock tank (γst)
4 Oil density (ρo) The molecular weight of the stock tank oil (MSTO); API; gas solubility (Rs); the specific gravity of the solution gas (γg); reservoir temperature (Tres); reservoir pressure (Pres)
5 Gas solubility (Rs) Reservoir temperature (Tres); reservoir pressure (Pres); oil density (ρo); API; the specific gravity of the solution gas (γg); oil specific gravity (γo); separator pressure (Psep); separator temperature (Tsep); oil formation volume factor (Bo)
6 Oil formation volume factor (Bo) Reservoir temperature (Tres); reservoir pressure (Pres); oil density (ρo); API; specific gravity of the solution gas (γg); oil specific gravity (γo); separator pressure (Psep); separator temperature (Tsep); isothermal compressibility (Co); bubble point pressure (Pb); gas solubility (Rs)
7 Isothermal compressibility (Co) Above bubble point (Pb) Reservoir pressure (Pres > Pb); oil density (ρo); API; specific gravity of the solution gas (γg); separator pressure (Psep); separator temperature (Tsep); bubble point pressure (Pb); gas solubility at the bubble point pressure (Rs)
Below bubble point (Pb) Reservoir temperature (Tres); reservoir pressure (Pres < Pb); API; specific gravity of the solution gas (γg); oil specific gravity (γo); oil formation volume factor (Bo); gas formation volume factor (Bg); bubble point pressure (Pb); gas solubility (Rs)
8 Oil viscosity (μo) Dead oil viscosity (μod) Reservoir temperature (Tres); API
Saturated oil viscosity (μob) Reservoir temperature (Tres); reservoir pressure (Pres); API; specific gravity of the solution gas (γg); oil specific gravity (γo); bubble point pressure (Pb); gas solubility (Rs); oil density at bubble point (ρob); dead oil viscosity (μod)
Under-saturated oil viscosity (μou) Reservoir pressure (Pres); API; bubble point pressure (Pb); gas solubility (Rs); dead oil viscosity (μod); saturated oil viscosity (μob)

Table 2.

Major input parameters for PVT equations.

Advertisement

3. Experimental PVT analysis

Crude oil samples can be obtained from downhole using a wireline sampling chamber, or at the separator, using separation facilities [19]. To carry a successful PVT analysis, the samples should be generated through the following procedure:

3.1 Sample handling

The pressurized liquid sample is heated to 200°F and stabilized in single phase prior to removal of any sample in order to ensure sample integrity and avoid any wax precipitation.

3.2 Sample validation

Carried out for both surface and bottom hole samples to ensure that no leakage occurs for the samples and the fluid virtually represents the reservoir.

3.2.1 Surface samples obtained from separator

They are subjected to assuring cylinders integrity during transportation through measurement of opening pressures of the separator oil cylinder and separator gas cylinder at ambient temperature and separator temperature. If the measured pressure equals the separator pressure, the sample is valid, otherwise, the sample is invalid. Moreover, quality control of separator samples carried out through measurement of the saturation pressure (bubble point) of the separator oil sample at the separator temperature and compared with the separator pressure. The sample was then transferred to the Lab cylinder, and the amount of water if present in the separator oil was determined, and discharged.

3.2.2 Bottom hole sample

Also, check cylinders’ integrity during transportation and sampling carried out by measurement of sample opening pressures at transferring temperature and reservoir temperature. The sample cylinder was then agitated and transferring to Lab cylinder at reservoir temperature and pressure of 1000 psi above reservoir pressure.

3.3 Primary study

3.3.1 Surface samples obtained from separator

A primary study carried only for surface samples is as follows:

3.3.1.1 Atmospheric flashing of separator oil

A portion of the separator oil sample was charged to a high-pressure visual PVT cell and thermally expanded to the separator temperature till equilibrium. The sample was subjected to flash separation from separator conditions to standard conditions of 14.73 psia and 60°F. Dissolved gas-oil ratio and oil formation volume factor were calculated using the volumetric data obtained so that the supplied stock-tank oil rate could be adjusted to separator oil rate for the well-stream calculations. The density and API gravity of the stock-tank oil were also measured using density meter. The dissolved gas composition was determined using extended gas chromatography and the flashed liquid by high-resolution capillary gas chromatography. The two analyses were then mathematically recombined according to the dissolved gas-oil ratio to get the composition of separator oil. Figure 3 illustrates the primary study procedure.

Figure 3.

Diagram of the primary study steps.

3.3.1.2 Composition analysis of stock-tank oil

The compositions of paraffins and aromatics in liquid samples were determined using Clarus 500 Perkin Elmer Gas Chromatograph connected with flame ionization detector (FID) using selective PIONA capillary column of 100 m in length and 0.25 mm internal diameter. Suitable sample capacity was injected into split / splitless injector through a microsyringe according to both response and linear range of FID detector. Nitrogen gas was used as carrier gas at its optimum flow rate. The system enables to detect the composition up to C36+. Further oil analysis includes density, viscosity, pour point, and distillation analysis, in addition to water, basic sediment (BS) and sulfur contents, in order to determine the oil quality.

3.3.1.3 Gas composition analysis

The analysis of both separator gas and dissolved gas were determined according to ASTM 1945, using Clarus 500 PerkinElmer Gas Chromatograph online with the condensate cell utilizing four packed columns connected in series with thermal conductivity detector and DP-1 capillary column matched with flame ionization detector. Helium was used as carrier gas and the oven isothermal temperature 100°C in order to provide highly efficient detection of aromatics and hydrocarbons paraffin (normal and iso) from methane (C1) up to undecane (C11) components and non-hydrocarbons including H2S, CO2, and N2. Once repeatability is verified, the validated bottle generated to proceed the complete PVT analysis, which, in the case of the reservoir fluid be a black oil, is composed by three main assays; constant mass depletion (CMD), differential liberation (DL), and separator test.

3.3.1.4 Measurement of density

The density of oil is measured by Anton Paar Density Meter. The density meter is calibrated with air and distilled water. The procedure reference is ASTM D-4052.

3.3.1.5 Measurement of molecular weight

The average molecular weight number is determined by using cryoscopy, (freezing point depression of benzene).

3.3.1.6 Corrected and adjustment of field separator GOR

Adjustment of field GOR carried out in the lab through the following equations, owing to higher deviations of the measured field gas deviation factor (Zf)

corrected separator GOR = field GOR × ρ f × Z f ρ l × Z l 0.5 E11
Zf = 1 fpv 0.5 E12
adjusted separator GOR = corrected separator GOR B o E13

3.3.1.7 Quality control of primary study

Carried out by applying Hoffmann method to data of separator oil & gas through plotting of (Ki Psep) against [1/Bi(1/Tbi − 1/Tsep)] on a Cartesian scale which results in a straight line, as shown in Figure 4.

Figure 4.

Hofmann QC-plot.

where,

Ki = yi xi E14
B i = log P ci / 14.7 / 1 / T bi 1 / T ci E15

3.3.1.8 Physical recombination of well stream

The samples were physically recombined at separator conditions (Psep & Tsep) using the corrected separator (GOR), (Bo) and (Bg) to form the wellstream. The wellstream composition can mathematically be determined by applying material balance equation to the separator product compositional analysis and their physical properties such as corrected separator GOR, dissolved GOR, Bo, the density of stock tank oil, separator oil and gas densities as well as Bg.

3.3.2 Bottom hole sample (BHS)

A known portion of the reservoir fluid was charged from (BHS) cylinder to a high-pressure visual cell and thermally expanded to the reported reservoir temperature and pressure. After establishing thermal equilibrium, the fluid sample was subjected to a flash separation at standard conditions (P = 14.7 psia & T = 60°F). The flashed data were then checked using Hoffmann method.

3.4 Constant mass depletion (CMD)

An integrated and automated mercury free system (PVT-cell) is used for PVT analysis of black-oil and volatile oil. A portion of the sample is pumped into the cell, and the volume and, therefore, pressure of the sample chamber are controlled with a piston. Most sample chambers are equipped with a magnetic stirrer to ensure that the system is well mixed. Constant mass depletion also known as flash liberation, in which a set of volume measurements as a function of pressure at fixed reservoir temperature are performed within the PVT cell. The recombined well stream was pressurized and heated to reservoir temperature and pressure. A portion of the sample (∼100 cc of live oil) was charged to a high-pressure visual PVT cell and thermally expanded to the reservoir temperature and the pressure was adjusted to a pressure of 1000 psi higher than the reservoir pressure. The measurements of pressures and volumes were recorded and then used to determine the bubble point pressure of the under-saturated reservoirs at reservoir temperature. Above bubble point pressure, the relative volume, oil formation volume factor, oil compressibility, density, and viscosity were determined, while below bubble point pressure the Y-Function is calculated at each pressure step.

3.5 Differential liberation (DL)

After the constant composition mass depletion, the sample was pressurized and equilibrated to a single phase condition. A differential vaporization procedure is then performed for the purpose of measuring two-phase properties as a function of reservoir pressure depletion. As with the flash expansion test, the volume of the sample chamber is increased incrementally; however, in this case, the evolved solution gas is withdrawn from the sample chamber after each incremental expansion. At each pressure depletion, the oil formation volume factor, gas formation volume factor, two-phase volume factor, gas gravity, gas solubility, density, and viscosity were determined. The viscosity of the evolved gas is usually calculated through Lee et al. correlations rather than measured. While the viscosity measurement for both live and dead oil is carried out over a wide range of high pressure and temperature by a rolling ball viscometer (RUSKA). Schematic representation of CMD and differential liberation are shown in Figure 5.

Figure 5.

Schematic illustration of (a) constant mass depletion and (b) differential liberation [6].

3.6 Separator test

Separator test is used to determine the effect of surface separation parameters on the gas-oil ratio, stock tank oil gravity and formation volume factor at separator conditions. About 100 cc of the recombined sample is raised to the saturation pressure at the reservoir temperature. Then the sample condition is lowered and adjusted to the separator conditions, where the oil and gas volumes are measured. In presence of the second separator, the sample is adjusted to the second separator temperature and pressure, where the oil and gas volumes are also measured. After testing all separator stages, the sample is flashed to stock tank conditions, where the volumes of stock tank oil and gas are monitored. The separator test volumetric data and the compositional analysis of the evolved gases from each pressure stage are reported. The separator test data were then used to adjust differential vaporization result to the surface conditions used in the separator test. Schematic representation of the multistage separator test is represented in Figure 6.

Figure 6.

Multistage separator test.

Advertisement

4. PVT correlations and software

Numerous correlations to calculate PVT properties for oils have been developed. The correlations are usually simple expressions (equations) that allow engineers to calculate the PVT properties from readily available inputs in case of absence of experimentally measured data. Correlations are generated from PVT Datasets associated with the particular geographical region. The most commonly known geographical correlations are Al-Marhoun [20] correlations for Saudi Arabian crudes; Standing [21] for California crudes; Petrosky [22] for Gulf of Mexico crudes; Hanafy et al. [23] for Egyptian crudes; and Omar and Todd [24] for Malaysian crudes. Some developed correlations use common oil samples of worldwide origin. The correlation is usually developed by fitting a predetermined form of an equation (or more) with the measured data and calculating the equation(s) parameters by use of multiple regression analysis. It is essential to note that the accuracy of correlation performance and applicability is associated with the range of data used in correlation development. PVT correlations developed as a function of PVT properties, or as groups of correlations developed from the same dataset for several PVT properties. Table 3 gives a detailed survey on the published PVT correlations. Numerous software were developed for prediction of PVT phase behavior and physical properties such as PVTp, PVTsim, Eclipse E100, Eclipse E300, Petrel, and Ecrin. A detailed description of this software will be described elsewhere.

Correlation name Pb Rs Cob Bob ρob μob Cou (GOR)STO Bou μou μod
Abdul-Majeed and Salman [25]
Abdul-Majeed et al. [26]
Abu-Khamsim and Al-Marhoun [27]
Ahmed [28]
Ahmed [29]
Al-Khafaji et al. [30]
Al-Marhoun [31]
Al-Marhoun [20]
Al-Marhoun [32]
Al-Marhoun [33]
Al-Marhoun [34]
Al-Mehaideb [35]
Al-Najjar et al. [36]
Al-Shammasi [37]
Andrade [38]
Frick [39]
Asgarpour et al. [40]
Asgarpour et al. [40]
Beal [41]
Beggs and Robinson [42]
Bennison [43]
Bergman [44]
Bergman and Sutton [45]
Sutton and Bergman [46]
Bolondarzadeh et al. [47]
Calhoun [48]
Casey and Cronquist [49]
Chew and Connally [50]
De Ghetto et al. [51]
Dindoruk and Christman [52]
Doklah and Osman [53]
Elam [54]
El-Banbi et al. [55]
Elmabrouk et al. [56]
Elsharkawy and Alikhan [57]
Elsharkawy and Alikhan [58]
Elsharkawy and Gharbi [59]
Farshad et al. [60]
Fitzgerlad [61]
Glaso [62]
Hanafy et al. [23]
Hassan [63]
Hemmati and Kharrat [64]
Hossain et al. [65]
Ikiensikimama et al. [66]
Ikiensikimama and Azubuike [67]
Ikiensikimama and Ogboja [68]
Karimnezhad et al. [69]
Kartoatmdjo and Schmidt [70]
Kartoatmdjo and Schmidt [71]
Kaye [72]
Khairy et al. [73]
Khamehchi et al. [74]
Khan et al. [75]
Khazam et al. [76]
Knopp and Ramsey [77]
Kouzel [78]
Labedi [79]
Labedi [80]
Labedi [81]
Lasater [82]
Levitan and Murtha [83]
Macary and El Batanony [84]
Mazandarani and Asghari [85]
McCain [14]
McCain et al. [86]
Mehran and Movagharnejad [87]
Moradi et al. [88]
Moradi et al. [89]
Mehran and Movagharnejad [87]
Naseri et al. [90]
Naseri et al. [91]
Ng and Egbogah [92]
Obomanu and Okpobiri [93]
Okeke and Sylvester [94]
Okoduwa and Ikiensikimama [95]
Omar and Todd [24]
Osorio [96]
Ostermann and Owolabi [97]
Owolabi [98]
Oyedeko and Ulaeto [99]
Petrosky [22]
Petrosky and Farshad [100]
Petrosky and Farshad [101]
Petrosky and Farshad [102]
Rollins et al. [103]
Sattarina et al. [104]
Spivey et al. [105]
Standing [21]
Standing [106]
Standing [107]
Sim [108]
Sulaimon et al. [109]
Twu [110]
Oyedeko and Ulaeto [99]
Valko and McCain [111]
Vasquez and Beggs [112]
Velarde et al. [113]
Whitson and Brule [4]

Table 3.

Summary of the published PVT correlations.

Advertisement

5. Conclusion

PVT studies are very critical and precious for reservoir engineers to monitor the phase behavior of petroleum fluids and generating material balance calculations. Traditionally, PVT experiments are carried out through high pressurized PVT-cell, but in case of unavailability of experimental facilities, we resort to empirically derived correlations. Several empirical correlations were published depending on geographical provinces. In this chapter, a brief description of the phase diagram for black and volatile oil was introduced as well as the importance of PVT properties in predicting petroleum fluids behavior. Then PVT physical parameters and its meanings were briefly illustrated and a summary of the input parameters in their predictions as well as a survey of the published correlations were provided. Moreover, a detailed description of the experimental PVT procedures and tests were illustrated to help engineers and analysts in conducting such tests.

Advertisement

Nomenclature

(GOR)STO

stock-tank gas/oil ratio

μob

saturated oil viscosity (at or below bubble point)

μod

dead oil viscosity correlations

μou

under-saturated oil viscosity

bbl/STB

barrel/stock tank barrel

Bg

gas formation volume factor

Bo

oil formation volume factor, bbL/STB

Bob

oil formation volume factor at the bubble point

Bou

under-saturated oil formation volume factor

Co

oil compressibility

Cob

saturated oil compressibility (at or below bubble point)

Cou

under-saturated oil compressibility (above bubble point)

GOR

producing gas-oil ratio

Ki

the equilibrium constant of component i in the recombined sample

N

number of separators

Pb

bubble point pressure

Pci

critical pressure of component i in the recombined sample

PVT

pressure-volume-temperature

Rs

gas solubility

Rsp

separator GOR, scf/STB

Rst

GOR from the stock-tank, SCF/STB

Scf/STB

standard cubic foot/stock tank barrel

Tbi

true boiling point temperature of component i in the recombined sample

Tci

the critical temperature of component i in the recombined sample

xi

the mole fraction of component i in the separator oil

yi

the mole fraction of component i in the separator gas

Z

gas compressibility factor

Zf

gas compressibility factor measured at the field

Zl

gas compressibility factor measured at PVT lab

γg

solution gas specific gravity

γsep

separator gas gravity

γst

gas gravity from the stock-tank

ρf

gas gravity measured at the field

ρl

gas gravity measured at PVT lab

ρob

oil density at the bubble point

References

  1. 1. Moses PL. Engineering applications of phase behavior of crude oil and condensate systems (includes associated papers 16046, 16177, 16390, 16440, 19214 and 19893). Journal of Petroleum Technology. 1986;38(07):715-723
  2. 2. Mukerji T, Jørstad A, Avseth P, Mavko G, Granli J. Mapping lithofacies and pore-fluid probabilities in a North Sea reservoir: Seismic inversions and statistical rock physics. Geophysics. 2001;66(4):988-1001
  3. 3. McCain WD Jr. Revised gas-oil ratio criteria key indicators of reservoir fluid type; part 5. Petroleum Engineer International; (United States). 1994;66:4
  4. 4. Whitson C, Brule M. Phase behavior. In: Monograph, Vol. 20 SPE Henry L. Doherty Series. Richardson, Texas: SPE Monograph series; 2000
  5. 5. Holditch S, Hill AD, Zhu D. Advanced Hydraulic Fracturing Technology for Unconventional Tight Gas Reservoirs. US department of energy; Texas Engineering Experiment Station; 2007
  6. 6. Danesh A. PVT and Phase Behaviour of Petroleum Reservoir Fluids. Scothland: Elsevier; 1998
  7. 7. Jarrahian A, Moghadasi J, Heidaryan E. Empirical estimating of black oils bubblepoint (saturation) pressure. Journal of Petroleum Science and Engineering. 2015;126:69-77
  8. 8. Ahmed T. Reservoir Engineering Handbook. Amestrdam: Elsevier; 2006
  9. 9. McCain WD, Spivey JP, Lenn CP. Petroleum Reservoir Fluid Property Correlations. Houston, Texas: PennWell Books; 2011
  10. 10. Di Primio R, Dieckmann V, Mills N. PVT and phase behaviour analysis in petroleum exploration. Organic Geochemistry. 1998;29(1–3):207-222
  11. 11. Mansour E, Farag A, El-Dars F, Desouky S, Batanoni M, Mahmoud M. Predicting PVT properties of Egyptian crude oils by a modified soave–Redlich–Kowng equation of state. Egyptian Journal of Petroleum. 2013;22(1):137-148
  12. 12. El-Hoshoudy A, Farag A, Ali O, El-Batanoney M, Desouky S, Ramzi M. New correlations for prediction of viscosity and density of Egyptian oil reservoirs. Fuel. 2013;112:277-282
  13. 13. El-Hoshoudy A, Desouky S, Betiha M, Abdelraouf S, Elaily M. Mathematical correlation predicting stock tank oil density for Egyptian crude oils at standard conditions. International Journal of Oil, Gas and Coal Engineering. 2014;2(4):55-58
  14. 14. McCain JW. Reservoir-fluid property correlations-state of the art (includes associated papers 23583 and 23594). SPE Reservoir Engineering. 1991;6(02):266-272
  15. 15. El-Hoshoudy AN, Desouky SM. An empirical correlation for estimation of formation volume factor of gas condensate reservoirs at separator conditions. Petroleum & Petrochemical Engineering Journal. 2018;2(2):9
  16. 16. El-Hoshoudy AN. Numerical prediction of oil formation volume factor at bubble point for black and volatile oil reservoirs using nonLinear regression models. Petroleum & Petrochemical Engineering Journal. 2018;2(2):11
  17. 17. El-Hoshoudy AN, Desouky SM, Gomaa S. Prediction of dew point pressure in gas condensate reservoirs based on a combination of gene expression programming (GEP) and multiple regression analysis. Petroleum & Petrochemical Engineering Journal. 2018;2(5)
  18. 18. Vazquez M, Beggs HD. Correlations for fluid physical property prediction. In: SPE Annual Fall Technical Conference and Exhibition. Society of Petroleum Engineers; 1977
  19. 19. Starling KE. Fluid Thermodynamic Properties for Light Petroleum Systems. Houston, Texas: Gulf Publishing Company; 1973
  20. 20. Al-Marhoun MA. PVT correlations for Middle East crude oils. Journal of Petroleum Technology. 1988;40(05):650-666
  21. 21. Standing M. A pressure-volume-temperature correlation for mixtures of California oils and gases. In: Drilling and Production Practice. Washington, United states: American Petroleum Institute; 1947
  22. 22. Petrosky GE. PVT Correlations for Gulf of Mexico Crude Oils. Lafayette, USA: University of Southwestern Louisiana; 1990
  23. 23. Hanafy H, Macary S, ElNady Y, Bayomi A, El Batanony M. Empirical PVT correlations applied to Egyptian crude oils exemplify significance of using regional correlations. In: International Symposium on Oilfield Chemistry; Society of Petroleum Engineers. 1997
  24. 24. Omar M, Todd A. Development of new modified black oil correlations for Malaysian crudes. In: SPE Asia Pacific Oil and Gas Conference; Society of Petroleum Engineers. 1993
  25. 25. Abdul-Majeed GH, Salman NH. An empirical correlation for oil FVF prediction. Journal of Canadian Petroleum Technology. 1988;27(06):118-122
  26. 26. Abdul-Majeed GH, Kattan RR, Salman NH. New correlation for estimating the viscosity of undersaturated crude oils. Journal of Canadian Petroleum Technology. 1990;29(03):80-85
  27. 27. Abu-Khamsin SA, Al-Marhoun MA. Development of a new correlation for bubble-point oil viscosity. Arabian Journal of Science and Engineering. 1991;16:99-106
  28. 28. Ahmed T. Hydrocarbon Phase Behavior (Contributions in Petroleum Geology and Engineering). Houston, Texas: Gulf Publishing Company; 1989
  29. 29. Ahmed T. Oil Recovery Mechanisms and the Material Balance Equation Reservoir Engineering Handbook. Houston, Texas: Gulf Professional Publishing; 2010. pp. 733-809
  30. 30. Al-Khafaji AH, Abdul-Majeed GH, Hassoon SF. Viscosity correlation for dead, live and undersaturated crude oils. Journal of Petroleum Research. 1987;6(2):1-16
  31. 31. Al-Marhoun M. Pressure-Volume-Temperature Correlations for Saudi Crude Oils. Saudi Arabia: University of Petroleum and Minerals; 1985
  32. 32. Al-Marhoun MA. New correlation for formation volume factor of oil and gas mixtures. Journal of Canadian Petroleum Technology. 1992;31(3):22-26
  33. 33. Al-Marhoun MA. The coefficient of isothermal compressibility of black oils. In: Middle East Oil Show. Bahrain: Society of Petroleum Engineers; 2003
  34. 34. Al-Marhoun MA. A New Correlation for Undersaturated Isothermal Oil Compressibility. 2006
  35. 35. Al-Mehaideb R. Improved PVT correlations for UAE offshore crudes. Journal of the Japan Petroleum Institute. 1997;40(3):232-235
  36. 36. Al-Najjar H, Al-Soof N, Al-Khalisy K. Correlations for bubble-point pressures, gas oil ratios and formation volume factors for Iraqi crude oils. Journal of Petroleum Research. June 1988;1988:13
  37. 37. Al-Shammasi A. Bubble point pressure and oil formation volume factor correlations. In: Middle East Oil Show and Conference; Society of Petroleum Engineers. 1999
  38. 38. Andrade EC. The Viscosity of Liquids. London, United Kingdom: Nature Publishing Group. 1930
  39. 39. Frick TC. Petroleum Production Handbook. Newyork, United states: McGraw-Hill; 1962
  40. 40. Asgarpour S, McLauchlin L, Wong D, Cheung V. Pressure-volume-temperature correlations for western Canadian gases and oils. Journal of Canadian Petroleum Technology. 1989;28(04):103-112
  41. 41. Beal C. The viscosity of air, water, natural gas, crude oil and its associated gases at oil field temperatures and pressures. Transactions of the AIME. 1946;165(01):94-115
  42. 42. Beggs HD, Robinson J. Estimating the viscosity of crude oil systems. Journal of Petroleum Technology. 1975;27(09):1140-1141
  43. 43. Bennison T. Prediction of heavy oil viscosity. Presented at the IBC Heavy Oil Field Development Conference. 2-4 December; 1998
  44. 44. Bergman D. Don’t forget viscosity. Petroleum Technology Transfer Council 2nd Annual Reservoir Engineering Symposium; Lafayette, Louisiana. 28 July; 2004
  45. 45. Bergman DF, Sutton RP. An update to viscosity correlations for gas-saturated crude oils. In: SPE Annual Technical Conference and Exhibition; Society of Petroleum Engineers. 2007
  46. 46. Sutton RP, Bergman DF. Undersaturated oil viscosity correlation for adverse conditions. In: SPE Annual Technical Conference and Exhibition; Society of Petroleum Engineers. 2006
  47. 47. Bolondarzadeh A, Hashemi S, Solgani B. The new PVT generated correlations of Iranian oil properties. In: 4th Iranian Petroleum Engineering Student Conference. 2006
  48. 48. Calhoun J Jr. New methods for determining permeability. The Petroleum Engineer. February 1947:103-108
  49. 49. Casey J. Estimate GOR and FVF using dimensionless PVT analysis. World Oil (United States). 1992;213(11):12-16
  50. 50. Chew J-N, Connally CA Jr. A Viscosity Correlation for Gas-Saturated Crude Oils. 1959
  51. 51. De Ghetto G, Paone F, Villa M. Pressure-volume-temperature correlations for heavy and extra heavy oils. In: SPE International Heavy Oil Symposium; Society of Petroleum Engineers. 1995
  52. 52. Dindoruk B, Christman PG. PVT properties and viscosity correlations for Gulf of Mexico oils. SPE Annual Technical Conference and Exhibition; Society of Petroleum Engineers; 2001
  53. 53. Dokla M, Osman M. Correlation of PVT properties for UAE crudes (includes associated papers 26135 and 26316). SPE Formation Evaluation. 1992;7(01):41-46
  54. 54. Elam FM. Prediction of Bubble Point Pressures and Formation Volume Factors from Field Data. University of Texas; 1957:53-79
  55. 55. El-Banbi AH, Fattah KA, Sayyouh H. New modified black-oil PVT correlations for gas condensate and volatile oil fluids. In: SPE Annual Technical Conference and Exhibition; Society of Petroleum Engineers. 2006
  56. 56. Elmabrouk SK, Zekri AY, Shirif EE. Prediction of bubblepoint pressure and bubblepoint oil formation volume factor in the absence of PVT analysis. In: SPE Latin American and Caribbean Petroleum Engineering Conference; Society of Petroleum Engineers. 2010
  57. 57. Elsharkawy AM, Alikhan AA. Correlations for predicting solution gas/oil ratio, oil formation volume factor, and undersaturated oil compressibility. Journal of Petroleum Science and Engineering. 1997;17(3–4):291-302
  58. 58. Elsharkawy A, Alikhan A. Models for predicting the viscosity of Middle East crude oils. Fuel. 1999;78(8):891-903
  59. 59. Elsharkwy A, Gharbi R. Comparing classical and neural regression techniques in modeling crude oil viscosity. Advances in Engineering Software. 2001;32(3):215-224
  60. 60. Frashad F, LeBlanc J, Garber J, Osorio J. Empirical PVT correlations for Colombian crude oils. In: SPE Latin America/Caribbean Petroleum Engineering Conference; Society of Petroleum Engineers. 1996
  61. 61. Fitzgerald D. A Predictive Method for Estimating the Viscosity of Undefined Hydrocarbon Liquid Mixtures. PA, United States: Pennsylvania State University, State College; 1994
  62. 62. Glaso O. Generalized pressure-volume-temperature correlations. Journal of Petroleum Technology. 1980;32(05):785-795
  63. 63. Hassan OF. Correlation for solution gas-oil ratio of Iraqi oils at pressures below the bubble point pressure. Iraqi Journal of Chemical and Petroleum Engineering. 2011;12(2):1-8
  64. 64. Hemmati MN, Kharrat R. A correlation approach for prediction of crude oil PVT properties. In: SPE Middle East Oil and Gas Show and Conference; Society of Petroleum Engineers. 2007
  65. 65. Hossain MS, Sarica C, Zhang H-Q, Rhyne L, Greenhill K. Assessment and development of heavy oil viscosity correlations. In: SPE International Thermal Operations and Heavy Oil Symposium; Society of Petroleum Engineers. 2005
  66. 66. Ikiensikimama S, Egbe T, Effiong E, Ogboja O. New screening technique for PVT correlations as applied to the Niger Delta crude. In: Refereed Proceedings. 2008
  67. 67. Ikiensikimama S, Azubuike I. Modeling approach for Niger-delta oil formation volume factor prediction using artificial neural network. In: Nigeria Annual International Conference and Exhibition; Society of Petroleum Engineers. 2012
  68. 68. Ikiensikimama SS, Ogboja O. New bubblepoint pressure empirical PVT correlation. In: Nigeria Annual International Conference and Exhibition; Society of Petroleum Engineers. 2009
  69. 69. Karimnezhad M, Heidarian M, Kamari M, Jalalifar H. A new empirical correlation for estimating bubble point oil formation volume factor. Journal of Natural Gas Science and Engineering. 2014;18:329-335
  70. 70. Kartoatmodjo TR, Schmidt Z. New Correlations for Crude Oil Physical Properties. 1991
  71. 71. Kartoatmodjo T, Schmidt Z. Large data bank improves crude physical property correlations. Oil and Gas Journal;(United States). 1994;92(27)
  72. 72. Kaye S. Offshore California viscosity correlations. COFRC, TS85000940; 1985
  73. 73. Khairy M, El-Tayeb S, Hamdallah M. PVT correlations developed for Egyptian crudes. Oil and Gas Journal. 1998;96(18):114-116
  74. 74. Khamehchi E, Rashidi F, Rasouli H, Ebrahimian A. Novel empirical correlations for estimation of bubble point pressure, saturated viscosity and gas solubility of crude oils. Petroleum Science. 2009;6(1):86-90
  75. 75. Khan S, Al-Marhoun M, Duffuaa S, Abu-Khamsin S. Viscosity correlations for Saudi Arabian crude oils. In: Middle East Oil Show. Society of Petroleum Engineers; 1987
  76. 76. Khazam M, Alkhaboli M, Shlak M. New PVT Correlations for Libyan Crudes
  77. 77. Knopp C, Ramsey L. Correlation of oil formation volume factor and solution gas-oil ratio. Journal of Petroleum Technology. 1960;12(08):27-29
  78. 78. Kouzel B. How pressure affects liquid viscosity. Hydrocarbon Processing. March 1965;1965:120
  79. 79. Labedi RM. PVT Correlations of the African Crudes. Colorado School of Mines; 1982
  80. 80. Labedi RM. Use of production data to estimate the saturation pressure, solution GOR, and chemical composition of reservoir fluids. In: SPE Latin America Petroleum Engineering Conference; Society of Petroleum Engineers. 1990
  81. 81. Labedi R. Improved correlations for predicting the viscosity of light crudes. Journal of Petroleum Science and Engineering. 1992;8(3):221-234
  82. 82. Lasater J. Bubble point pressure correlation. Journal of Petroleum Technology. 1958;10(05):65-67
  83. 83. Levitan LL, Murtha M. New correlations estimate PH, FVF. Oil & Gas Journal. 1999;97(10):10
  84. 84. Macary S, El-Batanoney M. Derivation of PVT correlations for the Gulf of Suez crude oils. Journal of the Japan Petroleum Institute. 1993;36(6):472-478
  85. 85. Mazandarani MT, Asghari SM. Correlations for predicting solution gas-oil ratio, bubblepoint pressure and oil formation volume factor at bubblepoint of Iran crude oils. In: European Congress of Chemical Engineering (ECCE-6). 2007. pp. 16-20
  86. 86. McCain WD Jr, Rollins JB, Lanzi AJV. The coefficient of isothermal compressibility of black oils at pressures below the bubblepoint. SPE Formation Evaluation. 1988;3(03):659-662
  87. 87. Mehran F, Movagharnejad K, Didanloo A. New correlation for estimation of formation volume factor and bubblepoint pressure for Iranian oil fields. In: 1st Iranian Petroleum Engineering Conference. 2006
  88. 88. Moradi B, Malekzadeh E, Mohammad A, Awang M, Moradie P. New oil formation volume factor empirical correlation for Middle East crude oils. International Journal of Petroleum and Geoscience Engineering (IJPGE);2013:12-23
  89. 89. Moradi B, Malekzadeh E, Amani M, Boukadi FH, Kharrat R. Bubble point pressure empirical correlation. In: Trinidad and Tobago Energy Resources Conference; Society of Petroleum Engineers. 2010
  90. 90. Naseri A, Nikazar M, Dehghani SM. A correlation approach for prediction of crude oil viscosities. Journal of Petroleum Science and Engineering. 2005;47(3–4):163-174
  91. 91. Naseri A, Yousefi S, Sanaei A, Gharesheikhlou A. A neural network model and an updated correlation for estimation of dead crude oil viscosity. Brazilian Journal of Petroleum and Gas. 2012;6(1)
  92. 92. Ng J, Egbogah E. An improved temperature-viscosity correlation for crude oil systems. In: Annual Technical Meeting; Petroleum Society of Canada. 1983
  93. 93. Obomanu D, Okpobiri G. Correlating the PVT properties of Nigerian crudes. Journal of Energy Resources Technology. 1987;109(4):214-217
  94. 94. Okeke H, Sylvester O. Improved correlation for predicting stock tank gas-oil ratio in Niger Delta. In: SPE Nigeria Annual International Conference and Exhibition; Society of Petroleum Engineers. 2016
  95. 95. Okoduwa I, Ikiensikimama S. Bubble point pressure correlations for Niger Delta crude oils. In: Nigeria Annual International Conference and Exhibition; Society of Petroleum Engineers. 2010
  96. 96. Osorio JG. Empirical PVT Correlations for Colombian Crude Oil. University of Southwestern Louisiana; 1990
  97. 97. Ostermann R, Owolabi O. Correlations for the reservoir fluid properties of Alaskan crudes. In: SPE California Regional Meeting; Society of Petroleum Engineers. 1983
  98. 98. Owolabi OO. Reservoir Fluid Properties of Alaskan Crudes. Fairbanks: University of Alaska; 1984
  99. 99. Oyedeko KF, Ulaeto UW. Predicting the dead oil viscosity of reservoir fluids: A case study of the Niger Delta. Journal of Energy Technology and Policy. 2011;3:1-7
  100. 100. Petrosky JG, Farshad F. Pressure-volume-temperature correlations for Gulf of Mexico crude oils. In: SPE Annual Technical Conference and Exhibition; Society of Petroleum Engineers. 1993
  101. 101. Petrosky G Jr, Farshad F. Viscosity correlations for Gulf of Mexico crude oils. In: SPE Production Operations Symposium; Society of Petroleum Engineers. 1995
  102. 102. Petrosky JG, Farshad F. Pressure-volume-temperature correlations for Gulf of Mexico crude oils. SPEREE. 1998;1(5):416-420. SPE-51395-PA. DOI: 10.2118/51395-PA
  103. 103. Rollins JB, McCain WD Jr, Creeger TJ. Estimation of solution GOR of black oils. Journal of Petroleum Technology. 1990;42(01):92-94
  104. 104. Sattarin M, Modarresi H, Bayat M, Teymori M. New viscosity correlations for dead crude oils. Petroleum & Coal. 2007;49(2):33-39
  105. 105. Spivey JP, Valkó PP, McCain WD. Applications of the coefficient of isothermal compressibility to various reservoir situations with new correlations for each situation. SPE Reservoir Evaluation & Engineering. 2007;10(01):43-49
  106. 106. Standing MB. Volumetric and Phase Behavior of Oil Field Hydrocarbon Systems: PVT for Engineers. California Research Corporation; 1951
  107. 107. Standing M. Volumetric and Phase Behavior of Oil Field Hidrocarbon Systems. Richardson, Texas: SPE; 1981
  108. 108. Sim SSK. Pressure-volume-temperature correlations for crude oils from the Illinois Basin. In: Illinois Petroleum. University of Illinois Urbana-Champaign. Vol. 140; 1993
  109. 109. Sulaimon A, Ramli N, Adeyemi B, Saaid I. New correlation for oil formation volume factor. In: SPE Nigeria Annual International Conference and Exhibition; Society of Petroleum Engineers. 2014
  110. 110. Twu CH. Internally consistent correlation for predicting liquid viscosities of petroleum fractions. Industrial & Engineering Chemistry Process Design and Development. 1985;24(4):1287-1293
  111. 111. Valko P, McCain W Jr. Reservoir oil bubblepoint pressures revisited; solution gas–oil ratios and surface gas specific gravities. Journal of Petroleum Science and Engineering. 2003;37(3–4):153-169
  112. 112. Vazquez M, Beggs H. Correlations for fluid physical property prediction. JPT. 1980;32(6):968-970. SPE-6719-PA. DOI: 10.2118/6719-PA
  113. 113. Velarde J, Blasingame T, McCain W Jr. Correlation of black oil properties at pressures below bubble point pressure—A new approach. In: Annual Technical Meeting; Petroleum Society of Canada. 1997

Written By

Abdelaziz El-Hoshoudy and Saad Desouky

Submitted: 28 June 2018 Reviewed: 26 October 2018 Published: 18 December 2019