Широкое распространение в диагностике электрооборудования получило тепловидение. Картина инфракрасного излучения весьма содержательна с точки зрения диагностической информации о техническом состоянии объекта. Возможность получения поверхностного распределения температурного поля позволяет выявлять многие дефекты, в частности, перегревы внешних контактных соединений, так называемые «тепловые мешки», засорение элементов систем охлаждения трансформаторов и многие другие. Однако продуктивный анализ этой информации требует от исследователя опыта и определенного искусства, что, очевидно, препятствует более широкому и более эффективному внедрению метода в практику. Одним из главных препятствий при этом видится отсутствие методик, которые позволили бы формализовать процесс анализа тепловизионных изображений.
Целью настоящей статьи является разработка подхода к такой формализации. Воспользуемся одним из наиболее универсальных и весьма эффективных исследовательских инструментов, каким зарекомендовал себя подход, основанный на электроаналогии, и, в частности, на электротепловой аналогии. Напомним некоторые электротепловые аналоги: температура T соответствует электрическому потенциалу U, тепловой поток Q — электрическому току I, процесс теплоотдачи (или теплопередачи) характеризуется коэффициентом теплоотдачи α (или коэффициентом теплопроводности λ) и моделируется электрической проводимостью G или сопротивлением R. Суть подхода такова. Полученное с помощью тепловизора распределение температуры на баке трансформатора может быть интерпретировано как распределение потенциала электрического поля, созданного источником, расположенным внутри трансформатора (рис. 1).
В качестве прообраза электрического источника (источника ЭДС или источника тока) будем рассматривать интегральный тепловой источник — активную часть трансформатора. Изменение внутреннего состояния трансформатора или реактора, очевидно, скажется на изменении его внутреннего тепловыделения и, как следствие, на изменении внешнего температурного распределения на баке.
Для дальнейшего важно следующее допущение. Будем рассматривать только стационарное тепловое состояние трансформатора или реактора, то есть такой тепловой режим, при котором выделение тепла равно его отводу. Реально практически всегда существует изменение теплового и, соответственно, температурного распределения, однако будем считать такое изменение достаточно медленным. Разобьем непрерывную поверхность бака трансформатора SB на n площадок. В пределах каждой площадки будем считать температуру постоянной и равной ее среднему значению T i (i =1, …, n). В целях упрощения модели будем считать тепловой источник заключенным внутри объема VА, ограниченного некоторой сглаженной поверхностью SА, след которой показан пунктиром на рис. 1. Интегральный тепловой поток Q разобьем на элементарные трубки qi: Q = Σqi (одна из трубок показана пунктиром на рис. 1). Для общности не будем вводить требования равенства сечений трубок на поверхности бака (а также на поверхности SА), поскольку, как будет показано далее, такое ограничение неудобно для практической реализации подхода. Предположим, что имеется некоторая условная замкнутая поверхность Sо , расположенная на определенном расстоянии от трансформатора, содержащая внутри себя трансформатор со всеми его элементами, температура которой равна температуре окружающей среды T о (см. рис. 1).
Элементом этой поверхности, например, является часть площадки, на которой установлен трансформатор. Если продолжить трубки теплового потока, рассеченные поверхностью бака, до пересечения с указанной условной поверхностью (обозначим эту часть i-ой трубки qнi) , то построение упрощенной тепловой модели трансформатора можно считать практически завершенной. Сопоставим каждую элементарную площадку с узлом электрической цепи, каждую элементарную трубку теплового потока — с ветвью цепи, источник теплового потока — с источником ЭДС. Свяжем узлы электрической цепи, соответствующие элементам поверхности бака, дополнительными ветвями, моделирующими тепловую связь указанных элементов между собой за счет теплопроводности стенки бака. Пример построения описанной цепной модели показан на рис. 1. Отдельно остановимся на способе моделирования системы охлаждения, а также крышки и дна бака. На данном этапе разработки представляется, что имеет смысл моделировать систему охлаждения источником тока или ЭДС.
Поскольку в общем случае интенсивность теплового потока в трубопроводах системы охлаждения не известна, а температуры входа и выхода масла, наоборот, как правило, определены, то удобно эквивалентировать этот элемент трансформатора идеальным источником ЭДС. При этом целесообразно помещать точки входа и выхода в узлах цепной модели. Что касается моделирования крышки и дна бака, то, учитывая практическую невозможность получения термограмм этих участков поверхности трансформатора, имеет смысл предварительно представить их изотермами с температурой, равной средней температуре граничных точек Tк = ΣTj/m (j = jк1, …, jкm), Tд = ΣTj/m (j = jд1, …, jдm), отразив затем соответствующие допущения в электрической модели. Здесь m — число площадок разбиения поверхностей крышки и дна бака (с учетом их геометрического подобия). Наличием различных дополнительных элементов конструкции (газоразводных труб, расширителя, вводов и др.) при этом пренебрегаем. Вполне возможно, что более удобным может оказаться способ моделирования системы охлаждения пассивными «теплоотдающими» ветвями. Критерием при этом, как обычно, будет практическое удобство модели.
Рассмотрим теперь особенности моделирования процессов теплоотдачи и теплопередачи в баке трансформатора. Будем считать коэффициент теплоотдачи постоянным на всей поверхности бака αB = const. Положим также одинаковым коэффициент теплопередачи вдоль стенки бака λB = const. Что касается теплового сопротивления стенки бака в направлении нормали к ее поверхности, то вследствие высокой теплопроводности металла стенки бака и очевидного соотношения PB, hB>> d (PB — периметр бака, hB — высота бака, d — толщина стенки) примем его равным нулю. После принятия этих допущений получаем RTo (αB) Ro = 1/Go; RTc(λB) ↔ Rc =1/Gc, где RTo — тепловое сопротивление теплоотдачи в окружающую среду элементарной трубки теплового потока; RTc — тепловое сопротивление теплопередачи вдоль стенки бака между центрами двух элементарных площадок; Ro, Gо — соответственно сопротивление и проводимость, моделирующие теплоотдачу с элемента поверхности бака; Rc, Gc — сопротивление и проводимость, моделирующие передачу тепла посредством теплопроводности вдоль стенки бака; ↔ — символ соответствия. Рассмотрим теперь модель теплопередачи внутри трансформатора.
Очевидно, что и тепловые сопротивления разных элементарных трубок и интенсивность теплового потока через элементы поверхности SАi существенно различаются, что сильно усложняет постановку задачи. Примем для определенности, что тепловые сопротивления трубок одинаковы, тогда все различие в параметрах теплопередачи будет отнесено на счет неравномерности распределения элементарных тепловых источников по поверхности SА . Обозначим тепловое сопротивление одной трубки RTм, а его электрические аналоги — Rм и Gм. С учетом всех принятых допущений и обозначений получаем электрическую модель теплового состояния трансформатора, изображенную на рис. 2. 
Запишем систему уравнений полученной цепной модели
Примем в качестве базового (нулевого) узел с потенциалом Uо(соответствующий температуре Tо). Порядок нумерации остальных узлов может быть выбран произвольно. Получаем:
Определение значений величин Gij (i, j =1, …, n+2) представляет собой самостоятельную задачу. Для ее решения необходимо иметь значения коэффициентов теплоотдачи и теплопроводности (αB, αМ, λB). Если значение λB можно считать известным (например, из справочника), то значения αB,αB должны быть найдены расчетным или экспериментальным способом. Будем считать, что для решения общей задачи не требуется высокая точность и можно обойтись приближенными оценками указанных величин (которые могут быть взяты из данных лабораторных или натурных экспериментов, выполненных ранее). Таким образом, значения величин Gij будем считать априори известными. Рассмотрим на примере шунтирующего реактора РОДЦ-110000/750 использование данного подхода.
Выбор типа оборудования в данном случае был обусловлен максимальной открытостью поверхности бака (отсутствием навесных конструкций), наличием подробной совокупности термограмм (рис. 3), а также возможностью осуществить сопоставление с состоянием другого аналогичного реактора. Разобьем поверхность бака на участки так, как показано на рис. 4.
Следуя указанной нумерации, составим вектор потенциалов (температур) и матрицу проводимостей (в роли нулевого узла выступает условная внешняя поверхность So; номера 41 и 42 соответствуют крышке и днищу). Последующее умножение матрицы проводимостей на столбец потенциалов дает нам столбец токов, то есть тепловых потоков. Разделив значения тепловых потоков на соответствующие площади, получаем столбец с удельными тепловыми потоками. Таким образом от распределения температур по поверхностям бака мы переходим к распределению удельных тепловых потоков, исходящих с внутренней условной поверхности, содержащей в себе активную часть. В табл. 1 приведены результаты расчетов и графики распределения температур по поверхности бака и плотностей теплового потока по поверхности So. Структура матрицы проводимостей показана на рис. 5.
В табл. 2 показан вид одного из блоков матрицы.
ВЫВОДЫ
1. Предложенный метод работоспособен, но требует совершенствования.
2. Профили распределения тепловых источников для рассмотренных реакторов близки (вне зависимости от внешних условий и режима работы системы охлаждения).
3. Режим работы системы охлаждения, то есть интенсивность массопереноса трансформаторного масла, может быть учтен изменением коэффициента теплопроводности масла.
Автор:Алпатов М.Е., докт. техн. наук, Голованов А.А. ОАО «Электрозавод»
Добавить комментарий