Učestanost i napon su osnovne mere kvaliteta rada i isporučene električne energije. Za razliku od napona, u sistemu nema pada učestanosti u stacionarnom stanju, bez obzira na nivo napona i udaljenost prenosa. To znači da svi delovi elektroenergetskog sistema (EES-a) koji su u sinhronizmu imaju istu učestanost (u stacionarnom stanju).
Pojedini EES-i povezuju se u interkonekcije, kako bi se lakše održavala učestanost u interkonekciji.
U Evropi postoji više interkonekcija:
•UCTE – zemlje zapadne Evrope (njoj pripada i Srbija)
•NORDEL – skandinavske zemlje
•СЭВ – istočnoevropske zemlje i druge.
Za regulaciju napona postoji više regulacionih resursa. Kod regulacije učestanosti i aktivnih snaga postoji samo jedan regulacioni resurs – sinhrona generatorska grupa.
Osnovna principijelna šema sistema za regulaciju učestanosti je:
Turbinski regulatori
Primenjuju se tri vrste turbinskih regulatora:
Tahometrijski (koriste se na termoelektranama)
Akcelerotahometrijski (koriste se najčešće na hidroelektranama)
Regulatori sa elastičnom povratnom spregom, ili drugačije nazvani regulatori sa prolaznim statizmom
Zavisno od konstrukcije mernog organa, sumatora i povratne sprege postoji dva tipa turbinskih regulatora:
Mehaničko-hidraulički
Elektro-hidraulički
Sistem za regulaciju učestanosti je znatno sporiji od sistema za regulaciju napona regulacionim transformatorima i za regulaciju pobude.
Akcelerotahometrijski turbinski regulatori
Organ osetljiv na promenu brzine je kombinacija tahometar-akcelerotahometar (prvi meri brzinu, a drugi ubrzanje). Referentna brzina određena je naprezanjem opruge tahometra koja se podešava uređajem za promenu brzine.
Njihova principijelna šema je data na slici:
Hidraulički servomotori imaju uljni razvodnik kod svih tipova turbinskih regulatora. Kretanjem klipa uljnog razvodnika (odnosno njegova promena oko referentne vrednosti) dobija se superpozicijom:
Δy=Δyf+Δya+Δyx−Δyfm;
gde su pojedini članovi:
1. Član koji zavisi od učestanosti
Δyf=Δf=Δωel=Δωm=F0F−F0,
gde se indeksi el i m odnose na električnu i mehaničku ugaonu brzinu obrtanja.
2. Član koji zavisi od akcelerotahometrije (ubrzanja)
Δya=Tadtd(Δf)=TaΔf˙,
gde je Ta vremenska konstanta akcelerotahometrije (kod tahometrijskih regulatora ovaj član ne postoji).
3. Član kretanja uljnog razvodnika koji zavisi od krute povratne sprege
Δyx=r4Δx,
gde je Δx promena kretanja klipa hidrauličkog servomotora.
4. Član koji daje referentnu vrednost uređaju za promenu brzine
Δyfm=Δfm.
Za mala odstupanja oko stacionarnog stanja, može se pretpostaviti da je jednačina promene kretanja klipa hidrauličkog motora (d(Δx)/dt) proporcionalna sa kretanjem ulja uljnog razvodnika (kojim se napaja servomotor):
dtd(Δx)=−Ts′1Δy,
gde su:
Δx – promena kretanja klipa servomotora;
Δy – promena kretanja klipa uljnog razvodnika.
Promena otvaranja sprovodnog aparata turbine (Δa) direktno je proporcionalna promeni kretanja servomotora:
Ulazno-izlazni model akcelerotahometrijskog regulatora sa stalnim statizmom je:
Stalni statizam ima vrednost 4–8 %, a tipična vrednost je r=5%, kada je pojačanje turbinskog regulatora k=1/0,05=20r.j.
Vremenska konstanta je Ts≈0,1−0,5s, što znači da je:
T0≈0,050,1−0,5=2−10s.
Turbinski regulator sa prolaznim statizmom (sa elastičnom povratnom spregom)
Jednačina promene kretanja klipa uljnog razvodnika je:
Δy=Δf+Δfp−Δfm+Δfx,
gde je Δfp hod prigušivača, sveden na promenu učestanosti (kod prethodnog tipa ovde je bio član koji je zavisio od akcelerotahometrije).
Jednačina prigušivača izvodi se pretpostavljajući da je relativni pomeraj klipa prigušivača u odnosu na položaj poluge stalnog statizma srazmeran protoku ulja iz jedne u drugu komoru. Protok ulja srazmeran je relativnoj brzini pomeranja klipa u odnosu na cilindar prigušivača.
Ako se uvedu sledeće oznake:
Te – vremenska konstanta prigušivača;
Δfc – promena pomeranja klipa prigušivača;
r4′ – faktor proporcionalnosti između položaja cilindra prigušivača i klipa servomotora,
dobijaju se sledeće jednačine:
Δfp=Tedtd(Δfc−Δfp);
Δfc=r4′Δx;
Δv=r′Δa−Δfp;
r′=r3r4′,
gde je Δx položaj klipa servomotora, a r′ prolazni statizam.
Konačno se dobija matematički model turbinskog regulatora sa elastičnom povratnom spregom:
Ulazno-izlazni model turbinskog regulatora sa prolaznim statizmom je:
Ovaj tip turbinskih regulatora primenjuje se kod hidroelektrana.
Tahometrijski turbinski regulator
Primenjuje se kod termoelektrana.
Oni su specijalan slučaj bilo kojeg od dva prethodna tipa, ako se zanemari akcelerotahometrijski član (akcelerotahometrijski tip), odnosno prolazni statizam (regulatori sa prolaznim statizmom).
Zamenska šema tahometrijskog regulatora je:
Funkcija prenosa tahometrijskog regulatora je:
Wf(s)=−Wfm(s)=−r+Tss1.
Ts je konstruktivna karakteristika servomotora i ne može se menjati.
Željeno ponašanje regulatora dobija se podešavanjem:
Promena referentne brzine na ulazu u turbinski regulator (Δfm) vrši se preko uređaja za promenu brzine, kojim se kod mehaničkih regulatora upravlja posredstvom malog jednosmernog motora. Ovaj uređaj za promenu brzine je vrlo važan pri stavljanju u pogon i zaustavljanju agregata, pri sinhronizaciji, a u paralelnom sinhronom radu sa mrežom on je ulazni element pomoću koga se u sistem za regulaciju uvode naredbe za promenu odate snage agregata.
Kada se želi ograničiti otvaranje sprovodnog kola turbine, posebno u slučaju hidrauličkih turbina, dejstvo organa za detekciju promene brzine i elemenata povratne sprege moraju se blokirati i zameniti dejstvom uređaja za ograničenje otvaranja, koji direktno komanduje položajem klipa uljnog razvodnika servomotora. Pri tome, dejstva uljnog razvodnika servomotora i kola za prenos signala od servomotora do sprovodnog kola turbine zadržavaju istu funkciju, kao i u slučaju upravljanja posredstvom uređaja za promenu brzine.
Po analogiji sa prethodnim izvođenjem, za jednačinu kretanja sprovodnog aparata turbine dobija se:
dtd(Δa)=−Tsr′′Δa+Ts1Δxm,
gde su:
r′′ – statizam sistema za ograničenje otvaranja;
Δxm – promena hoda uređaja za ograničenje otvaranja.
U praksi je: r′′>>r.
Prenosna funkcija sistema za ograničenje otvaranja je:
Wxm(s)=Δxm(s)Δa(s)=r′′+Tss1.
Odziv ove funkcije je mnogo brži nego kod sistema za promenu brzine.
Pošto je r′′ vrlo veliko, onda je i ova prava vrlo strma. Teorijski, r′′→∞, odnosno k′′→0.
Regulatori sa više stepeni pojačanja
Kada se zahtevaju veće sile za pokretanje sprovodnog aparata, primenjuje se kaskadna sprega sa dva i više stepeni pojačanja, tako što se pomeraj klipa prvog hidrauličkog motora ne usmerava direktno na otvaranje turbine, već na uljni razvodnik drugog hidrauličkog servomotora (drugi stepen pojačanja):
Funkcija prenosa sistema za ograničavanje brzine sa gornje slike je:
Spoljašnja povratna sprega je uvek odlučujuća i određuje statizam regulatora.
Regulacija turbina sa dvojnim regulisanjem
Turbine sa dvojnim regulisanjem su one turbine kod kojih se snaga na vratilu reguliše pomoću dva posebna organa. U ovu grupu spadaju Kaplanove i Peltonove turbine.
Kod Kaplanovih turbina snaga se reguliše promenom otvaranja sprovodnog kola i promenom položaja, odnosno ugla zakretanja lopatica rotora.
Kod Peltonovih turbina komanduje se sekačem mlaza i promenom položaja igle mlaza.
Regulacija ovakvih turbina vrši se pomoću dva posebna regulaciona kola:
Komanda na lopatice statora servomotora.
Komanda za lopatice rotora.
Osnovna regulacija je komandovanjem lopaticama sprovodnog kola, dok je komandovanje položajem lopatica rotora prateće (pomoćno).
Ulazno-izlazni model regulatora sa dvojnim regulisanjem je dat na slici:
Modelovanje parnih turbina
Proizvodnja električne energije u termoelektranama na paru vrši se kao rezultat konverzije toplotne energije. Međutim, ova konverzija nije direktna, već se obavlja pretvaranjem toplotne energije u mehaničku, a zatim mehaničke u električnu. Konverzija toplotne energije u latentnu energiju pare vrši se u parnim kotlovima. Vodena para je medijum za prenos toplotne energije iz parnog kotla u turbinu, tako da se može reći da se turboagregat sastoji iz:
Instalacije za proizvodnju pare (parni kotao, nuklearni reaktor sa izmenjivačem toplote).
Parna turbina.
Sinhroni generator sa pomoćnim uređajima.
Para se u turbinu dovodi iz generatora pare. Za njega se može pretpostaviti da ima veliku akumulacionu sposobnost i može se smatrati neograničenim izvorom radnog fluida (pregrejane pare).
U modelovanju se može zanemariti dinamika kotla i posmatrati samo dinamika parne turbine. Ovo je moguće pošto su fizički procesi u parnom kotlu mnogo sporiji od procesa u parnoj turbini. Procesi u kotlu karakterišu se vremenskom konstantom koja je reda nekoliko minuta (vrlo spora dinamika). Parne turbine karakterišu se brzom vremenskom dinamikom, čije su vremenske konstante reda nekoliko sekundi.
Regulacija odate snage turboagregata odvija se promenom dotoka radnog fluida u turbinu, posredstvom regulacionih ventila koji kontrolišu upuštanje pare. Obično ima 4 do 8 regulacionih ventila, koji se zatvaraju i otvaraju po određenom redosledu.
Postoje dva osnovna tipa turbina:
Kondenzacione, kod kojih se iskorišćena voda (para iz turbine) vodi u kondenzator i dalje u zatvorenom ciklusu ide nazad. Ove turbine rade se kao višeoklopne, pa se na istom vratilu razlikuju turbina visokog, srednjeg i niskog pritiska.
Turbine sa međupregrevanjem pare, kod kojih se vodena para posle prolaska kroz turbinu visokog pritiska opet vodi u kotao (međupregrejač), gde se dogreva i dalje vodi u turbinu srednjeg i niskog pritiska, da bi potpuno iskorišćena stigla u kondenzator.
Kondenzacione turbine
Principijelna šema kondenzacione turbine bez međupregrevanja pare je:
Ako se pretpostavi da su pritisak i temperatura konstantni, odata snaga turbine može se izraziti kao:
Pt=Ψ(A,F),
gde je A otvaranje regulacionog ventila, a F učestanost.
Promena snage turbine u relativnim jedinicama u okolini neke radne tačke je:
Δpt=(∂a∂pt)0Δa+(∂f∂pt)0Δf=Δpta+Δptf,
gde su:
Δa – promena položaja (otvaranja) regulacionog ventila (u r.j.);
Δptf=(∂f∂pt)0Δf=eTΔf;
eT=(∂f∂pt)0 – koeficijent samoregulacije turbine (faktor osetljivosti promene snage u odnosu na promenu učestanosti pri konstantnom otvaranju).
Ukoliko su regulacioni ventili sasvim blizu mlaznika, komponenta Δpta zavisi od položaja regulacionog ventila. Međutim, ako između regulacionog ventila i ulaza u turbinu postoji parni prostor konačne zapremine, on igra ulogu akumulatora toplote, pa postoji vremenska zavisnost između promene položaja ventila i odate snage turboagregata. Ta zavisnost je:
dtd(Δpta)=TuktΔa−Tu1Δpta,
gde su:
kt – faktor proporcionalnosti između promene položaja regulacionog ventila i promene odate snage turboagregata, ili pojačanje turbine (u stacionarnom stanju za većinu turbina u okolini nominalnog stanja je kt≈1);
Tu – vremensko kašnjenje u parnom prostoru, a vrednost je 0,2−0,4s.
Prenosna funkcija kondenzacionog turboagregata bez međupregrevanja pare je:
WTc(s)=Δa(s)Δpta(s)=1+Tuskt.
Ulazno-izlazni model kondenzacionog agregata bez međupregrevanja pare je:
Kondenzacioni turboagregat sa međupregrevanjem pare
Kondenzacione turbine uglavnom se grade za manje snage, dok se za veće snage koriste turbine sa međupregrevanjem pare.
Kondenzacione turbine mogu imati turbine visokog, srednjeg i niskog pritiska.
Principijelna šema termoenergetskog postrojenja sa turbinom sa jednim međupregrevanjem pare data je na donjoj slici:
Dinamičke karakteristike pojedinih delova termoenergetskog postrojenja su:
•Dinamika turbine i parovoda visokog pritiska (praktično ista jednačina kao prethodna)
Ulazno-izlazni model kondenzacionog agregata sa međupregrevanjem pare je:
Ako je parovod između turbine srednjeg pritiska i turbine niskog pritiska vrlo kratak, onda je Tn≈0, i umesto modela trećeg reda ima se model drugog reda. Ovaj slučaj ima se kod nuklearnih termoelektrana.
Dakle, matematički model nuklearnog turboagregata je:
dtd(Δpt1)=Tukt1Δa−Tu1Δpt1;
dtd(Δpt2)=Trkt2Δpt1−Tr1Δpt2.
Δpt=Δpta+Δptf=Δpta+eTΔf;
Δpta=cvΔpt1+(1−cv)Δpt2.
Ulazno-izlazni model nuklearnog turboagregata je:
Modelovanje hidroelektrana
Njihovo modelovanje je prilično složeno.
Modeluju se preko sledećih efekata:
1. Jednačina gubitaka u dovodnom tunelu
δhg=2Δqg,
gde su:
δ – oznaka za gubitke;
Δ – oznaka za priraštaje;
hg – pad;
qg – protok kroz tunel;
g – oznaka za dovodni tunel.
2. Jednačina gubitaka pada usled kinetičke energije vode
δhw=2Δqg.
3. Jednačina protoka na ulazu u vodostan
qv=Δqg−Δqc.
4. Jednačina dovodnog tunela
Tgdtd(Δqg)+2(cg+cw)Δqg+xv+(1+cg)Δhb=0,
gde su:
Tg – vremenska konstanta tunela;
qg – protok kroz tunel;
xv – nivo vode u vodostanu;
hb – bruto pad;
cg, cw – koeficijenti gubitaka pada u dovodnom tunelu u stacionarnom stanju.
Hidroelektrane su izrazito nelinearan element. Njihovo ponašanje opisuje se kao složena funkcija od osnovnih promenljivih turbine:
•Protok vode kroz turbinu (Qt)
•Koeficijent korisnog dejstva (ηt)
•Snaga turbine (pt)
•Brzina obrtanja (N)
•Otvaranje turbine (A)
•Neto pad (Hn)
Opšta jednačina turbine je oblika:
Φ=f(Qt,ηt,pt,N,A,Hn)=0.
Analiza ponašanja hidro turbine bazira se na eksperimentalno dobijenim krivama na fizičkom modelu. Sam fizički model pravi se pre konstrukcije same turbine. Obično se te krive daju tako da se ηt, pt i A uzimaju kao parametri, pri N=Const.
To znači da se one predstavljaju kao:
Hn=f(Qt), gde su ηt, pt i A parametri, a N=Const.
Njihov izgled je dat na slici u prilogu.
Dakle, osnovne izlazne promenljive su Qt, ηt i pt, pošto se može smatrati da je u sinhronom radu brzina konstantna:
Za mala odstupanja oko stacionarnog stanja (označenog sa “0”), posle linearizacije, umesto nelinearnih jednačina dobijaju se linearizovane jednačine. Sve promenljive izražavaju se kao priraštaji. Tada je u sistemu relativnih jedinica:
pošto je promena brzine obrtanja jednaka promeni učestanosti (Δn=Δf).
U prethodnim izrazima e1−e11 su koeficijenti osetljivosti, dobijeni parcijalnim izvodima za određeno stacionarno stanje. Pri tome prava četiri koeficijenta su nezavisna:
Kada se jednačine turbine pridruže jednačinama hidroenergetskog postrojenja, dobija se model hidroelektrane.
Pošto je model izuzetno složen, ovde će se obrađivati samo pribranske hidroelektrane, koje nemaju dovodni tunel i vodostan. To znači da u opštem slučaju izloženog modela ne postoji prvih pet od osam jednačina. Ostaju samo jednačine cevovoda i neto pada.
Pribranske hidroelektrane
Kod njih ne postoji uticaj dovodnog tunela i vodostana. Ovde će se zanemariti i promena neto pada (jednačina 8). U tom slučaju za analizu ostaju samo jednačine cevovoda i promene snage turbine.
Tada je promena protoka kroz cevovod (samim tim i kroz turbinu) data jednačinom:
Ako se primeni Laplace-ova transformacija i eliminiše promena protoka kroz turbinu, za funkciju prenosa pribranske hidroelektrane pri promeni otvaranja turbine dobija se:
Ulazno izlazni model pribranske hidroelektrane je:
Može se uvesti i dalje uprošćenje ovog modela, čime se dobija uprošćeni model pribranske hidroelektrane.
On se izvodi pod pretpostavkom da se turbina posmatra kao običan otvor kroz koji ističe voda pod dejstvom zemljine teže.
U tom slučaju je:
e1=1/2;e2=1;e3=e4=0.
Onda se ostali koeficijenti izračunavaju kao:
e5=3/2;e6=1;e7=0;e8=1;e9=e10=e11=0.
Jednačina protoka kroz turbinu i snaga turbine tada su:
dtd(Δqt)=dtd(Δqc)=−Tc2Δqc+Tc2Δa.
Δpt=3Δqc−2Δa.
Pojedine prenosne funkcije su:
WaHPU(s)=Δa(s)Δpt(s)Δf=0=1+0,5Tcs1−Tcs.
WfHPU(s)=Δf(s)Δpt(s)Δa=0=0.
Može se zaključiti da promena snage agregata ne zavisi od promene učestanosti.
Model potrošača
Potrošnja u EES-u sastoji se od velikog broja individualnih potrošača. Ako se posmatra kao skoncentrisana, u stacionarnom stanju se može izraziti kao nelinearna funkcija:
Pp=f(F,V,λ(F,V))=0,
gde su:
λ – karakteriše određeni radni režim;
F – učestanost;
V – napon.
U okolini nekog stacionarnog stanja, posle linearizacije i relativiziranja promenljivih, uz uvažavanje inercije obrtnih masa potrošača, priraštaj snage potrošača je:
Δvm – dejstvo lokalnih uređaja za regulaciju napona;
Ω – mehanička brzina obrtanja potrošača.
Ako se u prethodnu jednačinu uvedu sledeći parametri:
Tp=pp01∑iJpiΩp0i2 – vremenska konstanta inercije obrtnih masa potrošača u sistemu;
epf=(∂f∂pp)0 – koeficijent elastičnosti snage potrošača pri promeni učestanosti za konstantan napon u [r.j.] (faktor samoregulacije potrošača koji zavisi od učestanosti);
epv=(∂v∂pp)0 – koeficijent elastičnosti snage potrošača pri promeni napona za konstantnu učestanost u [r.j.];
β=ΔfΔv – koeficijent osetljivosti, koji karakteriše uticaj malih promena učestanosti na male promene napona u [r.j.];
Δpλ=(∂λ∂pp)0Δλ – promena potrošnje sistema zbog promene radnih uslova (za konstantne napon i učestanost);
ep=epf+βepv – globalni koeficijent samoregulacije potrošača pri promeni učestanosti i napona u sistemu.
Tada je:
Δpp=Tpdtd(Δf)+epΔf+epvΔvm+Δpλ.
Ulazno-izlazni model potrošača je:
Efekat inercije obrtnih masa
U složenim višemašinskim EES-ima jednačina kretanja obrtnih masa nekog (i-tog) generatora može se pisati kao:
Ji∂t∂ωmi=Mti−Mgi;i=1,2,⋯,m;
Ji=Jti+Jgi+(Jei);i=1,2,⋯,m;
gde su:
ωm – mehanička brzina obrtanja;
Mt – motorni momenat (momenat turbine);
Mg – otporni momenat (momenat generatora);
Jt – momenat inercije turbine;
Jg – momenat inercije generatora;
Je – momenat inercije pobudnice (ako je obrtna).
Ako se umesto momenata uvedu snage i ako se uvede vremenska konstanta mehaničke inercije, dobija se sledeća diferencijalna jednačina:
Tgidtd(Δfi)=Δpti−Δpgi;i=1,2,⋯,m,
gde su:
Tgi=Jip0iωmiωm0i≈p0iJiωmi2;
p0i – snaga inercije u stacionarnom stanju (ili snaga turbine);
p0i=pt0i=pg0i – jednake snage turbine i generatora;
Δpti=p0ipti−pt0i;Δpgi=p0ipgi−pg0i;
Δf=F0Fi−F0=ωm0iωmi−ωm0i=Δωmi.
Sada se jednačina obrtnih masa može napisati u obliku:
Tgidtd(Δfi)=Δpti−fi(Δpp);i=1,2,⋯,m.
Ako se pretpostavi da je učestanost u sistemu zajednička:
εi=pp0pi0=p0pi0 – udeo i-tog agregata u ukupnoj snazi sistema u stacionarnom stanju;
i=1∑mεi=1,0;
ΔpL=Δpλ+epvΔvm – ekvivalentna promena potrošnje u sistemu, u [r.j.];
i=1∑mεiΔpti=Δpt.
Zamenska šema koja odgovara jednačini obrtnih masa je:
Izolovani rad elektrana
Izolovani rad je kada pojedine elektrane ili generatori napajaju sopstveno potrošačko područje. Ovaj slučaj karakteriše početni stadijum razvoja EES-a, ili poremećene režime.
Izolovani rad turboagregata
Zamenska šema u ovom slučaju je:
Prenosna funkcija WT(s) može biti:
WTc(s) – kondenzacioni turboagregat
WTr(s) – kondenzacioni turboagregat sa međupregrevanjem pare
WTN(s) – nuklearni turboagregat
Model u prostoru stanja dobija se udruživanjem jednačina stanja regulatora, turboagregata i inercije obrtnih masa:
Ovaj sistem jednačina može se napisati u formi linearnog modela u prostoru stanja kao:
dtdx=Ax+Bu+Fz,
gde su:
x=[ΔaΔvΔqcΔf]T – vektor stanja
u=[Δfm] – vektor upravljanja
z=[ΔpL] – vektor poremećaja
U prethodnom modelu podesivi parametri su: Ta, r, r′ i Te.
Najkritičniji slučaj sa stanovišta ep (faktor samoregulacije potrošača) jeste kada potrošač ima čisto omsko opterećenje, odnosno kada je ep=0.
Modelovanje elektroenergetskog sistema (EES-a)
U EES-u više elektrana radi paralelno na zajedničku mrežu, tako da se model dobija generalizacijom izvedenih modela elektrana koje rade na sopstvenu mrežu. Pri tome, uvode se sledeće pretpostavke:
Brzina obrtanja (u relativnim jedinicama) svih paralelno spregnutih generatorskih jedinica u EES-u je ista u stacionarnom stanju i u prelaznom procesu. Ova pretpostavka nameće sledeću:
Uzajamne sprege između raznih generatorskih jedinica u sistemu su krute (zanemareni ugaoni pomeraji između rotora generatora i sinhronizacione sprege).
Zanemaruju se varijacije napona u sistemu, što znači da se potpuno rasprežu P-f i Q-V potproblemi.
Potrošnja u sistemu se posmatra kao skoncentrisana, a gubici aktivne snage su pridodati potrošnji.
Pri modelovanju povezanih EES-a pretpostavlja se da su veze između njih (spojni vodovi) elastične.
Šema EES-a koji se sastoji od po jednog hidro i termo agregata je:
Na prethodnoj slici εT=p0pT0 i εH=p0pH0 su udeli turbo i hidroagregata u ukupnoj snazi sistema u stacionarnom stanju, pri čemu uvek važi εT+εH=1.
Model u prostoru stanja biće dat pod pretpostavkom da je termoelektrana predstavljena agregatom sa međupregrevanjem pare (na fosilna goriva), a hidroelektrana pribranskim agregatom. Model je:
Razlika u odnosu na napred izvedene jednačine za pojedinačne agregate je gde je potrebno dodati indeks T za termički agregat, a indeks H za hidro agregat. Takođe, razlika je i u načinu formiranja jednačine inercije obrtnih masa (jednačina 8), gde sa svojim udelima učestvuju i termički i hidro agregat.
Ako u jednom mešovitom EES-u u opštem slučaju ima mT termoelektrana, onda one unose 4mT diferencijalnih jednačina (pod pretpostavkom da su turbine kondenzacione sa međupregrevanjem pare).
Ako u EES-u u opštem slučaju ima mH hidroelektrana, onda one unose 3mH diferencijalnih jednačina (pod pretpostavkom da je regulator sa prolaznim statizmom).
Za inerciju obrtnih masa uvek postoji jedna diferencijalna jednačina.
Znači, pod napred navedenim pretpostavkama maksimalan red sistema je:
n=4mT+3mH+1.
Model u prostoru stanja je:
dtdx=Ax+Bu+Fz,
gde su:
x=[ΔaTΔpt1Δpt2Δpt3ΔaHΔvΔqcΔf]T – vektor stanja
u=[ΔfmTΔfmH] – vektor upravljanja
z=[ΔpL] – vektor poremećaja.
Modelovanje povezanih elektroenergetskih sistema
Izolovani rad EES-a danas je izuzetak. Iz ekonomskih i tehničkih razloga teži se ukrupnjavanju EES-a i stvaranju velikih interkonekcija. Istorijski su se prvo stvorili regionalni EES-i, zatim nacionalne interkonekcije i na kraju međunarodne interkonekcije, pa čak i kontinentalne interkonekcije.
U Evropi postoji više interkonekcija, od kojih su najpoznatije: UCTE, СЭВ i NORDEL.
I ove interkonekcije su povezane, ali asinhronim spojnim vodovima (DC vodovi), tako da svaka interkonekcija ima svoju sopstvenu učestanost. Interesantno je da se sad čak povezuju sistemi 50 i 60 Hz u zalivskim zemljama.
Da bi se modelovali povezani EES-i, mora se prvo izvesti model spojnog voda.
Snaga razmene između ova dva EES-a je (izraz izveden u predmetu Analiza EES-a):
V1, V2 – moduli napona na sabirnicama postrojenja odakle polaze spojni vodovi;
V1=V1/θ1; V2=V2/θ2; θ12=θ1−θ2.
Varijacija snage razmene, ili sinhronizaciona snaga, između sistema 1 i 2 je:
ΔP12=(X12′V1′V2′cosθ12)Δθ12,
gde su sa ‘ označene tranzijentne vrednosti napona i reaktansi.
Ova snaga razmene data je u stvarnim jedinicama. Međutim, da bi se izabrala bazna snaga, mora se opredeliti između sistema 1 i 2, ili zbira njihovih snaga.
Ako bi se za baznu snagu izabrala snaga sistema 1, onda je:
gde su m1T=[00⋯ms(1)]T i m2T=[00⋯ms(1)]T, pri čemu ovi vektori imaju različit broj elemenata.
Konačno, može se napisati jedinstven sistem prostora stanja, kao:
dtdx=Ax+Bu+Fzy=CTx,
gde su:
x=x1Δp12x2;u=[u1u2];z=[z1z2];
A=A1m1T0a120−α12a210−m2TA2;
B=B10000B2;F=F10000F2;
y=[Δg1Δg2]=[Δp12+bs1Δf1−α12Δp12+bs2Δf2] – regulacione greške (ulazi u sekundarne regulatore)
bs1, bs2 – faktori strmine (bias factors)
CT=[bs1C1T01−α120bs2C2T];
C1T=[00⋯Δf11]T;C2T=[00⋯Δf21]T.
Sekundarna regulacija učestanosti
Sekundarna regulacija predstavlja superponirano dejstvo na primarnu regulaciju, koje potiče od centralizovanog ili lokalnog regulatora u nekom EES-u. Pri tome, pod primarnom regulacijom podrazumeva se spontano dejstvo primarnih regulatora mašina sinhronih generatora (parne, gasne ili hidrauličke turbine) koji su osetljivi na promenu učestanosti.
Primarna regulacija je brza regulacija statičkog tipa (prelazna pojava iščezava za 10 s), sa statičkom greškom učestanosti posle pojave poremećaja. Ta stacionarna greška određena je veličinom koeficijenta globalne samoregulacije sistema etot. Njena korekcija je moguća jedino naknadnim delovanjem na turbinske regulatore. To znači da je cilj sekundarne regulacije da se eliminiše stacionarno odstupanje učestanosti ili snage razmene. Proizvodne jedinice koje učestvuju u sekundarnoj regulaciji nazivaju se regulacione jedinice.
U sekundarnoj regulaciji primenjuje se Darius-ov princip neintervencije. On kaže da se regulaciono dejstvo primenjuje samo u oblastima koje su odgovorne za promenu ravnoteže proizvodnje i potrošnje, što je izazvalo odstupanje učestanosti i snage razmene sa susedima.
Sekundarna regulacija ima tri moda: osnovni, korekcioni i havarijski.
U osnovnom regulacionom modu postoje dva aspekta:
Održavanje učestanosti i aktivnih snaga razmene po međusistemskim spojnim vodovima na željenim vrednostima, i
Kontrola raspodele aktivnih opterećenja između generatorskih jedinica u pogonu (ekonomski dispečing).
Šema regulacije u osnovnom modu je:
Osnovni cilj regulacije učestanosti i snage razmene je održavanje energetske ravnoteže između proizvodnje i potrošnje aktivnih snaga u EES-u. Realizuje se putem spore regulacije odatih generatorskih snaga, u zavisnosti od promena snaga razmene i učestanosti.
Raspodela opterećenja između generatorskih grupa (ekonomski dispečing) određuje se na bazi nekog ekonomskog kriterijuma. Taj deo spada u tercijernu regulaciju.
U havarijskom modu osnovni cilj je da se očuva integritet interkonekcije i da se pomogne sistemu sa teškoćama. Početna pomoć daje se prirodnim delovanjem regulatora. Kasnije se vrši prepodešavanje referentnog ulaza, što znači da se sistem sam ispomaže korišćenjem sopstvene obrtne rezerve.
U trećem, korektivnom modu vrši se korekcija grešaka sinhronog vremena i ugovorenih energija razmene.
Centralizovani (mrežni) regulator jedne regulacione oblasti radi prema zakonu:
ω(t)=−kfPΔf(t)−kfI0∫tΔf(τ)dτ+ω0,
gde su:
ω – izlazni signal iz mrežnog regulatora;
kfP, kfI – proporcionalno i integralno dejstvo mrežnog regulatora;
ω0 – početni odziv.
Ovaj regulator najčešće je smešten u dispečerskom centru, pa se za prenos koristi telemetrija ili neki drugi savremeni medijum.
Primer sekundarnog regulatora u sistemu jedne termo i hidroelektrane je:
Raspodela upravljačkog signala iz sekundarnog regulatora u slučaju postojanja više regulacionih agregata:
– pored učestanosti koja se mora održavati konstantnom, mora se održavati i snaga razmene. Zato je sekundarni regulator istovremeno i regulator snage razmene. To se postiže dovođenjem na ulaz učestanosti i snage razmene.
Regulator sa integralnim članom je sposoban da kompenzuje poremećaje. Integralni član je obavezan kod sekundarnog regulatora. Proporcionalno dejstvo nije obavezno. Dodaje se radi poboljšanja dinamičkih performansi sistema.
Kod praktične realizacije sistema za sekundarnu regulaciju postoji niz problema: izbor parametara, prenos signala i slično.
Obično se bira: 1/kI=60−100 s.
Sekundarna regulacija je za red veličine sporija od primarne regulacije.
Izbor ν vrši se na bazi ekonomskog dispečinga, ili na bazi heurističkih pretpostavki o željenom učešću agregata u regulaciji.