Образование и распространение паровоздушных облаков сжиженных углеводородных газов при техногенных авариях тема диссертации и автореферата по ВАК РФ 05.26.03, кандидат технических наук Исламхузин, Дамир Ягфарович

• Исламхузин, Дамир Ягфарович
• кандидат технических наук
• 2003

Оглавление диссертации кандидат технических наук Исламхузин, Дамир Ягфарович

Введение.

Глава 1. Анализ причин и последствий возникновения и развития аварийных ситуаций при транспортировке СУГов.

1.1. Обзор типичных аварийных ситуаций

1.2. Теория атмосферной турбулентности. 11 1.2.1. Уравнения баланса кинетической энергии турбулентности и дисперсии температуры. 11 1.2.3. Теория Монина - Обухова для приземного слоя.

1.3. Современное состояние теорий, описывающих распространение вредных веществ, попадающих в атмосферу в условиях аварийных разгерметизации. Нормативные методики.

1.3.1. Модель диффузии примесей в атмосфере.

1.3.2. Граничные условия.

1.3.3. Коэффициенты уравнения атмосферной диффузии.

1.3.4. Обзор решений уравнения (1.9). Нормативные методики.

1.4. Постановка задачи исследования. '

Глава 2.Теоретическое исследование процесса истечения СУГов из аварийных отверстий характерных размеров.

2.1. Моделирование закономерностей формирования профиля давлений по длине продуктопровода при перекачке СУГов.

2.4. Полная схема расчета развития аварийной ситуации.

Глава 3. Математическое моделирование распространения ПВО в атмосфере.

3.1. Выбор и формулировка исходной системы уравнений.

3.2. Выбор приемов решения выбранного математического описания. 44 3.2.1. Обзор существующих методов решения.

3.3. Расчетная схема.

3.3.1. Матричная запись исходной системы уравнений.

3.3.2. Схема расчета.

3.3.3. Аппроксимация пространственных производных.

3.4. Коэффициенты турбулентной диффузии.

3.5. Описание численного решения математической модели.

3.5.1. Граничные и начальные условия.

3.5.2. Конечноразностный аналог математической модели изменения скорости ветра по высоте.

3.5.3. Конечноразностный аналог математической модели гидродинамического потока тяжелого газа.

3.6. Проверка адекватности математической модели.

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

4.1. Утечка ШФЛУ через свищ диаметром 5 мм.

4.1.1. Расчет размеров и интенсивности испарения.

4.1.2. Расчет пространственно-временного распределения концентрации взрывоопасного паровоздушного облака.

4.2. Утечка ШФЛУ через свищ диаметром 50 мм.

4.2.1. Расчет размеров и интенсивности испарения.

4.2.2. Расчет пространственно-временного распределения концентрации взрывоопасного паровоздушного облака.

4.3. Утечка ШФЛУ через свищ диаметром 100 мм. 95 4.3.1. Расчет размеров и интенсивности испарения. 95 Выводы. 96 Исследование температурной стратификации атмосферы на характеристики турбулентности. 98 Заключение. ЮЗ Литература. Ю<>

Введение диссертации (часть автореферата) на тему «Образование и распространение паровоздушных облаков сжиженных углеводородных газов при техногенных авариях»

Любое предприятие химической промышленности в той или иной мере представляет потенциальную опасность для экологии и населения. Особую опасность представляют предприятия, осуществляющие переработку, хранение и транспортировку нефтехимических продуктов, обладающих взрывоопасными и/или токсическими характеристиками. Особое место в списке потенциально опасных веществ занимают вещества с низкими температурами кипения. К ним относятся сжиженные углеводородные газы (этан, пропан, бутан), хлор и его соединения, оксид этилена и многие другие вещества. Процессы переработки этих веществ, а также их хранение и транспортировка проводятся под давлением (0,6-^10 МПа) и, как правило, в жидком состоянии. В основном транспортировка сжиженных газов на большие расстояния осуществляется по подземным магистральным трубопроводам, и, в случае повреждения трубопровода, возникает аварийная ситуация, отличная от аварий при разрушении оборудования на технологических установках. Во время таких аварий возникает опасность выброса большого количества сжиженного газа в окружающую среду. Основной причиной этого является большая протяженность магистралей, вследствие чего, бывает очень трудно обнаружить утечку опасного вещества. При этом утечка может происходить в опасной близости от населенных пунктов, авто-, железнодорожных магистралей, сель-хоз. угодий и т.д.

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

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

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

Для оценки последствий вероятных техногенных аварий на опасных производственных объектах с целью принятия компенсационных мер (выбора площадки застройки и т. д.), оценки риска от эксплуатации опасных производственных объектов, например, при разработке декларации промышленной безопасности согласно федеральному закону ФЗ 116 [1], необходима количественная оценка распространения паровоздушных облаков в пространстве. Декларация промышленной безопасности уточняется или разрабатывается в случае изменения сведений, содержащихся в декларации промышленной безопасности, или в случае изменения требований промышленной безопасности [2].

Приемов для полноценной количественной оценки риска от эксплуатации магистральных продуктопроводов недостаточно. Наиболее применяемые на сегодняшний день методики расчетов истечения и распространения СУГов в атмосфере дают противоречивые результаты. Это обусловлено тем, что процессы истечения СУГов из емкостей и трубопроводов, а также их рассеивания в нижних слоях атмосферы протекают в условиях, которые сильно зависят от множества природных факторов. Например, образование и рассеивание ПВО в атмосфере зависит от температуры окружающей среды, которая изменяется в пределах от-30°С до +30°С, от скорости ветра, от степени турбулизованности атмосферы и т.д. Для учета влияния всех этих факторов необходимо решать точные уравнения гидродинамики, записанные в дифференциальной форме. Существующие методики, как правило, основаны на определенных допущениях, которые игнорируют влияние тех или иных факторов.

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

Заключение диссертации по теме «Пожарная и промышленная безопасность (по отраслям)», Исламхузин, Дамир Ягфарович

Выводы.

Как видно из рис. 4.1-4.8, наиболее высокие концентрации углеводородов создаются непосредственно около псевдоледяной каверны, из которой испаряется ШФЛУ. На рис. 4.1, 4.5 показаны распределения концентраций углеводородов вдоль оси факела выброса при скоростях ветра на высоте 10 миа = 8 иЗ м/с для диаметров аварийных свищей 5 мм и 50 мм. На рис. 4.2, 4.6 изображены распределения концентрации поперек оси факела выброса.

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

Эти результаты позволяют дать оценку опасности образующегося при утечке ШФЛУ паровоздушного облака. Очевидно, наиболее опасным случаем является превышение локальной концентрации ШФЛУ в факеле выброса нижнего предела взрываемости газовоздушной смеси. Последняя имеет порядок 2.4 % объемных для пропана [71], что соответствует примерно 48 г/м3 при переходе к относительным массовым концентрациям.

Данные численного анализа показывают, что эта концентрация при диаметре коррозионного отверстия, равном 5 мм достигается на расстоянии 61м от источника на оси факела выброса при скорости ветра на высоте 10 м, равной 8 м/с, и 108 м при скорости ветра на высоте 10 м, равной 3 м/с. Т.е. распределение концентраций возрастает с уменьшением скорости ветра. Аналогичное явление наблюдается и при диаметре свища 50 мм. Это объясняется тем, что в приземном слое атмосферы турбулентность главным образом генерируется за счет градиента скорости, следовательно, с ее уменьшением падают значения турбулентных пульсаций, что приводит к снижению значений коэффициентов турбулентной диффузии. При этом происходит более медленное рассеивание облака и сохранение больших концентраций на значительных расстояниях от источника испарения.

Исследование температурной стратификации атмосферы на характеристики турбулентности.

U , м/с а'

Рис. 4.9. Влияние атмосферной стратификации на формирование профиля скорости ветра по высоте.

Z - высота над уровнем земли; иа - скорость ветра; L - характерный масштаб Обухова. t, °c i

-100

-150

-200

-50 0 50 100 150 200 250 300 350 400 450 500

X, м

Рис. 4.10. Температура на оси облака на разных высотах над уровнем земли в зависимости от продольной координаты. t - температура облака; X - продольная координата; Z - высота над уровнем земли.

Z = / 3,7 м •

1,4 M и--

-r iiz= 0,5 m •

----- i 1 z = ^- 0,2 м ----г 0 м. -1100

-600

-100

400

900 L, M

Рис. 4.11. Зависимость динамической скорости от атмосферной стратификации, и* динамическая скорость в турбулентном потоке; L — характерный масштаб Обухова.

Z, м

L = 0 м

20

L = 228,84 м L = 482,39 м

Kz,m7C

