SIMULACIÓN DE LA DISPERSION DE CONTAMINANTE MEDIANTE UN

MODELAMIENTO DEL PROCESO DE DISPERSION DE CONTAMINANTES

DEL AIRE

Quispe Quispe Gian Marco Adrian

Modelamiento y simulación ambiental, Facultad de Ingeniería Ambiental y Recursos Naturales, Universidad Nacional Del Callao,
Lima, Perú 2019

Resumen: El presente trabajo consiste en la simulación de la dispersión de contaminantes en el aire, el cual

hallaremos a través de dos ecuaciones, la primera que está dada por la ecuación general de dispersión de
contaminantes en el aire (ecuación de Gauss) el cual estima las concentraciones a diferentes distancias, y
la segunda que está dada por la ecuación de M. Embaby et al.(Cmax) el cual nos ayuda a obtener la
concentración máxima y el punto de localización de dicha concentración. Por lo tanto, se estableció una
situación en el cual tomaremos como contaminante disperso en el aire al dióxido de azufre (SO 2) para poder
hacer las comparaciones con dichas ecuaciones mencionadas se considerara los siguientes datos: un
caudal de SO2: Q g/s=26.7, una altura física de la chimenea de 59(m), el diámetro de la chimenea de 1.98
m, la velocidad de salida del gas Vgas(a=4, b=7)(m/s), temperatura de salida del gas T°gas(205°C),
presión(1.25 atm), Vviento(v1=1.5 y v5=5)(m/s), Temperatura del ambiente(21°C), con una dirección del
viento predominante(oeste a este). Se le añadirá la Concentración de fondo:26.3 ug/m 3, una constante
cinética de reacción química en la atmosfera 1.64*10 -4s-1.
Se calculó la concentración para 3 poblaciones de las cuales 3 sobrepasan el ECA para SO 2= 80 ug/m3 en
las estabilidades B y D, también se halló las distancias que separan a cada una de las poblaciones al punto
de máxima concentración. Finalmente se halló la altura mínima de la chimenea para que las 3 poblaciones
cumplan con el ECA.

Palabras claves: Ecuación de Gauss, Ecuación de M Embaby, Concentración Máxima, Dióxido de Azufre,
Calidad del Aire, Simulación, Altura de la Chimenea, Estabilidad.

Abstract: The present work consists of the simulation of the dispersion of pollutants in the air, which we will
find through two equations, the first one given by the general equation of dispersion of pollutants in the air
(Gaussian equation) which It estimates the concentrations at different distances, and the second is given by
the equation of M. Embaby et al. (Cmax) which helps us to obtain the maximum concentration and the
location point of said concentration. Therefore, a situation was established in which we will take sulfur
dioxide (SO2) as a dispersed pollutant in the air to be able to make comparisons with these mentioned
equations, the following data will be considered: a flow of SO2: Q g / s = 26.7, a physical chimney height of
59 (m), the chimney diameter of 1.98 m, the gas outlet speed Vgas (a = 4, b = 7) (m / s), gas outlet
temperature T ° gas (205 ° C), pressure (1.25 atm), Wind (v1 = 1.5 and v5 = 5) (m / s), Ambient temperature
(21 ° C), with a predominant wind direction (west to east ). The background concentration will be added: 26.3
ug / m3, a kinetic constant of chemical reaction in the atmosphere 1.64 * 10-4s-1.
The concentration was calculated for 3 populations of which 3 exceed the RCT for SO2 = 80 ug / m3 in the
stabilities B and D, the distances that separate each of the populations to the point of maximum
concentration were also found. Finally, the minimum height of the chimney was found so that the 3
populations comply with the ECA.

Keywords: Gaussian equation, M Embaby equation, Maximum concentration, Sulfur dioxide, Air quality,
Simulation, Chimney height, Stability.

1. INTRODUCCIÓN ecológico general, cuyas consecuencias a largo

plazo no son fácilmente previsibles.
La explotación intensiva de los recursos naturales Los modelos de concentración de contaminantes
y el desarrollo de grandes concentraciones son protocolos matemáticos que proporcionan
industriales y urbanas en determinadas zonas, una estimación de la concentración de un
son fenómenos que, por incontrolados, han dado contaminante.
lugar a la saturación de la capacidad asimiladora En el presente trabajo haremos una comparación
y regeneradora de la Naturaleza y pueden llevar a utilizando la ecuación general de Gauss para
perturbaciones irreversibles del equilibrio dispersión de contaminantes en el aire con la
 ecuación de M. Embaby et al. Utilizando el SO 2 Como el contaminante del aire se acumula justo
