Tema 3

Calidad del aire

Mireia Udina

Calidad del aire

MIREIA UDINA

3. Modelos de dispersin de
contaminantes en la atmsfera
3.1 Los modelos de dispersin de contaminantes.................................................................. 1
3.1.1 Qu son............................................................................................................................ 1
3.1.2 Componentes de un modelo ......................................................................................... 2
3.1.3 Validacin de los modelos ............................................................................................ 3
3.1.4 Tipo y clasificacin de los modelos .............................................................................. 4
3.2 Estadsticos de dispersin de contaminantes ............................................................... 5
3.2.1 Plumas instantneas y plumas promedio ........................................................................ 5
3.2.2 Centro de la masa de contaminante y desviacin estndar ............................................ 5
3.3 Modelos de pronstico ..................................................................................................... 6
3.3.1 Modelo de caja o columna ............................................................................................ 6
3.3.2 Modelos de malla .......................................................................................................... 7
3.4 Modelos de emisiones ...................................................................................................... 7
3.4.1 Generalidades de las fuentes emisoras ........................................................................... 8
3.4.2 Clculos en un modelo de emisiones .............................................................................. 9
3.5 Modelo de diagnstico gaussiano ................................................................................. 11
3.5.1 Formulacin del modelo gaussiano ............................................................................... 12
3.5.2 Determinacin de los principales parmetros ................................................................ 15
a) Coeficientes de dispersin .................................................................................................. 15
b) Velocidad del viento ............................................................................................................ 16
c) Sobreelevacin del penacho ............................................................................................... 16
c.1) Para estratificacin inestable o neutral (PGT = A, B, C, D) ............................................. 18
c.2) Para estratificacin estable (PGT = E, F) ........................................................................ 18
c.3) Para estratificacin estable sin viento (calmas) (PGT = E, F) ......................................... 18
Problema 1 ejemplo .................................................................................................................. 19

3.1 Los modelos de dispersin de contaminantes

3.1.1 Qu son
Un modelo numrico es, en general, un conjunto de ecuaciones matemticas que relacionan
las variables atmosfricas que se quieren conocer (viento, temperatura, concentracin de un
contaminante, etc.) y los procesos fsicos y qumicos que producen los cambios de estas
variables. Ms concretamente, la finalidad de un modelo de dispersin atmosfrica es
diagnosticar o prever la concentracin C de un determinado contaminante i en un punto del
espacio x, y, z y en un determinado tiempo t, es decir: Ci (x,y,z,t). Un modelo matemtico de
dispersin atmosfrica debe tratar de simular el comportamiento en conjunto de las plumas
emitidas desde fuentes a nivel del terreno o desde la altura de la chimenea.

Los parmetros de entrada que necesitan los modelos son:

La cantidad de contaminante emitida por unidad de tiempo, la posicin y la altura de
emisin

1
 Tema 3
Calidad del aire
Mireia Udina

La velocidad y direccin de los vientos dominantes, la estabilidad atmosfrica y la

altura de mezcla
El comportamiento qumico del contaminante: posibles reacciones, vida media, etc.

Los modelos nos permiten explorar el estado de la atmsfera en reas para las que no
tenemos datos o que hay dificultades para obtenerlos. Los resultados se pueden obtener con
antelacin a que se presente el problema de contaminacin atmosfrica. Adems, podemos
obtener informacin de propiedades que no se pueden medir. Cualquier simulacin matemtica
de un fenmeno tan complejo como es la dispersin atmosfrica no es nunca exacta pero los
resultados de un modelo son el instrumento ms vlido en la planificacin de decisiones y en la
adopcin de medidas correctoras. Por eso, la modelizacin mejora nuestra comprensin del
funcionamiento del sistema ambiental estudiado.

La previsin y cuantificacin del impacto ambiental atmosfrico slo es posible cuando se ha

conseguido la modelizacin con suficiente representatividad de las caractersticas bsicas de
los medios emisor, difusor y receptor en su interrelacin espacial y temporal.

Algunas de las aplicaciones de los modelos son:

- Evaluaciones de impacto ambiental de uno o varios focos de contaminacin
atmosfrica de carcter puntual, lineal o superficial existentes o previstos
- Optimizacin de alturas de chimeneas para instalaciones industriales
- Estudios de contaminacin de fondo
- Planificacin urbana e industrial
- Diseo de redes de calidad del aire
- Predicciones de contaminacin del aire
- Programas de prevencin

3.1.2 Componentes de un modelo

La mayora de clasificaciones de modelos se basan en:
1. Las tcnicas utilizadas para construir el modelo
2. El grado de conocimiento e interpretacin del sistema que el modelo representa

Los componentes de un modelo aplicado a la atmsfera se pueden separar en:

Modelo conceptual: esquema bsico de los procesos y fenmenos que rigen el
sistema que se est estudiando, elaborado a partir de la recopilacin de informacin in
situ (experiencia de sistemas similares), que permiten tener una idea general de lo que
est pasando y de los componentes que participan en l.
Modelo matemtico: puede ser estadstico o determinista. Un modelo estadstico est
basado en el anlisis estadstico de series de datos obtenidas por sensores diversos a
partir de los cuales de infieren relaciones y tendencias. Un modelo determinista, en
cambio, simula los procesos atmosfricos relacionados con la meteorologa y la
contaminacin del aire y, generalmente, estn basados en las ecuaciones
fundamentales que describen los fenmenos que se tratan.
Modelo computacional: solucionan las ecuaciones del modelo determinista a travs de
tcnicas numricas cuando stas no se pueden resolver analticamente.

