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

1275

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

s,c

Рис. 77 График распределения на­ сыщенности водной фазой и концен­ трации активной примеси в пласте в случае значительной сорбции ак­ тивной примеси породой (пунктиром показано распределение водонасыщенности при обычном заводнении)

вытеснения нефти водой Яф, т. е. 5+>5ф, ггде Яф определяется так же, как и по теории Баклея—Леверетта:

5ф=«0+ 4 Т Г Г ’

(6-38)

то решение для зоны II аналогично решению, данному Баклеем—

Левереттом, а именно:

 

Z(s, т) = / ' (s, 0) т

(6.39)

для любых s, удовлетворяющих неравенству

 

5ф < s < S+

 

Для этого случая характерное распределение насыщенности водной фазой s и концентрации примеси с будет таким, как пока­ зано на рис. 77.

Если я+СЯф, по пласту перед скачком концентрации будет распространяться «ступенька» насыщенности водой при s = s+ (в противном случае в пласте имелись бы точки с неоднозначным значением насыщенности). Найдем скорость продвижения фронта вытеснения нефти водой в этом случае, для чего рассмотрим урав­ нение материального баланса воды в зоне II. Рассмотрим моменты времени, при которых фронт вытеснения еще не дошел до линии отбора.

За время t в зону II войдет объем воды

/ *С «> t * с <*>

Ув = f v{t)dt— m b o

I (s—s0) dx = \v (t) dt + ms0xc {t)— m

j sdx,

o

о

который займет поровый объем зоны II, равный, очевидно, VB, т. е. *ф W

Ув = т j (s—sQ)dx = m{s+— s0) [*ф(0 xc(t)].

xc (i)

Вводя безразмерные переменные и приравнивая оба эти выра­ жения, получим

£с <*)

(6.40)

т+*!.(т)— f sdl = (s+ —So) [1ф(x)—1c(x)].

6

 

8*

231

Вычислим интеграл, входящий в это выражение. Заменяя пе­ ременную интегрирования £ на s и используя формулу (6.37), получим

d l = f" (s, с°) Tds.

При этом, очевидно, пределы интегрирования изменятся сле­

дующим образом: нижний £ =

0 будет иметь вид s = s°, а верхний

i = £с 00 — s =

sr. Тогда

 

Бс <т>

S-

 

| sdh, =

т

( sf" (s, с°) ds.

 

6

 

 

 

Вычисляя

последний интеграл по частям, получим

S ~

 

 

 

J sf" (s, с°) ds = т [s- /' (s_, c°) —s°f' (s°, c°) —f(s~, c°) + f (s°, c0)] =

 

 

 

= T [1 + S-/'(S-,

c°)— f(s-, c0)],

 

где / (s°, c°) =

1; f' (s°, c°) =

0.

Учитывая, что из формулы (6.37) следует

Ес(т) = /'(s"

с°)т,

(6.41)

и подставляя найденное выражение для интеграла в уравнение (6.40), получим

(s^ —s°) Еф (т) —т [f (s—, c°) — (s-— s+)f/ (s-, с0)].

(6.42)

Упростим это соотношение.

 

Из теории фильтрации известно, что на границах области, где насыщенность изменяется скачкообразно, должно соблюдаться следующее:

f(s~,

c°) -f(s+, 0) ,

(6.43)

dx

 

S~ — S+

 

 

/(s-

c°) +

K

 

1- К

(6.44)

Vc

a + К

s

 

1-

К

 

 

 

Здесь К и а — константы растворимости соответственно активной примеси в нефти и адсорбции примеси породой.

Дифференцируя уравнение (6.41) по т, найдем еще одно уравне­

ние для определения безразмерной скорости

vc:

4 с

= vc

с°)

(6.45)

dx

Из (6.43) и (6.45) следует:

 

 

fis-,

c°)-/(s+, 0)

 

 

с°)

s—— s+

 

 

 

 

232

Рис. 78. График распределения на­ сыщенности водной фазой и концен­ трации активной примеси в пласте в случае слабой сорбции примеси

Рис. 79. Зависимость доли водного раствора в ПАВ в потоке от водонасыщенности s при слабой абсорбции.

Концентрация ПАВ:

I — с = 0; 2 с = &'

