Размещено на http://www.allbest.ru/

МОДЕЛИРОВАНИЕ АВАРИЙНЫХ РЕЖИМОВ УЧАСТКА ГАЗОПРОВОДА С УЧЕТОМ СВОЙСТВ РЕАЛЬНОГО ГАЗА

Н.Н. Елин, С.В. Крюков, Т.В. Корюкина

ФГБОУВО «Ивановский государственный политехнический университет», г. Иваново, Российская Федерация

Авторское резюме

газопровод прорыв потеря моделирование

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

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

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

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

Ключевые слова: аварийный режим, динамика истечения газа, массовый расход, процесс опорожнения, локализация порыва, давление газа, гидравлические потери.

Background: Ensuring reliable and uninterrupted supply of gas fuel is one of the most urgent tasks for energy generation systems. Design and construction of gas pipelines require a technical and environmental expert review to evaluate the possible after-effects of an emergency gas escape through pipeline breakages. Emptying the emergency shutdown section of the pipeline is a complex non-stationary process, in which all major parameters are linked. It is very important to know the dynamics of gas outflow to estimate the technological and ecological damage. This will allow reliable planning of the pipeline repair. Most of the known works examine the process as an isothermal outflow of a perfect gas, which may lead to considerable errors in prognosis. Hence, it is necessary to build a mathematical model of the process under study and its computer implementation.

Введение

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

Задачи аварийного истечения газа рассматривались в [2, 3], однако расчетные прогнозы часто не совпадали с практикой. На наш взгляд, основной причиной этого были допущения об изотермичности истечения и об идеальности газа, которые далеко не всегда соответствуют действительности, так как при этом игнорируется изменение в процессе истечения некоторых важных параметров газа - теплоемкости и показателя адиабаты. Кроме того, не учитываются изменения давления вблизи отверстия порыва вследствие движения к нему газа и температуры при расширении газа, обусловленные его неидеальностью.

В работах М.В. Лурье [4, 5] показано, что учет сверхсжимаемости газа при больших давлениях может привести к увеличению расчетного расхода утечки на 25 % и более. Однако автор задает коэффициент сверхсжимаемости в исходных данных и принимает его постоянным в течение всего процесса. На самом деле коэффициент сверхсжимаемости зависит от давления и температуры газа, которые могут сильно изменяться в процессе опорожнения аварийного участка, а также от компонентного состава газа. Учет компонентного состава особенно актуален для генераторов энергии, работающих на попутном газе нефтяных месторождений, например газотурбинных приводов.

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

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

1) на аварийно отключенном участке;

2) на работающем газопроводе.

Знание динамики этих процессов необходимо для принятия технически обоснованных решений о сроках проведения ремонтных работ и прогнозирования последствий аварийных истечений газа [6-8].

Методы исследования

Истечение через порыв на аварийно отключенном участке газопровода представляет собой нестационарный процесс с непрерывно изменяющимися расходными и теплофизическими параметрами. Для его расчета организуется цикл по времени = n•, где n = 0,1,2,… - номер шага. Величина шага подбирается исходя из опыта (чем меньше эквивалентный диаметр отверстия и чем больше начальная масса газа в трубе, тем больше ).

На каждом шаге рассчитывается массовый расход истечения m_n через порыв эквивалентным диаметром d [9]:

(1)

где k_q- коэффициент расхода, зависящий от формы отверстия (например, для отверстия с острой кромкой можно принять k_q? 0,62); у - отношение давления в окружающей среде к давлению в трубе в сечении локализации порыва p_a; _а - плотность газа в этом же сечении; k - показатель адиабаты.

Если величина у меньше критической, т.е.

(2)

то в расчете по формуле (1) принимается y = y_c.

Масса газа в участке составляет M_n = M_n-1 - m_n•, его плотность определяется как отношение этой массы к объему участка.

Температура газа рассчитывается по формуле

(3)

где Kt - линейный (отнесенный к единице длины трубопровода) коэффициент теплопередачи от газа в окружающую среду; Т₀ - температура окружающей среды. Теплоемкость с_р и показатель адиабаты рассчитываются по эмпирическим формулам [10].

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

Давление газа находится из уравнения

(4)

Так как величина z_n зависит от p_n, то следует применить метод последовательных приближений, положив в качестве первого приближения z_n = z_n-1.

Величину критического отношения давлений (2) следует пересчитывать на каждом временном шаге.

Для расчетов по уравнениям (1)-(4) необходимо знать теплофизические параметры газа. Для этого используем кубическое относительно коэффициента сверхсжимаемости уравнение состояния реального газа Пенга-Робинсона, коэффициенты которого рассчитываются по громоздким формулам, построенным в виде функций компонентного состава газа [11].

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

