Equations of State and

PVT Analysis
Single-Component Systems

 A una determinada temperatura, la

presión de rocío y la presión de burbuja
son iguales.
 La curva de punto de rocío representa los
estados de gas saturado.
 La curva de punto de burbuja representa
los estados de líquido saturado.
 El punto donde se unen ambas curvas se
conoce como punto crítico (Tc, Pc, Vc)
 Por fuera de la envolvente de fase, solo
una fase puede existir.
 El punto crítico es el punto donde todas
las propiedades intensivas de las fases
líquida y gaseosa son iguales.
 𝑝𝑣 = 𝑝𝑐 exp 𝐴 + 𝜔𝐵

6,09648 6
𝐴 = 5,92714 − − 1,2886ln 𝑇𝑟 + 0,16934 𝑇𝑟
𝑇𝑟

15,6875 6
𝐵 = 15,2518 − − 13,4721ln 𝑇𝑟 + 0,4357 𝑇𝑟
𝑇𝑟

𝑝𝑣
𝜔 = − log −1
𝑝𝑐 𝑇=0,7𝑇𝑐

Fig 2. Diagrama Presión Temperatura para un sistema single-component.

• El factor acéntrico, , es un parámetro de
correlación empleado para caracterizar la
La curva de presión de vapor representa las centricidad o no-esfericidad de una molécula.
condiciones de presión y temperatura a las
cuales dos fases, vapor y líquido, pueden
coexistir en equilibrio.
Two-Component Systems

 Para un sistema single-component dos fases

pueden coexistir en equilibrio, a determinada
temperatura, a solo una presión, la presión de
vapor.
 En un sistema binario, dos fases pueden existir
en equilibrio a varias presiones a la misma
temperatura.
 Una característica de sistemas binarios es la
variación de sus propiedades termodinámicas y
físicas con la composición.
 En el punto de rocío, la composición total del
sistema es igual a la composición de la fase
gaseosa.
 En el punto de burbuja, la composición total del
sistema es igual a la composición de la fase
líquida.
 La líneas de calidad describen condiciones de
presión y temperatura de igual volumen de
líquido.
 El punto A de la curva de presión vs
composición, representa la presión de vapor
del componente más volátil.

 Las composiciones de líquido y vapor

coexistentes se hayan proyectando los
extremos de la línea recta a través de la
región bifásica hasta los ejes de composición.
Clasificación de Yacimientos y Fluidos de Yacimiento

 Oil Reservoirs: Si la temperatura del yacimiento es menor a la temperatura crítica del fluido,
se clasifica como un yacimiento de aceite.
 Gas Reservoirs: Si la temperatura del yacimiento es mayor a la temperatura crítica del fluido,
entonces se clasifica como un yacimiento de gas.

Oil Reservoirs

 De acuerdo a la presión de yacimiento pueden:

ser subsaturados, saturados y con capa de gas.
 Con respecto a sus propiedades físicas, se
clasifican en:
 Ordinary Black Oil
 Low-Shrinkage oil
 High Shrinkage (volátil) oil
 Near critical crude oil
Gas Reservoirs
 Retrograde gas Reservoirs
 Near-critical gas-condensate reservoirs
 Wet gas Reservoirs
 Dry gas Reservoirs

 Comportamiento de yacimientos de Gas Retrógrado

 A medida que la presión de yacimiento decrece durante la producción desde el punto inicial
hasta el punto de rocío, la atracción entre las moléculas de componentes livianos y pesados se
debilita.
 Al mismo tiempo, la atracción entre moléculas de componentes pesados se hace más efectiva,
por tanto, el líquido comienza a condensarse.
 Este comportamiento continúa hasta que el “liquid dropout” alcance un punto máximo.
 Una adicional reducción de la presión permite a las moléculas pesada comenzar un proceso de
vaporización normal.
 En la mayoría de yacimientos el condensado rara vez excede el 15-19% del volumen poroso.
 Natural Gas Properties
Comportamiento de Gases Ideales
 El volumen de las moléculas es insignificante comparado con el volumen total ocupado.
 Las moléculas no presentan fuerzas de atracción o repulsión entre ellas.
 Las colisiones entre moléculas son perfectamente elásticas.

Ley de gases ideales

P, presión absoluta, psia
V, volumen, ft3
𝑝𝑉 = 𝑛𝑅𝑇 T, temperatura absoluta, °R
n, número de moles de gas, lb-mol
R, constante de gases, (10,73 psia ft3/lb-mol °R)

