Разработка компактной кинетической модели пиролиза пропана методами анализа чувствительности тема диссертации и автореферата по ВАК РФ 02.00.04, кандидат наук Нурисламова, Лиана Фануровна

• Нурисламова, Лиана Фануровна
• кандидат наук
• 2015

Оглавление диссертации кандидат наук Нурисламова, Лиана Фануровна

Оглавление

Введение

Глава 1. Литературный обзор

1.1 Актуальность проблемы моделирования пиролиза пропана

1.2 Исследование пиролиза пропана в лабораторном реакторе

1.3 Обзор кинетических моделей пиролиза пропана

1.4 Анализ чувствительности как метод получения компактной кинетической модели

1.5. Выводы по главе 1: задачи диссертационной работы

Глава 2. Методы анализа чувствительности математической модели реакции

2.1 Уравнения формальной химической кинетики

2.2 Общее уравнение динамики и кинетики реакции в режиме идеального вытеснения

2.3 Методы анализа чувствительности кинетических моделей

2.3.1 Локальный анализ чувствительности

2.3.2 Глобальный анализ чувствительности Соболя И.М

2.4. Анализ чувствительности кинетических кривых к изменению констант скоростей реакции на примере реакции гидроалюминирования олефинов алкилаланами

2.5. Выводы по главе 2

Глава 3. Методика анализа чувствительности функционала модели

3.1 Описание методики

3.2 Выбор вида функционала математической модели

3.3. Апробация методики

3.3.1 Анализ чувствительности математической модели реакции окисления

формальдегида

3.3.2 Анализ чувствительности математической модели реакции окисления водорода

3.3.3. Анализ чувствительности математической модели пиролиза этана

3.4. Выводы по главе 3

Глава 4. Компактная кинетическая модель пиролиза пропана: разработка и исследование

4.1 Экспериментальные данные по пиролиза пропана

4.2 Разработка компактной кинетической модели пиролиза пропана

4.2.1 Локальный анализ чувствительности функционала модели

4.2.2 Глобальный анализ чувствительности функционала

4.3 Исследование компактной кинетической модели пиролиза пропана

4.3.1 Влияние температуры и времени контакта на состав продуктов

4.3.2 Анализ чувствительности концентраций веществ к параметрам компактной кинетической модели

4.4 Расчет динамики химически реагирующего газа в реакторе

4.5 Выводы по главе 4

Заключение

Литература

Введение диссертации (часть автореферата) на тему «Разработка компактной кинетической модели пиролиза пропана методами анализа чувствительности»

Введение

Актуальность

В мировой нефтехимической промышленности наметилась тенденция к увеличению потребности в низших олефинах. В частности, Россия согласно плану развития газо- и нефтехимии на период до 2030 года, обладает потенциалом увеличения спроса почти в 4 раза к 2030 году по сравнению с 2010 годом. Рост потребления низших олефинов, получаемых пиролизом насыщенных углеводородов, влечёт за собой необходимость в создании принципиально новых технологий, отличающихся высокой производительностью и сниженными энергетическими затратами. Одним из перспективных направлений, позволяющим увеличить эффективность пиролиза является комбинирование термического и лазерного воздействия на алкан. Данный способ получения низших олефинов исследуется в Институте катализа имени Г.К. Борескова (группа к.ф.-м.н. Снытникова В.Н.), преимуществом которого является возможность снижения пороговой температуры и температуры выхода целевых продуктов реакции за счет дополнительной генерации радикалов в поле лазерного излучения. Такое снижение рабочей температуры заметно уменьшает образование кокса и побочных продуктов реакции. А инициация процесса лазерным излучением позволяет увеличить выход газообразных ненасыщенных углеводородов за счет повышения реакционной способности смеси. Чтобы понять эффективность подачи лазерного излучения, необходимо прежде детально проанализировать процесс, протекающий только за счет внешнего нагрева стенок реактора.

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

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

В настоящее время при анализе кинетических моделей широкое распространение получил метод анализа чувствительности численных моделей к вариации ее параметров (А.8. ТошНп, А. 8акеШ). Существующие подходы, как правило, основаны на анализе чувствительности температуры, концентраций всех либо одного из веществ, скорости стадий, скорости тепловыделения и пр. к вариации кинетических параметров модели. Однако при таком подходе интерпретация получаемых результатов анализа не проста: требуется анализировать матрицы коэффициентов чувствительности большой размерности, построенные для различных условий проведения процесса (в частности, при различных температурах и временах контакта). В связи с этим, актуальной научной задачей является разработка методики анализа сложных кинетических схем и редуцирования кинетических моделей до размеров, приемлемых с точки зрения точности описания и простоты практического использования. В настоящей работе метод анализа чувствительности применен для разработки компактной схемы химических

превращений при пиролизе пропана в условиях внешнего нагрева стенок реактора.

Цель работы заключается в разработке компактной кинетической модели, позволяющей описать процесс низкотемпературного (820-980 К) пиролиза пропана при атмосферном давлении.

Задачи исследования:

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

2. Разработка программы для проведения локального и глобального анализа чувствительности выходных параметров модели к ее входным параметрам.

3. Построение компактной кинетической модели газофазного пиролиза пропана с применением методики анализа чувствительности, адекватно описывающей экспериментальные данные низкотемпературного пиролиза пропана (820-980 К).

4. Исследование кинетики реакции газофазного пиролиза пропана при различных температурах проведения реакции и разных расходах смеси с использованием компактной кинетической модели.

5. Проведение численного моделирования трехмерной динамики газового потока в реакторе в программном пакете ANSYS Fluent с использованием компактной кинетической модели.

Научная новизна результатов исследования:

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

2. Построена кинетическая модель реакции газофазного пиролиза пропана, позволяющая описать процесс низкотемпературного (820-980 К) пиролиза пропана при атмосферном давлении, и приемлемая для проведения газодинамических расчетов в программных комплексах ЗО-моделирования.

3. Проведено комплексное исследование пиролиза пропана в неравновесных условиях проточного реактора: исследована кинетика реакции и проведены расчеты динамики газового потока пиролиза пропана в лабораторном реакторе с уютом процессов диффузии, химических реакций и их тепловых эффектов в программном пакете ANSYS Fluent с включением компактной схемы реакции.

Практическая значимость работы

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

Личный вклад автора

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

Апробация работы

Основные положения диссертационной работы докладывались и обсуждались на следующих Международных и Всероссийских конференциях:

~~ Международная суперкомпьютерная конференции «Научный сервис в сети Интернет» (Абрау-Дюрсо, 2013);

II Всероссийская научно-практическая конференция с международным участием "Математическое моделирование процессов и систем" (Стерлитамак, 2013);

