УДК 519.217
DOI: 10.38028/Е81.2022.28.4.005
Полумарковская и скрытая марковская модели функционирования автономного ветродизельного комплекса
Обжерин Юрий Евгеньевич, Сидоров Станислав Михайлович, Федоренко Сергей Николаевич
Севастопольский государственный университет, Россия, Севастополь, xaevec@mail.ru
Аннотация. Рассматривается задача оценки надежности и анализа функционирования автономного ветродизельного комплекса. Построены его полумарковская и скрытая марковская модели. Скрытая марковская модель используется для нахождения оценок характеристик функционирования автономного ветродизельного комплекса и прогнозирования его состояний на основе полученного в результате функционирования вектора сигналов. Результаты исследования получены в общем виде и инвариантны относительно законов распределения случайных величин, описывающих элементы, позволяют прогнозировать режимы работы автономного ветродизельного комплекса.
Ключевые слова: ветродизельный комплекс, полумарковская модель, скрытая марковская модель, анализ функционирования, прогнозирование
Цитирование: Обжерин Ю.Е. Полумарковская и скрытая марковская модели функционирования автономного ветродизельного комплекса / Обжерин Ю.Е., Сидоров С.М., Федоренко С.Н. // Информационные и математические технологии в науке и управлении. - 2022. - № 4(28). - С. 70-82. Б01: 10.38028/ЕБ1.2022.28.4.005.
Введение. В настоящее время возобновляемые источники электроэнергии повсеместно внедряются в системы электроснабжения. Совместное использование ветроэнергетических установок и дизельных генераторов способствует не только получению экономической выгоды, но и сокращению выбросов продуктов переработки топлива. [1, 2] Однако, эффективность работы ветроэнергетических установок зависит от многих факторов, в частности, от ветровой нагрузки, которая является стохастической. В связи с этим возникает задача построения адекватных математических моделей, учитывающих стохастический характер работы и используемых природных ресурсов (например, ветра) и наличие дизель-генератора (резервного источника энергии).
В статье [3] рассматривается приближенная модель ветродизельного комплекса в Ма1!аЬ^ти1тк. В статьях [4, 5] приводятся результаты исследований имитационного моделирования. В [6, 7] исследуется оценка режимов работы комплекса электроснабжения с вет-родизельной электростанцией. В [1, 2] используются марковские модели для анализа функционирования ветродизельного комплекса. Следует отметить, что большинство результатов, полученных другими авторами, рассматривают случаи, когда случайные величины, описывающие систему, имеют экспоненциальное распределение (марковская модель). Полумарковская модель, предлагаемая в данной статье, позволяет использовать распределения произвольного вида. Она также позволяет находить временные характеристики надежности и эффективности системы, в отличие от марковской модели. Во-вторых, для территориально удаленных систем (или систем, управляемых дистанционно) оператору не всегда удается при изменении их состояний полностью получить информацию, содержащуюся в кодировке состояний полумарковской модели, но всегда имеется возможность получить некоторый сигнал (информацию), связанный с состояниями вложенной цепи Маркова, которые можно считать ненаблюдаемыми (скрытыми). В таких случаях использование теории скрытых марковских моделей позволяет решать задачи нахождения и оценки характеристик на основе наблюдаемого вектора сигналов.
Цель нашей работы - применение теории полумарковских процессов с общим фазовым пространством состояний и теории скрытых марковских моделей для анализа надежности и эффективности автономного ветродизельного комплекса (ВДК).
В данной статье рассматривается задача оценки надежности и анализа функционирования автономного ВДК, состоящего из ветроэнергетической установки, рабочего и резервного дизельных генераторов, инвертора и аккумуляторной батареи.
В работе сначала строится полумарковская модель [8-12] автономного ВДК, позволяющая рассчитать стационарные и временные характеристики надежности. Затем, следуя предложенной авторами методике [18], на ее основе разрабатывается скрытая марковская модель [15-17], которая используется для решения задач прогнозирования и оценки ее характеристик с учетом заданных параметров.
1. Полумарковская модель автономного ВДК. Рассмотрим систему S, представляющую собой автономный ветродизельный комплекс. Будем считать инвертор и аккумуляторную батарею абсолютно надежными в том смысле, что вероятности их отказа намного меньше вероятностей отказа других элементов. Следовательно, не будем их учитывать при построении полумарковской модели. Тогда систему S можно представить трехэлементной системой с резервом времени: элемент 1 - ветроустановка, элемент 2 - дизельный генератор (ДГ), элемент 3 - резервный дизельный генератор (РДГ) (резерв времени). Под отказом элемента 1 следует понимать любое событие, приводящее к нарушению нормального функционирования (например, отсутствие ветра, поломка одной из составляющих и т.п.). Система функционирует следующим образом: если элемент 1 (2) выходит из строя, то элемент 2 (1) отключается, элемент 3 включается, и система функционирует за счет резерва времени (элемент 3). Как только элемент 1 (2) восстанавливается, то элемент 3 отключается (при этом считаем, что к следующему моменту времени его характеристики полностью восстанавливаются), а элемент 2 (1) приступает к работе с прежним уровнем наработки. Отказ системы наступает в том случае, когда элемент 1 (2) находится на восстановлении, элемент 2 (1) отключен, и заканчивается резерв времени (отказ элемента 3). Резерв времени заканчивается в случае, если либо кончается топливо резервного генератора (время работы за счет запаса топлива), либо он выходит из строя.
Предположим, что времена безотказной работы элементов 1 и 2 описываются случайными величинами (СВ) а1 и а2, соответственно, имеющими функции распределения Р^х) = Р(а 1 < х), Р2(х) = Р(а2 < х) и плотности распределения /1(х), /2(х), а время восстановления - случайными величинами и имеющими функции распределения С1(х) = Р(^1 < х), С2(х) = Р(^2 < х) и плотности распределения ^1(х), ^2(х),. Рассматривается случайный мгновенно пополняемый резерв времени, который описывается СВ т , имеющей функцию распределения Р(х) = Р(т < х) и плотности распределения г(х). Случайные величины а1, а2, ^1, т предполагаются независимыми, имеющими конечные математические ожидания. Отметим, что СВ т представляет собой минимум из времени работы за счет запаса топлива резервного генератора и времени его безотказной работы.
На рис. 1 представлена структурная схема рассматриваемой системы.
Рис. 1. Структурная схема рассматриваемой система Рассмотрим дискретно-непрерывное фазовое пространство состояний рассматриваемой системы:
элемент 1
элемент 2
РДГ
элемент 3
Е = (1и2х2, 2112х1, 1021х2, 2201х1, 3020х1х2, 3200х1х2}..
Фазовые состояния обозначаются согласно следующей схеме: первая цифра кода означает номер элемента, в котором произошло изменение состояний; остальные цифры - что происходит с каждым элементом (1 - работа, 2 - отключен, 0 - восстановление). Непрерывные компоненты показывают время, оставшееся до следующей смены состояний. Расшифруем значения кодов состояний:
- 1112х2 (2112х1) - элемент 1(2) восстановился, элемент 2(1) включается, элемент 3 отключен, х2 (х-^ - время до конца работы элемента 2(1).
- 1021х2 (2201х1) - элемент 1(2) выходит из строя и начинает восстановление, элемент 2(1) отключается, включается элемент 3, х2 (х-^ - время до конца работы элемента 2(1).
- 3020х1х2 (3200х1х2) - отказ системы: х1 (х2) - время до конца восстановления элемента 1(2), х2 (хх) - время до конца работы элемента 2(1).
Запишем вероятности переходов между состояниями системы:
^пт? = /1(х2 + > Р^ = Л(*2 - У2)<0 < У2 < *2; Р211°121хХ12 = /2(х1 + х2)< *2 > ^т1;? = /2(х1 - УД 0 < У1 <
то
>3°2°^2=! гф^ + ОЛ; °
1021Х
п œ
pïoïï? = I ^l(t)/?(t)rft = Р(т > ^l); pi3
0 то
^¿J = I 02(t)ä(t)dt = Р(т >
10
р
3200х1х2
= /"KO^fe + Orft; ^MO^ = 1'
р2112Х1 _ г3200х1х2 1
22°1х1
Выпишем систему уравнений для нахождения стационарного распределения ВЦМ. г р(1112х2) = р(1021х2)/0Ю#1(ОР(О^ + /отр(3020х1х2)^х1, р(2112х1) = р(2201х1) /0т#2(0Р(0^ + /отр(3200х1х2)^х2, р(1021х2) = /Х2/1(у х2) Р(1112У2)^У2 + /° /2(х2 + Х1)р(2112Х1)^Х1 , Р(2201Х1) = / /2(у
(2)
(3)
р(3020х1х2) = р(1021х2) г(0,д1(х1 + р(3200х1х2) = р(2201х1) г(0^2(ж2 + 0^,
/ р(е)^е = 1 (условие нормировки). Подставим уравнение 5 в уравнение 1 системы (3). Получим
- ТО г> ТО г> то
р(1112х2) = р(1021х2) I ^(0^(0^ + + р(1021х2) I ^х1 | г(0^1(ж1 + =
г> ТО г> ТО г> то
= р(1021х2) I Д[1(0Р?(0^+ I ^х1 | г(0^1(ж1 + =
г> ТО г> ТО г> то
= р(1021х2) I Д[1(0Р?(0^+ I г(0^ I ^1(х1 + 0^х1 =
г> ТО г> то
= р(1021х2) I I г(0<^(0^ = р(1021х2).
°°
Аналогично, подставляя уравнение 6 в уравнение 2 системы (3), получаем:
р(2112х1) = р(2201х1). Следовательно, система (3) приводится к виду:
Гр(1021Х2) = /;/1(У2 - Х2)Р(1021У2)^У2 + /°°°/2(Х2 + Х^р^О^^ ,
(р(2201х1) = /™/2(У1 - Х1)Р(2201У1)^У1 + /°°°/1(Х1 + Х2)Р(1021Х2)^Х2 . Решение системы (4) найдено в [12] и имеет вид:
р(И12Х2) = р(1021Х2) = Р°^(Х2),
(4)
р(2112х1) = р(2201х1) = р0Р1(х1),
г. ж
р(3020х^) = Ро^гМ I +
0
Р(3200Х1Х2) = Р0Р1М + (5)
константа р0 находится из условия нормировки. Времена пребывания в состояниях системы:
$1112х2 = а1 Лх2, ^2112%! = а2 Л х1, ^1021х2 = Л т,
^2201%! = ^2 Л т, ^3020хгх2 = х1, ^3200хгх2 = х2, (6)
где Л — знак минимума.
Тогда средние времена пребывания в состояниях равны:
ГХ2 _ ГХ1_
м0ц12х2 = I М02112Х, = I Р^ж
00
М0Ю21Х2 = М022О1Х! = /^ёКОВД^. (7)
^^3020хгх2 = х1, ^^3200хгх2 = х2.
Разобьем пространство состояний Е на два непересекающихся подмножества: работоспособных (£+) и отказовых состояний (£-):
£+ = (1112х2, 2112х1, 1021х2, 2201х1}, £- = (3020х1х2, 3200х1х2} Найдем стационарные характеристики надежности системы.
Используя вероятности переходов (2), стационарное распределение (5), средние времена пребывания в состояниях (7), получаем:
Ж рЖ
+
^Г Г Ж п Ж
ш(х)р(^х) = Ро I /22(х2)^х2 I
£+ -'О -'О
р Ж рЖ рЖ р
+Ро I /Г1(х1)^х1 I ё2(0й(0^ + ро I /22(х2)^х2 I +
0 0 0 0 /ч Ж /ч
+р0 I ^1(х1)^х1 I Р2(0^ = р0Ма2М(^1 Л т) + р0Ма1М(^2 Л т) + 00
+р0Ма1Ма2.
/ч Ж /ч Ж /ч Ж
Р(х,£-)р(^х) = р0 I ^х1 I Р^(х2)^х2 I г^)^^ + +
£+ -'о -'о -'о
Ж Ж Ж
+Ро I ^х2 I р!(х1)^х1 I (х2 + =
7о ./о -'о
ЖЖ
= роМа2 I г(0С2(0^ + роМа1 I г(0£22(0^ =
оо
= роМ«2Р(^1 > т) + РоМ«1Р(^2 > т).
« Ж Ж Ж
ш(х)р(^х) = ро I х1^х1 I Р:2(х2)^х2 I г^)^^ + +
-'о -'о -'о
Ж Ж Ж
+Ро I х2^х2 I р1(х1)^х1 I г(0#2(х2 + =
о о о
Ж Ж Ж Ж
роМа2 I х1^х1 I г^)^^ + + р^^ I х2^х2 I г(0$2(х2 +
о о о о
юльзуя результаты, полученные в [18, с. 61], получаем
I ш(х)р(^х) = РоМ«2[М^1 — М(^1 Лт)] +РоМ«1[М^2 — Лт)].
^ш(х)р(^х) = I ш(х)р(^х) + I ш(х)р(^х) =
Р ¿К Зр.,
— рсм«2[М& - Л т)] + РсМ«1 [М& — Л т)] + Лт) +
+р0Ма1М(^2 Лт) + р0Ма1Ма2 = р0[Ма1Ма2 + Ма2М^1 + Ма1М^2], где выше обозначено:
ш(х) - средние времена пребывания в состояниях, Р(х, £—) - вероятности перехода из работоспособных состояний в отказовые, Л т) - математическое ожидание минимума двух СВ, Р(^ > т) - математическое ожидание того, что СВ ^ больше СВ т.
«то - то
Р(^>Т)=| г(0ё(0^ , М(^Лт)— | ё(0й(0^ . 00 Используя формулы, представленные в [12, 13], и выражения, найденные выше, получаем:
- среднее стационарное время безотказной работы 7+:
/ ш(х)р(^х) Ма2М(& Л т) + Ма1М(^2 Л т) + Ма1Ма2
Р+ —
J P(x,P-)p(dx) Mo^CÄ > т)+ M«iP(^2 > т)
- среднее стационарное время восстановления 7- :
J£_ rn(x)p(dx) - м(^1 Л т)] + Ма1 [Mß2 - M(ß2 Л т)]
Т_ —
/ Р(х,Р-)р(^х) Ма^ОД > т)+ Ма^ > т)
- стационарный коэффициент готовности
/£+ ш(х)р(^х) Ма2М(& Л т) + Ма1М(^2 Л т) + Ма1Ма2
— - — -.
г 4 ш(х)р(^х) М«1М«2 + М«2М& + М«^
В случае неслучайного резерва времени, т.е. когда Р(х) — 1(х — Л), где И - величина резерва времени, стационарные характеристики надежности примут вид:
- среднее стационарное время безотказной работы 7+:
+ = Ма^^+Ма-^ф) ' (8)
- среднее стационарное время восстановления 7-:
- = Ма2ё1;(Л)+Ма1С2(й) ' ( )
- стационарный коэффициент готовности
—-^-::-■ (10)
В качестве иллюстративного примера использования формул (8)-(9), рассмотрим систему для которой перед началом её функционирования принято, что СВ а1, а2, имеют распределение Эрланга IV порядка и Ма1 — 150 ч., Ма2 — 100 ч., — 24 ч., М^2 — 20 ч., а резерв времени И изменяется от 10 до 20 часов. Полученные результаты изображены на рис. 2 (а, б, в).
в)
Рис. 2. Графики стационарных характеристик надежности при различных значениях резерва времени И: а - среднее время безотказной работы, б - среднее время восстановления, в - стационарный коэффициент готовности 2. Нахождение вероятностей перехода укрупненной полумарковской модели. Для построения скрытой марковской модели (СММ) укрупним по непрерывным компонентам стационарные распределения (5), т.е. применим алгоритм стационарного фазового укрупнения к каждому состоянию. Отметим, что фазовое пространство полумарковских состояний получается добавлением к кодам физических состояний совокупности непрерывных компонент, фиксирующих остаточные времена действия факторов, изменяющих состояния системы. [13] Эти непрерывные компоненты и желательно укрупнять, оставляя только дискретное множество физических состояний. Получим
р ж ^ж ^ж
I р(1112Х2)^Х2 = I р(1021X2)^X2 = р° I Р2(Х2)^Х2 = р°М«2, ° ° °
пЖ пЖ пЖ
I р(2112х1)^х1 = I р(2201х1)^х1 = р° I Р1(х1)^х1 = р°Ма1, ° ° °
Ж Ж Ж Ж Ж
I I р(3020х1х2)^х1 ^х2 = р° I ^(х2)^х2 I х1^х1 I г^)^^ + = ° ° ° ° °
Ж
= р°М«2 I = р°М«2Р(^1 > т).
°
Ж Ж Ж Ж Ж
I I р(3200х1х2)^х1 ^х2 = р° I Р1(х1)^х1 I х2^х2 I г(0,д2(х2 + = ° ° ° ° °
Ж
= р°М«1 I г(0С2(0<^ = р°М«1Р(^2 > т). °
Следовательно, дискретное фазовое пространство укрупненных состояний рассматриваемой системы имеет вид:
Ё = {1112, 2112, 1021, 2201, 3020, 3200}. Используя формулы, представленные в [12, стр. 36],
— /£ р(^е)Р(е,£'г) / р02е)т(е)
=-—-, тк = -——■ (11)
найдем вероятности переходов р] между укрупненными состояниями и средние времена пребывания т] в них.
Вычислим знаменатели формул (11).
р(Ё1112) = р(Ё1°21) = р°М«2, р(Ё2112) = р(Ё22°1) = р°М«1, р(Ёз°2°) = р°М«2Р(01 > Т), р(Ёз2°°) = р°М«1Р(^2 > т). Средние времена пребывания тк равны:
т1112 = ~ м^^ , т2112 = ,
гп1°21 = --- = Л т), ГП22°1 = М(^2 Л т),
_М01-М(01Лт) _М02-М(02Лт)
?тз°2° = Р^ >^) , ?тз2°° = Р(& > т) ■ Найдем вероятности переходов р]:
л1021 _ /°Ж^22(Х2)^Х2 /°2 /1(Х2 - № _ м^2 - м(«1 Л «2)
р1°12 = м«2 = м«2 ,
л22°^ м(«1 Л «2) л1°21 _ м(«1 Л «2) л22°^ м«1 - м(«1 Л «2)
р1112= м«2 , р2112= м^ ,р2112= ма1 ,
РИИ = Р(Т> &), р!3°°221° = Р(Т < &), р|111 = Р(Т > &),
р22°1 = Р(т < ^2), рlql1 = 1, р2112 = 1.
Тогда матрица переходных вероятностей Р^ между укрупненными состояниями системы равна:
Ма9 -М(а л а) М(а л а) 0 0 -2-^-^ —^-^ 0 0
Ма2 Ма2
М (а л а) Ма, - М (ала) 0 0 —^-И -1-^-И 0 0
PJ =
Ma Ma
Р(г>Р) 0 0 0 РрР) 0
0 Р(*>р2) 0 0 0 РрР2)
1 0 0 0 0 0 0 1 0 0 0 0
3. Скрытая марковская модель рассматриваемой системы. Для полного описания СММ [16, 17] необходимо определить:
1. Множество состояний модели соответствует множеству состояний укрупненной модели Ё.
2. Алфавит наблюдаемой последовательности (множество сигналов). Предположим, что при функционировании системы £ состояния ВЦМ укрупненной модели не наблюдаются (скрытые состояния), а наблюдаются только количество работоспособных элементов во время смены состояний ВЦМ. Введём следующее множество сигналов:
/ = {0,1,2},
где
- 0 - отказ системы;
- 1 - элементы 1 и 2 системы функционируют;
- 2 - система работоспособна за счет резервного дизельного генератора.
Множество сигналов можно выбирать по-разному. Множество сигналов J выбрано в таком виде, т.к. «точную» информацию о количестве работоспособных основных и резервных элементов можно получить практически для любой системы.
3. Матрицу переходных вероятностей между состояниями системы.
Для нашей модели матрица переходных вероятностей состоит из переходных вероятностей Р^ укрупненной полумарковской модели.
4. Связь состояний модели с сигналами.
Рассмотрим связь между состояниями вложенной цепи Маркова (ВЦМ) укрупненной модели и сигналами, т.е. определим функцию связи Р(5|ж) [16, 18]:
Вдж) = Р(5П = 5|ХП = е £,5 е/,у Вдж) = 1, где 5П - п-ый сигнал.
Будем считать, что сигналы 0 и 2 испускаются верно с вероятностью 0,99 и с вероятностью 0,01 ошибочно.
Функция Р(5|ж) связи состояний ВЦМ укрупненной модели с сигналами представлена в таблице 1.
Таблица 1. Функция связи Р(5|х) состояний ВЦМ укрупненной модели с сигналами
' ' ■——^^^^ Сигнал, s Состояние,х ——____ Б=0 Б=1 б=2
1112 0 1 0
2112 0 1 0
1021 0,01 0 0,99
2201 0,01 0 0,99
3020 0,99 0 0,01
3200 0,99 0 0,01
5. Начальное распределение вероятностей модели.
Будем считать, что в начальный момент времени укрупненная модель равновероятно может находится в состоянии 1112 или 2112.
СММ на основе укрупненной полумарковской модели построена.
4. Оценка характеристик и прогнозирование состояний. В качестве иллюстративного примера, рассмотрим систему £, для которой перед началом её функционирования принято, что СВ «2,^1,^2 имеют распределение Эрланга IV порядка и = 150 ч., Ма2 = 100 ч., М^! = 24 ч., М^2 = 20 ч., Мт = 15 ч. Предполагается, что в начальный момент времени система находится в состояниях 1112 и 2112 с вероятностями 0,5.
Предположим, что в результате функционирования системы £ получен следующий вектор сигналов:
5зо = (1,2,1,2,1,2,1,2,0,1,2,1,2,1,2,0,1,2,1,2,0,1,2,1,2,1,2,1,2,1), п=30.
Рассмотрим задачи по оценке характеристик скрытой марковской модели с учетом введённых параметров, следуя [18, 19].
1. Определим вероятности состояний скрытой модели в момент испускания 30-го сигнала (последнего сигнала в цепочке).
Т.к. 30-й сигнал равен 1 (элементы 1 и 2 системы функционируют), то стоит задача определения вероятности нахождения укрупненной системы в состояниях, приводящих к испуска-
нию указанного сигнала. В результате получаем, что на 30-ом шаге укрупнённая модель с вероятностью 0,3553 находилась в состоянии 1112, с вероятностью 0,6447 в состоянии 2112. Для остальных состояний эта вероятность равна нулю.
2. Найдём вероятности, с которыми скрытая модель осуществит переход в состояния на следующем 31-ом шаге.
Вероятности перехода скрытой модели на 31-ом шаге: в состояние 1021 с вероятностью 0,419, в состояние 2201 - 0,581, во все остальные - с нулевой.
3. Определим вероятности появления сигналов на следующем 31-ом шаге.
Получаем, что вероятность появления сигнала 2 на 31-м шаге равна 0,99, сигнала 0 -
0,01, сигнала 1 - 0.
4. Найдем вероятность появления (испускания) полученного вектора сигналов ?30.
Вероятность появления полученного вектора сигналов ?30 равна 0,00000243.
5. Прогнозирование состояний скрытой модели по полученному вектору сигналов.
В таблице 2 указаны наиболее вероятные состояния скрытой модели на указанных в ней переходах и вероятности этих состояний.
Применяя алгоритм алгоритм Баума-Велша [15, 16], получаем переоцененную матрицу переходных вероятностей для рассматриваемой системы. Исходная матрица переходных вероятностей и матрица переходных вероятностей после переоценки представлены на рисунке 3.
Таблица 2. Наиболее вероятные состояния скрытой модели на переходах
Номер перехода 1 7 11 17 21 2б 29
Наиболее вероятное состояние 1112 2112 2201 2112 3200 2112 2201
Вероятность состояния 0,5233 0,б737 0,б520 0,5190 0,513б 0,б132 0,б447
Применяя алгоритм Витерби [15, 16] к переоцененной модели, определяем наиболее вероятную цепочку состояний для полученного вектора сигналов: 1112, 2201, 2112, 1021, 1112, 2201, 2112, 1021, 3020, 1112, 2201, 2112, 1021, 2112, 2201, 3200, 2112, 2201, 2112, 1021, 3020, 1112, 2201, 2112, 1021, 2112, 2201, 2112, 1021, 2112.
о о
0,2606
0
1 о
о о о
о,3469
0
1
0,1454 0,8546 0,5697 0,4303
о о
о,7394
о о о
о о о
о,6531
о о
о о
0,7075
0
1 о
о о о
0,8066
0
1
о,1283 о,8717 о,5369 о,4631
о о
о,2925
о о о
о о о
о,1934
о о
А
/
а)
б)
Рис. 3. Матрицы переходных вероятностей: а) исходная матрица переходных вероятностей Р/, б) переоцененная матрица переходных вероятностей РР
Обсуждение. Результаты данной статьи открывают большие возможности для оценки функционирования автономных ветродизельных комплексов. В статье показано, как применить аппарат теории скрытых марковских моделей (используемый для систем с конечным множеством состояний) для полумарковских моделей, у которых общее фазовое пространство состояний, что расширяет возможности применения теории СММ.
В статье приведены лишь иллюстративные примеры, показывающие расширение возможностей моделирования автономных ветродизельных комплексов при использовании
разработанной авторами методики. Однако, очевидно, что для специалистов-энергетиков не составит труда использовать параметры реальных ветродизельных комплексов в качестве входных данных моделей с последующим расчетом характеристик.
Следует отметить, что полученные результаты при должной переформулировке могут быть использованы для анализа функционирования технических систем различного назначения.
Рис. 4. Решетчатая диаграмма функционирования системы Решетчатая диаграмма для первых 12 полученных сигналов изображена на рис. 4. Для простоты используется следующая перекодировка состояний модели: 1112 ^ 1, 2112 ^ 2, 1021 ^ 3, 2201 ^ 4, 3020 ^ 5, 3200 ^ 6. Наиболее вероятные состояния (по алгоритму Ви-терби для переоцененной модели) для вектора полученных сигналов выделены жирной линией. Тонкими линиями соединены другие возможные состояния, которые система могла бы достичь без информации о векторе сигналов.
Заключение. В данной работе построена полумарковская модель автономного ВДК, позволяющая рассчитать стационарные и временные характеристики надежности. Затем на ее основе разработана скрытая марковская модель, которая используется для нахождения оценок характеристик автономного ВДК и прогнозирования его состояний на основе полученного в результате функционирования вектора сигналов. Результаты исследования позволят прогнозировать режимы работы автономного ВДК. Полученная СММ позволяет переоценивать параметры построенной модели (обучать) для повышения ее адекватности, согласно получаемым в процессе функционирования сигналам.
Результаты исследования получены в общем виде и инвариантны относительно законов распределения случайных величин, описывающих элементы автономного ВДК. Они позволяют моделировать функционирование системы при различных законах распределения, основываясь на статистических данных, не видоизменяя саму модель. В дальнейшем планируется обобщить результаты на случай парка (состоящего из N штук) ветродизельных комплексов.
Благодарности. Исследование выполнено при поддержке гранта Президента Российской Федерации для государственной поддержки молодых российских ученых - кандидатов наук № МК-329.2022.4.
Список источников
1. Ветродизельные комплексы в децентрализованном электроснабжении: монография / А. В. Бобров, В.А. Тремясов // Красноярск: Изд-во СФУ, 2012. - 216 с.
2. Григорьева О.А. Анализ надежности автономного ветродизельного комплекса / О.А. Григорьева, Т.В Кри-венко, В.А. Тремясов // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета, 2016. - № 2(243). - С. 45-52.
3. Salima Meziane, Omar Feddaoui, Riad Toufouti Modeling and simulation of hybrid wind-diesel power generation system. International journal of renewable energy, 2013, vol. 8, no. 2, pp. 49-58.
4. Sebastián R., Nevado A. Study and simulation of a wind hydro isolated microgrid. Energies, 2020, vol. 13(22), DOI: 10.3390/en13225937.
5. Sebastián R. Modeling, simulation and control of wind diesel power systems. Energies, 2022, vol. 15, no. 1712, DOI: 10.3390/en15051712.
6. Shklyarskiy Y.E., Batueva D.E. Operation mode selection algorithm development of a wind-diesel power plant supply complex. Journal of mining institute, 2022, vol. 253, pp. 115-126, DOI: 10.31897/PMI.2022.7.
7. Baikov A.I., Dar'enkov A.B., Sosnina E.N. Simulation modeling of a wind-diesel power plant. Russian electrical engineering, 2018, vol. 89, pp. 161-167, DOI: 10.3103/S1068371218030033.
8. Limnios N., Nikulin M. Recent advances in reliability theory: methodology, practice, and inference. New York, Springer science+business media, 2000, 514 p.
9. Руденко Ю.Н. Надёжность систем энергетики / Ю.Н. Руденко, И.А. Ушаков // Новосибирск: Наука, 1989. - 328 с.
10. Обжерин Ю.Е. Полумарковские и скрытые марковские и полумарковские модели систем энергетики / Ю.Е. Обжерин // Известия РАН. Энергетика, 2019. - № 5. - С. 26-32.
11. Обжерин Ю.Е. Анализ надежности и фазовое укрупнение систем с поэлементными накопителями / Ю.Е. Обжерин, С.М. Сидоров, М.М. Никитин // Известия РАН. Энергетика, 2019. - № 6. - С. 66-77.
12. Королюк В.С. Процессы марковского восстановления в задачах надёжности систем / В.С. Королюк, А.Ф. Турбин // Киев: Наук. думка, 1982. - 236 с.
13. Korolyuk V.S. Stochastic models of systems. Dordrecht: Springer science+business media, 1999, 185 p.
14. Копп В.Я. Стохастические модели автоматизированных производственных систем с временным резервированием / В.Я. Копп, Ю.Е. Обжерин, А.И. Песчанский // Севастополь: Изд-во СевГТУ, 2000. - 284 с.
15. Ross S.M. Introduction to probability models. USA, Elsevier academic press, 2006, 800 p.
16. Rabiner L.R. A tutorial on hidden Markov models and selected applications in speech recognition. Proceedings of the IEEE, 1989, vol. 77, no.2, pp. 257-286, DOI: 10.1109/5.18626.
17. Cappe' O., Moulines E., Ryde'n T. Inference in hidden markov models. New York, Springer-Verlag, 2005, 653 p.
18. Obzherin Y.E., Nikitin M.M., Sidorov S.M. Hidden Markov model based on signals from blocks of semi-markov system's elements and its application for dynamics analysis energy systems. Smart innovation, systems and technologies. Singapore: Springer, 2021, vol. 187, p. 477-486.
19. Обжерин Ю.Е. Скрытая марковская модель технического обслуживания по возрасту элементов системы / Ю.Е. Обжерин, М.М. Никитин, С.М. Сидоров // Методические вопросы исследования надежности больших систем энергетики. 92-е заседание семинара, учрежденного при ИСЭМ СО РАН. В 3-х книгах. - Иркутск: Изд-во ИСЕМ СО РАН, 2020. - С. 89-98.
Обжерин Юрий Евгеньевич. Доктор технических наук, профессор, кафедра высшей математики Севастопольского государственного университета. Исследовательские интересы связаны с теорией случайных процессов, теорией надежности, теорией массового обслуживания, математическим моделированием систем, теорию скрытых марковских моделей. AuthorID: 766997, SPIN: 7589-4167, ORCID: 0000-0003-1180-1084. objsev@mail.ru, 299053, Российская Федерация, г. Севастополь, ул. Университетская, 33, кафедра «Высшая математика».
Сидоров Станислав Михайлович. Кандидат технических наук, доцент кафедр «Высшая математика» и «Интеллектуальные сети энергоснабжения» Севастопольского государственного университета, ученый секретарь Института природно-технических систем. Исследовательские интересы включают прикладную математику, теорию надежности, математическое моделирование систем, теорию скрытых марковских моделей. AuthorID: 889130, SPIN: 2766-5222, ORCID: 0000-0002-9785-918. xaevec@mail.ru, 299011, Российская Федерация
г. Севастополь, ул. Ленина, д. 28.
Федоренко Сергей Николаевич. Старший преподаватель кафедры Приборные системы и автоматизация технологических процессов Севастопольского государственного университета. Исследовательские интересы включают автоматизацию и математическое моделирование систем. AuthorID: 930695, SPIN: 4100-2953, ORCID: 0000-0003-3080-7668, sergey.fedor@mail.ru, 299055, г. Севастополь, проспект генерала Острякова,
д.167, кв.32.
UDC 519.217
DOI: 10.38028/ESI.2022.28.4.005
Semi-Markov and hidden Markov models of autonomous wind diesel complex functioning
Yuriy E. Obzherin, Stanislav M. Sidorov, Sergey N. Fedorenko
Sevastopol State University, Russia, Sevastopol, xaevec@mail.ru
Abstract. The problem of assessing the reliability and analyzing the functioning of an autonomous wind-diesel complex is considered. Its semi-Markov and hidden Markov models are constructed. A hidden Markov model is used to find performance estimates for an autonomous wind-diesel complex and predict its states based on the resulting signal vector. The results of the study are obtained in a general form and are invariant with respect to the laws of distribution of random variables describing the elements. They make it possible to predict the operating modes of an autonomous wind-diesel complex.
Keywords: wind diesel complex, semi-Markov model, hidden Markov model, performance analysis, forecasting
Acknowledgements: The study was supported by the grant of the President of the Russian Federation for state support of young Russian scientists - Candidates of Sciences No. MK-329.2022.4.
References
1. Bobrov A.V., Tremjasov V.A. Vetrodizel'nye kompleksy v decentralizovannom jelektrosnabzhenii: monografija [Wind-diesel complexes in decentralized power supply: monograph]. Krasnojarsk, Izd-vo SFU [SFU Publishing House], 2012, 216 p.
2. Grigor'eva O.A., Krivenko T.V., Tremjasov V.A. Analiz nadezhnosti avtonomnogo vetrodizel'nogo kompleksa [Reliability analysis of an autonomous wind-diesel complex]. Nauchno-tekhnicheskie vedomosti Sankt-Peterburg-skogo gosudarstvennogo politekhnicheskogo universiteta [Scientific and technical statements of the St. Petersburg State Polytechnic University], 2016, no. 2 (243), pp. 45-52.
3. Salima Meziane, Omar Feddaoui, Riad Toufouti. Modeling and simulation of hybrid wind-diesel power generation system. International journal of renewable energy, 2013, vol. 8, no. 2, pp. 49-58.
4. Sebastián R., Nevado A. Study and simulation of a wind hydro isolated microgrid. Energies, 2020, vol. 13(22), DOI: 10.3390/en13225937.
5. Sebastián R. Modeling, simulation and control of wind diesel power systems. Energies, 2022, no. 15, DOI: 10.3390/en15051712.
6. Shklyarskiy Y.E., Batueva D.E. Operation mode selection algorithm development of a wind-diesel power plant supply complex. Journal of mining institute, 2022, no. 253, pp. 115-126, DOI: 10.31897/PMI.2022.7.
7. Baikov A.I., Dar'enkov A.B., Sosnina E.N. Simulation modeling of a wind-diesel power plant. Russian electrical engineering, 2018, vol. 89, pp. 161-167, DOI: 10.3103/S1068371218030033.
8. Limnios N., Nikulin M. Recent advances in reliability theory: methodology, practice, and inference. New York: Springer science+business media, 2000, 514 p.
9. Rudenko Y.N., Ushakov I.A. Nadjozhnost' sistem jenergetiki, 2-e izd [Reliability of power systems, 2nd edition]. Novosibirsk, Nauka [The Science], 1989, 328 p.
10. Obzherin Y.E. Polumarkovskie i skrytye markovskie i polumarkovskie modeli sistem jenergetiki [Semi-Markov and hidden Markov and semi-Markov models of energy systems]. Izvestiya RAN. Energetika [Proceedings of the Russian Academy of Sciences. Power Engineering], 2019, no 5, pp. 26-32.
11. Obzherin Y.E., Sidorov S.M., Nikitin M.M. Analiz nadezhnosti i fazovoe ukrupnenie sistem s pojelementnymi nakopiteljami [Analysis of reliability and phase merging of systems with component-wise storages]. Izvestiya RAN. Energetika [Proceedings of the Russian Academy of Sciences. Power Engineering], 2019, no. 6, pp. 66-77.
12. Koroljuk V.S., Turbin A.F. Processy markovskogo vosstanovlenija v zadachah nadjozhnosti system [Markov recovery processes in problems of system reliability]. Kiev, Nauk. Dumka [Scientific Thought], 1982, 236 p.
13. Korolyuk V.S., Stochastic models of systems. Dordrecht, Springer science+business media, 1999, 185 p.
14. Kopp V.Ja., Obzherin Yu.E., Peschanskij A.I. Stohasticheskie modeli avtomatizirovannyh proizvodstvennyh sistem s vremennym rezervirovaniem [Stochastic models of automated production systems with time redundancy]. Sevastopol', Izd-vo SevGTU [Publishing House of SevSTU], 2000, 284 p.
15. Ross S.M. Introduction to probability models. USA: Elsevier academic press, 2006, 800 p.
16. Rabiner L.R. A tutorial on hidden Markov models and selected applications in speech recognition. Proceedings of the IEEE, 1989, vol. 77, no.2, pp. 257-286, DOI: 10.1109/5.18626.
17. Cappe' O., Moulines E., Ryde'n T. Inference in hidden Markov models. New York, Springer-Verlag, 2005, 653 p.
18. Obzherin Yu.E., Nikitin M.M., Sidorov S.M. Hidden Markov model based on signals from blocks of semi-markov system's elements and its application for dynamics analysis energy systems. Smart innovation, systems and technologies, Singapore, Springer, 2021, vol. 187, pp. 477-486.
19. Obzherin Yu.E., Nikitin M.M., Sidorov S.M. Skrytaja markovskaja model' tehnicheskogo obsluzhivanija po voz-rastu jelementov sistemy [Hidden Markov model of maintenance by age of system elements]. Metodicheskie vo-prosy issledovanija nadezhnosti bol'shih sistem jenergetiki. 92-e zasedanie seminara, uchrezhdennogo pri ISJeM SO RAN. V 3-h knigah [Methodological issues of reliability research of large energy systems. 92nd meeting of the seminar established at ISEM SB RAS. In 3 books], 2020, pp. 89-98.
Obzherin Yuri Evgenievich. Doctor of Technical Sciences, Professor, Department of Higher Mathematics, Sevastopol State University. Research interests are related to the theory of random processes, reliability theory, queuing theory, mathematical modeling of systems, the theory of hidden Markov models. AuthorlD: 766997, SPIN: 7589-4167, ORCID: 0000-0003-1180-1084. objsev@mail.ru, 299053, Russian Federation, Sevastopol, st. Universitetskaya, 33, Department of Higher Mathematics.
Sidorov Stanislav Mikhailovich. Candidate of Technical Sciences, Associate Professor of the Departments of Higher Mathematics and Intelligent Power Supply Networks of the Sevastopol State University, Scientific Secretary of the Institute of Natural and Technical Systems. Research interests include applied mathematics, reliability theory, mathematical modeling of systems, the theory of hidden Markov models. AuthorID: 889130, SPIN: 2766-5222, ORCID: 00000002-9785-918. xaevec@mail.ru, 299011, Russian Federation, Sevastopol, st. Lenina, d. 28.
Fedorenko Sergey Nikolaevich. Senior Lecturer, Department of Instrument Systems and Automation of Technological Processes, Sevastopol State University. Research interests include automation and mathematical modeling of systems. AuthorID: 930695 SPIN: 4100-2953 ORCID: 0000-0003-3080-7668 sergey.fedor@mail.ru Статья поступила в редакцию 26.08.2022; одобрена после рецензирования 06.10.2022; принята к публикации 25.11.2022.
The article was submitted 08/26/2022; approved after reviewing 10/06/2022; accepted for publication 11/25/2022.