El proceso de modelizacin consistira en el desarrollo de los modelos anteriores como se

indica en la Figura 1.

2
 Tema 3
Calidad del aire
Mireia Udina

Figura 1. El proceso de la modelizacin numrica

3.1.3 Validacin de los modelos

La validacin y verificacin de un modelo son necesarias para asegurarse que los resultados
generados por el modelo son correctos. La validacin consiste en comprobar si el modelo
funciona de la manera esperada sobre un rango dado de condiciones. La verificacin trata de
comprobar si el resultado obtenido es comparable a las medidas directas de variables
utilizadas en el modelo.

Desafortunadamente, no es siempre posible validar un modelo, por ejemplo, cuando no es

posible medir todas la variables especificadas en el modelo acuradamente o cuando la escala
temporal o espacial de las medidas es muy grande. En estos casos ser necesario validar o
verificar el modelo contra otro modelo o mediante un anlisis de sensibilidad.

A menudo, se utilizan parmetros estadsticos para la verificacin y validacin. Algunos estn

resumidos en la Tabla 1.

Parmetro estadstico Definicin matemtica

Error medio (Mean 1 N
Bias) MB = (C m Co )
N 1
Error medio
1 N C m Co
normalizado (Mean MNBE = 100%
Normalized Bias Error) N 1 Co
Error medio absoluto 1 N
(Mean Absolute Gross MAGE = C m Co
Error) N 1
Error cuadrtico medio
1 N
(Root Mean Square RMSE = (Cm Co )2
Error) N 1
Tabla 1. Algunos parmetros estadsticos para la evaluacin de las previsiones de los modelos de
dispersin de contaminantes, dnde Cm es la concentracin dada por el modelo y Co la concentracin
medida con observaciones.

A continuacin se describe el significado de los estadsticos:

3
 Tema 3
Calidad del aire
Mireia Udina

Error medio (MB): el error medio indica, en promedio, si las previsiones estn por
encima o por debajo de las observaciones y en qu medida. La comparacin se realiza
para los mismos instantes de tiempo, espacio o tiempo y espacio.

Error medio normalizado (MNBE): representa el promedio de las diferencias entre los
resultados del modelo y las observaciones, normalizado segn las observaciones,
agrupndolas en tiempo o espacio. Un valor de 0 indicara que las sobreestimaciones e
subestimaciones del modelo frente a las observaciones se cancelan exactamente unas
con las otras.

Error medio absoluto (MAGE): representa una medida de lo cerca que estn en
promedio las medidas y las previsiones. La comparacin se realiza para los mismos
instantes de tiempo, espacio o tiempo y espacio.

Error cuadrtico medio (RMSE): mide la diferencia cuadrtica media entre las
observaciones y la previsin del modelo.

3.1.4 Tipo y clasificacin de los modelos

Podemos distinguir dos grandes tipos de modelos que simulan procesos fsicos similares: son
los modelos de diagnstico y los modelos de pronstico.

Los modelos de diagnstico proporcionan un diagnstico de las condiciones

presentes en un determinado instante a partir de una cierta cantidad de datos,
simulando las condiciones globales de la atmsfera. Los resultados de estos modelos
son campos de diferentes parmetros obtenidos a partir de interpolaciones y/o
extrapolaciones de medidas disponibles (meteorolgicas, de emisiones, etc.).

Los modelos de pronstico proporcionan una prediccin de la evolucin de la

atmsfera a partir de unas condiciones iniciales y de contorno. Nos dan una prediccin
para el futuro. Son modelos de mayor complejidad y en sus ecuaciones se incluyen
derivadas temporales.

Todos los modelos pretenden determinar la concentracin de un contaminante en un punto del

espacio y en un tiempo concreto o determinado. Segn el tratamiento de las ecuaciones de
transporte podemos distinguir entre modelos gaussianos, modelos eulerianos y modelos
lagrangianos. Los modelos gausianos son modelos de diagnstico. Los modelos eulerianos y
lagrangianos son modelos de pronstico.

Dependiendo de la forma de mirar el movimiento, podemos distinguir entre los modelos

eulerianos y los modelos lagrangianos.

Los modelos eulerianos usan un sistema de coordenadas fijo con respecto a la tierra.
Representan la clase ms sofisticada de modelos atmosfricos. En ellos las ecuaciones
que describen el movimiento y la transformacin qumica de los contaminantes en la
atmsfera se resuelven considerando un sistema fijo de coordenadas, y la regin a modelar
se puede dividir en cajas o celdas, tanto horizontal como verticalmente. Cuando el modelo
se simplifica y considera una sola caja, hablamos de los modelos de caja (o columna) y
cuando divide el espacio en mltiples celdas hablamos de los modelos de malla.

Los modelos lagrangianos usan un sistema de coordenadas que sigue el movimiento de

la atmsfera.

