Научный журнал
Международный журнал прикладных и фундаментальных исследований
ISSN 1996-3955
ИФ РИНЦ = 0,564

ИДЕНТИФИКАЦИЯ КОЛЕБАНИЙ НАПРЯЖЕНИЯ В ЭЛЕКТРОЭНЕРГЕТИЧЕСКИХ СИСТЕМАХ НА ОСНОВЕ ПАКЕТНОГО ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЯ

Долгих Н.Н. 1 Киселёв Б.Ю. 1 Осипов Д.С. 1 Горовой С.А. 1
1 ГОУ ВО «Омский государственный технический университет»
В статье приводится метод анализа колебаний напряжения, основанный на вейвлет преобразовании (ВП), который имеет более высокую точность при локализации сигнала по времени. Вейвлет анализ хорошо зарекомендовал себя для исследования периодических сигналов, как для стационарных режимов, так и при переходных процессах. Колебания напряжения могут быть изучены с применением различных типов вейвлет функций для повышения точности результатов. В рассматриваемом методе предлагается применение частотной декомпозиции сигнала с применением вейвлет фильтров для идентификации времени начала и окончания колебаний напряжения по высокочастотным детализирующим коэффициентам, а также определить частоту и амплитуду искажающего сигнала. В статье приводится численный эксперимент анализа колебания напряжения с применением пакетного вейвлет преобразования.
качество электроэнергии
колебания напряжения
вейвлет преобразование
1. ГОСТ 32144–2013. Электрическая энергия. Совместимость технических средств электромагнитная. Нормы качества электрической энергии в системах электроснабжения общего назначения. Введ. с 01.07.2014. – М.: Стандартинформ, 2014. – 16 с.
2. Дмитриева Е.Н. Совершенствование цифровой фликер модели для оценки колебаний напряжения / Е.Н. Дмитриева, В.А. Топчий // Электричество. – 2010. – № 7. – С. 20–25.
3. Долгих Н.Н. Применение дискретного вейвлет преобразования для расчета нестационарных режимов систем электроснабжения // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. – 2015. – № 113. – С. 286–296.
4. Жежеленко И.В. Электромагнитная совместимость потребителей: моногр. / И.В. Жежеленко, А.К. Шидловский, Г.Г. Пивняк и др. ‒ М.: Машиностроение, 2012. ‒ 351 с.
5. Салтыков А.В. Метод расчета колебаний напряжения на шинах питания одиночной и параллельных дуговых сталеплавильных печей // Вестник Самарского государственного технического университета. Серия: Технические науки. – 2006. – № 41. – С. 161–165.
6. Черепанов В.В. Расчет колебаний напряжения электрической сети при работе асинхронных двигателей с резкопеременной нагрузкой / В.В. Черепанов, Е.А. Калинина // Электротехнические системы и комплексы. – 2016. – № 3 (32). – С. 20–26.
7. Barros J. Rapid voltage changes in power system networks and their effect on flicker // IEEE Trans. on Power Delivery. ‒ 2016. Vol. 31, № 1. ‒ P. 262–270. DOI: 10.1109/TPWRD.2015.2452967.
8. Eghtedarpour N. Effective voltage flicker calculation based on multiresolution S-Transform / N. Eghtedarpour, Farjah E, Khayatian A. // IEEE Trans. on Power Delivery. ‒ 2012. – Vol. 27, № 2. ‒ P. 521–530. DOI: 10.1109/TPWRD.2011.2182064.
9. Morsi W.G. Effect of distributed generation on voltage flicker in distributed systems: A case study // Electrical and Computer Engineering CCECE Canadian Conference (May 5–7 2008 Niagara Falls). – 2008. – P. 65–70. DOI: 10.1109/CCECE.2008.4564497.
10. Osipov D.S. Applications of wavelet transform for analysis of electrical transients in power systems: the review // Przeglad Elektrotechniczny. ‒ 2016. Vol. 92, № 4. ‒ P. 162–165. DOI: 10.15199/48.2016.04.35.

Колебания напряжения в узлах нагрузок систем электроснабжения (СЭС) возникают при наличии электроприёмников, имеющих повторно-кратковременные режимы работы с быстрым изменением режима и потребляемой мощности. К таким электроприемникам относятся [4]: дуговые сталеплавильные печи (ДСП), прокатные станы, установки дуговой и контактной электросварки, нагнетательные насосы и др.

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

Колебания напряжения возникают в сетях различных классов напряжения. В работе [7] на конкретных примерах сетей низкого и среднего напряжения Испании авторы приводят методики анализа колебаний напряжения в соответствии с международным стандартом IEC 61000–2–12.

dol1.tiff

Рис. 1. Колебания напряжения

На рис. 1 представлена осциллограмма напряжения, где размах колебания dolgi1.wmf может быть определён из разности максимального dolgi2.wmf и минимального dolgi3.wmf значениями огибающей напряжения:

dolgi4.wmf, (1)

Колебания напряжения в системах электроснабжения следует рассматривать как модулированное колебание синусоидальной величины при изменении амплитуды сигнала. В этом случае простейшее амплитудно-модулированное колебание может быть записано по формуле

dolgi5.wmf, (2)

где Um – амплитудное значение напряжения, В; m – глубина модуляции, В; w – угловая частота питающей сети, рад/с; W – угловая частота огибающей, рад/с.

В то же время колебания напряжения можно представить в виде суммы синусоиды с частотой питающей сети а=50 Гц и двух дополнительных синусоид с частотами близкими, но не равными промышленной частоте [9]:

dolgi6.wmf, (3)

где ω±W называют боковыми частотами.

Детальный анализ причин и последствий возникновения колебаний напряжения в системах ограниченной мощности (при питании от автономных источников энергии) представлен в работе [9].

Теория и моделирование

Каждый электроприёмник, являющийся причиной возникновения колебаний напряжения имеет свой специфический график электрических нагрузок, определяемый технологическими условиями производства. Для решения задач моделирования, расчета и анализа режимов работы систем электроснабжения при условии возникновения колебания напряжения в узлах нагрузок следует учитывать особенности каждого типа электроприёмников. В работе [6] авторами исследуется проблема колебания напряжения в узлах с асинхронными двигателями, работающими с резкопеременной нагрузкой. При этом освещается применение теории вероятности в качестве математического аппарата, поскольку величины потребляемой двигателем активной P(t) и реактивной Q(t) мощностей могут быть учтены как случайные величины. В таком случае характеристиками режима работы узла нагрузок выступают функции математического ожидания M[P(t)], M[Q(t)] и дисперсии D[P(t)] и D[Q(t)]. На базе представленного анализа автор приходит к выводу, что для расчета колебания напряжения, вызванной работой асинхронных двигателей с резкопеременной нагрузкой целесообразно «использовать решетчатые модели графиков электрических нагрузок и напряжения сети, интерполированные сплайн-функциями первой степени» [6].

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

Одним из существующих подходов анализа колебаний напряжения является применение преобразования Стоквелла (S-transform), которое является некоторой модификацией оконного преобразования Фурье (ПФ). Отличие в данном случае заключается в дополнительном масштабировании оконной функции по частоте [7].

В последнее время для решения широкого круга задач цифровой обработки сигналов находит вейвлет преобразование [3].

Вейвлет преобразование (ВП) функции напряжения dolgi7.wmf при изменении параметров масштаба и сдвига по времени может быть реализовано путём свёртки dolgi8.wmf и вейвлет функции dolgi9.wmf [10]:

dolgi10.wmf, (4)

где k и j сдвиг по времени, и масштабирующий коэффициент, соответственно, а dolgi11.wmf нормирующий показатель. Преимущество ВП заключается в отсутствии эффекта «растекания спектра», а также в возможности трёхмерной интерпретации выходных данных (амплитуда-частота-время).

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

Предлагаемый алгоритм. Численный эксперимент

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

1. Исходный сигнал напряжения u(t) возводится в квадрат. При этом получаем сигнал, где представлены различные частотные составляющие:

dolgi12.wmf(5)

2. Производится ДВП сигнала u2(t) с такой глубиной разложения [3], чтобы разделить по частотным составляющим в соответствии с алгоритмом Малла.

3. По детализирующим коэффициентам уровня разложения d1 и d2 определяются интервалы времени, на которых происходит изменение режима – возникновение колебаний напряжения.

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

4. По вейвлет коэффициентам, имеющих значение энергии спектра выше Δε производится обратное ДВП (восстановление) сигнала.

5. На основе быстрого ПФ определяются частоты и амплитуды дополнительных частот, определяющих характер колебания напряжения. При этом ПФ производится только в отрезке времени, полученном в п. 3 настоящего алгоритма.

Проиллюстрируем предложенный алгоритм численным экспериментом. Пусть сигнал напряжения (рис. 1) задан функцией

dolgi13.wmf (6)

где на интервале времени от 0,04 (с) до 0,24 (с) (короткие промежутки времени выбраны для наглядности представляемого алгоритма) будет присутствовать колебание напряжения с двумя группами боковых частот (ω±Ω1)= (50±3) и (ω±Ω1)= (50±8). Глубина модуляции принята равной для обеих частот и составляет

dolgi14.wmf (В).

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

dol2.tiff

Рис. 2. График детализирующего вейвлет коэффициента d1

Характерный всплеск амплитуды высокочастотного вейвлет коэффициента d1 дает информацию о начале отрезка появления колебания напряжения в сети

dolgi15.wmf

и завершения данного отрезка

dolgi16.wmf,

что с достаточной степенью точности соответствует заданной функции (6).

Далее определяются вейвлет коэффициенты, чья энергия спектра превышает заданный пороговый уровень (в относительных величинах к основной частоте). В нашем случае такими коэффициентами определены d7 и d9 (рис. 3, 4).

dol3.tiff

Рис. 3. Вейвлет коэффициент d9

dol4.tiff

Рис. 4. Вейвлет коэффициент d7

По матрицам значений выделенных (имеющих существенный вклад в энергию спектра) вейвлет коэффициентов производится обратное вейвлет преобразование, другими словами – восстановление сигнала определённого частотного диапазона. На завершающем этапе производится быстрое преобразование Фурье (БПФ) восстановленных по вейвлет коэффициентам (рис. 3, 4) частотных компонент. Итогом БПФ становится определение амплитуды и частоты, присутствующих в сигнале боковых частот, которые характеризуют размах колебаний напряжения.

dol5.tiff

Рис. 5. Спектр частот, характеризующих колебания напряжения

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

Выводы

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


Библиографическая ссылка

Долгих Н.Н., Киселёв Б.Ю., Осипов Д.С., Горовой С.А. ИДЕНТИФИКАЦИЯ КОЛЕБАНИЙ НАПРЯЖЕНИЯ В ЭЛЕКТРОЭНЕРГЕТИЧЕСКИХ СИСТЕМАХ НА ОСНОВЕ ПАКЕТНОГО ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЯ // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 12-2. – С. 190-195;
URL: https://applied-research.ru/ru/article/view?id=10804 (дата обращения: 20.10.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074