Skip to content

Zadatak 1

Tropolni kratak spoj na krajevima generatora: rešen zadatak

Zadatak prof. dr Andrije T. Sarića. Rešenje: Miloš Sarić

Analiza elektroenergetskih sistema, TF Čačak

Zadatak

Primena matrične metode pri proračunu kratkih spojeva

Na slici je prikazan prenosni sistem sa svim podacima potrebnim za proračun. Potrebno je izračunati

  1. Vremenski tok struje tropolnog kratkog spoja za kvar na krajevima generatora K1\mathrm{K}_1, i to:
  • 1° prema potpunom izrazu,
  • 2° prema praktičnoj formuli,
  • 3° prema uprošćenoj praktičnoj formuli,

b) Vremenski tok struje tropolnog kratkog spoja za kvar na mestu K2\mathrm{K}_2.

GK1YNd, T1l12 = 200 kmzv = (0,08 + j0,4) Ω/kmK22mreža

Poznato

  • SnG=100S_{nG} = 100MVA
  • UnG=15,75U_{nG} = 15{,}75kV
  • xd=1,2x_d = 1{,}2r.j.
  • xd′=0,3x_d' = 0{,}3r.j.
  • xd′′=0,2x_d'' = 0{,}2r.j.
  • xq=0,75x_q = 0{,}75r.j.
  • xq′′=0,24x_q'' = 0{,}24r.j.
  • Td′′=0,04T_d'' = 0{,}04s
  • Td′=1T_d' = 1s
  • Ta=0,15T_a = 0{,}15s
  • SnT=100S_{nT} = 100MVA
  • Un1T=15,75U_{n1}^T = 15{,}75kV
  • Un2T=231U_{n2}^T = 231kV
  • xT=0,1x_T = 0{,}1r.j.
  • rT/xT=0,05r_T / x_T = 0{,}05
  • l12=200l_{12} = 200km
  • rv=0,08r_v = 0{,}08Ω/km
  • xv=0,4x_v = 0{,}4Ω/km

Rešenje

Šta se traži i kojim postupkom

Traži se kako se u vremenu menja struja tropolnog kratkog spoja, prvo na krajevima generatora (K1\mathrm{K}_1), a zatim na kraju voda (K2\mathrm{K}_2). Tropolni kratak spoj je simetričan, pa se računa po jednoj fazi. Struja koju sinhroni generator daje u kvar ima tri dela: naizmeničnu komponentu, koja opada od suptranzijentne preko tranzijentne do ustaljene vrednosti jer reaktansa generatora raste od xd′′x_d'' preko xd′x_d' do xdx_d, jednosmernu (aperiodičnu) komponentu, koja se prigušuje sa TaT_a, i, kod mašine kojoj xd′′x_d'' i xq′′x_q'' nisu jednaki, mali član dvostruke učestanosti. Za mrežu iza sabirnice 2 podaci nisu dati, pa računamo struju koju u kvar daje generator.

Pre kvara generator je u praznom hodu, na nazivnom naponu, pa su sve tri njegove elektromotorne sile jednake naponu: E′′=E′=E=1E'' = E' = E = 1 r.j. Reaktansa xqx_q zato nije potrebna: ona bi ušla u proračun tek preko elektromotorne sile opterećenog generatora.

Bazne veličine

Uzimamo SB=100S_B = 100 MVA, a za bazne napone nazivne napone transformatora, 15,75 kV na strani generatora i 231 kV na strani voda. Podaci generatora i transformatora dati su na snazi od 100 MVA, pa ih ne treba preračunavati. Bazne struje su:

IBG=SB3 UnG=1003⋅15,75=3,666 kA,IBV=1003⋅231=0,2499 kA.\begin{aligned} I_B^{G} &= \frac{S_B}{\sqrt{3}\, U_{nG}} = \frac{100}{\sqrt{3} \cdot 15{,}75} = 3{,}666\ \text{kA}, \\ I_B^{V} &= \frac{100}{\sqrt{3} \cdot 231} = 0{,}2499\ \text{kA}. \end{aligned}

Karakteristične vrednosti naizmenične komponente za kvar K1

Efektivne vrednosti naizmenične komponente u suptranzijentnom, tranzijentnom i ustaljenom periodu su:

Ik′′=E′′xd′′=10,2=5 r.j.=18,33 kA,Ik′=E′xd′=10,3=3,333 r.j.=12,22 kA,Ik=Exd=11,2=0,8333 r.j.=3,055 kA.\begin{aligned} I_k'' &= \frac{E''}{x_d''} = \frac{1}{0{,}2} = 5\ \rj = 18{,}33\ \text{kA}, \\ I_k' &= \frac{E'}{x_d'} = \frac{1}{0{,}3} = 3{,}333\ \rj = 12{,}22\ \text{kA}, \\ I_k &= \frac{E}{x_d} = \frac{1}{1{,}2} = 0{,}8333\ \rj = 3{,}055\ \text{kA}. \end{aligned}

1° Potpuni izraz

Za generator u praznom hodu, sa uglom α\alpha koji određuje trenutak nastanka kvara (fazni stav napona faze R u tom trenutku), struja faze R je:

iR(t)=2 E[(1xd′′−1xd′)e−t/Td′′+(1xd′−1xd)e−t/Td′+1xd]sin⁡(ωt−α)+2 E2(1xd′′+1xq′′)e−t/Tasin⁡α−2 E2(1xd′′−1xq′′)e−t/Ta⋅sin⁡(2ωt−α).\begin{aligned} i_R(t) = \sqrt{2}\, E \Big[ &\Big( \frac{1}{x_d''} - \frac{1}{x_d'} \Big) e^{-t/T_d''} \\ &+ \Big( \frac{1}{x_d'} - \frac{1}{x_d} \Big) e^{-t/T_d'} + \frac{1}{x_d} \Big] \sin(\omega t - \alpha) \\ &+ \frac{\sqrt{2}\, E}{2} \Big( \frac{1}{x_d''} + \frac{1}{x_q''} \Big) e^{-t/T_a} \sin\alpha \\ &- \frac{\sqrt{2}\, E}{2} \Big( \frac{1}{x_d''} - \frac{1}{x_q''} \Big) e^{-t/T_a} \\ &\quad \cdot \sin(2\omega t - \alpha) . \end{aligned}

Prvi red je naizmenična komponenta, drugi jednosmerna, a treći član dvostruke učestanosti, koji postoji zato što rotor nije isti u d i q osi. U trenutku t=0t = 0 zbir sva tri je nula, kako i mora biti, jer struja kroz induktivnost ne može skokom da se promeni. Sa brojevima (ω=314,2\omega = 314{,}2 rad/s):

iR(t)=2[1,667 e−t/0,04+2,5 e−t/1+0,8333]sin⁡(ωt−α)+2⋅4,583 e−t/0,15sin⁡α−2⋅0,4167 e−t/0,15sin⁡(2ωt−α)  r.j.\begin{aligned} i_R(t) = \sqrt{2} \Big[ &1{,}667\, e^{-t/0{,}04} + 2{,}5\, e^{-t/1} \\ &+ 0{,}8333 \Big] \sin(\omega t - \alpha) \\ &+ \sqrt{2} \cdot 4{,}583\, e^{-t/0{,}15} \sin\alpha \\ &- \sqrt{2} \cdot 0{,}4167\, e^{-t/0{,}15} \sin(2\omega t - \alpha)\ \ \rj \end{aligned}

Jednosmerna komponenta je najveća za α=90∘\alpha = 90^\circ, kada kvar nastane u prolasku napona faze R kroz maksimum. Tada je struja posle pola periode najveća, a najveća vrednost, udarna struja, iznosi ip=13,13i_p = 13{,}13 r.j. =48,15= 48{,}15 kA, u t=9,891t = 9{,}891 ms.

2° Praktična formula

Praktična formula iz predavanja o kratkim spojevima u distributivnim mrežama zanemaruje razliku xd′′x_d'' i xq′′x_q'': jednosmernu komponentu računa sa Ik′′I_k'', a član dvostruke učestanosti izostavlja. Ukupna struja je zbir dve komponente, ik(t)=iAC(t)+iDC(t)i_k(t) = i_{AC}(t) + i_{DC}(t):

iAC(t)=2[(Ik′′−Ik′)e−t/Td′′+(Ik′−Ik)e−t/Td′+Ik]sin⁡(ωt−α),iDC(t)=2 Ik′′e−t/Tasin⁡α.\begin{aligned} i_{AC}(t) = \sqrt{2} \Big[ &\left( I_k'' - I_k' \right) e^{-t/T_d''} \\ &+ \left( I_k' - I_k \right) e^{-t/T_d'} + I_k \Big] \sin(\omega t - \alpha), \\ i_{DC}(t) = \sqrt{2}\, &I_k'' e^{-t/T_a} \sin\alpha . \end{aligned}
iAC(t)=2[1,667 e−t/0,04+2,5 e−t/1+0,8333]sin⁡(ωt−α) r.j.,iDC(t)=2⋅5 e−t/0,15sin⁡α r.j.\begin{aligned} i_{AC}(t) = \sqrt{2} \Big[ &1{,}667\, e^{-t/0{,}04} + 2{,}5\, e^{-t/1} \\ &+ 0{,}8333 \Big] \sin(\omega t - \alpha)\ \rj, \\ i_{DC}(t) = \sqrt{2} \cdot &5\, e^{-t/0{,}15} \sin\alpha\ \rj \end{aligned}

Za α=90∘\alpha = 90^\circ udarna struja je ip=13,14i_p = 13{,}14 r.j. =48,15= 48{,}15 kA: od potpunog izraza se razlikuje tek u četvrtoj cifri. Jednosmerna komponenta potpunog izraza je manja (4,583 4{,}583 umesto 5 5), ali je posle pola periode član dvostruke učestanosti dodaje skoro tačno nazad, jer je 12(1xd′′+1xq′′)+12(1xd′′−1xq′′)=1xd′′\tfrac{1}{2}\big(\tfrac{1}{x_d''} + \tfrac{1}{x_q''}\big) + \tfrac{1}{2}\big(\tfrac{1}{x_d''} - \tfrac{1}{x_q''}\big) = \tfrac{1}{x_d''}. Zato je praktična formula dovoljna za udarnu struju.

3° Uprošćena praktična formula

Uprošćena formula iz predavanja o nesimetričnim stanjima daje samo efektivnu vrednost naizmenične komponente, anvelopu struje bez jednosmerne komponente:

IAC(t)=(Ik′′−Ik′)e−t/Td′′+(Ik′−Ik)e−t/Td′+Ik=1,667 e−t/0,04+2,5 e−t/1+0,8333 r.j.\begin{aligned} I_{AC}(t) &= \left( I_k'' - I_k' \right) e^{-t/T_d''} \\ &\quad + \left( I_k' - I_k \right) e^{-t/T_d'} + I_k \\ &= 1{,}667\, e^{-t/0{,}04} + 2{,}5\, e^{-t/1} + 0{,}8333\ \rj \end{aligned}

Na primer, posle 0,1 s, kada se otvaraju kontakti brzog prekidača, IAC=3,232I_{AC} = 3{,}232 r.j. =11,85= 11{,}85 kA. Suptranzijentni deo je tada već skoro nestao (Td′′=0,04T_d'' = 0{,}04 s), a tranzijentni traje još nekoliko sekundi.

b) Kvar na mestu K2

Između generatora i mesta kvara sada su transformator i vod. Njihove impedanse na baznom naponu 231 kV, gde je ZB=UB2/SB=533,6Z_B = U_B^2/S_B = 533{,}6 Ω:

xV=xvl12ZB=0,4⋅200533,6=0,1499 r.j.,rV=0,08⋅200533,6=0,02998 r.j.,rT=0,05⋅0,1=0,005 r.j.,\begin{aligned} x_V &= \frac{x_v l_{12}}{Z_B} = \frac{0{,}4 \cdot 200}{533{,}6} = 0{,}1499\ \rj, \\ r_V &= \frac{0{,}08 \cdot 200}{533{,}6} = 0{,}02998\ \rj, \\ r_T &= 0{,}05 \cdot 0{,}1 = 0{,}005\ \rj, \end{aligned}
Xe=xT+xV=0,1+0,1499=0,2499 r.j.,Re=rT+rV=0,03498 r.j.\begin{aligned} X_e &= x_T + x_V \\ &= 0{,}1 + 0{,}1499 = 0{,}2499\ \rj, \\ R_e &= r_T + r_V = 0{,}03498\ \rj \end{aligned}

Spoljna reaktansa XeX_e dodaje se reaktansi generatora u svakom periodu (otpornosti se u naizmeničnoj komponenti zanemaruju):

Ik′′=E′′xd′′+Xe=10,4499=2,223 r.j.=0,5555 kA,Ik′=E′xd′+Xe=10,5499=1,818 r.j.=0,4545 kA,Ik=Exd+Xe=11,45=0,6897 r.j.=0,1724 kA.\begin{aligned} I_k'' &= \frac{E''}{x_d'' + X_e} = \frac{1}{0{,}4499} \\ &= 2{,}223\ \rj = 0{,}5555\ \text{kA}, \\ I_k' &= \frac{E'}{x_d' + X_e} = \frac{1}{0{,}5499} \\ &= 1{,}818\ \rj = 0{,}4545\ \text{kA}, \\ I_k &= \frac{E}{x_d + X_e} = \frac{1}{1{,}45} \\ &= 0{,}6897\ \rj = 0{,}1724\ \text{kA}. \end{aligned}

To su struje na strani od 231 kV; na strani generatora je Ik′′=8,147I_k'' = 8{,}147 kA. Menjaju se i vremenske konstante. Date Td′T_d' i Td′′T_d'' važe za kratak spoj na krajevima generatora, pa se iz njih prvo nalaze vremenske konstante praznog hoda, koje su svojstvo mašine:

Td0′=Td′xdxd′=1⋅1,20,3=4 s,Td0′′=Td′′xd′xd′′=0,04⋅0,30,2=0,06 s,\begin{aligned} T_{d0}' &= T_d' \frac{x_d}{x_d'} = 1 \cdot \frac{1{,}2}{0{,}3} = 4\ \text{s}, \\ T_{d0}'' &= T_d'' \frac{x_d'}{x_d''} = 0{,}04 \cdot \frac{0{,}3}{0{,}2} = 0{,}06\ \text{s}, \end{aligned}

a sa spoljnom reaktansom:

Td′=Td0′xd′+Xexd+Xe=4⋅0,54991,45=1,517 s,Td′′=Td0′′xd′′+Xexd′+Xe=0,06⋅0,44990,5499=0,04909 s.\begin{aligned} T_d' &= T_{d0}' \frac{x_d' + X_e}{x_d + X_e} = 4 \cdot \frac{0{,}5499}{1{,}45} \\ &= 1{,}517\ \text{s}, \\ T_d'' &= T_{d0}'' \frac{x_d'' + X_e}{x_d' + X_e} = 0{,}06 \cdot \frac{0{,}4499}{0{,}5499} \\ &= 0{,}04909\ \text{s}. \end{aligned}

Jednosmerna komponenta se prigušuje sa odnosom reaktanse i otpornosti kola. Na krajevima generatora Ta=xi/(ωra)T_a = x_i/(\omega r_a), gde je xi=(xd′′+xq′′)/2=0,22x_i = (x_d'' + x_q'')/2 = 0{,}22 r.j. inverzna reaktansa, pa je otpornost statora ra=xi/(ωTa)=0,004669r_a = x_i/(\omega T_a) = 0{,}004669 r.j. Iza transformatora i voda:

Ta=xi+Xeω(ra+Re)=0,22+0,2499314,2⋅(0,004669+0,03498)=0,03772 s.\begin{aligned} T_a &= \frac{x_i + X_e}{\omega \left( r_a + R_e \right)} \\ &= \frac{0{,}22 + 0{,}2499}{314{,}2 \cdot (0{,}004669 + 0{,}03498)} \\ &= 0{,}03772\ \text{s}. \end{aligned}

Po praktičnoj formuli, za kvar na mestu K2:

ik(t)=2[(2,223−1,818)e−t/0,04909+(1,818−0,6897)e−t/1,517+0,6897]sin⁡(ωt−α)+2⋅2,223 e−t/0,03772sin⁡α  r.j.\begin{aligned} i_k(t) = \sqrt{2} \Big[ &\left( 2{,}223 - 1{,}818 \right) e^{-t/0{,}04909} \\ &+ \left( 1{,}818 - 0{,}6897 \right) e^{-t/1{,}517} \\ &+ 0{,}6897 \Big] \sin(\omega t - \alpha) \\ &+ \sqrt{2} \cdot 2{,}223\, e^{-t/0{,}03772} \sin\alpha\ \ \rj \end{aligned}

Suptranzijentna struja je više nego dvostruko manja nego na krajevima generatora, tranzijentni deo traje duže, a jednosmerna komponenta, zbog otpornosti voda, nestaje četiri puta brže. Udarna struja (α=90∘\alpha = 90^\circ) je ip=5,448i_p = 5{,}448 r.j. =1,362= 1{,}362 kA, 1,733 1{,}733 puta temena vrednost 2 Ik′′\sqrt{2}\, I_k''.

Vremenski tok u brojevima

Efektivna vrednost naizmenične komponente IACI_{AC} za oba kvara i najveća moguća jednosmerna komponenta (α=90∘\alpha = 90^\circ) za kvar K1:

tt [s]K1: IACI_{AC} [r.j.]K1: IACI_{AC} [kA]K1: iDCi_{DC} [r.j.]K2: IACI_{AC} [r.j.]K2: IACI_{AC} [kA]
0518,337,0712,2230,556
0,014,60616,896,6152,1410,535
0,053,68913,525,0671,9280,482
0,13,23211,853,631,7990,450
0,22,89110,601,8641,6860,421
0,52,358,610,25231,5020,375
11,7536,430,0089991,2740,318
21,1724,300,000011450,99170,248
50,85023,1200,73150,183

Rezultat

K1: Ik′′=18,33I_k'' = 18{,}33 kA, Ik′=12,22I_k' = 12{,}22 kA, Ik=3,055I_k = 3{,}055 kA, udarna struja ip=48,15i_p = 48{,}15 kA.

K2: Ik′′=0,5555I_k'' = 0{,}5555 kA, Ik′=0,4545I_k' = 0{,}4545 kA, Ik=0,1724I_k = 0{,}1724 kA (na 231 kV), Td′=1,517T_d' = 1{,}517 s, Td′′=0,04909T_d'' = 0{,}04909 s, Ta=0,03772T_a = 0{,}03772 s.

Metod u predavanju

Slični zadaci

Teme

Kratki spojevi · Sinhrone mašine