4
 Tema 3
Calidad del aire
Mireia Udina

3.2 Estadsticos de dispersin de contaminantes

3.2.1 Plumas instantneas y plumas promedio

La fotografa de plumas de contaminantes instantneas sera parecido a lo que vemos a simple

vista, con lmites poco definidos y movimientos arriba y abajo y de izquierda a derecha (Figura
2a). Cuando promediamos para un intervalo de tiempo la pluma emitida desde una chimenea
podemos aproximar que la mayora de contaminantes se encuentran cerca de la lnea central
(centerline) y que la concentracin decrece suavemente con la distancia a medida que nos
alejamos de la lnea central. El resultado del perfil de concentracin es una forma gaussiana (
Figura 2b).

Figura 2. (a) Fotografa de la pluma instantnea en distintos tiempos y su perfil de concentracin c en

altura, para la pluma ms oscura a la derecha. (b) Promedio de muchas plumas de contaminantes con el
perfil de concentracin promedio cavg a su derecha

3.2.2 Centro de la masa de contaminante y desviacin estndar

La lnea central de la pluma (zCL) se puede definir como la localizacin del centro de masa del
contaminante. Se define como el promedio pesado de la altura (zk) en cada altura k y la
concentracin correspondiente a cada altura (ck):
K

c z k k
k=1
zCL = K

c k
k=1
dnde K es el nmero total de alturas y k es el ndice de cada altura.

Para plumas fras que salen a poca velocidad de la chimenea la lnea central se puede
aproximar a la misma altura de la chimenea. Para plumas que salen a altas temperaturas, los
contaminantes sufren un empuje, una sobreelevacin (ver seccin 3.4) de manera que la lnea
central queda por encima de la altura de la chimenea.

Una vez promediada la pluma de contaminante, podemos calcular la desviacin estndar del
contaminante ( z ) que nos sirve para medir la dispersin de la pluma (Figura 3):

5
 Tema 3
Calidad del aire
Mireia Udina

1/2
K
ck (zk z )
2

z = k=1 K

ck
k=1
donde z = zCL , definido en la ecuacin anterior.
La funcin gaussiana en una dimensin se puede representar como en la Figura 3.

Figura 3. Curva gaussiana con unidades correspondientes a desviaciones estndar desde el centro. La
letra d corresponde a la amplitud nominal 1-D de la pluma, Q1 corresponde a la cantidad total de
contaminante emitido.

3.3 Modelos de pronstico

3.3.1 Modelo de caja o columna

Es un modelo euleriano de pronstico. Supone que la regin donde se quiere hacer la previsin
es una caja de dimensiones fijas y altura variable, en funcin de la altura de la capa de lmite (o
capa de mezcla) (Figura 4).

6
 Tema 3
Calidad del aire
Mireia Udina

Figura 4. Ilustracin de un modelo de caja

Se basa en el principio de conservacin de la masa expresado por la ecuacin:

c u 1 h Q uc 1 h
+ + c = i + b + c a + Ri (c1 , c 2 ,........ci .....c n )
t L h t h L h t
dnde:
- c es la concentracin media de las especies contaminantes
- u es la velocidad horizontal media
- h es la altura de la capa de mezcla
- L es la dimensin de la caja
- Qi es la emisin de la especie contaminante
- ca y cb son las concentraciones de fondo por encima de la capa de mezcla (para z>h) y
fuera del dominio respectivamente
- Ri es el trmino de reacciones qumicas que da cuenta del grado de produccin y/o
destruccin de las especies.

3.3.2 Modelos de malla

La concentracin de los contaminantes en cada celda se estima a intervalos especficos de

tiempo, teniendo en cuenta la informacin sobre campos meteorolgicos tridimensionales, as
como las concentraciones iniciales de los contaminantes, las emisiones, el transporte, la
dilucin y las transformaciones qumicas. La aplicacin de estos modelos resulta ms
conveniente cuando existen patrones complejos de emisin (por ejemplo, numerosas y
diversas fuentes de emisin, dispersas en un rea geogrfica amplia) y/o cuando los procesos
de transformacin qumica juegan un papel relevante en la generacin y destino de los
contaminantes (por ejemplo, para contaminantes secundarios). En general, cuando se cuenta
con informacin suficientemente detallada sobre el inventario de emisiones, la calidad del aire y
la meteorologa, estos modelos pueden aplicarse para evaluaciones detalladas de calidad del
aire a nivel urbano o regional.

3.4 Modelos de emisiones

7
 Tema 3
Calidad del aire
Mireia Udina

3.4.1 Generalidades de las fuentes emisoras

A igualdad de otros factores, la concentracin de contaminantes en un punto a sotavento del

foco es tanto mayor cuanto mayor sea la emisin o masa de contaminantes lanzadas en la
atmsfera. Los focos de emisin de contaminantes se llaman fuentes emisoras o simplemente
fuentes. Se llama intensidad de una fuente respecto a un contaminante a la masa de
contaminante que emite la fuente por unidad de tiempo, si bien la definicin debe matizarse
segn las caractersticas geomtricas de la fuente.

En un punto situado en el suelo a sotavento de la fuente, la concentracin de contaminantes es

