FACTORES DE EMISIÓN

Un factor de emisión es una relación entre la cantidad de contaminante emitido a la atmósfera

con una unidad de actividad asociada a dicha emisión, como por ejemplo la cantidad de material
procesado o la cantidad de combustible usado.

En este caso la emisión se obtiene multiplicando el factor de emisión dado por la cantidad de
material o combustible procesado:

E=AxF

Donde:
E = Estimado de emisión para la fuente (a nivel de proceso)
A = Nivel de actividad (por ejemplo material producido)
F = Factor de emisiones controladas (por ejemplo, Kg de contaminantes emitidos/t de material
procesado)
FACTORES DE EMISIÓN

Caldera de es ng asador Manufactu a de

g s natural e vapor baterías

1
kg de contaminante/ kg de contaminante/
•••
••
kg de contaminante/
1O m de gas quemado hr/ n de area superficial
del desengrasador 1if baterias

Ejemplos de factores de emisión

 FACTORES DE EMISIÓN
Ejemplo:

Calcular las emisiones anuales de NO2 de una caldera sin sistema de control en una planta
termoeléctrica (consumo de calor>100 millones de Btu/h) que quema gas natural.

El consumo anual de gas natural es de 40 millones de metros cúbicos.

Desarrollo:

FNOX = 8800 Kg/106 m3 (Tomado del AP-42. Tabla 1.4-2)

Qcon = 40 x 106m3/año
ENOX = FNOX x Qcon
= 8800 x 40
= 352000 Kg/año
 MODELAMIENTO DE LA DISPERSIÓN
ATMOSFÉRICA DE CONTAMINANTES
CONTAMINANTES EMITIDOS
DIFUSIÓN/
DISPERSIÓN

• REDUCCIÓN DE LA
CONCENTRACIÓN

• ARRASTRE
 MODELAMIENTO DE LA DISPERSIÓN
ATMOSFÉRICA DE CONTAMINANTES
La Difusión/Dispersión de contaminantes no representa en realidad su
eliminación, sino un reparto en volúmenes de aire cada vez mayores.
PERFIL VERTICAL DE TEMPERATURA
La dispersión de los contaminantes está estrechamente
relacionada con el movimiento turbulento vertical de la
atmósfera, la cual depende del perfil vertical de temperatura.

Aire caliente

N,o rmal lniVersión térmica

En Meteoro/agio Aplicada a la Contamlnac/6n Atmosf'rlca, las alturas de lnteris

raramente superan los 2 km
7
Formulación matemática de la dispersión de Gauss

