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

книги / Теория химических реакторов введение в основной курс

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
5.34 Mб
Скачать

Рис. 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 = aVi.

 

 

 

(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

Соседние файлы в папке книги