Monografía Factor de Compresibilidad Z PDF

UNIVERSIDAD DE ORIENTE
NÚCLEO DE MONAGAS
ESCUELA DE INGENIERÍA Y CIENCIAS APLICADAS
DEPARTAMENTO DE INGENIERÍA DE PETRÓLEO
CURSOS ESPECIALES DE GRADO
INGENIERÍA DE GAS NATURAL I-2018
ANÁLISIS COMPARATIVO DE LAS DIFERENTES METODOLOGÍAS

PARA LA DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD (Z)
DEL GAS APLICADO A MUESTRAS DE GAS NATURAL
REALIZADO POR:
YÉRIKA GLECMAR CARRERA MARCANO
JAVIER EDUARDO CASTILLO LA ROSA
FÉLIX DAVID ROQUE PÉREZ
Seminario de Área de Grado presentado como requisito parcial para

optar al Título de:
INGENIERO DE PETRÓLEO
MATURÍN, MAYO DE 2018

NÚCLEO DE MONAGAS

REALIZADO POR:
YÉRIKA GLECMAR CARRERA MARCANO, C.I. 21.051.876
JAVIER EDUARDO CASTILLO LA ROSA, C.I. 19.403.248
FÉLIX DAVID ROQUE PÉREZ, C.I. 19.804.888
REVISADO POR:
__________________________
Msc. Tomás Marín
ii
NÚCLEO DE MONAGAS

REALIZADO POR:
YÉRIKA GLECMAR CARRERA MARCANO, C.I. 21.051.876
JAVIER EDUARDO CASTILLO LA ROSA, C.I. 19.403.248
FÉLIX DAVID ROQUE PÉREZ, C.I. 19.804.888
APROBADO POR:
____________________
Msc. Tomás Marín
Asesor Académico
____________________ ____________________
Msc. Pedro Tineo Msc. Rubén Vega
Jurado Principal ____________________ Jurado Principal
Ing. Miguel Flores
Jurado Principal
iii
RESOLUCIÓN
DE ACUERDO AL ARTÍCULO 41 DEL REGLAMENTO DE TRABAJOS DE

GRADO:
“LOS TRABAJOS DE GRADO SON PROPIEDAD EXCLUSIVA DE LA

UNIVERSIDAD DE ORIENTE Y SOLO PODRÁN SER UTILIZADOS A
OTROS FINES CON EL CONSENTIMIENTO DEL CONSEJO DE NÚCLEO
RESPECTIVO, QUIEN LO PARTICIPARÁ AL CONSEJO UNIVERSITARIO.”
iv
AGRADECIMIENTOS
A Dios primeramente por permitirme llegar hasta este punto de vida y

alcanzar algo que siempre vi tan lejano. Hoy por hoy, ya con esta meta
cumplida siento que valió la pena el esfuerzo.
A mi Madre, que desde el primer día me acompañó en este transitar

académico y que a pesar de no tenerla físicamente conmigo, se que siempre
me acompaña y orienta a cada paso del camino. A mi Padre, por su crianza y
orientación que me permitieron no caer en malos pasos y ser hoy una
persona de bien.
A mi hermana, mi concu, mi compañera de aventuras Virginia Roque,

por su paciencia, ayuda, guía, consejos, que me diste en la realización de
este trabajo, ten presente que eres parte importante de él.
A mis compañeros de clase, que hicieron más llevadero y ameno

este camino, ya que si todo fuera responsabilidades y academia, saliéramos
directo al psicólogo y no al ámbito empresarial. A Uds. gracias, Nore, Paty,
Gaby, Lesmari, Thalia, María, Alfredo, Yilson, Saúl, Cesar, espero la vida les
tenga preparada cosas buenas siempre.
A las Áreas de Gas, por permitirme cursar este último tramo que
faltaba para terminar mi carrera, así como también conocer buenas personas
y excelentes profesores. Mención especial a los de Tomás, team Z y team
Coco, fuimos los mejores muchachos. Finalmente, a la Universidad de
Oriente. Casa de estudios a la que llegué con un montón de expectativas y a
pesar de la relación amor y odio que tuvimos, hoy puedo decir: soy Ingeniero
de Petróleo. Gracias totales UDO Monagas.
Félix Roque
v
AGRADECIMIENTOS
A DIOS todopoderoso, por haberme dado salud y sabiduría para

alcanzar esta importante meta. Toda la gloria sea para ti padre creador.
Agradezco también la crianza, consejos y apoyo brindado por parte de

mis padres, sin ustedes esto no hubiese sido posible. Gran parte de este
logro les pertenece.
A mis hermanos, en especial a mi hermana Daicelia Zapata por

brindarme su hogar en todos estos años y por compartir sus conocimientos
sobre esta hermosa carrera.
A mis abuelos, tíos, primos y sobrinas, por siempre contar con ustedes
cuando los necesitaba.
A mis compañeros de estudio, por todos los momentos vividos dentro

y fuera de la Universidad, ustedes saben lo difícil que es este camino pero
con esfuerzo, empeño y dedicación lo podemos lograr.
Gracias a todas aquellas personas que en algún momento de mi vida

aparecieron como estrellas para iluminar mi camino, por los buenos o malos
consejos que me dieron, por darme momentos inolvidables tanto felices
como tristes, pero que al final son parte de la vida, a todos ellos por compartir
de su tiempo, de su vida, sus molestias, sufrimientos y logros.
Javier Castillo
vi
AGRADECIMIENTOS
A Dios por ser mi guía, por darme todo lo necesario para seguir
adelante y por regalarme la dicha de cumplir uno más de mis objetivos.
A mis padres, Carmen Marcano y Julio Carrera, por todos sus consejos,
por la educación que me han brindado, por creer en mí y estar presente en
los momentos más importantes de mi vida. ¡Gracias por tanto! Su inmenso
amor, esfuerzo y sacrificio han sido mi motor y fortaleza.
A mis hermanos Junior y Yunaika Carrera, en especial a ti mi negra, por

darme siempre lo mejor de ti y por acompañarme en todo momento.
A ti José Reyes por respetar mis decisiones, por alentarme siempre con
tus sabias palabras a luchar por alcanzar mis sueños, por tus buenos
deseos, por el amor que me das y todo lo que haces por mí.
A mis compañeros de este trabajo, especialmente a Félix Roque por

tanta dedicación. Todo el esfuerzo hoy se ve reflejado en nuestros logros.
A mis amigos, en especial a María Pereda, Yulimar Reyes y Angélica

Patiño. A mi suegra, familiares y todas aquellas personas que de alguna
manera ayudaron y contribuyeron al cumplimiento de esta meta.
Yérika Carrera
vii
DEDICATORIA
A mis Padres, especialmente a mi Madre Bethsy Pérez a quien

siempre tengo presente y espero siempre este orgullosa de mi donde quiera
que este. Te amo Ma.
A mi Abuela, Berta Roque, esa que me abrió las puertas de su casa

desde que llegue a Maturín y siempre estuvo allí con su bondad, atenciones
y regaños que me permitieron de una u otra forma alcanzar esta meta. Te
quiero Bertalina.
A mis hermanos, Carlos y Virginia, especialmente a Carlos espero

que la vida te permita alcanzar este logro también, ten presente que a pesar
de las dificultades siempre es posible. Cuentas con nosotros.
A la princesa Miranda Victoria, la sobrina que llegó a esta vida para

enamorarme y llenar de felicidad cualquier rincón que visite. Espero seas
siempre una señorita de bien princesa de ojos claros.
Félix Roque
viii
DEDICATORIA
Dedico este trabajo primeramente a Dios por permitirme llegar a este

momento tan importante en mi vida, por los triunfos y momentos difíciles que
me han enseñado a valorar la vida cada día más, por cuidar de mi y por
otorgarme esta familia tan hermosa que tengo.
A mis padres, que han sabido formarme con buenos principios y

valores, lo cual me ha me ayudado a enfrentar todas las adversidades que
se me han presentado en el camino.
A mis hermanos, por compartir alegrías y tropiezos de los cuales

supimos salir adelante, por su confianza y por permitirme estar en sus vidas.
A mi familia en general, por brindarme su apoyo incondicional y por

todos los momentos compartidos.
A mi grupo de amigos, que gracias al equipo que hicimos logramos

llegar al final de esta importante meta, espero seguir contando con ustedes
por el resto de mi vida.
A mis profesores, por su tiempo, apoyo y así como por los

conocimientos que me transmitieron en el desarrollo de mi formación
profesional.
Javier Castillo
ix
DEDICATORIA
A Dios Todopoderoso por todas las bendiciones que ha derramado

sobre mí.
A mis padres y hermanos por su infinito amor y por su apoyo

incondicional.
A ti José Reyes que te has convertido en pieza fundamental en mi vida

ocupando ese lugar especial en mi corazón. Llegaste en el momento
indicado, fortaleces mis ganas de continuar.
A mi tía Jenny Marcano, por brindarme parte de su tiempo, por sus

enseñanzas, su paciencia y por todo su cariño.
A mis abuelos por ser parte importante de mi vida y por haber

contribuido con su apoyo al cumplimiento de esta meta.
A ustedes y todas esas personas q creyeron en mí.
Yérika Carrera
x
ÍNDICE
Pág.
RESOLUCIÓN ............................................................................................... iv
AGRADECIMIENTOS .................................................................................... v
DEDICATORIA ............................................................................................ viii
LISTA DE TABLAS ..................................................................................... xiii
RESUMEN ................................................................................................... xiv
INTRODUCCIÓN ............................................................................................ 1
CAPÍTULO I.................................................................................................... 2
EL PROBLEMA Y SUS GENERALIDADES .................................................. 2
1.1. PLANTEAMIENTO DEL PROBLEMA ................................................ 2
1.2. OBJETIVOS DE LA INVESTIGACIÓN .............................................. 3
1.2.1 Objetivo general ....................................................................... 3
1.2.2 Objetivos específicos ............................................................... 3
1.3. JUSTIFICACIÓN DE LA INVESTIGACIÓN ....................................... 4
CAPÍTULO II................................................................................................... 5
MARCO TEÓRICO ......................................................................................... 5
2.1. ANTECEDENTES DE LA INVESTIGACIÓN..................................... 5
2.2. BASES TEÓRICAS ............................................................................ 6
2.2.1 Estado gaseoso ....................................................................... 6
2.2.2 Gas ideal .................................................................................. 7
2.2.3 Gas real ................................................................................. 13
2.2.4 Gas natural ............................................................................ 13
2.2.5 Clasificación del gas natural según su origen ........................ 14
2.2.6 Clasificación del gas natural según su composición .............. 14
2.2.7 Mezcla de Gases ................................................................... 16
2.2.8 Ecuaciones de estado ............................................................ 18
2.2.9 Factor de compresibilidad (Z) ................................................ 26
2.2.10 Significado del valor de Z ....................................................... 28
2.2.11 Determinación del factor de compresibilidad Z ...................... 28
2.2.12 Análisis no paramétrico de la varianza................................... 29
2.2.13 Prueba de Friedman .............................................................. 29
2.3. PROCEDIMIENTO METODOLÓGICO ............................................ 31
2.3.1 DESCRIPCIÓN DE LOS MÉTODOS MÁS UTILIZADOS PARA LA
DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z ..................... 31
xi
2.3.2 DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z A
CUATRO (4) MUESTRAS DE GAS NATURAL, UTILIZANDO NUEVE
(9) MÉTODOS REPRESENTATIVOS PARA SU CÁLCULO ................ 32
2.3.3 COMPARACIÓN A TRAVÉS DE UN ANÁLISIS ESTADÍSTICO
DE LAS DIFERENCIAS EXISTENTES ENTRE LOS RESULTADOS
OBTENIDOS ........................................................................................ 32
2.4. DEFINICIÓN DE TÉRMINOS BÁSICOS ......................................... 33
CAPÍTULO III................................................................................................ 35
DESARROLLO DE OBJETIVOS ................................................................. 35
3.1 DESCRIPCIÓN DE LOS MÉTODOS MÁS UTILIZADOS PARA LA
DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z ................. 35
3.1.1 Métodos basados en ecuaciones de estado para el cálculo del
factor de compresibilidad Z................................................................... 35
3.1.2 Métodos gráficos para el cálculo del factor de compresibilidad
Z 50
3.2 DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z A
CUATRO (4) MUESTRAS DE GAS NATURAL, UTILIZANDO NUEVE (9)
MÉTODOS REPRESENTATIVOS PARA SU CÁLCULO ......................... 57
3.2.1 Ecuación de Van Der Waals .................................................. 58
3.2.2 Ecuación de Peng Robinson .................................................. 59
3.2.3 Ecuación de Redlich-Kwong .................................................. 60
3.2.4 Ecuación de Soave-Redlich-Kwong ....................................... 61
3.2.5 Ecuación de Dranchuk Abou Kassen ..................................... 62
3.2.6 Ecuación de Dranchuk-Pulvis-Robinson ................................ 63
3.2.7 Método de Standing-Katz....................................................... 64
3.2.8 Método de Refracción Molecular de Eykman ......................... 65
3.2.9 Método de Pitzer .................................................................... 66
3.3. COMPARACIÓN A TRAVÉS DE UN ANÁLISIS ESTADÍSTICO DE
LAS DIFERENCIAS EXISTENTES ENTRE LOS RESULTADOS
OBTENIDOS............................................................................................. 67
3.3.1 Resumen estadístico.............................................................. 68
3.3.2 Cálculo para el estadístico de prueba .................................... 69
CAPÍTULO IV ............................................................................................... 72
CONCLUSIONES Y RECOMENDACIONES................................................ 72
4.1. Conclusiones ................................................................................... 72
4.2. Recomendaciones ........................................................................... 73
REFERENCIAS BIBLIOGRÁFICAS ............................................................ 74
APÉNDICES ................................................................................................. 77
xii
LISTA DE TABLAS
Pág.
Tabla 2.1–Valores de la constante general de los gases para diferentes
unidades de presión (P), volumen (V) y temperatura (T) .............................. 12
Tabla 2.2-Datos para la prueba de Friedman ............................................... 30
Tabla 3.1-Valores de los coeficientes Aij usados en la Ec. 3.32.................... 55
Tabla 3.2- Resultados del factor de compresibilidad Z utilizando la ecuación
de estado de Van Der Waals ........................................................................ 59
de estado de Peng Robinson ........................................................................ 60
de estado de Redlich-Kwong ........................................................................ 61
de estado de Soave-Redlich-Kwong ............................................................. 61
de estado de Dranchuk Abou Kassen ........................................................... 62
de estado de Dranchuk-Pulvis-Robinson ...................................................... 63
Tabla 3.8- Resultados del factor de compresibilidad Z utilizando el método
gráfico de Standing y Katz ............................................................................ 64
gráfico de Eykman ........................................................................................ 66
Pitzer............................................................................................................. 67
Tabla 3.11- Análisis descriptivo de los valores de Z calculados ................... 68
Tabla 3.12- Resultados del factor de compresibilidad Z para cada método
definido ......................................................................................................... 69
Tabla 3.13- Rangos asignados a cada método y resultados de Rj ............... 69
xiii
NÚCLEO DE MONAGAS
RESUMEN
AUTORES ASESOR ACADÉMICO
Yérika Glecmar Carrera Marcano Msc. Tomás Marín

Javier Eduardo Castillo La Rosa
Félix David Roque Pérez
Mayo, 2018
El factor de compresibilidad Z, es uno de los parámetros que, con mayor

precisión diferencia el comportamiento de los fluidos en estado líquido del
estado gaseoso. Además de esto, define el comportamiento de los gases a
determinadas condiciones de presión y temperatura. Lo que hace que sea de
vital importancia, para todos los diseños e instalaciones que trabajan con
fluidos compresibles. El objetivo principal de esta monografía de
investigación fue realizar un análisis comparativo entre nueve métodos
seleccionados para el cálculo del factor de compresibilidad Z del gas natural,
mediante la aplicación de la prueba estadística de Friedman, con la finalidad
de establecer cuál de ellos es el más apropiado para un resultado óptimo y
eficaz tomando en cuenta cuatro muestras de gas obtenidas de tesis de
investigación con sus respectivas condiciones de operación. El análisis
estadístico indicó que de acuerdo a los resultados alcanzados, el método
más apropiado para el cálculo del factor de compresibilidad Z, fue el de
Soave-Redlich-Kwong y el menos apropiado el de Peng-Robinson, tomando
en cuenta las cromatografías mencionadas anteriormente.
Palabras claves: Factor de compresibilidad, ecuaciones de estado, prueba

de Friedman, análisis estadístico.
xiv
INTRODUCCIÓN
El gas natural es una mezcla homogénea de hidrocarburos parafínicos

de la forma CnH2n+2 y cantidades menores de gases inorgánicos como el N 2,
CO2, H2S y otros, de composición variable y acotada dentro de rangos
específicos (Irragorry, 2004). De acuerdo a este autor: Un gas ideal es aquel
que tiene sus partículas en movimiento libre y no tienen atracción entre ellas.
Este obedece exactamente a la ecuación de estado PV=n.R.T.
Ahora bien el factor de compresibilidad Z, es un factor de corrección,

que se introduce en la ecuación de estado de gas ideal descrita
anteriormente, para modelar el comportamiento de los gases reales, los
cuales se pueden comportar como gases ideales para condiciones de baja
presión y alta temperatura, tomando como referencia los valores del punto
crítico. A nivel bibliográfico podemos encontrar numerosos estudios que se
han hecho sobre el factor de compresibilidad, algunos de ellos son: Standing-
Katz, Pitzer, Van der Waals, Peng-Robinson, Redlich-Kwong, entre otros. La
referida literatura presenta también cálculos utilizando los métodos anteriores
y demuestra que existe diferencia entre uno y otro.
La presente monografía de investigación tiene por objeto realizar un