Всероссийская научно-практическая конференция «Математическое моделирование на основе методов Монте-Карло» (Бирск, 2013); Международная конференция "Актуальные проблемы вычислительной и прикладной математики (Новосибирск, 2014);

XIV Международная конференция "Высокопроизводительные параллельные вычисления на кластерных системах" (Пермь, 2014). Связь с научными программами

Работа выполнена при поддержке грантов РФФИ № 12-07-31029 "Идентификация механизма реакции гидроалюминирования олефинов параллельными методами", № 12-07-00324 "Структурная и параметрическая идентификация кинетических моделей реакций нейтрального металлокомплексного катализа ", № 15-07-01764 "Оптимальное управление химическими реакциями металлокомплексного катализа", № 14-38-50639 мол нр «Исследование методами анализа нелинейных систем кинетической схемы газофазного пиролиза пропана». Публикации

По теме диссертации опубликовано 13 статей, из них 3 - в центральных научных журналах, входящих в перечень изданий, рекомендуемых ВАК РФ, 2 - в изданиях, индексируемые в Scopus. Получены 3 свидетельства о государственной регистрации программ для ЭВМ. Структура и объем работы

Материал диссертационной работы изложен на 118 страницах машинописного текста. Состоит из введения, 4 глав, заключения, списка литературы, включающего 124 наименования и содержит 14 таблиц и 36 рисунков.

Глава 1. Литературный обзор

1.1 Актуальность проблемы моделирования пиролиза пропана

В мировой нефтехимической промышленности наметилась тенденция к увеличению потребности в низших олефинах. В частности, Россия согласно плану развития газо- и нефтехимии на период до 2030 года, обладает потенциалом увеличения спроса почти в 4 раза к 2030 году по сравнению с 2010 годом [1]. Однако в настоящее время не хватает мощностей по производству базовых мономеров (этилен, пропилен, бутадиен), прежде всего, пиролизов. С другой стороны, в России существует достаточный объем сырья, который в будущем продолжит увеличиваться, и из-за отсутствия перерабатывающих мощностей большая часть сырья экспортируется. Уже сейчас спрос не удовлетворяется внутренним производством: доля импорта в российском потреблении основных видов пластиков в 2009 году составляла 10%, а по отдельным видам — около 30% (поливинилхлорид, полистирол и сополимеры стирола) [1].

В настоящее время в промышленных условиях пиролиз углеводородов осуществляют при температурах 800-900 °С и при давлениях, близких к атмосферному [2, 3]. Процессы термического пиролиза с "водяным паром" в трубчатых печах являются основными источниками получения этилена и пропилена, которые используются в качестве сырья для получения полиэтилена, полипропилена, фенола, ацетона, спиртов, лаков, растворителей в разных отраслях народного хозяйства, а также промежуточных продуктов для синтеза других веществ [2, 4]. Существующие технологии ориентированы на высокие производительности свыше 100 тыс тонн в год и более. В настоящее время стоит задача в задействовании низкодебитных скважин и создании технологии переработки углеводородов с небольшой до 50 тыс. тонн в год производительностью по сырью. Для таких производителыюстей существующие технологии

экономически неэффективны, поэтому требуются принципиально новые решения.

В Институте катализа имени Г.К. Борескова СО РАН рассматриваются способы создания экономически эффективной технологии переработки для невысоких производительностей путем управления химическими реакциями пиролиза через дополнительную генерацию радикалов в реагирующем газе посредством лазерного излучения, поглощаемого самой смесью реагентов [58].

Благодаря излучению С02-лазера в трубах змеевиков создаются локальные зоны высоких температур, которые служат дополнительным источником радикалов, что приводит к снижению пороговой температуры реакции и температуры выхода целевых продуктов приблизительно на 150°С (в пристеночной зоне). По данным экспериментов пиролиза этана [7], в таком процессе можно достичь конверсии этана порядка 85% при относительно невысоких температурах (600-700 °С) и атмосферном давлении. Такое снижение рабочей температуры весьма заметно уменьшает образование кокса, а также тяжелых продуктов пиролиза. За счет этого процесс можно вести без разбавления исходного этана водяным паром.

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

Преимуществом лазерной термохимии [10, 11] является возможность контролируемого вложения энергии таким образом, что излучение может стать управляющим элементом для генерации радикалов в объеме предварительно разогретого реакционного газа. Совмещение процесса

инициирования реакции с использованием лазерного излучения и получения высоко реакционно-способных соединений с последующим проведением каталитических реакций в условиях малых времен контакта и быстрой остановки протекания реакции может приводить к получению продуктов селективного превращения с высокими выходами и селективностью по целевым продуктам [6, 12].

1.2 Исследование пиролиза пропана в лабораторном реакторе

Экспериментальная установка по пиролизу пропана [5, 6] включала непрерывный С02-лазер, регистратор мощности лазерного излучения, металлический реактор, формирователь газовых потоков для задания состава исходной смеси, а также газовые хроматографы для анализа исходной газовой смеси и продуктов пиролиза. В проточном металлическом реакторе предусмотрен смешанный способ разогрева реагентов: внешний обогрев реакционной зоны и введение лазерного излучения непосредственно в газ. Схема экспериментальной установки (группа аэрозольного катализа ИК СО РАН) для исследования газофазного пиролиза С2-Сз углеводородов представлена на рис. 1.1.

Рис. 1.1. Схема установки с лазерным реактором для исследования газофазного пиролиза С2-С3 углеводородов. 1. С02-лазер; 2. Оптические окна (КС1, ЫаС1, 2п8е, ОаАв); 3. Измеритель мощности излучения — ЬМ-2; 4. Диффузное окно (ЫаС1); 5. Фурье-спектрометр «Инфралюм ФТ-801»; 6. Зеркала оптической трассы; 7. Ректор; 8. А1, 7п8е, ОаАв - зеркала реактора; 9. Измеритель мощности излучения - М154; 10. Хроматографы

«Кристалл 2000 М» и «ЛХМ-80 МД»; 11. Формирователь потока гаха УФГП-4,

Измеритель давления газа в реакторе

Источником излучения служил непрерывный СС^-лазер (1) с длиной волны >.=10.6 мкм, с максимальной мощностью до 130 Вт. Контроль мощности излучения производился с помощью светоделителыюй пластины из 7п8е (2) и измерителя мощности типа ЬМ-2 (3). Измерение спектральных характеристик генерации лазера осуществлялось с помощью Фурьеспектрометра «Инфралюм ФТ-801» (5). Основная часть излучения через систему зеркал (6) подавалась на вход реактора (7) под небольшим углом к их продольной оси. Входными окнами служили или просветленные пластины из 2п8е или кристаллы КС1. На выходе реакторов устанавливались плоские зеркала из А1, 2п8е или Азва (8), непоглощенное лазерное излучение в реакторе регистрировалось с помощью измерителя мощности М154 (9). Химический анализ состава исходной смеси газов и продуктов реакции производили на хроматографах «Кристалл 2000 М» и «ЛХМ-80 МД» (10). Смешение исходных газов и задание их расхода через реакторы производилось с помощью 4-х канального формирователя потоков типа УФПГ-4 (11). Для измерения температуры в пристеночной области реактора использовались хромель-алюмелевые термопары, а в объеме реактора и в поле излучения лазера была использована термопара с электродами из молибдена и меди.

