Определение коэффициентов искажения синусоидальности формы кривой тока по вейвлет коэффициентам

Размещено на http://www.allbest.ru/

ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ИСКАЖЕНИЯ СИНУСОИДАЛЬНСТИ ФОРМЫ КРИВОЙ ТОКА ПО ВЕЙВЛЕТ КОЭФФИЦИЕНТАМ

Для расчета несинусоидальных режимов систем электроснабжения применяется метод наложения, при котором на основании преобразования Фурье определяются отдельные частотные компоненты сигнала тока (напряжения). В условиях стационарного режима (неизменного во времени) данный подход обеспечивает приемлемые результаты моделирования и анализа. В случае изменения режима работы электроустановки (изменения числа, потребляемой мощности, напряжения в узлах нагрузок) использование преобразования Фурье может исказить представление об исследуемом сигнале. Связано это, прежде всего, с отсутствием локализации базовых функций преобразования Фурье в пространстве времени. Проиллюстрируем данное утверждение осциллограммами тока в случае неизменного и нестационарного режима работы электроустановки. В первом случае (рисунок 1) на протяжении исследуемого интервала времени спектральный состав тока не изменяется и является суммой трех частотных компонент.

Рисунок 1- График и частотный спектр стационарного сигнала.

При анализе нестационарных процессов возникает эффект «растекания спектра» (рисунок 2), т.е. вместо узкого пика получается сложный спектр, в котором в общем случае могут содержаться все возможные частоты. Причина растекания спектра состоит в том, что ДПФ неявно подразумевает периодическое продолжение анализируемого фрагмента сигнала. Если на рассматриваемом промежутке укладывается целое число периодов синусоиды (это эквивалентно условию совпадения ее частоты с одной из частот анализа), периодически продолженный сигнал также будет непрерывной синусоидой, в спектре которой содержится единственная частота. Если же число периодов на интервале анализа не является целым, при периодическом продолжении сигнала непрерывность синусоиды окажется нарушенной и спектр «растечется» [1].

Рисунок 2. График и частотный спектр нестационарного сигнала

Для уменьшения влияния эффекта «растекания» используются «весовые» или оконные функции. При этом исходный сигнал перед преобразованием Фурье умножается на некоторую функцию, убывающую от середины к краям. Это ослабляет влияние разрывов, образующихся на стыке фрагментов сигнала при его периодизации [1]. Для работы с нестационарными сигналами так же часто пользуются модифицированной версией ПФ, которая носит название Оконное преобразование Фурье (ОПФ). Основная проблема ОПФ заключается в сложности подбора ширины используемой оконной функции (окна).

В настоящее время все более широкое применение для анализа нестационарных процессов в электроэнергетике и других областях техники получает вейвлет преобразование, которое «даёт лучшую, чем оконное преобразование Фурье, возможность, рассмотреть высокочастотные явления с коротким сроком жизни» [2].

Дискретное вейвлет преобразование определяется скалярным произведением двух кусочно-непрерывных функций - исследуемой функции тока и вейвлет фукнции на отрезке , при :

Система базисных функций , впервые предложенных А. Хааром (названных впоследствии вейвлетами), обладала главной особенностью вейвлетов: локальной областью определения (ограниченными носителями), ортогональностью и единичной нормой, нулевым средним и автомодельностью [3]. Семейство вейвлет функций образуется путём изменения масштаба и сдвига базисной функции , что в общем виде может быть представлено:

Изменение масштаба и сдвига базисной функции иллюстрирует алгоритм разложения, представленный на рисунке 4 для уровня разложения и на рисунке 5 для уровня разложения. Вейвлет преобразование позволяет получить временную и частотную информацию, сужая окно для выделения коротких высокочастотных участков или расширяя его для анализа длительных низкочастотных колебаний [4]. Различают непрерывное (НВП) и дискретное (ДВП) вейвлет преобразование. Так, представленное на рисунке 3 НВП исследуемого нестационарного сигнала даёт четкое представление в трёхмерном пространстве об изменении амплитуды, частоты и времени.

Рисунок 3. Непрерывное вейвлет преобразование исследуемого нестационарного сигнала тока

Расчет вейвлет-спектра при НВП сопряжен со значительными вычислительными трудностями при непрерывном изменении параметров j и k. Множество функций шj,k избыточно, поэтому для сохранения возможности восстановления исходного сигнала из его преобразования необходима дискретизация этих параметров. Зачастую дискретизация реализуется через степени двойки.

В настоящее время Хаара были вытеснены другими типами вейвлетов, но вместе с тем, вейвлеты Хаара являются достаточно удобными для формирования теоретического базиса исследования сигналов с применением вейвлет преобразования. Вейвлет Хаара можно задать функционально следующим образом:

При обработке числовых массивов данных дискретные вейвлеты используются, как правило, в паре со связанными с ними дискретными масштабирующими функциями [3].

Численный эксперимент дискретного вейвлет преобразования провеём для нестационарного сигнала тока, представленного на рисунке 2. На основании теоремы Котельникова (теоремы отсчетов) непрерывный сигнал тока i(t), спектр которого не содержит частот выше fm , полностью определяется дискретной последовательностью своих мгновенных значений {ik}, k=0,1,…, N-1, отсчитываемых через интервалы времени ?t:

,,

где ?t и fд - интервал (шаг) и частота дискретизации.

Таким образом, дискретизированный с шагом ?t сигнал можно определить выражением:

,

где д(t) - дельта функция.

Приняв частоту дискретизации

(Гц)

получаем дискретный сигнал тока, представляющий собой вектор-строку из 16 элементов (рисунок 4). Перемножая полученный дискретный сигнал с масштабирующей функцией , сдвигаемой по оси времени, так чтобы полностью покрыла исследуемый сигнал, получаем аппроксимирующие коэффициенты первого уровня разложения . Аналогичным образом, проводя перемножение с вейвлетом Хаара получаем детализирующие коэффициенты . Для второго (рисунок 5) и последующего уровней разложения необходимо по формуле (2) получить новое семейство отмасштабированных и сдвинутых функций и .

Рис. 4. Первый уровень вейвлет разложения исследуемого дискретного сигнала

На практике для анализа сигналов часто пользуются быстрым вейвлет преобразованием или алгоритмом Малла из-за возможности вычислять вейвлет коэффициенты без интегрирования, при помощи алгебраических операций свертки. Фиксированное количество коэффициентов h0(k) и h1(k) связывают масштабируемую функцию и вейвлет одного разрешения с масштабирующей функцией на следующем более низком разрешении [4]:

Рисунок 5. Второй уровень вейвлет разложения исследуемого дискретного сигнала

Представленная на рисунке 6 схема декомпозиции сигнала позволяет выделить аппроксимирующие и детализирующие вейвлет коэффициенты, характеризующие определенный спектр частот. А сигнал тока (напряжения, мощности) может быть представлен через вейвлет коэффициенты: синусоидальность вейвлет электрический ток

Следует отметить, что полученные при разложении сигнала (при выбранной частоте дискретизации) вейвлет коэффициенты будут характеризовать различные частотные компоненты исследуемого сигнала в соответствии с таблицей 1.

Таблица 1. Частотные диапазоны вейвлет коэффициентов

Рисунок 5. Исходный сигнал и вейвлет разложение

На основании равенства Парсеваля может быть вычислено действующее значение тока в заданном диапазоне частот. В нашей работ, действующее значение основной частоты (50 Гц) может быть получено через аппроксимирующие и детализирующие вейвлет коэффициенты третьего (в нашем случае максимального) уровня разложения:

Действующее значение токов высших гармоник определяется через детализирующие коэффициенты первых двух уровней разложения:

где, N - количество отсчётов исследуемого дискретного сигнала

Полученные таким образом действующие значения могут быть подставлены в формулу для расчёта коэффициента искажения синусоидальности:

Задача применения вейвлет анализа для идентификации искажений синусоидальных сигналов напряжения и тока в системах промышленного электроснабжения с обоснованием выбора вейвлет носителя решена в работе [5]. Автор исследует применимость более десятка различных типов вейвлетов в задачах определения показателей качества электроэнергии.

В нашей работе ограничимся четырьмя типами вейвлетов, представленных в таблице 2. Определим с помощью указанных вейвлетов действующие значения токов основной частоты и высших гармоник с последующих нахождением коэффициента искажения синусоидальности.

Таблица 2. Сравнение точности расчета по некоторым вейвлетам

Из таблицы видно, что тип используемого вейвлета во многом определяет точность расчета коэффициента искажения. Наиболее приемлемые результаты показали вейвлеты Добеши 24 и 42 порядков (db24, db42), погрешность не превысила 2%. Предлагаемые методы определения гармонического состава тока или напряжения совместно с современными способами активной фильтрации высших гармоник [6,7] помогут решать задачи обеспечения качества электроэнергии, используя преимущества вейвлет преобразования.

Библиографический список

1. Cергиенко А. Б. Цифровая обработка сигналов/АБ Сергиенко //СПб: ПИТЕР. - 2002.

2. Добеши И. Десять лекций по вейвлетам. - Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001. - 464 с.

3. Нагорнов, О.В. Вейвлет-анализ в примерах: учеб. пособие / О.В. Нагорнов, [и др.]. - М.:НИЯУ МИФИ, 2010. - 120 с.