tanto menor cuanto mayor sea la altura de la fuente. La relacin entre ambas, sin embargo, no
es proporcionalmente inversa y se ver en detalle en la seccin 3.5.

La concentracin de contaminantes en un punto dado no es solamente funcin de la emisin

sino tambin de la naturaleza y configuracin de las fuentes. A este respecto, las fuentes se
pueden clasificar en:

Puntuales. Las dimensiones de la fuente son despreciables frente a las del penacho
emitido. Por lo general, una chimenea aislada puede considerarse como fuente
puntual.
Lineales. Estn constituidas por una sucesin de focos emisores dispuestos a lo largo
de una lnea. Por ejemplo, una serie de chimeneas idnticas alineadas entre s y con el
viento es una primera aproximacin a una fuente lineal. Tambin lo es una autopista o
carretera concurrida.
La intensidad de esta clase de fuentes (puntuales y lineales) se define como la masa de
contaminantes emitida por unidad de longitud de la fuente y de tiempo.
Superficiales. Son fuentes de dimensiones no despreciables. Un ejemplo es un
vertedero, un estanque de depuracin de aguas residuales o un rea urbana o
industrial, en las que la densidad de fuentes individuales es suficientemente grande
para que no se puedan discernir los efectos propios de cada una sino, slo el efecto
conjunto.
En este caso la intensidad de la fuente superficial es la masa de contaminantes emitida por
unidad de tiempo y de superficie.
Volumtricas. Son emisiones producidas por fuentes tales como un edificio en cuyo
interior tiene lugar un proceso industrial que desprende gases a travs de chimeneas,
ventanas, puertas y aperturas de ventilacin. La intensidad es ahora la masa de
contaminantes emitida por unidad de tiempo y de volumen.
Difusas. Son fuentes de distribucin e intensidad irregulares y variables. El trfico
urbano es un ejemplo caracterstico. Su tratamiento es complicado y requiere tcnicas
especiales.

En cuanto al modo de emisin, una fuente puede ser continua (intensidad de emisin
constante) o discreta, en cuyo caso la fuente emite burbujas o bocanadas de contaminantes a
intervalos de tiempo regulares o no.
El concepto de fuente puntual continua se usar en muchos de los razonamientos y los
resultados obtenidos se aplicarn a otras clases de fuentes mediante artificios adecuados.

Se define la altura efectiva de la chimenea (He) como la altura sobre el terreno a la que el eje
del penacho se hace esencialmente horizontal.
A efectos del clculo de concentracin de contaminantes en un punto a sotavento de la fuente
interesa la altura efectiva adems de su altura geomtrica (H) (Figura 5). La diferencia entre
ambas es la llamada sobreelevacin (H) que se debe a que:

a) Los gases expulsados salen a una temperatura superior a la del aire ambiente, la cual
cosa les da un empuje aerosttico ascensional
b) En todas las chimeneas medianas o grandes los gases son emitidos a velocidades
-1
tpicas entre 10 y 20 ms

8
 Tema 3
Calidad del aire
Mireia Udina

La sobreelevacin no es, en general, despreciable frente a la altura geomtrica y puede llegar a

multiplicarla por factores entre 2 y 4.

Figura 5. Altura geomtrica (H) y altura efectiva (He = H+ h) de una chimenea

3.4.2 Clculos en un modelo de emisiones

Las emisiones son la entrada fundamental de un modelo acoplado de calidad del aire (un
modelo compuesto por varios modelos relacionados) y son totalmente necesarias para
representar bien la cantidad de contaminantes emitidos en la atmsfera. Las variaciones en las
emisiones suponen importantes variaciones en los campos escalares de concentracin de
varios contaminantes. Es necesaria una renovacin y mejora continua de los datos base del
modelo de emisiones para conseguir la mxima calidad posible en los resultados del modelo de
prediccin de calidad del aire.

Normalmente los modelos de emisiones suponen un gran esfuerzo inicial para conseguir y
tratar toda la informacin necesaria para el modelo. Una vez superado este proceso, el trabajo
es automatizarlo para extrapolar las emisiones a otras resoluciones o intervalos temporales.

Hay dos metodologas para calcular las emisiones en un modelo euleriano:

a) Metodologa Top & down: consiste en la adaptacin de las emisiones calculadas a
partir de otros modelos con menor resolucin espacial al rea que se quiere modelizar.
Habitualmente se recurre al inventario anual de emisiones EMEP/CORINAIR
(European Monitoring and Evaluation Programme /CORe INventory of AIR emissions)
(Figura 6) si no se dispone de una base de datos a nivel local.
Las fuentes emisoras presentes en este inventario se dividen en 11 sectores diferentes
segn las caractersticas y procedencias, teniendo en cuenta las principales fuentes
antropognicas, naturales y de trnsito martimo. Las emisiones recogidas en este
inventario se distribuyen por especies contaminantes y sectores de un dominio que
cubre toda Europa con una resolucin horizontal de 50 km. Para adaptar estas
emisiones a resoluciones mayores se utiliza una disgregacin espacial, que consiste en
dividir la emisin en las celdas que queden contenidas dentro de la celda mayor en
funcin del rea u otros criterios.