Lih
C{x,y,z) o eti o del d
he ¿e . ,,.
la
e t acmn e
hs u p o (x.y.z .

y
Formulación matemática de la dispersión de Gauss

Se ha tomado como base el algoritmo Gaussiano del modelo ISC3,

desarrollado en los años ’60 a partir de investigaciones sucesivas de
Turner (1964 y 1970), Gifford (1960) y Pasquill (1974) entre otros,
siendo la ecuación general del modelo la siguiente:

𝐶𝑒𝑠𝑡 =
𝑄× 106
×𝑒
[ (-)] × [[
𝑦 2
−0.5 𝜎
𝑦
𝑒
𝑧−ℎ𝑒 2
(-)] + 𝑒 [ (-)]]
−0.5 𝜎
𝑧
𝑧+ℎ𝑒 2
−0.5 𝜎
𝑧
2𝜋𝜎𝑦 𝜎𝑧 𝑢𝑠
Formulación matemática de la dispersión de Gauss

𝐶𝑒𝑠𝑡 =
𝑄 × 10
×𝑒
6
[
−0.5 (-)] × [[
𝑦 2
𝜎𝑦
𝑒
−0.5
𝑧−ℎ𝑒 2
(-)] + 𝑒 [ (-)]]
𝜎𝑧
−0.5
𝑧+ℎ𝑒 2
𝜎𝑧
2𝜋𝜎𝑦 𝜎𝑧 𝑢𝑠

Donde:
𝜇𝑔
𝐶𝑒𝑠𝑡 : 𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 𝑒𝑠𝑡𝑖𝑚𝑎𝑑𝑎 (𝑎𝑝𝑜𝑟𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑓𝑢𝑒𝑛𝑡𝑒 𝑝𝑢𝑛𝑡𝑢𝑎𝑙) 𝑑𝑒𝑙 𝑐𝑜𝑛𝑡𝑎𝑚𝑖𝑛𝑎𝑛𝑡𝑒; 𝑒𝑛 3
𝑚
𝑔
𝑄: 𝑡𝑎𝑠𝑎 𝑑𝑒 𝑒𝑚𝑖𝑠𝑖ó𝑛 𝑑𝑒𝑙 𝑐𝑜𝑛𝑡𝑎𝑚𝑖𝑛𝑎𝑛𝑡𝑒; 𝑒𝑛
𝑠
𝜎𝑦 : 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑑𝑖𝑠𝑝𝑒𝑟𝑠𝑖ó𝑛 𝑙𝑎𝑡𝑒𝑟𝑎𝑙; 𝑒𝑛 𝑚
𝜎𝑧 : 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑑𝑖𝑠𝑝𝑒𝑟𝑠𝑖ó𝑛 𝑣𝑒𝑟𝑡𝑖𝑐𝑎𝑙; 𝑒𝑛 𝑚
𝑚
𝑢𝑠 : 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑣𝑖𝑒𝑛𝑡𝑜 𝑎 𝑙𝑎 𝑎𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎 𝑒𝑚𝑖𝑠𝑖ó𝑛; 𝑒𝑛
𝑠
𝑥 ∶ 𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎 𝑑𝑒𝑙 𝑟𝑒𝑐𝑒𝑝𝑡𝑜𝑟 𝑐𝑜𝑛 𝑟𝑒𝑠𝑝𝑒𝑐𝑡𝑜 𝑎 𝑙𝑎 𝑓𝑢𝑒𝑛𝑡𝑒 (𝑥 = 0) 𝑒𝑛 𝑑𝑖𝑟𝑒𝑐𝑐𝑖ó𝑛 𝑎 𝑠𝑜𝑡𝑎𝑣𝑒𝑛𝑡𝑜; 𝑒𝑛 𝑚
𝑦 ∶ 𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎 𝑑𝑒𝑙 𝑟𝑒𝑐𝑒𝑝𝑡𝑜𝑟 𝑐𝑜𝑛 𝑟𝑒𝑠𝑝𝑒𝑐𝑡𝑜 𝑎 𝑙𝑎 𝑓𝑢𝑒𝑛𝑡𝑒 (𝑦 = 0) 𝑝𝑒𝑟𝑝𝑒𝑛𝑑𝑖𝑐𝑢𝑙𝑎𝑟 𝑎𝑙 𝑒𝑗𝑒 𝑑𝑒 𝑙𝑎 𝑝𝑙𝑢𝑚𝑎; 𝑒𝑛 𝑚
𝑧 ∶ 𝑎𝑙𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑟𝑒𝑐𝑒𝑝𝑡𝑜𝑟 𝑐𝑜𝑛 𝑟𝑒𝑠𝑝𝑒𝑐𝑡𝑜 𝑎𝑙 𝑛𝑖𝑣𝑒𝑙 𝑑𝑒𝑙 𝑠𝑢𝑒𝑙𝑜 (𝑏𝑎𝑠𝑒 𝑑𝑒 𝑙𝑎 𝑐ℎ𝑖𝑚𝑒𝑛𝑒𝑎); 𝑒𝑛 𝑚
ℎ𝑒 : 𝑎𝑙𝑡𝑢𝑟𝑎 𝑒𝑓𝑒𝑐𝑡𝑖𝑣𝑎 𝑑𝑒 𝑙𝑎 𝑝𝑙𝑢𝑚𝑎; 𝑒𝑛 𝑚 (𝑑𝑒𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑎𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎 𝑐ℎ𝑖𝑚𝑒𝑛𝑒𝑎 (ℎ𝑠 ) 𝑦 𝑙𝑎 𝑒𝑙𝑒𝑣𝑎𝑐𝑖ó𝑛
𝑑𝑒 𝑙𝑎 𝑝𝑙𝑢𝑚𝑎 (∆ℎ) )
Formulación matemática de la dispersión de Gauss
x: dato de campo
y: dato de campo
z: dato de campo
Q: dato de campo (medición) o estimado según factores de
emisión¿¿¿???
us: medido o estimado ¿¿¿¿????
he: calculado ¿¿¿???
y: calculado¿¿¿??
z: calculado¿¿¿??
Formulación matemática de la dispersión de Gauss
Algunos supuestos del Modelo:

• Los contaminantes emitidos son continuos o, al menos, las

emisiones ocurren por un intervalo de tiempo que es mayor que
el tiempo que demora un contaminante en llegar al receptor
desde punto de emisión.

• No hay ocurrencia de reacciones químicas.

• Se considera que los gases son de baja densidad como los gases
de combustión, o en el caso de partículas, se asumen como más
pequeñas que 20 µm de diámetro aerodinámico
(aproximadamente), pues estas suelen permanecer mayor
tiempo en suspensión.
Formulación matemática de la dispersión de Gauss
Algunos supuestos del Modelo:

Considera un estado estacionario, es decir no depende ni es

función del tiempo.
- La difusión de la masa es despreciable en el eje x.
- La velocidad del viento se considera constante y las variaciones
de esta en los tres ejes coordenados se consideran despreciables.
- El modelo no puede trabajar bien en condiciones de calma donde
la velocidad del viento es cercana a cero.
Coeficientes de dispersión
 Coeficientes de dispersión
2.-------,,-------.,----,-----,---,--------,,------,----,-----.----,
...
~

P, Difer en cia de t emp er atura

t;

'-4-- - - - - , M
..
1,
t;
E
.Y "' ~
'le
e' \ Gradiente v erti ca l
•o 1
·o
n:i
>
·\ r adiabá ti co seco

g = 1K/100 m
(L1 ~-
't;
w ~

.
~

Gradi ente \
v erti ca l amb iental D Estable
z \ F Isoterma
10 20 30 4D 50
Inestable \ E
Te mpe ratura oc \ Inversión
C
B

T
Gradiente de temperaturas según la clase de
atmósfera
 Coeficientes de dispersión
Clases de estabilidad atmosférica (Pasquill-Gifford-Turner)
Día Noche*
Velocidad de Radiación solar entrante Cobertura de nubes**
Viento superficial (a
Fuerte Moderado Leve 𝟒 𝟑
10 m), m/s ≥ ≤
(>700 W/m2) (350-700 W/m2) (<350 W/m2) 𝟖 𝟖
<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 después de la Categoría Estabilidad

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

16
 Coeficientes de dispersión
Una de las primeras expresiones matemáticas fue dada aproximadamente
en 1967 (Smith), quien mostró mediciones horarias a distancias de hasta
10 km desde una fuente de una altura de 108 metros. De acuerdo con
Smith, las mediciones correspondían a las siguientes expresiones:
Coeficientes de dispersión
• -i,,4 ~-~---r---.--~~~---.------.----r-"T""""T"........-,r-r------r-----r-----r-""T"""T""T""T"TII

A - : .· tet-em . ly uns:ttable
i8 - M:ode[·at ly 1n1sta b~.e
C - Sligb-Uy unstable
D - N eutra]
1: ·-· Sl.ighny stab~.e

.' - M:oderateJy stabl

Di,stanc from souro 1n

Coeficientes de dispersión

UJ)

E
ti,
~·
....
_ if'.J,

m
''¡.J ,
:r;-
~ UJ:l
~.
f,J•

§
''Gj

t
·~-
"CI
'¡j
y,
't:;
:. rn . - Ex· ,e r,etuely unstabl
~
B - Moderalely unstabl
C -· SUghtJy unstabl.e
D - Neutral
E - ligbUy stab]e
F - Moderalely table
m ~-~~-~~~--~~~~~-~~~~~
rn~ 103 10,
DistaH1ce frou1 · ouro m
 Coeficientes de dispersión

