UNIVERSIDAD NACIONAL DEL

CENTRO DEL PER

FACULTAD DE INGENIERA QUMICA
E.A.P. ING. QUMICA DEL GAS NATURAL Y ENERGA

CLCULO DEL PUNTO DE BURBUJA Y DE

ROCO EN SISTEMA MULTICOMPONENTES

Alumnos: Cabezas Snchez, Salustria

Catedra: Ingeniera del Gas Natural II
Docente: Ing. Felix Albert Villavicencio
Fecha de entrega de informe: 25 de mayo del 2015

Huancayo - Per

RESUMEN

El presente informe denominado Clculo de punto de burbuja y punto

de roco en sistemas multicomponentes tiene la finalidad de dar a
conocer los mtodos de clculo y las leyes que los rigen. Estos
sistemas

son

muchos

ms

complejos.

El

comportamiento

termodinmico de una mezcla de hidrocarburos depende de las

condiciones de presin y temperatura en que se encuentre, por lo que
un diagrama P vs T.
El informe se divide en cuatro partes: La primera explica la definicin
de sistemas multicomponentes de hidrocarburos, la segunda sobre el
coeficiente de equilibrio de estos sistemas y los dos ltimos, acerca
de los mtodos para calcular la temperatura y presin en el punto de
burbuja y de roco.
Finalmente espero que esta recopilacin de informacin se de gran
utilidad para el lector.

MARCO TERICO
COMPORTAMIENTO DE FASES DE LOS HIDROCARBUROS
El comportamiento de fases de los hidrocarburos es muy utilizado para mostrar
los efectos de presin y temperatura sobre los estados fsicos de un sistema de
hidrocarburos, adems, el ingeniero puede determinar desde un punto de vista
tcnico, el tipo de fluidos que contiene un yacimiento y posteriormente predecir
su comportamiento en el yacimiento.
Fase: Una porcin fsicamente distinta de materia que es uniforme en su
composicin y propiedades se llama fase. Las fases presentes en la materia son
el slido, lquido y gas; Fig. 1 Como se sabe en la fase slida los cuerpos poseen
volumen definido y forma propia a ciertas condiciones de presin y temperatura,
los lquidos poseen volumen propio pero no forma propia y los gases carecen de
ambos. Los lquidos y gases se consideran como fluidos. Un lquido en la medida
que llene un recipiente adoptar la forma de este, pero retendr su volumen,
mientras que un gas siempre llenar cualquier recipiente que lo confine.

Figura N1 Diagrama de fase P vs T para una sustancia pura

Los hidrocarburos presentan un comportamiento multifsico en altos rangos de

presin y temperatura. Las fases ms importantes para el estudio de

comportamiento de fases que ocurre en las reservas de Petrleo son la fase

lquida y de gas. Las condiciones bajo las cuales estas fases existen son de gran
importancia prctica para la clasificacin de yacimientos, las determinaciones
experimentales

matemticas

de

estas

condiciones

se

expresan

convenientemente en diferentes tipos de diagrama del tipo redondo llamados

diagramas de fase.
SISTEMAS BINARIOS
Cuando se agrega un segundo componente a un sistema, su comportamiento se
hace ms complejo, debido a la introduccin de una nueva variable: la
composicin. El efecto de esta nueva variable puede observarse comparando
los diagramas para uno y dos componentes, para el caso de diagramas P-T. En
un sistema binario o de ms componentes y para una composicin dada, las
lneas de punto de burbujeo y de roco ya no coinciden, sino que forman una
lnea envolvente en el diagrama de fases. Adems para cada composicin
diferente del sistema, se obtiene un diagrama de fases diferentes.

Figura N Diagrama de P-T para diferentes mezclas de metano-propano

SISTEMAS MULTICOMPONENTES
El comportamiento de fases de un sistema de hidrocarburos multicomponente
en la regin de dos fases es muy similar al comportamiento que se presenta en
los sistemas binarios. Pero los sistemas multicomponentes son mucho ms
complejos, con un gran nmero de diferentes componentes hidrocarburos y no
hidrocarburos y en ocasiones slidos, con rangos de presin y temperaturas
diferentes. Esto hace que el comportamiento de las isotermas cambie de una
forma perceptible a una imperceptible.
El comportamiento termodinmico de una mezcla de hidrocarburos depende de
las condiciones de presin y temperatura en que se encuentre, por lo que un
diagrama p vs. T del fluido nos dar una clara idea de qu tipo de fluido se est
manejando, y el comportamiento que podra tener a diferentes condiciones de
presin y temperatura.

I.

EQUILIBRIO DE FASES

Una mezcla de multicomponentes de dos fases se dice que est en equilibrio si

