Lady Julianalexandra Vásquez Delgado

MODELOS DE DISPERSIÓN

La dispersión es el transporte de contaminantes en el aire, la cual depende del estado de la atmósfera y

de las condiciones atmosféricas.

Modelos atmosféricos de Dispersión de Contaminantes

Los modelos atmosféricos de dispersión de contaminantes son herramientas físico-matemáticas que
permiten simular las condiciones reales de transporte y dispersión de los contaminantes como producto
de la interacción de las condiciones meteorológicas y las propias de las fuentes de emisión.

La información requerida considera la ubicación geográfica de las fuentes, la concentración de los

contaminantes emitidos, y factores meteorológicos tales como la turbulencia atmosférica, velocidad del
viento, temperatura y presión atmosférica. Por otro lado, los datos de salida de los modelos
(concentraciones diarias, horarias, etc) varían dependiendo de la complejidad de los modelos empleados.

Los modelos no describen con exactitud el problema de contaminación de una zona determinada. Estos
sólo constituyen una representación, cercana o no, de lo que sucede; por lo tanto, no hay que olvidar que
sólo son modelos y que la calidad de los resultados depende de diversos factores tal como:
 Exactitud del inventario de emisiones,
 Incertidumbre de los procesos físico químicos que sufren los contaminantes durante su transporte,
 Calidad y representatividad de los datos meteorológicos,
 Validación del modelo en una situación determinada. El modelo puede ser adecuado en una situación
particular pero no en otra.
Como se ha mencionado, los modelos de dispersión incluyen información tanto de la fuente emisora
como de las características particulares del lugar donde se ubica dicha fuente tal como se indica en la
Figura 1. Una vez ingresada toda la información, es posible determinar las concentraciones de los
contaminantes en un punto del espacio y en un tiempo determinado. De aquí que la concentración de un
contaminante dado, Ci, sea función del tiempo y del espacio:

Ci = Ci(x,y,z,t) donde Ci: es la concentración del contaminante i

1
 DATOS DE LA FUENTE
- Proceso: materias primas, productos,
subproductos, consumos, producciones.
- Chimenea: altura, diámetro.
- Gases: caudales, temperatura y
humedad.
- Contaminantes: caudales, medidas de
control.
NIVELES DE
CONCENTRACION
DATOS METEOROLOGICOS MODELO DE
- Máximos.
- Condiciones de estabilidad. DISPERSION - Promedio.
- Capas de inversión. - Probabilidad.
- Régimen de viento.
- Humedad relativa.
- Temperatura.
- Insolación.
- Precipitación.
- Nubosidad.
- Otros datos meteorológicos

DATOS TOPOGRAFICOS

Figura 1: Entradas y salidas de los modelos de dispersión.

Ecuación General de Difusión de Contaminantes

La concentración de un determinado tipo de contaminante en un volumen pequeño en el espacio y en un

tiempo dado se determina básicamente aplicando la ecuación de conservación de masa en un volumen
pequeño, que se considera como volumen de control.

C
 FMadvectiv o  FMdifusivo  Rxn Qca  Aporte fuentes  Re moción
t
El flujo másico advectivo (FMadvectivo) tiene relación con el movimiento de los contaminantes producto
del movimiento del fluido (en este caso de la atmósfera). Tiene dos términos en el balance de masa, uno
positivo (ingreso de flujo al volumen de control a través de una de sus superficies) y otro negativo (salida
de flujo a través de una de las superficies del volumen de control).

El flujo másico difusivo (FMdifusivo) es el paso de los contaminantes a través del volumen de control
debido a su propio movimiento. También tiene un término positivo y negativo en la ecuación de balance
de masa debido a la entrada y salida del contaminante por las superficies del volumen de control.

El término relacionado con reacción química (Rxn Qca) tiene que ver con la generación o eliminación de
contaminantes debido a procesos de reacción química que ocurren dentro del volumen de control.
Dependiendo del caso, el término debe ser considerado positivo o negativo en la ecuación de balance.

El aporte de fuentes corresponde a las emisiones de contaminantes que generan diversos tipos de
fuentes emisoras dentro del volumen de control considerado. Por último, el término correspondiente a la

