UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

DECANA DE AMÉRICA

Trabajo de Investigación

SOFTWARE CALPUFF VIEW 9.0

 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

ÍNDICE
ÍNDICE
Introducción Aplicaciones

Objetivos
Conclusiones

Descripciónn

Bibliografía
Características
s

Alcancese

Métodos
utilizados
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

INTRODUCCIÓN
Los modelos de dispersión de contaminantes atmosféricos usan una serie de técnicas matemáticas
para simular los procesos físicos y químicos que afectan a los contaminantes cuando se dispersan
como consecuencia de los factores meteorológicos y topográficos del emplazamiento; o reaccionan
en la atmósfera.
Los modelos de dispersión ayudan a predecir el comportamiento de los contaminantes y mediante
simulación prevenir que estos puedan dañar a las personas y el medio ambiente.
El objetivo de los modelos de dispersión de contaminantes son encontrar una relación entre las
sustancias que son emitidas a la atmosfera, las concentraciones de la sustancia y otras que se
originan a partir de las mismas en la atmosfera.
A través de modelos basados en experimentos de laboratorio se ha podido identificar mecanismos
de transporte atmosférico, pero el problema global a sido abordado por modelos matemáticos que
describen aspectos físicos y químicos los que requieren una calibración específica sobre cada
problema. Estos modelos se tratan de manera aproximada resolver las ecuaciones de difusión que
gobiernan el movimiento de los contaminantes en el flujo turbulento de la atmosfera.
Los métodos de solución se subdividen en: Modelos de caja, Modelos gaussianos y modelos de
diferencias finitas que enfocan la solución a la ecuación de difusión atmosférica utilizando diferentes
planteamientos físicos: Euleriano y lagrangiano.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

CALLPUFF VIEW 9.0

Modelizador medioambiental para la evaluación de la dispersión tipo Puff
(bocanada o soplo).
CALPUFF View es una interfaz basada en GIS diseñada para el modelo de
referencia CALPUFF aprobado por la US EPA. Se compone de tres módulos:
•CALMET. Un modelo de diagnóstico meteorológico en 3-D que acepta campos
de viento preprocesados.
•CALPUFF. Es uno de los Modelos de dispersión de aire más avanzados que
existen en la actualidad. Un modelo tipo soplo (PUFF) que ofrece información
completa hora a hora de la variación espacial del viento y su estabilidad.
Contiene algoritmos adicionales que le permiten emular al potente AERMOD en
distancias cortas donde los modelos PUFF son generalmente menos fidedignos.
Es razonablemente preciso en distancias no superiores a los 300 Km
•CALPOST. Programa de Postprocesado que compila resultados de CALPUFF
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Objetivos:
• Conocer la metodología langrangíana de dispersión
atmosférica CALPUFF.
• Conocer el rango de la aplicación de la metodología.
• Conocer tiempo de aplicación y áreas de cobertura de la
metodología
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Descripción
CALPUFF View es una interfaz basada en GIS diseñada para el modelo de referencia
CALPUFF aprobado por la US EPA. Se compone de tres módulos:
 CALMET. Un modelo de diagnóstico meteorológico en 3-D que acepta campos de viento preprocesados.
 CALPUFF. Es uno de los Modelos de dispersión de aire más avanzados que existen en la actualidad. Un
modelo tipo soplo (PUFF) que ofrece información completa hora a hora de la variación espacial del viento y
su estabilidad. Contiene algoritmos adicionales que le permiten emular al potente AERMOD en distancias
cortas donde los modelos PUFF son generalmente menos fidedignos. Es razonablemente precision en
distancias no superiores a los 300 Km.
 CALPOST. Programa de Postprocesado que compila resultados de CALPUFF
CALPUFF View versión 9.0 ahora es compatible con los siguientes ejecutables del sistema
BETA Versión 7 lanzados el 6 de febrero de 2020:
CALPUFF Version 7.3.2, Level 200113
• CALPOST Version 7.2.0, Level 150720
• POSTUTIL Version 7.1.0, Level 190930
• CALSUM Version 7.1.0, Level 190131
UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Características
Primera

