ВведениеИсследование систем искусственного интеллекта в настоящее время является актуальным направлением в области ComputerScience. К искусственному интеллекту относятся такие техно-логии, как машинное обучение, компьютерное зрение, методы оптимизации, в том числе гене-тический алгоритм. Такие технологии широко применяются для решения различных экономи-ческих задач, например в задачах моделирования высокотехнологических производственных предприятий [1]. Методы оптимизации применяются в решении разнообразных задач, таких как машинное обучение, составление расписаний, планирование. Наиболее актуально применение методов оптимизации в экономических задачах, где требуется минимизировать или, напротив, максимизировать определенные экономические показатели, например издержки, размер прибы-ли, время простоя, производительность. Аналогичные задачи существуют в области энергосбе-режения [2], в частности, актуальным является решение оптимизационных задач в системах централизованного отопления, т. к. это положительно влияет на энергосбережение и экологиче-скую обстановку. К таким задачам относится устранение температурного дисбаланса. Этот дис-баланс выражается в том, что есть «теплые» здания и «холодные» здания. То же самое относится и к помещениям внутри здания. Таким образом, проблема существует на двух уровнях. Она решается путем увеличения расходов теплоносителя или повышения его температуры в зданиях, где есть «холодные» комнаты. В результате в «теплых» помещениях температура становится еще выше. Такой температурный режим очень некомфортен для жильцов. В этом случае они просто открывают окна. По сути, жильцы «теплых» квартир «отапливают» улицу, что приводит к нерациональному расходу ресурсов, в частности угля, которым топят котельные, что, кроме всего прочего, увеличивает выбросы в атмосферу парниковых газов. На этапе проектирования зданий проблему «теплых» и «холодных» помещений можно решить путем подбора отопительных приборов оптимальной мощности, которая обеспечивает нормативную температуру в помещении. В действующих зданиях можно установить дополни-тельные радиаторы в «холодных комнатах», однако это не всегда приводит к устранению тем-пературного дисбаланса. Второй способ заключается в регулировании гидравлических сопро-тивлений клапанов стояков или отопительных приборов. Однако применение данного способа влечет за собой проблему системной сложности и нелинейности. Дело в том, что если мы путем регулирования клапана изменяем гидравлическое сопротивление в одной комнате, это влечет за собой изменение температуры не только в этой комнате, но также в других. Данная зависимость нелинейная. Следовательно, данная задача сводится к решению системы нелинейных уравнений размерностью порядка сотен (по количеству регулировочных клапанов во всех комнатах отап-ливаемого здания). В настоящее время аналитического способа решения данной математической задачи не существует, поэтому она решается приближенным итерационным методом. К таким методам относятся математические методы оптимизации. Анализ существующих решенийКомпьютерная программа, в которой можно выполнять расчет тепловых сетей разных уровней, представлена в [3]. При эксплуатации данного ПО используются методики и алгорит-мы, предложенные авторами. Данный программный продукт содержит функционал по расчету гидравлических сетей и равновесных температур помещений, однако возможность оптимизации температурных режимов в нем отсутствует. Решение задачи оперативного регулирования теплоснабжения было предложено в [4], где авторы применили имитационное моделирование зависимости температуры в помещениях от температуры наружного воздуха и других факторов. Предложенный ими метод основывался на штрафных функциях , которые зависят от температуры в помещении и количества потребляе-мой энергии. Данный метод позволяет экономить тепловую энергию. Следует отметить, что его применение не приводит к устранению температурного дисбаланса, т. к. нельзя просто запре-тить открывать форточку жильцам «теплых» квартир. Задачи многокритериальной оптимизации рассмотрены в [5] – они, как правило, возника-ют перед началом отопительного сезона. Авторы рассмотрели задачи оптимизации распредели-тельных тепловых сетей, имеющих древовидную и однолинейную структуру. Критерии оптими-зации: минимизация мест приложения усилий, потери теплоносителя и риски возникновения критических ситуаций.Другая идея представлена в работе [6]. Ее авторы предложили оптимизировать структуру тепловой сети. Метод состоит в применении алгоритма перебора деревьев и алгоритма имитации отжига. Однако предложенную методику можно применить только для вновь создаваемых тепловых сетей, т. к. экономически невыгодно перестроить существующую. В [7] авторы решали задачу оптимального отопления зданий. Их цель – оптимизировать энергопотребление в системе центрального отопления путем определения достижения минимума потерь энергии обогреваемого здания. Решалась двухкритериальная задача: с одной стороны, необходимо поддерживать комфортную температуру, но одновременно минимизировать потери тепловой энергии. Данные критерии являются противоречивыми. Авторы нашли простейшие частные решения и предложили расчетные формулы. В результате было предложено снизить коэффициент теплоотдачи через наружные стены (это весьма дорогостоящее мероприятие) и строго соблюдать температурный график, индивидуальный для каждого здания. Данные решения не всегда приемлемы. Кроме того, авторами не было предложено регулировать температуру внутри помещений.Отдельно можно отметить ряд зарубежных исследований, посвященных использованию методов искусственного интеллекта для решения задач энергосбережения. В работе [8] авторы представили иерархическую оптимизацию сети станций зарядки электромобилей. Была пред-ложена трехслойная модель, включающая (на первом слое) оптимизацию местоположений за-правочных станций, далее – модель очередей (второй слой) и, наконец, третий слой – оптимиза-ция динамики зарядки аккумулятора. В работе [9] была обнародована крикинг-модель, назначение которой – снижать потери энергии в распределительных сетях. Для ускорения процесса поиска решения авторами был предложен алгоритм аппроксимации сложных распределительных сетей. Применение искусственного интеллекта актуально во всех сферах энергосбережения, включая оптимизацию сетей центрального отопления и оптимизацию в сфере потребления электроэнергии. В частности, применение метода Монте-Карло для моделирования надежности оборудования сетевой инфраструктуры электрических сетей представлено в статье [10], авторы предлагают оптимизацию с помощью моделирования для поиска оптимального набора сетевых активов для повышения надежности системы. В частности, алгоритм Industrial Strength COMPASS реализуется для поиска оптимального набора новых линий передачи, который мак-симизирует надежность системы с учетом бюджетных ограничений. В [11] описывается текущий режим подключения для распределенной генерации в Вели-кобритании и представлено новое применение метода оптимального потока мощности для ав-томатического управления потоком энергии для управления тепловыми ограничениями в рас-пределительных сетях. Авторы оценили и обосновали пригодность данного подхода для управ-ления распределительной сетью в рамках онлайн-схемы и доказали, что такая схема обеспечи-вает устойчивость управления.Таким образом, несмотря на многообразие решений в области энергосбережения, проблема устранения температурного дисбаланса рассмотрена недостаточно и в настоящий момент является актуальной, в связи с чем исследование алгоритмов оптимизации температурных ре-жимов отапливаемых помещений является актуальной задачей. Постановка задачиЦелевой функцией данной задачи оптимизации является среднеквадратичное отклонение от целевой температуры [2] , (1)где z – число балансировочных клапанов стояков или регулировочных клапанов отопительных приборов; i – порядковый номер помещения; N – число помещений; ti – расчетная равновесная температура в помещении; tопт – целевая температура в помещении(в рамках данного исследо-вания это нормативная температура, принятая 20 oC); sj – сопротивление j-го клапана; j – номер клапана, принимает значение от 1 до z.Оптимизируемые параметры – это гидравлическое сопротивление регулировочных клапа-нов отопительных приборов и/или балансировочных клапанов стояков (можно оптимизировать и то, и другое). Размерность аргументов целевой функции – порядка сотен (количество клапанов во всем здании). Ограничение – диапазон регулирования клапанов: , (2)где и – минимальное и максимальное сопротивление соответственно.Сопротивление клапанов влияет на расход теплоносителя xпр в данном тепловом приборе, а следовательно, на поток тепловой энергии от радиатора , который пропорционален его нагревательной площади , и зависит от температуры воздуха в помещении tв. Теплоотдача отопительного прибора рассчитывается по формуле , (3)где – номинальная мощность отопительного прибора; n, p – эмпирические коэффициенты, подбираемые экспериментальным путем; cпр – коэффициент, учитывающий схему подключения отопительного прибора; – коэффициент учета дополнительного теплового потока, который получается из-за округления расчетной величины; – коэффициент учета тепловых потерь отопительного прибора у наружных ограждений; , – температура теплоносителя на входе и выходе из отопительного прибора [3]. Следует отметить, что тепловой поток рассчитывается при стандартных условиях, когда средний температурный напор равен 70 оС, поэтому в форму-лу введен коэффициент 70. Коэффициент 360 – перевод единиц измерения. Расходы теплоносителя в отопительных приборах xпр и трубах xтр образуют вектор расхода теплоносителя системы отопления здания x, который зависит от структуры теплосети, гидравлических сопротивлений и определяется путем решения системы нелинейных уравнений Кирхгофа (4)где A – матрица соединений узлов и ветвей; B матрица контуров; неизвестный вектор пе-репада давлений; S – диагональная матрица гидравлических сопротивлений (включая сопротив-ления регулировочных клапанов радиаторов ); X – диагональная матрица модулей расходов; – вектор гидравлических напоров; – вектор расходов теплоносителя в узлах. От теплового потока, в свою очередь, зависит температура в отапливаемом помещении. При расчете решается уравнение теплового баланса, учитываются потери тепла в помещении через ограждающие конструкции, а также решается гидравлическое уравнение тепловой сети здания.Температура в помещении определяется путем решения системы нелинейных урав-нений теплового баланса, отдельно для каждого i-го помещения [3]: (5)где – теплообмен со стенами; – теплообмен с полом; – теплообмен с потолком; – теплообмен с дверями; – инфильтрация через двери; – теплообмен с окнами; – инфильтрация через окна; – инфильтрация воздуха; , – приход тепла от отопительных приборов и труб. Более подробно данный расчет описан в [12].Описание генетического алгоритма, применяемого для целей данного исследования В [12] нами проведено исследование генетического алгоритма, который решает задачу оп-тимизации гидравлических сопротивлений регулировочных клапанов отопительных приборов. Оптимизируется значение целевой функции (1). Генетическим алгоритмом является алгоритм эвристического поиска, который используется для решения оптимизационных задач и моделирования посредством случайного комбинирования, подбора и вариации искомых параметров. При этом используются механизмы, аналогичные естественному отбору в природе. Задача формализована таким образом, что ее решение может быть закодировано в «генотипе» – специального вида векторе генов. Каждый ген может быть либо битом, либо числом, либо неким другим объектом. Такой набор параметров часто называется «особь», а множество особей в пределах одной итерации оптимизации – «популяцией». Хромосомой обычно является набор параметров, которые подлежат оптимизации в рамках решаемой алгоритмом задачи, это могут быть числа, строки и другие данные, в том числе объекты.В рамках решаемой задачи хромосомой является массив чисел из диапазона (2). Количество чисел в массиве порядка сотен.Генетический алгоритм, используемый в настоящем исследовании, можно описать следующими шагами:1. Создается начальная выборка, в которой гидравлические сопротивления клапанов задаются случайным образом из диапазона (2). 2. Чтение из базы данных информации, необходимой для расчета целевой функции (1).3. Рассчитывается целевая функция (1), включая расчет теплоотдачи от отопительных приборов (3) и решение уравнений (4), (5).4. Выборка сортируется по возрастанию целевой функции (т. к. целевую функцию необходимо минимизировать).5. Путем скрещивания, вероятность которого вычисляется по формуле (6)где n – число «особей» в «популяции» (наборов в выборке); k1 и k2 – первый и второй набор из проверяемой на возможность скрещивания пары соответственно, и мутации, вероятность кото-рой вычисляется по формуле (7)где m – номер набора, подлежащего удалению (начиная с нуля), создаются новые элементы и добавляются в выборку. Мутацией является изменение гидравлических сопротивлений согласно формулам (8)где – среднеквадратическое отклонение случайной величины, имеющей гауссовское распреде-ление (выбрано в размере 0,1); – выбираемый эмпирическим путем шаг мутации; x1 и x2 – случайные числа (могут быть в пределах [0;1], распределение этих случайных величин равномерно (во многие языки программирования встроены подобные генераторы случайных чисел)), и (9)где – максимально допустимое отклонение температуры от нормативной, при котором в любом случае происходит мутация (вероятность мутации равна 1).После мутации выполняется проверка условия (2) и корректировка, если произошел выход за этот диапазон (устанавливается граничное значение). 6. Рассчитываются целевые функции у новых элементов.7. Выборка сортируется по убыванию целевой функции (т. к. целевую функцию необходимо минимизировать).8. Селекция выборки (случайным образом удаляются элементы, чем хуже значение целевой функции, тем больше вероятность удаления (7)).9. Если достигнут критерий остановки, выходим из алгоритма, иначе переходим к шагу 5.Теперь рассмотрим более подробно некоторые моменты генетического алгоритма.Скрещивание происходит по следующей схеме:1. Из родительской популяции выбираются пары хромосом. 2. Для каждой такой пары случайно выбирается точка скрещивания: позиция гена (локус) в хромосоме, которая делит ее на две части.3. Создается новая хромосома. В нее берется часть из первой хромосомы (та, что до точки скрещивания) и из второй (та, что после точки скрещивания).Вероятность скрещивания определется по формуле (6). Данная величина принимает значения от 0 до n – 1. Вероятность удаления элемента из выборки при селекции определяется по формуле (7).Для «особи» (набора), имеющей лучшую целевую функцию, данный номер равен 0, иными словами, он никогда не будет удален, деление происходит в цикле, который может быть повторен несколько раз до тех пор, пока количество наборов не достигнет значения меньшего, чем максимально допустимое количество. Если число наборов становится меньше, чем минимально допустимое количество, селекция не производится на данной итерации.Мутация происходит путем изменения гидравлического сопротивления на случайную ве-личину, вычисленную по формуле (8).Вероятность мутации 0,5, т. е. примерно половина всех гидравлических сопротивлений в наборе, который получен скрещиванием, подвергаются мутации. Модификация «Направленная мутация». Наряду с классической схемой мутации, опи-санной выше, в программном продукте, который был использован для настоящего исследования, предусмотрена модификация алгоритма «Направленная мутация». Идея состоит в том, чтобы объединить генетический алгоритм с другими видами алгоритмов, создав гибридный алгоритм. Из результатов вычислительных экспериментов следует, что гибридные алгоритмы позволяют достичь большей глубины оптимизации [13]. В данной модификации вероятность мутации зависит от температуры в помещении, где установлен прибор с клапаном, гидравлическое сопротивление которого подвергается мутации: чем больше отклонение температуры от требу-емой, тем больше вероятность мутации. Данная вероятность вычисляется по формуле (9). Изме-няется также сам способ вычисления размера мутации: где σ – среднеквадратическое отклонение случайной величины, имеющей гауссовское распре-деление (выбрано в размере 0,1); х1, х2 – случайные числа (могут быть в пределах [0;1)); b – смещение, зависящее от температуры в помещении, для которого происходит мутация гид-равлического сопротивления клапана, вычисляется по формуле где ti – температура в помещении; tтр – требуемая температура в помещении.Описание алгоритма направленного поиска, применяемого для целей данного исследованияДанный алгоритм решает задачу оптимизации гидравлических сопротивлений регулиро-вочных клапанов отопительных приборов. Описать его можно следующими шагами:1. Начальная инициализация. Создаются два набора оптимизируемых параметров: про-шлый и текущий, в которых гидравлические сопротивления клапанов задаются случайным об-разом из диапазона (2). 2. Чтение из базы данных информации, необходимой для расчета целевой функции (1). Рассчитывается целевая функция для каждого из наборов, включая расчет теплоотдачи от ото-пительных приборов (3) и решение уравнений (4), (5).3. Вычисляются приращения для следующей итерации: (10)где k – номер итерации; j = 1, …, z – номер клапана; – требуемая температура в помещении; – расчетная температура на k-ой итерации в помещении, где установлен прибор с клапаном j.4. Создается новый набор оптимизируемых параметров за счет прибавления приращений к текущему набору параметров (приращение может быть и отрицательным). 5. Вычисляется целевая функция для нового набора.6. Набор, названный «предыдущим», заменяется на «текущий».7. Набор, названный «текущим», заменяется «новым».8. Если выполняется условие останова, то выход из алгоритма, иначе переход к п. 4.Целевая функция рассчитывается точно так же, как и для генетического алгоритма. При-ращение сопротивлений клапанов вычисляется по формуле (10).Результаты экспериментовВ целях сравнения минимумов достижения целевой функции в рамках данного исследо-вания были проведены следующие вычислительные эксперименты:– оптимизация гидравлических сопротивлений клапанов отопительных приборов;– оптимизация гидравлических сопротивлений балансировочных клапанов отопительных стояков.В рамках данных вычислительных экспериментов был использован генетический алго-ритм, в котором шаг мутации был взят 20 и 50 в единицах Па/(Кг/час)2, и алгоритм направ-ленного поиска.В первом эксперименте исследовалась зависимость глубины оптимизации (значение целевой функции, которую удалось получить в ходе оптимизации) от следующих характеристик алгоритма:– шаг мутации (величина , которая подставляется в формулу (8));– минимальный размер популяции (минимальное допустимое количество, после которого селекция элементов прекращается);– максимальный размер популяции (количество элементов, больше которого не может быть в популяции, данный предел ограничивается селекцией, которая происходит до тех пор, пока размер популяции выше этого значения);– начальный размер популяции (количество элементов, которое создается случайным об-разом на первой итерации алгоритма).Размерность целевой функции:– в случае оптимизации клапанов отопительных приборов (количество регулировочных клапанов в отапливаемом здании) – 254;– в случае оптимизации клапанов стояков (количество балансировочных клапанов) – 44.Результаты эксперимента приведены в табл. 1. Таблица 1Результаты работы генетического алгоритма при разных значениях характеристик алгоритмаЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функции*минимальныймаксимальныйначальныйСтояки20530102,51550202,4532072,985032072,311550202,26Батареи20530103,301550202,9932073,425032073,301550202,8815200202,81* Значения больше 5 – «неудовлетворительно» (значительный дисбаланс); 3–4 – «удовлетворительно» (имеется дисбаланс); 1,5–3 – «хорошо» (небольшой дисбаланс); 0,5–1 – «отлично» (дисбаланс незначительный); менее 0,5 – «идеально» (дисбаланс прак-тически отсутствует).В табл. 2 отдельно выведены значения характеристик алгоритма, для которых получились наилучшие результаты, в табл. 3 – наихудшие.Таблица 2Наилучшие результаты оптимизацииЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйСтояки501550202,26Батареи2,882002,81Таблица 3Наихудшие результаты оптимизацииЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйСтояки2032072,98Батареи3,42Таким образом, наилучшие результаты были получены при наибольших значениях харак-теристик «Шаг мутации», «Начальное количество элементов», «Максимальное количество эле-ментов» и «Минимальное количество элементов». Наглядно разница прослеживается на графике истории оптимизации (рис. 1) для варианта с худшими и лучшими параметрами (для варианта оптимизации клапанов отопительных приборов). Рис. 1. История оптимизации гидравлических сопротивлений клапанов отопительных приборов генетическим алгоритмомДля выявления наиболее важной характеристики (которая оказывает наибольшее влияние на сходимость алгоритма оптимизации) были проведены эксперименты, когда только одна ха-рактеристика менялась, остальные были фиксированные. В табл. 4 приведены результаты эксперимента, когда менялось минимальное количество элементов в популяции. Таблица 4Результаты оптимизации в зависимости от минимального количества элементов в популяцииЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйБатареи2032073,4245681015Таким образом, глубина оптимизации (столбец «Значение целевой функции») не менялась при изменении минимального размера популяции. Это можно объяснить тем, что размер популяции никогда не достигает минимума или достигает очень редко. В табл. 5 приведены результаты эксперимента, когда менялось максимальное количество элементов в популяции. Таблица 5Результаты оптимизации в зависимости от максимального количества элементов в популяцииЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйБатареи2032073,42213,412325303,36403,32503,39Таким образом, результат оптимизации почти всегда (за исключением последней строки) улучшается с увеличением максимального количества элементов в популяции. Такая законо-мерность связана с тем, что количество элементов в популяции после селекции сохраняется близким к верхней границе (максимальному значению); отодвигая верхнюю границу, мы увели-чиваем количество элементов, участвующих в создании новых элементов (увеличиваем генети-ческое разнообразие), что ведет к увеличению вероятности найти лучшие решения. Однако из-за ограничения вычислительных ресурсов невозможно бесконечно увеличивать максимальное количество элементов в популяции. Следует также отметить, что время работы алгоритма пропорционально квадрату максимального количества элементов.Исключение из закономерности связано с тем, что генетический алгоритм основан на факторе случайности. В табл. 6 приведены результаты эксперимента, когда менялось начальное количество эле-ментов в популяции.Таблица 6Результаты оптимизации в зависимости от начального количества элементов в популяцииЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйБатареи2032073,4283,43103,32153,30203,10303,1150Здесь с ростом начального количества элементов улучшается глубина оптимизации. Это объясняется тем, что чем больше генетическое разнообразие в начале, тем больше вероятность того, что лучшее решение попадет в дальнейшие итерации. В табл. 7 приведены результаты эксперимента, когда менялся шаг мутации.Таблица 7Результаты оптимизации в зависимости от шага оптимизации элементов в популяцииЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйБатареи2032073,42213,37232530403,34503,30Согласно табл. 7 с увеличением шага мутации (изменчивости) увеличивается глубина оп-тимизации. Это объясняется тем, что с увеличением изменчивости растет охват пространства параметров оптимизации и повышается вероятность найти лучшее решение. Улучшает качество оптимизации также модификация «Направленная мутация». В табл. 8 приведено сравнение данной модификации с классическим вариантом алгоритма.Таблица 8Сравнение модификации «Направленная мутация» с классическим вариантом генетического алгоритмаЧто оптимизируемШаг мутацииРазмер популяцииЗначение целевой функцииминимальныймаксимальныйначальныйбез модификациис модификациейБатареи20530103,32,271550202,992,1432073,422,263,32,351550202,882,08Таким образом, применение модификации «Направленная мутация» во всех случаях уменьшает значение достигнутой целевой функции. Качество оптимизации с применением дан-ной модификации, как и в случае с классическим вариантом генетического алгоритма, увеличи-вается с увеличением шага мутации и количества элементов в популяции. На рис. 2 приведен график истории целевой функции при данной модификации в сравне-нии с генетическим алгоритмом без модификаций и с алгоритмом направленного поиска. Рис. 2. История оптимизации гидравлических сопротивлений клапанов отопительных приборов генетическим алгоритмом с модификацией «Направленная мутация» в сравнении с генетическим алгоритмом без модификации и с направленным поиском Для определения экономической эффективности данной разработки сравним потребление тепловой энергии до и после оптимизации. Данные для сравнения сведены в табл. 9. Таблица 9Сравнение потребления тепловой энергии до оптимизации и после оптимизацииЧто оптимизировалиШаг мутацииРазмер популяцииПотребление тепловой энергии, Гкал/чминимальныймаксимальныйначальныйСтояки501550200,110БатареиБез оптимизации0,126Из таблицы видно, что потребление энергии снизилось примерно на 14 %, или на 0,016 Гкал/ч в абсолютном значении с одного дома. По данным Федеральной службы государственной статистики, на 2018 г. число многоквартирных домов составляет примерно 2,7 млн [14]. Таким образом, в масштабах страны в отопительный сезон экономия тепловой энергии составляет 43 200 Гкал/ч. Стоимость 1 Гкал тепла примерно 2 000 руб. Итого в масштабах страны экономия составляет примерно 360 млн руб. только за один час во время отопительного сезона.ЗаключениеВ ходе эксперимента мы пришли к следующим выводам:1. Генетический алгоритм позволяет добиться наилучших результатов при наибольшем значении максимального размера популяции, начального размера популяции и шага мутации. 2. Минимальный размер популяции практически не влияет на качество оптимизации.3. Так как значительное увеличение количества элементов в популяции невозможно, необходимо провести исследование влияния на качество оптимизации больших значений шага мутации. 4. Введение модификации «Направленная мутация» позволяет увеличить глубину оптимизации.В результате проведенной оптимизации гидравлических сопротивлений клапанов отопи-тельных приборов получены благоприятные эффекты:– уменьшение тепловых выбросов в атмосферу (экологический эффект);– благоприятный температурный режим помещения (оздоровительный эффект);– уменьшение перерасхода тепловой энергии примерно на 14 % (экономический эффект).