CC BY
86
DOI 10.5862/JEST.243.5 УДК 621.311.245
О.А. Григорьева, Т.В. Кривенко, В.А. Тремясов
АНАЛИЗ НАДЕЖНОСТИ АВТОНОМНОГО ВЕТРОДИЗЕЛЬНОГО КОМПЛЕКСА
O.A. Grigorieva, T.V. Krivenko, V.A. Tremyasov RELIABILITY ANALYSIS OF AN AUTONOMOUS WIND-DIESEL COMPLEX
Использование ветродизельных комплексов в автономных системах электроснабжения позволяет заместить часть дорогостоящего дизельного топлива возобновляемым источником энергии. Основной недостаток возобновляемых источников энергии — нестабильность, зависимость их функционирования от природных ресурсов, что влияет на надежность работы системы электроснабжения. Актуален анализ надежности работы систем генерации на основе возобновляемых источников энергии. В работе предложен логико-вероятностный метод для оценки надежности автономного ветродизельного комплекса, учитывающий его особенности и условия функционирования. Произведен расчет надежности ветродизельного комплекса, расположенного на севере Красноярского края.
ВОЗОБНОВЛЯЕМЫЕ ИСТОЧНИКИ ЭНЕРГИИ; ВЕТРОДИЗЕЛЬНЫЙ КОМПЛЕКС; НАДЕЖНОСТЬ; ДЕРЕВО ОТКАЗОВ; ОТКАЗ; ПОКАЗАТЕЛИ НАДЕЖНОСТИ
Using wind-diesel complexes in autonomous electricity supply systems allows to substitute the expensive part of diesel fuel by renewable energy sources. The renewable energy sources, however, are unstable energy sources whose behavior differs much from the conventional ones. They rely in their functioning on nature resources, which affects the stable operation of the electric supply system. Therefore, analyzing the reliability of generation systems based on renewable energy sources is of immediate interest. The logical and probabilistic method has been suggested for evaluating the reliability of an autonomous winddiesel complex considering the specifics and conditions of its functioning. The reliability of a wind-diesel complex located at the north of the Krasnoyarsk Krai is computed.
RENEWABLE ENERGY SOURCE; WIND-DIESEL COMPLEX; RELIABILITY; FAULT TREE; FAILURE; RELIABILITY INDICATORS.
Большая часть обширной территории России с малой плотностью населения до сих пор не присоединена к централизованным энергетическим системам: по статистическим данным не подключены к электрическим сетям около 10 млн населения, живущего на Северных территориях, Дальнем Востоке и в некоторых других регионах. Население этих регионов получает энергию в основном от дизельных генераторов (ДГ). Однако электроснабжение от ДГ децентрализованных областей обладает рядом недостатков. К ним относят малый ресурс и высокий износ энергооборудования, загрязнение окружающей
среды. Кроме того, необходимое топливо для ДГ завозится в отдаленные населенные пункты по зимникам, водными путями, а иногда даже авиацией, что делает производство электроэнергии очень дорогим [1]. Решением этой проблемы является применение автономных ветродизельных комплексов (ВДК). Применение в децентрализованных областях ветроэнергетических установок (ВЭУ) совместно с ДГ позволяет частично заменить дорогостоящее дизельное топливо возобновляемым источником энергии (ВИЭ).
Однако ВИЭ — это источники, мощность которых зависит от потенциала природных ре-
сурсов. Так, функционирование ВЭУ определяется режимом поступающего ветрового потока и зависит от его энергетических характеристик. Основными энергетическими характеристиками ветрового потока являются скорость ветра и годовой ход среднемесячных скоростей ветра [2]. Из-за изменчивости скорости ветра нестабильность выработки электроэнергии такими установками влияет на надежность системы электроснабжения. Анализ надежности систем генерации на основе ВИЭ — новая и актуальная задача.
Цель нашей работы — развитие логико-вероятностного метода на основе дерева отказов (ДО) для анализа надежности автономного ветроди-зельного комплекса (ВДК).
Для анализа надежности ВИЭ применяются различные методы [3—9]: аналитические; метод пространства состояний (марковские процессы); метод статистического моделирования (метод Монте-Карло).
Однако эти методы имеют недостатки применительно к оценке показателей надежности систем генерации на основе ВИЭ. Аналитические методы не учитывают стохастический характер природных ресурсов. Основная сложность при использовании метода пространства состояний связана с высокой размерностью математической модели. В общем случае число R возможных состояний, в которых может находиться система, равно
к
R = Ъм^,
I=1
где k — число групп с различными возможными
состояниями элементов; М — число возможных
' 1
состояний элементов ^й группы; N. — число элементов в ^й группе. Поэтому при анализе надежности систем генерации с большим количеством элементов применение данного метода требует составления сложного графа с большим количеством состояний, что существенно усложняет процедуру анализа.
Метод статистического моделирования (метод Монте-Карло) требует специального программного обеспечения. В большинстве случаев этот метод может быть использован в качестве эталонного.
Логико-вероятностный метод расчета надежности с использованием ДО является дедуктив-
ным методом и применяется в тех случаях, когда число различных видов отказов системы невелико [10]. Метод заключается в построении логической блок-схемы (дерево отказов), в которой используются элементы, подразделяющие и связывающие большое число событий и условий: логические символы и символы событий. События и условия формулируются специалистами на основе технического проекта системы генерации с помощью анализа ее поведения (при возникновении различных отказов и режимов) и записи в виде конъюнкций (П =.) и дизъюнкций (и=+).
Количественная оценка показателей надежности систем генерации на основе ВИЭ с помощью ДО может быть произведена различными способами. Для анализа надежности ВДК предлагается использовать логико-вероятностную модель на основе формирования логической функции отказа (ФО). Модель учитывает восстановление и профилактику элементов системы [10].
Метод ДО на основе формирования логической ФО позволяет дифференцировать отказы на кратковременные и длительные. Кратковременные отказы (КО) — это кратковременное отключение на время оперативных переключений, а длительные отказы (ДО) — это отключение на время аварийного ремонта или замены оборудования. К кратковременным отключениям приводят отказы коммутационных аппаратов (КА), отходящих линий и элементов, которые могут быть отсоединены от распределительного устройства (РУ) с помощью рубильников и разъединителей. К длительным отключениям приводят отказы генерирующего оборудования, преобразовательных и аккумулирующих установок, коммутационной аппаратуры, а также элементов, не имеющих разъединителей или рубильников между ними и РУ.
Метод ДО учитывает элементарные события, приводящие к отказу ВДК, а также такие сложные события, как отказ одного элемента во время аварийного ремонта другого или отказ элемента при аварийном ремонте другого и событии отказа в срабатывании коммутационного аппарата.
С помощью предложенного метода учитываются ветровые условия, влияющие на работу ВДК. При построении ДО для ВДК учитывают-
ся состояния, когда наблюдается расчетная скорость ветра ут|п < V < vтax , при которой ВЭУ развивает мощность от нуля до номинальной, и случаи нерасчетной скорости ветра vmax < V; V < vmin, при которой установка остановлена и не вырабатывает электроэнергию. Отказы по причине ветровых условий обозначаются символом (V).
Чтобы получить выражение ФО для конечного события необходимо построить ДО. Формирование ФО происходит путем выполнения указанных в ДО операций логического сложения («или») и умножения («и») над символами первичных элементарных событий и состояний начиная с нижнего уровня. В результате ФО для конечного события будет представлена логической суммой (дизъюнкцией) конъюнкций первого порядка (элементарные события), конъюнкций второго порядка (произведение события отказа и события отказа в срабатывании, произведения события и состояния отказа) и конъюнкций третьего порядка (произведение состояний, событий отказов и событий отказов в срабатывании):
ФО = и я и (я П *ос *) и (я п Уу) и (уу п у )х
i к у у
х и (уП у} П *ос к)
ук
где у, у,, у — элементарные события (отказы
1 к j
элементов 1, j, к); хо с к — события отказов в срабатывании КА, устройств релейной защиты и автоматики (РЗиА) для элемента к; у1, у — состояние элементов (аварийные ремонты), в обозначении которых используется тильда («волна») над буквой.
Здесь каждая комбинация символов в дизъюнкции представляет собой сечение отказов или условие его возникновения в виде конъюнкции первого, второго или третьего (и более высоких) порядков.
Полученная ФО подвергается сокращению по правилам алгебры логики. Так сокращаются одинаковые конъюнкции и конъюнкции, содержащие в себе конъюнкции более низкого порядка. Сокращаются также конъюнкции четвертого и пятого порядка, если в ФО содержатся конъюнкции третьего, второго и первого порядков.
Логическая ФО позволяет получить выражения для аналитического расчета частоты Л и вероятности q отказа системы по известным частотам и вероятностям отказов ее элементов. Для перехода от логической ФО к выражениям для расчета частоты и относительной длительности отказа, необходимо заменить знаки логического умножения П и сложения и на знаки умножения (•) и суммирования (+) соответственно, а символы событий у заменить символами частоты Х(у), символы событий хо с — символами условной вероятности Q(xос К), символы состояний у — символами вероятности состояний д (у). При этом вероятность состояния простоя элемента системы равна
д (У ) = Ч У )т( у),
где Х(у) — интенсивность отказа элемента системы; т( у) — среднее время восстановления функционирования элемента системы.
Выражение для оценки частоты отказа имеет вид
лвдк=ХЧ Уi )+ХЧ Ук )2 (к)+
i к
+ ХЧУ )я(уу ) + ХЧУ: >(у) +
1 А
+ХЧ у)? (у )з (*ос *)
ук
где 1, j, к — элементы системы; Х(у.), ^(ук), — интенсивность отказов элементов системы; Q(xo с к) — условная вероятность отказа срабатывания КА и устройств РЗиА.
Выражение для оценки вероятности (относительная длительность) состояния отказа имеет вид
?ВДК = X Ч У)т (У г) + X Ч у к)Q (С к) т(У к П *о с к)
г к
+ХЧ у, 9 (У у )т( У П Уу ) + ХЧ у у > (у, )т( у у ПУг)+
у уг
+ХЧу) 9(Уу) Q (Хоск) т(У,),
гук
где т (ук П хос к) — среднее время восстановления функционирования при отказе в отключении повреждения на присоединении; т (у П У у) — среднее время восстановления функционирования при совпадении отказа одного элемента системы с аварийным ремонтом другого элемента.
0,4 кВ
Ж
X
Р
G
Z
I
"X
Ж АБ
Рис. 1. Расчетная схема автономного ветродизельного комплекса: Ж — ветроэнергетическая установка; G — дизельный генератор; I — инвертор; АБ — аккумуляторная батарея; X, Y, Z — выключатели; Р — нагрузка
В качестве примера применения предложенного метода рассмотрим схему ВДК (рис.1), функционирующего в северных районах Красноярского края, с учетом характеристик ветра [2].
В состав автономного энергокомплекса входят ветроэнергетическая установка и дизельный генератор. Также в ВДК включена аккумуляторная батарея (АБ), которая позволяет сглаживать
возможные колебания мощности ВЭУ, связанные с переменным характером скорости ветра. ДО для оценки надежности электроснабжения ВДК приведено на рис. 2. ДО построено для конечного события — отказа функционирования всей системы ВЭУ—ДГ—АБ.
При расчетной скорости ветра формулируются элементарные события, которые приводят к отказу системы в зависимости от того, какой отказ (длительный или кратковременный), а также сложные события (отказ одного элемента во время аварийного ремонта другого и при учете отказа в срабатывании КА).
Как указывалось выше, при анализе надежности ВИЭ необходимо учитывать зависимость их функционирования от поступления природных ресурсов. В ДО (рис. 2) отражено влияние ветровых условий на работу ветроустановки посредством логического знака «запрет» (шестиугольник). Применение в ДО логического знака «запрет» позволяет учесть событие, которое происходит с какой-то определенной вероятностью [10]. На рис. 2 входным событием V, помещенным под знаком «запрет», является простой ВЭУ по причине ветровых условий (отсутствие ветра, ураган и др.). Условное событие Q, расположен-
Отказ всей системы ВЭУ—ДГ—АБ
Рис. 2. Дерево отказов для оценки надежности автономного ВДК
ное на рис. 2 сбоку от логического знака, — это ограничительное условие, которое представляет собой вероятностный весовой фактор. Ограничительное условие характеризуется величиной q условной вероятности нерасчетной скорости ветра. Величина вероятности нерасчетной скорости ветра определяется с помощью гистограммы распределения скоростей ветра в месте расположения ВДК. Гистограмма распределения скоростей ветра характеризует повторяемости скоростей ветра за исследуемый период [2]. Она показывает, какую часть времени в течение рассматриваемого периода ветер имел определенную скорость.
Функция отказа для события «Отказ системы ВЭУ—ДГ—АБ» имеет следующий вид:
ФО(ВДК)ДО = WGZОС + вШОС +
ЮС
+УШ ОС + W (А + />ос.;
ФО(ВДК)КО = X + Y + Z + ШОС + GYОС + ИОС;
ФО(ВДК)к = (В + У)^ с, где G, Ж, А, I — коды отказавшего элемента в соответствии со схемой ВДК; G, Ж, А, I — коды элемента, находящегося в аварийном ремонте; X, У, Z — отказ автоматических выключателей; ХО С, УО С, ZО С — отказ выключателей системы при появлении заявки на срабатывание; V — отказ ветроустановки по причине ветровых условий.
Для учета отказа ветроустановки по причине ветровых условий необходимо определить вероятность
4* = Чу Ч,
где qr — коэффициент простоя ВЭУ по причине погодных условий; q — условная вероятность появления нерасчетной скорости ветра.
Коэффициент простоя ВЭУ по причине погодных условий определяется по формуле
qv =
V +М
\-е~(Ху +Му У
где Ху = 1/70, год-1, — интенсивность появления нерасчетной скорости ветра; = 1/Тм, год1, — интенсивность восстановления расчетной скорости ветра [11]; Т0 и Ты — периоды погоды соответственно с нерасчетной и расчетной скоростью ветра [2].
По ФО составлены расчетные выражения для показателей надежности: для частоты отказов ВДК
Л(ВДК )д о ЯGQOCz +^gqwq ocz + + ^xqgqocz ocz +^w (яа + Я1 )QOCY;
Л(ВДК)к о = Хх Ч Ху Ч X г Ч
Ч К@о С У Ч с ?
Л(ВДК )у ={х6 ;
для относительной длительности отказов ВДК
Я (ВДК )до = V V ЯоQ ocZ +^0 TGЯWQOCZ +
+ TXЯGQ ocz + ocz +
+ tw(ЯА +Я1У2OCY; V (ВДК)КО = ТХ + V% + Т +
я (вдк)у =(*<?Т + 4TY )яV0OCZ. Суммарная частота отказов ВДК
Л(ВДК) = лд о ч лк 0 ч л,. Суммарная относительная длительность состояния отказа ВДК
q(ВДК) = qд о Ч qК 0 Ч ду. Чтобы получить численное значение показателей надежности ВДК, необходимо в полученные формулы подставить исходные данные для ВДК, представленные в табл. 1, 2 [10, 12]. Необходимые данные для учета в ДО влияния ветровых условий на работу ВДК для северных районов Красноярского края следующие:
д............................................................0,05
Т...........................................................0,42
о '
Тн..........................................................0,58
V 1/год................................................2,38
йг, 1/год................................................ 1,72
Чу...........................................................0,57
Чу........................................................... 0,0285
Результаты расчета показателей надежности по вышеприведенным формулам: Л(ВДК)д о = 3,110-7 1/год; Л(ВДК)К О = 3,59-10-3 1/год; Л(ВДК)у = 3,1410-5 1/год; Ч(ВДК)д о = 1,4940-9; Ч(ВДК)к о = 2,5Ы0-6; д(ВДК)у = 5,7-10-9. Окончательно получаем Л(ВДК) = 3,6240-3 1/год; д(ВДК) = 2,5240-6 .
Таблица 1
Показатели надежности элементов ВДК
Код элементов Значения показателей надежности
Элементы ВДК X(y), 1/год т( У), год )
W Ветроустановка Vestas V27 225 кВт 1,1 3,28-10-3 3,61-10-3
G Дизельный генератор (2±0,5)-10-2 (1±0,5)-10-2 (2±0,5)-10-4
I Инвертор (5±4)-10-2 (2±1)-10-4 Ы0-5
A Аккумуляторная батарея (2±1)-10-3 (1±0,5)-10-3 2-10-6
X, Y, Z Автоматический выключатель (1±0,5)-10-3 (2±1)-10-4 2-10-7
Таблица 2
Условная вероятность отказа срабатывания элементов
Код элемента Элемент ВДК Отказы автоматических выключателей: Условная вероятность отказа срабатывания (Q к)
X, Y, Z Автоматический выключатель В отключении КЗ Во включении (5±2)-10-4 (1±0,5>10-3
Выводы
Предложен логико-вероятностный метод дерева отказов для анализа надежности ВДК, который позволяет учесть особенности функционирования комплекса и погодные условия.
Сформулирована целесообразность дифференциации отказов ВДК на кратковременные, длительные и отказы по причине изменения ветровых условий.
Метод учитывает такие сложные события, происходящие в ВДК, как отказ одного элемента энергокомплекса во время аварийного ремонта другого с учетом отказа в срабатывании КА.
Надежность ВДК определяется с помощью следующих показателей: интенсивности отказа (Л) и относительной длительности отказа Эти показатели могут быть использованы при выборе оптимальною варианта структуры ВДК.
СПИСОК ЛИТЕРАТУРЫ
1. Куликова Л.В., Меновщиков Ю.А. Нетрадиционные и возобновляемые источники энергии / Алт. гос. техн. ун-т им. И.И. Ползунова. Барнаул: Изд-во АлтГУ, 2005. 365 с.
2. Старков А.Н., Ландберг Л., Безруких П.П., Бо-рисенко М.М. Атлас ветров России. М.: Можайск— Терра, 2000. 560 с.
3. Abouzahr I., Ramakumar R. Loss of power supply probability of stand-alone wind electric conversion systems: A closed form solution approach // IEEE Transactions on Energy Conversion. 1990. Vol. 5, № 3. P. 445-451.
4. Abouzahr I., ramakumar г. An approach to assess the performance of utility-interactive wind electric conversion systems // IEEE Transactions on Energy Conversion. 1991. Vol. 6, № 4. P. 627-638.
5. Bakirtzis A.G. A probabilistic method for the evaluation of the reliability of stand-alone wind energy systems // IEEE Transactions on Energy Conversion. 1992. Vol. 7, № 1. P. 99-107.
6. Leite A.P., borges C.L.T., Falcao D.M. Probabilistic wind farms generation model for reliability studies applied to Brazilian sites // IEEE Transactions on Power Systems. 2006. Vol. 21, № 4. P. 493-501.
7. billinton R., bagen, Cue Yu. Reliability evaluation of small stand-alone wind energy conversion systems using a time series simulation model. IEE Proceedings Generation, Transmission and Distribution. 2003. Vol. 150, № 1. P. 96-100.
8. billinton R., bagen B. A sequential simulation method for the generating capacity adequacy evaluation of small stand-alone wind energy conversion systems //
IEEE Canadian Conference on Electrical and Computer Engineering CCECE. 2002. Vol. 1. P. 72-77.
9. Li W.Y, Bagen B. Reliability evaluation of integrated wind/diesel/storage systems for remote locations // IEEE 11th International Conference on Probabilistic Methods Applied to Power Systems (PMAPS). 2010. P. 791-795.
10. Гук Ю.Б., Синенко М.М., Тремясов В.А. Расчет надежности схем электроснабжения. Л.: Энер-
гоатомиздат, Ленингр. отд-ние, 1990. 208 с.
11. Бобров А.В., Тремясов В.А. Ветродизельные комплексы в децентрализованном электроснабжении. Красноярск: Изд-во Сиб. федер. ун-та, 2012. 216 с.
12. Sienergy [Электронный ресурс] // URL: http: //www.sienergy. co.uk/WT/ downloads/Reliability-and-Condition-Monitoring-of-Wind-Turbines.pdf (дата обращения 07.02.2016).
REFERENCES
1. Kulikova L.V., Menovshchikov Yu.A. Netraditsionnye i vozobnovlyayemye istochniki energii [Nonconventional and renewable energy sources] / Altai state technical university name of I. I. Polzunov. Barnaul: ASU Publ.,
2005. 365 s. (rus)
2. Starkov A.N., Landberg L., bezrukikh P.P., borisenko M.M. Atlas vetrov Rossii [Wind atlas of Russia] Moscow: Mozhaisk-Terra, 2000. 560 s. (rus)
3. Abouzahr I., Ramakumar R. Loss of power supply probability of stand-alone wind electric conversion systems: A closed form solution approach. IEEE Transactions on Energy Conversion. 1990. Vol. 5, № 3. P. 445-451.
4. Abouzahr I., Ramakumar R. An approach to assess the performance of utility-interactive wind electric conversion systems. IEEE Transactions on Energy Conversion. 1991. Vol. 6, № 4. P. 627-638.
5. bakirtzis A.G. A probabilistic method for the evaluation of the reliability of stand-alone wind energy systems. IEEE Transactions on Energy Conversion. 1992. Vol. 7, № 1. P. 99-107.
6. Leite A.P., borges C.L.T., Falcao D.M. Probabilistic wind farms generation model for reliability studies applied to Brazilian sites. IEEE Transactions on Power Systems,
2006. Vol. 21, № 4. P. 493-501.
7. billinton R., bagen, Cue Yu. Reliability evaluation of small stand-alone wind energy conversion systems using a time series simulation model. IEEE Proceedings Generation, Transmission and Distribution, 2003. Vol. 150, № 1. P 96-100.
8. billinton R., bagen b. A sequential simulation method for the generating capacity adequacy evaluation of small stand-alone wind energy conversion systems, IEEE Canadian Conference on Electrical and Computer Engineering CCECE. 2002. Vol. 1. P. 72-77.
9. Li W.Y, bagen B. Reliability evaluation of integrated wind/diesel/storage systems for remote locations. IEEE 11th International Conference on Probabilistic Methods Applied to Power Systems (PMAPS), 2010. P. 791-795.
10. Guk Yu. B., Sinenko M.M., Tremyasov V.A. Raschet nadezhnosti skhem elektrosnabzheniia [Reliability calculation of schemes power supply]. Leningrad: Energoatomizdat, 1990. 208 s. (rus)
11. bobrov A.V., Tremyasov V.A. Vetrodizelnye kompleksy v detsentralizovannom elektrosnabzhenii [Winddiesel complexes in the decentralized power supply]. Krasnoyarsk: Siberian Federal University, 2012. 216 s. (rus)
12. Sienergy. Available: http://www.sienergy.co.uk/WT/ downloads/Reliability-and-Condition-Monitoring-of-Wind-Turbines.pdf (Accessed 07.02.2016)
СВЕДЕНИЯ ОБ АВТОРАХ/AUTHORS
ГРИГОРЬЕВА Ольга Анатольевна — кандидат технических наук доцент, начальник управления
аспирантуры, докторантуры и АНПК Сибирского федерального университета.
660041, Россия, г. Красноярск, пр. Свободный, 79.
E-mail: Grigorieva.00@mail.ru
GRIGORIEVA Olga A. - Siberian Federal University.
79, Svobodny Prospect, Krasnoyarsk, Russia, 660041.
E-mail: Grigorieva.00@mail.ru
КРИВЕНКО Татьяна Витальевна — аспирант Сибирского федерального университета.
660041, Россия, г. Красноярск, пр. Свободный, 79.
E-mail: tanya-1991mir@mail.ru
KRIVENKO Tatiana V. — Siberian Federal University.
79, Svobodny Prospect, Krasnoyarsk, Russia, 660041..
E-mail: tanya-1991mir@mail.ru
ТРЕМЯСОВ Владимир Анатольевич — кандидат технических наук профессор Сибирского федерального университета.
660041, Россия, г. Красноярск, пр. Свободный, 79. E-mail: emf_tva@mail.ru
TREMYASOV Vladimir A. - Siberian Federal University. 79, Svobodny Prospect, Krasnoyarsk, Russia, 660041. E-mail: emf tva@mail.ru
