Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Расчет геометрии трещины при проведении ГРП. Методические указания

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

L

p(x,t)dx

 

 

 

k

 

 

 

 

 

 

 

 

 

 

L

x

2

 

2

H

2L

.

0

 

 

 

 

2

 

 

 

 

 

 

(4)

Уравнение (3) предполагает одностороннюю трещину, взятую при полной производительности q0. Как и в первоначальной теории РК сделано допущение, что

 

 

q

= 0.

 

 

 

 

x

 

 

 

Желтов

и

Христианович

[2],

предложили

распределение давления жидкости, которое удовлетворяет уравнению (4) аппроксимировать при:

рf = p0; для 0 L0/L;

рf = 0; для L0/L 1;

с отношением L0/L замыкающемся на единице. Рассчитаем "спорную" длину трещины:

 

 

 

H

 

K c

 

 

 

0 = L0/L = sin

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

p f

 

p f

 

 

 

 

 

 

L

За исключением непосредственной близости к вершине трещины, форма трещины в горизонтальной плоскости имеет вид эллипса аксиальной шириной трещины у скважины:

(0,t) 2(1 )L( p f H ) .

G

Это соотношение подобно уравнению (1) в теории РК, если принять L = hf/2.

Определим уровень Pf, обусловленного быстрым снижением давления жидкости у вершины трещины, которая заключается в форме клина. Хорошее

11

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

L0

 

3

(0,t) dL

7

1 L20

1

 

 

 

 

 

 

2 .

 

3

4

0

 

(x,t)

 

 

 

После подстановки этого выражения в уравнение (3) реализуя дальше то, что объем односторонней трещины примерно равен:

 

1

 

 

 

 

 

 

V hf

2

)

1 / 2

d

hf

