Расчет промывки после ГРП

Расчет технологических параметров при проведении ГРП

Расчет проведения ГРП.

Расчет гидравлического разрыва пласта представляет собой достаточно сложную задачу, которая состоит из двух частей: основных характеристик процесса и выбор необходимой техники для его осуществления; определение вида трещины и расчет ее размеров.

Исходные данные:

Эксплуатационная скважина №17 Фроловского месторождения.

Таблица 2.7.1

Показатель Обозначение Величина Размерность
Глубина скважины L м
Диаметр по долоту D 0,168 м
Вскрытая толщина пласта H 22,7 м
Средняя проницаемость K 2,7*10-10 м2
Модуль упругости пород E 2,0*1010 Па
Коэффициент Пуассона v 0,25  
Средняя плотность пород над продуктивным горизонтом rп 2385,2 кг/м3
Плотность жидкости разрыва rн кг/м3
Вязкость жидкости разрыва m 0,2 Па*с
Концентрация песка C кг/м3
Темп закачки Q 1,2*10-2 м3

 

1. Рассчитываем вертикальную составляющую горного давления:

РГВ = rn . g . L . 10-6 = 2385,2 . 9,81 . 2887 . 10-6 = 67,55 МПа

2. Рассчитываем горизонтальную составляющую горного давления:

v 0,25

РГГ = РГВ . ___________ = 67,55 . _______________ = 22,52 МПа

( 1 – v ) 1- 0,25

 

В данных условиях при ГРП следует ожидать образование вертикальной трещины.

Запроектируем гидроразрыв не фильтрующейся жидкостью. В качестве жидкости разрыва и жидкости песконосителя используем загущенную нефть с добавкой асфальтина, плотность и вязкость в таблице исходных данных. Содержание песка по данным 300 кг / м3, для расклинивания трещины планируем закачку примерно 8 тонн кварцевого песка фракции 0,8 – 1,2 мм, темп закачки по данным 1,2 * 10 –2, что больше минимально допустимого при создании вертикальных трещин.

При ГРП непрерывно закачивают жидкость-песконоситель в объеме 10,08 м3 , которая одновременно является и жидкостью разрыва.

3.

5,25 . . (1 – v2)2
РЗАБ РГГ  
Рассчитаем забойное давление разрыва пласта:

.
=
Е Р ГГ
=
- 1
.
 
РЗАБ 3 3 2 2 Q .m

РГГ Р ГГ

0,000502
=
=
5,25 . ( 2 . 1010 )2 . 0,012 . 0,2 0,0504 . 10 20

(1 – 0,252)2 . (22,52 . 106)2 . 22,52 . 106 10038,012 . 1018

 
 
  РЗАБ = 24,3 МПа


= 1,079
РЗАБ

22,52

4. При закачке жидкости –песконосителя давление на устье скважины

Р У = РЗАБ - r жп . g . L + РТР ,

Где r жп - плотность жидкости - песконосителя, кг/м3

r жп = r н . (1 - b П) + r . b П ,

где r н - плотность жидкости разрыва, кг/м3

r - плотность песка, кг/м3 (2500 кг/м3)

b П - объемная концентрация песка в смеси

С / r 300 / 2500

b П = = = 0,107

С / r + 1 300/2500 + 1

rжп = 930 . (1-0,107) + 2500 . 0,107 = 1098 кг / м3

 

Потери давления на трение жидкости-песконосителя:

8 . l . Q 2 . L . r П

РТР = ,

p 2 . d 5ВН

где l - коэффициент гидравлических сопротивлений

64 4 . Q . r жп

l = Re = ,

Re p 2 . d вн . m жп

m жп вязкость жидкости с песком,Па*с

m жп = m . exp . ( 3.18 . bП )

m жп = 0,2. exp . (3,18 . 0,107) = 0,281

4 . 0,012 . 1098 52,704

Re = = = 964

3,14 . 0,062 . 0.281 0,0547

 

Коэффициент гидравлического сопротивления

l = = 0,066

Тогда,

8 . 0,066 . 0,012 2 . 2887 . 1098 241,0158

РТР = = = 27 МПа

3,14 2 . 0,062 5 0,000009

 

Учитывая, что Re = 964 >200, потери на трение составят:

Р ТР = 1,52 . РТР = 1,52 . 27 = 41,04 МПа

Давление на устье скважины при закачке жидкости-песконосителя:

РY = 24,3 – 1098 . 9,81 . 2887 . 10-6 + 41,04 = 34,24 МПа

