УДК 65.011.56
ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ДЛЯ ПОСТРОЕНИЯ АЛГОРИТМОВ ПУСКА И ПРОТИВОАВАРИЙНОЙ ЗАЩИТЫ ХИМИЧЕСКИХ РЕАКТОРОВ
А.Э. Софиев, И.А. Янкина
Кафедра «Техническая кибернетика и автоматика»,
ГОУ ВПО «Московский государственный университет инженерной экологии»; yankina.irina@gmail.com
Представлена членом редколлегии профессором В. И. Коноваловым
Ключевые слова и фразы: полимеризационный реактор; пуск; устойчивость.
Аннотация: Рассмотрено применение анализа устойчивости объекта для осуществления противоаварийной защиты при проведении пуска химического реактора непрерывного действия в производстве полиэтилена высокого давления.
Обозначения
Л, - предэкспоненциальный множитель реакции инициирования, м3/(кг-с);
Л0 - предэкспоненциальный множитель реакции обрыва цепи, м3/(кг-с);
Лр - предэкспоненциальный множитель
реакции роста цепи, м3/(кг-с);
ср - удельная теплоемкость, Дж/(кг-К);
Е - общая энергия активации, Дж/моль;
Е, - энергия активации реакции инициирования, Дж/моль;
Е0 - энергия активации реакции обрыва цепи, Дж/моль;
Ер - энергия активации реакции роста цепи, Дж/моль;
10 - начальная концентрация инициатора, кг/м3;
М0 - начальная концентрация мономера,
кг/м3;
Р - давление, Па;
Qp - тепловой эффект реакции, Дж/моль;
К - универсальная газовая постоянная, Дж/(моль-К);
Т - температура в реакторе, °С;
Т0 - температура на входе в полимериза-ционный реактор, °С;
Д V - объем активации, Дж/(Па-моль);
АУа' - объем активации реакции инициирования, Дж/(Па-моль); р - плотность реакционной смеси, кг/м3; т - время пребывания в реакторе, ч.
Химические реакторы являются основным оборудованием в технологических схемах получения многих веществ. Особенностью экзотермических химических реакторов непрерывного действия является наличие неустойчивого состояния системы [1]. Эта особенность является определяющей для процессов, протекающих с наибольшей эффективностью именно в неустойчивых режимах. Возможность нормального функционирования подобных систем в неустойчивом режиме обеспечивают системы автоматического регулирования.
Однако перед включением системы автоматического регулирования необходимо вывести объект в рабочий режим, что осуществляется в процессе пуска (или останова) реактора.
Задача исследования устойчивости экзотермических химических реакторов появляется на различных стадиях создания как технологического процесса, так и системы управления реакторами (рис. 1). К таким задачам относятся: конструиро-
ISSN 0136-5835. Вестник ТГТУ. 2011. Том 17. № 1. Transactions TSTU.
м
З
>.
с
о
■=t
идеальное смешение физ-хим. нрамеїрьі иосюянны мі иовеииое виранииваиие всех параметров (Т, С)
допущении нринимаюіси в зависимое і и от модели линеаризация в окресі носій сосі ом II им равновесия
идеальное смешение физ-хим. нрамеїрм носюннны мгновенное выравнивание всех параметров (Т, С)
Статика
Статика, динамика
Динамика
- :
О ft.
ОДУ. ДУ в час і мы \
11|ШМ «НО 1НЬ1\, ИН IV! |М1-інффг|і«н мнн.п.нмс у ракнгнмн
ОДУ
РИС, РИВ, реакморы промежуточною режима смешении (нолнаи сеі реї ации, максимум смешения, нолусеї регаиия и гл.)
'ч
РИС
Рис. 1. Применение анализа устойчивости при исследовании химических реакторов
вание химических реакторов, построение систем автоматического управления и противоаварийной защиты. При проектировании реакторов с помощью анализа устойчивости решается задача выбора оптимальных геометрических размеров реактора [2]. Анализ устойчивости позволяет решить задачу выбора структуры системы автоматического управления: организовать управление по принципу обратной связи или выбрать более подходящий для данного объекта тип регулирования, например, стабилизацию неустойчивого режима в полимеризационном реакторе с помощью периодических пульсаций управляющего параметра [3].
Во всех указанных случаях при исследовании устойчивости используется математическая модель химического реактора. В зависимости от вида решаемой задачи применяются модели, получаемые при разных допущениях.
Так как основной режим функционирования химических реакторов - установившийся, поэтому для решения задач, связанных с оптимальным проектированием реактора, используют модели статики. Данные модели описываются системой нелинейных алгебраических уравнений в случае допущения о реакторе идеального смешения или системой обыкновенных дифференциальных уравнений - в случае реактора идеального вытеснения. При построении систем автоматического управления и противоаварийной защиты используют как модели статики, так и модели динамики, например, при построении систем управления работой периодических реакторов. С помощью моделей статики исследуют область устойчивости системы, что также необходимо при построении систем противоаварийной защиты химического реактора, так как помогает выявить граничные значения параметров процесса, при достижении которых система выйдет в неустойчивое состояние [4].
Процесс полимеризации этилена под высоким давлением относится к числу наиболее опасных технологических процессов. Данный процесс протекает в автоклавном реакторе непрерывного действия и характеризуется высоким рабочим давлением, достигающим 150 МПа, и температурой, рабочие значения которой находятся в интервале от 240...270 °С.
Поэтому задача построения алгоритмов управления пуском и противоаварийной защиты полимеризационного реактора является актуальной.
Ввиду сложности рассматриваемого в качестве объекта управления реактора полимеризации этилена под высоким давлением используются методы математического моделирования для его исследования. Процесс полимеризации этилена под высоким давлением в автоклавном реакторе является хорошо изученным объектом, важнейшие результаты которого приведены в работах Р.1 Нойу7ег, Т.К Zwiete-пп§, Б.В. Вольтер, М.Г. Слинько.
Ранее была предложена последовательность пуска полимеризационного реактора, показанная на рис. 2 [5]. Данная последовательность пуска основана на применении математической модели реактора полимеризации этилена под высоким давлением.
Рис. 2. Блок-схема пуска полимеризационного реактора
Уравнения материального баланса по инициатору (1) и теплового баланса (2) описывают модель идеального смешения, получаемую при допущениях о постоянстве физико-химических параметров процесса (таких, как теплоемкость и плотность), квазистационарном протекании процесса и равенстве начальной и текущей концентраций продукта в реакторе вследствие малой степени превращения этилена в полиэтилен (примерно 10 %), а также отсутствии теплообмена с рубашкой реактора (вследствие толщины стенки порядка 20 см).
Существенной особенностью пускового режима непрерывного полимериза-ционного реактора является наличие гистерезиса - скачкообразного изменения температуры в реакторе, приводящего, в большинстве случаев, к срабатыванию аварийной защиты или взрыву реактора или, наоборот, к резкому снижению температуры и, как следствие, гашению реакции.
Наличие гистерезиса было выявлено с помощью статической характеристики полимеризационного реактора [6], приведенной на рис. 3.
Для своевременного обнаружения и предупреждения аварийных ситуаций во время пуска реактора предлагается дополнить алгоритм пуска проверкой устойчивости системы (рис. 4).
Текущие данные о состоянии процесса, а именно: температура, давление, расход инициатора в реактор (концентрация инициатора в реакторе), поступают на вход алгоритма и используются для идентификации статической математической моделим реактора Т = /То), получаемой из уравнений (1), (2) при = 0
Т, °С 300 т
dt
(E+av-p) ,
'AfMgLe RT + -(Ig _ L);
т
dt Pcp I Ao
_TE+AVaP) mUlG,5e RT + - To _ T).
T
(2)
25G
2GG
15G ■
120 140 160 180 200 220 240 260 T0, °C
Рис. 3. Статическая характеристика полимеризационного реактора [4, 6] при различных значениях начальной концентрации инициатора, I0, кг/м3:
1 - 0,08; 2 - 0,03; 3 - 0,01
Рис. 4. Алгоритм пуска реактора полимеризации этилена на основе анализа устойчивости объекта
На следующем шаге рассчитываются функции чувствительности по основным управляющим параметрам (начальной температуре Т0, начальной концентрации инициатора 10, давлению Р): £ .■ = То ^; где а - вектор управляющих парада у
метров. После чего осуществляется сравнение значений производных, полученных на текущем и предыдущем шаге с целью выявления бифуркационного со-
стояния системы (точки А или В на рис. 3). Если
дТ0
>>
о V
кдТ0 V
следова-
к-І
тельно система находится в окрестности точки бифуркации, что характеризуется как предаварийное состояние вследствие возможной потери устойчивости системы. В этом случае включается процедура противоаварийной защиты реактора, в ходе которой выполняются следующие действия: определение расхода инициатора на входе в реактор (кривые 1, 3, рис. 3), изменяющего концентрацию инициатора в реакторе; перевода системы в дремлющий режим путем снижения давления в реакторе до 100 МПа или выброс содержимого реактора в атмосферу и аварийный останов реактора.
k
В случае неравенства нулю значения производной по температуре от статической характеристики реактора и отсутствия предаварийного состояния реактора осуществляется расчет пусковых управлений с учетом заданных технологических ограничений и реализация этих управлений на объекте.
Таким образом, на основе использования математической модели был предложен алгоритм управления пуском полимеризационного реактора в производстве полиэтилена высокого давления. Использование математической модели позволило учесть особенности переходных режимов автоматизируемого объекта.
Применение системы управления пуском реактора, построенной на основе приведенного алгоритма, позволяет снизить время вывода реактора на рабочие значения параметров процесса за счет уменьшения времени простоя реакторного блока в результате попадания в неустойчивые режимы по время пуска реактора. Кроме этого, использование данного алгоритма пуска и противоаварийной защиты реактора позволяет уменьшить потери сырья (сброс в атмосферу в результате срабатывания аварийной защиты) и исключить аварийные ситуации в процессе пуска.
Список литературы
1. Арис, Р. Анализ процессов в химических реакторах / Р. Арис. - М. : Химия, 1967. - 328 с.
2. Янкина, И.А. Конструирование нефтехимического реактора на основе математического моделирования (на примере реактора полимеризации этилена) / И. А. Янкина, М.Г. Беренгартен, А.Э. Софиев // Нефтехим. и газ. машиностроение. -2010. - № 7. - С. 18-20.
3. Софиев, А.Э. Вибрационная стабилизация неустойчивого режима химического реактора с перемешиванием / А.Э. Софиев, А.М. Трахтенберг // Теорет. основы хим. технологии. - 1989. - Т. 23, № 5. - С. 644-650.
4. Вольтер, Б.В. Устойчивость режимов работы химических реакторов / Б.В. Вольтер, И.Е. Сальников. - 2-е изд., перераб. и доп. - М. : Химия, 1981 - 200 с.
5. Вольфсон, С.А. Расчеты высокоэффективных полимеризационных процессов / С. А. Вольфсон, И.С. Ениколопян. - М. : Химия, 1980. - 312 с.
6. Голованов, О.В. О гистерезисе в статической характеристике процесса полимеризации этилена / О.В. Голованов, А.Э. Софиев // Тр. Центр. науч.-исследов. ин-та комплексной автоматизации. - 1964. - Вып. 10. - С. 232-236.
Application of Mathematical Modeling for Start-up Algorithm and Emergency Protection of Chemical Reactors
A.E. Sofiev, LA. Yankina
Department of Technical Cybernetics and Automation,
Moscow State University of Environmental Engineering; yankina.irina@gmail.com
Key words and phrases: polymerization reactor; start-up; stability.
Abstract: The paper considers the application of stability analysis of the object for the emergency protection during start-up of a chemical reactor of continuous action in the production of high-pressure polyethylene.
Anwendung der matematischen Modellierung ftir den Bau der Algorithmen der Inbetriebnahme und des Notschutzes der chemischen Reaktoren
Zusammenfassung: Es ist die Anwendung der Analyse der Standfestigkeit des Objektes fur die Verwirklichung des Notschutzes bei der Durchfuhrung der Inbetriebnahme des chemischen Reaktors des ununterbrochenen Funktionierens bei der Produktion des Polyathylens des Hochdruckes betrachtet.
Application du modelage mathematique pour l’algorithme du lancement et de la securite anti-avarie des reacteurs chimiques
Resume: Est examinee l’application de l’analyse de la stabilite de l’objet pour la realisation de la securite anti-avarie lors du lancement du reacteur chimique de l’action continue dans la production du polyethilene de la haute pression.
Авторы: Софиев Александр Эльхананович - доктор технических наук, профессор, заведующий кафедрой «Техническая кибернетика и автоматика»; Янкина Ирина Аркадьевна - аспирант кафедры «Техническая кибернетика и автоматика», ГОУ ВПО «Московский государственный университет инженерной экологии», г. Москва.
Рецензент: Гатапова Наталья Цибиковна - доктор технических наук, профессор, заведующая кафедрой «Химическая инженерия», ГОУ ВПО «ТГТУ».