Подставляя это соотношение в уравнение (6.42), окончательно

получим

 

(s+—s0) M

T) = / ( S+, 0)т,

откуда скорость фронта вытеснения нефти водой

dx

(6.46)

s+ — s„

График характерного распределения насыщенности и концентра­ ции для этого случая показан на рис. 78.

Для определения ос, s~ и s+ имеем три уравнения (6.43)— (6.45), которые удобно записать в виде, позволяющем решить каж­

дое из них

независимо

от других:

 

 

 

 

/( s - , с°) +

 

(6.47)

/ ' ( S - , С»)

=

---------------------- b l T ~

 

 

 

s - + а+К

 

 

 

 

1

К

 

(6.48)

vc=f'(s~,

 

с°);

 

 

/ ( S - , C ° ) - / ( S + , 0)

 

 

(6.49)

 

 

С/с.

 

S—— S +

 

 

 

Как следует из этих уравнений, из первого определяем s~, из

второго — ос,

из третьего — значение

насыщенности s+ за

скач­

ком концентрации.

 

решать, применяя

графо­

Уравнения

(6.47) — (6.49) удобно

аналитическую технику. На диаграмме s, f (s), называемой фазовой,

построим функции f (s, с°) и f (s, 0)

(см. рис. 79) и проведем каса-

тельную из точки с координатами

f

а 4- К

К

\

к кривой

\

—------ — , — —---- —

/

 

1 — А

1 — Л

 

233

f (s, c°). На абсциссе в точке касания получим величину s~, а в точке пересечения касательной с кривой / (s, 0) найдем s+. Далее, оче­ видно, по тангенсу угла наклона касательной к оси абсцисс опреде­ лим скорость vc.

Две разобранные возможные ситуации (s+;>S(}, и s+ < 5ф) со­ ответствуют случаям сильной и слабой адсорбции закачиваемой активной примеси пористой средой. В обоих случаях фронту вы­ теснения нефти примесью предшествует фронт вытеснения нефти водой. Однако при значительной адсорбции фронт вытеснения нефти водой и передняя часть профиля водонасыщенности пол­ ностью совпадают с фронтом вытеснения нефти и профилем водонасыщенности при обычном заводнении (показано пунктиром на рис. 77).

Таким образом, сильно сорбирующаяся примесь не изменяет момента начала обводнения добываемой продукции. Прирост нефте­ отдачи получают лишь на поздней стадии разработки (к моменту подхода фронта концентрации к линии отбора). При закачке слабо сорбирующейся примеси (s+<s<jj) происходит отставание фронта вытеснения нефти водой по сравнению с обычным заводнением, т. е. увеличивается безводный период добычи, а прирост нефтеот­ дачи получаем уже на ранней стадии добычи нефти.

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

Для этого построим функцию доли водного раствора ПАВ с кон­ центрацией / (s, с°) и доли воды / (s, 0) в потоке фильтрующихся жидкостей в зависимости от насыщенности водной фазой s, исполь­ зуя заданные соотношения для относительных фазовых проницае­ мостей.

Построим вначале зависимость / (s, 0), придавая водонасыщенности s значения от s„ до 1 — s„ ост в равноотстоящих точках с ша­ гом

А___ (1 so 5н ост)

где М — общее число точек, для которых следует определить иско­ мую функцию / (s, 0). Аналогично находим функцию / (s, с°). Дан­ ные расчетов приведены в табл. 53 и на рис. 79.

Далее, воспользовавшись методикой, описанной выше, опреде­ ляем водонасыщенности s~, s+ и s$.

Учитывая, что по условию задачи растворимость ПАВ в нефти отсутствует, и полагая в уравнении для нахождения s- значение

К = 0,

получим

 

И * ".

С*)

I ( У . С») _

s- +

а

 

 

Проводя касательную к кривой / (s, с°) из точки с координатами (— а, 0), находим значение s~ (см. рис. 79):

s~ = 0,625.

234

Таблица 53

 

 

 

 

 

 

 

 

f is, 0)

 

f (s, c°)

 

f (s, 0)

 

f (s. d>)

0,200

0,000

0,200

0,0

0,463

0,707

0,611

0,840

0,226

0,006

0,241

0,013

0,489