5. Рабочие жидкости в скважину закачивают насосными агрегатами 4АН-700.

Таблица 2.7.1.1

Скорость Подача, м3 Давление, МПа
0,0063 71,9
0,0085 52,9
0,012 37,4
0,015 29,8

 

При работе агрегата 4АН-700 на III скорости РР = 37,4 Мпа, а Q p = 0,012 м3/с.

Необходимое число насосных агрегатов:

PY . Q 34,2 . 0,012

N = + 1 = + 1 = 2,83

РР . QP . K TC 37,4 . 0,012 . 0,5

где РР - рабочее давление агрегата;

QP – подача агрегата при данном РР;

K TC - коэффициент технического состояния агрегата K TC = 0,5-0,8

6. Необходимый объем продавочной жидкости (при закачке в НКТ):

Vп = 0,785 . d2ВН . L = 0,785 . 0,062 2 . 2887 . 1,2 = 10,08 м3

7. Объем жидкости для осуществления гидроразрыва (жидкость разрыва и жидкость-песконоситель):

QП 8000

Vж = = = 26,7 м 3

С 300

8. Суммарное время работы одного агрегата 4АН-700 на III скорости:

Vж + Vп 26,7 + 10,08

t = = = 3065 c. или 51 мин.

0,012

 

Расчет размеров трещин:

1. Вычисляем длину вертикальной трещины.

V ж . Е 5,6 . (1-v)2 . Н . ЗАБ - РГГ)  

l = =

 

       
   


= = = 64,8 м.

 

2. Рассчитываем раскрытость трещины:

4 . (1-v)2 . 1 . (РЗАБ – РГГ)

wo = =

Е

4 . (1-0,25)2 . 64,8 . (24,3-22,52) 259,524

= = = 0,013 м или 1,3 см.

2,0 . 10 4 20000

 

Таким образом, в результате проведения гидроразрыва в данной скважине образуется вертикальная трещина длиной 64,8 метров и шириной на стенке скважины 1,3 см.

 

Расчет промывки после ГРП.

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

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

Исходные данные:

Эксплуатационная скважина №17 Фроловского месторождения.

Таблица 2.7.2

Показатель Обозначение Величина Размерность
Глубина скважины Н м
Диаметр эксплуатационной колонны DК 168 х 8,9 мм
Диаметр промывочных труб dH 73 х 5,0 мм
Максимальный размер зерен а 1,2 мм
Кинематическая вязкость нефти n 0,4 * 10-4 м2
Плотность нефти rн кг/м3
Пористость песчаной пробки m 0,3  
Высота песчаной пробки l м
Плотность песчаной пробки rп кг/м3
Плотность промывочной жидкости rж кг/м3
Длина нагнетательной линии L м

 

Промывка ведется промывочным агрегатом АзИНМаш – 32.

Прямая промывка водой.

Гидравлические потери давления при движении жидкости внутри промывочных труб находятся по формуле (4) для всех скоростей агрегата.

, (Па) (4)

где: l - коэффициент гидравлического сопротивления (табл. 2.7.2.5);

H - глубина скважины, м ;

dв - внутренний диаметр промывочных труб, м;

vH - скорость нисходящего потока жидкости, м/с (табл. 2.7.2.1);

rж – плотность жидкости, кг/м3 .

На первой скорости:

2887 1,1852 . 1000

p1I = 0,035 . . = 1,13 МПа;

0,063 2

На второй скорости:

2887 1,8402 . 1000

p1II = 0,035 . . = 2,7 МПа;

0,063 2

На третьей скорости:

2887 3,032 . 1000

p1III = 0,035 . . = 7,4 МПа;

0,063 2

На четвертой скорости:

2887 4,72 . 1000

p1IV = 0,035 . . = 17,7 МПа.

0,063 2

Аналогично по формуле (5), находятся гидравлические потери давления при движении смеси жидкости с песком в кольцевом пространстве скважины:

, (Па) (5)

где φ = 1,1-1,2 -коэффициент, учитывающий повышение гидравлического сопротивления от содержания песка в жидкости;

λ - коэффициент гидравлического сопротивления при движении воды в кольцевом пространстве (определяется по диаметру труб, эквивалентному разности диаметров D и dн);

D – наружный диаметр эксплуатационной колонны, м;

dн - наружный диаметр промывочных труб, м;

υв - скорость восходящего потока жидкости в кольцевом пространстве, м/с (табл. 2.7.2.2) .

