CC BY
143
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Моделирование переходных процессов и анализ динамической
устойчивости синхронных генераторов при воздействии
больших возмущений
# 10, октябрь 2012
Б01: 10.7463/1012.0465566
Енин В. Н., Степанов А. В.
УДК [621.31:63]+631.001.57
Россия, МГТУ им. Н.Э. Баумана [email protected]
Обеспечение качественного управления элементами электроэнергетических систем, в частности синхронными генераторами, при больших возмущениях является необходимым условием их надежного функционирования и живучести. Большие возмущения в сочетании с несовершенством системы противоаварийной автоматики и конструкции генераторов, гидротурбин и другого оборудования могут приводить к авариям с катастрофическими последствиями. Вследствие этого для надежной работы и настройки противоаварийной автоматики требуется анализ динамической устойчивости синхронных генераторов при воздействии различных возмущений. В настоящее время разрабатываются методы анализа динамической устойчивости, основанные на методе функций Ляпунова. Обзору работ по методам анализа динамической устойчивости с использованием функций Ляпунова посвящена работа [1]. В работе основное внимание уделяется методу основанному на сведении многомашинной системе к одномашинному эквиваленту. Однако правомочность такого упрощения, сведения многомашинной системы к одномашинному эквиваленту не достаточно обосновано.
Более обоснованными представляются методы, использующие более полную модель системы: функцию Ляпунова энергетического типа для консервативной модели многомашинной системы [1, 2, 3]. В работах [2, 3] предложена методика анализа запаса динамической устойчивости использующая переходные траектории движения, полученные численным интегрированием и метод функций Ляпунова. Такой подход позволяет расчет переходных траекторий производить с помощью численных методов и достаточно подробных математических моделей системы [4], а для оценки запасов динамической устойчивости упрощенные консервативные модели и метод функций
Ляпунова. После расчета траекторий движения производится аппроксимация функции Ляпунова полиномом для заданных возмущений и определяется критическое значение.
1. Математическая модель электроэнергетической системы
Рассмотрим модель системы, которая включает синхронные генераторы, работающие на электрическую сеть. К узлам электрической сети с номерами к = 1,
2, 3,......, п - подключены синхронные генераторы. Упрощенная модель генератора
представляет ЭДС Ек (действующее значение) с фазой 6к включенной за реактивным сопротивлением 1/Ьк. Узлы примыкания генераторов к сети имеют номера к+п, где к -номер узла подключения генератора. К узлам с номерами 2п+1, 2п+2, ... , 2п+т подключаются нагрузки. В 5-ом узле сети напряжение задается действующим значением и и фазой е,. К узлу сети может быть подключена комплексная нагрузка Рн¡i+jQн¡i. Уравнения для такой системы имеют вид
( 20 —
Тк—к = = рк -ЬкЕк^к+п в1и(9к - 0к+п), к = 1, п, (1)
Ж2
п+т
Рн* + Ь^иЕ™ ^п(9.,я - е,) + X Ькимвтфк - е,) = 0, (2)
к=п+1
п+т
Qнs - ь,-пи, - и,Е,-пСов(е,-п - е,)) - X Ь,1(и - ии/^е/ - е,) = о, , = п+1,п+т, (3)
I =п+1
где Тк - постоянная времени к-го генератора; Рк - мощность турбины к-го генератора; Ь,к -реактивная проводимость ветви электрической сети между узлами , и к.
Для исследования взаимных возмущенных движений роторов синхронных
генераторов производим в уравнениях (1 - 3) замену переменных 5к = ек - ес, к = 1, п + т (где угол ес = Тс-1 Xп^Т^ь - Центр инерции системы) получаем
Тк = = Рк - Тк ^ - ЬкЕкик+пЫп(&к - бк+и), к = 1, п, (4)
Ж2 Тс
п +т
и^п (б,-„ - 5,) + X Ьии^Цб/ - б,) = 0, (5)
/ =п+1
п+т _
Qнs - Ь.„(и2 - и,Е,-„81п(б,-„ - б,)) - X (Ьи2 - Ь,/и,и/$т(б/ - б,)) = 0, , = п +1, п + т. (6)
/ =п+1
Нетрудно показать, что для системы (4) - (6) выполняется равенство
А
(Рк - Тк
к=1
Е (Рк - Тк-) = Е
Т V= п+1
и она является консервативной системой, поскольку мощность, вводимая в систему, равна мощности потребляемой нагрузкой.
Функция Ляпунова для консервативной модели системы (4) - (6) имеет вид
К(5,ю, и) = К(ю) + Ж(5, и), (7)
где
К(ю) = ^(ш) - К1(ю ), Ж(6,Ц) = ^1(5, и) - ^(5 *,и*).
Г1(б,ю,и) = ^(Ю) + ^1(5, и),
п лг
К:(®) = 0,5]ЕТ Юк = ^
^1(5, и) = + Жь + Жс,
(8) (9) (10)
(11) (12)
п Т \ \ п п+т
жа = -Е Рк-Т- §к-ЕЬкЕкип+ксов(5к- +п) + Е Рн*г,
к=1 ^ Тс у к=1 V =п+1
^=о,5Е \и1+п -ЕЕ кии со^г - ^),
к=1
V = п+11 = V+1
п+т п+т
Жс = о,5Е Еьи-Е I
О, (и*)
С ' * } Т Т
V = п+1 1 = V+1 V = п+1 0 V
-ли .
А =Еп Р + Еп+т р . (13)
¿—¡к=1 к ¿—¡ь =п+1 4 7
В уравнениях [7 - 10] обозначено: (5,и) - вектор переменных (51, 52, 53, ..., 5п, ип+1, Цп+2, ..., ип+т, 5п+1, 5п+2,, ..., 5п+т); 51, 52, 53, ..., 5п - фазы ЭДС генераторов; 5п+1, 5п+2,, ..., 5п+т - фазы напряжений в узлах сети относительно угла 9с; (5 ,ш ,и) - положение равновесия послеаварийной системы. В консервативной системе Д = 0.
2. Вычисление критического значения функции Ляпунова и запасов
устойчивости
При анализе динамической устойчивости вначале определяется критическое значение функции Ляпунова Жкр. Запас динамической устойчивости вычисляется по формуле
Ж, - г
К Д =
кр
ж
кр
где величина Гвозм вычисляется в момент устранения возмущения (¿кз) и равна V возм = V (^(¿кз ),®(^кз (¿кз)). Для поиска критического значения функции Ляпунова на
траектории движения определяется точка максимума функции Ж(5,Ц). Точки максимума
*
обозначим Жт. Через точку 5 и точку максимума в многомерном пространстве углов проводится прямая на которой в вычисляются значения функции Ж(5,Ц) [2]. Для более быстрого нахождения критического значения целесообразно использовать полиномиальную аппроксимацию функции Ж(5($,),и($>)). В качестве параметра £ можно использовать расстояние точек локальных максимумов функции Ж от точки положения равновесия (или время отключения возмущения)
£
V
I(ы-8* )2.
Пусть заданы значения параметра ¿а и соответствующие им значения функции Ж на прямой в многомерном пространстве
£0, £1, £2, £3, ....., £К+1,
Ж(£о), Ж(£1), Ж(£2), Ж(£з), ..., Ж(£к+1). Представим интерполяционный полином, проходящий через К+1 точку функции У(£), в виде
Р (£)=!/адИ ^^. (14)
к =0 г =0 (£к £г)
г Фк
Аппроксимируем функцию Ж(5(£)) локальных максимумов заданную в К+2 точках полиномом вида (14). Тогда для наилучшей аппроксимации должно выполняться
I(ж £)+(-1)к*)П= (ж (£к+,)+(-1)К
к =0 У ' г =0 (£к £г) У '
г Фк
Погрешность аппроксимации в заданных точках будет равна
¿Ж Су,)П^^-Ж (Ук+1)
к=о г=0 (Ук уг)
_ i *к
¿» „ ., К , I
(-1)к ±'-:с (-1)к ±1Л у-у
к =0 г =0 (Ук Уг)
i Фк
Полином, аппроксимирующий функцию Щ(у), тогда будет иметь вид
жл )(Ж (у ,)+(-1ЙП ^^ . (15)
к-И ' г =0 (Ук -
г Фк
Критическое значение функции Ляпунова определим из условия
СЩ.(у) _ 0
су 0
Дифференцируя выражение (15) по у и учитывая, что
СПП (У - Уг) ТГ (у - У)
гФк 1фк гФк
,ф1
получаем уравнение для определения Укр при котором функция Щ.(у) принимает максимальное значение:
К , К К ( )
ЕЖ^ЕПу-^ _ 0, (16)
к = 0 1 =0 г_0 (У к У г)
1 Фк г Фк
1Ф1
где Щ.(ук) _ Щ Су к)+(-1)*.
Для повышения точности при определении критического значения следует выбирать полином более высокой степени.
В качестве примера рассматривалась трехмашинная электроэнергетическая система, в которой три генератора через линию электропередачи работают на общую нагрузку. Аварийный процесс вызывается трехфазным КЗ на линии 2. Возмущение вызывает переходный процесс, который может привести к нарушению синхронной работы генераторов. Возмущение деблокируется отключением КЗ и одной цепи линии. По результатам расчетов при различных длительностях КЗ были получены следующие значения максимумов функции ЩЩ/а) на траектории
У0 = 0.0, У1 = 0.365, У2 = 0.945, У3 = 0.9625, У4 = 0.979, У5 = 1.095; Жт(У0) = 0, Щт(У1) = 0.0115, Щт(У2) = 0.04, Щт(Уз) = 0.0491, Щт(У4) = 0.0495, Щ^) = 0.029.
С использованием выражений (15), (16) (использовался полинома четвертой степени) было вычислено критическое значение функции Ляпунова равное Щкрит = 0.0548
и критическое время деблокирования ^крит=0.1648. Результаты расчета запаса динамической устойчивости по этим данным приведены в таблице.
Таблица. Запасы устойчивости при различных длительностях КЗ
¿КЗ, с ^озм Гкр Кд, %
0,05 0,0040 0,0548 92
0,10 0,0115 -» 79
0,15 0,0400 -» 27
0,16 0,0495 -» 9
0,17 0,0605 -» -10
В таблице обозначено: ¿КЗ - длительность короткого трехфазного замыкания; Квозм -значение функции Ляпунова после деблокирования возмущения; Жкрит - критическое значение функции Ляпунова; Кд - коэффициент запаса устойчивости (в процентах).
Особо опасными режимами для мощных синхронных гидрогенераторов считаются короткие трехфазные замыкания на шинах, которые могут приводить к выделению большой мощности и разрушению гидроагрегата. Далее рассматривается динамическая устойчивость одного из мощных синхронных генераторов 6400 МВт работающих параллельно на линию передач 525 кВ. Динамическая устойчивость исследовалась при воздействии возмущения в виде КЗ на шинах генератора. Исследование проводилось с использованием системы моделирования динамических систем MATLAB. Запасы устойчивости и критическое время отключения возмущения определялось с использованием предлагаемой методики. Параметры синхронного генератора СВФ 1285/275-42 УХЛ4:
£н = 711 МВА; ин =15.75 кВ; щ = 142.8 об/мин; J = 25400 тм
xs =0.184; Xd = 1.58; xd = 0.43; xd' = 0.295; xq = 0.97; xq = 0.31; (о. е.)
Tdo = 8.21 с; Td0 = 0.143 с.
Параметры трансформатора:
& = 711 МВА; инн =15.75 кВ; иш =525 кВ;
ДРх = 500 кВт; Ix = 0.35 %; APK = 1300 кВт; uK =14 %.
Структурная схема в Simulink с использованием блоков расширения SimPoweSystems приведена на рис. 1.
Рис 1. Структурная схема в системе МЛТЬЛБ мощного синхронного генератора
работающего на линию 525 кВ
Проводилось моделирование системы (рис. 2) при различных длительностях КЗ на шинах синхронного генератора. Результаты моделирования: напряжение на шинах генератора (фаза С); фаза ЭДС генератора для случая нарушения динамической устойчивости приведены на рис. 2 - 3.
г Ус(ри) □0 X
ц
(
1
0 1 2 3 4 5 Т ¡гше оИ:5е!:: 0
Рис. 2. Напряжение статора фазы С (КЗ на шинах генератора, о. е.).
-20-1-1-1-1-
0 1 2 3 4 5
Time offset: 0
Рис. 3. Изменение фазы ЭДС генератора относительно синхронной оси при КЗ на шинах
генератора (потеря устойчивости), о. е.
При различных длительностях КЗ на шинах проводился расчет переходных процессов и анализ динамической устойчивости синхронного генератора. Анализ проводился с использованием функций Ляпунова энергетического типа. Проведенные расчеты были использованы для аппроксимации функции Wm в зависимости от длительности КЗ. Затем было найдено критическое время отключения возмущения (¿крит~0.105 с).
Предлагаемая и развиваемая методика анализа динамической устойчивости синхронных генераторов совмещает преимущества методов численного моделирования и качественных методов (методов функций Ляпунова). Применение этой методики предполагает вычисление функции Ляпунова для всех точек многомерной траектории переходного процесса, получаемой численным интегрированием уравнений динамики системы. Это позволяет получать не только траектории движения при воздействии возмущения, но и производить некоторую интегральную оценку динамической устойчивости, для более эффективной настройки системы протиаварийного автоматического управления синхронными генераторами.
Список литературы
1. Павелла М. От общей теории Ляпунова к практическому прямому методу анализа
динамической устойчивости энергосистем // Электричество. 2000. № 6. C. 14-26.
2. Авраменко В.Н., Степанов А.В. Определение запаса динамической устойчивости энергосистем прямым методом Ляпунова // Техническая электродинамика. 1999. № 5. С. 55-58.
3. Кузовкин В.А., Степанов А.В. Оценка запаса динамической устойчивости энергосистем прямым методом Ляпунова // Электричество. 2002. № 1. С. 2-8.
4. Строев В.А., Шульженко С.В. Математическое моделирование элементов электрических систем : курс лекций. М.: Издательство МЭИ, 2002. 56 с.
SCIENTIFIC PERIODICAL OF THE RAIJMAN MS TU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Simulation of transients and analysis of dynamic stability of
synchronous generators under the influence of large perturbations
# 10, October 2012
DOI: 10.7463/1012.0465566
Enin V., N., Stepanov A.V.
Russia, Bauman Moscow State Technical University
To ensure reliable operation of electric power systems it's necessary to analyze and evaluate the inventory of dynamic stability of synchronous generators under large perturbations. In this paper the authors consider a methodology for analysis and evaluation of dynamic stability a feature of which is sharing high quality (the method of Lyapunov functions) and numerical methods. The authors propose an approximation method with the use of the interpolation polynomial to reduce computational costs and to carry out more accurate calculation of the critical value of the Lyapunov function. This method can improve accuracy of results to assess dynamic stability and make a more grounded choice of states of power systems and configuration of emergency automation systems.
Publications with keywords:step response process, synchronous generator, simulate power systems, transient stability, Lyapunov function
Publications with words:step response process, synchronous generator, simulate power systems, transient stability, Lyapunov function
References
1. Pavella M. Ot obshchei teorii Liapunova k prakticheskomu priamomu metodu analiza dinamicheskoi ustoichivosti energosistem [From the general theory of Lyapunov to practical direct method of analysis of the dynamic stability of power systems]. Elektrichestvo [Electricity], 2000, no. 6, pp. 14-26.
2. Avramenko V.N., Stepanov A.V. Opredelenie zapasa dinamicheskoi ustoichivosti energosistem priamym metodom Liapunova [Determination of reserve of dynamic stability of power systems by use the direct method of Lyapunov]. Tekhnicheskaia elektrodinamika [Technical electrodynamics], 1999, no. 5, pp. 55-58.
3. Kuzovkin V.A., Stepanov A.V. Otsenka zapasa dinamicheskoi ustoichivosti energosistem priamym metodom Liapunova [Reserve assessment of dynamic stability of power systems by use the direct method of Lyapunov]. Elektrichestvo [Electricity], 2002, no. 1, pp. 2-8.
4. Stroev V.A., Shul'zhenko S.V. Matematicheskoe modelirovanie elementov elektricheskikh sistem [Mathematical modeling of the elements of electrical systems]. Moscow, MEI Publ., 2002. 56 p.
