Comparative evaluation of SO₂ in less favorable wind direction from a thermal power plant

Cezar García, Analí Machado, Neyma García, Carmen Cárdenas, William González, Jaime Martínez y Suher Yabroudi

Escuela de Ingeniería Química, División de Postgrado, Facultad de ingeniería, Universidad del Zulia, Maracaibo, Estado Zulia, Venezuela. E-mail: garcía@falcon.ing.luz.ve

Abstract

There is a thermal power plant located at Maracaibo State of Zulia, Venezuela. This plant represents an important source of sulphur dioxide emissions in the zone because of combustion of fuel oil. The SO₂ emissions from the chimneys are dispersed into atmosphere, generating high levels of SO₂ in the areas near to the plant that surpassed the limits for air quality.This work analyzes the distribution of SO₂ concentration toward the northen zone of the plant and in the less favorable wind direction. For the experiment, samples were collected by absorption through the RAC 5 system (model all weather shelter). For computational purposes the atmospheric dispersion models presented by Briggs and Pasquill-Gifford were applied. To apply the model hourly values of meteorological conditions were recorded at Maracaibo such as: wind speed, wind direction, cloudiness, insolation and temperature during December 2000 and January 2001. The results indicated that the Briggs model for urban dispersion fit quite well experimental values from the sampling zone, which takes into account local meteorological conditions and data from the pollution source.

Key words: SO₂ dispersion, air pollution, dispersion model, atmospheric stability.

Evaluación comparativa de la concentración de SO₂ vientos arriba de una planta termoeléctrica

Resumen

En la zona urbana de la ciudad de Maracaibo existe una planta termoeléctrica, la cual es una fuente de contaminación atmosférica, produce entre otros poluentes, Dióxido de Azufre (SO₂) formado de la combustión de los componentes fósiles utilizados. Estas emisiones se dispersan en la atmósfera, encontrándose concentraciones en diferentes puntos que inclusive pueden sobrepasar estándares legales de calidad de aire. En este trabajo de investigación se analiza la distribución de las concentraciones del SO₂, hacia la zona norte de la planta termoeléctrica y en la dirección menos preferencial de los vientos. A nivel experimental las muestras fueron recolectadas por absorción a través del sistema RAC 5 modelo all weather shelter, mientras a nivel computacional se aplicaron los modelos de dispersión atmosféricos de Briggs y Pasquill-Gifford. Para la aplicación de los modelos se determinaron datos horarios de las condiciones metereologicas de la ciudad de Maracaibo tales como: velocidad del viento, dirección del viento, nubosidad, insolación y temperatura de los días de muestreo del mes de Diciembre de 2000 y Enero de 2001, al igual que los datos de la planta termoeléctrica tales como: tasa de emisión, altura física de la chimenea, diámetro de la chimenea, y temperatura de salida de los gases. El modelo de Briggs par la dispersión urbana ofreció la mejor aproximación con los valores experimentales medidos en la zona de muestreo teniendo en cuenta condiciones metereologicas locales y datos de la fuente puntual.

Recibido: 04 de Julio de 2002

Revisado: 14 de Julio de 2003

I. Introducción

El diagnostico de la calidad del aire por si solo representa un aspecto fundamental en el proceso de definición, legislación, control y evaluación de la contaminación atmosférica en las diversas regiones. Los factores que participan en la problemática de la contaminación de las atmósferas urbanas son múltiples, entre otros características físico-químicas del efluente gaseoso, configuración de la fuente, condiciones meteorológicas del área, y perfil topográfico del área [1].

La calidad del aire de la atmósfera de la ciudad de Maracaibo esta directamente influenciada por la climatología de la cuenca del lago de Maracaibo, la cual comprende una extensa superficie de agua, rodeada en su mayor parte por cordilleras, con un solo lado abierto hacia el norte; esta topografía particular hace que en la cuenca ocurra fenómenos locales cuya intensidad depende en gran medida de los vientos reinantes sobre el Caribe y Venezuela [2].

