Skip to content

Zadatak 13

DC model i linearizovani tokovi snaga na 400 kV: rešen zadatak

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

Analiza elektroenergetskih sistema, TF Čačak

Zadatak

Za dati trofazni jednopolno prikazani EES sastavljen od homogenih vodova, naponskog nivoa 400 kV, izračunati:

  1. Nepoznate θ2\theta_2 i θ3\theta_3 prethodno linearizujući odgovarajuće izraze za snage približnom zamenom sin funkcija njihovim argumentima.
  2. Tokove aktivnih snaga po granama mreže.
  3. Koliki je stepen korisnosti sistema.
  4. Nepoznate θ2\theta_2, θ3\theta_3 i P1P_1 primenom DC modela.
SP2 = (500 + j200) MVASG2 = (1000 + jQG2) MVASP3 = (400 + j200) MVASG3 = (800 + jQG3) MVA23L23 = 200 kmU2 = 420 kVU3 = 415 kVrv = 0,0295 Ω/kmxv = 0,327 Ω/kmbv = 3,5·10−6 S/kmL12 = 200 kmL13 = 200 km1SG1 = ?SP1 = (1800 + j600) MVAU1 = 400 kV

Poznato

  • U1=400U_1 = 400kV
  • U2=420U_2 = 420kV
  • U3=415U_3 = 415kV
  • L12=L13=L23=200L_{12} = L_{13} = L_{23} = 200km
  • rv=0,0295r_v = 0{,}0295Ω/km
  • xv=0,327x_v = 0{,}327Ω/km
  • bv=3,5b_v = 3{,}5·10⁻⁶ S/km
  • PP1=1800P_{P1} = 1800MW
  • QP1=600Q_{P1} = 600MVAr
  • PG2=1000P_{G2} = 1000MW
  • PP2=500P_{P2} = 500MW
  • QP2=200Q_{P2} = 200MVAr
  • PG3=800P_{G3} = 800MW
  • PP3=400P_{P3} = 400MW
  • QP3=200Q_{P3} = 200MVAr

Rešenje

Šta se traži i kojim postupkom

Čvor 1 je referentno-balansni (zadati su U1U_1 i θ1=0\theta_1 = 0), a čvorovi 2 i 3 su PV čvorovi: zadati su im moduli napona i aktivne snage injektiranja. Svi moduli napona su dakle poznati, a nepoznati su samo uglovi θ2\theta_2 i θ3\theta_3, pa su za njih dovoljne jednačine aktivnih injektiranja čvorova 2 i 3. Zadatak traži da se one linearizuju (sin⁡x≈x\sin x \approx x, cos⁡x≈1\cos x \approx 1), čime postaju sistem dve linearne jednačine. Iz dobijenih uglova računaju se tokovi snaga, gubici i stepen korisnosti, a na kraju se isti uglovi i snaga čvora 1 računaju i DC modelom.

Računa se sa trofaznim snagama u MW i MVAr, linijskim naponima u kV i admitansama u S, pa su snage injektiranja S=U⋅(YU)∗S = U \cdot (Y U)^* direktno u MVA.

Parametri voda

Z‾g=(rv+jxv)L=(5,9+j65,4) ΩY‾g=1Z‾g=(0,001368−j0,01517) SB0g2=bvL2=0,35⋅10−3 S\begin{aligned} \ph{Z}^g &= (r_v + \mathrm{j} x_v) L = (5{,}9 + \mathrm{j}65{,}4)\ \Omega \\ \ph{Y}^g &= \frac{1}{\ph{Z}^g} = (0{,}001368 - \mathrm{j}0{,}01517)\ \text{S} \\ \frac{B^{0g}}{2} &= \frac{b_v L}{2} = 0{,}35 \cdot 10^{-3}\ \text{S} \end{aligned}

Sva tri voda su ista, pa je Y‾12g=Y‾13g=Y‾23g=Y‾g\ph{Y}_{12}^g = \ph{Y}_{13}^g = \ph{Y}_{23}^g = \ph{Y}^g, sa Gg=0,001368G^g = 0{,}001368 S i Bg=−0,01517B^g = -0{,}01517 S.

Matrica admitansi

Svaki čvor je vezan sa dva voda, pa su svi dijagonalni elementi isti, kao i svi vandijagonalni:

Y‾11=Y‾22=Y‾33=2Y‾g+2⋅jB0g2=(0,002737−j0,02963) SY‾12=Y‾13=Y‾23=−Y‾g=(−0,001368+j0,01517) S\begin{gathered} \ph{Y}_{11} = \ph{Y}_{22} = \ph{Y}_{33} = 2 \ph{Y}^g + 2 \cdot \mathrm{j} \frac{B^{0g}}{2} \\ = (0{,}002737 - \mathrm{j}0{,}02963)\ \text{S} \\[1ex] \ph{Y}_{12} = \ph{Y}_{13} = \ph{Y}_{23} = -\ph{Y}^g \\ = (-0{,}001368 + \mathrm{j}0{,}01517)\ \text{S} \end{gathered}

Snage injektiranja i vrste čvorova

S‾1=(PG1−1800)+j(QG1−600)S‾2=(1000−500)+j(QG2−200)=500+j(QG2−200)S‾3=(800−400)+j(QG3−200)=400+j(QG3−200)\begin{aligned} \ph{S}_1 &= (P_{G1} - 1800) + \mathrm{j}(Q_{G1} - 600) \\ \ph{S}_2 &= (1000 - 500) + \mathrm{j}(Q_{G2} - 200) \\ &= 500 + \mathrm{j}(Q_{G2} - 200) \\ \ph{S}_3 &= (800 - 400) + \mathrm{j}(Q_{G3} - 200) \\ &= 400 + \mathrm{j}(Q_{G3} - 200) \end{aligned}
Čvor 1:U1=400U_1 = 400 kV, θ1=0\theta_1 = 0SL čvor
Čvor 2:U2=420U_2 = 420 kV, P2=500P_2 = 500 MWPV čvor
Čvor 3:U3=415U_3 = 415 kV, P3=400P_3 = 400 MWPV čvor

a) Linearizovane jednačine i uglovi

Jednačina aktivnog injektiranja čvora ii je

Pi=Ui2Gii+Ui∑j≠iUj[Gijcos⁡(θi−θj)+Bijsin⁡(θi−θj)].\begin{aligned} P_i = U_i^2 G_{ii} &+ U_i \sum_{j \neq i} U_j \Big[ G_{ij} \cos(\theta_i - \theta_j) \\ &+ B_{ij} \sin(\theta_i - \theta_j) \Big] . \end{aligned}

Sa cos⁡(θi−θj)≈1\cos(\theta_i - \theta_j) \approx 1 i sin⁡(θi−θj)≈θi−θj\sin(\theta_i - \theta_j) \approx \theta_i - \theta_j ona postaje linearna po uglovima, a sve što ne zavisi od uglova prelazi na desnu stranu:

∑j≠iUiUjBij(θi−θj)=Pi−Ui2Gii−Ui∑j≠iUjGij.\begin{gathered} \sum_{j \neq i} U_i U_j B_{ij} (\theta_i - \theta_j) \\ = P_i - U_i^2 G_{ii} - U_i \sum_{j \neq i} U_j G_{ij} . \end{gathered}

Za čvorove 2 i 3, sa θ1=0\theta_1 = 0, to su dve jednačine sa dve nepoznate. Na primer, desna strana jednačine čvora 2 je 500−4202⋅0,002737+420⋅(400+415)⋅0,001368=485,6500 - 420^2 \cdot 0{,}002737 + 420 \cdot (400 + 415) \cdot 0{,}001368 = 485{,}6 MW, a koeficijent uz θ2\theta_2 je 420⋅400⋅0,01517+420⋅415⋅0,01517=5192420 \cdot 400 \cdot 0{,}01517 + 420 \cdot 415 \cdot 0{,}01517 = 5192 MW/rad:

[5192−2644−26445161][θ2θ3]=[485,6394,3]⇒θ2=0,1792 rad=10,27∘,θ3=0,1682 rad=9,635∘\begin{gathered} \begin{bmatrix} 5192 & -2644 \\ -2644 & 5161 \end{bmatrix} \begin{bmatrix} \theta_2 \\ \theta_3 \end{bmatrix} = \begin{bmatrix} 485{,}6 \\ 394{,}3 \end{bmatrix} \\ \Rightarrow \quad \theta_2 = 0{,}1792 \text{ rad} = 10{,}27^\circ , \\ \theta_3 = 0{,}1682 \text{ rad} = 9{,}635^\circ \end{gathered}

Snaga čvora 1 i gubici

Sa dobijenim uglovima, iz tačne jednačine injektiranja čvora 1 (bez linearizacije):

P1=U12G11+U1∑j=23Uj[G1jcos⁡(−θj)+B1jsin⁡(−θj)]=−887,8 MW\begin{aligned} P_1 &= U_1^2 G_{11} + U_1 \sum_{j = 2}^{3} U_j \Big[ G_{1j} \cos(-\theta_j) \\ &\qquad + B_{1j} \sin(-\theta_j) \Big] \\ &= -887{,}8 \text{ MW} \end{aligned}

Generator u čvoru 1 daje PG1=P1+PP1=912,2P_{G1} = P_1 + P_{P1} = 912{,}2 MW. Gubici aktivne snage u mreži su zbir snaga injektiranja:

ΔPg=P1+P2+P3=−887,8+500+400=12,21 MW\begin{aligned} \Delta P_g &= P_1 + P_2 + P_3 \\ &= -887{,}8 + 500 + 400 = 12{,}21 \text{ MW} \end{aligned}

Linearizacija je približna: tačne jednačine čvorova 2 i 3 sa ovim uglovima daju P2=501,3P_2 = 501{,}3 MW i P3=401,2P_3 = 401{,}2 MW umesto zadatih 500 i 400 MW. Tačno rešenje (Newton-Raphson) je θ2=10,24∘\theta_2 = 10{,}24^\circ i θ3=9,607∘\theta_3 = 9{,}607^\circ, sa P1=−885,4P_1 = -885{,}4 MW i gubicima od 14,61 MW, pa je linearizovani model dovoljno tačan za ovakvu procenu.

c) Stepen korisnosti

Mreža prima P2+P3=900P_2 + P_3 = 900 MW iz čvorova 2 i 3, a u čvoru 1 predaje −P1=887,8-P_1 = 887{,}8 MW, pa je stepen korisnosti:

η=1−ΔPgP2+P3=1−12,21900=98,64 %\eta = 1 - \frac{\Delta P_g}{P_2 + P_3} = 1 - \frac{12{,}21}{900} = 98{,}64 \ \%

Reaktivne snage injektiranja

Qi=−Ui2Bii+Ui∑j≠iUj[Gijsin⁡(θi−θj)−Bijcos⁡(θi−θj)]\begin{aligned} Q_i = -U_i^2 B_{ii} &+ U_i \sum_{j \neq i} U_j \Big[ G_{ij} \sin(\theta_i - \theta_j) \\ &- B_{ij} \cos(\theta_i - \theta_j) \Big] \end{aligned}
Q1=−169 MVArQ2=33,13 MVArQ3=−57,33 MVAr\begin{aligned} Q_1 &= -169 \text{ MVAr} \\ Q_2 &= 33{,}13 \text{ MVAr} \\ Q_3 &= -57{,}33 \text{ MVAr} \end{aligned}

Generatori daju QGi=Qi+QPiQ_{Gi} = Q_i + Q_{Pi}: QG1=431Q_{G1} = 431 MVAr, QG2=233,1Q_{G2} = 233{,}1 MVAr i QG3=142,7Q_{G3} = 142{,}7 MVAr. Zbir Q1+Q2+Q3=−193,2Q_1 + Q_2 + Q_3 = -193{,}2 MVAr je negativan: kapacitivnosti vodova na 400 kV proizvode više reaktivne snage nego što se troši na njihovim reaktansama.

b) Tokovi snaga po granama

Tok snage od čvora ii ka čvoru jj, na početku grane ijij, je

Pij=Ui2Gg−UiUj(Ggcos⁡θij+Bgsin⁡θij)Qij=−Ui2Bg−Ui2B0g2−UiUj(Ggsin⁡θij−Bgcos⁡θij)\begin{aligned} P_{ij} &= U_i^2 G^g \\ &\quad - U_i U_j \left( G^g \cos \theta_{ij} + B^g \sin \theta_{ij} \right) \\ Q_{ij} &= -U_i^2 B^g - U_i^2 \frac{B^{0g}}{2} \\ &\quad - U_i U_j \left( G^g \sin \theta_{ij} - B^g \cos \theta_{ij} \right) \end{aligned}

sa θij=θi−θj\theta_{ij} = \theta_i - \theta_j. Na oba kraja svake grane:

GranaPijP_{ij} [MW]QijQ_{ij} [MVAr]PjiP_{ji} [MW]QjiQ_{ji} [MVAr]
1–2−461,37−95,58469,2865,49
1–3−426,42−73,47433,1431,64
2–331,97−32,35−31,91−88,97

Aktivna snaga teče iz čvorova 2 i 3 ka čvoru 1, a razlika snaga na dva kraja grane su gubici u njoj.

d) DC model

DC model zanemaruje otpornosti i otočne grane, uzima sve module napona jednake 1 r.j. i sin⁡θ≈θ\sin \theta \approx \theta, pa ostaju samo reaktanse grana. Za bazne veličine uzimamo SB=100S_B = 100 MVA i UB=400U_B = 400 kV:

ZB=UB2SB=4002100=1600 Ω,x=xvLZB=0,327⋅2001600=0,04088 r.j.\begin{gathered} Z_B = \frac{U_B^2}{S_B} = \frac{400^2}{100} = 1600\ \Omega , \\ x = \frac{x_v L}{Z_B} = \frac{0{,}327 \cdot 200}{1600} = 0{,}04088 \ \rj \end{gathered}

Matrica DC modela bez referentnog čvora i snage injektiranja u r.j. (P2=5P_2 = 5, P3=4P_3 = 4):

[−2/x1/x1/x−2/x][θ2θ3]=−[P2P3][−48,9324,4624,46−48,93][θ2θ3]=−[54]⇒θ2=0,1907 rad=10,93∘,θ3=0,1771 rad=10,15∘\begin{gathered} \begin{bmatrix} -2/x & 1/x \\ 1/x & -2/x \end{bmatrix} \begin{bmatrix} \theta_2 \\ \theta_3 \end{bmatrix} = -\begin{bmatrix} P_2 \\ P_3 \end{bmatrix} \\ \begin{bmatrix} -48{,}93 & 24{,}46 \\ 24{,}46 & -48{,}93 \end{bmatrix} \begin{bmatrix} \theta_2 \\ \theta_3 \end{bmatrix} = -\begin{bmatrix} 5 \\ 4 \end{bmatrix} \\ \Rightarrow \quad \theta_2 = 0{,}1907 \text{ rad} = 10{,}93^\circ , \\ \theta_3 = 0{,}1771 \text{ rad} = 10{,}15^\circ \end{gathered}

Gubitaka u DC modelu nema, pa čvor 1 prima svu snagu: P1=−(P2+P3)=−900P_1 = -(P_2 + P_3) = -900 MW. Tokovi su Pij=(θi−θj)/xP_{ij} = (\theta_i - \theta_j)/x:

P12=0−0,19070,04088=−4,667 r.j.=−466,7 MWP13=0−0,17710,04088=−4,333 r.j.=−433,3 MWP23=0,1907−0,17710,04088=0,3333 r.j.=33,33 MW\begin{aligned} P_{12} &= \frac{0 - 0{,}1907}{0{,}04088} = -4{,}667 \ \rj \\ &= -466{,}7 \text{ MW} \\ P_{13} &= \frac{0 - 0{,}1771}{0{,}04088} = -4{,}333 \ \rj \\ &= -433{,}3 \text{ MW} \\ P_{23} &= \frac{0{,}1907 - 0{,}1771}{0{,}04088} = 0{,}3333 \ \rj \\ &= 33{,}33 \text{ MW} \end{aligned}

DC model daje nešto veće uglove (θ2=10,93∘\theta_2 = 10{,}93^\circ prema 10,27° u linearizovanom modelu) i tokove ka čvoru 1 bez gubitaka, a raspodelu snage po granama gotovo istu kao linearizovani model.

Rezultat

a) θ2=0,1792\theta_2 = 0{,}1792 rad =10,27∘= 10{,}27^\circ, θ3=0,1682\theta_3 = 0{,}1682 rad =9,635∘= 9{,}635^\circ; P1=−887,8P_1 = -887{,}8 MW.

b) P12=−461,37P_{12} = -461{,}37 MW, P13=−426,42P_{13} = -426{,}42 MW, P23=31,97P_{23} = 31{,}97 MW (na početku grane); gubici 12,21 MW.

c) η=98,64\eta = 98{,}64 %.

d) DC model: θ2=10,93∘\theta_2 = 10{,}93^\circ, θ3=10,15∘\theta_3 = 10{,}15^\circ, P1=−900P_1 = -900 MW.

Metod u predavanju

Slični zadaci

Teme

Tokovi snaga · Gubici