Full Text

Введение Тепловое потребление является одной из основных статей топливно-энергетического баланса нашей страны. По данным Единой межведомственной информационно-статистической системы (ЕМИСС) [1], за 2020 г. расходы топливно-энергетических ресурсов на производство пара и горячей воды, отпущенной котельными, составили более 40 млн т условного топлива. Чтобы организовать теплофикацию, необходимо наличие следующей инфраструктуры [2]: · Источник генерации тепла (котельная, ТЭЦ, АТЭС). · Разветвленная сеть теплоснабжения. Представляет собой теплоизолированные теплопроводы, по которым осуществляется доставка теплоносителя абонентам (конечным потребителям тепла). · Тепловые пункты (ТП). Назначение ТП - подготовка теплоносителя согласно температурным графикам для конечных потребителей, происходит с помощью сетевой воды (пара). · Оборудование конечных потребителей. Предназначено для использования горячей воды (пара) в системах отопления, горячего водоснабжения (ГВС), вентиляции, кондиционирования воздуха. Температурный график - это неэмпирическая зависимость температуры сетевой воды от температуры наружного воздуха [3]. Значения температурных графиков принято рассчитывать по формулам. В основе этих формул заложены уравнения теплопередачи. Результаты расчетов обычно представляются в табличной форме. 1. Особенности теплоснабжения Основной задачей теплоснабжающей организации (ТСО) является поставка теплоносителя в многоквартирные дома (МКД) с соблюдением установленного нормативного температурного графика. Допускается отклонение в размере 5 % по температуре, поступающей в МКД, в меньшую сторону от температуры, обозначенной в температурном графике. Следовательно, актуальной задачей для ТСО является задача качественного (оптимального) регулирования тепловой нагрузки на котельной. В качестве базовой теплосети рассматривается закрытая сеть теплоснабжения, упрощенная структура которой представлена на рис. 1. Основными элементами, образующими теплосеть, являются: · Котельная - поставщик теплоносителя. · ТК - тепловые камеры. · МКД - многоквартирные дома, потребители теплоносителя. · Соединительные участки теплосети - трубы, соединяющие основные узлы теплосети. Регулирование на котельной, которая питает тепловую сеть с неизвестными теплофизическими характеристиками, представляет собой задачу, относящуюся к классу слабоформализуемых проблем. Это не в последнюю очередь связано с характерными особенностями эксплуатации сетей теплоснабжения: · инерционность - теплоноситель с измененными характеристиками поступает к потребителю через определенный промежуток времени; · износ сети теплоснабжения, который может быть как прогнозируемым, так и непрогнозируемым; · сложность топологии сети теплоснабжения, которая оказывает влияние на распределение и изменение характеристик теплоносителя; · неполнота данных с общедомовых приборов учета (ОДПУ), порой снимается только часть данных (например, только температура) или вообще ОДПУ может отсутствовать. Рис. 1. Упрощенная структура закрытой сети теплоснабжения В случае полноты данных по обозначенным выше характерным особенностям при разных температурах окружающей среды и в разные моменты времени задача относительно легко может быть решена. К сожалению, сбор такой информации практически невозможен. Приходится ставить и решать задачу в условиях неопределенности. Ввиду этого было решено разработать виртуальный стенд городской сети теплоснабжения в условиях неполноты данных, который позволит создавать модели сетей теплоснабжения с разнообразными условиями тепловых потерь, с потребителями тепла разных категорий и разным потреблением тепловой энергии, а также моделировать вносимые изменения без внесения фактических изменений в структуру городской теплосети. Таким образом, в данной работе объектом исследования является закрытая городская сеть теплоснабжения со сложной топологией, а предметом - расчетные модели теплосети, позволяющие определять потери температуры теплоносителя на участках теплосети на основе обработки данных с ОДПУ. Цель работы - повышение эффективности управления городской сетью теплоснабжения на основе данных, получаемых с помощью виртуального стенда. Для достижения обозначенной цели необходимо решить следующий перечень задач: 1. Сформулировать функциональные требования к виртуальному стенду и реализовать в виде конечного программного продукта. 2. Разработать алгоритм, позволяющий вычислять параметры сети теплоснабжения на основе данных с ОДПУ. 3. Провести апробацию стенда на реальной сети теплоснабжения: · построить топологию сети теплоснабжения; · вычислить параметры сети теплоснабжения на основе обработки данных с ОДПУ; · с учетом заданных и вычисленных параметров сети теплоснабжения провести ряд численных экспериментов, направленных на определение отклонений значений температуры по МКД, получаемых на стенде, от значений температуры с ОДПУ (фактических). 2. Расчет потерь тепла Расчет теплопотерь на участках теплосети осуществляется по следующей формуле [4]: (1) где τ1, τ2 - температуры теплоносителя в начале и в конце участка, °С; t0 - температура окружающей среды, °С; l - длина теплопровода, м; β - коэффициент местных потерь теплоты, зависит от диаметра труб[1]; R - суммарное термическое сопротивление тепловой изоляции, стенки канала и грунта, (м·°С)/Вт; G - расход теплоносителя, кг/с; - теплоемкость теплоносителя, Дж/(кг·°С). Неизвестным параметром в формуле (1) является суммарное термическое сопротивление, которое сложно точно определить по приведенным в источниках формулам [5-10]. Проблема определения точного значения суммарного термического сопротивления связана с неравномерным для всего участка трубопровода ухудшением изоляционных характеристик как самого теплоизоляционного материала, так и стенок канала, в том числе грунта [11]. Таким образом, необходимо реализовать в рамках виртуального стенда сети теплоснабжения алгоритм вычисления суммарного термического сопротивления в зависимости от топологии и показаний ОДПУ, с учетом погодных условий и характеристик теплоносителя на выходе из котельной. 3. Вычисление значений суммарных термических сопротивлений Для того чтобы определить значения суммарных термических сопротивлений на каждом участке теплотрассы, необходимо решить задачу оптимизации с заданными ограничениями. В качестве минимизируемой была выбрана функция следующего вида: где - фактически замеренная температура на входе в i-й МКД; - вычисленное на виртуальном стенде значение температуры на входе в i-й МКД; N - количество МКД. При этом вводятся ограничения по значениям суммарных термических сопротивлений: где - суммарное термическое сопротивление тепловой изоляции, стенки канала и грунта на участке теплосети, ведущем к i-му узлу теплосети от узла, с которым соединен; - количество узлов (ТК и МКД) теплосети. Для учета обозначенных ограничений сформированы следующие функции: (2) (3) (4) (5) (6) (7) (8) где λ - коэффициент, определяющий вид штрафной функции. С учетом формул (2)-(8) общий вид штрафной функции следующий: Таким образом, минимизируемая функция с учетом ограничений примет следующий вид: Для поиска минимума функции решено было использовать метод Нелдера - Мида (метод деформируемого многогранника [12, 13]). Основная идея метода заключается в последовательном перемещении и деформировании симплекса в районе точки экстремума. Метод находит локальный экстремум и может «застрять» в одном из них. В случае когда требуется найти глобальный экстремум, можно пробовать выбирать другой начальный симплекс [12, 13]. В реализации поиска глобального экстремума в рамках виртуального стенда был разработан алгоритм, в рамках которого происходит перестроение симплекса из текущей найденной точки экстремума путем увеличения граней симплекса на определенно заданную величину. Процедура поиска глобального экстремума происходит до тех пор, пока не будет найдена новая точка, где значение минимизируемой функции меньше наперед заданного значения, или пока новая точка экстремума не совпадет с точкой, найденной на предыдущем шаге. Кроме того, процедура поиска глобального экстремума ограничена количеством попыток поиска. 4. Функциональные требования Для реализации будущего виртуального стенда городской сети теплоснабжения были сформулированы следующие функциональные требования: 1. Генерация структуры упрощенной сети теплоснабжения случайно или вручную из набора заданных параметров/элементов. 2. Сохранение в файл и загрузка из файла конфигурации и текущего состояния сети теплоснабжения. 3. Управление параметрами сети теплоснабжения. 4. Расчет сети теплоснабжения: · расчет состояния теплосети через заданный временной интервал; · расчет установившегося состояния. 5. Определение среднего суммарного термического сопротивления тепловой изоляции для каждого МКД на основе: · данных по температуре и расходу теплоносителя на выходе из бойлерной; · данных по температуре теплоносителя на входе в МКД; · данных по температуре окружающей среды. 6. Визуализация расчетных результатов. 7. Имитация распределения температуры внутри сети теплоснабжения при изменяющихся параметрах: погодные условия, характеристики теплоносителя, поступающего с котельной. 5. Виртуальный стенд сети теплоснабжения В рамках виртуального стенда реализована возможность генерации случайной топологии сети. Случайная генерация происходит на основе следующих параметров: · Общая протяженность теплосети. Суммарная длина участков теплосети от котельной до каждого узла теплосети (ТК, МКД). · Количество ТК. Общее количество в генерируемой теплосети тепловых камер. Значение должно быть больше 1. · Количество МКД. Общее количество в генерируемой теплосети тепловых камер. Значение должно быть больше 1. · Минимальное расстояние между ТК. Минимально допустимое расстояние (длина трубопровода) от ТК до узла, к которому он подключен. · Максимальное/минимальное расстояние от ТК до МКД. Максимально/минимально допустимое расстояние (длина трубопровода) от МКД до ТК, к которому он подключен. · Максимальное суммарное термическое сопротивление. Максимально возможное допустимое значение суммарного термического сопротивления должно быть обязательно больше 1Е-4 - нижняя граница, используемая при генерации параметров теплосети. При случайной генерации теплосети учтены следующие моменты: 1. Диаметр труб, распределяющих теплоноситель от ТК, должен быть меньше диаметра трубы, подводящей теплоноситель к данному ТК. 2. Суммарное максимально возможное расстояние до каждого узла теплосети (ТК, МКД) должно быть меньше заданной общей протяженности сети. В случае некорректного ввода данных выводится сообщение с рекомендацией изменить один из заданных параметров. В случае ручного создания теплосети необходимо добавлять каждый узел в отдельности и задавать минимально необходимые параметры узла: диаметр и длину трубы, проводящей теплоноситель, а также величину суммарных термических потерь, а для МКД обязательным параметром является величина расхода (рис. 2). Сгенерированная теплосеть отображается в виде древовидной структуры, пример которой представлен на рис. 3. Для работы с ранее созданными моделями теплосетей в рамках виртуального стенда реализована возможность сохранения в файл и загрузки конфигурации и текущего состояния теплосети. Данные по теплосети в файлах представлены в формате JSON [14]. Предусмотрены следующие виды расчетов: · Расчет состояния теплосети через заданный временной интервал. С учетом скорости движения теплоносителя производится расчет распределения теплоносителя в теплосети. · Расчет установившегося состояния. Расчет распределения теплоносителя в теплосети без учета временной составляющей. Рис. 2. Генерация топологии теплосети вручную Рис. 3. Отображение топологии теплосети Для всех обозначенных расчетов необходимо задать исходные данные - погодные условия (рис. 4) и характеристики теплоносителя на выходе с котельной (рис. 5), остальные характеристики теплоносителя, необходимые для расчета, берутся из справочной литературы [15]. Рис. 4. Форма для задания погодных условий Рис. 5. Форма для задания характеристики теплоносителя на выходе с котельной 6. Апробация виртуального стенда В качестве объекта исследования была взята одна из действующих теплосетей закрытого типа г. Лысьвы (рис. 6). Теплосеть состоит из одной котельной и подключенных к ней 25 потребителей тепла (МКД) через соединительные узлы теплосети (ТК) в количестве 26 штук. Следовательно, для построения полнофункциональной модели данной теплосети нужно определить значения суммарного термического сопротивления для 51 участка. Рис. 6. Рассматриваемая теплосеть Апробация работы стенда была проведена на базе данных, представленных ООО «СофтМ», по теплосети г. Лысьвы Пермского края. Вычисление суммарных термических сопротивлений предполагалось проводить на основе данных, предоставляемых с ОДПУ. Однако не во всех домах, подключенных к теплосети, установлены ОДПУ. Исходя из этого, ввиду отсутствия данных по ряду МКД было решено задать примерные данные по расходу, схожие с аналогичными у МКД. В качестве снимаемых значений температуры теплоносителя задавалось нормативное значение из температурного графика с учетом допустимого отклонения не более 5 %. Были проведены расчеты по определению суммарных термических сопротивлений на основе осредненных данных, полученных 05.02.2019 г. за период времени с 1:06 до 1:36 ч. Данный временной интервал выбран с учетом инерционности сети, составляющей 18 мин, при учете неизменности погодных условий в заданный промежуток времени и до начала периода осреднения. На основе вычисленных значений суммарных термических сопротивлений были проведены расчеты по определению температуры в МКД с учетом усредненных данных, полученных 05.03.2019 г. с 6:40 до 8:20 ч. Временной интервал был также выбран с учетом вышесказанных соображений. Для сравнительного анализа результаты расчета сведены в таблицу. Сравнение результатов МКД Температура с ОДПУ Температура вычисленная Разница Градусы Проценты Гайдара, 1 77,9625 78,233 434 09 2,27 2,91 Гайдара, 22 76,511 25 78,538 450 9 2,03 2,65 Гайдара, 26 76,494 285 71 78,654 557 06 2,16 2,82 Гайдара, 28 75,631 25 77,846 215 7 2,21 2,93 Гайдара, 30 76,195 78,122 501 53 1,93 2,53 Гайдара, 32 75,988 75 77,897 294 94 1,91 2,51 Куйбышева, 8 76,271 25 78,261 401 09 1,99 2,61 Куйбышева, 10 76,566 25 78,742 516 93 2,18 2,84 Куйбышева, 14 76,508 75 78,802 990 73 2,29 3,00 Куйбышева, 16 76,66 79,057 376 49 2,40 3,13 Энгельса, 18 76,48 78,600 364 67 2,12 2,77 Энгельса, 22 75,546 25 78,573 796 12 3,03 4,01 Энгельса, 23 76,507 5 78,876 773 82 2,37 3,10 Энгельса, 25 76,275 78,554 994 16 2,28 2,99 Энгельса, 29 76,467 142 86 78,432 361 12 1,97 2,57 Энгельса, 33 76,48625 78,652 376 28 2,17 2,83 Энгельса, 35 76,41625 79,008 915 39 2,59 3,39 Энгельса, 37 75,75 78,238 607 13 2,49 3,29 Как видно из таблицы, наибольшее отклонение вычисленных значений температуры от фактически замеренных составило 4,01, или 3,03 °С. Данные показатели говорят о достоверности созданной виртуальной модели теплосети в рамках допустимых ограничений в 5 %. Заключение Таким образом, в рамках выполненной работы реализован виртуальный стенд городской сети теплоснабжения в условиях неполноты данных с обозначенным функционалом. Разработан алгоритм поиска суммарных термических сопротивлений на основе метода Нелдера - Мида. Проведена верификация получаемых значений с помощью созданного виртуального стенда с данными, получаемыми с ОДПУ. Погрешность вычислений по виртуальному стенду составила менее 5 %. Созданный стенд позволяет прогнозировать поведение сети теплоснабжения, что положительным образом может сказаться на эффективности управления конкретной городской сетью теплоснабжения. В дальнейшем планируется проведение дополнительных численных экспериментов и усовершенствование методики определения суммарных термических потерь с целью снижения погрешности вычислений.

About the authors

R. Yu Bannikov

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 96

PDF (Russian) - 57

Refbacks

• There are currently no refbacks.