Разработка программы для расчета гидратообразования в мг на программе Борланд дельфи 7. 0 Текст научной статьи по специальности «Химические технологии»

ДОБЫЧА

УДК 622.276 63

Разработка программы для расчета гидратообразования в мГ на программе борланд Дельфи 7.0

А.А. Паранук

аспирант1, см. инженер 2 кат.2 rambi.paranuk@gmail.com

1кафедра Машины и оборудование нефтяных и газовых промыслов, Кубанский государственный технологический университет (КубГТУ), Краснодар, Россия

2ПДС «Газпром трансгаз Краснодар» , Краснодар, Россия

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

материалы и методы

Использован язык программирования Дельфи 7.0.

Ключевые слова

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

Из-за низкого качества осушки газа на промыслах влага конденсируется в МГ, в результате чего снижается их пропускная способность и возникают условия для образования кристаллогидратов напоминающих внешним видом снег или лед. Кристаллогидраты были открыты английским химиком Дэви Х. 1810 году. По своей структуре газовые гидраты — соединения, которые образуется путем внедрения в пустоты кристаллических структур, восстановленных из молекул воды и газа. Обща формула газовых гидратов МпН20, где значение «п» изменяется от 5,75 до 17 в зависимости от состава газа и условий образования гидратов. Максимальное значение влагосодер-жания «при полном насыщении» зависят, от состава газа возрастая с увеличением содержания Н2Б, С02, и снижаясь с повышением содержания 1М2.

На диаграмме представлено условие образование гидратов. Гидраты образуются при наличии влаги в газе определенных давлениях и температурах. Гидраты образуются в областях расположенных влево от кривых 2, 5. При пересечении кривых 1, 2, 5 образуются критическая точка расположения гидратов точке С. В точке В образованной при пересечений кривых 2 и 3 , существует система гидрат — лед — вода — газа. При нагреве падение давления ниже критической температуры это приводит к разложению на воду и газ.

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

Способы борьбы с гидратами — понижение давление в системе ниже равновесного; повышения температуры газа выше равновесного; осушка газа для предупреждения конденсаций паров воды, ввод ингибиторов

(метанол, раствор ДЭГ, ТЭГ и др.)

Предупреждение образование гидратов подогревом газа заключается в том, что притом же давления в газопроводе температура газа поддерживается выше равновесной температуры образования гидратов. В условиях транспорта газа по МГ этот метод неприемлем, так как связан с большими затратами энергии. Как показываю расчеты, при больших объемах транспортируемого газа целесообразнее охлаждать его «с учетом затрат на более глубокую осушку газа», поскольку это позволяет заметно увеличить пропускную способность газопроводов, особенно газопроводов с большим числом компрессорных станции. Метод подогрева применяется на газораспределительных станция, где при больших перепадах давления вследствие дроссельного эффекта температура газа может значительно снижаться, в результате чего обмерзают редуцирующие клапаны, краны, и диафрагмы др.

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

Очевидно, что метод пригоден только для ликвидации гидратных пробок при положительных температурах. Иначе гидрат-ная пробка перейдет в ледяную фракцию. Поскольку минимальная температура газа в магистральных газопроводах близка к нулю, равновесное давления при этом для природного газа находиться в пределах от 1-1,5 МПа. Применение данного метода для

Рис. 1 — Диаграмма фазового состояния гидратов. 1 — упругость паров гидратообразователя; 2 — равновесные условия образования гидратов; 3 — равновесная граница между гидратом и льдом; 4 — понижение температуры замерзания воды; 5 — зависимость критической температуры разложения гидратов от давления.

1 ТЛН1

\ гга 283 \

V у

х

V ч.

т Ч

Ч - --

Т«268 1 -

«ли «дав о дм „.МП!

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

предупреждения гидратообразования в МГ оказывается не эффективным «оптимальное давление транспортируемого газа» от 5-7 МПа. Метод снижения давления применяется в аварийных случаях, для разложения гидратов в газопроводе в сочетаний с ингибиторами, так как в противном случай после повышения давления гидраты появляются вновь. Ингибитор, введенный в насыщенный водяными парами поток природного газа, частично поглощают водяные пары и переводят их вместо со свободной водой в раствор, который совсем не образует гидратов или образует их, но при более низких температурах. В качестве ингибиторов применяю метиловый спирт метанол, растворы этиленгликоля (ЭГ), диэти-ленгликоля (ДЭГ), триэтеленгликоля (ТЭГ), хлористого кальция, этилкарбитола (ЭК) и др. Удельный расход ингибитора для предупреждения процесса гидратообразования:

