close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
ПРОЦЕССЫ И АППАРАТЫ ХИМИЧЕСКОЙ ТЕХНОЛОГИИ
УДК 543.4:544.2
С. П. Плохотников, В. А. Богомолов, Е. Н. Белова,
О. И. Богомолова, Д. С. Плохотников, О. Р. Булгакова
ОСРЕДНЕНИЕ МОДЕЛЕЙ ТРЕХФАЗНОЙ ФИЛЬТРАЦИИ В НЕОДНОРОДНЫХ СЛОЯХ,
ПОДЧИНЯЮЩИХСЯ РАВНОМЕРНОМУ РАСПРЕДЕЛЕНИЮ
Ключевые слова: фильтрация, фазовые проницаемости.
В работе построены осредненные модели трехфазной фильтрации для неоднородных пластов. Неоднородность задавалась с помощью равномерного распределения, частного случая β–распределения при коэффициентах γ=1, η=1. Проведен сравнительный анализ построенных моделей с эталонными трехмерными моделями.
Keywords: filtration, phase permeability.
The averaged three-phase filtration models of inhomogeneous layers have been built in this work. The heterogeneousness was specified by beta– distribution law with γ=2, η=1. The compared analysis of filtration models with etalon 3D
models has been done.
Введение
ность, подчиняющуюся равномерному распределению.
При многовариантных гидродинамических
расчетах можно использовать трехмерные численные гидродинамические модели многофазной
фильтрации. Гидродинамические модели создаются на основе геологических моделей, которые
содержат до 100 млн. ячеек разностной сетки, время расчетов которых неприемлемо для практического использования. Поэтому делают ремасштабирование (up’scaling), объединяя мелкие ячейки
геологической сетки в крупные ячейки гидродинамической сетки. В работе с помощью вычислительного эксперимента исследованы возможности
применимости известных [1,2] осредненных моделей, в том числе с модифицированными относительными фазовыми проницаемостями (ОФП), исследованных ранее при двухфазной фильтрации [36]. Также создана методика осреднения для трехфазной фильтрации [7]. В данной работе рассматривалось трехфазное изотермическое вытеснение нефти и газа водой при площадном заводнении в слоистом пласте – пятиточечной и девятиточечной системах заводнения. Пятиточечная система заводнения – одна нагнетательная скважина в центре квадрата, а вокруг – 4 добывающих скважины. Девятиточечная система – одна нагнетающая скважина в
центре, а вокруг – 8 добывающих скважин.
Построение математической модели
Трехфазная фильтрация в слоистом пласте с
проницаемостью k(z) и исходными ОФП, обозначим
их k rw (S g , S w ), k ro (S g , S w ), k rg (S g , S w ) , была сведена к фильтрации в однородном пласте со средней
проницаемостью k* и модифицированными ОФП
вида [7] (Модель B1):
k (S )  k (S )
k mrw (S w )  rw w * w w ,
k
k mro (S g , S w ) 
k rw (S gfix , S w )  k o (S w )
k mrg (S gfix , S w ) 
k*
k rg (S gfix )  k o (S w )
k*
,
.
где k w (Sw ) , k o (S w ) – средние значения абсолютной проницаемости в зонах воды и нефти в каждом
вертикальном сечении пласта.
Модель B1 – осредненная модель, которая
получается путем осреднения по толщине пласта
проводимости для каждой из фаз. В ней учитывается
наличие газа только в зоне нефти. При осреднении
вводим следующие физические допущения: независимость давления Р от z , в каждом вертикальном
разрезе пласта существуют две зоны – зона воды
мощностью H w и в ней S w = S *w – максимальное
Цель работы
Для моделей трехфазной фильтрации построить двумерные модели со средними ОФП и с
модифицированными ОФП, полученными с учетом
струйности течения по пропласткам для равномерного
распределения,
частный
случай
βраспределения при γ=1, η=1. Провести сравнительный анализ с помощью ВЭ трехмерных и двумерных моделей трехфазной фильтрации. При этом для
трехмерных моделей задать слоистую неоднород-
значение S w , S g  S g * – минимальное значение
Sg , и зона нефти мощностью Hoil и в ней S w  S w * ,
S g  S gfix , исходные параболы относительных проницаемостей заменяем на секущие прямые. Но целесообразнее подправить их с помощью поправочных
99
Уравнение Ошибка! Источник ссылки не
коэффициентов по аналогии с работами [3-6], так же
как со случаем двухфазной фильтрации.
По аналогии с двухфазной фильтрацией [36], получим модифицированные ОФП для трехфазной модели.
k S   k *  J k
k mrw (S w )  rw *w
S w (S w )  k *
,
fix
(1)
k
(
S
,
S
)
J k
ro
w
g
k mro (S g , S w ) 
*
*
1 S w (S w ) k

найден. будет иметь вид:
1
dk  0
a ba
k
So  
 
ab
.
2
Решая уравнение относительно k , получаем
корень уравнения – k  b  (b  a)  S *w
.
Из (1) получим модифицированные относительные проницаемости для линейного случая имеют вид:
 S *  (b  a) 
,
k rom  k ro  1  w
b  a 

Средняя проницаемость будет k * 

k mrg (S w , S gfix ) 

k rg (S gfix )  J k
1  S
*
w
(S w ) k *
k
   k  f k dk ,
где J k 
а
0 *
krw(Sw )  krw
Sw , krg(Sgfix )  krg0 Sgfix ,



(2)
m
k rw
 k rw  1 
kro(Sw,Sgfix)  kro0 1 S*w  Sgfix ;

0
ro
    1,2,3 , k – максимальная ОФП нефти;

k rgm  k rg  1 

0
k rw
– максимальная ОФП воды; k rg – максималь-
k
 1  Sw*  Sgfix , k *   k  f (k )  dk – средняя
0
проницаемость, S *w 
(S w  S wc )
–
(1  S wc  S or  S gc )
вижная
S gfix (S g , S w ) 
под-
водонасыщенность,
(S g  S gc S )
*
w
– фиксированное значе(1  S *w )
ние газонасыщенности .
При этом значение k (Sw ) находим, решая
уравнение
k S 
(3)
S wc  S w
  f (k )dk .
1  S or  S wc
0
Получим модифицированные ОФП для для
трехфазной модели в случае равномерного распределения.
Плотность вероятности обобщенного β – распределения имеет вид:
Г(   )
1

f (k ) 

b  a Г(  )Г()
(4)
 k  a   1  k  a   1

 1 


ba
 ba
α≤k≤b, γ>0, η>0, где Г(η) – гамма функция.
Равномерное распределение - частный случай
β-распределения при коэффициентах γ=1, η=1.
Подставляя их в Ошибка! Источник ссылки
не найден. получим
(5)
1
k mrw  k rw  [1    3  (1  S *w )] ,
k mro  k ro  [1    3  S *w ] ,
(6)
k mrg  k rg  [1    3  S *w ]
w
f (k ) 
(1  S *w )  (b  a) 

ba

Можно получить коэффициент вариации слоистой неоднородности для равномерного распределения
b  a
(k ) 
a  b  3
и записать соотношение:
b  a  .
(k )  3 
a  b 
И произвести замену в формулах. Тогда для
трехфазной модели модифицированные ОФП примут вид:
ная ОФП газа; Sor – остаточная нефтенасыщенность;
Swc – насыщенность связанной воды; Sw, So – водои
нефтенасыщенности,
Swc≤Sw≤1Sor, So
(1  S *w )  (b  a) 
 ,
ba

На рис.1 показаны графики зависимости моди-
Рис. 1 - График зависимости модифицированных ОФП
(нефть-вода) от водонасыщенности Sw для равномерного распределения при линейных исходных ОФП
ba
100
фицированных ОФП (нефть-вода) от водонасыщенности Sw для равномерного распределения при линейных исходных ОФП. Графики модифицированных ОФП для системы нефть-газ выглядят аналогично.
другом. Пропластки располагались относительно
друг друга следующим образом: «лучший» (максимальное значение абсолютной проницаемости) рядом с «худшим» (минимальное значение абсолютной проницаемости), «лучший» из оставшихся рядом с «худшим» из оставшихся и т.д. снизу-вверх.
Вычислительный эксперимент
Задача была решена при заданном перепаде
давлений между нагнетательной и добывающими
скважинами,
внешняя
граница
задавалась
непроницаемой. При расчетах использовались сетки
из блоков: 11x11x10 (x,y,z) для эталонного
трехмерного случая; и 11x11x1 для двумерного
случая (ремаштабированная). В
расчетах
использовалась полностью не явная схема.
В модели были заданы следующие физические
параметры: 128 – начальное пластовое давление,
атм; 22 – температура пласта, С; 100 – температура
закачиваемой воды, С; 55 – забойное давление на
добывающей скважине, атм.; 170 - забойное давление на нагнетательной скважине, атм.; k ro0 =0.5 –
0
=0.7 – максимальная
максимальная ОФП воды; k rw
ОФП нефти; k rg0 =0.8 – максимальная ОФП газа;
Sor=0.2 – остаточная нефтенасыщенность; Swc=0.2 –
насыщенность связанной воды; Sgc=0.1 – насыщенность защемленного газа; Sw,So,Sg – водо-, нефте- и
газонасыщенности, Swc ≤Sw≤ 1- Sor.- Sgc.
Расчеты для пятиточечной и девятиточечной систем проводились для четырех вариантов:
Двумерные решения (осредненные):
1. C – решение – осредненное решение с
линейными исходными ОФП krw(Sw), krg(Sg,Sw),
kro(Sg,Sw), вида (2) и средней абсолютной проницаемости k  0.5 дарси . Задавался 1 пропласток,
высотой H = 10м
*
Рис. 2 - Зависимость суммарного объема добытых нефти и газа Vp от времени Date для пятиточечной системы заводнения, и соответственно для решений –
A7,A8,C,B1
Трехмерные решения (эталонные):
2. В1 – решение – осредненное решение, газ
находился только в связанном состоянии в зоне
нефти, задача решалась с модифицированными
ОФП k mrw (S w ) , k mro (S w , S gfix ) , k mrg (S w , S gfix ) вида (6) и
На рис.2 хорошо видно, что осредненное решение B1 с модифицированными ОФП имеет меньшую
погрешность относительно эталонного решения A8,
а решение C с исходными ОФП меньшую погрешность относительно эталонного решения A7. Такие
же результаты получены для девятиточечной системы заводнения, и для других показателей разработки, как для газа, так и нефти. Решения B1 и C ограничивают все трехмерные эталонные решения Ai для
которых проводились вычисления. Это говорит об
обоснованности применения этих двух осредненных
моделей в совокупности при площадном заводнении
в слоистых пластах при трехфазной фильтрации.
Выводы
Сравнительный анализ зависимости показателей разработки для различных трехфазных моделей,
показал правильность построения осредненных
трехфазных моделей. Осредненные трехфазных мо-
средней k  0.5 дарси для равномерного закона
распределения задания абсолютной проницаемости
k(z) исходного слоистого пласта. Задавался 1 пропласток, высотой H = 10м.
3. A8 –решение (с изолированными пропластками) – эталонное численное трёхмерное решение
задачи для десятислойного пласта с изолированными пропластками (отсутствуют перетоки), распределение абсолютной проницаемости которых подчинялось равномерному закону. Задавались 10 пропластков каждый высотой H1 = H2 = H3 = H4 … H10 = 1м,
изолированных друг от друга непроницаемыми перемычками;
4. A7 – решение (с неизолированными пропластками) – тоже, что и предыдущая модель, но с
неизолированными пропластками. Задавались 10
пропластков гидродинамически связанные друг с
*
101
Марвин, Р.Х. Фатыхов // Вестник Казан. технол. ун-та. 2005. № 1. – С.121–124.
5. Плохотников С.П. Математическое моделирование
трёхфазной фильтрации в слоистых пластах с учётом
схемы струй / С.П. Плохотников, Д.С. Плохотников,
В.В. Елисеенков, А.С. Климова // Вестник Казан. технол. ун-та. - 2005. - № 2. – С.173–178.
6. Плохотников С.П. Осредненные модели двухфазной
трехкомпонентной фильтрации при закачке в нефтяной
пласт химических реагентов – полимеров, водных растворов ПАВ. / С.П. Плохотников, В. А. Богомолов, О.Р.
Булгакова// Вестник Казан. технол. ун-та. – 2010. - № 10.
- С.350–356.
7. Bogomolov V.A. "Mathematical simulation of three-phase
filtration in stratified beds with account for the scheme of
jets" V.A. Bogomolov, S.P. Plokhotnikov, O.R. Bulgakova,
D.S. Plokhotnikov // Journal of Engineering Physics and
Thermophysics, - Springer, 2011, Vol. 84, No. 5, pp. 975–
979.
дели показали хорошие результаты, что говорит об
обоснованности применения этой методики осреднения и для трехфазных моделей фильтрации.
Литература
1. Булыгин Д.В. Геология и имитация разработки залежей
нефти / Д.В. Булыгин, В.Я. Булыгин. – М.: Недра, 1996.
– 382 с.
2. Методические указания по созданию постоянно действующих геолого-технологических моделей нефтяных и
газонефтяных месторождений (Часть 2. Фильтрационные модели). – М.: ВНИИОЭНГ, 2003. – 228с.
3. Плохотников С.П., Елисеенков В.В. Гидродинамические
расчеты в слоистых пластах на основе модифицированных относительных проницаемостей// ж. «Прикладная
механика и техническая физика», (ПМТФ), Новосибирск, РАН СО, т.42, №5, 2001, с. 115–121.
4. Плохотников С.П. Модифицированные фазовые проницаемости в задачах площадного заводнения слоистых
пластов / С.П. Плохотников, Д.С. Плохотников, О.Б.
___________________________________________________
© С. П. Плохотников – д-р техн. наук, проф. каф. ИПМ КНИТУ, plokhotnikov@kstu.ru; В. А. Богомолов – ст. препод. той
же кафедры, bogomolov@kfti.knc.ru; Е. Н. Белова –асп. той же кафедры; О. И. Богомолова – асс. той же кафедры;
Д. С. Плохотников – асп. той же кафедры; О. Р. Булгакова – асп. той же кафедры.
102
1/--страниц
Пожаловаться на содержимое документа