Экспериментальные исследования движения двухфазных систем в газовых скважинах Текст научной статьи по специальности «Физика»

УДК: 532.546 : 622.279.5

С.Н. Бузинов1, С.А. Бородин1, В.М. Пищухин1, А.Н. Харитонов2,

О.В. Николаев1, С.А. Шулепин1

'ООО «Газпром ВНИИГАЗ», Москва 2000 «Газпром добыча Надым», Надым [email protected], [email protected]

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

СКВАЖИНАХ

Одна из проблем истощенных газовых месторождений - накопление жидкости в стволе скважины, и как следствие остановка ее работы. Традиционно применяемые методы расчета оптимального режима работы скважин при осложненных условиях по эмпирическим формулам имеют ряд ограничений и допущений. В статье рассмотрен оптимизированный алгоритм расчета, опирающийся на результаты практических экспериментов, проведенных на экспериментальной установке ООО«ГазпромВНИИГАЗ» - стенд по отработке эксплуатации газовых скважин на поздней стадии разработки.

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

Эксплуатация газовых скважин сопровождается поступлением в ствол скважины пластовой и конденсационной воды. Гидравлические процессы при движении сухого газа и газа, содержащего жидкость, во многих случаях существенно отличаются между собой. Жидкость, находящаяся в стволе скважины, как правило, усложняет условия работы скважины, вызывая дополнительные потери сопротивления (Точигин, 1979). Как показывает практика разработки газовых месторождений, а также экспериментальные исследования, наибольшее отрицательное воздействие оказывает жидкость на работу скважины при малых расходах газа что, как правило, наблюдается на поздней стадии разработки. Это актуализирует проблему изучения механизма двухфазного течения газа и жидкости в трубах (Ли и др., 2008).

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

В настоящей статье ставится более скромная задача -

Г

0,04

0,03

0,02

0,01

~~| у = 0,0101х

« №314 (15 ата) • №323 (10 ата) №325 (15 ата) • №326 (6.5 ата) • №334 (5-9 ата)-Сухой воздух — Линейный (Сухой воздух )

0,5 1

1,5

2,5

3,5

4

Fr ■

Рис. 1. Характеристики лифтовой колонны 114 мм при расходе воды 68 л/ч и давлениях 5 ата, 6.5 ата, 9 ата, 10 ата, 15 ата и сухом воздухе.

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

-*-а -•- б -"-в — г —д —е -

1000 2000 3000 4000 Fr*, Фруд приведенный

Рис. 2. Зависимость индикаторныгх характеристик лифта < = 62 мм от расхода водыг. Расход водыг: а - 210 л/ч, б -130 л/ч, в - 67л/ч, г - 4 л/ч. д - сухой воздух, е - линейный (сухой воздух).

-•- а -■- б В

-*- г -ж- д -е

3000 4000 Fr*, Фруд приведенный

I

Рис. 3. Зависимость индикаторныгх характеристик лифта < = 76 мм от расхода водыг. Расход водыг: а - 350 л/ч, б -120 л/ч, в - 80 л/ч, г - 2.5 л/ч. д - сухой воздух, е - линейный (сухой воздух).

4 (36) 2010

^научно-техническим журнал

Георесурсы

Нами на установке (изложено в работе (Тер-Сарки-сов, 2001)), а также в настоящем сборнике (Бузинов и др., 2010), были проведены экспериментальные исследования. Кроме того, были использованы ранее проведенные в ООО «Газпром ВНИИГАЗ» экспериментальные исследования, проведенные Б.Г. Ахмедовым (Ахмедов, 1982), а также, к сожалению, неопубликованные исследования О.В. Николаева и А.А. Михайлова.

Исследования проводились применительно к условиям разработки сеноманских газовых залежей Надым-Пу-ровского газоносного района при сниженных пластовых давлениях.

Опыты проводились на вертикальных трубах с внутренним диаметром 62 мм, 76 мм, 100 мм, 152 мм. Давление при этом варьировалось 0,3 МПа до 1,5 МПа. Максимальный расход газа при этом составлял 350 м3/ч. Водный фактор изменялся от 0 до 1300 л/ч. В процессе исследования замеряется расход газа и воды, а также перепад давления на трубе длинной 29 м.

Результаты экспериментальных исследований обрабатывались графически в координатах i=Ap/(pgH),

pr = . _Рг . * «*> Р

где i - относительные потери давления,

Ар - потери давления в трубе на расстоянии Н, g - ускорение свободного падения, Н - высота колонны - 29 м, .г - приведенный параметр Фруда, V - скорость газа при рабочих условиях, Б - внутренний диаметр лифтовой колонны, рг - плотность газа при рабочих условиях, рж -плотность жидкости.

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

Здесь иллюстрируется зависимость г от .г, полученные на трубе ё =100 мм при различных значениях давлений в трубе. Как видно, наблюдается хорошее совпадение результатов исследования для разных значений давлений и других параметров.

