Приказ Ростехнадзора от 10.01.2023 № 4

«Об утверждении Руководства по безопасности „Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи“»

ФЕДЕРАЛЬНАЯ СЛУЖБА ПО ЭКОЛОГИЧЕСКОМУ, ТЕХНОЛОГИЧЕСКОМУ И АТОМНОМУ НАДЗОРУ

ПРИКАЗ

от 10 января 2023 г. № 4

ОБ УТВЕРЖДЕНИИ РУКОВОДСТВА ПО БЕЗОПАСНОСТИ «МЕТОДИКА АНАЛИЗА РИСКА АВАРИЙ НА ОПАСНЫХ ПРОИЗВОДСТВЕННЫХ ОБЪЕКТАХ НЕФТЕГАЗОДОБЫЧИ»

В соответствии с пунктом 5 статьи 3 Федерального закона от 21 июля 1997 г. № 116-ФЗ «О промышленной безопасности опасных производственных объектов», пунктом 1 Положения о Федеральной службе по экологическому, технологическому и атомному надзору, утвержденного постановлением Правительства Российской Федерации от 30 июля 2004 г. № 401, приказываю:

1. Утвердить прилагаемое Руководство по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи».

2. Признать утратившим силу приказ Федеральной службы по экологическому, технологическому и атомному надзору от 17 августа 2015 г. № 317 «Об утверждении Руководства по безопасности „Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи“.»

Руководитель А. В. Трембицкий

Утверждено приказом Федеральной службы по экологическому, технологическому и атомному надзору от 10 января 2023 г. № 4

РУКОВОДСТВО ПО БЕЗОПАСНОСТИ «МЕТОДИКА АНАЛИЗА РИСКА АВАРИЙ НА ОПАСНЫХ ПРОИЗВОДСТВЕННЫХ ОБЪЕКТАХ НЕФТЕГАЗОДОБЫЧИ»

I. ОБЩИЕ ПОЛОЖЕНИЯ

1. Руководство по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи» (далее — Руководство) разработано в целях содействия соблюдению требований федеральных норм и правил в области промышленной безопасности «Правила безопасности в нефтяной и газовой промышленности», утвержденных приказом Ростехнадзора от 15 декабря 2020 г. № 534 (зарегистрирован Минюстом России 29 декабря 2020 г., регистрационный № 61888).

2. Руководство содержит методические рекомендации по проведению анализа риска аварий на опасных производственных объектах нефтегазодобычи (далее — ОПО НГД) и не является нормативным правовым актом.

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

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

5. Руководство не устанавливает необходимость, периодичность проведения анализа риска аварий на ОПО НГД.

II. РЕКОМЕНДАЦИИ ПО АНАЛИЗУ РИСКА АВАРИЙ НА ОПАСНЫХ ПРОИЗВОДСТВЕННЫХ ОБЪЕКТАХ НЕФТЕГАЗОДОБЫЧИ

6. Общая процедура анализа риска аварий на ОПО НГД включает:

планирование и организацию работ на ОПО НГД;

идентификацию опасностей;

оценку риска аварий на ОПО НГД;

определение степени опасности ОПО НГД и (или) их участков;

разработку рекомендаций по уменьшению рисков аварий на ОПО НГД.

7. Выбор методов и оформление результатов оценки риска аварий на ОПО НГД проводится в соответствии с Руководством по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах», утвержденным приказом Ростехнадзора от 3 ноября 2022 г. № 387, с учетом настоящего Руководства.

8. Расчет пожарного риска на объектах ОПО НГД и сравнение его с законодательно установленными допустимыми значениями пожарного риска осуществляются в соответствии с Федеральным законом от 22 июля 2008 г. № 123-ФЗ «Технический регламент о требованиях пожарной безопасности».

9. Исходные данные, сделанные допущения и предположения, результаты оценки риска аварий на ОПО НГД обосновываются и документально фиксируются в объеме, достаточном для того, чтобы выполненные расчеты и выводы могли быть повторены и проверены в ходе независимой оценки риска аварий (например, страховщиком в рамках выполнения Федерального закона от 27 июля 2010 г. № 225-ФЗ «Об обязательном страховании гражданской ответственности владельца опасного объекта за причинение вреда в результате аварии на опасном объекте») или экспертизы промышленной безопасности.

10. Рекомендации к оформлению результатов оценки риска аварий на ОПО НГД приведены в Руководстве по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах», утвержденном приказом Ростехнадзора от 3 ноября 2022 г. № 387.

11. Исходными данными для количественной оценки риска аварий на ОПО НГД являются:

сведения о природно-климатических условиях района и площадки строительства (трассы);

данные о технологии, оборудовании и технических решениях по обеспечению безопасности ОПО НГД;

сведения о работниках ОПО НГД, населении и других объектах на прилегающей территории;

графические материалы (ситуационный и генеральный планы, иные графические материалы).

В случае недостаточности исходных данных для качественной и количественной оценки риска аварий на ОПО НГД рекомендуется обосновывать необходимые допущения и учитывать их при анализе неопределенностей результатов оценки риска аварий на ОПО НГД.

12. Сведения о природно-климатических условиях района и площадки строительства (трассы) включают в себя:

характеристику рельефа местности (например, ровный, холмистый, пересеченный, горный) с упоминанием наличия балок, оврагов, возвышенностей, естественных и искусственных подземных полостей, горных выработок;

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

сведения о сейсмичности площадки строительства (трассы) и наличии других опасных природных воздействий, отнесенных согласно своду правил СП 115.13330.2016 «Геофизика опасных природных воздействий. Актуализированная редакция СНиП 22-01-95», утвержденному и введенному в действие приказом Минстроя России от 16 декабря 2016 г. № 956/пр, к категориям «весьма опасные» и «чрезвычайно опасные (катастрофические)»;

данные о находящихся вблизи объекта (по трассе) водных объектах.

13. Данные о технологии, оборудовании и технических решениях по обеспечению безопасности ОПО НГД включают в себя:

описание технологических процессов и технологические схемы;

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

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

показатели пожаровзрывоопасности и токсичной опасности веществ и материалов, используемых в технологическом процессе;

описание объемно-планировочных и конструктивных решений зданий и сооружений;

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

14. Сведения о работниках ОПО НГД, населении и других объектах на прилегающей территории включают в себя:

данные о режиме работы опасного производственного объекта, количестве смен, списочной численности работников эксплуатирующей организации, численности наибольшей смены, распределении людей по территории объекта;

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

перечень населенных пунктов и мест массового скопления людей, находящихся в зоне действия поражающих факторов при максимальной аварии на ОПО НГД (с указанием характера застройки, максимальной и средней численности людей с учетом их суточной миграции);

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

15. Графические материалы включают в себя:

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

генеральный план площадочного ОПО НГД с экспликацией всех зданий и сооружений;

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

16. Количественная оценка риска аварий на ОПО НГД включает в себя следующие процедуры:

идентификация опасностей, которые могут привести к авариям, и определение вероятностей (частот) возникновения аварий (с использованием метода анализа «деревьев отказов»);

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

оценка количества опасных веществ, участвующих в аварии и создании поражающих факторов по каждому сценарию;

расчет зон действия поражающих факторов по каждому сценарию;

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

расчет и представление показателей риска аварий на ОПО НГД;

определение при необходимости степени опасности аварий на ОПО НГД и (или) наиболее опасных составных элементов ОПО НГД.

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

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

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

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

Рекомендуемые частоты разгерметизации емкостного оборудования и трубопроводов приведены в приложении № 4 к Руководству по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах», утвержденному приказом Ростехнадзора от 3 ноября 2022 г. № 387. Рекомендованные частоты разгерметизации и фонтанирования скважин на разных этапах их производственного цикла приведены в приложении № 5 к Руководству.

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

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

