Генератор наноразмерных аэрозолей для тушения пожаров на газовых скважинах и на разливах нефти Текст научной статьи по специальности «Физика»

Вычислительные технологии

Том 12, № 6, 2007

ГЕНЕРАТОР НАНОРАЗМЕРНЫХ АЭРОЗОЛЕЙ ДЛЯ ТУШЕНИЯ ПОЖАРОВ НА ГАЗОВЫХ СКВАЖИНАХ И НА РАЗЛИВАХ НЕФТИ*

А. Д. Рычков, Ю.И. Шокин Институт вычислительных технологий СО РАН, Новосибирск, Россия

e-mail: [email protected]

Numerical modeling is used to investigate an opportunity of application of the pulse aerosol system of fire fighting (PASFF) in the cases of flare smoke jets and oil spills. On the basis of numerical experiment it was established that gasdynamic structure of the flow field can provide an effective transport of aerosol particles of the flame retardant substance and their vapors into the combustion zone for reliable fire extinguishing of heat sources.

Введение

В настоящее время для эффективной борьбы с интенсивными очагами возгорания на больших площадях широкое применение находят системы порошкового пожаротушения на основе различного рода пламегасителей, в качестве которых используются фос-форорганические, металлоорганические, а также некоторые галогенсодержащие соединения. Попадая в зону пламени, эти пламегасящие вещества (ПГВ) разлагаются с образованием активных соединений, которые вступают в химические реакции с радикалами (частицами, ответственными за горение и распространение пламени), переводя их в неактивные частицы. При этом обрываются цепные разветвленные реакции горения, в результате чего и происходит гашение пламени. Эффективность тушения очагов горения при таком способе определяется главным образом показателем концентрации порошка пламегасителя в зоне горения [1], поэтому разработка эффективных способов его доставки в очаг пламени является чрезвычайно актуальной задачей. Заметим, что в соответствии с описанным выше механизмом гашения пламени большую роль в повышении эффективности пожаротушения играет и величина суммарной площади поверхности всех частиц пламегасителя, в связи с чем доставка их в зону горения в виде аэрозольного облака наноразмерных частиц наиболее предпочтительна. Однако получение аэрозолей с такими размерами частиц в существующих импульсных системах порошкового пожаротушения не представляется возможным. Такие системы обычно представляют собой газогенераторы, использующие унитарные аэрозольгенерирующие

* Работа выполнена при финансовой поддержке программы междисциплинарных интеграционных исследований СО РАН (проекты 2006-28, 2006-113), программы фундаментальных исследований ОЭММПУ РАН (проект 4.1.1) и научной школы НШ-9886.2006.9.

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2007.

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

В данной работе для этих целей предлагается использовать импульсную аэрозольную систему пожаротушения (ИАСП) [3, 4], представляющую собой устройство раздельного снаряжения, в котором распыление пламегасящего вещества и транспортировка его в зону горения в виде мелкодисперсного аэрозольного облака осуществляется продуктами сгорания обычного твердотопливного газогенератора кратковременного (импульсного) действия, не содержащими окислитель.

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

По-видимому, наиболее эффективным в этом случае при большой площади возгорания будет кассетное (множественное) применение таких ИАСП, которые можно сбрасывать с вертолета или метать иными способами внутрь очага горения, в результате чего пожар можно быстрее локализовать и потушить.

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

Струя газа Система ИАСП

ф

\ Труба

Газогенератор

Щ Горящий слой нефти

■> / / ? ? /;

о

У

Рис. 1. Схема тушения газового факела

Рис. 2. Схема тушения разлива нефти

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

1. Описание модели и основные уравнения

При работе твердотопливного газогенератора его газообразные продукты сгорания создают давление в несколько десятков атмосфер и их истечение в окружающую среду происходит в турбулентном режиме. Поэтому в качестве модели движения несущего газа использовалась осредненная по Фавру система уравнений Навье—Стокса, замыкаемая ^ — ^)-моделью турбулентности [5]. Движение полидисперной твердой фазы описывалось в рамках стохастического подхода [6], учитывающего влияние турбулентной природы поля течения несущего газа на движение частиц.

Основной целью данной работы являлось изучение газодинамической структуры поля течения при воздействии ИАСП на очаг возгорания и определение полей концентраций частиц ПГВ. Поэтому при построении математической модели принимались следующие допущения:

— совокупность ИАСП (см. рис. 1), воздействующих на газовый факел, заменялась тором прямоугольного сечения, что позволяет заменить пространственное течение осе-симметричным без существенного искажения сути протекающих процессов;

— несущий газ представлял собой химически не реагирующую многокомпонентную смесь, состоящую из воздуха, природного газа, продуктов сгорания газогенератора и продуктов разложения ПГВ;

— течение является турбулентным, двухфазным, осесимметричным и нестационарным;

— частицы ПГВ — полидисперсные, имеют сферическую форму и состоят из N фракций, внутри каждой из фракций частицы имеют одинаковый размер. Столкновениями между частицами на данном этапе моделирования пренебрегалось;

— распределение температуры по объему частицы полагалось однородным. Когда температура достигает заданной величины Т,„ар, происходит разложение ПГВ с образованием паров, причем линейная скорость такого разложения гр предполагалась постоянной;

— газоприход от гранул твердого топлива и от частиц ПГВ при их разложении моделировался источниковыми членами в уравнениях сохранения массы и энергии.

Система уравнений, описывающая такое течение в цилиндрической системе координат, имеет следующий вид:

г= 1

dpYO2 д 1 д

~дГ + дХ (PuYO2 + qx,O2) + rdrr(PvYO2 + qr,O2) = 0; (2)

друпгв , д , _ . s . 1 д , _ s ^ ^

дг + дХ(РиУпгв + Qx,M) + ГдГг(Р^Упгв + 5г,пгв) = — ^ < > ni; (3)

i=1

дрУм д 1 д

—д^--+ дХ (РиУм + qx,M) + ГдГг(Р^Ум + qr,M) = Mgen; (4)

дри д 2 1 д 2 дрк v^ du

^г + тг(ри + p + Txx) + -7T"r(puv + Txr= — > ni < mi — >; (5) дг дх гдг 3 дх ' dt

i=i

др^ д , . 1 д , 2 . 2 дрк ^ dvi

+ Trx) + r(pv + p + Trr) + + 2— = p — > Hi < mi— >; (6)

дГ дх r дг 3 дг r2 ' dt

i= 1

дрЕ д , , ^ . ч 1 д//7^ч .

-7— + — (и(рЕ + p) + UTxx + VTxr + <£x) + "7T"r (^(рЕ + p) + UTxr + VTrr + ) =

дг дх r дr

N

-J2n

i= 1

i

dUi dVi _ . . ,

u < m> +v < mi^- > + < ndiA ■ Nu (T — Ti) > dt dt

дрд д 1 д

+ Q gen-M-geni (7)

дг + дх (puq+Tx'g) + гзтг(р^+Tr,q ) = s; (8) др^ д , 1 д

"дГ + дх (ри^ + ^) + r w + ) = ; (9)

p = рТ^Л + ^ + 1м + -упгв ); YO2 + YN2 + Ум + Упгв = 1; (10)

MO2 Mn2 Мм Мпгв, Sq = Cqi(CMS/^2 — 2D/(3w) — 1)р^д, = [C^ (CMS/w2 — C^D/w) — C^w2,

где q = Vk, w = e/k — "турбулентные" переменные, связанные с кинетической энергией турбулентности k и со скоростью ее диссипации e; Е = CvТ + (u2 + v2)/2 — полная удельная энергия; ^ и — молекулярная и турбулентная вязкости соответственно; p, Т — статическое давление и температура; R0 — газовая постоянная. Величины Mgen и Qgen есть источники массы и энергии, моделирующие приход высокотемпературных продуктов сгорания от газогенератора. Члены в угловых скобках означают осреднение по объему ячейки разностной сетки. Константы, входившие в описание модели турбулентности, имели следующие значения:

Сш1 = 0.55, Сш2 = 0.833, Сшз = 0.666, CM = 0.09, Cq1 = 0.5, ^ = C^/w,

д^ /дм

Txx , Trr , Txr Trx ^e I TJ + TJ ) ,

дх дr \ дr дх /

дТ дТ

0x = — Ae^—, фг = — Ae^— , ^e = ^ + ^t, Ae = A + At, дх дr

_ , ^ дu 1 дrv At = ^tCp/Prtt, D + ,

дх r дr

дУк дУк

