Cornils Peter
Cornils Peter
Cornils Peter
DEPARTAMENTO DE INDUSTRIAS
VALPARAÍSO – CHILE
NOVIEMBRE 2014
Universidad Técnica Federico Santa María
Departamento de Industrias
Índice de Contenido
1. Agradecimientos .................................................................................................. 18
2. Resumen .............................................................................................................. 19
3. Abstract................................................................................................................ 20
4. Introducción ......................................................................................................... 21
5. Objetivos.............................................................................................................. 22
6. Alcance ................................................................................................................ 23
2
Universidad Técnica Federico Santa María
Departamento de Industrias
7.2.3. Colectores solares cilíndrico parabólicos (CCP) ...................................... 34
8.3. Energía transportada y parámetros de los fluidos en cada punto de interés ..... 57
3
Universidad Técnica Federico Santa María
Departamento de Industrias
8.4.2. Patio de estanques ................................................................................... 67
9.4. Diseño Básico del Sistema de Calentamiento Solar de Agua ........................ 147
9.5. Dimensionamiento de equipos del sistema de calentamiento solar de agua ... 168
4
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.4. Análisis de Costos ................................................................................ 194
5
Universidad Técnica Federico Santa María
Departamento de Industrias
13.1.1. Balances de Masa .............................................................................. 264
6
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.4. Método de la Utilizabilidad ............................................................... 292
7
Universidad Técnica Federico Santa María
Departamento de Industrias
Índice de Tablas
Tabla 16 - Energía disipada por efecto Joule por componentes de Nave EW ................ 84
8
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 17 - Entalpías o calores de formación de cada compuesto involucrado en al
reacción de electroobtención ........................................................................................ 85
Tabla 23 - Densidad del agua en el aire sobre las celdas de electroobtención .............. 106
Tabla 25 - Densidad del agua en el aire sobre el estanque de recirculación ................. 123
Tabla 29 - Factor de espaciamiento y Área bruta corregida de cada colector solar ...... 140
9
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 34 - Precios de la Electricidad Medios en Atacama (SING) .............................. 190
Tabla 36 - Proyecciones de los Precios del Gas Natural y Diesel desde el 2014 hasta el
2024 [24] ................................................................................................................... 192
Tabla 37 - Costo Total de Combustibles para los casos base evaluados ...................... 193
Tabla 40- Costo total por Viajes y Alojamiento por la Inspección de Sitio ................. 201
Tabla 41 - Costo total por Viajes y Alojamiento por el Estudio de Impacto Ambiental
.................................................................................................................................. 201
Tabla 42 - Costo total por Viajes y Alojamiento por la Topografía del Sitio ............... 205
Tabla 43 - Costo total por Viajes y Alojamiento por la Validación y Registro de GEI 205
Tabla 44 - Costo total por Viajes y Alojamiento por la Gerencia ................................ 205
Tabla 46 - Precio unitario del área de apertura de los Colectores Solares Evaluados .. 210
Tabla 49 - Resumen de Resultados: Calentamiento del Agua del Circuito de Calderas 236
10
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 51 - Resumen de Resultados: Calentamiento de Agua de Osmosis – Normal .... 238
Tabla 54 - Variables más influyentes en el VAN por caso base - Fracción solar 100% 250
11
Universidad Técnica Federico Santa María
Departamento de Industrias
Índice de Figuras
Figura 2 - Componentes básicos de un tubo de vacío típico – Imagen: Estado del arte
colectores [5] ............................................................................................................... 33
Figura 8 - Vista Satelital de la División Radomiro Tomic en Enero del 2014 (Extraída de
Google Earth) .............................................................................................................. 41
12
Universidad Técnica Federico Santa María
Departamento de Industrias
Figura 15- Volumen de Control: Intercambiador de Calor - Electrolito Rico y Electrolito
Pobre ........................................................................................................................... 70
Figura 18 - Perfil tentativo de la Convección del aire sobre celdas de la Nave Ew ........ 88
Figura 25 - Cara Lateral Ancha de Celda (Plano perpendicular a eje X) ..................... 107
Figura 28 - Cara Lateral Larga de Celda (Plano perpendicular a eje Z) ....................... 112
13
Universidad Técnica Federico Santa María
Departamento de Industrias
Figura 33 - Pérdidas de energía de los Estanques ....................................................... 120
Figura 37 - Diseño Básico de los Sistemas de Calentamiento Solar de Agua .............. 148
Figura 39 - Alternativa: Sistema de calentamiento solar de agua para pre calentar agua
de osmosis ................................................................................................................. 153
Figura 41 - Alternativa: Sistema de calentamiento solar de agua para pre calentar agua
de calderas ................................................................................................................. 155
Figura 42 - Perfil de Temperaturas: Intercambiador de Calor – Loop Solar / Carga .... 157
Figura 10 - Transferencia de calor entre dos fluidos separados por una placa ............. 269
14
Universidad Técnica Federico Santa María
Departamento de Industrias
Índice de Gráficos
Gráfico 1 - Evolución del Producto Interno Bruto de Chile, elaboración propia con
información de The World Bank. ................................................................................. 24
Gráfico 12 - Eficiencia Instantánea del Colector Solar versus Temperatura de entrada 144
Gráfico 16 - Eficiencia estacional del sistema de calderas según fracción solar - Gas
Natural....................................................................................................................... 166
15
Universidad Técnica Federico Santa María
Departamento de Industrias
Gráfico 17- Eficiencia estacional del sistema de calderas según fracción solar - Diesel
N° 2 ........................................................................................................................... 166
Gráfico 18 - Radiación solar en el plano del Colector Solar – Fijado .......................... 184
Gráfico 21 - Estructura Relativa de Costos de Inversión Solar (Según IREC) ............. 194
Gráfico 22 - Función de costo de cañerías de 8'' - Acero inoxidable 316L Sch. 40 ...... 212
Gráfico 25 - Función de costos de bombas de acero inoxidable - hasta 1035 kPa........ 216
Gráfico 28 - Evolución del precio de los EUA entre los años 2009 y 2014 [29] .......... 233
Gráfico 30 - Costo de inversión por área de parque solar versus área total del parque
solar ........................................................................................................................... 241
Gráfico 34 - Contribución a la varianza del VAN - Caso Base: Diesel – f : 7,7%........ 248
Gráfico 35 - Contribución a la varianza del VAN - Caso Base: Diesel – f: 100%........ 248
Gráfico 36 - Contribución a la varianza del VAN - Caso Base: Gas Natural – f : 7,7%
.................................................................................................................................. 249
Gráfico 37 - Contribución a la varianza del VAN - Caso Base: Gas Natural – f: 100%
.................................................................................................................................. 250
Gráfico 38 - VAN versus Precio del Gas Natural y Tipo de Cambio Euro/Dolar ........ 252
Gráfico 39 - VAN versus Precio del Gas Natural y Precio del Colector Solar ............. 253
Gráfico 40 - VAN versus Precio del Colector solar y Tipo de Cambio Euro/Dolar ..... 254
Gráfico 42 - VAN versus vida útil del proyecto de precalentar agua de calderas usando
en el caso base de combustible gas natural ................................................................. 256
Gráfico 43 - VAN versus vida útil del proyecto de precalentar agua de calderas usando
en el caso base de combustible diesel ......................................................................... 257
17
Universidad Técnica Federico Santa María
Departamento de Industrias
1. Agradecimientos
Quiero agradecer en primer lugar a mis padres Sven y Ruth, por todo el apoyo
incondicional y amor que me han hecho sentir a lo largo de toda una vida, por las
grandes experiencias que gracias a ustedes he podido vivir. Ustedes son la base del
edificio que hace 27 años llevo construyendo y solo gracias a eso tengo la certeza que
algún día llegará a ser muy alto.
Quiero agradecer a mi hermana Carol, por ser mi ejemplo y mi patrón a imitar desde
pequeño; por demostrarles a todos que con dedicación y perseverancia puedes lograr
cualquier objetivo que te propongas. Agradecerte también porque a pesar de la distancia,
siempre has hecho notar tu cariño y preocupación.
Quiero agradecer a todos mis amigos que conocí durante este ciclo, que ha sido una
etapa increíble en mi vida, y ustedes son una de las principales razones de ello. A
Francisco, Claudio, Antonella, Andrés, Aníbal y Eileen, porque además de todas las
buenas vibras entregadas, y que a pesar de que uno está al otro lado del mundo, ustedes
son de esos amigos que cuesta encontrar y que pase lo que pase hacen sentir que siempre
se puede contar con ustedes. Ustedes más que amigos, son familia; son hermanos.
Me voy con una sonrisa en la cara, sintiéndome pleno y listo para disfrutar de una vida
que ha sido y seguirá siendo maravillosa, y eso se lo debo a todos ustedes, por lo que lo
mínimo que puedo decir es:
Muchísimas Gracias.
18
Universidad Técnica Federico Santa María
Departamento de Industrias
2. Resumen
En el presente estudio se realizó una serie de evaluaciones de proyectos para la
incorporación de colectores solares térmicos como sistema de calentamiento de fluidos
para los procesos de hidrometalurgia. Para dicho propósito, se evaluó precalentar en un
caso el agua del circuito de calderas, en otro caso precalentar directamente el electrolito
rico previo a la etapa de electroobtención y en otro caso precalentar el agua de osmosis
para el lavado de cátodos, riego de celdas u otras demandas esporádicas de energía.
19
Universidad Técnica Federico Santa María
Departamento de Industrias
3. Abstract
In the present study, series of project evaluations incorporating solar thermal collectors
as a fluid heating system for hydrometallurgical processes were performed. For
achieving this objective, an evaluation of preheating process fluids in differents circuits
of the hydrometallurgy was held. These were the boilers system water circuit, the rich
electrolyte before entering to the electrowinning stage and the osmosis water principally
used for cathodes washing, irrigating electrolyte cells and other small demands.
Each of the project evaluations was performed using 26 models of solar collectors with 5
different technology types; flat solar collectors, evacuated tube solar collectors, a fresnel
colector, parabolic trough collectors and a concentrating with stationary reflector
collector. Their respective levelized energy costs were 51, 54, 63, 76 USD/MWh.
The most profitable project, with an project life cycle equal to 5 years, was considered to
preheat the boilers system water, and obtained an VPN equal to 1,57 M USD, and IVAN
equal to 31%, an simple capital return on the third year of the project and an negative
TIR equal to -23,8%. The mean capacity of this preheating system is up to 10.868
MWh/year, which with an solar fraction equal to 8,1%, contemplates an solar field gross
area equal to 21.064 m2 using 4.457 flat solar collectors. The same proyect with an
project life cycle equal to 10 years, improves his profits with an VPN equal to 5,69 M
USD, an IVAN equal to 111% and an TIR equal to 13,5%.
An sensibility analysis was performed, which proved the importance of the fossil
combustibles price against the projects profit. This is so because all of the evaluated
projects stopped being profitable when used natural gas as the boiler system combustible
of the base project line.
20
Universidad Técnica Federico Santa María
Departamento de Industrias
4. Introducción
Entre los años 2001 y 2010 la producción de cobre en la industria minera de Chile
aumentó en un 14%, desde 4,74 TMF hasta 5,42 TMF, observándose un importante
aumento de 50% en su consumo total de energía, desde 86.234 TJ hasta 129.583 TJ [1].
Algunos de los factores más significativos que explican el continuo incremento en los
costos medios del cobre son: la disminución de la ley promedio del mineral, que en total
en Chile ha disminuido desde 1,29% hasta 0,86% entre los años 2000 y 2012, el
incremento de la dureza del mineral, que entre los mismos años ha incrementado desde
15,4 kWh/T hasta 18,7 kWh/T y por último el incremento de las distancias de acarreo de
material debido a la expansión de las excavaciones. Acorde al envejecimiento de las
minas estos factores se ven incrementados, contribuyendo a que la demanda de insumos
(agua, energía, entre otros), necesarios para producir la misma cantidad de cobre fino,
aumente significativamente [2].
También hay que considerar que el impacto ambiental que conlleva el masivo consumo
de energía en la industria de la minería del cobre de Chile no es menor. Las emisiones
directas de gases de efecto invernadero (GEI) para el año 2010 resultaron ser iguales a
4,3 millones de TM CO2 equivalente, significando un incremento de 49% en una década
[1].
El norte de Chile posee un potencial enorme para la producción de energía solar, ya que
posee uno de los índices de radiación más altos del mundo. Considerando esto, las
presiones al alza en el precio de la energía y los costos medios de producción del cobre,
la producción de GEI, es que resulta interesante buscar alternativas para producir ERNC
para la minería en el norte de Chile.
21
Universidad Técnica Federico Santa María
Departamento de Industrias
5. Objetivos
Objetivo General
Objetivos Específicos
22
Universidad Técnica Federico Santa María
Departamento de Industrias
6. Alcance
Dado que existe un alto grado de similitud entre los parámetros de operación
térmicos de las líneas de procesos entre la mayoría de las plantas de hidrometalurgia del
cobre, es que se decidió analizar en detalle a una planta de hidrometalúrgica real que
permita servir de ejemplo para una posible adaptación a otras plantas. Para ello se
escogió la División Radomiro Tomic de CODELCO.
23
Universidad Técnica Federico Santa María
Departamento de Industrias
300 25
250
200
15
150
10
100
5
50
0 0
1970 1975 1980 1985 1990 1995 2000 2005 2010
Gráfico 1 - Evolución del Producto Interno Bruto de Chile, elaboración propia con información de The World
Bank.
A su vez, como se puede observar en el Gráfico 1, el aporte al PIB por las rentas de
la minería incrementa su tasa de crecimiento fuertemente a partir del año 2003 para
mostrar una fuerte caída el año 2007, lo que se puede explicar por una importante
disminución de la cantidad demandada total del cobre de Europa, la cual cayó en
aproximadamente 500 kTM/año dentro de un período de 3 años (sumado la continua
contracción de la demanda de EEUU que se ha visto progresivamente disminuyendo
desde el año 2000). Sin embargo, ésta baja en la demanda las rentas por la minería
Chilena se recuperan fuertemente el 2009 al verse compensada por la creciente demanda
de China, que incrementó en aproximadamente 4.000 kTM/año desde el año 2007 al
24
Universidad Técnica Federico Santa María
Departamento de Industrias
2011. El 2012 esto transformó a China como el mayor demandante de cobre del mundo,
alcanzando una demanda total de 8.840 kTM, seguido por la Eurozona con 2.601 kTM y
EEUU con 1.758 kTM, siendo el total mundial igual a 20.183 kTM [3].
Por otro lado, la tasa de inflación se ha mantenido positiva los últimos 5 años, y
considerando que la variación del IPC para el año 2014 sea igual a 5,5%, la variación
media del IPC de los últimos 5 años rondaría un valor igual a 3,3% (Gráfico 2).
6,00%
5,00%
Variación IPC (%)
4,00%
Media = 3,3%
3,00%
2,00%
1,00%
0,00%
2009 2010 2011 2012 2013 2014 2015
La producción histórica de cobre al año 2010 fue igual a 5,4 M MTF, siendo un 38,5%
de esta correspondiente a cátodos de cobre producidos por la hidrometalurgia (Gráfico
3).
25
Universidad Técnica Federico Santa María
Departamento de Industrias
6,0
5,5
5,0
El año 2014 la producción de cobre alcanzaría un valor aproximado de 5,8 M TMF [5],
proyectándose que la producción para el año 2020 (considerando los proyectos en
estados de prefactibilidad, factibilidad, ejecución y operativos) alcanzaría un valor
aproximado de 8,1 M TMF, significando un crecimiento del 6,8% anual de la industria
minera (Gráfico 4).
9,0
8,0
7,0
Producción Cu (M TMF)
6,0
Prefactibilidad
5,0
Factibilidad
4,0
Ejecución
3,0 Operaciones
2,0
1,0
0,0
2012 2013 2014 2015 2016 2017 2018 2019 2020
26
Universidad Técnica Federico Santa María
Departamento de Industrias
Con estos antecedentes, se puede afirmar que la industria minera chilena presenta un
constante crecimiento a la par del PIB, y considerando la estabilización de la demanda
internacional del cobre, se pronostica que la industria minera del cobre de Chile
aproximadamente aumentará su producción en un 40 % para el año 2020.
Al año 2012 la ley del cobre en la minería del cobre chileno promedia un valor del
0,86% (Gráfico 5).
1,25
1,20 Concentradora /
Concentrating Plant
1,15
Lixiviación / Leaching
1,10
Ley del Cobre (%)
0,95
0,90
0,85
0,80
2003 2004 2005 2006 2007 2008 2009 2010 2011 2012
27
Universidad Técnica Federico Santa María
Departamento de Industrias
7.1.4. Consumo de Energía
9.000
8.000
7.000
Consumo (MJ/TMF)
6.000
5.000
4.000
3.000
2.000
1.000
0
2003 2004 2005 2006 2007 2008 2009 2010 2011 2012
Mina Rajo / Open Pit Mina Subterránea / Underground Mine
Concentradora / Concentrating Plant Fundición / Smelter
Refinería / Electrolytic Refining LX / SX / EW
Servicios / Services
Gráfico 6 - Valores históricos del Consumo de combustible por unidad de cobre producida
28
Universidad Técnica Federico Santa María
Departamento de Industrias
45
35
30
25
20
15
2012 2013 2014 2015 2016 2017 2018 2019 2020
Según se proyecta, la minería del cobre en Chile demandará un total de energía cercano
a 39 TWh para el año 2020, significando un incremento del 52%.
Según la Comisión Nacional de Energía de Chile (CNE), el precio del diesel rondará por
los 62,4 USD/MWh para el año 2024 y el del gas natural valdría alrededor de los 37,8
USD/MWh para el mismo año. Esto significaría que el diesel experimentaría una tasa de
crecimiento anual promedio igual a 1,2% y el gas natural de -1,3% (Gráfico 8).
65,0
Precio Energía (USD/MWh)
60,0
55,0
50,0
45,0 Diesel
40,0
35,0 Gas Natural
30,0
25,0
20,0
2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024
29
Universidad Técnica Federico Santa María
Departamento de Industrias
Si bien es de esperar que en el mercado existan presiones al alza para todos los
combustibles fósiles, las nuevas tecnologías para extraer shale gas han permitido un
incremento en la oferta de gas natural de E.E.U.U. en un 30 % [4]. Como resultado, en el
corto plazo el precio del gas natural tiene presiones tanto al alza como a la baja, por lo
que en el corto-mediano plazo se especula una estabilización de su valor. Sin embargo,
luego del 2020 se pronostica que el precio del gas natural tenderá nuevamente a
recuperar su tendencia al alza.
El tipo de cambio del peso chileno con el dólar se ha mantenido en los últimos 10 años
oscilando dentro de un rango de los 600 a los 470 CLP/USD, alcanzando su máximo
valor el año 2014.
700 0,8200
0,8000
650 Tipo de Cambio (EUR/USD)
Tipo de Cambio (CLP/USD)
0,7800
0,7600
600
0,7400
550 0,7200
0,7000
500
0,6800
0,6600
450
0,6400
400 0,6200
2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014
30
Universidad Técnica Federico Santa María
Departamento de Industrias
El tipo de cambio del euro con el dólar se ha comportado de una manera más irregular,
situándose dentro de los 0,8 EUR/USD para el año 2004, para caer a 0,67 EUR/USD en
el año 2008. Sin embargo, desde ese año se ha presentado una tendencia al alza hasta la
fecha, recuperándose hasta alcanzar un valor de 0,78 EUR/USD el año 2014.
7.1.8. Resumen
Considerando todos estos antecedentes, se puede inferir que la industria minera hasta el
año 2020 experimentará un crecimiento importante, bajo un entorno económico estable.
Junto a esto y la continua disminución de la ley del cobre, las presiones al alza del
precio de los combustibles fósiles, la inminente expansión de la oferta del cobre de la
industria chilena (40% más al año 2020) y su inevitable incremento en la demanda de
energía (53% más para el año 2020), permiten que evaluar alternativas para mejorar el
rendimiento económico (o bien reducir los costos medios de operación) de los proyectos
y operaciones mineros sea atractivo para los inversionistas y/o el estado de Chile.
31
Universidad Técnica Federico Santa María
Departamento de Industrias
7.2. Colectores Solares en el mercado
Los colectores solares de placa plana son económicos, mecánicamente simples, tienen
un bajo costo de mantenimiento y sirven para colectar tanto como radiación directa
como radiación difusa.
Este tipo de colectores solares tienen la característica de que su área de apertura es muy
similar a su área bruta, por lo que las pérdidas de energía serán mayores a otros
colectores solares que presentan menos área expuesta. Por ello, la desventaja de los
colectores solares de placa plana es que a temperaturas de operación muy elevadas, su
eficacia tiende a ser baja.
32
Universidad Técnica Federico Santa María
Departamento de Industrias
7.2.2. Colectores solares de tubo de vacío (CTV)
Los colectores solares de tubo de vació consiste en varias filas de tubos situados en
paralelo, típicamente construidos de cobre y recubiertos con vidrios, dejando un espacio
con vacío (<10-2 Pa). El tubo de cobre que pasa por dentro transporta en su interior el
fluido a calentar, y tiene una forma de “U” recubierto en un material que transfiere el
calor. El vidrio exterior se pinta con una capa selectiva para mejorar la captura de la
radiación.
Figura 2 - Componentes básicos de un tubo de vacío típico – Imagen: Estado del arte colectores [5]
33
Universidad Técnica Federico Santa María
Departamento de Industrias
7.2.3. Colectores solares cilíndrico parabólicos (CCP)
Figura 3 - Esquema de un colector solar cilíndrico parabólico (CCP) – Fuente de la imagen: AEE INTEC,
Austria
Estos colectores solares al tener una superficie de absorbedor muy pequeña en relación a
su área de apertura, experimentan bajas pérdidas al operar a altas temperaturas,
permitiéndoles calentar fluidos sin problemas sobre 100 °C y llegar a temperaturas de
hasta 400 °C.
Esta tecnología cuenta con un sistema de seguimiento solar en un eje, por el cual el
colector solar será capaz de seguir el desplazamiento del sol de este a oeste. Por ello
deben contar con una estructura de soporte sofisticada y un flexible transportador del
fluido calo portador capaz de soportar el movimiento mecánico del colector solar.
34
Universidad Técnica Federico Santa María
Departamento de Industrias
7.2.4. Colector solar fresnel de concentración lineal (CCF)
Figura 4 - Esquema de un colector solar fresnel de concentración lineal (CCF) – Imagen: Estado del arte
colectores [5]
Estos colectores solares ocupan superficies muy grandes y están diseñados para usos en
parques de colectores de gran escala.
35
Universidad Técnica Federico Santa María
Departamento de Industrias
7.2.5. Colector de concentración con reflector estacionario (MAC)
Figura 5 - Esquema funcionamiento con reflector estacionario (MAC) – Imagen: Estado del arte colectores [5]
Entonces estos colectores solares se construyen situando una brazo con su eje situado en
el centro de la circunferencia imaginaria, el cual podrá rotar desplazando al tubo
absorbedor a lo largo del perímetro del círculo por el cual se va concentrando la
irradiación durante el transcurso del día (Figura 6).
36
Universidad Técnica Federico Santa María
Departamento de Industrias
Figura 6 - Esquema de un colector solar con reflector estacionario (MAC) – Imagen: Estado del arte colectores
[5]
37
Universidad Técnica Federico Santa María
Departamento de Industrias
7.3. Evolución de los sistemas de colectores solares
2.000
Colectores
1.500
1.000
500
0
1995 2005 2015 2025 2035
El norte de Chile presenta uno de los índices de irradiación en el plano horizontal mejor
evaluados del mundo, alcanzando en Calama un valor de 6,49 kWh/m2/d, con un índice
de transparencia atmosférica mayor a 0,7 el 92% de los días del año y un 1% menor a
0,6 [8]. La irradiación a lo largo del norte de Chile se puede observar a continuación
(Figura 7).
38
Universidad Técnica Federico Santa María
Departamento de Industrias
700
600
500
Chuquicamata
400 Radomiro Tomic
kTMF
300 Salvador
Andina
200
El Teniente
100
Gabriea Mistral
0
1993 1998 2003 2008 2013
Año
Gráfico 10 - Producción Anual de Cobre Fino por División de Codelco, Elaboración propia con datos de
CODELCO.
Radomiro Tomic se ubica a 3.000 metros sobre el nivel del mar, en la Región de
Antofagasta, en la provincia el Loa y a 45 km de Calama. Los yacimientos fueron
descubiertos el año 1950 y CODELCO comenzó las operaciones en la faena el año 1955,
dónde 3 años después se iniciaron las operaciones en régimen [6], con una capacidad
productiva nominal de 180.000 TMF. Hoy en día esta alcanza un valor igual 320.000
TMF y la producción anual de óxidos se ha mantenido en torno a los 300.000 TMF
desde el año 2002 y fijándose en 305.000 ton/año programadas para el año 2012 [7] [8].
Sin embargo, las proyecciones de producción de cátodos de cobre para el año 2013 son
de 240.000 TMF y estas continuarán disminuyendo paulatinamente hasta que el
procesamiento de óxidos deje de ser factible. Debido a esta abrupta disminución en la
producción de cobre fino, para fines del año 2016 se verá completada la construcción y
montaje del proyecto “Radomiro Tomic Sulfuros”, que significará la continuidad de la
40
Universidad Técnica Federico Santa María
Departamento de Industrias
explotación de la mina a rajo abierto de Radomiro Tomic mediante la explotación de sus
reservas de sulfuros, estimadas en unas aproximadas 2.6000 MT con una ley igual a
0,47% [9]. Mientras se construye la planta de sulfuros, el mineral sulfurado extraído es
enviado a procesar a la División Chuquicamata.
Figura 8 - Vista Satelital de la División Radomiro Tomic en Enero del 2014 (Extraída de Google Earth)
Para poder extraer el mineral enriquecido con cobre, primero se debe tronar la roca y
retirar el material estéril. Una vez que los bancos de mineral enriquecido con cobre del
rajo de la mina están expuestos, estos son extraídos y depositados mediante palas
eléctricas sobre camiones de alto tonelaje, que transportan el material hacia el proceso
41
Universidad Técnica Federico Santa María
Departamento de Industrias
de chancado. Acorde a la planificación de explotación del yacimiento del año 2012, se
extraen alrededor de 580 kton/día de material, con una relación media de estéril/mineral
igual a 2,7.
El mineral extraído del rajo (del cual un 2% proviene de material extraído de Mina
Sur) es depositado y procesado por un chancador giradorio Fuller de 60”x110” OSS 10”,
el cual busca reducir el diámetro de las rocas a menos de 10 pulgadas. Este equipo es
capaz de recibir una descarga simultánea de hasta 600 ton, con una capacidad de
procesado igual a 10 kton/h. La roca procesada es almacenada en un stock con una
capacidad de almacenaje de 80 kton. Este material es enviado mediante correas
transportadoras al proceso de chancado secundario.
42
Universidad Técnica Federico Santa María
Departamento de Industrias
7.7.2.3. Chancado Terciario
43
Universidad Técnica Federico Santa María
Departamento de Industrias
El mineral procesado del proceso de chancado terciario, es curado con una solución
concentrada de H2SO4 a una tasa de 2 a 7 kg ácido/T Mineral sobre correas transportadoras,
44
Universidad Técnica Federico Santa María
Departamento de Industrias
mediante un sistema de regadores dispuestos sobre la correa transportadora, para
disolver las especies mineralógicas contenedoras de cobre (atacamita, crisocola y
alterado de mina sur), al capturarlo mediante reacciones de óxido-reducción en forma de
CU2SO4, solución que es drenada y posteriormente recuperada. El proceso de lixiviación
se separa en las etapas de lixiviación primaria, secundaria (o de ripios) y de OBL, lsa
cuales se describen a continuación.
Lixiviación Primaria
45
Universidad Técnica Federico Santa María
Departamento de Industrias
Lixiviación Secundaria o Lixiviación de Ripios
Lixiviación de OBL
46
Universidad Técnica Federico Santa María
Departamento de Industrias
7.7.3.2. Planta de Extracción por Solventes (SX)
Los procesos de cada tren se puede dividir en tres etapas: Etapas de extracción o
carga del reactivo, otra de lavado y otras de re extracción o descarga del reactivo. Cada
una de las etapas cuenta con equipos mezcladores de fases que luego alimentan a un
estanque de decantación. Cada tren tiene una capacidad de diseño de 1.750 m3/h y
trabajan con una razón de orgánico/acuoso igual a 1,1. En los trenes A, B y C existen
dos etapas de extracción en serie, seguido de una etapa de lavado y luego una etapa de re
extracción. En el caso del tren D, existen dos etapas de extracción, una de lavado y dos
etapas de re extracción.
47
Universidad Técnica Federico Santa María
Departamento de Industrias
por sobre la acuosa debido a su menor peso específico. La fase acuosa resultante es
enviada a la segunda etapa de extracción, dónde pasa previamente por un mezclador que
recibe la solución orgánica descargada (con un contenido muy bajo de cobre). Las fases
se separan de la misma manera que en la primera etapa de extracción. La fase acuosa
resultante de esta etapa es enviada a la piscina de refino 1, mientras que en el caso de los
trenes A, B y C, la fase orgánica cargada resultante de la primera etapa de extracción es
enviada al patio de estanques, donde pasa por 3 coalescedores Chuquicamata en paralelo
con el fin de remover arrastres de electrolito. Luego pasa a un estanque de orgánico el
cual alimenta a la etapa de limpieza de los respectivos trenes. En el caso del tren D, la
fase orgánica es enviada de manera independiente al patio de estanques, donde se hace
pasar por 6 coalescedores y luego a un estanque de orgánico, para luego alimentar a la
etapa de limpieza del mismo tren.
48
Universidad Técnica Federico Santa María
Departamento de Industrias
7.7.3.3. Patio de Estanques
49
Universidad Técnica Federico Santa María
Departamento de Industrias
7.7.3.4. Planta de Electro Obtención o Electro Winning (EW)
Reacción Global:
La reacción global indica que se consume un mol de agua y de sulfato de cobre por cada
mol de cobre depositado en un cátodo, produciéndose un mol de ácido sulfúrico y medio
mol de oxígeno molecular.
50
Universidad Técnica Federico Santa María
Departamento de Industrias
efecto es minimizando mediante el uso, en la superficie de cada celda, de dos a tres
capas de esferas de polipropileno de 19 mm de diámetro en la y aceite de quillay para
disminuir la tensión superficial de la interface líquido-gas y así evitar la explosión
violenta de burbujas. La neblina ácida generada es controlada por un sistema de
ventilación forzada DESOM modificado. Esto se logra mediante sopladores situados
sobre la superficie de las celdas electrolíticas que redirigen el flujo de aire hacia el
centro de la nave de electro obtención, donde es evacuado hacia el exterior mediante
chimeneas situadas a lo largo del techo de la planta.
Las celdas están configuradas en dos circuitos, donde 148 de las celdas son para
limpieza (Celdas Scavenger) de la solución de electrolito y 852 celdas comerciales. Cada
celda cuenta con 61 ánodos y 60 cátodos, con una separación promedio entre cátodos
igual a 10 cm.
51
Universidad Técnica Federico Santa María
Departamento de Industrias
La sustancia de electrolito rico proveniente de los intercambiadores de calor llega
directamente al banco 1, en el cual se encuentran las Celdas Scavenger. En estas celdas
se busca atrapar todas las contaminaciones posibles que pudiese acarrear el electrolito
rico, de tal manera que el producto obtenido por estas celdas posiblemente no cumpla
con la calidad de cátodos de cobre grado A. El electrolito después de pasar por estas
celdas, es dirigido al estanque de recirculación.
El manejo de la cosecha de los electrodos es realizado por cuatro puentes grúa, que
operan automáticamente y se encuentran programados con tres lavadoras y máquinas
despegadoras de cátodos. Estas retiran 20 cátodos simultáneamente por celda y los
montan en las lavadoras, que consisten en cubas especiales donde se lavan con agua de
procesos, eliminando impurezas de la superficie de los cátodos. Luego son llevado a las
máquinas despegadoras de cátodos y luego éstos son pesados, corrugados, muestreados,
para finalmente elaborar paquetes contenedores de aproximadamente 26 placas de cobre,
que pesan entre un rango de 2 a 2,5 TMF. Luego estos son sometidos a un análisis
químico, para después ser cargados en camiones y transportados hasta el puerto de
Antofagasta y enviarlos a los clientes.
52
Universidad Técnica Federico Santa María
Departamento de Industrias
Es por ello que los procesos de interés para este trabajo corresponderán a los que
demandan energía directamente de un sistema de calentamiento externo.
En esta sección se presentan las demandas de energía históricas el año 2008 recopiladas
de la bibliografía para las etapas de lixiviación, extracción por solvente y electro
obtención. Se escogió dicho año ya que este es el único año del cual se logró conseguir
bases de datos específicas de los parámetros operacionales de los procesos. Las
demandas energéticas se separan según su fuente de obtención, que puede ser energía
eléctrica o calórica.
53
Universidad Técnica Federico Santa María
Departamento de Industrias
Para la generación de energía térmica Radomiro Tomic cuenta con 6 calentadores
Cleaver Brooks modelo CB-100-700-150, que el año 2008 combustionó mezclas de
aproximadamente un 96,6% de diesel y 3,6% gas natural. Sin embargo esta relación
cambia drásticamente el año 2009; siendo un 58,7% para el diesel y un 41,3% para el
gas natural respectivamente [16] y debido al elevado y creciente costo del diesel frente
al gas natural, es de esperar que en el futuro se priorice el uso de éste último.
54
Universidad Técnica Federico Santa María
Departamento de Industrias
El año 2008 la planta de hidrometalurgia de Radomiro Tomic presentó los siguientes
consumos energéticos (Tabla 1):
Tabla 1 - Consumo energético total por combustibles LX-SX-EW Radomiro Tomic (2008)
Según la información sobre control de demanda del SING, el año 2008 la División
Radomiro Tomic facturó un consumo eléctrico igual a 798,8 GWh [17]. Según las bases
de datos de [16] aproximadamente un 76% (610,3 GWh) de la energía eléctrica
corresponde a los procesos de hidrometalurgia a analizar.
Los antecedentes históricos reflejan que del total de energía eléctrica empleada durante
el año 2008 en procesos hidrometalúrgicos, solo un 6,69% (40.837 MWh) corresponde
al consumo energético realizado por operaciones en la etapa de lixiviación.
La energía consumida por combustible diesel equivale a un 0,02% (27 MWh) del
total de energía proveniente de la combustión de diesel.
55
Universidad Técnica Federico Santa María
Departamento de Industrias
8.2.3. Demanda energética de la etapa de Electro Obtención (EW)
56
Universidad Técnica Federico Santa María
Departamento de Industrias
8.3. Energía transportada y parámetros de los fluidos en cada punto de
interés
Se cuenta con una base de datos con parámetros operacionales de flujos desde la
lixiviación hasta la electroobtención. A continuación se determina el calor transportado
por cada fluido en cada punto presentado en la Figura 10, la cual se determinó en cada
caso utilizando la temperatura de referencia estándar de 25 °C.
Esta información es necesaria para construir más adelante los balances de energía para
los volúmenes de control que se quieran analizar.
57
Universidad Técnica Federico Santa María
Departamento de Industrias
Figura 10 – Flujos de los procesos Hidrometalúrgicos Considerados con sus Temperaturas y Caudales
volumétricos.
58
Universidad Técnica Federico Santa María
Departamento de Industrias
La energía transportada por cada fluido en cada punto “ i ” se calculó como:
( ) ̇
Tabla 3 – Energía transportada y parámetros medios de fluidos en cada punto – Basado en datos reales del
trabajo de Palacios, 2011
59
Universidad Técnica Federico Santa María
Departamento de Industrias
8.4. Balances
El volumen de control que globaliza los análisis de este estudio encierra a todos los
flujos y procesos involucrados con los intercambiadores de calor y la etapa de
electroobtención. Esto se debe a que estas son las dos etapas en las cuales se suministra
la mayor cantidad de energía en los procesos Hidrometalúrgicos; energía de origen
térmico en los intercambiadores de calor y energía de origen eléctrico en la etapa de
electroobtención.
60
Universidad Técnica Federico Santa María
Departamento de Industrias
8.4.1. Calderas: Calentadores Diesel/Gas Natural
Los calentadores analizados son 6 generadores marca Cleaver Brooks, modelo CB-
100-700-150, con una capacidad de 5,2*106 kcal/h cada uno. Generan calor al
combustionar Diesel, Gas o una combinación entre ambos. El quemador de aceite
funciona con aire atomizado a baja presión y el quemador de gas es de tipo orificio y no
requiere pre mezcla. La compuerta de aire de estos equipos es tipo rotatoria modulada
eléctricamente. Operan calentando agua desde una temperatura media de 54°C hasta una
61
Universidad Técnica Federico Santa María
Departamento de Industrias
temperatura media de 105 °C. El caudal medio total circulante es de 244 m3/h, operando
continuamente 4 calentadores con un caudal aproximado de 61 m3/h cada uno. Durante
el año 2008 los calentadores consumieron 157.789 MWh de Diesel y 5.733 MWh de
Gas Natural, por lo que los calentadores demandaron en total 163.522 MWh para fines
de calentar agua de los intercambiadores de calor.
∑ ∑
A. Consumo de Combustible
Por lo tanto, los consumos de combustibles estimados para el año 2008 serán:
62
Universidad Técnica Federico Santa María
Departamento de Industrias
B. Gases de combustión y aire consumido
Para simplificar este cálculo, se considerará como que todo mol de combustible es
quemado completamente, por lo que los gases de combustión estarán compuestos
principalmente por agua, dióxido de carbono y nitrógeno (proveniente del aire
consumido). Los pesos moleculares del diesel y gas natural son 170,1 g/mol y 16,95
g/mol respectivamente y sus estructuras moleculares promedio se muestran a
continuación:
Con esta información se puede determinar los moles de combustibles quemado y los
moles de oxígeno necesarios para llevar a cabo la reacción, y los moles de agua y
dióxido de carbono producidos por la combustión. Para ello se debe respetar la siguiente
ecuación de estequiometria:
63
Universidad Técnica Federico Santa María
Departamento de Industrias
Aire (Consumido) 6.926 329
N2 (Calentado) 5.408 257
Entradas Salidas
Flujoi ṁi (T/año) Flujoj ṁj (T/año)
Diesel 13.822 H2O (Combustión) 15.986
Gas Natural 557 CO2 (Combustión) 41.997
O2 (Combustión) 48.613 N2 (Calentado) 146.039
N2 (aire) 157.708 12 1.525.504
13 1.566.783 14 538.846
15 554.137
Total 2.341.621 2.276.050
M (T/año) -65.571
64
Universidad Técnica Federico Santa María
Departamento de Industrias
8.4.1.2. Balance de Energía
∑ ∑
Para el cálculo de los calores específicos de los gases de entrada y salida del calentador,
se consideró que estos ingresan a la temperatura ambiente media de 13 °C (286 K) y
salen a 500°C (773 K) [18]. Para ello se empleó el modelo propuesto por Gordon y
McBride mostrado a continuación:
Donde
( )
Cp (kJ/kg K)
Entrada (286 K) Salida (773 K)
H2O 1,850 2,123
CO2 0,833 1,160
N2 1,037 1,114
O2 0,913 1,048
65
Universidad Técnica Federico Santa María
Departamento de Industrias
Con esta información se determinaron las entradas y salidas de energía, junto con el
balance de energía, cuyos resultados se presentan en la Tabla 7:
Entradas Salidas
Flujoi Hi (MWh/año) Flujoj Hj (MWh/año)
Diesel 162.886 H2O (Combustión) 4.893
Gas Natural 5.733 CO2 (Combustión) 6.969
O2 (Combustión) -145 N2 (Calentado) 23.318
N2 (aire) -538 12 142.338
13 36.329 14 50.609
15 19.294
Total 223.558 228.128
Q (MWh/año) 4.569
66
Universidad Técnica Federico Santa María
Departamento de Industrias
8.4.2. Patio de estanques
67
Universidad Técnica Federico Santa María
Departamento de Industrias
8.4.2.1. Sistema de intercambiadores de calor.
A. Balance de Masa
La ecuación que describe el balance de masa para cada entrada y salida del
sistema de intercambiadores de calor, se muestra a continuación.
∑ ∑
68
Universidad Técnica Federico Santa María
Departamento de Industrias
Resolviendo la ecuación, las entradas y salidas de masa, junto con el balance de masa
resuelto se muestra a continuación en la Tabla 8:
Entradas Salidas
Flujoi ṁi (T/año) Flujoj ṁj (T/año)
4 2.937.120 6 2.901.600
10 2.866.107 11 2.866.107
12 172.530 13 178.254
Total 5.975.757 5.945.961
dm/dt (T/año) -29.796
Como se puede observar, el balance de masa refleja una pérdida de 29.796 T/año. Si
bien este valor es a simple vista bastante elevado, relativamente a los flujos de entrada
no lo es, ya que significaría una pérdida (o error en el balance) del 0,5%. Por lo tanto, se
considerará al balance de masa en equilibrio.
B. Balance de Energía
∑ ∑
Entradas Salidas
Flujoi Hi (MWh/año) Flujoj Hj (MWh/año)
4 11.588 6 412.885
10 577.941 11 271.784
12 142.338 13 36.329
Total 731.867 720.997
Q (MWh/año) -10.870
69
Universidad Técnica Federico Santa María
Departamento de Industrias
Como se puede observar, el balance de energía refleja una pérdida de energía
equivalente a -10.870 MWh/año. Esta pérdida equivale a 1,5% del flujo de energía de
entrada y pese a la correcta aislación térmica, puede deberse a pequeñas fugas naturales
por radiación y convección.
Figura 14- Volumen de Control: Intercambiador de Calor - Electrolito Rico y Electrolito Pobre
70
Universidad Técnica Federico Santa María
Departamento de Industrias
Entradas Salidas
Flujoi ṁi (T/año) Flujoj ṁj (T/año)
4 2.937.120 5 2.919.360
10 2.866.107 11 2.866.107
Total 5.803.227 5.785.467
M (T/año) -17.760
71
Universidad Técnica Federico Santa María
Departamento de Industrias
B. Balance de Energía - IC Electrolito Rico/Pobre
[ ] [ ] [ ] [ ] [ ] [ ]
[ ]
Entradas Salidas
Flujoi ṁi (T/año) Flujoj ṁj (T/año)
5 2.919.360 6 2.901.600
12 172.530 13 178.254
Total 3.091.890 3.079.854
M (T/año) -12.036
72
Universidad Técnica Federico Santa María
Departamento de Industrias
[ ] [ ] [ ] [ ] [ ] [ ]
Y resolviendo la ecuación:
[ ]
73
Universidad Técnica Federico Santa María
Departamento de Industrias
Nave de Electroobtención
74
Universidad Técnica Federico Santa María
Departamento de Industrias
existen otras dos salidas y una entrada para el volumen de control. Estas salidas son la
masa de cobre retirada del sistema por la cosecha de cátodos de cobre y el oxígeno (O2
gaseoso) liberado producto de la hidrólisis de agua. Por lo que para mantener el
equilibrio se considera un caudal másico de entrada de agua de reposición al sistema.
Debido a que el caudal volumétrico del electrolito rico es igual al del electrolito pobre
(2.400 m3/h) y existe un cambio en la composición química de los flujos de electrolito,
es de esperar que exista una variación en la densidad del fluido. La densidad del
electrolito pobre se determinó teóricamente al considerar la variación de su masa
producto de la reacción de electrodeposición.
∑ ̇ ∑ ̇
Dónde:
{ }
{ }
Se conoce que el año 2008 se produjo 240.954 T de cátodos comerciales de cobre puro,
por ende al considerar los pesos moleculares de cada componente químico involucrado y
respetando la estequiometria de la ecuación, se pueden determinar los moles y masas
para cada compuesto o elemento. Considerando que las 240.954 T de cobre corresponde
75
Universidad Técnica Federico Santa María
Departamento de Industrias
solo a la producción de cátodos comerciales, esto permitirá determinar cuántos moles de
agua y sulfato de cobre se consumieron, y cuántos moles de ácido sulfúrico y oxígeno se
produjeron durante la reacción, sin embargo, se desconoce cuántos de cada uno
estuvieron involucrados en la reacción de las celdas de limpieza.
Dónde:
( )
( )
[ ]
Dónde
76
Universidad Técnica Federico Santa María
Departamento de Industrias
( )
Por lo tanto, según esta ecuación la producción teórica de cobre en las celdas
comerciales sería:
[ ]
[ ] [ ] [ ]
[ ] [ ]
[ ] [ ] [ ]
Asumiendo que las celdas scavenger trabajan con la misma eficiencia de corriente, se
puede estimar la producción real de cátodos depositados (cátodos que no alcanzarían el
mismo nivel de pureza de los cátodos comerciales). Sin embargo para simplificar el
cálculo se considerarán como cátodos 100% puros de cobre, por lo que la masa teórica
de cobre depositada en los cátodos de limpieza equivale a:
[ ]
[ ] [ ] [ ]
[ ] [ ]
[ ] [ ] [ ]
[ ] [ ]
77
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo tanto la masa total de cobre depositado en el año 2008 (cátodos comerciales y de
limpieza) sería igual a:
Conociendo la masa total de cobre y su peso molecular, se obtienen los moles de cobre
involucrados en la reacción de óxido reducción, y realizando un balance estequiométrico
se obtienen los moles de los otros compuestos involucrados, según la ecuación
previamente citada:
Tabla 12 - Masa de los componentes químicos involucrados en la reacción de electro deposición durante el año
2008.
La masa de CuSO4 y H2O son disminuciones de la masa total del electrolito rico ya que
se consumen y el H2SO4 significa un incremento en su masa total. Sin embargo el Cu0
depositado es una disminución en la masa total del electrolito, ya que es retirado
mecánicamente del sistema durante la cosecha de cátodos. El O2 también es una
disminución, ya que al ser un gas es expulsado del sistema de forma natural.
78
Universidad Técnica Federico Santa María
Departamento de Industrias
C. Densidad Electrolito Pobre
Volumen (m3)
Cu0 Depositado 31.659
H2O Hidrolizado 80.227
O2 Liberado 58.782.561
H2SO4 Producido 237.544
CuSO4 Disociado 197.329
Por lo tanto, el volumen de agua de reposición para mantener el volumen del flujo será:
[ ] [ ]
̇ [ ] [ ] [ ] [ ]
∑ ̇ ̇ ̇ ̇ ̇
̇ ̇
79
Universidad Técnica Federico Santa María
Departamento de Industrias
∑ ̇ ̇ ̇ ̇ ̇ ̇
Dónde:
̇ ̇ ̇
̇ ̇
̇ ̇
Junto a esta información, ya se han determinado todos los flujos másicos ( ̇ ) menos el
flujo másico 10, porque aún se desconoce su densidad. Por lo tanto, asumiendo que el
balance de masa está en equilibrio se tendrá:
∑ ̇ ∑ ̇
Por lo que al reemplazar las sumatorias en la ecuación se podrá despejar el flujo másico
del electrolito pobre (10), resumiéndose a:
̇ ( ̇ ̇ ̇ ) ( ̇ ̇ )
Resolviendo:
̇ [ ]
̇ [ ]
̇ [ ]
[ ]
[ ] [ ]
80
Universidad Técnica Federico Santa María
Departamento de Industrias
D. Balance de masa
Entradas Salidas
Flujoi ṁi (T/año) Flujoj ṁj (T/año)
6 25.487.654 10 25.175.887
*H2O Reposición 156.572 Cu0 Cosechado 282.810
H2SO4 Producido 437.081 O2 Expulsado 71.313
O2 Producido 71.313 H2O Hidrolizado 80.227
0
Cu Depositado 282.810 CuSO4 Disociado 710.976
*H2O Evaporado 124.172
Total 26.430.995 26.430.995
M (T/año) 0
( ∫( ) )
81
Universidad Técnica Federico Santa María
Departamento de Industrias
̇
A. Electrolito Rico
[ ] [ ] [ ] [ ] [ ]
[ ]
B. Electrolito Pobre
[ ] [ ] [ ] [ ] [ ]
[ ]
[ ] [ ] [ ] [ ] [ ]
[ ]
Como se puede observar, existe una diferencia de 5.645 MWh al año, por lo que el
efecto por la variación de densidad puede ser despreciado.
[ ] [ ] [ ] [ ] [ ]
[ ]
82
Universidad Técnica Federico Santa María
Departamento de Industrias
D. Cambios de entalpía del electrolito producto de la reacción
de electroobtención
Por lo tanto el cobre, oxígeno y ácido sulfúrico generado son entalpías de entrada:
[ ] [ ] [ ] [ ] [ ]
[ ] [ ] [ ] [ ] [ ]
[ ] [ ] [ ] [ ] [ ]
[ ] [ ] [ ] [ ] [ ]
[ ] [ ] [ ] [ ] [ ]
La energía disipada por efecto Joule se produce principalmente en los contactos ánodo y
cátodo, electrodo anódico, electrolito y electrodo catódico [7]. Para determinar estas
pérdidas se debe conocer la caída de potencial operacional disipada por la resistencia a
83
Universidad Técnica Federico Santa María
Departamento de Industrias
la corriente eléctrica de cada componente ( ). Conociendo la intensidad de corriente ( )
se puede determinar el calor disipado mediante la siguiente ecuación:
( )
Donde
( )
Las entalpías de formación para cada compuesto se pueden observar en la Tabla 17:
ΔH0f (kJ/mol)
Cu0 0
H2O -285,84
O2 0
H2SO4 -810,94
CuSO4 -745,33
85
Universidad Técnica Federico Santa María
Departamento de Industrias
Reemplazando los valores en la ecuación
Por lo tanto
( )
Por lo tanto, el consumo energético por unidad de masa de cobre depositada será:
86
Universidad Técnica Federico Santa María
Departamento de Industrias
Las celdas de la nave de electroobtención están ubicadas en filas de a par por banco
y los bancos se distribuyen de a par en forma paralela. Por cómo están distribuidas las
celdas, se consideró que las pérdidas relevantes de energía a analizar son las pérdidas
tanto de convección como radiación desde los muros exteriores de las celdas, la
superficie de polipropileno y pérdidas por evaporación del electrolito. En la Figura 17 se
pueden observar los perfiles de temperatura de las celdas que se emplearán en los
análisis de convección y radiación, y el perfil de densidad de agua en el aire utilizado
para la determinación de las pérdidas por evaporación.
87
Universidad Técnica Federico Santa María
Departamento de Industrias
88
Universidad Técnica Federico Santa María
Departamento de Industrias
89
Universidad Técnica Federico Santa María
Departamento de Industrias
90
Universidad Técnica Federico Santa María
Departamento de Industrias
i. Pérdidas en la dimensión Y
a. Conducción y Convección
91
Universidad Técnica Federico Santa María
Departamento de Industrias
1) Superficie de la celda
92
Universidad Técnica Federico Santa María
Departamento de Industrias
Como no se cuenta con un modelo para determinar el número de Nusselt que
describa el comportamiento de la convección en la superficie de una celda de
electroobtención, se asumirá que el coeficiente de convección térmico promedio entre
los electrodos es similar al coeficiente de convección en todas las superficies del
electrolito, por lo tanto:
Dónde:
( )
( )
La velocidad vertical (dimensión “Y”) del electrolito que fluye en las celdas se
determinó como el caudal de electrolito que pasa por una celda dividido el área de la
sección perpendicular del flujo (plano “XZ”), es decir:
( )
93
Universidad Técnica Federico Santa María
Departamento de Industrias
El largo equivalente de los electrodos de una celda de electroobtención es:
El área de polipropileno corresponde al área superior total de la celda menos el área que
ocupan los electrodos en la superficie (0,4 m2 para los cátodos y 0,2 m2 para los ánodos).
Por lo tanto esta será:
En la nave de electroobtención las celdas comerciales operan con una temperatura media
de 48,5 °C, mientras que las celdas de limpieza operan a una temperatura media de 44
°C. Por lo tanto, se estima una temperatura media ponderada para todas las celdas, esta
es:
94
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo tanto
Lo que evidencia que los efectos de convección forzada son despreciables y predomina
un régimen de convección natural turbulento.
95
Universidad Técnica Federico Santa María
Departamento de Industrias
El ancho de la capa de polipropileno se determinó al considerar que según [7] estas se
sitúan de manera ordenada una sobre otra formando 3 capas de bolas de diámetro igual a
19 mm.
( )
Interface polipropileno/aire
Se consideró que la velocidad del aire que circula sobre las celdas es igual a 1 m/s [7].
Las propiedades termodinámicas del aire atmosférico (a una altura de 2.900 m y una
presión de 70.121 Pa) empleadas son:
96
Universidad Técnica Federico Santa María
Departamento de Industrias
ν (m2/s) 1,60E-05
-1
β (K ) 3,49E-03
2
DAB, Ts (m /s) 2,64E-05
Pr aire 0,71
Sc 0,61
( )
( )
( )
̅̅̅̅
̅̅̅̅
97
Universidad Técnica Federico Santa María
Departamento de Industrias
̅̅̅̅ ( )
Por lo tanto, el coeficiente de transferencia de calor del aire en la superficie de las celdas
de electroobtención es:
Y para el total de la Nave EW en el año 2008 (considerando las 1.000 celdas) es:
98
Universidad Técnica Federico Santa María
Departamento de Industrias
2) Fondo de la celda
Para este caso la velocidad del electrolito en el fondo de cada celda será:
( )
Por lo tanto
Como era de esperar, al igual que en el caso anterior, los efectos de convección forzada
pueden ser despreciados.
99
Universidad Técnica Federico Santa María
Departamento de Industrias
El régimen de convección natural sobre el fondo de la celda electrolítica se representó
mediante el modelo de convección natural de placa plana mostrado a continuación:
Por lo tanto, el coeficiente de convección de la cara interna del concreto del fondo de la
celda equivale a:
( )
100
Universidad Técnica Federico Santa María
Departamento de Industrias
Interface concreto polimérico/aire
Como tanto los costados como el fondo de las celdas electrolíticas se encuentran bajo el
nivel del suelo, se considera que no hay convección forzada alrededor de éstas.
El largo equivalente para la superficie exterior del fondo de las celdas equivale a:
Entonces, empleando las propiedades del aire previamente presentadas (Tabla 21), el
número de Grashof es:
( )
̅̅̅̅
Por lo tanto, el coeficiente de transferencia de calor del aire en el fondo de las celdas de
electroobtención es:
101
Universidad Técnica Federico Santa María
Departamento de Industrias
Y el Total Nave EW en el año 2008 (considerando las 1.000 celdas) es:
b. Radiación
El esquema que describe las pérdidas de energía debido a la radiación por la dimensión
Y se muestra a continuación en la Figura 22:
102
Universidad Técnica Federico Santa María
Departamento de Industrias
determinada por la temperatura superficial , el coeficiente de emisividad del
concreto polimérico y la temperatura cielo considerada en la nave de electroobtención
1) Superficie de polipropileno
103
Universidad Técnica Federico Santa María
Departamento de Industrias
previamente calculada ( ) la energía perdida por radiación por la
superficie de polipropileno será:
Por lo tanto la pérdida de energía por radiación de la Nave EW por la dimensión Y es:
c. Evaporación
104
Universidad Técnica Federico Santa María
Departamento de Industrias
Por analogía de capa límites, el número de Sherwood para convección forzada sobre una
superficie plana isotérmica equivale a:
̅̅̅
̅̅̅
̅̅̅ ( )
105
Universidad Técnica Federico Santa María
Departamento de Industrias
Las esferas de polipropileno sobre las celdas cumplen además la función de adicionar un
factor de recubrimiento a las celdas, evitando la evaporación de aproximadamente un
90% del agua desde la superficie de las celdas.
Las densidades del agua en el aire para una humedad relativa ambiental del 29% en
función de las temperaturas a considerar se muestran a continuación en la Tabla 23.
Ambiente Superficie
T (K) 286 298
ρwv (kg/m3) 0,247 0,818
̇ ( )
Y considerando el calor de vaporización del agua, la energía total perdida de las celdas
por evaporación el año 2008 fue:
106
Universidad Técnica Federico Santa María
Departamento de Industrias
de la energía calculada en esta sección está representada en rojo en la Figura 24 a
continuación:
a. Conducción y Convección
107
Universidad Técnica Federico Santa María
Departamento de Industrias
El flujo de calor está determinado por la temperatura interna del electrolito , el
coeficiente de convección en la superficie lateral ancha del electrolito , la
temperatura de la interface electrolito/muro del fondo , el coeficiente de
conductividad del concreto polimérico , la temperatura exterior del muro lateral ancho
, el coeficiente de convección del aire al fondo de las celdas y la temperatura
promedio del aire de la nave de electroobtención .
El largo equivalente para la cara interior del muro de la celda paralelo al plano ZY es:
La temperatura de la cara interna del concreto en el plano paralelo “ZY” se iteró hasta
obtener:
( )
Como tanto los costados como el fondo de las celdas electrolíticas se encuentran bajo el
nivel del suelo, se considera que no hay convección forzada alrededor de éstas.
Entonces, empleando las propiedades del aire previamente presentadas (Tabla 21), el
número de Grashof es:
109
Universidad Técnica Federico Santa María
Departamento de Industrias
( )
Por lo tanto, considerando el número de Nusselt para una superficie plana vertical
isotérmica para un régimen turbulento de convección natural, tendremos:
̅̅̅̅
Por lo tanto, el coeficiente de transferencia de calor del aire en el muro externo de las
celdas de electroobtención paralelo al plano ZY es:
b. Radiación
El esquema que describe las pérdidas de energía debido a la radiación por la dimensión
X se muestra a continuación en la Figura 26:
110
Universidad Técnica Federico Santa María
Departamento de Industrias
111
Universidad Técnica Federico Santa María
Departamento de Industrias
a. Conducción y Convección
112
Universidad Técnica Federico Santa María
Departamento de Industrias
La temperatura de la cara interna del concreto en el plano paralelo “XY” se iteró hasta
obtener:
113
Universidad Técnica Federico Santa María
Departamento de Industrias
Simplificando y asumiendo régimen de convección natural, el número de Rayleigh
queda:
( )
114
Universidad Técnica Federico Santa María
Departamento de Industrias
Como tanto los costados como el fondo de las celdas electrolíticas se encuentran bajo el
nivel del suelo, se considera que no hay convección forzada alrededor de éstas.
( )
Por lo tanto, considerando el número de Nusselt para una superficie plana vertical
isotérmica para un régimen turbulento de convección natural, tendremos:
̅̅̅̅
Por lo tanto, el coeficiente de transferencia de calor del aire en el muro externo de las
celdas de electroobtención paralelo al plano XY es:
Y el Total Nave EW en el año 2008, considerando que el área expuesta analizada solo
existe en 4 celdas por banco (recordar que cada banco tiene solo dos columnas de
celdas) , corresponde a:
115
Universidad Técnica Federico Santa María
Departamento de Industrias
Como se puede observar, al existir un área expuesta tan pequeña este análisis de pérdida
de energía podría ser perfectamente despreciado.
b. Radiación
El esquema que describe las pérdidas de energía debido a la radiación por la dimensión
Z se muestra a continuación en la Figura 29:
Convección y Conducción
Radiación
Evaporación
Por lo tanto, el total de pérdidas para las celdas de electroobtención será igual a:
117
Universidad Técnica Federico Santa María
Departamento de Industrias
G. Pérdidas de energía en Estanque de Recirculación
Todos los estanques están construidos con concreto de 0,35 m de espesor, recubiertos
internamente por una capa de aislante protectora de HDPE, la cual tiene 0,005 m de
espesor. Estos se ubican en el patio de estanques en la intemperie y la temperatura media
ambiental es de 14 °C.
Según la evaluación del recurso eólico a 5 metros de altura (Informe creado por Internet
por el departamento de Geofísica de la Universidad de Chile), la velocidad promedio del
viento en el patio de estanques es de 5,4 m/s, y como se observa en la rosa del viento en
la Figura 30, con una dirección de oeste a este.
Como se puede observar en la rosa del viento, este predomina de oeste a este siendo
prácticamente despreciable su magnitud en las otras direcciones. Por este motivo se
118
Universidad Técnica Federico Santa María
Departamento de Industrias
considerará que existe convección forzada en los lados de los muro por donde llega
directamente el viento, como se ve en la figura a continuación:
La manera en que están distribuidos estos estanques se puede observar en la Figura 32:
119
Universidad Técnica Federico Santa María
Departamento de Industrias
i. Pérdidas en la dimensión Y
a. Conducción y Convección
Interface electrolito/superficie
( )
Se realizó una iteración realizando los cálculos presentados a continuación hasta obtener
una temperatura de la superficie del electrolito ( ) que respete el flujo de energía
en todas partes..
Por lo tanto
121
Universidad Técnica Federico Santa María
Departamento de Industrias
Interface electrolito/aire
El área superficial y el perímetro en esta interface son iguales al anterior, por lo que el
largo equivalente es:
Por lo tanto, considerando el número de Nusselt para una superficie horizontal con
convección natural turbulenta:
̅̅̅̅
Por lo tanto, el coeficiente de transferencia de calor del aire en la superficie del estanque
de recirculación es:
122
Universidad Técnica Federico Santa María
Departamento de Industrias
b. Evaporación
̅̅̅
Las densidades del agua en el aire para una humedad relativa ambiental del 29% en
función de las temperaturas a considerar se muestran a continuación en la Tabla 25.
Ambiente Superficie
T (K) 286 320
ρwv (kg/m3) 0,247 0,749
̇ ( )
123
Universidad Técnica Federico Santa María
Departamento de Industrias
Y considerando el calor de vaporización del agua, la energía total perdida del estanque
de recirculación será:
c. Radiación
O sea
( )
Conociendo el factor de emisividad del electrolito (se consideró igual al del agua
), el área de emisión ( ) y su temperatura superficial
previamente calculada ( ) la energía perdida por radiación por la
superficie del estanque de recirculación será:
124
Universidad Técnica Federico Santa María
Departamento de Industrias
125
Universidad Técnica Federico Santa María
Departamento de Industrias
a. Conducción y Convección – Muro Oeste
La temperatura de la cara interna del HDPE en el plano paralelo “ZY” se iteró hasta
obtener:
126
Universidad Técnica Federico Santa María
Departamento de Industrias
( )
Para este caso se considera que fluye viento a una velocidad de 5,4 m/s. El largo
equivalente de esta interface es la misma que el anterior ( ).
Por lo tanto
127
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo tanto, considerando el número de Nusselt para una superficie plana vertical
isotérmica para un régimen turbulento de convección forzada, tendremos:
̅̅̅̅
Por lo tanto, el coeficiente de transferencia de calor del aire en el muro externo oeste del
estanque de recirculación es:
Por lo tanto, la pérdida estimada de energía por convección en el muro oeste es:
De la misma manera que en el caso anterior se determinó las pérdidas por radiación
conociendo el factor de emisividad del concreto ( ), el área de emisión
( ), su temperatura superficial previamente calculada (
) y la temperatura cielo previamente estimada ( ), la energía
perdida por radiación por la superficie del estanque de recirculación será:
128
Universidad Técnica Federico Santa María
Departamento de Industrias
Las pérdidas por convección por el muro sur ( ) dependen de la temperatura del
electrolito en el estanque de recirculación , del coeficiente de convección térmico
del electrolito en el muro sur , de los coeficiente de conductividad térmica del
HDPE y el concreto, de la temperatura del muro interno del estanque de alimentación de
electrolito , del coeficiente de convección térmico del electrolito del estanque
de alimentación y de la temperatura media de electrolito en el estanque de
alimentación .
Al igual que en el caso anterior, se desprendió que las pérdidas por el muro sur son
diferencias significativas en el muro norte. Por ende se considerarán como iguales las
pérdidas .
La temperatura de la cara interna Sur del HDPE en el plano paralelo “XY” se iteró hasta
obtener:
130
Universidad Técnica Federico Santa María
Departamento de Industrias
El régimen de convección natural la superficie interior del muro del estanque de
recirculación se representó mediante el modelo de convección natural turbulenta para
placa plana vertical, detallado a continuación:
Por lo tanto, el coeficiente de convección térmico del electrolito en el muro sur equivale
a:
( )
131
Universidad Técnica Federico Santa María
Departamento de Industrias
Interface HDPE / Electrolito de alimentación
( )
Por lo tanto
Por lo tanto, considerando el número de Nusselt para una superficie plana vertical
isotérmica para un régimen turbulento de convección natural, tendremos:
132
Universidad Técnica Federico Santa María
Departamento de Industrias
Convección y Conducción
Radiación
133
Universidad Técnica Federico Santa María
Departamento de Industrias
Evaporación
134
Universidad Técnica Federico Santa María
Departamento de Industrias
H. Balance de energía
Entradas Salidas
Flujoi Hi (MWh/año) Flujoj Hj (MWh/año)
6 412.885 10 577.941
Q Efecto Joule 275.079 CuSO4 + H2O + Cu0 + O2 4.920
H2O Reposición 6.768 Pérdidas Celdas EW 72.305
0
H2SO4 + Cu + O2 2.775 Pérdidas Estanque RC 47.529
Como se puede observar el balance de energía térmica refleja que las salidas se exceden
entrada, por lo que se considerará que este valor es un error aceptable y el sistema está
equilibrado.
135
Universidad Técnica Federico Santa María
Departamento de Industrias
Como se puede observar el balance de energía eléctrica refleja que las salidas se exceden
entrada, por lo que se considerará que este valor es un error aceptable y el sistema está
equilibrado.
136
Universidad Técnica Federico Santa María
Departamento de Industrias
9. Evaluación de Proyectos
En esta sección se describen y evalúan económicamente (mediante el software
RETSCreen) las alternativas de calentamiento de fluidos de procesos mediante distintas
configuraciones de sistemas de calentamiento solar de agua, variando desde las
temperaturas y flujos másicos en los circuitos de los colectores solares y fluido de
almacenaje, hasta las tecnologías de colectores solares empleadas.
137
Universidad Técnica Federico Santa María
Departamento de Industrias
9.1. Estructura base de Evaluación de Proyectos de Calentamiento Solar
de Agua
La IREC (Catalonia Institute for Energy Research) el año 2012 realizó una investigación
de mercado para colectores solares para aplicaciones de gran escala, en la cual presenta
una base de datos conteniendo 25 modelos de colectores solares característicos del
mercado europeo, con sus especificaciones técnicas y sus respectivos costos de inversión
por proveedor al año 2012 [19].
Por lo tanto, se cuenta con todos los parámetros técnicos (coeficientes de eficiencia
, coeficientes de pérdida de calor de primer orden ó , coeficientes
de pérdida de calor de segundo orden ó , áreas brutas del colector solar y áreas
de apertura del colector solar).
A continuación se muestra una tabla con los modelos de los colectores solares a estudiar
junto a su tecnología y especificaciones técnicas (Tabla 28).
138
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 28 - Información Técnica de Colectores Solares Evaluados (Según Área de Apertura)
Modelo Tecnología η0 UL1 (W/m2 K) UL2 (W/m2 K2) AApertura (m2) ABruta (m2)
HT-A 28/10 CPP 85,2% 3,481 1,30E-02 12,52 14
KAIROS XP 2,5 V CPP 80,8% 3,131 1,60E-02 2,26 3
SOL250 CPP 81,0% 3,480 2,00E-02 2,37 3
FKC-2S CPP 76,6% 3,216 1,50E-02 2,25 2
XL USB CPP 70,0% 3,560 1,00E-02 2,65 3
gigaSol S CPP 73,5% 3,653 1,10E-02 11,06 12
SKS 4,0 CPP 85,1% 4,036 1,08E-02 2,09 2
GS CPP 79,0% 3,979 1,40E-02 13,10 14
HT CPP 81,1% 2,710 1,00E-02 13,10 14
AL10 CPP 78,9% 3,834 1,10E-02 9,42 10
FA(3/6) CPP 74,8% 2,869 1,20E-02 16,58 20
VKF 125 CPP 75,3% 3,936 1,70E-02 2,35 3
VKF 145V CPP 79,0% 2,414 4,90E-02 2,35 3
XL10 CPP 78,9% 3,834 1,10E-02 9,43 10
10 HTF CPP 83,0% 3,523 1,50E-02 10,10 11
GJ140A CPP 80,3% 3,430 1,45E-02 13,9 15
DF100 CTV 77,3% 1,430 6,00E-03 3,23 4
HP250 CTV 76,1% 1,360 7,40E-03 3,23 4
DF30 - 300 CTV 77,6% 1,194 6,40E-03 3,05 5
HP30-300 CTV 70,3% 2,224 5,00E-03 3,05 5
APSE-30 CTV 68,7% 1,505 1,10E-02 2,98 4
CPC 45 CTV 64,4% 0,749 5,00E-04 4,50 5
X10T CCP 54,8% 0,849 3,00E-03 10,37 11
Polytrough 1800 CCP 68,5% 0,300 1,00E-03 36,90 37
CCStaR 2V MAC 71,8% 0,181 1,56E-03 37,10 44
LF-11 CCF 63,5% 0,000 4,32E-04 22,00 30
RETScreen realiza todos sus cálculos considerando que los coeficientes de la función de
eficiencia del colector solar a evaluar se encuentran en base al área bruta, sin embargo
los coeficientes de la función de eficiencia de los colectes previamente presentados se
encuentran en base a su respectiva área de apertura. Para ingresar estos parámetros a
RETScreen es necesario adaptar los coeficientes ( , , ) presentados en esta tabla
(Tabla 28) a su área bruta. En la práctica el área física que ocupa un colector solar será
mayor que el área bruta, ya que se deben ubicar dejando una distancia entre una fila de
colectores paralelas con otra, para evitar la generación de sombras, y según las
139
Universidad Técnica Federico Santa María
Departamento de Industrias
características físicas y técnicas del colector puede requerir espacio adicional. Para ello
se cuenta con un factor de espaciamiento ( ), el cual fue determinado para cada
uno de estos colectores solares por la IREC de la siguiente manera:
Dónde
Factor que considera la distancia mínima requerida entre cada colector solar.
Los factores de espaciamiento y las áreas brutas corregidas a considerar en este estudio
se presentan a continuación (Tabla 29):
140
Universidad Técnica Federico Santa María
Departamento de Industrias
VKF 125 CPP 1,5 3,77
VKF 145V CPP 1,5 3,77
XL10 CPP 1,2 12,06
10 HTF CPP 1,3 14,89
GJ140A CPP 1,5 22,53
DF100 CTV 1,4 5,95
HP250 CTV 1,4 5,96
DF30 - 300 CTV 1,3 6,10
HP30-300 CTV 1,3 6,30
APSE-30 CTV 1,4 5,90
CPC 45 CTV 1,3 6,41
X10T CCP 2,6 28,40
Polytrough 1800 CCP 2,6 95,90
CCStaR 2V MAC 2,6 113,60
LF-11 CCF 2,6 79,20
Por lo tanto, adaptando los coeficientes para el software RETScreen, cada coeficiente se
multiplicó por . Los valores de los coeficientes utilizados durante las
evaluaciones se muestran a continuación (Tabla 30):
Tabla 30 - Coeficientes de la Función de Eficiencia de los Colectores Solares (base área bruta)
141
Universidad Técnica Federico Santa María
Departamento de Industrias
HP250 CTV 41,2% 0,73655 4,01E-03
DF30 - 300 CTV 38,8% 0,59700 3,20E-03
HP30-300 CTV 34,0% 1,07705 2,42E-03
APSE-30 CTV 34,7% 0,76015 5,56E-03
CPC 45 CTV 45,2% 0,52590 3,51E-04
X10T CCP 20,0% 0,31000 1,10E-03
Polytrough 1800 CCP 26,4% 0,11543 3,85E-04
CCStaR 2V MAC 23,4% 0,05911 5,09E-04
LF-11 CCF 17,6% - 1,20E-04
Todos los colectores solares de tecnología de concentración (CCP, MAC y CCF) tienen
un sistema de seguimiento axial de oeste a este del sol, mientras los demás colectores
son fijos.
Considerando una irradiación solar en Calama sobre la superficie horizontal igual a 1065
W/m2 ( ), una temperatura del aire promedio igual a 13,2 °C ( ) y una temperatura de
entrada del agua a calentar ( ), las curvas de eficiencia de todos estos modelos se
pueden observar en el Gráfico 11:
142
Universidad Técnica Federico Santa María
Departamento de Industrias
40%
30%
20%
10%
0%
0,00 0,02 0,04 0,06 0,08 0,10 0,12 0,14
Te-Ta/Gh (°C m2 / W)
Como se puede observar, todos los colectores solares disminuyen su eficiencia a medida
que aumenta, y si se considera una temperatura ambiental ( ) y radiación
143
Universidad Técnica Federico Santa María
Departamento de Industrias
40%
30%
20%
10%
0%
0 20 40 60 80 100 120 140 160 180
Te (°C)
La selección del colector solar óptimo no es trivial; no basta con escoger simplemente el
colector solar con mayor eficiencia, ni tampoco el con el menor costo. Este análisis
considera una radiación solar constante y emula el mejor escenario. En Calama un
sistema de colectores solares fijo recibe una irradiación promedio anual de 282 W/m2,
mientras que los colectores solares con un sistema de seguimiento de oeste a este pueden
recibir una irradiación promedio anual de 375 W/m2 (pueden recibir casi un 33% más de
energía solar que un sistema fijo). Otro factor que puede ser decisivo es que operar con
144
Universidad Técnica Federico Santa María
Departamento de Industrias
temperaturas muy bajas puede significar la necesidad de operar con caudales muy
elevados, incrementando el costo por el dimensionamiento de cañerías y bombas frente a
un sistema que opere con temperaturas más elevadas y un caudal menor.
145
Universidad Técnica Federico Santa María
Departamento de Industrias
9.3. Selección de Procesos de Análisis
Los dos primeros procesos mencionados, reciben el calor entregado por flujos de agua
de calderas al pasar cada uno por un intercambiador de calor propio. El diagrama que
resume la alimentación de energía actual de estos dos procesos por el sistema de agua
caliente de calderas se puede ver a continuación, en la Figura 35.
146
Universidad Técnica Federico Santa María
Departamento de Industrias
En este estudio se evaluó la incorporación de energía en cualquiera de estos tres
procesos, ya que dependiendo de las temperaturas y caudales de operación se pueden
obtener rendimientos muy distintos para un mismo colector solar, y por lo tanto, generar
la diferencia entre que un proyecto evaluado sea o no rentable.
El diseño básico del sistema de calentamiento solar de agua para todos los proyectos
evaluados considera cinco circuitos:
147
Universidad Técnica Federico Santa María
Departamento de Industrias
Luego esta mezcla se hace pasar por un intercambiador de calor, donde el calor atrapado
por los colectores solares es entregado al fluido del circuito de carga (F), el cual eleva su
temperatura desde hasta .
148
Universidad Técnica Federico Santa María
Departamento de Industrias
Si la temperatura del estanque de reservorio ( ) es menor a la de carga ( ), el sistema
automático alimentará el estanque de reservorio a un flujo másico ̇ con una
temperatura hasta que la temperatura de la superficie del estanque de reservorio ( )
se iguale a , y la del fondo ( ) se iguale a . Si se dan esas condiciones, el estanque
de reservorio se considerará cargado a su máxima capacidad ( ).
149
Universidad Técnica Federico Santa María
Departamento de Industrias
adiciona calor auxiliar ( ) proveniente de un sistema de calderas en la medida justa
y necesaria para que el fluido de proceso alcance una temperatura igual a .
150
Universidad Técnica Federico Santa María
Departamento de Industrias
9.4.2. Alternativas de los sistemas evaluados
Como los fluidos de procesos cambiarán sus parámetros de operación (temperaturas y/o
caudales) dependiendo de dónde se analicen, y para tener un espectro más amplio sobre
como rinden económicamente los mismos colectores solares en diferentes condiciones,
se consideraron diferentes alternativas para satisfacer la misma demanda de energía.
En este caso el fluido de proceso a calentar es electrolito rico, el cual requiere que su
temperatura media se eleve desde 37,3 °C ( ) hasta 41,5 °C ( ). Lo
que significa que con un flujo medio de 2.400 m3/h ( ̇ ), una
3
densidad de 1.209 kg/m y un calor específico de 3,56 kJ/kg °C, requiere la continua
adición de 100.575 MWh al año ( ).
Para ello se diseñó un sistema de calentamiento solar con el cual se pretende pre calentar
el electrolito rico en un intercambiador de calor (IC – E.Rico / Descarga) que se
instalaría previo a el ingreso del electrolito rico al intercambiador de calor de las
calderas, de manera que este último opere como sistema de calentamiento auxiliar.
151
Universidad Técnica Federico Santa María
Departamento de Industrias
Figura 37 - Alternativa: Sistema de calentamiento solar de agua para pre calentar elelectrolito rico
En este caso el fluido de proceso a calentar es agua de osmosis, la cual requiere que su
temperatura media se eleve desde 9,0 °C ( ) hasta 60,0 °C ( ). Lo
que significa que con un flujo de 60 m3/h ( ̇ ), una densidad de
1.003 kg/m3 y un calor específico de 4,21 kJ/kg °C, requiere la continua la continua
adición de 31.256 MWh al año ( ).
Para ello se diseñaron dos alternativas de para el calentamiento solar de agua para
precalentar el agua de osmosis previo a su ingreso al sistema de calentamiento auxiliar
(intercambiador de calor agua de osmosis / agua de caldera), los cuales se presentan a
continuación junto a un diagrama que describe la alternativa y sus respectivos
parámetros de operación a evaluar.
152
Universidad Técnica Federico Santa María
Departamento de Industrias
Aportar energía solar de manera al agua de osmosis mediante intercambiador de
calor de agua de osmosis y agua de reserva calentada por energía solar (Figura
38).
Figura 38 - Alternativa: Sistema de calentamiento solar de agua para pre calentar agua de osmosis
153
Universidad Técnica Federico Santa María
Departamento de Industrias
Figura 39 - Alternativa: Sistema de calentamiento solar de agua modificado para pre calentar agua de osmosis
154
Universidad Técnica Federico Santa María
Departamento de Industrias
específico de 4,226 kJ/kg °C, para lograr el incremento de temperatura deseado
significaría la continua adición de 137.324 MWh al año ( ).
Para ello se diseñó un sistema de calentamiento solar de agua para pre calentar el agua
de caldera previo a su ingreso a las calderas, para lo cual se evalúa la incorporación de
un intercambiador de calor que le entregará al agua de caldera el calor capturado por el
sistema de calentamiento solar de agua. Esto junto a los parámetros operacionales a
evaluar se puede observar a continuación (Figura 40).
Figura 40 - Alternativa: Sistema de calentamiento solar de agua para pre calentar agua de calderas
155
Universidad Técnica Federico Santa María
Departamento de Industrias
9.4.3. Parámetros de Operación de los Proyectos
Para el diseño de equipos, los parámetros de operación del loop solar se determinaron en
base a su potencia nominal de captura solar. La potencia nominal de captura solar
( ) se puede determinar mediante la función de eficiencia del colector solar
evaluado, el área bruta total del campo solar, la temperatura de entrada al parque solar y
la radiación máxima incidente sobre la superficie de los colectores solares, o sea:
Donde
( )
( )
( )
156
Universidad Técnica Federico Santa María
Departamento de Industrias
minera Gaby. Este intercambiador de calor opera con las siguientes diferencias de
temperaturas en el lado de flujos calientes ( ) y el lado de flujos fríos ( ):
Y despejando tendremos:
157
Universidad Técnica Federico Santa María
Departamento de Industrias
9.4.3.2. Circuito de Carga
Para los diseños, se considerará que la temperatura de operación del fluido frío de carga
es la misma que la temperatura de operación del fluido de descarga, y que a la vez esta
es igual a la temperatura fría ( ) del fluido en el fondo del reservorio, o sea:
Donde
Para lograr esto, el flujo másico nominal de este circuito deberá ajustarse como:
158
Universidad Técnica Federico Santa María
Departamento de Industrias
Si bien el aporte solar medio ( ) será menor que el máximo ( ), los equipos
deben ser dimensionados para ser capaces de operar bajo las condiciones extremas, o sea
debe conocerse el flujo de masa para dicha condición ( ̇ ). Por lo tanto, este
corresponderá a:
Como la temperatura caliente del fluido del circuito de descarga ( ) siempre tendrá una
temperatura mayor a la temperatura del fluido de proceso pre calentado ( ), la
temperatura del fluido de descarga (considerando un mínimo igual a 3
°C) será al menos , pero si la diferencia entre las temperaturas del
fluido de proceso es muy pequeña, el sistema se forzará a mantener un flujo másico muy
elevado. Por lo tanto, al temperatura se ajustará como:
Donde K es una variable dicotómica multiplicada por un factor de 20 °C, cuya función
es activar una diferencia mínima de 20 °C entre y de la siguiente manera:
Esto es para que en todo caso se respete una holgura entre la temperatura caliente y fría
de al menos 20 °C en el fluido de reserva, evitando forzar la operación con caudales
másicos muy elevados.
159
Universidad Técnica Federico Santa María
Departamento de Industrias
Y corresponderá a incrementos realizados durante la optimización, para
determinar la temperatura de operación que permita maximizar los retornos del
proyecto.
Por otra parte, elevar la temperatura caliente de los circuitos y manteniendo constante la
temperatura fría, producirá un aumento en la diferencia de temperaturas, por lo que para
que el balance de energía se encuentre en equilibrio, significará una disminución de los
flujos másico. Una disminución de los flujos másicos permite operar con bombas más
chicas y de un costo reducido.
160
Universidad Técnica Federico Santa María
Departamento de Industrias
50% 1,10
45% 1,05
35% 0,95
30% 0,90
25% 0,85
20% 0,80
60 70 80 90 100 110 120
Temperatur de Salida de Colectores (°C)
Gráfico 13 - Fracción Solar y Variación del Costo de los Intercambiadores de Calor vs Temperatura de entrada
del Parque Solar
161
Universidad Técnica Federico Santa María
Departamento de Industrias
9.4.3.4. Circuito de Calderas
Dónde:
162
Universidad Técnica Federico Santa María
Departamento de Industrias
86,3%
y = 0,064x3 - 0,152x2 + 0,106x + 0,839
86,2% R² = 1
86,1%
Eficiencia (%)
86,0%
85,9%
85,8%
85,7%
85,6%
0% 20% 40% 60% 80% 100%
Carga (%)
89,4%
89,3%
89,3%
89,2%
89,2%
89,1%
89,1%
0% 20% 40% 60% 80% 100%
Carga (%)
De los ajustes, se calculó que la máxima eficiencia alcanzada por las calderas se alcanza
a un 51,8% de carga, con un 86,2% para Gas Natural y con un 89,5% para Diesel N° 2.
La potencia máxima de una de estas calderas es de 6.048 kW, por lo que se puede
determinar la potencia a la cual debe operar cada caldera para maximizar su eficiencia
para gas natural ( ) y para, esta será:
163
Universidad Técnica Federico Santa María
Departamento de Industrias
Sin embargo, para determinar la eficiencia estacional considerar que las calderas
siempre operan en su carga óptima sería un error que sobrevaloraría la capacidad del
sistema de calderas. Por ello se determinó un factor de ajuste de la eficiencia empleando
los datos reales de eficiencia estacional del caso base y relacionándolos con la eficiencia
teórica.
Si durante el año 2008 las calderas combustionaron 137.324 MWh, la carga del sistema
sería:
Y su carga porcentual:
164
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo tanto, la eficiencia instantánea del sistema de calderas será (para cada caso de
combustibles), corresponderá a la suma ponderada de la eficiencia instantánea de cada
caldera incluida en el sistema, o sea:
Entonces, como la eficiencia estacional del gas natural para el año 2008 fue igual a
83,7% y la de diesel N° 2 fue igual a 86,3%, los factores de ajuste serán:
165
Universidad Técnica Federico Santa María
Departamento de Industrias
Con esta información, se construyó un modelo para cada caso de combustible que
describe la eficiencia instantánea y el número de calderas operando dependiendo del
aporte de energía solar relativo al consumo base de combustible ( ).
Cabe destacar que cada modelo se corrigió con el factor de ajuste, transformando la
eficiencia instantánea a un estimado de eficiencia estacional. Estos modelos graficados
se pueden observar a continuación (Gráfico 16 y Gráfico 17).
86,0% 5
N° de calderas Activas
85,5% 4
85,0% 3
84,5% 2
84,0% 1
83,5% 0
0% 20% 40% 60% 80% 100%
Fracción Solar Total (%)
Gráfico 16 - Eficiencia estacional del sistema de calderas según fracción solar - Gas Natural
89,5% 5
N° de calderas Activas
89,0% 4
88,5% 3
88,0% 2
87,5% 1
87,0% 0
0% 20% 40% 60% 80% 100%
Fracción Solar Total (%)
Gráfico 17- Eficiencia estacional del sistema de calderas según fracción solar - Diesel N° 2
166
Universidad Técnica Federico Santa María
Departamento de Industrias
Considerando que durante 10 horas del día el sistema de calentamiento solar con
almacenaje entregue un 20% más de calor al fluido de proceso que la media diaria
( ), la eficiencia de las calderas media diaria ponderada variaría para un
combustible ( ) de la siguiente manera:
( ) ( )
Si bien el glicol puede operar desde un rango de los -28°C hasta los 188 °C, el fluido
para almacenar calor en el estanque de reservorio es agua, la cual se considerará que
puede ser calentada dentro de las iteraciones hasta un máximo de 100 °C y la presión de
vapor debe ser tal que el cambio de fase del agua se encuentre al menos 10 °C sobre la
temperatura caliente del agua ( ). La caída de presión necesaria para que el agua
experimente un cambio de fase a una temperatura determinada dada, según la ecuación
de Antonie, 1988:
( )
Donde
167
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo tanto y considerando los 10 °C de seguridad sobre la temperatura de evaporación,
la temperatura máxima permitida para el agua de reserva será:
L (m) 1,473
W (m) 0,908
s (m) 0,003
Dónde
168
Universidad Técnica Federico Santa María
Departamento de Industrias
Los coeficientes de transferencia de calor considerados empleados, corresponden a
valores medios de pares de fluidos circulantes en intercambiadores de calor de placa
típicos [21], cuyas magnitudes se pueden ver en la Tabla 32.
( )
Donde
169
Universidad Técnica Federico Santa María
Departamento de Industrias
: Coeficiente de transferencia de calor global del intercambiador de calor a diseñar
Y con
( )
( )
9.5.2. Cañerías
El diámetro de las cañerías usadas se basó en el diámetro estándar para cañerías de acero
inoxidable 316L Sch.40 de 8 pulgadas de diámetro (0,2032 m), diseñada para presiones
de trabajo de 2.362 kPa [22].
La longitud de la cañería del circuito de carga se consideró igual a 430 m, mientras que
la cañería del circuito de descarga se consideró igual a 85 m. En ambos casos se
requieren dos cañerías para llevar a cabo la recirculación de los fluidos de los circuitos.
170
Universidad Técnica Federico Santa María
Departamento de Industrias
9.5.3. Bombas
171
Universidad Técnica Federico Santa María
Departamento de Industrias
Donde es el área de apertura total del parque solar evaluado. Dentro de este
valor de caída de presión, también se considera la caída de presión del intercambiador de
calor, ya que la relación obtenida del campo solar de la minera Gaby considera todo el
circuito de colectores solares.
( )
( )
( )
Por lo tanto, las pérdidas de la cañería (carga o descarga según su caso respectivo) serán:
Donde
: Aceleración de gravedad ( )
172
Universidad Técnica Federico Santa María
Departamento de Industrias
Caída de presión en Intercambiadores de Calor
Con
Y dónde
173
Universidad Técnica Federico Santa María
Departamento de Industrias
: Largo de las placas del intercambiador de calor (m)
La potencia nominal de las bombas será la requerida para vencer la caída del fluido en el
circuito operando bajo su caudal nominal.
174
Universidad Técnica Federico Santa María
Departamento de Industrias
La caída de presión total a considerar en las bombas ( ) corresponderá a la
caída de presión de los circuitos ( ) más la presión requerida para evitar la
evaporación del fluido ( ), considerando un factor de seguridad de 10 °C por sobre
la temperatura máxima de operación, por lo tanto esta será:
Con
( )
Donde
Por lo tanto la potencia nominal requerida por el sistema de bombas de cada circuito
será:
175
Universidad Técnica Federico Santa María
Departamento de Industrias
( )
176
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6. Evaluación de proyectos con RETScreen
Para cada modelo de colector solar, se iteraron 10 valores de temperatura del circuito de
descarga ( ), y dentro de cada una de esas iteraciones, se realizó una iteración de 20
valores de número de colectores a considerar ( ) en el parque solar.
177
Universidad Técnica Federico Santa María
Departamento de Industrias
Mediante la programación de macros con Visual Basic en Microsoft Excel y en
combinación con el software RETScreen, se construyó un programa para evaluar los
proyectos, cuyo algoritmo de optimización se muestra a continuación (Figura 43).
178
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.1. Adaptación al Circuito de agua de RETScreen
179
Universidad Técnica Federico Santa María
Departamento de Industrias
Donde el calor total demandado por el fluido de proceso se puede determinar en función
del flujo másico del fluido de proceso ( ̇ ), su calor específico ( ), su
temperatura fría inicial ( ) y su temperatura de operación deseada ( ), como:
Vale considerar que el caudal evaluado por RETScreen irá disminuyendo a medida que
el valor de aumenta durante las iteraciones, de tal manera de respetar la misma
cantidad de energía.
Una vez que RETScreen determina el área bruta del parque solar, entrega el valor de la
fracción solar del sistema. Es decir, la magnitud relativa de energía que aporta el parque
solar al circuito de agua. O sea:
180
Universidad Técnica Federico Santa María
Departamento de Industrias
THOT (°C) 105,0 41,5 60,0
TCOLD (°C) 47,6 37,3 8,5
ṁProceso (kg/h) 232.046 2.337.335 59.558
Una vez ingresada toda la información sobre la ubicación se deben ingresar los datos
para que RETScreen alimente el modelo de energía que optimiza el tamaño del parque
solar, esta información se ingresa en la pestaña “Modelo de Energía”.
181
Universidad Técnica Federico Santa María
Departamento de Industrias
A. Temperatura óptima de calentamiento
Con
Con
182
Universidad Técnica Federico Santa María
Departamento de Industrias
El caudal que se debe configurar en RETScreen para el ajuste de la iteración
( ), deberá respetar el balance de energía del calor total requerido por el
fluido de proceso, o sea:
En este ítem se deben incluir la información sobre la orientación física de los colectores
solares. RETScreen calcula automáticamente la radiación solar incidente sobre la
superficie de estos, dependiendo de su inclinación, azimut y sistema de rastreo solar.
Los colectores solares de placa plana (CPP) y de tubo de vacío (CTV) considerados en
este estudio no tienen sistema de rastreo solar, por lo que se considerarán fijos.
183
Universidad Técnica Federico Santa María
Departamento de Industrias
2,50
1,50
2,00-2,50
1,00
1,50-2,00
0,50 1,00-1,50
0,000 0,50-1,00
103 0,00-0,50
206
Azimut °
309 310
210
110
10 Inclinación °
Para determinar la inclinación y azimut óptimos de los colectores solares con sistema de
rastreo uniaxial, se hicieron variar hasta determinar la combinación óptima que
maximice la radiación solar anual en el plano de los colectores (Gráfico 19)
184
Universidad Técnica Federico Santa María
Departamento de Industrias
.
3,5 3-3,5
Radiación (MWh/año m2) 2,5-3
3
2,5 2-2,5
2 1,5-2
1,5 1-1,5
1 0,5-1
0,5 0-0,5
0
0 310
103 210
206 Inclinación °
Azimut ° 110
309
10
Gráfico 19 - Radiación solar en el plano del Colector Solar – Sistema de Rastreo Uniaxial
185
Universidad Técnica Federico Santa María
Departamento de Industrias
Los datos de radiación solar diaria media para cada mes ( ) considerados por
RETScreen para todas las evaluaciones, según tipo de rastreo solar, se muestran a
continuación en el Gráfico 20:
12,00
11,00
10,00
9,00
kWh/m2 día
8,00
Horizontal
7,00
6,00 Fijado
5,00 Uniaxial
4,00
3,00
Gráfico 20 - Radiación Solar en Calama: Horizontal e Inclinada (Superficie de los Colectores Solares)
La radiación solar diaria promedio anual en el plano horizontal sería de 6,49 kWh/m2/d,
sobre la superficie de los colectores solares fijos sería de 6,80 kWh/m2/d y sobre la
superficie de los colectores solares con sistema de seguimiento uniaxial sería de 9,03
kWh/m2/d.
En este ítem deben incluirse los parámetros técnicos del colector solar a evaluar en el
proyecto. Se puede incluir solo un modelo de colector solar por evaluación a la vez.
Coeficiente Fr UL: ( )
186
Universidad Técnica Federico Santa María
Departamento de Industrias
Número de colectores:
Pérdidas varias: 3,5%
o Se utilizó la recomendación del panel de ayuda de RETScreen, que dice
que este valor puede oscilar de 2 a 5% para colectores con un buen
mantenimiento. Por lo cual se considerará el promedio de este rango.
El número de colectores óptimos se determinó para cada colector iterando su valor con
el fin de maximizar el VPN (Valor Presente Neto) del proyecto. La iteración se realizó
usando como primer valor un 10% del tamaño óptimo de colectores sugerido por
RETScreen ( ) y se continuó incrementando este valor de manera progresiva
con esta misma magnitud, o sea:
Iteración 1:
187
Universidad Técnica Federico Santa María
Departamento de Industrias
Iteración 2:
… … … …
Iteración “i”:
Se iteró hasta i=20, de tal manera de alcanzar el 200% del valor del RETScreen y
asegurarse de encerrar al óptimo (Entre y ). Luego se escogió el valor de VPN
máximo ( ) y se tomó su tamaño de parque solar respectivo que correspondería
al tamaño óptimo de parque solar ( )
Con esta información RETScreen entrega el valor del área bruta total ( ) y la
capacidad nominal de la planta ( ) considerando un factor de conversión igual a
0,7 kWth/m2.
Almacenamiento: Sí
Capacidad de almacenamiento/Área Col: 119
o Este valor se escogió usando como referencia el estanque de reservorio
del parque solar de la minera Gaby, el cual maneja un factor de 119 L/m2
(4.300 m3 divididos en un área de apertura de 36.025 m2). La
recomendación de RETScreen que dice que normalmente las
instalaciones de calentamiento solar de agua consideran un rango entre 25
a 100 L/m2. Por lo tanto se considerará un valor de 119 L/m2 para ser
consistente con los niveles de redundancia exigidos en la minería.
Intercambiador de calor: Sí.
o Este ítem RETScreen lo consulta para considerar parámetros de eficiencia
en el sistema de intercambio de calor.
188
Universidad Técnica Federico Santa María
Departamento de Industrias
Eficiencia del intercambiador de calor:
o La eficiencia para cada intercambiador de calor se determinó para cada
caso como se muestra en la sección de metodología, y para la eficiencia
global se consideró como la multiplicación de la eficiencia de los
intercambiadores de calor involucrados (si en un proyecto no se considera
intercambiador, este parámetro se deja igual a 1):
Pérdidas varias: 5%
o RETScreen recomienda que para sistemas con almacenamiento de calor,
un valor entre 5 y 10% de pérdidas debe ser escogido. Considerando que
los diseños contemplan una aislación óptima, para el dimensionamiento
de los proyectos se consideró el mínimo como razonable.
Con
̇
∑
Donde corresponde a la potencia total a consumir por las bombas del sistema
de calentamiento solar evaluado y se determina como la sumatoria de la potencia de
cada circuito. Y corresponde al área de apertura total de los colectores
solares en el parque solar.
189
Universidad Técnica Federico Santa María
Departamento de Industrias
Tarifa de Electricidad ( ): 0,096 USD/kWh
Tabla 35 - Pronósticos del Precio de la Electricidad en Atacama desde el 2014 hasta el 2024
190
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.3.5. Resumen
Las eficiencias estacionales del caso base a considerar, para cada combustible son:
191
Universidad Técnica Federico Santa María
Departamento de Industrias
Diesel: 86,3%
Gas Natural: 83,7%
Precio de Combustible:
Se extrajo de la CNE los pronósticos de los precios del gas natural licuado ( ) y del
diesel ( ), desde el año 2014 hasta 2024 se muestran en la Tabla 36.
Tabla 36 - Proyecciones de los Precios del Gas Natural y Diesel desde el 2014 hasta el 2024 [24]
192
Universidad Técnica Federico Santa María
Departamento de Industrias
Según estos pronósticos, el Diesel experimentaría un incremento promedio en su precio
de 1,2% anual. Mientras que el Gas Natural tendería a disminuir 1,3% anual.
Considerando una disminución en el precio de un combustible fósil hoy en día puede
resultar un tanto polémico, dado que debido a la polémica mundial del agotamiento de
los combustibles fósiles es de esperar que este aumente. Por este motivo, considerando
que los pronósticos son confiables, adicionalmente se evaluarán los proyectos
considerando un escenario que el precio del gas natural se mantiene constante.
Al precio del Gas Natural se le incluyó un costo por regasificación del gas licuado, que
según la CNE equivale a 0,41 USD/MWh.
Los proyectos se evaluaron utilizando dos casos bases para cada alternativa; los cuales
consistieron en que las calderas operaran combustionando solo diesel en un caso y solo
gas natural en otro. En todos los casos propuestos se consideró que las calderas utilizan
gas natural para producir energía (por su menor costo).
El costo total de combustible para los dos casos bases se muestran a continuación en la
Tabla 37.
193
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.4. Análisis de Costos
Según el IREC el desglose de los costos de la inversión solar para un campo solar de
10.000 m2 cada uno de estos colectores solares presenta la siguiente distribución relativa
[19] (Gráfico 21):
18,0% Colectores
27,0%
Estructuras
Accesorios
8,0%
Almacenaje
Tuberías/IC/Bombas
2,0%
8,0% Sistema de Control
Obras Civiles
5,0%
2,0% Ingeniería
Otros
9,0%
11,0% Comercialización
3,0% 7,0%
Algunos ítems de los proyectos se basaron en esta estructura relativa de costos, como
son los costos asociados a las estructuras, accesorios, almacenaje y sistemas de control.
194
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 38 - Indice de Marshall & Swift actualizados al año 2013
El costo de la mano de obra propia de la minería en Chile el año 2012 fue de 3.790 M
USD [25] y para el mismo año la mano de obra propia de las mineras contemplo
145.470.568 h-persona [26]. Por lo tanto el costo de la mano de obra promedio por día
( ), considerando jornadas de 8,5 h el día, sería:
Para fines de limpieza de los colectores solares se consideró que un valor prudente para
ese tipo de trabajos sería un 60% del costo de mano de obra promedio de la minería, o
sea:
Para obtener un costo de mano de obra para los trabajos de ingeniería de diseño y
gestión cercano al real, se determinó en base a la remuneración mensual esperada de un
ingeniero civil en minas de 10 años de experiencia, que según la base de datos de la
página web del gobierno de chile www.mifuturo.cl este valor equivale a 4.860.839
CLP/mes, que considerando un tipo de cambio de 550 CLP/USD, 1.981 h de trabajo en
un año, tendremos:
195
Universidad Técnica Federico Santa María
Departamento de Industrias
Según RETScreen el precio de la mano de obra para evaluaciones de proyectos oscila
entre los 300 y 1.000 USD/día persona, por lo que el valor obtenido cae dentro de lo
esperado.
196
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.4.2. Estudio de Factibilidad
Entonces:
197
Universidad Técnica Federico Santa María
Departamento de Industrias
C. Diseño preliminar
El estimado de costos detallado se basa en los resultados del diseño preliminar y otras
investigaciones realizadas durante el estudio de factibilidad. Completar este tipo de
trabajos puede requerir entre 3 a 100 días-persona. Sin embargo, el máximo considera
proyectos de mayor complejidad, como son por ejemplo, las centrales hidroeléctricas.
Según RETScreen un proyecto chico de colectores solares puede demandar simplemente
un par de horas para completar este estudio. Por dicho motivo, y considerando que se
evaluará un proyecto de gran escala, es que se aceptó que 30 días de trabajo para 2
ingenieros deberían ser suficientes para completar esta etapa. Por lo tanto:
La línea base de los gases de efecto invernadero (GEI) debe contener toda la
información detallada sobre los consumos de energía; como por ejemplo, planes de
expansión de la red energética, modelos de despacho, consumo de combustible
198
Universidad Técnica Federico Santa María
Departamento de Industrias
marginal, patrones de consumos de energía y sus factores de emisión. Este corresponde a
la evaluación del proyecto pero sin alterar el estado del arte (o sea sin implementar el
proyecto de energía limpia), de tal manera de contar con las emisiones de gases de
efecto invernadero que se producirían. El plan de monitoreo (MP) debe asegurarse de
que existe toda la información necesaria para asegurar la reducción de emisiones y debe
describir una metodología para cuantificar estas emisiones, de tal manera poder
comparar el proyecto con la línea base y demostrar la reducción de GEI. Según
RETScreen los costos para desarrollar estos estudios en proyectos de gran escala oscilan
entre 30.000 y 55.000 USD. Por lo tanto, asumiendo el escenario de mayor costo:
F. Preparación de informes
Estos corresponden a los informes del estudio de factibilidad del proyecto, debiendo
incluir estudio detallado de mercado, técnico, legal, organizacional y financiero;
describiendo claramente todo el proyecto y su desempeño y riesgo, con el fin de que los
inversionistas puedan evaluar y tomar una decisión. Según RETScreen preparar un
estudio de factibilidad puede requerir de 2 a 15 días-persona. Por el tamaño de este
proyecto, se considerarán 15 días-persona. Por ende:
Este costo corresponde al de administrar todas las etapas del proyecto de factibilidad y
debe considerar reuniones con los stakeholders para conseguir su apoyo y colaboración.
Este costo debería involucrar entre 2 a 8 días-persona y considerar no más de 3 días-
persona para presentar el proyecto a los stakeholders. Sin embargo, se deberá el análisis
de costos detallado es el estudio limitante que más tarda. Por lo tanto se considerará el
tiempo para administrar todas las etapas como 25 días y 3 días para presentar el proyecto
a los stakeholders, y considerando que la responsabilidad es tomada por un ingeniero
experimentado, el costo será:
199
Universidad Técnica Federico Santa María
Departamento de Industrias
H. Viajes y alojamiento
Este costo debe incluir todos los costos relacionados con los viajes y alojamiento de los
integrantes del equipo del estudio de factibilidad que así lo requieran. Por este motivo se
considerará que se deben realizar viajes a terreno, comprendiendo vuelos desde Santiago
a Calama, el arriendo de vehículos para transportar al personal a terreno y el costo de
alojamiento y alimentación.
Para conocer los costos de viajes y alojamiento se realizaron las siguientes cotizaciones:
A continuación se presenta una tabla con el resumen de costos unitarios cotizados (Tabla
39).
Costo Unitario
Vuelo Santiago - Calama (LATAM) (USD) 273
SUV Toyota Runner 4x4 (USD/día) 186
Domos Apart Hotel (USD/noche*persona) 84
Alimentación (USD/día) 20
200
Universidad Técnica Federico Santa María
Departamento de Industrias
el estudio de impacto ambiental requeriría dos personas y cuatro días para realizar el
trabajo, se considerarán 5 días; uno para llegar y cuatro para realizar el trabajo y volver.
Los costos totales por sub ítem se muestran a continuación (Tabla 40 y Tabla 41
respectivamente).
Tabla 40- Costo total por Viajes y Alojamiento por la Inspección de Sitio
Tabla 41 - Costo total por Viajes y Alojamiento por el Estudio de Impacto Ambiental
Por lo tanto el costo total por Viajes y Alojamiento del Estudio de Factibilidad
correspondería a:
9.6.4.3. Desarrollo
Este ítem contempla los procedimientos que se deben llevar a cabo una vez que el
estudio de factibilidad está listo, de tal manera de comenzar a poner en marcha el
proyecto.
201
Universidad Técnica Federico Santa María
Departamento de Industrias
A. Negociaciones del contrato
B. Permisos y autorizaciones
Para que el proyecto pueda comenzar a construirse, se requiere una serie de permisos y
autorizaciones Principalmente estos son:
Autorizaciones ambientales
Autorización para el uso del terreno
Permisos de edificación
Permiso para el uso del recurso hídrico
Otros acuerdos de operaciones
Autorización de inspectores eléctricos, de calderas, de seguridad y de la agencia
reguladora ambiental.
202
Universidad Técnica Federico Santa María
Departamento de Industrias
C. Topografía del sitio y derechos de servidumbre
Según RETScreen realizar una topografía del sitio dependerá principalmente del estado
de la propiedad del sitio, planes que se tengan para uso del lugar, de la ubicación y su
tamaño y posibles problemas legales y de seguros. Se requieren los derechos del terreno
donde el proyecto será implementado, incluyendo caminos, líneas eléctricas y de
colección, la red de agua a calentar, subestaciones y construcciones para la operación y
mantenimiento. Según RETScreen el costo de realizar una topografía de 1 a 10 hectáreas
puede costar alrededor de 750 USD. Para este estudio consideraremos un costo de 750
USD/He (0,750 USD/m2). Entonces el costo de la evaluación topográfica a considerar:
Este sub ítem hace referencia al costo de conseguir el financiamiento el cual dependerá
del esfuerzo de los expertos por hacer todos los arreglos, identificar a los inversionistas y
solicitar los fondos. Según RETScreen este puede tardar de 3 a 100 días-persona. Por lo
tanto, considerando que se emplea una mano de obra calidad negociante ( )
tendremos:
203
Universidad Técnica Federico Santa María
Departamento de Industrias
F. Legal y contabilidad
Durante el desarrollo de las distintas etapas del proyecto será necesario llevar una
contabilidad. La necesidad de contar con soporte legal dependerá de los arreglos
financieros, de los dueños, contrato de seguros, asunción de responsabilidades y la
complejidad de los contratos y acuerdos. Según RETScreen esto involucraría entre 3 a
200 días-persona dependiendo de la complejidad y dimensionamiento del proyecto. Se
consideró que 150 días-persona debería ser suficiente para la complejidad de este
proyecto, por lo tanto:
Este ítem corresponde al costo de administrar todos los sub ítems de la etapa de
desarrollo, la cual según RETScreen normalmente requeriría unos 100 días-persona para
proyectos de gran escala. Sin embargo, el proceso cuello de botella vendría siendo de
permisos y negocios, el cual requeriría 60 días. Por lo que considerando 2 personas para
abarcar las funciones de gerencia tendremos:
H. Viajes y alojamiento
Se deberán realizar numerosas visitas a terreno durante la etapa de desarrollo, las cuales
deberán ser principalmente destinadas a reuniones. Dentro de la etapa de desarrollo será
completamente necesario viajar a terreno durante la topografía del sitio y altamente
204
Universidad Técnica Federico Santa María
Departamento de Industrias
probable que sea necesario durante la validación y registro de GEI. Para la topografía del
sitio se considerará un equipo de 20 personas para un trabajo de 5 días, y para la
validación y registro de GEI se considerarán 5 días de visita a terreno para 3 personas.
También se contemplarán 10 días hábiles y un equipo de 3 personas para viajes de la
gerencia del proyecto, por potenciales reuniones y consultas en terreno. En todos los
casos se contempla que se llega a Calama un día antes del trabajo a realizar.
Por lo tanto, empleando las cotizaciones presentadas en la Tabla 39, los costos para
viajes y alojamiento de estos sub ítems se muestran a continuación en la Tabla 42, Tabla
43 y Tabla 44.
Tabla 42 - Costo total por Viajes y Alojamiento por la Topografía del Sitio
Tabla 43 - Costo total por Viajes y Alojamiento por la Validación y Registro de GEI
205
Universidad Técnica Federico Santa María
Departamento de Industrias
Alimentación (USD/día) 11 660
Total 7.670
9.6.4.4. Ingeniería
Este ítem de RETScreen coincide con el ítem “Ingeniería” del IREC. Sin embargo
RETScreen proporciona una estructura más detallada. Dado que el nivel de complejidad
de los estudios de esta sección está directamente relacionado con el dimensionamiento
del parque solar, es lógico esperar que el costo de estos sean también proporcionales.
Por ello se respetará el valor relativo para este ítem presentado en la estructura de costos
del IREC y se utilizará la distribución relativa de RETScreen para comparar y conocer
en detalle la distribución de costos de los sub ítems de Ingeniería.
Este costo depende de los recursos necesarios para completar todos los dibujos técnicos
y especificaciones, incluyendo la configuración final del sistema de calentamiento solar,
el dimensionamiento de todos los equipos del sistema, el lugar exacto de la colocación
de los colectores solares, los estanques de almacenaje, intercambiadores de calor, unidad
de control automático, tuberías y las unidades de manejo de líquidos. También debe
incluir el cómo este sistema se integrará con la fontanería existente. Se deberá trabajar
junto al diseñador del sistema de calentamiento de agua, de las edificaciones y de la
fontanería existente para alcanzar un diseño óptimo.
206
Universidad Técnica Federico Santa María
Departamento de Industrias
B. Diseño mecánico
Este corresponde al diseño de la estructura que soportará a los colectores solares. Según
RETScreen, el diseño para estructuras inclinadas en el suelo requiere de unas 22 a 44
horas de ingeniería y unas 30 a 60 horas en el dibujo de planos. Considerando que estos
los realiza un equipo de 2 ingenieros ( ), demandando la cota máxima de
tiempos, los costos estimados serán:
( )
C. Diseño Eléctrico
RETScreen no recomienda valores en este ítem, ya que se asume que un diseño eléctrico
es bastante básico en este tipo de instalaciones. Dado que debe existir un diseño para
alimentación de bombas, iluminación y sistema control automático de las líneas de agua,
se consideró que este trabajo podría ser realizado en una semana por un ingeniero
diseñador ( ), por lo tanto:
D. Diseño Civil
207
Universidad Técnica Federico Santa María
Departamento de Industrias
E. Licitaciones y contratos
F. Supervisión de la construcción
Este costo busca asegurar que el proyecto se está construyendo como se diseñó. La
supervisión de la construcción para proyectos de larga escala implica presencia de
tiempo completo en la obra y la deberá realizar un administrador del proyecto. Se estima
que para proyectos de larga escala 500 días-persona debería ser suficiente para
completar la obra. Asumiendo que las obras tarden 167 días, se necesitarían
aproximadamente 3 supervisores a costo de ingeniero de terreno ( ).
Como los supervisores deben estar presente de principio a fin durante la construcción de
toda la obra, se consideró el arriendo de un departamento en Calama a un costo de 840
USD/mes durante 17 meses y 10 USD/día para el transporte a la faena. Los costos de
viaje y alojamiento para el supervisor se muestran en la Tabla 45:
208
Universidad Técnica Federico Santa María
Departamento de Industrias
*Transporte (10 USD/día) 501 15.030
Arriendo Departamento (840 USD/mes) 17 42.840
Alimentación (20 USD/día) 501 30.060
Total 88.750
Los costos presentados en este ítem que están basados en la estructura del IREC, se
determinan de manera proporcional al costo del material del colector solar ( ) según su
peso relativo.
209
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 46 - Precio unitario del área de apertura de los Colectores Solares Evaluados
Como estos precios corresponden a los ofrecidos por fabricantes de Europa, se estimó un
factor por importación de colectores solares, basándose en la estimación indirecta del
precio de los colectores solares del fabricante Sunmark, basándose en los costos de
inversión del parque solar de la minera Gaby y del parque solar de la cervecera austríaca
Goess Brewery, ya que cuentan con colectores solares con las mismas áreas de apertura
y del mismo fabricante.
210
Universidad Técnica Federico Santa María
Departamento de Industrias
inversión total, entonces el valor correspondiente al colector solar importado en Chile
(precio estimado de Gaby):
211
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo tanto, el costo de la estructura solar será:
( )
C. Cañerías (Tubo)
Estos corresponden a los costos asociados a las cañerías del circuito de carga, descarga y
de los colectores solares. Para todas las cañerías se consideraron líneas de 8 pulgadas de
diámetro de acero inoxidable 316L Sch.40.
Del complemento online del libro “Plant Design and Economics for Chemical Engineers
– 5th Edition” se extrajo la function que describe el comportamiento del costo para este
tipo de cañerías, dependiendo de su largo. Esta función se muestra a continuación
(Gráfico 22).
250.000
R² = 1
150.000
100.000
50.000
0
0 2.000 4.000 6.000 8.000 10.000 12.000
Largo de Matriz (m)
212
Universidad Técnica Federico Santa María
Departamento de Industrias
fluido frío y la otra mitad a la cañería que transporta el fluido caliente. Por lo tanto, el
costo será:
( )
Donde corresponde al índice de Marshall & Swift para actualizar el año 2012
al 2013.
Con
( )
(( ) )
(( ) )
213
Universidad Técnica Federico Santa María
Departamento de Industrias
1.200.000
600.000
400.000
200.000
0
0 2.000 4.000 6.000 8.000 10.000 12.000
Largo de Matriz (m)
Por lo tanto, el costo total por aislar térmicamente las cañerías será:
( )
D. Intercambiadores de Calor
El costo de los intercambiadores de calor depende del área de transferencia de calor con
el cual se dimensionaron y del tipo de intercambiador de calor.
214
Universidad Técnica Federico Santa María
Departamento de Industrias
300.000
250.000
y = 19313x0,3616
Costo IC (USD/IC)
200.000 R² = 1
150.000
100.000
50.000
0
0 500 1.000 1.500 2.000
Área (m2)
215
Universidad Técnica Federico Santa María
Departamento de Industrias
E. Bombas
Las bombas a considerar son del tipo centrífugas, y la función de costos a utilizar
considera implícitamente el costo del motor electro. El costo asociado se encuentra en
función de la potencia de la bomba en (kW) y se puede ver a continuación en el Gráfico
25 (en azul). Cada bomba que se dimensiona en esta ecuación, considera una diferencia
de presión de hasta 1035 kPa.
90.000 4,00E-03
70.000 3,00E-03
60.000
2,50E-03
50.000
2,00E-03
40.000
1,50E-03
30.000 y = 9042,8e0,0123x
R² = 0,9978 1,00E-03
20.000
10.000 5,00E-04
0 0,00E+00
0 50 100 150 200
Potencia (kW)
Como se puede observar, se graficó otra curva (en rojo) que describe la potencia que
puede entregar una bomba por dólar gastado. De esta manera se consideró que cualquier
circuito que exceda por al menos el doble de la capacidad máxima de esta última
función, empleará otra bomba de la capacidad que maximiza la potencia para el mismo
costo.
216
Universidad Técnica Federico Santa María
Departamento de Industrias
Si , la potencia a emplear por bomba en ese circuito será:
Y el costo total de bombas será la sumatoria del costo de bombas de cada circuito, o sea:
Para estos costos se consideró la proporción de la estructura de costos del IREC, para lo
cual considera un 2% para accesorios y un 3% para el sistema de control automático. Por
lo tanto:
( )
G. Tanque de reservorio
( )
217
Universidad Técnica Federico Santa María
Departamento de Industrias
H. Instalación
Los costos de instalación están divididos en los costos de instalación de los colectores
solares, la instalación del loop solar y la instalación de equipamiento auxiliar. Según el
IREC estos costos corresponden al 18% de la inversión total. Por lo tanto, el costo de
instalación total a considerar en los proyectos es:
Este costo corresponde a las obras civiles que se deben hacer durante la obra, como
construcción de caminos, edificaciones, excavaciones, etc. Este costo idealmente debe
ser cotizado con los contratistas locales, sin embargo el estudio del IREC incluye dentro
de las cotizaciones de cada colector solar un ítem correspondiente a obras civiles. Por
ende, se utilizará este valor que es relativo al costo de la inversión total de colectores
solares el cual es equivalente a un 9%. Por lo tanto, el costo de la construcción de las
obras civiles a considerar equivale a:
( )
J. Repuestos
RETScreen recomienda incluir dentro del costo de inversión total, un costo asociado a
repuestos que en la mayoría de los casos son requeridos luego de haber realizado las
compras pertinentes. Estos corresponderían a repuestos para equipos del sistema y para
proyectos de gran escala, considerar un 3% de la inversión del total de equipos debería
ser suficiente. Por lo tanto, el costo asociado a los repuestos a considerar equivale a:
218
Universidad Técnica Federico Santa María
Departamento de Industrias
K. Transporte
Este se refiere al costo de transportar todos los materiales y equipos a instalar. Este costo
está considerado por el IREC dentro del costo de instalación. Por lo tanto este sub ítem
se deja en blanco dentro del software RETScreen.
Según RETScreen una sesión de 4 a 8 horas debería ser suficiente para que un experto
aprenda a operar de manera adecuada una planta de calentamiento solar de agua. Por lo
tanto, considerando que un ingeniero capacitado le enseñe a un equipo de 10 personas, el
costo asociado al entrenamiento será:
219
Universidad Técnica Federico Santa María
Departamento de Industrias
o Colector solar (Material del colector solar)
o Estructuras de soporte del colector solar
Para incluir los efectos del segundo descuento por economía de escala en los proyectos,
se construyó un modelo con la información proporcionada por el IREC. Con ella se
construyó la curva de economías de escala, mostrada a continuación (Gráfico 26).
100,00%
99,00%
98,00%
Economía de escala (%)
97,00%
96,00%
95,00%
94,00%
93,00%
92,00%
91,00%
90,00%
500 2.500 4.500 6.500 8.500 10.500
Área de apertura (m2)
Cabe considerar, que un colector solar cilíndrico parabólico (CCP) Polytrough 1800 del
proveedor NEP Solar, no se adapta a este modelo. Por lo que se determinó otra modelo
de economías de escala, mostrada en el gráfico a continuación (Gráfico 27).
220
Universidad Técnica Federico Santa María
Departamento de Industrias
100,00%
95,00%
85,00% y = 1,7193x-0,084
R² = 0,9518
80,00%
75,00%
500 2.500 4.500 6.500 8.500 10.500
Área de apertura (m2)
Gráfico 27 - Economías de Escala en la Inversión Total del colector solar Polytrough 1800 (CCP)
El modelo que describe estas economías de escala para el colector solar Polytrough 1800
se muestra a continuación, con un R2=0,9518.
N. Contingencias
221
Universidad Técnica Federico Santa María
Departamento de Industrias
En el presente trabajo se realizaron las evaluaciones considerando contingencias de un
20%.
Este costo dependerá del tiempo que tarde la construcción de la obra, la tasa de interés.
RETScreen calcula este valor asumiendo que la deuda promedio durante la duración del
proyecto (en meses) equivale a un 50% del subtotal de todos los costos del proyecto.
Para ello se considera una tasa de interés de la deuda de un 15%.
Este ítem debe considerarse en el caso que no se cuente con la propiedad el terreno a
utilizar para construir el parque de colectores solares. En este estudio se considera que
ya se cuentan con los derechos del terreno, por lo que los costos asociados a arriendo no
dependen de ese proyecto.
Para este estudio, se considerará que el costo anual del seguro ( ) equivale a un
1% del costo de materiales de la inversión total, o sea:
222
Universidad Técnica Federico Santa María
Departamento de Industrias
C. Partes y labor
Este costo hace referencia a los gastos anuales en partes de repuestos y el trabajo
requerido en rutinas de emergencia y de operación y mantenimiento del proyecto. Los
colectores solares deben limpiarse de la suciedad (polvo), realizar inspecciones,
monitoreo y mantenciones programadas.
( )
Dónde
Instalar un colector solar de placa plana típico de 2,1 m2 de área de apertura requiere
aproximadamente 3,45 L de fluido para su operación (1,64 L/m2). Entonces el volumen
total contenido en el loop solar ( ) se considerará como la suma del volumen por
área de apertura requerido más el volumen de la cañería principal ( ), o sea:
223
Universidad Técnica Federico Santa María
Departamento de Industrias
La limpieza de los colectores solares en la mina Gaby se realiza a una tasa aproximada
de 0,36 min/m2 .Si se considera que se limpian los colectores solares una vez cada dos
semanas, el tiempo total requerido para limpieza anual sería:
Por lo que el personal de limpieza necesaria para llevar a cabo este trabajo, considerando
8,5 horas laborales por día por persona, sería (en día*persona/año):
Por lo tanto el costo anual de limpieza de los colectores solares, considerando un costo
de mano de obra de limpieza ( ) de 133 USD/día-persona, equivale a:
O sea:
224
Universidad Técnica Federico Santa María
Departamento de Industrias
E. General y Administrativo
F. Contingencias
G. Costo de Combustible
225
Universidad Técnica Federico Santa María
Departamento de Industrias
También debe incluirse el costo asociado al consumo de energía eléctrica del sistema de
bombas del caso propuesto ( ), el cual será:
226
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.4.7. Ahorros Anuales (Ingresos)
RETScreen realiza una comparación entre el caso base y el caso propuesto para la
determinación de ingresos por ahorros. Para ello considera que una vez puesto en
marcha el caso propuesto, el caso base “se deja de hacer”, por lo que se considerarán
como ingresos todos los gastos asociados al consumo anual de combustible del caso
base.
Donde
227
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.5. Análisis de Emisiones
Sistema eléctrico
228
Universidad Técnica Federico Santa María
Departamento de Industrias
Emisiones del caso base
Análogamente al caso anterior, RETScreen determina las emisiones GEI para el caso
base ( ), que corresponderán a las emisiones que se dejarían de producir al
implementar el caso propuesto. Considerando el consumo de energía ( ) y el
respectivo factor de emisión del combustible a evaluar como base para contrastar el caso
propuesto ( ), las emisiones del caso base serán:
229
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.6. Análisis Financiero
Esta corresponde a la tasa promedio con la cual se espera que incremente el costo del
combustible de manera anual. Se evaluaron los proyectos considerando la tasa de
escalamiento anual pronosticada por la CNE (comisión nacional de energía), que
corresponderían en el caso respectivo de su consumo, para el diesel ( ) y el gas
natural ( ):
B. Tasa de inflación
Según el INE (Instituto Nacional de Estadísticas) la variación anual del IPC en Chile
desde diciembre del año 2009 hasta el año 2013 corresponde a:
Entre los meses de enero y septiembre del 2014 el IPC varió en un 3,85%, por lo que
considerando que el resto del año seguirá aumentando a una tasa de 0,48% mensual, la
variación total del año 2014 sería igual a 5,47%.
230
Universidad Técnica Federico Santa María
Departamento de Industrias
C. Tasa de Descuento (WACC)
( [ ] )
O sea
( )
231
Universidad Técnica Federico Santa María
Departamento de Industrias
O sea:
Para fines de realizar una evaluación con un riesgo reducido, se optó por considerar la
vida útil del proyecto igual a 5 años. Sin embargo los proyectos de estas dimensiones
pueden tener una vida útil real de 20 años.
E. Financiamiento
Para todos los proyectos se consideró una relación del 100% de deuda, con una tasa de
interés igual al costo de la deuda del mercado minero ( ). La duración de la
deuda se fijó para cada proyecto igual a su vida útil (5 años).
En este ítem se considera el impuesto a la renta que se deberá rendir por los ingresos
generados en los casos propuestos, el cual se consideró igual a un 27%.
232
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.6.2. Renta Anual por reducción de GEI
La evolución del precio de los derechos de emisión europeos (EUA) no ha sido buena
los últimos años; el año 2010 el valor de un EUA rondaba los 14,62 EUR/tCO2eq y el
año 2013 alcanzó su valor más bajo llegando a los 4,58 EUR/tCO2eq. Sin embargo, entre
el año 2013 y 2014 se produjo un alza, alcanzando un valor de 6,55 EUR/tCO2eq. Lo
mencionado previamente se puede observar en el Gráfico 28.
16,00
14,00
12,00
EUA's (EUR/tCO2eq)
10,00
8,00
6,00
4,00
2,00
0,00
2008 2009 2010 2011 2012 2013 2014 2015
Año
Gráfico 28 - Evolución del precio de los EUA entre los años 2009 y 2014 [29]
233
Universidad Técnica Federico Santa María
Departamento de Industrias
9.6.6.3. Viabilidad Financiera
∑
∑
Donde
Si bien un requisito para que el proyecto sea factible es que el costo nivelado de la
energía sea igual o menor que el costo del combustible del caso base, se considerará
viablemente financiero un proyecto cuyo VAN sea mayor o igual a cero.
234
Universidad Técnica Federico Santa María
Departamento de Industrias
10. Resultados
10.1. Resultados de las evaluaciones de proyecto
Para cada proyecto y alternativa se construyó una tabla con los resultados de la
optimización y sus respectivos rendimientos financieros, presentado para cada
tecnología el modelo de colector solar que logró obtener el menor costo nivelado de
energía (LEC). Las variables y parámetros presentados en las tablas de resultados se
muestran a continuación.
: Tecnología del colector solar evaluado (CPP, CTV, CCP, MAC o CCF).
235
Universidad Técnica Federico Santa María
Departamento de Industrias
10.1.1. Calentamiento del Agua del Circuito de Calderas
Los resultados de este proyecto considerando el uso de gas natural o diesel en el caso base, se muestran en la Tabla 49.
236
Universidad Técnica Federico Santa María
Departamento de Industrias
10.1.2. Calentamiento de Electrolito Rico
Los resultados de este proyecto considerando el uso de gas natural o diesel en el caso base, se muestran en la Tabla 50.
237
Universidad Técnica Federico Santa María
Departamento de Industrias
10.1.3. Calentamiento de Agua de Osmosis
Configuración normal (agua de osmosis se calienta a través de un intercambiador de “agua de osmosis y circuito de descarga”.
Los resultados de este proyecto considerando el uso de gas natural o diesel en el caso base, se muestran en la Tabla 51.
238
Universidad Técnica Federico Santa María
Departamento de Industrias
Configuración modificada (agua de osmosis se hace pasar directamente como fluido de carga y descarga).
Los resultados de este proyecto considerando el uso de gas natural o diesel en el caso base, se muestran en la Tabla 52.
239
Universidad Técnica Federico Santa María
Departamento de Industrias
Se lograron resultados con un VAN positivo utilizando colectores solares de placa plana
(CPP) y te tubo de vacío (CTV), para pre calentar el agua de calderas de la faena. Y
también se logró un resultado de VAN positivo solo al usar CPP para precalentar el
electrolito rico. Sin embargo, todos los escenarios donde se obtuvo un valor de VAN
positivo, se realizaron utilizando combustible Diesel en el caso base.
El modelo de colector solar que permitió obtener el mayor VAN fue el colector de placa
plana XL USB, del fabricante EZINC Metal San Tic. A.S., utilizado para precalentar el
agua de calderas. Obteniendo un VAN igual a 1,57 M USD, un IVAN del 31%, un LEC
igual a 51 USD/MWh y considerando un retorno de capital simple al 2 año de vida útil
del proyecto.
El modelo con otra tecnología que lo sigue es el colector de tubo de vacío APSE-30, con
un VAN igual a 1,02 M USD, un IVAN de 20%, un LEC igual a 54 USD/MWh y para
un retorno de capital simple a los 3 años.
Si estas operaran funcionando solo con gas natural, el proyecto no se hace factible para
ningún caso, con un VAN igual a -3,58 M USD para el modelo XL USB en el mejor
caso (precalentamiento de agua de calderas). Sin embargo, hay que considerar que los
proyectos con gas natural fueron evaluados con una tasa de escalamiento de combustible
igual a -1,3%, escenario que claramente perjudica los ingresos del proyecto.
En todos los proyectos evaluados los colectores solares que rindieron según su
tecnología fueron siempre los mismos y en el mismo orden. El colector solar que mejor
rindió en todos los proyectos fue el colector solar de placa plana (CPP) modelo XL USB
del fabricante EZINC Metal San Tic. A.S.
De los resultados se puede ver que el costo nivelado de la energía disminuye a medida
que el tamaño óptimo de la planta es mayor (Gráfico 29).
240
Universidad Técnica Federico Santa María
Departamento de Industrias
100
90
LEC (USD/MWh)
80
CPP
70 CTV
CCP
60
MAC
50 CCF
40
0 5.000 10.000 15.000 20.000 25.000
Área Total del Parque Solar (m2)
Esto se puede explicar por los elevados costos fijos que existen en la minería, que a
medida que el tamaño de la planta aumenta su peso frente al costo variable disminuye.
Esto se puede evidenciar implícitamente, ya que el costo de inversión por área de parque
solar alcanza su valor mínimo cuando el tamaño de la planta solar es aproximadamente
20.000 m2 (Gráfico 30).
1.400
1.200
IT/AT.BRUTA (USD/m2)
1.000
CPP
800
CTV
600
CCP
400 MAC
200 CCF
0
0 5.000 10.000 15.000 20.000 25.000
Área Total del Parque Solar (m2)
Gráfico 30 - Costo de inversión por área de parque solar versus área total del parque solar
241
Universidad Técnica Federico Santa María
Departamento de Industrias
Costo Nivelado de la Energía
90
80
Costo Energía (USD/MWh)
70
60 Diesel
50
Gas Natural
40
30
20
10
0
CPP CTV CCP MAC CCF
Gráfico 31 - Costo Nivelado de la Energía para todos los colectores - Calentamiento de Agua de Calderas
Como se puede observar solo colectores solares de placa plana (CPP) y de tubo de vacío
(CTV) obtuvieron costos de energía (LEC) que compiten con el costo de energía del
diesel, y todos se encuentran por sobre el costo de energía del gas natural. Esto es
concordante y tiene sentido con los resultados, ya que ningún proyecto evaluado resultó
rentable para el gas natural como combustible del caso base.
242
Universidad Técnica Federico Santa María
Departamento de Industrias
10.3. Análisis de Riesgo y Sensibilidad
243
Universidad Técnica Federico Santa María
Departamento de Industrias
Se realizaron 4.115 pruebas cambiando los valores de estas variables mediante
distribuciones beta en el Crystal Ball. Los valores de esperanza, mínimo y máximo de
las distribuciones de cada variable se muestran a continuación (Tabla 53).
244
Universidad Técnica Federico Santa María
Departamento de Industrias
El gráfico de distribución probabilística del VAN resultante de la simulación, se muestra
a continuación (Gráfico 32).
Gráfico 32 - Distribución probabilística del VAN - Proyecto calentamiento de agua de calderas mediante
colector solar plano XL USB
Como se puede observar, bajo las consideraciones realizadas para la simulación, existiría
una probabilidad del 83,09% de que el proyecto obtenga un VAN mayor o igual a cero.
Por lo que el riesgo de realizar un proyecto de este tipo se considerará como:
245
Universidad Técnica Federico Santa María
Departamento de Industrias
El gráfico de contribución a la varianza del VAN por cada variable participante en la
simulación se muestra a continuación (Gráfico 33).
0,80
0,60
Contribución a la Varianza del VAN
0,40
0,20
0,00
P. Diesel
EUR/USD
T.E. Diesel
T.E. GN
T.E. EUA
CLP/USD
T.E. Elect.
P. EUA
Seguro
C. Operarios
P. Elect
C. Limpieza
C. Cambio Glicol
F.Importacion
C. Negociante
C. Ingeniero
T. Interés
Contingencias
P. GN
-0,20
-0,40
-0,60
-0,80
Las variables que más influyen sobre la rentabilidad del proyecto, corresponden a la
variación del precio del diesel con un 68% y el precio del gas natural con un -63%. Esto
significa que un aumento significativo en el precio del diesel impacta positivamente y de
manera importante en el valor del VAN, mientras que el gas natural afecta
negativamente al VAN.
246
Universidad Técnica Federico Santa María
Departamento de Industrias
10.4. Pruebas de Sensibilidad
Para conocer una contribución a la varianza proporcional entre cada variable, se optó por
realizar otra simulación Montecarlo, pero esta vez considerando cada variable con una
distribución normal y una desviación estándar del 10% de su esperanza.
Se realizó el análisis de sensibilidad al proyecto con fracción solar igual a 7,7% y luego
al proyecto con fracción solar igual a 100%.
247
Universidad Técnica Federico Santa María
Departamento de Industrias
100,0%
80,0%
40,0%
20,0%
0,0%
η0
P. Diesel
Años Proyecto
EUR/USD
Deuda (%)
T.E. Diesel
P. EUA
T.E. GN
Seguro
T.E. EUA
T.E. Elect.
CLP/USD
C. Operarios
P. Elect
C. Obrero
U2 (W/m2 K2)
U1 (W/m2 K)
Contingencias (%)
C. Limpieza
C. Negociante
C. Ingeniero
F.Importacion
PU
P. GN
-20,0%
-40,0%
-60,0%
-80,0%
0,60
Contribución a la Varianza del VAN
0,40
0,20
0,00
EUR/USD
P. Diesel
U1 (W/m2 K)
Deuda (%)
U2 (W/m2 K2)
P. EUA
T.E. EUA
CLP/USD
T.E. Elect
Años Proyecto
P. GN
C. Operarios
T.E. GN
T.E. Diesel
C. Negociante
C. Ingeniero
Seguro
C. Limpieza
C. Cambio Glicol
P. Elect
Contingencias (%)
T. Interés (%)
η0
F.Importacion
PU
-0,20
-0,40
-0,60
248
Universidad Técnica Federico Santa María
Departamento de Industrias
Análisis de sensibilidad – Caso Base: Gas Natural (fracción solar = 7,7 %, n =
4.485)
0,50
0,40
Contribución a la Varianza del VAN
0,30
0,20
0,10
0,00
EUR/USD
η0
P. GN
Deuda (%)
Años Proyecto
P. EUA
C. Operarios
T.E. Elect.
T.E. EUA
T.E. GN
P. Diesel
-28 / 200
T.E. Diesel
CLP/USD
P. Elect
Seguro
U2 (W/m2 K2)
C. Limpieza
C. Cambio Glicol
C. Negociante
C. Ingeniero
Contingencias (%)
U1 (W/m2 K)
F.Importacion
PU
-0,10
-0,20
-0,30
-0,40
-0,50
Gráfico 36 - Contribución a la varianza del VAN - Caso Base: Gas Natural – f : 7,7%
249
Universidad Técnica Federico Santa María
Departamento de Industrias
0,60
0,20
0,00
EUR/USD
P. GN
U1 (W/m2 K)
Deuda (%)
T.E. GN
U2 (W/m2 K2)
CLP/USD
P. EUA
T.E. Elect.
C. Ingeniero
C. Negociante
T.E. EUA
T.E. Diesel
P. Diesel
C. Cambio Glicol
C. Operarios
Seguro
C. Limpieza
Años Proyecto
P. Elect
Contingencias (%)
T. Interés (%)
η0
F.Importacion
PU
-0,20
-0,40
-0,60
Gráfico 37 - Contribución a la varianza del VAN - Caso Base: Gas Natural – f: 100%
Se realizó una simulación considerando varios escenarios para las variables más
influyentes de cada proyecto y se hicieron variar hasta lograr que el proyecto resultase
rentable, considerando el número de colectores óptimo y el número mínimo necesario
para que la fracción solar del sistema sea de un 50% y de un 100%.
Para el caso de diesel y gas natural, la fracción solar se hace igual a 100% cuando el
tamaño del parque solar es igual a 128.178 colectores solares planos XL USB.
Para este análisis se seleccionaron las variables más influyentes para el proyecto
evaluado con una fracción solar igual al 100%, y sus pesos relativos de contribución a la
varianza del VAN se muestran a continuación (Tabla 54).
Tabla 54 - Variables más influyentes en el VAN por caso base - Fracción solar 100%
250
Universidad Técnica Federico Santa María
Departamento de Industrias
Precio Diesel 22% Precio GN 17%
Precio Colector -52% Precio Colector -52%
F. Importación -33% F. Importación -34%
η0 -29% η0 -27%
Como se puede observar, las variables más influyentes para los dos casos bases son las
mismas, considerando que el precio del combustible diesel y gas natural son variables
análogas en cada caso base.
Se forzó cada una de estas variables para que el proyecto de rentable (VAN = 0), y
considerando todas las demás variables constantes, los únicos escenarios que lograron
rentabilizar el proyecto, para una fracción solar igual al 100% (n = 128.178), fueron los
siguientes (Tabla 55):
Como el gas natural como combustible del sistema auxiliar de calentamiento en el caso
base es el escenario más probable a enfrentar en los próximos años, se analizaron
distintas combinaciones de variación de estas variables para encontrar los escenarios que
rentabilicen el proyecto.
Se analizará el proyecto con una fracción solar igual al 50%., construyendo un escenario
con cada variable dentro de los rangos mostrados a continuación (Tabla 57).
251
Universidad Técnica Federico Santa María
Departamento de Industrias
Tabla 57 - Rango de simulación de variables del análisis de rentabilidad
Los resultados se presentan de manera gráfica a continuación, tomando para cada caso
dos variables y graficándolas versus el VAN del proyecto, considerando en el caso
respectivo todas las demás variables constantes y en su valor normal del proyecto para
un tamaño de parque solar igual a 38.159 colectores solares (fracción solar
aproximadamente igual a 50%).
El VAN del proyecto se torna igual o mayor a cero en la mayoría de los casos cuando el
precio del combustible es mayor a 116 USD/MWh, y en todos los casos cuando el tipo
de cambio es mayor a 1 EUR/USD (Gráfico 38).
20
VAN (M USD)
10
0
3
10-20
2,5456
-10
2,0912
1,6368
0-10
-20
1,1824
0,728
-30 -10-0
43,3
-20--10
61,44
79,58
97,72
-30--20
115,86
134
Gráfico 38 - VAN versus Precio del Gas Natural y Tipo de Cambio Euro/Dolar
252
Universidad Técnica Federico Santa María
Departamento de Industrias
El VAN del proyecto se torna igual o mayor a cero en todos los casos que el precio
del colector solar con las características del modelo XL USB es menor a 24,4
EUR/m2 y el precio del gas natural es mayor a 97 USD/MWh (Gráfico 39).
20
VAN (M USD)
10
0
2
10-20
13,2
-10
24,4
35,6
-20 0-10
46,8
58
-30 -10-0
43,3
-20--10
61,44
79,58
97,72
-30--20
115,86
134
Gráfico 39 - VAN versus Precio del Gas Natural y Precio del Colector Solar
El VAN del proyecto se torna igual o mayor a cero algunos casos que el precio del
colector solar con las características del modelo XL USB es menor a 10 EUR/m2
para todos todos los valores de eficiencia (Gráfico 39).
253
Universidad Técnica Federico Santa María
Departamento de Industrias
5
0
VAN (M USD)
-5 0-5
-10 -5-0
-15
-10--5
-20
-15--10
-25 2
24,4 -20--15
46,8 -25--20
Precio de
Colector Solar
(USD/m2)
η0 (-)
El VAN del proyecto se torna igual a cero solo en los casos que el precio del colector
solar con las características del modelo XL USB es menor que 2 EUR/m2 y los casos
en el que el tipo de cambio es muy cercano a 1,18 EUR/USD (Gráfico 40).
5
0 0-5
VAN (M USD)
-5 -5-0
-10 -10--5
-15 -15--10
Gráfico 40 - VAN versus Precio del Colector solar y Tipo de Cambio Euro/Dolar
-5-0
-10--5
0 -15--10
-5 -20--15
VAN (M USD)
-25--20
-10
-25 1
0,728 1,1824 1,6368 2,0912 2,5456 3
Tipo de Cambio (EUR/USD)
Gráfico 41 - VAN versus Tipo de Cambio Euro/Dolar y Factor de importación de Colectores Solares
255
Universidad Técnica Federico Santa María
Departamento de Industrias
10.4.3. Variación de vida del proyecto
0,00
0 5 10 15 20 25
-2,00
-4,00
CPP
VAN (M USD)
-6,00 CTV
CCP
-8,00
MAC
-10,00
CCF
-12,00
-14,00
Vida útil del proyecto (años)
Gráfico 42 - VAN versus vida útil del proyecto de precalentar agua de calderas usando en el caso base de
combustible gas natural
256
Universidad Técnica Federico Santa María
Departamento de Industrias
12,00
10,00
8,00
6,00
CPP
VAN (M USD)
4,00
CTV
2,00
CCP
0,00
MAC
-2,00 0 5 10 15 20 25
CCF
-4,00
-6,00
-8,00
Vida útil del proyecto (años)
Gráfico 43 - VAN versus vida útil del proyecto de precalentar agua de calderas usando en el caso base de
combustible diesel
Como se puede observar, en ningún caso al evaluar con combustible gas natural
se logró obtener un VAN positivo, a diferencia de las evaluaciones con diesel,
que a partir del año 12 se obtiene un VAN positivo para todas las tecnologías.
257
Universidad Técnica Federico Santa María
Departamento de Industrias
11. Conclusiones
11.1. De los Balances de energía
De todos los proyectos evaluados, queda en evidencia que entre colectores solares de
placa plana (CPP), tubo de vacío (CTV), cilíndrico parabólico (CCP), concentración con
reflector estacionario (MAC) y fresnel de concentración lineal, los CCP son la
tecnología que más se ajusta a los proyectos de minería. La principal causa de esto
radica en el costo de inversión de las otras tecnologías, que sigue siendo muy elevado
con respecto a los colectores solares planos, ya que la mayoría de las tecnologías de
concentración requieren una gran cantidad de partes móviles y estructuras sofisticadas,
para poder realizar un seguimiento en un eje del sol y a la vez sostener una gran cantidad
de espejos concentradores de radiación solar.
258
Universidad Técnica Federico Santa María
Departamento de Industrias
La rentabilidad de los proyectos de esta naturaleza, dependerá fuertemente de las
exigencias, objetivos y tolerancia al riesgo de los inversionistas, ya que los proyectos
evaluados pueden pasar fácilmente de rentables a no rentables solo cambiando el tipo de
combustible a emplear como caso base o tan simple como cambiar la vida útil del
proyecto.
259
Universidad Técnica Federico Santa María
Departamento de Industrias
12. Bibliografía
[3] COCHILCO, «Informe tendencias del mercado del cobre: Julio - Spetiembre de
2013.,» 2013. [En línea]. Available: www.cochilco.cl.
[5] M. Rommel, «Process Heat Collectors - State of the Art within Task 33/IV: Solar
Heat for Industrial Processes,» Werner Weiss, Alemania, 2008.
260
Universidad Técnica Federico Santa María
Departamento de Industrias
[9] CODELCO, «Radomiro Tomic Fase II,» 2011. [En línea]. Available:
http://www.codelco.com/explotacion-sulfuros-radomiro-tomic-fase-
ii/prontus_codelco/2011-07-06/132737.html.
[10] CODELCO, «Codelco Educa: Celdas electrolíticas,» 2014. [En línea]. Available:
https://www.codelcoeduca.cl/procesos_productivos/tecnicos_electroobtencion_cel
daselectroliticas.asp.
[12] A. Picard, «Revised formula for the density of moist air (CIPM-2007),» 2008.
[13] Incropera, Fundamentals of Heat and Mass Transfer, 6th ed., 2006.
[16] J. Bustos, «Análisis de la huella de carbono en una empresa minera del cobre en
Chile,» 2011.
261
Universidad Técnica Federico Santa María
Departamento de Industrias
[19] IREC, «Suitable Solar Thermal Collectors for Large Scale Industrial
Applications,» Catalunia, 2012.
[21] W. Green, Perry's Chemical Engineers' Handbook, 8th ed., The McGraw-Hill
Companies, Inc., 2008.
[24] CNE, «Comisión Nacional de Energía de Chile,» 2014. [En línea]. Available:
www.cne.gov.cl.
262
Universidad Técnica Federico Santa María
Departamento de Industrias
[27] Investing.com, «US 30 YR T-Bond Futuros - Dec 14 (USZ4),» 2014. [En línea].
Available: http://es.investing.com/rates-bonds/us-30-yr-t-bond.
263
Universidad Técnica Federico Santa María
Departamento de Industrias
13. Anexos
En este capítulo se describen las fórmulas y ecuaciones utilizadas para desarrollar tanto
los balances de energía de los procesos como cálculos y optimizaciones de la evaluación
de proyectos con el software RETScreen.
También se presentan las tablas con los resultados de todos los colectores solares para el
proyecto mejor evaluado.
( ∫ ) ∫ ⃑
∫ ∫ ⃑
∫ ∑ ∑
264
Universidad Técnica Federico Santa María
Departamento de Industrias
Por lo que para que el sistema esté perfectamente balanceado se deberá cumplir que:
∑ ∑
| |
∑
( ∫ ) ∫ ⃑
Y la energía para los flujos másicos “e” viene dada significativamente por:
265
Universidad Técnica Federico Santa María
Departamento de Industrias
( ∫( ) )
( ∫ ) ∑ ̇ ∑ ̇
Para la determinación del calor transportado por cada flujo, se consideró la temperatura
estándar de referencia igual a 25°C o 298,15 K, como se muestra a continuación:
∑ ̇ ∑ ( ) ̇
266
Universidad Técnica Federico Santa María
Departamento de Industrias
| |
∑ ( ) ̇
267
Universidad Técnica Federico Santa María
Departamento de Industrias
13.2. Ecuaciones de Transferencia de Calor
13.2.1. Convección
Las pérdidas de energía por convección vienen dadas por la transferencia de calor que
experimenta la superficie de un cuerpo hacia un fluido que circula por ella. La expresión
que representa esta transferencia de calor está dada por:
Donde
( )
13.2.2. Conducción
La conducción se puede definir como el flujo de calor que pasa a través de un material
con un espesor “L” y con una superficie “A”. Asumiendo un estado estacionario,
transferencia de calor unidimensional, que no existe generación de calor en el material y
que las temperaturas se encuentran en función de solo una dirección, la ecuación que
describe el flujo de calor por conducción vendrá dada por:
268
Universidad Técnica Federico Santa María
Departamento de Industrias
conducción respectivos. Por ejemplo si un fluido transfiere calor a otro fluido por una
placa que los separa (Figura 45);
Figura 45 - Transferencia de calor entre dos fluidos separados por una placa
Dónde:
13.2.4. Radiación
269
Universidad Técnica Federico Santa María
Departamento de Industrias
13.3. Ecuaciones de Transferencia de Masa
( )
Esta misma expresión si se multiplica por el peso molecular de los compuestos, se puede
dejar en función de sus densidades, permitiendo determinar directamente la transferencia
de masa, como se muestra a continuación:
̇ ( )
270
Universidad Técnica Federico Santa María
Departamento de Industrias
Dónde:
̇ ( )
Dónde:
̇ ( )
( )
( )
( )
271
Universidad Técnica Federico Santa María
Departamento de Industrias
Dónde:
( )
( )
La densidad del vapor saturado del aire se determinó mediante la siguiente función [12]:
[ ( )]
[ ] ( )
Donde:
[ ]
i [ ]
272
Universidad Técnica Federico Santa María
Departamento de Industrias
[ ]
[ ]
[ ]
[ ]
[ ]
[ ]
[ ]
[ ]
[ ]
[ ]
[ ]
273
Universidad Técnica Federico Santa María
Departamento de Industrias
[ ]
[ ]
[ ]
274
Universidad Técnica Federico Santa María
Departamento de Industrias
13.4. Modelos para Coeficientes de Convección
Convección Libre
Convección Forzada
Mezcla de Convección Libre y Forzada
275
Universidad Técnica Federico Santa María
Departamento de Industrias
13.4.2. Analogía de capas límite
De esta manera luego se pueden relacionar los números de Nusselt y Sherwood para
llegar a siguiente la expresión:
276
Universidad Técnica Federico Santa María
Departamento de Industrias
13.4.3. Coeficientes de transferencia de Calor
En esta sección se presentan los modelos para la determinación del número de Nusselt
según el tipo de convección.
̅̅̅̅ [ ]
277
Universidad Técnica Federico Santa María
Departamento de Industrias
Si se cumple que entonces el flujo es turbulento, con:
̅̅̅̅ ( )
̅̅̅̅ ( )
̅̅̅̅ ( )
278
Universidad Técnica Federico Santa María
Departamento de Industrias
̅̅̅̅ ( )
̅̅̅̅ ( )
̅̅̅̅
̅̅̅̅
279
Universidad Técnica Federico Santa María
Departamento de Industrias
̅̅̅̅
Válido para
Donde el símbolo permite identificar si los flujos se oponen en el sistema y n=7/2 para
placas horizontales.
̅̅̅̅
[ ( ) ]
{ }
280
Universidad Técnica Federico Santa María
Departamento de Industrias
13.5. Coeficientes de transferencia de Masa
En esta sección se presentan los modelos para la determinación del número de Sherwood
según el tipo de convección.
Se propone el siguiente modelo para electrodos paralelos (Válido para electrodos con
):
( )
281
Universidad Técnica Federico Santa María
Departamento de Industrias
̅̅̅ ( )
Este modelo describe el proceso de evaporación natural del agua en un estanque [11].
282
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6. Energía Solar
13.6.1.1. Declinación
( )
Dónde:
Corresponde al desplazamiento angular del sol al este u oeste del meridiano local y se
determina como:
Dónde:
283
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.1.3. Radiación extraterrestre e Índice de Claridad
( ( ))
Dónde:
̅
̅
̅
Dónde:
̅ ̅ ̅ ̅ ( ) ̅ ( )
Dónde:
284
Universidad Técnica Federico Santa María
Departamento de Industrias
̅
Si ̅ ̅( ̅ ̅ ̅ )
Si ̅ ̅( ̅ ̅ ̅ )
̅ ̅ ̅
Asumiendo que el cielo es un cuerpo negro perfecto, la radiación emitida por el cielo
puede expresarse como:
Esta radiación puede estimarse también como la suma de la radiación emitida por el
cielo con nubes ( ) y la radiación emitida por el cielo completamente despejado
( ), considerando la fracción del cielo cubierto con nubes (“c”):
285
Universidad Técnica Federico Santa María
Departamento de Industrias
Con
̅
̅
̅ ̅ ̅ ̅ ̅
̅ ̅
{
286
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.2. Colectores Solares
Dónde:
Para estos colectores solares tanto como son independientes del viento.
287
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.2.1. Pérdidas por cañerías y estanques
288
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.3. Método: Carta “f”
̅̅̅
Dónde
̅̅̅
( )
289
Universidad Técnica Federico Santa María
Departamento de Industrias
̇
[ ( )( )]
̇
Dónde
Las capacidades calóricas que circulan por los intercambiadores de calor “Colector –
Reservorio” son:
Por lo que
Dónde
290
Universidad Técnica Federico Santa María
Departamento de Industrias
Si el intercambiador de calor es de flujos cruzados:
( )
Válida para
Con
291
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.4. Método de la Utilizabilidad
13.6.4.1. Principio
̅̅̅
Dónde:
∑ ∑ ̅̅̅
Donde “ + “ denota que solo se consideran los valores positivos entre los soportes.
∑ ∑
̅
̅̅̅̅
̅̅̅ ̅̅̅̅ ̅
̅
̅
292
Universidad Técnica Federico Santa María
Departamento de Industrias
La proporción entre la radiación de la hora de mediodía de la superficie inclinada y la
horizontal para un día promedio del mes será:
( ) ( )( ) ( )
Donde
̅
̅
̅ [( )( ̅ ̅ )]
Dónde:
̅ ̅
̅ ̅
̅ ̅
293
Universidad Técnica Federico Santa María
Departamento de Industrias
13.6.5. Resultados Precalentamiento de Agua de Calderas – Diesel
Modelo TA (°C) nC AT.BRUTA (m2) ƒSOLAR (%) QCOL (MWh/año) IT (M USD) GEIRED (tCO2) LEC (USD/MWh) TIR (%) RCAP (año) VAN (M USD)
XL USB 76,6 4.457 21.064 8,1% 10.868 5,2 12.232 51 -23,8% 2 1,57
HT-A 28/10 76,2 724 12.770 8,1% 10.928 5,8 12.282 53 -33,0% 3 1,12
APSE-30 71,6 2.288 13.500 8,4% 11.328 6,0 12.352 54 -35,3% 3 1,02
DF30 - 300 71,6 1.831 11.167 8,3% 11.193 6,3 12.321 54 -40,8% 3 0,78
GJ140A 76,4 710 16.001 8,1% 10.918 6,5 12.259 55 -49,8% 3 0,44
KAIROS XP 2,5 V 76,2 4.228 18.113 8,2% 11.027 6,6 12.252 56 -53,1% 3 0,35
HP30-300 71,6 2.252 14.188 8,5% 11.373 7,2 12.357 57 - 3 0,01
HT 75,6 677 12.182 8,3% 11.138 7,3 12.323 57 - 3 -0,11
GS 76,9 815 14.668 7,9% 10.679 7,3 12.237 58 - 3 -0,26
XL10 76,6 1.098 13.236 8,0% 10.760 7,9 12.261 59 - 4 -0,77
AL10 76,5 1.096 15.229 8,0% 10.819 8,2 12.254 60 - 4 -0,99
CPC 45 71,6 1.417 9.083 8,3% 11.116 8,7 12.349 61 - 4 -1,29
gigaSol S 76,6 1.014 16.969 8,0% 10.812 9,1 12.248 63 - 4 -1,79
FKC-2S 76,4 4.599 16.348 8,1% 10.918 9,2 12.262 63 - 4 -1,87
Polytrough 1800 71,6 156 14.932 8,3% 11.114 9,3 12.243 63 - 4 -1,84
SKS 4,0 76,4 4.506 15.160 8,1% 10.872 9,4 12.253 64 - 4 -2,02
FA(3/6) 76,0 610 16.807 8,2% 11.074 9,6 12.281 64 - 5 -2,17
HP250 71,6 1.805 10.764 8,3% 11.210 9,9 12.339 65 - 5 -2,34
CCStaR 2V 71,6 147 16.642 8,3% 11.108 10,5 12.201 67 - 5 -2,92
VKF 125 77,2 4.940 18.601 7,9% 10.619 10,6 12.204 67 - 5 -3,09
SOL250 76,7 4.262 18.186 8,1% 10.838 10,7 12.228 67 - 5 -3,12
DF100 71,6 1.770 10.532 8,3% 11.200 11,3 12.337 69 - 5 -3,46
10 HTF 76,4 944 14.066 8,1% 10.876 11,5 12.265 70 - 5 -3,77
294
Universidad Técnica Federico Santa María
Departamento de Industrias
VKF 145V 78,3 4.807 18.098 7,8% 10.433 11,4 12.169 70 - 6 -3,82
LF-11 71,6 274 21.693 8,2% 11.088 13,4 12.139 76 - 6 -5,47
X10T 71,6 753 21.385 8,4% 11.249 14,9 12.254 80 - 7 -6,74
295
Universidad Técnica Federico Santa María
Departamento de Industrias
FA(3/6) 76,0 610 16.807 8,2% 11.074 9,6 12.281 59 -15,4% 5 2,25
HP250 71,6 1.805 10.764 8,3% 11.210 9,9 12.339 60 -16,3% 5 2,20
CCStaR 2V 71,6 147 16.642 8,3% 11.108 10,5 12.201 61 -22,0% 5 1,60
VKF 125 77,2 4.940 18.601 7,9% 10.619 10,6 12.204 62 -25,3% 5 1,32
SOL250 76,7 4.262 18.186 8,1% 10.838 10,7 12.228 62 -25,3% 5 1,33
DF100 71,6 1.770 10.532 8,3% 11.200 11,3 12.337 63 -28,8% 5 1,15
10 HTF 76,4 944 14.066 8,1% 10.876 11,5 12.265 63 -41,7% 5 0,76
VKF 145V 78,3 4.807 18.098 7,8% 10.433 11,4 12.169 63 positivo 6 0,63
LF-11 71,6 274 21.693 8,2% 11.088 13,4 12.139 67 negativo 6 -0,83
X10T 71,6 753 21.385 8,4% 11.249 14,9 12.254 70 negativo 7 -2,00
296