Paramd.ers
Stability conditions a b e d
A OAO 0,91 0.41 0.91
B1 0.36 0.86 0.33 0.86
e 0.32 0.78 0.22 0.78
< D 0.31 0.71 0.06 0.71
 Coeficientes de dispersión
Una de las maneras más simplificadas y largamente estudiadas y
aceptadas por la EPA (sobre todo para áreas urbanas) son las funciones
Sigma establecidas por Briggs en 1973; cuentan con una gran validez para
distancias a sotavento (x) de entre 100 m y 10 km. Las funciones Sigma de
Briggs para áreas urbanas también son llamadas las funciones Sigma
McElroy-Pooler, pues derivaron de una serie de experimentos de la
dispersión urbana realizadas por el investigador del mismo nombre en
1968.
 Coeficientes de dispersión

FUNCIONES SIGMA DE LA DISPERSIÓN PARA ÁREAS URBANAS

Grados de Estabilidad 𝝈𝒚 (m) 𝝈𝒛 (m)
de Pasquill
A–B 0.32(x)(1+0.0004x)-0.5 0.24(x)(1+0.001x)0.5
C 0.22(x)(1+0.0004x)-0.5 0.20x
D 0.16(x)(1+0.0004x)-0.5 0.14(x)(1+0.0003x)-0.5
E–F 0.11(x)(1+0.0004x)-0.5 0.08(x)(1+0.00015x)-0.5
 Coeficientes de dispersión

