CC BY
35
ТЕОР1Я I МЕТОДИ АВТОМАТИЧНОГО УПРАВЛ1ННЯ
ТЕОРИЯ И МЕТОДЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ
THEORY AND METHODS _OF AUTOMATIC CONTROL_
УДК 681.5.01
Дорофеев Ю. И.1, Любчик Л. М.2
1Канд. техн. наук, доцент кафедры системного анализа и управления, Национальный технический университет
«Харьковский политехнический институт», Украина 2Д-р техн. наук, профессор, заведующий кафедрой компьютерной математики и математического моделирования, Национальный технический университет «Харьковский политехнический институт», Украина
РОБАСТНОЕ ПОДАВЛЕНИЕ ВОЗМУЩЕНИЙ ПРИ УПРАВЛЕНИИ НАСОСНЫМИ СТАНЦИЯМИ В СИСТЕМЕ ЦЕНТРАЛИЗОВАННОГО _ВОДОСНАБЖЕНИЯ_
Предложен подход к решению задачи синтеза стабилизирующего робастного управления запасами воды в системе подачи и распределения воды крупного города. Математическая модель системы подачи и распределения воды представлена в виде нелинейной дискретной модели в пространстве состояний с запаздыванием. Предложена методика факторизации матриц модели, описывающих влияние нелинейных термов, которая позволила представить заданные структурные ограничения в виде линейных матричных неравенств. Для подавления влияния возмущений, моделирующих изменения неизвестного, но ограниченного внешнего спроса, одновременно с обеспечением устойчивости замкнутой системы, применена методика инвариантных эллипсоидов, которая позволила сформулировать задачу в терминах линейных матричных неравенств, а синтез управления свести к последовательности задач одномерной выпуклой оптимизации и полуопределенного программирования. В качестве примера рассмотрен фрагмент системы подачи и распределения воды города Харьков.
Ключевые слова: система подачи и распределения воды, управление запасами воды, метод инвариантных эллипсоидов, прямой метод Ляпунова, линейное матричное неравенство, задача полуопределенного программирования.
НОМЕНКЛАТУРА
ЛМН - линейное матричное неравенство; СПРВ - система подачи и распределения воды; А1 - матрица динамики; В{ - матрицы влияния управлений; С - матрица выходов;
О - допустимое множество значений внешних возмущений;
О, - диаметр трубы, исходящей из резервуара г; Е - матрица влияния внешних возмущений; Я, - радиус резервуара г;
и - допустимое множество значений управляющих воздействий;
X - допустимое множество значений состояний;
© Дорофеев Ю. И., Любчик Л. М., 2014
d (к) - вектор внешних возмущений в момент времени к;
Стах - вектор верхних граничных значений объемов потребляемой воды и утечек;
Ст1П - вектор нижних граничных значений объемов потребляемой воды и утечек;
ё - ускорение свободного падения; кг (к) - напор (давление) воды в узле г в момент времени к ;
к - номер дискретного момента времени; т - количество управляемых гидравлических элементов (насосов и вентилей);
п - количество резервуаров;
q - количество секторов потребления;
Ягп,г (к) - г -й входящий поток воды в момент времени к;
164 DOI 10.15588/1607-3274-2014-2-21
Чош, / (к) - }-й исходящий поток воды в момент времени к;
и (к) - вектор управляющих воздействий в момент времени к;
и шах - вектор, определяющий максимальные мощности насосов;
V, (к) - объем воды в резервуаре , в момент времени к;
хшш - вектор, определяющий минимально допустимые объемы воды в резервуарах;
Хшах - вектор, определяющий объемы резервуаров; - высота расположения центра тяжести узла ,;
А/ - период дискретизации;
Л - целочисленная переменная, кратная А/, которая определяет максимальное значение интервала времени, необходимого для транспортировки воды между резервуарами сети;
аК - коэффициент Кориолиса.
ВВЕДЕНИЕ
Система водоснабжения представляет собой жизненно необходимую составляющую существования любого города. Управление системами водоснабжения крупных городов по мере их роста становится предметом пристального внимания специалистов в области автоматического управления. Проблема устойчивого развития и управления системой водоснабжения для города Харьков стоит особенно остро. Предпосылками являются территориальная расположенность города в относительно маловодном районе, большое водопотребление населением и промышленностью, групповое водоснабжение целого ряда прилежащих к мегаполису городов, высокая изношенность инженерных коммуникаций, возросшие требования к ка -честву водоснабжения со стороны населения.
Современные системы водоснабжения, представляющие собой комплекс сооружений для подъема, очистки и подачи потребителям воды, с точки зрения их функционирования можно представить в виде трех относительно независимых по характеру и критериям управления подсистем [1]: источники водоснабжения и водоприемные сооружения, очистные сооружения, СПРВ.
Наиболее дорогостоящей и сложной среди указанных подсистем является СПРВ, которая имеет специфические характеристики, усложняющие ее анализ и управление:
- СПРВ имеет распределенную крупномасштабную архитектуру;
- резервуары и исполнительные элементы имеют эксплуатационные ограничения;
- для подачи воды в высотные здания необходимо обеспечить определенное давление воды, значение которого нелинейно зависит от объема перекачиваемой воды, что приводит к необходимости использования нелинейных моделей;
- транспортирование воды требует определенных временных затрат, что приводит к необходимости учета запаздываний.
Гидравлические элементы, из которых состоят СПРВ, делятся на две категории: управляемые и неуправляемые. Управляемыми являются насосные станции и вентили. Неуправляемыми элементами являются трубы, резервуары, а также узлы, которые соответствуют точкам сети, где водные потоки объединяются или разделяются. Еще одним элементом моделей СПРВ являются секторы потребления, которые позволяют моделировать объемы спроса на воду со стороны потребителей.
Для обеспечения непрерывного и эффективного функционирования СПРВ необходимо создание запасов питьевой воды. Управление запасами воды заключается в определении объемов воды, которую необходимо закачать в резервуары, и моментов времени, когда необходимо это делать. В системах управления запасами воды предусматривается периодический контроль состояния через равные промежутки времени, а момент принятия решения о пополнении запаса обычно совпадает с началом каждого периода.
В настоящее время в большинстве случаев управление СПРВ осуществляется с помощью эвристических правил и «исторически сложившихся» стратегий, которые являются результатом многолетнего эмпирического опыта. Безусловно, эти стратегии вцелом являются адекватными, однако они не гарантируют оптимальной политики управления в современных крупномасштабных сетях. В связи с этим весьма актуальной является задача разработки оптимальных стратегий управления системами водоснабжения, в основе которых лежат математические модели работы сети и методы современной теории автоматического управления.
Целью настоящей работы является синтез робастной стабилизирующей стратегии управления запасами воды в системе централизованного водоснабжения, гарантирующей полное и своевременное удовлетворение неизвестных, но ограниченных объемов потребления воды при условии минимизации затрат на производство и транспортировку воды с учетом заданных эксплуатационных ограничений.
1 ПОСТАНОВКА ЗАДАЧИ
В работе [2] представлены различные подходы к построению математических моделей СПРВ. В настоящей работе в качестве переменных состояний модели рассматриваются объемы воды в резервуарах, измеряемые в м3, и величины давления (напора) воды в исходящих из резервуаров потоках, измеряемые в м. Предполагается, что структура СПРВ известна, а состояния доступны непосредственному измерению. В качестве управляющих переменных выступают объемы воды, перекачиваемые через управляемые гидравлические элементы в текущем периоде, измеряемые в м3. Внешними возмущениями являются объемы потребляемой воды и объемы утечек за период времени А/, также измеряемые в м3.
Для любого узла сети можно записать уравнение баланса массы в виде разностного уравнения, учитывающего все входящие и исходящие потоки:
V,(к +1) = V, (к) + £ д1П,,(к) - £ дои(,} (кX I = 1П.
Связь между объемом перекачиваемой воды и напором воды в потоке определяется уравнением Бернулли, которое выражает закон сохранения механической энергии в жидкости и может быть записано для любых двух произвольных узлов сети г и ] по ходу потока:
Ьк + (к)2 + Р^ = к!(к) + (к)2 + , (!)
п О п и ^
где и (к) е Ят; а (к) е Я
Напор воды в потоке, исходящем из резервуара, имеющего форму цилиндра, зависит от объема воды и радиуса резервуара:
h (k) = -Ц- V, (k).
nR,2
(2)
Чтобы избавиться в уравнениях (1) от постоянных величин г введем переменные hj(к) = hj(к) + pg(- ) и запишем уравнения с учетом (2) в следующем виде:
hj (k) = -^rVi (k) + Fu2 (k) + Hfd2 (k), (3) nR,2
где F, 6 R 8a k p
nXffl т-р _ 8aKp
Fl = • diag(1.0,..,0),
n Di
Fm =_TKT• diag(0,...,0,l),H} 6 Rnxq
^n 4
n Dm
Hi = -
8a K P
n2 Di4
• diag(l,0,...,0), ..., Hq =^-80KPr• diag(0,...,0,i).
^n 4 n Dq
Представим вектор и 2(k) = {и,- 2(k)}
u 2(k) = \uf (k)f i = 1, m , элементами которого являются значения управляющих воздействий в квадрате, в следующем виде:
u2(k) = eUuT(k) • El1 • u(k) +... + eumuT(k) • Eum • u(k),
где eu e Rmx1 : e:u = (l,0,...,0)T,..., em = (0,...,0,l)T ,
Eu e Rmxm : Eu = diag(l,0,...,0),..., Em = diag(0,...,0,l). Аналогичным образом представим вектор
d2(k) = {dj2(k)} j = 1,q , элементами которого являются значения внешних возмущений в квадрате:
d2(k) = ef dT (k) • E:d • d(k) +... + edqdT (k) • Ed • d(k),
где
ef 6 Rqxi : ef = (l,0,...,0)T,..., ed = (0,...,0,l)T
Еа е к. = ахав(1,0,...,0),..., Е^ = ^(,...,0,1).
Тогда нелинейные слагаемые в уравнениях (3) могут быть представлены в виде кусочно-непрерывных нелинейных функций: первая зависит от дискретного времени и вектора управлений Ьи (к,и): Ят+1 ^ Яп, вторая -
от дискретного времени и вектора возмущении hd (k, d): R q+i ^ R n:
m m
hu (k, и) = £ F, • и 2(k) = £ F •Su (k),
i=i
i=i
hd (k, d) = £ Hj • d2(k) = £ Hj • Sf (k),
j=i
j=i
где Еи (к) = £ е> т(к) Е?и(к), (к) = £ т(к )Е*а (к).
г=1 j=1
Для моделирования запаздываний, связанных с транспортировкой воды от источников до резервуаров, а затем до потребителей, используется модель дискретной задержки, поскольку предполагается, что значения интервалов времени, требуемых для транспортировки воды по каждому участку сети, известны и не меняются в процессе функционирования.
Тогда математическая модель СПРВ может быть представлена в виде нелинейной дискретной модели в пространстве состояний с запаздыванием следующего вида: л
х(к +1) = А1х(к) + £ Вги(к - г) + Ей (к) + ЬХ (к, и) + ЬХ (к, а), (4)
г=0
где x(k)=
— Т
v(k) , h(k) J 6 R n - составной вектор со-
стоянии модели;
> 2nxm
4 6 R
«m+l
2nx2n
E 6 R
2nxq.
Bt 6 R2nxm, t = 0,Л ; h„x : Rm+i ^ R2n, : Rq+i ^ Rn,
0nxm 0nxq
N m £ Fr •Su (k ^ hl = q £ Hj •S d (k).
_i=1 _ _j=i _
Очевидно, что структура сети определяется матрицами Ai, Bt, E, методика построения которых изложена в работе [3].
В процессе функционирования системы должны выполняться эксплуатационные ограничения:
x(k) 6 X = { x 6 Rn : xmin < x < xmax } ,
u(k) 6U = {u 6 Rm : 0 < и < иmax }, (5)
min max „ max
где векторы x , x и и считаются заданными.
Предполагается, что векторы внешних возмущений удовлетворяют ограничениям:
d(k) 6 D = {d 6 Rq : dmin < d < dmax}
где векторы dmin и dmax предполагаются известн^1ми.
Преобразуем модель (4) к стандартному виду без запаздываний с помощью расширения вектора состояний [4] путем включения в него векторов предыдущих управляющих воздействий, определяющих объемы воды,
находящейся в процессе транспортировки. Вектор состояний расширенной модели будет иметь вид:
^(к) = [(к)т, И(к)т, и (к - 1)т, и(к - 2)т,..., и(к - Л)т].
Тогда уравнения расширенной модели сети примут
вид:
£(к +1) = А^(к) + Ви (к) + Оа (к) + иех (к, и) + (к, а), (6) х(к) = С^(к),
где Неих : Ят+1 ^ Я *, ке/ : Яя+1 ^ Я*, N = 2п + Лт;
ии^ (к, и) = Р£и (к), Иех (к, е) = ЯХа (к); матрицы модели имеют соответствующую блочную структуру:
А =
А1 В • ■ ВЛ-1 Вл В0
тх2п 0тхт 0тхт 0тхт 1тхт
тх2п 1тхт 0тхт 0тхт , В = 0тхт
тх2п 0тхт 1тхт 0тхт _ _°тхт
О =
£
утхд ^тхд
-^тхд
сТ
12пх2п 0пхт 0пхя
0тх2п , р = т I Р ¡=1 , я = я IЯ7 7=1
_ 0тх2п _ _0Лтхт _ 0Лтхя
Предполагается, что пара матриц (А, В) является стабилизируемой, а нелинейные функции являются ограниченными и удовлетворяют квадратичным неравенствам [5]:
Иех (к, и)Кх (к, и) < а2 ^ (к)РТ^ (к) <
< а2ииТ(к)
/ \Т
т
I е> Т(к) Е?
V г=1
Р Т Р
т А
I е> Т(к) Е V г=1
и (к) <
< о^и Т(к)
(т \
I еии шах Е?
V г=1
Т
р Т Р
< а ии (к )2и Р Р^и и (к),
, 2 >Т %т.шах^ ^огшах
^ т Л
I егии шах Е? V ¡=1
и(к) <
Т
ке/ (к, а)ке/ (к, а) < а ее Т(к)£шах я1 я£ 7х а (к),
я
где = I е.ии шахЕ.и, ЕШаХ = I еаа шахЕа,
¡=1 7=1
аи > 0, а а > 0 - некоторые скаляры.
Для системы (6) рассматривается задача синтеза ро-бастной относительно неопределенных, но ограниченных внешних воздействий а (к) е Б, стратегии управле-
ния запасами, обеспечивающей для любого начального состояния х(0) е X минимизацию критерия качества, который гарантирует безопасность работы системы, а также асимптотическую робастную устойчивость замкнутой системы при условии выполнения заданных ограничений на состояния и управления (5).
2 ОБЗОР ЛИТЕРАТУРЫ
Выбор модели управления запасами определяется характером спроса со стороны внешних потребителей. В настоящее время для синтеза стратегии управления запасами с заданной моделью спроса широко применяется метод прогнозирующего управления [6], основанный на решении оптимизационной задачи, в которой последовательность управлений вычисляется в каждый момент времени так, чтобы минимизировать некоторый критерий качества, определяемый прогнозируемыми состояниями системы на некотором временном горизонте. После определения последовательности управляющих воздействий только первый ее элемент используется для управления в текущий момент времени. Затем измеряется новое состояние системы и процедура повторяется, т.е. реализуется принцип отступающего горизонта. Указанный подход позволяет строить управление как в программном виде, так и в виде обратной связи, а также учитывать ограничения на состояния и управления.
Однако, на практике, как правило, отсутствует информация для построения адекватной модели внешнего спроса, которая необходима для синтеза прогнозирующего управления. Одним из подходов к решению задачи управления запасами в условиях неопределенности спроса является использование концепции «неизвестных, но ограниченных» воздействий [7]. При этом соответствующая модель спроса характеризуется интервальной неопределенностью.
В последнее десятилетие сформировался новый подход к рассматриваемой проблематике, основанный на концепции инвариантных множеств. Инвариантные множества широко используются в различных задачах теории автоматического управления в динамических системах при наличии неопределенностей [8]. Среди различных форм инвариантных множеств особо выделяются эллипсоиды вследствие их простой структуры и прямой связи с квадратичными функциями Ляпунова.
В рамках метода инвариантных эллипсоидов [9] в качестве технического средства используется математический аппарат ЛМН [10]. После того, как были развиты вычислительные методы, основанные на идеях выпуклой оптимизации, и для их реализации были разработаны соответствующие алгоритмы и программное обеспечение, ЛМН стали рассматриваться в качестве общего метода анализа и синтеза линейных систем как в непрерывном, так и в дискретном случае. Однако в большинстве работ, посвященных задаче подавления ограниченных внешних возмущений, техника ЛМН применяется для подавления возмущений, ограниченных в какой-либо норме. Тогда как спецификой задач управления запасами воды является неотрицательность значений
Т
переменных и наличие эксплуатационных ограничений, что приводит к необходимости учета несимметричных ограничений на значения внешнего спроса, состояний и управляющих воздействий.
3 МАТЕРИАЛЫ И МЕТОДЫ
Согласно теореме об аппроксимации произвольных выпуклых множеств, не обладающих свойством симметрии относительно начала координат [11], множество О может быть аппроксимировано эллипсоидом:
Е(/, Ра) = е Я1: (((к) - а*)ТРа_1(((к) - /)< ),
матрица которого Ра и вектор а *, определяющий координаты центра, определяются на основании известных граничных значений объемов потребляемой воды:
( 2 2 А
Ра = ахав ^(агах - аг1)2,..,^^ - С")2
а * = 2 (а т1п + а тах).
Аналогично выполним аппроксимацию множества X эллипсоидом максимального объема
Е(х\ Рх) = !х е Я2п : (х(к) - х*)Рх_1(<к) - х*)< Л (7)
у которого вектор х , определяющий координаты центра, определим позднее, а матрица эллипсоида Рх вычисляется на основании векторов хт1п и хтах:
Рх = аЦ±-хт1п, хтах -х*} 1 („,■ < * т1п тах * и2 |
...-^ШШх - хп , хп - хпЦ I .
Построим закон управления в виде нестационарной линейной обратной связи по сигналу рассогласования между наличным и страховым уровнями запаса воды:
и(к) = К (к )(^(к) - )
(8)
где К (к) е ЯтхМ - нестационарная матрица коэффициентов обратной связи в момент к .
Значения элементов вектора определяющего уровень страховых запасов воды, вычисляются на основании средних значений объемов потребления с помощью продуктивной модели Леонтьева:
Г =
*т *т
V ,01хп, V ,
= (/ - п )-1 а
-1 лтеап
= 1 2
1 (а гх
. лтт
+ а ]
), j = 1,1,
0, j = 1 +1, п,
(9)
где П - матрица смежности ориентированного графа, описывающего модель сети.
Значения элементов вектора х*, определяющего координаты центра эллипсоида (7) будут равны -|Т
*
х =
*т
V ,0
1хп
Расширенную модель замкнутой системы для управления (8) представим в виде:
^(к+1) = Аг (к )(^(к) - )+А^* + о(а (к) - а *) +оа + ьи^ (к, и)+(к, а),
х(к) = С^(к), Аг (к) = А + ВК (к).
(10)
Оптимальное управление запасами воды в СПРВ сводится к задаче построения стратегии управления потоками воды от источников к потребителям с целью полного и своевременного удовлетворения спроса на воду при условии оптимизации критериев качества, обеспечивающих безопасность и устойчивость работы системы.
Для безопасного функционирования системы водоснабжения и минимизации утечек, а также поддержания высокого качества воды ее объемы в резервуарах и величины напора не должны превышать максимально возможных и не опускаться ниже минимально допустимых, то есть быть приближенными к значениям, определяющим размер страховых запасов. Для математического выражения данного критерия в виде квадратичной функции в случае бесконечного временного горизонта используется следующее выражение:
(к) =£|(к) - )Т% ((к) - )
к=0
х N
где Я^ е Я *' - положительно определенная диагональная весовая матрица, элементы которой определяют штрафы за отклонение текущих значений уровня воды от страховых.
Экономические затраты обусловлены затратами на химическую обработку воды и на электроэнергию, необходимую для перекачивания воды с помощью насосов. Критерий, учитывающий стоимость производства и транспортировки воды, может быть описан следующим выражением:
да
• 2°° (к) = £ и т(к )Яии(к), к=0
где Яи е Ятхт - положительно определенная диагональная весовая матрица.
*
V
т
т
л
Регулирование мощности насосов должно осуществляться плавно, чтобы избежать нежелательных переходных процессов в герметичных трубах, что может привести к их повреждению. Для обеспечения эффекта сглаживания скачков управляющих воздействий предлагается критерий следующего вида:
J 3° (k) = £Au T(k )RAAu(k), k=1
Ди(к) = и(к) - и(к -1),
где Яд е Ятхт - положительно определенная диагональная весовая матрица.
В результате критерий качества синтезируемой системы управления будет иметь вид:
j f (k) = jf (k)+jf (k)+j3° (k).
(11)
Стабилизирующие алгоритмы управления, как правило, основаны на оценивании верхнего граничного значения критерия качества с помощью функции Ляпунова. Определим квадратичную функцию Ляпунова, построенную на решениях системы (10):
V ((k) - )= ((k) - )T P(k )((k) - ), (12)
где P(k) = PT (k) e RN xN - симметрическая положительно определенная матрица.
В соответствии с прямым методом Ляпунова вычислим первую разность функции (12) по k в силу системы (10) и потребуем выполнения следующего неравенства:
V($(k +1) - )-V($(k) - )<-J"(k). (13) Если неравенство (13) выполняется Vk > 0, то функция V ($(k) - ) определяет верхнее граничное значение критерия (11) [12]:
max
d(k)eE(d*, Pd)
J f (k) < V
fe(k) - ).
(14)
Управляющие воздействия и (к) будем искать из условия минимизации верхней оценки критерия (11), что позволяет упростить решение задачи. Тогда в соответствии с (14) найдем управляющие воздействия и (к) из условия минимизации функции Ляпунова:
u(k) = arg min V\t(k) - )
u (k )eU
В соответствии с методом инвариантных эллипсоидов [8] их можно рассматривать как характеристику влияния внешних возмущений на траектории динамической системы. Для оценки степени влияния возмущений d (k) на выходы x (k) удобно использовать инвариантные эллипсоиды вида:
E *, P (k)) = e R N : | - $ * )T P ~\k)(( - ) < 1}, (15)
которые аппроксимируют множества достижимости
замкнутой системы (10) при действии возмущений
*
d(k) e E(d ,Pd). В соответствии с [12] введем матричную переменную:
Q(k) = Y(k) P "1(k),
(16)
где у (к) - положительный скаляр.
Тогда задача синтеза робастного стабилизирующего управления заключается в вычислении в каждый дискретный момент времени к матрицы коэффициентов обратной связи К (к), которая стабилизирует замкнутую систему (10) и обеспечивает минимизацию по критерию следа эллипсоида (15). Соответствующий результат представлен следующей теоремой.
Т е о р е м а. Рассмотрим систему (6), замкнутую с помощью закона управления (8) при ограничениях (5), и
пусть нестационарные матрицы Р(к) = у(к )0> _1(к) и
К (к) = 7 (к)0 "1(к) (17)
получены в результате решения оптимизационной задачи:
Т(к) (18)
min
Q(k ),Y (k ),au ,ad ,T1,T2,T3
при ограничениях на матричные переменные 0(к) = 0 Т(к) е Я *х *, 7 (к) е Ятх* и скалярные параметры у (к), аи, а а, т1, т 2, Т3:
у (к) > 0, т1 > 0, т 2 > 0, т 3 > 0, (19)
Q(k) ^ 0,
, 1 , ( - )T ,(=(k) - ) Q(k)
^ 0,
(20)
" Ок) 7Т(к)ЯдК 0 0 0 0 0 ©Т(к) ©Т(к) 0 ОкЯ 7Т(Щ/2 7Т(к)^шахГ
К%7(к) -у(к)К%К 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 ((4-I)T ((4-I)T 0 0 0 0
0 0 0 -Т(к)т1/ 0 0 0 I I 0 0 0 0
0 0 0 0 -1(кК1 0 0 I I 0 0 0 0
0 0 0 0 0 0 ^ша^я^Т 0 0 0 0 0 0
0 0 0 0 0 яг V) I та2, 0 0 0 0 0 0
©(к) 0 А-1 I I 0 0 I ОТ 0 0 0 0
©(к) 0 А-1 I I 0 0 О Ок) да/2 0 0 0
0 0 0 0 0 0 0 0 ■»(кт,/ 0 0 0
ЯШ 0 0 0 0 0 0 0 0 0 #У 0 0
Я/Чк) 0 0 0 0 0 0 0 0 0 0 #У 0
0 0 0 0 0 0 0 0 0 0 0
(21)
где ©(к) = А(к)0(к) + В7 (к), К = К (к -1), ЯЕ = Яи + ЯД,
^ 0,
у(к)Рх у(к)С у(к )С Т 0(к)
7(к)(^(к) - ^ ) У 0 :
7 Т(к)1 (ш
((к) -((к) - Г )шах )+7(к)
7 Т (к )1 (и
((к) -
(22)
(23)
0, (24)
0(к)
где «+» - псевдообращение Мура-Пенроуза.
Если задача (18)-(24) имеет решение, то система (6), замкнутая с помощью закона управления (8), для любого начального состояния х(0) е X и неопределенного, но ограниченного, внешнего возмущения
а (к) е Е (а *, Ра) является асимптотически робастно устойчивой при ограничениях (5).
Доказательство теоремы аналогично приведенному в работе [13].
Таким образом, задача синтеза робастного стабилизирующего управления запасами воды в СПРВ, описываемой нелинейной дискретной моделью (4) в пространстве состояний, при наличии транспортных запаздываний в условиях неизвестного, но ограниченного внешнего спроса сводится к задаче минимизации линейной функции (18) при ограничениях в виде ЛМН (19)-(24), которая при фиксированных значениях скалярных параметров аи, а а, ть т 2, т3 является задачей полуопределенного программирования.
Если оптимизационная задача (18)-(24) имеет решение, то найденное в соответствии с (8), (17) управляющее воздействие применяется для управления СПРВ в текущий момент времени. Затем измеряется новое значение вектора состояний и, в соответствии с принципом отступающего горизонта, задачи полуопределенного программирования и одномерной выпуклой оптимизации
решаются в реальном времени для вычисления нового значения матрицы коэффициентов обратной связи и, соответственно, нового управляющего воздействия.
4 ЭКСПЕРИМЕНТЫ
В качестве примера рассмотрим фрагмент СПРВ города Харьков, изображенный на рис. 1. Он состоит из наземного источника, трех насосных станций М4, и 5), двух подкачивающих насосов (, М3), пяти резервуаров
(у15..., У5) и четырех секторов потребления (а!,..., а4).
В табл. 1, 2 и 3 приведены значения параметров рассматриваемого фрагмента СПРВ.
Технологическая матрица П совпадает с матрицей смежности графа О = ({1,2,3,4,5}, {(2,1), (4,3), (5,2), (5,4)}), описывающего модель сети.
Выберем период дискретизации равным 1 час. Пусть заданы значения интервалов, требуемых для перекачивания воды Л1 = Л 2 = Л3 = Л 4 = 1, Л 5 = 3, а также граничные значения объемов потребления и утечек воды
а™п =[90,60,70,60]Т, ашах=[150,100,110,120]. Моделирование осуществлялось в течение 24 периодов, что соответствует суткам. С целью учета суточных колеба -
Площадка
[ТУВХ Роганская резервуаров "Донец" насЛм Л ос ев о'
и-, V
КТ 1 (ул. Столетова)
О".
КТ 2 (ул. Халтурина)
Рисунок 1 - Фрагмент системы подачи и распределения воды города Харьков
Таблица 1 - Характеристики резервуаров
Название Переменная Минимальный объем воды (м3) Максимальный объем воды (м3) Относительная высота г (м)
КТ 1 50 500 156,40
КТ 2 50 500 120,98
КТ 3 50 500 144,91
КТ 4 50 500 146,97
«Лосево» 300 3000 180,08
Таблица 2 - Граничные значения напора воды в системе
Название Переменная Минимальный напор воды (м) Максимальный напор воды (м)
КТ 1 К 3,5 5,5
КТ 2 3,0 5,0
КТ 3 3,5 5,5
КТ 4 ь4 3,0 7,0
Таблица 3 - Характеристики насосов
Название Переменная Минимальная мощность (м3/час) Максимальная мощность (м3/час) Относительная высота г (м)
Насос КТ 1 и1 0 150 156,40
Нас/ст 27 А и2 0 200 160,50
Насос КТ 3 и3 0 150 144,91
Нас/ст 27 Б и4 0 200 162,38
Роганская нас/ст и5 0 400 178,08
6 ОБСУЖДЕНИЕ
Значения весовых матриц Я^, Яи, Яд были выбраны из условия разрешимости оптимизационной задачи (18)-(24) и определяют размер эллипсоида, аппроксимирующего мно-жество достижимости замкнутой системы (10)
*
при действии возмущений а (к) е Е(а , Ра).
Анализ полученных результатов позволяет сделать вывод об эффективности предложенного подхода к решению задачи синтеза робастной стабилизирующей стратегии управления запасами воды в СПРВ, поскольку значения объемов воды в резервуарах, а также напора воды в трубах в условиях стохастически изменяющихся объемов потребления воды находятся в пределах заданных ограничений, а объемы перекачиваемой воды не превышают максимальных возможностей насосов.
ВЫВОДЫ
В работе предложен подход к решению задачи построения робастного стабилизирующего управления запасами воды в централизованной системе подачи и распределения воды. Основными особенностями рассмотренной задачи являются: нелинейность математической модели; неопределенность, но ограниченность внешнего спроса на воду; наличие запаздываний, связанных с транспортировкой воды; наличие эксплуатационных ограничений на значения состояний и управлений.
ний значения объемов потребляемой воды моделировались как нормально распределенные случайные величины с изменяющимся средним значением:
а1 ^ N(т1, 7), т1 е [60,120], а2 ^ N(т2,4),
т2 е[30,70], а3 ^ N(m3,6), т3 е [50,100],
а4 ^ N(т4, 3), т4 е [20, 60].
Значения страховых уровней запаса воды, вычисленные в соответствии с (9), равны
V* =[364,588,276,446,1034]Т Выберем начальные условия v(0) = [400,500,400,500,2000]т и значения элемен-
-7
тов диагональных весовых матриц ^ = 1,5 10 , ги = 1,0-10-5, Гд = 1,5-10-6. 5 РЕЗУЛЬТАТЫ
Моделирование выполнено в среде МЛТЬЛБ. Численное решение задачи (18)-(24) получено с помощью свободно распространяемого пакета СУХ [14]. Графики переходных процессов для узла 1 при аи = а а = 0,15,
-3
Т1 = т 2 = Т3 = 5,0 • 10 и стохастическом внешнем спросе приведены на рис. 2. На рис. 3 представлена фазовая траектория в пространстве состояний, определенном для резервуаров 1, 2 и 3.
Рисунок 2 - Графики переходных процессов в системе подачи и распределения воды: значения объема воды в резервуаре 1; б - значения напора воды в потоке, исходящем из резервуара 1; в - значения объема потребленной воды в секторе 1 и объема воды, перекачанной насосом 1
дослщження та шформацшш технологи. - 2013. - № 1. -
C. 16-27.
4. Blanchini F. Feedback control of production-distribution systems with unknown demand and delays / F. Blanchini, R. Pesenti, F. Rinaldi and W. Ukovich // IEEE Transaction on Robotics and Automation, Special Issue on Automation of Manufacturing Systems. - 2000. - Vol. RA-16. - No. 3. -P. 313-317.
5. Stipanovic D. M. Robust stability and stabilization of discrete-time non-linear systems: the LMI approach /
D. M. Stipanovic, D. D. Siljak // International Journal of Control. - 2001. - Vol. 74. - P. 873-979.
6. Bemporad A. Robust model predictive control: a survey /
A. Bemporad, M. Morari // Lecture Notes in Control and Information Sciences. - 1999. - Vol. 245. - P. 207-226.
7. Bertsekas D. P. Recursive state estimation for a set-membership description of uncertainty / D. P. Bertsekas, I. Rhodes // IEEE Trans. Automat. Control. - 1971. - Vol. 16. -P. 117-128.
8. Blanchini F. Set theoretic methods in control / F. Blanchini, S. Miani. - Boston : Birkhauser, 2007. - 504 p.
9. Хлебников М. В. Оптимизация линейных систем при ограниченных внешних возмущениях (техника инвариантных эллипсоидов) / М. В. Хлебников, Б. Т. Поляк,
B. М. Кунцевич // Автоматика и телемеханика. - 2011. -№ 11. - С. 9-59.
10. Linear matrix inequalities in system and control theory / [S. Boyd, El. Ghaoui, E. Feron, V. Balakrishnan]. -Philadelphia: SIAM, 1994. - 187 p.
11. Черноусьиэ Ф. Л. Оценивание фазового состояния динамических систем. Метод эллипсоидов / Ф. Л. Черноусько. -М. : Наука, Главная редакция физ.-мат. литературы, 1988. -320 с.
12. Kothare M. V. Robust constrained model predictive control using linear matrix inequalities / M. V. Kothare, V. Balakrishnan, M. Morari // Automatica. - 1996. - Vol. 32 № 10. - P. 1361-1379.
13. Lyubchyk L. M. Robust model predictive control of constrained supply networks via invariant ellipsoids technique / L. M. Lyubchyk, Y. I. Dorofieiev, A. A. Nikulchenko // Proc. IFAC Conf. on Manufacturing Modelling, Management and Control MIM'2013. - 2013. -P. 1618-1623.
14. Grant M. CVX: Matlab software for disciplined convex programming, version 2.0 beta [Electronic resource] / M. Grant, S. Boyd // Acces mode: http://cvxr.com/cvx.
Огатья поступила в редакцию 09.09.2014.
После дороботки 23.09.2014.
Рисунок 3 - Фазовая траектория и инвариантный эллипсоид, полученный на последнем шаге, для резервуаров 1, 2 и 3
Для подавления влияния возмущений, моделирующих изменения внешнего спроса, одновременно с обеспечением устойчивости замкнутой системы, применена методика инвариантных эллипсоидов, которая позволила сформулировать задачу в терминах линейных матричных неравенств, а синтез управления свести к последовательности задач одномерной выпуклой оптимизации и полуопределенного программирования.
Полученное управление зависит от выбранного значения уровня страховых запасов воды. В рамках предложенного подхода возможен выбор оптимальных значений страховых запасов, поскольку полученное решение задачи синтеза робастного управления задает, фактически, алгоритмическую зависимость между уровнем страховых запасов и оптимальным значением критерия качества.
СПИСОК ЛИТЕРАТУРЫ
1. Петросов В. А. Управление региональными системами водоснабжения / В. А. Петросов. - Харьков : Основа, 1999. - 320 с.
2. Brdys M. A. Operational control of water systems: structures, algorithms, and applications / M. A. Brdys, B. Ulanicki. - London: Prentice Hall, 1994. - 364 p.
3. Дорофеев Ю. И. Построение математических моделей управляемых сетей поставок с учетом запаздываний потоков / Ю. И. Дорофеев, А. А. Никульченко // Системш
в
а
Дорофеев Ю. I.1, Любчик Л. М.2
'Канд. техн. наук, доцент, Нацюнальний техшчний ушверситет «Харкiвський полiтехнiчний iнститут», Укра'на
2Д-р техн. наук, професор, Нацiональний техшчний ушверситет «Харювський полiтехнiчний iнститут», Укра'на
РОБАСТНЕ ПРИДУШЕННЯ ЗБУРЕНЬ ПРИ УПРАВЛШШ НАСОСНИМИ СТАЦ1ЯМИ В СИСТЕМ1 ЦЕНТРАЛ1ЗО-ВАНОГО ВОДОПОСТАЧАННЯ
Запропоновано пiдхiд до виршення задачi побудови стабiлiзуючого робастного управлшня запасами води в системi подачi i розподiлу води великого мюта. Математична модель системи подачi i розподiлу води подана у виглядi нелшшно'' дискретно'' моделi в просторi станiв i3 запiзненням. Запропоновано методику факторизаци матриць моделi, що описують вплив нелiнiйних термiв, яка дозволила подати задаш структурнi обмеження у виглядi лiнiйних матричних нерiвностей. Для придушення впливу збурень, що моделюють змiни невiдомого, але обмеженого зовшшнього попиту, одночасно iз забезпеченням стшкосп замкнуто'' системи, застосована методика iнварiантних елшсо'дов, яка дозволила сформулювати задачу в термiнах лiнiйних матричних нерiвностей, а синтез управлiння звести до послщовност задач одновимiрноi опуклою оптимiзацii та напiввизначеного програ-мування. Як приклад розглянуто фрагмент системи водопостачання мюта Харюв.
Ключовi слова: система подачi i розподiлу води, управлшня запасами води, метод iнварiантних елшсо'дов, прямий метод Ляпунова, лшшна матрична нерiвнiсть, задача напiввизначеного програмування.
Dorofieiev Yu. I.1, Lyubchyk L. M.2
'PhD, Associate Professor, National Technical University «Kharkiv Polytechnic Institute», Ukraine
2Dr. Sc., Professor, National Technical University «Kharkiv Polytechnic Institute», Ukraine
ROBUST DISTURBANCES REJECTION IN PUMPING STATIONS CONTROL FOR CENTRALIZED DRINKING WATER DISTRIBUTION SYSTEM
An approach to solving the problem of stabilizing robust water inventory control synthesis for drinking water distribution system of a large city is proposed. The mathematical model of drinking water distribution system is presented in the form of a nonlinear discrete state-space model with time-delay. The technique of model matrix factorization describing the influence of nonlinear terms, which allowed to introduce structural constraints in the form of linear matrix inequalities, is proposed. To suppress the disturbances influence simulating unknown but bounded external demand, while ensuring robust stability of the closed-loop system, is used the invariant ellipsoids technique, which allowed to formulate the control problem in terms of linear matrix inequalities. As a result the control synthesis problem is reduced to a sequence of one-dimensional convex optimization problems and semi-definite programming. As an example, a fragment of the Kharkiv drinking water distribution system is consider.
Keywords: drinking water distribution system, water inventory control, invariant ellipsoids method, Lyapunov's direct method, linear matrix inequality, semi-definite programming.
REFERENCES
1. Petrosov V. A. Upravlenie regional'nimi sistemami vodosnabzhenija. Khar'kov, Osnova, 1999, 320 p.
2. Brdys M. A., Ulanicki B. Operational control of water systems: structures, algorithms, and applications. London, Prentice Hall, 1994, 364 p.
3. Dorofieiev Yu. I., Nikulchenko A. A. Postroenie matematicheskih modelej upravliaemih setej postavok s uchetom zapazdivanij potokov, System Research & Information Technologies, 2013, No. 1, pp. 16-27.
4. Blanchini F., Pesenti R., Rinaldi F., Ukovich W. Feedback control of production-distribution systems with unknown demand and delays, IEEE Transaction on Robotics and Automation, Special Issue on Automation of Manufacturing Systems, 2000, Vol. RA-16, No. 3, pp. 313-317.
5. Stipanovic D. M., Siljak D. D. Robust stability and stabilization of discrete-time non-linear systems: the LMI approach, International Journal of Control, 2001, Vol. 74, pp. 873-979.
6. Bemporad A., Morari M. Robust model predictive control: a survey, Lecture Notes in Control and Information Sciences, 1999, Vol. 245, pp. 207-226.
7. Bertsekas D. P., Rhodes I. Recursive state estimation for a set-membership description of uncertainty, IEEE Trans. Automat. Control, 1971, Vol. 16, pp. 117-128.
8. Blanchini F., Miani S. Set theoretic methods in control. Boston, Birkhauser, 2007, 504 p.
9. Hlebnikov M. V., Poljak B. T., Kuncevich V. M. Optimizacija linejnih sistem pri ogranichennih vneshnih vozmuschenijah (tehnika invariantnih ellipsoidov), Automation and Remote Control, 2011, № 11, pp. 9-59.
10. Boyd S., Ghaoui El., Feron E., Balakrishnan V. Linear matrix inequalities in system and control theory, Philadelphia, SIAM, 1994, 187 p.
11. Chernous'ko F. L. Ocenivanie fazovogo sostojanija dinamicheskih sistem. Metod ellipsoidov. Moscow, Nauka, 1988, 320 p.
12. Kothare M. V., Balakrishnan V., Morari M. Robust constrained model predictive control using linear matrix inequalities, Automatica, 1996, Vol. 32(10), pp. 1361-1379.
13. Lyubchyk L. M., Dorofieiev Y. I., Nikulchenko A. A. Robust model predictive control of constrained supply networks via invariant ellipsoids technique, Proc. IFAC Conf. on Manufacturing Modelling, Management and Control MIM'2013, 2013, pp. 1618-1623.
14. Grant M., Boyd S. CVX: Matlab software for disciplined convex programming, version 2.0 beta. Acces mode: http:// cvxr.com/cvx.
