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

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

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

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

Добавлен: 11.04.2024

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

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

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

Математическое описание процесса. Процесс дегидрирования

бутана можно описать

тремя группами уравнений.

С т е х и о м е т р и

ч е с к и е у р а в н е н и я . При дегидри­

ровании бутана на промышленном катализаторе в условиях, близ­ ких к производственным, процесс описывается следующими уравне­ ниями реакций:

н-С4Н10

я-C.jH8 -f Н2,

 

 

 

(а)

я-С4Н8 -ч=±

С4Н0 +

Н2,

 

 

 

(б)

н-С4Н10

ызо-С4Н10,

 

 

 

(в)

я-С4Н8

изо-С4Н8,

 

 

 

(г)

я-С4Н10

0,62СН4+ 0,13С2Н4 +

0,44С2Н„ +

 

 

+ 0,04С3Н„ + 0,2С3Н8 +

0,ЗС5,

 

 

(д)

С4Н10 ->■ч-С4Н8 ->

С4Не

С8Н12 (4-винилциклогексан-1)

 

поликонденсация, дегидрирование, крекинг -*■

 

полициклические углеводороды (смолы или «уголь»)

 

 

пли упрощенно:

^

_ 4С +

 

 

(е)

У р а в н е н и я

с к о р о с т е й

р е а к ц и й . В

эту группу

входят уравнения скоростей реакций дегидрирования

бутана, кре­

кинга бутана

и бутилена, а также углеобразования

(см. гл. IV).

Дегидрирование бутана.

 

 

 

 

 

ЛХ

_ и рс4н10

P Ct H„P H ,

\

(IX, 1)

 

d(W/F)

к

р0.5

ХрРС4Н,0

/

 

 

 

 

 

С>4г1в

 

Константа скорости реакции находится из уравнения *

 

 

\gk =

— 41 400

+ 9,0.

 

(IX,2)

 

 

 

4,5757

 

 

 

Константа равновесия Кр рассчитывается по уравнению [354]

 

lg^p =

— 30 500

+ 7,574.

 

(IX,3)

 

4.575Т

 

Крекинг бутана и бутилена. Приближенное уравнение имеет

вид

 

 

dxk

 

 

 

 

 

= К (Рс4нІ0 + p ctH,).

(IX,4)

 

 

 

d(W/F)

Уравнение (IX,4) в интегральной форме при

Р0бщ = 1 атм

имеет следующий вид:

(W/F)0,5-----ln (W /F f,b0,5+

 

In 1

1

хк

2kK

1

* Уравнение

(IX ,1) справедливо до температур 550° С; выше 550° С

 

 

 

'

+ 2'2'

 

221


Для промышленного катализатора К-5

 

lgv4 =

— 2300

+

1,727,

 

 

т

 

 

 

lg*K =

— 14 300

 

1,035.

(IX,5)

4,5757

 

Количество отлагающегося на катализаторе «угля» (% к весу

катализатора) [168]:

 

С = Л / - 68,

(IX,6)

где т — время, мин\ А х — температурный коэффициент, определяе­ мый из уравнения

lg Л = 20 + 6,57. (IX,7)

Выход бутадиена примерно в десять раз меньше выхода бутиле­ на; скорости реакций (в) и (г) незначительны, приблизительно 0,05 скорости реакции (а), и при расчетах состава продуктов реакции могут не учитываться.

Значение констант скорости дегидрирования, рассчитываемое по уравнению (IX,2), справедливо для размера частиц катализатора 0,7 мм и меньше; в случае использования частиц большего размера (например, шарикового катализатора) в найденную по этому урав­ нению константу необходимо ввести поправку, определяемую из соотношения [354]:

где kd — константа скорости при размере частиц катализатора d (мм).

Этой поправкой учитывается степень использования внутренней поверхности шарикового катализатора. Скорость остальных реак­ ций не зависит от размера частиц катализатора.

Уравнение изменения активности катализатора за счет длитель­ ной эксплуатации в систему уравнений на включено, так как сред­ няя активность поддерживается на постоянном уровне за счет перио­ дического ввода в систему некоторого количества свежего катализа­ тора. Изменение активности катализатора при прохождении из реак­ тора вследствие удаления избыточного кислорода может быть учтено на основании данных, приведенных в главе V.

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

скоростью

циркуляции

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

Примем

следующие

обозначения:

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

катализатора соответственно на входе в ре­

актор и на

выходе из него, °С; tH— температура бутана на входе

в реактор,

°С; tK— температура контактного газа на выходе из

222


реактора, °С; Q — скорость циркуляции

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

кг/кг бу­

тана; q — суммарный тепловой эффект

дегидрирования

бутана,

ккал/кг; ср — соответствующие удельные теплоемкости, ккал/кг х

X град — для катализатора или ккалікг-моль

■град — для

газа;

Хк — конечная степень превращения бутана,

мольные доли;

ин­

дексы при ср: кт — катализатор; к. г — контактный газ.

 

Тепловой баланс 1 кг-моль бутана (58 кг) при конверсии Хк за­ пишется так:

58Q (Сркт • Т п Срк Т к) = EcpiirtK cpCiiitJu + 58^ХК.

Найдено [354], что с погрешностью меньше 5% мольная тепло­ емкость бутана равна сумме мольных теплоемкостей продуктов реак­ ции, т. е.

СРС,Ніа ~ 2 с Рк.г=

ср -

Тогда

 

58Q (сРктТи сРктТк) = Ср (/к

tH) -j- 58qXK,

откуда

 

т -

1к

Ркт

СР ^н)