2
remoción es negativo dada la eliminación de contaminantes producto de procesos de transformación u
otro.

Cada uno de estos términos representa fenómenos complejos donde se requiere mucha información
sobre los movimientos de la atmósfera, niveles de turbulencia, reacciones químicas, características de las
fuentes emisoras, etc. A pesar de la complejidad, es posible realizar una serie de simplificaciones
llegando a soluciones que quedan descritas por fórmulas relativamente simples.

Tipos de Modelos
La agencia de protección ambiental de Estados Unidos (EPA, por sus siglas en inglés) agrupa a los
modelos de calidad de aire (o de difusión) en las siguientes categorías:

a) Gaussianos:
Es la técnica más ampliamente usada para contaminantes no reactivos. Proporciona una excelente
aproximación matemática a la dispersión de contaminantes. Esta basado en la distribución Gaussiana o
distribución normal de estadísticas (se verá con mayor detalle más adelante en este capítulo).

El modelo Gaussiano asume que el material proveniente de una fuente, continuamente es transportado
en dirección del vector de velocidad de viento, estando las concentraciones más altas en el centro de la
pluma y las mas bajas en los extremos. En resumen, asume que la concentración de la pluma, en cada
distancia vientos abajo, posee una distribución Gaussiana o normal independiente horizontal y
verticalmente.

La mayoría de modelos recomendados por la Agencia de protección ambiental de los Estados Unidos son
Gaussianos.

b) Numéricos:
Este tipo de modelos es usado cuando los contaminantes analizados están reaccionando y formando
contaminantes secundarios. Son mucho más apropiados que los modelos Gaussianos para análisis de
una fuente de área urbana cuando están involucrados contaminantes reactivos; pero requieren unas
bases de datos mucho más extensas y complejas, por lo cual no son usados ampliamente.

Sus aplicaciones más comunes son para la formación de ozono y en el estudio de la conversión de SO 2 y
NOx para formar sulfatos y nitratos que contribuyen a la lluvia ácida. Un modelo numérico ampliamente
usado es el UAM (Urban Airshed Model) para la simulación de áreas urbanas calculando concentraciones
de ozono.

c) Estadísticos
Las técnicas estadísticas o empíricas son frecuentemente empleadas donde los procesos químicos y
físicos no poseen una clara interpretación científica, o en situaciones donde no se cuente con una fiable y
completa base de datos.

Están basados en técnicas estadísticas o semi-empíricas para analizar tendencias, relaciones de la

calidad del aire, las mediciones atmosféricas y para predecir la evolución de situaciones de
contaminación de corto plazo. Lo anterior se hace usando análisis de distribución de frecuencias, análisis
de series de tiempo y otros métodos estadísticos.

Poseen limitaciones por cuanto no se pueden establecer relaciones causa-efecto; sin embargo son muy
útiles para realizar predicciones de corto plazo en tiempo real en donde se cuenta con información de
concentraciones y sus tendencias.

d) De caja:
Este modelo asume que los contaminantes emitidos a la atmósfera se mezclan uniformemente en un
volumen o caja de aire de dimensiones finitas. Se deben establecer algunos parámetros como las

3
dimensiones con viento a favor, viento de costado y las dimensiones verticales de la caja. También se
debe establecer el periodo de tiempo de emisión de contaminantes.

Los modelos de caja suponen también que las emisiones se mezclan totalmente de una manera
inmediata con el aire disponible para la dilución. Se considera que los contaminantes son químicamente
estables y que permanecen en el aire (se verá con mayor detalle más adelante en este capítulo).

e) Físicos:
Involucra el uso de túneles de viento, canales de agua u otros medios para modelar fluidos. Esta clase de
modelo requiere un alto nivel técnico especializado, así como el acceso a los recursos necesarios.
Pueden ser muy útiles para simular situaciones de flujo complejo, con terrenos muy irregulares, donde se
presentan edificaciones o sobre alturas del terreno sobresalientes. Es aplicable a un grupo de fuentes
agrupadas en unos pocos kilómetros cuadrados.

Esta es la mejor técnica y la más aproximada a la realidad, sin embargo es poco utilizado por los recursos
que se deben poseer para utilizarla.