L (0, t) (1

 

4

 

0

 

 

 

 

 

 

 

 

 

 

 

 

окончательно получаем: L(t) = 0,68

 

 

1 q

3

 

(0,t) 1,87

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

Ghf

 

 

 

 

 

 

 

 

 

 

 

L (0, t)

 

Gq

3

 

0

 

 

 

 

(1 )h

 

 

 

 

 

t

1 / 3

.

 

 

 

 

 

 

 

q

t

0

 

 

 

3

t

 

f

.

2 / 3

;

Обобщения результатов по двум моделям

Таблицы 1, 2 обобщают результаты для длины трещины, максимальной ширины трещины и давления для моделей РКN и GdK. Сравнения для длин трещин, максимальной ширины и давления нагнетания при постоянной производительности нагнетания.

 

 

 

 

 

 

Таблица 1

 

 

 

 

 

 

 

 

РКN

ОДНО

ДВА

 

GDK

ОДНО

 

ДВА

 

КРЫЛО

КРЫЛА

 

 

КРЫЛО

 

КРЫЛА

C1

0,68

0,45

 

C4

0,68

 

0,48

C2

2,50

1,89

 

C5

1,87

 

1,32

C3

2,75

2,31

 

C6

2,27

 

1,19

 

 

 

 

 

 

Таблица 2

 

 

 

12

 

 

 

 

Модель

РКN

GdK

Длина трещины

L(t)

 

 

 

 

 

 

 

Gq

3

 

 

1 / 5

 

 

 

 

 

 

 

 

 

 

 

 

C

 

0

 

 

 

t

4 / 5

(1 )h

4

 

1

 

 

 

 

f

 

 

 

 

 

 

 

 

L(t)

 

 

 

 

 

 

 

Gq

3

 

 

1 / 6

 

 

 

 

 

 

 

 

 

 

 

 

C4

0

 

 

 

t

2 / 3

(1 )h

3

 

 

 

f

 

 

 

 

 

 

 

 

Ширина трещины

(0,t)

 

 

 

 

 

 

 

 

1 q

2

 

1 / 5

 

C2

 

Ghf

0

 

t1 / 5

 

 

 

 

 

 

 

(0,t)

 

 

 

 

 

 

 

1 q

 

3

 

1 / 6

 

 

 

 

 

 

C5

 

 

0

 

 

 

t

1 / 3

Ghf

 

 

 

 

 

 

 

 

 

 

 

 

Избыточное

давление

P(0, t)

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

Gq

 

3

L

1 / 4

 

 

 

 

 

3

 

 

 

0

 

 

 

 

 

2H

 

(1 )

3

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

P(0,t)

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

Gq h

3

 

1 / 4

 

 

 

6

 

f

 

 

 

 

 

0

 

 

 

 

 

2H

 

(1 )

3

 

2

 

f

 

 

L

 

 

 

 

 

 

 

 

 

 

 

Влияние поглощения на размеры трещины

Основные элементы описания влияния потерь жидкости на распространение трещины взяты по Картеру[4]. В его теории принято, что ширина и высота трещины должны быть постоянными, переменной является только длина трещины. Производительность нагнетания также принята постоянной.

Функция скорости поглощения жидкости, нормальной к поверхности трещины, принята в виде:

v

 

 

K

L

 

 

 

 

 

 

 

L

 

t

 

 

 

,

(4)

где КL коэффициент общего поглощения (потерь) жидкости, измеренный при лабораторных испытаниях фильтрации.

отражает время, с которого фильтрация началась. Для геометрии материального баланса, имеем:

dVdt hf dLdt q qt ,

где

13

 

Afl

Afl

dA

Afl

dA

d

 

 

 

 

qL

vL dAfl

vL

fl

d KL

fl

 

d

d

1

 

0

0

0

и Afl=2hfL.

Эта формула дает:

q

tl

dL

d

 

dL

 

 

 

 

 

 

h

 

2K L

d

1

 

dt

f

0

 

 

 

 

 

 

 

которое приводит к

 

qt

 

2

 

 

 

 

2

 

 

 

 

L , где

2K

 

t

V =

 

 

e

 

erfc

 

1 hf

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для > 4

 

e

 

2

erf

 

1

 

(см. таблицу 5),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

атакже для больших значений: hf L 2q0t ,

.

которое приводит к: L

q0

 

 

t

 

 

 

 

 

h

f

K

L

 

 

 

 

 

 

 

 

 

 

 

Если эффективность жидкости определена как

отношение объема разрыва

hf

L

 

к нагнетаемому объему

q0t, то для больших значений :

 

 

 

 

 

 

 

 

nf =

h

 

L

 

2

.

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

t

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Это уравнение подчеркивает важность удержания как можно меньшим. Для =10 эффективность составляет только 11%.

Таким образом, распределение скоростей в трещине зависит от эффективности жидкости, посредством этого воздействуя на распределение давления.

14

Обозначим безразмерный параметр

Al

hf

LKl

/ q0 .

Для Аl> 1 теория показывает, что скорость течения резко снижается в соответствии с

q(x)

 

v(x)

q

 

v(0)

0

 

 

 

 

1

2

arcsin

 

 

 

.

Для этих условий коэффициент Аl соотносится с эффективностью жидкости согласно зависимости:

 

h L K

 

2

t

 

 

 

2

 

1

 

A

 

 

 

 

 

 

 

 

 

 

 

.

f

 

l

 

 

 

4

 

 

 

 

 

 

 

 

 

l

q

t

 

2

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

f

(5)

Тогда потерю струи представим в виде:

 

 

dV

q ql

Vsp

dA

 

 

 

 

 

fl

;

 

 

dt

dt

 

 

 

 

 

 

и тогда

 

 

 

dL

 

 

dL

d

q

 

 

 

t

 

 

 

 

 

 

 

 

h

 

( 2Vsp )

dt

2Kl

d

t

f

 

 

 

0

 

 

 

 

 

 

 

 

.

Сравнение уравнений показывает, может быть принята в расчет при эффективную ширину трещины ( 2Vsp

что потеря струи

замене

 

через

), таким образом:

 

2K

l

 

t

 

 

 

 

 

 

 

( 2V

sp

)

 

 

 

 

.

Дополнение в модель РКN.

Перкинс и Керн не принимали во внимание поглощение (утечки) жидкости в породу недостаток возмещенный Нордгреном [1]. Уравнение местной непрерывности теперь получило вид:

q

 

hf

 

ql

0 ,

x

 

 

t

4

 

 

 

 

 

 

 

15

 

 

 

(для сравнения:

q

 

h

f

 

 

 

x

 

4

 

t

),

где ql потери на единицу длины трещины:

q

l

 

 

2h K

 

f

l

 

 

t (x)

.

Расширенная

форма

уравнения

 

 

G

 

 

 

 

 

 

64(1

)h

 

 

 

f

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

0

x

2

 

t

 

 

 

 

 

 

 

 

 

 

 

 

G

 

64(1 )h

 

f

 

теперь принимает вид:

 

 

 

 

 

8K

 

 

2

 

4

 

 

 

 

 

0

x

 

 

t

 

 

l

 

 

 

 

 

2

 

t (x)

 

 

 

 

 

 

 

 

При начальном условии (х,0)=0 и граничных условиях

(x,0)=0 для х > L .

 

4

(0,t)

 

256 (1 )q

 

 

 

x

G

 

 

для односторонней трещины.

qx

Система уравнений может быть решена, принимая

ql

0

, получаем

1

d

 

q

 

, где зависит от

 

 

 

L

 

 

0

 

 

 

 

 

 

 

 

 

 

 

0

t (x)

 

2h

K

l

 

 

 

 

 

 

f

 

=x/L, так что = (L(t )) = t для 0<t <t

Это решение подобно решению Картера:

L

1 q

0

t

 

 

 

 

 

 

 

 

h

K

l

 

 

f

 

и

q( )

q

0

 

1 arcsin

.

Это уравнение является хорошим приближением для

h

LK

2

q

(0)

.

 

 

 

 

 

 

l

 

 

 

 

 

 

f

 

 

0

 

 

 

 

 

 

 

 

 

Для этих условий:

 

 

 

 

 

 

 

4

 

G3q L

 

 

1

2

 

p4 (0)

 

 

 

 

 

0

 

 

 

(1

 

arcsin )d ,

 

4

(1

 

3

)

 

 

 

 

 

 

h

 

 

0

 

f

что дает в итоге

16

(0) 4( 2 )3 (1 )q2 1/ 4 t1/ 8 и

Gh K 0

f l

p(0) 4(

2

 

3

 

 

 

 

)

 

 

 

 

 

h

 

 

 

 

 

f

 

 

 

 

 

Gq

 

2

 

0

 

(1 )K

l

 

 

1/ 4

 

 

1/ 8

 

t

 

 

 

 

.

Дополнение в модель GdK.

В отличие от модели РК первоначальная модель GdK принята для моделирования условий месторождения более подходящих для малых значений т.е. для малых значений коэффициента поглощения и малого времени обработки. Материальный баланс учтен в общей форме также, как и у Картера, включая потери струи:

dV

q q

 

2V

 

h

 

dt

l

sp

f

 

 

 

 

 

 

 

 

 

dL dt

.

Объем трещины отличается от объема по Картеру изза эллиптической формы поперечного сечения в горизонтальной плоскости, таким образом:

 

dV

 

 

 

 

 

dL

L

d (0)

 

 

 

 

 

 

 

 

 

 

L

d (0) dL

 

 

 

 

 

 

 

hf

(0)

 

 

 

 

 

 

 

 

 

 

 

hf

 

(0)

 

 

 

 

.

 

dt

 

4

 

 

 

dt

 

 

 

dt

 

 

4

 

 

 

 

 

 

 

 

 

dL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

Для модели GdК было найдено, что L

d (0)

 

(0)

 

 

dL

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0)

