Количественный критерий классификации аварийных газовых выбросов Текст научной статьи по специальности «Математика»

ПОЖАРОВЗРЫВОБЕЗОПАСНОСТЬ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ

С. Е. Якуш

д-р физ.-мат. наук, ведущий научный сотрудник, Институт проблем механики им. А. Ю. Ишлинского РАН, г. Москва, Россия

Г. М. Махвиладзе

профессор, руководитель Центра по исследованиям пожаров Универ-си тета центрального Ланкашира, г. Престон, Великобритания

УДК 532.525.2:614.841.2:614.841.12

КОЛИЧЕСТВЕННЫМ КРИТЕРИИ КЛАССИФИКАЦИИ АВАРИЙНЫХ ГАЗОВЫХ ВЫБРОСОВ

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

Ключевые слова: разгерметизация, аварийные выбросы, горение, огненный шар.

Введение

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

К настоящему времени разработаны методики инженерного расчета перечисленных явлений, позволяющие определять зоны действия поражающих факторов [1-5]. Однако критерии реализации того или иного типа истечения в принятых методиках отсутствуют либо являются качественными. Так, при практическом применении инженерных формул для оценки последствий выбросов сжатых газов возникает вопрос о том, при каких размерах выходного отверстия выброс следует относить к мгновенному типу, а при каких — к струйному. Расследование происшедших аварий может привести к обратной задаче: по известным размерам отверстия в оболочке резервуара определить, какой вид горения реализовался в ходе аварии — факел или огненный шар.

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

Критерий классификации газовых выбросов

Для получения количественного критерия классификации выбросов используем физическую модель [6-8], основанную на сопоставлении времени истечения и характерного времени турбулентного смешения выбрасываемого газа с атмосферным воздухом. Рассмотрим истечение заданной массы газа М0 в атмосферу через круглое отверстие диаметром В. Если время истечения ^ значительно меньше времени смешения, то результатом выброса будет переобогащенное топливом облако, сгорающее при зажигании в режиме огненного шара. В другом предельном случае, когда время выброса значительно превосходит время смешения, в атмосфере образуется струя газа, сгорающая при зажигании в виде факела. С целью нахождения параметров, при которых реализуется тот или иной тип истечения,

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

Время истечения газа

Мп

Мп

кР/ 4

(1)

сd Р ъиъ

А я

Бг

Диаметр отверстия

Рис. 1. Схема классификации выбросов конечной продолжительности с разделением их критическими диаметрами Рс и Ру на выбросы с образованием облака (а), струи (в) и промежуточные истечения (б)

где О — массовый расход газа в единицу времени; Cd — коэффициент истечения; ръ и иъ — плотность и скорость газа в выходном отверстии соответственно. Время турбулентного смешения в облаке направленного мгновенного выброса tс определяется по формуле

tс =

0,5Р

иъ

-1/3/

т,

л

4/3

V т

X

-4/3 ВКПВ,

(2)

где а — индекс, обозначающий величины, относящиеся к окружающей атмосфере;

— молекулярная масса истекающего газа. Вывод формулы (2) на основе интегральной модели и определение констант из экспериментальных данных подробно рассмотрены в [6, 7].

Время смешения в стартующей турбулентной струе tJ вычисляется по формуле

ч 2

3,73Р

и,

-1/2/

X

-2

ВКПВ ,

(3)

е У

где числовой множитель вычислен в [6, 7] из условия, что передний фронт за время у достигает точки, в которой осевая относительная объемная концентрация газа равнаХВКПВ.

Из формул (2) и (3) следует, что

^ = 0,134 I ^ У V Ръ

1/6/.

Х2/3 = 1 Х ВКПВ = 1,

где для типичных углеводородовХВКПВ = 0,05^0,20, т. е. разбавление газа воздухом при мгновенном выбросе происходит значительно интенсивнее.

Времена смешения по (2) и (3) пропорциональны Р, тогда как время истечения (1) убывает с ростом диаметра как Р 2 (рис. 1). Точки пересечения

