1/1 0 /2 0 2 4

CALIDAD DEL AIRE

Semana 07

Docente:
MSc. Ing. Jorge Alfonso Ugarte Díaz

Modelo Gaussiano de Dispersión de

Contaminantes Atmosféricos

1
 1/1 0 /2 0 2 4

PROPOSITO

Identifica r y comprender la
ecuación del modelo
gaussiano de dispersión de
contaminantes del aire

¿ Como definirías la dispersión de contaminantes

atmosféricos?

Actividad

• Indaga y compártenos un
concepto.

tiempo
2
5 min

plenaria
 1/1 0 /2 0 2 4

¿ Como definirías: modelo de dispersión de

contaminantes atmosféricos?

Actividad

• Indaga y compártenos un
concepto.

tiempo

5 min

plenaria

Modelos de dispersión

Son modelos
matemáticos que
permiten estimar la
concentración de un
contaminante en base
a un conjunto de
variables
3
 1/1 0 /2 0 2 4

Modelos de dispersión

¿Qué variables considerarías en un modelo de

dispersión de contaminantes atmosféricos?

Actividad Grupal

• Indaga y determina algunas de las variables que crees

que se debe considerar.

tiempo plenaria
4
5 min
 1/1 0 /2 0 2 4

Factores que determinan la dispersión de

contaminantes atmosféricos

Variables del orígen

• Flujo másico de emisión
• Altura de la emission (chimenea)
• Velocidad de salida del gas
• Temperatura de salida del gas
• Area de salida del gas (diametro
chimenea)

Variables meteorológicas
• Velocidad y dirección del viento
• Temperatura ambiental
• Presión barométrica
• Estabilidad atmosferica

Hipótesis del Modelo Gaussiano

Parámetros:
Q:Cantidad de contaminante
H: altura efectiva de chimenea (h +Δh )
u: Velocidad del viento a la altura efectiva
σz,σy: coeficientes de dispersión turbulenta

Z
C(x,y,z)

z X 5

y
x

Y
 1/1 0 /2 0 2 4

Ecuación General de Modelo Gaussiano

Ecuación Advección-Difusión Consideraciones

• Emisión constante
• No hay transformación química
𝜕𝐶 = −𝑢 𝜕𝐶 . −𝑦 𝜕𝐶 . −𝑧 𝜕𝐶 + 𝐾 𝜕 2 𝐶 + 𝐾 𝜕 2 𝐶 + 𝐾 𝜕 2𝐶 • Velocidad de viento constante
𝑥 𝑦 𝑧 • Terreno simple
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥2 𝜕𝑦 2 𝜕𝑧2 • Escala de tiempo: 10 – 20 min
• Reflexión turbulenta cerca
al suelo

Ecuación del modelo en (x,y,z) y en estado estacionario

−𝑦 2 − 𝑧−𝐻 2 − 𝑧+𝐻 2
𝑄 1 2𝛿 2 2𝛿 2 2𝛿 2
𝐶(𝑥, 𝑦, 𝑧) = 𝑒 𝑒 +𝑒
𝑢 2𝜋𝛿𝑦𝛿𝑧

−𝑦 2 − 𝑧−𝐻 2 − 𝑧+𝐻 2
𝑄 1 2𝛿 2 2𝛿 2 2𝛿 2
𝐶(𝑥, 𝑦, 𝑧) = 𝑒 𝑒 +𝑒
𝑢 2𝜋𝛿𝑦𝛿𝑧

Donde:
C :Concentración de contaminantes (g/m3)
Q:Cantidad de contaminante (g/s)

u: Velocidad del viento a la altura de la emisión (m/s)

6
δy, δz:Coeficientes de dispersión turbulenta (m)
H:Altura efectiva a la q es emitida el contaminante (m)
y:Distancia en el eje Y (m)
z:Distancia en el eje Z(m)
 1/1 0 /2 0 2 4