На первой скорости:

 

2887 0,2652 . 1000

p2I = 0,034 . 1,2 . . = 0,05457 МПа;

0,1502 – 0,073 2

На второй скорости:

 

2887 0,4122 . 1000

p2II = 0,034 . 1,2 . . = 0,1295 МПа;

0,1502 – 0,073 2

На третьей скорости:

 

2887 0,6872 . 1000

p2III = 0,034 . 1,2 . . = 0,360 МПа;

0,1502 – 0,073 2

На четвертой скорости:

 

2887 1.0522 . 1000

p2IV = 0,034 . 1,2 . . = 0,8443 МПа.

0,1502 – 0,073 2

Потери давления для уравновешивания разности удельных весов жидкости в трубах и в затрубном пространстве находится по формуле К. А. Апресова (6).

, (6)

где m – пористость песчаной пробки (m = 0,3÷0,45);

F-площадь сечения обсадной колонны, м2;

l - высота промытой песчаной пробки по длине одной трубы или одного колена труб, м;

f - площадь сечения кольцевого пространства при прямой промывке и площадь сечения промывочных труб при обратной промывке, м2;

ρн - плотность песка, кг/м3 (для кварцевого песка ρн = 2650÷2700кг/м3);

ρж - плотность промывочной жидкости, кг/м3;

ω - скорость свободного падения частиц песка и жидкости, м/с (табл. 2.7.2.3);

vв - скорость восходящего потока, м/с.

 

 

На первой скорости:

 

р3I = = 0,0559 МПа;

 

На второй скорости:

 

р3II = = 0,099 МПа;

 

На третьей скорости:

 

р3III = = 0,1265 МПа;

На четвертой скорости:

 

р3IV = = 0,143 МПа.

Гидравлические потери давления в шланге и вертлюге находятся по опытным данным (табл. 2.7.2.4).

На первой скорости: р4 I = 0,0572 МПа;

На второй скорости: р4 II = 0,145 МПа;

На третьей скорости: р4 III = 0,407 МПа;

На четвертой скорости: р4 IV = 0,844 МПа.

Гидравлическое сопротивление в нагнетательной линии p5 от насоса до шланга определяется аналогично сопротивлению в промывочных трубах диаметр нагнетательной принимается равным диаметру промывочных труб. (формула 4).

, (Па) (4)

На первой скорости:

40 1,185 2 . 1000

р5 I = 0,035 . . = 0,0156 МПа;

0,063 2

 

На второй скорости:

40 1,840 2 . 1000

р5 II = 0,035 . . = 0,0376 МПа;

0,063 2

На третьей скорости:

40 3,03 2 . 1000

р5 III = 0,035 . . = 0,1020 МПа;

0,063 2

На четвертой скорости:

40 4.7 2 . 1000

р5 IV = 0,035 . . = 0,2454 МПа

0,063 2

 

Давление на выкиде насоса на различных скоростях находится по формуле (7).

Рн=pобщ=p1 + p2+ p3+p4+p5. (7)

РнI = 1,13 + 0,05457 + 0,0559 + 0,0572 +0,0156 = 1,31 Мпа;

РнII = 2,7 + 0,1295 + 0,099 + 0,145 +0,0376 = 3,11 МПа;

РнIII = 7,4 + 0,360 + 0,1265 + 0,407 + 0,1020 = 8,39 МПа;

РнIV = 17,7 + 0,8443 + 0,143 + 0,844 + 0,2454 = 19,78 МПа.

Давление на забой скважины определяется по формуле (8):

Рзаб= . (8)

 

2887 . 10 3 . 9,81

РзабI = + 0,05457 + 0,0559 = 28,43 МПа;

10 6

2887 . 10 3 . 9,81

РзабII = + 0,1295 + 0,099 = 28,55 МПа;

10 6

2887 . 10 3 . 9,81

РзабIII = + 0,360 + 0,1265 = 28,81 МПа;

10 6

2887 . 10 3 . 9,81

РзабIV = + 0,8443 + 0,143 = 29,31 МПа.

10 6

Мощность, необходимая для промывки песчаной пробки вычислим по формуле (9):

N =Pн Q/ ηа, (9)

где Q - производительность насоса в л/с;

ηа- общий механический к. п. д. промывочного агрегата.

 

1,31 . 106 . 3,58 . 10 -3

NI = = 6,89 кВт.;

10 3 . 0,68

3,11 . 106 . 5,56 . 10 -3