como contaminante en el aire. Asimismo, en la dirección de la velocidad del viento, es decir,
utilizaremos una concentración de fondo y como a lo largo del eje x, donde y = 0, entonces el
actúa con reacción o sin reacción química. punto en el eje x donde se da la Cmax, será:

I. ECUACIÓN GENERAL DE GAUSS

Es una de las expresiones matemáticas muy

utilizadas en la actualidad para predecir cómo
varía la concentración de un contaminante
atmosférico en la atmósfera a una distancia Y el máximo nivel de concentración Xmax se
determinada, partiendo de un foco fijo de emisión. obtendrá con la siguiente ecuación:
(Birds, 1995)

Q
C ( x , y , z , H )= .¿
π .u . σy .σz
Donde:
En la siguiente tabla se muestran los siguientes
● C: concentración del contaminante µg/m3
contantes a diferentes estabilidades:
● x: Distancia a la chimenea en dirección del
viento(m)
● y: Distancia en dirección horizontal, del foco Estabilidad\ a b c d
fijo constantes
Muy 0.00022 2.1 0.59 0.85
● z: distancia en dirección vertical, del foco fijo inestable (A)
de emisión Inestable (B) 0.056 1.1 0.41 0.86
● Q: flujo másico del contaminante µg/seg. Ligeramente 0.12 0.1 0.24 0.88
● u : Velocidad del viento m/sec inestable (C)
● δy: Desviación estándar el contaminante en la Neutro (D) 0.73 0.5 0.14 0.89
dirección y, 5
Ligeramente 0.82 0.4 0.11 0.89
● δz: Desviación estándar del contaminante en
estable (E) 8
la dirección z.
Estable (F) 0.63 0.4 0.075 0.89
● H: Al tura eficaz de la chimenea. 5
Tabla 1 Valores de las constantes a, b, c y d en diferentes
Para el presente trabajo usaremos la ecuación de clases de estabilidades
gauss de primer orden:
III. DIOXIDO DE AZUFRE (SO2).
Q
C= exp ¿
π . u . σy . σz Se emite a la atmósfera en forma de SO 2 durante
la quema de combustibles y el procesamiento de
A medida que la pluma viaja los coeficientes de los minerales. Es un gas incoloro, irritante, no
dispersión del contaminante están en función a la inflamable y con un olor penetrante que consiste
distancia de ocurrencia en un momento dado, en un átomo de azufre y dos de oxigeno u vida
estos coeficientes es posibles estimarlos media en la atmósfera es corta, de unos 2 a 4
gráficamente y matemáticamente a través de días.

los valores de dispersiónσy y σz se calcularon El SO2 también se emplea en la industria del

matemáticamente según Martin, a partir de las papel y la pulpa de madera
ecuaciones como agente blanqueador.

2. OBJETIVOS
σy=a . x b σz=c . x d
● Comprobación del cálculo del modelo
Donde “x” es la distancia recorrida del gas desde
desarrollado.
la salida de la chimenea y “a, b, c, d” son
● Analizar el proceso de dispersión de
constantes que dependen de la estabilidad dada
contaminantes en el aire
que ya se encuentran tabuladas.
3. MATERIALES Y METODOS
II. ECUACIÓN DE EMBABY
3.1. Materiales
Para el presente trabajo se utilizó una laptop
Ahora, consideramos el caso donde la altura
personal con el programa de hojas de cálculo
efectiva H es constante, es decir, no depende
Excel. Y se emplearon las fórmulas de
de x e y.
Gauss y Embaby (“Máximum Ground Level
 Concentration of Air Pollutant”) para 1. El comportamiento de las
determinar el comportamiento del SO2 en el concentraciones de SO2 en las
aire. estabilidades B, C y D.

3.2. Métodos a) Sin considerar Cfondo ni reacción química

Para poder aplicar las fórmulas de dispersión en la atmosfera
de contaminantes en el aire y analizar el
proceso de dispersión de contaminantes
Velocidad del gas 4m/s
estableceremos las siguientes para el
cálculo:
a) Fuente: QgSO2/s=26.7; hf=59(m);
diámetro de la chimenea (1.98 m);
Vgas(a=4, b=7)(m/s); Tgas(205°C);
presión(1.25 atm).
b) Condiciones atmosféricas: Vviento(v1=1.5
y v5=5)(m/s); Tambiente(21°C); dirección
del viento predominante(oeste a este).
c) ECA del SO2= 80 ug/m3
d) Cfondo=(26.3 ug/m3)
e) Cte cinética de reacción química en la
atmosfera (1.64*10-4)s-1.
f) Localización; hay tres poblaciones (P-1,P-
2 Y P-3) en dirección al este. Grafico n°1: estimaciones de las
P-1: x1=1270m e y1=82m. concentraciones de SO2 a diferentes
P-2: x2=1435m e y2=77m. estabilidades (B, C Y D) Y diferentes
P-3: x3=1751m e y3=80m. velocidades de viento.
e) aplicar las Ecuaciones a 3 estabilidades:
“B, C, D”.
Velocidad del gas 7m/s
DETERMINAR

1. Estimar las concentraciones de SO2 en las

estabilidades B, C y D cada 150m hasta los
1,5 Km y cada 500 m desde los 1,5Km hasta
los 12Km considerando que sigue la dirección
predominante del viento, para:
a) Sin considerar Cfondo ni reacción química
en la atmosfera.
b) Considerando Cfondo ni reacción química
en la atmosfera.
2. Estimar la Cmax del SO2 aplicando la formula
general y comparar con el obtenido por la
ecuación M. Embaby et. al. (2002). Para:
a) Sin considerar Cfondo ni reacción química Grafico n°2: estimaciones de las
en la atmosfera concentraciones de SO2 a diferentes
b) Considerando Cfondo ni reacción química estabilidades (B, C Y D) Y diferentes
en la atmosfera
velocidades de viento.
3. Estimar la concentración en c/u de las 3
poblaciones a las 3 estabilidades; así mismo
calcular la distancia en (m) que los separa de
la CMAX considerando solo viento de 2m/s, el b) Considerando Cfondo ni reacción química en la
valor de la concentración de fondo y reacción atmosfera
química atmosférica. Velocidad del gas 4m/s
4. Indica en cual población o poblaciones
sobrepasara el ECA del SO2, de sobrepasar
el ECA estimar cual sería la altura que la
chimenea mínima garantice no sobrepasar el
ECA.

4. RESULTADOS
 Cuadro n°2: comparación de la concentración
Grafico n°3: estimaciones de las concentraciones máxima Gauss-Embaby a diferentes
de SO2 a diferentes estabilidades (B, C Y D), estabilidades.
diferentes velocidades de viento y concentración
de fondo.

Concentración máxima (ug/m3), con Vgas=7m/s sin

Velocidad del gas 7m/s reacción química y Viento=1.5m/s.

Grafico n°4: estimaciones de las concentraciones Cuadro n°3: comparación de la concentración

de SO2 a diferentes estabilidades (B, C Y D), máxima Gauss-Embaby a diferentes
diferentes velocidades de viento y concentración estabilidades.
de fondo.

2. Estimación de la concentración máxima y

comparación de la ecuación general con la Concentración máxima (ug/m3), con vgas=7m/s sin
ecuación de Embaby. reacción química y Viento=5m/s.
a) Sin considerar Cfondo ni reacción química en la
atmosfera
Concentración máxima (ug/m3), con Vgas=4m/s sin
reacción química y Viento=1.5m/s.

Cuadro n°4: comparación de la concentración

máxima Gauss-Embaby a diferentes
estabilidades.
Cuadro n°1: comparación de la concentración
máxima Gauss-Embaby a diferentes
estabilidades.
Concentración máxima (ug/m3), con Vgas=4m/s sin
reacción química y Viento=5m/s.
b) Considerando Cfondo ni reacción química en la
atmosfera

Concentración máxima + Cfondo (ug/m3), con

Vgas=4m/s sin reacción química y Viento=1.5m/s.
Cuadro n°5: comparación de la concentración Cuadro n°8: comparación de la concentración
máxima Gauss-Embaby a diferentes máxima Gauss-Embaby a diferentes estabilidades
estabilidades.
3. Estimación de la concentración de SO 2
considerando reacción química,
Concentración máxima + Cfondo (ug/m3), con concentración de fondo y viento de 2m/s en 3
Vgas=4m/s, sin reacción química y Viento=5m/s. poblaciones y la distancia que los separa de la
Cmax.
● Concentración de SO2 en las poblaciones
(ug/m3) y Viento=2m/s a diferentes estabilidades.

Población n°1 concentración (ug/m3) con

Vgas=4m/s y Vviento=2m/s.
Estabilida concentració
c+cf con Rx
d n
B 126.16 152.46 137.38
C 0.00 26.30 23.69
D 183.77 210.07 189.29
Cuadro n°6: comparación de la concentración Cuadro n°9: estimación de concentración
máxima Gauss-Embaby a diferentes considerando reacción química (Rx) y la
estabilidades. concentración de fondo (c+cf).

Población n°1 concentración (ug/m3) con

3
Concentración máxima + Cfondo (ug/m ), con Vgas=7m/s y Vviento=2m/s.
Vgas=7m/s, sin reacción química y Viento=1.5m/s Estabilida concentració
c+cf con Rx
d n
B 124.38 150.68 135.77
C 0.00 26.30 23.69
D 147.86 174.16 156.93
Cuadro n°10: estimación de concentración
considerando reacción química (Rx) y la
concentración de fondo (c+cf).

Población n°2 concentración (ug/m3) con

Vgas=4m/s y Vviento=2m/s.
Estabilida concentració
Cuadro n°7: comparación de la concentración c+cf con Rx
máxima Gauss-Embaby a diferentes estabilidades d n
B 103.45 129.75 115.34
C 0.00 26.30 23.38
D 206.46 232.76 206.91
Cuadro n°11: estimación de concentración
Concentración máxima + Cfondo (ug/m3), con considerando reacción química (Rx) y la
Vgas=7m/s sin reacción química y Viento=5m/s. concentración de fondo (c+cf).
 D 1.32E+03 1647.55 0
Población n°2 concentración (ug/m3) con
Cuadro n°16: estimación de concentración
Vgas=7m/s y Vviento=2m/s.
máxima solo con Vviento=2m/s.
Estabilida concentració
c+cf con Rx
d n
B 102.33 128.63 114.35 ● Distancia entre el punto de concentración
máxima con las tres poblaciones.
C 0.00 26.30 23.38
D 170.73 197.03 175.15 Sabiendo que los puntos de ubicaciones de las 3
Cuadro n°12: estimación de concentración poblaciones son:
considerando reacción química (Rx) y la X(m) Y(m)
concentración de fondo (c+cf). P-1 1270 82
P-2 1435 77
Población n°3 concentración (ug/m3) con P-3 1751 80
Vgas=4m/s y Vviento=2m/s.
Estabilida concentració Població Distancia desde el punto de
c+cf con Rx
d n n concentración máxima (m)
B 73.52 99.82 86.47
P-1 804.45
C 0.00 26.30 22.78
P-2 2.39E+22
D 227.88 254.18 220.17
Cuadro n°13: estimación de concentración P-3 302.06
considerando reacción química (Rx) y la Cuadro n°17: cálculo de la distancia entre el
concentración de fondo (c+cf). punto de concentración máxima con las 3
poblaciones, considerando una altura efectiva
Población n°3 concentración (ug/m3) con de 64.94 m y altura física de la chimenea de 59
Vgas=7m/s y Vviento=2m/s. m.
Estabilida concentració
c+cf con Rx
d n Població Distancia desde el punto de
B 73.01 99.31 86.02 n concentración máxima (m)
C 0.00 26.30 22.78 P-1 775.30
D 195.62 221.92 192.23 P-2 4.65E+22
Cuadro n°14: estimación de concentración P-3 130.77
considerando reacción química (Rx) y la Cuadro n°17: cálculo de la distancia entre el
concentración de fondo (c+cf). punto de concentración máxima con las 3
poblaciones, considerando una altura efectiva
de 69.493 m y altura física de la chimenea de 59
● Concentración máxima (C-max) (ug/m3) y m.
ubicación del punto de concentración máxima
en (m) 4. Estimación de la altura de la chimenea para
que no sobrepase el ECA. Considerando
Concentración máxima a diferentes estabilidades estabilidades (B, C, D) y viento de 2m/s en 3
Con Vgas=4m/s en (ug/m3) Xmax Ymax poblaciones
ESTABILIDAD C-MAX (m) (m)

B 9.62E+02 469.74 0 Población 1

Altura
2.3921E+2 ECA mínima de la
C 8.54E-13 0 Estabilidad C con Rx
2 80(ug/m3) chimenea
(m)
D 1.57E+03 1459.73 0
Cuadro n°15: estimación de concentración B 137.38 sobrepasa 147.21
máxima solo con Vviento=2m/s.
C 23.69 no sobrepasa *
Concentración máxima a diferentes estabilidades D 189.29 sobrepasa 74.78
Con Vgas=7m/s en (ug/m3) Xmax Ymax Cuadro n°18: cálculo de la altura mínima de la
ESTABILIDAD C-MAX (m) (m) chimenea para que no sobrepase el ECA para
B 8.55E+02 499.05 0 SO2, considerando una altura efectiva de 64.94
4.6548E+2 m y altura física de la chimenea de 59 m inicial.
C 4.45E-13 0
2
 SO2, considerando una altura efectiva de 64.94 m
