Материал: Модели управления функционированием систем подачи и распределения воды. Панов М.Я., Петров Ю.Ф

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Zj H j q j , j J

 

 

 

 

*

r

H j Zj

H j , j J

R

 

 

 

, j Jf (p) ,

 

Zj

 

 

 

 

r

f

J

r

f

;

( )

J ( )

( )

J ( )

Jf R

Jr R

Jf R ;

 

 

здесь Hj q j – характеристика насоса, регулятора в составе РФ и метасистемы; H*j – фиксированный пьезометрический напор в резервуарных энергоузлах РФ

и метасистемы; Zj – фиксированный узловой потенциал в энергоузлах-стоках метасистемы (геодезическая высота расположения потребителя).

В замкнутости систем уравнений (1.28) – (1.30), (1.31) – (1.33) легко убедиться, поскольку сумма (е + μ) выражает полное число узлов БРС, а вместе с r (количество элементарных контуров) равно без единицы общему количеству участков (nl + n2 + nЗ = n). То есть число уравнений соответствует числу неизвестных участковых расходов БРЗ. А объединенные матрицы (1.28) – (1.30), (1.31) – (1.33) являются квадратными.

Рассмотрим «конструкцию» этих матриц применительно к схеме водопроводной сети (рис. 1.1). Упорядоченная нумерация реальных и фиктивных участков БРЗ позволяет сформировать клеточную структуру матриц

(1.28) – (1.33).

Знак слагаемого в матрице цепей определяется совпадением направления течения потока на участке с выбранной положительной ориентацией по цепи. Предельное число независимых цепей равно 11, при 11 узлах с фиксированным потенциалом и 1 узле с фиксированной характеристикой НС, 4 участка не охвачены цепями, причем это участки, замыкающие цепи в контуры. Любая новая цепь сверх 11 образует контур, то есть является линейно зависимой. Матрица независимых контуров повторяет аналогичную матрицу в составе модели невозмущенного состояния, поскольку фиктивные участки не образуют новых контуров.

При переходе от модели невозмущенного к модели возмущенного состояния, вследствие присоединения фиктивных участков, энергоузлы РФ утрачивают свой статус ЭУ, превращаясь в промежуточные узлы ветвления (поз. 1 – 6, 8, рис. 1.1), кроме узлов НС, 7, 9, поскольку к этим ЭУ – источникам вода поступает из окружающей среды. В новых ЭУ – стоках (поз. 10 – 18, рис. 1.1) сформированы ГУ I рода, в форме узловых потенциалов, равных соответствующим геодезическим уровням, что является существенным признаком модели возмущенного состояния СПРВ.

36

Рис. 1.4. Цепной подграф бинарного структурного графа модели полноразмерной СПРВ: (1–11) – предельное число независимых цепей; 7, 9 – 18 – узлы с фиксированным потенциалом (ГУ I рода); НС – узел с фиксированной характеристикой (ГУ II рода)

Цепной подграф (рис. 1.4), построенный на вершинах (ЭУ с фиксированным потенциалом) и дугах (независимых цепях), также представляет собой плоский (планарный) структурный орграф (вариант "дерева"), к которому можно применить условие Л.Эйлера о соотношении дуг, вершин и циклов с учетом того, что по определению независимой цепи число циклов r = 0 [22, 51, 63]:

Pц

е 1.

(1.34)

 

 

С целью уяснения правой части подсистемы уравнений независимых цепей (1.31), применительно к рис. 1.1, ниже приведем структуру векторматрицы, матрицы смежности и цепного подграфа бинарного СТГ модели полноразмерной СПРВ на рис. 1.4.

37

Структура вектора–матрицы (рис. 1.1):

 

 

ZHC HHC

(qHC )

 

 

Z

 

H

 

*

 

 

 

 

 

 

7

 

7

 

 

 

 

 

Z9 H9 *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

Z11

 

 

 

 

 

H

Z

 

 

 

 

 

.

 

e 1

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

Z13

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

Z15

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

17

 

 

 

 

 

 

 

Z18

 

 

 

 

 

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

Нетрудно убедиться в том, что объединенная топологическая матрица

A K C

бинарной сетевой структуры в (1.31) – (1.33) является

квадратной, а задача моделирования потокораспределения – определенной.

Структура топологической матрицы смежности БРС (рис. 1.1):

38

1.4. Реструктуризация системы водоснабжения на основе модели возмущенного состояния

1.4.1. Имитационное моделирование варианта реструктуризации кольцевой системы водоснабжения

Проиллюстрируем на примере реальной водопроводной сети, (рис. 1.5), моделирование одной из возможных задач реконструкции системы, а именно реструктуризацию сети, то есть уменьшение числа элементарных колец путем постепенного перекрытия участка 2–11 с помощью дистанционно управляемого дросселя (УД) (рис. 1.5, а). Система оснащена стальными трубами с двумя источниками питания: насосной станцией с параллельным подключением центробежных насосов Д200/95 и Д630/95 в составе НС г.Губкина, Белгородской области и водонапорной башней (поз. НС и ВБ, рис. 1.5). Исходная информация по системе представлена в табл. П.5, П.6.

