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

УДК 681.5

СИНТЕЗ СИСТЕМЫ УПРЕЖДАЮЩЕГО УПРАВЛЕНИЯ ПРОЦЕССОМ ПОДАЧИ ТЕПЛА НА ОТОПЛЕНИЕ ЗДАНИЯ

С.А. Добротин, Е.Л. Прокопчук Дзержинский политехнический институт,

(филиал Нижегородского государственного технического университета

им. Р.Е. Алексеева)

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

Ключевые слова: упреждающее управление, алгоритм, прогнозирование, отопление, синтез, здание, тепло.

SINTEZA SISTEMULUI DE DIRIJARE A PROCESULUI DE LIVRARE A CALDURII LA

fNCALZIRE A BLOCULUI S.A. Dobrotin, E.L.Procopciuc Institutl Politehnic din Dzerjinsc,

(filiala a Univesitatii Tehnice din Nijnii Novgorod , R.E. Alexeev)

Rezumat. Este sintezat un sistem combinat de dirijare cu prognozare a incalzirii blocului. Sunt prezentati algoritm de prognozare de durata foarete scurta al temperaturii §i algoritm de dirijare cu prognozare a livrarii caldurii pentru incalzirea blocului tinand cont de efectele diferite ale influentei climei din exterior la fatadele blocului.

Cuvinte-cheie: dirijare cu prognozare, algoritm, pronosticare, incalzire, sinteza, bloc, caldura.

PREDICTIVE CONTROL SYSTEM SYNTHESIS OF DISTRICT HEATING

SYSTEM S.A. Dobrotin, E.L. Prokopchuk Dzerjinsk Politechnical Institute,

(branch of the Nijniy Novgorod Technical University, R.E. Alexeev)

Abstract. A combined system of building heating predictive control has been synthesized in the research. Following algorithms are described: adaptive algorithm of very short-term temperature forecast and algorithm of anticipatory control of heat supply to building heating taking into account different outside climate influence on elevations of a building.

Key words: predictive control, algorithm, forecasting, heating, synthesis, building, heat.

Структурный синтез системы упреждающего управления (ССУ) процессом подачи тепла на отопление здания

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

53

В связи с вышесказанным предложена прогнозно-компенсационная схема управления подачей тепла на отопление здания [3-5] (рис. 1). Управление по данной схеме строится на основе прогнозирования параметров наружного климата, то есть на вход динамического компенсатора в разомкнутом контуре системы управления подается прогнозная информация о возмущении, формируемая алгоритмами, заложенными в блоке прогнозирования.

1 - блок прогнозирования; 2 - динамический компенсатор; 3 - регулятор; 4,5 -сумматоры; / - возмущающее воздействие; и - управляющее воздействие /рг -прогнозируемое возмущение; їрг, - заданная температура в контрольном помещении; є - отклонение текущей температуры в контрольном помещении от заданной (рассогласование).

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

отопление здания

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

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

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

Динамический компенсатор (поз. 2, рисунок 1) введен в структурную схему для компенсации инерционности объекта по каналу возмущения. Передаточная функция динамического компенсатора (Wk(s)) находится из условия инвариантности системы отопления по отношению к возмущающему воздействию [6] и имеет вид:

uni \ /kl' е s(r1 г2) . (Top ■ s +1) ч -s-(-z2 К

