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