Расчет концентрации поллютанта с помощью барометрической модели с учетом влияния барьера растительности

Введение . Из источника выброса загрязняющее вещество попадает в сложно организованное пространство наземно-воздушной среды, компоненты которого обладают барьерными свойствами [3]. К числу важнейших барьеров относится растительный покров, на внешней границе которого в системе «листовая пластинка – воздух» идут процессы сорбции поллютантов [1, 9–11], заканчивающиеся их поступлением внутрь листа и аккумуляцией в тканях растений. Однако барьерная роль растительности как фактора, изменяющего содержание загрязняющего вещества, не нашла должного отражения в моделях рассеяния и оседания поллютантов [4, 5]. Цель работы – изучить характер изменения содержания диоксида серы в замкнутой среде при наличии на пути диффузии SO 2 растительного барьера и использовать полученные количественные оценки в барометрической модели расчета его концентрации.

Объекты и методы исследований . Для изучения изменений содержания диоксида серы была выполнена серия опытов в стеклянных эксикаторах объемом около 7 л, герметично закрывающихся притертой крышкой с отверстием, через которое по окончании времени экспозиции шприцем брали пробу воздуха объемом 100 мл для анализа. Определение концентрации диоксида серы проводили по стандартной методике, измеряя оптическую плотность на ФЭК-56 при длине волны 400 нм [7]. Источником SO 2 служила химическая реакция сульфита натрия с серной кислотой [8]. Масса навески соли в опытах составляла 10 и 30 мг. На пути диффузии диоксида серы размещали побеги растений, которые крепили специальным образом внутри эксикатора. Срезанные побеги выдерживали не менее 1 ч, чтобы уменьшить влияние газообмена листьев на процессы сорбции поллютанта. Площадь листьев определяли весовым методом, хвои – через удельную линейную плотность [8]. Использовали облиственные побеги основных средообразующих древесных пород селитебной и санитарно-защитной зоны промышленных предприятий г. Ангарска (сосна, береза, тополь).

Результаты и обсуждение . На рисунке 1 представлена схема проведения и наиболее типичный результат опытов с 10-минутной экспозицией. Концентрации SO 2 в воздухе эксикатора с растительным барьером на пути его диффузии (рис. 1, б, г) были достоверно ниже, чем в воздухе эксикатора без облиственных побегов (рис. 1, а, в). Подобное изменение концентрации диоксида серы было обусловлено барьерными свойствами внутреннего пространства эксикатора и находящегося в нем побега. Отношение количества поллютанта в воздухе (экспериментальная концентрация, ЭК) к общему его количеству, выделившемуся в результате реакции (теоретическая концентрация, ТК), характеризовали коэффициентом К h . Разность «1 – К h »

соответствовала доле SO 2 , которая была адсорбирована на поверхности его барьеров ( К адс ). Влияние адсорбции поллютанта на изменение его концентрации ( ΔС ) оценивали формулой

ΔС = С i – C j = С i – К адс · С i = С i · (1 – К адс ) = С i ·К h ,

где С i , C j – концентрации поллютанта в воздухе с учетом и без учета растительного барьера соответственно.

По результатам опытов была рассчитана средняя величина коэффициента адсорбции ( К адс ). В данной статье при апробации модели расчета концентрации поллютанта использовали величины коэффициента адсорбции, полученные в опытах эксперимента, представленного на рисунке 1, б, г.

Рис. 1. Концентрация диоксида серы в воздухе эксикатора при отсутствии и наличии растительного барьера: ТК, ЭК – теоретическая и экспериментальная концентрация поллютанта. Теоретическая концентрация (ТК) рассчитана по уравнению химической реакции образования диоксида серы, экспериментальная концентрация (ЭК) определена по методике [7];

«Навеска 10 мг», «Навеска 30 мг» – навески сульфита натрия, которые были использованы для получения диоксида серы; РП – растительный побег

Повторность каждого из опытов «а – г» в одном эксперименте составляла 3. На рисунке приведены средние величины концентрации диоксида серы, ошибка среднего арифметического не превышала 10%. По схеме, представленной на рисунке, было проведено более 50 экспериментов.

