Алгоритм расчета потерь мощности, обусловленных высшими гармониками и интергармониками на основе вейвлет-преобразования
Автор: Осипов Дмитрий Сергеевич, Лютаревич Александр Геннадьевич, Ткаченко Всеволод Андреевич, Логунова Яна Юрьевна
Журнал: Вестник Южно-Уральского государственного университета. Серия: Энергетика @vestnik-susu-power
Рубрика: Электроэнергетика
Статья в выпуске: 1 т.23, 2023 года.
Бесплатный доступ
В статье приводится модернизация алгоритма расчета дополнительных потерь мощности в электрических сетях в несинусоидальных нестационарных режимах. Основой метода является математический аппарат вейвлет-преобразования. Отличием от ранее существующих методов является предварительный выбор оптимальной вейвлет-функции по условию энергии спектра гармонической группы. В работе предложено совместить гармонические и интергармонические группы, определенные в действующей нормативно-технической документации с полосой пропускания вейвлет-функции. Вторым определяющим фактором при выборе оптимальной вейвлет-функции является минимальное среднеквадратичное отклонение исходного и восстановленного сигнала. Предложенный алгоритм позволяет для нестационарных режимов определять действующие значения высших гармоник и интергармоник, а также длительность (время) присутствия высокочастотных компонент в сигнале тока. Полученные величины необходимы для расчета дополнительных потерь мощности и энергии в токоведущих частях. Повышение точности расчета потерь при несинусоидальных нестационарных режимах позволит корректировать технико-экономические обоснования внедрения устройств для фильтрации высших гармоник и интергармоник.
Вейвлет-преобразование в электроэнергетике, высшие гармоники, потери мощности, интергармоники, качество электрической энергии
Короткий адрес: https://sciup.org/147240391
IDR: 147240391 | DOI: 10.14529/power230104
Текст научной статьи Алгоритм расчета потерь мощности, обусловленных высшими гармониками и интергармониками на основе вейвлет-преобразования
Неизбежный рост доли электроприемников, имеющих нелинейную вольт-амперную характеристику, приводит к искажению синусоидальности формы кривой тока и напряжения в электрических сетях и системах электроснабжения. Как правило, эти режимы являются не экономичными, вызывают дополнительные потери в токоведущих частях, могут приводить к преждевременному выходу из строя электрооборудования. Потери, обусловленные несинусоидальностью, являются причиной дополнительного нагрева токоведущих частей. Действующий ГОСТ 32144–2013 в части несину-соидальности определяет допустимые значения гармонических составляющих напряжения в зависимости от класса напряжения и порядка высшей гармоники. Кроме высших гармоник в электрических сетях возможно возникновение интергармоник, однако их уровень в ГОСТ 32144–2013 «находится на рассмотрении» (п. 4.2.4.2). В странах Европейского Союза уровень интергармоник определяется стандартом IEC 61000-4-7:2, где предельно допустимый уровень интергармоник определен как 0,2 % для каждого номинального напряжения. В США уровень интергармоник регламентируется специальной номограммой стандарта IEEE Std 519 для каждого класса напряжения.
Одной из причин появления в сетях интергармоник являются несинхронные переключения устройств силовой электроники [1]. Источниками интергармоник в системах электроснабжения могут быть преобразователи частоты, дуговые печи, частотно-регулируемый электропривод и др. [1].
Вопросы анализа, моделирования и учета интергармоник вызывают существенный интерес в научном сообществе. В статье [2] приводится имитационное моделирование асинхронного электродвигателя, работающего на переменную периодическую нагрузку. При этом авторы определяют возникающие в обмотках статора высшие гармоники и интергармоники. В результате исследования авторы приходят к выводу о целесообразности применения активных фильтров [2, с. 43]
Зависимость уровня интергармоник в питающей сети системы освещения и кратковременной дозы фликера приведена в работе [3]. В ходе исследований было установлено, что интергармони- ческие искажения влияют на колебания напряжения в сети, что негативно отражается на системах освещения.
Основные расчетные выражения для определения частот интергармоник при работе различных типов преобразователей, включая циклоконверторы, приведены в исследовании [4].
Как отмечают авторы [5], под определение интергармоник «попадают процессы, имеющие совершенно разный генезис», при этом подобные явления пока мало изучены.
Анализ российских и зарубежных нормативных документов, регламентирующих как показатели качества электроэнергии в целом, так и не-синусоидальность в частности, приводится в статье [6]. Авторы считают целесообразным разработку отечественных нормативных значений допустимых пределов гармонических и интергармонических составляющих тока [6, с. 31]
Таким образом, возникает задача оценки негативного влияния интергармоник с целью внесения рекомендаций об их нормировании в электрических сетях России.
-
1. Основные положения исследования
Математический аппарат вейвлет-преобразо-вания получает все более широкое применение для анализа динамики сложных нестационарных несинусоидальных процессов. Теория вейвлетов применяется для модуляции напряжения на выходе автономных инверторов напряжения [7]. Вейвлет-анализ токов асинхронных двигателей и уровней вибрации позволил авторам [8] произвести оценку остаточного ресурса машины с достаточно высокой (до 96 %) достоверностью. Вейвлет-анализ находит применение для анализа несинусоидальных режимов, в том числе расчета гармоник и интергармоник [9, 10].
В настоящей работе для реализации алгоритма расчета потерь мощности в токоведущих частях от высших гармоник и интергармоник будет использован метод вейвлет-анализа.
В статье [11] авторы применяют частотновременной контур быстрого преобразования Фурье для анализа интергармоник. Применение вейвлет-преобразования позволит упростить решение, предложенное в [11]. В статье [12] представлены разработки метода идентификации интергармоник в сети промышленного предприятия при наличии частотно-регулируемого привода.
Как известно, любая вейвлет-функция представляет собой цифровой фильтр с определенной амплитудно-частотной характеристикой, которая в общем случае определяется типом вейвлета и глубиной разложения при реализации алгоритма вейвлет-преобразования. Таким образом, каждая вейвлет-функция может быть представлена характеристиками цифрового фильтра. Как и любой реальный фильтр, вейвлеты вносят некоторые искажения в полосе пропускания и не полностью подавляют частоты в полосе заграждения. На рис. 1 представлены амплитудно-частотные характеристики вейвлет-функций Ингрид Добе-ши 10 (обозначено пунктирной линией) и 44 порядков (обозначено сплошной линией).
В настоящий момент для решения задач спектрального анализа токов и напряжений при несинусоидальных режимах электрических сетей широко применяется метод преобразования Фурье. В случае анализа нестационарных режимов возникает эффект «растекания спектра», который может быть минимизирован с помощью оконного преобразования Фурье. В соответствии с действующей нормативной документацией [13] для устранения эффекта «растекания спектра» используется принцип группирования гармоник и интергармоник в соответствии с формулой
2 1 2 N* 2-1 2
Y g , h = T Y C ,( Nh ) - N / 2 + t Y C ,( Nh ) + k +
2 k = ( - N / 2) + 1
+ 2 Y C ,( Nh ) + N 2 , (1)
где YC ,( Nh ) _ N2 - среднеквадратичное значение частотной компоненты;
( Nh ) + k - номер частотной компоненты;
Yg , h – результирующее среднеквадратичное значение гармонической группы.
Принцип формирования групп и подгрупп гармоник и интергармоник представлен на рис. 2 и 3.
Для выбора оптимальной вейвлет-функции предлагается совместить гармонические и интергармонические группы, определенные в [13] с полосой пропускания вейвлет-функции таким образом, чтобы узлу aj , k дерева вейвлет-разложения
соответствовала полоса пропускания нижних час-
тот, а следующему узлу a j k + 4 дерева полоса пропускания верхних частот вейвлет-функции (рис. 4).
Если в качестве одного из условий выбора оптимальной вейвлет-функции для анализа несинусоидальных режимов электрических сетей принять минимальное среднеквадратичное отклонение исходного сигнала ik тока (напряжения) от функции тока, восстановленного в результате обратного вейвлет-преобразования ik , то это позволит говорить о минимальной потере информации о сигнале.
Характерным признаком совпадения полос пропускания и заграждения вейвлет-функций при совмещении гармонических и интергармонических групп в соответствии с рис. 4 может быть энергия спектра, определяемая для каждого узла дерева вейвлет-разложения. Энергия спектра гармонической группы En должна быть максимальной именно в границах узла вейвлет-разло-жения aj , k . В таком случае критерий для выбора оптимальной вейвлет-функции может быть записан в следующем виде:
g =
<
Zt (kk - k )2
k =1-- > min;
n
tk
!■ k |2+ Z ЬГ,
E n n = E ;
k = n
k = n
^ ^ ^
zb-k| '-1 к k =-^ j=J k =-^
Dj , k 2
max.

