|
73 Это облако в модели имеет цилиндрическую форму (высота Н, радиус R) (рис. 1.8). Под действием сил тяжести происходит увеличение горизонтальных размеров облака, а сверху часть газа уносится ветром. Скорость расширения облака моделируется зависимостью U где р плотность газа в облаке; ра плотность воздуха; g ускорение свободного падения; at постоянная воздуха. Скорость уноса массы из облака (интенсивность вторичного источника) £ зависит от скорости ветра, плотности газа в облаке и его размеров и определяется соотношением где £0 скорость ветра на высоте Zo; Мо масса газа в облаке; я, постоянная воздуха. Изменение массы облака по времени описывается формулой где Es -■ интенсивность первичного источника. Масса газа в облаке к моменту действия первичного источника Г/ определяется по формуле 1-ехр / Изменение площади основания облака по времени описывается формулой / . . \д1/2 1 1 Учитывая вышеизложенное, для потенциально опасных технологических установок с учетом вероятного сценария аварии (рис. 1.6) проведены вычислительные |
106 * (для Московской области Л = 140); коэффициент 77 (77= 1); средняя температура окружающего воздуха за самый жаркий месяц по данным Гидрометцентра о 25 С. Кроме того, учитываются параметры, характеризующие точечные, линейные и плоские источники загазованности промышленных территорий. Основным модулем программы является программа расчета параметров поля концентраций, создаваемого одним точечным источником поллютанта. На рис. 1.32 в качестве примера приведено распределение концентраций метана, образованных в результате разгерметизации адсорбера К-401 на установке получения водорода. Рассеивание большого количества тяжелого отличается дисперсии малых примесей. Поэтому для оценки рассеивания примеси важно знать поведение облака, образующегося вблизи первичного источника при мощном выбросе. Это облако в модели имеет цилиндрическую форму (высота Н, радиус К) (рис. 1.33). Под действием сил тяжести происходит увеличение горизонтальных размеров облака, а сверху часть газа уносится ветром. Скорость расширения облака моделируется зависимостью: Uf=atA\g (Р-Ра) Ра где р плотность газа в облаке; ра плотность воздуха; g ускорение свободного падения; at постоянная воздуха. Скорость уноса массы из облака (интенсивность вторичного источника) Е зависит от скорости ветра, плотности газа в облаке и его размеров и определяется соотношением: р 7 2 Ра U О лЕ р = а t -7 Р Zo м где Uo скорость ветра на высоте Zo; Мо масса газа в облаке; at постоянная воздуха. 108 г Изменение массы облака по времени описывается формулой dt Р Z о где Es интенсивность первичного источника. Масса газа в облаке к моменту действия первичного источника ti определяется по формуле: Z М 1 EsZop аРоРа \-exp atPaU pz0 \ о t1 / Изменение площади основания облака по времени описывается как dS_ dt ng, Ра 1/2 2 аt г 1 1 РУ Учитывая вышеизложенное, для потенциально опасных технологических установок с учетом вероятного сценария аварии (рис. 1.31) проведены вычислительные эксперименты рассеивания выбросов тяжелых газов и паров разлившейся жидкости. Для каждого объекта рассчитывались зависимости осевой концентрации примеси С на уровне земли, параметров вертикальной и горизонтальной дисперсии Sy полуширины зоны постоянной концентрации В (м), эффективной полуширины Вер (м) и высоты Нер(м) облака, плотности облака и времени прихода примеси г (с) от осевой координаты X (м) по направлению ветра. Кроме этого, определялись координаты линии постоянных концентрации на различных уровнях по высоте Z В качестве граничных выбирались концентрации, соответствующие 5, 25, 50 и 100 % от уровня НКПР. Для технологических установок в табл. 1.36 приведены исходные данные для расчета параметров рассеивания выбросов углеводородных газов, результаты которого даны в табл. 1.37. |