Методы математического моделирования ЭС. Прогнозирование электропотребления. Оптимизация краткосрочных и оперативных режимов

Содержание

Слайд 2

Общие принципы прогнозирования физических процессов


Общие принципы прогнозирования
физических процессов

Слайд 3

Интуитивные (экспертные) методы прогнозирования Базируются на интуитивно-логическом мышлении. В основе разработки

Интуитивные (экспертные) методы прогнозирования

Базируются на интуитивно-логическом мышлении. В основе разработки прогноза

лежат мнение и профессиональный научный опыт экспертов.
Используются:
когда невозможно учесть влияние многих факторов из-за значительной сложности объекта прогнозирования;
при наличии высокой степени неопределенности информации.
Среди интуитивных методов широкое распространение получил метод экспертных оценок.
Сущность метода экспертных оценок заключается в том, что в основу прогноза закладывается мнение специалиста или коллектива специалистов, основанное на профессиональном, научном и практическом опыте в сочетании с количественными методами оценки и обработки получаемых результатов.
Слайд 4

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

Формализованные (фактографические) методы прогнозирования

Сущность экстраполяционного метода заключается в изучении устойчивых тенденций

развития объекта прогноза и переноса их на будущее.
Различают формальную и прогнозную экстраполяцию
Формальная базируется на сохранении в будущем прошлых и настоящих тенденций развития объекта прогноза.
При прогнозной фактическое развитие связано с гипотезами об исследуемом объекте и о влиянии различных факторов в перспективе на него.
Моделирование предполагает конструирование модели на основе предварительного изучения объекта или процесса, выделения его существенных характеристик или признаков.
Предполагают использование в процессе прогнозирования экономико-математических моделей, которые описывают исследуемый объект в виде математических зависимостей и отношений.
В электроэнергетике чаще всего применяется различное сочетание указанных методов.
Слайд 5

Общая математическая модель нагрузки ЭЭС Графики нагрузки в ЭЭС являются последовательностями

Общая математическая модель нагрузки ЭЭС

Графики нагрузки в ЭЭС являются последовательностями наблюдений