−𝑦 2 − 𝑧−𝐻 2 − 𝑧+𝐻 2
𝑄 1 2𝛿 2 2𝛿 2 2𝛿 2
𝐶(𝑥, 𝑦, 𝑧) = 𝑒 𝑒 +𝑒
𝑢 2𝜋𝛿𝑦𝛿𝑧

Donde:
C :C oncentración de contaminantes(g/m3)
Q:C antidad de contaminante (g/s)
u: Velocidad del viento a la altura de la emisión (m/s)
δy, δz: Coeficientes de dispersión turbulenta (m)
H: Altura efectiva a la q es emitida el contaminante (m)
y: Distanc ia en el eje Y (m)
z: Distanc ia en el eje Z(m) Z=0

Ecuación del área de concentración a nivel del

suelo (isolineas)

• Si solo son requeridos las concentraciones a nivel de

suelo (z=0)
• Genera curvas de isoconcentración en el plano xy

Donde:
2 C :g/m3
− −𝐻 2
Q:g/s
2 2 2𝛿 2
u:m/s
7
𝑦 𝑧 δy, δz:m
H: m
y:m
 1/1 0 /2 0 2 4

Inversión Térmica

Estabilidad Atmosférica

8
 1/1 0 /2 0 2 4

Inversión Térmica

Espejismos

Se trata de un fenómeno producido por la diferente temperatura a la cual se encuentran las

distintas capas de la atmósfera. Esta diferencia de temperatura genera variaciones en el índice de
refracción de las diversas capas atmosféricas, de modo que la luz no se propaga en línea recta,
sino que se va curvando gradualmente. El espejismo se produce, finalmente, porque el ojo sitúa la
imagen en la prolongación rectilínea de los rayos que inciden en él.

Inversión Térmica

9
 1/1 0 /2 0 2 4

Categorías de la atmosfera y efectos sobre la dispersión

Espiral Neutral

A D
Cono Fumigación

B E
Abanico Fumigación

C F

Clases de atmósfera según su estabilidad (Turner)

Viento de Día Noche*
superficie Radiación solar (W m-2)
(m s-1) Fuerte Moderada Débil Nubes** 4/8 Nubes**< 3/8
(a 10 m de altura) > 580 (1) 290 - 580 (2) <290 (3) (4) (5)
0-2 A A-B B E F
2-3 A-B B C E F
3-5 B B-C C D E
5-6 C C-D D D D
6 C D D D D

* La noche se define como el período desde una hora Categoría Estabilidad

después de la puesta de sol, hasta una hora antes de
A Extremadamente inestable
la salida del mismo
B Moderadamente inestable 10
** Nubosidad expresada en octavos de bóveda celeste C Ligeramente inestable
cubierta por nubes
D Neutra
Para cielos totalmente cubiertos, tanto para el día
E Ligeramente estable
como para la noche, debe asumirse clase de
estabilidad D F Moderadamente estable
 1/1 0 /2 0 2 4

1 2 3
Δh

Hh

6
4 5

H Δh

Coeficientes de dispersión turbulenta

δy , δ z

−𝑦 2 −𝐻 2
11
2𝛿 2 2𝛿 2
𝑦 𝑧
 1/1 0 /2 0 2 4

Coeficientes de dispersión - Ecuación de O. Martin

𝝈𝒚 = 𝒂. 𝒙𝒃 𝝈𝒛 = 𝒄. 𝒙𝒅 + 𝒇
- x (kilometros) y  (metros)!!!
- Las constantes a, c, d y f dependen de la categoría de estabilidad de Pasquill
- La constante b = 0.8941
x < 1 km x >1 km
a c d f c d f
A 213 440.8 1.941 9.27 459.7 2.094 -9.6
B 156 106.6 1.149 3.3 108.2 1.098 2
C 104 61 0.911 0 61 0.911 0 12

D 68 33.2 0.725 -1.7 44.5 0.516 -13

E 50.5 22.8 0.678 -1.3 55.4 0.305 -34
F 34 14.35 0.740 -0.35 62.6 0.180 -48.6
 1/1 0 /2 0 2 4