Схема реактора приведена на рис. 1.2. Реактор (1) изготовлен из нержавеющей стали 12X18Н9Т с внутренним диаметром 21 мм, длиной 190 мм и толщиной стенок 2 мм. Характерная длина реакционной зоны составляла 50 мм. С внешней стороны реакционной зоны располагался теплоизолированный нагреватель (6) с максимальной тепловой мощностью 350 Вт. Оптические окна лазерного реактора (2) выполнены из просветленного в области 10 мкм 2п8е. Два ввода (3) служили для подачи в реактор пропановой или пропан-этиленовой смеси. Выход смеси защитного газа и продуктов пиролиза производится через канал (5). Вводы (4)

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

—' JJ | VL ±J

Рис. 1.2. Реактор пиролиза с внешним нагревом. 1. Реактор, диаметр—21 мм, длина - 190 мм ; 2. Окна, Si02; 3. Ввод исходной смеси газов, диаметр - 4 мм; 4. Входы защитного газа СНЦ, диаметр - 4 мм; 5. Выход продуктов пиролиза, диаметр - 6мм; 6.

Подогреватели; 7. Защитный конус из полированной меди; 8. Реакционная зона.

1.3 Обзор кинетических моделей пиролиза пропана

Моделирование химического процесса является необходимым инструментом как на этапе проектирования реактора (особенно для конструкций сложной геометрии) с целью изучения процесса в широком спектре различных условий протекания реакции и технологических режимов и для оптимального управления процессом. Решающее значение при выборе условий проведения химико-технологических процессов имеют вопросы определения скорости химических превращений, изучаемые химической кинетикой. Кинетические уравнения, содержащие необходимую информацию об основных закономерностях химических превращений, являются первоосновой математической модели химического процесса и реактора [13]. Количество стадий и веществ, включенных в кинетическую модель, должно быть ограничено. Это мотивируется необходимостью интегрирования кинетической модели в программный комплекс 3D-моделирования FLUENT для расчета трехмерной динамики газа. Брутто-схемы не дают удовлетворительного описания пространственного

распределения по реактору протекания химических реакций с радикально-цепными механизмами. Использование детальных кинетических механизмов при моделировании процесса пиролиза пропана важно для полного и точного описания процессов в широком диапазоне температур, однако их использование в трехмерных численных расчетах реакторов приводит к неприемлемым временам расчетов даже на современных суперкомпьютерах. Поэтому при CFD-моделировании (computational fluid dynamics), как правило, используют компактные схемы, что обусловлено большим временем счета по детализированным моделям [14, 15]. Таким образом, необходимо найти баланс между точностью описания модели и временем вычислений [16-18].

Термическое разложение предельных углеводородов сопровождается разнообразными и многочисленными параллельными процессами, что приводит к образованию десятка веществ и возможности протекания сотни реакций. Состав продуктов пиролиза весьма многообразен и может варьироваться в широких пределах в зависимости от типа сырья и технологических условий проведения реакции. Именно поэтому исследователями предложены многочисленные кинетические схемы для описания процессов пиролиза углеводородов на различных лабораторных и промышленных установках и накоплен большой объем экспериментальных данных [2, 19, 20].

Первые кинетические модели пиролиза пропана строились вручную и основывались на теории цепных реакций [21, 22]. Пиролиз пропана является радикально-цепным процессом и протекает согласно механизму Райса-Герфельда -Косякова, который складывается из трех стадий: инициирование цепи, продолжение цепи и обрыв цепи [2, 21]. Общепринятый механизм пиролиза пропана включает следующие основные реакции [4]:

Стадия инициирования цепи

С3Н8->С2Н5- + СН3-Стадии продолжения цепи

СН3« + С3н8 -> СН4+ п-СъН7« СН3« + С3Н8^ СН4 + ш?-С3Н7-п-С3Н7« С2Н4 + СН3» 150-С3Н7- С3Н6 + Н« Н« + С3Н8 Н2+ п-С3Н7« Н« + С3Н8 -> Н2 + ¿?о-С3Н7 С2Н5- -> С2Н4 + н-С2н5. + с3н8 С2Н6 + л-С3Н7. С2Н5- + С3Н8 -> С2Н6 + ио-С3Н7-Н« + «о-С3Н7« -> н2 + с3н6

Стадии обрыва цепи

н« + н« н2 с2н5- + С2н5. -> С2Н4 + С2Н6 н- + С2Н5- -> С2Н6 СНз« + СНз* С2Н6 СНз' + С2Н5« -> СН4 + С2Н4 В работе Мухиной Т.Н. [2] представлена аналогичная схема с небольшими изменениями и содержит следующие стадии обрыва цепи:

Н- + и-С3Н7« -»• Н2 + С3Н6 Н« + я-С3Н7« С3Н8 СНз- + ¿уо-С3Н7« -»• СН4 + С3Н6 СН3« + мо-С3Н7- -> С4Ню Для выявления наиболее существенных реакций обрыва цепи часто проводят предварительный кинетический анализ.

Анализ кинетической модели пиролиза пропана по этим двум схемам с относительно небольшим числом стадий показал, что расчетные кинетические зависимости расходятся с экспериментальными данными для низкотемпературного пиролиза при атмосферном давлении и высоких конверсиях пропана (рис. 1.3). Кинетические расчеты пиролиза пропана по [2, 4] для температуры 970 К показали, что основными продуктами пиролиза

пропана являются С3Н6 и Н2, количество которых нарабатывалось в равных долях. Другая пара продуктов СН4 и С2Н4 нарабатывалась в меньших количествах, но также в равных долях.

Длина реакционной зоны, мм Рис. 1.3. Временные зависимости концентраций пропана и основных продуктов пиролиза, полученных на основе схемы, представленной в [4]

В работе [8] была предложена автокаталитическая схема пиролиза пропана в смеси с этиленом при атмосферном давлении для температуры 1000 К. Согласно этой схеме ключевой реакцией зарождения радикалов является стадия взаимодействия пропана и этилена с получением н-пропильного и этильного радикалов по двухстадийному маршруту: образование бирадикала этилена и его взаимодействие с пропаном. Кинетические расчеты температуры в районе 970 К показали удовлетворительное согласие с экспериментальными данными. Однако при более высоких и низких температурах расчеты по представленной в [8] схеме показали расхождение с экспериментальными данными.

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

Примерами использования схем дегидрирования пропана с большим набором элементарных стадий могут служить работы [23-25], которые представляют результаты численных исследований пиролиза пропана в адиабатическом реакторе и реакторе быстрого смешения для температур 1000 и 1250 К. Предложенный в работе [26] механизм состоит из 827 элементарных реакций и 227 веществ.

Практически используется в настоящее время модель «ЗРУЯО» нидерландской фирмы КТ1 и модель «Терасуг» ВНИИОС (Россия). Модель «Терасуг» включает около 650 реакций и 77 молекулярных веществ, 72 радикала. Программа «БРУЯО» позволяет моделировать распад углеводородов в пирозмеевиках. Полная кинетическая модель этой программы включает 2000 реакций радикально-цепного механизма термического распада углеводородов.

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

А.Б. ТошНп и др. [27] представили анализ кинетической схемы высокотемпературного пиролиза пропана [28], включающей 422 химические реакции, и показали, что уменьшение числа элементарных стадий вплоть до 50 наиболее важных, дает схожие численные кинетические зависимости по продуктам пиролиза.

Были рассмотрены различные детальные схемы химических превращений, предназначенные для описания пиролиза углеводородов с различными типами сырья и технологических условий проведения реакции:

схемы, приведенные в работах Raseev S.D. [29], Мухиной Т.Н. [2], Saeys М. [30], Marinov N.M. [31].

В работе Raseev S.D. [29] приведены некоторые наиболее важные стадии из модели «SPYRO», в работе Мухиной Т.Н. [2] - модели «Терасуг». Sayes М. исследует процесс промышленного парового крекинга углеводородов и в его работе [30] приведена детальная схема процесса. Наибольшее применение среди исследователей получила работа [31]. Предложенный механизм был разработан для моделирования ароматических и полициклических ароматических углеводородов, образованных в пламени предварительно перемешанной бутан-аргон-кислородной смеси при атмосферном давлении.

В основу расчетов была положена детальная схема, полученная на основе анализа и комбинирования схем, приведенных в работах Raseev S.D., Saeys М., Marinov N.M, Мухиной Т.Н. [2, 29, 30, 31]. Расход пропана, образование основных (СН4, С2Н4, Н2, С3Н6)и побочных (С2Н6, С4Н6, С4Н8, С2Н2) продуктов реакции наилучшим образом соответствовало экспериментальным данным с использованием кинетических параметров, приведенных в работе Marinova N.M. (рис. 1.4).

Полученная схема состоит 157 стадий и включает в себя следующие вещества, наблюдаемые в ходе реакции: С3Н8, СН4, С2Н4, Н2, С3Н6, С2Н6, С4Н6, С4Н8, С2Н2; а также С4Ню, /¿о-С^Ню и радикалы (Н% С2Н5% СН3», С2Н3% С3Н5*, h-C3N7«, «о-С3Н7% С4Н7*, n-C4H9», iso-C4H9*). Углеводороды выше, чем С4 не были включены, поскольку при низкотемпературном пиролизе пропана появление С5 и С6 углеводородов, если и возможно, то в очень малых количествах. А их включение значительно бы увеличило число стадий реакции.

Список литературы диссертационного исследования кандидат наук Нурисламова, Лиана Фануровна, 2015 год

Литература

1. Приказ Минэнерго России от 01.03.2012 N 79 (с изм. от 03.12.2014) «О Плане развития газо- и нефтехимии России на период до 2030 года». -СПС «КонсультантПлюс».

2. Мухина, Т.Н. Пиролиз углеводородного сырья / Т.Н.Мухина H.J1. Барабанов, С.Е. Бабаш и др. - Москва: Химия, 1987. - 240 с.

3. Богомолов, А.И. Химия нефти и газа / А.И. Богомолов; под ред. В.А. Проскурякова. - 3-е изд., доп. и исправл. - СПб: Химия, 1995. - 446 с.

4. Жоров, Ю.М. Кинетика промышленных органических реакций: справочник / Ю. М. Жоров. - Москва : Химия, 1989. - 384 с.

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. Снытников, B.H. Автокаталитическое газофазное дегидрирование этана в «бесстеночном» реакторе / В.Н. Снытников, Т.П. Мищенко, Вл.Н. Снытников, О.П. Стояновская, В.Н. Пармон // Кинетика и катализ. - 2010. -Т. 51,№ 1.-С. 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, Т. I. Mishchenko, VI. N. Snytnikov, I. G. Chernykh // Research on Chemical Intermediates. - 2014. - Vol. 40. - P. 345-356.

9. Козлов, Г.И. Сильное воздействие мощного С02 лазера на горение пропано-воздушных смесей / Г.И. Козлов, В.А. Кузнецов, А.Д. Сокуренко // Письма в ЖТФ. - 1990. - Т. 18, № 9. - С. 55-59.

Ю.Баранов, В.Ю. Воздействие лазерного излучения с Я=308 нм на пиролиз 1,2-дихлорэтана / В.Ю.Баранов В.М.Борисов, А.В.Виноходов,

Ф.И.Высикайло, Ю.Б.Кирюхн, В.Д.Письменный, О.А.Зайдман, Ю.А.Трегер, К.А.Чагир, И.С.Заслонко, Ю.К.Мукосеев, В.Н.Смирнов // Квантовая электроника. - 1983. - Т. 10, № 7. - С. 1406-1412. 11 .Манташян, А.А. Сложные газофазные реакции превращения углеводород-кислородных смесей в условиях лазерного нагрева / А.А. Манташян, JI.A. Хачатрян, А.А. Сукиасян // Кинетика и катализ. - 1989. -Т. 2, № 30. - С. 272-276.

12.Сукиасян, А.А. Окислительный крекинг пропана под воздействием лазерного излучения / А.А. Сукиасян, JI.A. Хачатрян, А.А. Манташян. // Кинетика и катализ. - 1989. - Т. 2, № 30. - С. 277-280.

13.Слинько, М.Г. Основы и принципы математического моделирования каталитических процессов / М.Г. Слинько. - Новосибирск : Наука, 2004. -488 с.

14.Gascoina, G.N. Methane pyrolysis: literature survey and comparisons of available data for use in numerical simulations / G.N. Gascoina, Ph. Gillarda, J. Steelant // Journal of Analytical and Applied Pyrolysis. - 2013. - Vol. 104. -P. 1-9.

15.Bedarev, I. Numerical study of methane pyrolysis in shock waves / I. Bedarev, V. Parmon, A. Fedorov, N. Fedorova, V. Fomin // Combustion, Explosion, and Shock Waves. - 2004. - Vol. 40. - P. 580-590.

16.Кондратьев, B.H. Кинетика и механизм газофазных реакций / В.Н. Кондратьев, Е.Е. Никитин. - Москва : Наука, 1974. - 558 с.

17.Варнатц, Ю. Горение. Физические и химические аспекты, моделирование, эксперименты, образование загрязняющих веществ / Ю. Варнатц, У. Маас, Р. Диббл. - Москва : Физматлит, 2003. - 352 с.

18.Губайдуллин, И.М. Редуктивный подход при моделировании сложных задач химической кинетики / И. М. Губайдуллин, В. Б. Маничев, JT. Ф. Нурисламова // Журнал Средневолжского математического общества. — 2012.-Т. 4, №4. _с. 26-33.

19. Khana, R.U. Pyrolysis of propane under vacuum carburizing conditions: an experimental and modeling study / R.U. Khana, S. Bajohra, et.al. // Journal of Analytical and Applied Pyrolysis. - 2008. - Vol. 81, № 2. - P. 148-156.

20.Westbrook, C.K. A comprehensive chemical kinetic reaction mechanism for oxidation and pyrolysis of propane and propene / C.K. Westbrook, W.J. Pitz // Combustion Science and Technology. - 1984. - Vol. 37. - P. 117-152.

21.Эмануэль, H.M. Курс химической кинетики: учебник для химических факультетов университетов / Н.М. Эмануэль, Д.Г. Кнорре. - 4-е изд., перераб. и доп. - Москва : Высшая школа, 1984. - 463 с.

22.Панченков, Г.М. Химическая кинетика и катализ: учебное пособие для вузов / Г.М. Панченков, В.ГТ. Лебедев. - 3-е изд., исправл. и доп. — Москва : Химия, 1985. - 592 с.

23.Starik, A.M. Kinetics of the oxidation of the products from the thermal destruction of C3H8 and п-СДТю in the mixtures with air / A.M. Starik, N.S. Titova, L.S. Yanovskii // Kinetics and catalysis. - 1999. - Vol. 40, № 1. - P. 722.

24.Ktalkherman, M. G. Investigation of high-temperature pyrolysis of propane in a fast-mixing reactor / M. G. Ktalkherman, I. G. Namyatov, V. A. Emel'kin, B. A. Pozdnyakov // High Temperature. - 2009. - Vol. 47, № 5. - P. 707-717.

25.Кталхерман, М.Г. Пиролиз углеводородов в потоке теплоносителя при быстром смешении компонентов / М.Г. Кталхерман, И.Г. Намятов // Физика горения и взрыва. - 2008. - Т. 44, № 5. - С. 38-44.

26.Norinaga, К. Detailed kinetic modeling of gas-phase reactions in the chemical vapor deposition of carbon from light hydrocarbons / K. Norinaga O. Deutschmann // Ind. Eng. Chem. Res. - 2007. - T. 46, № 11. - P. 3547-3557.

27.Tomlin, A.S. Reduced mechanisms for propane pyrolysis / A.S. Tomlin, M.J. Pilling, J.H. Merkin, J. Brindley, N. Burgess, A. Gough // Ind. Eng. Chem. Res. - 1995. - T. 34. - P. 3749-3760.

28.Dente, M. E., Ranzi, E. M. Mathematical modeling of hydrocarbon pyrolysis reactions // Pyrolysis: Theory and industrial practice В кн.: Pyrolysis: Theory

and industrial practice / L. F. Albright, B. L. Crynes, W. H. Corcoran. - New York, USA : Academic Press, 1983. - P. 133-175. 29.Raseev, S.D. Thermal and catalytic processes in petroleum refining / S.D.

Raseev. - New York, USA, 2003. - 920 p. 30.Saeys, W.S. Construction of an ab initio kinetic model for industrial ethane pyrolysis / W.S. Saeys Wenjie // AIChE Journal. - 2011. - Vol. 57, № 9. - P. 2458-2471.

31.Marinov, N.M. Aromatic and polycyclic aromatic hydrocarbon formation in a laminar premixed n-butane flame / N.M. Marinov, W.J. Pitz, C.K. Westbrook, et.al. // Combustion and Flame. - 1998. - Vol. 114, № 1-2. - P. 192-213.

32.Database NIST Chemical Kinetics [Электронный ресурс]. - Режим доступа: http://kinetics.nist.gov/kinetics/.

33.Tsang, W. Chemical kinetic data base for combustion chemistry. Part I. Methane and related compounds / W. Tsang, R.F. Hampson // J. Phys. Chem. Ref. Data. - 1986. -Vol. 15.-P. 1087-1279.

34.Tsang, W. Chemical kinetic data base for combustion chemistry. Part 3. Propane / W. Tsang // J. Phys. Chem. Ref. Data. - 1988. - Vol. 17. - P. 887951.

35.Baulch, D.L. Evaluated kinetic data for combusion modelling. Supplement I. / D.L. Baulch, C.J. Cobos, et. al. // J. Phys. Chem. Ref. Data. - 1994. - Vol. 23. -P. 847-1033.

36.Tsang, W. Chemical kinetic data base for combustion chemistry. Part 3. Propane / W. Tsang // J. Phys. Chem. Ref. Data. - 1988. - Vol. 17. - P. 887951.

37.Harding, L.B. Predictive theory for hydrogen atom - hydrocarbon radical association kinetics / L.B. Harding, Y. Georgievskii, S.J. Klippenstein // J. Phys. Chem. A. - 2005. - Vol. 109. - P. 4646 - 4656.

38.Stewart, P.H. Pressure and temperature dependence of reactions proceeding via a bound complex. 2. Application to 2CH3 —► C2H5 + H / P.H. Stewart, C.W.

Larson, D.M. Golden 11 Combustion and Flame. - 1989. - Vol. 75. - P. 25 -31.

39.Dean, A.M. Predictions of pressure and temperature effects upon radical addition and recombination reactions / A.M. Dean // J. Phys. Chem. - 1985. -Vol. 89. - P. 4600 - 4608.

40.Harding, L.B. On the combination reactions of hydrogen atoms with resonance-stabilized hydrocarbon radicals / L.B. Harding, S.J. Klippenstein, Y. Georgievskii // J. Phys. Chem. A. - 2007. - Vol. 111. - P. 3789 - 3801.

41.Tsang, W. Chemical kinetic data base for combustion chemistry. Part V. Propene / W. Tsang // J. Phys. Chem. Ref. Data. - 1991. - Vol. 20. - P. 221 -273.

42.Rao, V.S. Formation of H and D atoms in pyrolysis of 1,3-butadiene and 1,3 butadiene-1,1,4,4,-d4 behind shock waves / V.S. Rao, K. Takeda, G.B. Skinner // Int. J. Chem. Kinet. - 1988. - Vol. 20. - P. 153-164.