mp_n = (1- M_n/M₀), (5)

где M₀ - начальная масса газа в участке до начала истечения.

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

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

Давление перед порывом p_r на шаге n можно найти из закона сохранения массы. Так как массовый расход истечения равен сумме массовых расходов от концов участка к отверстию, то

m_1n + m_2n = m_n, (6)

где массовые расходы газа от концов участка к отверстию (l₁ + l₂ = L) составляют:

(7)

(8)

где - коэффициент гидравлического сопротивления участков.

Массовый расход истечения рассчитывается по формуле (1), в которую подставляются давление и плотность газа перед порывом p_r и _r.

На работающем газопроводе давление в сечении локализации порыва определяется по известной «газопроводной» формуле [1] при обычных допущениях о том, что течение газа изотермическое и происходит в квадратичной зоне сопротивления.

Массовый расход истечения через порыв эквивалентным диаметром dр рассчитывается по формуле (1).

Давление в конце газопровода длиной L с учетом порыва составляет

(9)

где p_b - давление в голове газопровода, в который подается массовый расход газа М; T_w - температура транспортируемого газа; L_a - расстояние от начала участка до места порыва.

Если задача заключается в определении места порыва при известном диаметре отверстия, то, предполагая, что истечение критическое (а это справедливо при давлении в трубе не менее 1,7-1,8 ата), из предыдущей формулы выразим

(10)

Результаты

На рис. 1-4 представлены некоторые результаты расчетов динамики процесса опорожнения аварийно отключенного участка газопровода диаметром 3256 мм, длиной 1000 м при начальных давлении и температуре газа в нем 55 ата и 303 К. Линейный коэффициент теплопередачи в окружающую среду составляет 3 Вт/(м•К). Расчеты выполнены для истечения в окружающую среду с атмосферным давлением и температурой 263 К через отверстия с эквивалентным диаметром 3, 5 и 7 мм.

Рис. 1. Динамика температуры газа: жирные линии - расчет по предлагаемой методике; тонкие - для изотермического истечения идеального газа при эквивалентном диаметре отверстия порыва: 1 - 3; 2 - 5; 3 - 7 мм; 4 - температура окружающей среды

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

Расчеты, выполненные при обычных допущениях об изотермическом течении идеального газа со средними теплофизическими свойствами (рис. 1-4, тонкие линии), показывают монотонное уменьшение температуры газа до температуры окружающей среды [2, 3].

Влияние адиабатного расширения увеличивается при увеличении эквивалентного диаметра отверстия порыва. Например, при диаметре 7 мм разница может достигать 19 К.

Так как величины, входящие в (1)-(4), зависят от температуры, то все параметры исследуемого процесса, рассчитанные по предлагаемой методике, будут отличаться от результатов, рассчитанных для изотермического истечения идеального газа.

На рис. 2 представлена динамика степени опорожнения аварийного участка газопровода.

Рис. 2. Динамика степени опорожнения аварийного участка: жирные линии - расчет по предлагаемой методике; тонкие - для изотермического истечения идеального газа, при эквивалентном диаметре отверстия порыва: 1 - 3; 2 - 5; 3 - 7 мм

Сравнение результатов показывает, что обычно используемые допущения [2, 3] могут привести к ошибкам до 15 %. Как и для динамики температуры, расхождения тем больше, чем больше сечение порыва.

На рис. 3, 4 представлены результаты расчетов динамики давления газа в отключенном участке и массового расхода истечения, выполненные по предлагаемой методике.

Рис. 3. Динамика давления газа в трубе, при эквивалентном диаметре отверстия порыва: 1 - 3; 2 - 5; 3 - 7 мм

Рис. 4. Динамика массового расхода утечки, при эквивалентном диаметре отверстия порыва: 1 - 3; 2 - 5; 3 - 7 мм

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

Тестовые расчеты показали, что размер порыва на процесс истечения оказывает заметное влияние при значениях d/D > 0,3 (рис. 5).

Рис. 5. Влияние размера порыва х = d/D на снижение давления перед отверстием y = p_r/p_b

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

На рис. 6-8 показаны некоторые результаты расчетов по предлагаемой методике, выполненные для газопровода диаметром 5308 мм, длиной 2850 м с давлением в голове 5,9 ата при средней температуре транспорта 303 К при подаче по нему газа относительной плотностью по воздуху 0,894 с массовым расходом 29,632 кг/с. При выполнении расчетов учитывалось отклонение свойств газа от идеального.

На рис. 6 показана зависимость давления в конце участка от эквивалентного диаметра порыва при заданной его локализации L_a = 1000 м.

Рис. 6. Зависимость давления в конце участка от диаметра порыва при заданной его локализации (1); давление в начале участка (2)

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