b) Metodologa Bottom-up: consiste en la adaptacin de las emisiones a partir de

inventarios pblicos que normalmente son ms detallados (mayor resolucin) que las
celdas del modelo de calidad del aire. Un ejemplo sera adaptar los datos de las
fuentes de emisin de las industrias de un polgono.

9
 Tema 3
Calidad del aire
Mireia Udina

Figura 6. Dominio de emisiones del EMEP CORINAIR

10
 Tema 3
Calidad del aire
Mireia Udina

3.5 Modelo de diagnstico gaussiano

Una fuente que emite continuamente en la atmsfera produce un penacho de contaminantes.

Las partculas o gases inyectados a la atmsfera en un mismo punto siguen trayectorias
ligeramente diferentes, separndose entre ellas a medida que se alejan de la fuente emisora.

Suponiendo que la tasa de emisin y las condiciones atmosfricas son constantes, se llega a
una condicin estacionaria, en la cual el penacho adquiere una forma constante en el tiempo.
La concentracin es mxima en el eje del penacho, disminuyendo hacia los bordes (Figura 7)
de tal forma que la distribucin de la concentracin se puede aproximar a una funcin
gaussiana (ver seccin 3.2).

Figura 7. Coordenadas en un penacho gaussiano

El modelo gaussiano es el ms usado en la prctica para estimar la dispersin de

contaminantes. El eje es el de la direccin del viento. El eje representa la direccin vertical.
El eje es la direccin ortogonal al plano . Si representamos la concentracin del
contaminantes vemos que se ajusta a la curva de Gauss en los dos planos y .

e
z = Direccin vertical
x = Distancia en direccin del viento
y = Direccin ortogonal al plano xz
h = Sobreelevacin del penacho
H = Altura de chimenea
He = Altura efectiva = h + h

11
 Tema 3
Calidad del aire
Mireia Udina

3.5.1 Formulacin del modelo gaussiano

La funcin de Gauss es la que ms se ajusta en los dos planos, de manera que la

concentracin en cualquier punto (C) sigue la frmula:

( 1)
C = A exp(By 2 )

donde y es la distancia al eje y A y B son constantes a determinar. A representa la

concentracin en el eje, cuando y = 0 , es decir C0

C = C0 exp(By 2 )

Para determinar el coeficiente B utilizamos propiedades de la funcin gaussiana. Por definicin

en la funcin gaussiana cuando B 0 , C 0 , de manera que la dispersin lateral del
penacho es infinita.
Segn una propiedad de la funcin gaussiana, el rea limitada por la curva, el eje de abscisas y
las ordenadas trazadas a distancias y = 2 y del eje (2 veces la desviacin estndar),
comprende el 95.45% del rea total bajo la curva (es decir, en nuestro caso el porcentaje de
masa total del contaminante). Entonces si tomamos el coeficiente C/C0 y para y = 2 y las
tablas de la funcin dan:

C / C0 = 0.1353 = exp(B4 2y )

Podemos aislar B y resolver entonces:

1
B=
2 y2

Entonces la funcin resultante es:

" y2 %
C = C0 exp $$ 2 ''
# 2 y &
En el plano haramos lo mismo, de manera que:

" z2 %
C = C0 exp $ 2 '
# 2 z &
Entonces, para un punto del espacio x,y,z e instante de tiempo t, la concentracin sera:

" y2 % " z2 %
C(x, y, z, t) = C0 exp $$ 2 '' exp $ 2 '
# 2 y & # 2 z &

Podemos estimar la concentracin C0 en funcin de la emisin (Q), la velocidad del viento

medio (u) y una constante (k), ya que sabemos que la concentracin es inversamente
proporcional a la velocidad del viento:

Q
C0 = k
u

12
 Tema 3
Calidad del aire
Mireia Udina

Entonces:

Q " y2 % " z2 %
C(x, y, z, t) = k exp $$ 2 '' exp $ 2 ' ( 2)
u # 2 y & # 2 z &

Adems, sabiendo que toda la masa del contaminante emitido tiene que atravesar un plano
indefinido contenido en un volumen elemental (dx, dy, dz), se deduce que:

+ +
Q= Cudydz ( 3)

Sustituyendo C por la expresin en (2):

