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

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

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

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

Добавлен: 29.02.2024

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

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

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

Здесь

/ — функция Баклея—Леверетта

(функция

распределения фаз);

п — функция, определяющая долю

воды в

потоке.

 

 

Скорость движения фронта вытеснения (скачка насыщенности) по слою /

drji

w

г t* / \

-* / vi

Qj fUs) — /b(so)

(XVI. 17)

dtГ =

®ск = m (s - s 0)

(s) -

(so)]

2nrm

s — So

где s — насыщенность непосредственно за скачком (фронтом вытеснения); s0 — насыщенность перед ним (связанная вода).

Когда фронт застывания (в связи со снижением температуры до критического значения) нагоняет фронт скачка насыщенности, то насыщенность на фронте мо­ жно определить из соотношения

Fн (s) Оц (k, 7\ s) = p,D(Т) w.

(XVI.18)

В этом случае режим будет поршневым.

Когда переход на поршневой

режим отсутствует, то скачок насыщенности

оценивается как точка, где достигается максимум функции

(So)

s — So

Здесь sQ— связанная водонасыщенность. При этом доля воды /£ определяется

из (XVI. 16) и зависит она как от структурных свойств, так и от суммарной ско­ рости фильтрации.

Уравнение для насыщенности в слое / (при m = const)

dsi ,

Q/ dfh

(XVI.19)

dt 2Jimr dr

Расчетный конечно-разностный аналог этого уравнения для ячейки ij, через которую прошел фронт вытеснения к моменту А/,

Qj&t (/в i_i, j

!в ij)

(XVI.20)

Asij =

 

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

Если фронт вытеснения расположен в t-й ячейке и не покинул ее за время А/, то приращение насыщенности в данной ячейке

Qj

,_i, /

[■(XVI.2I)

As*/ =

 

n m { r] - rU ) ’

 

гДе f! /_i берется

в промежуточный момент времени.

 

Если фронт вытеснения за взятый отрезок времени пересекает границу ячейки, то определяется время Atp и Atq, когда фронт вытеснения находится в t-й

или (i +

1)-й ячейках,

тогда

соответствующие приращения насыщенности

Мр =

Jim (г? — т\

)

(XVI.22)

' Л

? , -Т - ’

 

Q/f, (sck)

 

 

tatq = М 4- Atp.

 

 

Приращение

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

в (i + 1)-й ячейке

 

Qj

(sck)

(XVI.23)

As,-+iлт (/■?_, -г* ) ‘

 

Для предыдущей ячейки изменение насыщенности определяется следующим образом. На первом этапе приращение определяется по формуле (XVI.21), в ко­ торой At заменяется на Atp. На втором этапе используется формула (XVI. 19) с уточненными на момент t + Atp значениями доли воды в потоке.

333


Рис. XVJ.1. Температурное поле в четырех­

Рис. XVI.2. Значение насыщенности вы­

слойном пласте при нагнетании горячей (/)

тесняющей фазой (водой) на различных

и холодной (2) воды на момент обводнен­

расстояниях г в пластах с проницаемостью

ности продукции на 93 % (начальная пла­

0,6 мкм2(1, 2) и 0,06 мкм2(3, 4) при на­

стовая температура Т0 = 65 °С)

гнетании горячей (Г = 90 °С) и холодной

 

 

(7=12 °С) воды на момент обводненности

 

 

продукции на 93

%.

 

 

1,3 —при нагнетании горячей воды; 2,

 

 

4 —при нагнетании холодной воды

Расход воды по пропластку

 

 

Q/

dr

1_ Цн (Гр) .

Rk

I 2nk (Гп/Цп (Т) + Гц/ц,, (Т)) г

2nkF„ (s„)

