Бесплатный автореферат и диссертация по наукам о земле на тему
Разработка методики расчета параметров течения сырого газа с малым содержанием жидкости в морских трубопроводах
ВАК РФ 25.00.18, Технология освоения морских месторождений полезных ископаемых

Автореферат диссертации по теме "Разработка методики расчета параметров течения сырого газа с малым содержанием жидкости в морских трубопроводах"

РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ НЕФТИ И ГАЗА ИМЕНИ И. М. ГУБКИНА

0034818Ь4

УДК 622.692.4 На правах рукописи

Бычкова Ольга Андреевна

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

Специальность 25.00.18 - Технология освоения морских месторождений полезных ископаемых (технические науки)

АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата технических наук

Москва-2009

003481854

Работа выполнена в Государственном образовательном учреждении высшего проффесионального образования «Российский государственный университет нефти и газа имени И. М. Губкина» и Обществе с ограниченной ответственностью "Газпром ВНИИГАЗ"

Научный руководитель:

доктор технических наук, профессор Мирзоев Дилижан Аллахверди оглы

Официальные оппоненты:

доктор технических наук, профессор Григорьев Леонид Иванович

кандидат технических наук Староконь Иван Викторович

Ведущая организация:

ОАО «Гипроспецгаз»

Защита состоится 25 ноября 2009г. в 1_5 час. 00 мин. на заседании диссертационного совета Д212.200.11 при Российском государственном университете нефти и газа им. И. М. Губкина по адресу: 119991, Москва, ГСП-1, Ленинский проспект, 65.

С диссертацией можно ознакомиться в библиотеке Российского государственного университета нефти и газа им. И. М. Губкина

Автореферат разослан «д^/ » 2009г.

Учёный секретарь диссертационного совета, доктор технических наук, доцент

И.Е. Литвин

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

Перспективы освоения месторождений и определения транспортных маршрутов связаны с территориями, расположенными в акватории шельфа.

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

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

Принимая во внимание вышесказанное и тот факт, что на шельфе РФ открыто множество газовых месторождений с малым содержанием конденсата и пластовой воды (Штокмановское, Северо-Каменномысское, Каменномысское-море, Семаковское и др.), точное проектирование гидравлических режимов в газопроводах от этих месторождений становится важной и актуальной задачей.

Цель работы

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

Основные задачи исследования

Поставленная цель достигается путем решения следующих основных задач:

■ Выбор критерия для оценки возможности накопления жидкости в трубопроводе сырого газа.

■ Разработка методики расчета параметров течения сырого газа с малым содержанием жидкости в морских газопроводах.

■ Апробация разработанного алгоритма на примере транспортной системы Штокмановского, Северо-Каменномысского месторождений и месторождения Каменномысское-море.

Научная новизна

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

Защищаемые положения

1. Алгоритм для оценки возможности накопления жидкости в морском трубопроводе сырого газа, основанный на вычислении выносной скорости газа и построении критериальной зависимости выносной скорости газа от давления, полученной из соответствующих зависимостей Стина-Уоллиса и Клапчука-Елина.

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

Практическая ценность результатов работы

В настоящее время ведутся активные работы по разработке проектов обустройства месторождений Баренцева и Карского морей, Обской и Тазовской губ, в частности Штокмановского газоконденсатного месторождения, Северо-Каменномысского газового месторождения и т.д. Результаты исследований данной работы использованы при разработке «Корректировки проекта разработки Штокмановского газоконденсатного месторождения», «Проекта разработки месторождения Каменномысское-море», «Проекта разработки Северо-Каменномысского месторождения», «Инвестиционного замысла освоения ресурсов Обской и Тазовской губ», «Корректировки «Программы освоения углеводородных ресурсов шельфа РФ на период до 2030 года».

Апробация работы

Основные положения диссертационной работы обсуждались на:

■ Международной конференции «Нефть и газ Арктического шельфа-2006» (2006 г., Мурманск);

■ Международной конференции «Освоение шельфа арктических морей России RAO/CIS Offshore 2007» (2007 г., Санкт-Петербург), где получен «Диплом» конференции RAO-07;

" Международной научно-практическая конференция

«Газотранспортные системы: настоящее и будущее (GTS - 2007)» (2007 г., Москва);

■ на Международной конференции «Освоение ресурсов нефти и газа российского шельфа: Арктика и Дальний Восток ROOGD-2008» (2008 г., Москва);

■ заседании секции «Освоение морских нефтегазовых месторождений» Ученого совета ООО «ВНИИГАЗ», Протокол № 3 от 9 апреля 2007 года;

■ совместных заседаниях секций «Техника и технология разработки морских месторождений» Научно-технического совета ОАО «Газпром» и секции «Освоение морских нефтегазовых месторождений» Ученого совета ООО «ВНИИГАЗ», от 4 октября 2006 года, 26 ноября 2008 года.

Объем и структура диссертационной работы

Диссертационная работа содержит введение, три главы, основные результаты с выводами, список использованной литературы из 57 наименований. Содержание изложено на 101 странице машинописного текста и включает 38 рисунков и 23 таблицы.

СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе проведен анализ развития теоретических исследований по движению газожидкостных смесей в трубах от простейших эмпирических зависимостей начала XX века (Тимерина, Костерина, Вильде, Мартинелли, Арманда, Телетова, Франкля и др.) и упрощенных стационарных моделей (Лоренца, Меликова, Гибсона, Лейбензона, Поэтманна-Карпентера, Баксендела-Томаса, Фэнчера-Брауна и др) до сложных механистических (Уоллиса, Хасана-Кабира, Ансари и др., Петаласа-Азиза, ВНИИГАЗ'а и др.) и современных динамических моделей (OLGA, TACITE, TUFF). Во второй части главы показано современное стояние применения технологии многофазного транспорта в мировой практике, а также приведены конкретные примеры по реализации проектов трубопроводных систем, транспортирующих сырой газ с малым содержанием жидкости (Brent, Yacheng, Frigg, Troll, Viking, Kvitebjom, Mensa и др.), которые могут послужить аналогами при проектировании подобных месторождений на российском шельфе, таких как Штокмановское, месторождения Обской и Тазовской губ, Карского моря.

Во второй главе описываются особенности течения сырого газа с малым содержанием жидкости (менее 56 см3 на 1 м3 газа или 10 баррелей на 1 млн. кубических стандартных футов газа):

" Результаты сравнения фактических и расчетных потерь давления для существующих трубопроводов с содержаниями жидкости от 5.6 до 45 см3/м3, показывают, что для потоков сырого газа с малым содержанием жидкой фазы разность между расчётными и фактическими значениями потерь давления

минимальна при использования гидравлических зависимостей АвА для сухого газа, а использование методик Беггса-Брилла и Даклера для расчёта двухфазного течения транспортируемой продукции приводит к большим погрешностям в оценке гидравлических потерь для широкого интервала значений производительности газопроводов.

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

Ж

Рисунок 1 - Структуры течения двухфазных потоков сырого газа с малым содержанием жидкости, характерных для протяжённых рельефных газопроводов: а - дисперсное течение; б - расслоенное течение с сегментной (наверху) и серповидной (внизу) поверхностью раздела

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

Формула Стина-Уоллиса выглядит следующим образом:

Кс= 0.246— (1)

Ра

Формула Клапчука-Елина имеет вид:

/ \0-25 ;----

Гс= 5.5

Рь-Ро

й-

где: - вязкость газовой фазы, сП;

а - поверхностное натяжение, Н/м;

рс- плотность газовой фазы, кг/м3;

Р1- плотность жидкой фазы, кг/м3.

Алгоритм оценки условий для накопления жидкости в полости трубопровода:

1. По формулам (1) и (2) рассчитывается скорость газа Ус для заданного диапазона давлений. При этом плотности жидкой фазы рплотности ро и вязкости ро газовой фазы, а также поверхностное натяжение а рассчитываются для заданного состава газа при заданной температуре сырого газа в программном продукте РУТБт (СаЬер).

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

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

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

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

* Раздельное сплошное движение газовой и жидкой компонент двухфазного потока подчиняется раздельным уравнениям сохранения импульса.

■ На границе раздела фаз действует сдвиговое напряжение, зависящее в основном от толщины слоя жидкости. Кроме того, оно зависит от плотности, скорости и вязкости газовой фазы. Сдвиговое напряжение вызывает замедление движения газового потока и ускорение слоя жидкости.

■ Геометрия слоя жидкости принята в виде симметричного кольцевого сегмента на нижней образующей трубы (Рисунок 1, 6). Выбор такой геометрии поверхности раздела фаз объясняется, в основном, соображениями упрощёния расчётных соотношений и не приводит к значительным погрешностям при определении величины истинного содержание жидкости при рассматриваемых в работе объёмах жидкости.

На основе предлагаемой модели можно:

■ рассчитывать толщину слоя жидкости в сечении трубы;

■ рассчитывать величину истинного содержания жидкости;

■ рассчитывать потери давления при транспортировке сырого газа;

■ определять границы существования раслоенно-волнового режима течения двухфазного потока.

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

| - + г,-5, + р1А1йпсс = 0, (3)

•■Л"* "

-4

с1х

- М ~ -гД. + рсАс эта = 0.

(4)

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

<РР1и1 + 0 - Ф)Рс,ио = соп^ • (5)

Исключая из уравнений (3) и (4) градиент давления, получаем:

\

11¥С

^с ^ V ^и

(6)

Касательные напряжения на стенках трубы вычисляются традиционным способом:

г -7 Р^к 'иг — А1 о

' па

Рсис

Т =Я ' ' 8

(7)

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

Г Г. V

Л?

(8)

Уравнение (6) с учётом уравнений (7), (8) и неравенства коэффициентов гидравлического сопротивления на поверхности раздела газ-жидкость и на границе газ-стенка трубы, имеет вид:

X'

А

уЛс

А

ЗА

^ у

-4У = 0

(9)

Для приведения переменных к безразмерной форме необходимо отнести их к: £) для длины; О1 для площади; ¡Д и ис для скорости жидкой и газовой фаз соответственно, где и - относительные скорости жидкой и газовой фаз при условии, что соответственно жидкость и газ заполняют весь периметр трубы 5, м/с. Для обозначения безразмерных величин используется символ (~).

В уравнении (6):

* "pv5Si- (10)

Y=(PL-PG)gS™a

(П)

\{dPldx)G\ '

SL=a- cos"' (2 hL -l), §G= cos"' (2A, -1),____= /l-(2Ä, -1)2 , (12)

= 0.25^ - cos"1 (2hL -1) + (2hL -1) Д- (lhL -1)2 ), (13)

2a = 0.25(cos_1 (2hL -1)-(lhL -1),/Г-(2^. -I)2), (14)

А А

аа = -ч~, *1 = Т> (15)

_ Л<

Хю ~ ] , (16)

Коэффициент гидравлического сопротивления на поверхности раздела фаз с учетом экспериментальных исследований Хамерсма и Харта по определению влияния толщины плёнки жидкости рассчитываем по формуле Эка (1973):

__ __0.25_

> .Г \Л2 • (18)

Log

15- + 0.619Я,

