Файл: Тюряев И.Я. Теоретические основы получения бутадиена и изопрена методами дегидрирования.pdf

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

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

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

Добавлен: 11.04.2024

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

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

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

И з б и р а т е л ь н о с т ь . Поскольку величина WIF, а также изменение температуры по реактору известны, избирательность лег­ ко можно рассчитать по уравнениям (IX,4) и (IX,8). Для этого в простейшем случае по средней температуре в реакторе по уравне­ нию (IX,5) находится константа скорости крекинга и далее по урав­ нению (IX,4) определяется величина хк — доля бутана, разлагаю-

Рнс. 50. Изотермы и адиабаты дегидрирования бутана по уравнениям (IX ,16) и (IX,8а).

щегося в результате крекинга. Согласно определению, избиратель­ ность

З а у г л е р о ж е н н о с т ь к а т а л и з а т о р а . По величи­ не W/F и Q определяется время пребывания (мин) катализатора в реакторе:

_ W/F • 60 58Q '

Как в случае расчета избирательности, по средней температуре в реакторе рассчитывается величина А г (IX,7) и далее по уравнению (IX,6) определяется зауглероженность в процентах к весу катали­ затора.

С о с т а в к о н т а к т н о г о г а з а . Объем получаемого кон­ тактного газа рассчитывается по коэффициенту увеличения объема газа а = 1 + 1,1 Х к. Содержание С4Н8 + С4Н0 в газе равно Ху Іа, а водорода — Xyß/a. Количество продуктов, получаемых в резуль­ тате крекинга, рассчитывается по величине хк и уравнению крекин­ га:

С4Н10

0,62СН4

0,13С2Н4 +

0,44С2І1в +

 

-f- 0,04C3He +

0,2С3Н8 +

0,ЗС5

226

Количество бутадиена с достаточной точностью можно принять рав­ ным 0,1 количества C4HS + C4HG

Гидродинамические расчеты (сопротивления слоя, высоты затво­ ров, распределительных устройств для катализатора) проводятся известными методами [349, 355].

П р о в е р к а р а с ч е т о в . При оптимальном режиме вьь ход по расчету 40 об. %, конверсия 50 об. %; выход на полупромыш­ ленной установке при тех же условиях 40—41 об.%, конверсия — 48—52 об.% [354].

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

ОСНОВЫ РАСЧЕТА РЕАКТОРА СО ВЗВЕШЕННЫМ СЛОЕМ КАТАЛИЗАТОРА [191, 354]

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

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

 

(IX,9)

Как было показано ранее (см. главу IV),

 

 

(IX, 10)

где ш,- — линейная скорость газа в реакторе (на

общее сечение),

шкр — скорость начала образования взвешенного

слоя (критиче­

ская).

 

Скорости крекинга и углеотложения в первом приближении не зависят от проскока газа; скорость углеотложения не зависит от объемной скорости, поэтому проскок бутана не влияет на зауглероженность катализатора. В связи с этим соотношение основной и побочных реакций при дегидрировании бутана во взвешенном слое катализатора. несколько отличается от соотношения скоростей в

15*

227


реакторе с плотным слоем катализатора. Это сказывается на величи­

не коэффициентов а

и ß; при дегидрировании бутана во взвешен­

ном слое а = 1 +

1,15 X и

ß =

1 + X.

У р а в н е н и я

т е п л о

в ы х

б а л а н с о в , число которых

равно количеству слоев, связывают температуру катализатора в слое с кратностью циркуляции катализатора, температурой ката­

 

 

Выход

лизатора

и газа на входе

в реактор и на

 

В х о д

выходе, глубиной конверсии бутана в слое

 

частиц

г а з а

 

 

U L

и эффективностью

секционирующих реше­

 

 

ток по перемешиванию твердой

фазы.

 

 

 

Для

упрощения

расчетов

принимаем,

 

\

 

что перемешивание катализатора в каждом

 

/

слое идеальное, т. е. температура катали­

 

 

 

затора в пределах одного слоя постоянная;

 

 

 

вследствие хорошей теплопередачи во взве­

 

 

 

шенном слое [215] температура газа на вы­

 

 

 

ходе из соответствующего

слоя равна тем­

 

 

 

пературе

этого слоя.

 

 

 

 

 

 

щшШ

Для п-слойного реактора, разделенного

 

решетками провального

типа

с

эффектив­

 

-ШЁШ-

ностью ц, в расчете на

подачу бутана F =

 

 

4-

1

 

= 1 кг-моль/ч = 58

кг/ч

можно

написать

1ы,%4

следующую

систему

уравнений

тепловых

I

балансов с учетом влияния перемешивания

 

-

катализатора

между

слоями

(нумерация

1

 

 

,,

v>,;-g5

слоев по ходу газа):

 

 

 

 

 

 

tun-VU R B

580 ск (Т2-

7\) =

сг (7\ -

g

+

58<7 • AXlt

 

Л

 

 

 

 

 

 

 

 

В х о д г а з а __V

 

ск (Т3

Т2) = сГ (Т2

 

7\) +

 

 

 

 

 

f В ы х о д ч а с т и и ,

Рис. 51. Схема реактора со взвешенным слоем катали­ затора, секционированного провальными решетками:

hp — расстоян и е меж ду

ре ­

ш еткам и, Лр з — вы сота

р а з ­

реж ен ной зоны .

 

+ 58Сск(Т2- Г 1) ( ^1-

1) +

58? • АХ2,

58О

 

 

 

С К ( Т п

Т п — \ ) -

С Г ( Т „ „ I —

Т п — 2 ) +

+ 58GcK„_і - Г

п_2) ^ -----1) +

