Рис. 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 (х) dx—bhmxBsCB. |
(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= Ло1'П2, |
где 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 |