Файл: Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Проектирование разработки.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.02.2024

Просмотров: 299

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

а

в

 

 

U V s

V

1

Рис. XVI.9. Схема выделения зон.

а —вокруг нагнетательных скважин; б —вокруг добывающих скважин; в —зона сек­ тора и прилегающей ячейки

При равномерном разбиении на кольца гi = 0,25/Дх.

I = 0,196s! +

0,589s2

0,215s3;

$i+it j ~ ^i-i}j ~

$it y+i ^

S/1/-1 = 0,19s3 4~ 0,2665^ -j- 0,266s5 -f- 0,2786sg.

 

 

(XVI.37)

Числовые коэффициенты сохраняются и при определении средних тем­ ператур по ячейкам. Здесь множителем является удельная (объемная) тепло­ емкость каждой из кольцевых зон. Например, для центральной ячейки

0,196^7! + 0,589С272 + 0,215С373 0,196С! + 0,589С2 + 0,215С3

В ы ч и с л е н и е д е б и т о в н а г н е т а т е л ь н о й с к в а ж и н ы

по

п р о п л а с т к а м

 

 

 

 

 

 

 

Расход

 

 

 

 

 

 

 

 

 

 

 

V2 Дх

 

 

 

 

 

 

 

 

Я=

РСК-

J

(1 -П С „ (А .Г )А --

 

 

 

 

 

Pi, 1 + Р_1,1 -4- р1, _1 -Ь Р-1,-1

«в.

 

 

 

(XV1.38)

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

1 ,

0,25 Дх

 

1

i

1

In

4 V~2

 

Лв =

2nk

— In---------

La

Xi

1- 1

Xe

 

 

Xi

rc

 

 

 

у =

---

]___ ------

 

 

 

 

 

 

 

x

Hb(T)

+ Рн(Г)

 

 

 

 

 

 

 

При определении дебитов следят за положением скачка насыщенности и зна­

чением температуры. Если до прохождения скачка

насыщенности температура

в слое достигла

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

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

чается от работы.

343


Р а з б и е н и е п р и з а б о й н ы х зон

д о б ы в а ю щ и х

с к в а

ж и н.

Призабойная зона добывающих скважин в пределах одного шага сетки

(гс <

г < А*) разбивается на четыре кольцевые

зоны (г* =

/Дх/4) и

четыре

сектора, соответствующие четырем соседним узлам сетки. В

каждом

секторе

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

Средняя насыщенность в узле, где расположена добывающая скважина,

соответствующая средним значениям по секторам и

по радиальным кольцам,

4

(XVI.39)

s, 7 = 0,25 2 (0,196sy + 0,589s^+0,215s^).

i= 1

 

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

„ _

S4 +

ss

*

_

^4 + *3

 

 

bi , j —

2

 

’ */»/—

2

9

 

 

s4 ==

j

j

S3 =

2S/, j

s3 -}- ds[t j == s4 4~ dsi,j,

 

dslf j = 0,5 (ds3 4- ds4),

 

$4 =

2sif / — s3.

 

(XV1-40)

Знаком

Д

обозначены

новые

значения

насыщенности.

Р а с х о д

по с е к т о р а м

д о б ы в а ю щ и х

с к в а ж и н

— Реек == (Pit / — Рек — Рпот)/#д»

 

 

(XVI-41)

где

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

Рпот -

0 ~

/сек) Gn (k' Т)

^д = ~nk 2

"хГ ^ ri-l

г/= 0,25/Ах;

г0 = гс.

 

 

i=1

 

 

 

 

 

 

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

4 м

* =

S

C i! Q =

Ё

qZ-

 

 

 

 

 

V=l

 

 

2=1

 

 

 

 

 

О б в о д н е н н о с т ь

 

 

 

 

 

h r

( s i

Т]) q lJ q .

 

 

 

 

 

V—1

 

 

 

 

 

 

 

 

 

В ы ч и с л е н и е д а в л е н и й в у з л а х с е т к и

 

Уравнение

(XVI.32)

заменяется

разностной аппроксимацией

 

J _ /

 

Pi+V j

Pi, j

 

Pit j~~ Pi-l,/\

 

 

Ax\ V

4 - . /

A*

- f /

A*

j

 

4-

 

X.

Pi, /+1 — P/»/

 

