Calculul puterilor activa si reactiva ale generatorului sincron cu poli inecati functionând în regim capacitiv
Se presupune puterea aparenta:
unde
Calculul tensiunii electromotoare Eq si trasarea diagramei fazoriale a tensiunilor si curentilor. Determinarea axelor d si q ale masinii
Argumentul t.e.m. Eq, unghiul , determina directia axei transversale a generatorului sincron fata de referinta sincrona care este constituita din t.e.m. a sistemului de putere infinita. Axa longitudinala a acestuia este decalata cu 90o electrice în urma, în sistemul în care axa q conduce axa d.
Se verifica valorile precedente obtinute pentru puterile din nodul 1 (se scad din aceste valori pierderile de putere din retea si trebuie sa se obtina valorile puterilor din nodul 2):
Pierderile din retea sunt:
Deci, puterile din nodul 2 s-au verificat.
În regim capacitiv de functionare a generatorului sincron scade tensiunea electromotoare Eq si creste unghiul intern al masinii , ceea ce conduce la înrautatirea conditiilor initiale, implicit înrautatirea stabilitatii.
Stabilitatea starilor de echilibru (stabilitatea statica)
in ipoteza modelarii generatorului sincron
printr-o t.e.m. in spatele reactantei sincrone in cazul
neactionarii RAT si in cazul actionarii RAT
Scopul acestei etape este trasarea grafica a puterii electrice posibila a fi produsa de generatorul cu poli inecati pentru o anumita valoare a tensiunii de excitatie in functie de unghiul rotoric la functionarea in regim inductiv si in regim capacitiv.
Un sistem este stabil intr-o anumita stare de functionare, daca dupa o mica perturbatie oarecare, ajunge intr-o stare de functionare identica, sau apropiata celei dinaintea pertubatiei.
Stabilitatea regimului permanent (a starilor de echilibru) pentru sistemul nereglat
Conditiile de legare in paralel ale unui generator la sistemul energetic (conditiile de sincronizare ale generatorului cu sistemul) sunt:
- turatia generatorului sa fie egala cu turatia sincrona;
- tensiunea la bornele generatorului sa fie egala cu tensiunea sistemului;
- succesiunea fazelor tensiunilor generatorului sa fie aceeasi cu succesiunea fazelor tensiunilor sistemului.
Dupa realizarea acestor conditii si inchiderea intreruptorului de legare la sistem a generatorului, se realizeaza o stare de regim permanent, un punct de echilibru caracterizat prin δ=0, puterea mecanica Pm transferata in Pel este zero si aceasta corespunde punctului de origine. Se presupune ca puterea mecanica creste incet, deci cresc turatia si unghiul rotoric si deci in mod corespunzator creste Pel astfel incat un punct nou de echilibru se realizeaza, punct in care Pel= Pm.
Sistemul este static stabil daca o crestere/descrestere corespunzatoare in putere mecanica cauzeaza o crestere/descrestere corespunzatoare in puterea electrica. Daca reactia sistemului se opune la aceasta, adica o crestere(descrestere) in Pm este insotita de o crestere (descrestere) a Pel, atunci nici un punct de echilibru nu poate fi atins.
Pentru o anumita valoare a Pm marcata ca “veche” sunt doua puncte de functionare 1 si 2. Daca Pm este crescuta cu o noua valoare (a) aceasta duce la un serplus de Pm in punctul 1. Acest surplus reprezinta a putere acceleratoare, va accelera rotorul astfel incat creste unghiul δ si deci Pel. Miscarea rezultata a punctului de functionare este reprezentata prin sageata catre noua stare de echilibru din punctul 5. O situatie opusa apare in punctul 2 de functionare, aici puterea acceleratoare egala cu segmentul 2-4 va accelera mai departe rotorul, va creste unghiul δ, dar aceasta crestere a unghiului δ duce la scaderea puterii electrice. Miscarea rezultatnta a punctului de functionare este aratata prin sageata. Un raspuns similar este obtinut daca se reduce Pm (b). Pentru punctele de echilibru de pe partea stanga a caracteristicii putere-unghi, miscarea rotorica este de la punctul 1 catre noul punct de echilibru 5. Pe de alta parte cand se pleaca de pe punctul de echilibru 2 pe portiunea descendenta a caracteristicii de putere, nu este posibil sa se ajunga in punctul de echilibru 6 si micsorarea rotorica continua pina in puntul de echilibru 5. Puterea mecanica posibila a fi transmisa de la generator spre sistem se noteaza cu PEq.cr.
Din cele prezentate rezulta ca generatorul, cu excitatia constanta, ce debiteaza in sistemul de putere infinita, este stabil in regim permanent numai pe portiunea ascendenta a caracteristicii de putere, adica pe portiunea pe care panta caracteristicii este pozitiva:
unde - puterea sincronizanta in regim permanent
- puterea maxima posibila a fi transmisa
Pentru o anumita valoare a excitatiei, puterea electrica posibila a fi debitata de generator reprezinta o sinusoida:
Generatorul, cu excitatia constanta, ce debiteaza in sistemul de putere infinita, este stabil in regim permanent numai pe portiunea ascendenta a caracteristicii de putere (pe portiunea pozitiva):
Valoarea PEq max se mai numeste limita stabilitatii RP si poate fi folosita pentru calculul marginii stabilitatii de RP (sau rezerva stabilitatii statice).
-Pentru regim inductiv
-Pentru regim capacitiv
Concluzie: Sistemul este stabil când punctul de functionare se afla pe partea ascendenta a caracteristicii de putere, respectiv instabil când se afla pe partea descendenta a caracteristicii de putere.
Atât timp cât punctului de functionare ii corespunde un coeficient al puterii sincronizante pozitiv, ne aflam intr-un sistem stabil de functionare. Când coeficientul este negativ, punctului de functionare ii corespunde un sistem instabil.
Stabiliatea RP pentru sistemul reglat
In aceasta parte se introduce actiunea RAT, iar influenta RAT se face in 2 stadii: a)In primul stadiu va fi dedusa ecuatia modificata a puterii electrice, datorita RAT; b)In stadiul al doilea se va arata posibilitatea de functionare dincolo de punctul critic.
Se considera un generator cu poli inecati, adica xd=xq si r=0. Se doreste obtinerea expresiei puterii electrice debitata de generator, atunci cand actioneaza RAT.
Acest lucru se realizeaza (Ug sa ramana constanta) prin modificarea excitatiei si deci a t.e.m. Eq, prin urmare va fi necesar sa se inlocuiasca Eq prin Ug si δ.
;
xTOT=xd+xext=0,6459
δ:=0,π/12,…π
δ
P1
P2
P3
P4
P5
0
0
0
0
0
0
0.2617
0.77
1.37
1.9736
2.5137
3.1737
0.5235
1.49
2.65
3.8128
4.9720
6.1311
0.7853
2.11
3.75
5.3922
7.0314
8.6705
1.0472
2.59
4.6
6.6040
8.6117
10.6194
1.3089
2.89
5.13
7.3659
9.6051
11.8444
1.5707
2.99
5.31
7.6257
9.9440
12.2683
1.8325
2.89
5.13
7.3659
9.6057
11.8444
2.0943
2.59
4.6
6.6040
8.6117
10.6194
2.3561
2.11
3.75
5.3922
7.0314
8.6707
2.6179
1.49
2.65
3.8128
4.9720
6.1311
2.8797
0.77
1.37
1.9736
2.5137
3.1737
Analiza cantitativa a dinamicii rotorului
la mici perturbatii in ipoteza utilizarii modelului clasic
pentru generator si a ecuatiilor de miscare
sub forma standard
Caracterul miscarii perturbate nu depinde nici de natura perturbatiilor, nici de valoarea lor concreta, cu conditia ca aceste sa fie suficient de mici.
Caracterul miscarii libere are scopul de a arata daca sistemul este capabil sau nu sa se intoarca la starea initiala, daca mica perturbatie a disparut. In aceasta lucrare, miscarea libera este descrisa de variatia lui in raport cu timpul.
Ecuatia de miscare rotorica la mici perturbatii in forma curenta este:
unde : Jws: [MWs2/rad];
Pm-puterea mecanica cu care este incarcat generatorul, [MW];
Pe-puterea electrica la bornele generatorului, [MW];
PD-puterea electrica de amortizare ce depinde de unghiul rotoric si de abaterea vitezei unghiularefata de viteza sincrona, [MW];
Relati de mai sus se mai poate scrie sub alta forma folosind timpul de lansare Ta care este definit:
[s]
[s]
rezulta:
Pentru un generator debitand pe barele de putere infinita, puterea de amortizare este:
rezulta:
Se presupune o mica perturbatie care are drept efect, faptul ca unghiul δ =δs+Δδ, cu δs cel din regimul de echilibru stabil. In urma liniarizarii, se obtine relatia:
(1)
unde: Pm,Pe-puterile in unitati relative, recalculate la noua putere de baza Sn
Conditiile initiale, din momentul in care cauza perturbatiilor a disparut dar sistemul a fost adus intr-o stare perturbata in care urmeaza sa oscileze liber sunt:
Se considera ca moment initial, momentul in care sistemul in cepe sa oscileze liber. Ecuatia de mai sus este o ecuatie diferentiala liniara de ordinul 2 a carei solutie e determinata de radacinile ecuatiei caracteristice:
(2)
si are radacinile
Solutia ecuatiei de miscare rotorica (1) va fi de forma:
(*)
in care constantele de integrare sunt: si se determina din conditiile initiale.
Stabilitatea sistemului depinde de valorile radacinilor p1 si p2 astfel:
- daca sunt valori reale, atunci este un mod neoscilatoriu; valoarea reala negativa corespunde unui mod amortizat, iar valoarea reala pozitiva corespunde unei instabilitati aperiodice;
- daca sunt valori complexe conjugate (), atunci duc la oscilatii de forma care reprezinta o sinusoida amortizata pentru si o sinusoida cu amplitudini crescatoare pentru .
Componenta reala a valorii proprii da amortizarea, iar componenta imaginara da frecventa de oscilatie.
O parte negativa reala reprezinta o oscilatie amortizata in timp ce o parte pozitiva reala reprezinta o oscilatie de amplitudine crescatoare. Deci, pentru o pereche complexa de radacini caracteristice, asa dupa cum se va vedea:
Frecventa de oscilatie este data prin . Aceasta relatie da frecventa curenta sau frecventa amortizata.
Raportul de amortizare (cantitatea de amortizare prezenta in raspunsul sistemului) este dat prin:
, unde - partea reala a radacinilor
- partea imaginara a radacinilor
Ecuatia de miscare rotorica la mici perturbatii sub forma standard
Ecuatia de miscare rotorica la mici perturbatii sub forma standard este:
(**)
Din relatiile anterioare rezulta:
Ecuatia caracteristica a ecuatiei de miscare rotorica (**) este:
si are radacinile
Solutia ecuatiei de miscare rotorica (**) va fi de forma:
unde: - frecventa naturala amortizata;
.
Raportul de amortizare ζ determina cantitatea de amortizare prezenta in raspunsul sistemului. Toate valorile radacinilor depinde de valorile curente ale lui psE’, D si Tm care determina tipul raspunsului. Coeficientul de inertie Tm este constanta, in timp ce D si psE’ depind de sracina generatorului. S-a constatat ca coeficientul puterii sincronizante psE’ descreste cu sarcina.
Ecuatia (*) atrata ca psE’ >0, in functie de valorile curente ale lui psE’ si D, radacinile ecuatiei carracteristice sunt fie reale, fie complexe.
Calculul parametrilor care intra in ecuatia oscilatiei libere a rotorului
Cazul
- Calculul t.e.m. tranzitorii
I=2.12-j0.3570
unde
- Calculul puterii sincronizante
- Calculul frecventei naturale
- Calculul raportului de amortizare si al frecventei naturale amortizate
pentru: D=0 (fara amortizare), ζ=0
D=1, ζ=0,0024
D=-1, ζ=-0,0024
Se creste Pm cu 0,3.
Pm(2)=2.12+0,3=2,4
Puterea mecanica trebuie sa fie egala cu cea electrica:
Noua putere reactiva corespunzatoare lui Pm=2,4 va fi:
- Calculul t.e.m. tranzitorii
- Calculul puterii electrice
(s-a verificat )
- Calculul puterii sincronizante
- Calculul frecventei naturale
- Calculul raportului de amortizare si al frecventei naturale amortizate
Pentru: D=0 (fara amortizare), ζ2=0
D=1, ζ2=0,0025
D=-1, ζ2=-0,0025
Se creste Pm cu 0,6.
Pm(2)=2.12+0,6=2,72
S3=2,72-j0,7478
Puterea mecanica trebuie sa fie egala cu cea electrica:
Noua putere reactiva corespunzatoare lui Pm=2,72 va fi:
- Calculul t.e.m. tranzitorii
- Calculul puterii electrice
(s-a verificat )
- Calculul puterii sincronizante
- Calculul frecventei naturale
- Calculul raportului de amortizare si al frecventei naturale amortizate
Pentru: D=0 (fara amortizare), ζ3=0
D=1, ζ3=0,0034
D=-1, ζ3=-0,0034
Caracteristici obtinute pentru D>0 , PsE>0, timp de studiu 10s
Caracteristici obtinute pentru D=0 , PsE>0, timp de studiu 10s
Caracteristici obtinute pentru D<0 , PsE>0, timp de studiu 10s
Studiul stabilitatii tranzitorii pentru diverse
puncte de scurtcircuit trifazat pe unul din circuitele liniei
dublu circuit, utilizand legea ariilor si
metoda Runge-Kutta de ordinul IV
Scurtcircuitul trifazat pe unul din circuitele liniei dublu circuit face ca generatorul sa oscileze fata de sistemul de putere infinita. Metoda arrilor egale poate determina stabilitatea generatorului aflat in RT fara a rezolva ecuatia de miscare. In aplicarea criteriului ariilor egale, pentru generator s-a folosit modelul clasic x’d=x’q. Puterea mecanica ramane constanta.
Nu actioneaza RAT sau RAV, t.e.m. E’ din spatele reactantei tranzitorii ramane constanta in modul. Secventa evenimentelor care a fost luata in considerare este urmatoarea:
- La timpul T=0 cind generatorul functioneaza in RP bine precizat, un scurtcircuit trifazat la pamant apare pe circuitul L2 al liniei dublucircuit. Reactanta liniei in partea stanga a defectului este λx1.
- Pentru unghiul δ=δdec respectiv la timpul t=tdec intreruptoarele la cele 2 capete ale circuitului avariat se deschid si circuitul este pierdut. Intre unghiurile δ=0 si δ=δdec, respectiv timpii t=0 si t=tdec dureaza regimul de avarie, iar pentru δ>δdec respectiv t=tdec dureaza regimul de postavarie.
- Celelalte elemente ale retelei au fost modelate ca si pina acum numai prin reactante inductive longitudinale, conform schemelor electrice de principiu si schemelor echivalente de functionare sunt prezentate in continuare.
Schema electrica de principiu
Schema electrica echivalenta pentru RPN
Pentru RPN: x’’=x’d+xT+xL/2
cu E’-t.e.m. ct in spatele reactantei tranzitorii
δ’-argumentul t.e.m. E’
U2-tensiunea in punctu 2 de putere infinita.
Sa presupunem ca se produce un scurtcircuit trifazat cu pamantul pe un ul din circuitele liniei dublu circuit, puterea debitata la bornele generatorului sincron se poate scrie in functie de reactanta xa din timpul avariei.
Conform teremei lei Thevenim, avem:
Presupunand ca se deconecteaza circuitul cu defect, noua reactanta devine:
Se observa ca puterea debitata dupa deconectarea defectului este mai mica decat in RPN (creste reactanta echivalenta a sistemului). Aceasta face ca pulsatia sa creasca, adica generatorul capata turatie suprasincrona=> Pm-Pe=Pacc<0, creste unghiul rotoric δ’. Pentru δ=δdec=δcr are loc deconectarea avariei si se face cu functionarea pe caracteristica post-avarie.
Dupa acest moment Pm-Pe=Pacc<0 si turatia devine negativa. Procesul se reia pina cand in functie de prezenta sau absenta amortizarii in modulul matematic, unghiul intern se stabilizeaza la valoarea corespunzatoare regimului post-avarie sau variaza la infinit cu aceeasi amplitudine.
Conform principiului ariilor egale, cind aria de accelerare este mai mica decait aria de frinare sistemul este stabil: Aacc≤Afranare. Daca nu se tine seama de amorizare, aceasta afirmatie este exprimata matematic prin relatia:
unde : Pmec- puterrea mecanica la bornele generatorului sincron;
Pav-puterea electrica debitata in timpul avariei;
Pe max av-maximul de putere electrica debitata in timpul avariei;
PPA= Pe max pasinδ - puterea electrica debitata dupa deconectarea liniei avariate;
Pe max pa- maximul de putere electrica debitata dupa deconectarea liniei avriate.
Locul scurtcircuitului trifazat se gaseste la barele de inalta tensiune al centralei electrice (λ=0)
In acest caz reactanta echivalenta de transfer in regim de avarie este:
In aceste conditii, puterea mecanica la arborele turbinei fiind aceeasi, generatorul se accelereaza si unghiul rotoric δ creste ceea ce conduce la pierderea stabilitatii dac nu se iau masuri de deconectare a liniei cu defect. Deconectarea liniei cu defect conduce la o noua configuratie a retelei, rezultand trecerea punctului de functionare de pe caracteristica de avrie pe cea de post-avarie. Reactanta echivalenta post-avarie este:
=0.2422
problema care se pune este cat de tarziu se poate deconecta defectul astfel incat generatorul sa nu-si piarda stabilitatea. Unghiul critic la care se face deconectarea este δcr in calculul caruia consideram ca aria de accelerare este cel mult egala cu aria de franare.
PN=Pmec = 2,2492
λ=0
pa=0 xpa=0,2422 δo=0.4454
δb=3,14-0,4454
Fig 6.1.1Pentru cazul 1 (λ=0), amortizare pozitiva (D>0), timpul de rulare 10 s, tdec=tcr.
Fig 6.1.0 Pentru cazul 1 (λ=0), amortizare D=0, timpul de rulare 10 s, tdec=tcr.
Fig 6.1.-1.Pentru cazul 1 (λ=0), amortizare D<0, timpul de rulare 10 s, tdec=tcr.
Locul scurtcircuitului trifazat se gaseste intre bornele de inalta tensiune ale centralei electrice si locul determinat de λ1(λ= λ1/2)
Punand conditia ca maximul puterii electrice debitata de generatorul sincron in timpul avariei sa fie cel mult egala cu Pm la arborele generatorului, se determina distanta pina la locul de defect notata cu λ1x. Valoarea lui λ1 se determina astfel:
λ1=0.4323
In acest caz reactanta echivalenta nu mai este infinita:
,expresia caracteristii de avarie
δcr=1,6401 tdec=0,2176
Unghiul critic de deconectare este mai mare decat in cazul 1, ceea ce inseamna ca stabilitatea se imbunatateste.
Pnmax=6,1172 u.r.
PPAmax=4,6882 u.r.
PN=2,1 u.r.
Pa=1,3006 u.r
a)
b)
c)
d)
Fig. 6.2.1 Caracteristicile putere-unghi de regim normal, de avarie, de postavarie si ariile de accelerare, respectiv franare pentru cazul 2, amortizare pozitiva (D>0)
Fig. 6.2.0 Caracteristicile putere-unghi de regim normal, de avarie, de postavarie si ariile de accelerare, respectiv franare pentru cazul 2, fara amortizare (D=0)
Fig. 6.2.-1 Caracteristicile putere-unghi pentru cazul 2, amortizare negativa (D<0)
Locul scurtcircuitului trifazat se gaseste in locul determinat de λ1(λ= λ1)
Punand conditia ca maximul puterii electrice debitata de generatorul sincron in timpul avariei sa fie cel mult egala cu Pm la arborele generatorului, se determina distanta pina la locul de defect notata cu λ1x. Valoarea lui λ1 se determina astfel:
λ1=0.4323
=0,5406
,expresia caracteristii de avarie
δcr=2,1 este mai mare decat in cazul 2
tdec=0,3581 este mai buna, creste implicit si timpul critic de deconectare.
Cu cat locul de defect este mai indepartat de barele centralei cu atat unghiul maxim necesar indepartarii defectului pentru pastrarea stabilitatii este mai mare.
Pnmax=6,1172 u.r.; PPAmax=4,6882 u.r.; PN=2,1 u.r.; Pa=2,1 u.r
Fig. 6.3.1 Caracteristicile putere-unghi de regim normal, de avarie, de postavarie si ariile de accelerare, respectiv franare pentru cazul 3, amortizare pozitiva (D>0)
Locul scurtcircuitului trifazat se gaseste intre punctul determinat de λ1 si locul determinat de λ2 ()
In aceasta situatie se observa prezenta ariei de decelerare in timpul avariei. Conditia cu care se determina acest interval se refera la egalitate ariei de accelerare si a ariei de decelerare in timpul avariei. Aceasta noua distanta este reflectat prin valoarea λ2 si se determina punand conditiile:
λ2 =0.6310
δcr=2,3 este mai mare decat in cazul 3
Cu cat locul de defect este mai indepartat de barele de inalta tensiune, scade reactanta echivalenta in timpul avariei, puterea electrica debitata in timpul avariei creste, creste si stabilitatea generatorului la mari perturbatii.
Pnmax=6,1172 u.r.; PPAmax=4,6882 u.r.; PN=2,1 u.r.; Pa=2,4431 u.r
Fig. 6.4.1 Caracteristicile putere-unghi de regim normal, de avarie, de postavarie si ariile de accelerare, respectiv franare pentru cazul 4, amortizare pozitiva (D>0)
Legea ariilor pentru λ=λ2
In cazul in care aria de accelerare cuprinsa intre caracteristica de putere mecanica si caracteristica de avarie este egala cu aria de franare, cuprinsa intre caracteristica de avarie si putere mecanica (cazul de fata), pentru a pastra stabilitatea generatorului nu mai este necesara deconectarea defectului de pe linia de transport – indepartarea defectului se va face din alte considerente cum ar fi: nivelul de tensiune; securitatea instlatiilor electrice si functionarea consumatorilor.
λ2=0.6310
=0,4466
Pm=max ,expresia caracteristii de avarie
=>δmax=2,2103
Pnmax=6,1172 u.r.; PPAmax=4,6882 u.r.; PN=2,1 u.r.; Pa=2,5110 u.r
Fig. 6.5.1 Caracteristicile putere-unghi de regim normal, de avarie, de postavarie si ariile de accelerare, respectiv franare pentru cazul 5, amortizare pozitiva (D>0)
Legea ariilor pentru λ>λ2
In acest caz avem situatia in care aria de accelerare cuprinsa sub caracteristica de putere mecanica si caracteristica de avarie este mai mica decat aria de franare, cuprinsa intre caracteristica de avarie si putere mecanica. Acest caz corespunde lui λ>λ2, =0,8155.
Nu se mai impune deconectarea deoarece generatorul este stabil, putind functiona in continuare cu defectul pe linie
λ2=0.6310
=0,4003
Pnmax=6,1172 u.r.; PPAmax=4,6882 u.r.; PN=2,1 u.r.; Pa=2,8012 u.r
a)
b)
c)
d)
Fig. 6.6.1 Caracteristicile putere-unghi de regim normal, de avarie, de postavarie si ariile de accelerare, respectiv franare pentru cazul 6, amortizare pozitiva (D>0)
In concluzie, legea ariilor ofera posibilitatea aprecierii caliative a fenomenului de stabilitate tranzitorie la schimbarea configuratiei retelei, daca generatorul este stabil cu defectul pe linie sau este instabil si atunci defectul trebuie eliminat, situatie in care se poate calcula unghiul critic de deconectare δdec=δcr.
|