Рис. 4.12. Зависимость коэффициента турбулентной диффузии атмосферного воздуха от высоты при различных значениях атмосферной стратификации.

Z - высота над уровнем земли; Kz - коэффициент турбулентной диффузии; L - характерный масштаб Обухова.

X, м I

Рис. 4.13. Изменение коэффициентов турбулентной диффузии в облаке на разных высотах над уровнем земли.

Kz - коэффициент турбулентной диффузии; X - продольная координата Z -высота над уровнем земли.

На рисунках 4.9-4.13 показаны графики, характеризующие влияние температурной стратификации на параметры турбулентности приземного слоя атмосферы. Как видно из рис. 4.9, температурная устойчивость практически не влияет на формирование профиля скорости в приземной области (до 20 м.). Согласно рис. 4.11, при стремлении характерного масштаба Обухова L к плюс или минус бесконечности, динамическая скорость и* асимптотически приближается к значению, которую принимает при нейтральной устойчивости атмосферы.

Read Save 1*"™*" Stall 17 Sl-fi Ц, DoseJ Картинка |

Г" AutoS ave I? Л»: W Показывать картинку

Данные I ГрФикХ-С Рисинок |

24% об 16.8% об j 7,2? об 2,4% об

2,4 4,8 7,2 9,6 12,0 14,6 17,4 20,4 23.7 27,2 31,0 35,1 39.6 44,4 49,5

55,1

Время 87с

Перемещение jo ^J

Масштаб X j7o~j£j |

Масштаб Y |l0 J

Высота j~

Ув/умен. мае. |1 -^j Координаты/Концентрация

X= 47,8 м

Z- 0,8 и

С 2,429%об

Разрез'

О ПяоскХ-Z

Г Плоек X-Y !

Read

Save

F Stop f JlS"*» | Даргинка | p Показывать картинку

Данные) ГрФикХ-С Рисунок |

I 24% об 116.8% об 7,2% об 2,4% об

На высоте 0,0 м

Перемещение jo jjj

МасштабХ jig

Масштаб Y jp3~jjij

Высота jo jj

Ув/умен. мае. p ^ Координэты/Концентрация