La Agencia Ambiental Europea (EEA) identifica, además de las anteriores, los siguientes tipos o técnicas
de modelamiento de la dispersión de contaminantes:

a) Modelos Eulerianos:
Es un modelo que hace una aproximación matemática de la dispersión de contaminantes a las
ecuaciones Eulerianas (conservación de la masa del contaminante). Estos modelos están usualmente
incluidos en modelos de pronóstico meteorológico. Los modelos de Eulerianos refinados incluyen
submodelos para simular turbulencias.

En estos modelos la dispersión de contaminantes se describe usando las ecuaciones básicas de la

hidrodinámica teniendo en cuenta:
· El balance del momentum
· El balance de calor
· El balance de humedad
· Las leyes de los gases
· Balance de masa del aire por debajo de la altura de mezcla
· El balance de los contaminantes estudiados

Estos modelos no funcionan bien en situaciones atmosféricas convectivas.

b) Modelos Lagrange:
Incluyen a todos los modelos en los que las plumas están “rotas” o segmentadas, emisiones puff o
partículas. Es decir cuando se trata de emisiones instantáneas no continuas.

Los modelos de Lagrange utilizan un cierto número de partículas ficticias para simular la dinámica de un
parámetro físico determinado. Estas partículas poseen las siguientes características:
· No tienen volumen
· Poseen una masa tal que la suma de todas las partículas emitidas en una unidad de tiempo es igual a la
tasa de emisión de la fuente
· No consideran flotación
· El promedio de las trayectorias de todas las partículas equivale a la pluma emitida por la fuente.

Este tipo de modelos se ajusta muy bien a la realidad en condiciones convectivas.

c) Modelos de receptor:
En contraste con modelos de dispersión (qué calculan la contribución de una fuente a un receptor como
el producto de la proporción de la emisión multiplicado por un coeficiente de dispersión), los modelos de
receptor parten de concentraciones observadas en un receptor y buscan repartir proporcionalmente las

4
concentraciones observadas en un punto de muestreo entre varios tipos de fuente. Están basados en
balances de masa.

En este curso se verán dos tipos de modelos, los modelos de caja y modelos gaussianos.

Modelo de Caja
Como se mencionó anteriormente este modelo asume que los contaminantes emitidos a la atmósfera se
mezclan uniformemente en un volumen o caja de aire de dimensiones finitas. La mezcla es prácticamente
instantánea. En general, las dimensiones en el sentido horizontal (eje x) de la caja son mucho mayores
que en el sentido vertical (eje z) tal como se muestra en la Figura 2. El tope vertical, en algunos casos,
está dado por el nivel de la capa de inversión.

Co, u u
C
z Q

Figura 2: Esquema del modelo de caja

Dado que el modelo asume que los contaminantes son químicamente estables y que permanecen en el
aire, los términos de reacción química y de remoción del balance de masa (ver ecuación de difusión) no
son considerados en este caso. Si además se considera que el movimiento advectivo de la atmósfera
predomina sobre el movimiento difusivo, éste último puede ser despreciado de la ecuación de difusión.
Por lo tanto, la ecuación de difusión considerará los siguientes términos:

Variación de masa = FlujoMásico advectivo que entra – FlujoMásico adevectivo que sale +
en la caja Aporte de Fuentes.

Si se considera que el origen de las coordenadas es uno de los vértices de la caja, las dimensiones de
ésta serán x, y, z. Los aportes de las fuentes, Q, son tal que se mezclan inmediatamente con el aire
circundante generando una concentración uniforme dentro de la caja, C. Por último, se considera que la
velocidad del viento es u y que antes de ingresar a la caja, el aire tiene una concentración del
contaminante estudiado Co (nivel base o background).

Por lo tanto, la ecuación de balance quedará:

d (Cxyz )
 uCo yz  Qxy  uCyz
dt
Si se alcanza el estado estacionario (no hay cambios en el tiempo de ninguna variable) se establece un
equilibrio entre los flujos de entrada y salida, donde la concentración de contaminantes en la caja y la que
sale de ella se puede estimar a partir de la siguiente expresión:

Qx
C  Co 
uz

