Propiedades Fisicoquimicas y Termodinamicas Del Gas Natural
Propiedades Fisicoquimicas y Termodinamicas Del Gas Natural
Propiedades Fisicoquimicas y Termodinamicas Del Gas Natural
"#$5!
,+%,$(1*1()&9$)$'%:"$4$'*)&;&
2(+4%1$.*4$'*)&1(#&/*)&.*2"+*#&
!"#$%&'()*+&,(+(-&*./"#%&
".$0(+)$1*1&$.1")2+$*#&1(&)*.2*.1(+&
()'"(#*&1(&$./(.$(+$*&1(&,(2+%#(%)&
(),('$*#$-*'$%.&(.&$./(.$(+$*&1(#&/*)&
3"'*+*4*./*51$'$(43+(&6786&
&
&
&
&
&
&
&
Página
INTRODUC CIÓN
1. C O N C E P T OS B ÁSI C OS «««««««««««««««««
1.1 TERMINOLOGÍA BÁSICA PARA CÁLCULOS DE INGENIERÍA
1.1.1 Elementos y Átomos
1.1.2 Compuestos Químicos y Moléculas
1.1.3 Masa Atómica Relativa
1.1.4 0DVD0ROHFXODU5HODWLYD««««««««««««««««««
1.1.5 Mol
1.1.6 Mezcla
1.1.7 Masa, Fuerza y Peso
1.1.8 Longitud, Área y Volumen ................................................................................... 4
1.1.9 Densidad .............................................................................................................. 5
1.1.10 Presión
1.1.11 Temperatura
1.1.12 Trabajo y Potencia
1.1.13 (QHUJtD&DORUtILFD«««««««««««««««««««««««
1.1.14 Viscosidad
1.1.15 Conversiones Masa ± Volumen
1.2 48Ë0,&$25*È1,&$««««««««««««««««««««
1.2.1 Términos Básicos
1.2.1.1 Teoría Estructural
1.2.1.2 Enlace Químico
1.2.2 1RPHQFODWXUDGH+LGURFDUEXURV«««««««««««««««««
1.3 '(),1,&,21(6«««««««««««««««««««««««
2. PR O PI E D A D ES D E L G AS N A T U R A L ......................................................... 32
6 E Q U I L I B R I O G AS - L Í Q U I D O .................................................................. 247
6.1 SOLUCIONES IDEALES
6.1.1 Ecuación de Raoult .......................................................................................... 248
6.1.2 Ecuación de Dalton
6.1.3 Composiciones y Cantidades de las Fases Gas y Líquido en Equilibrio de
una Solución Ideal
6.1.4 Cálculo de la Presión del Punto de Burbuja de una Solución Líquida Ideal ..... 251
6.1.5 Cálculo de la Presión del Punto de Rocío de una Solución de Gas Ideal ........ 252
INTRODUC CION
La industria del gas en Colombia tiene sus inicios en los años setenta en regiones como
la costa Atlántica y Santander, en aplicaciones para los sectores industrial,
termoeléctrico y doméstico. Los descubrimientos sucesivos de campos de gas en la
Guajira, Huila y Piedemonte Llanero, marcan un nuevo rumbo en el desarrollo de dicha
industria, que empieza a consolidarse con la puesta en marcha del Plan Nacional de
Masificación, cuyo principal objetivo es contribuir a la solución de las necesidades
energéticas del país. La ingeniería del gas en Colombia está en proceso de
fortalecimiento, con miras a responder al nuevo reto de una industria en formación.
La materia que nos ocupa se refiere a la esencia misma del gas, es decir, sus
propiedades fisico-químicas y termodinámicas. Del conocimiento que se tenga en este
aspecto, se deriva el comprender su comportamiento, el entender de los procesos y
transformaciones que se llevan a cabo, así como la seguridad en la toma de decisiones al
diseñar los diferentes procedimientos que sean necesarios para obtener un producto de
la mejor calidad que responda de manera eficiente a las condiciones requeridas.
Este documento esta organizado de la siguiente manera: un primer capítulo que trata
sobre los conceptos básicos sobre el gas natural. Un segundo capítulo presenta las
diferentes maneras de calcular las propiedades el gas natural y los líquidos producidos
en su manejo. En la tercera parte, se muestra las ecuaciones de estado más utilizadas
para el manejo de los hidrocarburos, teniéndose en cuenta posteriormente una sección
para definir la caracterización de las fracciones pesadas. Finalmente, presenta los
conceptos básicos de comportamiento de fase y el equilibrio líquido - vapor utilizando
Ecuaciones de Estado.
1. C O N C E P T OS B ASI C OS
A lo largo del libro utilizaremos terminologías químicas y físicas básicas, para lo cual a
continuación haremos una breve revisión.
1.1 T E R M I N O L O G Í A B ÁSI C A PA R A C Á L C U L OS D E I N G E N I E R Í A .
Toda la materia en el universo está compuesta por elementos los cuales no pueden ser
desglosados o subdivididos en partes más pequeñas de forma ordinaria. Más de 100
materiales han ido encontrados los cuales son clasificados como elementos. Entre estos
están carbono, hidrogeno, azufre, oxigeno, nitrógeno y helio.
El átomo es la unidad básica de cada elemento que puede ser combinado con sí mismo
o con los átomos de otros elementos para formar un compuesto.
Un compuesto es una sustancia de más de un átomo que satisface las siguientes dos
condiciones:
También podemos tener la molécula biatómica, la cual está formada por la combinación
de dos átomos del mismo elemento. Los ejemplos más comunes son Nitrógeno ( ) y
Oxigeno ( ).
CARBONO C 12
NITROGENO N 14
OXIGENO O 16
AZUFRE S 32
El peso molecular relativo de una molécula es la suma de todos los pesos atómicos
relativos de los átomos que componen la molécula. Por ejemplo, el agua tiene una masa
molecular relativa de 18 ( ଶ = 2 + 16 = 18). Una molécula biatómica como el
oxigeno ( ଶ ) tiene un peso molecular relativo de 32.
1.1.5 Mol.
1.1.6 Mezcla.
Una mezcla es una combinación de elementos y compuestos los cuales pueden ser
separados por medios físicos.
Las propiedades de las mezclas son un reflejo de las propiedades de sus componentes.
El gas natural y el crudo son mezclas, y son analizados separando la mezcla en sus
componentes e identificando cada una de sus propiedades.
Una de las mayores diferencias entre el sistema métrico internacional y los sistemas
métricos anteriores es la forma en que la relación fuerza-masa es usada. La tabla 1.3 nos
muestra cinco sistemas que han sido usados.
El peso es un término que debe evitarse para denotar masa. El peso es una medida de la
atracción de la gravedad sobre una masa dada, por lo tanto, varía con la magnitud de la
gravedad, el peso es realmente una fuerza.
La diferenciación entre peso y fuerza ha sido confundida por la definición de una libra
fuerza (lbf) como una libra masa (lbm) sujeta a la aceleración gravitacional de la tierra
de 32.2 ft/ଶ a 45º latitud y a nivel del mar. Estableciendo la constante de
proporcionalidad igual a la aceleración, pero una libra fuerza era numéricamente igual a
una libra masa. En el sistema métrico internacional no existe esta ambigüedad. La
unidad básica de la fuerza es el Newton (N) y la unidad básica de la masa es el
kilogramo (Kg).
Tabla 1.3
(1) (2) (3) (4) (5)
Masa lbm lbm kg g kg
Longitud ft ft m cm m
Tiempo s s s s s
Fuerza lbf ft-lb/ ଶ kgf dina N
a 32.17 ft/ ଶ 1.0 ft/ ଶ 9.81 m/ ଶ 981 cm/ ଶ 1.0 m/ ଶ
୪ୠ୫ି୲ ୪ୠ୫ି୲ ୩ି୫ ିୡ୫ ୩ି୫
gc 32.17 ୪ୠି୲మ 1.0 1.0 ୩ିୱమ 1.0 ୢ୧୬ୟିୱమ 1.0
୮ୢ୪ିୱమ ିୱమ
La unidad básica del volumen es el metro cubico. El litro (L) también lo llamamos
decímetro cubico. El barril API (bbl), el cual es 42 galones U.S., no hace parte del
sistema métrico internacional.
1.1.9 Densidad.
Las unidades más comunes son kg/ଷ o lbm/ f ଷ . La densidad varía con la presión y
con la temperatura. La densidad de referencia o la estándar usualmente está dada a una
presión atmosférica estándar y a una temperatura como 0ºC, 15ºC, 25ºC, 60ºF, 77ºF. En
la mayoría de los cálculos ingenieriles con participación neta de agua se utiliza la
medida estándar de la densidad. Pero para sistemas de transferencia o de gas, se tiene
que calcular la densidad a la determinada T y P que se esté manejando.
1.1.10 Presión.
Presión es la fuerza ejercida por unidad de área. El termino estándar utilizado para la
presión es el pascal (Pa) = 1 N/ଶ . Una unidad pequeña de presión es 1 lbf/ଶ
(psi)=6895 Pa. Una presión de 500 psi seria 3447500 Pa, un numero difícil de manejar y
tratar. Así que, para estos casos el kilopascal (kPa) y el megapascal (MPa) son usados
más comúnmente. En algunas partes el termino bar es usado como nombre especial para
100 kPa.
Cabeza de fluido, una columna vertical de gas o líquido ejercerá una presión sobre la
parte inferior. (VWD ³FDEH]D´ H[SUHVDGD HQ SLHV R PHWURV SXHGH VHU FRQYHUWLGD a
términos de presión por la ecuación
1.1.11 T emperatura.
En el sistema internacional usamos los grados Celsius (ºC) y Kelvin (K). La palabra
³FHQWtJUDGR´HVWiREVROHWD
Como se observa en los últimos capítulos, la energía calorífica debe tener unidades
equivalentes al trabajo y la potencia, para hacer los cálculos que implican la conversión
entre ellos. En el sistema métrico internacional, el Joule es usado también en la energía
calorífica; el Btu es la contraparte tradicional. Hay varias definiciones para Btu, pero la
tabla de valores internacional es usada para todas las conversiones. Para nuestros
propósitos las diferencias son triviales.
1.1.14 V iscosidad.
En todos los cálculos en este libro se usa la viscosidad dinámica. Sin embargo, la
viscosidad dinámica dividida por la densidad es conocida como la viscosidad
cinemática, viscosidad con unidad primaria centistoke (cSt).
ୢୣ୪୪୧୯୳୧ୢ୭
ɀ = (1.3)
ୢୣ୪ୟ୳ୟ
ୢୣ୪ୟୱ ୢୣ୪ୟୱ
ɀୋ = =
ୢୣ୪ୟ୧୰ୣ ୢୣ୪ୟ୧୰ୣ
Resolviendo la ecuación 1.3, la densidad del agua normalmente puede ser tomada como
1000 kg/ ଷ , 1.0 g/cଷ , o 1.0 kg/L; M del aire es tomada como 29.
ଵସଵǤହ
ɀ = (1.4)
ሺଵଷଵǤହା୍͑ሻ
ଵସଵǤହ
API = - 131.5
ஓై
Los pies cúbicos estándar y los metros cúbicos estándar son datos erróneos que
históricamente se han utilizado para la medición y transferencia de gas natural. En
cálculos de ingeniería, es necesario hacer la conversión de este volumen a términos
apropiados de masa, esto puede llegar a ser una rutina, pero también una molestia.
Este volumen estándar fue designado a una presión y temperatura base. En U.S la
temperatura base es de 60ºF y la presión base es de 14.4-15.025 psia, dependiendo de la
organización. En años recientes 14.696 psia (presión atmosférica estándar), ha sido la
más usada.
Ha habido gran controversia con las antiguas unidades inglesas acerca de las unidades
del sistema métrico internacional. Por supuesto, hay una pequeña diferencia con la
precisión en los cálculos cuando se hace la conversión.
Los siguientes son los factores básicos de conversión aplicados a estas condiciones.
S.Ingles: 379.3 f ଷ /lb-mol, 2636 lb-mol/MMscf Densidad del aire = 0.0764 lb/ f ଷ
S. Metrico: 23.96 ଷ /kmol, 41740 kmol/ͳͲ std ଷ Densidad del aire = 1.21 kg/ଷ
El subíndice sc, siempre significa condiciones estándar. Así que ୱୡ = 1.6 ൈ ͳͲ ଷ Ȁ
Significa tasa de flujo de 1.6 millones de metros cúbicos estándar por día.
Para la conversión más común, la tasa de flujo volumétrico por día es convertida a la
tasa de masa por hora.
୫య ୢ ଵ୩ୌమ ɀ
Kmol/h = ቚ ቚ ቚ ቚ ቚ = 41.7 (ଷ /d) (Ȗ /MW)
ୢ ଶସ୦୰ ୫య
1.2 Q U Í M I C A O R G Á N I C A .
Los compuestos derivados de fuentes orgánicas tienen una cosa en común: todos
contienen el elemento carbono. Hoy la mayoría de los compuestos del carbón son
sintetizados, en lugar de ser obtenidos de fuentes animales o vegetales.
Dos de las grandes fuentes de material orgánico son el petróleo y el carbón. Los
compuestos de estas fuentes se emplean para formar compuestos orgánicos más
complicados los cuales son fundamentales para la civilización actual.
Hay otros elementos (tales como el boro, silicio y fósforo) que pueden formar cadenas
de átomos enlazados unos con otros. El carbono es único debido a que no solamente
forma fuertes enlaces carbono-carbono, sino que también estos enlaces permanecen
fuertes cuando los átomos de carbono están enlazados con otros elementos. Los
compuestos del carbono son estables y químicamente no- reactivos.
Las moléculas que tienen estructuras similares exhiben propiedades físicas y químicas
similares. De esta forma, los químicos orgánicos han sido agrupados en familias
OODPDGDV³VHULHVKRPRORJDV´HQODVFXDOHVODVHVWUXFWXUDVPROHFXODUHV\ODVSURSLHGDGHV
son similares.
Recordemos que cada átomo consiste de un núcleo cargado positivamente rodeado por
electrones cargados negativamente ubicados en órbitas concéntricas. Hay un máximo
número de electrones que pueden ser acomodados en cada órbita: dos en la primera
órbita, ocho en la segunda, ocho o dieciocho en la tercera órbita.
La mayor estabilidad se alcanza cuando la órbita exterior esta llena, como en el caso del
helio, que tiene dos electrones en su única órbita. El helio es no-reactivo.
E nlace iónico.
E nlace covalente.
H + H H H
F + F F F
H
H + H + O H O
Los enlaces covalentes son los de interés en nuestro estudio de química orgánica. El
elemento fundamental de la química orgánica estructural es el átomo de carbono
tetravalente.
El enlace que resulta de compartir dos electrones, como el que se presenta en el enlace
carbono - hidrógeno y el enlace carbono - carbono, en compuestos como el etano, es
llamado enlace simple. El átomo de carbono es capaz de formar cuatro enlaces simples.
Un átomo de carbono puede compartir dos o tres electrones con otros átomos diferentes
del carbón. Ejemplo de esto ocurre en los enlaces dobles del dióxido de carbono y el
enlace triple del cianuro de hidrógeno.
O C O H C N
Dióxido de carbono Cianuro de hidrógeno
Normalmente una línea recta que conecta a los símbolos atómicos representa a un
enlace simple; dos líneas representan a un enlace doble ; y tres líneas representan a un
enlace triple. Las fórmulas estructurales de algunos compuestos del carbono se dan a
continuación:
O C O H C N
H
H H
H C H C C H C C H
H H
H
O
C H3 C
OH
Acido acético
C O2 HCN C H4
Dióxido de Carbono Cianuro de hidrógeno Metano
Con base en la estructura, los hidrocarburos son divididos en dos clases principales,
alifáticos y aromáticos. Los hidrocarburos alifáticos posteriormente se dividen en las
siguientes familias: alcanos, alquenos, alquinos y sus análogos cíclicos.
Una familia de químicos orgánicos es conocida como serie homóloga. Los miembros de
una serie homóloga tienen estructuras moleculares similares y tienen propiedades físicas
que difieren de otro miembro de acuerdo al número de átomos de carbono en la
estructura.
%&'()*+,!
Los compuestos de estas familias algunas veces son llamados hidrocarburos saturados
debido a que los átomos de carbono están unidos a tantos átomos de hidrógeno como
sea posible. Estos alcanos también son llamados hidrocarburos parafínicos. Los
Ingenieros de Petróleos normalmente los llaman parafinas. La tabla 1.1 da algunos
ejemplos.
T abla 1.1
Alcanos (CnH2n+2)
Número de átomos
De carbono, n ! Nombre
1 Metano
2 Etano
3 Propano
4 Butano
5 Pentano
6 Hexano
7 Heptano
8 Octano
9 Nonano
10 Decano
20 Eicosano
30 Triacontano
H H H H H H
H C H H C C H H C C C H
H H H H H H
Nótese que los átomos de carbono forman cadenas continuas. Cuando el número de
átomos de carbono incrementa en el compuesto, los átomos de carbono pueden estar
conectados en cadenas continuas o pueden estar conectados como ramales unidos entre
sí. Por ejemplo el butano puede tomar cualquiera de las dos estructuras dadas a
continuación:
n-butano iso-butano
H H H H H C H
H H
H C C C C H H C C C H
H H H H H H H
H H H H H H H H H
H C C C C C H H C C C C H
H H H H H H H H
H C H
C H3 C H2 C H2 C H2 C H3 H
C H3 C H2 C H C H3
n-pentano
C H3
Isopentano
H
H H C H H
H C C C H
H H C H H
H
Neopentano
C H3
C H3 C C H3
C H3
Los puntos de ebullición y fusión de los alcanos son bastante bajos, particularmente en
el caso de los miembros más pequeños de la familia. Esto se debe a que las moléculas
son altamente simétricas, y por lo tanto las atracciones entre las moléculas, conocidas
como fuerzas intermoleculares, son bastante pequeñas. Ya que los procesos de
ebullición y fusión involucran vencer estas fuerzas intermoleculares, los puntos de
fusión y ebullición son más altos para moléculas más grandes debido a que las fuerzas
intermoleculares son más grandes entre moléculas de mayor tamaño.
Todos los hidrocarburos son atacados por el oxígeno a temperaturas elevadas. Si hay
suficiente oxígeno presente, ocurre combustión completa para producir dióxido de
carbono y agua. Esto es particularmente importante ya que los alcanos son los mayores
constituyentes del gas natural, y el principal uso del gas natural radica en la combustión.
A lquenos.
Nomenclatura de los alquenos. Los nombres de los alquenos se forman con los
mismos prefijos que se emplean al nombrar los alcanos. Estos prefijos corresponden al
número de átomos de carbono en el compuesto. El sufijo -eno, indica que el compuesto
pertenece a la familia alqueno. Así, el miembro más simple de la familia alqueno,
C2H2, se debe llamar eteno. El eteno se conoce comúnmente como etileno. El siguiente
miembro de la familia es el propeno, comúnmente conocido como propileno.
H H
H H H
C C H C C C
H H H
H
C H2 C H2 C H3 C H C H2
Los enlaces dobles son menos estables que los enlaces simples y por lo tanto son más
susceptibles a ser atacados por otros químicos. De esta forma los alquenos, debido a sus
dobles enlaces, son más reactivos que los alcanos. La mayoría de las reacciones de los
alquenos involucran la eliminación del doble enlace y la formación de dos enlaces
simples más fuertes.
A lquinos.
Nomenclatura de los alquinos. Las reglas IUPAC (International Union of Pure and
Applied Chemistry) son exactamente las mismas para nombrar alquinos que para
nombrar alquenos, excepto que el sufijo -ino reemplaza al sufijo -eno.
De acuerdo con este sistema, el alquino más simple se debe llamar etino; sin embargo
este compuesto comúnmente es llamado acetileno.
Por razones que aún se desconocen, el triple enlace carbono-carbono es menos activo
que el doble enlace carbono-carbono hacía algunos reactivos y más activo hacía otros
reactivos.
Al igual que los alquenos, la mayoría de las reacciones químicas de los alquinos
involucran la eliminación del triple enlace a favor de un enlace simple o un enlace
doble. Normalmente las reacciones continúan de tal manera que el enlace doble es
eliminado para formar dos enlaces simples.
Serie aromática.
La mayoría de los crudos contienen solamente una traza de aromáticos. Algunos pueden
contener cantidades significativas. Cualquier análisis de crudo o gas natural debería
incluir aromáticos.
C
H C C H
Benceno
H C C H
C
Serie nafténica.
Radicales.
Un radical representa a un grupo de átomos que actúan como una sola unidad en la
formación de muchos compuestos comunes.
El paréntesis indica el grupo radical. El radical alquil normalmente tiene una valencia
de +1. En muchos casos el radical alquil se indica por el símbolo `R´. La fórmula del
metanol es CH3OH; para el etanol es C2H5OH. Ambos se pueden escribir como ROH.
Cuando se utiliza R, no es posible identificar el radical específico. La letra R se emplea
para mostrar las características generales de una reacción.
Radical hidroxilo (O H). Esta combinación ocurre entre muchos compuestos comunes.
El radical OH se combina con el hidrógeno para formar agua (H2O), con sales metálicas
como el sodio, calcio, y magnesio para formar hidróxidos (bases); y con radicales alquil
para formar alcoholes, tales como el metanol, etanol, etc.
(SO 4), (C O 3). Si radicales como estos se combinan con hidrógeno, forman ácidos.
Cuando se combinan con sales metálicas como sodio, calcio y magnesio se forman
sales (las cuales ocurren comúnmente en sistemas de agua). Las escamas formadas en
sistemas de agua son causadas por la precipitación de sales como éstas. Los nombres
comunes para algunos radicales de este tipo son:
SO4 - sulfato
SO3 - sulfito
A lcoholes.
Los alcoholes más comunes son formados por la adición de un solo radical hidroxilo a
un radical alquil. El nombre de los alcoholes termina en -ol, o el nombre del radical
alquil seguido por la palabra alcohol.
H
CH3OH Metanol o metil alcohol
H C OH
H
H H
H H
Tanto el C2H5OH como el CH3OH se pueden escribir como ROH al denotar la reacción
general de un alcohol.
Mercaptanos.
Los compuestos con la fórmula general RSH son conocidos como mercaptanos.
También se les conoce como alcoholes de azufre, ya que la fórmula es la misma si se
reemplaza el átomo de oxígeno en el radical (OH) por un átomo de azufre.
Hay varios otros compuestos de azufre- carbono presentes en crudos agrios. Algunos
son:
Sulfuro de carbonilo
Disulfuro de carbono
tiofeno
El azufre es un elemento muy reactivo que se combina químicamente con muchos otros
elementos y compuestos. Sus compuestos reaccionan con el acero al carbón para
formar sulfuros y óxidos de hierro. Muchos compuestos polimerizan y forman el fango,
tan común en los sistemas de crudos agrios. Este fango es muy corrosivo y debería
removerse por filtración.
Las alcanolaminas son comúnmente usadas para tratar líquidos y gases agrios,
particularmente la monoetanolamina y la dietanolamina. Como sus nombres lo indican,
las alcanolaminas pueden ser consideradas como una combinación de alcohol y
amoníaco.
H H H
Monoetanolamina (MEA) HO C C N
H H H
H H
HO C C
Dietanolamina (DEA)
H H
N H
H H
HO C C
H H
Julio César Pérez Angulo
PROPIEDADES FISICOQUIMICAS Y TERMODINAMICAS DEL GAS NATURAL 21
Notar que la única diferencia en los compuestos anteriores, es la manera cómo los
átomos de hidrógeno, pertenecientes al amoníaco, son reemplazados por el radical
(C2H4OH).
G licoles.
Los glicoles son una familia de químicos, algunas veces conocidos como dioles. Ellos
pueden considerarse como alcoholes complejos, ya que contienen radicales hidroxilos y
radicales alquil. Los glicoles empleados para la deshidratación están basados en el
radical etil.
Al igual que la mayoría de compuestos que contienen grupos hidroxilos los glicoles
reaccionan rápidamente con otros compuestos y elementos. En el DEG y el TEG el
átomo de oxígeno es muy reactivo.
H H
C H2C H2 OH
HO C C OH O
H H C H2C H2 OH
C H2 O C H2C H2 OH
C H2 O C H2C H2 OH
Trietilenglicol (TEG)
1.3 D E F I N I C I O N ES
Sistema. Cuerpo de materia con límites finitos (definidos) que puede considerarse, si se
desea, aislado del medio que lo rodea. Representa el material bajo consideración para un
proceso dado. Ejemplo: hidrocarburos en un yacimiento.
Sistema heterogéneo. Está compuesto por dos o más sistemas homogéneos y cuyas
propiedades intensivas cambian bruscamente en las superficies de contacto. Ejemplo: un
yacimiento por debajo de la presión de saturación; petróleo más agua.
F ase. Cuerpo de material homogéneo diferente en sus propiedades a los demás que le
rodean. Las zonas homogéneas, discutidas en sistemas heterogéneos, se refieren a fases.
Existe sólo una fase gaseosa, pero en líquidos y sólidos se suelen llamar fases a los
estratos diferenciados. Por ejemplo, en separadores se habla de tres fases para referirse
al gas, petróleo y agua.
F=C+2-P
Donde:
C = número de componentes constitutivos del sistema.
P = número de fases.
F = número de grados de libertad o número de variables requeridas para determinar el
Estado de un sistema en equilibrio.
Región de dos fases. Zona encerrada por las curvas de punto de rocío y de burbujeo, en
un diagrama presión-temperatura del sistema, donde el gas y líquido coexisten en
equilibrio.
Estas definiciones de propiedades críticas no son válidas para sistemas con más de un
componente. En el caso de mezclas, el punto crítico es el punto en cual todas las
propiedades en el líquido y el gas son idénticas.
C ricondembárico. Presión máxima a la cual las fases líquido y vapor pueden coexistir.
Para componentes puros, el cricondembárico es igual a la presión crítica del
componente.
Punto triple. Condiciones a las cuales coexisten: sólido, líquido y vapor. Es un punto
único para un componente puro.
F actor volumétrico del gas (gas formation volume factor). Es un factor que
representa el volumen de gas libre, a presión y temperatura del yacimiento, por unidad
volumétrica de gas libre a condiciones normales. Se expresa por Bg y sus unidades
generalmente son barriles en el yacimiento (BY) por pie cúbico de gas a condiciones
normales (PCN).
F actor de merma (shrinkage factor). Es el inverso del factor volumétrico del petróleo,
es decir, barriles normales por barril de petróleo a condiciones de yacimiento.
G as natural. El Gas Natural se define como un gas que se obtiene del subsuelo en
forma natural. Casi siempre contiene una gran cantidad de metano acompañado de
hidrocarburos más pesados como etano, propano, isobutano, butano normal, etc. En su
estado natural a menudo contiene una cantidad considerable de sustancias que no son
hidrocarburos como el nitrógeno, dióxido de carbono y sulfuro de hidrógeno. Además
puede contener trazas de compuestos como el helio, sulfuro de carbonilo y varios
mercaptanos. En su estado natural casi siempre está saturado con agua.
G as residual o gas para la venta. Un gas que tiene la calidad requerida para usarse
como combustible industrial o doméstico.
G as ácido o agrio. Gas que contiene más de 1 gramo de H2S por 100 pies cúbicos
estándar de gas. Casi siempre mucho mayor.
G as dulce. Gas que contiene menos de 1 gramo de H2S por 100 pies cúbicos estándar
de gas.
G as rico. Gas que contiene una cantidad de compuestos más pesados que el etano,
alrededor de 0.7 galones de C3+ por 1000 pies cúbicos estándar de alimento a una torre
absorbedora.
G as pobre. Un gas que contiene muy poca cantidad de propano y de compuestos más
pesados que el propano, o el gas efluente de una unidad de absorción.
G as seco (dry gas o lean gas). Son hidrocarburos en estado gaseoso compuestos casi
exclusivamente por metano (generalmente más del 90 por ciento). Puede provenir
directamente de yacimientos de gas, caso en el cual se le denomina también Gas no-
asociado (nonassociated gas), o sea hidrocarburos gaseosos que ocurren como gas libre
en el yacimiento, o también puede provenir de plantas de gasolina natural, donde el gas
húmedo (condensado de gas) ha sido despojado (stripped) de sus productos más pesados
en forma líquida.
G as asociado (associated gas). Hidrocarburos gaseosos que ocurren como gas libre en
un yacimiento a condiciones iniciales, en contacto con petróleo crudo comercialmente
explotable.
Punto de fusión. Temperatura a la cual una sustancia pasa de la fase sólida a la fase
liquida.
2. PR O PI E D A D ES D E L G AS N A T U R A L
Un gas es definido como un fluido homogéneo de baja densidad y viscosidad,el cual no
tiene un volumen definido pero que se expande hasta ocupar todo el recipiente donde
este se encuentre .Las propiedades de presión-volumen-temperatura (PVT) y otras
propiedades físicas y químicas de los gases son escenciales para resolver problemas
ingenieriles que se tienen en yacimientos de gas natural. Estas propiedades se deben
obtener directamente de pruebas de laboratorio o por predicciones de la composición
química conocida del gas. En este último caso los cálculos deben basarse en las
propiedades físicas de cada componente del gas.
2.1 I M PO R T A N C I A D E L C A L C U L O D E PR O PI E D A D ES
El ingeniero estructural no puede diseñar puentes sin conocer las densidades del acero y
del concreto. De una manera similar los científicos y los ingenieros necesitan
frecuentemente de las propiedades de gases y líquidos. En particular, los ingenieros
químicos, de petróleos o de procesos encuentran que el conocimiento de las propiedades
físicas es esencial para el diseño de muchas clases de equipo industrial. Aun el físico
teórico tiene que comparar la teoría con cantidades medidas.
2.2 PR O PI E D A D ES F ISI C O Q U Í M I C AS D E L G AS N A T U R A L
Las propiedades fisicoquímicas del gas natural pueden ser obtenidas directamente ya
sea por mediciones en laboratorio o por predicción a partir de la composición química
conocida del gas. En este caso, los cálculos están basados en las propiedades físicas de
componentes individuales del gas y en las leyes físicas, frecuentemente conocidas como
leyes de mezclas, las cuales relacionan las propiedades de los componentes a la de la
mezcla de gas.
Un gas ideal es un modelo útil porque se describe con una ecuación sencilla que con
frecuencia se aplica a gases reales con una buena aproximación.
Evidencia experimental.
Teoría cinética.
En este texto nos ocuparemos del desarrollo de la ecuación de estado de un gas ideal a
partir de la evidencia experimental.
E cuación de Boyle.
1
Vv ó PV cons tan te 2.1
P
E cuación de C harles.
V
V vT ó cons tan te 2.2
T
L ey de A vogadro.
Las ecuaciones de Boyle, Charles y Avogadro se pueden combinar para dar la ecuación
de estado de un gas ideal.
El proceso se inicia con una masa de gas dada con volumen V1, a presión P1 y
temperatura T 1, y termina con volumen V2, a presión P2 y temperatura T 2.
T 1 = constante P2 = constante
P1V1
P1V1 P2V ó V 2.4
P2
La ecuación de Charles se aplica al cambio en volumen del gas durante el segundo paso
ya que la presión y la masa del gas se mantienen constantes. Por lo tanto,
V V2 V2 T1
ó V 2.5
T1 T2 T2
P1V1 V2 T1
2.6
P2 T2
P1V1 P2V2
2.7
T1 T2
Así, para una masa dada de gas, PV/T es una constante que designaremos con el
símbolo de R cuando la cantidad de gas sea igual a un peso molecular,
PV M
R 2.8
T
Hasta este punto no sabemos si R es la misma para todos los gases ideales. De modo
que necesitaremos considerar la ley de Avogadro, la cual establece que un peso
molecular de un gas ideal ocupa el mismo volumen que el peso molecular de otro gas
ideal a la misma presión y temperatura. Esto es:
V MA V MB 2.9
PV MA RAT
RA ó V MA 2.10
T P
PV MB RB T
RB ó V MB 2.11
T P
R AT RB T
ó RA RB 2.12
P P
De esta forma la constante R es la misma para todos los gases ideales y es llamada
"Constante Universal del Gas".
Por lo tanto, la ecuación de estado para un peso molecular de cualquier gas ideal es,
PV M RT 2.13
Donde V M es el volumen de un peso molecular del gas, el volumen molar. Para n mole
de gas ideal esta ecuación se convierte en:
PV nRT 2.14
m RT
PV RT ó PX 2.15
M M
Esta expresión es conocida con varios nombres tal como la "Ley del Gas Ideal", "Ley
General del Gas", o "Ley del Gas Perfecto". Aquí la llamaremos "Ecuación de Estado
de un Gas Ideal o la Ecuación del Gas Ideal".
El valor numérico de la constante R depende de las unidades usadas para expresar las
variables. La tabla 2.1 da los valores numéricos de R para varios sistemas de unidades.
Las unidades de campo usadas en la industria del gas y petróleo para la ecuación de gas
ideal son:
Solución:
m PVM
PV nRT PV RT m
M RT
¸ 3.20 ft
lb ·
1,000 psia §¨16.04 3
© lb m ol ¹
m
§ psia ft ·
3
¨¨10.732 ¸ 528 q R
© lb mol q R ¸¹
m 9.1 lb
Ya que la densidad está definida como la masa por unidad de volumen de sustancia, la
ecuación del gas ideal puede ser resuelta para la densidad:
m PM
Ug 2.16
V RT
Al tratar con gases a muy baja presión, la ecuación del gas ideal es una herramienta
conveniente y satisfactoria. Para el cálculo de propiedades físicas de gases naturales a
elevadas presiones, el uso de la ecuación de estado del gas ideal puede conducir a
errores mayores de los 500%, comparados con errores de 2 a 3% a presión atmosférica.
Solución :
lb ·
20 psia §¨ 44.097 ¸
© lb mol ¹
Ug 0.1468 lb ft 3
§ psia ft 3
·
¨¨ 10.732 ¸¸ 100 460
© lb m ol q R ¹
La presión total ejercida por una mezcla de gases es igual a la suma de las presiones
ejercida por sus componentes. La presión ejercida por cada uno de los componentes
gaseosos es conocida como su presión parcial.
Dalton postuló que la presión parcial ejercida por cada componente en una mezcla de
gases es igual a la presión que ejercería cada componente si estuviera sólo en el
volumen ocupado por la mezcla de gases. Esto es válido cuando la mezcla y cada
componente de la mezcla actúan como gases ideales.
La presión parcial ejercida por cada componente de la mezcla de gases se puede calcular
usando la ecuación de gas ideal. Considerar una mezcla que contiene nA moles del
componente A, nB moles del componente B, n C moles del componente C , y así
RT RT RT
PA nA , PB nB , PC nC , 2.17
V V V
P PA PB PC 2.18
De esta forma:
RT RT RT
P nA nB nC 2.19
V V V
RT RT
P ¦nj n 2.20
V j V
Pj nj nj
yj 2.21
P ¦nj n
j
Donde yj está definida como la fracción molar del componente j en la mezcla de gases.
Por lo tanto, la presión parcial de un componente de una mezcla de gases ideales es el
producto de su fracción molar y la presión total.
Pj yjP 2.22
Recuerde que esto es solamente válido para mezclas ideales de gases ideales.
E jemplo 2.3. Calcular la presión parcial ejercida por el metano en el siguiente gas
cuando el gas está a una presión de 750 psia. Asumir que el gas es una mezcla de gases
ideales.
Componente Composición;
F racción molar
Metano 0.85
Etano 0.10
Propano 0.05
1.00
Solución:
PC 1 P yC1
Amagat postuló que el volumen total ocupado por una mezcla de gases es igual a la
suma de los volúmenes que los componentes puros ocuparían a la misma presión y
temperatura. Esto se conoce como la ley de los volúmenes aditivos. La ley de Amagat
de volúmenes parciales es análoga a la ley de Dalton de presiones parciales.
Los volúmenes ocupados por los componentes individuales son conocidos como
volúmenes parciales. Esta ecuación es correcta si la mezcla y cada uno de los
componentes obedecen a la ecuación de gas ideal.
RT RT RT
VA nA , VB nB , VC nC , 2.23
P P P
V V A V B VC 2.24
Así,
RT RT RT
V nA nB nC 2.25
P P P
RT RT
V ¦nj n 2.26
P j P
La razón del volumen parcial del componente j, al volumen total de la mezcla es:
RT
Vj nj nj
P yj 2.27
V RT n
n
P
Vj y jV
Las propiedades básicas de gases ideales se expresan comúnmente en términos del peso
molecular aparente, volumen estándar, densidad, volumen específico y gravedad
específica.
Puesto que una mezcla de gases está compuesta de moléculas de varios tamaños, decir
que una mezcla de gases tiene un peso molecular no es estrictamente correcto. Sin
embargo, una mezcla de gases se comporta como si tuviera un peso molecular definido.
Este peso molecular es conocido como peso molecular aparente:
Ma ¦ yjM j 2.28
j
E jemplo 2.4. El aire seco es una mezcla de gases que consiste esencialmente de
nitrógeno, oxigeno, y pequeñas cantidades de otros gases. Calcular el peso molecular
aparente del aire dada su composición aproximada.
Solución:
Ma y N 2 M N 2 yO 2 M O 2 y A M A
Volumen estándar.
Dada la ecuación de gas ideal: PV = nRT y resolviendo para volumen estándar tenemos:
1 R TS C
VS C 2.29
PS C
VS C 379.6 S C F lb mol
Donde,
Densidad.
PMa
Ug 2.30
RT
Donde,
Ug = densidad de la mezcla de gases, ft3/lb,
M a = peso molecular aparente.
Volumen específico.
El volumen específico está definido como el volumen ocupado por una unidad de masa
de gas. Para un gas ideal, esta propiedad se puede calcular aplicando la ecuación:
§ m · RT
PV nRT V ¨¨ ¸¸
©Ma ¹ P
V RT 1
X 2.31
m P Ma Ug
Donde,
X = volumen específico, ft3/lb,
Ug = densidad del gas, lb/ft3.
La gravedad específica de un gas está definida como la razón entre la densidad del gas y
la densidad del aire seco, ambas medidas a la misma presión y temperatura.
Simbólicamente,
Ug
Jg 2.32
U a ire
Suponiendo que el comportamiento del gas y del aire puede ser representado por la
ecuación de un gas ideal, la gravedad específica puede estar dada como:
PMg
RT Mg Mg
Jg 2.33
P M a ire M a ire 29
RT
Donde, M aire es el peso molecular aparente del aire y Mg es el peso molecular del gas.
Si el gas es una mezcla, esta ecuación se convierte en:
Ma Ma
Jg 2.34
M a ire 29
Donde M a es el peso molecular aparente de la mezcla de gases. Note que esta ecuación
es estrictamente válida si el gas y el aire se comportan como gases ideales.
Componente Composición
F racción molar
Metano 0.85
Etano 0.09
Propano 0.04
n-Butano 0.02
1.00
Solución:
Ma 19.26
Jg 0.664
29 29
ni ni
yi 2.35
n ¦ ni
i
Donde:
mi mi
wi 2.36
m ¦ mi
i
Donde:
Vi Vi
vi 2.37
V ¦ Vi
i
Donde:
n i RT
PV i n i RT Vi
P
n i RT n i RT
Vi P P ni
vi
¦ Vi n RT
¦i i P
RT
¦ ni ¦ ni
i
P i i
vi yi
Paso 1: Asumir que el número total de moles de la fase gaseosa es uno, n =1.
Paso 2: A partir de la ecuación 2.35 es aparente que:
ni = yi
Paso 3: Debido a que el número de moles de un componente es igual al peso del
componente dividido por el peso molecular del componente tenemos:
mi mi
ni yi mi yi M i
Mi Mi
y
m ¦ mi ¦ yi M i
i i
yi M i
wi
¦ yi M i
wi M i
yi
¦ wi M i
Solución:
Componente yi Mi m i = yi M i wi = m i/m
C1 0.65 16.04 10.4260 0.3824
C2 0.10 30.07 3.0070 0.1103
C3 0.10 44.10 4.4100 0.1618
C4 0.10 58.12 5.8120 0.2132
C5 0.05 72.15 3.6075 0.1323
m = 27.2625
Solución:
Los investigadores han propuesto centenares de ecuaciones de estado para gases reales,
las cuales podemos reunir en dos grandes grupos: ecuaciones de compresibilidad y
ecuaciones de estado analíticas. De estas últimas hablaremos con detalle en el
capítulo 3.
ZRT PM
PV ZnRT , PV M ZRT , PQ , Ug 2.38
M ZRT
V rea l
Z 2.39
Videa l
El factor Z no es constante. Varía con los cambios de composición del gas, presión
y temperatura. Para un fluido determinado, este factor tiene la forma mostrada en
la figura 2.1
Los factores Z para varios hidrocarburos se muestran en las figuras 2.2, 2.3 y 2.4
E jemplo 2.8. Calcular la masa de gas metano contenida a 1000 psia y 68°F en un
cilindro con volumen de 3.20 ft3. No asuma que el metano es un gas ideal.
Solución:
PMV
m
ZRT
¸ 3.20 ft
lb ·
1,000 psia §¨16.04 3
© lb mol ¹
m
§ psia ft 3 ·
0.89 ¨10.732
¨ ¸¸ 528q R
© lbm ol q R ¹
m 10.2 lb
Nota: la suposición que se hizo en el ejemplo 2.1 de que el metano actúa como un gas
ideal resultó en un error en la masa de aproximadamente 11%.
La verificación de que esta similitud es cierta para casi todos los gases reales llevó al
desarrollo de la ley de los estados correspondientes y la definición de los términos de
temperatura reducida y presión reducida. La temperatura reducida y la presión reducida
se definen como:
T P
Tr y Pr
Tc Pc
Donde,
Pr = presión reducida.
Pc = presión crítica.
Tr = temperatura reducida.
Tc = temperatura crítica.
La ley de los estados correspondientes dice que todos los gases puros tienen el mismo
factor Z a los mismos valores de presión reducida y temperatura reducida.
La figura 2.5 da una prueba de esta teoría para datos de compresibilidad del etano,
propano, n-pentano y n-hexano. La desviación entre las líneas a temperaturas reducidas
constantes puede ser debida a errores experimentales y a las inexactitudes de la teoría.
a. Metano.
b. Etano.
c. Propano.
Solución:
P T
Pr y Tr
Pc Tc
Entonces,
El ejemplo anterior muestra que a valores iguales de Tr y Pr los factores Z para las tres
sustancias son muy similares, indicando el claro poder del principio de los estados
correspondientes.
PPC ¦ y i PC i 2.40
i
T PC ¦ y i TC i 2.41
i
Donde,
Las ecuaciones 2.40 y 2.41 son conocidas como reglas de mezclas de Kay.
Los estados reducidos para mezclas gaseosas son conocidos como temperatura y presión
pseudo-reducidas, y se expresan por las siguientes relaciones:
P
Ppr 2.42
PPC
T
T pr 2.43
T PC
Donde,
Los estudios de factores de compresibilidad para gases naturales han mostrado que éstos
pueden ser generalizados con suficiente exactitud para la mayoría de los propósitos de
ingeniería introduciendo el concepto de presiones pseudo-reducidas y temperatura
pseudo-reducidas.
§ m ·
PV Z ¨¨ ¸¸ RT
©Ma ¹
Entonces,
V Z R T
X 2.44
m PMa
1 PMa
Ug 2.45
X Z R T
Donde,
Componente yi
C1 0.75
C2 0.07
C3 0.05
n - C4 0.04
n - C5 0.04
C6 0.03
C7 0.02
Asumiendo un comportamiento de gas real, calcular la densidad del gas a 1000 psia y
100°F.
Solución:
Componente yi Mi yi M i T Ci (°R) yi T Ci P Ci yi P Ci
C1 0.75 16.04 12.030 343.5 257.6 673 504.7
C2 0.07 30.07 2.105 550.1 38.5 708 49.6
C3 0.05 44.10 2.205 666.2 33.3 617 30.9
NC4 0.04 58.12 2.325 765.6 30.6 551 22.0
NC5 0.04 72.15 2.886 847.0 33.9 485 19.4
C6 0.03 86.18 2.585 914.6 27.4 434 13.0
C7 0.02 100.21 2.004 972.8 19.5 397 7.9
M a = 26.14 Tpc = 440.8 Ppc = 647.5
P 1000
Ppr 1.54
Ppc 647.5
T pr
T 100 460 12.7
T pc 440.8
Ug
PMa 1,000 26.14 6.0 lb ft 3
Z R T 0.725 10.73 560
La gráfica del factor de compresibilidad, mostrada en la figura 2.6, es aplicable a la
mayoría de los gases encontrados en los yacimientos de petróleo y suministra una
predicción satisfactoria para todos los cálculos de ingeniería.
En los casos donde la composición del gas natural no está disponible, las propiedades
pseudocríticas Ppc y T pc, se pueden predecir solamente a partir de la gravedad específica
del gas. Brown y colaboradores (1948) presentaron un método gráfico para una
aproximación conveniente de la presión y temperatura pseudocrítica de gases cuando
sólo se dispone de la gravedad específica. La correlación es presentada en la figura 2.9.
Standing (1977) expresó esta correlación gráfica en una forma matemática:
Donde,
T pc = temperatura pseudo-crítica, °R
Ppc = presión pseudo-crítica, psia.
Jg = gravedad específica de mezcla de gases.
E jemplo 2.11. Usando los datos dados en el ejemplo 2.10, calcular nuevamente la
densidad del gas estimando las propiedades pseudocríticas a partir de las ecuaciones
2.46 y 2.47.
Solución:
Ma 26.14
J 0.903
29 29
560
T pr 1.24
451.28
1,000
Ppr 1.52
659.97
Z 0.69
Ug
PMa 1,000 26.14 6.3 lb ft 3
Z R T 0.69 1073 560
Errores en los cálculos del factor de compresibilidad mayores al 10% pueden ocurrir en
concentraciones más altas de componentes no - hidrocarburos en mezclas de gases.
Los gases naturales que contienen H2S y/o CO2 frecuentemente exhiben diferente
comportamiento del factor de compresibilidad en comparación con los gases dulces.
Wichert y Aziz (1972) desarrollaron un simple procedimiento de cálculo que tiene en
cuenta estas diferencias. Este método permite usar la gráfica de Standing- Katz, figura
2.6, empleando un factor de ajuste a la temperatura pseudo-crítica, el cual es una
función de la concentración de CO2 y H2S en el gas agrio. Este factor de corrección
luego es usado para ajustar la presión y temperatura pseudo-crítica de acuerdo a las
siguientes expresiones:
T pc ' T pc 2.50
'
' Ppc T pc
Ppc 2.51
T pc B 1 B
Donde,
120 A 0.9 A1.6 15 B 0.5 B 4.0 2.52
A y H 2S y C O2
Los pasos de cálculo para incorporar el factor de ajuste en los cálculos del factor Z se
resumen a continuación:
Paso 3: Ajustar la Ppc y T pc calculada en el paso 1 aplicando las ecuaciones 2.50 y 2.51.
Componente yi
CO2 0.10
H2S 0.20
N2 0.05
C1 0.63
C2 0.02
Solución:
Componente yi Mi yi M i P ci yi P ci T ci yi T ci
CO2 0.10 44.01 4.401 1071 107.1 547.57 54.757
H2S 0.20 34.08 6.816 1306 261.2 672.37 134.474
N2 0.05 28.01 1.401 493 24.65 227.29 11.3645
C1 0.63 16.04 10.105 667.8 420.714 343.06 216.128
C2 0.02 30.07 0.601 707.8 14.156 549.78 11.00
M a = 23.324 827.82 427.72
23.324
Jg 0.8054
28.96
Ppc 827.82
T pc 427.72
1,000
Ppr 1.208
827.82
570
T pr 1.333
427.72
Ug
PMa 1,000 23.324 4.651 lb ft 3
Z R T 0.82 10.73 570
B YH 2S 0. 2
29.86
Ppc'
827.82 397.86 727.07 psia
427.72 0.21 0.2 29.86
1,000
Ppr 1.375
727.07
570
T pr 1.433
397.86
Ug
1,000 23.324 4.56 lb ft 3
0.837 10.73 570
Donde,
E jemplo 2.13. Calcular de nuevo la densidad del gas del ejemplo 2.12 a partir de su
gravedad específica.
Solución:
Debido a que la gravedad específica del gas es 0.8054, calcular T pc y Ppc aplicando
las ecuaciones 2.46 y 2.47.
Ajustar las Ppc y T pc calculadas en el paso anterior, empleando las ecuaciones 2.53 y
2.54.
Ppc' 650.9 4400.10 6000.20 1700.05 763.9 psia
1000 570
Ppr 1.31 , T pr 1.33
763.9 427.15
Z 0.808
Ug
PMa 1,000 23.324 4.72 lb ft 3
Z R T 0.808 10.73 570
Sutton propuso que esta desviación puede ser minimizada utilizando las reglas de
mezcla desarrolladas por Stewart (1959), junto con factores de ajuste empíricos
relacionados con la presencia de la formación de Heptano - plus en la mezcla de gases.
2
§ Tc · º 2 ª § Tc · º
0.5
1ª
J «¦ y i ¨ ¸ » «¦ y i ¨ ¸ » 2.55
3 ¬ i © Pc ¹ i ¼ 3 «¬ i © Pc ¹ i »¼
ª § Tc · º
K ¦ « y i ¨ 0.5 ¸ » 2.56
i ¬ © Pc ¹ i ¼
Donde,
2
2 ª § Tc · º
0.5
1 ª § Tc ·º
FJ y¨ ¸ « y¨ ¸ » 2.57
3 «¬ © Pc ¹»¼ C 7 3 «¬ © Pc ¹ »¼
C7
EK
§ Tc ·
>
¨ 0.5 ¸ 0.3129 y C7 4.8156 y C7 27.3751y C7
2 3
@ 2.59
© Pc ¹ C7
Donde,
J' J EJ 2.60
K' K EK 2.61
Donde,
T PC
K '2
2.62
J'
T Pc
PPC 2.63
J'
Las reglas de mezcla propuestas por Sutton para calcular las propiedades pseudo-
críticas para gases de yacimientos con alto peso molecular (Jg > 1.25) deberían mejorar
la exactitud del factor de compresibilidad.
x Papay.
x Hall ± Yarborough
x Dranchuk y Abu - Kassem
x Dranchuk ± Purvis - Robinson.
x Hankinson ± Thomas - Phillips.
Método de Papay.
Ppr ª § Ppr ·º
Z 1 «0.36748758 0.04188423¨¨ ¸»
¸»
2.64
T pr «¬ © T pr ¹¼
Donde,
E jemplo 2.14. Usando los datos dados en el ejemplo 2.12, evaluar nuevamente el
factor de compresibilidad usando la corrección de Papay, y calcular la densidad del gas.
Solución:
A partir del ejemplo 2.12, Ppr =1.375 y T pr =1.433, resolver la ecuación 2.64 para el
factor Z.
1.375 ª § 1.375 ·º
Z 1 «0.36748758 0.04188423¨ ¸» 0.687
1.433 ¬ © 1.433 ¹¼
Ug
1,000 23.324 5.56 lb ft 3
0.687 10.73 570
Z
ª 0.06125 Ppr t º
«
Y
>
» Exp 1.21 t
2
@ 2.65
¬ ¼
Donde,
>
F Y 0.06125 Ppr t Exp 1.21 t
2
@
Y Y 2 Y 3 Y 4
1 Y 3
14.76 t 9.76 t 2 4.58 t 3 Y 2
90.7 t 242.2 t 2 42.4 t 3 Y 2.182.82t 0 2.66
Esta ecuación no - lineal puede ser convenientemente resuelta para la densidad reducida
Y, empleando la técnica de iteración de Newton Raphson. El procedimiento de cálculo
para resolver la ecuación 2.66 se resume en los siguientes pasos:
YK >
0.06125 Ppr t Exp 1.21 t
2
@
Paso 2: Sustituir este valor inicial en la ecuación 2.66 y evaluar la función no - lineal.
A menos que se halla seleccionado un valor correcto de Y, la ecuación 2.66
tendrá un valor F(Y K ), diferente de cero.
F Y K
Y K 1 YK
F ' Y K
1 4Y 4Y 2 4Y 3 Y 4
F ' Y
29.52t 19.52t 2 9.16t 3 Y
1 Y 4
2.67
2.18 2.82t 90.7t 242.2t 2 42.2t 3 Y 1.18 2.82t
Paso 4: Los pasos 1 al 3 se repiten n veces, hasta que el error (Yn-Yn-1) llegue a ser
más pequeño que una tolerancia dada, por ejemplo, 10-12.
Paso 5: El valor corregido de Y es usado luego para evaluar la ecuación 2.65 para el
factor de compresibilidad Z .
E jemplo 2.15. Calcular el factor Z para el siguiente gas a 3,000 psia y 150°F.
T (°F) = 150
P (psia) = 3,000
Tc (°R) = 491.3837
Pc (psia) = 816.7643
Fp (°R) = 31.1365
Tc' (°R) = 460.2472
Pc' (psia) = 756.2226
Z = 0.6695
ª A2 A A A º ª A A º ªA A º
Z « A1 3 3 4 4 5 5 » U r « A6 7 8 2 » U r2 A9 « 7 8 2 » U r5
¬ Tpr Tpr Tpr Tpr ¼ ¬ Tpr Tpr ¼ ¬ Tpr Tpr ¼
U r2
A10 1 A11 U r2
Tpr 3
Exp A11 U r2 1 2.68
Donde Ur = densidad reducida del gas y está definida por la siguiente expresión:
0.27 Ppr
Ur 2.69
Z Tpr
0.27 Ppr
Ur 2.70
Tpr
ª A A º ª A º §A A ·
Z 1 « A1 2 3 3 » U r « A4 5 » U 2 ¨¨ 5 6 ¸¸ U r5
¬ Tpr Tpr ¼ ¬ Tpr ¼ © Tpr ¹
2.71
ª A º
« 7 3 U r2 1 A8 U r2 Exp A8 U r2 »
¬ Tpr ¼
donde Ur está definido por la ecuación 2.69 y los coeficientes del A1 al A8 tienen los
siguientes valores:
A1 = 0.31506237 A5 = -0.61232032
A2 = -1.0467099 A6 = -0.10488813
A3 = -0.57832729 A7 = 0.68157001
A4 = 0.53530771 A8 = 0.68446549
1 ª A º § Ppr · § Ppr 2 ·
1 « A4 Tpr A2 6 2 » ¨¨ 2 ¸¸ A3Tpr A1 ¨¨ 3 ¸
3 ¸
Z ¬ Tpr ¼ © Z Tpr 2 ¹ © Z Tpr ¹
Tabla 2.2
Coeficientes para la ecuación 2.72
1 § wV ·
Cg ¨ ¸ 2.73
V © wP ¹ T
Donde:
C g = Compresibilidad isotérmica del gas, psi-1
nRTZ
V
P
§ wV · ª 1 § wZ · Z º
¨ ¸ nRT « ¨ ¸ 2 »
© wP ¹ T ¬ P © wP ¹ P ¼
1 1 § wZ ·
Cg ¨ ¸ 2.74
P Z © wP ¹ T
§ wZ ·
Para un gas ideal, Z =1 y ¨ ¸ 0 , por lo tanto,
© wP ¹ T
1
Cg 2.75
P
1 1 § wZ ·¸
Cg ¨
Ppr Ppc Z ¨© w Ppr Ppc ¸¹ T
pr
1 1 ª wZ º
C g Ppc Cr « » 2.76
Ppr Z «¬ wPpr »¼ T
pr
Cr C g Ppc 2.77
Donde:
C r = Compresibilidad isotérmica pseudo-reducida.
C g = Compresibilidad isotérmica del gas, en psia-1.
Ppc = Presión pseudo-reducida, psia.
§ wZ ·
Los valores de ¨ ¸
¨ wP ¸ pueden ser calculados a partir de la pendiente de la isoterma T pr
© pr ¹
sobre la gráfica del factor Z de Standing y Katz en el factor Z de interés.
Componente yi
C1 0.75
C2 0.07
C3 0.05
n-C4 0.04
n-C5 0.04
C6 0.03
C7 0.02
Solución:
1
Cg 1,000 *10 6 psia 1
1,000
1000
Ppr 1.54
647.5
560
T pr 1.27
440.8
§ wZ ·
¨ ¸ 0.1678
¨ wP ¸
© pr ¹ T pr
1 § 1 ·
Cr ¨ ¸ 0.1678 0.8808
1.54 © 0.725 ¹
Cr 0.8808
Cg 1,361.1 *10 6 psia 1
Ppc 647.5
E jemplo 2.17. Empleando las gráficas generalizadas de Trube, repetir el ejemplo 2.16
Solución:
C r = 0.90
0.90
Cg 1390 *10 6 psia 1
647.5
Mattar, Brar y Aziz (1975) presentaron una técnica analítica para calcular la
compresibilidad isotérmica del gas. Los autores expresaron C r como una función de
§ wZ · ª wZ º
¨¨ ¸¸ En lugar de utilizar « » .
© wU r ¹ T pr ¬« wPpr ¼» T pr
ª º
ª wZ º 0.27 « wZ wU r T pr »
« » « » 2.78
«¬ wPpr »¼ ZT pr « U r
«¬
1 wZ wU r T pr »
»¼
Z
ª º
1 0.27 « wZ wU r T pr »
Cr 2 « » 2.79
Ppr Z T´ pr « U
1 r wZ wU r T pr »
«¬ Z »¼
Donde:
Ur = densidad psuedo-reducida del gas.
§ · § ·
¨ A1 A2 A3 ¸ 2 ¨ A4 A5 ¸ U r 5 A5 A6 U r
ª wZ º 4
« » ¨
¬ wU r ¼ T pr
¨
© T pr T pr3 ¸¹ © T pr ¸¹ T pr
2.80
2A U
73 r 1 A8 U r2 A8 U r4 EXP A8 U r2
T pr
2
Introducir
P , T , yi
Calcular
Ppc y Tpc
El factor volumétrico de formación del gas es usado para relacionar el volumen de gas,
medido a condiciones estándar (60 °F y 14.7 psia). Esta propiedad del gas es entonces
definida como el volumen actual ocupado por cierta cantidad de gas a una presión y
temperatura determinada, dividido por el volumen ocupado por la misma cantidad de
gas a condiciones estándar.
V P ,T
Bg
V sc
Donde:
Bg = factor volumétrico de formación del gas, ft3/scf
V P,T = volumen del gas a presión P y temperatura T , ft3
Vsc = volumen del gas a condiciones estándar, scf.
Aplicando la ecuación de estado del gas natural (ecuación 2.38) a la anterior relación
tenemos:
ZnRT
P PS C ZT
Bg
Z S C nRT S C TS C P
PS C
Donde:
Z sc = factor Z a condiciones estándar = 1.0
Psc , T sc= presión y temperatura estándar, 14.7 psi y 60 °F
ZT
Bg 0.02827 , ft 3 scf 2.81
P
En otras unidades de campo, el factor volumétrico de formación del gas puede ser
expresado en bbl/scf, para dar:
ZT
Bg 0.005035 , bbl scf 2.82
P
P
Eg 35.37 . scf ft 3 2.83
ZT
P
Eg 198.6 . scf bbl 2.84
ZT
E jemplo 2.18. Un pozo de gas está produciendo a una rata de 15,000 ft3/día a partir de
un yacimiento de gas con unas condiciones de 1,000 psi y 100 °F. La gravedad
específica del gas es 0.903. Calcular el flujo de gas a scf/día.
Solución:
T pc = 427 °R
Ppc = 650 psia
Calcular T pr , Ppr y Z .
T pr = 1.3
Ppr = 1.54
Z = 0.748
1,000
Bg 35.37 84.44 scf ft 3
0.748560
Calcular la rata de flujo de gas en scf/día
Cuando la viscosidad incrementa, cada capa de fluido ejerce un gran arrastre friccional
sobre las capas adyacentes y el gradiente de velocidad disminuye.
Donde:
E jemplo 2.19. Dada la siguiente composición del gas, calcular la viscosidad del gas a
3,000 psia y 150 °F.
Composición yi
C1 0.850
C2 0.055
C3 0.035
C4 0.010
Solución:
Composición yi
Mi y i Mi
C1 0.850 16.04 13.634
C2 0.055 30.07 1.654
C3 0.035 44.09 1.543
C4 0.010 58.12 0.5812
M a = 17.412
17.412
Jg 0.6
28.9
T pc = 358.5
Ppc = 672.5
Paso 2. Calcular la viscosidad del gas natural a una atmósfera y 150 °F a partir de la
figura 2.13
P1 0.0119 cp
3,000
Ppr 4.46
672.5
610
T pr 1.70
358.5
Pg
1.7
P1
Pg
Pg P1 1.70.0119 0.0202 cp
P1
E jemplo 2.20. Dada la siguiente composición del gas, calcular la viscosidad a 200 °F y
3,500 psia
Componente yi
N2 0.05
CO2 0.05
H2S 0.02
C1 0.80
C2 0.05
C3 0.03
Solución:
Componente yi Mi yi M i
N2 0.05 28.02 1.4010
CO2 0.05 44.01 2.2005
H2S 0.02 34.08 0.6816
C1 0.80 16.07 12.8560
C2 0.05 30.07 1.5035
C3 0.03 44.09 1.3227
M a = 19.97
19.97
Jg 0.689
28.97
T pc = 380.43 °R
Ppc' 665.11 4400.05 6000.02 1700.05 690.61 psia
P1 0.0123 cp
Paso 4. Emplear los insertos en la figura 2.13, estimar ('P)N2, ('P) C O2, y ('P)H2S
'P N2 0.00042 cp
'P C O 2 0.00026 cp
'P H 2S 0.00005 cp
3500
Ppr 5.21
690.61
660
T pr 1.80
366.53
Pg
1.52
P1
6.1510 log J
3
g
2.86
Donde:
P1 = Viscosidad del gas a presión atmosférica y temperatura del yacimiento, cp.
T = Temperatura del yacimiento, °R.
Jg = Gravedad del gas.
y N 2 , y C O2 , y H 2S = Fracción molar de N2, CO2 y H2S, respectivamente.
ª § Pg ·º
ln «T pr ¨¨ ¸¸» a 0 a1 Ppr a 2 Ppr2 a 3 Ppr3
¬« © P1 ¹¼»
T pr a 4 a 5 Ppr a 6 Ppr2 a 7 Ppr3 2.90
T pr2 a 8 a 9 Ppr a10 Ppr2 a11 Ppr3
T pr3 a 12 a13 Ppr a14 Ppr2 a15 Ppr3
Donde:
T pr = Temperatura pseudo-reducida de la mezcla gaseosa.
Ppr = Presión pseudo-reducida de la mezcla gaseosa.
a0 - a15 = Los coeficientes de las ecuaciones se dan a continuación:
Lee, González y Eakin (1966) presentaron una relación semi-empírica para calcular la
viscosidad de gases naturales. Los autores expresaron la viscosidad del gas en términos
de la temperatura del yacimiento, densidad, y peso molecular del gas. La ecuación tiene
la siguiente forma:
ª § U g ·Y º
Pg 10 K EXP « X ¨¨
4
¸¸ » 2.91
«¬ © 62.4 ¹ »¼
Donde:
9.4 0.02 M a T 1.5
K 2.92
209 19 M a T
986
X 3.5 0.01 M a 2.93
T
La correlación propuesta puede predecir los valores de viscosidad con una desviación
estándar de 2.7% y una desviación máxima de 8.99%. Los autores afirman que el
método puede ser usado para gases agrios.
E jemplo 2.21. Repetir el ejemplo 2.19 y calcular la viscosidad usando el método Lee -
González - Eakin.
Solución:
Z = 0.89
Ug
3,00017.412 8.97 lb / ft 3
61010.730.89
Calcular los parámetros K , X y Y usando las ecuaciones 2.92, 2.93 y 2.94
respectivamente.
K
>9.4 0.02 17.412@6101.5 127.73
209 19 17.412 610
986
X 3.5 0.0117.412 5.291
610
ª § 8.97 · º
1.342
T 8
9
pr
P1 34 10 5
, para T pr d 1.5 2.95
[m
P1
>
166.8 10 5 0.1338 T pr 0.0932 @5
9
, para T pr > 1.5 2.96
[m
T pc
1
6
[m 5.4402 2.97
M P 0.5 2
3
a pc
Pg P1
>
10.8 10 5 EXP 1.439 U r EXP 1.111U r 1.888
@ 2.98
[m
Donde:
Pg = Viscosidad del gas a la presión y temperatura del yacimiento, cp.
P1 = Viscosidad del gas a la presión atmosférica y temperatura del yacimiento, cp.
Ur = Densidad reducida del gas definida por la ecuación 2.104.
Componente yi
C1 0.850
C2 0.053
C3 0.035
C4 0.010
Solución:
M a = 17.412
Ppc = 672.5 psia
T pc = 358.5 °R
Ppr = 4.46
T pr = 1.70
Z = 0.89
[m 5.4402
358.5 6
1
0.4527
17.4120.5 672.52 3
Debido a que T pr > 1.5 aplicamos la ecuación 2.96 para calcular P1.
P1
166.8 10 5 >0.1338 1.7 0.0932@
5
9
0.0121 cp
0.04527
Ur
0.274.46 0.7959
0.891.70
Obtener la viscosidad del gas usando la ecuación 2.98 para dar:
Pg 0.01845 cp
Solución
T pr = (300+460)/440 = 1.73
De la figura 2.16
km
¦ yi k i 3 M i 2.99
¦ yi 3 M i
La regla de la raíz cúbica fue probada en cerca de 17 sistemas con una desviación
promedio de 2.7%.
km
¦ yi k i 3 M i 0.05774
¦ yi M i
3 2.822
-./0.1(231.!45!
La industria petrolera usa cuatro adjetivos para describir poder calorífico, humedo,
seco, bruto y neto, Cuando se usan las palabras humedo y seco, se refieren a la
condición del gas antes de ser quemado. Humedo, significa el gas saturado con vapor de
agua y seco significa que el gas no contiene vapor de agua.
Las palabras bruto y neto se refieren a la condición del agua de combustión después de
ser quemado. El poder calorífico bruto es el calor producido en la combustion completa
bajo presión constante con los productos de combustión enfriados a condiciones
estándar y el agua condensada en estado liquido, esta cantidad también es llamada poder
calorífico total.
El poder calorífico neto es parecido al anterior, excepto que el agua permanece como
vapor a condiciones estándar; la diferencia entre los valores brutos y netos es
precisamente el calor de vaporización del agua de combustión.
La anterior ecuación es usada para calcular el poder calorífico bruto o poder calorífico
neto. Los valores caloríficos para gases reales los podemos obtener dividiendo el poder
calorífico ideal por el factor de compresibilidad del gas a condiciones estándar.
ܮௗ
ܮ ൌ
ܼ
ଶ
ܼ ൌ ͳ െ ൬ ݕ ටͳ െ ܼ ൰
Componentes Composición
F racción molar
CO2 0.0167
N2 0.0032
C1 0.7102
C2 0.1574
C3 0.0751
i-C4 0.0089
n-C4 0.0194
i-C5 0.0034
n-C5 0.0027
C6 0.0027
C7+ 0.0003
1.0000
Solucion:
ଶ
ܼ ൌ ͳ െ ൬ ݕ ටͳ െ ܼ ൰
ܼ ൌ ͳ െ ሺͲǤͲͶͷͻሻଶ
ܼ ൌ ͲǤͻͻͷͺ
Segundo, calcule el poder calorífico bruto del gas como un gas real.
ܮௗ
ܮ ൌ
ܼ
ܷܶܤ
ͳ͵ͳͷǤͻ ܷܶܤ
݂ܿݏ
ܮ ൌ ൌ ͳ͵ʹͳ
ͲǤͻͻͷͺ ݂ܿݏ
El cambio del gas ideal al real, usualmente es menor que el 0.5 % y normalmente es
ignorado.
El poder calorífico neto (seco) puede ser convertido en poder calorífico neto (humedo)
asi:
ܮሺ௨ௗሻ ൌ ሺͳ െ ͲǤͲͳͷሻܮሺ௦ሻ
Donde 0.0175 es la fracción molar de vapor de agua en el gas cuando esta saturado a
condiciones estándar.
Del mismo modo el poder calorífico bruto (seco) puede ser llevado apoder calorífico
bruto (humedo) asi :
Las unidades de la fracción molar son lb-mol del componente j por lb-mol de gas. La
fracción molar se puede convertir en gal/Mscf de la siguiente manera:
yjM j gal
= 19.65 2.100
U oj Mscf
yjM j gal
GPM j 0.3151 2.101
J oj Mscf
Componentes Composición
F racción molar
CO2 0.0167
N2 0.0032
C1 0.7102
C2 0.1574
C3 0.0751
i-C4 0.0089
n-C4 0.0194
i-C5 0.0034
n-C5 0.0027
C6 0.0027
C7+ 0.0003
1.0000
Solución:
yjM j gal
GPM j 0.3151
J oj Mscf
Etano: de 5 a 25 % de recobro
Propano: de 80 a 90 % de recobro
Butanos: 95 % o más
Componentes más pesados: 100%
La figura 2.21 está basada en datos experimentales y muestra la solubilidad del agua en
hidrocarburos líquidos dulces. En hidrocarburos ácidos la solubilidad del agua puede
ser substancialmente alta.
E jemplo 2.26. Determinar el contenido de agua que satura un gas hidrocarburo dulce a
150ºF y 1,000 psia.
Solución:
!
De la figura 2.22
W = C g (220)
!
1. El contenido de agua en CO2 y H2S puro, puede ser significativamente mas alto que
el de un gas natural dulce, particularmente a presiones de alrededor de 700 psia y a
temperatura ambiente.
Correcciones para H2S y CO2 podrán ser aplicadas cuando las mezclas de gas
contengan más del 5% de H2S y/o CO2 a presiones por encima de 700 psia. Esas
correcciones comienzan a incrementarse significativamente a altas concentraciones
y a altas presiones.
2. La adición de pequeñas cantidades de CH4 o N2 a CO2 o H2S pueden reducir
dramáticamente el contenido de agua comparado con un gas ácido puro.
Para un gas con una cantidad de componentes ácidos por debajo del 40%, se puede
estimar el contenido de agua utilizando la ecuación 2.102 y las figuras 2.22, 2.27, y
2.28.
Note que las figuras 2.27 y 2.28 proporcionan valores para los cuales esta determinado
el contenido efectivo de agua para el CO2 y H2S en mezclas de gases naturales, para
usar únicamente en la ecuación 2.102.
Otro método para la estimación del contenido de agua en mezclas de gases ácidos hasta
6,000 psia utiliza las figuras 2.29 y 2.30.
!
Método 1.
Método 2.
Primero la composición debe ser convertida para ser usada con la figura 2.29
F igura 2.23 Contenido de agua en el CO2. F igura 2.24 Contenido de agua en H2S
F igura 2.25 Valores experimentales del contenido de agua F igura 2.26 Contenido de agua saturada en mezclas ricas
en gases ácidos. en CO2 a 100 ºF.
F igura 2.29 Contenido de agua en mezclas de F igura 2.30 Contenido de agua en mezclas
gases ácidos a 2000 psia. de gases ácidos a 6000 psia.
2.3 PR O PI E D A D ES T E R M O D I N A M I C AS
Una vez propuesta su expresión matemática, estas leyes dan lugar a un sistema de
ecuaciones que permiten obtener numerosos resultados y conclusiones prácticas. La
termodinámica tiene una aplicación universal, dado que físicos, químicos e ingenieros la
emplean por igual. Los principios básicos son siempre los mismos, pero sus
aplicaciones difieren. El ingeniero de procesos debe estar capacitado para resolver una
amplia variedad de problemas; entre los más importantes se pueden citar el cálculo de
los requerimientos de calor y trabajo en el procesamiento del gas, así como en la
determinación de las condiciones de equilibrio para la transferencia de los componentes
hidrocarburos entre fases.
Al aplicar la primera ley a un proceso dado, la esfera de influencia del proceso se divide
en dos partes: el sistema y sus alrededores. El primero se refiere a la parte en la que
ocurre el proceso; todo lo demás, que no está incluido en el sistema, constituye los
alrededores.
La elección del signo empleado con Q y W depende del sentido de la transferencia que
se tome como positivo.
El primer término de la ecuación 2.103 puede ampliarse para mostrar los cambios de
energía en diferentes formas. Si la masa del sistema es constante y sólo intervienen
cambios en las energías interna, cinética y potencial,
La elección convencional de los signos en el lado derecho de la ecuación 2.104 hace que
el valor numérico del calor sea positivo cuando se transfiere de los alrededores al
sistema, y que el valor numérico del trabajo sea positivo en el sentido opuesto. Con este
entendimiento, la ecuación 2.104 queda:
En palabras, la ecuación 2.105 establece que la energía total del sistema es igual al calor
agregado al sistema menos el trabajo realizado por éste. Esta ecuación se aplica a los
cambios que ocurren durante un cierto tiempo en sistemas cerrados.
Con frecuencia, en estos sistemas los procesos no presentan cambios en las energías
cinética y potencial, sino solamente variaciones en la energía interna. En este caso, la
ecuación 2.105 se reduce a:
'U = Q - W 2.106
En un sistema fluyendo ( abierto) está involucrada una cantidad adicional de energía PV.
En la entrada del sistema ingresa una cantidad ( V1) de fluido a una presión P1. De esta
forma hay involucrada una energía adicional P1 V1. En forma similar las condiciones del
fluido a la salida del sistema son P2 y V2, lo cual conlleva a una cantidad de energía P2 V2
a la salida.
H U PV 2.107
De esta manera un sistema abierto puede ganar o perder energía en cuatro formas:
E ntonces,
Hay varios métodos básicos que pueden ser usados para calcular los cambios de entalpía
entre dos puntos en un sistema.
1 Calcular los cambios de calor sensible para un proceso en el cual no ocurre un
cambio de fase.
Calor sensible. Es aquel que no involucra cambio de fases sino cambio en la
temperatura del sistema.
2 Cuando hay un cambio de fase calcular el calor latente con una correlación
apropiada.
Calor latente. El calor latente es la cantidad de entalpía (energía) necesaria para
cambiar de fase de líquido a fase vapor. En otras palabras, es la energía que se
tiene que suministrar para vaporizar el líquido.
3 Calcular directamente el ' H a partir de una correlación para componentes puros
o para mezclas, basada en las propiedades del fluido.
4 Calcular el ' H a partir de una ecuación de estado PVT y ecuaciones
termodinámicas.
9 Componentes puros.
H TP H 00 H 0
T
H 00 H T0 H TP 2.109
Donde:
H 0
T
H 00 = La entalpía en el estado de gas ideal con respecto al datum, H 0, a la
Temperatura deseada, Btu / lb-mol.
P
H T H T = Efecto de la presión sobre la entalpía, dado por la diferencia de
0
H TP
H T0 H T0 H TP 2.110
H
H0 H0 H 2.111
La ecuación 2.112, proporciona la relación que fue usada para correlacionar los datos.
Esta correlación consiste de dos partes. Un parte da el efecto de la presión sobre la
entalpía para un fluido simple (un fluido simple es aquel que tiene un factor acéntrico
igual a cero). La segunda parte es una corrección por la desviación que tiene un fluido
real del comportamiento de fluido simple.
H 0
H RT c ®«
°ª H 0 H º 0
» Z«
ª H 0 H º ½°
´
» ¾
2.112
°̄¬ RT c ¼ ¬ RT c ¼ °¿
Donde:
ª H0 H º
«
0
ª H0 H º
«
´
En ambas figuras, si el fluido es un líquido, las funciones son leídas a partir de las
isotermas en el tope de la gráfica. Las presiones críticas ( P c), las temperaturas críticas
(T c), y los factores acéntricos de los componentes puros se dan en la tabla de constantes
físicas.
Las entalpías en el estado de gas ideal están en unidades de Btu/lb; por lo tanto, el peso
molecular del componente se debe usar para convertir la entalpía en el estado de gas
ideal o el efecto de la presión sobre la entalpía antes de sustituirlo en la ecuación 2.111.
9 Mezclas definidas.
H m0 ¦ yi H i0 2.113
i
Tc m ¦ y i Tc i 2.114a
i
Pc m ¦ yi Pc i 2.114b
i
Zm ¦ y i Zi
i
2.115
0 0 c½
°ª H H º ªH0 H º °
H Hm » Zm «
0
RT cm ®« » ¾ 2.116
°̄¬ RT c ¼ ¬ RT c ¼ °
¿
La entalpía de una mezcla de dos fases puede calcularse primero realizando un cálculo
flash, para obtener las moles y la composición de cada fase a las condiciones deseadas.
Entonces la entalpía de cada fase podrá ser calculada como se describió arriba. La
entalpía molar de la mezcla total se calcula por la combinación de las entalpías de la
fase gas y líquida en base molar.
Las cartas de entalpía total se incluyen en las figuras 2.37 a la 2.45 para proporcionar
un método rápido de calcular los balances de calor. Estos pueden reemplazar el cálculo
detallado por componentes para las entalpías del vapor, líquido o entalpías de mezclas
vapor - líquido. Las cartas cubren un amplio rango de composición, temperatura, y
presión encontradas en la mayoría de los sistemas de gas natural desde los separadores
de cabeza de pozo hasta los sistemas de gas natural licuado.
La entalpía en estado de gas ideal fue corregida por el efecto de la presión en la entalpía
del fluido simple y del fluido real interpolando los datos usados en la obtención de las
figuras 2.35 y 2.36. Cálculos de presión se realizaron a partir de presiones reducidas de
0.2 a 3,000 psia en incrementos de 100 psia. El rango de temperaturas usado fue de ±
300 ºF o T r = 0.35 hasta 600 ºF en incrementos de 50 ºF.
Las cartas de entalpía total pueden ser usadas para cálculos rápidos al chequear balances
de calor en diseños de plantas. Aunque el valor absoluto de la entalpía para una
corriente leído de las cartas puede ser diferente al valor calculado por métodos más
rigurosos, los balances en intercambiadores de calor, torres, etc., realizados por ambos
métodos serán similares. Para cálculos de diseño, se recomienda el método más riguroso
por componentes.
Se requirió una ligera extrapolación de las entalpías (en el estado de gas ideal) para las
fracciones para valores de ±20 ºF. Las entalpías en el estado de gas ideal para dos
temperaturas se muestran en las columnas 4 y 10 de la figura 2.31. Note que la entalpía,
en el estado de gas ideal, para cada componente fue multiplicada por el peso molecular
de cada componente para obtener unidades de Btu/lb-mol (columnas 5 y 11). Las
entalpías en el estado de gas ideal para una mezcla de fluidos fue calculada usando la
ecuación 2.113, por ejemplo, para el gas tenemos:
Las temperaturas y presiones reducidas de las mezclas fueron calculadas por T/T cm y
P/P cm y están dadas en las filas 19 y 20. Los valores del efecto de la presión sobre la
entalpía del fluido simple [(H 0 - H)/RT c] 0, y la desviación del efecto de fluido simple
[(H 0 - H)/RT c]´, fueron leídos de las figuras 2.35 y 2.36.
Estos valores se encuentran en las filas 22 y 23 de la figura 2.31 para las tres corrientes.
La ecuación 2.116 fue usada para calcular el efecto de la presión sobre las entalpías de
las tres corrientes, las cuales se dan en la fila 25. Los cálculos para el gas del separador
se muestran a continuación:
H 0
H
m RTcm 0.730 0.0260.020 1.986371.00.7305 538.4 Btu / lb mol
La entalpía de la mezcla para cada corriente fue calculada con la ecuación 2.111,
restando el efecto de la presión, a la entalpía en el estado de gas ideal. Los resultados
para cada corriente se muestran en la fila 27. La entalpía total de la corriente de salida
fue calculada combinando las entalpías de las fases líquidas y vapor en base molar:
Usando las cartas de entalpía total, la entalpía de cada corriente puede ser calculada
como sigue:
2. Leer la entalpía de cada fase en Btu/lb a partir de las cartas de entalpía total con el
peso molecular, la temperatura y la presión.
La entalpía total de la corriente de salida por mol de mezcla a la entrada del separador
(ver figura 2.31, fila 33) es:
La entalpía total de la corriente de entrada (ver figura 2.31, fila 33) es:
Dado que el calor y el trabajo son dos formas de energía, es un hecho que existe una
diferencia en su calidad. Parece ser que el trabajo es de una calidad superior al calor.
Los siguientes hechos lo confirman:
Se han hecho esfuerzos para tratar de convertir toda la energía transferida como
calor en trabajo, y no se ha logrado, a pesar de todas las precauciones técnicas
tomadas. Las eficiencias son sólo del orden del 40 %
2.3.2.2 E ntropía
La temperatura cambia cuando la presión se reduce debido a una corriente de gas que
pasa a través de una restricción, válvula, choke, o perforaciones en el casing. Esto se
conoce como el efecto de Joule ± Thompson. El cambio en la temperatura está
directamente relacionado con la atracción de las moléculas.
La ecuación
§ wV ·
T¨ ¸ V
© wT ¹ P
'T 'p 2.122
CP
esta ecuación es válida cuando el cambio en presión es adiabático, es decir que no hay
entrada ni salida de calor del sistema, y cuando el cambio en presión es pequeño.
PV ZnRT
VT § wZ ·
¨ ¸
Z © wT ¹ P
'T 'P 2.123
CP
Los términos V,T ,Z y C p son siempre positivos. De esta forma, la dirección del cambio
en temperatura depende del signo del término 'P y de la derivada wZ / wT P . Los
ingenieros de petróleos usualmente se enfrentan a situaciones en las cuales la presión
disminuye. El flujo de gas a través de una restricción como una válvula o un choke son
ejemplos comunes.
Los volúmenes específicos de líquidos son muy fáciles de medir. Para la mayoría de los
líquidos orgánicos comunes, hay por lo menos un dato experimental. Hay un número de
referencias donde se pueden encontrar estos volúmenes experimentales tabulados (o
densidades) o en las cuales se dan constantes que permiten calcularlos por medio de una
ecuación empírica.
Se pueden usar métodos aditivos, donde a cada elemento y a ciertos enlaces se les
asignan valores numéricos de modo que se puede calcular el volumen molar en el punto
normal de ebullición por medio de la adición de estos valores. El método más práctico
es sin embargo el método de Tyn y Calus. Vb, el volumen en el punto normal de
ebullición se relaciona al volumen crítico, V c así:
1.048
Vb 0.285 Vc 2.126
Vb y V c, tienen unidades de centímetros cúbicos/gr. mol. Esta sencilla relación tiene una
exactitud general del 3%, excepto para los gases permanentes con bajo punto de
ebullición (He, H2, Ne, Ar, Kr) y algunos compuestos polares con fósforo o nitrógeno
(HCN o PH3), el error promedio es de solo el 2%.
E jemplo 2.30: Estimar la densidad del metano líquido en el punto normal de ebullición
(111.6 K) y compararla con el valor experimental de 26.47 lb/ft3.
Solución:
ft 3 ª 30.483 cm 3 º ª 1 lb º ª gr º cm 3
Vc 0.0988 « » 16.043 = 99.1
lb ¬ 1 ft 3 ¼ «¬ 453 gr »¼ «¬ gr . mol »¼ gr . mol
cm 3
Vc 99.1
gr . mol
T abla 2.3
Propiedades físicas
T abla 2.3
Propiedades físicas
(Continuación)
T abla 2.3
Propiedades físicas
(Continuación)
T abla.2.4.diagrama ph-propano
cm 3
0.28599.1
1.048
Vb 35.22
gr . mol
m asa gr 1 gr . mol
Densidad ; U b 16.043
volum en gr . mol 35.22 cm 3
gr lb
U b 0.4555 ; U b 28.42
cm 3 ft 3
Donde:
R = Constante del gas, 10.73 ft3 psia/lb-mol ºR
T r = Temperatura reducida, T/T c
T c = Temperatura crítica
P c = Presión crítica
Z RA = Constante para cada componente y sus valores se encuentran tabulados. Si
no se conocen los valores de Z RA se pueden estimar por medio de
Donde:
Z : es el factor acéntrico de cada componente
I
VS VSR Z RA 2.129
Donde:
VSR = volumen de referencia
I Está dado por:
I 1 Tr 2 7 1 TrR 2 7 2.130
Note que la ecuación 2.127 no predice el volumen correcto en el punto crítico a menos
que Z RA = Z c
Donde:
Z RAm ¦ xi Z RAi
i
T r = T/T cm , se recomiendan las reglas de Prausnitz y Chueh para T cm:
Tcm ¦ ¦ I i I j Tc ij 2.131 a
i j
Donde:
x i V ci
Ii y Tc ij 1 k T
ij ci Tcj 1 2
¦ x i V ci
Para lo cual se tiene:
8V ci V cj
12
1 k ij
V 13
ci V cj1 3
3
9 Densidades de H idrocarburos.
Kw
punto de ebulliciónprom edio, º R 1 3
gra vedad específica a 60 º F
1. Hacer una tabla de cálculo como la que se muestra en el siguiente ejemplo (tabla
2.4).
Densidad del C 3 (ó H 2 S )
Peso de los componente s de C 3 ( ó H 2 S )
2.132
Volumen de los componente s de C 3 ( ó H 2 S )
Peso (N 2 C 2 )
% en Peso de (N 2 C 2 ) 100 2.133
Peso (N 2 C 2 )
5. Si no hay C O 2 vaya al paso 6. Si lo hay, se tiene en cuenta en una base aditiva del
volumen, como se muestra a continuación
Peso de CO 2 Peso de (N 2 C 2 )
Densidad del CO 2 y (N 2 C 2 )
Volumen de CO 2 Volumen de (N 2 C 2 )
Donde:
peso de (N 2 C 2 )
Volumen de (N 2 C 2 ) 2.134
densidad de (N 2 C 2 )
Peso de metano
% de Metano 100 2.135
Peso total
7. Entre a la parte superior de la figura 2.54 con la seudodensidad del paso 4 ó 5, y baje
verticalmente hasta la línea que representa el porcentaje en peso del metano. Lea la
seudodensidad de la mezcla (a 60 ºF y 14.7 psia) en la escala a la derecha de la
carta.
Volumen del C 3+ = Vol. (C 3) + Vol. (n-C 4) +... Vol. (n-C 14) = 1.0128 ft3
0.567
Paso 3: % en Peso del C 2 en C 2 * 100 1.25%
0.567 44.836
3.353
Paso 6: % en Peso del CH 4 * 100 5.1%
66.24
Error = (0.6/41.2)100=1.5%
T abla 2.4 Cálculo de la Densidad Líquida de una Mezcla a 120 ºF y 1760 psia
Puntos de ebullición
% Destilado T emperatura ºF
IBP 100
5 130
10 153
20 191
30 217
40 244
50 280
60 319
70 384
80 464
90 592
EP 640
F igura 2.57 Puntos de Ebullición que Caracterizan a las Fracciones del Petróleo.
Los datos correspondientes se llevan a la figura 2.57. Con la pendiente de 5.49 baje e
intercepte la curva interpolada para 325 ºF en la serie dibujada con líneas punteadas
(MeABP ). Lea una corrección de ±54 ºF en la ordenada y entonces:
E jemplo 2.33 Cálculo del peso molecular para la mezcla del ejemplo anterior.
M promedio 204.387310.1180.75151.88 e>0.00218 731 3.075 0.7515@ 127
9 Propiedades C riticas
x i V ci
Ti
¦ xi V ci
i
TcT ¦i T i Tc i
Donde :
Chueh y Prausnitz propusieron una relación similar. Se define una fracción superficial
Ti por medio de:
xi V ci2 3
Ti
¦ xi Vci2 3
i
El método de Prausnitz y Chueh no es tan sencillo como el método de Li.
E jemplo 2.34 Estimar, por el método de Li, la temperatura crítica verdadera de una
mezcla con 19.3 % molar de metano 47 % de etano y 33.7 % de n ± butano. El valor
experimental de Li es de 354 K. Los valores de T ci y V ci se dan en la siguiente tabla
T ci K V ci , cm3/mol Ti
Metano 190.4 99.2 0.11
Etano 305.4 148.3 0.398
n - Butano 425.2 255 0.492
Error = 2 K
No hay métodos adecuados para predecir el volumen crítico de una mezcla o la presión
crítica de una mezcla. Los métodos existentes son solo aproximados y complejos
La figura 2.60 es una correlación de presión seudocrítica versus pesos moleculares hasta
100 y gravedad específica del líquido. La figura 2.61 es una correlación similar para
MeABP hasta de 800.
Como alternativas de estas figuras, pueden usarse las siguientes ecuaciones tomadas del
API Technical Data Book:
Ambas correlaciones se han evaluado para pesos moleculares entre 80 a 690; punto
normal de ebullición de 70 a 295 ºF y 6.6 a 95 ºAPI.
VABP = 325 ºF
MeABP = 271 ºF
API = 56.8 º
Solución:
De la figura 2.57 con la pendiente de 5.49 se lee una corrección del VABP de
aproximadamente -85 ºF (extrapolando del grupo a mano izquierda).
En la figura 2.59, se lee con la línea de MABP = 240 hasta encontrar la curva
correspondiente a un API de 56.8 º
T pc = 560 ºF
Esto representa una discrepancia del 7% con respecto al resultado gráfico. Para este
fluido de 56.8 º API utilice la figura 2.61 para estimar la presión seudocrítica. El
MeABP es de 271 ºF y entonces:
Ppc
3.12281 10 9 271 460 0.75152.3201
2.3125
386 psia
F actor A céntrico
Z ¦ xi Z i 2.142
E jemplo 2.36 Estimación del factor acéntrico.Una fracción de petróleo tiene un VABP
de 418 ºF, una pendiente de ASTM de 0.75 y una gravedad API de 41 º. Estime su
factor acéntrico.
Para usar la figura 2.62 o la ecuación 2.143 se necesita el punto de ebullición promedio
MeABP; la temperatura seudocrítica (función del MABP); y la presión seudocrítica
(función del MeABP).
Solución:
Para una fracción del petróleo (mezcla indefinida). El factor acéntrico puede estimarse
si se conocen el MeABP, T pc y Ppc por medio de la ecuación derivada del trabajo de Lee
-Kessler
§ Ppc ·
ln ¨¨ ¸¸ 5.92714 6.09648T 1 1.28862 ln T 0.169347T 6
Z © 14.7 ¹ 2.144
15.2518 15.6875T 1 13.4721 ln T 0.43577T 6
Donde T = MeABP/T pc
E jemplo 2.37 Calcule el factor acéntrico para una mezcla de compuestos definidos y
una fracción del petróleo.
Solución:
Del ejemplo 2.35, para una fracción del petróleo con una gravedad API de 56.8 º
MeABP = 271 ºF
T pc = 602 ºF
F igura 2.62 Factor Acéntrico a Partir del Punto Crítico y del Punto de Burbuja.
MABP = 240 ºF
M = 127
Note que los valores seleccionados para Ppc y T pc se basan en las ecuaciones 2.139 y
2.140.
T = (271+460)/(602+460) = 0.6883
Si se sustituye
Z = 0.34
E jemplo 2.38 Estimación de las propiedades físicas de una mezcla de un gas con una
fracción del Petróleo. El cálculo de las propiedades para una mezcla gas - fracción del
petróleo se da en la tabla 2.5. La mezcla consiste de 1.222 moles de la fracción de
petróleo anterior y 1.0 mol de gas.
M = 79.04
Z = 0.202
Gravedad Especifica = 0.7034 = 69.7 ºAPI
MABP = 30.5
T abla 2.5 Cálculo de las Propiedades de una mezcla de fracción Petróleo - Gas.
Presión de V apor
Esta sección cubre los métodos para estimar y correlacionar las presiones de vapor de
líquidos puros cuando la fase vapor de un fluido puro está en equilibrio con la fase
líquida. La igualdad de potencial químico, temperatura y presión en ambas fases
conduce a la ecuación de Clausius ± Clapeyron
d ln Pvp 'H v
2.146
d 1 T R'Z v
B
ln Pvp A 2.147
T
Esta ecuación solo es válida para pequeños intervalos de temperatura y para condiciones
alejadas del punto crítico.
§ ·
ln Pvp
1
h¨¨1 ¸¸ 2.148
© Tr ¹
ln Pc 1.01325
h Tbr 2.149
1 Tbr
La forma lineal de la ecuación anterior ln( Pvp) vs. 1/T r no es satisfactoria para
materiales que se asocian y en general las ecuaciones 2.148 y 2.149 sobrepredicen para
presiones de vapor a temperaturas menores que T b.
Para lograr más exactitud, varios investigadores propusieron formas de tres parámetros.
La expresión de Pitzer es una de las más exitosas:
ln pvpr f 0 Tr Z f 1 Tr 2.150
Lee y Kessler expresaron los parámetros f(0) y f(1) en la forma analítica que se da a
continuación:
6.09648
f 0 5.921714 1.22862 ln Tr 0.169347Tr 2.151
Tr
15.6875
f 1 15.2518 13.4721ln Tr 0.43577Tr 2.152
Tr
Los valores del factor acéntrico se encuentran tabulados para sustancias puras en la
Tabla 2.3 pero si se utiliza la ecuación 2.150, se recomienda utilizar Z de la ecuación
2.144. La forma de Lee - Kessler de la ecuación de Pitzer predice las presiones de
vapor dentro de un 1 a 2% entre T b y T c. Por debajo de T b la exactitud es mucho menor y
en general el valor predicho es menor que el valor real.
B
ln Pvp A 2.153
T C
G ráfico de Cox
Z R 2 Z R1
2.154
Los superíndices Rl y R2 se refieren a las dos sustancias de referencia, que pueden ser
propano y octano o Benzeno y pentafluorotolueno. En la ecuación anterior todas las
presiones de vapor se calculan a la temperatura reducida de la sustancia cuya presión de
vapor se va a predecir. Las presiones de vapor reducidas de las sustancias de referencia
se calculan por medio de la ecuación de Wagner:
aW bW 1.5 cW 3 dW 6
ln Pvpr 2.155
Tr
Z R1 Z Z R 2 2.156
La mayoría de las ecuaciones que se han presentado son para estimar la presión de
vapor. Es decir, si se tienen los datos del punto de ebullición y las propiedades críticas
es posible desarrollar las constantes para poder estimar las presiones de vapor como
función de la temperatura.
Sin embargo, hay ocasiones donde se tienen datos experimentales de la presión de vapor
sobre un amplio rango de temperaturas y puede ser útil almacenar esta información en
forma analítica. Por medio de técnicas estándar de regresión, se pueden determinar los
mejores valores de las constantes en cualquiera de las ecuaciones anteriores para estimar
la presión de vapor.
aW bW 1.5 cW 3 dW 6
ln Pvpr 2.157
Tr
B DPvp
ln Pvp A C ln T 2 2.158
T C T
ln Pvpr
A
Tr
>
1 Tr2 B 3 Tr 1 Tr
3
@ 2.159
Ajusta mejor los datos experimentales que la ecuación de Antoine. Las tres constantes
de ajuste son: A, B y T c
E jemplo 2.39 Calcule el factor acéntrico a partir de la presión de vapor de la mezcla del
Ejemplo 2.37
Solución:
Para una fracción con un MeABP 271 ºF, T pc = 602 ºF y Ppc = 386 psia, estimar la
presión de vapor reducida a T r = 0.7.
M(C 9) = 128
Mm (mezcla) = 127
De la ecuación 2.141
Si este valor se compara con el valor más exacto de 0.34, que se obtuvo a partir de la
ecuación 2.144, se concluye que la caracterización de una mezcla con amplios puntos de
ebullición por medio de un solo compuesto parafínico no es confiable.
3. E C U A C I O N ES D E EST A D O .
Una ecuación de estado (EOS) es una expresión analítica que relaciona la presión P,
temperatura T y volumen V. Una descripción apropiada de esta relación PVT para
fluidos hidrocarburos reales es esencial para determinar el comportamiento de fase y
volumétrico de los fluidos del yacimiento y en predecir el desempeño de las facilidades
de separación en superficie.
El ejemplo más simple y conocido de una ecuación de estado, es la ecuación del gas
ideal, expresada matemáticamente por la siguiente expresión:
RT
P 3.1
VM
Al desarrollar la EOS para el gas ideal, ecuación 3.1, se realizaron dos suposiciones:
Primera suposición: el volumen de las moléculas de gas es insignificante comparado
con el volumen del recipiente y la distancia entre las moléculas.
Segunda suposición: no hay fuerzas de atracción o repulsión entre las moléculas o
paredes del recipiente.
Van der Waals intentó eliminar estas dos suposiciones al desarrollar una ecuación de
estado empírica para gases reales. En su intento de eliminar la primera suposición, Van
der Waals estableció que las moléculas de gas ocupan una fracción significante de
volumen a elevadas presiones y propuso que el volumen de las moléculas, denotado por
el parámetro b, sea sustraído del volumen molar real V M, en la ecuación 3.1:
RT
P 3.2
VM b
RT a
P 2 3.3
VM b VM
Donde
Los parámetros a y b son constantes que caracterizan las propiedades moleculares de los
compuestos individuales. El símbolo a es una medida de las fuerzas de atracción
intermolecular entre las moléculas. La ecuación 3.3 muestra las siguientes
características importantes:
A elevadas presiones ( P o f), el volumen V M llega a ser muy pequeño a tal punto
que se aproxima al valor de b, que representa el volumen real molecular.
La EOS de Van der Waals o cualquier otra ecuación de estado se puede expresar en una
forma más generalizada así:
Al determinar los valores de las dos constantes a y b para cualquier sustancia, Van der
Waals observó que la isoterma crítica tiene una pendiente horizontal y un punto de
inflexión en el punto crítico, como se muestra en la figura 3.1.
ª wP º ª w2 P º
« wV » 0 , « 2» 0 3.5
¬ ¼ TC ¬ w V ¼ TC
ª wP º RT 2a
« wV » 0 3.6
¬ ¼ TC Vc b 2
V c3
ª w2 P º 2 RT c 6a
« 2» 0 3.7
¬ w V ¼ TC Vc b 3
V c4
§1·
b ¨ ¸V c 3.8
©3¹
§8·
a ¨ ¸ RT c V c 3.9
©9¹
La ecuación 3.8 sugiere que el volumen de las moléculas b sea aproximadamente 0.333
del volumen crítico de la sustancia. Estudios experimentales revelan que b está en el
rango de 0.24 - 0.28 del volumen crítico.
La ecuación 3.10 muestra que sin importar el tipo de sustancia, la EOS de Van der
Waals produce un factor de compresibilidad crítico universal Z c de 0.375. Estudios
experimentales muestran que los valores de Z c para sustancias varían en el rango de 0.23
- 0.31.
Para obtener una expresión más conveniente que calcule los parámetros a y b, la
ecuación 3.10 puede combinarse con las ecuaciones 3.8 y 3.9 de esta forma:
R 2 Tc2
a :a 3.11
Pc
RT c
b :b 3.12
Pc
Donde
R = constante del gas, 10.73 psia-ft3 /lbmol-ºR
P c = Presión crítica, psia
T c = temperatura crítica, ºR
:a = 0.421875
:b = 0.125
3 § RT · 2 § a · § ab ·
VM ¨ b ¸V M ¨ ¸V M ¨ ¸ 0 3.13
© P ¹ ©P¹ ©P¹
La ecuación 3.13 usualmente se conoce como la ecuación cúbica de Van der Waals de
dos parámetros. Los dos parámetros se refieren a las constantes a y b. El término
ecuación cúbica implica una ecuación que al ser expandida, contiene términos de
volumen elevados a la primera , segunda y tercera potencia.
Tal vez la característica más significante de la ecuación 3.13 es que describe el
fenómeno de condensación de líquido y el paso desde la fase gas a la fase líquida
cuando el gas es comprimido.
La ecuación 3.13 se puede expresar en una forma más práctica en términos del factor de
compresibilidad Z reemplazando el volumen molar V M en la ecuación 3.13 por ZRT/P
para dar:
Z 3 1 B Z 2 AZ AB 0 3.14
Donde:
aP
A 3.15
R 2T 2
bP
B 3.16
RT
Z = factor de compresibilidad
P = presión del sistema, psia
T = temperatura del sistema, ºR
La ecuación 3.14 da una raíz real en la región de una fase y tres raíces reales en la
región de dos fases (donde la presión del sistema es igual a la presión de vapor de la
sustancia ). En el último caso, la raíz más grande corresponde a factor de
compresibilidad de la fase vapor Z V, mientras que la raíz positiva más pequeña
corresponde al Z L del líquido.
E jemplo 3.1. Propano puro es mantenido en un recipiente cerrado a 100 ºF. Tanto el gas
como el líquido están presentes. Calcular, empleando la EOS de Van der Waals, la
densidad de las fases líquida y gaseosa.
Solución.
Paso 1. Determinar la presión de vapor P V del propano a partir del gráfico de Cox
(sección 2.4). Esta es la única presión a la cual las dos fases pueden existir a la
temperatura especificada.
P V = 185 psia
a 0.421875
10.73 666
2 2
34,957.4
616.3
10.73666
b 0.125 1.4494
616.3
A
34,957.4185 0.179122
10.732 5602
B
1.4494185 0.044625
10.73560
Paso 4. Sustituir los valores de A y B en la ecuación 3.14 para obtener:
Paso 5. Resolver el polinomio de tercer grado sacando la raíz más grande y más
pequeña del polinomio empleando un método numérico conveniente, tenemos:
Z V = 0.843504
Z L = 0.07534
Ug
18544.0 1.61 lb/ft3
0.84350410.73560
UL
18544 17.98 lb/ft
3
0.0753410.73560
Estas ecuaciones, muchas de las cuales son una modificación de la ecuación de estado
de Van der Waals, varían en complejidad desde expresiones simples que contienen dos
o tres parámetros, hasta formas complicadas que contienen más de cinco parámetros.
Aunque la complejidad de cualquier ecuación de estado no representa un problema
desde el punto de vista computacional, la mayoría de los autores prefieren mantener la
simplicidad que presenta la ecuación cúbica de Van der Waals, mientras mejoran la
exactitud a través de modificaciones.
Todas las ecuaciones de estado generalmente se desarrollan primero para fluidos puros,
y luego para mezclas mediante el uso de las reglas de mezclas.
Estas reglas de mezclas constituyen un medio simple de calcular parámetros de mezclas
equivalentes a los que presentan las sustancias puras.
Redlich y Kwong (1948) demostraron que por medio de un simple ajuste, el término de
presión de atracción, a/V M2, en la ecuación de Van der Waals, puede considerablemente
mejorar la predicción de las propiedades físicas y volumétricas de la fase vapor. Los
autores reemplazaron el término de presión de atracción por un término dependiente de
la temperatura. La ecuación tiene la siguiente forma:
RT a
P 3.17
V M b V M V M b T 0.5
b = 0.26V c 3.18
Imponiendo las condiciones del punto crítico, (como se expresa en la ecuación 3.5), en
la ecuación 3.17, y resolviendo las ecuaciones resultantes en forma simultánea,
tenemos:
R 2 Tc2.5
a :a 3.19
Pc
RT c
b :b 3.20
Pc
Donde
:a = 0.42747
:b = 0.08664
P c V c = 0.333RT c 3.21
La ecuación 3.21 muestra que la EOS de Redlich - Kwong produce un factor crítico de
compresibilidad Z c de 0.333 para todas las sustancias.
Como se indicó al principio, el factor compresibilidad crítico del gas varía entre 0.23 y
0.31 para la mayoría de las sustancias.
Reemplazando el volumen molar, V M en la ecuación 3.17 por ZRT/P tenemos,
Z 3 - Z2 + (A - B - B2)Z - AB = 0 3.22
Donde
aP
A 3.23
R T 2.5
2
bP
B 3.24
RT
Como sucede en la EOS de Van der Waals , la ecuación 3.22 da una raíz real en la
región de una fase (región de la fase gaseosa o líquida), y las tres raíces reales en la
región de dos fases. En el último caso, la raíz más grande corresponde al factor de
compresibilidad de la fase gas Z V mientras que la raíz positiva más pequeña corresponde
a la del líquido Z L.
Solución
a 0.42747
10.73 666
2 2.5
914,110.1
616.3
b 0.08664
10.73666 1.0046
616.3
A
914,110.1185 0.197925
10.732 5602.5
B
1.0046185 0.03093
10.73560
Paso 2. Substituir los parámetros A y B en la ecuación 3.22 y sacar la raíz más grande y
más pequeña.
Z 3 Z 2 0.1660384Z 0.0061218 0
UL
18544.0 25.7 lb/ft
3
0.052737710.73560
UV
18544 1.688 lb/ft
3
0.80264110.73560
Redlich y Kwong extendieron la aplicación de su ecuación a mezclas de gases y de
líquidos empleando las siguientes reglas de mezcla:
2
am ª¦n x a 0.5 º
«¬i 1 i i »¼
3.25
n
bm ¦ xi bi 3.26
i 1
Donde
Para calcular am y bm en una mezcla de gas con una composición de yi, se deben usar las
ecuaciones 3.25 y 3.26 y reemplazar xi por yi . El factor de compresibilidad de la fase
gaseosa y de la fase líquida está dado por la ecuación 3.22 con los coeficientes A y B
definidos por medio de las ecuaciones 3.23 y 3.24. La aplicación de la ecuación
Redlich y Kwong para mezclas de hidrocarburos se puede ilustrar por medio del
siguiente ejemplo:
E jemplo 3.3 Calcular la densidad de una fase gaseosa con la siguiente composición a
4,000 psia y 160 º F. Emplear la ecuación de estado de R-K.
Componente yi Mi Pc Tc
C1 0.86
C2 0.05
C3 0.05
C4 0.02
C5 0.01
C6 0.005
C7+ 0.005 215 285 825
Solución
am = 241.118
bm = 0.5701225
Paso 2. Calcular A y B
A = 0.8750
B = 0.3428
Paso 3. Determinar Z V
Z V = 0.907
UV
4,00020.89 13.85 lb/ft3
10.736200.907
RT aD
P 3.27
V M b V M V M b
D 1 m1 T
r
0.5 2
3.28
donde:
T r = temperatura reducida, T/Tc
Z = factor acéntrico de la sustancia
Para cualquier componente puro, las constantes a y b en la ecuación 3.27, considera las
mismas restricciones clásicas del punto crítico de Van der Waals (ecuación 3.5).
Resolviendo las ecuaciones resultantes, para obtener lo siguiente:
R 2 Tc2
a :a 3.30
Pc
RT c
b :b 3.31
Pc
Donde :a y :b son los parámetros adimensionales del componente puro para la
ecuación de Soave ± Redlich - Kwong, y tienen los siguientes valores:
:a = 0.42747
:b = 0.08664
Edmister y Lee (1984) mostraron que los parámetros a y b se pueden determinar por un
método más conveniente. Para la isoterma crítica tenemos:
VM Vc
3
V M 3 3Vc V M 2 3Vc2V M Vc3 0 3.32
3 ª RT º 2 ª aD bRT º ª aDb º
V M « »V M « b 2 »V M « 0 3.33
¬ P ¼ ¬ P P ¼ ¬ P »¼
En el punto crítico, las ecuaciones 3.32 y 3.33 son idénticas, y D = 1. Igualando los
términos correspondientes en las ecuaciones 3.32 y 3.33, tenemos:
RT c
3V c 3.34
Pc
a bRT c
3V c2 b2 3.35
Pc Pc
ab
V c3 3.36
Pc
b = 0.26 V c
Z 3 Z 2 ( A B B 2 ) Z AB 0 3.37
con
aC P
A 3.38
RT 2
bP
B 3.39
RT
Donde
E jemplo 3.4 Repetir el ejemplo 3.1 y calcular la densidad de las dos fases usando la
ecuación SRK.
Solución.
T c = 666.01 ºR
P c = 616.3 psia
Z = 0.1524
T r = 560/666.01 =0.8408
m = 0.7051
a = 35,427.6
b = 1.00471
Paso 6. Calcular los coeficientes A y B por medio de las ecuaciones 3.38 y 3.39
A = 0.203365
B = 0.034658
Z L = 0.06729
Z V = 0.80212
UV = 1.6887 lb/ft3
UL = 20.13 lb/ft3
Para usar la ecuación 3.37 en mezclas, se requieren las reglas de mezcla para determinar
los términos a y b de la mezcla. Soave adoptó las siguientes reglas de mezclas:
aD m
i j
>
¦ ¦ xi x j a i a j D iD j
0.5
k ij @
1 3.40
bm ¦ >xi bi @ 3.41
i
Con
a D m P
A 3.42
RT 2
bm P
B 3.43
RT
k i , j 1 ! k i , j
k i, j 0
k i, j k j ,i
Slot - Petersen (1987) y Vidal - Daubert (1978) presentaron un estudio del significado
de los coeficientes de interacción y las técnicas para determinar sus valores.
E jemplo 3.5 Un sistema de hidrocarburo de dos fases existe en equilibrio a 4,000 psia y
160 º F. El sistema tiene la siguiente composición:
Componente xi yi
C1 0.45 0.86
C2 0.05 0.05
C3 0.05 0.05
C4 0.03 0.02
C5 0.01 0.01
C6 0.01 0.005
C7+ 0.40 0.0005
Solución .
Paso 1. Calcular los parámetrosD, a y b aplicando las ecuaciones 3.28, 3.30 y 3.31
Componente Di ai bi
Cl 0.6869 8,689.3 0.4780
C2 0.9248 21,040.8 0.7225
C3 1.0502 35,422.1 1.0046
C4 1.1616 52,390.3 1.2925
C5 1.2639 72,041.7 1.6091
C6 1.3547 94,108.4 1.9455
C7+ 1.7859 232,367.9 3.7838
Paso 2. Calcular los parámetros de mezcla (a D)m y bm para la fase líquida y gaseosa
aplicando las ecuaciones 3.40 y 3.41
(a D)m = 9,219.3
bm = 0.5680
(a D)m = 104,362.9
bm = 1.8893
Paso 3. Calcular los coeficientes A y B para cada fase aplicando las ecuaciones 3.42 y
3.43
A = 0.8332
B = 0.3415
A = 9.4324
B = 1.136
Paso 4. Resolver la ecuación 3.37 para el factor de compresibilidad de la fase gaseosa.
Z V = 0.9267
Z L = 1.41211
Paso 6. Calcular el peso molecular aparente de la fase gaseosa y de fase líquida a partir
de su composición.
M a = 20.89
M a = 100.25
UV
4,00020.89 13.556 lb/ft
3
10.736200.9267
UL
4,000100.25 42.68 lb/ft3
10.736201.4121
Sim y Daubert (1980) señalaron que debido a que los coeficientes de la ecuación 3.44 se
determinaron analizando los datos de presión de vapor de hidrocarburos de bajo peso
molecular, es poco probable que la ecuación 3.44 se aplique para fracciones de petróleo
con alto peso molecular.
Teniendo en cuenta que los factores acéntricos de las fracciones pesadas del petróleo se
pueden calcular empleando la correlación de Edmister (ecuación 2.143) o la correlación
de Lee-Kesler (ecuación 2.144), los autores propusieron las siguientes expresiones para
determinar el parámetro m:
Donde
k ijf > H i
2
@
H j / 2H i H j 3.51
L
Vcorr V M L ¦ xi c i 3.53
i
v
Vcorr VM v ¦ yi ci 3.54
i
Donde:
Los autores propusieron seis diferentes esquemas para calcular el factor de corrección c i
de cada componente. Para fluidos del petróleo e hidrocarburos pesados, Peneloux y sus
colaboradores sugirieron que el mejor parámetro de correlación para c i es el factor de
compresibilidad de Rackett Z RA. El factor de corrección es definido matemáticamente
por la siguiente relación:
Donde:
Donde:
Zi = factor acéntrico del componente i.
E jemplo 3.6 Trabaje de nuevo el ejemplo 3.5 incorporando la corrección al volumen de
Peneloux en la solución.
Solución
Paso 1. Calcular el factor de corrección c i utilizando la ecuación 3.52.
VM v
10.736200.9267 1.54119 ft3/mol
4000
VM L
10.736201.4121 2.3485 ft3/mol.
4000
Paso 3. Calcular los volúmenes de gas y líquido corregidos utilizando las ecuaciones
3.53 y 3.54
L
Vcorr 2.3485 0.38564 1.962927 ft3/mol
v
Vcorr 1.54119 0.02349 1.5177 ft3/mol.
v
Z corr
40001.5177 0.91254
10.73620
L
Z corr
40001.962927 1.18025
10.73620
Paso 5. Determine las densidades corregidas de ambas fases.
Uv
400020.89 13.767 lb/ft
3
10.736200.91254
UL
4000100.25 51.07 lb/ft
3
10.736201.18025
Peng y Robinson (1975) llevaron a cabo un estudio general para evaluar el uso de la
ecuación SRK en la predicción del comportamiento de sistemas de hidrocarburos. Ellos
mostraron la necesidad de mejorar la habilidad de la ecuación de estado para predecir la
densidad y otras propiedades del fluido particularmente en las cercanías de la región
crítica. Como una base en la creación de un modelo mejorado, Peng y Robinson dieron
la siguiente expresión:
RT aD
P
V b V b 2 cb 2
RT aD
P 3.57
V b V V b bV b
Asumiendo las condiciones clásicas del punto crítico (ecuación 3.5) en la ecuación 3.57
y resolviendo para los parámetros a y b, se obtiene lo siguiente:
R 2 Tc2
a :a 3.58
Pc
RT c
b :b 3.59
Pc
Donde :
:a = 0.45724
:b = 0.07780
D 1 m1 T
r
0.5 2
3.60
Donde:
m 0.3746 1.5423Z 0.2699Z 2
Esta expresión posteriormente fue expandida por los investigadores (1978) para dar la
siguiente relación:
Z 3 B 1Z 2 A 3B 2 2 B Z AB B 2 B 3 0 3.62
Donde A y B se dan en las ecuaciones 3.38 y 3.39 para componentes puros y por las
ecuaciones 3.42 y 3.43 para mezclas.
Solución.
Paso 1. Calcular los parámetros de la mezcla, (a D)m y bm para la fase gaseosa y la fase
líquida,
(a D)m = 10423.54
bm = 0.862528
(a D)m = 107325.4
bm = 1.696543
A = 0.94209
B = 0.30669
Z v = 0.8625
Z L = 1.2645
Uv
400020.89 14.566 lb/ft
3
10.736200.8625
UL
4000100.25 47.67 lb/ft
3
10.736201.2645
Donde i se refiere a los componentes principales N2, CO2 ó CH4; y j se refiere al otro
componente hidrocarburo de la pareja binaria. Los coeficientes G0, G1 y G2 que
dependen del factor acéntrico se determinan para cada par binario aplicando las
siguientes expresiones:
G1 >
0.584474 1.328 logZ j 2.035767 logZ j @2
3.65
y
G2 > @
2.257079 7.869765 logZ j 13.50466 logZ j 8.3864 logZ j
2
> @3
3.66
G1 >
0.48147 3.35342 logZ j 1.0783 logZ j @
2
3.69
G2 >
0.4114 3.5072 logZ j 0.78798 logZ j @
2
3.70
Para los parámetros de interacción del CO2, la siguiente corrección por presión es
sugerida:
La aplicación de las anteriores correlaciones se debe realizar con precaución para casos
donde T rj excede el valor de 1.
Para reproducir las presiones de vapor a temperaturas reducidas por debajo de 0.7,
Stryjek y Vera posteriormente modificaron el parámetro m en la ecuación PR
introduciendo un parámetro ajustable m1 característico de cada compuesto, en la
ecuación 3.60. Ellos propusieron la siguiente relación generalizada para el parámetro
m:
Donde:
Tr = temperatura reducida del componente puro
Para todos los componentes con una temperatura reducida por debajo de 0.7, Stryjek y
Vera recomendaron situar m1 = 0. Para componentes con una temperatura reducida
mayor de 0.7, los valores óptimos de m1 para compuestos de interés industrial se dan en
la tabla 3.1.
Debido a la naturaleza totalmente empírica de m1, Stryjek y Vera no pudieron encontrar
una correlación generalizada para el parámetro m1 en términos de los parámetros del
componente puro. Ellos establecieron que los valores de m1 dados en la tabla 3.1 no se
deben cambiar.
ci S i bi 3.77
Donde:
Si = parámetro adimensional que también es llamado el parámetro traslación
bi = co - volumen de Peng - Robinson: esta dado por la ecuación 3.59
T abla 3.1
Parámetro " m1 " de compuestos puros
Compuesto m1 Compuesto m1
Nitrógeno 0.01996 Nonano 0.04104
Dióxido de Carbono 0.04285 Decano 0.04510
Agua -0.06635 Undecano 0.02919
Metano -0.00159 Dodecano 0.05426
Etano 0.02669 Tridecano 0.04157
Propano 0.03136 Tetradecano 0.02686
Butano 0.03443 Pentadecano 0.01892
Pentano 0.03946 Hexadecano 0.02665
Hexano 0.05104 Heptadecano 0.04048
Heptano 0.04648 Octadecano 0.08291
Octano 0.04464
L
Vcorr V M L ¦ xi c i
i 1
v
Vcorr VM v ¦ yi ci
i 1
Donde:
T abla 3.2
Parámetro de translación para hidrocarburos
Componente Si
C1 -0.1540
C2 -0.1002
C3 -0.08501
i-C4 -0.07935
n-C4 -0.06413
i-C5 -0.04350
n-C5 -0.04183
n-C6 -0.01478
La tabla 3.2 proporciona los valores de los parámetros de translación para los
hidrocarburos livianos bien definidos, como los presentados por los autores.
Para los componentes no definidos como el C7+, los autores propusieron la siguiente
expresión para calcular el parámetro de translación:
d
Si 1
M i e
Donde:
M i = peso molecular del componente i
d, e = coeficientes positivos de correlación
El enfoque usual para tales VLWXDFLRQHV HV HO GH ³VLQWRQL]DU´ ORV SDUiPHWURV GH OD
ecuación de estado para mejorar la exactitud de la predicción. Durante el proceso de
³VLQWRQL]DFLyQ´ ODV SURSLHGDGHV FUtWLFDV GH ODV IUDFFLRQHV &7+ y los coeficientes de
interacción binarios son ajustados para obtener una correspondencia razonable con los
datos experimentales disponibles para mezclas de hidrocarburos.
2
ª ª § 520 · 0.5 º º
D «1 m «1 ¨ ¸ »» 3.78
«¬ «¬ © T ¹ »¼ »¼
m
ª M J c7 º a Mc a Mc a Mc a Jc
« »
¬ a1 a 2 M J c 7 ¼
3 7 4 7 5 7 6 7
3.79
a 7J 2 c 7 a8 Jc 7
Donde:
T abla 3.3
Coeficientes a1 - a8 para ecuaciones 3.79 y 3.80
Coeficiente m a b
a1 -36.91776 -2.433525E7 -6.8453198
a2 -5.2393763E-2 8.3201587E3 1.730243E-2
a3 1.7316235E-2 -18.444102 -6.2055064E-6
a4 -1.3743308E-5 3.6003101E-2 9.0910383E-9
a5 12.718844 3.4992796E7 13.378898
a6 10.246122 2.838756E7 7.9492922
a7 -7.6697942 -1.1325365E7 -3.1779077
a8 -2.6078099 6.418828E6 1.7190311
aób a1 a 2 M J c 7 a 3 M J 2 c 7 a 4 M J 3 c 7 a 5 J C M c 7
7
3.80
a 6J a 7J 2
a8 J
C 7 C 7 C 7
Los coeficientes a1 - a8 para los parámetros a y b se dan en la tabla 3.3. Los coeficientes
a, b y m para el dióxido de carbono, nitrógeno, metano y etano también fueron
optimizados usando la metodología anterior. Los valores optimizados son:
Coeficiente
Componente a b m
en la ecuación 3.79
CO2 1.499914*104 0.41503575 -0.73605717
N2 4.5693589*103 0.4682582 -0.97962859
C1 7.709708*103 0.46749727 -0.549765
C2 2.416260*104 0.6690577 -0.6952108
Schmidt y Wenzel (1980) han mostrado que la mayoría de las ecuaciones cúbicas de
estado pueden ser expresadas en una forma generalizada por la siguiente ecuación de
cuatro constantes:
RT a
p 2 3.81
V b V ubV wb 2
Cuando a los parámetros u y w se les asigna ciertos valores, la ecuación 3.81 se reduce a
una ecuación de estado específica. La relación entre u y w para un número de
ecuaciones de estado, presentada por Yu y colaboradores se presenta a continuación:
Z 3 UB B 1Z 2 WB 2 UB 2UB A Z WB 3 WB 2 AB 0 3.82
4. C A R A C T E R I Z A C I Ó N D E L A F R A C C I Ó N P ESA D A D E
H I D R O C A R B U R OS.
Casi todos los sistemas de hidrocarburos que ocurren naturalmente contienen una
cantidad de fracciones pesadas que no están bien definidas y no corresponden a mezclas
de componentes discretamente identificados.
Estas fracciones pesadas se identifican como una "fracción más", por ejemplo C7+.
Una descripción adecuada de las propiedades físicas de las fracciones pesadas y otras
fracciones del petróleo, que no están definidas en mezclas de hidrocarburos, es esencial
para realizar cálculos del comportamiento de fases y estudios de modelamiento
composicional.
4.1 M É T O D OS B ASA D OS E N L A D E T E R M I N A C I Ó N D E L PN A
Una gran cantidad de compuestos hidrocarburos que constituyen los crudos se han
agrupado químicamente en varias series de compuestos. Como se anotó en secciones
anteriores, cada serie consiste de aquellos compuestos similares en sus características y
pesos moleculares. Dentro de una serie, existen compuestos que van desde los más
livianos o químicamente simples, hasta los más pesados o químicamente complejos. En
general, se asume que las fracciones pesadas (no definidas) de hidrocarburos están
compuestas de tres grupos de hidrocarburos:
Parafinas (P)
Naftenos (N)
Aromáticos (A)
Donde:
Donde:
wi = fracción en peso del corte de destilación i.
Donde:
xi = fracción molar del corte de destilación i.
CABP >¦ x T @
i
13 3
bi 4.4
MABP CABP
M e ABP 4.5
2
Los tres parámetros (peso molecular, gravedad específica y VABP /WABP ) son
empleados para estimar el contenido de PNA de la fracción pesada de hidrocarburo, la
cual a su vez se usa para predecir las propiedades críticas y factor acéntrico de la
fracción. Hopke - Lin (1974), Erbar (1977), Bergman (1977) y Robinson - Peng han
empleando el concepto de PNA para caracterizar las fracciones no definidas de
hidrocarburos.
¦ Xi 1 4.6
i P,N , A
¦ >M i X i @ M 4.8
i P, N , A
Donde:
Las ecuaciones 4.6 hasta la 4.8 se pueden escribir en forma matricial, así:
ª 1 1 1 ºªX pº ª 1 º
«>M T @ >M Tb @N >M Tb @A »» «« X N »» « M WABP » 4.9
« b P « »
«¬ >M @P >M @N >M @A »¼ «¬ X A »¼ «¬ M »¼
Robinson y Peng afirmaron que es posible obtener valores negativos para los contenidos
de PNA. Para prevenir estos valores negativos, los autores impusieron las siguientes
restricciones:
0 d X p d 0.90
XN t0
XA t0
Para resolver la ecuación 4.9, y de esta forma hallar el contenido de PNA, se requiere el
punto de ebullición promedio en peso y el peso molecular para el corte de la fracción no
definida de hidrocarburo. Si no se disponen los valores experimentales de estos cortes,
se pueden usar las siguientes correlaciones propuestas por Robinson y Peng.
i 1
>
Grupo parafínico: ln Tb ln 1.8 ¦ a i n 6
6 i 1
@ 4.10
Grupo nafténico:
i 1
>
ln Tb ln 1.8 ¦ a i n 7
6
@i 1
4.11
Donde:
n= número de átomos de carbono en la fracción no definida de hidrocarburo.
ai = Coeficientes de las ecuaciones (ver tabla 4.1)
T abla 4.1
Coeficientes a i en las ecuaciones 4.10 a la 4.12
Pc X P Pc P X N Pc N X A Pc A 4.16
Donde
P c = presión crítica de la fracción pesada de hidrocarburo, psia.
La presión crítica para cada corte de fracción pesado es calculada de acuerdo a las
siguientes ecuaciones:
206.126096 n 29.67136
Grupo parafínico: Pc P 4.17
0.227 n 0.3402
206.126096 n 206.126096
Grupo nafténico: Pc N 4.18
0.227 n 0.1372
206.126096 n 295.007504
Grupo aromático: Pc A 4.19
0.227 n 0.3252
Paso 3. Calcular el factor acéntrico de cada corte de fracción no definida usando las
siguientes expresiones:
Tc X P Tc P X N Tc N X A Tc A 4.23
Donde
S 0.996704 0.00043155 n
S1 0.99627245 0.00043155 n
3 >log Pc 14.7 @
Z 1 4.27
7 >Tc Tb 1@
Donde:
E jemplo 4.1 Calcular la presión crítica, temperatura crítica y factor acéntrico de una
fracción no definida de hidrocarburos con una peso molecular promedio de 94 y un
punto de ebullición promedio en peso de 655 °R. El número de átomos de carbono del
componente es de 7.
Solución:
Paso 1. Calcular el punto de ebullición de cada corte aplicando las ecuaciones 4.10
hasta la 4.12:
Paso 2. Calcular el peso molecular de varios cortes usando las ecuaciones 4.13
hasta la 4.15.
Pc 471.7 psia
Paso 6. Calcular el factor acéntrico para cada corte en la fracción usando las
ecuaciones 4.20 hasta 4.22
Paso 7. Resolver (T c)P , (T c)N y (T c)A usando las ecuaciones 4.24 hasta la 4.26
Tc 964.1 qR
Z 0.3680
wA 8.47 K w 4.28
Donde:
wA = fracción en peso de los aromáticos
K w = fractor de caracterización de Watson, definido matemáticamente por la
siguientes expresión.
Kw Tb 13 J 4.29
Donde:
J = gravedad específica de la fracción no definida.
Tb = Punto de ebullición promedio en peso, °R.
Paso 2. Con la estimación del contenido de aromáticos, las fracciones en peso de los
cortes parafínicos y nafténicos se calculan resolviendo el siguiente sistema de
ecuaciones lineales:
wP w N 1 wA 4.30
wP w N 1 wA
4.31
JP JN J JA
Donde:
wP = fracción en peso del corte de parafina.
wN = fracción en peso del corte de naftenos.
J = gravedad específica de la fracción no definida.
JP , JN , JA = gravedad específica de los tres grupos en el punto de ebullición
promedio en peso de la fracción no definida. Estas gravedades se calculan a
partir de las siguientes relaciones:
4.32
0.3207736 10 9 Tb 460 3
4.33
0.330966 10 9 Tb 460
3
4.34
0.166318 10 9 Tb 460
3
- Parafinas
4.36
0.85103 10 6 Tb 460
3
Z p
0.14 0.0009 Tb 460 0.233 10 6 Tb 460
2
4.37
- Naftenos.
4.38
0.2544 10 5 Tb 460
3
Pc N
726.414 1.3275 Tb 460 0.9846 10 3 Tb 460
2
4.39
0.45169 10 6 Tb 460
3
Z N Z P 0.075 4.40
Bergman asignó los siguientes valores especiales de factor acéntrico a los naftenos C 8,
C 9 y C 10.
C 8 : Z N 0.26
C 9 : Z N 0.27
C 10 : Z N 0.35
- Aromáticos.
4.41
0.82358 10 6 Tb 460
3
4.42
0.23416 10 5 Tb 460
3
Z A Z P 0.1 4.43
Pc X P Pc P X N Pc N X A Pc A 4.44
Tc X P Tc P X N Tc N X A Tc A 4.45
Z X P Z P X N Z N X A Z A 4.46
4.2 O T R OS M É T O D OS PA R A C A R A C T E R I Z A R L AS F R A C C I O N ES
P ESA D AS D E H I D R O C A R B U R OS.
Tc C 7
1.8 961 10 a 4.47
Donde:
(T c)C7+ = temperatura crítica del C7+, °R.
a = coeficiente de la ecuación, el cual se da por la siguiente expresión:
2
a 2.95597 0.090597 n 3
n M C 7 2.0 14 4.48
Pc C 7 10 5Y Tc C 7 4.49
Tb C 7 0.0004347 Tc C 7 265
2
4.50
1. Correlaciones de Rowe.
2. Correlaciones de Standing.
Solución:
n = 15.29
Tc C 7 1,279.8 qR
Pc C 7 230.4 psia
Tc C 7 1,269.3 qR
Pc C 7 270. psia
Whitson (1983) encontró inconsistencia en los datos de peso molecular tabulados por
Katz - Firoozabadi después de analizar y comparar éstos datos con las fuentes a partir de
las cuales se desarrollaron.
Whitson modificó las propiedades físicas originales tabuladas para hacer su uso más
consistente. La modificación se realizó empleando la forma de la correlación de Riazi y
Daubert: para extrapolar los datos de peso molecular desde el C22 al C45.
T aTbbJ c
Donde:
% Desviación
T a b c promedio Máxima
M 4.5673*10-5 2.1962 -1.0164 2.6 11.8
T c ºR 24.2787 0.58848 0.3596 1.3 10.6
P c, psia 3.12281*109 -2.3125 2.3201 3.1 9.3
V c, ft3/lb 7.5214*10-3 0.2896 -0.7666 2.3 9.1
Ahmed (1985) correlacionó las propiedades físicas tabuladas por Katz ± Firoozabadi -
Whitson con el número de átomos de carbono de la fracción usando un modelo de
regresión. La ecuación generalizada tiene la siguiente forma:
a5
T a1 a 2 n a3 n 2 a 4 n 3 4.53
n
Donde
T abla 4.2
T a1 a2 a3 a4 a5 MA E rro
D% r%
M -131.11375 24.96156 -0.34079022 2.4941184*10-3 468.32575 0.418 1.31
T c ºR 915.53747 41.421337 -0.7586859 5.8675351*10-3 -1.3028779*103 0.126 0.5
P c, psia 275.56275 -12.522269 0.29926384 -2.8452129*10-3 1.7117226*103 2.071 5.3
Tb, ºR 434.38878 50.125279 -0.9027283 7.0280657*10-3 -601.85651 0.178 0.7
Z -0.50862704 8.700211*10-2 -1.8484814*10-3 1.4663890*10-5 1.8518106 1.152 3.7
J 0.86714949 3.4143408*10-3 -2.839627*10-5 2.4943308*10-8 -1.1627984 0.103 0.47
V c, ft3/lb 5.223458*10-2 7.8709139*10-4 -1.9324432*10-5 1.7547264*10-7 4.4017952*10-2 0.099 0.49
Tc >
Tb 1 1.25127 0.137242n0.884540633@ 4.54
Donde
n= número de átomos de carbono
Tb = punto de ebullición promedio de la fracción no definida, °R.
E jemplo 4.3. Calcular las propiedades críticas y el factor acéntrico de C 7+ con un peso
molecular medido de 198.71 y gravedad específica de 0.8527. emplear los siguientes
métodos:
1. Correlación de Rowe.
2. Correlación de Standing.
Solución:
n = 14.0507
T c = 1,247.06 °R
P c = 245.089 psia
T b = 941.03 °R
(T c)C7+ = 1,247.73 °R
Paso 2. Determinar la presión crítica de la ecuación 4.52
5. C O N C E P T OS B ÁSI C OS D E C O M PO R T A M I E N T O D E F ASES
Las condiciones bajo las cuales existen estas fases es una materia de considerable
importancia práctica. Las determinaciones experimentales o matemáticas de estas
condiciones se expresan convenientemente en diferentes tipos de diagramas,
comúnmente llamados Diagramas de Fases.
El objetivo de este capítulo es revisar los principios básicos del comportamiento de fase
de los hidrocarburos e ilustrar el uso de los diagramas de fase para describir y
caracterizar el comportamiento volumétrico de los sistemas de un solo componente, dos
componentes y multi-componentes.
5.1 SUST A N C I AS PU R AS
Primero se considerará un sistema que consiste tan solo de una sustancia pura. Estos
sistemas se comportan en forma diferente a como lo hacen los sistemas de dos o más
componentes. En particular, el principal interés en el estudio del comportamiento de
fases, es determinar las condiciones de temperatura y presión para las cuales las
diferentes fases pueden existir.
La línea TC de la figura 5.1 es llamada línea de presión de vapor. Esta línea separa las
condiciones de presión - temperatura para las cuales una sustancia es líquida, de las
condiciones para las cuales es un gas. Los puntos de presión - temperatura que están
sobre esta línea, indican las condiciones para las cuales la sustancia es líquida. De
forma similar, los puntos que se encuentran debajo de ésta línea, representan las
condiciones para las cuales la sustancia es un gas. Los puntos que caen exactamente
sobre la línea indican las condiciones para las cuales coexisten tanto el líquido como el
gas.
E l punto crítico.
Para una sustancia pura, la temperatura crítica puede definirse como la temperatura por
encima de la cual el gas no puede ser licuado a pesar de la presión aplicada. De forma
similar, la presión crítica de una sustancia pura es la presión por encima de la cual el
líquido y el gas no pueden coexistir a pesar de la temperatura.
A temperatura por debajo de la temperatura del punto triple, la línea de presión de vapor
divide las condiciones para las cuales la sustancia es sólida, de las condiciones para las
cuales es gas. Esta línea es llamada línea de presión de sublimación. Teóricamente esta
línea se extiende a la temperatura y a la presión del cero absoluto.
La línea del punto de fusión es la línea casi vertical sobre el punto triple. La línea
separa las condiciones sólidas de las condiciones líquidas. De nuevo los puntos de
presión y de temperatura que caen exactamente en esta línea indican un sistema de
doble fase, en este caso, coexisten el sólido y el líquido.
Para lograr tener un mejor entendimiento de la utilidad del diagrama de fases, considere
una celda en la cual la temperatura puede ser controlada y el volumen varía por la
inyección o remoción de mercurio como se muestra en la figura 5.2. La figura 5.2A
muestra una sustancia pura que ha sido colocada en una celda a una presión P1, y a
alguna temperatura por debajo de la temperatura crítica de la sustancia. La temperatura
es mantenida constante y el volumen de la celda aumentará por la remoción de mercurio
causando una disminución en la presión.
que aquellas moléculas con la energía cinética más alta pueden escapar del líquido y
formar el gas.
Si el mismo proceso es mostrado para una temperatura mayor que la temperatura crítica,
por ejemplo la línea 45 de la figura 5.3, la remoción de mercurio causa que la presión
disminuya. Sin embargo, no habrá un cambio brusco en la densidad de la sustancia. La
línea presión de vapor no será atravesada y no se observará un cambio abrupto de fase.
Note que los términos gas y líquido simplemente designan, la fase menos densa y la
más densa respectivamente. Estas fases se unen y pierden sus identidades en el punto
crítico.
F igura 5.3 Diagrama de fase típico de una sustancia pura, con dos líneas de expansión
isotérmica: 123 debajo de T c, 45 encima de T c .
El resultado del proceso descrito en la figura 5.2 puede ser representado en forma de un
diagrama presión - volumen. La figura 5.4 muestra dos isotermas de un típico diagrama
presión - volumen para una sustancia pura. Los procesos 1-3 y 4-5 corresponden a los
procesos mostrados en la figura 5-3.
F igura 5.4 Diagrama típico presión volumen de una sustancia pura mostrando dos
isotermas: 13 por debajo de T c, 45 por encima de T c.
Considere un proceso que comienza en el punto 1 con una sustancia en la fase líquida.
La temperatura se mantiene constante y el volumen es aumentado por la remoción de
mercurio. Lo anterior causará una reducción en la presión desde P1, hasta P V. Un
cambio relativamente grande en la presión causará un pequeño cambio en el volumen;
debido a que los líquidos son relativamente incompresibles.
La línea 45 en la figura 5.4, ilustra el mismo proceso a una temperatura por encima de la
temperatura crítica de la sustancia. La figura muestra que simplemente hay una
expansión de la sustancia y que no ocurre un cambio abrupto de fase.
Domo de saturación.
La figura 5.5 muestra un diagrama presión - volumen más completo. La línea punteada
muestra el contorno de todos los puntos de burbuja y los puntos de rocío. El área dentro
de la línea punteada indica las condiciones para las cuales el líquido y el gas coexisten.
A menudo esta área es llamada domo de saturación.
La línea del punto de burbuja y la línea del punto de rocío coinciden en el punto crítico.
Nótese que la isoterma en la temperatura crítica muestra un punto de inflexión
horizontal como si ésta pasara a través de la presión crítica.
5.2 M E Z C L AS D E D OS C O M P O N E N T ES
A medida que la presión se disminuye por debajo de la presión de burbuja aparece gas
adicional. Finalmente, sólo una minúscula cantidad de líquido permanece. Este es el
punto de rocío, Pd. Una reducción adicional de la presión, hasta el punto dos
simplemente provoca una expansión del gas.
F igura 5.6 Diagrama típico de una mezcla de dos componentes con una línea expansión
E l punto crítico.
La definición de punto crítico como se usa para una sustancia pura, no se puede aplicar
a una mezcla de dos componentes. En una mezcla de dos componentes, el líquido y el
gas pueden coexistir, a temperaturas y presiones por encima del punto crítico. Observe
que el domo de saturación existe a temperaturas más altas que la temperatura crítica y
presiones más altas que la presión crítica. Podemos ver ahora que la definición de
punto crítico es simplemente aquel en el cual la línea de punto de rocío y de burbuja se
juntan.
Una definición más rigurosa de punto crítico es el punto en el cual todas las propiedades
del líquido y del gas llegan a ser idénticas.
F igura 5.7 Diagrama de fase típico de una mezcla de dos componentes, con sus
respectivas líneas de presión de vapor
La figura 5.7 muestra las líneas de presión de vapor de los 2 componentes de una
mezcla sobre puestas en el diagrama de fases de la mezcla. El domo de saturación para
la mezcla se encuentra entre las líneas de presión de vapor de los dos componentes. La
temperatura crítica de la mezcla se encuentra entre la temperatura crítica de los dos
componentes puros. Sin embargo, la presión crítica de la mezcla usualmente está por
encima de las presiones críticas de ambos componentes. La presión crítica de una
mezcla de dos componentes, usualmente será más alta que la presión crítica de
cualquiera de los dos componentes. La figura 5.8 muestra datos de fase para ocho
mezclas de metano y etano, entre las líneas de presión de vapor del etano puro y metano
puro. De nuevo observe, que el domo de saturación de cada una de las mezclas se
encuentra entre las líneas de presión de vapor de las dos sustancias puras y que las
presiones críticas de la mezcla se encuentran por encima de las presiones críticas de los
dos componentes puros. La línea a trazos representa los puntos críticos del domo de
mezclas de metano y etano.
C ricondentérmica y cricondenbárica.
Condensación retrógrada.
F igura 5.9 Diagrama de fases típico de una mezcla con dos componentes y su
definición de cricondenbárica y cricondertérmica
La figura 5.11 muestra una sola isoterma en un diagrama de presión - volumen de una
mezcla de dos componentes con una, composición constante. La característica que
distingue este diagrama de un diagrama presión - volumen de una sustancia pura, figura
5.4, es el descenso en presión a medida que el proceso pasa desde el punto de burbuja al
punto de rocío. La isoterma dentro del domo de saturación no es horizontal y tampoco
necesariamente recta. El descenso en la presión es causado por los cambios en las
composiciones del líquido y el gas a medida que el proceso atraviesa la región de dos
fases.
F igura 5.11 Diagrama típico de una mezcla de dos componentes mostrando una
isoterma por debajo de la temperatura crítica.
F igura 5.12 Diagrama de presión volumen para una mezcla de 47.6 de porcentaje en
peso de n-pentano y 52.4 de porcentaje en peso de n-heptano
Una serie de experimentos de laboratorio con una sustancia pura (Figura 5.2) arroja
datos para presión, volumen, temperatura. Una serie similar de experimentos con un
sistema de dos componentes, genera datos para variables adicionales (la composición
de la mezcla, del líquido en equilibrio y del gas en equilibrio son todas importantes).
Por lo tanto, además de graficar combinaciones de temperatura, presión y volumen,
gráficos adicionales de estas variaciones contra la composición pueden ser obtenidos y
son de gran utilidad.
La figura 5.13, muestra un diagrama típico de presión - composición para una mezcla de
dos componentes a una temperatura dada. Las combinaciones de composición y
presión las cuales caen encima del sobre indican condiciones para las cuales la mezcla
es completamente líquida. Cuando estos caen por debajo del sobre indican condiciones
a las cuales la mezcla es gaseosa. Cualquier combinación de presión y composición que
caiga dentro del sobre indica que la mezcla existe en dos fases, gas y líquido. El punto
de burbuja y el punto de rocío tienen las mismas definiciones anteriores. La línea de
punto de burbuja es también la localización de composiciones del líquido cuando las
dos fases están presentes.
La línea de punto de rocío es la localización de las composiciones de gas cuando gas y
líquido están en equilibrio.
La línea que une la composición del líquido y la composición del gas en equilibrio es
conocida como línea enlazante de equilibrio. Las líneas enlazantes son siempre
horizontales para mezclas de dos componentes.
La línea más baja del sobre de saturación es la línea de punto de burbuja y la más alta es
la línea de punto de rocío. Las condiciones de temperatura y composición, que caen por
debajo del sobre de saturación indican que la mezcla es completamente líquida. Cuando
la presión es menor que las presiones críticas de ambos componentes las líneas de punto
de burbuja y de punto de rocío se juntan en las presiones de vapor de los componentes
puros a cada lado del diagrama. Cuando la presión excede la presión crítica de uno de
los componentes, la línea de punto de burbuja y punto de rocío se unen en el punto
crítico. Por ejemplo una mezcla de 98% molar de metano y 2% molar de etano tiene
una temperatura crítica de menos de 110qF a una presión crítica de 700 psia. Cuando la
presión de interés excede las presiones críticas de ambos componentes el sobre de fases
presenta dos puntos críticos. Por ejemplo mezclas de metano y etano exhiben puntos
críticos a 900 psia y -62qF y a 900 psia y 46qF. Las líneas enlazantes que dan la
composición del líquido y el gas en equilibrio son siempre horizontales. La línea de
punto de burbuja da la composición del líquido en equilibrio y la línea del punto de
rocío no da la composición del gas en equilibrio. Las longitudes de las líneas enlazantes
representan las cantidades de gas y líquido en el equilibrio, de la misma manera como
en el diagrama presión - composición.
6. E Q U I L I B R I O G AS L Í Q U I D O
El área limitada por las curvas de punto de burbuja y punto de rocío en el diagrama de
fases de una mezcla multicomponente define las condiciones para que el gas y el
líquido existan en equilibrio. Las cantidades y composiciones de las dos fases varían
en diferentes puntos dentro de los límites del domo de saturación.
9 Cálculo de las condiciones para las cuales una mezcla exhibe un punto de rocío.
9 Cálculo de las condiciones para las cuales una mezcla exhibe un punto de burbuja.
6.1 SO L U C I O N ES I D E A L ES
Como en el caso de los gases ideales, las soluciones ideales de líquidos, no existen.
Realmente, las únicas soluciones que se aproximan al comportamiento de una solución
Ideal, son las mezclas de gases a bajas presiones. Mezclas líquidas de componentes de
las mismas series homólogas, se aproximan al comportamiento de la solución ideal
solamente a bajas presiones. Sin embargo, estudios del comportamiento de fase
de las soluciones ideales nos ayudan a entender el comportamiento de las soluciones
reales.
Pj x j Pvj 6.1
Como se trató anteriormente, la ecuación de Dalton puede ser usada para calcular la
presión parcial ejercida por un componente de una mezcla de gas ideal.
Pj yjP
y j Pvj
xj P
La ecuación 6.2 establece una relación entre la fracción molar de un componente
en el gas y su fracción molar en el líquido. Para determinar los valores de xj e yj , esta
ecuación debe ser combinada con otra ecuación que relacione estas dos cantidades.
Tal ecuación puede ser desarrollada considerando un balance de materia sobre el j-
ésimo componente de la mezcla.
Entonces:
z jn x j nL y j ng 6.3
Combinando las ecuaciones 6.2 y 6.3, para eliminar yj resulta la siguiente expresión:
Pvj
z jn x j nL x j ng 6.4
P
z jn
xj 6.5
Pvj
nL ng
P
z jn
¦xj ¦ 1 6.6
j j Pvj
nL ng
P
Las ecuaciones 6.2 y 6.3 pueden ser combinadas para eliminar xj, con el resultado:
z jn
¦ yj ¦ 1 6.7
j j P
ng nL
Pvj
La ecuación 6.6 ó la ecuación 6.7, son utilizadas para calcular las composiciones de las
fases líquida y gaseosa de una mezcla en equilibrio. En cualquiera de los casos, se
requiere una solución por prueba y error. El cálculo es simplificado si un mol de mezcla
total es tomado como base, de tal manera que: ñL = nL/n, ñg = ng/n y ñL + ñg = 1.
Basados en estas definiciones las ecuaciones 6.6 y 6.7 pueden ser reagrupadas y
simplificadas así:
zj
¦xj ¦ 1 6.8
j j § Pvj ·
1 ñg ¨¨ 1¸¸
© P ¹
y
zj
¦ yj ¦ 1 6.9
j j § P ·
1 ñL ¨ 1¸
¨P ¸
© vj ¹
Si la ecuación 6.9 es usada para los cálculos, se deben seleccionar valores de prueba de
ñL. Cuando se halla un valor de ñL que hace que la sumatoria total sea igual a 1.0, los
términos resultantes en la suma representan los valores de yj y ñL representa el número
de moles en el líquido por mol de mezcla total.
E jemplo 6.1. Calcular las composiciones y cantidades de gas y de liquido cuando 1.0
lb-mol de la siguiente mezcla es llevada al equilibrio a 150 ºF y 200 psia.
Suponga comportamiento de solución ideal.
Solución:
zj
¦ yj ¦ 1
j j § P ·
1 ñL ¨ 1¸
¨P ¸
© vj ¹
Para ahorrar espacio, solo se muestra el cálculo con el valor de prueba final de ñ L =
0.487.
La sumatoria es igual a 1.0 para ñL = 0.487; de esta manera hay 0.487 moles de líquido
y 0.513 moles de vapor por cada mol de mezcla total y las composiciones del gas y del
líquido en equilibrio están dadas en las columnas 4 y 5 de la tabla anterior.
6.1.4 C álculo de la Presión del Punto de Burbuja de una Solución L íquida Ideal
zj
¦ 1 6.10
j Pb
Pvj
Pb ¦ z j Pvj 6.11
j
Por lo tanto, la presión del punto de burbuja de una solución líquida ideal a una
determinada temperatura, es simplemente la sumatoria de la fracción molar de cada
componente por su correspondiente presión de vapor.
E jemplo 6.2. Calcular la presión del punto de burbuja a una temperatura de 150 ºF para
la mezcla presentada en el ejemplo 6.1. Asuma comportamiento de solución ideal.
Solución:
Resolver la ecuación
Pb ¦ z j Pvj
j
6.1.5 C álculo de la Presión del Punto de Rocío de una Solución de gas ideal
zj
¦ 1 6.12
j Pvj Pd
ó
1
Pd 6.13
§ zj ·
¦ ¨¨ ¸
¸
j
© Pvj ¹
El punto de rocío de una mezcla de gas ideal a una temperatura dada, es simplemente el
recíproco de la sumatoria de la fracción molar dividida por la presión de vapor para
cada componente.
E jemplo 6.3. Calcular la presión del punto de rocío a 150 ºF para la mezcla dada en el
ejemplo 6.1. Suponga comportamiento de solución ideal.
Solución:
Resolver la ecuación:
1
Pd
§ zj ·
¦ ¨¨ ¸
¸
j
© Pvj ¹
1
Pd 136 psia
0.00738
6.2 SO L U C I O N ES N O I D E A L ES
El uso de las ecuaciones 6.6 ó 6.7 para predecir el equilibrio gas líquido, está
restringido severamente por varias razones. Primero, la ecuación de Dalton está basada
en la suposición de que el gas se comporta como una solución ideal de gases ideales.
Para propósitos prácticos, la suposición de gas ideal limita las ecuaciones a
presiones por debajo de los 100 psia y temperaturas moderadas.
Varios métodos teóricos para superar estos problemas han sido tratados. Pero el único
camino razonablemente exacto para predecir el equilibrio gas - líquido, es a través
de correlaciones basadas en observaciones experimentales del comportamiento del
equilibrio gas - líquido. Estas correlaciones involucran la razón de equilibrio K, la
cual es definida como:
yj
Kj , 6.14
xj
Un análisis de la ecuación 6.2, indica que las ecuaciones 6.6 y 6.7 pueden ser
modificadas, al tener en cuenta el comportamiento de una solución real; simplemente
reemplazando las relaciones de presión por las relaciones de equilibrio
experimentalmente determinadas. Por lo tanto,
z jn
¦xj ¦ 1 6.15
j nL ng K j
z jn
¦ yj ¦ 1 6.16
j n
ng L
Kj
zj
¦xj ¦ 1 6.17
j 1 ñg K j 1
zj
¦ yj ¦ 1 6.18
j § 1 ·
1 ñL ¨ 1¸
¨K ¸
© j ¹
La elección entre las ecuaciones 6.17 y 6.18 es complemente arbitraria; pues ambas
trabajan con igual eficiencia y requieren una solución de prueba y error. Normalmente
se seleccionan valores sucesivos de ñg ó ñL para prueba y error hasta que la sumatoria
sea igual a 1.0.
El valor de ñg o ñL que hace que la suma sea igual a 1.0 es el valor correcto de ñg o ñL.
Por consiguiente, los términos en la suma representan la composición bien sea del
líquido o del gas dependiendo de la ecuación usada. La exactitud de los resultados,
depende estrictamente de los valores de relaciones de equilibrio utilizados.
E jemplo 6.4. Repita el ejemplo 6.1, sin suponer comportamiento de solución ideal.
Suponga que las cartas de relación de equilibrio dadas en el capítulo 25 del libro de la
GPSA pueden usarse para esta mezcla.(Nota: utilice las cartas que tiene una presión de
convergencia de 5000 psia)
Solución:
zj
¦ yj ¦ 1
j § 1 ·
1 ñL ¨ 1¸
¨K ¸
© j ¹
Determine los valores de K j a 150 ºF y 200 psia, use las cartas del capítulo 25 del libro
de la GPSA. Determine un primer valor de prueba de 0.5 para ñL. Esto dará como
resultado una sumatoria de 0.992; entonces pruebe con un valor de 0.6 para ñL. Este
dará como resultado una suma de 1.010.
La sumatoria se hace igual a uno cuando ñL = 0.547, los cálculos son mostrados en la
siguiente tabla.
Por lo tanto, hay 0.547 moles de líquido por cada mol de mezcla total, y las
composiciones de las fases gas y líquido son:
El análisis de la forma como las ecuaciones 6.17 y 6.18 son afectadas por los valores
incorrectos de prueba de ñg o ñL, proporciona una guía para determinar rápidamente la
convergencia a los valores correctos. La figura 6.1 muestra los resultados de los
cálculos realizados con la ecuación 6.17.
ido
íqu
tod
ol
sf
s
a se
lo
p do
r eje
m
plo,
ejem
o
b, p
T< T r
P> Pbó
Td, po
T<
<
Pd<P<PbyTb
P
<P
d
ó T
>T
d
,p
or
ej oo
to
ds
a
g
eml
p
Valor de prueba de ñg
Por otro lado, la mezcla será todo gas si la temperatura elegida está por encima de la
temperatura de punto de Rocío, o si la presión escogida es menor que la presión de
punto de Rocío. La línea a trazos inferior en la figura 6.1, muestra que no se obtiene
zj
solución para estas condiciones. Esta situación ocurre cuando ¦ es menor que 1.0.
j K j
Por lo tanto, la línea continua, la cual representa condiciones en la región de dos fases,
se presenta solo si ambas cantidades son mayores que 1.0. Un simple chequeo para
asegurar que las condiciones están en la región de dos fases es:
¦ z j K j ! 1.0 6.19
j
zj
¦ ! 1.0 6.20
j Kj
Este chequeo debe realizarse antes de efectuar cualquier cálculo de equilibrio gas -
líquido.
E jemplo 6.5. Realice el chequeo que debería haberse desarrollado antes de hacer los
cálculos en el ejemplo 6.1.
Solución:
¦ z j K j ! 1.0
j
zj
¦ ! 1.0
j Kj
Ambas sumas son mayores de 1.0, así que la mezcla está en la región de dos fases a 150
ºF y 200 psia.
Los significados de las curvas a trazos son similares a los de las curvas a trazos de la
Figura 6.1. La línea continua representa la solución al ejemplo 6.4, en el cual un valor
correcto de 0.547 moles es obtenido para nL.
¦zjK j 1 6.21
s
ga
todo
,
Td
T>
<P
dó ses
fa
P dos
< Td,
b <T
by T
d<P<P
P
P q
í
l
ouo
d
i
>
Pó
b<
Tt
o
,
b
Td
Valor de prueba de ñL
F igura 6.2. Comportamiento de la ecuación 6.18 con valores de prueba de ñL
incorrectos y cálculos a presiones o temperaturas por encima de la región de
dos fases
Un valor aproximado de presión de punto de burbuja para el uso como valor inicial de
prueba, puede ser determinado usando la ecuación 6.11. Cuando la presión correcta de
punto de burbuja ha sido establecida por este proceso de prueba y error, la composición
del líquido en el punto de burbuja está dada por cada uno de los términos de la
sumatoria.
Solución:
¦zjK j 1
Segundo, trate con un nuevo valor de prueba de presión de punto de burbuja de 220
psia. Determine valores de K j a 150 ºF y 220 psia.
Tercero, graficar los valores de prueba contra la suma resultante como en la figura 6.3,
e interpolar para determinar un tercer valor de prueba de presión de punto de burbuja de
237 psia.
Valor resultante de 6 zj K j
En el punto de rocío, la mezcla está enteramente en la fase gas excepto por una cantidad
infinitesimal de líquido. Por lo tanto nL = 0 y ng = n así que la ecuación 6.15 llega a ser:
zj
¦ 1 6.22
j Kj
E jemplo 6.7. Estime la presión de punto de rocío de la mezcla dada en el ejemplo 6.1 a
una temperatura de 150 ºF. Asuma que las cartas de relación de equilibrio dadas en el
capítulo 25 del libro de la GPSA pueden usarse para esta mezcla.
Solución:
zj
¦ 1
j Kj
Para ahorrar espacio, sólo se muestra el último cálculo con un valor final de
prueba de
Pd = 141 psia.
Determinar los valores de K j a 141 psia y 150 ºF a partir de las cartas de relación de
equilibrio dadas en el capítulo 25 del libro de la GPSA
zj
¦ 1 , Entonces Pd = 141 psia
j Kj
6.3 V A PO R I Z A C I Ó N F L ASH
Los cálculos del equilibrio gas- liquido usando ya sea la ecuación 6.17 ó 6.18, son a
menudo llamados cálculos de vaporización flash.
Los resultados de los cálculos flash están en términos de libra-mol. Las composiciones
del gas y del liquido en equilibrio, son usadas para calcular los pesos moleculares
y densidades.
6.4 V A P O R I Z A C I Ó N D I F E R E N C I A L
La figura 6.4A muestra nL moles de una mezcla líquida contenida en una celda a
temperatura constante. La presión se coloca en la presión de punto de burbuja del
liquido. Cuando una cantidad de mercurio se remueve de la celda (como se muestra en
la figura 6.4B), la presión se reduce ligeramente, y dnL moles de liquido se vaporizan.
El mercurio es luego reinyectado a la celda, y la presión se mantiene constante
permitiendo que el gas se escape a través de una válvula.
El procedimiento anterior deja la celda llena de un líquido, el cual tiene una presión de
punto de burbuja igual a la presión en la celda. Ahora, solamente quedan nL - dnL moles
de líquido, puesto que se han removido dnL moles del gas. La repetición de este proceso
trae como resultado la vaporización diferencial del líquido.
Además, xjnL representa las moles del j-ésimo componente en el líquido al comienzo de
la vaporización diferencial,
(xj - dxj)(nL - dnL) representa las moles del j-ésimo componente en el liquido al final del
incremento, y yjdnL representa las moles del j-ésimo componente en el gas removido.
El número de moles del j-ésimo componente en el gas, se puede igualar con el número
de moles del j-ésimo componente que se pierde en el líquido. Entonces:
y j dn j x j n j x j dx j n L dn L 6.23
y j dn j x j n L x j n L n L dx j x j dn L dx j dn L 6.24
Puesto que dxj es pequeña, el término dxjdnL se puede despreciar. También K jxj se
puede sustituir por yj para obtener:
dn L 1 §¨ dx j ·
¸ 6.26
nL K j 1 ¨© x j ¸
¹
n Lf 1 x jf
ln ln , 6.28
n Li K j 1 x ji
x jf § nLf ·
K j 1 6.29
¨¨ ¸¸
x ji © nLi ¹
Puesto que las ecuaciones 6.28 ó 6.29 involucran dos incógnitas, xjf y nLf, se requiere
otra ecuación para completar el cálculo. La ecuación 6.30 es la apropiada.
¦ x jf 1.0 6.30
j
Las composiciones calculadas del líquido se suman como se indica en la ecuación 6.30.
Si la sumatoria es igual a 1.0, el valor de prueba de nLf, es la cantidad correcta de moles
de líquido que queda después de la vaporización llevada a la presión final. Si la
sumatoria de fracciones molares de los componentes del líquido no son iguales a 1.0,
se debe seleccionar otro valor de prueba de nLf y repetirse el procedimiento hasta que la
sumatoria xjf sea igual a 1.0.
Valor resultante de 6 xjf
Note que un valor de prueba de nLf / nLi igual a 1.0 resulta en ¦ x jf 1.0 ,
j
prescindiendo de los valores de K j. Esta es una solución trivial. También, a medida que
el valor de nLf / nLi se aproxima a cero, ¦ x jf llega a ser muy grande. Por lo tanto la
j
curva que relaciona los valores de prueba de nLf/nLi, a los valores que resultan de ¦ x jf
j
toma la forma que se muestra en la figura 6.5. Un análisis del comportamiento de esta
curva, ayuda a la selección de valores de prueba de nLf/nLi.
Componente Composición,
fracción molar
Metano 0.160
Propano 0.380
n-Pentano 0.460
1.000
Solución:
Determine por prueba y error, un valor de nLf/nLi, que satisfaga las ecuaciones
K j 1
x jf § n Lf ·
¨¨ ¸¸ 6.29
x ji © n Li ¹
¦ x jf 1.0 6.30
j
Comience con un valor de nLf/nLi = 0.90. Use la forma de la figura 6.5 como una guía
para valores sucesivos de las pruebas.
Puesto que, nLf/nLi = 0.93 es el valor correcto. Por lo tanto el 7% en mol de liquido se
ha vaporizado, y la composición del líquido que queda, está dada en la última columna
de la tabla anterior.
La sumatoria de xjf debe ser igual a 1.0; si no, se debe seleccionar un nuevo valor de
prueba de la presión final y repetirse el cálculo.
Los resultados de este cálculo de prueba y error no tienen una curva característica tal
como la observada previamente en la Figura 6.5. Sin embargo, la relación entre ¦ x jf y
j
Solución:
Determine por prueba y error un valor de P f el cual dará valores de K j que satisfaga las
ecuaciones
K j 1
x jf § n Lf ·
¨¨ ¸¸
x ji © n Li ¹
¦ x jf 1.0
j
Comience con un valor de prueba de P f =300 psia, es decir, Ppromedio = 450 psia.
La suposición clave en la derivación de estas ecuaciones que se usan para predecir los
resultados de la vaporización diferencial se da en la ecuación 6.27; K j permanece
constante para cambios razonables de presión. Las relaciones de equilibrio no
permanecen constantes con cambios de presión.
6.5 R E L A C I O N D E E Q U I L I B R I O - C O R R E L A C I O N ES
En esta sección se examinaran las correlaciones del factor K . La ecuación que define el
factor K es:
yj
Kj 6.14
xj
yj Pvj
6.2
xj P
Esta ecuación puede ser combinada, obteniéndose una ecuación para calcular el factor K
de un componente, j, en una mezcla la cual se comporta como una solución ideal.
Pvj
Kj 6.31
P
Por lo tanto, el factor K para un componente de una solución real no depende solamente
de la presión y la temperatura sino también del tipo y cantidades de otras sustancias
presentes. Esto significa que cualquier correlación del factor K debe estar basada en al
menos tres propiedades: presión, temperatura y una tercera propiedad, la cual describe
el comportamiento de una solución no ideal. Esta propiedad debe representar tanto los
tipos de moléculas presentes como sus cantidades en el gas y en el líquido.
Para definir esta tercera propiedad se han realizado varios tipos de pruebas. Este texto
tratará solo la presión de convergencia, la cual parece ser la más conveniente para los
tipos de cálculos requeridos por los ingenieros de petróleos.
Un gran número de correlaciones del factor K han sido propuestas. Presentaremos sólo
dos: una gráfica y otra en forma de ecuación. Un conjunto parcial de curvas del factor K
es dado en el capítulo 25 del libro de la GPSA. Note que en la GPSA se presentan cartas
del factor K con varias presiones de convergencia. Un valor de presión de convergencia
aplicable al fluido de interés debe ser estimado a fin de seleccionar el conjunto correcto
de gráficas.
Las formas de las curvas en estas figuras son características de la mayoría de mezclas
multicomponentes. A bajas presiones, la pendiente de cada curva es aproximadamente -
1.0. Una pendiente de -1.0, para un factor K ideal es predecida por la ecuación 6.31.
Cada curva pasa a través del valor unidad a una presión muy cercana a la presión de
vapor del componente a la temperatura de la gráfica.
A elevadas presiones, se observan los efectos del comportamiento real. La curva para
cada componente parte de una pendiente de -1.0. Las curvas tienden a converger hacia
un factor K de 1.0.
Varios métodos han sido propuestos para estimar la presión de convergencia. Estos
métodos han sido evaluados usando datos de laboratorio de muestras de yacimientos
petroleros. El método descrito abajo es tan seguro como cualquiera de los otros y es el
más fácil de aplicar.
PK 60 M C 7 4,200 6.32
Las propiedades críticas del metano y las propiedades críticas promedio del componente
pesado hipotético se encuentran en la figura 6.9; y la ubicación de la presión de
convergencia para este par, es interpolada usando el domo crítico adyacente como guía.
E jemplo 6.10.
Las composiciones del gas y del liquido de un aceite negro a 1,300 psia y 160 ºF
calculado con la ecuación 6.17 son dadas en la siguiente tabla. Una presión de
convergencia de 5,000 psia fue usada para obtener el factor K . Qué valor de presión de
convergencia debería usarse para esta mezcla a 160 ºF?
Solución:
Tercero, grafique el punto crítico del peso promedio, en la figura 6.9, interpolando un
domo de presión de convergencia, y lea P K a 160 ºF.
Una presión de convergencia de 10,000 psia debería haber sido usada para esta mezcla.
Ver figura 6.10. Repita el cálculo de equilibrio gas - líquido, utilizando un factor
K a P K = 10,000.
Existen cartas para la presión de convergencia de 800, 1,000, 1,500, 2,000, 3,000, 5,000
y 10,000 psia. Cuando la presión de convergencia para la mezcla, está entre los valores
para los cuales las cartas son dadas interpole entre las cartas. La interpolación es
necesaria cuando la presión de operación está cerca de la presión de convergencia. A
bajas presiones, simplemente se usa la carta con presión de convergencia más cercanas
al valor de la mezcla.
La figura 6.11, muestra los resultados de los cálculos para un aceite negro a condiciones
de yacimiento. En esta situación la relación de equilibrio para el metano es obviamente
la única que debe ser conocida con extrema exactitud.
La figura 6.12, muestra los resultados de los cálculos para un aceite negro a condiciones
de superficie. En esta situación, las relaciones de equilibrio para todos los
componentes más livianos contribuyen a errores en el volumen del líquido calculado.
Así, un error sustancial en la relación de equilibrio es requerido para producir un error
significativo en el volumen el líquido calculado.
La figura 6.13 muestra los resultados de los cálculos para un gas retrógrado a
condiciones de yacimiento. Inexactitudes en la relación de equilibrio, tanto para el
componente más liviano, como para el componente más pesado representado por el
heptano plus, da como resultado un drástico error en el volumen del líquido calculado.
De hecho un error de solo el 8% en la relación de equilibrio para el metano resulta en el
cálculo de un punto de rocío en vez del volumen de liquido correcto. Esto indica que
los cálculos cercanos al punto de rocío de una mezcla dependen altamente de la
exactitud con la cual la relación de equilibrio es conocida.
calculado
F 6*1'.)2(7.!8.!.11*1!.)!&(+!'*)+2()2.+!8.!.93:&:;1:*!
igura 6.11. Efectos de los errores en las relaciones de equilibrio sobre el
cálculo de la cantidad de líquido de un petróleo negro a condiciones de
yacimiento.
La mejor fuente de factores K , son los datos de laboratorio de una misma mezcla a de
una mezcla similar de petróleo.
§ 1 1·
bj ¨ ¸ 6.33
¨T ¸
© Bj T ¹
PK, psia
b[1/TB - 1/T]
T abla 6.1
V alores de bj
H idrocarburo V alores A tenuados
Metano 805
Etano 1412
Propano 1799
i-Butano 2037
n-Butano 2153
i-Pentano 2368
n-Pentano 2480
n-Hexano 2780
n-Heptano 3068
n-Octano 3335
n-Nonano 3590
n-Decano 3828
n-Undecano 4055
n-Dodecano 4291
n-Tridecano 4500
n-Tetradecano 4715
n-Pentadecano 4919
n-Hexadecano 5105
n-Heptadecano 5290
n-Octadecano 5470
n-Nonadecano 5630
n-Eicosano 5790
n-Heneicosano 5945
n-Docosano 6095
Los datos del factor K son atenuados graficando una línea recta de los datos y moviendo
verticalmente cada punto hasta la línea. La línea recta puede no ser un "buen ajuste"
lineal en un sentido estadístico. Uno o dos puntos dados pueden estar obviamente en
error deberían ser ignorados al dibujar la línea.
E jemplo 6.11: Los datos del factor K para el metano mediante la línea 2 del hexano a
2,900 psig, de la figura 6.14 son dados a continuación. Atenuar los datos.
Componente F actor K
Metano 1.97
Etano 0.956
Propano 0.545
i-Butano 0.355
n-Butano 0.376
i-Pentano 0.178
n-Pentano 0.292
Hexano 0.121
Solución:
Primero calcule:
§ 1 1·
PK j y bj ¨ ¸
¨T T ¸¹
© Bj
Componente Psia bj ºR § 1 1·
P Kj T Bj bj ¨ ¸
¨T ¸
© Bj T ¹
C1 5742 805 201.0 2.79
C2 2786 1412 332.2 2.11
C3 1589 1799 416.0 1.60
i-C4 1035 2037 470.5 1.25
n-C4 1096 2153 490.8 1.13
i-C5 519 2368 541.8 0.79
n-C5 851 2480 556.6 0.70
C6 353 2780 615.4 0.31
> @
Segundo. grafique P Kj contra b j 1 TBj 1 T dibuje una línea con buen ajuste. y
ajuste los puntos dados a la línea. Ver la figura 6.15.
6000 C1
C2
PK, psia
C3
n-C4
1000 i-C4
n-C5
i-C5
300 C6
b/(1/TB-1/T)
F igura 6.15. PK contra b/(1/TB-1/T) (parte de la solución
del ejemplo 6.11
Componente P Kj K j.
A justado psia A justado
C1 5844 2.00
C2 2797 0.960
C3 1610 0.552
i-C4 1102 0.378
n-C4 968 0.332
i-C5 669 0.230
n-C5 607 0.208
C6 398 0.137
La pendiente y los interceptos de las líneas rectas trazadas por los datos del factor K a
varias presiones diferentes. tales como las líneas 2 a 7 de la figura 6.14 pueden ser
linealmente correlacionados contra la presión. Esto permite la interpolación y alguna
extrapolación de los datos. Por ejemplo el límite superior de los datos es una línea
horizontal donde P Kj. es igual a la presión de convergencia.
E jemplo 6.12. Los datos de pendientes e interceptos del conjunto de líneas rectas de la
2 a la 7 de la figura 6.14 son dados a continuación. Correlacionadas las pendientes y los
interceptos contra la presión es posible una interpolación y una extrapolación limitada
de los datos.
§ 1 1·
log pK j pendiente b j ¨¨ ¸ intercepto
¸ 6.35
© Bj T ¹
T
Solución:
Primero. Grafique la pendiente y el intercepto contra la presión. Ver figuras 6.16 y 6.17.
La forma del gráfico del factor ecuación 6.33. indica la posibilidad de extrapolación a
temperaturas diferentes a la temperatura a la cual los datos fueron obtenidos.
6.)8:.)2.!
Presión, psia
F igura 6.16. Pendientes de P K contra presión (parte
de la solución del ejercicio 6.12.
>)2.1'.02*!!
61.+:<)=!0+:(!
F igura 6.17. Interceptos de PK contra presión (parte de la
solución del ejemplo 6.12.)
6.6 E Q U I L I B R I O G AS L Í Q U I D O C O N E C U A C I O N ES D E EST A D O
Las ecuaciones de estado pueden ser utilizadas para calcular el equilibrio gas ± líquido
como una alternativa al uso de las correlaciones para el factor K. Se puede asumir que
las ecuaciones presentadas en el capítulo 3 predicen las relaciones presión ± volumen ±
temperatura para líquidos como para gases.
Una ecuación de estado de la forma de Van der Waals podrá producir isotermas como
las mostradas en la figura 6.18 para una sustancia pura. Nótese que las isotermas por
encima y en el punto crítico se parecen mucho más a las isotermas experimentales
correspondientes a la figura 5.5 correspondiente al diagrama presión ± volumen para el
etano
El rizo de Van der Waals es usado para determinar los volúmenes molares del equilibrio
gas y líquido y luego es reemplazado por una línea conectora entre esos dos volúmenes.
La base de este cálculo está en que el potencial químico del gas en equilibrio es igual al
potencial químico del líquido en equilibrio.
Presión, P
Volumen molar VM
dG V M dP 6.38
La figura 6.19A muestra una isoterma calculada por debajo de la temperatura crítica. El
punto a es seleccionado arbitrariamente a lo largo de la parte líquida de la isoterma.
Los puntos b y c son puntos máximos y mínimos en la curva de Van der Waals. Los
puntos e a lo largo de la curva, representan los puntos para los cuales el potencial
químico es igual. El punto f es el punto en la línea bc, el cual tiene la misma presión
que los puntos e.
A B a
a
Presión, P
Presión, P
c e c
e e d
d
f f
A B
Presión, P
Presión, P
0 Volumen molar, VM 0 Potencial químico, G
El cálculo de las propiedades del equilibrio gas ± líquido se resuelve con el cálculo de
los potenciales químicos de las dos fases.
Se requiere una ecuación de estado para evaluar de la ecuación 6.38. Para un gas ideal
se tiene:
RT
VM
P
dG RT P dP RTd ln P 6.39
!
!
!
Presión, P
!
!
!
! 0 Volumen molar, VM
!
! F igura 6.21. Valores de volumen molar, calculados con una ecuación de
estado de dos constantes, "el rizo de Van der Waals" ha sido removido de
! la isoterma
6.6.1.2 L a fugacidad
dG RT d ln f 6.40
lim f P 6.41
P o0
Nótese que la fugacidad simplemente reemplaza la presión en la ecuación del gas ideal
y la convierte en una ecuación de un gas real. La fugacidad tiene unidades de presión.
La ecuación 6.41 simplemente establece que a bajas presiones el fluido se comporta
como un fluido ideal.
Combinando las ecuaciones 6.38 con la 6.40 se define la ecuación para la fugacidad de
una substancia pura.
d ln f V M RT dP 6.42
Recuérdese que el potencial químico para el líquido debe ser igual al potencial químico
del gas en el equilibrio. Para una substancia pura, esto significa que en cualquier punto
a lo largo de la línea de presión de vapor, el potencial químico del líquido será igual al
del gas.
De esta forma la ecuación 6.40 muestra que la fugacidad del líquido debe ser igual a la
fugacidad del gas en equilibrio sobre la línea de presión de vapor. Así el equilibrio
puede ser calculado bajo la siguiente condición:
fg fL 6.43
La siguiente expresión puede ser derivada a partir de la ecuación 6.42 bajo la restricción
de la ecuación 6.41
f 1 V M § RT ·
ln Z 1 ln Z ³ ¨¨ P ¸¸ dV M 6.44
P RT f © V M ¹
donde Z es el factor de compresibilidad como está definido por las ecuaciones 2.38 y
2.39.
Para una substancia pura, la razón de fugacidad a presión, f/P, es llamada coeficiente de
fugacidad.
La ecuación de estado de Peng ± Robinson, ecuación 3.57 podrá ser utilizada con la
ecuación 6.44 con el fin de desarrollar un procedimiento para el cálculo de la presión de
vapor de una sustancia pura. La presión de vapor es simplemente la presión en los
puntos e de la figura 6.20, para los cuales la fugacidad del líquido es igual a la
fugacidad del gas.
ª aT º
«P V M b RT 3.57
¬ V M V M b bV M b »¼
§ RT · § a 2bRT · §a ·
V M3 ¨ b ¸ V M2 ¨ 3b 2 ¸ V M b¨ b b 2 ¸ 0 6.45
© P ¹ ©P P ¹ ©P ¹
RT C
b 0.07780 6.46
PC
y
R 2 TC2
aC 0.45724 6.47
PC
aT a CD 6.48
D
1
2
©
1 0.37464 1.54226Z 0.26992Z 2 §¨1 Tr 2 ·¸
1
¹
6.49
El factor acéntrico, Z, es una constante. Cada substancia pura tiene un valor de factor
acéntrico.
A los coeficientes en las ecuaciones 6.46 y 6.47 también les son asignados símbolos
como sigue:
VM ZRT P
Z 3 1 B Z 2 A 2 B 3B 2 Z AB B 2 B 3 0 6.51
Donde:
aT P
A 6.52
R 2T 2
bP
B 6.53
RT
La ecuación 6.51 es una ecuación cúbica con coeficientes reales. Así, tres valores del
factor Z hacen que la ecuación sea igual a cero. Esas tres raíces son todas reales
cuando la presión y la temperatura están en la línea de presión de vapor, que es cuando,
el líquido y el gas están presentes. Una raíz real y dos raíces complejas existen cuando
la temperatura está por encima de la temperatura crítica.
La figura 6.22 presenta la forma de una isoterma calculada con la ecuación 6.51 a la
temperatura por debajo de la temperatura crítica. Los puntos aebf son equivalentes a
los puntos aebf de las figuras 6.19 y 6.20. Los puntos e son valores del factor Z que
pueden ser medidos experimentalmente. El punto f no es una solución física.
Como antes, la curva ebfce es eliminada. Ver la figura 6.23. El punto e superior es el
factor Z del gas en equilibrio y el punto e inferior es el factor Z del líquido en equilibrio.
La línea punteada que conecta ambos puntos no tiene un significado físico. La curva a
trazos representa la envolvente de fases completa. Nótese la similitud de la isoterma de
la figura 6.23 con la isoterma experimental a 104ºF de la figura 2.4.
Así, la ecuación 6.51 es resuelta para sus tres raíces. Si se tiene sólo una raíz real, la
temperatura está por encima de la temperatura crítica. Si las tres raíces son reales, la
mayor corresponde al factor Z de equilibrio del gas y la menor al factor Z de equilibrio
del líquido. Estos son los puntos e en la figura 6.22. La raíz del medio, representada por
el punto f, es descartada.
La ecuación 6.51 puede ser usada para completar la integración de la ecuación 6.44
obteniendo como resultado:
ª Z § 2 12 1· B º
§ fg · A « g ¨© ¸ »
ln ¨¨ ¸¸ Z g 1 ln Z g B 1.5 ln « ¹
» 6.54
© P ¹ 2 B « Z §¨ 2 2 1·¸ B »
1
g
¬ © ¹ ¼
y
ª Z § 2 12 1· B º
§f · A « L ¨© ¸ »
¹
ln ¨ L ¸ Z L 1 ln Z L B 1.5 ln « » 6.54
© P¹ 2 B « Z §¨ 2 2 1·¸ B »
1
L
¬ © ¹ ¼
La ecuación 6.54 se puede aplicar dos veces: una con el factor Z del líquido para
calcular la fugacidad del líquido y otra con el factor z del gas a para calcular la
fugacidad del gas.
El procedimiento para calcular la presión de vapor de una substancia pura involucra las
ecuaciones 6.46 a la 6.54. Una vez que se selecciona una temperatura, los resultados de
las ecuaciones 6.46 a la 6.54, son fijos. El problema consiste en encontrar una presión
para utilizarla en las ecuaciones 6.51 a la 6.53 la cual dará los valores de factores Z para
los cuales puedan ser iguales los valores de las fugacidades del gas y el líquido.
Factor Z
Presión, P
Punto crítico
Factor Z
Presión, P
F igura 6.23. Valores del factor Z a una temperatura por debajo de la temperatura crítica,
calculados con una ecuación de estado de dos constantes, "el rizo de Van der Waals" ha
sido removido.
E jemplo 6.13 Calcular la presión de vapor del iso ± butano a 190 ºF usando la ecuación
de estado de Peng ± Robinson.
Solución
T 190 459.67
Tr 0.88495
Tc 274.46 459.67
b :b
RT c
0.07780 10.732734.13 1.1611
Pc 527.9
aC :a
R 2 Tc2
0.45724 10.732 734.13
2 2
53,765
Pc 527.9
D
1
2
©
1
1 0.37464 1.54226Z 0.26992Z 2 §¨1 Tr 2 ·¸
¹
Paso 2. Por prueba y error hallar la presión que ocasiona que la fugacidad del líquido
calculada con las ecuaciones 6.51, 6.52, 6.53 y 6.54 sea igual a la fugacidad del gas
calculada con las mismas ecuaciones. Únicamente el cálculo con un valor final de
prueba de P = 228.79 psia se mostrará.
A
aT P 57995228.79 0.27295
R 2T 2 10.7322 649.672
B
bP 1.1611228.79 0.038101
RT 10.732649.67
Las tres raíces se obtienen de,
Z 3 1 B Z 2 A 2 B 3B 2 Z AB B 2 B 3 0
ZL 0.067258 y Z g 0.70786
Entonces
ª Z § 2 12 1· B º
§ fg · A « g ¨© ¸ »
ln ¨¨ ¸¸
Z g 1 ln Z g B 1.5 ln « ¹
»
© P ¹ 2 B « Z §¨ 2 2 1·¸ B »
1
g
¬ © ¹ ¼
y
fg 176.79 psia
ª Z § 2 12 1· B º
§f · A « L ¨© ¸ »
¹
ln ¨ L ¸ Z L 1 ln Z L B 1.5 ln « »
© P¹ 2 B « Z §¨ 2 2 1·¸ B »
1
L
¬ © ¹ ¼
fL 176.79 psia
f L = fg, de esta forma el valor de prueba de P, 228.79 psia, es la presión de vapor del
isobutano a 190 ºF.
Nótese que el valor de los volúmenes molares puede ser calculado fácilmente.
V Mg
Z g RT 0.7078610.732649.67 21.57 ft 3 lb mol
P 228.79
V ML
Z L RT 0.06725810.732649.67 2.050 ft 3 lb mol
P 228.79
Compare los resultados con los datos experimentales.
C alculado E xperimental
Presión de vapor, psia 228.8 228.3
Volumen molar del líquido saturado, ft3/lb mol 2.050 2.035
Volumen molar del vapor saturado, ft3/lb mol 21.57 21.68
Factor- Z del líquido 0.0673 0.0666
Factor- Z del gas 0.7078 0.7101
Esos métodos pueden ser agrupados en dos tipos generales: sustituciones sucesivas y
Newton ± Raphson. Estas técnicas no son discutidas en este texto.
6.6.2 Mezclas
La situación concerniente a las mezclas algunas veces es más difícil de visualizar. Sin
embargo, el equilibrio se obtiene cuando el potencial de cada componente en el líquido
es igual al potencial químico de cada componente en el gas.
dG j RT d ln f j 6.55
Esto significa que cuando la presión tiende a cero, el fluido se acerca al comportamiento
ideal y la fugacidad de un componente se acerca a la presión parcial de dicho
componente.
fg j fL j 6.57
Los valores de fugacidad para cada componente son calculados con una ecuación de
estado. Cualquier ecuación de estado puede ser usada para dichos cálculos. Más
adelante desarrollaremos las ecuaciones en un ejemplo, utilizando la ecuación de Peng ±
Robinson.
fj
Ij 6.58
yjP
ª º
1 V « RT §¨ wP ·
¸ » dV ln Z
ln I j ³ 6.59
RT f « V ¨© wn j ¸
¹ T ,V ,n j »¼
¬
f Lj
I Lj xjP yj
Kj 6.60
I gj f gj xj
yjP
RT aT
P
V M b V M V M b b V M b
b ¦ y jbj
j
y
¦ ¦ y i y j a T i a T j 2 1 k i j
1
aT
i j
k ii k jj 0 6.61
ki j k ji 6.62
Los valores de los coeficientes para los componentes individuales son calculados de
esta forma:
RT c j
bj 0.07780 6.46
Pc j
y
aT j ac j D j 6.48
Donde:
R 2 Tc2j
ac j 0.45724 6.47
Pc j
D
1
2
©
1 0.37464 1.54226Z j 0.26992Z 2j §¨1 Trj 2 ·¸
1
¹
6.49
z 3 1 B z 2 A 2 B 3B 2 z AB B 2 B 3 0 6.51
Donde:
aT P
A 6.52
R 2T 2
bP
B 6.53
RT
Cuando se obtienen tres raíces, la menor corresponde al factor Z del líquido, la mayor al
factor Z del gas y la del medio es descartada. Esto es análogo a la eliminación del punto
f en la figura 6.22
.
Combinando la ecuación de estado de Peng ± Robinson y la ecuación 6.59 se obtiene
una ecuación para el coeficiente de fugacidad de cada componente.
ª Z § 2 12 1· B º
A « ¨ ¸ »
ln I j ln Z B Z 1 B ' j 1.5 A' j B ' j ln « © ¹
» 6.63
2 B « Z §¨ 2 2 1·¸ B »
1
¬ © ¹ ¼
Donde:
bj
B' j 6.64
b
2 a Tj ¦ y i a Ti2 1 k ij º
1 ª 12 1
A' j 6.65
a T «¬ i »¼
Las ecuaciones de la 6.63 a la 6.65 están escritas dos veces: una para valores de Z , b y
a T del gas y otra para valores de Z , b y a T del líquido. El procedimiento para el cálculo
del equilibrio gas ± líquido a una temperatura y presión dada se da a continuación.
Los valores a Tj y bj para cada componente de la mezcla se obtienen con las ecuaciones
6.46 a 6.49 a partir de las propiedades críticas conocidas y de los factores acéntricos de
los componentes puros.
Se obtiene una primera serie de valores de prueba del factor K . Por ejemplo, la
siguiente ecuación para el factor K puede utilizarse como primera estimación.
Luego estos valores de K j son usados en los cálculos de equilibrio gas ± líquido, como
se escribe en la sección 6.2, para determinar las composiciones del gas y el líquido. Las
restantes ecuaciones se resuelven tanto para el líquido como para el gas.
También se debe calcular B´j y A´j para cada componente j. B´jL se obtiene cuando se
usa bL en la ecuación 6.64 y B´jg si se utiliza bg. A´jL y A´jg se obtienen de igual forma de
la ecuación 6.65.
La ecuación 6.63 se resuelve para los coeficientes de fugacidad de los componentes del
líquido, ) Lj, usando valores de A L , B L , zL , A´jL y B´jL .
Los valores de )gi resultan cuando se usan en la ecuación 6.63 el factor Z y los
coeficientes correspondientes al gas.
Posteriormente los valores de la fugacidad del líquido y el gas para cada componente se
conocen a partir de:
fL j x j PI L j 6.58
fg j y j PI g j
6.58
fg j fL j 6.57
j fL j fg j 6.66
¦ j
2
6.67
j
j
K T
j K Cj
2
6.68
K Tj K Cj
Donde K j C son los factores K calculados y K j T son los valores de prueba de los factores
K.
La convergencia a la solución correcta se obtiene cuando la suma de las funciones error
es menor que alguna tolerancia seleccionada. Si la suma de las funciones error es
mayor que la tolerancia, el factor K j C se usa como nuevo valor de prueba de K j y se
repite el proceso.
E jemplo 6.14. Calcular las composiciones del gas y el líquido cuando la mezcla con la
composición dada es llevada al equilibrio a 160 ºF y 1,000 psia. Utilice la ecuación de
estado de Peng ± Robinson. Use los coeficientes de interacción binaria de 0.02 para el
metano - n-butano, 0.04 para metano - n-decano, y 0.0 para n-butano - n-decano.
Componente Composición,
F racción molar
Metano 0.5301
n-Butano 0.1055
n-Decano 0.3644
1.0000
Solución:
Paso 1. Calcule los coeficientes de los componentes de la mezcla.
Dj
1
2 ©
1 0.37464 1.54226Z j 0.26992Z 2j §¨1 Trj 2 ·¸
1
R 2 Tc2j
ac j 0.45724
Pc j
aT j a c jD j
RT c j
bj 0.07780
Pc j
Componente T cj P cj Zcj Dj a cj a Tj bj
°R psia
C1 343.0 666.4 0.0104 0.7481 9,297 6,956 0.4297
n-C4 765.3 550.6 0.1995 0.1394 56,017 63,827 1.1604
n-C10 1111.7 305.2 0.4898 0.6139 213,240 344,149 3.0411
Paso 2. Seleccione valores de prueba para los factores K y calcule las composiciones de
prueba del gas y líquido en equilibrio. Únicamente la prueba final, con los valores de K
obtenidos se muestra a continuación:
zj
¦xj ¦
j j 1 ñ g K j 1
Este cálculo requiere prueba y error; solamente la prueba final con ñg = 0.4025 es
mostrada.
Componente Kj zj xj yj
C1 3.992 0.5301 0.2408 0.9613
n-C4 0.2413 0.1055 0.1517 0.0366
n-C10 0.00340 0.3644 0.6075 0.0021
1.0000 1.0000 1.0000
¦ ¦ y i y j a T i a T j 2 1 k i j
1
aT
i j
b ¦ y jbj
j
aT P
A
R 2T 2
bP
B
RT
F ase aT b A B
Líquido 171,446 2.1270 3.8766 0.31983
Gas 8,177 0.4619 0.1849 0.06945
Z 3 1 B Z 2 A 2 B 3B 2 Z AB B 2 B 3 0 6.50
ZL 0.3922 y Z g 0.9051
A' j
1 ª2 a 12 ¦ y a 12 1 k º
aT «¬ Tj i i T i ij »
¼
bj
B' j
b
L íquido G as
Paso 6. Calcular los coeficientes de fugacidad de los componentes del líquido y del gas.
ª Z § 2 12 1· B º
« ¨ ¸ »
ln Z B Z 1 B ' j 1.5 A' j B ' j ln «
A © ¹
ln I j »
2 B « Z §¨ 2 2 1·¸ B »
1
¬ © ¹ ¼
Componente I Lj Igj
C1 3.67552 0.92065
n-C4 0.12878 0.53373
n-C10 0.000699 0.20600
IL j
Kj
Ig j
j
K T
j K Cj
2
K Tj K Cj
Componente K jC j
C1 3.992 0.000
n-C4 0.2413 0.000
n-C10 0.00340 0.000
¦ j 0.000
j
La suma de las funciones error es menor que una tolerancia de 0.001, de modo que la
serie de valores de prueba de los factores K fue correcta y los valores calculados de
composiciones de líquido y de gas son correctos.
C alculado E xperimental
Componente xj yj xj yj
ZRT
VM
P
PM
U
ZRT
F ase VM Ma U
Z ft3/lb mol lb/lb-mol lb/ft3
Líquido 0.3922 2.61 99.12 38.00
Gas 0.9051 6.02 17.84 2.96