НЕФТЕХИМИЯ, 2021, том 61, № 3, с. 359-366
УДК 66.011+665.644.
ДИСКРЕТНАЯ КИНЕТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА M-DSO
ДЛЯ ГИДРООЧИСТКИ БЕНЗИНА КАТАЛИТИЧЕСКОГО КРЕКИНГА
© 2021 г. Ying Xiang1, Jiaying Shen1, Fusheng Ouyang1,*
1 International Joint Research Center of Green Energy Chemical Engineering,
East China University of Science and Technology, Shanghai, 200237 China
*E-mail: ouyfsh@ecust.edu.cn
Поступила в редакцию 03 сентября 2020 г.
После доработки 27 декабря 2020 г.
Принята к публикации 18 марта 2021 г.
Селективное гидрирование бензина каталитического флюид-крекинга (FCC) является эффективным
методом снижения содержания в нем серы и олефинов для получения чистого бензина. На основе харак-
теристик процессов DSO-M и M-DSO для гидроочистки тяжелой фракции бензина FCC и определения
вклада различных углеводородов в октановое число разработаны две пятикомпонентные реакционные
схемы для стадии обогащения M и стадии гидродесульфурации DSO, соответственно. По эксперимен-
тальным данным, полученным при использовании микротрубчатого реактора с неподвижным слоем
для гидроочистки бензина FCC, с помощью алгоритма Рунге-Кутта и генетического алгоритма были
рассчитаны кинетические параметры 5-компонентных моделей для стадий M и DSO, соответственно.
Кинетический анализ показывает, что процесс M-DSO лучше, чем процесс DSO-M. Более того, проверка
показывает, что установленные модели обладают хорошей надежностью и пригодны для экстраполяции.
Ключевые слова: бензин FCC, гидроочистка, процесс M-DSO, дискретная кинетическая модель
DOI: 10.31857/S0028242121030072
В связи с постоянным ужесточением требова-
ния олефинов и восстановлении октанового чис-
ний по охране окружающей среды почти во всех
ла [5]. Функция катализатора DSO заключается в
странах предъявляются более высокие требования
гидродесульфурации бензина FCC [6]. Промыш-
к содержанию серы и олефинов в бензине. Меж-
ленные результаты показали, что процесс M-DSO
ду тем, бензин FCC составляет более 30% объема
имел такое же влияние на десульфурацию и сни-
мирового производства бензина. Следовательно,
жение содержания олефинов в бензине FCC, при
очистка бензина FCC является ключом к повыше-
этом октановое число было увеличено примерно на
нию качества бензина. Был разработан целый ряд
0.8 единицы, по сравнению с процессом DSO-M.
процессов гидроочистки бензина FCC, в частно-
Существующие кинетические модели для ги-
сти, процессы RSDS [1], Prime-G+ [2], SCANfining
дроочистки бензина FCC [7] не учитывают в до-
[3], DSO-M и т. д. Среди них процесс DSO-M [4],
статочной степени все реакции, поэтому они не
разработанный компанией PetroChina, сыграл важ-
играют большой роли при моделировании и оп-
ную роль в десульфуризации и восстановлении
тимизации процессов гидроочистки бензина FCC.
олефинов. Однако после первого использования
Хотя гидроочистка бензина FCC была улучшена
этого процесса потеря октанового числа бензина
за счет преобразования процесса DSO-M в про-
составила около 1.7 единицы, что серьезно повлия-
цесс M-DSO, мало что известно об изменениях в
ло на экономический эффект.
протекании этих реакций и влиянии условий реак-
Процесс M-DSO развивается путем реверсиро-
ции на процесс. В этом исследовании, основанном
вания стадий DSO и M. Функции катализатора М
на всестороннем рассмотрении различных реак-
заключаются в ароматизации, снижении содержа-
ций бензина FCC и метода дискретной кинетики
359
360
YING XIANG и др.
Рис. 1. Дискретные схемы реакций стадий (а) M и (б) DSO. kj - константа скорости реакции на стадии M; kj - константа
скорости реакции на стадии DSO; A, O, N, nP, iP - ароматические соединения, олефины, нафтены, н-парафины и изопара-
фин, соответственно.
[8], ожидается прояснение того, почему процесс
активные центры всех реакций имеют одинаковую
M-DSO демонстрирует более высокую эффектив-
активность; используемый трубчатый реактор с
ность при очистке бензина по сравнению с процес-
неподвижным слоем можно приблизительно рас-
сом DSO-M.
сматривать как реактор идеального вытеснения,
поэтому диффузия внутри частиц не учитывается;
благодаря защите циркулирующего водорода деак-
РАЗРАБОТКА ДИСКРЕТНОЙ РЕАКЦИОННОЙ
тивацией катализатора можно пренебречь.
МОДЕЛИ И СОСТАВЛЕНИЕ КИНЕТИЧЕСКОГО
На основании сделанных выше предположений
УРАВНЕНИЯ
кинетическое уравнение реакции может быть вы-
Поскольку выход бензина в обоих процессах
ражено формулой (1) [9-11]:
DSO-M и M-DSO составляет около 99%, газы от
крекинга бензина в этом исследовании не учиты-
(1)
вались. Для углеводородов с одинаковым числом
атомов углерода порядок распределения их окта-
где K - матрица констант скорости реакции, состо-
нового числа от высокого к низкому следующий:
ящая из kj; a - вектор массовых долей всех компо-
ароматические углеводороды > изопарафины > на-
нентов (рис. 2, 3); X - безразмерное относительное
фтены ≈ олефины > н-парафины. В соответствии
расстояние в сечении X в слое пласта; SWH - объем-
с вкладом различных углеводородов в октановое
ная скорость (WHSV). Пояснения ко всем симво-
число реакции стадий DSO и M были разделены
лам приводятся в работе [12].
на пять групп, а именно: ароматические углеводо-
роды (A), олефины (O), нафтены (N), н-парафины
ЭКСПЕРИМЕНТЫ ПО
(nP) и изопарафины (iP). На основе механизма ги-
СОВЕРШЕНТСВОВАНИЮ ГИДРООЧИСТКИ
дроочистки и характеристик стадий DSO и M были
созданы две 5-компонентные схемы реакций для
Эксперименты по гидроочистке бензина FCC
стадий M и DSO, соответственно (как показано на
проводили с тяжелыми фракциями бензина FCC из
рис. 1).
ректификационной колонны (далее именуемыми
сырьем) в лабораторном устройстве, представля-
Кинетическое уравнение было составлено на ос-
нове следующих допущений: согласно существую-
ющем собой микротрубчатый реактор с неподвиж-
ным слоем.
щей литературе по дискретной кинетике, все реак-
ции обычно рассматриваются как реакции первого
Для процесса M-DSO сначала сырьем стадии M
порядка, что подходит и для данного исследования;
было исходное сырье (номер 0#), а катализатором -
НЕФТЕХИМИЯ том 61 № 3 2021
ДИСКРЕТНАЯ КИНЕТИЧЕСКАЯ МОДЕЛЬ
361
Рис. 2. Матрица констант скорости для кинетиче-
Рис. 3. Матрица констант скорости для кинетиче-
ской модели стадии M. [kN = -(k5+k6+k7+k10), kO =
ской модели стадии DSO. [kN = - (k5+k6+k7+k10), kO =
-(k1+k2+k8+k9), kA = -k11, kiP = -k4, knP = -k3).
-(k1+k2+k8+k9), kA = 0, kiP = -k4, knP = -k3).
M при давлении 2.0 МПа, объемной скорости
рьем стадии DSO был продукт 3#, эксперименты
1.0 ч-1 (WHSV), температурах 350, 370 и 390°C,
проводились при температуре 240, 260 и 280°C,
соответственно; полученные продукты были по-
соответственно; соответствующие продукты были
следовательно пронумерованы как 1#, 2#, 3#. Для
последовательно пронумерованы как 10#, 11#, 12#.
стадии DSO сырьем были продукты 1#, 2#, 3#, а ка-
Точно также для процесса DSO-M сначала сы-
тализатором - M при 2.3 МПа, WHSV 2.0 ч-1, соот-
рьем стадии DSO было исходное сырье (номер
ветственно. Когда сырьем стадии DSO был продукт
0#), а катализатором - DSO при 2.3 МПа, WHSV
1#, эксперименты проводили при температурах
2.0 ч-1 и 240, 260, 280°C, соответственно; эти про-
240, 260 и 280°C, соответственно; соответствую-
дукты были последовательно пронумерованы как
щие продукты реакции были последовательно про-
13#, 14#, 15#. Сырье на стадии M представляло
нумерованы как 4#, 5#, 6#. Когда сырьем стадии
собой экспериментальные продукты 13#, 14#, 15#,
DSO был продукт 2#, эксперименты проводились
соответственно, а катализатором был M при давле-
при температурах 240, 260 и 280°C, соответствен-
нии 2.0 МПа, объемной скорости 1.0 ч-1, темпера-
но; соответствующие продукты были последова-
турах реакции 350, 370, 390°C, соответственно; эти
тельно пронумерованы как 7#, 8#, 9#. Когда сы-
продукты были последовательно пронумерованы
Рис. 4. Схема расчета кинетических параметров.
НЕФТЕХИМИЯ том 61 № 3 2021
362
YING XIANG и др.
Таблица 1. Состав продуктов экспериментаa
No.
Продукт реакции
A, мас. %
O, мас. %
N, мас. %
iP, мас. %
nP, мас. %
ОЧИ
0#
Сырье
34.41
20.23
16.78
23.89
4.69
85.72
1#
M350
33.86
14.58
17.28
28.10
6.18
86.92
2#
M370
34.26
16.68
16.40
26.91
5.75
87.06
3#
M390
34.11
15.8
16.61
27.17
6.32
88.09
4#
M350-D240
33.94
11.87
17.87
29.82
6.51
86.28
5#
M350-D260
33.47
9.40
18.76
31.43
6.94
86.02
6#
M350-D280
33.33
6.96
19.63
32.92
7.16
85.96
7#
M370-D240
33.7
11.46
17.92
30.26
6.66
86.46
8#
M370-D260
33.82
9.91
19.20
31.06
6.01
86.67
9#
M370-D280
33.48
6.92
19.85
32.83
6.92
85.97
10#
M390-D240
34.13
12.81
17.31
28.88
6.88
87.27
11#
M390-D260
34.60
10.11
18.10
30.06
7.13
86.87
12#
M390-D280
34.61
8.19
19.12
31.05
7.03
86.29
13#
D240
33.02
14.23
17.35
28.72
6.68
84.14
14#
D260
33.02
12.72
17.83
29.42
7.01
84.06
15#
D280
32.68
10.28
18.50
31.28
7.25
83.84
16#
D240-M350
33.11
12.43
17.39
29.84
7.23
85.66
17#
D240-M370
33.24
12.15
17.27
29.80
7.54
86.05
18#
D240-M390
33.39
12.95
17.18
29.12
7.36
86.33
19#
D260-M350
33.18
11.32
17.67
30.59
7.26
85.25
20#
D260-M370
32.86
11.38
17.78
30.50
7.48
85.50
21#
D260-M390
33.45
10.07
17.99
30.88
7.61
85.69
22#
D280-M350
32.80
8.24
18.90
32.36
7.71
84.75
23#
D280-M370
32.94
9.31
18.53
31.51
7.71
84.90
24#
D280-M390
32.98
9.33
18.35
31.49
7.85
85.29
a ОЧИ: октановое число по исследовательскому методу; D240-M350 - продукт реакции после стадии DSO при температуре 240°C
и стадии M при температуре 350°C.
как 16#...24#. В табл. 1 приведены результаты ана-
чение. Процесс расчета кинетических параметров
лиза экспериментальных данных по углеводород-
показан на рис. 4.
ному составу для процессов M-DSO и DSO-M.
(2)
РАСЧЕТ КИНЕТИЧЕСКИХ ПАРАМЕТРОВ
СРАВНЕНИЕ КИНЕТИЧЕСКИХ ПАРАМЕТРОВ
Расчет кинетических параметров производили
ПРОЦЕССОВ DSO-M И M-DSO
путем программирования в MATLAB. С одной сто-
С использованием продуктов 0#, 1#, 2#, 3#, 4#,
роны, для решения дифференциальных уравнений
5#, 6#, 7#, 8# 9# из табл. 1 в качестве исходных дан-
использовали метод Рунге-Кутта четвертого поряд-
ных были получены кинетические параметры ста-
ка (из-за его высокой точности). С другой стороны,
дий M и DSO в процессе M-DSO, как показано в
для оптимизации целевой функции использовали
табл. 2 и 3. С использованием продуктов 0#, 13#,
генетический алгоритм (GA), поскольку он пред-
14#, 15#, 18#, 19#, 20#, 21#, 22#, 24# из табл. 1 в
ставляет собой алгоритм глобальной оптимизации,
качестве исходных данных, были получены кине-
позволяющий найти оптимальное решение, имити-
тические параметры стадий DSO и M в процессе
DSO-M, как показано в табл. 2 и 3.
руя естественный отбор и генетический механизм,
чтобы избежать попадания в цикл локальных зна-
Анализ данных табл. 2 показывает, что:
чений [13]. Целевая функция (func) - формула (2),
1) У констант скорости реакции значения k1
где yic - расчетное значение, yir - фактическое зна-
(превращение олефина в н-парафин), k2 (превраще-
НЕФТЕХИМИЯ том 61 № 3 2021
ДИСКРЕТНАЯ КИНЕТИЧЕСКАЯ МОДЕЛЬ
363
Таблица 2. Сравнение кинетических параметров стадии M процессов DSO-M и M-DSOа
k (константа скорости реакции), см-3·г·ч-1
стадия post-M
стадия pre-M
стадия pre-M
стадия post-M
(DSO260-M)
350°C
370°C
390°C
350°C
370°C
390°C
1
O→nP
15.64
3.42E+01
55.99
2.44E+04
1.67
1.84
2.01
0.49
0.69
0.95
2
O→iP
4.90
1.16E+01
50.22
2.20E+04
4.52
4.66
4.79
1.36
1.84
2.44
3
nP→iP
51.53
3.07E+04
46.86
1.27E+04
1.47
2.00
2.68
1.50
1.99
2.59
4
iP→nP
60.99
2.91E+04
57.17
1.90E+04
0.22
0.32
0.46
0.31
0.43
0.60
5
N→nP
68.92
4.35E+04
65.98
9.72E+04
0.07
0.11
0.16
0.29
0.43
0.62
6
N→iP
58.93
5.54E+04
59.77
5.26E+04
0.64
0.91
1.26
0.51
0.74
1.03
7
N→O
86.10
1.21E+06
86.24
1.12E+06
0.07
0.12
0.20
0.07
0.11
0.18
8
O→N
36.59
4.96E+02
58.75
4.66E+04
0.43
0.53
0.65
0.55
0.79
1.10
9
O→A
134.41
3.13E+09
130.97
9.29E+09
0.02
0.04
0.08
0.10
0.21
0.45
10
A→N
120.68
1.92E+09
120.57
2.03E+09
0.15
0.30
0.60
0.16
0.33
0.65
11
O→nP
230.70
9.98E+17
230.39
8.11E+17
0.05
0.18
0.67
0.04
0.16
0.58
а Направление реакции (A → B): реакция превращения дискретного вещества A в дискретное вещество B; стадия Pre-M: стадия M
в процессе M-DSO; стадия Post-M : стадия M в процессе DSO-M.
ние олефина в изопарафин) и k3 (изомеризация н-
пользован для снижения потери октанового числа
парафинов) намного больше, чем у других констант
на стадии pre-M.
на стадии M; это указывает на то, что эти три реак-
2) Константа скорости насыщения олефинов
ции являются основными реакциями, влияющими
(k1+k2+k8+k9) стадии pre-M намного выше, чем ста-
на стадию М в обоих процессах DSO-M и M-DSO.
дии post-M. Таким образом, функция уменьшения
Между тем, значение k2 на стадии pre-M (стадия M
содержания олефинов может быть улучшена для
в процессе M-DSO) намного больше, чем на стадии
стадии pre-M.
post-M (стадия M в процессе DSO-M). Причина
3) Что касается энергии активации, то энергия
заключается в том, что при поддержании концен-
активации превращения олефина в изопарафин
трации олефина на высоком уровне нормальный
на стадии pre-M намного меньше, чем на стадии
(изомерный) олефин с большей вероятностью по-
post-M, так что реакции изомеризации протекают
лучит протоны водорода из кислотных центров ка-
легче и могут улучшить качество бензина на ста-
тализатора с образованием вторичных (или трет-)
дии pre-M.
карбокатионов, а вторичные карбокатионы будут
4) Ароматические углеводороды имеют высокое
легко превращаться в трет-карбокатионы с более
октановое число, но реакции превращения олефи-
высокой стабильностью. Если трет-карбокати-
на и нафтена в ароматические углеводороды на
оны теряют протоны, они образуют изоолефин, а
стадии М протекают нелегко из-за больших значе-
затем изоолефин насыщается путем гидрирования
ний энергии активации и малых значений констант
в изопарафин с высоким октановым числом. Таким
скорости реакции. Следовательно, функция арома-
образом, функция изомеризации может быть улуч-
тизации катализаторов М нуждается в дальнейшем
шена, а катализатор M может быть полностью ис-
улучшении.
НЕФТЕХИМИЯ том 61 № 3 2021
364
YING XIANG и др.
Таблица 3. Сравнение кинетических параметров стадии DSO процессов DSO-M и M-DSOa
k (константа скорости реакции), см-3·г·ч-1
стадия post-DSO
стадия pre-DSO
стадия pre-DSO
стадия post-DSO
(M350-DSO)
240°C
260°C
280°C
240°C
260°C
280°C
1
O→nP
34.98
9.80E+03
44.8
5.89E+04
2.69
3.66
4.87
1.62
2.40
3.46
2
O→iP
30.92
1.20E+04
65.18
1.80E+07
8.54
11.21
14.43
4.17
7.40
12.60
3
nP→iP
51.44
2.15E+05
53.99
5.50E+05
1.25
1.96
2.98
1.40
2.28
3.57
4
iP→nP
55.95
2.49E+05
57.32
1.35E+05
0.50
0.82
1.30
0.20
0.33
0.52
5
N→nP
70.99
4.51E+06
59.67
3.07E+05
0.27
0.50
0.89
0.26
0.44
0.71
6
N→iP
68.99
6.90E+06
55.99
2.24E+05
0.65
1.20
2.11
0.45
0.73
1.16
7
N→O
87.40
5.90E+07
73.94
2.10E+06
0.07
0.16
0.33
0.06
0.12
0.22
8
O→N
56.18
8.59E+05
58.12
2.21E+06
1.64
2.69
4.25
2.68
4.47
7.18
9
O→A
142.7
4.90E+11
133.38
1.00E+13
0.00
0.01
0.02
0.27
0.86
2.54
10
A→N
95.90
3.05E+09
68.15
1.86E+06
0.53
1.23
2.68
0.22
0.39
0.68
a Стадия Pre-DSO: Стадия DSO в процессе DSO-M; Стадия Post-DSO: Стадия DSO в процессе M-DSO.
Анализ данных табл. 3 показывает, что:
3. Для энергии активации реакций: E2 < E1 <
< E3 < E4 < E7 < E6 < E5 < E8 < E10 < E9. Можно ви-
1. У констант скорости реакции значения k1, k2,
деть, что реакция превращения олефинов в арома-
k3 и k8 (превращение олефина в нафтен) намного
тические соединения является наиболее трудной.
больше, чем у других констант на стадии DSO, что
указывает на то, что насыщение олефина и изоме-
4. Стадия pre-DSO имеет немного большую
ризация н-парафина являются основными реакция-
константу скорости насыщения олефина гидриро-
ми, влияющими на стадию DSO. Между тем, неко-
ванием (k1+k2+k8), чем стадия post-DSO. Это ука-
торые олефины насыщаются до н-алканов с самым
зывает на то, что ее воздействие на уменьшение со-
низким октановым числом, что вызывает потерю
держания олефина немного сильнее, чем у стадии
октанового числа.
pre-DSO. Однако константы скорости насыщения
олефинов на стадии pre-М намного больше, чем
2. Стадия рre-DSO (стадия DSO в процессе
на стадии post-М. Таким образом, процесс M-DSO
DSO-M) имеет гораздо меньшую константу ско-
сильнее влияет на снижение содержания олефинов
рости реакции превращения олефина в аромати-
ческие углеводороды и гораздо более высокую
в бензине FCC.
константу скорости реакции превращения арома-
В целом, сравнивая кинетические параметры
тических углеводородов в нафтен, чем стадия post-
стадий M и DSO, можно видеть, что необходимо
DSO (стадия DSO в процессе M-DSO), что ускоря-
иметь высокое содержание олефинов в сырье для
ет расход ароматических углеводородов с высоким
изомеризации олефинов на стадии M. Для стадии
октановым числом, в то же время стадия pre-DSO
post-M олефины в сырье будут сильно насыщены
имеет гораздо меньшую энергию активации пре-
после стадии DSO, что не позволяет очистить бен-
вращения олефина в изопарафин, чем стадия post-
зин FCC на стадии M и повысить октановое число
DSO. В целом, процесс DSO-M ведет к потере ок-
бензина. Между тем, для стадии pre-M содержание
танового числа бензина FCC.
олефинов уменьшилось, так что процесс M-DSO
НЕФТЕХИМИЯ том 61 № 3 2021
ДИСКРЕТНАЯ КИНЕТИЧЕСКАЯ МОДЕЛЬ
365
Таблица 4. Сравнение фактических и расчетных значений для процесса M-DSO
10#
11#
12#
A
34.13
34.00
-0.37
34.60
34.04
-1.61
34.61
34.26
-1.01
O
12.81
12.22
-4.60
10.11
10.31
1.95
8.19
7.85
-4.15
N
17.31
17.54
1.36
18.10
17.99
-0.60
19.12
18.48
-3.34
iP
28.88
29.3
1.46
30.06
30.52
1.54
31.05
32.07
3.29
nP
6.88
6.94
0.88
7.13
7.14
0.20
7.03
7.35
4.49
Таблица 5. Сравнение фактических значений и расчетных значений для процесса DSO-M
16#
17#
23#
A
33.11
32.88
-0.69
33.24
32.81
-1.31
32.94
32.44
-1.50
O
12.43
12.84
3.30
12.15
12.38
1.89
9.31
8.97
-3.67
N
17.39
17.23
-0.95
17.27
17.21
-0.36
18.53
18.17
-1.93
iP
29.84
29.93
0.32
29.80
30.29
1.66
31.51
32.61
3.48
nP
7.23
7.12
-1.56
7.54
7.31
-3.00
7.71
7.80
1.14
в целом превосходит процесс DSO-M в этом от-
ВЫВОДЫ
ношении. Кроме того, по сравнению с процессом
1. На основе механизма гидрирования нефти и
DSO-M, в процессе M-DSO образуется больше
характеристик процессов M-DSO и DSO-M были
ароматических углеводородов и повышается окта-
разработаны 5-компонентные реакционные схемы
новое число бензина.
стадии M и стадии DSO, соответственно.
2. Кинетические параметры стадий M и DSO
ПРОВЕРКА МОДЕЛИ
в процессах M-DSO и DSO-M были рассчитаны,
соответственно, по экспериментальным данным
С использованием полученных выше кинетиче-
с тяжелой фракцией бензина FCC в микрореакто-
ских параметров была проведена проверка моделей
ре с неподвижным слоем. Показано, что процесс
процессов M-DSO (табл. 1, 10#, 11#, 12#) и DSO-M
M-DSO лучше, чем DSO-M в плане уменьшения
(табл. 1, 16#, 17#, 23#). Результаты проверок пред-
содержания олефинов и повышения октанового
ставлены в табл. 4 и 5. Относительные погрешно-
числа с точки зрения динамики.
сти между расчетным и фактическим значениями
для моделей находятся в пределах 5%, что указыва-
3. Оценка параметров моделей показывает, что
ет на превосходную экстраполяцию и надежность
относительные погрешности между расчетным и
моделей. Кроме того, эти модели позволяют дать
фактическим значениями составляют менее 5%.
дополнительные рекомендации по моделированию
Это говорит о том, что дискретные кинетические
и оптимизации процесса гидроочистки бензина FCC.
модели, установленные для процессов M-DSO и
НЕФТЕХИМИЯ том 61 № 3 2021
366
YING XIANG и др.
DSO-M, имеют хорошую экстраполяцию и надеж-
5.
Wang B.C., Huo D.L., Cui D.Q., Cui J.W., Liu Y.L.,
ность и позволяют давать рекомендации по моде-
Zhang X.J. Industrial regeneration of hydro-upgrading
лированию и оптимизации процессов M-DSO и
catalyst TM2-DSO process // Mod. Chem. Ind. 2010.
DSO-M.
V. 30. P. 59-61. https://doi.org/10.16606/j.cnki.
issn0253-4320.2010.08.012
6.
Lan L., Ju Y.N. Development and commercial
ИНФОРМАЦИЯ ОБ АВТОРАХ
trial of DSO technology for FCC gasoline hydro
Ying Xiang, ORCID: http://orcid.org/0000-0002-
desulfurization // China Pet. Process. Petrochem.
5619-9720
Technol. 2010. V. 41. P. 53-56. https://doi.org/10.3969/j.
Jiaying Shen, ORCID: http://orcid.org/0000-0002-
issn.1005-2399.2010.11.010
6906-6484
7.
Ouyang F.S., Ling Q., Yu Z.K. Reaction kinetic model for
flexible dual-riser catalytic cracking process // J. Chem.
Fusheng Ouyang, ORCID: http://orcid.org/0000-
Eng. Chin. Univ. 2015. V. 29. P. 1106-1113. https://doi.
0001-7523-6673
org/10.3969/j.issn.1003-9015.2015.00.025
8.
Weekman V.W.J., Nace D.M. Kinetics of catalytic
КОНФЛИКТ ИНТЕРЕСОВ
cracking selectivity in fixed, moving, and fluid bed
Авторы заявляют об отсутствии конфликта
reactors // AIChE J. 1970. V. 16. P. 397-404. https://doi.
интересов, требующего раскрытия в данной статье.
org/10.1002/aic.690160316
9.
Zhu B.C., Weng H.X., Zhu Z.B. Catalytic reaction
engineering. Ixueshu, Beijing, 1999. ISBN:
СПИСОК ЛИТЕРАТУРЫ
9787800438547.
1. Lappas A.A., Valla J.A., Vasalos I.A., Kuehler C., Francis J.,
10.
Jiang H., Huang S. Eight-Lump reaction kinetic model
O’Connor P., Gudde N.J. The effect of catalyst properties
for the maximizing isoparaffin process for cleaning
on the in situ reduction of sulfur in FCC gasoline // Appl.
Catal. A. 2004.V. 262. P. 31-41. https://doi.org/10.1016/j.
gasoline and enhancing propylene yield // Energy Fuels.
apcata.2003.11.014
2016. V. 30. P. 10770-10776. https://doi.org/10.1021/acs.
2. Dong H.M., Qu Y., Sun L.L. Application of Prime-G+
energyfuels.6b02208
technology in FCC naphtha hydrodesulfurization
11.
Fusheng O., Yongqian W., Qiao L. A lumped kinetic
unit // China Pet. Process. Petrochem. Technol.
model for heavy oil catalytic cracking FDFCC process //
2012. V. 43. P. 27-30. https://doi.org/10.3969/j.
Pet. Sci. Technol. 2016. V. 34, P. 335-342. https://doi.org/
issn.1005-2399.2012.11.006
10.1080/10916466.2015.1132236
3. Ghosh P., Andrews A.T., Quann R.J., Halbert T.R.
12.
Zong G., Ning H., Jiang H., Ouyang F.S. The lumping
Detailed kinetic model for the hydro-desulfurization
kinetic model for the heavy oil catalytic cracking MIP
of FCC Naphtha // Energy Fuels. 2009. V. 23. P. 5743-
process // Pet. Sci. Technol. 2010. V. 28. P. 1778-1787.
5759. https://doi.org/10.1021/ef900632v
https://doi.org/10.1080/10916460903261749
4. Yang Q.J., Pan D., Zhang N.N. Industrial application
13.
Ge J.K., Qiu Y.H., Wu C.M., Pu G.L. Summary of
of catalytic gasoline M-DSO combined with solvent
extraction and desulfurization // Petrochem. Technol.
genetic algorithms research // Appl. Res. Comput.
Appl. 2018. V. 36. P. 343-346. https://doi.org/10.3969/j.
2008. V. 25. P. 2911-2916. https://doi.org/10.3969/j.
issn.1009-0045.2018.05.016
issn.1001-3695.2008.10.008
НЕФТЕХИМИЯ том 61 № 3 2021