análisis comparativo de las metodologías existentes para la determinación
del factor de compresibilidad del gas (Z) aplicado a muestras de gas natural.
Para ello se desarrollará todo lo concerniente al fundamento teórico del tema
en estudio, seguidamente se hará la revisión de los modelos existentes más
utilizados para la determinación del factor de compresibilidad del gas y por
último se determinará el factor de compresibilidad a muestras de gas
aplicando los diferentes métodos seleccionados y se analizaran los
resultados estadísticamente.
1
CAPÍTULO I
EL PROBLEMA Y SUS GENERALIDADES
1.1. PLANTEAMIENTO DEL PROBLEMA
El comportamiento de presión, volumen, temperatura y moles para un

gas natural puede ser fácilmente representado por las leyes de los gases
ideales. No obstante, esto es posible si hablamos de presiones bajas y
temperaturas altas. A medida que la presión se va incrementando los
parámetros que se conjugan en la ecuación de los gases ideales (PV=n.R.T)
comienza a presentar cierta desviación. Esta desviación puede presentarse
por encima o por debajo al valor que arroja las dimensiones del gas en
cuestión. (Romero, 2003)
El factor de desviación o factor de compresibilidad (Z) de un gas natural
se define como la relación que existe entre el volumen real de un gas medido
a condiciones de presión y temperatura y el volumen ideal medido a las
mismas condiciones (Gonzalez, 2006). Debido a la necesidad de determinar
este valor, diversos autores han presentado sus modelos específicos. El libro
titulado “Características y comportamiento de los hidrocarburos” (Pérez y
Martínez 1994) presentan los siguientes: Standing–Katz, Standing–Katz
modificado, Pitzer, Papay, Van Der Waals, Redlich–Kwong, Soave–Redlich–
Kwong, Peng–Robinson, Starling, Sarem, Berthelot, Wohl, Esferas Duras,
Hall Y Yarborough, Beattie-Bridgeman, Benedict–Webb–Rubin, Refracción
Molecular de Eykman, Dranchuk–Purvis–Robinson, Dranchuk–Abou-Kassen,
Gray–Sims, Carlile–Gillett, Brill. Haciendo una comparación de los diferentes
métodos, se pueden visualizar diferencias considerables entre cada uno de
ellos.
2
3
Los gases naturales de Venezuela presentan características típicas que

deben ser considerados en el desarrollo de estas correlaciones. Es por ello
que la presente monografía plantea un análisis comparativo de los diferentes
métodos utilizados para la determinación del factor de compresibilidad del
gas a varias muestras, y de este modo analizar de forma estadística las
posibles diferencias entre los resultados.
1.2. OBJETIVOS DE LA INVESTIGACIÓN
1.2.1 Objetivo general
Analizar comparativamente las diferentes metodologías para la

determinación del factor de compresibilidad (Z) del gas aplicado a muestras
de gas natural.
1.2.2 Objetivos específicos
 Describir los métodos más utilizados para la determinación del factor
de compresibilidad Z.
 Determinar el factor de compresibilidad Z a cuatro (4) muestras de gas
natural, utilizando nueve (9) métodos representativos para su cálculo.
 Comparar a través de un análisis estadístico las diferencias existentes
entre los resultados obtenidos.

4
1.3. JUSTIFICACIÓN DE LA INVESTIGACIÓN
El factor de compresibilidad Z, es uno de los parámetros que, con

mayor precisión diferencia el comportamiento de los fluidos en estado líquido
del estado gaseoso. Además de esto, define el comportamiento de los gases
a determinadas condiciones de presión y temperatura. Lo que hace que sea
de vital importancia, para todos los diseños e instalaciones que trabajan con
fluidos compresibles.
Ahora bien, cualquier estudio o investigación que se desarrolle para

producir una herramienta que facilite o simplifique los cálculos de ingeniería,
es por sí misma justificada, sobre todo si en su ejecución, los recursos
invertidos serán la intelectualidad, la práctica y la ampliación de
conocimientos y técnicas que al adquirirlas producirían un significativo
crecimiento de las personas que lo realizan, en este caso nosotros como
futuros ingenieros.
Es por ello que esta investigación busca entre un universo de

metodologías existentes para el cálculo del factor de compresibilidad Z,
analizar un número finito de dichos métodos y que esto permita establecer en
base a criterios estadísticos, cuál de ellos es el más adecuado para su
cálculo, dependiendo claro está, de las condiciones o limitantes que
establezca cada método.
CAPÍTULO II
MARCO TEÓRICO
2.1. ANTECEDENTES DE LA INVESTIGACIÓN
Barrios, R. y Martínez, K. (2011) Universidad de Oriente. Realizaron el

trabajo de grado “Validación del factor de compresibilidad Z, a través del
gasoducto ACOGAS de la planta de inyección de gas El Furrial
mediante ecuaciones de estado”. Para ello utilizaron la prueba estadística
de Friedman. Según la metodología y los resultados obtenidos se pudo
concluir que el método más preciso para el cálculo de Z, tomando en cuenta
las condiciones de operación del gasoducto ACOGAS es la ecuación de
estado de Van der Waals, seguidamente por la ecuación de Soave-Redlich-
Kwong.
El aporte de este trabajo a nuestra investigación se basó principalmente
en las cromatografías obtenidas a través de 2 reportes de inyección diario
del cromatografo ACOGAS, que sirvieron como muestras de gas, para la
aplicación de los métodos desarrollados.
Rodríguez, R. (2014) Universidad de Carabobo. Realizó el trabajo de

ascenso titulado: “Estimación de propiedades termodinámicas de
mezclas de gases en condiciones reales utilizando una herramienta
computacional”. Desarrolló una herramienta computacional para determinar
las propiedades termodinámicas de gases y sus mezclas. El algoritmo de
cálculo fue codificado en Visual Basic 6.0. Según la metodología y los
resultados obtenidos se pudo concluir que el programa desarrollado permite
analizar de manera eficiente la relación funcional que hay entre las diversas
propiedades de los gases y sus mezclas. Al comparar con resultados
5
6
experimentales, permite discernir sobre la precisión y la aplicabilidad entre

los diversos modelos propuestos.
Su aporte a nuestra investigación radica en los cálculos que realizaron

utilizando los modelos de Van der Waals, Redlich-Kwong, Soave-Redlich-
Kwong y Peng-Robinson, lo cual sirvió de guía para desarrollar las
ecuaciones de estado propuestas.
2.2. BASES TEÓRICAS
2.2.1 Estado gaseoso
El estudio de los gases ha proporcionado claves para comprender tanto

fenómenos físicos como químicos. Históricamente los gases ocupan un lugar
especial, debido a que su estudio llevó a desarrollar gran cantidad de
conocimientos, de ahí la importancia de las características que distinguen al
estado gaseoso. La cantidad de un gas se determina midiendo sus
propiedades como lo son el volumen, la presión y su temperatura (Triana,
2012), por lo que es importante hablar de estas variables:
Volumen
Una de las propiedades características de los gases es su carencia de

forma y su expansión ilimitada debido a su estructura interna, es decir, que
todos los gases se expanden indefinidamente hasta llenar el espacio dentro
del cual están contenidos; por eso una determinada muestra de gas no tiene
ni forma, ni volumen definido. El estado gaseoso también se caracteriza por
su alto grado de compresibilidad. Para producir una disminución mínima en
el volumen de un líquido o un sólido se requiere aplicar una presión enorme,
7
mientras que un gas fácilmente puede ser reducido a una pequeña fracción
de su volumen original (Cotton, 1986).
Presión
La razón por la cual los gases ejercen presión, es que estos en una
especie de bombardeo continuo (Movimiento Browniano), golpean las
paredes del recipiente que los contienen. Clásicamente, se define la presión,
como la fuerza (F) por unidad de área (A), por lo que la presión de un gas, es
la fuerza que éste ejerce sobre el recipiente que lo contiene, dividida por el
área de superficie del recipiente:
𝐹𝑢𝑒𝑟𝑧𝑎
𝑃𝑟𝑒𝑠𝑖ó𝑛 = (2.1)
Á𝑟𝑒𝑎
Temperatura
La temperatura es una propiedad de un sistema que representa la

existencia de una condición del equilibrio térmico que es independiente de la
composición y del tamaño del sistema. Esta propiedad en un cuerpo es la
misma que determina el flujo de calor hacia otro cuerpo o de otros cuerpos
hacia él (Petrucci, 1988).
2.2.2 Gas ideal
Se ha dado el nombre de gas ideal o gas perfecto a un gas teórico

cuyas moléculas no ocupan espacio ni ejercen atracción alguna entre ellas,
por lo que en él se cumplen todas las leyes que aplican a los gases. La
principal ecuación que define el estado gaseoso en condiciones ideales es
8
PV=n.R.T (Pino, 2008). Adicionalmente, los choques entre las moléculas son
completamente elásticos. A bajas presiones, muchos gases se comportan
como un gas ideal. Cuando la presión del gas se incrementa se observan
desviaciones de su comportamiento respecto a un comportamiento ideal.
Ley de Boyle: efecto de la presión sobre el volumen
Este fenómeno fue estudiado independientemente por Robert Boyle

(1627-1691) y por Edmund Mariotte (1620-1684), quienes comprobaron
experimentalmente que el volumen de una determinada cantidad de gas
ideal, cuando la temperatura se mantiene constante, es inversamente
proporcional a la presión que se ejerce sobre el gas. Matemáticamente, esto
se puede representar como:
(2.2)
P=Vα 1
Boyle llevó a cabo experimentos con el aire para encontrar la relación

(2.2)
que existe entre el volumen y la presión. Tomó un tubo en forma de “U” y
vertió mercurio hasta que los niveles en ambas ramas fuesen iguales. De
esta forma, se tiene el volumen del aire V a la presión atmosférica P, ya que
las ramas del mercurio tienen el mismo nivel. Enseguida, vertió más mercurio
hasta lograr un desnivel. Midiendo la columna de aire encerrado, obtuvo el
nuevo volumen V1. Después de varias lecturas con este dispositivo Boyle
encontró que al incrementar la presión sobre un gas confinado, éste reducía
su volumen en la misma proporción. Con base en estas observaciones Boyle
estableció que: “El volumen de un gas seco varía inversamente con la
presión ejercida sobre él, si la temperatura permanece constante”. La
siguiente figura muestra la relación entre volumen y presión estudiada por
Boyle.
9
Escrita en forma de ecuación, resulta:
𝑃1 𝑉
= 𝑉2 𝑃1 𝑉1 = 𝑃2 𝑉2 (2.3)
𝑃2 1
PV = constante
Ley de Charles: efecto de la temperatura sobre el volumen
El físico francés Jacques Charles (1763-1823) a principios del siglo XIX

utilizó el mismo diseño empleado un siglo antes por Boyle, pero ahora
variando la temperatura y manteniendo constante la presión. La siguiente
figura muestra la relación entre volumen y temperatura analizada por Charles
(Triana, 2012).
Charles estableció que, si la presión ejercida sobre una cantidad

determinada de gas se mantiene constante, para cualquier cambio de estado
que sufra el gas: el volumen del mismo varía directamente proporcional a su
temperatura absoluta, escrita en forma de ecuación es:
10
𝑉1 𝑇1 𝑇1 𝑇2 𝑇
= = = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 (2.4)
𝑉2 𝑇2 𝑉1 𝑉2 𝑉
Ley de Gay Lussac: relación entre presión y temperatura
La relación existente entre los cambios de temperatura y la presión a

volumen constante se le atribuyen a Joseph Louis Gay-Lussac y a Guillaume
Amontons (1663-1705) (Triana, 2012). Sus observaciones se resumen en la
ley de Gay-Lussac, la cual establece que: si el volumen de una cantidad
determinada de gas se mantiene constante, para cualquier cambio de estado
que sufra el gas, la presión absoluta a la cual es sometida el gas varía
directamente proporcional a su temperatura absoluta, escrita en forma de
ecuación es:
𝑃1 𝑇1 𝑃1 𝑃2 𝑃 (2.5)
= = = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝑃2 𝑇2 𝑇1 𝑇2 𝑇
Ley combinada
Existen tres variables que describen una determinada cantidad de gas:

P, V y T; es por ello que se estudia la dependencia entre estas propiedades,
cuando todas las demás se conservan constantes. La ley de Boyle y la ley de
Charles han sido verificadas y confirmadas mediante cuidadosas
realizaciones experimentales. En ciertas condiciones de temperatura y/o
presión, las propiedades de la mayor parte de los gases reales se desvían
por completo de las de un gas ideal (Sabino, 2010). Las dos primeras leyes
pueden utilizarse, para averiguar el nuevo volumen que adquiere un gas
cuando se modifica su temperatura y su presión, pero no cuando ambas
variables lo hacen de manera simultánea. Debido a que normalmente los
11
gases están sometidos a cambios simultáneos de presión, volumen y

temperatura, las leyes de Boyle y Charles se combinan para cuantificar esos
efectos como se muestra a continuación: Si se supone que una cantidad
cualquiera de gas está a una temperatura inicial (T 1), presión inicial (P1),
volumen inicial (V1) y queremos saber que volumen final (V2) ocupará cuando
se lleva a una presión final (P2) y temperatura final (T2), es decir:
V1, P1, T1 Vx, P2, T1 V2, P2, T2 (2.6)
En el primer paso de V1 a Vx, la temperatura se mantiene constante.

Según Boyle, en el segundo paso, de Vx a V2, la presión se mantiene
constante. Luego según Charles y como los dos volúmenes V x son iguales,
las ecuaciones también lo son. Ordenando la expresión:
P2 V2 T 1 = P1 V1 T 1 (2.7)
Ley de Avogadro
Amadeo Avogadro propuso una ley en el siglo XIX, la cual establece

que bajo las mismas condiciones de presión y temperatura, volúmenes
iguales de cualquier gas ideal contienen el mismo número de moléculas. Se
ha demostrado que hay 733x1026 moléculas en cada lb-mol de cualquier gas
ideal. De esta ley, se puede concluir que el peso de un volumen de gas es
una función de los pesos de las moléculas y que existe un volumen en el cual
el gas pesa en valor numérico igual a su peso molecular.
12
Ley de los gases ideales
Esta Ley relaciona la temperatura, presión y volumen para gases

perfectos. Es una combinación de las leyes de Boyle, Charles (conocida
también como la Ley de Gay-Lussac) y de Avogadro (Pérez y Martínez,
1994). En forma de ecuación, se expresa como:
PV = n.R.T (2.8)
Donde:
P: presión absoluta del gas
R: Constante del gas o constante universal de los gases en las unidades
apropiadas
T: temperatura absoluta del gas
V: volumen total ocupado por el gas
n: número de moles del gas
El valor de la constante R depende de las unidades empleadas de

presión, volumen y temperatura. La tabla 2-1 incluye diferentes valores de R
(GPSA, 1987).
Tabla 2.1–Valores de la constante general de los gases para diferentes

unidades de presión (P), volumen (V) y temperatura (T)
P V T n R
Atm. Litros ºK (273,16 + ºC) g/M 0,082057
Atm. pie3 ºR (460 + ºF) lb/M 0,73024
lpca pulg3 ºR (460 + ºF) lb/M 18.510,0
lpca pie3 ºR (460 + ºF) lb/M 10,732
Fuente: Pérez y Martínez, 1994
13
2.2.3 Gas real
Son aquellos que se considera están formados por partículas con

volumen entre las cuales existen fuerzas de atracción y repulsión. En
realidad, en la naturaleza, sólo existen gases reales. La denominación de
gas ideal o real dependerá de las condiciones en que sea estudiado el gas.
Sin embargo, a medida que aumenta la presión y/o la temperatura hay que
incluir ciertas modificaciones en la ecuación de los gases, por tanto, es de
gran importancia tener una noción de las desviaciones que sufren éstos bajo
determinadas condiciones de temperatura, presión y volumen (Pino, 2008).
2.2.4 Gas natural
El gas natural es una mezcla homogénea de hidrocarburos gaseosos

que se encuentra frecuentemente en yacimientos fósiles, solo o
acompañando al petróleo o a los depósitos de carbón; formado en
proporciones variables de hidrocarburos parafínicos, denominados alcanos,
éstos debido al tipo de hibridación (sp3), indica que debe existir un enlace
simple entre los carbonos, por lo que reciben el nombre de hidrocarburos
saturados, para diferenciarlos de los alquenos (sp2) y alquinos (sp) que se
denominan insaturados. El componente principal del gas natural es el
metano (CH4), cuyo contenido varía generalmente entre 60 y 90% molar.
Posee cantidad de gases inorgánicos, como el nitrógeno (N2), el dióxido de
carbono (CO2), sulfuro de hidrógeno (H2S), monóxido de carbono (CO),
oxígeno (O2), vapor de agua (H2O), entre otros; todos estos componentes
son considerados impurezas del gas natural (Martínez y Ramos, 2010).
14
2.2.5 Clasificación del gas natural según su origen
Gas asociado: es un gas que se encuentra en contacto y/o disuelto en el

petróleo crudo del yacimiento. Puede ser clasificado como gas de casquete
(libre) o gas en solución (disuelto). El gas natural asociado debe pasar por
distintos procesos de separación (extracción por absorción, fraccionamiento,
entre otros) para obtener las diferentes fases (líquidas, sólidas y gaseosas).
Más del 90% de las reservas de gas natural del país es de gas asociado.
Gas no asociado o libre: este gas se encuentra en yacimientos que no

contienen hidrocarburos líquidos. El gas natural no asociado está unido con
agua en yacimientos de gas seco, permaneciendo en fase gaseosa a
condiciones de yacimiento y superficie, este puede ser tratado y distribuido
sin necesidad del proceso de separación.
2.2.6 Clasificación del gas natural según su composición
Para determinar la composición del gas natural se utilizan

