Задаемся значением объемного расходного газосодержания у входа в насос в пределах 0,15…0,25, т.к. bвсу < 0,5 и определяем по кривой 2 рисунка 1 расстояние Lн от устья скважины до сечения эксплуатационной колонны, в котором газосодержание равно принятой величине, а по кривой 1 - давление у входа в насос в стволе скважины на найденной глубине. Пусть bгвх = 0,15. Тогда Lн = 1050 м и Рвх = 5,5 МПа.
Вычисляем обводненность жидкости у входа в насос, найдя предварительно значение объемного коэффициента нефти при Рвх = 5,5 МПа:
.
Проверяем, выполняется ли неравенство (93) то есть условие бескавитационной
работы насоса. Для этого вычисляем по (93) значение (bгвх)н,
поскольку bввх<
0.5 и газожидкостная смесь относится к типу (Г+В)/Н:
сопоставляем найденное значение с bгвх = 0,277. Так как (bгвх)н > bгвх, приходим к заключению, что насос в скважине не будет кавитировать и газосепаратор перед насосом ставить нет необходимости.
. Вычисляем по (74) значение коэффициента сепарации свободного газа перед входом продукции в насос при работе его на глубине Lн = 1050 м, принимая Кс=0. Так как bввх<0,5, берем wдр.г= 0,02 м/с.
Принимаем, что для отбора заданного дебита жидкости из скважины диаметром 0,152 м надо использовать насос группы 5А. Тогда диаметр всасывающей сетки насоса, согласно таблице в п.4.7 [1, стр.28], будет Dсн=0,103 м.
Вычисляем значения приведенной скорости жидкости в зазоре между эксплуатационной колонной скважины и насосом перед всасывающей сеткой его:
Вычисляем
значение Кск:
Кс = Кск = 0,186.
.
Вычисляем по (75) [5] действительное давление насыщения жидкости в колонне НКТ,
приняв Кфн = Кфв = 1:
методом последовательной итерации находим Рд.нас=10,562 с погрешностью 10-5
. Рассчитываем методом сверху низ кривую Р(Lнкт) изменения давления вдоль колонны НКТ в интервале от устьевого сечения ее (Lнкт= 0) до глубины Lн=1050 м, найденной в п. 3.4.
. Расчет Р(Lнкт) в основном аналогичен расчету кривой Р(Lэк) и отличается от него необходимостью учета потерь давления на преодоление гидравлического трения в НКТ, то есть ведется на базе уравнения (92), но с учетом второго слагаемого в знаменателе его правой части, а также нагрева продукции, поступающей в колонну НКТ, теплом, выделяемым двигателем и насосом УЭЦН.
.
Разбиваем перепад давлений Рд.нас - Ру = 10,562 - 1,5 =
9,62 МПа на 4 ступени: DР1 = 1,9; DР2
= 2,1; DР3 =2,3; DР4
=2,5 и находим значения среднего давления для каждой ступени: Рср1=2,45;
Рср2= 4,45; Рср3= 6,65; Рср4= 9,05. Вычисляем
значения rн1 по
(12), bн1 по (11) и bвж1 по (70) для 1-го участка колонны НКТ, примыкающего к
устью скважины:
,
вычисляем средние значения объемных расходов и приведенных скоростей нефти и воды для 1-го участка НКТ:н2=0,0015*(1-0,30)*1,12 = 0,00118 м3/с; QВ2=0,3*0,0015 = 0,00045м3/сН2 =0,00114 / 0,003018 =0,377 (м/с); wВ2 =0,00045/ 0,003018 =0,0149
(м/с);
вычисляем приближенно длину первого участка колонны НКТ, соответствующего перепаду DР1, положив rн1 =757 кг/м3; rв1=1011 кг/м3, bвж1=0,276, bн1 = 1-bвж1 = 1-0,276 = 0,724; wг1 = 0, wсм1 = 0, то есть допустив, что колонна НКТ на первом участке заполнена неподвижной смесью нефти и воды с водосодержанием bвж1=0,276. Подставляем перечисленные данные в (92) и получаем:
вычисляем
расстояние от устья до середины участка DL1:
Определяем приращение температуры потока продукции за счет нагрева ее теплом двигателя и насоса по (67). Для этого предварительно оцениваем значения входящих в (67) величин Н, Сср, hд, hн, а также bвжн и rжн, используемых при вычислении Н.
Находим приближенно водосодержание в насосе по (70) при bн = bн.нас:
Вычисляем
приближенно напор насоса при работе его в скважине по (68):
вычисляем
приближенно среднюю теплоемкость жидкости в насосе по (71):
,
где Сн -средняя теплоемкость нефти, равная » 2000 Дж/(кг*К), Свсу - средняя теплоемкость пластовой воды, равная » 4380 Дж/(кг*К).
Значение hд принимают равным номинальному КПД двигателя, который должен быть спущен в скважину вместе с насосом (двигатели диаметром 117 мм, комплектуемые с насосами группы 5А, имеют hд= 0,81.
Для оценки значения КПД насоса при работе в скважине сначала определяем значение номинального КПД насоса группы 5А, номинальная подача которого не меньше (равна или несколько больше) среднего расхода продукции через насос, равного приближенно величине:жн =130*(1,103*120,0199*(1-0,27)+1*0,27) = 145,1 м3/сут.
Из справочника [2] находим ближайшую по подаче установку группы 5А - ЭЦН5А - 250 с КПД насоса 0,6. Затем находим приближенно кажущуюся вязкость продукции в насосе. Для этого сначала определяем приближенно вязкость нефти в насосе, являющейся внешней фазой проходящей через него продукции, при температуре пласта по (13):
Па*с
Но
поскольку температура продукции в насосе ниже Тпл и равна
приближенно температуре в стволе скважины перед входом в насос:
вносим поправку на вязкость нефти по номограмме Льюиса и Сквайрса [1, рисунок 4, стр.22] Вязкость нефти в насосе при Т = 318,75 К будет: mнн » 0,052 Па*с.
Так
как bвжн=0,270
<0,5, то значение кажущейся вязкости определяем по (40):
Па*с
Находим
по (73) значение параметра Вm, учитывающего влияние вязкости жидкости на КПД насоса:
Так
как Вm <
47950, КПД насоса при работе в скважине, согласно (72), будет:
Теперь по (67) находим:
Вычисляем по (65) температуру потока в НКТ на середине 1-го участка, то есть на глубине L1 = 120,3 м:
Вычисляем
значение коэффициента сверхсжимаемости попутного газа в НКТ на глубине Lнкт1=120,3
м, аналогично как в п.1.2.3 [5]:
Находим
zу по соответствующей формуле из (50) по Рпр1 и Тпр1:
Далее вычисляем значение коэффициента сжимаемости азотной части попутного газа по (62):
,
а
по (58) - значение z2:
Вычисляем объемный расход газа через среднее сечение 1-го участка колонны НКТ по (79) без слагаемого с сомножителем aг (поскольку bвсу< 0,65), положив Ккф=Кфв=1, Кс=0,186, Рвх=5,5 МПа:
Вычисляем значение приведенной скорости газа, скорости жидкости и скорости ГЖС в среднем сечении 1-го участка НКТ:
wпр.г1= 0,0011 / 0,003018 =0,368 (м/с);
wпр.ж1= 0,391 + 0,2290,149 = 0,540 (м/с);
wсм1 = 0,368 + 0,540 = 0,909 (м/с);
Вычисляем
значения 1-й и 2-й критических скоростей потока в среднем сечении 1-го участка:
Определяем
по приложению 2 тип и структуру потока нефтеводогазовой смеси через среднее
сечение 1-го участка НКТ. Так как bвж1< 0,5, wсм1 > wкр2, Рср1>
0,7МПа, смесь относится к типу (В+Г)/Н и имеет эмульсионно-пузырьковую
структуру. Вычисляем значение поверхностного натяжения между фазами ГЖС
Вычисляем по (13) значение вязкости нефти при Рср1=2,45 МПа и Тпл=357К:
Пересчитываем это значение на Т1=295,1 К по номограмме Льюиса и Сквайрса [1, рисунок 4, стр.22]. Так как снижение температуры нефти DТ1 = 357-295,1 = 61,9 К, то вязкость нефти при 295,1 К будет mн1 = 0,068 Па*с.
Вычисляем
значение параметра А по (42) и (25):
,
где
, тогда
Находим
кажущуюся вязкость жидкой части ГЖС по (41) т.к. А > 1:
Вычисляем истинное газосодержание jг1 по (36):
Вычисляем
истинную долю в жидкой части ГЖС на 1-ом участке колонны НКТ по (30) [5],
поскольку внешней фазой потока является нефть:
Находим долю нефти в жидкости по (32): jнж1 = 1 - 0,266 = 0,734.
Вычисляем истинное водосодержание по (43) и нефтесодержание по (44) в ГЖС на участке 1:
jв = jвж*(1-jг) = 0,266*(1-0,361) = 0,170;
jн = jнж*(1-jг) = 0,734*(1-0,361) = 0,469;
делаем проверку результатов оценки значений истинных долей фаз в трехфазном потоке: сумма долей фаз должна быть равна 1:
,361 + 0,170 + 0,469 = 1,000.
Вычисляем
значение плотности попутного газа при Рср1=2,45 МПа и Т1=295,1
К по (56)
Оцениваем кажущуюся вязкость ГЖС в среднем сечении 1-го участка НКТ, принимая ее равной кажущейся вязкости смеси нефти и воды в том же сечении, то есть mгжс1=mж1=0,227 (Па*с).
Вычисляем
значение числа Рейнольдса потока ГЖС по (48):
Определяем значение lсм1 по (49), поскольку полученное значение Reсм1 меньше 2000:
Вычисляем значение DL1 по (92):
12. Рассчитываем значения DL2¼DL4 колонны НКТ аналогично расчету DL1 и определяем расстояние по оси скважины от ее устья до сечения НКТ, в котором давление равно Рд.нас. Эта длина оказывается 680,5 м.
. Определяем длину участка DL5 колонны
НКТ от сечения, где давление равно Рд.нас, до глубины спуска насоса Lн = 1050 м. DL5 = 1050 -
680,5 = 369,5 м.
. Вычисляем перепад давления на длине DL5 НКТ, учитывая, что на этом участке течет водонефтяная смесь, не содержащая свободного газа, jг5 = 0, что можно принять: bн5 = bн.нас, rн5 = r н.нас и что вязкость нефти mн5 отличается от вязкости mнпл при Тпл. Расчет выполняем аналогично расчету участка DL1 эксплуатационной колонны.
Результаты расчетов кривой Р(Lнкт) представлены в приложении 2, в которой Li - расстояние по оси скважины от устья до нижнего сечения i-го участка НКТ; Рi - давление в этом сечении.
. Строим кривую Р(Lнкт) - линия 3 на рисунка 1 по значениям Рi, Li таблицы 3 и экстраполируем ее в область L > Lн = 1050 м в расчете на возможность спуска насоса в процессе дальнейшего подбора УЭЦН к скважине на глубину, большую 1050 м.
. Определяем давление в НКТ на выходе из насоса по кривой 3 рисунка 1 и давление Рс, которое требуется для работы системы скважина - УЭЦН с заданным дебитом жидкости:
Рвых = 13,43, тогда Рс = Рвых - Рвх = 13,43 - 5,5 = 7,93 МПа.
. Вычисляем среднюю температуру продукции в насосе по (64):
Вычисляем среднеинтегральный расход жидкой части продукции через насос по (88), принимая Кфн = 0,9; Кфв = 0,1; Рд.нас = 10,562 МПа:
. Вычисляем по (89) среднеинтегральный расход свободного газа через насос.
Сначала
находим значения А,В и zcр в
насосе:
,
Значение
zcр
определяем по (58) при Тср.н= 299,07 и Рвх = 5,5 МПа и rу.отн =
0,996.
Находим
zу по
соответствующей формуле из (50) по Рпр1 и Тпр1:
Далее вычисляем значение коэффициента сжимаемости азотной части попутного газа по (62):
а
по (58) - значение z2:
Подставив
значения А, В и zср, в (89)
получаем:
. Вычисляем среднеинтегральный расход ГЖС через насос по (85):
Qср = 0,00188 + 0,000028 = 0,0019 (м3/с).
.
Вычисляем массовый расход через насос по (76):
.
Вычисляем среднеинтегральную плотность продукции в насосе по (90):
. Вычисляем напор, который необходим для работы системы скважина - УЭЦН с заданным дебитом Qжсу = 0,0015 м3/с по (91):
.
Вычисляем среднеинтегральное газосодержание в насосе:
. Определяем кажущуюся вязкость жидкости и ГЖС в насосе при Тср.н=299,07 К.
Поскольку вязкость нефти, являющейся внешней фазой продукции в насосе, Тпл=357 К равна 0,0329 Па*с (согласно п.1.9.1), то при Тср.н=299,07 К, пользуясь графиком Льюиса и Сквайрса, находим mнн » 0,067 Па*с. Кажущаяся вязкость жидкой части так же как и ГЖС в насосе, будет: