CC BY
49А.В. Завгороднев, к.т.н., генеральный директор ООО «Газпром трансгаз Ставрополь»;
Г.С. Акопова, к.т.н., начальник лаборатории охраны окружающей среды и ресурсосбережения,
ООО «Газпром ВНИИГАЗ»; А.В. Мельников, заместитель начальника отдела хранения СПГ и газового конденсата, ООО «Газпром добыча шельф»
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ РАССЕИВАНИЯ В АТМОСФЕРЕ ОРГАНИЗОВАННЫХ НЕСТАЦИОНАРНЫХ ВЫБРОСОВ ГАЗА НА ОБЪЕКТАХ ГАЗОТРАНСПОРТНЫХ ПРЕДПРИЯТИЙ
Статья содержит анализ действующих моделей и методик по рассеиванию загрязняющих веществ в атмосфере и возможности их адаптации для расчетов загрязнения атмосферы природным газом от залповых выбросов для газотранспортных предприятий при нестационарных условиях. Приведены выводы о целесообразности применения методики PHAST DNV для расчета нестационарных выбросов.
Процессы переноса загрязняющих примесей от источников естественного и антропогенного происхождения являются существенным фактором, влияющим на экономическое состояние промышленных и агропромышленных регионов. Для их описания и прогнозирования требуются математические модели и методики, адекватно отражающие рассеивание загрязняющих веществ в атмосфере.
Характерной особенностью функционирования газотранспортных предприятий являются организованные непродолжительные «залповые» выбросы в атмосферу различной мощности в результате регламентных технологических операций. При этом выбросы природного газа поступают в атмосферу при стравливании через свечи со скоростями, близкими к звуковым.
В настоящее время при проведении расчета рассеивания выбросов загрязняющих веществ в атмосферном воздухе широко используется нормативный документ ОНД-86 [1]. Методика утверждена на государственном уровне и рекомендована для использования
при разработке природоохранных мероприятий и осуществлении государственного контроля охраны воздуха. Существенными недостатками методики ОНД-86 являются неучет скорости истечения газа и атмосферной диффузии, приводящие к тому, что по результатам расчетов газовое облако имеет концентрацию газа существенно выше предельно допустимой. Последнее обстоятельство не подтверждается результатами инструментального контроля качества атмосферного воздуха. Это приводит к увеличению платы за его выброс в 5-кратном размере (по временно согласованному выбросу - ВСВ).
За рубежом получили распространение ряд методик, входящих в программные пакеты: ANSYS CFX, FLOTRAN, РИАБТ DNV, ТИО и другие, позволяющие рассчитать нестационарные выбросы. К недостаткам данных методик относятся большая стоимость программного продукта (от 30 до 150 тыс. евро), трудность задания начальных и граничных условий, длительность расчета и необходимость наличия сертифицированного персонала для работы с ней.
В России и за рубежом известны модели, позволяющие прогнозировать загрязнение атмосферы при разрывах газо- и конденсатопроводов. На их основе созданы методики [2-16], которые используются при проектировании трубопроводных систем. Так, методика (5) применима для условий сравнительно ровной и малопересеченной местности. Имеется возможность учета двух состояний устойчивости атмосферы - устойчивого и неустойчивого. Методика носит частный характер и может использоваться лишь при стационарном выбросе (истечение через малое отверстие в трубе) для некоторых метеоусловий.
Большинство известных методик и программных комплексов предназначено для решения узкого класса задач.
В рамках настоящей статьи рассматривается возможность применения для расчетов загрязнения атмосферного воздуха природным газом от «залповых» выбросов для предприятий газовой отрасли методики, разработанной компанией DNV, и реализующего эту методику программного комплекса РИДБТ .
Одной из наиболее важных метеорологических характеристик атмосферы является степень ее устойчивости, характеризующая «рассеивающую способность» атмосферы.
В литературе нет общепринятого критерия определения категорий стабильности, однако большинство исследователей пользуются простейшей классификацией Паскуилла по вертикальному градиенту температур (табл. 1).
Поскольку состояние устойчивости атмосферы, по сути, определяется интенсивностью вертикальных конвективных течений, оно может существенно измениться и в течение суток. Характерное распределение стабильности атмосферы в течение дня для широт объектов расположения ООО «Газпром трансгаз Ставрополь» представлено на рисунке 1 для различных времен года. Как видно, распределение состояний устойчивости атмосферы в холодное и теплое время года весьма отличаются. Так, в холодное время года доминируют нейтральное и стабильное состояние атмосферы, а в теплое время года: ночью - стабильное, а днем - неустойчивое.
Для описания зависимости реализации того или иного класса устойчивости атмосферы от скорости ветра и были проанализированы данные станций наблюдений в США и в РФ (Обнинский ИЭМ). Результаты представлены в таблице 2 и на рисунке 2. Определенное их различие свидетельствует о необходимости использования только местной метеоинформации, строго соответствующей определенному типу земной поверхности (равнинная, холмистая и т.д.). Влияние состояния земной поверхностью проявляется также и в эффекте торможения воздушных масс, который характеризуется, как правило, так называемым параметром шероховатости (Д0) в зависимости от структуры ландшафта и
Таблица 1. Классификация устойчивости атмосферы
Вертикальный градиент температуры (ДT/Дz) Класс устойчивости
<-1,9 А — сильная конвекция
— (1,9 1,7) В — конвекция
-(1,7.. .1,5) С — умеренная конвекция
— (1,5....0,5) 0 — нейтральная
—0,5...+1,5 Е — инверсия
+(1,5...4,0) и более F — сильная инверсия
Вероятность,% 100%
90%
80%
70%
60%
50%
40%
30%
20%
10%
0%
11 13 15 17 19 21 23
Час суток ноябрь-апрель
Вероятность, %
■ класс Е.Р
□ класс С,Б
□ класс А,В
9 11 13 15
Час суток
май-октябрь
Рис. 1. Распределение частоты реализации различных классов устойчивости атмосферы в зависимости от времени суток
сезона в разрезе года. Наиболее характерные значения Д0 даны в таблице 3. Таким образом, для выбора параметров модели рассеивания в условиях
конкретного региона расположения объекта следует принимать во внимание характерные особенности этого региона.
Таблица 2. Значения класса устойчивости в зависимости от скорости ветра по сезонам года
и„ м/с Апрель - сентябрь Октябрь - март
IV V VI VII VIII IX X XI XII I II III
А В С 0 Е F А В С 0 Е F
0—1 1,4 7,0 29,8 23,7 4,0 34,1 0,0 0,1 5,3 17,2 20,8 56,6
2—3 3,2 10,5 36,5 28,7 3,8 17,3 0,0 0,0 5,8 23,0 26,0 45,2
4—5 2,2 12,4 39,6 31,0 3,7 11,1 0,0 0,1 6,4 26,0 28,1 39,4
6—7 4,3 11,9 39,4 31,5 3,4 9,5 0,1 0,2 6,7 27,3 29,0 36,7
8—10 4,3 12,2 40,4 32,8 3,7 6,6 0,0 0,2 5,9 26,0 27,2 49,7
Рис. 2. Распределение классов устойчивости в зависимости от скорости ветра (А. Макнаут, В. Бызова)
Таблица 3. Значения параметра шероховатости для различных типов ландшафта
Ландшафт Сезон До, м
Городской все сезоны 1
Сельскохозяйственный весна 0,03
лето 0,25
зима 0,001
Пастбище весна 0,02
лето, осень 0,05
зима 0,001
Лиственные и хвойные леса все сезоны 1,0
Болото весна 0,1
лето 0,15
осень 0,1
зима 0,001
Смесь сельскохозяйственных культур с пастбищами весна 0,03
лето 0,1
осень 0,08
зима 0,001
Таблица 4. Значения показателя р от параметра шероховатости и стабильности атмосферы
Категория стабильности атмосферы Параметр шероховатости Д0, м
0,01 0,1 1 3
А 0,05 0,08 0,17 0,27
В 0,06 0,09 0,17 0,28
С 0,06 0,11 0,20 0,31
0 0,12 0,16 0,27 0,37
Е 0,34 0,32 0,38 0,47
F 0,53 0,54 0,61 0,69
Известно, что скорость ветра существенно изменяется с высотой. В прикладных исследованиях наиболее часто используется степенная зависимость вида и^)=110^Д0)р, где и0 - скорость на «стандартной» высоте z0 (обычно z0 =10 м). Значения показателя р также зависят от класса устойчивости атмосферы и шероховатости поверхности Д0 (табл. 4).
МЕТОДИКИ РАСЧЕТА КОНЦЕНТРАЦИЙ НЕЙТРАЛЬНОЙ ПРИМЕСИ В АТМОСФЕРНОМ ВОЗДУХЕ
Среди относительно простых методик расчета полей рассивания следует выделить две, нашедшие широкое применение на практике ввиду своей сравнительной простоты - это методика ОНД-86 и Гауссова модель.
ОНД-86 является основным документом для расчета концентраций вредных веществ в атмосферном воздухе в России. Методика предназначена для определения приземных концентраций (в двухметровом слое над поверхностью земли), а также вертикального распределения концентраций для стационарных источников. Методика не распространяется на расчет концентраций на дальних расстояниях (более 100 км от источника выброса). ОНД-86 нельзя использовать для сценарных и прогнозных расчетов при конкретных метеоусловиях и специфических источников выбросов, характерных для объектов ОАО «Газпром». Он может использоваться лишь для оценки максимально возможной концентрации газовоздушной примеси при наихудших условиях рассеивания.
Наибольшая концентрация загрязняющего вещества в воздухе См для заданного источника выбросов (высота, скорость и температура выбросов и другие параметры которого зафиксированы) может достигаться только при одной характерной для данного источника выбросов скорости ветра им (при скоростях ветра, меньших им или больших им, максимальное значение концентрации на уровне земли оказывается меньшим, чем См). Более того, концентрация См, при ветре им может достигаться только на определенном удалении от источника хм (на меньших и на больших расстояниях от источника, чем хм ,при опасной скорости ветра концентрация будет меньше, чем См).
Рассмотрим кратко основные положения модели. Расчет пространственного распределения концентраций базируется на расчете трех основных параметров: См, им и хм. Максимальное значение приземной концентрации вредного вещества См (мг/м3) при выбросе газовоздушной смеси из одиночного точечного источника с круглым устьем достигается при неблагоприятных метеорологических условиях на расстоянии хм (м) от источника и определяется по формуле:
_ АМРтт1
где А - коэффициент, зависящий от температурной стратификации атмосферы (коэффициент определяется географическим местоположением места, и значение задается для различных регионов: Средняя Азия, Бурятия, Читинская область, европейская территория южнее 50° с. ш., Кавказ, Молдавия; Дальний Восток, европейская территория от 50 до 52° с. ш. и др.);
М (г/с) - масса вредного вещества, выбрасываемого в атмосферу в единицу времени;
F - безразмерный коэффициент, учитывающий скорость оседания вредных веществ в атмосферном воздухе; т и п - коэффициенты, учитывающие условия выхода газовоздушной смеси из устья источника выброса; значения этих коэффициентов определяются в зависимости от ряда параметров, характеризующих источник выброса;
Н (м) - высота источника выброса над уровнем земли;
ц - безразмерный коэффициент, учитывающий влияние рельефа местности (в случае ровной или слабопересеченной местности с перепадом высот, не превышающим 50 м на 1 км, г) = 1);
ДТ (°С) - разность между температурой выбрасываемой газовоздушной смеси Тг и температурой окружающего атмосферного воздуха Тв;
V1 (м3/с) - расход газовоздушной смеси. Значение безразмерного коэффициента F для газообразных вредных веществ принимается равным 1.
Расстояние хм (м) от источника выбросов, на котором приземная концентрация с (мг/м3) при неблагоприятных метеорологических условиях достигает максимального значения см, определяется по формуле:
ли
Х„=^(1Н ,
где безразмерный коэффициент d находится по весьма громоздким формулам и здесь не приводится / 1/. Максимальное значение приземной концентрации вредного вещества сми (мг/м3) при неблагоприятных метеорологических условиях и скорости ветра и (м/с), отличающейся от «опасной» скорости ветра им (м/с), определяется по формуле с = ГС
'-ми 1'-м '
где г - безразмерная величина, определяемая по формулам через величины и и им.
Расстояние от источника выброса хми (м), на котором при скорости ветра и и неблагоприятных метеорологических условиях приземная концентрация вредных веществ достигает максимального значения сми (мг/м3), определяется по формуле: хми = РХм ,
где р - безразмерный коэффициент, определяемый величиной отношения
и/им.
При «опасной» скорости ветра им приземная концентрация вредных веществ с (мг/м3) в атмосфере по оси факела выброса на различных расстояниях х (м) от источника выброса определяется по формуле: с = с(х) = Sl(x)Смu
где - безразмерный коэффициент, определяемый в зависимости от отношения х/хми и коэффициента Р. Значение приземной концентрации вредных веществ в атмосфере су (мг/м3) на расстоянии у (м) по перпендикуляру к оси факела выброса определяется по формуле:
су = су(х,у) = s2(y/x,u)c,
где зависимость s2 от х/у, и задана.
РАССМОТРИМ ДАЛЕЕ ОСНОВНЫЕ ПОЛОЖЕНИЯ МОДЕЛИ ГАУССА
При выбросе из трубы высотой И газовоздушная смесь имеет некоторую вертикальную скорость истечения из трубы. Кроме того, выбрасываемая газовоздушная смесь, как правило, имеет более высокую температуру, чем окружающий воздух. Это приводит к тому, что струя выбрасываемого газа поднимается над источником на некоторую дополнительную высоту и даль-
ше распространяется так, как если бы оно изначально было выпущено на этой высоте. Чтобы учесть это, в качестве высоты источника берется не реальная высота источника И, а так называемая эффективная высота Не№
Пусть задана система координат, начало которой совпадает с проекцией источника на поверхность земли,ось Ъ направлена вертикально вверх, ось X лежит в горизонтальной плоскости и направлена по ветру, ось Y лежит в горизонтальной плоскости и направлена перпендикулярно направлению ветра. В этой системе координат концентрация загрязняющего вещества, выбрасываемого непрерывно действующим источником в точке с координатами х, у, z, определяется соотношением:
С(Х,У,2)=
2.1ы5ехр('2?;)(иф<
(г-Не„):
2 о2,
-)+
где
оу = оу (х), о, = о, (х) - функции, характеризующие расширение факела в перпендикулярных к ветру направлениях при удалении от источника выброса: сгу, - в горизонтальном направлении, с, - в вертикальном направлении; обе функции зависят только от удаления х;
НеН = ($, т) - эффективная высо-
та источника ($ - задает температуру газовоздушной смеси, расход, скорость истечения, высоту устья трубы, т - задает температуру атмосферного воздуха и скорость ветра на высоте устья источника);
и„ = и„ (т, Н^) - скорость ветра на высоте 10 м, при Н^ « 10 м; либо на высоте НеИ при Н^ > 10 м;
S - мощность выброса загрязняющего вещества.
В методике ОНД-86 введено избыточное число параметров (четыре), тогда как достаточно двух. Избыточность параметров усложняет расчет, а иногда из-за противоречивых условий приводит к ошибкам. Некоторые формулы методики ОНД-86 имеют недостаточную точность.
В методике ОНД-86 отсутствует четкое определение проектной и проверочной задач установления предельно допустимых выбросов (ПДВ), вследствие чего задача ошибочно сводится к расчету
Таблица 5. Значения коэффициентов
Z0 = 0,01 Zo = 0,1 N о = 1, 0
Класс устойчивости а Ь с d с d с d
А 0.527 0.856 0.154 0.940 0.28 0.90 0.615 0.830
В 0.371 0.866 0.133 0.890 0.23 0.85 0.539 0.770
С 0.209 0.897 0.121 0.850 0.22 0.80 0.533 0.720
D 0.128 0.905 0.108 0.810 0.20 0.76 0.456 0.680
Е 0.098 0.902 0.078 0.780 0.15 0.73 0.348 0.650
F 0.065 0.902 0.062 0.720 0.12 0.67 0.309 0.580
концентрации загрязняющего вещества в приземном слое атмосферы. Фактически же необходимо определять расход атмосферного воздуха, участвующего в рассеивании выброса из данного источника, по отношению заданной точки при определенном направлении и скорости ветра.
Указанные ошибки и другие неточности методики ОНД-86 делают сомнительными результаты расчетов, на основе которых разрабатываются проекты ПДВ.
Модель Гаусса базируется на следующих предположениях:
• метеоусловия и характеристики поверхности, определяющие распространение и рассеивание облака на всей территории, одинаковы и постоянны по времени;
• источник выбросов действует в течение длительного времени;
• вследствие влияния турбулентной диффузии концентрация загрязняющего вещества распределяется в горизонтальном и вертикальном направлениях по закону Гаусса, с дисперсиями су, а2, зависящими от расстояния вниз по ветру от источника, состояния устойчивости атмосферы и характера ландшафта территории;
• выбрасываемое источником в атмосферу загрязняющее вещество консервативно;
• поглощение загрязняющего вещества на поверхности земли отсутствует. Модель Гаусса работает при скоростях ветра, превышающих 1 м/с.
Величины дисперсий гауссовского распределения концентрации су,, о2, входящие в формулу расчета разовых концентраций, описываются следующими эмпирическими зависимостями:
ау(т)=ахь
а2(т)=схС,
где х - удаление от источника по ветру.
Значения коэффициентов а, с, Ь, С при различных значениях коэффициента шероховатости Z0 в формулах определяются по таблице 5.
Коэффициент шероховатости поверхности Z0 здесь измеряется в метрах. Для расчета скорости ветра на заданной высоте Н используется зависимость:
С(ХУ2) =_____-____ехрГ-^^2)
Ц ,У, ; 2я3/2с,аца7 2о\}
2о2,
)+
У(Н)=У1*
ід№)
ч>
где V! - скорость ветра, измеренная на высоте флюгера Н^ Н = 10 м,
Z0 - коэффициент шероховатости поверхности для данной территории. Эффективная высота источника определяется по формуле
Н^ = Ьі+ДЬі
При выбросе из труб величина ДЬі определяется по формуле
где
Уд - скорость истечения газовоздушной смеси из сопла источника, [м/с],
Тд - температура выбрасываемой газовоздушной смеси, [град С°],
Та - температура атмосферного воздуха [град С°],
d - диаметр сопла трубы источника [м],
Ьі - реальная высота трубы источника [м],
Ук - скорость ветра на высоте Ьі [м/с]. В случае мгновенного залпового выброса ЗВ веществ массы М мгновенные значения концентрации (период усреднения - 3-5 мин.) могут быть определены на основе нестационарной модели Гаусса:
6Хр(- 2о2 )(еХр(
+еХр(-(22о2Г)2))
где и - скорость ветра на высоте Не№ t - время, ох = ох (х) - функция, характеризующая расплывание облака вдоль направления ветра при удалении от источника выброса.
Известным фактом является значительное расхождение между собой результатов расчетов концентраций по Гауссовой модели и методике ОНД-86, в частности:
• ширина подфакельной зоны по фронту. Методика ОНД-86 дает более широкий по фронту факел, чем методика Гаусса (даже при неустойчивом состоянии атмосферы, соответствующем наиболее активному ее перемешиванию);
• глубина подфакельной зоны по ветру. При устойчивом состоянии атмосферы Гауссова методика дает большие значения глубины зоны повышенных концентраций по ветру, чем методика ОНД-86, а при неустойчивом состоянии атмосферы - меньшие. В отдельных случаях (низкие холодные источники высотой 1-5 м, при выбросе сильно перегретой ГВС, при выбросе ГВС через сопло трубы со скоростями более 20 м/ сек. и др.) указанное соотношение глубины зон повышенных концентраций нарушается, но в целом оно является верным.
Существующие в настоящее время зарубежные и российские методики по расчету рассеивания загрязняющих веществ в атмосфере:
• имеют большое количество ограничений либо ориентированы на расчет рассеивания определенных загрязнителей;
• методики, применяемые для расчетов параметров нестационарных источников, представляют собой сложные системы дифференциальных уравнений в частных производных, которые требуют задания граничных условий и большого комплекса параметров внешней среды. Реализуются они в виде сложных комплексов программ, что требует высококвалифицированных специалистов. Существенным недостатком этих методик является невозможность задания на практике всего комплекса граничных и начальных условий и существенные временные затраты для расчета даже одного варианта, который соответствует одному состоянию атмосферы (конкретные значения скорости и направления ветра, стабильность атмосферы).
Известны различные экспериментальные установки по изучению стационарных процессов, позволяющие установить влияние отдельных факторов на рассеивание вредных примесей, а также
исследовать поведение пассивной примеси в потоках турбулентности различной интенсивности от точечного и линейного источника, получить высоту начального подъема струи, расположение точки максимума концентрации и т.д.
С целью изучения процесса рассеивания струй при давлениях до 1 МПа для случая залповых выбросов газа, имеющих место на газотранспортных системах, и для обоснования возможности применения наиболее продвинутой методики DNV PHAST выполнены специальные исследования на экспериментальной установке, состоящей из горизонтальной аэродинамической трубы с трактами подвода воздуха и газа. Проведены также натурные исследования по определению фактического рассеивания природного газа (метана) при продувках и стравливании из технологического оборудования объектов ООО «Газпром трансгаз Ставрополь» в полевых условиях.
Таким образом, для выбора параметров модели рассеивания в условиях конкретного региона расположения объекта следует принимать во внимание характерные особенности этого региона.
ВЫВОДЫ
Приведена сравнительная оценка существующих в России и за рубежом моделей и методик, описывающих нестационарный процесс, и показано, что они имеют массу ограничений, ориентированы на расчет рассеивания определенных загрязнителей.
С целью получения достоверных данных при рассеивании природного газа в атмосфере от залповых выбросов, при продувках и стравливании газа из технологического оборудования целесообразно использовать методику, разработанную компанией DNV, и реализующий ее программный комплекс РНАБТ.
Литература:
1. ОНД-86 «Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий».
2. Пененко В.В., Цветова Е.А. Моделирование процессов переноса в прямых и обратных задачах климатоэкологического мониторинга и прогнозирования. Оптика атмосферы и океана, 1999,12, № 6, с. 482.
3. Шварц К.Г., Шкляев В.А. Моделирование процессов переноса примеси в свободной атмосфере. Метеорология и гидрология, 2000, № 8, с. 44-54.
4. Методика прогнозирования масштабов заражения сильнодействующими ядовитыми веществами при авариях (разрушениях) на химически опасных объектах и транспорте.
5. Методика расчета концентрации вредных веществ в атмосфере промышленных площадок от низких источников с учетом особенностей объемно-планировочных решений.
6. Practical Guide to Atmospheric Dispersion Modeling// Trimity Consultants, Inc., Dallas, Texas. 1990, с. 75-80.
7. Texas Air Control Board Quality Modelling Procedures// Texas Air Control Board, Austin, Texas, 1988.
8. Practical Guide to Atmospheric Dispersion Moiling// Trinity Consultans, Inc., Dallas, Texas, 1988.
9. Guideline on Air Quality Models (Revised) // U.S. Environmental Protection Agency, EPA-450/2-78-027R, 1986.
10. Клименко Е.Т., Максименко А.Ф. (РГУНГ), Стативко В.Л. (ОАО «Газпром»). Метод расчетной оценки приземных концентраций выбросов загрязняющих веществ от нестационарных источников. Научно-технический сборник Проблемы газовой промышленности. ИРЦ Газпром, 2000, № 1. Раздел: Охрана воздушного бассейна.
11. Акопова Г.С., Ганага С.В., Толстова Н.С. Проблемы моделирования рассеивания залповых выбросов загрязняющих веществ (природного газа) в атмосфере // Трубопроводный транспорт. Теория и практика. 2010. №1 (17).
12. Онищенко Г.Г., Новиков С.М., Рахманин Ю.А., Авалиани С.Л., Буштуева К.А. Основы оценки риска для здоровья населения при воздействии химических веществ, загрязняющих окружающую среду / Под. ред. Ю.А. Рахманина, Г. Г. Онищенко. -М. НИИ ЭЧ и ГОС, 2002. с. 168-171.
13. Филатова Е.Н. Моделирование загрязнения атмосферы по оперативным данным. Дис. канд. физ.-мат. наук. - Санкт-Петербург, 2005. с. 100.
14. Гриценко А.И., Максимов В.М., Самсонов Р.О., Акопова Г.С. Экология: нефть и газ. - М. ИКЦ «Академкнига», 2009.
15. Максимов В.М., Розенберг Г.Д., Исаев В.И., Лимар Е.Е. Теоретические основы прогнозирования последствий аварий на газо- и конденсатопроводах// Защита от коррозии и охрана окружающей среды. 1995. № 5. с. 8-15.
16. Максимов В.М., Розенберг Г.Д., Исаев В.И. и др. Методические рекомендации по расчету параметров выброса газовой смеси и ее рассеивания в атмосфере при аварийных разрывах газопроводов. - М. ВНИИГаз, 1992. с. 64.
Ключевые слова: залповый выброс, нестационарный процесс, продувка, модели, рассеивание, атмосфера.
