CC BY
100
УДК 519.17:6 Дата подачи статьи: 03.10.17
DOI: 10.15827/0236-235X.031.2.387-395 2018. Т. 31. № 2. С. 387-395
ПРИМЕНЕНИЕ АЛГОРИТМА ПЕРЕБОРА ДЕРЕВЬЕВ И МЕТОДА ИМИТАЦИИ ОТЖИГА ДЛЯ СХЕМНО-СТРУКТУРНОЙ ОПТИМИЗАЦИИ ТЕПЛОВЫХ СЕТЕЙ
В.А. Стенников 1, чл.-корр. РАН, ВРИО директора, sva@isem.irk..ru А.А. Чемезов 1, аспирант, a.a..chemezov@gmail.com
1 Институт систем энергетики им. Л.А. Мелентьева СО РАН, ул. Лермонтова, 130, г. Иркутск, 664033, Россия
В статье рассматривается математически сложная задача схемно-структурной оптимизации тепловых сетей, излагаются ее постановка, методические подходы и алгоритмы решения. В практике проектирования эта задача традиционно решается сопоставлением 2-3 заранее намеченных вариантов схемы. Значительное усложнение схемы и масштабов систем обусловило то, что получаемые решения далеки от оптимальных и нередко приводят к слабой загруженности участков сети и даже их неработоспособности. Развитие методов математического программирования и оптимизации открывает новые широкие возможности для решения практических задач. Вместе с тем это требует понимания энергетических, математических особенностей и соответствующих связей между ними с целью правильного применения математических методов. Накоплен большой арсенал методов, выбор из их числа эффективных требует проведения специальных исследований, связанных с возможностями и ограниченностью применения.
Сложность задачи схемно-структурной оптимизации заключается в том, что целевая функция в виде приведенных затрат является выпуклой по расходам и вогнутой по напорам. Фиксируя напоры путем преобразований, удается свести ее к вогнутой функции. Задача по своей сути является многоэкстремальной, что характерно для задач вогнутого программирования. Оптимальное решение будет иметь вид дерева. Каждому возможному варианту дерева, соответствующего вершине многогранника ограничений, будет отвечать локальный минимум целевой функции. Именно это делает данную задачу трудноформализуемой, что не позволяет найти аналитическое решение.
В статье рассматривается исторически применяемый метод перебора деревьев, предлагаются его модификации, а также другие методы, ранее не применявшиеся для оптимизации тепловых сетей. Осуществляется сравнительный анализ предлагаемых методов и алгоритмов, дается оценка эффективности их практического применения для расчета сетей различных масштабов и сложности.
Реализации алгоритмов имеют универсальный характер и могут применяться для различных типов энергетических систем, имеющих сетевую структуру.
Ключевые слова: методы теории графов, случайное дерево, алгоритм отжига, избыточная схема, оптимизация, системы теплоснабжения, тепловые сети.
Современные теплоснабжающие системы представлены сложными пространственно распределенными тепловыми сетями (ТС) многоуровневой иерархической структуры, связывающей источники тепловой энергии с потребителями. Задача минимизации затрат на строительство и эксплуатацию ТС продолжает оставаться актуальной из-за увеличивающихся масштабов городов и других инфраструктурных систем, для которых теплоснабжение является неотъемлемой частью нормального функционирования [1].
Для решения сложных комплексных прикладных задач, заключающихся в синтезе ТС, имеющих нелинейное математическое представление, используется методический подход, основанный на их декомпозиции на структурное и параметрическое решения. Наиболее сложной и трудноформа-лизуемой является задача схемно-структурной оптимизации ТС по экономическому критерию.
Разработкой подходов к решению задачи структурной оптимизации ТС занимались многие отечественные ученые [1-3]. В исследованиях последних лет для решения такой задачи активное применение находят различные эвристические алгоритмы [4-8].
Для решения задачи схемно-структурной оптимизации ТС успешно применяется развиваемая в ИСЭМ СО РАН (г. Иркутск) методика избыточных проектных схем [9]. Избыточная схема представляет собой схему системы теплоснабжения, в которую включены всевозможные способы и направления прокладки сети трубопроводов от источников тепловой энергии до потребителей.
В соответствии с этой методикой задача схемно-структурной оптимизации ТС в общем виде формализуется следующим образом: на заданной избыточной схеме ТС требуется выделить сеть, соединяющую узлы производства и потребления тепловой энергии таким образом, чтобы обеспечить минимум расчетных затрат на ее строительство и эксплуатацию, при соблюдении физико-технических ограничений. По своей сути она сводится к задаче оптимального потокораспределения в ТС.
Постановка задачи
Структура избыточной схемы задается в виде графа О* = (V, Е*), где V = ^ |] = 0, 1, ..., т} -множество узлов графа; Е* = {е, = {ум, Ук}| и ф к, , = 0, 1, ..., п} - множество ветвей, соответствую-
щих участкам схемы, т - число узлов, п - число участков. Множество V определяется как объединение непересекающихся множеств V = V0 и V1 и V2, где V0 - множество потребителей; V — множество источников; V2 - множество узлов разветвления. Избыточная схема дополняется следующими параметрами элементов схемы: Qj - расход теплоносителя в узле j; р , р - верхнее и нижнее ограничения на давление в узле у; ¡1 - длина участка ,;
, щ - верхнее и нижнее ограничения на скорость движения теплоносителя по участку ,.
Структура выделяемой сети теплоснабжения определяется подграфом О = (V, Е), где Е с Е*. Оптимальная топологическая структура сети теплоснабжения с учетом экономического критерия находится среди множества древовидных схем избыточной схемы [3], поэтому известные алгоритмы ее поиска относительно минимальной функции затрат Е(О) ориентированы на топологию сетей древовидной структуры. В общем виде алгоритм оптимальной сети включает несколько взаимосвязанных этапов, в числе которых нормализация структуры (определение и исключение ветвей с нулевым расходом), расчет потокораспределения сети, определение суммарных приведенных затрат. Содержание и последовательность реализации этого алгоритма могут быть представлены следующим образом.
1. Граф О = (V, Е), где Е с Е, ориентируется относительно источников теплоносителя е, = (уи, Ук), е, е Е, где Уи - начальная вершина; у/ - конечная вершина.
2. Для каждого узла Vj определяется М,- = = {т | ет = (Vу, Ук)} - множество ветвей, направленных от узла Vу.
3. По заданным нагрузкам, начиная с концевых узлов и исходя из материальных балансов в узлах сети, определяются расходы на участках:
х=а + Xхт' е= (vJ' V).
(1)
4. Производится нормализация сети, для каждой ветви с нулевым расходом V ei е В = { ei = (уу,
Ук) е Е| Xi = 0} положим hi = 0, Н = 0, Мк = 0, Wi = = щ .
5. Из-за ограничения на скорость течения теплоносителя принято потери давления на ветвях рассчитывать относительно рекомендуемой скорости.
Скорость потока wi жидкости на участке i рассчитывается по формуле 4 • х.
f
3,6 ■ ж- p ■
d.
I
1000
Y
e. e E / B
(2)
где V - скорость течения жидкости, м/с; х, - расход жидкости (объемная скорость), т/ч; р - плотность жидкости, кг/м3; й, - внутренний диаметр трубы, мм; л - число Пи.
Потери давления И, на участке i [4] рассчитываются по формуле
ю. ■ x2 ■I -(1 + а.)
h = -d,e. eE/B
i d i
(3)
где ф,- - коэффициент, зависящий от эквивалентной шероховатости труб; а, - коэффициент местных потерь.
Выразив диаметр из (2) и заменив его в выражении (3), получим зависимость потерь давления от величины расхода теплоносителя:
ф • х2 • I (1 + а)• (3,6• л• р• щ)2,625 ... к=^' ' '(-' ) ( , 2 62р ')-, е'Е / В .(4)
' / л \2,625 '
(4 • х •ЮОО2)
Потери давления И, на участке , при заданной (рекомендуемой) скорости течения теплоносителя V рассчитываются по формуле (4) как функции от переменной расхода теплоносителя на участке сети.
6. Рассчитываются значения давления в узлах и действующего напора на участках.
Для концевых узлов давление принимается равным среднему между верхним и нижним ограничениями давления в этом узле:
P + h Р =—-- ,М. = !
(5)
Для остальных узлов V,- давление Ру берется равным максимальному рассчитанному среди ветвей
из Му\
р = р + к - Н
] кг ' '
Vm е М. :(Р + к - Н * > Р, + к - Н *,
J ^ а т т к ' ' '
ет = (^^ , ^а ) , е = (^^ , ^к ) ^ т = ') . (6)
Производится проверка соблюдения ограничений на давление [ р < р < р ] в узлах, и, если требуется, осуществляется установка насосных (дросселирующих) станций (НС, ДС).
Правила выбора действующего напора НС на участке е, = (уу-, Ук) е ЕВ включают
P - P при p + h > P, H * =\p - P при P, + h < P ,
1 —j j r k 1 —j'
(7)
0 при р < р + к < р.
7. Осуществляется увязка давления в узлах, общих для нескольких участков. Производится корректировка давлений таким образом, чтобы для всех ветвей, выходящих из него, значение давления было равно максимальному из рассчитанных:
Иг = Рк + И, - Ру-, е, = (У-, Ук) е Е/В. (8)
8. Производится расчет общих суммарных затрат на основе массивов переменных х, И, Н для сети О:
п
Р (°) = X (рС (х,к ) + FN (х,Н ,) + ГЕ (х,,к )) ,
'=1
где (х, к) - расчетные затраты на сооружение и эксплуатацию трубопровода на участке с но-
w
мером К" (х , Н) - расчетные затраты на сооружение и эксплуатацию НС на ,-м участке; гЕ (х, И) - расчетные затраты на электроэнергию, расходуемую на перекачку теплоносителя по ,-му участку и подачу его потребителю.
Математическая постановка задачи с учетом вышеизложенного формулируется следующим образом: требуется минимизировать функцию приведенных затрат
п
£ (К? (х, И ) + Р,М (х, Н) + (X,, И)) ^ шт (9)
I=1
при соблюдении следующих ограничений:
Р < Р < Р , V. е V, (10)
—11 Р 1 ' 4 '
щ < ш < щ, е. е Е, в = (V,Е). (11)
Эта задача относится к задачам дискретного нелинейного программирования. Среди алгоритмов ее решения хорошо зарекомендовал себя алгоритм целенаправленного перебора деревьев (ПД) [6].
Целенаправленный ПД
Алгоритм ПД заключается в поэтапном улучшении дерева начального приближения путем замены входящих в него участков хордами (участками, не входящими в текущее решение). При включении хорды в дерево образуется цикл. При исключении любого участка, входящего в цикл, будет определено новое дерево [10].
В качестве дерева начального приближения (О0) берутся деревья минимальной суммарной длины, максимальной суммарной длины, минимальных расстояний от потребителей до источника, максимальных расстояний от потребителей до источника, случайно сгенерированные остовные деревья. На каждой г-й итерации алгоритма целенаправленного ПД на базе принятого дерева начального приближения генерируется новое дерево (О,).
Алгоритм целенаправленного ПД представляет последовательность следующих операций.
1. Определяется начальное приближение О0 = = (V, Ео).
2. Формируется массив хорд и = {и е Е*/Ео}.
3. Берется первая хорда и = (у-, Ук) из массива хорд и е и. Хорда и исключается из и. Определяется образовавшийся в результате включения хорды и контур С = РИ- Д РИк, где РИ- с Е - путь от у- до источника, РИк с Е - путь от ук до источника.
4. Генерируются новые решения О е А, А = = {О = (V, Е) | Е = {и} + Ег / {с}} путем поочередного исключения одного из участков выделенного контура с е С и включения хорды и.
5. Дерево О,+1 определяется как дерево с наименьшими расчетными затратами:
\в,, К (в )< К (о ш1"),
в+1 = (12) в ш,",к (о )> к (ошш),
V О е А: (Е(О)< Е(От1П) ^ О = От1П).
6. Если массив хорд пуст (и = 0), работа алгоритма завершается, иначе осуществляется переход к шагу 3.
Повышение эффективности алгоритма ПД возможно путем его модификации, заключающейся в динамическом изменении массива хорд, или на основе перехода к новому решению аналогично алгоритму имитации отжига [11].
Модификация динамического массива хорд
В качестве модификации алгоритма ПД предлагается участки, исключенные из структуры, и хорды, рассмотренные, но не включенные в решение, добавлять в конец массива хорд. Предположение об эффективности данного подхода основано на том, что исключенный или не включенный в схему участок в случае дальнейшего изменения структуры решения на следующих итерациях алгоритма может быть снова рассмотрен и включен, тем самым улучшая решение. Чтобы вычислительный процесс не зацикливался, вводится счетчик количества итераций без изменения текущего решения. Условием завершения алгоритма принимается равенство значения счетчика числу хорд в сети. Данное условие гарантирует обнаружение локального минимума, так как в окрестности текущего решения отсутствуют сети, для которых значение функции затрат меньше текущего решения.
Модифицированный алгоритм перебора деревьев.
1. Принимается, что пО соответствует количеству хорд, а п1 представляет счетчик итераций без изменения решения.
2. Определяется начальное приближение О0 = = (V, Ео).
3. Формируется массив хорд и = {и е Е*/Ео}.
4. Определяется пО = |Ц|, п1 = 0.
5. Берется первая хорда и = (у-, Ук) из массива хорд и е и. Хорда и исключается из и. Определяется образовавшийся в результате включения хорды и контур С = РИ- Д РИк, где РИ- с Е - путь от у- до источника, РИк с Е - путь от ук до источника.
6. Генерируются новые решения О е А , А = = {О = (V, Е) | Е = {и} + Ег / {с}} путем поочередного исключения одного из участков выделенного контура с е С и включения хорды и.
7. Дерево Ож определяется как дерево с наименьшими расчетными затратами (12).
8. Если Ог+1 ф Ог, то исключенная ветвь с добавляется в конец массива хорд и, обнуляется счетчик п1 = 0, иначе рассмотренная хорда и добавляется в конец массива и и увеличивается значение счетчика п1 + 1.
9. Если п1 = п0, то нужно завершить алгоритм, иначе перейти к шагу 5.
Изложенный алгоритм сохраняет рассмотренные хорды для последующих расчетов, что позволяет расширить область исследований и повысить точность получаемых результатов.
Метод отжига
Метод отжига основывается на имитации процесса кристаллизации вещества при понижении температуры [9]. Предполагается, что атомы уже выстроились в кристаллическую решетку, но еще допустимы переходы отдельных атомов из одной ячейки в другую. Процесс протекает при постепенно понижающейся температуре. Переход атома из одной ячейки в другую происходит с некоторой вероятностью, причем вероятность уменьшается с понижением температуры. Устойчивая кристаллическая решетка соответствует минимуму энергии атомов, поэтому атом либо переходит в состояние с меньшим уровнем энергии, либо остается на месте, что может быть записано в следующем виде: 1, если Г ( вп)- Г ( вь )< 0,
p ( Gn, Gb ) =
exp
F ( G")-F ( Gb)
•100
(13)
р{сь)-1
если Г ( вп)-Г ( вь )> 0.
Начальное приближение определяется как Оь, и производится выбор плана Оп из окрестности текущего плана Оъ. Вероятность перехода (13) от плана Оь к Оп регулируется коэффициентом 1°, имитирующим падение температуры. Чем выше показатель, тем больше вероятность перехода к плану, хуже текущего Е(Оп) > Е(Оъ). Такой подход позволяет выйти из локального минимума функции Е. Если р(Оп, Оъ) больше случайного числа из интервала (0, 1), то на следующей итерации рассматривается Оъ = Оп. С течением времени значение 1° убывает, что обеспечивает сходимость алгоритма к локальному оптимуму.
Необходимо экспериментально выбрать начальное ?1аг1 и конечное значения температуры и закон Т его изменения. В качестве Т предлагается использовать один из следующих подходов:
^аИ
- отжиг Коши Т (И) =-, где И - номер ите-
и
рации алгоритма;
- метод тушения Т(1°) = с • 1°, где с - некоторая константа, 0,7 < с < 1.
Алгоритм ПД и колец с имитацией отжига
Алгоритм представляет собой комбинацию модифицированного алгоритма ПД и алгоритма имитации отжига. Переход от одного решения к другому в алгоритме ПД осуществляется в соответствии с законом перехода алгоритма имитации от-
жига, который включает следующие основные этапы.
1. Принимается, что n0 соответствует количеству хорд, а n 1 представляет счетчик итераций без изменения решения.
2. Определяется начальное приближение G0 =
= (V, Eo).
3. Формируется массив хорд U = {u е E'/Eo}.
4. Определяется n0 = |U|, n1 = 0, f = ftart.
5. Берется первая хорда u = (Vj, vk) из массива хорд u е U. Хорда u исключается из U. Определяется образовавшийся в результате включения хорды u контур C = Phj Д Phk, где Phj с E - путь от Vj до источника, Phk с E - путь от vk до источника.
6. Генерируются новые решения G е A, A = = {G = (V, E) | E = {u} + Et/{c}} путем поочередного исключения одного из участков выделенного контура c е C и включения хорды u, Gt+1 = Gt.
7. V G е A, еслиp(G, Gt+\) > r при (0 < r < 1), где r - случайное число, то Gt+\ = G.
8. Если Gt+1 Ф Gt, то исключенная ветвь c добавляется в конец массива хорд U, обнуляется счетчик n1 = 0, иначе рассмотренная хорда u добавляется в конец массива U и увеличивается значение счетчика n1 + 1.
9. Если n1 > n0 и t° < tstop, то завершить алгоритм, иначе t° = T(t°, t), перейти к шагу 5.
Использование принципов имитации отжига определяет вероятностный характер процесса поиска решения. Однако ввиду того, что с увеличением числа выполненных итераций алгоритма вероятность переходов к менее оптимальным решениям стремится к нулю, процесс поиска решения приобретает детерминированный характер и вырождается в предложенную выше модификацию алгоритма ПД. Вероятностный характер на первых итерациях алгоритма позволяет снизить зависимость получаемого решения от начального приближения.
Сравнительный анализ методов схемно-структурной оптимизации ТС
С целью более компактного изложения результатов исследования введем сокращенные названия алгоритмов:
ПД - исходный алгоритм целенаправленного ПД;
ПДХ - модификация алгоритма ПД с динамическим добавлением хорд в массив хорд;
МОК - модификация ПДХ с законом перехода метода имитации отжига, в качестве закона изменения параметра температуры принят отжиг Коши;
МОТ - модификация ПДХ с законом перехода метода имитации отжига, в качестве закона изменения параметра температуры принят метод тушения;
СП - случайный поиск (генерация и расчет случайных деревьев).
В рамках вычислительного эксперимента сравнивались результаты работы алгоритмов для четырех случайно сгенерированных избыточных и различных по масштабам и параметрам ТС (табл. 1).
Число возможных вариантов решений определяется числом остовных деревьев, связывающих все вершины графа. Структура сети возможного решения является древовидной, но необязательно может быть представлена в виде единого дерева (это может быть лес деревьев). Вместе с тем, если узлы-источники с помощью фиктивных ветвей свести в один узел, то такой граф будет деревом. Для определения количества возможных решений избыточной схемы воспользуемся матричной теоремой о деревьях. Для этого найдем матрицу Кирхгофа K для графа с совмещенными узлами-источниками G*. Число вариантов остовных деревьев равно определителю матрицы, получающейся из матрицы K путем удаления 1 -й строки и 1 -го столбца.
Характеристики рассчитываемых сетей, включая количество возможных различных деревьев, сведены в таблице 1.
Таблица 1
Характеристики расчетных сетей
Table 1
Characteristics of calculated networks
Характеристики сети Количество воз-
Сеть Участ- Уз- Источ- Узлов раз- можных вариан-
ков лов ников ветвления тов деревьев
1 332 200 1 30 2,20159E+75
2 506 250 2 0 2,0796E+104
3 600 100 10 10 3,7534E+101
4 2016 1000 3 10 _^ да
Каждому расчету соответствует определенное дерево начального приближения. Последовательность расчетов на заранее определенных начальных приближениях представлена в таблице 2.
Таблица 2
Деревья начального приближения, принятые в расчетах, и последовательность их использования в предложенных алгоритмах
Table 2
Initial approximation trees in the calculations and a sequence of their use in the proposed algorithms
Начальное приближение Номер расчета
Дерево минимальной суммарной длины 1
Дерево максимальной суммарной длины 2
Дерево минимальных расстояний от потребителей до источника 3
Дерево максимальных расстояний от потребителей до источника 4
Случайное дерево 5, 6, ..., да
Для каждой сети выполнялась серия расчетов по каждому алгоритму. Одному проходу алгоритма соответствует набор входных параметров, включая заранее определенное дерево начального приближения. В результате одного расчета алгоритма определяются структура сети, суммарные затраты и время, потраченное на поиск решения. Серия расчетов алгоритма для сети 4 ограничивалась заданным временем вычислений t, для остальных сетей - количеством расчетов алгоритма r . Параметры алгоритмов для тестовых сетей приведены в таблице 3.
Таблица 3
Параметры и ограничения алгоритмов для расчетных сетей
Table 3
Parameters and constraints of the algorithms for calculated networks
Алгоритм Сеть 1 Сеть 2 Сеть 3 Сеть 4
ПД - - - -
ПДХ - - - -
МОК tstart = 200 ftop = 0,1 tstart = 100 tstop = 0,05 /start = 1 tstop = 0,001 ftart = 100 tstop = 0,01
МОТ tstart = 10 с = 0,99 ftop = 0 tstart = 10 с = 0,99 ftop = 0 tstart = 1,0 с = 0,99 ftoP = 0 ftart = 100 с = 0,99 ftoP = 0
СП Не выполнялся Не выполнялся Не выполнялся -
Ограничения r=1000 t=1800с
Для сопоставления результатов работы алгоритмов аналогично работе [9] определялись следующие показатели для каждой серии расчетов:
minF(G), midF(G), maxF(G) - соответственно минимальное, среднее арифметическое и максимальное значения F(G) среди всех результатов расчетов по всему набору алгоритмов для данной сети (меньше - лучше), тыс. руб.;
mid/ - среднее арифметическое значение затраченного времени на один расчет рассматриваемого алгоритма (меньше - лучше), сек.;
midG - среднее арифметическое количество рассмотренных сетей за один расчет рассматриваемого алгоритма;
midp - среднее арифметическое отклонение найденных решений от лучшего найденного решения среди всех алгоритмов (меньше - лучше), %;
X - отношение количества расчетов (best) по соответствующему алгоритму, в которых обнаружено лучшее решение F(G) = minF(G), к общему количеству расчетов (all) алгоритма (больше -
. „. . best
лучше), %: X =-;
all
tmin - расчетное время, необходимое для поиска наилучшей сети с минимальным значением F(G) =
= minF(G) (меньше - лучше), сек.: tmn = midt
rmin - расчетное количество расчетов алгоритма, необходимое для нахождения сети с F(G) =
= minF(G) (меньше - лучше): rmm = 1.
В таблице 4 отражены результаты работы рассмотренных выше алгоритмов. Экспериментальные расчеты выполнялась на компьютере на базе процессора Intel® Core™ i5-4670K.
Визуально результаты всех проведенных исследований согласно рассматриваемым схемам систем теплоснабжения представлены на рисунках 1-4 в виде точечных графиков. Ось абсцисс соответствует значениям функции затрат F(G) для полученных решений, ось ординат отражает требуемое на их поиск время. Каждая точка представляет полученное решение работы алгоритма, а именно значение функции затрат найденной сети и время, затраченное на поиск. Чем точка ближе к началу координат, тем меньше времени потрачено на вычисление и выгоднее решение. Выделенные области решений соответствуют распределению результатов, полученных определенным алгоритмом.
Представленная на рисунках 1-4 область решений, полученная с помощью ПД, имеет наибольший разброс результатов относительно других рассматриваемых алгоритмов. Это объясняется сильной зависимостью получаемого результата от начального приближения. Множества локальных минимумов для первой и второй сетей, найденных алгоритмами МОК, МОД и ПДХ, как видно из рисунков 1 и 2, принимают форму столбцов из точек. Их графическое сопоставление отражает наличие небольшого количества локальных минимумов,
при этом, чем столбец плотнее, тем больше начальных приближений приводят к решениям с равными значениями функции затрат в процессе оптимизации. Крайний левый столбец определяет решение, соответствующее условно-глобальному минимуму функции затрат. Высота столбца отражает распределение затраченного времени на поиск решения. Различие затрат времени на поиск решения объясняется отличием начальных приближений. Для каждого начального приближения алгоритмом для поиска оптимального решения выполняется различное количество итераций, на выполнение которых затрачивается различное время расчета. Для сетей 3 и 4 локальные минимумы представлены размытыми областями решений. Это свидетельствует о наличии большого числа минимумов, из которых исследуемые алгоритмы не способны выйти, используя предложенную для них логику вычислений. Вместе с тем распределение точек области решений для алгоритмов МОК, МОД и ПДХ смещено в сторону оптимума функции затрат. Наличие значительного количества обнаруженных локальных минимумов функции затрат свидетельствует о высокой сложности задачи оптимизации для данной сети в отличие от случаев, когда число обнаруживаемых минимумов незначительное, например, для сетей 1 и 2.
Анализ результатов экспериментальных исследований (табл. 4, рис. 1 -4) показывает, что предложенная модификация динамического массива хорд обеспечивает прирост эффективности работы алгоритма ПД в 2-5 раз в процентном отношении (midp), но при этом затрачивается больше времени
Сопоставление результатов работы алгоритмов Comparison of the algorithm output
Таблица 4 Table 4
Сеть Алгоритм minF(G) midF(G) maxF(G) mid t mid G mid p k tmin rmin
ПД 71126,0 72787,9 74018,0 0,26c 2617,4 2,63% 0,10% 260c 1000
1 ПДХ 70922,4 71746,1 72575,4 0,78c 7962,9 1,16% 9,70% 9c 11
МОК 70922,4 71191,3 72209,6 1,55c 15626,3 0,38% 29,70% 6c 4
МОТ 70922,4 71255,2 72575,4 1,12c 11072,1 0,47% 32,20% 4c 4
ПД 150662,5 153829,5 160643,6 0,42c 1281,1 2,13% 0,10% 416c 1000
2 ПДХ 150620,9 151576,9 152485,5 1,68c 5154,8 0,63% 7,60% 24c 14
МОК 150620,9 150785,8 152092,3 3,27c 9998,9 0,11% 74,60% 7c 2
МОТ 150620,9 151235,9 152834,0 2,58c 7668,8 0,41% 29,70% 10c 4
ПД 950636,6 968046,9 1017854,2 0,30c 1209,3 1,90% 0,10% 298c 1000
3 ПДХ 949990,3 955862,6 974523,4 1,00c 4153,4 0,62% 1,20% 84c 84
МОК 949990,3 955528,0 973488,5 1,09c 4525,9 0,58% 1,30% 84c 77
МОТ 949990,3 954288,0 966761,3 1,58c 6619,1 0,45% 1,80% 88c 56
ПД 365134,3 375475,6 388603,2 10,68c 6383,3 3,39% 0,60% 1784c 167
ПДХ 364019,4 365170,9 369091,8 66,75c 38610,2 0,55% 3,85% 1736c 26
4 МОК 363158,0 364429,6 366246,1 105,27c 66334,1 0,35% 6,25% 1684c 16
МОТ 363240,8 364794,2 367187,6 77,64c 44681,9 0,45% 4,55% 1708c 22
СП 436955,2 615916,3 3794069,7 0,002c 1,0 69,60% 6,3E-6% 1800с 159109
о
70500 71000 71500 72000 72500 73000 73500 74000 74500
Расчетные затраты, т.р.
Рис. 1. Результаты работы алгоритмов для сети 1 Fig. 1. The algorithm output for network 1
Рис. 2. Результаты работы алгоритмов для сети 2 Fig. 2. The algorithm output for network 2
Рис. 3. Результаты работы алгоритмов для сети 3 Fig. 3. The algorithm output for network 3
0
360000 365000 370000 375000 380000 385000 390000
Расчетные затраты, т.р.
Рис. 4. Результаты работы алгоритмов для сети 4 Fig. 4. The algorithm output for network 4
(в 2-3 раза) на один расчет (т1&). Однако ПДХ при достаточном количестве расчетов (гтт) находит условно-глобальный минимум (сети 1 и 2) или точки, близкие к нему (сеть 4), в то время как ис-
ходный алгоритм ПД не находит глобальные опти-мумы для рассчитываемых сетей.
Предложенные модификации имитации отжига имеют большую вероятность попадания в глобаль-
ный оптимум (X) по сравнению с ПДХ при достаточном количестве расчетов согласно предложенному алгоритму и подходящих параметрах.
Для больших сетей (сеть 4) при ограниченном времени расчета (? ) модификации имитации отжига находят более выгодные решения, чем ПД и ПДХ, несмотря на меньшее количество расчетов. В общем случае для МОК и МОТ время 1тт, обеспечивающее нахождение глобального минимума функции (9) или близкого к нему, меньше чем у ПДХ.
Сравнительный анализ решений, полученных с помощью алгоритмов имитации отжига, показывает, что лучшие результаты имеет подход, использующий отжиг Коши.
Проведенные вычислительные исследования показали, что предложенные модификации алгоритма ПД являются более эффективными для решения задачи схемно-структурной оптимизации ТС.
Следует отметить, что предложенные алгоритмы могут быть использованы для решения оптимизационных задач в других отраслях экономики, где решение может быть представлено в виде древовидной структуры, например, для задач транспортного планирования, газо-, нефте-, водоснабжения, оптимизации систем водоотведения и вентиляции, проектирования электрических сетей и др.
Заключение
Изложенный в статье методический и вычислительный инструментарий ориентирован на решение сложных задач схемно-структурной оптимизации теплоснабжающих систем на уровне предпро-ектных и проектных документов. Предложенные методические и алгоритмические разработки отличаются высокой вычислительной эффективностью и позволяют получать решения, близкие к глобальному оптимуму, и, что является не менее важным, эти решения работоспособны, реализуемы и применимы на практике. Предлагаемый инструментарий может быть представлен в виде меню с целью выбора наиболее эффективного из имеющихся методов получения решения при расчете систем различной размерности и сложности.
В рамках решения задачи схемно-структурной оптимизации следует отметить, что предлагаемые алгоритм ПД и его модификации предназначены для оптимизации теплоснабжающих сетей древовидной структуры. Однако в последнее время предъявляются повышенные требования к надежности и эффективности работы ТС. С точки зрения структуры ТС это приводит к необходимости структурного резервирования, что подразумевает наличие в сети кольцевых структур. Кроме того, кольцевание сети может быть экономически целесообразным при реконструкции существующей сети, когда новые и старые участки в определенной конфигурации образуют кольцо, которое выгоднее древовидной структуры. Вышеупомянутые проблемы структурной оптимизации ТС с кольцевыми участками задают направление дальнейших исследований.
Литература
1. Меренков А.П., Хасилев В.Я. Теория гидравлических цепей. М.: Наука, 1985. 175 с.
2. Некрасова О.А., Хасилев В.Я. Оптимальное дерево трубопроводной системы // Экономика и мат. методы. 1970. Т. 4. № 3. С. 427-432.
3. Чупин В.Р., Майзель И.В., Чупин Р.В., Житов А.В. Схемно-структурная оптимизация систем водоотведения поверхностного водостока и ливневой канализации // Изв. вузов. Инвестиции. Строительство. Недвижимость. 2014. N° 4. С. 97-105.
4. Keedwell E., & Khu S.T. A hybrid genetic algorithm for the design of water distribution networks. Eng. Appl. Artif. Intell., 2005, vol. 18, no. 4, pp. 461-472.
5. Kadu M.S., Gupta R., & Bhave P.R. Optimal design of water networks using a modified genetic algorithm with reduction in search space. Jour. of Water Resources Planning and Management, 2008, vol. 134, no. 2, pp. 147-160.
6. Zecchin A.C., Simpson A.R., Maier H.R., Marchi A., & Nixon J.B. Improved understanding of the searching behavior of ant colony optimization algorithms applied to the water distribution design problem. Water Resources Research, 2012, vol. 48, no. 9.
7. Johns M.B., Keedwell E., & Savic D.A. Pipe smoothing genetic algorithm for least cost water distribution network design. Proc. GECCO, 2013, pp. 1309-1316.
8. McClymont K., Keedwell E., & Savic D., Randall-Smith M. A general multiobjective hyper-heuristic for water distribution network design with discolouration risk. Jour. Hydroinform, 2013, vol. 15, no. 3, pp. 700-716.
9. Системные исследования в энергетике: ретроспектива научных направлений СЭИ-ИСЭМ; [отв. ред. Н.И. Воропай]. Н.: Наука, 2010. 686 с.
10. Наумов И.В., Ямщикова И.В. Математическое обоснование выбора оптимизационной модели трассировки электрической сети // Евразийский союз ученых. 2015. № 7-2. С. 123-127.
11. Лопатин А.С. Метод отжига // Стохастическая оптимизация в информатике. 2005. Т. 1. № 1-1. С. 133-149.
Software & Systems Received 03.10.17
DOI: 10.15827/0236-235X.031.2.387-395 2018, vol. 31, no. 2, pp. 387-395
APPLICATION OF A TREE SEARCH ALGORITHM AND A ANNEALING SIMULATED METHOD IN OPTIMIZATION OF HEAT NETWORK CONFIGURATION AND STRUCTURE
V.A. Stennikov1, Corresponding Member of RAS, Acting Director, sva@isem.irk.ru A.A. Chemezov1, Postgraduate Student, a.a.chemezov@gmail.com
1 Melentiev Energy Systems Institute SB RAS, Lermontov St. 130, Irkutsk, 664033, Russian Federation
Abstract. The paper considers a mathematically complex problem of heat network configuration and structure optimization. It presents a problem statement and methodological approaches and algorithms to solve it.
Traditionally in design practice, comparing 2 or 3 pre-planned alternate schemes solve this problem. Significant complication of the scheme and system scale led to the fact that the obtained solutions are far from optimal ones and often lead to insufficient load of network sections and even their inoperability. The development of mathematical programming and optimization methods gives new possibilities for solving practical problems. At the same time, it requires understanding of energy and mathematical features, as well as corresponding connections between them in order to apply mathematical methods correctly.
There is a large arsenal of methods. Choosing effective ones requires special research related to capabilities and limited application.
The complexity of scheme-structural optimization problem is in the fact that the objective function in the form of reduced costs is convex in terms of flow rate and concave on the heads. Fixing the heads by transformations, it can be reduced to a concave function. The problem is multiextremal as concave programming problems. An optimal solution will look like a tree. Each possible variant of the tree corresponding to the vertex of a constraint polyhedron will correspond to the local minimum of the objective function. This makes this task difficult to formalize, which does not allow finding an analytical solution.
The paper considers a historically applied method of tree search, suggests its modifications, as well as other methods that were not previously used to optimize heat networks. The authors compare the proposed methods and algorithms, and assesse their performance in calculation of networks of different sizes and complexity.
The algorithm realizations are universal and might be applied to different types of energy systems with a network structure.
Keywords: graph theory methods, random tree, simulated annealing algorithm, redundant diagram, optimization, heat supply systems, heat networks.
References
1. Merenkov A.P., Khasilev V.Ya. Teoriya gidravlicheskikh tsepey [Theory of Hydraulic Circuits]. Nauka Publ., 1985, 175 p.
2. Nekrasova O.A., Khasilev V.Ya. Optimal tree of a pipeline system. Ekonomika i mat. metody [Economics and Mathematical Methods]. 1970, vol. 4, no. 3, pp. 427-432 (in Russ.).
3. Chupin V.R., Maizel I.V., Chupin R.V., Zhitov A.V. Scheme and structural optimization of water discharge systems of surface watercource and a storm-sewage system. Izvestiya vuzov. Investitsii. Stroitelstvo. Nedvizhimost [News of Higher Educational Institutions. Investments. Building. Property]. 2014, no. 4 (9), pp. 97-105 (in Russ.).
4. Keedwell E., Khu S.T. A hybrid genetic algorithm for the design of water distribution networks. Eng. Appl. Artif. Intell. 2005, vol. 18, no. 4, pp. 461-472.
5. Kadu M.S., Gupta R., Bhave P.R. Optimal design of water networks using a modified genetic algorithm with reduction in search space. Jour. of Water Resources Planning and Management. 2008, vol. 134, no. 2, pp. 147-160.
6. Zecchin A.C., Simpson A.R., Maier H.R., Marchi A., Nixon J.B. Improved understanding of the searching behavior of ant colony optimization algorithms applied to the water distribution design problem. Water Resources Research. 2012, vol. 48, no. 9.
7. Johns M.B., Keedwell E., Savic D.A. Pipe smoothing genetic algorithm for least cost water distribution network design. Proc. GECCO, 2013, pp. 1309-1316.
8. McClymont K., Keedwell E., Savic D., Randall-Smith M. A general multiobjective hyper-heuristic for water distribution network design with discolouration risk. J. Hydroinform. 2013, vol. 15, no. 3, pp. 700-716.
9. Sistemnye issledovaniya v energetike: retrospektiva nauchnykh napravleny CEI-ISEM [System research in the energy sector: a retrospective of SEI-ISEM scientific directions]. N.I. Voropay [Ed.]. Novosibirsk, Nauka Publ., 2010, 686 p.
10. Naumov I.V., Yamshchikova I.V. Mathematical justification of choosing an optimization model of electric network tracing. Evrazysky soyuzuchenykh [Eurasian Union of Scientists]. 2015, no. 7-2 (16), pp. 123-127 (in Russ.).
11. Lopatin A.S. Annealing method. Stokhasticheskaya optimizatsiya v informatike [Stochastic Optimization in Computer Science]. 2005, vol. 1, no. 1-1, pp. 133-149 (in Russ.).
Примеры библиографического описания статьи
1. Стенников В.А., Чемезов А.А. Применение алгоритма перебора деревьев и метода имитации отжига для схемно-структурной оптимизации тепловых сетей // Программные продукты и системы. 2018. Т. 31. № 2. С. 387-395. DOI: 10.15827/0236-235X.031.2.387-395.
2. Stennikov V.A., Chemezov A.A. Application of a tree search algorithm and a annealing simulated method in optimization of heat network configuration and structure. Programmnye produkty i sistemy [Software & Systems]. 2018, vol. 31, no. 2, pp. 387-395 (in Russ.). DOI: 10.15827/0236-235X.031.2. 387-395.