La existencia de una planta termoeléctrica en la ciudad de Maracaibo constituye una importante fuente emisora de dióxido de azufre (SO₂), oxidos de nitrógeno (NO_X), y cenizas en forma de partículas suspendidas, que afectan adversamente la calidad del aire de zonas cercanas a ella. La planta termoeléctrica está ubicada a 10°36^’15” de latitud norte y 71°36’50” de longitud oeste en el municipio Maracaibo; la capacidad instalada de producción es de 804 Mw., el combustible de uso frecuentemente es el fuel oil residual que contiene alrededor de 3% de azufre, 450ppm de vanadio, 0.1% de cenizas y un poder calorífico de 18055 BTU/lbm [3]

Araujo y otros [4] realizaron en 1987 un muestreo de SO₂ en áreas alrededor de la planta termoeléctrica con los siguientes valores extremos encontrados en las diferentes estaciones: en la dirección preferencial de los vientos, A (sur, 650 m de la planta) con 3.70- 22.9µg/m³, B (sur, 1200 m de la planta) con 41.60 – 117.6 µg/m³; vientos arriba de la planta, E (norte, a 3700 m), con 7.8 - 23.4µg/m³; hacia el oeste de la planta, C (oeste a 500 m) con 39.8 - 53.8µg/m³, D (oeste a 800 m) con 29.4 – 247.4µg/m^3. Las concentraciones de SO₂ resultaron por encima del standard de calidad del aire y en relación directa con la continua emisión de SO₂ desde la planta termoeléctrica.

Rossell y otros [5] efectuaron en 1989 un monitoreo de SO₂ y partículas suspendidas totales(PST) en áreas aledañas a la planta termoeléctrica y reportaron los siguientes valores extremos en las diferentes estaciones, en la dirección preferencial de los vientos, A (sur a 1200 m) con SO₂ 17.1 – 189.9µg/m³ y PST 53.9 – 192.1µg/m³, B (sur a 1800 m) con 5.7 – 72.8µg/m³ y PST 46.3 – 76.2 µg/m³, vientos arriba de la planta, E (norte, a 3000 m) con SO₂ 17.3 – 23.3µg/m³, y PST 22.5 – 79.7 µg/m³; hacia el oeste de la planta, C (oeste a 800 m) con SO₂ 11.4 – 30.2µg/m³, y PST43.4 - 101.2µg/m³; D (oeste, a 1800 m) con SO₂ 5.6 – 66.5µg/m³, y PST 40.9 – 66.0µg/m³. De todos los sitios de muestreo la estación A resultó la más influenciada por la planta termoeléctrica con los niveles de SO₂ y PST por encima del estándar anual de calidad del aire.

García y otros [6] acoplaron el modelo de dispersión Gaussiana desarrollado por Turner [7] a las condiciones operacionales de la planta termoeléctrica estudiada y a las estabilidades atmosféricas de la cuidad de Maracaibo en la predicción de las concentraciones de SO₂ y su comparación con los valores experimentales encontrados en los alrededores de la planta termoeléctrica. Los resultados indican que para la clase de estabilidad C que predomina en la atmósfera de la ciudad de Maracaibo y con los coeficientes de dispersión dados por Pasquill-Gifford se obtiene una alta correlación entre las concentraciones estimadas por el modelo y las medidas de SO₂ en el área.

II. Fundamentos Teóricos

Para predecir las concentraciones de SO₂ en cualquier posición (X, Y, Z) tomando como origen de las coordenadas cartesianas la planta termoeléctrica Ramón Laguna, siendo X el eje trazado por la dirección preferencial de los vientos, Y el eje perpendicular al X sobre el nivel del suelo, Z el eje Normal al eje X hacia arriba del suelo, se utiliza la ecuación Gaussiana de dispersión desarrollada por Turner [7]:

                                    Q                             - Y²                       -(Z - H_e)²                   -(Z + H_e)²

C_(X,Y,Z,He) = ————————  exp[ ———— ]{ exp[ ————— ] exp[ —————]}     (1)

                      2* π * µ * s_Y * s_Z              2 * s_Y²                       2 * s_Z²                        2 * s_Z²      

donde la altura efectiva involucra la altura física mas la elevación por los efectos combinados de temperaturas y momento, He = H +DH, el programa de computación a partir de los datos operacionales de descarga de las emisiones calcula con las ecuaciones de Briggs[8] la altura adicional DH,

                 F^1/3l ^2/3

