Расчет процесса фонтанирования с помощью кривых распределения давления
Умение рассчитывать при любых заданных условиях кривую распределения давления вдоль НКТ при движении по ним газожидкостной смеси позволяет по-новому подойти к расчету процесса фонтанирования, выбора диаметра труб и режима в целом. Использование кривых распределения давления Р(х) при проектировании и анализе фонтанной эксплуатации (а также других способов эксплуатации скважин) позволяет решить ряд новых задач, недоступных при использовании прежних расчетных методов. Далее будем исходить из того, что при любых заданных условиях кривая распределения давления Р(х) в НКТ может быть определена и построена любыми возможными методами.
Заметим, что для проектирования или для анализа фонтанной эксплуатации не требуется распределение давления Р(х) вдоль всей длины НКТ. Достаточно знать забойное или башмачное давление, соответствующее данному забойному давлению, давление на устье при заданных параметрах работы скважины или наоборот, устьевое давление и соответствующее давление на забое при заданных параметрах работы скважины.
Однако поскольку простых и надежных формул (кроме формул А. П. Крылова), связывающих устьевое и забойное давления при прочих заданных условиях, нет, то приходится прибегать к численному интегрированию процесса движения ГЖС по трубе, т. е. расчету по шагам. При таком решении неизбежно получаются значения давлений в промежуточных точках между устьем и забоем, использование которых необязательно. Рассмотрим для начала простейший случай, когда задан дебит скважины Q и соответствующее этому дебиту забойное давление Рс.Отметим, что во всех случаях проектирования процесса эксплуатации скважины любым способом знание уравнения притока или индикаторной линии обязательно. В противном случае любой инженерный расчет становится невозможным, если не говорить о предположительных оценках возможных показателей работы скважины. Итак, если задан дебит, то по индикаторной линии или по уравнению притока определяется соответствующее этому дебиту давление на забое скважины.
В отношении фонтанных труб уже указывалось, что их диаметр выбирается из соображений технологических условий и возможности спуска в скважину глубинных приборов для различных исследований. Можно сказать, что для подавляющего числа случаев это будут либо трубы диаметром d = 60 мм, либо d = 73 мм. Лишь для редких случаев, когда ожидаемые отборы могут достигать нескольких сот м3/сут, можно говорить о целесообразности использования труб d = 89 мм. Во всяком случае для последующего расчета диаметром НКТ задаемся.
Зная дебит, газовый фактор, плотность нефти, воды и обводненность продукции, а также другие данные, такие как температура и ее распределение по стволу скважины, объемный коэффициент нефти (жидкости), необходимые для расчета, строим кривую распределения давления Р(х), начиная от точки с известным давлением Рс на забое скважины (рис. 8.5).
Рис. 8.5. Построение кривой распределения давления в фонтанных трубах
по методу «снизу вверх» и определение давления на устье
При этом могут возникнуть разные условия расчета, которые необходимо учитывать.
а. Башмак НКТ находится непосредственно на забое скважины, так что Рс = Рб.
б. Башмак НКТ находится выше забоя на некотором расстоянии а = Н - L, так что Рб < Рс.
в. Давление на забое или у башмака больше давления насыщения, т. е. Рс = Рб > Рнас.
г. Давление на забое меньше давления насыщения, т. е. Рс < Рнас.
Возможны также сочетание условия "а" с условиями "в" или "г", а также условия "б" с теми же "в" или "г". Предположим простейший случай: действуют условия "а" и "г". В этом случае ГЖС движется от башмака до устья, и расчет ведется по соответствующим формулам для газожидкостной смеси по шагам, начиная от башмака НКТ от точки с давлением Рс и до устья. Давление на устье получаем путем суммирования элементарных перепадов давления на n шагах:
. (8.55)
Если действуют условия "а" и "в", т.е. выделение газа начинается выше забоя в НКТ, то до точки Рнас расчет ведется по обычным формулам трубной гидравлики, с помощью которых определяются потери давления на трение.
Обозначим длину участка НКТ от забоя до точки с давлением Рнас, на котором будет двигаться однородная жидкость, через h (см. рис. 8.5). Тогда для этого участка запишется очевидное равенство давлений:
, (8.56)
где - гидростатическое давление столба жидкости высотою h и плотностью ρж;
- потери давления на трение при скорости жидкости С, м/с.
Подставляя значения Рг и Ртр в (8.56) и решая относительно h, получим
. (8.57)
Обычно второе слагаемое в круглых скобках знаменателя мало, поэтому им часто можно пренебречь.
На остальной длине НКТ, равной L - h, т. е. от точки давления насыщения и выше, будет происходить движение ГЖС, поэтому давление на устье будет равно
. (8.58)
Если действует условие "б", т. е. когда башмак НКТ выше забоя на величину a = H - L, то на этом участке при расчете распределения давления вместо диаметра трубы подставляется диаметр обсадной колонны.
Поскольку потери давления на трение из-за большого диаметра на этом участке малы, то ими всегда можно пренебречь. Давление на устье Ру определяется либо по формуле (8.55), либо по (8.58) в зависимости от того, выделяется ли газ с самого забоя (8.55) или НКТ (8.58).
Рассчитав кривую распределения давления и определив давление на устье скважины при заданном режиме ее работы, сопоставим вычисленную величину Ру с возможным давлением в выкидной линии Рл, по которому продукция скважины поступает в систему нефтегазосбора промысла. Если Ру > Рл, то работа скважины на рассчитанном режиме возможна, а избыточное давление на устье ΔРшт = Ру - Рл должно быть понижено созданием в арматуре устья дополнительного гидравлического сопротивления в виде регулируемого или нерегулируемого штуцера, в котором поток ГЖС дросселируется с давления Ру до давления Рл. Если при расчете окажется, что Ру < Рл, то фонтанирование скважины на проектируемом режиме будет невозможно. В таком случае необходимо задаться меньшим отбором Q, при котором давление на забое возрастает. Это в свою очередь приведет к более высокому давлению на устье скважины.
Изменяя отбор, а следовательно, и давление на забое, можно подобрать такие соотношения, при которых окажется Ру > Рл, когда фонтанирование будет возможно. Если ни одна комбинация Q и соответствующего Рс при построении кривой распределения давления Р(х) не дает давление на устье Ру > Рл, то фонтанирование такой скважины вообще невозможно.
Изложенная система расчета процесса фонтанирования может быть повторена для труб меньшего или большего диаметра для определения возможных режимов фонтанирования и дебита скважины при других диаметрах фонтанных труб.
Рассмотрим другой, наиболее общий случай, когда возникает необходимость определения всего комплекса возможных и невозможных условий фонтанирования скважины. При этом будем считать, что все проектируемые отборы жидкости из пласта допустимы и не противоречат принципам рациональной разработки залежи.
а. Задаемся несколькими забойными давлениями Рсi, лежащими в пределах Рmin < Рсi < Pпл, где Pпл - пластовое давление, a Pmin - минимальное давление на забое, при котором фонтанирование скважины заведомо неосуществимо.
б. Для принятых значений Pci определяем приток жидкости в скважину Qi по уравнению притока или по индикаторной линии.
в. Задавшись диаметром НКТ, рассчитываем распределение давления P(х) по методу снизу вверх для принятых значений забойных давлений Рci и соответствующих им дебитов Qi. В результате получаем i кривых Р(х) (рис. 8.6).
г. По полученным кривым Р(х) определяем i значений устьевых давлений Рyi.
д. Получаем систему данных, состоящих из нескольких забойных давлений Рсi, дебитов скважины Qi, и устьевых давлений Рyi.
Причем каждому конкретному давлению на забое Рci соответствует конкретный дебит и вычисленное давление на устье Рyi. Поскольку увеличение давления на забое Рci сопровождается уменьшением притока Qi и, как правило, увеличением давления на устье Рyi, то полученная система данных будет находиться в следующих соотношениях:
. (8.59
Рис. 8.6. Кривые распределения давления в фонтанном
подъемнике при нескольких (четырех) режимах работы
По полученным данным (8.59) можно построить две графические зависимости Q = f 1(Рc) и Ру = f 2 (Рс) (рис. 8.7). Графики отражают совместную, согласованную работу пласта и газожидкостного подъемника, общей точкой для которых является давление на забое скважины Рс. Отметим, что понижению давления Рс не всегда должно соответствовать уменьшение давления на устье Ру, как это показано на рис. 8.7.
Рис. 8.7. Согласование индикаторной линии (1) с зависимостью устьевого давления Ру от давления
на забое скважины Рс (2). Точки а - b разделяют возможные и невозможные режимы фонтанирования
Изменение Рс, сопровождаемое соответствующим изменением притока Q, приводит к изменению режима работы самого газожидкостного подъемника, который при определенных условиях может совпадать с режимами оптимальной или максимальной подач или иметь какой-то промежуточный режим. К.п.д. при этих режимах различный. Это может привести к различным зависимостям давления на устье от давления на забое и, в частности, к зависимостям, имеющим максимум или минимум. Это выявляется при расчете кривых Р(х).
На оси Ру можно отложить давление в выкидной линии Рл, по которой продукция скважины поступает в систему промыслового нефтегазосбора. Эта величина отсечет на графике (см. рис. 8.7) возможные режимы фонтанирования для условий данной скважины. Точка а соответствует минимально допустимому давлению на устье ( Ру = Рл), а ее проекция на ось абсцисс определит соответствующее этому режиму работы критическое забойное давление Ркр. Пересечение вертикали с кривой Q (рс) (точка b) дает критический дебит скважины Qкр, превышение которого приведет к давлению Ру < Рл. Таким образом, область режимов фонтанирования скважины, лежащая влево от вертикали, проходящей через точки а и b, нереальная, а область режимов, лежащая вправо от той же вертикали, осуществима, так как при условиях Рс; Q; Ру пластовая энергия превышает необходимую для подъема жидкости. Избыток энергии обусловливает устьевое давление Ру, превышающее давление в выкидной линии Рл. Для поглощения этой энергии применяется штуцер, в котором создается перепад давлений ΔРшт = Ру - Рл.