Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

CALCULATION OF HEAT LOSSES IN TANK OF POWER TRANSFORMER AT GEOMAGNETIC STORMS

Kretov D.A. 1
1 Togliatti State University, Togliatti, Russia
In article the developed design procedure of additional thermal losses in a tank of the power transformer is presented at course of the geoinduced currents at geomagnetic storms. At course on windings of the power transformer of the geoinduced current there is a saturation of magnetic system of the power transformer and increase in a stream of dispersion that leads to additional losses of active capacity and heating of windings of the power transformer and massive metal constructive elements. For calculation of additional thermal losses in a tank the TDTs real power transformer model – 400000/220 is used. At considerable saturation of magnetic system calculation of a magnetic field near walls of a tank is executed without magnetic system of the transformer. Windings of each phase were replaced with one infinitely thin winding on which the current of magnetization brought to a winding of a high voltage flows. At calculation influence of all three phases was considered. For an assessment of thermal losses in a tank of the power transformer the currents received as a result of calculations, were displayed on harmonious components. Losses paid off for each harmonica separately, and then the received losses were summarized. As a result of calculations values of specific and total thermal losses in a tank of the power transformer are received at course on its windings of the geoinduced current. The developed technique can be put into practice for calculation of heating of walls of a tank of the power transformer, and as for an assessment of extremely excess of temperature according to a standard at geomagnetic storms.
thermal losses.
geomagnetic storms
power transformer
Протекание геоиндуцированных токов (ГИТ), возникающих при геомагнитных бурях (ГМБ), по обмоткам силового трансформатора (СТ) приводит к насыщению магнитной системы и увеличению потоков рассеяния. Увеличение потока рассеяния вызывает дополнительные потери активной мощности от вихревых токов в обмотках силового трансформатора и в массивных металлических конструкциях. Электрическая энергия, идущая на увеличение дополнительных потерь активной мощности, полностью преобразуется в тепловую энергию. Дополнительно выделившееся тепло в обмотках и массивных металлических конструкциях увеличивает их температуру, которая при определенных условиях может превышать предельно допустимую. В стандартах МЭК 354-91 и ГОСТ 14209-97 [2] указаны предельные значения температуры при перегрузках силовых трансформаторов. Предельная температура обмоток мощных силовых трансформаторов при продолжительных аварийных перегрузках - 130 °С, а при кратковременных - 160 °С, а масла в верхних слоях - 115 °С.

Главным риском во время кратковременной (менее 30 минут) аварийной перегрузки является снижение электрической прочности вследствие образования пузырьков газа, которые могут возникать в бумажной изоляции, когда температура наиболее нагретой точки обмотки силового трансформатора достигает 140 - 160 °С при нормальном содержании влаги в изоляции [2].

Расчёт потерь в баке при насыщении магнитной системы проведён по методике, изложенной в [3]. Для расчёта потерь в баке была использована модель реального СТ типа ТДЦ - 400000/220 (рис.1). Бак представлен в виде прямоугольного параллелепипеда с толщиной стенок, превышающей глубину проникновения магнитного поля. При значительном насыщении магнитной системы расчёт магнитного поля вблизи стенок бака проводился без учета магнитной системы трансформатора. Обмотки каждой фазы заменялись одной бесконечно тонкой обмоткой, по которой течёт ток намагничивания, приведённый к обмотке высокого напряжения. При расчёте учитывалось влияние всех трёх фаз A, B, C.

Рис. 1. Расчетная модель реального трехфазного силового трансформатора ТДЦ - 400000/220 (размеры даны в метрах, число витков обмотки ВН равно w = 364)

Методика, по которой проводились расчеты, состоит из четырёх этапов: 1) по конструкции трансформатора и токам обмоток рассчитывается напряженность магнитного поля в воздухе для трансформатора со снятым баком Н0; 2) по Н0 и коэффициенту kH находится касательная составляющая напряжённости магнитного поля на внутренней поверхности бака ; 3) по  и активному поверхностному сопротивлению материала бака rп находятся потери на единицу поверхности pп; 4) общие потери определяются суммированием удельных потерь pп по поверхности бака.

В изложенную методику внесено изменение: вместо интегрирования по поверхности бака рассчитывались среднеквадратичные значения напряжённости магнитного поля на трёх площадках, расположенных так же, как и по поверхности бака, - nx (внешняя нормаль , направленная вдоль оси x), ny (внешняя нормаль ), nz (внешняя нормаль ).

Учитывая, что по своему влиянию на магнитное поле стенка бака силового трансформатора является электромагнитным экраном, коэффициент kH был принят равным двум (kH = 2) [4].