(д^ - а), к.

-■НО '«г.

(1)

с, - с,

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

Величина с2 определяется по графику (рис. 3). В зависимости от требуемого

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

йТ=Т-Т

Р г

(2)

где Тр — равновесная температура гидратообразования газа; Тг — температура газа в газопроводе.

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

При этом предпочтение следует отдать диэтиленгликолю, так как возможность его

регенерации и сравнительно небольшие потери большинстве случаев делают этот ингибитор наиболее экономичным.

При больших объемах транспортируемого газа его осушка является наиболее эффективным и экономичным способами предупреждения образования кристаллогидратов в МГ. При промысловой подготовке газа к дальнейшему транспорту его осушают сорбционными способами или охлаждением газового потока. В результате осушки газа точка росы паров воды должна быть снижена ниже минимальной температуры при транспортировке газа (влажность должна составлять не более 0,05-0,1 г/м3).

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

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

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

/

/

/

У

_^ У

£

Рис. з — График снижения температуры гидратообразования при вводе в газовый поток метанола

Рис. 4 — График образования гидратов природного газов в зависимости от температуры и давления

/ г /

I г

I 1 /

у / /

/ У

Рис. 5 — График для определения содержания метанола в жидкости

Рис. 6 — Зависимость содержания метанола Кж в паровой в жидкой фазах от давления и температуры газа Тж

-ям-' if-nnnh >u i :о:-гган-i in lurini ып^шл реки i^rfra^nK

jih ujh. F -.xi у ivyt- шары cjt i -auuipriu^iaMr|iiuu! >mviid 1»пнч гаакпЛ01||. r

il'h ЯЯР f+tf>H H<rtJTMPK*IMI

. >^flIIHT1 I ГЯ1 I --.».4.1 ■ I Tl Г Ы i'l!

ILtr^HlT МП |фИИ I HI НПНЦКНВН»

I

Вили литр »jrtbqmnji i ы

иг ЫШк »»»i кч II» raj»in eji «к мч1*(>пу|ф i+ii

Рис. 7

rii | ■ 11- ' i г.- n ■ . ^-.--И i nL ГЧ1Щ Ihll UmJ-HI.»*"' ■

Рис. 9

Для абсорбционной осушки природного газа в основном используется ДЭГ и ТЭГ (таб. 1).

При образований гидратных пробок в газопроводе их разрушают: путем ввода в трубопровод ингибиторов, снижение давление или подогревом газа. Место нахождения гидратной пробки определяю замером давлений на трассе «повышенный перепад давления, на каком либо участке свидетельствует об образовании гидратной пробки», при помощи радиолокационной антенны и передвижной радиолокационной станций, а также просвечиванием труб гамма излучением с помощью радиоизотопного прибора РИК-6М.

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

Методика расчета при гидратообразова-ний в газопроводе

(3)

где k — коэффициент теплопередачи в окружающую среду; D — наружный диаметр газопровода; плотность газа; Q — объемный расход газа, Ср — удельная теплоемкость газа.

Алгоритм программы для расчета формулы (з; 4) выглядит таким образом var

Formi: TFormi;

k, d, q, Cp, nu, T2, To, L, pi, p2, po, G:Real — здесь мы объявляем все переменные real — вещественный тип данных (диапазон значение 2.9* 10-39 до 1.7*10з8)[5]. begin-^ачало программы) программа просит ввести данные k: = StrToFloat(Editi.Text); // коэффициент теплопроводности d: = StrToFloat(Edit2.Text)/iooo; // диаметр трубопровода q: = StrToFloat(Edit3.Text); // Расход газа Cp: = StrToFloat(Edit4.Text); // теплоемкость газа To: = StrToFloat(Edit5.Text); // точка росы природного газа Ti: = StrToFloat(Edit6.Text); // Начальная температура газ T2: = StrToFloat(Edit7.Text); // Температура окружающей среды pi: = StrToFloat(Edit8.Text); // Давление в начале участка газопровода p2: = StrToFloat(Edit9.Text); // Давление в конце участка газопровода L: = StrToFloat(Editio.Text); // Длина участка газопровода nu: = StrToFloat(Editii.Text); // Джоуля-Томсона

