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

Ип пппапг гтк описи

003456317

Лежнев Максим Владимирович

СТРУКТУРНОЕ МОДЕЛИРОВАНИЕ И АВТОМАТИЧЕСКОЕ УПРАВЛЕНИЕ ТЕМПЕРАТУРОЙ

АБСОРБЕНТА В ТЕПЛООБМЕННОМ АППАРАТЕ УСТАНОВКИ КОМПЛЕКСНОЙ ПОДГОТОВКИ ГАЗА

Специальность 05.13.06 - Автоматизация и управление технологическими процессами и производствами (промышленность)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Самара-2008

003456317

Работа выполнена на кафедре «Автоматика и управление в технических системах» ГОУ ВПО «Самарский государственный технический университет»

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

- Заслуженный деятель науки и техники РФ, доктор технических наук, профессор Рапопорт Эдгар Яковлевич

Официальные оппоненты: - доктор технических наук, профессор

Абакумов Александр Михайлович

- кандидат технических наук, доцент Масляницын Александр Петрович

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

ГОУ ВПО "Саратовский государственный технический университет", г. Саратов

Защита диссертации состоится «24» декабря 2008 г. в 11 часов на заседании диссертационного совета Д 212.217.03 в ГОУ ВПО "Самарский государственный технический университет" по адресу: 443010, г. Самара, ул. Галактионовская, 141, ауд. 28.

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО "Самарский государственный технический университет" по адресу: 443100, г. Самара, ул. Первомайская, 18, корп. №1 и на официальном сайте www.samgtu.ru.

Отзывы на автореферат в двух экземплярах, заверенные печатью просим направлять по адресу: 443100, г. Самара, ул. Молодогвардейская 244, СамГТУ, Главный корпус, ученому секретарю диссертационного совета Д

212.217.03.

Автореферат разослан ноября 2008 г.

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

Н.Г. Губанов

ОБЩАЯ ХАРАКТЕРИСТИКА ДИССЕРТАЦИОННОЙ

РАБОТЫ

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

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

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

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

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

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

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

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

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

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

3. Компьютерное моделирование и исследование динамики системы управления температурным нолем противоточного теплообменного аппарата в пакете компьютерного моделирования сосредоточенных динамических систем.

4. Структурно-параметрический синтез и промышленное внедрение системы автоматического управления температурой абсорбента в теплообменном аппарате установки комплексной подготовки газа.

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

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

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

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

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

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

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

1. Разработано математическое, алгоритмическое, программное обеспечение и инженерная методика структурно-параметрического синтеза систем автоматического управления противоточным теплообменным аппаратом установки комплексной подготовки газа.

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

3. Методом компьютерного моделирования обоснована эффективность предложенной структуры системы автоматического управления.

Реализация результатов работы. Результаты работы использованы в проектных разработках и при промышленном внедрении перспективных систем управления технологическими процессами в ООО Научно-Внедренческая Фирма «СМС» (г. Новокуйбышевск).

Апробация работы. Основные положения и результаты работы докладывались и обсуждались на международных и всероссийских конференциях: Тринадцатой международной научно-технической конференции студентов и аспирантов «Радиоэлектроника, электротехника и энергетика» (Москва, 2007); Всероссийской научной конференции с международным участием «Математическое моделирование и краевые задачи» (Самара 2007, 2008); XIII Международная научно-практическая конференция студентов, аспирантов и молодых ученых «Современные техника и технологии» (Томск, 2007).

Публикации. По теме диссертации опубликовано 6 печатных работ, в том числе 2 публикации в журналах из перечня, рекомендованного ВАК.

Структура и объем работы. Диссертация состоит из введения, 5 глав и заключения, изложенных на 110 страницах машинописного текста, содержит 45 рисунков, 3 таблицы, список литературы го 81 наименования и 2 приложения.

На защиту выносятся следующие положения:

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

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

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

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

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

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

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

В настоящее время существуют три широко распространенные технологии осушки природного газа, при подготовке его к транспортировке:

- низкотемпературная сепарация; .

- адсорбция;

- абсорбция.

В России наиболее широкое распространение получила технология абсорбционной осушки газа.

Абсорбционные установки применяются для осушки природного газа до точки росы -15 -30 °С. В качестве абсорбента используют в основном растворы диэтиленгликоля (ДЭГ) и триэтиленгликоля (ТЭГ). Степень осушки газа на абсорбционных установках определяется главным образом концентрацией подаваемого в абсорбер раствора, а концентрация раствора, в свою очередь, зависит от используемого на установке метода регенерации отработанного абсорбента. Для глубокой регенерации раствора и получения низких (-20 -30°С) точек росы осушенного газа регенерацию диэтиленгликоля проводят под вакуумом.

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

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

