Skip to content

Predavanje 05 i 06

Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora

Predavanje prof. dr Andrije T. Sarića

Planiranje elektroenergetskih sistema, TF Čačak i FTN

Posle planiranja razvoja izvora zna se veličina, vrsta i vremenski trenutak stavljanja u pogon novih agregata. Lokacija zavisi od rasporeda energetskih resursa (hidroelektrane, termoelektrane, obnovljivi izvori energije i drugi), i da li je model jednotačkasti.

Zatim se povezuju novi izvori na mrežu, odnosno planira prenosna mreža. Neophodno je usklađivati ova dva planiranja, pošto jedno na drugo utiču.

Osnovni zahtev u planiranju prenosne mreže je da ona ne sme biti ograničavajući faktor prenosu energije. To je zato što su elementi prenosne mreže znatno jeftiniji.

Formulacija problema planiranja prenosnih mreža je slična kao kod izvora: naći strategiju razvoja prenosne mreže, tako da se minimiziraju troškovi i zadovolje sva tehnička ograničenja.

Kod planiranja prenosne mreže ne može se koristiti jednotačkasti model, već se koristi celokupni graf modela mreže.

Vrlo je važno da su planiranje izvora i prenosne mreže međusobno povezani. Sprovodi se korektivno planiranje mreže za optimalnu i više suboptimalnih planova razvoja izvora.

U planiranju razvoja prenosne mreže važni su sledeći problemi:

  • Studija dogradnje prenosne mreže.
  • Studija pouzdanosti prenosne mreže.
  • Studija dogradnje izvora reaktivne snage.
  • Provera planova razvoja prenosne mreže.

Generalno se planiraju tri tipa mreža:

  • Prenosne mreže najviših napona (iznad 220 kV).
  • Prenosne mreže visokih napona (35 do 220 kV).
  • Distributivne mreže.

Osnovni cilj planiranja prenosnih mreža jeste odrediti optimalnu konfiguraciju i parametre višenaponske mreže. Pri tome se uvažavaju porast opterećenja i razvoj izvora. Mora se zadovoljiti pouzdanost, sigurnost i ekonomičnost sistema, uz potreban kvalitet električne energije.

Planiranje razvoja izvora i prenosne mreže su međusobno uslovljeni. Kada je planiranje izvora jednotačkasto, onda je planiranje razvoja prenosne mreže korektivni faktor. Dakle, primenjuju se principi dekompozicije i agregacije.

Planiranje prenosne mreže se najčešće sprovodi iz ugla proširenja mreže. Pri tome, postoje različiti aspekti problema:

  • Obezbediti plasman energije iz novog izvora.
  • Obezbediti tranzit i povezivanje sa susednim sistemima.
  • Obezbediti dozvoljeni nivo struja kratkog spoja.
  • Bolje naponske prilike (manji padovi napona).
  • Obezbediti zahtevani nivo pouzdanosti i sigurnosti.
  • Zamena amortizovane opreme.

Zahtevi pri planiranju definišu se za:

  • Normalan pogon.
  • Poremećena stanja, odnosno (n−1)(n - 1)-koncept sigurnosti.

Troškovi vazdušnih vodova

Investicioni troškovi – troškovi izgradnje.

Zavise od naponskog nivoa, prenosne snage i prirode terena. Čine ih:

  • Troškovi projektovanja, otkupa zemljišta, montaže, nadzora …
  • Troškovi stubova, temelja i zatega, koji zavise od naponskog nivoa.
  • Troškovi provodnika, koji zavise od poprečnog preseka provodnika.

Moguće opcije realizacije voda na određenom rastojanju prikazane su na slici:

Dve opcije voda na istom rastojanju: gore pet raspona a1 sa manjim ugibom, dole tri raspona a2 sa većim ugibom.

Vidi se da je ovo klasičan optimizacioni problem.

Ukupna dužina provodnika praktično ne zavisi od raspona.

To znači da je:

  • cpc_p – cena provodnika;
  • cic_i – cena izolatora (smanjuje se sa povećanjem raspona);
  • csc_s – cena stuba (nelinearno zavisi od visine stuba);
  • cvc_v – cena voda.

Prethodne cene prikazane su na slici:

Cene c u funkciji raspona a: cena provodnika cp je stalna, cena izolatora ci blago opada, a cena stuba cs i cena voda cv, koja je zbir, prvo opadaju pa rastu. Strelica sa krive cv spušta se do raspona aE.

Troškovi izgradnje voda su:

CI′=A′+B′Vn+C′qv[NJ/km],C'_I = A' + B'V_n + C'q_v \quad [\text{NJ/km}],

gde su:

  • A′A', B′B', C′C' – konstante u usaglašenim jedinicama;
  • VnV_n – nominalni napon;
  • qvq_v – poprečni presek voda.

Struktura troškova vodova je:

  • Stubovi: 25 %.
  • Provodnici, zaštitna užad i uzemljenje stubova: 25 %.
  • Izolatori sa armaturom za vešanje: 14 %.
  • Temelji stubova: 10 %.
  • Montaža: 10 %.

Investicioni troškovi su koncentrisani na sam period izgradnje, ali se raspoređuju na čitav životni vek:

CIg=cIgCI′[NJ/km.god],C_I^g = c_I^gC'_I \quad [\text{NJ/km.god}],

gde je, na primer, cIg=3,33 %c_I^g = 3{,}33\ \% za životni vek od 30 godina.

Eksploatacioni troškovi vazdušnih vodova

Postoje dve vrste troškova:

  • Stalni troškovi pogona i održavanja.
  • Promenljivi troškovi gubitaka.

Stalni troškovi pogona i održavanja

Podrazumevaju tekuće investiciono održavanje, plate, poreze i drugo.

Obično se računaju kao neki deo (cOMog=1−3 %/godc_{OM}^{og} = 1 - 3\ \%/\text{god}) od ukupnih investicionih troškova.

Ako se ovi troškovi aktualizuju, onda je:

COMaog=cOMogCI′∑k=1n1(1+i)k=cOMog(PWF)CI′,C_{OMa}^{og} = c_{OM}^{og}C'_I\sum_{k=1}^{n}\frac{1}{(1 + i)^k} = c_{OM}^{og}(PWF)C'_I,

gde je (PWF)(PWF) faktor svođenja uniformnog niza na sadašnju ekvivalentnu vrednost.

To znači da su fiksni troškovi voda:

CF=CIg+COMaog=cIgCI′+cOMog(PWF)CI′=cFgCI′,\begin{aligned} C_F &= C_I^g + C_{OMa}^{og} = c_I^gC'_I + c_{OM}^{og}(PWF)C'_I \\ &= c_F^gC'_I, \end{aligned}

gde je cFg=cIg+cOMog(PWF)c_F^g = c_I^g + c_{OM}^{og}(PWF).

Godišnji troškovi Joule-ovih gubitaka su:

ΔPJM=3ρqvIM2[W/km];IM=PM3Vncos⁡φ[A].\Delta P_{JM} = 3\frac{\rho}{q_v}I_M^2 \quad [\text{W/km}];\quad I_M = \frac{P_M}{\sqrt{3}V_n\cos\varphi} \quad [\text{A}].

Specifični godišnji gubici energije zavise od ekvivalentnog vremena trajanja gubitaka pri maksimalnom opterećenju (τgg\tau_g^g):

ΔWJg=ΔPJMτgg;CJg=ΔWJgcwg=3ρqvIM2τggcwg,\begin{gathered} \Delta W_J^g = \Delta P_{JM}\tau_g^g; \\ C_J^g = \Delta W_J^gc_{wg} = 3\frac{\rho}{q_v}I_M^2\tau_g^gc_{wg}, \end{gathered}

gde je cwgc_{wg} prosečna cena izgubljene energije (različita od cene energije).

Aktualizovani godišnji troškovi specifičnih Joule-ovih gubitaka su:

CJag=CJg∑k=1n1(1+i)k=CJg(PWF)=3ρqvIM2τggcwg(PWF),\begin{aligned} C_{Ja}^g &= C_J^g\sum_{k=1}^{n}\frac{1}{(1 + i)^k} = C_J^g(PWF) \\ &= 3\frac{\rho}{q_v}I_M^2\tau_g^gc_{wg}(PWF), \end{aligned}

odnosno:

CJag=jΔIM2qv,\boxed{C_{Ja}^g = j_\Delta\frac{I_M^2}{q_v}},

gde je jΔ=3ρτggcwg(PWF)j_\Delta = 3\rho\tau_g^gc_{wg}(PWF).

Gubici na koronu su puno manji od Joule-ovih gubitaka. Važni su kod vazdušnih vodova vrlo visokih napona. Računaju se kao:

Ccorag=kVM(VM−Vkrcor)[NJ/km⋅god],C_{cora}^g = kV_M\left(V_M - V_{kr}^{cor}\right) \quad [\text{NJ/km}\cdot\text{god}],

gde su:

  • Vkrcor=80cvnq1v⋅10−9V_{kr}^{cor} = \dfrac{80}{c_v}\sqrt{nq_{1v}}\cdot 10^{-9} – kritični napon korone;
  • cvc_v – podužni kapacitet voda;
  • kk – faktor aktualizacije.

Dakle, ukupni eksploatacioni troškovi su:

CExag=COMaog+CJag+Ccorag0[NJ/km⋅god],C_{Exa}^g = C_{OMa}^{og} + C_{Ja}^g + \overset{0}{\cancel{C_{cora}^g}} \quad [\text{NJ/km}\cdot\text{god}],

a ukupni troškovi voda su:

Cag=CIg+CExag=CIg+COMaog+CJag=A+BVn+Cqv+jΔIM2qv.\boxed{\begin{aligned} C_a^g &= C_I^g + C_{Exa}^g = C_I^g + C_{OMa}^{og} + C_{Ja}^g \\ &= A + BV_n + Cq_v + j_\Delta\frac{I_M^2}{q_v} \end{aligned}}.

Optimizacija prenosa između dve tačke

Cilj je odrediti optimalne karakteristike prenosnog voda nominalne prenosne snage Pn=Sn=3VnInP_n = S_n = \sqrt{3}V_nI_n (faktor snage približno jednak jedinici), poznatog godišnjeg vremena iskorišćenja maksimalne snage (TMgT_M^g), koja teče između te dve tačke.

Minimizira se funkcija aktualizovanih godišnjih troškova prenosnog voda po tri promenljive: VnV_n, qvq_v i InI_n.

Koristi se Lagrange-ov metod. Proširena funkcija je:

L=Cag+λ(Pn−3VnIn)=A+BVn+Cqv+jΔIM2qv+λ(Pn−3VnIn).\begin{aligned} L &= C_a^g + \lambda\left(P_n - \sqrt{3}V_nI_n\right) \\ &= A + BV_n + Cq_v + j_\Delta\frac{I_M^2}{q_v} + \lambda\left(P_n - \sqrt{3}V_nI_n\right). \end{aligned}

Potrebni uslovi za minimum su:

∂L∂Vn=B−λ3In0=0⇒B=λ3In0.∂L∂qv=C−jΔIn02qv02=0⇒C=jΔIn02qv02.∂L∂In=2jΔIn0qv0−λ3Vn0=0⇒λ3Vn0=2jΔIn0qv0.\begin{gathered} \frac{\partial L}{\partial V_n} = B - \lambda\sqrt{3}I_{n0} = 0 \quad\Rightarrow\quad B = \lambda\sqrt{3}I_{n0}. \\[0.5em] \frac{\partial L}{\partial q_v} = C - j_\Delta\frac{I_{n0}^2}{q_{v0}^2} = 0 \quad\Rightarrow\quad C = j_\Delta\frac{I_{n0}^2}{q_{v0}^2}. \\[0.5em] \frac{\partial L}{\partial I_n} = 2j_\Delta\frac{I_{n0}}{q_{v0}} - \lambda\sqrt{3}V_{n0} = 0 \quad\Rightarrow\quad \lambda\sqrt{3}V_{n0} = 2j_\Delta\frac{I_{n0}}{q_{v0}}. \end{gathered}

Optimizacija poprečnog preseka provodnika – Kelvin-ova pravila

Polazi se od napred izvedene jednačine:

C=jΔIn02qv02,C = j_\Delta\frac{I_{n0}^2}{q_{v0}^2},

odakle je:

qv0=jΔCIn0=jΔCIM.q_{v0} = \sqrt{\frac{j_\Delta}{C}}I_{n0} = \sqrt{\frac{j_\Delta}{C}}I_M.

Za ovako izabrani poprečni presek provodnika, optimalni (minimalni) aktualizovani godišnji troškovi prenosnog voda su:

Camg=Ca0g=A+BVn0+CjΔCIn0+jΔIn02In0CjΔ;Camg=A+BVn0+CjΔIn0+CjΔIn0;Camg=A+BVn0+2CjΔIn0.\begin{gathered} C_{am}^g = C_{a0}^g = A + BV_{n0} + C\sqrt{\frac{j_\Delta}{C}}I_{n0} + j_\Delta\frac{I_{n0}^2}{I_{n0}}\sqrt{\frac{C}{j_\Delta}}; \\[0.4em] C_{am}^g = A + BV_{n0} + \sqrt{Cj_\Delta}I_{n0} + \sqrt{Cj_\Delta}I_{n0}; \\[0.4em] \boxed{C_{am}^g = A + BV_{n0} + 2\sqrt{Cj_\Delta}I_{n0}}. \end{gathered}

Iz prethodnog izraza se vidi da su pri optimumu aktualizovanih godišnjih troškova koji zavise od poprečnog preseka i gubitaka u prenosu jednaki.

Gornji zaključak je poznat kao Osnovno Kelvin-ovo pravilo:

Pri optimalnom poprečnom preseku prenosnog voda izjednačuju se godišnji troškovi srazmerni poprečnom preseku provodnika i troškovi gubitaka, aktualizovani na ceo životni vek prenosnog voda.

Izraz za optimalan poprečni presek često se daje u formi optimalne gustine struje provodnika:

J0=In0qv0=CjΔ.J_0 = \frac{I_{n0}}{q_{v0}} = \sqrt{\frac{C}{j_\Delta}}.

Vidi se da optimalna gustina struje ne zavisi od prenosne snage i naponskog nivoa voda. To je konstanta u koju su preko jΔj_\Delta uključeni podužna otpornost i godišnje vreme trajanja gubitaka.

Primer:

J0(vaz)=CjΔ≈0,6 Amm2– vazdusˇni vod.J0(kab)=CjΔ≈1,0 Amm2– kablovski vod.\begin{gathered} J_{0(vaz)} = \sqrt{\frac{C}{j_\Delta}} \approx 0{,}6\ \frac{\text{A}}{\text{mm}^2} \quad \text{– vazdušni vod.} \\[0.4em] J_{0(kab)} = \sqrt{\frac{C}{j_\Delta}} \approx 1{,}0\ \frac{\text{A}}{\text{mm}^2} \quad \text{– kablovski vod.} \end{gathered}

Problem izbora optimalnog poprečnog preseka može se postaviti i kao:

Pri kojoj struji (ili gustini struje) se za neki dati poprečni presek provodnika minimizuju aktualizovani godišnji troškovi?

CagIn=A+BVn+CqvIn+jΔInqv.\frac{C_a^g}{I_n} = \frac{A + BV_n + Cq_v}{I_n} + \frac{j_\Delta I_n}{q_v}.

Onda je:

∂(Cag/In)∂In=−A+BVn+CqvIn02+jΔqv=0,\frac{\partial\left(C_a^g/I_n\right)}{\partial I_n} = -\frac{A + BV_n + Cq_v}{I_{n0}^2} + \frac{j_\Delta}{q_v} = 0,

odakle je:

A+BVn+CqvIn02=jΔqv.\frac{A + BV_n + Cq_v}{I_{n0}^2} = \frac{j_\Delta}{q_v}.

Onda su fiksni aktualizovani troškovi voda:

CFag=A+BVn+Cqv=jΔqvIn02.\boxed{C_{Fa}^g = A + BV_n + Cq_v = \frac{j_\Delta}{q_v}I_{n0}^2}.

Ova formula definiše Modifikovano Kelvin-ovo pravilo:

Ukupni fiksni aktualizovani godišnji troškovi prenosnog voda CFagC_{Fa}^g (a ne samo deo koji je srazmeran poprečnom preseku qvq_v) treba pri optimumu da budu jednaki aktualizovanim godišnjim troškovima aktivnih gubitaka.

Ova druga gustina struje je veća od prethodne za oko 15 %, pošto zavisi od poprečnog preseka qvq_v i naponskog nivoa VnV_n.

Koji je praktičan smisao ova dva različita izraza?

Objašnjenje:

Pri određivanju optimalnog poprečnog preseka polazi se od pretpostavke da je struja maksimalna, na osnovu koje se računa optimalna gustina struje.

Kada se vod izgradi i stavi u pogon, više ne postoji ekonomski interes da se ta gustina struje pri promenljivim snagama održava, jer novi optimum daje veće gustine struje.

Optimizacija izbora prenosnog voda

Polazi se od poslednjeg napred izvedenog uslova optimalnosti:

2jΔIn0qv0=λ3Vn0=2jΔIn0jΔCIn0=2CjΔ(qv0=jΔCIn0).\begin{gathered} 2j_\Delta\frac{I_{n0}}{q_{v0}} = \lambda\sqrt{3}V_{n0} = \frac{2j_\Delta I_{n0}}{\sqrt{\dfrac{j_\Delta}{C}}I_{n0}} = 2\sqrt{Cj_\Delta} \\[0.4em] \left(q_{v0} = \sqrt{\frac{j_\Delta}{C}}I_{n0}\right). \end{gathered}

Takođe, izvedeno je i:

λ=B3In0iPn0=3Vn0In0.\lambda = \frac{B}{\sqrt{3}I_{n0}} \quad \text{i} \quad P_{n0} = \sqrt{3}V_{n0}I_{n0}.

Iz prvog izraza je:

Vn0=2CjΔλ3=2CjΔBIn0=2CjΔB⋅Pn3Vn0;Vn0=2CjΔB3⋅Pn=(CjΔ3)1/4(2PnB)1/2.\begin{gathered} V_{n0} = \frac{2\sqrt{Cj_\Delta}}{\lambda\sqrt{3}} = \frac{2\sqrt{Cj_\Delta}}{B}I_{n0} = \frac{2\sqrt{Cj_\Delta}}{B}\cdot\frac{P_n}{\sqrt{3}V_{n0}}; \\[0.4em] V_{n0} = \sqrt{\frac{2\sqrt{Cj_\Delta}}{B\sqrt{3}}}\cdot\sqrt{P_n} = \left(\frac{Cj_\Delta}{3}\right)^{1/4}\left(\frac{2P_n}{B}\right)^{1/2}. \end{gathered}

Optimalan nominalni napon prenosnog voda je srazmeran korenu nominalne prenosne snage.

Ako je poznat prenosni napon, onda je ekonomska snaga prenosa:

Pev=(32⋅BCjΔ)Vn2,\boxed{P_{ev} = \left(\frac{\sqrt{3}}{2}\cdot\frac{B}{\sqrt{Cj_\Delta}}\right)V_n^2},

odnosno ekonomska snaga prenosa je srazmerna kvadratu nominalnog napona prenosa.

Pod prethodnim uslovima, optimalni aktualizovani godišnji troškovi voda su:

Ca0g=A+BVn0+B2Vn0+B2Vn0=A+2BVn0,C_{a0}^g = A + BV_{n0} + \frac{B}{2}V_{n0} + \frac{B}{2}V_{n0} = A + 2BV_{n0},

pošto je: In0=B2CjΔVn0I_{n0} = \dfrac{B}{2\sqrt{Cj_\Delta}}V_{n0} i qv0=B2CVn0q_{v0} = \dfrac{B}{2C}V_{n0}.

Na sličan način vrši se i optimizacija prenosnih postrojenja. Samim tim, u praksi se optimizacija mora sprovesti jedinstveno za vodove i prenosna postrojenja (V + TS). Izrazi su znatno komplikovaniji, tako da se neće izvoditi (videti knjigu!).

Opteretivost prenosnih vodova

Maksimalno opterećenje prenosnih vodova određeno je sa tri faktora:

  • Dozvoljeno zagrevanje provodnika.
  • Dozvoljeni padovi napona.
  • Granica statičke stabilnosti sistema.

Termičke granice određene svojstvima materijala. Ovo se nalazi u priručnicima za standardne materijale i uslove.

Dozvoljeni padovi napona definišu granice naponske i ugaone stabilnosti (standardni postupak je kompenzacija reaktivnih snaga).

Granice statičke stabilnosti računaju se pod pretpostavkom da se napon na prijemnom kraju može održavati kompenzacijom na nominalnoj vrednosti.

Posmatra se sistem sa zanemarenim gubicima aktivne snage:

Zamenska π-šema voda između napona V1 i V2: redna reaktansa Xv i otočne susceptanse Bv/2 na oba kraja.

Za prethodnu šemu je:

Pp1=Pp2=Pp=V1V2Xvsin⁡θ12;Qp1=V1V2Xvcos⁡θ12−(1+BvXv2)V12Xv=V1Xv(V2cos⁡θ12−AV1);Qp2=(1+BvXv2)V22Xv−V1V2Xvcos⁡θ12=V2Xv(AV2−V1cos⁡θ12),\begin{aligned} P_{p1} &= P_{p2} = P_p = \frac{V_1V_2}{X_v}\sin\theta_{12}; \\ Q_{p1} &= \frac{V_1V_2}{X_v}\cos\theta_{12} - \left(1 + \frac{B_vX_v}{2}\right)\frac{V_1^2}{X_v} \\ &= \frac{V_1}{X_v}\left(V_2\cos\theta_{12} - AV_1\right); \\ Q_{p2} &= \left(1 + \frac{B_vX_v}{2}\right)\frac{V_2^2}{X_v} - \frac{V_1V_2}{X_v}\cos\theta_{12} \\ &= \frac{V_2}{X_v}\left(AV_2 - V_1\cos\theta_{12}\right), \end{aligned}

gde je: A=1+BvXv2A = 1 + \dfrac{B_vX_v}{2}

Ako se definiše ugao faktora snage na prijemnom kraju voda kao φ2=arctg⁡(Qp2Pp)\varphi_2 = \operatorname{arctg}\left(\dfrac{Q_{p2}}{P_p}\right) i ako se u ovaj izraz zamene PpP_p i Qp2Q_{p2}, dobija se:

tg⁡φ2=V2Xv(AV2−V1cos⁡θ12)V1V2Xvsin⁡θ12=AV2−V1cos⁡θ12V1sin⁡θ12,\operatorname{tg}\varphi_2 = \frac{\dfrac{V_2}{X_v}\left(AV_2 - V_1\cos\theta_{12}\right)}{\dfrac{V_1V_2}{X_v}\sin\theta_{12}} = \frac{AV_2 - V_1\cos\theta_{12}}{V_1\sin\theta_{12}},

odakle je:

V2=V1A(cos⁡θ12+sin⁡θ12tg⁡φ2).V_2 = \frac{V_1}{A}\left(\cos\theta_{12} + \sin\theta_{12}\operatorname{tg}\varphi_2\right).

Onda je aktivno opterećenje prenosnog voda:

Pp=V1V2Xvsin⁡θ12=V12AXvsin⁡θ12(cos⁡θ12+sin⁡θ12tg⁡φ2).\boxed{\begin{aligned} P_p &= \frac{V_1V_2}{X_v}\sin\theta_{12} \\ &= \frac{V_1^2}{AX_v}\sin\theta_{12}\left(\cos\theta_{12} + \sin\theta_{12}\operatorname{tg}\varphi_2\right) \end{aligned}}.

Na primer, ako se napon na početku voda održava na V1=1,05V_1 = 1{,}05 r.j., a dozvoli pad napona od 10 %, odnosno V2=0,95V_2 = 0{,}95 r.j., onda je:

V2A=V1(cos⁡θ12+sin⁡θ12tg⁡φ2);1,05⋅sin⁡θ12⋅tg⁡φ2=0,95A−1,05⋅cos⁡θ12.\begin{gathered} V_2A = V_1\left(\cos\theta_{12} + \sin\theta_{12}\operatorname{tg}\varphi_2\right); \\ 1{,}05\cdot\sin\theta_{12}\cdot\operatorname{tg}\varphi_2 = 0{,}95A - 1{,}05\cdot\cos\theta_{12}. \end{gathered}

Teorijski granični ugao statičke stabilnosti iz izraza Pp=V1V2Xvsin⁡θ12P_p = \dfrac{V_1V_2}{X_v}\sin\theta_{12}, može se dobiti iz uslova ∂Pp∂θ12=V1V2Xvcos⁡θ12=0\dfrac{\partial P_p}{\partial\theta_{12}} = \dfrac{V_1V_2}{X_v}\cos\theta_{12} = 0 i može se grafički predstaviti kao:

Snaga Pp u funkciji ugla θ12: kriva (V1V2/Xv) sin θ12 sa maksimumom PpM = V1V2/Xv pri uglu π/2.

Sa prethodne slike se vidi da je teorijska granica statičke stabilnosti θ12stgr=π/2\theta_{12st}^{gr} = \pi/2.

Praktični granični ugao statičke stabilnosti zavisi od zahtevane rezerve stabilnosti (RS):

sin⁡θ12stgr=1−(RS),\sin\theta_{12st}^{gr} = 1 - (RS),

gde su pojedine vrednosti date u donjoj tabeli:

(RS)(RS) [r.j.]0,10,150,200,250,30
θ12stgr\theta_{12st}^{gr} [stepeni]64,1658,2153,1348,5944,43

Dozvoljeni pad napona još više smanjuje taj granični ugao stabilnosti.

Ovo je naročito izraženo kod vodova srednje dužine.

VnV_n [kV]60110220400
Pnat=Vn2ZcP_{nat} = \dfrac{V_n^2}{Z_c} [MW]930,2121457
ZcZ_c [Ω\Omega]400400400350
Broj provodnika po fazi1112
Odnos Pp/Pnat u funkciji dužine voda Lv [km], od 200 do 800 km: kriva pada od blizu 3 pri najkraćim vodovima do ispod 1 pri najdužim.

Generalni zaključci:

  • Za kratke vodove (do 100 km) merodavna je termička granica opteretivosti.
  • Za vodove srednje dužine (od 100 do 500 km) merodavna je granica određena padom napona (i granicom naponske stabilnosti).
  • Za duge vodove (preko 500 km) granica moguće snage prenosa određena je rezervom statičke stabilnosti.

Studija dogradnje prenosne mreže

Mora se posmatrati ceo EES. Međutim, za razliku od izvora, model je složeniji, pošto se uvažavaju: topologija mreže, koridori vodova, lokacije postrojenja i slično.

Ograničenja koja se javljaju u optimizaciji mogu biti vrlo složena i često nelinearna: naponska ograničenja, tranzijentna stabilnost (diferencijalne jednačine).

Osnovni principi pri planiranju prenosnih mreža su:

  • Mora se uvažiti simultanost proizvodnje i potrošnje.
  • Moraju se unapred definisati norme kvaliteta električne energije (određeno propisima).
  • Koordinacija rada proizvodnih i prenosnih kapaciteta.
  • Dobra analiza ekonomskih efekata (inženjerska ekonomija).

