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

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

F

Z11 S2 11 Qfz2 11 2

 

 

Z11

S2 11 Qf2 11

2

2 Z12 S4 12 Q4fz 12

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

2

 

2

 

 

 

 

 

fz

 

 

2

 

 

 

 

f

 

2

 

 

2

 

Z12

S4 12 Q4 12

 

 

 

 

Z13

S6 13 Q6 13

 

 

Z13

S6 13 Q6 13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z14 S8 14 Q8fz 14 2

Z14 S8 14 Q8f 14

2

2

Z15 S10 15 Q10fz 15

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z S

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qf

 

 

2

 

Qfz

Qfz

 

Qfz

 

Qfz

 

Qfz

Qfz

 

;

 

 

 

 

10 15

10 15

 

 

 

 

 

2 11

4 12

 

 

 

6 13

8 14

10 15

 

Б

 

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

или в общем случае для системы водоснабжения любой конфигурации, с расположением управляемых дросселей по схеме рис. 2.1, 2.2:

 

 

 

 

 

2

 

 

2

 

2

 

 

 

 

 

 

 

 

Sf

Qfz

 

 

Qf

 

 

 

 

 

 

 

 

F

 

 

Sf

 

 

Qfz

Qfz ,

 

 

 

j

j

 

j

j

 

 

 

 

j

j

 

 

j JH

 

 

 

 

 

 

 

 

 

j J

 

j J

 

