Файл: Кафаров, В. В. Принципы математического моделирования химико-технологических систем (введение в системотехнику химических производств) учеб. пособие.pdf

перед нагревом в печи с целью доведения температуры смеси почти до температур

в которой разделяются и распределяются крекинг-продукты.

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

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

всепаратор, в котором отделяется

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

Задача оптимизации заключается в том, чтобы максимизировать текущую прибыль ХТС. Прибыль является явной функцией следующих переменных: стоимости продукта, стоимости перерабатываемой нефти и стоимости перера­

Ограничения на зависимые переменные ХТС имеют следующий вид: выход

ками (ТКЗ) sg 691 ~ С; мольная доля 0 2 в газе нижней зоны регенератора (Х 0 2) ^ ^ 0,05; массовая доля остаточного кокса или катализатора (ХКК) ^ 0,01;

Температура катализатора в сепараторе........................Х12; Z22

Декомпозиция ХТС каталитического крекинг-процесса на две подсистемы (реактор и печь регенератора) осуществлена посредством разрыва потоков между

Рис. ѴІ-13. Декомпозиция задачи оптимизации ХТС каталитического крекинга.

реактором и печью, печью и нижним бункером катализаторного пневмоподъем­ ника (рис. ѴІ-13). Фактически этой декомпозицией достигнуто немного: исходная полная задача имеет восемь, а подзадачи I и II — пять и семь оптимизирующих переменных соответственно (рис. ѴІ-13). Таким образом, размерность задачи уменьшается незначительно.

Большое преимущество декомпозиции обычно заключается в том, что избе­ гают расчета рецикла. Однако в данном случае это преимущество не было ис­ пользовано, так как примененная для исследования описываемой ХТС регрес­ сионная модель уже включала расчет рецикла.

Две подзадачи оптимизации, связанные с двумя подсистемами, будут:

322

Решение этих подзадач составляет первую ступень двухуровневого алго­ ритма, являющегося результатом декомпозиции. Вторая ступень предназначена для согласования решений подзадач с целью решения полной задачи оптимиза­ ции ХТС.

Задача согласования заключается в минимизации двойственной функции, которая дана выражением

Y =/? + /; + Pi ( Z U ~ *Й) + Рг (ZUXU) +

+ Рз (zb-xU)+P4 (ZU-XU)

Верхние индексы «звездочка» у переменных указывают оптимальные зна­ чения, полученные в подзадачах. Задачей согласования*является изменение Pj, Р 2, Р3 и Р4 для минимизации у. Градиенты двойственной функции определяются следующим образом:

ду

дР1 — zu - X;

21*

Г Л А В А VII

СПЕЦИАЛЬНЫЕ ПРОГРАММЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ХИМИКО-ТЕХНОЛОГИЧЕСКИХ СИСТЕМ

Рассмотренные ранее принципы математического моделирования и алгоритмы оптимизации стратегии исследования ХТС являются теоретической базой для разработки специальных программ мате­ матического моделирования их с помощью ЦВМ.

1. ПРИН СПЕЦИАЛЬНЫХ ПРОГРАММ МОДЕЛИРОВАНИЯ ХТС

Возможная блок-схема организации решения задачи математи­ ческого моделирования ХТС представлена на рис. ѴІІ-1.

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

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

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

Возможная организация работы научно-исследовательского кол­ лектива, занимающегося решением задачи математического модели­ рования ХТС, приведена на рис. ѴІІ-2.

В обобщенную специальную программу моделирования ХТС входят: подпрограмма ввода исходной информации; подпрограмма математических моделей элементов системы; основная исполнитель­ ная подпрограмма; подпрограмма массива информации о физико-хи­ мических константах и физических свойствах компонентов и смесей; подпрограмма оптимизации и прогнозирования возможных техноло­ гических режимов; подпрограмма обеспечения сходимости вычисли­ тельных операций; подпрограмма вывода результатов.

324

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

Исходную информацию вводят в ЦВМ либо с помощью специаль­ ных карт, либо посредством свободного формата. Оба метода имеют достоинства и недостатки.

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

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

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

В топологической матрице каждой г-ой вершине параметриче­ ского потокового графа соответствует г'-ая строка, элементы которой имеют следующие значения: іі — номер вершины графа; і2 — наи­ менование отвечающего этой вершине графа типа элемента; /3+ -f- і/г — номера ветвей графа, соответствующих входным технологи­ ческим потокам данного элемента, со знаком плюс; г (п -f- 1) -f- іт— номера ветвей графа, отвечающих выходным технологическим пото­ кам данного элемента, со знаком минус.

Наиболее удобный и перспективный метод представления техно­ логической и информационной топологии ХТС состоит в примене­ нии информационно-потоковых мультиграфов. Информацию о топо­ логических характеристиках параметрического потокового графа или информационно-потокового мультиграфа представляют с по­ мощью матрицы ветвей графа [LL

326

Информацию о параметрах технологических потоков сырья и про­ дуктов XTG представляют в виде таблицы параметров технологиче­ ских потоков. Информацию о технологических и конструкционных параметрах каждого элемента ХТС представляют в форме таблицы параметров определенного типа элементов.

Подпрограмма математических моделей элементов ХТС строится по модульному принципу, сущность которого заключается в следу­ ющем. Математическую модель каждого элемента получают в виде совокупности математических моделей типовых технологических

Некоторые модули, используемые при моделировании ХТС, не соответствуют реально существующим элементам системы. К такого типа модулям относятся, например, модуль-калькулятор стоимости; модуль расчета материальных и тепловых балансов системы; модуль эквивалентного преобразования единиц измерения физико-химиче­ ских величин и др.

Как правило, один элемент ХТС может быть описан совокупно­ стью нескольких модулей. Так, многослойный контактный реактор при моделировании ХТС производства серной кислоты представляет­ ся математической моделью в виде совокупности нескольких модулей химического превращения, нагрева и смешения — разделения.

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

Какие именно модули выбираются для моделирования отдельных элементов, зависит от поставленных целей исследования системы, глубины понимания физико-химических основ технологических про­ цессов и точности исходных данных. Основой для разработки под­ программы математических моделей элементов ХТС по модульному принципу является библиотека стандартных программ математиче­ ских моделей типовых технологических операторов и операторная схема системы.

Общая стратегия разработки модулей включает следующие этапы: а) точность результатов моделирования ХТС в целом; б) точность индивидуальных модулей для каждого элемента; в) уровень стандар­

о точности модулей зависит от того, изучается отдельный элемент или система в целом. Например, можно получить точные результаты моделирования ХТС для стоимостных и экономических оценок при

327