PNG 520
Phase Relations in Reservoir Engineering

Peng-Robinson EOS (1976)

PrintPrint

The Peng-Robinson EOS has become the most popular equation of state for natural gas systems in the petroleum industry. During the decade of the 1970’s, D. Peng was a PhD student of Prof. D.B. Robinson at the University of University of Alberta (Edmonton, Canada). The Canadian Energy Board sponsored them to develop an EOS specifically focused on natural gas systems. When you compare the performance of the PR EOS and the SRK EOS, they are pretty close to a tie; they are “neck to neck,” except for a slightly better behavior by the PR EOS at the critical point. A slightly better performance around critical conditions makes the PR EOS somewhat better suited to gas/condensate systems.

Peng and Robinson introduced the following modified vdW EOS:

( P+ αa v ˜ 2 +2b v ˜ − b 2 )( v ˜ −b )=RT This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.1a)
 

or, explicitly in pressure,

P= RT v ˜ −b − αa v ˜ 2 +2b v ˜ − b 2 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.1b)
 

where:

Peng and Robinson conserved the temperature dependency of the attractive term and the acentric factor introduced by Soave. However, they presented different fitting parameters to describe this dependency (equation 4.11c), and further manipulated the denominator of the pressure correction (attractive) term. As we have seen before, coefficients “a” and “b” are made functions of the critical properties by imposing the criticality conditions. This yields:

α= ⌊ 1+( 0.37464+1.54226ω−0.26992 ω 2 )( 1− T r ) ⌋ 2 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.1c)
 
a=0.45724 R 2 T c 2 P c This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.2a)
 

The PR cubic expression in Z becomes:

b=0.07780 R T c P c This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.2b)
 

where:

Z 3 −( 1−B ) Z 2 +( A−2B−3 B 2 )Z−( AB− B 2 − B 3 )=0 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.3a)
 
A= αaP R 2 T 2 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.3b)
 
B= bP RT This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.3c)
 

Similar to SRK, the PR mixing rules are:

( αa ) m = ∑ ​ ∑ ​ y i y j ( αa ) ij  ;  ( αa ) ij = ( αa ) i ( αa ) j ( 1− k ij ) This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.4a)
 
b m = ∑ i y i b i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(11.4b)
 

where binary interaction parameters (kij) again play the important empirical role of helping to better fit experimental data. Due to the empirical character of these interaction parameters, kij’s calculated for PR EOS are unlikely to be the same as the kij’s calculated for SRK EOS for the same pair of molecules.