n está definido como la masa del gas dividida entre el peso molecular. Combinando con ley
de gases, se puede calcular un valor de densidad, lb/ft3
Peso Molecular Aparente, 𝑴𝒂 𝑀𝑎 = ෍ 𝑦𝑖 𝑀𝑖
𝑖=1

1 𝑅𝑇 (1)(10,73)(520)
𝑉𝑠𝑐 = =
𝑝𝑠𝑐 14,7
Volumen Estándar, 𝑽𝒔𝒄
𝑉𝑠𝑐 = 379,4 𝑠𝑐𝑓/𝑙𝑏𝑚𝑜𝑙𝑒

𝑝𝑀𝑎
Densidad del Gas, 𝝆𝒈 𝜌𝑔 =
𝑅𝑇

𝑀𝑎 𝑀𝑎
Gravedad específica, 𝜸𝒈 𝛾𝑔 = =
𝑀𝑎𝑖𝑟 28,96
Comportamiento de Gases Reales
 La magnitud de la desviación de los gases reales del comportamiento ideal aumenta con el aumento en la presión y
la temperatura.
𝑉𝑟𝑒𝑎𝑙
𝑝𝑉 = 𝑍𝑛𝑅𝑇 𝑍=
𝑉𝑖𝑑𝑒𝑎𝑙
 El factor de compresibilidad puede generalizarse cuando se expresa en
términos de las propiedades pseudo-reducidas, con las cuales se ingresa
a la Carta de Factor de Compresibilidad de Standing and Katz.
𝑝
𝑝𝑝𝑟 = 𝑝𝑝𝑐 = ෍ 𝑦𝑖 𝑝𝑐𝑖
𝑝𝑝𝑐
𝑖=1

𝑇
𝑇𝑝𝑟 = 𝑇𝑝𝑐 = ෍ 𝑦𝑖 𝑇𝑐𝑖
𝑇𝑝𝑐
𝑖=1

𝑝𝑐 𝐶7+ = 1188 − 431 log 𝑀𝐶7+ − 61,1 + 2319 − 852 log 𝑀𝐶7+ − 53,7 𝛾𝐶7+ − 0,8

𝑇𝑐 𝐶7+ = 608 − 364 log 𝑀𝐶7+ − 71,2 + 2450 log 𝑀𝐶7+ − 3800 𝑙𝑜𝑔 𝛾𝐶7+
 En casos donde la composición del gas natural no está disponible, las propiedades pseudo-críticas
pueden ser predichas a partir de la gravedad específica

Caso 1: Gas natural

𝑝𝑝𝑐 = 677 + 15𝛾𝑔 − 37,5𝛾𝑔2

𝑇𝑝𝑐 = 168 + 325𝛾𝑔 − 12,5𝛾𝑔2

Caso 2: Gas húmedo

𝑝𝑝𝑐 = 706 − 51,7𝛾𝑔 − 11,1𝛾𝑔2

𝑇𝑝𝑐 = 187 + 330𝛾𝑔 − 71,5𝛾𝑔2

Desde la página 155 se proponen varias correlaciones para

calcular el factor de compresibilidad.
Factor Volumétrico de Formación, 𝑩𝒈 : Utilizado para relacionar el volume de gas, medido a condiciones
de yacimiento, con el volume de gas medido a condiciones estándar, 60°F y 14,7psia

𝑍𝑇 𝑓𝑡 3൘
𝐵𝑔 = 0,02827 , 𝑠𝑐𝑓
𝑍𝑛𝑅𝑇 𝑝
𝑉 𝑝,𝑇 𝑝 𝑝𝑠𝑐 𝑍𝑇
𝐵𝑔 = = =
𝑉𝑠𝑐 𝑍𝑠𝑐 𝑛𝑅𝑇𝑠𝑐 𝑇𝑠𝑐 𝑝
𝑍𝑇 𝑏𝑏𝑙ൗ
𝑝𝑠𝑐 𝐵𝑔 = 0,005035 ,
𝑝 𝑠𝑐𝑓

Factor de Expansión del Gas, 𝑬𝒈

𝑝 𝑠𝑐𝑓
𝐸𝑔 = 35,37 , ൘𝑓𝑡 3
𝑍𝑇
1
𝐸𝑔 =
𝐵𝑔
𝑝 𝑠𝑐𝑓ൗ
𝐸𝑔 = 198,6 , 𝑏𝑏𝑙
𝑍𝑇
 PVT Properties of Crude Oils
