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

1209

.pdf
Скачиваний:
2
Добавлен:
15.11.2022
Размер:
11.79 Mб
Скачать

Рис. 77. Схема элемента пласта при непоршневом вы­ теснении нефти водой

ЬЬиъ

h

Т7777Т77Т7Т77Т7ТГП ТП

А х

Рассматривая непоршневое вытеснение нефти водой в пря­ молинейном пласте, выделим элемент длиной Ах, высотой h и

шириной

b в направлении,

перпендикулярном к

плоскости

(рис. 77).

В общем случае

слева в элемент пласта

поступают,

а справа вытекают нефть и вода. При этом расход воды слева равен bhvB, а справа — bh(vв-{-Д^-в Ах).

Количество

накопленной воды в элементе пласта состав-

д8

(v — скорость фильтрации воды; s — водонасы-

ляет Ь/гт— Ах

щенность пласта; t — время). Согласно закону сохранения мас­ сы вещества, разность между скоростями входящей в элемент пласта воды и выходящей из него равна скорости накопления объема воды в элементе пласта. Выражая сказанное в матема­ тической форме, получим

— bh fvB+

Ахj + bhvB= bhm -^- Д*.

 

После сокращения соответствующих членов при устремле­

нии Ах— >-0 имеем

 

dvB .

ds

л

(IV.43)

а Г + 'п Ж = 0 -

 

Поскольку в пористой среде содержатся только нефть и вода, то насыщенность пористой среды нефтью sH= l —s. Рас­ сматривая аналогично предыдущему скорости проникновения нефти в элемент пласта и выхода из него, получим

fog

ds_

 

(IV.44)

дх

т dt = 0.

 

 

Складывая уравнения (IV.43) и (IV.44), имеем

 

 

(^н + О = °;

»ii+0B= 0(O-

(IV.45)

Таким образом, суммарная скорость фильтрации нефти и воды не изменяется по координате х, что и следовало ожидать,, так как нефть и воду принимают за несжимаемые жидкости.

Следовательно, режим пласта жесткий водонапорный.

141

Скорости фильтрации воды и нефти подчиняются обобщен­ ному закону Дарси, так что

kkB (s)

dp

 

kk„(s)

др

 

(IV.46)

рв

дх

н

рн

дх

*

 

где kB и &н, цв и ||н — относительные проницаемости, зависящие от водонасыщенности s и вязкости воды и нефти.

Рассмотрим

функцию f(s),

называемую функцией Бакли —

Леверетта. При этом

 

 

 

/(*) =

ув

 

М$)

 

(IV.47)

УВ+ Ун

£ в (s) +

рв

(s)

 

 

Рн

 

 

 

 

 

 

или

 

 

 

 

 

/(«) =

Ув

 

 

 

(IV.48)

у (0 ‘

 

 

 

Из

(IV.48), дифференцируя

vB по х, получим

 

 

 

 

 

(IV.49)

После подстановки (IV.49) в (IV.43) получим одно диффе­ ренциальное уравнение первого порядка для определения s, т. е.

v(t)f'(S) - ^ - + m ^ - = 0.

 

(IV. 50)

По мере вытеснения

нефти водой из прямолинейного

пласта

фронт вытесняющей нефть воды продвигается к концу

пласта

и водонасыщенность в каждом сечении заводненной области

непрерывно

увеличивается. Процесс

вытеснения

нефти

водой

из прямолинейного пласта можно представить и

иным

обра­

зом, следя

за изменением по пласту

некоторой

водонасыщен­

ности. Если, например, в какой-то момент времени в некотором сечении пласта водонасыщенность составляла s=si, то спустя определенное время эта водонасыщенность будет и в конце пласта, так как нефть постепенно извлекается из него и ее место занимает вода. Для указанного s = const можно принять

ds = - ^ - dx-\--^-dt = 0.

или

ds

dx i

ds

(IV. 51)

дх

dt

dt

 

 

Сравним (IV.50) и (IV.51). Они будут идентичными, если

положить

 

 

dx___ Г (*) у (0

(IV. 52)

dt

 

т

 

 

142

Умножим и разделим (IV.52) на bh и проинтегрируем, полу­ чим

 

t

 

bhmx = f (s)Qm (t);

<?„(/) = f bhv (f) dt.

(IV.53)

 

О

 

Обозначим

 

 

bhmx

 

(IV. 54)

QBS(0 ’

 

 

 

тогда

 

 

l =

 

(IV. 55)

Задавая s в

формуле (IV.55), можно

определить расстоя­

ние от входа в пласт для данного значения водонасыщенности Однако в период безводной эксплуатации закачиваемая вода еще не достигает конца пласта. Чтобы установить положение фронта вытеснения нефти водой и водонасыщенность на фрон­ те вытеснения, рассмотрим материальный баланс закачанной в пласт воды. Если к моменту времени t в пласт закачан объем

воды, равной QB3(t),

длина фронта вытеснения

составит хв,

насыщенность пласта связанной водой s = sCB> то

 

<2вз (0 = bhm J s (х) dxbhmxBsCB.

(IV.56)

 

о

 

 

 

 

 

Введем следующие обозначения:

 

х

QВЗ t.

„ _

QB3 £ .

 

= bhm s’

 

 

bhm Sd’

 

 

 

 

(IV.57)

 

 

 

 

 

d x = bhm dl.

 

 

 

 

Тогда, подставляя (IV.57) в (IV.56), получим

 

J s ( 9 d S - s CBl, =

l.

 

 

(IV. 58)

О

 

 

 

 

 

Поскольку £,= f'(s),

то

 

dl = f" (s) ds.

 

 

 

 

Следовательно, из

(IV.58)

 

SB

 

 

 

 

 

\ s r ( s ) d s = l + s CBf' (s,).

(IV.59)

В выражении (IV.59) принято, что при х=0 и £= 0, т. е. на

143

Рис. 78.

График зависимости f(s)

от s Рис.

79.

График зависимости f'(s)

 

 

 

 

от s

 

 

 

 

входе

в

пласт,

мгновенно

устанавливается

водонасыщенность

s*, при

которой

kB=0 (см. рис. 40),

а

на

фронте

вытеснения

значение ее в течение всего процесса составит sB.

по частям.

Выполним интегрирование в левой части

(IV.59)

Имеем

 

 

 

 

 

 

 

 

о

 

В

В

 

 

 

 

 

J s/" (s) ds = s/' (s)— |' f' (s)ds =

 

 

 

 

 

 

S ,

S i

 

 

 

 

 

= SB/ ' (Sb)—S*f' (S*)—/ (SB)+

/ (S*)•

 

 

 

(IV. 60)

В соответствии со сказанным водонасыщенность s* устанавли­ вается в сечении £=0. Следовательно, /'(s*) = 0, поэтому и вто­ рой член в формуле (IV.60) равен нулю. Далее, поскольку &H(s*)=0, т о , согласно формуле (IV.47), f(s*) = l. Таким обра­ зом, из (IV.59) и (IV.60) получим

sBf' (sB)— f(sB) = scJ ' (sB), откуда

/'(*.) =

7 ^

- -

 

(IV. 61)

 

 

’ *св

 

 

На рис. 78 приведен график, построенный

с учетом

кривых

относительных проницаемостей, данных на

рис. 40, при

рв/цн =

= 0,5.

кривой f(s) можно найти значение

sB графическим пу­

По

тем. В самом деле, согласно рис. 78

 

 

п * в ) =

t g a

/Ев)

 

 

SB--- SQ[J

 

 

 

 

 

 

144

Проведя

касательную к кривой f(s) из точки s = sCB, по точ­

ке касания

(см. рис. 78) определяем f(sB) и sB.

Для того же, чтобы найти распределение водонасыщенности

по длине пласта, необходимо построить кривую f'(s) (рис. 79).

Это можно сделать методом графического дифференцирования кривой f(s) или, представив кривые относительных проницае­ мостей аналитически, выполнить дифференцирование аналити­ ческим путем, сделав соответствующее построение.

Определим теперь длительность безводного периода добычи

нефти, т. е. момент времени £ = £*, когда

фронт

вытеснения до­

стигнет конца пласта и, следовательно, хв будет равен I. Будем

считать, что к этому моменту времени в

пласт

закачано QBa =

= Q*(t*) воды. Имеем из (IV.57)

 

 

 

 

 

о т = / '& )•

 

 

 

 

 

 

(lv -62>

Из (IV.62) определим Q*(7*) и, следовательно,

 

Величина

bhml равна объему Vn пор пласта. Так

как

режим

жесткий

водонапорный,

объем

закачанной в

пласт воды

к

моменту

времени t = t*

равен

объему добытой

из

пласта

нефти QH* к

этому же моменту времени, т. е. Q*(t*) = QH*. Безводная нефте­

отдача rjo= Ло12,

где r]oi — коэффициент вытеснения

нефти во­

дой, достигнутый в безводный период. Поэтому

 

у , ______ _____________________ % ______

(IV.63)

l o - Vn ( l - s CB)

f'(sB) ( l - s CB)

 

Заметим, что

распределение водонасыщенности

в пласте

изменяется по мере продвижения в глубь пласта фронта вы­ теснения нефти водой таким образом, что значения sBна фрон­ те вытеснения хв и s* на входе в пласт остаются неизменными. Таким образом, кривая распределения водонасыщенности как бы «растягивается», оставаясь подобной себе. Такое распреде­

ление некоторого

параметра,

будь

то водонасыщенность или

какой-либо другой

параметр,

называется

а в т о м о д е л ь н ы м .

Соответствующие

решения задач

также

именуются автомо­

дельными.

Полученные формулы позволяют рассчитать распределение водонасыщенности к моменту подхода воды к линии добываю­ щих скважин, т. е. в безводный период разработки пласта.

Однако добыча нефти из пласта продолжается и после про­ рыва фронта вытеснения к концу пласта при х=1.

Для определения текущей нефтеотдачи и обводненности продукции при t> t*, т. е. в водный период разработки пласта, поступим следующим образом. Будем считать, что продвиже­ ние фронта вытеснения происходит и в водный период разра­ ботки пласта, но этот фронт распространяется вправо за пре­ делы пласта (рис. 80). Водонасыщенность на таком фиктивном фронте вытеснения и в этом случае остается постоянной, рав­

ной sB, а водонасыщенность при х = 1 уже составит 5. Пусть в

10 Ю. П. Желтов

145

1

 

 

 

Рис. 80. Схема вытеснения

 

 

S

2

нефти водой из прямолиней­

----------------\

/

ного пласта в водный пери­

 

/

од разработки. Распределе­

S*

■-

- ‘— г

ние водонасыщенности:

*7777777777777777777777777777777777777777777^

 

__ __

1 — истинное; 2 — фиктивное

 

 

 

■Тцш

^

 

 

 

О

 

 

X

 

некоторый момент времени на расстоянии хВф от входа ствии с формулами (IV.54) сать

t > t * фиктивный фронт находится

впласт (см. рис. 80). В соответ­

и(IV.55) при t> t* можно напи­

bhml

= /

' (s).

(IV.64)

QBS(0

 

Из (IV.62)

и (IV.64) получим

 

Г (SB)

QD3 (0

(IV.65)

 

По формуле (IV.65) находим s для различных значений вре­ мени t. Так, зная Q*(t*), f'(sB) и QB3(X)> определим вначале /'(s), а затем по графику функции f'(s) — значение s.

Дебаты нефти и воды в водный период разработки пласта составят

_

 

(

д р \

(IV.66)

рв

\

дх Jx-l

 

Отсюда для определения текущей обводненности продукции v получим формулу

V —

(IV.67)

9ю+ 9н

*-<5> + 7 Г

Текущую нефтеотдачу в водный период разработки пласта мож­ но определить в принципе следующим образом:

1) установлением объема накопленной добычи нефти по фор­ муле

t

Q„= ffc (<)<";

2) отнесением этого объема накопленной добычи нефти к первоначальному объему нефти в пласте, равному bhm( 1—sCB).

Однако во втором случае можно определять объем добытой

из пласта нефти по изменению в нем водонасыщенности, учи­ тывая опять-таки то, что режим разработки пласта жесткий во­ донапорный. Так, на основе равенства объема вошедшей в пласт воды объему вытесненной из него нефти имеем

Qn= bhm[[s (х) dx— sCB/J = qt\^sd

—sCB

j =

о

0

 

 

5W

= т ш ~ [Isf" <s)

 

 

= qt |J sdl—sCB/' (s)

Щ

 

о

 

isf’ ( s ) - s j ' ( S j - f ( s ) + f ( s * ) - s j ' (s)].

(IV. 68)

Формула (IV.68)

должна быть справедлива для всех моментов

времени, когда

При t— »-оо, вообще говоря, водонасыщен-

ность должна стать равной s* во всем пласте. Однако при лю­ бом другом значении времени водонасыщенность s=s* только на входе в пласт, т. е. при £= 0. Тогда, как следует из формулы

(IV.55), /'(s*)= 0 .

Следовательно, из (IV.68) получим

 

<2н = ^п

1—т

(IV.69)

Г

(s)

 

 

Из (IV.69) вытекает, что текущая нефтеотдача пласта в период водной его эксплуатации

 

 

i - т

ла

 

S C B

Г (s)

*1 =

 

(IV.70)

SCB)

 

M l

 

 

Таким

образом, мы определили основные технологические

показатели разработки элемента пласта — текущую нефтеотда­ чу и обводненность добываемой продукции.

Рассмотрим непоршневое вытеснение нефти водой в радиаль­ ном направлении, например при разработке элемента семито­ чечной системы с использованием заводнения. Схема элементар­ ного объема пласта для такого случая показана на рис. 81. Уравнение неразрывности фильтрующейся воды в таком объеме получим с учетом баланса втекающей и вытекающей воды за

Рис. 81. Схема элементар­ ного объема радиального пласта

10*

147

время dt в виде

2nrdQhvBdt— (r-\-dr) dQh ^ в+ -^г- drj dt—

2nrdrdQmds = 0.

(IV.71

Раскрывая скобки в выражении (IV.71), сокращая в нем соот­ ветствующие члены и заменяя обозначения обыкновенных про­ изводных на частные, имеем

dvя

ds Л

 

дг

■тж = 0

 

 

 

ИЛИ

 

 

- г Ч

г - + т Ж = 0-

(IV-72)

Вполне аналогичным образом, но с учетом того, что насыщен­ ность пористой среды нефтью sH= l —s, установим соответствую­ щее уравнение неразрывности для фильтрующейся в пласте нефти в следующем виде:

г

дМ

ds

п

(IV.73)

дг

m w = ° -

 

Складывая уравнения (IV.72) и (IV.73), получим

 

0 = 0"+

t’B = . |l r -

 

(IV-74)

Вводя, как и в случае прямолинейного вытеснения нефти водой, функцию f(s), определяемую формулой (IV.47), и подставляя ее в (IV.72) с учетом (IV.74), будем иметь одно дифференци­ альное уравнение для определения водонасыщенности s в виде

ds

qV (s)

ds

n

(IV.75)

т dt

2nrh

dr

 

Так же, как и в прямолинейном случае, рассматриваем переме­

щение со временем в пласте линий s = const.

В этом случае

ds= - ^ dr+

^ d t = 0 -

(l v -76>

Из

(IV.75) и

(IV.76)

 

dr

_ q(l) f' (s)

 

dt

2nrhm

 

 

Отсюда

mnkr3

(IV.77)

~QZ~

 

Qm = \q (t)d t.

 

148

Рассмотрим баланс закачанной в пласт и извлеченной из него воды. Устремляя для простоты радиус скважины к нулю (чс— >-0), имеем

j hmsrdrnmsCBhrB2= QB3.

 

(IV.78)

о

 

 

 

 

 

 

Учитывая из

(IV.71), что

 

 

f"(s)ds

2 n m h r d r

Г (SB)

nmhrB2

 

Q.B3

Q B3 *

 

 

 

 

и подставляя эти выражения в (IV.78), приходим

к интеграль­

ному соотношению

 

 

SB

 

 

 

 

 

 

j >

( s ) d s

= l

+ S cl/ ( S B ) .

 

 

S*

 

 

 

 

 

 

в

точности

совпадающему

с соответствующим

соотношением

(IV.59) для случая вытеснения нефти водой из прямолинейного пласта. Можно поэтому утверждать, что и при вытеснении неф­ ти водой из радиального пласта справедливы соотношение (IV.60) и все последующие рассуждения, включая формулу (IV.61), пригодную для нахождения водонасыщенности на фрон­ те вытеснения нефти водой, а также описанный графический метод определения sB.

Время t* безводной разработки пласта радиусом гк опреде­

лим из

(IV.77). Если полагать, что Q*z = qt, имеем

 

лhrK2m

(IV.79)

t *

Я

 

 

Аналогично по формулам (IV.66) и (IV.67) находим теку­ щую обводненность v продукции, добываемой из пласта при /> /* . Соответственно текущую нефтеотдачу г| вычислим по фор­ муле (IV.70). Таким образом определяем все важнейшие тех­ нологические показатели процесса вытеснения нефти водой.

П р и м е р IV.2. Из элемента однорядной системы разработки осущест­ вляется вытеснение нефти водой. Будем условно принимать процесс движения нефти и воды во всем элементе пласта одномерным, прямолинейным.

Длина пласта /=600 м, его ширина 6= 600 м, общая эффективная толщи­

на Л0= 20 м. Коэффициент охвата пласта

заводнением по толщине г)2= 0,75,

так что охваченная заводнением толщина

пласта составляет Л= 15 м. Прони­

цаемость пласта 6=0,5 мкм2, пористость т = 0,2. Вязкость нефти в пластовых условиях рн = 4-10-3 Па-с, вязкость воды |лв=10-3 Па-с. Через границу пласта при х = 0 закачивается вода с расходом q= 200 м3/сут и столько же жидкости добывается с конца элемента при х=1 в течение всего рассматриваемого периода разработки.

149

Относительные проницаемости заданы следующим образом:

Ли (s) = ( 1 7 = ^ " ) ПРИ «св < « < s*;

К (S) ms( - £ Е г £ )

ПР"

5св ^

S ^ ®х;

*в(s) — 0 ,8 (

1

ПРИ

51 < « ^ V

При этом 5св“ 0,1; s*=0,8. Значение s 1 определяем из условия равенства

относительных проницаемостей для нефти и воды при s=s\.

Требуется рассчитать на основе изложенной теории изменение во времени текущей нефтеотдачи и зависимость обводненности продукции v от текущей

нефтеотдачи ц.

Приступая к решению задачи, определим прежде всего эь Имеем

ljL = h a A='

0 8

——

\i/«

/Д =*а_\:

— «е. )

 

^

/

Отсюда

 

 

 

«1 —«ев \Va

0 ,8;

 

( s'.-t* . У = ° '64; s ,- s c .= (s.-%.)0.641/*:

Таким образом, при sc*<s<0,7032 функция Бакли — Леверетта

м«)

(

V

/(*>=

 

(«—«ед)а

 

” Cs— s»)» + « ,2 5 ( s . - s ) ' -

 

При 0,7©32«s«0„8

 

« — «СВ У

« * — «СВ /

т

Фуакщя Ms), построенная по приведенным формулам, представлена на графике (рпе. Щ , Определим по формуле (1\Ш ) водонаеыщенность па фрон­ те вытеснения, для чего проведем касательную к кривой f(s) из точки s=Sc*. Мз рис. Ш подучаем, чпго $®=0,413; Ш®) =0,723. Теперь необходимо построить кривую f (s). Так как формулы для }(§) в рассматриваемом примере срашштеаьно простые, функпвю f (s) можно получить путем обычного, а не графиче­

ского дщфферевдпроваппя функвдию Ц$).

Прн имеем

ISO

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