автоматизация
Реализация синхронного включения при БАВР
В статье рассмотрена проблема реализации синхронного включения секционного выключателя (СВ) в комплексе быстродействующего автоматического ввода резерва. Для наглядности приведены экспериментальные результаты и описание алгоритма для точного определения момента выдачи команды на включение секционного выключателя.
Одно из требований к быстродействующему автоматическому вводу резерва (далее — БАВР) — включение секционного выключателя в момент, когда вектора напряжения на «здоровой» секции шин и ЭДС обобщенного двигателя, выбегающего на секцию шин, потерявшую питание, имеют небольшое расхождение по фазе (в [1] указан диапазон до 30°). Это позволяет минимизировать значение переходного тока при включении, однако методические погрешности в определении угла и частоты в цифровых блоках и игнорирование скорости изменения частоты могут привести к нарушению этого требования.
Рассмотрим режим потери питания на секции шин №1 (СШ1), к которой подключен синхронный двигатель (далее — СД). Примем, что напряжение на «здоровой» секции шин №2 (СШ2) неизменно по модулю и частоте (рис. 1).
Уравнение движения вектора напряжения СШ2 имеет следующий вид:
<РСШ2 = шном ' Ь + <РоСШ2, (1)
где фсш2 - текущая фаза, рад; юном - номинальная круговая частота (угловая скорость), рад/с;
Фо.
- начальная фаза, рад,
и представляет собой равномерное движение по окружности.
Закон движения вектора напряжения СШ1 зависит от закона выбега двигателя (см. ниже). При потере питания напряжение на СШ1 индуцируется только полем возбуждения двигателя. В этом случае частота и фаза напряжения однозначно связаны с частотой и фазой вращения ротора.
В общем виде закон выбега двигателя есть второй закон Ньютона для вращающихся тел:
(2)
где J - момент инерции вращающихся масс, кгм2; Ф - угол поворота ротора относительно заранее выбранной неподвижной оси, рад; М-суммарныймоментсил,действующийнатело,Нм.
При двукратном интегрировании уравнения (2), в процессе нахождения зависимости Ф() возникает две константы, определяемые начальными условиями:
• Ф0 - начальный угол (угол в начальный момент времени), рад;
• ^фМ)1=0=ю0 - начальная угловая скорость, рад/с.
Для выдачи команды на включение СВ необходимо спрогнозировать значение угла в момент времени включения выключателя, что требует знать величину угла между напряжениями — ф, значение угловой скорости — ю и закон ее изменения во времени.
Проанализируем каждый из этих параметров.
Частота
Расчет частоты в устройствах РЗА может реализовываться различными способами [2, 3]. Примем, что алгоритм без погрешности рассчитывает мгновенную частоту, но имеет некоторое окно усреднения для повышения надежности результата (стандарт [4] регламентирует расчет частоты не более, чем за 60 мс).
Это приводит к тому, что при усреднении значений за время Т , при линейно меняющейся частоте, алгоритм будет давать погрешность, равную:
(3)
где / — скорость изменения частоты, Гц/с; Туср — время усреднения, с.
Использование частоты с погрешностью (3) при оценке угла, на который повернется вектор напряжения за время включения выключателя, приведет к погрешности по углу, равной:
Д= 2л ■ Д/ • ГВ1
(4)
При этом действительный вектор успеет описать большую дугу, чем это следует из расчета.
-от погрешности по частоте —от погрешности по углу -от неучетэ линейного изм частоты —суммарная погрешность
Твкл = 22 мс
Начальное ускорение е0, рад/сг
Рис. 2 — Зависимости угловых погрешностей от начального ускорения при использовании быстродействующего выключателя
78
Экспозиция НЕфть газ апрель 2 (62) 2018
Угол
Ортогональные составляющие будем выделять с помощью разложения сигнала в ряд Фурье:
С = Ст1У ■ е->шЧг. (5)
Так как частота входного сигнала отличается от частоты, используемой для вычисления ортогональных составляющих (см. п. Частота), функция под интегралом равна:
у = А- 5т((&) + Дсо)1 + <р). (6)
Интуитивно можно ожидать что вектор уже не будет неподвижен, а начнет совершать вращение с разностной частотой Аю.
Точное выражение после взятия интеграла имеет вид:
С = А-зШ- [1вЛД«"ч»-*> -
__—еЧ((2о+Л <и)Г+ф-А)1 (7)
27Г+Я J'
где Х=п (Аю/ю).
Предположение оказалось верным, вектор имеет низкочастотную составляющую Аю. Однако, в дополнение к этому (в соответствии с АЧХ ДПФ), изменилась амплитуда вектора (в sinX/X раз) и добавилась высокочастотная составляющая 2ю+Аю.
Для темы статьи важно, что вектор в среднем всегда отстает от «реального» вектора на угол X:
Д <р2
п-
Аш
(8)
То есть, при выполнении прогноза действительный вектор уже находится чуть «дальше», чем это предполагает блок.
Закон выбега
Выбег двигателя описывается уравнением [5, 6, 7]:
+ ,)-йП (9)
где кзг — коэффициент загрузки двигателя по активной мощности;
т — постоянная инерции системы «двигатель - приводной агрегат», с;
т0' — приведенный начальный момент сопротивления, Ы/(кг-Ы м);
& — приведенная угловая скорость, ю/ю .
Темп выбега определяется начальным ускорением (dw/dt)t0=-k Jt=e0. При реальных начальных ускорениях от 40 до 400 рад/с2 [5] отклонение закона изменения угловой скорости (частоты) от линейного за время расчетов и усреднений в блоках РЗА мало. Но также может использоваться предположение о том, что за время включения выключателя частота не меняется. В этом случае погрешность по углу равна:
Д<р3
2 '
(10)
где в — угловое ускорение (производная частоты), рад/с2.
Все погрешности имеют одинаковый знак и могут привести к тому, что действительная фаза вектора в момент включения СВ будет отличаться от фазы напряжения «здоровой» СШ больше, чем допускается.
На рис. 2 приведены погрешности по углу в зависимости от значения начального ускорения при использовании быстродействующего выключателя 1БМ15_ЗИе11_РТ2 с электронным модулем управления СМ_1501_01(4) с временем включения 22 мс [8]. На рис. 3 — то же при использовании небыстродействующего выключателя с временем включения 50 мс.
Постоянная инерции для газоперекачивающих агрегатов с электродвигателями СТД-4000-2 и нагнетателем 280-12-7 равна 3,25 с [9], что позволяет вычислить начальное ускорение при выбеге — 96,6 рад/с2. По графикам видно, что только методическая погрешность будет составлять 7 и 18 градусов соответственно. Учитывая цикличность расчетов блоках, погрешность по углу будет еще увеличена, так как «правильный» момент выдачи команды может попасть на середину программного цикла. При больших значениях начального ускорения и выключателях, включающихся за 70—100 мс, методическая погрешность будетсоставлять десятки градусов!
В устройстве БмРЗ-БАВР производства ООО «НТЦ «механотроника» применяются следующие решения для точного определения момента выдачи команды на включение секционного выключателя: 1. Устройство анализирует не только скорость, но и ускорение вращения векторов секций шин относительно друг друга.
2. Учитывается линейное изменение частоты за время включения выключателя, так как погрешность без учета этого факта зависит квадратично от времени включения.
3. Используется адаптивное изменение порядка фильтра для более быстрой реакции на изменение режима сети. Максимальная величина окна наблюдения равна половине периода, что позволяет не учитывать погрешности по углу при расчете ортогональных составляющих.
4. Предпочтительнее использование быстродействующих выключателей, так как точность прогноза на длительный промежуток времени снижается.
Список литературы
1. Устройства быстродействующего автоматического ввода резерва с тиристорным коммутатором напряжением более 1000 В. Общие технические требования. ПАО «Транснефть», 2017.
2. Иванов Н.Г., Солдатов А.В., Наумов В.А., Антонов В.И. Оценка частоты сети в цифровых системах РЗА по переходу через нуль: характеристики точности // Релейная защита и автоматизация. 2013. №4. С. 22-25.
3. Антонов В.И., Наумов В.А., Шевцов В.М. Оценка частоты электрической сети: теоретические основы и практические алгоритмы. Цифровая электротехника: проблемы и достижения.
Выпуск 1. Чебоксары: РИЦ СРЗАУ, 2012. 120 с.
4. Релейная защита и автоматика. Автоматическое противоаварийное управление режимами энергосистем. Микропроцессорные устройства автоматической частотной разгрузки. Нормы и требования. М.: Системный оператор Единой энергетической системы, 2016. 19 с.
5. Крышнев Ю.В. Исследование динамики выбега синхронных двигателей с учетом характеристик приводных механизмов // Вестник Гомельского государственного технического университета
им. П.О. Сухого. 2002. №3-4. С. 74-81
6. Сыромятников И.А. Режимы работы асинхронных и синхронных двигателей. М.: Энергоатомиздат, 1984. 240 с.
7. Курганов В.В., Крышнев Ю.В., Верига Б.А. Учет механической характеристики обобщенного привода в алгоритме быстродействующего самозапуска высоковольтных синхронных электродвигателей // Энергетика. Известия высших учебных заведений и энергетических объединений. 2007. №2. С. 17-33.
8. Комплекс быстродействующего АВР БАВР10_ШеллП"2_1. Техническая информация. Таврида электрик, 50 с.
9. Беляев А.В. Автоматика и защита на подстанциях с синхронными и частотно регулируемыми электродвигателями большой мощности. СПб: ПЭИПК, 2012. Ч.1. 72 с.
Uli
Рис. 3 — Зависимости угловых погрешностей от начального ускорения при использовании небыстродействующего выключателя
МЕХАНОТРОНИКА
интеллектуальные устройства релейной защиты
ООО «НТЦ «механотроника» г. Санкт-Петербург тел. 8-800-250-63-60 www.mtrele.ru