Рис. 1 - Технологическая схема абсорбционной установки осушки газа

Осушенный природный газ со всех установок поступает в промысловый газосборный коллектор 5.

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

Установка регенерации представляет собой систему взаимосвязанных технологических аппаратов тепловой подготовки и массообмена, в состав которой входят: теплообменник 8, регенерационная колонна 9, испаритель 10, насос откачки регенерированного раствора 11, холодильники 12 и 15 (воздушные или водяные), сборник флегмы 16, насос откачки флегмы 17 и вакуум-насос 18.

Абсорбент регенерируется в регенерационной тарельчатой колонне 9 под вакуумом, который создается в результате откачки из нее паров вакуум-

5

насосом 18. Нагретый в теплообменнике 8 насыщенный раствор ДЭГ поступает на тарелку питания. Паровая фаза нагретого насыщенного раствора суммируется с паром, поступающим из испарителя, проходит через тарелки верхней (концентрационной) секции колонны и выводится из нее в холодильник 15.

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

Поток насыщенного ДЭГ смешивается на тарелке питания с жидкостью, поступающей из верхней части колонны, и стекает по тарелкам нижней (отгонной) секции колонны в испаритель 10. Здесь вода окончательно выпаривается. В качестве теплоносителя используется пар (иногда применяется огневой подогрев природным газом). Из испарителя паровая фаза направляется в колонну, а регенерированный раствор, нагретый до температуры кипения, насосом 11 через теплообменник 8 и холодильник 12 подается в емкость 13.

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

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

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

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

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

(1)

= о<*</

(2)

д1 1 дх с граничными условиями на левом конце

и начальным условием

№<» = &,(*) (3)

Здесь 0,(а',/), в зависимости от рассматриваемого процесса, имеет смысл температуры или концентрации, изменяющейся в зависимости от времени г и пространственной координаты л: в направлении движения управляемого потока со скоростью /¡(х,/) - функция, описывающая пространственное и временное распределение источников тепла или вещества. Коэффициент а1 имеет смысл коэффициента теплообмена нагреваемого агента с внешней средой. Для уравнения массопереноса коэффициент а, принимается равным нулю в виду отсутствия потерь вещества. В дальнейшем исследовалась модель более общая модель процесса теплообмена, которая в частном случае при а = 0 описывает процессы массопереноса.

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

а, 0(х,1)

Нагреваемый агент

0,(0,/) Греющий

0,(0,0

Рис. 2 - Схема противоточпого теплообменного аппарата

Температурное распределение нагреваемого агента описывается уравнениями (1)—(3). Аналогичные уравнения, записанные для греющего агента, будут иметь вид:

= О**,/

с4 дх

(4)

62<Х0) = &>(*)■. (6)

При конвективном характере теплообмена между нагреваемым и греющим агентами, связи между ними описываются соотношениями:

/дх.0=в„Сасх,о

Структурная схема противоточного теплообменного аппарата, описываемого системой уравнений (1)-(7) и полученная с помощью структурной теории распределенных систем, приведена на рис.3., где в роли входных сигналов фигурируют внешние воздействия g^(p), g2(р) и О0.(х), (2ю(х) по граничным и начальным условиям в (2), (3) и (5),(6), а в качестве управляемых величин рассматриваются температуры (^¡{ир) и Я2(0,р) нагреваемого и греющего агентов на выходе го теплообменника.

Рис. 3 - Структурная схема объекта управления

Здесь и №2(х,£,р) - известные передаточные функции (изо-

бражения по Лапласу функций Грина) типовых распределенных блоков, описываемых уравнениями (1)-(3) и (4)-(6) соответственно:

1 ——-(т>+1>и-;> 1Г,(х,£,р) = , V, > О

1 —, ' .(Г,р-И)(,'чО

Щх,£,р) = 1(£-х) — е'™ ,у2>0 v2

Где:

7>—, Г3=—; (9)

а1 + агх аг+аа

Е, и х -пространственные аргументы входных воздействий и выходных переменных блоков; р - комплексная переменная; - единичная функция переменной х с запаздыванием £; 5(г) - дельта - функции пространственного аргумента г, в терминах которых формируются передаточные функции переходных блоков структурной схемы; ®х{х,р) , сог{х,р) - изображения частей стандартизирующих функций, определяемых начальными и граничными условиями в (2), (3) и (5), (6). Структурная схема рис. 3 приводит к интегральным уравнениям Фредгольма 2-го рода относительно передаточных функций объекта от управлений £,(/?), 82{р) по граничным условиям к каждому из выходов и (2г(0,р) ■ В диссертации получены точные решения этих уравнений в следующей форме:

V,

(10)

I ауаЛ-е^'^-ЬХрУ-ЬЛрУУе™ V, V,у2 (6,(Р) + Й2(7))"' + аааи{1 - е-*1'™» - Ь, (р)1 - Ь2 (/,)/)

Шг1 (р) = ШШ. =-^^- . [] - е-МР)<А(,У ]+

+_Е»_х • (11)

(1 — е"*-"'»'-^1''»' _ А. - (/?>/)-в"*-«'" V, ^ (6, (р) + Ь2 (р))2 + аааг1 (I - - 6, О?)/ - 62 (р)Г)

&0>) ЧъЫрУ + ЬгСр)) 1 ]

(12)

У,У,(6,00 + 6,0?))

а31а,Д2/(ВД + Чр)} + е- О

где:

Ь1{р) = —{р + ах +аи), Ь2(р) = —{р + аг + ап).

V, V,

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

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

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

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

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

Для вывода выражений, описывающих распределение температурного поля во времени по длине теплообменного аппарата, используются дифференциальные уравнения (1)-(6) с заменой соотношений (7) на равенства:

/(х,о=(х,о -е, (*,о] (м)

/2(*,о=«»ЬГ(*>')-а(*,о]- • (15)

Здесь 0;4Ш'(х,1) и О'^Ш) - аппроксимированные функции пространственного распределения температур в нагреваемом и греющем потоках соответственно. В условиях возможности реализации температурного контроля

12

только на входе и выходе теплообменника используется соответствующая двухпараметрическая аппроксимация зависимостей и Qf{x,t) в

форме линейной функции пространственной координаты в целях последующей параметрической идентификации этих величин по наблюдаемым переменным:

0;1п\х,1) = к«)х + к&) (16)

= *,(*)* +*4(0. (17)

Здесь А,(0, кг{О, А,(0> коэффициенты аппроксимации. При фиксированных значениях параметров аппроксимации в (16), (17) краевые задачи (1)-(3), (14), (16) и (4)-(6), (15), (17) для нагреваемой и греющей сред могут быть решены независимо друг от друга. Соответствующие выражения для изображений управляемых функций состояния, описывающие пространственно-временное распределение температурок полей (х, р) и Q1 (х, р) в зависимости от измеряемых граничных воздействий g¡(p), Ег(р) и параметров аппроксимации А,, к2, ку! /с, при нулевых начальных условиях найдены в диссертации в следующем виде:

<2Х*,р)=011(х,р)+д11(х,р)+оа(х,р);

где:

йг(х> Р) = <2гх Р) + вп (х> Р) + й>з (х- Р) .

в„(х,р) = Иг„(х,р)

- —-—Лт,р* Ое-о 1Г11(х,р) = е«™1 ■

( (Г|р+1Х/-х)

(Г^х-ч А

(18)

(19)

(20) (21)

вп(х,Р) = ИГп(х,р)-к1(р)

■аЛа,+0Си)

КАХ>Р)

I

I', (<г,+ап) 4

(22)

Qn{x,p) = Wn{x,p)-k1(p)-

an(a2+ccn)

g "lCjfOu) _ g »i(<«j+«4»)

Wn{x,p) =

ftp + l) Q„(x,p) = fFa(x,p)-kt{p);

f (r„tl)(.t)

al2(a2+al2)

-(гд/4-ixz-.q >

(25)

(Г2р + 1)

где Г; и Г2 - определяются согласно (9).

Уравнения (18) и (19), полученные для модели противоточного теплооб-менного аппарата, описывают поведение управляемых величин при постоянных скоростях течения обоих потоков. Однако, в большинстве случаев, основным возмущающим воздействием является именно изменение этих скоростей.

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

гас*,*)

dt

dQ2(x, Q

dt

+ v,(í)

дх

■(a,+a2¡)Ql(x,t) = a2lQ2(x,t);

- v> w dQltJ)+(«i+fe ')=a«a o •

&C

(26)

(27)

Переход к линеаризованной системе уравнений в приращениях относительно установившегося режима работы объекта с постоянными скоростями потоков Vio и \го'.

3AQM dt

8Щ(х;/)

+ А v,(í)

90о(*) , . SAQ(a,í)

-Av2(Í)

дх

¿Шх)

- + v,

&

дЮг(х,Г)

+ (at + а2) )Ag, (х, í) = 0; (28)

+ (а1+а11)Да(х,0 = 0, (29)

Зс ™. дх

и последующее решение этой системы в изображениях по Лапласу относительно Д<2: и Лб2 приводит к следующим результатам:

Ав1(х,р) = !¥и(х,р)-Ау1(р); (30)

( г (Г,р+1) Л

1-е

V

Л,

+ — Р

1 { (Г,я+1) , .1

К V -Г ' —

и Ггт */» 1 1 ы

у\

А()1(х,р) = 1Уи(х,р)-Ауг(р);

...(я.+о,,)

е ,я 1 "- е

1

<7*20 (Р)

А

Р

/ « \ е ""(<г,"'"> - е

(Т2р +1)

с Пя*' Л

V У

г

(31)

р(«2+а|2)

\\

Л/

где к}, к6, к^, - постоянные, определяемые по характеристикам номинального режима работы теплообменного аппарата. На рис. 4 представлена структурная схема объекта управления, построенная по выражениям (20)-(25), (30), (31), где IV,,, IV,2, 1У13, \У14, 1У„, 1¥23, \У21 - передаточные функции распределенных х - блоков, имеющих сосредоточенный вход и распределенный выход.

ГТТ:—Л

*,0>)

ш -

Щр)

Рис. 4 - Структурная схема объекта управления

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

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

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

Рис. 5 - Экспериментальные данные, использованные для идентификации объекта

управления

На рис. 6. представлена соответствующая рис. 4 структурная схема, используемая для идентификации объекта управления.

Зужсп(р)

Лр)

сп(р)

г~Н8Н 1/1

81-зксп(Р)

Дуг

Чр)

Ы

М }-

йэкга ('•/>)

^Ь-жсп(Р)

1/1

к2р)

^(М

&(0,Р)

м -I

1/р

Ар)

йгжсп^Р)

Рис. 6 - Структурная схема, используемая для параметрической идентификации объекта управления.

Здесь: - коэффициенты зависимости скорости от степени открытия

регулирующего клапана; Спкс,/р), Огжс/р) - степени открытия регулирую-

4 4

щих клапанов в процентах; (1,р)(1,Р); й7^ (О,/?) .

Коэффициенты аппроксимации в (16) и (17) рассчитываются по измеряемым температурам на входе и выходе теплообменника на каждом шаге итерационной процедуры процесса по простейшим зависимостям:

*г(0=у(еД0-Й(0,0); ^(0 = 0(0,0 -

(32)

(33)

В качестве параметров идентификации выбираются постоянные времени 7) и Т2 в выражениях (20)-(25), (30), (31), и скорости у^,\-2о номинального режима работы аппарата, поиск оптимальных значений которых, из условий минимизации интегральной величины суммарной абсолютной ошибки

J

->пип;

производится в пакете структурного моделирования ЭгшиНпк в составе пакета МАТЬАВ.

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

Таблица 1. Конструктивные параметры теплообмениого аппарата.

Номинальный расход 14 м3/ч

Общая длина рабочей области 12 м

Объем в трубной части 6м3

Объем в корпусной части 11.7 м3

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

Таблица 2. Расчетные и найденные путем параметрической идентификации па____раметры теплообмениого аппарата.__

Наименование параметра Рассчитанное по характеристикам значение Значение, полученное по результатам идентификации

Постоянная времени (Т¡) 501 с 630.72 с

Постоянная времени (Тг) 503 с 629.07 с

Номинальная скорость потока (у,0) 0.0078065 м/с 0.0062452 м/с

Номинальная скорость потока (у ,„) 0.00400068 м/с 0.0033339 м/с

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

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

В качестве регулятора был выбран ПИ-регулятор, широко используемый в инженерной практике для управления объектами подобного рода.

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

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

3 =

=

гг —» Ш1П .

(34)

Для параметрической оптимизации регулятора использовался пакет структурного моделирования БтшПпк, в котором реализовывалась структура системы автоматического управления объектом, представленная на рис. 7, где !-](р) и Р\(р) - возмущения по температурам на входах теплообменника; К3 -коэффициент зависимости скорости от степени открытия клапана; 0(р)- степень открытия регулирующего клапана.

я,(я)_

ОиЛр) 1-

Ш

к,

Мр)

^ (Ш

ОАКр)

1 //

: К(Р)

\ (0,р)

0г(0,р)

глр)

Рис. 7 - Структурная схема системы автоматического управления температурой на выходе теплообменного аппарата.

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

Т.иР .

ад:

Т.ПР +1

Т,пР + 1

параметры которых, наряду с параметрами ПИ-регулятора, находятся в процессе минимизации критерия (34). На рис. 8 представлена структурная схема системы управления температурой на выходе теплообменного аппарата с компенсаторами, отрабатывающими возмущения /\(р) и Ег{р).

Рис. 8 - Структурная схема системы управления с двумя компенсирующими звеньями.

Переходный процесс изменения температуры нагреваемого потока на выходе теплообменного аппарата в комбинированной системе управления при воздействии 10% возмущений Г, и относительно входных температур представлен на рис. 9.

О 05 1

4

х10 .С

Рис. 9 20

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

Для практической реализации алгоритмов управления противоточным тепдообменным аппаратом был использован программируемый логический контроллер фирмы Siemens серии S7-417H. Данная серия обладает высоким быстродействием (несколько млн. инструкций в секунду), а кроме того является программируемым логическим контроллером повышенной надежности за счет использования аппаратного резервирования ПЛК и шин передачи данных.

Каждый вычислительный блок реализуется в виде отдельного функционального блока (Function Block) программы контроллера. Это позволяет модифицировать алгоритм любого блока, не изменяя остальные блоки. Большой объем памяти контроллера и его быстродействие позволяют реализовать на контроллере не только систему автоматического управления противоточным теплообмешшм аппаратом, но и управление всеми системами, участвующими в процессе абсорбционной осушки газа: систему управления абсорберами; систему управления процессом регенерации; систему контроля загазованности в цехах; систему блокировок и противоаварийных защит; а также систему экстренного останова производства осушки, связанную с другими системами управления установки комплексной подготовки газа через резервированную оптическую сеть Industrial Ethernet.

Супервизорное управление процессом комплексной подготовки газа осуществляется с нескольких автоматизированных рабочих мест оператора, реализованных на промышленных ЭВМ, причем для обеспечения непрерывности наблюдения и управления за технологическим процессом используется клиент-серверная архитектура с резервированными серверами, обеспечивающая повышенную надежность системы управления. Станции операторов и сервера работают под управлением операционной системы Microsoft Windows NT 4.0 и SCADA-системы Siemens SIMATIC PCS7 v5.2 Service Pack 3. Мнемосхемы SCADA-системы WinCC обеспечивают удобный, интуитивно понятный пользовательский интерфейс позволяющий в процессе работы отслеживать изменения всех контролируемых параметров, изменять параметры регуляторов и различные настройки алгоритмов управления.

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

1. Методами структурной теории распределенных систем разработана структурная схема противоточного теплообменного аппарата как объекта автоматического управления. Путем решения интегральных уравнений

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

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

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

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

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

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

1. Данилушкин И.А., Лежнев М.В. Разработка математической модели процесса абсорбционной осушки газа.// Вестник Самарского государственного технического университета. Сер. "Технические науки". Самара, 2006. Вып. 41. С. 184-188

2. Данилушкин И.А., Лежнев М.В. Структурное представление процесса теплообмена при встречном направлении взаимодействующих потоков.// Вестник Самарского государственного технического университета. Сер. "Технические науки". Самара, 2007. Вып. №1(19). С. 16-22.

3. Лежнев М.В. Построение математической модели процесса абсорбционной осушки газа.//Тринадцатая международная научно-техническая конференция студентов и аспирантов "Радиоэлектроника, электротехника и энергетика", 1-2 марта 2007г. : Тезисы докладов. В 3-х т. -М: Издательский дом МЭИ, 2007 - Т. 1., с.340-341.

4. Лежнев М.В. Разработка упрощенной модели противоточного теплообменника как объекта управления с распределенными параметрами.// XIII Международная научно-практическая конференция студентов, аспирантов и молодых ученых «Современные техника и технологии», 26 -30 марта 2007 г. Труды в 3-х т. - Томск: Изд-во ТПУ, 2007 - Т.2., с.375-377.

5. Лежнев М.В., Данилушкин И.А. Моделирование процессов тепло-массопереноса с использованием структурной теории распределенных систем.//Труды четвертой Всероссийской научной конференции с международным участием "Математическое моделирование и краевые задачи", 29-31 мая 2007г.: Ч. 2: Моделирование и оптимизация динамических систем и систем с распределенными параметрами. -Самара: СамГТУ, 2007.,

6. Лежнев М.В., Рапопорт Э.Я., Данилушкин И.А. Численное моделирование процессов теплопереноса в противоточном теплообменном аппарате. // Труды пятой Всероссийской научной конференции с международным участием "Математическое моделирование и краевые задачи", 29-31 мая 2008г.: Ч. 2: Моделирование и оптимизация динамических систем и систем с распределенными параметрами. -Самара: СамГТУ, 2008., с.66-69.

Разрешено к печати диссертационным советом Д 212.217.03 ГОУВПО «Самарский государственный технический университет» (протокол № 9 от 14.11. 2008)

Заказ № 773. Формат 60x84 1/16. Усл. печ. л. 1. Тираж 100 экз. Отпечатано на ризографе. ГОУ ВПО «Самарский государственный технический университет» Отдел типографии и оперативной печати 443100, г. Самара, ул. Молодогвардейская, 244.

с.90-92.

Содержание.

Введение.

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

1.1 Современные технологии осушки природного газа.

1.2 Современное состояние проблемы управления процессом теплообмена как сложной физической системой.

Выводы по первой главе.

2. Структурное моделирование противоточного теплообменпого аппарата как объекта управления с распределенными параметрами.