Crude Oil Gravity

𝜌𝑜 𝜌𝑜 141,5 141,5
𝛾𝑜 = = 𝐴𝑃𝐼 = − 131,5 𝛾𝑜 =
𝜌𝑤 62,4 𝛾𝑜 𝐴𝑃𝐼 + 131,5

Crude Oil Density

 Con base en la composición, se propone el desarrollo de Standing and Katz (página 184)
 Basado en datos PVT, Standing (1981) propone la siguiente correlación para calcular la
densidad, donde la temperatura está en °R.

62,4𝛾𝑜 + 0,0136𝑅𝑠 𝛾𝑔
𝜌𝑜 = 1,175
𝛾𝑔 0,5
0,972 + 0,000147 𝑅𝑠 + 1,25(𝑇 − 460)
𝛾𝑜
 Ahmed expresó el peso molecular aparente del petróleo crudo con su gas en solución disuelto a
través de la siguiente expresión.

 𝑀𝑠𝑡 , peso molecular del STO y puede ser tomado

como el peso molecular de la fracción heptano-plus.
0,0763𝑅𝑠 𝛾𝑔 𝑀𝑠𝑡 + 350,376𝛾𝑜 𝑀𝑠𝑡
𝑀𝑎 =
0,0026537𝑅𝑠 𝑀𝑠𝑡 + 350,376𝛾𝑜
 𝛾𝑜 , gravedad específica del STO o el 𝐶7+

 La densidad del aceite a condiciones estándar se puede determinar, cuando el peso molecular no
está disponible, a partir de la siguiente correlación

0,0763𝑅𝑠 𝛾𝑔 + 350,4𝛾𝑜
𝜌𝑠𝑐 =
0,0027𝑅𝑠 + 2,4893𝛾𝑜 + 3,491

 A partir de la página 200 se presentan algunas correlaciones para calcular la solubilidad del gas.
Oil Formation Volume Factor

 Siempre es mayor o igual a la unidad

 Relación entre volumen a condiciones de
yacimiento y volumen a condiciones estándar.
 Correlaciones a partir de la página 213 del
libro.

Total Formation Volume Factor

𝑉𝑜 𝑝,𝑇 + 𝑉𝑔
𝑝,𝑇
𝐵𝑡 = = 𝐵𝑜 + 𝑅𝑠𝑏 − 𝑅𝑠 𝐵𝑔
𝑉𝑜 𝑠𝑐

 Correlaciones a partir de la página 234.

 A partir de la página 218 se presentan correlaciones para calcular la compresibilidad isotérmica, definida
como la tasa de cambio en volumen respecto a un incremento en la presión.

 A partir de la página 237 se presentan correlaciones para calcular la Viscosidad del Petróleo Crudo, antes y
después del punto de burbuja.
Laboratory Analysis of Reservoir Fluids

 Routine laboratory tests

 Análisis composicional del sistema
 Constant Composition Expansion (CCE)
 Liberación Diferencial
 Prueba de Separador
 Constant Volume Depletion (CVD)

 Special laboratory tests

 Slim-tube test
 Swelling test
 Aplicadas si el yacimiento se va a
depletar bajo inyección miscible de gas.
 Constant Composition Expansion Test

 Conducida para simular la relación presión/volumen de un sistema hidrocarburo.

 Permite determinar: presión de saturación, coeficiente de compresibilidad isotérmica, factor de
compresibilidad de la fase gaseosa y el volumen total de hidrocarburo.
 El volumen de hidrocarburo como función de la presión de la celda se reporta como una relación del
volumen de referencia:

𝑉𝑡 𝑝𝑠𝑎𝑡 − 𝑝 𝑝𝑠𝑎𝑡 − 𝑝
𝑉𝑟𝑒𝑙 = 𝑌= 𝑉𝑟𝑒𝑙 = 1 +
𝑉𝑠𝑎𝑡 𝑝 𝑉𝑟𝑒𝑙 − 1 𝑝(𝑎 + 𝑏𝑝)

𝑌 = 𝑎 + 𝑏𝑝

 La función Y es utilizada para suavizar los valores de

volumen relativo.
 Cuando se grafica presenta la forma de una línea recta o
una pequeña curvatura.
 Differential Liberation (Vaporization) Test

 El gas en solución liberado durante cada reducción de

