автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.06, диссертация на тему:Автоматизация сложных химико-технологических производств

На правах рукописи

Шауро Василий Сергеевич

АВТОМАТИЗАЦИЯ СЛОЖНЫХ ХИМИКО-ТЕХНОЛОГИЧЕСКИХ ПРОИЗВОДСТВ (на примере производств фенола-ацетона и полиэтилена в реакторах смешения)

05.13.06- Автоматизация и управление технологическими процессами и производствами (химическая технология, нефтехимия и нефтепереработка, биотехнология)

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Москва - 2006.

Работа выполнена на кафедре Технической кибернетики и автоматизации в Московском государственном университете инженерной экологии

Научный руководитель:

доктор технических наук, профессор Софиев Александр Эльхананович

Официальные оппоненты:

Ведущая организация:

доктор технических наук, профессор

Ромаш Эдуард Михайлович

доктор технических наук, профессор

Бессарабов Аркадий Маркович

Институт программных систем РАН г.Переславль-

Залесский

Защита диссертации состоится » К Р-* с/^Э 2006 г в часов на заседании диссертационного совета Д 212.145.02 при Московском государственном университете инженерной экологии, 105066, г. Москва, ул. Старая Басманная, д. 21/4.

С диссертацией можно ознакомиться в библиотеке МГУИЭ.

Автореферат разослан «¿О» ОиЛЗ^с^ ^ 2006 г

Ученый секретарь диссертационного совета Д 212.145.02 к.т.н., доцент

Мокрова Н.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность работы. Среди комплекса основных задач, решаемых на химико-технологических производствах для обеспечения их эффективной и бесперебойной работы, особое место занимают задачи автоматизации, как производства в целом, так и отдельных его участков. Эффективное решение этих задач - один из наиболее значимых факторов, определяющих общую эффективность управления химическим предприятием. Автоматизированные системы управления производством, являются основой для рентабельного функционирования производства, служат фундаментом для построения систем управления предприятием в целом.

Многие российские химические производства не оснащены, в достаточной степени, высокоточной информационно-измерительной техникой и АСУ ТП. .В последнее время задача автоматизации химических производств становится приоритетной для большинства компаний работающих в химической промышленности.

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

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

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

Научная новизна. В результате проведенной работы были получены следующие основные научные результаты:

1. Впервые разработана комбинированная система управления автоклавным реактором полимеризации на основе вибрационного регулирования и регулирования по обратной связи.

2. Разработана и предложена инженерная методика применения типовых решений при создании АСУТП.

Практическая ценность и реализация результатов работы.

1. Принят к внедрению проект АСУТП фенола-ацетона

разработанный на основе методики применения типовых решений.

Апробация работы. Результаты работы докладывались и обсуждались на научно-технических советах ОАО "ЦНИИКА", Международной конференции «Программные системы: теория и приложения», ИПС РАН (г. Переславль-Залесский ,■ май 2004), Международной научно-практической конференции «Интеллектуальные системы принятия решений и информационные технологии» (Украина, г. Черновцы 19-21 мая 2004 г.).

Публикация результатов исследования. Основные результаты научных исследований по теме диссертации содержатся в 4 печатных работах.

Структура и объем работы. Диссертация состоит из введения, 6-и глав, заключения и 4-х приложений. Работа изложена на 122 страницах машинописного текста, содержит 50 рисунков и 4 таблицы и приложения объемом 17 страниц. Список литературы включает 72 наименования.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В главе 1 приведена общая характеристика комплексного подхода к построению АСУТП крупного производства, приводятся примеры существующих АСУТП в различных отраслях промышленности, спроектированных в соответствии с таким подходом.

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

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

Смысл комплексного подхода к автоматизации производства - обеспечить руководство предприятия средствами контроля работы и рычагами воздействия для обеспечения максимальной прибыли от функционирования объекта. Основные положения теории комплексного подхода к автоматизации приведены по работам О.В. Голованова, Э.Л. Ицковича, Е.П. Стефани, Ю.З. Шапиро и др.

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

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

В представленной работе рассматриваются задача унифицированного подхода к построению АСУ крупнотоннажным производством, на примере производств фенола-ацетона и полиэтилена высокого давления, и возникающая при решении общей задачи, подзадача управления отдельным сложным неустойчивым объектом, на примере реактора полимеризации этилена. Другими словами задачи уровней 1, 2 и 3 на рисунке 1.

Автоматизированные системы управления химическим производством являются частным случаем автоматизированной системы управления технологическим процессом. Ее основными составляющими являются программно технический комплекс (ПТК), набор документации и обслуживающий персонал.

Современная АСУ 111 решает комплекс задач:

A) Информационные, включающие в себя централизованный сбор информации о параметрах технологического процесса и состояния оборудования, контроль отклонений основных технологических параметров, контроль состояния оборудования, контроль ТЭП, подготовку информации для передачи на более высокие уровни MES, ERP, OLAP.

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

B) Вспомогательные, включающие в себя протоколирование действий обслуживающего персонала, диагностирование оборудования.

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

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

Выделяется три класса элементов подлежащих унификации:

- Техника - проектные решения по выбору и применению НТК.

- Задача - комплекс алгоритмов и программ, реализующих функции системы.

- Персонал - инструкции, регламентирующие действия персонала в условиях функционирования АСУ.

Применение такой унификации для всех подзадач, полученных после декомпозиции общей задачи, позволяет сделать разработку АСУТП в значительной степени независимой от выбора технических и программных средств.

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

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

Задача построения АСУТП при унифицированном подходе разбивается на три подзадачи: организация системы сбора информации, организация системы обработки и хранения информации, организация представления информации человеку — задача человеко-машинного интерфейса HMI. Каждая из этих задач решается отдельно, независимо от других.

Обработка и хранение информации -структурированные кабельные сети на основе серверов

Система сбора информации - контроллеры, датчики

Рис.2 Уровни АСУТП

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

Задача решается в соответствии с унифицированным подходом. На каждом уровне предлагаются свои типовые решения.

Уровень датчиков и исполнительных механизмов.

Датчики предоставляют информацию на верхние уровни автоматической системы. Качество и надежность работы АСУТП невозможно обеспечить без достоверной и полной информации о текущем состоянии объекта, за которую отвечают датчики.

Для решения задачи предложено опустить некоторый интеллект на уровень датчиков, заменив интерфейс 4...20 мА на промышленную сеть общепринятых стандартов HART и Profibus-PA. В результате вместо большого количества 2-х, 3-х и 4-х проводных линий связи, идущих от множества датчиков и исполнительных механизмов к каналам ввода-вывода PLC, предлагается один "малопроводный" кабель, подключаемый к соответствующему интерфейсному модулю программируемого логического контроллера (PLC).

Уровень сбора информации.

Система сбора информации строится исходя из требований надежного приема сигналов с уровня датчиков и передачи их в микроконтроллер. Главным ограничением здесь является территориальная распределенность автоматизируемых объектов и места куда сводятся кабельные трассы от датчиков и исполнительных механизмов.

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

Для производства фенола-ацетона выделяется три удаленных объекта: 101 103 и 602. Отдельно выделяется система противоаварийной защиты, которая по существующим нормам должна быть независима от АСУТП и построена со строгим соблюдением свойства отказоустойчивости.

Для производства полиэтилена выделяется пять удаленных друг.от друга объектов: корпус грануляции линий А и Б, корпус реактора линии А, корпус реактора линии Б, корпус компрессии линий А и Б, центральный пульт управления (ЦПУ). Отдельно выделяется система противоаварийной защиты (ПАЗ), которая по существующим нормам должна быть независима от АСУТП. Расстояние между корпусами достигает нескольких километров.

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

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

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

Исходя из технологических особенностей рассматриваемых производств, а также требований надежности и поэтапности внедрения, техническая структура IIIК для АСУТП разделена на подсистемы следующим образом:

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

Для производства полиэтилена: прием энергоресурсов, компрессия этилена, реакторный блок, возвратный газ высокого давления, возвратный газ низкого давления, грануляция, станция горячей воды, противоаварийная защита производства.

Чтобы связать сигналы от подсистем в единую систему ввода-вывода необходимо выбрать технологию полевой шины. В рассматриваемых системах выбор оборудования фирмы Siemens обуславливает применение технологии Profibus.

Пример такой сети для АСУТП производства фенола-ацетона приведен на рисунке 3, для АСУТП производства полиэтилена - на рисунке 4.

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

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

В АСУТП производства полиэтилена сеть Profibus АСУ аналогична по исполнению сети производства фенола-ацетона. Основное отличие двух систем в построении сети ПАЗ. Станции ввода-вывода ПАЗ на производстве полиэтилена располагаются на удалении нескольких сот метров от контроллеров, что предъявляет особые требования к построению сети.

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

о. о

I £

J

Рис. 3 Структура сети Profibus производства фенола-ацетона

Контроллеры отделения 423

Контроллеры отделения 424

Оптическая Сеть Profibus DP топология кольцо

Оптическая Сеть Profibus DP топология кольцо

Удаленные станции ввода вывода

Удаленные станции ввода вывода

Рис. 4 Структура сети Profibus для стадии полимеризации

Уровень обработки и хранения информации.

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

Для производства фенола-ацетона сеть имеет топологию кольцо, построенную с помощью модулей OSM (Optical Switch Module) для усиления сигнала при передаче его на большие расстояния и реализации кольца. Кольцо имеет протяженность около километра и частично проходит по открытой местности, что обуславливает применение оптической среды передачи данных. Контроллеры подключены к модулям OSM с помощью электрической сети. Для питания модулей OSMXESM использована встроенная в модули схема резервированного питания.

Производство полимеризации имеет более сложную структуру сети контроллеров. Сеть Ethernet имеет структуру двух резервировано-соединенных

оптоволоконных колец, кольца ядра системы и кольца полимеризации, связь осуществляется через специальные порты OSM ITP53, предназначенные для соединения колец. При этом, в кольце ядра резервированы OSM, к которым подключается кольцо полимеризации.

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

Для этого проводится исследование объекта. Заметим, что методика исследования также не зависит от объекта и применима к любому производству.

Исследование объекта заключается в:

A) Определении общего количества сигналов, которые должна обрабатывать АСУТП, и количества регуляторов, с выделением одноконтурных и каскадных систем регулирования.

Б) Определении уровней сигналов и циклов их обработки.

B) Определении территориальной распределенности системы.

Г) Алгоритмизации стандартных задач первичной обработки информации, регулирования, создания интерфейсов оператора-технолога.

Д) Выделении нестандартных задач, требующих отдельной алгоритмизации.

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

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

На основе примененных решений строится общая структура системы.

Операторские станции АСУ

Операторские станции ПАЗ

Сеть Ethernet контроллеров

Контроллеры АСУ

Сеть РгоАЬиз АСУ -

Станции ввода/вывода АСУ Сети датчиков

Резервированный сервер

Контроллеры ПАЗ

- Сеть Profibus ПАЗ

Станции ввода/вывода ПАЗ

Рис.5 Техническая структура производства фенола-ацетона

Уровень НМ1

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

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

A) Верхний уровень: базовые обзорные видеограммы объектов, отражающие технологическую структуру данного объекта.

Б) Средний уровень: базовые обзорные видеограммы - мнемосхемы для отдельных технологических стадий, узлов и отделений.

B) Нижний уровень: к видеограммам нижнего уровня относятся: панели контроля и регулирования, группы значений взаимосвязанных контролируемых и регулируемых технологических параметров, тренды, видеопротоколы.

Рассмотренный подход применен к двум различным производствам: фенола-ацетона и полиэтилена.

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

В главе 3 рассматривается задача управления реактором полимеризации этилена с перемешивающим устройством в неустойчивом режиме.

Приведен краткий обзор существующих типов реакторов полимеризации этилена, дано краткое описание кинетики реакции.

Управление реактором полимеризации сложная задача, реактор является сложным нелинейным объектом, работающим в неустойчивом режиме.

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

Задачу управления реактором предложено решать с применением метода вибрационного регулирования. Этот метод впервые был предложен С. М. Мерковым, в дальнейшем применение этого метода к химическим реакторам рассматривалось в работах В.В. Кафарова, А.Э. Софиева, Р. Беллмана. Дан литературный обзор основных положений теории вибрационного регулирования и применения вибрационного регулирования для промышленных установок. Приведено обоснование целесообразности использования метода вибрационного регулирования по сравнению с традиционными системами управления по обратной связи. Метод вибрационного регулирования позволяет стабилизировать объект, в случае, когда измерение параметров объекта затруднено. Широкого применения на практике этот метод не нашел, из-за сложности настройки системы регулирования и колебаний значения стабилизируемого параметра, иногда достаточно

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

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

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

продукт

Рис.6 Принципиальная схема вибрационного управления реактором

полимеризации

1 - Реактор, 2 - Датчик температуры, 3 - Регулятор температуры, 4 -дозировочный насос, 5 - выгрузной клапан, 6 - датчик давления, 7 — регулятор давления, 8 — генератор колебаний.

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

РЮ | Ртах

▲

Оп(Т

О о Л} У» / • ' 1 1 ---ь.

Ре Рт1п

Те т

Рис.7 Стационарное состояние реактора СЦТ,Р<;) - скорость тепловыделения, СЬ(Т) — скорость теплоотвода

1

1

Рис.8 Введение колебаний по давлению

0 - к

Оп ^^^ЬпСГз.Р)

1 1 Опэ | 1

1 1 -1-

Рггйп Рэ Ртах

Рт'т ре Ртах Р

Рис.9 Зависимость скорости тепловыделения от давления

Рис.10 Зависимость наклона кривой тепловыделения от Р

Стационарное состояние реактора, показанное на рис. 7 и соответствующее точке Т5 является неустойчивым. Зависимость скорости тепловыделения от параметра Р, при стационарной температуре Т^ в некоторых случаях имеет вид показанный на рисунке 9. А зависимость наклона кривой тепловыделения от параметра Р - вид представленный на рисунке 10.

При введении периодических колебаний параметра Р, с неизменным средним за период 9, с периодом малым по сравнению с характеристическим временем переходных процессов в самом реакторе (рис.8), среднее значение кривой тепловыделения изменяется.

Функция §(Т5,Р) выпукла вверх в окрестности Р8, а потому, в искусственно создаваемом периодическом режиме, среднее за период значение наклона кривой меньше своего стационарного значения и находится на выпуклой оболочке У^.Р) и множестве точек {§,Р}.

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

В главе 5 приведены основные результаты вычислительного эксперимента по моделированию системы управления реактором полимеризации с мешалкой.

