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

книги / Моделирование физико- химических процессов нефтепереработки и нефтехимии

..pdf
Скачиваний:
29
Добавлен:
19.11.2023
Размер:
31.95 Mб
Скачать

Рис. IX-11. Влияние концентрации кислорода в газе* поступающем в пер­ вую зону регенератора, на результраты процесса в ней:

I — температура; 11 — содержание кокса на катализаторе (для режимов 1, 711, охарактеризованных в табл. IX -6),

Рис. IX -12. Влияние температуры кислородсодержащего газа, поступающего в первую вону, на результаты процесса в регенераторе:

1 — температура; I I — содержание кокса на катализаторе (для режимов I, 12—16, охарактеризованных в табл. IX -6).

6,4 т/ч при температуре

27 °С, плотность поверхности охлажда­

ющих

змеевиков в III

зоне 8,05 м2/м. Расчеты выполнены на

ЭВМ

М-20, результаты

расчетов представлены на рис. IX -10 —

IX -14 [27, 28].

Изменение расхода кислородсодержащего газа в I Бону реге­ нератора в пределах 3,4 —8,4 т/ч (рис. IX -10) или концентрации

зоножйароэ ознон он H jB u t/D acn 'adom % )('зэон

JeMnepamypat

Длина по оси регенератора, м

Рис. IX-13. Влияние расхода воздуха, поступающего во вторую зону, на результаты процесса в регенераторе:

1 — температура; I I — содержание кокса на катализаторе (длл режимов 1, 1721% охарактеризованных в табл. IX -6).

кислорода в нем от 15 до 22,5% (масс.) (рис. IX -11) при неизмен­ ных в каждом случае прочих показателях существенно влияет лишь на температуры в I и II зонах. При уменьшении расхода газа в I зону регенератора или концентрации кислорода в нем «горячая точка» регенератора и зоны с равной закоксованностыо катализатора смещаются к выходу из аппарата. При этом темпе­ ратура в «горячей точке» понижается. Максимальное смещение «горячей точки» и зон с равной закоксованностыо, для выбранных условий (режим 1 и 6, а также 7 и 11) в обоих случаях составляет около 2 м (высота одной секции), а понижение температуры — 20 °С. Однако такое снижение интенсивности процесса не приво­

з а

дит к нарушению установленного технологическим регламентом содержания остаточного кокса на катализаторе. Так, для режи­ мов б и 11 эта величина составляет 0,03% .

Изменение температуры кислородеодержащего газа на входе в I зону регенератора на 100 °С при неизменных остальных пока­ зателях существенно не влияет на поля концентраций и темпера­ тур в регенераторе (рис. IX -12). При расходе воздуха во II зопу 8,2 т/ч (режим 17) содержание кокса на катализаторе на выходе

Длина по оси регенераторам

Рис. IX-14. Влияние плотности поверхности включенных во вторую зону змеевиков на результаты процесса в регенераторе:

1 — температура; I I — содержание кокса на катализаторе (для режимов 19 22—26, ■охарактеризованных в табл. IX -6).

из этой

зоны составляет (рис. IX-13)

около 0,5% (масс.), а при

расходе

16,0 т/ч

(режим 21) — 0,14%

(масс.).

Температура

в «горячей точке» в

регенераторе повышается

на 20 °С и «горячая точка» перемещается от выхода из регенера­ тора ко входу в III зону.

Расчеты показывают, что регулирование расхода воздуха во II вонуэффективно для управления процессом, так как позволяет 'Снизить остаточное содержание кокса на катализаторе в случае роста коксовой нагрузки регенератора.

Использование газов и8 III зоны регенератора, богатых кисло­ родом, в верхних секциях позволяет увеличить нагрузку на

II зону регенератора и выжигать больше кокса в этой зоне при неизменном общем расходе воздуха. В этом заключается преиму­ щество регенератора с линией перехода газов. Влияние интенсив­ ности теплоотвода к охлаждающим змеевикам (режимы 22— 26), включенным в средней зоне регенератора, показано на рис. IX -14.

