NGPO
.pdfПостоим паспортную характеристику насоса. Паспортная характеристика |
||||||||||||
насоса |
построена |
для |
эталонной |
жидкости |
(плотность |
1000 |
кг/м3, |
вязкость |
||||
1мПа∙с). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Паспортная характеристика УЭЦН |
|
|
|
||||||
|
|
3000 |
|
|
|
|
|
|
|
|
|
|
|
|
2500 |
|
|
|
|
|
|
|
|
|
|
|
м |
2000 |
|
|
|
|
|
|
|
|
|
|
|
Н, |
|
|
|
|
|
|
|
|
|
|
|
|
Напор |
1500 |
|
|
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
500 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
10 |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
90 |
|
|
|
|
|
|
|
Дебит Q, м3/сут |
|
|
|
|
|
3. Пересчет паспортной характеристики насоса в реальную с учетом
влияния газа, вязкости, наличия твердых взвешанных частиц и
конструктивных особенностей ЭЦН
Пересчитываем полученную характеристику для продукции скважины.
Пересчет будем вести по методике изложенной в [4, с. 683].
Откладываем Qоптв=60 м3/сут на шкале подач номограммы для определения параметров центробежных насосов на вязких жидкостях - графики П. Д. Ляпкова.
Из полученной точки проводим вертикальную линию вверх до пересечения с соответствующей линией числа оборотов n=2890 об/мин. Из точки пересечения линий Qопт в и n проводим горизонтальную линию до пересечения с линией,
соответствующей средней кинематической вязкости жидкости в насосе,
определяемой по формуле:
н пл
Из точки пересечения этих линий проводим вертикальную прямую до пересечения ее с соответствующими кривыми пересчетных коэффициентов: КQ и
КН. Получаем КQ =0,97 и КН=0,96.
11
Пересчитываем Qопт ж для продукции скважины:
опт К |
|
м |
|
опт |
|
||
сут |
|||
|
|
К
м
сут
К
м
сут
Пересчитаем оптимальный напор Ноптж для продукции скважины:
опт Кн |
Нопт |
м |
Кн |
Н |
м |
Кн |
Н |
м |
4.Проведение корректировки режима работы УЭЦН и скважины
4.1Определение диаметра штуцера. Построение зависимости потерь
давления на штуцере от его диаметра
Расчет штуцера состоит из трех областей: потери при внезапном сужении канала, при внезапном расширении канала и гидравлические потери в штуцере.
При внезапном сужении канала поток жидкости отрывается от стенок входного участка и затем касается стенок канала меньшего сечения. Образуются две зоны интенсивного вихреобразования. Коэффициент потерь напора при гидравлическом сопротивлении внезапного сужения потока определяется по следующей формуле:
где S1, S2 – площади внутренних поперечных сечений выкидной линии до штуцера и штуцера.
Потери на трение из-за внезапного сужения трубы:
где v1, v2 – скорости движения потока до штуцера и в штуцере.
12
Выразив одно через другое получим:
При внезапном расширении канала образуется зона, в которой жидкость практически не участвует в движении по трубам, образуя локальный вихревой поток, где претерпевает деформацию. Тогда коэффициент потерь напора определяется по следующей формуле:
где S3 – площадь внутреннего поперечного сечения трубы после штуцера.
Потери жидкости на трение во время движения по внутренней полости штуцера:
где l – длина штуцера, м;
d2 – внутренний диаметр штуцера, м.
Тогда суммарные потери напора на штуцере:
Учитывая что d1=d3 перепишем:
Приняв, что в устьевой арматуре УЭЦН на выкидных линиях используются трубы 89×6 мм, получим:
d1=89 – 2∙6 = 77мм.
Число Рейнольдса определим по формуле:
где - линейная скорость потока ГЖС в штуцере, м/с, определяется:
13
где - вязкость жидкости, Паּс;
ρгжс - плотность ГЖС при Ратм Дебит скважины: QСКВ = 103 м3/сут;
Режим течения жидкости - турбулентный, следовательно:
Потери при d1=3 мм:
Тогда потери давления на штуцере:
Аналогично рассчитываем потери и для других диаметров штуцера.
Результаты расчетов представлены в таблице 4.1.
|
|
|
|
|
Таблица 4.1 |
|
|
|
|
|
|
|
|
Внутренний |
ν, м/с |
Re |
λ |
Потери |
Потери |
|
диаметр |
|
|
|
напора на |
давлении на |
|
штуцера |
|
|
|
штуцере ∆h, |
штуцере |
|
dшт, мм |
|
|
|
м |
∆P, МПа |
|
3 |
168,4 |
289387,4 |
0,014 |
1403,5 |
9,622 |
|
6 |
42,1 |
144693,7 |
0,016 |
79,5 |
0,545 |
|
7 |
30,9 |
123900,4 |
0,017 |
41,6 |
0,285 |
|
8 |
23,7 |
108606,2 |
0,017 |
23,7 |
0,162 |
|
10 |
15,2 |
87068,3 |
0,018 |
9,2 |
0,063 |
|
14
|
|
|
Потери давления в штуцере |
|
|
|||
, МПа |
12 |
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
в штуцере |
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
давления |
6 |
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Перепад |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
0 |
2 |
4 |
6 |
8 |
10 |
12 |
14 |
|
|
|
|
Диаметр штуцера, мм |
|
|
|
Рис. 4.1 Потери давления на штуцере Из рисунка видно, что наибольшие потери давления достигаются при
минимальном диаметре штуцера, поэтому при необходимости возможна установка штуцеров больших диаметров.
4.2 Регулирование глубины подвески насоса. Построение кривой
распределения давления, температуры и газосодержания по стволу
скважины
Используя метод Ф. Поэтмана – П. Карпентера. Расчет ведем «сверху-вниз».
1. Задаем шаг а, нас а и определяем число расчетных точек:
2. Рассчитываем температурный градиент потока
ж ст |
Т |
м
15
где Тпл м - средний геотермический градиент скважины, Qж ст – дебит скважины по жидкости при стандартных условиях; DТ – внутренний диаметр колонны НКТ, м.
3. Определяем температуру ГЖС на устье скважины
Ту Тс пН
Рассчитываем температуру потока в рассматриваемых сечениях (точках)
потока. Например, в сечении, где , температура будет:
5. Используя данные исследования проб пластовой нефти, определяем
физические параметры, соответствующие заданным давлениям
Г н
6. Вычислим коэффициент сверхсжимаемости газа Z, для этого определим приведенные параметры смеси газов:
гу – относительная по воздуху плотность смеси газов,
пр
пр
Коэффициент сверхсжимаемости газовой смеси при Р = 2,25 МПа: При и
7. Определяем удельный объем ГЖС при Р=2,25 |
а |
|||||
|
|
в |
|
|
в |
|
с |
н |
|
|
|
в |
|
|
|
|
|
|
16
см |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
м м |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
8. |
Определяем удельную массу смеси при стандартных условиях |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
во |
в |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с |
|
н |
|
|
о |
|
|
в |
в |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
см |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
кг м |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
9. |
Рассчитываем идеальную плотность ГЖСм при Р=2,25 а |
|||||||||||||||||
|
|
|
|
см |
|
|
|
|
|
|
|
|
кг м |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
см |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
см
10. Определяем корреляционный коэффициент :
всм
жнкт
11. Вычисляем полный градиент давления при Р=2,25 |
|
а |
|||||||
|
|
|
|
|
|
ст |
в |
с |
|
|
|
|
|
с и |
|
|
|
|
|
|
|
|
|
|
с |
и н т |
|||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
Вычисляем - обратные расчетным градиентам давления. 12. Распределение давления (см. табл.)
Нi=(Ркон-Ру)/N*(((dH/dP)заб+(dH/dP)i)/2+( dH/dP)1+( dH/dP)2+……+( dH/dP)i-1)
Проводим численное интегрирование зависимости dH/dp = f(p), в результате чего получаем распределение давления на участке НКТ, где происходит течение газожидкостного потока.
Расчет ас е е ения ав ения в э с атационной о онне
Аналогично рассчитываем давление по методу Поэтмана-Карпентера, но
«снизу вверх», при этом в качестве берем , и учитываем, что расчет ведем в кольцевом зазоре ( экв).
экв |
эк нкт |
м |
17
Результаты расчета распределения давления в колонне НКТ
|
|
Vгр, |
Vгв, |
|
|
|
|
Vсм, |
Mсм, |
ρсм, |
|
dP/dH, |
dH/dP, |
|
P, МПа |
Т, К |
м3/м3 |
м3/м3 |
bн |
Pпр |
Tпр |
Z |
м3/м3 |
кг/м3 |
кг/м3 |
f |
МПа/м |
м/МПа |
H, м |
1,25 |
277,9 |
1,06 |
86,54 |
0,94 |
0,270 |
1,81 |
0,988 |
7,94 |
823 |
104 |
0,039 |
0,0022 |
454,4 |
0 |
2,25 |
279,2 |
3,32 |
84,28 |
1,01 |
0,487 |
1,82 |
0,980 |
4,79 |
823 |
172 |
0,039 |
0,0024 |
417,0 |
428 |
3,25 |
280,5 |
17,87 |
69,73 |
1,05 |
0,703 |
1,82 |
0,972 |
3,23 |
823 |
255 |
0,039 |
0,0030 |
335,3 |
798 |
4,25 |
281,8 |
28,91 |
58,69 |
1,08 |
0,920 |
1,83 |
0,966 |
2,49 |
823 |
330 |
0,039 |
0,0036 |
276,9 |
1099 |
5,25 |
283,1 |
37,60 |
50,00 |
1,10 |
1,136 |
1,84 |
0,960 |
2,09 |
823 |
394 |
0,039 |
0,0042 |
239,3 |
1353 |
6,25 |
284,4 |
44,78 |
42,82 |
1,12 |
1,352 |
1,85 |
0,955 |
1,84 |
823 |
447 |
0,039 |
0,0047 |
214,6 |
1576 |
7,25 |
285,7 |
50,89 |
36,71 |
1,14 |
1,569 |
1,86 |
0,951 |
1,68 |
823 |
490 |
0,039 |
0,0051 |
197,7 |
1779 |
8,25 |
287,0 |
56,21 |
31,39 |
1,15 |
1,785 |
1,87 |
0,948 |
1,57 |
823 |
524 |
0,039 |
0,0054 |
185,9 |
1968 |
9,25 |
288,2 |
60,92 |
26,68 |
1,17 |
2,001 |
1,87 |
0,945 |
1,49 |
823 |
552 |
0,039 |
0,0058 |
171,8 |
2143 |
10,25 |
289,5 |
65,14 |
22,46 |
1,18 |
2,218 |
1,88 |
0,942 |
1,43 |
823 |
574 |
0,039 |
0,0061 |
165,0 |
2309 |
11,25 |
290,8 |
68,97 |
18,63 |
1,19 |
2,434 |
1,89 |
0,941 |
1,39 |
823 |
591 |
0,039 |
0,0063 |
159,9 |
2469 |
12,25 |
292,1 |
72,48 |
15,12 |
1,20 |
2,650 |
1,90 |
0,939 |
1,36 |
823 |
605 |
0,039 |
0,0064 |
156,1 |
2624 |
13,25 |
293,4 |
75,71 |
11,89 |
1,21 |
2,867 |
1,91 |
0,938 |
1,33 |
823 |
617 |
0,039 |
0,007 |
153,1 |
2776 |
14,25 |
294,7 |
78,70 |
8,90 |
1,22 |
3,083 |
1,92 |
0,938 |
1,31 |
823 |
626 |
0,039 |
0,0066 |
150,8 |
2926 |
15,25 |
296,0 |
81,49 |
6,11 |
1,22 |
3,299 |
1,92 |
0,937 |
1,30 |
823 |
633 |
0,039 |
0,0067 |
148,9 |
3073 |
16,25 |
297,3 |
84,11 |
0 |
1,23 |
3,516 |
1,93 |
0,938 |
1,27 |
823 |
649 |
0,039 |
0,0069 |
145,0 |
3217 |
17,25 |
298,6 |
86,56 |
0 |
1,24 |
3,732 |
1,94 |
0,938 |
1,27 |
823 |
646 |
0,039 |
0,0069 |
145,6 |
3360 |
18,25 |
299,9 |
88,88 |
0 |
1,24 |
3,949 |
1,95 |
0,939 |
1,28 |
823 |
643 |
0,039 |
0,0068 |
146,3 |
3504 |
18
Результаты расчёта распределения давления в эксплуатационной колонне
|
|
Vгр, |
Vгв, |
|
|
|
|
Vсм, |
Mсм, |
ρсм, |
|
dP/dH, |
dH/dP, |
|
P, МПа |
Т, К |
м3/м3 |
м3/м3 |
bн |
Pпр |
Tпр |
Z |
м3/м3 |
кг/м3 |
кг/м3 |
f |
МПа/м |
м/МПа |
H, м |
1,25 |
277,9 |
1,06 |
86,54 |
0,94 |
0,270 |
1,807 |
0,988 |
7,94 |
823 |
104 |
0,118 |
0,0046 |
216,7 |
1383 |
2,25 |
279,2 |
3,32 |
84,28 |
1,01 |
0,487 |
1,815 |
0,980 |
4,79 |
823 |
172 |
0,118 |
0,0039 |
259,3 |
1621 |
3,25 |
280,5 |
17,87 |
69,73 |
1,05 |
0,703 |
1,824 |
0,972 |
3,23 |
823 |
255 |
0,118 |
0,0040 |
252,3 |
1877 |
4,25 |
281,8 |
28,91 |
58,69 |
1,08 |
0,920 |
1,832 |
0,966 |
2,49 |
823 |
330 |
0,118 |
0,0044 |
228,9 |
2118 |
5,25 |
283,1 |
37,60 |
50,00 |
1,10 |
1,136 |
1,841 |
0,960 |
2,09 |
823 |
394 |
0,118 |
0,0048 |
207,8 |
2336 |
6,25 |
284,4 |
44,78 |
42,82 |
1,12 |
1,352 |
1,849 |
0,955 |
1,84 |
823 |
447 |
0,118 |
0,0052 |
191,6 |
2536 |
7,25 |
285,7 |
50,89 |
36,71 |
1,14 |
1,569 |
1,857 |
0,951 |
1,68 |
823 |
490 |
0,118 |
0,0056 |
179,6 |
2721 |
8,25 |
287,0 |
56,21 |
31,39 |
1,15 |
1,785 |
1,866 |
0,948 |
1,57 |
823 |
524 |
0,118 |
0,0059 |
170,8 |
2896 |
9,25 |
288,2 |
60,92 |
26,68 |
1,17 |
2,001 |
1,874 |
0,945 |
1,49 |
823 |
552 |
0,118 |
0,0061 |
164,2 |
3064 |
10,25 |
289,5 |
65,14 |
22,46 |
1,18 |
2,218 |
1,883 |
0,942 |
1,43 |
823 |
574 |
0,118 |
0,0063 |
159,2 |
3226 |
11,25 |
290,8 |
68,97 |
18,63 |
1,19 |
2,434 |
1,891 |
0,941 |
1,39 |
823 |
591 |
0,118 |
0,0076 |
130,9 |
3371 |
12,25 |
292,1 |
72,48 |
15,12 |
1,20 |
2,650 |
1,900 |
0,939 |
1,36 |
823 |
605 |
0,118 |
0,0078 |
127,6 |
3500 |
19
Строим график распределения давления:
Распределение давлений
Н, м
Р, МПа
0 |
2 |
4 |
6 |
8 |
10 |
12 |
14 |
16 |
18 |
0
500
1000
1500
2000
2500
3000
3500
4000
НКТ ЭК
Газонасыщенность нефти
G, м3/т
140,0 |
|
|
|
|
120,0 |
|
|
|
|
100,0 |
|
|
|
|
80,0 |
|
|
|
|
60,0 |
|
|
|
|
40,0 |
|
|
|
|
20,0 |
|
|
|
|
0,0 |
|
|
|
|
0,00 |
5,00 |
10,00 |
15,00 |
20,00 |
Р, МПа
20