procedimientos de análisis, para establecer los porcentajes de los distintos
hidrocarburos y de los demás componentes en el gas natural y productos
afines. Entre los métodos o aparatos comúnmente usados para la
determinación de los hidrocarburos se encuentran el fraccionamiento a baja
temperatura, espectrómetro de masa, espectrómetro de infracción infrarroja y
el análisis cromatográfico de gases. (Martínez y Ramos, 2010). A través de
estos procedimientos, se puede establecer la siguiente clasificación según la
composición de los gases:
Gas pobre o seco: no contiene hidrocarburos que se licuan a temperatura y

presión ambiente. Tiene un contenido bajo de propano hasta heptano. En el
15
gas seco, el componente principal es el metano, cuya composición

porcentual alcanza valores mayores al noventa por ciento en la relación
molar (C1>90%), el gas natural por lo general tiene pequeñas cantidades de
pentano y componentes más pesados (C5<1%).
Gas rico (húmedo): gas predominantemente con metano y otros

hidrocarburos (propano y compuestos más pesados). Una corriente de gas
con unos 2-3 galones de hidrocarburos líquidos por mil pies cúbicos
normales de gas (GPM), se considera como riqueza líquida, por lo que es
muy rentable a los propósitos de la inversión, dependiendo del caudal
disponible para este propósito.
Gas condensado: se define como un gas con líquido disuelto. El contenido

de metano es mayor al 60% molar, mientras que el contenido de heptano y
fracciones más pesadas (C7+) deben ser mayores al 12,5% molar. Aquí la
mezcla de hidrocarburos a las condiciones iniciales de presión y temperatura
se encuentra en fase gaseosa o en el punto de rocío.
Gas ácido: contiene cantidades de sulfuro de hidrógeno (H2S) mayor a

cuatro (4) partes por millón en base a volumen (4ppm, V), más de 3% de
dióxido de carbono (CO2) y otros componente ácidos, razón por la cual se
vuelve corrosivo en presencia de agua libre. El gas natural ácido es
independiente si el gas es seco o pobre, rico o húmedo y condensado.
Gas dulce: éste es un gas que contiene cantidades de sulfuro de hidrógeno

(H2S) menores a cuatro (4) partes por millón en base a volumen (4 ppm, V) y
menos de 3% en base molar de dióxido de carbono (CO2).
16
Gas hidratado: este es un gas que tiene más de siete (7) libras de agua por
cada millón de pies cúbicos normales de gas (7 lb de Agua/MMPCN), lo que
indica que el gas deberá ser sometido al proceso de deshidratación, para
poder comercializarlo.
Gas anhidro: es un gas que tiene menos de siete (7) libras de agua para
cada millón de pies cúbicos normales de gas (7lb de agua / MMPCN).
2.2.7 Mezcla de Gases
La composición de una mezcla de gases, se expresa generalmente en

porcentajes por mol, porcentaje (o fracción) por volumen o porcentaje por
peso. Si una mezcla de gases se compone de 1, 2, 3,..i,..n componentes, el
porcentaje por mol, por volumen o por peso es la razón del número de moles,
ni, del volumen Vi, o del peso W i de un componente i en la mezcla, al número
total de moles, Σni, volumen total, ΣVi, o peso total, ΣW i, del gas
respectivamente, multiplicado por 100. Analíticamente:
𝑊𝑖
% 𝑊𝑖 = ∑𝑁
𝑥 100 (2.9)
𝑗=1 𝑊𝑗
𝑉𝑖
% 𝑉𝑖 = ∑𝑁
𝑥 100 (2.10)
𝑗=1 𝑉𝑗
𝑛𝑖 (2.11)
% 𝑛𝑖 = 𝑁 𝑥 100
∑𝑗=1 𝑛𝑗
17
Peso molecular aparente
Martínez, (1994) plantea que, cuando se trata de mezcla de gases,

generalmente no se habla de peso molecular en el mismo sentido de
componentes puros, ya que una mezcla no se puede representar como una
formula sencilla. Sin embargo, su comportamiento puede estudiarse
considerando el denominado “peso molecular aparente” o “peso molecular
promedio molar”. En forma similar a componentes puros, se define como el
peso de 379,6 pie3 de la mezcla gaseosa a 60º F y 14,7 lpca. Usando este
valor, su comportamiento es igual que para gases ideales. El peso molecular
aparente puede calcularse a partir de la composición molecular (por mol) de
la mezcla y de los pesos moleculares individuales de los componentes. Si Yi
es la fracción molecular de un componente i, de peso molecular Mi en una
mezcla de n componentes, el peso molecular aparente será:
𝑀𝑎 = ∑𝑁
𝑖=1 𝑦𝑖 𝑀𝑖 (2.12)
Propiedades pseudocríticas
Se conoce como temperatura crítica (Tc) aquella temperatura a partir de

la cual la sustancia sólo puede existir en estado gaseoso y la presión crítica
(Pc) aquella necesaria para provocar la licuefacción a la temperatura crítica.
Cada mezcla de gases tendrá su propia temperatura y presión crítica
verdaderas y para conocerlas habría que determinarlas experimentalmente.
Para el estudio del comportamiento de la mezcla, Kay incluyo la definición de
la temperatura pseudocrítica, lo que indica que una mezcla gaseosa de n
componentes, Yi es la fracción molar del componente i de presión crítica Pci
y temperatura crítica Tci, se tiene:
18
𝑃𝑠𝑐 = ∑ 𝑌𝑖 ∗ 𝑃𝑐𝑖 (2.13)
𝑇𝑠𝑐 = ∑ 𝑌𝑖 ∗ 𝑇𝑐𝑖 (2.14)
Cuando se conoce la composición del gas, fácilmente se calculan los

valores de Psc y Tsc, ya que los valores de Pc y Tc para los componentes
puros son datos tabulados (Barrios y Martínez, 2011).
2.2.8 Ecuaciones de estado
Las ecuaciones de estado se usan para correlacionar datos PVT y

calcular diferentes propiedades físicas y termodinámicas de sistemas de
hidrocarburos en un amplio intervalo de valores de presión y temperatura.
Con el auxilio de los computadores, tales ecuaciones tienen un uso y una
aplicación cada vez mayores. Son ecuaciones más o menos complicadas,
pero, en general, bastante exactas. Las principales ecuaciones de estado
son las siguientes:
Ecuación de Van Der Waals
Puede escribirse en la siguiente forma:
𝑎 (2.15)
(𝑃 + 𝑉 2 ) (𝑉 − 𝑏) = 𝑅𝑇
𝑅𝑇 𝑎 (2.16)
𝑃= − 2
𝑉−𝑃 𝑉
Ya que por definición, Z = PV/RT. Además, puede formularse:

19
𝑃𝑉 𝑏 𝑎
=𝑍 =1+ − (2.17)
𝑅𝑇 𝑉−𝑏 𝑅𝑇𝑉
Donde:
27𝑅2 𝑥𝑇𝑐 2 𝑅𝑥𝑇𝑐 8𝑃𝑐𝑥𝑉𝑐

𝑎= 𝑏= 𝑅=
64𝑃𝑐 8𝑃𝑐 3𝑇𝑐
Introduciendo estos valores en la ecuación anterior y expresándola en

términos de cantidades reducidas, resulta:
(2.18)
3
( 𝑃𝑟 + 2 ) (3𝑉𝑟 − 1) = 8𝑇𝑟
𝑉𝑟
La ecuación de Van der Waals posee dos constantes, y expresada en

función de Z tenemos (GPSA, 1987)
𝑍 3 − (1 + 𝐵)𝑍 2 + 𝐴𝑍 − 𝐴𝐵 = 0 (2.19)
Donde:
𝑎𝑃 𝑏𝑃
𝐴= (2.20) 𝐵= (2.21)
𝑅2 𝑇 2 𝑅𝑇
Ecuación de Berthelot
Es muy similar a la ecuación de Van der Waals. Sólo varía en que usa
a/TV2 en lugar de a/V2. Los valores de a, b y R son los mismos, sin embargo,
sugiere una mejora y usa el valor de b = 9RTc/128 Pc.
20
Ecuación de Redlich-Kwong, (R-K)
Esta ecuación también contiene dos constantes y pueden determinarse

a partir de la presión y temperatura crítica. No obstante, la ecuación R-K es
más exacta y tiene una mayor aplicación. Se escribe de la siguiente manera:
𝑎
(𝑃 + ) (𝑉 − 𝑏) = 𝑅𝑇 (2.22)
√𝑇 𝑉 (𝑣+𝑏)
𝑅𝑇 𝑎
𝑃= − (2.23)
𝑉−𝑏 √𝑇 𝑉 (𝑉+𝑏)
Donde:
0,4278𝑅2 𝑥𝑇𝑐 2,5 0,0867𝑅𝑥𝑇𝑐
𝑎= 𝑏=
𝑃𝑐 𝑃𝑐
Para mezclas de gases, pueden usarse varias reglas de combinación,

como se presentan a continuación:
Regla A:
𝑛 2 𝑛
𝑎𝑚 = (∑ 𝑦𝑖 𝑎𝑖 ) 𝑏𝑚 = ∑ 𝑦ᵢ 𝑏ᵢ
𝑖=1 𝑖=1
Regla B:
𝑛 𝑛 𝑛
𝑎𝑚 = ∑ ∑ 𝑦ᵢ 𝑦𝑗 𝑎𝑖𝑗 𝑏𝑚 = ∑ 𝑦ᵢ 𝑏ᵢ
𝑖=1 𝑗=1 𝑖=1
21
El valor de aij está dado por:

𝑎𝑖𝑗 = √𝑎𝑖 𝑎𝑗
La ecuación R-K es comúnmente escrita en la siguiente forma:
𝑃𝑉 1 𝐷2 ℎ
=𝑍= − 𝑥 (2.24)
𝑅𝑇 1−ℎ 𝐸 1+ℎ9
Donde:
𝑎 𝑏 0,0867 𝐸𝑃 𝑏𝑃 𝑏
𝐷2 = 𝐸= = ℎ= = =
𝑅2 𝑥𝑇 2,5 𝑅𝑇 𝑃𝐶 𝑇𝑟 𝑍 𝑍𝑅𝑇 𝑉
Ya que Z aparece en ambos lados de la ecuación R-K (h contiene el

término Z), ésta puede expresarse en función de Z, lo que la hace muy
apropiada para programas en computadores.
𝑓(𝑍) = 𝑍 3 − 𝑍 2 + (𝐷2 𝑃 − 𝐸𝑃 − 𝐸 2 𝑃2 )𝑍 − (𝐷2 𝐸𝑃2 ) = 0 (2.25)
La ecuación original de Redlich-Kwong fue modificada en el sentido de

que los parámetros "a" y "b" se determinan incluyendo una serie de datos
experimentales. En esta forma, se habla de la ecuación de Redlich-Kwong
modificada.
Ecuación de Soave Redlich Kwong (SRK)
𝑍 3 − 𝑍 2 + (𝐴 − 𝐵 − 𝐵2 )𝑍 − 𝐴𝐵 = 𝑜 (2.26)
22
Donde:
𝑎𝑃 𝑏𝑃
𝐴= 𝐵=
𝑅2 𝑇𝐶2 𝑅𝑇𝐶
𝑎 = 0,42747 ( )𝛼 𝑏 = 0,08664 ( )
𝑃𝐶 𝑃𝐶
𝑇
√𝛼 = 1 + 𝑚 [1 − √ ]
𝑇𝐶
𝑚 = 0,48 + 1,574𝑤 − 0,176𝑤 2
Ecuación de Peng Robinson (PR)
𝑍 3 − (1 − 𝐵)𝑍 2 + (𝐴 − 3𝐵2 − 2𝐵)𝑍 − (𝐴𝐵 − 𝐵2 − 𝐵3 ) = 0 (2.27)
Donde:
𝑎𝑃 𝑏𝑃
𝐴= 𝐵=
𝑅2 𝑥𝑇𝐶2 𝑅𝑇𝑐
𝑎 = 0,45724 ( )𝛼 𝑏 = 0,0778 ( )
𝑃𝐶 𝑃𝑐
𝑇
√𝛼 = 1 + 𝑚 [1 − √ ]
𝑇𝑐
𝑚 = 0,37464 + 1,54226𝑤 − 0,26992𝑤 2

23
Ecuación de Wohl
Esta ecuación contiene tres constantes, también obtenidas a partir de

las propiedades críticas (presión, temperatura, volumen). La ecuación es:
𝑎 𝐶 (2.28)
(𝑃 + − 2 3 ) (𝑉 − 𝑏) = 𝑅𝑇
𝑇𝑉 (𝑉 − 𝑏) 𝑇 𝑉
𝑅𝑇 𝑎 𝐶
𝑃= − + 2 3 (2.29)
𝑉 − 𝑏 𝑇𝑉(𝑉 − 𝑏) 𝑇 𝑉
𝑃𝑉 𝑉 𝑎 𝐶
=𝑍= − + (2.30)
𝑅𝑇 𝑉−𝑏 𝑅𝑇 2 (𝑉−𝑏) 𝑅𝑉 2 𝑇3
Donde:
𝑉𝐶
𝑎 = 6𝑃𝑐 𝑇𝑐 𝑉𝑐2 𝑏= 𝑐 = 4𝑃𝑐 𝑇𝑐2 𝑉𝑐3
4
𝑅𝑇𝐶 15
=
𝑃𝐶 𝑉𝐶 4
Ecuación de esferas duras
Las ecuaciones de estado, denominadas de esferas duras,

consideran que un fluido está compuesto de esferas inelásticas rígidas, en
lugar de moléculas reales. Aunque a primera vista este concepto no tiene
mucho sentido, presenta una serie de ventajas, ya que a altas presiones y
temperaturas, los fluidos reales se comportan muy similares a fluidos
compuestos de esferas duras. Además, este sistema lleva a un desarrollo
teórico simple, característica de una ecuación de estado práctica.
24
Hall y Yarborough
Presentan, en forma muy práctica, un tipo de ecuación de estado de

esferas duras, para su uso en la industria del gas natural.
𝑃 1+𝑦+𝑦2 −𝑦3
Z= = − (14,76 − 9,76𝑡 + 4,58𝑡 2 ) 𝑡𝑦 +
𝜌𝑅𝑇 (1−𝑦)3
(90,7 − 242,2𝑡 + 42,4𝑡 2 )𝑡𝑦 (1,18+2,82𝑡) (2.31)
Además: (
𝑏𝑃𝐶 2] 𝑏𝜌
= 2,45. 𝑒 [−1,2.(1−𝑡) 𝑦= 2
𝑅𝑇𝐶 4
.
2
Donde: 9
3
b = constante de Van der Waals, pie /lbmol, y t = 1/Tr = Tc/T = Tsc/T, inverso )
de la temperatura reducida. Esta ecuación es muy práctica para programas
de computadoras.
Ecuación de Beattie-Bridgeman
𝑎
𝑅𝑇 𝐶 𝑏 𝐴𝑂 (1−𝑉)
𝑃= (1 − ) [𝑉 + 𝐵𝑂 (1 − )] − (2.32)
𝑉2 𝑉𝑇 3 𝑉 𝑉2
En ella Z = PV/RT, de donde:
𝑎
𝑃𝑉 1 𝐶 𝑏 𝐴𝑂 (1− )
𝑉
= 𝑍 = (1 − 3
) [𝑉 + 𝐵𝑂 (1 − )] − (2.33)
𝑅𝑇 𝑉 𝑉𝑇 𝑉 𝑇𝑅𝑉
25
Ecuación de de Starling
Es una ecuación similar a la de BWR. Igualmente, la presión (P) está

expresada en función de la densidad molar (ρ) la temperatura absoluta y una
serie de constantes. En este caso, son once constantes, en lugar de ocho,
las que contiene la ecuación de BWR. Ellas son: Ao, Bo, Co, Do, Eo, a, b, c, d,
α y Ƴ. La ecuación es la siguiente:
𝐶𝑂 𝐷𝑂 𝐸𝑂 2 𝑑
𝑃 = 𝜌𝑅𝑇 + (𝐵𝑂 𝑅𝑇 − 2
− 3 − 4 ) 𝜌 + (𝑏𝑅𝑇 − 𝑎 − ) 𝜌3 +
𝑇 𝑇 𝑇 𝑇
𝑑 1+Ƴ𝜌2
𝛼 (𝑎 + 3
) 𝜌6 + 𝑐𝜌3 2 (2.34)
𝑇 𝑇 2 𝑒 (Ƴ 𝜌 )
𝑃𝑉 𝐴𝑂 𝐶𝑂 𝐷𝑂 𝐸𝑂
= 𝑍 = 1 + (𝐵𝑂 − − 3
− 4
− )𝜌 +
𝑅𝑇 𝑅𝑇 𝑅𝑇 𝑅𝑇 𝑅𝑇 5
𝑎 𝑑 𝑎 𝑑 1+Ƴ𝜌2
(𝑏 − − 2
) 𝜌2 + 𝛼 ( + 2
) 𝜌5 + 𝑐𝜌2 2 (2.35)
𝑅𝑇 𝑅𝑇 𝑅𝑇 𝑅𝑇 𝑅𝑇 3 𝑒 (Ƴ 𝜌 )
La anterior ecuación ha sido muy usada para el cálculo de propiedades

termodinámicas de hidrocarburos (entalpía, entropía, fugacidad, densidad,
etc.) y constantes de equilibrio. Se han calculado los diagramas de presión-
entalpía para el metano, etano, propano, n-butano, isobutano, n-pentano,
isopentano, hexano, heptano, octano, etileno, propileno, sulfuro de
hidrógeno, dióxido de carbono y nitrógeno. Estos diagramas cubren la fase
gaseosa, la fase liquida y la región de dos fases. La GPSA usa las
propiedades termodinámicas calculadas con la ecuación de Starling y
enumera los principales modelos matemáticos utilizados actualmente
(Redlich-Kwong, Soave Redlich Kwong (SRK), Peng Robinson y Benedict-
Webb-Rubin-Starling o BWRS). (GPSA, 1987).
26
Ecuación de de Martin-Hou
Se han desarrollado otras ecuaciones de estado, como la de Martin-

