1275
.pdfпарциальное давление |
кислорода, Па; |
Е — энергия активации, |
||
кДж/( моль-К); |
Rr— |
универсальная |
газовая |
постоянная, |
кДж/(моль-К); Т — абсолютная температура, К; К — коэффициент, учитывающий влияние степени окисленности нефти на скорость
реакции; / — степень окисленности нефти, доли единицы; п — по казатель степени при парциальном давлении кислорода, доли еди ницы.
Энергию активации необходимо определять с использованием метода наименьших квадратов. Для этого зависимость (5.11) не обходимо привести к линейному виду. Прологарифмировав правую и левую части (5.11), получим
1пау=1пЛ — В (1/Г), |
(5.12) |
где
Л = Л0р£2ехр ( — Kf)\
B = EIRr.
Так как в опытах f = 0 и все величины скорости реакции при ведены к одному значению парциального давления кислорода, то
Л = Л0ро2= const.
Тогда значение коэффициента В вычисляют по формуле
где кЭг/— |
число замеров скорости |
реакции; |
winp — скорость ре |
||
акции |
нефти с кислородом воздуха |
в i-u |
опыте, кг/(кг-с); 7\- — |
||
температура кернодержателя в i-м опыте, |
К; |
||||
£ _ |
6-0,351678— 113,74071-0,0184258 |
_ |
_ |
8133 |
|
_ |
6-56,87815-10-°— (18,4258)М0-С |
~ |
|
|
Е = B R r = 8133 -8,3736 = 68 102 кДж/моль.
З а д а ч а 5.24. Вычислить энергию активации Е по результа там лабораторных исследований кинетики окисления нефти.
Эксперименты проводились по методике, которая описана в за даче 5.20. Результаты опытов приведены в табл. 51. По результа там предыдущих исследований было определено, что показатель степени при парциальном давлении кислорода в формуле (5.11)
п = 1,0, а коэффициент, учитывающий влияние степени окислен ности нефти на скорость реакции, К = 5.
181
Таблица 51
1
Номер |
Скорость реак |
Температура |
Среднее парциаль |
Степень |
экспери |
ции, кг/(кг-с) |
кернодержателя. |
ное давление |
окисленности |
мента |
|
К |
кислорода, 105 Па |
нефти |
1 |
3,34-10—9 |
294,1 |
4,98 |
0,003 |
2 |
1,41 • 10- 8 |
306,8 |
5,01 |
0,008 |
3 |
3,01-10-8 |
318,6 |
4,70 |
0,013 |
4 |
1,49-10 - 7 |
337,1 |
4,93 |
0,029 |
5 |
3,35-10—7 |
352,0 |
4,80 |
0,056 |
6 |
3,41 • 10—7 |
357,3 |
4,83 |
0,098 |
При решении задачи необходимо замеренные скорости реакции привести к средним значениям парциального давления кислорода и степени окисленности нефти.
О т в е т : |
£ = 64 510 кДж/(моль-К). |
З а д а ч а |
5.25. В нагнетательную скважину с целью прогрева |
призабойной зоны пласта при инициировании процесса внутрипла-
стового горения закачивается горячий воздух. |
температура |
|||||||
Темп |
нагнетания |
воздуха |
^Воз = 6-104 |
м3/сут; |
||||
воздуха Т = |
300 °С; |
теплоемкость |
воздуха |
с^з = 1 |
кДж/(кг-°С); |
|||
плотность |
воздуха |
рв03 = 1,3 |
кг/м3; толщина пласта |
h = 10 м; |
||||
коэффициент |
охвата |
пласта |
по |
толщине |
т|о = 0,8; |
начальная |
пластовая температура Т 0 = 20 °С; теплопроводность пласта и ок
ружающих |
его |
пород |
А,пл = А,оп = |
8 Вт/(м-К); |
теплоемкость |
|
пласта и |
окружающих |
его |
пород |
спл = соп = 1,1 |
кДж/(кг-К); |
|
плотность |
пласта |
рпл = |
2600 |
кг/м3. |
|
|
Определить температуру пласта Т на расстоянии 0,5 м от нагне
тательной скважины через 1 день после начала прогрева. |
|
|||||
Р е ш е н и е . |
Температуру пласта можно рассчитать по за |
|||||
висимости, полученной А. Б. Золотухиным, |
|
|||||
0 =■ |
О*)2 |
(5.13) |
||||
12v(l — V и) |
||||||
|
|
|
||||
где |
|
|
|
|
|
|
|
|
hr\2 |
|
|
|
|
v = ? В О З Р в о з ^ в о з |
|
|||||
|
|
4л/АплП2 |
|
|||
и |
= |
Т - |
т0 |
|
||
т Б О З |
|
T Q |
|
|||
|
|
------ |
|
|||
0 |
|
^-ПЛt |
. |
|
||
|
|
|
|
|
СплРпл^ ^2
г — расстояние от оси нагнетательной скважины, м; t — время нагнетания агента, с.
182
Для условий задачи преобразуем формулу (5.13) к виду
Г - Г „ + ( Г „ — |
|
|
|
|
|
|
|||
Подставив |
исходные данные, получим |
|
|
|
|||||
г* = 0,5/(10-0,8) = 0,0625; |
|
|
|
|
|
||||
v = |
0,07-1,3-1 |
= 0,1132; |
|
|
|
|
|||
|
|
|
|
|
|
||||
|
4-3,14-Ю -0,88 -10—3 |
|
|
|
|
|
|||
|
8 -10~3-8,64-104 |
0,003776. |
|
|
|
|
|||
|
1,1-2600- 102 (0,8)2 |
|
|
|
|
||||
|
|
|
|
|
|
||||
Тогда |
|
|
|
|
|
|
|
||
Т = 2 0 + (3000- 20) [ 1- |
(0,0625)2 |
|
88 |
°С. |
|||||
12-0,1132-3,776-10= rj |
|||||||||
|
|
|
|
|
|
||||
З а д а ч а |
5.26. Вычислить время t, через которое температура |
||||||||
пласта достигнет 329,6 °С на расстоянии |
1 м от нагнетательной сква |
||||||||
жины при закачке в пласт воздуха. |
|
температура |
нагнетае |
||||||
Расход воздуха |
qB03 = |
1000 м3/сут; |
|||||||
мого воздуха |
Т воз = 500 сС; |
объемная |
теплоемкость |
нагнетаемого |
|||||
воздуха Своз = |
1,507 |
кДж/(м3-К); толщина |
пласта |
Д = |
10 м; на |
чальная пластовая температура Т 0 = 20 °С; объемная теплоемкость
пласта |
с'пл = 2860 |
кДж/(м3-К); |
теплопроводность |
пласта |
Хпл = |
|||||
= 2,4 Вт/(м- К). |
|
|
|
|
|
|
|
|||
|
О т в е т : |
t = 8,64-105 с = 10 сут. |
qB03, |
при котором тем |
||||||
|
З а д а ч а |
5.27. Найти расход |
воздуха |
|||||||
пература пласта будет равна 200 °С на расстоянии |
3,3 м |
через |
||||||||
100 сут после начала закачки горячего воздуха. |
|
|
|
|||||||
|
Толщина |
пласта |
h = 10 м; начальная |
пластовая температура |
||||||
Т = 20 °С; плотность пласта рпл = |
2600 кг/м3; теплоемкость пласта |
|||||||||
сПл=1,1 |
кДж/(кг-К); теплопроводность пласта Хпл = |
2,4 Вт/(м-К); |
||||||||
плотность воздуха |
рв03 = 1,3 кг/м3; теплоемкость |
воздуха С = |
||||||||
= |
1,1 |
кДж/(кг-К); |
температура |
нагнетаемого |
воздуха |
7+3 = |
||||
= |
500 |
°С. |
qB03 = |
5,787-10~2 |
м3/с = 5000 м3/сут. |
|
|
|||
|
О т в е т : |
|
|
§ 3. РАСЧЕТ ТЕХНОЛОГИЧЕСКИХ ПОКАЗАТЕЛЕЙ РАЗРАБОТКИ НЕФТЯНЫХ МЕСТОРОЖДЕНИЙ С ПРИМЕНЕНИЕМ ВЛАЖНОГО ВНУТРИПЛАСТОВОГО ГОРЕНИЯ И ЗАКАЧКИ В ПЛАСТ ПАРА.1 ЗАДАЧИ ДЛЯ УЧЕБНОЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКОЙ РАБОТЫ
При разработке нефтяных месторождений с применением] влаж ного внутрипластового горения (ВВГ) в пласт вместе с воздухом закачивается вода. В результате этого в пласте образуется высоко температурная зона, заполненная продуктами горения, азотом и во дяным паром. В передней части этой зоны выпариваются легкие фракции нефти и остается кокс, который затем сгорает.
183
Изменение отношения количества закачиваемой в пласт воды к количеству закачиваемого воздуха (водо-воздушного отношения) может приводить к существенному изменению температурного поля и скорости окислительной реакции. При проектировании ВВГ необходимо уметь рассчитывать его технологические параметры.
Срасчетом этих параметров связаны задачи 5.28К—5.31 К. Приближенный способ оценки показателей разработки нефтя
ного месторождения при нагнетании в него водяного пара дан в за даче 5.32К, вариант которой — задача 5.33К.
Задачи 5.34Н—5.38Н предназначены для учебной научно-иссле- довательной работы.
З а д а ч а 5.28К. Нефтяное месторождение, толщина продук тивного пласта которого составляет 17 м, насыщенное нефтью вяз костью р,н = ЮО мПа с, разрабатывается с самого начала, т. е. после ввода из промышленной разведки в разработку, с примене нием метода влажного внутри пластового горения (ВВГ).
Система разработки — однорядная, |
при этом длина одного эле |
|||||||||
мента |
(«полуэлемента») I = |
400 м, |
а |
ширина Ь0 = |
200 м (общий |
|||||
sc = 8-104 м2/скв). |
Коэффициент |
охвата |
пласта |
процессом |
ВВГ |
|||||
по толщине г)2т = 0,882, а по площади т]2п =0,9. |
|
|
|
|||||||
Абсолютная проницаемость пласта k = |
0,8 -10~12 м2; пористость |
|||||||||
т = |
0,22; плотность |
нефти |
р„ = |
0,85-103 |
кг/м3; |
плотность |
гор |
|||
ных |
пород |
рт = |
2,5-103 |
кг/м3; |
теплоемкость |
пород |
сТ = |
=1,257 кДж/(кг-°С); коэффициент теплопроводности пород кровли
иподошвы пласта Хт = 3-(Вт/(м-°С); температуропроводность пород хт = 9-10-7 м2-с = 0,07776 м2-сут. Среднее пластовое
давление в элементе р = 102 Па.
При осуществлении влажного горения в нагнетательную сква жину закачивается воздух с расходом 0,926 м3/с = 80-103 м3/сут при Т = 30 °С, так что в каждый элемент поступает 40-103 м3/сут воздуха. По данным лабораторных исследований процесса ВВГ
содержание кокса в породе пласта составляет |
zT = 28 |
кг/м3, теп |
|||||
лота |
сгорания |
кокса |
А = 25,14-103 кДж/кг; |
стехиометриче |
|||
ский |
коэффициент а = |
11,2 м3/кг. Пластовая |
температура Тпл = |
||||
- 30 °С. |
|
|
|
|
|
|
|
Начальная нефтенасыщенность |
пласта s„0 = |
0,93; |
насыщен |
||||
ность |
связанной |
водой |
sCB= 0,07. |
Относительные |
проницаемости |
пласта для нефти, воды и газов зависят линейно от соответствую щих насыщенностей. Остальные исходные данные, необходимые для расчета ВВГ, указаны в решении задачи.
Технология осуществления ВВГ в элементе системы разработки такова, что в начальной стадии ВВГ в пласте создаются высоко температурные зоны (паровое плато, зона горения) таким образом, что фронт конвекции хт(рис. 67) движется быстрее, чем фронт го рения, в 1,5 раза. Начальная стадия ВВГ продолжается до момента
времени |
t = tlt когда фронт горения достигнет расстояния Хф = |
||
= 50 м. |
Затем наступает стадия установившегося ВВГ, |
когда ско |
|
рость |
фронта конвекции vT = dxT/dt равна скорости |
фронта го |
|
рения |
Уф = d ^/d t. Эта последняя стадия продолжается до конца |
184
1 Z |
3 4 |
Рис. 67. Схема влажного внутрипластового горения (ВВГ). Зоны:
1 — /; 2 — горения / / ; 3 — парового плато I I I ; 4 — I V
разработки элемента пласта методом ВВГ, т. е. до момента времени, когда х? = I.
Требуется определить для рассматриваемого элемента пласта
температуру Т ш в паровом |
плато в зоне / / / , водо- и газонасыщен- |
ности в зонах I, III и IV, |
водовоздушное отношение Хвв на входе |
в пласт, дебиты нефти и воды, обводненности продукции в среднем за начальную стадию и в течение всего остального времени ВВГ.
На основе экспериментальных и теоретических исследований процесса принимается, что средняя температура в зоне/ / постоян ная и составляет Т% = 330 °С, а длина зоны / / Д£ = 1 м.
Р е ш е н и е . 1. Расчет начальной стадии создания высокотем пературной зоны в пласте.
При определении технологических показателей разработки пласта методом ВВГ, а также параметров самого процесса ВВГ применим одну из разновидностей так называемых зонных моделей, описание которых можно найти, например, в книге [5].
Согласно модели, используемой при решении настоящей задачи, поле температуры и насыщенностей в пласте при ВВГ разбивается на несколько характерных зон. На рис. 67 сплошной линией по казано распределение температуры в пласте, наблюдаемое, напри мер, в лабораторных экспериментах, а пунктирной — расчетное распределение температуры, принимаемое в «зонной» модели про-
185
цесса, на основе которой будут рассчитываться технологические показатели ВВГ
Всоответствии со схемой распределения веществ и температуры
впласте, показанной на рис. 67, в зоне I в пласте совместно дви жутся воздух и вода, в узкой зоне II происходит основная часть внутрипластовой окислительной реакции. В этой зоне находятся водяной пар, кислород, реагирующий с нефтяным остатком — кок сом, газообразные вещества, остающиеся или получающиеся в ре зультате реакции горения (азот, углекислый газ, окись углерода),
именуемые в дальнейшем как «газы горения». Водяной пар в зоне
/ / во многих случаях, по-видимому, бывает перегретым. |
|
||||||||
|
Зона I I I (паровое плато) содержит в основной своей части газы |
||||||||
горения и водяной пар, а в части, прилегающей к зоне IV, |
где про |
||||||||
исходят сложные |
процессы дистилляции |
легких |
фракций нефти |
||||||
и образование кокса,— смесь всех движущихся |
в |
пласте веществ. |
|||||||
Однако приближенно принимается, что |
в зоне |
I I I |
содержатся |
||||||
только газы горения и водяной пар. В зоне IV происходит совмест |
|||||||||
ное движение газов горения, нефти и воды. |
|
|
|
|
|||||
|
Граница между III и IV зонами считается фронтом конвекции |
||||||||
с координатой хт, |
а за |
координату фронта горения |
Хф, |
условно |
|||||
для |
удобства расчетов, |
принимается |
граница |
между зонами II |
|||||
и / / / . |
|
|
оф, исходя из стехиометрии |
||||||
|
Определим скорость фронта горения |
||||||||
горения. |
|
|
|
|
|
|
|
|
|
|
Имеем следующую формулу: |
|
|
|
|
|
|
||
|
Уф — Ц0г/-^В03> |
Уот = Явоз^г |
|
|
|
|
|
(5»14) |
|
где |
<7во3 — расход |
закачиваемого в пласт воздуха |
в |
нормальных |
|||||
условиях, м3; 5 — площадь поперечного |
сечения |
пласта; |
3 — |
объем воздуха, требующийся для выжигания кокса из 1 м3 пласта в целом, м3.
Соответственно для # в03 |
имеем |
формулу |
R воз — |
|
(5.15) |
где zx — содержание кокса |
в 1 м3 |
пласта, кг/м3; а — стехиомет |
рический коэффициент. |
процесса ВВГ, проведенным примени |
|
Согласно исследованиям |
тельно к конкретным условиям рассматриваемого месторождения, а = 11,2 м3/кг. По условию задачи zr = 28 кг/м3. Следовательно, Яв03 = 11,2-28 = 313,6 м3/м3.
Поскольку по условию задачи коэффициент охвата пласта про
цессом |
по |
площади т)2П = |
0,9, а по |
толщине |
т)2Т = |
0,882, |
то за |
|||
эффективную |
ширину пласта |
примем |
b = |
Ь0 0,9 = 200 0,9 = |
||||||
— 180 |
м, |
а |
за эффективную толщину пласта h = /г0 -0,882 = |
17X |
||||||
X0,882 = |
15 м. Тогда 5 = |
bh = |
180-15 = |
2700 м2. |
м3/с = |
40X |
||||
По |
формуле (5.14) с учетом того, |
что |
<7воэ= 0,46 |
|||||||
Х103 м3/сут |
|
|
|
|
|
|
|
|
||
|
|
|
40 -10s |
= 0,547-10~6 м/с = 0,04724 |
м/сут- |
|
||||
|
2700-313,6 0,864-10» |
|
||||||||
|
|
|
|
|
|
|
|
1S6
По условию задачи в начальной стадии ВВГ скорость фронта конвекции ут = 1,5 Уф = 1,5 0,04724 = 0,07086 м/сут.
Перейдем к определению температуры Т щ в зоне III, учитывая, что по условию задачи температура в зоне окислительной реакции II задана и остается в течение всего процесса ВВГ постоянной, равной Т* = 330 °С.
Вначале получим формулу для текущего ухода тепла из зоны I I I в кровлю и подошву пласта. В соответствии с условием задачи
и ее схематизацией будем считать, |
что температуры |
в зонах I и |
||
IV равны |
пластовой, так что |
= |
7 IV = Тпл. Будем |
также при |
ближенно |
принимать, что средняя температура в зоне I I I остается |
неизменной в течение всей начальной стадии процесса ВВГ. Если бы температура в зоне / была равной температуре в зоне III, т. е. при 0 < х < хтзначение Т = Тщ, то тепловые потери в кровлю и подошву пласта за счет теплопроводности можно определить по формуле
t
1 |
|
(* 2Хт (Гщ |
Т пл) vTbdx |
|
|
|
|
|
|
j |
|||||
* Ш“ |
1 ----- [п*, (< -*)]'2 |
• |
|
|
|
|
|
(5Лб) |
|||||||
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
Однако согласно схеме, представленной на рис. 67, температура |
|||||||||||||||
при х < |
Хф снижается |
скачком до пластового значения. Если же |
|||||||||||||
рассматривать только уход тепла из зоны III, то необходимо при |
|||||||||||||||
нимать, что уже при х < |
|
Хф температура снижается до пластовой. |
|||||||||||||
При этом нужно учитывать, что граница |
зон I I I |
и IV движется |
|||||||||||||
со скоростью |
ут, |
а |
граница |
зон II |
и |
I I I — со скоростью Уф. |
|||||||||
Учитывая |
сказанное, |
получаем |
|
|
|
|
|
|
|||||||
|
|
t 2А»т |
( 7*H |
I |
ТПл) vTbdx |
|
i |
2ХТ (Т HI — Тпл) ифbdx |
|||||||
?тШ = ^ |
[яхт (t —т)]12 |
|
|
[лхт (f — т)]1,2 |
|||||||||||
|
|
о |
|
|
|||||||||||
|
|
|
г |
(^т |
|
Уф) bdx |
|
л т", |
'г |
|
т |
|
|||
= |
2Arj-АТ п I f |
|
’ |
|
(5.17) |
||||||||||
(п*т)"2 |
'1 |
(f- т ) ' 2 |
4 Г - Г ш |
- Г п - |
|||||||||||
По |
условию |
задачи |
|
ут = |
1,5, |
|
Уф = |
const. |
Следовательно, |
||||||
ут— Уф = |
0,5 уф. В этом случае |
|
|
|
|
|
|
||||||||
йтП1 = |
2ХтД'Гц10,5иф6 |
Г |
dx |
|
|
|
|
|
|
||||||
|
|
(лхт) ь |
|
оS |
( t - |
X)12 |
|
|
|
|
|
||||
|
|
|
|
|
t |
\1:2 |
|
|
|
|
|
|
(5.18) |
||
= 2ЯтДТи1&Цф ( ------^ |
|
|
|
|
|
|
|
||||||||
В более общем же случае, когда существует другое, но постоян |
|||||||||||||||
ное отношение |
ут/Уф, |
имеем |
|
|
|
|
|
|
|
||||||
? Т Ш = г ^ Д Г ш Н ( |
V ф |
— 1 ) ( — — )'я |
|
( 5. 19) |
|||||||||||
|
|
|
|
|
\ |
|
|
/ |
\ ЯХт |
/ |
|
|
|
187
Всоответствии с условием задачи зона окислительной реакции
сне изменяющейся во времени шириной Д£ и температурой Т * движется в пласте также со скоростью Оф.
в |
Обозначим через |
t* время, за которое фронт горения проходит |
|
пласте расстояния |
Д£. В этом случае t%= А|/«ф. |
пласта |
|
в |
Выражение для текущих тепловых потерь из зоны II |
||
кровлю и подошву |
получаем аналогично формуле (5.17). |
Имеем |
„ |
_ 2Хт6иф (Т% |
Тпл) |
|
^11 |
/ |
N|/2 |
|
|
(лхт)1'2 |
|
|
_ |
2\тЬифАТ* |
г' |
|
|
(лхт)Ь2 |
Lo |
|
- |
(п*т) |
[2 |
+ |
|
|
|
г t |
|
dx |
t-U |
dx |
Г |
|
Г |
||
) |
V - ' ) 112 ) |
(t-т;)112 |
||
Lo |
|
|
t* |
1 Н и1 |
------ |
— |
1 ьо |
||
- т |
|
1 |
|
|
|
|
U |
1 |
|
|
|
- 2 |
V t - U |
]• |
При t > t* имеем
ЧП1 = 4 )ф > * Ь Т ,( - ± - ') . |
(5.20) |
V лхт /
Поскольку согласно условию задачи из зон I и IV тепло не уходит в кровлю и подошву пласта (TY = T1V = Тпл), то при / > полные тепловые потери в кровлю и подошву qT можно оп ределить следующим образом:
= ^7тц+ <7т1ц = |
2Я,тЬиф j^A r in |
^ -f- |
+ 2ДТ, ( ———У 2] . |
(5.21) |
|
V лхт / |
J |
|
Тепловые потери в кровлю и подошву пласта QT, накопленные ко времени t , получают интегрированием (5.21). Имеем
Q T = Q T II + Q T I I I — |
4Хт6Д7'.1.Оф ( — 5— ^ |
t + |
|
V лхт / |
|
+ -?-М щ ,А Г „,- |
^ v . |
(5.22) |
3(лхт)
Всоответствии со схемой процесса ВВГ, представленной на рис. 67, для количества тепла Qnjl, накопленного к моменту вре мени t в самом пласте, имеем следующую формулу:
<2пл = Стрт(1 —m) bhA Т ш (от—оф) t +
н - m C epaSni ATu l(vT — Оф) bht - f
+ |
m (срРп + сгРг) (1 — 5 щ ) Ы гА Т щ (ит — Цф) t -}- |
|
+ |
СтРт (1 — т) Д Т *& Л Д £ , |
(5.23) |
188
где sn |
и sm — соответственно насыщенности |
водой |
в |
зонах II |
и III. |
Количество тепла Q*, выделившегося в |
пласте |
к |
моменту |
времени t в результате реакции горения, выражается следующим образом:
Q* ~ z TAbh(v^t — А£). |
(5.24) |
Согласно балансу тепла имеем |
|
Q* = QT+ Q™. |
(5.25) |
Как видно из структуры формул, с помощью которых опреде |
|
ляли Q*, QT и <2пл, члены, входящие в указанные выше |
величины, |
по-разному зависят от времени t. Это означает, что либо в случае постоянства АТ* и ATlU в зоне I I I будет изменяться насыщенность sIH пласта жидкой водной фазой, либо в случае существования в зоне I I I только перегретого пара АТщ не останется постоянной. Однако при сравнительно небольшом продвижении в пласте фронта горения АТ 1п может изменяться в сравнительно небольших преде лах, так что приближенно в этих случаях можно считать, что АТ т
равна среднему значению ATjn в зоне I I I и постоянна в течение рас
сматриваемого времени протекания ВВГ. |
с зоны |
IV, |
|
Расчет |
технологических показателей ВВГ начнем |
||
в которой |
необходимо определить насыщенности зоны |
нефтью, |
га |
зами горения и водой и, как наиболее важный показатель, скорость
фильтрации воды uBiv в зоне 1У-
Для этого рассмотрим совместное движение газов горения, нефти и воды в зоне IV
Обозначая насыщенность зоны IV газами горения srIV водой
sIV, |
нефтью sHiv, имеем |
|
|
|
||
|
•Sriv _1_SHIV “Г siv = |
1* |
|
|
(5.26) |
|
Для скоростей фильтрации соответствующих веществ известны |
||||||
следующие выражения: |
|
|
|
|||
|
Vr — |
kkv (sr) |
dp |
_ |
kkH(sH) |
dp |
|
-------------------------1 , |
VH — |
■ |
, |
||
|
|
(Xr |
dx |
|
|i„ |
dx |
|
Vв |
kkB(sB) |
dp |
|
|
(5.27) |
|
Ив |
dx |
|
|
||
|
|
|
|
|
||
где |
kr, |
kn и /гв — относительные |
проницаемости соответственно |
для газа, нефти и воды.
Примем, исходя из условия задачи, что относительные прони цаемости в зоне IV для газов, нефти и воды линейно зависят от со ответствующих насыщенностей, а именно
К |
Sp — Spo |
kB |
S — |
S e e |
(5.28) |
||
1 |
|
||||||
1 |
SCB |
SCB |
|||||
|
|
|
|||||
где sro — насыщенность пористой |
среды газом, при которой про |
||||||
ницаемость для газа равна нулю; |
sCB— насыщенность пласта свя |
||||||
занной водой; s* — водогазонасыщенность, при |
которой проницае |
||||||
мость для |
нефти становится равной нулю. |
|
|
|
189
Из зоны I I I |
в зону IV входит количество нефти, |
пропорцио |
|||||
нальное скорости |
ыт> а также первоначальному содержанию нефти |
||||||
в пласте за вычетом количества сгоревшего кокса, так что |
|||||||
и„ = ^ms„0---- vTt |
|
|
|
(5.29) |
|||
где рк — плотность |
кокса. |
0,95* 103 кг/м3 |
и |
учитывая, |
|||
Вычислим |
UHIV- |
Принимая рк = |
|||||
что т = 0,22; sH0 = |
0,93; sCB— 0,07; |
vr = 1,5 иф = |
0,07086 м/сут, |
||||
получим |
|
|
|
|
|
|
|
uHiV = |
Г0,22 -0,93---------------- 0,07086 - 0,0124 м/сут. |
|
|||||
|
V |
|
|
0,95-10_3 J |
|
J |
|
Составим |
отношение |
|
|
|
|||
и г |
________ Р н |
(s rIV — s r o ) |
|
|
(5.30) |
||
Z>HIV |
Н-г (s* — SIV — Sriv) |
|
|
||||
|
|
|
Будем полагать, согласно закону идеальных газов, при среднем дав лении в пласте р = 107 Па, что
yriv — ^го Ро 40-10М05 = 0,1481 м/сут.
р180-15-107
При sr0 = 0,05 |
значение |
s* = 0,95. |
Учитывая, что в зоне IV |
||
рн = |
ЮО |
мПа-с, |
рг = 0,02 |
мПа-с, в |
соответствии с формулой |
(5.30) |
получаем |
|
|
|
|
0,1481_________ 100 (SrIV — 0,05) |
|
||||
0,0124 |
0,02(0,95 — sjy — sriy) |
|
|||
Выпишем отношение |
|
|
|||
a BIV _ |
Рн (SIV — $св) |
|
(5.31) |
||
«HIV |
MB (S* — «IV — Sriv) |
|
|||
|
|
Из этого отношения получаем выражение для uBlv:
l>uiv |
1,24 (Siv — 0,07) |
(5.32) |
|
0,9 — siv |
|||
|
|
Таким образом, на основе теории установившегося трехфазного течения в зоне IV получены уравнения (5.30) и (5.31) для опреде ления трех неизвестных uBlv, sIV и sriv. Недостающее третье уравнение получаем на основе выражения для скорости конвек
тивного переноса тепла на границе зон II I |
и IV в виде |
|
||
vT |
_______________СгРг^Г ~Т Снр нЦц 4~ СврвРв_______________ |
(5.33) |
||
стрт (1 — т) + т [Сцр„ (1 — S — sr) + cBpDs + crprsr] |
||||
|
|
Систему уравнений (5.30), (5.31) и (5.33) будем решать методом последовательных приближений.
190