или расчетных значений, показывающих изменения мощности в течение определенного периода времени. В суточных, недельных и годовых графиках отображается периодичность процесса изменения мощности нагрузки, связанная с режимом работы людей, сменой дня и ночи, недельными циклами и сезонными изменениями в течение года.
В энергосистемах имеются фактические данные почасовой (или по-получасовой) нагрузки с разной дискретностью и глубиной архива, причем за много лет. Поэтому чаще других применяются методы статистической экстраполяции.
В общем случае, не выделяя интервалы прогнозирования, модель нагрузки можно представить в виде трех составляющих временного ряда
где Q(t) – тренд – устойчивые систематические изменения;
S(t  – периодическая составляющая – колебания относительно тренда;
U(t) – нерегулярная (случайная) составляющая.
Слайд 6

Общая математическая модель нагрузки ЭЭС Ретроспектива временного ряда Можно выделить характерные

Общая математическая модель нагрузки ЭЭС

Ретроспектива временного ряда

Можно выделить характерные периодичности потребления:
Суточные

максимумы нагрузки в неделе;
Недельные максимумы нагрузки в месяце;
Месячные максимумы нагрузки в году.
Такие периодичности характерны не только для максимальных, но и для любых однотипных точек графика.
Слайд 7

Общая математическая модель нагрузки ЭЭС Для выделения трендовой составляющей часто используют

Общая математическая модель нагрузки ЭЭС

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

модель до третьего порядка включительно. Иногда при выделении тренда предварительно применяют процедуру сглаживания, которая устраняет периодическую и случайную составляющие.
Слайд 8

Краткосрочное прогнозирование электропотребления в ЭЭС Учитывая, что в энергосистемах имеются фактические

Краткосрочное прогнозирование электропотребления в ЭЭС

Учитывая, что в энергосистемах имеются фактические данные

нагрузки с разной дискретностью и глубиной архива, причем за много лет, чаще других применяется метод статистической экстраполяции, иногда называемый календарный метод.
Для имеющейся на интервале ретроспективы (назад) статистической информации находится некоторая аппроксимирующая зависимость
, называемая трендом, которая используется для экстраполяции вперед на интервал упреждения. Аналитическое выражение тренда обычно имеет вид полинома (линейного, квадратичного, трехстепенного). Пусть необходимо на некоторый час предстоящих суток спрогнозировать значение нагрузки . Имеются значения нагрузки в тот же час за несколько предшествующих однотипных суток. Например, известны нагрузки в 12 часов Т предшествующих вторников и требуется спрогнозировать нагрузку на 12 часов нового вторника. Расположим эти данные в статистический ряд:
где Т – количество наблюдений:
Слайд 9

Краткосрочное прогнозирование электропотребления в ЭЭС Для простоты составим линейное уравнение регрессии

Краткосрочное прогнозирование электропотребления в ЭЭС
Для простоты составим линейное уравнение регрессии (или

тренда):
,
где а и b – искомые коэффициенты (параметры) уравнения регрессии; t – независимая переменная (время, календарная дата, номер равноудаленного измерения, зависит от интервала упреждения и дискретности архива данных). Эта зависимость называется статистической, так как она неоднозначна и зависит от объема статистической информации – числа наблюдений Т.
Слайд 10

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

Краткосрочное прогнозирование электропотребления в ЭЭС

Необходимо найти такую прямую, чтобы точки лежали

как можно ближе к ней. Решаем эту задачу методом наименьших квадратов, т.е. определением таких значений а и b, при которых сумма квадратов отклонений была бы минимальной. Для этого отыскиваем минимум функции
,
которую можно представить как
Чтобы найти, при каких значениях а и b функция S достигает минимума, продифференцируем ее и приравняем частные производные нулю:
Слайд 11

Краткосрочное прогнозирование электропотребления в ЭЭС В результате получим два уравнения решив

Краткосрочное прогнозирование электропотребления в ЭЭС

В результате получим два уравнения
решив которые

относительно а и b , получаем параметры уравнения регрессии. После достоверизации его путем вычисления среднего значения и дисперсии (формулы не приводятся) и сравнения их с критериями математической статистики, можно вычислить прогнозное значение нагрузки подстановкой в уравнение момента упреждения
Ниже показаны примеры определения прогнозного значения активной мощности методом линейной и квадратичной аппроксимации тренда:
Слайд 12

Общая математическая модель нагрузки ЭЭС После вычитания из X(t) трендовой составляющей

Общая математическая модель нагрузки ЭЭС

После вычитания из X(t) трендовой составляющей получается

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

Временной ряд без трендовой составляющей

Слайд 13

Общая математическая модель нагрузки ЭЭС Если имеется N результатов наблюдений за

Общая математическая модель нагрузки ЭЭС

Если имеется N результатов наблюдений за период

T (N = 12 в годовом цикле, N = 7 – в недельном и N = 24 – в суточном), то периодическая модель процесса может быть представлена рядом Фурье
где n – количество частот, включенных в модель, а
основная частота гармонического ряда.
Коэффициенты полигармонического ряда a(k) и b(k) находятся известными методами спектрального анализа.
Слайд 14

Краткосрочное прогнозирование электропотребления в ЭЭС - основная частота гармонического ряда. Коэффициенты

Краткосрочное прогнозирование электропотребления в ЭЭС
- основная частота гармонического ряда.
Коэффициенты полигармонического

ряда a(k) и b(k) находятся известными методами спектрального анализа.
В общем случае наивысшая частота гармонического разложения дискретного ряда, называемая частотой Найквиста, определяется половиной интервала между наблюдениями, в случае годового статистического ряда N = 12, n = 24.
Дисперсия, учитываемая i-й гармоникой:
Суммарная дисперсия:
Как правило, первые три гармоники описывают до 90 % всей дисперсии.
Слайд 15

Система рынков в электроэнергетике Рынок электроэнергии: обеспечивает краткосрочную надежность и экономическую

Система рынков в электроэнергетике

Рынок электроэнергии:
обеспечивает краткосрочную надежность и экономическую эффективность
формирует наиболее

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

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

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

Слайд 16

Модель оптового рынка: секторы РЫНОК ЭЛЕКТРОЭНЕРГИИ Свободные двусторонние договоры (СДД) Договоры

Модель оптового рынка: секторы

РЫНОК ЭЛЕКТРОЭНЕРГИИ

Свободные двусторонние договоры (СДД)
Договоры купли/продажи электроэнергии между

поставщиком и покупателем, условия которых, включая цены, определяются по соглашению сторон
Рынок на сутки вперед (РСВ)
Система отношений в рамках ценовых зон оптового рынка электроэнергии между участниками оптового рынка и ФСК, связанная с поставкой/потреблением электроэнергии в объемах, определенных по результатам конкурентного отбора ценовых заявок на сутки вперед
Балансирующий рынок (БР)
Сфера обращения отклонений от плановых объемов поставки электроэнергии, определенных в результате конкурентного отбора ценовых заявок для балансирования системы и (или) определенных по факту производства/потребления электрической энергии на основе данных коммерческого учета

13

Слайд 17

Прогнозирование электропотребления для оптового рынка электроэнергии «на сутки вперед» РСВ (краткосрочное прогнозирование)


Прогнозирование электропотребления для оптового рынка электроэнергии
«на сутки вперед» РСВ
(краткосрочное прогнозирование)

Слайд 18

Краткосрочное прогнозирование электропотребления в ЭЭС Типовые (характерные) суточные графики электропотребления: а) б) в) г) д) е)