Wk(s) = -(-------- v *2---------- +k2■ (Top■ s +1)■ e ' )), (1)

(Tst ■ s +1)

где kl, k2 - обобщенные коэффициенты усиления теплоемких и нетеплоемких ограждающих конструкций соответственно, т1 - транспортное запаздывание по каналу возмущения, т2- время прохождения теплоносителя по трубопроводу, Тор - постоянная

54

времени радиатора, - постоянная времени теплоемких ограждающих

конструкций (стенка).

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

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

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

Жк (я) =

7 1 _ 1

к1 * е

+ к 2

(Тяі * я +1)2

(2)

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

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

1 - СУ, регулирующая подачу теплоносителя по текущему значению наружной температуры; 2 - СУ, регулирующая подачу теплоносителя с учетом динамических характеристик объекта по каналу возмущения; тгеё1 - время регулирования СУ1; тгеё2 - время регулирования СУ2;

&лоР - допустимый диапазон температуры; ЖорЛ - оптимальный диапазон температуры.

Рис. 2. Графики изменения температуры в отапливаемом помещении

1 - = 8.3 часа; 2 - = 20 часов.

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

На рис. 4 для сравнения приведены графики изменения температуры в отапливаемом помещении при работе разных СУ. Из рисунка видно, что снижение негативного влияния динамических характеристик по каналу управления можно достичь регулированием с упреждением по прогнозируемым метеорологическим параметрам.

т, ч

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

Рис. 4. Графики изменения температуры в отапливаемом помещении

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

Таблица 1 - Сравнение систем управления

СУ ІМт °С Тгв8, ч 1т, °Оч

1 3.30 34 70.45

2 2.47 4 10.12

3 0.92 0 3.0

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

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

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

Имитационное моделирование производилось с помощью пакета визуального программирования Simulink [7].

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

57

На рисунке 5 приведены переходные процессы изменения температуры в помещении для следующих параметров канала управления: т2 = 0,7 ч, Тор = 0,8 ч, Т\о = 1,35 ч (Т\о

- постоянная времени воздушного объема).

т, ч

1 - трг = 0 мин; 2 - трг = 30 мин; 3 - трг = 45 мин; 4 - трг = 60 мин; 5 - трг = 75 мин;

6 - трг = 90 мин; 7 - трг = 105 мин.

Рис. 5. Графики изменения температуры в отапливаемом помещении, при разных

интервалах прогнозирования

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

(3)

19 < 1 < 21

Для решения поставленной задачи оптимизации использовался метод сканирования, результаты приведены на рисунке 6 (на рисунке 6, N - номер эксперимента).

Т

0

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

Рис. 7. Зависимость интервала прогнозирования от инерционности помещения

Экспериментально полученные данные были аппроксимированы по методу наименьших квадратов полиномом первого порядка:

где к = 1.

1рГ

Ъг + к'тор ,

(4)

Адаптивный алгоритм сверхкраткосрочного прогнозирования температуры наружного воздуха

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

59

/" г = к - р • ^(х)

^(г) = -0,0012 т3 + 0,062 т2 - 0,94 т - 3,5 Г

<

Р =

-4,89; -2,14 -2,40; -1,20 -2,40; -0,68 -2,40; -0,76 -3,76; -1,16 -6,49; -3,55 ,-9,77; -6,23_у

(5)

Полученная на ее основе формула для прогнозирования имеет вид:

грГ = к) - р- g(x + Хрг)

(6)

Алгоритм осуществляется в несколько этапов:

1. Задание допустимой ошибки прогнозирования и интервала прогнозирования. Допустимая ошибка прогнозирования и интервал прогнозирования задаются либо

со станции оператора, либо непосредственно вводятся в контроллер.

2. Определение текущего месяца.

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

3. Определение максимального значения амплитуды суточного хода температуры.

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

Месяц Р

солнечный пасмурный

Октябрь -4,89 -2,14

Ноябрь -2,4 -1,2

Декабря -2,4 -0,68

Январь -2,4 -0,76

Февраль -3,76 -1,16

Март -6,49 -3,55

Апрель -9,77 -6,23

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

60

Расчетные значения текущей температуры определяются по разработанной модели

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

Выбор р осуществляется в зависимости от полученной на текущей момент ошибки прогнозирования, которая определяется как разность между текущим действительным значением температуры и текущим расчетным значением температуры (погрешности будет две: для солнечного и облачного дня):

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

Если условие не выполняется, из полученных ранее Р^ и ро, выбирается Р, соответствующая минимальной А.

Если условие выполняется, происходит адаптация модели (3.10), которая заключается в определении скорректированного ро. Скорректированное Ра определяется по модели (3.10) из условия совпадения текущей расчетной и действительной температур - г = Ьвк-

6. Определение прогнозного значения температуры.

Прогнозное значение температуры - 1ргг определяется по формуле (6).

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

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

Алгоритм упреждающего управления процессом подачи тепла на отопление здания

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