хк.

(IX,8)

 

58Qc

Ркт

Qc

 

Ркт

Ркт

 

Разность температур в реакторе между катализатором и газом не превышает 5°, поэтому расчет реактора будем вести по температу­ ре катализатора, которую для любой конверсии X можно найти по видоизмененному уравнению (IX,8):

т =

Т„-

 

-----(Хк - X).

(IX,8а)

Ч

т

58QCPKT

QCP KT

 

В этом уравнении при X =

О Т =

Тк, при X = Хк Т =

Тк, Тн >

>>ік\ для упрощения можно принять t = tK.

Расчет реактора. Покажем, как на основе этого математического

описания можно рассчитать реактор.

И с х о д н ы е д а н н ы е . Из помещенных в гл. IV данных сле­ дует, что при дегидрировании бутана на катализаторе К-5 в лабора­ торных условиях максимальная конверсия достигает 75% равновес­ ной. Примем, что при оптимальных условиях избирательность де­ гидрирования должна быть порядка 80—85 мол.%.

Проверочными расчетами (после определения объема реактора) по уравнениям (IX,4) и (IX,5) можно показать, что температура контактного газа на выходе из реактора в этом случае не должна превышать 600° С, а следовательно, температура катализатора на входе в реактор должна быть не больше 610° С. При 600° С равновес­ ная глубина дегидрирования составляет около 67% и, следователь­ но, конверсия может достигать 50%. Максимальная температура перегрева исходного бутана (перед поступлением в реактор) опреде-

2 2 Д


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

О б ъ е м р е а к т о р а . Количество катализатора (величина W/F), требуемое для достижения заданной конверсии, определяется уравнениями (IX,1) и (IX,2), в которых температура (величина пере­ менная) находится из уравнения (IX,8а).

Покажем, как найти объем реактора простейшим графическим способом. При принятых оптимальных условиях определение объе­ ма реактора сводится к совместному решению уравнений (IX, 1), (IX,2) и (IX,8а) относительно величины W/F. Это можно выполнить следующим образом. Выразим парциальные давления реагентов в уравнении (IX,I) через величину конверсии бутана за проход X. Пусть общее давление в реакторе Р, атм. При конверсии X и проте­ кании только реакции дегидрирования в контактном газе будет со­ держаться следующее количество молей:

С4Н10 1 - Х

С4Н8 X

Н2 X Всего 1 + X.

Парциальные давления в данном случае определяются уравнени­ ями:

^с,н10 = -[ + х р-> р с,н, = Рнt = Tq“ x~ •

При протекании побочных реакций общее увеличение объема газа будет равно а, бутилена получается Ху, а водорода — Xyß. Напомним, что у — избирательность процесса, а ß — коэффициент избытка водорода, равный P HJ P C<H,- Тогда парциальные давления

р с,н,а

■—X Рс0\ Рс.„, = ^ - Р с » ; Рн, = ^ * -

Рсп

где с„ — содержание С4Н 10 в исходном бутане, подаваемом на де­ гидрирование.

Подставляя эти выражения в уравнение (IX, 1), получим

dX

_ и

(с0Р)°-5

1 - X I ,

с0у2ßP

X*

(IX, Іа)

d(W/F)

к

(вт)о.5 '

х о,5

Кра '

І - Х

 

Если обозначить через cp (X, t) величину, обратную выражению

(срР)0-5

1 - Х

л

с„у«рр

ха

\

 

(ау)0-5 '

X0-5

V

КР*

l - X

j ’

 

dX

k

 

W/F =■ J

l it ' Ü- . dX.

(IX, 16)

d{WIF) -

cp(X, t)

 

 

 

 

о

224


Необходимые для расчетов значения коэффициентов а, ß и у определяются из'уравнений, составленных по данным, приведенным

в главе

IV.

 

 

Для

реактора идеального вытеснения

а = 1 + 1 , 1 Х ,

ß =

1 Ч- 0,45ЛГ,

у = 0,5595 + 0,074Х + 0,00216/ — 0,8Х2 — 2,4 • ІО-6/2,

где t — температура,

°С.

 

Выражения для сс

и ß

справедливы и для дифференциального

реактора, поскольку они определяются только стехиометрией сум­ марного разложения, не за­

висящей от глубины превра­

 

 

 

 

 

 

 

щения. Избирательность у,

 

 

 

 

 

 

 

определяемая по приведен­

 

 

 

 

 

 

 

ному

выше

уравнению —

 

 

 

 

 

 

 

величина

интегральная.

 

 

 

 

 

 

 

Использование ее в уравне­

 

 

 

 

 

 

 

нии (IX,Іа)

вносит некото­

 

 

 

 

 

 

 

рую погрешность в опреде­

 

 

 

 

 

 

 

ление

W/F,

однако это

 

 

 

 

 

 

 

оправдывается тем, что при­

 

 

 

 

 

 

 

нимаемая в начале расчетов

о

4 8 0

5 0 0 5 2 0 5 4 0

5 6 0

5 8 0 6 0 0

общая величина у проверя­

 

 

 

Теппврашура,°С

 

 

ется затем при расчетах глу­

Рис. 49.

Изменение температуры

по высоте

бины крекинга

по уравне­

реактора (по ходу газа) при разных глубинах

ниям (IX,4) и (IX,5).

превращения

X

и скоростях

циркуляции

Q:

Графическое

решение

/ — <3 =

6; 2

7; 3 — 8: 4 — 9;

5 —

10; 5 —

12.

уравнений (IX ,16) и (IX,8а)

 

Прежде всего следует постро­

выполняется в следующем порядке.

ить график зависимости изменения температуры по высоте реактора

(по ходу газа) при

разных

глубинах превращения бутана X и ско­

ростях циркуляции катализатора Q по уравнению (IX,8а). Полу­

ченные адиабаты

представлены на рис. 49 (при Т я = 610° С,

іи —

550° С, tK=

600° С, <7 = 460 ккал/кг, Х к =

0,5).

 

Далее строим

изотермы

по уравнению

(IX, 16), т. е. получаем

зависимость cp (X,

t)/k от X

(рис. 50, пунктирные линии). Затем при

постоянном

Q, задаваясь

последовательно значениями X

(0,05;

0,1 и т. д. до 0,50), по рис. 49 находим величины Т, а по получен­ ным точкам проводим кривые (см. рис. 50, сплошные линии). Пло­ щадь, ограниченная кривой для заданного значения Q, осью абс­ цисс и ординатой Х к, равна W/F. При этих построениях константа скорости берется с учетом среднего размера частиц катализатора. Кроме.того, в найденную величину объема катализатора нужно вне­ сти поправку на изменение активности катализатора при прохожде­ нии им реактора (т. е. поправку на пониженную активность из-за избыточного кислорода)’ Соответствующий поправочный коэффи­ циент определяется на основании данных главы V.

1 5 3 — 1 3 1 8

225