REDUCCIÓN DEL TIEMPO EN EL PROYECTO

Descarga y procesa de forma
automática los datos
Meteorológicos y geográficos
de la web (gratuitamente)

• Permite la
introducción
gráfica de
mallados ,
fuentes,
receptores y
líneas de costa.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Segunda

ADMITE TODO TIPO DE FUENTES

Puntuales, de línea, volumen y
área. Además de fuentes de
emisión variable o
constante. Emisiones inestables
y condiciones meteorológicas:

• Campos de variables meteorológicas (vientos,

temperatura) en enmallados en 3-D
• Campos de variables espaciales combinando:
altura , fricción de la velocidad, cociente de
convección y de precipitación.
• Cocientes de turbulencia de variación vertical y
horizontal Rayos de dispersión.
• Datos de emisiones y fuentes variables con el
tiempo.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Tercera SOPORTE PARA DATOS DE USO DE SUELO

SHAPELIFE
• Los datos de uso de la tierra del
archivo Shapefile, incluidos los
datos de uso de la tierra de Geo
Base disponibles en el archivo
FTP Geogratis de Natural
Resources Canada, ahora se
pueden importar a Creador de
uso de la tierra utilidad para
generar archivos de datos de uso
del suelo compatibles con el
procesador geofísico de
CALPUFF View. Esto incluye una
utilidad para que los modeladores
definan códigos de uso de la
tierra para cada atributo shapefile.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Cuarta
CALCULA EL OZONO MENSUAL A PARTIR DEL
FONDO POR HORA
• Cuando se proporciona un
archivo de ozono de fondo
cada hora, los datos en la
pestaña O3 mensual de las
opciones de
Transformación química
se pueden calcular
utilizando los datos del
propio archivo. En este
caso, los datos de la
pestaña O3 mensual se
utilizarán para llenar los
vacíos en el archivo de
datos por hora. .
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Alcances
Alcance indicado en el trabajo
 Este método se aplica a pequeñas y medianas empresa además de aplicar en
ciudades y regiones cubriendo un área de 5x5 km2 a 3360x3360 km2.
 Este método a sido aplicado en:
 Villasenor et al. (2003) estudiaron la dispersión del polvo procedente de ciudad
México.
 Cohen et al. (2005), describieron las emisiones provenientes de vehículos
automotores en Portland, Oregón.
 Taseko (2008) lo emplearon en una zona minera en Canadá.
 Song et al, (2010) analizaron las emisiones de buques en un puerto italiano de la
costa del mar jónico.
 Maclontosh et al.( 2010) lo utilizaron en la industria metalúrgica.
 Poplawski et al. (2011) Investigaron cruceros en West Virginia.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Métodos utilizados
Método Lagrangiano
• Existe un gran número y variedad de alternativas en la representación del flujo de contaminantes en la atmosfera, lo que
sugiere que no existe un único problema, sino una combinación de fenómenos, y que la elección de una u otra solución
depende fundamentalmente de la escala del espacio y tiempo escogido, que condicionan la consideración de estos
fenómenos.
• El planteamiento general sigue el esquema de un modelo de procesos o fenómenos, es decir, consisten en aplicar y sumar las
contribuciones de todos los modelos de proceso en una ecuación de balance o conservación para cada propiedad considerada.

• Donde el primer término representa la variación de concentración del contaminante en el tiempo.

• El segundo, la entrada neta del contaminante asociada al flujo atmosférico (advección y turbulencia).
• El tercero, la variación de concentración de contaminante debido a la difusión molecular.
• El cuarto, la velocidad de remoción del contaminante por reacción química.
• El quinto, la aparición de contaminante en el entorno considerado, proveniente de focos emisores en el propio entorno.

• Sin embargo, la resolución analítica de estas ecuaciones no es posible, puesto que no se conoce forma de función ni de la difusión de
turbulencia asociada. Además hay que tener en cuenta que es necesario resolver las N ecuaciones diferenciales simultáneamente, para
tener en cuenta la influencia de la concentración de un contaminante en la velocidad de reacción de los otros contaminantes en su
balance.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
DECANA DE AMÉRICA

