УДК: [628.2+620.1](035)
МЕТОДИКА ОБРАБОТКИ РЕЗУЛЬТАТОВ ТЕПЛОВИЗИОННЫХ ОБСЛЕДОВАНИЙ ВЫСОКОВОЛЬТНОГО ОБОРУДОВАНИЯ
О. АЛЬ-АОМАРИ, Ю.В. ВАНЬКОВ, Е.Е. КОСТЫЛЕВА, Р.Н. ВАЛИЕВ Казанский государственный энергетический университет
В работе исследуются термограммы высоковольтного энергетического оборудования. Для ухода от субъективности при анализе термограмм, во время проведения тепловизионных исследований силового оборудования объектов энергетики, предложено использовать термографические информационные функции. Метод анализа реализован средствами математической среды ЫаЛСай. В качестве критерия оценки технического состояния силового оборудования используется коэффициент дефектности.
Ключевые слова: теплограмма, матрица температуры, методика анализа, неразрушающий контроль, коэффициент дефектности.
Современные тенденции в развитии промышленности и производства диктуют необходимость высоких требований к надежности оборудования. В России более 40 % силового оборудования, эксплуатируемого на объектах энергетики, выработали свой ресурс. Только слаженная и взаимоувязанная работа всех элементов единого энергетического комплекса обеспечивает безопасность работы энергетических объектов [1]. В связи с этим все более актуальными становятся вопросы, связанные с надежностью работы силовых энергетических объектов (Федеральный закон от 21.07.1997 N 116-ФЗ (ред. от 13.07.2015) "О промышленной безопасности опасных производственных объектов") [2,3]. Необходимый уровень эксплуатационной надежности силового оборудования можно поддерживать, применяя методы и средства неразрушающего контроля [4].
На сегодняшний день одним из перспективных направлений развития высокоэффективной системы технической диагностики в энергетике является внедрение приборов инфракрасной техники. Своевременное проведение тепловизионного обследования силового энергетического оборудования позволяет на ранних стадиях выявить возникающие в нем дефекты, прогнозировать их развитие [5,6]. Однако существует проблема, когда на анализ результатов тепловизионного обследования оборудования влияют субъективные факторы [7].
Все элементы работающего энергетического оборудования имеют определенную температуру. Элементы оборудования, имеющие дефекты, отличаются по температуре от нормально функционирующих элементов. Тепловизионное исследование позволяет видеть тепловое излучение и получать тепловое изображение (теплограмму) объекта, где определенному значению температуры соответствует определенный цвет.
Для ухода от субъективности при анализе теплограмм, полученных во время проведения тепловизионного исследования различных силовых объектов энергетики, предложено разработать методику компьютерного анализа теплограмм. Предложено использовать методику анализа термографических информационных функций (ТИФ).
Анализ теплограмм ведется в математической среде ЫаЛСАБ. Средствами ЫаЛСАБ рисунок превращается в матрицу [8], состоящую из чисел от 0 до 255. Числу 0 соответствует пиксель с минимальной, а числу 255 - с максимальной яркостью.
© О. Аль-Аомари, Ю.В. Ваньков, Е.Е. Костылева, Р.Н. Валиев Проблемы энергетики, 2015, № 11-12
Поскольку на теплограмме видно распределение температурного поля исследуемого объекта, то было установлено соответствие между значениями полученной матрицы и температуры в соответствующей точке исследуемого объекта. То есть, выполнено преобразование матрицы яркости теплограммы со значениями от 0 до 255 в матрицу распределения температур по поверхности объекта со значениями от ^т до tmax.
Анализу были подвергнуты три термограммы одного объекта (рис.1), полученные в разное время проведения регламентных работ.
1 зо б Sru« 47
'••«•»Pcml SFW 23 2012-0*41 10:00:44 40 ■ «120 —С 9«
TiTO .Sör.ru-МСЧ-.! ЛГОМ 23 гОЧМв-13 %Г 4?:5» -М- мго »-О» "С
(А) (Б) (С)
Рис. 1. Термограммы объекта: (А) - первая; (Б) - вторая; (С) - третья
Документ обработки первой термограммы (А) в программном продукте MathCAD представлен ниже.
ORIGIN = 1 - устанавливаем отсчет индексов от 1 (по умолчанию от 0) A = READGB (Рис. 1(A)) - чтение цветного изображения из файла. Получаем матрицу А: m = rows(A) = 356 число строк матрицы А n= cols(A) = 1728 число столбцов матрицы А Amin= min(A) = 0 Amax= max(A) = 255 N1 = Amax - Amin + 1=256
Преобразование матрицы яркости (со значениями от 0 до 255) в матрицу температур (со значениями от tmin до tmax):
™ , , t1max— t1min „„„„„., , , timax = 42, timin = 30, h1 =-—-= 0.047 K1 = 1..N1+1
N1
t1ki = ti min+ hi.(ki — i)
i = 1..N1 +1 , x = i- 1 At (m,n,M) = for i e 1..m for j e 1..n for k e 1..N1+1 s ^ k-1
A11,J ^tlk if My=s. Статистическая обработка матриц температуры дала возможность получить функцию распределения температуры (ТИФ) [9].
Определяем оптимальное число интервалов изменения температуры: N = round (1+3,323.1n(m.n)) = 45.
nA red =
hist(x, A_red)
m.n
- относительная частота значений х в матрицы
(рис. 2).
Рис. 2. ТИФ красного цвета
. hist(x,A green)
nA _ green =- - относительная частота значении х в матрицы
m.n
A_green (рис. 3).
Рис. 3. ТИФ зеленого цвета
. ,, hist(x, A blue) nA _ blue =--относительная частота значении х в матрице
m.n
A_blue (рис. 4).
Рис. 4. ТИФ синего цвета
. , ы81(х,а) „ . , ..
пА_оЬ =--относительная частота значений х в матрице А (рис. 5).
т.п
Рис. 5. Общая ТИФ обследуемого объекта (термограмма А)
Следующие две термограммы обработаны по этой же схеме. Результаты обработки представленны на рис. 6 и 7.
Рис. 6. Общая ТИФ обследуемого объекта (термограмма В)
Рис. 7. Общая ТИФ обследуемого объекта (термограмма С)
На графиках показано распределение относительной частоты регистрации различных температур на термограмме исследуемого объекта.
При недостатке статистических данных за «эталонную» ТИФ принимается ТИФ объекта с минимальным тепловыделением из нескольких (не менее двух) объектов в тождественных условиях.
В качестве иллюстрации на рис. 8 показано смещение графиков ТИФ трех теплограмм в сторону более высокой температуры. Оценка степени развитости локального дефекта производится в соответствии с табл.1
град. С
Рис. 8. Смещение графиков ТИФ трех теплограмм Оценка технического состояния объекта
Таблица 1
Оценка технического состояния Норма Норма с отклонениями Норма со значительными отклонениями Ухудшенное Предаварийное
Кдис До 1,2 1,2-1,4 1,4-1,6 1,6-2 Более 2
В качестве критерия для оценки мощности диссипационных явлений в энергетическом оборудовании используем коэффициент дефектности [10], определяемый соотношением
Кдис = ?2 р Р , С1)
м
где Р; - уровни мощности диссипации.
Вычислялась мощность диссипации для каждого цвета всех трех теплограмм по соотношению:
Р1 = {^ВД^Л; Р2 = {^2(1)^1. В среде ЫаЛСАБ этот расчет выглядит так:
вычисление коэффициентов дефектности - Красный цвет:
P1 = Jfl- FA _ red (t).tdt = 1,126
P2 = Jt2"" FA _ red (t).tdt = 2,442
P3 = J^"" FA _ red (t ).tdt = 2,836
Р 2 - Р1 Р3 - Р1 Р 3 - Р 2
К12 red =-= 1,169,К13 red =-= 1,519 , К23 red =-= 0,162 .
p1 p1 ' - Р2
Вычисление коэффициентов дефектности - Зеленый цвет:
P1 = J *W fa _ green(t ).tdt = 1,453 ti ■ 1min
P2 = J'2max FA _ green(t).tdt = 2,092
Jtj . 2min
P3 = Jt3max fa _ green(t).tdt = 2,299
t3 ■
3min
Р 2 - Р1 Р 3 - Р1
К12 green =-= 0,44, К13 red =-= 0,583,
Р1 ~Р1
Р3 - Р2
К23 red =-= 0,099 .
Р2
Вычисление коэффициентов дефектности - Синий цвет:
P1 = Jt1max fa _ blue(t).tdt = 1,601 t1
1min
P2 = J*2max fa _ blue(t).tdt = 2,691
t2min
P3 = Jt3max fa _ blue(t).tdt = 2,697
J to
3
Р 2 - Р1 Р 3 - Р1
К12 blue =-= 0,681, К13 red =-= 0,685,
Р1 ~Р1
Р 3 - Р 2 К23 red =-= 2.327 .
Р2
Вычисление общего коэффициента дефектности:
P1 = Jf1max FA _ ob(t).tdt = 1,393
P2 = Jt2max FA _ ob(t).tdt = 2,408
Jfo .
2min
P3 = Jt3max fa _ ob(t).tdt = 2,611.
3min
Р 2 - Р1 Р3 - Р1
К12 blue =-= 0,729, К13 red =-= 0,874,
Р1 Р1
Р3 - Р 2 К23 red =-= 0,084 .
Р2
Оценка технического состояния по тепловым условиям производится по уровню коэффициента Кдис в соответствии с табл. 1.
Коэффициент дефектности K13_red=1,532 соответствует норме со значительным отклонениями, коэффициент дефектности K23_red= 2,327 соответствует предаварийному состоянию. Все остальные коэффициенты в норме.
Заключение
Распределение температур по поверхности объектов несет информацию следующего характера: о наличии распределенных источников тепловыделения, эффективности системы охлаждения, наличии локальных температурных аномалий, которые обусловлены скрытым дефектом термического характера.
Для ухода от субъективности при анализе термограмм, полученных во время проведения тепловизионного исследования различных силовых объектов энергетики, предложено разработать методику компьютерного анализа термограмм. Предложено использовать методику анализа термографических информационных функций. Метод анализа реализован средствами математической среды MathCad. Проведена апробация метода при анализе термограмм реальных объектов энергетики.
В качестве критерия оценки технического состояния трансформатора используется коэффициент дефектности.
Summary
The conducted researches allowed to reveal shortcomings of the existing techniques of thermo vision control of an assessment of a condition of transformers one of which is low extent of automation of the analysis of thermo grams.
Keywords: Power transformer, thermo grams, Mathcad programs, thermo graphic functions, analysis of pictures.
Литература
1. Аль-Аомари О.М. Результаты тепловизионного обследования силовых трансформаторов ОАО «ТАТНЕФТЬ» // Научному прогремссу - творчество молодых: материалы Х международной молодежной научной конференции по естественным и
техническим дисциплинам. (Йошкар-Ола, 17-18 апреля 2015г.): в 2 ч. / редкол.: Д.В. Иванов [и др.]. Йошкар-Ола: Поволжский государственный техноолгический университет, 2015. Ч.1. С.225.
2. Ванин Б.В., Львов Ю.Н., Львов М.Ю. и др. О повреждении силовых трансформаторов напряжением 110-500 кВ в эксплуатации // Электрические станции. Ежемесячный производственно-технический журнал. 2001. № 9. С.53-58.
3. Лила В.Б., Костюков А.В. Экспертная система диагностики силовых трансформаторов // «Инженерный вестник Дона». 2013. №1. URL: ivdon.ru/ru/magazine/archive/n1y2013/1504.
4. Долин А.П., Першина Н.Ф., Смекалов В.В. Опыт проведения комплексных обследований силовых трансформаторов // Электрические станции. Ежемесячный производственно-технический журнал. 2000. №6. С. 46-52.
5. РАО ЕЭС РФ «Основные положения методики инфракрасной диагностики электрооборудования и ВЛ». РД153-34.0-20.363-99. М.: Изд-во НЦ «ЭНАС», 2004. 85 с.
6. Объем и нормы испытаний электрооборудования. РД 34.45-51.300-97. М.: Изд-во НЦ «ЭНАС», 2000. 256 с.
7. Вавилов В.П. Инфракрасная термография и тепловой контроль. М.: Изд-во Спектр, 2009. 562 с.
8. Аль-Аомари О., Ваньков Ю.В., Костылева Е.Е., Валиев Р.Н. Методика обработки термограмм силовых объектов энергетики // Инженерный вестник Дона. 2015. №3. URL: http://www.ivdon.ru/ru/magazine/archive/n3y2015/3261.
9. Yuri,V,Vankov and Osamah, M, Al-aomari. Technique of Processing of Thermo grams of Power Transformers // International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 4, Issue 5, September 2015. рр.134-138.
10. Козлов В.К., Гарифуллин М.Ш., Валиуллина Д.М. Перспективы использования спектральных методов при диагностике состояния маслонаполненных трансформаторов // Энергетика Татарстана. 2005. №2. С. 36-40.
Поступила в редакцию 14 декабря 2015 г.
Осама Аль-Аомари - аспирант кафедры «Промышленная теплоэнергетика и системы теплоснабжения» (ПТЭ) Казанского государственного энергетического университета (КГЭУ). Тел: 8(927)248-22-61. E-mail: al_aomari2010@mail.ru.
Ваньков Юрий Витальевич - д-р техн. наук, профессор, заведующий кафедрой «Промышленная теплоэнергетика и системы теплоснабжения» (ПТЭ) Казанского государственного энергетического университета (КГЭУ). Тел: 8(903)305-62-51. E-mail: yvankov@mail.ru.
Костылева Елена Евгеньевна - канд. техн. наук, доцент кафедры «Промышленная теплоэнергетика и системы теплоснабжения» (ПТЭ) Казанского государственного энергетического университета (КГЭУ). Тел: 8(919)699-43-90. E-mail: dunysha@bk.ru.
Валиев Радик Нурттинович - канд. техн. наук, доцент кафедры «Промышленная теплоэнергетика и системы теплоснабжения» (ПТЭ) Казанского государственного энергетического университета (КГЭУ). Тел: 8(903)306-57-63. E-mail: valievkgeu@mail.ru.