/-С1(

 

И

1

 

=

Др— J GM(k, Т) dr — j [1

 

1 ~Ь ^вМ^н (T)/Fнр,в (Т)]о„ (k,T)dr.

(XVI.24)

До прорыва воды движение по слою возможно лишь при условии Ар >

> J (7Н(/?, Т) dr, в противном случае пропласток считается застывшим.

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

Вычисления показали, что при таком выборе шага расчетная схема устой­ чива.

В качестве примера приведем результаты расчета процесса вытеснения вязко-пластичной нефти в неоднородном слоистом пласте при нагнетании холод­

ной и горячей вод.

д ан н ы е .

Четырехслойный пласт со значениями прони­

И с х о д н ы е

цаемостей:

а) 0,6,

0,040,

0,06, 0,6 мкм2 и б) 0,4, 0,1, 0,05, 0,2 мкм2, толщина

каждого слоя И= 4 м;

RK= 200

м, m = 0,23.

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

7\ = 12 °С,

при нагнетании

горячей воды Г1=90°С. Начальная пластовая

температура

Г„ =65°С.

Перепад

давления Ар = 4 МПа.

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

Для

вязкости

воды

и

нефти

при 7’>20°С

Кв (Т)

 

46J

 

 

m

_

151

 

Т +15,7°’

^

'

'

Т — 20°'

 

 

334


Рис. XVI.3. Значение общей нефтеотдачиц

Рис. XV1.4. Показатели текущей нефтеот­

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

дачи л в четырехслойном пласте (с прони­

горячей (/) и холодной (2) воды в различ­

цаемостью слоев 0,4; 0,1;

0,03;

0,2 мкм2)

ные периоды времени

 

 

 

 

при закачке горячей (Г = 90 °С)

и холод­

 

 

 

 

 

 

 

ной (Т = 12 °С) воды.

 

 

 

 

 

 

 

 

 

Вода:

 

 

 

 

 

 

 

 

 

1

—горячая. 2 —холодная

 

При Т <20°С

нефть считалась малоподвижной, что обеспечивалось принятием

рн (Т <,20°)=

1500

мПа-с.

Начальный

градиент давления

сдвига

^

^

0,00158

/ с

Г

\

f мПа 1

 

 

Gu iT,k) = - v

r ^

( 5 - w

)

[ П Т

J-

 

 

Фазовые проницаемости

 

 

 

 

 

 

 

 

FB=

( f E T i) 3 01 ’

f " =

 

) 2 2 nP"

si = °-20 + °-001

~ 65°) 11s? =

= 0,8+ 0,002 (T — 65°). Шаг по координате г был равен 5 м. Результаты расчетов представлены на рис. XVI. 1—XVI.4 и в табл. XVI. 1.

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

Так, в случае нагнетания горячей воды обводненность продукции до 93 % на расстоянии 200 м была достигнута через 22,5 месяца разработки, при нагнета­ нии холодной воды — только через 42 месяца. Причем общая нефтеотдача в первом случае была значительно выше, чем во втором, так как при нагнетании

ТАБЛИЦА XVI.1

ЗНАЧЕНИЯ НАСЫЩЕННОСТИ su И ТЕМПЕРАТУРЫ НА РАЗЛИЧНЫХ

РАССТОЯНИЯХ В

ПЕРВОМ

СЛОЕ ПРИ

ОБВОДНЕННОСТИ

ПРОДУКЦИИ

НА 95 %

 

 

 

 

 

 

 

 

30

40

50

60

70

80

90

Г, Сс

15,6

18,9

23,5

29,9

37,8

46,3

53,83

SB

0,675

0,718

0,729

0,735

0,739

0,741

0,739

м

100

 

120

130

140

160

200

т, °с

59,3

64,42

64,9

64,97

65,0

65,0

sB

0,732

0,712

0,701

0,691

0,673

0,644

335


холодной воды вследствие массового выпадения парафина «отключались» низко­ проницаемые слои.

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

Расчеты по данному методу могут выполняться на ЭВМ средней мощности.

§ 2. РАСЧЕТ ПРОЦЕССА

ВЫТЕСНЕНИЯ НЕФТИ ВОДОЙ

В НЕИЗОТЕРМИЧЕСКИХ

УСЛОВИЯХ

ПРИ ПЛОЩАДНЫХ СИСТЕМАХ ЗАВОДНЕНИЯ

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

Расчет показателей неизотермического вытеснения нефти водой в пласте при пятиточечной схеме расположения скважин