На рисунках 2 - 4 помещены графики зависимости величин .г и г для труб с тремя внутренними диаметрами 62 мм, 76 мм, 100 мм. Общим является хорошо прослеживаемая зависимость увеличения потерь давления в любой трубе от увеличения количества воды (расхода) в ней. Черная линия - линия тренда, экстраполирующая линию, описывающую потери сопротивления в трубе для сухого газа. Следует отметить, что при малых расходах воды наблюдаются потери давления меньшие, чем при движении сухого газа, вплоть до расхода воды равного 10 л/ч для трубы 114 мм.

Как видно при некотором значении Бг* функция 1(Бг*, Бг ) имеет минимум. Зная зависимость г от Бг*, Бг , а так-

ж' * 7 ж7

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

Перепад давления в лифте (разница давлений на забое и на устье) можно определить путем численного решения дифференциального уравнения:

dP

aZ ,

где ось Ъ - направлена вверх.

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

pL~pL.=AQ+BQ2

(2)

можно построить зависимость забойного давления от дебита (Одишария, Точигин, 1998). Эта зависимость называется индикаторной характеристикой скважины. Пример индикаторной характеристики скважины и газожидкостного подъемника приведен на рис. 5.

Пересечение индикаторных характеристик позволяет определить рабочий дебит скважины. Устойчивая работа скважины на левой ветви не возможна. В связи с этим минимально-возможный устойчивый дебит скважины может быть, когда обеспечивается значение приведенного числа Фруда, соответствующее минимуму индикаторной характеристики лифта. Обозначим это значение через .г . Тогда минимально-возможный дебит скважины будет равен значению, соответствующему точке пересечения правой ветви индикаторной характеристики лифта с индикаторной характеристикой пласта. Например, на рис. 5 это значение приблизительно равно 112 м3/сут.

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

* кр

Как видно из рисунков 2 - 4, при значениях ¥г >¥т

пр. кр.

наблюдается линейная связь между величиной г и Бг*.

Экспериментально было также установлена зависимость величины 1 от параметра Фруда для жидкости, ¥т , для труб диаметром ё = 62 мм, 76 мм, 100 мм. Для значений ¥т >¥т была получена эмпирическая зависимость

пр. кр. * г

г от Бг* и приведенного расхода для жидкости. Под приведенным расходом воды, д, понимается расход, при котором значение числа Фруда для жидкости будет равно числу Фруда для трубы ё =62 мм и определяется формулой:

= Яв

62

А у

(3)

где д - реальный расход жидкости, в л/ч; Б - внутренний диаметр трубы в мм.

Эмпирическая (аппроксимационная) формула имеет вид:

I = -0.000877 + 0.000276 ■ д - 0,547 • д2 ■ 10^ • дв

+ 0A0\4-Frnp -10

-в

(4)

(1)

Эту формулу можно рекомендовать для проведения ориентировочных расчетов для труб, внутренний диаметр которых изменяется от 62 мм до 100 мм. Более правильно необходимо ориентироваться на экспериментальные исследования, которые должны быть проведены на сертифицированном стенде.

|— научно-технический журнал

I еоресурсы 4 (36) 2010

к 0,06

5 X

а>

5 0,05

со

Ч

g. 0,04

I

■= 0,03

5 0,02 0,01

1 /

\

— а б в г Д — е

1000 2000 3000

Fr*, Фруд приведенный

Рис. 4. Зависимость индикаторных характеристик лифта ^ = 100 мм от расхода водыг. Расход водыг: а - 132 л/ч, б - 97 л/ ч, в - 68 л/ч, г - 9 л/ч. д - сухой воздух, е - линейный (сухой воздух).