Протекание в электроэнергетической системе (ЭЭС)  ГИТ приводит к насыщению стали силовых трансформаторов и появлению высших гармоник тока и напряжения. Ток намагничивания трансформатора при этом является несинусоидальным и имеет выраженную однополупериодную несимметрию.

Для оценки потерь в баке токи в обмотках, полученные в расчетах, раскладывались на гармонические составляющие. Потери рассчитывались для каждой гармоники в отдельности, а затем полученные потери суммировались. Следует отметить, что предложенный подход является приближённым, т.к. не учитывает намагничивание стали бака.

Удельные потери pп рассчитывались по формуле:

   (1)

где pп - удельные потери на единицу поверхности;  - действующее значение тангенциальной составляющей напряжённости магнитного поля у поверхности бака; rп - удельное поверхностное сопротивление; kφ - коэффициент, равный 1,1-1,2 для конструкционной стали (при расчётах полагали kφ=1,15).

Удельное поверхностное сопротивление rп определялось как проводимость скин-слоя:

 (2)

где  толщина скин-слоя, ρ - удельное электрическое сопротивление конструкционной стали, ω - частота переменного электромагнитного поля; μэк - эквивалентная магнитная проницаемость. В [5] показано, что при частоте f0=50 Гц для конструкционных сталей толщина скин-слоя δ лежит в диапазоне δ=0,7 - 1,5 мм и зависит от напряжённости магнитного поля у поверхности. При расчётах полагали ρ=0,144·10-6 Ом·м, а δ = 1,1 мм, что соответствует эквивалентной магнитной проницаемости конструкционной стали μэк=603.

Таким образом, на частоте 50 Гц поверхностное сопротивление rП,1 составляет rП,1=1,31·10-4 Ом. На частотах высших гармоник n поверхностное сопротивление рассчитывалось по формуле:

       (3)

где n - номер гармоники.

С учётом сказанного формулу (1) для удельных потерь можно переписать в виде:

 (4)

где Н0,т - амплитудное значение напряжённости магнитного поля на заданной площадке, совпадающей с поверхностью бака, которое получается после среднеквадратичного усреднения соответствующих компонент поля. Для каждой поверхности nx, ny, nz результирующие удельные потери равны сумме удельных потерь от каждой гармоники. Полные потери в баке получали путем умножения удельных потерь на соответствующие площади поверхностей nx, ny, nz и удваивая полученный результат (вклад противоположных поверхностей - nx, - ny, - nz).

Расчёт магнитного поля кругового ленточного тока в воздухе был протестирован методом, изложенным в [4] для случая бесконечно тонкого кольца радиуса r0.

 (5)

 (6)

где I1 и K1 - модифицированные функции Бесселя первого и второго рода, а постоянная c0 равна:

 (7)

Для воздуха полагаем μ=1, в случае многовитковой катушки ток I=w·I0 (w - число витков, I0 - ток одного витка).

Напряжённость магнитного поля в любой точке пространства находилась по формулам:

 (8)

При численных расчётах использовалась квадратная плоская сетка с шагом .

Для гармонического анализа удобно использовать кривую тока намагничивания, изображённую на рис.2. В силу значительной симметрии в этом случае достаточно ограничиться только косинус-членами разложения. Кривая тока намагничивания на рис.2 соответствует величине геомагнитного тока равного iГИТ=65А.

Рис.2. Кривая тока намагничивания

Разложение в ряд Фурье тока намагничивания по высшим гармоникам определялось по следующим выражениям:

 (9)

 (10)

Из полученных результатов разложения кривой тока намагничивания (рис.2) выявлено, что коэффициенты bn много меньше коэффициентов аn первой гармоники и ими можно пренебречь. Таким образом, разложение будет иметь вид:

  (11)

В трёх обмотках токи сдвинуты по фазе. Поэтому разложение в ряд Фурье для функции тока , имеющей временной сдвиг h (если принять за i0(t) ток фазы A, то для фаз B и C имеем ). Учитывая, что bn=0, на основе формулы (10), разложение для сдвинутого тока по формуле (9) будет иметь вид:

 (12)

При  величина nωh соответствует фазовому сдвигу  соответственно для фаз B и C; при  - фазовому сдвигу ; при  - фазовому сдвигу Δφ=0. На комплексной плоскости первый случай соответствует расположению векторов фаз A, B и C в виде прямой последовательности, второй случай - в виде обратной последовательности, а третий - в виде нулевой последовательности.

Расчёт напряжённости магнитного поля проводился с помощью комплексных амплитуд:    и т.д. - верхний индекс обозначает площадку, первый нижний индекс обозначает компоненту поля, второй нижний индекс - обмотку (1 - фаза A, 2 - фаза B, 3 - фаза C), третий нижний индекс - номер гармоники. Сдвиг по фазе даёт разную ориентацию векторов комплексных амплитуд на комплексной плоскости, а пространственное размещение обмоток даёт разную величину этих векторов. Для комплексных амплитуд его можно записать следующим выражением (для краткости опущен верхний индекс площадки):

  (13)

где hx1, hx2, hx3 и т.д. - удельные компоненты напряжённостей, создаваемых обмотками 1, 2, 3 (напряженности, отнесённые к величине тока; верхняя строка - для прямой последовательности ( ), средняя - для обратной ( ), нижняя - для нулевой ( ); фазовые коэффициенты определены формулами:

 (14)

Сложение полей обмоток даёт:

 (15)

Анализ (15) показывает, что выражения в скобках не зависят от коэффициентов разложения тока аn, но зависят от типа последовательности. Обозначив их через  и т.д., где последний индекс означает тип последовательности, получим:

 (16)

При этом вместо 27 действительных удельных компонент напряжённостей   и т.д. для трёх обмоток (на трёх площадках) получаем 27 комплексных величин  и т.д. для трёх типов последовательностей (на трёх площадках).

Формула (16) показывает, что усреднение комплексных амплитуд поля можно заменить усреднением удельных напряжённостей, т.к. между среднеквадратичными значениями их модулей будет связь, аналогичная (16). Среднеквадратичные значения модулей комплексных удельных напряжённостей находится усреднением по соответствующим площадкам:

  (17)

Поскольку удельные напряжённости прямой и обратной последовательностей отличаются только фазой, а не величиной, то их среднеквадратичные значения будут равны:  Поэтому для каждой площадки имеем 6 среднеквадратичных значений, а всего - 18. Умножение этих значений на соответствующие коэффициенты гармоник аn даёт среднеквадратичные значения компонент напряжённостей , создаваемых каждой гармоникой.

Для каждой площадки существуют две тангенциальные компоненты магнитного поля, которые получаются после усреднения по формулам:

  (18)

Результаты расчета по (4) удельных потерь и суммарные потери бака представлены в табл. 1. Из табл. 1 видно, что полученная оценка потерь в баке силового трансформатора ТДЦ - 400000/220 составляет 223,9 кВт при условии, что стенки бака являются электромагнитным экраном по своему влиянию на магнитное поле.

Следует отметить, что при сильном насыщении магнитной системы силового трансформатора магнитная проницаемость электротехнической стали близка к магнитной постоянной, а относительная магнитная проницаемость стремится к единице [3], т.е. можно полагать, что основная часть рабочего магнитного потока будет замыкаться по баку. В этом случае дополнительные потери активной мощности в баке СТ, устанавливаемых на электрических станциях и трансформаторных подстанций электроэнергетических систем, при геомагнитных бурях, могут значительно превысить расчетные значения.

Таблица 1. Удельные потери по гармоникам (Вт/м 2) и суммарные потери

n

nx

ny

nz

n

nx

ny

nz

1        

53,39

384,28

268,45

11     

0,28

2,02

1,41

2        

61,16

440,18

307,50

12     

0,03

0,19

0,08

3        

98,85

681,09

263,49

13     

0,02

0,16

0,11

4        

43,45

312,72

218,46

14     

0,09

0,66

0,46

5        

30,49

219,47

153,31

15     

0,24

1,63

0,63

6        

34,68

238,99

92,46

16     

0,13

0,94

0,66

7        

11,61

83,54

58,36

17     

0,10

0,72

0,50

8        

6,27

45,11

31,51

18     

0,11

0,77

0,30

9        

5,26

36,23

14,01

19     

0,03

0,23

0,16

10        

1,13

8,13

5,68

20     

0,01

0,08

0,06

Суммарные удельные потери, Вт/м 2

347

2457

1417

Площадь поверхностей, м 2

2·9,7

2·33,6

2·18,2

Суммарные потери, кВт

6,8

165,3

51,8

Суммарные потери в баке, кВт

223,9

Выводы

В результате проведенных расчетов были получены значения удельных потерь и суммарных тепловых потерь в баке силового трансформатора при протекании по его обмоткам геоиндуцированного тока. Рассмотрен реальный силовой трансформатор, устанавливаемый на крупных объектах электроэнергетических систем, таких как электростанции и крупные подстанции. Разработанная методика может быть применена на практике для расчета нагрева стенок бака и обмоток силового трансформатора, а также для оценки предельно допустимого приращения температуры согласно стандарту при геомагнитных бурях.

Рецензенты:

  • Кувшинов Алексей Алексеевич, д.т.н., профессор, доцент кафедры «Электроснабжение и электротехника», федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Тольяттинский государственный университет», г. Тольятти.
  • Абрамов Геннадий Николаевич, д.т.н., профессор, профессор кафедры «Промышленная электроника», федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Тольяттинский государственный университет», г. Тольятти.