УДК 517.977.5
МОДЕЛИРОВАНИЕ ПРОЦЕССА ПИРОЛИЗА ПРОПАНА И ПРОПАН-ПРОПИЛЕНОВОЙ ФРАКЦИИ В СРЕДЕ МА^АВ
© И. М. Губайдуллин1*, А. Е. Мусина2, Р. В. Узянбаев3
'Институт нефтехимии и катализа УФИЦ РАН Россия, Республика Башкортостан, 450075 г. Уфа, пр. Октября, 141.
2Уфимский государственный нефтяной технический университет Россия, Республика Башкортостан, 450062 г. Уфа, ул. Космонавтов, 1.
3Национальный исследовательский Мордовский государственный университет им. Огарева Россия, Республика Мордовия, 430005 г. Саранск, ул. Большевистская 38.
Тел.: +7 (347) 242 0712.
*Етай: irekmars@mail.ru
В статье описана методика моделирования процесса пиролиза пропана, основанная на изменении концентраций компонентов во времени. Модель представляет собой систему обыкновенных нелинейных дифференциальных уравнений, решение которых позволяет получить полную картину о протекании процесса. С помощью программы, написанной в среде МАПАВ, был произведен поиск оптимального температурного режима и начального состава сырья процесса пиролиза.
Ключевые слова: пиролиз, математическое моделирование, пропан, пропилен, компактная модель.
Введение
Процесс пиролиза в современной нефтепереработке занимает одно из главных положений, как по масштабам перерабатываемого сырья, так и по объемам получаемых продуктов. Одними из главных продуктов процесса являются олефины (этилен, пропилен, ацетилен, бутадиен, бутилен), которые являются ценным сырьем для многих процессов нефтехимии.
В связи с доминирующим положением процесса пиролиза в современной нефтепереработке, существует необходимость его дальнейшего развития и совершенствования. Наиболее главными путями являются расширение сырьевой базы и изменение конструктивного оформления процесса.
Теоретический анализ
В настоящее время более половины перерабатываемого сырья составляет прямогонный бензин (54%) и газойль (6%), остальное — газы с НПЗ и конденсатных месторождений (этан — 29%, пропан — 6%, бутан — 5%) [1—4]. Таким образом, в связи с большими запасами природного газа, стоит вопрос о вовлечении в процесс пиролиза большего количества газов, таких как этан, пропан, бутан, а также их смесей.
Совершенствование конструктивного оформления процесса является путем, позволяющим повысить селективность процесса. При изучении литературы по данному направлению, был найден способ интенсификации процесса, позволяющий проводить процесс не в привычной печи пиролиза. В институте катализа имени Г. К. Борескова СО РАН рассматриваются способы создания экономически эффективной технологии переработки для невысоких производительностей путем управления химическими реакциями пиролиза через дополнительную генерацию радикалов в реагирующем газе посредством лазерного излучения, поглощаемого самой смесью реагентов [5—8]. Для этого используется реактор проточного типа с использованием микроволнового излучения. Микроволновое излучение позволяет проводить равномерный нагрев во всем объеме вещества, в результате проведения эксперимента в реакторе с микроволновым излучением были получены экспериментальные данные выхода целевых продуктов при различной температуре. [9—10]. Микроволновое излучение позволило повысить селективность процесса, за счет того, что процесс стал проходить при более высоких температурах.
Методика расчетов
В рамках поиска перспективного пути развития процесса, был исследован процесс пиролиза пропана. Расширенная модель составляет около 160 стадий и десятки веществ, в то время как для достаточно точного описания процесса достаточно гораздо меньшего количества стадий. За основу была взята компактная 30 стадийная модель пиролиза пропана, полученная Л. Ф. Нурисламовой методом локального анализа чувствительности [11]. 30-стадийная модель пиролиза пропана приведена в табл. 1. При исследовании механизмов сложных процессов нефтепереработки и нефтехимии на начальном этапе моделируются только химические превращения. Предполагается, что химическая реакция протекает в закрытой системе, без изменения объема. Такая система адекватно описывается законом действующих масс в виде системы обыкновенных нелинейных дифференциальных уравнений, задающихся задачей Коши.
¿с 0
— = / (с), с(0) = с0 ,
¿т
где с — концентрация компонента.
Таблица 1
30-стадийная модель пиролиза пропана
А, 1/с или Е, n
№ Стадия см3/(моль-с) кДж/моль
1 C3H —► C2H5^ + СНз^ 2.78E+18 376 -1.8
2 C2H5^ + СНз^—► C3H 2.83E+13 0 -0.5
3 C2H5^ — C2H4 + H^ 4.31E+12 155 1.19
4 C2H4 + C2^ 4.09E+12 4.15 1.49
5 C3H8+CH3 • — CH4 + n-C3H,^ 2.98E+08 29.93 3.65
6 C3H8+CH3 • — CH4 + iÄO-C3H7^ 5.48E+08 22.95 3.46
7 C3H8 +H — H2 + n-C3H7^ 2.55E+12 28.27 2.54
8 C3H8 +H — H2 + iÄO-C3H7^ 1.13E+12 18.71 2.4
9 n-C3H^ — C2H4 + CH3^ 1.2E+13 126 0
10 iro-C3H,^ — C3H + H 1.6E+13 150 0
11 n-C3H,^ — C3H + H 1.09E+13 149 0.17
12 C3H + n-C3H7^ 1.3E+13 13.64 0
13 C3H8 + C2^ — C2H6 + n-C3H,^ 9.7E+08 38.25 3.65
14 C3H8 + C2H3^ — C2H4 + iÄO-C3H7^ 4.79E+10 36.92 3.1
15 C3H + № — H2 +C3H5^ 2.61E+11 10.39 2.5
16 H2 +C3^ — C3H + H 8.37E+10 79.49 2.38
17 C3H8 + C3H5^ — C3H + n-C3H,^ 3.44E+10 83.06 3.3
18 C2H4+CH3 • — CH4 + C2H3^ 9.45E+09 39.74 3.7
19 CH4 + C2H3^ — C2H4+CH3 • 1.28E+10 22.86 4.02
20 H2+CH3 • — CH4 + H 1.52E+10 36.42 3.12
21 C2H4+ C2H3^ — C4H7 9.21E+11 19 0
22 C3H6 + C2H5^ — C2H + C3H5^ 1.02E+09 27.77 3.5
23 C4H7 —> C4H6 + H 6.4E+12 144.2 0
24 C4H7 — C2H4+ C2H3^ 2.1E+13 149.2 0
25 C2H + H^ — H2 +C2H5^ 2.48E+12 35.34 1.5
26 C3H5^ — C2H2 + C^ 3E+13 151 0
27 C4H8 —► C3^ + CH3^ 1E+16 305 0
28 CH3^ + CH3^ — C2H6 2.64E+13 0 0
29 C3H5^ + CH3^ —► C4H8 1.64E+13 -0.55 -0.32
30 C2H3^ + C^ — C3H 9.56E+14 0.57 -0.54
В процессе пиролиза принимают участие 17 веществ, из которых 10 — индивидуальные компоненты, 7 — радикалы, которые приведены в табл. 2. Для удобства принята нумерация всех участвующих в реакции компонентов.
Таблица 2
Индивидуальные вещества и радикалы, участвующие в процессе пиролиза
| Индивидуальные вещества | Радикалы |
1 -ЗД,2 - С2-Н4,3 - СН4,4 - Н2,5 - СЪНЬ,6 - С2Н6, П - с2Н5',12 - СН3',13 - И',
7 - С4Н7,8 - С4Н6,9 - С2Н2 ДО - С4Н8 14 - „ - ^,15 - 1 - С3Н7',
16 - С2Н3' ,17 - С3Н5'
На следующем этапе прописывалась система дифференциальных уравнений, которые учитывали изменение концентрации веществ во времени с учетом модифицированного уравнения Аррениуса:
- E
K = A • (—)П • eRT , 298
где А — предэкспоненциальный множитель, см3/(мольс); Е — энергия активации, кДж/моль; R — универсальная газовая постоянная, кДж/(мольК); Т — температура, К; n — температурная экспонента, экспериментально определяемый параметр.
В результате была получена следующая система дифференциальных уравнений:
dci
— = -K1 ■ c1 + K2 ■ c11 ■ c12 - K5 ■ c1 ■ c12 - K6 ■ c1 ■ c12 - K7 ■ c1 ■ c13 - K8 ' c1 ■ c13 - K13 ■ c1 ■ c11 - K14 ■ c1 ■ c16 - K17 c1 ■ c17 dc 2
— = K3 ■ У11 - K4 ■ c2 ■ c13 + K9 ■ c14 + K14 ■ C1 ■ c16 - K7 ' C1 ' c13 - K18 ' c2 ' c12 + K19 ' c3 ' c16 -K21 ' c2 ' c16 + K24 ' c7
dc3
-= K5 ■ c1 ' c12 + K6 ■ c1 ' c12 + K18 ' c2 ■ c12 - K19 ' c3 ' c16 + K20 ' c4 ■ c12
dt
DC4
-= Xy ■ ci ■ ci3 + K8 ■ c7 ■ ci3 + K15 ■ c5 ■ ci3 - K16 ■ c4 ■ c17 - K20 ' c4 ■ c12 + K25 ' c6 ' c13
dt
de 5
—— = K10 ■ e15 + K11 ■ e14 - K12 ■ e5 ■ e13 - K15 ■ e5 ' e13 + K16 'e4 ' e17 + K17 ' e1 ' e17 - K22 ' e5 ' e11 + K30 ' e16 dt
dc6
— = K13 ■ C1 ■ C11 + K22 ■ c5 ■ C11 - K25 ■ % ■ c13 + K28 ' c12 ■ c12
dt
dey
— = K21 ■ e2 ■ e16 - K23 ■ e7 - K24 ■ e7 ■ e12 dt
dc§
— = К23 ■ c7i dt
dcg
— = K26 ' c17
dt
dcio
—— = k29 ' c17 ' c12 - K27 ' c10
dt
den
-= ■ q - K2 ■ cn ■ ci2 - K3 ■ cn + К4 ■ с2 ■ c13 - K13 ■ q ■ cn - K21 ■ c5 ■ cn + K25 ' c6 ' c13
dt
dc12
—— = Ki ■ q - К2 • en ■ C12 - K5 • ci ■ C12 - K6 ■ c1 ■ c12 + K9 ■ c14 - K18 ' c2 ' c12 + K19 ' c3 ' c16
dt
— K20 ■ c4 ■ c12 + K26 ■ c17 + K27 ■ c10 - K28 ■ c12 ■ «12 - K29 ■ c12 ■ "17 - K30 ■ c12 ■ "16 dcl3
=K3 ■ c11-K4 'c2'C13-K7 ■ c1 ■ "13-K8 ■ c1 ■ c13+K10'c15 + К1Г C14-K12'c5 'C13-K15'c5 'c13+Кб c4 'c17+K20'c4 'c12 + K23'c7 -K25'c6 'c13
- K
^ ^^ т ' i " i iu W ü i-r W i.. 16 4 i/ 20 4 12 /
dci4
-= K5 • ñj.cji + K7 ■ ci ■ CJ3 -K9 ■ CJ4 -Kjj • C4 + Kj2 • C5 • ñj3 + Kj3 • cj • ñjj + Kj^ • cj • cj^
dt
de 15
—— = K6 ' ñl-e12 + K8 ' cl' e13 - K10 ' c15 + K14 ' el ' "6 dt
dc\6
—— = -Ki4 ■ ñ1.c16 + K18 • c2 • ci2 - K19 ■ c3 ' ñ16 - K21' c2 ' "16 + K24 ' c7 - K30 ' c12 ' "16 dt
dc\i
~~ = K15 ' "5-c13 - K16 ' c4 ' c17 - K17 ' cl' "17 + K22 ' c4 ' "11 - K26 ' c17 + K27 ' c10 - K29 ' c12 ' c17 dt
На основе полученной системы дифференциальных уравнений, описывающих процесс пиролиза, была написана программа в среде MATLAB, позволяющая рассчитать изменение концентрации компонентов во времени, а также константы скорости каждой реакции по формуле 2.
Целью данной работы являлось установление оптимальной температуры ведения процесса, при которой наблюдался максимальный выход этилена. Система дифференциальных уравнений решалась с помощью метода Гира.
По заводским данным рассматривались три варианта состава исходных реагентов:
1) ППФ с большим количеством пропилена — 75% и с содержанием пропана 25%;
2) ППФ с равным содержанием пропилена и пропана;
На основе разработанной математической модели вычислительный эксперимент проводили в температурном интервале 750—900 °С с шагом 50°С. Такой интервал и шаг варьирования температуры был принят на основе технологических возможностей работы промышленных реакторов.
Результаты и обсуждение
В табл. 3 приведены значения констант скорости реакций, полученные при расчете.
Таблица 3
Вычисленные значения констант скорости реакций
Стадия Значение, см3/с Стадия Значение, см3/с
1 2 3 4
1 1.2723e+010
2 9.4333e+008
3 1.3836e+013
4 1.0311e+014
5 1.8008e+016
6 9.1018e+015
7 8.1035e+016
8 1.3873e+016
9 1.1800e+010
10 1.5682e+010
11 3.3962e+010
12 1.2976e+010
13 5.8551e+016
14 6.8603e+016
15 6.3335e+015
16 8.8962e+014
17 1.9087e+017
18 8.0135e+017
19 9.5926e+018
20 2.4944e+016
21 9.1866e+008
22 2.2225e+016
23 6.2778e+009
24 2.0585e+010
25 6.6644e+013
26 2.9400e+010
27 9.6004e+012
28 2.6400e+010
29 1.8600e+009
30 2.4274e+010
На рис. 1—3 приведены зависимости выхода целевых и побочных продуктов от времени при проведения процесса пиролиза 100% пропана при различных температурах. Их анализ позволил сделать следующие выводы:
• Выход этилена при 850 °С составляет 18%масс, при 900 °С 38%масс. Повышение температуры проведения процесса позволяет повысить выход целевого продукта 20% масс.
• Выход пропилена при 850 °С составляет 13% масс., при 900°С 38%масс. Повышение температуры процесса позволяет повысить выход пропилена на 15% масс.
• Из полученных результатов видно, как при повышении температуры проведения процесса пиролиза выход метана уменьшается с 33% масс до 21%масс.
• По графику, представленному на рис. 3, видно, что пропан расходуется в первые 0.15 секунд, дальнейшее образование этилена происходит за счет свободно-радикальных реакций, протекающих в оставшееся время.
Рис. 1. График зависимости выхода продуктов при пиролизе пропана с течением времени при температуре 750°С.
Рис. 2. График зависимости выхода продуктов при пиролизе пропана с течением времени при температуре 850°С.
Рис. 3. График зависимости выхода продуктов при пиролизе пропана с течением времени при температуре 900°С.
Рис. 4. График зависимости выхода целевых продуктов от времени протекания процесса для ППФ фракции с содержанием пропана 25% об.
Рис. 5. График зависимости выхода целевых продуктов от времени протекания процесса для ППФ фракции с содержанием пропана 50% об.
При проведении вычислительного эксперимента с ППФ фракцией, содержащей 75% пропилена, были получены результаты, представленные на рис. 4. Для данного варианта процесса была так же подобрана оптимальная температура, составившая 850°С.
Затем был проведен процесс пиролиза с ППФ фракцией, в которой равное содержание пропана и пропилена. Результаты представлены на рис. 5. Для данного варианта процесса так же подбиралась температура, позволяющая достичь наибольшего выхода этилена, которая составила 850°С.
По представленным графикам можно сделать выводы:
• Оптимальный температурный режим, при котором протекает процесс пиролиза ППФ фракции любого состава — 850 °С. При данной температуре достигается наибольший выход целевого продукта — этилена.
• При проведении процесса с фракцией с большим содержанием пропилена, выход этилена составляет 12% масс. Повышение содержания пропана в ППФ на 25% масс. и снижение содержания пропилена на 25% масс. позволило повысить выход этилена при той же самой температуре процесса до 32.5% масс.
• Из графика, изображенного на рисунке 4, видно, что процесс протекает 0,5 секунд и останавливается из-за того, что пропан полностью израсходован. Исходя из этого, можно сделать вывод, что небольших концентраций пропана не хватает для развития свободно-радикального механизма.
Таким образом, данная компактная модель позволяет судить о протекании процесса, избегая громоздких вычислений и проведения сложных экспериментов, протекающих при условиях, которые трудно достичь в лаборатории. Было установлено, что проведение процесса пиролиза чистого пропана позволяет достичь высокого выхода этилена. Однако, использование в качестве сырья ППФ с одинаковым содержанием пропилена и пропана позволяет снизить выход метана на 17%, а также оптимальная температура, которая позволяет достичь высокий выход пропилена, ниже на 50°С. Это позволяет снизить экономические затраты, связанные с нагреванием при проведении процесса.
Работа выполнена в рамках реализации программы фундаментальных научных исследований государственных академий наук на 2013—2020 гг. (№ гос. Регистрации АААА-А17-117012610061-4) и частичной поддержки проекта РФФИ №1807-00341 А.
ЛИТЕРАТУРА
1. Хафизоф И. Ф. Современные тенденции развития процесса пиролиза / И.Ф. Хафизов, Р.Р. Мусин // Вестник Казанского технологического университета. 2015. С. 231-234.
2. Морозов А. Ю. Получение непредельных углеводородов в процессе пиролиза / А. Ю. Морозов, О. Н. Каратун, А. З. Саушин // Вестник Астраханского государственного технического университета. 2008. С. 161-163.
3. Беренц А. Д. Переработка жидких продуктов пиролиза / А.Д. Беренц [и др.]. М.: Химия, 1985. 216 с.
4. Приказ Минэнерго России от 01.03.2012 N 79 (с изм. от 03.12.2014) «О Плане развития газо- и нефтехимии России на период до 2030 года». СПС «КонсультантПлюс».
5. Snytnikov V. N. A reactor for the study of homogeneous processes using laser radiation energy / V.N. Snytnikov T. I. Mischenko, Vl. N. Snytnikov, I. G. Chernykh // Chemical Engineering Journal. 2009. Vol. 150. P. 231-236.
6. Снытников, В. Н. Автокаталитическое газофазное дегидрирование этана в «бесстеночном» реакторе / В. Н. Снытников, Т. И. Мищенко, Вл. Н. Снытников, О. П. Стояновская, В. Н. Пармон // Кинетика и катализ. 2010. Т. 51, №»1. C. 12-20.
7. Snytnikov V. N. Autocatalytic gas-phase dehydrogenation of ethane / V. N. Snytnikov, T. I. Mishchenko, Vl. N Snytnikov, S. E. Malykhin, V. I. Avdeev, V. N. Parmon // Research on Chemical Intermediates. 2012. Vol. 38. P. 1133-1147.
8. Snytnikov V. N. Autocatalytic dehydrogenation of propane / V. N. Snytnikov, T. I. Mishchenko, Vl. N. Snytnikov, I. G. Chernykh // Research on Chemical Intermediates. 2014. Vol. 40. P. 345-356.
9. Полищук Т. С., Череватюк Г. В., Патрушева О. В. Использование микроволнового излучения в нефтехимии // Молодой ученый. 2017. №2.1. С. 23-27.
10. Рахманкулов Д. Л., Шавшукова С. Ю., Даминев Р. Р. и др. Применение микроволнового излучения в нефтехимических процессах // Российский хим. журнал об-ва им. Д. И. Менделеева. 2008. №4. С. 136-141.
11. Нурисламова Л. Ф. Разработка компактной кинетической модели пиролиза пропана методами анализа чувствительности: дисс. канд. ... физ-мат. наук. ФГБОУН Институт нефтехимии и катализа РАН, Уфа, 2015.
Поступила в редакцию 01.11.2018 г.
SIMULATION ON THE PYROLYSIS OF PROPANE AND PROPANE-PROPYLENE FRACTION WITHIN MATLAB INTERFACE
© I. M. Gubaydullin1*, A. E. Musina2, R. V. Uzyanbaev3
1 Institute of Petrochemistry and Catalysis, Ufa Federal Research Center of RAS 141 Oktyabrya Avenue, 450075 Ufa, Republic of Bashkortostan, Russia.
2 Ufa State Petroleum Technological University 1 Kosmonavtov Street, 450044 Ufa, Republic of Bashkortostan, Russia.
3Ogarev Mordovia State University 68 Bolshevistskaya Street, 430005 Saransk, Russia.
Phone: +7 (347) 242 0712.
*Email: irekmars@mail.ru
In the present work, a computational methodology of the simulation on the pyrolysis of propane and propane-propylene fraction are described. It is based on the changes in concentrations of the chemicals in time and made up with the system of simple nonlinear differential equations. The solution of these equations enables constructing the whole picture of the process. Using the program written in the MATLAB program suite, the authors elucidated the optimal temperature regimes of the process and initial contents of the raw materials.
Keywords: pyrolysis, mathematical modelling, propane, propylene, compact model. Published in Russian. Do not hesitate to contact us at bulletin_bsu@mail.ru if you need translation of the article.
REFERENCES
1. Khafizof I. F. Vestnik Kazanskogo tekhnologicheskogo universiteta. 2015. Pp. 231-234.
2. Morozov A. Yu. Vestnik Astrakhanskogo gosudarstvennogo tekhnicheskogo universiteta. 2008. Pp. 161-163.
3. Berents A. D. Pererabotka zhidkikh produktov piroliza [Processing of liquid pyrolysis products] / A.D. Berents [i dr.]. Moscow: Khimi-ya, 1985.
4. Prikaz Minenergo Rossii ot 01.03.2012 N 79 (s izm. ot 03.12.2014) «O Plane razvitiya gazo- i neftekhimii Rossii na period do 2030 goda». SPS «Konsul'tantPlyus».
5. Snytnikov V N. Chemical Engineering Journal. 2009. Vol. 150. Pp. 231-236.
6. Snytnikov, V. N. Kinetika i kataliz. 2010. Vol. 51, No. 1. Pp. 12-20.
7. Snytnikov V N. Research on Chemical Intermediates. 2012. Vol. 38. Pp. 1133-1147.
8. Snytnikov V N. Research on Chemical Intermediates. 2014. Vol. 40. Pp. 345-356.
9. Polishchuk T. S., Cherevatyuk G V, Patrusheva O. V. Molodoi uchenyi. 2017. No. 2.1. Pp. 23-27.
10. Rakhmankulov D. L., Shavshukova S. Yu., Daminev R. R. i dr. Primenenie mikrovolnovogo izlucheniya v neftekhimicheskikh protsessakh. Rossiiskii khim. zhurnal ob-va im. D. I. Mendeleeva. 2008. No. 4. Pp. 136-141.
11. Nurislamova L. F. Razrabotka kompaktnoi kineticheskoi modeli piroliza propana metodami analiza chuvstvitel'nosti: diss. kand. ... fiz-mat. nauk. FGBOUN Institut neftekhimii i kataliza RAN, Ufa, 2015.
Received 01.11.2018.