FUNCIONES SIGMA DE LA DISPERSIÓN PARA ÁREAS RURALES

Grados de Estabilidad 𝝈𝒚 (m) 𝝈𝒛 (m)
de Pasquill
A 0.22(x)(1+0.0001x)-0.5 0.20x
B 0.16(x)(1+0.0001x)-0.5 0.12x
C 0.11(x)(1+0.0001x)-0.5 0.08(x)(1+0.0002x)-0.5
D 0.08(x)(1+0.0001x)-0.5 0.06(x)(1+0.0015x)-0.5
E 0.06(x)(1+0.0001x)-0.5 0.03(x)(1+0.0003x)-1
F 0.04(x)(1+0.0001x)-0.5 0.016(x)(1+0.0003x)-1
Velocidad de Viento a la Altura de la Emisión
p ,

h,s
l s == U ref

Donde:

Altura de la ch'ime1 nea (m)

p Exp,orierite die la formula de potencia que estima la velocidad
del viento (adimens.·onal)
Us La velocidad del v·ento (mis) a la alltura die la ch·menea
u.-nef Ve ocidladl del v·ento medida a la altura die referencia (mis)
-
,_, l'€,J/' Altura de reierencia (m)
Velocidad de Viento a la Altura de la Emisión
Los valores de P pueden ser establecidos por el u1suario en función de la clase de estabilidad y die
la clase debido al terreno, el modelo tiene los sig1uientes valores predeterminados:

Valores dell exponente P para calcular la velocidad del viento

Categorías de Exponente Exponente

estabilidad (condiciones (condiciones
111rales) urbanas)
A 0 .07 0.15
B 0.07 0.15
e 0.10 0.20
D 0.1~ 0.25
E 0.3 '" 0.30
F 0.55 0.30
 CÁLCULO DE ELEVACIÓN DE LA PLUMA
ℎ𝑒 = ℎ𝑠 + ∆ℎ

Una de las principales variables relacionadas a la elevación de la pluma y la

consecución de una buena dispersión (máxima elevación de la pluma).

Existen órdenes de magnitud ya más o menos establecidas para las

velocidades de salida de los gases (Vs):

1. Sistemas de Tiro Natural Vs = 6 m/s.

2. Sistemas de Corriente Forzada:

• Chimeneas de hasta 20 m de altura.- Vs = 6 m/s.
• Chimeneas de altura comprendida entre 20 y 40 m.- Vs = 9 m/s.
• Chimeneas de altura superior a 40 m de altura.- Vs = 12 m/s.
 CÁLCULO DE ELEVACIÓN DE LA PLUMA

