8. Setkání uživatelů programu MODES ČVUT FEL Praha 19.2. 2007 Program 1. seznámení s inovacemi v nové verzi 2.2/11 programu MODES (inovované modely GASA, alternativní modelu regulátoru buzení a systémového stabilizátoru, třívstupový externí regulátor) 2. uživatelské rozhraní MODMAN v.2.11 (automatická tvorba scénářů) 3. verifikace programu 1. INOVACE PROGRAMU MODES 1.1. 1.
Inovovaný model paroplynového cyklu - PPC Zatímco předchozí model PPC (implementovaný ve verzi 2.2.10) měl omezení v tom, že. a) b)
jedna plynová turbína (model GASA) zásobuje jeden kotel na zbytkové teplo (model HRSG), generátory obou turbín jsou připojeny na jednu přípojnici (uzlu).
Nový model tato omezení částečně překračuje tím, že: A) více plynových turbín (modely GASA) může zásobuje jeden kotel na zbytkové teplo (model HRSG), B) generátory jednotlivých turbín jsou připojeny na různé přípojnice, přičemž parní turbína musí pracovat do samostatné přípojnice. Následující obrázek ukazuje alternativní způsob modelování PPC pro případ více plynových turbín zásobujících jeden kotel: Pdod
HRSG GASA 1
Tabulka bloků GASA 1
ST_S
GASA2 ….
NS
TE 1 W1
GASA n
GASA n
TE 2 W2
ST_S HRSG
TE n Wn
Obr. 1 Uspořádání rozšířeného modelu PPC Oproti stávajícímu modelu může být výkon jednotlivých plynových turbín vyveden do samostatných a/nebo společných přípojnic. Naopak výkon parní turbíny musí být vyveden do samostatné přípojnice. Model plynové turbíny GASA je upraven tak, aby na výstupu dával nejen teplotu výfukových plynů na výstupu z turbíny, ale i jejich množství v kg/s. Ekvivalentní teplota výfukových plynů vstupující do vztahu pro výpočet využitelného zbytkového tepla se určí jako vážený průměr teplot jednotlivých plynových turbín. Blíže o modelování PPC viz také [ 1 ] . 1.2. 2.
Alternativní modely buzení Je možné modelovat alternativní struktura regulátoru buzení odpovídající ruským regulátorům tzv. silnogo dejstvija, zobrazené na následujícím obrázku. Příznak použití tohoto modelu je jeden parametr navíc na řádku typových parametrů regulátoru buzení. Předposlední čtyři parametry mají odlišný význam od zobecnělého regulátoru.
Uzmin
Uzmax
Unec
Kstat
Kp(-)
Kse(-)
TI(s)
Ts(s)
T1(s)
T3(s)
v(%/s)
Urmax
Urmin
UImin
Uimax
T1f T1u K1f K1u K_f
1
. U ST Systémový Ug abc Obvody Přídavný PI regulátor stabilizátor statiky Ig ac signál kP Sekund. U Zmax U Rmax STAB + regulátor - U nec + + + + + v K_f Σ Σ Σ Σ Σ T IQ UZ 1 + Σ U Rmin + RB US DZV U Zmin pT I Necitlivost KOR B Pevná zpětná vazba Derivační zpětná vazba K SE pT S UB K_f(K_f-1) IB Σ 1+pT S + Kanál derivace napětí K 1U pT 1U 1+pT 1U Kanál derivace frekvence K 1f pT 1f sU 1+pT 1f
UG
Obr. 2 Alternativní struktura regulátoru buzení Volbou k_f=1 bude struktura alternativního i zobecnělého regulátoru obdobná (včetně omezení integrátoru hodnotami URminURmax) s tím rozdílem, že signály DZV a STAB jsou připojeny až na výstup PI regulátoru a ne do součtového bodu (na vstup regulátoru). Následující obrázek ukazuje časový průběh veličin generátoru při blízkém třífázovém zkratu. UGEN [p.j.]
UBUD
PG [p.j.]
SG
[
1.5 1.0 0.5 0.0 0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0
-0.5 -1.0 -1.5 -2.0
t[s]
. Obr. 3 Svorkové a budící napětí, výkon a skluz generátoru s regulátorem tzv. silnogo dejstvija Zatímco rychlá tyristorová nezávislé budící soupravy po zkratu nárazově přibudí generátor, je v obrázku přibuzení pomalé –viz zelený průběh UB. Jedná se totiž o bezkroužkové buzení s točivými diodami a navíc vlivem silné pevné zpětné vazby k_f je nárůst zastaven ještě během trvání zkratu. Po vypnutí zkratu způsobí derivační zpětná vazba od svorkového napětí další odbuzování. Tyto vlivy jsou příčinou pozdějšího zotavení napětí, zmenšení elektrického výkonu po vypnutí zkratu a tím i zmenšení brzdící plochy, která je v obrázku pro názornost vybarvena oranžově. 3. Pro modelování systémových stabilizátorů PSS2A (implementovaný v budících soupravách Unitrol 5000) je od verze 2.2.11 zaveden alternativní model podle následujícího obrázku: Kanál otáček sG
Σ Kanál výkonu PG
pT S 1+pT S
Dvojitý Lead -lag
Dolní propust
pT S 1+pT S + K SP 1+pT S
+
1 + Σ [1+pT 9]5 -
U Smax K SS
1+pT S1 1+pT S2 -U Smax
U Smax 1+pT S3 1+pT S4
STAB
-U Smax
Obr. 4 Schéma alternativního modelu systémového stabilizátoru Model má oproti výchozímu modelu PSS2A zjednodušenou strukturu (vynechání dvojitých členů „washout“ na vstupech), ale zachovává hlavní vlastnosti PSS2A, kterými jsou: 1.
tlumení zesílení kanálu výkonu pro frekvence vyšší jak fZS =1/(2πTS)) způsobené členem zpoždění 1.řádu 2
odfiltrování frekvencí vyšších jak fZ9=1/(2πT9) dolní propustí, čímž se pro vyšší frekvence (odpovídající elektromechanickým kyvům 1-2 Hz) eliminuje vliv kanálu otáček 3. naopak pro nižší frekvence jak fZ9=1/(2πT9) (odpovídající systémovým kyvům 0.1 -1 Hz) eliminace vlivu kanálu výkonu. Příznakem použití alternativního modelu jsou dva přídavné parametry TS3, TS4 a T9 v úseku přídavných automatik regulátoru buzení katalogu typových parametrů bloku: 2.
KSs
KSP
TS
TS1
TS2
USmax
UOMmax
UOMzad
TiOM
TOM
KdH
KHMP
IHMPzad
Alfa
TiHMP
UHmax
TdHMP
kQ
TIQ
Qmin
Qmax
eQnec
eQom
eQmez
IBmax
IGmax
TS3 TS4 T9
Alternativní model systémového stabilizátoru byl verifikován porovnáním s měřením frekvenčních charakteristik. Tento tzv. frekvenční test spočívá v tom, že do součtového bodu regulátoru napětí se přivádí harmonický o definované amplitudě ∆U a frekvenci a sleduje se odezva veličin generátoru jako činného výkonu PG , jalového výkonu QG a svorkového napětí UG a to jednou s vypnutým a jednou ze zapnutým stabilizátorem. Následující obrázky ukazují odezvu činného výkonu a napětí (v poměrných hodnotách) generátoru na harmonický signál o amplitudě 1% při sedmi různých frekvencích od 0.25 do 1.75 Hz. u [pj] 0.987
0.5 Hz
0.25 Hz
0.986
0.75 Hz
1 Hz
0.79
p
u
0.985
p [pj] 0.8
1.25Hz
0.78
0.984
0.77
0.983
0.76
0.982
0.75
0.981
0.74
0.98
0.73
0.979
0.72
0.978
0.71
0.977
0.7
0.976
0.69
0.975
0.68
0.974
0.67
0.973
1.5Hz
0.972 0
10
20
30
40
50
60
0.66
1.75Hz 70
0.65 t80 [s]
Obr. 5 Časový průběh výkonu (modře) a napětí (fialově) při frekvenčním testu bez PSS u [pj] 0.986
p0.795 [pj]
u
0.985
0.785
0.984
0.775
0.983
0.765 p
0.982
0.755
0.981
0.745
0.98
0.735
0.979
0.725
0.978
0.715
0.977
0.705
0.976
0.695
0.975
0.685
0.974 0
10
20
30
40
50
60
70
t [s]80
0.675
Obr. 6 Časový průběh výkonu (modře) a napětí (fialově) při frekvenčním testu s PSS Časové průběhy budícího signálu i odezvy jednotlivých provozních veličin jsou ukládány do výstupního souboru a zpracovány přídavným programem HARM_ANA ( je součástí balíku programů MODES), který vypočítá zesílení zP=∆PG/∆U a zU=∆UG/∆U 3
přímo v decibelech (20*log(z)). Na následujících obrázcích jsou vypočítaná zesílení porovnány s měřením [ 2 ] - vpravo pro tlumení činného výkonu a vlevo pro tlumení kyvů napětí.
Mod 0[dB]
500
1000
1500 mHz
0
Mod 0[dB] 20
bez PSS
-5
500
1500 mHz
1000
bez PSS
15 10
-10
5
-15
0 -5
-20 -25
-15
uG
-30
s PSS
-10
s PSS
pG
-20
Obr. 7 Frekvenční charakteristiky simulačního výpočtu (tence) a měření (silně) Červeně jsou kresleny průběhy s vypnutým stabilizátorem a modře se zapnutým. Porovnání obou charakteristik ukazuje kvalitativně na schopnost tlumit kyvy dané frekvence. Charakteristika se zapnutým stabilizátorem (PSS) by měla být vždy pod charakteristikou s vypnutým PSS, aby PSS nezhoršoval tlumení. V obrázcích jsou tlustou čarou (včetně značek) kresleny i měřené charakteristiky. Je vidět vcelku zanedbatelný rozdíl mezi měřením a simulací u tlumení kyvů činného výkonu pro frekvence odpovídající systémovým kyvům (0.1-1Hz). U napětí jsou rozdíly o něco větší, ale nepřesahují 3 dB. V oblasti lokálních elektromechanických kyvů (1-2 Hz) je přesnost simulace horší a to vlivem o něco větší rezonanční frekvence pro elektromechanické kyvy. Přesto je přesnost simulačního výpočtu dostatečná pro posouzení tlumících schopností systémového stabilizátoru i pro zkoumání vlivu nastavení jednotlivých parametrů na tlumení vnucených kyvů. 1.3.
Třívstupový externí regulátor
4. V nové verzi je možno používat třívstupový externí regulátor (dosavadní model měl dva vstupy). X1
k1 1+pT 1
X2
k2 1+pT 2
X3
k3 1+pT 3
ZAD R
+ Σ
+
EPS R
+ -
Σ
NEC
R max 1 pT I R min
Y i=K RTR pro Kod=3
KP
REG R
+ Σ
+
R max R min
Y i=KOR B pro Kod=1 Y i=KOR T pro Kod=2
K1
Y1
K2
Y2
K3
Y3
Kn
Yn
Obr. 8 Blokové schéma externího stabilizátoru Upravena suma těchto proměnných tvoří regulovanou veličinu. Zadaná hodnota ZADR se zadává pro výchozí stav a lze jej měnit během výpočtu zásahem 'SETR'. Regulační odchylka EPSR se zpracuje v PI regulátoru. Výstup omezovače REGR slouží ke generací korekčních signálů Y1 až Yn. Regulátor tedy může mít až n výstupů, které se přičítají jako korekční signál k zadané hodnotě regulátoru buzení nebo regulátoru pohonu nebo hladinovému regulátoru trafa. Funkčnost nového regulátoru byla testována na případu korekce zadané hodnoty napětí pro tři paralelně pracující transformátory tak, aby jalový výkon se rozložil rovnoměrně. Následující obrázek schématicky popisuje tento případ. Do uzlu KRA4 je vyveden blok EDLS1, který reguluje pomocí externího regulátoru napětí v uzlu. Je znázorněn i sousední pilotní uzel TYN4, jehož napětí regulují bloky ECHVA. Z uzlu KRA4 jsou napájeny tři trojvinuťové transformátory vybavené hladinovými regulátory regulujícími napětí na straně 110 kV. Zadané hodnoty těchto regulátoru jsou korigovány výstupy externích regulátorů HRT401 – HRT403, jejichž parametry jsou nastaveny tak, aby zatížení jalovým výkonem všech traf bylo co nejrovnoměrnější (v poměru jmenovitých výkonů traf).
4
Obr. 9 Schéma modelu pilotního uzlu Následující obrázky ukazují časové průběhy přechodných dějů po sepnutí tří přípojnic v čase 20 s, kdy transformátory T401, T402 a T403 přejdou do paralelního provozu se stejným napětím u na sekundární straně 110 kV. V okamžiky sepnutí se aktivují externí regulátory, které změní zadané hodnoty napětí (kreslené tenkou čarou) z původní hodnoty 1.08 o hodnotu KRTR. Jelikož pouze u trafa T401 bude rozdíl zadané a skutečné hodnoty větší než necitlivost 0.009, dojde jen u tohoto trafa k přepnutí odbočky, jak ukazuje zelený tečkovaný průběh převodu trafa T401. Převody zbylých traf (kreslené tečkovaně modře a červeně) zůstanou nezměněny: [p.j. ] 1.1 T401
1.09 KRTR=0.0144
1.08
u
KRTR=-0.0035 T403
KRTR=-0.0094
T402
1.07 sepnutí spínačů Přepnutí odbočky T401
1.06 0
10
20
30
40
50
60
70
80
90
100 t [s]
Obr. 10 Časové průběhy zadaných hodnot a převodů (tečkovaně) traf a napětí na straně 110 kV Následující obrázek ukazuje časové průběhy jalových výkonů měřených na sekundární straně. Před sepnutím přípojnic tyto výkony odpovídaly odběru zátěže (mají zápornou hodnotu podle konvence, podle které je výkon vytékající z uzlu kladný). Po sepnutí se jalový výkon přerozdělí podle jmenovitého výkonu, uk a převodu jednotlivých traf. To způsobí odlehčení trafa T401 a naopak zatížení zbylých dvou. Ideální zatížení (v poměru jmenovitých výkonů) je v obrázku kresleno čerchovaně –zeleně pro trafa T401 a T402 s jmenovitým výkonem 250 MVA a modře pro T403 s výkonem 350 MVA. Korekce zadané hodnoty zapůsobí správně, protože zrovnoměrní zatížení T401 a T402 (přiblíží zatížení ideální hodnotě). Zatížení T403 pak bude takřka ideální.
5
Q 0 [MVAr] -10 T401 -20 -30 -40 -50 -60
T402
T403
-70 -80 sepnutí spínačů
Přepnutí odbočky T401
-90 0
10
20
30
40
50
60
70
80
90
100 t [s]
Obr. 11 Časové průběhy toků jalového výkonu 2. UŽIVATELSKÉ ROZHRANÍ MODMAN 1. V menu Nástroje je možno nalézt nový příkaz Tvorba scénářů. V případě, že je zadán ve stávajícím scénáři zkrat na vedení po spuštění příkazu se vygeneruje nový scénář, který obsahuje posloupnost zkratů a vypnutí všech vedení, které jsou připojeny do uzlu, v jehož blízkosti byl zkrat zadán. Všechny zkraty se spočítají na jedno spuštění programu MODES (využívá se mechanismus zapamatování stavu soustavy a návrat do tohoto stavu).
3. VERIFIKACE PROGRAMU Během vývoje a praktického používání dochází k neustálému ověřování správnosti modelů a metod implementovaných v programech MODES. K tomu slouží porovnání výsledků simulačních výpočtů: • s měřením • s výsledky v jiných programech, v [ 3 ] jsou např. porovnány průběhy zkratových proudů spočítané programy MODES a ATP. • se záznamy při systémových poruchách. Pro verifikaci může být použita (kromě již uvedeného frekvenčního testu) i skoková změna regulační odchylky. Následující obrázek ukazuje měřený (tlustě - převzato z [ 2 ] ) a simulovaný průběh svorkového napětí při o skoku regulační odchylky o 4%. UG[V] 16800 16600 16400 16200 16000 15800 15600 15400 0
1
2
3
4
5
6
7
8
9
t [s] 10
Obr. 12 Časové průběhy svorkového napětí při skokové změně regulační odchylky Simulační výpočet dobře sleduje začátek přechodného děje, dochází při něm k určitém překyvu, který však není velký. V loňském roce došlo k řadě velkým systémových poruch (viz např. [ 4 ] a [ 5 ] . Síťový simulátor tvoří důležitý prostředek pro analýzu příčin a následků takových poruch. Zároveň může taková analýza sloužit i pro verifikaci použitého dynamického modelu. Ukážeme to na příkladu dvou poruch z 25.7. a 3.8.2006.
6
Následující schéma ukazuje situaci v přenosové soustavě 25.7.2006 krátce po poledni , kdy vedení V402, V415 a V420 byla vypnuta pro přetížení a poruchu.
Chotějovice
Babylón
EPRU1, ETU
~
~
~
ETI2, PVR
EME3
~
Hradec
Výškov
Vítkov
~EPO
Bezděčín Neznášov
Čechy Střed
ECHV
~ Týnec
EDS
~ G1,2
V41 Krasíkov
V420 Přeštice
V402
Opočínek
Milín
~ EORK ~
Tábor
V203
~
Obr. 13 Schéma přenosové soustavy –stav z 25.7.2006 v 12:07 Po vypnutí vedení V402 zůstala část přenosové soustavy (v obrázku vyznačená modře) propojena se zbytkem jen jedním vedením V203. Toto vedení během necelé sekundy po výpadku vedení V402 vypnuto ochranami a modrá část přenosové soustavy přešla do ostrovního provozu, čímž se zabránilo šíření poruchy a zachoval se stabilní provoz. Vypnutí zajistila distanční ochrana měřící zdánlivou impedanci. ukazuje trajektorii zdánlivé impedance měřené v rozvodně Opočínek. 20 15 10
X/Ohm(Secondary)
5 0 -5 -10 -15 -20 -25 -30 -35 -40 0
25
50
75
100
125
150
175
200
R/Ohm(Secondary)
Obr. 14 Porovnání měření a výpočtu zdánlivé impedance vedení V203 v r. Opočínek Měřená trajektorie je kreslena tmavější barvou. Koncový bod zdánlivé impedance se v komplexní rovině RX pohybuje ve směru šipky z hodnot s velkým odporem R směrem k ose imaginárních hodnot X. V určitém čase trajektorie proniká do 1.zóny distanční ochrany (její charakteristika je dokreslena v horní části obrázku). Jelikož 1.zóna nemá závoru proti kývání, dává ochrana povel na vypnutí okamžitě, jakmile trajektorie dosáhne hranice charakteristiky. K vypnutí došlo asi 60 ms po překročení hranice 1.zóny, což představuje vlastní dobu vypínače. Tečkovaně světlejší fialovou je kreslena trajektorie spočítaná simulačním výpočtem na dynamickém modelu. Je vidět dobrá shoda měření se simulací zvláště v okolí nulové reálné složky. Relevantní pro srovnání je hustě trajektorie značkovaná čtverečky. Řídce značkovaná měřená trajektorie ukazuje kývání po vypnutí vedení. Je vidět dobrá shoda měřeného a simulovaného průběhu – síťový simulátor je schopen postihnout i činnost distančních ochran měřících zdánlivou impedanci.
7
Jako další příklad přechodné stability ukážeme chování elektrárenských bloků při zkratu, ke kterému došlo 3.8.2006. Situaci ukazuje následující obrázek:
DE
PL
400 kV 220 kV
~
~
~
~
-- vypnuto
~ V209
~
V432
~~ ETE ~
~
Čebín
V422 V476
V423
EDU
Slavětice
V433
~
Sokolnice
V417
V435 V436
APG
V437
V424
SK
V497
Obr. 15 Schéma poruchy v přenosové soustavě – 3.8.2006 v 14:51 Ke zkratu došlo na přípojnicích rozvodny 400 kV Sokolnice za probíhající dlouhodobé rekonstrukce, kdy není možné provozovat rozdílovou ochranu přípojnic. Jelikož se jednalo o přípojnicový zkrat, vypnuly 2.stupně ochran v delších časech (zkrat na vedení je vypínán 1.stupni v časech kratších než 100 ms) všechna vedení a vazební transformátor připojené do této rozvodny (v obrázku vyznačeno čárkovaně). Porucha měla následující časový průběh (za počátek časové osy budeme pro zjednodušení považovat čas 100 ms před vznikem zkratu) : 1. 2. 3. 4. 5. 6. 7.
t=0.1 s: dvoufázový zemní zkrat v r. Sokolnice t=0.32 s: dvoufázový zemní zkrat v r. Sokolnice se změní na třífázový t=0.58 s: distanční ochrany vypínají vedení V417, V435, V436 a transformátor v r. Sokolnice t=0.62 s: distanční ochrana vypíná vedení V424 t=0.65 s: distanční ochrana vypíná vedení V497 t=0.68 s: distanční ochrana vypíná vedení V423 t=0.71 s: distanční ochrany vypínají vedení V433 a V422.
Vypnutím vedení V423 v čase t=0.68 s byl zkrat eliminován (trval celkově 580 ms) a krátce na to vypnutím vedení V433 a V422 se dostaly rozvodny 400 kV Slavětice a Čebín do ostrovního provozu, jak je v naznačeno fialovou barvou. Následující obrázek ukazuje porovnání měřených a vypočítaných průběhy činného a jalového výkonu generátoru ( P a Q) vzdálenější elektrárny ETE, jejíž výkon byl do zbytku soustavy vyveden dvěma vedeními V432 a V476.
8
p.j. 1.6 P
1.4 1.2 měření
1 0.8 výpočet
0.6 NT
0.4 Q
YRV
0.2 0 0
1
2
3
4
5
6
7
t [s] 8
-0.2
Obr. 16 Porovnání měřených (tučná čára) a vypočítaných průběhů výkonů P a Q Na první části průběhu (od 0.1 do 0.68 s) jsou patrny vidět jednotlivé komutace v síti (zkraty a vypnutí vedení). Simulační výpočet vcelku věrně reprodukuje měřené průběhy (pro zjednodušení vstupních dat byl počáteční dvoufázový zemní zkrat modelován jako dvoufázový). Po vypnutí zkratu dojde k stabilnímu kývavému přechodnému ději. Stabilní přechod je umožněn i tím, že v regulačním systému turbíny byl od překročení nastaveného zrychlení rotoru turbíny aktivován tzv. elektrický urychlovač, který způsobil rychlé zavírání regulačních ventilů turbíny (v obrázku je vypočítaný průběh polohy regulačních ventilů kreslen tečkovaně jako proměnná YRV) a odpovídající snížení výkonu turbíny (v obrázku kresleno čárkovaně jako proměnná NT). V poněkud odlišné situaci byla elektrárna EDU, která byla blíže místu poruchy a která se po vypnutí vedení V433 a V422 octla v ostrově. Následující obrázky ukazují porovnání měřených a vypočítaných průběhů činného a jalového výkonu generátoru ( P a Q). p.j.
1
Tyristorový budič závislý ST_1 0.8 Q
0.6 P
výpočet
0.4
0.2 měření
0 0
1
2
3
4
5
6
7
t [s]8
Obr. 17 Měřené (tučná čára) a vypočítané průběhů výkonů P a Q pro blok se závislým budičem
9
p.j. 1
Nezávislý budič AC_1 0.8 Q 0.6 P
výpočet
0.4
0.2 měření 0 0
1
2
3
4
5
6
7
t [s]8
Obr. 18 Měřené (tučná čára) a vypočítané průběhů výkonů P a Q pro blok s nezávislým budičem Jsou kresleny průběhy na dvou blocích, z nichž každý měl jiný typ budící soupravy. Závislé buzení je realizováno řízeným tyristorovým můstkem napájeným ze svorkového napětí generátoru. Nezávislé buzení představuje pomocný alternátor (na hřídeli s hlavním generátorem) jehož napětí je usměrněno v neřízeném můstkovém usměrňovači. Více podrobností je uvedeno v [ 6 ] Porovnání měřených a vypočítaných průběhů ukazuje dobrou shodu – síťový simulátor je schopen věrně reprodukovat přechodné děje i u složitých poruch jako jsou zkraty (nesymetrické i symetrické) a přechody do ostrovního provozu. Dynamický model musí být ověřovaný porovnáním vypočítaných průběhů elektrických a mechanických veličin (např. napětí, proudů, výkonů a otáček) s dostupnými měřeními. Takto prověřený a komplexní model lze pak využít jak pro kontrolu a návrh ochran, tak i pro implementaci v dalších aplikací, např. v dispečerském trenažéru
Literatura [ 1 ] K. Máslo : Model a testování ostrovního provozu paroplynového cyklu, sborník III. Mezinárodního vědeckého symposia Elektroenergetika, St.Lesná, září 2005 [ 2 ] Ověření funkce systémového stabilizátoru TG12 v EDU, Technická zpráva ORGREZ z 8.12.2006 [ 3 ] V.Mach, K. Máslo: Model asynchronního generátoru, 7. mezinárodní konference Automatizace energetických procesů, Zlín květen 2006 [ 4 ] K. Máslo: Popis poruch v přenosové soustavě ČR, vzniklých 25. 7. 2006, Energetika č.9/2006 [ 5 ] K. Máslo: Zpráva o systémové poruše v propojení UCTE, k níž došlo 4.11.2006, Energetika č.12/2006 [ 6 ] K. Máslo: Dynamický model ES - použití v dispečerském tréninkovém simulátoru, seminář Aktuální otázky a vybrané problémy řízení ES, Poděbrady listopad 2006
10