0,787

0,652

0,889

0,253

0,026

0,282

0,055

0,516

0,852

0,693

0,926

0,279

0,064

0,323

0,130

0,542

0,902

0,735

0,952

0,305

0,122

0,364

0,232

0,568

0,939

0,776

0,971

0,331

0,199

0,406

0,351

0,595

0,965

0,818

0,984

0,358

0,294

0,447

0,475

0,621

0,982

0,858

0,992

0,384

0,399

0,488

0,591

0,647

0,993

0,899

0,997

0,411

0,508

0,529

0,692

0,674

0,998

0,940

0,999

0,437

0,613

0,570

0,774

0,700

1,000

0,981

1,000

По рис. 79 находим абсциссу точки пересечения построенной касательной с кривой f (s, 0), после чего определим значение s+:

s+ = 0,440.

Графическим построением находим точку Яф пересечения ка­ сательной к кривой f (s, 0), проведенной из точки с координатами, (s0, 0):

s.p= 0,475.

 

 

 

 

Отметим, что

sф есть не что иное,

как

корень уравнения (6.38 )

найденный графическим способом.

что

s+< ^ .

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

Из

сравнения

Бф и s+ заключаем,

фронт

вытеснения нефти будет отставать от фронта вытеснения

нефти

при обычном заводнении (случай малой

сорбируемости

ПАВ

пористой средой). Характерный

график распределения насы­

щенности водной фазой и концентрации ПАВ приведен на рис. 78. Перейдем к расчету нефтеотдачи. Для этого разобьем весь срок разработки залежи на два этапа, соответствующих безводному

периоду и периоду добычи обводненной продукции.

Очевидно, что в безводный период добычи нефти вследствие предположения о несжимаемости жидкостей объем добытой нефти равен объему закачанного в пласт раствора ПАВ, т. е.

Ун — Vзак в ~

Из формулы для определения коэффициента нефтеотдачи сле­

дует,

что

 

 

qt

Т1

= ----------—

= -------------------" --------------------

 

Vзап

tnbhl (\ — SQ)

или в безразмерных переменных

т

(6.50)

 

где т — безразмерное время; т]2т — коэффициент охвата пласта воздействием.

235

Определим момент окончания безводного периода нефтеотдачи. Скорость движения фронта вытеснения нефти дается выражением (6.46), интегрируя которое по т, получаем

?Ф(^) =

/(s+. 0) т<

S+ — So

Очевидно, безводный период заканчивается в момент подхода фронта вытеснения к добывающей галерее (линии отбора).

Полагая хф (t) = / (Нф (т) = 1 в безразмерных переменных), получим

 

S+ — %

,

 

(6.51)

 

/(*+. 0)

'

 

 

 

 

Определяя значение / (s+, 0) с помощью рис. 79, находим

 

т*

0,440 — 0,200

0,372

 

0,645

 

 

 

 

или в размерных величинах тЫг01 0,2-200-10-400 0,372 = 297,6 сут.

q

200

Вычислим коэффициент нефтеотдачи на конец безводного пе­ риода добычи нефти, подставив значение т* в формулу (6.50):

_Vbr_

0,372-0,8

= 0,372.

1 — %

0,8

 

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