'Ч; LJJ

Коэффициент гидравлического сопротивления для газовой фазы lG определяется по формуле Колбрука-Уайта. Напомним, что Тайтель и Даклер решали уравнение (6) при условии xuf* 1 •

В (3) - (18) используются следующие обозначения:

Al и Ас - площадь поперечного сечения, занимаемая соответственно жидкостью и газом, м2;

Sü Sq и S, - части периметра трубы, занимаемые соответственно жидкостью, газом и границей раздела фаз, м;

Тц<1, Tjкс и т, - касательные напряжения соответственно между стенкой трубы и жидкостью, между стенкой трубы и газом и на поверхности раздела жидкости и газа, Н/м2;

Pi и Рс - плотность жидкости и газа соответственно, кг/м3;

а - угол наклона рассматриваемого участка газопровода, град;

Аь Ха и Я, - коэффициенты гидравлического сопротивления жидкости, газа и на поверхности раздела жидкости и газа соответственно;

UiHUa- скорости жидкости и газа соответственно, м/с;

Dl и Dg - гидравлические диаметры жидкости и газа соответственно, м;

ClhCc- параметрические коэффициенты;

п, т- показатели степени;

Хи Y- комплексы Локкарта-Мартинелли;

(dP/dX)c и {dP/dX)L - дифференциальные потери давления на трение, рассчитанные для газовой (жидкой) фазы при условии её перекачки в однофазном состоянии на участке трубопровода диаметра D с весовым расходом, равным весовому расходу газа (жидкости) в исходном двухфазном потоке.

При выводе (6) использовалось предположение о том, что m¿ « нетипичные значения показателей степеней пят для условий магистрального транспорта сырого газа составляют 0.25 и 0.1 соответственно.

Все величины, входящие в уравнение Тайтеля-Даклера (6), кроме комплексов X и У, определяются приведенной величиной высоты слоя жидкой фазы = h¡/D, где hL- геометрическая высота слоя жидкой фазы, м.

Значение истинного содержания жидкой фазы определяется по классической формуле:

'"íV (19>

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

1. Трасса газопровода делится на расчётные участки, характеризующиеся углами их наклона к горизонтали.

2. Термобарические и расходные параметры транспортировки сырого газа на каждом расчётном участке оцениваются с помощью гидравлических опций Беггса-Брилла или Мукерджи-Брилла, которые входят в состав известных программных комплексов PipePhase и PipeSim или программы OLGA 5.

3. На каждом расчётном участке по полученным на шаге 2 локальным значениям термобарических и расходных параметров рассчитываются значения X и К

4. С помощью уравнения (6) рассчитываются значения относительной высоты жидкой фазы, площади внутренней поверхности трубы и истинного содержания жидкой фазы на рассматриваемом расчётном участке. Решение осуществляется перебором величин hL; поиск искомого значения hL завершается при условии, что левая часть уравнения (6) становится меньше заранее заданной малой величины. Кроме того, рассчитывается истинное содержание жидкости и градиент давления на рассматриваемом участке с помощью уравнений (19) и (4) соответственно.

5. Кроме того, на каждом расчётном участке газопровода проверяется предположение о расслоено-волновой структуре течения сырого газа:

■ нижняя граница (переход в дисперсный режим) определяется с помощью значения отношения коэффициентов гидравлического сопротивления на поверхности раздела фаз и газовой фазы ххз с помощью уравнения (16). Если x¡c .< 1, то режим течения дисперсный, если же хю 1> то режим течения расслоенный;

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

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

Ргм<Ргт, (20)

где:

(21)

Гг =

1 ' »V

\Pl~PG ¡Dgcosa, С

1

'в'

(22)

Ла

0 = 1-^. (23)

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

Описанный алгоритм реализован с помощью алгоритмического языка Фортран.

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

Так как условия для С-КГМ и К-МГМ практически одинаковы, то в работе представлен расчет только для транспортной системы Северо-Каменномысского месторождения.

В начале главы подробно приведены исходные данные для расчета:

■ профили трасс газопроводов сырого газа;

" составы транспортируемых флюидов (два состава для ШГКМ с содержаниями жидкости 8.93 и 23.1 г/м3 и сеноманский газ С-КГМ с содержанием жидкости менее 1 г/м3);

■ производительность трубопроводов (23.1 млрд. м3/год для ШГКМ и 15.3 млрд. м3/год для С-КГМ);