Uno de los modelos lagrangianos más versátiles y mas aplicado en años recientes es el
CALPUFF que es un modelo no estacionario, multicapas y multiespecies de dispersión de pufs,
que puede simular los efectos, temporales y espaciales, de las condiciones meteorológicas en el
transporte, la transformación y la deposición de contaminantes en la atmósfera.
El modelo CALPUFF puede utilizar los campos meteorológicos generados por el modelo
CALMET o en su defecto, las medidas de estaciones meteorológicas.
El modelo CALPUFF contiene algoritmos que consideran diversos fenómenos.
• Efectos en la cercanía de los focos emisores.
• El flujo descendente tras obstáculos.
• Penetración parcial del penacho sobre la capa de mezcla.
• La interacciones del terreno a escala submalla.
• Los efectos a largo plazo.
• La deposición de contaminantes (lavado húmedo y deposición seca)
• La transformación química.
• La cizalladura vertical del viento, las interacciones costeras y las singularidades del transporte
sobre masas de agua.
• También puede simular variaciones arbitrarias de las emisiones y contiene diferentes opciones
para el tratamiento de los procesos físicos de los contaminante en la atmosfera, en diferente nivel
de detalle en función de la aplicación del modelo.
El modelo CALPUFF, ha ido recomendado por el grupo interinstitucional sobre modelos de calidad del aire
(IWAQM), para la caracterización del transporte y dispersión de contaminantes.

El modelo CALPUFF es útil para diversos campos mostrados a continuación

Industria
Química

Metalúrgica Agricultura

CALPUFF

Industria
Minería automotriz e
Hidrocarburos

Plantas de
mantenimiento
 Aplicaciones
MODELADO CALPUFF-WRF DE LA DISPERSIÓN DE PMX, NOX Y SO2 EMITIDOS
POR LA TERMOELÉCTRICA DE MIRAFLORES EN EL CANAL DE PANAMÁ
Tipo de Artículo: Original
Fecha de Recepción: 27 de febrero de 2017
Fecha de Aceptación: 25 de septiembre de 2017

Se escogió el modelo CALPUFF por ser el recomendado por la EPA para recorridos mayores de 50 km y
donde se presentan cambios de cobertura: agua, suelo. Este modelo, además de ser complejo, requiere
el uso de datos continuos en tiempo y espacio para generar resultados válidos. A falta de una
cobertura total en estaciones meteorológicas superficiales, recurrimos al uso del modelo meteorológico
WRF con el cual se pudo simular el campo de vientos tridimensional para un dominio de 50x50km2 con
resolución de 4x4km2. Los resultados obtenidos que se presentan en este artículo muestran los
puntos de mayor concentración promedio de 1h de contaminantes, cuyos valores aplicando la
metodología EPA para Puerto Rico como referencia, arrojan buenos índices de calidad del aire.

• Esto conlleva, un método semiempírico donde las ecuaciones matemáticas de la dispersión se calibran con
observaciones para que el modelo genere resultados. En esta investigación se utilizó el modelo de dispersión de
licencia libre más reciente aprobado por la EPA: CALPUFF v.5.8.5
• El objetivo fue modelar la dispersión y el transporte de material particulado (PMx), óxidos de nitrógeno (NOx) y dióxido
de azufre (SO2) que son los contaminantes típicos de la combustión fósil.
• CALPUFF requiere un modelo meteorológico CALMET, la caracterización de la fuente de emisión para calcular la
dispersión y el transporte con CALPUFF y finalmente, el cálculo de concentraciones en la zona de estudio con CALPOST.
 Caracterización de la Planta Termoeléctrica de Miraflores como fuente fija de emisión