43.Zhang N. CFD simulation of propane cracking tube using detailed radical kinetic mechanism / N. Zhang, T. Qiu, B. Chen // Chinese Journal of Chemical Engineering.-2013.-Vol. 21, № 12.-P. 1319-1331.

44.Warnatz, J. Rate coefficients in the C/H/O system. B kh.: Combustion Chemistry / W.C. Gardiner. - New York : Springer, 1984. - P. 197-360.

45.Baldwin, R.F. Rate constants for hydrogen + oxygen system, and for H atoms and OH radicals + alkanes / R.F. Baldwin, R.W. Walker // J. Chem. Soc. Faraday Trans. 1. - 1979. - Vol. 75.-P. 146-154.

46.Tsang, W. Chemical kinetic data base for combustion chemistry. Part 4. Isobutane / W. Tsang // J. Phys. Chem. Ref. Data. - 1990. - Vol. 19. - P. 1-68.

47.Curran, H.J. Rate constant estimation for C-l to C-4 alkyl and alkoxyl radical decomposition / H.J. Curran // Int. J. Chem. Kinet. - 2006. - Vol. 38. - P. 250275.

48.Kungwan, N. Kinetics of the hydrogen abstraction center dot CH3 plus alkane -> CH4 plus alkyl reaction class: an application of the reaction class transition

state theory / N. Kungwan, T.N. Truong // J. Phys. Chem. A. - 2005. - Vol. 109.-P. 7742-7750.

49.Lin, M.C. The thermal decomposition of ethane. Part III. Secondary reactions / M.C. Lin, M.H. Back // Can. J. Chem. - 1966. - Vol. 44, № 20. - P. 23692380.

50.Matheu, D.M. Capturing pressure-dependence in automated mechanism generation: reactions through cycloalkyl intermediates / D.M. Matheu, W.H. Green, J.M. Grenda // Int. J. Chem. Kinet. - 2006. - Vol. 35. - P. 95 - 119.

51.Lovas T. Model reduction techniques for chemical mechanisms. В кн.: Chemical Kinetics / V. Patel. - InTech, 2012. - P. 79-114.

52.Розенвассер, E.H. Чувствительность систем автоматического управления / Е.Н. Розенвассер, P.M. Юсупов. - Москва: Наука, 1981. -464 с.

53.Brown N. J. Mechanism reduction via principal component analysis / N. J. Brown, G. Li, M. L. Koszykowski // Int. J. Chem. Kinet. - 1997. - Vol. 29, № 6.-P. 393-414.

54.Pepiot, P. Systematic Reduction of Large Chemical Mechanisms / P. Pepiot, H. Pitsch // Proc. 4th Joint Meeting of the U.S. Sections of the Combustion Institute. - Philadelphia, USA, 2005. - 6 p.

55.Miles, S.O. Simplification of mathematical models of chemical reaction systems / S.O. Miles, M.L. Mavrovouniotis // Chemicals Reviews. - 1998. -Vol. 98, №2.-P. 391-408.

56.Shi, Y. Automatic chemistry mechanism reduction of hydrocarbon fuels for HCCI engines based on DRGEP and PCA methods with error control / Y. Shi, H.-W. Ge , J. L. Brakora, R.D. Reitz // Energy Fuels. - 2010. - Vol. 24. - P. 1646-1654.

57.Tomlin, A.S. Mathematical tools for the construction, investigation and reduction of combustion mechanisms / A.S. Tomlin, T. Turanyi, M.J. Pilling // Comprehensive chemical kinetics. - 1997. - Vol. 35. - P. 293-437.

58.Lu, T. F. A directed relation graph method for mechanism reduction / T.F. Lu, C.K. Law // P. Combust. Inst. - 2005. - Vol. 30. - P. 1333-1341.

59.Lu, Т. F. Linear time reduction of large kinetic mechanisms with directed relation graph: N-heptane and isooctane / T.F. Lu, C.K. Law // Combust Flame. - 2006. - Vol. 144. - P. 24-36.

60.Tomlin, S. Mechanism reduction for the oscillatory oxidation of hydrogen: sensitivity and quasi-steady-state analyse / S. Tomlin, M. Pillinc, T. Turanyi, J. Merkin, J. Brindley//Combustion and Flame. - 1992,-Vol. 91.-P. 107-130.

61.T. Turanyi Sensitivity analysis of complex kinetic systems. Tools and applications // Journal of Mathematical Chemistry. - 1990. - Vol. 5, 3. - P. 203-248.

62-Saltelli, A. Sensitivity analysis practices: strategies for model-based inference / A. Saltelli, M. Ratto, S. Tarantola, F. Campolongo // Reliability Engineering & System Safety. - 2006. - Vol. 91, № 10-11. - P. 1109-1125.

63.Применение вычислительной математики в химической и физической кинетике / под ред. Л.С.Полака. - М.: Наука, 1969. - 280 с.

64.Полак, JI.C. Вычислительные методы в химической кинетике / JI.C. Полак, М.Я. Гольденберг, А. А. Левицкий. - Москва : Наука, 1985. - 280 с.

65.Saltelli, A. Sensitivity Analysis for Chemical Models / A. Saltelli, M. Ratto, S. Tarantola, F. Campolongo // Chem Rev. - 2005. - Vol. 205, № 7. - P. 28112828.

66.Xia, A. G. Mechanism reduction for the formation of secondary organic aerosol for integration into a 3-dimensional regional air quality model: a-pinene oxidation system / A. G. Xia, D. V. Michelangeli // Atmos. Chem. Phys. - 2009. - Vol. 9. - P. 4341-4362.

67.Марчук, Г.И. Методы вычислительной математики / Г.И. Марчук — Москва : Наука, 1980. - 536 с.

68.Волкова, Е.В. Методы анализа чувствительности для моделей фильтрации и массопереноса в подземной гидросфере: дис. ... канд. тех. наук: 05.13.18 / Е.В. Волкова. - Казань, 2009. - 144 с.

69.Кокотовича, П. В. Чувствительность систем автоматического управления / П. В. Кокотовича, Р. С. Рутман // Автомат, и телемех. - 1965. - Т. 6, № 4. -С. 730-750.

70.Соболь, И.М. Об оценке чувствительности нелинейных математических моделей // Математическое моделирование. — 1990. — Т. 2, № 1. - С. 112118.

71.Соболь, И.М. Анализ чувствительности нелинейных математических моделей: численные опыты / И.М. Соболь, А. Сальтелли // Математическое Моделирование. - 1995. - Т. 7, № 11. - С. 16-28.

