Producción de Éter Dimetílico A Partir de Metanol: Universidad Mayor de San Simón Facultad de Ciencias Y Tecnología
Producción de Éter Dimetílico A Partir de Metanol: Universidad Mayor de San Simón Facultad de Ciencias Y Tecnología
Producción de Éter Dimetílico A Partir de Metanol: Universidad Mayor de San Simón Facultad de Ciencias Y Tecnología
Cochabamba Bolivia
1
INDICE2
1. INTRODUCCION...3
2. OBJETIVOS...3
3. ESTUDIO DE LA DEMANDA Y CAPACIDAD DE PLANTA..4
3.1. Introduccin.4
3.2. Dimetil ter y sus usos.4
3.3. Gas Licuado de Petrleo (GLP)..5
3.4. GLP en Bolivia.7
3.5. Comercializacin de DME13
3.6. Determinacin de la capacidad y localizacin de la planta de DME...14
4. PROCESO DE OBTENCION DE ETER DIMETILICO...14
4.1. Introduccin.14
4.2. Produccin de DME por Sntesis Indirecta15
4.3. Deshidratacin cataltica del metanol..17
4.4. Propiedades del Dimetil ter..19
5. BALANCES DE MATERIA Y ENERGIA..20
5.1. Descripcin general del proceso..21
5.2. Esquema del proceso..21
5.3. Balance macroscpico de materia y energa al proceso completo22
5.4. Balances macroscpicos de materia y energa del reactor.29
6. DISEO DE EQUIPOS...34
6.1. Diseo del Tanques de Almacenamiento.34
6.2. Diseo del Mezclador35
6.3. Diseo de los Intercambiadores de calor40
6.4. Diseo del Reactor Cataltico.54
6.5. Diseo de las Columnas de Destilacin..65
6.6. Diseo de las Bombas..75
7. ANALISIS FINANCIERO DEL PROYECTO...79
7.1. Estructura orgnica de la planta.82
7.2. Plan de ejecucin del proyecto...82
7.3. Costos de inversin.83
7.4. Costos de operacin..84
7.5. Servicio a la deuda..84
7.6. Resumen costos de operacin...85
7.7. Estudio financiero86
8. SIMULACIN EN HYSYS...96
9. CONCLUSIONES98
10. BIBLIOGRAFIA...98
11. ANEXOS...99
2
1. INTRODUCCION
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
3
3. ESTUDIO DE LA DEMANDA Y CAPACIDAD DE PLANTA
3.1. Introduccin
4
3.2.1. Produccin y demanda de DME a Nivel Mundial
La produccin mundial total de DME est entre 100.000 y 150.000 toneladas mtricas
por ao; una lista de los pases productores de DME se presenta en la Tabla 3.1.
5
Existen dos formas de obtener GLP, una es en plantas de separacin de gas natural y
la segunda consiste en la produccin del GLP derivado del petrleo, tal como se
obtiene de las refineras del pas.
3.3.1.1. Produccin
3.3.1.2. Demanda
La tasa de crecimiento promedio mundial total para el GLP desde el ao 2000 es del 2
%/ao, sin embargo las tasas de crecimiento en las diferentes regiones del mundo han
variado.
Domstico
Qumica
Industria
Agricultura
Refinera
Transporte
3.3.2.1. Produccin
6
Brasil, Venezuela y Argentina. Juntos, estos cuatro pases representaron ms del 80 %
de la produccin de GLP en Latinoamrica. El procesamiento de gas representa ms
del 60 % de la produccin de GLP, se prev que este porcentaje se incremente en el
futuro.
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
demanda de GLP fue de 28,3 millones de toneladas. Se estima que la demanda para el
2010 ser de 30 millones de toneladas.
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.
7
Total 490 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.
Tabla 3.5. Produccin Nacional Total de GLP 2000 2007 [YPFB, 2008]
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
Ao Promedio (TM/Da)
2000 731
2001 745
8
2002 791
2003 834
2004 893
2005 933
2006 973
2007 997
2008 1.100
Por otro lado para la produccin de GLP de plantas separadoras se estima que estas
trabajarn al mximo de su capacidad de operacin y que probablemente se realizarn
ampliaciones de las mismas. En esta estimacin no se tom en cuenta los proyectos de
instalacin de nuevas plantas de separacin de GLP.
La Tabla 3.7. muestra el resumen de estas proyecciones.
9
3.4.3.2. Proyeccin de la Demanda
Para realizar la proyeccin anual del GLP segn los datos histricos se tom en cuenta
los siguientes aspectos principales:
Ao Promedio Ao Promedio
(TM/Da) (TM/Da)
2009 1.151 2019 1.804
2010 1.204 2020 1.887
2011 1.259 2021 1.974
2012 1.317 2022 2.065
2013 1.377 2023 2.160
2014 1.441 2024 2.259
2015 1.507 2025 2.363
2016 1.576 2026 2.472
2017 1.649 2027 2.585
2018 1.725 2028 2.704
10
3.4.3.3. Demanda Insatisfecha Pronosticada de GLP
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.
11
3.4.4. Precio de GLP
Segn datos histricos internacionales, para los ltimos 5 aos evaluados, el precio de
venta de DME grado combustible oscila entre 430-630 $us/TM como se muestra en la
Tabla 3.11.
Ao Precio ($us/TM)
2004 433
2005 428
2006 534
2007 507
2008 636
( ) ( )
12
Ao Precio ($us/TM)
2008 665
2009 696
2010 728
2011 761
2012 796
2013 833
Segn datos histricos internacionales, para los ltimos aos, el precio de venta de
Metanol grado comercial oscila entre 246-552 $us/TM como se muestra en la Tabla
3.13.
Ao Precio ($us/TM)
2002 175
2003 249
2004 270
2005 314
2006 417
2007 472
2008 550
2009 253
2010 372
2011 442
2012 452
2013 513
13
La similitud de las propiedades de manipulacin y almacenamiento que tiene ambos
combustibles permite realizar una mezcla de DME/GLP hasta 30/70 (30 %, DME; 70 %
GLP, en masa). Debido a la mezcla 30/70 la prdida de valor calorfico es de 7%. Una
prdida del 7 % en el poder calorfico, es lo suficientemente pequea para pasar
desapercibida por el consumidor final. Es notable que, en una base independiente, el
DME tenga un menor poder calorfico que el GLP. En la mezcla el DME mejora el
rendimiento del GLP. Cuando se mezcla DME/GLP, el GLP se quema completamente
evitando prdidas. Esto se debe al alto contenido de oxgeno que tiene el DME.
4.1. Introduccin
14
El proceso de produccin de metanol est tcnicamente probado y es
altamente confiable.
El proceso de deshidratacin de metanol es relativamente simple y est
comercialmente establecido.
El escalamiento del proceso de produccin de DME se puede lograr con poco
riesgo tcnico por la integracin de los dos procesos mencionados
anteriormente.
La tasa de produccin de metanol y DME puede ser ajustado a cualquier
relacin dependiendo de la demanda de mercado.
Los parmetros del proceso (temperatura y presin) son ms flexibles y el
rango de conversin de metanol es ms alto.
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
Se ha determinado que las mejores opciones para producir DME por deshidratacin
de metanol son las desarrolladas por las empresas Toyo Engineering Corporation
16
y Haldor Topsoe debido a su alto grado de conversin y los intervalos de operacin
en los que se efecta la reaccin.
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.
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:
( )
Metanol DME
( )
DME Metanol
17
Bandiera y Naccache (1991), hasta modelos complejos con mecanismos
catalticos de varios para la deshidratacin cataltica del metanol.
Un ejemplo de la velocidad de reaccin de la deshidratacin cataltica del metanol,
emplendose especficamente un catalizador zeoltico, viene expresada por la
ecuacin de Bandiera y Naccache (1991):
[ ]
{ } , - , - , -
18
De modo general se ha encontrado que la reaccin segn el modelo de Bercic y
Levec controlada por la velocidad superficial de reaccin y con adsorcin
disociativa del metanol se correlaciona bien con resultados experimentales.
19
En estado gaseoso su poder calorfico neto alcanza un valor de 28,8 MJ/Kg, que
es ms alto que del metano pero ms bajo que del propano y butano.
Considerando las propiedades de combustin su lmite de explosin es ms
amplio que del propano y butano, pero casi idntico al del metano y muy cercano
al del metanol.
El DME es un compuesto muy estable que solo reacciona o se descompone en
condiciones muy crticas de presin y temperatura. La velocidad de
descomposicin estimada a 493 K es alrededor de uno por ciento por ao y es
inerte qumicamente. Tiene una alta solubilidad en sustancias polares y no polares.
20
5.2. Esquema del proceso
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.
21
5.3. Balance macroscpico de materia y energa al proceso completo
5.3.1. Introduccin
22
Figura 5.2. Balance de materia del proceso completo
23
La reaccin que va a tener lugar en el reactor, y que va a ser el centro de todo
nuestro proceso, va a ser la siguiente:
Por ltimo se debe definir la denominacin que se dar a los tres componentes
presentes para la resolucin del balance de materia. El componente 1 ser el
metanol, componente 2 ser el agua y, finalmente, componente 3 ser el DME.
Todas las corrientes estarn expresadas en mol/h.
24
acumulacin, generacin y consumo tienen valor cero, y dicho balance ser de la
forma: (Entrada) = (Salida).
(Por cada mol de metanol, se consume 0,8 y se produce 0,4 mol de DME y 0,4 mol
de agua; de modo que queda sin reaccionar 0,2 mol de metanol).
1. Metanol: R1 = C1 + D1.
2. Agua: R2 = 0,009 D D1 + 0,995 C C1.
3. DME: R3 = 0,005 C + 0,991 D.
1. Metanol: C1 = B1 + 0,750 A.
2. Agua: 0,995 C C1 = 0,790 B + 0,250 A A3.
3. DME: 0,005 C = 0,210 B B1 + A3.
25
Contabilidad y resolucin del balance
Por tanto este balance va a tener dos grados de libertad, a los cuales daremos
valores coherentes para poder resolver el balance de materia. Se elegir como
grados de libertad las composiciones de entrada al reactor de metanol y agua, es
decir E1 y E2.
La resolucin de este balance se llevar a cabo mediante el empleo de una hoja de
clculo, en la cual se definirn los valores ya conocidos de la base de clculo y
composiciones y se introducirn las ecuaciones de los balances. Dndole valores a
los dos grados de libertad establecidos, resolveremos el balance, y podremos
conocer el valor de todas las variables que intervienen en el balance de materia del
proceso completo. Dndole a los dos grados de libertad los valores de E 1 = 1100
mol/h y E2 = 36 mol/h, se resuelve el balance de materia para dar los siguientes
valores de las incgnitas (se puede comprobar que con los valores obtenidos las
ecuaciones se cumplen perfectamente y todos los balances cuadran):
26
Figura 5.3. Resolucin del balance de materia al proceso completo
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
27
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.
28
5.3.3. Balance macroscpico de energa del proceso completo.
5.4.1. Introduccin
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.
29
Figura 5.5. Balances de materia al reactor
Para poder realizar el balance energtico del reactor, ser necesario conocer la
expresin que relaciona las capacidades calorficas a presin constante de los tres
componentes con la temperatura. Dicha expresin viene dada por la ecuacin:
( )
( )
30
donde la temperatura viene expresada en C, y los valores de a, b, c y d estn
recogidos en la siguiente tabla:
Tabla 5.2. Valores de las capacidades calorficas promedio para los diversos
componentes
Cp Metanol Agua ter dimetlico
(kJ/molC)102
25C-250C 5,2532 3,5940 8,0896
25C-400C 5,7865 3,7277 8,9478
En primer lugar, se plantear el balance energtico del proceso reactivo, que viene
descrito en la figura anterior. En este caso, la corriente de entrada slo estar
constituida por el reactivo metanol, mientras que la corriente de salida estar
constituida por los componentes agua y ter dimetlico, productos de la reaccin,
as como por el metanol que no ha reaccionado.
31
Figura 5.6. Balance macroscpico de energa del reactor. Proceso reactivo.
* +
Hs = 62071,0 mol/h0,057865(T 25) kJ/mol + 124142 mol/h0,037277(T 25)
kJ/mol + 124142 mol/h0,089478(T 25) kJ mol
Hs = 3591,74(T 25) kJ/h + 4627,64(T 25) kJ/h + 11108(T 25) kJ/h
Hs = 19327,36(T 25) kJ/h
Trmino de reaccin:
,( ) -
32
Figura 5.7. Balance macroscpico de energa del reactor. Proceso no reactivo.
He = {FaHa + FdHd}e
He = 10157,07 mol/h 8,0862 kJ/mol + 262,39 kJ/mol 18,2015 kJ/mol
He = 86908 kJ/h
Hs = {FaHa + FdHd}e
Hs = 10157,07 mol/h3,727710-2(T 25) kJ/mol + 262,39 kJ/mol8,947810-2(T
25) kJ/mol
Hs = 402,10(T 25) kJ/h
33
6. DISEO DE EQUIPOS
Dimetro = 7,06 m
Altura = 7,06 m
34
Ahora se calcula el volumen de ter dimetlico a almacenar al poseer el valor de la
densidad del DME, 604 kg/m3. Por lo tanto, el volumen til del depsito se
determina:
Dimetro = 6,61 m
Altura = 6,61 m
Dimetro = 2,05 m
Altura = 2,05 m
35
6.2.2. Caractersticas del mezclador:
36
anchura de las hlices y la anchura de las cuatro placas deflectoras o buffles se
calculan mediante las expresiones que se pueden observar en la Figura 6.1.
Dado que cada una de las corrientes entrantes a la mezcladora se encuentra a una
diferente temperatura hace falta realizar el clculo de la temperatura de salida de la
mezcla. Asumiendo que las prdidas de calor durante el mezclado sean
despreciables se obtiene el balance siguiente:
* ( ) ( ) ( )+
* ( ) ( ) ( )+
* ( ) ( ) ( )+
Donde:
Te1 = 25 C
Te2 = 152,8 C
* ( ) ( )
( )+
* ( ) ( )
( )+
* ( ) ( )
( )+
6.3.1.1. Introduccin
37
desde 25C a 250C y en cuyo interior tendr lugar la vaporizacin del metanol y del
agua. Este intercambio trmico ser llevado a cabo en un cambiador de carcasa y
tubos.
( ) ( ) ( )
Donde A B, y C son constantes para el DME, MeOH y agua, estas constantes son:
38
A B C
DME 7,400 1120,00 241,96
MeOH 7,878 1473,11 230,00
Agua 7,967 1668,21 228,00
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.
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:
2 {[ ( )] [ ( )]}
{[ ( )] [ ( )]} {[ ( )]
[ ( )]}3
{ 20 ( ) 1
0 ( ) 13
39
20 ( ) 1
0 ( ) 13
20 ( ) 1
0 ( ) 13}
= 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:
( )
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
Procesos Qumicos (Richard M. Felder & Ronald W. Rousseau).
( )
( )
( )
( )
40
Las caractersticas del intercambiador utilizado para esta operacin de intercambio
de calor son las que se muestran a continuacin:
Para el clculo del coeficiente de pelcula tanto interior (hi) como exterior (he) se
emplea la expresion para flujo laminar:
{ } { } { } { }
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.
{ } { } { }
41
Llamando nt al nmero de pasos, por el lado de los tubos, L a la longitud de
los tubos, Nt al nmero total de tubos del haz y M al caudal msico del fluido
Considerando los siguientes valores:
M = 6103,5 kg/h.
Nt = 4136.
nt = 2.
di = 0,027 m.
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:
Tendremos:
Dc = 2,286 m.
B = 0,3 m.
Pt = 1 9/16 plg = 0,040 m.
dt = 1 plg. = 0,032 m.
42
[. / ( *]
{ }
[ ]
Resultando: Deq = 0,137 m
{ } { }
{ } { }
{ }
Es decir, que:
43
( * ( * ( * ( *
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.
t = espesor del tubo.
{. / . / . /}
{[. / . /] [. / ] . /}
El coeficiente global U se determinar a partir de la ecuacin anterior.
Siendo:
hi = 23,17 kcal/hm2K.
he = 76,30 kcal/hm2K .
de = 0,032 m.
di = 0,027 m.
F = factor de ensuciamiento, se asignar el valor de 0,010 hm 2K/kcal tanto
para la superficie exterior como para la interior.
Siendo:
44
Q = 15294696,2 kJ/h = 3655520,12 kcal/h.
U = 11,82 kcal/hm2K.
{( ) ( )}
( )
8 9
( )
*( ) ( )+
( )
{ }
( )
Tml = 105,5K.
Despejando el valor de A:
A = 2932,61 m2
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.
45
Para el exterior de los tubos, la determinacin de la prdida de carga se realizar
segn el mtodo:
( )
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:
6.3.2.1. Introduccin
46
Caudal molar = 320774,50 mol/h
Caudal msico = 10126,26 kg/h
Temperatura de entrada = 331C (604 K)
Temperatura de salida = 100C (373 K)
Temperatura media = 215,5C (488,5 K)
Composicin molar:
19,35% Metanol
41,87% Agua
38,78% Dimetil ter
Las propiedades a la temperatura media son las siguientes:
o Cp (media) = 0,574 kcal/kgK
o (media) = 0,051 kg/hm
o (media) = 52,543 kg/m3
( ) ( ) ( )
Donde A B, y C son constantes para el DME, MeOH y agua, estas constantes son:
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.
47
Para llevar a cabo el clculo del calor retirado por la corriente fra deber tenerse en
cuenta los calores latentes de vaporizacin que son, respectivamente, 36,13; 40,60 y
21,51 kJ/mol para cada compuesto. 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 331C y 173,3C as como entre la temperatura
de roco y 100C. Dichos valores se recogen en las siguientes tablas:
2 {[ ( )] [ ( )]}
{[ ( )] [ ( )]} {[ ( )]
[ ( )]}3
{ 20 ( ) 1
0 ( ) 13 20 (
) 1 0 ( ) 13
20 ( ) 1 , ( )-3}
= 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:
( )
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.
( )
48
( )
( )
( )
49
= viscosidad a la temperatura del fluido.
w = viscosidad a la temperatura de la pared del tubo.
CP = capacidad calorfica.
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:
Tendremos:
Dc = 1,676 m.
B = 0,3 m.
Pt = 1 9/16 plg = 0,040 m.
dt = 1 plg. = 0,032 m.
[. / ( *]
{ }
[ ]
Resultando: Deq = 0,137 m
51
6.3.2.7. Clculo del coeficiente global de transmisin de calor
{ } { }
{ } { }
{ }
Es decir, que:
( * ( * ( * ( *
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.
t = espesor del tubo.
{. / . / . /}
{[. / . /] [. / ] . /}
El coeficiente global U se determinar a partir de la ecuacin anterior.
52
Siendo:
hi = 40,92 kcal/hm2K.
he = 109,31 kcal/hm2K .
de = 0,032 m.
di = 0,027 m.
F = factor de ensuciamiento, se asignar el valor de 0,010 hm 2K/kcal tanto
para la superficie exterior como para la interior.
Siendo:
{( ) ( )}
( )
8 9
( )
*( ) ( )+
( )
{ }
( )
Tml = 138,7K.
Despejando el valor de A:
A = 1162,08 m2
53
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.
( )
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:
54
alcanzar valores de conversin ptimos de metanol a ter dimetlico, asegurando las
mejores condiciones de velocidad de reaccin, de prdidas de carga y menores
cantidades de catalizador.
En el diseo del reactor tambin tendrn una gran importancia las caractersticas
fsico-qumicas del catalizador, tales como su densidad real y aparente, el rea
superficial y su porosidad, entre otras. Todos estos factores repercutirn en la
obtencin de mejores rendimientos globales del proceso.
4 ( *5
( ( ) )
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
( )
( *
( )
( ( ) )
( ) ( )
( )
( ) ( )
( )
55
( ) ( )
( )
6.4.2.1. Introduccin
Aunque, a priori, para que fuera mayor el rendimiento interesara eliminar el calor
desprendido durante la deshidratacin (operacin isotrmica), las dificultades que
presenta la eliminacin del calor de la reaccin de la masa cataltica hace que se haya
optado por la operacin en un rgimen aproximadamente adiabtico.
Balance de materia
La expresin de la conservacin de la masa para este tipo de reactor est dada por la
expresin:
que ser la ecuacin del balance de materia para el reactor una vez sustituida la
expresin de la velocidad de reaccin. Entonces, slo es necesario integrar la
ecuacin y as determinar la conversin para cualquier masa de catalizador.
Balance de energa
La expresin que nos relaciona la entalpa de la reaccin (expresada en kJ/mol) con
la temperatura (expresada en K) de la misma es la siguiente:
56
( ) ( ){ } ( )8 9
( )8 9 ( )8 9
La expresin que nos relaciona la cada de presin con la longitud del lecho poroso es
la ecuacin de Ergun:
( )
{ } { } 86 7 9
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.
( )
{ } { } 86 7 9
57
La masa de catalizador, se calcular tomando mP0 = 0 kg, la temperatura T0 =
523K, y la presin P0 = 1610 kPa que hay a la entrada del lecho, aplicando el
mtodo de Runge-Kutta a las ecuaciones con intervalos de longitud de mP = 1 kg
desde el valor inicial de mP = 0 hasta el valor final de mP = 140.
Los valores de las constantes que aparecen en las ecuaciones anteriores son los
siguientes:
T0 = 523 K
P0 = 1610 kPa
M0 = 10126,26 kg/h
r = 0,375 m
Ac = 0,4118 m2
G = 22921,16 kg/m2h
NA0 = 310,3552 kmol/h
NB0 = 10,15707 kmol/h
NC0 = 0,26238 kmol/h
= 57,86 kJ/kmolK
= 37,28 kJ/molK
= 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.
58
140 0,0840 1595,51 542,91
150 0,0960 1595,22 544,32
160 0,1030 1595,05 545,72
170 0,1100 1594,67 547,12
180 0,1200 1594,56 548,51
190 0,1310 1594,23 549,91
200 0,1390 1594,16 551,30
210 0,1510 1594,03 552,68
220 0,1620 1593,85 554,07
230 0,1750 1593,67 555,45
240 0,1890 1593,52 556,84
250 0,2030 1593,33 558,21
260 0,2170 1593,04 559,59
270 0,2370 1592,95 560,97
280 0,2510 1592,82 562,34
290 0,2720 1592,66 563,71
300 0,2930 1592,45 565,07
310 0,3100 1592,33 566,44
320 0,3310 1592,27 567,80
330 0,3550 1592,16 569,16
340 0,3720 1592,04 570,52
350 0,4000 1591,97 571,87
360 0,4210 1591,77 573,22
370 0,4380 1591,56 574,57
380 0,4550 1591,48 575,92
390 0,4760 1591,35 577,27
400 0,4970 1591,25 578,61
410 0,5140 1591,10 579,95
420 0,5280 1591,02 581,29
430 0,5450 1590,95 582,63
440 0,5620 1590,88 583,96
450 0,5760 159076 585,29
460 0,5890 1590,62 586,62
470 0,6000 1590,55 587,95
480 0,6170 1590,44 589,26
490 0,6280 1590,37 590,60
500 0,6400 1590,27 592,05
510 0,6520 1590,19 594,34
520 0,6620 1590,07 595,55
530 0,6720 1590,00 595,87
540 0,6830 1589,95 598,18
550 0,6890 1589,92 598,49
560 0,7000 1589,84 600,80
570 0,7050 1589,75 602,10
580 0,7140 1589,67 604,40
59
590 0,7240 1589,59 606,70
600 0,7280 1589,50 608,00
610 0,7340 1589,42 609,30
620 0,7390 1589,34 607,59
630 0,7480 1589,26 608,89
640 0,7520 1589,18 610,17
650 0,7590 1589,10 611,46
660 0,7630 1589,02 612,75
670 0,7690 1588,95 614,03
680 0,7750 1588,87 615,31
690 0,7795 1588,79 616,59
700 0,7840 1588,71 617,87
710 0,7870 1588,63 619,14
720 0,7910 1588,55 620,42
730 0,7930 1588,47 621,69
740 0,7960 1588,40 622,96
750 0,7970 1588,32 624,22
760 0,7980 1588,24 625,49
770 0,7985 1588,15 626,75
780 0,7990 1588,07 628,01
790 0,7995 1587,99 629,27
800 0,8000 1587,90 630,52
Tabla 6.6. Variacin de las composiciones molares de los tres componentes con la
altura del lecho cataltico.
60
180 0,8475 0,0917 0,0608
190 0,8365 0,0972 0,0663
200 0,8285 0,1012 0,0703
210 0,8165 0,1072 0,0763
220 0,8055 0,1127 0,0818
230 0,7925 0,1192 0,0883
240 0,7785 0,1262 0,0953
250 0,7645 0,1332 0,1023
260 0,7505 0,1402 0,1093
270 0,7305 0,1502 0,1193
280 0,7165 0,1572 0,1263
290 0,6955 0,1677 0,1368
300 0,6745 0,1782 0,1473
310 0,6575 0,1867 0,1558
320 0,6365 0,1972 0,1663
330 0,6125 0,2092 0,1783
340 0,5955 0,2177 0,1868
350 0,5675 0,2317 0,2008
360 0,5465 0,2422 0,2113
370 0,5295 0,2507 0,2198
380 0,5125 0,2592 0,2283
390 0,4915 0,2697 0,2388
400 0,4705 0,2802 0,2493
410 0,4535 0,2887 0,2578
420 0,4395 0,2957 0,2648
430 0,4225 0,3042 0,2733
440 0,4055 0,3127 0,2818
450 0,3915 0,3197 0,2888
460 0,3785 0,3262 0,2953
470 0,3675 0,3317 0,3008
480 0,3505 0,3402 0,3093
490 0,3395 0,3457 0,3148
500 0,3275 0,3517 0,3208
510 0,3155 0,3577 0,3268
520 0,3055 0,3627 0,3318
530 0,2955 0,3677 0,3368
540 0,2845 0,3732 0,3423
550 0,2785 0,3762 0,3453
560 0,2675 0,3817 0,3508
570 0,2625 0,3842 0,3533
580 0,2535 0,3887 0,3578
590 0,2435 0,3937 0,3628
600 0,2395 0,3957 0,3648
610 0,2335 0,3987 0,3678
620 0,2285 0,4012 0,3703
61
630 0,2195 0,4057 0,3748
640 0,2155 0,4077 0,3768
650 0,2085 0,4112 0,3803
660 0,2045 0,4132 0,3823
670 0,1985 0,4162 0,3853
680 0,1925 0,4192 0,3883
690 0,1880 0,42145 0,39055
700 0,1835 0,4237 0,3928
710 0,1805 0,4252 0,3943
720 0,1765 0,4272 0,3963
730 0,1745 0,4282 0,3973
740 0,1715 0,4297 0,3988
750 0,1705 0,4302 0,3993
760 0,1695 0,4307 0,3998
770 0,1690 0,43095 0,40005
780 0,1685 0,4312 0,4003
790 0,1680 0,43145 0,40055
800 0,1675 0,4317 0,4008
Tabla 6.7. Variacin de los caudales molares de los tres componentes con la altura
del lecho cataltico.
62
200 8389,61 1024,78 711,88
210 8268,09 1085,54 772,63
220 8156,70 1141,23 828,33
230 8025,06 1207,05 894,15
240 7883,29 1277,93 965,03
250 7741,53 1348,82 1035,92
260 7599,76 1419,70 1106,80
270 7397,23 1520,96 1208,06
280 7255,47 1591,85 1278,95
290 7042,81 1698,17 1385,27
300 6830,16 1804,50 1491,60
310 6658,02 1890,57 1577,67
320 6445,36 1996,90 1684,00
330 6202,33 2118,41 1805,51
340 6030,19 2204,49 1891,59
350 5746,65 2346,25 2033,35
360 5534,00 2452,58 2139,68
370 5361,85 2538,65 2225,75
380 5189,71 2624,73 2311,83
390 4977,06 2731,05 2418,15
400 4764,41 2837,38 2524,48
410 4592,26 2923,45 2610,55
420 4450,49 2994,34 2681,43
430 4278,34 3080,41 2767,51
440 4106,20 3166,48 2853,58
450 3964,43 3237,37 2924,46
460 3832,79 3303,19 2990,28
470 3721,40 3358,88 3045,98
480 3549,25 3444,95 3132,05
490 3437,87 3500,65 3187,75
500 3316,35 3561,41 3248,50
510 3194,84 3622,16 3309,26
520 3093,57 3672,79 3359,89
530 2992,31 3723,43 3410,52
540 2880,92 3779,12 3466,22
550 2820,16 3809,50 3496,60
560 2708,77 3865,19 3552,29
570 2658,14 3890,51 3577,61
580 2567,01 3936,08 3623,18
590 2465,74 3986,71 3673,81
600 2425,24 4006,96 3694,06
610 2364,48 4037,34 3724,44
620 2313,85 4062,66 3749,75
63
630 2222,71 4108,22 3795,32
640 2182,21 4128,48 3815,57
650 2111,33 4163,92 3851,02
660 2070,82 4184,17 3871,27
670 2010,06 4214,55 3901,65
680 1949,31 4244,93 3932,03
690 1903,74 4267,71 3954,81
700 1858,17 4290,50 3977,59
710 1827,79 4305,69 3992,78
720 1787,28 4325,94 4013,04
730 1767,03 4336,06 4023,16
740 1736,65 4351,25 4038,35
750 1726,53 4356,32 4043,42
760 1716,40 4361,38 4048,48
770 1711,34 4363,91 4051,01
780 1706,27 4366,44 4053,54
790 1701,21 4368,97 4056,07
800 1696,15 4371,51 4058,61
Para estimar el volumen total del reactor, se va a considerar el lecho del reactor como
un cilindro de base AC y altura L, que se encuentra relleno por el catalizador poroso,
de modo que todo el cilindro estar ocupado tato por el volumen del catalizador como
por el volumen de huecos. El volumen del cilindro considerado ser:
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:
64
6.5. Diseo de las Columnas de Destilacin
El DME, MeOH y agua son los principales componentes que ingresan a la torre
de destilacin. A continuacin se define los componentes clave:
Clave liviano: DME.
Clave pesado: Metanol.
Donde A B, y C son constantes para el DME, MeOH y agua, estas constantes son:
65
A B C
MeOH 7,878 1473,11 230,00
Agua 7,967 1668,21 228,00
DME 7,400 1120,00 241,96
( ) ( ) ( )
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.
( ) ( ) ( )
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.
66
El clculo de la volatilidad relativa se realiz con la temperatura promedio
(TProm) de las temperaturas de burbuja y roco.
TProm = 133,60 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.
2 3
67
6.5.1.6. Clculo de la Relacin de Reflujo Mnimo
( )
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.
q = 1,0
XDME,Destilado = 0,991
XMeoH,Fondo = 0,313
= 0,9719
Vmin/D = 1,530
Lmin/D = 0,530
Rmin = 0,530
68
que los mnimos, estos parmetros pueden ser calculados si uno de ellos es
fijado como mltiplo, obteniendo el otro por el medio de la correlacin;
( )
( *
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
( ) ( )
69
{ }
Donde:
V = 223560,80 mol/h
vap,m = 36,13 kJ/mol
vap,a = 40,50 kJ/mol
vap,d = 21,51 kJ/mol
{ }
( ) ( )
{ }
Donde:
V = 223560,91 mol/h
{ }
70
Metanol 0,313
Agua 0,682
Dimetil ter 0,005
Total 1,000
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.
( ) ( ) ( )
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.
71
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.
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.
2 3
( )
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.
73
XH2O,Fondo = 0,790
= 0,7443
Vmin/D = 2,6928
Lmin/D = 1,6928
Rmin = 1,6928
( )
( *
La columna fue designada para operar con un reflujo actual de 10% 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.
(L/D)min = 1,6928
L/D = 1,5(L/D)min = 2,5392
Abscisa = 0,2391
[N-Nmin]/(N+1) = 0,522033724
Nmin = 2,533
o N = 6,40
74
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.
( ) ( )
{ }
Donde:
V = 136731,89 mol/h
vap,m = 36,13 kJ/mol
vap,a = 40,50 kJ/mol
vap,d = 21,51 kJ/mol
{ }
( ) ( )
{ }
Donde:
V = 223560,91 mol/h
{ }
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.
75
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.
Resultando as:
hL = 1,389 m.
76
vA = 0 m/s.
Y por lo tanto:
( )
( ) [ ] ( )
Donde:
Sustituyendo los valores en la expresin anterior, se obtiene que las bombas deben
tener una NPSH de 11,24 m.
77
6.6.2. Diseo de las Bombas 3 y 4
( * 4 5 ( * 4 5
Y por lo tanto:
78
7. ANLISIS FINANCIERO DEL PROYECTO
79
Gerencia General
Secretaria
Asesora Legal
Ejecutiva
Gerencia de
Gerencia Gerencia de
Aministracin y Gerencia Tcnica
Comercial Mantenimiento
Finanzas
Departamento de
Secretaria Secretaria Secretaria
Mantenimiento
Departamento de Departamento de
Departamento de Departamento de Departamento de Seguridad
Recursos Adquisicin de Planta Laboratorio Supervisin
Contabilidad Ventas Publicidad Industrial
Humanos Materia Prima
80
7.1.3. Requerimiento de personal y costo de personal
El rea tcnica operacional contar con tres turnos de trabajo de 8 horas por turno
con el siguiente personal:
Supervisor
Operador de sala
Operador de campo
81
Ayudante de campo
Se considera un operador extra para cubrir vacaciones y otras contingencias. En la
Tabla 7.2. se muestra el nmero aproximado de operarios para el rea tcnica de la
planta tambin se detalla el salario del personal de campo y jefes de la planta.
El rea tcnica de la planta requiera un total de 99 empleados, y tiene un costo total
de MM$us/Ao.
82
que la produccin efectiva sea igual al 75% de la capacidad total proyectada para este
proyecto.
La Tabla 7.4. muestra los resultados de los clculos efectuados del costo de
inversin total estimado bajo las condiciones supuestas anteriormente:
83
Para el monto de transporte de equipos se tom el 10 % del costo de los
mismos, dando un total de 0,3 MM$us que es el monto total de compra de
equipo puesto en Patacamaya.
La instalacin y montaje se asumieron como 40 % de la compra de equipo.
Para las tuberas, vlvulas y accesorios se tom un 60 % de la compra de
equipo.
La instrumentacin para la automatizacin de la planta representa un 8 % de
la compra de equipo.
El sistema elctrico se considera como 10 % de la compra de equipo.
Los servicios auxiliares son 30 % dela compra de equipo.
La infraestructura y obras civiles se asumieron como un 5 % de la compra de
equipo.
El capital de trabajo representa el 15% del activo fijo [Peters et al, 2002].
84
Metanol 78.998 0,41574 32,84
Total 78.998 0,41574 32,84
Donde:
frc = Factor de recuperacin de capital
P = Monto de prstamo
El factor de recuperacin de capital se calcula con la ecuacin [Sapag, 2007]:
85
( )
( )
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.
Los costos de operacin se clasifican en costos fijos y costos variables. Los costos
variables calculados se detallan en la Tabla 7.8., y en la Tabla 7.9. se detallan los
principales costos fijos de operacin.
86
Descripcin Monto (Miles de $us/Ao)
Costo de Produccin
Mano de Obra Indirecta 0,30729984
Mantenimiento 0,100
Laboratorio 0,032
Directorio 0,064
Seguros 0,083
Subtotal 0,586
Costos Generales
Gerencia 0,812
Investigacin y Servicios 0,271
Subtotal 1,082
Total 1,668
Los costos fijos de la planta de DME fueron calculados en base a los siguientes
criterios:
La Tabla 7.10. muestra un resumen de los costos fijos y variables. La Tabla 7.11.
muestra los costos del proyecto incluidos los costos financieros, la depreciacin y la
amortizacin de los activos diferidos.
87
Tabla 7.11. Costos del Proyecto en MM$us/Ao
DESCRIPCION 1 2 3 4 5 6 7 8 9 10
Programa de
75 100 100 100 100 100 100 100 100 100
Produccin
Produccin 37500 50000 50000 50000 50000 50000 50000 50000 50000 50000
Costos Variables 26,648 34,858 34,862 34,866 34,870 34,874 34,881 34,889 34,897 34,905
Costo de Materia
24,632 32,843 32,843 32,843 32,843 32,843 32,843 32,843 32,843 32,843
Prima
Costo de Insumos 0,794 0,794 0,794 0,794 0,794 0,794 0,794 0,794 0,794 0,794
Costo de MOD 0,739 0,739 0,743 0,747 0,751 0,755 0,762 0,770 0,778 0,786
Gastos de
0,483 0,483 0,483 0,483 0,483 0,483 0,483 0,483 0,483 0,483
Comercializacin
Costos Fijos 1,668 1,668 1,672 1,676 1,680 1,684 1,691 1,699 1,707 1,715
Mano de obra
0,307 0,307 0,311 0,315 0,319 0,323 0,331 0,338 0,346 0,354
indirecta
Directorio 0,064 0,064 0,064 0,064 0,064 0,064 0,064 0,064 0,064 0,064
Laboratorio 0,032 0,032 0,032 0,032 0,032 0,032 0,032 0,032 0,032 0,032
Seguros 0,083 0,083 0,083 0,083 0,083 0,083 0,083 0,083 0,083 0,083
Investigacin y
0,271 0,271 0,271 0,271 0,271 0,271 0,271 0,271 0,271 0,271
servicios
Gerencia 0,812 0,812 0,812 0,812 0,812 0,812 0,812 0,812 0,812 0,812
Mantenimiento 0,100 0,100 0,100 0,100 0,100 0,100 0,100 0,100 0,100 0,100
COSTO DE
28,316 36,526 36,534 36,542 36,549 36,557 36,573 36,588 36,604 36,619
OPERACION
Costos financieros 0,420 0,387 0,352 0,315 0,276 0,235 0,193 0,148 0,101 0,052
Depreciacin Activos
0,282 0,282 0,282 0,282 0,282 0,282 0,282 0,282 0,282 0,282
Fijos
Amortizacin activos
1,872 1,872 1,872
diferidos
COSTO DE
30,89 39,07 39,04 37,14 37,11 37,07 37,05 37,02 36,99 36,95
PRODUCCION
88
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.
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.
89
Tabla 7.13. Estado de Resultados del Proyecto en MM$us/Ao
DESCRIPCION 1 2 3 4 5 6 7 8 9 10
1. INGRESOS 31,2375 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65
Ventas 31,2375 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65
Valor Residual
2. EGRESOS 31,558 39,768 39,776 37,912 37,919 37,927 37,943 37,958 37,974 37,989
Costo de Operacin 28,316 36,526 36,534 36,542 36,549 36,557 36,573 36,588 36,604 36,619
Intereses 0,420 0,387 0,352 0,315 0,276 0,235 0,193 0,148 0,101 0,052
Amortizacin 0,668 0,701 0,736 0,773 0,812 0,852 0,895 0,940 0,987 1,036
Depreciacin 0,282 0,282 0,282 0,282 0,282 0,282 0,282 0,282 0,282 0,282
Amortizacin Diferida 1,872 1,872 1,872 0,000 0,000 0,000 0,000 0,000 0,000 0,000
UTILIDAD BRUTA -0,320 1,882 1,874 3,738 3,731 3,723 3,707 3,692 3,676 3,661
Impuestos IVA (13%)
Impuestos IT (2,5%) 0,781 1,041 1,041 1,041 1,041 1,041 1,041 1,041 1,041 1,041
UTILIDAD IMPONIBLE 0,000 0,841 0,833 2,697 2,689 2,682 2,666 2,650 2,635 2,619
Impuesto Utilidad IUE
0,000 0,210 0,208 0,674 0,672 0,670 0,666 0,663 0,659 0,655
(25%)
UTILIDAD NETA -1,101 0,630 0,625 2,023 2,017 2,011 1,999 1,988 1,976 1,965
90
Tabla 7.14. Corrientes de Liquidez para la Planificacin Financiera en MM$us/Ao
DESCRIPCION 0 1 2 3 4 5 6 7 8 9 10
1. FUENTES 8,4 31,2375 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65 42,443
Aporte Propio
Crdito 8,4
Ventas 31,2375 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65 41,65
Valor Residual
Recuperacin de K
0,793
de Trabajo
2. USOS 8,296 30,184 38,865 38,871 39,345 39,351 39,357 39,368 39,380 39,392 39,403
Inversiones 8,296
Costo de
28,316 36,526 36,534 36,542 36,549 36,557 36,573 36,588 36,604 36,619
Operacin
Intereses 0,420 0,387 0,352 0,315 0,276 0,235 0,193 0,148 0,101 0,052
Amortizacin 0,668 0,701 0,736 0,773 0,812 0,852 0,895 0,940 0,987 1,036
Impuestos 0 0,781 1,251 1,249 1,716 1,714 1,712 1,708 1,704 1,700 1,696
EXCEDENTE /
0,104 1,053 2,785 2,779 2,305 2,299 2,293 2,282 2,270 2,258 3,040
DEFICIT
ACUMULADO 0,104 1,157 3,941 6,720 9,025 11,324 13,618 15,899 18,169 20,428 23,467
91
Tabla 7.15. Flujo de Fondos Neto en MM$us/Ao
DESCRIPCION 0 1 2 3 4 5 6 7 8 9 10
1. INGRESOS 0 31,238 42 42 42 42 42 42 42 42 42
Ventas 31,238 42 42 42 42 42 42 42 42 42
Valor Residual
Recuperacin de
K de Trabajo
2. EGRESOS 8,296 30,184 38,865 38,871 39,345 39,351 39,357 39,368 39,380 39,392 39,403
Inversiones 8,296 0 0 0 0 0 0 0 0 0 0
Costo de
28,316 36,526 36,534 36,542 36,549 36,557 36,573 36,588 36,604 36,619
Operacin
Intereses 0,420 0,387 0,352 0,315 0,276 0,235 0,193 0,148 0,101 0,052
Amortizacin 0,668 0,701 0,736 0,773 0,812 0,852 0,895 0,940 0,987 1,036
Impuestos 0,781 1,251 1,249 1,716 1,714 1,712 1,708 1,704 1,700 1,696
FLUJO NETO -8,296 1,053 2,785 2,779 2,305 2,299 2,293 2,282 2,270 2,258 2,247
FLUJO NETO
-8,296 -7,243 -4,459 -1,680 0,625 2,924 5,218 7,499 9,769 12,028 14,274
ACUMULADO
92
7.7.3. Clculo de indicadores financieros para el proyecto
( )
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.
( )
Donde:
i1 = Tasa de actualizacin del ltimo VAN positivo
93
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 %.
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%.
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.
( )
( )
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
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 %.
94
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.
Donde:
I0 = Inversin inicial
FC = Flujos de caja
t = Periodo de tiempo 1,2,..., n.
Descripcin 1 2 3 4 5 6 7 8 9 10
Flujo Neto -8,296 1,053 2,785 2,779 2,305 2,299 2,293 2,282 2,270 2,258
Flujo Neto
-8,296 -7,243 -4,459 -1,680 0,625 2,924 5,218 7,499 9,769 12,028
Acumulado
95
8. SIMULACIN EN HYSYS
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.
96
97
9. CONCLUSIONES
10. BIBLIOGRAFIA
98
ANEXOS
99
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).
100
Grfica A.1. Factor de transferencia de calor en los tubos
101