Tabla 1. Motores operativos en la Planta Termoeléctrica
Al momento de esta investigación, la misma contaba con 8 motores (ver tabla
de Miraflores al momento de la investigación
1), unos a diésel y otros a búnker, de los cuales se contaba con el informe de
emisiones de los motores 6, 7 y 8 medidas en febrero de 2012 (ver tabla 3).
El informe de 2013 estaba bajo revisión aún. Las unidades 6, 7 y 8 pertenecen
a un motor de combustión interna que utiliza búnker para su operación y diésel
para el arranque y parada de la unidad La ubicación geográfica de las
chimeneas de estos motores se muestra en la figura 1 y en la tabla 2 .
Figura 1. Localización de las chimeneas de los motores 6, 7 y 8 de la Planta Termoeléctrica de
Miraflores. Imagen de Google Earth.
 Caracterización de la Planta Termoeléctrica de Miraflores como fuente fija de emisión
Tabla 2. Coordenadas UTM WGS84 de las chimeneas de los motores 6 Tabla 4. Datos de la Chimenea Equivalente como Fuente Fija de
,7 y 8 de la Planta Termoeléctrica de Miraflores Emisión para el Modelado CALPUFF

Tabla 3. Emisiones de PMx, NOx y SO2 en los motores 6, 7 y 8 de la

Planta Termoeléctrica de Miraflores medidas del 27 al 29 de febrero de
2012

*estimación en sitio, de la base respecto al NMM.

**valor calculado por el programa CALPUFF utilizando el MDT.
Para la modelación CALPUFF se reemplazaron las tres chimeneas por
una sola equivalente cuya posición se calculó en el centroide del
triángulo formado por las chimeneas 6, 7 y 8. Esta chimenea
equivalente tiene la misma configuración geométrica de las
chimeneas originales. El factor de ponderación para cada motor será
igual a su caudal medido entre el caudal total emitido por los tres
motores. La suma de los productos de cada factor de ponderación
por la concentración medida del contaminante dará como resultado
su emisión ponderada.
 Modelo Meteorológico CALMET
DATOS DE ENTRADA
Para la investigación se escogió un dominio de 50x50 km2, Como Modelo Digital del Terreno se utilizó el SRTM1 (Shuttle
cuyo centro se ubicó aproximadamente en el centro del Canal Radar Topography Mission Global Coverage) Versión 3, cuya
de Panamá (9°06’44”N, 79°45’55”W, WGS84) con un precisión es de 30 m (aproximadamente 1”). Este modelo
cuadriculado de 4x4 Km2, 35 niveles verticales con el nivel más contiene información sobre las elevaciones utilizando el
bajo a 20 m del suelo. El modelo WRF generó, a partir de WGS84 como referencia geodésica con una cobertura global
datos reales (medidos y de análisis) una meteorología entre 56° de Latitud Sur a 60° de Latitud Norte (ver figura 3).
simulada tridimensional para el dominio de estudio,
correspondiente al año 2014, del 1 de enero a las 0:00 h hasta
el 31 de diciembre a las 23:00 h. (ver figura 2).
Figura 2. Dominio para el modelo meteorológico WRF.

Figura 3. Modelo del Terreno extraído del SRTM1 para el dominio

de estudio en una cuadrícula de 1x1 km2 de resolución espacial.
 Modelo Meteorológico CALMET
Para el uso del suelo y cobertura se utilizó el GLCC (Global
Land Cover Characterization), específicamente los datos de
Uso del Suelo y Cobertura para América del Norte. Este modelo
de uso del suelo y cobertura global tiene una resolución de 1
km (30” aproximadamente), (ver figura 4).

Figura 5. Modelo Meteorológico CALMET donde se muestra la

dirección del viento en cada celda para el dominio de estudio en una
cuadrícula de 1x1 km2 de resolución espacial.

