Skip to content

Zadatak 4

Topološko-matrična metoda, zvezda radijalnih grana: 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

Za slučaj kvara u šemi iz prethodnog zadatka primenom topološko – matrične metode odrediti:

  1. Ekvivalentna impedansa kvara za mrežu za sva tri redosleda,
  2. Impedanse ekvivalentne zvezde radijalnih grana,
  3. Približni vremenski tok struje u granama mreže (naizmenične komponente).
SnS_nxd′′=xix_d'' = x^ixd′x_d'xdx_dx0x^0Td′′T_d''Td0′T_{d0}'TaT_a
G1, G3200 MVA0,12 r.j.0,18 r.j.1,6 r.j.0,05 r.j.0,03 s10 s0,15 s
G2500 MVA0,1 r.j.0,15 r.j.1,2 r.j.0,03 r.j.0,02 s8 s0,1 s
SnS_nxTx_TrTr_T
T1, T3200 MVA0,12 r.j.0,05⋅xT 0{,}05 \cdot x_T r.j.
T2500 MVA0,1 r.j.0,05⋅xT 0{,}05 \cdot x_T r.j.

Vodovi: lAB=150l_{AB} = 150 km, lAC=100l_{AC} = 100 km, lBC=100l_{BC} = 100 km, z‾v=(0,1+j0,4)\ph{z}_v = (0{,}1 + \mathrm{j}0{,}4) Ω/km, x0=3⋅xvx^0 = 3 \cdot x_v.

G1YNd, T1AlAB = 150 kmzv = (0,1 + j0,4) Ω/kmx0 = 3·xvBYNd, T2G2lAC = 100 kmlBC = 100 kmCYNd, T3G3

Poznato

  • SnG1=SnG3=200S_{nG1} = S_{nG3} = 200MVA
  • xd1′′=x1i=0,12x_{d1}'' = x_1^i = 0{,}12r.j.
  • xd1′=0,18x_{d1}' = 0{,}18r.j.
  • xd1=1,6x_{d1} = 1{,}6r.j.
  • x10=0,05x_1^0 = 0{,}05r.j.
  • Td1′′=0,03T_{d1}'' = 0{,}03s
  • Td01′=10T_{d01}' = 10s
  • Ta1=0,15T_{a1} = 0{,}15s
  • SnT1=SnT3=200S_{nT1} = S_{nT3} = 200MVA
  • xT1=0,12x_{T1} = 0{,}12r.j.
  • rT1/xT1=0,05r_{T1}/x_{T1} = 0{,}05
  • SnG2=500S_{nG2} = 500MVA
  • xd2′′=x2i=0,1x_{d2}'' = x_2^i = 0{,}1r.j.
  • xd2′=0,15x_{d2}' = 0{,}15r.j.
  • xd2=1,2x_{d2} = 1{,}2r.j.
  • x20=0,03x_2^0 = 0{,}03r.j.
  • Td2′′=0,02T_{d2}'' = 0{,}02s
  • Td02′=8T_{d02}' = 8s
  • Ta2=0,1T_{a2} = 0{,}1s
  • SnT2=500S_{nT2} = 500MVA
  • xT2=0,1x_{T2} = 0{,}1r.j.
  • rT2/xT2=0,05r_{T2}/x_{T2} = 0{,}05
  • lAB=150l_{AB} = 150km
  • lAC=100l_{AC} = 100km
  • lBC=100l_{BC} = 100km
  • rv=0,1r_v = 0{,}1Ω/km
  • xv=0,4x_v = 0{,}4Ω/km
  • x0/xv=3x^0/x_v = 3

Rešenje

Šta se traži i kojim postupkom

Kvar je jednopolni kratak spoj u čvoru C mreže iz zadatka 3, sa istim podacima i istim usvojenim pretpostavkama, rešen simetričnim komponentama (mreža 220 kV, prazan hod sa E=1E = 1 r.j., otpornosti zanemarene, reaktanse u r.j. na SB=100S_B = 100 MVA). Topološko-matrična metoda iz predavanja o nesimetričnim stanjima polazi od jednačine injektiranja I‾=Y‾BUSV‾\ph{\boldsymbol{I}} = \ph{\boldsymbol{Y}}_{BUS} \ph{\boldsymbol{V}}. Čvorovi se dele na karakteristične (pp) i nekarakteristične (qq, bez injektiranja), a nekarakteristični se eliminišu:

Y‾p=Y‾pp−Y‾pqY‾qq−1Y‾qp.\ph{\boldsymbol{Y}}_p = \ph{\boldsymbol{Y}}_{pp} - \ph{\boldsymbol{Y}}_{pq} \ph{\boldsymbol{Y}}_{qq}^{-1} \ph{\boldsymbol{Y}}_{qp} .

Šta je karakterističan čvor, a šta referentni, zavisi od toga šta se traži: mesto kvara za impedanse kvara, tačke iza generatora za zvezdu radijalnih grana.

a) Impedanse koje se vide sa mesta kvara

Mesto kvara C je jedini karakteristični čvor, a zajednička tačka elektromotornih sila (zemlja) referentni. Tada Y‾p\ph{\boldsymbol{Y}}_p ima jedan element, a njegova recipročna vrednost je impedansa mreže viđena iz C. Matrice admitansi su one iz zadatka 3: u direktnom redosledu generatorske grane xd′′+xTx_d'' + x_T, xd′+xTx_d' + x_T, odnosno xd+xTx_d + x_T prema zemlji, u inverznom iste kao u suptranzijentnom periodu, a u nultom samo xTx_T uzemljenih zvezdišta i vodovi sa X0=3XX^0 = 3X. Za suptranzijentni period, na primer, Yp=−j25,62Y_p = -\mathrm{j}25{,}62 r.j., pa je:

Z‾d′′=j0,03904 r.j.,Z‾d′=j0,04553 r.j.,Z‾d=j0,1801 r.j.,Z‾i=j0,03904 r.j.,Z‾0=j0,03722 r.j.\begin{aligned} \ph{Z}^{d''} &= \mathrm{j}0{,}03904\ \rj, \\ \ph{Z}^{d'} &= \mathrm{j}0{,}04553\ \rj, \\ \ph{Z}^{d} &= \mathrm{j}0{,}1801\ \rj, \\ \ph{Z}^i &= \mathrm{j}0{,}03904\ \rj, \\ \ph{Z}^0 &= \mathrm{j}0{,}03722\ \rj \end{aligned}

Ekvivalentna impedansa kvara, kojom se uticaj inverznog i nultog redosleda unosi u šemu direktnog, kod jednopolnog kratkog spoja je:

Z‾(k)=Z‾i+Z‾0=j(0,03904+0,03722)=j0,07626 r.j.\begin{aligned} \ph{Z}_{(k)} &= \ph{Z}^i + \ph{Z}^0 = \mathrm{j}(0{,}03904 + 0{,}03722) \\ &= \mathrm{j}0{,}07626\ \rj \end{aligned}

b) Ekvivalentna zvezda radijalnih grana

Svaka generatorska grana deli se na deo j(X−xd′′)\mathrm{j}(X - x_d'') uz elektromotornu silu i deo jxd′′\mathrm{j}x_d'' do mreže; tačka između njih je Mi\mathrm{M}_i. Desno od tačaka M je pasivna mreža, ista za sve periode: xd′′+xTx_d'' + x_T do čvorova A, B i C, vodovi i, iz čvora C, impedansa kvara Z‾(k)\ph{Z}_{(k)} do fiktivne tačke B iza nje, koja je referentni čvor. Karakteristični čvorovi su M1\mathrm{M}_1, M2\mathrm{M}_2 i M3\mathrm{M}_3, a A, B i C se eliminišu:

Y‾p=j[−5,222,6041,0162,604−9,562,7031,0162,703−6,54]r.j.\ph{\boldsymbol{Y}}_p = \mathrm{j} \begin{bmatrix} -5{,}22 & 2{,}604 & 1{,}016 \\ 2{,}604 & -9{,}56 & 2{,}703 \\ 1{,}016 & 2{,}703 & -6{,}54 \end{bmatrix} \rj

Admitanse radijalnih grana, od svake tačke M do tačke B, dobijaju se kao Y‾B=Y‾p[1]\ph{\boldsymbol{Y}}_B = \ph{\boldsymbol{Y}}_p [\boldsymbol{1}], gde je [1][\boldsymbol{1}] kolona jedinica, to jest kao zbirovi vrsta matrice Y‾p\ph{\boldsymbol{Y}}_p. Admitanse između tačaka M se zanemaruju: u suptranzijentnom periodu sve tačke M su na istom potencijalu EE i kroz njih struja ne teče.

Y‾B1=−j1,599,Z‾B1=j0,6253 r.j.,Y‾B2=−j4,253,Z‾B2=j0,2351 r.j.,Y‾B3=−j2,822,Z‾B3=j0,3544 r.j.\begin{aligned} \ph{Y}_{B1} &= -\mathrm{j}1{,}599, & \ph{Z}_{B1} &= \mathrm{j}0{,}6253\ \rj, \\ \ph{Y}_{B2} &= -\mathrm{j}4{,}253, & \ph{Z}_{B2} &= \mathrm{j}0{,}2351\ \rj, \\ \ph{Y}_{B3} &= -\mathrm{j}2{,}822, & \ph{Z}_{B3} &= \mathrm{j}0{,}3544\ \rj \end{aligned}

Provera: zbir admitansi zvezde je 8,674 r.j., što je tačno 1/(Xd′′+X(k))1/(X^{d''} + X_{(k)}) iz tačke a), pa radijalne grane u suptranzijentnom periodu daju struju kvara bez greške. Najkraću radijalnu granu ima G2, najjači generator, vezan na čvor B dvostrukim vodom: kraću i od G3, koji je u samom čvoru C.

c) Približni vremenski tok struja u granama

Svaka radijalna grana je sada jedan generator koji preko svoje spoljne reaktanse Xe,i=XBi−xd,i′′X_{e,i} = X_{Bi} - x_{d,i}'' napaja kvar, pa se na nju primenjuju formule za radijalnu granu. U tranzijentnom i ustaljenom periodu na red dolazi još xd′−xd′′x_d' - x_d'', odnosno xd−xd′′x_d - x_d'':

Ii′′=EXBi,Ii′=EXBi+xd′−xd′′,Ii=EXBi+xd−xd′′,Ti′′=Td0′′xd′′+Xe,ixd′+Xe,i,Ti′=Td0′xd′+Xe,ixd+Xe,i,\begin{aligned} I_i'' &= \frac{E}{X_{Bi}}, \qquad I_i' = \frac{E}{X_{Bi} + x_d' - x_d''}, \\ I_i &= \frac{E}{X_{Bi} + x_d - x_d''}, \\ T_i'' &= T_{d0}'' \frac{x_d'' + X_{e,i}}{x_d' + X_{e,i}}, \\ T_i' &= T_{d0}' \frac{x_d' + X_{e,i}}{x_d + X_{e,i}}, \end{aligned}

gde je Td0′′=Td′′xd′/xd′′T_{d0}'' = T_d'' x_d'/x_d'' (0,045 s za G1 i G3, 0,03 s za G2), a Td0′T_{d0}' dato. Naizmenična komponenta struje u grani ii (direktni redosled) je:

Ii(t)=(Ii′′−Ii′)e−t/Ti′′+(Ii′−Ii)e−t/Ti′+Ii.\begin{aligned} I_i(t) &= \left( I_i'' - I_i' \right) e^{-t/T_i''} \\ &\quad + \left( I_i' - I_i \right) e^{-t/T_i'} + I_i . \end{aligned}
GranaXe,iX_{e,i} [r.j.]Ii′′I_i'' [kA]Ii′I_i' [kA]IiI_i [kA]Ti′′T_i'' [s]Ti′T_i' [s]
G1–T10,56530,4200,4000,1920,042944,8
G2–T20,21511,1161,0710,5770,028784,309
G3–T30,29440,7400,6830,2400,041493,513

Struje u vodovima dobijaju se iz struja radijalnih grana matricom raspodele iz predavanja, računatom jednom, za pasivnu mrežu:

I‾ℓ=C‾I‾p,C‾=Y‾g(ApT−AqTY‾qq−1Y‾qp)Y‾p−1.\begin{gathered} \ph{\boldsymbol{I}}_\ell = \ph{\boldsymbol{C}} \ph{\boldsymbol{I}}_p , \\ \ph{\boldsymbol{C}} = \ph{\boldsymbol{Y}}_g \left( \boldsymbol{A}_p^{\mathrm{T}} - \boldsymbol{A}_q^{\mathrm{T}} \ph{\boldsymbol{Y}}_{qq}^{-1} \ph{\boldsymbol{Y}}_{qp} \right) \ph{\boldsymbol{Y}}_p^{-1} . \end{gathered}

Kako su sve impedanse reaktanse, elementi matrice C‾\ph{\boldsymbol{C}} su realni brojevi, udeli struje svake grane generatora u struji voda (smer od prvog slova ka drugom). Generator G3 je u čvoru C, pa njegova struja ide pravo u kvar i ni u jednom vodu nema udela:

IAB(t)=0,3333 I1(t)−0,1667 I2(t),IAC(t)=0,6667 I1(t)+0,1667 I2(t),IBC(t)=0,3333 I1(t)+0,8333 I2(t).\begin{aligned} I_{AB}(t) &= 0{,}3333\, I_1(t) - 0{,}1667\, I_2(t), \\ I_{AC}(t) &= 0{,}6667\, I_1(t) + 0{,}1667\, I_2(t), \\ I_{BC}(t) &= 0{,}3333\, I_1(t) + 0{,}8333\, I_2(t). \end{aligned}

Struja voda A–B je mala i teče od B ka A, pa je u tabeli negativna: deo struje G2 ide preko čvora A i voda A–C.

Struja kvara je zbir struja radijalnih grana, puta tri: Ik1(t)≈3[I1(t)+I2(t)+I3(t)]I_{k1}(t) \approx 3\left[ I_1(t) + I_2(t) + I_3(t) \right]. Efektivne vrednosti naizmeničnih komponenti, u kA na strani od 220 kV:

tt [s]I1I_1I2I_2I3I_3IABI_{AB}IACI_{AC}IBCI_{BC}Ik1I_{k1}
00,4201,1160,740−0,0460,4661,0706,83
0,050,4041,0730,694−0,0440,4481,0296,51
0,10,3981,0610,675−0,0440,4421,0176,40
0,20,3921,0480,659−0,0440,4361,0046,30
0,50,3801,0160,624−0,0430,4230,9746,06
10,3610,9680,573−0,0410,4020,9275,71
20,3290,8870,490−0,0380,3680,8495,12
50,2660,7310,346−0,0330,2990,6984,03
100,2180,6250,265−0,0310,2500,5943,33

Poređenje sa tačnim vrednostima iz zadatka 3 pokazuje koliko je približenje dobro. U suptranzijentnom periodu radijalne grane daju tačno Ik1′′=26,02I_{k1}'' = 26{,}02 r.j. U tranzijentnom periodu zbir je 24,62 r.j. prema tačnih 24,63 r.j. (greška 0,05 %), a u ustaljenom 11,53 r.j. prema 11,7 r.j. (1,48 %): u tim periodima tačke M više nisu na istom potencijalu, pa zanemarene grane između njih unose malu grešku, kako i predavanje kaže. Zato se metoda koristi tamo gde je potreban vremenski tok u granama, za izbor prekidača i strujnih transformatora. Inverzne i nulte komponente struja u granama raspodeljuju se po svojim šemama na isti način; kod jednopolnog kratkog spoja one su u mestu kvara jednake direktnoj, pa prate njen vremenski tok.

Rezultat

a) Z‾d′′=Z‾i=j0,03904\ph{Z}^{d''} = \ph{Z}^i = \mathrm{j}0{,}03904, Z‾0=j0,03722\ph{Z}^0 = \mathrm{j}0{,}03722, Z‾(k)=j0,07626\ph{Z}_{(k)} = \mathrm{j}0{,}07626 r.j.

b) Z‾B1=j0,6253\ph{Z}_{B1} = \mathrm{j}0{,}6253, Z‾B2=j0,2351\ph{Z}_{B2} = \mathrm{j}0{,}2351, Z‾B3=j0,3544\ph{Z}_{B3} = \mathrm{j}0{,}3544 r.j.

c) Ik1′′=6,829I_{k1}'' = 6{,}829 kA na početku kvara, Ik1≈3,026I_{k1} \approx 3{,}026 kA u ustaljenom stanju; najveća struja grane je u grani G2–T2, 1,116 kA na početku kvara.

Metod u predavanju

Slični zadaci

Teme

Kratki spojevi · Simetrične komponente