72.Соболь, И.М. Глобальные показатели чувствительности для изучения нелинейных математических моделей / И.М. Соболь // Математическое моделирование. - 2005. - Т. 7, № 9. - С. 43-52.

73.Соболь, И.М. Глобальные показатели чувствительности для изучения нелинейных математических моделей(дополнение) / И.М. Соболь // Математическое моделирование. - 2007. - Т. 19, № 11. - С. 23-24.

74.Крянев, A.B. Математические методы обработки неопределенных данных / A.B. Крянев, Г.В. Лукин. - Москва : Физматлит, 2006. - 216 с.

75.Быков, В.И. Моделирование критических явлений в химической кинетике / В.И. Быков. - Москва : КомКнига, 2006. - 328 с.

76.Прямые и обратные задачи в химической кинетике / Под ред. В.И. Быкова. - Новосибирск: Наука, 1993. - 288 с.

77.Брин, Э.Ф. Обратные задачи химической кинетики как метод исследования механизмов сложных реакций / Э.Ф. Брин // Успехи химии. - 1987.-№3.-С. 428-446.

78.Яблонский, Г.С. Кинетические модели каталитических реакций / Г.С. Яблонский, В.И. Быков, А.Н. Горбань. - Новосибирск : Наука, 1983. — 254 с.

79.Димитров, В.М. Простая кинетика / В.М. Димитров. - Новосибирск : Наука, 1982. - 379 с.

80.Яблонский, Г.С. Математические модели химической кинетики / Г.С. Яблонский, С.И. Спивак. - Москва : Знание, 1977. - 64 с.

81.Штиллер, В. Уравнение Аррениуса и неравновесная кинетика / В. Штиллер. - Москва : Мир, 2000. - 176 с.

82.Семиохин, И.А. Физическая химия / И.А. Семиохин. - Москва : МГУ, 2001.-272 с.

83.Байназарова Н.М. Анализ чувствительности кинетических кривых к изменению констант скоростей реакции модели реакции гидроалюминирования олефинов / Н. М. Байназарова, JI. Ф. Нурисламова, И. М. Губайдуллин // Журнал Средневолжского математического общества. - 2013. - Т. 15, № 1. - С. 34-40.

84.Байназарова, Н. М. Анализ чувствительности кинетических параметров модели реакции гидроалюминирования олефинов с диизобутилалюминийхлоридом / Н. М. Байназарова, К. Ф. Коледина, Л. Ф. Нурисламова // Научный сервис в сети Интернет: все грани параллелизма: труды международной суперкомпьютерной конференции. - Новороссийск, 2013. - С. 45-50.

85.Zador, J. Local and global uncertainty analysis of complex chemical kinetic systems / J. Zador, I.G. Zsely, T. Turanyi // Rel Eng & Sys Safety. - 2006. -Vol. 91, № 10-11.-P. 1232-1240.

86.Sobol' I.M. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates / I.M. Sobol' // Mathematics and Computers in Simulation.-2001.-Vol. 55, № 1-3.-P. 271-280.

87.Канторович, Л.В. О некоторых новых подходах к вычислительным методам и обработке наблюдений / Л.В. Канторович // Сибирский математический журнал. - 1962. - Т. 3, № 5. - С. 701-709.

8 8. Соболь, И.М. Равномерно распределенные последовательности с дополнительным свойством равномерности / И.М. Соболь // Ж. вычисл. матем. и матем. физ. - 1976. - Т. 16, № 5. - С. 1332-1337.

89.Губайдуллин, И.М. Методология распараллеливания при решении многопараметрических задач химической кинетики / И.М. Губайдуллин, Ю.Б. Линд, К.Ф. Коледина //Вычислительные методы и программирование. - 2012. - Т. 13, №1. - С. 236-244.

90.A\vrejcewicz, J. Modern information technologies in construction of kinetic models for reactions of metal complex catalysis / J. Awrejcewicz, K.F. Koledina, Y.B. Lind, I.M. Gubaidullin // Theoretical and Applied Mechanics Letters. - 2012. - Vol. 2, №4.-4 p.

91.Нурисламова, Л. Ф. Двухкритериальная идентификация кинетических параметров реакции гидроалюминирования олефинов алкилаланами / И. М. Губайдуллин, А .П. Карпенко, Л. Ф. Нурисламова, А. С. Савелов // Наука и образование: электронное научно- техническое издание. — 2013. — № 12. URL: http://technomag.bmstu.ru/doc/645511.html (дата обращения: 12.12.2013).

92.Коледина, К.Ф. Последовательно-параллельное определение кинетических параметров при моделировании детального механизма гидроалюминирования олефинов: дне. ... канд. физ.-мат. наук: 02.00.04 / К.Ф. Коледина. - Уфа, 2011. - 111 с .

93.Nurislamova, L. F. Kinetic Model of Isolated Reactions of the Catalytic Hydroalumination of Olefins / I. M. Gubaydullin, K.F. Koledina, L. F. Nurislamova // Reaction Kinetics, Mechanisms and Catalysis. - 2015. DOI: 10.1007/s 11144-015-0876-6.

94.Boivin, P. An explicit reduced mechanism for H2-air combustion / P. Boivin, C. Jimenez, A.L. Sanchez, F.A. Williams // Proceedings of the Combustion Institute. - 2011. - Vol. 33, № 1. - P. 517-523.

95.Нурисламова, Л. Ф. Двухкритериальная идентификация кинетических параметров реакции гидроалюминирования олефинов алкилаланами / И. М. Губайдуллин, А .П. Карпенко, Л. Ф. Нурисламова, А. С. Савелов // Наука и образование: электронное научно- техническое издание. — 2013. —

№ 12. URL: http://tcchnomag.bmstu.ru/doc/645511.html (дата обращения: 12.12.2013).

96.Turanyi, Т. Applications of sensitivity analysis to combustion chemistry / T. Turanyi // Reliability Engineering and System Safety. - 1997. - Vol. 57. - P. 41-48.

97.Turanyi, T. Sensitivity analysis of complex kinetic systems. Tools and application / T. Turanyi // Journal of mathematical chemistry. - 1990. - Vol. 5. - P. 203-248.

98.Нурисламова, Л.Ф. Методика получения редуцированной математической модели химической реакции / Л. Ф. Нурисламова, И. М. Губайдуллин // Системы управления и информационные технологии. - 2014. - Т. 3, № 57. -С. 266-271.

99.Нурисламова, Л.Ф. Редукция детальных схем химических превращений окислительных реакций формальдегида и водорода на основании результатов анализа чувствительности математической модели / Л. Ф. Нурисламова, И. М. Губайдуллин // Вычислительные методы и программирование. - 2014. - Т. 15. - С. 685-696.

100. Спивак С.И. Обратные задачи химической кинетики / С.И. Спивак, И.М. Губайдуллин, Е.В. Вайман. - Уфа: РИО БашГУ, 2003. - 110 с.