2.1 Базовая аналитическая модель температурного поля противоточного теплообменника.

2.2 Передаточные функции теплообменника.

Выводы по второй главе.

3. Приближенное моделирование объекта управления.

3.1 Структурное представление объекта управления с линеаризованным пространственным распределением взаимодействующих источников тепла.

3.2 Структурное представление объекта управления с переменными скоростями движения встречных потоков в теплообменном аппарате

Выводы по третьей главе.

4. Идентификация объекта управления.

4.1 Исходные данные для проведения параметрической идентификации объекта управления.

4.2. Параметрическая идентификация объекта управления.

Выводы по четвертой главе.

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

5.1 Параметрическая оптимизация настроек ПИ-регулятора.

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

Выводы по пятой главе.

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

Актуальность работы.

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

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

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

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

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

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

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

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

3. Компьютерное моделирование и исследование динамики системы управления температурным полем противоточного теплообменного аппарата в пакете компьютерного моделирования сосредоточенных динамических систем.

4. Структурно-параметрический синтез и промышленное внедрение системы автоматического управления температурой абсорбента в теплообменном аппарате установки комплексной подготовки газа.

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

Научная новизна.

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

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

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

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

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

Практическая полезность работы. Прикладная значимость проведенных исследовании определяется следующими результатами:

1. Разработано математическое, алгоритмическое, программное обеспечение и инженерная методика структурно-параметрического синтеза систем автоматического управления противоточным теплообменным аппаратом установки комплексной подготовки газа.

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

3. Методом компьютерного моделирования обоснована эффективность предложенной структуры системы автоматического управления.

Реализация результатов работы. Результаты работы использованы в проектных разработках и при промышленном внедрении перспективных систем управления технологическими процессами в ООО Научно-Внедренческая Фирма «СМС» (г. Новокуйбышевск). Апробация работы.

Основные положения и результаты работы докладывались и обсуждались на международных и всероссийских конференциях: Тринадцатой международной научно-технической конференции студентов и аспирантов «Радиоэлектроника, электротехника и энергетика» (Москва, 2007); Всероссийской научной конференции с международным участием «Математическое моделирование и краевые задачи» (Самара 2007, 2008); XIII Международная научно-практическая конференция студентов, аспирантов и молодых ученых «Современные техника и технологии» (Томск, 2007).

Публикации. По теме диссертации опубликовано 6 печатных работ, в том числе 2 в журналах из перечня, рекомендованного ВАК.

Структура и объем работы. Диссертация состоит из введения, 5 глав и заключения, изложенных на 110 страницах машинописного текста, содержит 45 рисунков, 3 таблицы, список литературы из 81 наименования и 2 приложения.

Выводы по пятой главе

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

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

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

Заключение

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

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

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

3. Методами структурной теории распределенных систем разработана структурная схема противоточного теплообменного аппарата как объекта автоматического управления.

4. Выполнено преобразование структурной схемы, обеспечившее нахождение передаточных функций независимо по каждому из каналов (вход-выход).

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

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

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

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

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

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

11. Синтезирована численная модель системы управления, обеспечивающая приемлемые результаты управления противоточным теплообменным аппаратом.

12. Проведена параметрическая оптимизация ПИ-регулятора системы управления температурным полем теплообменника.

13. Выполнен параметрический синтез системы управления с компенсирующими звеньями

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

1. Тараненко Ф.Б., Герман В.Т. Автоматическое управление газопромысловыми объектами. М.: Недра, 1976.

2. Коротаев Ю.П. и др. Добыча, подготовка и транспорт природного газа и конденсата. Справочное руководство в 2-х т. М.: Недра, 1984.

3. Корчажкин М. Т., Технологическая схема сепарации высоконапорного конденсатного газа, в кн.: Добыча газа М., 1961.

4. Базлов М. Н., Жуков А. И., Алексеев Т. С., Подготовка природного газа и конденсата к транспорту. —М., 1968.

5. Фролов Ю. Г. Курс коллоидной химии. Поверхностные явления и дисперсные системы. М.: Химия, 1989.

6. Кельцев Н. В. Основы адсорбционной техники. М.: Химия, 1984.

7. Грег С., Синг К. Адсорбция, удельная поверхность, пористость. М.: Мир, 1984.

8. Коротаев Ю.П., Тагиев В.Г., Самородкин В.Д. Оптимизация режимов эксплуатации объектов добычи природного газа. М.: Недра, 1982.

9. Репин Н.Н., Тагиев В.Г. Оптимальное управление установками комплексной подготовки природного газа. -М.:Недра, 1992.

10. Гриценко А.И., Истомин В.А, Кульков А.Н., Сулейманов Р.С. Сбор и промысловая подготовка газа на северных месторождениях России . -М.: Недра, 1999.