Hou y las denominadas ecuaciones de estado viriales que se expresan así:
𝑃𝑉 𝐴2 𝐴3 𝐴3
=𝑍 =1+ + + + ⋯, (2.36)
𝑅𝑇 𝑉 𝑉2 𝑉3
Donde A2, es el segundo coeficiente virial, A3, es el tercer coeficiente

virial, etc. Estos coeficientes son funciones sólo de temperatura. Tanto la
ecuación de Martin-Hou como las viriales han sido poco empleadas en
cálculos de sistemas de hidrocarburos.
2.2.9 Factor de compresibilidad (Z)
El factor de compresibilidad Z, es un factor de corrección, que se

introduce en la ecuación de estado de gas ideal para modelar el
comportamiento de los gases reales, los cuales se pueden comportar como
gases ideales para condiciones de baja presión y alta temperatura, tomando
como referencia los valores del punto crítico, es decir, si la temperatura es
mucho más alta que la del punto crítico, el gas puede tomarse como ideal, y
si la presión es mucho más baja que la del punto crítico el gas también se
puede tomar como ideal.
Este factor es uno de los parámetros que, con mayor precisión

diferencia el comportamiento de los fluidos en estado liquido del estado
gaseoso. Indica la falta de idealidad de un gas y ademas define el
comportamiento de los gases a determinadas condiciones de presión y
27
temperatura. Lo que hace que sea de vital importancia para todos los
diseños e instalaciones que trabajan con fluidos compresibles.
Las ecuaciones aplicadas a gases ideales se cumplen solo para
temperaturas y presiones relativamente bajas, y con el objetivo de poder
utilizar la ecuación general de los gases a altas presiones y temperaturas, es
necesario agregarle a dicha ecuación el factor de compresibilidad Z,
obteniendo la ecuación general del estado gaseoso en condiciones reales.
(Pino, 2008)
(2.37)
𝑃𝑉 = 𝑍𝑛𝑅𝑇
En la ecuación (2.38) (Z) es el factor de compresibilidad y representa

un factor de desviación entre el comportamiento ideal de los gases y el
comportamiento real. Es un parámetro con el cual se determina el efecto de
comprimir un gas para llevarlo a sus condiciones reales de operación. Este
se representa a través de la siguiente ecuación:
𝑉(𝑟𝑒𝑎𝑙)
𝑍= (2.38)
𝑉(𝑖𝑑𝑒𝑎𝑙)
Es decir que Z expresa la relación entre un volumen real de un gas a

una determinada presión y temperatura de operación, con respecto al
volumen del mismo gas en condiciones ideales. Aplicando el principio de los
estados correspondientes, el cual establece que el factor Z es
aproximadamente el mismo para todos los gases si estos tienen la misma
presión y temperatura reducidas; la presión reducida Pr y la temperatura
reducida Tr se definen como:
𝑃 (2.39) 𝑇
𝑃𝑟 = 𝑇𝑟 = (2.40)
𝑃𝑐 𝑇𝑐
28
En estas ecuaciones deben utilizarse presiones y temperaturas

absolutas. Por lo tanto la presión y temperatura crítica de una sustancia se
usa para definir un estado reducido. En realidad todos los métodos sencillos
para calcular Z, se fundamentan en el principio de los estados
correspondientes, desarrollados por Van Der Waals, y que permite
determinar el factor de compresibilidad de los gases.
2.2.10 Significado del valor de Z
Si el valor de Z = 1 esto indica que el gas se comporta como ideal. Si el

valor de Z es mayor o menor que 1 el gas se comporta como un gas real.
Z > 1 gases como el hidrógeno y neón, que son difícilmente compresibles.
Z < 1 gases como el oxigeno, argón y metano, fácilmente compresibles.
Mientras más grande sea la desviación del valor de Z con respecto a 1,
mayor es la desviación respecto al comportamiento ideal del gas. (Pérez y
Martínez 1994).
2.2.11 Determinación del factor de compresibilidad Z
Hoy en día se disponen de muchos métodos para la determinación del

factor de compresibilidad Z de un gas, a continuación se nombran los más
comunes: Standing–Katz, Pitzer, Van Der Waals, Redlich–Kwong, Soave–
Redlich–Kwong, Peng–Robinson, Starling, Sarem, Berthelot, Wohl, Esferas
Duras, Hall Y Yarborough, Beattie-Bridgeman, Benedict–Webb–Rubin,
Refracción Molecular de Eykman, Dranchuk–Purvis–Robinson, Dranchuk–
Abou-Kassen, Carlile–Gillett, entre otros.
La selección del método a utilizar para determinar Z depende de los
requerimientos y las condiciones de operación, ya que cada uno de ellos
considera diferentes criterios para su aplicación. Con el avance de la
29
tecnología se han creado distintos programas que permiten agilizar los

cálculos de algunas de estas ecuaciones, por ejemplo, Gray–Sims, Papay y
Brill, lo que permite obtener valores de Z de manera más sencilla.
2.2.12 Análisis no paramétrico de la varianza
Los procedimientos inferenciales usuales denominados “métodos

paramétricos”, se fundamentan en el comportamiento o supuestos de la
distribución de la población subyacente de los datos. Por lo que pueden
hacerse ciertas suposiciones o aseveraciones de los parámetros que la
especifican.
En la práctica frecuentemente nos encontramos con datos, los cuales
por su naturaleza o por la manera de obtenerlo no son capaces de llenar los
requerimientos necesarios para analizarlos por métodos paramétricos. Para
analizar estos datos se hacen necesarios métodos que no requieran el
conocimiento previo acerca de la distribución de la población subyacente, es
decir, necesitamos la aplicación de “métodos libres de distribución”. A estos
métodos se le denomina no paramétricos. (Barrios y Martínez, 2011).
2.2.13 Prueba de Friedman
Esta prueba puede utilizarse en aquellas situaciones en las que se

seleccionan n grupos de k elementos de forma que los elementos de cada
grupo sean lo más parecidos posible entre sí, y a cada uno de los elementos
del grupo se le aplica uno de entre k ‘’tratamientos’’, o bien cuando a cada
uno de los elementos de una muestra de tamaño n se le aplican los k
‘’tratamientos’’.
30
La hipótesis nula que se contrasta es que las respuestas asociadas a

cada uno de los ‘’tratamientos’’ tienen la misma distribución de probabilidad o
distribuciones con la misma mediana, frente a la hipótesis alternativa de que
por lo menos la distribución de una de las respuestas difiere de las demás.
Para poder utilizar esta prueba las respuestas deben ser variables continuas
y estar medidas por lo menos en una escala ordinal. Los datos se disponen
en una tabla en la que en cada fila se recogen las respuestas de los k
elementos de cada grupo a los k tratamientos:
Tabla 2.2-Datos para la prueba de Friedman

Grupo\ Tratamiento 1 2 … Ἥ … K
1 x11 x12 … x1Ἥ … x1k
… …. …. … …. … ….
i xi1 xi2 … xij … xik
… … … … … … …
N xn1 xn2 … Xnj … Xnk
A las observaciones de cada fila se les asignan rangos de menor a

mayor desde 1 hasta k; a continuación se suman los rangos
correspondientes a cada columna, siendo RἭ la suma correspondiente a la
columna Ἥ-ésima. Si la hipótesis nula es cierta, la distribución de los rangos
en cada fila se debe al azar, y es de esperar que la suma de los rangos
correspondientes a cada columna sea aproximadamente igual a n (k + 1)/2.
La prueba de Friedman determina si las RἭ observadas difieren
significativamente del valor esperado bajo la hipótesis nula.
El estadístico de prueba es:
12
𝑋𝑟 2 = (∑ 𝑅𝑗 2 ) − 3𝑛(𝐾 + 1) (2.41)
𝑛𝑘 (𝑘+1)
31
Si la hipótesis nula Ho es cierta y el número de columnas y/o de filas es

2
moderadamente grande la distribución de Xr se aproxima a una chi-
cuadrado con k – 1 grados de libertad; de forma que se rechaza la hipótesis
nula para valores de X superior al valor crítico para el nivel de significación
fijado.
2.3. PROCEDIMIENTO METODOLÓGICO
2.3.1 DESCRIPCIÓN DE LOS MÉTODOS MÁS UTILIZADOS PARA LA

DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z
Esta etapa se desarrolló a través de la revisión de tesis, trabajos de

ascenso, investigación y documentos en internet, con el fin de obtener la
información requerida acerca de los métodos más utilizados para el cálculo
del factor de compresibilidad Z, tomando en consideración los criterios de
aplicación, rangos de uso y limitantes en cada caso.
Se tomaron en cuenta tanto los métodos basados en ecuaciones de

estado, así como los métodos gráficos empleados para el cálculo de Z. En
base a lo descrito anteriormente se describieron los siguientes métodos:
 Métodos basados en ecuaciones de estado: Van Der Waals, Peng-

Robinson, Hall–Yarborough, Redlich-Kwong, Soave-Redlich-
Kwong, Virial, Benedict-Webb-Rubín, Dranchuk-Abou-Kassen,
ajuste de Dranchuk-Purvis-Robinson.
 Métodos gráficos: Standing-Katz, Eykman, Sarem y Pitzer.
32
2.3.2 DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z A

MÉTODOS REPRESENTATIVOS PARA SU CÁLCULO
Para cumplir con este objetivo se realizó una serie de cálculos del factor
de comprensibilidad Z, mediante la selección de nueve (9) métodos de los
descritos en la etapa anterior, empleando los valores de presión, volumen y
temperatura, así como la composición del gas natural. Para la selección de
dichos métodos se tomó en cuenta el rango de aplicación, facilidad de
cálculo y frecuencia de uso en la industria petrolera. En base a ello se trabajó
con las siguientes ecuaciones de estado: Van Der Waals, Peng Robinson,
Redlich-Kwong, Soave-Redlich-Kwong, Dranchuk-Abou-Kassen, ajuste de
Dranchuk-Purvis-Robinson y los métodos gráficos: Standing-Katz, Eykman y
Pitzer. Por otra parte se utilizaron cuatro (4) muestras de gas provenientes de
cromatografías encontradas en tesis de investigación de la Universidad de
Oriente, núcleo Monagas. Para facilitar los cálculos se utilizó un programa
bajo el entorno de Excel realizado por los autores, de esta forma se
obtuvieron los valores de los parámetros requeridos para ser incluidos en la
ecuación general de cada uno de los métodos.
2.3.3 COMPARACIÓN A TRAVÉS DE UN ANÁLISIS ESTADÍSTICO DE

OBTENIDOS
Esta etapa se realizó a través de la aplicación del método estadístico

prueba de Friedman, siendo esta prueba un método no paramétrico para el
análisis de la varianza de doble entrada de rangos, el cual se desarrolló bajo
el entorno del programa estadístico Statgraphic. Esto permitió la comparación
33
de los valores del factor de compresibilidad Z obtenidos para cada método

establecido, con respecto al número de muestras que se seleccionaron.
2.4. DEFINICIÓN DE TÉRMINOS BÁSICOS
Compresibilidad: es una propiedad de la materia a la cual se debe que

todos los cuerpos disminuyan de volumen al someterlos a una presión o
compresión determinada manteniendo constantes otros parámetros. (Pino,
2008).
Cromatografía de gases: se define como un método de separación de los

diferentes componentes de una mezcla. Este método logra la separación de
los mismos a través del paso de una muestra por una fase estacionaria con
la ayuda de la fase móvil, cada componente de la muestra tiene propiedades
particulares que permitirá su interacción en forma diferente entre la fase
estacionaria y móvil. (Skoog y Leary, 2008).
Estadística: La estadística es la ciencia de los datos, la cual implica su

recolección, clasificación, síntesis, organización, análisis e interpretación,
para la toma de decisiones frente a la incertidumbre. (Ángel, 1998).
Hipótesis nula: es una hipótesis construida para anular o refutar, con el

objetivo de apoyar una hipótesis alternativa. Cuando se la utiliza, la hipótesis
nula se presume verdadera hasta que una evidencia estadística en la forma
de una prueba empírica de la hipótesis indique lo contrario. (Ángel, 1998).
Presión crítica: La presión crítica es, a temperatura crítica, la más alta

presión a la que una sustancia pura puede existir como líquido en equilibrio
con su gas. Además es la mínima presión que debe ejercerse para licuar un
gas a temperatura crítica. (Rojas, 2003).
34
Presión reducida: La presión reducida se define como la presión actual 𝑃

dividida por la presión crítica 𝑃𝑐 (Cengel, 2002).
Punto crítico: es aquel donde que las propiedades intensivas (aquellas que
no dependen de la masa: densidad, viscosidad, etc.) del gas y líquidos son
idénticas. (Rojas, 2003).
Temperatura crítica: Toda sustancia tiene una temperatura crítica, por

arriba de la cual la fase gaseosa no se puede licuar, independientemente de
la magnitud de la presión ejercida. Ésta es también, la temperatura más alta
a la que puede existir una sustancia en estado líquido. Esto quiere decir, que
por encima de la temperatura crítica no puede distinguirse un líquido de un
gas; se tiene simplemente, un fluido. (Rojas, 2003).
Temperatura reducida: La temperatura reducida se define como la

temperatura actual 𝑇 dividida por la temperatura crítica 𝑇𝑐. (Cengel, 2002).
CAPÍTULO III
DESARROLLO DE OBJETIVOS
3.1 DESCRIPCIÓN DE LOS MÉTODOS MÁS UTILIZADOS PARA LA

DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z
3.1.1 Métodos basados en ecuaciones de estado para el cálculo del

factor de compresibilidad Z
Ecuación virial
La ecuación virial de estado, es la única que tiene una base teórica

firme y se basa en el comportamiento intermolecular previsto por la mecánica
estadística.
Expresa el factor de compresibilidad como una serie de potencia del

inverso del volumen molar de la siguiente manera:
𝑃𝑣̃ 𝐵 𝐶 𝐷
Z= =1+ ̃ + + +⋯ (3.1)
𝑅𝑇 𝑣 𝑣̃ 2 𝑣̃ 3
En donde los coeficientes B, C, D… se llaman coeficientes viriales, los cuales

son independientes de la presión o de la densidad y para los componentes
puros ellos son funciones únicamente de la temperatura y del compuesto.
Estos coeficientes están relacionados directamente con las fuerzas
intermoleculares que existen entre grupos de moléculas; por ejemplo el
coeficiente B describe interacciones entre pares de moléculas, C entre
grupos de tres moléculas y así sucesivamente (González, 2006).
El factor de compresibilidad en algunas ocasiones se expresa como

una serie de potencias en presión (más fácil de usar):
35
36
𝑃𝑣̃
Z= = 1 + 𝐵′ 𝑃 + 𝐶 ′ 𝑃2 + 𝐷 ′ 𝑃3 + ⋯ (3.2)
𝑅𝑇
donde los coeficientes B’, C’, D’, dependen de la temperatura pero son
independientes de la presión o densidad. Para mezclas, sin embargo, estos
coeficientes dependen de la composición de una manera más complicada
que los de la ecuación (3.1). La relación entre los coeficientes de las
ecuaciones (3.1) y (3.2) son las siguientes:
′
𝐵 ′
𝑐 − 𝐵2 ′
𝐷 − 3𝐵𝐶 + 2𝐵3
𝐵 = 𝐶 = 𝐷 =
𝑅𝑇 (𝑅𝑇)2 (𝑅𝑇)3
Las dos formas de la ecuación virial son series infinitas; su uso práctico
en ingeniería resulta cuando converge rápidamente, es decir cuando no se
necesita más de dos o tres términos para que las series proporcionen valores
razonablemente aproximados. Esto es posible para gases y vapores a
presiones bajas o moderadas.
La ecuación virial truncada después del segundo término es de gran

utilidad y a densidades moderadas cerca de la mitad de la densidad crítica
proporciona muy buenos resultados.
𝐵
𝑍 =1+ (𝑣𝑖𝑟𝑖𝑎𝑙 𝑡𝑟𝑢𝑛𝑐𝑎𝑑𝑎) (3.3)
𝑣̃
𝐵𝑃
𝑍 = 1 + 𝐵′ 𝑃 = 1 + (3.4)
𝑅𝑇
Esta ecuación representa satisfactoriamente el comportamiento para

muchos vapores a temperatura sub crítica y hasta presiones de 15 bar. A
temperaturas mayores es apropiada para gases dentro de un intervalo de
presiones creciente a medida que aumenta la temperatura.
37
Las desventajas de usar la ecuación virial es la dificultad para calcular

los coeficientes sin recurrir a datos experimentales. Además por la forma que
presenta el ajuste a datos experimentales tampoco es tarea sencilla, por lo
tanto la ecuación virial es poco usada en ingeniería. Sin embargo, muchas de
las otras ecuaciones comúnmente utilizadas pueden ser reducidas a una
forma similar a la ecuación virial, ya que de esta forma las ecuaciones son
más útiles.
Ecuación de Benedict-Webb-Rubin (BWR)
Benedict-Webb-Rubin propusieron una ecuación de estado inspirada en

la ecuación virial la cual ha sido extremadamente valiosa en la correlación de
datos volumétricos, tanto líquidos como gaseosos para hidrocarburos livianos
y mezclas, con más aplicabilidad en el análisis de hidrocarburos ligeros,
debido a su exactitud y a su alto rango de acción.
El desarrollo de esta ecuación fue inicialmente empírico debido a que

los factores constantes que la componían se determinaron de esta forma,
pero posteriormente surgieron variaciones que permitieron un tratamiento
matemático más simple, ya que estas constantes se presentaron como
valores dependientes en función de otros factores como la densidad o el
volumen, a su vez estas extensiones de la ecuación permitieron reducir los
condicionantes de su aplicación y así se logro ampliar el rango de acción de
la misma hasta el posterior desarrollo de la ecuación de Lee-Kesler.
En esta ecuación, Ao, Bo, a, b y c son constantes para cada