101. Померанцев, А.Л. Методы нелинейного регрессионного анализа для моделирования кинетики химических и физических процессов: дис. ... д-ра физ-мат. наук.: 01.04.17 / А.Л. Померанцев. - Москва, 2003. - 304 с.

102. Безменов, В.М. Фотограмметрия. Построение и уравнение аналитической фототриангуляции / В.М. Безменов. - Казань : КГУ, 2009. -87 с.

103. Raymond, F.M. Evaluation of objective functions for estimation of kinetic parameters / F.M. Raymond, T.C. Bradley // Medical Physics. - 2006. - Vol. 33, № 2. - C. 342-353.

104. Turanyi, T. Determination of rate parameters based on both direct and indirect measurements / T. Turanyi, T. Nagy, I.Gy Zsely, M. Cserhati, T.

Varga, B.T. Szabo, I. Sedyo, P.T. Kiss, A. Zcmpleni, H.J. Curran // Int. J. Chem. Kinet. - 2012. - Vol. 44, № 5. - P. 284-302.

105. Губайдуллин, И.М. Информационно-аналитическая система решения многопараметрических обратных задач химической кинетики: дис. ... д-ра физ.-мат. наук: 02.00.04/ И.М. Губайдуллин. -Уфа, 2012. - 243с.

106. Нурисламова, JI. Ф. Редукция математической модели реакции окисления водорода / Л. Ф. Нурисламова, И. М. Губайдуллин // Материалы XIV Международной конференции «Высокопроизводительные параллельные вычисления на кластерных системах (НРС 2014)». - Пермь, 2014. - С. 318-324.

107. Байназарова, Н. М. Исследование чувствительности результатов кинетического моделирования к константам скорости реакции / Н. М. Байназарова, Л. Ф. Нурисламова, И. М. Губайдуллин // Журнал Средневолжского математического общества. - 2014. - Т. 16, № 4. - С. 23-32.

108. Анализ чувствительности кинетической системы к варьированию входных параметров / Н. М. Байназарова, Л. Ф. Нурисламова, И.М. Губайдуллин // Свидетельство о регистрации электронного ресурса № 19919 от 07.02.2014. - Москва: ИНИПИ РАО ОФЭРНиО. № 19919.

109. Глобальный анализ чувствительности и неопределенности кинетической системы: / Н. М. Байназарова, Л. Ф. Нурисламова, И. М. Губайдуллин // Свидетельство о регистрации электронного ресурса № 2014617100 от 11.07.2014. - Москва: Реестр программ для ЭВМ.

110. Вычисление локальных показателей чувствительности кинетической системы / Н. М. Байназарова, Л. Ф. Нурисламова, И. М. Губайдуллин // Свидетельство о регистрации электронного ресурса № 2014617098 от 11.07.2014. - Москва: Реестр программ для ЭВМ.

111. The Kintecus simulation software [Электронный ресурс]. - Режим доступа: http://www.kintecus.com/.

112. Hidaka, Y. Shock-tube study of CH20 pyrolysis and oxidation / Y. Hidaka, T. Taniguchi, H. Tanaka, et al. // Combustion and Flame. - 1993. - Vol. 92, № 4.-P. 365-376.

113. O'Connaire, M. A Comprehensive modeling study of hydrogen oxidation / M. O'Connaire, H. J. Curran, J. M. Simmie, W. J. Pitz, С. K. Westbrook // Int. J. Chem. Kinet. - 2004. - Vol. 36. - P. 603-622.

114. NIST Chemistry WebBook [Электронный ресурс]. - Режим доступа: http://webbook.nist.gov/chemistry/.

115. Masel, R.I. Chemical Kinetics and Catalysis / R.I. Masel - New York: Wiley Interscience, 2001. - 952 p.

116. Nurislamova, L. F. Few-Step Kinetic Model of Gaseous Autocatalytic Ethane Pyrolysis and Its Evaluation by Means of Uncertainty and Sensitivity Analysis / L. F. Nurislamova, O. P. Stoyanovskaya, O. A. Stadnichenko, I. M. Gubaidullin, V. N. Snytnikov, A. V. Novichkova // Chemical Product and Process Modeling. - 2014. - 9(2) - P. 143-154.

117. Nurislamova, L. F. Sensitivity analysis of kinetic model parameters of the gas-phase ethane pyrolysis reactions / L. F. Nurislamova O. P.Stoyanovskaya, N. M. Bainazarova, O. A.Stadnichenko, I. M. Gubaydullin, V. N. Snytnikov // Proc. the international conference «Advanced Mathematics, computations and applications». - Novosibirsk, 2014. - P. 82.

118. Стадниченко, O.A. Математическое моделирование потоков многокомпонентного газа с энергоемкими химическими процессами на примере пиролиза этана / О.А. Стадниченко, В.Н. Снытников // Вычислительные методы и программирование. - 2014. - Т. 15. - С. 658668.

119. Нурисламова, Л.Ф. Математическое моделирование реакции пиролиза этана Л. Ф. Нурисламова Н. М. Байназарова // Материалы всероссийской научно-практической конференции «Математическое моделирование на основе методов Монте-Карло». — Бирск, 2013. - С. 34 - 42..

120. Нурисламова Л.Ф. Математическое моделирование процесса газофазного дегидрирования этана / Л. Ф. Нурисламова Н. М. Байназарова // Сборник трудов II Всероссийской научно-практической конференции с международным участием «Математическое моделирование процессов и систем». - Стерлитамак. - 2013. - С. 292-296.

121. Katare, S. An intelligent system for reaction kinetic modeling and catalyst / S. Katare, J.M. Caruthers, W.N. Delgass, V. Venkatasubramanian // Ind. Eng. Chem. Res. - 2004. - Vol. 43. - P. 3484-3512.

122. Keyvanloo, K. Genetic algorithm model development for prediction of main products in thermal cracking of naphtha: comparison with kinetic modeling / K. Keyvanloo, M. Sedighi, J. Towfighi // Chem. Eng. J.- 2012. - Vol. 209. -P. 255-262.

123. Tsuchiya, T. Advantages of external periodic events to the evolution of biochemical oscillatory reactions / T. Tsuchiya, J. Ross // J. Phys. Chem. -2001.-Vol. 105.-P. 4052-4058.

124. Khan, U. Pyrolysis of propane under vacuum carburizing conditions: an experimental and modeling study / U. Khan, S. Bajohr, D. Buchholz, R. Reimert, H.D. Minh, K. Norinaga, V.M. Janardhanan, S. Tischer, O. Deutschmann // J. Anal. Appl. Pyrolysis. - 2008. - Vol. 81. - P. 148-156.