Вадим Романов - Прикладные аспекты аварийных выбросов в атмосферу

δ 2 z= Кt

где К = К х= К у= K z— эффективное значение коэффициента диффузии, определяемое состоянием атмосферы.

Кроме того, учитывается зависимость дисперсий клубов в виде облаков от скорости ветра с помощью соотношений:

σ 2= σ 2 T+ σ 2 дол

где σ 2 T— дисперсия, определяемая только коэффициентом горизонтальной турбулентной диффузии; σ 2 дол— дисперсия, описывающая рассеяние примеси по горизонтали в потоках с изменяющейся по высоте скоростью за счет вертикальной турбулентности.

При сильном перемешивании примеси по вертикали и большом градиенте скорости ветра вклад слагаемого σ 2 долв суммарную дисперсию может оказаться значительно большим, чем σ 2 T.

Такой метод позволяет учитывать реальную метеорологическую обстановку в месте инцидента, однако его вычислительный аппарат очень громоздок и сложен. Некоторые получаемые результаты, учитывая заложенную в расчет гауссову модель распределения концентраций, носят иллюстративный характер.

В работе [135], проведенной в районе испытаний ядерных энергетических установок, получена полуэмпирическая формула для расчетов концентраций радиоактивной примеси. Она записывается так:

В этой формуле предполагается использование следующих выражений для характеристик дисперсий загрязнений по координатным осям:

где U — скорость ветра в слое распространения струи ядерной энергетической установки; [U]=км/час; t — время процесса диффузии, час; г — расстояние выброса от места проведения работ, км; U x,U y— составляющие скорости ветра по осям х и у.

Недостатком этой формулы является узкий диапазон ее возможного применения — в месте испытаний. Кроме того, она не позволяет провести анализ влияния исходных параметров установки и среды на результаты расчетов.

Более универсальный характер имеет формула для расчета приземной концентрации на оси прохождения радиоактивного облака. Она получена путем обобщения экспериментальных после подземных ядерных взрывов [154]. В случае мгновенного выброса примеси получают:

где Q — общее количество выброшенной из источника примеси; U — как и ранее, скорость ветра в слое распространения загрязнений, считающаяся постоянной.

Дисперсии примеси описываются следующими временными зависимостями:

σ 2~ t 2при t малых;

σ 2~ 2Kt при t больших.

Распределение примеси вдоль направления распространения облака можно получить, интегрируя записанное выше уравнение по времени после подстановки в него значений метеопараметров и диффузии.

Приведем широко используемую для инженерных оценок эмпирическую формулу для расчета концентраций загрязняющих веществ при выбросах примеси из мощных источников типа дымовых труб тепловых электростанций или химических предприятий. Она имеет следующий вид:

где α и β — некоторые постоянные; U — скорость ветра на высоте флюгера; Q — мощность выброса.

Высота источника Н, входящая в эту формулу, складывается из высоты трубы и начального (динамичного) подъема струи ΔН:

где W 0,R 0, и ΔТ 0— начальные значения скорости газа струи, ее радиуса и перегрева; g — ускорение силы тяжести; θ — температура окружающего воздуха в абсолютной шкале.

Недостатком приведенных выше формул является отсутствие универсальности в выборе коэффициентов α и β, а также некорректность при U → 0 С уменьшением скорости ветра до нуля динамический подъем струи и концентрация загрязнений неограниченно возрастают. Вместе с тем известно, что при инверсионных состояниях атмосферы эти условия заведомо не выполняются, так как существует некоторый «потолок» для начального подъема примеси.

В заключение этого раздела приведем формулы для оценок влияния параметров диффузии на максимальную концентрацию примесей С mи расстояние х mот источника до этого максимума [150]. Зависящий от устойчивости атмосферы режим распространения описывается сигма — значениями σ уи σ z, входящими в расчетную формулу гауссовой модели дымового факела: