Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методическое пособие 733

.pdf
Скачиваний:
5
Добавлен:
30.04.2022
Размер:
6.36 Mб
Скачать

узла 2) являются энергоузлами, чем и определяется количество независимых цепей.

Рис. 3.6. Расчетная схема разветвленной системы водоснабжения с водонапорной башней: - дистанционно управляемый дроссель

Математическая модель оперативного управления для разветвленного варианта системы водоснабжения (рис.3.6), II версия ЦФ:

 

Сn1 p

T

 

hn1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cn1D p

 

hn1D 1

 

M

 

H

 

;

 

 

 

 

 

 

 

 

 

C

 

 

h

 

 

n 2 p

 

n 2 1

 

p e

 

e 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.40)

Cn 2D p

 

hn 2D 1

 

 

 

 

 

An1 m

T

 

Qn1 1

 

 

 

 

 

 

 

 

 

 

 

 

An1D m

 

Qn1D 1

 

 

 

 

 

 

0 ;

 

 

 

A

n 2 m

 

 

 

Q

n 2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.41)

 

 

 

 

 

 

 

 

 

An 2D m

 

 

Qn 2D 1

 

 

131

 

A

n1 m

T

Qz

 

 

 

 

 

 

n1 1

 

 

An1D m

z

 

 

 

 

Qn1D 1

 

0 ;

 

 

 

A

n 2 m

 

Qz

 

 

 

 

 

 

n 2 1

 

 

 

 

z

 

 

 

 

 

 

 

 

 

An 2D m

Qn 2D 1

 

 

Cn1 p

Cn1D p

Cn 2 p

Cn 2D p

T

Qz

 

 

 

n1 1

 

z

n1D 1

Qz

Qzn 2 1

n 2D 1Q

 

Cn1 p T

Qn1 1

 

 

 

 

 

 

 

 

 

Cn1D p

Qn1D 1

 

 

 

.

 

 

 

C

n 2 p

 

 

Q

n 2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cn 2D p

Qn 2D 1

 

(3.42)

(3.43)

Объединенная матрица (3.40) – (3.43) имеет квадратную конфигурацию, то есть соответствующая система уравнений является определенной, а ее решение – однозначным. В данном случае число уравнений: 11 (цепных) + 13 (узловых балансовых относительно Qi) + 13 (узловых балансовых относительно QiZ ) + 11 (оптимизационных цепных уравнений) = 48. Число неизвестных:

24 Qi 13(Qiz ) 11 SiD 48 . В систему уравнений вместо 11 задаваемых Qiz включены 11 неизвестных коэффициентов SiD УД, призванных обеспечить заданный режим водопотребления. С результатами моделирования (рис. 3.6) можно ознакомиться из табл. 3.7, 3.8.

Как следует из табл. 3.7, 3.8, погрешность исполнения прогноза не зависит от структуры системы, а определяется алгоритмом реализации модели оперативного управления, в основе которой лежит условие (3.33). То есть усовершенствованием алгоритма можно добиться выполнения теоретического прогноза водопотребления.

Условие квадратичности конфигурации объединенной матрицы модели оперативного управления для принятой системы функциональных ограничений является инвариантным по отношению к структуре водопроводной сети, что подтверждается результатами проведенного системного анализа.

Возможность размещения управляемых дросселей на фиктивных ветвях БСТГ вынуждает переходить на более высокий уровень детализации абонентских подсистем с целью разработки специального алгоритма идентификации гидравлических характеристик этих подсистем.

132

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

водоснабжения с водонапорной башней (рис.3.6)

 

Обозначение

 

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

 

 

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

 

 

 

участка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q(QZ), л/с

S·104

 

 

h, м

 

Q (QZ), л/с

 

S·104

 

 

 

h, м

 

 

 

 

 

ВБ-1

 

187,0

0,3660385

 

1,28

 

153,096

0,3660385

 

0,85794

 

 

 

 

1-5

 

10,0

 

268,0

 

 

2,68

 

8,18696

1805,288

 

12,1

 

 

 

 

 

1-9

 

38,0

 

29,293629

 

4,23

 

31,1104

29,293629

 

2,835

 

 

 

 

 

1-2

 

124,0

1,580385

 

2,43

 

101,518

1,580385

 

1,629

 

 

 

 

 

9-8

 

18,0

 

32,098765

 

1,04

 

14,736

421,545

 

9,154

 

 

 

 

 

9-10

 

20,0

 

51,25

 

 

2,05

 

16,374

433,112

 

11,612

 

 

участки

 

2-11

 

32,0

 

37,10937

 

3,8

 

26,198

37,10937

 

2,547

 

 

 

 

2-7

 

35,0

 

40,163266

 

4,92

 

28,654

40,163266

 

3,297

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7-6

 

15,0

 

158,6666

 

3,57

 

12,280

894,378

 

13,488

 

 

Реальные

 

11-12

 

16,0

 

148,828

 

3,81

 

13,099

799,293

 

13,715

 

 

 

3-13

 

20,0

 

140,0

 

 

5,6

 

16,374

530,471

 

14,222

 

 

 

 

2-3

 

57,0

 

7,35611

 

2,39

 

46,665

7,35611

 

1,602

 

 

 

 

 

3-4

 

12,0

 

332,6389

 

4,79

 

9,824

 

1526,599

 

14,734

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1-14

 

15,0

 

1256,444

 

28,27

 

12,280

1902,515

 

28,692

 

 

 

 

5-24

 

10,0

 

2729,0

 

27,29

 

8,187

 

2729,0

 

 

18,291

 

 

 

 

8-15

 

18,0

 

589,506

 

19,1

 

14,736

589,506

 

12,802

 

 

 

 

10-16

 

20,0

 

587,25

 

23,49

 

16,374

587,25

 

 

15,744

 

 

 

 

11-19

 

16,0

 

948,828

 

24,29

 

13,099

1559,9

 

 

26,766

 

 

участки

 

12-17

 

16,0

 

880,078

 

22,53

 

13,099

880,078

 

15,101

 

 

 

6-23

 

15,0

 

953,333

 

21,45

 

12,280

953,333

 

14,377

 

 

 

 

7-22

 

20,0

 

640,5

 

 

25,62

 

16,374

1061,717

 

28,465

 

 

Фиктивные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3-20

 

25,0

 

466,4

 

 

29,15

 

20,467

743,85

 

 

31,161

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4-18

 

12,0

 

1670,833

 

24,06

 

9,824

 

1670,833

 

16,126

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13-21

 

20,0

 

501,25

 

20,05

 

16,374

501,25

 

 

13,439

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.8

 

 

Прогноз и исполнение прогноза режима водопотребления по результатам

 

 

 

 

моделирования системы водоснабжения (рис. 3.6)

 

 

 

 

 

Обозначение

1-14

5-24

8-15

10-16

 

11-19

12-17

7-22

 

6-23

3-20

 

 

4-18

13-21

участка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q Z(O) ,

 

л/с, до

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

внес.

 

 

15,0

10,0

18,0

20,0

 

16,0

16,0

20,0

 

15,0

25,0

 

12,0

20,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

возмущния

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q Z(К) ,

 

л/с,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

прогно

 

 

12,0

8,0

14,4

16,0

 

12,8

12,8

16,0

 

12,0

20,0

 

9,6

16,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

водопотреб.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q Z(К) ,

 

л/с,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

исполнен.прогн

12,28

8,187

14,736

16,374

13,099

13,099

16,374

12,28

20,467

9,824

16,374

оза

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

133

3.4. Алгоритм идентификации гидравлических характеристик управляемых дросселей на ветвях структурного графа абонентских подсистем

Результаты системного анализа процесса оперативного управления СПРВ иллюстрируют функциональную эффективность установки управляемых дросселей на реальных или фиктивных ответвлениях от кольцевых структур БСТГ. Естественно, что процедура идентификации гидравлических характеристик УД относится к дросселям, устанавливаемым на фиктивных ветвях, эквивалентирующих реальные сетевые структуры абонентских подсистем.

Очевидным является то, что режимом водопотребления непосредственно управляет дроссель, присоединенный к реальной ветви АП, который по гидравлическим характеристикам несовместим с УД, установленным на фиктивной линии. То есть необходимо определить область рабочей характеристики реального дросселя, адекватную пределам изменения SD фиктивного УД, полученным из результатов моделирования процесса оперативного управления.

Таким образом, неизбежен обратный переход от фиктивных к реальным сетевым структурам АП, с разработкой специального алгоритма идентификации коэффициентов гидравлического сопротивления УД, установленного на реальных и фиктивных ветвях абонентских подсистем.

Содержательная постановка задачи сводится к следующему: при известных в энергоузле j потенциале, расходе воды через него и коэффициенте SDj фиктивного участка, присоединенного к ЭУj, для двух состояний системы (до и после внесения возмущения) определить значения (для этих же состояний) коэффициента гидравлического сопротивления УД, присоединенного к реальной ветви структурного графа АП. Исходная информация по энергоузлу j считается известной из результатов моделирования системы. Иллюстрация «свертывания» абонентских подсистем, в рамках перехода к более высокому уровню детализации этой задачи, приведена на рис.3.7.

Рис. 3.7. Иллюстрация «свертывания» абонентских подсистем с управляемым дросселем: а – гидравлический эквивалент АП; б – реальная сеть АП; управляемый дроссель

134

Более общим случаем формирования абонентской подсистемы представляется гидравлически связанная сетевая структура, запитываемая от множества энергоузлов РФ. Реструктуризация ее на отдельные независимые фрагменты, по аналогии с рассмотренным выше вариантом, довольно проблематична (рис. 3.8).

Если в предыдущем случае формирования эквивалента АП рекомендуется использовать МПГС (рис. 1.2, б), то в данном случае целесообразным представляется модель ПГС, структурированная более гибкой схемой абонентской метасистемы, согласно рис. 1.2,в.

Абонентская сеть, как и в предыдущем случае, будучи ограниченной энергоузлами с ГУ I рода, является автономным объектом для моделирования. Определение необходимых для идентификации расходов на участках АП возможно решением прямой задачи анализа невозмущенного состояния (1.14) – (1.16). В качестве источников выступают энергоузлы PЗ с известными (из результатов моделирования задачи оперативного управления) узловыми потенциалами. В энергоузлах - стоках также сформированы ГУ I рода в форме эквивалентных геодезических уровней, определенных с учетом разноэтажности коммунально-бытовых и промышленных потребителей. В этой части просматривается аналогия с алгоритмом формирования МВС. Таким образом, алгоритм идентификации гидравлических характеристик УД строится на основе моделирования невозмущенного состояния абонентской сети.

Как было установлено ранее, эта задача имеет однозначное решение, хотя нетрадиционным в ее постановке является отсутствие информации по гидравлическим характеристикам ветвей АП с присоединенными УД, подлежащими идентификации. Это влечет за собой необходимость в дополнительных аналитических связях. Именно в этой части оказываются востребованными частные условия энергетического экивалентирования, обстоятельно изученные во второй главе [55].

На примере абонентской подсети, запитываемой от четырех энергоузлов исследуемого фрагмента системы (рис.3.8), рассмотрим часть алгоритма идентификации, которая строится на моделировании невозмущенного состояния АП. Поскольку решение приводится без привязки к конкретному информационному материалу, оно может быть реализовано для любого состояния (начального или конечного) МПГС.

Фрагмент полноразмерной гидравлической системы (рис.3.8, а) отображается физической моделью (рис.3.8, б). Для любого состояния МПГС, из результатов моделирования процесса оперативного управления,известны потенциалы в энергоузлах РФ 1,2,3,4, узле ветвления 26 и узле – стоке 27, а также известны расходы воды на всех ветвях (рис. 3.8, б).

Потеря энергии на транспортировку воды на фиктивных участках модели АП:

W*

1

 

S1 26Q1 261

S2 26Q2 261

S3 26Q3 261

S4 26Q4 261

S26 27Q26 127 const.

1

 

 

 

 

 

 

135

Условие постоянства W* справедливо для конкретного состояния МПГС, в рамках итерационного моделирования прямой задачи анализа невозмущенного состояния АП.

Рис. 3.8, Идентификации гидравлических характеристик УД:

 

 

 

 

 

а – фрагмент ПГС; б – фрагмент МПГС; В – абонентская сеть;

 

 

 

 

-

 

 

 

управляемый дроссель

 

 

 

 

 

136

 

 

 

 

 

Исходя из этого запишем для абонентской сети (рис. 3.8, в), частное условие энергетического эквивалентирования (2.23) [55]:

1 W* S

Q 1 S

Q 1

S

Q 1

S

Q 1

S Q 1

S

Q 1 S Q 1

 

 

 

 

 

 

 

 

1 5

 

 

1 5

2 6

 

 

2 6

 

 

 

 

 

3 7

 

3 7

 

 

4 8

 

 

 

 

4 8

 

 

 

6 5 6 5

 

 

 

 

 

7 6

 

 

7 6

7 8

 

7 8

 

S

 

Q 1

 

S

 

 

Q 1

S

 

 

 

Q 1

S

Q 1

 

S

 

 

Q 1

S

 

Q 1

S

 

 

Q 1

 

5 14

5 14

 

 

 

 

15 14 15 14

 

16 15

 

16 15

 

 

17 16

17 16

 

 

 

14 21

14 21

 

 

 

 

15 21

 

15 21

16 21

16 21

 

S

Q 1 S

 

Q 1 S

 

 

Q 1 S

Q 1

 

S

Q 1

S

 

Q 1

 

S

Q 1

 

 

5 9

 

5 9

 

 

14 18

 

14 18

6 10 6 10

 

7 11

7 11

 

 

 

8 12

 

 

8 12

15 19 15 19

 

 

 

17 13 17 13

 

 

 

 

S8 17Q8 171

S16 20Q16 120

S21 22Q21 122 ; (*)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или в относительных отклонениях, с учетом постоянства W* и переменности SD

участка АП с присоединенным дросселем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1)(Q1 5h1 5 Q1 5 Q2 6h2 6 Q2 6

Q3 7h3 7 Q3 7

Q4 8h4 8 Q4 8

 

Q6 5h6 5 Q6 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q7 6h7 6 Q7 6 Q7 8h7 8 Q7 8 Q5 14h5 14 Q5 14 Q15 14h15 14 Q15 14 Q16 15h16 15 Q16 15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q14 21h14 21 Q14 21 Q15 21h15 21 Q15 21 Q16 21h16 21 Q16 21 Q17 16h17 16 Q17 16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q5 9h5 9 Q5 9 Q14 18h14 18 Q14 18 Q6 10h6 10 Q6 10 Q7 11h7 11 Q7 11 Q8 12h8 12 Q8 12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q15 19h15 19 Q15 19 Q17 13h17 13 Q17 13

Q8 17h8 17 Q8 17 Q16 20h16 20 Q16 20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q21 22h21 22 Q21 22 ) Q14 18h14 18 Q14 18 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнее уравнение и является той дополнительной аналитической связью, которой не хватает для определения коэффициента гидравлического сопротивления участка с присоединенным УД (S14-18, рис. 3.8. в), полученной из соответствующего частного условия ЭЭ.

Приведем полную систему уравнений идентификации для квадратичного режима течения в трубах (=2,0):

I. Цепные уравнения для независимых цепей

1.S

Q2

 

 

 

S

 

 

 

Q2

(Z H ) (Z H

 

);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 5

 

 

 

1 5

 

 

5 9

 

 

 

 

5 9

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

5

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.S

 

Q2

 

 

 

S

 

 

 

 

Q2

 

S

 

 

Q2

 

 

 

 

 

(Z

 

H ) (Z

H

 

);

 

 

 

 

 

 

 

 

 

 

 

 

1 5

 

 

1 5

 

 

5 14

 

 

 

 

5 14

 

 

 

 

14 18

 

14 18

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

18

 

 

 

18

 

 

 

 

 

 

 

 

 

 

 

 

 

3.S

 

 

Q2

 

 

 

S

 

 

 

 

 

 

Q2

 

 

 

(Z

H

2

) (Z

 

 

H

 

);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 6

 

2 6

 

6 10

 

 

6 10

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.S

 

 

Q2

 

 

 

S

 

 

 

 

 

 

Q2

 

 

 

(Z

H ) (Z

 

 

H

 

);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 7

3 7

 

7 11

 

7 11

 

 

 

 

 

3

 

 

 

 

3

 

 

 

 

 

11

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.S

 

Q2

 

 

 

S

 

 

 

 

 

Q2

 

 

(Z

H

4

) (Z

 

H

 

);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 8

 

 

4 8

 

8 12

 

 

 

8 12

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.S

 

 

Q2

 

 

 

S

 

 

 

 

 

Q2

 

 

S

 

 

Q2

 

 

 

 

 

(Z

4

H

 

) (Z

 

H );

 

 

 

 

 

 

 

 

 

 

 

4 8

 

4 8

 

8 17

 

 

 

8 17

 

 

 

 

 

17 13

 

17 13

 

 

 

 

 

 

 

 

 

4

 

 

13

 

 

 

13

 

 

 

 

 

 

 

 

 

 

 

 

7.S

 

 

Q2

 

 

 

S

 

 

 

 

 

Q2

 

 

S

 

 

Q2

 

 

 

 

S

 

 

 

 

 

Q2

 

 

 

(Z

 

H

) (Z

 

H

20

);

 

 

 

 

4 8

 

4 8

 

8 17

 

 

 

8 17

 

 

 

 

 

17 16

 

17 16

 

 

 

16 20

 

16 20

 

 

4

 

 

4

 

20

 

 

 

 

 

 

 

 

8. S

 

 

 

 

Q2

 

 

S

 

 

Q2

 

 

S Q2

 

 

 

(Z H

20

) (Z H

22

);

 

 

 

 

 

 

 

20 16

 

 

20 16

 

 

 

 

 

 

16 21

 

 

16 21

 

 

21 22

 

 

21 22

 

 

 

 

 

 

 

20

 

 

 

 

 

 

22

 

 

 

 

 

 

 

 

 

 

9.S

Q2

 

 

 

S

 

 

 

 

Q2

 

S

 

 

Q2

 

 

 

 

S

 

 

 

Q2

 

 

 

(Z H ) (Z

H

22

);

 

 

 

 

 

1 5

 

 

 

1 5

 

 

5 14

 

 

 

 

5 14

 

 

 

 

14 21

 

14 21

 

 

 

 

21 22

 

 

 

21 22

 

 

1

 

 

 

1

 

22

 

 

 

 

 

 

 

 

 

10. S

 

 

 

 

 

Q2

 

 

 

 

S

 

 

 

 

Q2

 

 

S

 

 

 

 

 

Q2

 

 

 

 

 

(Z

 

H

20

) (Z

H

 

);

 

 

 

 

 

 

 

 

 

 

16 20

 

16 20

 

 

 

 

16 15

 

 

16 15

 

 

 

15 19

 

15 19

 

 

 

 

 

20

 

 

 

 

 

19

 

19

 

 

 

 

 

 

 

 

11.S

 

 

 

Q2

 

 

S

 

 

 

 

Q2

 

 

S

 

Q2

 

 

 

(Z H

) (Z H );

 

 

 

 

 

 

 

 

 

 

 

 

4 8

 

 

4 8

 

 

 

8 7

 

8 7

 

 

 

 

 

7 11

 

7 11

 

 

 

 

 

4

 

 

 

 

 

 

4

 

 

 

 

11

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

12. S

 

 

 

Q2

 

 

S

Q2

 

 

S

 

 

Q2

 

 

 

 

(Z

 

 

H

 

) (Z

H

);

 

 

 

 

 

 

 

 

 

 

 

 

 

7 11

 

 

7 11

 

 

 

 

 

 

 

7 6

 

 

7 6

 

 

6 10

 

 

 

6 10

 

 

 

 

11

 

 

11

 

 

 

 

10

 

 

10

 

 

 

 

 

 

 

 

 

 

 

II. Контурные уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13.S

 

 

 

Q2

 

 

S

 

 

 

 

 

Q2

 

 

 

S

 

 

 

Q2

 

 

 

 

S

 

 

 

 

 

 

Q2

 

 

 

S

 

 

 

Q2

 

S

 

 

Q2

 

S Q2

S

Q2

0;

6 5

 

 

6 5

 

 

 

5 14

5 14

 

 

 

15 14

 

15 14

 

 

16 15

 

16 15

17 16

17 16

8 17

8 17

7 8

7 8

7 6

7 6

 

14.S

 

 

 

 

Q2

 

S

 

Q2

 

 

S

 

 

 

Q2

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15 14

 

 

 

15 14

 

 

 

 

 

 

14 21

 

 

14 21

 

 

 

15 21

 

 

15 21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15.S

 

 

 

 

Q2

 

S

 

Q2

 

 

S

 

 

 

Q2

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16 15

 

 

 

16 15

 

 

 

 

 

 

15 21

 

 

15 21

 

 

 

16 21

 

 

16 21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

III. Узловые балансовые уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16.Q6 5 Q1 5

 

Q5 9

Q5 14

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

137

17.Q7 6

Q2 6

Q6 10

 

Q6 5

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18.Q3 7

Q7 6

Q7 11

Q7 8

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19.Q4 8

Q8 12

Q8 17

 

Q7 8

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20.Q5 14

Q14 21

Q15 14 Q14 18

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21.Q16 15 Q15 14 Q15 21

Q15 19

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22.Q17 16

Q16 15

Q16 21

Q16 20

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23.Q8 17

Q17 16

Q17 13

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24. Q14 21 Q15 21 Q16 21

 

Q21 22

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IV. Уравнение энергетического эквивалента:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25.S

 

 

Q3

S

 

Q3

S

 

Q3

S

Q3

 

S

 

Q3

 

S

Q3

S

Q3

 

S

 

Q3

 

 

 

1 5

 

1 5

 

2 6 2 6

 

 

3 7

3 7

 

 

4 8

4 8

6 5

6 5

7 6

 

7 6

 

7 8

 

7 8

 

5 14

5 14

 

 

S

 

Q3

 

S

Q3

 

 

 

S

 

Q3

S

Q3

 

S Q3

S

Q3

 

S

Q3

 

 

15 14

 

 

15 14

 

16 15 16 15

 

17 16

 

 

17 16

 

14 21

14 21

 

15 21

 

15 21

 

16 21

 

16 21

 

 

5 9

 

5 9

 

S

Q3

 

S

 

Q3

 

S Q3

 

 

S

 

Q3

S Q3

 

S

 

Q3

 

S

 

Q3

 

 

 

6 10

 

 

 

6 10

 

7 11 7 11

 

 

 

8 12 8 12

15 19

 

15 19

 

17 13 17 13

8 17

8 17

 

 

16 20

 

16 20

 

 

 

S

 

 

Q3

 

S

Q3

 

 

3W*.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21 22

 

21 22

 

14 18 14 18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующая система уравнений в относительных отклонениях, решаемая на отдельной итерации, в рамках алгоритма идентификации:

1. 2h(k 1)

 

 

(k)

2h(k 1)

 

 

 

(k)

 

 

 

 

 

 

 

 

 

 

Q

 

 

Q

0 ;

 

 

 

 

 

 

 

 

1 5

 

 

 

 

 

1 5

5 9

 

 

 

 

5 9

 

 

 

 

 

 

 

 

 

 

2. 2h(k 1)

 

Q

(k)

2h(k 1)

 

Q

(k)

 

2h(k 1)

 

Q

(k)

h(k 1)

 

S

(k)

0 ;

1 5

 

 

 

 

 

1 5

5 14

 

 

 

 

 

5 14

14 18

 

 

14 18

14 18

 

14 18

 

3.2h(k2 61) Q(k2 )6 2h(k6 101) Q(k6 )10 0 ;

4.2h3(k 71) Q3(k)7 2h(k7 111) Q(k7 )11 0 ;

5.2h(k4 81) Q(k4 )8 2h8(k121) Q8(k12) 0 ;

6. 2h(k 1)

 

 

 

 

(k)

2h(k 1)

 

 

 

 

(k)

2h(k 1)

 

 

 

(k)

 

 

 

 

 

 

 

 

 

 

 

Q

 

Q

 

Q

0 ;

 

 

 

 

 

 

 

4 8

 

 

 

 

 

 

4 8

 

 

 

8 17

 

 

 

 

 

 

8 17

 

 

 

 

17 13

 

 

 

 

 

17 13

 

 

 

 

 

 

 

 

 

 

7. 2h(k 1)

 

Q

(k)

2h(k 1)

Q

(k)

 

2h(k 1)

 

Q

(k)

2h(k 1)

 

Q

(k)

0 ;

4 8

 

 

 

 

 

 

4 8

 

 

 

8 17

 

 

 

 

 

 

8 17

 

 

 

 

 

17 16

 

 

 

 

 

 

17 16

 

 

16 20

 

 

 

 

16 20

 

8. 2h(k 1)

Q

(k)

2h(k 1)

Q

(k)

 

2h(k 1)

Q

(k)

0 ;

 

20 16

20 16

 

 

16 21

 

 

16 21

 

 

 

21 22

 

21 22

 

 

 

 

 

 

 

9. 2h(k 1)

 

Q

(k)

2h(k 1)

 

Q

(k)

2h(k 1)

Q

(k)

2h(k 1)

Q

(k)

0 ;

1 5

 

 

 

 

 

 

1 5

 

 

5 14

 

 

 

 

 

5 14

 

 

 

 

14 21

 

 

 

 

 

14 21

 

 

21 22

 

 

 

21 22

 

10. 2h(k 1)

Q

(k)

 

2h(k 1)

Q

(k)

 

 

2h(k 1)

Q

(k)

0 ;

 

16 20

 

16 20

 

 

 

16 15

 

 

 

16 15

 

 

 

 

15 19

 

 

15 19

 

 

 

 

 

 

 

11.2h(k4 81) Q(k4 )8 2h8(k 71) Q8(k)7 2h(k7 111) Q(k7 )11 0 ;

12.2h(k7 111) Q(k7 )11 2h(k7 61) Q(k7 )6 2h(k6 101) Q(k6 )10 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13. 2h(k 1)

Q

 

(k)

 

2h(k 1)

Q

 

 

(k) 2h(k 1)

Q

(k)

2h(k 1)

 

Q

(k)

 

6 5

 

 

 

 

 

6 5

 

5 14

 

 

 

 

5 14

15 14

 

15 14

 

16 15

 

 

 

16 15

 

2h(k 1)

Q

(k)

2h(k 1)

Q

(k) 2h(k 1)

Q

(k)

2h(k 1)

Q

(k) 0 ;

17 16

 

17 16

 

 

8 17

 

8 17

7 8

7 8

 

 

7 6

 

7 6

 

14. 2h(k 1)

 

Q

(k)

 

 

2h(k 1)

Q

(k)

2h(k 1)

 

Q

(k)

 

0 ;

 

 

 

 

 

15 14

 

 

 

 

 

15 14

14 21

 

 

 

 

 

 

14 21

15 21

15 21

 

 

 

 

 

 

 

15. 2h(k 1)

 

Q

(k)

 

 

2h(k 1)

Q

(k)

2h(k 1)

 

Q

(k)

 

0 ;

 

 

 

 

 

16 15

 

 

 

 

 

16 15

15 21

 

 

 

 

 

 

15 21

16 21

16 21

 

 

 

 

 

 

 

16.Q(k6 51) Q(k6 )5 Q1(k51) Q1(k5) Q5(k 91) Q5(k)9 Q5(k 141) Q5(k)14 0 ;

17.Q(k7 61) Q7(k)6 Q(k2 61) Q(k2 )6 Q(k6 101) Q(k6 )10 Q(k6 51) Q(k6 )5 0 ;

18.Q3(k 71) Q3(k)7 Q(k7 61) Q(k7 )6 Q(k7 111) Q(k7 )11 Q(k7 81) Q(k7 )8 0 ;

19.Q(k4 81) Q(k4 )8 Q8(k121) Q8(k12) Q8(k171) Q8(k17) Q(k7 81) Q(k7 )8 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20. Q(k 1)

 

Q

 

 

(k)

Q(k 1)

Q

(k)

Q(k 1)

Q

(k)

Q(k 1)

Q

(k)

0 ;

5 14

 

 

 

 

 

 

 

 

 

5 14

 

14 21

 

 

 

 

 

 

14 21

 

15 14

 

 

 

 

15 14

 

14 18

 

 

 

 

 

14 18

 

 

21. Q(k 1)

 

 

Q

 

(k)

 

Q(k 1)

Q

 

(k)

 

Q(k 1)

Q

 

(k)

 

Q(k 1)

 

Q

(k)

 

0 ;

16 15

 

 

 

 

 

 

 

 

16 15

15 14

 

 

 

 

 

 

 

 

15 14

15 21

 

 

 

 

 

15 21

15 19

 

 

 

 

 

 

15 19

 

 

22. Q(k 1)

 

 

Q

(k)

 

Q(k 1)

 

Q

(k)

 

Q(k 1)

Q

(k)

 

Q(k 1)

 

 

Q

(k)

 

0 ;

17 16

 

 

 

 

 

 

 

 

17 16

16 15

 

 

 

 

 

 

 

 

16 15

16 21

 

 

 

 

 

 

16 21

16 20

 

 

 

 

 

 

16 20

 

23. Q(k 1)

 

Q

(k)

Q(k 1)

Q

(k)

Q(k 1)

Q

(k)

0 ;

 

 

 

 

 

 

 

 

 

 

8 17

 

 

 

 

 

 

 

 

8 17

 

17 16

 

 

 

 

 

17 16

 

17 13

 

 

 

17 13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

138

 

 

 

 

 

 

 

 

 

 

