Автоматизированное оперативное управление техногенными химико-технологическими объектами при возникновении запроектных аварийных ситуаций Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

Труды ИСП РАН, том 27, вып. 6, 2015 г..

Автоматизированное оперативное управление техногенными химико-технологическими объектами при возникновении запроектных аварийных ситуаций1

Ю.Н. Матвеев <[email protected]>

Я.А. Стукалова <[email protected]>

Тверской государственный технический университет,

170026, Россия, наб. А. Никитина, 22

Аннотация. В статье обосновывается необходимость применения распределенных вычислительных систем для математического моделирования процесса образования облака зараженного воздуха при возникновении запроектных аварийных ситуаций на техногенных химико-технологических объектах. Проанализирована зависимость скорости изменения концентрации опасной примеси в произвольной точке пространства и приведены математические модели процессов образования облака зараженного воздуха при высокотемпературных выбросах (взрывах, пожарах) и проливах больших количеств токсичных химических веществ на различные поверхности с последующим их испарением. Доказана целесообразность декомпозиции задачи моделирования процесса образования облака зараженного воздуха.

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

1.Введение

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

1 Работа выполнена при финансовой поддержке РФФИ грант №15-29-07970

395

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

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

2. Классификация аварийных ситуаций

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

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

К «проектным авариям» относятся такие ситуации, которые в случае их возникновения не приводят к аварии т.к. заранее запланированы дополнительные организационные и технические мероприятия по их нейтрализации. Как правило, это аварийные ситуации, причиной которых являются различного рода отказы оборудования. Для снижения вероятности и возможных масштабов таких аварий в типовую проектную и технологическую документацию вносят определённые дополнения, предусматривается установка различного рода блокирующих, сигнальных и других систем.

К «запроектным» относятся аварии, не вошедшие в первую группу. Причинами таких аварий служат в основном внешние непредсказуемые заранее события - различного рода стихийные бедствия (землетрясение, тайфуны, ураганы и т.д.) или непрогнозируемые последствия человеческой деятельности (взрывы, крупные пожары и пр.). Возможность снижения вероятности их возникновения учитывается в основном за счёт таких организационных мероприятий, как выбор площадки для размещения объекта,

396

Труды ИСП РАН, том 27, вып. 6, 2015 г..

например, в сейсмически неактивной зоне. В проектной и технологической документации возможность возникновения таких аварий не находит своего отражения, поэтому такие аварии принято называть «запроектными». Вероятность возникновения проектных аварий гораздо выше, чем аварий запроектных. Однако масштабы запроектных аварий в случае их возникновения гораздо больше, чем масштабы проектных, так как при малой частоте их появления они обладают гораздо большей разрушительной силой. Причём, по упомянутым выше причинам, на техногенных объектах не предусматривается введение элементов и систем, предназначенных для противостояния или ограничения этой разрушительной силы.

Анализ возможных сценариев аварийных ситуаций на техногенных объектах химических производств, приводящих к выбросам токсичных химических веществ (ТХВ) в атмосферу с образованием облака зараженного воздуха (ОЗВ), показывает, что основными вариантами сценариев таких аварий являются[1]:

высокотемпературные выбросы ТХВ, которые по времени протекания могут быть кратковременными и продолжительными (взрывы, пожары);

пролив больших количеств ТХВ на различные поверхности с последующим их испарением.

3. Процессы образования ОЗВ как объекты управления.

Зависимость скорости изменения концентрации примеси в произвольной точке пространства - dCa/dt определяется расположением в пространстве источников примеси и рядом параметров, в число которых входят составляющие скорости ветра - и, v, w вдоль осей Ox, Оу, Oz трехмерного пространства, коэффициенты атмосферной турбулентности - к и др. [2]

В общем виде эта зависимость определяется дифференциальным уравнением

баланса примеси или уравнением переноса примеси:

дСа дСа дСа дСа д j дСа д т дСа.

---= -U------V------W------+ — (К---)+ — (К-----) +

dt дх ду dz дх дх ду ду

д т дСа