DH = 160 ———               (2)

donde el parámetro l depende de la clase de estabilidad atmosférica, si es neutral o inestable.

l = 49F^5/8    si F < 55          (3)

l = 119F^2/5   si F ≥ 55         (4)

Mientras que para la clase de estabilidad atmosférica estable:

                µ

l = 207 ——                       (5)

               S^1/2

donde el parámetro flotante F esta definido en función de las condiciones de descarga desde el tope de la chimenea, así:

          T_s - T_a           D 

             T_a               2

Y el parámetro de estabilidad S definido a parte del gradiente vertical de temperatura, de esta manera:

         g          dT            dT

S = —— [ —— + ( —— ) adiabatica ]        (7)

        T_a         dz             dz

S= 0, atmósfera adiabática o estabilidad neutra; S>0 atmósfera estable; S<0 atmósfera inestable.

El modelo Gaussiano de dispersión esta orientado por los coeficientes de dispersión s_y, s_x como función de la dirección preferencial de los vientos (x) y según la estabilidad atmosférica de Pasquill [9] mostrada en la Tabla 1. En este trabajo de investigación se utilizan tres aproximaciones para los coeficientes de dispersión:

i. Ecuaciones de Briggs [10], para dispersión urbana, Tabla 2.

ii. Ecuaciones de Briggs [10], para dispersión rural, Tabla 2.

Las siguientes ecuaciones (8 y 9) son utilizadas para estimar los parámetros de dispersión vertical y horizontal de Pasquill-Gifford:

                                              T

          [1000* X * Tan(—————)]

                                        57.29569

s_y = ——————————————             (8)

                               215

donde:

s_y: es el parámetro de dispersión horizontal

T: es determinada para cada estabilidad como función de X, Tabla 3

s_z = a X^b        (9)

donde:

s_z: es el parámetro de dispersión vertical.

X: es la distancia a favor del viento en Km.

a y b: es determinado para cada estabilidad como función de X, Tabla 4.

Tabla 1
Clave para la clase de estabilidad Pasquill

Tabla 2
Coeficiente de dispersión rurales y urbanos de Briggs

Tabla 3
Parámetros de Pasquill Gifford

Tabla 4
Parámetros a y b de Pasquill Gifford

III. Metodología y Análisis de Variable

Se seleccionaron tres estaciones de muestreo localizadas hacia la zona norte de la planta termoeléctrica, en la dirección menos preferencial de los vientos, denominadas de la siguiente manera: Estación N° 1 ubicada a una distancia de 3000 m, Estación N° 2 a 2400 m y Estación N° 3 a 700 m, tal cual se muestra en la Figura 1. La recolección de las muestras de aire se realizo con un absorbedor tipo RAC 5 modelo All Weather Shelter, utilizando solución de peroxido de hidrogeno al 0.6% y posterior medición por cromatografía iónica. Las condiciones de operación del RAC 5 se calibraron para un caudal de 1.80 ft³ /min., 29 pulgadas de Hg. de vació y una temperatura de 32.2°C; mientras que para el cromatógrafo iónico se emplearon columnas Dionex Ion Pac-AS4A y como eluyente la mezcla de Na₂CO₃ 1.8 nM - NaHCO₃ 1.7nM, los resultados del monitoreo del SO₂ para las estaciones 1-3 se presentan en las Tablas 5 y 6 respectivamente.

Figura 1. Mapa donde se ubican las estaciones seleccionadas

Tabla 5
Concentraciones de SO₂ en las Estaciones Nº 1 y 2

Tabla 6
Concentraciones de SO₂ (µg/m³)
en la Estación Nº 3

Las variables meteorológicas correspondientes al periodo de muestreo fueron reportadas por la estación metereologica, la Cañada de la Fuerza Aérea de Venezuela. En la Figura 2 se presenta la variación de la velocidad de los vientos en la atmósfera de la ciudad de Maracaibo como promedio horario para los meses de Diciembre y Enero, con valores entre 2 y 5 m/s la velocidad de los vientos va disminuyendo como avanza la mañana entre 1:00 am y 12:00 am, luego del medio día, la velocidad empieza a aumentar hasta alcanzar los valores máximos durante la noche entre 8:00 pm y 10:00 pm. En la Figura 3 se presenta la variación de la temperatura del aire atmosférico como promedio horario para los meses de diciembre y enero con los valores máximos entre 35°C y 36°C hacia el mediodía y transcurso de la tarde entre 12:00m y 4:00 pm; mientras los mínimos entre 23°C y 25°C ocurren entre 4:00 am y 8:00 am.

