Файл: Баясанов, Д. Б. Автоматизированные системы управления трубопроводными объектами коммунального хозяйства.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.10.2024

Просмотров: 127

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Выражение функционала (3.177) может быть с успехом ис­ пользовано для решения задач оптимизации при стацио­

нарных и нестационарных режимах газоснабжения в АСУ

этими процессами. В случае динамических процессов ис­

пользуют аппроксимации нестационарных режимов газо-

Рис. 31. Блочная схема специализированного вычислительного уст­ ройства

а — блоки вычислений; б — блоки задания постоянных величин

снабжения в виде смены стационарных

состояний

этих

процессов.

Вышеуказанные задачи

решаются

при

следу­

ющих ограничениях: Plmax (t) > Р х

>

Р 3

(t)\ Р г (t) >

Р ц ( 0 ;

Р 3 ( 0 ^ Р 3 min (^)-

 

 

 

 

Пока не разработаны критерии, позволяющие оцени­

вать непосредственно как в информационном, так и в режи­ ме замкнутого контура автоматизированного управления

213

функционирование таких сложных и больших систем, как

системы газоснабжения. Поэтому предложенные решения

являются начальным этапом к постановке и реализации

более сложных задач оптимизации в описываемой области. Для реализации задач управления в АСУ может быть ре­

комендовано специализированное вычислительное устрой­

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

управления режимами газоснабжения в общем комплексе АСУ как элемент вышеприведенного функционала. Блоч­

ная схема этого устройства показана на рис. 31. Эта схема

состоит из блоков вычислений z, T cv, S, QBbIx. газ, 3 и С.

При вычислении такие параметры, как протяженность га­

зопроводных линий, коэффициенты гидравлических сопро­

тивлений, диаметры труб и т. п., могут быть приняты по­

стоянными. Другие параметры, как температуры tlt ts, frp, давления газа P lt Р 2, 3 3 Должны вводиться в устройство

непосредственным измерением и передачей показаний на вход вычислительного устройства. Эти данные могут вво­

диться в устройство дискретно диспетчером, который мо­

жет получить их и по телефону с мест контроля. Вычисли­

тельное устройство может эксплуатироваться совместно

сЭВМ при внедрении АСУ.

§7. ЭЛЕМЕНТЫ с т а т и с т и ч е с к о г о

МОДЕЛИРОВАНИЯ в АСУ

Внастоящее время при рассмотрении технологических

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

сания происходящих в них процессов, которые, как уже

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

тора во времени. Этот накопленный статистический мате­

риал может быть использован для составления статистиче­

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

систем, так'и коммунальных хозяйств в комплексе. Каж­ дая отрасль коммунальных хозяйств и входящие в нее

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

214


еПолным основанием сгруппировать их в отдельные группы.

Кпримеру, все технологические трубопроводные системы

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

каждого отдельного случая представлена на рис. 32. Эта

схема построена на основе обобщения опыта по управлению

объектами газораспределения и раскрывает характер ос­

новных взаимосвязей, существующих при организации

АСУ этими объектами. В элемент сравнения (или в один из

w

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

/ _ объект управления; 2 — математическая

модель

объекта; 3 — блок

прог­

нозирования добычи и подачи газа; 4 — блок выбора

оптимальных отклонений

регулируемых параметров-, оптимизированных

в технологическом смысле;

5 ^

блок выбора оптимальных отклонений регулируемых параметров, оптимизи­ рованных по материальным затратам; 6 — сумматор

отделов предприятия) поступает плановое задание на пла­

нируемый

период

времени в

виде

временной функции

Y 3

( 0

= у\, у\

,..., уп, где у*

(г =-

1

— задание

на

подачу

газа на i-й месяц. В

блоке сравнения функция

Y 3 (/)

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

газа в

систему

У*пр (t) = у[пр,

у\пр

 

УпР-

 

 

Если соответствующие значения

Ayt =

(г/? — г/гпр) Ф О,

то принимаются

решения по

выработке таких регулиру-

 

 

 

 

 

—Ф

 

^

^

 

ющих

воздействий

X (/) = х * (/),

х\ (/),..., Хт (t),

где

i =

1 ,....,

т — количество регулируемых параметров, кото­

рые обеспечивали хотя бы условие Ayt =

0. Для этого не­

обходимо знать не только величину Ауи но и качественную

215


взаимосвязь подачи газа с определяющими ее параметрами

xt (t) при известных затратах Hj, приходящихся на измене­

ние единицы каждого фактора. Выработанные в блоке 4,

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

лируемых параметров Але,- (t), реализующие условия Ay t ==

=0 , минимизируются, например, по материальным затра­

там в блоке 5. В результате этого получаются оптимальные,

в экономическом смысле, характеристики отклонений регу-

лируемых параметров Ах, (t). Эти отклонения суммируются

с соответствующими значениями хг (t), выдавая командную

информацию в виде вектора функции хг ( t) на объект управ­

ления.

С выхода объекта управления получается функция

подачи газа Y (t) = у 1} г/2,..., уп, где y t (i = 1, ...., п)

фактическая подача газа в t-м месяце, которая в силу слу­ чайных возмущений w, действующих на объектах управле­

ния, отличается от плановой подачи газа у*.

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

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

газа Y (t) с определяющими ее параметрами х (t), т. е,

Y (t) =■ Ах (/).

Необходимо построить математическую модель объекта управления и рассмотреть близость этой модели к объекту.

Будем

искать

оптимальный

оператор А * математической

модели

Y* (t)

= A *x{t),

где Y * (t) = yt, yt ,..., у*

рассчитанные значения функции подачи газа в систему. Статистические данные можно представить в виде таблицы:

Yi

А Д

'^21- ••

A j i ........

Х т ±\

 

 

Х 1Л

Х „ . . .

X j i ............

Х т 2>

 

Уп

X ln

Х 2П. . .

X j n ........... х тп,

 

где столбцы представляют собой показатели подачи

газа,

а строки — месячные данные,

характеризующие эти

пока-

216


затели. Предполагая, что элементы первого столбца приведенной таблицы взаимосвязаны линейно с элементами

других столбцов, можно записать:

Y

= b0 -f- fe1x 1 +

+ b2x 2 + ... + Ьтхт, где bj+1 (/ =

1 ,

.... т) — коэф­

фициенты регрессии, удовлетворяющие принятому крите­ рию оптимальности. Эти коэффициенты определяются для

каждой конкретной системы газоснабжения. Для нахожде­

ния статистической зависимости можно использовать метод

наименьших квадратов. Основываясь на этом методе, мож­

но определить коэффициенты bJ+1, выражающие уровень

подачи газа как функции регулируемых факторов для раз­

личных объемов обрабатываемого статистического материа­

ла: по ряду лет эксплуатации системы газоснабжения. По­

лучение таких зависимостей, учитывающих различную

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

рактер влияния каждого учтенного фактора на показатели

подачи газа в систему газоснабжения. Полученные числен­

ные значения коэффициентов показывают различную сте­

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

если коэффициент bj+1 имеет положительное значение, то

соответствующий фактор имеет тенденцию повышать подачу

газа и наоборот, если bj+l имеет знак минус, то соответ­

ствующий фактор имеет тенденцию понижать подачу газа

в систему. Полученные статистические модели процессов

необходимо проанализировать. Оценка их связана со сле­ дующими погрешностями при составлении моделей. Это,

во-первых, погрешность, связанная с выбором вида операто­

ра А * и учитываемых факторов, т. е. погрешность, связан­ ная с идентичностью модели реальному процессу, во-вторых, погрешность определения постоянных коэффициентов, в-третьих, погрешность, связанная с ошибками измерений отдельных факторов.

В качестве примера рассмотрим погрешность выбран­

ного оператора А * математической модели относительно

оператора реального объекта управления А. Ясно, что

требование близости А * к А будет разумным при условии выполнения требования близости функции Y * (t) на выходе

модели к функции Y (t) на выходе объекта. Близость А * к А связана с выбранным критерием, зависящим от конкрет­

ной задачи. Когда критерий задан, то говорят о близости

Y * (t) к Y (t) в смысле этого критерия, или об оптималь­

ности этого критерия определения модели объекта. Сравним значения динамического ряда tji (i — 1 , ..., п) фактической

2J7


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

Y*

(i =

1,

я). Естественно, что рассчитанные значения

У*

во

многом

не совпадают 'с фактическими

значениями

Y t. Эти расхождения, очевидно, вызываются

причинами,

действие которых уже не имеет отношения к изменениям

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

значениями

параметров называются остаточными zt =

= yl yi,

а величина s =

— стандартной

 

V

п — 1

ошибкой оценки. Величина s будет являться критерием иден­

тичности оператора А * к А

и измеряется теми же единица­

ми, что и значение уровня

подачи газа в систему газоснаб­

жения, и может сравниваться с этими значениями. Можно

рассчитать и относительные ошибки по соотношению s0TH =

= - 1 0 0 %, где у

— среднемесячная подача газа в систему.

У

что величина ошибки зависит от периода

Следует отметить,

времени, за который строится математическая модель объек­

та, т. е. она возрастает с увеличением продолжительности

динамического ряда. Это объясняется в основном тем, что с увеличением динамического ряда все существеннее сказы­

вается нелинейность взаимосвязей, которая и приводит к увеличению ошибки.

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

при ориентировочных расчетах. Существует метод Брауна,

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

стикой этой взаимосвязи может служить коэффициент кор­

реляции между этими величинами. Корреляционный анализ

статистического материала нередко является первым эта­

пом извлечения полезной информации об объекте управле­

ния статистическими методами. При расчете выборочного

парного коэффициента парной корреляции, осуществляе-

'мой по формуле;

218