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

1275

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

парциальное давление

кислорода, Па;

Е — энергия активации,

кДж/( моль-К);

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

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