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

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

При анализе и синтезе БГС возникает проблема учета действия большого количества факторов и необходимости быстрого получения надежного результата, что составляет основное противоречие современного моделирования. В этом случае БГС традиционно подвергается структурной декомпозиции на РФС и метасистему. РФС – ограниченный в пространстве фрагмент, представляющий интерес с точки зрения анализа протекающих в нем процессов. Сечение по внутренним вершинам графа (как процедура выделения РФС) переводит последние в разряд энергоузлов (ЭУ), через которые осуществляется энерго- и массообмен между РФС и метасистемой, сохраняющей при этом статус БГС.

Для адекватного описания гидравлических процессов в РФС в энергоузлах формируются граничные условия. Для анализа невозмущенного состояния традиционно используется определенная форма ГУ: а) в виде априорно заданных (достоверных) значений узловых потоков qj или потенциалов Hj (ГУ I рода); б) в виде функционально и параметрически определенной взаимосвязи Hj=H(qj) (ГУ II рода – паспортные характеристики регуляторов, насосов и т.д.). Достоверность определенных форм ГУ обеспечивает корректность постановки задач. Анализ возмущенного состояния, то есть прогноз последствий любых структурных или параметрических возмущений в РФС, например аварийное отключение участков, присоединение новых потребителей, замена труб других диаметров и т.д., относится к области задач реконструкции. В этом случае ГУ II рода, являясь заданной характеристикой элемента, обладает "поглощающей" (то есть устойчивой к возмущениям) способностью, в то время как параметры в составе ГУ I рода, в силу связности СТГ, отклоняются от фиксированных значений, что равносильно их информационной утрате. В итоге задача анализа потокораспределения переводится в разряд неопределенных, состояние системы описывается уже прямоугольной матрицей, что является признаком ее вырожденности.

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

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

21

и т.п.), нагруженное и ненагруженное резервирование, аппроксимационнотопологические методы оптимального синтеза и так далее. То есть область моделирования возмущенного состояния систем значительно обширнее, и уместнее было бы рассматривать невозмущенное состояние как частный случай возмущённого [13, 31, 38, 40, 66, 68 и др.].

Приступая к решению поставленной задачи развития декомпозиционного подхода, рассмотрим механизм выделения РФС из состава гидравлической системы, полагая, что на эту процедуру не накладываются ограничения. Известно [2], что традиционной является схема отдачи ЦП потребителям исключительно через узлы, причем это касается и различных вариантов упрощения расчетных схем с использованием участковых путевых отборов. Поэтому естественной границей РФС должен считаться узел расчетной схемы, который в этом случае приобретает статус ЭУ [40, 54], то есть узла, через который осуществляется обмен ЦП между РФС и метасистемой.

Для стационарных состояний РФС (установившееся

потокораспределение) в составе условий однозначности начальные условия отсутствуют и остаются лишь граничные условия (ГУ), которые для j-гo ЭУ можно в общем случае представить зависимостью

где Hj ,qj – узловой потенциал и расход через ЭУ j;. ξj – функция, задающая взаимосвязь между потенциалом и расходом.

Введем следующую классификацию граничных условий, исходя из функции ξj, обобщающую накопленный опыт математического моделирования ГС:

1. ГУ I рода – фиксируется (задается) один из параметров (1.17) в узле, Hj Hj или q j q j второй определяется по результатам анализа

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

2.ГУ II рода – определён вид функции ξj, a Hj ,qj считаются неизвестными, но если один из них найден, то по соотношению (1.17) определяется и второй. Примером ГУ II рода являются напорные характеристики нагнетающего оборудования (насосов, компрессоров), паспортные характеристики регуляторов. Специфической формой граничных условий II рода можно считать функциональные характеристики подсистем при "кибернетическом" моделировании.

3.ГУ III рода предполагаются известными по форме, но не определены

