Электронный научный журнал
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,791

МОДЕЛИРОВАНИЕ ОДНОМАШИННОЙ ЭНЕРГОСИСТЕМЫ С ПОМОЩЬЮ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ В MATLAB SIMULINK

Кулешова Е.О. 1 Шишка Н.В. 1
1 ФГБОУ ВПО «НИ ТПУ» («Национальный исследовательский Томский политехнический университет»)
Проведено исследование статической устойчивости одномашинной энергосистемы, синхронный генератор которой оборудован автоматическим регулятором возбуждения сильного действия (АРВ СД). В качестве возмущения, вызывающего переходный процесс в энергосистеме, было принято отклонение активной мощности генератора от номинального. Параметрами регулирования являются генераторное напряжение и частота вращения ротора генератора. Энергосистема описывается системой нелинейных дифференциально-алгебраических уравнений, для упрощения которой используется метод линеаризации. В среде Matlab Simulink было смоделировано решение линеаризованной системы дифференциальных уравнений, описывающих энергосистему, с помощью передаточных функций. Правильность моделирования была подтверждена расчетами корней характеристического уравнения. В результате моделирования были получены графики по отклонению частоты и напряжения от номинальных значений. Выбраны коэффициенты усиления АРВ по отношению напряжения, по первой и второй производным напряжения.
моделирование
коэффициенты усиления
статическая устойчивость
частота
генератор
напряжение
мощность
энергосистема
переходный процесс
1. Веников В.А. Теория подобия и моделирования (применительно к задачам электроэнергетики). - Изд. 2-е, доп. и перераб. – М. : Высшая школа, 1976. – 479 с.
2. Ерофеев А.А. Теория автоматического управления : учебник для вузов. – 2-е изд., перераб. и доп. – СПб. : Политехника, 2005. – 302 с. : ил.
3. Лазарев Ю. Моделирование процессов и систем в MATLAB. Учебный курс. – СПб. : Питер ; Киев : Издательская группа BHV, 2005. – 512 с. : ил.
4. Справочник по проектированию электрических сетей / под ред. Д.Л. Файбисовича. – 2-е изд., перераб. и доп. – М. : Изд-во НЦ ЭНАС, 2006. – 352 с. : ил.
5. Хрущев Ю.В. Методы расчета устойчивости энергосистем : учебное пособие. – Томск : Изд-во ТПУ, 2005. – 176 с.
6. Хрущев Ю.В., Заподовников К.И., Юшков А.Ю. Электромеханические переходные процессы в электроэнергетических системах : учебное пособие. – Томск : Изд-во ТПУ, 2010. – 160 с.
7. Shishka N.V., Kuleshova E.O. The description of the one-machine power supply system by means of transfer function // Vedecky Pokrok na Prelomu Tysyachalety : materialy VIII Mezinarodni Vedecko-Prakticka Konference, Prague, 27 May-5 June 2012. – Praha : Education and Science, 2012. - Vol. 29. Technicke vedy. Telovychova a sport. - Р. 8-10.

Введение

Устойчивость системы характеризует способность системы возвращаться в состояние равновесия после исчезновения внешних сил, которые вывели ее из этого состояния. Следовательно, только устойчивая система является работоспособной.

На все переходные процессы в энергосистемах действуют возмущения – начальные отклонения от номинальных параметров режима. Это могут быть отклонения значений токов, напряжений, мощностей и других параметров.

Для исследования статической устойчивости энергосистемы необходимо составить систему дифференциальных уравнений переходных процессов. Решение такой системы является сложной задачей. Для упрощения дифференциальных уравнений, описывающих энергосистему, был использован метод линеаризации «в малом» [5; 6].

В среде Matlab Simulink было смоделировано решение линеаризованной системы дифференциальных уравнений.

Для наглядности расчета коэффициентов характеристического уравнения зададимся исходным данными для схемы, представленной на рис. 1 [5].

Математическая модель энергосистемы

Для описания энергосистемы с помощью передаточных функций рассмотрим одномашинную энергосистему (рис. 1), синхронный генератор которой оборудован автоматическим регулятором возбуждения сильного действия (АРВ СД). Система дифференциальных уравнений переходных процессов для данной энергосистемы имеет вид [1]:

,

,

,

.

Линеаризация функции независимой переменной уравнений (1-4) приводит к производной малого приращения этой функции [5]. В результате получим линейную систему вида:

.

Линеаризация уравнения (5) проводится относительно угла δ и синхронной ЭДС Eq.

Приращение угла δ пропорционально частоте вращения ротора генератора:

,

,

,

,

.

В совокупности уравнения (5-10) составляют математическую модель энергосистемы.

 

Рис. 1. Схема одномашинной энергосистемы

Для рассматриваемой одномашинной энергосистемы в качестве варьируемого параметра режима была принята активная мощность генератора. Параметрами регулирования являются модуль UГ вектора генераторного напряжения UГ и частоту f вращения ротора генератора.

Параметры схемы

Примем следующие типы электрооборудования и параметры [4]:

Генератор Г является эквивалентом двух турбогенераторов типа ТГВ-200-2УЗ. Трансформатор Т1 является эквивалентом двух трансформаторов типа ТДЦ-250000/220.

Трансформатор Т2 является эквивалентом двух групп однофазных трансформаторов типа АТДЦТН-250000/220/110.

Линия ВЛ – двухцепная, выполнена проводом марки АСО-300 с параметрами: длина , .

Возбудитель генератора: Те=0,04 с. Регулятор возбуждения: Тр = 0,04 с.

Система С: Uc=115 кВ.

Турбина эквивалентного генератора: режим: Рт(1) = 0,5 Рmax.

Коэффициенты усиления АРВ: К0U=50 ед. возб. хх*с/ед. напряж., К1U=5 ед. возб. хх*с/ед. напряж.

Расчет коэффициентов

Для исследования установившегося режима энергосистемы запишем формулы частных производных в уравнениях:

 

 

.

Моделирование в среде Matlab Simulink

На основе выше полученных расчетов построим одномашинную энергосистему в среде Matlab Simulink [3; 7].

 

Рис. 2. Схема одномашинной энергосистемы, представленная в среде Matlab Simulink

Фиксируем коэффициенты усиления АРВ К0U и К1U по отклонению напряжения и по первой производной напряжения соответственно. Коэффициенты усиления АРВ К0f и К1f по отклонению частоты и по первой производной частоты должны попадать в область устойчивости моделируемой энергосистемы. Обычно для определения области устойчивости при проектировании энергосистем с двумя параметрами, влияющими на устойчивость энергосистемы, используют метод D-разбиений [2]. Процедура построения области устойчивости является трудоемкой. Более того, необходимо построение нескольких областей устойчивости при различных режимах работы генератора.

Используя программные средства MatLab Simulink, настройку коэффициентов усиления АРВ К0f и К1f так, чтобы система была устойчивой, можно выполнить, не прибегая к построению области устойчивости. Выполнить это можно с помощью построения корневого годографа (рис. 3).

Рис. 3. Использование блока SISO Design Task при настройке коэффициентов К0f и К1f

В результате моделирования коэффициенты усиления были приняты равными К0f =0,58 и К1f =2,726. В качестве проверки нами была построена область устойчивости для данной энергосистемы методом D-разбиений. Найденные коэффициенты попадают в полученную область устойчивости.

При полученных коэффициентах усиления АРВ К0f и К1f корни характеристического уравнения располагаются в левой комплексной полуплоскости, как это видно на рис. 5. В качестве проверки сравним полюса передаточных функций с расчетами корней характеристического уравнения в MathCad при тех же коэффициентах усиления АРВ.

 

Рис. 4. График линейного анализа, показывающий корни характеристического уравнения

Из рис. 3 видно, что отклонение напряжения от заданного значения стремится к нулю, следовательно, система устойчива. Время переходного процесса в среднем составляет 10 с.

Заключение

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

Рецензенты:

Лукути Б.В., д.т.н., профессор, зав. кафедрой ЭПП ЭНИН ФГБОУ ВПО «НИ ТПУ», Национальный исследовательский Томский политехнический университет, г. Томск.

Новиков С.А., д.ф.-м.н., профессор кафедры ЭСиЭ ЭНИН ФГБОУ ВПО «НИ ТПУ», Национальный исследовательский Томский политехнический университет, г. Томск.


Библиографическая ссылка

Кулешова Е.О., Шишка Н.В. МОДЕЛИРОВАНИЕ ОДНОМАШИННОЙ ЭНЕРГОСИСТЕМЫ С ПОМОЩЬЮ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ В MATLAB SIMULINK // Современные проблемы науки и образования. – 2013. – № 2.;
URL: http://science-education.ru/ru/article/view?id=8744 (дата обращения: 22.11.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074