(U

с:

39 34 29 24 19 14 9

--Инд. кривая плас -Инд. кривая лиф

та

| 0 50 100 150 200 250

Дебит газа, куб.м/сут. I

Рис. 5. Индикаторныге характеристики лифтовой колонныг и пласта.

На рисунке 6 приведено сопоставление графика зависимости, (г = ¿(Уг^, двИй,), полученного в результате проведенных экспериментов на трубах ё =76 мм и ё =100 мм при одинаковые значениях приведенного расхода воды. Как видно, получается удовлетворительное сопоставление зависимостей, отображающих процесс при одинаковых значениях приведенного расхода воды. Отклонения составили не более 5 - 6 %. По графикам видно, что индикаторная кривая лифта ё =76 мм (зеленого цвета) очень близка с кривой ё =100 мм (малинового цвета) при одинаковых значениях приведенного расхода воды 27 л/ч.

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

На рисунке 7 приведено сопоставление экспериментальных данных, полученных нами с данными расчета по аппроксимационной формуле. Как видно получается удовлетворительное сопоставление. Отклонения от расчетной формулы (общей для всех перечисленных в статье диаметров труб) составило не более 8 %.

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

Следует отметить, что необходимое сопоставление с промысловыми данными имеет одну особенность: необходимо не только (а может быть и не столько) точ-

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

В заключении отметим еще два момента. О.В. Николаевым и A.A. Михайловым проведены экспериментальные исследования, когда в качестве агента при газлифте применяли гелий, плотность которого в 7,5 раз меньше плотности воздуха. Зависимость гидравлических потерь от приведенного параметра Фруда практически полностью совпала с аналогичной зависимостью, полученной в экспериментах с воздухом.

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

Констатируется, ссылаясь на исследования Данса и Роса, о наличии 8 безразмерных параметров, из которых 4 являются определяющими. Приведем их (в интерпретации авторов).

- Показатель скорости жидкости

N =V -4

lylv r sl i

Pl

0,027

0,023

S.®

? I

0,019

tu s 31 s

8 0,015

0,011

/

-20.4/100.4 -27.6/100.4 -27/76 21/76

300

700

1100

1500 1900 Число Фруда

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

—270—270теор—210- 210теор — 69 — 69теор —48 — 48теор 27 27теор 21 21тбор 2 2теор

Рис 7. Графики сравнения сопоставления расчетным формул с эксперименталыныгми данныгми.

4 (36) 2010

^научно-техническим журнал

Георесурсы

- Показатель скорости газа

- Показатель диаметра трубы

Pl

Nd=d- 4 М^-

\g°l

- Показатель вязкости жидкости

Pl

где V и - средние скорости (расходные) движения газа и жидкости в трубе, р - плотность жидкости, аь - поверхностное натяжение, р - вязкость жидкости, g - ускорение свободного падения.

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

N2 V2

Fr = LV = sl

ж Nd gD

Fr =

N:

gv

N,

-iL gD

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

V2

метром является не параметр Фруда, —— , а приведенный

' & '

параметр Фруда, Fr* = ■ , где р - плотность газа.

V2

Ф Рь

Данс и Рос проводили исследования при условиях, близких к атмосферному. В связи с этим, параметр р^рь был постоянным и в этих условиях определяющими параметрами были параметр Фруда для жидкости и параметр Фруда

V2

для газа, Рг* = .

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

представить в виде Fr

Ks Р-Р-Тст.

тельная плотность газа по отношению к воздуху.

Представление определяющего параметра Фруда в

У2 Р

форме рг* = ---— более естественно. Подъемная спо-

2 Ф Рь

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

Литература

Точигин A.A. Гидродинамика газожидкостных смесей в разработке и эксплуатации газоконденсатных месторождений. Авто-реф. дисс. уч. ст. д.техн.н. М. 1979.

Ли Дж., Никенс Г.В., Уэллс М. Эксплуатация обводняющихся газовых скважин. М. ООО «Премиум Инжинириг». 2008.

Тер-Саркисов Р.М. Газовая промышленность. 2001.

Бузинов С.Н., Бородин С.А., Пищухин В.М., Шулепин С.А. Стенд по отработке технологии эксплуатации скважин на поздней стадии разработки месторождений. 2010.

Ахмедов Б.Г. Оптимизация режимов работы газовых скважин с помощью ПАВ. Дисс. на соискание уч. ст. к.г.-м.н. 1982.

Одишария Г.Э., Точигин A.A. Прикладная гидродинамика газожидкостных смесей. М.: Всероссийский научно-исследовательский институт природных газов и газовых технологий. Ивановский государственный энергетический университет. 1998.

Мукерджи X. Производительность скважин (Руководство). М. 2001.

S.N. Buzinov, S.A. Borodin, V.M. Pischukhin, A.N. Kharitonov, O.V. Nikolaev, S.A. Shulepin. Experimental researches of diphasic systems movement in gas wells.

One of the problems of depleted gas fields is liquid accumulation in a well trunk and as consequence a stop of its work. Traditionally applied methods of calculation of an optimum operating mode of wells under the complicated conditions on empirical formulas have a number of restrictions and assumptions. In article the optimized algorithm of calculation leaning against results of practical experiments, «GazpromVNIIGAZ» Ltd spent on experimental installation - the stand on working off of operation of gas wells at a late stage of working out is considered.

Keywords: dynamics of fluids and gases, the diphasic current, the resulted (brought) parameter of Fruda, the exhausted deposits, is minimum-admissible wells debit, wells itself rate termination.

Олег Валерьевич Николаев ведущий инженер. Научные интересы: физическое моделирование процессов многофазовой фильтрации.

ООО «Газпром ВНИИГАЗ»

142717, РФ, Пос. Развилка, Ленинский р-н, Московская область. Тел.: (495)355-96-46.

gD

Z-Pm, Т

ст.

где р - относи-

Андрей Николаевич Харитонов к.тех.н., Зам. начальника ИТЦ Надым ООО «Газпром добыча Надым». Научные интересы: проблематика эксплуатации газоконденсатных месторождений на поздней стадии разработки.

629730, РФ, Ямало-Ненецкий автономный округ, г.Надым, ул. Полярная, 1/1 (инженерно-технический центр -ИТЦ). Тел./факс: 8 (3499) 56-70-08.

^^р«» •— научно-технический журнал

I еоресурсы 4 во; 2010