по содержанию. Иными словами, Hj ,qj считаются неизвестными, а для функции ξj установлен только ее вид. Впервые граничные условия III рода рассмотрены в работах [54, 55].

4.ГУ IV рода заранее не определены ни по форме, ни по содержанию и

генерирования граничных условий IV рода является метод диакоптики [38]. Здесь в узлах "разрезания" больших гидравлических цепей (БГЦ) на совокупность малых (МГЦ) граничные условия формируются посредством итеративного подбора характеристик источников давления и расхода вместо элементов цепи пересечений. С точки зрения ГУ диакоптику можно рассматривать как метод преобразования граничных условий от IV рода ко II роду.

Следует подробнее остановиться на модели невозмущенного состояния, поскольку она обстоятельно изучалась на протяжении длительного времени, ещё до появления компьютерных технологий и именно с этой модели начиналось развитие процесса моделирования ГС. Структурный граф модели невозмущенного состояния сформирован реальными линиями; в энергоузлах, ограничивающих РФС, заданы ГУ I и II рода. То есть РФС, являющийся автономным объектом для моделирования, называется расчётной зоной. Иными словами, границы РФС и РЗ невозмущенного состояния совпадают.

Прежде всего рассмотрим “конструкцию” топологических матриц на примере водопроводной сети (рис.1.1).

Рис. 1.1. Бинарный структурный граф водопроводной сети:

нс – насосная станция; Rl, R2 – контррезервуары; г – участки с реальными трубопроводами в составе РЗ, f – фиктивные трубопроводные линии в составе эквивалента АП;

НС, 1,2,…,

9 – энергоузлы РФС (РЗ невозмущенного состояния); НС, 10, 11,…,18 – энергоузлы РЗ возмущенного состояния

Напомним, что РЗ, содержащая в своём составе только реальные трубопроводные линии, является автономным объектом моделирования и отображается структурным графом, построенным на реальных дугах (позиции r, рис.1.1). Резервуарные узлы (поз. R1, R2, рис.1.1), водонапорные башни могут ,в зависимости от режима водопотребления, функционировать как в

23

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

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

Все участки РЗ, кроме НС-1, несут путевую нагрузку, следовательно, согласно методу перевода переменного путевого расхода в постоянный к каждому узлу присоединяется фиктивный (фиксированный по определению) узловой отбор, то есть все узлы, включая и промежуточные, являются энергоузлами. Однако последние относятся к множеству ЭУ с незаданным узловым потенциалом, так как в них неизвестен пьезометрический напор, относительно которых и сформирована топологическая матрица [А]. Нейтральные узлы ветвления (χ) также относятся к узлам с фиксированным (нулевым) отбором и неизвестным потенциалом.

К узлам с фиксированным потенциалом, соответствующим ГУ I и II рода, относятся узлы НС, 6, 7, 8, 9, причем в узле НС сформированы ГУ II рода в форме характеристики насоса (эквивалентной характеристики насосной станции), в остальных узлах – ГУ I рода.

Число независимых цепей определяется из условия Эйлера для плоских графов, р = е – 1 (при r = 0 в составе цепного подграфа); в данном случае е = (НС, 6, 7, 8, 9). В этой задаче размер объединенной (квадратной), топологической матрицы в левой части (1.14) – (1.16) составляет 13 х 13, что обусловлено определенным характером граничной информации в составе модели потокораспределения.

Граничные условия I рода в виде фиксированных узловых потенциалов H отборов (притоков) q включены в состав матриц-столбцов к узлам с

фиксированным отбором отнесены также и узлы ветвления с нулевым отбором (здесь Hj Zj H*j , где H*j – фиксированный пьезометрический напор узла j).

Граничные условия II рода содержат в своей основе напорные характеристики насосов в составе насосных станций, и эти (заданные) характеристики насосов в составе насосных станций формируют дополнительные функциональные связи, не включенные в систему уравнений

(1.14) – (1.16).

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

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

24

участки

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

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

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

25