componente puro. En el caso de mezclas, se han desarrollado diferentes
reglas de combinación para determinar las constantes respectivas. Tiene
38
ocho constantes para cada componente. En los constituyentes puros, las

constantes se encuentran definidas.
𝐶 𝑐𝜌 3 (1+Ƴ 𝜌 2)
𝑃 = 𝑅𝑇𝜌 + (𝐵𝑂 𝑅𝑇 − 𝐴𝑂 − 𝑇𝑂2 ) 𝜌2 + (𝑏𝑅𝑇 − 𝑎)𝜌3 + 𝑎 𝛼 𝜌6 + 2 (3.5)
𝑇 2𝑒 (Ƴ 𝜌 )
𝑃𝑉 𝐴𝑂 𝐶𝑂 𝑎 𝑎𝛼 𝑐𝜌2 (1+Ƴ𝜌2 )
= 𝑍 = 1 + (𝐵𝑂 − − ) 𝜌 + (𝑏 − ) 𝜌2 + + 2 (3.6)
𝑅𝑇 𝑅𝑇 𝑅𝑇 3 𝑅𝑇 𝑅𝑇 𝑅𝑇 3 𝑒 (Ƴ𝜌 )
Donde ρ = densidad molar y Ao, Bo, Co, a, b, c, α, Ƴ son constantes

para cada componente. El uso de la ecuación BWR requiere el empleo de
computadores digitales, en vista de la cantidad de cálculos que implica. Ha
sido muy usada para la determinación de fases (constantes de equilibrio),
factores de compresibilidad, Z, fugacidades, entalpía, entropía, etc. También
ha sido utilizada para cálculos a bajas temperaturas, variando los
coeficientes con temperatura.
Ecuación de Van Der Waals (VDW)
Para una descripción exacta del comportamiento PVT de los fluidos

sobre intervalos amplios de temperatura y presión, se requiere una ecuación
de estado más completa que la ecuación virial. Tal ecuación debe ser
aplicable tanto para líquidos, como gases y vapores. Aun así, no debe ser
tan compleja como para presentar grandes dificultades numéricas o
analíticas en su aplicación.
Las ecuaciones polinomiales que son cúbicas en el volumen molar

ofrecen un compromiso entre generalidad y simplicidad apropiado para
muchos fines. Las ecuaciones cúbicas son, de hecho, las ecuaciones más
sencillas capaces de representar el comportamiento de líquidos y vapores.
39
La primera ecuación cúbica de estado práctica, fue desarrollada por J.

D. Van Der Waals. Este derivo una ecuación capaz de describir la
continuidad entre las fases líquido y vapor. Propuesta en 1873, significó una
modificación de la ecuación de gases ideales, a través de la inclusión de dos
correcciones importantes con respecto a ese modelo, según el cual un fluido
está constituido por partículas puntuales que no interactúan entre sí. En
primer lugar, Van Der Waals establece que las moléculas tienen un tamaño
finito y real, de manera que parte del volumen no está disponible al
movimiento molecular. Esto aumenta el número de colisiones con las
paredes del recipiente que contiene el gas, lo que implica un aumento en la
presión (Malanowski y Anderko, 1992). En segundo lugar, existen
interacciones entre moléculas en forma de atracciones mutuas (cohesión),
cuyo efecto es una disminución de la presión de la sustancia. De esta
manera, llega a una expresión de la presión con dos contribuciones:
𝑃 = 𝑃𝑟𝑒𝑝𝑢𝑙𝑠𝑖ó𝑛 − 𝑃𝑎𝑡𝑟𝑎𝑐𝑐𝑖ó𝑛 (3.7)
Donde:
𝑅𝑇 (3.8)
𝑃𝑟𝑒𝑝𝑢𝑙𝑠𝑖ó𝑛 =
𝑣−𝑏
Esta presión de repulsión es equivalente a la expresión para un gas

ideal pero con el volumen disponible para el movimiento molecular corregido
a través del covolumen (cuatro veces el volumen propio de las moléculas), b
una constante característica de cada sustancia. El término atractivo por su
parte viene expresado como:
𝑎
𝑃𝑎𝑡𝑟𝑎𝑐𝑐𝑖ó𝑛 = (3.9)
𝑣2
40
Lo que supone que la disminución en la presión es proporcional al

número de moléculas por unidad de volumen e inversamente proporcional al
volumen.
La ecuación de Van Der Waals a pesar de que es muy superior a la ley

del gas ideal, siendo cualitativamente correcta, es incapaz de predecir con
exactitud las propiedades de las sustancias puras. Esto se debe a que tanto
el término repulsivo como el atractivo son teóricamente incorrectos. El
término repulsivo de la ecuación, es incapaz de presentar las repulsiones de
cuerpo rígido de manera adecuada, por lo que el término atractivo debe
compensar esta deficiencia.
Esta ecuación permite en forma relativamente fácil determinar el factor

de compresibilidad Z para los gases reales, su facilidad se sustenta por el
desarrollo de la computación, y la utilización de los modelos de simulación,
con los que se facilitan muchos de los cálculos de este parámetro
operacional; la ecuación (3.10) es utilizada para gases puros, en esta
ecuación (a) se denomina el parámetro de atracción y (b) el parámetro de
repulsión o el volumen molar efectivo.
𝑏𝑃 𝑎𝑃 𝑎𝑏𝑃 2
𝑍 3 − (𝑅𝑇 + 1) 𝑍 2 + ((𝑅𝑇)2) 𝑍 − ( 𝑅𝑇 3 ) = 0 (3.10)
Donde:
27𝑅2 𝑥𝑇𝑐 2 𝑅𝑥𝑇𝑐 8𝑃𝑐𝑥𝑉𝑐

𝑎= 𝑏= 𝑅=
64𝑃𝑐 8𝑃𝑐 3𝑇𝑐
Las constantes en esta ecuación de estado pueden evaluarse mediante

un ajuste de los datos PVT disponibles. Sin embargo, al tratarse de una
41
ecuación cúbica sencilla, las estimaciones adecuadas provienen de las

constantes críticas Tc y Pc.
Aunque esta ecuación tal vez no proporcione los mejores valores

posibles, ellos son razonables y casi siempre pueden determinarse, debido a
que a menudo se conocen las temperaturas y presiones críticas (en
contraste con los datos PVT extensos) o pueden estimarse de manera
confiable. Para mezclas es utilizada la ecuación (2.19).
Las constantes A y B pueden ser para sustancias simples o mezclas,
aunque a y b son conocidas como las constantes de Van Der Waals no son
realmente constantes, ya que varían de un gas a otro, es decir, que son
constantes para el gas que se trate.
Ecuación de Redlich-Kwong (RK)
Se deriva de la ecuación de Van Der Waals, fue introducida en 1949 y

pasó a ser la primera ecuación de estado cúbica que encontró una aplicación
industrial exitosa. Motivó además una gran cantidad de estudios que
exploraron su versatilidad y limitaciones para la determinación de
propiedades termodinámicas de compuestos puros y sus mezclas. Redlich y
Kwong aseguraban que su ecuación proporcionaba resultados
“satisfactorios” por encima de la temperatura crítica para cualquier presión
(Figueira, 2005).
No se utiliza con tanta frecuencia por la mayor dificultad en sus cálculos

y el uso general, ya que no da buenos resultados sobre la fase líquida y por
ello no puede usarse para calcular precisamente los equilibrios líquido-vapor.
Este método es empleado tanto para gases dulces como para gases
agrios. Para este último caso, cuando el gas contiene CO 2, y H2S o ambos,
42
Wichert y Azis proponen corregir la presión y la temperatura crítica de cada

componente de la mezcla, por medio de un factor de ajuste (Frk), en la
siguiente forma:
𝑇′𝑐𝑖 = 𝑇𝑐𝑖 − 𝐹𝑟𝑘 (3.11)
𝑃𝑐𝑖 − 𝑇′𝑐𝑖
𝑃′𝑐𝑖 = 𝑇𝑐𝑖
(3.12)
Donde: T’ci y P'ci son la temperatura y la presión críticas corregidas, T ci y Pci

son los respectivos valores sin corregir (de tablas).
El valor de Frk se calcula de la ecuación (3.13)
𝐹𝑟𝑘 = 27𝐴(1 − 𝐴) + 7,5(√𝐵 − 𝐵2 ) (3.13)
Donde: A = fracción molar de (CO2 + H2S) y B = fracción molar de H2S.
Para determinar el factor de compresibilidad (Z), este método se basa en la

solución de la siguiente ecuación:
𝑍 3 − 𝑍 2 + (𝐷2 𝑃 − 𝐸𝑃 − 𝐸 2 𝑃2 )𝑍 − (𝐷2 𝐸𝑃2 ) = 0 (3.14)

Donde:
𝑎𝑥𝑃 0,42748𝑃𝑅 𝑏𝑥𝑃 0,08664𝑃𝑅

𝐴=( = ) 𝐵= ( = )
𝑅2 𝑥𝑇𝑐2 𝑇𝑅2,5 𝑅𝑇 𝑇𝑅
𝑎 𝑏
𝐷= √ 𝐸=
𝑅2 𝑇 2,5 𝑅𝑇
P= Presión a condiciones de operación.
TR= Temperatura reducida.

43
R= Constante de los gases, en el Sistema Británico de Unidades tiene un

valor de 10,732 (lpca) (pie3)/(lbM )(ºR)
Los parámetros D y E son ecuaciones resultantes de la simplificación

de los parámetros A y B originalmente usado en la ecuación de Redlich-
Kwong, lo cual permite facilitar los cálculos.
Actualmente se conoce que la ecuación tiene limitaciones en la

predicción de propiedades volumétricas de los líquidos, debido a que su
factor de compresibilidad crítico Zc = 0,333 es demasiado alto y por lo tanto
no es muy exacta cerca del punto crítico. Sin embargo, es mucho más exacta
que la ecuación de Van Der Waals y la ecuación de los gases ideales,
además da buenos resultados a presiones altas y temperaturas sobre la
temperatura crítica. Por lo tanto esta ecuación es muy útil, ya que combina la
simplicidad de una ecuación de dos parámetros con un alto grado de
exactitud.
Ecuación de Soave Redlich-Kwong (SRK)
La Ecuación de Soave-Redlich-Kwong (SRK), propuesta en 1972, fue

la primera modificación de la forma simple de la ecuación de Redlich-Kwong
(RK). Es apropiada para el cálculo del equilibrio líquido-vapor, pero su rango
de aplicabilidad es más limitado y menos confiable en sistemas no ideales
que las ecuaciones de Peng-Robinson y Lee-Kesler-Plocker. Su uso típico es
para refinación, petroquímica y procesamiento de gas. Rango de aplicación:
todos (menos exacto cerca de los valores de las propiedades críticas).
Soave, al modificar la ecuación de estado Redlich-Kwong obtuvo el siguiente

resultado:
44
𝑅𝑇 𝑎𝑐 𝛼(𝑇𝑟 ,𝑤)
𝑃= − (3.15)
𝑣−𝑏 𝑣(𝑣+𝑏)
Donde:
1 𝑅2 𝑇𝑐2
𝑎𝑐 = 3 − 2
9(√2 − 1) 𝑃𝑐
𝑇
√𝛼 = 1 + 𝑚 [1 − √ ]
𝑇𝑐
𝑚 = 0,48 + 1,574𝑤 − 0,176𝑤 2 (Para 0 ≤ w ≤ 0,5)
Esta ecuación, al introducir otro parámetro, permite mejorar aún más los
resultados entregados por la ecuación de Redlich-Kwong, sin embargo aún
no logra representar fielmente el comportamiento de un fluido para grandes
intervalos de presión y temperatura. Para determinar el factor de
compresibilidad (Z) se usa la ecuación (2.26)
La ecuación de SRK tuvo una rápida aceptación y se convirtió en una

de las más populares en la industria de los hidrocarburos. Sin embargo,
continuó el interés en mejorar la capacidad de predicción de la ecuación de
estado de diferentes propiedades de sustancias puras y mezclas. Se
propusieran diversas modificaciones, no solo de la función de cohesión, sino
también de la funcionalidad en volumen del término atractivo (Figueira,
2005).
45
Ecuación de Peng- Robinson (PR)
Otra modificación de la ecuación de estado de Redlich-Kwong fue la

publicada por Peng-Robinson (1976-1977), quienes desarrollaron una nueva
ecuación ante la necesidad de obtener predicciones más precisas del
comportamiento volumétrico de las fases coexistentes en el cálculo del
equilibrio líquido-vapor, puesto que la ecuación SRK, generaba en muchos
casos volúmenes de líquido mayores que los reales, como consecuencia de
su alto factor de compresibilidad crítico.
De esta manera, propusieron una nueva funcionalidad en volumen del
término atractivo de la ecuación de estado bajo la idea de que la selección de
una funcionalidad adecuada puede hacer que el factor de compresibilidad
crítico predicho por la ecuación de estado sea más cercano a sus valores
reales.
𝑅𝑇 𝑎𝑐 ∝(𝑇𝑟 ,𝑤)
𝑃= − (3.16)
𝑣−𝑏 𝑣(𝑣+𝑏)+𝑏(𝑣−𝑏)
Generalmente la ecuación de Peng Robinson da unos resultados

similares a la de Soave, aunque es bastante mejor para predecir las
densidades de muchos compuestos en fase líquida.
Es una de las ecuaciones de mayor aplicabilidad para determinar Z,

debido al poco tiempo de computación que emplea. Está representada de la
siguiente manera:
𝑍 3 − (1 − 𝐵)𝑍 2 + (𝐴 − 3𝐵2 − 2𝐵)𝑍 − (𝐴𝐵 − 𝐵2 − 𝐵3 ) = 0
Donde:
46
𝑎𝑥𝑃 𝑏𝑥𝑃
𝐴= (3.17) 𝐵= (3.18)
𝑅2 𝑥𝑇 2 𝑅𝑇
𝑅2 𝑥𝑇𝐶2 𝑅𝑥𝑇𝑐
𝑎 = 0,45724 ( )𝛼 𝑏 = 0,0778 ( )
𝑃𝐶 𝑃𝑐
𝑇
√𝛼 = 1 + 𝑚 [1 − √ ]
𝑇𝑐
𝑚 = 0,37464 + 1,54226𝑤 − 0,26992𝑤 2
Los parámetros A y B se determinan por las ecuaciones (3.17) y (3.18)
Donde:
√𝛼 = Mejora la correlación de la presión de vapor de puros.
w = Factor acéntrico de Pitzer.
Tc = Temperatura crítica. (Adimensional)
Pc = Presión crítica. (Adimensional)
Muchas ecuaciones modificadas basadas en el modelo básico de Peng-

Robinson están disponibles con rangos extendidos de aplicación no ideal,
incluyendo sistemas ácidos. Su uso típico es en la refinación, petroquímica y
procesamiento de gas. Aplica para cualquier rango, sin embargo es menos
exacto cerca de los valores de las propiedades críticas.
47
Ecuación de Hall-Yarborough
Hall y Yarborough basados en la ecuación de estado de Starling-

Carnaban desarrollaron las siguientes ecuaciones para el cálculo de Z:
2
0.06125𝑃𝑠𝑟 𝑡𝑒 −1,2(1−𝑡) (3.19)
𝑍=
𝑦
Donde:
t: recíproco de la temperatura pseudorreducida, t = Tsc/T
y: densidad reducida, la cual se obtiene a partir de la siguiente ecuación:
𝑦+ 𝑦 2 + 𝑦 3−𝑦 4
−𝐴𝑃𝑠𝑟 + (1−𝑦)3
− 𝐵𝑦 2 + 𝐶𝑦 𝐷 = 0 (3.20)
Donde:
2
𝐴 = 0.06125𝑡𝑒 −1,2(1−𝑡) 𝐶 = 90.7𝑡 − 242.2𝑡 2 + 42.4𝑡 3
𝐵 = 14.76𝑡 − 9.76𝑡 2 + 4.58𝑡 3 𝐷 = 2.18 + 2.82𝑡
Esta es una ecuación no lineal, su solución requiere de la aplicación de

técnicas numéricas y métodos de aceleración de convergencia como el
método de Newton –Raphson.
Para obtener buenos resultados con este método, se recomienda que la

Psr y Tsr estén entre los siguientes rangos:
0.1 ≤ Psr ≤ 24.0 y 1.2 ≤ Tsr ≤ 3.0

48
Para gases naturales con impurezas se recomienda corregir P scM y

TscM por el método de Wichert y Aziz.
La ecuación de Hall- Yarborough tiene una precisión del 0.3 % para

temperaturas reducidas mayores a 1.0
Ecuación de Dranchuk, Purvis y Robinson (DPR)
Estos autores hicieron un ajuste de la ecuación de estado BWR para

correlacionar la grafica de Standing-Katz. Escribieron la ecuación en la
siguiente forma:
𝐴 𝐴 𝐴 𝐴5 𝐴6 𝜌𝑟5 𝜌2 2
𝑍 = 1 + (𝐴1 + 𝑇 2 + 𝑇33 ) 𝜌𝑟 + (𝐴4 + 𝑇 5 ) 𝜌𝑟2 + 𝑇𝑠𝑟
+ 𝐴7 (1 + 𝐴8 𝜌𝑟2 ) 𝑇3𝑟 𝑒 −𝐴8𝜌𝑟 (3.21)
𝑠𝑟 𝑠𝑟 𝑠𝑟 𝑠𝑟
La densidad reducida ρr se calcula con la siguiente expresión (obtenida de la

ley de los gases):
𝜌 𝑍𝑐 𝑃𝑠𝑟 𝑃𝑠𝑟 (3.22)

𝜌𝑟 = = = 0.27
𝜌𝑐 𝑍𝑇𝑠𝑟 𝑍𝑇𝑠𝑟
En esta ecuación se tomó como factor de compresibilidad en el punto

crítico Zc = 0,27, considerado como un valor apropiado para mezclas
compuestas principalmente por metano. Además, determinaron los
siguientes valores para las constantes A1-A8 como se muestra a
continuación:
A1 = 0,31506237 A5= -0,61232032

A2 = -1,0467099 A6 = -0,10488813
A3 = -0,57832729 A7¨= 0,68157001
A4 = 0,53530771 A8 = 0,68446549
49
Esta correlación reproduce con buena aproximación la grafica de

Standing-Katz. Los rangos de aplicación del método son:
0,2 ≤ Psr ≤ 30,0 y 1,05 ≤ Tsr ≤ 3,0
Para gases naturales con impurezas se recomienda corregir P scM y

TscM por el método de Wichert y Aziz.
Ecuación de Dranchuk-Abou-Kassen (DAK)
Este método es muy similar al método anterior de DPR. En lugar de la

ecuación de BWR, utilizaron la ecuación de estado de Starling la cual
escribieron en la siguiente forma:
𝐴 𝐴 𝐴 𝐴 𝐴 𝐴 𝐴 𝐴
𝑍 = 1 + (𝐴1 + 𝑇 2 + 𝑇33 + 𝑇44 + 𝑇55 ) 𝜌𝑟 + (𝐴6 + 𝑇 7 + 𝑇28 ) 𝜌𝑟2 − 𝐴9 (𝑇 7 + 𝑇28 ) 𝜌𝑟3 +
𝑠𝑟 𝑠𝑟 𝑠𝑟 𝑠𝑟 𝑠𝑟 𝑠𝑟 𝑠𝑟 𝑠𝑟
𝜌𝑟2 2
𝐴10 (1 + 𝐴11 𝜌𝑟2 ) 3 𝑒 −𝐴11𝜌𝑟 (3.23)
𝑇𝑠𝑟
Igualmente, emplearon la ecuación (3.22) para la densidad reducida y

por un procedimiento similar al método anterior, determinaron los siguientes
valores de las constantes A1-A11:
A1 = 0,3255 A2 = -1,07 A3 = -0,5339
A4 =0,01569 A5 = -0.05165 A6 = 0.5475
A7 = -0,7361 A8 = 0,1844 A9 = 0,1056
A10 = 0,6134 A11= 0,721

50
Los rangos de aplicación del método son:
0,2 ≤ Psr ≤ 30,0 y 1,0 ≤ Tsr ≤ 3,0 y para P < 1,0 con 0,7 < Tsr ≤ 1,0
Se efectuó un estudio comparativo de los métodos de Y-H, DPR y

DAK y se concluyó que estos tres métodos producen resultados
satisfactorios en el cálculo de Z, en las regiones:
0,2 ≤ Pr ≤ 30
1,0 < Tr ≤ 3,0
Pr < 1,0; 0,7 < Tr ≤ 1,0
Pero conducen a resultados inaceptables en la región:
Tr = 1,0 y Pr ≥ 1,0. Sin embargo, para Pr < 1,0 y 0,7 < Tr ≤ 1,0 se
recomienda el método DAK.
3.1.2 Métodos gráficos para el cálculo del factor de compresibilidad Z
Método de Standing-Katz
Este método es probablemente el más popular para calcular el factor de

compresibilidad de gases dulces y que contengan pequeñas cantidades de
gases no hidrocarburos. El método se basa en la figura 1B (ver apéndice B),
donde Z está representado gráficamente en función de la temperatura y
presión reducida o más comúnmente pseudorreducidas. Estos valores se
calculan de diferentes maneras, de acuerdo con la forma de combinar la
composición y propiedades críticas de los componentes puros de un sistema,
la figura 1B puede emplearse hasta 10.000,0 lpca, con desviaciones del 1%.
51
El método Standing-Katz para gases dulces ha sido modificado con el

fin de aplicarlo a gases agrios (gases que contienen H 2S y/o CO2), Wichert y
Aziz utilizaron rangos de temperatura entre 40°F y 300ºF, presiones variables
entre 0 y 7000 lpcm y gases ácidos (CO2 y H2S) con una concentración del
80% o menos. Los valores de Psc y Tsc se calculan por la regla de Kay y se
modifican por un factor de ajuste (Fsk), calculado con:
𝐹𝑠𝑘 = 120(𝐴0.9 − 𝐴1.6 ) + 15(√𝐵 − 𝐵4 )
La temperatura y la presión pseudorreducidas se ajustan en la forma

siguiente:
𝑃𝑐𝑖 − 𝑇′𝑐𝑖
𝑇′𝑠𝑐 = 𝑇𝑠𝑐 − 𝐹𝑠𝑘 𝑃′𝑐𝑖 =
𝑇𝑐𝑖
Con estos valores, se calculan la presión y temperatura pseudoreducidas. El

factor Z se calcula luego de la figura 1B. (ver Apéndice B)
Donde:
A = fracción molar de (CO2 + H2S) y B= fracción molar de H2S

Fsk = factor de ajuste para Standing y Katz, ºR
Psc = Presión pseudocrítica de la mezcla, lpca
Tsc = Temperatura pseudocrítica de la mezcla, ºR
Método de Refracción Molecular de Eykman (EMR)
La ecuación de Eykman de Refracción Molecular es:

52
(𝑛 2 −1)𝑀
𝐸𝑀𝑅 = (𝑛+0,4)𝜌 (3.24)
Donde, EMR = refracción molecular de Eykman y n = índice de refracción del

gas o líquido, usando en el refractómetro luz amarilla de la línea D del sodio.
Para hidrocarburos parafínicos normales se ha encontrado, que:
𝐸𝑀𝑅 = 2,4079 + 0,7293𝑀 + 0,00003268𝑀2 (3.25)
Esta ecuación presenta una correlación bastante exacta y puede usarse

para determinar el EMR de la fracción más pesada de un sistema, a partir de
su peso molecular. Si se conoce el peso molecular y la gravedad específica
de la fracción más pesada, este último valor es generalmente el más exacto.
Para una serie homóloga, la relación entre el índice de Refracción Molecular
de Eykman, (EMRI = EMR/M) y la densidad (ρ) es lineal.
Cuando se usa el EMR para obtener compresibilidad del gas, es posible

utilizar los valores críticos obtenidos con base en la refracción molecular de
Eykman, en lugar de los pseudovalores.
I. Para CO2, H2S, N2 y C1 (Grupo I):
(𝐴)𝑇𝑐 (3.26)
= 0,4238 + 0,005117𝐸𝑀𝑅
𝑃𝑐
II. Para los demás hidrocarburos, C2 (Grupo II):
(𝐴)𝑇𝑐
= 0,1419 + 0,02437𝐸𝑀𝑅 + 0,00007911𝐸𝑀𝑅2 (3.27)
𝑃𝑐
53
Los valores A y B se utilizan cuando se esté usando el sistema métrico. Por

ello cada uno de estos parámetros es igual a la unidad en el sistema ingles.
Los siguientes pasos resumen el método para determinar el factor Z:
a. Se divide la composición en dos grupos:

Grupo I: CO2, H2S, N2 y C1
Grupo II: C2 y demás hidrocarburos
b. Si se encuentra presente C7+, calcular su EMR a partir de la ecuación

(3.25)
c. Para cada grupo del Paso a, se calcula:
𝐸𝑀𝑅𝐼 = ∑ 𝑦𝑖 ∗ 𝐸𝑀𝑅𝑖 (3.28)
Donde yi* es la fracción molar de cada componente en el Grupo I
𝐸𝑀𝑅2 = ∑ 𝑦𝑖𝑖 ∗ 𝐸𝑀𝑅𝑖𝑖 (3.29)
Donde yii* es la fracción molar de cada componente en el Grupo II
d. Con EMR1 se obtiene [(𝐴)𝑇𝑐/𝑃𝑐 ]1 para el Grupo I y con EMR2 se

obtiene [(𝐴)𝑇𝑐/𝑃𝑐 ]2 para el grupo II. Luego se calcula (𝐴)𝑇𝑐/𝑃𝑐 por:
𝑇𝐶 𝑇𝑐 𝑇𝑐 (3.30)
= [(∑ 𝑦𝑖 ) ( )] + [(∑ 𝑦𝑖 ) ( )]
𝑃𝐶 𝑃𝑐 1 𝑃𝑐 2
54
e. Se calcula el EMR de la mezcla:
𝑛
(3.31)
𝐸𝑀𝑅𝑚 = ∑ 𝑦𝑖 𝐸𝑀𝑅𝑖
𝑖=1
f. Con el valor anterior, se obtiene (𝐵)𝑇𝑐/√(𝑃𝑐 )
g. Con los valores de (𝐴)𝑇𝑐/𝑃𝑐 y (𝐵)𝑇𝑐/√(𝑃𝑐 ) se calcula Tc y P
h. Se calcula Pr = P/Pc y Tr = T/Tc
i. Del gráfico de la figura 2B (ver apéndice B) se determina Z.
Método de Sarem
Este método para la determinación de Z se basa en los polinomios de

Legendre de gado 0 a 5. La ecuación básica es:
6 6
(3.32)
𝑍 = ∑ ∑ 𝐴𝑖𝑗 𝑃𝑖 (𝑥 ) 𝑃𝑗 (𝑦)
𝑖=1 𝑗=1
Donde:
2𝑃𝑟 − 15 (𝑃⁄𝑃𝑐 ) − 7,5

𝑥= =
14,8 7,4
2𝑇𝑟 − 4 (𝑇⁄𝑇𝑐 ) − 2
𝑦= =
1,9 0,95
55
Los Polinomios de Legendre de grado 0 - 5, Pi y Pj son los siguientes:
P1(a) = 0,7071068
P2(a) = 1,224745 (a)
P3 (a) = 0,7905695 (3a2 - 1)
P4 (a) = 0,9354145 (5a3 - 3a)
P5(a) = 0,265165 (35a4 - 30a2 + 3)
P6(a) = 0,293151 (63a5 - 70a3 15a)
Donde "a" se reemplaza por "x" e "y", al efectuar la sumatoria de la

ecuación (3.32). La constante Aij, se obtiene de la tabla 3.1
Tabla 3.1-Valores de los coeficientes Aij usados en la Ec. 3.32
i j=1 j=2 j=3 j=4 j=5 j=6

1 2,1433504 0,0831762 -0.0214670 -0,0008714 0,0042846 -0,0016595
2 0,3312352 -0,1340361 0,0668810 -0,0271743 0,0088512 -0,0021521
3 0,1057287 -0,0503937 0,0050925 0,0105513 -0,0073182 0,0026960
4 -0,0521840 0,0443121 -0,0193294 0,0058973 0,0015367 -0,0028327
5 0,0197040 -0,0263834 0,0192621 -0,0115354 0,0042910 -0,0008130
6 0,0053096 0,0089178 -0,0108948 0,0095594 -0,0060114 0,0031175

Fuente: Pérez y Martínez, 1994
El método de Sarem puede ser utilizado cuando el gas natural contiene

impurezas corrigiendo por Wichert y Aziz.
Método de Pitzer
En este caso, el factor de compresibilidad se calcula con la ecuación:

56
𝑍 = 𝑍 0 + 𝑤. 𝑍 ʹ (3.33)
Donde Z0 y Z’ son funciones de la presión y la temperatura reducidas,

Pr y Tr, del sistema y pueden obtenerse de las figuras 3B y 4B (ver apéndice
B), donde "w" es el denominado factor acéntrico de Pitzer.
El factor acéntrico de Pitzer está definido por la siguiente ecuación:
𝑤 = −(𝑙𝑜𝑔𝑃𝑟0 + 1)𝑇𝑟 = 0,7
Donde Pr0 = presión de vapor reducida = P0/Pc para una temperatura

reducida, Tr = T/Tc = 0,7
En este método se establece que Z no es función exclusiva de T r y Pr;

se requiere más información, por esta razón incluyeron el factor acéntrico en
su ecuación. El factor acéntrico de cualquier sistema puede considerarse
como una medida de la desviación del comportamiento de un fluido simple, o
sea, la desviación del sistema del principio de estados correspondientes.
Para un fluido simple, el factor acéntrico w = 0.
Anteriormente se han discutido los métodos para encontrar los valores

de Psc y Tsc, de acuerdo con diferentes reglas de combinación de
propiedades y composición. En el método de Pitzer, estos valores se
calculan según la regla de combinación de Stewart Burkhardt-Voo (SBV), y
se ajustan en la forma sugerida por Wichert y Aziz, a saber:
𝑇𝑠𝑐 1 𝑇𝑐𝑖 2 𝑇𝑐𝑖 (3.34)

𝐽= = ∑ 𝑦𝑖 𝑥 + ∑ 𝑦𝑖 𝑥 √
𝑃𝑠𝑐 3 𝑃𝑐𝑖 3 𝑃𝑐𝑖
57
𝑇𝑐𝑖 𝑇𝑐𝑖
𝐾= = ∑ 𝑦𝑖 𝑥 (3.35)
√𝑃𝑐𝑖 √𝑃𝑐𝑖
Por lo tanto, se puede deducir que:
𝐾2 𝑇𝑠𝑐
𝑇𝑠𝑐 = 𝑃𝑠𝑐 =
𝐽 𝐽
𝑃𝑠𝑐 − 𝑇′𝑠𝑐
𝑇′𝑠𝑐 = 𝑇𝑠𝑐 − 𝐹𝑝 𝑃′𝑠𝑐 =
𝑇𝑠𝑐
Fp es el denominado factor de ajuste de Pitzer y está dado por:
𝐹𝑝 = 950𝑤(𝐴 − 𝐴3 ) + 9(√𝐵 − 𝐵2 )
Donde:
A = fracción molar de (CO2 + H2S)
B = fracción molar de H2S.
El método de Pitzer para determinar Z, aunque es muy exacto en la

forma descrita, no es muy apropiado para su uso con computadores, debido
a la dificultad para programar los valore de Z 0 y Z’.
3.2 DETERMINACIÓN DEL FACTOR DE COMPRESIBILIDAD Z A

MÉTODOS REPRESENTATIVOS PARA SU CÁLCULO
Considerando cuatro (4) muestras de gas natural obtenidas de dos (2)

tesis de grado de la Universidad de Oriente Núcleo Monagas, y al tener
definidos nueve (9) métodos basados en ecuaciones de estado y métodos
58
gráficos, permitió obtener cuatro (4) valores de Z para cada uno de los
métodos seleccionados. Ahora bien, al usar los datos de las cuatro (4)
cromatografías mencionadas para determinar Z, las operaciones de cálculo
se vuelven repetitivas, por esta razón se detalla a continuación una
descripción del cálculo que se realizó, haciendo referencia a la muestra 1A
(ver apéndice A) para cada método definido.
3.2.1 Ecuación de Van Der Waals
Para determinar el factor de compresibilidad Z utilizando los datos de

presión, temperatura y composición del gas de la muestra 1A (ver apéndice
A), se utilizó un programa de hoja de cálculo bajo el entorno de Microsoft
Excel realizado por los autores, lo que permitió desarrollar las operaciones
de una forma más sencilla.
En la figura 1C (ver apéndice C), considerando los valores de la

composición del gas de la muestra 1A (ver apéndice A), se observan los
datos de operación de dicha muestra (presión y temperatura) y la constante
general de los gases (R=10.732 Psia ft3 / lbM °R). Con estos datos se
procedió al cálculo de las variables A y B, cuyos resultados se aprecian en
las casillas de color verde. Posteriormente se desarrolló la ecuación cúbica
que requiere este método para el cálculo del factor de compresibilidad,
sustituyendo los valores de las constantes A y B en la ecuación general de
Van Der Waals para mezclas (ecuación 2.19) y obteniendo como resultado
Z=0.7835, que se observa en la celda de color amarillo.
Los resultados de Z con los valores de composición de las tres (3)

muestras de gas restantes para este método se muestran a continuación:
59
Tabla 3.2- Resultados del factor de compresibilidad Z utilizando la ecuación de

estado de Van Der Waals
Método Muestra 1 Muestra 2 Muestra 3 Muestra 4

Van Der
Waals 0.7835 0.8015 0.9192 0.9207
Fuente: Los autores
Nótese que para la muestra 1 y 2, provenientes de los reportes del

cromatografo ACOGAS (ver apéndice A), se aprecia una diferencia
significativa con respecto a las muestras 3 y 4, provenientes de un análisis
PVT del campo Santa Bárbara (ver apéndice A). Esto se debió básicamente
a las condiciones de operación de cada una de las cromatografías, a saber
presión y temperatura. Del mismo modo se observó un comportamiento
ascendente del valor de Z para cada una de las muestras seleccionadas.
Este comportamiento se repite para los métodos de Peng Robinson, Redlich-
Kwong y Soave-Redlich-Kwong respectivamente.
3.2.2 Ecuación de Peng Robinson

procedió al cálculo de las variables A y B, cuyos resultados se aprecian en
sustituyendo los valores de las constantes A y B en la ecuación general de
60
Peng Robinson para mezclas (ecuación 2.27) y obteniendo como resultado

Z=0.7183, que se aprecia en la celda de color amarillo.
Los resultados de Z con los valores de composición de las tres (3)


estado de Peng Robinson

Peng
Robinson 0.7183 0.7387 0.8833 0.8853
Fuente: Los autores
3.2.3 Ecuación de Redlich-Kwong

procedió al cálculo de las variables D y E, cuyos resultados se aprecian en
sustituyendo los valores de las constantes D y E, así como la presión de
operación en unidades absolutas (lpca), en la ecuación general de Redlich-
Kwong para mezclas (ecuación 2.25) y obteniendo como resultado Z=0.8093,
que se aprecia en la celda de color amarillo. Los resultados de Z con los
valores de composición de las tres (3) muestras de gas restantes para este
método se muestran a continuación:
61