Figura 2. Variación horaria promedio de la velocidad del viento en la ciudad de Maracaibo.

Figura 3. Variación horaria promedio de la temperatura en la ciudad de Maracaibo.

En relación a la dirección de los vientos, el análisis del promedio horario indica que en general durante las primeras horas del día entre 1:00 am y 9:00 am y las ultimas horas de la tarde y noche entre 6:00 pm y 12:00 pm la dirección predominante de los vientos es nor-este; lo anterior implica que las emisiones de SO₂ por parte de la fuente generadora de energía eléctrica estarían siendo distribuidas hacia lo que se conoce como vientos abajo durante un periodo de tiempo de 14 horas. Por el contrario durante el horario cercano al medio día y primeras horas de la tarde, la dirección sur-este es la prevaleciente; similarmente, ello significa que las emisiones de SO₂ de la fuente bajo estudio, estarían dispersándose en la dirección catalogada como vientos arriba por un periodo de tiempo de 6 horas. En la Tabla 7 se presenta los valores diarios predominantes de la dirección de los vientos y de la estabilidad atmosférica según el criterio de Pasquill (9) y los valores promedios diarios de la velocidad de los vientos y la temperatura atmosférica.

Tabla 7 ]]> Parámetros meteorológicos registrados los días de muestreo

IV. Discusión de Resultados

Del análisis de los resultados del Estación N° 1, los valores de SO₂ se ubican entre 66.20µg/m³, y 85.20µg/m³, con 76.07µg/m³ como concentración promedio para el periodo de muestreo del 07 al 13 de Diciembre de 2000 con velocidades de vientos entre 2.08 m/s y 3.82 m/s y temperaturas entre 27°C y 29°C. En principio las concentraciones de SO₂ registradas en la Estación N° 1 fluctúan en un rango considerado moderado por la uniformidad de las condiciones meteorológicas predominante para el periodo de muestreo y señala al mismo tiempo la posibilidad de la existencia de una línea base de SO₂ en las adyacencias de la planta termoeléctrica.

En la Estación N° 2, las concentraciones de SO₂ resultaron entre 58.90µg/m³ y 75.52µg/m³ con un valor promedio de 66.26µg/m³ para el periodo de muestreo correspondiente del 14 al 20 de Diciembre de 2000, con variaciones meteorológicas de velocidad del viento entre 2.82 m/s y 4.59 m/s y temperaturas entre 27°C y 28°C. Del análisis comparativo de las Tablas 5 y 7 (meteorología del lapso) se observa la tendencia hacia la disminución de la concentración de SO₂ con la disminución de la velocidad de los vientos; 72.52µg/m³, 70.57µg/m³, 64.73µg/m³, mientras la velocidad de los vientos varia así 3.74m/s, 3.13m/s y 2.84m/s, esto se corresponde con las premisas del modelo Gaussiano de dispersión, donde la concentración es inversamente proporcional a la velocidad del viento.

En la Estación N° 3, los valores extremos de la concentración de SO₂ se ubicaron entre 58.89µg/m³ y 65.71 µg/m³ con un valor promedio de 63.97 µg/m³ para el periodo de muestreo del 12 al 18 de Enero de 2001, bajo las velocidades del viento entre 3.09 m/s y 4.79 m/s, y temperaturas entre 27.49°C y 36.01°C.

Del análisis global del muestreo experimental de SO₂ en las Estaciones 1-3, se destaca las siguientes observaciones:

i. Las concentraciones de SO₂ aumentan con el incremento de la posición horizontal( trazada en el eje preferencial de los vientos), esto es los valores de la Estación N° 1 (3000 m) resultaron mayores que los de la Estación N° 2 (2400 m) y estos mayores que los de la Estación N° 3 (700 m).