tr = tс и tr = у дают два критических диаметра — Рс и Ру разделяющих выбросы различного типа. Критерий классификации выбросов формулируется следующим образом [6, 7]:

• tr < tс при Р > Рс — "мгновенный" выброс с образованием переобогащенного топливом облака с концентрацией топлива выше ХВКПВ;

• tr > tJ при Р < Ру — "непрерывное" истечение с образованием квазистационарной топливной струи;

• tс < tr < у при Рс > Р > Ру — промежуточный случай, развивающееся топливовоздушное облако, в котором за время истечения существенная часть газа разбавляется до концентраций ниже хвкпв.

Конкретный вид формул для критических диаметров Рс и Ру зависит от условий истечения, при этом, как следует из (1)-(3), все характерные времена пропорциональны и-1, поэтому скорость газа не входит в окончательные выражения.

При низком давлении в резервуаре истечение можно считать изотермическим, при этом плотности газа в резервуаре р0 и в выходном отверстии ръ можно принять равными плотности газа при наруж-

ных температуре и давлении: р0 ; случае из (1)-(3) получаем:

Ръ ~ Ре,а. В этом

РЬхр = 3

та

V а У

Х4/3 у .

Х ВКПВ у 0 ;

р3,ьр = 0,4

/ \ 3/2

т„

(4)

та V а У

Х ВКПВ у0,

где у0 — объем выбрасываемого газа при условиях хранения; У0 = М0/р0; константы вычислены при Cd = 0,85; индекс "ЬР" означает "низкое давление".

3 =

а

tJ =

а

При высоком давлении в резервуаре (Р0 > Р») происходит запирание потока, при этом давление и плотность газа в выходном сечении описываются формулами [9]:

Рь = Ро