Figura 4. Modelo de uso del suelo y cobertura extraído del GLCC para
el dominio de estudio en la cuadrícula de 1x1 km2 de resolución espacial.
3.2 PROCESADO Y RESULTADOS
Con los datos de entrada geodésicos (modelo digital del
terreno, uso del suelo y cobertura) y los datos del modelo WRF
postprocesados para que fueran compatibles con CALMET, se
especificó una cuadrícula de 1x1 km2 para generar el campo
de vientos tridimensional que se ajuste a la forma del terreno Figura 6. Modelo Meteorológico CALMET muestra la variación
y al tipo de cobertura. Los resultados del procesado CALMET tridimensional de los vientos tanto en dirección como en velocidad,
se muestran en las figuras 5 y 6. entre otras variables meteorológicas.
 Modelo de dispersión CALPUFF
Tabla 2. Coordenadas UTM WGS84 de las chimeneas de los motores
ÁREA DE MUESTREO Y CUADRÍCULA DISCRETA 6 ,7 y 8 de la Planta Termoeléctrica de Miraflores
Para ello, CALPUFF pidió el rango de pixeles a considerar en
las direcciones X e Y. El dominio total fue de 50 pixeles de 1x1
km2 en cada dirección, se le indicó al programa que
considerara como cuadrícula de muestreo del 2 al 50 en
ambas direcciones. Una vez definida la cuadrícula de
muestreo, ésta se transformó a una red de puntos llamados
receptores. Para ello, CALPUFF calculó las coordenadas del Tabla 3. Emisiones de PMx, NOx y SO2 en los motores 6, 7 y 8 de la
centro de cada pixel transformando la cuadrícula a discreta. Planta Termoeléctrica de Miraflores medidas del 27 al 29 de febrero
de 2012
FUENTE DE EMISIÓN DE CONTAMINANTES A
MODELAR
Como fuente puntual fija se escogió la chimenea equivalente a
los motores 6, 7 y 8. El programa CALPUFF requiere como
datos de entrada: las coordenadas de la chimenea
equivalente, la elevación de la base respecto al Nivel Medio
del Mar (NMM), la altura de la chimenea, la temperatura
promedio del gas en la chimenea en Kelvin, el diámetro
interior en m y la velocidad del gas en la salida de la chimenea
en m/s.
CONTAMINANTES A MODELAR
Finalmente, se escogieron las especies a modelar: PMx, NOx,
SO2 y se introdujeron las emisiones estimadas en ton/año.
 Procesado y resultados CALPOST
CALPOST es el tercer componente modelador de CALPUFF, el Tabla 6. Concentraciones máximas promedio en 1h calculadas por
cual calculó la concentración promedio a nivel de terreno y el CALPUFF para la dispersión y transporte de PMx, NOx y SO2
flujo de deposición promedio de cada contaminante modelado emitidos por la Planta Termoeléctrica de Miraflores para el año 2014
para periodos específicos de tiempo (1 h, 3 h, 24 h, 1 año), a y su correspondiente Índice de Calidad de Aire calculado con los
partir de los resultados CALPUFF (ver tabla 5). parámetros de la EPA
CALPOST arrojó los resultados en tablas y gráficos (ver figuras
7 a 10) con el visor CALPUFF View v 8.4.0 de Lakes
Environmental Software.
Tabla 5. Concentraciones máximas calculadas por CALPUFF para la
dispersión y transporte de PMx, NOx y SO2 emitidos por la Planta
Termoeléctrica de Miraflores para el año 2014

*Calculado con el convertidor de μg/Nm3 a ppb para gases.

** Air Quality Index calculado a partir de las concentraciones en ppb
para los gases y en μg/m3 para las PMx.
 Procesado y resultados CALPOST

Figura 7. Concentraciones promedios en 1h para PMx, Nox y SO2 en el dominio de estudio en μg/Nm3. Gráficos de CALPUFF View v.8.4.0
montados sobre una imagen de Google Earth del área.

Figura 8. Concentraciones promedios en 3h para PMx, Nox y SO2 en el dominio de estudio en μg/Nm3. Gráficos de CALPUFF View v.8.4.0
montados sobre una imagen de Google Earth del área.
 Procesado y resultados CALPOST

Figura 9. Concentraciones promedios en 24h para PMx, Nox y SO2 en el dominio de estudio en μg/Nm3. Gráficos de CALPUFF View v.8.4.0
montados sobre una imagen de Google Earth del área.

Figura 10. Concentraciones promedios en 1 año para PMx, Nox y SO2 en el dominio de estudio en μg/Nm3. Gráficos de CALPUFF View
v.8.4.0 montados sobre una imagen de Google Earth del área.
 Procesado y resultados CALPOST

Figura 11. Pluma de dispersión de contaminantes provenientes de la Figura 12. Rosa de los vientos localizada en Planta Termoeléctrica de
Planta Termoeléctrica de Miraflores en el momento de mayor Miraflores. Muestra la distribución estadística de los vientos durante el
concentración promedio de 1h: 12 de octubre de 2014 a las 4:00 h. año 2014. Gráfico de CALPUFF View v.8.4.0 montado sobre una
Gráfico de CALPUFF View v.8.4.0 montado sobre una imagen de imagen de Google Earth del área.
Google Earth del área.
 Resultados Conclusiones

Para los tres contaminantes modelados, el comportamiento de • Los resultados obtenidos de la distribución de
dispersión fue similar. Las concentraciones promedio mayores concentraciones han sido de utilidad para determinar el
se presentaron en un período de 1h y fueron disminuyendo al área real de influencia de la Planta Termoeléctrica de
aumentar el periodo de modelado como era de esperarse. Miraflores e identificar los puntos receptores donde el
Siendo los valores mayores para las concentraciones de SO2 y programa calculó las mayores concentraciones de
los menores para el material particulado PMx. contaminantes.

• La fortaleza de este primer modelado radica en la

utilización de los datos del modelo meteorológico regional
WRF con el cual se logra un modelo meteorológico
tridimensional más completo, con la cual CALMET trabajó
para ajustar los datos a una Cuadrícula. Sin embargo, para
contaminantes que sufren transformaciones químicas, como
el NOx, será interesante ajustar los valores calculados por
CALPUFF tomando en cuenta esas transformaciones

Figura 13. Localización de la Fuente de Emisión Fija y los receptores

para mayor concentración de contaminantes modelados en periodos
de 1h, 3h, 24h y 1 año en 2014. Imagen de Google Earth del área de
estudio.
Conclusión
 CALPUFF es un software que facilita el trabajo de evaluación ambiental con el objetivo de conocer el
comportamiento de la nube de gases su dirección y la distancia que esta llegaría según las condiciones
geográficas de la zona donde se instale la industria.
 En los distintos estudios a demostrado ser el más eficiente y es por esto que es uno de los más utilizados a nivel
mundial, junto a esto también este método es aplicable a diferentes fuentes emisoras su rango de aplicación
también juega un papel importando su buen funcionamiento.
 CALPUFF ha sido empleado en un amplio rango de tamaño de región, desde aplicaciones locales de 5x5 km2
hasta regionales de 3360x33600 km2.
 Esto indica que permite abordar problemas más allá de la escala espacial para la cual fue diseñado.
 Existe una gran variedad para configurar CALPUFF. Se reportan ensayos desde los de alta resolución de 0.1 km
hasta alcanzar los 28 km en estudios regionales. En cuanto a los periodos de simulación evaluados, estos varían
desde 24 h hasta 4 años. Los niveles de altura por su parte se extienden desde 8 hasta 12 niveles.
 Debido a que CALPUFF tiene en cuenta los efectos del terreno y los flujos atmosféricos en zonas costeras, ha
sido ampliamente utilizado en dominios de terreno complejo y costeros
BIBLIOGRAFÍA
• https://www.addlink.es/images/pdf/2020/CALPUFF_View_Releas
e_Notes_V.9.0.pdf
• http://www.src.com/
• Guía de uso de modelo CALPUFF (5 Edición), 2000.
• Aplicación a escala urbana y validación del modelo
CALPUFF, Varsovia-Polonia, 2015.
• A distributed computing solution for CALPUFF
Environment-Canada, 2007.
• Validation of CALMET/CALPUFF models simulations around
a large power plant stack- España, 2015