МОДЕЛИРОВАНИЕ И ОПТИМИЗАЦИЯ В РЕАЛЬНОМ ВРЕМЕНИ УМНОЙ МИКРОСЕТИ ЭЛЕКТРОПИТАНИЯ РАЙОНА С РАСПРЕДЕЛЕННЫМИ ИСТОЧНИКАМИ ЭНЕРГИИ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

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

энергии

П.М. Ведель, Н.Ю. Овчинникова НГТУ, Новосибирск, Россия

Аннотация. В системах электроснабжения все более важную роль играют возобновляемые генераторы энергии, для которых разрабатываются способы управления, учитывающие изменение цены на электроэнергию. В данной работе представлен вариант управления потреблением энергии на основе прогнозирующей экономической модели (Economic MPC), которая минимизирует стоимость эксплуатации микросети с распределенным производством электроэнергии. Контроллер для координации и оптимизации работы микросети в реальном времени, который включает в себя блок прогнозирования потребления электроэнергии и профилей цен на электроэнергию на основе моделей ARIMA, разработан на базе прогнозирующей экономической модели с учетом смешанного целочисленного программирования и протестирован в среде Matlab / Simulink. Было исследовано изменение параметров умной микросети на общую экономию затрат.

Ключевые слова: умная микросеть, Economic MPC, модель ARIMA, ценообразование в реальном времени.

ВВЕДЕНИЕ

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

В связи с этим, в 2015 году в Париже выработано соглашение по объединению усилий всех мировых держав по сдерживанию климатических изменений и регулированию методов генерации энергии. Цели данного соглашения достигаются путем отказа от невозобновляемых источников энергии и переходу к возобновляемым источникам энергии.

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

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

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

1. СОСТАВЛЕНИЕ МОДЕЛИ ОБЪЕКТА УПРАВЛЕНИЯ

Для работы была выбрана эталонная микросеть низкого напряжения касающаяся жилой зоны из CIGRE (международный совет по крупным электрическим сетям) [1]. Эта сеть включает в себя 5 конечных пользователей (частные дома и многоквартирные дома), подключение к главной энергосети, хранилище энергии в виде аккумуляторов, возобновляемые генераторы (солнечные панели и ветровые турбины) и управляемые генераторы (микротурбины и топливная ячейка). Данная микросеть представлена на Рис. 1.

Модель умной микросети была создана в среде Matlab с использованием Simulink / Simscape при помощи решателя, основанного на комплексных амплитудах (phasor).

Подключение к главной сети

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

Рис. 1. Эталонная микросеть низкого напряжения СЮЯЕ

Типовая связь блоков разных библиотек Основная часть, которая связывает БтиНпк (управляющий сигнал мощности) и Simscape (трехфазная сеть) можно описать через формулу:

панели (СП) солнечного

управляются излучения

£ = 3• I 1 (УаЬ -УЬс • г-'120)• 1А

(1)

источника

мощность

/„, - комплексная

К -

где: £ - полная электричества [ВА]; + рА

амплитуда тока фазы А [А]; УаЬ и

линейный напряжения [В] .

Комплексные амплитуды токов других фаз

получаются с помощью умножения IрА на

поворотные коэффициенты е±'120.

Описание неконтролируемых элементов Данная микросеть имеет неконтролируемые элементы, такие как конечные потребители и возобновляемые источники энергии. Данные элементы управляются по реальным историческим профилям за 2018 год для Чешской Республики [2].

Для каждого конечного пользователя был выбран типовой профиль потребления, соответствующий его конфигурации [3, 4].

Солнечные значением

К (/) [Вт / м2 ] :

Рр (') = К К ('). (2)

где: РСП ^) - мощность блока солнечных панелей в данный момент времени [Вт]; КСП -

коэффициент, характеризующий площадь СП, КПД СП и КПД инвертора.

Ветровая турбина (ВТ) управляется значением скорости ветра V (() [м / с]:

Рвт ^) = КВТ (г)3, (3)

где: РВТ (/) - мощность ветровой турбины

в

времени

[Вт]

К

ВТ

параметры

данный момент

коэффициент, характеризующий воздуха и характеристики ВТ.

Описание контролируемых генераторов

Микротурбина (МТ) и топливный элемент (ТЭ) в выбранной микросети - это управляемые генераторы, работающие на натуральном газе. В данной модели рассматривается только электрические характеристики МТ и ТЭ. Стоит учитывать, что электрический КПД ТЭ намного