В соответствии с рисунком 1 алгоритм [8] будет состоять из 3 основных блоков:

- блок 1 - блок прогнозирования, на основании разработанного адаптивного алгоритма сверхкраткосрочного прогнозирования температуры (рисунок 8) формирует прогнозное значение возмущающего воздействия;

- блок 2 - динамический компенсатор, по прогнозному значению возмущения формирует управляющее воздействие;

- блок 3 - блок корректировки, корректирует управляющее воздействие для каждого фасада.

Ьвк'* = к - Р/ £-(т)

іівкг’° = к + ро' £(т)

(7)

(8)

(9)

(10)

А^ > Е и Д0 >Е

(11)

Ра = ^0 - Ы) / £(Т)

(12)

61

наружного воздуха

62

PROBLEMELE ENERGETICII REGIONALE 2(1в) 2011

_____________________________________________в3

© Добротин С. А., Прокопчук Е. Л., 2G11 г.

Рис. 9. Алгоритм упреждающего управления процессом подачи тепла на отопление

здания

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

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

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

В случае получения дополнительного тепла, например, от солнечного излучения одним из фасадов и, как следствие, увеличения температуры в помещениях данного фасада (гр1 > гр1г, гр2 > гР2г) управляющее воздействие на этот фасад корректируется в блоке 3, который реализует ПИ закон регулирования.

Алгоритм упреждающего управления процессом подачи тепла на отопление здания приведен на рис. 9.

Выводы

1. Получено математическое описание динамического компенсатора,

использованного как звено разомкнутого контура СУУ.

2. Установлена зависимость интервала прогнозирования от параметров

инерционности помещения.

3. Выполнен структурный и параметрический синтез комбинированной

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

Литература

1. Прокопчук, Е.Л. Комбинированная система управления отоплением здания / Е.Л. Прокопчук, С.А. Добротин // Тезисы докладов IV Международной молодежной научно-технической конференции. Н. Новгород, 2005. - С. 258.

2. Прокопчук, Е.Л. Диагностический подход к управлению отоплением здания / Е.Л. Прокопчук // VI Международная конференция «Неразрушающий контроль и техническая диагностика в промышленности» Москва, 2007. - С. 100-101.

3. Прокопчук, Е.Л. Регулирование отпуска тепловой энергии по прогнознокомпенсационной схеме / Е.Л. Прокопчук, С.А. Добротин // Тезисы доклада VI Международной молодежной научно-технической конференции «Будущее технической науки» Н. Новгород, 2007. - С. 268.

4. Прокопчук, Е.Л. Прогнозно-компенсационная схема управления отоплением здания / Е.Л. Прокопчук // Материалы XII Нижегородской сессии молодых ученых (технические науки) Н.Новгород, 2007. - С. 125.

5. Пат. 73509 Российская Федерация, МПК51 О 05 Б 23/19 Система управления процессом подачи тепла на отопление здания / Прокопчук Е.Л., Добротин С.А.;

64

заявитель и потентообладатель Государственное образовательное учреждение высшего профессионального образования Нижегородский государственный технический университет им. Р.Е.Алексеева. - № 2008100508/22; заявл. 09.01.2008; опубл. 20.05.2008.

6. Автоматическое управление в химической промышленности: Учебник для вузов. Под ред. Е. Г. Дудникова. - М.; Химия, 1987. - 368 с.

7. Прокопчук, Е.Л. Применение имитационного моделирования на этапе синтеза системы автоматического регулирования подачи тепла на отопление здания / Е.Л. Прокопчук, С.А. Добротин // Успехи современного естествознания. 2008, №3. - С. 50-51.

8. Прокопчук, Е.Л. Алгоритм упреждающего управления процессом подачи тепла на отопление здания / Е.Л. Прокопчук // Успехи современного естествознания. 2008, №3. -С. 60-62.

Сведения об авторах:

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

Прокопчук Елена Леонидовна, магистр техники и технологии, старший преподаватель кафедры "Автоматизация и информационные системы". Область научных интересов: автоматизация

технологических процессов, моделирование систем управления. тел. (8313) 26-59-16, факс (8313) 23-04-10, e-mail: llp@hotmail.ru