У доб а (0 —^Ш в(0—У доб в (0-

Деля обе части этого уравнения на геологические запасы нефти получим текущий коэффициент нефтеотдачи на момент времени /

,, _

^зак в Щ —^даб в (О

 

1 ~

у

 

 

гзап

 

Далее очевидно, что

 

 

i

t

 

I)d ! =

i qf{s*. 0 ) d l= q f( s + , 0)(<— /.).

так как добыча воды начинается с момента времени t = /*. Таким образом,

4 t — 4f { s b , 0 ) ( f - U

rnsSsU(1 —

или в безразмерных величинах

„ _

II

 

0) II „

1 ~

к '

" и - й *

' ь ‘

ш ,

Подставляя выражение для т* (6.51), окончательно находим

Л =

T l ‘ - f ( s + ' 0) l + s + ~ s " n . T .

(6 . 52)

 

1 — S0

 

Из этого соотношения следует, что нефтеотдача в водный период при незначительной сорбируемости ПАВ нарастает по линейному закону.

Найдем время окончания разработки залежи, которое по ус­ ловию задачи заканчивается в момент подхода фронта ПАВ к ли­ нии отбора. Из уравнения (6.41) следует, что положение фронта ПАВ

в момент времени т определяется следующим образом:

 

ic (T )= r:(s-, с°)т.

 

 

 

Полагая в

этом

соотношении

 

хс (/) = I (|с (т) = 1),

получаем

Т * = ------'------.

 

 

 

 

 

/' (S-. с")

 

 

 

 

Из рис. 79 следует, что

 

 

 

/ ' (s-,

с°)

- ---- ,

 

 

 

 

где sm — координата на оси абсцисс точки пересечения

касатель­

ной к кривой / (s, с°) с прямой / =

1, откуда

 

т;*=5т + а-

 

 

 

(6.53)

Подставляя в это соотношение найденное из рис. 79 значение

s,n — 0,740 и а = 0,1, получим

 

 

 

т* =0,7404-0,1 =0,840

 

 

 

или в размерных величинах

 

 

 

=

mbh0l т.3 _

0,2-200-10-400

0,840 = 672 сут.

 

 

 

 

200

 

 

 

 

 

 

 

 

Конечный коэффициент нефтеотдачи можно определить по фор­

муле (6.52),

подставив значение т = т*:

 

,* __

[1 — f (s+,

0)] т* -L s+ s0

Лат-

(6.54)

ту5

 

1— s0

 

 

 

 

 

Это выражение можно упростить. Из рис. 79 следует, что

1—/(s+, 0)______ 1

 

 

 

s,n — s +

Sm - f а

 

 

 

Используя это соотношение и (6.53), находим

 

[\— f{s+,

0)]т* = sm— s+

 

 

 

После подстановки этого выражения в формулу (6.54) оконча-

тель

 

им

 

 

 

 

Т|*

Sm

Sp

 

 

 

(6.55)

1 Sos0 Л2т

 

 

 

 

 

 

237

Рис. 80. График зависимости теку­ щей нефтеотдачи от количества прокачанного раствора ПАВ в долях порового объема:

I — период безводной добычи; 2 — период обводненной продукции

Подставляя в (6.55) исход­ ные данные, находим конечный коэффициент нефтеотдачи

(0,740 — 0,2)0,8

Г)1' = —--------------------

0,8

-- 0,540.

Наиболее удобно результаты вычисления коэффициента нефте отдачи представить в виде графика зависимости его от количества

прокачанного раствора ПАВ т = тт]2тРезультаты таких расчетов по данным примера приведены на рис. 80.

Определим общий объем добытой нефти Vn (t*), Из формулы для определения коэффициента нефтеотдачи следует, что

Ун (^ * )= У за п 'П * (**)•

Подставляя исходные данные, получаем

У„ {t*) = mbhl (1 —s0) rf = 0,2 • 200 ■12,5 • 400-0,8 • 0,540 = 86400 м3.

По распределению насыщенности (см. рис. 78) обводненность добываемой продукции будет постоянной и равной f (s+, 0), т е

v3 = /(s+, 0)-0,645.

Таким образом, определены все технологические показатели разработки полосообразной залежи при закачке водного раствора ПАВ, которые требовалось найти по условию задачи.

З а д а ч а

6.9К. Из круговой залежи радиусом гк = 200 м,

толщиной h =

12,5 м, имеющей начальную водонасыщенность s0,

равную 0,2, и

пористость т = 0,2, вытесняется нефть вязкостью

р„ = 10 мПа-с при закачке в нагнетательную скважину, располо­ женную в центре залежи, водного раствора полиакриламида (ПАА) с концентрацией с° = 0,001 и расходом q = 250 м3/сут.

Вязкость водной фазы при закачке раствора ПАА в зависимости от его концентрации можно определить следующим образом:

рВф = M l + 0,93 • 103с + 0,74 • 106с2),

где рв — вязкость чистой воды (|iB= 1 мПа-с). При этом содержа­

ние

остаточной нефти

 

 

ост (^) = SHост.(1

100с)3,

где

S], ост — остаточная

нефтенасыщенность при вытеснении нефти

водой ( S O C T — 0,3)

238

ПАА не растворяется в нефти и сорбируется скелетом пористой среды по закону (изотерма сорбции Генри), формула которого имеет вид

а (с) = ас,

где а = 0,7.

Зависимости относительных фазовых проницаемостей нефти и воды (водной фазы) от насыщенности принять такими же, как и в задаче 6.8К.

Определить текущую и конечную нефтеотдачу при вытеснении нефти раствором ПАА в зависимости от прокачанных поровых объемов, время разработки залежи и общий объем добытой нефти. Процесс считать законченным к моменту подхода фронта концен­ трации ПАА к линии отбора, за которую принимаем круговую галерею, расположенную на расстоянии гк = 200 м от центра за­ лежи. Движение жидкостей считать плоскорадиальным. Коэффи­ циент охвата пласта по толщине воздействием г|от принять равным

0, 8.

Р е ш е н и е. Прежде чем переходить к решению поставленной задачи, выведем уравнения неразрывности для водной фазы и опре­ деления концентрации ПАВ при плоско-радиальном течении жидко­ стей.

Выделим элементарный объем пласта

AV = 2nh0rAr

ирассмотрим баланс объема водной фазы в элементе за время At. За это время количество водной фазы, поступившей в элемент

пласта, составило

 

qB\rAt=qf(s,

c)\rAt.

 

 

За то же время количество водной фазы, вытекающей из этого

элемента,

будет

 

 

 

 

йв|л + д г = qf [s,

с) |Л^.дгАГ

 

 

Очевидно, что разница между ними равна количеству водной

фазы, накопленной в элементе пласта за время At, т. е.

 

 

[qf{s,

c)\-r— qf{s,

c)\r + Ar]At = m AV[s{r, t + A t) — s{r, Q],

 

где

г — радиус

некоторой точки элемента AV пласта (т.

е.

г <

г

<

г + Аг).

 

 

Подставляя в последнее уравнение выражение для AV и деля

обе части на ArAt, получим

 

 

[qf (s,

с) 1г — qf (s,

с) \r+ л,]

 

 

 

 

Дл

 

 

 

Устремляя Аг и At к нулю и учитывая, что q (г) = q (г +

Лг)>

получим

 

 

 

 

 

m

ds

2 nh0r

(6 .5 6 )

 

 

dt

dr

 

23Э

Уравнение неразрывности для раствора ПАА получаем анало­ гично (6.56). Выделив тот -же элементарный объем А У, запишем балансовые соотношения.

За время At количество раствора ПАА, поступившего в элемент

объемом АУ,

составило

 

qBc\r At.

 

 

(Здесь

учтен

тот факт, что

ПАА не растворяется в нефти, т. е.

ср (с) =

0, и,

следовательно,

нефтяной фазой не переносится).

За то же время из элемента АУ количество вышедшего раствора ПАА составит

ЯъС|л + Д г At.

В момент времени t в элементе объема А У пласта находилось ПАА

mAV (sc + ac) |*.

В момент времени t + At

mAV (sc + ac)|,+A,.

Составляя уравнение баланса, получим

(двс\г — двС\г+Аг) At = mAV [(sc4 ас) |/+Д/—(sc + ac) |,]

или, деля

последнее соотношение

на ArAt и

устремляя Дг и Л/

к нулю,

получим

 

 

 

 

 

т ---- (sc-f ас)

h0r

дг

СЯъС) = 0 •

 

dt

 

 

 

 

Учитывая, что

qB = qf (s, с), окончательно

имеем

т “T r(sc + ac) + — -------7~[/(s,

с) с] —0.

(6.57)

dt

 

2пп0г

дг

 

 

 

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

показатели разработки залежи.

 

В соответствии с

условиями задачи сформулируем начальные

и граничные условия для уравнений (6.56), (6.57).

нефть

В начальный момент времени t = 0 в пласте находились

и связанная

вода. Следовательно,

 

s(r, 0) =

sCB= so;

с (г, 0) = 0-

(6.58)

В это же время через центральную нагнетательную скважину радиусом гс начала осуществляться непрерывная закачка водного раствора полиакриламида с концентрацией с°, т. е.

s (r c, /) = 1 — s„ осТ( с°) = S°; c (rCt t) = c

(6 .5 9 )

240

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