выше, чем КПД МТ, однако ТЭ имеет меньшую номинальную электрическую мощность [5, 6]:

Рмт (t) = Гш^мт (t), (4)

Ртэ (t) = Гтэ -Qts (t), (5)

где: РМТ (t), РТЭ (t) -мощности, выделяемые МТ и ТО в момент времени t [Вт] ; гМТ, ГТЭ -электрические КПД МТ и ТО при постоянной температуре окружающей среды; QMT (t), QT3 (t) - полное количество энергий подаваемых на МТ и ТО через натуральный газ [Дж / ч].

Хранилище энергии представлено в виде аккумуляторов [7]. В рассматриваемом диапазоне безопасной зарядки/разрядки, состояние заряда (Soc) будет меняться по линейному закону:

( р (t)^ пТ

S°ck+1 = Sock + Pchar (t) - гЛаг I • — -100, (6)

I Г J Cap

при условии:

SoCmin < S0Ck+1 < S0Cmax , (7)

где: Pchar (t) , Pdis (t) - мощности

зарядки/разрядки аккумулятора в момент времени t, представленные как положительные компоненты полной мощности Pbatt (t)[Вт] ; Sock+1, Sock - состояние заряда аккумулятора в момент времени k+1 и k [%]; Socm¡n, Socmax -минимальное и максимальное состояние заряда аккумулятора (соответственно 20% и 90%); Г char, Г ¿и* - КПД зарядки и разрядки блока аккумулятора соответственно; ПТ - шаг дискретизации по времени [ч] ; Cap - ёмкость аккумулятора [Вда • ч].

Ток зарядки/разрядки аккумулятора Im (t)[А] в момент времени t находится как:

I,

Irms (t) =

Г char ' • Pbatt (t ) , Pat, (t )< 0

P

— ■ PT-Pat, (t), Pbau (t0

hdis dis

(8)

Подсистема подсчета затрат микросети Модель микросети подсчитывает общий эксплуатационные расходы по формуле: Costmg (t) = Costel (t)• PMG (t) +

+ |CostdL-PMe (t)| + Costgas-P;TG (t), (9)

где: Costmg (t) - цена за энергоснабжение

микросети в момент времени t полученная с рынка электроэнергии на день вперед (DAM)

[Kc / (МВт • ч)] ; Costl = 1000 [Kc / (МВт • ч)] - цена на дистрибьюцию электричества принятая как константа; Costgas = 800 [Kc / (МВт-ч)] -полная цена за натуральный газ принятая как константа; PMG (t) - мощность, используемая микросетью из главной сети электроэнергии в

течение часа t [МВт-ч] ; PMG ( t ) - мощность, используемая микросетью из сети газопровода в течение часа t [МВт - ч].

2. УПРАВЛЕНИЕ МИКРОСЕТЬЮ НА ОСНОВЕ ОПТИМИТИЗАЦИИ ЗАТРАТ НА ЭНЕРГОПОТРЕБЛЕНИЕ

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

Был выбран экономический контроллер с прогнозирующей моделью (EMPC). Главное отличие EMPC от стандартного MPC, это то, что EMPC использует линейную весовую функцию, а не квадратную.

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

Таким образом, структуру EMPC можно представить из двух блоков (Рис. 2):

- участок генерации будущих горизонтов событий, на основе предсказаний;

- участок минимизации затрат, на основе матриц ограничений и представленных горизонтов.

Рассмотрим синтез каждого из участков ниже.

Блок минимизации затрат

Основа блока минимизации затрат — это функция минимизации. В его идее лежит итерационная модель оптимитизации, старающаяся минимизировать функцию затрат J каждый момент времени [8].

В данной работе все функциональные формулы были описаны в линейном виде, поэтому используется смешано-целочисленное линейное программирование (MILP). В Matlab это функция intlinprog, которую можно представить как:

minx JTx при условии

integer : x (intcon ) A • x < B

Aeq-x = Beq

lb < x < ub.

(10)

где: 3 - вектор, представляющий весовую функция на всю длину горизонтов, построенный при модификации формулы (9); х - решение М1ЬР на всю длину горизонтов, которое включает:

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

- значение заряда хранилищ на каждый час горизонтов;

- переключающая переменная для описания попеременного режима работы

(зарядки/разрядки) аккумулятора на каждый час горизонтов.

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

Рассмотрим более детально формулировку у^у

проблемы МГЬР. Критерий М1ЬР на всю длину горизонтов можно представить как:

Длина

Кол-во. упр.

j = £ гори— £ генераторов^Ц ,,. ■ PûWer, ,,., (11)

где: Costki - будущие цены, полученные из

горизонтов (ARIMA); Powerk ; - часть вектора х

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

Мы усложняем данный критерий строгими

требованиями ( Aeq ■ х = Beq ) для суммарной

генерации энергии для покрытия нагрузки на каждый будущий час горизонтов:

Кол-во. упр.

£ генераторовPower, i = Loadk , (12)

где: Loadk - общая предсказанная нагрузка микросети на каждый час горизонтов (ARIMA). Также с помощью строгих ограничений, мы устанавливаем начальные условия для Soc в данный момент времени:

SOC = SOC initial.

(13)

Рис. 3. Принцип работы выбранной ЛЯ1МЛ модели (20), где у (t) - временная серия, е( t) - ошибка предсказания.

Также вводится дополнительная

целочисленная переменная Ьк е[0;1], которая

является одной из частей вектора х и ограничивает одновременную зарядку/разрядку аккумуляторов. Ограничения описываются в виде:

0 < Ра^егкМ5 < Рапег^-(1 -Ьк), (15) 0 < Pawer.cn < РапегМН ■ К, (16)

С помощью не строгих ограничений ( A ■ x < B) устанавливается область оптимального

заряда/разряда аккумулятора на основе формул (6, 7):

Длина

Socm% < SOC0 + X Z~V SOCk < Soc9m0% , (14) где: □ SOCk - оптимальное изменение Soc на час k.

где: Power,

Power,

часть решения

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

горизонтов; Power,

функциональные аккумулятора.

-MAs, Powert -

лимиты разрядки/зарядки

500 к_L

-150 -100 -50 О

Бремя [ч]

Рис. 4. Профили нагрузки и цены за электричество за прошлую неделю и будущую неделю. Предсказания с помощью модели ARIMA представлены красным цветом

С помощью двухсторонних фиксированных ограничений (lb < x < ub) на решение в виде

вектора x, задаются эксплуатационные пределы искомых переменных на каждый час горизонтов: PowerMN < Power,J < Powerk^f , (17)

Soc™ < SOCk<SocZ°, (18)

0 < bk < 1. (19)

Выходом контроллера является решение MILP в виде вектора, с указаниями для блоков на текущий час (Powerli, SOC1,b1 ). Через час задача

оптимитизации MILP повторяется снова, при новых полученных горизонтах системы. Блок генерации горизонтов Для корректной работы блока минимизации затрат необходимо учитывать будущие значения системы, для того чтобы находить оптимальное решение оптимитизации. Например, учитывать будущие значения нагрузок микросетей и цен. Эти значения известны не точно, но могут быть предсказаны.

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

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

ARIMA ([1,24,168], 1, [1,24,168]) :

(1 -jL-jL24 -jL16S)(1 -L) yt =

= c + ( 1 + exL + q4 L24 + 0168 L168 ) e, (20)

где: L - лаговый оператор (L • xt = xt-1); yt -значение временной серии в момент времени t □ yt =(1 - L)1 yt); yt+1 - предсказание модели

(yt+1 = yt +□ yt); j, q, c - параметры модели (AR коэффициенты, MA коэффициенты, константа модели); et - ошибка предсказания в момент времени t.

Принцип работы ARIMA модели можно рассмотреть на Рис. 3. Для трансформации в стационарный ряд (во избежание влияния дрифта) используется предсказание разностей 1 порядка, а не абсолютных значений. В свою очередь расчет предсказанных разностей происходит с использованием сезонностей: рассматриваем вклад предыдущих разностей и предыдущих ошибок предсказаний час назад, день назад и неделю назад.

В связи со сложным расчетом структуры ARIMA модели, параметры были получены с использованием Econometric Toolbox в Matlab.

Генерация горизонтов в нашей модели работает, как показано на Рис. 4. Текущий час 0. Система смотрит на предыдущие значения разностей и ошибок предсказаний и шаг за шагом (или час за часом) конструирует новые горизонты с длиной в 7 дней. Для предсказания цен за электричество мы также используем данные с DAM.

320 330 34 0 350 360 370 330

Время [ч]

Рис. 5. Графики, отражающие поведения модели с EMPC с горизонтами в 1 день и оригинальными параметрами

0,5 1 1.5

Коэффициент параметра Рис. 6. Графики, отражающие поведения модели с EMPC при оригинальных

3. ИССЛЕДОВАНИЕ РАЗРАБОТАННОЙ УМНОЙ МИКРОСЕТИ

Проверка работоспособности модели умной микросети

Объединив модель микросети с распределенными генераторами и EMPC, включающий ARIMA модели и MILP, была получена модель умной микросети в Matlab.

Пример работоспособности данной модели умной микросети с оригинальными параметрами системы и горизонтами длиной 1 день можно увидеть на Рис. 5. Когда цена за электричество достигает максимума (например, в момент времени 330-го часа), система использует © AUTOMATICS & SOFTWARE ENGINERY. 2019

■0,02 0 0,02 Дельта параметра изменении параметров умной микросети относительно

энергию из аккумулятора и альтернативные малоэффективные источники энергии. Когда цена достигает минимума на горизонте (например, в момент времени 340-го часа), то энергия в аккумуляторах восполняется, и система отключает альтернативные малоэффективные источники энергии. Также контролер рассчитывает момент, когда стоит начать зарядку/разрядку аккумулятора.

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

Различные конфигурации умной микросети приводят к разным общим эксплуатационным затратам. В работе было изучено, как изменения таких параметров как длина горизонтов, номинальные мощности генераторов, емкость № 2 (28) 42

аккумулятора и КПД блоков влияют на общую экономию затрат £С [%], которая рассчитывается по формуле:

SC = ^ - Cost- .100%,

Cost,,

где:

Cost„

(21)

на

- суммарные затраты эксплуатацию микросети при условии, что вся микросеть питается только от возобновляемых источников энергии и главной энергосети;

CostE

суммарные затраты на эксплуатацию

микросети при использовании EMPC и всех указанных генераторов на Рис. 1.

Средняя экономия затрат умной микросети при стандартных параметрах и длинам горизонтов 3 дня равна 5.17%. Изменяя параметры умной микросети, можно увеличить экономию затрат до 21%.

ЗАКЛЮЧЕНИЕ

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

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

Для этой модели был синтезирован EMPC, который включает в себя:

- участок генерации будущих горизонтов, работающий с помощью ARIMA моделей;

- участок минимизации затрат микросети, работающий на основе смешано-целочисленного линейного программирования (MILP).

Работоспособность модели с EMPC была проверена.

Также был изучен вклад каждого из параметров системы в общую экономию затрат умной микросети. Общая экономия затрат достигает значения 5.17%, а экономия затрат измененной системы достигает значения 21%.

ЛИТЕРАТУРА

[1] Ayaz M.S., Azizipanah-Abarghooee R., Terzija V. European LV microgrid benchmark network: Development and frequency response analysis // 2018 IEEE International Energy Conference (ENERGYCON). 06.2018. С.1-6. D0I:10.1109/ENERGYC0N.2018.8398759.

[2] OTE. Statistics. Normalized Load Profile. 2019. URL: https://www.ote-cr.cz/en/statistics/electricity-loadprofiles/normalized-

lp?set language=en&date=201903-30; [Дата

обращения 9-04-2019].

[3] Carl C. Calculating Solar Photovoltaic Potential on Residential Rooftops in Kailua Kona, Hawaii: дис. маг. / Carl Caroline. Faculty of the USC graduate school. University of Southern California: University of Southern California, 05.2014.

[4] Johnson G.L. Wind Energy Systems. Prentice Hall, 1985. ISBN 0139577548.

[5] Capstone. Datashit C30 MicroTurbine. Natural Gas. 2019. URL: https://horizonpowersystems.com/wp-content/uploads/2016/03/C30-

NatGas 331031E lowres.pdf; [Дата обращения 1004-2019].

[6] Barone L., Kharitonov Z., Roth N. Solid Oxide Fuel Cells Powered by Natural Gas as The Solution for Clean Consumer Energy // University of Pittsburgh Swanson School of Engineering. 2016. Март. С.1-7.

[7] Rodriguez D.I.H., Myrzik J.M. Economic Model Predictive Control for Optimal Operation of Home Microgrid with Photovoltaic-Combined Heat and Power Storage Systems**This work was supported by the German Federal Ministry of Education and Research (BMBF) under the grant number 03EK3547 // IFAC-PapersOnLine. 2017. Т. 50, № 1. С. 10027 10032. ISSN2405-8963. DOI: https://doi.org/10.1016/j.ifacol.2017.08.2039.

[8] Halvgaard R. Model Predictive Control for Smart Energy Systems: дис. канд. / Halvgaard Rasmus. Kongens Lyngby: Technical University of Denmark. Department of Applied Mathematics, Computer Science, 2014.

[9] Asteriou D., Hall S.G. Applied econometrics. 2-е изд. Basingstoke: Palgrave Macmillan, 2011.

Павел Михайлович Ведель -

магистрант кафедры

автоматики.

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

автоматизации. E-mail: [email protected] 630073, Новосибирск, просп. К.Маркса, д. 20

Надежда Юрьевна

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

E-mail: [email protected] 630073, Новосибирск, просп. К. Маркса, д. 20

Статья получена 13.04.2019.

Modeling and Real Time Optimization of a Smart Microgrid of Power Supply of a District

with Distributed Energy Sources

P. M. Vedel, N.Y. Ovchinnikova

Abstract: In power supply systems, renewable energy generators play an increasingly important role. Different control methods are developed for such generators andt these methods take into account changes in the price of electricity. This paper presents a variant of energy management based on an economic model predictive control (Economic MPC), which minimizes the cost of microgrid operating with distributed power generation. The controller for coordinating and optimizing the microgrid operation in real time was developed based on a predictive economic model taking into account mixed-integer linear programming and it was then tested in Matlab / Simulink. The controller includes a block for predicting electricity consumption and price profiles of electricity based on ARIMA models. The change in the parameters of the smart microgrid was investigated for overall cost savings.

Keywords: Smart grid, Economic MPC, ARIMA models, Real time pricing.

REFERENCES

[1] AYAZ, M. S.; AZIZIPANAH-ABARGHOOEE, R.; TERZIJA, V. European LV microgrid benchmark network: Development and frequency response analysis. In: 2018 IEEE International Energy Conference (ENERGYCON). 2018, pp. 1-6. Available from DOI: 10.1109/ENERGYCON.2018.8398759.

[2] OTE. Statistics. Normalized Load Profile. 2019. URL: https://www.ote-cr.cz/en/statistics/electricity-loadprofiles/normalized-

lp?set language=en&date=201903-30; [Online;

accessed 9-04-2019].

[3] Carl, Caroline. Calculating Solar Photovoltaic Potential on Residential Rooftops in Kailua Kona, Hawaii. Faculty of the UCS graduate school. University of Southern California, 2014. Master's thesis. University of Southern California.

[4] Johnson G.L. Wind Energy Systems. Prentice Hall, 1985. ISBN 0139577548.

[5] Capstone. Datashit C30 MicroTurbine. Natural Gas. 2019. URL: https://horizonpowersystems.com/wp-content/uploads/2016/03/C30-

NatGas 331031E lowres.pdf; [Online; accessed 10-042019].

[6] BARONE, Luke; KHARITONOV, Zane; ROTH, Nicholas. Solid Oxide Fuel Cells Powered by Natural Gas as The Solution for Clean Consumer Energy. University of Pittsburgh Swanson School of Engineering. 2016, pp. 1-7.

[7] Diego I. Hidalgo Rodriguez, Johanna M.A. Myrzik, Economic Model Predictive Control for Optimal Operation of Home Microgrid with Photovoltaic-Combined Heat and Power Storage Systems**This work was supported by the German Federal Ministry of

Education and Research (BMBF) under the grant number 03EK3547, IFAC-PapersOnLine, Volume 50, Issue 1, 2017, Pages 10027-10032, ISSN 2405-8963, https://doi.Org/10.1016/j.ifacol.2017.08.2039.

[8] HALVGAARD, Rasmus. Model Predictive Control for Smart Energy Systems. Kongens Lyngby, 2014. PhD thesis. Technical University of Denmark. Department of Applied Mathematics and Computer Science.

[9] ASTERIOU, Dimitrios; HALL, Stephen G.

Applied econometrics. 2n Macmillan, 2011.

ed. Basingstoke: Palgrave

Pavel Vedel - master's student of Department of Automatic. Research interests: development, design, modeling and research of real automation objects. Email: [email protected] 630073, Novosibirsk, st. Prosp. K. Marksa, h. 20

Nadezhda Ovchinnikova -

master's student of Department of Automatic.

Research interests: development and design of automated process control systems. Email: [email protected] 630073, Novosibirsk, st. Prosp. K. Marksa, h. 20

The paper has been recieved on 13.04.2019.