CC BY
184
УДК 648.3.048.1
М.А. Поляков, В.В. Василевский
ОЦЕНКА ОСТАТОЧНОГО РЕСУРСА ИЗОЛЯЦИИ НА ОСНОВЕ УЧЕТА ИНДИВИДУАЛЬНЫХ ОСОБЕННОСТЕЙ ЖИЗНЕННОГО ЦИКЛА СИЛОВОГО ТРАНСФОРМАТОРА
Запропонована методика оцінки залишкового ресурсу ізоляції силового трансформатору шляхом імітаційного моделювання прогнозованих результатів впливів зовнішніх факторів на стан ізоляції.
Предложена методика оценки остаточного ресурса изоляции силового трансформатора путем имитационного моделирования прогнозируемых результатов воздействий внешних факторов на состояние изоляции.
ВВЕДЕНИЕ Силовой трансформатор (СТ) - важнейший элемент систем электроснабжения. Продление фазы эксплуатации его жизненного цикла (ЖЦ) всегда было актуальной задачей науки и практики. Оценка такой возможности производится на основе анализа остаточных ресурсов СТ, среди которых определяющим является ресурс изоляции активной части обмоток (далее "изоляции"). Стандартами в области "надежность в технике" ресурс определяется как временная характеристика объекта - время, за которое техническое состояние объекта с определенной степенью уверенности достигнет своего критического значения, при котором дальнейшая эксплуатация объекта недопустима. Критическое состояние изоляции связывается с определенной степенью полимеризации целлюлозы - материала изоляции.
АНАЛИЗ ИССЛЕДОВАНИЙ Расход ресурса изоляции характеризуется интегралом старения [1]:
1 ?п +Т ,
х0 + т) - т |
к
\ 2,05
к
БАЗ
W
л
1.493
(1)
Wl
БАЗ
• ехр[1п2 •
Фи - 98) д
]йХ,
где Ь(^п+Т) - износ изоляции в интервале времени, который начинается в момент ^ и имеет длительность Т; К, КБАЗ - текущий и базовый показатель кислотности масла, соответственно; W, WБАЗ - текущий и базовый показатель влагосодержания масла, соответственно; 0. -температура наиболее нагретой точки изоляции; Д -интервал температур 0*, на котором происходит, при прочих равных условиях, удвоение износа изоляции, Д = 6 °С (в соответствии с рекомендациями [2]) или 7 °С [3]. Таким образом, выражение (1) определяет износ изоляции как функцию показателей К, W, 0. узлов СТ, значения которых в свою очередь зависят от времени.
В работах [1, 3] утверждается, что расход ресурса изоляции это случайный процесс, а остаточный ресурс - случайная величина. В соответствии с этим подходом на основании данных о расходе ресурса по нескольким годам работы СТ средствами математической статистики вычисляется у - процентный износ. Вместе с тем, в известных работах не раскрыты механизмы и характер случайности расхода ресурса, его взаимосвязь с факторами эксплуатации СТ, что снижает достоверность оценки у - процентного износа.
ЗАДАЧА ИССЛЕДОВАНИЯ Задача работы состоит в разработке методики оценки остаточного ресурса изоляции СТ как случайной величины, включая:
• выбор/разработку моделей зависимостей значений показателей входящих в "интеграл старения" от параметров эксплуатации и технического состояния СТ;
• выбор/разработку методики прогнозирования параметров эксплуатации технического состояния СТ;
• определение путем компьютерного моделирования значений показателей СТ во временном интервале прогноза износа изоляции и оценку статистических результатов расчета по этим данным "интеграла старения", включая у - процентный износ.
ОСНОВНАЯ ЧАСТЬ Фазы жизненного цикла СТ и их роль в формировании и расходовании ресурса изоляции приведены на рис. 1.
С позиций поставленной задачи оценки остаточного ресурса изоляции, фазу проектирования, включая проектирование технологических процессов (ТП) изготовления, охарактеризуем как процесс принятия решений, относительно конструкции и технологии изготовления, ограничивающих потенциальный Я" ресурс. А также как процесс (конечно непреднамеренного) внесения скрытых конструктивных и технологических дефектов в СТ.
Так, принимая решение о выборе материала изоляции обмоток, проектант фактически задает параметры относительной скорости расхода ресурса этой изоляции вследствие термического износа [2]:
• для обычной бумаги:
(6и-98)
V - 2 д ,
• для термически улучшенной бумаги: 15000 15000
(2)
„ , 110+273 9. +273 , ,,ч
V = ек й ' . (3)
Входными данными для проектирования (ВДП) являются цели проектирования, требования, представления об идеальных решениях и решениях-аналогах, а также качество персонала, включая его профессионализм и ответственность. Результатами проектирования являются комплект конструкторской (КД) и технологической (ТД) документации на СТ, а также модель потенциального ресурса Я", с параметрами, которые отражают принятые решения относительно конструкции и технологии изготовления.
© М.А. Поляков, В.В. Василевский
X
X
0
X
Рис. 1. Фазы жизненного цикла СТ
В фазу изготовления включим, кроме собственно изготовления в условиях завода, транспортировку и монтаж изделия на объекте заказчика. Эта фаза характеризуется внесением явных и скрытых дефектов в ходе технологических (формообразующих, сборочномонтажных и др.) операций ТП и коррекцию качества изделия в ходе контрольных операций.
Входными данными для оценки ресурса на стадии изготовления являются качество КД, ТД, материалов, комплектующих, оборудования и персонала (МКОП). Результатом изготовления является изделие с реальным Яр ресурсом относительно некоторого диапазона ожидаемых условий эксплуатации. Техническое состояние СТ определяется вероятным наличием в нем скрытых дефектов, определяемых качеством КД, ТД и МКОП. Вследствие чего, как правило, Яр < Я". Модель Яр опишем как нечеткую функцию ресурса зависящую от переменных отклонений входных факторов фазы изготовления от нормы. Разработка этой модели выходит за рамки данной работы и будет рассмотрена отдельно.
В фазе эксплуатации происходят изменения технического состояния (ТС) и ресурса СТ. Входными данными для учета этих изменений в подфазе, собственно, эксплуатации являются данные мониторинга параметров технического состояния, дестабилизирующих факторов эксплуатации и изменений нагрузки СТ (ФЭиН).
Результатами процесса эксплуатации, кроме эффектов у потребителей электроэнергии, являются расход ресурса, выявленные дефекты и отклонения от нормы параметров ТС а также параметры оперативной (среднесрочной) модели прогноза (ОП) динамики факторов внешней среды и нагрузки.
В подфазе технического обслуживания (ТО) происходит устранение выявленных дефектов и восстановление ТС изоляции, например путем регенерации или замены масла. Эффективность этой работы зависит от моментов времени проведения и качества ТО, включая качество МКОП. В результате ТО в экс-
плуатацию возвращается СТ, который может находиться в эксплуатации больший срок, чем до ТО. Оценка этого срока должна проводиться на основе ОП и СП значений параметров внешних воздействий на СТ, моделей параметров СТ как функций этих воздействий и моделей расхода ресурса. Параметры ОП должны учитывать характер "быстрых" изменений внешних факторов, длительности которых соизмеримы со значениями тепловых постоянных обмотки и масла, а также параметры суточных, недельных, сезонных циклов температуры, влажности и нагрузки.
Горизонт СП соизмерим со сроками продления эксплуатации СТ и может превышать десять лет. Примером параметра ОП может служить амплитуда бросков нагрузки СТ, а СП - динамика среднесуточного, среднегодового тока нагрузки конкретного СТ в ближайшем десятилетии. При оценке ресурса Я° должна учитываться также запланированная периодичность ТО и ее влияние на износ. Модель Я° базируется на моделях реального ресурса Яр и модели расхода Ь. Модель Я° должна учитывать вероятностный характер ОП и СП.
Фаза эксплуатации переходит в фазу утилизации при возникновении финальной аварии СТ. Эта фаза используется при построении модели финального ресурса ЯФ.
Рассмотрим более подробно процесс определения ресурса Я°. С этой целью определим следующие моменты времени: €° - завершения изготовления; €' -последнего экспериментального определения ресурса Я°; 1 - прогнозирования; !х+Тг - прогнозируемого с у
- процентной степенью уверенности завершения эксплуатации.
Задача определения в момент 1 ресурса Я° сводится к нахождению V, т.е. Я°= V, для которого:
(4)
или
Яр - Ьр (1р°,1х) = Ьп (1х,1х + Т1)
Яр> - Ьр (Гр<, Iх) = Ьп (1х,Iх + Т1)
(5)
где Ьр(1р° ,1х), Ьр(1р' ,Г) - измеренный ("реальный")
расход ресурса на интервалах [1р° ,1х] и [1р1 ,1х], соответственно; Ьп(Гс, 1х + V) - прогнозируемый расход ресурса на интервале [1х, f + V].
Под измерением расхода будем понимать его расчет по формуле (1), куда подставляются измеренные в процессе мониторинга параметров реального СТ значения 0/(0, Ж“(1), К“ (1), 1е[1р°, 1х] или 1е[1рг, 1е].
Под прогнозированием расхода будем понимать его расчет по формуле (1), по прогнозируемым с помощью компьютерной модели значениям 0/(0, Ж"(1), К"(0, [1х, 1х + V] и параметров технического состоя-
ния СТ.
Компьютерные модели К, Ж, 0* могут быть построены на основе уравнений соответствующих балансов: теплового - для 0*; баланса влажности масла -для Ж; баланса кислотосодержания - для К. Уравнения соответствующих балансов - это дифференциальные уравнения, позволяющие определить изменение параметра, влияющего на расход ресурса изоляции, в зависимости от значений факторов эксплуатации. Параметры этих уравнений могут быть определены через конструктивные характеристики СТ, результаты испытаний или мониторинга параметров. Так, известны термодинамические модели СТ [2, 5, 6], позволяющие при известных тепловых параметрах СТ определить изменение 0* как реакцию на уровень тока нагрузки I, и температуру окружающей среды 0а В работах [7, 8] приведены примеры моделирования изменений технического состояния элементов СТ.
Решение задачи прогнозирования расхода ресурса связано с последовательным решением этих уравнений для всех значений 1,0) и ва(0 в интервале времени 1е [1х, 1х + V] с интервалом изменения времени Д. Неопределенность прогнозов нагрузки и температуры окружающей среды заставляет рассматривать I, (1) и 0а(1) как случайные функции прогноза (СФП). Из-за этого уравнения балансов приобретают свойства стохастических дифференциальных уравнений и могут быть решены относительно некоторого количества реализаций СФП с оценкой вероятности появления этих реализаций. Совокупность этих оценок формирует функцию вероятностей реализаций СФП (ФВР СФП).
Таким образом, каждой сочетанию реализаций СФП I/, 0/ имеющему вероятность появления Р(1/, вау), соответствует полученное в результате моделирования 0* и рассчитанный по формуле (1) износ Ь. Следовательно, вероятность износа - это вероятность появления соответствующего сочетания реализации СФП. Совокупность вероятностей износа для различных сочетаний реализаций СФП I, и 0а задает функцию вероятностей износа, по которой может быть определен у - процентный износ.
Прогнозирование индивидуального энергопотребления множества потребителей на срок продления эксплуатации СТ назначаемый на многие годы является сложной задачей. Учитывая необходимую дискретность определения I, соизмеримую с тепловыми постоянными обмоток - единицы минут.
Для целей определения остаточного ресурса параметры прогноза факторов эксплуатации должны включать параметры быстрых изменений (бросков)
нагрузки, параметры суточных, недельных, сезонных, годовых циклов изменения нагрузки, температуры 0a, параметры среднегодовых токов потребления:
1d = Рг '1 гаг ' Nг + A/ds ' Pc ' (Nдг ) + A/dc ' Рд (Nдн) + (7)
+ A[d4 • Рч (Nч ) + A/d; • K, (N,) где Рг - коэффициент изменения среднегодового тока; Isav - среднегодовой ток нагрузки на момент прогноза; N - номер годового интервала прогноза; AIds - среднегодовое изменение тока нагрузки; Рс - коэффициент среднесуточных изменений; Nds - номер дня в году; A/dc - среднесуточные изменения тока нагрузки; Pd(NdH) - коэффициент зависимости от дня недели; NdH
- номер дня в неделе; AId4 - среднечасовое изменение тока нагрузки; Рч - коэффициент, зависящий от номера часа; N4 - номер часа в сутках; A/di - амплитуда броска тока нагрузки; Ki - интервал броска тока нагрузки; N, - номер интервала броска тока нагрузки.
Для формирования значений переменных Ns, Nds, NdH, N4, N в компьютерную модель расчета расхода ресурса вводится блок календаря. Значения коэффициентов Рс, Рч, Рд, ток Isav, изменения токов AIds, AId4 определяются по данным мониторинга до момента прогнозирования. Значение Рг - из прогнозов роста электропотребления в зоне нагрузок СТ. Для определения AIdi и K, воспользуемся методикой [9].
Результаты моделирования по предложенной методике существенно отличаются от результатов статистической оценки по методике [1, 3].
ВЫВОДЫ
Проведенный анализ показывает, что все фазы жизненного цикла СТ оказывают существенное влияние на ресурс изоляции. Основным инструментом определения остаточного ресурса является имитационное моделирование процессов старения. Определяющим фактором повышения точности оценки является повышение точности прогнозирования динамики факторов изоляции за счет одновременного учета факторов из разных временных диапазонов - диапазонов скачков нагрузки, суточных, недельных и сезонных циклов, долголетних прогнозов нагрузки. Дальнейшие исследования будут направлены на разработку и исследование компьютерной модели СТ более полно отражающей взаимосвязи дестабилизирующих факторов эксплуатации и параметров изоляции, а также разработку моделей формирования реального ресурса.
СПИСОК ЛИТЕРАТУРЫ
1. Васин В.П., Долин А.П. К задаче оценки остаточного ресурса изоляции силовых маслонаполненных трансформаторов // Новое в российской электроэнергетике. - 2008. -№3. - С. 42-55.
2. IEC 60076-7 Ed. 1: Power transformers - Part 1: Loading guide for oil-immersed power transformers. Vol. 14/512/FDIS, Sept. 2005.
3. Васин В.П., Долин А.П. Ресурс изоляции силовых маслонаполненных трансформаторов // ЭЛЕКТРО. Электротехника, электроэнергетика, электротехническая промышленность. - 2008. - №3. - С. 12-17.
4. Поляков М.А., Василевский В.В. Моделирование динамики распределений дефектов в изделиях в ходе технологического процесса // Системные технологии. - 2013. -№4(87). - С. 73-81.
5. Поляков МА., Климов С.И. Методы и информационные технологии обработки данных мониторных параметров силового трансформатора // Вісник національного університету "Львівська політехніка". - 2009. - №637. - С. 70-73.
6. Конограй С.П. Применение модели старения твердой изоляции силовых маслонаполненных трансформаторов для их диагностики в режиме эксплуатации // Електротехніка і електромеханіка. - 2010. - №1. - С. 43-45.
7. Поляков МА., Мирошниченко A.r. Исследование процессов идентификации тепловых параметров силового масляного трансформатора по данным мониторинга // Вісник східноукраїнського національного університету. - 2007. -№11. Ч.1(117). - С. 1б0-1б7.
В. Поляков МА. От мониторинга параметров - к мониторингу состояний силового трансформатора // Електротехніка і електромеханіка. - 2011. - №1. - C. 49-52.
9. Поляков МА. Конструирование и выбор вейвлетов для анализа тепловых процессов в силовом трансформаторе // Технічна електродинаміка. - 2012. - №3. - C. 119-120.
REFERENCES: 1. Vasin V.P., Dolin A.P. K zadache ocenki ostatoch-nogo resursa izoljacii silovyh maslonapolnennyh transformatorov. No-voe v rossijskoj elektroenergetike, 200В, no.3, pp. 42-55. 2. IEC б007б-7 Ed. 1: Power transformers - Part 1: Loading guide for oil-immersed power transformers. Vol. 14/512/FDIS, Sept. 2005. S. Vasin V.P., Dolin A.P. Resurs izoljacii silovyh maslonapolnennyh transformatorov. ELEKTRO. Elektrotehnika, elektroenergetika, elektrotehnicheskaja promyshlennost', 200В, no.3, pp. 12-17. 4. Poliakov M.A., Vasilevskij V.V. Modelirovanie dinamiki raspredelenij defektov v izdelijah v hode tehnologicheskogo processa. System technologies, 2013, no.4(87), pp. 73-В1. 5. Poliakov M.A., Klimov S.I. Metody i informacionnye tehnologii obrabotki dannyh monitornyh parametrov silovogo transfor-matora. Bulletin of Lviv Polytechnic National University, 2009, п0.637, pp. 70-73. б. Konogray S.P. A solid insulation ageing model
application to oil-filled power transformer for operation diagnostics.
Electrical engineering & electromechanics, 2010, no.1, pp. 43-45. 7. Poliakov M.A., Miroshnichenko A.G. Issledovanie processov identifi-kacii teplovyh parametrov silovogo masljanogo transformatora po dannym monitoringa. Bulletin of V. Dahl East Ukrainian National University, 2007, no.11, part 1(117), pp. 160-167. 8. Poliakov M.A. From monitoring of parameters - to monitoring of the states of power transformer. Electrical engineering & electromechanics, 2011, no.1, pp. 49-52. 9. Poliakov M.A. Konstruirovanie i vybor vejvletov dlja analiza teplovyh processov v silovom transformatore. Technical electrodynamics, 2012, no.3, pp. 119-120.
Поступила (received) 25.10.2013
Поляков Михаил Алексеевич1, доц., к.т.н.,
Василевский Владимир Валентинович1, аспирант,
1 Запорожский национальный технический университет, 69063, Запорожье, ул. Жуковского, 64, тел/phone +38 066 794 60 37, e-mail: polyakov@zntu.edu.ua, Lisses@ukr.net
M.A. Poliakov1, V.V. Vasilevskij1 1 Zaporizhzhya National Technical University 64, Zhukovsky Str., Zaporizhzhya, 69063, Ukraine Evaluation of power transformer insulation residual life based on its individual life cycle characteristics.
Analysis reveals that all phases of a power transformer life cycle have a significant effect on the isolation resource. The paper intro-duces a technique for evaluating residual life of power transformer insulation via simulation of predicted results of external factor actions on the transformer insulation condition. Key words - insulation, power transformer, technical maintenance, failure, energy consumption.