Л (К-----1-WaCa) + F + R - Р-W = f(u,v,-w,k,w.,F,R,P,W)

dz dz

где:

Wa-собственная вертикальная скорость примеси,

R и Р - скорость образования и уничтожения примеси в результате химических реакций,

W - скорость выпадения примеси на подстилающую поверхность;

397

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

Ca -количество примеси a, содержащейся в единице объема воздуха (объемная концентрация примеси);

F - скорость поступления примеси в воздух от того или иного источника (интенсивность источника).

Уравнение (1) показывает, что перенос примеси в движущейся среде обусловлен двумя различными физическими факторами: во-первых, при наличии разности концентраций идет процесс молекулярной диффузии; во-вторых, частицы примеси увлекаются движущейся средой и переносятся вместе с ней.

Кроме параметров, входящих в уравнение (1) должны быть заданы условия, ограничивающие область распространения ОЗВ. Если область распространения ОЗВ ограничена сверху поверхностью z=H, а снизу земной поверхностью z=0, то условия на этих границах обычно задаются в виде:

при z=H Са=0 или к ° =0, dz (2)

при z=0 к дСа +waCa=pCa. dz (3)

При рассмотрении примеси в слое толщиной Н порядка 2-5 км условия (2) определяют исчезновение примеси или ее вертикального потока на этой границе. Условие (3) на земной поверхности имеет смысл равенства вертикального турбулентного потока и потока примеси при ее гравитационном оседании на поверхность со скоростью wa количеству примеси, поглощаемой поверхностью - рСа, где Р - коэффициент аккомодации, зависящий от физических свойств подстилающей поверхности, наличия на ней растительности или застроек. Тогда задачу оптимизации управления процессом образования ОЗВ в общем виде можно записать:

V0 = F\X, U, Ca, t) —> min (4)

т. е. минимизировать объем ОЗВ при связях и ограничениях:

dCaldt=f(u,v,w,k,Wa,F,R,P,W)=f( X ,U ), (5)

Ca(x,y,z,t;ho) < ПДК, (6)

^min — % ^Х-тах ? (7)

У min —У ~Уmax-’ (8)

0< z <гф, (9)

398

Труды ИСП РАН, том 27, вып. 6, 2015 г..

о < t<tdup,

(10)

где

Ca(x,y,z,t;ho) - концентрация ТВХ в точке с координатами x,y,z в момент времени t от источника, расположенного на высоте ho;

X, U - векторы входных и управляющих параметров;

2ф - высота флюгера;

tdup - директивное время локализации аварии.

4. Математические модели процессов образования облака зараженного воздуха

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

При взрыве процессом поступления ТХВ в атмосферу будет испарение жидких частиц и капель с поверхности образовавшегося пароаэрозольного облака. Поступление ТХВ в атмосферу определяется производительностью источника [3]:

Ло „ q==EoQX р

то

^ akEot^

V

то У

(П)

где:

q - производительность источника, кг/с;

- начальная плотность заражения; то - масса частицы среднемассовым диаметром dk, кг;

Е0 - начальная скорость испарения, кг/(м2-с):

Eq = К p27T^ + l,2~2RQ2/3)coD dny/o/D , (12)

где:

Кр - коэффициент, зависящий от типа ТХВ и вида почвы (КР= 0,9 для большинства ТХВ);

399

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

Re - число Рейнольдса:

_ 0,25 midn . va

C0— максимальная концентрация насыщенного пара, % ; D - коэффициент диффузии вещества; d„ - диаметр пятна от частицы вещества, м :

(13)

dn~^>dk (14)

D - коэффициент кинематической вязкости воздуха;

cik - параметр, зависящий от доли свободной или связанной жидкости;

t - время, прошедшее от начала заражения, с.

При проливе процессом поступления ТХВ в атмосферу будет только испарение с поверхности образовавшейся «лужи»:

q = S пр Е (15)

где:

q - производительность источника, кг/с; Е - скорость испарения, кг/(м 2-с);

Snp - площадь зеркала пролива, м2.

Испарение с гладкой поверхности пролитого вещества рассчитывается по методу Братсерта [4].