24. Q(k 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

Q(k 1)

 

(k)

 

Q(k 1)

 

(k)

 

Q(k 1)

 

 

 

 

 

(k)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

(k)

 

 

 

 

 

 

 

Q

 

Q

 

Q

 

0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14 21

14 21

 

 

15 21

 

 

15 21

 

 

 

 

16 21

 

 

 

16 21

 

21 22

 

 

 

 

 

 

 

 

 

21 22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25. 3Q(k 1)h(k 1)

 

Q

(k)

3Q(k 1) h(k 1)

 

Q

(k)

 

3Q(k 1) h(k 1)

Q

 

(k) 3Q(k 1) h(k 1)

Q

(k)

 

 

 

 

 

1 5

 

1 5

 

 

 

 

 

 

 

 

1 5

 

 

2 6

 

2 6

 

 

 

 

 

 

 

 

 

 

2 6

 

 

 

 

3 7

3 7

 

 

 

 

 

 

 

 

 

3 7

 

4 8

4 8

 

 

4 8

 

 

 

 

3Q(k 1)h(k 1)

Q

 

(k) 3Q(k 1) h(k 1)

Q

(k) 3Q(k 1) h(k 1)

Q

(k)

 

3Q(k 1) h(k 1)

Q

(k)

 

 

 

 

 

 

 

 

 

6 5

6 5

 

 

 

 

 

 

 

 

 

 

 

6 5

7 6

7 6

 

 

 

7 6

 

 

 

7 8

 

7 8

 

 

 

7 8

 

 

 

 

 

5 14

5 14

 

 

5 14

 

 

 

 

 

 

 

 

 

3Q(k 1) h(k 1)

 

Q

 

 

 

(k)

3Q(k 1) h(k 1)

 

Q

 

(k)

 

 

3Q(k 1) h(k 1)

Q

(k)

 

3Q(k 1) h(k 1)

Q

(k)

 

 

15 14

15 14

 

 

 

 

 

 

 

 

 

 

 

 

 

15 14

 

16 15

16 15

 

 

 

 

 

 

 

 

 

16 15

 

 

 

17 16

17 16

 

 

 

17 16

 

14 21 14 21

 

 

 

 

 

14 21

 

3Q(k 1) h(k 1)

 

Q

 

(k)

3Q(k 1) h(k 1)

 

Q

(k)

 

 

3Q(k 1) h(k 1)

Q

(k)

3Q(k 1) h(k 1)

Q

(k)

 

15 21

15 21

 

 

 

 

 

 

 

 

 

 

 

 

15 21

 

16 21

16 21

 

 

 

 

 

 

 

 

 

 

16 21

 

 

 

5 9

5 9

 

 

 

 

 

 

 

 

 

 

 

5 9

 

 

14 18

14 18

 

 

 

 

14 18

 

 

3Q(k 1)h(k 1)

Q

 

(k) 3Q(k 1) h(k 1)

Q

(k) 3Q(k 1) h(k 1)

Q

 

(k) 3Q(k 1) h(k 1)

Q

(k)

 

 

 

 

 

6 10

6 10

 

 

 

 

 

 

 

 

 

 

 

 

6 10

7 11

7 11

 

 

 

 

 

7 11

 

 

 

8 12 8 12

 

 

 

 

8 12

 

15 19

15 19

15 19

 

 

 

3Q(k 1) h(k 1)

 

Q

 

(k)

3Q(k 1) h(k 1)

Q

(k)

3Q(k 1) h(k 1)

 

Q

(k)

3Q(k 1)

h(k 1)

 

Q

(k)

 

 

17 13

17 13

 

 

 

 

 

 

 

 

 

 

 

 

 

17 13

 

8 17

 

8 17

 

 

 

 

 

 

 

8 17

 

 

 

16 20

16 20

 

 

 

 

 

 

 

 

 

16 20

 

 

21 22

21 22

 

 

 

 

 

 

 

21 22

 

3Q(k 1) h(k 1)

 

 

Q

(k)

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14 18

14 18

 

 

 

 

 

 

 

 

 

 

 

 

14 18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача моделирования процесса управления решается для двух состояний (до и после исполнения прогноза), а излагаемый ниже алгоритм идентификации гидравлических характеристик УД может быть реализован для любого из этих состояний. Для каждого из них известны по результатам моделирования узловые потенциалы в энергоузлах РФ, расходы воды и потери энергии W* абонентской сети, запитываемой от этих ЭУ.

По своей сути алгоритм идентификации сводится к решению системы нелинейных уравнений прямой задачи анализа (включающей условие энергетического эквивалентирования), например уравнений (1) – (25), рис. 3.8.

1. Для формирования начального приближения означенной системы необходимо задаться значением SDi участка i с присоединенным УД (например

S14 – 18, рис. 3.8в).

2. Произвольное значение SDi противоречит условию ЭЭ (уравнению (25)). Поэтому это условие должно быть исключено на данном этапе.

3. Для оставшейся системы уравнений (уравнения (1)–(24), рис.3.8, в), отражающей модель невозмущенного состояния АП с граничными условиями в форме узловых потенциалов, необходимо сформировать начальное приближение. Одним из возможных вариантов решения данной задачи является топологическо-узловой метод с предварительным заданием режимов водопотребления через отдельные энергоузлы-стоки АП, с соблюдением

условий сплошности водопотоков (Q1-26 + Q2-26 + Q3-26 + Q4-26 = Q5-9 + Q14-18 + Q6-10

+ Q7-11 + Q8-12 + Q15-19 + Q16-20 + Q17-13 + Q21-22, рис. 3.8, б,в).

4.Производится решение системы уравнений, сформированной согласно п.2, «стартующее» от начального приближения, определенного согласно п.3. Результат решения этой системы уравнений формирует новое начальное приближение по расчетным расходам в сочетании с заданным значением SDi.

5.От начального приближения п.4 стартует итерационное решение полной системы уравнений (уравнения (1) – (25), рис. 3.8, в), содержащей уравнение энергетического эквивалентирования (25). Результаты решения помимо расходов содержат искомую величину коэффициента гидравлического сопротивления реального УД.

6.По результатам решения идентификационной задачи (см. п.5) для двух состояний МПГС (до и после исполнения прогноза водопотребления)

определяется предел изменений ΔSD (гидравлическая перенастройка) реального УД, обеспечивающий заданный режим водопотребления.

139

Построение алгоритма идентификации с помощью «вложенных» итерационных процедур вызвано необходимостью выполнения на всех его этапах сетевых законов, что гарантирует качество сходимости решения.

В случае установки нескольких дросселей на реальных ветвях АП должно быть задано, в рамках режима водопотребления отношения между расходами воды QD2/QD1, QD3/QD1 и т.д. Это позволяет выразить расходы воды по дроссельным линиям через один из них, не меняя алгоритма идентификации.

3.5. Приближенное решение задачи идентификации гидравлических характеристик управляемых дросселей

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

n 2

SQ1 SD1QD11 SiQi 1,

i 2

где первое слагаемое в первой части равенства отражает потерю энергии на реальном участке АП с присоединенным дросселем; вторая группа слагаемых – общие потери энергии на остальных (реальных) участках АП.

Запишем уравнение энергетического баланса в отклонениях:

n 2

SSQQ1 SD1 SD1 QD1 QD1 1 Si Qi Qi 1 .

i 2

После разложения в ряд Тейлора с удержанием линейных слагаемых:

 

 

 

 

1

 

 

 

1

 

 

 

SSQ

1 QQ

SD1 SD1 QD1

1 QD1 QD1

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

Si Q 1

( 1)Q Q .

 

 

 

 

 

 

 

i

 

 

i i

 

 

 

 

 

 

i 2

После раскрытия скобок и использования более удобной записи в относительных отклонениях получаем:

n 2

1 SQ1 QSQ1 S1 SD1QD11 QD1 SD1QD11 SD1 ( 1) SiQi 1Qi ;

i 2

Исходя из условия энергетического баланса допускается, что изменения потерь энергии, обусловленное изменением участковых расходов, приближенно равны на эквиваленте и реальных трубопроводных линиях АП, то есть

n 2

1 SQ1 Q1 SD1QD11 QD1 1 SiQi 1Qi ;

i 2

140