НЕЧЁТКОЕ  УПРАВЛЕНИЕ НАСОСНОЙ СТАНЦИЕЙ

 

Кудинов Ю.И., Уварова Л.В., Кудинов И.Ю.

(Липецкий государственный технический университет)

В настоящей работе рассматривается методология построения нечёткого алгоритма управления насосной станцией (НС), поддерживающей уровень жидкости в городском резервуаре для промышленных, природных и бытовых стоков.

НС состоит из трех насосов  производительностью м3/ч, 700 м3/ч,  1000 м3/ч, обеспечивающих путём сброса стабилизацию уровня воды в резервуаре (рис. 1). Изменение уровня возможно за счёт притока природных осадков, канализационных и промышленных стоков.

Рис.1. Схема стабилизации уровня в резервуаре

 

Основными характеристиками объекта управления являются уровень жидкости  и скорость его изменения. Уровень жидкости не может быть выше зоны перелива = 8 м и ниже зоны опорожнения = 0.2 м.

На основании многолетних наблюдений были установлены нижний = -4 м (отрицательный) и верхний  = 4 м (положительный) пределы скорости изменения уровня в резервуаре. НС функционирует в условиях неопределенности, обусловленных неконтролируемыми значениями притока и состава жидкости, меняющего производительность насосов, а также зависимостью их производительности от уровня жидкости.  С достаточной степенью достоверности параметры процесса стабилизации (уровень и его скорость) можно охарактеризовать соответствующими интервалами их изменения и степенями принадлежности к этим интервалам, т.е. использовать нечёткие множества как для оценки параметров, так и для создания нечеткого алгоритма управления НС.

На рис. 2 а показано разбиение высоты резервуара у = 8 м, а на рис. 2 б - интервала [] изменения скорости уровня  на нечеткие множества  с лингвистическими значениями: NB отрицательное большое; NM отрицательное среднее; NSотрицательное малое; ZEблизкое нулю; PSположительное малое; PMположительное среднее; PB – положительное большое.

 

 

 

 

 


                                                         а)

 

 

 

 

 

 

Рис. 2. Функции принадлежности

Подпись:                                          Таблица 1
          Характеристики насосов
l	Pl	 	u1	u2	u3
1	400	0,99	-1	0	0
2	500	1,24	1	0	0
3	600	1,49	0	-1	0
4	700	1,74	0	1	0
5	850	2,11	0	0	-1
6	1000	2,48	0	0	1
7	1100	2,73	-1	1	0
8	1200	2,97	1	1	0
9	1350	3,32	-1	0	1
10	1500	3,69	1	0	1
11	1600	3,94	0	-1	1
12	1700	4,22	0	1	1
13	1950	4,80	-1	1	1
14	2200	5,46	1	1	1

Алгоритм управления насосами можно представить совокупностью нечетких правил

 если  есть  и  есть, то ,                (7)

где  нечеткое множество, характеризующее уровень жидкости в резервуаре  в момент времени , ; нечеткое множество, характеризующее скорость изменения уровня жидкости в резервуаре  в момент времени , ; u(t) вектор управлений. Цифрой (-1) обозначим импульсный режим работы, цифрой 1 включение, а цифрой 0 – выключение насоса или совокупности насосов, вектор управления которой  состоит из одного ; двух  , или трех  элементов. В столбце Pl вычисляется производительность включенных насосов и работающих в импульсной режиме из соотношений

Pl = 0.7, l =1;

, l = 2, 4, 6, 8, 10, 14;

Pl = Pl-1 + 0.5(Pl+1 - Pl-1), l = 3, 5, 7, 9, 11, 13.

Элементы столбца  рассчитываются по формуле (3), записанной в виде

.

Задание пределов изменения скоростей  дает возможность сформировать некоторую логическую функцию f q(·), которая определяет l-ую строку с управлениями u1,l, u2,l, u3,l, позволяющими стабилизировать уровень в резервуаре.

В правой части продукционных правил (7) присутствует четыре типа функции f θ(·) = (f1(·), …, f4(·)). Первый тип функции определяется выражением

,

которое означает отключение всех насосов, когда уровень жидкости достигает нижней отметки и характеризуется лингвистическим значением NS.

Второй тип функции  определяется выражением

,

согласно которому управления u1(t), u2(t), u3(t) в текущий момент времени t остаются неизменными и равными соответствующим управлениям u1(t-D), u2(t-D), u3(t-D) в предыдущий момент времени t-D.

Функция третьего типа  находит управления

в строке l (табл. 1), удовлетворяющей условию

.

Четвёртый тип функции  определяется выражением

,

которое соответствует включению всех насосов, если уровень достигает верхней отметки и характеризуется лингвистическим значением PS. Тогда правила (6) можно записать в унифицированной форме для каждого лингвистического значения уровня жидкости у(t) и скорости изменения ее уровня .

Интервал NS. 1. Если скорость уровня  отрицательная NB, NM, NS, т.е. он понижается, то отключаются все насосы, что соответствует правилам

R1: если у(t) есть NS,  есть NB, то u(t) = f1(·);

R2: если у(t) есть NS,  есть NМ, то u(t) = f1(·);

R3: если у(t) есть NS,  есть NS, то u(t) = f1(·).

2. Если скорость уровня  нулевая (ZE) или положительная, т.е. он повышается, то состояние насосов не меняется

R4: если у(t) есть NS,  есть ZE, то u(t) = f2(·);

R5: если у(t) есть NS,  есть PS, то u(t) = f2(·);

…………………………………………………………………..

Интервал PS. 1. Если скорость изменения уровня  отрицательная и равная NB, NM, NS, ZE, то состояние насосов остается неизменным

R29: если у(t) есть РS,  есть NB, то u(t) = f2(·);

R30: если у(t) есть РS,  есть NM, то u(t) = f2(·);

R31: если у(t) есть РS,  есть NS, то u(t) = f2(·).

2. Если скорость уровня нулевая или положительная ZE, PS, то следует включить дополнительный насос j

R32: если у(t) есть РS,  есть ZE, то u(t) = f 0(·) È ui(t) = 1;

R33: если у(t) есть РS,  есть PS, то u(t) = f 0(·) È ui(t) = 1.

3. Если скорость средняя и большая положительная, то следует включить все насосы

R34: если у(t) есть РS,  есть PМ, uj(t) = 0, то u(t) = f4(·);

R35: если у(t) есть РS,  есть PB, uj(t) = 0, то u(t) = f4(·).

Алгоритм импульсного управления состоит в таком изменении интервалов включения  и выключения  насоса i, при котором скорость изменения уровня  стала близкой нулю или равной ZE. Запишем нечеткие правила, реализующие импульсное управление. При средней и малой отрицательной (положительной) скорости необходимо увеличить на dt = D интервал отключения  (увеличить интервал включения )

: если  есть NM, то (t) = (t – 1) + dt;

: если  есть NS, то (t) = (t – 1) + dt;

: если  есть PS, то (t) = (t – 1) + dt;

: если  есть PM, то (t) = (t – 1) + dt.

Таким образом с привлечением алгоритма импульсного управления удаётся стабилизировать уровень стоков в городском резервуаре с минимальными энергетическими  затратами.