Краткосрочное прогнозирование электропотребления в ЭЭС

Типовые (характерные) суточные графики электропотребления:
а) б) в)

г)
д) е)
Слайд 19

Краткосрочное прогнозирование электропотребления в ЭЭС Недельный график мощности нагрузки ЕЭС РФ

Краткосрочное прогнозирование электропотребления в ЭЭС

Недельный график мощности нагрузки ЕЭС РФ

Слайд 20

Краткосрочное прогнозирование электропотребления в ЭЭС Для суточных графиков нагрузки различают следующие

Краткосрочное прогнозирование электропотребления в ЭЭС

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

активной и реактивной нагрузок Pmax и Qmax;
коэффициент мощности максимума нагрузки cos φmax;
коэффициент неравномерности графика Кн = Pmax / Pmin ;
суточный расход активной и реактивной энергии Wа.сут и Wр.сут;
средневзвешенный за сутки коэффициент реактивной мощности
коэффициенты заполнения суточного графика активной и реактивной энергии
Эти показатели отображают конфигурацию графика потребления. Предполагают, что для однотипных дней они являются условно постоянными.
Слайд 21

Краткосрочное прогнозирование электропотребления в ЭЭС Задача планирования режимов и, в частности,

Краткосрочное прогнозирование электропотребления в ЭЭС

Задача планирования режимов и, в частности, задача

расчета установившегося режима энергосистемы начинается с прогнозирования суточного графика электропотребления. При этом непрерывная кривая мощности аппроксимируется кусочно-линейно с временным шагом дискретности :
В отечественной практике = 1 час, хотя в зарубежных энергосистемах рассматриваются шаги дискретности 30 и 15 мин. На рисунке выделен только интервал текущих суток
Тем самым интегрирование заменяется суммирование часовых расходов энергоносителей - топлива, воды и т.д.
Слайд 22

Краткосрочное прогнозирование электропотребления в ЭЭС Применение календарного метода, в котором для

Краткосрочное прогнозирование электропотребления в ЭЭС

Применение календарного метода, в котором для прогноза

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

Сезонная кривая для 10=00 часа
каждых суток планируемого года

Слайд 23

Краткосрочное прогнозирование электропотребления в ЭЭС Другая модель прогнозирования потребления, также учитывающая

Краткосрочное прогнозирование электропотребления в ЭЭС

Другая модель прогнозирования потребления, также учитывающая предположение

о стабильности приведенных на слайдах 15 - 17 характерных конфигураций и величин для каждых типичных суток года, основана на явной недельной периодичности потребления электроэнергии
Для понедельного потребления электроэнергии с учетом линейного тренда запишем
Здесь годовой цикл разбит на 52 недели, а, b, c – параметры аппроксимации кривой регрессии (тренда), определяемые методом наименьших квадратов, n – количество недель, предшествующих началу ретроспективных данных Э(i)ф (I = 1…N) по отношению к первой неделе периода наблюдений.
Для определения параметров а, b, c минимизируется функция
куда подставляются выражения для Э(i) и соответствующие значения Э(i)ф
Слайд 24