11. Гвоздев Б.П., Гриценко А.И., Корнилов А.Е. Эксплуатация газовых и газоконденсатных месторождений: Спр. пос. М.:Недра, 1988.

12. Михеев М.А., Михеева И.М. Основы теплопередачи. М.: Энергия, 1979.

13. Кутателадзе С.С. Основы теории теплообмена. М.: Наука, 1970.

14. Кафаров В.В., Мешалкин В.П., Гурьева J1.B. Оптимизация теплообменных процессов и систем. М.:Энергоатомиздат, 1988.

15. Манысовский О.Н., Толчинский А.Р., Александров A.M. Теплообменная аппаратура химических производств. Л.:Химия, 1976.

16. Романков П.Г., Фролов В.Ф. Теплообменные процессы химической технологии. Л.: Химия, 1982.

17. Миртов Н.Н., Белякова И.Г. Теплоотдача при конденсации пара в термосифоне//ИФЖ. 1982. Т.43. №3.

18. Карслоу Г., Егер Д. Теплопроводность твердых тел. М.:Наука, 1964.

19. Эрриот П. Регулирование производственных процессов: Пер. с англ. М.: Энергия, 1967.

20. Полоцкий Л.М., Лапшенков Г.И. Автоматизация химических производств. Теория, расчет и проектирование систем автоматизации. М.: Химия, 1982.

21. Казаков А.В. Кулаков М.В., Мелюшев Ю.К., Основы автоматики и автоматизации химических производств. М.: Машиностроение, 1970.

22. Алифанов О.М. Обратные задачи теплообмена. -М.Машиностроение, 1988.

23. Пантакар С. Численные методы решения задач теплообмена и динамики жидкости. — М.: Энергоатомиздат, 1984.

24. Матвейкин В.Г., Муромцев Д-Ю. Теоретические основы энергосберегающего управления динамическими режимами установок производственно-технического назначения. М.: "Издательство Машиностроение-1", 2007.

25. Рэй, Д. Экономия энергии в промышленности / Д. Рэй ; пер. с англ. -М„ 1985.

26. Аджиев, М.Э. Энергосберегающие технологии / М.Э. Аджиев. М., 1990.

27. Манысовский О.Н., Толчинский А.Р., Александров М.В. Теплообменная аппаратура химических производств. Л.: Химия, 1976.

28. Мигай В.К. Повышение эффективности современных теплообменников. Л.: Энергия, 1980.

29. Бажан П.И., Каневец Г.Е., Селиверстов В.М. Справочник по теплообменным аппаратам. М.: Машиностроение, 1989.

30. Mozley J.M. Predicting Dynamics of Concentric Pipe Heat Excangers // Ind. Eng Chem. 1956. Vol.48. No.6. P. 1035-1041.

31. Stremole F.J., Larson M.A. The Dynamics of Flow Forced Distributed Parameter Heat Exchangers // AIChE Journal. 1964. 10. No. 5.

32. Девятов Б.Н., Демиденко Н.Д., Охорзин B.A. Динамика распределенных процессов в технологических аппаратах, распределенное конструирование и управление. Красноярск.'.Красноярское книжное издательство, 1976.

33. Девятов Б.Н, Демиденко Н.Д. Теория и методы анализа управляемых распределенных процессов. — Новосибирск.:Наука, 1983.

34. Дорри М.Х., Поляк Г.Л. Математическое моделирование теплообменников при переменном расходе теплоносителя.//Исследование тепло и массообмена в технологических процессах и аппаратах. Под ред Лыкова А.В. — Минск.:Наука и техника, 1966.

35. Бутковский А.Г. Структурная теория распределенных систем. М., Наука, 1977.

36. Рапопорт Э.Я. Анализ и синтез систем автоматического управления с распределенными параметрами: Учеб. пособие. М.: Высш. шк., 2005.

37. Рапопорт Э.Я. Структурное моделирование объектов и систем управления с распределенными параметрами. М.:Высш. шк., 2003.

38. Бутковский А.Г. Методы управления системами управления системами с распределенными параметрами. М.:Наука., 1975.

39. Першин И.М. Синтез систем с распределенными параметрами. — Пятигорск, 2002.

40. Шевяков А. А., Яковлева Р. В. Управление тепловыми объектами с распределенными параметрами. — М. Энергоатомиздат, 1986.

41. Девятов Б.Н. Теория переходных процессов в технологических аппаратах с точки зрения задач управления. Новосибирск.: АН СССР, 1964.

42. Кэмпбелл Д.П. Динамика процессов химической технологии. -М.:Госхимиздат, 1962.

43. Сиразетдинов Т.К. Устойчивость систем с распределенными параметрами. — Казань. :КАИ, 1971.

