CC BY
26
Гидравлика и инженерная гидрология
УДК 502/504:628.14
С. Н. Карамбиров, доктор техн. наук Контактная информация: тел. 8 (495) 153-97-66
С. А Трикозюк, аспирант
Контактная информация: 8 (495) 976-49-39
Федеральное государственное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства»
МНОГОРЕЖИМНАЯ СТОХАСТИЧЕСКАЯ ОПТИМИЗАЦИЯ СИСТЕМ ПОДАЧИ И РАСПРЕДЕЛЕНИЯ ВОДЫ
Определены новые подходы к задаче выбора оптимальных параметров (диаметров участков и напоров водопитателей) систем подачи и распределения воды. Расчет сделали одновременно для нескольких режимов с учетом стохастической неопределенности. При этом целевая функция является аддитивной по расчетным периодам (кроме пожара), а задача оптимизации предполагает процесс управления и должна удовлетворять всем ограничениям.
New approaches are determined to the task of optimal parameters choice (diameters of parts and pressures of water feeders) of the water supply and distribution system. The calculation was made simultaneously for several regimes taking into account stochastic uncertainty. The target function is additive on calculated periods (except fire) and the task of optimization presumes a process of control and must meet all restrictions.
Для нахождения экономически оптимальных диаметров линий водопроводных систем при заданном пото-кораспределении существует значительное число методов и их модификаций [1]. Практически все они основаны на детерминированных методах расчета для одного из характерных режимов
работы. В большинстве случаев выбор диаметров производится на час наибольшего транзита или на час наибольшего водопотребления [2]. Чаще используется второй случай — на час наибольшего водопотребления.
Остальные условия работы учитываются единственным коэффициентом
неравномерности расходования электроэнергии, обоснованность которого без проведения дополнительных гидравлических расчетов является весьма сомнительной. Кроме того, для схемы с контррезервуаром трудно учесть затраты на транзит воды в башню.
Системы подачи и распределения воды (ПРВ) по своей сути являются ресурсоемкими предприятиями системы водопользования, рациональное использование которых существенно улучшает экономические показатели работы всего водопровода. Удельный вес капитальных и эксплуатационных затрат, приходящихся на них, составляет до 60...80 % общей стоимости системы водоснабжения.
С появлением современных быстродействующих вычислительных машин существовавшие методы были автоматизированы, однако это не привело к реализации новых возможностей и качественному изменению постановки задач. Во многом это связанно с разде-
лением труда между специалистами-водоснабженцами, ставящими задачи, и профессиональными программистами, их реализующими. Даже в условиях одного расчетного случая величина во-допотребления является интервальной, что зафиксировано соответствующими нормами проектирования. Рассматривая истинное водопотребление как случайную величину, а отмеченные интервалы — доверительными, можно, с некоторыми допущениями, составить стохастическую модель многорежимной стохастической оптимизации.
Целью работы является создание новых методов оптимизации системы подачи и распределения воды, более адекватно описывающих реальные режимы работы систем водоснабжения.
Постановка задачи заключается в определении оптимальных диаметров участков и напоров водопитателей (для насосных станций зависящих от времени), доставляющих экстремальное значение целевой функции с учетом
Рис. 1. Расчетная схема сети
заданных ограничений, определяемых нормативными документами. Рассмотрим сеть, представленную на рис. 1.
Так как система является стохастической, случайной будет и целевая функция. В качестве оптимизирующей величины примем математическое ожидание дисконтированных затрат на строительство и эксплуатацию системы:
тп Г тг п \
I
t в
К,
+ ■
(1 + ЕУ (1 + Ef
(1)
где К — капитальные затраты на ¿-м интервале времени, величина которого принята равной одному году; С — текущие затраты на ¿-м интервале, за исключением капиталовложений; Е — норма дисконта (требуемая инвестором норма доходности на единицу авансированного капитала); Т — продолжительность расчетного периода (для систем водоснабжения 20-25 лет).
Будем считать, что капитальные затраты приведены к началу строительства, т. е. - = 0 и продолжительность строительства равна одному году. Тогда
Z = К0 + С1
1 + Е 1 + Е 1 -
г)
Т„+1
+ 1,
1 + Е
(2)
где К0 — капитальные затраты; С1 — ежегодные эксплуатационные затраты (средние за расчетный период).
Для решения задачи используем метод линейного программирования. Он базируется на выборе искомых переменных — протяженности труб стандартного сечения, укладываемых на каждом участке. Таким образом, участок разбивается на несколько последовательных телескопических секций (по-дучастков) из числа возможных к прокладке, среди которых впоследствии выбираются оптимальные.
Такое рассмотрение основано на том, что длины подучастков (при заданном расходе) являются единственными параметрами, от которых линейно зависят потери напора. Обозначим: / — номера сечений, укладка которых возможна на системе (/ = 1...та); Ji — под-многочленство / для ¿-го участка; х.. — протяженность труб /-го диаметра на ¿-м участке, которые являются искомыми переменными. Длина участка — это сумма длин его подучастков:
(i = 1-n), (3)
je Ji
где n — число всех участков.
Это задает первую систему ограничений. Для дальнейшего рассмотрения необходимо конкретизировать задачу.
Известно, что расходование воды из сети для расчетных суток определяется суммарными ступенчатыми графиками водопотребления 24-часовой размерности. Такие графики имеют несколько многочасовых ступеней (будем называть их обобщенными ступенчатыми графиками), вокруг которых происходят малые отклонения расходов. Обобщенные графики имеют меньшее число ступеней, определяемое пользователем по его внешнему виду. Обычно оно не превышает четырех. Мерой близости исходного и обобщенного графиков можно выбрать сумму квадратов их разностей для каждого часа суток. Обобщенный график формируется минимизацией указанной суммы методом динамического программирования [3]: k ti + 1
ЕЕ (Q06i - Q)2 min> (4)
i = lj = t i
Vi Q.
где ®o6i = ¿j 1 Т'' k — число ступеней обоб-j - ti h +1 _ h щенного графика; t e 1...24, ч.
Задачу можно записать так: fk(S) = min [Q(S, U) — fk_i(S')], (5)
где S — состояние системы (время переключения, ч); U — управление (ширина интервала обобщенного графика, ч); Q(S, U) — значение целевой функции (сумма квадратов разностей) при U-м варианте управления на k-м шаге; S' — состояние системы за k — 1 ступеней до конца графика; fk(S) — сумма квадратов разностей при оптимальном управлении за k ступеней до конца графика для S-го варианта состояния.
Расчет начинается с последней ступени. Полученный обобщенный график будем рассматривать как математические ожидания случайного водо-потребления, а усредненные квадраты разностей — как соответствующие дисперсии (рис. 2). Подобная информация может быть получена непосредственно по результатам измерений на
Рис. 2. Графики водопотребления: | | — исходный; | | — обобщенный
действующих объектах.
От стохастических характеристик водопотребления перейдем к соответствующим характеристикам потокорас-пределения (вероятностным распределениям расходов участков напоров в узлах и подач водопитателей).
Основой расчета является система уравнений, описывающих первый и второй законы Киргофа, а также напорно-расходные характеристики водопитате-лей, станций подкачки и регулирующих задвижек. В общем виде эта система может быть записана так:
F(X) = Q,
(6)
где Х = Н-"Нт' ?НС1-?ИСк)ТГ,1; ?1 — рас-
ходы участков; Н. — напоры узлов; дНС. —
подачи водопитателей; Q — нормально распределенный случайный вектор с математическим ожиданием Q и ковариационной матрицей
к„ =
к,
(7)
Используя линеаризацию исходной системы' получим' что вектор решения распределен нормально с мате-
матическим ожиданием
х = р~\Я)
и ковариационной матрицей
к = ^тк^]-1*'
где звездочкой показана операция транспонирования' а
(8) (9)
F' =
К
ôx1
К
дх1
дА дхг
дх„
(10)
х=х
Полученные выражения можно использовать для оценки вероятностей событий. Так, отказом узла Ь будем называть событие' при котором свободный напор в нем Н. меньше минимально необходимого Н .
Вероятность отказа:
Рот=Р(Н1<Яс1) = 0,б-Ф
i— > н, - н:
1 f -J-
где Ф(м) = _ е 2 dt — функция Лапласа.
л/2я Jo
№ 5' 2008
Вернемся к нашей задаче оптимизации. Пусть на ¿-м участке допускается п. типов труб с диаметрами ^, / = 1.п..
Математическое ожидание потерь напора на /-м подучастке для режима ¿.:
1Щ]
(12)
ЕЕ
<е;
1ЁЦ=1
(13)
^ (14)
я:: - я;
РОТ»
ш
^ы - ЕЕ
* + р}
ы
РОТ| аш'
(16)
что приводит к следующему ограничению:
Ы ^ доп 1
(18)
где к.., р.., т.. — зависящие от материала труб коэффициенты.
Тогда гидравлические ограничения второго закона Киргофа для кольца М и условия внешней увязки напоров ¿-го и /-го водопитателей для режима -можно записать в следующем виде:
где Ныи — математическое ожидание напора у Ы-го водопитателя для режима ^ — отметка Ы-го водопитателя; I — путь (участки), соединяющий водопитатели; е — точность увязки.
Условие (11) на требуемую вероятность безотказной работы можно записать так:
где и — процентная точка нормального распределения, соответствующая вероятности 1 - Р .
При оптимизации, как правило, не учитываются некоторые требования, предъявляемые к системе, например, намеченный график водоподачи в случае нескольких водопитателей. Даже если в начальном приближении требуемые подачи были выдержаны, технико-экономический расчет не накладывает на них ограничений, в результате чего подачи перераспределяются, что может привести к удорожанию сооружений. Чтобы учесть условия для подачи, зафиксируем диаметры смежных с водо-питателями участков и обозначим их параметры с помощью индекса 1Ь. Чтобы математическое ожидание подачи Ы водопитателя находилось для режима t¿ в интервале ?2Й], необходимо задать такие условия:
(19)
(15) \(К + гыу(нЦ + г1Ъ)\<
Чь
,т1Ь Ы1Ь
V (20)
где ирот — процентная точка нормального распределения, соответствующая вероятности РОт; тти и — соответственно математическое ожидание и среднеквадратичное отклонение напора ¿-го узла для режима
При малых значениях РОт величина ирот меньше нуля. Обозначим ее как (- | и I), тогда (15) запишется в виде
где правые части задают трехчленную формулу потерь напора.
Расшифруем значения К0 и С1 в выражении (2):
*о = ЕЕсА (21)
где С.. — стоимость единицы длины подучаст-ка /-го диаметра;
где Иыи — математическое ожидание напора водопитателя, питающего ¿-й узел в режиме (для разных режимов могут быть разными), а внешнее суммирование производится по любой последовательности участков, соединяющих ¿-й узел с его водопитателем.
У водопитателей вероятность превышения напора Н й над предельно допустимыми значением Н й не долж-
дот
на превышать заданного значения Р :
(17)
1=1 ; = 1
Кв
+2:1
О** .ст —г ,.
пер1 е ттЪ.т/-Ъ1 МЫКЫ>
(22)
и и 102 с Ь1
где Р. — амортизационные отчисления для ¿-го участка; ст — стоимость 1 кВт-ч электроэнергии; - . — длительность периода в году;
СпеР 1 _ _ «.,1
ь. — коэффициент полезного действия насосных установок водопитателей Ы; Кв — число насосных станций с математическим ожиданием подач фьи, л/с, и напоров Ньй, м, в период
¿ : -К"ы - 1 + ¿^и; К и — ковариация между ■"ьМЙ 4
подачей и напором ¿-го водопитателя в период (стоимость водопитателей для рассматриваемых вариантов принимается постоянной).
Минимизация функционала с ограничениями (3), (13), (14), (16), (18), (19), (20) составляют задачу линейного программирования. При этом общее число ограничений может составлять несколько тысяч.
Решение находим итерационно при заданных вероятностных характеристиках потокораспределения для каждого расчетного момента, а на случай пожара — с соответствующими корректировками требуемых значений напоров до совпадения начальных и конечных диаметров.
Вначале находим первое приближение параметров системы, которые обеспечивают ее работоспособность (возможно, не оптимально). Выбираем возможные диаметры стандартного сортамента в окрестности начальных. Затем определяем расчетные режимы системы и последовательно вызываем процедуры стохастической увязки и формирования
целевой функции системы ограничений (аддитивно для каждого периода) при заданном потокораспределении. В расчете учитываем процесс управления, т. е. определяем оптимальные напоры водопитателей для каждого периода. После этого с помощью симплекс-метода решаем задачу линейного программирования. Для разветвленных сетей в виде дерева на этом процесс оптимизации заканчивается. Для кольцевых сетей расходы участков зависят от их диаметров. Поэтому если они поменялись, изменится и потокораспределение, которое по условиям задачи считается заданным. Это вынуждает провести следующую итерацию и далее до совпадения диаметров на двух последних.
Результаты расчета приведены на рис. 3. Для сравнения на рис. 4 показана однорежимная оптимизация. Здесь каждый вертикальный скачок целевой функции означает новый вероятностно-гидравлический расчет, при котором система снова становится неоптимальной.
/" Ход оптимизации - П
1
Число итераций
Рис. 3. Ход процесса многорежимной стохастической оптимизации
№ 5' 2008
68
Рис. 4. Ход процесса однорежимной оптимизации
Отсюда видно, что применение многорежимной стохастической оптимизации позволяет снизить затраты на систему водоснабжения почти на 15 % .
Выводы
В статье приведено решение задачи многорежимной стохастической оптимизации системы по экономическому фактору с приближением подач водопитателей в характерные периоды к заданным значениям. Показано, что решение не требует округления и корректировки, оно является оптимальным для характерных режимов в целом и позволяет получить значительный экономический эффект.
Ключевые слова: водопроводная сеть, графики водопотребления, целевая функция оптимизации затрат.
Список литературы
1. Николадзе, Г. Н. Водоснабжение [Текст] : учеб. для вузов / Г. Н. Николадзе, М. А. Сомов. — М. : Стройиздат, 1995. — 688 с.
2. Расчет водопроводных сетей [Текст] : учеб. пособие для вузов / Н. Н. Абрамов [и др.]. — М. : Стройиздат, 1983. — 278 с.
3. Карамбиров, С. Н. Математическое моделирование систем подачи и распределения воды в условиях многорежим-ности и неопределенности [Текст] : монография / С. Н. Карамбиров. — М. : ФГОУ ВПО МГУП, 2004. — 196 с.