При моделировании изменений концентрации диоксида серы с учетом растительной поверхности на пути его диффузии использовали схему, показанную на рисунке 2. При этом было принято, что массопере-нос поллютанта осуществлялся по оси «Z». Зависимость концентрации газа (n z ) от высоты (Z) определялась по барометрической интегральной формуле для неизотермических условий [5] с учетом местных особенностей температурной стратификации атмосферы, обусловленной теплоотдачей промышленной зоны урбанизированной территории [2]:

In n z = In n о - (m^g/k) j dz / T ,                                     (2)

o где   nо – концентрация молекул поллютанта при Zо;

g – гравитационная постоянная;

k – константа Больцмана;

m – масса частиц газа.

Рис. 2. Схема изменения концентрации поллютанта, попавшего с выбросом в атмосферу: C i – концентрации загрязняющих веществ; H – высота источника; h – высота пространства, занятого барьером растительности; h 1 – высота приземного слоя; h 2 = h - h 1 – высота слоя адсорбции поллютантов (кроны деревьев); сплошные стрелки показывают направление изменения концентрации поллютанта, а пунктирные – «слой адсорбции» поллютанта

При рассеивании в атмосфере и оседании в направлении земли на поверхностях пространства поллютант соприкасается со слоем растительности h 2 (см. рис. 2). При движении поллютанта через этот слой часть его сорбируется листовой поверхностью, это изменяет его концентрацию. Принимали, что до отметки «h» (см. рис. 2) концентрация поллютанта изменялась по барометрической зависимости, но на отрезке высоты «h 2 » помимо барометрического фактора вступил в действие «барьерный» фактор, вносящий изменения в концентрацию поллютанта в соответствии с выражением, представленным формулой (1).

Для расчета концентрации диоксида серы по барометрической модели использовали блок-схему, представленную на рисунке 3. Программа написана в среде Borland Delphi. При расчете концентрации диоксида серы по барометрической модели местоположение барьера растительности отнесли к отметкам высоты 10–20 м. Результаты расчетов изменения концентрации SO 2 с высотой с помощью барометрической модели даны в таблице.

Результаты расчета концентрации диоксида серы с помощью неизотермической барометрической модели без учета и с учетом адсорбции растительным барьером

N	Отметки высоты, м	Температура воздуха, K	Концентрация, кг/м3		Концентрация, кг/м3
			без учета К адс	с учетом К адс	без учета К адс	с учетом К адс
			Вариант 1		Вариант 2
1	1,1000000000E+02	-1,2088575000E+03	5,5870463813E-04	5,5870463813E-04	1,6162527032E-03	1,6162527032E-03
2	1,0000000000E+02	-2,7107150000E+03	5,5867706291E-04	5,5867706291E-04	1,6161729320E-03	1,6161729320E-03
3	9,0000000000E+01	-4,2125725000E+03	5,5864648684E-04	5,5864648684E-04	1,6160844798E-03	1,6160844798E-03
4	8,0000000000E+01	-5,7144300000E+03	5,5861217574E-04	5,5861217574E-04	1,6159852227E-03	1,6159852227E-03
5	7,0000000000E+01	-7,2162875000E+03	5,5857308784E-04	5,5857308784E-04	1,6158721470E-03	1,6158721470E-03
6	6,0000000000E+01	-8,7181450000E+03	5,5852767361E-04	5,5852767361E-04	1,6157407701E-03	1,6157407701E-03
7	5,0000000000E+01	-1,0220002500E+04	5,5847348036E-04	5,5847348036E-04	1,6155839968E-03	1,6155839968E-03
8	4,0000000000E+01	-1,1721860000E+04	5,5840627667E-04	5,5840627667E-04	1,6153895861E-03	1,6153895861E-03
9	3,0000000000E+01	-1,3223717500E+04	5,5831776476E-04	5,5831776476E-04	1,6151335338E-03	1,6151335338E-03
10	2,0000000000E+01	-1,4725575000T+04	5,5818778143E-04	3,2039978654E-04	1,6147575106E-03	1,3321749462E-03
11	1,0000000000E+01	-1,6227432500E+04	5,5793861529E-04	3,2025676517E-04	1,6140367085E-03	1,3315802845E-03