В сценариях развития возможных аварий учитываются следующие явления, связанные с условиями обращения и выброса опасных веществ на ОПО НГД:

пролив (разлив) жидкости;

испарение пролива (разлива) жидкости;

пожар пролива (разлива), в том числе горение в резервуаре;

разрыв сосудов (аппаратов), трубопроводов с образованием воздушной волны сжатия (ударной волны) и осколков;

дрейф облака топливно-воздушных смесей (далее — ТВС);

распространение токсического облака;

пожар-вспышка;

вскипание и выброс горящей жидкости из резервуара при пожаре;

истечение газа при фонтанировании скважин;

струйное горение (факел);

взрыв облака ТВС (дефлаграция или детонация) как в замкнутом помещении, так и на открытом или частично загроможденном пространстве;

взрыв расширяющихся паров вскипающей жидкости (BLEVE);

диффузионное горение облака ТВС (огненный шар);

пожар в помещении (в том числе с учетом действий по тушению).

Примеры типовых сценариев развития возможных аварий приведены в приложении № 1 к Руководству.

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

20. Массу аварийного выброса опасных веществ рекомендуется определять как массу вещества в аппарате (трубопроводе) с учетом перетоков от соседних аппаратов (участков) в течение времени обнаружения выброса и перекрытия запорной арматуры (задвижек) с учетом массы стока вещества из отсеченного блока (трубопровода).

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

22. Рекомендуемый порядок расчета истечения опасных веществ из технологических трубопроводов на площадочных объектах приведен в приложении № 4 к Руководству по безопасности «Методика оценки риска аварий на опасных производственных объектах нефтегазоперерабатывающей, нефте- и газохимической промышленности», утвержденному приказом Ростехнадзора от 28 ноября 2022 г. № 414.

23. Количество веществ, участвующих в аварии на промысловых трубопроводах ОПО НГД, определяется с учетом данных, указанных в приложении № 2 к Руководству.

24. При расчете истечения газа при фонтанировании скважин рекомендуется использовать зависимости, указанные в приложении № 3 к Руководству.

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

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

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

28. Для сценариев с пожаром пролива (разлива) в случае примерно равных площадей пролива (разлива) форму пламени при горении рекомендуется аппроксимировать наклонным цилиндром с радиусом, равным эффективному радиусу пролива (разлива). Для этого цилиндра определяются параметры теплового излучения в соответствии с пунктом 23 приложения № 3 к Методике определения расчетных величин пожарного риска на производственных объектах, утвержденной ведомственным приказом МЧС России (далее — Методика определения расчетных величин пожарного риска на производственных объектах).

29. Для расчета сценариев с образованием огненного шара рекомендуется использовать пункт 24 приложения № 3 к Методике определения расчетных величин пожарного риска на производственных объектах.

30. Для расчета концентрационных полей при рассеивании и дрейфе облака рекомендуется использовать Руководство по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденное приказом Ростехнадзора от 2 ноября 2022 г. № 385.

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

32. Массу во взрывоопасных пределах, способную участвовать во взрыве, рекомендуется определять согласно приложению № 3 к федеральным нормам и правилам в области промышленной безопасности «Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств», утвержденным приказом Ростехнадзора от 15 декабря 2020 г. № 533 (зарегистрирован Минюстом России 25 декабря 2020 г., регистрационный № 61808) (далее — федеральные нормы и правила в области промышленной безопасности «Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств»).

33. При отсутствии сведений о распределении источников воспламенения и вероятности зажигания облака ТВС расчет зон поражения при взрыве облаков ТВС рекомендуется выполнять из условия воспламенения облака ТВС в момент времени, когда облако ТВС достигает наибольшей массы, способной к воспламенению.

34. Для расчета параметров волн давления (давление P и импульс I), образующихся при сгорании или взрыве облаков ТВС, и зон поражения рекомендуется использовать Руководство по безопасности «Методика оценки последствий аварийных взрывов топливно-воздушных смесей», утвержденное приказом Ростехнадзора от 28 ноября 2022 г. № 412.

35. Последствия сценария со струйным горением и расчета размеров зон поражения термическим излучением рекомендуется определять в соответствии с Методикой определения расчетных величин пожарного риска на производственных объектах.

36. Для расчета размеров зон поражения ударными волнами и расчета вероятности гибели людей, находящихся в зданиях, при взрыве рекомендуется использовать пробит-функцию в соответствии с подпунктами 2.2 и 2.3 приложения № 3 к федеральным нормам и правилам в области промышленной безопасности «Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств».

37. Для оценки гибели людей при пожарах в здании с учетом их эвакуации рекомендуется использовать формулы в соответствии с приложением № 4 к Методике определения расчетных величин пожарного риска на производственных объектах.

38. Для расчета вероятности гибели людей от поражения токсичными опасными веществами рекомендуется применять формулы согласно Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному приказом Ростехнадзора от 2 ноября 2022 г. № 385.

39. Для расчета вероятности гибели людей от поражения токсичными продуктами горения в помещениях рекомендуется применять формулы согласно приложению № 4 к Методике определения расчетных величин пожарного риска на производственных объектах.

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

41. При оценке зоны разлета осколков оборудования под давлением рекомендуется руководствоваться положениями Руководства по безопасности «Методика оценки риска аварий на опасных производственных объектах магистрального трубопроводного транспорта газа», утвержденного приказом Ростехнадзора от 22 декабря 2022 г. № 454 (далее — Руководство по безопасности «Методика оценки риска аварий на опасных производственных объектах магистрального трубопроводного транспорта газа».

42. Для проведения оценки возможного числа пострадавших и материального ущерба от аварий на ОПО НГД используются критерии и законы поражения людей и разрушения (повреждения) оборудования, зданий и сооружений, приведенные в приложении № 3 к федеральным нормам и правилам в области промышленной безопасности «Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств», пункте 40 Руководства по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденного приказом Ростехнадзора от 2 ноября 2022 г. № 385, в Руководстве по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах», утвержденном приказом Ростехнадзора от 3 ноября 2022 г. № 387.

43. Оценка возможного числа пострадавших от аварий на ОПО НГД проводится на основе:

полученных зон действия поражающих факторов;

критериев и законов поражения людей на открытой местности и в зданиях;

распределения и времени нахождения людей в зоне действия поражающих факторов.

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

44. Оценка материального ущерба от аварий на ОПО НГД производится на основе расчета вреда (экономического ущерба и вреда окружающей среде) от аварии, а также справки о размере причиненного вреда и оценке экономического ущерба и вреда окружающей среде от аварии.

45. Для определения основных показателей риска аварий на ОПО НГД, приведенных в приложении № 4 к Руководству, рекомендуется использовать соотношения, приведенные в Руководстве по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах», утвержденном приказом Ростехнадзора от 3 ноября 2022 г. № 387.

46. Качественная оценка риска аварий на ОПО НГД включает следующие процедуры:

оценка качественных признаков возможности возникновения инцидентов и (или) аварий;

оценка качественных признаков тяжести возможных последствий анализируемых инцидентов и (или) аварий;

определение качественных характеристик степени опасности аварий на ОПО НГД и (или) составных элементах ОПО НГД.

В процедуре качественной оценки риска аварий на ОПО НГД используются экспертные специальные знания в области аварийности и травматизма на ОПО НГД.

Приложение № 1 к Руководству по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи», утвержденному приказом Федеральной службы по экологическому, технологическому и атомному надзору от 10 января 2023 г. № 4

ТИПОВЫЕ СЦЕНАРИИ РАЗВИТИЯ ВОЗМОЖНЫХ АВАРИЙ

Примеры «деревьев событий» для сценариев A₁ - A₁₂, связанных с выбросом нефти и нефтепродуктов на площадочных объектах, приведены на рисунках 1 - 12 приложения № 1 к Руководству.

На рисунке 1 приложения № 1 к Руководству принимают следующие условные вероятности событий:

резервуар теряет целостность после появления разрушения (a) — условная вероятность определяется согласно данным Руководства по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах»;

разрушение соседних (находящихся в одном обваловании) резервуаров и дополнительный выброс нефти, нефтепродуктов (b):

для длительных выбросов — 0;

для залповых выбросов — определяется на основе расчета зон действия опасных факторов;

возможность пролива за пределы обвалования определяется по результатам моделирования распространения жидкости по поверхности, и условная вероятность (c) его составляет:

при длительном выбросе — 1, если приподнятая струя, образующаяся при истечении из резервуара, выпадает за пределы обвалования, в противном случае — 0;

при залповом выбросе — 0, если конструкция обвалования вмещает всю выброшенную нефть, нефтепродукты, исключает перехлест нефти, нефтепродукта через обвалование и его разрушение/размыв, в противном случае — 1;

мгновенное воспламенение и образование горящих проливов (d) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы;

образование дрейфующего облака ТВС (e) для взрывопожароопасных жидкостей с давлением насыщенных паров в условиях содержания в емкости менее 10 кПа — 0, в остальных случаях — 1;

образование капельной взвеси нефти, нефтепродукта в атмосфере (f): для бензинов и керосинов при высоте выброса более 5 м — 1; в остальных случаях — 0;

отсроченное воспламенение (g) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы.

Рис. 1. «Дерево событий» разрушения/перелива наземного резервуара (сценарий A1)

Рис. 2. «Дерево событий» при выходе газовой фазы с наземного резервуара (сценарий А2)

На рисунке 2 приложения № 1 к Руководству приняты следующие условные вероятности событий:

воспламенение шлейфа паров нефти, нефтепродукта (a) — 0,2;

прекращение горения (b) при наличии на дыхательной арматуре исправного огнепреградителя — 1, на зазоре — 0,75, на люке — 0,2;

зажигание нефти, нефтепродукта в резервуаре (c) — 0,2, в зависимости от надежности огнепреградителей или с учетом способности потушить пожар на зазоре;

при переходе горения на резервуар в резервуаре происходит взрыв (d) — 0,2 для резервуаров со стационарной крышей, для резервуаров с плавающей крышей — 0;

взрыв вызывает разрушение резервуара (в том числе обрушение крыши с переливом горящего продукта) (e) — 0,5;

при проведении пенной атаки произошел перелив нефти, нефтепродукта — 0,2 (без пенной атаки не задается);

выброс горящей нефти, нефтепродукта при вскипании воды в резервуаре (g) — в зависимости от обстоятельств;

потеря устойчивости резервуара при пожаре в нем (h) — в зависимости от обстоятельств.

На рисунке 3 приложения № 1 к Руководству приняты следующие условные вероятности событий:

при взрыве внутри резервуара образуются разлетающиеся элементы резервуара (a) — 0,02;

условная вероятность выхода нефти из резервуара (b) определяется согласно данным Руководства по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах» из отношения частоты полного разрушения к частоте разгерметизации;

зажигание нефти, нефтепродукта в резервуаре при отсутствии выброса из резервуара (c) — 0,2.

Рис. 3. «Дерево событий» при взрыве внутри наземного резервуара (сценарий A3)

Рис. 4. «Дерево событий» разрушения/перелива подземного резервуара (сценарий A4)

На рисунке 4 приложения № 1 к Руководству приняты следующие условные вероятности событий:

мгновенное воспламенение и образование горящих проливов (d) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы;

образование дрейфующего облака ТВС (e) для взрывопожароопасных жидкостей с давлением насыщенных паров в условиях содержания в емкости менее 10 кПа — 0, в остальных случаях — 1;

отсроченное воспламенение (g) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы.

Рис. 5. «Дерево событий» при выходе газовой фазы из подземного резервуара (типа ЖБР) (сценарий A₅)

На рисунке 5 приложения № 1 к Руководству приняты следующие условные вероятности событий:

воспламенение шлейфа паров нефти, нефтепродукта (a) — 0,2;

прекращение горения (b) — 0,75;

зажигание нефти, нефтепродукта в резервуаре (c) — 0,2 или в зависимости от надежности огнепреградителей;

при переходе горения на резервуар в резервуаре происходит взрыв (d) — 0,2;

взрыв вызывает разрушение резервуара (разрушение крыши) (e) — 0,5;

при проведении пенной атаки произошел перелив нефти, нефтепродукта — 0,2 (в случае отсутствия данного варианта — пенной атаки — не задается);

выброс горящей нефти, нефтепродукта при вскипании воды в резервуаре (g) — в зависимости от обстоятельств.

Рис. 6. «Дерево событий» при взрыве внутри подземного резервуара (типа ЖБР) (сценарий A6)

На рисунке 6 приложения № 1 к Руководству по безопасности приняты следующие условные вероятности событий:

при взрыве внутри резервуара образуются разлетающиеся элементы крыши резервуара (a) — 0,02;

зажигание нефти, нефтепродукта в резервуаре при отсутствии выброса из резервуара (c) — 0,2.

Рис. 7. «Дерево событий» при разрушении емкости под давлением (сценарий A7)

На рисунке 7 приложения № 1 к Руководству приняты следующие условные вероятности событий:

емкость сохраняет целостность после появления разрушения (a) — условная вероятность определяется согласно данным Руководства по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах»;

разрушение ниже уровня жидкости (b) — пропорционально отношению средней высоты уровня жидкости (взлива) к высоте емкости (если нет данных, принимают 0,8);

мгновенное воспламенение и образование горящих проливов/факелов (c) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы или газовой фазы;

образование дрейфующего облака ТВС (d) — для взрывопожароопасных жидкостей с давлением насыщенных паров в условиях содержания в емкости менее 10 кПа — 0, в остальных случаях — 1;

отсроченное воспламенение (e) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы или газовой фазы.

образование капельной взвеси (диспергированной струи) (h) — 0,7 (только в случае свищей).

образование BLEVE/огненного шара (j) — условная вероятность определяется согласно данным Руководства по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах» по полному и частичному разрушению емкости с учетом возможности образования BLEVE/огненного шара для конкретного продукта в емкости.

На рисунке 8 приложения № 1 к Руководству приняты следующие условные вероятности событий:

емкость сохраняет целостность после появления разрушения (a) — условная вероятность определяется согласно данным Руководства по безопасности «Методические основы анализа опасностей и оценки риска аварий на опасных производственных объектах»;

разрушение ниже уровня жидкости (b) — пропорционально отношению средней высоты уровня жидкости (взлива) к высоте емкости (если нет данных, принимают 0,8);

мгновенное воспламенение и образование горящих проливов/факелов (c) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы или газовой фазы;

образование дрейфующего облака ТВС (d) — для взрывопожароопасных жидкостей с давлением насыщенных паров в условиях содержания в емкости менее 10 кПа — 0, в остальных случаях — 1;

отсроченное воспламенение (e) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы или газовой фазы.

Рис. 8. «Дерево событий» при разрушении подземной емкости под давлением (сценарий A8)

Рис. 9. «Дерево событий» при аварии в насосных (сценарий A9)

На рисунке 9 приложения № 1 к Руководству приняты следующие условные вероятности событий:

возможность образования капельной смеси (a) — 0,7; мгновенное воспламенение и образование горящих проливов/факелов (b) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы;

образование дрейфующего облака ТВС (c) для взрывопожароопасных жидкостей с давлением насыщенных паров в условиях содержания в емкости менее 10 кПа (насосы на открытой площадке) или менее 3 кПа (насосы в помещении) — 0, в остальных случаях — 1;

отсроченное воспламенение (d) — 0,05 для истечения стабильной жидкой фазы; 0,2 для истечения нестабильной жидкой фазы.

Рис. 10. «Дерево событий» при разрушении устьевой арматуры нефтяных фонтанирующих скважин (сценарий A10)

Рис. 11. «Дерево событий» при разрушении устьевой арматуры нефтяных фонтанирующих скважин (сценарий A10)

Рис. 12. «Дерево событий» при авариях с оборудованием печей (сценарий A₁₂)

Таблица № 1

Вероятность возникновения источника воспламенения для сценариев на рисунках 10, 11, 12 приложения № 1 к Руководству

Таблица № 2

Вероятность воспламенения опасного вещества от появившегося энергетического (теплового) источника (различные типы опасного вещества) для сценариев на рисунках 10, 11, 12 приложения № 1 к Руководству

Таблица № 3

Вероятность эффекта BLEVE и огненного шара для сценариев на рисунках 10, 11, 12 приложения № 1 к Руководству

Приложение № 2 к Руководству по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи», утвержденному приказом Федеральной службы по экологическому, технологическому и атомному надзору от 10 января 2023 г. № 4

РАСЧЕТ КОЛИЧЕСТВА ВЕЩЕСТВ, УЧАСТВУЮЩИХ В АВАРИИ НА ЛИНЕЙНОЙ ЧАСТИ ПРОМЫСЛОВЫХ ТРУБОПРОВОДОВ

Расчет объема (массы) разлившейся нефти (пластовой жидкости, нефтепродукта)

Рассмотрим участок линейной части промыслового трубопровода (ЛЧ ПТ) протяженностью LН между насосными станциями (далее — НС1, НС2), на котором на расстоянии x от станций произошла аварийная утечка через эффективную площадь аварийного отверстия Sэфф.

Отметим, что как на участке LН, так и за его пределами (до НС1 и после НС2) участок может и не представлять собой изолированную систему, то есть в отдельных точках ЛЧ ПТ к нему могут подходить/отходить другие линейные участки. Эти подходящие/отходящие участки могут либо замыкаться на рассматриваемый участок (лупинги), либо осуществлять транспортирование жидкости в не связанные с данным линейным участком места (через отводы, ответвления).

Для штатного режима функционирования рассматриваемого участка расход жидкости составляет Q0. Также известны давления на входе и выходе отдельных линейных участков.

Общий объем вытекшей жидкости определяется процессами во всей разветвленной трубопроводной системе. Общий объем V определяют по формуле:

V = V₁ + V₂ + V₃, (2-1)

где: V₁ - объем жидкости, вытекшей в напорном режиме, с момента повреждения до остановки перекачки, м3; V₂ - объем жидкости, вытекшей в безнапорном режиме, с момента остановки перекачки до закрытия трубопроводной арматуры, м3; V3 — объем жидкости, вытекшей с момента закрытия трубопроводной арматуры до прекращения утечки (до момента прибытия аварийно-восстановительной бригады и ликвидации утечки или до полного опорожнения отсеченной части трубопровода), м3.

Объем V1, вытекший из участка ЛЧ ПТ за интервал времени \tau_{1} с момента возникновения аварии до остановки перекачки, определяют численным решением системы дифференциальных уравнений в частных производных, включающей законы сохранения массы, импульса и энергии потока ньютоновской жидкости:

уравнение неразрывности (уравнение сохранения массы):

A \cdot \frac{\partial\rho}{\partial\tau} + \frac{\partial(\rho \cdot A \cdot u)}{\partial x} = - M_{0} ; (2-2)

уравнение сохранения импульса:

A \cdot \frac{\partial(\rho \cdot u)}{\partial\tau} = - I_{0} - \frac{\partial(A \cdot \rho \cdot u^{2})}{\partial x} - A \cdot \frac{\partial P_{c}}{\partial x} - A \cdot \frac{\lambda_{\text{тр}}(Re)}{2 \cdot D} \cdot \rho \cdot u \cdot |u| - A \cdot g \cdot \rho \cdot \beta ; (2-3)

уравнение сохранения энергии:

A \cdot \frac{\partial}{\partial\tau} (\rho \cdot \epsilon) = - \frac{\partial}{\partial x} (\rho \cdot A \cdot u \cdot \epsilon) - E_{0} - \frac{A}{2 \cdot D} \cdot \lambda_{\text{тр}}(Re) \cdot \rho \cdot u^{3} + Q + A \cdot \Theta(T,T_{\text{ср}}) ; (2-4)

связь давления, плотности и температуры (уравнение состояния):

P_{c} - P_{0} = c^{2}(\rho - \rho_{0}) - c^{2}V(T_{0} - T) , (2-5)

где: x — расстояние от начала ПТ, м; Pc — осредненное по сечению давление жидкости, Па; P₀ - давление при нормальных условиях, Па (101325 Па); T0 — температура при нормальных условиях, К (293,15 К); \rho — осредненная по сечению плотность, кг/м3; \rho_{0} — плотность жидкости при нормальных условиях, кг/м3; u — осредненная по сечению скорость жидкости, м/с; \lambda_{\text{тр}}(Re) — коэффициент трения, зависящий от режима течения в трубе (от числа Рейнольдса Re = u \cdot D/\upsilon ), при необходимости в эту величину включаются и местные сопротивления на различных элементах (задвижках, клапанах и т.д.); A(x) — площадь поперечного сечения трубопровода, в общем случае переменная по трассе, м2; D(x) — номинальный диаметр ПТ, в общем случае переменный по трассе; M0(x) — удельная (на единицу длины трубы) интенсивность выброса жидкости из трубы на месте разрушения, кг/с/м; I0(x) — удельная (на единицу длины трубы) интенсивность потери импульса при выбросе жидкости из трубы на месте разрушения, кг/с2/м2; E0(x) — удельная (на единицу длины трубы) интенсивность потери внутренней энергии при выбросе жидкости из трубы на месте разрушения, Дж/с/м; g — ускорение свободного падения, м/с2; \beta — локальный угловой коэффициент трассы ПТ ( \beta = dz / dx); \epsilon — удельная внутренняя энергия, Дж/кг; \varsigma — коэффициент теплового объемного расширения, 1/К; Q(x) — удельная (на единицу длины трубы) интенсивность энергии, поступающей к транспортируемому продукту при его нагревании в нагревателях; \Theta(T,T_{\text{ср}}) — удельная (на единицу длины трубы) интенсивность теплообмена с окружающей средой; c — скорость распространения звука в транспортируемой продукции, м/с; z — нивелирная отметка трассы, м; \nu — кинематический коэффициент вязкости (v = \mu/\rho) , м2/с; \mu — динамический коэффициент вязкости жидкости (в общем случае зависящий от температуры транспортируемой среды), Н·с/м2.

В случае если температура в трубе остается постоянной (или меняется незначительно) на всем протяжении ПТ (T(x) — const, изотермическое течение), допускается решение только системы уравнений (2-2), (2-3) без уравнения (2-4).

Уравнение состояния (2-5) используют для сжатой среды (при Pc > P₀), при иных условиях (при растяжении) плотность полагают равной начальной плотности.

В случае отсутствия ветвлений участка систему уравнений (2-2) - (2-5) записывают только для одного участка. В случае наличия нескольких линейных участков, соединенных в единую систему (ответвления, лупинги и т.д.), систему уравнений (2-2) - (2-5) записывают для всех линейных участков, составляющих разветвленную трубопроводную систему и влияющих на массу выброса на месте аварии. При этом в уравнения (2-2) - (2-4) в правую часть добавляются слагаемые, описывающие дополнительное поступление или забор массы, импульса и энергии из отдельно взятого линейного участка в смежные участки; эти слагаемые аналогичны величинам M0(x), I0(x), E0(x).

Систему уравнений (2-2) - (2-5) дополняют начальными и граничными условиями.

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

Граничные условия выбирают следующим образом:

на входе ПТ производная давления полагается равной нулю, а скорость потока определяется с учетом этого давления по характеристике насоса Н — Q0 «напор-расход», также задается температура жидкости на входе;

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

После срабатывания запорной арматуры граничные условия на входе/выходе ПТ изменяются. Граничные условия соответствуют условию «жесткой стенки»: равенство нулю скорости на границах и равенство нулю первых производных по давлению и температуре.

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

Для определения величины используется зависимость Коулбрука-Уайта, связывающая коэффициент трения с числом Рейнольдса и характеристиками ПТ:

\frac{1}{\sqrt{\lambda_{\text{тр}}}} = - 2lg[\frac{2,51}{Re \sqrt{\lambda_{\text{тр}}}} + \frac{z_{\text{шер}}}{3,71D}] ; (2-6)

где Zшер — шероховатость внутренней поверхности ПТ.

Соотношение (2-6) представляет собой трансцендентное уравнение, решая которое можно определить \lambda_{\text{тр}}(Re) .

Помимо соотношения (2-6) для определения величины могут использоваться иные обоснованные соотношения:

\lambda_{\text{тр}} = - \frac{64}{Re} \text{при }Re<2000 ; (2-7)

\lambda_{\text{тр}} = (0,16 \cdot Re - 13) \cdot 10^{- 4}\text{при }2000 \leq Re \leq 2800 ; (2-8)

\lambda_{\text{тр}} = \frac{0,3164}{Re \cdot 0,25} \text{при }2800 \leq Re \leq Re_{1} ; (2-9)

\lambda_{\text{тр}} = B + \frac{1,7}{Re \cdot 0,5} \text{при }Re_{1} \leq Re \leq Re_{2} . (2-10)

Предельные значения Re₁, Re₂ и значения B приведены в таблице № 1 приложения № 2 к Руководству.

Таблица № 1

Предельные значения Re₁, Re₂ и значения B

При числах Рейнольдса, больших, чем указанные в таблице № 1 приложения № 2 к Руководству значений Re2 (в квадратичной зоне), коэффициент гидравлического сопротивления остается постоянным и равным значению \lambda_{\text{тр}} , рассчитанному по формуле (2-10) при Re = Re2.

Для определения величины Q используют зависимость:

\Theta(T,T_{\text{ср}}) = \frac{4\alpha_{\text{т}}}{D} (T_{\text{ср}} - T) , (2-11)

где: T — температура нефти (нефтепродукта), °C; Tср — температура окружающей среды, °C; \alpha_{\text{т}} — коэффициент теплопередачи нефти (нефтепродукта) с окружающей средой, определяемый по формуле:

(2-12)

где: Hпр — приведенная толщина, мм; Hпр — H + Hэ; Hэ — эффективная толщина, мм; Hэ = Hсн \lambda_{\text{гр}}/\lambda_{\text{сн}} ; Hсн — толщина стенки, мм; Nuгр — число Нуссельта (Nu_{\text{гр}} = \alpha_{\text{ввод}}D_{N}/\lambda_{\text{гр}}) ; \lambda_{\text{гр}},\lambda_{\text{сн}} — коэффициенты теплопроводности; \lambda_{\text{гр}} = 1,4\text{Вт} \cdot \text{К}/\text{м} ; \lambda_{\text{сн}} = 0,465\text{Вт} \cdot \text{К}/\text{м} ; \alpha_{\text{ввод}} — коэффициент теплопередачи воздуха (\alpha_{\text{ввод}} = 11,6\text{Вт} \cdot \text{К}/\text{м}^{2}) .

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

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

U_{0} = \sqrt{2 \frac{P - P_{\text{нар}}}{\rho}} , (2-13)

где Pнар — давление снаружи ПТ, Па.

Для сухопутных участков Pнар = 101325 Па, для подводных ПТ величину Pнар определяют как сумму атмосферного давления и давления столба жидкости над отверстием разгерметизации.

Соответственно, поток массы через отверстие задают выражением:

M_{0} = \alpha \cdot S_{j} \cdot U_{0} \cdot \rho , (2-14)

где: \alpha — коэффициент, который принимает максимально возможное значение, равное 0,6; Sj — площадь отверстия разгерметизации, м2.

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

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

Для вывода интегральных напорно-расходных характеристик НС используют формулу:

H = a - b \cdot Q_{0}^{2} , (2-15)

где a, b — экспериментально определенные коэффициенты штатного режима работы насосов.

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

до спада давления в трубопроводе (в частности, до установления вакуумметрических давлений в самых высоких точках трассы в каждом из прилегающих к месту аварии участков трубопровода) решением системы уравнений (2-2) - (2-14). В этом случае расчет объемов вытекшей жидкости V2 является продолжением расчета объемов вытекшей жидкости V1 с изменением граничных условий (остановкой насосов на входе трубопроводной системы);

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

Алгоритм расчета объема жидкости V2 в этом случае зависит от размеров отверстия разгерметизации. При свищах размер отверстия настолько мал, что существенного движения среды в трубе не наблюдается. Поэтому при расчете интенсивности истечения можно, пренебрегая столь малым движением, жидкость в трубопроводе считать покоящейся, а зеркало жидкости в каждом из Nст участков трубопроводной системы будет находиться на одном уровне (zзер). Давление в трубопроводе будет определяться гидростатикой:

\text{р}_{\text{внутр}} = \text{р}_{\text{у}} + \rho g(z_{\text{зер}} - z) , (2-16)

где: zзер — уровень (нивелирная отметка трассы), на котором находится нефть, нефтепродукт (зеркало жидкости), м; z — нивелирная отметка трассы, м; pвнутр — внутреннее давление в трубопроводе, Па; pу — вакуумметрическое давление паров нефти, Па; \rho — осредненная по сечению плотность, кг/м3; g — ускорение свободного падения, м/с2.

Давление на месте разрушения (высотная отметка z*) составит:

\text{р}_{\text{внутр}} = \text{р}_{\text{у}} + \rho g(z_{\text{зер}} - z*) , (2-17)

где z* — уровень (нивелирная отметка трассы), на котором находится место разрушения, м.

Расход жидкости через свищ определяют следующим образом:

M = \alpha S_{1}\rho \sqrt{2 \frac{p_{\text{внутр}} - p_{\text{нар}}}{\rho}} , (2-18)

где: \alpha — коэффициент, который принимает максимально возможное значение, равное 0,6; Si — площадь свища, м2; \rho — осредненная по сечению плотность, кг/м3; pнар — наружное давление в окружающей среде на месте разрушения, Па.

Положение высотной отметки зеркала жидкости zзер по мере истечения жидкости меняется (zзер(t)). В начальный момент времени площадь зеркала находится на уровне максимальной высоты из всех трубопроводов. В последующие моменты времени высота зеркала жидкости уменьшается, при этом по достижении высоты максимально высоких точек в других трубопроводных ответвлениях в этих ответвлениях также будут появляться свои поверхности жидкости. В конечном итоге в системе сформируется столько поверхностей, сколько в ней ответвлений. Во всех ответвлениях высота зеркала нефти (нефтепродукта) совпадает. По мере стока жидкости происходит не только постепенное снижение высоты зеркал нефти (нефтепродукта) zзер(t), но и перемещение их вдоль трубопровода (x1(t), x2(t) и т.д.). Такое перемещение по длине происходит как непрерывно, так и скачками. Скачкообразное изменение xi(t) происходит, когда на пути зеркала встречается V-образный спуск-подъем и высота зеркала сравнивается с высотой лежащего по ходу слива жидкости локальным максимумом. Скачок происходит на величину расстояния, которое разделяет локальный максимум и точку на спуске с той же высотой, что и локальный максимум. V-образный профиль между этими точками остается заполнен жидкостью, и слив далее будет происходить из участка, расположенного после локального максимума.

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

(\frac{p_{\text{у}}}{\rho g} + z_{1}(t)) - (\frac{p_{\text{внутр}}}{\rho g} + z*) = \lambda_{\text{тр}}(Re_{1}) \frac{(x* - x_{1})}{D} \frac{u_{1}^{2}}{2g} ; (2-19)

(\frac{p_{\text{у}}}{\rho g} + z_{2}(t)) - (\frac{p_{\text{внутр}}}{\rho g} + z*) = \lambda_{\text{тр}}(Re_{2}) \frac{(x_{2} - x*)}{D} \frac{u_{2}^{2}}{2g} ; (2-20)

\frac{u_{1}}{4} \pi D_{1}^{2} - \frac{u_{2}}{4} \pi D_{2}^{2} = \alpha S_{j} \sqrt{2 \frac{p_{\text{внутр}} - p_{\text{нар}}}{\rho}} ; (2-21)

\frac{dx_{1}}{dt} = u_{1} , \frac{dx_{2}}{dt} = - u_{2} , (2-22)

где: t — время, с; z* — уровень (нивелирная отметка трассы), на котором находится место разрушения, м; z1(t) — уровень (нивелирная отметка трассы), на котором находится перемещающееся зеркало жидкости на участке до места разрушения, м; z2(t) — уровень (нивелирная отметка трассы), на котором находится перемещающееся зеркало жидкости на участке после места разрушения, м; pвнутр — внутреннее давление в трубопроводе на месте разрушения, Па; pнар — наружное давление в окружающей среде на месте разрушения, Па; pу — вакуумметрическое давление паров жидкости, Па; \rho — осредненная по сечению плотность, кг/м3; g — ускорение свободного падения, м/с2. \lambda_{\text{тр}}(Re_{1}) — коэффициент трения, зависящий от режима течения в трубе, для участка до места разрушения, где число Рейнольдса равно Re1; \lambda_{\text{тр}}(Re_{2}) — коэффициент трения, зависящий от режима течения в трубе, для участка после места разрушения, где число Рейнольдса равно Re2; u1 — скорость движения жидкости к месту аварии (участок до места разрушения) на стадии самотечного истечения, м/с; u2 — скорость движения жидкости к месту аварии (участок после места разрушения) на стадии самотечного истечения, м/с; x* — координата по трассе места разрушения, м; x1 — координата перемещающейся поверхности (зеркала жидкости) в участке до места разрушения, м; x2 — координата перемещающейся поверхности (зеркала жидкости) в участке после места разрушения, м; D1 — внутренний диаметр ПТ до места разрушения, м; D2 — внутренний диаметр ПТ после места разрушения, м.

Когда на месте разрушения избыточное давление отсутствует (pвнутр = p0, например, при гильотинном разрыве), скорость выброса будет определяться потоками жидкости к месту аварии. При этом также используют условия непрерывности давления и сохранения потока массы в местах изменения скорости потока (места разрыва, изменения диаметра, ветвления и т.д.). Например, для стока из двух участков необходимо решить систему из следующих уравнений:

(\frac{p_{\text{у}}}{\rho g} + z_{1}(t)) - (\frac{p_{\text{внурт}}}{\rho g} + z*) = \lambda_{\text{тр}}(Re_{1}) \frac{(x* - x_{1})}{D} \frac{u_{1}^{2}}{2g} ; (2-23)

(\frac{p_{\text{у}}}{\rho g} + z_{2}(t)) - (\frac{p_{\text{внурт}}}{\rho g} + z*) = \lambda_{\text{тр}}(Re_{2}) \frac{(x_{2} - x*)}{D} \frac{u_{2}^{2}}{2g} ; (2-24)

\frac{u_{1}}{4} \pi D_{1}^{2} - \frac{u_{2}}{4} \pi D_{2}^{2} = \alpha S_{j} \sqrt{2 \frac{p_{\text{внутр}} - p_{\text{нар}}}{\rho}} ; (2-25)

\frac{dx_{1}}{dt} = u_{1} , \frac{dx_{2}}{dt} = - u_{2} . (2-26)

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

Расчет аварийных выбросов на промысловых трубопроводах со сжиженным углеводородным газом

При разгерметизации ПТ с газовым конденсатом, продукцией скважин, сжиженными углеводородными газами (далее — СУГ) и иными продуктами, способными образовывать при аварийном падении давления облака ТВС, по трубопроводу от места разрыва распространяются волны разгрузки. На участке до места разрыва (по потоку) среда ускоряется, а на участке после места разрыва — замедляется, вплоть до образования обратного потока перекачиваемой среды.

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

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

Описанные эффекты имеют место как при транспортировании СУГ, так и при транспортировании изначально газокапельной среды.

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

Таким образом, при исследовании аварийного истечения жидких продуктов из ПТ с СУГ можно выделить две основные стадии:

истечение в напорном режиме, когда в трубопровод подается и из трубопровода отбирается продукт;

самотечный режим истечения.

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

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

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

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

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

Движение среды по каждому участку описывается следующими одномерными нестационарными уравнениями для осредненного по сечению течения:

закон сохранения массы:

A \cdot \frac{\partial\rho}{\partial\tau} + \frac{\partial(\rho \cdot A \cdot w)}{\partial x} = - M_{0} , (2-27)

где: x — пространственная координата вдоль оси трубопровода, м; \tau — время, с; \rho — общая плотность смеси, кг/м3; w — скорость движения среды, м/с; A — площадь сечения трубопровода, м2; M0(x,t) — расход выброса, кг/с/м;

закон сохранения отдельных компонент/фаз:

A \cdot \frac{\partial(\rho \cdot Y_{m})}{\partial\tau} + \frac{\partial(\rho \cdot A \cdot Y_{m} \cdot w)}{\partial x} = - Y_{m}M_{0} - Y_{m} , (2-28)

где Y_{m} — скорость исчезновения/появления m-й компоненты/фазы в результате испарения/конденсации;

закон сохранения импульса:

A \cdot \frac{\partial(\rho \cdot w)}{\partial\tau} + \frac{\partial(\rho \cdot A \cdot w^{2})}{\partial x} = - I_{0} - A \cdot \frac{\partial p}{\partial x} - A \cdot \frac{\lambda_{\text{тр}}(Re)}{2D} \cdot \rho \cdot w \cdot |w| - A \cdot \rho \cdot g \cdot \beta , (2-29)

где: p — давление в системе, Па; I0(x,t) — потери импульса в системе при выбросе среды, кг/с2; g — ускорение свободного падения (9,81 м/с2); z — высотная отметка трубопровода, м; \pi — число Пифагора (3,14159...); \lambda — коэффициент трения;

закон сохранения энергии:

A \cdot \frac{\partial}{\partial\tau} (\rho \cdot \epsilon) + \frac{\partial}{\partial x} (\rho \cdot A \cdot u \cdot \epsilon) = - E_{0} - \frac{A}{2 \cdot D} \cdot \lambda_{\text{тр}}(Re) \cdot \rho \cdot w^{3} + Q + A \cdot \Theta(T,T_{sur}) ; (2-30)

где: \epsilon — удельная внутренняя энергия (на единицу массы), Дж/кг; E0(x,t) — потери энергии в системе при выбросе среды, Дж/(м·с); Q — удельная (на единицу массы) скорость энерговыделения/энергополощения за счет фазовых переходов, протекающих в системе, либо за счет подогрева смеси в отдельных точках трассы, Дж/(м3·с); \Theta(T,T_{sur}) — теплообмен с окружающей средой через стенки трубы, Дж/(м·с).

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

В случае отсутствия ветвлений систему уравнений (2-27) - (2-30) решают только для одного участка. В случае наличия нескольких линейных участков, соединенных в единую систему (ответвления, лупинги и т.д.), систему уравнений (2-27) - (2-30) решают для всех линейных участков, составляющих разветвленную трубопроводную систему и влияющих на массу выброса на месте аварии. При этом в уравнения (2-27) - (2-30) в правую часть добавляются слагаемые, описывающие дополнительное поступление/забор массы, массы отдельных компонент/фаз, импульса и энергии из отдельно взятого линейного участка в смежные участки; эти слагаемые аналогичны величинам M0(x,t), I0(x,t), E0(x,t).

Представленная система уравнений учитывает широкий перечень факторов и процессов:

изменение сечения по длине трубопровода;

изменение параметров среды (плотность, концентрация, скорость и энергия) в каждой точке пространства (первые слагаемые левых частей в уравнениях (2-27) - (2-30);

перенос в поле скорости (конвективный перенос) массы в целом, массы отдельных компонент смеси, импульса и энергии (вторые слагаемые левых частей в уравнениях (2-27) - (2-30);

потери массы (в том числе отдельных компонент), импульса и энергии в системе при выбросе среды на местах разрушения (первые слагаемые правых частей в уравнениях (2-27) - (2-30);

действие сил давления (второе слагаемое правой части (2-29);

действие силы тяжести (четвертое слагаемое правой части (2-29);

потери на трение на стенках трубы и на фитинге трубопровода (третье слагаемое правой части (2-29) и второе слагаемое правой части (2-30);

энерговыделение/энергопоглощение за счет реакций/фазовых переходов, протекающих в транспортируемой среде, в том числе за счет принудительного внешнего нагрева транспортируемого продукта (третье слагаемое правой части (2-30);

теплообмен продукта, транспортируемого по трубопроводу, с окружающей средой через стенки трубопровода (четвертое слагаемое правой части (2-30).

Система уравнений (2-27) - (2-30) должна быть дополнена следующими соотношениями:

уравнениями состояния жидкой и газокапельной сред; для случая течения только жидкости в качестве уравнения состояния используется соотношение (2-5); для случая наличия двух фаз для газовой фазы используется уравнение, описывающее зависимость давления газа от плотности и температуры (например, уравнение состояния идеального газа), а жидкая фаза может считаться несжимаемой;

откорректированным соотношением для \lambda_{\text{тр}} , учитывающим двухфазность потока (например, зависимость Локкарта-Мартинелли);

соотношением для расчета скоростей фазовых переходов Y_{m} (долей газовых/жидких фракций); в простейшем случае допускается использовать предположение о термодинамическом давлении фаз; в этом случае в качестве уравнения состояния используются результаты термодинамического расчета задачи об отыскании параметров в системе при известных плотности (берется из решения уравнения (2-27)), энергии (берется из решения уравнения (2-30)) и составе (берется из решения уравнения (2-28)). Именно в ходе термодинамического расчета отыскивается давление, температура и доля вскипания в соответствующем сечении трубопровода.

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

Система уравнений (2-27) - (2-30) решается численно на ЭВМ с использованием разностных методов необходимой точности.

Решение системы (2-27) - (2-30) позволяет отыскать для каждого места утечки величину M0(t), а также пространственно-временное распределение всех параметров, обусловливающих M0(t): скорость продукта в трубопроводе, давление, температуру, вскипевшую долю и состав СУГ и т.д. Именно расход продукта M0(t) и его параметры на месте выброса позволяют проводить в дальнейшем необходимые расчеты по последствиям выброса.

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

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

Расчет аварийных выбросов на промысловых газопроводах

Для инженерных расчетов аварийных выбросов на промысловых газопроводах применяются соотношения, изложенные в приложении № 7 к Руководству по безопасности «Методика оценки риска аварий на опасных производственных объектах магистрального трубопроводного транспорта газа».

Приложение № 3 к Руководству по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи», утвержденному приказом Федеральной службы по экологическому, технологическому и атомному надзору от 10 января 2023 г. № 4

МЕТОДИКИ РАСЧЕТА ИНТЕНСИВНОСТИ ИСТЕЧЕНИЯ ПРИ ФОНТАНИРОВАНИИ СКВАЖИН

Методика расчета интенсивности истечения газа при фонтанировании скважин

Исходные данные:

\lambda — коэффициент гидравлического сопротивления, б/р;

g — ускорение свободного падения, м/с2;

\rho_{\text{н}} , Tн, Pн — плотность, температура, давление газа при нормальных условиях;

Tср — средняя температура в скважине, К;

Zср — среднее значение коэффициента сжимаемости;

Pз — давление в скважине напротив работающего интервала;

a, b — коэффициенты фильтрационного сопротивления;

Pпл — пластовое давление;

Pа — атмосферное давление (рассматривается открытый фонтан).

Геометрия секции скважины:

li — длина секции, м;

dнар i — наружный диаметр кольцевого пространства, м;

dвн i — внутренний диаметр кольцевого пространства, м;

\alpha_{zi} — зенитный угол, град;

\eta — вязкость газа в пластовых условиях, сП;

с — скорость звука, м/с;

к — проницаемость пласта в окрестности скважины, Д;

m — пористость пласта в окрестности скважины, б/р;

hпл — эффективная толщина пласта, м;

zср — среднее значение коэффициента сжимаемости.

Искомые параметры:

G1(t) — массовый расход из скважины, кг/с.

Порядок расчета

Модель установившегося истечения из скважины

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

\rho \cdot Q = const , (3-1)

где: Q — объемный расход через сечение; \rho — средняя плотность газа в сечении.

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

l_{i},d_{\text{нар }i},d_{\text{вн }i},\alpha_{zi},i = 1,2,3,...,N , (3-2)

где: li — длина секции; dнар i — наружный диаметр кольцевого сечения; dвн i — внутренний диаметр; \alpha_{zi} — зенитный угол (угол между направлением оси секции и вертикалью).

В пределах секции справедливо уравнение количества движения:

\frac{d(\rho \cdot W^{2})}{dl} + \frac{dP}{dl} = - \frac{\lambda_{i}\rho \cdot W|W|}{2d_{i}} + \rho \cdot g \cdot cos\alpha_{zi} , (3-3)

где: dj — гидравлический диаметр, вычисляемый по формуле:

di = dнар i — dвн i; (3-4)

l — расстояние от устья (при фонтанировании скорость в выражении (3-3) отрицательна); W — скорость газа, м/с; P — давление газа, Па.

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

Уравнение состояния газа записывается в обычном виде:

\rho = \rho_{\text{н}} \cdot \frac{\text{Т}_{\text{н}}}{\text{Т}} \cdot \frac{1}{Z} \cdot \frac{P}{P_{\text{н}}} , (3-5)

где \rho_{\text{н}} , Tн, Pн — плотность, температура, давление газа при нормальных условиях.

Потери давления в пласте при стационарной фильтрации описываются уравнением:

P_{\text{з}}^{2} = \text{Р}_{\text{пл}}^{2} - \alpha \cdot Q_{\text{н}} - \alpha \cdot Q_{\text{н}}^{2} , (3-6)

где: Pз — давление в скважине напротив работающего интервала; a, b — коэффициенты фильтрационного сопротивления; Pпл — пластовое давление.

Уравнение (3-6) можно рассматривать как граничное условие для системы (3-1) - (3-3). Условие на устье имеет вид:

где: Pу, Wу, C — давление, скорость газа и скорость звука на устье; Pа — атмосферное давление (рассматривается открытый фонтан).

Задача решается в предположении постоянства температуры и коэффициента сжимаемости:

Z = const; T = const. (3-8)

Задача состоит в определении дебита фонтана при заданных параметрах пласта (пластовое давление, коэффициенты фильтрационного сопротивления), геометрии ствола и параметров уравнения состояния (3-1), (3-8). Рекомендуется использовать в качестве средних значений в (3-8) среднеарифметические значения температуры и коэффициента сжимаемости для пластовых и устьевых условий. Решение проводится методом деления отрезка пополам. В качестве нижней границы корня принимается нулевое значение дебита. Верхняя граница определяется путем расчета забойного давления для нескольких последовательно возрастающих значений QН.

Модель залпового выброса из скважины

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

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

\frac{\partial\rho}{\partial t} + \frac{\partial(\rho W)}{\partial t} = 0 ; (3-9)

\rho \frac{dW}{dt} + \frac{\partial P}{dl} + gp[- 1 + \frac{|W| \cdot W}{W_{gid}^{2}}] = 0 ; (3-10)

W_{gid} = \sqrt{\frac{2gD_{gid}}{\lambda}} . (3-11)

Уравнение состояния принимается в виде

P = \rho \cdot z_{\text{ср}} \cdot R \cdot T_{\text{ср}} , (3-12)

где: \frac{d}{dt} = \frac{\partial}{\partial t} + W \frac{\partial}{\partial l} , t — время; l — длина вдоль оси ствола; g — ускорение свободного падения; P, W, \rho — давление, скорость и плотность газа; \lambda — коэффициент гидравлического сопротивления (принимается постоянным).

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

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

- \frac{\partial P_{0}}{\partial l} + g\rho = 0 . (3-13)

На забое давление в скважине совпадает с пластовым давлением:

P0(lb) = Pпл, (3-14)

где: lb — координата кровли проявляющего интервала; Pпл — пластовое давление в местах расположения скважин.

Пусть в момент времени t = 0 происходит мгновенная разгерметизация устья. Скорость в выходном сечении будет равна местной скорости звука, а текущий дебит вычисляется по формуле:

Q_{h} = C_{D} \cdot W_{S} \cdot F_{h} \cdot \rho_{h} , (3-15)

где: Qh — дебит фонтана; Ws — местная скорость звука; Fh — площадь выходного сечения; \rho_{h} — плотность флюида в выходном сечении; CD — коэффициент расхода, зависящий от формы выходного сечения.

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

Таким образом, в приустьевой части ствола выполняются уравнения:

\rho FW = Q_{h}(t) = const ; (3-16)

- \frac{dP}{dl} - \frac{d(\rho W^{2})}{dl} + g\rho[1 + (\frac{W}{W_{gid}})^{2}] = 0 . (3-17)

На подвижной границе выполняется условие:

P(lf, t) = P0(lf), (3-18)

где lf — текущее положение фронта.

Зная распределение давления по стволу, можно определить массу газа, находящегося в момент t в скважине. Из сказанного следует, что масса полностью определяется положением фронта: M(t) = M(lf). Из условия материального баланса, примененного ко всему стволу, следует уравнение перемещения фронта:

M^{l} \cdot \frac{dl_{f}}{dt} = Q_{h} . (3-19)

После того как волна достигла забоя, приходит в движение флюид в пласте. Предполагая течение симметричным относительно оси скважины, обозначим через Rf радиус границы (радиус депрессионной воронки), отделяющей область неподвижного флюида от прискважинной области, в которой поток стационарен, и его дебит равен мгновенному дебиту фонтана. Для вычисления Rf используется уравнение, аналогичное (3-19), в котором под M понимается масса газа в стволе и круговой области пласта, радиус которой Ro выбран так, что на рассматриваемом интервале времени Rf < Ro. Чтобы определить M(Rf), рассматривается задача о стационарном течении в системе «скважина — пласт», удовлетворяющем условию на устье (3-15) и условию на подвижном контуре:

P(Rf, t) = Pпл. (3-20)

Кроме того, выполняется условие сопряжения — непрерывное изменение давления и массового расхода при переходе от пласта к скважине.

Методика расчета интенсивности истечения жидкости при фонтанировании скважин

Скважина представляется как совокупность цилиндрических каналов переменного диаметра, состыкованных последовательно торец к торцу.

Скважина может иметь произвольный угол наклона к вертикали на различных своих участках.

Рассматривается случай стационарного истечения.

Движение флюида в такой скважине описывается уравнением сохранения импульса в предположении изотермичности потока:

\frac{d(\rho \cdot U^{2})}{dx} + \frac{dP}{dx} = - \frac{\lambda\rho \cdot U|U|}{2d(x)} + \rho \cdot g \cdot cos(\alpha(x)) , (3-21)

где: \rho — плотность флюида; U — скорость флюида; x — расстояние от конечной точки скважины; P — давление по длине скважины; \alpha — угол отклонения скважины от горизонтали; \lambda — коэффициент трения; d(x) — внутренний диаметр канала.

Граничными условиями для данного уравнения служат давления в пласте и окружающей среде (1 атм).

Уравнение движения флюида замыкается двумя соотношениями:

уравнением состояния;

соотношением для определения трения.

Приложение № 4 к Руководству по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи», утвержденному приказом Федеральной службы по экологическому, технологическому и атомному надзору от 10 января 2023 г. № 4

ОСНОВНЫЕ ПОКАЗАТЕЛИ РИСКА АВАРИЙ НА ОПО НГД

Таблица № 1

Показатели риска аварий на промысловых трубопроводах

Таблица № 2

Показатели риска аварий на площадочных ОПО и их составляющих

Приложение № 5 к Руководству по безопасности «Методика анализа риска аварий на опасных производственных объектах нефтегазодобычи», утвержденному приказом Федеральной службы по экологическому, технологическому и атомному надзору от 10 января 2023 г. № 4

РЕКОМЕНДОВАННЫЕ ЧАСТОТЫ РАЗГЕРМЕТИЗАЦИИ И ФОНТАНИРОВАНИЯ СКВАЖИН

Частоты аварийной разгерметизации газовых, нефтяных и газоконденсатных скважин, включая морские, рекомендуется принимать согласно таблице № 1 приложения № 5 к Руководству.

Частоты аварийной разгерметизации промысловых трубопроводов рекомендуется принимать согласно сведениям Руководства по безопасности «Методика оценки риска аварий на технологических трубопроводах, связанных с перемещением взрывопожароопасных газов», утвержденного приказом Ростехнадзора от 28 ноября 2022 г. № 410 и Руководства по безопасности «Методика оценки риска аварий на технологических трубопроводах, связанных с перемещением взрывопожароопасных жидкостей», утвержденному приказом Ростехнадзора от 28 ноября 2022 г. № 410.

Таблица № 1

Частоты разгерметизации газовых, нефтяных и газоконденсатных скважин, включая морские

Пример:

На проектируемой кустовой площадке предусматривается размещение трех газовых скважин с аномально высоким пластовым давлением (АВПД). Предполагаемый период эксплуатации скважин 20 лет.

Общая продолжительность строительных работ на кустовой площадке, включая бурение скважин, — 3 месяца.

Решение:

Используя данные таблицы № 1 приложения № 5 к Руководству, получаем:

Частота выбросов (включая фонтанирование) для добывающих газовых скважин:

на стадии эксплуатации — 7,2 · 10-5 + 8,8 · 10-5 = 1,6 ·10-4 на одну скважину в год;

на стадии капитального ремонта скважины — 4,0 · 10-4 + 7,6 · 10-4 = 1,16 · 10-3 на операцию;

при бурении (для скважин с АВПД) — 2,6 · 10-4 + 2,5 · 10-3 = 2,76 · 10-3 на пробуренную скважину;

на стадии заканчивания скважины — 2,8 · 10-4 + 4,8 · 10-4 = 7,6 · 10-4 на операцию.

Частота выброса из скважин в период их строительства с учетом принятой продолжительности строительных работ:

(2,76 · 10-3 + 7,6 · 10-4) · 3 / 12 = 8,8 · 10-4 год-1.

С учетом периода эксплуатации скважин 20 лет и допущения о периодичности проведения КРС и его средней продолжительности расчетная частота выброса из всех скважин куста в период их эксплуатации может быть оценена величиной:

(1,6 · 10-4 · 20 + 4 · (30 / 365) · 1,16 · 10-4) / 20 = 1,62 · 10-4 год-1.