NII = = 25,43 кВт.;

10 3 . 0,68

 

8,39 . 106 . 9,15 . 10 -3

NIII = = 112,89 кВт.;

10 3 . 0,68

19,78 . 106 . 14,2 . 10 -3

NIV = = 413,05 кВт.

10 3 . 0,68

 

Агрегат АзИНМаш – 32 имеет максимальную мощность 110 кВт., поэтому работа на третьей и четвертой скоростях невозможна.

Коэффициент использования максимальной мощности промывочного агрегата (%) рассчитывается по формуле (10)

Kа= . (10)

6,89 . 100

КаI = = 6,3 %;

25,43 . 100

КаII = = 23,1 %.

 

Скорость подъема размытого песка определим оп формуле (1):

υп = υв – ω, (1)

где υп - скорость подъема песчинок

υв - скорость восходящего потока жидкости;

ω - скорость падения наиболее крупных частиц в жидкости.

vПI = 0,265 – 0,1102 = 0,155 м/с;

vПII = 0,412 – 0,1102 = 0,302 м/с.

Время t необходимое для подъема размытой песчаной пробки на поверхность с глубины Н определим по формуле (2):

t=H/ υп (2)

tI = = 310 мин.;

0,155 . 60

tII = = 159 мин.

0,302 . 60

Обратная промывка водой.

Гидравлические потери давления при движении смеси жидкости с песком в кольцевом пространстве скважины определим по формуле (11):

. (Па) (11)

 

2887 0,265 2 . 1000

p1I = 0,034 . . = 0,0446 МПа;

0,1502 – 0,073 2

2887 0,412 2 . 1000

p1II = 0,034 . . = 0,107 МПа;

0,1502 – 0,073 2

2887 0,687 2 . 1000

p1III = 0,034 . . = 0,300 МПа;

0,1502 – 0,073 2

2887 1,052 2 . 1000

p1IV = 0,034 . . = 0,703 МПа.

0,1502 – 0,073 2

 

Гидравлические потери давления при движении смеси жидкости с песком внутри промывочных труб определим по формуле (12):

, (Па) (12)

 

2887 1,185 2 . 1000

p2I = 0,035 . 1,2 . . = 1,3513 МПа;

0,063 2

 

2887 1,840 2 . 1000

p2II = 0,035 . 1,2 . . = 3,2580 МПа;

0,063 2

 

2887 3,03 2 . 1000

p2III = 0,035 . 1,2 . . = 8,8351 МПа;

0,063 2

 

2887 4,7 2 . 1000

p2IV = 0,035 . 1,2 . . = 21,2579 МПа.

0,063 2

Потери давления для уравновешивания разности удельных весов жидкости в трубах и в затрубном пространстве найдем по формуле К.А. Апресова (6):

р3I = = 0,2472 МПа;

 

р3II = = 0,4310 МПа;

р3III = = 0,5637 МПа;

 

р3IV = = 0,6327 МПа.

Гидравлическими потерями в шланге и вертлюге при обратной промывке можно пренебрегать.

Потери давления в нагнетательной линии будут такие же, что и при прямой промывке.

p5I = 0,0156 МПа;

p5II = 0,0376 МПа;

p5III = 0,1020 МПа;

p5IV = 0,2454 МПа.

Давление на выкиде насоса на различных скоростях определяем по формуле (7):

РнI = 0,0446 + 1,3513 + 0,2472 + 0,0156 = 1,66 МПа;

РнII = 0,107 + 3,2580 + 0,4310 + 0,0376 = 3,83 МПа;

РнIII = 0,300 + 8,8351 + 0,5637 + 0,1020 = 9,80 МПа;

РнIV = 0,703 + 21,2579 + 0,6327 + 0,2454 = 22,84 МПа.

Давление на забой скважины определяется по формуле (8):

Рзаб= . (8)

 

2887 . 10 3 . 9,81

РзабI = + 1,3513 + 0,2472 = 29,92 МПа;

10 6

2887 . 10 3 . 9,81

РзабII = + 3,2580 + 0,4310 = 32,01 МПа;

10 6

2887 . 10 3 . 9,81

РзабIII = + 8,8351 + 0,5637 = 37,72 МПа;

10 6

2887 . 10 3 . 9,81

РзабIV = + 21,2579 + 0,6327 = 50,21 МПа.

10 6

Мощность, необходимую для промывки песчаной пробки вычислим по формуле (9):

1,66 . 106 . 3,58 . 10 -3