■ характеристики труб (диаметры 42", шероховатость 25 мкм, толщины стенок 37 мм для ШГКМ и 14 мм для С-КГМ, условия окружающей среды);

■ граничные условия для расчета в программе (выходные давления 60 атм для ШГКМ и 28 атм для С-КГМ и входные температуры 40 С для ШГКМ и 20 С для С-КГМ).

Согласно описанным в главе 2 алгоритмам, сначала строятся кривые Стина-Уоллиса и Клапчука-Елина в заданном диапазоне давлений (Рисунок 2, Рисунок 3, Рисунок 4).

\ 1 1 ! ¡ i ! ~ — по Стииу-Уоллвсу — по Клялчуку-Елаяу ;

-— —

— -

— — — .....

20 30 40 SO (О 70 80 90 100 110 120 130 140 ISO 160 170 180 190 ZOO Давление, атм

Рисунок 2- Зависимость скорости газа от давления для сырого газа для сырого газа ШГКМ с содержанием жидкости 8.93 г/м3

Рисунок 3- Зависимость скорости газа от давления для сырого газа для сырого газа ШГКМ с содержанием жидкости 23.1 г/м3

Рисунок 4- Зависимость скорости газа от давления для сырого газа С-КГМ с содержанием жидкости менее 1 г/м3

Строится критериальная кривая Клапчука-Стина путем соединения верхних половин полученных кривых Стина-Уоллиса и Клапчука-Елина, если они пересекаются (Рисунок 5, Рисунок 6). В противном случае (Рисунок 7) выбираем кривую с большими значениями скоростей газа с целью получения более надежных результатов. Таким образом, мы считаем, что в области скоростей газа, лежащих выше данных кривых, жидкость не будет накапливаться в полости газопровода.

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

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

Далее проводится моделирование процесса течения сырого газа по трубопроводу. В результате компьютерного моделирования процесса транспортировки сырого газа получаем зависимости скорости газа от давления по длине трубопроводов для различной загрузки трубопровода от 100% до 50% с интервалом в 10% и сравниваем их с соответствующими критериальными кривыми (Рисунок 5, Рисунок 6, Рисунок 7).

12 ------1 ,— I "I----■_. |__

10

i « 8

w et u

§ 4 ----

U ____________

2------------------

0 -----------

60 70 80 90 100 110 120 130 140 150 160 170 Давление, атм

Рисунок 5- Сравнение критериальной кривой со скоростями газа для трубопровода сырого газа ШГКМ с содержанием жидкости 8.93 г/м3

--- -100% загрузка -----90% загрузка

--------- •70% загрузка - - 60% загрузка — 50% загрузка -критериальная кривая

••О-

1

------- ------------

60 70 80 90 100 110 120 130

Давление, атм

140 150 160 170

Рисунок 6- Сравнение критериальной кривой со скоростями газа для трубопровода сырого газа ШГКМ с содержанием жидкости 23.1 г/м3

Рисунок 7- Сравнение критериальной кривой со скоростями газа для трубопровода сырого газа С-КГМ с содержанием жидкости менее 1 г/м3

Результаты сравнения показывают, что для обеспечения уноса жидкости из полости трубы трубопроводы необходимо загружать минимум на 60% от проектного уровня загрузки, то есть производительность должна быть не менее 40.6 млн. м /сут. для трубопровода от Штокмановского месторождения и 27.9

млн. м3/сут. для трубопровода от Северо-Каменномысского месторождения, а давление на входе в трубопровод нужно держать на уровне 110 атм и 34 атм соответственно. При меньших значениях расхода сырого газа режимы эксплуатации газопровода отличаются существенной нестабильностью, что будет приводить к технологическим осложнениям в их работе и отрицательно сказываться на качестве подготовки и обработки газа.

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

Таблица 1 - Результаты расчета толщины слоя жидкости и истинного содержания жидкости для различных участков трассы трубопровода от ШГКМ

Конденсатный фактор =8.93 Конденсатный фактор =23.1

Расстояние от г/м3 г/м3

начала Толщина Истинное Толщина Истинное

газопровода, слоя содержание слоя содержание

км жидкости, жидкости, доли жидкости, жидкости, доли

мм ед. мм ед.

50 4.1 0.0004 0 0

100 7.5 0.0011 0 0

150 10.8 0.0019 0 0

200 12.0 0.0022 0 0

250 14.5 0.0029 21.2 0.0063

300 15.1 0.0031 27.3 0.0075

350 13.2 0.0025 26.7 0.0073

400 16.3 0.0035 35.3 0.0111

450 19.8 0.0047 37 0.0119

500 19.9 0.0047 37.1 0.0119

550 20.5 0.0049 36.5 0.0116

573 16.7 0.0036 28.5 0.0081

Таблица 2 - Результаты расчета толщины слоя жидкости и истинного содержания жидкости для различных участков трассы трубопровода от С-КГМ

Расстояние от начала Толщина слоя Истинное содержание

газопровода, км жидкости, мм жидкости, доли ед.

10 4.4149 0.00047

20 4.7473 0.00052

30 5.0070 0.00057

40 5.1836 0.00059

48 5.2978 0.00061

50 5.2979 0.00061

Таким образом, мы видим, что в обоих рассматриваемых трубопроводах от ШГКМ наблюдается накопление жидкой фазы (толщина слоя жидкости для лёгкого состава сырого газа ШГКМ увеличивается от 0 в начале трубопровода до 17 мм в конце, для тяжёлого состава газа ШГКМ - от 0 до 29 мм соответственно). Что касается транспортной системы С-КГМ, то здесь накопление жидкости выражено слабо (толщина слоя жидкости колеблется в пределах от 4 до 5 мм).

Далее определяем нижнюю границу существования расслоенного режима (то есть границу перехода его в дисперсный режим) с помощью критерия Хю по формуле (16) и верхнюю границу существования расслоенного режима (то есть границу перехода его в кольцевой или пробковый режим) с помощью критерия Кельвина-Гельмгольца по формуле (20) (Таблица 3, Таблица 4).

Таблица 3 - Значения критериев в рассматриваемых сечениях трубопровода от ШГКМ

Расстояние от начала газопровода, км Конденсатный фактор 8.93 г/м3 Конденсатный фактор 23.1 г/м1

FrM Frrtv JCiG FrM Frrev XiG

50 0.749 2.464 0.74 0 0 0

100 0.736 2.119 0.9 0 0 0

150 0.723 1.992 1.05 0 0 0

200 0.739 1.871 1.1 0 0 0

250 0.751 1.779 1.19 0.729 1.604 1.41

300 0.765 1.761 1.2 0.724 1.495 1.54

350 0.785 1.827 1.14 0.762 1.504 1.56

400 0.812 1.723 1.25 0.789 1.385 1.81

450 0.852 1.636 1.36 0.829 1.365 1.85

500 0.913 1.634 1.37 0.889 1.364 1.86

550 1.011 1.619 1.39 0.988 1.371 1.84

573 1.075 1.713 1.26 1.053 1.475 1.62

Таблица 4- Значения критериев в рассматриваемых сечениях трубопровода от С-КГМ

Расстояние от

начала FrM Frrev XiG

газопровода, км

10 0.470 2.444 0.305

20 0.473 2.399 0.306

30 0.478 2.368 0.307

40 0.483 2.346 0.307

48 0.549 2.340 0.307

50 0.507 2.334 0.308

Принятые критерии существования расслоенного режима показывают, что практически вдоль всей трассы газопровода ШГКМ (где выпадает жидкая фаза,

кроме начального участка) наблюдается расслоенная структура течения двухфазного потока с высотой слоя жидкой фазы в интервале от 0 до 17 мм для лёгкого состава сырого газа ШГКМ, от 0 до 29 мм для тяжёлого состава газа ШГКМ. Что касается трубопровода от С-КГМ, то критерий хю определяет существующий режим как дисперсный с очень малым слоем жидкой фазы (от О до 5 мм).

Здесь следует отметить, что предложенный нами критерий перехода расслоенного режима течения в дисперсный показывает более точные результаты по сравнению с существующими представлениями о гидравлике сырого газа с малым содержанием жидкости, согласно которым течение на участке газопровода характеризуется чисто дисперсной структурой, если в результате расчётов значение истинного содержания жидкости <р составляет 0.005-^0.006 и меньше. Данные значения <р основаны на обобщении экспериментальных наблюдений, и при установлении границы смены режимов не принималась во внимание незначительная по толщине пленка жидкости. Таким образом, граница перехода дисперсного режима в расслоенный смещается в область меньших истинных содержаний жидкости по сравнению с экспериментальными наблюдениями и допущениями, принятыми при этом.

Таблица 5- Результаты расчета толщины слоя жидкости для различных участков трассы трубопровода сырого газа от ШГКМ с содержанием жидкости 8.93 г/м3

Без учета коэффициента С учетом коэффициента

гидравлического гидравлического

Расстояние сопротивления на границе сопротивления на границе

от начала раздела фаз раздела фаз

газопровода, Толщина Истинное Толщина Истинное

км слоя содержание слоя содержание

жидкости, жидкости, жидкости, жидкости,

мм доли ед. мм доли ед.

50 3.75 0.00038 4.13 0.00044

100 7.2 0.001 7.45 0.0011

150 11.03 0.002 10.8 0.0019

200 12.4 0.0023 12.0 0.0022

250 15.4 0.0032 14.5 0.0029

300 16.05 0.0034 15.0 0.0031

350 13.7 0.0027 13.4 0.0025

400 17.6 0.0039 16.3 0.0035

450 21.9 0.0054 19.7 0.0047

500 22.03 0.0055 19.9 0.0047

550 22.9 0.0058 20.5 0.005

573 18.1 0.004 16.7 0.0036

Рассчитанный по формуле (18) коэффициент гидравлического сопротивления на межфазной границе Я, для типичных условий транспортировки

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

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

Факт смещения границ существования расслоено-волнового течения сырого газа по сравнению с рассчитанными по методике Тайтела-Даклера подтверждается результатами экспериментальных исследований Вуда.

Основные результаты и выводы

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

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

3. Разработан алгоритм для оценки возможности накопления жидкости в морском трубопроводе сырого газа, основанный на вычислении выносной скорости газа и построении критериальной кривой, полученной из кривых Стина-Уоллиса и Клапчука-Елина

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

5. Данная методика позволила провести расчет названных показателей для газопроводов сырого газа от Штокмановского, Северо-Каменномысского месторождения и месторождения Каменномысское-море.

Перечень опубликованных работ

1. Бычкова, O.A. Некоторые проблемы транспорта газожидкостной смеси Штокман-берег [Электронный ресурс] / O.A. Бычкова // Науч. сб. тр. IV международной конференции «Нефть и газ Арктического шельфа», Мурманск, 15-17 ноября 2006. - 1 эл. опт. диск (CD-ROM).

2. Бычкова, O.A. Оценка рабочих параметров безопасной транспортировки сырого газа по морскому трубопроводу / O.A. Бычкова, Д.А. Мирзоев, В.А. Сулейманов // Технологии ТЭК. - 2007. - № 1. - С. 54-58.

3. Бычкова, О. А. Исследование гидравлического режима транспортировки углеводородного сырья в многофазном состоянии (на примере Штокмановского газоконденсатного месторождения) / O.A. Бычкова //Сб. мат. 7-ой всерос. науч.-тех. конф. «Актуальные проблемы состояния и развитая нефтегазового комплекса России», Москва, 29-30 января 2007. - М. : РГУ нефти и газа им. И.М. Губкина, 2007. - С. 63-64.

4. Бычкова, O.A. Технологические требования к режимам транспортировки сырого газа по протяжённым морским трубопроводам [Электронный ресурс] / O.A. Бычкова, Д.А. Мирзоев, В.А. Сулейманов // Сб. тр. межд. конф. RAO/CIS Offshore, Санкт-Петербург, 11-13 сентября 2007. - 1 эл. опт. диск (CD-ROM).

5. Бычкова, O.A. Технологические требования к режимам транспортировки сырого газа по протяжённым морским трубопроводам / O.A. Бычкова, В.А. Сулейманов // Сб. тез. 2-ой межд. науч.-пр. конф. «Газотранспортные системы: настоящее и будущее (GTS-2007)», Развилка, 6-7 декабря 2007. - М. : ВНИИГАЗ, 2007. - С. 55.

6. Бычкова, O.A. Транспортировка сырого газа по морским трубопроводам: результаты исследования режимов течения / O.A. Бычкова, В.А. Сулейманов // Oil&Gas Journal Russia. - 2009. - № 6. - С. 36-42.

7. Бычкова О. А. Исследование режимов течения сырого газа с малым содержанием жидкости в протяженных морских трубопроводах / O.A. Бычкова, В.А. Сулейманов // Трубопроводный транспорт: теория и практика. — 2009. - № 2

(14).-С. 12-15.

Соискатель

Бычкова O.A.

Подписано к печати /¿Г /£ 09 Формат 60x90/16 Бумага офсетная Усл. п. л.

Тираж /00 экз. Заказ

Издательский центр РГУ нефти и газа имени И.М. Губкина 119991, Москва, Ленинский проспект, 65 Тел. (499) 233-93-49

Содержание диссертации, кандидата технических наук, Бычкова, Ольга Андреевна

Специальность — 25.00.18 ~ «Технология освоения морских месторождений полезных ископаемых»

Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель: д.т.н., профессор Д. А. Мирзоев

Москва, 2009 г.

СОДЕРЖАНИЕ

ВВЕДЕНИЕ.б

ГЛАВА 1 ИСТОРИЯ РАЗВИТИЯ ТРАНСПОРТИРОВКИ МНОГОФАЗНОЙ

ПРОДУКЦИИ ПО МОРСКИМ ТРУБОПРОВОДАМ.

1.1 История развития теоретических представлений о течении газожидкостной смеси по трубопроводам.

1.1.1 Эмпирические соотношения.

1.1.2 Упрощенные стационарные модели.

1.1.3 Механистические модели.

1.1.4 Динамические модели (программные продукты).

1.2 Примеры морских многофазных трубопроводных систем с малым содержанием жидкости.

ГЛАВА 2 МЕТОДИКА РАСЧЁТА ПАРАМЕТРОВ ТЕЧЕНИЯ СЫРОГО ГАЗА С

МАЛЫМ СОДЕРЖАНИЕМ ЖИДКОСТИ В МОРСКИХ ГАЗОПРОВОДАХ.

2.1 Сравнение фактических и расчетных потерь давления в газопроводах сырого газа с малым содержанием жидкости.

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

2.3 Разработка методики расчёта параметров течения сырого газа с малым содержанием жидкости в морских газопроводах.

2.3.1 Определение условий для накопления жидкости в трубопроводе.46%

2.3.2 Механистическая модель течения сырого газа по трубопроводу.

2.3.3 Алгоритм расчета.

ГЛАВА 3 РАСЧЁТ ПАРАМЕТРОВ ТЕЧЕНИЯ СЫРОГО ГАЗА С МАЛЫМ СОДЕРЖАНИЕМ ЖИДКОСТИ НА ПРИМЕРЕ ТРАНСПОРТНОЙ СИСТЕМЫ ОТ ШТОКМАНОВСКОГО ГАЗОКОНДЕНСАТНОГО МЕСТОРОЖДЕНИЯ.

3.1 Исходные данные для расчета.

3.2 Результаты расчета.

3.2.1 Проверка условия накопления жидкости в трубопроводе.76 ■

3.2.2 Расчёт параметров двухфазного потока.

3.3 Выводы.

ВЫВОДЫ

Введение Диссертация по наукам о земле, на тему "Разработка методики расчета параметров течения сырого газа с малым содержанием жидкости в морских трубопроводах"

Перспективы освоения месторождений и определения транспортных маршрутов связаны с территориями, расположенными в акватории шельфа.

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

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

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

За последние 70 лет было разработано множество теорий и экспериментальных зависимостей- для предсказания падения давления и истинного содержания жидкости при течении двухфазных потоков. Тем не менее, ни одно из них не является универсальным и не позволяет достаточно точно рассчитывать термобарические и расходные параметры потоков для всех структур течения, характерных для двухфазных потоков углеводородных систем. Что касается вопросов транспортировки сырого газа с малым 6 содержанием жидкости, то они относятся к наименее изученным в истории исследований многофазных течений, и поэтому информация по этим системам весьма ограничена.

На шельфе РФ открыто множество газовых месторождений с малым содержанием конденсата и пластовой воды (Штокмановское, Северо-Каменномысское, Каменномысское-море, Семаковское, Киринское и др.). Проектирование газопроводов от таких месторождений является важной и актуальной задачей. Необходимо использовать такую методику расчётов двухфазных потоков, которая позволила бы прогнозировать режимы течения двухфазных потоков и обеспечивать приемлемую с точки зрения инженерных приложений оценку значений потерь давления в трубопроводе и истинного содержания жидкости в потоке.

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

Поставленная цель достигается путем решения следующих основных задач:

Выбор критерия для оценки возможности накопления жидкости в трубопроводе сырого газа.

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

Апробация разработанного алгоритма на примере транспортной системы Штокмановского, Северо-Каменномысского месторождений и месторождения Каменномысское-море.

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

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

Защищаемые положения

1. Алгоритм для оценки возможности накопления жидкости в морском трубопроводе сырого газа, основанный на вычислении выносной скорости газа и построении критериальной зависимости выносной скорости газа от давления, полученной из соответствующих зависимостей Стина-Уоллиса и Клапчука-Елина.

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

Практическая ценность результатов работы

В настоящее время ведутся активные работы по разработке проектов обустройства месторождений Баренцева и Карского морей, Обской и Тазовской губ, в частности Штокмановского газоконденсатного месторождения, Северо-Каменномысского газового месторождения и т.д. Результаты исследований данной работы использованы при разработке «Корректировки проекта разработки Штокмановского газоконденсатного месторождения», «Проекта разработки месторождения Каменномысское-море», «Проекта разработки Северо-Каменномысского месторождения», «Инвестиционного замысла освоения ресурсов Обской и Тазовской губ», «Корректировки «Программы освоения углеводородных ресурсов шельфа РФ на период до 2030 года».

Апробация работы

Основные положения диссертационной работы обсуждались на:

Международной конференции «Нефть и газ Арктического шельфа-2006» (2006 г., Мурманск);

Международной конференции «Освоение шельфа арктических морей России RAO/CIS Offshore 2007» (2007 г., Санкт-Петербург), где получен «Диплом» конференции RAO-07;

Международной научно-практической конференции «Газотранспортные системы: настоящее и будущее (GTS - 2007)» (2007 г., Москва);

Международной конференции «Освоение ресурсов нефти и газа российского шельфа: Арктика и Дальний Восток ROOGD-2008» (2008 г., Москва); заседании секции «Освоение морских нефтегазовых месторождений» Ученого совета ООО «ВНИИГАЗ», Протокол № 3 от 9 апреля 2007 года; совместных заседаниях секций «Техника и технология разработки морских месторождений» Научно-технического совета ОАО «Газпром» и секции «Освоение морских нефтегазовых месторождений» Ученого совета ООО «ВНИИГАЗ», от 4 октября 2006 года, 26 ноября 2008 года.

Объем и структура диссертационной работы. Диссертационная работа содержит введение, три главы, основные результаты с выводами, список использованной литературы из 62 наименований. Содержание изложено на 102 страницах машинописного текста и включает 38 рисунков и 23 таблицы.

Заключение Диссертация по теме "Технология освоения морских месторождений полезных ископаемых", Бычкова, Ольга Андреевна

выводы

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

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

3. Разработан алгоритм для оценки возможности накопления жидкости в морском трубопроводе сырого газа, основанный на вычислении выносной скорости газа и построении критериальной кривой, полученной из кривых Стина-Уоллиса и Клапчука-Елина

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

5. Данная методика позволила провести расчет названных показателей для газопроводов сырого газа от Штокмановского, Северо-Каменномысского месторождения и месторождения Каменномысское-море.

Библиография Диссертация по наукам о земле, кандидата технических наук, Бычкова, Ольга Андреевна, Москва

1. Арманд, А. Сопротивление двухфазного потока в горизонтальных трубах / А. Арманд // Известия. - 1946. -№ 15. - С. 6.

2. Брилл, Дж.П. Многофазный поток в скважинах / Дж.П. Брил, X. Мукерджи. — Москва-Ижевск : Институт компьютерных исследований, 2006. — 384 с.

3. Еловацкий, И.П. Нефтяная и газовая промышленность капиталистических и развивающихся стран / И.П. Еловацкий. — М. : Просвещение, 1977. — 104 с.

4. Ершов, Т.Б. Математическое моделирование нестационарных газожидкостных потоков в системе пласт-скважина: дис. . канд. техн. наук : защищена 31.10.2006 : утв. 20.01.2007 / Т.Б. Ершов. М. : РГУ нефти и газа им. И.М. Губкина, 2006. - 150 с.

5. Гутман, Э.М. Защита газопроводов нефтяных промыслов от сероводородной коррозии / Э.М. Гутман, М.Д. Гетманский, О.В. Клапчук, JI.E. Кригман. М. : Недра, 1988.-200 с.

6. Костерин, С.И. Исследование структур потока двухфазной среды в горизонтальных трубах / С.И. Костерин // Изв. ОТН АН СССР. 1943. - № 7. - С. 37-43.

7. Мамаев, В.А., Гидродинамика газожидкостных смесей в трубах / В.А Мамаев, Г.Э Одишария, Н.И. Семенов, А.А. Точигин. М. : Недра, 1969. - 208 с.

8. Медведский, Р.И. Метод расчета движения газонефтяной смеси в вертикальных трубах / Р.И Медведский, Ф.Г. Аржанов, В.А. Попов // Труды Гипротюменненфтегаза. 1974. - Вып. 41. - С. 26-31.

9. Мирзоев, Д.А. Оценка рабочих параметров безопасной транспортировки сырого газа по морскому трубопроводу / Д.А. Мирзоев, В.А. Сулейманов, О.А. Бычкова // Технологии ТЭК. 2007. - № 1. - С. 54-58.

10. Сулейманов, В.А. Гидравлика в газопроводах с внутренним гладкостным покрытием труб / В.А. Сулейманов // Наука и техника в газовой промышленности. 2006. - № 1.-С. 65-71.

11. Сулейманов, В.А., Транспортировка сырого газа по морским трубопроводам: результаты исследования режимов течения / В.А. Сулейманов, О.А. Бычкова // Oil&Gas Journal Russia. 2009. - № 5. - С. 36-42.

12. Телетов, С.Г. Уравнения гидродинамики двухфазных жидкостей / С.Г. Телетов // ДАН СССР. 1945. - Т. 50. - С. 99-102.

13. Транспортировка многофазных смесей по трубопроводам: пер. с англ. // Technoscoop. 1997. -№ 14 (1). - 84 с.

14. Уоллис, Г. Одномерные двухфазные течения / Г. Уоллис ; пер. с англ. B.C. Данилина, Ю.А. Зайгарника. М.: Мир, 1972. - 440 с.

15. Франкль, Ф. И. Уравнения энергии для движения жидкостей с взвешенными наносами / Ф.И. Франкль // ДАН СССР. 1955. - Т. 102. - № 5. - С. 903-906. ,

16. Чисхолм, Д. Двухфазные течения в трубопроводах и теплообменниках / Д. Чисхолм ; пер. с англ. Б.Л Кривошеина. М. : Недра, 1986. - 204 с.

17. Agrawal,. S. An analysis of horizontal stratified two-phase flow in pipes / S. Agrawal, G.A. Gregory, G.W. Govier // Canadian Journal of Chemistry Engineering. 1973. - № 51.-PP. 280-286.

18. Ansari, A. M. A comprehensive mechanistic model for upward two-phase flow in wallboards / A.M. Ansari // SPE Production Engineering 1994.-№ 5. - PP. 143-147.

19. Baxendell, P.B. The calculation of pressure gradients in high-rate flowing wells / P.B. Baxendell, R. Thomas // Journal of Petroleum Technology. 1961. - № 10. - PP. 10231028.

20. Beggs, H.D. A study of two-phase flow in inclined pipes / H.D. Beggs, J.P. Brill // Journal of Petroleum Technology. 1973. - № 4. - PP. 607-617.

21. Bendiksen, К. The dynamic two-fluid model OLGA: theory and application / K. Bendiksen, D. Malnes, R. Мое, S. Nuland // SPE Production Engineering. 1991. - № 6.-PP. 171-180.

22. Benzoni-Gavage, S. Multi-Dimensional Hyperbolic Partial Differential Equations: First-Order Systems and Applications / S. Benzoni-Gavage, D. Serre. Oxford University Press, 2007. - 536 pp.

23. Caetano, E. F. Upward vertical two-phase flow through an annulus / E.F. Caetano, O. Shoham, J.P. Brill // Journal of Energy Resource Technology. 1992. - V. 114. — № 1. — PP. 1-31.

24. Fancher, G.H. Prediction of pressure gradients for multiphase flow in tubing / G.H. Fancher., K.E. Brown // SPE Journal. 1963. - V. 3. - № 1. - PP. 59-69.

25. Fernandes, R.S. Hydrodynamic model for gas-liquid slug flow in vertical tubes / R.S. Fernandes., T. Semait, A.E. Dukler // AIChE Journal. 1983. - V. 29. - № 6. - PP. 981989.

26. Gould, T.L. Design of offshore gas pipelines accounting for two-phase flow / T.L. Gould, E.L. Ramsey // Journal of Petroleum Technology. 1975. - V. 27. - № 1. - PP. 366-374.

27. Hasan, A.R. Performance of a two-phase gas-liquid flow model in vertical wells / A.R Hasan., C.S. Kabir // Journal of Petroleum Science and Engineering. 1990. - № 4. -PP. 273-289.

28. Hewitt, G.F. Annular two-phase flow / G.F. Hewitt, N.S. Hall-Taylor. Oxford : Pergamon Press, 1970. - 310 pp.

29. Hope, P.M. A new approach to the design of wet-gas pipelines / P.M. Hope, R.G Nelson // Proceedings of ASME Conference (Houston, USA, September, 1977). Houston, 1977.-PP. 136-142.

30. Lerat, A. Three Dimensional Calculation of Transonic Viscous Flows by an Implicit Method / A. Lerat, R. Peyret // The American Institute of Aeronautics and Astronautics Journal.-1985.-V. 23.-№ 11.-PP. 1670- 1678.

31. Lockhart, R.W. Proposed Correlation of Data for Isothermal Two-Phase, Two-Component Flow in Pipes / R.W. Lockhart, R.C. Martinelli // Chemical Engineering

32. Progress. 1949. - V. 45. - № 1. - PP. 39-45.

33. Lopez, D. Performances of transient two-phase flow models / D. Lopez, P. Duchet-Suchaux // Proceedings of SPE International Petroleum Conference (Villahermosa, Mexico, 3-5 March 1998). Villahermosa, 1998. - PP. 227-242.

34. Martinelly, R.C. Isothermal pressure drop for two-phase two-component flow in a horizontal pipe / R.C. Martinelly, L.M.K. Boelter, T.H.M. Taylor, E.G. Tomson, Moen R. H. // Trans. Amer. Soc. Mech. Engrs. 1944. -V. 66. -№ 2. - PP.139-151.

35. Minami, K. Transient two-phase flow behavior in pipelines: experiment and modeling / K. Minami, O. Shoham // International Journal of Multiphase Flow. 1994. - V. 20. - № 4. - PP.739-752.

36. Moore, T.V. Experimental measurement of slippage in flow trough vertical pipes / T.V. Moore, H.D. Wilde // Trans. Amer. Inst. Mining and Met. Engrs. 1931. - № 92. - PP. 296-313.

37. Norwegian Petroleum Technology. A success story : collected articles / Edit, by H. Keilen. Trondheim : NTVA, 2005. - 89 pp.

38. Petalas, N. Mechanistic Model for Multiphase Flow in Pipes / N. Petalas, K, Aziz // Journal of Canadian Petroleum Technology. 2000. - V. 39. - № 6. - PP. 43-55.

39. Sylvester, N.D. A mechanistic model for two-phase vertical slug flow in pipes / N.D. Sylvester // Journal of Energy Resources Technology. 1987. - V.109. - №4. - PP. 206213.

40. Taitel, Y. A model for predicting flow regime transitions in horizontal and near horizontal gas-liquid flow / Y.A. Taitel, A.E. Dukler, // AIChE Journal. 1976. - V. 22. -№ 1. - PP. 47-55.

41. Taitel, Y. Simplified transient solution and simulation of two-phase flow in pipelines / Y. Taitel, O.Shoham, J.P. Brill // Chemical Engineering Science. 1989. - V. 44. - № 6. -PP.1353-1359.50.