Файл: Резниковский, А. Ш. Управление режимами водохранилищ гидроэлектростанций.pdf

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

ожидание состояний системы по данным наблюдений. Так, следуя теории оптимальных фильтров Калмана—Бюси, можно модифициро­ вать исходную систему уравнений каскада с учетом погрешности в наблюдениях и затем производить регулирование по этой новой системе уравнений. Таким образом устраняется влияние искаженных измерений на определение правил управления работой каскада водо­ хранилищ ГЭС.

Целью управления является нахождение такого оптимального марковского управления U 0 (x, t ) , которое доставляет максимум ма­ тематическому ожиданию интеграла экономии затрат в системе:

Оптимальное управление ищется в предположении о его линей­

де К (/) и К (0 матрица и вектор управления.

Функцию экономии затрат также примем линейной функцией своих аргументов (суммы среднеиптервальных мощностей ГЭС), т. е.

U . x . f ( x . U )

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

Ц ? ш . налагаемые на управляющие функции t/t неэнергетическими участниками комплекса.

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

Выбор же порядка штрафной функции ставится в зависимость от уровня обеспеченности того или иного ограничения. Именно, за­ давшись некоторыми значениями Р ша* и P mln, можно подобрать необходимый порядок штрафа так, чтобы эти обеспеченности огра­ ничений выполнялись.

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

Таким образом, управление U 0 понимается оптимальным еще и в смысле непревышения заданных уровней обеспеченности ограниче­ ний в течение периода регулирования (0; Г,,).

Если учесть все ограничения штрафными функциями второй степени, то в силу сделанных выше предположений относительно

Здесь квадратичная и линейная формы (U , GU ) и ( U , G,) есть результаты совместного разложения по степеням //< в функции эконо­ мии затрат З эк, а также в штрафных функциях, налагаемых на управ-

ные функции по состояниям системы (объемам водохранилищ). Кро­ ме того, матрицы G, Е, С обеспечивают отрицательную определен­ ность своим формам и имеют соответственно размерность п х п , n x U

/X/; Gi, Ci — n и /-мерные векторы; © (/)— некоторая известная функция времени.

Далее, используя аппарат, разработанный в [Л. 86], можно найти нелинейное матричное дифференциальное уравнение Риккати,

Матрица Р (/) и вектор R (/) находятся из матричного уравнения Риккати, которое для краткости запишем в общем виде

164

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

де (Т’ч) = (*.—х т , s ( х — х т ) ) .

° к к

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

матрица, обеспечивающая отрицательную полуопределенность штра-

фа 'Ы Гк).

Заметим, что выбор периода регулирования, т. е. момента Гк, является самостоятельной задачей в проблеме долгосрочного регу­ лирования. Так, при сезонном регулировании стока всегда можно указать по времени точку, в которой водохранилище будет при всех обстоятельствах наполнено до отметки ИПУ.

Для случая же многолетнего управления можно, например, за­

дальнейшее изменение перестанет заметно влиять на оптимальное управление за интересующие пас ближайшие годы периода развития

циальноеуравнение рассматриваемой системы и решая для него обратное уравнение Колмогорова [Л. 69], можно найти плотность распределения вероятностей состояний каскада и проверить выпол­ нение заданного уровня обеспеченности ограничений по х. При не­ обходимости производится корректировка порядка штрафной функ­ ции. То же, хотя и несколько сложнее, можно проделать и с другими ограничениями.

Численное решение уравнений (15) производится на ЦВМ по стандартным программам для решения систем обыкновенных диф­ ференциальных уравнений. Размерность системы (15) в общем слу­ чае равна /г+ ((/+ 1)1/2). Возможности же такой ЦВМ, как БЭСМ-4, позволяют решать системы размерности порядка 250.

Рассмотренная методика применялась для отыскания оптималь­ ного (в вышеуказанном смысле) управления каскадом из трех гидро­ узлов (ГЭС № 12, 13 и 16 на рис. 5-3). Боковая приточность аппроксимировалась марковскими процессами, имеющими логнор­ мальное распределение вероятностей. Соответствующие уравне­ ния (15) решались'на ЦВМ. БЭСМ-4 методом Рунге—Кутта. Время расчета на ЦВМ не превышает нескольких минут, что дает надеж­

165

ды на возможность использования данного метода для сложных си­ стем.

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

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

ПРИЛОЖЕНИЕ II

ОСНОВЫ ОТБОРА РАСЧЕТНЫХ ПЕРИОДОВ РЕГУЛИРОВАНИЯ ПРИ ИСПОЛЬЗОВАНИИ НАБЛЮДЕННЫХ И ИСКУССТВЕННЫХ ГИДРОЛОГИЧЕСКИХ РЯДОВ

Отбор расчетных периодов производится: 1) в целях уменьше­ ния объема вычислительных работ, 2) по условию задачи, когда она должна решаться в некоторой определенной гидрологической обста­ новке.

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

При любой расчетной обеспеченности гарантированной отдачи

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

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

Для оценки среднемноголетней выработки энергии ' гидроузла после определения величины гарантированной отдачи и построения

166

Сй йравила распределения избытков; 2) выполняются расчеты регу­ лирования отдачи по всему имеющемуся или отобранному более короткому ряду периодов и подсчитывается среднемноголетняя вы­ работка энергии. Правила распределения избытков могут быть по­

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

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

Отбор расчетных гидрологических периодов производится из всего имеющегося наблюденного или искусственного ряда. Для этого из ряда выделяется некоторая представительная совокупность пе-

167