ii. Los valores promedios de SO₂ en función de la posición horizontal así: 76.07µg/m³ a 3000 m, 66.26µg/m³a 2400 m, 63.97µg/m³ a 700 m desde la base de la Planta termoeléctrica. Esta distribución de concentración de SO₂ con la distancia X se corresponde con lo predicho en el modelo Gaussiano de dispersión.

iii. El rango de variación de la concentración de SO₂ en la medida que se aleja de la Planta termoeléctrica tiende a una mayor abertura entre los niveles máximos y mínimos, así DC= 19.00µg/m³ a 3000 m, DC= 13.62µg/m³a 2400 m, DC= 6.82µg/m³ a 700 m.

iv. La existencia de un rango de variación moderada entre los niveles de SO₂ aconseja sobre la presencia de una línea base de SO₂ en la adyacencia de la Planta Termoeléctrica como consecuencia de la emisión continua de SO₂ por más de 50 años.

1. La tasa de emisión de 535.09 g/s correspondiente a la generación mínima de energía eléctrica en la planta.

2. Los coeficientes variables de dispersión en las direcciones, Y, Z, esto es, los parámetros s_y y s_z hasta encontrar el orden de Magnitud apreciable entre las concentraciones experimentales y calculadas.

Los resultados de esta tercera aproximación se presentan en los Tablas 8, 9 y 10, para el Modelo de Briggs de dispersión urbana; en las Tablas 11, 12 y 13 para el Modelo de Briggs de dispersión rural; y en las Tablas 14, 15 y 16 para el Modelo d Pasquill-Gifford; del análisis global de los anteriores resultados se desprende las siguientes anotaciones:

Tabla 8
Ajuste de las concentraciones calculadas para el modelo Briggs urbano en la Estación Nº 1

Tabla 9
Ajuste de las concentraciones calculadas para el modelo Briggs urbano en la Estación Nº 2

Tabla 10
Ajuste de las concentraciones calculadas para el modelo de Briggs urbano en la Estación Nº 3

Tabla 11
Ajuste de las concentraciones calculadas para el modelo de Briggs rural en la Estación Nº 1

Tabla 13
Ajuste de las concentraciones calculadas para el modelo de Briggs rural en la Estación Nº 3

Tabla 14
Ajuste de las concentraciones calculadas para el modelo de Pasquill-Gifford en la Estación Nº 1

Tabla 15
Ajuste de las concentraciones calculadas para el modelo de Pasquill-Gifford en la Estación Nº 2

Tabla 16
Ajuste de las concentraciones calculadas para el modelo de Pasquill-Gifford en la Estación Nº 3

i. Indudablemente que la distribución Gausiana o distribución normal de la concentración varia con los parámetros de dispersión s_y y s_z que según Pasquill la relación binaria (s_y y s_z ) conjuntamente con la dirección preferencial de los vientos define igualmente la Estabilidad Atmosférica. De tal modo que al variar la relación binaria (s_y y s_z) se estarán analizando las estabilidades A, B, C, D, E, y F, en la determinación de las concentraciones.

ii. El Modelo de Briggs de dispersión urbana ofrece la mejor aproximación para la primera estación para la relación binaria s_y= 373.99 metros y s_z = 354.99 metros equivalente a una estabilidad atmosférica D, con rango de concentración entre 51.84µg/m³ y 62.36µg/m³ referido al rango experimental entre 85.18µg/m³ y 66.19µg/m³. Para la segunda estación de muestreo, el mejor ajuste corresponde al par de valores sy= 335.64 metros _y s_z = 316.59 metros equivalente a una estabilidad atmosférica D, con un rango de concentración entre 38.16µg/m³ y 51.48µg/m³ comparados con el rango experimental entre 58.89µg/m³ y 72.52µg/m³. Y para la tercera estación el mejor ajuste se logra para la relación binaria s_y = 216.68 metros y s_z = 247.81 metros correspondiente a una estabilidad atmosférica A o B con el rango de concentración entre 54.36µg/m³ y 67.68µg/m³ en comparación con el rango experimental entre 58.89µg/m³ y 66. 68µg/m³.

