Пределы
выкипания, ˚С
|
Средняя
t, ˚С
|
Молекулярная
масса, Мi
|
%
выход на нефть
|
Массовое
кол-во, кг/ч
|
Массовая
доля, Yi
|
Число
молей, кмоль/ч
|
Мольная
доля, Yi’
|
120-130˚С
130-140˚С 140-160˚С 160-180˚С
|
125
135 150 170
|
113,1
118,7 127,5 139,9
|
1,5
2,0 3,5 3,5
|
9682,6
12910,2 22592,8 22592,8
|
0,14 0,20 0,33 0,33
|
85,6
108,8 177,2 161,5
|
0,16
0,20 0,33 0,31
|
Сумма
|
-
|
-
|
19,5
|
677778,4
|
1
|
533,1
|
1
|
Температура паров
на входе в конденсатор - холодильник определяется методом постепенного
приближения по уровню изотермы паровой фазы:
Σ
= 1 , (2.5)
где Ki - константы фазового равновесия при давлении П = 103000 Па и
температуре начала конденсации, определяется по формуле [4,c. 124]:
где Pi - парциальное давление компонентов при температуре Тн.к.
Примем Тн.к. = 428
К и по диаграмме Кокса [ 5, с.45,46 ] найдем давления
насыщенных паров всех компонентов:
1 = 3,7 кг/см2 = 360500 Па;
Р2 =
2,95 кг/см2 = 288400 Па;
Р3 = 0,7
кг/см2 = 69010 Па;
Р4 = 0,3
кг/см2 = 29870 Па.
Аналогично найдем:
К2 = 2,8; К3 = 1,5; К4 = 0,44;
Подставив в
уравнение изотермы паровой фазы все известные величины, получим:
Полученная сумма
несколько больше 1, однако расхождение составляет 0.06%, поэтому пересчета не
делаем.
Следовательно,
температура Тн.к. = 428 К.
Температуру конца
конденсации определяем методом постепенного приближения по уравнению изотермы
жидкой фазы [1.с.226]:
Σxi’*Ki = 1 , (2.7)
где Ki - константа фазового равновесия , при давлении П = 103000 Па и
температуре конца конденсации Тк.к., определяемая по формуле [ 4.с. 124]:
где Pi - давление насыщенных
паров компонентов, при температуре Тк.к. Примем Тк.к. = 403 К и но диаграмме Кокса найдем давления
насыщенных паров всех компонентов:
1 = 2,6 кг/см2 = 254974,2 Па;
Р2 = 1,9
кг/см2 = 186327,3 Па;
Р3 = 0,8
кг/см2 = 78453,6 Па;
Р4 =
0,15 кг/см2 = 14710,1 Па;
Аналогично найдем :
К2 = 1,81, К3= 0,76, К4= 0,143
Подставив в
уравнение изотермы жидкой фазы все известные величины, получим
Полученная сумма
несколько меньше 1, однако расхождение составляет примерно 0,98% , поэтому
пересчета не делаем.
Следовательно,
температура Тк.к. = 403 К.
2.4 Тепловая нагрузка
конденсатора
2.4.1 Тепловая
нагрузка зоны конденсации паров
Тепловую нагрузку
определяем по формуле [4, с.106]:
1 = G (JТн.к. - iТн.к.), (2.9)
где G - количество конденсируемого
бензина, кг/ч;
JТн.к.
- энтальпия паров бензина, при Тн.к. и его плотности ρ293277 = 0,7526 [2.с.83];
iТн.к. - энтальпия жидкого
бензина, при Тн.к. и ρ293277.
Пользуясь
приложением 2; 3 [ 1,c.328;332], интерполяцией находим JТн.к.
и iТн.к. JТн.к = 648,7 кДж/кг ; iТн.к. = 350,6 кДж/кг .
Подставив известные
величины в формулу 2.9 , получим:
1= 67778,4 ( 648,7 - 350,6 ) = 20204741,04 кДж/ч = 5612,4 кВт.
2.4.2 Тепловая нагрузка зоны
охлаждения конденсата
Тепловую нагрузку
определяем по формуле [4,с.133]:
воздушный
охлаждение конденсация бензин
Q2 = G (iТн.к.- iТохл), (2.10)
Тохл = 80°С = 353 К
и ρ293277 = 0,7526
Пользуясь
Приложением 3 [ 4, с.332 ], интерполяцией находим iТохл: iТохл = 167,5 кДж/кг.
Подставив известные
величины в формулу 2.10, получим:
2= 67778,4 ( 350,6 - 167,5) = 12410225,04кДж/ч = 3447,3 кВт.
Общая тепловая
нагрузка конденсатора:
= Q1 + Q2
Q = 5612,4 + 3447,3 = 9059,7 кВт.
2.5 Определение
расхода воздуха
Не учитывая потерь
тепла в окружающую среду и принимая температуру воздуха на выходе из аппарата T1" = 328 К, определим расход воздуха.
Уравнение теплового
баланса конденсатора [4, c.134]:
Q = WСв(T1" - T1'),
(2.12)
откуда ,
(2.13)
где W - расход воздуха,
кг/ч;
Q - тепловая нагрузка аппарата, кВт;
T1' = 291 К - начальная
температура воздуха;
Св = 1,005 кДж
/(кг*К) - теплоемкость воздуха [4, e. 107],
Получим:
кг/ч
2.6 Определение
конечной температуры воздуха
Тепловой баланс
первой зоны [4, с.135] :
1
= W ( T1"-T1 )СВ
, (2.14)
откуда:
(2.15)
.
2.7 Расчет зоны
конденсации
2.7.1 Коэффициенты
теплоотдачи со стороны бензина
Коэффициент
теплоотдачи от конденсирующих паров бензина к наружной поверхности
горизонтальной трубки аппарата рассчитываем по формуле [4,с.137]:
где r - теплота конденсации углеводородных паров. Дж/кг;
ρ - плотность
конденсата, кг/м3;
λ
- коэффициент теплопроводности конденсата, Вт/(м • к) ;
μ
- динамическая вязкость конденсата. Па • с ;
∆Т - разность
температур конденсирующего пара и стенки, К;
dН
= 0,028 м - наружный диаметр трубок аппарата. dН
берется по чертежу.
Значения ρ,λ и μ берутся при средней температуре пограничного слоя конденсата Тm,
r - при средней температуре конденсации Ts.
= (Tн.к. + Тк.к.)/2 ; [4,с.138] (2.17) Ts = (428 + 403)/2 = 415,5
К.
Принимаем, что
температура стенки трубы Тw = 326 К. Тогда средняя температура пограничного слоя
конденсата:
= 0,5(Ts + Tw) (2.18)
Tm = 0.5 (415,5 + 326) = 371 К.
Коэффициент
теплопроводности находим по формуле [4. c 138] :
, (2.19)
где р288288
- относительная плотность конденсата.
Ρ288288 = ρ420 + 5α, (2.20)
где α - средняя температурная поправка на 1К определяемая по таблице.
Зная плотность
конденсата в конце первой зоны р293277 = 0,7526 и
пользуясь таблицей средних температурных поправок находим по
формуле Менделеева [1, c.20]:
Tm277
= p293277
- α (Tm - 293); (2.21)
P367277
= 0,7526 - 0.000874 (371 - 293) = 0,684.
Значение r находим, как разность энтальпий конденсата в паровой и жидкой
фазах, при:
= 415,5 К и р277293
= 0,7526;
г = Jтs
- i тs
[4, с.138] (2.22)
Пользуясь Приложениями
2;3 [4, с.328: 332], интерполяцией находим
Jts
и iTs: J393
= 1314,8 кДж/кг; i379
= 1127,9 кДж/кг.
Подставив полученные
величины в формулу 2.22 получим:
= 1314,8- 1127,9 =
186,9 кДж/кг.
Зная кинематические
вязкости бензина при 293 К и 323 К, находим значение νTm, пользуясь формулой Гросса [4, с 99]:
, (2.23)
где n - коэффициент.
В нашем случае:
ν1
= ν293 = 1,105*10-6 м2/с , ν2
= ν323 = 0,78*10-6 м2/с [2, c. 216]
T1
= 293 К; Т2 = 323 К.
Решив формулу
относительно n при известных ν1,
и ν2, получим:
;
Тогда кинематическая
вязкость для бензина, при 367 К определится из уравнения:
(2.25)
;
Откуда ν371 = 0,875*10-6 м2/с.
Для перевода
кинематической вязкости в динамическую воспользуемся формулой [1, с.23]:
μ = 0,875*10-6
* 684 = 5,98*10-4 Па·с;
Величина ∆T будет равна:
∆Т = Ts - Tw (2.39)
∆Т = 379 - 326 =
53 K.
Подставляя найденные
величины в формулу 2.16 получим:
Вт/м2*К
В соответствии с
конструкцией аппарата [7], среднее число трубок, расположенных в одном
вертикальном ряду, равно n
= 15. Поэтому коэффициент теплоотдачи α1
, от конденсирующего пара к поверхности пучка горизонтальных трубок найдем по
формуле [4, с.138]:
α1
= εn · α1’ , (2.40)
где En - усредненный поправочный коэффициент, при различном размещении
труб в пучке, - определяется по графику [8, с.305].
В нашем случае n = 15 получим:
α1
= 0,6 * 1447 = 868,2 Вт/м2*К
2.7.2 Коэффициент теплоотдачи со стороны воздуха в случае применения гладких труб
В целях правильного
выбора расчетной формулы для α2,
следует определить значения критерия Рейнольдса для воздуха, при поперечном
обтекании шахматного пучка труб аппарата.
Примем, что фронтальное
к потоку воздуха сечение аппарата L × В = 4 × 4 м2 [4 ,с.110] с шагом труб по ширине пучка S1 = 0,052 м. Шаг труб по глубине пучка S2
найдем следующим образом:
(2.41)
мм
Определим число n труб в одном горизонтальном ряду пучка из формулы:
Примем число рядов труб
по вертикали одной секции nв
= 8.
Определим площадь
сжатого (наименьшего) сечения в пучке труб, через которое проходит воздух [4,c.111]:
= L(B - nd3); (2.44)
Fc = 4 (4 - 76 · 0,028) = 7.5 м2
Скорость воздушного
потока в сжатом сечении
(2.45)
где Vd - действительный
секундный расход воздуха (берется из паспорта данных на вентилятор), м3/с
Средняя температура
воздуха в конце первой зоны [4, c.139]
По таблице [4,
с.107] находим интерполяцией кинематическую вязкость воздуха при его средней
температуре, а также коэффициент теплопроводимости воздуха при его средней
температуре Тср: ν = 16.55*10-6 м2; λ =
0,0.0268 Вт/м*К. Теперь
определим величину критерия Рейнольдса
Re = (W0*d3)/ν [4, c. 111] ; (2.47)
Re
= =
14380.7
Коэффициент теплоотдачи α
определим из уравнения [9, с.191], справедливого при:
= 2*102 -
0.2*106
Получим:
α2 =
0,37*(λ/d3)*Еат*Re0.6
, (2.48)
где Еат = 1 -
поправочный коэффициент, учитывающий угол атаки, принимается по таблице [9, с.
190].
α2
= 0,37*(0,0268/0,028)*14380.70,6*1 = 110.6 Вт/м2*К.
2.7.3 Коэффициент теплопередачи для пучка гладких труб
Для биметаллических труб
(латунь-алюминий) и загрязненной поверхности теплообмена (внутренней и
наружной) этот коэффициент определяется пи формуле [4,с.111]:
, (2.49)
где -
тепловое сопротивление внутреннего слоя загрязнения (принимаем для бензина
равным 0.00035 (м2.к)/Вт (см.Приложение 5)
(м2. К) / Вт
- тепловое сопротивление латунной стенки т трубы при δ
= 0,002 м. и λ
= 91,9 Вт / (м.К) [4 .с. 140];
=
0,0015/205 = 0,000073 (м2. К) / Вт - тепловое сопротивление
алюминиевого слоя трубы, при δ = 0,0015
м и λ = 205 Вт/(м.К) [4,c.112];
- тепловое
сопротивление наружного слоя загрязнения - выбираем в пределах
0.00017 - 0.00086, для дальнейшего расчета эта величина принята равной 0,00060
(м2К)/Вт. [4 ,с. 112].
Подставив эти
значения в формулу 2.38, получим
Вт/(м2
К) .
2.7.4 Коэффициент теплоотдачи αк , при поперечном
обтекании воздухом пучка оребренных труб
При спиральном
оребрении труб, расположенных в шахматном порядке, для определения коэффициента
теплоотдачи [в Вт/(м2.К)], воспользуемся формулой [4, c.113]
, (2.50)
где λ = 0,02619Bt/(m.K) -
коэффициент теплопроводности воздуха, при его средней температуре [4,с. 107];
ρВ - плотность воздуха при
Тср, кг/м3, ρВ = 1,17727 кг/м3
;
W0’ - скорость воздуха в сжатом (узком) сечении одного ряда
труб оребренного пучка, м/с;
μ - динамическая вязкость воздуха, при Тср, Па • с;
Pr - критерий Прандтля, при Тср Pr = 0.7123[4.с. 107];
δР - средняя толщина
ребра, м.
Величину Wo' определяем по
формуле [4,с. 113]:
,
(2.51)
где WН - скорость набегающего воздушного потока, при входе в трубный
пучок, то есть в свободном сечении перед секциями оребренных труб,
δ = Sn/d3 , (2.52)
где Sn - поперечный шаг
оребренных труб, принятый ранее равным 0,052 м
δ = 0,052/0,028 =
1,86;
Нр = 0,0105 м
высота ребра;
Sp = 0,0035 м - шаг ребер;
Скорость
набегающего воздушного потока [4,с. 113]:
, (2.53)
где VД - действительный секундный
расход воздуха, м/с; Fcb = L×B = 4*4 = 16m2
- фронтальное
к потоку воздуха сечение аппарата.
Таким образом: WН = 64/16 = 4 м/с
Средняя толщина
ребра:
δр = (δ1 + δ2) / 2, (2.54)
где δ1 = 0,0006 м - толщина ребра в его вершине;
Имеем:
δр = (0,0006 + 0,0011) / 2 = 0,00085 м
Скорость воздушного
потока в сжатом сечении:
м/с.
Динамическая вязкость
воздуха при Тср:
μ = νρВ
[4, c.114] , (2.55) где ν
- кинематическая вязкость воздуха;
ν = 15,296*10-6
м2/с;
μ = 16,55*10-6*1,1442
= 1,89*10-5 Па·с.
Подставив в формулу 2.39
значение всех величин, получим:
Коэффициент теплоотдачи αК
сильно зависит от диаметра трубы d3
, несколько меньше от толщины ребра δр,
и почти не зависит от диаметра ребра d4.
2.7.5 Приведенный коэффициент
теплоотдачи αпр со стороны воздуха в
случае пучка оребренных труб
Приведенный
коэффициент теплоотдачи αпр учитывает конвективный
теплообмен между оребренной поверхностью и потоком воздуха и передачу тепл
теплопроводностью через металл ребер. Его величину необходимо знать, чтобы определить
коэффициент теплопередачи К.
В литературе
отсутствуют надежные данные , позволяющие подсчитать αпр для случая ребристой спиральной поверхности. В виду малого шага
спирали определяем приведенный коэффициент теплоотдачи по формуле для дисковых
(круглых) ребер [4,с. 115]
, (2.56)
где Fp - поверхность ребер, приходящихся на 1 м длины трубы, м2/м;
Е - коэффициент
эффективности ребра, учитывающий понижение его температуры по мере удаления от
основания [4,с. 115];
ε∆
- коэффициент, учитывающий трапециевидную форму сечения ребра (рисунок 2.7
[1,сП5]);
ψ
= 0,85 - экспериментально найденный коэффициент учитывающий неравномерность
теплоотдачи по поверхности ребра;
β3 =
(δ/λ)З.Н - тепловое сопротивление загрязнения наружной поверхности трубы,
принимаемое (для сравнения) равным тепловому сопротивлению наружного
загрязнения поверхности гладких труб (смотри выше). Находим поверхность ребер,
приходящуюся на 1 метр длины трубы
, (2.57)
где Х = 286 - число
спиральных витков ребер, приходящихся на 1 м длины трубы.
м2/м.
Определяем наружную
поверхность участков трубы между ребрам приходящуюся на 1 м длины трубы [4 ,с.
116]:
Fтр = πd3(1-хδ2) , (2.58) где δ2 = 0.0011 м - ширина ребер у основания Fтр = 3.14 • 0,028 ( 1- 286 • 0,0011
) - 0,06 м2/м. Полная наружная поверхность 1м оребренной трубы будет
равна: Fn = Fр + Fтр [4,с. 116]; (2.59) Fn = 0,761 + 0,06 = 0,821
м2/м.
Предварительно
вычисляем соотношения, необходимые для пользования рисунками 2.6 и 2.7 [4,с.
115]:
4
/ d3 = 49 / 28 = 1.75;
√(δ1/δ2) = √(0.6/1.1) = 0.738.
Чтобы учесть наружное загрязнение труб, необходимо вычислить
подкоренные выражения аргументов графических зависимостей (смотри рисунки 2.6 и
2.7 [4,с.115] ) и поделить их на величину ( 1 + β3
ψ Lк ):
, (2.60)
где αК = 205 Вт/ (м • к) - коэффициент теплопроводности алюминиевого
ребра:
Тогда по рисункам 2.6 и
2.7 [4, c.115] Е = 0,96; ε∆
= 1,02
По формуле 2.45
определяем приведенный коэффициент теплоотдачи:
.
2.7.6 Коэффициент теплопередачи для пучка оребренных труб
Ведем расчет на единицу
гладкой поверхности трубы по формуле:
, (2.61)
где Fст
- поверхность гладкой трубы (по наружному диаметру), приходящаяся на
1м ее длины :
т = πd3 • 1 [4,с. 116];
(2.62)
Fст = 3,14 •
0,028 • 1 = 0,088 м2/м.
Все остальные
величины и обозначения - прежние. Получим:
Вт/(м2 *К),
Следовательно, при
прочих равных условиях оребрения гладкой трубы со стороны воздуха приводит к
значительному увеличению коэффициента теплопередачи ( в 314,5 / 92,78 = 3,4
раза ).
2.7.7 Средний температурный напор зоны
конденсации и ее поверхность
При многоходовом
потоке теплоносителя в трубном пространстве аппарата и одноходовом потоке
теплоносителя в межтрубном пространстве определяется по методу Белоконя [4,с.
112]:
, (2.63)
где ∆Тср -
средний температурный напор первой зоны, К;
Тmax1 , Tmin1 - соответственно
большая и меньшая разность температур, определяется по формулам:
Тmax1 = Ө1 + 0,5∆TI ; (2.64)
Tmin1 = Ө1 -
0,5∆TI , (2.65)
где Ө -
разность среднеарифметических температур горячего и холодного теплоносителей
, [4, c.112] (2.66)
Расчитаем ∆ТI по формуле:
∆ТI =
√(∆Т1 + ∆Т2)2 - 4Р∆Т1∆Т2
, (2.67)
где ∆Т =
Т н.к - Т к.к - перепад температур в горячем потоке ; (2.68)
∆Т2 = T1 - T'1’
- перепад температур в холодном потоке ; (2.69)
Р - индекс
противоточности. В нашем случае Р = 0,98. Имеем:
∆Т1 =
408 - 354 = 54 К
∆Т2 =
319 - 298 = 21 К
∆Т =
√(54 + 21)2 - 4*0,98*54*21 = 34,4
К
Tmax1
= 72,5 + 0.5 * 34,34=90 К
Тmin1 =72,5-0,5*34,4 = 55 К
К.
Поверхность теплообмена
в первой зоне определим по формуле:
; (2.70)
м2 .
2.3.8 Расчет зоны охлаждения
Средний температурный
напор зоны охлаждения определяем по методу Белоконя [4,с. 112]:
где ∆Тср2
- средний температурный напор второй зоны, К;
Т max2, Тmin2
- большая и меньшая разность температур, определяется
по формулам:
Тmax2 = Ө2 + 0,5∆TII
; (2.72)
Tmin2
= Ө2 - 0,5∆TII
, (2.73)
где Ө2 -
разность среднеарифметических температур горячего и холодного теплоносителей.
, (2.74)
где Тв - конечная
температура бензина (на выходе), К;
Характеристическую
разность температур ∆TII
рассчитаем по формуле (2.75), где ∆T1
= Т кк - Тв - перепад температур в горячем потоке; (2.76)
∆T2 = T1"
- Т1 - перепад температур в холодном потоке; (2.77)
Индекс противоточности
Р= 0.98 [4,с. 112] оставляем прежним:
∆Т1 =
354-333 = 21 К;
∆Т2 =
328 -319 = 9 К;
∆ТII =
√(21+9)2 - 4*0,98*21*9 = 12,6 К;
К;
Tmax2
=20+0,5*12,6= 26,3 К;
Тmin2 = 20-0,5*12,6 = 13,7 К;
К .
Принимаем коэффициент
теплоотдачи:
1
= 150 Вт/(м2. к),
Тогда поверхность
теплообмена второй зоны определяем по формуле:
; (2.78)
м2 .
2.7.9 Компановка и аэродинамическое сопротивление пучка труб
Находим общую
поверхность теплообмена аппарата:
= F1 + F2 ;
(2.79)
F
= 1134+369=1503 м2 .
Количество труб равно
[4,с. 117]:
= F / Fo, (2.80)
где Fo = 3,14*0,028*8 = 0,704
м2 - поверхность теплообмена одной трубы
N = 590 / 0,704 =
2135 шт.
Определяем число
труб для одного хода бензина, при скорости движения W = 0,74м/с [4,с. 109], по формуле
[4, с. 117]:
; (2.81)
трубы .
Устанавливаем 4
аппарата с общим числом труб в каждом - 564, Пучок распределяем на 3 секции по
188 труб. Nд = 564*4=2256 труб. Аэродинамическое сопротивление пучка труб (в
Па) определяем по формуле:
, (2.82)
где ρВ
= 1,18505кг/м3 - плотность воздуха при его начальной температуре;
W0’ = 10,6 м/c - скорость воздуха в сжатом (узком) сечении оребренного
пучка труб;
nВ = 8 - число
горизонтальных рядов труб в пучке (по вертикали);
d3 = 0,028 м - наружный
диаметр трубы .
Критерий
Рейнольдса, отнесенный к диаметру труб d3 [4,c. 118]:
; (2,83)
.
Подставляя указанные
величины в формулу 2.70, получим:
Па
Принимаем коэффициент
оребрения равным 9 [3, c.263],
тогда истинное количество труб равно:
д
= 564*4 = 2256 шт.
Выбираем аппарат по ГОСТ
26-02-1522-77 : число труб в аппарате 188 шт., поверхность теплообмена 376 м2
, длина труб 8000 мм. АВО - малопоточного типа.
Для проектируемого
аппарата выбираем осевой вентилятор ГАЦ-28-8М2 с регулируемым углом установки
лопастей.
Основные технические
данные вентилятора
диаметр рабочего колеса
-2.8 м;
число лопастей -8шт.
частота вращения
428об./мин.
мощность приводного
электродвигателя -37кВт;
номинальный напор
-320Па;
диапазон рабочих температур
-от -60 до 60 0С.
Список использованных
источников
1.
Молоканов Ю.К; Процессы
и аппараты нефтегазопереработки. М.: Химия. 1987, с.368.
2.
Нефти СССР (справочник) Том I. Нефти северных районов европейской части СССР и Урала.
Издательство Химия, М., 1971г. с.504.
3.
Рудин М. Г., Драбкин А. Е. Краткий справочник нефтепереработчика. - Л.: Химия,
1980. - 328с.:ил.
.
Кузнецов А.А; Кагерманов С.И.; Судаков Е.Н; Расчеты процессов и аппаратов нефтеперерабатывающей
промышленности. Л:
Химия. 1974.с.344.
5.
Эмирджанов Р.Т. Основы
технологических расчетов в нефтепереработке. М - Л.: Химия. 1965г. с.544.
6.
Эмирджанов Р.Т. Примеры
расчетов нефтезаводских процессов и аппаратов. Баку, Азнефтеиздат.,
1957г. с.404.
7.
Павлов. Новое нефтяное оборудование, М., Гостоптехиздат, 1961г. с. 155.
.
Касаткин А.Г. Основные процессы и аппараты химической технологии. М., Химия
1973г. с.784.