Pi* i ~ Pit i-1

= 0.

(XVI.42,)

 

1 '

Ap

— X

Ap

*■'4- —

 

 

AP

 

 

344


При

Да: =

А// для

центральной

точки

 

 

 

X. .

1

Р .,,

. +

Х.

1

Р.

.

. + Х.

1 Р.

 

 

 

 

 

 

i .

р

 

 

 

 

 

 

 

 

+

X.

2

 

 

 

 

 

 

 

 

 

.

1 р.

.

,

 

(XVI.43)

Л/

1C. ,

1 ,+ЗС.

1

. + Х.

 

1 +Х. . 1

 

 

 

 

1+ Т ’>

‘" Т

1'

‘• '+ т

• . / - J -

 

где х — гидропроводность с учетом структурных свойств нефти и воды,

■ F n (s,T )f ,

(XVI.44)

X соответствует данным между узлами сетки, средние значения которых опреде­ ляют по закону среднего гармонического:

 

4 b /Xi+i, / .

 

 

j%j+1> i

i'-h*2”. /

X/» / “b Xi+1’ /

i, /4—2”

Xi* 7

Xl>j+1

Для узлов, соседних с добывающими скважинами, разностная аппроксима­ ция уравнения для давления записывается так:

{

i

I х - Т* 1

/+т

Pi+1» j — Pi* j

1

<

Pi’ J+l Pi’ j

,

Ду

+

h - 1

= °-

Pi» j

Pi-1» j У

Др +

,

Да:

 

(XVI.45)

Здесь <7сек — расход по третьему (верхнему)

сектору добывающей скважины.

В данном случае

(XVI.46)

Ясен — (Pi* f-i — Pi. ,)//?д,

г. .

,

= * д.

 

‘• '

- т

д

 

Для узлов, соседних с нагнетательными скважинами, пользуются этой же формулой, но вместо ^Сек подставляют 0,25^Скв-

Расчет насыщенности в призабойной зоне нагнетательных скважин вы­ полняется по явной схеме с выделением фронта насыщенности до момента его выхода за пределы радиальной зоны. Насыщенность при этом определяют из условия, что дробь /* (s) — /* (s0)/(s — s„) принимает наибольшее значение при фиксированных температуре, скорости фильтрации и проницаемости.

Приращение насыщенности в кольцевой зоне

А

^скв {fi—1

fi)

А.

(XVI.47)

 

= -----71-----2“^Г ДЛ

 

 

ЯШ

j)

 

 

Здесь

/* — доля воды

в

потоке с учетом ее структурных

свойств.

Временной шаг выбирают по наиболее проницаемому слою из условия ус­ тойчивости разностной схемы с учетом того, чтобы за этот период фронт вытес­

нения не проходил более

одного

координатного шага,

.

пт (f*+, -

г»)

(XVI.48)

at = min

------—£ --------.

 

^скв^в (sck)

Здесь v — номер временного слоя; sCk — скачок насыщенности.

345


 

Приращение насыщенности в ячейках сетки

 

 

 

 

 

As = г

о

 

 

; -

NZ , / -

NZ , i + NZ ,

/->),

 

(XVI.49)

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wi =

+

U[_i, j Ax;

N2 = -f ui, j Ax-,

 

 

 

 

 

 

N3 = + v i,j Ay;

Ni = + Vitj_1Ay;

 

 

 

 

 

 

 

 

 

, (

PQ

A

 

 

t' \

 

 

Pi+1, j — PL j

 

 

uH = ~ k Vhb(T)

^ „ ( Л

6 )

 

 

 

 

**

 

 

 

 

 

 

 

 

 

 

 

 

 

'+ T - '

 

 

 

 

 

Насыщенности в призабойной зоне добывающих скважин определяют по

методу

разностного аналога типа уголок

 

 

 

 

 

 

m

Si — St

,

4б/ск

/ p* '

—si+i

 

к/

T’i — T’i+l

 

(XVI.50)

 

 

M

 

1

я

 

S r2 Y-r9

■+ /r

r2 — r2

= 0.

 

 

 

 

 

 

 

 

 

ri

r rt+l

 

 

ri

ri+ 1

 

 

 

В

качестве производных

принимают

 

 

 

 

 

 

fs

=fs

(si+ 1»

 

 

fa

= fr

(si+ ь

 

 

 

 

 

Sj означает

значения

на

предыдущем

временном

слое.

имеет вид

 

 

 

Разрешенное

относительно s*

уравнение

(XVI.50)

 

 

..

 

5г + [ / ; ч + . + / г '( 7 ,,+1-

 

г .) ] ф

 

 

 

(XVI.51)

 

 

 

 

 

 

 

1+ /,ф

 

 

 

 

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4(/ск А/

 

 

 

 

 

 

 

 

 

 

 

 

Ф = —

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Производные в каждом из элементов призабойной зоны рассчитывают по

формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/5*' =

50 [/*(s 4- 0,01; Т) -

/* (s -

0,01; Г)],

 

 

 

 

 

/*' =/*(s, T + 0 ,5 ))-/* (s;T — 0,5).

 

 

 

 

(XVI.52)

 

Если температурный фронт не достиг забоя добывающей скважины, то /f' =

=

0 и формулы упрощаются.

 

 

 

 

 

 

 

 

скважин

по

При расчете температуры в призабойных зонах нагнетательных

формуле

Ловерье

теплоемкость всего

продуктивного

разреза

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(rv ‘

 

 

НС. =

С (Н - Н Эф) +

2

hv

 

 

 

 

>V'ck;

(XVI.53)

 

 

~Ь (С11В— Сп,|) Sj--

 

где /=

1,5Ддг;

— скачок

насыщенности на

фронте.

 

 

по

Т е м п е р а т у р у

в

у з л а х

 

к в а д р а т н о й

с е т к и вычисляют

разностному

уравнению

 

 

 

 

 

 

 

 

 

 

. .

 

 

 

 

|

Ят

I

, 1

~

ЯТ

|

 

 

 

 

 

CzATi,i — Cz^Tij

‘+т

 

‘~ т - ‘

 

+

 

 

 

 

 

 

 

At

 

 

 

 

 

 

Ах

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЯТ

 

, - q T

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

•'+ т

ду

~

 

+

V

w AT>'i =

 

 

 

(XVI.54)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

340


Здесь

А Г = (Г — Г0); q = J]/tv (^в^в г СнТ’нУ^

^о) — конвективный

 

v

 

поток тепла.

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

Новое значение теплоемкости пласта, отнесенной -к единице площади,

Сг = С (Н — Ндф) +

+ 2 Y ^ +

2

 

ICb*V+

С“ (1 -

sV)l-

(XVI.55)

Новое значение температуры

 

 

 

< 7 -T ,)w -

 

 

 

 

 

 

С дг

I

/

т

-<77

Я?

- q T

+ ^

 

г

 

Г

 

«'+ -j-. /

./+ т )

 

 

ч . <

 

 

 

 

С2 4- А/ г

TtC

 

 

 

 

 

д/

 

(XVI.56)

Если узловые точки — скважины, необходимо учитывать потоки тепла, вносимые или уносимые жидкостью. Так, например, характер потоков вдоль оси х определяют по четырем вариантам, соответствующих случаям: нагнетатель­ ная скважина, узловая точка; добывающая скважина, узловая точка; узловая точка, нагнетательная скважина; узловая точка, добывающая скважина.

ят

1 ‘

L

. =

2

AV<г ‘ • / -

т °> W

(Ѱð +

• /»

 

 

2

' 1

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qV

 

 

 

 

ят .. J_

. =

2 * v (т '+>. / -

го ) ^ Г

(с -£ + с вА*).+1. /.

 

^

 

,

. =

-

2

ftV (г ж

. / - го) T Z j

(CJI +

c / r) t+I.

 

 

i H——, /

 

V

 

 

 

 

^

 

 

 

 

~ 2

 

 

 

 

 

 

 

 

 

 

 

 

^

. =

-

2

ftV (r <. / -

r o)

 

 

(Cȣ +

,.

(XVI.57)

Для

тепловых

потоков вдоль

оси

у /qT

, \

формулы

аналогичны.

При

в ы ч и с л е н и и

 

 

V

 

/+ - )

в п р и з а б о й н ы х

т е м п е р а т у р ы

зо н а х

д о б ы в а ю щ и х

с к в а ж и н

уравнение записывается для сек­

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

347