Diseño Inverso de Chimeneas Industriales
Diseño Inverso de Chimeneas Industriales
Diseño Inverso de Chimeneas Industriales
Nomenclatura
Caracteres latinos:
AG(s) algoritmo(s) gentico(s)
CF
CT
CV
DA
deterioro ambiental
inters anual
presin [Pa]
PM
Q&
Caracteres griegos:
h
Introduccin
1.1
El Problema Inverso
La definicin del problema inverso admite otras formulaciones alternativas, que estn muy
relacionadas con las anteriores. stas incluyen:
Encontrar las constantes de un modelo que mejor se ajusten a los resultados de las
observaciones.
De la anterior descripcin, se puede concluir que la mayora de las resoluciones directas que se
hacen de un problema de ingeniera son en realidad uno de los pasos de un problema inverso que
debe ser resuelto.
Por ejemplo, cuando se calcula un flujo determinado de una instalacin, es frecuente hacerlo con
el objetivo de mejorar el proceso de un producto. A raz de los resultados del clculo, el
ingeniero, la mayora de las veces, introduce cambios en el diseo o en las condiciones de
operacin con la esperanza de hacer el proceso ms eficiente. Este ingeniero est, sin saberlo,
intentando resolver un problema inverso.
Se puede acudir a varias tcnicas para resolver el problema inverso, entre las que se encuentran
gran parte de los mtodos de optimizacin existentes. En este caso se emplearn como base para
la resolucin de este problema los Algoritmos Genticos (AGs), ya que presentan la ventaja de
poderse desacoplar del problema a resolver y de ser muy robustos, con lo cual es posible su
aplicacin a una gran variedad de problemas.
1.1.1 Problema inverso, Algoritmos Genticos y Combustin
Debido a que el problema inverso, en cualquiera de las formas presentadas anteriormente, es
muy relevante en la generacin de electricidad a partir de combustibles fsiles, este trabajo
pretende dedicarse al estudio de su resolucin, especialmente a la problemtica entorno a la
generacin de contaminantes en la combustin.
Los Algoritmos Genticos, desde su aparicin en la dcada de los 70, desarrollados por John
Holland (1975), han sido empleados hasta ahora, en multitud de campos y aplicaciones, entre las
que destacan: la optimizacin, la programacin automtica, el aprendizaje mquina, la
economa, el estudio de sistemas sociales o la asignacin de recursos. Su extensin a todas estas
reas de conocimiento se debe a la robustez que presentan en su estructura ms sencilla y a su
correcta adaptacin para resolver problemas reales, que ha sido posible gracias al esfuerzo
realizado por varios grupos de investigacin en el desarrollo de tcnicas avanzadas para AGs.
Un Algoritmo Gentico (AG) consigue resolver problemas de optimizacin de una manera
rpida y eficaz, mediante la emulacin del proceso natural de supervivencia de los mejor
adaptados. Para ello, genera poblaciones cuyos individuos son las posibles soluciones a un
problema y las hace evolucionar, consiguiendo, en el transcurso de un cierto nmero de
3
generaciones, conocer la solucin ptima para dicho problema. Como no requiere la derivacin
de ninguna funcin, puede dar solucin a problemas irresolubles con otros mtodos. Adems,
los AGs pueden ser masivamente paralelizados, reduciendo enormemente el tiempo de
computacin frente a otros mtodos de optimizacin.
La utilizacin de combustibles fsiles en calderas produce una serie de contaminantes (CO2,
CO, SO2 y NOx entre otros), que tras haberse realizado los procesos correspondientes de
limpieza de gases, deben ser expulsados a la atmsfera garantizando la mxima seguridad para
las personas y el medio ambiente circundante. Los dispositivos empleados para expulsar los
gases son las chimeneas, y de su correcto diseo depende la calidad del servicio.
En este trabajo se propone una metodologa para el diseo de chimeneas que consigue
minimizar el impacto ambiental que stas suponen, debido a la emisin de gases contaminantes,
a la vez que se minimizan los costes de construccin y explotacin.
1.2
Diseo de chimeneas
El diseo de una chimenea industrial radica en calcular su dimetro interno (D) y su altura (H).
Por otra parte, debe comprobarse que la chimenea diseada crea suficiente tiro.
Antes de comenzar con la discusin acerca de cmo afrontar el diseo de chimeneas, se definen
los siguientes conceptos:
La misin de una chimenea en una industria que emite gases a la atmsfera es crear el tiro
necesario para conseguir que la emisin de dichos gases se haga de tal forma que la
concentracin de contaminantes en inmisin no sobrepase en ningn caso los lmites permitidos
por la legislacin. Esto supone que se tenga que conocer, antes de la construccin de la
chimenea, que concentraciones se van a dar en inmisin a partir de ciertos valores de emisin
que es posible calcular en el proceso de diseo del proceso industrial correspondiente.
Para disear una determinada chimenea se deben conocen los lmites para cada uno de los
contaminante producidos, as como los valores de concentraciones de emisin para cada uno de
ellos. Se necesita algn medio para relacionar estos valores, de tal forma que se pueda asegurar
4
que la chimenea diseada es apta, o por el contrario, genera concentraciones en inmisin por
encima de los valores establecidos en la legislacin. El mtodo buscado debe cumplir las
siguientes caractersticas:
Ser universal, es decir, que sea posible aplicarlo a todo tipo de chimeneas, sin restriccin
de dimensiones ni de situacin geogrfica.
Ser fiable, de tal forma que garantice un valor correcto en inmisin a partir de un valor
de emisin dado.
Objetivos
Realizar el diseo de tal forma que se consiga cumplir simultneamente los dos
siguientes propsitos:
o Minimizar el impacto ambiental producido por los gases emitidos.
o Minimizar el coste econmico de la construccin de chimenea y de su posterior
explotacin.
El Algoritmo Gentico
El AG empleado constar de un AG bsico ms una serie de tcnicas avanzadas que mejoran sus
capacidades hacindolo ms rpido. Adicionalmente se le incluyen una serie de tcnicas
avanzadas, como son la capacidad de realizar anlisis multiobjetivo y la paralelizacin del
clculo.
3.1
AG bsico
Ciega, es decir, no disponen de ningn conocimiento especfico del problema (no hay
ninguna funcin analtica que lo defina a priori) de manera que la bsqueda se basa
exclusivamente en los valores de la funcin objetivo.
Codificada, puesto que no trabajan directamente sobre el dominio del problema, sino
sobre representaciones de sus elementos.
Estocstica referida tanto a las fases de seleccin como a las de transformacin. Esto
proporciona control sobre la explotacin y exploracin del espacio de bsqueda.
Generar
poblacin incial
(n individuos)
Funcin de
coste (Fitness)
Seleccin
4
generacin =
= generacin +1
Reproduccin
(Operadores
genticos)
Reemplazo
n
NO
6
Criterio de
terminacin?
S
Individuo ptimo
10
La seleccin de los individuos que pasarn de una generacin a otra debe hacerse considerando
la aptitud del individuo y la aptitud promedio de la poblacin actual. De modo que los
individuos con una aptitud mayor que el promedio tienen una alta probabilidad de sobrevivir y
los individuos con una aptitud inferior al promedio tienen una alta probabilidad de morir.
Los criterios con que se seleccionan los criadores no necesariamente han de ser los mismos con
que se seleccionan los supervivientes, de ah la necesidad de especificarlos por separado.
[6] Iteracin
Se deben concretar las condiciones con las que se considera que el AG ha dado con una solucin
aceptable o, en su defecto, ha fracasado en la bsqueda y no tiene sentido continuar.
3.2
Anlisis Multiobjetivo
consideran del mismo nicho y se rebaja su aptitud para favorecer la diversidad evitando su
reproduccin.
3.3
Paralelizacin
12
AG en serie
Poblacin/
Variables AG
Poblacin /
Seleccin
Entrecruzamiento
Mutacin
Nueva
generacin
Fitness
CPU 1
Evaluacin
del problema
Scripts de
distribucin
(CPU 0)
CPU i
Individuo
Variables AG
Evaluacin
del problema
Individuo
Fitness
Scripts de
recoleccin
(CPU 0)
CPU n
Evaluacin
del problema
Coste
H, D
Diseo
chimenea
H, D, v
ISCST
Evala
inmisiones
Valor
funcin
objetivo
DA
De las tres posibles vas planteadas para el diseo de chimeneas industriales (seccin 1.2), tan
slo la tercera va cumple correctamente las dos caractersticas que permiten aplicar algoritmos
genticos y realizar el diseo inverso deseado.
En este trabajo se ha optado, de entre los modelos de dispersin de contaminantes existentes y
atendiendo a las consideraciones de la Tabla 1 (Zannetti, 1993), por hacer uso de un programa
de dispersin de contaminantes de penacho gaussiano.
El programa de penacho gaussiano a emplear es Industrial Source Complex (ISC) (US EPA,
1995).
13
Escala de
procesos
atmosfrico
Microescala
Mesoescala
Macroescala
s
Escala de
fenmenos
de
Local
Local -
Regional -
Regional
Continental
Global
dispersin
Modelo
Sobreelevac
in
Gaussiano
Semiempri
co
Euleriano
Lagrangian
o
Qumico
1,2,4
1,2,4
1,2
1,2,3,4
1,2,4
1,2,4
2,3,4
2,4
2,4
(1,2,)4
2,3,4
2,4
Recepcin
2,4
Estocsticos
2,4
2,4
2,4
Diseo de chimeneas
El diagrama de bloques de la Figura 2 muestra la relacin entre los distintos mdulos de clculo
que intervienen en el proceso de diseo de la chimenea.
4.2.1 Diseo de la chimenea
El dimetro y la altura son parmetros necesarios para disear una chimenea, pero se precisa
conocer tambin la energa empleada para inducir el tiro y la velocidad media de los gases
dentro de ella.
14
El flujo a analizar es compresible; se considera que la chimenea es adiabtica y que los gases
pueden tratarse como ideales.
Tomando los puntos de referencia de la Figura 3 y realizando un balance de fuerzas a los gases
dentro de la chimenea, se obtiene la ecuacin 1.
P1 = Patm + r gH + DPf 1 2 - ( r a - r ) gH
(1)
Donde:
Patm
r gH
DPf 1 2
( r a - r ) gH
2
x
(x,y,z)
(x,0,0)
Gases de
chimenea
(x,y,0)
y
Hs
(a)
(b)
Figura 3 (a) Chimenea y ventilador para inducir el tiro. (b) Representacin del penacho gaussiano
calculado por ISC.
La prdida de carga entre el punto 1 y el 2 de la chimenea (Sandler, 1985) se calcula como suma
de dos trminos (Ec. 2), uno debido a la energa cintica y otro debido a la friccin.
DPf 12 = DPK + DPF
(2)
15
P
1
DPK = r 1 v12 1
P2
2
DPF = K FA f
k
- 1
(3)
L r 12 v12
D 2r
(4)
El factor de friccin para las prdidas de carga se calcula a partir de la ecuacin de Coolebrok
(Ec. 5). El sistema de ecuaciones resultante se resuelve de forma iterativa.
e
2,51
= -2 log
+
3
,
7
D
f
Re f
(5)
Los parmetros de clculo empleado son los mostrados en la Tabla 2, y el clculo de las
propiedades del flujo en los puntos 1 y 2 se realizan de acuerdo a las expresiones de la Tabla 3.
Valor
P0
Tamb
Patm
PM
[Pa]
[k]
[Pa]
[g/mol]
298
101325
29,92
10132
5
R
[J/mol
k]
m&
[kg/s] =Cp/Cv
8,31
140
1,4
0,03
T [k]
400
2
P
T2 = T1 2
P1
k -1
k
16
Ecuacin (1)
P [Pa]
3
r1 =
[kg/m ]
P1 PM
RT1
Patm
k
P
v2 = v1 1
P2
4m& RT1
v1 =
P1 PMpD 2
v [m/s]
Tabla 3 - Valores de T, P,
P
r 2 = r 1 2
P1
(6)
(7)
El coste de una chimenea puede variar dependiendo del material de construccin, que a su vez
depende de la temperatura de los gases de salida. La ecuacin 7 puede ser sustituida por aquella
ms ajustada al tipo de chimenea que se vaya a construir.
4.2.2.2 Costes variables
Los costes variables se calculan mediante una estimacin de la energa necesaria para
proporcionar el tiro inducido necesario para la correcta expulsin de los gases por la chimenea.
En el apartado de Diseo de la chimenea, se calcula la presin en el punto 1. La energa
necesaria para tener esa presin es igual a la proporcionada por el ventilador menos la presin en
el punto cero.
17
Los costes variables para n aos de funcionamiento (en Euros), con crecimiento de precios i y
con un coste C de la electricidad (Tabla 4), es:
CV = 8,76C
( P1 - P0 )Q 1 - (1 + i ) - n
h
i
C
[Euros/kW
[aos]
h]
Valor
0,02
3%
(8)
0,9
15
18
Respecto a la topografa: malla del terreno, con la altura del terreno en cada nodo de la
malla.
Los resultados del programa son:
Resumen de todos los parmetros de control empleados para la ejecucin del programa.
Tablas con los mayores valores en inmisin calculados para cada receptor.
La ecuacin de difusin gaussiana relaciona los niveles de inmisin en un punto con la cantidad
de contaminante vertido a la atmsfera desde la fuente emisora, teniendo en cuenta las
condiciones de emisin de los gases y las caractersticas atmosfricas y topogrficas del medio
receptor.
4.2.4 Estimacin del Impacto Ambiental
Una vez aplicado ISC, los resultados de concentracin en inmisin en cada uno de los receptores
son conocidos. Es necesario establecer un criterio para transformar toda la informacin obtenida
con el modelo y reflejarla en uno o varios parmetros homogneos y adimensionales que
permitan comparar unos casos frente a otros.
19
Los criterios que se muestran en la siguiente lista cumplen las caractersticas necesarias, luego
podran ser tiles para realizar esta evaluacin:
Los ndices de calidad del aire se pueden emplear cuando se traten simultneamente las
emisiones de varios contaminantes.
El criterio que se propone usar es el de ndices de Calidad Ambiental (CA), en particular los
ndices propuestos por el instituto Battelle y que son aplicados por la EPA.
4.2.4.1 Indice de calidad ambiental Battelle-EPA
El uso de estos ndices permiten obtener un parmetro, que se denominar Deterioro Ambiental
(DA), con el que poder comparar concentraciones de distintos contaminantes medidos para
distintos promedios de tiempo. En este caso slo se trabaja con concentraciones de NOx, por lo
tanto slo se emplea para comparar concentraciones en distintas escalas temporales.
El valor de DA queda delimitado entre 0 y 10, tanto para los promedio horarios como para los
promedios anuales, que son, los promedios cuyos valores de concentracin estn definidos en la
legislacin.
Basndose en los valores mostrados en la Tabla 5 se obtiene la Ecuacin 9. Para los dos casos,
un valor de DA correspondiente a 7 es el mximo permitido por la legislacin.
Los mismos conceptos presentados para las concentraciones horarias se aplican a los valores de
concentracin promedio anual, obtenindose la expresin (10).
Promed Unidad
NO2
io
1 hr.
es
g/m3 02 45 6
5
8 11 15 20 27 40 40
0*
20
Escala adimensional
10
Tabla 5 - Escala para transformar los valores de inmisin de NO2. * Valor lmite de NO2 para
promedio horario.
2]
2]
(9)
(10)
Escenario
La central trmica propuesta para el escenario, produce energa elctrica a partir de un ciclo
combinado, empleando como nico combustible gas natural.
La tecnologa de ciclo combinado consiste en el aprovechamiento de los gases residuales de una
turbina de gas, los cuales se encuentran a temperaturas elevadas, como fuente de energa para un
ciclo de vapor. Las ventajas de un ciclo combinado respecto a otros tipos de ciclos trmico para
la produccin de energa elctrica es su elevado rendimiento, pudiendo considerarse valores
entorno al 55% en condiciones ISO.
En la central trmica de turbinas de gas existen 4 turbinas de gas y 4 calderas de recuperacin,
que producen vapor para dos turbinas de vapor, produciendo un potencia total de 250 MW al
100% de carga. A la salida de estas calderas de recuperacin de calor se sitan las 4 chimeneas
que son objeto del proceso de diseo que se quiere llevar a cabo.
21
Legislacin
22
350
50(3),(4),(5)
Tabla 6 Niveles de emisin permitidos para NOx. (1) Niveles de emisin para NOx se
refiere al 3 % de O2 en base seca de los gases residuales, calculando el volumen de
los gases en condiciones normales (101.3 KPa y 0 C) Niveles de emisin para NOx,
(2) Concentracin referida al 15 % de O2 en los gases residuales en base seca,
calculando el volumen de los gases en condiciones normales (101.3 KPa y 0 C). (3)
Los valores lmites se aplican nicamente por encima de una carga del 70 %. (4) Para
turbinas de gas utilizadas en ciclos combinados cuyo rendimiento elctrico global
medio anual sea superior al 55 %, el lmite se establece en 75 mg/Nm3, medido en las
mismas condiciones (referido al 15 % de O2 en base seca).
5.2.2 Niveles de inmisin
Los lmites de inmisin de los xidos de nitrgeno se fijaron en el R.D.833/1975 y
posteriormente se actualizaron en el R.D.717/1987. En el ao 2002 se traspuso la Directiva
99/30/CE, la Directiva 2000/69/CE y la Decisin de la Comisin del 17 de Octubre de 2001 en
el R.D.1073/2002, que actualmente es la legislacin vigente referida a los niveles de inmisin de
las principales sustancias contaminantes.
En la tabla 12 se presentan los valores de los niveles de inmisin para los xidos de nitrgeno
que marcan el Real Decreto 1073/2002.
R.D. 1073/2002
Periodo de
promedio
Valor lmite
Margen de
tolerancia
Fecha de
cumplimiento de
valor lmite
50 % a la entrada
200 mg/m3 de NO2,
Valor lmite
horario para la
proteccin de la
salud humana
superarse ms de 18
ocasiones por ao
civil
en vigor de la
directiva, con una
reduccin lineal a
partir del 1 de
1 de enero de
2010
enero de 2001,
hasta alcanzar el 0
% el 1 de enero de
23
2010.
50 % a la entrada
en vigor de la
directiva, con una
Valor lmite
anual para la
proteccin de la
reduccin lineal a
40 mg/m3 de NO2
1 ao civil
1 de enero de
partir del 1 de
2010
enero de 2001,
salud humana
hasta alcanzar el 0
% el 1 de enero de
2010.
Valor lmite
anual para la
proteccin de la
30 mg/m3 de NOx
1 ao civil
19 de Junio de
Ninguno
2001
vegetacin
Tabla 7 - Lmites de inmisin para NOx en el R.D. 1073/2002. Los valores que aparecen en esta
tabla se miden en condiciones estndar (101,3 KPa y 20 C)
5.3
La informacin necesaria para ejecutar ISC se describe a continuacin. Los datos son
introducidos mediante ficheros externos que incluyen toda la informacin necesaria.
5.3.1 Fuentes de contaminantes
Como se ha indicado en la descripcin del escenario se consideran cuatro fuentes puntuales
cuyos parmetros se dan en la Tabla 8.
Tipo de
fuente
[gNOx/s]
H [m]
T [K]
V [m/s]
D [m]
Para cada
una de las
Puntual
4
chimenea
13,53
Variable
del AG
Calculada
400
(ver Tabla
3)
Variable
del AG
5.3.2 Topografa
El terreno se ha modelado sobre la base del modelo digital del terreno MDT200 del Instituto
Geogrfico Nacional (IGN). Cubre una extensin de 1.600 km2. Este modelo digital proporciona
informacin para simular el terreno con un paso de malla de 200 m. Se considera que la esquina
inferior izquierda, del cuadrado de 40 x 40 km simulado, es el origen del sistema de coordenadas
de la malla. En la Figura 5 se muestra un plano del terreno incluido en la simulacin.
5.3.3 Receptores
La evaluacin de la concentracin en inmisin se realiza en puntos discretos sobre el terreno
simulado. En este caso se han empleado receptores dispuestos en la misma disposicin que la
malla del terreno pero con un paso de 400 m. Disponiendo de unos 6 receptores por cada
kilmetro cuadrado.
5.3.4 Meteorologa
Los datos meteorolgicos se obtienen a partir del INM. En la Figura 4 se presenta la rosa de los
vientos elaborada a partir de los datos disponibles para el ao 2001.
5.3.5 Calidad del aire
El entorno propuesto para la central es tpicamente rural, alejado de grandes vas de circulacin.
Por ello, la contaminacin de fondo derivada de focos mviles (automocin) o de actividad
industrial debe ser, muy baja. Se fijan los siguientes valores de concentracin de fondo para
NOx:
Concentracin de
Tiempo promedio
NOx
[ g/Nm3]
1 hora
1 ao
25
Ao 2001
25
NNW
Calmas
4,7%
NNE
20
NW
NE
15
WNW
10
ENE
5
W
WSW
ESE
SW
SE
SSW
SSE
S
Figura 4 - Rosa de los vientos anual correspondiente a datos de 2001 de la estacin meteorolgica
representativa de la zona de estudio.
340
220
40
30
0
30
22 0
260
300
30 0
30
2
26 0 2
22 0
<Ce ntral>
o
1 80
180
1 80
22
0
18
26
300
26
22
220
180
20
30
22
38
380
30
26
260
300
30
0
34
0
34 38 0
10
34
10
0
22
300
3 40
38
42
340
30
42
20
30
40
5.4
26
Parmetros
Valor
Generaciones
23
24
Elitismo
Seleccin
Torneo
Niche
No
Paralelizacin
Multiobjetivo
No
---
share
Mutaciones
Jump mutation
0,01
Creep mutation
No
Puntual
0,65
Uniforme
No
Cruce
Caracterstica
Valor
Lmites a variables
H [m]
10 350
D [m]
1,5 15
Precisin
H [m]
11 bits
D [m]
0,1660
9 bits
0,02637
27
1
max
CT + 10 max( DA1max
h , DAanual )
10 5
(11)
Se explica en secciones posteriores como se calculan, tanto los costes como los valores de
Deterioro Ambiental (DA).
6
6.1
Resultados
Resultados del AG
La Figura 6 muestra la evolucin del fitness (Ec.11) medio y mximo a lo largo de las
generaciones. Los valores del fitness son negativos porque el objetivo es minimizar el minimizar
el coste y el deterioro ambiental simultneamente.
Se puede observar que el comportamiento del GA es bueno, ya que los valores medios y
mximos convergen rpidamente.
El tiempo de ejecucin de un caso de ISCST junto con el diseo de la chimenea es de 25
minutos en un procesador Pentium II a 500 MHz. Se ha realizado la paralelizacin empleando 8
CPUs, con lo cual el tiempo total para completar las 23 generaciones ha sido de algo menos de
29 horas.
28
-50
-100
-150
Fitness
-200
-250
-300
-350
-400
-450
0
10
15
20
25
Ge ne raciones
Figura 6 - Evolucin del valor de la funcin objetivo, tanto para el mejor individuo de cada
generacin como para la media.
6.2
Chimenea ptima
El diseo del caso ptimo supone obtener la altura y el dimetro de las cuatro chimeneas (las
cuatro son iguales) que minimizan, tanto el coste ambiental que implica realizar la emisin de
NOx como el coste econmico de su instalacin y operacin.
La Tabla 12 muestra las dimensiones ptimas de las chimeneas encontradas por el GA; as
como, la velocidad de los gases dentro de la chimenea.
H [m]
D [m]
v [m/s]
47,54
2,47
31,0
Valor
Tabla 12- Altura, dimetro y velocidad de los gases de las chimeneas ptimas diseadas.
La Tabla 13 muestra los costes que implica instalar y operar con las cuatro chimeneas ptimas.
CV [miles
CF[miles
CT[miles
230,8
454,6
685,4
Tabla 13- Costes fijos, variables y totales de la suma de las cuatro chimeneas del caso ptimo.
29
30
0
34
260
22
0
30
3 00
260
30
18
220
22
22
Central
Central
22
18
220
30
2 20
30
30
26
18
26
300
34 0
30
380
0
38
34 0
260
300
34 0
380
340
30
0
38
34
300
26 0
22 0
26 0
220
300
22
180
18
2 60
Central
Central
30
220
300
220
2 60
30 0
34
260
30
38
38
30 0
300
26 0220
0
30
0
34
30
30
0
34
30
38 0
30
380
34
4 60
microg/Nm3
149
139
129
119
109
99
89
79
69
59
50
40
30
20
10
microg/Nm3
2.99
2.87
2.76
2.64
2.52
2.41
2.29
2.17
2.06
1.94
1.83
1.71
1.59
1.48
1.36
1.24
1.13
1.01
0.90
0.78
0.66
0.55
0.43
0.32
0.20
Figura 7. (a) Concentraciones promedio horarias mximas junto con el terreno simulado. El mximo
permitido por la concentracin es de 200 microgramos de NOx por metro cbico
normal. Se aprecia claramente en la figura como los contornos de concentraciones
mximas horarias se ajustan a la orografa del terreno. (b) Concentraciones promedio
anuales mximas junto con el terreno simulado. El mximo permitido por la legislacin
es 40 microgramos de NOx por metro cbico normal. Se aprecia en la figura como las
31
6.3
Anlisis multiobjetivo
Parmetros
Valor
Generaciones
40
28
Elitismo
Seleccin
Torneo
Niche
Paralelizacin
Multiobjetivo
S
0,8
share
Mutaciones
Jump mutation
0,01
Creep mutation
0,02
Puntual
No
Uniforme
Cruce
32
La chimenea diseada anteriormente, en base a una funcin objetivo, se situara en el punto con
coste 0,68 y DA 6,3. Esto implica que se ha diseado una chimenea en la que predomina el
objetivo de bajo coste.
Hay que tener en cuenta, en el anlisis de la Figura 7, que el clculo del Deterioro Ambiental
realizado es conservativo, es decir, eligiendo un diseo de chimenea en el que el DA sea 6,99
(menor que 7), ya se est cumpliendo la legislacin, por lo tanto no es necesario invertir ms
dinero en chimeneas ms altas. Ahora bien, siempre sera correcto invertir ms, con la idea de
causar el menor impacto ambiental posible.
En la Figura 7.b se representan mediante lnea discontinua los diseos que pertenecen al ptimo
de Pareto, acotando el espacio de representacin entre 4,3 y 7 para el DA. Del anlisis de esta
grfica, junto con la anterior, se aprecia que llega un punto, en el que es necesario realizar
inversiones muy altas para poder reducir levemente el impacto ambiental producido por los
gases emitidos.
Otra representacin que ofrece informacin til a la hora de tomar la decisin sobre qu diseo
escoger se muestra en la Figura 7.c. En ella se puede ver el valor de DA mximo que se obtiene
en funcin de la altura de la chimenea escogida. No se representa el valor de DA en funcin del
dimetro de la chimenea, porque ste no afecta tan decisivamente en la concentracin en
inmisin.
Los puntos en negrita de la figura son los pertenecientes al frente ptimo de Pareto, esto quiere
decir que los puntos situados en la lnea que une todos los diseos del frente proporcionan
informacin sobre el Deterioro Ambiental mnimo que se puede conseguir con una determinada
altura de chimenea.
Al igual que en la representacin del plano objetivo (Figura 7.a) se puede apreciar que llega un
punto, entorno a un valor de deterioro ambiental de 5,3, que para poder disminuir el impacto
producido hay que aumentar mucho la altura de la chimenea. Esto probablemente se deba a la
orografa del terreno, ya que las elevaciones existentes en el entorno de la central, aunque no de
gran altura, s provocan concentraciones elevadas a pesar de aumentar la altura de la chimenea.
33
920
10
720
Coste [miles ]
Coste [Millones de ]
820
620
520
420
320
220
0
3
10
5,8
6,1
6,4
6,7
DA
DA
(a)
(b)
310
260
H[m]
210
160
110
60
10
3
10
DA
(c)
Figura 7 - (a) Plano objetivo de las soluciones para el diseo de las chimeneas. (b) Frente ptimo de
Pareto para el diseo de las chimeneas. (c) Representacin de la variacin del impacto
ambiental mximo con la altura de la chimenea. Los puntos en negrita corresponden al
frente ptimo de Pareto.
34
Conclusiones
Se ha presentado una metodologa para el diseo de chimeneas que se basa en la resolucin del
problema inverso mediante la aplicacin de Algoritmos Genticos. El objetivo, para el problema
presentado, es el diseo de la chimenea de mnimo coste que cumpla con la legislacin
medioambiental aplicable. Para ello se utilizan modelos para el coste de la chimenea y el
impacto ambiental de los efluentes como funcin del dimensionado de la misma. Aunque el
modelo de dispersin de contaminantes es relativamente caro en trminos computacionales, la
paralelizacin del mtodo entre varios ordenadores permite obtener la solucin ptima en un
tiempo razonable.
Al analizarse un problema en el que existe una confrontacin de dos objetivos, en este caso el
impacto ambiental y el coste econmico, es posible la aplicacin del anlisis multiobjetivo. Este
mtodo genera gran cantidad de informacin, resultando una buena herramienta que ayuda a
tomar la decisin sobre las dimensiones ms apropiadas para la chimenea. Por esto, la
metodologa propuesta de diseo inverso de chimeneas resulta ms apropiada que la propuesta
por el Ministerio de Industria y Energa (1987), por ser capaz de generar ms informacin. Esto
no excluye tener que realizar un anlisis estadstico de la meteorologa de la zona objeto de
estudio.
La generalidad exigida al mtodo empleado para poder realizar el diseo inverso de chimeneas
industriales, hace que su uso se pueda extender a todo tipo de fuentes de contaminantes. As
pues, la metodologa de trabajo propuesta sera vlida tambin para emisiones lineales, como
carreteras; o de rea, como excavaciones al aire libre.
La metodologa empleada es muy flexible, y permite no slo el cambio sencillo de los
elementos del mtodo (por ejemplo, la funcin de coste), sino la adaptacin a problemas
enteramente distintos. Por ejemplo, con un modelo adecuado puede ser aplicado a la
optimizacin del proceso de fabricacin para reducir las emisiones de contaminantes. En la
referencia (Izquierdo 2003) se exploran algunas de estas aplicaciones alternativas.
35
Bibliografa
ESPAA. Decreto 833/75, que desarrolla la Ley 38/72 de proteccin del Ambiente Atmosfrico.
BOE nm. 96, 22 Abril 1975.
ESPAA. MINISTERIO DE INDUSTRIA Y ENERGA. Manual de Clculo de Altura de
Chimeneas Industriales. 1987.
36
Methods for Design, Optimization and Control with Applications to Industrial Problems.
EUROGEN 2003, CIMNE, Barcelona, Septiembre 2003.
NATALIE M. DSIDRI, J.A., LANTERI, S. Multi-Objective Optimization in CFD by Genetic
Algorithms. Rapport de recherch n 3686, INRIA. 1999.
OEI CK, GOLDBERG DE, CHANG SJ. Tournament Selection, Niching and the Preservation of
Diversity. Urbana, University of Illinois at Urbana-Champaign, Illinois Genetic
37
38