Научная статья на тему 'Математическое моделирование неявнополюсного синхронного генератора мини-ТЭЦ'

Математическое моделирование неявнополюсного синхронного генератора мини-ТЭЦ Текст научной статьи по специальности «Энергетика»

CC BY
199
27
Поделиться
Область наук
Ключевые слова
НЕЯВНОПОЛЮСНЫЙ СИНХРОННЫЙ ГЕНЕРАТОР / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / УРАВНЕНИЯ ПАРКА-ГОРЕВА / ЧАСТИЧНЫЕ ХАРАКТЕРИСТИКИ НАМАГНИЧИВАНИЯ / ПЕРЕХОДНЫЕ ПРОЦЕССЫ / NONSALIENT-POLE SYNCHRONOUS GENERATOR / IMPLICIT-POLE GENERATOR / MATHEMATICAL MODEL / PARK-GOREV DIFFERENTIAL EQUATION / MAGNETIC CHARACTERISTICS / TRANSIENTS

Аннотация научной статьи по энергетике, автор научной работы — Вагин Г.Я., Солнцев Е.Б., Мамонов А.М., Петров А.А.

Предложена математическая модель неявнополюсного синхронного генератора мини-ТЭЦ на базе полных уравнений Парка-Горева, представленных во взаимной системе относительных единиц. Учет насыщения осуществляется с помощью метода частичных характеристик намагничивания неявнополюсной машины. Данный метод позволяет наиболее точно оценить изменение параметров схемы замещения неявнополюсного генератора при отсутствии достоверных сведений о магнитных свойствах материалов, применяемых при проектировании и изготовлении машины. Реализация математической модели осуществляется в среде Simulink. Применение модели возможно при расчете установившихся и переходных режимов в распределительных сетях, содержащих неявнополюсные синхронные генераторы.

Похожие темы научных работ по энергетике , автор научной работы — Вагин Г.Я., Солнцев Е.Б., Мамонов А.М., Петров А.А.,

Mathematical model of implicit-pole synchronous generator of mini-CHP

In this paper we propose mathematical model of implicit-pole synchronous generator of mini-CHP. The model is based on a full Park-Gorev differential equation. The saturation of the machine is taken by using the method of partial magnetic characteristics. This method allows most accurate evaluation of the change of parameters of an equivalent circuit of the nonsalient-pole generator in the absence of information about the magnetic properties of the machine. The model is Implemented in Simulink. Application of the model is possible when calculating the steady-state and transients in distribution electrical net wich contains implicit-pole synchronous generators.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Текст научной работы на тему «Математическое моделирование неявнополюсного синхронного генератора мини-ТЭЦ»

Математическое моделирование неявнополюсного синхронного генератора мини-ТЭЦ

Г.Я. Вагин, Е.Б. Солнцев, А.М. Мамонов, А.А. Петров Нижегородский государственный технический университет им. Р. Е. Алексеева

Аннотация: Предложена математическая модель неявнополюсного синхронного генератора мини-ТЭЦ на базе полных уравнений Парка-Горева, представленных во взаимной системе относительных единиц. Учет насыщения осуществляется с помощью метода частичных характеристик намагничивания неявнополюсной машины. Данный метод позволяет наиболее точно оценить изменение параметров схемы замещения неявнополюсного генератора при отсутствии достоверных сведений о магнитных свойствах материалов, применяемых при проектировании и изготовлении машины. Реализация математической модели осуществляется в среде Simulink. Применение модели возможно при расчете установившихся и переходных режимов в распределительных сетях, содержащих неявнополюсные синхронные генераторы.

Ключевые слова: неявнополюсный синхронный генератор, математическая модель, уравнения Парка-Горева, частичные характеристики намагничивания, переходные процессы.

Скорость вращения роторов агрегатов мини-ТЭЦ на базе газотурбинных установок обычно составляет 3000 об/мин, поэтому они комплектуются неявнополюсными синхронными генераторами [1]. Особенности регулирования скорости вращения газовых турбин [2], а также зависимость их КПД от загрузки делают целесообразным применение ГТУ в базовой части суточного графика нагрузки системы электроснабжения.

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

При моделировании синхронных генераторов предпочтение следует

отдавать моделям, основанным на системе дифференциальных уравнений Парка-Горева [4]. Во взаимной системе относительных единиц данная система имеет вид (символ «*» здесь и далее опущен):