C2

 

L1/ 2

Из L C1t

 

и

(0) C2t

 

 

 

следует, что

 

 

 

2 / 3

1/ 3

 

 

 

 

 

C1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d (0)

 

 

 

C2

 

 

1/ 2

 

(0)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dL

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Материальный баланс запишется:

 

 

 

 

 

 

 

 

 

 

 

 

 

q

0

 

 

3

 

 

 

 

 

dL

 

 

 

t

 

dL

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0,t)

2Vsp

 

 

 

 

 

2Kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hf

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

dt

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это решение может быть использовано потому, что в

произведении

dL dt

изменения в виде функции времени

намного меньше, чем в выражении

dL dt

.

Однако, постоянное значение должно быть выбрано таким образом, чтобы правильное решение не «перекликалось» с поглощением жидкости: (0,t) в уравнении должно быть заменено на 2/3 (0,tp), где tp время остановки насосного оборудования. Коэффициент 2/3 делает приблизительное решение численно равным единственному точному для Кl=Vsp=0:

dL

 

q

dt

h

 

0 f

 

 

4

 

 

(0,t

 

) 8V

 

 

p

sp

 

 

 

2

L

erfc L

e

 

 

 

 

 

 

 

, где

 

 

 

8K

L

 

t

 

 

 

 

 

) 8

 

