Риформинг

Принятая для моделирования кинетическая схема риформинга

Для принятой кинетической схемы риформинга (рис. 8.1) уравне­ние (8.2а) трансформируется в систему уравнений, левые части ко­торых описывают накопление или расходование каждого класса углеводородов по координате реактора, а вид правых частей опреде­ляется реакциями образования и расходования данного класса углеводородов согласно кинетической схеме. При этом расходом углеводородов на образование кокса пренебрегают, а водород рас­сматривается как отдельный класс углеводородов из одного предста­вителя.

 

Принятая для моделирования кинетическая схема риформинга

 

Рис. 8.1. Принятая для моделирования кинетическая схема риформинга.

 

При проведении конкретных численных расчетов надо раскрыть структуры уравнений системы (8.2) и иметь значения констант.

В принципе кинетика реакций процесса риформинга должна опи­сываться уравнениями лэнгмюровского типа, однако поскольку ад­сорбционные коэффициенты различных углеводородов, участвую­щих в риформинге отдельных нефтяных фракций близки, знаменатель в указанных уравнениях существенно не меняется по ходу реакций и скорости их достаточно хорошо описываются уравнениями первого порядка по реагирующему углеводороду с эффективными констан­тами. Парциальное давление водорода по-разному влияет на отдель­ные реакции, составляющие процесс риформинга. Так, его увеличе­ние ускоряет реакцию гидрокрекинга и тормозит реакцию дегидроциклизации. Точно так же водород за счет гидрирования жидких полимеров уменьшает скорость коксоотложения. Принимается, что влияние парциального давления водорода достаточно удовлетвори­тельно описывается степенной функцией. По данным [316] показа­тель в ней равен для реакции дегидроциклизации 0,6, гидрокрекинга 0,6, изомеризации н-парафинов 0.   

Ход превращения смеси углеводородов можно представить как марковскую цепь. Это означает, что состояние системы в данный момент времени или, что то же, в данной точке по координате реак­тора полностью определяет вероятность ее состояния в следующий момент времени или следующей точке реактора. В соответствии со стохастической трактовкой кинетики химических реакций соотно­шение вероятностей взаимных превращений углеводородов идентично соотношению констант скорости этих превращений при равной струк­туре описывающих их уравнений. Тогда для решения поставленной задачи значения переходных вероятностей (вероятностей взаимных превращений углеводородов) должны описываться функцией от мо­лекулярной массы в каждом классе углеводородов, а для машинных расчетов такая зависимость может быть дана в виде таблицы.

Таблица 8.1. Таблица переходных вероятностей разрыва С—С-связей в н-парафинах (расчетные данные)

X

X

l

2

3

4

5

6

7

8

9

10

11

12

1

0

1

0,5

0,1

0,1

0,05

0,01

0,001

0

0

0

0

2

0

0

0,5

0,8

0,4

0,2

0,1

0,049

0,004

0,003

0,002

0,001

3

0

0

0

0,1

0,4

0,5

0,39

0,2

0,126

0,1

0,07

0,059

4

0

0

0

0

0,1

0,2

0,39

0,5

0,37

0,23

0,15

0,12

5

0

0

0

0

0

0,05

0,1

0,2

0,37

0,334

0,278

0,21

6

0

0

0

0

0

0

0,01

0,49

0,126

0,23

0,278

0,23

7

0

0

0

0

0

0

0

0,001

0,004

0,1

0,15

0,21

8

0

0

0

0

0

0

0

0

0

0,003

0,07

0,12

9

0

0

0

0

0

0

0

0

0

0

0,002

0,059

10

0

0

0

0

0

0

0

0

0

0

0

0,001

11

0

0

0

0

0

0

0

0

0

0

0

0

12

0

0

0

0

0

0

0

0

0

0

0

0

 

Таблица 8.2. Таблица переходных вероятностей разрыва С—С-связей в н-парафинах по данным [316]

X

X

1

2

3

4

5

6

7

8

9

10

11

12

1

0

1

0,5

0,4

0,2

0,1

0,05

0,001

0,001

0,001

0,001

0,001

2

0

0

0,5

0,4

0,3

0,15

0,1

0,005

0,069

0,078

0,048

0,018

3

0

0

0

0,1

0,3

0,3

0,12

0,12

0,1

0,1

0,1

0,10

4

0

0

0

0

0,2

0,2

0,15

0,15

0,15

0,15

0,12

0,1

5

0

0

0

0

0

0,25

0,35

0,37

0,4

0,35

0,27

0,25

6

0

0

0

0

0

0

0,23

0,3

0,25

0,23

0,23

0,22

7

0

0

0

0

0

0

0

0,001

0,029

0,08

0,15

0,1

8

0

0

0

0

0

0

0

0

0,001

0,01

0,06

0,09

9

0

0

0

0

0

0

0

0

0

0,001

0,01

0,07

10

0

0

0

0

0

0

0

0

0

0

0,001

0,05

11

0

0

0

0

0

0

0

0

0

0

0

0,001

12

0

0

0

0

0

0

0

0

0

0

0

0

 

В табл. 8.1 и 8.2 представлены значения переходных вероятностей, вычисленные на основе проведенных экспериментов, общетеоретиче­ских соображений органической химии и данных работы [316]. Существенное преимущество изложенного подхода то, что при по­строении кинетических моделей для конкретных катализаторов не требуется определять константы скорости реакций всех углево­дородов, а только отдельных, реперных.

Для конкретных расчетов по кинетической модели с учетом де­зактивации катализатора должны быть раскрыты структуры уравне­ний (8.26) и (8.2в).