Модель обмотки силового трансформатора для учета влияния квазипостоянного тока на режим работы силового трансформатора

Аннотация

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

Ключевые слова: силовой трансформатор, насыщение, моделирование, обмотка силового трансформатора, система электроснабжения, квазипостоянный ток, система электроснабжения, линии электропередачи, геомагнитные бури, космическая погода.

Изменения космической погоды вызывают относительно быстрые изменения геомагнитного поля Земли, что приводит к появлению в системах электроснабжения квазипостоянных токов (КПТ), которые попадают с СЭС через заземленные нейтрали силовых трансформаторов (СТ). Амплитудные значения КПТ в протяженных высоковольтных линиях электропередачи (ЛЭП) могут быть сравнимы со значениями рабочих токов ЛЭП. Квазипостоянные токи вызывают изменение режима работы СТ за счет насыщения ферромагнитного сердечника. Для учета влияния квазипостоянных токов на режим работы силового трансформатора необходимо учитывать изменение индуктивности ветви намагничивания от величины протекающего квазипостоянного тока.

Описание модели обмотки силового трансформатора

Для получения зависимости индуктивности ветви намагничивания СТ от величины КПТ рассмотрена модель обмотки СТ с бесконечным ферромагнитным сердечником и цилиндрической магнитной оболочкой. Обмотка СТ заменена круговым поверхностным током I, текущим по поверхности круговой ленты радиуса r0 и шириной 2z0. Расчёт выполнен в цилиндрической системе координат (r, ц, z) (рисунок 1).

а) — система координат; б) — расчетные размеры.

На рисунке 1 обозначены:

I — ферромагнитный сердечник с относительной магнитной проницаемостью µ1;

II — область, включающая обмотку СТ с относительной магнитной проницаемостью µ2 = 1;

III — цилиндрическая оболочка из конструкционной стали с относительной магнитной проницаемостью µ3;

IV — область учитывающая бак СТ с магнитной проницаемостью µ4 = 1; r1, r0, r2, r3 — радиусы областей I, II, III и IV.

Относительная магнитная проницаемость µ1 в модели принята постоянной при заданном токе I, т. е. все точки сердечника имеют одно и то же значение µ1. Это означает, что не учитывается неоднородность намагничивания сердечника током I.

Под током I понимается полный ток, текущий по ленте.

.

где w — число витков в обмотке.

Для моделирования процессов насыщения СТ в соответствии с [5] рассчитаем векторный потенциал кругового ленточного тока, расположенного в свободном пространстве (сердечник отсутствует). Поскольку плотность тока имеет только ц — составляющую, то и векторный потенциал имеет единственную составляющую Aц(r, z), которая не зависит от ц.

Система уравнений для компонент Aц векторного потенциала исследуемой системы: ток трансформатор намагничивание.

(1).

Коэффициенты С 1 и С 2, определяются из уравнений (2):

(2).

Коэффициенты A, B, D, E, F, G могут быть найдены из граничных условий при r = r1, r = r2, r = r3 (по 2 на каждой границе):

при r = r1, Aц1 = Aц21,.

; (3).

при r = r2, Aц22 = Aц3,.

; (4).

при r = r3, Aц3 = Aц4,.

. (5).

Обозначив.

.

.

.

.

.

.

.

подставив (2) в граничные условия (3)-(5), получим:

(6).

Из системы (6) получены выражения для определения коэффициентов А, В, D, E, F, G:

(7).

где

(8).

(9).

(10).

(11).

С учетом выражений (7) — (11) составляющая Aц(r, z) векторного потенциала исследуемой системы:

(12).

Для нахождения индуктивности обмотки силового трансформатора применим выражения векторного потенциала Aц(r0, z) для вычисления энергии магнитного поля WМ и Aц(r1 = r0, 0) — для вычисления потока Ф 1 через ферромагнитный сердечник в плоскости z = 0, а также Aц (r2, 0) и Aц (r3, 0) для вычисления потока Ф 3 через сечение магнитной оболочки в плоскости z = 0:

(13).

Для частного случая r1 = r0, когда поверхностный ток находится на поверхности сердечника:

(14).

При отсутствии оболочки (a = 1, b = 0) выражение (14) преобразуется в вид:

(15).

или при z = 0.

. (16).

Для вычисления энергии магнитного поля и потоков Ф 1 и Ф 3 введем обозначения Q3,Q4:

(17).

причём Q3 = Q3(r0, z0, м1, r2, r3, м3,) и Q4 = Q4(r0, z0, м1, r2, r3, м3,), т. е. Q3 и Q4 зависят также от параметров цилиндрической магнитной оболочки r2, r3, м3.

Из (13) с учетом (12) получим выражение для вычисления потока Ф 3.

.(18).

Введём обозначение для интеграла в выражении (18),.

. (19).

Аналогично для Q3, Q4 интеграл Q5 = Q5(r0, z0, м1, r2, r3, м3,).

Поток Ф 3 через сечение оболочки в плоскости z=0 из выражения (18):

(20).

откуда.

. (21).

Таким образом, получены выражения для магнитной индукции и индуктивности обмотки от тока намагничивания:

(22).

(23).

. (24).

В формулах (21) — (23) м1 = м1(B), м3 = м3(B), B1 = B1(I), B3 = B3(I).

Зависимость м1 = м1(B1) = f4,1(B1) для сердечника силового трансформатора получена на основе кривой намагничивания электротехнической стали. Для получения зависимости м3 = м3(B) = f4,3(B3) в случае цилиндрической магнитной оболочки рассмотрена кривая намагничивания конструкционной стали, которая отличается от кривой для электротехнической стали.

Из (21), (22) ток в ленте:

(24).

где ц3 = ц3(B1)= B3(B1).

Зависимость B1(I) может быть получена путём нахождения обратной функции к ц3:

(25).

В результате зависимость индуктивности от тока в ленте:

. (26).

Индуктивность L, отнесённая к току одного витка равна:

(27).

Результаты моделирования

Для рассмотренной модели обмотки СТ с бесконечным ферромагнитным сердечником и цилиндрической магнитной оболочкой был проведен расчет зависимости индуктивности обмотки от величины тока намагничивания. Расчет проводился для обмотки с числом витков w = 200 и радиусом круговой ленты r0 = 0,5 м. На рисунке 2 приведены две зависимости индуктивности обмотки СТ от тока намагничивания. Кривая 1 учитывает влияния бака СТ (с оболочкой), кривая 2 не учитывает влияние бака силового трансформатора на индуктивность обмотки.

Анализ полученных зависимостей показывает, что в области токов намагничивания от 0 до 100 А наблюдается различия в значениях индуктивностей обмотки до 20% без учета и с учетом оболочки из конструкционной стали. В области токов намагничивания выше 100 А различия практически не наблюдаются.

Рис. 2. — Зависимость L = L(i0) для обмотки с бесконечным сердечником и цилиндрической магнитной оболочкой при радиусе круговой ленты r0 = 0,5 и числе витков обмотки w = 200 (1 — с оболочкой, 2 — без оболочки)

Заключение

Рассматриваемая модель обмотки силового трансформатора позволяет учитывать влияние КПТ на режим его работы. Сравнение полученных зависимостей индуктивностей обмоток от тока намагничивания для модели с учетом и без учета бака показало, что наибольшие отличия, порядка 20%, наблюдаются при значениях тока намагничивания меньше 100 А. При значениях тока намагничивания выше 100 А отличия в характеристиках L(i0) не значительны, что объясняется насыщением цилиндрической магнитной оболочки и уменьшением её магнитной проницаемости. Поэтому при моделировании процессов насыщения в присутствии квазипостоянных токов бак силового трансформатора необходимо учитывать только при токах намагничивания величиной меньше 100 А.

• 1. Вахнина, В.В., Черненко, А.Н., Кузнецов, В. А. Влияние геоиндуцированных токов на насыщение магнитной системы силовых трансформаторов // Вектор науки Тольяттинского государственного университета. — 2012. — № 3(21). — С. 65−66.
• 2. Вахнина, В.В., Кретов, Д. А. Определение допустимых уровней геоиндуцированных токов для обеспечения работоспособности силовых трансформаторов при геомагнитных бурях // Интернет-журнал «НАУКОВЕДЕНИЕ». — 2012. — № 3. URL: naukovedenie.ru/sbornik12/12−93.pdf.
• 3. Бадамс, А.М., Сомов, В.А., Шмидт, А. О. Трансформаторы и стабилизаторы, регулируемые подмагничиванием шунтов. — М.: Госэнергоатомиздат, 1959. — 159 с.
• 4. Дружинин, В. В. Магнитные свойства электротехнической стали. — М.: Энергия, 1974. — 240 с.
• 5. Пановский, В., Филипс, М. Классическая электродинамика. — М.: Физматизд, 1963. — 432 с.
• 6. Зихерман, М.Х., Кузьмин, Г. П., Лейтес, Л. В. Магнитная характеристика электротехнической стали при сильном насыщении // Электротехническая промышленность. Аппараты высокого напряжения, трансформаторы, силовые конденсаторы. — 1972. № 7(16). — С.3−5.
• 7. Лейтес, Л. В. Электромагнитные расчёты трансформаторов и реакторов. — М.: Энергия, 1981. — 392 с.
• 8. Розенблат, М. А. Магнитные элементы автоматики и вычислительной техники. — М.: Наука, 1974. — 768 с.
• 9. Aubin, J., 1992. Effect of geomagnetically induced currents of power transformers. Electra, 141: pp. 24 — 33.
• 10. Lahtinen, M. and J. Elovaara, 2002. GIC occurrences and GIC tests for 400 kV system transformer. IEEE Transactions on Power Delivery, 17: pp. 555−561.
• 11. Кралин А. А., Охотников М. Н. Моделирование режимов параллельной работы трехфазных распределительных трансформаторов // Инженерный вестник Дона, 2014, № 4 URL: ivdon.ru/ru/magazine/archive/N4y2014/2602.
• 12. Папков Б. В. Основа безопасности государства — электроэнергетика // Инженерный вестник Дона, 2013, № 4 URL: ivdon.ru/ru/magazine/archive/n4y2013/2023.

References

• 1. Vakhnina, V.V., Chernenko, A.N., Kuznetsov, V.A. Vliyaniye geo indutsirovannykh tokov na nasyshcheniye magnitnoy sistemy silovykh transformatorov [Influence of geo induced currents at saturation magnetic system of power transformers]. Vektor nauki of Togliatti State University. 2012. № 3(21). pp. 65−66.
• 2. Vakhnina, V.V., Kretov, D.A. Opredeleniye dopustimykh urovney geo indutsirovannykh tokov dlya obespecheniya rabotosposobnosti silovykh transformatorov pri geomagnitnykh buryakh [Definition of admissible levels of the geoinduced currents for ensuring operability of power transformers at geomagnetic storms]. Internet — zhurnal «NAUKOVEDENIYE». 2012. № 3. URL: naukovedenie.ru/sbornik12/12−93.pdf
• 3. Badams, A.M., Somov, V.A., Shmidt, A.A. Transformatory i stabilizatory, reguliruyemyye podmagnichivaniyem shuntov (Rus) [Transformers and stabilizers, adjustable magnetization shunts.]. M.: Gosenergoatomizdat, 1959. pp. 159 .
• 4. Druzhinin, V.V. Magnitnyye svoystva elektrotekhnicheskoy stali (Rus) [The magnetic properties of electrical steel]. M .: Energiya, 1974. pp. 240
• 5. Panovskiy, V., Filips, M. Klassicheskaya elektrodinamika (Rus) [Classical lectrodynamics]. M.: FizmatIzd 1963 g. 432 p.
• 6. Zikherman, M.H., Kuz’min, G.P., Leytes, L.V. Magnitnaya kharakteristika elektrotekhnicheskoy stali pri sil’nom nasyshchenii (Rus). Elektrotekhnicheskaya promyshlennost'. Apparaty vysokogo napryazheniya, transformatory, silovyye kondensatory. 1972. № 7 (16). pp. 3−5.
• 7. Leytes, L.V. Elektromagnitnyye raschoty transformatorov i reaktorov (Rus) [The electromagnetic calculations of transformers and reactors]. M.: Energiya, 1981. pp. 392
• 8. Rozenblat, M.A. Magnitnyye elementy avtomatiki i vychislitel’noy tekhniki (Rus) [Magnetic elements of automation and computer technology]. M.: Nauka, 1974. pp.768.
• 9. Aubin, J., 1992. Effect of geomagnetically induced currents of power transformers. Electra, 141: pp. 24 — 33.
• 10. Lahtinen, M. and J. Elovaara, 2002. GIC occurrences and GIC tests for 400 kV system transformer. IEEE Transactions on Power Delivery, 17: pp. 555−561.
• 11. Kralin A.A., Okhotnikov M.N. Inћenernyj vestnik Dona (Rus), 2014, № 4 URL: ivdon.ru/ru/magazine/archive/N4y2014/2602.
• 12. Papkov B.V. Inћenernyj vestnik Dona (Rus), 2013, № 4 URL: ivdon.ru/ru/magazine/archive/n4y2013/2023.