1209
.pdfРис. 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 |
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 — |
9а |
(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— 2л (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 |
q(О V (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 2лhmsrdr—nmsCBhrB2= 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