В качестве гипотетического запроектного пожара рассмотривается аварию в хранилище отравляющих веществ. На основании анализа данных по техническому состоянию многих складских помещений и ёмкостей, в которых хранится отравляющий запас России, при любом из указанных выше инициирующих событий наиболее вероятно катастрофическое развитие аварии. Такой ход событий объясняется тем,что сами здания и ёмкости с токсикантами находятся в неудовлетворительном состоянии. Поэтому практически любое из описанных выше внешних воздействий приведёт к максимально возможной аварии. При такой аварии в пожаре будет задействовано всё токсическое содержимое помещений.
Наиболее уязвимым в пожарном отношении звеном в складском помещении, безусловно, является деревянная крыша, загорание которой, как указывалось выше, может инициировать развитие инцидента. Предполагается, что при пожаре ёмкости с отравляющим веществом будут последовательно терять герметичность и разрушаться после обрушения частей горящей кровли. Разлившийся токсикант будет испаряться и воспламеняться, поддерживая пожар примерно на одном уровне. Таким образом, предполагается, что пожар имеет каскадный самоподдерживающийся характер, пока не выгорит всё рабочее тело.
Кроме того, считаем, что при нагревании разрушение больших ёмкостей и цистерн происходит по наиболее слабым их частям — люкам и верхним полуцилиндрам. Через эти разрывы происходит интенсивный выброс паров отравляющих веществ и их загорание под воздействием открытого пламени. Кирпичные стены здания относительно локализуют загорание, не давая ему распространиться в горизонтальном направлении. Высокотемпературные продукты горения формируются при выходе из верхнего проёма хранилища в виде мощного восходящего потока.
При запроектной аварии на арсенале создаются физические условия, при которых горящие и частично разложившиеся отравляющие вещества восходящим потоком выносятся из помещения и доставляются в виде конвективной струи на большие высоты. По мере подъёма струи внутри пограничного слоя атмосферы она многократно разбавляется воздухом и сносится ветром в направлении его движения с учётом высотного разворота. Скорость потока уменьшается^ начиная с некоторой продольной координаты, струя динамически становится неразличимой на фоне пульсационного движения окружающего воздуха.
В этом случае возникает приподнятый объёмный источник токсической опасности. Рассеивание токсикантов из него происходит под действием атмосферной диффузии. Для вычисления полей рассеивания токсикантов из высотного источника необходимо знание его геометрических,массовых, динамических и концентрационных характеристик. Сложностью расчёта высокотемпературного потока при мощном пожаре является то, что он может подняться на высоты большие, чем пограничный слой атмосферы, где характеристики турбулентного обмена и диффузии претерпевают резкие изменения.
Особенность высокотемпературного восходящего потока смеси продуктов горения и разложения отравляющих веществ с воздухом состоит в том, что в формирующийся при пожаре крупномасштабной струе будет присутствовать заметное количество твёрдой и жидкой фаз, пренебречь которыми было бы не корректно. Относительно крупные частицы и аэрозоли появляются при пожаре из-за неполного сгорания отравляющих веществ, элементов конструкции здания и оборудования склада. Они присутствуют в виде сажи, капельной влаги, а также аэрозольных окислов органических и неорганических соединений.
Кроме того, мощные восходящие потоки вовлекаемого в пожар воздуха захватывают атмосферную и поверхностную пыль широкого диапазона дисперсности. В процессе движения внутри струйного потока частицы адсорбируют на свою поверхность газообразные вещества, становясь тоже токсичными.
После подъёма струи на высоты в несколько километров её вещество охлаждается, и транспортируемые потоком частицы под действием силы тяжести выпадают на земную поверхность на значительных расстояниях от места пожара. Таким образом, на удалениях в десятки, а при определённых условиях и в сотни километров от места инцидента, совершенно неожиданно для местного населения возникают токсичные дожди, туманы, смоги или пылевые осадки, приводящие к гибели всё живое.
Количество токсичных веществ,поднимающихся в единицу времени в атмосферу на высоту эффективного источника Но, зависит как от динамики температурного поля очага пожара, так и от метеопараметров в высотном слое 0 < z < Но. Для получения физических характеристик источника загрязнений проводят аэродинамический расчёт струи с учётом жидкой и твёрдой фаз в ней. Умение расчитывать распространение загрязнений в приземном слое атмосферы и в тропосфере совершенно необходимо не только для прогноза опасности персонала, обслуживающего указанные объекты и находящегося в непосредственном контакте с загрязнённой средой, но и для предсказания опасности региона в масштабах ~ 100 км.
Проблемами возникновения и развития больших открытых пожаров, при которых в атмосфере появляются крупномасштабные конвективные движения, переносящие газоаэрозольные продукты горения и дымления на значительные расстояния, занимаются во всём мире. Такие атмосферные образования называют конвективными колонками; они являются разновидностью струйнх потоков в свободной атмосфере [73 — 81]. Отмечается,что конвективные колонки приводят к загрязнению верхних слоев атмосферы большим коли-чеством мелкодисперсного оптически активного аэрозоля и могут вызвать как региональные, так и глобальные климатические изменения.
Распространение продуктов горения от крупных пожаров с помощью метеорологических моделей дождевых облаков исследовалось в [76, 77],в приближении Буссинеска в [78, 79], с использованием уравнений Навье-Стокса для сжимаемого газа с постоянными коэффициентами турбулентного переноса — в [80]. Формирование конвективной колонки над пожарами исследовалось, например, в [81], струи метеотрона — в [74].
Построение математических моделей тепловых колонок основывается на базе уравнений сохранения макроскопических характеристик вещества струйного потока, в результате чего удаётся численно исследовать динамику формирования струи продуктов горения с учётом фазовых переходов, обусловленных наличием влаги в атмосфере. Очаг пожара моделируется объёмным источником тепла Q(Bt/m3) и массы в(кг/с/м3) мелкодисперсного инертного аэрозоля с заданными законами их изменения по времени. Предполагается,что известны начальные величины Q и S внутри зоны тепловыделения с характерным размером R и высотой h и равны нулю вне этой зоны. При рассмотрении развития турбулентных конвективных движений пылегазовой смеси над очагом пожара во влажной стратифицированной атмосфере учитывается, что влажный вовлекаемый воздух в процессе подъёма и расширения охлаждается. При достижении условий насыщения водяной пар конденсируется с выделением тепла.
Например, в работе [73] формирование и подъём тепловой колонки дымового аэрозоля рассматрива-
Пожар на нефтехимическом производстве.
ется при штилевых условиях в рамках односкорост-ной и однотемпературной модели дисперсной среды, применение которой правомерно, так как размеры дисперсных частиц (дыма, пара, водяных капель) намного превышают характерные молекулярнокине-тические пробеги, а время их скоростной и температурной релаксации значительно меньше времени развития конвективных движений. Кроме того, предполагается малое объёмное содержание дисперсной фазы, не учитываются эффекты столкновения частиц и коагуляция. Конвективные движения среды описываются осесимметричными уравнениями Навье-Стокса в цилиндрической системе координат, тепловыделение пожара имеет площадной характер.
Выполненные по этой модели расчёты показали, что при тепловой мощности источника 5 • 105 Вт/м2 наибольшая вертикальная скорость ( 43 м/с ) наблюдается на оси симметрии;при этом максимальная радиальная скорость у границ очага горения не превышает 17 м/с.
Вовлечение холодного воздуха в восходящую струю наблюдается до высоты 4 км; на высотах от 7 до 11 км образуется зона зависания, в которой дымовой аэрозоль и вовлечённый в струю воздух растекаются в горизонтальном направлении от оси симметрии течения. Через час дымовое облако растекается на площади 700 км2, что почти на порядок больше площади очага горения. Вода, выделяющаяся при конденсации влажного воздуха, может выпасть в виде осадков — дождевых капель, снега и льда.
В работе [74] используется более реалистический подход к моделированию пожара — рассматривается высокотемпературная и высокоскоростная струя в сносящем ветровом потоке. Уравнения сохранения макроскопических характеристик среды записываются в квазиодномерном приближении, не требующем задания детальной информации о начальных и граничных параметрах струи и ветрового потока. Расширение струйного потока интегрально учитывается введением функции вовлечения, имеющей физический смысл массы воздуха, вовлекаемой в колонку в единицу времени через единичную боковую поверхность потока. Газодинамические параметры потока осреднены по поперечному сечению струи. Расчёты,выполненные для мощностей тепловых источников ~ 5 • 10 5
— 5 • 10 6 кВт, показали, что скорости потока достигали 8 — 9 м/с, высоты подъёма до 3 км, а диаметры потоков на завершающих участках подъёма — нескольких сот метров. Причём расчётные значения высоты подъёма струй существенно зависели от степени устойчивости атмосферы и силы ветра и слабо зависели от тепловой мощности источника.
Полученные в этих работах результаты показывают, что динамика формирования конвективной колонки, высота подъёма аэрозоля и характер его распределения в выбросе зависят не только от мощности пожара, но и от влажности атмосферы. Фазовые переходы, вызванные присутствием влаги в атмосфере, существенно влияют на подъём струи и её физические характеристики.
Дальнейшее распространение газообразных и парообразных токсикантов происходит из высотного объёмного источника, которым можно мысленно заменить струйный выброс на завершающем этапе его подъёма. В зависимости от турбулизации атмосферы вторичный высотный источник будет иметь непрерывный или дискретный характер. В любом случае -при непрерывном завершении струи (в виде однородного потока) или в случае её дробления на клубы газообразные и дисперсные составляющие струйного течения рассеиваются в сносящем ветровом потоке под действием атмосферной турбулентности и достигают поверхности земли многократно разбавленные воздухом на больших удалениях как от очага пожара,так и от высотного источника.
Рассмотрению этих достаточно сложных и важных вопросов посвящено большое количество работ, например, [14, 38, 49, 51, 52, 94, 113, 114].
Предлагаемые разными авторами методики вычисления полей распределения концентраций загрязняющих примесей базируются на решении полуэмпирического уравнения диффузии или исходят из рассмотрения гауссовой модели рассеяния. Максимальные значения концентраций примесей у поверхности земли и расстояния от источника загрязнений определяются типом устойчивости атмосферы в слое перемешивания и скоростью ветра.
Использование для расчётов той или иной методики должно приводить к одинаковым результатам при условии корректного задания параметров дисперсии, класса устойчивости атмосферы и периода осреднения используемого набора метеопараметров.
В настоящее время наибольшее распространение получили методика Главной геофизической обсерватории (ГТО) [114] и методика Пасквилла, развитая Ми-дом, Гиффордом и Брайантом [94] и рекомендованная МАГАТЭ для оценки загрязнений атмосферы от предприятий атомной промышленности и энергетики.
Методика Пасквилла позволяет сравнительно просто с использованием таблиц и графиков для дисперсий бокового ау и вертикального az расширения струйного потока определять концентрацию загрязняющей примеси в приземном воздушном слое при мгновенном, разовом и непрерывном выбросах для различных категорий устойчивости и разных расстояний от источника загрязнений в направлении ветрового переноса. Подходы этой методики основаны на гауссовой модели рассеивания загрязняющей примеси и рассматривают три вида точечных источников, отличающихся по продолжительности выбросагмгновенный (одномоментный), кратковременный (от 20 минут до нескольких часов) и постоянно действующий . В соответствии с этой моделью концентрация токсикантов в струе q при выбросе из мгновенного точечного источника мощностью Q в ветровом потоке со средней скоростью и в слое распространения записывается в виде [94]:
где:
tfx> ау> az> — дисперсии распределения концентраций токсиканта в направлениях осей х, у, z — соответственно.
Дисперсии зависят от метеоданных и стратификации атмосферы и являются экспериментально определяемыми функциями времени.
Для кратковременного или постоянно действующего точечного источника пренебрегается диффузией в направлении х по сравнению с ветровым потоком в том же направлении, и выражение для концентрации загрязняющей примеси имеет такой вид:
при получении этой формулы использовался принцип суперпозиции дискретным суммированием результатов, полученных для мгновенного источника.
Причём |
Дисперсии ау
и gz при применении этой методики предлагается находить, используя экспериментальные данные по расширению струи дыма — в вертикальном направлении h и в горизонтальном — углу 0.
для источника продолжительностью от 20 минут до нескольких часов;
для постоянно действующего источника при круговой розе ветров.
В этих формулах:
Н — высота источника загрязнений,м;
0 и 0р — боковое (горизонтальное) расширение струи при мгновенном и продолжительном выбросе, град; h — вертикальное расширение струи, м.
Отметим, что в настоящее время имеется большое количество ведомственных и отраслевых методик расчёта полей загрязнений от приподнятых источников, разработанных как в Росии, так и за рубежом. Однако сопоставить результаты расчётов по ним затруднительно из-за разных наборов используемой входной информации. Обязательными являются разработанные на основе наиболее представительных методик единые Указания по расчёту рассеивания в атмосфере вредных веществ, содержащихся в выбросах предприятий, которые периодически обновляются и модернизируются (например, [68,114]).
В заключение отметим, что основным недостатком рассматриваемых в итературе математических моделей крупного пожара является предположение о химической пассивности загрязняющих примесей. Применительно к проблеме переноса и дисперсии отравляющих веществ неучёт химической деструкции продуктов в нагретом восходящем потоке и при атмосферном его рассеивании многократно увеличивает токсическую опасность последствий пожара [118].