Torre
Torre
FACULTAD DE QUÍMICA
Tesis
INGENIERO QUÍMICO
PRESENTA
1
UNAM – Dirección General de Bibliotecas
Tesis Digitales
Restricciones de uso
DERECHOS RESERVADOS ©
PROHIBIDA SU REPRODUCCIÓN TOTAL O PARCIAL
Todo el material contenido en esta tesis esta protegido por la Ley Federal
del Derecho de Autor (LFDA) de los Estados Unidos Mexicanos (México).
Dr. Néstor Noé López Catillo Ing. Gerardo Bonifacio Lemus Morales
Sustentante
2
3
4
Índice
5
Glosario…………………………………………………………………..……………..…8
Capítulo I. Introducción………………………………...…………………………..….15
Objetivos………………………………………………………………………………....18
6
2.7.2. Desventajas de la simulación de procesos………………..…………56
Capítulo V. Conclusiones……………………………………………………………106
Bibliografía………………………………………………………………………..……109
Anexos……………...……………………………………………………………..……113
Anexo A. Índice de Tablas……………...…………………………………………114
Anexo B. Índice de Figuras…………………………………………………...…...118
Anexo C. Índice de Gráficas……………………………………………………….122
7
Glosario
8
Agua amarga “Sour Water”. Agua residual con contenido de ácido sulfhídrico
(H2S), amoníaco (NH3) y compuestos de azufre de baja masa molecular que
proviene de los procesos de refinación. Se obtiene de los domos (parte superior) de
las torres de fraccionamiento de hidrocarburos, tales como primarias, catalíticas,
hidrodesulfuradoras, etc. El agua recolectada en los acumuladores es
posteriormente enviada a plantas de tratamiento de aguas amargas para eliminarle
el ácido sulfhídrico y amoníaco.
Barril de petróleo “oil barrel”. Unidad de volumen basada en la medida del barril
utilizado en la industria del petróleo. Equivale a 158.9873 litros (42 galones de
Estados Unidos).
Para los diferentes tipos de petróleo crudo, cada uno tiene una curva de destilación
única que ayuda a caracterizar los compuestos químicos que se encuentran en él.
Generalmente los compuestos que contienen más átomos de carbono tienen más
alta temperatura de ebullición.
Las pruebas de destilación más usuales para el crudo, sus fracciones o productos,
son la ASTM atmosférica (ASTM D86), la ASTM al vacío (ASTM D1160), la TBP
atmosférica (ASTM D2892).
9
Curva ASTM D86. Esta prueba se lleva a cabo en un aparto simple diseñado para
hervir el líquido a analizar y condensar vapores al mismo tiempo que son generados,
consta de una sola etapa en equilibrio y sin reflujo. La separación de los
componentes es muy pobre.
Curva TBP “true boiling point”. Es una gráfica de los puntos de ebullición de los
componentes puros o de las fracciones de crudo. Se realiza usando un complejo
aparato para destilación por lotes con un gran número de etapas en equilibrio y con
un alto reflujo.
Crudo base aromática “Aromatic base crude oil”. Crudo que contiene grandes
cantidades de compuestos aromáticos de bajo peso molecular y naftenos, junto con
cantidades más pequeñas de asfaltos y aceites lubricantes.
Crudo base asfáltica “Asphalt base crude oil”. Crudos que producen altos
rendimientos de asfalto y aceite combustible pesado.
Crudo base nafténica “Naphtene base crude oil”. Crudo que contiene
principalmente naftenos, esto es, compuestos cíclicos saturados con cadenas
laterales nafténicas y parafínicas, pueden contener mucho material asfáltico. Al
refinarse estos crudos producen aceites lubricantes que se diferencian de los
obtenidos de crudos parafínicos por ser de más baja gravedad y viscosidad, así
como presentar un menor contenido de carbón.
Crudo base parafínica “Paraffin base crude oil”. Crudo de alto contenido en
ceras y fracciones de aceites lubricantes, conteniendo pequeñas cantidades de
naftenos, asfaltos, azufre, nitrógeno y oxígeno.
10
Se consume principalmente en máquinas de combustión interna de alto
aprovechamiento de energía, con elevado rendimiento y eficiencia mecánica. Su
uso se orienta, fundamentalmente, como energético en el parque vehicular
equipado con motores diseñados para combustible diesel, tales como camiones de
carga de servicio ligero y pesado, autobuses de servicio urbano y de transporte
foráneo, locomotoras, embarcaciones, maquinaria agrícola, industrial y de la
construcción.
11
Un gap muestra como debería ser una separación típica en una columna de
destilación atmosférica y esta diferencia numérica debe ser positiva de lo contrario
es un overlap.
Gasóleo “gas oil”. Producto refinado del petróleo cuya densidad es mayor que la
nafta y turbosina, pero menor que la del crudo reducido; generalmente comprende
los hidrocarburos destilados entre 200 y 370 ºC, cuyo rango de pesos específicos
es de 0.820 a 0.890. Esta mezcla de hidrocarburos tiene dos usos principales:
Combustible para pequeñas máquinas diesel y para hornos o calentadores, de
donde toma sus nombres populares, diesel y aceite para hornos (furnace oil).
Nafta “naphtha”. Primer corte o fracción de la destilación del petróleo, cuyo rango
de destilación comprende de los 30 a los 232 ºC. En este tipo de clasificación se
denomina nafta ligera a la fracción de la destilación comprendida hasta los 101 ºC,
y nafta pesada a la fracción comprendida entre 101 ºC y 232 ºC.
12
En la destilación atmosférica del crudo, la fracción considerada en este rango
comprende la gasolina primaria la cual es estabilizada, hidrodesulfurizada y
después de estos procesos sirve como carga a las plantas reformadoras de nafta;
con el objetivo principal de incrementar su octanaje, utilizándose como componente
básico de las gasolinas comerciales.
Punto de corte “TBP cut point”. Temperatura en la curva TBP del crudo que
representa los límites (superior e inferior) de una fracción a ser producida.
13
Punto de inflamación “flash point”. Temperatura a la cual un hidrocarburo libera
vapores en cantidad suficiente para permitir la combustión.
14
Capítulo I
Introducción
15
Introducción
El objetivo de esta tesis es establecer los criterios de diseño para llevar acabo el
diseño de la torre de destilación atmosférica que procesará 250,000 barriles por día
de crudo Istmo y Maya, enfocada particularmente a la ingeniería de proceso, que
es parte de la ingeniería básica, relacionada con los balances de materia y energía,
las especificaciones en los productos finales, el grado de separación entre las
diferentes fracciones del crudo Istmo y Maya, entre otros. Posteriormente se evalúa
de forma aproximada las dimensiones de la torre de destilación con base a las
especificaciones de los internos.
Una vez establecidos los criterios de diseño, según lo reportado en la literatura como
Petroleum Refinery Distillation de Watkins[2] y Handbook of Petroleum Processing
de David Jones y Peter Pujadó[3] se genera el diseño de la torre de destilación
atmosférica en el simulador de procesos Aspen HYSYS de la compañía AspenTech.
16
Aspen HYSYS es hoy en día una herramienta computacional indispensable para la
simulación de cualquier planta química que desee llevarse a operación con la
ventaja de tener un mejor conocimiento del proceso a través de las variables críticas
y así poder optimizar el proceso.
Los capítulos tratados a lo largo de esta tesis se enfocan en aspectos del proceso
de destilación atmosférica del petróleo. Estos mencionan conceptos básicos para
una mejor comprensión del caso de estudio, abarcan desde qué es el petróleo hasta
la descripción del proceso de destilación atmosférica, y explica cómo se llevó a cabo
la simulación del proceso para finalmente presentar y evaluar los resultados.
17
Objetivos
Objetivo General
Establecer los criterios para el diseño de una torre de destilación atmosférica que
procesará 250,000 barriles por día de crudo Istmo y Maya con la ayuda del
simulador Aspen HYSYS.
Objetivos Específicos
1. Corroborar que la curva TBP y de gravedad API del crudo Istmo y Maya
generada solo con pseudocomponentes sea igual a las curvas de la técnica
assay.
2. Hacer la simulación de la torre de destilación atmosférica bajo los criterios de
diseño según lo reportado en Petroleum Refinery Distillation de Watkins[2] y
Handbook of Petroleum Processing de David Jones y Peter Pujadó[3] en el
simulador Aspen HYSYS.
3. Lograr que se cumpla el balance materia y energía en la torre de destilación
atmosférica.
4. Obtener el grado de separación entre las fracciones adyacentes para que los
crudos Istmo y Maya se encuentren en valores permisibles según lo
reportado en Petroleum Refinery Distillation de Watkins[2] y Handbook of
Petroleum Processing de David Jones y Peter Pujadó[3].
5. Conseguir que los productos del crudo Istmo y Maya estén dentro de las
especificaciones que se requieren en México hoy en día.
6. Que las curvas TBP de cada fracción del crudo Istmo y Maya se traslapen
sobre toda la curva TBP generada a partir de pseudocomponentes.
7. Lograr el equilibrio entre la eficiencia de separación y recuperación de
energía en la torre de destilación atmosférica.
18
Capítulo II
Marco Teórico
En este capítulo se establecerá de manera general los fundamentos teóricos y solo
está incluida la información necesaria para entender aspectos de los que se
hablarán en capítulos posteriores.
19
2.1. Definición de petróleo
Las moléculas de hidrocarburos son más complejas a medida que aumenta el punto
de ebullición de las mismas. En un crudo están presentes en mayor o menor
proporción hidrocarburos parafínicos, aromáticos y nafténicos. Los crudos tienen
diferente composición en cuanto al tipo y cantidad de familias químicas
hidrocarbonadas de las que están formadas, las características físicas variarán en
función de esta composición[5].
Son miles los compuestos químicos que constituyen el petróleo, y entre muchas
otras propiedades, estos compuestos se diferencian por su volatilidad (dependiendo
de la temperatura de ebullición). Al calentarse el petróleo, se evaporan
preferentemente los compuestos ligeros (de estructura química sencilla y baja masa
molecular), de tal manera que conforme aumenta la temperatura los componentes
más pesados van incorporándose al vapor[5].
20
ligero y viceversa[6]. A continuación se mencionan los dos crudos utilizados en este
trabajo con sus respectivas características:
En la destilación TBP del crudo Istmo tiene buen rendimiento de destilados hasta
538 °C (81.33% Vol.), por su viscosidad el crudo reducido hasta 538 ºC, puede ser
utilizado como materia prima para la obtención de combustóleo pesado con
contenido de azufre menor a 4.0% peso[7].
21
600
Temperatura (°C) a 760 mm Hg
500
400
300
200
100
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
En la destilación TBP del crudo Maya tiene bajo rendimiento de destilados hasta
538 °C (64.04% Vol.), por su viscosidad el crudo reducido hasta 538 ºC, puede ser
utilizado como materia prima para la obtención de combustóleo pesado con alto
contenido de azufre menor a 5.5% peso[7].
22
El crudo Maya presenta las características que se muestran Tabla 2 y Figura 2.
600
500
Temperatura (°C) a 760 mmHg
400
300
200
100
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Voluemen Vaporizado
23
Para cierto valor de temperatura llamado T1 (Figura 3) está definido el corte para
la fracción de nafta, la curva TBP muestra el rendimiento que se obtendría; para el
crudo Istmo se obtiene un rendimiento del 24% de nafta mientras que para el crudo
Maya sólo se obtiene el 16%; al mismo tiempo, la comparación de volumen
recuperado para un mismo valor de temperatura muestra las variaciones en el
rendimiento entre el crudo Istmo y Maya, las temperaturas a las que se debe
someter un crudo pesado son muy altas debido a que está compuesta por
hidrocarburos de gran masa molecular y la cantidad recuperada de productos de
alto valor agregado (gasolinas y diesel) son menores con respecto al crudo Istmo;
El crudo Istmo tiene mayor cantidad de productos de alto valor agregado y no
requiere ser sometido a altas temperaturas[8].
24
600
550 Istmo
500 Maya
450
Temperatura (°C) a 760 mmHg
400
350
300
250
200
T1
150
100
50
0
0 10 20 30 40 50 60 70 80 90 100
% Volumen Vaporizado
2.3. Caracterización
25
petróleo de composición desconocida, el punto de ebullición puede representarse
mediante una curva de temperatura frente al porcentaje en volumen de mezcla
vaporizada a continuación se muestran dos métodos muy frecuentes para la
caracterización del petróleo:
26
de la curva ASTM D86. Por tanto, el rango de puntos de ebullición basado en la
curva ASTM D86 es menor que el rango de la curva TBP. Esta diferencia se debe
a las distintas condiciones de operación en las que se dan los dos ensayos, una de
ellas es porque la curva ASTM D86 produce un gran fraccionamiento, separando
mejor los componentes ligeros es decir con menor masa molecular, en el primer
tramo. En el segundo tramo aparece por arriba la curva TBP por la eficiencia en
fraccionar los componentes de mayor masa molecular[6].
260
TBP
ASTM D86
240
Temperatura (°C) a 760 mmHg
220
200
180
160
140
120
0 10 20 30 40 50 60 70 80 90 100
% Volumen Vaporizado
27
temperatura aumenta, pero el efecto de la presión en la densidad de los líquidos es
generalmente insignificante[6].
Dado que las condiciones estándar adoptadas por la industria del petróleo son 15
°C y 1 atm, el peso relativo de los hidrocarburos líquidos se da generalmente en
esas condiciones. A esa temperatura de referencia, la densidad del agua líquida es
de 0.999 g/cm3. Por tanto, para un hidrocarburo o para una fracción de petróleo, el
peso específico se define como lo muestra la ecuación 2:
En los primeros años de la industria, el Instituto Americano del Petróleo (API) definió
la densidad API (grados API) para cuantificar la calidad de los productos del petróleo
y del crudo. La densidad API se define como lo muestra la ecuación 3:
141.5
°𝑨𝑷𝑰 = − 131.5 (3)
𝑆𝐺 (𝑎 60 °𝐹)
28
Los hidrocarburos líquidos con un peso relativo bajo tienen una densidad API
elevada. Los hidrocarburos aromáticos tienen mayor peso específico (menor
densidad API) que los de las parafinas [6].
3
√𝑇𝑀
𝑲𝑼𝑶𝑷 = ( ) (4)
𝑆𝐺
Donde:
Los valores de KUOP de los crudos oscilan entre 12-13 para crudos parafínicos y
entre 10.5-11 para crudos nafténicos y aromáticos, permite estimar la composición
química del petróleo y las diferentes fracciones[5].
2.3.4. Pseudocomponentes
29
Estos compuestos hipotéticos constituyen una poderosa herramienta que mediante
métodos predictivos permite, con pocos datos conocidos (tales como masa
molecular, temperatura, densidad, etc.), predecir las propiedades de las fracciones
del crudo tales como la presión crítica (Pc), temperatura crítica (Tc) y factor
acéntrico (w) y así tener una aproximación del comportamiento del sistema, pueden
ser definidos a partir de la curvas de destilación TBP o ASTM D86[3].
Mientras que el punto de corte es una temperatura ideal utilizada para definir el
rendimiento de una fracción, los puntos finales son las temperaturas reales de una
fracción. Ningún proceso tiene la capacidad de separar perfectamente los
componentes de una fracción, cuando dos fracciones se separan en un proceso
algunos de los componentes más ligeros permanecen en la fracción más ligera
adyacente y del mismo modo para los componentes pesados, por lo tanto, la
temperatura inicial de ebullición real de la fracción será más baja que el punto de
corte inicial y su temperatura final de ebullición será más alta que su punto de corte
final[3].
30
Figura 5. Ejemplo del punto de ebullición medio para un pseudocomponente sobre
la curva de destilación TBP para un crudo[3].
31
Temperaturas de ebullición
32
2.4. Aspectos generales sobre destilación
33
localizada arriba del plato de alimentación. La sección de agotamiento permite
mejorar la recuperación del producto más pesado y se encuentra localizada debajo
del plato de alimentación.
Son aquellas que definirán el diseño y las condiciones de operación del proceso;
para el diseño de la torre de destilación atmosférica debe tomarse en cuenta todas
las variables involucradas, principalmente las de mayor impacto, ya que son las que
servirán para el diseño tales como:
34
2.6. Criterios y variables del proceso para el diseño de la torre de
destilación atmosférica
35
a que, desde el punto de vista de equilibrio, el líquido de la extracción de calor es
exterior a la zona en la que se emplea[2].
Este método es práctico para remover calor y utilizado para el precalentamiento del
crudo antes de entrar al horno pero si se extrae demasiado calor de la torre afecta
de manera negativa al fraccionamiento[2].
Se observa que la turbosina, cuyos puntos típicos de corte son entre 180 y 230 °C
en ASTM D86[5], está en realidad compuesta por hidrocarburos cuyas temperaturas
de ebullición varían entre aproximadamente 130 y 280 °C, es decir que la turbosina
contiene un porcentaje de nafta y diesel.
36
Si T>0; existe una gap
1150
1100 Nafta
1050 Turbosina
1000 Diesel
GOPA
950
Crudo Reducido
900
Curva ASTM D86, Temperatura (°C)
850
800
750
700
650
600
550
500
450
400
350
300
250
200
150
100
50
0
0 10 20 30 40 50 60 70 80 90 100
% Volumen Vaporizado
37
280
275
270
265
260
255
250
245
240
Temperatura (°C)
235
230
gap
225
220
215
210
205
200
195
190 Turbosina
185 Diesel
180
22.0 23.0 24.0 25.0 26.0 27.0 28.0 29.0 30.0 31.0 32.0 33.0 34.0
% Volumen Vaporizado
Figura 10. Eficiencia de fraccionamiento a través del gap y overlap entre la fracción de turbosina y diesel sobre la curva de destilación
ATMS D86[2].
38
Un ejemplo de un fraccionamiento eficiente es como el que se muestra en la Figura
10 entre la turbosina y diesel ya que existe un gap de 28 °C sobre la destilación en
ASTM D86, en la literatura, Petroleum Refinery Distillation[2] o Handbook of
Petroleum Processing[3] mencionan el grado de separación permisible entre las
fracciones adyacentes, mencionados a continuación:
39
Figura 11. Curva de Packie que relaciona la dificultad de separación y la
capacidad de separación del fraccionamiento entre el domo y el producto lateral
superior[3].
40
Figura 12. Curva de Packie que relaciona la dificultad de separación y la
capacidad de separación del fraccionamiento entre los productos laterales
adyacentes[3].
41
2.6.3. Tipo de condensador en la sección del domo
La corriente vapor de domos de la torre atmosférica está formada por nafta, reflujo
y vapor de agotamiento adicionado en los fondos de la torre y en los agotadores
laterales. El arreglo de condensación es en dos etapas (Figura 13) y servirá para
tener una máxima recuperación de calor en el intercambio de calor con el crudo en
el tren de precalentamiento. El primer condensador es un condensador parcial, el
líquido condensado se separa de la fase vapor, en un tanque acumulador y retorna
como reflujo al plato inmediato del domo (Plato 1) de la torre de destilación[10].
La fase vapor de la primera etapa de condensación se enfría con una parte del
crudo, aire y agua de enfriamiento. La nafta y agua condensadas se separan en otro
tanque acumulador y la temperatura del producto destilado tiene valores típicos
entre 32 y 38 °C[10].
Esta configuración es adecuada debido al gran flujo que se tiene en la parte del
domo; su propósito es mantener un diseño en los intercambiadores de calor y
tanques acumuladores en dimensiones manejables, en esta configuración el primer
tanque acumulador que pertenece al reflujo es considerado como un plato teórico
adicional debido a que el reflujo está en equilibrio con la fase vapor del producto [3].
Las caídas de presión entre el domo y el condensador son entre 5-7 psi tomando
en cuenta las caídas de presión por accesorios cuando el sistema está limpio pero
la caída de presión normal está en el orden de 3 psi[2].
42
2.6.4. Presión de operación de la columna de destilación
La presión de operación de la torre atmosférica debe ser la más baja, limitada por
el tanque acumulador del reflujo de la torre ya que en muchos casos, la corriente
del domo no puede ser totalmente condensada a una baja presión[2].
Muchas de las torres tienen de 25 a 35 platos reales entre la zona flash y el domo,
usando el número de platos y sabiendo la caída de presión por plato (0.1 a 0.2 psi),
entonces se conoce la caída de presión total entre la zona flash y el domo[2].
Trabajar a alta presión requiere una mayor temperatura a la salida del horno, y por
lo tanto un mayor consumo de energía lo que conduce a la posibilidad de formación
de coque en los tubos del horno si se excede de la temperatura permisible. De
hecho, trabajar a altas presiones hace el fraccionamiento más difícil[2].
Esta variable debe ser controlada con precisión, ya que determina el grado de
vaporización del crudo y de la cantidad de calor suministrada al calentador. La
temperatura de entrada del crudo en la columna está relacionada con las
características del propio crudo, además de la cantidad y calidad de las fracciones
que se producen; porque toda la energía se proporciona a través del
precalentamiento del crudo, cualquier variación de temperatura modifica el equilibrio
entre las fases y el perfil de la temperatura a lo largo de la columna, tal como para
causar una variación en las características de los productos[11].
Una baja temperatura a la salida del horno produce una cantidad más pequeña de
material vaporizado y en consecuencia una reducción de la cantidad de destilados,
mientras que una temperatura demasiado elevada provoca reacciones de craqueo
con la formación de gas y depósitos de carbono, en particular en los tubos del horno,
lo que hace necesario el paro de la planta con mayor frecuencia para limpiar el
horno[11].
43
La temperatura máxima permisible del crudo es en un intervalo de 343 a 357 °C ya
que la mayoría de los crudo de base parafínica y nafténica comienza a sufrir
descomposición térmica como es el caso del crudo Istmo y Maya[2].
2.6.6. Overflash
44
Tabla 5. Especificaciones de los diferentes productos en la torre de destilación atmosférica para la
nueva refinería (Tula, Hidalgo)[12].
Gasóleo
Residuo
Propiedad Unidades Método Nafta Turbosina Diesel Pesado
Atmosférico
Atmosférico
Densidad API °API 64.27 45.76 32.96 20.90 7.09
ASTM
Peso Específico 0.7228 0.7983 0.8604 0.9284 1.021
D1298
Factor de
12.150 11.838 11.664 11.450 11.621
Caracterización
Azufre % Peso 0.191 0.587 1.85 2.89 5.014
Nitrógeno Básico % Peso ---- ---- 0.012 0.038 0.0966
Nitrógeno Total % Peso 0.04 0.115 0.496
Níquel ppm ---- ---- ---- ---- 89
Vanadio ppm ---- ---- ---- ---- 444
Sodio ppm ---- ---- ---- ---- 26.663
Fierro Ppm ---- ---- ---- ---- 14.166
Carbón Conradson % Peso ---- ---- ---- 0.0032 19.581
Temperatura de 63.7 81.8
°C ---- 139.8 ----
Inflamación (38 mín.) (41 mín.)
Temperatura de ASTM
°C ---- -47.3 ---- ---- ----
Congelamiento D2386
Temperatura de
°C -44.4 -9.1
Niebla
Viscosidad ASTM
cSt ---- ---- ---- ---- 788,070
Cinemática a 50°C D445
Viscosidad ASTM
cSt ---- ---- ---- ---- 2,177
Cinemática a 100°C D445
Destilación, % Vol. ASTM
°C ASTM D86 ASTM D86 ASTM D86 ASTM D 1160
a TIE D86
TIE 44.1 125.7 147.4 289.5 297.7
5 70.3 157.1 224.6 331.2 374.9
244.8
10 86.3 170.2 350.6 413.8
(275 máx.)
50 127.7 196.3 281.4 409.1 622.2
343.3
90 165.42 217.3 447.5 ----
(345 máx.)
95 169.7 222.5 354 452.5 ----
TFE 175.0 máx. 230.0 máx. 365 457.8 ----
45
2.6.8. Reflujo
El calor que no se elimina por medio de los productos que salen de la columna se
debe extraer mediante la extracción del líquido en las secciones apropiadas de la
columna y el envío de vuelta a la columna después de intercambiar calor con otras
corrientes o sean enfriadas por medio de agua o aire[11].
46
La distribución del calor extraído a lo largo de la columna está determinada por el
grado de fraccionamiento necesario, así como por la cantidad de calor que se retira.
Además de la temperatura, el flujo de las diversas extracciones de calor en la
columna también tiene un efecto en el fraccionamiento[11].
47
En cada punto a lo largo de la columna atmosférica, y especialmente en el punto de
extracción en la columna de una corriente o un corte lateral, la temperatura del
líquido refleja la temperatura de ebullición del líquido en este punto exacto. Sin
embargo, un aumento en el flujo de líquido extraído de la columna conduce a una
reducción del reflujo interno, y como resultado, el líquido se vuelve más rico en
componentes menos volátiles y, en consecuencia, eleva su temperatura de
ebullición en la sección en cuestión. A la inversa, una reducción del flujo de la
extracción causa que más líquido permanezca en la sección de la columna
inmediatamente debajo del punto de extracción, con la consiguiente existe una
reducción de la temperatura[11].
Como regla general se puede afirmar que, cuando el flujo de extracción aumenta,
el punto de ebullición de los productos extraídos se eleva, mientras si se reduce,
disminuye el punto de ebullición, por lo que hace que se vuelva más rica en los
componentes más volátiles[11].
Las torres de platos son cilindros donde el líquido y el gas se ponen en contacto en
forma de pasos sobre los platos. El líquido entra en la parte superior y fluye en forma
descendente por gravedad y fluye en el plato a través de orificios que pueden ser
de distintos tipos; entonces el vapor se desplaza en dirección contraria, es decir
hacia arriba, tal como lo indica la Figura 7. Cada plato en la torre es una etapa,
puesto que sobre el plato se ponen en contacto íntimo los fluidos donde ocurre la
difusión interfacial y los fluidos se separaran[9].
48
Tabla 6. Platos requeridos en las diferentes secciones de la torre de
destilación atmosférica[2].
Sección Número de Platos
Nafta Ligera a Nafta Pesada 6-8
Nafta Pesada a Destilado Ligero 6-8
Destilado Ligero a Destilado Pesado 4-6
Destilado Pesado a Aceite Atmosférico 4-6
Zona Flash a primer plato extraído 3-4
Vapor y Reboiler en sección de Agotamiento 4
Existe una gran variedad de tipos de platos para las torres de destilación, como
válvula, perforados, burbujeo o tipo chimenea, la selección de cada plato es
determinada por varios factores tales como: el proceso, costo, diseño mecánico,
mantenimiento, capacidad, eficiencia, caída de presión, corrosión, etc.[13].
Solamente se mencionará el plato tipo válvula ya que fue utilizado para el diseño de
la torre de destilación en la simulación.
49
2.6.12.1.1. Plato tipo válvula
Los platos válvula (Figura 15) son muy flexibles con un rango de operación típico
del 30 al 85% de inundación; además son de 30 a 50% más económicos que los de
burbujeo de un mismo diámetro a causa de su construcción más simple. Presentan
capacidades 20% más altas para el vapor que de los platos de cachucha[13].
2.6.12.2. Eficiencia
En la operación real las corrientes de vapor y líquido que salen de un plato no están
completamente en equilibrio, para lograr el mismo grado de separación, más platos
50
tendrán que ser añadidos para compensar esta deficiencia, el número real de platos
necesarios es ajustado con el concepto de eficiencia, una aproximación general
como lo muestra la ecuación 6:
Debido a que los simuladores trabajan con etapas ideales y el diseño quiere llevarse
a condiciones reales se considera una eficiencia en cada una de las diferentes
secciones de la columna (Tabla 7), estos valores son lo que se presentan hoy en
día en cualquier torre de destilación atmosférica[15]:
51
Figura 16. Diámetro y altura de una torre de destilación[16].
Los platos deben estar separados a una distancia tal que permita la separación del
líquido y el vapor, tanto en el bajante como en el área activa. El espaciamiento entre
platos tiene un efecto importante, tanto como en la inundación como en el arrastre.
Bajo las condiciones limitantes de inundación, hay una relación recíproca entre el
espaciamiento y el diámetro, de tal forma que se requiere modificar la inundación
dejando fijo el espaciamiento, tendrá que modificarse el diámetro. El espaciamiento
más común es de 45.72 cm[13].
2.6.12.6. Inundación
Ocurre cuando los flujos del líquido o del vapor resulten mayores que la capacidad
del equipo (Figura 17). No puede decirse que la torre está inundada por exceso de
líquido o por demasiado vapor, la inundación es una combinación de los dos flujos.
Una torre se inunda principalmente cuando se exceden los valores máximos
52
permitidos de la retención del líquido en el bajante y del arrastre, la inundación de
una torre se caracteriza porque se presenta un súbito incremento en la caída de
presión y una marcada disminución de la eficiencia[13].
53
Figura 18. Bajante “Downcomer backup” en torre de destilación atmosférica[14].
54
2.7. Simulación de procesos
55
tipo de análisis “What if” son muy usados en la industria para el diseño y
el rediseño de sistemas, por tanto se convierte en una herramienta muy
importante para el ingeniero de proceso.
Dentro de la industria, los simuladores son útiles en investigación y
desarrollo al agilizar los ensayos en laboratorios y en plantas piloto
porque permite predecir resultados o rangos de trabajo óptimo lo cual
ahorra mucho tiempo en la experimentación; y en la etapa de diseño, al
ofrecer una manera rápida de diseñar un proceso para establecer
comparaciones entre diferentes alternativas[20].
Los simuladores son muy útiles en el entrenamiento de personal. Los
operadores de las plantas pueden aprender de sus errores y de este
modo, aprender a operar mejor[18].
56
2.8. Simulador Aspen HYSYS
Las aplicaciones de Aspen HYSYS son muy variadas e importantes, tales como en
los procesos de refinación de petróleo, el transporte de fluidos por tubería, procesos
de endulzamiento de gas, procesos químicos, etc.
57
Capítulo III
Simulación de la Torre de
Destilación Atmosférica en
Aspen HYSYS
En este capítulo se describe de manera detallada como se llevó a cabo la simulación
de la torre de destilación atmosférica en Aspen HYSYS, la cual incluye: descripción
del proceso, caracterización del crudo Istmo y Maya, especificaciones de equipos y
los diferentes productos, etc.
58
3.1. Descripción general del proceso de destilación atmosférica
Figura 19. Diagrama de flujo de proceso para el proceso de destilación atmosférica implementado en
Aspen HYSYS.
La destilación del crudo se efectúa en una columna de platos que cuenta con una
zona de rectificación, formada por varias secciones, de las que se extraen diferentes
fracciones según su masa molecular: gas y nafta que se obtienen como productos
de domos, turbosina, diesel y gasóleo pesado atmosférico que se obtienen como
extracciones laterales, contando los tres primeros, con columnas de agotamiento
con vapor de agua sobrecalentado (2.5 kg/cm2 man. y 343 °C), para obtener las
especificaciones requeridas. Además de las extracciones de productos, la torre de
59
destilación cuenta con dos extracciones de calor “Pumparounds” para la extracción
de calor“[10] (Figura 20).
60
especificaciones de producto de la sección correspondiente. El arreglo de
condensación es en dos etapas para tener una máxima recuperación de calor en
favor del crudo. El primer condensador es un condensador parcial, el líquido
condensado se separa de la fase vapor a 118 °C y 0.7 kg/cm2 man., en el
acumulador de la primera etapa de condensación y se retorna como reflujo al plato
número 1 de la torre de destilación, la segunda etapa de condensación es total, la
nafta y agua condensados se separan a 38 °C y 0.7 kg/cm2 man., en el acumulador
de la segunda etapa de condensación[10].
61
Cantidad de líquido de las extracciones de calor (Pumparounds).
Cantidad de vapor en los agotadores.
Con ayuda del simulador Aspen HYSYS, se pudo simular la torre de destilación
atmosférica, para ello se realizaron unas series de pasos descritos a continuación:
Cada fracción del crudo está formado por una gran cantidad de componentes por
consiguiente su caracterización resultaría muy dificíl si se hiciera tomándolos como
componentes puros por lo que es necesario definir pseudocomponentes que son
mezclas de hidrocarburos en función del punto de ebullición medio (mid boiling
point) y la gravedad API; por lo tanto, las curvas de destilación TBP y de gravedad
API son usadas para definir dichos pseudocomponentes y de esta manera poder
estimar sus propiedades.
Se utilizaron datos del Instituto Mexicano del Petróleo que son curvas elaboradas
con la técnica assay que definen las propiedades de ambos crudos como la curva
de destilación TBP (Figura 21 y 23), gravedad API y contenido de azufre (Figura
22 y 24)[7]; una mayor cantidad de información para la caracterización del crudo, se
verá reflejada en mayor confiabilidad de los resulatdos.
62
540
520
500
480
460
440
420
400
380
360
340
320
Temperatura (°C)
300
280
260
240
220
200
180
160
140
120
100
80
60
40
20
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Figura 21. Curva de destilación TBP de crudo 100% Istmo elaborada con la
técnica assay[7].
63
100 2.5
90
80 2
70
60 1.5
Gravedad API
50
40 1
30
20 0.5
10
Gravedad API
Contenido de Azufre
0 0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80
%Volumen Vaporizado
Figura 22. Curva de gravedad API y contenido de azufre de crudo 100% Istmo elaboradas con la técnica assay[7].
64
540
520
500
480
460
440
420
400
380
360
340
Temperatura (°C)
320
300
280
260
240
220
200
180
160
140
120
100
80
60
40
20
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Figura 23. Curva de destilación TBP de crudo 100% Maya elaborada con la
técnica assay[7].
65
100.0 3.5
90.0
3
80.0
70.0 2.5
60.0
Gravedad API
50.0
1.5
40.0
30.0 1
20.0
0.5
10.0 Gravedad API
Contenido de Azufre
0.0 0
0 5 10 15 20 25 30 35 40 45 50 55 60 65
% Volumen Vaporizado
Figura 24. Curva de gravedad API y contenido de azufre de crudo 100% Maya elaboradas con la técnica assay[7].
66
El primer paso antes de la definición de los pseudocomponentes en Aspen HYSYS,
fue extrapolar la curva de destilación TBP, la de gravedad API y la del contenido de
azufre de la técnica assay debido a que ninguna de ellas llega al 100% de volumen
vaporizado a causa de las limitaciones que se tienen en la caracterización de los
crudos en el laboratorio; la extrapolación de la curva TBP de ambos crudos fue
elaborada con el método númerico spline cúbico en el simulador de procesos PRO/II
debido a que este método realiza una mejor extrapolación de la curva TBP.
67
780
760
740
720
700
680
660
640
620
600
580
560
540
520
Temperatura (°C)
500
480
460
440
420
400
380
360
340
320
300
280
260
240
220
200
180
160
140
120
100
80 Mid Boiling Point de
60 pseudocomponente Hypo2006*
40
20
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50 52 54 56 58 60 62 64 66 68 70 72 74 76 78 80 82 84 86 88 90 92 94 96 98100
%Volumen Vaporizado
Figura 25. Ejemplo de la caracterización del crudo por pseudocomponentes con base a la curva de destilación TBP[7].
68
Una vez que todas las curvas para ambos crudos fueron divididas en
pseudocomponentes se introdujeron los valores leídos de las curvas (TBP y
gravedad API) en la sección de Hypotheticals del simulador, donde solo requiere
una mínima cantidad de información (masa molecular, propiedades críticas,
densidad o temperatura) para poder definir cada pseudocomponente y así poder
estimar las demás propiedades de cada uno, al menos 2 datos requiere el simulador
para calcular las propiedades restantes (Figura 26).
69
Este método de caracterización del crudo solo consiste en dividir la curva de
destilación TBP y gravedad API en rangos estrechos, siendo cada rango la
representación de un pseudocomponente, en general, mientras más
pseudocomponentes se utilicen, mayor será la precisión y la confiabilidad de la
estimación de las propiedades y por ende los rendimientos de los productos.
70
Figura 27. Definición del paquete termodinámico en Aspen HYSYS.
71
Tabla 9. Corriente de alimentación para crudo Istmo y
Maya en Aspen HYSYS[10].
Temperatura (°C) 28
Presión (kg/cm2 man.) 3.97
Flujo Volumétrico (MBPD Std.) 250
Tabla 11. Condiciones de entrada de corrientes de vapor sobrecalentado para agotadores y torre
atmosférica para crudo Istmo y Maya en Aspen HYSYS.
Vapor de
343 2.47 972 856 29
Agotamiento
Vapor de
343 2.47 82 84 3
Turbosina
Vapor de Diesel 343 2.47 259 173 3
Vapor de GOPA 343 2.47 56 35 3
Vapor Total
--- --- 1,369 1,148 ---
Requerido
72
3.4.2. Torre de destilación atmosférica
La torre de destilación atmosférica, agotadores y extracciones de calor
“Pumparounds” fueron construidos con las siguientes especificaciones (Tabla 12,
13 y 14) provenientes de la recopilación de los criterios de diseño según lo reportado
en la literatura[2] y [3]:
Tabla 13. Especificación de agotadores para crudo Istmo y Maya en Aspen HSYSY.
Número de Etapas Etapa de Etapa de
Teóricas salida Retorno
Agotador de Turbosina 3 10 9
Agotador de Diesel 3 17 16
Agotador de GOPA 3 24 23
73
Tabla 14. Especificación de las extracciones de calor para
crudo Istmo y Maya en Aspen HYSYS.
Etapa de Etapa de
salida Retorno
PA 1 (Nafta-Turbosina) 10 9
PA 2 (Turbosina-Diesel) 17 16
74
3.4.3. Especificaciones para la convergencia de la torre de
destilación atmosférica
Se añadió una cantidad de especificaciones con el fin de tener los grados de libertad
necesarios para logar la convergencia de la torre de destilación (Figura 30), además
de proponer buenos estimados iniciales como la temperatura de retorno y cantidad
de líquido en las extracciones de calor, la temperatura final de nafta y turbosina,
etc.; para que fuese rápida y el simulador no tuviera problemas en la convergencia,
ya que si no se tenían buenos estimados el simulador empieza a generar ciclos y
no lograría resolver la columna. En la Tabla 15 se muestran las especificaciones
dadas al simulador para que se tuvieran los grados de libertad necesarios para
ambos crudos:
Tabla 15. Especificaciones para la convergencia de la torre de destilación para crudo Istmo y
Maya en Aspen HYSYS.
Especificación Variable Unidad Crudo Istmo Crudo Maya
Temperatura de
PA1 °C 140 140
Retorno
(Nafta-Turbosina)
Flujo MBPD Std. 172 144
Temperatura de
PA2 °C 160 155
Retorno
(Turbosina-Diesel)
Flujo MBPD Std. 113 99
Temperatura de
TFE en ASTM D86 de Nafta °C 174 177
Producto
TFE en ASTM D86 de Temperatura de
°C 232 234
Turbosina Producto
Temperatura al 90% en ASTM Temperatura de
°C 345 345
D86 de Diesel Producto
Overflash Flujo MBPD Std. 12.5 12.5
75
Figura 30. Especificaciones para la convergencia de la torre de destilación en
Aspen HYSYS.
Finalmente en Tools/Utilities/Tray sizing se especificaron las características de los
internos dela torre de destilación atmosférica (Tabla 16) donde la operación se
realizó con platos, con el fin de tener un valor aproximado de la dimensiones de la
columna y así conocer, a rasgos muy generales, las limitaciones desde el punto de
vista de la hidráulica de la torre de destilación, haciendo énfasis al parámetro de
inundación.
Tabla 16. Especificaciones de internos de la torre
de destilación atmosférica[13].
Máxima inundación (%) 85
Tipo de plato Válvula
Tipo de rejilla Glitsch
Número de pasos del flujo 4
Espaciamiento en cada plato (m) 0.6096
Máximo downcomer backup (%) 60
Espaciamiento del bajante (mm) 50.80
Tabla 16. Especificaciones de internos de la torre de destilación atmosférica [13].
76
Capítulo IV
Análisis de Resultados
77
4.1. Resultados comparativos para ambos crudos
78
740
Assay
690
HYPO
640
590
540
490
Temperatura (°C)
440
390
340
290
240
190
140
90
40
-10
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
1400
1300 Assay
1200 HYPO
1100
1000
Temperatura (°C)
900
800
700
600
500
400
300
200
100
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
79
900
Assay
800
Extrapolación con el método de probabilidad
700
Temperatura (°C)
600
500
400
300
200
100
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Gráfico 3. Extrapolación de la curva TBP de la técnica assay para crudo Maya con
el método de probabilidad en Aspen HYSYS.
80
1600
Spline Cúbico en PRO/II
1400 Probabilidad en HYSYS
1200
Temperatura (°C)
1000
800
600
400
200
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Gráfico 4. Curva TBP generada por el método spline cúbico en PRO/II a partir de
la curva TBP de la técnica assay comparada contra la curva TBP generada a partir
de la extrapolación por el método de probabilidad.
La Gráfica 5 y 6 muestran las curvas de gravedad API para el crudo Istmo y Maya,
el crudo Maya reproduce de manera satisfactoria cada punto de la curva, caso
contrario al crudo Istmo ya que hay desviaciones en los puntos iníciales debido a
que se hizo un ajuste de la gravedad API porque el factor de caracterización (KUOP)
para la fracción de la nafta era superior 11.96 y no podía ser mayor a dicho valor,
así que a partir de la ecuación del factor de caracterización se ajustó la gravedad
API solo en los puntos iniciales que pertenecen a la fracción de la nafta hasta que
estuviera en el límite de 11.96, es el valor que representa a todo el crudo Istmo, el
ajuste se realizó una vez que la torre de destilación hubo convergido y se conociera
la curva de destilación TBP, ya que es necesario conocer la temperatura al 50%
para realizar el ajuste, es la causa de las desviaciones de la Gráfica 5 y el flujo
volumétrico de la fracción de la nafta se modifica por el ajuste de la gravedad API
pero es una desviación mínima.
81
95
90
Assay
85
HYPO
80
75
70
65
60
Gravedad API
55
50
45
40
35
30
25
20
15
10
5
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
85
80 Assay
75 HYPO
70
65
60
55
Gravedad API
50
45
40
35
30
25
20
15
10
5
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Gráfico 6. Curva de gravedad API generada a partir de pseudocomponentes en
Aspen HYSYS contra curva TBP de la técnica assay del crudo 100% Maya.
82
En resumen, se puede realizar la caracterización del crudo Istmo y Maya generando
solo pseudocomponentes en las curvas de la técnica assay ya que es un método
fácil y no generará desviaciones en dichas curvas; la buena reproducción de estas
curvas en el simulador fueron un indicador y el más importante para el diseño de la
torre de destilación que el balance de materia se hizo de forma correcta.
83
4.1.2. Perfil de flujo, presión y temperatura de torre de destilación
atmosférica
Tabla 17. Perfil de flujos, presión y temperatura de la torre de destilación atmosférica cuando se
procesa crudo Istmo y Maya.
Etapa Presión Temperatura Líquido Vapor
Istmo Maya Istmo Maya Istmo Maya
kg/cm2 man. °C MBPD Std.
Condensador 0.707 117 118 81.7 57.4 116.3 56.7
1 1.057 157 152 98.1 65.9 158.1 113.2
2 1.071 166 161 100.2 66.7 168.2 121.5
3 1.085 171 166 99.9 66.5 172.9 122.3
4 1.099 173 169 99.0 65.9 176.0 122.1
5 1.113 176 171 97.9 65.0 178.5 121.4
6 1.127 178 174 96.5 63.9 181.0 120.5
7 1.141 180 176 94.9 62.4 183.6 119.4
8 1.155 183 179 93.0 60.9 186.6 117.9
PA1Retorno 1.169 186 182 328.7 247.7 190.0 116.4
PA1Extracción 1.183 198 194 150.0 96.1 201.2 154.4
11 1.197 208 204 152.0 96.6 209.1 168.9
12 1.212 214 209 150.3 94.4 214.7 169.3
13 1.226 219 214 146.8 91.1 219.7 167.1
14 1.240 224 219 140.3 85.9 225.0 163.8
15 1.254 231 226 126.6 75.8 231.7 158.5
PA2Retorno 1.268 242 238 301.7 227.8 242.7 148.5
PA2Extracción 1.282 264 263 89.0 68.3 264.3 190.1
18 1.296 299 298 77.7 61.2 297.9 186.3
19 1.310 316 316 75.3 58.4 315.3 179.3
20 1.324 324 324 73.8 56.3 322.5 176.4
21 1.338 327 328 71,5 54.0 325.9 174.4
22 1.352 329 330 67.7 51.5 328.1 172.1
23 1.366 331 332 61.9 48.6 330.4 169.7
24 1.380 334 334 45.3 39.3 333.3 165.4
25 1.394 338 337 34.3 33.1 337.3 161.4
26 1.409 343 341 12.7 12.6 342.5 155.3
Alimentación 1.423 349 351 102.9 141.8 347.2 22.7
28 1.437 346 349 98.4 137.8 344.0 15.2
Fondo 1.451 342 346 91.4 130.7 339.6 11.6
84
350
Crudo Istmo
330
Crudo Maya
310
290
Temperatura (°C)
270
250
230
210
190
170
150
130
110
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
Etapa
85
Finalmente, si se extrae de la torre de destilación una gran cantidad de calor
perjudica el perfil de temperatura y la eficiencia de separación (gap’s, Tabla 18) de
los diferentes productos del crudo Istmo y Maya, es una variable importante para el
diseño de la torre de destilación atmosférica y debe ser una cantidad que no
comprometa la eficiencia de separación porque ya se demostró el impacto negativo
que tiene si solo se enfoca el diseño hacía la recuperación de energía.
350
Sin pumparounds
330
310
290
270
Temperatura (°C)
250
230
210
190
170
150
130
110
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
Etapa
86
La distribución del líquido y de vapor a lo largo de la torre también se muestra en la
Tabla 17, la mayor acumulación de líquido en la columna se presenta en los platos
donde regresa el líquido de los pumparounds (Platos 9 y 16), y la mínima cantidad
de líquido se encuentra encima de la alimentación (Plato 26) un valor especificado
para la convergencia en HYSYS de la torre de destilación que es el overflash con
un valor del 5% con base volumétrica con respecto a la alimentación; la gran
cantidad de líquido y vapor en la sección de rectificación cuando se procesa crudo
Istmo es debido a que hay mayor cantidad de componentes ligeros y como
consecuencia requiere mayor cantidad de vapor para remover estos componentes
ligeros de las fracciones más pesadas, contrario cuando se procesa crudo Maya,
hay menor cantidad de componentes ligeros en la sección de rectificación y requiere
menor cantidad de vapor de agotamiento.
El punto de corte (TBP cut point) de cada fracción del crudo Istmo y Maya también
se pueden ver en la Gráfica 9 y 10, solo con colocarse sobre la curva TBP del crudo
Istmo o Maya y ver la temperatura que está en el punto medio de la temperatura
final de ebullición de la fracción ligera y la temperatura inicial de ebullición de la
fracción pesada se encuentra dicho valor, los TBP cut point para cada fracción del
crudo Istmo y Maya son valores muy similares debido a que estuvieron en función
de la temperatura final de ebullición de cada fracción que es la misma especificación
cuando se procesa crudo Istmo y Maya, que sean similares los TBP cut point de
cada fracción no significa que tuvieran el mismo rendimiento porque el crudo Istmo
y Maya difieren en la cantidad de componentes ligeros.
87
Tabla 19. Punto de corte “TBP cut point” para
las diferentes fracciones del crudo Istmo y Maya.
Punto de Corte (TBP cut point)
°C
Maya Istmo
Nafta 183 180
Turbosina 238 229
Diesel 376 373
GOPA 384 391
88
780
Crudo Istmo
730 Nafta
Turbosina
680 Diesel
GOPA
630 Crudo Reducido
580
530
T (°C ) a 760 mm Hg
480
430
380
330
280
230
180
130
80
30
-20 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Gráfico 9. Traslape de las diferentes fracciones del crudo 100% Istmo sobre la curva TBP.
89
1,797 MMkcal/d Ligeros
1,163 MMkcal/d
Agua Amarga
1 1,594 t/d
Reflujo 82 MBPD Std.
Nafta
60.92 MBPD Std.
Carga Térmica = 791 MMkcal/d
Flujo Volumétrico = 172 MBPD Std. 9
Temperatura de Extracción = 198°C
Temperatura de Retorno = 140°C 1
10 82 t/d de Vapor
Agotador 3
Turbosina Turbosina
Carga Térmica = 1,007 MMkcal/d 22.48 MBPD Std.
Flujo Volumétrico = 113 MBPD Std.
16
Temperatura de Extracción = 264°C
Temperatura de Retorno = 160°C 1 259 t/d de Vapor
17 Agotador
3
Diesel
Diesel
68.16 MBPD Std.
23
1
24 56 t/d de Vapor
Agotador 3
Crudo + Vapor de Arrastre TsH=360 °C GOPA
250 MBPD Std. + 238 t/d
GOPA
27 6.58 MBPD Std.
29
TeH=270 °C* 2,570 MMkcal/h*
Figura 31. Configuración de la torre de destilación atmosférica procesando crudo 100% Istmo.
90
1350
1300 Crudo Maya
1250 Nafta
1200 Turbosina
1150 Diesel
1100 GOPA
1050 Crudo Reducido
1000
950
900
850
T (°C ) a 760 mm Hg
800
750
700
650
600
550
500
450
400
350
300
250
200
150
100
50
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
% Volumen Vaporizado
Gráfico 10. Traslape de las diferentes fracciones de crudo 100% Maya sobre la curva TBP.
91
1,545 MMkcal/d Ligeros
773 MMkcal/d
Agua Amarga
1 1,382 t/d
Reflujo 57 MBPD Std.
Nafta
50.69 MBPD Std.
Carga Térmica = 585 MMkcal/d
Flujo Volumétrico = 144 MBPD Std. 9
Temperatura de Extracción = 194°C
Temperatura de Retorno = 140°C 1
10 84 t/d de Vapor
Agotador 3
Turbosina Turbosina
Carga Térmica = 895 MMkcal/d 18.13 MBPD Std.
Flujo Volumétrico = 99 MBPD Std.
16
Temperatura de Extracción = 263°C
Temperatura de Retorno = 155°C 1 173 t/d de Vapor
17 Agotador
3
Diesel
Diesel
46.37 MBPD Std.
23
1
24 35 t/d de Vapor
Agotador 3
Crudo + Vapor de Arrastre TsH=360 °C GOPA
250 MBPD Std. + 238 Ton/día
GOPA
27 3.99 MBPD Std.
29
TeH=270 °C* 2570 MMkcal/h*
Figura 32. Configuración de la torre de destilación atmosférica procesando crudo 100% Maya.
92
4.2. Parámetros de operación en torre de destilación atmosférica
Tabla 20. Resultados comparativos de los parámetros de operación para crudo Istmo y Maya.
Crudo Crudo
Unidades
100% Istmo 100% Maya
Q 1ra. Etapa de condensación MMkcal/d 1,063 773
Q 2da. Etapa de condensación MMkcal/d 1,797 1,545
Temperatura de tanque acumulador °C 38 38
Presión de tanque acumulador kg/cm2 man. 0.707 0.707
Reflujo MBPD Std. 82 57
Temperatura del domo °C 157 152
Presión del domo kg/cm2 man. 1.057 1.057
Q PA1 (Nafta-Turbosina) MMkcal/d 791 585
Flujo PA1 MBPD 172 144
Temperatura de extracción PA1 °C 198 194
Temperatura de retorno PA1 °C 140 140
ΔT PA1 °C 58 54
Q PA2 (Turbosina-Diesel) MMkcal/d 1,007 895
Flujo PA2 MBPD Std. 113 99
Temperatura de extracción PA2 °C 264 263
Temperatura de retorno PA2 °C 160 155
ΔT PA2 entre TRenorto - TPlato de Retorno °C 82 83
MBPD Std. 12.67 12.60
Overflash
% Vol. 5.07 5.04
Temperatura de la zona flash °C 349 351
Presión de la zona flash kg/cm2 man. 1.423 1.423
Temperatura a la salida del horno °C 360 360
Temperatura de alimentación °C 354.4 355.5
Vaporización a la entrada de la torre % mol. 81 76
Temperatura de fondo °C 342 346
Presión de fondo kg/cm2 man. 1.451 1.451
Vapor sección Turbosina t/d 82 84
Vapor sección Diesel t/d 259 173
Vapor sección GOPA t/d 56 35
Vapor sección Agotamiento t/d 972 856
Vapor Total t/d 1,369 1,148
Relación de Vapor, sección de Turbosina lb/BBL 8 10
Relación de Vapor, sección de Diesel lb/BBL 8 8
Relación de Vapor, sección de GOPA lb/BBL 19 19
Relación de Vapor , sección de agotamiento lb/BBL 23 15
Factor de inundación % 76 58
Calor extraído de la columna por PA´s % 39 39
93
La Tabla 20 muestra los parámetros de operación cuando se procesa crudo Istmo
y Maya, se discuten los resultados más relevantes de la torre de destilación.
La temperatura a la salida del horno pudo haber tenido valores entre 343 a 357 °C[2];
las Gráfica 11 y 12 muestran el efecto que tiene sobre el consumo de vapor de agua
en la sección de agotamiento y la vaporización a la entrada de la torre de destilación
cuando se modifica la temperatura a la salida del horno cuando se procesa crudo
Istmo y Maya, con una diferencia de 5% en mol en la vaporización entre 347 a 377
°C a la salida del horno para ambos crudos, a 347 °C hay buena vaporización pero
se requiere una gran cantidad de vapor en la sección de agotamiento para evitar
que la torre de destilación se enfrié a causa del calor extraído en las extracciones
laterales “Pumparounds”, si se quisiera operar a 347 °C a la salida del horno cuando
se procese crudo Maya se tendría que alimentar a la sección de agotamiento 8,500
t/d de vapor de agua, que es una cantidad excesiva de vapor.
Una temperatura de 377 °C a la salida del horno tiene la ventaja de que se puede
recuperar una gran cantidad de componentes de alto valor por la alta vaporización
a la entrada de la torre de destilación comparado con una temperatura de 347 °C a
la salida del horno además requiere menor consumo de vapor pero con la
desventaja de que habrá problemas de craqueo por la alta temperatura, por lo tanto
para evitar un alto consumo de vapor y posibles problemas de craqueo en la torre
de destilación se encontró que es conveniente operar a una temperatura a la salida
del horno de 360 °C, una temperatura que está en equilibrio con estos factores.
94
2000 84
Vaporización a la entrada de la
Vapor de Agotamiento (t/d) Vaporización
1600
1400 82
columna (% mol)
1200
81
1000
80
800
600 79
400
78
200
0 77
345 350 355 360 365 370 375 380
Temperatura a la salida del horno (°C)
Gráfico 11. Efecto sobre la vaporización a la entrada de la torre y vapor en la
sección de agotamiento variando la temperatura a la salida del horno para crudo
100% Istmo.
9000 79
8500 Vapor de Agotamiento
8000 Vaporización
7500
Vapor de Agotamiento (t/d)
7000 78
Vaporización a la entrada
6500
de la columna (% mol)
6000
5500 77
5000
4500
4000
3500 76
3000
2500
2000 75
1500
1000
500
0 74
345 350 355 360 365 370 375 380
Temperatura a la salida del horno (°C)
Gráfico 12. Efecto sobre la vaporización a la entrada de la torre y vapor en la
sección de agotamiento variando la temperatura a la salida del horno para crudo
100% Maya.
95
4.2.2. Extracciones de calor “Pumparounds”
96
350
330
310
290
Temperatura (°C)
270
250
230
210
190
170
150 PA1=Nafta-Turbosina, PA2=Turbosina-Diesel
130 PA1=Turbosina-Diesel, PA2=Diesel-GOPA
110
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
Etapa
97
En la Tabla 22 se muestra el efecto en la eficiencia de separación (gap’s), el gap
con mayor repercusión, como era de esperarse, es en la sección de diesel y GOPA,
el gap de estas fracciones esta fuera del criterio de diseño según lo reportado en la
literatura[2] y [3].
4.2.3. Overflash
98
92 16
Crudo Reducido
91
GOPA 14
Crudo Reducido (MBPD Std.)
90
89
12
88
87 10
86
8
85
84
6
83
82 4
0% 1% 2% 3% 4% 5% 6%
Overflash (% Vol.)
99
conclusión, para evitar este problema es conveniente mayor porcentaje de overflash
aunque no haya un gran rendimiento de GOPA, pero cabe mencionar que lo que no
se recuperó de GOPA en la torre de destilación atmosférica se recuperará en la
destilación al vacío donde se obtiene un gasóleo ligero (GOLV) y uno pesado
(GOPV), además estos productos de la sección de vacío se mezclan junto con el
GOPA y van directamente a un reactor FCC con el fin de obtener gasolina.
100
4.2.5. Hidráulica de torre de destilación atmosférica
Las propiedades físicas de las fracciones más pesadas como el GOPA y el crudo
reducido no son valores totalmente confiables porque no reproducen con precisión
la densidad, peso específico, etc.; esto podría deberse a que las correlaciones del
simulador no son muy buenas a altas temperaturas [24], además de que no hay una
cantidad suficiente de pseudocomponentes en la sección final de la curva de
destilación TBP para ambos crudos, también puede verificarse que el balance de
materia en dichas tablas se cumple para ambos casos, por lo tanto se cumple el
balance de energía.
101
Tabla 24. Resultados comparativos de los diferentes productos de la torre de destilación
atmosférica para crudo Istmo y Maya (Primera parte).
Unidad Crudo 100% Istmo Crudo 100% Maya
Nafta
Flujo MBPD Std. 60.92 50.69
Rendimiento % Vol. 24.37 20.28
Densidad °API 61.45 58.53
Peso específico 0.7333 0.7446
Factor de caracterización 11.94 11.79
Contenido de azufre % peso 0.04 0.20
TBP a 760 mmHg
5% °C 12 24
10% °C 36 38
90% °C 171 171
95% °C 180 183
100% °C 191 193
ASTM D86 a 760 mmHg
5% °C 46 57
10% °C 63 65
90% °C 160 159
95% °C 167 168
100% °C 175 175
Turbosina
Flujo MBPD Std. 22.48 18.13
Rendimiento % Vol. 8.99 7.25
Densidad °API 45.83 43.17
Peso específico 0.7980 0.8101
Factor de caracterización 11.90 11.73
Contenido de azufre % peso 0.23 0.72
Temperatura de niebla °C -46 -47
Temperatura de inflamación °C 64 64
TBP a 760 mmHg
5% °C 166 165
10% °C 179 180
90% °C 229 230
95% °C 235 237
100% °C 250 252
ASTM D86 a 760 mmHg
5% °C 184 184
10% °C 190 191
90% °C 217 217
95% °C 221 222
100% °C 230 230
gap Nafta-Turbosina °C 17 16
102
Tabla 25. Resultados comparativos de los diferentes productos de la torre de destilación
atmosférica para crudo Istmo y Maya (Segunda parte).
Unidad Crudo 100% Istmo Crudo 100% Maya
Diesel
Flujo MBPD Std. 68.16 46.37
Rendimiento % Vol. 27.26 18.55
Densidad °API 34.78 32.28
Peso específico 0.8510 0.8640
Factor de caracterización 11.81 11.67
Contenido de azufre % peso 0.95 1.95
Índice de cetano 53 50
Temperatura de niebla °C -7 -6
Temperatura de inflamación °C 101 105
TBP a 760 mmHg
5% °C 226 226
10% °C 237 244
90% °C 362 364
95% °C 377 378
100% °C 397 393
ASTM D86 a 760 mmHg
5% °C 246 249
10% °C 253 259
90% °C 345 345
95% °C 352 352
100% °C 363 361
gap Turbosina-Diesel °C 26 27
Gasóleo atmosférico pesado (GOPA)
Flujo MBPD Std. 6.58 3.99
Rendimiento % Vol. 2.63 1.59
Densidad °API 25.87 24.09
Peso específico 0.8992 0.9095
Factor de caracterización 11.80 11.68
Contenido de azufre % peso 1.65 2.70
TBP a 760 mmHg
5% °C 337 342
10% °C 356 360
90% °C 439 428
95% °C 450 438
100% °C 474 461
ASTM D86 a 760 mmHg
5% °C 352 355
10% °C 362 365
90% °C 415 405
95% °C 420 409
100% °C 430 419
103
Tabla 26. Resultados comparativos de los diferentes productos de la torre de destilación
atmosférica para crudo Istmo y Maya (Tercera parte).
Unidad Crudo 100% Istmo Crudo 100% Maya
D1160 a 760 mmHg
5% °C 348 352
10% °C 363 367
90% °C 439 428
95% °C 450 438
100% °C 474 461
gap Diesel-GOPA °C -1 2
Crudo Reducido
Flujo MBPD Std. 91.45 130.75
Rendimiento % Vol. 36.58 52.30
Densidad °API 17.01 13.92
Peso específico 0.9528 0.9730
Factor de caracterización 11.82 11.96
Contenido de azufre % peso 2.28 3.62
TBP a 760 mmHg
5% °C 382 396
10% °C 409 423
90% °C 672 991
95% °C 393 1092
100% °C 749 1282
ASTM D86 a 760 mmHg
5% °C 396 395
10% °C 417 418
90% °C 659 1106
95% °C 666 1126
100% °C 684 1164
D1160 a 760 mmHg
5% °C 391 397
10% °C 414 422
90% °C 673 991
95% °C 693 1092
100% °C 749 1282
104
4.4. Criterios para el diseño de la torre de destilación atmosférica
Finalmente para procesar 250,000 barriles por día de crudo Istmo y Maya deben
seguirse los siguientes criterios de diseño que están dentro los valores según lo
reportado en la literatura[2] y [3]:
Tabla 27. Criterios para el diseño de la torre de destilación atmosférica para procesar
250,000 barriles por día de crudo Istmo y Maya.
Unidad
Proceso de destilación Tipo A
Número de platos reales 34
Sistema de condensación Dos Etapas
Presión del tanque acumulador kg/cm2 man. 0.707
Caída de presión por plato kg/cm2 0.014
Overflash % Vol. 5
Temperatura a la salida del horno °C 360
Relación de vapor, sección de turbosina lb/BBL 8 a 10
Relación de vapor, sección de diesel lb/BBL 8
Relación de vapor, sección de GOPA lb/BBL 19
Relación de vapor , sección de agotamiento lb/BBL 15 a 23
Número de extracciones de calor 2
Cantidad de extracción de calor de la columna % 39
Temperatura de retorno de las extracciones de calor °C ≥140
105
Capítulo V
Conclusiones
106
Los criterios de diseño seleccionados para la torre de destilación atmosférica para
procesar crudo Istmo y Maya son los siguientes: La temperatura a la salida del horno
es de 360 °C, el vapor de agotamiento para los productos (turbosina, diesel y
gasóleo pesado atmosférico) y sección de agotamiento se establece en el rango de
8 a 23 lb/BBL. Los traslapes, (5-95) gap, en grados Celsius para crudo Istmo y Maya
son: Nafta-Turbosina, +16 y +17, Turbosina-Diesel, +26 y +27 y Diesel-Gasóleo
Pesado Atmosférico, -1 a +2. El overflash es de 5% con respecto a la carga en la
sección atmosférica.
El objetivo principal se cumplió al lograr establecer los criterios de diseño con base
a lo reportado en Petroleum Refinery Distillation de Watkins[2] y Handbook of
Petroleum Processing de David Jones y Peter Pujadó[3] para llevar a cabo el diseño
de la torre de destilación atmosférica que procesará 250,000 barriles por día de
crudo Istmo y Maya en el simulador de procesos Aspen HYSYS, así mismo se
cumplieron las especificaciones de los productos tal como se requieren en México
hoy en día.
107
agotadores secundarios y la sección del fondo, solo ajustando adecuadamente
estas variables se cumplieron las especificaciones de los productos para el crudo
Istmo y Maya habiendo un equilibrio entre la cantidad de calor que se puede extraer
de la torre de destilación mediante las extracciones laterales y la eficiencia de
separación (gap’s).
108
Bibliografía
109
[1]. David. A., Paul E. Davis, “Revisión de las Operaciones de Refinación: Una
Mirada detrás del cerco”, Houston Texas, 2007.
[3]. Jones S., David J. y Pujadó, Peter R., “Handbook of Petroleum Processing”,
Springer, 2006.
[10]. Instituto Mexicano del Petróleo, “Planta Combinada U-10000 Refinería Gral.
Lázaro Cárdenas, Minatitlán, Veracruz”, (Paquete de Ingeniería Básica, Agosto
2003).
110
[11]. Romolo, Paolo, “Atmospheric Distillation of Crude Oil”, Volume II, Refining
and Petrochemicals, 2010.
[14]. http://procesosbio.wikispaces.com/
[15]. Gerald Kaes, “Refinery Process Modeling”, Athens Priting Company., 2000.
[16]. Fuentes A., Erica y Moreno C., Argenis, “Estudio de las Propiedades
Coligativas Aplicadas a una Torre Despropanizadora” (Tesis de Licenciatura en
Ingeniería Química, Universidad de Oriente Núcleo de Anzoátegui, 2011).
[19]. Espinoza N. Ruth de las N y Ricci Z., Rita, “Estudio y Evaluación del Proceso
de Alquilación de la Refinería Palito”, (Tesis de Licenciatura en Ingeniería
Química, Universidad Central de Venezuela, 2002).
111
[21]. Fogler, H. “Elementos de Ingeniería de las Reacciones Químicas”, Cuarta
Edición, Pearson, 2008.
112
Anexos
113
Anexo A
Índice de Tablas
114
Tabla 1. Propiedades físicas del crudo 100% Istmo…………………………………21
Tabla 3. Rangos de corte en torre de destilación atmosférica para cada fracción del
crudo…………………………………………………………………………………..24
115
Tabla 14. Especificación de las extracciones de calor para crudo Istmo y Maya en
Aspen HYSYS………………………………………………………………….…………74
Tabla 19. Punto de corte “TBP cut point” para las diferentes fracciones del crudo
Istmo y Maya……………………………………………………………………………..88
Tabla 23. Temperatura final de ebullición en ASTM D86 del GOPA modificando el
overflash en torre de destilación atmosférica…………………………………………99
116
Tabla 25. Resultados comparativos de los diferentes productos de la torre de
destilación atmosférica para crudo Istmo y Maya (Segunda parte)………………103
117
Anexo B
Índice de Figuras
118
Figura 1. Curva de destilación TBP que relaciona la temperatura y la cantidad de
volumen vaporizado a presión atmosférica para crudo 100% Istmo……………….22
Figura 6. Ejemplo del típico del desglose de los pseudocomponentes sobre la curva
de destilación TBP para un crudo………………………………………………32
Figura 10. Eficiencia de fraccionamiento a través del gap y overlap entre la fracción
de turbosina y diesel sobre la curva de destilación ATMS D86…………..38
119
Figura 12. Curva de Packie que relaciona la dificultad de separación y la capacidad
de separación del fraccionamiento entre los productos laterales
adyacentes………………………………………………………………………………..41
Figura 21. Curva de destilación TBP de crudo 100% Istmo elaborada con la técnica
assay………………………………………………………………………………….…..63
Figura 22. Curva de gravedad API y contenido de azufre de crudo 100% Istmo
elaboradas con la técnica assay……………………………………….……………….64
Figura 23. Curva de destilación TBP de crudo 100% Maya elaborada con la técnica
assay…………………………………………………………………………….………..65
120
Figura 24. Curva de gravedad API y contenido de azufre de crudo 100% Maya
elaboradas con la técnica assay……………………………………….……………….66
121
Anexo C
Índice de Gráficos
122
Gráfico 1. Curva TBP generada a partir de pseudocomponentes en Aspen HYSYS
comparada con curva TBP de la técnica assay del crudo Istmo…………………...79
Gráfico 3. Extrapolación de la curva TBP de la técnica assay para crudo Maya con
el método de probabilidad en Aspen HYSYS…………………………………………80
Gráfico 4. Curva TBP generada por el método spline cúbico en PRO/II a partir de
la curva TBP de la técnica assay comparada contra la curva TBP generada a partir
de la extrapolación por el método de probabilidad……………………………..……..81
Gráfico 9. Traslape de las diferentes fracciones del crudo 100% Istmo sobre la
curva de destilación TBP………………………………………………………………..89
Gráfico 10. Traslape de las diferentes fracciones de crudo 100% Maya sobre la
curva de destilación TBP………………………………………………………………..91
123
Gráfico 11. Efecto sobre la vaporización a la entrada de la torre y vapor en la
sección de agotamiento variando la temperatura a la salida del horno para crudo
100% Istmo……………………………………………………………………………….95
124