5
Modelo Gaussiano
Como se mencionó anteriormente, los principios básicos de este modelo son los más utilizados para
describir el movimiento de los contaminantes en el atmófera. Supone que la contaminación se produce
“vientos abajo” de acuerdo a la dirección del viento (eje x), la que es contante y continua, y que se
produce a partir de una fuente puntual que emite en forma continua un flujo de contaminantes denotado
con la letra “Q”. Por otro lado, el modelo supone que el contaminante estudiado no es reactivo.

La distribución de contaminantes es tal que puede ser representada por una distribución gaussiana en el
eje z e y, donde la mayor concentración se produce en el eje del penacho tal como lo muestra la Figura
3. Hay que notar, que los coeficientes de dispersión en el sentido z e y no son los mismos, en ambas
direcciones existe una distribución gaussiana pero su dispersión cambia.

La concentración en un punto del espacio, C(x,y,z), suponiendo que no existe límite en la dispersión de
contaminantes, está dada por la siguiente expresión básica:
Q  y2 z2 
C ( x, y , z )  exp  2  2 
2u y z  2 2 z 
 y

Donde:
Q. es el flujo de emisión de la fuente contaminante.
u: es la velocidad promedio del viento
y: parámetro de dispersión en el sentido y
z: parámetro de dispersión en el sentido z

Los parámetros de dispersión dependen de la distancia de la fuente emisora (distancia x) y del nivel de
turbulencia de la atmósfera.

La expresión anterior es válida cuando la emisión se produce en un punto del espacio, localizado en el
origen del sistema de coordenandas (0,0,0). Pero dado que las emisiones que considera este modelo se
producen, por lo general, a través de una chimenea (ver Figura 3) la expresión para determinar la
concentración de una sustancia contaminante en un punto cualquiera del espacio, cambia de la siguiente
forma, considerando que la emisión se produce a una distancia (0,0,H ef):

Q  y 2     z  H ef  2    z  H ef  2 
C ( x, y , z )  exp   exp    exp  
2u y z  2 2    2 2   2 2 
 y   z   z 

En esta ecuación se supone que los gases contaminantes se “reflejan” en el suelo, lo cual se maneja
suponiendo que existe una chimenea virtual que tiene las mismas características que la chimenea real
pero dirigida hacia el plano negativo del eje z. La reflexión total de un gas en el suelo no existe, salvo en
determinadas superficies y en condiciones de temperatura y humedad específicas.

La altura efectiva de la chimenea está dada por la siguiente expresión:

Hef = Hchim + z

6
Donde Hchim corresponde a la altura geométrica de la chimenea y z es la elevación del penacho
producto de la temperatura caliente de los gases de salida con respecto a la temperatura de la atmósfera.
Esta elevación adicional depende de parámetros como la temperatura de los gases de salida, las
dimensiones de la chimenea, la temperatura de la atmósfera, la velocidad del viento, velocidad de salidad
de los gases a través de la chimenea, entre otros.

Figura 3: Esquema del modelo gaussiano

Los parámetros de dispersión se determinan una vez que se ha definido el nivel de turbulencia de la
atmósfera. Este nivel de turbulencia define la capacidad de la atmósfera para diluir más o menos la
contaminación que está siendo generada por una fuente puntual en un momento dado. Para definir el
nivel de turbulencia, lo que se hace en la práctica es trabajar con la “estabilidad” de la atmósfera a través
de un concepto cualitativo denominado “categorías de estabilidad”. El concepto más utilizado son las
categorías de estabilidad de Pasquill, que se muestran en la Tabla 1.

Una vez definida la categoría de estabilidad se determinan los parámetros de dispersión, los cuales
también dependen de la distancia a la fuente emisora. En resumen estos parámetros dependen de dos
aspectos básicos, que son los mismos en ambos casos pero la función es distinta:

x = f(distancia a la fuente emisora, categoría de estabilidad de la atmósfera)

z = g(distancia a la fuente emisora, categoría de estabilidad de la atmósfera)

7
Tabla 1: Categorías de Pasquill de estabilidad de la atmósfera.

En la Tabla 2 se muestran las expresiones de los parámetros de dispersión, tanto para zona rural como
para zona urbana.

Tabla 2: Expresión para los parámetros de dispersión, para 102< x <104 m