Разрез""".

С ПлоскХ-Z

•" Плоек X-Y

Время 83с

14,6 17.4 20,4 23,7

Рис. 4.14. Интерфейс программы.

Заключение.

1. Разработана методика расчета интенсивности испарения СУГов с поверхности разлива и размеры источника выброса при аварийных ситуациях на магистральных продуктопроводах. Искомыми параметрами для математического описания процесса испарения являются массовая скорость испарения СУГов с единицы поверхности источника и эквивалентный диаметр разлива.

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

3. Выполнены сравнения результатов расчетов по разработанной в диссертации методике и нормативным методикам ОНД-86 и "Токси" с экспериментальными данными, которые показали, что предлагаемая методика лучше согласуется с экспериментом.

4. Результаты научно-исследовательской работы позволили выбрать безопасное месторасположение центрального пункта управления (ЦПУ) вновь строящейся газофракционирующей установки ГФУ-300 управления «Татнефтегазпереработка» ОАО «Татнефть», определить расстояния от трассы продуктопровода от д. Абдрахманово и других производственных объектов с учетом возможности дрейфа облака ШФЛУ в случае аварии на производстве.

Разработанные в диссертации математические модели процесса истечения СУГов из аварийных отверстий и модель распространения паровоздушного облака в атмосфере представляют собой законченную методику моделирования развития аварийной ситуации.

Полный расчет развития аварийной ситуации состоит из двух этапов.

1. Расчет интенсивности выброса СУГов из аварийного отверстия и инте-сивности испарения с поверхности разлива. Процесс выброса СУГов из трубопровода черз свищи происходит вследствие перепада давления в зоне аварийного отверстия. Перепад давления возникает из-за разности давлений внутри и снаружи трубопровода. Вещество, истекающее из аварийного отверстия, условно считается жидким. Одновременно с истечением СУГов из трубопровода происходит адиабатическое расширение части жидкости, попавшей в окружающую среду (в исследовавшемся случае в грунт.). Частично СУГи остаются в жидком состоянии. Оставшаяся неиспаренной жидкость будет иметь температуру кипения СУГов при атмосферном давлении, которая, как правило, существенно ниже температуры окружающего грунта. Вследствие этого между жидкостью и грунтом начинается интенсивный теплообмен с испарением оставшейся жидкости. Параллельно жидкость и пар выталкиваются грунтовыми водами на поверхность земли, образуя лужу неиспаренной жидкости. Здесь возникает дополнительный подвод тепла к жидкости от турбулизовнного атмосферного воздуха. При этом натупает стационарный режим испарения жидкости с поверхности лужи. Искомым параметром для математичаского описания процесса испарения является ма-совая скорость испарения СУГов с единицы поверхности разлива (кг/м2сек.). Для описания всех перечисленных процессов предложены математические модели, реализованные на ЭВМ. Достоверность получаемых по ним результатов подтверждена экспериментальными данными.

2. Нахождение пространственно-временного распределения концентраций. Динамика эволюции паровоздушного облака очень сложна и не может быть с достаточной степенью точности описана простыми алгебраическими уравнениями. Аналитические формулы, предлагаемые в различных методиках ([66, 67]) дают противоречивые результаты. Процесс турбулентного переноса неизотропного неизотермического газа в неограниченном пространстве может быть записан в дифференциальной форме с максимальным учетом всех факторов, влияющих на процесс массопереноса. К ним относятся скорость ветра, атмосферная стратификация, рельеф местности (в настоящей работе рассматривается развитие аварийной ситуации лишь на ровной поверхности) и т.д. Существует множество способов нахождения характеристик турбулентности, представляющие собой хаотические пульсации скорости, температуры, плотности и т.д. В настоящей работе решение уравнений Рейнольдса основано на использовании классического понятия "коэффициента турбулентной диффузии", т. к. имеется множество экспериментальных данных по этим коэффициентам. Численное решение системы уравнений Рейнольдса была осуществлена на ЭВМ в виде программного продукта. Сравнение результатов численных расчетов по выбранной математической модели с экспериментальными данными, взятыми из литературных источников, показало их хорошее согласование. В то же время результаты расчетов, выполненных по нормативным методикам для условий экспериментов, дали заниженные результаты. Проведенные исследования показали, что разработанная в диссертации методика более точно согласуется с экспериментальными данными по сравнению с существующими нормативными методиками. Следовательно, она может быть рекомендована к применению в практических расчетах.

Список литературы диссертационного исследования кандидат технических наук Исламхузин, Дамир Ягфарович, 2003 год

1. Федеральный закон «О промышленной безопасности опасных производственных объектов». - 2-е изд., с изм. — М.: Государственное унитарное предприятие «Научно-технический центр по безопасности в промышленности Госгортехнадзора России», 2003. - 28 с.

2. Порядок разработки декларации безопасности промышленного объекта Российской Федерации. Утвержден приказом МЧС России и Гос-гортехнадзором России 40496 №222/59.

3. Бесчастнов М. В. Промышленные взрывы. Оценка и предупреждение. М.: Химия, 1991. 432 с.

4. Бесчастнов М. В., Соколов В.М, Предупреждение аварий в химических производствах. — М.: Химия, 1979. 392 с.

5. Бесчастнов М. В. Взрывобезопасность и противоаварийная защита химико-технологических процессов. -М.: Химия, 1983.-471 с.

6. Мокроусрв С.Н. Состояние технической безопасности объектов магистральных трубопроводов// Безопасность труда в промышленности. 1998. №9. с. 2-5.

7. Бережковский М.И. Трубопроводный транспорт химических продуктов. Л.: Химия, 1979. 240 с.

8. Кравец В.А. Системный анализ безопасности в нефтяной и газовой промышленности. -М.: Недра, 1984. 117 с.

9. Сафонов B.C., Одишария Г.Э., Швыряев А.А. Теория и практика анализа риска в газовой промышленности. — М.: НУМЦ Минприроды России, 1996, 208 с.

10. Едигаров А.С. Прогноз токсической опасности при аварийных выбросах сероводорода из скважин и промысловых трубопроводов. Сб. на-учн. трудов ВНИИГАЗа "Перспективы развития экологическогострахования в газовой промышленности". М.: ВНИИГАЗ, 1998, с. 49-70.

11. Черкасских В.Г. Проблемы безопасной и безаварийной эксплуатации объектов магистрального, внутри- и межпромыслового транспорта нефти и газа.// Безопасность труда в промышленности. 1997. №5. с. 711.

12. Эльтерман В.М. Охрана воздушной среды на химических и нефтехимических предприятиях. М.: Химия, 1985. 160 с.

13. Бард. В.Д., Кузин А.В. Предупреждение аварий в нефтеперерабатывающей и нефтехимической промышленности. — М.: Химия, 1984. — 315с.

14. Пененко В.В., Алоян А.Я. Модели и методы для задач охраны окружающей среды. Неб.: Наука, 1985,256 с.

15. Легасов В.А., Проблемы безопасного развития техносферы. Коммунист, 1987, №8, с. 92-101.

16. Савельев П.С. Пожары-катастрофы. М.: Стройиздат, 1983. 315 с.

17. Атмосферная турбулентность и моделирование распространения примесей: Пер. с англ. Л.: Гидрометеоиздат, 1985. 351 с.

18. Бетчов Р., Яглом A.M. Замечания о теории подобия для турбулентности в неустойчиво стратифицированной жидкости. Изв. АН СССР. Физ. атмосф. и океана, 1971, т. 7, №12, с. 1270 - 1279.

19. Волковицкий З.И., Иванов В.Н. Диссипация турбулентной энергии в пограничном слое атмосферы. Изв. АН СССР. Сер. Физика, атмосферы и океана, 1970, №5, с. 435-444.

20. Зубковский С.Л., Копров Б.М. О балансе турбулентной энергии в пограничном слое атмосферы. Изв. АН СССР. Сер. Физика, атмосферы и океана, 1970, №10, с. 1307-1310.

21. Kai К. The budget of turbulent energy measured at the ERC 30-m meteorological tower. J. Met. Soc. Japan, 1982, Ser. 11, v 60, No. 5, p. 11171131.

22. Лайхтман Д.Л. Физика пограничного слоя атмосферы. Л.: Гидроме-теоиздат, 1970, 342 с.

23. Кухарец В.П., Цванг Л.Р. О скорости диссипации турбулентной энергии в неустойчиво стратифицированном атмосферном пограничном слое. Изв. АН СССР. Физ. атмосф. и океана, 1977, т. 13, №6, с. 620 -628.

24. Кухарец В.П., Цванг Л.Р. Спектры турбулентных потоков в пограничном слое атмосферы. — Изв. АН СССР. Физ. атмосф. и океана, 1969, т. 5, №11, с. 1132-1142.

25. Кухарец В.П. Спектры пульсаций вертикальной компоненты скорости в пограничном слое атмосферы. Изв. АН СССР. Сер. Физика атмосферы и океана, 1974, №6, с. 613 - 618.

26. Яглом A.M. Закономерности мелкомасштабной турбулентности в атмосфере и океане. Изв. АН СССР. Физ. атмосф. и океана, 1981, т. 17, №12, с. 1235-1257.

27. Обухов A.M. Турбулентность в температурно-неоднородной атмосфере. Труды Института теоретической геофизики АН СССР, 1946, № 1, с. 95-115.

28. Обухов A.M. О структуре температурного поля и поля скоростей в условиях свободной конвекции. Изв. АН СССР. Сер. геофиз., 1960. № 9, с. 1392- 1396.

29. Обухов A.M. Турбулентность и динамика атмосферы. Л.: Гидроме-теоиздат, 1985. 408 с.

30. Роди В. Модели турбулентности окружающей среды. Сб. Методы расчета турбулентных течений. — М.: Мир, 1984.

31. Краснов Н.Ф. Аэродинамика. М. Мир, 1980, Т. 1. 353 е.; Т. 2. 416 с.

32. Монин А.С., Яглом A.M. Статистическая гидромеханика. Ч. 1. М.: Наука, 1965. 640 с.

33. Монин А.С. О температурно неоднородном пограничном слое атмосферы. - Изв. АН СССР. Физ. атмосф. и океана, 1965. т. 1, № 5, с. 490 -500.

34. Монин А.С., Озмидов Р.В. Океанская турбулентность. JL: Гидроме-теоиздат, 1981. 319 с.

35. Монин А.С., Обухов A.M. Основные закономерности турбулентного перемешивания в приземном слое атмосферы. Труды Геофизического института АН СССР, 1954, № 24 (151), с. 163-187.

36. Businger, J.A., J.C. Wyngaard, Y. Izumi and E.F. Brandley, 1971: Fluxpro-file relationship in the atmospheric surface layer. J. Atmos. Sci., 28, 181189.

37. Businger, J.A. and J.W. Deardorff, 1968: On the distinction between "total" heat and eddy heat flux. J. Atmos. Sci., 25, 521-522.

38. Businger, J.A. and S.P.S. Arya, 1974: Height of the Mixed Layer in the Stably Stratified Planetary Boundary Layer. Advances in Geophysics, 18A, Academic Press, N.Y., 73-92.

39. Taylor C.I. Proc. Lond. Math., 1922, Ser. 2, v. 20, p. 3 - 18.

40. Schmidt W. Der Massen austausch in freien Luft. etc. Hamburg, 1925. 20 S.

41. Берлянд M.E. В кн: Современные проблемы климатологии. JL: Гид-рометеоиздат, 1979, с. 280-291.

42. Яглом A.M. О турбулентной диффузии в приземном слое атмосферы. Изв. АН СССР. Физ. атмосф. и океана, 1972. т. 8, № 6, с. 579 - 593.

43. Roberts O.F.T. -Proc. Roy. Soc., 1923, Ser. A, v. 104, p. 640.

44. Сэттон О.Г. Микрометеорология. Пер. с англ. М.: Гидрометеоиздат, 1958. 230 с.

45. Сэттон О.Г. Микрометеорология. Исследование физических процессов в нижних слоях атмосферы. Пер. с англ. Под ред. Д.Л. Лайхтмана. -Л.: Гидрометеоиздат, 1958. 356 с.

46. Sutton O.G. A Theory of Eddy Diffusion in the Atmosphere. Proc. Roy. Soc. (London), Ser. A. 1932. v. 135. p. 143.

47. Frenkiel F.N. Turbulent Diffusion. Adv. Appl. Mech., 1953. v. 90. p. 61.

48. Колмогоров A.H. ДАН СССР, 1941, т. 31, №6, c.538-541.

49. Обухов A.M. Изв. АН СССР. Сер. геогр. и геофиз., 1941, №4-5, с. 512-522.

50. Calder K.L. Quart. J. Mech. Appl. Math., 1943, v.2, p. 153.

51. Лайхтман Д.Л. Физика пограничного слоя атмосферы. Л.: Гидромет-издат, 1970. 180 с.

52. Берлянд М.Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы. Л.: Гидрометеоиздат, 1975. 448 с.

53. Вызова Н.Л. Методическое пособие по расчету рассеяния примесей в пограничном слое атмосферы. М.: Гидрометеоиздат, 1973. 32 с.

54. Вызова Н.Л. Рассеяние примеси в пограничном слое атмосферы. М.: Гидрометеоиздат, 1974,202 с.

55. Вызова Н.Л., Махонько К.П. О взаимодействии аэрозоля с подстилающей поверхностью. Изв. АН СССР. Сер. Физика атмосферы и океана, 1968, т. 4, №9, с. 1000-1003.

56. Вызова Н.Л., Гаргер Е.К., Иванов В.Н. экспериментальные исследования атмосферной диффузии и расчеты загрязнения примеси. Л.: Гидрометеоиздат, 1991. 256 с.

57. Безуглая Э.Ю. Тр. ГГО, вып. 243, 1970. с. 69-79.

58. Скарер Р. Аэродинамика окружающей среды.// Пер. с англ. М.: Мир, 1980. 550 с.

59. Cramer Н.Е. In: Proc. Inst. Nath. Conf. Appl. Meteorol., Amer. Meteorol. Soc., 1957, p. 63-72.

60. Pasquill F. Meteorol. Mag., 1963, v. 90, p. 25-37.

61. Pasquill F. The Estimation of the Dispersion of Windborne Materials. Meteorol. Mag., 1961. v. 90. p. 33-49.

62. Андреев П.И. Рассеяние в атмосферном воздухе газов, выбрасываемых промышленными предприятиями. М.: Госстройиздат, 1952. 88 с.

63. Теверовский Е.Н. В кн.: Новые идеи в области изучения аэрозолей. Изд-во АН СССР, 1949. с. 108-115.

64. Хайкина А.Е. Тр. НИИСТ, 1969. №30, с. 120-125.

65. Recommended Guide for prediction of the dispersion of au borne effluent (sec. ed.). New York, 180 p. The American Society of Mechanical Engineers, 1973, №4.

66. Методика оценки последствий химических аварий (Методика "ТОКСИ". Вторая редакция). М.: НТЦ "Промышленная безопасность", 1999, 83 с.

67. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий (ОНД-86). JL: Гидро-метеоиздат, 1987. 94 с.

68. Атмосферная диффузия. Труды ОРДЕНА ТРУДОВОГО КРАССНОГО ЗНАМЕНИ ГЛАВНОЙ ГЕОФИЗИЧЕСКОЙ ОБСЕРВОТОРИИ им. А.И. Войейкова под ред. д-ра физ-мат. наук М.Е. Берлянда. Л.: Гидро-метеоиздат, 1987. Выпуск 511.

69. Атмосферная диффузия. Труды ОРДЕНА ТРУДОВОГО КРАССНОГО ЗНАМЕНИ ГЛАВНОЙ ГЕОФИЗИЧЕСКОЙ ОБСЕРВОТОРИИ им.

70. А.И. Войейкова под ред. д-ра физ-мат. наук М.Е. Берлянда. Д.: Гидро-метеоиздат, 1985. Выпуск 495.

71. Атмосферная диффузия. Труды ОРДЕНА ТРУДОВОГО КРАССНОГО ЗНАМЕНИ ГЛАВНОЙ ГЕОФИЗИЧЕСКОЙ ОБСЕРВОТОРИИ им. А.И. Войейкова под ред. д-ра физ-мат. наук М.Е. Берлянда. Л.: Гидро-метеоиздат, 1985. Выпуск 467.

72. Маршалл В. Основные опасности химических производств: Пер. сангл.//Под ред. Б.Б. Чайнова, А.Н. Черноплекова. М.: Мир, 1989. 678 с.

73. Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств. ПБ 09-540-03. Утв. постановлением Госгортехнадзора России от 05.05.03г. №29

74. Едигаров А.С. Численный анализ различных моделей турбулентного переноса в задаче диффузии тяжелого газа// ИФЖ, 1991, т. 61, № 3, с. 501-503. Деп. в ВИНИТИ 24.04.91 г., № 1726-В91.

75. Методы расчета турбулентных течений: Пер. с англ.// Под ред. В. Колльмона. М.: Мир, 1984. 464 е., ил.

76. Анучина Н.Н. О методах расчета течений сжимаемой жидкости с большими деформациями. Численные методы механики сплошной среды. Новосибирск: 1970, т. 1, № 4, с.3-84.

77. Фарлоу С. Уравнения с частными производными для научных работников и инженеров. Пер. с англ. Под ред. С.И. Похожаева. М.: Мир. 1985.384 с.

78. Бахвалов Н.С. Численные методы. — М.: Высш. шк., 1973. 631 с.

79. Белоцерковский О.М. Численное исследование современных задач газовой динамики. М.: Высш. шк., 1974. - 397 с.

80. Белоцерковский О.М., Давыдов Ю.М. Метод крупных частиц в газовой динамике. М.: Высш. шк., 1982. - 392 с.

81. Вазов В.Р., Форсайт Дж. Е. Разностные методы для уравнений в частных производных. — М.: Наука, 1963. — 487 с.

82. Годунов С.К. Забродин А.В., Иванов М.Я., Крайко А.Н., Проколов Г.А. Численное решение многомерных задач газовой динамики. М.: Наука, 1976.-400 с.

83. Годунов С.К., Рябенький B.C. Разностные схемы. М.: Наука., 1973. — 400 с.

84. Демидович Б.П., Марон И.А. Основы вычислительной математики. — М.: Высш. шк., 1970. 664 с.

85. Калиткин Н.Н. Численные методы. М.; Высш. шк., 1978. - 512 с.

86. Пасконов В.М., Полежаев В.И., Чудов JI.A. Численное моделирование тепло- массообмена. М.: Наука, 1983. - 288 с.

87. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. М.: Высш. шк., 1970. - 368 с.

88. Самарский А.А. Введение в численные методы. М.: Наука, 1982. -271 с.

89. Самарский А.А. Теория разностных схем. М.: Наука, 1977. - 656 с.

90. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. М.: Наука, 1978, - 592 с.

91. Самарский А.А., Попов Ю.П. Разностные методы решения задач газовой динамики. М.: Наука, 1980. - 352 с.

92. Едигаров А.С. Численный расчет турбулентного течения холодного тяжелого газа в атмосфере// Журнал вычисл. математики и мат. физики. 1991, т.31, № 9, с.1369-1380.

93. Доброчеев О.В. рассеяние тяжелых газов в атмосфере. Физический механизм. Математические модели. М.: РНЦ "Курчатовский институт", 1993, 113 с.

94. Доброчеев О.В., Кулешов А.А., Черноплеков А.Н., Киселев М.Н. Математическая модель рассеяния инертных облаков тяжелых газов при авариях на промышленных предприятиях. М.: Препринт ИАЭ-4789/1,1989,36 с.

95. Shang J.S. An assessment of numerical solutions of the compressible Navier Stokes equations// J. Aircraft. 1985. v. 22. № 5. p. 353-370.

96. Havens J.S., Spicer Т.О., Sheurs P.J. Evaluation of 3-D hydrodynamic computer models for prediction of LNG vapor dispersion in the atmosphere// Internat. Conf. Liquefied Natur. Gas. France. 1985. Session III. № 2. p. 320.

97. Ковеня B.M., Яненко H.H. Метод расщепления в задачах газовой динамики. Новосибирск: Наука, 1981. 304 с.

98. Яненко Н.Н. Численное решение задач механики жидкости. Труды III всесоюз. семинара по моделям механики сплошной среды. Новосибирск: ВЦСО АН СССР, 1976, с. 177-199.

99. Яненко Н.Н. Метод дробных шагов решения многомерных задач математической физики. Новосибирск: Наука, 1967. - 197 с.

100. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1970. 904 с.

101. Абрамович Г.Н. Прикладная газовая динамика. М.: Наука, 1969. 824 с.

102. Spicer Т.О., Havens J.A.// J. hazard, materials. 1987. vol. 16. p. 231-245.

103. Kntha L.H., Phillips O.M., Azard R.S.// J. Flid mech. 1977. vol. 79. P.753. -168.

104. Lofquist K.// Phys. fluids. 1960. vol. 3. № 2, p.158-175.

105. Colebrook C.F., Turbulent flow in pipes with particular reference to the transition region between the smooth and rough pipe laws. J. Institution Civil Engineers 1939.

106. Goldstein S., The similarity theory of turbulence, and flow between planes and trough pipes. Proc. Roy. Soc. A 159,473 (1937).

107. Seinfeld J.H. Air Pollution Physical and Chemical Fundamentals. N.-Y.: McGraw-Hill Co., 1975, 523 p.

108. Seinfeld J.H. Atmospheric Chemistry and Physics of Air Pollution. N.-Y.: J.Wiley, 1986, 738 p.

109. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоатомиздат, 1984. 152 с.

110. Patankar S.V., Liu С.Н., and Sparrow E.M. (1977). Fully Developed Flow and Heat Transfer in Ducts Having Streamwise-Periodic Variations of Cross-Sectional Area, J. Heat Transfer, vol. 99, p. 180.

111. Patankar S.V. and Spalding D.B. (1972a). A Calculation Procedure for the Transient and Steady-State Behavior of Shell-and-Tube heat Exchangers,chap. 7 in Heat Exchangers: Design and Theory Sourcebook, Hemisphere, Washington, D.C.

112. Новиков Ф.А. Дискретная математика для программистов. СПБ.: Питер, 2002. 304 с.

113. Delphi 5. Руководство разработчика, том 1,2. Основные методы и технологии программирования: Пер. с англ. М.: Издательский дом "Вильяме", 2001.

114. Colenbrander G.W., PuttokJ.S. Maplin Sands experiments 1980: interpretation and modeling of liquefied gas spills onto the sea// IUTAM Symp. Atmospheric Dispersion Heavy Gases and Small Particles. Netherlands, 1983. P. 277-295.

115. Chan S.T., Ermak D.L., Morris L.K. //J. Hazardous Materials. 1987. Vol. 16. P. 267-292.1. В.Д. Шашин исемендэге

116. ТАТНЕФТЬ» АЧЫК АКЦИОНЕРЛЫК /&ЭМГЫЯТЕ

117. ТА ТНЕФТЕГАЗПЕРЕРАБОТКА» ИДАРЭСЕ

118. ОТКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО «ТАТНЕФТЬ» им. В.Д. Шашина

119. АБ «Девон-Кредит» г. Альметьевск1. БИК 049202792к/с 301018104000000007921. УТВЕРЖДАЮ»

120. ГТервьш'Заместитель начальника управлениял. /, ч. ^ . ^ - *t / » '

121. Технико-экономический эффект использования результатов работы заключается в обеспечении промышленной безопасности при эксплуатации опасных производственных объектов управления «Татнефтегазпереработка» и в сокращении сроков проектирования.

122. Заместитель главного инженера по промышленной безопасности и охране труда: Начальник технического отдела:1. А.А. Попов1. Г.М. Ахметзянов