44. Шевяков А. А., Яковлева Р. В. Инженерные методы расчета теплообменных аппаратов М.: Машиностроение, 1968.

45. Голицейский Б.М., Дрейцер Г.А., Изосимов В.Г. и др. М.: Машиностроение, 1973.

46. Кудряшев Л.И., Меньших H.JT. Приближенные методы решения нелинейных задач теплопроводности. М.Машиностроение, 1979.

47. Гудмэн Т. Применение интегральных методов в нелинейных задачах нестационарного теплообмена.//Проблемы теплообмена. — М.:Атомиздат, 1967.

48. А.В. Лыков, Ю.А. Михайлов Теория тепло- и массопереноса. М.-Л., Государственное энергетическое издательство, 1963.

49. Лыков А.В. Теория теплопроводности. М.:Высшая школа, 1967.

50. Кошляков Н.С. Уравнения в частных производных математической физики. -М.:Высшая школа, 1970.

51. Бутковский А.Г. Характеристики систем с распределенными параметрами (справочное пособие). — М.: Наука, 1979.

52. Данилушкин И.А., Лежнев М.В. Разработка математической модели процесса абсорбционной осушки газа.// Вестник Самарского государственного технического университета. Сер. "Технические науки". Самара, 2006. Вып. 41.

53. Манжиров А.В., Полянин А.Д. Справочник по интегральным уравнениям: Методы решения. — М.: Изд-во «Факториал Пресс», 2000.

54. Данилушкин И.А., Лежнев М.В. Структурное представление процесса теплообмена при встречном направлении взаимодействующих потоков.// Вестник Самарского государственного технического университета. Сер. "Технические науки". Самара, 2007. Вып. №1(19).

55. Амелысин С.А., Цирлин A.M. Предельные возможности теплообменников при различных моделях потоков теплоносителей.//Теплоэнергетика. 2001. №5.

56. Рэй У. Методы управления технологическими процессами. Пер. с англ. —М.:Мир, 1982.

57. Янг В. Дж., Кларк И.А., Арпаси B.C. Труды Американского общества инженеров-механиков, пер. с англ. М.:Теплопередача, №3, 1961.

58. Эннс М. Труды Американского общества инженеров-механиков, пер. с англ. М.:Теплопередача, №4, 1963.

59. Рапопорт Э.Я., Данилушкин А.И. Алгоритмы функционирования процесса непрерывно-последовательного индукционного нагрева. //Алгоритмизация и автоматизация технологических процессов и промышленных установок. Куйбышев.: Изд. КПтИ, 1976, вып.7.

60. Полянин А.Д. Справочник по линейным уравнениям математической физики. -М.: ФИЗМАТ ЛИТ, 2001.

61. Бугров Я.С. Высшая математика: Учеб для вузов: В 3 т./Я.С. Бугров, С.М. Никольский; Под ред. В.А. Садовничего. -5-е изд., стереотип. — М.:Дрофа, 2003.

62. Бусленко Н.П. Математическое моделирование производственных процессов. М.:Наука, 1964.

63. Сейдж Э.П., Мелса Д.Л. Идентификация объектов управления. — М.: Наука, 1974.

64. Дьяконов В.П. MATLAB 6: Учебный курс. -СПб.: Питер, 2001. — 592 с.

65. Бакластов A.M., Горбатенко В.А., Данилов O.JI. и др. Под ред Бакластова A.M. Промышленные теплообменные процессы и установки. Учебник для вузов. М.: Энергоатомиздат, 1986

66. Гуров A.M., Починкин С.М. Автоматизация технологических процессов. Учебное пособие. М.:Высш.шк. 1979.

67. Беннет К.О., Мойерс Дж. Е. Гидродинамика, теплообмен и массообмен. пер. с англ. М.: Недра, 1966.

68. Борисов Г. С., Брыков В.П., Дытперский Ю.И. Основные процессы и аппараты химической технологии: Пособие по проектированию 2-е изд. М.: Химия, 1991.

69. Эйкхофф П. Основы идентификации систем управления. — М.:Наука, 1975.

70. Молчанов А.А. Моделирование и проектирование сложных систем. -К.:Выща шк. Головное изд-во, 1988.

71. Романов В.А. Автоматическое управление в сложных объектах и системах. Конспект лекций. JI.:Jla6. полигр. маш. политехи, инст., 1970.

72. Фельдбаум А.А. Основы теории оптимальных автоматических систем. Физматгиз, 1963.

73. Бесекерский В.А., Попов Е.П., Теория систем автоматического регулирования. -М.:Наука, 1975.

74. Иващенко Н.Н. Автоматическое регулирование. Теория и элементы систем. М.'.Машиностроение, 1973.

75. Изерман Р. Цифровые системы управления: Пер. с англ. М.: Мир, 1984.