Краткосрочное прогнозирование электропотребления в ЭЭС

Краткосрочное прогнозирование электропотребления в ЭЭС

 

Слайд 25

Краткосрочное прогнозирование электропотребления в ЭЭС Учитывая гипотезу, указанную на слайдах 15

Краткосрочное прогнозирование электропотребления в ЭЭС

Учитывая гипотезу, указанную на слайдах 15 -

17, далее полученная
величина суточной электроэнергии потребления конфигурируется
в соответствии с графиком типового (характерного) дня
прогнозируемой недели. Данная процедура называется «вливание
в график».
Слайд 26

Прогнозирование электропотребления для балансирующего оптового рынка электроэнергии БР (оперативное планирование)


Прогнозирование электропотребления для балансирующего оптового рынка электроэнергии БР
(оперативное планирование)

Слайд 27

Оперативное планирование электропотребления в ЭЭС Все изменения уровня потребления, сетевые ограничения,

Оперативное планирование электропотребления в ЭЭС

Все изменения уровня потребления, сетевые ограничения, состояние

генерирующего оборудования невозможно спрогнозировать заранее с абсолютной точностью. При управлении режимом работы энергосистемы в реальном времени необходимо компенсировать возникающие отклонения от режима, запланированного сутки назад. Расчет и доведение до объектов управления графиков генерации при фактическом управлении ЕЭС в режиме реального времени реализуется в рамках балансирующего рынка.
В течение суток Системный оператор многократно проводит формирование прогнозов спроса (потребления), актуализацию расчетной модели с учетом изменившихся системных условий. На основании данной информации проводятся конкурентные отборы ценовых заявок поставщиков, обеспечивающие экономическую эффективность загрузки станций и требования к надежности. Критерием конкурентного отбора в реальном времени является минимизация стоимости поставки электроэнергии, а также меры, делающие невыгодными самовольные отклонения от плановых значений и стимулирующие выполнение команд СО по управлению режимом
Слайд 28

Оперативное планирование электропотребления в ЭЭС Учитывая природу формирования нагрузки потребителями в

Оперативное планирование электропотребления в ЭЭС

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

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

Оперативное планирование электропотребления в ЭЭС При этом задача оперативного прогноза нагрузок

Оперативное планирование электропотребления в ЭЭС

При этом задача оперативного прогноза нагрузок решается

методом аппроксимации значений нагрузки на интервале моделирования τmo нелинейной аддитивной функцией с учетом данных краткосрочного прогноза. Интервал моделирования τmo включает в себя интервал ретроспективных данных τro и интервал упреждения оперативного прогноза τuo. Учет текущих данных в суточном разрезе позволяет также уточнить краткосрочный (суточный!) прогноз и, как следствие, улучшить оперативный прогноз.
Для определения параметров прогнозной модели в момент времени t0 используются фактические значения нагрузки текущего ретроспективного интервала τro, а также часовые значения краткосрочного прогноза нагрузки на интервале упреждения τuo. Ожидаемое значение нагрузки, взятое из результатов краткосрочного прогноза нагрузки, назовем “точкой притяжения”.
Слайд 30

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

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

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

Оперативный прогноз нагрузки

Слайд 31

«Скользящий» прогноз для нужд БР Учет ретроспективных данных за прошедшие сутки

«Скользящий» прогноз для нужд БР

Учет ретроспективных данных за прошедшие сутки от

t0 позволяет уточнить текущее значение «точки притяжения». Для этого применяется метод скользящих суток - циклическая аппроксимация значений нагрузки на суточном интервале ретроспективы относительно текущего момента времени t0. В качестве модели нагрузки на суточном интервале τmc используется так же разложение в ряд Фурье.
Метод скользящих суток обеспечивает расчет “точек притяжения” и коррекцию краткосрочного прогноза с использованием текущей информации о поведении нагрузки на суточном интервале.
Слайд 32

Математические модели комплексной оптимизации режимов электроэнергетических систем


Математические модели комплексной оптимизации режимов электроэнергетических систем

Слайд 33

Постановка задачи оптимизации режима ЭЭС Под задачей оптимизации текущего режима энергосистемы

Постановка задачи оптимизации режима ЭЭС

Под задачей оптимизации текущего режима энергосистемы или