Рис. 1. АЧХ вейвлетов Добеши db10, db44
Fig. 1. Frequency response of Daubechies wavelets db10, db44

Рис. 2. Принцип образования гармонических и интергармонических групп [13, с. 15] Fig. 2. The formation of harmonic and interharmonic groups [13, p. 15]

Рис. 3. Принцип образования интергармонических подгрупп [13, с. 17]
Fig. 3. The formation of interharmonic subgroups [13, p. 17]
Частотная локализация узла aj,k
Частотная локализация узла aj,k+i

Рис. 4. Совмещение полос пропускания и заграждения вейвлет-функций с гармоническими группами по ГОСТ 30804.4.30–2013 [14]
Fig. 4. The combination of the passbands and the blocking of wavelet functions with harmonic groups according to GOST 30804.4.30–2013 [14]
Предложенный критерий (выражение (2)) позволит на предварительном этапе выбрать оптимальную вейвлет-функцию для расчета действующих значений токов высших гармоник и интергармоник.
-
2 Полученные результаты
Как было отмечено выше, несинусоидальные режимы могут приводить к росту потерь мощности и электроэнергии, а также к повышению температуры токоведущих частей. Необходимость учета интергармоник в расчете потерь мощности определяется в статье [15].
В статье [16] представлен метод, в основе которого лежат преобразования Гильберта, который позволяет проводить спектральный анализ гармоник и интергармоник, при этом производится сдвиг оконной функции и интерполяция.
В зависимости от выходной частоты преобразователя частоты f 0 спектральный состав тока (напряжения) может быть различным. В общем случае спектр частот может быть представлен в виде основной гармоники и двух боковых частот:
i ( t ) = I m • sin(2 n ft ) + I m • m • sin((2 n f + Q ) t ), (3) где I m – амплитуда основной частоты, А;
I m • m - амплитуда боковых частот, А;
f – частота питающей сети, Гц;
Q - угловая частота огибающей, рад/с.
Наличие боковых частот вызывает эффект колебания напряжения.
В качестве исходного сигнала в предлагаемом алгоритме используются мгновенные значения токов, которые формирует измерительный прибор. Следует отметить, что частота дискретизации измерительного прибора должна соответствовать теореме Котельникова, в противном случае возможны искажения в процессе цифровой обработки сигналов [17].
В настоящей работе производится модернизация метода, представленного в работе [18] с учетом дополнительных условий выбора оптимального типа вейвлет-функции по выражению (2).
На первом этапе производится вейвлет-преобразование сигнала тока и на основании детализирующих вейвлет-коэффициентов высоких порядков определяются интервалы присутствия той или иной гармонической (интергармонической) составляющей. Участок времени, на котором происходит существенное изменение сигнала, характеризуется значительным увеличением детализирующих коэффициентов (рис. 5).
Аналогичный подход использования детализирующих вейвлет-коэффициентов высоких порядков для определения границ изменения режима при анализе отклонения напряжения был использован в статье [19, c. 133] и в работе [20, c. 276–277].
Далее каждый элемент матрицы тока возводится в квадрат, в результате чего можно записать следующее выражение i2 (t) = Im • sin2 (tot) + Im • m2 • sin2 ((to+Q)t) =
_ 2 I 1 2 I - , 2
= I m '[ 2 + m J + 2 I m
• m • cos Q • t + I m • m • cos 2 Q t +
, 2 I 2 1 I 2
+ Im - \m — — I • cos 2 to t - Im • m • cos( to - П ) -
2 1 • m
— Im * m • cos( to + П ) —------ • cos 2 • ( to — П ) —
-
— — cos2 • ( to + Q ). (4)
Из выражения (4) видно, что при возведении исследуемой функции в квадрат в сигнале появляется постоянная составляющая тока и частоты Ω, 2Ω, 2ω и др. Практический интерес представляет частота огибающей - Q , что позволит определить частоту интергармоники.
В результате частотной декомпозиции сигнала тока (выражение (4)) и оценки спектральной плотности (метод Уэлча, Барлетта и др.) получаем частоту огибающей Ω, из которой определяется частота интергармоники (рис. 6).
Модернизированный алгоритм расчета потерь мощности от высших гармоник и интергармоник представлен на рис. 7.

