CC BY
6
Технические науки — от теории к практике № 4 (52), 2016г_'_
СЕКЦИЯ
«ЭНЕРГЕТИКА И ЭНЕРГЕТИЧЕСКИЕ ТЕХНИКА И ТЕХНОЛОГИИ»
ЧАСТОТНАЯ ДЕКОМПОЗИЦИЯ НА ОСНОВЕ ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЯ ДЛЯ ИССЛЕДОВАНИЯ НЕСИНУСОИДАЛЬНЫХ РЕЖИМОВ СИСТЕМ ЭЛЕКТРОСНАБЖЕНИЯ
Осипов Дмитрий Сергеевич
канд. техн. наук,
доц. кафедры «Электроснабжение промышленных предприятий» Омского государственного технического университета,
РФ, г. Омск E-mail: ossipovdmitriy@list. ru
Долгих Надежда Николаевна
аспирант кафедры «Электроснабжение промышленных предприятий» Омского государственного технического университета,
РФ, г. Омск E-mail: nabal2006@list. ru
Горовой Сергей Анатольевич
студент кафедры «Электроснабжение промышленных предприятий» Омского государственного технического университета,
РФ, г. Омск E-mail: serega1593728@mail.ru
Технические науки — от теории к практике _№ 4 (52), 2016 г
FREQUENCY DECOMPOSITION BASED ON WAVELET TRANSFORM FOR STUDY OF NON-SINUSOIDAL MODES OF POWER SUPPLY SYSTEMS
Dmitriy Osipov
candidate of Science, assistant professor of Omsk State Technical University, Russia, Omsk
Nadezhda Dolgikh
postgraduate student of Omsk State Technical University,
Russia, Omsk
Sergey Gorovoy
student of Omsk State Technical University,
Russia, Omsk
АННОТАЦИЯ
Вейвлет преобразование и преобразование Фурье находят широкое применение для анализа показателей качества электроэнергии. Главная задача определения спектрального состава, действующих значений и коэффициента гармонического искажения возникает при исследовании нестационарных режимов. Данная статья представляет результаты численного моделирования применения пакетного вейвлет преобразования для частотной декомпозиции сигналов. Показано применение для расчета действующего значения тока в каждом частотном диапазоне.
ABSTRACT
The wavelet transform and Fourier transform are popular techniques for power quality analysis in electrical power systems. They are mainly used to identify power harmonics and power disturbances and to estimate power quantities in the presence of nonstationary power components such as root-mean-square (RMS) values and total harmonic distortions. This paper presents a numerical example of application of wavelet packet transform (WPT) for estimation of standard harmonic groups. The advantage of WPT may be useful in measuring the RMS value of current in each harmonic component.
Ключевые слова: пакетное вейвлет преобразование; высшие гармоники; действующее значение гармоник.
Технические науки — от теории к практике № 4 (52), 2016г._
Keywords: wavelet packet transform; power system harmonics; root-mean square values.
Развитие цифровой обработки сигналов в электроэнергетике предопределило разработку и внедрение нового поколения устройств автоматического управления, релейной защиты и диагностики, в основе которых лежат новые математические методы анализа динамических систем. В последние годы большое количество исследователей-энергетиков обращают внимание на методы вейвлет преобразования (ВП) для решения актуальных задач электроэнергетической отрасли [1].
Идея применения вейвлет преобразования для исследования показателей качества электроэнергии находит всё более широкое применение благодаря качественному отличию при анализе поведения сигнала не только в частотной, но и в временной области [2]. Благодаря локализации вейвлет функции появляется возможность не только идентифицировать гармонический состав параметров режима систем электроснабжения (токов, напряжений, мощностей), но и исследовать время изменения указанных несинусоидальных режимов.
Вейвлет преобразование заданной функции f(t) по составляющим в различных масштабах и частотных диапазонах может быть реализовано путём свёртки f (t) и вейвлет функции *¥(t) [3]:
где: k и j сдвиг по времени и масштабирующий коэффициент,
Принято различать непрерывное ВП и дискретное ВП. В данной работе продемонстрируем частотную декомпозицию несинусоидального тока на основе разновидности дискретного ВП, называемого пакетным ВП.
Если сигнал тока, содержащий три частотные компоненты (50, 150, 250 Гц) задан непрерывной величиной представлен функцией:
i(t) = Im1 sin (®lt + Ф1 ) + 4з sin (®3t + Ф3 ) + 45 sin (®5t + Ф5 ) (2)
(1)
соответственно, а l нормирующий показатель.
Технические науки — от теории к практике _№ 4 (52), 2016 г
где: 1т1 = л/2; 1т3 = л/2/2; 1т5 = л/2/2 - амплитудное значение гармоник тока, (А);
о - угловая частота, соответствующей гармоники. По теореме Котельникова, определим частоту дискретизации для исходного сигнала:
F, > F
d m
(3)
где: F - максимальная частота сигнала, Гц.
^ m ' '
По условию (3) принимаем частоту дискретизации Fd = 800 Гц. Непрерывный и дискретный сигналы тока, представлены на рисунке 1.
Рисунок 1. Непрерывный и дискретный сигналы тока
Схема пакетного вейвлет разложения сигнала до уровня ] = 2 представлена на рисунке 2. Индексом т на схеме обозначены узлы вейвлетного дерева разложения. В общем случае каждый узел содержит некоторое количество элементов (вейвлет коэффициентов), каждый из которых имеет индексы к, диапазон изменения которых внутри каждого узла т может быть определён:
2-jnm < к.т < 2-Jn(m + 1)-1
к < к < к
min j ,m m
где: n - количество элементов (отсчётов) исходного сигнала
Технические науки — от теории к практике № 4 (52), 2016г_
n = T ■ F,
(6)
Для рассматриваемого численного примера n = 0.02 ■ 800 = 16 отсчётов. Численный эксперимент проведем с применением вейвлета Ингрид Добеши 24 порядка (db42). От типа выбранного вейвлета во многом зависит точность полученных результатов [4], так, наиболее удобным для пояснения принципа вейвлет преобразования является вейвлет Хаара, в то же время, вейвлеты Добеши показывают в задачах исследования несинусоидальности токов и напряжений в системах электроснабжения более высокую точность.
г= о-юо гц г = 100-200 Гц г=2оо-зоогц с=зоо-40«гц
Рисунок 2. Схема пакетного вейвлет разложения сигнала тока
В качестве фильтров низких (ФНЧ) и фильтров высоких (ФВЧ) частот применяются коэффициенты, определяемые для каждого типа вейвлета. Для выбранной глубины разложения до уровня ] = 2 можно
определить, что в каждом узле т будет представлено по 4 вейвлет коэффициента:
Технические науки — от теории к практике _№ 4 (52), 2016 г
/2„(к) =
/2,(к) =
-1.7619 2.2123 1.7619 -2.2123 -0.7928 1.1510 0.7928 -1.1510
Т т
1 1
1 2 3 4
1 1 т
1 *
/ (к) =
1.3042 0.5875 -1.3042 -0.5875
Рисунок 3. Величины вейвлет коэффициентов второго уровня
разложения
Действующее значение тока в каждом частотном диапазоне, определяемом по схеме (рисунок 2) можно определить по формуле:
1
=, - I ^ (к)
(7)
Для нахождения действующего значения тока третьей гармоники (частота 150 Гц) по формуле (7) необходимо просуммировать квадраты коэффициентов в узле т = 1:
1 -==161 ы)=,/т^
(1,3042 )2 +(0,5875)2 + + (-1,3042)2 +(-0,5875)2
= 0,5057 (А)
к
к
Технические науки — от теории к практике № 4 (52), 2016г._
Результаты расчета сведём в таблицу, кроме того целесообразно сравнение с действительными (истинными) значениями, каждой частотной компоненты, определяемое выражением (2).
Частотный диапазон 50 Гц 150 Гц 250 Гц
Действующее значение аналогового сигнала, А 1 0,5 0,5
Действующее значение через вейвлет коэффициенты, А 0,9999 0,4941 0,5057
Погрешность вычисления, % 0,01 1,18 1,14
Следует отметить, что в моделируемом сигнале отсутствует сигнал в коридоре частот 300 < f < 400 Гц, что хорошо отражается в равенстве нулю коэффициентов вейвлет разложения в узле j = 2, m = 3 .
Список литературы:
1. Аббакумов А.А. Разработка методики и алгоритмов идентификации отклонений от нормативов параметров качества электроэнергии в системах электроснабжения: дис. ... канд. техн. наук: 05.13.18 -математическое моделирование, численные методы и комплексы программ / Саранск: Мордовский гос. ун-т им. Н.П. Огарева, 2005. - 180 с.
2. Мисриханов А.М. Применение методов вейвлет-преобразования в электроэнергетике // Автоматика и телемеханика. - 2006. - № 5. -С. 5-23.
3. Определение коэффициентов искажения синусоидальности формы кривой тока по вейвлет коэффициентам / Н.Н. Долгих, Ф.В. Катеров, Д.С. Осипов, Е.В. Птицына, Л.А. Файфер // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. - 2015. - № 113. - С. 814-825.
4. Яковлев А.Н. Введение в вейвлет-преобразования: учеб. пособие / -А. Новосибирск: Изд-во НГТУ, 2003. - 104 с.