y altura física de la chimenea de 59 m inicial.

Población 1
Altura
mínima de
ECA
Estabilidad C con Rx la Población 3
80(ug/m3) chimenea
Altura mínima
(m) Estabilida ECA
C con Rx de la
B 135.77 sobrepasa 142.74 d 80(ug/m3) chimenea (m)
no
C 23.69 * B 86.02 sobrepasa #¡NUM!
sobrepasa
D 156.93 sobrepasa 70.31 no
C 22.78  *
sobrepasa
Cuadro n°19: cálculo de la altura mínima de la
chimenea para que no sobrepase el ECA para D 192.23 sobrepasa 81.19
SO2, considerando una altura efectiva de 69.493 Cuadro n°23: cálculo de la altura mínima de la
m y altura física de la chimenea de 59 m inicial. chimenea para que no sobrepase el ECA para
SO2, considerando una altura efectiva de 69.493
m y altura física de la chimenea de 59 m inicial.
Población 2
Estabilida ECA
Altura mínima 5. CONCLUSIONES
C con Rx de la
d 80(ug/m3) chimenea (m) ● Podemos concluir que para la estabilidad “C”
B 115.34 sobrepasa 131.5 no se pudieron obtener concentraciones y solo
no se obtuvo más que la concentración de fondo
C 23.38 * 26.7 ug/m3.
sobrepasa
● Las poblaciones con mayor concentración de
D 206.91 sobrepasa 79.91 SO2 en el aire fue la población 3 en la
Cuadro n°20: cálculo de la altura mínima de la estabilidad “D” con 220.17 ug/m 3 el cual se
chimenea para que no sobrepase el ECA para consideró una altura efectiva de 64.94 m y
SO2, considerando una altura efectiva de 64.94 m altura física de la chimenea de 59 m inicial. E
y altura física de la chimenea de 59 m inicial. igual excedió el ECA para SO2.
● La población más alejada del punto de
concentración máxima fue la población 2, el
Población 2 cual se consideró una altura efectiva de
Estabilida ECA Altura mínima 69.493 m y altura física de la chimenea de 59
C con Rx de la
m.
d 80(ug/m3) chimenea (m)
● En las 3 poblaciones se excedió el valor del
B 114.35 sobrepasa 127.08 ECA para SO2, en las estabilidades B y D.
no ● No se pudo determinar la altura mínima de la
C 23.38 * chimenea 3 para la estabilidad “B”.
sobrepasa
D 175.15 sobrepasa 75.43
6. REFERENCIAS BIBLIOGRAFICAS
Cuadro n°21: cálculo de la altura mínima de la
chimenea para que no sobrepase el ECA para ● Daniel Pereira-Peláez (2017) Simulación de la
SO2, considerando una altura efectiva de 69.493 dispersión de contaminantes en la atmósfera de
m y altura física de la chimenea de 59 m inicial. una planta de generación de electricidad a
biomasa. Empresa GREENCON S.R.L.
Consultora en Desarrollo y Medio Ambiente
Población 3 Calle Andrés Muñoz N° 2405, Sopocachi,
Estabilida ECA Altura mínima ciudad de La Paz, Bolivia.
C con Rx de la
d 80(ug/m3) chimenea (m) ● M. EMBABY, A. B. MAYHOUB, K. S. M. ESSA
B 86.47 sobrepasa #¡NUM! and S. ETMAN.(2002) Maximum ground level
no concentration of air pollutant.
C 22.78 * 
sobrepasa
● Jorge A. Cabrera Fernández (2012)
D 220.17 sobrepasa 85.67 APLICACIÓN DE UN MODELO DE
Cuadro n°22: cálculo de la altura mínima de la DISPERSIÓN ATMOSFÉRICA. PONTIFICIA
chimenea para que no sobrepase el ECA para
 UNIVERSIDAD CATÓLICA DE VALPARAÍSO
FACULTAD DE INGENIERÍA ESCUELA DE
INGENIERÍA QUÍMICA.

● Dispersiónde contaminantes atmosféricos:

Modelo gaussiano de la columna de humo.
Universidad Pablo de Olavide.
● María
E Manzur * a , Graciela Benzal b , Silvia N
González c(2012) MODELO DE DISPERSIÓN
DE CONTAMINANTES ATMOSFÉRICOS.