qx,k = —(рDfc + setqr,k = —(рDk + set

где индекс k принимает следующие значения: k = O2, N2, M, ПГВ.

Уравнения движения г-й частицы ПГВ вдоль ее траектории записывались в виде

(т.

(Ь

(и + и' — щ) ; (у + у' — Уг);

(щ _ 3 (Ь 4 рр(2

(Уг _ 3 (Ь 4 рр(2

(Т _ 6 Л • Ми (Ь С г Рр (2

(Ь г

(Уг

_ Уг,

(Т — Тг);

11) 12)

13)

14)

15)

16)

где Сг — удельная теплоемкость частицы; хг, уг — ее координаты; пг — концентрация частиц г-й фракции в ячейке в данный момент времени; и', у' — компоненты случайного вектора возмущения скорости газа V' _ (и' ,у ); V' _ \j2kj3 вт/с-1(Мг)У/\У|; V — вектор осредненной скорости; Мг — случайное число из диапазона [—1, +1]; eтfc-1( ) — обратная функция ошибок.

Коэффициенты сопротивления в уравнениях движения частиц рассчитывались по формуле

24 { . Керг/3

Сдг _ ^ ^ + — , еслиКеРг ^ 103' (гР

V — V

Р\ / К>ерг

0.424, если И,ерг > 103, V

(17)

Для расчета течения несущего газа применялась противопотоковая ЬИ-разностная схема второго порядка точности, обладающая ТУБ-свойствами, близкая к схеме из работы [7]. Расчет движения частиц проводился с помощью А-устойчивой разностной схемы также второго порядка точности [8]. Область решения системы уравнений (1)—(10) снизу ограничивалась непроницаемой поверхностью (подстилающая поверхность, ось ОУ направлена в горизонтальном направлении, ось ОХ — в вертикальном), слева — осью симметрии течения, справа и сверху не ограничивалась ничем (безграничное пространство). Граничные условия для несущего газа задавались следующим образом. На оси ОХ задавались условия симметрии течения, на поверхностях ИАСП и на подстилающей поверхности — условия прилипания, на остальных границах — неотражающие граничные условия. Для системы (11)—(16) на подстилающей поверхности задавались условия неупругого отражения:

5

у- _ —(0.993 — 1.760 + 1.5602 — О.4903)Уг+га, V- _ 7УУ+,

где Угп, УЦ~п, Уг"т, Щт — нормальная и тангенциальная компоненты вектора скорости г-й частицы, символ (+) относится к падающей и символ (—) — к отраженной частицам; 0 — угол падения.

Величина массоприхода Мдеп продуктов сгорания газогенератора в (1) определялась решением уравнения

М = —рдг^ • г'ьд^дд'

где рдг, S, гЬд — плотность материала гранулы твердого топлива, ее текущая поверхность и скорость горения; N0 — число гранул в единице объема заряда твердого топлива, которое остается постоянным в процессе их горения и определяется из начальных условий.

2. Некоторые результаты расчетов

Во всех приведенных ниже расчетах скорость испарения частиц гр = 10 мм/с, число фракций полидисперсной среды равнялось 10, их среднемассовый размер (143 = 20 мкм, минимальный размер частиц равен 0.1 мкм, температура испарения частиц ПГВ Туар = 400 К, теплотворность твердого топлива газогенератора Qgen = 1780 кДж/кг.

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

Рис. 3. Структура поля течения и изотермы Рис. 4. Структура поля течения и изотермы в холодной струе (Ь =10 мс) в холодной струе (Ь =40 мс)

такой струи моделировался сверхзвуковой холодной струей природного газа с числом Маха, равным 3. Диаметр струи был равен 0.5 м, температура в струе Т = 320 К.

Внутренний и внешний диаметры тора ИАСП равны 0.9 ми 1.06 м соответственно, масса пламегасителя — 60 кг, масса заряда твердого топлива — 2 кг.

На рис. 3 и 4 показаны распределения температуры в поле течения и положения частиц пламегасящего вещества в два различных момента времени. Видно, что на начальном этапе воздействия ИАСП на факел сначала происходят "развал" струи, что способствует быстрому отсеканию поступления газа в горящую часть факела, и формирование в центре струи зоны частиц ПГВ с достаточно высокой их концентрацией. После прекращения работы ИАСП структура струи частично восстанавливается и осуществляется транспорт частиц пламегасителя в зону горения.

Горящий участок факела моделировался горячей дозвуковой струей нереагирующе-го природного газа с числом Маха = 0.4. В этом случае из-за теплового расширения струи ее диаметр полагался равным 1 м, температура в струе Т = 1200 К. Внутренний и внешний диаметры тора ИАСП были равны 1.4 м и 1.56 м соответственно, масса пламегасителя — 90 кг, масса заряда твердого топлива — 3 кг. На рис. 5 и 6 представлены распределения температуры в поле течения и положения частиц ПГВ на начальном и конечном этапах воздействия ИАСП на струю. Видно, что на начальном этапе воздействия дозвуковая струя, в отличие от сверхзвуковой, теряет свою устойчивость, что также приводит к ее "развалу" и к существенной перестройке течения в ней, а также к затягиванию процесса ее "релаксации", что положительно влияет на процесс гашения пламени.

Как уже говорилось, для тушения разливов горящей нефти предлагалась схема ИАСП в виде цилиндра, в центральной части которого располагался заряд твердого топлива газогенератора (см. рис. 2). В этом случае очаг горения моделировался воздухом, нагретым до температуры Т = 800 К, ИАСП располагался на расстоянии Н = 0.1 м от подстилающей поверхности (разлива нефти) и основной целью исследования явля-

1ПЯЗ ЧйП й/н 77й к/н ЧЯП 1й» "шк 1186. 1083. 980. 878. 778. 678. 580. 182. 306

Рис. 5. Течение в дозвуковой струе (г = 20 мс)

Рис. 6. Течение в дозвуковой струе (г = 50 мс)

2388. 2215. 2065. 1922. 1784. 1650. 1518. 1389. 1261. 1135. 1011. 887. 765. 644. 523. 307. 2.00 1.83 1.69 1.55 1.42 1.29 1.17 1.04 0.92 0.80 0.68 0.56 0.45 0.33 0.21 0.01

Рис. 7. Поле течения при тушении разлива Рис. 8. Концентрация паров ПГВ

нефти

лось определение структуры аэрозольного облака частиц и распределения концентрации паров ПГВ. Диаметр цилиндра ИАСП был равен 0.36 м, высота цилиндра — 0.2 м, масса заряда твердого топлива — 0.1 кг, масса ПГВ — 10 кг. На рис. 7 показаны структура поля течения и положения частиц ПГВ в момент времени £ = 10 мс. Распределение температуры свидетельствует о ее понижении в области, занятой частицами ПГВ, из-за поглощения ими тепла, начальная температура которых была равна 300 К, а температура испарения — ниже температуры окружающей среды. Вихревая структура за облаком частиц обеспечивает устойчивую компактную форму облака аэрозоля в форме тора, генерирующую достаточно высокую концентрацию паров ПГВ (рис. 8). Видно, что предлагаемое устройство способно осуществлять транспортировку достаточно большой массы аэрозольного облака ПГВ в зону горения без заметного разрушения его структуры продуктами сгорания газогенератора.

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

[1] Баратов А.Н., ВогмАн Л.П. Огнетушащие порошковые составы. М.: Стройиздат, 1982. 72 с.

[2] ИсАвнин Н.В. Средства порошкового пожаротушения. М.: Стройиздат, 1983. 156 с.

[3] Рычков А.Д., Шокин Ю.И. Моделирование работы генератора аэрозолей в качестве пламегасящего устройства // Матер. XV междунар. конф. по вычисл. механике и современным прикладным программным системам (ВМСППС-2007), 25-31 мая. Алушта, 2007. С. 446-448.

[4] Рычков А.Д., Шокин Ю.И., Милошевич Х. Исследование возможности применения твердотопливных газогенераторов для тушения пожаров на газовых скважинах // Математическое моделирование научно-технологических и экологических проблем в нефтегазодобывающей промышленности (Матер. VI Казахстанско-российской междунар. научно-практ. конф. Астана: Изд-во Евразийского ун-та, 2007. С. 274-278.

[5] Cqakley T.J. Turbulence modeling for high speed flows // AIAA Paper. 1992. 92-0436.

[6] Crow C.T. Review — Numerical models for dilute gas-particles flows // Trans. of ASME J. of Fluid Eng. 1982. Vol. 104. P. 297-303.

[7] Yoon S., Jameson A. An LU-SSOR scheme for the Euler and Navier-Stokes equations // AIAA Paper. 1987. 87-600.

[8] Рычков А.Д. Математическое моделирование газодинамических процессов в каналах и соплах. Новосибирск: Наука, 1988. 222 с.

Поступила в редакцию 12 ноября 2007 г.