(2.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Sfj – коэффициент гидравлического сопротивления фиктивного участка j;

λ – неопределенный множитель Лагранжа; Jη, Jπ – множество стоков и источников МПГС соответственно; JH – множество узлов-стоков (потребителей) с фиксированным потенциалом.

Целевая функция (2.20) сформирована на основе сопоставления полных пьезометрических напоров в энергоузлах расчетной зоны, определяемых по соответствующим двум расходам Qfj и Qfzj . Вторая группа слагаемых (2.20)

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

Принципиальным является вопрос выбора независимых переменных в (2.20), каковыми являются априорно заданные значения расходов через фиктивные линии БСТГ ( Qfzj , j JH ), хотя и подверженные влиянию

субъективных факторов пользователя. Фактические расходы Qfj , j JH

независимыми не являются, будучи "связанными" моделью потокораспределения. Такой выбор упрощает структуру экстремума (минимума) целевой функции, позволяя сформировать систему нормальных уравнений:

F

 

 

f fz

 

2

f f 2

 

 

f fz

 

 

 

2

Sj Qj

 

 

Sj Qj

 

Sj Qj

 

0, j JH

fz

 

Qj

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

(2.21)

86

Исключение λ приводит к условию

f

2

fz

 

fz

 

2

f

2

const, j JH

 

 

Sj

 

Qj

Qj

 

Qj

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(2.22)

 

 

 

 

 

 

 

 

 

 

Благодаря (2.22) удается синтезировать mн – 1 дополнительные независимые связи, которые позволяют сформировать систему нормальных уравнений, размерностью mн – 1:

f

 

2

fz

 

fz

 

2

f

 

2

f

 

2

 

fz

 

 

fz

 

2

 

f

 

 

2

 

 

 

 

 

S1

 

Q1

Q1

 

Q1

 

S2

 

Q2

Q2

 

Q2

 

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

2

fz

 

fz

 

2

f

 

2

f

 

2

 

fz

 

 

fz

 

2

 

f

 

 

2

 

 

 

 

 

S1

 

Q1

Q1

 

Q1

 

S3

 

Q3

Q3

 

Q3

 

 

0;

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.................................................................................

 

 

 

 

 

 

 

 

S1f 2 Q1fz Q1fz 2 Q1f 2 Sfm

 

2

Qfzm

 

Qfzm

 

2

Qfm

 

2

 

 

H

H

H

H

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Оптимизационная подмодель, построенная на основе МНК и содержащая систему нормальных уравнений (2.23), представлена ниже в матричном виде:

где

SrDj(k)

 

 

 

Ef

f

Qf

 

Ef f

Qfz

;

 

 

 

 

 

 

 

(d)

 

 

(d)

 

 

 

 

 

(2.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fj

Sfj 2 Qfzj

Qfj ; fj Sfj 2

Qfzj 2 . При ограничениях:

если

SDjr(k) SDj0r , то

Sr

; если

Sf (k)

Sf

,

то

 

Sf (k) Sf

,

где

 

S

 

коэффициент

 

Dj0

 

Dj

Dj0

 

 

 

Dj

Dj0

 

 

 

Dj0

 

 

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

Оптимизационная (нормальная) матрица

Ef

, составленная из

 

 

 

 

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

(mН - 1)× mН.

87

Ниже приведен один из вариантов нормальной матрицы для МПГС

(рис. 2.1).

Вследствие большой нелинейности нормальных уравнений могут возникнуть проблемы вычислительного характера, поэтому рассматривалась вторая версия целевой функции, построенная не на законе Бернулли, а на сопоставлении расходов воды через фиктивный участок, что дало возможность существенно уменьшить нелинейность:

F Qfzj Qfj

2

 

Qfzj

 

 

 

Qfzj

.

 

 

 

 

 

 

 

j J

 

 

 

 

j JH

 

j J

 

(2.25)

 

 

 

 

 

 

 

 

Условия минимума целевой функции, обозначаемой нами как функция

ошибок:

 

 

 

 

 

 

 

 

F

2 Qfzj Qfj 0 .

 

 

 

fz

 

 

 

Q

j

j JH

 

 

 

 

(2.26)

 

 

 

 

 

 

 

 

После исключения λ получаем условие:

 

 

 

 

 

 

Qfzj Qfj const, j JH

,

 

(2.27)

из которого вытекает система нормальных уравнений, размерностью mн - 1:

Q1fz Q1f Qfz2

Qf2 0;

 

 

Q1fz Q1f Q3fz Q3f 0;

 

 

 

 

.......................................

.

 

 

 

Q1

Q1

QmH QmH 0;

 

fz

f

fz

f

 

 

 

 

 

 

 

(2.28)

 

 

 

 

 

 

 

 

88

 

 

Как и в предыдущем случае, допускается любое другое, не повторяющееся сочетание фиктивных участков в составе (2.28).

Оптимизационная подмодель, построенная на новой версии целевой

функции (функции ошибок):

 

 

 

 

 

 

Ef Qf

Ef Qfz .

 

 

 

 

 

 

 

(2.29)

 

 

 

 

 

 

В дальнейших исследованиях, проводившихся на основе вычислительных экспериментов, изучалось влияние (2.24) и (2.29) на конфигурацию дроссельных характеристик. Более глубокое различие между этими подмоделями не исследовалось, поскольку такая цель в работе не ставилась.

Оптимизационно-топологическая модель РСВ для рассматриваемого случая схемы замещения АП может быть сформирована объединением топологической (1.31) – (1.33) и оптимизационной (2.24) моделей (для первой версии целевой функции) в единую модель. Подобное "механическое слияние" двух моделей обусловлено, как отмечалось ранее, переводом определяющего вектора SD в разряд определяемых, с возникшим при этом дефицитом уравнений и изменением конфигурации матрицы МВС. Дополнение традиционного числа уравнений теперь уже топологической подматрицы нормальными уравнениями позволяет восстановить квадратную конфигурацию матрицы МВС, но уже на уровне оптимизационно-топологической модели, представленную ниже в матричном виде (рис. 2.1):

 

С

 

T

 

 

h

n1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1 p

 

 

 

 

 

 

 

 

M

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

n1D p

 

 

h

n1D 1

 

;

 

 

 

 

 

 

 

 

 

 

p e

 

e 1

 

 

С

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.30)

 

 

 

n 2 p

 

 

 

 

 

 

n 2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

n1 r

 

T

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1 1

0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

On1D r

 

 

 

0

 

 

 

 

 

 

 

 

O

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

(2.31)

 

 

 

 

 

n 2 r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

n1 m

T

 

 

Q

n1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

An1D m

Qn1D 1

 

;

 

 

 

 

 

A

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

(2.32)

 

 

 

 

 

n 2 m

 

 

 

 

 

n 2 1

 

 

 

 

 

 

89

произошел переход к

 

On1 (m 1)

T

0

 

0

 

0

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

On1D (mH 1)

 

0

 

0

 

0

 

 

 

 

 

 

 

 

0

(m

 

E

n 2 (mH 1)

 

H

 

 

 

 

 

 

 

 

0

 

1)(d)

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Qf

 

 

 

n 2 1

 

 

 

 

 

 

 

 

On1 (m 1)

T

0

 

0

 

0

 

 

0

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

On1D (mH 1)

 

0

 

0

 

0

 

 

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

(m 1)(d)

Qfzn 2 1

 

E

 

 

 

 

 

 

n 2 (mH 1)

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Число участков с дросселями строго соответствует числу без единицы ЭУ – стоков. Возможно уменьшение числа дроссельных элементов при соответствующем снижении ранга матрицы, но с сохранением ее квадратной конфигурации. В этом случае часть потоков воды через ЭУ, не подконтрольная дросселям, остается плохо прогнозируемой. Матрицы (2.30) – (2.33), благодаря упорядоченной нумерации реальных, фиктивных участков и участков с дросселями, имеют блочную структуру.

Для проведения линейных преобразований матрицы (2.30) – (2.33) необходимо выполнить процедуру линеаризации системы нормальных уравнений (2.33), с целью чего перепишем нормальное уравнение из состава (2.23) в форме удобной для преобразований, рассматривая общий случай режима течения. Решение ищется для двух фиктивных участков i и j в составе отдельного уравнения:

Sif 2 Qifz 2 1 Sif 2 Qifz 1 Qif Sfj 2 Qfzj 2 1 Sfj 2 Qfzj 1 Qfj .

Дальнейшие выкладки производятся для варианта расположения дросселей на реальных участках БРЗ. Внесем возмущение в гидравлическую настройку дросселей, в результате чего нормальное уравнение перепишется с

учетом отклонений

расходов, при неизменных коэффициентах Sf ,

Sf

,

 

i

j

 

относящихся к фиктивной части МПГС:

 

 

Sif 2 Qifz Qifz 2 1 Sif 2 Qifz Qifz 1 Qif Qif

 

 

Sfj 2 Qfzj

Qfzj 2 1 Sfj 2 Qfzj Qfzj 1 Qfj Qfj .

 

 

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

сущости возмущающего воздействия: вместо Sif Sfj возмущениям Qifz Qfzj :

90

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