ПРОГНОЗИРОВАНИЕ ПОСЛЕДСТВИЙ АВАРИЙНЫХ ЗАЛПОВЫХ ВЫБРОСОВ СЖИЖЕННЫХ ГАЗОВ (26.09.2012)

Автор: Старовойтова Евгения Валерьевна

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

в написании статей и тезисов, участии в конференциях.

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

Практическая значимость работы состоит в том, что предлагаемый комплекс моделей может использоваться для решения задач оценки последствий аварий на объектах хранения, переработки и использования сжиженных газов при разработке Планов локализации и ликвидации аварийных ситуаций (ПЛАС), разделов «Инженерно-технические мероприятия Гражданской обороны. Мероприятия по предупреждению чрезвычайных ситуаций» (ИТМ ГО и ЧС) и «Мероприятий по обеспечению пожарной безопасности» проектной документации, Деклараций промышленной безопасности опасных производственных объектов, Деклараций пожарной безопасности, паспортов безопасности опасных объектов, при выборе пассивных мер защиты по ограничению распространения опасных веществ на стадии проектирования.

Методика и программная система использовались при разработке планов локализации и ликвидации аварийных ситуаций для ОАО «Татнефтегазпереработка» и аммиачно-холодильных станций ОАО «Нижнекамскнефтехим», при разработке паспорта безопасности и декларации промышленной безопасности для хлораторных МУП «Водоканал», паспорта безопасности для ОАО «Татнефтегазпереработка».

Апробация работы. Основные положения и результаты диссертационной работы докладывались и обсуждались на Всероссийской студенческой научно-технической конференции «Интенсификация тепло-массообменных процессов, промышленная безопасность и экология» (Казань, 2005, 2008, 2012); Всероссийской научно-практической конференции «Инновации и высокие технологии XXI века» (Нижнекамск, 2009); Международной научно-практической конференции «Устойчивое развитие территорий: управление природными, техногенными, пожарными, биолого-социальными и экологическими рисками: материалы» (Оренбург, 2011); XXIV Международной научной конференции «Математические методы в технике и технологиях –24» (Саратов, 2011); Всероссийской молодежной конференции с элементами научной школы «Нефть и нефтехимия» (Казань, 2011); Международной научно-практической конференции «XL Неделя науки СПбГПУ» (Санкт-Петербург, 2011); на ежегодных научных сессиях КГТУ (г. Казань).

Публикации. Основное содержание диссертации отражено в 19-и публикациях, в том числе в 10-и статьях, опубликованных в научных изданиях, входящих в перечень ВАК Минобразования и науки РФ для соискателей ученых степеней доктора и кандидата наук, а также в материалах 9-ти российских и Международных научных конференций.

Структура и объем работы. Диссертация состоит из введения, пяти глав, выводов, библиографического списка, включающего 143 наименования. Общий объем работы составляет 195 страниц, включая 11 таблиц, 68 рисунков.

Основное содержание работы

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

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

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

Рисунок 1 – Схема установки измерения интенсивности парообразования

1– электронные весы;

2 – подложка; 3 – чаша;

4 – сжиженный газ; 5 – песок;

6 – слой фольги; 7 – секундомер

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

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

где CP,liq – удельная теплоемкость жидкости, Дж/(кг(К); T0 – температура сжиженного газа в емкости, К; Tb – температура кипения жидкости, К; (Hg – удельная теплота парообразования при температуре кипения и атмосферном давлении, Дж/кг.

Принималось, что в начальный момент времени (на месте выброса) первичное облако опасного вещества имеет форму цилиндра с радиусом R, равным высоте H:

где mliq — масса жидкости в оборудовании, кг; (выб — плотность опасного вещества в первичном облаке в начальный момент, кг/м3.

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

Для определения движения аэрозольных частиц использовалось уравнение:

где up – проекция скорости движения капли на ось x, м/с; ul – проекция вектора скорости несущей фазы на ось x, м/с; g1 – проекция вектора ускорения свободного падения на ось x, м/с2; ?p – плотность капли, кг/м3; ? – плотность несущей фазы, кг/м3; dt – шаг по времени, с; FD(ul–up) – удельная сила сопротивления на движению капли, Н/кг.

Уравнение, определяющее интенсивность испарения аэрозольных частиц:

где Ap – площадь поверхности капли, м2; Mg – молекулярная масса, кг/моль; Cg,s – концентрация компонента у поверхности капли, моль/м3; Cg,? – концентрация компонента в окружающем газе, моль/м3; ? – коэффициент массоотдачи, м/с, рассчитывался из соотношения Ранца-Маршалла.

Изменение температуры капли определялось уравнением:

где mp – масса капли, кг; CP,p – теплоемкость капли, Дж/(кг(К); (s – коэффициент конвективной теплопередачи, Вт/(м2(К); T( – локальная температура несущей фазы, К.

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

Изменение импульса несущего потока F, Н/м3, обусловленное гидродинамическим сопротивлением аэрозольных частиц, рассчитывалось из следующего соотношения:

– текущий массовый расход частиц данной фракции через грань контрольного объема, кг/с. СD – коэффициент сопротивления для частиц сферической формы; dp – диаметр капли, м; Red – число Рейнольдса.

Источник массы S, кг/(м3·с) в уравнении неразрывности и уравнении переноса примеси, возникающий в результате испарения капель, определялся через соотношение:

Сток энергии в газовой фазе Q, Дж/(м3(с), обусловленной теплообменом с каплями, рассчитывался следующим образом:

 – начальный массовый расход частиц данной фракции через грань контрольного объема, кг/с.

Распределение капель по размерам в источнике задавалось с использованием соотношения Розин-Раммлера.

Модель кипения (испарения) пролива. Оставшаяся после мгновенного вскипания жидкая фаза разливается на подстилающей поверхности и переходит в состояние кипения. Интенсивность парообразования при этом пропорциональна скорости теплопритока из окружающей среды. Схема теплообмена пролива приведена на рис. 2.

Рисунок 2 – Схема теплообмена пролива сжиженного газа

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

При разработке модели парообразования из пролива сжиженного газа были сделаны следующие допущения:


загрузка...