iii. Del análisis anterior se reafirma que para la estación mas cercana a la Planta Termoeléctrica, Estación N° 3 los valores de concentración de SO₂ pueden expresarse por la distribución Gaussiana para el Modelo de Briggs de dispersión urbana con los coeficientes binarios s_y y s_z (216.68 m; 247.81 m) para la estabilidad atmosférica B en el periodo de muestreo del presente estudio, y para la tasa mínima de emisión.

iv. Para las Estaciones más lejanas de la Planta Termoeléctrica, Estación N° 2 y N° 3, las predicciones del modelo de Briggs de Dispersión Urbana resultaron ligeramente menores a las concentraciones experimentales para la tasa mínima de emisión de SO₂. Si el valor real de descarga de SO₂ desde las chimeneas coincidiera con la tasa mínima de emisión, entonces habría que considerar otras fuentes emisoras de SO₂ que contribuyan a mantener los niveles registrados en las Estaciones N° 2 y N° 3. En otras palabras, la Planta Termoeléctrica aporta la concentración de 62.36µg/m³ mientras el diferencial DC = 22.82µg/m³ tendrá que ser aportado por otras fuentes de emisión (Por ejemplo: El flujo vehicular).

v. Las concentraciones calculadas por el modelo de Briggs de Dispersión Rural resultan menores a las reportadas por el modelo de Briggs de Dispersión Urbana en referencia a las mediciones experimentales. Sin embargo los valores calculados se aproximan en orden de magnitud a los experimentales para la estabilidad atmosférica C y los coeficientes binarios s_y= 345.53 metros y s_z= 223.11 metros para la Estación N° 1; para la estabilidad atmosférica B y los coeficientes binarios s_y= 164.25 m y s_z = 155 metros para la Estación N° 3.

vi. Las mejores aproximaciones del modelo de Pasquill-Gifford se resumen así: Estabilidad atmosférica C y coeficientes binarios s_y= 334.78 metros y s_z= 200.97 metros para la estación N° 1, estabilidad atmosférica B y coeficientes binarios s_y= 427.83 metros y sz= 385.54 metros para la estación N° 2; y estabilidad A y coeficientes binarios, s_y= 166.66 metros y s_z= 264.61 metros para la estación N° 3. Los valores de concentración calculados resultaron menores a los valores experimentales y menores a las predicciones del modelo de Briggs de dispersión urbana. En las Tablas 17, 18 y 19, se observan los mejores ajuste de los modelos de Briggs y Pasquill-Gifford para la tasa mínima de emisión de SO₂ para las correspondientes estaciones 1, 2 y 3. Tales Tablas son de particular importancia para adecuar la predicción y adecuación del modelo a las mediciones experimentales de SO₂ en las referidas estaciones. Lo anterior se explica con los siguientes ejemplos:

Tabla 17
Concentraciones Ajustadas Estación N° 1 (Hielo el Toro)

Tabla 18
Concentraciones Ajustadas Estación N° 2 (Extintores del Zulia)

Tabla 19
Concentraciones Ajustadas estación N° 3 (Carrocerías Buena Vista)

1. Si la tasa de emisión de SO₂ corresponde a la generación mínima de energía eléctrica de 100 Mw., esto es, 535.09g de SO₂ /s la Estación N° 3 es entonces representada por el modelo de Briggs para dispersión urbana.

3. Si la tasa de emisión de SO₂ resulta ser el 52 % mayor a la tasa mínima de emisión, esto es, 813.34 g de SO₂/s, entonces las Estaciones N° 1 y N° 2 pueden expresarse por el Modelo de Pasquill-Gifford.

Tabla 20
Concentraciones Calculadas de SO₂ ajustando flujo masico (m) a 730,86g/s

V. Conclusiones

Las principales conclusiones de este estudio enfocado en la medición experimental y de la aplicación de tres modelos de dispersión atmosférica, basándose en el principio Gaussiano para estimar las concentraciones de Dióxido de Azufre emitidas por la planta termoeléctrica Ramón Laguna son las siguientes:

Del análisis global del muestreo experimental del SO₂ se destaca:

1. Las concentraciones de SO₂ aumentan con el incremento de la posición horizontal (trazada en el eje preferencial de los vientos), esto es los valores de la estación N° 1 (3000 metros) resultaron mayores que los de la estación N° 2 (2400 metros) y estos mayores que los de la estación N° 3 (700 metros).