Postupci planiranja:

  • Statički (rešavaju problem za neku godinu iz horizonta planiranja, ne vodeći računa o načinu prelaska u nju) – nema vremenskog plana izrade.
  • Dinamički (uključuju se i vremenski planovi razvoja; generalno se koristi statički model za horizont godinu, a onda se proveravaju prethodne godine). Ako bi se postupilo suprotno (proširenje od prve do horizont godine), može se dogoditi da se mreža dograđuje sa više malih elemenata, što je neekonomično. Na primer, izgradnja dva voda 220 kV, umesto jednog 400 kV.

Statički kriterijumi dogradnje prenosne mreže:

  • Ne sme biti vodova opterećenih iznad termičke granice, i drugih elemenata iznad nominalnog opterećenja.
  • Dozvoljeni naponi u opsegu ±5 % (izuzetno −10 %).
  • Nekompenzovani ugao prenosne snage vodova ne sme biti veći od 45°.
  • Mreža mora zadovoljiti za (n−1)(n - 1) ispade, uz smanjenje napona od 5 do 10 %.

Postupci koji se koriste u planiranju:

  • Provera opterećenja: U testiranju topološke šeme prvo se mora proveriti da li grane imaju dovoljan prenosni kapacitet. Zato se sprovodi provera tokova snaga za osnovno stanje i (n−1)(n - 1) ispad. Najčešće se koriste paketi za analizu sigurnosti, zasnovani na DC modelu.
  • Analiza osetljivosti: Ako se utvrdi preopterećenje neke grane, mora se odrediti dogradnja kog elementa najbolje doprinosi otklanjanju tog preopterećenja. Zato se koristi analiza osetljivosti. Kriterijum je minimum specifičnih investicija (po km izgrađenog voda) za dogradnju novog voda. Kod transformatora mora se dograditi nova jedinica.
  • Formiranje topološke šeme mreže: Pravi se lista opcija dogradnje prenosne mreže sa dinamikom dogradnje. Bira se najbolja opcija.

Dakle, osnovni metod je analiza sigurnosti i DC model proračuna tokova snaga, koji se izlažu u nastavku.

DC model tokova snaga

Opšti izraz za jednačinu aktivnih injektiranja je:

Pi=Pgi−Ppi=Vi2Gii+Vi∑j∈αi[Vj(Gijcos⁡θij+Bijsin⁡θij)],P_i = P_{gi} - P_{pi} = V_i^2G_{ii} + V_i\sum_{j \in \alpha_i}\left[V_j\left(G_{ij}\cos\theta_{ij} + B_{ij}\sin\theta_{ij}\right)\right],

gde su GijG_{ij} i BijB_{ij} elementi matrice admitansi Y‾BUS\ph{\boldsymbol{Y}}_{BUS}:

Y‾ij={Y‾i0+∑j∈αiY‾ijg,i=j−Y‾ijg,i≠j;\ph{Y}_{ij} = \begin{cases} \ph{Y}_i^0 + \sum_{j \in \alpha_i}\ph{Y}_{ij}^g, & i = j \\ -\ph{Y}_{ij}^g, & i \ne j \end{cases};
  • αi\alpha_i – skup čvorova incidentnih i-tom čvoru.

Osnovna uprošćenja koja se uvode u gornju jednačinu za dobijanje DC modela tokova snaga su:

Gij≈0; Vi=Vj≈1,00 r.j.; – nema potrebe za Q-V delom,sin⁡θij≈θij; cos⁡θij≈1,00 r.j.\begin{gathered} G_{ij} \approx 0; \ V_i = V_j \approx 1{,}00\ \rj; \ \text{– nema potrebe za Q-V delom}, \\ \sin\theta_{ij} \approx \theta_{ij}; \ \cos\theta_{ij} \approx 1{,}00\ \rj \end{gathered}

Sa napred uvedenim aproksimacijama jednačina injektiranja je:

Pi=∑j∈αiθijXijg=−∑j∈αiBijgθij=−∑j∈αiBijg(θi−θj)=−θi∑j∈αiBijg+∑j∈αiBijgθj=−Biiθi−∑j∈αiBijθj;i=1,2,⋯ ,N;i∉αSL,\begin{gathered} \begin{aligned} P_i &= \sum_{j \in \alpha_i}\frac{\theta_{ij}}{X_{ij}^g} = -\sum_{j \in \alpha_i}B_{ij}^g\theta_{ij} = -\sum_{j \in \alpha_i}B_{ij}^g\left(\theta_i - \theta_j\right) \\ &= -\theta_i\sum_{j \in \alpha_i}B_{ij}^g + \sum_{j \in \alpha_i}B_{ij}^g\theta_j = -B_{ii}\theta_i - \sum_{j \in \alpha_i}B_{ij}\theta_j; \end{aligned} \\ i = 1, 2, \cdots, N; \quad i \notin \alpha_{SL}, \end{gathered}

jer je: Bii=∑j∈αiBijgB_{ii} = \sum_{j \in \alpha_i}B_{ij}^g, Bij=−BijgB_{ij} = -B_{ij}^g.

Dakle, redukovana [(N−1)(N - 1)-dimenziona, bez balansnog čvora] matrična jednačina injektiranja je (pošto je jednačina sa balansnim čvorom singularna):

Pr=−Brθr,\boldsymbol{P_r} = -\boldsymbol{B_r}\boldsymbol{\theta_r},

gde je:

Bij={∑j∈αiBijg,i=j−Bijg,i≠jB_{ij} = \begin{cases} \sum_{j \in \alpha_i}B_{ij}^g, & i = j \\ -B_{ij}^g, & i \ne j \end{cases}

(uz pretpostavku da je Bi0=0B_i^0 = 0, što ne utiče na generalnost).

Za referentno-balansni čvor je:

PSL=−∑i=1i≠αSLNPi.P_{SL} = -\sum_{\substack{i=1 \\ i \ne \alpha_{SL}}}^N P_i.

Iz prethodne jednačine injektiranja je:

θr=−Br−1Pr,\boldsymbol{\theta_r} = -\boldsymbol{B_r}^{-1}\boldsymbol{P_r},

što znači da se rešenje dobija u jednom koraku (bez iterativnog procesa), pri čemu se usvaja θSL=0\theta_{SL} = 0.

Iz dobijenih uglova fazora napona mogu se izračunati tokovi snaga po granama:

Pij=θi−θjXijg=−Bijg(θi−θj);i=1,2,⋯ ,N;j∈αi,P_{ij} = \frac{\theta_i - \theta_j}{X_{ij}^g} = -B_{ij}^g(\theta_i - \theta_j); \quad i = 1, 2, \cdots, N; \quad j \in \alpha_i,

ili u matričnoj formi:

Pℓ=−BℓgArTθr=BℓgArTBr−1Pr,\boldsymbol{P}_\ell = -\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{\theta_r} = \boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{B_r}^{-1}\boldsymbol{P_r},

gde je:

  • Pℓ\boldsymbol{P}_\ell – LL-dimenzioni vektor tokova aktivnih snaga PℓP_\ell po granama (LL je ukupan broj grana);
  • Bℓg=diag{Bℓg}\boldsymbol{B}_\ell^g = \mathrm{diag}\left\{B_\ell^g\right\} – (L×L)(L \times L)-dimenziona dijagonalna matrica susceptansi grana;
  • θr\boldsymbol{\theta_r} – redukovani (N−1)(N - 1)-dimenzioni vektor faznih uglova napona u čvorovima;
  • Ar={Aij}\boldsymbol{A_r} = \left\{A_{ij}\right\} – [(N−1)×L][(N - 1) \times L]-dimenziona matrica incidencije čvorova i grana (Aij=±1A_{ij} = \pm 1 ako je ℓ\ell-ta grana povezana sa čvorom, odnosno Aij=0A_{ij} = 0 ako ℓ\ell-ta grana nije povezana sa čvorom).

Analiza osetljivosti

DC model je linearan model, pa je pogodan za linearnu analizu osetljivosti.

Faktori preraspodele tokova snaga po granama pri promeni injektiranja

ΔPr=−BrΔθr – ne menja se topologija i matrica Br;ΔPr=Pr(1)−Pr(0);Δθr=θr(1)−θr(0).\begin{gathered} \Delta\boldsymbol{P_r} = -\boldsymbol{B_r}\Delta\boldsymbol{\theta_r} \ \text{– ne menja se topologija i matrica}\ \boldsymbol{B_r}; \\ \Delta\boldsymbol{P_r} = \boldsymbol{P_r}^{(1)} - \boldsymbol{P_r}^{(0)}; \\ \Delta\boldsymbol{\theta_r} = \boldsymbol{\theta_r}^{(1)} - \boldsymbol{\theta_r}^{(0)}. \end{gathered}

Faktor preraspodele akia_k^i toka aktivne snage u k-toj grani usled promene injektiranja u i-tom čvoru je:

aki=ΔPkΔPi=1Xkg(Δθm−ΔθnΔPi),a_k^i = \frac{\Delta P_k}{\Delta P_i} = \frac{1}{X_k^g}\left(\frac{\Delta\theta_m - \Delta\theta_n}{\Delta P_i}\right),

gde je k-ta grana između m-tog i n-tog čvora.

U prethodnom izrazu ΔPi\Delta P_i je, na primer, ispad generatorske ili potrošačke snage.

Tada je novi tok snage po k-toj grani:

Pk(1)=Pk(0)+ΔPk=Pk(0)+akiΔPi.P_k^{(1)} = P_k^{(0)} + \Delta P_k = P_k^{(0)} + a_k^i\Delta P_i.

Pošto je model linearan, onda važi princip superpozicije, pa se model može lako primeniti na slučaj višestrukih promena u vektoru injektiranja:

Pk(1)=Pk(0)+∑i∈αpakiΔPi,P_k^{(1)} = P_k^{(0)} + \sum_{i \in \alpha_p}a_k^i\Delta P_i,

gde je αp\alpha_p skup čvorova mreže gde je došlo do promene injektiranja.

Izvođenje faktora osetljivosti

ΔPr=−BrΔθr;Δθr=−Br−1ΔPr,ΔPℓ=−BℓgΔθℓ=−BℓgArTΔθr=BℓgArTBr−1ΔPr,\begin{gathered} \Delta\boldsymbol{P_r} = -\boldsymbol{B_r}\Delta\boldsymbol{\theta_r}; \\ \Delta\boldsymbol{\theta_r} = -\boldsymbol{B_r}^{-1}\Delta\boldsymbol{P_r}, \\ \Delta\boldsymbol{P}_\ell = -\boldsymbol{B}_\ell^g\Delta\boldsymbol{\theta}_\ell = -\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\Delta\boldsymbol{\theta_r} = \boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{B_r}^{-1}\Delta\boldsymbol{P_r}, \end{gathered}

odnosno:

aki={ΔPkΔPi}=BℓgArTBr−1,\boxed{a_k^i = \left\{\frac{\Delta P_k}{\Delta P_i}\right\} = \boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{B_r}^{-1}},

gde se može definisati matrica BℓgArTBr−1\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{B_r}^{-1} kao matrica osetljivosti promene tokova snaga po granama na promene injektiranja u čvorovima.

Distribucioni faktori preraspodele tokova snaga po granama pri uključenju/isključenju grane

Ako dolazi do ispada grane između m-tog i n-tog čvora, onda se menja matrica Br\boldsymbol{B_r}. Tada je promena uglova fazora napona:

Δθr=−ΔBr−1Pr;ΔBr=Br(1)−Br(0)=−BmngekekT,\begin{gathered} \Delta\boldsymbol{\theta_r} = -\Delta\boldsymbol{B_r}^{-1}\boldsymbol{P_r}; \\ \Delta\boldsymbol{B_r} = \boldsymbol{B_r}^{(1)} - \boldsymbol{B_r}^{(0)} = -B_{mn}^g\boldsymbol{e}_k\boldsymbol{e}_k^{\mathrm{T}}, \end{gathered}

gde je:

ek=[01⋮−1⋮0]←m-ta pozicija⋮←n-ta pozicija⋮.\boldsymbol{e}_k = \begin{bmatrix} 0 \\ 1 \\ \vdots \\ -1 \\ \vdots \\ 0 \end{bmatrix} \begin{matrix} \\ \leftarrow m\text{-ta pozicija} \\ \vphantom{\vdots} \\ \leftarrow n\text{-ta pozicija} \\ \vphantom{\vdots} \\ \\ \end{matrix}.
Δθr=θr(1)−θr(0);θr(1)=θr(0)+Δθr=−(Br(1))−1Pr=−(Br(0)+ΔBr)−1Pr,\begin{gathered} \Delta\boldsymbol{\theta_r} = \boldsymbol{\theta_r}^{(1)} - \boldsymbol{\theta_r}^{(0)}; \\ \boldsymbol{\theta_r}^{(1)} = \boldsymbol{\theta_r}^{(0)} + \Delta\boldsymbol{\theta_r} = -\left(\boldsymbol{B_r}^{(1)}\right)^{-1}\boldsymbol{P_r} = -\left(\boldsymbol{B_r}^{(0)} + \Delta\boldsymbol{B_r}\right)^{-1}\boldsymbol{P_r}, \end{gathered}

ili:

θr(1)=−(Br(0)−BmngekekT)−1Pr.\boxed{\boldsymbol{\theta_r}^{(1)} = -\left(\boldsymbol{B_r}^{(0)} - B_{mn}^g\boldsymbol{e}_k\boldsymbol{e}_k^{\mathrm{T}}\right)^{-1}\boldsymbol{P_r}}.

Međutim, prethodnu formulu je neophodno ponovo invertovati za svaki ispad (promenu topologije), što je numerički nepovoljno.

Zato se koristi Lema o inverziji matrice (Šerman-Morisonova lema):

(A+CDT)−1=A−1−A−1C(Im+DTA−1C)−1DTA−1.\left(\boldsymbol{A} + \boldsymbol{C}\boldsymbol{D}^{\mathrm{T}}\right)^{-1} = \boldsymbol{A}^{-1} - \boldsymbol{A}^{-1}\boldsymbol{C}\left(\boldsymbol{I_m} + \boldsymbol{D}^{\mathrm{T}}\boldsymbol{A}^{-1}\boldsymbol{C}\right)^{-1}\boldsymbol{D}^{\mathrm{T}}\boldsymbol{A}^{-1}.

Ovde je poznato (ranije rešeno za stanje bez ispada) A−1\boldsymbol{A}^{-1} (u našim uslovima (Br(0))−1\left(\boldsymbol{B_r}^{(0)}\right)^{-1}), dok je izraz u zagradi dijagonalna matrica, pa je lak za inverziju:

C=c=−Bmngek;D=d=ek;A=Br(0).\begin{gathered} \boldsymbol{C} = \boldsymbol{c} = -B_{mn}^g\boldsymbol{e}_k; \\ \boldsymbol{D} = \boldsymbol{d} = \boldsymbol{e}_k; \\ \boldsymbol{A} = \boldsymbol{B_r}^{(0)}. \end{gathered}

Kada se prethodni izraz sredi, može se napisati kao:

θr(1)=θr(0)+(Br(0))−1ek(θm(0)−θn(0))1Bmng−ekT(Br(0))−1ek.\boxed{\boldsymbol{\theta_r}^{(1)} = \boldsymbol{\theta_r}^{(0)} + \frac{\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k\left(\theta_m^{(0)} - \theta_n^{(0)}\right)}{\dfrac{1}{B_{mn}^g} - \boldsymbol{e}_k^{\mathrm{T}}\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k}}.

Koristeći prethodno izvedeni izraz za uglove fazora napona posle ispada, dobijaju se tokovi snaga kao:

Pℓ(1)=−BℓgArTθr(1)=BℓgArT(Br(0)−BmngekekT)−1Pr,\boldsymbol{P}_\ell^{(1)} = -\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{\theta_r}^{(1)} = \boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\left(\boldsymbol{B_r}^{(0)} - B_{mn}^g\boldsymbol{e}_k\boldsymbol{e}_k^{\mathrm{T}}\right)^{-1}\boldsymbol{P_r},

odnosno:

ΔPℓk=−BℓgArTΔθrk=−BℓgArT(Br(0))−1ek(θm(0)−θn(0))1Bmng−ekT(Br(0))−1ek,\Delta\boldsymbol{P}_\ell^k = -\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\Delta\boldsymbol{\theta_r}^k = -\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\frac{\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k\left(\theta_m^{(0)} - \theta_n^{(0)}\right)}{\dfrac{1}{B_{mn}^g} - \boldsymbol{e}_k^{\mathrm{T}}\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k},

pa je:

dℓk={ΔPℓkPk(0)},Pℓ(1)=Pℓ(0)+dℓkPk(0),\boxed{d_\ell^k = \left\{\frac{\Delta P_\ell^k}{P_k^{(0)}}\right\}}, \quad \boldsymbol{P}_\ell^{(1)} = \boldsymbol{P}_\ell^{(0)} + d_\ell^kP_k^{(0)},

gde je A\boldsymbol{A} (N×L)(N \times L)-dimenziona matrica incidencije čvorova i grana (Aij=±1A_{ij} = \pm 1 ako je ℓ\ell-ta grana povezana sa čvorom, odnosno Aij=0A_{ij} = 0 ako ℓ\ell-ta grana nije povezana sa čvorom), odnosno:

A=[⋮⋯1⋯⋮⋯−1⋯⋮]→ Grane (L)↓cˇvorovi.\boldsymbol{A} = \overset{\rightarrow\ \text{Grane}\ (L)}{\begin{bmatrix} & \vdots & \\ \cdots & 1 & \cdots \\ & \vdots & \\ \cdots & -1 & \cdots \\ & \vdots & \end{bmatrix}} \downarrow \text{čvorovi}.

Osetljivost promene tokova snaga po granama pri simultanom ispadu grane i promeni injektiranja u čvorovima

Ovo je najteži slučaj, ali ujedno i najvažniji. Teškoća u primeni postupka je što je veći poremećaj, veća je i greška linearizacije.

U izvođenju osetljivosti polazi se od:

ΔPℓ=Pℓ(1)−Pℓ(0)=aℓiΔPi+dℓkPk(1),\Delta P_\ell = P_\ell^{(1)} - P_\ell^{(0)} = a_\ell^i\Delta P_i + d_\ell^kP_k^{(1)},

gde je snaga koja bi tekla po grani k≠ℓk \ne \ell posle ispada generatora u i-tom čvoru (ΔPi=−Pgi\Delta P_i = -P_{gi}):

Pk(1)=Pk(0)+ΔPk=Pk(0)+akiΔPi=Pk(0)−akiPgi.P_k^{(1)} = P_k^{(0)} + \Delta P_k = P_k^{(0)} + a_k^i\Delta P_i = P_k^{(0)} - a_k^iP_{gi}.

Sada je:

ΔPℓ=dℓkPk(0)+(aℓi+dℓkaki)ΔPi;ΔPℓ=dℓkPk(0)−(aℓi+dℓkaki)Pgi=aℓi,kPgi.\begin{gathered} \Delta P_\ell = d_\ell^kP_k^{(0)} + \left(a_\ell^i + d_\ell^ka_k^i\right)\Delta P_i; \\ \Delta P_\ell = d_\ell^kP_k^{(0)} - \left(a_\ell^i + d_\ell^ka_k^i\right)P_{gi} = a_\ell^{i,k}P_{gi}. \end{gathered}

Tok snage posle simultanog ispada je:

Pℓ(1)=Pℓ(0)+ΔPℓ=Pℓ(0)+aℓi,kPgi,P_\ell^{(1)} = P_\ell^{(0)} + \Delta P_\ell = P_\ell^{(0)} + a_\ell^{i,k}P_{gi},

gde je, dakle, za simultani ispad i-tog generatora i k-te grane:

aℓi,k=dℓkPk(0)Pgi−(aℓi+dℓkaki);ℓ=1,2,⋯ ,L≠k.\boxed{a_\ell^{i,k} = d_\ell^k\frac{P_k^{(0)}}{P_{gi}} - \left(a_\ell^i + d_\ell^ka_k^i\right)}; \quad \ell = 1, 2, \cdots, L \ne k.

Primena metoda osetljivosti za planiranje mreže

Vrlo je važna osetljivost tokova snaga na dodavanje novih vodova. Prvi metod koji se obično primenjuje je metod sukcesivne ekspanzije, koji sudi o efektivnosti novog voda-kandidata za izgradnju na osnovu rasterećenja ostalih vodova. Vodovi se dodaju po svojoj efektivnosti.

Metod sukcesivne redukcije je suprotan od prvog. Dodaju se svi potencijalni vodovi, a onda se eliminišu jedan po jedan po kriterijumu najmanje efektivnosti.

Metod sukcesivne ekspanzije

Indeks efektivnosti za metod sukcesivne ekspanzije dobija se iz sledećeg izvođenja (ranije izveden izraz pri dodavanju k-te grane, između m-tog i n-tog čvora):

Δθr=−(Br(0))−1ek(θm(0)−θn(0))1Bmng+ekT(Br(0))−1ek,\Delta\boldsymbol{\theta_r} = -\frac{\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k\left(\theta_m^{(0)} - \theta_n^{(0)}\right)}{\dfrac{1}{B_{mn}^g} + \boldsymbol{e}_k^{\mathrm{T}}\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k},

gde se u prethodnom izrazu mogu uvesti sledeće oznake:

Xkg=−Xmng=1Bmng;Xkd=ekT(Br(0))−1ek;βk=−1Xkg+Xkd;θk(0)=θm(0)−θn(0),\begin{aligned} X_k^g &= -X_{mn}^g = \frac{1}{B_{mn}^g}; \\ X_k^d &= \boldsymbol{e}_k^{\mathrm{T}}\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{e}_k; \\ \beta_k &= -\frac{1}{X_k^g + X_k^d}; \\ \theta_k^{(0)} &= \theta_m^{(0)} - \theta_n^{(0)}, \end{aligned}

pa je:

Δθr=βkXr(0)ekθk(0).\Delta\boldsymbol{\theta_r} = \beta_k\boldsymbol{X_r}^{(0)}\boldsymbol{e}_k\theta_k^{(0)}.

Elementi matrice Xr(0)\boldsymbol{X_r}^{(0)} se teško računarski određuju (zbog inverzije). Ovaj pristup se može primeniti za sukcesivnu ekspanziju, ali se ne može primeniti za sukcesivnu redukciju, pošto tada βk→∞\beta_k \to \infty.

Indeks efektivnosti definiše se kao:

Uticaj dodavanja k-te grane (između m-tog i n-tog čvora) na ugao fazora napona ℓ\ell-te grane:

Eℓk=ΔθℓkCIvk=eℓTΔθrCIvk=βkeℓTXr(0)ekθk(0)CIvk;k∈αs.E_\ell^k = \frac{\Delta\theta_\ell^k}{C_{Iv}^k} = \frac{\boldsymbol{e}_\ell^{\mathrm{T}}\Delta\boldsymbol{\theta_r}}{C_{Iv}^k} = \frac{\beta_k\boldsymbol{e}_\ell^{\mathrm{T}}\boldsymbol{X_r}^{(0)}\boldsymbol{e}_k\theta_k^{(0)}}{C_{Iv}^k}; \quad k \in \alpha_s.

Prethodni izraz kaže da se analizira odnos rasterećenja vodova (tokova snaga, odnosno razlike uglova) i investicija.

Posle testiranja svih vodova – kandidata za izgradnju bira se najpovoljniji vod iz ugla rasterećenja ℓ\ell-te grane:

max⁡k∈αs{Eℓk}; k∈αs.\max_{k \in \alpha_s}\left\{E_\ell^k\right\}; \ k \in \alpha_s.

Ako u sistemu postoji više preopterećenih vodova (αℓ\alpha_\ell), onda se računa globalni indeks efektivnosti:

Ek=∑ℓ∈αℓEℓk=βkθk(0)CIvk∑ℓ∈αℓeℓTXr(0)ek.E^k = \sum_{\ell \in \alpha_\ell}E_\ell^k = \frac{\beta_k\theta_k^{(0)}}{C_{Iv}^k}\sum_{\ell \in \alpha_\ell}\boldsymbol{e}_\ell^{\mathrm{T}}\boldsymbol{X_r}^{(0)}\boldsymbol{e}_k.

Čitav algoritam proširenja prenosne mreže sukcesivnim izborom najefikasnijeg voda od svih kandidata za izgradnju sprovodi se kroz sledeće korake:

  1. Ulazni podaci za DC proračun tokova snaga (topologija, reaktanse grana, injektiranja...).
  2. Iz Br\boldsymbol{B_r} dobija se vektor faznih uglova θ=[θSLθr]\boldsymbol{\theta} = \begin{bmatrix} \theta_{SL} & \boldsymbol{\theta_r} \end{bmatrix}, gde je θr=−(Br(0))−1Pr\boldsymbol{\theta_r} = -\left(\boldsymbol{B_r}^{(0)}\right)^{-1}\boldsymbol{P_r}.
  3. Proračun tokova aktivnih snaga po granama Pℓ=−BℓgArTθr\boldsymbol{P}_\ell = -\boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{\theta_r}, ili Pℓ=BℓgArTBr−1Pr\boldsymbol{P}_\ell = \boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{B_r}^{-1}\boldsymbol{P_r}.
  4. Identifikacija preopterećenja grana ∣Pℓ∣≤PℓM\left|\boldsymbol{P}_\ell\right| \le \boldsymbol{P}_{\ell M}.
  5. Vodovi koji ne zadovolje uslov iz Tačke 4, svrstavaju se u skup αℓ\alpha_\ell.
  6. Proračun globalnog indeksa efektivnosti za svaki vod koji je kandidat za izgradnju EkE^k, k∈αsk \in \alpha_s (αs\alpha_s je skup vodova kandidata za izgradnju).
  7. Bira se najefektivniji vod po kriterijumu max⁡k∈αs{Ek}\max\limits_{k \in \alpha_s}\left\{E^k\right\}.

Metod sukcesivne redukcije

Suprotno od prethodnog metoda, prvo se formira fiktivna mreža, koja obuhvata postojeću mrežu i dodate sve vodove-kandidate za izgradnju. Ova fiktivna mreža je jako isprepletana, redundantna i neekonomična. Onda se sukcesivno eliminišu vodovi po kriterijumu:

Ek=∣Pk∣CIvk;k∈αs,E^k = \frac{\left|P_k\right|}{C_{Iv}^k}; \quad k \in \alpha_s,

što znači da kriterijum valorizuje odnos maksimalnih prenosnih snaga i investicionih troškova za skup analiziranih vodova αs\alpha_s.

Vidi se da se ovde koristi različit indeks efektivnosti u odnosu na metod sukcesivne ekspanzije.

Kriterijum eliminacije pojedinih opcija je:

min⁡k∈αs{Ek}.\min_{k \in \alpha_s}\left\{E^k\right\}.

Programski paketi za planiranje prenosnih mreža

Kod nas se najčešće koristi programski paket TransPlan (Transmission Planning). Razvijen je u belgijskoj firmi System Europe i isporučen bivšoj SFRJ, pa se zato i danas koristi u Srbiji (sa znatnim unapređenjima).

Povezan je sa CLF-OPF (Constrained Load Flow – Optimal Power Flow) programom, čime je umesto DC modela tokova snaga omogućeno korišćenje AC modela.

Koristi linearan planerski model, primenom dualnog simpleks algoritma.

Primena linearnog programiranja za planiranje prenosnih mreža

Prethodna dva metoda su korak po korak (sukcesivni), pa su praktično delimično heuristični, odnosno ne garantuju globalni optimum. Optimum po koracima ne garantuje globalni optimum. Egzaktnije rešenje je primena optimizacionih tehnika, od kojih je najveću primenu našlo linearno programiranje.

Model postojeće mreže (označen indeksom P)

Osnovna jednačina je linearizovana jednačina tokova snaga:

P(P)=−B(P)θ(P), ili Pr(P)=−Br(P)θr(P).\boldsymbol{P}_{(P)} = -\boldsymbol{B}_{(P)}\boldsymbol{\theta}_{(P)}, \ \text{ili} \ \boldsymbol{P}_{\boldsymbol{r}(P)} = -\boldsymbol{B}_{\boldsymbol{r}(P)}\boldsymbol{\theta}_{\boldsymbol{r}(P)}.

uz uvažavanje ograničenja tokova snaga po granama:

∣Pℓ(P)∣≤Pℓ(P)M,\left|\boldsymbol{P}_{\ell(P)}\right| \le \boldsymbol{P}_{\ell(P)}^M,

ili u matričnoj formi:

∣Bℓ(P)gA(P)Tθ(P)∣≤Pℓ(P)M.\left|\boldsymbol{B}_{\ell(P)}^g\boldsymbol{A}_{(P)}^{\mathrm{T}}\boldsymbol{\theta}_{(P)}\right| \le \boldsymbol{P}_{\ell(P)}^M.

Iz prethodnog se vidi da se ne mora koristiti redukovani sistem jednačina DC tokova snaga, pošto nema inverzije matrice.

Model buduće proširene mreže (postojeći i novododati vodovi)

P(D)=A(D)Pℓ(D).\boldsymbol{P}_{(D)} = \boldsymbol{A}_{(D)}\boldsymbol{P}_{\ell(D)}.

U prethodnom izrazu indeks D označava skup novih, dodatih čvorova proširene mreže (novi izvori).

A(D)\boldsymbol{A}_{(D)} je (ND×L)(N_D \times L)-dimenziona matrica incidencije dodatnih čvorova i svih grana (L=LP+LDL = L_P + L_D) proširene mreže.

Novi vodovi (kandidati za izgradnju) obično nemaju prenosnih ograničenja.

Ako su investicioni troškovi srazmerni tokovima aktivnih snaga, onda je optimizacioni kriterijum:

min⁡Pℓ(P+D){Z}=CIv(P+D)T∣Pℓ(P+D)∣,\min_{\boldsymbol{P}_{\ell(P+D)}}\{Z\} = \boldsymbol{C}_{Iv(P+D)}^{\mathrm{T}}\left|\boldsymbol{P}_{\ell(P+D)}\right|,

uz uvažavanje ograničenja:

A(P)Pℓ(P+D)⏟P(P)+A(D)Pℓ(P+D)⏟P(D)=P(P+D).∣Pℓ(P+D)∣≤Pℓ(P+D)M.\begin{gathered} \underbrace{\boldsymbol{A}_{(P)}\boldsymbol{P}_{\ell(P+D)}}_{\boldsymbol{P}_{(P)}} + \underbrace{\boldsymbol{A}_{(D)}\boldsymbol{P}_{\ell(P+D)}}_{\boldsymbol{P}_{(D)}} = \boldsymbol{P}_{(P+D)}. \\ \left|\boldsymbol{P}_{\ell(P+D)}\right| \le \boldsymbol{P}_{\ell(P+D)}^M. \end{gathered}

U prethodnom optimizacionom problemu je CIv(P)=0\boldsymbol{C}_{Iv(P)} = 0, pošto nema investicija na postojeće vodove.

Ograničenje tipa jednakosti u prethodnoj formulaciji predstavlja jednačinu bilansa snaga injektiranja po metodu superpozicije, dok je ograničenje tipa nejednakosti ograničenje tokova snaga po granama (novi vodovi obično nisu kritični sa stanovišta ovog ograničenja).

Kada se rešava prethodni problem može doći do dve karakteristične situacije:

  • Z=0Z = 0, što znači da nema preopterećenja vodova, pa nisu potrebni novi vodovi (investicije su nula). Ovaj uslov praktično drže ograničenja.
  • Z>0Z > 0, što znači da postoje preopterećenja nekih vodova. Tada treba izdvojiti tokove snaga po dodatnim vodovima Pℓ(D)\boldsymbol{P}_{\ell(D)} iz ukupnog vektora Pℓ(P+D)\boldsymbol{P}_{\ell(P+D)}. U ovom slučaju najefektivniji vod je onaj sa najvećom prenosnom snagom. Njegovo dodavanje na postojeću mrežu daje najveće smanjenje preopterećenja (najčešće).

Primena drugih metoda matematičke optimizacije za planiranje prenosnih mreža

Svaki optimizacioni problem čine četiri osnovna elementa:

  • Optimizacioni kriterijum (funkcija cilja).
  • Promenljive.
  • Ograničenja.
  • Tehnika rešavanja.

Optimizacioni kriterijum je funkcija promenljivih odlučivanja. Obično u sebe uključuje investicije za proširenje mreže i pogonske troškove. Kod dugoročnog planiranja mora se izvršiti aktualizacija troškova.

Promenljive mogu biti: odlučivanja i stanja.

Promenljive odlučivanja su celobrojne i određuju da li će se neki vod graditi, čime utiču na topologiju.

Promenljive stanja određuju režim eksploatacije EES-a (naponi, uglovi fazora napona, tokovi snaga po granama, injektiranja i druge).

Ograničenja uvažavaju konstrukcione i pogonske karakteristike opreme, granice promenljivih stanja i odlučivanja i sl. Najčešće su to preopterećenja vodova i jednačine tokova snaga. Ređe su to naponska ograničenja i ograničenja stabilnosti.

Tehnika rešavanja zavisi od formulacije problema. To su najčešće: nelinearno programiranje, linearno programiranje, dinamičko programiranje, mešovito-celobrojno programiranje, tehnika grananja i ograničavanja, metode konveksne optimizacije, metode veštačke inteligencije i druge.

Standardizacija elemenata prenosne mreže

Do sada je rađena pojedinačna optimizacija: vodova, transformatora, poprečnog preseka, napona i sl. Međutim, optimizacija pojedinih elemenata ne daje globalni optimum čitavog EES-a. Ako bi ostavili optimalne pojedinačne parametre imali bi više (mnogo) naponskih nivoa i poprečnih preseka. To je nepovoljno sa stanovišta održavanja i rezervnih delova. Zbog serijske proizvodnje, povoljno je da postoji manje standardnih elemenata.

Postoji više međunarodnih i nacionalnih standarda: JUS, VDE, IEC, NEMA, IEEE (ANSI) i drugi.

Standardni poprečni preseci (postoji najčešće 1 do 2 po naponskom nivou):

  • 110 kV: 150/25 Ač i 240/40 Ač.
  • 220 kV: 360/57 Ač i 490/65 Ač.
  • 400 kV: 2×490/65.

Standardne snage transformatora 110 kV/SN su:

20 MVA, 31,5 MVA i 63 MVA.

Standardne snage autotransformatora:

  • 220 kV/110 kV: 150 MVA i 250 MVA.
  • 400 kV/110 kV: 300 MVA.
  • 400 kV/220 kV: 400 MVA.

Čini se da praksa ruši ekonomske efekte optimizacije. U praksi to obično nije slučaj, pošto elementi nisu previše osetljivi na promenu parametara. Efekat udaljavanja od optimuma ne menja previše rezultate.

Provera izabranog plana razvoja sistema

Kada se izabere optimalan plan razvoja treba pristupiti proveri tehničke adekvatnosti. Zato se vrše detaljniji proračuni za karakteristična rešenja:

  • Proračun tokova snaga i naponskih stanja uz korišćenje kompletnog AC modela (najčešće se koriste brzi raspregnuti algoritmi). Ovako se proverava adekvatnost topologije i parametara mreže.
  • Proračun nivoa struja i snaga simetričnih i nesimetričnih kratkih spojeva i zemljospojeva. Ovako se proverava adekvatnost izbora rasklopnih aparata, relejne zaštite, napona koraka i dodira i uticaja energetskih na telekomunikacione vodove.
  • Analiza statičke sigurnosti pri jednostrukim (n−1)(n - 1) i ređe dvostrukim (n−2)(n - 2) poremećajima.
  • Studija naponsko-reaktivnih prilika, kompenzacije i naponskih regulacionih resursa.
  • Studija statičke i tranzijentne stabilnosti. Određuje tipove i karakteristike zaštite vodova i automatskog ponovnog uključenja.
  • Studija regulacije učestanosti i aktivnih snaga. Izbor tipova primarne i sekundarne regulacije, karakteristike turbinskih regulatora i sekundarnih regulatora učestanosti i snage razmene.
  • Izbor sistema relejne zaštite generatora, vodova i transformatora.
  • Studija sistema veza i informatike (SCADA).
  • Studija pouzdanosti. Provera kriterijuma isporuke (LOLP, LOEP).

Studija dogradnje reaktivnih izvora

Realizuje se za optimalnu strategiju razvoja prenosne mreže. Ova situacija sprovodi se ako su naponsko-reaktivne prilike nezadovoljavajuće.

Ako samo naponska ograničenja nisu zadovoljena, onda je jeftinije ugraditi izvore reaktivne snage nego prenosne elemente. Postoji jaka sprega između napona i reaktivnih snaga.

Po pravilu se koriste statičke metode (za jednu godinu). Pošto se proveravaju naponska ograničenja, onda se više ne može koristiti DC model tokova snaga, već mora AC model.

Kod nas se koristi najviše CLF-OPF (Constrained Load Flow – Optimum Power Flow). Bazira se na optimalnim tokovima snaga.

Kod planiranja izvora reaktivne snage period planiranja je puno kraći (do 10 godina). Vrlo su velike neizvesnosti u pogledu reaktivnih snaga.

Izvori reaktivne snage mogu biti induktivni i kapacitivni, čak i u istom EES-u (za režime maksimalnog (kapacitivni) i minimalnog (induktivni) opterećenja).

Ulazni podaci za ovu studiju su:

  1. Podaci o prognoziranim potrošnjama.
  2. Podaci o proizvodnim kapacitetima.
  3. Podaci o prenosnoj mreži.
  1. Prognoza mora biti za aktivnu i reaktivnu snagu. Najčešće se prognoza reaktivne snage direktno vezuje za aktivnu snagu preko konstantnog faktora snage. Ovo se mora uraditi za minimalan i maksimalan režim u svakoj godini planiranja.
  2. Od tehničkih karakteristika proizvodnih agregata treba znati minimalne i maksimalne proizvodnje aktivne snage, redosled angažovanja (funkciju troškova) i što je najvažnije karakteristike proizvodnje. One se daju preko pogonskih dijagrama sinhrone mašine. Ovo je detaljno rađeno prilikom učenja sinhronih generatora.
  3. Od tehničkih karakteristika prenosne mreže treba znati sve elemente ekvivalentnih šema, topološka struktura, zatim eksploataciona ograničenja: prenosne snage vodova i transformatora, minimalni i maksimalni naponi, regulacione mogućnosti transformatora.

Vrlo važno je da se koristi AC model tokova snaga za proveru naponskih ograničenja.

Za dogradnju novih izvora postoji više različitih metoda. Ovde će se izložiti metod potencijala čvorova.

Polazi se od linearnog modela injektiranja, koji je osnov za model tokova snaga:

I‾=Y‾BUSV‾,\ph{\boldsymbol{I}} = \ph{\boldsymbol{Y}}_{BUS}\ph{\boldsymbol{V}},

gde su elementi matrice admitanse nezavisnih čvorova Y‾BUS\ph{\boldsymbol{Y}}_{BUS}:

Y‾ij={Y‾i0+∑j∈αiY‾ijg,i=j−Y‾ijg,i≠j.\boxed{\ph{Y}_{ij} = \begin{cases} \ph{Y}_i^0 + \sum_{j \in \alpha_i}\ph{Y}_{ij}^g, & i = j \\ -\ph{Y}_{ij}^g, & i \ne j \end{cases}}.

Redukovana forma prethodne jednačine (bez referentno-balansnog čvora) je:

I‾r=Y‾rV‾r,\ph{\boldsymbol{I}}_{\boldsymbol{r}} = \ph{\boldsymbol{Y}}_{\boldsymbol{r}}\ph{\boldsymbol{V}}_{\boldsymbol{r}},

čije formalno rešenje je:

V‾r=(Y‾r)−1I‾r.\ph{\boldsymbol{V}}_{\boldsymbol{r}} = \left(\ph{\boldsymbol{Y}}_{\boldsymbol{r}}\right)^{-1}\ph{\boldsymbol{I}}_{\boldsymbol{r}}.

Ako se u i-tom čvoru dogradi izvor reaktivne snage, to se vidi u članu Y‾i0\ph{Y}_i^0, tako da se menja matrica Y‾BUS\ph{\boldsymbol{Y}}_{BUS}. Dakle, menjaju se samo dijagonalni elementi, tako da se ne mora invertovati matrica kompletno, već se može primeniti formula:

(Y‾r′)−1=(Y‾r)−1−1C‾Y‾i(1)AiT(Y‾r)−1;\left(\ph{\boldsymbol{Y}}'_{\boldsymbol{r}}\right)^{-1} = \left(\ph{\boldsymbol{Y}}_{\boldsymbol{r}}\right)^{-1} - \frac{1}{\ph{\boldsymbol{C}}}\ph{\boldsymbol{Y}}_i^{(1)}\boldsymbol{A}_i^{\mathrm{T}}\left(\ph{\boldsymbol{Y}}_{\boldsymbol{r}}\right)^{-1};

gde su:

  • C‾=Z‾i+(Y‾ii)−1\ph{\boldsymbol{C}} = \ph{Z}_i + \left(\ph{Y}_{ii}\right)^{-1};
  • Y‾i(1)\ph{\boldsymbol{Y}}_i^{(1)} – (N−1)(N - 1)-dimenziona i-ta kolona inverzne polazne (redukovane) matrice admitansi (Y‾r)−1\left(\ph{\boldsymbol{Y}}_{\boldsymbol{r}}\right)^{-1};
  • Ai=[0⋯1↓ i-ta⋯0]T\boldsymbol{A}_i = \begin{bmatrix} 0 & \cdots & \overset{\downarrow\, i\text{-ta}}{1} & \cdots & 0 \end{bmatrix}^{\mathrm{T}};
  • (Y‾ii)−1\left(\ph{Y}_{ii}\right)^{-1} – i-ti dijagonalni element matrice (Y‾r)−1\left(\ph{\boldsymbol{Y}}_{\boldsymbol{r}}\right)^{-1}.

Ako se prethodna jednačina pomnoži sa vektorom strujnih injektiranja, dobija se:

V‾r′=V‾r−V‾iZ‾i+(Y‾ii)−1Y‾i(1).\ph{\boldsymbol{V}}'_{\boldsymbol{r}} = \ph{\boldsymbol{V}}_{\boldsymbol{r}} - \frac{\ph{V}_i}{\ph{Z}_i + \left(\ph{Y}_{ii}\right)^{-1}}\ph{\boldsymbol{Y}}_i^{(1)}.

Iz prethodnog izraza dobija se fazor napona u i-tom čvoru posle dodavanja otočne impedanse.

V‾i′=V‾i−V‾i(Y‾ii)−1Z‾i+(Y‾ii)−1=Z‾iZ‾i+(Y‾ii)−1V‾i.\ph{V}'_i = \ph{V}_i - \frac{\ph{V}_i\left(\ph{Y}_{ii}\right)^{-1}}{\ph{Z}_i + \left(\ph{Y}_{ii}\right)^{-1}} = \frac{\ph{Z}_i}{\ph{Z}_i + \left(\ph{Y}_{ii}\right)^{-1}}\ph{V}_i.

Dakle, promena napona je direktno srazmerna priključenoj impedansi.

Priključena impedansa je Z‾i=jXi\ph{Z}_i = \mathrm{j}X_i, dok je redukovana matrica admitansi Y‾r\ph{\boldsymbol{Y}}_{\boldsymbol{r}} približno imaginarna (Y‾ii)−1=−jXii\left(\ph{Y}_{ii}\right)^{-1} = -\mathrm{j}X_{ii}, pa je čisto realan izraz:

Vi′=−XiXii−XiVi.V'_i = \frac{-X_i}{X_{ii} - X_i}V_i.

Iz prethodnog izraza mogu se izvesti uslovi za XiX_i da bi se postigli minimalni/maksimalni naponi u i-tom čvoru:

Xi=Vimin,maxVimin,max−ViXii.\boxed{X_i = \frac{V_i^{min,max}}{V_i^{min,max} - V_i}X_{ii}}.

Kao i ranije, može se definisati indeks efektivnosti dogradnje izvora reaktivne snage kao:

Fi=∑j∈αi∣ΔVj∣Ci(Xi);i∈γ,F_i = \frac{\sum\limits_{j \in \alpha_i}\left|\Delta V_j\right|}{C_i(X_i)}; \quad i \in \gamma,

gde su:

  • ΔVj\Delta V_j – promena napona u čvoru sa narušenim ograničenjem;
  • CiC_i – investicioni troškovi u i-tom čvoru;
  • αi\alpha_i – skup indeksa čvorova gde postoje naponska prekoračenja;
  • γ\gamma – skup čvorova gde se postavljaju izvori reaktivne snage.

Indeks efektivnosti predstavlja odnos poboljšanja naponskih prilika sa starim (“lošim”) naponima. Međutim, otežan je troškovima te priključene reaktanse.

Kriterijum izbora je:

max⁡i∈γ{Fi}.\max_{i \in \gamma}\left\{F_i\right\}.

Drugi (složeniji) pristup je primenom metoda optimalnih tokova snaga.

Teme

Planiranje · Tokovi snaga · Kompenzacija reaktivne snage

Zadaci uz ovo predavanje