Данная схема размещения скважин широко применяется в промысловой практике [29].

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

Рис. XVI.5. Схемы расположения нагнетательных и добывающих скважин при различных системах заводнения.

а —радиальная фильтрации; б —нагнетание и линейные ряды (добывающие ряды ап­ проксимированы добывающей галереей); « питнточечная система заводнения; г —семн- точечнаи система; д —семиточечная обращенная система

336

Учитывая

симметричность

расположения

 

скважин, показатели вытеснения

определяются

 

для одного из элементов схемы

(рис.

XVI.5,

 

XVI.6).

 

ниже

метод показал свою эф­

 

Излагаемый

 

фективность в расчетной практике.

 

по­

 

Рассматриваемая

область фильтрации

 

крывается

прямоугольной

сеткой. Нагнетатель­

 

ная скв. 1

расположена

в

точке

(0, 0),

добы­

 

вающая скв. 2

в точке с координатами а, Ь. В

 

скважину нагнетается вода

с температурой, от­

 

личной от начальной

пластовой.

 

 

 

 

Пусть шаг сетки по координате х равен Д*,

 

поу — Д*/;

 

0

у

Ь.

 

 

Рис. XVI.6. Схема расположе­

Обозначив число шагов по х через пх, а по

у — через пи, имеем а =

t±xnx\ b —купу*

i

ния эксплуатационных и нагне­

тательной скважины при пяти­

= 0, 1, 2,

пх\ j =

0,

1,

2,

 

 

 

точечной схеме заводнения.

Влияние логарифмической особенности поля

1 —нагнетательная скважина;

давления

вблизи скважины

можно учитывать

2 — добывающая скважина

двояким способом. Г. Г. Вахитовым было показа­ но, что решение разностной задачи для поля давления совпадает с аналитиче­

ским, если скважину поместить в узле и шаг сетки взять в 5 раз большим. Можно ввести дополнительное сопротивление, тогда дебит скважины опре

деляется уравнением

2nHkFH

pc — рц d= Он (k, T) Дх

Q = И*н (^)

In (Дх/5rc)

где pc — давление на забое скважины; Pij — давление в узле сетки, в которой

расположена

скважина.

Было предложено считать в призабойной зоне проницаемость условной, т. е.

kx = Kk\

К = In (Дх/гс)

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

профиля давления.

Процесс неизотермической фильтрации в этом случае описывается системой дифференциальных уравнений

(XVI.25)

(XVI.26)

{c . + i V - т г ) - ж + 1с -,: + с " ( | - , *)1

(XVI.27)

337


Здесь

с* wti = fuw;

GK(T,k)

I Vp | '

u — проекция скорости w на ось Ox, v — проекция скорости w на ось 0у.

В (XVI.27) предусмотрен теплообмен пласта с кровлей и подошвой по закону (асимптотическое приближение уравнения Ловерье [29])

(XVI.28)

Согласно данному методу параметры давления, температуры и насыщенности определяют в узловых точках сетки, а компоненты скорости — в промежуточных, средних между соседними узлами точках: компоненты и (х, у) — между сосед­ ними по горизонтали точками, v (х, у) — между соседними по вертикали точ­

ками.

Т — от 0 до пх и от 0 до пу,

Граничные пары индексов: для массивов р, s,

для массива и — от 0 до пх — 1 и от 0 до пу\

для массива v — от 0 до п[

и от 0 до пу — 1.

 

Уравнения для давления во всех узлах сетки (кроме узлов скважин) аппрок­ симируются разностной схемой

 

 

(XVI.29)

Для узлов, расположенных на границе

с координатами х =

0, х = а\

у = 0, у = Ь, применяется метод зеркального

отражения значений

давления

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

С учетом структурных свойств нефти процесс итерации проводится дважды. Первый раз фиксируются подвижности системы. Когда давления станут близкими по значению, они снова пересчитываются в связи с влиянием градиента давления. При пересчете требуется всего 2—5 итераций (при точности 10-2).

Уравнение для насыщенности аппроксимируется явной разностной схемой типа уголок

(XVI.30)

338