2. De la comparación de las concentraciones obtenidas experimentalmente con los valores estipulados por las normas de calidad del aire, se concluye que en la estación N° 1 tan solo un 30% de las muestras superan el limite de 80µg/m³ cumpliendo así con el porcentaje de excedencia permitido; para el caso de las Estaciones N° 2 y N° 3, en la totalidad de los días muestreados los valores obtenidos resultaron inferiores a la norma.

Del análisis global de las concentraciones calculadas del SO₂ se refiere:

3. Similarmente al comportamiento observado para los valores experimentales, las predicciones de los modelos de Briggs (dispersión urbana y dispersión rural) y Pasquill-Gifford describen adecuadamente el orden de magnitud para las diferentes Estaciones de monitoreo sujeto a la tasa de emisión de SO₂ y de la estabilidad atmosférica, sin embargo, el modelo de Briggs de dispersión urbana ofrece la mayor aproximación respecto a los valores experimentales de SO₂.

Referencias Bibliográfícas

1. García C.: “Evaluación cuantitativa de Polución del Aire”, edición Septiembre 1987, Facultad de Ingeniería, LUZ, paginas 138-160, Maracaibo.        [ Links ]

2. Rossell N.: “Diagnostico de la Química Atmosférica y Procesos de Remoción de los Principales Compuestos de Azufre en Maracaibo”, Trabajo de Ascenso categoría asociado, Facultad de Ingeniería, LUZ, 100 pagina, Octubre 1990.        [ Links ]

3. González I. y Tremont C. “Efectos Socio-Económicos y Ambientales Ocasionados por el Funcionamiento de Planta Termoeléctrica Ramón Laguna, Sector la Arreaga -Maracaibo”, Trabajo Especial de Grado, Facultad de Ciencias Econrómicas y Sociales, Universidad del Zulia. Maracaibo, Venezuela pp 120, 1987.        [ Links ]

4. Araujo J., Petit E. y Flores C.: “Determinación Cuantitativa de los Contaminantes Gaseosos Emitidos por la Planta Termoeléctrica Ramón Laguna”, Trabajo Especial de Grado, Escuela de Ingeniería Civil, Universidad del Zulia. Maracaibo, Venezuela, pp 115, 1987.        [ Links ]

5. Rossell N., Velásquez H. y Linares T.: “Niveles Atmosféricos de SO2, PST y Sulfatos alrededor de una Planta termoeléctrica en Época Seca”, Rev. Tec. Univ. Zulia, Vol. 15 N° 2, 79-85, 1992.        [ Links ]

6. García C., Rojano R. y Rossell N.: “Efectos de las Variables Operacionales y Meteriologicas en la Dispersión Atmoferica del SO₂. Rev. Tec. Univ. Zulia, Vol. 23, No. 2, 150-160, 2000.        [ Links ]

7. Turner B.: “Worbook atmospheric dispersion estimates: an introduction to dispersion modeling”, Lewis Publishers, CRC Press Inc., London, pp 144, 1994.        [ Links ]

8. Briggs G.A.: Some resents analyses of plume rise observations, pp 1029 – 1032 in Proceedings of the second international Clear Air Congress, England H.M. and Berry W.T. Editors, Academic Press, New York.        [ Links ]

9. Pasquill F.: The estimations of the dispersion of windborne material, Meteorology Magazine, 90 (1063), 33 - 49, 1961.        [ Links ]

10. Briggs G.A.: Discussion on chimney plumes in neutral and stable surrounding, Atmos, Enviran. 6, 507 - 510, 1972.        [ Links ]

11. Gifford F.A.: Atmospheric dispersion calculations using the generalized Gaussian plume model, Nuclear safety, 2 (2), 56 – 68, 1960.        [ Links ] ]]> 1 LUZ^dFacultad de Ingeniería Sept ie mb 138-160 2 LUZ^dFacultad de Ingeniería Octu br e 100 3 1987 120 4 1987 115 5 1992 15 2 2 79-85 6 2000 23 2 2 150-160 7 1994 144 8 1029 - 1032 9 1961 90 1063 1063 33 - 49 10 1972 6 507 - 510 11 1960 2 2 2 56 - 68