estado de Redlich-Kwong

Redlich-
Kwong 0.8093 0.8270 0.9252 0.9270
Fuente: Los autores
3.2.4 Ecuación de Soave-Redlich-Kwong

datos de operación de dicha muestra (presión y temperatura), la constante
general de los gases (R=10.732 Psia ft3 / lbM °R) y los valores obtenidos de
las variables a y b (resaltados en color morado) que son requeridos por este
método. Con estos datos se procedió al cálculo de las variables A y B, cuyos
resultados se aprecian en las casillas de color verde. Posteriormente se
desarrolló la ecuación cúbica que requiere este método para el cálculo del
factor de compresibilidad, sustituyendo los valores de las variables A y B, en
la ecuación general de Soave-Redlich-Kwong para mezclas (ecuación 2.26) y
obteniendo como resultado Z=0.8638, que se aprecia en la celda de color
amarillo. Los resultados de Z con los valores de composición de las tres (3)

estado de Soave-Redlich-Kwong

Soave-Redlich-
Kwong 0.8638 0.8723 0.9615 0.9631
Fuente: Los autores
62
3.2.5 Ecuación de Dranchuk Abou Kassen

datos de operación de dicha muestra (presión y temperatura), al igual que el
valor de Z calculado por Standing y Katz que se aprecia en la celda de color
verde y es requerido para el cálculo de la densidad reducida que utiliza este
método para su desarrollo. La teoría indica que se requiere de corrección por
impurezas cuando están presentes componentes ácidos, por lo que se utilizó
el método de Wichert y Azís. Del mismo modo se muestran en las celdas de
color morado, las once (11) constantes que utiliza este ajuste y su respectiva
agrupación para finalmente obtener el resultado de Z=0.8433, que resalta en
la celda de color amarillo. Los resultados de Z con los valores de
composición de las tres (3) muestras de gas restantes para este método se
muestran a continuación:

estado de Dranchuk Abou Kassen

Dranchuk
Abou Kassen 0.8433 0.8471 0.8564 0.8575
Fuente: Los autores
Nótese que para este método siguió manteniéndose la tendencia de

aumento del factor de compresibilidad para cada una de las muestras, pero
esta vez la diferencia entre cada uno de los valores no es tan significativa
como se observó en las ecuaciones de estado estudiadas anteriormente.
Esto se debió a que este ajuste ademas de los datos de operación (presión y
63
temperatura), toma en cuenta otras variables para el cálculo de Z, a saber las

once (11) constantes propias del método y la densidad reducida que
depende tanto del Z hallado por Standing, así como también de la presión y
temperatura pseudoreducidas respectivamente.
3.2.6 Ecuación de Dranchuk-Pulvis-Robinson

datos de operación de dicha muestra (presión y temperatura), al igual que el
valor de Z calculado por Standing y Katz que se aprecia en la celda de color
verde y es requerido para el cálculo de la densidad reducida que utiliza este
método para su desarrollo. La teoría indica que se requiere de corrección por
impurezas cuando están presentes componentes ácidos, por lo que se utilizó
el método de Wichert y Azís. Del mismo modo se muestran en las celdas de
color morado, las ocho (8) constantes que utiliza este ajuste. Finalmente al
desarrollar la formula se obtiene el resultado de Z=0.8778, que resalta en
color amarillo. Los resultados de Z con los valores de composición de las tres
(3) muestras de gas restantes para este método se muestran a continuación:

estado de Dranchuk-Pulvis-Robinson

Dranchuk-
Pulvis- 0.8778 0.8846 0.9541 0.9549
Robinson
Fuente: Los autores
64
Nótese que si se aprecian diferencias significativas entre las muestras

1 y 2 provenientes de los reportes del cromatografo ACOGAS (ver apéndice
A), con respecto a las muestras 3 y 4, provenientes de un análisis PVT del
campo Santa Bárbara (ver apéndice A), si lo comparamos con el método
anterior. Se mantuvo el comportamiento ascendente del valor de Z.
3.2.7 Método de Standing-Katz

datos de operación de dicha muestra (presión y temperatura). Al aplicar este
método gráfico es necesario obtener primeramente los valores de la presión
pseudorreducida (Psr) y temperatura pseudorreducida (Tsr) respectivamente,
considerando que en la composición del gas están presente impurezas se
aplicó el método de Wichert y Azis para hacer las correcciones necesarias.
Los valores de Psr y Tsr resaltan en la celda color turquesa, el valor de
Z=0.84 es obtenido de la figura 1B (ver apéndice B) y se aprecia en la celda
de color amarillo. Los resultados de Z con los valores de composición de las
tres (3) muestras de gas restantes para este método se muestran a
continuación:

gráfico de Standing y Katz

Standing-Katz 0.84 0.87 0.95 0.96
Fuente: Los autores
65
Nótese que se mantuvo la tendencia ascendente de los valores de Z.

pero existe una diferencia considerable entre los resultados para la muestra
1 y 2 provenientes de los reportes del cromatografo ACOGAS (ver apéndice
A), con respecto a las muestras 3 y 4, provenientes de un análisis PVT del
campo Santa Bárbara (ver apéndice A). Es importante resaltar que a pesar
de que este método es uno de los más utilizados para el cálculo de Z, tiene
un porcentaje de error considerable debido a que el resultado se toma de
una gráfica, por lo que tiende a tener errores de lectura o apreciación. Este
comportamiento se repitió para el método gráfico de Eykman descrito a
continuación.
3.2.8 Método de Refracción Molecular de Eykman

datos de operación de dicha muestra (presión y temperatura). Para aplicar
este método gráfico fue necesario obtener primeramente los valores de la
presión crítica (Pc) y temperatura crítica (Tc) respectivamente. Ambos
valores se observan en las celdas color verde y dependen de las variables
X1, X2, X3, X4 presentes en la tabla. Posteriormente se calcularon los valores
de temperatura reducida (Tr) y presión reducida (Pr) mostrados en las celdas
color turquesa. Finalmente se procedió a leer de la figura 2B (ver apéndice B)
el valor de Z=0.82 representado en la celda de color amarillo.
Si la cromatografía presenta componentes pesados (C 3+, C5+, C7+)

como es el caso de las muestras tres (3) y cuatro (4), provenientes de un
análisis PVT del campo Santa Bárbara (ver apéndice A), la hoja de cálculo
permite calcular su EMR, conocido previamente su peso molecular. Esto se
muestra en las celdas color morado. Los resultados de Z con los valores de
66

gráfico de Eykman

Eykman 0.82 0.84 0.93 0.94
Fuente: Los autores
3.2.9 Método de Pitzer

datos de operación de dicha muestra (presión y temperatura). Para aplicar
este método gráfico fue necesario obtener primeramente los valores de a, k y
b que se obtuvieron de la tabla y están representados en color rojo. Esto se
hizo con el fin de aplicar la regla de combinación de Stewart-Burkhardt-Voo
(SBV) para mezclas y determinar el valor de los parámetros J, temperatura
pseudocritica sin corregir y presión pseudocritica sin corregir
respectivamente. Estos resultados se muestran en las celdas color morado.
Considerando que en la composición del gas están presentes impurezas se
aplicó el método de Wichert y Azis, pero en este caso con el factor de ajuste
de Pitzer (Fp), para hacer las correcciones necesarias. Una vez hallados los
valores de temperatura reducida y presión reducida que se aprecian en las
celdas color azul claro, se ubicaron los valores de Z 0= 0.85 y Z’=0.17 en las
figuras 3B y 4B respectivamente (ver apéndice B). Finalmente se aplicó la
fórmula definida para este método y se obtiene el resultado de Z=0.8549, que
resalta en la celda de color amarillo. Los resultados de Z con los valores de
67

Pitzer

Pitzer 0.8549 0.8656 0.9635 0.9736
Fuente: Los autores
Nótese que al igual que en los casos anteriores las diferencias

significativas se apreciaron entre las muestras provenientes de los reportes
del cromatografo ACOGAS, a saber muestras 1 y 2 (ver apéndice A), con
respecto a las muestras 3 y 4, provenientes de un análisis PVT del campo
Santa Bárbara (ver apéndice A). Esto se debió a la cantidad de
consideraciones que toma este método para el cálculo de Z, como se aprecia
en la figura 9C (ver apéndice C). Al igual que en los casos anteriores se
mantuvo la tendencia ascendente de los valores de Z.
3.3. COMPARACIÓN A TRAVÉS DE UN ANÁLISIS ESTADÍSTICO DE

OBTENIDOS
Luego de la determinación de 36 valores del factor de compresibilidad

Z, se realizó la comparación de los valores obtenidos con el uso de la prueba
de Friedman, siendo esta una prueba estadística no paramétrica, la cual se
desarrolló bajo el entorno del programa estadístico Statgraphic, considerando
que las muestras utilizadas para el cálculo de Z fueron escogidas al azar y
que los métodos aplicados no tienen relación uno con el otro en cuanto a los
68
criterios y metodología de uso.
3.3.1 Resumen estadístico
En la tabla 3.11, se muestran los estadísticos descriptivos obtenidos

mediante el programa Statgraphic, en base a los resultados del factor de
compresibilidad Z (ver tabla 3.12) calculados por nueve (9) métodos
diferentes.
Tabla 3.11- Análisis descriptivo de los valores de Z calculados
Métodos Recuento Promedio Desviación Coeficiente Mínimo Máximo Rango

Estándar de Variación
VDW 4 0,856225 0,0739518 8,63696% 0,7835 0,9207 0,1372
PR 4 0,8064 0,0903396 11,2028% 0,7183 0,8853 0,167
RK 4 0,872125 0,0627468 7,1947% 0,8093 0,927 0,1177
SRK 4 0,915175 0,0545297 5,95839% 0,8638 0,9631 0,0993
DAK 4 0,851075 0,00697346 0,819371% 0,8433 0,8575 0,0142
DPR 4 0,91785 0,042412 4,6208% 0,8778 0,9549 0,0771
S-K 4 0,905 0,0591608 6,5371% 0,84 0,96 0,12
EMR 4 0,8825 0,0613052 6,94677% 0,82 0,94 0,12
PITZER 4 0,9144 0,0628149 6,86952% 0,8549 0,9736 0,1187
Total 36 0,880083 0,0646392 7,34467% 0,7183 0,9736 0,2553
Fuente: Statgraphic
Según los resultados arrojados por el programa, el método que

presentó mayor homogeneidad fue el basado en la ecuación de Dranchuk
Abou Kassen, ya que su coeficiente de variación fue CV=0.819371%.
Mientras que el método con menor homogeneidad fue la ecuación de estado
de Peng Robinson con un CV=11.2028%. El coeficiente de variación se
refiere básicamente a que tanto cambia el valor de Z para un mismo método,
con respecto a la muestra que se estudió.
69
Tabla 3.12- Resultados del factor de compresibilidad Z para cada método

definido
Método VDW PR RK SRK DAK DPR S-K EMR PITZER
Muestra1 0.7835 0.7183 0.8093 0.8638 0.8433 0.8778 0.84 0.82 0.8549
Muestra2 0.8015 0.7387 0.8270 0.8723 0.8471 0.8846 0.87 0.84 0.8656
Muestra3 0.9192 0.8833 0.9252 0.9615 0.8564 0.9541 0.95 0.93 0.9635
Muestra4 0.9207 0.8853 0.9270 0.9631 0.8575 0.9549 0.96 0.94 0.9736
Fuente: Los autores
3.3.2 Cálculo para el estadístico de prueba
Considerando los valores de factor de compresibilidad Z calculado

para cada método definido como se muestra en la tabla 3.12, se aplicó el
procedimiento para la prueba de Friedman, asignándole un rango a cada
valor y obteniendo como resultado lo mostrado en la siguiente tabla:
Tabla 3.13- Rangos asignados a cada método y resultados de Rj
Método VDW PR RK SRK DAK DPR S-K EMR Pitzer
Muestra 1 2 1 3 8 6 9 5 4 7
Muestra 2 2 1 3 8 5 9 7 4 6
Muestra 3 3 2 4 8 1 7 6 5 9
Muestra 4 3 2 4 8 1 6 7 5 9
Rj 10 6 14 32 13 31 25 18 31
2
Rj 100 36 196 1024 169 961 625 324 961
Fuente: Los autores
70
La tabla 3.13 representa el procedimiento para la aplicación de la

prueba de estadística de Friedman. De ella se obtuvo el valor de chi-
cuadrado aplicando la ecuación 2.41, que representa el estadístico de
prueba. Para ello se consideró n=4 y k=9. Del mismo modo se plantearon las
siguientes consideraciones:
 Hipótesis nula (Hl o): no hay diferencias estadísticamente significativas
entre los nueve (9) métodos estudiados.
 Hipótesis alternativa (Hl): si existen diferencias estadísticamente
significativas entre los nueve (9) métodos estudiados.
 Nivel de significancia= 5%
 Grados de libertad (gl)= k-1 = 9-1= 8
Entonces el valor de chi-cuadrado (ecuación 2.41) viene dado por:
12
𝑋2𝑟 = (∑ 4396) − 3 ∗ 4(9 + 1) = 26.533
4 ∗ 9 (9 + 1)
Este valor calculado por el programa se comparó con el valor tabulado

de la distribución chi-cuadrado que se muestra en la figura 1D (ver apéndice
D), entrando con los grados de libertad (8) y cortando con la columna de
nivel de significancia (0.05) y se obtuvo el valor Xr2 (8, 0.05) = 15.51. Como el
valor tabulado es menor que el valor calculado (26.533), se rechazó la
hipótesis nula y se aceptó la hipótesis alternativa, en consecuencia se
concluyó que, estadísticamente existen diferencias significativas entre los
métodos aplicados con un 95% de confianza.
Según el análisis estadístico realizado el método más relevante, fue el

basado en la ecuación de estado de Soave-Redlich-Kwong, por presentar el
mayor rango (Rj=32), mientras que el menos relevante fue el método basado
71
en la ecuación de estado de Peng Robinson (Rj=6). Finalmente se pudo

establecer, que considerando los métodos definidos para el cálculo del factor
de compresibilidad Z y basándose en los resultados obtenidos con la prueba
estadística de Friedman, el método que resultó más apropiado para el
cálculo de Z, tomando en cuenta las cuatro muestras de gas seleccionadas
con sus respectivas condiciones de operación (presión y temperatura), fue el
método basado en la ecuación de estado de Soave-Redlich-Kwong,
seguidamente el ajuste Dranchuk-Purvis-Robinson y el método de Pitzer
respectivamente.
CAPÍTULO IV
CONCLUSIONES Y RECOMENDACIONES
4.1. Conclusiones
 Se describieron 13 de los métodos más utilizados para el cálculo del factor

de compresibilidad Z.
 Se determinó el factor de compresibilidad Z a 4 muestras de gas natural en
base a 9 métodos representativos, de los cuales 6 se basaron en ecuaciones
de estado y 3 en métodos gráficos.
 Se realizaron hojas de cálculo en Microsoft Excel que permitieron en base a
las condiciones de operación y a la composición conocida de cada muestra
en particular, determinar el factor de compresibilidad del gas para cada uno
de los 9 métodos definidos.
 Se obtuvieron 36 valores del cálculo de Z, los cuales fueron comparados con
el uso de la prueba de Friedman bajo el entorno del programa estadístico
Statgraphic.
 El uso del programa estadístico Statgraphic permitió obtener un análisis
descriptivo de los valores de Z calculados en el cual se pudo apreciar lo
siguiente: el método que presentó mayor homogeneidad fue el basado en la
ecuación de estado de Dranchuk Abou Kassen y el de menor homogeneidad
fue el método de Peng-Robinson, considerando para ello sus respectivos
coeficientes de variación.
 Mediante la prueba estadística de Friedman se demostró que existen
diferencias significativas entre los métodos estudiados con un 95% de
confianza.
 La prueba estadística de Friedman indicó que de acuerdo a los resultados
obtenidos, el método más apropiado para el cálculo del factor de
72
73
compresibilidad Z, fue el de Soave-Redlich-Kwong y el menos apropiado el

de Peng-Robinson. Estas consideraciones se basan principalmente en las 4
muestras de gas natural, con sus respectivas condiciones de operación
(presión y temperatura) que inicialmente fueron seleccionadas.
4.2. Recomendaciones
 Se recomienda el uso computacional bajo el entorno de Microsoft

Excel y programas estadísticos, ya que los mismos permiten la simplificación
de los cálculos de ingeniería.
 Ampliar la investigación utilizando una mayor cantidad de muestras y
métodos respectivamente, para evaluar el comportamiento del factor de
compresibilidad Z.
 Se recomienda validar los resultados obtenidos por medio de métodos
gráficos como Standing-Katz, Eykman y Pitzer, ya que al ser valores
obtenidos de gráficas, es posible una mala lectura o apreciación y por ende
dicho resultado tenga un porcentaje de error significativo.
 Para las muestras de gas natural utilizadas, calcular el factor de
compresibilidad Z por el método de Soave-Redlich-Kwong, seguidamente por
el ajuste de Dranchuk-Abou-Kassen y el método de Pitzer, los cuales
permitirán obtener mejores resultados.
 Evitar el uso de la ecuación de estado desarrollada por Peng
Robinson para el cálculo del factor de compresibilidad Z utilizando las
muestras de gas seleccionadas, ya que no arroja buenos resultados.
REFERENCIAS BIBLIOGRÁFICAS
Ángel, J. (1998). Estadística general aplicada. Universidad Eafit.
Banzer, C. (1996). Correlaciones numéricas PVT. Edición preliminar.

Maracaibo, Venezuela.
Barrios, R. y Martínez, K. (2011). Validación del factor de compresibilidad Z,