Вычислительный эксперимент проводится с помощью математической модели химического реактора, предложенной в ряде работ Р. Ариса, Д. Перлмуттера, Б. В. Вольтера и др. Модель включает в себя уравнения материального и теплового балансов и имеет общий вид.

¿С, V , _ _.

1Г=7(С«>-С<)-Й,<

срУ^ = -МР-Тс)-срд(Т-Т.)

О)

Где V - объемная скорость подачи реагирующей смеси; V- объем реактора; Сю - концентрация ¡-го вещества на входе в реактор; С, - его концентрация внутри реактора и на выходе из него, Ц - тепловые эффекты реакций обладающих

скоростями Тс - температура стенки; Ь - коэффициент теплопередачи; 8 -поверхность теплообмена реактора с окружающий средой. Модель построена при следующих допущениях:

Эффективная скорость потока газа существенно больше скорости передачи тепла вдоль потока. Это предположение дает возможность пренебречь

слагаемыми вида —г и рассматривать гиперболическое уравнение теплового

дх

баланса, вместо параболического.

Пренебрегаем, вязкостью газа и связанными с ней потерями на внутреннее трение, а также весом газа, поскольку его движение обусловлено в основном перепадом давления, а не силой тяжести.

Горизонтальную скорость движения газа считаем малой по сравнению с вертикальной.

Тогда математическая модель реактора будет иметь имеет вид:

¿с 1

— = —к2(с,Т) +—(С0 -С) = Р(г,С0,С5 Т3)

Л Г (2)

^¡- = 1к1(с,Р,Т) + -(Т.-Т) +—(ТС -Т) = д(т,Тв,Тс,С5Т5) тех с

где С,С0- концентрация реакционной смеси соответственно на входе в реактор и в

реакторе; Тв Т,Те- температура на входе в реактор в реакторе и стенке реактора; Р —

давление в реакторе, г1(с,Г) = Л:1е"£'/яг[с|',,г2(£г5Г) = ^е'^'^^-скорости химических

реакций, Кь К2 - предэкспоненциальные множители; <7 - тепловой эффект реакции; Е ,Е2(Р) - £ ° - а (Р - Р0) - энергии активации а - удельный объем активации; п и

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

Далее проводится исследование устойчивости реактора, для этого чего находятся стационарные состояния.

Координаты стационарных состояний можно получить, решая систему уравнений

Р(г,Св,Св,Г8) = 0

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

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

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

Для стабилизации этого стационарного состояния предлагается применить к рассматриваемому реактору принцип вибрационного управления. Стабилизация достигается путем периодического изменения давления вида Р = Р3 + /(ох), где

- 1 Г°

а> = 27г/Т0 - частота колебания давления, причем Р = —

То о

Здесь Р1 - некоторое стационарное давление в реакторе. Как видно из исходных уравнений давление входит в тепловой баланс в степени экспоненты.

Е,{Р) = Е°-а{Р-Рй) (4)

Или, после введения колебаний:

£,(/>) = -а(Р!+Аая)-Р0) (5)

Введение колебаний приводит к следующим эффектам:

А) Скорость теплоотвода остается неизменной. Средняя скорость тепловыделения изменяется следующим образом: при малых значениях амплитуды скорость возрастает, что ухудшает поведение системы, а с ростом амплитуды средняя скорость тепловыделения начинает падать, и система приближается к области устойчивости. График изменения скорости представлен на рисунке 11.

Рис.11 Изменение скорости тепловыделения при росте амплитуды колебаний Для проверки поведения системы в окрестностях состояния равновесия, необходимо решить систему балансовых уравнений. При создании промышленной АСУ необходимо, чтобы решение системы проводилось на ЭВМ достаточно быстро, и позволяло рассчитывать и корректировать параметры колебаний в реальном времени. Оптимальным, с точки зрения быстродействия и легкости алгоритмизации, являются методы Рунге-Кутта 4-го порядка.

Решение системы балансовых уравнений с помощью метода Рунге-Кутта выполнено в пакете МАТКАД 11.

Получены следующие результаты:

Т(х)

-"150

0 50 100

х

Рис. 12 Решение системы для второй стационарной точки с введением

колебаний

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

В главе 6 рассматривается разработка приложения для проведения исследования эффективности различных систем регулирования.

Моделирование реактора и систем регулирования выполнено в среде МАТЪАВ с использованием специального пакета БшпЦтк.

Рис.13 Реализация модели реактора с системами управления в пакете БшшПпк

Исследуются три различных типа систем регулирования:

- Система регулирования по обратной связи с использованием ПИД-регулятора.

- Система регулирования с вибрационным регулированием.

- Комбинированная система регулирования, состоящая из ПИД-регулятора и вибрационного управления.

Результаты вычислительного эксперимента по управлению реактором для различных настроек ПИД-регулятора приведены на рисунке 14.

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