Рис. 5. График детализирующего вейвлет-коэффициента d1
Fig. 5. Detail of the wavelet coefficient d1

Рис. 6. Поиск огибающей частоты методом Уэлча
Fig. 6. Finding the frequency envelope using the Welch method

Рис. 7. Предлагаемый алгоритм Fig. 7. The proposed algorithm
Результаты численного эксперимента Numerical experiment results
f , Гц |
50 |
130 |
150 |
230 |
250 |
350 |
370 |
450 |
470 |
I m , А |
85,685 |
14,985 |
18,23 |
8,89 |
17,94 |
11,927 |
8,87 |
7,47 |
4,48 |
£ , % |
0,02 |
0,1 |
1,27 |
3,42 |
0,33 |
0,61 |
2,64 |
0,4 |
0,44 |
A P , кВт |
4,23 |
0,13 |
0,19 |
0,05 |
0,19 |
0,08 |
0,05 |
0,03 |
0,01 |
Список литературы Алгоритм расчета потерь мощности, обусловленных высшими гармониками и интергармониками на основе вейвлет-преобразования
- Hanzelka Z., Bien A. “Power Quality Application Guide: Harmonics & Interharmonics”, A guide material by Leonardo Power Quality Initiative, Copper Development Association, 2004.
- Определение интергармоник тока асинхронного двигателя с переменной периодической нагрузкой / Б.А. Авдеев, С.Г. Черный, И.С. Моисеев, А.А. Жиленков // Электротехника. 2022. № 6. С. 39–44. DOI: 10.53891/00135860_2022_6_39
- Eidson B., Halpin М. An evaluation of the extent of correlation between interharmonic and voltage fluctuation measurements // IEEE Transactions on Power Delivery. 2016. Vol. 31, iss. 2. P. 753–760. DOI: 10.1109/tpwrd.2015.2480715
- Interharmonics: theory and modeling / A. Testa [et al.] // IEEE Transactions on Power Delivery. 2007. Vol. 22, iss. 4. P. 2335–2348. DOI: 10.1109/TPWRD.2007.905505
- Сивоконь В.П., Лапшов Д.В., Белов О.А. Диагностические признаки нестандартного проявления нелинейности в электрических сетях // Вестник КамчатГТУ. 2019. № 48. С. 18–27. DOI: 10.17217/2079-0333-2019-48-18-27
- Макашева С.И., Пинчуков П.С. Качество тока: аспекты оценки и нормирования // Вестник ЮУрГУ. Серия «Энергетика». 2020. Т. 20, №4. С. 23–35. DOI: 10.14529/power200403
- Доброскок Н.А. Лавриновский В.С. Спектральный анализ базовых алгоритмов широтно-импульсного управления без обратной связи для двухуровневых преобразователей частоты // Электротехника. 2021. № 3. С. 21–26.
- Степанов П.И. Закураев В.В. Модель и алгоритм оценки остаточного ресурса электромеханического оборудования на основе комплексного анализа тока и вибрации // Южно-Сибирский научный вестник. 2021. № 2 (36). С. 122–127. DOI: 10.25699/SSSB.2021.36.2.015
- Barros J., Diego R.I., de Apráiz M. Applications of Wavelet Transform for Analysis of Harmonic Distortion in Power Systems: A Review // IEEE Transactions on Instrumentation and Measurement. 2012. Vol. 61, iss. 10. P. 2604–2611. DOI: 10.1109/TIM.2012.2199194
- Morsi W., El-Hawary M. Time-frequency single-phase power components measurements for harmonics and inter-harmonics distortion based on wavelet packet transform; Part II: case studies and results // Canadian Journal of Electrical and Computer Engineering. 2013. Vol. 35 (1). P. 8–14. DOI: 10.1109/cjece.2010.5783379
- Hui J., Xu W., Yang H. A Method to determine the existence of genuine interharmonics // IEEE Transactions on Power Delivery. 2012. Vol. 27 (3). P. 1690–1692. DOI: 10.1109/tpwrd.2012.2198313
- Sheshyekani K., Fallahi G., Hamzeh M., Kheradmandi M. A General Noise-Resilient Technique Based on the Matrix Pencil Method for the Assessment of Harmonics and Interharmonics in Power Systems // IEEE Transactions on Power Delivery. 2017. Vol. 32, no. 5. p. 2179–2188. DOI: 10.1109/TPWRD.2016.2625329
- ГОСТ 30804.4.7–2013. Совместимость технических средств электромагнитная. Общее руководство по средствам измерений и измерениям гармоник и интергармоник для систем электроснабжения и подключаемых к ним технических средств. Введ. 01.01.2014. М.: Стандартинформ, 2013. 40 с.
- ГОСТ 30804.4.30–2013. Электрическая энергия. Совместимость технических средств электромагнитная. Методы измерений показателей качества электрической энергии. Введ. 01.01.2014. М.: Стандартинформ. 2014. 52 с.
- Бирюлин В.И., Куделина Д.В., Горлов А.Н. Анализ нагрева кабельных линий токами высших гармоник и интергармоник // Вестник Казанского государственного энергетического университета. 2020. № 2 (46). С. 61–67.
- Xiong J., Wang B. Measuring power system harmonics and interharmonics by envelope spectrum analysis» // Przeglad Elektrotechniczny. 2010. Vol. 86, no. 12. P. 319–324.
- Обеспечение когерентности цифровой обработки сигналов токов и напряжений электроэнергетических систем при снижении частоты дискретизации / А.Л. Куликов, А.Б. Лоскутов, П.В. Ильюшин, А.А. Севостьянов // Электричество. 2020. № 8. С. 5–16. DOI: 10.24160/0013538020208516
- Разработка метода расчета потерь мощности в токоведущих частях при наличии интергармоник / Д.С. Осипов, Д.В. Коваленко, Л.А. Файфер и др. // Омский научный вестник. 2017. № 4 (154). С. 60–65.
- Latran M.B., Teke A. A novel wavelet transform based voltage sag/swell detection algorithm // International journal of electric power and energy systems. 2015. Vol. 71. P. 131–139. DOI: ijepes.2015.02.040
- Voltage sag and swell detection and segmentation based on independent component analysis / E.A. Nagata, D.D. Ferreira, C.A. Duque, A.S. Cequeira // Electric Power Systems Research. 2018. Vol. 155. P. 274–280. DOI: 10.1016/j.epsr.2017.10.029
- Жежеленко И.В., Саенко. Ю.Л. Показатели качества электроэнергии и их контроль на промышленных предприятиях. М.: Энергоатомиздат, 2000. 252 с.
- Манусов В.З., Хрипков В.В., Фролова В.В. Сравнительный анализ математических моделей для определения коэффициента увеличения активного сопротивления проводников от высших гармоник // Научные проблемы транспорта Сибири и Дальнего Востока. 2018. № 1. С. 184–188.
- Влияние высших гармоник тока на режимы работы кабелей распределительной сети 380 В / В.Н. Тульский, И.И. Карташев, Р.Р. Насыров, М.Г. Симуткин // Промышленная энергетика. 2013. № 5. С. 39–44.
- Экспериментальные исследования электромагнитной совместимости современных электроприводов в системе электроснабжения металлургического предприятия / А.А. Николаев, Г.П. Корнилов, Т.Р. Храмшин и др. // Вестник Магнитогорского государственного технического университета им. Г.И. Носова. 2016. Т. 14, № 4. С. 96–105. DOI: 10.18503/1995-2732-2016-14-4-96-105