2 ^у-1 ( 2 ^у-1

; Р ь = Р о I

у + 1

Р = Р„

у + 1

У + 1V-1

"Г

где у — отношение удельных теплоемкостей газа.

Дорасширение газа учтем введением эквивалентного диаметра Вед = В (Рь/Ра)1/2 и плотности ред = р(Ра/Рь), после подстановки которых вместо истинных величин получим (индекс "НР " означает "высокое давление"):

В

с, НР

= В

у + 8

3 (у + 1 ^ 6(у-1) ( Ра

с, ЬР

Р

В

У, НР

= В

3 ( у + 1 ) 2(у -1) ( Р

У, ЬР

Р

1/2

1/2

(5)

Важно отметить, что входящие в (5) величины ВсЬР и ВУЬР следует вычислять по формулам (4), при этом под ¥0 понимается объем выброшенного газа при условиях хранения, т. е. при высоком давлении Ро.

Формулы (5) получены в предположении о постоянстве расхода газа в течение всего времени истечения. Реализация такого случая возможна, например, если давление в резервуаре поддерживается за счет подачи газа из каких-либо сторонних источников. Если же происходит разгерметизация и опорожнение резервуара конечного объема, то давление и расход падают со временем. Чтобы учесть переменность расхода, заменим переменный расход его средним за время истечения значением Оау = цО0, где О0 — начальный расход; ^ — коэффициент пропорциональности, зависящий от начального давления, как показано в [7]; с приемлемой точностью можно принять ^ = 0,6(Ра /Р0)1/6. В результате приходим к выражениям для критических диаметров при истечении из резервуаров конечного объема (индекс "ГУ"):

В

3

с, ГУ

= 2,15В

у+ 8

3 ( у + 1 ^ 6(у-1) ( Р С, ЬР

В

3

У, ГУ

= 2,15В

3

У, ЬР

у + 1 ^ 2(у -1) ( Ра

1/4

1/4

(6)

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

Низкое

Конечный Высокое объем

2,0

1,5

1,0

0,5

Рис. 2. Критерий классификации выбросов конечной продолжительности на плоскости параметров 5 - 3. Граничные линии 5с А) и 5У А) разделяют выбросы с образованием облака (а), струи (в) и промежуточные истечения (б). Вертикальными линиями показаны значения параметра 5 для веществ, перечисленных в таблице

мерный диаметр отверстия 3 и параметр свойств вещества 5:

5 =

В (Р.

У

1/3

РI; 5 =

а

Ш„

та

,1/2

2/3

вкпв,

(7)

где п = 0 — для выбросов низкого давления; п = 1/6 — для выбросов высокого давления; п = 1/12 — для истечений из резервуаров конечного объема.

В переменных (7) границы между различными типами выбросов описываются простыми соотношениями:

5 с = ^ с 5

2/3 .

5 у =Х У 5.

(8)

Приведем значения коэффициентов формул (8), вычисленные при у = 1,4, для выбросов из резервуаров: низкого давления — Хс = 1,44, = 0,74; высокого давления — Хс = 1,83, = 0,93; конечного объема— Хс = 2,36,= 1,20. Граничные кривые (8) представлены на рис. 2. Значения критических диаметров для ряда пожароопасных газов приведены в таблице [7].

Горение газовых выбросов

Если выброс горючего газа сопровождается его зажиганием, то возникает задача оценки зон теплового поражения. В то время как в предельных случаях непрерывного и мгновенного истечения можно применять существующие формулы для факельного горения и огненного шара соответственно [1-5], для оценки истечений, попадающих согласно предложенной выше классификации в категорию

3

п

0

3

0

Критические диаметры для некоторых пожароопасных веществ [7]

Вещество % Истечение при низком давлении Истечение при высоком давлении Истечение из резервуара конечного объема

sc s^ sc s^ sc

Водород Н2 0,219 0,162 0,522 0,204 0,664 0,263 0,857

Метан СН4 0,210 0,155 0,509 0,195 0,647 0,252 0,834

Этан С2Н6 0,287 0,212 0,627 0,267 0,796 0,344 1,030

Пропан С3Н8 0,257 0,190 0,582 0,239 0,739 0,308 0,953

Бутан С4Н10 0,271 0,201 0,603 0,252 0,767 0,325 0,989

Бутен-1 С4Н8 0,286 0,212 0,625 0,266 0,794 0,343 1,020

Пропилен С3Н6 0,289 0,214 0,629 0,269 0,800 0,347 1,030

Аммиак КН3 0,327 0,242 0,684 0,304 0,869 0,392 1,120

Этилен С2Н4 0,487 0,360 0,880 0,453 1,120 0,585 1,440

промежуточных, инженерные формулы отсутствуют. В этом случае для оценки опасности можно воспользоваться подходом, предложенным в [6, 7]. Считается, что в результате зажигания возникает огненный шар, однако масса топлива, вовлеченного в горение, уменьшается с ростом времени задержки зажигания ti из-за разбавления воздухом части топлива, выброшенной в атмосферу до момента зажигания, на величину ДМ. Коэффициент участия топлива в горении огненного шара

X =

M 0 -AM Mn

получен в [7] в безразмерных переменных как функция параметра вещества времени задержки зажигания х = ti /tr (рассматривается только случай ti < tr) и безразмерного диаметра отверстия 5. Пока-

зано, что

X= 1

3/2

(9)

где а — коэффициент, равный 0,42 для выбросов при низком давлении, 0,59 — для выбросов при высоком давлении и, наконец, 0,87 — для истечений из резервуаров конечного объема.

Как видно, коэффициент участия топлива падает с уменьшением диаметра отверстия 5 вследствие более сильного разбавления газа. Нижняя граница для коэффициента участия может быть получена подстановкой минимального диаметра, соответствующего границе промежуточных и непрерывных истечений: при 5 = 5у для всех типов истечений получаем х(5у) = 1-0, откуда следует, что минимальный коэффициент участия топлива соответствует зажиганию в момент окончания истечения (х = 1): %тт = 0,33. Полученные результаты дают нижнюю границу для массы топлива, вовлеченного в диффузионное горение в режиме огненного шара, таким образом Мрв = (0,33^1)М0.

Примеры использования количественного критерия

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

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

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

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

истечения метана при низком давлении хранения составляют: 5. =0,155 и 5С = 0,509 (см. таблицу). Для большей конкретности проведем дальнейшие оценки для типичного газгольдера диаметром В0 = 25 м и высотой Н0 = 20 м (полный начальный объем У0 = 9800 м3, масса газа М0 = 6525 кг). Критический размер пробоины, отвечающий границе между струйными и промежуточными выбросами, составляет I. = 5.Г01/3 = 3,3 м, что дает общую площадь пробоины 8,64 м2 (1,75 % площади крыши). Диаметр отверстия, соответствующий переходу от промежуточных истечений к выбросам компактного облака, Вс = 5сГ01/3 = 10,9 м, в этом случае площадь повреждения равна 93 м2 (около 20 % общей площади крыши). Минимальная масса топлива, которая может участвовать в горении огненного шара, может быть оценена исходя из минимального коэффициента участия хтщ = 0,33, найденного из (9): мтп = X тт М 0 а 2000 кг. Соответственно, максимальный диаметр огненного шара и его время выгорания согласно эмпирическим формулам [10] (см. также [5])

Брв = 5,33М00,327; 1РВ = 0,92 М00,303 (10)

оцениваются как БРВ = 64^94 м, 1РВ = 9,2^13,2 с. Эти данные согласуются с характеристиками разрушения газгольдера и параметрами огненного шара, зафиксированными камерами видеонаблюдения в упоминавшемся инциденте 1993 г.

В качестве второго примера рассмотрим истечение сжатого метана при постоянном высоком давлении в резервуаре. Для конкретности принимаем внутреннее давление равным 50 бар, а диаметр отверстия I = 1 м. В данном случае цель состоит в оценке времен истечения, соответствующих различным режимам истечения. Массовый расход газа через выходное отверстие может быть получен на основе стандартных соотношений с учетом запирания потока [9]:

О = Си

пБ/

2 ) 2(у-1)

4 ^у + 1

у +1

уЯ <

1/2

Р 0.

е У

Полная выброшенная масса газа М0 = Gtr, а его объем при начальных условиях ¥0 = М0/р0. Приравнивая диаметр I каждому из критических диаметров Бс и I. и разрешая соответствующие уравнения Бс, . = 5 с, .Л,1/3(Ра/Р0)1/6 (см. (7)), найдем критические времена истечения, соответствующие границам между различными режимами выброса:

у+ 1

С . = ГУ + 1 ^ 2(У"1) I Р

г' = п с,

1/2/

,1/2

У Я 0Г,

5С3,..

0 у

Подставляя наружную температуру Т0 = 293 К и критические безразмерные диаметры отверстия

для метана (при высоком внутреннем давлении

5. = 0 ,195 и 5С = 0,647 согласно таблице) , получим

с

размерные значения критических времен tr а а 0,02(Р0/Р)1/2В, ^ а 0,5(Р0/Ра)1/2I (где I — в метрах, tr — в секундах). Для указанных выше значений диаметра и внутреннего давления получим ^ а 0,14 си tJr а 3,5 с. Эти времена весьма малы, следовательно, если после открытия в резервуаре поддерживается высокое давление, в течение нескольких первых секунд устанавливается струйное истечение газа, и струйный выброс может рассматриваться в качестве наиболее вероятного исхода разгерметизации. Только при истечениях, длящихся менее секунды, можно ожидать образования в атмосфере переобогащенного топливного облака.

Наконец, оценим возможные последствия истечения из резервуара конечного объема. Рассмотрим сосуд объемом 200 м3 со сжатым природным газом при двух начальных давлениях 10 и 100 бар, при этом масса газа составит соответственно 1350 и 13500 кг. Безразмерные критические диаметры при истечении из резервуара конечного объема при высоком начальном давлении для метана 5. = 0,252 и 5С = 0,834 (см. таблицу), что приводит к следующим размерным величинам: при меньшем начальном давлении Б. =5.А01/3(Р^Р0)1/12 а 1,17 м, Вс =5СГ01/3(Ра/Р0)1/12 а 3,9 м, тогда как при большем — I. = 1,0 м, Вс = 3,2 м. Сопоставив полученные данные, можно увидеть, что критические диаметры довольно слабо зависят от давления хранения. Кроме того, очевидно, что истинно "мгновенным" истечение можно считать только в том

случае, если размер отверстия сравним с характер-

1/3

ным линейным размером сосуда Г0 = 5,8 м. Приведем оценки параметров огненных шаров (см. формулу (10)), которые могут возникнуть при разгерметизации резервуара: при меньшем внутреннем давлении МРВ = 420^1350 кг, БРВ = 38^56 м, tрв = = 5,7^8,0 с; при большем — МРВ = 4200^13500 кг, БРВ = 82^120 м, tрв = 11^16 с.

Заключение

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

т

е

I

38

0869-7493 ООЖАРООЗРЫЮОБЕЗООАСООСТЬ 2010 ТОМ 19 №2

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

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

СПИСОК ЛИТЕРАТУРЫ

1. AlChE/CCPS. Guidelines for Chemical Process Quantitative Risk Assessment. — New York: Amer. Inst. Chem. Engineers, 1989.

2. AlChE/CCPS. Guidelines for Evaluating the Characteristics of Vapor Cloud Explosions, Flash Fires, and BLEVEs. — NewYork : Amer. Inst. Chem. Engineers, 1995.

3. TNO. Methods for the Calculation of Physical Effects Resulting from Releases of Hazardous Materials (Liquids and Gases) / TNO "Yellow Book". 2nd Ed. — Voorburg, TNO, 1992.

4. ГОСТ P 12.3.047-98. ССБТ. Пожарная безопасность технологических процессов. Общие требования. Методы контроля. — Введ. 2000-01-01. — М. : Изд-во стандартов, 1998.

5. Методика определения расчетных величин пожарного риска на производственных объектах // Пожарная безопасность. — 2009. — № 3.

6. Makhviladze, G. М. Burning regimes for the finite-duration releases of fuel gases / G. M. Makhvi-ladze, J. P. Roberts, S. E. Yakush // Twenty Sixth Int. Symp. on Combustion. — The Combustion Institute, Pittsburgh, PA, USA. — 1996. — Vol. 1. — P. 1549-1555.

7. Махвиладзе, Г. M. Образование и горение газовых облаков при аварийных выбросах в атмосферу / Г. М. Махвиладзе, Дж. П. Роберте, С. Е. Якуш// Физика горения и взрыва. — 1997. — Т. 33, № 2. — С. 23-38.

8. Makhviladze, G. М. Large-scale unconfined fires and explosions. Invited Paper for a special memorial session devoted to the September 11th fire / G. M. Makhviladze, S. E. Yakush // 29th Symposium (Int.) on Combustion. — Proceedings of the Combustion Institute. — 2003. — Vol. 29. — P.195-210.

9. Лойцянский, Л. Г. Механика жидкости и газа / Л. Г. Лойцянекий. — М.: Наука, 1987. — 840 с.

10. Moorhouse, J. Thermal radiation hazards from large pool fires and fireballs : a literature review /

J. Moorhouse, M. J. Pritchard // IChemE Symp. Series. — 1982. — No. 71. — P. 397-428.

Материал поступил в редакцию 24 ноября 2009 г.

©Якуш С. Е., Махвиладзе Г. М, 2010 г. (e-mail: yakush@ipmnet.ru, gmakhviladze@uclan.ac.uk).

А.Н. ЧЛЕНОВ, ГЛ. БУЦЫНСНАЯ, И,Г. ДР0ВНИК08А. -Ч.1.- 316 с.

В,П. БАБУРОВ, В.В, БАБУРИН, В.И. ФОМИН. — Ч, 2. - 300 С,

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

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

р*- 4JJ и I -ЧЩ

СИСТЕМ ДАННОЙ И ДОЖАРНОЙ СИГНАЛИЗАЦИИ

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

WEB-САЙТ: www.flrepress.ru

ЭЛ. ПОЧТА:

mai l@fi repress, ru; izdat pozhnauka @ mail.ru

Телефон: (495) 228-09-03,

тел./факс: (495)445-42-34