На рисунке 15 показан переходный процесс с использованием ПИД-регулятора, когда регулирование начинается при значении температуры близкой к заданию. Эксперимент показал, что при использовании рабочей точки, отстоящей от заданного значения температуры на величину порядка 5-10 градусов, замкнутая система регулирования позволяет достичь устойчивости состояния равновесия и при этом значения динамической ошибки и времени переходного процесса невелики.

Основная трудность при управлении реальным реактором состоит в том, каким образом достичь рабочей точки.';::'"'"-" ■ • V '

Рис.14 Переходный процесс замкнутой системы регулирования с ПИД-регулятором.

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

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

На рисунках 16 и 17 показаны переходные процессы при запаздываниях г1 и г2 соответственно. Вычислительный эксперимент показал, что качество регулирования ухудшается с ростом величины запаздывания, Существует критическая величина запаздывания, после которой система с обратной связью не способна стабилизировать систему. Величина динамической ошибки остается большой, даже при использовании рабочей точки.

Рис. 16 Переходный процесс с с запаздыванием г 1

Рис. 17 Переходный процесс запаздыванием г2

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

На рисунке 18 показано поведение системы при использовании вибрационного регулирования.

Рис.18 Изменение значения температуры при использовании вибрационного регулятора

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

При построении комбинированной системы для корректной работы ПИД-регулятора необходимо подавать на регулятор усредненное значение регулируемого параметра. Для этого в систему введен оператор скользящего среднего.

На рисунке 19 показан результат вычислительного эксперимента для комбинированной системы регулирования.

Рис.19 Изменение значения температуры при использовании комбинированной системы..

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Предложена методика создания АСУ 111 с использованием типового подхода.

2. Разработаны и прошли промышленную апробацию техническое задание, информационное обеспечение, техническое обеспечение, математическое обеспечение технорабочего проекта производства, фенола-ацетона на основе предложенного подхода к созданию АСУ 111.

3. Проведено сравнение различных типов АСР для управления реактором полимеризации с перемешивающим устройством.

4. Проведено исследование влияния запаздывания в измерениях на различные типы АСР.

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

6. Результаты работы используются в промышленном производстве фенола-ацетона и полиэтилена.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Шауро B.C. Софиев А.Э. Моделирование АСР реактором полимеризации с перемешивающим устройством в среде MATLAB // Приборы и автоматизация. -2006г. №7

2. Шауро B.C. Канищев П.А. Краев В.П. Софиев А.Э. Савельев A.M. Система управления установкой производства фенола-ацетона на Уфимском ОАО УФАОРГСИНТЕЗ // Приборы и автоматизация. - 2006г. №6

3. Софиев А.Э. Шауро B.C. Применение вибрационного управления для химико-технологических объектов. // Труды международной конференции «Программные системы : теория и приложения», ИПС РАН, г. Переславль-Залесский, май 2004 / под редакцией С.М. Абрамова. -М.: Физматлит, 2004. -Е.2 стр 475

4. Софиев А.Э. Гончаренко М.В. Шауро B.C. Моделирование трубчатого реактора с учетом баланса импульса. // Сборник трудов 15 международной конференции ММТТ15, Т10 стр 10-13, Тамбов 2002

Принято к исполнению 19/10/2006 Исполнено 19/10/2006

Заказ № 775 Тираж: 100экэ.

Типография «11-й ФОРМАТ» ИНН 7726330900 Москва, Варшавское ш., 36 (495) 975-78-56 www.autoreferat.ru

Введение

Глава 1. Подход к проектированию АСУТП

1.1 Понятие комплексного подхода

1.2 Постановка задачи построения АСУТП с использованием унифицированного подхода

1.3 Примеры применения комплексного подхода в создании АСУТП

Глава 2. Построение АСУТП крупнотоннажным химическим производством на примере производств фенола-ацетона и полиэтилена высокого давления

2.1 Уровень датчиков и исполнительных механизмов

2.2 Уровень сбора информации.

2.3 Уровень обработки и хранения информации

2.4 Уровень человеко-машинного интерфейса

2.5 Техническая структура системы.

2.6 Разработка унифицированного приложения для протоколирования информации для производства фенола-ацетона

Глава 3. Система управления реактором полимеризации.

3.1 Краткие сведения о конструкционных особенностях реакторов полимеризации этилена с перемешивающим устройством

3.2 Механизм и кинетика полимеризации

3.3 Управление реактором полимеризации в неустойчивом режиме

3.4 Применение вибрационного управления к различным системам управления

3.4.1. Линейные системы

3.4.2 Нелинейные системы

3.4.3 Оптимальная Форма Вибрации

3.4.4 Вибрационная стабилизация систем с распределенными параметрами

Глава 4. Применение вибрационного управления для промышленных установок

Глава 5. Применение вибрационного управления для реактора полимеризации этилена с перемешивающим устройством

5.1 Моделирование реактора

5.1.1 Материальный баланс реактора

5.1.2 Тепловой баланс реактора

5.1.3 Исходные данные для моделирования

5.1.4 Исследование устойчивости

5.2 Проведение вычислительных экспериментов

5.3 Влияние колебаний на поведение системы в окрестностях второго стационарного состояния

Глава 6. Сравнение различных систем управления реактором

6.1 Сравнение различных АСР в среде MATLAB

6.2 Реализация комбинированной системы регулирования на средствах программирования контроллеров Siemens . 111 Заключение .116 Библиографический список использованной литературы

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

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

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

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

Сегодня системы диспетчерского управления (Manufacturing Execution System - MES системы) из экзотики превращаются в непременный атрибут любой системы управления. Более высокий уровень автоматизации АСУП или ERP системы (Enterprise Resource Planning), в России пока остается мало распространенным и встречается в основном в нефтеперерабатывающей промышленности и энергетике. Более того, в настоящее время современная техника позволяет реализовывать более высокий уровень автоматизации АСУ Отрасли (On-Line Analytical Processing OLAP системы)

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

Исследование, постановка и решение данных задач проводились на ОАО «Казаньоргсинтез» и ОАО «Уфаоргсинтез». Это крупные химические производства, что позволяет глубоко рассмотреть проблемы решения задач создания АСУТП с помощью унифицированного подхода и задач, решение которых в рамках унифицированного подхода невозможно.

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

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

Научная новизна. В результате проведенной работы были получены следующие основные научные результаты:

1. Впервые разработана комбинированная система управления автоклавным реактором полимеризации на основе вибрационного регулирования и регулирования по обратной связи.

2. Разработана и предложена инженерная методика применения типовых решений при создании АСУТП.

3. Проведено исследование влияния запаздывания в измерениях на различные типы Автоматизированных систем регулирования (АСР).

Практическая ценность и реализация результатов работы.

1. Принят к внедрению проект АСУТП фенола-ацетона разработанный на основе методики применения типовых решений.

2. Разработаны и прошли промышленную апробацию техническое задание, информационное обеспечение, техническое обеспечение, математическое обеспечение технорабочего проекта производства фенола-ацетона на основе предложенного подхода к созданию АСУТП.

В соответствии с поставленной целью работы в первой главе приведена общая характеристика комплексного подхода к построению АСУТП крупного производства. Приводятся примеры существующих АСУТП в различных отраслях промышленности, спроектированных в соответствии с таким подходом.

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

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

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

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

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

Приложения содержат в себе листинг программ созданных для реализации вычислительных экспериментов.

6. Результаты работы используются в промышленном производстве фенола-ацетона и полиэтилена.

Заключение.

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

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

1. Дудников Е.Г., Левин А.А., Промышленные автоматизированные системы управления, М., Энергия, 1973. 192 с.

2. Беляков В.К., Голованов О.В. Автоматизация управления химической промышленностью. М.: Химия, 1977. - 279 с.

3. Голованов О.В., Шкарупа В.В. Разработка программного обеспечения отраслевой АСУ. М.: Статистика, 1978. - 199 с.

4. Борозденков Н.В. Плюсы и минусы комплексного подхода к автоматизации производства / http://www.plcsystems.ru/] (13.04.2005)

5. Власов В.А., Комисарчук С.Ю., Лебедев В.О., Обносов А.В. Системно-ориентированное программирование как средство решения задач реального времени. //Промышленные АСУ и контроллеры. 1999. - №1.

6. Шапиро Ю.З. АСУ химическими производствами. Унифицированные решения. -М.: Химия, 1983. 224 с.

7. Ю.Б. Типовой проект автоматизации технологических процессов на базе технологии Industrial Ethernet //Промышленные АСУ и контроллеры. -2005.-№1-с. 14-22.

8. АСУ ТП литейного отделения №3 ОАО "КрАЗ". / http://www.sinetic.ru/proiects/allum/kraz/lo3/] (12.10.2004).

9. Чичелов В.А Мостовой А.В Бернер Л.И Илюшин С.А. Интегрированная многоуровневая АСУ Предприятия "Пермтрансгаз". //Приборы и Системы. Управление. Контроль. Диагностика. 2002г. - №5 -с. 12-15.

10. Прокопчик В.Ф., Сапожников Г.Н., Терлецкий М.Ю. Опыт внедрения и эксплуатации АСУТП Пинских электрических сетей на базе SCADA-пакета FIX32. // Промышленные АСУ и контроллеры. 2002 - №3 -с. 41-46.

11. Автоматизированная система управления технологическим процессом производства полистиролов ОАО "Нижнекамскнефтехим. / http://www.szma.org/himprom3.plip] (11.07.2004).

12. Онищенко А.Г., Терлецкий М.Ю. Новые решения для эффективного управления технологическими процессами // RM Magazine. 2001.-№2.

13. Дульнева В.В., Кухтин С.Г. Новые возможности SCADA-пакетов FIX и iFIX // Промышленные АСУ и Контроллеры. 2001. №7. - с.36-38.

14. Альперович И.В. iFix "крупноблочное" построение диспетчерских систем АСУТП // PCWeek/RE. - 2001. - №30.

15. Терлецкий М. Ю. Универсальное средство построения производственного архива iHistorian // Промышленные АСУ и Контроллеры. - 2001. - №11 - с. 48-52.

16. Гучек С.И., Вронко Ю.Ф., Гелейша Ю.А., Терлецкий М.Ю. Применение SCADA-пакета iFIX на гродненском производственном объединении "Азот" // Промышленные АСУ и Контроллеры. 2001. - № 9.

17. Капица П.Л. ЖЭТФ 1951, т.21. №5.

18. Боголюбов Н.Н. "О некоторых статистических методах в математической физике" Киев 1945. 137 с.

19. S. М. Meerkov. Averaging of trajectories of slow dynamics system Differential Equation, // IEEE Trans Automat Control, vol. 9, noil, pp. 12391245, 1973.

20. S. M. Meerkov. Principle of Vibrational Control: Theory and Applications // IEEE Trans Automat Control, vol. AC-25, pp 755-762, 1980.

21. R. Bellman, J.Bentsman, S.M. Meerkov. Vibrational control of system with Arrhenius dynamics // J. Math, Anal. Appl. Vol. 90, 1983.

22. S.M. Meerkov. Conditions of vibrational stabilizability for a class of nonlinear systems // IEEE Trans Automat Control, vol. AC-27 pp 485-487, 1982.

23. S.M. Meerkov. Vibrational Stabilizability of Distributed Parameter System // J. Math Anal. Appl vol. 98, 1984.

24. Челомей B.H. О возможности повышения устойчивости упругих систем при помощи вибраций //ДАН СССР, 1956, т.110, №3. с. 345-348.

25. G.F. Wolf. The dynamic stabilization of the Rayleigh-Taylor instability and the corresponding dynamics equilibrium // Z. Phys, vol. 227, pp 291-300, 1969.

26. J.P. Boris. Dynamic Stabilization of the imploding-shell Rayleigh-Taylor instability // Comments Plasma Phys Cont Fusion, vol. 3, pp 1-13, 1977.

27. G.Berge. Equilibrium and stability of MHD-fluids by dynamic techniques // Nuclear Fusion, vol 12, pp 99-117, 1972.

28. P.M. Gresho, R.L. Sani. The effects of gravity modulation on the stability of a heated fluid layer, // J. Fluid Mech vol 40, pp 783-806, 1970.

29. Боголюбов H.H, Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний, -М.: Наука. 1974,- 408 с.

30. S.M. Meerkov. Averaging of trajectories of slow dynamic systems // J. Differential Equations vol. 9, pp 1239-1245, 1973.

31. R. Bellman, J. Bentsman S.M. Meerkov. On vibrational Stabilizability of nonlinear system // J. Optim. Theory ahd Appl. vol.4, 1985.

32. R. Bellman, J. Bentsman S.M. Meerkov. Vibrational control of nonlinear systems: Vibrational Stabilizability // IEEE Transactions on automatic control, vol. AC-11, №8, 1986.

33. R.Bellman, J.Bentsman S.M. Meerkov. Vibrational control of nonlinear systems: Vibrational controllability and transient behavior // IEEE Transactions on automatic control vol. AC-11, №8, 1986.

34. S.M. Meerkov. Condition of Vibrational Stabilizability for a Class of Nonlinear Systems // IEEE Transactions on Automatic Control, vol. AC-27, pp. 485-487, 1982.

35. Демидович Б.П. Лекции по Математической теории устойчивости, -М.: Наука,- 1971.

36. Габасов Р. Кирилова Ф. Качественная теория оптимальных процессов, М.: Наука.- 1971.

37. D. Graupe. Identification of systems // Huntington, NY: Krieger, 1976.

38. S. M. Meerkov and M. Yu. Tsitkin. The effectiveness of the method of vibrational control for the dynamic systems of the order n // Automat. Remote Contr, vol. 36, pp. 525-529, 1975.

39. Beranek Leo L., Ver Istvan L. Noise and vibration control engineering -Principles and applications. John Wiley & Sons, Inc., 814 p. 1992.

40. A. Cinar, J. Deng, S.M. Meerkov, and X. Shu, Vibration Control of an Exothermic Reaction in a CSTR: Theory Simulation, Experiment. // Annual Meeting, San Francisco, CA. Paper No. 3, 1984.

41. A. Cinar, J. Deng, S.M. Meerkov, and X. Shu, Vibrational Control of an Exothermic Reaction in a CSTR: Theory and Experiments. //J. AlChE, vol. 33, pp. 353-365, 1987.

42. J. E. Bailey, and F.M. Horn, Comparisons Between two Sufficient Conditions for Improvement of un Optimal Steady State Process by Periodic Optimization. // J. of Optimization Theory and Applications, vol. 7, №5, 5pp. 378384,1971.

43. S. Bittanti, G. Fronza, and G. Guardabassi, Periodic Control: A Frequency Domain Approach. // IEEE Transactions on automatic control, AC-18, pp. 33-38, 1973.

44. N. Watanabe, К. Onogi, and M. Matsubara. Periodic Control of Continuous Stirred Tank Reactors // Chem. Eng. Sci., vol. 36, pp. 809-818, 1981

45. N. Watanabe, K. Kurimoto, M. Matsubara, and K. Onoqi, Periodic Control of Continuous Stirred Tank Reactors -II. // Chem. Eng. Sci., vol. 37, pp. 745-752, 1982.

46. M. Fjeld. Asynchronous Quenching of Nonlinear Systems Exhibiting Limit Cycles. // IEEE Transactions on automatic control, AC-13, 201-202, 1968

47. M. Fjeld. Relaxed Controls in Asynchronous quenching and Dynamical Optimization. // Chem. Eng.Sci.,vol. 29, pp. 921-933, 1974.

48. D.D. Bruns, and J.E. Bailey. Process Operation Near an Unstable Steady State Using Nonlinear Feedback Control. // Chem. Eng. Sci., vol. 30, pp. 755-762, 1975.

49. D.D. Bruns, and J.E. Bailey. Nonlinear Feedback Control for Operating a Nonisothermal CSTR Near an Unstable Steady Stale.// Chem. Eng. Sci., vol. 32, pp. 257-264, 1977.

50. Кафаров В.В. Софиев А.Э. Трахтенберг A.M. Стабилизация неустойчивых режимов в нелинейных химико-технологических системах с помощью параметрических колебаний. // ДАН СССР, 1983, т. 268, №1 - с 137-142.

51. A. Cinar, J. Deng, S.M. Meerkov, and X. Shu. Vibrational Stabilization of a Chemical Reactor: An Experimental Study // IEEE Trans. Automat. Control, vol. 32, pp.348-352, 1987.

52. Софиев А.Э. Трахтенберг A.M. Геометрическая интерпретация эффекта вибрационного управления // Нестационарные процессы в катализе 3-я всесоюзная конференция. —Т. 1. ■—Новосибирск, 1986. с. 42-43.

53. Дубинина О.В. Дышлевич И.И. Туник О.А. Складановский Е.Н. Ярошевский C.JI. Работа доменной печи с пульсирующем давлением газа на колошнике // Метал и литье Украины. №7. 2001.

54. Трахтенберг A.M. Исследование устойчивости и разработка метода стабилизации неустойчивых режимов в ядерных реакторах типа БВЭР-1000 с помощью параметрических колебаний // Отчет о научно-исследовательской работе ЦНИИКА, Москва, 1984.

55. Софиев А.Э. Трахтенберг A.M. Вибрационная стабилизация неустойчивого режима химического реактора с перемешиванием // ТОХТ, 1989, т. XXII, №5 - с. 644-650.

56. Поляков А.В. Дунтов Ф.И. Софиев А.Э. Полиэтилен высокого давления/ JL: Химия, 1988.

57. Вольтер Б.В. Сальников И.Э. Устойчивость режимов работы химических реакторов. -М.: Химия, 1972 г. 192с.

58. Левич В.Г. Физико-химическая гидродинамика. -М.: ГИФМЛ, 1959 г.- 700с.

59. Кафаров В.В. «Методы кибернетики в химии и химической технологии».-М.: Химия 1968г.

60. Эльсгольц Л.Э. Дифференциальные уравнения. М.: Комкнига, 2006 г.- 312 с.

61. Степанов В.В. Курс дифференциальных уравнений М.: Комкнига, 2006 г. - 472 с.

62. Ляпунов А. М. Общая задача об устойчивости движения. -М.: Гостехиздат, 1950 г. -472 с.

63. Андронов А.А. Витт А.А. Хайкин С.Э. Теория колебаний. М.: Физматгиз, 1959 г. 916 с.

64. Айзерман М.А. Теория автоматического регулирования. М.: Наука, 1966 г.-427 с.

65. Солодовников В.В., Плотников В.Н., Яковлев А.В. Основы теории и элементы систем автоматического регулирования: Учеб. пособие.- М.: Машиностроение, 1985 г.

66. Софиев А.Э. Софиева Ю.Н. Теория управления М.: МИХМ 1975 г. -166 с.

67. С. Van Heerden Autothermic Processes: Properties and Reactor Design. // Ind. Eng. Chem., vol.45, pp. 1242, 1953.

68. Андронов А.А., Леонтович E.A., Гордон И.И., Майер А.Г. Теория бифуркаций динамических систем на плоскости. М.: Наука, 1967 г. - 488 с.

69. Авторское свидетельство. Способ управления реактором полимеризации этилена. А.С. 1159927 СССР.- С08 F 110/02/ А.Э. Софиев, A.M. Трахтенберг. -J 3643037, заявлено 19.09.85.