Diagramas y tablas termodinámicas de crudos pesados de la
faja del Orinoco
Diagrams and thermodynamic tables of heavy oils from
Orinoco belt
G. Vielma., R. Varillas y J.I. Belandria
Escuela de Ingeniería Química, Facultad de Ingeniería, Universidad de Los Andes,
Mérida- Venezuela.
[email protected]
Resumen
Los diagramas y tablas termodinámicas de varios crudos pesados de la faja del Orinoco como Cerro Negro,
Hamaca, Jobo, Merey y Zuata son presentados en este trabajo. La información termodinámica fue obtenida
empleando el simulador de procesos Aspen Hysys con modelos termodinámicos adecuados para examinar
los datos experimentales de estos crudos como Redlich Kwong Soave, Peng Robinson y Peng Robinson
Stryek Vera. En particular, se muestran análisis y comparaciones de los diagramas presión-volumen,
temperatura-volumen, presión-temperatura, entropía-presión, temperatura–entalpía, presión-entalpía,
temperatura-entropía y presión-entropía. Asimismo, las tablas termodinámicas presentan información sobre
volumen específico, entalpía y entropía a diferentes presiones y temperaturas en las regiones de líquido
comprimido, líquido saturado, vapor saturado y vapor sobrecalentado de los crudos mencionados. Esta
información es útil para ayudar a comprender el comportamiento de estas sustancias y diseñar las unidades
de fraccionamiento y conversión para transformarlos en productos de importancia energética y
petroquímica.
Palabras clave: Diagramas y tablas termodinámicas, crudos pesados, faja Orinoco.
Abstract
The thermodynamic diagrams and tables of some typical heavy oils from Orinoco belt like Cerro Negro,
Hamaca, Jobo, Merey, and Zuata are presented in this work. The thermodynamic information was got by
using the process simulation software Aspen Hysys with suitables thermodynamic models to examine the
experimental data of these oils like Redlich Kwong Soave, Peng Robinson, and Peng Robinson Stryek Vera. In
particular, it is shown the analysis and comparison of the following thermodynamic diagrams: Pressure-
volume, pressure-temperature, temperature-volume, entropy-pressure, temperature-enthalpy, pressure-
enthalpy, temperature-entropy and pressure-entropy. Also, the thermodynamic tables present information
about volume, entropy and enthalpy at different pressures and temperatures in the regions of compressed
liquid, superheated vapor, saturated liquid and saturated vapor of these crude oils. This information is
helpful to understand the behavior of these substances and to design the fractionating and conversion units to
convert them in products of energetic and petrochemical importance.
Keywords: Diagrams and thermodynamic tables, heavy oils, Orinoco belt.
1 Introducción fundamentada, principalmente, en los crudos
livianos y medianos explotados antes de 1950
Venezuela posee las reservas de crudos pesados (Maxwell, 1950). Por lo tanto, sus
y extra pesados más grandes del mundo en la faja extrapolaciones al entorno de los crudos pesados
del Orinoco por lo cual es pertinente desarrollar son restringidas y aproximadas.
tecnologías para su conversión en productos
útiles para el sector energético e industrial. Ante esta circunstancia, queda la alternativa de
utilizar, los simuladores de procesos modernos
Para lograr este objetivo es necesario crear utilizando ecuaciones de estado convenientes para
tecnologías capaces de transformar estos operar con crudos pesados como los presentes en
materiales, los cuales debido a su naturaleza física la faja del Orinoco. Conviene mencionar que
y química son difíciles de extraer, transportar y algunas ecuaciones de estado como Peng-
procesar. En general, la mayoría de las tecnologías Robinson (PR), Peng-Robinson-Stryek–Vera
propuestas todavía no están completamente (PRSV), Soave-Redlich-Kwong (SRK), entre
desarrolladas y muchas de ellas se encuentran otras, han sido usadas para modelar hidrocarburos
bajo estudios a diferentes escalas en INTEVEP, de distintas características (Firoozabadi, 1999;
en las universidades nacionales y en otros países PRO II KEY WORD MANUAL, 1991; ASPEN
del mundo (Belandria, 1987; Vielma, 2015). PLUS Manual, 2004; FAQ Aspen HYSYS,
2010).
Ante esta situación y reconociendo que la
información acerca de las propiedades Al respecto, los creadores del simulador Aspen
termodinámicas son útiles y fundamentales en el Hysys sugieren emplear para la simulación de
diseño de cualquier proceso físico y químico, es crudos pesados la ecuación de estado de PR ó
evidente que la estimación de las propiedades PRSV (Aspen HYSYS, 2010). En su opinión,
termodinámicas de los crudos pesados y extra dejan entrever, que la ecuación PRSV es una
pesados constituye un apoyo para comprender mejor alternativa, porque los crudos pesados
mejor su naturaleza y diseñar racionalmente los presentan metales capaces de crear interacciones
esquemas de procesamiento, instalaciones y intermoleculares, relativamente importantes. En
equipos adecuados para transformar estos este caso, la ventaja de PRSV sobre PR, es la
materiales en sustancias de interés energético e presencia de un parámetro adicional en la
industrial, entre otras posibilidades. ecuación PRSV, el cual toma en cuenta posibles
interacciones entre las moléculas metálicas y
En este sentido, en el presente trabajo se muestran
moléculas vecinas. La ecuación PR tiene un solo
diagramas y tablas termodinámicas de los crudos
parámetro que cuantifica las interacciones
pesados Cerro Negro, Hamaca, Jobo, Merey y
moleculares comunes en los hidrocarburos no
Zuata. Estos diagramas y tablas fueron
polares, mientras que la ecuación PRSV posee
construidos con el simulador Hysys utilizando
dos parámetros para ponderar el impacto de las
modelos termodinámicos adecuados para los
fuerzas intermoleculares polares existentes en los
crudos pesados de la faja del Orinoco (Vielma,
asfaltenos del crudo (Peng y col., 1976; Stryek y
2015). Es pertinente señalar, que hasta el nivel del
col., 1986).
conocimiento de los autores, no se han realizado
estudios previos de esta naturaleza, lo cual Basado en una investigación previa se compararon
justifica la intención de este artículo.
las predicciones de los modelos PRSV, PR y
SRK, y se concluyó que el modelo PRSV tiende
a ajustar adecuadamente los datos experimentales
2 Estimación de propiedades termodinámicas de los crudos pesados de la faja del Orinoco
con el simulador Hysys (Vielma 2015). Por esta razón se escogió este
modelo para estimar las propiedades
Debido a la complejidad de la estructura termodinámicas de los crudos pesados y construir
molecular de los crudos pesados, los métodos los correspondientes diagramas y tablas
clásicos para estimar las propiedades
termodinámicas de los hidrocarburos como los
basados en el famoso libro Data Book on
Hydrocarbons son limitados, porque la
información presentada en este libro está
3 Simulación de los crudos pesados. 4 Resultados y discusión
Para lograr estimar las propiedades 4.1 Obtención de la curvas de destilación TBP de
termodinámicas de un crudo cualquiera como un los crudos pesados
crudo pesado, empleando una ecuación de estado
especifica como PRSV, PR, SRK, entre otras, los
simuladores de procesos como el Aspen HYSYS Para construir las curvas de destilación TBP de los
utilizan las curvas TBP del crudo para generar un crudos pesados de la faja del Orinoco estudiados
conjunto de pseudocomponentes que simulan el en el presente proyecto, se introducen al
crudo de interés (Chang y col., 2012). simulador los datos experimentales de
destilación disponibles en trabajos publicados en
También, es necesario suministrar al simulador la INTEVEP y el Laboratorio de Petróleo y Catálisis
gravedad API de las fracciones del crudo o en su de la ULA (Montoya, 1985; Bellet, 1998).
defecto la gravedad API del crudo global. Si se Luego, el simulador genera las correspondiente
proporciona la gravedad API del crudo global, curvas TBP siguiendo el procedimiento descrito
como en el presente caso, el simulador estima el en (Vielma, 2015). Las curvas TBP obtenidas se
factor de caracterización K del crudo global presentan en la Fig. 1
conocido el punto de ebullición medio del crudo
y su densidad SG (Stryek y col., 1986).
En seguida, el simulador asume que todos los 1000 Jobo
pseudo componentes tienen un factor de Hamaca 2X
caracterización constante e igual al del crudo Cerro Negro
800
Temperatura (°C)
global. Conocido el factor de caracterización y Zuata
el punto de ebullición normal de cada pseudo Merey
componente, el simulador calcula su densidad 600
usando un procedimiento rutinario.
Con esta información básica y un algoritmo 400
iterativo, el simulador define y caracteriza todos
los pseudo componentes del crudo asignándole a 200
cada a uno ellos la densidad y el punto de
ebullición normal correspondiente.
0
Si se dispone de la viscosidad, el simulador asigna 0 % volumen 100
también una viscosidad a cada pseudo
componente; si no se posee esta información, el
simulador no toma en cuenta esta caracterización
adicional, y usa solamente la densidad y el punto Fig.1. Curvas TBP de los crudos
de ebullición normal como en la situación actual.
Una vez definidos los pseudo componentes, según
En la Tabla 1 se muestran el peso molecular
el esquema anterior, el simulador puede generar
MW , la presión crítica Pc, la temperatura
todas las propiedades termo físicas requeridas
crítica Tc, la presión cricondentermica Pcrin y la
para la simulación de cualquier proceso, como el
temperatura cricondentermica Tcrin de los crudos
peso molecular MW, temperatura crítica Tc,
Merey , Zuata, Jobo, Cerro Negro y Hamaca .
presión crítica Pc, factor acéntrico ω , densidad
También se presenta la gravedad API
SG , viscosidad σ , conductividad térmica π ,
experimental y estimada
entalpía H , entropía S, y otras que puedan ser
requeridas por el usuario. Todas estas propiedades
se estiman tanto para los pseudo componentes
como para el crudo global. Las del crudo global se
evalúan combinando las de los pseudo
componentes usando reglas de mezcla adecuadas.
Tabla 1. Propiedades críticas y cricondentérmicas de los crudos pesados
Merey Zuata Jobo Cerro negro Hamaca
°API 16 8,1 12 8 9,4
°API (modelos) 16,316 8,406 11,927 8,286 9,365
MW ( kg / k mol) 303,8 371,7 476,6 514,1 601,5
Tc (°C) 661,3 727,6 783,1 814,0 837
Pc (kPa) 3537,9 2661 2326 2054,0 1601
Vc (m³ / K mol) 1,1 1,2 1,9 2,0 2,3
T Crin (°C) 720,0 757,4 854,9 857,7 876,5
P Crin (k Pa) 3560,7 2668 2330 2072,0 1610
En la tabla se detecta como los pesos de los De manera similar, es posible explicar el
crudos estudiados oscilan entre 303.8 y 601.5 comportamiento presentado por la temperatura y
Kg/Kmol sugiriendo una alta complejidad presión cricondentermicas.
molecular. También, es posible detectar como el
volumen crítico, la temperatura crítica y la En general, no se observa una tendencia definida
temperatura cricondentérmica crecen al aumentar entre las propiedades mencionadas y la gravedad
el peso molecular. En cambio, la presión crítica y API de los crudos.
cricondentérmica decrecen con el aumento del
peso molecular.
El comportamiento crítico observado es 4.3 Propiedades termodinámicas de los crudos
explicable, porque en general, la temperatura pesados a condiciones estándar
crítica y el volumen crítico tienden a crecer al
incrementarse el peso molecular, como se En la Tabla 2 se presentan las propiedades
aprecia en las comparaciones anteriores (Maxwell, termodinámicas de los crudos pesados a las
1950; Smith y Van Ness, 2001). Esto ocurre condiciones estándar de 101.3 kPa y 15.56 °C. En
porque la temperatura crítica está asociada con particular se presentan los valores del volumen V,
parámetros energéticos relacionados con la entalpía H, entropía S, Exergía Ex, factor de
intensidad de las fuerzas intermoleculares, las comprensibilidad Z, capacidad calorífica Cp y Cv.
cuales tienden a intensificarse con la masa Del mismo modo, se muestran el factor de
molecular de la sustancias (Praustnitz, 1969). caracterización K, el peso molecular MW y la
Asimismo, el volumen crítico está asociado con gravedad API experimental y estimada para
parámetros vinculados con el tamaño molecular interpretar los datos en términos de estos
(Praustnitz, 1969). parámetros.
Con respecto a la presión crítica, las deducciones
estadísticas indican que ésta depende de un
compromiso entre la intensidad de las fuerzas
intermoleculares y del tamaño molecular, y tiende
a disminuir cuando aumenta el peso molecular,
como se observa en nuestro caso.
Tabla 2. Propiedades termodinámicas de los crudos pesados a condiciones estándar
Merey Zuata Jobo Cerro negro Hamaca
°API 16 8,1 12 8 9,4
MW ( kg / k mol) 303,8 371,7 476,6 514,1 601,5
K 11,07 10,87 11,47 11,33 11,63
V (m³ / K mol) 0,3174 0,3676 0,4832 0,5079 0,5990
H (K J / Kg) -2152,23 -2158 -2128 -2131,61 -2116
S (K J / Kg °C) 1,0761 1,089 1.518 1,5084 2.128
Cp (K J / K mol°C) 495,166 582,3 772,7 836,243 1001
Cv (K J/ Kmol°C) 486,852 537.9 724,86 827,928 992,6
Ex (K J / Kg) 0,255 0,241 0,256 0,251 0,257
Z 0,0134 0,0155 0,0204 0,0215 0,0253
Los datos indican como el volumen V, las vapor saturado, la región líquido-vapor, la zona
capacidades caloríficas Cp y Cv, y el factor de de líquido, y la zona de vapor.
comprensibilidad Z crecen al aumentar el peso
molecular de los crudos. Tal comportamiento se
debe a que estas propiedades generalmente crecen
a medida que aumenta el peso molecular de las
sustancias. En cuanto a la entalpia H, entropía S,
Exergía Ex y el factor de caracterización K no se
observa una tendencia definida con el peso
molecular. Tampoco se observa una correlación
definida entre las propiedades de la tabla y la
gravedad API de los crudos.
4.4 Diagramas termodinámicos de los crudos
pesados
Debido a limitaciones de espacio, solamente se
visualizarán los diagramas correspondientes al
crudo Cerro Negro, precisando que los diagramas Fig. 2. Diagrama T-V del crudo Cerro Negro.
de todos los crudos examinados en este trabajo
están reportados en (Vielma, 2015).
Con esta aclaración, en las Figs. 2, 3 y 4 se
muestran los diagramas presión - volumen P-V,
temperatura –volumen T-V y presión-temperatura
P-T, respectivamente. La Fig. 2 describe el
diagrama P-V señalando algunas isotermas por
encima y por debajo de la isoterma crítica. Allí
se muestra la línea de líquido saturado, la línea de
Fig.3. Diagrama P-V del crudo Cerro Negro Fig.5. Diagrama P-H del crudo Cerro Negro.
Fig.6. Diagrama T-H del crudo Cerro Negro
Fig.4. Diagrama P-T del crudo Cerro Negro. También, en las Figs. 7 y 8 aparecen los
diagramas presión-entropía P-S y temperatura –
En la Fig.3 se detalla el diagrama T-V mostrando entropía T-S, respectivamente. La Fig. 7
varias isóbaras. También se detecta la línea de representa el diagrama P-S con isotermas
líquido saturado, la línea de vapor saturado, la específicas, y la Fig.8 describe el diagrama T-S
región líquido-vapor, la zona de líquido, y la con la isóbara de 500 kPa. En ambos diagramas
zona de vapor. se aprecia la zona líquido-vapor, la zona de
líquido y la zona de vapor. También se percibe la
En la Fig.4 se muestra el diagrama P-T. Aquí,
línea de líquido saturado y la línea de vapor
también se visualiza la línea de líquido saturado,
saturado
la línea de vapor saturado, la región líquido-vapor,
la zona de líquido, la zona de vapor, y se puede
observar el punto crítico y las condiciones
cricondentérmicas del crudo.
Análogamente, en las Figs. 5 y 6 se presentan los
diagramas presión-entalpía P-H y temperatura-
entalpía T-H, respectivamente. En la Fig.5 se
ilustra el diagrama P-H detallando algunas
isotermas, y en la Fig.6 se describe el diagrama T-
H mostrando la isóbara de 1000 kPa. En ambos
diagramas se visualizan la línea de líquido
saturado, la línea de vapor saturado, la región
líquido-vapor, la zona de líquido y la zona de Fig.7. Diagrama P-S del crudo Cerro Negro
vapor
4000
Merey
Cerro Negro
3000 Jobo
Hamaca 2X
Zuara
2000
1000
Fig.8. Diagrama T-S del crudo Cerro Negro
P (kPa)
0
-7 -2 3 8
Entalpía (kJ7kmol)
Siguiendo un procedimiento similar al utilizado
con el crudo Cerro Negro se construyeron los Fig.10. Diagrama P-H de los crudos Cerro
diagramas termodinámicos de los crudos Hamaca Negro, Zuata, Hamaca 2X, Jobo y Merey.
2X, Jobo, Merey y Zuata los cuales aparecen en
(Vielma , 2015).
4000
Merey
4.5 Diagramas termodinámicos del líquido y
Cerro Negro
vapor saturado de los crudos pesados 3000 Jobo
Los diagramas termodinámicos presión-volumen Hamaca 2X
P-V, presión entalpía P-H y presión-entropía P-S 2000 Zuata
del líquido y vapor saturado de todos los crudos
pesados investigados son ilustrados en las Figs.
1000
9,10 y 11
P (kPa)
4000 0
Merey 500 1500 2500 3500 4500
Cerro Negro Entropía (kJ/kmol °C)
3000 Jobo
Hamaca 2X
Zuata Fig.11. Diagrama P-S de los crudos Cerro
2000 Negro, Zuata, Hamaca 2X, Jobo y Merey.
En el diagrama P-V, mostrado en la Fig.9, la
1000 diferencia más resaltante está relacionada con la
presión crítica de ambos crudos. Allí se detecta
que la presión el crudo Merey duplica,
P (kPa)
0 prácticamente, la presión crítica del crudo Hamaca
0 10 20 30 40 50 2X. Los demás crudos presentan valores de
volumen (m³/kmol) presión críticas intermedios entre los valores
anteriores.
Fig.9. Diagrama P-V de los crudos Cerro Como se explicó previamente, valor de la presión
Negro, Zuata, Hamaca 2X, Jobo y Merey. crítica depende de un compromiso entre las
fuerzas intermoleculares y el tamaño o volumen
de las moléculas del crudo y, en general, salvo
excepciones, en una serie homologa de
hidrocarburos, la presión crítica decrece al
aumentar el peso molecular, como observamos en
nuestro caso (Maxwell,1950; Praustnitz,1969). El
tamaño molecular juega un papel preponderante Considerando el diagrama P-S presentado en la
en este comportamiento. Fig.11 se ve como la entropía del líquido y vapor
saturado de los crudos aumenta y sus curvas se
También, se puede apreciar como a presiones separan notablemente, según la tendencia
menores de 500 Kpa y volúmenes mayores a 15 creciente del peso molecular. También, este efecto
m3 /mol, aproximadamente, las curvas del vapor puede asociarse al peso molecular de los crudos,
saturado de todos los crudos se aproximan en la que al aumentar, hay más campo para el
región de gas ideal. Pero, a presiones mayores a crecimiento entrópico y desorden molecular
500 Kpa, se notan diferencias en las curvas
representativas del vapor y liquido saturado, de tal
manera que el volumen del vapor saturado tiende
a crecer a medida que el peso molecular decrece. 4.6 Tablas termodinámicas de los crudos pesados
En cambio, el volumen del líquido saturado
aumenta a medida que el peso molecular Manejando adecuadamente la información
incrementa. La interacción entre la complejidad suministrada por el simulador es posible
molecular y las fuerzas intermoleculares presentarla en forma tabulada como se detalla a
presentes en el líquido y vapor saturado del crudo continuación para el caso Cerro Negro. En este
explican este comportamiento. sentido, en las Tablas 3,4,5 y 6 se muestran las
propiedades termodinámicas del líquido sub
Con respecto al diagrama P-H ilustrado en la enfriado, del líquido saturado, del vapor saturado
Fig.10 se puede señalar que a una presión y del vapor sobrecalentado del crudo mencionado.
específica, tanto la entalpía del líquido saturado Debido a limitaciones de espacio las tablas
de los crudos como la entalpía del vapor saturado correspondiente a todos los crudo estudiados
de los crudos crecen al aumentar el peso aparecen detalladas en (Vielma, 2015).
molecular de estos, lo cual es razonable porque la
entalpía depende de la capacidad calorífica Cp, la
cual incrementa al aumentar el peso molecular,
como se ve en la Tabla 2.
3. Propiedades termodinámicas del líquido sub enfriado del crudo Cerro Negro
440 °C 475 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
239,143 0,699 -576416,818 1820,523 339,242 0,726 -521221,951 1895,967
826,992 0,698 -576187,870 1820,229 902,067 0,725 -521043,829 1895,630
1414,842 0,697 -575953,309 1819,945 1464,891 0,723 -520857,799 1895,306
2002,691 0,696 -575713,530 1819,670 2027,716 0,722 -520664,497 1894,993
2590,541 0,694 -575468,888 1819,403 2590,541 0,720 -520464,481 1894,690
500 °C 525 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
425,812 0,747 -480801,295 1949,021 526,056 0,770 -439570,704 2001,401
966,994 0,746 -480667,008 1948,648 1042,177 0,769 -439487,485 2000,989
1508,176 0,744 -480522,612 1948,291 1558,298 0,767 -439391,334 2000,596
2049,358 0,742 -480368,998 1947,948 2074,419 0,765 -439283,503 2000,220
2590,541 0,740 -480206,940 1947,617 2590,541 0,763 -439165,067 1999,859
550 °C 575 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
640,828 0,796 -397544,605 2053,130 769,379 0,826 -354733,372 2104,235
1128,256 0,796 -397520,575 2052,675 1224,669 0,826 -354777,606 2103,730
1615,684 0,793 -397479,997 2052,242 1679,960 0,823 -354800,643 2103,253
2103,112 0,791 -397424,636 2051,829 2135,250 0,820 -354804,972 2102,801
2590,541 0,788 -397355,981 2051,435 2590,541 0,817 -354792,658 2102,371
600 °C 625 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
911,249 0,860 -311142,611 2154,743 1065,872 0,901 -266774,225 2204,686
1331,072 0,860 -311265,280 2154,183 1447,039 0,901 -266986,612 2204,065
1750,895 0,858 -311360,753 2153,658 1828,206 0,901 -267164,070 2203,486
2170,718 0,854 -311432,555 2153,163 2209,374 0,896 -267311,603 2202,944
2590,541 0,850 -311483,552 2152,694 2590,541 0,891 -267433,184 2202,433
650 °C 675 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
1230,675 0,951 -221619,733 2254,103 1402,403 1,017 -175663,284 2303,042
1570,641 0,951 -221934,551 2253,413 1699,438 1,017 -176095,093 2302,273
1910,608 0,951 -222204,407 2252,775 1996,472 1,017 -176468,826 2301,570
2250,574 0,950 -222436,426 2252,182 2293,506 1,017 -176794,672 2300,920
2590,541 0,942 -222636,116 2251,627 2590,541 1,015 -177080,283 2300,317
700 °C 725 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
1576,199 1,117 -128873,893 2351,561 1744,781 1,352 -81200,670 2399,738
1829,784 1,117 -129440,070 2350,701 1956,221 1,352 -81924,024 2398,769
2083,370 1,117 -129930,837 2349,924 2167,661 1,352 -82548,584 2397,904
2336,955 1,117 -130360,849 2349,212 2379,101 1,352 -83095,583 2397,121
2590,541 1,117 -130740,748 2348,555 2590,541 1,352 -83579,891 2396,404
750 °C 775 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C) P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
1898,140 1,284 -32565,936 2447,675 2015,369 1,460 17186,304 2495,561
2071,240 1,252 -33479,892 2446,568 2159,162 1,412 16006,353 2494,238
2244,340 1,225 -34262,843 2445,593 2302,955 1,372 15007,464 2493,094
2417,441 1,202 -34945,677 2444,720 2446,748 1,339 14143,553 2492,084
2590,541 1,182 -35549,262 2443,929 2590,541 1,310 13384,327 2491,178
800 °C
P (kPa) v (m3/kmol) h (kJ/Kmol) s (kJ/kmol*°C)
2072,108 1,733 68232,425 2543,604
2201,716 1,650 66632,405 2541,909
2331,324 1,584 65291,182 2540,464
2460,932 1,531 64143,005 2539,207
2590,541 1,486 63143,418 2538,093
Tabla 4. Propiedades termodinámicas del líquido saturado del crudo Cerro Negro
Presión Temperatura Volumen Entalpía Entropía
(K Pa) (°C) (m³ / k mol) (K J / K mol) (k J / K mol*°C)
176,900 412,895 0,684 -618409,975 1772,278
221,213 432,977 0,698 -587547,841 1816,569
275,285 453,921 0,713 -554763,578 1862,263
340,720 475,726 0,730 -520001,200 1909,302
419,168 498,377 0,750 -483219,463 1957,607
512,230 521,851 0,772 -444395,171 2007,081
621,319 546,111 0,797 -403527,805 2057,606
747,443 571,105 0,826 -360645,621 2109,037
890,924 596,759 0,860 -315815,527 2161,200
1052,016 623,131 0,903 -268870,853 2214,197
1225,249 649,591 0,956 -220885,293 2266,780
1405,703 675,949 1,028 -172186,844 2318,624
1586,137 701,987 1,138 -123162,733 2369,382
1756,634 727,473 1,451 -74239,588 2418,710
1864,307 745,041 1,268 -39928,199 2452,572
1954,966 762,143 1,373 -6020,787 2485,480
2021,881 778,711 1,507 27362,510 2517,377
2049,343 789,428 1,618 49281,800 2538,070
2058,388 796,414 1,705 63719,093 2551,598
2059,795 800,279 1,759 71760,786 2559,101
2059,745 800,996 1,769 73257,905 2560,495
2056,999 805,514 1,840 82725,418 2569,295
2049,813 809,965 1,918 92114,361 2577,994
2041,640 813,141 1,979 98854,339 2584,222
Tabla 5. Propiedades termodinámicas del vapor saturado del crudo Cerro Negro.
Presión Temperatura Volumen Entalpía Entropía
(K Pa) (°C) (m³ / k mol) (K J / K mol) (k J / K mol*°C)
202,650 804,369 41,608 121678,635 2619,839
294,854 819,778 28,356 147773,464 2640,982
423,668 833,688 19,381 170837,471 2659,219
598,953 845,251 13,287 189162,343 2673,172
828,394 853,381 9,145 200567,513 2681,086
1111,189 856,796 6,336 202494,790 2680,889
1142,773 856,828 6,105 202052,009 2680,324
1354,338 855,393 4,813 195898,885 2673,859
1555,518 851,249 3,874 184539,983 2663,004
1716,618 845,440 3,244 170367,788 2649,857
1841,253 838,565 2,796 154461,926 2635,257
1933,607 830,968 2,463 137455,495 2619,688
1997,843 822,853 2,206 119700,748 2603,412
2037,823 814,343 2,004 101414,243 2586,585
2041,640 813,141 1,979 98854,339 2584,222
Tabla 6. Propiedades termodinámicas del vapor sobrecalentado del crudo Cerro Negro.
200 kPa 465,6 kPa 731,2 kPa 996,8 kPa 1262,5 kPa 1528,1 kPa 1793,7 kPa
v (m³/kmol) 45,538 18,425 11,012 7,560 5,579 4,317 3,471
870 °C h (kJ/kmol) 235370,072 232208,126 228818,222 225183,346 221310,622 217260,500 213176,954
s (kJ/kmol°C) 2710,086 2700,756 2694,499 2689,204 2684,308 2679,626 2675,155
v (m³/kmol) 47,024 19,217 11,628 8,103 6,085 4,797 3,921
900 °C h (kJ/kmol) 288950,066 286052,823 283012,170 279834,784 276544,951 273192,880 269858,264
s (kJ/kmol°C) 2756,352 2747,250 2741,295 2736,395 2732,003 2727,924 2724,100
v (m³/kmol) 49,454 20,478 12,580 8,916 6,819 5,474 4,551
950 °C h (kJ/kmol) 378563,242 376047,471 373474,827 370861,777 368232,280 365617,717 363054,923
s (kJ/kmol°C) 2831,155 2822,372 2816,808 2812,379 2808,539 2805,077 2801,899
v (m³/kmol) 51,834 21,679 13,463 9,650 7,464 6,058 5,086
1000 °C h (kJ/kmol) 468315,135 466122,353 463920,569 461724,627 459551,750 457420,631 455350,107
s (kJ/kmol°C) 2903,072 2894,548 2889,282 2885,188 2881,714 2878,639 2875,856
v (m³/kmol) 54,173 22,834 14,293 10,327 8,049 6,580 5,559
1050 °C h (kJ/kmol) 557914,189 555999,304 554102,032 552233,331 550404,384 548626,002 546908,013
s (kJ/kmol°C) 2972,102 2963,792 2958,761 2954,919 2951,710 2948,907 2946,396
v (m³/kmol) 56,477 23,950 15,082 10,961 8,590 7,057 5,990
1100 °C h (kJ/kmol) 647068,549 645395,979 643755,309 642153,809 640598,218 639094,473 637647,492
s (kJ/kmol°C) 3038,242 3030,112 3025,271 3021,627 3018,622 3016,023 3013,713
v (m³/kmol) 58,750 25,035 15,838 11,561 9,097 7,501 6,387
1150 °C h (kJ/kmol) 735485,564 734026,156 732605,553 731228,106 729897,493 728616,633 727387,654
s (kJ/kmol°C) 3101,489 3093,512 3088,828 3085,345 3082,500 3080,061 3077,907
v (m³/kmol) 60,997 26,092 16,567 12,133 9,577 7,918 6,759
1200 °C h (kJ/kmol) 822872,476 821602,154 820373,040 819187,360 818046,755 816952,290 815904,480
s (kJ/kmol°C) 3161,841 3153,994 3149,443 3146,092 3143,379 3141,068 3139,040
v (m³/kmol) 63,221 27,127 17,273 12,683 10,034 8,314 7,109
1250 °C h (kJ/kmol) 908935,875 907834,025 906773,042 905753,692 904776,307 903840,832 902946,867
s (kJ/kmol°C) 3219,296 3211,561 3207,122 3203,883 3201,278 3199,074 3197,148
v (m³/kmol) 65,426 28,141 17,959 13,214 10,473 8,691 7,442
1300 °C h (kJ/kmol) 993382,555 992431,782 991519,890 990646,685 989811,699 989014,228 988253,381
s (kJ/kmol°C) 3273,851 3266,214 3261,871 3258,726 3256,213 3254,098 3252,258
5 Conclusiones Firoozabadi,A., Thermodynamics of hydrocarbon
reservoirs, Mc Graw Hill, New York, 1999.
Stryek,R., Vera,J.H., J.Can.Chem.Eng .,
1. Empleando la información experimental 64,April,1986, p.334
obtenida en INTEVEP y en el Laboratorio de
petróleo y catálisis de la ULA se logró obtener Peng,D.Y. and Robinson,D.B., A two Constant
con el simulador Hysys las curvas de destilación equation of State, I.E.C. Fundamentals, (1976)
TBP de los crudos Cerro Negro, Merey, Zuata, ,p.p. 59-64
Hamaca y Jobo las cuales son necesarias para la
simulación de crudos señalados Chang A.F., Pashikanti K., Liu Y.A., Refinery
engineering: integrated process modeling and
2. De los modelos de ecuaciones de estado optimization, Wiley-VCH, New york, 2012.
estudiadas en este proyecto, consideramos que la
ecuación de estado de Peng-Robinson-Strjek-Vera PRO II . KEY WORD MANUAL, 1991.
(PRSV) es la más adecuada para describir el
ASPEN PLUS , Manual, 2004
comportamiento de los crudos pesados de la faja
del Orinoco, porque posee un parámetro de Prausnitz, J.M., Molecular Thermodynamics of
interacción molecular ajustable capaz de Fluid Phase Equilibria, Prentice Hall, Inc., New
considerar las interacciones moleculares de Jersey,1969.
compuestos polares existentes en la matriz
asfalténica de estos crudos. Smith,J.M.,Van Ness, N.C., Abbott, M.M.,
Termodinámica en Ingeniería Química, McGraw
3. Con las herramientas del simulador Hysys y la Hill, Mexico, 2001.
ecuación de estado de Peng-Robinson-Strjek-Vera
PRSV se construyeron los diagramas Vielma, G., Tesis de grado, Facultad de
termodinámicos más comunes P-V, P-T, T-V, P- Ingeniería, Universidad de Los Andes, Venezuela
H, T-H. P-S y T-S, de los crudos Cerro Negro, 2015
Merey, Jobo, Hamaca y Zuata.
4. Utilizando creativamente el simulador se
consiguió elaborar las tablas termodinámicas de
los crudos Cerro Negro, Merey, Jobo, Hamaca y
Zuata con el modelo PRSV
Referencias
Belandria J.I, Procesamiento de los crudos
pesados Ciencia e Ingeniería. Vol, 19, N° 1 (1987)
Montoya D., Rubén A., Destilación al vacío de
crudos pesados. Trabajo de Ascenso. Facultad de
Ingeniería, Universidad de Los Andes, Venezuela
1985
Bellet, A. Molecular distillation: an analytical tool
for heavy and extra heavy crude oil assays,
PDVSA, INTEVEP, Informe Tecnico N°
1998088, 7Th Unitar Conference on heavy crudes
and tar sands, 1998.
Maxwell J.B., Data book on Hydrocarbons, Van
Nostrand, Princeton, N.J., 1950.
FAQ: Modeling heavy oils using Aspen HYSYS,
2010
