Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information

# Difference between revisions of "Calculating PVT properties"

This is an example of calculating PVT properties. The specific correlations that should be used for a specific crude oil or reservoir may vary, as discussed in the referenced pages focusing on specific properties.

Determine the PVT properties for a United States midcontinental crude oil and natural gas system with properties listed in Table 1. Table 2 lists the correlations to be used. Measured data are provided for comparison with the calculated results. For correlations that rely on other correlations, these data illustrate the effects of error propagation in the calculations.

## Gravity and molecular weight

Determine the crude oil specific gravity, ....................(1)

and molecular weight, ....................(2)

## Bubblepoint pressure

Use the Lasater correlation to estimate bubblepoint pressure. Calculate the gas mole fraction in the oil, ....................(3)

and the Lasater bubblepoint pressure factor, ....................(4)

with Lasater’s relationship between bubblepoint pressure factor and bubblepoint pressure, ....................(5)

For comparison, Standing = 2,316 psia, Glasø = 2,725 psia, Al-Shammasi = 2,421 psia, and Velardi = 2,411 psia.

Modify the calculated bubblepoint pressure to account for the effects of nitrogen in the surface gas with Jacobson’s equation. ....................(6)

Therefore, the bubblepoint pressure should be increased by 9.8% to 2,251 psia. The measured bubblepoint pressure was reported to be 2,479 psia.

## Bubblepoint oil formation volume factor

Calculate the bubblepoint oil formation volume factor (FVF) using the correlation from Al-Shammasi. ....................(7)

For comparison (in bbl/STB), Standing = 1.410, Glasø = 1.386, Al-Marhoun = 1.364, Farshad = 1.364, and Kartoatmodjo = 1.358. The measured bubblepoint oil FVF is 1.398 bbl/STB.

## Isothermal compressibility

Calculate the isothermal compressibility of oil using the Farshad correlation. ....................(8) ....................(9)

The measured isothermal compressibility is 11.06 × 10-6psi-1.

## Undersaturated oil formation volume factor

With the results from Lasater’s method for bubblepoint pressure, use Al-Shammasi’s method for bubblepoint oil FVF, and Farshad’s equation for isothermal compressibility, the undersaturated oil FVF is given by ....................(10)

which compares to a measured value of 1.367 bbl/STB. Because this calculation uses the results from multiple correlations, individual correlation error compounds and propagates through to the final result. The calculated value is 1.367 bbl/STB with the actual bubblepoint value of 1.398 bbl/STB; therefore, the accuracy of the bubblepoint FVF is primarily affected by the accuracy of the undersaturated FVF.

## Oil density

Calculate the oil density. ....................(11)

Calculate the dead oil viscosity using the correlation from Glasø. ....................(12)

For comparison, Fitzgerald = 1.808 cp, and Bergman = 2.851 cp. The measured dead oil viscosity is 1.67 cp.

## Bubblepoint oil viscosity

Calculate the bubblepoint oil viscosity using the method developed by Chew and Connally. ....................(13) ....................(14) ....................(15)

For comparison, Beggs and Robinson = 0.515 cp. The measured viscosity at bubblepoint is 0.401 cp.

## Undersaturated oil viscosity

Calculate the undersaturated oil viscosity by applying the Vazquez and Beggs correlation. ....................(16)

For comparison, Beal = 0.730 cp and Kouzel = 0.778 cp. The measured value is 0.475 cp. This example illustrates the steps necessary to calculate oil viscosity requiring correlations for dead oil viscosity, bubblepoint viscosity, undersaturated viscosity, and bubblepoint pressure/solution GOR. Errors in individual correlations can compound and propagate through to the resulting answer. For instance, if the measured bubblepoint viscosity is used in Eq. 16, the result is 0.52 cp—much closer to the measured value. Therefore, care should be exercised in the selection of accurate correlations for individual properties.

## Gas/oil interfacial tension

Estimate the gas/oil surface tension using the method developed by Abdul-Majeed. Calculate the dead oil surface tension. ....................(17)

Determine the live oil adjustment factor. ....................(18)

Calculate the live gas/oil surface tension. ....................(19)

For comparison, Baker and Swerdloff = 4.73 dynes/cm.

## Water/oil interfacial tension

Estimate the water/oil surface tension using Firoozabadi and Ramey. Calculate the pseudocritical temperature of the dead oil. ....................(20)

Calculate the pseudocritical temperature of the gas. ....................(21)

Calculate the pseudocritical temperature of the live gas/oil mixture. ....................(22)

Convert oil density units from lbm/ft3 to g/cm3. ....................(23)

Calculate the surface tension between the oil and water phases. ....................(24)

## Nomenclature

 Bg = gas FVF, ft3/scf Bo = oil FVF, bbl/STB Bob = oil formation volume at bubblepoint pressure, bbl/STB co = oil isothermal compressibility, Lt2/m, psi-1 cob = oil isothermal compressibility at bubblepoint, Lt2/m, psi-1 Kw = Watson characterization factor, °R1/3 Mg = gas molecular weight, m, lbm/lbm mol Mgo = gas/oil mixture molecular weight, m, lbm/lbm mol Mo = oil molecular weight, m, lbm/lbm mol Mog = oil-gas mixture molecular weight, m, lbm/lbm mol p = pressure, m/Lt2, psia pb = bubblepoint pressure, m/Lt2, psia = bubblepoint pressure of oil with N2 present in surface gas, m/Lt2, psia pbh = bubblepoint pressure of oil without nonhydrocarbons, m/Lt2, psia pf = bubblepoint pressure factor, psia/°R pr = pressure ratio (fraction of bubblepoint pressure) Rs = solution GOR, scf/STB T = temperature, T, °F Tabs = temperature, T, °R Tb = mean average boiling point temperature, T, °R Tcg = gas pseudocritical temperature, T, °R Tcm = mixture pseudocritical temperature, T, °R Tco = oil pseudocritical temperature, T, °R Tr = reduced temperature, T Tsc = temperature at standard conditions, T, °F V = volume, L3 Vo = volume of crude oil, L3 Wg = weight of dissolved gas, m Wo = weight of crude oil, m xg = gas "component" mole fraction in oil xo = oil "component" mole fraction in oil yg = gas "component" mole fraction in gas = mole fraction N2 in surface gas yo = oil "component" mole fraction in gas Z = gas compressibility factor γAPI = oil API gravity γg = gas specific gravity, air=1 γgc = gas specific gravity adjusted for separator conditions, air=1 γghc = gas specific gravity of hydrocarbon components in a gas mixture, air=1 γgs = separator gas specific gravity, air=1 γo = oil specific gravity μo = oil viscosity, m/Lt, cp μob = bubblepoint oil viscosity, m/Lt, cp μod = dead oil viscosity, m/Lt, cp ρg = gas density, m/L3, lbm/ft3 ρo = oil density, m/L3, lbm/ft3 ρob = bubblepoint oil density, m/L3, lbm/ft3 ρw = water density, m/L3, g/cm3 σhw = hydrocarbon/water surface tension, m/t2, dynes/cm σgo = gas/oil surface tension, m/t2, dynes/cm σod = dead oil surface tension, m/t2, dynes/cm