Расчет геометрии трещины при проведении ГРП. Методические указания
.pdf
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