Asimismo, usualmente se considera por temas estructurales y visuales que

las chimeneas cumplan la siguiente relación entre su altura (hs) y su diámetro
externo medio (Dem)

ℎ𝑠
= 𝑒𝑛𝑡𝑟𝑒 15 𝑦 20
𝐷𝑒𝑚
 CÁLCULO DE ELEVACIÓN DE LA PLUMA

Fases de Evolución del Penacho

a) Primera fase:
El penacho es prácticamente vertical, ya que inicialmente posee un
movimiento ascendente debido a.- la velocidad de salida y al efecto de
flotación cuando la temperatura de salida es mayor a la del aire.
El final de esta fase vendrá marcado, consecuentemente, cuando las
velocidades y temperaturas del penacho se asemejen a las del aire.

b) Segunda fase:
Al perder su impulso ascensional, el penacho va curvando progresivamente
dominado por el viento horizontal. Concluye esta fase cuando la componente
vertical de la velocidad del penacho tiene el mismo orden que la de la
turbulencia atmosférica.
 CÁLCULO DE ELEVACIÓN DE LA PLUMA

Fases de Evolución del Penacho

3) Tercera fase:
La emisión tiene una densidad y una cantidad de movimiento semejantes a
las del aire atmosférico. El penacho ya no tiene movimiento propio , sino que
es arrastrado y su movimiento responde al viento.
 CÁLCULO DE ELEVACIÓN DE LA PLUMA

La altura física modificada de la emisión (hs') (Efecto downwash)

En prim er luga r, deberemos hall ar una variab le !llamada: La alt ura fís ica modifi cada de la em isión (h s' }. Co n el
f in de tener en cuenta el efect o de tu rbulencia (l avado ) a la altura de la ch im enea (hs), la m od ificació n de la
altu ra física de la ch imenea se rea li za sigu iendo la for mulac ión brindada por Bri ggs (1 974 ).

h~ = hs + 2ds e:-1.5) para Vs < 1. Sus

 CÁLCULO DE ELEVACIÓN DE LA PLUMA

Para hallar el Parámetro de Flujo de Bouyancia 'Fb' (Factor de Flotabilidad)

Para efectos de l cálculo de la altura efect iva de la pluma (he) será necesa ri o ca lcu lar previamente este
parámetro:

Para efectos de la presente gu ía, cuando se calcu le 'he' obviaremos los pos ibles casos en los que Ts es me no r o
igual que Ta. Es decir la aplicab illidad del algoritmo constru ido estará dada so lamente pa ra casos cuando
Ts > Ta .
 CÁLCULO DE ELEVACIÓN DE LA PLUMA

• Pa ra los casos en que la at1

m ósfera se encuentr a inestable o neutr a (A, IB,. C, ID):

Cuando Fb < SS
(Xm = 49F 5/8) F o.75)
he = h~ + 21.425 ( b
Us
..·

Cuando Fb > SS

(Xm = 119F 2/5)

 CÁLCULO DE ELEVACIÓN DE LA PLUMA
• Pa ra los casos en que la atmósfera se encuentra estab le (E, F) :

Para l as situaciones de atmósfera estab le, pri mero se debe calcu lar el pará1metro de estabil idad ('sp')

ae [Xm = 3.14us/(sp1/2)]
sp == g
az
Ta

ae
-az ➔ es la gradiente de temperatura, la cua l se considera como una ap roximación por defecto:

ae
Para grado de estab ilidad E (o 5) ➔ -=0.020 K/m
az
ae
Para grado de estab ilidad F (o 6) ➔ "a;'=0.035 K/m

Luego,
CÁLCULO DE ELEVACIÓN DE LA PLUMA
 OTROS ASPECTOS

Sorne of the fe atures of ADMS

fluctuations

plume
visibility

r a di oacti ve
decay and

time
varying
t changes in
gamma dose

emission
dispersion sutface
around '\ roughne~
buildinqs~