CC BY
164
Ключевые слова:
эксплуатация
скважин,
поздняя стадия
разработки,
обводнение
скважин,
минимальный
дебит,
газожидкостные
потоки,
забой,
глубина спуска лифтовых труб.
Keywords:
well operation, later stage of development, drowning gas wells, the minimum gas flow rate, gas-liquid flows, borehole bottom, setting depth of tubing.
УДК 622.279.23
О.В. Николаев, С.А. Шулепин
Влияние глубины спуска лифтовых труб на работу обводненной газовой скважины
На поздней стадии разработки газовых месторождений имеет место обводнение скважин, что приводит к необходимости создания и применения расчетных моделей, учитывающих влияние жидкости в продукции на технологические показатели работы скважин, кустов, сборных пунктов (УКПГ) и промысла в целом. В работах [1, 2] рассмотрены две основные группы методик расчетов техрежимов скважин, содержащих в продукции жидкость. Первая группа основана на использовании математических выражений для минимальной скорости газа, необходимой для выноса жидкости в вертикальном потоке. Существует ряд таких соотношений, теоретически выведенных разными авторами исходя из рассмотрения различных физических условий, имеющих место в газожидкостных потоках и отражающих определенные аспекты взаимодействия газовой и жидкой фаз между собой и твердой поверхностью стенки трубы. Основой второй группы является метод сопряженных элементов, рассматривающий совместную работу двух соседствующих звеньев цепи - пласта и лифтовой трубы, соединенных забоем в качестве общей точки, в которой определенные физические параметры газожидкостного потока принимаются одинаковыми для обоих элементов. Методики второй группы требуют знания характеристики лифтовой трубы для газожидкостных потоков и отличаются между собой в основном выбором способа задания в ней потерь давления.
Оба подхода имеют общий недостаток - они игнорируют физические явления на забое. Однако при определенных условиях технологический режим газовой скважины может в значительной степени определяться процессами, происходящими на забое скважины. Дело в том, что даже в устойчивом режиме работы скважины, т.е. при полном выносе жидкости потоком газа по лифтовой трубе на поверхность, пространство забоя заполнено газожидкостной смесью с изменяющимся по высоте соотношением фаз, что приводит к наличию на участке скважины значительного вертикального градиента давления. В связи с этим при больших расстояниях между башмаком лифтовой трубы и нижней границей интервала перфорации локальная депрессия с глубиной может заметно уменьшаться, что в условиях поздней стадии разработки, характеризуемой пониженным пластовым давлением, может приводить к существенному снижению рабочего дебита. Отметим, что газожидкостная смесь на забое может формироваться как пластовой, так и конденсационной жидкостью по мере ее выпадения в стволе в связи с изменением термобарических условий вдоль направления движения извлекаемого газа.
Рассмотрим газожидкостные процессы, которые имеют место на забое обводненной скважины, и их влияние на технологический режим. В расчетах используем корреляции для потерь давления в вертикальных газожидкостных потоках, полученные на специализированном стенде при давлениях от 0,2 до 3,0 МПа для труб диаметром от 0,62 до 0,15 м [3-5]. Результаты экспериментов позволили получить зависимости потерь давления от диаметра лифтовой трубы, что дало возможность путем экстраполяции рассчитывать процессы в трубах диаметром 21 см (типичный внутренний размер эксплуатационной колонны на западно-сибирских месторождениях).
Вначале рассмотрим вариант, когда башмак лифтовой трубы располагается на уровне кровли пласта, затем перейдем к изучению вопроса о том, как глубина спуска насосно-компрессорной трубы (НКТ) влияет на технологический режим работы скважины. Будем рассматривать совершенную по характеру и степени вскрытия скважину.
№ 4 (15) / 2013
Проблемы эксплуатации газовых, газоконденсатных и нефтегазоконденсатных месторождений
47
Вариант 1
Башмак лифтовой трубы расположен на уровне кровли продуктивного пласта (рис. 1а).
Для расчетов технологического режима работы скважины выделим на забое два возможных характерных интервала. В интервале h12 находится однофазная жидкость, газ в этом интервале из пласта в скважину не поступает. В данной работе будем считать, что обратной фильтрации жидкости в пласт не происходит. В интервале h11 находится жидкость, барбо-тируемая газом, поступающим из пласта. Эта двухфазная смесь поднимается вверх до баш-
а
мака лифтовой трубы, поступает в лифтовую трубу и течет к устью скважины.
Разобьем перфорированный интервал высоты H на N одинаковых ячеек достаточно малого по высоте размера АН = H/N. Занумеруем ячейки (начиная снизу) и рассмотрим j-ю ячейку (рис. 2). Газ в нее поступает из двух источников: из нижней ячейки в количестве Q(j _ ц и из пласта в количестве Qj а утекает в верхнюю ячейку в количестве Q(J). Примем, что расходное количество жидкости определяется водогазовым фактором W = qx2\/Q.
б
Рис. 1. Схематичное изображение забоя скважины: а - вариант 1 - башмак лифтовой трубы расположен на уровне кровли; б - вариант 2 - башмак лифтовой трубы спущен ниже кровли. Rc - радиус скважины на забое; d - диаметр лифтовой трубы; Н - интервал вскрытия пласта; h11 - интервал, заполненный газожидкостной смесью; h12 - интервал, заполненный жидкостью; h1 - интервал от башмака до нижней точки забоя; h2 - расстояние от кровли до башмака
Рис. 2. Схема расчета восходящего газожидкостного потока на забое скважины
№ 4 (15) / 2013
48
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
Определим следующие соотношения для каждой ячейки. Для притока газа в ячейку j из пласта:
Рт( j) Рзаб (j) a(j ')Qm(j') + b(j )Qm( j)> (1)
где a(j) = aN, b(f) = bN2; a, b - коэффициенты фильтрационного сопротивления пласта; pmj) - пластовое давление напротив j-й ячейки; p3a6j) - давление в j-й ячейке.
Для притока газа в ячейку j из нижележащей ячейки (j- 1) воспользуемся формулами, полученными для вертикальных газожидкостных потоков [4, 5]:
2 2 _
Рзаб (j ) Рзаб (j-1)
= р(j-1) S АН(j-1) + +(j-1) рж &Н(j-1), (2)
где - ц - плотность газа в (j - 1)-й ячейке; g - ускорение свободного падения.
Безразмерные потери давления i = i(Fr*, Bu) определяются по эмпирическим корреляциям, в которых в качестве аргументов выступают модифицированный параметр Фруда и параметр Бузинова [5]:
Fr"
2
. _Ц_
Р* gDc
(3)
Bu =
( 2 V'
О w
V РжgDc У l gDc J
1/2
(4)
где Dc - диаметр забоя; рг, рж - плотности газа и жидкости; u, w - средние по сечению потока скорости газа и жидкости; с - величина поверхностного натяжения на границе «газ - жидкость».
Для перетока газа в верхнюю ячейку справедливо соотношение
Рзаб (J+1) Рзаб (j ) P(j) 8 j) + '( у)Рж ^Я( J) ■ (5)
Соотношение материального баланса для газа в каждой ячейке выражается как
Q( j) Q( 1 -1) + Qrm( j) • (6)
Материальный баланс жидкости в ячейке определяется аналогичным образом:
рами: давлением, плотностью газа, расходом жидкости, расходом газа, модифицированным параметром Фруда, параметром Бузинова, потерями на трение, объемным водосодержанием и т.д. Задаваясь величиной давления в нижней точке забоя, по соотношениям (1)-(7) можно рассчитать давление на башмаке лифтовой трубы и дебит скважины. Расчет потерь давления в лифтовой трубе осуществляется по методике, изложенной в работе [4].
Поскольку в скважине между башмаком лифтовых труб и нижней точкой забоя находится не газ, а барботируемая жидкость, градиент давления по вертикальной координате больше по величине, чем в случае для однофазного газа. Это приводит к тому, что «локальная» депрессия
^P(J) Роли) Рзаб (j)
(8)
в каждом из интервалов AHj уменьшается с глубиной, и в нижнем интервале АН1 она минимальна. Возможны режимы, когда в нижней части забоя находится однофазная жидкость (интервал hl2 на рис. 2). Такие режимы могут реализовываться при достаточно большом промежутке между башмаком лифта и нижней точкой забоя (более 50 м). В этом случае работающий интервал кп меньше перфорированного, что позволяет рассматривать скважину как несовершенную по степени вскрытия.
Воспользуемся известными соотношениями для расчета коэффициентов фильтрационного сопротивления a* и b* в несовершенных по степени вскрытия скважинах[6]:
1 + -
ln( Rk / Rc)
(9)
b* = b(l + C2),
(10)
где С1 и С2 - коэффициенты, учитывающие степень вскрытия пласта.
Обозначим у = ки1И. Поскольку у - аналог относительного вскрытия пласта, коэффициенты C1 и C2 можно определить из соотношений
Cj = —ln y + -—У ln
У У
Г Rk_ " V RC J
(11)
^ж (j) qж(j-1) + qжm(j)
(7) c2 =
У
(12)
В каждой ячейке устанавливается динамическое равновесие с присущими ему парамет-
где Rk - радиус контура питания; Rc - радиус скважины в интервале перфорации.
№ 4 (15) / 2013
Проблемы эксплуатации газовых, газоконденсатных и нефтегазоконденсатных месторождений
49
Расчетные коэффициенты a(J) и в (1) определим по а* и Ь*. Будем считать, что давление на контуре питания распределяется по толщине пласта в соответствии с барометрической формулой, и средневзвешенная ее величина определяет значение пластового давления.
Проведенные расчеты показывают, что, как и следовало ожидать, при одинаковых характеристиках пласта скважина с большим интервалом перфорации на большом расстоянии от башмака лифтовой трубы до нижней точки забоя работает в менее благоприятных условиях; предельное значение устьевого давления у такой скважины выше, а минимальный дебит больше.
Параметры, принятые при расчете, отражены в табл. 1.
На рис. 3 изображены кривые продуктивности скважины для разных величин интервала перфорации. На рис. 4 и 5 приведены результаты расчетов параметров предельного режима от протяженности интервала перфорации. Из анализа графиков следует, что высокий подъем башмака относительно забоя в обводненной скважине может приводить к существенному ухудшению условий эксплуатации скважины. Особенно значительно это отрицательное влияние сказывается на поздней стадии разработки, когда необходимое для устойчивой эксплуатации скважины снижение устьевого давления требует дополнительных затрат, связанных, например, с вводом следующей очереди дожимной компрессорной станции (ДКС).
Таблица 1
Пластовое давление, рпл, МПа Размеры лифтовой трубы, м Радиус забоя, м Радиус зоны дренирования, м Коэффициенты фильтрационного сопротивления пласта
L D а, 10-2 (МПа2 • сут)/тыс. м3 Ь, 10-2 ((МПа • сут)/тыс. м3)2
2,0 1000 0,15 0,21 500 0,1 0,001
Рис. 3. Влияние высоты башмака лифтовой трубы относительно забоя H на продуктивность скважины
Рис. 4. Зависимость минимального дебита от высоты башмака лифтовой трубы над нижней границей интервала перфорации
№ 4 (15) / 2013
50
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
Рис. 5. Зависимость максимальных устьевого давления и давления на башмаке от протяженности интервала перфорации
х
Рис. 6. Зависимость дебита скважины от глубины спуска лифтовой трубы при разных значениях устьевого давления
Вариант 2
Башмак лифтовой трубы расположен ниже уровня кровли продуктивного пласта (рис. 1б).
Введем обозначения:
х = h2IH - степень спуска башмака лифтовой трубы ниже кровли пласта (или безразмерный интервал перекрытия лифтовой трубой интервала перфорации);
y = (h11 + h2)IH - безразмерный работающий интервал (аналог относительного вскрытия пласта в теории несовершенных скважин);
z = h11Ih1 - безразмерный работающий интервал в пространстве ниже башмака лифтовой трубы;
yg = h2I(h2 + h11) - безразмерный интервал с нисходящим потоком.
Для исследования влияния глубины спуска лифтовой трубы на технологический режим работы скважины рассмотрим три характерных интервала на забое. Как и в предыдущем варианте, имеются два интервала h12 и h11, располо-
женные ниже башмака лифтовой трубы. Эти интервалы описаны выше. Дополнительно к ним появляется третий промежуток h2 от башмака лифтовой трубы до верхних дыр перфорированного интервала, в котором газ (возможно с жидкостью) поступает из пласта в кольцевое пространство и течет вниз до башмака, где смешивается с восходящим из нижнего участка забоя газожидкостным потоком, заходит в лифтовую трубу и направляется к устью скважины. Расчеты потерь давления в нисходящем потоке кольцевого пространства интервала h2 проводятся по известным гидродинамическим соотношениям для однофазного потока. Расчеты потерь давления в восходящем газожидкостном потоке интервала h11 проводятся аналогично варианту 1.
Исходные данные для расчетов представлены в табл. 1. Интервал перфорации был принят равным 100 м. На рис. 6 представлены значе-
№ 4 (15) / 2013
Проблемы эксплуатации газовых, газоконденсатных и нефтегазоконденсатных месторождений
51
ния дебитов скважины при различных значениях устьевого давления в зависимости от глубины спуска лифтовой трубы. Согласно рисунку, при устьевых давлениях выше 1,55 МПа и башмаке лифтовой трубы, установленной на высоте верхних дыр перфорации, скважина устойчиво работать не будет. Для работы скважины башмак лифтовой трубы должен быть спущен ниже; например, при устьевом давлении 1,72 МПа скважина будет работать только в случае, если башмак располагается от забоя не выше, чем на 20 м.
На рис. 7 отражены зависимости значений предельного устьевого давления от глубины расположения башмака в случае разных значений водогазового фактора. Глубина спуска лифтовой трубы существенно влияет на режим работы обводненной газовой скважины. Расчет методом сопряженных элементов, т.е. без учета процессов на забое, для условий рассмотренного примера дает значения параметров предельного режима, приведенные в табл. 2. Как видно из сравнения с рис. 6, эти величины соответствуют вариантам значительного спуска лифтовой трубы к нижней точке забоя (X > 0,8).
Таблица 2
w, м3/м3 Qmin’ тыс. м3/сут Руст max’ МПа
1 166 1,695
10 173 1,627
Как отмечено выше, рассматриваемые процессы на забое могут сопровождаться так называемой обратной фильтрацией жидкости в пласт [7]. В данной работе этот эффект не учитывался, поскольку целью авторов статьи являлась общая демонстрация методики, основанной на последних экспериментальных исследованиях газожидкостных потоков. Очевидно, в предложенной расчетной схеме (см. рис. 2) можно задать условия оттока жидкости в пласт, что и предстоит сделать по мере развития методики. Для дальнейшего развития теории эксплуатации скважин явление обратной фильтрации, несомненно, должно быть изучено экспериментально и теоретически и учтено в расчетных моделях, поскольку обратная фильтрация может существенно влиять на величину водосодержания на забое в сторону ее уменьшения и, как следствие, приводить к значительному улучшению условий извлечения газа.
Таким образом, предложена новая методика расчета режима работы обводненных скважин, учитывающая наличие газожидкостной смеси на забое и основанная на результатах экспериментальных исследований газожидкостных потоков. Разработанная расчетная схема позволяет получать более точные значения рабочих и предельных параметров работы скважин по сравнению с существующими методиками, поскольку учитывает дополнительные важные факторы, определяющие эксплуатационные характеристики.
X
Рис. 7. Зависимость предельно максимальной величины устьевого давления от глубины спуска лифтовой трубы при разных значениях водогазового фактора
№ 4 (15) / 2013
52
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
На поздней стадии разработки месторождений режим работы скважины зависит от глубины спуска лифтовых труб. Эта зависимость обусловлена повышенным (по сравнению с пластом и лифтовой трубой) содержанием подвижной жидкой фазы на забое, что создает дополнительное сопротивление газовому потоку и приводит к перераспределению давлений по высоте интервала перфорации.
Глубина спуска лифтовой трубы влияет не только на текущие параметры рабочего ре-
Список литературы
1. Брилл Дж.П. Многофазный поток в скважинах / Дж.П. Брилл, Х. Мукерджи. - М.; Ижевск:
Ин-т комп. иссл., 2006. - 384 с.
2. Изюмченко Д.В. Газожидкостные потоки в вертикальных трубах: парадоксы гидродинамики / Д.В. Изюмченко,
О.В. Николаев, С.А. Шулепин // Вести газовой науки: Проблемы эксплуатации газовых, газоконденсатных и нефтегазоконденсатных месторождений. - М.: Газпром ВНИИГАЗ,
2013. - С. 36-45.
3. Бузинов С.Н. Экспериментальные исследования движения двухфазных систем в газовых скважинах / С.Н. Бузинов, С.А. Бородин, В.М. Пищухин и др. // Георесурсы. - 2010. - № 4. - С. 55-58.
4. Бузинов С.Н. Расчет потерь давления
в газовых скважинах на поздней стадии разработки месторождений / С.Н. Бузинов,
Г.М. Гереш, С.А. Бородин и др. // Газовая промышленность. - 2011. - № 12. - С. 18-21.
жима, но и на границы устойчивого режима, в частности на величину минимального дебита скважины.
Влияние глубины спуска лифтовых труб наиболее заметно при значительном интервале вскрытия пласта (более 100 м), поскольку наличие на забое протяженного столба газожидкостной смеси с повышенным содержанием жидкой фазы создает предпосылку для ухудшения условий поступления газа из пласта в скважину.
5. Николаев О.В. Регулирование работы газовых скважин на завершающей стадии разработки залежей по результатам экспериментальных исследований газожидкостных потоков
в вертикальных трубах: дис. канд. техн. наук / Николаев Олег Валерьевич. - М., 2012. - 128 с.
6. Зотов ГА. Расчет фильтрационных сопротивлений скважины, несовершенной по степени вскрытия пласта, при нелинейном законе фильтрации / Г. А. Зотов // Подземная гидродинамика: тр. ВНИИГАЗа. - М.: ГТТИ, 1963. - Вып. 18/26. - С. 64-70.
7. Плосков А.А. Фильтрация конденсационной воды в пласт в работающих скважинах сеноманских залежей / А.А. Плосков,
И. В. Шулятиков, В. З. Минликаев и др. // Газовая промышленность. - 2013. - № 5. -С. 62-66.
№ 4 (15) / 2013