На рис. 7 показана зависимость давления в конце участка от его локализации при эквивалентном диаметре порыва 100 мм.

Рис. 7. Зависимость давления в конце участка от его локализации при заданном диаметре порыва

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

На рис. 8 показана зависимость локализации порыва (величины L_a) от диаметра порыва при давлении в конце участка p_e = 4,5 ата. Чем больше эквивалентный диаметр порыва, тем ближе он должен располагаться к концу участка для достижения заданной величины p_e. Для данного случая минимальная величина d равна 58 мм при расположении порыва в самом начале участка (L_a = 0). Максимальная величина d равна 170 мм, так как при этом весь газ уходит в порыв.

Рис. 8. Зависимость локализации порыва от диаметра порыва при заданном давлении в конце участка (1); общая длина участка(2)

На рис. 9 показана зависимость локализации порыва (величины L_a) от величины утечки, измеренной как разница между массовыми расходами на входе и выходе участка при разных давлениях в его конце p_e = 4,2; 4,5 и 4,8 ата.

Рис. 9. Зависимость локализации порыва от величины утечки при разных давлениях в конце участка: 1 -p_e = 4,2; 2 - 4,5; 3 - 4,8 ата; 4 - общая длина участка

Выводы

Проведенные исследования показали, что игнорирование реальных свойств газа при расчетах аварийных режимов газопроводов может привести к существенным ошибкам. Расход истечения на отдельных этапах процесса может быть занижен на 10-15 %, а температура завышена на 5-10 К.

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

Использование в предлагаемой математической модели сложного уравнения состояния Пенга-Робинсона позволяет, в отличие от известных методик [1-5], прогнозировать возможное изменение фазового состояния газа, т. е. частичную конденсацию его бензиновых фракций, обусловленную резким падением температуры (рис. 1), при которой протекание процесса значительно ускоряется [6]. Это особенно актуально для газотурбинных генераторов, работающих на попутном нефтяном газе.

Представленные математические модели использованы при разработке расчетных модулей аварийных режимов в составе программного комплекса OISPipe [6-8]. Они позволяют для аварийно отключенного участка газопровода прогнозировать динамику опорожнения, а для работающего участка - локализацию порыва по изменению давления и расхода газа в конце участка. Это необходимо для обоснованного планирования производства ремонтных работ и мероприятий по обеспечению безопасности оборудования и персонала.

Список литературы

1. Лурье М.В., Некляев А.В. Об одном опасном явлении, сопровождающем истечение газа из газопровода // Газовая промышленность. - 2008. - № 1. - С. 82-83.

2. Лурье М.В. Математическое моделирование процессов трубопроводного транспорта нефти, нефтепродуктов и газа. - М.: Изд. центр РГУ нефти и газа им. И.М. Губкина, 2012. - 456 с.

3. Куцова Е.В., Сердюков С.Г., Васильев Е.М. Математическое моделирование аварийных режимов магистральных газопроводов // Вестник Воронежского государственного технического университета. - 2011. - Т. 7, № 9. - С. 17-21.

4. Лурье М.В. Экспертиза утечек газа из резервуаров с высоким давлением // Территория нефтегаз. - 2014. - № 4. - С. 52-57.

5. Лурье М.В., Найденов Р.А. Уточненный расчет утечек газа через отверстия в стенках газопроводов высокого давления // Газовая промышленность. - 2014. - № 8(710). - С. 82-85.

6. Программный комплекс OISPipe для мониторинга и оптимизации систем сбора газа месторождений различных типов / Н.Н. Елин, А.В. Бардин, Д.В. Загинайко, А.П. Попов // Нефтяное хозяйство. - 2008. - № 5. - С. 65-69.

7. Программный комплекс OIS PIPE для математического моделирования сложных трубопроводных систем промыслового обустройства / Н.Н. Елин, Ю.В. Нассонов, О.В. Белоусов и др. // Нефтяное хозяйство. - 2002. - №12. - С. 91-93.

8. Моделирование и оптимизация сетей с периодически работающими источниками / Н.Н. Елин, А.П. Попов, Д.В. Загинайко, М.Г. Королев // Вестник ИГЭУ. - 2016. - Вып. 3. - С. 73-77.

9. Кириллин В.А., Сычев В.В., Шейндлин А.Е. Техническая термодинамика. - М.: Энергоатомиздат, 1983. - 512 с. 

10. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей: справочное пособие. - Л.: Химия, 1982. - 592 с.

11. Brill J.P., Mukherjee H. Multiphase flowin wells. Henry L. Doherty Fund of AIME Society of Petroleum Engineers Inc. Richardson. - Texas, 1999. - 157 p.

Размещено на Allbest.ru