Actividad
Ejercicio
C alcule los coefic ientes de
dispersión turbulenta a una
distancia de 3 km de la
fuente de emisión en
dirección del viento
considerando un día nublado

Resolución
Datos necesarios:
Clase de Estabilidad: D
X=3 km
Empleando la ecuación de Martin

𝜎𝑦 = 68 (3)0.8941 = 181.59 𝑚
13
𝜎𝑧 = 44.5(3)0.516 + −13 = 65.44 𝑚
 1/1 0 /2 0 2 4

Altura efectiva (H) y del penacho (Δh)

−𝑦 2 −𝐻 2
2𝛿 2 2𝛿 2
𝑦 𝑧

Altura del penacho (Δh) - Ecuación de Holland

𝑣𝑠 . 𝑑 𝑇𝑠 − 𝑇𝑎
𝛥ℎ = 1.5 + 2.68x10−3. 𝑃. 𝑑
𝑢 𝑇𝑠
Δh
Donde:
H
vs: Velocidad de salida del gas (m/s)
d: Diámetro de salida de la emisión (m)
h
P: Presión barométrica (mbar)
Ts: Temperatura de salida del gas (K) 14

Ta: Temperatura del aire (K)

u: velocidad del viento a la altura efectiva (m/s)

 1/1 0 /2 0 2 4

Actividad

Ejercicio
Calcule la altura efectiva de la emisión:
Δh
Datos de la fuente:
Altura de la chimenea:80 m
Diámetro de la chimenea: 2.5 m H
Velocidad de salida del gas: 5.2 m/s h
Temperatura de salida del gas: 85°C
Met:
Temperatura ambiental:18°C
Presión barométrica :688 mbar
Velocidad del viento a la altura efectiva: 4.9 m/s

Resolución
Empleando la ecuación de
Holland con los datos del ejercicio

5.2(2.5) 358 − 291

𝛥ℎ = 1.5 + 2.68x10−3(688)(2.5)
4.9 358

𝛥ℎ = 6.27 m
finalmente
15

𝐻 = 80 m + 6.27 m = 86.27 m
 1/1 0 /2 0 2 4

Velocidad del viento a la altura efectiva(u)

−𝑦 2 −𝐻 2
2𝛿 2 2𝛿 2
𝑦 𝑧

Velocidad del viento

Medio Urbano
Terreno co mplejo

Medio Rural Terreno simple

16
 1/1 0 /2 0 2 4

Estación Meteorológica:
h
href = 10 m
uref
Ta
Pa
Rad

Velocidades del viento a la altura efectiva

𝐻
𝑢 = 𝑢𝑟𝑒𝑓
ℎ 𝑟𝑒
Los valores del coeficiente p como función de la clase de estabilidad y el entorno en que se
mueve el viento son los siguientes:

Categoría de Exponente para Exponente para

Estabilidad Medio rural Medio Urbano
A, B 0.07 0.15
17
C 0.10 0.20
D 0.15 0.25
E 0.35 0.40
F 0.55 0.60
30
 1/1 0 /2 0 2 4

Actividad
Calcule la velocidad del viento a la altura de
la emisión:
Datos de la fuente:
Ubicación: Medio Rural
Altura de la chimenea: 70 m
Diámetro de la chimenea: 2 m
Velocidad de salida del gas: 5.1 m/s
Temperatura de salida del gas: 98°C
Met:
Altura de la estación: 10 m
Temperatura ambiental: 17°C
Presión barométrica: 770 mbar
Velocidad del viento: 3.2 m/s
Radiación: 600 W/m2

Resolución
i.Clase de estabilidad atmosférica .
Velocidad del viento: 3.2 m/s
Radiación:600 W/m2
Clase de Estabilidad: B

ii. Exponente p
Clase de Estabilidad: B
Medio: rural
p=0.07
0.07
iii. Altura efectiva 70
𝑢 = 3.2
asumiendo de H =h 10 18
Por lo tanto H =70 m
𝑢 = 3.7 m/s
iii. Finalmente calculamos la velocidad
1/1 0 /2 0 2 4

19