NI = = 8,74 кВт.;

10 3 . 0,68

 

3,83 . 106 . 5,56 . 10 -3

NII = = 31,32 кВт.;

10 3 . 0,68

 

9,80 . 106 . 9,15 . 10 -3

NIII = = 131,87 кВт.;

10 3 . 0,68

 

22,84 . 106 . 14,2 . 10 -3

NIV = = 476,95 кВт.

10 3 . 0,68

 

Работа на третьей и четвертой скоростях невозможна.

Коэффициент использования максимальной мощности промывочного агрегата (%) рассчитаем по формуле (10):

8,74 . 100

КаI = = 7,9 %;

 

31,32 . 100

КаII = = 28,5 %.

Скорость подъема размытого песка рассчитаем оп формуле (1):

vПI = 1,185 – 0,1102 = 1,07 м/с;

vПII = 1,840 – 0,1102 = 1,73 м/с.

Время t необходимое для подъема размытой песчаной пробки на поверхность с глубины Н определим по формуле (2):

tI = = 45 мин.;

1,07 . 60

 

tII = = 28 мин.

1,73 . 60

 

Справочные таблицы.

 

Скорость движения жидкости в промывочных трубах в м/с.

Таблица 2.7.2.1

Расход жидкости, л/с Наружный диаметр промывочных труб, мм
0,495 0,331 0,22 0,12
0,999 0,662 0,44 0,252
1,485 0,993 0,66 0,378
1,985 1,324 0,88 0,504
2,475 1,655 1,1 0,66
2,97 1,986 1,32 0,756
3,465 2,317 1,54 0,882
3,96 2,648 1,76 1,008
4,95 3,31 2,2 1,26
7,426 4,966 3,3 1,89

 

 

Скорость движения жидкости в кольцевом пространстве υв в м/с.

Таблица 2.7.2.2

Расход жидкости, л/с Наружный диаметр эксплуатационной колонны, мм.
Наружный диаметр промывочных труб, мм.
 
0,12 0,27 0,10 0,12 0,07 0,07 0,09 0,05 0,05 0,06 0,07 0,03 0,04 0,04 0,05
0,39 0,53 0,20 0,23 0,14 0,15 0,18 0,09 0,10 0,11 0,15 0,07 0,07 0,08 0,09
0,59 0,8 0,30 0,35 0,2 0,22 0,26 0,14 0,15 0,17 0,22 0,10 0,10 0,11 0,14
0,79 1,06 0,40 0,46 0,27 0,3 0,35 0,2 0,20 0,22 0,29 0,14 0,14 0,15 0,18
0,98 1,33 0,50 0,58 0,34 0,37 0,44 0,24 0,25 0,28 0,37 0,17 0,18 0,19 0,23
1,18 1,6 0,60 0,69 0,41 0,45 0,52 0,28 0,30 0,34 0,44 0,20 0,21 0,23 0,27
1,38 1,86 0,70 0,81 0,47 0,52 0,61 0,33 0,35 0,39 0,51 0,24 0,25 0,27 0,32
1,58 2,13 0,80 0,92 0,54 0,59 0,7 0,38 0,40 0,45 0,58 0,27 0,28 0,30 0,36
1,97 2,66 1,0 1,15 0,68 0,74 0,87 0,47 0,50 0,56 0,73 0,34 0,35 0,38 0,45
2,95 3,99 1,50 1,73 1,01 1,11 1,31 0,71 0,75 0,84 1,1 0,51 0,53 0,57 0,68
                                 

 

Критическая скорость падения частичек ω.

Таблица 2.7.2.3

Максимальный размер зерен, мм Скорость свободного падения, м/с Максимальный размер зерен, мм Скорость свободного падения, м/с · 104 Максимальный размер зерен, мм Скорость свободного падения, м/с · 104
0,01 0,0001 0,23 0,028 0,095
0,03 0,0007 0,25 0,03 1,2 0,1102
0,05 0,0019 0,3 0,035 1,4 0,1254
0,07 0,0036 0,35 0,0397 1,6 0,14
0,09 0,006 0,4 0,0444 1,8 0,149
0,11 0,009 0,45 0,049 0,157
0,13 0,0126 0,5 0,0535 2,2 0,165
0,15 0,0167 0,6 0,0625 2,4 0,172
0,17 0,0214 0,7 0,0707 2,6 0,179
0,19 0,0239 0,8 0,0789 2,8 0,186
0,21 0,026 0,9 0,087 0,192