satisface las siguientes condiciones:
1) La temperatura de la fase vapor es igual a la temperatura de la fase
lquida
2) La presin total a travs de la fase de vapor es igual a la presin total a
travs de la fase lquida.
3) La tendencia de cada componente a escapar de la fase liquida a la fase
vapor es exactamente igual a la tendencia de escapar de la fase de vapor
a la fase liquida.
Para poder determinar la composicin de

las fases gaseosa y lquida en

sistemas multicomponentes en la regin de dos fases, se consider que el

sistema, segua un comportamiento de soluciones ideales. En este tipo de
soluciones no se producen fuerzas de atraccin de ninguna especie entre los
constituyentes moleculares y por lo tanto pueden aplicarse la layes de Raoult y
de Dalton.

Ley de Raoult

En una solucin ideal, la presin parcial, Pi, de un componente, i, en la fase de

vapor es igual al producto de la fraccin molar del componente en la fase lquida,
xi, y la presin de vapor, (P)i, del componente i puro, es decir:
= . ( )
En este caso de soluciones ideales, la suma de las presiones parciales de los
componentes, ser la presin total y tambin la presin de vapor del sistema.
= . ( )
Ley de Dalton
En una mezcla de gases perfectos cada gas ejerce una presin parcial igual a
la que se ejercera si estuvieran solo en el volumen ocupado por la mezcla.
Si una mezcla gaseosa contiene n1, n2moles de cada uno de los
componentes 1,2...n, para un componente, i, segn la ley de Dalton, la presin
Pi, ser:
= .

Y la presin total de la mezcla ser,

= = .

Las ecuaciones anteriores se dividen:

=

Pero ni/ni, por definicin, es la fraccin molar del componente, i, en la mezcla,
o sea por definicin, yi. Por lo tanto:
Pi=yi.Pt
Igualando esta ecuacin, con la ecuacin de Raoult

 . ( ) = .
( )
=
=

Si se tiene una mezcla de dos componentes:

= .
= .

= /

= .

= /

Siguiendo con la siguiente ecuacin de la Ley de Raoult se tiene:

= . ( )
= . ( )

= . ( )

Luego se tiene:
. = . ( )

. = . ( )

. ( )

. ( )

Adems: Pt=PA+PB, XA+XB=1; YA+YB=1

Por lo tanto:
. ( ) . ( )
+
=1

. ( ) + . ( ) =

+ = 1

Estas dos ltimas ecuaciones se resuelven y se obtiene las siguientes:

=

( )
( ) ( )

( )
=
( ) ( )

Equilibrio Lquido-Vapor en una solucin ideal

A continuacin se estudiarn las consecuencias de la ley de Raoult en los
sistemas binarios en los cuales ambos componentes son voltiles.

En una solucin binaria, por ejemplo,

, y se tiene:
..(1)

Si se denomina P a la presin total de la solucin, entonces:

..(2)
La anterior expresin relaciona la presin total sobre la mezcla con la fraccin
molar del componente 1 en el lquido, es decir establece que P es una funcin
lineal de x1.
La presin total puede tambin expresarse en funcin de y1, la fraccin molar del
componente 1 en la fase de vapor, as:

(3)
Reemplazando los valores de P1 y de P, se obtiene:

(4)
Resolviendo esta ecuacin en funcin de x1, se halla:

.(5)

Empleando este valor de x1 en la ecuacin

, y despus de

reagrupar, se obtiene:

.(6)
Segn la regla de las fases, resulta que

. Debido a que P es mayor

o igual que 1, entonces F es igual o menor que 3, por lo que se debe especificar
como mximo 3 variables, si una es la temperatura, las otras debern ser (P,x1) o
(P,y1) .
De modo que la presin total, P, slo podra representarse grficamente con
respecto a la fraccin molar del componente 1 en el lquido (x1), la cual describe
la composicin del lquido y de todo el sistema si slo hay lquido presente o con
respecto a la fraccin molar del componente 1 en la fase vapor (y1), que describe
la composicin del vapor y de todo el sistema si slo hay vapor presente.
Surge ahora el problema de cmo representar grficamente los estados lquido
y vapor en un solo diagrama. Para ello es razonable representar de forma
esquemtica, la presin en funcin de x1, la fraccin molar del componente 1 en
todo el sistema (fase lquida + fase vapor).
En la figura 3 se presenta un diagrama en el que se representa P en funcin
de x1 para un sistema cualquiera en el que estn en equilibrio las fases lquido y
vapor. La curva superior (cuando el sistema se comporta idealmente) es una
lnea recta que se denomina curva del lquido, e indica el lmite entre las regiones
de lquido y de la mezcla lquido-vapor y es simplemente una grfica de la presin
de vapor total de la mezcla lquida contra x1. Por encima de esta lnea, el vapor
se condensa formando un lquido.

Figura N2 Diagrama de fase P2 vs T para una sustancia pura

La curva inferior se denomina curva del vapor, es una grfica de P contra y1, e
indica el lmite de fases entre las regiones de lquido-vapor y de vapor. Por
debajo de esta lnea, el lquido no puede existir en equilibrio.
Figura N3 Esquema de un diagrama P vs. x1

Fuente: Ayuso L. (1991). Fisicoqumica. Santaf de Bogot, Unisur.

Considrese un punto cualquiera (b) dentro de la regin lquido + vapor; el
valor x1 correspondiente a b es la fraccin molar del componente 1 en todo el
sistema, lquido + vapor.

Hasta ahora, se ha empleado x1 o y1 para describir el sistema, sin embargo, se

sabe que

, y

, y por tanto tambin se hubiese podido

escoger x2 y y2.
Para conocer la composicin del sistema en el punto b, se procede as: La
interseccin de una lnea horizontal (lnea de unin o de enlace), a presin
constante, con la curva de lquido en (l) da el valor x1 de que describe la
composicin del lquido, mientras que la interseccin con la curva de vapor en
(v) da el valor de y1, que describe la composicin del vapor.
Si x1 es la variable independiente, la interseccin de la lnea vertical en x1 con la
curva de lquido establece el valor deP. A partir de P puede obtenerse de
inmediato el valor de y1. As se define el sistema totalmente, especificando tan
soloP y x1 o P y y1, en la regin de dos fases.

Equilibrio Lquido-Vapor en soluciones reales

El tratamiento antes descrito es para un sistema que tiene comportamiento ideal
o que presenta desviaciones pequeas a la ley de Raoult.
Cuando se observan grandes desviaciones positivas con respecto a la ley de
Raoult, habr un mximo en la curva de presin de vapor total (P) en funcin de
la composicin.
Entonces, se establecer un equilibrio lquido- vapor que se esquematiza en la
figura 17. (Por ejemplo cloruro de etileno-etanol, para el cual se presenta un
mximo a 38% molar de etanol a una temperatura constante de 313.1 K).

De manera similar, habr un mnimo en la curva de presin de vapor total para

las desviaciones negativas de la ley de Raoult. El equilibrio lquido-vapor se
establece como se indica en la figura 18, la cual esquematiza el sistema acetona-

cloroformo. En este caso se produce un mnimo a 58% molar de cloroformo a

308 K.
Figura N4 Esquema del equilibrio binario lquido-vapor Acetona- cloroformo

Fuente: Maron & Prutton. (2002). Fundamentos de Fisicoqumica. Mxico.

Limusa S.A.
Estos puntos mximos o mnimos en las curvas Presin de vapor vs.
Composicin reciben el nombre de azetropos.

Clculo del equilibrio en funcin de volatilidades

El clculo del equilibrio en funcin de volatilidades es especialmente til cuando
las volatilidades son independientes de la temperatura, aunque tambin puede
usarse en aquellos casos en que stas varan con la temperatura.
La volatilidad relativa se define:

(7)

Donde el subndice b identifica al componente base con respecto al que se

calcula la volatilidad, que normalmente es uno de los de la mezcla.
Para un componente i:

.(8)
que puede reordenarse:

..(9)
y realizando un sumatorio para todo i, se obtendra cualquiera de las dos
expresiones siguientes que corresponden a la ecuacin para el clculo de
puntos de burbuja:

(10)
Anlogamente, de la ecuacin (8) se llega a:

(11)
y el sumatorio para todo i permite obtener las siguientes expresiones, que
permiten el clculo del punto de roco en funcin de volatilidades:

.(12)
Por tanto, las ecuaciones del punto de burbuja y de roco en funcin de las
volatilidades relativas se reducen a las ecuaciones (7) y (8), respectivamente. Si
las volatilidades son independientes de la temperatura, el valor de
Kb correspondiente a un conjunto de valores de xi o yi se determina directamente
de las ecuaciones (10) o (12) y, una vez calculado, la temperatura

correspondiente se obtiene fcilmente leyendo en grficos de K b = f(T), o de las

relaciones algebraicas entre ambas magnitudes.
Cuando las volatilidades relativas varan con la temperatura se pueden utilizar
esquemas similares a los descritos en el apartado 2.3.2, pero sustituyendo la
expresin del

i=

ixi

= 1 por la ecuacin (10) y el xi =

i/Ki

= 1 por la

ecuacin (12).

II.

CLCULO DEL PUNTO DE BURBUJA

