CC BY
0УДК 519.873+621.311 DOI:10.25729/ESI.2023.31.3.007 Моделирование автономной микросети
Еделев Алексей Владимирович1, Карамов Дмитрий Николаевич2, Башарина Ольга Юрьевна3' 4
1Институт систем энергетики им. Л.А. Мелентьева СО РАН, 2Иркутский национальный исследовательский технический университет, 3Институт динамики систем и теории управления им. В.М. Матросова СО РАН, Россия, Иркутск, flower@isem. irk. ru,
4Уральский государственный экономический университет, Россия, Екатеринбург
Аннотация. Данная работа является первой из серии статьей, посвященных применению разработанной технологии цифрового двойника сложной технической системы для исследования ее свойств. Одним из таких свойств является живучесть, под которой понимается свойство системы адаптироваться к крупным возмущениям и восстанавливать свое исходное состояние после их воздействия. В качестве сложных технических систем в статье рассматриваются автономные микросети. Исследование живучести этих энергетических комплексов локального уровня обычно строится на многовариантных вычислительных экспериментах, однако при существовании обратной связи цифрового двойника с микросетью или испытательным стендом можно также использовать натурные опыты. Ключевым компонентом цифрового двойника микросети должен быть комплекс разнообразных моделей, с высокой точностью описывающих все аспекты поведения этого энергетического комплекса. В данной статье описывается оптимизационная модель микросети на основе концепции энергетического хаба. Энергетический хаб концентрирует производственные мощности по генерации, преобразованию и аккумулированию энергоресурсов, которые соединяются между собой и с потребителями с помощью сетей транспорта энергоресурсов. В статье представлена программная реализация оптимизационной модели микросети на основе концепции энергетического хаба. В заключении производится сравнение с принципиально близкой двухуровневой моделью энергетического комплекса города. Моделирование микросети на основе концепции энергетического хаба обеспечивает более высокую вычислительную эффективность. С другой стороны, модель энергетического комплекса города более точно, чем представленная в данной статье модель микросети, отражает физические особенности некоторых процессов, например транспорта тепла и электроэнергии.
Ключевые слова: цифровой двойник, испытательный стенд, микросеть, живучесть, математическая модель, энергетический хаб
Цитирование: Еделев А.В. Моделирование автономной микросети / А.В. Еделев, Д.Н. Карамов, О.Ю. Башарина // Информационные и математические технологии в науке и управлении. - 2023. - № 3(31). -С.74-85. - DOI: 10.25729/ESI.2023.31.3.007.
Введение. Автономные системы энергоснабжения являются разновидностью энергетических комплексов локального уровня (микросетей), в которых совместное использование возобновляемых источников энергии с традиционными энергетическими установками и накопителями электрической энергии является экономически эффективным способом энергоснабжения потребителей. Особую значимость и большое распространение возобновляемые источники энергии получили при электрификации труднодоступных территорий, не имеющих связи с традиционными топливно-энергетическими комплексами [1].
Изолированность автономных микросетей повышает значимость исследования их функционирования в следующих экстремальных условиях:
- природные бедствия, такие как ураганы, наводнения, землетрясения, резкие похолодания и т.д.;
- отказ элементов и подсистем вследствие износа, ошибок персонала и т.д.;
- преднамеренные (умышленные) нарушения, такие как кибератаки и т.д. Способность противостоять вышеперечисленным крупным возмущениям, не допуская
каскадного развития аварий с массовым нарушением режима энергоснабжения потребителей,
и восстанавливать исходное состояние или близкое к нему характеризует одно из свойств энергетических комплексов - живучесть [2].
Исследование живучести обычно строится на проведении крупномасштабных вычислительных экспериментов, в которых разыгрывается множество сценариев крупных возмущений. Для исследования живучести автономных микросетей вследствие их малого размера помимо вычислительных экспериментов также доступны натурные опыты на реальном оборудовании или на испытательных стендах, позволяющих моделировать функционирование микросетей в экстремальных условиях. Эти вопросы изучались в рамках научного проекта РФФИ и Правительства Иркутской области № 20-47-380002-р_а, который был посвящен разработке фундаментальных основ математического и информационного моделирования инфраструктурных объектов Байкальской природной территории, использующих природосберегающие технологии. В процессе выполнения проекта получены следующие основные методические результаты [3, 4]:
- разработана система мониторинга условий функционирования инфраструктурных объектов и их внешнего окружения, обеспечивающая сбор, хранение, обработку, агрегирование, структурирование и анализ ретроспективных и текущих предметных данных большого объема;
- созданы модели, методы и инструментальные средства спецификации и моделирования процессов функционирования исследуемых объектов на основе многокритериальных методов с использованием высокопроизводительных вычислительных систем в процессе проведения многовариантных расчетов;
- разработаны методы и инструментальные средства конструирования наборов цифровых двойников, представляющих различные подсистемы объекта, а также системы интеллектуальных агентов, которым делегируются права и обязанности субъектов, связанных с обеспечением эксплуатации объекта и потреблением предоставляемых им ресурсов и услуг;
- реализованы методы и средства микро-сервисного взаимодействия между цифровыми двойниками, а также двойниками и их конечными пользователями разных категорий. Цифровой двойник используется для исследования свойств моделируемого инфраструктурного объекта. Инфраструктурным объектом в данном статье выступает автономная микросеть, соответственно, предметной областью является её живучесть, а решаемым классом задач - анализ уязвимости этих энергетических комплексов. Традиционно исследование живучести строится на разыгрывании сценариев крупных возмущений с применением математических моделей вследствие утверждения, что проведение необходимых натурных экспериментов на реальных энергетических комплексах недопустимо [2]. В настоящее время испытательные стенды позволяют разворачивать полноценные реалистичные схемы микросетей и имитировать воздействие на них крупных возмущений. Таким образом, цифровой двойник микросети, соединенный с испытательным стендом, может для разыгрывания заданного множества крупных возмущений использовать не только вычислительные эксперименты, но и натурные опыты. Главным преимуществом цифрового двойника по сравнению математической моделью является наличие обратной связи с объектом исследования, что позволяет в процессе эксперимента подстраивать управление микросетью в экстремальных условиях, исходя из разницы между прогнозируемыми (полученными в результате расчета на математической модели) и реальными значениями технологических параметров оборудования.
Согласно [5], в современной энергетике наибольшее распространение получили три основных типа цифровых двойников:
1. Двойник-прототип - виртуальный аналог реального объекта, содержащий информацию, которая представляет объект в целом на всех стадиях его жизненного цикла.
2. Двойник-экземпляр - содержит информацию по описанию оборудования объекта, то есть данные о материалах, комплектующих, информацию от системы мониторинга оборудования.
3. Агрегированный двойник - объединяет два вышеперечисленных типа.
Основным требованием к цифровому двойнику микросети в [6] определена его ориентация на моделирование всевозможных воздействий (не только физических), прогнозирование последствий таких воздействий, выработка мер по предотвращению их негативного техногенного влияния. Ключевым компонентом цифрового двойника микросети должен быть комплекс оптимизационных, имитационных, нейросетевых моделей, с высокой точностью описывающих все аспекты поведения этого энергетического комплекса. Данная статья представляет оптимизационную модель, которая используется в двойнике-экземпляре для исследования живучести автономной микросети.
1. Математическая модель микросети. В целом в энергетических исследованиях применяется иерархический подход к моделированию, использующий принципы декомпозиции и агрегирования при формировании моделей энергетических систем на различных уровнях иерархии и взаимосвязей моделей разных уровней. Так, например, в модели энергетического комплекса современного города [7] выделяются два уровня иерархии: уровень централизованных систем энергоснабжения и уровень децентрализованно управляемых потребителей электро- и теплоэнергии. Формирование подобных моделей и расчёты на них в настоящее время могут осуществляться с помощью инструментария, представленного в [5, 8].
Сущность предлагаемого в данной работе подхода к моделированию автономных микросетей заключается в следующем: энергетический комплекс локального уровня состоит из произвольного числа взаимосвязанных компонентов (хабов) [9], в которых концентрируются производственные мощности по генерации, преобразованию и аккумулированию энергоресурсов [10-14]. Данный подход поддерживается специализированным инструментарием [15].
Концепция энергетического хаба позволяет детально учитывать факторы совместной работы накопителей электрической энергии в микросетях, влияющие как на техническую, так и экономическую эффективность использования тех или иных типов накопителей. В большей степени эти факторы касаются неравномерности генерации источниками, использующими энергию ветра и солнца в качестве первичных двигателей.
Накопители электрической энергии в микросетях, как правило, представлены разного рода аккумуляторными батареями. Накопители, в отличие от элементов генерации и трансформации электрической энергии, имеют определенный ряд специфических особенностей, которые необходимо учитывать при моделировании микросетей. Исходя из условий их эксплуатации [16], к данным особенностям можно отнести внутренние процессы старения активной массы аккумуляторных батарей, включая расслоение электролита, разрушение положительного электрода, тяжелую/необратимую сульфацию, кристаллизацию и т.п.
Стохастическая составляющая при генерировании электрической энергии существенно сказывается на основных эксплуатационных показателях аккумуляторных батарей. Например, в автономных фотоэлектрических системах эксплуатация аккумуляторных батарей сопровождается процессами с ярко выраженным циклическим характером в течение суток.
1.1. Общая схема энергетического хаба. Большинство схем хабов может быть построено на базе пяти типов элементарных блоков: входы, входные накопители энергии, преобразователи энергии, выходные накопители энергии и выходы [14] (рис. 1).
Входные накопители
Выходные накопители
Рис. 1. Общая схема энергетического хаба
Согласно данной схеме моделирование потоков энергии в хабе от его входов Р к его выходам Ь происходит в несколько этапов. Сначала необходимо описать, как входные потоки разбиваются на потоки Qвх-к входным накопителям с запасом энергии Евх и потоки F к преобразователям энергии. Преобразователи могут конвертировать один вид энергии в несколько различных. Затем преобразованные потоки суммируются и вместе с потоками выходных накопителей Qвblх с запасом энергии Евьа, образуют выходные потоки [14].
Входные и выходные потоки. Пусть рассматриваемый интервал времени состоит из Т периодов длительностью Ж, энергоресурсы представлены множеством К, а автономный энергетический комплекс представлен набором из взаимосвязанных хабов в количестве Н.
Каждый хаб Л = {1, ...,Н} имеет множество входов Ми, множество выходов Ыи, множество входных накопителей I™, множество выходных накопителей I. Каждый элемент из множеств Ми, , связан с одним определённым энергоресурсом к Е К. Также хаб Л
имеет набор преобразователей в количестве .
Как показано на рис. 2, в хабе Л = {1, ...,Н] в период времени t = {1,...,Т} поток энергоресурса Рьт на входе т Е Ми разделяется на поток Q{li входного накопителя энергии / Е I™ и потоки к преобразователям] = {1, ...,]ь}, связанным со входом т Е М^
РИт = ^¡Ев^т) ^Ит] + Qhi,
С Р ГИт] — гИт],
(1) (2)
где - пропускная способность связи преобразователя] со входом т хаба Л; &и(т) - множество, содержащее преобразователи хаба Л, связанные со входом т.
Р
кт
О
р
кт]
Рис. 2. Переток энергии в энергетическом хабе от входа к преобразователям
Образование потока энергоресурса на выходе п Е хаба Л = {1, ...,Н] показано на рис. 3. В период времени t = {1,...,Т} преобразователь энергии} = {1,...,]ъ} хаба Л конвертирует поток РИт] с коэффициентом преобразования и направляет его к выходу п. Принимая во внимание поток Qhi выходного накопителя I Е , поток Ь\1П на выходе п Е Ыи определяется как
Цгп = X
тЕМъ
Х]'Евп(т,п) Qhi,
(3)
где (ш, п) - множество преобразователей, посредством которых вход т Е Ми соединяется с выходом п хаба Л.
Рис. 3. Переток энергии в энергетическом хабе от преобразователей к выходу
Накопители энергии. В формулах (1) и (3) поток О,£ накопителя I Е I™ и I Ци11 хаба к = [1, ...,Н] формируется потоками заряда ОС^3'1 и разряда Qdd¡'t в период t = [1, ...,Т}:
^ = - ОН?'1- (4)
Запас энергии Е^ в накопителе I Е I™ и 1%и1 хаба к = [1, ...,Н} на конец периода времени t = [1,... ,Т} зависит от его емкости Ем, запаса ЕЫ-1 на конец предыдущего периода t — 1, потоков заряда и разряда Qddls't в период Ь:
Е*ы = ЕI-1 + (Лыв Оы3'11 — ^ Ь (5)
0<Е*Ы<_ Ем, (6)
Еы = Еы, (7)
где г/1^3 - коэффициент эффективности зарядки накопителя / хаба к; - коэффициент эффективности разрядки накопителя / хаба к. Уравнение (7) предполагает, что перед началом работы ^ = 0) все накопители содержат полный заряд.
Так как накопитель I Е I™ и I Ци1 хаба к = [1, ...,Н} в период t = [1, ...,Т} может находиться только в состоянии заряда или разряда, то вводится дополнительная бинарная переменная г,I Е [0,1} и уравнения, связанные с ней:
0 < О^1 < (1 — , (8)
0<QchЫ3'l<zíli(QchЫ3, (9)
где - максимально возможное значение потока разряда накопителя / хаба к; О^3 - максимально возможное значение потока заряда накопителя / хаба к. Если г,£ = 0, то согласно уравнению (8) накопитель / хаба к находится в состоянии разряда в период Если г,I = 1, то согласно уравнению (9) накопитель / хаба к находится в состоянии заряда в период
Уравнения (6), (8), (9) задают технические ограничения на работу накопителя / в период .
1.2. Сети транспорта энергоресурсов. Энергетические хабы соединяются между собой и с потребителями с помощью сетей транспорта энергоресурсов. Каждая сеть представляет собой граф Ск = (Ук,Ек), где Ук - множество узлов, Ек - множество ориентированных дуг, к - энергоресурс из множества К. На граф Ск могут накладываться дополнительные ограничения, связанные с физическими законами перемещения энергоресурса к.
В сеть йк может входить множество узлов V3 сУк, которые представляют обособленные источники энергоресурса к Е К, не относящиеся ни к одному из Н хабов. Например, в V3 могут быть выделены источники топлива.
Через Vк
с Ук обозначается множество узлов потребления энергоресурса к & К
Остальные узлы из множества Ук = Ук\(У3 и Ук ) называются промежуточными и служат местами пересечения потоков энергоресурса в сети транспорта. Каждый узел V сети ^ должен удовлетворять одному из условий баланса потоков энергоресурса к Е К в период времени =
[1.....Т}:
+ %^'У)ЕУк(у) — = Е Ук, (10)
^^'Р)ЕУк(р) акqvfкqv — ^(уд)ЕПк(у) + = Е Ук > (11)
^(q,v)EVk(v)akqvfk.qv ^(v,q)Enk(v) fkvq = 0, v Е Vk , (12)
0<utkv<ätkv,vEV£, (13)
0<fLq<fkvq,(v,q)EEk, (14)
0<gkv<gtkv,veVke, (15)
где Wki,v) — Ek - множество дуг, входящих в узел v; Ok(v) — Ek - множество дуг, исходящих из узла v; fkqv - поток энергоресурса к по дуге (q, v) Е Ek в период времени t; gkv - производство энергоресурса к в узле v EV^ в период времени t; dkv - потребность узла v Е Vk в энергоресурсе к в период времени t; utkv - недопоставка энергоресурса к в узел v EVk в период времени t; akqv - коэффициент, характеризующий потери при транспорте энергоресурса к по дуге(q, v) Е Ek; fkvq, (v, q) - пропускная способность дуги (q, v) Е Ek.
Если учитывать энергетические хабы, соединённые с Gk в узле v, то уравнения (10)-(12) для периода времени t = {1, ...,Т} запишутся как:
gkv + ^L(q,v)EWk(v) akqvfkqv + ^hEA3(v) LLhn — ^hEA^(v) Phm — ^(v,q)EDk(v) fkvq = 0,
v E V/, (16)
^(q,v)EVk(v) akqvfkqv + ^hEA3(v) Ltin — ^hEA.k(v)Phm — ^(v,q)EDk(v) fkvq + ukv = ^kv,
v E V£, ^ ^ (17)
,v)EWk(v)akqvflcqv + ^hEA3(v)^hn ^hEA%(v)Phm ^(v,q)Enk(v)fkvq = 0,
v E V£, (18)
где Ak(v) - множество хабов, входы которых соединены с сетью Gk в узле v; P^k - поток энергоресурса к Е К, поступающего на вход m Е Mh хаба h Е Ak(v) в период t; Ak(v) - множество хабов, выходы которых соединены с сетью Gk в узле v; Lthn - поток энергоресурса к Е К, поступающего на выход п Е Nh хаба h Е Ak(v) в период t.
1.3. Целевая функция. Целью оптимизации функционирования автономного энергетического комплекса является наиболее полное удовлетворение потребностей в энергоресурсах при минимальных затратах на эксплуатацию хабов и транспортных сетей на протяжении всех периодов времени t = {1,... ,Т}:
Et=1 ZmEMh Ch.kPhm + ^kEK ^vEvg ykvukv + ^kЕК ^(v,q)EEk ckvqfkvq\ ^ тп, (19)
где C^k - цена на энергоресурс к Е К, поступающего на вход m Е Mh хаба h = {1,..., Н} в период t; Phm - поток энергоресурса к Е К, поступающего на вход m Е Mh хаба h = {1, ...,Н} в период t; ykv - штраф за недопоставку энергоресурса к Е К в узел v EVk в период t; ckvq -удельная стоимость транспорта энергоресурса к Е К по дуге (v, q) Е Ek; fkkqv - поток энергоресурса к Е К по дуге (q, v) Е Ek в период t.
Ограничениями для целевой функции (19) выступают уравнения (1)-(9), (13)-(18).
2. Описание конфигурации микросети. Конструирование цифрового двойника реальной микросети [17, 18] невозможно без её спецификации или описания конфигурации этого энергетического комплекса. В целом, понятие конфигурации энергетического комплекса в литературе трактуется весьма широко. В рамках научного проекта РФФИ и Правительства Иркутской области № 20-47-380002-р_а конфигурацией был назван определённый вариант инфраструктуры энергетического комплекса, которая включает в себя перечень объектов энергетики, входящих в состав комплекса, а также список природно-климатических, социальных, экономических и прочих параметров энергетических технологий, применяемых на этих объектах. Например, для детального моделирования возобновляемых источников энергии требуется подробная природно-климатическая информация. К ней относятся температура воздуха,
атмосферное давление, влажность, плотность, солнечная радиация, скорость и направление ветра, а также некоторые другие метеорологические показатели [ 19, 20].
Математическая постановка (1)-(9), (13)-(19) предоставляет достаточный набор элементов для описания возможных конфигураций микросети:
- узлы сетей транспорта: источники и потребители энергоресурсов, промежуточные узлы;
- дуги сетей транспорта;
- преобразователи;
- накопители энергии (хранилища энергоресурсов).
Таким образом, описание конфигурации микросети в виде набора из Н взаимодействующих энергетических хабов состоит из двух частей. Первая часть описывает инфраструктуру передачи энергоресурсов к Е К в виде множества узлов Ук и дуг Ек транспортных сетей ^. Вторая часть описывает инфраструктуру генерации, преобразования и аккумулирования энергоресурсов в виде информации о входах, входных накопителях, преобразователях, выходных накопителях и выходах каждого из Н энергетических хабов. Данный способ спецификации энергетического комплекса достаточен для описания конфигураций автономных микросетей. Следует отметить, что в литературе представлены более сложные и универсальные способы спецификации микросетей на базе концепции энергетического хаба [10-13].
В рамках проекта № 20-47-380002-р_а были разработаны инфологическая и физическая модели данных для хранения конфигураций микросетей [21].
3. Программная реализация. После того, как описание конфигурации микросети считано из базы данных, оно трансформируется в вектора и матрицы согласно математической модели (1)-(9), (13)-(19). Этот шаг достаточно легко автоматизируется, так как каждому элементу конфигурации соответствует свой набор уравнений, описывающий функционирование этого элемента. Затем вектора и матрицы заполняются данными, характеризующими функционирование микросети на протяжении Т периодов времени. В конце происходит расчёт пото-кораспределения энергоресурсов К в микросети. Этот шаг состоит из этапов предварительной оптимизации и решения задачи (1)-(9), (13)-(19).
Вышеописанные шаги реализованы в виде программного модуля, который использует один из пакетов математической оптимизации: GLPK [22], С1р [23], СЬс [24], СРЬБХ [25] и 8оР1ех [26]. Пакеты CPLEX и СЬс ориентированы на параллельные и распределённые вычисления. Алгоритм работы реализованного программного модуля показан на рис. 4.
Заключение. В качестве выводов можно выделить следующие основные отличия математической постановки (1)-(9), (13)-(19) от аналогичной модели городского энергетического комплекса, описанной в [7]:
- более низкий уровень территориальной иерархии рассматриваемых энергетических комплексов и меньшая на три порядка величина нагрузок на стороне потребителя;
- использование концепции энергетического хаба, жестко заданная структура которого позволяет достаточно просто автоматизировать перевод текстового описания конфигурации энергетического комплекса в целевые структуры данных конкретного пакета математической оптимизации;
- возможность оптимизации показателей функционирования автономного энергетического комплекса по данным типичного метеорологического года или заданного многолетнего периода метеорологических наблюдений [16,17,27];
- ориентация на более детализированное моделирование режимов работы возобновляемых источников энергии и аккумуляторов электроэнергии [1].
Рис. 4. Алгоритм расчёта потокораспределения энергоресурсов в микросети
Особенностями программной реализации модели микросети (1)-(9), (13)-(19) по сравнению с подходом, представленным в [8], являются:
- применение готовых пакетов математической оптимизации, ориентированных на параллельные и распределённые вычисления;
- применение концепции параллельных и распределённых вычислений, что требует специализированной вычислительной среды для сборки, тестирования и выполнения разрабатываемого программного обеспечения.
Таким образом, сравнение с похожей двухуровневой моделью энергетического комплекса города [7] показывает, что благодаря тому, что все процессы производства, преобразования, аккумулирования и передачи энергоресурсов описываются линейными зависимостями, моделирование микросети на основе концепции энергетического хаба обеспечивает более высокую вычислительную эффективность. С другой стороны, модель энергетического комплекса города более точно и качественно, чем представленная в статье модель микросети, отражает физические особенности некоторых процессов, например транспорта тепла и электроэнергии. Дальнейшее совершенствование модели микросети связано с добавлением новых элементарных блоков, более детально описывающих работу электрохимических накопителей энергии и возобновляемых источников энергии, а также с формированием различных вариантов компоновки аккумулирующего звена на основе методики категоризации электрохимических накопителей энергии [1].
Благодарности. Исследование выполнено при поддержке Министерства науки и высшего образования Российской Федерации, проект № FNEW-2021-0005 «Технологии разработки и анализа предметно-ориентированных интеллектуальных систем группового управления в недетерминированных распределенных средах».
Список источников
1. Карамов Д.Н. Интеграция процесса категоризации электрохимических накопителей энергии в задачу оптимизации состава оборудования автономных энергетических комплексов, использующих возобновляемые источники энергии / Д.Н. Карамов // Известия Томского политехнического университета. Инжиниринг георесурсов, 2019. - Т. 330. - № 5. - С. 113-130.
2. Надежность систем энергетики и их оборудования. В 4-х томах. Т. 1. Справочник по общим моделям анализа и синтеза надежности систем энергетики / ред. Ю.Н. Руденко. - М.: Энергоатомиздат, 1994.- 480 с.
3. Bychkov I.V., Gorsky S.A., Edelev A.V. et al. Support for managing the resilience of energy systems based on a combinatorial approach. Journal of computer and systems sciences international, 2021, v. 60, no. 6. pp. 981-994, D0I:10.1134/S1064230721060071.
4. Kostromin R., Basharina O., Feoktistov A., Sidorov I. Microservice-based approach to simulating environmentally-friendly equipment of infrastructure objects taking into account meteorological data. Atmosphere, 2021, v. 12, no. 9, 1217, pp. 1-24, DOI: 10.3390/atmos12091217.
5. Hosseini S.H.R., Allahham A., Walker S.L., Taylor, P. Optimal planning and operation of multi-vector energy networks: A systematic review. Renewable and sustainable energy reviews, 2020, 133, p.110216.
6. Тягунов М.Г. Использование цифрового двойника гибридных энергетических комплексов для оптимизации их параметров и режимов / М.Г. Тягунов, Р.П. Шевердиев // Вестник МЭИ, 2023. - № 3. - С. 109-118. - D0I:10.24160/1993-6982-2023-3-109-118.
7. Воропай Н.И. Методика управления спросом на электро-и теплоэнергию в интегрированной энергосистеме c активными потребителями / Н.И. Воропай, В.А. Стенников, Е.А. Барахтенко и др. // Известия РАН. Энергетика, 2020. - Т. 4. - С. 11-23.
8. Barakhtenko E., Sokolov D. An architecture of the technology platform for computer modeling, design, and optimization of intelligent integrated energy systems. Proceedings of International multi-conference on industrial engineering and modern technologies (FarEastCon), 2019, pp. 1-5.
9. Eladl A.A., El-Afifi M.I., El-Saadawi M.M., Sedhom B.E. A review on energy hubs: Models, methods, classification, applications, and future trends. Alexandria engineering journal, 2023, vol. 68, pp. 315-342.
10. Liu T., Zhang D., Wang S., Wu T. Standardized modelling and economic optimization of multi-carrier energy systems considering energy storage and demand response. Energy conversion and management, 2019, vol. 182, pp.126-142.
11. Wang Y., Cheng J., Zhang N., Kang C. Automatic and linearized modeling of energy hub and its flexibility analysis. Appl energy feb., 2018, pp. 211:705-14.
12. Wang Y., Zhang N., Kang C., Kirschen D.S., Yang J., Xia Q. Standardized matrix modeling of multiple energy systems. IEEE Trans smart grid, 2019, vol. 10, pp. 257-70.
13. Huang W., Zhang N., Wang Y., Capuder T., Kuzle I., Kang C. Matrix modeling of energy hub with variable energy efficiencies. International journal of electrical power & energy systems, 2020, 119, p. 105876.
14. Almassalkhi M., Hiskens I. Optimization framework for the analysis of large-scale networks of energy hubs. Proceedings of 17th power system computation, 2011, p. 7.
15. Kriechbaum L., Scheiber G., Kienberger T. Grid-based multi-energy systems - modelling, assessment, open source modelling frameworks and challenges. Energy, sustainability and society, 2018, 8, p.35.
16. Карамов Д.Н. Формирование исходных метеорологических массивов с использованием многолетних рядов fm 12 Synop и metar в системных энергетических исследованиях / Д.Н. Карамов // Известия Томского политехнического университета. Инжиниринг георесурсов, 2018. -Т. 329. - № 1. - С. 69-88.
17. Kostromin R., Feoktistov A., Voskoboinikov M. Service-oriented tools for automating digital twin development. Proceedings of the 4th Scientific-practical workshop on information technologies: algorithms, models, systems (ITAMS 2021), CEUR-WS Proceedings, 2021, vol. 2984, pp. 95-100, DOI:10.47350/ITAMS.2021.12.
18. Kostromin R., Feoktistov A. Agent-based DevOps of software and hardware resources for digital twins of infra-structural objects. Proceedings of the 4th International conference on future networks and distributed systems (ICFNDS 2020), ACM, 2020, pp. 1-6, DOI:10.1145/3440749.
19. Карамов Д.Н. Математическое моделирование солнечной радиации с использованием многолетних метеорологических рядов, находящихся в открытом доступе / Д.Н. Карамов // Известия Томского политехнического университета. Инжиниринг георесурсов, 2017. - Т. 328. - № 6. - С. 28-37.
20. Karamov D.N., Naumov I.V. Modeling a solar power plant with regard to changes in environmental parameters. Proceedings of Power technology and engineering, 2020, 4, pp. 448-454.
21. Еделев А.В. База данных для моделирования автономных энергетических комплексов / А.В. Еделев, А.Г. Феоктистов // Техническая и технологическая модернизация России. Проблемы, приоритеты, перспективы: Тр. Международной научно-практической конф. - Казань: Аэтерна, 2021. - С. 3-9.
22. GLPK - GNU Linear Programming Kit. Available at: http://www.gnu.org/software/glpk/ (accessed: 05/04/2023).
23. Clp - Coin-Or Linear Programming. Available at: http://projects.coin-or.org/Clp (accessed: 05/04/2023).
24. Cbc - Coin-Or Branch and Cut. Available at: http://projects.coin-or.org/Cbc (accessed: 05/04/2023).
25. Cplex V12. 1: User's Manual for CPLEX. International Business Machines Corporation, 2009, 46, p.157.
26. SoPlex - The Sequential Object-Oriented Simplex. Available at: http://soplex.zib.de (accessed: 05/04/2023).
27. Svoboda V., Wenzl H., Kaiser R., Jossen A. et al. Operating conditions of batteries in off-grid renewable energy
systems. Solar energy, 2007, 81, pp. 1409-1425.
Еделев Алексей Владимирович. К.т.н., старший научный сотрудник Института систем энергетики им. Л.А. Мелентьева. Область научных интересов - живучесть энергетических комплексов, математическое моделирование, распределенные вычисления. AuthorlD: 710350, SPIN: 8447-9522, ORCID: 0000-0003-2219-9754, flower@isem.irk.ru.
Карамов Дмитрий Николаевич. К.т.н., доцент, Иркутский национальный исследовательский технический университет, Байкальский институт БРИКС. Область научных интересов - энергетические системы и комплексы, возобновляемые источники энергии, математическое моделирование. AuthorlD: 905773, SPIN: 99771303, ORCID: 0000-0001-5360-4826, dmitriy.karamov@mail.ru.
Башарина Ольга Юрьевна. К.т.н., доцент, научный сотрудник Института динамики систем и теории управления им. В.М. Матросова СО РАН, доцент Уральского государственного экономического университета. Область научных интересов - системный анализ и моделирование сложных систем. AuthorID: 702499, SPIN: 2612-2891, ORCID: 0000-0002-7151-782X, basharinaolga@mail.ru.
UDC 519.873+621.311 DOI:10.25729/ESI.2023.31.3.007
Modelling autonomous microgrids
Alexey V. Edelev1, Dmitriy N. Karamov2, Olga Yu. Basharina3' 4
1Melentiev Energy Systems Institute SB RAS,
2National Research Irkutsk State Technical University,
3Matrosov Institute for System Dynamics and Control Theory of SB RAS,
Russi a, Irkutsk, flower@isem.irk.ru,
4Ural State University of Economics, Russia, Yekaterinburg
Abstract. This work is the first of a series of articles devoted to the application of the developed digital twin technology of a complex technical system for the study of its properties. One of these properties is resilience, which is understood as the property of these energy complexes to adapt to large disturbances and restore their initial state after their impact. Autonomous micro-networks are considered as complex technical systems in the article. The study of the resilience of these local-level energy complexes is usually based on multivariate computational experiments, however, if there is a feedback of a digital twin with a micro-grid or a test bench, field experiments can also be used. The key component of the digital twin of the microgrid should be a complex of various models that accurately describe all aspects of the behavior of this energy complex. This article describes an optimization model of a microgrid based on the concept of an energy hub. The energy hub concentrates production capacities for generating, converting and accumulating energy resources, which are connected to each other and to consumers using energy transport networks. The article presents a software implementation of an optimization model of a microgrid based on the concept of an energy hub. In conclusion, a comparison is made with a fundamentally similar two-level model of the city's energy complex, showing that the modeling of a microgrid based on the concept of an energy hub provides higher computational efficiency. On the other hand, the model of the energy complex of the city more accurately than the microgrid model presented in this article reflects the physical features of some processes, for example, the transport of heat and electricity. Keywords: digital twin, test-bed, micro-grid, resilience, mathematical model, energy hub
Acknowledgements: The study was supported by the Ministry of Science and Higher Education of the Russian
Federation, project no. № FWEW-2021-0005 «Technologies for the development and analysis of subject-oriented
intelligent group control systems in non-deterministic distributed environments».
References
1. Karamov D.N. Integratsiya protsessa kategorizatsii elektrokhimicheskikh istochnikov energii v unikal'nom sostave oborudovaniya avtonomnykh energeticheskikh kompleksov, ispol'zuyushchikh istochniki energii [Integration of the storage battery categorization process into the task of optimizing the equipment of stand-alone energy systems with renewable energy sources]. Izvestiya Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov [Bulletin of the Tomsk Polytechnic University. Geo assets engineering], 2019, v. 330, no. 5, pp. 113-130.
2. Nadezhnost' sistem energetiki i ikh oborudovaniya. V 4-kh tomakh. T. 1. Spravochnik po obshchim modelyam analiza i sinteza nadezhnosti sistem energetiki [Reliability of energy systems and their equipment. In 4 volumes. T. 1. Handbook of general models of analysis and synthesis of reliability of energy systems], vol. 1., Ed Rudenko Yu.N., Moscow, Energoatomizdat, 1994, 480 p.
3. Bychkov I.V., Gorsky S.A., Edelev A.V., Kostromin R.O., Sidorov I.A., Feoktistov A.G., Fereferov E.S., Fedo-rov R.K. Support for managing the resilience of energy systems based on a combinatorial approach. Journal of computer and systems sciences international, 2021, v. 60, no. 6. pp. 981-994, DOI: 10.1134/S1064230721060071.
4. Kostromin R., Basharina O., Feoktistov A., Sidorov I. Microservice-based approach to simulating environmentally-friendly equipment of infrastructure objects taking into account meteorological data. Atmosphere, 2021, v. 12, no. 9, 1217, pp. 1-24, DOI: 10.3390/atmos12091217.
5. Hosseini S.H.R., Allahham A., Walker S.L., Taylor, P. Optimal planning and operation of multi-vector energy networks: A systematic review. Renewable and sustainable energy reviews, 2020, 133, p.110216.
6. Tyagunov M.G., Sheverdiyev R.P. Ispol'zovaniye dvoynika gibridnykh energeticheskikh kompleksov dlya opti-mizatsii ikh parametrov i rezhimov [Using a digital twin of hybrid energy complexes to optimize their parameters and modes]. Vestnik MEI [Bulletin of MPEI], no 3, pp. 109-118, DOI: 10.24160/1993-6982-2023-3-109-118.
7. Voropay N.I., Stennikov V.A., Barakhtenko E.A. et al. Metodika upravleniya sprosom na elektro- i teploenergiyu vo vstroyennoy energosisteme c aktivnymi pol'zovatelyami [Demand management technique for electricity and heat in an integrated energy system with active consumers]. Izvestiya RAN. Energetika [Izvestiya RAS. Energy], 2020, v. 4, pp. 11-23.
8. Barakhtenko E., Sokolov D. An architecture of the technology platform for computer modeling, design, and optimization of intelligent integrated energy systems. Proceedings of International multi-conference on industrial engineering and modern technologies (FarEastCon), 2019, pp. 1-5.
9. Eladl A.A., El-Afifi M.I., El-Saadawi M.M., Sedhom B.E. A review on energy hubs: Models, methods, classification, applications, and future trends. Alexandria engineering journal, 2023, vol. 68, pp. 315-342.
10. Liu T., Zhang D., Wang S., Wu T. Standardized modelling and economic optimization of multi-carrier energy systems considering energy storage and demand response. Energy conversion and management, 2019, vol. 182, pp.126-142.
11. Wang Y., Cheng J., Zhang N., Kang C. Automatic and linearized modeling of energy hub and its flexibility analysis. Appl energy feb., 2018, pp. 211:705-14.
12. Wang Y., Zhang N., Kang C., Kirschen D.S., Yang J., Xia Q. Standardized matrix modeling of multiple energy systems. IEEE Trans smart grid, 2019, 10, pp. 257-70.
13. Huang W., Zhang N., Wang Y., Capuder T., Kuzle I., Kang C. Matrix modeling of energy hub with variable energy efficiencies. International journal of electrical power & energy systems, 2020, vol. 119, p. 105876.
14. Almassalkhi M., Hiskens I. Optimization framework for the analysis of large-scale networks of energy hubs. Proceedings of 17th power system computation, 2011, p. 7.
15. Kriechbaum L., Scheiber G., Kienberger T. Grid-based multi-energy systems - modelling, assessment, open source modelling frameworks and challenges. Energy, sustainability and society, 2018, 8, p.35.
16. Karamov D.N. Formirovaniye iskhodnykh meteorologicheskikh massivov s ispol'zovaniyem mnogoletnikh rya-dov fm 12 Synop i metar v sistemnykh energeticheskikh proizvodstvakh [Formation of initial meteorological arrays using long-term series fm 12 Synop and metar in system energy studies]. Izvestiya Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov [Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering], 2018, vol. 329, no. 1, pp. 69-88.
17. Kostromin R., Feoktistov A., Voskoboinikov M. Service-oriented tools for automating digital twin development. Proceedings of the 4th Scientific-practical workshop on information technologies: algorithms, models, systems (ITAMS 2021), CEUR-WS Proceedings, 2021, v. 2984, pp. 95-100, DOI:10.47350/ITAMS.2021.12.
18. Kostromin R., Feoktistov A. Agent-based DevOps of software and hardware resources for digital twins of infra-structural objects. Proceedings of the 4th International conference on future networks and distributed systems (ICFNDS 2020), ACM, 2020, pp. 1-6, DOI:10.1145/3440749.
19. Karamov D.N. Matematicheskoye modelirovaniye sluchaynykh izlucheniy s ispol'zovaniyem mnogoletnikh me-teorologicheskikh ryadov, nakhodyashchikhsya v dostupe [Mathematical modeling of solar radiation using publicly available long-term meteorological series]. Izvestiya Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov [Bulletin of the Tomsk Polytechnic University. Geo Аssets Engineering], 2017, vol. 328, no. 6, pp. 28-37.
20. Karamov D.N., Naumov I.V. Modeling a solar power plant with regard to changes in environmental parameters. Proceedings of Power technology and engineering, 2020, 4, pp. 448-454.
21. Edelev A.V., Feoktistov A.G. Baza dannykh dlya modelirovaniya avtonomnykh energeticheskikh kompleksov [Database for modeling autonomous energy complexes]. Trudy mezhdunarodnoy nauchno-prakticheskoy konfer-entsii «Tekhnicheskaya i tekhnologicheskaya modernizatsiya Roscii. Problemy, prioritety, perspektivy» [Proc. of International scientific and practical conference "Technical and technological modernization of Russia. Problems, priorities, prospects"], Kazan', Aeterna,2021, pp. 3-9.
22. GLPK - GNU Linear Programming Kit. Available at: http://www.gnu.org/software/glpk/ (accessed: 05/04/2023).
23. Clp - Coin-Or Linear Programming. Available at: http://projects.coin-or.org/Clp (accessed: 05/04/2023).
24. Cbc - Coin-Or Branch and Cut. Available at: http://projects.coin-or.org/Cbc (accessed: 05/04/2023).
25. Cplex V12. 1: User's Manual for CPLEX. International Business Machines Corporation, 2009, 46, p.157.
26. SoPlex - The Sequential Object-Oriented Simplex. Available at: http://soplex.zib.de (accessed: 05/04/2023).
27. Svoboda V., Wenzl H., Kaiser R., Jossen A. et al. Operating conditions of batteries in off-grid renewable energy systems. Solar energy, 2007, 81, pp. 1409-1425.
Edelev Alexey Vladimirovich. Candidate of Technical Sciences, Researcher, Melentiev energy systems institute. The main direction of research - resilience of energy systems, mathematical modeling, distributed computing. AuthorID: 710350, SPIN: 8447-9522, ORCID: 0000-0003-2219-9754, flower@isem.irk.ru.
Karamov Dmitriy Nikolaevich. Candidate of Technical Sciences, Associate Professor, National research Irkutsk state technical university, Baikal School ofBRICS. The main direction of research - energy systems and power engineering complexes, renewable energy sources, mathematical modelling. AuthorID: 905773, SPIN: 9977-1303, ORCID: 00000001-5360-4826, dmitriy.karamov@mail.ru.
Basharina Olga Yurievna. Candidate of Technical Sciences, Associate Professor, researcher at the Matrosov institute for system dynamics and control theory of SB RAS, Associate Professor, Ural state university of economics. The main direction of research - system analysis and modeling of complex systems. AuthorID: 702499, SPIN: 2612-2891, ORCID: 0000-0002-7151-782X, basharinaolga@mail.ru.
Статья поступила в редакцию 30.08.2023; одобрена после рецензирования 12.09.2023; принята к публикации 25.09.2023.
The article was submitted 08/30/2023; approved after reviewing 09/12/2023; accepted for publication 09/25/2023.
