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

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

соответствует схеме, когда дроссельные элементы установлены на ответвлениях к энергоулам. Поскольку от каждого энергоузла РЗ запитывается множество разнородных потребителей, их режим потребления определяется конкретной дроссельной характеристикой.

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

QDj

j SDj , j ID

,

(2.1)

 

 

 

где QDj – расход воды через дроссель j; SDj – коэффициент гидравлического сопротивления дросселя j.

Если множество ID дросселей контролирует потребление воды через множество ЭУ – стоков, то можно утверждать, что дроссельные характеристики являются явной формой взаимодействия между режимом водопотребления и гидравлической настройкой дросселей. Как отмечалось выше, это зависит от схемы расположения дросселей.

Дроссельные характеристики (2.1) могут быть синтезированы из модели возмущенного состояния (1.31) – (1.33), то есть (2.1) является объектом моделирования, а ее конфигурация определяется схемой размещения дросселей, структурным графом водоснабжающей системы и не зависит от режима водопотребления конкретной сетевой системы.

2.2. Формализация дроссельной характеристики в рамках прямого анализа системы водоснабжения

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

Моделирование дроссельных характеристик водораспределительной системы основывается на линейной модели потокораспределения возмущенного состояния РСВ. С этой целью проведем линейные преобразования МВС (1.31) – (1.33) с учетом установленных на ряде участков МПГС дроссельных элементов. Обозначим через Ij множество участков в составе независимой цепи j, контура j, а также инцидентных узлу j, причем j Irj Ifj IrjD IfjD . Здесь Irj , Ifj – множество реальных и фиктивных участков,

исключая участки с присоединенными дросселями соответственно; IrjD , IfjD

- множество участков с дросселями, установленными на реальных (РФС) и фиктивных (АП) линиях соответственно.

Вначале проведем линейные преобразования цепного уравнения для наиболее общего случая.

66

Уравнение независимой цепи в составе МПГС с учетом возможной установки дросселей как на реальной, так и фиктивной ее части:

ZjH H*jH ZjK H*jK

sgnSiQi

sgnSiQi ,

 

i Irj Ifj

i IrjD IfjD

где ZjH , H*jH , ZjK , H*jK – геодезический и пьезометрический напоры в начальном

и конечном узлах независимой цепи j; sgn – оператор присвоения знака слагаемому участка: (+) – если направление потока на участке совпадает с положительным направлением по цепи, (–) – в противоположном случае.

Введем отклонения параметров, обусловленные изменением настройки дроссельных элементов в составе цепи j:

ZjH

H*jH H*jH ZjK

H*jK H*jK sgnSi Qi Qi

 

 

 

i Irj Ifj

 

sgn Si Si Qi

Qi .

 

i IrjD

IfjD

 

Используя разложение в ряд Тейлора, при условии Qi << Qi , Si << Si ,

получаем после неучета слагаемых высоких порядков малости, кроме линейных:

ZjH H*jH H*jH ZjK H*jK H*jK

sgn SiQi SiQi 1 Qi

 

i Irj Ifj

sgn SiQi SiQi 1 Qi Qi Si .

 

i IrjD IfjD

 

В соответствии с постановкой задачи возмущенного состояния РСВ

пьезометрический напор H*jK 0 и H*jK

0 . Отклонение H*jH может быть

обусловлено изменением расхода через насос в соответствии с его напорной характеристикой из-за изменения сопротивления дросселей и общего сопротивления сети. Но в данном случае мы анализируем поведение системы с источниками - резервуарами (водонапорными башнями), для которых можно принять с достаточной степенью точности, что H*jH 0 .

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

h*j

sgn SiQi SiQi 1 Qi

sgn SiQi SiQi 1 Qi Qi Si ;

 

i Irj Ifj

i IrjD IfjD

или с учетом уравнения независимой цепи:

67

sgn SiQi 1 Qi

sgn SiQi 1 Qi Qi Si 0.

i Irj Ifj

i IrjD IfjD

Разделим и умножим каждое из слагаемых на выражение для потери напора соответствующего участка цепи j: hij SijQij , после чего поучим цепное

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

sgn hi

 

i

sgn hi

 

i hi

 

i 0.

 

Q

Q

S

 

i Irj Ifj

i IrjD IfjD

(2.2)

 

 

 

 

 

 

 

 

Контурное уравнение для контура j с учетом того, что в составе контуров отсутствуют фиктивные участки, согласно постановке задачи:

sgnSiQi sgnSiQi 0.

i Irj

i IrjD

Контурное уравнение после внесения возмущений, обусловленных изменением сопротивления дроссельных элементов:

sgnSi Qi Qi sgn Si Si Qi Qi 0.

i Irj

i IrjD

Разложение в ряд Тейлора с удержанием линейных членов приводит к линейной форме записи контурного уравнения:

sgn SiQi SiQi 1 Qi sgn SiQi SiQi 1 Qi Qi Si 0.

i Irj

i IrjD

Перепишем последнее уравнение с учетом уравнения независимого контура j, приведенного выше:

sgn SiQi 1 Qi sgn SiQi 1 Qi

Qi Si 0.

i Irj

i IrjD

 

Разделим и

умножим каждое слагаемое

последнего равенства на

выражение hij SijQij , в результате чего получается контурное уравнение в относительных отклонениях:

sgn hi Qi sgn hi Qi hi Si 0.

i Irj

i IrjD

(2.3)

 

 

Узловое балансовое уравнение из состава МВC:

68

sgn Qi

sgn Qi sgn q*j Qj 0,

i Irj Ifj

i IrjD IfjD

где Q j – фиктивный узловой отбор, обусловленный наличием путевой нагрузки на участках; q*j – сосредоточенный (фиксированный) узловой приток (отбор):

при этом ( +) в случае притока к узлу, (–) – в случае оттока от узла.

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

sgn Qi Qi

sgn Qi Qi sgn q* Qj 0.

i Irj Ifj

i IrjD IfjD

Узловое балансовое уравнение в относительных отклонениях:

sgn Qi

 

i

sgn Qi

 

i

sgn Qi

 

i 0.

 

Q

Q

Q

 

i Ir

If

i Ir

If

i I

 

 

 

(2.4)

j

j

jD

jD

 

j

 

 

Уравнение (2.4) получено в предположении, что возмущения в настройке дроссельных элементов не отражаются на значениях путевых составляющих расчетных расходов участков, а влияют на отклонения только транзитной составляющей.

Линейная модель возмущенного состояния РСВ с выделением в качестве определяющих факторов коэффициентов гидравлического сопротивления дроссельных элементов, полученная на основе (2.2) – (2.4), в матричном виде приведена ниже:

Сn1 p

T

 

hn1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

0

 

Q

n1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hn1D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qn1D 1

 

 

Сn1D p

 

 

0

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hn 2

 

 

 

 

 

 

 

 

 

 

 

Сn 2 p

 

 

 

 

 

 

 

 

 

 

Qn 2 1

 

 

 

0

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сn 2D p

 

 

0

 

0

 

 

0

 

hn 2D Qn 2D 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hn1

0

0

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

hn1D

 

 

0

 

0

 

Sn1D 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 ;

 

 

 

 

 

 

 

 

0

 

0

 

 

hn 2

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

0

 

hn 2D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2D 1

 

 

 

 

 

 

 

 

(2.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

69

 

 

 

T

hn1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kn1 r

0

 

 

 

 

0

 

 

 

 

0

 

 

 

 

Q

n1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hn1D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qn1D 1

 

Kn1D r

 

 

0

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

 

 

 

0

 

0

 

 

h

 

 

 

 

 

0

 

 

 

Q

 

 

 

n 2 r

 

 

 

 

n 2

 

 

 

 

 

 

n 2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

On 2D r

 

 

0

 

0

 

 

 

 

0

 

 

 

hn 2D

Qn 2D 1

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hn1(d)

 

 

0

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

hn1D

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sn1D

1

 

 

0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

hn 2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

0

hn 2D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

An1 m

 

 

Qn1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

0

 

Q

n1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

An1D m

 

0

 

Qn1D

 

 

0

 

 

 

 

0

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1D 1

 

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

An 2 m

 

 

 

0

 

0

 

Qn 2

 

 

 

 

0

 

 

 

Qn 2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2D 1

 

 

 

An 2D m

 

0

 

0

 

 

0

Qn 2D

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.6)

(2.7)

где n = nl + n1D + n2 + n2D – общее число участков БСТГ; nl, n1D, n2, n2D –

число реальных участков, исключая участки с управляемыми дросселями;

число

реальных участков с присоединенными управляемыми дросселями;

число

фиктивных участков, число фиктивных участков БСТГ с

 

 

 

 

присоединенными УД соответственно;

hi – потери напора участка i;

S

Dr ,

 

 

 

 

 

 

S

f

– относительные

отклонения

коэффициентов

гидравлического

 

D

 

 

 

 

 

 

сопротивления дроссельных

элементов,

расположенных

на реальных и

фиктивных участках БСТГ соответственно; Q QQ, S SS – определяемый и определяющий факторы соответственно; Q , S –

абсолютные отклонения соответствующих параметров; Q, S – текущие значения тех же параметров.

Линейная модель (2.5) – (2.7) учитывает структуру БСТГ, в общем случае закольцованную в рамках реальной части БСТГ и разветвленную (незакольцованную) в пределах МАП, условия невырожденности матрицы (2.5)

(2.7), вытекающие из аналогичных условий невырожденности матрицы (1.31)

(1.33), позволяют получить решение этой системы уравнений, которое в матричном виде приведено ниже:

70

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