INSTITUTO POLITÉCNICO NACIONAL - Alexander Gelbukh

A pesar de esta variedad de modelos y paquetes ninguno se utiliza para resolver
problemas inversos como sería la determinación de posibles fuentes
contaminantes a partir de los datos de dispersión de estos sobre un área definida.
2. Modelos atmosféricos aplicados a la ciudad de México
Para la ciudad de México se han aplicado varios modelos como son los
meteorológicos, los de dispersión, los químicos etc. Su aplicación ha dependido
del tipo de estudio requerido y de los contaminantes de interés.
Los modelos HOTMAC y RAPTAD desarrollados en los Laboratorios Nacionales
de Los Álamos, Nuevo México, se aplicaron para conocer la distribución de las
partículas suspendidas en la atmósfera de la zona metropolitana del valle de
México. Se usó también el modelo tridimensional CIT utilizado en los Estados
Unidos de Norte América, para simular la dispersión de los contaminantes
químicos atmosféricos sobre el Distrito Federal.
El modelo ISCST3 que significa Industrial Source Complex Short Term se usó en
el Instituto Mexicano del Petróleo [Technical manual, 1995a, 1995b]. La base del
modelo es la ecuación de la pluma gaussiana de estado estable y lineal que se
usa con algunas modificaciones para modelar emisiones provenientes de
chimeneas.
Este modelo acepta datos meteorológicos horarios para definir las condiciones de
elevación de la pluma, transporte, difusión y deposición. El modelo estima la
concentración o el valor de deposición para cada combinación de fuente y receptor
en forma horaria y calcula los promedios a corto plazo. El modelo tiene opciones
para modelar emisiones provenientes de un amplio rango de fuentes que puedan
estar presentes en un complejo industrial.
3. Modelo matemático utilizado en este trabajo
El modelo usa la ecuación de estado estable de pluma gaussiana para una fuente
continua elevada. Para cada fuente y cada hora, el origen del sistema coordenado
de las fuentes se coloca a nivel de la superficie en la base de la chimenea. El eje
X es positivo en la dirección viento abajo, el eje Y está normal al eje X y el eje Z se
extiende verticalmente. Las ubicaciones fijas de los receptores son convertidas
para cada sistema coordenado de la fuente, para cada cálculo horario de
concentración.
Las concentraciones horarias calculadas para cada fuente en cada receptor son
sumadas para obtener la concentración total producida en cada receptor por la
combinación de las emisiones de las fuentes.
14