ий = № - Уд (1 + б) - ; Юс

ид =--РУд + Уd (1 + )-^д

Ю 1

с

и1 = — РУ/ + 1

с

^ =—pwdi + ЪМ, г=;

Юс

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(1)

Пд. =- РУд. + ЯдХ

д.

ЮГ

д{дг

г = 0...п

д

т

мт - Мэ =^рю = ^ рБ,

—с

где па, пд - количество демпферных контуров по соответствующей оси; TJ -инерционная постоянная; б - скольжение ротора генератора относительно синхронно вращающихся осей; Мт - механический момент на валу генератора, создаваемый первичным двигателем; Мэ - электромагнитный момент сопротивления на валу двигателя.

Пренебрегая влиянием насыщения на сопротивления рассеяния контуров машины [5], потокосцепления генератора можно представить:

у d = +У8d;

V д = хагд + Убд;

VI = ХМ/ + У^; (2)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

у di = xdюidi +y8d: у дг = хдю1д. +У5д,

г = 0..п^ г = 0...п

д

где у^, ^5д - проекции результирующего потокосцепления воздушного зазора на продольную и поперечную оси соответственно; х0, х/а, хдга -сопротивления рассеяния соответствующего контура.

Проекции потокосцепления воздушного зазора определяются:

г п \

ч§д = пдхад

и +'/ +Ъйг

'"=1 У

л

п

и

V '=1 У

(3)

где Пи Пд - эквивалентные коэффициенты, учитывающие насыщение стали неявнополюсной машины вдоль соответствующих осей.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Выразив из системы (2) токи и подставив их в (3), получим зависимости у^, у5д для неявнополюсного генератора:

чи +Ч/Щ

хс х/с 1=1х<Из

Г

(4)

х

Чд

+ Е-

=1Х

с '

дю

гд

где Ги, Гд - проводимости по соответствующим осям:

1

Г

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

и

1 1 п 1

+— +—+1—

пихаи хс х/с ¡=1хис

Гд

1

1

у

+1-

1

(5)

пдхад хс г=1Хдю

Наиболее полно учесть насыщение стали неявнополюсного генератора можно, используя его частичные характеристики намагничивания [5], причем, при отсутствии информации о магнитных свойствах материалов, возможно применение нормальных частичных характеристик намагничивания неявнополюсной машины [6]: магнитопровода статора и зазора Фг = /(Р1г); потока рассеяния Ф/0 = /Ру); магнитопровода ротора Ф2 = /(^2), представленных на рис. 1, а.

Ввиду линейной зависимости потока рассеяния Ф/а от МДС Г1/

п

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ф/а = к/а F1fi

обмотки возбуждения за вычетом магнитного напряжения ротора, запишем:

(6)

где к/а - коэффициент взаимосвязи между Ф/а и F\f; к/а ~ 0,0457.

к2 А к1гД

0 0,5 1,0 1,5 2,0 2,5 3,0 а)

0,80,60,40,20 -

1,61,20,8 0,4 0

А

у А

к\г =/ Ф1Г) ✓

\

h = /Ф2) _

X

0 0,4 0,8 1,2 1,6 б)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ф

>

Рис. 1. Характеристики намагничивания синхронной неявнополюсной машины: а - частичные кривые намагничивания; б - коэффициенты насыщения Взаимосвязь между магнитным напряжением F и магнитным потоком Ф характеристик магнитопровода статора и зазора и магнитопровода ротора можно представить:

F = к Ф, (7)

где к - коэффициент взаимосвязи между Ф и F (рис. 1, б), полученный из соответствующих характеристик намагничивания Ф = /F).

МДС, создаваемая обмоткой возбуждения машины, определяется: F/m = F2d + Fr1d - Fadm'; (8)

где F2d - магнитное напряжение ротора машины по оси d машины; Fr1d -результирующее магнитное напряжение статора и зазора по оси d машины; Fadm - МДС реакции якоря и демпферных контуров машины по оси d. МДС обмотки возбуждения направлена по продольной оси машины, поэтому сумма МДС контуров и магнитных напряжений по поперечной оси:

0 = F2q + Fr1q - Famq. (9)

Продольная и поперечная составляющие магнитного потока ротора:

Фъи = Фги + Фм Ф2д = Фгд + Ф/ад.

Составляющие выражений (8) и (9) определяются:

к =

айт

х

ас1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

tg ас

пи

Ь +

Л

= У

К =

адт

х.

ад

tg а0

I +

д

д

ъ

г=1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(10)

(11)

где tg а0 - тангенс угла наклона касательной к нормальной характеристике холостого хода неявнополюсной синхронной машины в нулевой точке [7]; tg а0 = 1,2.

Ф и

К = К гй Г1Ы Г1г ф

г

К = К Ф 2й 1 2 й ~ 1 2

К

1гд

Ф

К гд ± 1,.

Фг

Ф,

Ф

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

к = К 2д К2д К Фп

(12)

(13)

22 Взаимосвязь между МДС и током обмотки возбуждения во взаимной системе относительных единиц выглядит:

К/т = '/

х

ай

(14)

Подход к проектированию неявнополюсной синхронной машины предполагает расчет ее параметров в соответствии с равномерным распределением магнитной индукции по каждому отдельному элементу генератора [8]. Это позволяет записать в относительных единицах:

Фг = Уб; Фгй = УбЛ Фгд = Убд, (15)

тогда решая совместно уравнения (6) - (15), получим выражения для

проекций потокосцепления уб на соответствующие оси:

Г п \

[1 + КкГ с]

х

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ай

/

Чбй =

1 + к^к -Г с

V 2 /с

+ 'и +

'=1

йг

/ .

(16)

^ а0 ( + к1г [1 + к2к/с])

I х I г + Vг

ад\ д / • <

дг

г=1 У

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

[1 + к2 кГо] ^ tg ао ( + К [1 + к2к/а])

При определении коэффициента к2 следует учесть, что:

(17)

Ф2 =л/ Ф2а + Ф2д^

Ф 2а = Ф а [1 + КкГ а]"

Ф2д = Фгд [1 + к1гк/а]

к

Га

tg ао

х

аа

ь + £ г

к

а г

г=1 У

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Га

tg ао

х

ад

па

\

г + , -

д ¿—I дг

г=1 У

Сопоставляя (16), (17) с (3), получим выражения для эквивалентных коэффициентов, учитывающих насыщение неявнополюсного генератора:

[1 + к2к/а ]

Пд =

tg ао (к2 + к1г [1 + к2к/а])

(18)

Па = П

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

к к га 1 _. 2 1

1 + к к га

2 1 а

па

гГ + +

л

V гаг

г=1 У

(19)

БШиПпк-модель блока насыщения неявнополюсного генератора (8а1игаиопБ1оск), реализующая зависимости (5) с учетом (18) и (19), приведена на рис. 2. Математическая модель синхронного неявнополюсного генератора, реализованная в операционной среде БтиНпк в соответствии с выражениями (1) - (5), приведена на рис. 3. Блоки $>у$1етТоМасЫпе и МасЫпеТо$>у$1ет переопределяют векторы напряжения и тока из синхронно вращающейся системы координат в жестко связанную с ротором генератора систему координат и обратно [4].

Моделирование переходных процессов осуществлялось на примере предварительно нагруженного синхронного неявнополюсного генератора

Т2-6-2, работающего на изменяющуюся автономную нагрузку [9]: статическую [3, 6] в виде заданных сопротивлений или динамическую [10, 11, 12] в виде асинхронного двигателя.

Рис. 2. Математическая модель блока насыщения SaturationBlock

синхронного неявнополюсного генератора Параметры генератора [13]:

• Рном = 6 МВт, cos фном = 0,8; ишм = 6,3 кВ;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

• Xad = Xaq= 1,541; хс = 0,147; Rs = 0,042; x/c = 0,023; R/ = 0,001; хЛо = 0,216; Rd1 = 0,058; хдкт = 0,216; Rq1 = 0,058.

Параметры статической (таблица № 1) и динамической (таблица № 2) нагрузки и результаты моделирования даны в относительных единицах при номинальных условиях генератора.

Рис. 3. Структурная схема математической модели синхронного

неявнополюсного генератора Учет влияния устройств регулирования на характер переходного процесса обеспечивается применением соответствующих моделей регуляторов напряжения и первичных двигателей [3, 14]. Постоянная времени первичного двигателя принимается равной 1 с.

Таблица № 1

Параметры статической нагрузки генератора

№ расчета Исходное состояние Конечное состояние

Я, о.е. X, о.е. 2, о.е. ооб ф Я, о.е. X, о.е. 2, о.е. ооб ф

1 1,6 1,2 2 0,8 1,6 1,2 2 0,8

2 3,2 2,4 4 0,8 1,066 0,8 1,333 0,8

Таблица № 2

Параметры динамической нагрузки генератора [9]

№ расчета Исходное состояние Конечное состояние

Sном, °.е. еОБ фном М, о.е. Sном, о.е. ООБ фном М, о.е.

1 0,5 0,8 0 0,5 0,8 1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2 0,5 0,8 0,5 0,5 0,8 0,5

Примечание. М - момент сопротивления на валу динамической нагрузки относительно ее номинального момента.

Результаты моделирование переходных режимов с использованием разработанной модели (рис. 4) соответствуют реальным представлениям о протекании процессов в неявнополюсной синхронной машине.

Выводы. В среде БШиШк разработана математическая модель неявнополюсного генератора, основанная на дифференциальных уравнениях Парка-Горева. Учет насыщения неявнополюсной машины осуществляется с помощью частичных характеристик намагничивания, что позволяет повысить точность при проведении расчетов переходных режимов в электрической сети с неявнополюсными генераторами. Использование модели возможно при расчетах переходных режимов, устойчивости, настройках регуляторов напряжения и скорости вращения первичных двигателей.

а) б)

Рис. 4. - Результаты моделирования переходных процессов: а - расчет №1; б - расчет №2; Ur - напряжение на зажимах генератора; 1г - ток генератора; 1д - ток динамической нагрузки; 1с - ток статической нагрузки; sT -скольжение ротора генератора относительно синхронно вращающихся осей; 5*д - скольжение ротора эквивалентного асинхронного двигателя динамической нагрузки относительно синхронно вращающейся оси

Литература

1. Белоусенко И.В., Шварц В.Г., Шпилевой В.А. Энергетика и электрификация газовых промыслов и месторождений, Тюмень: Тюмень. 2000. 273 с.

2. Цанев С.В., Буров В.Д., Ремезов А.Н. Газотурбинные и парогазовые установки тепловых электростанций. Москва: Издательский дом МЭИ, 2009. 584 с.

3. Вагин Г.Я., Солнцев Е.Б., Мамонов А.М., Петров А. А. Математическая модель электроагрегата мини-ТЭЦ на базе явнополюсного

синхронного генератора // Известия Томского политехнического института. 2015. №8 (Т. 326). С 92-101.

4. Груздев И.А., Кадомская К.П., Кучумов Л.А., Лугинский Я.Н., Портной М.Г., Соколов Н.И. Применение аналоговых вычислительных машин в энергетических системах. Методы исследования переходных процессов. Изд. 2-е, перераб. и доп. Москва: Энергия, 1970. 400 с.

5. Иванов-Смоленский А.В. Электрические машины: учебник для вузов. В двух томах. Т. 2. Москва: Издательский дом МЭИ, 2006. 532 с.

6. Вагин Г.Я., Солнцев Е.Б., Мамонов А.М., Петров А.А. Математическая модель явнополюсного синхронного генератора мини-ТЭЦ // Инженерный вестник Дона. 2015. №2. (ч.2). URL: ivdon.ru/ru/magazine/ archive/n2p2y2015/2950.

7. Вольдек А.И. Электрические машины. Учебник для студентов высш. техн. учебн. заведений. 3-е изд. перераб. Ленинград: Энергия, 1978. 832 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

8. Извеков В.И., Серихин Н.А., Абрамов А.И. Проектирование турбогенераторов. Москва: Издательство МЭИ, 2005. 440 с.

9. Granovskiy A. V. Analysis of physical processes and operating behavior of gas turbine units with high loading of the gas turbine stage // Power Technology and Engineering. 2011. Vol 44. №6. pp. 467-470.

10. Петров А.А. Математическое моделирование динамической асинхронной нагрузки // Сборник статей II Международной научно-практической конференция «Молодой Ученый: технические и математические науки». 2015. №1(1). С. 58-65.

11. Гуревич Ю.Е., Либова Л.Е., Окин А.А. Расчеты устойчивости и противоаварийной автоматики в энергосистемах. Москва: Энергоатомиздат, 1990. 390 с.