+ +
Q $ y2 ' $ z2 '
Q= k u u exp && 2 2)
) exp & 2 ) dydz ( 4)
% y ( % 2 z (

La emisin Q que se supone constante y la velocidad del viento se eliminan de manera que la
expresin queda:

1 + + " y2 % " z2 %
= exp $$ 2 '' exp $ 2 ' dy dz
k # 2 y & # 2 z &
Y utilizamos de nuevo una propiedad de la gaussiana en que

+

exp (By ) dy =
2

Donde nuestra B:
1
B=
2 y2
Entonces

1 1
= ; =
k 2 y2 2 z2 k 1
2 y z
y nos queda:
1
k=
2 y z

Y finalmente logramos la ecuacin bsica de concentracin gaussiana:

" y2 % " (z H )2 % ( 5)
Q e
C(x, y, z, t) = exp $$ 2 '' exp $ 2
'
2 y z u # 2 y & # 2 z &

13
 Tema 3
Calidad del aire
Mireia Udina

Vemos que directamente la concentracin no depende de la distancia a la fuente (x) pero

(5
veremos que en el clculo de los parmetros que intervienen en ) s que aparece la variable
x.

Teniendo en cuenta que z es la altura que medimos sobre el suelo y por tanto hay que restarle
la altura efectiva (He), e H = H + h , donde H es la altura de la chimenea y h la
sobreelevacin del penacho.

Con la expresin ( 5) podemos calcular la concentracin del contaminante en cualquier punto

conociendo la emisin (Q), el viento medio a la altura efectiva u(H e ) y las sigmas y , z .

Hay, sin embargo, una contribucin adicional a C no incluida en ( 5) y debida al efecto del
suelo. Cuando a una cierta distancia de la chimenea el penacho alcanza el suelo, la mayor
parte de los contaminantes son reflejados por l y vuelven a la atmsfera, dando una
concentracin superior a la calculada con ( 5). Se introduce el efecto de la pluma imagen, una
chimenea ficticia, simtrica a la real que emite un penacho imagen especular del real respecto
al suelo. La concentracin de contaminantes correspondiente a la zona sombreada en la Figura
8 es la suma de la obtenida con la expresin ( 5) ms la debida a la chimenea imagen que se
simula sumando a la segunda exponencial:

" y 2 %( " (z H )2 % " (z + H )2 %+ ( 6)

Q e e
C(x, y, z, t) = exp $$ 2 ''*exp $ 2
' + exp $ 2
'-
2 y z u # 2 y &*
) # 2 z & # 2 z &-,

z=h

z=0

Figura 8. Penacho real y penacho imagen

Para un sensor al nivel del suelo, z=0, y situado a lo largo de la proyeccin del eje del penacho
sobre el suelo, y=0, tenemos:

Q ( " H 2 %+ ( 7)
C(x, 0, 0) = *exp $ e2 '-
y z u *) # 2 z &-,

Vemos que en las expresiones ( 6) y ( 7) intervienen los siguientes datos:

a) Datos de la fuente: se trata de todos los datos de la instalacin de la fuente como las
coordenadas del punto, la intensidad de emisin, la altura de la chimenea, la
temperatura y velocidad de salida de los gases, etc.
b) Datos meteorolgicos: seran la direccin y velocidad del viento, las variables ligadas a
la intensidad de turbulencia, etc. Pueden ser datos ficticios o bien resultado de una
medida o datos climatolgicos.

14
 Tema 3
Calidad del aire
Mireia Udina

c) Parmetros: son valores conocidos empricamente, para los cuales existen mltiples
opciones en la bibliografa. Los ms importantes son los coeficientes de dispersin (las
sigmas), la sobreelevacin del penacho y la altura efectiva.

La eleccin de las frmulas de las sigmas y la sobreelevacin es una cuestin delicada. La

velocidad del viento se necesita a la altura del penacho (en rigor, se requerira la media integral
de dicha velocidad desde el borde inferior al borde superior del penacho), mientras que el dato
disponible suele ser el viento a 10 m y, en algunos casos en torres de 30 o 50 m.

Examinando la ecuacin ( 7), se aprecia que en el denominador del primer trmino figura el
producto de las sigmas que aumentan con la distancia. El numerador Q, en cambio, se
mantiene constante; el trmino, pues, disminuye al aumentar la distancia. Por otra parte la
exponencial aumenta con la distancia, tendiendo a la unidad. Por lo tanto, el producto, es decir
la concentracin C(x, 0, 0) crece al principio con la distancia, pasa por un mximo y disminuye
despus, tendiendo a 0 cuando x tiende a infinito.

Mientras que las fuentes a nivel del suelo producen concentraciones tanto ms altas cuanto
ms flojo es el viento, con fuentes elevadas no es obvio a priori cules son las condiciones que
conducirn a la mxima concentracin en tierra. El viento favorece la dispersin pero reduce la
sobreelevacin y ambos efectos se contraponen. Entonces, hay que recurrir a un mtodo de
clculo por aproximaciones sucesivas.

3.5.2 Determinacin de los principales parmetros

Los parmetros que necesitamos calcular para resolver la ecuacin de concentracin de
contaminantes son: los coeficientes de dispersin o sigmas ( y , z ), la velocidad del viento a

la altura efectiva ( u(H e ) ) y la sobreelevacin (( h )) del penacho (para calcular la altura

efectiva ( H e )).

a) Coeficientes de dispersin
Para obtener los coeficientes de dispersin o sigmas ( y , z ) lo ideal sera disponer de
medidas de la intensidad de turbulencia. Entonces, los coeficientes se calcularan a partir de
las fluctuaciones de la direccin de la componente horizontal y del mdulo de la componente
vertical del viento. No obstante, en la prctica no disponemos de estas medidas y se recurre a
otros procedimientos. El ms frecuente consiste en parametrizar la intensidad de la turbulencia
en trminos de variables meteorolgicas corrientes como la direccin y velocidad del viento, la
temperatura y su gradiente vertical y la irradiancia. A menudo, se clasifica la turbulencia en
categoras, como por ejemplo, usando las Categoras de Dispersin de Pasquill-Gifford-Turner
(PGT) (ver tema 2, seccin 2.3.3, CA_lectura_2.3.pdf). Identificada la categora de dispersin,
los coeficientes y , z se pueden obtener mediante frmulas, tablas o grficas que en general
se obtuvieron empricamente a partir de experimentos u observaciones. Unas de las ms
conocidas son las frmulas de Briggs que son curvas de ajuste de y , z en funcin de la
distancia a la fuente y que son varan segn las categoras PGT y de si el medio es rural o
urbano. La frmula principal es:

y,z = mx(1+ nx) p ( 8)

donde x es la distancia a la fuente (en m) y los coeficientes y exponentes m, n, p varan

segn la clase de dispersin (Tabla 2).

15
 Tema 3
Calidad del aire
Mireia Udina

Tabla 2. Coeficientes m, n y exponentes p en la expresin para distancias entre 100 m y 10 km.

b) Velocidad del viento

Como vimos en el tema 2, seccin 2.3.3 (B), y explicado en la lectura 2.3, el viento se puede
aproximar a una variacin logartmica con la altura. En este caso, para calcular el viento a la
altura efectiva utilizaremos el ajuste potencial emprico del viento que nos permite calcular el
viento a una cierta altura sabiendo el dato del viento a otra altura de referencia y con el
coeficiente n que depende de la categora PGT y de si el medio es rural o urbano (ver tabla 6.1
de la CA_lectura_2.3.pdf).

Ejemplo: queremos saber el viento a 100 m y sabemos que una estacin de medida registra un
viento a 10 m de 5 m/s en un medio urbano y una categora de PGT clasificada como A.

Solucin:
Viendo la tabla 6.1 de CA_lectura_2.3.pdf, n=0.15. Calculamos el viento a 100 m
n 0.15
! 100 $ ! 100 $
u(100m) = u(10m) # & = 5* # & = 7.06m / s
" 10 % " 10 %

c) Sobreelevacin del penacho

La observacin de los efluentes que emergen de una chimenea demuestra que el penacho en
general sigue ascendiendo durante un cierto tiempo antes de propagarse horizontalmente, es
la llamada sobreelevacin (h), que sera la diferencia entre la altura efectiva He y la altura
geomtrica H (Figura 5), de manera que: H e = H + h
Esta sobreelevacin depende de:

16
 Tema 3
Calidad del aire
Mireia Udina

Pujanza o fuerza ascensional: cuando se trata de un penacho caliente hay una

diferencia de temperatura de los gases (Tg) y la temperatura del aire (Ta)
Condiciones meteorolgicas: si la situacin es de estabilidad o inestabilidad
Caractersticas de la fuente: radio de la fuente y velocidad de salida de los gases.

Se definen los siguientes parmetros para penachos calientes:

P es el parmetro de pujanza
"T T %
P = gVg r 2 $ g a '
# Ta &
2
donde g es la constante de gravedad g = 9.8m / s , Vg es la velocidad de salida de los
gases, r es el radio de la fuente, Tg es la temperatura de salida de los gases y Ta la
temperatura del aire a la altura geomtrica.
Para calcular P la Ta(H) se tiene que calcular a la altura geomtrica. Teniendo la
temperatura en superficie, utilizamos el clculo de la variacin de la temperatura con la
altura aplicando el gradiente trmico de la atmsfera.
Ta (H ) = Ta H
Para estimar el gradiente trmico miramos la Tabla 3 y sabiendo la categora PGT podemos
estimar el gradiente . Para una estratificacin indiferente normalmente supondremos que el
gradiente es adiabtico, es decir 0.98 K/100 m.

Tabla 3. Categoras de estratificacin de PGT segn la velocidad del viento y el gradiente vertical de
temperatura. Fuente: Safety Series (1980): Atmospheric dispersin in nuclear in nuclear power plant
siting, N.

x f es el valor de x a partir del cual la sobreelevacin no crece, tiene valores distintos:

x f = 3.5x * para casos inestables o neutrales
*
x representa el valor de x por el cual la turbulencia atmosfrica empieza a destruir el
penacho y adquiere los siguientes valores
x * = 14P 5/8 para P < 55
x * = 34P 2/5 para P 55
3.14u
xf = para casos estables
S1/2
El parmetro de estabilidad de Hesselberg S es:

17
 Tema 3
Calidad del aire
Mireia Udina

g
S= ( )
Ta
Con estos datos, podemos utilizar las frmulas de Briggs para el clculo de la sobreelevacin
en diferentes casos.

c.1) Para estratificacin inestable o neutral (PGT = A, B, C, D)

Se usa una de las expresiones:

1.6P1/3 x 2/3
h = para x < x f o x < 10H ( 9)
u
para distancias inferiores a 10 veces la altura geomtrica

2/3
1.6P1/3 (10H )
h = para x x f o x 10H ( 10)
u
para distancias superiores o iguales a 10 veces la altura geomtrica

c.2) Para estratificacin estable (PGT = E, F)

1.6P1/3 x 2/3
h = para x < x f o x < 10H ( 11)
u
para distancias inferiores a 10 veces la altura geomtrica

1/3
"P%
h = 2.4 $ ' para x x f o x 10H ( 12)
# uS &
para distancias superiores o iguales a 10 veces la altura geomtrica