presión es continuamente removido del contacto con el
aceite.
 De la prueba se pueden obtener los siguientes datos:
 Cantidad de gas en solución como función de la
presión
 El “shrinkage” del aceite como función de la presión
 Composición, Z-factor y gravedad específica del gas
liberado
 Densidad del aceite como función de la presión.

 Separator Test
 Conducidas para determinar cambios en el comportamiento
volumétrico del fluido a medida que pasa a través de
separadores y el tanque.
 El comportamiento volumétrico resultante se ve influenciado
en gran medida por las condiciones de operación (presión y
temperatura) de las facilidades de superficie.
 Slim Tube Test
 Diseñada para medir la Mínima Presión de Miscibilidad
 Permite evaluar datos para estudios específicos, como vaporización de aceite por inyección de gas.
 Para examinar la mezcla de fluidos durante un proceso de desplazamiento miscible.
 Los resultados de la prueba se interpretan graficando el recobro acumulado de aceite versus el volumen
poroso de gas inyectado.

 MMP incrementa con la temperatura de yacimiento.

 MMP es menor si el porcentaje de intermedios en la
fase aceite es mayor.
 MMP alto para aceites con alta concentración de C1
y N2.
 Equation of State and Phase Equilibria
 Equilibrium Ratios
Definidos como la relación entre la fracción molar del componente en la fase gaseosa y la fracción del
componente en la fase líquida.

𝑦𝑖 𝑦𝑖 𝑝𝑣𝑖
𝐾𝑖 = 𝑝𝑖 = 𝑥𝑖 𝑝𝑣𝑖 𝑝𝑖 = 𝑦𝑖 𝑝 = = 𝐾𝑖
𝑥𝑖 𝑥𝑖 𝑝
Raoult’s Law Dalton’s Law

𝑧𝑖
𝑥𝑖 = ෍ 𝑥𝑖 = 1
𝑛𝐿 + 𝑛𝑉 𝐾𝑖 𝑖
𝑧𝑖 𝐾𝑖 − 1 𝑧𝑖 𝐾𝑖 − 1
෍ =0 𝑓 𝑛𝑉 = ෍ =0
෍ 𝑦𝑖 = 1 𝑛𝐿 + 𝑛𝑉 𝐾𝑖 𝑛𝑉 𝐾𝑖 − 1 + 1
𝑖 𝑖
𝑖
𝑧𝑖 𝐾𝑖
𝑦𝑖 =
𝑛𝐿 + 𝑛𝑉 𝐾𝑖 ෍ 𝑦𝑖 − ෍ 𝑥𝑖 = 0
𝑖 𝑖
 Flash Calculations
A continuación se resumen los pasos para determinar las composiciones en líquido y vapor a
partir de una composición general a partir de un conjunto de constantes de equilibrio
 Asuma un valor arbitrario de 𝑛𝑉 entre 0 y 1. Un buen valor se puede calcular de la
siguiente manera:

𝐴 1
𝑛𝑉 = 𝐴 = ෍ 𝑧𝑖 𝐾𝑖 − 1 𝐵 = ෍ 𝑧𝑖 −1
𝐴+𝐵 𝐾𝑖
𝑖 𝑖

Evalue la función f 𝑛𝑉 . Si el resultado es menor que una tolerancia establecida, como

10−6 , entonces 𝑛𝑉 asumido es correcto. De lo contrario, calculi un nuevo valor de 𝑛𝑉 a
partir de la siguiente correlación:

𝑓 𝑛𝑉 𝑧𝑖 𝐾𝑖 − 1 2
𝑛𝑉 𝑛𝑒𝑤 = 𝑛𝑉 − 𝑓′ 𝑛𝑉 = −෍ 2
𝑓′ 𝑛𝑉 𝑛𝑉 𝐾𝑖 − 1 + 1
𝑖

Repetir proceso hasta que se alcance la convergencia, f 𝑛𝑉 < 10−6

  Calcule el número de moles de líquido teniendo en cuenta que 𝑛𝑉 + 𝑛𝐿 = 1.
 Calcule la composición de la fase líquida y la fase gaseosa aplicando las ecuaciones descritas
anteriormente.

 Equilibrium Ratios for Real Solutions

Se presentan correlaciones para calcular las constantes de equilibrio, como la de Standing (página 340)
y la de Wilson, descrita como sigue:

𝑝𝑐𝑖 𝑇𝑐𝑖
𝐾𝑖 = 𝑒𝑥𝑝 5,37 1 + 𝜔𝑖 1−
𝑝 𝑇

Los valores de las constantes de equilibrio de todos los componentes del sistema, convergen a un valor
común de uno, a una cierta presión, denominada presión de convergencia.

Para componentes No-Hidrocarburos se pueden calcular las constantes de equilibrio como sigue:

0,8
𝑝 1399,2204 18,215052 1112446,2
𝑙𝑛 𝐾𝐻2𝑆 = 1− 6,3992127 + − ln 𝑝 0,76885112 + −
𝑝𝑘 𝑇 𝑇 𝑇2
 0,4
𝑝 1184,2409
𝑙𝑛 𝐾𝑁2 = 1− 11,294748 + − 0,90459907 ln 𝑝
𝑝𝑘 𝑇

0,6
𝑝 152,7291 1719,2956 644740,69
𝑙𝑛 𝐾𝐶𝑂2 = 1− 7,0201913 + − ln 𝑝 1,8896974 + +
𝑝𝑘 𝑇 𝑇 𝑇2

 Vapour-Liquid Equilibrium Calculations

 Dew Point Pressure

Tener en cuenta que la composición del sistema es igual a la composición de la fase gaseosa,
por tanto, 𝑛𝑙 = 0 𝑛𝑣 = 1

 Asumir un valor de 𝑝𝑑 . Un buen valor inicial se puede calcular empleando la siguiente

ecuación.
1
𝑝𝑑 =
𝑧
σ𝑖 𝑖
𝑝𝑣𝑖
 Usando la presión de rocío asumida, calcule las constantes de equilibrio para cada
componente, a la temperatura del sistema. Luego calcule la suma, la cual se obtiene aplicando
condiciones a la ecuación de composición.
𝑧𝑖
෍
𝐾𝑖
𝑖

 Si la suma es menor que 1, repita procedimiento asumiendo un valor más alto de presión
inicial; si la suma es mayor que 1, repita asumiendo un valor menor de presión inicial. El valor
correcto de punto de rocío se obtiene cuando la sumatoria da igual a la unidad.

 Mirar y desarrollar ejemplo presentado en página 353.

 Bubble Point Pressure

 Procedimiento similar

𝑝𝑏 = ෍ 𝑧𝑖 𝑝𝑣𝑖 ෍ 𝑧𝑖 𝐾𝑖 = 1
𝑖 𝑖
 Separator Calculations

 Cada condición de presión y temperatura a la cual las fases de hidrocarburo son flasheadas, se
denomina etapa de separación.
 Existen dos tipos de separación gas/aceite: diferencial y flash.
 En la separación diferencial el gas liberado es removido del contacto con el aceite a medida que la
presión es reducida. Cuando el gas es separado en esta manera, la máxima cantidad de
componentes intermedios y pesados permanecen en la fase líquida, por tanto, hay un menor
“shrinkage” del aceite y una mayor recuperación de líquido en el tanque.
 En la separación flash, el gas liberado permanece en contacto con el aceite hasta que es
instantáneamente removido en la etapa final de separación. Una máxima proporción de
componentes intermedios y pesados son atraídos a la fase gaseosa con este procedimiento.
 Si la presión de separador es alta, gran cantidad de componentes livianos permanece en la fase
líquida y se pierde junto a otros componentes valiosos en el tanque.
 Si la presión es muy baja, gran cantidad de componentes livianos se separan y llevan consigo
cantidades substanciales de componentes intermedios y pesados.
 Ecuaciones de Estado (EOS)

 Redlich-Kwong
𝑅𝑇 𝑎 𝑅2 𝑇𝑐2,5 𝑅𝑇𝑐
𝑝= − 𝑎 = 0,42747 𝑏 = 0,08664
𝑉 − 𝑏 𝑉(𝑉 + 𝑏) 𝑇 𝑝𝑐 𝑝𝑐

𝑎𝑝 𝑏𝑝
𝑍3 − 𝑍2 + 𝐴−𝐵− 𝐵2 𝑍 − 𝐴𝐵 = 0 𝐴 = 2 2,5 𝐵=
𝑅 𝑇 𝑅𝑇

 La ecuación arroja una raíz real en la región de una fase, y tres raíces reales en la región
bifásica. En el segundo caso, la raíz de mayor valor corresponde al factor de compresibilidad
de la fase gas y la de menor valor corresponde al factor de compresibilidad de la fase líquida.