L

 

(0,t

 

 

 

 

 

 

 

 

 

 

 

 

p

 

Интегрирование уравнения показывает, что для односторонней трещины:

L

q0

 

(0,t p ) 8Vsp (

2 L

 

1 e 2 L erfc L ) и

16 h

K 2

 

 

 

 

 

 

 

 

 

 

 

 

 

f

L

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 ) q L

 

1/ 4

 

 

 

(0,t p ) 2,27

 

 

 

2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ghf

 

 

 

Соотношение между L и для определения обеих неизвестных можно найти как расчетным, так и графическим путем.

Основываясь на этом, рассчитаем геометрию трещины при проведении гидроразрыва пласта без учета потерь жидкости и с учетом потерь жидкости [3].

Примем известными:

18

Qi (производительность течения жидкости гидроразрыва на входе в трещину) = 2 м3/мин;

μ (вязкость жидкости разрыва) = 1,67·10-6 КПа·мин; h (мощность пласта) = 22 м;

Е (модуль Юнга) = 2400000 КПа;

Ke (полный коэффициент потерь жидкости)

Ke = 0,0026 м/мин0,5;

ν (коэффициент Пуассона) = 0,2 ед.; t (время работ по закачке) = 120 мин.

В таблице 3 приведены значения модуля Юнга для различных образцов горных пород.

Таблицы 4 и 5 обобщают результаты расчетов геометрии трещины для моделей РКN и GdK при постоянной производительности нагнетания.

19

 

 

 

 

 

Таблица №3

 

 

 

 

 

 

 

 

Описание

 

 

Модуль

 

 

 

 

 

 

Юнга

 

 

 

 

 

 

(МПа)

 

 

 

Песчаник Ойл-Грик, массивный, очень твердый,

 

251

весьма тонкозернистый, хроно сцементированный

 

 

 

 

 

Песчаник Бартелессвил, массивный, тонкозернистый,

 

293

хорошо сцементированный

 

 

 

 

 

 

 

Песчаник массивный, однородный, среднезернистый

 

223

 

 

 

 

 

 

 

Известняк

 

 

 

 

 

629

Известняк

ЛЕС

Невада-сайт,

плотный,

 

764

тонкозернистый

 

 

 

 

 

 

 

 

Доломит Блайр, однородный, тонкозернистый

 

443

 

 

 

Доломит Фусильмен, однородный, тонкозернистый

 

410

Доломит Хасмарк, однородный, грубозернистый,

 

245

насыщенный

 

 

 

 

 

 

 

 

 

 

 

 

 

Доломит

Стоктон

и

брекчия

доломита,

 

267

известковистый, от средне до тонкозернистого

 

 

 

 

 

Доломит со стилолитами, с накоплением глин

 

619

 

 

 

 

 

Ангидрид, тонкозернистый

 

 

 

172

Девонский известняк, однородный, грубозернистый

 

252

 

 

 

Известняк Вольф Кемп, однородный, тонкозернистый

 

413

 

 

 

 

 

 

Известняк

Вемс

Стейшн,

однородный,

 

468

мелкозернистый

 

 

 

 

 

 

 

 

 

 

 

 

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]