то есть мы вводим нужное значение Edit.Text; StrToFloat-переводит изображение вещественное;

для формулы (i) — программа присваивает переменной f нашу расчетную формулу для f:=(24*k*pi*D)/(s*q*Cp);// расчет коэффициента после этого определяем x^^

Начало зоны конденсации x^. в газопроводе при Т=Тт.р. определяется по формулам (i;2) в программе дельфи это формула выглядит таким образом

-la

LXP) ~Pi) + <XTf J

(4)

где м — коэффициент Джоуля Томсона, р1 — начальное давление, р2 — конечное давление Тн — начальная температура газа, Т0 — температура окружающей среды, а — расчетный коэффициент, 1_ — длина участка газопровода, Тт.р. — температура точки росы газа.

^:=(пи*(р1-р2)+(Т1-Т2)*1_*0;

delta:=i/f*Ln(log/logi);// Начало зоны концесации Х.н

1 . [аЦТп-Т0) Л

= —hl —!-ü—llL+ \

a IviPi-P:)

(5)

где м — коэффициент Джоуля Томсона, р1 — начальное давление, р2 — конечное давление Тн — начальная температура газа, Т0 — температура окружающей среды, а — расчетный коэффициент, L — длина участка газопровода; [4]

По формуле (3) происходит расчет конца зоны конденсаций

Г2:=Ц*1:*(Т1-Т2)/(пи*р1-пи*р2);

l/f*Ln(f2+l); ае[1аТ:=(Т1-Т2); de[tа2:=(nu*pl-nu*p2); Определяем фактическую температуру в газопроводе, считая, что Т=Тт.р.

(6)

где Т, Т0 — температура соответственно газа в газопроводе и окружающей среды; Тн — начальная температура газа; а - расчетный коэффициент; x — расстояние от начала газопровода до рассматриваемой точки; ц — коэффициент Джоуля Томсона, L — длина участка газопровода, р1 — начальное давление, р2 — конечное давление.

Tni := (i-Exp(-f*delta))/f;

Tn2:=(T2 + (deltaT*Exp(-f*delta)))-(delta2/L)*Tni;

Tn3:=(i-Exp(-f*f3))/f;

Tn4 := (T2 + (deltaT*Exp(-f*f3)))-(delta2/L)*Tn3; Tn5:=Tn2+ Tn4;

Tn:= Round(Tn5/2); case Tn of 233: df:=(o.l45l); 234: df:=(o.i6i6); 235: df:=(o.i78o); 236: df:=(o.i985); 237: df:=(o.2i89); 238: df:=(o.243o); 239: df:=(o.267o); 24o: df:=(o.2953); 24i: df:=(o.3235); 242: df:=(o.3583); 243: df:=(o.393o); 244: df:=(o.4323); end;

путем подбора выбирается коэффициент

(7)

Параметр ДЭГ ТЭГ

Относительная молекулярная масса 106,12 150,17 Плотность (г/см3) при температуре, К:

293 1,118 1,125

288 1,119 1,127 Температура кипения (К) при давлении, МПа

0,1013 518 558

0,00133 401 435 Температура К,

начала разложения 437,5 479 начала замерзания

воспламенения на воздухе 350,5 173,9 Скрытая теплота парообразования (Дж/кг) при давлении

0,1 МПа 0,0357 0,0478 Вязкость динамическая (Па*с) при температуре 293К

Таб. 1

po:=2/3*(pi+(p2*p2)/(pi+p2)); // расчет среднего давления

Описывать весь алгоритм мы не будем в рамках это статьи, таким образом мы получаем

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

Коэффициент теплопроводности k = 1740 Диаметр трубопровода D = 1220 Расход газа Q = 12700000 Теплоемкость газа Cp = 2512 Точка росы природного газа Тт.р. = 255 Начальная температура газ Тнач. = 263 Температура окружающей среды Токр. = 238 Давление в начале участка газопровода Рнач. = 5,5

Давление в конце участка газопровода Рконц. = 4

Длина расчетного участка газопровода L = 120000

Коэффициент Джоуля-Томсона nu = 3,3 Плотность газа = 0,672

Расчетные данные:

Расчетный коэффициент a = 0,000007466

Начало зоны конденсации

Х.н = 40708метров = 40,708 км

Конец зоны конденсации

Х.к = 228941 метров = 228,940 км

Среднее давление газа в газопроводе

Рср = 4,79 МПа

Температура в начале зоны конденсации Т = 255,0 K

Температура в конце зоны конденсации Т = 238,0 K

Средняя температура в газопроводе на рассматриваемом участке Т = 247,0 K Влагосодержание газа 0,117 г/м3 Влагосодержание холодного газа 0,48 г/м3 Изменение влагосодержания газа 0,363 г/м3 Количество жидкой фазы, выделившейся из газа 4607 кг/сут

Температура образования гидрата (Тгидр) определяем на основании среднего давления (Рср) и плотности газа Тгидр = 285 К Определяем разность между температурой гидратообразования и начальной температурой газа (Тгидр-Тн) = 22 К Содержание метанола в жидкости Мж определяется на оснований (Тгидр-Тн)+273 К Содержание метанола Муж = 34% Концентрация метанола в газе Км = 0,108 г/м3 Количество метанола, насыщающего

жидкость вм.ж.=2373,4 кг/сут

Количество метанола, насыщающего газ

вм.г.=1365 кг/сут

Количество метанола вводимого

в газопровод вм=3739 кг/сут

Удельный расход метанола

рм=0,0002940 кг/м3

На рис. 1 мы видим рабочее поле программы где необходимо ввести следующие данные производительность газопровода Q (м3/сут.), диаметр газопровода (м), начальное давление, конечное давление (МПа), начальную температуру газа, температуру окружающей среды (К), длину расчетного участка (м), коэффициент Джоуля-Томсона (К/МПа), плотность газа (кг/м3), температуру точки росы (К), теплоемкость газа (Дж/К) и коэффициент теплопроводности (Вт/м2К).

После этого надо подвести курсор мыши надави кнопку «произвести расчет» и программа строить график гидратообра-зования на основе своего алгоритма и введенных данных как на рис. 2.

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

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

Итоги

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

Выводы

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

что требуется это ввести необходимые данные для расчета. Она сама выводит графики

необходимые для расчета в отдельном окне.

Список используемой литературы

1. Архангельский А.Я. Интегрированная среда разработки DELPHI. М.: Бином, i999. С. 34-54.

2. Буч Г. Объектно-ориентированное проектирование с примерами применения. Иркутск: Конкорд, i992. С. 54-76.

3. Кенту М. DELPHI-4 для профессионалов. СПб: Питер, i999. С. 63-78.

4. Орлик С.В. Секреты DELPHI на примерах. М.: Бином, i996. С. i3-24.

5. Стефен Моррис. Объектно-ориенти-

Abstract

Designed program allows make a payment, points of condensation, that is the beginning and the end of condensation determines the actual temperature at these points determines the average pressure of the gas in the pipeline, the moisture content of the gas, changes in moisture content, the amount of methanol injected into the pipeline, the specific consumption of methanol, all of these data written to a text file for later analysis.

рованное программирование. Серия «Enter». Ростов-на-Дону: Феникс, i997. С. 28-45.

6. Стефен Моррис. Объектно-ориентированное программирование. Серия «Enter». Ростов-на-Дону: Феникс, i997. С. 65-74.

7. Фаронов В.В. DELPHI 5. Москва: Но-лидж, i999. С. 56-65.

8. Фаронов В.В. DELPHI 6. Москва: Но-лидж, 2ooi. С. 76-85.

9. Гриценко А.И., Истомин В.А., Кульков А.Н., Сулейманов Р.С.

Сбор и промысловая подготовка газа на северных месторождениях России. М.: Недра, i999. С. 473.

10. Истомин В.А. Предупреждение образования газовых гидратов в системах

Materials and methods

Used programming language Delphi 7.0.

Results

This program allows you to simplify the calculation, make a detailed analysis of different modes of the pipeline.

Conclusions

It creates an array to determine the coefficients and reads from it the necessary given, based on the entered starting temperature and the estimated

сбора и промысловой подготовки газа.

М.: ИРЦ Газпром, 1996. №12. С. 23-31.

11. Катц Д.Л., Корнелл Д., Кобаяши Р., Поеттманн Ф.Х., Вери Дж. А, Еленбаас Дж., Уайнауг Ч.Ф. Руководство по добыче, транспорту и переработке природного газа. М.: Недра, 1965. С. 531.

12. Мустафин Ф.М., Коновалов Н.И., Гильметдинов Р.Ф. и др.

Машины и оборудование газонефтепроводов: Учеб. пособие М38 для вузов 2-е изд., перераб. и доп. Уфа. С. 60-65.

13. Макогон Ю.Ф. Газовые гидраты, предупреждение их образования и использование. М.: Недра, 1985. С. 232.

UDC 622.276

actual temperature of the gas in the pipeline. This program simplifies the calculation to the extent that is required is to enter the necessary data for the calculation. The program displays graphs required for the calculation in a separate window.

Keywords

hydrate formation, specific expense of methanol, moisture content, actual temperature, material type of data, range of values

References 6. Stefen Morris. Ob"ektno- of commercial gas]. Moscow: IRTs

1. Arkhangelsk A.Ya. Integrirovannaya orientirovannoe programmirovanie. Gazprom, 1996, issue 12, pp. 23-31.

sreda razrabotki DELPHI Seriya «Enter» [Object-oriented 11. Katz D.L., Cornell D., Kobayashi R.,

[Integrirovannaya sreda razrabotki programming. Enter series]. Rostov- Poettmann F.H., Verri A.J., Elenbaas J.,

DELPHI]. Moscow: Binom, 1999, pp. on-Don: Phoenix, 1997, pp. 65-74. Uaynaug C.F. Rukovodstvo po dobyche,

34-54. 7. Faronov V.V. DELPHI 5. Moscow: transportu i pererabotke prirodnogo

2. Buch G. Ob"ektno-orientirovannoe Nolidzh, 1999, pp. 56-65. gaza [Guidelines for the extraction,

proektirovanie s primerami 8. Faronov V.V. DELPHI 6. Moscow: transport and processing of natural

primeneniya [Object-oriented design Nolidzh, 200, pp. 76-85. gas]. Moscow: Nedra, 1965, p. 531

with application examples]. Irkutsk: 9. Gritsenko A.I., Istomin V.A., Kulkov 12. Makogon Y.F. Mashiny i oborudovanie

Concorde, 1992, pp. 54-76. A.N., Suleimanov R.S. Sbor i gazonefteprovodov: Ucheb. posobie

3. Kentu M. DELPHI-4 dlya professionalov promyslovaya podgotovka gaza na M38 dlya vuzov2-e izd., pererab. i

[DELPHI-4 M for professionals]. SPb: severnykh mestorozhdeniyakh Rossii dop. [Hydrates of Natural Gases].

Piter, 1999, pp. 63-78. [Collection and preparation of fishing Moscow: Nedra, 1974, p. 232

4. Orlik S.V. [DELPHI secrets on on the northern gas fields in Russia]. 13. Mustafin. F.M Konovalov N.I.

examples]. Moscow: Binom: 1996, pp. Moscow: Nedra, 1999, pp. 473 Gilmetdinov R.F. Gazovye gidraty,

13-24. 10. Istomin V.A. Preduprezhdenie preduprezhdenie ikh obrazovaniya

5. Stefen Morris. Ob"ektno- obrazovaniya gazovykh gidratov i ispol'zovanie [Machinery and

orientirovannoe programmirovanie. vsistemakh sbora i promyslovoy equipment, oil and gas pipelines:

Seriya «Enter» [Object-oriented podgotovkigaza [Preventing the Textbook. M38 manual for schools.

programming. Enter series]. Rostov- formation of gas hydrates in the etc. 2e ed., Rev. and add]. Ufa, pp.

on-Don: Phoenix, 1997, pp. 28-45. systems of collection and preparation 60-65.

ENGLISH OIL PRODUCTION

Development program for the calculation of hydrate

in oil and gas pipelines with using Bordland Delphi 7.0 software

Authors:

Arambiy A. Paranuk — graduate student1, engineer of 2nd category2; rambi.paranuk@gmail.com

•Machinery and equipment of oil and gas fields , Kuban State Technological University (KubGTU), Krasnodar, Russian Federation 2Gazprom Transgaz Krasnodar, Krasnodar, Russian Federation