Актуальность проблемы. Топливно-энергетический комплекс (ТЭК) является важнейшей отраслью промышленности — базисом экономики промышленно развитых стран. Одним из основных компонентов ТЭК являются системы газопроводов высокого давления, значительную часть которых составляют магистральные газопроводы (МГ), транспортирующие природный газ от места добычи до потребителя.
Единая система газоснабжения Российской Федерации (ЕСГ) имеет более 150 тысяч километров МГ и отводов, из них диаметром 1420мм — около 50 тысяч километров, природный газ по которым транспортируется 251 компрессорной станицей (КС) с 689 компрессорными цехами (КЦ) общей установленной мощностью более 42млн.кВт. Протяженность МГ от промыслов северных районов Тюменской области до наиболее удаленных импортеров, таких как Франция и Италия, составляет более 5 тысяч километров [1]. Все это говорит об огромных размерах и сложности системы ЕСГ.
В последнее время в мире отмечается снижение надежности работы МГ. В первую очередь это обусловлено тем, что основной парк МГ составляют трубопроводы, имеющие срок эксплуатации свыше 25 лет. Кроме этого, параметры режимов работы МГ очень часто отличаются от проектных значений. Все вышесказанное обуславливает участившиеся случаи разрушения отдельных участков МГ, сопровождающиеся массовыми выбросами природного газа и сопутствующих фракций в окружающую среду. В процессе аварийного выброса природного газа в месте разрушения газопровода происходит его интенсивное перемешивание с атмосферным воздухом. Образовавшаяся метановоздушная смесь является легко воспламеняемой. По статистике аварий на российских МГ свыше 80% процентов разрывов МГ сопровождаются интенсивными пожарами [2].
Многие газопроводы высокого давления проходят по территории или в непосредственной близости от населенных пунктов. Пожары, сопровождающие разрывы таких трубопроводов, угрожают, как промышленным и жилым зданиям и сооружениям, так и жизни людей. Так, например, при аварии на МГ «Уренгой-Петровск» ООО «Баштрансгаз», произошедшей 09 апреля 2003 года, разрушился участок МГ с возгоранием транспортируемого природного газа. В результате аварии в близлежащей деревне сгорело три дома [3].
К сожалению, как констатируют руководители Госгортехнадзора Российской Федерации в статье [2], только на магистральных трубопроводах с 1992 по 2001 годы произошло 545 аварий. Таким образом, среднегодовой уровень аварийности составляет 50−60 аварий в год и в целом не имеет устойчивой тенденции к снижению [2]. На фоне общего роста числа аварий возрастает и число разрывов МГ. В период с октября 2001 года по февраль 2002 года на предприятиях ОАО «Газпром» зарегистрировано 5 разрывов газопроводов высокого давления. Из них 4 аварии сопровождались возгоранием транспортируемого природного газа. В период с мая 2003 года по июнь 2004 года на МГ, транспортирующих природный газ, произошло 18 аварий, из них 15 аварий сопровождались возгоранием метановоздушной смеси [3−10] (рис.В.1 — В. З).
Помимо пожаров, выбросы в атмосферу газовых смесей (транспортируемых по трубопроводам высокого давления) наносят окружающей среде ощутимый экологический урон. Так, например, при аварии на МГ «Ямбург-Западная граница» ООО «Тюментрансгаз» 11 апреля 2003 года произошло разрушение трубопровода без возгорания транспортируемого газа, сопровождавшееся большим выбросом природного газа в атмосферу [2].
Рис.В.1, Фотографии пожара на месте разрыва МГ в черте города.
Рис.В.2. Фотографии места разрыва МГ в Западной Европе, сопровождавшегося интенсивным пожаром.
Кроме экологического фактора, авария на МГ сопровождается тяжелыми экономическими последствиями. Это не только выплата штрафов и компенсаций за причиненный ущерб природе и социальной сфере жизнедеятельности людей, но и большие затраты на восстановительные работы.
Рис.В. З. Фотография места разрыва российского МГ.
Одним из важных направлений повышения пожарной и промышленной безопасности МГ является уменьшение продолжительности и интенсивности пожара, сопровождающего разрыв на МГ. Продолжительность такого пожара напрямую связана с оперативностью обнаружения разрыва МГ. Кроме того, от оперативности часто зависит и сама возможность локализации аварии, т.к. линейные краны, отсекающие аварийную нитку многониточной линейной части магистрального газопровода (ЛЧМГ), при длительном развитии аварии могут потерять работоспособность. Этим обусловлена актуальность проблемы оперативного обнаружения разрывов МГ.
Современный научный подход к повышению безопасности, экологичности и эффективности сложных газопроводных систем ТЭК заключается в построении и численном анализе комплекса взаимосвязанных математических моделей трубопроводных сетей, в основу которых положены базовые уравнения механики сплошных сред [11−17]. Одно из центральных мест в данном комплексе занимает математическая модель течения гомогенных газовых смесей по системам длинных разветвленных трубопроводов и компрессорных станций, лежащая в основе построения высокоточных газодинамических симуляторов (ГДС) [11,12,15]. ГДС позволяют:
• проводить численный анализ и оптимизацию процессов транспортирования газа по трубопроводным системам;
• получать исходные данные для численного прочностного анализа трубопроводных систем [18,19];
• получать исходные данные для проведения численного моделирования пожаров на МГ [20−26] и т. д.
Проведение углубленного анализа физических процессов транспортирования газа по трубопроводным сетям промышленных энергообъектов в номинальных режимах и при авариях с использованием ГДС предоставляет возможности для решения задачи оперативного обнаружения разрывов МГ.
Данная диссертационная работа направлена на решение проблемы оперативного обнаружения разрывов МГ с использованием ГДС.
Научная новизна работы состоит в следующем:
1. Разработана новая технология повышения пожарной и промышленной безопасности многониточных магистральных газопроводов, основанная на методе оперативного обнаружения и локализации их разрывов с использованием высокоточных газодинамических симуляторов. Оперативное обнаружение разрывов требует построения ГДС по технологии, впервые предложенной В. Е. Селезневым [11] и развитой С. Н. Пряловым [12], В. В. Киселевым [13] при непосредственном участии автора диссертации. Построенный ГДС используется для высокоточного численного анализа изменения основных параметров транспортирования природного газа (давлений, температур и массовых расходов) по длине трубопроводов при их эксплуатации в номинальных и аварийных режимах. Обнаружение аварийного разрыва многониточной ЛЧМГ производится в результате автоматического сравнительного анализа невязок между рассчитанными (с использованием ГДС) и измеренными градиентами давлений в транспортируемой среде на границах контролируемых участков трубопроводной сети. Локализация аварийного разрыва трубопровода производится в результате сравнительного анализа измеренных и рассчитанных с помощью ГДС нестационарных значений перепада давления на границах поврежденного трубопровода.
2. Впервые предложен и обоснован метод моделирования крановых площадок (крановых узлов), оборудованных межниточными перемычками, для обеспечения оперативности обнаружения разрывов многониточных линейных частей магистральных газопроводов. Он направлен на повышение скорости при одновременном сохранении точности численного анализа газодинамических режимов функционирования ЛЧМГ с использованием ГДС. Метод подразумевает использование (в зависимости от требований к точности и скорости проводимых расчетов):
• композиции узлов сочленения трубопроводов для моделирования крановых площадок ЛЧМГ (при рассмотрении трубопроводной системы газотранспортного предприятия);
• подробных моделей межниточных перемычек (при рассмотрении небольших участков трубопроводной системы газотранспортного предприятия).
Композиция узлов сочленения трубопроводов представляет собой набор совокупностей трубопроводов, газодинамически связанных между собой. Газодинамически связанными считаются трубопроводы, между которыми может происходить переток газа в пределах моделируемой крановой площадки. В каждый момент времени композиция узлов сочленения трубопроводов должна отражать состояние кранов на межниточных перемычках. Для реализации динамического построения композиции узлов сочленения трубопроводов предложен алгоритм, в котором в качестве базовой модели узла сочленения трубопроводов используется модель, разработанная С. Н. Пряловым [12,14].
3. Разработан и научно обоснован новый метод автоматической настройки газодинамических симуляторов на реальные параметры трубопроводных сетей промышленных энергообъектов для повышения адекватности расчетных оценок основных параметров транспортирования газовых смесей при аварийных ситуациях. Он требует формализации технологической задачи подбора обобщенных параметров газопроводов (например, шероховатостей участков труб, местные гидравлические сопротивления и т. д.) в виде постановки и решения классической дискретной минимаксной задачи оптимизации. Решение минимаксной задачи проводится с использованием широко известного метода модифицированных функций Лагранжа.
4. С использованием разработанных методов получены новые результаты, расширяющие и углубляющие представления о:
• функционировании трубопроводных систем газотранспортных предприятий;
• причинах возникновения и механизмах протекания аварий на ЛЧМГ.
К таким результатам относятся: расчетные оценки параметров полного спектра безопасных режимов транспортирования природного газа, как через отдельные участки ЛЧМГ и технологические газопроводы (ТГ) КС, так и по сетям газотранспортных предприятий в целомрасчетные оценки параметров режимов безаварийного функционирования сетей КС газотранспортных предприятий при проведении ремонта и реконструкциианализ причин возникновения и сценариев развития аварий на объектах трубопроводного транспорта.
Предложенные новые методы и результаты их применения, направленные на повышение пожарной и промышленной безопасности объектов трубопроводных систем ТЭК, выносятся на защиту в виде научных положений диссертации, принадлежащих лично автору.
Практическая ценность работы. Выносимые на защиту методы реализованы в ряде производственных ГДС (построенных на базе программно-математического комплекса (ПМК) «CorNet»), предназначенных для повышения пожарной и промышленной безопасности и эффективности трубопроводных систем конкретных газотранспортных предприятий [15−17]. Эти ГДС активно используются для решения практических задач, возникающих при проектировании, эксплуатации и реконструкции трубопроводных сетей газотранспортных предприятий, как в России, так и за рубежом. Так, например, на базе ПМК «CorNet» при активном участии автора диссертации была разработана комплексная компьютерная аналитическая система «AMADEUS» для обеспечения требований безопасности, экономической эффективности и экологии в Международной газотранспортной компании «SPP» [27−29] (см. Приложение 1).
С помощью ПМК «CorNet» установлены причины и механизмы развития свыше десяти аварий (см., например, [15,30,31]). Он успешно применялся при решении задач ОАО «ГАЗПРОМ», Минатома РФ, Госгортехнадзора РФ, Международной газотранспортной компании «SPP», Сандийских национальных лабораторий (США) и др. (см., например, [27−35]).
Достоверность изложенных в диссертации основных научных положений обеспечивается [11,12]:
• научным обоснованием использования современных методов вычислительной механики газов и жидкостей, математических моделей и методов их численного анализа;
• результатами натурных и численных экспериментов;
• многолетней практикой успешного применения рассматриваемых в диссертации методов на реальных объектах ТЭК, как в России, так и за рубежом.
Благодарности. Автор диссертации выражает искреннюю признательность и глубокую благодарность своему научному руководителю, заместителю главного конструктора ФГУП «РФЯЦ-ВНИИЭФ» — начальнику отделения, доктору технических наук Селезневу Вадиму Евгеньевичу за постоянное внимание к работам автора диссертации, поддержку, консультации, научное руководство и практическую помощь при определении направлений научных исследований, положенных в основу диссертации.
Автор диссертации выражает признательность главному конструктору ФГУП «РФЯЦ-ВНИИЭФ», директору Научно-производственного комплекса ФГУП «РФЯЦ-ВНИИЭФ» Клишину Геннадию Семеновичу, заместителю начальника отделения ФГУП «РФЯЦ-ВНИИЭФ», кандидату технических наук Алешину Владимиру Васильевичу, начальнику научно-исследовательской группы, кандидату технических наук Прялову Сергею Николаевичу, начальнику научно-исследовательской лаборатории, кандидату технических наук Киселеву Владимиру Владимировичу и начальнику научно-исследовательской лаборатории Фотину Сергею Валентиновичу за научные консультации по теме диссертации и поддержку его научных исследований и разработок.
Автор благодарит за научные дискуссии по теме диссертации своих иностранных коллег из Математического института Словацкой Академии Наук, Братиславского государственного университета им. Комениуса и Международной газотранспортной компании «SPP».
Автор выражает признательность и благодарность ученым и ведущим специалистам Академии ГПС МЧС РФ за плодотворные дискуссии по теме диссертации.
Автор выражает искреннюю благодарность и признательность своим коллегам Мотлохову Владиславу Владимировичу и Комиссарову Алексею Сергеевичу за сотрудничество и поддержку.
СПИСОК ОСНОВНЫХ ИСПОЛЬЗУЕМЫХ СОКРАЩЕНИЙ.
ГДС — высокоточный газодинамический симулятор режимов транспорта природного газа через трубопроводные сети газотранспортных предприятий (или фрагменты трубопроводных сетей газотранспортных предприятий) — ГИС — географическая информационная системаГПА — газоперекачивающий агрегат на КС (КЦ) — ГРС — газораспределительная станцияГУ — граничные условия;
КС — компрессорная газоперекачивающая станция;
КЦ — компрессорный газоперекачивающий цех;
ЛЧМГлинейная часть магистральных газопроводов;
МГ — магистральный газопровод;
ПМК — программно-математический комплекс;
ТГ — технологические газопроводы на КС (КЦ);
ТЭК — топливно-энергетический комплекс;
ЦВТМ — Центр вычислительных технологий механики ООО «НПО.
ВНИИЭФ-ВОЛГОГАЗ" — ЦН — центробежный нагнетатель ГПА;
SCADA-система — (Supervisory Control And Data Acquisition) — система.
Символьные обозначения, применяемые в формулах, оговариваются особо в каждой Главе диссертации.
ВЫВОДЫ.
1) Разработана новая технология повышения пожарной и промышленной безопасности эксплуатации газопроводных сетей промышленных энергообъектов, основанная на методе обнаружения и локализации аварийных разрывов многониточных трубопроводов с использованием высокоточных газодинамических симуляторов транспортирования газов по трубопроводным системам.
2) Впервые предложен и верифицирован на тестовых и производственных задачах критерий наличия аварийного истечения газа из разрушившейся нитки.
3) Разработан и верифицирован на тестовых и производственных задачах метод моделирования крановых площадок (крановых узлов), оборудованных межниточными перемычками, для обеспечения оперативности обнаружения разрывов многониточных линейных частей магистральных газопроводов с использованием газодинамических симуляторов.
4) Впервые разработан, научно обоснован и верифицирован на тестовых производственных задачах метод автоматической настройки газодинамических симуляторов на реальные параметры трубопроводных сетей промышленных энергообъектов для повышения адекватности расчетных оценок основных параметров транспортирования газовых смесей в номинальных режимах и при аварийных ситуациях, включая пожары на магистральных газопроводах.
5) С помощью разработанных методов получены новые результаты:
• по расчетным оценкам параметров безопасных динамических режимов транспорта природного газа по сетям газотранспортных предприятий;
• по оценкам возможности реконструкции объектов трубопроводных систем;
• по сценариям развития аварий на магистральных газопроводах. Полученные результаты расширяют и углубляют представления о функционировании трубопроводных систем газотранспортных предприятий.
6) По результатам численных экспериментов определены рекомендации по разработке системы автоматического обнаружения разрывов для конкретных ЛЧМГ. При использовании датчиков давления с точностью ±-0,25атм расстояние между датчиками, размещаемыми на концах контролируемого трубопровода, должно составлять не менее 1 км. При этом время от возникновения разрыва до момента его идентификации и определения аварийной нитки составит ЗО+бОс.