a través del gasoducto ACOGAS de la planta de inyección de gas El Furrial
mediante ecuaciones de estado. Trabajo de grado, Universidad de Oriente,
Núcleo de Monagas. Maturín, Venezuela.
Cengel, A. y Boles, M. (2002). Thermodynamics: an engineering approach.

Boston: McGraw-Hill. pp. 91-93.
Cotton, F. (1996). Química general superior. Interamericana, México.
Figueira, F. (2005). Desarrollo de ecuaciones de estado del tipo Van Der

Waals para fluidos polares y no polares. Trabajo de ascenso. Universidad
Simón Bolívar.
González, A. (2006). Análisis regresional de modelos matemáticos

polinómicos para relacionar propiedades termodinámicas del gas natural
venezolano. Trabajo de grado, Universidad Rafael Urdaneta. Maracaibo,
Venezuela.
Engineering Data Book (1982), Natural Gas Processors Supliers

Asociation, GPSA, 1982 First Place, 5th & Boston, Tulsa, Oklahoma 74103.
Tenth Edition, 1987.
Irragorry, M. (2004). Elaboración de un programa de simulación para el

cálculo de equilibrio líquido-vapor en columnas de separación usadas en el
74
75
procesamiento del gas natural. Trabajo de grado, Universidad Rafael

Urdaneta. Maracaibo, Venezuela.
Malanowski, S. y Anderko, A. (1992). Modeling phase equilibrium.
Thermodynamic background and practical tools. Wiley. New Jersey.
Martínez, A. y Ramos, J. (2010). Estudio del impacto ambiental generado

por el uso del gas metano y la gasolina como combustible vehicular. Trabajo
de grado, Universidad de Oriente, Núcleo Monagas. Maturín, Venezuela.
Pérez, R. y Martínez, M. (1994). Ingeniería de gas natural: características y

comportamiento de los hidrocarburos. Ingenieros consultores. Maracaibo,
Venezuela.
Petrucci, R. (1988). Química general. Sistemas técnicos de edición.

Addison- Wesley Iberoamericana. México.
Piersol, Allan. y Bendat Julius. (2010). Analysis and Measurement

Procedures, Fourth edition. Wiley Editorial
Pino, F. (2008a). Características y propiedades del gas natural y sus

derivados (Documento en formato PDF). Universidad de Oriente, Núcleo de
Monagas. Maturín, Venezuela.
Pino, F. (2008b). Mediciones y medidores del gas natural (Documento en

formato PDF). Universidad de Oriente, Núcleo de Monagas. Maturín,
Venezuela.
Pino, F. (2008c). Proceso de compresión y expansión del gas natural.

Escuela de Ingeniería de Petróleo. Universidad de Oriente, Núcleo de
Monagas. Maturín, Venezuela.
Rodríguez, R. (2014). Estimación de propiedades termodinámicas de

mezclas de gases en condiciones reales utilizando una herramienta
76
computacional. Trabajo de ascenso, Universidad de Carabobo. Carabobo,

Venezuela.
Rojas, G. (2003). Ingeniería de yacimientos de gas condensado. Universidad

de Oriente. Puerto la Cruz, Venezuela.
Romero, H. (2003). Procedimiento para la determinación experimental del

factor de desviación de un gas natural. Trabajo de ascenso, Universidad del
Zulia. Maracaibo, Venezuela.
Sabino, M. (2010). Química general aplicada a las ciencias de la salud. 1ª

edición. Mérida, Venezuela.
Skoog, D. y Leary, J. (2008). Principios del análisis instrumental. Editorial

Saunders College Pub. Universidad de Michigan.
Triana, M. (2012). Propuesta experimental aplicada al aula para la

enseñanza del tema de gases. Trabajo de ascenso, Universidad nacional de
Colombia. Bogotá. DC, Colombia.
APÉNDICES
77
APÉNDICE A
“Muestras de gas natural utilizadas para el cálculo del factor
de compresibilidad Z”
Tabla 1A- Cromatografías utilizadas para la determinación del factor de
compresibilidad Z por nueve métodos distintos
Muestra 1A Muestra 2A
P=1232.59 psi P=1233.57 psi
T=107.07 °F T= 104.74 °F
Componente Fracción molar Fracción molar
Yi Yi
(%) (%)
Nitrógeno 0.20 0.0020 0.24 0.0024
Dióxido de
3.40 0.0340 4.32 0.0432
Carbono
Metano 86.08 0.8608 84.94 0.8494
Etano 7.48 0.0748 9.19 0.0919
Propano 1.87 0.0187 0.82 0.0082
I-Butano 0.22 0.0022 0.12 0.0012
N-Butano 0.38 0.0038 0.22 0.0022
I-Pentano 0.15 0.0015 0.09 0.0009
N-Pentano 0.12 0.0012 0 0
Hexano 0.10 0.0010 0.06 0.0006
Total 100 1 100 1
Fuente: Reporte de inyección Cromatografo ACOGAS

Tabla 2A- Cromatografías utilizadas para la determinación del factor de
compresibilidad Z por nueve métodos distintos
Muestra 3A Muestra 4A
P= 681.8 lpca P=682.6 lpca
Componente T=693.5 °R T= 694 °R
Fracción molar Yi Fracción molar Yi
(%) (%)
Nitrógeno 0.46 0.0046 0.22 0.0022
Dióxido de
4.81 0.0481 4.91 0.0491
Carbono
Metano 80.52 0.8052 81.12 0.8112

Etano 7.95 0.0795 7.62 0.0762
Propano 3.34 0.0334 3.31 0.0331
I-Butano 0.72 0.0072 0.69 0.0069
N-Butano 0.98 0.0098 0.92 0.0092
I-Pentano 0.36 0.0036 0.35 0.0035
N-Pentano 0.29 0.0029 0.29 0.0029
Hexano 0.27 0.0027 0.27 0.0027
Heptano+ 0.19 0.0019 0.19 0.0019
Octano 0.09 0.0009 0.09 0.0009
Nonano 0.02 0.0002 0.02 0.0002
Total 100 1 100 1
Fuente: Análisis PVT del campo Santa Bárbara, edo Anzoátegui.

APÉNDICE B
“Gráficas utilizadas para el cálculo del factor de
compresibilidad Z”
Figura 1B- Factor de compresibilidad de gases naturales como función
de Psr y T sr por el método de Standing-Katz.
Figura 2B- Factor de corrección Z, como función de Pr y Tr por el
método de Refracción Molecular de Eykman.
Figura 3B- Factor de compresibilidad Z 0 como función de Pr y T r por el
método de Pitzer.
Figura 4B- Factor de compresibilidad Z’ como función de Pr y Tr por el
método de Pitzer.
APÉNDICE C
“Hojas de cálculo en Microsoft Excel utilizadas para
desarrollar cada uno de los métodos”
Figura 1C- Programa en hoja de cálculo para determinar el factor de
compresibilidad Z por la ecuación de estado Van Der Waals. (Fuente:
los autores)
compresibilidad Z por la ecuación de estado Peng-Robinson. (Fuente:
los autores)
compresibilidad Z por la ecuación de estado Redlich-Kwong. (Fuente:
los autores)
compresibilidad Z por la ecuación de estado Soave-Redlich-Kwong.
(Fuente: los autores)
compresibilidad Z por la ecuación de estado Dranchuk-Abou-Kassen.
(Fuente: los autores)
compresibilidad Z por la ecuación de estado Dranchuk-Purvis-
Robinson. (Fuente: los autores)
compresibilidad Z con el método gráfico de Standing-Katz. (Fuente: los
autores)
compresibilidad Z con el método gráfico de Eykman. (Fuente: los
autores)
compresibilidad Z con el método gráfico de Pitzer. (Fuente: los autores)
APÉNDICE D
“Datos tabulados para desarrollar la prueba estadística de
Friedman”
Figura 1D- Distribución probabilística chi-cuadrado para n y k grandes
Fuente: Random Data: Analysis and Measurement Procedures, Fourth

edition. Allan G. Piersol.
HOJAS METADATOS
Hoja de Metadatos para Tesis y Trabajos de Ascenso - 1/6
ANÁLISIS COMPARATIVO DE LAS DIFERENTES
METODOLOGÍAS PARA LA DETERMINACIÓN DEL
Título
FACTOR DE COMPRESIBILIDAD (Z) DEL GAS
APLICADO A MUESTRAS DE GAS NATURAL
El Título es requerido. El subtítulo o título alternativo es opcional.
Autor(es)
Apellidos y Nombres Código CVLAC / e-mail
Carrera Marcano Yérika CVLAC C.I: 21.051.876
Glecmar
e-mail
Castillo La Rosa Javier CVLAC C.I: 19.403.248
Eduardo
e-mail
CVLAC C.I: 19.804.888
Roque Pérez Félix David
e-mail
Se requiere por lo menos los apellidos y nombres de un autor. El formato
para escribir los apellidos y nombres es: “Apellido1 InicialApellido2., Nombre1
InicialNombre2”. Si el autor esta registrado en el sistema CVLAC, se anota el
código respectivo (para ciudadanos venezolanos dicho código coincide con
el numero de la Cedula de Identidad). El campo e-mail es completamente
opcional y depende de la voluntad de los autores.
Palabras o frases claves:
Factor de compresibilidad
Ecuaciones de estado
Gas
Prueba de Friedman
Análisis estadístico
Líneas y sublíneas de investigación:
Área Sub-área
Ingeniería de Petróleo
Tecnología y Ciencias Aplicadas
Debe indicarse por lo menos una línea o área de investigación y por cada
área por lo menos un subárea. El representante de la subcomisión solicitará
esta información a los miembros del jurado.
Resumen (Abstract):
El factor de compresibilidad Z, es uno de los parámetros que, con
mayor precisión diferencia el comportamiento de los fluidos en estado
líquido del estado gaseoso. Además de esto, define el comportamiento
de los gases a determinadas condiciones de presión y temperatura. Lo
que hace que sea de vital importancia, para todos los diseños e
instalaciones que trabajan con fluidos compresibles. El objetivo
principal de esta monografía de investigación fue realizar un análisis
comparativo entre nueve métodos seleccionados para el cálculo del
factor de compresibilidad Z del gas natural, mediante la aplicación de la
prueba estadística de Friedman, con la finalidad de establecer cuál de
ellos es el más apropiado para un resultado óptimo y eficaz tomando en
cuenta cuatro muestras de gas obtenidas de tesis de investigación con
sus respectivas condiciones de operación. El análisis estadístico indicó
que de acuerdo a los resultados alcanzados, el método más apropiado
para el cálculo del factor de compresibilidad Z, fue el de Soave-Redlich-
Kwong y el menos apropiado el de Peng-Robinson, tomando en cuenta
las cromatografías mencionadas anteriormente.
Contribuidores:
Apellidos y Nombres Código CVLAC / e-mail
ROL CA AS TU JU
Msc. Tomás Marín CVLAC C.I: 11.538.773
e-mail [email protected]
ROL CA AS TU JU
Ing. Miguel Flores CVLAC C.I: 10.467.516
ROL CA AS TU JU
Msc. Pedro Tineo CVLAC C.I: 10.100.632
ROL CA AS TU JU
Msc. Rubén Vega CVLAC C.I: 12.155.280
Se requiere por lo menos los apellidos y nombres del tutor y los otros dos (2)
jurados. El formato para escribir los apellidos y nombres es: “Apellido1
InicialApellido2., Nombre1 InicialNombre2”. Si el autor esta registrado en el
sistema CVLAC, se anota el código respectivo (para ciudadanos
venezolanos dicho código coincide con el numero de la Cedula de Identidad).
La codificación del Rol es: CA = Coautor, AS = Asesor, TU = Tutor, JU =
Jurado.
Fecha de discusión y aprobación:

Año Mes Día
2018 05 10
Fecha en formato ISO (AAAA-MM-DD). Ej.: 2005-03-18. El dato fecha es requerido.
Lenguaje: spa Requerido. Lenguaje del texto discutido y aprobado, codificado usando
ISO 639-2. El código para español o castellano es spa. El código para
ingles en. Si el lenguaje se especifica, se asume que es el inglés (en).
Archivo(s):
Nombre de archivo
TESIS_DIGITAL_Z_FELIX_YERIKA_JAVIER.docx
Caracteres permitidos en los nombres de los archivos: A B C D E F G H I J K
LMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwx
yz0123456789_-.
Alcance:
Espacial: __________________ (opcional)
Temporal: __________________ (opcional)
Título o Grado asociado con el trabajo:

Ingeniero de Petróleo
Dato requerido. Ejemplo: Licenciado en Matemáticas, Magister
Scientiarium en Biología Pesquera, Profesor Asociado, Administrativo III,
etc.
Nivel Asociado con el trabajo: Ingeniería

Dato requerido. Ej.: Licenciatura, Magister, Doctorado, Post-doctorado,
etc.
Área de Estudio:
Tecnología y Ciencias Aplicadas
Usualmente es el nombre del programa o departamento.
Institución(es) que garantiza(n) el Título o grado:

Universidad de Oriente Núcleo Monagas
Si como producto de convenciones, otras instituciones además de la
Universidad de Oriente, avalan el título o grado obtenido, el nombre de estas
instituciones debe incluirse aquí.
Hoja de metadatos para tesis y trabajos de Ascenso- 5/6
Hoja de metadatos para tesis y trabajos de Ascenso- 6/6
De acuerdo al Artículo 41 del reglamento de Trabajos de Grado:
Los Trabajos de Grado son de la exclusiva propiedad de la

Universidad de Oriente, y sólo podrán ser utilizados a otros fines con el
consentimiento del Consejo de Núcleo respectivo, quién deberá
participarlo previamente al Consejo Universitario, para su autorización.
Javier Castillo Yérika Carrera

C.I. 19.403.248 C.I. 21.051.876
Félix Roque
C.I. 19.804.888
Autores
Msc. Tomás Marín
Asesor

Monografía Factor de Compresibilidad Z PDF

Cargado por

Copyright:

Formatos disponibles

Monografía Factor de Compresibilidad Z PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Monografía Factor de Compresibilidad Z PDF

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE ORIENTE

ANÁLISIS COMPARATIVO DE LAS DIFERENTES METODOLOGÍAS

Seminario de Área de Grado presentado como requisito parcial para

MATURÍN, MAYO DE 2018

ANÁLISIS COMPARATIVO DE LAS DIFERENTES METODOLOGÍAS

MATURÍN, MAYO DE 2018

ANÁLISIS COMPARATIVO DE LAS DIFERENTES METODOLOGÍAS

MATURÍN, MAYO DE 2018

DE ACUERDO AL ARTÍCULO 41 DEL REGLAMENTO DE TRABAJOS DE

“LOS TRABAJOS DE GRADO SON PROPIEDAD EXCLUSIVA DE LA

A Dios primeramente por permitirme llegar hasta este punto de vida y

A mi Madre, que desde el primer día me acompañó en este transitar

A mi hermana, mi concu, mi compañera de aventuras Virginia Roque,

A mis compañeros de clase, que hicieron más llevadero y ameno

A DIOS todopoderoso, por haberme dado salud y sabiduría para

Agradezco también la crianza, consejos y apoyo brindado por parte de

A mis hermanos, en especial a mi hermana Daicelia Zapata por

A mis compañeros de estudio, por todos los momentos vividos dentro

Gracias a todas aquellas personas que en algún momento de mi vida

A mis hermanos Junior y Yunaika Carrera, en especial a ti mi negra, por

A mis compañeros de este trabajo, especialmente a Félix Roque por

A mis amigos, en especial a María Pereda, Yulimar Reyes y Angélica

A mis Padres, especialmente a mi Madre Bethsy Pérez a quien

A mi Abuela, Berta Roque, esa que me abrió las puertas de su casa

A mis hermanos, Carlos y Virginia, especialmente a Carlos espero

A la princesa Miranda Victoria, la sobrina que llegó a esta vida para

Dedico este trabajo primeramente a Dios por permitirme llegar a este

A mis padres, que han sabido formarme con buenos principios y

A mis hermanos, por compartir alegrías y tropiezos de los cuales

A mi familia en general, por brindarme su apoyo incondicional y por

A mi grupo de amigos, que gracias al equipo que hicimos logramos

A mis profesores, por su tiempo, apoyo y así como por los

A Dios Todopoderoso por todas las bendiciones que ha derramado

A mis padres y hermanos por su infinito amor y por su apoyo

A ti José Reyes que te has convertido en pieza fundamental en mi vida

A mi tía Jenny Marcano, por brindarme parte de su tiempo, por sus

A mis abuelos por ser parte importante de mi vida y por haber

A ustedes y todas esas personas q creyeron en mí.

AUTORES ASESOR ACADÉMICO

Yérika Glecmar Carrera Marcano Msc. Tomás Marín

El factor de compresibilidad Z, es uno de los parámetros que, con mayor

Palabras claves: Factor de compresibilidad, ecuaciones de estado, prueba

El gas natural es una mezcla homogénea de hidrocarburos parafínicos

Ahora bien el factor de compresibilidad Z, es un factor de corrección,

La presente monografía de investigación tiene por objeto realizar un

1.1. PLANTEAMIENTO DEL PROBLEMA

El comportamiento de presión, volumen, temperatura y moles para un

Los gases naturales de Venezuela presentan características típicas que

1.2. OBJETIVOS DE LA INVESTIGACIÓN

1.2.1 Objetivo general

Analizar comparativamente las diferentes metodologías para la

1.2.2 Objetivos específicos

 Describir los métodos más utilizados para la determinación del factor

 Determinar el factor de compresibilidad Z a cuatro (4) muestras de gas

natural, utilizando nueve (9) métodos representativos para su cálculo.

 Comparar a través de un análisis estadístico las diferencias existentes

entre los resultados obtenidos.

1.3. JUSTIFICACIÓN DE LA INVESTIGACIÓN

El factor de compresibilidad Z, es uno de los parámetros que, con