En el caso del clculo del punto de burbuja (condiciones correspondientes al

caso en que la mezcla se encuentra como lquido saturado, en equilibrio con el
vapor) se especifican la composicin del lquido y la presin (o la temperatura) y
se trata de determinar el valor de la temperatura (o de la presin, P) corresponde
equilibrio, as como la composicin del vapor en equilibrio. La ecuacin que debe
resolverse es la siguiente:

Para el clculo del punto de burbuja existen tres formas de resolver, el primero
es usando un software y los dos ltimos mediante tablas y clculos manuales.
2.1 Algoritmo de clculo para hallar el punto de burbuja (temperatura y
presin)
En los siguientes esquemas se muestra el algoritmo de clculo para resolver un
problema de clculo de temperatura y presin de burbuja:

Figura N5 Algoritmo de clculo de presin en el punto de burbuja

|

Figura N6 Algoritmo de clculo de la temperatura en el punto de burbuja

2.2 Mtodo de Dodge:

Este mtodo se lleva a cabo mediante los grficos del coeficiente de
equilibrio de cada componente. Para ello se siguen los siguientes pasos
1) Determinar el componente clave (Debe ser el ms voltil y el
menos pesado)
2) Se debe suponer una Temperatura
3) Se obtiene o se halla Ki para cada componente a la T supuesta.
4) Se efecta yi=Ki.Xi para cada componente.
5) Se suma yi para calcular la sumatoria, yi=(Ki.Xi)
6) Si yi 1, se calcula en yi modificado en:

7) Determinar Ki* modificado:

, donde i es el componente clave escogido

8) Con el Ki modificado en el grfico de coeficientes de distribucin de

equilibrio K, determinar las nuevas temperaturas.
9) Con la nueva temperatura se halla en grficos los Ki de los otros
componentes.

Observacin: El rango aceptable para hallar la T es: 0.9999yi1,005

2.3 Mtodo de interpolacin-extrapolacin

Para poder aplicar este mtodo se siguen los siguientes pasos:
1) Asumir una T
2) Determinar Ki
3) Efectuar: yi=Ki.Xi para cada componente.
4) Realizar la sumatoria: yi=(Ki.Xi)
5) Si yi 1 entonces se asume como punto en el grfico Si yi vs T.

6) Asumir otra T
7) Repetir los pasos del 2 a 5
8) Se une los dos puntos hallados y luego interpolar o extrapolar la
yi =1 para hallar la T de burbuja.
9) Comprobar para la temperatura de punto de burbuja hallada:
0.9999yi1,005

III.

CLCULO DEL PUNTO DE ROCO

3.1 Algoritmo de clculo para hallar el punto de roco (temperatura y

presin)
En los siguientes esquemas se muestra el algoritmo de clculo para resolver un
problema de clculo de temperatura y presin de roco:

Figura N7 Algoritmo de clculo de presin en el punto de roco

Figura N8 Algoritmo de clculo de temperatura en el punto de roco

3.2 Mtodo de Dodge:
Este mtodo se lleva a cabo mediante los grficos del coeficiente de
equilibrio de cada componente. Para ello se siguen los siguientes pasos
1) Determinar el componente clave (Debe ser el ms voltil y el menos
pesado).
2) Se debe suponer una Temperatura.
3) Se obtiene o se halla Ki para cada componente a la T supuesta y
presin de trabajo.
4) Se efecta xi=yi/ki para cada componente.
5) Se suma xi para calcular la sumatoria, xi=(yi/ki)
6) Si xi 1, se calcula en xi* modificado en:

7) Determinar Ki* modificado:

, donde i es el componente clave escogido

8) En el diagrama K vs P se obtiene la nueva temperatura a Ki

modificado y presin trabajo.
9) Con la nueva temperatura se obtiene los Ki

del resto de

componentes y se repite los pasos 4 al 8.

BIBLIOGRAFA

IQ. (2 de Junio de 2010). MetAprox. Obtenido de

http://iq.ua.es/MetAprox/231_puntos_de_burbuja_y_de_roco.html
Perez, R., & Martinez, M. (2011). Ingeniera del Gas Natural, Cracatersticas y
comportamiento de hidrocarburos. Marcaibo-Venezuela: Edwards Nava.
Torres, J., & Turrubiate, A. (2007). Ecuacin General de Balance de Materia y
una Herramienta Computacional para sus clculos. Mxico: Universidad
Nacional Autnoma de Mxico.
UNAD. (21 de Julio de 2012). DATATECA. Obtenido de
http://datateca.unad.edu.co/contenidos/201604/201604_Modulo_exe/exe
fiscoq28enero/leccion_no_25_regla_de_la_palanca_o_regla_de_la_nivel
acion.html