С изменением плотности поверхности охлаждающих змеевиков изменяются только температуры в средней и нижней зонах, поло­ жение «горячей точки» не меняется. Влияние на количество выго­ рающего кокса в этом случае почти отсутствует, что объясняется диффузионным характером процесса в этих зонах.

В заключение следует отметить, что моделирование процессов регенерации может быть эффективно осуществлено и для других конструкций аппаратов, например, с «кипящим» или восходя­ щим слоями контактного материала. Математические описание и моделирование возможны после изучения режима перемешива­ ния в этих аппаратах методами, описанными в главе III.

Кроме того, существенное влияние на процессы регенерации может оказать и неустойчивость режима, которая отмечена экспе­ риментально при разогревах, превышающих 30 °С [10, 29]. При высоких разогревах необходим переход к нестационарному описанию или проведение процесса в устойчивой области, для чего следует использовать методы, рассмотренные в главе V.

ЛИТЕРАТУРА

1 .

Панченков Г. М ., Голованов

Н. В . И з в .

А Н

С С С Р .

О Т Н , 1 9 5 1 , № 1 0 ,

 

с . 1 9 1 3 - 1 9 1 8 ; 1 9 5 2 , № 8 , с . 1 0 3 1 - 1 0 3 7 .

 

 

 

 

7 , с .

 

 

2 .

Добычин Д . И ., Клибанова Ц . М. Ш

Ф Х ,

1 9 5 8 ,

т . 8 ,

1 0 2 3 — 1 0 5 9 .

3 .

Prater С. D. C h e m . E n g . Sci.,

1 9 5 8 ,

v . 8 , №

2 ,

p .

2 4 8 — 2 5 5 .

 

 

4 .

Bondi C. e . a . I n d . E n g . C h e m .

P r o c .

D e s .

D e v .,

1 9 6 2 ,

v . 1 ,

1 , p .

1 9 6 —

 

202.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 .

Лавров И. В ., Шурыгин А. П. В в е д е н и е в

т е о р и ю г о р е н и я

и г а з и ф и к а ц и и

6.

т о п л и в а . М ., И а д -в о А Н С С С Р , 1 9 6 2 . 2 1 5 с .

 

 

 

 

 

 

 

Ж аров Ю. М. И с п о л ь зо в а н и е

м е т о д о в х и м и ч е с к о й

к и н е т и к и д л я

р а з ­

 

р а б о т к и , м о д е л и р о в а н и я и о п т и м и за ц и и п р о ц е с с о в н е ф т е п е р е р а б а т ы в а ­

 

ю щ е й п н е ф т е х и м и ч е ск о й п р о м ы ш л е н н о ст и . А в т о р е ф . д о к т . д и с с . М .,

 

М И Н Х и Г П , 1 9 6 8 . 3 5 с .

 

 

 

 

 

 

 

 

 

 

 

 

7 .

П ат енков Г . М ., Жоров Ю. М. « Н е ф т е х и м и я , п е р е р а б о т к а

н е ф т и и

г а з а » .

8.

Т р у д ы М И Н Х и Г П . В ь ш ., 2 8 . М ., Г о с т о п т е х и з д а т , 1 9 6 0 , с . 3 — 1 0 .

 

Hoge W. А. а .

о . P e t r o l . E n g . ,

1 9 4 9 ,

1 1 ,

р .

1 4 9 —

1 5 6 .

 

 

 

9 .

Галинин Е . Т., Левинтер М. Е .

В

с б . : « Т е к у щ а я и н ф о р м а ц и я п о н е ф т е ­

 

п е р е р а б о т к е

и н еф тех и м и и » .

М .,

Ц Н И И Т Э н е ф т ь ,

1 9 6 4 , №

2 4 . 7 с .

1 0 .

Sam path В ., Ram achandrah Р. , Hughes R . C a n a d .

J . C h e m .

E n g . ,

