
Лекции Сумин А Н 5
.pdf
РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
предыдущему скорости проникновения нефти в элемент пласта и выхода из него,
получим
|
|
|
|
|
|
vн m |
s |
0 , |
|
|
|
(5.44) |
|
|
|
|
|
|
|
x |
t |
|
|
|
|
|
|
Складывая уравнения (5.43) и (5.44), имеем |
|
|
|
|
|
||||||||
|
|
v |
|
v |
|
0 ; |
|
v |
|
v |
|
v t , |
(5.45) |
|
|
н |
в |
|
н |
в |
|||||||
|
x |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, суммарная скорость фильтрации нефти и воды не изменяется по координате x , что и следовало ожидать, так как нефть и воду принимают за несжимаемые жидкости.
Следовательно, режим пласта жесткий водонапорный.
Скорости фильтрации воды и нефти подчиняются обобщенному закону Дарси, так что
vв |
|
kkв s |
|
p |
; |
|
|
vн |
|
kkн s |
|
p |
, |
(5.46) |
||||||||||
в |
|
|
x |
|
|
|
н |
|
x |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где k в и k н , в и |
н - |
относительные |
|
проницаемости, |
зависящие от |
|||||||||||||||||||
водонасыщенности s и вязкости воды и нефти. |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Рассмотрим функцию f s , |
называемую функцией Бакли - |
Леверетта. При |
||||||||||||||||||||||
этом |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f s |
|
vв |
|
|
|
|
|
|
kв s |
|
|
|
|
, |
|
|
|
|
(5.47) |
|||||
vв vн |
kв s |
в |
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
н |
kн s |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f s |
|
vв |
|
. |
|
|
|
|
|
|
|
|
|
|
|
(5.48) |
||||
|
|
|
|
|
v t |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Из (6.48), дифференцируя vв по x , получим: |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
vв |
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
s x , |
|
|
|
|
|
|
|
|||||
|
|
|
|
x v t f |
|
|
|
|
|
|
(5.49) |
После подстановки (5.49) в (5.43) получим одно дифференциальное уравнение первого порядка для определения s , т.е.
108

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
|
|
s |
|
s |
|
|
v t f |
s |
|
m |
|
0, |
(5.50) |
x |
t |
|||||
По мере вытеснения нефти водой |
из |
прямолинейного пласта |
фронт |
вытесняющей нефть воды продвигается к концу пласта и водонасыщенность в каждом сечении заводненной области непрерывно увеличивается. Процесс вытеснения нефти водой из прямолинейного пласта можно представить и иным образом, следя за изменением по пласту некоторой водонасыщенности. Если,
например, в какой-то момент времени в некотором сечении пласта водонасыщенность составляла s s1, то спустя определенное время эта водонасыщенность будет и в конце пласта, так как нефть постепенно извлекается из него и ее место занимает вода. Для указанного s const можно принять
ds |
s |
|
dx |
s |
dt 0 |
|
|||||||||||
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
x |
|
|
|
t |
|
|||||||
или |
|
|
|
|
|
|
|
|
|||||||||
|
s |
|
|
dx |
|
s |
0 , |
(5.51) |
|||||||||
|
|
|
|
t |
|||||||||||||
|
x dt |
|
|
|
|
|
|
||||||||||
Сравним (5.50) и (5.51). Они будут идентичными, если положить |
|
||||||||||||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
f s v t |
, |
(5.52) |
|||||||||||||
|
|
|
|
||||||||||||||
|
dt |
m |
|
|
|
|
|
|
|||||||||
Умножим и разделим (5.52) на bh и проинтегрируем, получим: |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
Qвз t bhv t dt , |
(5.53) |
|||||||||||
bhmx f s Qвз t ; |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
||||
Обозначим |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
bhmx |
|
|||||||||||
|
|
|
|
|
|
|
, |
(5.54) |
|||||||||
|
|
|
|
Q |
вз |
t |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
тогда
(5.55)
Задавая s в формуле (5.55), можно определить расстояние от входа в пласт для данного значения водонасыщенности. Однако в период безводной эксплуатации закачиваемая вода еще не достигает конца пласта. Чтобы установить положение фронта вытеснения нефти водой и водонасыщенность на
109

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
фронте вытеснения, рассмотрим материальный баланс закачанной в пласт воды.
Если к моменту времени t в пласт закачан объем воды, равный Qвз t , длина фронта вытеснения составит xв , насыщенность пласта связанной водой s sсв ,
то
|
|
|
|
|
|
|
|
xв |
|
|
|
|
|
|
|
|
|
|
|
Qвз t bhm s x dx bhmxвsсв . |
|
(5.56) |
|||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
Введем следующие обозначения: |
|
|
|
|
|
|
|
|
||||||||
|
x |
Qвз |
; |
x |
в |
|
Qвз |
|
в |
; |
dx |
Qвз |
d . |
(5.57) |
||
|
|
|||||||||||||||
|
bhm |
bhm |
bhm |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Тогда, подставляя (5.57) в (5.56), получим |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
в |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s d sсв в |
1 . |
|
|
(5.58) |
|||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Поскольку f |
s ds , то |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d f s ds . |
|
|
|
|
|
||||
Следовательно, из (5.58) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
sв |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sf s ds 1 sсвf sсв . |
|
|
(5.59) |
||||||||
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
В выражении (5.59) |
принято, что при |
x 0 |
и 0 , |
т.е. на входе в пласт, |
||||||||||||
мгновенно устанавливается водонасыщенность |
s , при |
которой kн 0 , а на |
||||||||||||||
фронте вытеснения значение ее в течение всего процесса составит sв . |
|
|||||||||||||||
Выполним интегрирование в левой части (6.59) по частям. Имеем |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sв |
sв |
|
|
sв |
|
|
sвf sв s f s f sв f s |
|
|
|||||||
sf s ds |
|
|
sf s |
f s ds |
|
(5.60) |
||||||||||
s |
s |
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
В соответствии |
со |
сказанным |
водонасыщенность |
s |
устанавливается в |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сечении 0. Следовательно, f |
s 0, поэтому и второй член в формуле (5.60) |
|||||||||||||||
равен нулю. Далее, поскольку kн s 0 , то, согласно формуле (5.47), |
f s 1 . |
Таким образом, из (5.59) и (5.60) получим:
sвf sв f sв sсвf sв ,
110

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
откуда
f sв |
f sв |
. |
(5.61) |
|
sв sсв |
||||
|
|
|
На рис.47 приведен график, построенный с учетом кривых относительных проницаемостей, данных на рис. 40, при в н 0.5.
По кривой f s можно найти значение sв , графическим путем. В самом деле,
согласно рис.47
f sв tg |
f sв |
. |
|
sв sсв |
|||
|
|
Проведя касательную к кривой f s из точки s sсв , по точке касания (см.
рис.47) определяем f sв и sв .
Для того же, чтобы найти распределение водонасыщенности по длине пласта,
необходимо построить кривую f s (рис.48). Это можно сделать методом графического дифференцирования кривой f s или, представив кривые относительных проницаемостей аналитически, выполнить дифференцирование аналитическим путем, сделав соответствующее построение.
Определим теперь длительность безводного периода добычи нефти, т.е.
момент времени t t , когда фронт вытеснения достигнет конца пласта и,
следовательно, xв будет равен l .
Рис.47. График зависимости f s от S |
|
s от S |
Рис.48. График зависимости f |
111

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
Будем считать, что к этому моменту времени в пласт закачано Qвз Q t
воды. Имеем из (5.57)
|
|
bhml |
|
f sв . |
|
|
|
|
(5.62) |
|||||
|
|
Q t |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
||||
Из (5.62) определим Q t и, |
следовательно, t . |
|
|
Величина |
bhml равна |
|||||||||
объему Vп пор пласта. Так как режим |
жесткий водонапорный, объем закачанной |
|||||||||||||
в пласт воды к моменту времени t t |
равен объему добытой из пласта нефти |
|||||||||||||
Qн к этому же моменту времени, |
т.е. Q t Qн . |
Безводная нефтеотдача |
||||||||||||
0 01 2 , где 01- коэффициент |
вытеснения |
нефти |
водой, достигнутый в |
|||||||||||
безводный период. Поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
Qн 2 |
|
|
|
2 |
|
|
. |
(5.63) |
||||
Vп 1 |
sсв |
f sв 1 sсв |
||||||||||||
|
|
|
|
Заметим, что распределение водонасыщенности в пласте изменяется по мере
продвижения в глубь пласта фронта вытеснения нефти водой таким образом, что
значения sв на фронте вытеснения xв и s на входе в пласт остаются неизменными. Таким образом, кривая распределения водонасыщенности как бы
«растягивается», оставаясь подобной себе. Такое распределение некоторого параметра, будь то водонасыщенность или какой-либо другой параметр,
называется автомодельным. Соответствующие решения задач также именуются автомодельными.
Полученные формулы позволяют рассчитать распределение
водонасыщенности к моменту подхода воды к линии добывающих скважин, т.е. в
безводный период разработки пласта.
Однако добыча нефти из пласта продолжается и после прорыва фронта вытеснения к концу пласта при x l .
Для определения текущей нефтеотдачи и обводненности продукции при t t , т.е. в водный период разработки пласта поступим следующим образом.
Будем считать, что продвижение фронта вытеснения происходит и в водный период разработки пласта, но этот фронт распространяется вправо за пределы пласта (рис. 49).
112

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
Рис.49. Схема вытеснения нефти водой из прямолинейного пласта в водный период разработки. Распределение водонасыщенности: 1-истинное; 2-фиктивное
Водонасыщенность на таком фиктивном фронте вытеснения и в этом случае остается постоянной, равной sв , а водонасыщенность при x l уже составит s .
Пусть в некоторый момент времени t t фиктивный фронт находится на расстоянии xвф от входа в пласт (см. рис.49). В соответствии с формулами (5.54)
и (5.55) при t t можно написать
bhml f s .
Qвз t
Из (5.62) и (5.64) получим
f s |
|
Q t |
|
|||
|
|
|
|
|
. |
|
f sв |
Qвз t |
|||||
|
|
(5.64)
(5.65)
|
По формуле (5.65) находим s для различных значений времени t . Так, зная |
||||||||
Q t , f sв и Qвз t , определим вначале |
f s , а затем по графику функции |
||||||||
|
s - значение s . |
|
|
|
|
|
|
|
|
f |
|
|
|
|
|
|
|
|
|
|
Дебиты нефти и воды в водный период разработки пласта составят |
|
|||||||
|
qн |
bhkkн s |
p |
qв |
bhkkв s |
p |
|
||
|
|
|
; |
|
|
. |
(5.66) |
||
|
н |
в |
|||||||
|
|
|
x x l |
|
|
x x l |
|
Отсюда для определения текущей обводненности продукции получим
формулу
|
qв |
|
kв s |
|
|
. |
(5.67) |
|
qв qн |
kв s |
в |
kн s |
|||||
|
|
|
|
н |
|
|
113

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
Текущую нефтеотдачу в водный период разработки пласта можно определить в принципе следующим образом:
1) установлением объема накопленной добычи нефти по формуле
t
Qн qн t dt ;
0
2) отнесением этого объема накопленной добычи нефти к первоначальному объему нефти в пласте, равному bhm 1 sсв .
Однако во втором случае можно определять объем добытой из пласта нефти по изменению в нем водонасыщенности, учитывая опять-таки то, что режим разработки пласта жесткий водонапорный. Так, на основе равенства объема вошедшей в пласт воды объему вытесненной из него нефти имеем
|
н |
|
|
|
l |
s x dx s |
св |
|
|
|
|
l |
|
|
bhmx |
|
св |
|
bhml |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Q |
|
bhm |
|
|
|
l |
|
|
qt |
|
|
sd |
|
|
qt |
|
s |
|
|
|
qt |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
bhml |
s |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
qt |
sd sсв f |
|
s |
|
|
sf |
|
|
|
|
|
|
|
|
(5.68) |
||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
f s |
|
|
s ds sсв f |
s |
||||||||||||||||||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
bhml |
sf s |
s f s f s f s sсв f s |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
f s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Формула (5.68) должна быть справедлива для всех моментов времени, когда |
||||||||||||||||||||||||||||||||
t t . При t , вообще говоря, |
водонасыщенность должна стать равной s во |
всем пласте. Однако при любом другом значении времени водонасыщенность
s s |
только на входе в пласт, |
т.е. при |
0 . Тогда, |
как следует из формулы |
|||||||||||||||
(5.67), |
|
. Следовательно, из (5.68) получим |
|
|
|
|
|
||||||||||||
f s 0 |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
1 f |
s |
|
|
|
|
|
||||
|
|
Qн |
s sсв |
|
|
|
|
|
|
. |
|
|
|
|
(5.69) |
||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
f s |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из (5.69) вытекает, что текущая нефтеотдача пласта в период водной его |
|||||||||||||||||||
эксплуатации |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 f s |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
s |
sсв |
|
|
2 |
|
|
|
|
|
|
|
Qн 2 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
f s |
. |
(5.70) |
||||||
|
|
V |
|
1 s |
св |
|
|
|
|
|
1 s |
св |
|
||||||
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
114

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
Таким образом, мы определили основные технологические указатели разработки элемента пласта - текущую нефтеотдачу и обводненность добываемой продукции.
Рассмотрим непоршневое вытеснение нефти водой в радиальном направлении, например, при разработке элемента семиточечной системы с использованием заводнения. Схема элементарного объема пласта для такого случая показана на рис.50. Уравнение неразрывности фильтрующейся воды в таком объеме получим с учетом баланса втекающей и вытекающей воды за время dt в виде
|
vв |
2 rd hvвdt 2 r dr d h vв |
r |
|
|
(5.71) |
dr dt 2 rdrd mds 0. |
|
|
|
Рис.50. Схема элементарного объема радиального пласта
Раскрывая скобки в выражении (5.71), сокращая в нем соответствующие члены и заменяя обозначения обыкновенных производных на частные, имеем
vв |
vв |
m s |
0 |
|
|||||
|
|
||||||||
|
r |
|
r |
|
|
t |
|
|
|
или |
|
|
|
|
|
|
|
||
|
1 vвr |
|
|
s |
|
|
|||
|
|
|
r |
m |
t 0 . |
(5.72) |
|||
|
r |
Вполне аналогичным образом, но с учетом того, что насыщенность пористой среды нефтью sн 1 s , установим соответствующее уравнение неразрывности для фильтрующейся в пласте нефти в следующем виде:
1 vнr |
m |
s |
0. |
(5.73) |
|
|
|
|
|||
r r |
t |
Складывая уравнения (5.72) и (5.73), получим
115

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
v vн vв |
q t |
. |
(5.74) |
|
|||
|
2 rh |
|
Вводя, как и в случае прямолинейного вытеснения нефти водой, функцию f s , определяемую формулой (5.47) (Бакли-Леверетта), и подставляя ее в (5.72) с
учетом (5.74), будем иметь одно дифференциальное уравнение для определения водонасыщенности s в виде
|
s |
|
|
s s |
|
|
||
m |
|
q t f |
0 . |
(5.75) |
||||
t |
2 rh |
|
r |
|||||
|
|
|
|
|
Так же, как и в прямолинейном случае, рассматриваем временем в пласте линий s const . В этом случае
ds rs dr st dt 0 .
Из (5.75) и (5.76)
dr q t f s . dt 2 rhm
Отсюда
f s mhkr2
Qвз
t
Qвз q t dt .
0
перемещение со
(5.76)
(5.77)
Рассмотрим баланс закачанной в пласт и извлеченной из него воды.
Устремляя для простоты радиус скважины к нулю ( rc 0 ), имеем
rв |
2 hmsrdr msсвhrв2 |
Qвз . |
|
|
|||||
|
|
|
(5.78) |
||||||
|
0 |
|
|
|
|
|
|
|
|
Учитывая, что |
|
|
|
|
|
|
|
|
|
f s ds |
|
2 mhrdr |
|
|
sв |
mhr2 |
|||
|
|
|
; |
f |
|
в |
, |
||
|
|
Qвз |
|
||||||
|
|
|
|
|
|
|
Qвз |
|
и подставляя эти выражения в (5.78), приходим к интегральному соотношению
sв
f s ds 1 sсвf sв ,
s
116

РАЗРАБОТКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
в точности совпадающему с соответствующим соотношением (5.59) для случая вытеснения нефти водой из прямолинейного пласта. Можно поэтому утверждать,
что и при вытеснении нефти водой из радиального пласта справедливы соотношение (5.60) и все последующие рассуждения, включая формулу (5.61),
пригодную для нахождения водонасыщенности на фронте вытеснения нефти водой, а также описанный графический метод определения sв .
Время t безводной разработки пласта радиусом rк определим из (5.77).
Если полагать, что Qвз qt , имеем
|
hr2m |
|
|
|
t |
k |
. |
(5.79) |
|
q |
||||
|
|
|
||
Аналогично по формулам (5.66) и (5.67) находим текущую обводненность |
||||
продукции, добываемой из пласта при |
t t . Соответственно |
текущую |
нефтеотдачу вычислим по формуле (5.70). Таким образом, определяем все важнейшие технологические показатели процесса вытеснения нефти водой.
5.4. РАСЧЕТ ПЛАСТОВОГО ДАВЛЕНИЯ И ДЕБИТОВ СКВАЖИН При определении забойного давления в скважинах с целью выбора способов
подъема жидкости из глубины на дневную поверхность, оценки фазового состояния нефти и воды, а также для вычисления градиентов пластового давления с целью определения скоростей перемещения фильтрующихся веществ, границ разделов между нефтью и водой необходимо знать поле пластового давления.
При решении задач фильтрации неоднородных жидкостей, в частности нефти и воды, наряду с вычислением поля водонасыщенности определяют и поле пластового давления. В случае вытеснения нефти водой из прямолинейного или радиального пласта при использовании модели поршневого вытеснения поле давления вычисляется просто по формулам, приведенным в предыдущем разделе.
В случае непоршневого вытеснения нефти водой даже из прямолинейного пласта распределение давления в нем устанавливать несколько сложнее. Поэтому рассмотрим последний случай более подробно.
117