Скорость стационарного испарения жидкого ТХВ со свободной поверхности пролива, когда необходимо оценивать общую массу испарившегося вещества, можно определить из выражения [5]:

/_ V/3

Е = 0,0735-Ст'U*

Рм V °в

кг

(16)

где:

Ст - максимальная концентрация паров ТХВ при данной температуре испарения, кг/м 5;

- динамическая скорость воздуха над поверхностью испарения, м/с,

- коэффициент молекулярной диффузии паров ТХВ в воздухе, м2/с\

400

Труды ИСП РАН, том 27, вып. 6, 2015 г..

Ув - кинематическая вязкость воздуха, м2 /с.

Основное уравнение для расчета Ст применительно к ТХВ имеет вид:

Cm = 16 • 109 Мов^ Рт , кг/М3 (17)

где:

Мое ~ молекулярная масса ТХВ, кг/моль;

Рт - давление насыщенного пара, мм рт.ст. при температуре Т.

Давление насыщенного пара ТХВ при температурах, характерных для процессов испарения ТХВ при авариях, может быть рассчитано по уравнению

LqЮ (рт) — A J,

(18)

где:

А, В - коэффициенты уравнения Антуана, рассчитываемые по экспериментальным данным для каждого типа ОВ в определенном интервале температур;

Т - температура жидкого ТХВ, град.К.

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

и

,/и