c.3) Para estratificacin estable sin viento (calmas) (PGT = E, F)

5P1/4
h =
S 3/8
Hemos aplicado los casos para penachos calientes
Penacho caliente: cuando Tg - Ta > 25 C.
Penacho fro: cuando Tg - Ta < 25 C.

La sobreelevacin para penachos fros viene dada por la expresin:

2/3
" %
$ '
V / u ' " x %1/3
h = 1.195$ g $ '
$ 3 ' #r&
1+
$ V /u'
# g &

El clculo de la sobreelevacin a menudo consiste en un proceso iterativo de clculo de la

altura efectiva, velocidad del viento a la altura efectiva, luego clculo de la sobreelevacin, y
volver a empezar: altura efectiva velocidad del viento a la altura efectiva, luego clculo de la
sobreelevacin, etc.

18
 Tema 3
Calidad del aire
Mireia Udina

Problema 1 ejemplo
8
Supongamos una fuente situada en un medio rural que emite 1.67x10 mg/s de SO2 por una
chimenea de 250 m de altura. En una noche de cielo nuboso la velocidad del viento de un
anemmetro situado a 10 m de altura mide 3.5 m/s. Un experto nos suministra los datos de los
coeficientes de dispersin a 2.5 km de distancia horizontal del pie de la chimenea, y = 250m
y z = 165m .
Se desea calcular la concentracin de SO2 en el suelo a 2.5 km de distancia horizontal a la
fuente,
a) sin suponer sobreelevacin y admitiendo el dato de viento ledo
b) suponiendo sobreelevacin y calculando el viento a la altura efectiva. La chimenea
tiene 5 m de dimetro, los gases salen a una velocidad de 15 m/s y a una temperatura
de 150 C. La temperatura del aire es de 17 C. Supondremos que hay una
estratificacin neutral.

Suponer tambin que el clculo se realiza suponiendo el punto bajo el eje del penacho.

Solucin

a) Sin sobreelevacin admitiendo el dato de viento ledo

Entonces la altura efectiva coincidira con la altura geomtrica. Aplicamos la ecuacin ( 7)
porque nos dicen que estamos en el suelo (z=0) y en el punto bajo el eje del penacho (y=0) de
manera que:

1.67 108 ) # 250 2 &, 3

C(x, 0, 0) = +exp % 2(
. = 116.89mg m
250 165 3.5 * $ 2(165) '-
La concentracin obtenida de SO2 sin tener en cuenta la sobreelevacin ni la variacin del
-3
viento con la altura, en el suelo, por debajo del eje del penacho es de 116.89 mgm .

b) Con sobreelevacin y calculando el viento a la altura efectiva

1. Calculamos el parmetro de empuje

Primero, la temperatura a la altura geomtrica, sabiendo que la estabilidad es neutral, es decir
el gradiente es el adiabtico.

Ta (H ) = Ta H = (17 + 273) (0.98 250 /100) = 287.55K = 14.33 !C

"T T % " 423 287.55 %
P = gVg r 2 $ g a ' = 9.815 2.52 $ 4 3
' = 432.78m s
# Ta & # 287.55 &

2. Calculamos el viento a la altura efectiva

Como la altura efectiva es lo que queremos calcular, hacemos un mtodo iterativo de

aproximaciones sucesivas utilizando el ajuste potencial del viento (el exponente n=0.15 para
medio rural y categora D).

n 0.15
u !z$ ! 250 $ 1
= # & u = 3.5# & = 5.67ms
u0 " z0 % " 10 %

19
 Tema 3
Calidad del aire
Mireia Udina

Seleccionamos la expresin para la sobreelevacin: caso estratificacin neutral para el caso en

que x 10H , ecuacin (9):
2/3
1.6P1/3 (10H ) 1.6 432.781/3 (10 250)2/3
h1 = = = 393.18m
u 5.67
Volvemos a calcular el viento a esta altura
0.15
! 643.18 $ 1
u = 3.5# & = 6.54ms
" 10 %
1.6 432.781/3 (10 250)2/3
h2 = = 340.87m
6.54
Y as sucesivamente hasta obtener valores que converjan

Podemos crear una tabla para los datos que vamos obteniendo en las iteraciones

Iteracin 0 1 2 3 4
-1
U (m s ) 3.5 6.54 6.45 6.46 6.46
H (m) 250 250 250 250 250
H (m) 0 393.2 340.9 345.6 345.1
He (m) 250 643.2 590.9 595.6 595.1

Vemos que de la cuarta a la quinta iteracin los valores convergen, y ya tenemos entonces la
altura efectiva He=595 m

Ahora ya podemos calcular la concentracin

Q ( " H 2 %+ 1.67 108 ( " 595.12 %+ 3 3

C(x, 0, 0) = *exp $ e2 '- = *exp $ '- = 0.299mg m = 299 g m
y z u *) # 2 z &-, 250 165 6.46 ) # 2(165)2 &,

Vemos que el valor de la concentracin con sobreelevacin es mucho menor (del orden de
microgramos) comparado con el valor sin sobreelevacin (del orden de miligramos), por lo
tanto, hay que considerar la sobreelevacin para un clculo fiable del campo de concentracin.

20