1 9 7 5 ,

 

v . 5 3 , № 2 , p . 1 8 4 — 1 9 1 .

 

 

 

 

 

 

 

 

 

 

 

 

1 1 .

Deemter J .

I n d .

E n g . C h e m .,

1 9 5 4 , v . 4 6 ,

1 1 , p .

2 3 0 0 — 2 3 0 9 .

 

1 2 .

Johnson В .

M. ,

Froment G.

F ., Watson

С.

C. C h e m . E n g . ,

1 9 6 2 ,

v . 1 7 ,

5 , p . 8 3 5 — 8 4 2 .

1 3 .

Shulm an B.

Z. I n d . E n g . C h e m .

I n t . E d . , 1 9 6 3 ,

v .

5 5 , № 1 , p . 4 4 — 5 3 .

1 4 .

Лдельчик И .

К . С п р а в о ч н и к п о

г и д р а в л и ч е с к и м

с о п р о т и в л е н и я м .

М . ,

 

Г о с э н е р г о и з д а т , 1 9 6 0 . 1 5 6 с .

 

 

 

 

1 5 .

Иоффе И . И ., Письмен Л. М. И н ж е н е р н а я х и м и я

г е т е р о г е н п о г о

к а т а ­

л и з а . И з д . 2 - е . Л . , « Х и м и я » , 1 9 7 2 . 4 6 2 с .

334

16.Панченков Г . М ., Лазъян Ю. И ., Жаров Ю. М. ЖФХ, 1967, т. 41, № 9,

с.2276—2280.

17.Жаров 10. М ., Панченков Г . М ., Лазъян 10 . И. ЖФХ, 1967, т. 41, № 7,

с.1574—1578.

18.Панченков Г . М ., Л азъян Ю. И ., Жоров Ю. М . — В кн.: Математические описания и оптимизация процессов нефтепереработки и нефтехимии.

Труды МИНХ и ГП, выл. 74. Л., «Химия», 1967, с. 77—83.

19. K aranth N ., H ughes R . Chiem. Eng. Sci., 1974, v. 29, № 2, p. 197—204.

20.M assoth F . E . Ind. End. Chem. Proc. Des. Dev., 1967, v. 6, № 2, p. 200— 206.

21.Smith R ., Dresser T. Petrol. Ref., 1957, № 7, p. 196—203.

22.Опришко А. А . и др., «Химия и технология топлив и масел», 1970, № 3,

с. 37—40.

23. Hagerbaumer W ., Lee R . Petrol. Process., 1947, № 2, p. 308—315.

24.Америк Б . К . и др. Труды ГрозНИИ. Вып. 4, М., «Недра», 1959, с. 16—29.

25.Аделъсон С. В ., Заитова А. Я . «Химия и технология топлив и масел», 1964, № 5.

26.Жоров Ю. М . и др. '— В кн.: Математическое описание и оптимизация процессов переработки нефти и нефтехимии. Труды МИНХ и ГП, вып. 74. Л ., «Химия», 1967, с. 88—96.

27.Панченков Г . М . и др. «Химия и технология топлив и масел», 1969, № 7,

с. 4—9.

28. Панченков Г . М . и др. «Нефтепереработка н нефтехимия», 1970, № 6,

с> д_11'92

29. Sam path В ., H ughes R . Process Technol., 1973, v. 18, № 1-2, p. 39—43.

Глава X

Моделирование

каталитических процессов

нефтепереработки

инефтехимии

Вэтой главе рассмотрены математические описания и резуль­ таты моделирования крупнотоннажных каталитических •процес­ сов. Выбор процессов связан как с их большим техническим значе­ нием, так и с их «представительностью»; при их моделировании удобно продемонстрировать применение математических описаний для аппаратов различных типов при переработке индивидуальных

веществ (или «групповых» компонентов) и технических смесей.

1. МОДЕЛИРОВАНИЕ КАТАЛИТИЧЕСКОГО РИФОРМИНГА БЕНЗИНОВЫХ ФРАКЦИЙ

