CC BY
56
УДК 621.6
Осипова Лилия Эдуардовна
кандидат технических наук, доцент
E-mail: osipovalilij [email protected]
Варсегова Евгения Владиславовна
кандидат физико-математических наук, доцент
E-mail: [email protected]
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зелёная, д. 1
Варсегов Александр Владиславович
инженер
E-mail: alexandervav@mail. ru
ОАО «Газпром трансгаз Казань»
Адрес организации: 420073, Россия, г. Казань, ул. А. Кутуя, д. 41 К расчету распространения поражающих факторов аварий на газопроводах
Аннотация
Постановка задачи. В Российской федерации единая система газоснабжения с диаметром трубопроводов до 2,6 м, которые в течение всего срока эксплуатации подвергаются значительным внутренним и внешним нагрузкам. Поскольку трассы магистральных газопроводов имеют в основном подземную прокладку, проходят по селитебным территориям, пересекают дороги или иные коммуникации, то возможные аварии с выбросом газа с последующим пожаром или взрывом могут привести к катастрофическим последствиям, как для окружающей среды, так и для населения. Расчет интенсивности истечения аварийных выбросов газа и определение зон возможного разлета фрагментов трубы и грунта необходим при расследовании аварии и разработке документов, связанных с безопасностью эксплуатации трубопроводного транспорта, а так же при определении размера компенсационных выплат.
Результаты. Для участка газопровода высокого давления с диаметром d = 1420 мм, рабочим давлением Рнаг = 7,5 МПа и протяженностью L = 120 км рассчитан график суммарной динамики истечения газа (разрыв газопровода - гильотинный, на полное сечение) и определено время полного опорожнения. Определены размеры зон (на поверхности земли) разлета обломков трубы и фрагментов грунта.
Выводы. Значимость полученных результатов для строительной отрасли состоит в том, что они позволяют оценить масштабы возможной аварии уже на стадии проектирования. Найденное пространственное распределение параметров разлета обломков трубы и фрагментов грунта (дальность, высота и время) дает возможность провести оценку последствия ударных воздействий.
Ключевые слова: авария, газопровод, разгерметизация, интенсивность истечения, ударное воздействие.
Введение
Известны работы, в которых рассмотрены методы поддержания эксплуатационной безопасности магистральных газопроводов (МГ) [1-3]. Оптимизация трассировки газовых сетей не только снижает потребление газа, но и повышает безаварийную эксплуатацию газопроводов уже на стадии проектирования [4]. Расчет распространения поражающих факторов аварий на газопроводах позволяет оценить с достаточной долей вероятности возможные зоны поражения и масштабы последствий аварии [5-8].
При оценке последствий возможной аварии в первую очередь рассчитывают интенсивность истечения газа из разрыва трубопровода. Основными параметрами, которые влияют на интенсивность истечения, являются внутренний диаметр газопровода, фактическое доаварийное давление газа в точке разрыва, местоположение точки разрыва и время отсечения аварийных участков. Время отсечения определяется уровнем автоматизации средств технологического контроля и управления. Если операторы с
помощью средств дистанционного контроля перекрывают задвижки газопровода, то время отсечения зависит от алгоритма идентификации факта аварии и реакции диспетчера. Если сработала автоматическая защита газопровода при возникновении аварийной ситуации, то время отсечения минимальное.
Расчет поражающих факторов рекомендует проводить ПАО «Газпром» по методике СТО Газпром 2-2.3-351-2009, где определены методические подходы, термины и определения, а также процедуры и требования к оформлению результатов количественного анализа техногенного риска, обусловленного авариями и инцидентами на опасных производственных объектах транспорта газа и газового конденсата, подземных хранилищ газа. Расчет интенсивности истечения при аварийной разгерметизации газопровода требует строгого определения характеристик транспортируемого газа. В работе [9] приведен расчет параметров природного газа в зависимости от температуры и давления, которые были приняты в дальнейшем при проведении расчета массы аварийного выброса газа.
Постановка задачи
Смоделирован разрыв участка газопровода - «гильотинный», на полное сечение (рис. 1). Разрыв делит газопровод на два аварийных участка Ь1 и Ь2, для каждого участка рассчитывается массовый расход газа 01(() и 02(^.
Рис. 1. Расчетная схема газопровода [9] График суммарной динамики истечения газа
Согласно методическим указаниям СТО Газпром 2-2.3-351-2009 проведен расчет расход газа. Динамика суммарного истечения газа из двух концов трубопровода, представлена на рис. 2. Результаты расчета приведены в табл. 1.
Рис. 2. Динамика суммарного истечения газа из двух концов газопровода (иллюстрация авторов)
Таблица 1
Исходные данные:
Коэффициент гидравлического сопротивления 1 °,°115
Диаметр трубопровода, мм 14°°
Длина трубопровода, км Ь 12°
Протяженность аварийных участков, км Ьи Ь 6°
Температура природного газа в трубопроводе, К —ср 11 281,7
Давление нагнетания, МПа Р 1 наг 7,5
Давление в конце трубопровода, МПа р 1 вс 5,1
Производительность газопровода, м31°9/год е 3°
Плотность природного газа, кг/м3 Рпг °,679
Время закрытия линейных кранов, сек 12°
Показатель адиабаты природного газа к 1,385
Расстояние от компрессорной станции (КС) до середины аварийных участков, км Х] — Х2 1°
Результаты расчета: первый второй участок участок
Давление в момент аварии в точке разрыва газопровода, МПа Р° 6,41
Масса газа, которая нагнетается в аварийный участок газопровода до момента отсечения, кг мкс 828°°,°
Масса газа, находящаяся в аварийном участке газопровода до аварии, кг мг 4186732,7 4°21123,7
Масса газа, истекающего в адиабатическом режиме, кг мн 341744,3 328227,4
Начальный критический массовый расход газа, кг/с о° 17774,1 17774,1
Параметр пь °,23 °,22
Постоянная времени, с £ь 2°38,7 2°38,7
Скорость звука в газе до разрыва, м/с а° 445,7 445,7
Масса газа, выброшенного из аварийной секции на первом этапе истечения, кг м1 448°55,2 43998°,4
Расход газа на момент закрытия линейного крана, кг/с О(Гл) 2821,1 2697,6
Постоянная времени, с £х 138,7 138,7
Расход газа после перекрытия линейного участка (" = 6° сек), кг/с ОО) 183°,3 175°,2
Масса газа, выброшенного из аварийной секции на втором этапе истечения, кг м2 39125°,9 374129,1
Полная масса газа, кг М 8393°6,6 8141°9,5
Расчет разлета фрагментов грунта и осколков трубы
При разрыве подземного газопровода количество образующихся фрагментов грунта и осколков трубы, их форма и направление полета являются величинами случайными. Принято считать, что направления движения осколков равновероятны, рассеивание осколков по площади зависит только от расстояния от места разрыва газопровода, а зона возможного поражения на поверхности земли имеет форму круга. Согласно методическим указаниям СТО Газпром 2-2.3-351-2009 проведен расчет разлета фрагментов грунта и осколков трубы.
Предполагается, что часть потенциальной энергии сжатого газа расходуется на фрагментацию оболочки трубы и слоя грунта, и приобретение ими начальной скорости. В дальнейшем рассчитывается динамика движения этих фрагментов в воздухе, и определяются размеры зон (на поверхности земли) их разлета.
Закон сохранения энергии на единицу длины траншеи имеет вид:
Мгр7'=л-А -мт, (1)
где Мг - масса сжатого газа для единичной длины газопровода, кг/м; Аг - работа расширения единицы массы газа в предположении адиабатического процесса, Дж/кг; Мгрунта - масса метаемого грунта, кг/м; у°° - скорость метаемой массы грунта, м/с.
От общей энергии взрыва конденсированного взрывчатого вещества на образование воронки расходуется значительная часть п энергии взрыва: ~ 20 % - для плотных суглинков и глин, до 30^35 % - для средних грунтов. В расчетах следует принимать значение (консервативная оценка) п = 0,35.
Уравнение для Мгрунта имеет вид:
^^ грунта Р гр
(H + DH )2tg ф-pR2-(DH - R0) tg jj, (2)
DH (3)
sin ф
Mr =Рг,тр -P-Ro2' (4)
P
Рг, тр = тпг R J" Z ' (5)
Rr "Ал- Zc
где ргр - плотность грунта, кг/м3; рг, ^ - плотность природного газа в трубопроводе, кг/м3.
Для адиабатического расширения значение удельной потенциальной энергии сжатого газа определяют по формуле:
1-AP"4 ' (6)
Рг,тр (к -1) L J
AP = РТр - Pc. (7)
Из уравнения (1) определяется скорость выбрасываемых фрагментов оболочки трубы, которая по принятому приближению совпадает со скоростью метаемой массой грунта v0.
После выброса из траншеи фрагменты трубы и грунта продолжают свое движение в атмосфере. Уравнения движения фрагментов грунта и оболочки трубы в воздухе задаются системой уравнений:
4 -■ DP
dy
— - v, dt
dv 1 ,2 Г2 2 mo ~T = ~ 2 VAVv + w -v - mog'
2 dx ^
— - w, dt
dw 1 ,2 г mo^7 1ф cxPWv
"о 2 »ф-^^+ ^ ,
где у, V - вертикальные значения расстояния и скорости; х, я - горизонтальные значения расстояния и скорости; /ф, 1х, 1у, /2 - эффективный и характерные размеры фрагмента по трем ортогональным осям; сх - коэффициент сопротивления фрагмента в воздухе (для плохо обтекаемого тела сх = 0,2); рв - плотность воздуха:
/ф =7 ¡2 + ¡2+¡2, (9)
¡х = 5Тр, (10)
¡у = я-Д), (11)
/2 = ^, (12)
Поск
где 5тр - толщина стенки, м; поск - общее количество одинаковых осколков,
образующихся при разрушении МГ (рекомендуемое значение поск = 3^5). Масса фрагмента трубы эффективного размера /ф:
т0 = Ртр1ф , (13)
где ртр - плотность оболочки трубы, кг/м .
Уравнения для скорости системы (8) примут следующий вид:
ЖV сх Рв нг 2 —=-22+ •V - я,
ж 2/ф Рф
дм> с рв г~2 2 — =--ху—VV + м> •я,
ж 2/ф Рф
где Рф - плотность фрагмента трубы или грунта, кг/м3.
(14)
Результаты расчета представлены в табл. 2.
Таблица 2
Исходные данные:
Давление природного газа в газопроводе, МПа р 1 тр 7,5
Атмосферное давление, кПа рс 1°1,3
Диаметр газопровода, мм о° 14°°
Толщина стенки газопровода, мм 5тр 21,5
Плотность природного газа, кг/м3 Рпг °,679
Плотность грунта, кг/м3 Ргрунта 2°°°
Плотность материала оболочки трубы, кг/м3 Ртр 78°°
Показатель адиабаты природного газа к 1,385
Длина разрыва, м Ьразр 1°
Число осколков поск
Расстояние от оси газопровода до поверхности грунта, м и 1,°
Результаты расчета:
и/я 2,43
Плотность газа в газопроводе, кг/м3 рг, тр 51,93
Разность давлений, Па АР 73987°°
Работа расширения, Дж/кг Аг 462528,75
Ф 3° 45 6°
АН 1,4° °,99 °,81
Масса грунта, кг ■^^грунта 97°5 11226 13775
Начальная скорость метания грунта (фрагмента трубы), м/с 58,94 48,°2 37,23
Уравнения для скорости движения фрагментов грунта и оболочки трубы в воздухе задаются линейной системой дифференциальных уравнений (14), которая была решена численным методом Рунге-Кутта четвертого порядка.
На рис. 3-4 представлены найденные графические зависимости вертикальных и горизонтальных составляющих скорости от времени полета фрагмента трубы и грунта при угле раскрытия ф — 3°о.
Зная зависимости составляющих скоростей, численным интегрированием находятся зависимости от времени дальности полета и высоты подъема вылетевшего фрагмента трубы (грунта):
х^) = | w(t )Л
с t
у а) = | )А .
(15)
Графические зависимости дальности полета, высоты подъема и времени полета фрагмента трубы и грунта при этом же угле раскрытия ф — 3°о приведены на рис. 5-6.
Рис. 3. Зависимость от времени вертикальных Рис. 4. Зависимость от времени горизонтальных составляющих скорости (иллюстрация авторов) составляющих скорости (иллюстрация авторов)
Рис. 5. Зависимость дальности разлета грунта (трубы) от времени (иллюстрация авторов)
Рис. 6. Зависимость высоты разлета грунта (трубы) от времени (иллюстрация авторов)
Анализируя рис. 3-6 можно сделать следующие выводы:
- по осям х и у скорость имеет линейный характер;
- при начальной скорости метания грунта (фрагмента трубы) у0 =58,94 м/с примерно за 12 сек. фрагменты грунта (трубы) достигнут поверхности земли;
- за счет высокого давления транспортируемого природного газа и развиваемой скорости, после разрыва, фрагменты грунта (трубы) могут вылететь в атмосферу на достаточно большие расстояния, так по оси х на 700 м и по оси у до 180 м.
Найденные значения позволяют получить пространственное распределение параметров разлета обломков и грунта (дальность, высота и время), что соответственно дает возможность провести оценку последствия ударных воздействий (рис. 7).
—---Раэ.чег
Рис. 7. Траектории разлета обломков и грунта при аварии (иллюстрация авторов)
Таблица 3
Pг 7,5 МПа
нт 2,43
ф 30° 45° 60°
т Фрагменты грунта
1 г 171 147 117
1 кг 507 370 242
1 т 679 457 278
Фрагменты оболочки
1 г 292 237 174
1 кг 607 423 265
1 т 696 465 281
В табл. 3 приведены расчеты дальности разлета осколков для масс фрагментов от 1 грамма до 1000 кг в предположении, что фрагменты имеют компактную форму (линейные размеры по трем ортогональным осям фрагментов равны).
Заключение
Для принятия обоснованных инженерных и управленческих решений по безопасной эксплуатации газопроводов необходимо знать характеристики транспортируемого природного газа, разбираться в природе поражающих факторов возможной аварии, а инженерные расчеты позволяют проводить экспресс оценки возможной аварийной ситуации.
Список библиографических ссылок
1. Халлыев Н. Х., Усманов Р. Р., Чучкалов М. В., Велиюлин И. И. Методология поддержания, продления и обновления линейной части магистральных газопроводов - основа надежной и безопасной эксплуатации газотранспортной системы // Газовая промышленность. 2017. № 9 (757). С. 34-40.
2. Ямаева Э. Г., Фомина Е. Е. Разработка балльной оценки факторов влияния на безопасную эксплуатацию объектов газораспределения на этапе проектирования // Безопасность жизнедеятельности. 2016. № 1 (181). С. 18-23.
3. Шебеко Ю. Н., Гордиенко Д. М., Пономаренко А. А. и др. Оценка пожарного риска линейной части магистральных трубопроводов // Пожарная безопасность. 2010. № 4. С. 47-58.
4. Барышева О. Б., Садыков Р. А., Батюшков Н. Ю. Оптимизация трассировки газораспределительной сети низкого давления // Современное строительство и архитектура. 2016. № 3 (03). С. 25-28.
5. Рыбаков А. В., Замана К. Ю. Модель оценки зоны возможного поражения от разлета осколков фрагментов оборудования при возникновении чрезвычайных ситуаций техногенного характера на пожаровзрывоопасных объектах // Научные и образовательные проблемы гражданской защиты. 2015. № 4. С. 56-64.
6. Metropolo P., Brown A. Natural gas pipeline accident consequence analysis // Proc. Safety Prog. 2004. 23: doi: 10.1002 / prs.10054. P. 307-310.
7. Сафонов В. С., Ковалев С. А., Овчаров С. В., Мельников А. В. Оценка показателей безопасности и риска перспективных газопроводов высокого давления // Вестник РУДН. Экология и безопасность жизнедеятельности. 2007. № 1. С. 54-62.
8. Rothwell B., Stephens M., Risk analysis of sweet natural gas pipelines: Benchmarking simple consequence models. In 2006 International Pipeline Conference. American Society of Mechanical Engineers, 2006. P. 911-919.
9. Быков А. И. Определение средней температуры газа на аварийном участке магистрального газопровода // Пожарная безопасность. 2015. Т. 24, № 6. С. 43-50.
Osipova Liliya Eduardovna
candidate of technical sciences, associate professor E-mail: osipovalilij [email protected] Varsegova Evgeniya Vladislavovna
candidate of physical-mathematical sciences, associate professor E-mail: [email protected]
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1 Varsegov Alexander Vladislavovich
engineer
E-mail: alexandervav@mail. ru OJSC «Gazprom transgaz Kazan»
The organization address: 420073, Russia, Kazan, A. Kutuya st., 41
To the calculation of the distribution of damaging factors of accidents on gas pipelines
Abstract
Problem statement. In the Russian Federation, a unified gas supply system with a pipe diameter of up to 2,6 m, which is subjected to significant internal and external loads during the whole period of operation. Since the main gas pipeline routes are mainly underground, they pass through residential areas, cross roads or other communications, possible accidents with gas emissions followed by fire or explosion can lead to disastrous consequences both for the environment and for the population. The calculation of the intensity of the outflow of emergency gas emissions and the determination of zones of possible dispersion of pipe and soil fragments is necessary in the investigation of the accident and the development of documents related to the safety of operation of pipeline transport, as well as in determining the amount of compensation payments.
Results. For the section of the high pressure gas pipeline with a diameter of d = 1420 mm, a working pressure PHar = 7,5 MPa and a length of L = 120 km, a graph of the total gas outflow dynamics was calculated (the gas pipeline was guillotine, for a full cross section) and the time for complete emptying was determined. The sizes of zones (on the surface of the earth) of scattering of pipe fragments and soil fragments have been determined.
Conclusions. The significance of the results obtained for the construction industry is that they allow us to estimate the scale of a possible accident already at the design stage. The found spatial distribution of the scattering parameters of pipe fragments and soil fragments (range, height and time) makes it possible to assess the impact of shock effects.
Keywords: accident, gas pipeline, depressurization, flow rate, shock effect.
References
1. Khallyev N. Kh., Usmanov R. R., Chuchkalov M. V., Veliyulin I. I. Methodology for the maintenance, renewal and renovation of the linear part of main gas pipelines is the foundation of reliable and safe operation of the gas transport system // Gazovaya promishlennost. 2017. № 9 (757). P. 34-40.
2. Yamaeva E. G., Fomina E. E. Development of scoring factors influencing the safe operation of the gas distribution in the design phase // Bezopasnost zhiznedeyatelnosti. 2016. № 1 (181). P. 18-23.
3. Shebeko Yu. N., Gordienko D. M., Ponomarenko A. A. Fire risk assessment of the linear part of trunk pipeline // Pozharnaya bezopasnost. 2010. № 4. P. 47-58.
4. Barysheva O. B., Sadykov R. A., Batyushkov N. Yu. Optimization of trace in gas distribution network of flow pressure // Sovremennoe stroitelstvo i arkhitektura. 2016. № 3 (03). P. 25-28.
5. Rybakov A. V., Zamana K. Yu. The model for assessing the area of possible damage from the fragmentation of fragments of equipment fragments in the event of man-made emergencies at fire and explosion hazard objects // Nauchnye i obrazovatelnye problem grazhdanskoi zashchity. 2015. № 4. P. 56-64.
6. Metropolo P., Brown A. Natural gas pipeline accident consequence analysis // Proc. Safety Prog. 2004. 23: doi: 10.1002 / prs.10054. P. 307-310.
7. Safonov V. S., Kovalev S. A., Ovcharov S. V., Melnikov A. V. Assessment of safety performance and risk of promising high-pressure gas pipelines // Vestnik RUDN. Ekologiya i bezopastnost zhiznedeyatelnosti. 2007. № 1. P. 54-62.
8. Rothwell B., Stephens M. Risk analysis of sweet natural gas pipelines: Benchmarking simple consequence models. In 2006 International Pipeline Conference. American Society of Mechanical Engineers, 2006. P. 911-919.
9. Bykov A. I. Determination of the average gas temperature in the emergency section of the main pipeline // Pozharnaya bezopasnost. 2015. Vol. 24. № 6. P. 43-50.
