student
Russian Federation
The article considers the modeling of the dispersion of a pollutant during toxicant emissions, depending on the characteristics of atmospheric transport, wind speed and the stability class of the atmosphere. To describe the dependence of the implementation of atmospheric stability clas-ses on wind speed, data corresponding to a certain type of Earth's surface (flat, hilly, etc.) were also analyzed
modeling, pollutants, the stability class of the atmosphere, the emission power of pollutants
Размеры зоны токсической опасности при выбросах токсиканта зависят как от мощности выброса, так и от характеристик атмосферного переноса, прежде всего от скорости ветра и от категории (класса) устойчивости (стабильности) атмосферы [1,2]. Категории различаются в основном интенсивностью вертикального перемещения. Наиболее неустойчивая категория «А» отмечается при слабом ветре и сильной солнечной радиации, когда воздух, нагретый теплом от земной поверхности, всплывает. Обычно это состояние возникает после полудня или несколько раньше. Категория «С» наблюдается при усилении ветра от умеренного до сильного и чаще всего вечером при ясном небе или днём при низких кучевых облаках, а также летними ясными днями при высоте солнца 15 – 30°. Нейтральная категория «D» соответствует условиям сплошной облачности, как днём, так и ночью, когда влияние прямых солнечных лучей незначительно. Устойчивые категории «Е» и «F» фиксируют обычно ночью при чистом небе или слабой облачности, когда земная поверхность выхолаживается и над ней устанавливается инверсионный слой. Естественная конвекция при этом подавляется.
Стабильность повышается с ростом скорости ветра и снижением интенсивности солнечного облучения. Облачность проявляется по-разному.
Ночью её снижение усиливает охлаждение Земли, образование температурной инверсии. Днём, наоборот, конвективные потоки и нестабильность усиливаются.
Для определения категории стабильности используют классификацию по вертикальному градиенту температур (классификация Паскуилла, таблица 1). Поскольку состояние устойчивости атмосферы по существу определяется интенсивностью вертикальных конвективных течений, оно может существенно изменятся и в течение суток. Характерное распределение стабильности атмосферы в течение дня средних широт представлено для различных времён года.
Таблица 1 – Классификация классов устойчивости по Паскуилу
Градиент температуры (ΔT/ΔZ), оС/100 м |
Класс устойчивости |
< 1,9 |
А - сильная конвекция |
-(1,9÷1,7) |
В - конвекция |
-(1,7÷1,5) |
С – умеренная конвекция |
-(1,5÷0,5) |
D – нейтральная |
-0,5÷+1,5 |
Е – инверсия |
Как видно, распределение состояний устойчивости атмосферы в холодное и тёплое время года весьма отличаются. Так, в холодное время года доминируют нейтральное и стабильное состояние атмосферы, а в тёплое время года: ночью - стабильное, а днём - неустойчивое. Для мощных вертикальных высокоскоростных выбросов или высокотемпературных выбросов, а также в случае расположения источника выброса на значительной высоте от поверхности земли весьма важным является учёт зависимости слоя перемешивания (приземный слой атмосферы с примерно постоянным касательным напряжением сдвига) от состояния атмосферы. По данным зарубежных исследований высота слоя перемешивания в среднем изменяется от 100 м ночью до 2000 м в дневное время. При этом максимальное значение высоты слоя перемешивания достигается через 3-4 часа после восхода солнца.
Для описания зависимости реализации того или иного класса устойчивости атмосферы от скорости ветра «U» были проанализированы данные соответствующей определенному типу земной поверхности (равнинная, холмистая и т.д.).
Влияние состояния земной поверхности проявляется также и в эффекте торможения движения воздушных масс, который характеризуется, как правило, так называемым «параметром шероховатости» (Д0) в зависимости от структуры ландшафта и сезона в разрезе года. Наиболее характерные значения Д0 даны в таблице 2.
Известно, что скорость ветра существенно изменяется с высотой [1]. В прикладных исследованиях наиболее часто используется степенная зависимость вида:
Значения показателя «р» также зависят от класса устойчивости атмосферы и "шероховатости" поверхности Д0 (таблица 2).
Таблица 2 – Зависимость параметра "р" от величины шероховатости
Категория стабильности атмосферы |
Параметры шероховатости Д0, м |
|||
0,01 |
0,1 |
1 |
3 |
|
А |
0,05 |
0,08 |
0,17 |
0,27 |
В |
0,06 |
0,09 |
0,17 |
0,28 |
С |
0,06 |
0,11 |
0,2 |
0,31 |
D |
0,12 |
0,16 |
0,27 |
0,37 |
Е |
0,34 |
0,32 |
0,38 |
0,47 |
F |
0,53 |
0,54 |
0,61 |
0,69 |
Также существует связь между категориями устойчивости атмосферы и характеристиками турбулентного переноса с масштабами переноса примеси.
Анализ для простейшего случая переноса «нейтральной» примеси от точечного источника постоянной мощности G0 показал, что распределение концентрации примеси на оси следа (у = 0) на поверхности земли (z = 0) равно:
Обычно коэффициенты дисперсии в горизонтальном и вертикальном направлении и
вычисляют по эмпирическим соотношениям. Наиболее известны номограммы «Гиффорда-Паскуилла», составленные по наблюдениям концентрации на равнинной местности:
,
где x – расстояние, м (102 ≤ x ≤ 104); i = y, z.
Коэффициенты аппроксимации ai, bi и сi, даны в таблице 3:
Таблица 3 – Константы формул, аппроксимирующих кривые Гиффорда-Паскуилла
Константы |
Категории стабильности атмосферы |
|||||
А |
В |
С |
D |
Е |
F |
|
ay |
-1,104 |
-1,634 |
-2,054 |
-2,555 |
-2,754 |
-3,143 |
by |
0,9878 |
1,035 |
1,0231 |
1,0423 |
1,0106 |
1,0148 |
cy |
0,0076 |
0,0096 |
0,0076 |
-0,0087 |
0,0064 |
-0,007 |
az |
4,679 |
-1,999 |
-2,341 |
-3,186 |
-3,783 |
-4,49 |
bz |
1,7172 |
0,8752 |
0,9477 |
1,1737 |
1,301 |
1,4024 |
cz |
0,277 |
0,0136 |
-0,002 |
-0,0311 |
-0,045 |
-0,054 |
Бриггс провёл аналогичные наблюдения в городской местности и поэтому его коэффициенты иногда называют «городскими». Для нестабильной атмосферы городские коэффициенты выше сельских приблизительно до 5 км, потом существенно снижаются. С учётом стабильности (для категории F) позиция превышения городских коэффициентов увеличивается до 40 км. Расхождение коэффициентов (
ещё более существенны.
Представленные выше данные об изменениях и корреляциях между основными параметрами атмосферного переноса свидетельствуют о необходимости построения предельно чётких логических схем различных вариантов (исходов) развития аварийного процесса в атмосфере, построенных по принципу «деревьев событий».
Очевидно также, что по целому ряду характеристик (шероховатость дневной поверхности, функции распределения скорости ветра и классов стабильности атмосферы в разрезе года и др.) «функции переноса» непосредственно связаны и с «функцией источника» (давление насыщенных паров, температура воздуха и поверхности грунта и др.), то есть с мощностью выброса.
Проиллюстрируем важность обоснования характерного сценария на примере аварийного струевого выброса токсического газа с постоянной мощностью. Возможная неопределённость – направление выброса. Варианты: 1 – вертикальная свободная струя в сносящем потоке ветра; 2 – «настильная» струя (то есть ориентированная вдоль земной поверхности) с характерной высотой источника - 1,0 м. Для корректности сравнения диаметр источника и масса выброса в обоих случаях приняты одинаковыми.
С использованием классических законов сохранения количества движения, массы и неразрывности можно получить для вертикального выброса следующие выражения для расчёта концентрации токсиканта в расширяющейся струе на высоте «z» от среза трубы радиусом «r0»:
а также «предельную» высоту подъёма струи - Z*, на которой её средняя по сечению скорость будет близка к скорости ветра на этой же высоте - :
где С0, р0 – концентрация токсиканта и плотность газа (смеси) на срезе; ра – плотность атмосферного воздуха; а – коэффициент «захвата» (эжекции) воздуха (в соответствии с многочисленными экспериментами для свободной турбулентной струи, а = 0,06 ÷ 0,08); w0 – скорость выброса.
На высоте Z = Z* струя практически погасит свою исходную кинетическую энергию, и дальнейший процесс рассеивания примеси будет подчиняться законам турбулентного обмена в атмосфере, то есть произойдёт смена инжекционного механизма на диффузионный.
Определив значения Z* = Нэк и C(Z*), где по смыслу Нэк является высотой эквивалентного источника с помощью токсиканта Q* можно перейти ко второй стадии анализа - к расчёту рассеивания по направлению ветра токсиканта в атмосфере под действием массопереноса и турбулентной диффузии.
Для прогноза распределения концентраций токсиканта вокруг источника используют простейшую модель Гаусса турбулентной диффузии.
Математическое выражение для концентрации вещества от точечного источника с постоянной мощностью – Q* записывается в следующем виде:
где Q* – мощность источника, кг/с; и
, – дисперсионные параметры, зависящие от устойчивости атмосферы и расстояния от источника «х», м; U – скорость ветра, м/с; Н – высота источника, м; x, y, z – осевая, поперечная и вертикальная координаты; f(A) – доля примеси в слое перемешивания (А – высота слоя перемешивания).
Эта зависимость применительно к реальным условиям выброса с концентрацией С0 корректируются введением понятия виртуального источника. В этом случае вычисляют х0 таким образом, чтобы в точке (х = 0, у = 0, z = Н) соблюдалось равенство:
А при дальнейших расчётах параметры дисперсии корректируются с учётом значения х0, то есть ;
. Как правило, для относительно небольших значений высот подъёма факела (Н <100 ÷ 200 м) f(A) = l.
1. Vorob'ev E.I., Prusakov V.M., Dushutin K.K. Ohrana atmosfery i neftehimiya. L.: Gidrometeoizdat, 1985. - 231 s.
2. Zamyshlyaev B.V. Vliyanie nachal'nyh razmerov vybrosa (istecheniya) neve-somoy primesi na ocenku koncentracionnyh poley, voznikayuschih pri rasprostranenii oblaka gaussova tipa // Problemy bezopasnosti i ChS. 2004. №4. S. 35-42.