12. Glazyrin M. V. Study of synchronous overload capacity in interloading system // Russian Electrical Engineering. 2009. Vol 82. №6. pp 324-327.

13. Неклепаев Б.Н., Крючков И.П. Электрическая часть станций и подстанций: справочные материалы для курсового и дипломного проектирования. М.: Энергоатомиздат, 1989. 608 с.

14. Хватов О.С., Дарьенков А.Б., Поляков И.С. Математическое описание алгоритма управления топливоподачей дизель-генераторной электростанции переменной скорости вращения // Инженерный вестник Дона, 2013, №3. URL: ivdon.ru/magazine/archive/n3y2013/1869.

References

1. Belousenko I.V., Shvarts V.G., Shpilevoy V.A. Energetika i elektrifikatsiya gazovykh promyslov i mestorozhdeniy [Energy and electrification of the gas fields and deposits]. Tyumen': Tyumen'. 2000. 273 p.

2. Tsanev S.V., Burov V.D., Remezov A.N. Gazoturbinnye i parogazovye ustanovki teplovykh elektrostantsiy[Gas turbine and combined cycle plant of thermal power plants]. Moskva: Izdatel'skiy dom MEI. 2009. 584 p.

3. Vagin G.Ya., Solntsev E.B., Mamonov A.M., Petrov A.A. Bulletin of the tomsk polytechnic university. 2015. №8 (T. 326). pp. 92-101.

4. Gruzdev I.A., Kadomskaya K.P., Kuchumov L.A., Luginskiy Ya.N., Portnoy M.G., Sokolov N.I. Primenenie analogovykh vychislitel'nykh mashin v energeticheskikh sistemakh. Metody issledovaniya perekhodnykh protsessov [The use of analog computers in power systems. Research of electrical transients]. 2-e izd., pererab. i dop. M.: Energiya, 1970. 400 p.

5. Ivanov-Smolenskiy A.V. Elektricheskie mashiny: uchebnik dlya vuzov [Electrical machines]. V dvukh tomakh. Tom 2. M.: Izdatel'skiy dom MEI, 2006. 532 p.

6. Vagin G.Ya., Solntsev E.B., Mamonov A.M., Petrov A.A. Inzenernyj vestnik Dona (Rus), 2015. №2 (p. 2). URL: ivdon.ru/ru/magazine/archive/ n2p2y2015/2950.

7. Vol'dek A.I. Elektricheskie mashiny. Uchebnik dlya studentov vyssh.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

tekhn. uchebn. Zavedeniy [The electric machine. Textbook for technical colleges]. 3-e izd. pererab. Leningrad: Energiya, 1978. 832 p.

8. Izvekov V.I., Serikhin N.A., Abramov A.I. Proektirovanie turbogeneratorov [ The design of turbogenerators]. Moskva: Izdatel'stvo MEI, 2005. 440 p.

9. Granovskiy A. V. Analysis of physical processes and operating behavior of gas turbine units with high loading of the gas turbine stage. Power Technology and Engineering. 2011. Vol 44. №6. pp. 467-470.

10. Petrov A.A. Sbornik statey II Mezhdunarodnoy nauchno-prakticheskoy konferentsiya «Molodoy Uchenyy: tekhnicheskie i matematicheskie nauki», 2015. №1 (1). pp. 58-65.

11. Gurevich Yu.E., Libova L.E., Okin A.A. Raschety ustoychivosti i protivoavariynoy avtomatiki v energosistemakh [ Calculations of the stability and emergency control in power systems]. Moskva: Energoatomizdat, 1990. 390 p.

12. Glazyrin M. V. Study of synchronous overload capacity in interloading system. Russian Electrical Engineering. 2009. Vol 82. №6. pp 324-327.

13. Neklepaev B.N., Kryuchkov I.P. Elektricheskaya chast' stantsiy i podstantsiy: spravochnye materialy dlya kursovogo i diplomnogo proektirovaniya [The electrical part of stations and substations]. M.: Energoatomizdat. 1989. 608 p.

14. Khvatov O.S., Dar'enkov A.B., Polyakov I.S. Inzenernyj vestnik Dona (Rus), 2013, №3. URL: ivdon.ru/magazine/archive/n3y2013/1869..