Процедура реструктуризации должна производиться с сохранением в силе всех сетевых законов и в итоге трехкольцевая сеть, рис. 1.5а должна преобразоваться в двухкольцевую, рис. 1.5б.

Для выделения в составе структурного графа водопроводной сети – узла питания, отождествляющего насосную станцию, необходимо провести

эквивалентирование двух параллельно

подключенных насосов: 1-й насос

Д200/95 с напорной характеристикой

H a q2

c

0, 0058q2

113 и 2-й насос

 

 

 

 

 

 

 

 

1

 

1

1

1

 

 

1

 

Д630/95 с характеристикой H

2

a

2

q2

b

q

2

c

2

0, 00271q2

0, 797q

2

33, 68 .

 

 

2

2

 

 

 

 

 

2

 

 

Напорная характеристика эквивалентного насоса согласно частным условиям эквивалентирования источников питания (2.18) [55] может быть представлена квадратичным трехчленом

Hэ aэ q1 q2 2 bэ q1 q2 cэ .

Рис. 1.5. Расчетная схема системы водоснабжения:

а) до реструктуризации; б) после реструктуризации; – дистанционно управляемый дроссель; – -неуправляемый дроссель

39

Используя расчетные зависимости [55, (2.41)] для определения коэффициентов напорной характеристики центробежного эквивалентного насоса запишем для рассматриваемого случая:

 

 

 

а q3

а

q3

0,001283 ; b

 

 

b

q2

 

0, 4855;c

 

 

c q c

q

2

46,97.

а

 

1 1

 

2

2

 

 

2

2

 

 

1 1 2

 

Э

q q

2

3

Э

q

q

2

2

Э

q1 q2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Таким образом, в узле питания, отождествляющем насосную станцию, реализуются граничные условия II рода с аналитической характеристикой

HНС 0, 001283QНС2 -8 0, 4855QНС-8 46, 97 .

В целом система водоснабжения (рис. 1.5, а), представляет собой плоский (планарный) структурный орграф, содержащий 12 энергоузлов, 11 – с ГУ I рода (10 узлов-стоков и 1 узел-источник, поз. ВБ) и 1 узел-источник (поз. НС) с ГУ II рода.

Для моделирования потокораспределения необходимо составить 11 цепных (по числу узлов без единицы с определенными формами ГУ), 3 контурных и 11 узловых (по числу узлов с незаданным узловым потенциалом) балансовых уравнений, соответствующих 25 неизвестным (по числу участков БСТГ) расчетным расходам. Для определения НС вводится дополнительное уравнение напорной характеристики, приведенное выше.

Для моделирования возмущенного состояния СПРВ (рис. 1.5, а) используется модель (1.31) – (1.33). Возмущения в систему вносятся с помощью дистанционно управляемого дросселя, размещенного на участке 2 – 11, коэффициенты сопротивлений остальных участков остаются неизменными.

Для выведения из эксплуатации участка 2 – 11 (например, с целью проведения ремонтных работ, нанесения противокоррозионного защитного покрытия и т.п.) задается сопротивление этого участка (вместе с сопротивлением УД) следующим алгоритмическим условием:

S(к)

S(к 1)

S(к 1)

 

 

 

 

(к)

S(к)

S(к)

 

 

 

S

(к)

, при этом

 

S

0,1

.

2 11

2 11

2 11

 

2 11

 

2 11

2 11

2 11

 

Задаваемое число итераций K = 106.

Таблица 1.1 Результаты моделирования процесса реструктуризации системы водоснабжения

(рис. 1.5)

Обозначение

до внесения возмущения

 

после внесения возмущения

 

участка

Q, л/с

S·105

h, м

Q, л/с

S·105

h, м

 

НС-8

230,5

1,1293

0,6

234,3243

1,1293

0,62007

 

ВБ-3

169,5

0,6961304

0,2

157,8920

0,69613,04

0,17354

 

10-1

25

560,0

3,5

15,71897

560,0

1,3837

участки

2-1

40

93,75

1,5

49,75517

93,75

2,32085

2-11

40

125,0

2,0

1,517·10-2654

1,14·104937

2,6·10-376

9-11

25

800,0

5,0

37,56602

800,0

11,2896

9-10

50

40,0

1,0

40,47656

40,0

0,65534

реальные

3-2

105

31,746032

3,5

75,09339

31,746032

1,79016

3-4

34,5

210,04

2,5

52,79273

210,04

5,85397

4-11

35

244,898

3,0

55,40327

244,898

7,5172

40

Источник: https://studfile.net/preview/16563865/