𝑛 2 𝑛

𝑎𝑚 = ෍ 𝑦𝑖 𝑎𝑖 𝑏𝑚 = ෍ 𝑦𝑖 𝑏𝑖
𝑖=1 𝑖=1

 Lo mismo para la fase líquida, cambiando 𝑦𝑖 por 𝑥𝑖

 Soave-Redlich-Kwong

𝑅𝑇 𝑎𝛼 𝑇 2
𝑝= − 𝛼 𝑇 = 1 + 𝑚 1 − 𝑇𝑟 𝑚 = 0,480 + 1,574𝜔 − 0,17𝜔2
𝑉 − 𝑏 𝑉(𝑉 + 𝑏)

𝑅2 𝑇𝑐2 𝑅𝑇𝑐
𝑎 = 0,42747 𝑏 = 0,08664
𝑝𝑐 𝑝𝑐

(𝑎𝛼)𝑝 𝑏𝑝
𝑍 3 − 𝑍 2 + 𝐴 − 𝐵 − 𝐵2 𝑍 − 𝐴𝐵 = 0 𝐴= 2 2 𝐵=
𝑅 𝑇 𝑅𝑇

(𝑎𝛼)𝑚 = ෍ ෍ 𝑥𝑖 𝑥𝑗 𝑎𝑖 𝑎𝑗 𝛼𝑖 𝛼𝑗 (1 − 𝑘𝑖𝑗 ) 𝑏𝑚 = ෍ 𝑥𝑖 𝑏𝑖
𝑖 𝑗 𝑖
Concepto de Fugacidad

 Es una medida de la energía molar de Gibbs de un gas real.

 Tiene unidades de presión y puede ser vista como la presión de vapor modificada para
correctamente representar la tendencia de las moléculas de una fase a escapar a la otra.
 Para un componente puro, Soave aplicó la siguiente correlación:

0 𝑓 0 , fugacidad de un componente puro, psia

𝑓 𝐴 𝑍+𝐵
𝑙𝑛 = 𝑙𝑛 Φ = 𝑍 − 1 − ln 𝑍 − 𝐵 − ln
𝑝 𝐵 𝑍
Φ, coeficiente de fugacidad

 La fase con la menor fugacidad del componente, acepta el componente de la fase con la mayor
fugacidad del componente. Fugacidades iguales de un componente en ambas fases resulta en cero
transferencia.
 Cero transferencia para todos lo componentes implica que el sistema se encuentra en equilibrio
termodinámico.

𝑓𝑖
𝑣
𝑓𝑖
𝐿 Φ𝑖𝐿
Φ𝑖𝑣 = Φ𝑖𝐿 = 𝐾𝑖 = 𝑣
𝑦𝑖 𝑝 𝑥𝑖 𝑝 Φ𝑖
 Soave propone la siguiente correlación para calcular el coeficiente de fugacidad el componente i en
la fase líquida.

𝑓𝑖𝐿 𝐿 𝑏𝑖 𝑍 𝐿−1
𝐿
𝐴 2Ψ𝑖 𝑏𝑖 𝐵
𝑙𝑛 = 𝑙𝑛 Φ𝑖 = − 𝑙𝑛 𝑍 − 𝐵 − − 𝑙𝑛 1 + 𝐿
𝑥𝑖 𝑝 𝑏𝑚 𝐵 (𝑎𝛼)𝑚 𝑏𝑚 𝑍

Ψ𝑖 = ෍ 𝑥𝑗 𝑎𝑖 𝑎𝑗 𝛼𝑖 𝛼𝑗 1 − 𝑘𝑖𝑗 (𝑎𝛼)𝑚 = ෍ ෍ 𝑥𝑖 𝑥𝑗 𝑎𝑖 𝑎𝑗 𝛼𝑖 𝛼𝑗 (1 − 𝑘𝑖𝑗 )

𝑗 𝑖 𝑗

 Mismo procedimiento para la fase vapor, cambiando 𝑥𝑖 por 𝑦𝑖 .

 En la página 384 se presenta un procedimiento para seleccionar la raíz correcta para una fase en caso
de que la se tengan tres raíces reales.
 Peng-Robinson

𝑅𝑇 𝑎𝛼 𝑅𝑇 𝑎𝛼 Ecuación optimizada para