Каталитический риформинг бензиновых фракций на платино­ вом катализаторе (платформинг) — ведущий технический процесс для получения высокооктановых бензинов и ароматических углеводородов. Сырьем являются обычно фракции прямогонных ■бензинов, содержащие парафиновые, нафтеновые и ароматические углеводороды и небольшое количество олефинов. В сырье присут­ ствуют также, как микропримеси, различные элементоорганиче- ■ские соединения и вода. Процесс проводится при температурах -около 500 °С и давлениях 1—4 МПа с разбавлением сырья водо­ родсодержащим газом до мольного соотношения «водород/сырье», равного 5 —8. Обычно его осуществляют в системе из трех после­ довательно соединенных адиабатических реакторов с неподвиж­ ными слоями катализатора. Между реакторами происходит подо­ грев продукта.

Разработка математического описания платформинга пред­ ставляет большой интерес для целей оптимального проектирования и управления процессом. Существенно отметить, что статисти­ ческие описания [1, 2] не могут быть эффективно использованы для решения задач проектирования. Однако статистические методы могут быть полезны при переходе от рассчитываемых на •основании дифференцированного описания физических характе­ ристик (состав продукта) к техническим (октановое число) [21.

336

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

В работе [3] процесс описывается 53 уравнениями первого порядка, учитывающими превращения 20 компонентов реакцион­ ной смеси. Однако точный анализ сырья и катализата с определе­ нием содержания 20 веществ чрезвычайно сложен и неудобен как в промышленных, так и в лабораторных условиях. Неточ­ ность или ошибка в определении хотя бы одного из компонентов существенно сказывается на результатах. Кроме того, анализ 20 компонентов требует значительного времени, в течение кото­ рого качество сырья может существенно измениться.

Математическая модель платформинга в виде системы пз 8 дифференциальных уравнений предложена Смитом [4]. Он объединяет химические превращения в 6 обобщенных реакций, определяющих состояние системы: нафтены р ароматические + -{-ЗН 2, парафины р нафтены + Н 2, гидрокрекинг нафтенов, гид­ рокрекинг парафинов. Недостатком работы Смита является то, что для обеспечения сходимости результатов расчета с эксперимен­ тальными данными используются произвольно выбранные значе­ ния теплот и констант равновесия, отличающиеся от точных. Описание Смита лишь качественно характеризует процесс, вслед­ ствие неточного определения постоянных коэффициентов и непра­ вильной записи кинетических уравнений (оно неверно учитывает влияние давления на результаты).

Описание платформинга, позволяющее решать технические

задачи, получено

в СССР [5—8].

В соответствии

с методикой, изложенной в главе V, сырье

рассматривается как смесь некоторых обобщенных углеводородов; естественно положить, что каждый из них содержит среднее для сырья число углеродных атомов. Такое допущение существенно упрощает кинетическую схему, сводя число рассматриваемых реакций к относительно небольшой величине. Вместе с тем, оно дает оценку процесса по выходу ароматических углеводородов, содержаншо в сырье нафтенов и парафинов, т. е. по групповому, а не индивидуальному составу сырья и катализата. Олефиновые углеводороды, содержание которых в сырье и катализате невелико, целесообразно при дальнейшем анализе объединить с парафинами.

Полную кинетическую схему платформинга можно записать следующей системой реакций (для краткости записи водород,

участвующей в реакциях, в схеме не указывается):

 

 

W1

 

Wo

Wz

 

 

 

NG

а

 

No

 

2 р '

N. - ^ 2 p'

 

 

 

 

WA

 

 

 

(X .l)

Р

Шв

г‘Р

Ws

N 5

oilо

Wl2

гР

2 р'

N6

5

 

 

ЬУ7

 

 

W9

N

 

 

 

 

 

U>i3

 

здесь P — нормальный парафиновый углеводород; iP — изопара­ финовый углеводород; Nc — шестичлепный нафтеновый углево-

3 3 7

дород; N5 — пятичленный нафтеновый углеводород; А — арома­

тический углеводород; 2 ^ — газообразные парафины; w — ско­ рость реакции.

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

Вместе с тем, можно предположить, что удовлетворительное соответствие с промышленными данными может быть получено и для более простых схем, представляющих частные случаи общей

схемы (Х .1).

Нами исследованы

схемы

[9]:

 

2 Р '

<

Р — ►А «<— N

2

Р ' ^—

Р

N — > А

2 Р '

«—

P ? ± N ? ± A

2

р/ «— - Р ^

N =<=* А

Как отмечено выше (стр. 79,80), исследование можно произвести, используя некоторые очевидные для каждой схемы связи, не требующие определения вида кинетических уравнении. Так, например, из системы уравнений материального баланса для пер­ вой из простых схем следует, что одно из уравнений можно про­ интегрировать независимо, после чего получим для изменения чисел мольных потоков щ Ещ\

1 -Д » р= ДиА-Д « к + т Д)г2 р#

Здесь б — число молей сухого газа, образующегося из одного моля парафина; оно должно быть постоянно для опытов с сырьем одного вида. Проверка постоянства б требует знания лишь точных материальных балансов работы установки (для расчета величин Ап) и выполняется сравнительно просто. Кроме того, независимо может быть проинтегрировано и проверено уравнение теплового баланса (см. стр. 66). Сопоставление рассчитываемых и экспери­ ментальных перепадов температур позволяет судить об обоснован­ ности схемы.

Мы убедились [9], что вполне приемлемая для расчетов точ­ ность достигается при использовании схемы:

N

&i

Wt

НУа

А + З Н ,

Р * = £ N+ На

Р + бИ2 —

Р ' • . • (Х .2)

 

 

Wt

 

 

Более простые схемы не обеспечивают удовлетворительной точ­ ности описания, более сложная схема не улучшает заметно совпа­ дения расчетных и экспериментальных значений. Поэтому мы пользовались преимущественно схемой (Х .2). Для ее применения необходимо определить теплоты реакций АН , вид кинетических зависимостей для w и кинетические коэффициенты к 0 и Е .

.338

Из данных [10] следует, что теплоты реакций для правращений сложных смесей при платформинге слабо зависят от молеку­ лярной массы сырья. Для практических расчетов в промышлен­ ных условиях можно пренебречь зависимостью теплоты от темпе­ ратуры. Процессы каталитического риформинга проводят при 450—550 °С, поэтому для расчетов использованы теплоты реак­ ций при 500 °С.

Усреднение мольных теплот реакций различных углеводородов приводит к следующим значениям (в кДж/моль):

A7/ N -V A = 214

АЯР^ = 4 1 ,8

ДЯр^2р ,= -56,2

Знание АН и кинетической схемы

процесса позволяет пред­

ложить простой метод расчета теплот, основанный на.независи­ мом интегрировании уравнения теплового баланса ([10], см.

также главу

 

II). Так для схемы

(Х .2)

имеем:

 

dnA

 

dnN

 

chip

 

w b

(W i— ,u>6) +

(w 2wi)

~dV~ ~

( !t'2 — wa

G c d T

~~ A ^ N - VA (*^1 wb) ~Ь A ^P->-N

wi) + Д Я р ^ 2

 

d V

p , •Ш3

 

 

 

 

 

Выразив Wiw6, w2— Wi и w3 из первых трех уравнений и подста­ вив их в четвертое, получим уравнение для определения теплоты процесса при 500 °С (в кДж/кг):

Az,

 

^ZN

<?5оо = 3 1 2 . ю з

 

м т

М Ак

 

+ 56,2-103

Ягр

Z O P \

м Рк

^Рк )

 