(0,55 - 0,52в°ш™) ехр (о, 52 + 0,283в°'™) In z„,В > 0;

< 0,55 • zT\B = 0 (19)

[(о, 55 +0,52(-В)ода)]ехр[(о, 52 + 0,0313(-в)0'155)] ]nz„,B <0,

где:

В - значение параметра Бызовой (табл. 1);

Z0 - шероховатость подстилающей поверхности, м (табл.2) и/— скорость ветра на высоте \м над слоем шероховатости, рассчитываемая по формуле:

mi = M^[(l+Zo)/x^] , (20)

где:

Z0 - высота флюгера, м:

401

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

иф - скорость ветра на высоте флюгера, м/с:

т - показатель профиля скорости ветра в приземном слое атмосферы, который рассчитывается по формулам:

т = <

(го/гф ~ l) 1п(го

А ф)

^(l-zo/гф)

{l-zo/гфУ

,В = 0; ,В*0.

(21)

Табл. 1 Значения коэффициента Бызовой

Характер истика устойчив ости атомосфе ры Очень неустой чивая Умерен но неустой чивая Слабо неустой чивая Безразл ичная Слабо неустой чивая Умерен но неустой чивая Очень неустой чивая

Класс устойчив ости по Паскуилу -Тернеру 1 2 3 4 5 6 7

В -0,1 -0,01 -0,001 0,0 0,001 0,01 0,05

у

Табл. 2 Шероховатость подстилающей поверхности ^ 0 в зависимости от типа местности, где происходит рассеяние паров ТХВ

Тип поверхности Z0, м

Ровная местность, покрытая снегом 0,001

Ровная местность с высотой травы до 0,001

0,1

Ровная местность с высотой травы до 0,01

0,15 м

Ровная местность с высотой травы до 0,05

402

Труды ИСП РАН, том 27, вып. 6, 2015 г..

0,6 м

Неровная местность, покрытая кустарником 0,1

Лес высотой до 10 м 0,4

Г ородская застройка 1,0

Коэффициент молекулярной диффузии DM может быть получен из уравнения Андрусова [4]:

Dou

17,2(l + -Jmb +Мт )

Ра(з,08+У1ру ■ у/мв-Мт

(22)

где:

Dqm - коэффициент диффузии при 0 °С, м 2/с:

Мв, Мг- молекулярная масса воздуха и ТХВ соответственно, г/моль; Ра— давление, атм.

Для температур, отличаюпщхся от 0°С, используется уравнение вида [4]:

Dm (/ ) ~ Dqm

(23)

где Т - температура.

Расчетное уравнение для оценки величины кинематической вязкости воздуха при атмосферном давлении имеет вид [4]:

ов = 2,84 • 10 6 • Т • ехр(4,36 + 2,844 • IQ-3 • Т) (24)

При пожаре образование ОЗВ будет осуществляется следующими процессами: выброс паров ТХВ в атмосферу в результате взрывного разрушения емкостей при термическом воздействии (через 20 минут после возникновения пожара) и конвективный подъем в очаге пожара; поступление паров ТХВ в атмосферу в результате конвективного подъема в очаге пожара и испарения с поверхности пролива при повышенной температуре.

Методика количественной оценки выброса паров ТХВ при взрывном разрушении емкости под воздействием тепловой нагрузки в [4] отсутствует. Поэтому можно ориентировочно принять массу ТХВ, поступившего в ОЗВ, равной 10 ... 20% от общего количества ТХВ в емкости.

Мв = (0,1...0,2)Мя, (25)

403

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

где:

MB - масса паров ТХВ, поступивших в ОЗВ при термическом взрыве в момент времени tB, кг;

МЕ - масса ТХВ в емкости, кг.

Производительность источника поступления ТХВ в ОЗВ рассчитывается по формулам (15)-(22).

Высота подъема конвективной колонки может быть определена из соотношения [6]:

hk -

Кг Rn

U

2,5 + 3,3 ■ g-Rn-{Tn-TB) Tb-U2

(26)

где:

hk- высота источника относительно поверхности земли, м:

Кг - размерный коэффициент (равный 1,5 м/с):

Rn - радиус площади зоны пожара, м:

U - скорость ветра на высоте 10 метров, м/с, g - ускорение силы тяжести;

Т

1 п - температура горения в абсолютной шкале, °К;

Тв - температура окружающего воздуха в абсолютной шкале, ° К.

Интенсивность выделения i-ro продукта горения находится по формуле

Qi = ГПг ■ са (27)

где

Qi - интенсивность выделения i-ro продукта горения, кг/{м 2-с);

да, - массовая скорость выгорания, кг/{м 2-с):

а , - весовая доля i-ro компонента продукта горения, выделяющегося

при горении единицы массы горючей нагрузки (справочные данные).

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

404

Труды ИСП РАН, том 27, вып. 6, 2015 г.

На основании проведенного анализа процессов перехода ТХВ из начального источника химического заражения в ОЗВ можно сделать вывод, что переход ТХВ осуществляется в результате следующих процессов:

1. мгновенный выброс паров ТХВ в атмосферный воздух;

2. непрерывное поступление ТХВ в атмосферу - испарение из образовавшейся на поверхности (на неограниченной поверхности, в обваловке, в поддоне) «лужи» с ТХВ и испарение с поверхности выседания жидких частиц и капель;

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

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

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

5. Заключение

Рассмотренный в статье подход к моделированию процесса образования ОЗВ предполагает введение существенных допущений и ограничений. Например, величины компонент вектора скорости ветра по координатным осям OX,OY приняты постоянными и независящими от высоты.

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

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

Список литературы

[ 1 ]. Методика выявления и оценки химической обстановки при разрушении (аварии) объектов, содержащих СДЯВ.-М: ГШВС.1989-116с.

[2] . Батырев В.В., Минько С.М. и др. Обоснование размеров зоны защитных

мероприятий вокруг объектов по хранению и уничтожению химического оружия. Отчет по НИР «Вагонетка-О», этап 2,- Новогорск: АГЗ, 1998.-184с.

[3] . Белов П.Н. Теория расчета распространения атмосферных примесей и ее

применение для оценки загрязнения природной среды // Географическое

405

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

прогнозирование и охрана природы: Сб.науч.тр./ Под ред. Звонковой Т.В., Касимова Н.С. -М.:Изд-во МГУ, 1990. - 176с.

[4] . Братсерт У.Х. Испарение в атмосферу. Теория, история, приложения./ Пер. с

англ. Под ред. А.С. Дубова. - Л.: ГИМЗ, 1985. - 350 с.

[5] . Методика прогнозирования развития и последствий аварийных ситуаций на

объектах уничтожения химического оружия. Шифр «Система» (первая редакция). М.: ГУП ГосНИИОХТ, 2001. - 123 с.

[6] . Разработка и обоснование перечня мероприятий, осуществляемых при

ликвидации последствий чрезвычайных ситуаций, применительно к возможным сценариям их развития. Промежуточный отчет по НИР «Конда», этап2, № 715 -Редкино: ОАО «РОКБА», 2001 - 142 с.

406

Труды ИСП РАН, том 27, вып. 6, 2015 г..

Computer-aided operational management technogenic chemical-technological objects at occurrence of beyond design basis emergency situations

Y. Matveev <matveev4700(a)ynail.ru>

N. Stukalova <nast77(a).mail.ru>

Tver state technical University, 170026, Russia, A. Nikitin emb., 22

Abstract. The article substantiates the necessity of using distributed computing systems for mathematical simulation of process of formation of cloud contaminated air at occurrence of beyond design basis emergency situations of technogenic chemical-technological objects. We analyzed the dependence of the rate of change of concentration of harmful impurities in the arbitrary point of the space and the mathematical model of processes of formation of cloud contaminated air at high temperature releases (explosions, fires) and the Straits of large quantities of toxic chemicals on various surfaces and their subsequent evaporation. Proved the feasibility of decomposing the problem of modeling of process of formation of the cloud of contaminated air.

Keywords: emergency; chemical engineering; cloud contaminated air, mathematical model; distributed computing system.

References

[1] . Metodika vyjavlenija i ocenki himicheskoj obstanovki pri razrushenii (avarii) obektov,

soderzhashhih SDJaV [Methods of identification, M: General Staff of the Armed Forces. 1989-116 p (in Russian)

[2] . Batyrev V.V., Minko S.M. i dr. Obosnovanie razmerov zony zashhitnyh meroprijatij

vokrug obektov po hraneniju i unichtozheniju himicheskogo oruzhija [Rationale for the size of the zone of protective actions around objects on storage and destruction of chemical weapons], Otchet po NIR «Vagonetka-О», jetap 2. [Scientific research reports "Trolley-On", step 2.] -Novogorsk: AGZ, 1998.-184 p. (in Russian)

[3] . Belov P.N. Teorija rascheta rasprostranenija atmosfemyh primesej i ее primenenie dlja

ocenki zagijaznenija prirodnoj sredy [The theory for calculating the distribution of atmospheric constituents and its application to evaluation of environmental pollution], Geograficheskoe prognozirovanie i ohrana prirody: Sb.nauch.tr. Pod red. Zvonkovoj T.V., Kasimova N.S. - M.:Izd-vo MGU[MSU Publishing house], 1990. - 176 p. (in Russian)

407

Trudy ISP RAN [The Proceedings of ISP RAS], vol. 27, issue 6, 2015.

[4] . Bratsert U.H. Isparenie v atmosferu. Teorija, istorija, prilozhenija [ Evaporation into the

atmosphere. Theory, history, applications]. Per. s angl. Pod red. A.S. Dubova. - L.: GIMZ [GIS], 1985. - 350 p. (in Russian)

[5] . Metodika prognozirovanija razvitija i posledstvij avarijnyh situacij na ob#ektah

unichtozhenija himicheskogo oruzhija. Shifr «Sistema» (pervaja redakcija) [Methods for prediction of development and consequences of emergency situations on objects of chemical weapons destruction. The code "System" (first edition)]. M.: GUP GosNIIOHT [sue GosNIIOKhT], 2001. - 123 p. (in Russian)

[6] . Razrabotka i obosnovanie perechnja meroprijatij, osushhestvljaemyh pri likvidacii

posledstvij chrezvychajnyh situacij, primenitel'no к vozmozhnym scenarijam ih razvitija. Promezhutochnyj otchet po NIR «Konda», jetap 2 [Development and validation of the list of activities undertaken during the liquidation of consequences of emergency situations, with respect to possible scenarios of their development. Interim report on research "Condat", phase 2], № 715 - Redkino: ОАО «ROKBA» [JSC "ROKBA"], 2001 - 142 p. (in Russian)

408