книги / Теория химических реакторов введение в основной курс
..pdfРис. 59. Зависимости расхода, температуры и концентраций от времени для реактора периодического действия с регулятором
171
Начальные условия при решении системы уравнений следующие:
172
12. ЧИСЛЕННАЯ МОДЕЛЬ РЕГЕНЕРАТОРА КАТАЛИЗАТОРА УСТАНОВКИ КАТАЛИТИЧЕСКОГО КРЕКИНГА
При работе установки каталитического крекинга газойлевых фракций на движущемся пористом гранулированном катализаторе последний загрязняется отложениями углерода, теряет активность и селективность. В связи с этим необходимо непрерывно выводить его из процесса и выжигать накопившийся углерод. При этом процесс крекинга, проходящий с выделением углерода, и процесс выжигания отложившегося углерода, должны быть согласованы по тепловому балансу, т.е. потеря температуры катализатора при протекании эндотермических реакций каталитического крекинга должна полностью компенсироваться тепловым эффектом выжигания коксовыхотложенийсучетомтепловыхпотерь(рис. 60).
Техническая реализация совмещенного процесса предусматривает, как правило, расположение реактора и регенератора либо друг над другом, либо рядом, с промежуточным подъемом катализатора элеватором(устаревшийвариант), либо газлифтнымаппаратом.
Рис. 60. Движение потоков в регенераторе
173
В регенераторе возможны два варианта выжигания кокса: путем подачи окислителя противотоком к движению катализатора либо путем подачи воздуха (или горячие дымовые газы с остатками кислорода) в регенератор в виде сопутствующего тока.
Поскольку основная часть катализатора двигается сплошным потоком и не имеет контакта со стенками аппарата для передачи тепла в окружающую среду, любой выделенный микрообъем катализатора этого потока находится в условиях, близких к адиабатным. Практически вся теплота, выделяющаяся в этом микрообъеме, в нем и остается, вызывая его разогрев, и только незначительная часть теплоты все же переносится движущимся газообразным окислителем, имеющим многократно меньшую теплоемкость, чем твердый катализатор.
Рассмотрим математическую модель противоточного регенератора.
12.1. Распределение воздуха по сечению регенератора
Высота слоя катализатора по сечению реактора-регенера- тора не является постоянной величиной в связи с конструкционными особенностями аппарата. Разобьем условно площадь поперечного сечения регенератора на квадратные площадки со стороной a (м).
Введем также следующие обозначения: i – номер квадратной площадки;
hi – высота столбика катализатора на i-й площадке, м;
B – общий суммарный расход воздуха через все сечение зоны регенератора, м3/с;
Vi – скорость движения воздуха через i-ю площадку, м/с;
ξ – коэффициент сопротивления слоя катализатора единичной длины;
F – общая площадь поперечного сечения катализатора, м2.
174
Учитывая наличие достаточно обширных буферных емкостей выше и ниже регенератора для выравнивания давления воздуха и дымовых газов, которое происходит практически мгновенно (со скоростью звука в газе), можно считать, что перепады давления ∆pi по любому i-му сечению являются одинаковыми. В этом случае
в первом приближении |
расход |
воздуха (газов) через |
любую |
|||||
i-ю площадку может быть вычислен следующим образом: |
|
|||||||
|
p h V 2 |
, |
|
|
(224) |
|||
|
|
i |
|
i |
|
|
|
|
отсюда |
|
|
|
|
|
|
|
|
|
V |
p |
. |
|
|
|
(225) |
|
|
|
|
|
|
||||
|
i |
hi |
|
|
|
|
||
|
|
|
|
|
|
|||
Следовательно, расход воздуха через i-ю ячейку |
|
|||||||
|
Bi = a2·Vi. |
|
|
|
(226) |
|||
Сумма всех расходов должна быть нормирована на общий |
||||||||
расход воздуха по всему сечению зоны, т.е. |
|
|||||||
Bi a2 |
Vi a2 |
|
|
p |
B. |
(227) |
||
|
|
|||||||
i |
i |
|
|
i |
|
hi |
|
Исходя из этого легко может быть найдено выражение для расхода воздуха через i-ю площадку регенератора в зависимости от фактической формы поверхности катализатора и наличия местных сопротивлений, или размера зерен катализатора в местном сечении реакционного объема.
12.2. Тепловой баланс элементарного объема
Разделим i-й столбец катализатора на n ячеек с номером k. Тогда для каждого элементарного объема катализатора в системе координат, связанных с движущимся катализатором, можно написать уравнение баланса тепла (Дж) без учета теплообмена между соседними ячейками (пренебрегаем):
175
qкат |
qвозд q |
рекц |
qкат qвозд ; |
(228) |
||||
нач |
вх |
|
кон |
вых |
|
|||
|
qкат |
Скат tкат |
mкат ; |
|
(229а) |
|||
|
нач |
Р |
|
нач |
|
|
|
|
qвозд Свозд tвозд |
mвозд; |
|
(229б) |
|||||
вх |
Р |
вх |
|
|
|
|
||
|
qкат |
Скат tкат |
mкат; |
|
(229в) |
|||
|
кон |
Р |
|
кон |
|
|
|
|
qвозд |
Свозд |
tвозд mвозд; |
(229г) |
|||||
|
вых |
Р |
|
вых |
|
|
|
|
qрекц |
С |
|
|
Qr |
, |
(230) |
||
100 mкат |
где СРвозд, СРкат – теплоемкость при постоянном давлении воздуха и
катализатора соответственно, Дж/(кг·К); t – температура воздуха и катализатора (°С) в начале и конце элементарного интервала време-
ни , с; mвозд, mкат – масса элементарного объема воздуха и
катализатора, кг; Qr – удельная теплота сгорания углерода в зависимости соотношения СО/СО2; при соотношении СО/СО2 = 1,
Qr = 22 190 кДж/кг; С – изменение концентрации углерода в выделенной ячейке за элементарный промежуток времени, мас. %;O2 – концентрациякислородавпотокегаза.
Кинетическое уравнение горения углерода в катализаторе имеет вид
W f C , O2 , t . |
(231) |
В рамках некоторых ограничений реальный вид этого уравнения может быть представлен экспоненциальной формой:
d C |
A e (t 600) , |
(232) |
|
d |
|||
|
|
где коэффициент может быть найден экспериментально или определен непосредственно по результатам обработки данных для реальной промышленной установки.
176
Из литературных данных известно, что порядок реакции горения кокса зависит от его содержания в катализаторе. При [C] > 2 % порядокнулевой, применьшихконцентрациях– первый.
На основании промышленных данных для регенератора установки КК-1 ООО «ЛУКОЙЛ-Пермнефтеоргсинтез» были получены следующие результаты:
А = 4,245·10–4 [O2] при [C] > 2 %; |
|
|
|
|
||||||||||||||||||
A |
4,245 10 4 O2 C 0,5 , |
|
|
|
|
|
(233) |
|||||||||||||||
|
|
|
2 0,5 |
C 2 %, |
|
|
|
|
|
|
|
|||||||||||
|
|
|
при 0,5 % |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
А = 0 при C 0,5 % . |
|
|
|
|
|
|
|
|
|||||||||||
Таким образом, для конечных приращений выражение (232) |
||||||||||||||||||||||
может быть записано в следующем виде: |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
нач |
|
|
|
|
|
|
|
|
|
|
|
||||
[C] A exp |
0,00768 tкат |
|
|
600 |
|
. |
|
|
|
(234) |
||||||||||||
|
|
|
|
|
|
|||||||||||||||||
Изменение концентрации кислорода |
[O2] определяется в |
|||||||||||||||||||||
зависимости от изменения концентрации углерода С |
с учетом |
|||||||||||||||||||||
стехиометрии реакции горения: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
C 1/2O2 CO 12 p1 |
|
16 p1 |
22 p1 |
|
mвозд О2 |
|
|
16 p |
32 p |
|||||||||||||
C O2 CO2 12 p2 |
32 p2 |
|
38 p2 |
|
|
|
|
|
1 |
|
2 , |
|||||||||||
|
|
mкат С |
|
|
12 |
|||||||||||||||||
O |
|
C 3p 8 p |
|
|
|
|
|
mкат |
|
, |
|
|
|
|
(235) |
|||||||
|
|
3 |
mвозд |
|
|
|
|
|||||||||||||||
|
2 |
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
где р1 и р2 – доли углерода, сгорающего до СО и СО2 соответственно.
В математической модели сделаем следующее допущение:
tвыхвозд tконкат , |
(236) |
обозначим эти температуры Тк.
177
Масса катализатора в элементарной ячейке |
mкат определя- |
|
ется выражением |
|
|
mкат hi |
a2 кат , |
(237) |
n |
|
|
где ρкат – плотность катализатора.
Элементарный промежуток времени, за который указанная масса катализатора проходит через элементарную площадку a2:
|
mкат 3,6 F |
, |
(238) |
|
К a2 |
||||
|
|
|
где F – площадь поперечного сечения регенератора, м2; К – циркуляция катализатора, т/ч.
Завремя черезсечениеа2 проходитобъемвоздухамассой
mвозд V a2 |
возд. |
(239) |
i |
|
|
Обозначим шаги интегрирования по времени как j и введем признак движения катализатора shift, который будет принимать значение 0 и 1 для неподвижного и движущегося катализатора соответственно.
Тогда переменными в уравнениях являются Tkj shift1 tначкат и Tkj 1 tвхвозд (k – номер слоя, для самого нижнего слоя k = 1, для самого верхнего k = n). Температуру катализатора, поступающего в регенератор сверху, обозначим Tnj 1 , а температуру воздуха, по-
ступающего снизу, – T0j .
Концентрация кислородав поступающем воздухе определяется внешними условиями, связанными с возможным незначительным подогревом его за счет частичного выжига кислорода в топке под давлением. Она может быть приближенно описана в зависимости от температурыпоступающего воздухапоуравнению
О2 21 21/1500 Т0j %. |
(240) |
178
Всоответствиисэтимможнопреобразоватьвыражение(234):
C 6,67 10 5 O2 k 1 exp 0,00768 Tkj shift1 600 . (241)
Для кислорода |
|
|
|
|
|
|
|
|
|
|
|
|
где |
∆[O2] = S·∆ [C], |
|
|
|
|
(242) |
||||||
8 p 16 p |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
h |
|
кат |
|
|
1 |
|
|
||
S |
1 |
|
|
1 |
|
|
|
|
|
, |
(243) |
|
3 |
|
|
возд |
V |
|
|||||||
|
|
|
|
n |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
i |
|
|
|
окончательно, после соответствующих подстановок, можно получить соотношение
T |
j |
|
A1 Tkj shift1 |
A2 Tkj 1 A3 |
C |
, |
(244) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
k |
|
|
|
|
|
A1 |
A2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A Cкат |
|
|
|
hi , |
|
|
|
|||||
|
|
|
1 |
Р |
|
|
кат |
|
n |
|
|
|
|||
|
|
|
A |
Cвозд V |
|
возд , |
|
|
(245) |
||||||
|
|
|
2 |
Р |
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
A3 100 |
|
|
A1 |
|
|
|
|
|
||||
|
|
|
Cкат Q |
|
|
|
|
||||||||
|
|
|
|
|
|
|
Р |
|
|
|
r |
|
|
|
|
Энергия, накапливающаяся в катализаторе в виде тепла, |
|
||||||||||||||
|
|
|
|
q j A a2 |
T j . |
|
|
(246) |
|||||||
|
|
|
|
k |
1 |
|
|
|
|
k |
|
|
|
В начале расчета задается начальный массив концентраций
углерода по высоте i-го столбца: [C]0k = [C]нач = const, и массив начальных температур: T0k = Tнач = const.
Изменение концентрации кислорода и углерода определяется пошагово по выражениям:
C j |
C j 1 |
C , |
(247а) |
|
|
k |
k shift |
|
|
O2 |
|
O2 |
O2 . |
(247б) |
|
k |
k 1 |
|
|
|
|
|
|
179 |
Для установки КК-1 учтено, что регенератор состоит из двух одинаковых секций высотой по 2,5 м, в которые подается одинаковое количество воздуха. Подача свежего воздуха осуществляется под (n/2+1)-ю ячейку.
На рис. 61–63 показаны переходные процессы в модели при реализации различных возмущений на входе в регенератор. На рисунках представлено изменение следующих параметров:
–температуры термонапряженного слоя – где была достигнута максимальная температура;
–температуры на термопаре I (II). Расположение термопар задано и отмечено на рисунках пунктиром (правая шкала);
–температуры нижнего слоя – выходящего из регенератора
снизу;
–коэффициента поглощения кислорода – выражается отношением концентрации кислорода в воздухе, поступающем в верхнюю секцию регенератора, к концентрации кислорода в газе, выходящего из секции, деленным на 10 для уменьшения масштаба;
–максимальной температуры – достигнутая в регенераторе в процессе работы;
–положения термонапряженного слоя – положение слоя, считаяснизу, вкоторомбыладостигнутамаксимальнаятемпература.
Модель максимально подробно описывает работу регенератора установки каталитического крекинга и может продемонстрировать поведение любого учтенного параметра как по высоте слоя катализатора, так и во времени.
Наиболее интересным эффектом, выявленным при математическом моделировании, является установление факта движущихся «высокотемпературных зон» («вспышек») в катализаторе. Установлено, что они появляются внезапно с резким возрастанием температуры, иногда до недопустимых величин. Потом они постепенно сдвигаются вместе с катализатором и исчезают (см. рис. 63, кривые для положения термонапряженного слоя).
180