CC BY
28
ПРИМЕНЕНИЕ АЛГОРИТМОВ ВЗВЕШИВАНИЯ ПРИ ЦИФРОВОЙ ОБРАБОТКЕ СИГНАЛОВ РЕЛЕЙНОЙ ЗАЩИТЫ
КУЛИКОВ А.Л., канд. техн. наук
Рассматриваются варианты взвешенной обработки сигналов релейной защиты с применением «окон» и обосновывается необходимость ее применения.
Ключевые слова: релейная защита, цифровая обработка сигналов, окно данных, алгоритм взвешивания, дискретизация непрерывной функции.
WEIGHTING ALGORITHM APPLICATION IN RELAY PROTECTION DIGITAL SIGNAL PROCESSING
A.L. KULIKOV, Ph.D.
The work represents the variants of relay protection weighted signal processing, using “windows”. The author tries to justify the necessity of its application.
Key words: relay protection, digital signal processing, data window, weighting algorithm, continuous function discretization.
В целях уменьшения ошибок аналогоцифрового представления сигналов в задачах релейной защиты применяют специальные адаптивные методы [1, 2]. К ним относятся методы интерполяции [1], сглаживания [1], специального Фурье-анализа [1, 2, 6], наименьших квадратов [1, 2, 6], Прони [2, 6] и ряд других [1].
Однако в известной литературе по цифровой релейной защите не рассматривались преимущества методов взвешенной обработки с применением «окон» [3-6].
Отметим, что имеющуюся конечную последовательность дискретных данных можно рассматривать как некоторую часть соответст-
Хт(п)
s(t)1 ‘1
U
1l_
вующей бесконечной последовательности, видимую через применяемое окно. Например, последовательность наблюдаемых данных хт(п) из N отсчетов (рис. 1,г) можно представить как произведение прямоугольной функции единичной амплитуды, или прямоугольного окна S(n) (рис. 1,в), и бесконечной последовательности хд(п) (рис. 1,б), причем
S(t) Л1' I'1 s T/2'
1°, |t| > T/2,
где Т - интервал наблюдения.
gx(f)
fo
9а(<)
-1 -fo 0
T
fo 1
2
T
Рис. 1. Дискретное преобразование синусоидальной функции
f
0
n
При таком подходе принимается допущение о том, что все ненаблюдаемые (находящиеся за пределами интервала наблюдения) отсчеты равны нулю независимо от того, так ли это на самом деле или нет. Дискретное преобразование Фурье (ДПФ) дт(0 взвешенной окном последовательности хт(п) равно свертке ДПФ последовательности хц(п) и ДПФ прямоугольного окна Б(п).
Дискретное преобразование Фурье наблюдаемой конечной последовательности хт(п) является искаженной версией ДПФ бесконечной последовательности хд(п). Анализ влияния прямоугольного окна на дискретную синусоиду с частотой 10 (рис. 1,д-ж) показывает, что острые спектральные пики ДПФ дд(Г бесконечной последовательности данных хд(п) расширились за счет воздействия ДПФ окна. Минимальная ширина спектральных пиков дт(1) ограничена шириной, определяемой главным лепестком ДПФ этого окна, и не зависит от исходных данных. Боковые лепестки ДПФ окна будут изменять амплитуды соседних спектральных пиков. Это явление называется «просачиванием» [3-6]. Поскольку ДПФ имеет периодический характер, то наложение боковых лепестков соседних спектральных периодов может привести к дополнительному смещению пиков в спектре. Просачивание приводит не только к появлению амплитудных искажений в спектрах дискретных сигналов, но и
может маскировать слабые сигналы на фоне сильных, поскольку боковые лепестки более сильного сигнала могут почти полностью подавить главный лепесток более слабого сигнала. Такая ситуация может, например возникнуть при дискретном спектральном анализе переходных процессов. Аналогичные искажения будут наблюдаться при изучении и других несинусоидальных сигналов.
Обработка с помощью окна используется для уменьшения эффектов, обусловленных наличием боковых лепестков в спектральных преобразованиях. «Окно» (весовая функция) может применяться как во временной области (окно данных), так и в частотной области (спектральное окно). Существуют различные варианты введения окна при цифровой обработке сигналов в устройствах релейной защиты (рис. 2).
Снижение уровня боковых лепестков возможно за счет перехода от прямоугольного окна Б© к окну другой формы. Однако это дается ценой расширения главного лепестка спектра окна, что, естественно, приводит к ухудшению разрешения. Следовательно, вид окна должен выбираться на основе компромисса между шириной главного лепестка и уровнем подавления боковых лепестков. Приведем характеристики наиболее употребительных окон (см. таблицу) и их аналитические описания во временной и частотной областях.
Рис. 2. Варианты применения весовой обработки
Параметры взвешивающих окон
Вид окна Формула Ь||_, дБ Рц Р 0,71 Р 0,99 V
Дирихле (1) -13,3 0,82 0,44 0,07 1/!
Барлета (2) -26,6 1,63 0,64 0,11 1/Г2
Парзена (3) -58,1 3,26 0,91 0,13 1/^
Ханна (4) -31,5 1,88 0,71 0,12 1/Г3
Хемминга (5) -42,2 1,92 0,65 0,11 1/!
Блекмана (6) -58,2 2,83 0,82 0,14 1/Г3
Блекмана-Хэрриса (7) -90,1 3,94 0,94 0,16 1/Г3
Плосковершинное (8) -81,1 4,84 1,86 0,75 1/Г
Рисса (9) -21,3 1,28 0,58 0,1 1/Г2
Коши, а = 3 (10) -30,8 3,25 0,67 0,11 1/Г
Римана (11) -26,5 1,51 0,62 0,1 1/Г2
Гаусса, а = 3 (12) -55,9 3,34 0,8 0,13 1Г
Примечание: - уровень боковых лепестков, дБ; Рь - приведенная ширина основного лепестка, равная произведению шири-
ны основного лепестка АЧХ на длительность окна; Р0,71 - ширина основного лепестка по уровню 0,71; Р0,99 - ширина основного лепестка по уровню 0,99; V - скорость спада амплитуды боковых лепестков с ростом частоты.
Все ниже перечисленные окна - непрерывные функции времени:
S(t) Л1’ ^ < ^2’ gs(f) = S'n(nft);
'1°, Ш > T/2; s nft
S(t) =
S(t) =
1 , It < t/2,
T/2 ' gs(f) =
o, |t| > t/2;
2 t3
S'n(n f T/2) n f T/2
1 - 24
48
|< T/4,
2 (1 - 2 o, Itl > t/2;
3, TI4 < It < T/2,
gs(f) =
S'n(n f T/4)
n fT/4
S(t) = 0,5 [1 + Cos(2n t/T)], |t < T/2;
g(f) J Sin(nfT) 1
|_ n fT 1 - (fT)2 S(t) = 0,54 + 0,46 Cos(2n VT), |t| < T/2,
gs(f) =
S'n(n fT) n fT
1
0,46 (fT)2 0,54 1 - (fT)2
S(t) = 0,42 + 0,5 Cos(2n VT) + 0,08 Cos(4n VT), |t| < T/2,
gs(f) =
S'n(n fT) n fT
0,5 (fT)2 0,08 (fT)2
0,42 1 - (fT)2 0,42 4 - (fT)2
S(t) = 0,3587 + 0,4883 Cos(2nVT) +
+0,1413 Cos(4n VT) + 0,017 Cos(6n VT), It < T/2;
gs(f) =
S'n(n fT)
1
0,4883 (fT)2
0,3587 1 - (fT)2
0,1413 (fT)2 0,0117 (fT)2
0,3587 4 - (fT)2 0,3587 9 - (fT)2
S(t) = 1 +1,9333 Cos(2nVT) +1,286 Cos(4nVT),
It < t/2;
S'n(n fT)
gs(f) = -
1
1,93 (fT)2 - 1,29 (fT)2 1 - (fT)2 4 - (fT)2
0,388 (fT)2 - 0,0322 (fT)2 9 - (fT)2 16 - (fT)2
S(t) = 1 - (2t/T)2, |t < t/2;
s(n)
0,5-
0
S(t) =
S(t) =
1
1 + (2aVT)2 ; S'n(2rct/T)
It < T2; |< T2;
2rct/T
S(t) = exp{-2(aVT)2}, |t < T/2.
Переход к дискретному (или ступенчатому) аналогу окна мало влияет на его параметры.
Весовые коэффициенты S(n) дискретной (или импульсной) весовой функции могут быть определены из соотношения S(n) = S(tn).
При определении моментов времени tn, как правило, учитывается два обстоятельства.
Во-первых, непрерывное окно определено на интервале от -Т/2 до Т/2, в то время как номер дискреты меняется от 0 до N-1.
Во-вторых, возможны два способа дискретизации непрерывной функции. При первом способе промежуток времени Т делится на N равных интервалов и берутся отсчеты непрерывной функции, соответствующие серединам этих интервалов. При такой дискретизации
S(n) = S(tn) = S | T 2n - N +1
2 N
При втором способе промежуток Т, на котором задана непрерывная функция, делится на (Ы - 1) равных интервалов и находятся значения функции на границах. При этом амплитуды крайних отсчетов, для которых п = 0 и п = N - 1, уменьшаются вдвое. В итоге -1',
при п = 0, п = N -1,
S(n) =
1 еf T 2n -N-
2 I 2 N -1
T 2n - N +1 2 N -1
при 0 < n < N -1.
В качестве примера на рис. 3 изображены окна данных Ханна (рис. 3,а) и Хемминга (рис. 3,б).
а)
N-1
Рис. 3. Окна данных: а -Ханна; б -Хемминга © ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»
S
Проиллюстрируем применение весовой обработки с помощью окна на примере исходных данных [2]. Так, в приборах определения места повреждения (ОМП) типа ИМФ используется ортогональное разложение для цифровой обработки дискретных отсчетов фазного напряжения:
N
иа = М £8,и(! -Д д),
'=1
N
ир = М £ Ь и(! -Д д),
!=1
где иа и ир - дискретные квадратурные составляющие; Д1д - шаг дискретизации; I - номер отсчета; N - число отсчетов за период; М - масштабный множитель; а| и Ь| - коэффициенты, соответствующие функциям э1п(ю1) и ооб(ю1) в моменты 1-Д1д.
С учетом выбора в приборах ИМФ N = 12 и принятого масштабирования коэффициенты в' и Ь' определяются множествами:
А = {в'} = {0, 4, 7, 8, 7, 4, 0, -4, -7, -8, -7, -4};
В = {Ь'} = {8, 7, 4, 0, -4, -7, -8, -7, -4, 0, 4, 7}.
Анализ амплитудно-частотной характеристики (АЧХ) такого ортогонального разложения (рис. 4) [2] показывает, что АЧХ соответствует цифровому фильтру в диапазоне 0-300 Гц с шириной полосы пропускания на уровне 0,7, оцениваемой 20-80 Гц. Боковые лепестки АЧХ, имеющей вид |Б1пх/х|, позволяют пропускать нежелательные частоты, искажающие комплексные дискретные оценки напряжения (тока).
Операция взвешивания с использованием окон данных, например, Ханна (рис. 5,а) или Хемминга (рис. 5,б), позволяет улучшить фильтрующие свойства ортогонального разложения особенно при наличии интенсивной апериодической составляющей. Эффект улучшения фильтрующих свойств можно оценить по результирующим АЧХ (рис. 5,а,б), которые свидетельствуют о существенном подавлении нежелательных боковых лепестков, либо по количественным оценкам (см. таблицу).
Увых
а)
^ Гц
^ Гц
б)
Рис. 4. Амплитудно-частотная характеристика ортогонального разложения: а - разложение за 1 период промышленной частоты; б - разложение за 2 периода промышленной частоты
Рис. 5. Результирующая АЧХ совместной фильтрации с использованием ортогонального разложения и окна данных: а - окна Хана; б - окна Хемминга
Таким образом, использование алгоритмов взвешивания во временной или частотной областях позволяет существенно повысить точность оценки фазных составляющих при цифровой обработке сигналов релейной защиты. Анализ состава операций алгоритмов взвешивания позволяет сделать вывод о незначительном увеличении сложности результирующих совокупных вычислений.
Список литературы
1. Johns A.T., Salman S.K. Digital Protection for Power Systems - IEE Power Series 15, Peter Peregrims Ltd., 1995.
2. Аржанников Е.А., Лукоянов В.Ю., Мисриха-нов М.Ш. Определение места короткого замыкания на высоковольтных линиях электропередачи / Под ред. В.А. Шуина. -М.: Энергоатомиздат, 2003.
3. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. - М.: Мир, 1978.
4. Каппелини А.В. и др. Цифровые фильтры и их применение: Пер. с англ. - М.: Энергоатомиздат, 1983.
5. Радиоэлектронные системы: основы построения и теория. Справочник / Под ред. Я.Д. Ширмана. - М.: ЗАО «МАКВИС», 1998.
6. Марпл-мл. С.Л. Цифровой спектральный анализ и его приложения: Пер. с англ. - М.: Мир, 1990.
Куликов Александр Леонидович,
Филиал ОАО «ФСК ЕЭС» - Нижегородское ПМЭС, кандидат технических наук, директор, телефон (831) 257-85-50, e-mail: info@npmes.elektra.ru