𝑝= − 𝑝= −
𝑉 − 𝑏 (𝑉 + 𝑏)2 −𝑐𝑏 2 𝑉 − 𝑏 𝑉 𝑉 + 𝑏 + 𝑏(𝑉 − 𝑏) un valor de c=2 y Zc=0,307
 Si 𝜔 > 0,49

𝑚 = 0,3796 + 1,5422𝜔 − 0,2699𝜔2 𝑚 = 0,379642 + 1,48503𝜔 − 0,1644𝜔2 + 0,016667𝜔3

𝑅2 𝑇𝑐2 𝑅𝑇𝑐 𝛼 = 1 + 𝑚 1 − 𝑇𝑟
2
𝑎 = 0,45724 𝑏 = 0,07780
𝑝𝑐 𝑝𝑐

(𝑎𝛼)𝑚 𝑝 𝑏𝑚 𝑝 (𝑎𝛼)𝑚 = ෍ ෍ 𝑥𝑖 𝑥𝑗 𝑎𝑖 𝑎𝑗 𝛼𝑖 𝛼𝑗 (1 − 𝑘𝑖𝑗 ) 𝑏𝑚 = ෍ 𝑥𝑖 𝑏𝑖

𝐴= 2 2 𝐵=
𝑅 𝑇 𝑅𝑇 𝑖 𝑗 𝑖

𝑍 3 − 𝐵 − 1 𝑍 2 + 𝐴 − 3𝐵2 − 2𝐵 𝑍 − (𝐴𝐵 − 𝐵2 − 𝐵3 ) = 0

 En la página 392 se presenta una tabla con coeficientes de interacción binaria.

 Equation of State Applications
 Equilibrium Ratios

 Asuma valores iniciales para las relaciones

de equilibrio, 𝐾𝑖 . Puede emplearse la
correlación de Wilson para tal fin.
 Use la composición general y las constantes
para determinar 𝑥𝑖 , 𝑦𝑖 , 𝑛𝑣 , 𝑛𝐿
 Determine los coeficientes de fugacidad
para cada componente en ambas fases y
calcule nuevos valores de relaciones de
equilibrio, 𝐾𝑖𝐴
𝑛 2
𝐾𝑖
෍ 𝐴−1 ≤ 𝜀 = 0,0001
𝐾𝑖
𝑖=1

 Si no se cumple lo anterior, repita

procedimiento haciendo 𝐾𝑖 =𝐾𝑖𝐴
 Dew Point Pressure, 𝑷𝒅
 Asuma un valor inicial de punto de rocío
 Usando el valor asumido calcule las constantes de equilibrio. Se puede emplear la correlación de
Wilson.
 Calcule la composición de la fase líquida usando la definición matemática de la constante de
equilibrio, teniendo en cuenta que 𝑦𝑖 = 𝑧𝑖
 Calcule 𝑝𝑑𝐴 usando la composición de la fase gaseosa, 𝑧𝑖 , y Φ𝑖𝐿 usando la composición de la fase
líquida a: 𝑝𝑑𝐴 , 𝑝𝑑𝐴 + ∆𝑝𝑑 , 𝑝𝑑𝐴 − ∆𝑝𝑑 . Donde ∆𝑝𝑑 = 5 𝑜 10𝑝𝑠𝑖
𝑛 𝑛
𝑓𝑖𝑣 𝑓𝑖𝑣
𝑝𝑑 = ෍ 𝐿 𝑓 𝑝𝑑 = ෍ 𝐿 − 𝑝𝑑 = 0
Φ𝑖 Φ𝑖
𝑖=1 𝑖=1

 Evalúe la función 𝑓 𝑝𝑑 y su derivada y su derivada:

𝜕𝑓𝑖𝑣 𝑓𝑖𝑣 𝑝𝑑 + ∆𝑝𝑑 − 𝑓𝑖𝑣 𝑝𝑑 − ∆𝑝𝑑
=
𝜕𝑝𝑑 2∆𝑝𝑑
𝜕𝑓 Φ𝑖𝐿 𝜕𝑓𝑖𝑣 Τ𝜕𝑝𝑑 − 𝑓𝑖𝑣 𝜕Φ𝑖𝐿 Τ𝜕𝑝𝑑
=෍ 2 −1
𝜕𝑝𝑑 Φ𝐿 𝑖 𝜕Φ𝑖𝐿 Φ𝑖𝐿 𝑝𝑑 + ∆𝑝𝑑 − Φ𝑖𝐿 𝑝𝑑 − ∆𝑝𝑑
=
𝜕𝑝𝑑 2∆𝑝𝑑
 𝑓𝑖𝑣 𝑝𝑑 + ∆𝑝𝑑 , fugacidad del componente i @ 𝑝𝑑 + ∆𝑝𝑑

