УДК 536.7 + 519.85 ББК 78.34
МЕТОДЫ ТЕПЛОГИДРАВЛИЧЕСКОЙ ОПТИМИЗАЦИИ И УПРАВЛЕНИЯ ТЕПЛОВЫМИ СЕТЯМИ
Майков И. Л.1, Директор Л. Б.2, Зайченко В. М.3
(Учреждение Российской академии наук Объединенный Институт высоких температур РАН,
Москва)
Представлен метод решения одной из практически важных задач, связанных с эффективным теплоснабжением абонентов тепловых сетей, как в случае сильной гидравлической разбалансировки сети (статическая задача), так и при изменениях тепловых нагрузок в процессе эксплуатации (динамическая задача). Разработан численный метод для анализа и оптимизации тепловых сетей, позволяющий решать задачи оптимизации и управления тепловой сетью с большим количеством абонентов с целью реализации оптимального режима теплопотребления. Показано, что и статическая, и динамическая задачи могут быть решены в рамках одной общей нелинейной оптимизационной задачи. Разработан алгоритм решения такой задачи, в котором минимизация целевой функции осуществляется в итерационном процессе с использованием метода покоординатного спуска.
Ключевые слова: математическое моделирование, нелинейная оптимизация, задачи минимизации, численные методы, тепловые сети.
1 Игорь Леонидович Майков, доктор физико-математических наук, в.н.с. ([email protected]).
2 Леонид Бенцианович Директор, доктор технических наук, в.н.с. ([email protected]).
3 Виктор Михайлович Зайченко, доктор технических наук, зав. лаб. ([email protected]).
1. Введение
Тепловые сети предприятий или городских микрорайонов, обеспечивающие подачу тепловой энергии конкретным потребителям (абонентам) для целей отопления, являются достаточно сложными теплогидравлическими системами, включающими, помимо собственно трубопроводов различного диаметра, регулирующую арматуру, насосы, смесительные устройства, отопительные приборы. При проектировании жилых зданий, промышленных объектов тепловая сеть рассчитывается таким образом, чтобы ее гидравлические характеристики обеспечивали заданную температуру воздуха внутри помещений в широком диапазоне сезонного изменения температуры наружного воздуха при качественном регулировании (регулировании соответствующим изменением температуры теплоносителя) [5].
В процессе развития или перепрофилирования предприятий их тепловые сети претерпевают значительные изменения и тепловые нагрузки, заложенные в первоначальный проект системы теплоснабжения, не всегда соответствуют фактическим потребностям абонентов тепловой сети. Кроме того, в процессе длительной эксплуатации значительно изменяются гидравлические характеристики сети. Возникающие при этом проблемы, связанные с недостаточным обеспечением теплом отдельных абонентов, решаются наименее эффективно - увеличением тепловой нагрузки всей тепловой сети.
Настройка тепловой сети, насчитывающей десятки абонентов, представляет собой сложную, многопараметрическую задачу, решить которую без наличия адекватной математической модели сети практически нереально.
В работе представлен метод решения одной из практически важных задач, связанных с эффективным теплоснабжением абонентов тепловых сетей, как в случае сильной гидравлической разбалансировки сети, так и при изменениях тепловых нагрузок в процессе эксплуатации. В основу теплогидравлической модели тепловой сети положена рекурсивная схема [3]. Модель сети строится «сверху вниз»: сначала в качестве основных элементов рассматриваются крупные участки сети, затем
они детализируются, впоследствии детализируются их абоненты и т.д. Таким образом, с помощью одной и той же программы может быть смоделирована сеть практически любого уровня сложности, и ее конфигурацию можно изменять в широких пределах.
В работе сформулировано несколько нелинейных задач оптимизации с ограничениями (статическая и динамическая задачи) с соответствующими критериями оптимизации, определены условия разрешимости этих задач. Сформулирована общая задачи оптимизации теплогидравлической сети, имеющая однозначное решение.
Разработан численный метод для анализа и оптимизации тепловых сетей, позволяющий решать задачи оптимизации и управления сетью с большим количеством абонентов с целью реализации оптимального режима теплопотребления.
2. Постановка задачи
При построении модели рассматривается формализованная схема теплотрассы с присоединенными к ней абонентами (рис. 1). Такой формальный подход позволяет описывать схемы тепловых сетей практически любой сложности и глубины детализации.
Обозначения на схеме: абонент (А) - подключенный между подающей и обратной магистралью элемент тепловой сети, в котором происходит полезное тепловыделение;
коллектор (К1, К2) - часть тепловой сети, содержащая параллельно включенные элементы (абоненты, коллекторы, трассы) с общими точками на входе и выходе;
трасса (Т1, Т2) - часть тепловой сети, включающая набор коллекторов, абонентов и трасс;
узел (У) - точки на прямой и обратной магистралях, между которыми подключен абонент (группа абонентов);
дроссельное устройство (Д) - регулирующий клапан.
К2
К1
У
Т2 -->
гтт
Д-
Рис. 1. Схема элементов теплотрассы
Исходными данными для расчета теплотрассы служат: проектные тепловые нагрузки абонентов Qi и соответствующие им расходы теплоносителя (воды) Gi (при заданных температурах в помещениях абонентов Д/)й); температура наружного воздуха Та; температура Tin и расход воды Gin на входе в теплотрассу; параметры участков теплотрассы и абонентских вводов; потери давления на абонентских вводах и выводах; потери давления на абонентах. Индекс i = 1, N, где N - количество абонентов.
Расчетные гидродинамические сопротивления на абонентах вычисляются из уравнения Бернулли dp(i)
Dz (і) =
(1)
Dzeta = ^ Dz (і).
В уравнениях (1) под dp(i) понимаются потери давления: DP1i - потери давления на абоненте (между точками 2 и 3); DP2i - потери давления на абонентском вводе (между точками 1
2
г =1
и 2); DP3i - потери давления на абонентском выводе (между точками 3 и 4); DP4i - потери давления на дроссельном устройстве (между точками 4 и 5); G(i) -расход теплоносителя через абонент i; р - плотность теплоносителя, Dzeta - общее сопротивление между точками 1 и 5.
В дальнейшем, не теряя общности, пренебрегаем потерями тепла в теплотрассе (их учет может быть произведен для любой конкретной задачи [5]). Температура воды на выходе из абонента вычисляется из условия теплового баланса
(2) Ti Ti_________'
(2) out in ss-^1 \ ,
(Gtcp)
где cp - теплоемкость теплоносителя.
Для узловой точки j используется средневзвешенное значение температуры
G Ti + G T і-1
(3) T j _ i out i-1 out
v s out '
Gi + Gi-1
Алгоритм расчета.
1. По располагаемым напорам на абонентах и потерям напора на абонентских вводах по уравнению (1) определяются их коэффициенты гидравлического сопротивления. Все DP4i предполагаются равными 0.
2. Рассчитывается температура воды на абонентах и во всех узловых точках теплотрассы по формулам (2) и (3).
Пример расчетной схемы тепловой сети с разбивкой на расчетные элементы в соответствии с рис. 1 приведен на рис. 2.
Полученные поля температур и расходов при этом будут обеспечивать работу теплотрассы в проектном режиме, т.е. все характеристики теплотрассы согласованы таким образом, что обеспечивается заданная температура воды на выходе из абонентов и, соответственно, заданная температура воздуха в помещениях абонентов. Тепловая сеть состоит из одной трассы, 25 коллекторов и 45 абонентов.
а)
-СХ-| 26 |-СХ5-
ЧХЪ| 27 |чхь
-Сх}-| 28 ІЧХ-
НХ-| 29
ЧХ-| 30 |-0<3-
чхь| 31 |чхз-
-Чх1-| 32 Нхь-
-СХ-| 33 Нхн
-СХ-| 34 |-СХ]-
ЧХЪ| 35 |чх-
-СХ-| 36 1—00-
-00-| 37
б) Трасса
Коллекгор1
Абонент1 Коллектор11
Абоншт2 Коллектор 12
АбоненгЗ
Абоненг4
Коллектор13
АбонентЗ
Абонентб
Абоненг7
Коллектор2
Абонент8 Коллектор14
Абонеит9
АбонштЮ
КоллекторЗ
Коллектор25
Абонент 11 Абонеит12 Абоншт13 Коллектор15
Коллектор 16
Абонент 14 Абонент 15 Коллектор 17
Абоненг16
Коллектор18
Абонент17 Коллектор 19
Абонент 18 Абонент 19
Коллектор4
Абонент20
Коллектор20
Абоншт21 Коллектор 21
Абонент22
Абонент23
Абоненг24
Абонент25
Коллектор5
Абонент26
Абонент27
Абонент28
Коллекторб
Абонент29 АбонентЗ О АбоненгЗ 1 Коллектор7
АбонентЗ 2 Коллектор22
АбонштЗЗ
Абоншт34
Коллектор8
АбонентЗ 5 АбонентЗ 6 АбонентЗ 7 Коллектор9
АбонентЗ 8 Коллектор23
Абонент39 Коллектор 24
Абоненг40
Абоненг41
КоллекторЮ
Абонент42
Абонент43
Абонент44
Абонент45
Рис. 2. Тепловая сеть: а) - расчетная схема, б) - разбивка на расчетные элементы
Результаты расчетов (выходная температура воды на абонентах) приведены на рис. 3.
Номер абонента
Рис. 3. Распределение температуры теплоносителя на выходе из абонентов
3. Процедуры оптимизации
3.1. СТАТИЧЕСКАЯ ЗАДАЧА
Рассмотрим случай, когда в результате расчетов (рис. 3) не удалось получить согласованные температуры обратной воды, соответствующие заданной температуре наружного воздуха. Это может быть связано с гидродинамической рассогласованностью сети.
Задача 1. Требуется найти такие DP4Z, чтобы температура воды на выходе из каждого абонента была равна Таи, = 70 0С (расчетная температура воды в обратной магистрали при расчетной температуре воздуха для проектирования. При этом температура воды в прямой магистрали равна Tin = 130 0С). Расход теплоносителя в тепловой сети согласован с теплопо-треблением Q, т.е.
N
IQ
(4) G =----^-------,
Ср (Тт - Taut )
где О - расход теплоносителя в теплосети; ср - теплоемкость теплоносителя. Все 0Р4г- > 0, так как
где 0г4г- - сопротивления дроссельных (регулирующих) устройств. Из технических соображений сопротивления дроссельных устройств ограничиваются и сверху, что приводит к ограничениям и потерь давления на дроссельных устройствах 0Р4г- < 0Р4тах, где 0Р4тах - максимальный перепад давления на дроссельных устройствах.
Рассмотрим функцию
Функция (6) является функцией N переменных - перепадов давления на дроссельных устройствах 0Р4Ь стоящих на каждом абоненте. В выражении (6) каждая выходная температура учитывается со своим весом ОI согласованно с (3).
3.2. МЕТОД РЕШЕНИЯ
Сформулируем следующую оптимизационную задачу: найти такие 0Р4Ь при которых целевая функция (6) при Тоиг = 70 0С достигает минимума:
причем 0 < 0Р4г- < 0Р4тах и выполняется условие (4).
Оптимизационная задача (7) является нелинейной задачей оптимизации с ограничениями 0 < 0Р4г- < 0Р4тах и с размерностью N, причем N велико. Большинство методов решения таких задач требует вычисления первых или даже вторых производных [2, 4]. На практике такие вычисления часто оказываются слишком трудоемкими либо их принципиально невозможно осуществить точно по причине отсутствия явных выражений не только для производных (как и в рассматриваемой задаче), но и для самой целевой функции (6). Наиболее перспективным методом, не связанным с (точным или приближенным) вычислением производных, является метод покоординатного спуска [2], также называемый еще релаксационным методом [1].
N
(6) L(DP41,..,DP4N) = ^0\ГШ -Тоиі
І=1
(7) L(DP41,..,DP4N) ® тіп,
Общая схема метода включает следующие шаги.
1. Задание нулевого приближения (к = 0 - номер приближения) БР4(к} = (ОР 4(к},.., ОР 4(к}).
2. Рассмотрение функции ¿(0Р4(к),.., 0Р4^)) при фиксированных 0Р4(2к),.., 0Р4N) как функции одной переменной 0Р4Ь
3. Нахождение минимума функции одной переменной Ь(0Р41,0Р4(2к}..,0Р4(к}), определение 0Р4(к+1).
4. Повторение пунктов 2 и 3 для остальных переменных 7 = 2, ..N для вычисления БР4(к+1).
5. Проверка условия ^(БР4(к+1))-ь(рР4('к) )< е (е - заданная
точность). Если условие выполняется - выход, в противном случае к = к + 1 и переход к п. 2.
Центральным звеном предложенного алгоритма является поиск минимума функции одной переменной (п. 3). В работе используется алгоритм Брента [6], в основе которого лежит комбинация методов золотого сечения и последовательной параболической интерполяции.
Результаты расчетов 0Р47 представлены на рис. 4, распределение температуры теплоносителя на выходе из абонентов -на рис. 5.
Максимальное отклонение расчетной температуры воды от заданной составляет менее 0,3 0С (рис. 5). При отсутствии регулирования отклонение температур составляло более 10 0С (рис. 3).
Таким образом, дополнительные сопротивления (сопротивления дроссельных устройств) на абонентах (рис. 4) обеспечивают согласованный (проектный) режим работы теплосети - все абоненты имеют проектную температуру воздуха внутри помещения (расчеты проводились при 7(7% = 16 0С).
Номер абонента
Рис. 4. Распределение перепадов давления на дроссельных устройствах абонентов
Номер абонента
Рис. 5. Распределение температуры теплоносителя на выходе из абонентов после процедуры оптимизации
3.3. ДИНАМИЧЕСКАЯ ЗАДАЧА
Рассмотрим динамическую задачу, когда произошло изменение тепловой нагрузки одного или нескольких абонентов.
Задача 2. Требуется найти DP4{ для каждого абонента и расход воды G в теплотрассе, обеспечивающие температуру воды на выходах всех абонентах Tout = 70 0С.
Оптимизационная задача (7) в стандартной постановке не дает решения поставленной задачи, так как не выполняется условие (4). Но задача (7) имеет сходящееся решение, если определить выходную температуру воды следующим образом
N
IGTL
(8) ToUt = ^n-----.
I Gi
i=1
Оптимизационная задача может быть сформулирована следующим образом:
найти такие DP4i, при которых целевая функция (6) достигает минимума:
(9) L( DP 41,., DP 4 N) ® min,
причем 0 < DP4i < DP4max и ToUt удовлетворяет условию (8).
Решение задачи (9) определяет DP4i для каждого абонента такие, что температура обратной воды ToUt удовлетворяет условию (8).
Рассмотрим функцию
N
(10) L(DP4i,..,DP4n,G) = IGiTiut -Tj .
i=l
Тогда общая оптимизационная задача формулируется следующим образом
найти такие DP4i, при которых целевая функция (10) при ToUt = 70 0C достигает минимума:
(11) L(DP41,..,DP4N,G) ® min,
причем 0 < DP4i < DP4max и выполняется условие (4).
Следует заметить, что задача (7) является частным случаем общей задачи (11) (при фиксированном G задача (11) переходит в (7)).
Минимизация функционала (10) осуществляется аналогично схеме, описанной в п. 3.2, с учетом, что вектор DP4 = (DP41, ... , DP4n, G) имеет размерность N + 1.
3.4. ЗАМЕЧАНИЕ
Рассмотрим случай, когда условие (4) выполняется, но решение задачи (7) не существует, т.е. невозможно получить Tout = 70 0C одновременно для всех абонентов. То есть можно получить только решение задачи (9), в котором Tout Ф 70 0C, а удовлетворяет условию (8). Математически это объясняется ограничениями на DP4i. По условиям задачи (7) они должны удовлетворять условию неотрицательности, и поэтому, достигнув нулевого значения, DP4t не меняются, при этом значение целевой функции (6) имеет некоторое конечное значение. Если снять условие неотрицательности, получится сходящееся решение задачи (7), и значение целевой функции будет стремиться к нулю. Отрицательные значения DP4t на дроссельных устройствах говорят о том, что для обеспечения оптимального режима необходимо уменьшить сопротивления самих абонентов. В постановке задачи сопротивления абонентов заданные величины (либо заданы потери напора на абонентах, связанные с сопротивлениями соотношением (1)) и меняться не могут.
Таким образом, даже выполнение условия (4) не гарантирует возможности получения решения Tout = 70 0C одновременно для всех абонентов без изменения расхода теплоносителя в тепловой сети, а значит необходимо решение задачи (11).
Общий алгоритм решения задач 1 и 2 состоит из следующих шагов:
1. M = 1 (M - номер внешней итерации).
2. m = 1 (m - номер внутренней итерации).
3. Задание Tout = 70 0C.
4. Решение задачи (9) (определение DP4, i = 1, ..., N).
5. Расчет нового значения Tout согласно (8).
6. Проверка условия сходимости (по изменениям L и DP4i и в итерациях).
7. При выполнении условий шага 6 - переход к п. 8, в противном случае - переход к п. 4 (задание m = m + 1).
8. Проверка условия Т(ли = 70 0С?
9. При выполнении условия шага 8 - выход, в противном случае - задание Тои1 = 70 0С и переход к п. 10.
10. Решение задачи (11) при фиксированных ВР4Ь г = 1, ..., N (определение О).
11. Переход к п. 2 (задание М = М + 1).
Итерации внутреннего цикла дают решение задачи (9). Выполнение условия п. 8 при М = 1 дает решение задачи (7) (задача 1 без корректировки расхода). Выполнение условия п. 8 при М > 1 (итерации внешнего цикла) дает решение задачи (11) (задача 1 с корректировкой расхода или динамическая задача -задача 2).
4. Результаты
Расчеты проводились для реальной тепловой сети одного из московских предприятий. Расчетная схема тепловой сети представлена на рис. 2, расчетная тепловая нагрузка теплосети Q = 25,497 Гкал/ч, расчетный расход теплоносителя в теплосети О = 424,95 т/ч (тепловая нагрузка и расход удовлетворяют условию (4)), Тгп = 130 0С.
Начальное распределение температуры воды на выходе из абонентов и результаты расчетов представлены на рис. 6.
Результаты расчетов показывают, что, несмотря на выполнение условия (4), сходящееся решение (решение задачи (9) получается только при Тои1 = 68,345 0С (кривая 2). Корректировка расхода (решение задачи (11)) приводит к Тои1 = 69,983 0С, при этом О = 438,96 т/ч.
Следующий пример показывает решение динамической задачи. Пусть по каким-то причинам (ухудшение теплоизоляционных свойств ограждающих конструкций, рост или снижение тепловыделения внутри зданий и т.п.) изменилось теплопотреб-ление абонентов.
Соответственно изменились температуры теплоносителя на выходах абонентов. Зная расходы и разности температур по формуле (2) можно определить новые Qi и решить задачу (11).
Исходные данные модельной задачи - решение предыдущей задачи (вычисленные новый расход и перепады давления на дроссельных устройствах, которые обеспечивают Tout = 69,983 0C ). Кроме того, введем изменения в абоненте №32 (рис. 2): тепловая нагрузка выросла в 2 раза и составила величину 11,02 Гкал/ч. Сходящееся решение при заданном расходе теплоносителя получилось при температуре воды в обратной магистрали ТоШ = 57,359 °С.
Номер абонента
Рис. 6. Распределение температуры теплоносителя на выходе из абонентов: 1 - исходное распределение;
2 - решение задачи (9); 3 - решение задачи (11)
Результаты решения задачи оптимизации (11) представлены на рис. 7. Новый расход составил 533,438 т/ч, температура воды в обратной магистрали Tout = 70,029 0C.
Как видно из рис. 7, на изменение нагрузки одного абонента отреагировала вся система: произошло изменение сопротивления всех дроссельных устройствах, максимальное изменение потери давления составило более 2 м вод. ст.
14-1
Номер абонента
Рис. 7. Распределение перепадов давления на дроссельных устройствах абонентов:1 - исходное распределение;
2 - реакция системы на изменение теплопотребления
5. Выводы
В работе показана возможность расчетной оптимизации тепловой сети как в случае сильной гидравлической разбаланси-ровки сети, так и при изменениях тепловых нагрузок в процессе эксплуатации. Показано, что даже согласованный расход теплоносителя в тепловой сети с теплопотреблением не гарантирует получения оптимального одновременно для всех абонентов без изменения расхода теплоносителя в тепловой сети, и необходимо решение общей оптимизационной задачи.
Разработанная оптимизационная процедура позволяет проводить расчет характеристик регулирующих элементов и приводить тепловую сеть к оптимальным (проектным) параметрам, а также отслеживать изменение параметров тепловой сети и
возвращать их в оптимальное состояние за счет изменения
характеристик регуляторов и общего расхода.
Литература
1. АОКИ М. Введение в методы оптимизации. - М.: Наука, 1977. - 344 с.
2. БОГЛАЕВ Ю.П. Вычислительная математика и программирование. - М.: Высшая школа, 1990. - 544 с.
3. ДИРЕКТОР Л.Б., ЗАЙЧЕНКО В.М., МАЙКОВ И.Л., ФРИД С.Е. Математическое моделирование и оптимизация разветвленных тепловых сетей // Новости теплоснабжения. - 2002. - №5(21). - С. 36-38.
4. ИЗМАИЛОВ А.Ф., СОЛОДОВ М.В. Численные методы оптимизации: Учеб. пособие. - М.: ФИЗМАТЛИТ, 2005. -304 с.
5. СОКОЛОВ Е.Л. Теплофикация и тепловые сети: Учебник для ВУЗов. - 5-е изд. - М.: Энергоатомиздат, 1982. - 360 с.
6. ФОРСАЙТ Д., МАЛЬКОЛЬМ М., МОУЛЕР К. Машинные методы математических вычислений. - М.: Мир, 1980. -280 с.
METHODS OF THERMO-HYDRAULIC OPTIMIZATION AND MANAGEMENT OF HEAT SUPPLY NETWORKS
Igor Maikov, Joint Institute for High Temperatures of RAS, Moscow, professor ([email protected]).
Leonid Director, Joint Institute for High Temperatures of RAS, Moscow, professor ([email protected]).
Victor Zaitchenko, Joint Institute for High Temperatures of RAS, Moscow, professor, head of a laboratory ([email protected]).
Abstract: The method is suggested to solve one of the most important problems in the field of effective heat supply both in the case of the strong hydraulic deregulation of a heating system (a static problem) and in the case of varying thermal loadings (a dynamic
problem). The numerical procedure is developed for heating system analysis and optimization, which allows solving heat-consumption regime optimization and control problems for considerably large heat networks. Both static and dynamic problems are shown to reduce to a single general nonlinear optimization problem. The optimization algorithm is built where goal function minimization is carried out in an iterative process using the method of coordinate-wise descent.
Keywords: mathematical modeling, nonlinear optimization, minimization problems, numerical techniques, heat supply networks.
Статья представлена к публикации членом редакционной коллегии Н. Н. Бахтадзе
Уважаемые читатели! Обращаем Ваше внимание на выпущенный в 2010 году Институтом проблем управления им. В.А. Трапезникова РАН тематический сборник статей «Сетевые модели в управлении». Подробная информация о Сборнике доступна на сайте www.mtas.ru.