CC BY
107
7. Клевцов С.И. Прогнозирование изменения состояния параметров технического объекта с помощью интеллектуального микропроцессорного модуля // Всероссийская научнотехническая конференция «Проблемы разработки перспективных микро- и наноэлек-тронных систем (МЭС)». Сборник трудов. - 2010. - № 1. - С. 619-623.
8. Лукашин Ю.П. Адаптивные методы краткосрочного прогнозирования временных рядов.
- М.: Финансы и статистика, 2003. - 416 с.
Статью рекомендовал к опубликованию д.т.н., профессор А.Е. Панич.
Клевцов Сергей Иванович - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: [email protected]; 347900, г. Таганрог, ул. Петровская, 81; тел.: 88634328025; к.т.н.; доцент.
Klevtsov Sergey Ivanovich - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: [email protected]; 81, Petrovsky street, Taganrog , 347900, Russia; phone: +78634328025; cand. of eng. sc.; associate professor.
УДК 621. 315.2/075.8
М.Н. Дубяго
ИССЛЕДОВАНИЕ ХАРАКТЕРИСТИК ЧАСТИЧНЫХ РАЗРЯДОВ В ЭЛЕКТРИЧЕСКОЙ ИЗОЛЯЦИИ СИСТЕМ ЭЛЕКТРОСНАБЖЕНИЯ*
Показана взаимосвязь характеристик частичных разрядов (ЧР) с состоянием электрической изоляции силовых кабельных линий (СКЛ). Рассмотрены вопросы старения и разрушения изоляции ЧР СКЛ в условиях эксплуатации. Отмечено, что регистрация характеристик ЧР позволяет оценить качество изоляции и выявить местные дефекты, которые практически невозможно определить обычными испытаниями высоким напряжением или измерениями каких-либо интегральных характеристик изоляции (тангенс угла диэлектрических потерь, сопротивление изоляции и др.). Показано, что для наблюдения за характеристиками ЧР необходимо, использовать модели MATLAB Simulink, позволяющие выявить максимальный ЧР, количество частичных значений разряда, содержание частичного импульсного разряда и другие параметры. Методика расчета может быть использована для моделирования старения и отказа электрической изоляции. Показано, что исследование появления включения, процесса старения и достижения предельного состояния изоляции необходимо для диагностики и оценки остаточного ресурса изоляции систем электроснабжения.
Электрическая изоляция; частичные разряды; моделирование старения.
M.N. Dubyago
RESEARCH OF CHARACTERISTICS OF PARTIAL DISCHARGE IN ELECTRIC INSULATION OF SYSTEMS OF POWER SUPPLY
The relationship ofpartial discharge (PD) with the state of the electrical insulation ofpower cable lines (PCL). The problems of aging and failure insulation PD PCL in the field. It is noted that the registration of PD to evaluate the quality of insulation and identify local defects that are almost impossible to conventional high voltage tests or measurements of any integral characteristics of the insulation (the dielectric loss tangent, insulation resistance, etc.). It is shown that for monitoring the characteristics of PD is necessary to use MATLAB Simulink models help identify
*
Работа поддержана грантом РФФИ №213-01.-24/2013_96:офи_м_РЖД «Разработка и исследование распределенной системы интеллектуального управления процессом производства, передачи и распределения энергии».
the maximum PD, the number of PD values, the contents of the PD pulse, and more. Method of calculation may be used for the simulation of aging and failure of electrical insulation. It is shown that the study of the appearance of inclusion of the aging process and achieve the ultimate state of isolation is necessary to diagnose and assess the residual life of insulation of electrical systems.
Electrical insulation; partial discharges; modeling aging.
Введение. Анализ современных методов и средств диагностики изоляции под рабочим напряжением показывает высокую эффективность метода, основанного на измерении характеристик ЧР [1, 5]. Взаимосвязь характеристик ЧР с состоянием электрической изоляции известна давно [2, 3]. Разработка более совершенной методики расчета электрической изоляции СКЛ по длительной электрической прочности является сложной задачей. Такая методика основывается на математическом моделировании старения и отказа электрической изоляции (ЭИ). Исследование процесса старения и дости-
жения предельного состояния изоляции СКЛ необходимо для разработки как методов электрического расчета, так и диагностики и оценки остаточного ресурса. Методику расчета составляет подход, основанный на соотношениях, - степенном законе старения
Епт = с ons t, (1)
и двух- или трехметрическом распределении Вейбула
F (Е) = 1-ехр( -|-|) &, (2)
где - вероятность отказа, b - параметр формы, - параметр масштаба, -параметр сдвига распределения.
Выражение (1), справедливо для весьма грубых дефектов типа больших газовых включений. При таких дефектах основным механизмом старения электрической изоляции являются частичные разряды и выражения (1, 2), определяющие методику электрического расчета. Разработка модели старения, учитывающая качество изоляции, может быть получена путем измерения совокупности характеристик ЧР.
Механизм явления ЧР. Пусть к изоляции находящейся между электродами
(рис. 1,а, где Св - емкость включения, Сд - емкость диэлектрика, расположенного последовательно с включением, Са - емкость диэлектрика) приложено синусоидальное напряжение u = Umsinrnt, то при отсутствии ЧР напряжение на емкости включения Св будет синусоидальным, u-R = U-Rmsinwt, где UEX - амплитуда напряжения на включении [4].
Ув m = Уш СД / ( С В + Сд) = Уш?7 - (3)
где п = Сд/(Сд + Св).
Представим включение в виде прослойки толщиной dв, вытянутой поперек силовых линий поля:
JJ Um£q(d-dB) £Дав
в m (£В+_5а_Л ш £В (d - йв)+£ДЙВ ' ( )
\dB d-dj
В момент воздействия на изоляционную структуру первого полупериода напряжения (рис. 2) ЧР возникнет тогда, когда напряжение на включении (емкость
Св, рис. 1) достигнет значения ^.з (напряжения зажигания разряда во включении). При пробое напряжение на включении падает до Цв.п (напряжения погасания разряда во включении), при котором разряд гаснет. После погасания разряда напряжение на включении начинает нарастать от значения Цв.п по кривой, соответствующей изменению приложенного напряжения, смещенной по вертикали на значение постоянной составляющей, возникшей вследствие появления зарядов на
поверхности включения (на емкости Св). Когда напряжение UB на емкости Св достигнет значения ^.з, процесс повторяется.
Отт
Ход |
"Т" Са_|_
СВЕ Т
о ^ I
-ш
б) - в)
Рис. 1. Изоляционная структура (а), эквивалентная схема (б), схема для обнаружения ЧР в изоляции (в)
Следовательно, разряды в рассматриваемой области диэлектрика повторяются через промежутки времени, соответствующие изменению напряжения на Св на ДUв=Uв.3-Uв.п. При прохождении напряжения через максимум включение находится под напряжением, значение которого лежит в интервале между и Ц^п. Затем происходит уменьшение напряжения на включении до 0 и дальнейший рост до -Uв.з. В первом приближении можно принять, что при обеих полярностях пробивное напряжение включения одинаковое. При достижении напряжением на включении значения -^^ происходит пробой включения и процесс продолжается (рис. 1, где Va - напряжение на целой части и V - напряжение на включении).
Рис. 2. Развитие во времени ЧР: а - напряжение и б - ток в полости
При таком механизме явления ЧР должны прекращаться при прохождении напряжения через максимум и вновь возникать, когда напряжение на выводах испытуемого объекта достигнет значения П\, соответствующего изменению напряжения на включении на
и _ Ув-3 ^в п _ ^В1 и 1 — — ■
V V
В дальнейшем, после возобновления ЧР, в рассматриваемой области диэлектрика, по-прежнему, повторяются через промежутки времени, соответствующие изменению напряжения на емкости СВ на ДUв:
Д и = (5)
V
Таким образом, число разрядов во включении за один полупериод
__ _ 2ивт-(иВ.З + иВп) , 1 _ 2(ивт~иВл) _ 2(ит~ип)
тЧ . р = ^в.з - 1/вл +1= А ^ =—йТ-’ (6)
... £/В.П ^ВЗ тг тт
где ип =-----. При с/в . з . ~ ив п . и —— = из = /4 .р . имеем:
2(Увт-Ув.з) _ 2(ит-иЧР) ^
тч . р .----ХГВ = ХБ . (7)
При ив.п << ив.з
т = 2 ^вт = 2 (8)
тч . р . = т,;— = тг—- (8)
иВ.З иЧ.Р.
При этом число разрядов во включении за 1с будет
Пч . р . = 2 /тч. р . , (9)
где f - частота приложенного напряжения.
Учитывая ИИч р. = ИИи — И— = (С в+сд) (^з ^в.п), рчр . = Пч.р.|/[/ч р. и (4), мощность ЧР в одном включении может быть представлена в виде
Рч р . = И4 р Пч р . = 2 / ( св + с д) ( ив . з + ив .п) ( ивт — СВ .п) - (1°)
если , то
Рч .р . = 4/ (с в + с д) /в . з ( /в т — ив . 3 ) = 4/ Сд /в . з ( /т — ич .р .) (11)
если , то
Рч.р . = 2 / (Св + Сд) /в .з ивт = 2 /Сд /в .з /т. (12)
Из (10-12) следует, что для одного включения или при фиксированном количестве одинаковых включений мощность ЧР линейно возрастает с ростом приложенного напряжения. При , учитывая и ,
мощность ЧР может быть также представлена в следующем виде:
Рч.р. = Чч.р. ич.р.Пч.р. = ^ч.р. ич.р.- (13)
соответственно при /в п « /в з, учитывая (1-17) и (1-23),
р _ ЯЧ.Р.иЧ.Р.пЧ.Р. _ ^ЧР.^ЧР. чр. 2 2 . () На основании Рчр . = 1 / (<7чр. х//х + <7чр 2/2 + —+ <?ч.р.тит) ^, мощность Рчр ., выделяемая в единичном включении, может быть также представлена в виде
рч . р . = (=1р. Чч . р . гиг- (15)
где ui - мгновенное значение напряжения на диэлектрике в момент времени разряда с кажущимся зарядом <7ч р . ;; тч р . - общее количество разрядов в единичном включении за полупериод Т/2. При использовании (15) необходимо учитывать знаки при <7ч р . I и иР
Рассмотрим один полупериод (рис. 2) от максимума напряжения одной полярности до максимума напряжения противоположной полярности, и определим энергию WT/2, подводимую к выводам изоляционной структуры (рис. 1,а) за полпериода 7У2, для компенсации энергии, выделившейся во включении вследствие ЧР.
Энергия ЧР, возникающая от момента времени ^ до момента времени 12, соответствующего переходу напряжения через нуль, должна быть взята с одним знаком, а
энергия ЧР, возникающих от момента времени 12 до момента времени 13, соответствующего переходу напряжения через максимум, - с другим.
Если принять <7ч р I = <7ч р = с о П5 £, то формула ■ ■ - х Рч р — ~? \
(Чч.Р.1^1 + Чч.Р.2^2 + I" Чч.Р.т^т)^! может быть представлена в следующем виде:
Рч.р. = 7<?ч.р.{-[ди + 2Ди + -+(т1- 1)ди + (ит - Щ)] + [Ди + 2Ди + ...+т2-1Д и+1/т=2Тдч.р. /т—/1+Д и + /т—/1+Ди + ...+и.т—1У1++т3-1Д и+ /т, (16)
где ; ; .
Формула (16) может быть преобразована к виду
п 77 I лтА ^ч.р./^1 (2^т ^1+Аи) /п\
Рч . р . = Чч .р .т 3 / ( 2 ит — и1 + А и ) = -------—-----------. (17)
Учитывая (16) и (17), имеем
р _ Яч.^.и1пЧ.Р. _ Ы.рРг
ч . р . = 2=2. ( )
При //в з = ив п имеем и1~2ЦчР и Рч р = /ч р «ч р (19)
При ив.п «« ив . зрч. р . = ^ . р . ич . р /2 = 2 /Чч.р . ^, (20)
Мощность, расходуемая во включении (13, 14) и подводимая к изоляционной структуре (19, 20), показывает их идентичность.
Прогнозирование появления включения. С ростом напряжения изоляции количество включений, в которых возникают ЧР, увеличивается, что приводит к более сильной зависимости числа разрядов в секунду и мощности ЧР от напряжения, чем по формулам (16, 19, 20). Увеличение с ростом напряжения числа включений, в которых возникают ЧР, происходит в соответствии со статистическим разбросом напряжений зажигания ЧР в отдельных включениях. Распределение напряжений зажигания ЧР подчиняется нормальному закону со средним значением иВЗср и среднеквадратичным отклонением свз. Таким образом, плотность вероятности появления включений с напряжением зажигания ЧР иВЗ. равна
(^В.З-^В.Зср)
°в.з
(21)
а вероятность появления включения с напряжением зажигания менее ивз составляет
Гив.з I Ов.З-Ь'в.Зср)
7В.З. -° У 2о-|3
Переходя к напряжениям на выводах испытуемого объекта, имеем
иЧР = ивз/т/; Уч.Рср. = ив Зср/т7; <тЧР = оъз/г\\
* («в. 3) = вз е*р ( —) ®в.з. (22)
* (ич.р ) = ч ( — ( " ч '^7 ) В ) Л ич.р (23)
При напряжении, меньшем напряжения начальных ЧР ин, вероятность возникновения ЧР во включении должна быть достаточно мала. В области малых вероятностей, зависимость F(U), а следовательно, и определяемые ею зависимости пч.р (и) и 1ч. р(и) имеют вид
Р(£7) = А^и/ипУ'щ. Р = пч.Р0(и/ин)а1 = А2иа* = А'2Еа'-,
/ч р = /ч .р о ( и / /„) “ 1 = А з« “ 1 = А 3 Е “ 1, (24)
где а1 = 4ии/ачр, пч. ро и 1ч. р0 - число ч. р. и ток ч. р. при напряжении ин.
Учет статистического разброса напряжения зажигания ЧР в отдельных включениях приводит к тому, что мощность ЧР зависит от напряжения более сильно, чем по (10). Так как в этом случае ЧР во включениях, с наибольшим напряжением зажигания, а следовательно, обладающие и наибольшей энергией, возникают вблизи амплитуды приложенного напряжения, то Рч р . * /ч р /т. Учитывая 1чр (24), имеем
Рч.р . = А^ “ ^ 1 = А 4/ “ = Рчр о (// /И) “ = А 4 Е“. (25)
Таким образом, учет статистического разброса напряжения зажигания ЧР в
отдельных включениях приводит к степенной зависимости мощности ЧР от напряжения. Значение кажущегося заряда связано с физическим механизмом раз-
12 -11
вития разрядов. Начальные ЧР с интенсивностью до 10 -10 Кл при длитель-
ном воздействии напряжения вызывают старение изоляции и разрушение. При дальнейшем повышении напряжения в изоляции резко возрастает кажущийся заряд единичного ЧР, обусловленный изменением физики развития ЧР, например, образование дендрита в твердой изоляции. Такие ЧР называются критиче-
и
8 -7
скими и характеризуются кажущимся зарядом ~10- -10 Кл - для аппаратной и кабельной изоляции. Интенсивность определяется конструкцией, структурой изоляции и объемом, в котором развиваются разряды. Критические ЧР приводят к более интенсивному разрушению изоляции и их возникновение резко сокращает срок службы изоляционной конструкции. Эти разряды так же, как и начальные ЧР имеют степенную зависимость характеристик от напряжения (напряженности электрического поля) в соответствии с (25), однако для них показатель степени а имеет значительно большие значения, чем для начальных ЧР, и обычно лежит в пределах 12-16.
Результаты моделирования и обсуждения. Непосредственно в высоковольтном энергетическом оборудовании характеристики ЧР не могут быть измерены. Следовательно для наблюдения за характеристиками ЧР необходимо использовать модели МАТЬАВ Simulink. Модель содержит область с пониженной электрической прочностью (включение). Модель для обнаружения ЧР с использованием МАТЪАВ состоит из трех конденсаторов, рис.1,б. Конденсатор С,; представляет емкость пустоты в исследуемом объекте, Сь - емкость здоровой части, соединенной последовательно с пустотой, Са - емкость здоровой части изолятора. Емкости значения трех конденсаторов могут быть рассчитаны:
Сс = (е0А) / С = 2 . 079 * 1 0-13 Р;
Сь = (е0егА ) / (<2 - С) = 5.42 0 * 1 0- 13Р;
Са = (е0егА)/й = 3 0 . 6 * 1 0 - 1 3 Р.
Параметры, используемые для моделирования: Высоковольтный измерительный конденсатор Ст =1000 пФ, разделительный конденсатор Ск = 1000, мкФ, диэлектрическая проницаемость Е 0 = 8,85 * 10-12 Ф/м, относительная диэлектрическая проницаемость Еа=3,5. Сопротивление R = 50,0м, Индуктивность L = 0,60 мГн, Емкость С=0,45 мкФ. Характерные осциллограммы ЧР приведены на рис. 5 при
5 кВ, где пики - это всплески активности ЧР. Необходимо проанализировать каждый из выявленных частичных разрядов с выделением явно выраженных (задача отделения сигналов ЧР, возникающих в изоляции от внешних и внутренних помех различного происхождения, рис. 1,в, решается способом селекции сигналов ЧР с использованием вейвлет-анализа [6]). После этого определяется расстояние на кабеле до его источника.
Рис. 5. Характерные осциллограммы ЧР при 5кВ
Таким образом, MatLab Simulink модели, необходимы для наблюдения активности ЧР внутри твердой изоляции. Установлено, что с ростом напряжения изоляции количество включений, в которых возникают ЧР, увеличивается. На основе модели Simulink исследованы характеристики ЧР, чтобы выявить максимальный ЧР, количество частичных значений разряда, содержание частичного импульсного разряда и другие параметры ЧР.
Заключение. Получены зависимости характеристик ЧР от напряжения испытания. С увеличением напряжения испытания уровень частичных разрядов в изоляции оборудования повышается. Исправная изоляция имеет значение максимального кажущегося заряда намного меньшей величины даже при установленном испытательном напряжении. Таким образом, взаимосвязь характеристик ЧР с состоянием электрической изоляции не вызывает сомнений. Регистрация импульсов ЧР и их последующий анализ позволяет не только осуществлять постоянный мониторинг состояния высоковольтной изоляции, но и, при определенных условиях, выявлять характер развивающегося дефекта и место его расположения.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК:
1. Привалов И.Н. Современные методы и технические средства для испытаний и диагностики силовых кабельных линий номинальным напряжением до 35 кВ включительно: Учеб. пособие, 2008.
2. Nagasaki S., Yoshiba N., Aihara M. Philosophy of design and experience on high voltage XLPE cables and accessories in Japan // CIGRE. - 1988. - № 21-01. - P. 8.
3. Sekki Y. Recent technical progress in extrahigh-voltage XLPE cables // Hitachi Cable Review.
- August 1993. - № 12. - P. 3-10.
4. Кучинский Г.С. Частичные разряды в высоковольтных конструкциях. - Л.: Энергия. Ленингр. отд-ние, 1979. - 224 с.
5. Дубяго М.Н., Полуянович Н.К. Метод амплитудного и фазового распределения импульсов частичных разрядов в задачах исследовании изоляции кабельных линий // Известия ЮФУ. Технические науки. - 2012. - № 7 (132). - C. 200-205.
6. Дубяго М.Н., Полуянович Н.К. Метод селекции сигналов ЧР с помощью Вейвлет преобразований // Известия ЮФУ. Технические науки. - 2013. - № 2 (139). - C. 99-104.
Статью рекомендовал к опубликованию д.т.н., профессор В.Э. Бурлакова.
Дубяго Марина Николаевна - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: [email protected]; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: 89508519304; кафедра электротехники и мехатроники.
Dubyago Marina Nikolaevna - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: [email protected]; 44, Nekrasovsky, Taganrog, 347928, Russia; phone: +79508519304; the department of electric technics and mechatronics.
