CC BY
164
УДК 681.518
СИСТЕМА ДИАГНОСТИРОВАНИЯ СИЛОВОГО МАСЛЯНОГО ТРАНСФОРМАТОРА ПО ТЕПЛОВЫМ ПАРАМЕТРАМ
Кандидаты техн. наук, доценты ЗАЛИЗНЫЙ Д. И., ШИРОКОВ О. Г.
Гомельский государственный технический университет имени П. О. Сухого
Из-за продолжающегося старения трансформаторного парка актуальность применения современных систем диагностики силовых трансформаторов в Республике Беларусь возрастает. К сожалению, отечественная промышленность практически не выпускает подобные системы, и предприятия вынуждены закупать их за рубежом, затрачивая значительные средства.
Системы диагностики силовых трансформаторов без отключения от сети, предлагаемые зарубежными производителями, такими как ABB [1], характеризуются значительным количеством цепей для измерения токов, напряжений, температур, влажности, химического состава масла, параметров устройств охлаждения, параметров РПН и др. В результате аппаратная часть этих систем имеет высокую стоимость, не говоря о стоимости программного обеспечения.
В данной статье рассматривается разработанный авторами аппаратно-программный комплекс, позволяющий при минимальных аппаратных затратах на ранней стадии выявлять неисправности, медленно развивающиеся в трансформаторе и приводящие к его дополнительному, анормальному нагреву [2].
Функциональная схема аппаратной части предлагаемой системы диагностирования [3] приведена на рис. 1.
Система состоит из прибора контроля и расчета температур силового трансформатора (ПКРТ), компьютера (К) и двух датчиков температуры, расположенных в одном корпусе.
Прибор ПКРТ имеет два канала измерения температуры: поверхность бака и воздух. Датчики температуры расположены в специальной кон-
струкции, монтируемой на поверхности бака трансформатора 1. Эта конструкция состоит из пластмассового корпуса 2, в котором сделаны отверстия для размещения датчиков. Датчик температуры поверхности бака 3 расположен в полом металлическом цилиндре 4, заполненном теплоизоляционным материалом 5. Металлический цилиндр закреплен в пластмассовом корпусе 2. Датчик температуры воздуха 6 расположен в аналогичном металлическом цилиндре 7. Для крепления конструкции к поверхности бака предусмотрен постоянный магнит 8. Датчики температуры подключаются к прибору ПКРТ посредством экранированных кабелей 9. Основная задача всей конструкции - обеспечить хороший тепловой контакт датчика 3 с поверхностью бака трансформатора и поместить датчик 6 в те же условия, в которых находится поверхность бака по отношению к окружающей среде.
Внешний вид всех предлагаемых устройств приведен на рис. 2.
Рис. 2. Внешний вид аппаратной части системы диагностирования силового масляного трансформатора по тепловым параметрам
Действующие макеты прибора ПКРТ и датчиков температуры собраны на кафедре «Электроснабжение» ГГТУ имени П. О. Сухого, а в качестве компьютера можно, например, использовать промышленный компьютер фирмы Advantech.
Прибор ПКРТ реализован на основе восьмиразрядного микроконтроллера фирмы Atmel AT90S8535 и оснащен жидкокристаллическим дисплеем, двумя кнопками, сигнальными светодиодами и разъемами для подключения датчиков температуры и компьютера. Связь с компьютером осуществляется по интерфейсу RS232.
В целом работа системы реализуется следующим образом. Прибор ПКРТ с помощью соответствующих датчиков температуры каждые 10 минут (значение периода дискретизации обосновано в [4]) измеряет температуры поверхности бака трансформатора и воздуха, передавая их значения в компьютер и отображая на своем дисплее. Программное обеспечение компьютера с помощью алгоритмов, о которых будет сказано ниже, рассчиты-
вает значения коэффициента нагрузки трансформатора, температуры его масла, обмоток и магнитопровода по отношению к высоте, на которой расположен датчик температуры поверхности его бака, а также диагностирует трансформатор на анормальный нагрев. Результаты расчетов компьютер передает обратно в прибор ПКРТ, который отображает необходимые значения и сообщения на дисплее. На лицевой панели прибора ПКРТ имеются сигнальные светодиоды «Перегрев масла», «Перегрев обмоток», «Связь с компьютером» и др.
Очевидно, что основной функциональной частью предлагаемой системы диагностирования являются расчетные алгоритмы и соответственно математический аппарат, на котором они базируются. Рассмотрим этот вопрос подробнее.
В [4] предложены математическая модель и алгоритмы, полностью соответствующие требованиям для предлагаемой системы диагностирования, но не прошедшие экспериментальной апробации. Для проверки их адекватности и возможной корректировки авторами данной статьи с помощью разработанных датчиков и прибора ПКРТ проведены экспериментальные исследования на автотрансформаторе АТДЦТН 63000/220/110/6. В процессе исследований измеряемые температуры записывали во встроенную Аа8Ь-память прибора ПКРТ. Также дежурными по подстанции несколько раз в сутки производилась запись значений температуры верхних слоев масла автотрансформатора со штатного термометра. После этого работу алгоритмов анализировали на компьютере.
Было проведено несколько экспериментов. В первом случае датчики устанавливали в верхней части бака на солнечной стороне вблизи от охладителей автотрансформатора.
Исследования показали плохую работоспособность алгоритмов, приведенных в [4]. Расчетный коэффициент нагрузки трансформатора претерпевал значительные колебания как в положительную, так и в отрицательную сторону, делая все дальнейшие вычисления некорректными.
Было предложено фильтровать расчетный коэффициент нагрузки через программное апериодическое звено первого порядка с постоянной времени, в семь раз превышающей период дискретизации. Алгоритм фильтрации записывается следующим образом:
К2 +—К2
Н.] гр н.фл.,7-1
^н.фл, =--> (!)
Т,
где Кн фл^ - расчетный коэффициент нагрузки после фильтрации; у - номер интервала расчета; Кн ^ - расчетный коэффициент нагрузки до фильтрации; - период дискретизации.
После применения алгоритма (1) величины К„ф ■ не превысили 0,5, но
наблюдались их отрицательные значения, что некорректно для расчета температур. Поэтому первый эксперимент был признан неудачным. Очевидно, здесь значительное влияние на датчик температуры оказало прямое солнечное излучение, а также сказалась близость охладителей.
Во втором эксперименте датчики температуры были установлены в верхней части бака на теневой стороне автотрансформатора на максимальном расстоянии от охладителей.
Значения температур воздуха, поверхности бака, масла исследуемого автотрансформатора во втором эксперименте и результаты расчетов температуры масла по алгоритмам, приведенным в [4] с учетом (1), показаны
Рис. 3. Экспериментальные и расчетная диаграммы для автотрансформатора АТДЦТН 63000/220/110/6: 1, 2, 3 - измеренные температуры воздуха, поверхности бака и масла;
4 - расчетная температура масла
В процессе расчетов отрицательных значений . не наблюдалось.
Как видно из рис. 3, отклонения расчетных температур масла от экспериментальных велики. Максимальная абсолютная погрешность расчета температуры масла по отношению к эксперименту составила 13,2 °С, что недопустимо. Коэффициент корреляции равен 0,83, что свидетельствует об определенной близости расчетных и экспериментальных значений в динамике.
По мнению авторов данной статьи, основной причиной высокой погрешности расчета является реальная зависимость коэффициентов теплоотдачи с масла в воздух от температуры. То есть тепловое сопротивление Л4 в тепловой схеме замещения трансформатора [4] должно быть нелинейным, как показано на рис. 4.
На рис. 4 обозначены [4]: С - С - соответственно теплоемкости маг-нитопровода, обмотки, масла и внутренней окружающей среды; ^ - Л4 -тепловые сопротивления магнитопровода, обмотки, масла и внутренней
окружающей среды; 9, - 9- - температуры магнитопровода, обмотки, масла, внутренней окружающей среды (слоя воздуха, температура которого зависит от температуры поверхности бака) и внешней окружающей среды (воздуха); АРХХ - активная мощность потерь холостого хода трансформатора; А1'и - активная мощность нагрузочных потерь трансформатора.
Тепловые сопротивления (кроме Д) можно рассчитать в соответствии с [5] и ГОСТ 14209-85:
^2 =
А
АР
х>
д
шг
АР.
(2)
й.
А1> + АР,
где г)с ,л - номинальное превышение температуры магнитопровода над температурой масла (г)с м = 25 °С); г)[|[п ,л - номинальное превышение температуры наиболее нагретой точки (ННТ) трансформатора над температурой масла (т31[1п ,л = 23 °С); О,, - номинальное превышение температуры масла над температурой окружающей среды (Фм = 55 °С для трансформаторов с системами охлаждения М и Д; Фм =40 °С для трансформаторов с системами охлаждения ДЦ и Ц).
Теплоемкости можно найти по формуле
С,=судД, /е[1;4], (3)
где суд ( - удельная теплоемкость 7-го однородного тела, Вт • с/(°С • кг); ^ - масса /-го однородного тела, кг.
Тепловое сопротивление Д определим следующим образом:
1
ам.к+«б.к+аб.и
(4)
где амк, абк, аби - коэффициенты теплоотдачи соответственно с масла конвекцией, с поверхности бака конвекцией, с поверхности бака излучением, Вт/(°С • м2); площадь поверхности бака трансформатора, м2.
Пренебрегая перепадом температуры по толщине стенки бака трансформатора, величины амк, абк, аби определим по эмпирическим формулам в соответствии с [6, с. 383-385]:
Гь Л0-5
^ ; (5)
=38- е4-е5
абк=2,5- е4-е5
,50 у
293
л 0,25
273 + 9,
(6)
0.25
О 1 9
аби =2,85- е4-е5 • -0,95 1 +
/
(7)
5
100
V
Таким образом, тепловое сопротивление Я4 - это функция температуры поверхности бака 04 и температуры воздуха 65. Наличие нелинейного сопротивления Я4 несколько усложняет вывод расчетных соотношений. Поэтому будем использовать пошаговый расчет, считая Я4 линейным, не зависящим от температуры на временном интервале каждого шага. Тогда, учитывая, что &]' .. ~ const, для первого шага расчета на временном интервале от 0 до периода дискретизации Td и нулевых начальных условиях операторное изображение температуры поверхности бака 04 j р запишется в следующем виде [4]:
где а, Ь и ё - коэффициенты, зависящие от тепловых сопротивлений и теплоемкостей [4]; Кн1 - значение коэффициента нагрузки трансформатора для первого шага расчета, которое считаем неизменным; 95л - значение температуры воздуха для первого шага расчета, которое также считаем неизменным.
Очевидно, что если учесть начальные условия, то соотношение (8) будет справедливо и для каждого /-го шага расчета при соответствующих значениях К 0- / и перерасчете теплового сопротивления Я4 по формулам (4)-(7).
Найдя оригинал для выражения (8), далее можно воспользоваться всеми алгоритмами, приведенными в [4]. При этом формулы (2), (3) достаточно будет применить только вначале расчета. Остальные расчеты необходимо вести циклически, начиная с (4)-(7) и далее по всем этапам алгоритмов, рассмотренных в [4], с учетом (1).
Экспериментальные диаграммы температур воздуха, поверхности бака и масла, а также расчетные диаграммы температур масла, обмотки и маг-нитопровода с учетом нелинейности теплового сопротивления Л4 приведены на рис. 5. Диаграмма расчетного эквивалентного квадрата коэффициента нагрузки показана на рис. 6. Эти значения были использованы для расчетных диаграмм на рис. 5.
Максимальная абсолютная погрешность расчета температуры масла по отношению к экспериментальным данным составила 5,3 °С, а коэффициент корреляции между расчетными и измеренными значениями этой температуры равен 0,94, что свидетельствует о хорошей адекватности математической модели. Расчетные значения температур масла, обмоток и магнито-провода можно использовать для защиты трансформатора от перегрева [3]. В случае задач диагностирования необходимы более подробный анализ и дополнительные алгоритмы. Рассмотрим возможные варианты.
е41 р -ДЗД +
к2
bxp + b0 —d3p3 + d2p2 + dxp + 1
_P_
о
e
(8)
5.1
4 3 2 i
a4p +a3p +a2p +a^p +1
, °с
56
6
4 5
.Г4 ^оос 3 >С^Ос > о (уЧ О л
к Л- \ А Л.
X д^4
288 г, ч 336
Рис. 5. Экспериментальные и расчетные диаграммы температур автотрансформатора АТДЦТН 63000/220/110: 1, 2, 3 - экспериментальные диаграммы температуры воздуха, поверхности бака и масла; 4, 5, 6 - расчетные диаграммы температуры масла, обмотки и магнитопровода
2,0
288 г, ч 336
Рис. 6. Диаграмма расчетного квадрата коэффициента нагрузки
При медленном развитии неисправностей в трансформаторе, например разбухании изоляции с перекрытием масляных каналов или увеличении диэлектрических потерь трансформаторного масла, возникает его общий анормальный нагрев [2]. Очевидно, что температура поверхности бака, подвергаемая измерению в предлагаемой системе диагностирования, также возрастет, что, в свою очередь, приведет к возрастанию расчетного коэффициента нагрузки и всех остальных расчетных температур.
В качестве диагностических параметров выберем средние значения расчетных температур за некоторые равные интервалы времени. Эти интервалы должны превышать величину максимальной из тепловых постоянных времени трансформатора, чтобы анормальный нагрев успел проявиться на фоне естественных нестационарных тепловых процессов. Примем длительность расчетных интервалов равной 24 ч. Тогда, например для температуры обмоток, диагностический параметр необходимо рассчитывать по формуле
2.диагл
17—1 17—1
_ ./=0
N.
2.Л./
(9)
где 62д1ШГ1 - значение диагностической температуры обмоток; / - номер расчетного интервала, соответствующий периоду дискретизации [4]; Т -период дискретизации (10 мин); / - номер расчетного интервала диаг-
42
28
14
0
48
96
44
192
240
0
48
96
144
192
240
ностирования; Тдиаг - расчетный интервал времени для диагностирования (24 ч); ^диаг - количество отсчетов температуры за расчетный интервал времени для диагностирования.
Диаграмма расчетных значений 02дпаг,, полученных на основе диаграммы температуры обмоток, приведенной на рис. 5, за последние 11 суток для исключения ошибки начальных условий, которая наблюдалась в течение первых трех суток расчета, показана на рис. 7.
60-—-
02.диаг.г ^^^^^^
50
40
30-
0 1 2 3 4 5 6 7 8 9 10 / 11
Рис. 7. Диаграмма расчетных значений диагностической температуры обмоток автотрансформатора АТДЦТН 63000/220/110
Очевидно, что за столь короткий интервал времени, который показан на рис. 7, нельзя судить о наличии или отсутствии анормального нагрева трансформатора. Колебания диагностической температуры обусловлены естественными причинами: изменением температуры воздуха, токов нагрузки трансформатора и других параметров. Только монотонное длительное возрастание диагностического параметра будет свидетельствовать об анормальном нагреве трансформатора.
В Ы В О Д
Предложенная система диагностирования силового трансформатора находится на этапе опытных образцов и требует длительных испытаний на силовых трансформаторах. Основное ее отличие от существующих систем - возможность раннего выявления неисправностей трансформатора при минимальных аппаратных затратах. Для установки системы не требуется отключать трансформатор от сети и вмешиваться в штатные схемы измерений. Ее внедрение позволит повысить надежность эксплуатации силовых масляных трансформаторов, а также будет способствовать переходу к их ремонтам по текущему состоянию.
Л И Т Е Р А Т У Р А
1. S t e n e s t a m, B. Electronic transformer control with monitoring & diagnostics / B. Ste-nestam, C. Bengtsson, P. Ydersten // CIGRE 2006, rep. A2-110.
2. Ш и р о к о в, О. Г. Основные теоретические положения по тепловой защите масляного трансформатора / О. Г. Широков, В. В. Курганов, Д. И. Зализный // Энергетика... (Изв. высш. учеб. заведений и энерг. объединений СНГ). - 2008. - № 2. - С. 16-23.
3. У с т р о й с т в о защиты силового трансформатора от перегрева: пат. Респ. Беларусь № 8204, МПК (2012) H02H6/00 / Д. И. Зализный, О. Г. Широков, Г. О. Широков. - Заявл. 31.10.2011; опубл. 30.04.2012 (www.belgospatent.org.by).
4. З а л и з н ы й, Д. И. Расчет температур основных элементов силового масляного трансформатора на основе анализа температуры поверхности его бака / Д. И. Зализный, О. Г. Широков // Энергетика. (Изв. высш. учебн. заведений и энерг. объединений СНГ). -2012. - № 4. - С. 18-28.
5. Б о д н а р, В. В. Нагрузочная способность силовых масляных трансформаторов / В. В. Боднар. - М.: Энергоатомиздат, 1983.
6. В а с ю т и н с к и й, С. Б. Вопросы теории и расчета трансформаторов / С. Б. Васю-тинский. - Л: Энергия, 1970. - 432 с.
Представлена кафедрой электроснабжения Поступила 13.03.2013
