УДК 681.513.5
ПОСТАНОВКА ЗАДАЧИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ ТЕПЛОВЫМ РЕЖИМОМ ЗДАНИЯ ПРИ КОМБИНИРОВАННОЙ СИСТЕМЕ ОТОПЛЕНИЯ
М.М. Тверской, Д.В. Румянцев
STATEMENT OF A PROBLEM OF BUILDING THERMAL CONDITION OPTIMAL CONTROL AT DUAL HEATING SYSTEM
M.M. Tverskoy, D.V. Rumyantsev
Рассмотрен принцип управления тепловым режимом здания при наличии в системе отопления радиаторного и воздушного приборов отопления. Сформулирована задача оптимального управления тепловым режимом здания в соответствии с предложенным критерием при наличии ограничений.
Ключевые слова: тепловой режим здания, комбинированное отопление, оптимальное управление.
The principle of thermal condition control of a building with air and radiator sources in the heating system is considered in the article. The problem of building thermal condition optimal control, in terms of given criteria and appropriate boundary conditions, is stated in the article.
Keywords: building thermal condition, dual heating system, optimal control.
Введение
Тепловой режим здания зависит от множества параметров, определяющих тепловую обстановку здания [1]. Для оценки этих параметров используются различные показатели. Показателем, характеризующим взаимодействие между людьми и окружающей тепловой обстановкой (тепловым режимом здания), является уровень теплового комфорта, определяемый интенсивностью теплообмена между человеком и окружающей его средой. Интенсивность отдачи тепла человеком зависит от тепловой обстановки в помещении, которая характеризуется температурой воздуха, радиационными условиями, а также влажностью и подвижностью внутреннего воздуха [1].
Американским обществом инженеров по отоплению, охлаждению и кондиционированию воздуха «ASHRAE» был предложен показатель теплового комфорта PMV (Predictive Mean Vote - прогнозируемая средняя оценка). Согласно [2] про-
гнозируемая средняя оценка (PMV) - показатель, с помощью которого прогнозируется среднее значение ощущения к температуре большой группы людей на основе баланса температуры тела человека. Баланс температуры достигается, когда вырабатываемое телом человека тепло равно потере телом тепла в окружающую среду. PMV введен международным стандартом 1807730 [3].
Однако выражение для расчета PMV является достаточно сложным и содержит параметры, на которые сложно и не предполагается воздействовать при управлении тепловым режимом здания. К таким параметрам, например, относятся скорость движения воздуха в помещении и его влажность.
Согласно ГОСТ 30494-96, регламентирующему параметры микроклимата в жилых и общественных зданиях, комплексным показателем, характеризующим тепловую обстановку в помещении, может являться результирующая температура, определяемая по следующей зависимости:
Тверской Михаил Михайлович - д-р техн. наук, профессор, заслуженный деятель науки и техники РФ, заведующий кафедрой автоматизации механосборочного производства, Южно-Уральский государственный университет; [email protected]
Румянцев Дмитрий Владимирович - аспирант кафедры автоматизации механосборочного производства, Южно-Уральский государственный университет; dimfirst@gmail. сот
Tverskoy Mikhail Mikhailovich - Doctor of Science (Engineering), Professor, Honored Worker of Science and Technology of the Russian Federation, Head of Mechanical Assembly Production Automation Department, South Ural State University; [email protected]
Rumyantsev Dmitry Vladimirovich - post-graduate student of Mechanical Assembly Production Automation Department, South Ural State University; [email protected]
Т3П = Ъ1Тг + Ъ2Та , (1)
где Та - температура воздуха в помещении, °С; Тг - радиационная температура помещения, °С; Ъ1, Ъ2 - коэффициенты, принимающие значения по
0,5 при скорости движения воздуха в помещении ниже 0,2 м/с или 0,4 и 0,6 соответственно при скорости движения воздуха от 0,2 до 0,6 м/с.
Радиационная температура помещения может быть рассчитана как:
I 3Т
т =■
(2)
где <%■ - площадь внутренней поверхности ограждений и наружной поверхности отопительных приборов, м2; Тг - температура внутренней поверхности ограждений и наружной поверхности отопительных приборов, °С.
Результирующая температура является более предпочтительным показателем, характеризующим тепловой комфорт в помещении, чем индекс РМУ, поскольку при использовании комбинированного отопления (с радиационными и конвективными источниками тепла) такие параметры Тзи, как температура воздуха в помещении и радиационная температура, являются контролируемыми.
Комбинированные системы отопления
В настоящее время довольно часто используются системы отопления, содержащие комбинацию нескольких различных приборов отопления, например, воздушного и радиаторного с жидкостным теплоносителем. Также встречаются приборы отопления, состоящие из радиатора, создающего радиационную составляющую, и воздушного вентилятора, при включении которого обеспечивается повышение конвективной составляющей прибора [4]. Кроме того, в некоторых случаях в качестве воздушного отопления можно использовать систему приточной вентиляции с нагревом воздуха от теплоносителя, при наличии возможности переключения ее в режим рециркуляции воздуха.
Причины, по которым используются комбинированные системы отопления, могут быть различными [5, 6]. Воздушное отопление обладает меньшей инерционностью и позволяет нагреть воздух в помещении за относительно короткий промежуток времени. Температура внутренней поверхности стен помещения может быть ниже, чем при радиаторном отоплении, что снижает трансмиссионные потери тепла (через ограждающие конструкции). Радиаторное отопление содержит радиационную составляющую, позволяющую предавать тепло за счет инфракрасного излучения. Поэтому при радиаторном отоплении температура воздуха помещения может иметь более низкое значение, чем при воздушном отоплении. Это снижает инфильтрационные тепловые потери, которые имеют значительную величину при высоком воздухообмене, например, при включенной принудительной вентиляции. В соответствии с [7]
конвективные приборы воздушного отопления могут быть предпочтительнее в помещениях с отсутствующей или отключенной принудительной вентиляцией, а также при низком уровне воздухообмена. Но поскольку воздушное отопление требует дополнительных затрат электроэнергии, вызванных наличием воздушного вентилятора, то эффективность его применения с учетом суммарной энергии будет зависеть от стоимости электроэнергии. В то же время согласно [7] применение отопления с радиационной составляющей, каким является радиаторное отопление, более выгодно с точки зрения потребления энергии при высоком воздухообмене.
Комбинация приборов отопления с различными свойствами позволяет использовать их индивидуальные преимущества в единой системе отопления и избавиться от основных недостатков.
В России комбинированное отопление наиболее часто применяется в выставочных залах, торговых центрах и других объектах с прерывистым графиком работы. Комбинированные системы отопления работают по следующему принципу. Радиаторное отопление рассчитывается как дежурное из условия обеспечения минимально допустимой температуры воздуха в помещении на уровне 12 °С [8]. В качестве основного отопления используется воздушное, что не всегда является выгодным с точки зрения потребления энергии.
Функционирование системы автоматического
управления тепловым режимом здания
при комбинированной системе отопления
Рассмотрим принцип, по которому должна работать система управления комбинированным отоплением с прерывистым графиком работы. При настройке такой системы управления отоплением пользователь, кроме уставок результирующей температуры в дежурном и рабочем режимах соответственно, задает время начала и окончания рабочего режима для помещения, а также стоимостные показатели тепловой и электрической энергии с учетом ночного и дневного тарифов.
Задача системы управления заключается в автоматическом поддержании заданной результирующей температуры в соответствующих режимах при минимальном потреблении суммарной приведенной энергии. Поскольку здание и система отопления имеют значительную инерционность, то при переходе из дежурного режима в рабочий тепловая нагрузка системы отопления должна быть увеличена заблаговременно таким образом, чтобы в заданный момент времени начала рабочего режима уровень результирующей температуры в помещении был не менее минимально допустимого уровня, необходимого для комфортного состояния людей. Также с целью экономии энергии тепловая нагрузка системы отопления должна быть снижена за такой промежуток времени до окончания рабочего режима, чтобы в заданный момент
окончания рабочего режима значение результирующий температуры в помещении соответствовало минимально допустимому уровню. Моменты времени начала увеличения и снижения тепловой нагрузки система управления должна определять исходя из заданного графика работы, статических и динамических параметров воздушной и радиаторной систем отопления и самого помещения, а также распределять нагрузку между источниками тепла таким образом, чтобы количество потребляемой приведенной энергии за цикл работы (24 часа) было минимальным.
Подобные задачи рассмотрены в работах [9, 10], но представленный в них алгоритм ориентирован на систему с однотипными приборами отопления и требует переработки для использования в комбинированной системе.
Постановка задачи
Целью оптимального управления тепловым режимом здания при комбинированной системе отопления является минимизация общих затрат энергии на отопление при поддержании заданного уровня теплового комфорта.
В связи с наличием различных видов потребляемой энергии в системе отопления (тепловой и электрической), а также различных тарифов на указанную энергию в зависимости от времени суток будем выражать суточные затраты на энергию в виде следующего стоимостного показателя:
п т
Се = ЕСЭ (( + Е?В) + £ СТ ( + Е]В), (3)
¿=1 у=1
Э ]
где с; , су - стоимостные коэффициенты для соответствующего тарифа электрической и тепловой
ЭР ЭВ
энергии соответственно; Ej , Ej - количество
электрической энергии, потребленное по ¿-му тарифу радиаторным и воздушным отоплением со-
Г-.ТР г-ТВ
ответственно; Еу , Еу - количество тепловой
энергии, потребленное по у-му тарифу радиаторным и воздушным отоплением соответственно; п, т - количество различных тарифов, действующих в течение цикла работы системы отопления по заданному графику (за 24 часа) для электрической и тепловой энергии соответственно.
Количество потребленной электрической и тепловой энергии за время определенного тарифа можно рассчитать по следующим формулам:
гг гг
еЭр = | д(г)ЭР аг; ЕгЭВ = } д(г)ЭВ аг; (4)
у ^
Е]Р = } д(г)ТР аг; Е]в = | д(г)ТВ аг, (5)
г0 у г0 у
ЭР ЭВ
где д , д - электрическая мощность, потребляемая радиаторным и воздушным отоплением
ТР ТВ
соответственно; д , д - тепловая мощность,
потребляемая радиаторным и воздушным отоплением соответственно; ¿0; , ^ - пределы интегрирования, представляющие собой время начала и окончания ;-го тарифа для электрической энергии соответственно; г0у , гу - пределы интегрирования, представляющие собой время начала и окончания у-го тарифа для тепловой энергии соответственно.
Значения тепловых и электрических мощностей, используемых в (4) и (5), могут быть получены из математической модели объекта управления.
Объектом управления в данном случае будет являться отапливаемое помещение. В результате исследования, проведенного кафедрой строительной физики Чалмерского технологического университета (Швеция), была получена математическая модель такого объекта [11]. Она включает в себя выражения, описывающие статические и динамические характеристики отдельных элементов и конструкций здания, но при этом имеет достаточно сложную структуру и требует упрощения. Модель, разработанная в диссертационной работе [12], имеет более простую структуру и в то же время позволяет получить динамику изменения радиационной температуры и температуры внутреннего воздуха помещения. Однако в [12] не учитывается инерционность приборов отопления, что может потребовать доработки данной модели.
Из [11] путем упрощения и преобразования была получена математическая модель, представляющая собой систему уравнений в переменных состояниях:
аТа
аг +я1,б б<
- = а1,1Та а1,2Тзиг/ + а1,50 +
+ а Т '
сот ^“1,11 а,уеп >
ат
Г = а2,1ТаДП а2,2Тшг/ + а2,3Т» + а2,70гаа ;
аг
ат
аг
™ = а3,2Тзиг/ а3,3Т» + а3,4Тш/,ехг;
ат.
зиг/ ,ехг
аг аОВ
= а4,3Т» а4,4Тзиг/,ех + а4,10Та,ехг ;
(б)
аг
aQrad
= -а5 5£>В + а5 8иВ;
аг ад
сопу
= аб 7б™а + аб9и ;
р р
а7 бдсопу + а7 9и ,
аг
где ащп - постоянный коэффициент, соответствующий номеру строки т и столбца п матрицы коэффициентов правой части системы уравнений (б); переменные состояния: Та ¡п( - температура внутреннего воздуха; Тзи/ - температура внутренней поверхности ограждения; Т.№ - температура слоя внутри ограждения; Тзи/,ехг - температура наруж-
ной поверхности ограждения; 0 - конвективные
теплопоступления от воздушного отопления; 0Р -радиационные теплопоступления от радиаторного отопления; 0Сропу - конвективные теплопоступления от радиаторного отопления; управляющие воздействия: иР - сигнал управления приборами радиаторного отопления; иВ - сигнал управления приборами воздушного отопления; возмущающие
воздействия: Та
температура наружного воз-
духа; Тауеп - температура воздуха, поступающего в помещение от приточной вентиляции.
Согласно (1) и (2) уравнение выхода будет иметь вид:
(7)
Пусть температуры внутренних поверхностей ограждений одинаковы и равны Тзи/, температуры
наружной поверхности приборов радиаторного ото-
т-Р
пления также примем одинаковыми и равными Т , суммарную площадь внутренних поверхностей ограждений и наружных поверхностей приборов радиаторного отопления обозначим как Бзиг/ и £Р соответственно. Тогда (7) примет вид:
%шггТшН + 5 РТ Р
Ти = ъ -----+ Ь2Та .
+ 5
(8)
Тепловую мощность радиаторного отопления можно рассчитать по следующей формуле:
дТР = ор„а + С,, = кг*г ( - Т«п, )■ (9)
где КР - коэффициент теплопередачи приборов радиаторного отопления.
Выразив из (9) Т Р и подставив в (8), после соответствующих преобразований получим уравнение выхода в следующем виде:
Тг = 21Та,1п1 + 22Тшг/ + 250Сопу + 2бОТаа , (10)
где 2п - постоянный коэффициент, с номером п, определяемый в процессе идентификации математической модели объекта управления.
Коэффициенты из системы (б) и уравнения (10) могут быть получены в процессе параметрической идентификации модели объекта и зависят от таких постоянных, как тепловая емкость воздуха и материала ограждающих конструкций; тепловое сопротивление материала ограждающих конструкций; площадь внутренних поверхностей ограждения и наружных поверхностей приборов радиаторного отопления; плотность воздуха; коэффициенты передачи приборов воздушного и радиаторного отопления; постоянные времени приборов воздушного и радиаторного отопления.
Ограничения, накладываемые переменные состояния, включают в себя ограничения на температуру воздуха и радиационную температуру в
помещении. Пределы изменения этих переменных ограничены по ГОСТ 30494-96 [13]. Поскольку помещение эксплуатируется в двух режимах (рабочем и дежурном), то целесообразно сформировать различные ограничения вектора переменных состояния для рабочего и дежурного режимов. Важно также учитывать ограничение тепловой мощности приборов отопления.
Введем ошибку поддержания результирующей температуры:
е(г) = Т£и (г)- Тда (г), (11)
где Т$и (г) и Т$и (г) - установленная и реальная
траектория изменения результирующей температуры во времени соответственно.
Представив ошибку результирующей температуры как ограничение, получим однокритериальную задачу оптимального управления с критерием оптимальности:
(12)
рдеж . '■'шах >
при ограничениях:
ераб < ераб (г) < ераб • едеж < едеж (г) < е1
Ш1п _ \ ) ~ шах > е V / _ з
Траб < траб (г)< Траб • Тдеж < т^еж (г)< Тдеж . Т а,ш1п < Т а V / < Т а,шах ; Т а,ш1п < Т а \г / < Т а,шах ;
Траб < траб (г) < траб • Тдеж < т^еж () < Тдеж . Тг,ш1п < Тг \г/ < Тг,шах ; Тг,ш1п < Тг V/< Тг,шах ;
0В(г)<оВах; 0р (г)<0шах,
где Се - суточные затраты на энергию, определяе мые по зависимости (3); е
раб ераб ераб едеж
ш1п ’ шах > -
едеж , ешаж - ошибка, ее нижняя и верхняя грани-
цы соответственно, в рабочем и дежурном режи-
а ; траб траб траб тдеж тдеж тдеж
мах; Т а , Т а,ш1п, а,шах, Т а , а,ш1п, а,шах -
температура воздуха, ее нижняя и верхняя границы соответственно, в рабочем и дежурном режи-
мах;
Т раб т раб т раб т Деж Т Деж Т Деж
, Т , -írшяx, , г,ш1п , г,шах —
ра-
диационная температура, ее нижняя и верхняя границы соответственно, в рабочем и дежурном
режимах; еВ , бЩ^ах , 0е, QШa
- тепловая мощ-
ность, ее верхняя граница соответственно, для приборов воздушного и радиаторного отопления.
Переключение между рабочим и дежурным режимами происходит за счет изменения уставок и допустимых границ ошибки, температуры воздуха и радиационной температуры в помещении.
Заключение
Несмотря на то что системы комбинированного отопления находят свое применение в зданиях различного назначения, на сегодняшний день существующие алгоритмы управления такими системами не используют весь потенциал для экономии затрат на энергию.
Принцип работы системы автоматического управления тепловым режимом здания при комби-
нированной системе отопления, рассмотренный в статье, позволяет более рационально использовать энергию при работе по заданному графику. При этом учитываются динамические и статические свойства как источников тепла, так и объекта управления, а также тарифы на электрическую и тепловую энергию, потребляемую системой отопления.
Решение поставленной задачи оптимального управления позволит наилучшим образом распределять тепловую нагрузку между приборами отопления, что снизит суммарные затраты энергии на систему отопления, сохраняя значение результирующей температуры, характеризующей уровень теплового комфорта в помещении, в заданных пределах.
Литература
1. Отопление и вентиляция: учеб. для вузов. В 2 ч. Ч. 1: Отопление / П.Н. Каменев, A.H. Скана-ви, В.Н. Богословский и др. - 3-е изд., перераб. и доп. - М.: Стройиздат, 1975. - 483 с.
2. ГОСТ Р ИСО 7730-2009. Эргономика термальной среды. Аналитическое определение и ин-терпертация комфортности теплового режима с использованием расчета показателей PMV и PPD и критериев локального теплового комфорта. -М.: Стандартинформ, 2011. - 48 с.
3. ISO 7730:2005(Е). Ergonomics of the thermal environment - Analytical determination and interpretation of thermal comfort using calculation of the PMV and PPD indices and local thermal comfort criteria. - Switzerland, 2005. - 60 с.
4. Табунщиков, Ю.А. Экспериментальное исследование оптимального управления расходом энергии /Ю.А. Табунщиков, М.М. Бродач //АВОК. -
2006. - № 1. - С. 32-36.
5. Бродач, М.М. Отопление соборов - практика альтернативных решений / М.М. Бродач // «АВОК». - 2004. - № 2. - http://www.abok.ru/ for_spec/articles.php?nid=2342
6. Грудинин, В. С. Адаптивная компьютерная система управления микроклиматом / В. С. Грудинин // Аграрная наука Евро-Северо-Востока. -
2007. - № 9. - С. 137-142.
7. Малявина, Е.Г. Тепловая нагрузка на системы лучистого отопления. Сравнительный анализ / Е.Г. Малявина // АВОК. - 2009. - № 7. -С. 48-58.
8. СНиП 41-01-2003. Отопление, вентиляция и кондиционирование. - М.: Госстрой России: ГУ ЦПП, 2008. - 58 с.
9. Панферов, В. И. Об экономии энергии при оптимальном управлении режимом прерывистого отопления / В.И. Панферов, Е.Ю. Анисимова // Инженерные системы. - 2009. - № 5. - С. 32-35.
10. Панферов, В.И. Об оптимальном управле-
нии отоплением зданий как процессом с распределенными параметрами / В.И. Панферов, Е.Ю. Анисимова //Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». -
2009. - Вып. 9. - № 3 (136). - С. 24-28.
11. Angela, S.K. HAM-Tools International Building Physics Toolbox Block documentation Report: R-02:6 / S.K. Angela. - Sweden: Department of Building Physics Chalmers Institute of Technology,
2008. - 58 с.
12. Нагорная, А. Н. Математическое моделирование и исследование нестационарного теплового режима зданий: дис. ... канд. техн. наук /А. Н. Нагорная. - Челябинск, 2008. - 150 с.
13. ГОСТ 30494-96. Здания жилые и общественные. Параметры микроклимата в помещениях. -М.: Госстрой России: ГУ ЦПП, 1999. - 7 с.
Поступила в редакцию 9 февраля 2012 г