Dispersión de Contaminantes en La Atmosfera
Dispersión de Contaminantes en La Atmosfera
Dispersión de Contaminantes en La Atmosfera
w w
𝜕<𝑐>
< 𝑐 𝑉 ™ > = −𝐾™ , 𝑖 = 𝑥, 𝑦, 𝑧
𝜕𝑖
Que mas interesa es el dispersión en el dirección vertical 𝐾•
Debido la naturaleza aleatoria de las fluctuaciones de la dirección horizontal del
viento, la distribución media de las concentraciones a lo largo plazo del eje (y) se
podrá describir mediante una distribución gaussiana , pero ahora es claro que el
ancho de dicha distribución (o su desviación estándar) va a depender de cuan lejos
este la fuente y de las condiciones locales de estabilidad de la atmósfera.
Sistema de coordenadas para describir el penacho de una
chimenea usando la ecuación gaussiana.
• Se toma perspectiva pseudoestacionaria, en cual se supone que el tiempo de
viaje del penacho hasta llegar a la posición x viene dado por x/<V>, donde <V>
es la velocidad media del viento medida a la altura de descarga de chimenea.
• Además, suponiendo que en la dirección x se puede despreciar la dispersión
advectiva, se deduce que en condiciones turbulentas el penacho sigue siendo
gaussiano, y la concentración queda descrita por la expresión:
V V V
1 𝑦 (𝑧 − ℎ• )
< 𝑐 𝑥, 𝑦, 𝑧 >= 𝐴𝑒𝑥𝑝 − V +
2 𝜎œ 𝜎•V
ℎ• - altura efectiva de chimenea o de estabilización del penacho, determinada
por el balance de la boyancia y momentum de la emisión.
𝑄
𝐴=
2𝜋 < 𝑉 > 𝜎œ 𝜎•
𝜎œ 𝜎• desviaciones estándar varían tanto con la distancia a la fuente como según
las condiciones meteorológicas
Q tasa de emisión de la fuente µg/s
A distancias tales que 2,15*𝜎 es mayor o igual a la altura efectiva del penacho (punto 1 en la fig.), hay que
considerar el efecto de superficie del suelo (z=0). En tal caso penacho comienza a tocar el suelo y hay que
conocer como interactúa el contaminante con el suelo (condición de borde de problema). Si el
contaminante no se absorbe en la superficie; la ecuación gaussiana queda como:
𝑄 −𝑦 V (𝑧 − ℎ• )V (𝑧 + ℎ• )V
< 𝑐 𝑥, 𝑦, 𝑧 >= 𝑒𝑥𝑝 V 𝑒𝑥𝑝 − V + 𝑒𝑥𝑝 −
2𝜋 < 𝑉 > 𝜎œ 𝜎• 2𝜎œ 2𝜎• 2𝜎•V
Ÿ cœ § (¨© )§
< 𝑐 𝑥, 𝑦, 0 >= ¡¢£¤¥¤¦
𝑒𝑥𝑝 V¤§ 𝑒𝑥𝑝 − V¤§ - para concentraciones superficiales (z=0)
¥ ¦
Ÿ (¨© )§
< 𝑐 𝑥, 0,0 >= 𝑒𝑥𝑝 − - para concentraciones superficiales en eje de emission (y=0)
¡¢£¤¥¤¦ V¤¦§
Ÿ
𝑠𝑖 ℎ• = 0: < 𝑐 𝑥, 0,0 >= - si fuente de emission es superficial
¡¢£¤¥¤¦
Variación de velocidad del
viento con elevación
• En la altura h vamos a tener viento con velocidad:
¬
ℎ
< 𝑉 > = 𝑉•
ℎ•
Clase de estabilidad n
atmosférica
A,B,C 0,2
D 0,25
E 0,33
F 0,5
Altura efectiva de chimenea ∆ℎ = ℎ• − ℎQ
𝑣® 𝑑 𝑄¨•,¯
• Para condiciones superadiabáPcas ∆ℎ = 3,47
𝑉
+ 5,15
𝑉
(inestables)
𝑣® 𝑑 𝑄¨•,¯
∆ℎ = 0,35 + 2,64
𝑉 𝑉
• Para la estabilidad neutral
𝑣® 𝑑 𝑄¨•,¯
∆ℎ = −1,04 + 2,24
𝑉 𝑉
• Para condiciones subadiabáPcas
(estabes)
𝑣® = velocidad de salida del gas (en m/s),
d = diámetro de la chimenea (en m), y
Qh = tasa de emisión de calor de la chimenea (en kJ/s).
Sensabilidad de formula
h-altura de chemenea
Ws-velocidad de emission
Ta-temperatura de ambiente
Ts-temperatura de emission
• Una planta de energía tiene una chimenea con un diámetro de 2 m y
emite gases con una velocidad de salida de la chimenea de 15 m/s y
una tasa de emisión de calor de 4.800 kJ/s. La velocidad del viento es
de 5 m/s. La estabilidad es neutral. Estime el aumento de la pluma. Si
la chimenea tiene una altura geométrica de 40 m, ¿cuál es la altura
efectiva de la chimenea?
𝑣® 𝑑 𝑄¨•,¯ 15 4 2 4800
∆ℎ = 0,35 + 2,64 = 0,35 + 2,64 = 38,7𝑚
𝑉 𝑉 5 5
ℎ• = 38,7 + 40 = 78,7m
• Se emite dióxido de azufre a una tasa de 160 g/s desde una chimenea
con altura efectiva de 60 m. La velocidad del viento a la altura de
emisión es de 6 m/s, y la clase de estabilidad atmosférica es D,
suponiendo un día nublado. Calcule la concentración máxima a nivel
del suelo a una distancia de 500 metros desde la chimenea.
1. Vamos a calcular dispersión horizontal para
condiciones rurales y clase de estabilidad atmosférica
D:
𝑥 c•,¯
𝜎œ = 0,08 4 𝑥 1 + a
10
500 c•,¯
= 0,08 4 500 1 + = 39,04𝑚
10000
Dispersión vertical:
1,5𝑥 c•,¯
𝜎• = 0,06 4 𝑥 1 + U = 22,7 𝑚
10
Considerando que el SO2 no se absorbe en el suelo
(estimación conservadora) podemos usar la ecuación:
𝑄 (ℎ• )V
< 𝑐 𝑥, 0,0 >= 𝑒𝑥𝑝 −
𝜋 < 𝑉 > 𝜎œ 𝜎• 2𝜎•V
V
𝑄 (ℎ• )V −𝑦 V 50
< 𝑐 𝑥, 𝑦, 0 >= 𝑒𝑥𝑝 − 𝑒𝑥𝑝 = 292,3 exp −0,5 = 128,7µ𝑔/𝑚U
𝜋 < 𝑉 > 𝜎œ 𝜎• 2𝜎•V 2𝜎œV 39,04
292,3 µ𝑔/𝑚U
𝜇𝑔/𝑠 𝑚V 𝑚V
V V
= 𝜇𝑔/𝑚U
𝑚/𝑠 4 𝑚 4 𝑚 𝑚 𝑚
• Incorporación de la cinéPca de primer
orden en el modelo gaussiano
La concentración obtenida por el modelo gaussiano en sus disPntas modificaciones debe corregirse en el caso de que el contamínante en estudio sufra una reacción
química o fotoquímica que provoque que su concentración disminuya con el Pempo. Para ello basta recurrir a las ecuaciones cinéPcas correspondientes e integrar la
ecuación respecto del Pempo. Por ejemplo, si la molécula se degrada según una ecuación de primer orden, la concentración respecto del Pempo es:
𝑑𝐶
= 𝐶• exp(−𝑘𝑡)
𝑑𝑡
donde k es la constante de velocidad. Si definimos el Pempo de vida medio t1/2 de la molécula como el Pempo necesario para que C(t)=C0/2, tenemos:
t1/2 =ln2/k - Pempo de vida media
El factor de 1/2 en el primer término surge porque las par“culas que caen no se reflejan en la superficie del
suelo.
• La tasa, ω, a la que se deposita la materia parPculada en el suelo, está
relacionada a la concentración ambiental como se muestra en
𝜔 = 𝑣® < 𝑐 𝑥, 𝑦, 0 >
donde ω = la tasa de sedimentación (en g/m²/s).
º ¼
» §
Ÿ cœ § (¨© c½¾¿)
𝜔 = 𝑣® < 𝑐 𝑥, 𝑦, 0 > = 𝑣® 𝑒𝑥𝑝 𝑒𝑥𝑝 −
V ¡¢£¤¥¤¦ V¤¥§ V¤¦§
• La planta generadora de electricidad de combustión emite 1,1 kg/min
de SO2 de una chimenea con una altura efectiva de 60 m. la emisión
consiste partículas de 10 μm de diámetro y con una densidad de 1
g/cm3, calcule la concentración del nivel del suelo a 500 m a lo largo
de la línea central de la pluma,y la tasa de sedimentación en ese
punto. La viscosidad del aire es 0.0185 g/m/s a 25 ° C.
• Estabilidad atmosférica D. 𝜎œ = 35m, 𝜎• = 19𝑚 velocidad del viento
5 m/s
𝜌 1𝑔/𝑐𝑚U
𝑣® = 𝑔𝑑V V c¯
= 0.0029𝑚/𝑠 = 9.8 𝑚/𝑠 (10 𝑚) V
𝑔
18𝜇 10c° 𝑚U /𝑐𝑚U 18 4 0.0185 𝑚 /𝑠
Q = 1.1 kg/min = 18 g/s.
𝑣® 𝑥 V
𝑄 −𝑦 V (ℎ• − < 𝑉 > )
< 𝑐 500,0,0 >= 𝑒𝑥𝑝 𝑒𝑥𝑝 −
2𝜋 < 𝑉 > 𝜎œ 𝜎• 2𝜎œV 2𝜎•V
0.0029𝑚/𝑠500𝑚 V
18𝑔/𝑠 (60𝑚 − ) 𝑔
5𝑚/𝑠
= 𝑒𝑥𝑝 − = 19,4 4 10c° = 19,4µ𝑔
2𝜋5𝑚/𝑠 4 35𝑚 4 19𝑚 2 4 (19 𝑚)V 𝑚U
x P x
Tasa de sedimentación: 𝜔 = 𝑣® < 𝑐 500,0,0 > = 19,4 4 10c° PÀ 0.0029 ®
= 5,6 4 10cd P§®
• Cuando tenemos el caso de una emisión lineal, la dispersión de
contaminantes también se puede modular en enfoque gaussiano.
Para esto se considera que la emisión proviene desde una fuente
lineal coincidente con eje Y. En tal caso que X es la dirección
perpendicular a la línea de emisión o de la carretera, la concentración
que sopla en forma perpendicular a la carretera viene dada por:
ÂV
2𝑞 (ℎ• )V exp(−0,5𝑝^2)
< 𝑐 𝑥, 0,0 >= 𝑒𝑥𝑝 − ” 𝑑𝑝
2𝜋 < 𝑉Á > 𝜎• 2𝜎•V ÂL 2𝜋
𝐿1 = 𝑦L /𝜎œ ; 𝐿2 = 𝑦V /𝜎œ extremos (estandartizados) de la carretera o línea de emisión
Integral del lado derecho corrige por efectos de tamaño de la línea de emisión. Si se aproxima
que la línea de emisión es muy larga, la integral de L1 a L2 tiende 1,0
q – intensidad de emisión por unidad de longitud (g/(s*m))
La dispersión horizontal no aparece en el primer factor del lado derecho, puesto que no puede
haber variabilidades a lo largo de la coordenada Y. Por la misma razón esta coordenada no
aparece en factor exponencial.
< 𝑉Á > - componente del viento perpendicular a la carretera
• Calculemos la concentración total de hidrocarburos en un punto
ubicado a 100 m viento debajo de una carretera. La velocidad del
viento perpendicular a la carretera es 4m/s; flujo del trafico es de
8000 vehículos por hora y la velocidad media de estos es 50 km/h. La
emisión media de cada vehículo en dichas condiciones es de 0,02 g/s.
Supongamos que se trata de una carretera muy larga
ÂV
2𝑞 (ℎ• )V exp(−0,5𝑝^2)
< 𝑐 𝑥, 0,0 >= 𝑒𝑥𝑝 − ” 𝑑𝑝
2𝜋 < 𝑉Á > 𝜎• 2𝜎•V ÂL 2𝜋
º©Ä
d••• Ä K•¨ K•¨ x x
Calculamos numero de autos por unidad de longitud: Å =0,16 P
, factor de emisión: 𝑞 = 0,16 P
0,02 ® K•¨
= 0,0032 ®P
¯•••• Ä
𝑔
2 4 0,0032 𝑠 𝑚 𝑔
< 𝑐 𝑥, 0,0 >= = 91,2 4 10c° U = 91,2𝜇𝑔/𝑚U
2𝜋 4 4𝑚/𝑠 4 7𝑚 𝑚
Modelos para Dispersión de Fuentes Lineales