Umss.. Ingenieria Quimica
Umss.. Ingenieria Quimica
Umss.. Ingenieria Quimica
Cochabamba Bolivia
1
INDICE2
1 INTRODUCCION...3
2 OBJETIVOS...3
3 ESTUDIO DE LA DEMANDA Y CAPACIDAD DE PLANTA..4
1 Introduccin.4
2 Dimetil ter y sus usos.4
3 Gas Licuado de Petrleo (GLP)..5
4 GLP en Bolivia.7
5 Comercializacin de DME13
6 Determinacin de la capacidad y localizacin de la planta de DME...14
4 PROCESO DE OBTENCION DE ETER DIMETILICO...14
1 Introduccin.14
2 Produccin de DME por Sntesis Indirecta15
3 Deshidratacin cataltica del metanol..17
4 Propiedades del Dimetil ter..19
5 BALANCES DE MATERIA Y ENERGIA..20
1 Descripcin general del proceso..21
2
3
4
8 SIMULACIN EN HYSYS...96
9 CONCLUSIONES98
10 BIBLIOGRAFIA...98
11 ANEXOS...99
2
1. INTRODUCCION
La produccin de GLP en los ltimos aos no ha sido suficiente para abastecer el
consumo del mercado interno. El subsidio que el gobierno asigna al precio del GLP
genera prdidas al Tesoro General de la Nacin. Por estas razones se plantea producir
dimetil ter a partir de gas natural para utilizarlo como combustible alternativo al GLP.
Una vez identificada la oportunidad que presenta el mercado para los prximos aos,
se pretende elaborar un estudio de pre factibilidad para una planta de dimetil ter a
partir del metanol proveniente del gas natural boliviano y de esta manera dar una
solucin al problema planteado.
El proyecto que se propone consiste en el diseo de un reactor cataltico en el que
tiene lugar la deshidratacin del metanol para la sntesis del ter dimetlico o DME.
El proceso incluye, adems de la deshidratacin cataltica del metanol, la separacin
de los componentes de la corriente de salida del reactor mediante destilacin (para
purificar el producto y recircular el reactante), siendo esta el proceso de mayor
importancia tras la etapa de reaccin qumica.
2. OBJETIVOS
Objetivo General
Llevar a cabo el anlisis y diseo de una planta de produccin de dimetil ter mediante
la deshidratacin cataltica de metanol comercial.
Objetivos Especficos
Produccin
Demanda
Produccin
Demanda de GLP
Amrica latina es un gran mercado para el GLP. Alrededor del 80 % del GLP es para
uso residencial. Mxico y Brasil representan 60 % del consumo de GLP. El ao 2007 la
7
demanda de GLP fue de 28,3 millones de toneladas. Se estima que la demanda para el
2010 ser de 30 millones de toneladas.
3.4.
GLP en Bolivia
En las Tablas 3.3. y 3.4. se muestra la capacidad de proceso de las refineras y las
plantas de separacin que producen GLP en Bolivia.
Tabla 3.3. Produccin de GLP en Refineras [Ministerio de Hidrocarburos y Energa,
2008]
Refinera
Gualberto
Villarroel
Guillermo Elder
Bell
Total
Capacidad
Produccin de
GLP (BBLD)
Instalada
(BBLD)
39.750
Proceso
(BBLD)
25.300
21.500
18.100
623
61.250
43.400
2.560
1.937
Rendimiento
GLP (%)
Rio Grande
Vuelta Grande
Carrasco
Kanata
Paloma
Colpa
Total
91
91
95
95
94
29
Capacidad
Proceso
(MMpcd)
180
100
70
50
40
50
490
Capacidad
produccin GLP
(TM/Da)
267
181
132
101
99
41
821
La Tabla 3.5. muestra el GLP de plantas que es producido por las diferentes
operadoras existentes en el pas y el GLP producido en las refineras de Cochabamba
y Santa Cruz.
El GLP producido por plantas de separacin de gas natural es mucho mayor al
producido por las refineras. En el ao 2008 aproximadamente el 74 % de la produccin
de GLP provino de las plantas de separacin mientras que el 26 % restante fue el
aporte de las refineras.
8
Tabla 3.5. Produccin Nacional Total de GLP 2000 2007 [YPFB, 2008]
Ao
2000
2001
2002
2003
2004
2005
2006
2007
2008
GLP Plantas de
Separacin
(TM/Da)
831
966
1.086
916
961
805
713
725
560
GLP de
Refineras
(TM/Da)
126
126
142
167
189
207
205
256
225
Total (TM/Da)
958
1.092
1.228
1.083
1.151
1.011
918
981
785
En la Tabla 3.5. se puede observar que la produccin de GLP hasta el ao 2002 se fue
incrementando hasta alcanzar un mximo de 1228 TM/Da con una tendencia a
disminuir hasta el ao 2008, con un ligero aumento el ao 2007.
Actualmente la produccin de GLP tanto en refinacin, como en plantas de extraccin
est llegando a su lmite de capacidad de proceso.
3.4.2. Anlisis de Demanda
El GLP, es el segundo producto hidrocarburfero de mayor consumo en el mercado
interno boliviano despus del Diesel Oil.
3.4.2.1.
Promedio (TM/Da)
731
745
791
834
893
933
973
997
1.100
Proyeccin de Oferta
GLP de
Refineras
(TM/Da)
276
317
336
416
497
497
497
497
497
497
GLP de Plantas
(TM/Da)
Total
(TM/Da)
821
821
821
821
821
821
821
821
821
821
1.097
1.138
1.157
1.237
1.138
1.138
1.138
1.138
1.138
1.138
Proyeccin de la Demanda
Para realizar la proyeccin anual del GLP segn los datos histricos se tom en cuenta
los siguientes aspectos principales:
Demanda para el ao n
Demanda para el ao cero (ao base, es conocido)
i= Tasa de crecimiento
n= Nmero de aos transcurridos entre el ao base (0) y el que se va a calcular
(n).
Utilizando la formula anterior se calcul la demanda de GLP pronosticada para el
periodo 2009 al 2028. En la Tabla 3.8. puede notarse que la demanda de este
combustible tiene una tendencia creciente, aumentando el consumo a nivel nacional
aproximadamente 100% en relacin al ao 2008.
Tabla 3.8. Demanda pronosticada de GLP periodo 2008 -2017
Ao
2009
2010
2011
2012
2013
2014
2015
2016
2017
2018
Promedio
(TM/Da)
1.151
1.204
1.259
1.317
1.377
1.441
1.507
1.576
1.649
1.725
Ao
2019
2020
2021
2022
2023
2024
2025
2026
2027
2028
11
Promedio
(TM/Da)
1.804
1.887
1.974
2.065
2.160
2.259
2.363
2.472
2.585
2.704
3.4.3.3.
Produccin
(TM/Da)
1.151
1.204
1.259
1.317
1.377
1.441
1.507
1.576
1.649
1.725
1.804
1.887
1.974
2.065
2.160
2.259
2.363
2.472
2.585
2.704
Demanda
(TM/Da)
1.097
1.138
1.157
1.237
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
1.138
Dficit
(TM/Da)
-54
-66
-102
-80
-59
-123
-189
-258
-331
-407
-486
-569
-656
-747
-842
-941
-1.045
-1.154
-1.267
-1.386
En la Tabla 3.9. se puede observar que habr un dficit pronosticado de GLP durante el
Periodo 2009-2028.
El dficit de GLP el ao 2008, provoc que el gobierno boliviano emitiera un decreto
para la importacin de GLP, lo que representa prdidas para el TGN por el precio
subsidiado que tiene este combustible.
El ingresar el DME como un sustituto al GLP representa un ahorro energtico y
econmico, por lo tanto es una alternativa interesante introducir el DME al mercado
nacional, y destinar el GLP para la exportacin lo que incrementara los beneficios para
el estado boliviano. El dficit de GLP representa un mercado potencial para
combustibles alternativos como es el caso del DME.
12
Precios ($us/Kg)
0,29
1,40
0,69
1,62
1,10
1,11
Precios (Bs/Kg)
2,25
10,91
5,38
12,62
8,57
8,65
Precio ($us/TM)
433
428
534
507
636
Ao
2008
2009
2010
2011
2012
2013
Precio ($us/TM)
665
696
728
761
796
833
Precio ($us/TM)
175
249
270
314
417
472
550
253
372
442
452
513
Comercializacin de DME
Introduccin
15
Antes de enviar esta corriente a las torres de destilacin para obtener el DME
puro, es necesario reducir la presin hasta 10 bares. El DME producido se
obtiene por cabeza de la primera columna de destilacin. El producto de fondo
de esta primera columna se lleva a una segunda torre donde se separa el
metanol que no ha reaccionado del agua, la cual tras ser enfriada se enva a un
proceso de tratamiento de agua residual para eliminar los restos de compuestos
orgnicos que pudiera contener. El metanol que se obtiene por cabeza de esta
columna se recicla a cabeza del proceso
En el mtodo indirecto el metanol es deshidratado con un slido cido a DME
mediante la siguiente ecuacin:
2C H 3 OH C H 3 OC H 3 + H 2 O
La deshidratacin es una simple reaccin exotrmica de deshidratacin
intermolecular y procesos industriales ya han sido establecidos.
El catalizador para esta reaccin es principalmente un slido acido como una
almina o una zeolita.
Cuando la acidez es fuerte, en el caso de la zeolita, la reaccin procede a baja
temperatura mientras que los productos secundarios se incrementan. Los
productos secundarios causan una baja selectividad de DME.
En la Tabla 4.1. se muestra una comparacin de las industrias que se desarrollan la
sntesis de DME por deshidratacin de metanol de manera comercial.
Alimentacin
Tipo de
reactor
Catalizador
Temperatura
(C)
Presin
(Psig)
Conversin
(%)
SK
Metanol
crudo
Lecho fijo
Almina +
Na40-HZSM5
230-330
145-165
164
145-362
145-290
147-162
80
80
70-80
70-85
80
Se ha determinado que las mejores opciones para producir DME por deshidratacin
de metanol son las desarrolladas por las empresas Toyo Engineering Corporation
17
Para poder realizar el anlisis correcto para el reactor en lo que ser la sntesis es
necesario conocer que tipo de reaccin se trata, porque se escoge esa conversin,
el efecto que tiene la temperatura en la reaccin junto con la presin adems de
saber que tipo de alimentacin se tiene para identificar las reacciones que ocurren
y la secuencia que tienen estas si hay ms de una. Para poder hacer esto hay que
conocer los datos cinticos y de equilibrio, lo que a su vez afecta la decisin de
emplear un proceso en batch o continuo, el catalizador usado particularmente para
el producto deseado y si el proceso necesita o no de calor.
4.3.1. Tipo y secuencia segn la alimentacin
La sntesis indirecta involucra la deshidratacin de metanol para la sntesis de
DME, esta reaccin es una combinacin de las reacciones. Esta sntesis se da
mediante la siguiente ecuacin:
2C H 3 OH C H 3 OC H 3 + H 2 O
Se sabe que esta reaccin es exotrmica y que se busca que la conversin sea del
80 % mol de metanol. Esto es porque se reporta que para metanol puro las
conversiones son mayores al 83% en un rango de temperatura de 250C a 400C.
Se tienen entonces una reaccin principal que es reversible y exotrmica que
puede considerarse est conformada por las siguientes reacciones:
2C H 3 OH k 1 (C H 3)2 O+ H 2 O
Metanol
DME
(C H 3 )2 O+ H 2 O k 2 2C H 3 OH
DME
Metanol
18
[ ]
Ea
Pmetanol
RT
K =[ H O ] [ DME ] / [ C H OH ]
{2708
T }
2
19
20
22
5.2.
A continuacin se presenta la figura 5.1., donde se muestra el proceso de produccin del ter dimetlico, y se
presenta una visin general de las principales etapas que tienen lugar, as como los principales equipos que son
necesarios para la operacin.
23
5.3.
5.3.1. Introduccin
Se va a exponer a continuacin, el balance macroscpico de materia al proceso
completo de produccin del DME, teniendo en cuenta que se deben cumplir las
especificaciones del proyecto, es decir, deben producirse 50.000 toneladas anuales de
DME de una riqueza en peso del 99,5%.
5.3.2. Balance macroscpico de materia del proceso completo
El balance de materia del proceso completo, que se presenta a continuacin, se va a
realizar para una base de clculo de 1000 mol/h de alimentacin fresca al proceso para,
una vez resuelto, redimensionar las corrientes para una corriente de cabeza del proceso
de 50.000 toneladas anuales (especificacin del proyecto). La razn por la cual se va a
realizar el balance en caudales molares y no en caudales msicos reside en una mayor
sencillez en la realizacin de los clculos cuando se emplean los caudales molares por
tratarse ste de un proceso reactivo.
A continuacin se muestra el esquema del balance de materia del proceso completo:
24
25
Valor calculado
1100,00 mol/h
36,00 mol/h
0,93 mol/h
136,93 mol/h
0,93 mol/h
220,00 mol/h
476,00 mol/h
440,93 mol/h
695,50 mol/h
217,46 mol/h
441,43 mol/h
2,54 mol/h
558,57 mol/h
114,76 mol/h
28
29
Una vez tengamos todas las corrientes expresadas en toneladas anuales y todas las composiciones
expresadas como composiciones molares, deberemos multiplicar todas las corrientes por un factor de escala
para que el balance quede resuelto para 50,000 toneladas anuales de producto, especificacin del proyecto.
Dicho factor de escala ser = 50,000/177,22 = 282,1352.
30
5.4.1. Introduccin
Como se ha indicado anteriormente, la etapa de mayor relevancia dentro del
proceso de fabricacin del ter dimetlico es la etapa de reaccin qumica, siendo
sta la etapa en la que se centra el objetivo de este proyecto. Por ello, se presenta
a continuacin el balance macroscpico de materia del reactor, teniendo en cuenta
que se pretende conseguir una produccin anual de 50.000 toneladas de ter
dimetlico de un 99,5% de pureza como cabeza del proceso, tal como se ha
planteado en el balance material del proceso completo.
Por otra parte, en el reactor, que opera de forma adiabtica, va a tener lugar una
reaccin exotrmica que va a producir un incremento en la temperatura en el
interior del mismo y que va a afectar por tanto a la velocidad de reaccin, de modo
que sta y, consecuentemente, el diseo del reactor, van a quedar condicionados
por los efectos trmicos asociados a la reaccin qumica. Por ello, ser necesario
resolver el balance de energa del reactor, para poder as conocer el incremento de
temperatura que se va a producir a lo largo del mismo como consecuencia de la
energa liberada por la reaccin de deshidratacin del metanol.
5.4.2. Balance macroscpico de materia del reactor
El balance macroscpico de materia del reactor ha quedado resuelto con el balance
material del proceso completo, planteado en el captulo anterior. Se ha calculado
dicho balance para una produccin de ter dimetlico conocida de 50.000 toneladas
anuales de una pureza del 99,5% w/w, y una conversin por paso del reactor de un
80%. Para estas especificaciones, el balance macroscpico de materia del reactor
queda resuelto tal como se indica a continuacin:
31
32
C P ( T )=0,001
( a+ bT +c T 2 +d T 3 ) kJ
mol
33
34
36
6. DISEO DE EQUIPOS
6.1. Diseo del Tanques de Almacenamiento
6.1.1. Diseo del Tanque de Almacenamiento 1 (Metanol)
6.1.1.1.
37
M =5707,76
kg
h
3 das 24
=4,11 105 kg
h
da
kg
1h
30 min
=5,064 103 kg
h
60 min
39
=8,1 hp=6,04 kW
3
1000 galUS
3,7854 L
1m
El nivel de lquido, como se puede observar, es el mismo que el dimetro del
depsito (en el caso en estudio son casi iguales la altura y el dimetro, siendo
ligeramente superior la altura del mezclador, por lo que es una suposicin vlida). El
40
41
Donde:
Te1 = 25 C
Te2 = 152,8 C
T
T
T
( es
( s25)+ 0,0
( s25)+0,0317 0,0345 =9017,84 { 0,9973 0,0438 (2525)+0,0027 0,0342 (2525)+0 0,0961(2
0,9675 0
10126
470,14 T s 11753,57=6465,67
T s=38,75 C
6.3. Diseo de los Intercambiadores de calor
6.3.1. Diseo del Intercambiador de calor 1 (Entrada del reactor)
6.3.1.1.
Introduccin
P =10
B
)
T +C
Donde A B, y C son constantes para el DME, MeOH y agua, estas constantes son:
DME
MeOH
Agua
Donde:
PTotal = Presin de operacin de la columna 1, en mmHg.
xi = Fraccin molar del componente i en la corriente de entrada.
Pi* (TB) = Presin parcial del componente i a la temperatura de burbuja, calculado
con la ecuacin de Antoine.
Mediante un mtodo de prueba y error, se obtiene el valor de la temperatura de
burbuja.
PTotal: 12000 mmHg
xDME = 0,001; xMeOH = 0,967; xH2O = 0,032
43
TBurbuja = 158,6 C
Para llevar a cabo el clculo del calor necesario por la corriente fra deber tenerse
en cuenta los calores latentes de vaporizacin del metanol y el agua que son,
respectivamente, 36,13 y 40,60 kJ/mol. Deber conocerse las capacidades
calorficas del metanol y del agua en estado lquido evaluadas entre 25C y la
temperatura de burbuja. Asimismo, ser necesario conocer las capacidades
calorficas del metanol, del agua y del DME evaluadas entre 25C y 158C as como
entre la temperatura de burbuja y 250C. Dichos valores se recogen en las
siguientes tablas:
Tabla 6.2. Capacidades calorficas promedio de los componentes en estado lquido y
gaseoso.
Cp Lquido
Cp Gas
(kJ/molC)
(kJ/molC)
Metanol (38,8-158,6C/158,6-250C)
0,049
0,058
Agua (38,8-158,6C/158,6-250C)
0,035
0,037
DME (38,8-158,6C 158,6-250C)
0,051
0,066
Por lo tanto el calor requerido por la corriente fra ser el siguiente:
Qabsorbido=320774,48
{ {[
mol
kJ
kJ
kJ
0,967 0,049
C ( 158,6 38,8 ) C +36,13
+ 0,058
C ( 250 158,
h
mol
mol
mol
Qabsorbido=15294696,2 kJ /h
= 3655520,12 kcal/h
El calor cedido por el vehculo calefactor (vapor de agua a 20 bar y 350C) debe ser
idntico al calor requerido por la corriente fra. La estimacin del caudal de vapor
necesario para lograr dicho grado de calentamiento puede calcularse de la siguiente
forma:
Qcedido =Qabsorbido =m vapor ( h e hs )
Donde:
mvapor: es el caudal msico requerido de vapor, kg/h.
he: entalpa especfica del vapor en las condiciones de entrada, kJ/kg.
hs: entalpa especfica del vapor en las condiciones de salida, kJ/kg.
Los valores de las entalpas especficas del vapor de agua sobrecalentado en las
condiciones de entrada y de salida se encuentran en Principios Elementales de los
44
Qabsorbido
(he hs )
(15294696,2)kJ /h
(3139 633,1)kJ /kg
45
Para el clculo del coeficiente de pelcula tanto interior (h i) como exterior (he) se
emplea la expresion para flujo laminar:
0,33
0,33
0,33
0,14
hi di
di G
di
Cp
=1,86
k
L
w
{ } { } {} { }
Donde:
hi = coeficiente de pelcula interior.
di = dimetro interior de los tubos.
G = velocidad msica.
k = conductividad trmica.
= viscosidad a la temperatura del fluido.
w = viscosidad a la temperatura de la pared del tubo.
CP = capacidad calorfica.
Para gases se puede aplicar la siguiente simplificacin de la ecuacin:
0,33
0,33
0,33
{ } { } {}
hi di
dG
=1,86 i
k
Cp
d
i
L
di N t
( ){
46
En lo que se refiere a los coeficientes de la pelcula exterior en los tubos, o del lado
de la carcasa, se puede utilizar una frmula similar a la usada para los coeficientes
interiores cuando se utiliza un rea de la superficie transversal del flujo y un dimetro
equivalente definidos como se indica a continuacin:
a) rea de la superficie transversal del flujo:
Se considera como tal, la mxima rea entre dos deflectores consecutivos
multiplicada por la relacin existente entre el espacio que dejan dos deflectores
consecutivos (el paso, menos dos radios del tubo) y el existente entre los centros de
estos tubos:
Pt d t
A=D c B
Pt
{ }
Siendo:
Dc = dimetro interior de la carcasa.
B = distancia entre deflectores. La distancia mnima recomendada es 0,3 m.
Pt = paso entre tubos.
dt = dimetro externo de los tubos.
Tendremos:
Dc = 2,286 m.
B = 0,3 m.
Pt = 1 9/16 plg = 0,040 m.
dt = 1 plg. = 0,032 m.
De modo que: A = 0,137 m2
b) Dimetro equivalente del lado de la carcasa:
Se define considerando que a pesar de los deflectores, el flujo se produce a lo largo
de los tubos, en consecuencia, ser cuatro veces el rea de la seccin transversal
de flujo, segn esta hiptesis, dividido por el permetro mojado y a travs del cual
tiene lugar la transmisin de calor. Considerando el espacio comprendido entre tres
tubos del tresbolillo:
Deq =4
[(
) ( )]
d 2t
4 2
P
3 t
8
[ ]
dt
2
=0,36
k
w
{ } { } { }
{ } { }
Sustituyendo:
Deq = 0,137 m.
G = 73828,07 kg/hm2.
k = 0,090 kcal/hmK.
= 0,801 kg/hm.
CP = 0,661 kcal/kgK.
Resulta: he = 76,30 kcal/hm2K
6.3.1.7.
{ }{
Resistencia
R= trmica +
del flujo
interior
}{ }{
}{
Resistencia
Resistencia
Resistencia
Resistencia
trmica
trmica
trmica
+ de la suciedad + trmica
de la suciedad +
del material
del flujo
sobre la
sobre la
del tubo
exterior
superficie interior
superficie exterior
Es decir, que:
R=
F
F
1
t
1
+ i +
+ e+
A i hi
Ai
kA m A e A e he
( )( )( )
( )
Donde:
Ae, Ai = rea de las superficies exterior e interior de los tubos respectivamente.
Am = rea media logartmica = (Ae Ai)/ln(Ae/Ai)
Fe, Fi = unidad de resistencia de suciedad en la superficie exterior e interior del tubo,
respectivamente.
he, hi = coeficientes de transferencia de calor de los fluidos exterior e interior
respectivamente.
48
Ae U
Ae R
Sustituyendo el valor de R y desprecian do la transferencia de calor por conduccin
en la expresin de U se obtiene:
1
U=
Ae
A F
1
+ e i + F e+
A i hi
Ai
he
{( ) ( )
( )}
U=
Resultando para U la siguiente expresin:
{[(
) ( ) ] [( )
( )}
de 1
d
1
+ e Fi + F e +
d
hi
di
he
49
T ml=
{( T ceT fs )( T cs T fe ) }
ln
T ml=
( T ceT fs )
( T csT fe )
{ ( 623523 )( 423301,8 ) }
ln
( 623523 )
( 423301,8 )
Tml = 105,5K.
Despejando el valor de A:
A=
Q
U T ml
A = 2932,61 m2
6.3.1.9. Clculo de las prdidas de carga.
Para el interior de los tubos, la prdida de carga se calcula mediante la ecuacin:
4 f G2 n L
Ptubos=
di
Siendo:
- f = coeficiente de friccin, funcin del nmero de Reynolds.
(el factor de friccin se calcula mediante la grfica A1)
- G = velocidad de flujo msico.
- P = prdida de carga.
- Di = dimetro interior del tubo.
- L = longitud del tubo.
Considerando los siguientes valores:
Re = 373,60.
f = 0,055.
L = 7,315 m.
di = 0,020 m.
= 455,827 kg/m3.
G = 4993,69 kg/hm2.
Resultando: Ptubos = 4,95 kg/ms2 = 0,005 kPa
50
4 f c G2c (N c +1) D c
D eq
Siendo:
f = coeficiente de friccin, funcin del nmero de Reynolds (0,055).
(el factor de friccin se calcula mediante la grfica A2)
Nc + 1 = L/B = 24.
Nc = nmero de deflectores (23).
fc = factor de rozamiento (0,055).
Dc = dimetro interior de la carcasa (2,286 m).
Deq = dimetro equivalente del lado de la carcasa (0,137 m).
GC = 7219,12 kg/hm2.
Resultando:
Pcarcasa = 96,83 kg/ms2 = 0,0968 kPa.
6.3.2. Diseo del intercambiador de calor 2
6.3.2.1.
Introduccin
x i=
P1 ( T R )
P2 ( T R )
P3 ( T R )
P =10
B
)
T +C
Donde A B, y C son constantes para el DME, MeOH y agua, estas constantes son:
DME
MeOH
Agua
Donde:
x = Fraccin molar del componente i en la primera gota del lquido producido.
yi = Fraccin molar del componente i en la corriente de entrada.
PTotal = Presin de operacin, en mmHg.
Pi* (TR) = Presin parcial del componente i a la temperatura de roco, calculado
con la ecuacin de Antoine.
52
Qcedido =320774,48
{ {[
mol
kJ
kJ
kJ
0,194 0,0555
C ( 331174 ) C +36,13
+ 0,0482
C ( 174100 ) C
h
mol
mol
mol
Qcedido =11426399,3 kJ /h
= 2730975 kcal/h
El calor absorbido por el vehculo refrigerante (agua lquida a 1 bar y 25C) debe ser
idntico al calor requerido por la corriente caliente. La estimacin del caudal de vapor
necesario para lograr dicho grado de calentamiento puede calcularse de la siguiente
forma:
Qabsorbido=Q cedido =m agua ( hshe )
Donde:
mvapor: es el caudal msico requerido de agua, kg/h.
he: entalpa especfica del agua en las condiciones de entrada, kJ/kg.
hs: entalpa especfica del vapor en las condiciones de salida, kJ/kg.
53
magua=
magua=
Qcedido
(hs he )
(11426399,3) kJ /h
(2673,2 104,8)kJ /kg
6.3.2.4.
54
0,8
1 /3
0,14
{ } { } { }
hi di
dG
=0,023 i
k
Cp
k
Donde:
hi = coeficiente de pelcula interior.
di = dimetro interior de los tubos.
G = velocidad msica.
k = conductividad trmica.
= viscosidad a la temperatura del fluido.
w = viscosidad a la temperatura de la pared del tubo.
CP = capacidad calorfica.
Para gases se puede aplicar la siguiente simplificacin de la ecuacin:
0,8
1 /3
{ } { }
hi di
dG
=0,023 i
k
Cp
k
di N t
( ){
En lo que se refiere a los coeficientes de la pelcula exterior en los tubos, o del lado
de la carcasa, se puede utilizar una frmula similar a la usada para los coeficientes
55
{ }
Siendo:
Dc = dimetro interior de la carcasa.
B = distancia entre deflectores. La distancia mnima recomendada es 0,3 m.
Pt = paso entre tubos.
dt = dimetro externo de los tubos.
Tendremos:
Dc = 1,676 m.
B = 0,3 m.
Pt = 1 9/16 plg = 0,040 m.
dt = 1 plg. = 0,032 m.
De modo que: A = 0,101 m2
e) Dimetro equivalente del lado de la carcasa:
Se define considerando que a pesar de los deflectores, el flujo se produce a lo largo
de los tubos, en consecuencia, ser cuatro veces el rea de la seccin transversal
de flujo, segn esta hiptesis, dividido por el permetro mojado y a travs del cual
tiene lugar la transmisin de calor. Considerando el espacio comprendido entre tres
tubos del tresbolillo:
Deq =4
[(
) ( )]
d 2t
4 2
P
3 t
8
[ ]
dt
2
0,55
1/ 3
0,14
{ } { } { }
he Deq
D G
=0,36 eq
k
Cp
k
{ } { }
Sustituyendo:
Deq = 0,137 m.
G = 4448,84 kg/hm2.
k = 0,273 kcal/hmK.
= 1,624 kg/hm.
CP = 1,004 kcal/kgK.
Resulta: he = 109,31 kcal/hm2K
6.3.2.7.
{ }{
Resistencia
R= trmica +
del flujo
interior
}{ }{
}{
Resistencia
Resistencia
Resistencia
Resistencia
trmica
trmica
trmica
+ de la suciedad + trmica
de la suciedad +
del material
del flujo
sobre la
sobre la
del tubo
exterior
superficie interior
superficie exterior
Es decir, que:
R=
F
F
1
t
1
+ i +
+ e+
A i hi
Ai
kA m A e A e he
( )( )( )
( )
Donde:
Ae, Ai = rea de las superficies exterior e interior de los tubos respectivamente.
Am = rea media logartmica = (Ae Ai)/ln(Ae/Ai)
Fe, Fi = unidad de resistencia de suciedad en la superficie exterior e interior del tubo,
respectivamente.
he, hi = coeficientes de transferencia de calor de los fluidos exterior e interior
respectivamente.
k = conductividad trmica del material del tubo.
R = resistencia trmica total desde el fluido interior hasta el exterior.
57
Ae U
Ae R
Sustituyendo el valor de R y desprecian do la transferencia de calor por conduccin
en la expresin de U se obtiene:
1
U=
Ae
A F
1
+ e i + F e+
A i hi
Ai
he
{( ) ( )
( )}
U=
Resultando para U la siguiente expresin:
{[(
) ( ) ] [( )
( )}
De 1
d
1
+ e Fi + F e +
Di hi
di
he
58
T ml=
{( T ceT fs )( T cs T fe ) }
ln
T ml=
( T ceT fs )
( T csT fe )
{ ( 604373 )( 373298 ) }
ln
( 604373 )
( 373298 )
Tml = 138,7K.
Despejando el valor de A:
A=
Q
U T ml
A = 1162,08 m2
6.3.2.9.
59
4 f c G2c (N c +1) D c
D eq
Siendo:
f = coeficiente de friccin, funcin del nmero de Reynolds (0,09).
(el factor de friccin se calcula mediante la grfica A2)
Nc + 1 = L/B = 24.
Nc = nmero de deflectores (23).
fc = factor de rozamiento (0,038).
Dc = dimetro interior de la carcasa (1,676 m).
Deq = dimetro equivalente del lado de la carcasa (0,137 m).
Resultando:
Pcarcasa = 32,26 kg/ms2 = 0,0323 kPa
6.4.
60
k s K CH 3 OH C CH 3 OH
R m=
C C2 H 6 O C H 2 O
K eq
))
donde:
CCH3OH,CH2O and Cf son concentraciones de metanol, agua y sitios libres
respectivamente
CMS y CWS, son las concentraciones de metanol y agua adsorbidos en los sitios
activos.
R es la velocidad de deshidratacin del metanol
ks, es la constante de velocidad
Keq es la constante de equilibrio la cual es una funcin de la temperatura
KCH3OH y KH2O, son los coeficientes de adsorcin de las especies respectivas
Reemplazando los valores de las constantes se tiene:
7
( (
1,554 10 C CH 3 OH
R m=
CC 2H 6OC H2O
4019
6,561 104
3
7
2
exp
+3,707 logT 2,783 10 T + 3,8 10 T
26,64
T
T3
( 1+2 ( K CH 3 OH C CH 3 OH )0,5 + K H2 O C H 2 O )
Donde a su vez las concentraciones son funciones de la presin reducida Y (Y=P/P 0),
de la conversin del metanol xm y de la temperatura T:
CCH 3OH =
CC 2 H 6 O =
C H 2 O=
N CH 3 OH N CH 3 OH , 0 ( 1 xm )
y (1x m)
=
=C CH 3 OH ,0 T 0
V
T
T P0
V 0
( 1+ x m )
T0 P
C 2 H 6 O , 0 ( 1xm )
N C 2 H 6 O N
y (1x m )
=
=C C 2 H 6 O ,0 T 0
T
V
T P0
V 0
( 1+ x m )
T0 P
H 2 O, 0 ( 1x m )
H 2O
N
N
y (1 xm )
=
=C H 2 O, 0 T 0
T
V
T P0
V 0
( 1+ x m )
T0 P
Introduccin
61
T
( H r ) =11,7 + ( 2 a1 a2a3 ) 298
1000
62
T
2000
+ ( 2 c 1c 2c 3 ) 2983
T3
T4
+ ( 2 d 1d 2d 3 ) 2984
3000
4000
La expresin que nos relaciona la cada de presin con la longitud del lecho poroso es
la ecuacin de Ergun:
150 ( 1 )
dP
G
1
=
+ 1,75G
3
dz
0 gC D P
DP
} { } {[
Donde:
G: velocidad msica superficial.
0: densidad de la mezcla gaseosa a la entrada del lecho.
DP: dimetro de las partculas del lecho (catalizador).
: porosidad del lecho cataltico.
: viscosidad de la mezcla gaseosa a travs del lecho.
Expresando la presin reducida en funcin de la masa de catalizador se tiene:
} { } {[
150 ( 1 )
dY
1
G
1
=
+1,75 G
3
d mP P0 P AC 0 gC DP
DP
= 37,28 kJ/molK
~
C PC
= 89,48 kJ/molK
= 0,5
DP = 410-3 m.
0 = 11,63 kg/m3
(promedio entre entrada y salida del reactor) = 1,626 cP = 0,0585 kg/hm.
Aplicando el mtodo anteriormente expuesto, conocidos los valores de todos los
parmetros y constantes y las ecuaciones que deben ser resueltas, obtenemos:
Tabla 6.5. Variacin de la conversin cataltica, temperatura, entalpa de reaccin
y presin total con la altura del lecho cataltico.
mP (kg)
0
10
20
30
40
50
60
70
80
90
xm
0,0000
0,0048
0,0068
0,0150
0,0200
0,0240
0,0290
0,0360
0,0420
0,0480
P (kPa)
1600,00
1599,44
1599,05
1598,73
1598,55
1598,01
1597,66
1597,25
1597,11
1596,75
64
T (K)
523,00
523,24
524,33
526,12
530,03
531,17
532,60
534,02
534,44
536,14
100
110
120
130
140
150
160
170
180
190
200
210
220
230
240
250
260
270
280
290
300
310
320
330
340
350
360
370
380
390
400
410
420
430
440
450
460
470
480
490
500
510
520
530
540
0,0560
0,0630
0,0690
0,0760
0,0840
0,0960
0,1030
0,1100
0,1200
0,1310
0,1390
0,1510
0,1620
0,1750
0,1890
0,2030
0,2170
0,2370
0,2510
0,2720
0,2930
0,3100
0,3310
0,3550
0,3720
0,4000
0,4210
0,4380
0,4550
0,4760
0,4970
0,5140
0,5280
0,5450
0,5620
0,5760
0,5890
0,6000
0,6170
0,6280
0,6400
0,6520
0,6620
0,6720
0,6830
1596,44
1596,21
1596,07
1595,79
1595,51
1595,22
1595,05
1594,67
1594,56
1594,23
1594,16
1594,03
1593,85
1593,67
1593,52
1593,33
1593,04
1592,95
1592,82
1592,66
1592,45
1592,33
1592,27
1592,16
1592,04
1591,97
1591,77
1591,56
1591,48
1591,35
1591,25
1591,10
1591,02
1590,95
1590,88
159076
1590,62
1590,55
1590,44
1590,37
1590,27
1590,19
1590,07
1590,00
1589,95
65
537,27
538,69
540,10
541,51
542,91
544,32
545,72
547,12
548,51
549,91
551,30
552,68
554,07
555,45
556,84
558,21
559,59
560,97
562,34
563,71
565,07
566,44
567,80
569,16
570,52
571,87
573,22
574,57
575,92
577,27
578,61
579,95
581,29
582,63
583,96
585,29
586,62
587,95
589,26
590,60
592,05
594,34
595,55
595,87
598,18
550
560
570
580
590
600
610
620
630
640
650
660
670
680
690
700
710
720
730
740
750
760
770
780
790
800
0,6890
0,7000
0,7050
0,7140
0,7240
0,7280
0,7340
0,7390
0,7480
0,7520
0,7590
0,7630
0,7690
0,7750
0,7795
0,7840
0,7870
0,7910
0,7930
0,7960
0,7970
0,7980
0,7985
0,7990
0,7995
0,8000
1589,92
1589,84
1589,75
1589,67
1589,59
1589,50
1589,42
1589,34
1589,26
1589,18
1589,10
1589,02
1588,95
1588,87
1588,79
1588,71
1588,63
1588,55
1588,47
1588,40
1588,32
1588,24
1588,15
1588,07
1587,99
1587,90
598,49
600,80
602,10
604,40
606,70
608,00
609,30
607,59
608,89
610,17
611,46
612,75
614,03
615,31
616,59
617,87
619,14
620,42
621,69
622,96
624,22
625,49
626,75
628,01
629,27
630,52
Tabla 6.6. Variacin de las composiciones molares de los tres componentes con la
altura del lecho cataltico.
mP (kg)
0
10
20
30
40
50
60
70
80
90
100
110
120
130
yMetanol
0,9675
0,9627
0,9607
0,9525
0,9475
0,9435
0,9385
0,9315
0,9255
0,9195
0,9115
0,9045
0,8985
0,8915
yAgua
0,0317
0,0341
0,0351
0,0392
0,0417
0,0437
0,0462
0,0497
0,0527
0,0557
0,0597
0,0632
0,0662
0,0697
66
yDME
0,0008
0,0032
0,0042
0,0083
0,0108
0,0128
0,0153
0,0188
0,0218
0,0248
0,0288
0,0323
0,0353
0,0388
140
150
160
170
180
190
200
210
220
230
240
250
260
270
280
290
300
310
320
330
340
350
360
370
380
390
400
410
420
430
440
450
460
470
480
490
500
510
520
530
540
550
560
570
580
0,8835
0,8715
0,8645
0,8575
0,8475
0,8365
0,8285
0,8165
0,8055
0,7925
0,7785
0,7645
0,7505
0,7305
0,7165
0,6955
0,6745
0,6575
0,6365
0,6125
0,5955
0,5675
0,5465
0,5295
0,5125
0,4915
0,4705
0,4535
0,4395
0,4225
0,4055
0,3915
0,3785
0,3675
0,3505
0,3395
0,3275
0,3155
0,3055
0,2955
0,2845
0,2785
0,2675
0,2625
0,2535
0,0737
0,0797
0,0832
0,0867
0,0917
0,0972
0,1012
0,1072
0,1127
0,1192
0,1262
0,1332
0,1402
0,1502
0,1572
0,1677
0,1782
0,1867
0,1972
0,2092
0,2177
0,2317
0,2422
0,2507
0,2592
0,2697
0,2802
0,2887
0,2957
0,3042
0,3127
0,3197
0,3262
0,3317
0,3402
0,3457
0,3517
0,3577
0,3627
0,3677
0,3732
0,3762
0,3817
0,3842
0,3887
67
0,0428
0,0488
0,0523
0,0558
0,0608
0,0663
0,0703
0,0763
0,0818
0,0883
0,0953
0,1023
0,1093
0,1193
0,1263
0,1368
0,1473
0,1558
0,1663
0,1783
0,1868
0,2008
0,2113
0,2198
0,2283
0,2388
0,2493
0,2578
0,2648
0,2733
0,2818
0,2888
0,2953
0,3008
0,3093
0,3148
0,3208
0,3268
0,3318
0,3368
0,3423
0,3453
0,3508
0,3533
0,3578
590
600
610
620
630
640
650
660
670
680
690
700
710
720
730
740
750
760
770
780
790
800
0,2435
0,2395
0,2335
0,2285
0,2195
0,2155
0,2085
0,2045
0,1985
0,1925
0,1880
0,1835
0,1805
0,1765
0,1745
0,1715
0,1705
0,1695
0,1690
0,1685
0,1680
0,1675
0,3937
0,3957
0,3987
0,4012
0,4057
0,4077
0,4112
0,4132
0,4162
0,4192
0,42145
0,4237
0,4252
0,4272
0,4282
0,4297
0,4302
0,4307
0,43095
0,4312
0,43145
0,4317
0,3628
0,3648
0,3678
0,3703
0,3748
0,3768
0,3803
0,3823
0,3853
0,3883
0,39055
0,3928
0,3943
0,3963
0,3973
0,3988
0,3993
0,3998
0,40005
0,4003
0,40055
0,4008
Tabla 6.7. Variacin de los caudales molares de los tres componentes con la altura del
lecho cataltico.
mP (kg)
0
10
20
30
40
50
60
70
80
90
100
110
120
130
140
150
FMetanol
9797,16
9748,55
9728,30
9645,26
9594,63
9554,13
9503,50
9432,61
9371,85
9311,10
9230,09
9159,20
9098,44
9027,56
8946,55
8825,04
FAgua
321,00
345,31
355,43
396,95
422,27
442,52
467,83
503,28
533,65
564,03
604,54
639,98
670,36
705,80
746,31
807,06
68
FDME
8,10
32,40
42,53
84,05
109,36
129,62
154,93
190,37
220,75
251,13
291,64
327,08
357,46
392,90
433,40
494,16
160
170
180
190
200
210
220
230
240
250
260
270
280
290
300
310
320
330
340
350
360
370
380
390
400
410
420
430
440
450
460
470
480
490
500
510
520
530
540
550
560
570
580
8754,15
8683,27
8582,01
8470,62
8389,61
8268,09
8156,70
8025,06
7883,29
7741,53
7599,76
7397,23
7255,47
7042,81
6830,16
6658,02
6445,36
6202,33
6030,19
5746,65
5534,00
5361,85
5189,71
4977,06
4764,41
4592,26
4450,49
4278,34
4106,20
3964,43
3832,79
3721,40
3549,25
3437,87
3316,35
3194,84
3093,57
2992,31
2880,92
2820,16
2708,77
2658,14
2567,01
842,50
877,95
928,58
984,27
1024,78
1085,54
1141,23
1207,05
1277,93
1348,82
1419,70
1520,96
1591,85
1698,17
1804,50
1890,57
1996,90
2118,41
2204,49
2346,25
2452,58
2538,65
2624,73
2731,05
2837,38
2923,45
2994,34
3080,41
3166,48
3237,37
3303,19
3358,88
3444,95
3500,65
3561,41
3622,16
3672,79
3723,43
3779,12
3809,50
3865,19
3890,51
3936,08
69
529,60
565,05
615,68
671,37
711,88
772,63
828,33
894,15
965,03
1035,92
1106,80
1208,06
1278,95
1385,27
1491,60
1577,67
1684,00
1805,51
1891,59
2033,35
2139,68
2225,75
2311,83
2418,15
2524,48
2610,55
2681,43
2767,51
2853,58
2924,46
2990,28
3045,98
3132,05
3187,75
3248,50
3309,26
3359,89
3410,52
3466,22
3496,60
3552,29
3577,61
3623,18
590
600
610
620
630
640
650
660
670
680
690
700
710
720
730
740
750
760
770
780
790
800
2465,74
2425,24
2364,48
2313,85
2222,71
2182,21
2111,33
2070,82
2010,06
1949,31
1903,74
1858,17
1827,79
1787,28
1767,03
1736,65
1726,53
1716,40
1711,34
1706,27
1701,21
1696,15
3986,71
4006,96
4037,34
4062,66
4108,22
4128,48
4163,92
4184,17
4214,55
4244,93
4267,71
4290,50
4305,69
4325,94
4336,06
4351,25
4356,32
4361,38
4363,91
4366,44
4368,97
4371,51
3673,81
3694,06
3724,44
3749,75
3795,32
3815,57
3851,02
3871,27
3901,65
3932,03
3954,81
3977,59
3992,78
4013,04
4023,16
4038,35
4043,42
4048,48
4051,01
4053,54
4056,07
4058,61
70
V t =2,06 m
Finalmente obtenemos la longitud del reactor que ser necesario para nuestras
condiciones de operacin, tomando como dato:
- A = 0,4418 m2
Por lo que la longitud del reactor viene a ser:
V
2,06
L= t =
=4,6 m
A 0,4418
6.5.
Fraccin Molar
(% mol)
0,194
0,419
0,387
1,000
Metanol
Agua
Dimetil ter
Total
P= 10400 KPa (7800 mmHg)
T= 100 C
El DME, MeOH y agua son los principales componentes que ingresan a la torre
de destilacin. A continuacin se define los componentes clave:
71
( A
P =10
B
)
T +C
Donde A B, y C son constantes para el DME, MeOH y agua, estas constantes son:
Tabla 6.9. Constantes de Antoine
A
B
C
7,878
1473,11
230,00
7,967
1668,21
228,00
7,400
1120,00
241,96
MeOH
Agua
DME
x i=
y 1 PTOTAL
P (T R)
y 2 PTOTAL
P (T R)
y 3 PTOTAL
P3 ( T R )
Donde:
72
=
y C / xC PC /P P C
Donde:
AB = Volatilidad relativa del componente A respecto al componente B.
PA*(Tprom) = Presin parcial del componente A a la temperatura promedio.
PB*(Tprom) = Presin parcial del componente B a la temperatura promedio.
Tabla 6.10. Volatilidades Relativas Calculadas
Presin Parcial
Volatilidad Relativa
PMeOH*(T) = 6679,72 mmHg
AC = 0,2553
PH2O*(T) = 2255,78 mmHg
BC = 0,0862
PDME*(T) = 26164,19 mmHg CC = 1,0000
6.5.1.5.
73
ln
N min =
x LD x HF
x HD x LF
ln AC
Rmin =
Lmin
D
Donde:
q = Condicin trmica de alimentacin.
Vmin = Flujo de vapor mnimo.
Lmin = Flujo de lquido condensado para ser devuelto a la columna.
74
q = 1,0
XDME,Destilado = 0,991
XMeoH,Fondo = 0,313
= 0,9719
Vmin/D = 1,530
Lmin/D = 0,530
Rmin = 0,530
6.5.1.7.
Y=
NN min
N +1
N=N F , min
( N1Y+Y )
min
abcisas
[N-Nmin]/(N+1)
0<x<0.01
0.01<x<0.9
0.9<x<1
0,252856964
0,522033724
0,159273728
La columna fue designada para operar con un reflujo actual de 50% por encima
del reflujo mnimo. Una vez fijada la relacin de reflujo se procede a estimar el
nmero de platos necesarios con la correlacin de Gilliland.
o N = 15,26
El nmero de platos reales requerido es de 16 incluyendo el rehervidor.
Tpicamente se adiciona un plato extra por cada alimento y por cada
intercambiador de calor lateral. En este caso hay un alimento, luego utilizando
este criterio el nmero de platos reales dentro la torre fue de 17.
6.5.1.8.
Q=223560,80
V = 223560,80 mol/h
vap,m = 36,13 kJ/mol
vap,a = 40,50 kJ/mol
vap,d = 21,51 kJ/mol
mol
kJ
kJ
kJ
0,3127 36,13
+0,6823 40,6
+ 0,005 21,51
h
mol
mol
mol
Q=8742743,09 kJ /h
76
Q=4842360,02kJ /h
77
( A
P =10
6.5.2.2.
B
)
T +C
x i=
y 1 PTOTAL
P (T R)
y 2 PTOTAL
P (T R)
y 3 PTOTAL
P3 ( T R )
Donde:
x = Fraccin molar del componente i en la primera gota del lquido producido.
yi = Fraccin molar del componente i en la corriente de entrada.
PTotal = Presin de operacin de la columna 1, en mmHg.
Pi* (TR) = Presin parcial del componente i a la temperatura de roco, calculado
con la ecuacin de Antoine.
Mediante un mtodo de prueba y error, se obtiene el valor de la temperatura de
roco.
78
=
y B / x B PB /P PB
Donde:
AB = Volatilidad relativa del componente A respecto al componente B.
PA*(Tprom) = Presin parcial del componente A a la temperatura promedio.
PB*(Tprom) = Presin parcial del componente B a la temperatura promedio.
Tabla 6.10. Volatilidades Relativas Calculadas
Presin Parcial
Volatilidad Relativa
PMeOH*(T) = 6679,72 mmHg
AA = 1,0000
PH2O*(T) = 2255,78 mmHg
BA = 0,3767
PDME*(T) = 26164,19 mmHg CA = 3,1561
6.5.2.5.
x LD x HF
x HD x LF
ln AB
xLD = 0,7500
xHD = 0,2432
xHF = 0,7900
xLF = 0,2055
AB = 2,6545
Rmin =
Lmin
D
Donde:
q = Condicin trmica de alimentacin.
Vmin = Flujo de vapor mnimo.
Lmin = Flujo de lquido condensado para ser devuelto a la columna.
D = Flujo total de destilado.
Rmin = Relacin de reflujo mnimo (L/D)min.
= Parmetro ajustable.
EI primer paso fue evaluar el parmetro , mediante un proceso iterativo comn
la ecuacin.
Posteriormente se determin el valor de Vmin con la ecuacin; y el valor de Lmin con
la ecuacin; despus se calcul la relacin de reflujo mnimo con la ecuacin;
finalmente se determin el plato mnimo de alimentacin con la ecuacin.
Estos resultados son:
q = 1,0
80
XMeOH,Destilado = 0,980
XH2O,Fondo = 0,790
= 0,7443
Vmin/D = 2,6928
Lmin/D = 1,6928
Rmin = 1,6928
6.5.2.7.
Y=
NN min
N +1
N F =N F , min
( NN )
min
(L/D)min = 1,6928
L/D = 1,5(L/D)min = 2,5392
Abscisa = 0,2391
81
[N-Nmin]/(N+1) = 0,522033724
Nmin = 2,533
o N = 6,40
El nmero de platos reales requerido es de 7 incluyendo el rehervidor. Tpicamente
se adiciona un plato extra por cada alimento y por cada intercambiador de calor
lateral. En este caso hay un alimento, luego utilizando este criterio el nmero de
platos reales dentro la torre fue de 8.
6.5.2.8.
Q=136731,89
V = 136731,89 mol/h
vap,m = 36,13 kJ/mol
vap,a = 40,50 kJ/mol
vap,d = 21,51 kJ/mol
mol
kJ
kJ
kJ
0,2055 36,13
+0,79 40,6
+ 0,0045 21,51
h
mol
mol
mol
Q=5413968,74 kJ /h
6.5.2.9.
V = 223560,91 mol/h
82
Q=136731,90
mol
kJ
kJ
kJ
0,75 36,13
+0,2432 40,6
+0,0068 21,51
h
mol
mol
mol
Q=5075175,16 kJ /h
6.6.
Para impulsar los fluidos a travs de las lneas de tuberas y equipos que
componen la instalacin, ser necesario emplear los adecuados sistemas de
impulsin de fluidos.
Para impulsar la corriente de alimentacin al proceso, que se encuentra en estado
lquido, se emplearn dos bombas centrfugas situadas en paralelo, de modo que
en caso de avera de una de ellas, la otra bomba pueda realizar su funcin sin
interferir en la continuidad del proceso. Estas bombas se encargarn de compensar
la prdida de energa mecnica e impulsar la corriente de alimentacin hasta el
equipo intercambiador de calor. La corriente gaseosa que sale de dicho cambiador
de calor, ser impulsada hasta el reactor por un compresor situado a la salida del
equipo intercambiador de calor.
A continuacin, se muestran las caractersticas esenciales que debern poseer
dichos equipos para que puedan suministrar la potencia necesaria para alcanzar la
presin requerida y superar las prdidas de carga que se producen en la
instalacin.
6.6.1. Diseo de las Bombas 1 y 2
Las bombas 1 y 2 son las encargadas de impulsar la corriente lquida de
alimentacin desde el tanque de almacenamiento hacia la entrada del mezclador.
Se ha dispuesto de dos bombas en paralelo para que en caso de avera, la
instalacin pueda funcionar con normalidad modificando el paso del flujo de la
bomba 1 a la bomba 2 sin interferir as en el proceso, y de modo que sea posible la
reparacin de la bomba averiada sin que sea necesario parar el proceso.
Aplicando la ecuacin de Bernoulli entre el punto A, situado a la entrada de la
bomba, y un punto B, situado a la entrada al reactor (ya que la potencia de la
bomba debe superar las prdidas de carga que tienen lugar en las lneas situadas
desde sta hasta la entrada al reactor):
2
2
PA
vA
PB
vB
+ z A+
+ H util =
+ zB+
+h L
g
2g
g
2g
( ) ( )
( ) ( )
( NSPH )=
(P1Pv )
+ ( z1 za ) h f
g
Donde:
P1: presin en el punto inicial del tramo de aspiracin; 95992,1 N/m 2.
Pv: presin de vapor en el tramo de aspiracin. Se ha tenido en cuenta la
temperatura del fluido, de modo que Pv = 12300 N/m2.
z1: altura en el punto inicial del tramo de aspiracin; 1 m.
za: altura en la aspiracin de la bomba; 0,1 m.
hf: prdida de carga entre el punto inicial y final del tramo de aspiracin:
0,497 m.
: densidad del lquido; 787,47 kg/m 3.
g: aceleracin de la gravedad; 9,81 m/s2.
Sustituyendo los valores en la expresin anterior, se obtiene que las bombas deben
tener una NPSH de 11,24 m.
6.6.2. Diseo de las Bombas 3 y 4
Las bombas 3 y 4 son las encargadas de impulsar la corriente lquida de reflujo
hacia la entrada del mezclador. Se ha dispuesto de dos bombas en paralelo, como
se ha dicho anteriormente, para que en caso de avera, la instalacin pueda
funcionar con normalidad modificando el paso del flujo de la bomba 3 a la bomba 4,
de modo que sea posible la reparacin sin que sea necesario parar el proceso.
Aplicando la ecuacin de Bernoulli entre el punto A, situado a la entrada de la
bomba, y un punto B, situado a la entrada al mezclador:
2
( ) ( )
( ) ( )
PA
v
P
v
+ z A + A + H util = B + z B + B +h L
g
2g
g
2g
B = 662,50 kg/m3.
vB = 0,75 m/s
zB = 0,10 m.
87
Gerencia
Gerencia
General
General
Secretaria
Ejecutiva
Gerencia de
Aministracin
y Finanzas
Secretaria
Secretaria
Departamento
de
de
Contabilidad
Contabilidad
Asesora
Legal
Asesora Legal
Gerencia
Comercial
Secretaria
Secretaria
Departamento
de
de Recursos
Recursos
Humanos
Humanos
Departamento
de Ventas
Departamento
de Adquisicin
de Materia
Prima
Gerencia de
Mantenimient
o
Gerencia
Tcnica
Departamento
de
Mantenimient
o
Secretaria
Secretaria
Departamento
de Publicidad
Planta
Planta
Laboratorio
Laboratorio
Seguridad
Industrial
Supervisin
Supervisin
Supervisin
Supervisin
Supervisin
Operadores
Operadores
Operadores
Operadores
Operadores
Operadores
Operadores
88
Cargo
Gerente General
Secretaria
Gerente Administrativo
Jefe de Departamento
Legal
Asesor Legal
Jefe de Departamento de
Gerencia
Administrativa Contabilidad
Asesor de Contabilidad
Jefe de Recursos Humanos
Secretaria
Personal de Limpieza
Gerente Comercial
Departamento de Ventas
Gerencia de
Encargado Adquisicin
Comercializacin Materia Prima
Encargado Publicidad
Secretaria
TOTAL
Gerencia General
Total
Personal Sueldo
Requerido $us/Mes $us/Mes
1
3500 4668,18
1
400
533,51
1
3000 4001,30
Total
$us/Ao
56018,20
6402,08
48015,60
1
1
1150
850
1533,83
1133,70
18405,98
13604,42
1
1
1
1
6
1
1
1150
850
1000
400
150
3000
800
1533,83
1133,70
1333,77
533,51
400,13
4001,30
2134,03
18405,98
13604,42
16005,20
6402,08
4801,56
48015,60
25608,32
1
1
1
20
1000
800
400
1333,77 16005,20
1067,01 12804,16
533,51
6402,08
25608,32 307299,84
DESCRIPCION
Programa de
Produccin
1
75
100
100
100
100
10
100
100
100
100
7.3.
Costos de inversin
El capital de inversin ha sido estimado en base a datos de referencia empleados en
el material bibliogrfico y algunos datos extrados de "Plant Design and Economics
for Chemical Engineers".
La estimacin se encuentra en moneda de dlares americanos, ninguna
contingencia se ha agregado por inflacin ni fluctuaciones del tipo de cambio
monetario, los montos estimados son para la planta de DME con capacidad de 140
TM/da.
La Tabla 7.4. muestra los resultados de los clculos efectuados del costo de
inversin total estimado bajo las condiciones supuestas anteriormente:
Tabla 7.4. Inversin Total en MM$us para la Planta de DME
ACTIVO FIJO
MM$us
Maquinaria y equipo
1,9
Transporte de FOB a CIF (10%)
0,19
Total compra de equipo
2,09
Instalacin y montaje (40%)
0,836
Tuberas, vlvulas y accesorios
1,254
Instrumentacin y control (8%)
0,1672
Sistema elctrico (10%)
0,209
Servicios auxiliares (30%)
0,627
Infraestructura y obras civiles
0,1045
Subtotal
5,288
ACTIVO DIFERIDO
Ingeniera y supervisin de la obra
0,627
Gastos de construccin (40%)
0,836
Gastos legales (4%)
0,0836
Pago a contratista (22%)
0,4598
Contingencia (10%)
0,209
Subtotal
2,2154
ACTIVO CIRCULANTE
Capital de trabajo (15%)
0,793
MONTO DE INVERSIN TOTAL
8,296
EI capital de inversin fue calculado en base a los siguientes criterios de Peters et
al, (2002).
El capital de trabajo representa el 15% del activo fijo [Peters et al, 2002].
7.4.
Costos de operacin
En esta seccin se realiza la estimacin del costo de operacin del proyecto, se
plantea que la planta de DME funciona por 24 horas cada da, realizndose procesos
peridicos de mantenimiento durante los cuales se desviarn los flujos hacia equipos
duplicados de reemplazo mientras se lleva a cabo el mantenimiento de los equipos.
Se operar la planta 365 das al ao. Los principales costos de operacin se
detallan a continuacin.
7.4.2. Costo de Materia Prima
Los requerimientos de materia prima de la planta de DME, con una capacidad de
50 Ton/Ao, ya fueron desarrollados en el captulo de Ingeniera de Proyecto y son
de 78,998 Ton/Ao de metanol.
El costo de materia prima abarca la mayor parte de los costos variables.
La Tabla 7.5. muestra un detalle de la cantidad y precio de materia prima para la
planta.
Descripcin
Metanol
Total
Costo
($us/Ao)
0,05
1458,25
12774270
638714
0,03
589,17
5161129,2
154834
793547
Servicio a la deuda
Se consider que las condiciones del prstamo fueron las siguientes:
- Monto: 100 % del monto de inversin total
- Plazo: 10 aos, sin periodo de gracia
- Inters: 5% anual
- Forma de Pago: Cuotas iguales de amortizacin
El monto total de inversin es de 8,4 MM$us.
La cuota anual "A" de reembolso constante est dado por la ecuacin [Sapag, 2007]
A=frcP
Donde:
frc = Factor de recuperacin de capital
P = Monto de prstamo
El factor de recuperacin de capital se calcula con la ecuacin [Sapag, 2007]:
93
i(1+i)n
frc=
(1+ i)n1
Donde:
i = Tasa de inters del prstamo
n = Periodo de pago del prstamo
Por lo tanto:
i = 0,05
n = 10
frc = 0,129
A = 8,4 MM$us/Ao
En la Tabla 7.7. se muestra el detalle del servicio a la deuda calculado para los 10
aos.
Perodo
1
2
3
4
5
6
7
8
9
10
7.6.
Monto (MM$us/Ao)
1,668
34,858
36,526
95
10
75
100
100
100
100
100
100
100
100
100
37500
26,648
50000
34,858
50000
34,862
50000
34,866
50000
34,870
50000
34,874
50000
34,881
50000
34,889
50000
34,897
50000
34,905
24,632
32,843
32,843
32,843
32,843
32,843
32,843
32,843
32,843
32,843
0,794
0,739
0,794
0,739
0,794
0,743
0,794
0,747
0,794
0,751
0,794
0,755
0,794
0,762
0,794
0,770
0,794
0,778
0,794
0,786
0,483
0,483
0,483
0,483
0,483
0,483
0,483
0,483
0,483
0,483
1,668
1,668
1,672
1,676
1,680
1,684
1,691
1,699
1,707
1,715
0,307
0,307
0,311
0,315
0,319
0,323
0,331
0,338
0,346
0,354
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,064
0,032
0,083
0,271
0,271
0,271
0,271
0,271
0,271
0,271
0,271
0,271
0,271
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
0,812
0,100
28,316
36,526
36,534
36,542
36,549
36,557
36,573
36,588
36,604
36,619
0,420
0,387
0,352
0,315
0,276
0,235
0,193
0,148
0,101
0,052
0,282
0,282
0,282
0,282
0,282
0,282
0,282
0,282
0,282
0,282
1,872
1,872
1,872
30,89
39,07
39,04
37,14
37,11
37,07
37,05
37,02
36,99
36,95
96
7.7.
Estudio financiero
Una vez determinados las inversiones y los costos en los que incurre el proyecto, se
realiza la evaluacin financiera para determinar la rentabilidad del proyecto. El
anlisis de corrientes de liquidez demostr que el proyecto cuenta con la liquidez
necesaria para cumplir con sus obligaciones.
Asimismo, se realiz el clculo de indicadores financieros como la tasa interna de
retomo (TR), valor actual neto (VAN), relacin Beneficio/Costo y periodo de
recuperacin de capital, finalmente se realiz el anlisis de sensibilidad del proyecto
en funcin del costo de materia prima y producto final principalmente.
7.7.2. Determinacin de ingresos por venta de productos
La planta de DME tiene una capacidad de produccin de 140 TM/Da (5830 Kg/h)
que representa una produccin anual de 50.000 TM/Ao. Se tiene que el precio de
venta de DME es de 833 $us/TM.
En la Tabla 7.12. se muestra la produccin de la planta y los ingresos anuales
percibidos por la venta de DME.
Descripcin
DME
Total
97
1
31,237
5
31,237
5
10
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
31,558
28,316
0,420
0,668
0,282
1,872
-0,320
39,768
36,526
0,387
0,701
0,282
1,872
1,882
39,776
36,534
0,352
0,736
0,282
1,872
1,874
37,912
36,542
0,315
0,773
0,282
0,000
3,738
37,919
36,549
0,276
0,812
0,282
0,000
3,731
37,927
36,557
0,235
0,852
0,282
0,000
3,723
37,943
36,573
0,193
0,895
0,282
0,000
3,707
37,958
36,588
0,148
0,940
0,282
0,000
3,692
37,974
36,604
0,101
0,987
0,282
0,000
3,676
37,989
36,619
0,052
1,036
0,282
0,000
3,661
0,781
0,000
1,041
0,841
1,041
0,833
1,041
2,697
1,041
2,689
1,041
2,682
1,041
2,666
1,041
2,650
1,041
2,635
1,041
2,619
0,000
0,210
0,208
0,674
0,672
0,670
0,666
0,663
0,659
0,655
-1,101
0,630
0,625
2,023
2,017
2,011
1,999
1,988
1,976
1,965
98
0
8,4
1
31,237
5
10
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
42,443
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
41,65
8,4
31,237
5
0,793
8,296
8,296
30,184
38,865
38,871
39,345
39,351
39,357
39,368
39,380
39,392
39,403
28,316
36,526
36,534
36,542
36,549
36,557
36,573
36,588
36,604
36,619
0,420
0,668
0,781
0,387
0,701
1,251
0,352
0,736
1,249
0,315
0,773
1,716
0,276
0,812
1,714
0,235
0,852
1,712
0,193
0,895
1,708
0,148
0,940
1,704
0,101
0,987
1,700
0,052
1,036
1,696
0,104
1,053
2,785
2,779
2,305
2,299
2,293
2,282
2,270
2,258
3,040
0,104
1,157
3,941
6,720
9,025
11,324
13,618
15,899
18,169
20,428
23,467
99
0
0
1
31,238
31,238
2
42
42
3
42
42
4
42
42
5
42
42
6
42
42
7
42
42
8
42
42
9
42
42
10
42
42
8,296
8,296
30,184
0
38,865
0
38,871
0
39,345
0
39,351
0
39,357
0
39,368
0
39,380
0
39,392
0
39,403
0
28,316
36,526
36,534
36,542
36,549
36,557
36,573
36,588
36,604
36,619
-8,296
0,420
0,668
0,781
1,053
0,387
0,701
1,251
2,785
0,352
0,736
1,249
2,779
0,315
0,773
1,716
2,305
0,276
0,812
1,714
2,299
0,235
0,852
1,712
2,293
0,193
0,895
1,708
2,282
0,148
0,940
1,704
2,270
0,101
0,987
1,700
2,258
0,052
1,036
1,696
2,247
-8,296
-7,243
-4,459
-1,680
0,625
2,924
5,218
7,499
9,769
12,028
14,274
100
VAN =I 0 +
i=1
Fi
(1+t)i
Donde:
I0 = Inversin inicial
Fi = Flujo neto anual
t = Tasa de rentabilidad mnima aceptable,
n = Aos de duracin del proyecto.
El valor actual neto puede resultar una cantidad positiva o negativa. Si el VAN es
positivo significa que la inversin tendr ganancias y el proyecto ser aceptable.
Por el contrario, si el VAN es negativo significa que la inversin tendr prdidas y
el proyecto se rechazar. Alternativamente un VAN de cero indica que la inversin
no producira ni ganancias ni prdidas debido a que el proyecto no agregara valor
monetario, la decisin de aceptar o rechazar el proyecto debera basarse en otros
criterios, como la obtencin de un mejor posicionamiento en el mercado u otros
factores.
El clculo de VAN se realiz tomando los datos del flujo de caja presentados en la
Tabla 7.13.; siendo los siguientes:
- I0 = 8,4 MM$us
- Fi = Flujo neto anual
- t = 10 %
- n = 10 aos
A partir de la realizacin de estos clculos, se establece el valor del VAN que es 9,0
MM$us.
7.7.5. Tasa Interna de Retorno (TIR)
Para el clculo de la TIR del proyecto de DME se aplic la siguiente formula de
interpolacin lineal [Sapag 2007]:
TIR=i1 +(i2 i1 )
VAN 1
VAN 1+VAN 2
101
Donde:
i1 = Tasa de actualizacin del ltimo VAN positivo
i2 = Tasa de actualizacin del primer VAN negativo
VAN1 = Valor actual neto obtenido con i1.
VAN2 = Valor actual neto obtenido con i2.
En cuanto a la conveniencia de realizar Ia inversin la decisin ser favorable
cuando la tasa interna de retorno sea mayor que la tasa de inters en el mercado,
esta tasa alcanza un valor de 10 %.
TIR < i Proyecto econmicamente no rentable
TIR = i Proyecto econmicamente no concluyente
TIR > i Proyecto econmicamente recomendable
El valor del TIR calculado para el proyecto es de 22,94%, como el inters es menor
al valor de la TIR, se determin que el proyecto estudiado es econmicamente
rentable debido a que este valor es mayor a tasa de actualizacin del 10%.
7.7.6. Relacin Beneficio/Costo
La relacin Beneficio/Costo (B/C) muestra la cantidad de dinero actualizado que
recibir el proyecto por cada unidad monetaria invertida. Se determina dividiendo
los ingresos brutos anualizados (beneficios) entre los costos anualizados [Sapag,
2007].
A continuacin se describen las situaciones que se pueden presentar en la relacin
Beneficio/Costo.
Relacin B/C > 1: Indica que los ingresos son mayores a los egresos, en
consecuencia el proyecto es aconsejable o debe aceptarse.
Relacin B/C = 1: Indica que los ingresos son iguales a los egresos, en
consecuencia el proyecto es indiferente. Y la decisin de aceptar o rechazar el
proyecto depende de otros factores.
Relacin B/C < 1: Indica que los ingresos son menores a los egresos, por lo tanto
el proyecto no es aconsejable o debe rechazarse.
La frmula utilizada es la siguiente [Sapag, 2007]:
n
Vi
(1+i)
n
B i=0
= n
C
Ci
(1+i)
n
i=0
Donde:
B/C = Relacin Beneficio/Costo
Vi = Valor de la produccin (beneficio bruto)
Ci = Egresos (i = 0, 2, 3, 4,, n)
i = Tasa de descuento
102
El clculo de la relacin B/C se realiz con los datos del flujo de fondos de la Tabla
7.15. Tomando una tasa de actualizacin del 10 %.
Relacin B/C calculada fue 1,03, como el valor calculado es superior a la unidad
indica que los ingresos del proyecto son mayores a los egresos, por lo tanto el
proyecto es aconsejable o debe adaptarse.
7.7.7. Periodo de Recuperacin de Capital (Payback)
Indica el nmero de aos requerido para recuperar el capital invertido en un
proyecto o el nmero de aos que han de transcurrir para que la acumulacin de
los flujos de previstos iguales a la inversin inicial. Si el payback es menor al
periodo definido por la empresa entonces se acepta el proyecto. En el estudio se
calcul mediante la siguiente frmula [SaPag, 2007]:
N
I 0 = FCt
t=1
Donde:
I0 = Inversin inicial
FC = Flujos de caja
t = Periodo de tiempo 1,2,..., n.
Los resultados se presentan en la Tabla 7.16.
Tabla 7.16. Periodo de Recuperacin de capital en MM$us/Ao
Descripcin
Flujo Neto
Flujo Neto
Acumulado
10
-8,296
1,053
2,785
2,779
2,305
2,299
2,293
2,282
2,270
2,258
-8,296
-7,243
-4,459 -1,680
0,625
2,924
5,218
7,499
9,769
12,028
103
8.
9.
SIMULACIN EN HYSYS
10. Finalmente, con base en los clculos realizados en las secciones previas de balance de masa y de diseo de
equipos se llev a cabo la simulacin del proceso en el programa HYSYS 8.0, obtenindose el diagrama grfico y las
tablas presentadas a continuacin.
11.
12.
13.
104
14.
15.
16.
105
17.
24.
25.
CONCLUSIONES
18.
19.
Concluido el desarrollo del presente trabajo se
concluye que se llev a cabo de la manera ms completa posible el anlisis y diseo
de una planta de produccin de dimetil ter mediante la deshidratacin cataltica de
metanol comercial.
20.
Para ello se condujo una investigacin respecto al
mtodo de produccin de dimetil ter a partir de metanol, as como los parmetros y
condiciones necesarios para conducir la reaccin en cuestin tal que se obtenga el
mejor rendimiento y la mejor calidad de producto posible.
21.
En el proceso se disearon los equipos bsicos
empleados en la produccin de dimetil ter a partir de metanol, clculos que luego nos
permitieron efectuar la simulacin del proceso en el software HYSYS.
22.
Finalmente, se realiz el correspondiente anlisis
financiero del proceso indispensable para evaluar la factibilidad econmica del
proyecto y as estar confirmar la existencia de ganancias y beneficios provenientes de
la inversin econmica en este proyecto.
23.
BIBLIOGRAFIA
Estudio de Prefactibilidad de una planta de dimetil ter a partir de Gas Natural
Boliviano. Castro Alba, Vanesa Ins. Tesis Universidad Mayor de San Simn. 2010.
Diseo de reactor cataltico para el proceso de obtencin de dimetil ter a partir de
metanol. Chica Prez, Carlos. Tesis Universidad de Cdiz. 2008.
G. Bercic, J. Levec, "Catalytic dehydration of methanol to dimethyl ether. Kinetic
investigation and reactor simulation," Industrial & Engineering Chemistry Research,
vol. 32, pp. 2478-2484, 1993.
Liang Zhang, Hai-Tao Zhang, Wei-Yong Ying, Ding-Ye Fang, Intrinsic Kinetics of
Methanol Dehydration over Al2O3 Catalyst, World Academy of Science, Engineering
and Technology vol. 59, pp. 1538-1543, 2011.
Nassir Sapag Chain, Reinaldo Sapag Chain, Preparacin y evaluacin de
proyectos, 1989.
Max S. Peters, Klaus D. Timmerhaus, Plant Design and Economics for Chemical
Engineers, 1991.
26.
106
27.
28.
TEMARIO
29.
1. GENERALIDADES
30.
1.1.INTRODUCCIN
31.
1.2.ANTECEDENTES
32.
33.
1.4.OBJETIVOS
34.
1.5.JUSTIFICACIN
35.
1.6.ALCANCE
36.
37.
2. FUNDAMENTACIN TEORICA
38.
2.1.INGENIERA BASICA
39.
40.
41.
42.
43.
44.
45.
46.
BIBLIOGRAFA
47.
107
ANEXOS
108
Tabla A.1.Densidad del metanol a diferentes concentraciones y temperaturas (Perry et al., 2001)
Tabla B.1. Tubos de 1 de dimetro exterior, con paso triangular de 1 9/16. Fuente:
Normas T.E.M.A. (Standards of Tubular Exchanger Manufacturers Association).
109
110
Nivel macro
montaje de equipos.
La EIBH considera a santa cruz como eje industrial al campo rio grande, por donde
pasan la mayor cantidad de gasoductos y distintas plantas.
111
es el segundo mayor
1.1.1.2.
Nivel micro
625449.59
latitud
181126.40
112
113