где Я — выход катализата; zLи zoi — содержание г'-го компонента

в катализате и сырье, масс,

доли; М {к и Mic — молекулярные

массы

г-го компонента в катализате и сырье.

По

имеющимся данным,

теплота платформинга различных

технических видов сырья может меняться от 380 до 1200 кДж/кг, в зависимости от содержания нафтенов в сырье.

Процесс платформинга проводится при сильном разбавлении реагирующей смеси водородом. Кроме того, платина активно сорбирует водород и, следовательно, адсорбционный коэффициент водорода больше адсорбционных коэффициентов других веществ.

Поэтому для реакций N А, А ->• N, Р

N, скорости

w/ запи­

сываются в виде:

 

 

 

Wj— lcQ/e-E jjR T К н с х

, ^

Г си сх

(Х.4)

ilKs

4 ПНг

 

Скорости реакций гидрокрекинга, протекающих с разрывом С—С связи, пропорциональны объемной концентрации исходных веществ.

3 3 9

Математические оппсапия технических процессов, используемые для их моделирования

Платформинг бензина

Г л д р о к р е к и п г бензина

Г и д р о к р е к и н г г а з о й л я

по схеме:

по схем е:

п о с х е м е :

A

■шв N ----

N — н-Рп

uao~Vn С — ч- v3Д + V4O

 

WA y

 

 

wt

Ws

 

«-P» — vltl( P i- P 3)+ 0 — V

vi,ir +

Wt

 

 

— p

V

P'

4 V2 , 1 P 4 + V3 , 1 P 5 +

+ V2,I B + V3,XA

v*

 

+ V4 , l P e

 

 

 

Д — -

Vl, 2 r + V 2.2I

 

 

 

Ws

 

 

 

 

 

и з о - V n --------* v 1, 2 X

Г и д р о о ч и с тк а д и зел ьн о й ф р а кц и и п о схем е :

С — - V3A + V40;

О- — ► V i i r + v 2i B + v 3i A

д -

-v V i2r -J- V 22^

 

 

v3Sfl + v 4 S o

S G " " U I ^ C + H 2S

dnt

 

 

2 v'iW‘

 

d V

 

 

 

 

 

 

0

 

 

m (0) =

•Щ-0

 

 

 

щ

о

 

 

d T

 

 

 

 

 

A

=

 

-

±

\

A E

d V

 

 

 

G c

L

 

x (u>i—W5) + AHp^ ffX

X (

u>i)~\~Aff-p

ур/М'з

T (0 ) = Zo

0 4 - 0 4 A*A + 0 4 N . N +

+ 0 4 p 2 p - b 22

X (P i - P 3)+ v 2>2P4+ + V3 oP5 + y4 ,2 в

d n ;

d V ■ 2 v/w nu (0) = niln

dr

 

=■ {AHytOy-j-

d V

G c

- j- Д Я 2Ш24 Д Я 3Ш3

4 -

Д Я 4Ш4)

G‘' f r = 2 v^,“’,

3

Уд (0) = V3

У о ( ° ) = v 4

рг (0) = рБ (°) = °

d T

1

ЙУ

Gc

Х[Д Я 2ш2+Д//з^з]

- Vl , l S r +

 

 

+

V2 , 1 S B +

S 0 "

 

+

V3 , 1 S 4

 

 

- 3

*

0 + H

2S

-

v i , 2 S r

+

 

 

+

V2 , 2 S B

J2

-+

Д + H o S

0 J E l * - Г + H o S

 

b r

 

 

 

 

O

 

B - b H 2 S

 

Ь Б

 

 

 

 

ТАБЛИЦА Х-1

К а т а л и т и ч е с к и й к р е ­ к и н г га з о й л я п о схе м е :

А— - v 11A r +

+ V 2 , 1 A B + V 3 , 1 А ф 4

+V 4 , 1 A K

А Б

Wt

v i , 2 A r +

4

V2 , 2 A <I> + V 3 2 A K

А Ф

 

v l f S A r +

+ V 2 , 3 A B 4 ' V 3 3 A K

A

 

W *

V J I 4A b 4

Ar

 

4

V 2 , 4 A K

 

T l (0 )= T io

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