Разработка нефтяных месторождений
.pdfvk.com/club152685050 | vk.com/id446425943
ти, воды и породы при снижении давления. Аналитически, согласно закону Гука, можно записать:
V |
V |
V |
V (m |
æ |
|
) p V p. |
(2.4) |
óç |
æ |
ï |
ï |
ñ |
ï |
|
ãäå æ, ñ коэффициент сжимаемости жидкости и скелета породы; Vï геометрический объем пласта; ð снижение пластового давления ( ðmax ðïë ðíàñ).
На начальном этапе создания теории упругого режима в работах М. Маскета, Т. Мура и др. при определении упругоемкости пласта сжимаемостью породы пренебрегали и полагали, что
|
|
ò æ. |
(2.5) |
|
Наиболее полно теория упругого режима разработана В.Н. Щелкачевым. Необходимость учета сжимаемости породы была им показана на основании результатов экспериментальных исследований Д.А. Антонова (1952 г.). К аналогичному выводу пришел и американский ученый Х. Холл в 1953 г.
Упругие деформации могут быть обратимыми и необратимыми. Экспериментально это было подтверждено исследованиями В.М. Добрынина. Теоретическое описание упругих деформаций в пласте было сделано Г.И. Баренблаттом и А.П. Крыловым, а режим назван упруго-пласти- ческим (по определению И.Н. Стрижова упруго-гру- зовой режим). Дальнейшее развитие теории этого процесса было осуществлено В.Н. Николаевским и А.Т. Горбуновым.
Изменение давления при упругом законе описывается дифференциальным уравнением пьезопроводности по аналогии с уравнением теплопорводности Фурье, известным в математической физике с начала 19 в. Уравнение пьезопроводности имеет вид
2p 1 |
p , |
(2.6) |
|
t |
|
где коэффициент пьезопроводности горной породы, характеризующий распределение давления в пласте.
41
vk.com/club152685050 | vk.com/id446425943
Коэффициент пьезопроводности был предложен В.Н. Щелкачевым по аналогии с коэффициентом температуропроводности в уравнении Фурье.
Для плоскорадиальной фильтрации уравнение пьезо-
проводности имеет вид |
|
||||||
2p 1 |
p 1 |
p ; |
|
|
(2.7) |
||
r2 |
r |
r |
|
t |
|
|
|
с учетом инерционных сил |
|
||||||
2p 1 ð |
æ æ |
2p ; |
(2.8) |
||||
|
t |
|
|
t2 |
|
|
|
и в случае неизотермической фильтрации |
|
||||||
k |
|
|
ð |
|
T |
, |
(2.9) |
|
p |
t |
a |
z |
|||
|
|
|
|
|
|
||
ãäå à коэффициент, учитывающий объемное расширение жидкости и среды.
Такое решение уравнения пьезопроводности при определенных начальных и граничных условиях было получе- но А. Ван-Эвердингеном и У. Херстом в виде интегрального выражения.
Введем понятие безразмерного времени
ñ |
|
kt |
|
. |
(2.10) |
|
|
|
|
|
|||
m |
|
|
r2 |
|||
|
|
æ |
|
æ c |
|
|
Ïðè ñ 100 решение уравнения пьезопроводности
можно представить в виде |
|
pc( c) 0,5 ln 0,80907 . |
(2.11) |
При разработке нефтяных месторождений на упругом режиме выделяют две фазы. Первая фаза существует до тех пор, пока возмущение от пуска скважины в работу не достигнет контура нефтеносности (т.е. когда ðê const). Длительность этой фазы равна
t |
Rê2 |
rc2 |
. |
(2.12) |
|
|
|||
1 |
4 |
|
||
|
|
|||
42 |
|
|
|
|
vk.com/club152685050 | vk.com/id446425943
Положение границы области пониженного давления оценивают по формуле Э.Б. Чекалюка:
R(t) rc |
t. |
(2.13) |
Распределение давления в приведенной области влияния скважины можно определить с использованием многих формул.
Формула Дюпюи:
p p |
|
q |
ln |
rïð |
, |
|
|
|
|
|
|
|
|
|
(2.14) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
0 |
c |
|
|
2 kh rc |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
ãäå rïð |
2 t, |
(по И.А. Чарному). |
|
|||||||||||||||
|
Формула А.М. Пирвердяна |
|
||||||||||||||||
p p |
|
q |
ln |
rïð |
|
q |
|
|
1 |
rïð |
|
, |
(2.15) |
|||||
|
|
|
|
|
|
|
||||||||||||
0 |
c |
|
|
2 kh |
|
r |
|
|
2 kh |
|
|
|
r |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
c |
|
|
|
|
|
|
c |
|
|
|||
ãäå rïð |
12 |
t. |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Формула Г.П. Гусейнова |
|
|
|
|
|||||||||||||
|
|
|
|
q |
|
rïð |
|
|
q |
|
|
|
|
r2 |
|
|
|
|
p0 pc |
|
|
|
ln |
|
|
|
|
|
|
1 |
|
c |
|
|
, |
(2.16) |
|
|
2 kh |
rc |
|
|
|
|
r2 |
|
||||||||||
|
|
|
|
|
|
|
4 kh |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ïð |
|
|
||
ãäå rïð |
8 t. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Упрощенная основная формула упругого режима |
|||||||||||||||||
|
|
|
|
q |
|
|
|
r2 |
|
|
|
|
|
|
|
|||
p p |
|
|
E |
|
|
, |
|
|
|
|
|
(2.17) |
||||||
|
2 kh |
4 t |
|
|
|
|
|
|||||||||||
0 |
c |
|
|
|
i |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
ãäå ð0 начальное пластовое давление; ðñ давление на забое добывающей скважины; rïð приведенный радиус скважины.
Сопоставление результатов расчетов по указанным формулам, сделанные В.Н. Щелкачевым, показало, что наилучшие результаты дает упрощенная основная формула упругого режима, которая справедлива для бесконечного пласта.
43
vk.com/club152685050 | vk.com/id446425943
Изменение давления во второй фазе упругого режима можно определить по нескольким по формулам:
Формула И.А. Чарного для определения давления на стенке скважины и на контуре нефтеносности (или непроницаемой границе):
|
|
|
q |
|
4 |
|
|
|
q |
|
t |
|
|
|
|
p |
p |
|
ln |
t |
1 |
|
|
|
; |
(2.18) |
|||||
4 kh |
|
kh |
R2 |
r2 |
|||||||||||
c |
0 |
|
r2 1 |
|
|
|
|
|
|||||||
|
|
|
|
|
c |
|
|
|
|
ê |
|
c |
|
|
|
pk p0 q t ,
kh Rê2
ãäå t1 продолжительность первой фазы.
Формула Г.И. Баренблатта для определения пласто-
вого давления в точке r в момент времени t ïðè t 1 :
Rê2 12
|
|
q |
|
|
R |
r |
|
7 |
|
t |
|
|
|
p p |
|
|
|
ln |
ê |
|
|
|
2 |
|
|
. |
(2.19) |
|
|
|
|
||||||||||
0 |
|
|
|
r |
Rê |
|
|
|
R2 |
|
|
||
|
|
2 kh |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
ê |
|
|
|
При расчетах пластового давления и дебита скважин для месторождения в целом пользуются понятием укрупненной скважины. Радиус такой скважины равен
Róñ NF ,
ãäå F площадь месторождения; N число скважин.
|
В расчетах для укрупненной скважины, когда |
Rê |
|
|
|
r |
|
|
|
c |
|
|
1,5 2, пользоваться строгой гидродинамической форму- |
||
лой нельзя и следует использовать формулы М. Маскета, А. Хупера, А. Ван-Эвердингена, У. Херста. Так, формула М. Маскета имеет вид
ð (t) ð |
|
Q í |
p( ), |
(2.20) |
|
|
|||||
c |
0 |
|
2 kh |
|
|
|
|
|
|
|
|
ãäå |
t . |
|
|
||
|
r2 |
|
|
|
|
|
c |
|
|
|
|
Функция p( ) подсчитана А. |
Ван-Эвердингеном и |
||||
У. Херстом для значений 1 . |
|
||||
44 |
|
|
|
|
|
vk.com/club152685050 | vk.com/id446425943
Для расчета дебита скважины в круговой залежи ограниченного пласта в период действия первой фазы упругого режима предпочтительнее пользоваться следующими формулами:
Э.Б. Чекалюка
q |
q |
|
|
|
|
1 |
|
|
, |
(2.21) |
2 kh (ð |
ð ) |
ln(1 |
|
f |
) |
|||||
|
0 |
c |
|
|
|
0 |
|
|
|
|
ãäå q è q соответственно безразмерный и размерный дебит
скважины; f |
t |
параметр Фурье; |
|
0 |
r2 |
|
|
|
c |
|
|
Б.А. Богачева |
|
|
|
q 4 kh (ln f0 1)(p0 |
pc) . |
(2.22) |
|
(ln f )2 |
|
|
|
|
0 |
|
|
Погрешность этой формулы не превышает 1,2 %.
Для расчета дебита в период действия второй фазы режима рекомендуется пользоваться формулой В.Н. Щелка- чева (при Rê rc с точностью до 3 8 %)
|
|
|
|
|
|
|
|
|
2 |
F |
|
|
|
|
q |
|
|
|
e |
|
D 0 |
, |
(2.23) |
||
q |
|
|
|
|
|
|
|
|
||||
2 kh(p |
p ) |
ln |
Rê |
|||||||||
|
|
0 |
c |
|
r |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||
c
ãäå D ln Rk 1.
rc
Упруговодонапорный режим (УВНР) наиболее распространенный режим в практике разработки нефтяных месторождений. Он может быть как искусственным, так и естественным. В России доля добычи нефти при разработке на естественном УВНР составляет около 6 %. На искусственном более 90 % всей добываемой нефти извлекается из терригенных коллекторов, оставшаяся часть приходится на карбонатный коллектор. Искусственный УВНР на 98 % обеспечивается внутриконтурным заводнением, в том числе 10 % за счет площадного и избирательного.
Подавляющее большинство объектов ( 70 %) разрабатывается при плотности сетки скважин 10 50 га/скв. На долю этих объектов приходиться около 90 % добываемой
45
vk.com/club152685050 | vk.com/id446425943
за счет заводнения нефти. Наиболее распространенные сетки скважин 20 30 га/скв.
Рассмотрим особенности естественного УВНР.
Пуск скважины с постоянным дебитом в круговой залежи неограниченного пласта. Перераспределение давления в пласте после пуска добывающей скважины можно рассчитать с помощью формулы, которая давно известна в математической физике, но для подземной гидродинамики впервые применена Тсейсом.
|
|
|
|
|
|
|
q |
|
|
|
|
4 t |
|
|
|
|
r |
2 |
|
1 |
|
r |
2 |
|
2 |
||
ð0 |
ð |
|
|
|
|
|
0,5772 |
|
|
|
|
||||||||||||||||
|
|
|
|
ln |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
4 kh |
|
|
r |
|
|
|
|
|
4 t |
|
4 |
|
4 t |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
r |
2 |
|
|
3 |
1 |
|
|
r |
2 |
|
4 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
... |
, |
|
|
|
|
|
|
|
|
(2.24) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
18 |
|
|
|
|
|
|
19 |
|
4 t |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
4 t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ãäå r расстояние от скважины до точки с давлением ð; q установившийся дебит скважины.
Используя понятие упрощенной интегральной экспоненциальной функции Ei( x), получим
|
|
|
|
|
|
q |
|
|
|
|
|
|
|
|
r2 |
|
|
|
|
|
|
|
|
|
|
|||||
ð |
ð |
|
|
|
|
|
E |
|
|
|
|
, |
|
|
|
|
|
|
|
|
(2.25) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
0 |
|
|
|
|
4 kh |
|
|
|
i |
|
4 t |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ãäå |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 t |
|
|
|
|
|
|
|
r |
2 |
|
1 |
|
r |
2 |
|
2 |
|||
Ei( x) |
|
|
|
0,5772 |
|
|
|
|
||||||||||||||||||||||
ln |
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
4 t |
|
4 |
|
4 t |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
r |
2 |
|
|
3 |
|
1 |
|
|
r |
2 |
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
... . |
|
|
|
|
|
|
|
|
(2.26) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
18 |
|
|
|
|
|
|
|
19 |
|
|
4 t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
4 t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Точное решение уравнения пьезопроводности для укрупненной скважины в общем виде в случае развития упругого режима в законтурной области и поступления воды в нефтяную зону было получено А. Ван-Эвердингеном и У. Херстом (при q const, ð ð0 ïðè t 0, R r ).
ð ð( , ) |
q |
p( , ), |
(2.27) |
|
|||
0 |
4 kh |
|
|
|
|
|
|
46 |
|
|
|
vk.com/club152685050 | vk.com/id446425943
ãäå
|
|
2 |
(1 e u2 ) (u)Y (u, ) Y (u)J (u, ) |
|
|
|||||||||
ð ( , ) |
|
|
|
|
1 |
|
0 |
1 0 |
|
du; |
(2.28) |
|||
|
|
|
u |
2 |
2 |
2 |
|
|
||||||
|
|
|
0 |
|
|
|
J |
(u) Y |
(u) |
|
|
|
||
|
|
|
|
|
|
1 |
1 |
|
|
|
|
|||
|
R ; xt ; |
|
|
|
|
|
|
|
|
|||||
|
r |
|
r |
|
|
|
|
|
|
|
|
|||
|
c |
|
|
|
c |
|
|
|
|
|
|
|
|
|
u |
R2 |
|
|
|
; |
|
|
|
|
|
|
|
|
|
4 (t t ) |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
ãäå J0(u, ); J1(u); Y0(u, ); Y1(u) функции Бесселя. Для 1 функция p ( , ) табулирована.
Для 1 можно воспользоваться аппроксимацией Ю.П. Желтова.
p (1, ) 0,5 |
|
3,81 |
|
0,487 ln(1 ). |
(2.29) |
1 |
(1 ) |
|
Имеются также решения М. Маскета, В.Н. Щелкачева, Г.И. Баренблатта, А.М. Пирвердяна, И.А. Чарного, М.А. Гусейн-заде, А.К. Колосовской и др.
Пуск скважины с постоянным забойным давлением. Дебит скважины и накопленную добычу можно определить по формулам М. Маскета (для галереи)
q |
2 kh |
|
R2(ð |
ð )q( ); |
|
|
(2.30) |
||||||
|
|
|
|
||||||||||
|
|
ã |
0 |
ã |
|
|
|
|
|
|
|
||
Q |
|
2 kh |
R2 |
(ð |
|
|
|
|
|
|
|
|
|
|
ð )Q |
( ), |
|
|
(2.31) |
||||||||
|
|
|
|
||||||||||
|
|
|
ã |
0 |
ã |
|
|
|
|
|
|
|
|
|
|
( ) è |
q ( ) |
безразмерные функ- |
|||||||||
Rã радиус галереи; Q |
|||||||||||||
öèè |
дебита |
и накопленной |
добычи (табулированы |
||||||||||
Ó.Херстом, А. Ван-Эвердингеном).
Учет различия гидропроводностей нефтяной и законтурной областей. По предложению А.П. Крылова для расчетов при УВНР пользуются поправочным коэффициентом:
Z |
pôàêò(t) |
. |
(2.32) |
|
|||
|
p |
|
|
|
òåîð |
|
|
|
|
|
47 |
vk.com/club152685050 | vk.com/id446425943
В.С. Орловым было предложено использовать два коэффициента: Z1, учитывающий различие гидропроводностей, и Z2, учитывающий различие пьезопроводностей нефтяной и водонасыщенной частей залежи.
p(t) |
|
|
|
|
r2 |
|
|
|
|
Z E |
|
Z |
. |
(2.33) |
|||
|
|
|||||||
|
4 kh |
1 |
i |
|
4 t |
2 |
|
|
|
|
|
|
|
|
|
||
Для расчета давления на контуре залежи при изменяющейся во времени добыче жидкости используют теоремы суперпозиции Дюамеля. Существует две теоремы суперпозиции Дюамеля. Согласно первой теореме, криволинейную зависимость дебита (рис. 2.1) во времени можно заменить ступенчатой ломаной линией, т.е. участками с постоянным дебитом.
После пуска скважины в момент с дебитом q изменение давления на контуре питания будет равно
pê(t) ð0 2 kh 1 (q1 q0)f(1, t)d
0
2 |
|
|
(q2 |
q1) f(1, t 1)d . |
(2.34) |
|
|
|
1 |
|
|
Согласно второй теореме, при работе нескольких скважин (рис. 2.2) общее изменение давления в точке Ì в момент времени t будет равно
ð(r, t) ðA(r1, t) ðB(r2, t) ðC(r3, t), |
(2.35) |
ãäå ðÀ (r1, t) p0 pA(r2, t) падение давления в т. Ì от работы скважины À; ð (r2, t) p0 pB(r2, t) падение давления в т. Ì от работы скважины Â; ðÑ
(r3, t) p0 pÑ(r3, t) падение давления в т. Ì от работы скважины Ñ.
Прогнозирование давления на контуре питания при упругом режиме в законтурной области. Предполагается, что добыча жидкости во времени изменяется в соответствии с графиком (рис. 2.3).
48
vk.com/club152685050 | vk.com/id446425943
Рис. 2.1. Зависимость дебита во Рис. 2.2. Схема залежи
времени:
текущее время; время ступенчатого изменения дебита
Согласно первой теореме Дюамеля, изменение давления на контуре питания будет определяться выражением (при
Rê 1)
rc
ð |
(t) ð |
|
|
|
q f (1, )d . |
(2.36) |
|
2 kh |
|||||||
ê |
0 |
|
|
|
|||
|
|
|
|
0 |
|
|
|
Характер внедрения воды в нефтяную залежь является типичным, и можно выделить четыре периода.
1. Увеличение поступления законтурной воды, qçâ , 0 1 , ãäå 1 конец разбуривания.
Рис. 2.3. Зависимость добычи
жидкости во времени:
qi const прирост добычи жидкости
49
vk.com/club152685050 | vk.com/id446425943
2. Стабилизация отбора и поступления жидкости. qçâ q1 const,
1 .
где начало законтурного заводнения.
3. Снижение поступления законтурной воды за счет законтурного заводнения.
qçâ q1 1 , ïðè .
4. Установившийся приток воды в залежь. qçâ q2 const, ïðè .
Интеграл Дюамеля для 1 вычисляется с помощью аппроксимации Ю.П. Желтова
J( ) 0,5 0,178 |
|
|
1 |
|
|
1 |
|
|
|||
2,81 |
|||||
|
|
|
(1 ) |
|
|
0,487[(1 )ln(1 ) ].
Тогда для 1-го периода
pê(t) p0 J( );
2 kh
для 2-го периода
pê(t) p0 [J( ) J( 1)];
2 kh
для 3-го периода |
|
|
|
|
|
|||
p (t) p |
|
[J( ) J( )] |
1 |
J( ) |
; |
|||
|
|
|||||||
ê |
0 |
|
|
1 |
2 kh |
* |
|
|
|
|
2 kh |
|
|
|
|
||
для 4-го периода
pê(t) p0 [J( ) J( 1)]
2 kh
1 [J( ) J( )].
2 kh
50
(2.37)
(2.38)
(2.39)
(2.40)
(2.41)