энергообъединения понимают наивыгоднейшее распределение генерируемых активных и реактивных мощностей между электростанциями, а также другими регулируемыми источниками реактивной мощности – синхронными компенсаторами, управляемыми шунтирующими реакторами, устройствами FACTS, которому отвечает минимум эксплуатационных издержек И на производство, передачу и распределение электрической и тепловой энергии в топливном (сегодня чаще в стоимостном) выражениях:
Слайд 34

Постановка задачи оптимизации режима ЭЭС . ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ Рассмотрим математическую основу одного

Постановка задачи оптимизации режима ЭЭС

.

ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ

Рассмотрим математическую основу одного из многочисленных методов

нахождения экстремума целевой функции одной и нескольких переменных
Слайд 35

Безусловный экстремум

Безусловный экстремум

Слайд 36

Безусловный экстремум

Безусловный экстремум

Слайд 37

Безусловный экстремум

Безусловный экстремум

Слайд 38

Безусловный экстремум

Безусловный экстремум

Слайд 39

Безусловный экстремум

Безусловный экстремум

Слайд 40

Относительный экстремум и методы его нахождения

Относительный экстремум и методы его нахождения

Слайд 41

Относительный экстремум и методы его нахождения

Относительный экстремум и методы его нахождения

Слайд 42

Распределение активных мощностей между ТЭС с приближенным учетом потерь Рассмотрим применение

Распределение активных мощностей между ТЭС с приближенным учетом потерь

Рассмотрим применение

метода неопределенных множителей Лагранжа для часто встречающейся задачи распределения активных мощностей между ТЭС с приближенным учетом потерь в электрической сети.
Слайд 43

Распределение активных мощностей между ТЭС с приближенным учетом потерь Функция Лагранжа

Распределение активных мощностей между ТЭС с приближенным учетом потерь

Функция Лагранжа

Слайд 44

Распределение активных мощностей между ТЭС с приближенным учетом потерь грузки. В

Распределение активных мощностей между ТЭС с приближенным учетом потерь

грузки. В то

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

Распределение активных мощностей между ТЭС с приближенным учетом потерь

Распределение активных мощностей между ТЭС с приближенным учетом потерь

Слайд 46

Распределение активных мощностей между ТЭС с приближенным учетом потерь

Распределение активных мощностей между ТЭС с приближенным учетом потерь

Слайд 47

Постановка задачи комплексной оптимизации режимов ЭЭС ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ

Постановка задачи комплексной оптимизации режимов ЭЭС

ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ

Слайд 48

Постановка задачи комплексной оптимизации режимов ЭЭС

Постановка задачи комплексной оптимизации режимов ЭЭС

Слайд 49

Постановка задачи комплексной оптимизации режимов ЭЭС ЖЖЖЖЖЖЖЖЖЖ

Постановка задачи комплексной оптимизации режимов ЭЭС

ЖЖЖЖЖЖЖЖЖЖ

Слайд 50

Постановка задачи комплексной оптимизации режимов ЭЭС

Постановка задачи комплексной оптимизации режимов ЭЭС

Слайд 51

Современные методы оптимизации Метод неопределенных множителей Лагранжа без дополнительных «ухищрений» не

Современные методы оптимизации

Метод неопределенных множителей Лагранжа без дополнительных «ухищрений» не позволяет

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

Современные методы оптимизации

Современные методы оптимизации

Слайд 53

Слайд 54

Современные методы оптимизации

Современные методы оптимизации

Слайд 55

Заключение Интересно отметить, что в градиентном методе в расчете производных от

Заключение

Интересно отметить, что в градиентном методе в расчете производных от потерь

мощности и в расчете градиента участвует та же самая матрица Якоби, что и в расчете установившегося режима методом Ньютона. Поэтому градиентный метод получил наиболее широкое распространение.
Кроме того, имея возможность формировать в одной программе различные целевые функции, с помощью такой программы, например, ПК «Bars-Lincor», можно решать разнообразные задачи: эквивалентирование, утяжеление, расчет режима при неноминальной частоте, оптимизация режима по реактивной мощности, оптимизация режима по активной мощности (по ХОП, по ХОПС, по тарифам и по рыночным ценовым характеристикам), ввод режима в допустимую область и т.д. При этом изменяются только целевые функции, а вычислительная схема остается такой же самой.