+

58^ ■АХп—і

(IX, 1 1 )

 

58GcK(Тн-

Тп) =

сг (Та -

Тп-0 +

 

 

 

+ 58GcK(Т„ - -7Ѵ_,) Ң -----і) +

58<7 АХп.

 

Кроме того,

АХх +

АХй+

— +

АХп = Хк.

 

 

 

 

В

этих уравнениях

G — кратность циркуляции катализатора,

кгікг

бутана;

Т ъ

Т%,

Тп — температура

катализатора

в со­

ответствующем

слое,

°С; ср,

ск — средние

теплоемкости

газа,

228


ккалімоль ■град, и катализатора, ккал/кг ■град; А Х 1г АХ2, ..., Х к

— приращение конверсии

по слоям и общая

конверсия бутана;

q — тепловой эффект реакции, ккал/кг бутана;

q зависит от изби­

рательности и изменяется в пределах 380—460 ккал/кг.

У р а в н е н и я г и д р о д и н а м и ч е с к и х х а р а к т е ­

р и с т и к р е а к т о р а .

В вышеприведенных уравнениях в яв­

ном виде фигурируют две величины: критическая скорость газа wKр и эффективность решеток по перемешиванию твердой фазы т].

Для катализатора К-5 [354] критическая скорость (м/сек)

шкр= 1,485ІО“ 4

d1'82

(IX, 12)

г"

а эффективность секционирующих решеток провального типа

Ч =

 

(IX, 13)

и р

 

 

где dcр — средний размер частиц катализатора, мм;

р — вязкость

контактного газа, спз; Gu — скорость

циркуляции

катализатора

через реактор, кг/сек; Gp — количество катализатора (кг/сек), кото­ рое решетка способна пропустить при данной скорости газа в ее отверстиях.

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

Р а с ч е т р е ш е т о к [264]. Литературные данные по иссле­ дованию ситчатых провальных решеток для секционирования реак­ торов со взвешенным слоем катализатора почти отсутствуют.

Свободное истечение мелкозернистого материала через одиночное отверстие изучалось в работах [357, 358]. Показано [357], что в этом случае скорость истечения q выражается уравнением

<7 = ц (х>У 2gr,

где ц — коэффициент, со — площадь отверстия, т— радиус отвер­ стия.

Условия, при которых частицы взвешенного слоя начинают про­ валиваться через отверстия поддерживающих (распределительных) решеток, исследованы в работах [358, 359]. По данным работы [358], скорость газа в отверстиях, при которой начинается проваливание частиц через распределительную решетку, зависит от конструкции решетки и определяется только свойствами газа и твердых частиц.

229


В работе [264] исследование ситчатых решеток проводилось на стеклянных моделях диаметром 56 и 105 мм. Характеристика приме­ нявшегося мелкозернистого материала приведена в табл. 76. Взве­ шенный слой создавался воздухом, углекислым газом или водоро­ дом.

 

Т а б л и ц а

76

 

Характеристика мелкозернистого материала

 

 

 

 

П лотность

Разм ер

 

М атериал

слоя,

 

частиц,

мкм

 

 

кг/л’

 

 

 

Катализатор К-5

1125

115-150

 

Алюмосиликатный катализатор

960

250—450

 

950

110—400

 

Песок

1750

250—450

 

Диатомовый кирпич

420

250-450

 

 

 

450—600

 

Использовались ситчатые решетки с диаметром

отверстий: 2;

3; 5; 5,2; 7,1; 10; 14,5; 18,6 мм и

площадью отверстий около 16%

площади свободного сечения модели.

 

 

 

Опыты были разделены на три серии: исследование пропускной

способности решеток при свободном истечении

частиц без

подачи

в модель газа; определение предельной скорости газа wnp

отвер­

стиях решеток), при которой просыпание частиц через отверстия прекращается; определение пропускной способности решеток во всем интервале скоростей подачи газа от w = 0 до w = wnр.

На основании общих соображений и анализа эксперименталь­ ных данных [357, 358] можно предположить, что скорость свободно­ го истечения мелкозернистого материала GMaKC— функция плот­ ности слоя ртв, диаметра отверстия решетки d0, числа решеток п и ускорения силы тяжести g, т. е.

Смаке = Ф (Ртв> ^оі ё)- (IX, 14)

В уравнении (IX,14) четыре размерные величины (я = 4) и три основные размерности = 3). По теории подобия, число безразмер­ ных групп, определяющих данное явление, равно п р — 4 —3 = = 1. В соответствии со сказанным, уравнение (IX,14) можно пере­ писать следующим образом:•

- т - = ф (")-

(IX, 15)

 

Ртв^об

 

в котором левая часть уравнения и есть искомая безразмерная груп­ па. Методом анализа размерностей найдено, что а — 1, b = 1, с = 2,5, d = 0,5; откуда

РтА V ëdo = Ф(л).

(IX, 16)

 

230