𝑓𝑖𝑣 𝑝𝑑 − ∆𝑝𝑑 , fugacidad del componente i @ 𝑝𝑑 − ∆𝑝𝑑

Φ𝑖𝐿 𝑝𝑑 + ∆𝑝𝑑 , coeficiente de fugacidad del componente i @ 𝑝𝑑 + ∆𝑝𝑑

Φ𝑖𝐿 𝑝𝑑 − ∆𝑝𝑑 , coeficiente de fugacidad del componente i @ 𝑝𝑑 − ∆𝑝𝑑

 Calcule una nueva presión de rocío aplicando la fórmula de Newton-Raphson.

𝑓 𝑝𝑑 Φ𝑖𝐿
𝑝𝑑 = 𝑝𝑑𝐴 − 𝑝𝑑 − 𝑝𝑑𝐴 ≤5 𝐾𝑖 = 𝑣
𝜕𝑓Τ𝜕𝑝𝑑 Φ𝑖

 Chequee convergencia, si no cumple, repita procedimiento calculando un nuevo set de relaciones

de equilibrio.
  Bubble Point Pressure, 𝑷𝒃

𝑛 𝑛
𝑓𝑖𝐿 𝑓𝑖𝐿 𝜕𝑓 Φ𝑖𝑣 𝜕𝑓𝑖𝐿 Τ𝜕𝑝𝑏 − 𝑓𝑖𝐿 𝜕Φ𝑖𝑣 Τ𝜕𝑝𝑏
𝑝𝑏 = ෍ 𝑣 𝑓 𝑝𝑏 = ෍ 𝑣 − 𝑝𝑏 = 0 =෍ 2 −1
Φ𝑖 Φ𝑖 𝜕𝑝𝑏 Φ𝑣
𝑖=1 𝑖=1 𝑖

 Simulation of Laboratory PVT Data by Equations of State

 Constant Volume Depletion
 Asuma que el sistema original con una composición 𝑧𝑖 ocupa un volume inicial de 1ft3
 Usando la composición inicial, calcule el factor de compresibilidad del gas a la presión de rocío.
 Calcule las moles iniciales aplicando la ecuación de gases reales:

(1)𝑝𝑑
𝑛𝑖 =
𝑍𝑑𝑒𝑤 𝑅𝑇

 Reduzca la presión del sistema y calcule constantes de equilibrio, moles y composición de líquido
y gas, factores de compresibilidad…
 Calcule las moles reales de las fases líquida y gaseosa:
𝑛𝐿 𝑎𝑐𝑡𝑢𝑎𝑙 = 𝑛𝑖 𝑛𝐿 𝑛𝑣 𝑎𝑐𝑡𝑢𝑎𝑙 = 𝑛𝑖 𝑛𝑣
 Calcule los volúmenes de líquido y gas aplicando las
siguientes expresiones.
𝑣
𝑛𝐿 𝑎𝑐𝑡𝑢𝑎𝑙 𝑍
𝐿
𝑅𝑇 𝑛𝑣 𝑎𝑐𝑡𝑢𝑎𝑙 𝑍 𝑅𝑇
𝑉𝐿 = 𝑉𝑔 =
𝑝 𝑝

 Remueva el exceso de gas en la celda y numero de

moles removidas.
𝑝 𝑉𝑔𝑝
𝑝,𝑇
𝑉𝑔𝑝 = 𝑉𝑡 − 1 𝑛𝑝 =
𝑝,𝑇
𝑍 𝑣 𝑅𝑇

𝑛𝑖 = 𝑛𝐿 𝑎𝑐𝑡𝑢𝑎𝑙 + 𝑛𝑣 𝑟 𝑛𝑣 𝑟 = 𝑛𝑣 𝑎𝑐𝑡𝑢𝑎𝑙 − 𝑛𝑝

𝑥𝑖 𝑛𝐿 𝑎𝑐𝑡𝑢𝑎𝑙 + 𝑦𝑖 𝑛𝑣 𝑟
𝑧𝑖 =
𝑛𝑖
 Encontrar procedimientos para realización de otras pruebas a partir de la página 415 del libro.