4. Захарова, Т.В. Вейвлет-анализ и его приложения: учеб. пособие / Т.В. Захарова, О.В. Шестаков. - М.:ИНФРА-М, 2012. - 150 с.

5. Аббакумов, А.А. Разработка методики и алгоритмов идентификации отклонений от нормативов параметров качества электроэнергии в системах электроснабжения: дис. … канд. техн. наук: 05.13.18 - математическое моделирование, численные методы и комплексы программ / Аббакумов Андрей Александрович ; рук. работы С.А. Федосин. - Саранск: Мордовский гос. ун-т им. Н.П. Огарева, 2005. - 180 с.

6. Долингер, С.Ю. Схематические решения активной фильтрации кривой тока в четырёхпроводной трехфазной сети для обеспечения качества электрической энергии / С.Ю. Долингер, В.Н. Горюнов, А.А.Планков, О.А. Сидоров // Омский научный вестник. - 2013. - №3(103). - С. 214-217.

7. Горюнов, В.Н. Определение управляющего воздействия активного фильтра гармоник / В.Н. Горюнов, Д.С. Осипов, А.Г. Лютаревич // Электро. Электротехника, электроэнергетика, электротехническая промышленность. - 2009. - №6 - С. 20-24.

Аннотация

ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ИСКАЖЕНИЯ СИНУСОИДАЛЬНСТИ ФОРМЫ КРИВОЙ ТОКА ПО ВЕЙВЛЕТ КОЭФФИЦИЕНТАМ

Долгих Надежда Николаевна

ФГБОУ ВПО «Омский Государственный технический университет», г. Омск, Россия Катеров Филипп Викторович

Магистрант ФГБОУ ВПО «Омский Государственный технический университет», г. Омск, Россия

Осипов Дмитрий Сергеевич к.т.н., доцент

ФГБОУ ВПО «Омский Государственный технический университет», г. Омск, Россия

Птицина Елена Витальевна д.т.н., профессор

ФГБОУ ВПО «Омский Государственный технический университет», г. Омск, Россия

Файфер Лилия Андреевна Магистрант

ФГБОУ ВПО «Омский Государственный технический университет», г. Омск, Россия

В статье представлен алгоритм определения действующего значения тока и коэффициента искажения синусоидальности формы кривой тока в электрических системах на основании дискретного вейвлет преобразования. Показатели качества электроэнергии имеет существенное значение в случае нестационарных режимов. В этом случае удовлетворительных результатов не могут обеспечить ни подходы, основанные на анализе частотного спектра, использующие быстрое преобразование Фурье, ни методы, основанные на локализации в временнуй области. Представленный алгоритм предполагает разделение сигнала тока по частотным диапазонам, отвечающих за определенный спектр гармоник. Предложенный алгоритм позволяет избежать эффекта растекания спектра. Эффективность алгоритма была проверена компьютерным моделированием

Ключевые слова: ВЕЙВЛЕТ - ПРЕОБРАЗОВАНИЕ, ДЕЙСТВУЮЩЕЕ ЗНАЧЕНИЕ, КОЭФФИЦИЕНТ ИСКАЖЕНИЯ СИНУСОИДАЛЬНОСТИ

USING DISCRETE WAVELET TRANSFORM FOR CALCULATION OF TRANSIENT REGIMES OF POWER SUPPLY SYSTEMS

Dolgikh Nadezhda Nikolaevna postgraduate student

Omsk State Technical University, Omsk, Russia

Katerov Philipp Viktorovich master student

Omsk State Technical University, Omsk, Russia

Osipov Dmitriy Sergeevich Cand.Tech.Sci., associate professor

Omsk State Technical University, Omsk, Russia

Pticina Elena Vital'evna Dr.Sci.Tech., professor Omsk State Technical University, Omsk, Russia

Faifer Liliya Andreevna master student

Omsk State Technical University, Omsk, Russia

This article presents an algorithm based on the discrete wavelet transform for the analysis of current root mean square (RMS) value and total harmonic distortion (THD) in power systems. Power quality indices play an important role in case of non-stationary distorted waveforms, where neither a frequency-domain-based approach using fast Fourier transform tools nor a time-domain-based approach using real time data give satisfactory results. The algorithm proposed decomposes the current waveforms into uniform frequency bands corresponding to the odd harmonic components of the signal. The proposed algorithm overcomes the spectra leakage problem. Computer simulations verified the effectiveness of the proposed algorithm

Keywords: WAVELET TRANSFORMS, ROOT MEAN, SQUARE, TOTAL HARMONIC DISTORTION VALUE

Размещено на Allbest.ru