Skip to content

Planiranje EES, TF Čačak i FTN: pitanja za usmeni sa odgovorima

Odgovori iz predavanja prof. dr Andrije T. Sarića

Svih 38 pitanja sa spiska za usmeni ispit iz Planiranja elektroenergetskih sistema, prvi deo (pitanja 1–19) i drugi deo (pitanja 20–38). Svaki odgovor je deo njegovog predavanja, skraćen na ono što pitanje traži, a ispod odgovora je link na ceo odeljak.

Prvi deo

1. Pojam, vrste, vremenska dekompozicija i metode prognoze potrošnje električne energije i snage.

Prognoza potrošnje je aktivnost predviđanja razvoja potreba potrošača električne energije u budućnosti.

Koriste se različiti matematički postupci. Osnovni cilj je objektivna procena. Prognoza je prvi korak u planiranju EES-a. Zbog njenog značaja po pravilu je rade iskusni inženjeri.

Postoje tri tipa osnovnih karakteristika za prognozu:

  • Prognoza potreba sistema u energiji.
  • Prognoza maksimalnih (vršnih) i minimalnih opterećenja.
  • Prognoza kompletne (normalizovane) krive trajanja opterećenja.

Po vremenskom horizontu podela je:

  • Kratkoročna prognoza, 1 do 7 dana i satna diskretizacija. Bitna je za dnevno i nedeljno planiranje rada EES-a (program rada elektrana, planiranje razmene električne energije, planiranje operativne rezerve, proračun pouzdanosti i drugi).
  • Srednjeročna prognoza, 1 mesec do 5 godina i nedeljna diskretizacija. Vezana za srednjeročno operativno planiranje (mesečni i godišnji operativni bilansi, izrada plana remonata, punjenje i pražnjenje akumulacija hidroelektrana, razmena energije, nabavka goriva). Ovo je spona planiranja i eksploatacije EES-a.
  • Dugoročna prognoza, 5 do 30 godina i godišnja diskretizacija. Globalno planiranje investicija, dogradnja i rekonstrukcija proizvodnih i prenosnih kapaciteta.

Podela po metodu prognoziranja i prisustvu neelektričnih veličina:

  • Kvalitativni metodi.
  • Ekstrapolacioni (nezavisni) metodi.
  • Korelacioni (zavisni) metodi.
  • Kombinovani metodi.

Podela po prisustvu neizvesnosti:

  • Deterministički metodi (egzaktni ulazni podaci).
  • Stohastički metodi (neizvesni ulazni podaci).

Iz predavanja: Prognoza potrošnje električne energije i snage; Prognoza potrošnje električne energije i snage: Podele prognoze potrošnje.

2. Ekstrapolacioni (nezavisni) metod prognoze potrošnje.

Oslanjaju se na ekstrapolaciju podataka iz prošlosti.

Mogu biti deterministički i stohastički.

Osnovna ideja je analiza vremenskog razvoja: dolazi se do vremenske zavisnosti potrošnje od vremena, koja se zatim ekstrapolira.

Prvo se crta disperzioni dijagram iz prošlosti, što je ilustrovano na slici:

Disperzioni dijagram: tačke X_t^0, ostvarenja iz prošlosti, kroz koje prolazi kriva X_t = f(t), a dalje se tačkasto nastavlja kao prognoza za budućnost

Generalno, prognoza je pouzdanija ako je period iz prošlosti puno duži od perioda prognoze.

Dakle, osnovni cilj je pronaći najbolju funkciju Xt=f(t)X_t = f(t).

Postoji veći broj analitičkih funkcija za analizu, od kojih su najvažnije:

  • Polinom: Xt=∑iβitiX_t = \sum\limits_i\beta_it^i
  • Parabola: Xt=atbX_t = at^b
  • Logaritamska prava: log⁡Xt=a+bt\log X_t = a + bt
  • Logaritamska parabola: log⁡Xt=a+bt+ct2\log X_t = a + bt + ct^2
  • Logistička funkcija: log⁡Xt=a1+be−ct\log X_t = \dfrac{a}{1 + be^{-ct}}

Za određivanje koeficijenata optimalna je primena kriterijuma minimuma odstupanja sume kvadrata:

Osnovna pretpostavka je da će se razvoj u budućnosti nastaviti istim trendom (istom godišnjom stopom porasta).

Ako je taj trend konstantan: pt=p=Const.p_t = p = \text{Const.}, onda je (najčešće u praksi):

Xt=X0(1+p)t.X_t = X_0(1 + p)^t.

Ako se prethodni izraz logaritmuje:

log⁡Xt=log⁡X0+tlog⁡(1+p)=a+bt,\log X_t = \log X_0 + t\log(1 + p) = a + bt,

gde je a=log⁡X0a = \log X_0 i b=log⁡(1+p)b = \log(1 + p).

Stopa porasta u dve susedne godine je:

1+p=XtXt−1.1 + p = \frac{X_t}{X_{t-1}}.

Iz predavanja: Prognoza potrošnje električne energije i snage: Ekstrapolacioni (nezavisni) metodi.

3. Izvesti potrebne uslove za određivanje parametara prognostičkog modela pomoću logističke funkcije.

Postoji veći broj analitičkih funkcija za analizu, od kojih su najvažnije:

  • Polinom: Xt=∑iβitiX_t = \sum\limits_i\beta_it^i
  • Parabola: Xt=atbX_t = at^b
  • Logaritamska prava: log⁡Xt=a+bt\log X_t = a + bt
  • Logaritamska parabola: log⁡Xt=a+bt+ct2\log X_t = a + bt + ct^2
  • Logistička funkcija: log⁡Xt=a1+be−ct\log X_t = \dfrac{a}{1 + be^{-ct}}

Za određivanje koeficijenata optimalna je primena kriterijuma minimuma odstupanja sume kvadrata:

Na primer:

εt(βi, a, b, c)=∑t=1N(Xt−Xt0)2,\varepsilon_t(\beta_i,\, a,\, b,\, c) = \sum_{t=1}^{N}\left(X_t - X_t^0\right)^2,

ili:

εt(βi, a, b, c)=∑t=1N(log⁡Xt−log⁡Xt0)2.\varepsilon_t(\beta_i,\, a,\, b,\, c) = \sum_{t=1}^{N}\left(\log X_t - \log X_t^0\right)^2.

Ovaj minimum ima se za:

∂εt(βi, a, b, c)∂βi=0;∂εt(βi, a, b, c)∂a=0;∂εt(βi, a, b, c)∂b=0;∂εt(βi, a, b, c)∂c=0.\begin{gathered} \frac{\partial\varepsilon_t(\beta_i,\, a,\, b,\, c)}{\partial\beta_i} = 0; \\[0.5em] \frac{\partial\varepsilon_t(\beta_i,\, a,\, b,\, c)}{\partial a} = 0; \\[0.5em] \frac{\partial\varepsilon_t(\beta_i,\, a,\, b,\, c)}{\partial b} = 0; \\[0.5em] \frac{\partial\varepsilon_t(\beta_i,\, a,\, b,\, c)}{\partial c} = 0. \end{gathered}

Ovako se dobija onoliko jednačina koliko ima nepoznatih parametara: βi\beta_i, aa, bb i cc.

Iz predavanja: Prognoza potrošnje električne energije i snage: Ekstrapolacioni (nezavisni) metodi.

4. Zavisni (korelacioni) metod prognoze potrošnje; način određivanja parametara modela.

Potrošnja energije zavisi od niza uticajnih faktora: veličine i strukture društvenog proizvoda, cene drugih vrsta energije (vrlo važan uticaj), nivoa elektrifikacije domaćinstava i industrije (upotreba gasa i struktura energetskih resursa) i drugih.

Cilj je da se ove vrednosti nekako kvantifikuju.

Zato se u proračun uvode egzogene (spoljašnje) i endogene (unutrašnje) promenljive.

Najčešće egzogene promenljive za opisivanje razvoja potrošnje su:

  • Bruto nacionalni proizvod (indeks društveno-ekonomskih aktivnosti).
  • Priraštaj stanovništva.
  • Cene i prodaje električne energije.
  • Razvoj industrije (poljoprivreda, građevinarstvo).

Za opis sezonskih varijacija potrošnje koriste se sledeće egzogene promenljive:

  • Temperatura.
  • Vlažnost.
  • Osvetljaj (oblačnost).
  • Brzina i pravac vetra.

Dakle, u opštem slučaju koriste se nelinearne relacije – višestruka nelinearna regresija:

Xt=F(Egzogene promenljive)+Slucˇajna promenljiva.X_t = F(\text{Egzogene promenljive}) + \text{Slučajna promenljiva}.

Detaljnija matematička forma je:

Xt=β0+∑i=1LβiYti+ηt.X_t = \beta_0 + \sum_{i=1}^{L}\beta_iY_{ti} + \eta_t.

Kao i kod nezavisnih modela prognoze potrošnje, ovi koeficijenti se određuju iz minimuma sume kvadrata odstupanja:

min⁡{εt(β0, β1, ⋯ βL)}=min⁡{∑t=1N(Xt−Xt0)2};\min\{\varepsilon_t(\beta_0,\, \beta_1,\, \cdots\, \beta_L)\} = \min\left\{\sum_{t=1}^{N}\left(X_t - X_t^0\right)^2\right\};
∂εt(β0, β1, ⋯ βL)∂β0=0;  …  ; ∂εt(β0, β1, ⋯ βL)∂βL=0.\frac{\partial\varepsilon_t(\beta_0,\, \beta_1,\, \cdots\, \beta_L)}{\partial\beta_0} = 0;\ \ \ldots\ \ ;\ \frac{\partial\varepsilon_t(\beta_0,\, \beta_1,\, \cdots\, \beta_L)}{\partial\beta_L} = 0.

Iz predavanja: Prognoza potrošnje električne energije i snage: Zavisni (korelacioni) modeli.

5. Zavisni (korelacioni) metod prognoze potrošnje u funkciji bruto nacionalnog proizvoda.

Poznato je da bruto nacionalni proizvod (BNP) utiče na potrošnju električne energije.

Pretpostavlja se linearni korelacioni model:

Xt=a+bYt;Xt=Wptg;Yt=(BNP)t,X_t = a + bY_t; \quad X_t = W_{pt}^g; \quad Y_t = (BNP)_t,

gde su:

  • XtX_t – ukupna godišnja proizvodnja električne energije (endogena promenljiva);
  • YtY_t – bruto nacionalni proizvod (egzogena promenljiva).

Promena BNP je:

Yt=Yt−1(1+p(BNP)t),Y_t = Y_{t-1}(1 + p_{(BNP)t}),

pa je:

XtXt−1=WptgWp(t−1)g=1+pt=a+bYta+bYt−1=a+bYt−1(1+p(BNP)t)a+bYt−1.\frac{X_t}{X_{t-1}} = \frac{W_{pt}^g}{W_{p(t-1)}^g} = 1 + p_t = \frac{a + bY_t}{a + bY_{t-1}} = \frac{a + bY_{t-1}(1 + p_{(BNP)t})}{a + bY_{t-1}}.

Onda je:

pt=p(BNP)t(1−aa+bYt−1).p_t = p_{(BNP)t}\left(1 - \frac{a}{a + bY_{t-1}}\right).

Ako se uvedu srednje vrednosti:

p(BNP)t=pBNP=Const.;Yt−1=Y0(1+pBNP)t−1.\begin{gathered} p_{(BNP)t} = p_{BNP} = \text{Const.}; \\ Y_{t-1} = Y_0(1 + p_{BNP})^{t-1}. \end{gathered}

Tada je stopa porasta potrošnje električne energije:

pt=pBNP(1−aa+bY0(1+pBNP)t−1).\boxed{p_t = p_{BNP}\left(1 - \frac{a}{a + bY_0(1 + p_{BNP})^{t-1}}\right)}.

Iz predavanja: Prognoza potrošnje električne energije i snage: Zavisni (korelacioni) model u funkciji bruto nacionalnog proizvoda.

6. Prognoza vršne (maksimalne) snage; primena ekstrapolacione metode.

Osnovna dilema je da li prognozu treba sprovesti na bazi prognoze potrošnje, ili nezavisno.

Ako se povezuju, onda se vezuju preko faktora opterećenja.

Prednost nezavisnog planiranja je što je to direktna metoda, ali postoji problem uticajnih faktora.

Neophodni podaci su: godišnje potrošnje, faktori opterećenja, odnos min/max opterećenja i drugi.

Veze između ovih veličina su obrađene ranije:

PMg=Wg8760⋅mg;Pmg=m0gPMg,\begin{gathered} P_M^g = \frac{W^g}{8760 \cdot m^g}; \\[0.5em] P_m^g = m_0^gP_M^g, \end{gathered}

gde su:

  • mgm^g – faktor godišnjeg opterećenja;
  • m0gm_0^g – odnos godišnjeg minimalnog i maksimalnog opterećenja.

Ako se vrši srednjeročna prognoza, vrlo je veliki uticaj sezonske komponente, tako da generalno postoji linija trenda na koju se superponira sezonska komponenta. Ovaj efekat je prikazan na slici:

Levo: snaga P_p kroz godine, talasasta kriva koja se ponavlja svake godine i raste. Desno: ista kriva razložena na liniju trenda i sezonsku komponentu

Iz predavanja: Prognoza potrošnje električne energije i snage: Prognoza vršne (maksimalne) snage; Prognoza potrošnje električne energije i snage: Ekstrapolaciona metoda prognoze vršne snage.

7. Troškovi u hidro i termoelektranama.

Troškovi u hidroelektranama

Postoje investicioni (CIC_I) i eksploatacioni, odnosno pogonski troškovi (CExC_{Ex}).

Godišnji investicioni troškovi su proporcionalni instalisanoj snazi (sa koeficijentom csc_s):

CIg=cIgCI=cIgcsPEi,C_I^g = c_I^gC_I = c_I^gc_sP_{Ei},

gde su:

  • cIgc_I^g – godišnja stopa investicionih troškova;
  • csc_s – koeficijent proporcionalnosti troškova sa instalisanom snagom;
  • PEiP_{Ei} – instalisana snaga.

Eksploatacioni troškovi sastoje se od stalnih i promenljivih troškova pogona i održavanja:

CEg=CIg+CExg=(cIg+cOMog)CI+cOMvWEg=CFg+Cvarg,C_E^g = C_I^g + C_{Ex}^g = \left(c_I^g + c_{OM}^{og}\right)C_I + c_{OM}^vW_E^g = C_F^g + C_{var}^g,

Troškovi u termoelektranama

Investicioni troškovi su proporcionalni instalisanoj snazi:

CIg=cIgCI=cIgcsPEi,C_I^g = c_I^gC_I = c_I^gc_sP_{Ei},

što je isto kao u slučaju hidroelektrana.

Eksploatacioni troškovi sastoje se od:

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

Promenljivi troškovi goriva predstavljaju najvažnije troškove. Na godišnjem nivou su reda veličine investicionih troškova. Da bi se izračunali mora se izračunati ukupna potrošnja toplote:

DT=∫0TD(Pg) dt;Pg=Pg(t).D_T = \int_0^T D(P_g)\,\mathrm{d}t; \quad P_g = P_g(t).

Ovaj integral nije moguće izračunati za mnogo godina unapred, jer se ne poznaje dijagram opterećenja. Zato se koriste linearne aproksimacije, koje su date ranije.

Iz predavanja: Energetsko-eksploatacione karakteristike hidroelektrana: Troškovi u hidroelektranama; Energetsko-eksploatacione karakteristike termoelektrana: Troškovi u termoelektranama.

8. Osnovni faktori svođenja kapitala.

Ako se uvedu sledeće oznake:

  • PP – sadašnja (present) ekvivalentna vrednost novčanog toka (prihodi, rashodi);
  • FF – buduća (future) ekvivalentna vrednost novčanog toka;
  • RR – uniformni niz novčanog toka (na primer, otplata duga u jednakim vremenskim intervalima), gde je RgR^g godišnja rata (anuitet);
  • nn – broj intervala (godina) posmatranja;
  • ii – stopa aktualizacije (može, a ne mora biti jednaka kamatnoj stopi);
  • aa – stopa inflacije (neće se dalje uvoditi, ali ona povećava stopu aktualizacije);
  • (x/y, i, n)(x/y,\, i,\, n) – faktor svođenja veličine yy na veličinu xx pri aktualizacionoj (kamatnoj) stopi ii i jednakim vremenskim intervalima.

Prost kamatni faktor (Simple Interest Factor – SIF)

Kod njega se kamata računa samo na glavnicu:

F=P+niP=P⋅(SIF);SIF=1+ni.\boxed{F = P + niP = P \cdot (SIF)}; \quad \boxed{SIF = 1 + ni}.

Složeni kamatni faktor (Compound Interest Factor – CIF)

Kod njega se računa složena kamata, odnosno kamata na kamatu:

F=P(1+i)n=P⋅(CIF);CIF=(1+i)n.\boxed{F = P(1 + i)^n = P \cdot (CIF)}; \quad \boxed{CIF = (1 + i)^n}.

Faktor svođenja na sadašnju vrednost (Present Value Factor – PVF)

Recipročna vrednost složenog kamatnog faktora (CIF). Svodi buduću vrednost kapitala na sadašnju vrednost:

P=F(1+i)−n=F⋅(PVF)=F(CIF);PVF=(1+i)−n.\boxed{P = F(1 + i)^{-n} = F \cdot (PVF) = \frac{F}{(CIF)}}; \quad \boxed{PVF = (1 + i)^{-n}}.

Kamatni faktori u uniformnom nizu

1. Faktor svođenja uniformnog niza na sadašnju/buduću vrednost (Present Worth Factor – PWF)

Specifikacija problema: Postoji niz godišnjih jednakih uplata ili isplata. Cilj je da se utvrdi njihova sadašnja vrednost. Vrlo aktuelno za kreditiranje (treba malo bolje objasniti na realnom primeru kredita):

P=Rg1+i+Rg(1+i)2+⋯+Rg(1+i)n;pomnozˇiti sa /(1+i).\boxed{P = \frac{R^g}{1 + i} + \frac{R^g}{(1 + i)^2} + \cdots + \frac{R^g}{(1 + i)^n}}; \quad \text{pomnožiti sa } /(1 + i).

To znači da je uniformna rata otplate:

Rg=Pi(1+i)n(1+i)n−1,R^g = P\frac{i(1 + i)^n}{(1 + i)^n - 1},

koja određuje faktor povraćaja kapitala (Capital Recovery Factor – CRF)

(CRF)=(Rg/P, i, n)=i(1+i)n(1+i)n−1.(CRF) = (R^g/P,\, i,\, n) = \frac{i(1 + i)^n}{(1 + i)^n - 1}.

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Dijagrami novčanih tokova; Planiranje elektroenergetskih sistema i inženjerska ekonomija: Prost kamatni faktor (Simple Interest Factor – SIF); Planiranje elektroenergetskih sistema i inženjerska ekonomija: Složeni kamatni faktor (Compound Interest Factor – CIF); Planiranje elektroenergetskih sistema i inženjerska ekonomija: Faktor svođenja na sadašnju vrednost (Present Value Factor – PVF); Planiranje elektroenergetskih sistema i inženjerska ekonomija: Kamatni faktori u uniformnom nizu.

9. Amortizacija i ostatak vrednosti.

Prilikom kupovine treba izdvojiti sredstva u visini nabavne vrednosti. Nakon puštanja u pogon, dolazi do istrošenosti, tako da je na kraju životnog veka on u nemogućnosti da obavlja svoju funkciju.

Amortizacija je knjigovodstvena procedura za smanjenje vrednosti nekog kapitalnog dobra sa vremenom, usled starenja i istrošenosti tokom ekonomskog veka. U principu, amortizacioni vek je različit od životnog.

Postoji više šema za proračun amortizacije:

Pravolinijska amortizacija (UNIF), koja podrazumeva da se ona ravnomerno raspodeljuje tokom životnog veka (jednostavna za primenu).

Ubrzana (usporena) amortizacija (SYD i ACRS šeme), koja podrazumeva veću (manju) stopu amortizacije u prvim godinama, a kasnije manju (veću). Jednostavan primer za ovaj efekat je amortizacija automobila (pa i većine drugih dobara), koji se znatno više amortizuju u početnim godinama.

Pravolinijska šema amortizacije (UNIF) opisuje se jednačinom:

[pag(UNIF)]k=1ℓ;k=1, 2, ⋯ , ℓ,\boxed{[p_a^g(UNIF)]_k = \frac{1}{\ell}}; \quad k = 1,\, 2,\, \cdots,\, \ell,

gde je ℓ\ell očekivani ekonomski vek objekta (može biti jednak životnom, ali ne mora).

Iz prethodnog izraza sledi da je:

∑k=1ℓ[pag(UNIF)]k=1.\boxed{\sum_{k=1}^{\ell}[p_a^g(UNIF)]_k = 1}.

Ako je vrednost objekta na početku V0V_0, a na kraju ekonomskog (životnog) veka VℓostV_\ell^{ost} (ostatak vrednosti, “staro gvožđe”), godišnji iznos amortizacije je:

[Qg(UNIF)]k=V0−Vℓostℓ[NJ/god].[Q^g(UNIF)]_k = \frac{V_0 - V_\ell^{ost}}{\ell} \quad [\text{NJ/god}].

Ukupna akumulisana amortizacija do k-te godine je:

Vkakum=kV0−Vℓostℓ,V_k^{akum} = k\frac{V_0 - V_\ell^{ost}}{\ell},

dok je preostala (ostatna) vrednost

Vkost=Vℓost+ℓ−kℓ(V0−Vℓost),V_k^{ost} = V_\ell^{ost} + \frac{\ell - k}{\ell}(V_0 - V_\ell^{ost}),

Ubrzana šema amortizacije (SYD) opisuje se jednačinom:

[pag(SYD)]k=2ℓ−k+1ℓ(ℓ+1);k=1, 2, ⋯ , ℓ,\boxed{[p_a^g(SYD)]_k = 2\frac{\ell - k + 1}{\ell(\ell + 1)}}; \quad k = 1,\, 2,\, \cdots,\, \ell,

gde je karakteristični koeficijent dobijen iz 1+2+⋯+ℓ=ℓ(ℓ+1)21 + 2 + \cdots + \ell = \dfrac{\ell(\ell + 1)}{2}.

Godišnji iznos amortizacije je:

[Qg(SYD)]k=2ℓ−k+1ℓ(ℓ+1)(V0−Vℓost)[NJ/god].[Q^g(SYD)]_k = 2\frac{\ell - k + 1}{\ell(\ell + 1)}\left(V_0 - V_\ell^{ost}\right) \quad [\text{NJ/god}].

Usporena šema amortizacije sprovodi se tako što se osniva poseban amortizacioni fond za prikupljanje sredstava. Ova sredstva se koriste i donose dobit po godišnjoj stopi aktualizacije i. Ovde se neće detaljnije izlagati (videti knjigu prof. Levi).

Napred opisane šeme amortizacije prikazane su na slici:

Vrednost objekta V u funkciji godine k od V_0 do V_ℓ^ost: usporena šema amortizacije iznad prave, jednosmerna šema je prava, a ubrzana šema ispod prave, sa ostatkom V_k^ost na svakoj

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Amortizacija i ostatak vrednosti.

10. Metodi ekvivalentnih vrednosti za ekonomsku ocenu investicija.

Metod sadašnje ekvivalentne vrednosti

Suština je pretvaranje svih novčanih tokova u ekvivalentnu sadašnju vrednost.

Za ovo svođenje koristi se stopa aktualizacije i:

  • Ako se radi o privatnim kompanijama i je jednaka minimalno atraktivnoj stopi povraćaja kapitala (MARR – Minimum Attractive Rate of Return). Ovu stopu određuju kompanije u cilju optimizacije svog poslovanja.
  • Za državna preduzeća ova se stopa propisuje od strane države. Obično je jednaka aktualizacionoj (diskontnoj) stopi, koja je obično 8–10 %.

Sadašnja ekvivalentna vrednost (PV) je:

PV=F1g(1+i)−1+⋯+Fng(1+i)−n=∑k=1nFkg(1+i)−k,PV = F_1^g(1 + i)^{-1} + \cdots + F_n^g(1 + i)^{-n} = \sum_{k=1}^{n}F_k^g(1 + i)^{-k},

gde je FkgF_k^g budući novčani tok u godini k=1, 2, ⋯ , nk = 1,\, 2,\, \cdots,\, n:

Fkg=Bkg−(CIkg+CExkg),F_k^g = B_k^g - \left(C_{Ik}^g + C_{Exk}^g\right),

gde su:

  • BkgB_k^g – godišnji prihod;
  • CIkgC_{Ik}^g – godišnji investicioni troškovi;
  • CExkgC_{Exk}^g – godišnji eksploatacioni troškovi.

Naravno, ekonomska dobit je veća što je PV veće.

Ako ima više investicionih šema, kriterijum izbora je:

max⁡j{(PV)j}=max⁡j{∑k=1nFkg(1+i)−k}j;j=1, 2, ⋯ , m.\max_j\left\{(PV)_j\right\} = \max_j\left\{\sum_{k=1}^{n}F_k^g(1 + i)^{-k}\right\}_j; \quad j = 1,\, 2,\, \cdots,\, m.

Svaka šema sa (PV)j>0(PV)_j > 0 je ostvariva (dobra).

Metod godišnjih ekvivalentnih vrednosti

Bazira se na proračunu kumulativnih godišnjih ekvivalentnih vrednosti (AV):

(AV)=(AVB)−(AVC)−[AV(CR)],(AV) = (AVB) - (AVC) - [AV(CR)],
  • BB – prihodi (benefit);
  • CC – troškovi (cost);
  • CRCR – povraćaj kapitala (Capital Recovery);
  • CRFCRF – faktor povraćaja kapitala (definisan ranije);
  • PVFPVF – faktor svođenja budućih novčanih tokova.

Projekat je ekonomski isplativ ako je (AV)j>0(AV)_j > 0.

Kriterijum izbora je max⁡j{(AV)j}\max\limits_j\left\{(AV)_j\right\}.

Na sličan način može se definisati i metod buduće ekvivalentne vrednosti.

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metod sadašnje ekvivalentne vrednosti; Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metod godišnjih ekvivalentnih vrednosti.

11. Metodi najkraćeg perioda otplate projekta za ekonomsku ocenu investicija.

Nekada se problem izbora može posmatrati i sa gledišta brzine realizacije i povraćaja uloženog kapitala. Cilj je što pre povratiti uloženi kapital (Payback Period). Dakle, umesto isplativosti, kao osnovni kriterijum se posmatra vreme povraćaja kapitala.

Period povraćaja investicije definiše se kao broj godina za koji je ušteda na eksploatacionim troškovima jednaka investicionim ulaganjima. Ušteda u eksploatacionim troškovima je konstantna sa vremenom.

Projekat je bolje rangiran što je vreme povraćaja kraće.

Metod je pogodan za inicijalno filtriranje i eliminaciju nekih opcija, i izbor užih kandidata.

Matematički izraženo, za broj godina θ\theta je:

∑k=1θ(Bkag−Ckag)>CI, odnosno ∑k=1θ(Bkg⋅(PVF)k−Ckg⋅(PVF)k)>CI\sum_{k=1}^{\theta}\left(B_{ka}^g - C_{ka}^g\right) > C_I\text{, odnosno } \sum_{k=1}^{\theta}\left(B_k^g \cdot (PVF)_k - C_k^g \cdot (PVF)_k\right) > C_I

gde su CIC_I ukupne investicije.

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metodi najkraćeg perioda otplate projekta.

12. Metodi stope povraćaja i poređenja koristi i troškova za ekonomsku ocenu investicija.

Kod metode ekvivalentnih vrednosti pretpostavlja se konstantna kamatna stopa, odnosno stopa aktualizacije i (na primer, jednaka MARR).

Kod ove metode mera ekonomičnosti je kamatna stopa, koja je nepoznata. Koriste se metode:

  • Interne stope povraćaja (IRR – Internal Rate of Return)
  • Eksterne stope povraćaja (ERR – External Rate of Return).

Po pravilu, primenjuje se kod ulaganja u preduzeća konkurentske privrede (aktuelno kod deregulacije).

Metod interne stope povraćaja (metod profitnog indeksa)

Suština je da se traži stopa aktualizacije i koja izjednačava sadašnje ekvivalentne vrednosti prihoda i rashoda u nekom vremenskom periodu.

Praktično se ispituju novčani tokovi za nekoliko stopa aktualizacije. Ona stopa koja rezultuje u buduće novčane tokove koji izjednačuju prihode i rashode, naziva se interna stopa povraćaja (IRR). IRR je pozitivno kada priliv novca nadmašuje izdatke.

Najbolje rešenje između više opcija je ono sa najvećim IRR.

Matematički se ovaj metod specificira kao:

∑k=1nBkg(P/Fkg, i′, k)=∑k=1nCkg(P/Fkg, i′, k),\sum_{k=1}^{n}B_k^g\left(P/F_k^g,\, i',\, k\right) = \sum_{k=1}^{n}C_k^g\left(P/F_k^g,\, i',\, k\right),

gde je BkgB_k^g godišnji priliv u k-toj godini, a CkgC_k^g godišnji trošak u k-toj godini.

Metod poređenja koristi i troškova

Suština je u poređenju koristi (dobiti) i troškova (Benefit to Cost (B-C) Method).

Cilj je da se proračuna odnos B/C = Korist/Troškovi, i to svođenjem novčanih tokova na sadašnju ekvivalentnu vrednost (PV):

(BC)=(PVBPVC)=Sadasˇnja ekvivalentna vrednost koristiSadasˇnja ekvivalentna vrednost trosˇkova\left(\frac{B}{C}\right) = \left(\frac{PVB}{PVC}\right) = \frac{\text{Sadašnja ekvivalentna vrednost koristi}}{\text{Sadašnja ekvivalentna vrednost troškova}}

ili na godišnju ekvivalentnu vrednost:

(BgCg)=(AVBAVC)=Godisˇnja ekvivalentna vrednost koristiGodisˇnja ekvivalentna vrednost trosˇkova\left(\frac{B^g}{C^g}\right) = \left(\frac{AVB}{AVC}\right) = \frac{\text{Godišnja ekvivalentna vrednost koristi}}{\text{Godišnja ekvivalentna vrednost troškova}}

Projekat je isplativ ako je (BC)≥1\left(\dfrac{B}{C}\right) \geq 1.

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metodi stope povraćaja; Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metod interne stope povraćaja (metod profitnog indeksa); Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metod poređenja koristi i troškova.

13. Metod minimalnih troškova za ekonomsku ocenu investicija.

Često je teško izračunati prihode elektroenergetskih objekata (na primer, vodovi i transformatori). U tom slučaju se polazi od pretpostavke da sve opcije imaju istu dobit, odnosno da se ona ne uzima u obzir. Znači, analiziraju se samo troškovi.

Znači, kriterijum je minimizacija sadašnjih ekvivalentnih troškova (PVC):

min⁡j{(PVC)j}=min⁡j{∑k=1nCkg(PVF)k}j;j=1, 2, ⋯ , m,\min_j\left\{(PVC)_j\right\} = \min_j\left\{\sum_{k=1}^{n}C_k^g(PVF)_k\right\}_j; \quad j = 1,\, 2,\, \cdots,\, m,

gde su Ckg=CIkg+CExkgC_k^g = C_{Ik}^g + C_{Exk}^g investicioni i eksploatacioni troškovi u k-toj godini.

Može se izračunati i preko izravnatih godišnjih ekvivalentnih troškova (L – Levelized):

min⁡j{(LAVC)j};j=1, 2, ⋯ , m.\min_j\left\{(LAVC)_j\right\}; \quad j = 1,\, 2,\, \cdots,\, m.

L-metod izravnava godišnje ekvivalentne vrednosti na jedinstvenu godinu. Ovaj pristup izravnavanja može se primeniti na većinu prethodnih metoda.

Međutim, kod prethodnog kriterijuma problemi nastaju ako opcije imaju različito trajanje. Tada se koristi svođenje na najmanji zajednički sadržilac, odnosno isti fiktivni životni vek objekta. Ovo je demonstrirano u zadacima.

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Metod minimalnih troškova.

14. Primeri svođenja investicionih troškova na sadašnju vrednost.

Analizira se izgradnja objekta u k-toj godini perioda planiranja od n godina, čiji su investicioni troškovi CIC_I. Životni vek objekta je ℓ>n−(k−1)=n−k+1\ell > n - (k - 1) = n - k + 1.

Analiziraju se dva karakteristična slučaja:

1. Postojanje sopstvenog kapitala

Vremenska osa godina od 1 do n: kapital C_I se ulaže odjednom na početku k-te godine (sopstveni kapital), a na kraju n-te godine ostaje ostatak vrednosti V_ℓ^ost

Sadašnja ekvivalentna vrednost je:

PV=CI(P/F, i, k−1)−Vℓost(P/F, i, n)PV=CI⋅1(1+i)k−1−Vℓost⋅1(1+i)n.\begin{gathered} PV = C_I(P/F,\, i,\, k - 1) - V_\ell^{ost}(P/F,\, i,\, n) \\[0.5em] PV = C_I \cdot \frac{1}{(1 + i)^{k-1}} - V_\ell^{ost} \cdot \frac{1}{(1 + i)^n}. \end{gathered}

VℓostV_\ell^{ost} se oduzima od investicija, pošto se umanjuje.

2. Korišćenje pozajmljenog kapitala (kredit)

Vremenska osa godina od 1 do n: kredit C_I na početku k-te godine i jednake godišnje rate R^g na kraju svake godine od k-te do n-te

Uniformni godišnji novčani tok (anuitet) je:

Rg=CI(Rg/P, i, n−k+1),R^g = C_I(R^g/P,\, i,\, n - k + 1),

što je određeno faktorom povraćaja kapitala (CRF)=i(1+i)n−k+1(1+i)n−k+1−1(CRF) = \dfrac{i(1 + i)^{n-k+1}}{(1 + i)^{n-k+1} - 1}.

Sadašnja ekvivalentna vrednost jednog anuiteta je (zanemarena ostatna vrednost VℓostV_\ell^{ost}):

PV(Rg)=Rg(P/Rg, i, n−k+1)⋅(P/F, i, k−1),PV(R^g) = R^g(P/R^g,\, i,\, n - k + 1) \cdot (P/F,\, i,\, k - 1),
(P/Rg, i, n−k+1)(P/R^g,\, i,\, n - k + 1)(P/F, i, k−1)(P/F,\, i,\, k - 1)
Postoji n−k+1n - k + 1 anuitetaPrvi anuitet počinje u k-toj godini
(PWF)=(1+i)n−k+1−1i(1+i)n−k+1(PWF) = \dfrac{(1 + i)^{n-k+1} - 1}{i(1 + i)^{n-k+1}}(PVF)=(1+i)−(k−1)(PVF) = (1 + i)^{-(k-1)}

Primer:

Postoji objekat bazne vrednosti 10000 NJ u sadašnjem trenutku, čija je ostatna vrednost na kraju životnog veka nula. Životni (amortizacioni) vek je 5 godina, dok je stopa aktualizacije 10 %. Primenjuje se pravolinijska šema amortizacije.

GodinaVrednost objekta na početku godineAmortizacija + Aktualizacija (kamata) i = 10 %Sadašnja vrednost povraćaja kapitala
1100002000+10003000⋅(P/F, 10 %, 1)=2727,33000 \cdot \left(P/F,\, 10\ \%,\, 1\right) = 2727{,}3
280002000+8002800⋅(P/F, 10 %, 2)=2313,92800 \cdot \left(P/F,\, 10\ \%,\, 2\right) = 2313{,}9
360002000+6002600⋅(P/F, 10 %, 3)=1953,42600 \cdot \left(P/F,\, 10\ \%,\, 3\right) = 1953{,}4
440002000+4002400⋅(P/F, 10 %, 4)=1639,22400 \cdot \left(P/F,\, 10\ \%,\, 4\right) = 1639{,}2
520002000+2002200⋅(P/F, 10 %, 5)=1365,92200 \cdot \left(P/F,\, 10\ \%,\, 5\right) = 1365{,}9
Σ=13000\Sigma = 13000Σ≈10000\Sigma \approx 10000 NJ
(PVF)=(1+i)−k(PVF) = (1 + i)^{-k}

Godišnja rata je Rg=10000⋅(Rg/P, 10 %, 5)=2638R^g = 10000 \cdot \left(R^g/P,\, 10\ \%,\, 5\right) = 2638 NJ, gde je

(CRF)=(Rg/P, i, n)=i(1+i)n(1+i)n−1.(CRF) = \left(R^g/P,\, i,\, n\right) = \frac{i(1 + i)^n}{(1 + i)^n - 1}.

Iz predavanja: Planiranje elektroenergetskih sistema i inženjerska ekonomija: Svođenje investicionih troškova na sadašnju vrednost; Planiranje elektroenergetskih sistema i inženjerska ekonomija: 1. Postojanje sopstvenog kapitala; Planiranje elektroenergetskih sistema i inženjerska ekonomija: 2. Korišćenje pozajmljenog kapitala (kredit).

15. Vetroelektrane (pojam, osnovne tehničke/konstrukcione karakteristike i proračun snage vetroturbine).

Vetroelektrane preko vetroturbina i električnih generatora pretvaraju kinetičku energiju vetra prvo u mehaničku, a zatim u električnu energiju. Vetroelektrane su danas najpopularniji tip novih i obnovljivih izvora, jer je vetar dostupan i obilan obnovljivi energent, imaju usavršenu tehnologiju, pa već sada igraju značajnu ulogu u proizvodnji električne energije. Proizvodnja vetrogeneratora hemijski ne zagađuje životnu sredinu, ne ostavlja rizičan otpad, niti zahteva velika infrastrukturna ulaganja.

Snaga vetroturbine izračunava se preko formule

PWT=PW0cp(λ,θ)=ARρv2cp(λ,θ)vW3,P_{WT} = P_W^0c_p(\lambda, \theta) = A_R\frac{\rho_v}{2}c_p(\lambda, \theta)v_W^3,

gde su:

  • ρv\rho_v – gustina vazduha (≈\approx 1,225 kg/m3^3);
  • ARA_R – površina zahvatanja lopatica rotora;
  • cp(λ,θ)=PWTPW0c_p(\lambda, \theta) = \dfrac{P_{WT}}{P_W^0} – aerodinamički koeficijent snage, za različite vrednosti ugla zahvatanja rotorskih lopatica (θ\theta);
  • λ=RωWTvW\lambda = \dfrac{R\omega_{WT}}{v_W} – odnos obimne brzine rotora turbine poluprečnika RR i brzine vetra vWv_W.

Na slici je prikazan tipičan izgled karakteristike aerodinamičkog koeficijenta snage cp(λ,θ)c_p(\lambda, \theta), za različite vrednosti ugla θ\theta.

Aerodinamički koeficijent snage c_p u funkciji odnosa brzina λ od 0 do 16, za uglove θ od 1°, 2°, 5°, 10°, 15° i 25°: najveći, oko 0,5, ima kriva za 1° pri λ blizu 7; što je ugao veći, to je vrh niži i bliže manjem λ (za 25° oko 0,1 pri λ oko 3).

Konstrukciono se vetrogeneratori (vetroturbina i generator) izvode na stubovima, nekoliko desetina metara (100 do 150 m) iznad tla, kao na Slici. Najčešće su potpuno automatizovani, bez stalne posade i upravljani daljinski.

Vetrogenerator na stubu: vetar pokreće lopatice rotora vetroturbine; u kućištu na vrhu stuba su prenosnik (kutija sa zupčanicima) i generator; energetski kabl se spušta kroz stub do transformatora za povezivanje sa mrežom, a od njega ide dalekovod.

Kako se potencijal vetra menja sa promenom visine, to se isti modeluje prema nekom poluiskustvenom nelinearnom zakonu i njegov model se uklapa u napred izloženi model vetroturbine. Za promene brzine vetra u radnom opsegu iznad brzine od 4 m/s, brzina vetra na visini hh iznad tla, nalazi se iz relacije

vW(h)=v10(hh10)a,v_W(h) = v_{10}\left(\frac{h}{h_{10}}\right)^a,

gde su:

  • a≈0,14−0,18a \approx 0{,}14 - 0{,}18 – Helmanov eksponent;
  • v10v_{10} – izmerena ili proračunata brzina vetra na visini h10=10h_{10} = 10 m iznad tla.

Iz predavanja: Novi i obnovljivi izvori električne energije: Vetroelektrane; Novi i obnovljivi izvori električne energije: Konstrukcione karakteristike vetrogeneratora.

16. Zavisnost odate snage vetroelektrane od brzine vetra i ekološka ograničenja u korišćenju vetroelektrana.

Kriva brzina vetra-odata snaga je osnovna karakteristika svakog vetroagregata, koja daje odatu snagu u funkciji brzine vetra: PWT=f(vW)P_{WT} = f(v_W). Tipičan izgled ove karakteristike dat je na Slici.

Kriva snage vetroagregata: relativna snaga P_WTG [%] u funkciji brzine vetra v_W [m/s]; snaga raste od startne brzine v_in (oko 1 m/s) do nazivne snage pri v_n (13 m/s), ostaje 100 % do najveće dozvoljene brzine v_out (25 m/s), kad pada na nulu; najveća startna brzina v_ret je 22 m/s.

Na Slici se uočavaju četiri karakteristične brzine vetra:

  • startna (vinv_{in}) – brzina pri kojoj turbina počinje da odaje snagu PWT>0P_{WT} > 0;
  • nazivna (vnv_n) – brzina pri kojoj turbina dostiže svoju nazivnu snagu (PWTnP_{WTn});
  • najveća startna (vretv_{ret}) – najveća brzina pri kojoj se vetroturbina može staviti u pogon;
  • najveća dozvoljena (voutv_{out}) – najveća brzina pri kojoj turbina još može da radi (obično oko 20–25 m/s).

Shodno Slici, neiskoristiv je energetski potencijal pri brzinama vetra vW<vinv_W < v_{in} i vW>voutv_W > v_{out}, a delimično iskoristiv u opsegu brzina vetra vin<vW<voutv_{in} < v_W < v_{out} i to na nivou vW=vnv_W = v_n, što se može zapisati kao:

PWTG={0za0≤vW≤vinf(vW)zavin≤vW≤vnPWTnzavn≤vW≤vout0zavW≥vout.P_{WTG} = \begin{cases} 0 & \text{za} \quad 0 \le v_W \le v_{in} \\ f(v_W) & \text{za} \quad v_{in} \le v_W \le v_n \\ P_{WTn} & \text{za} \quad v_n \le v_W \le v_{out} \\ 0 & \text{za} \quad v_W \ge v_{out} \end{cases}.

Ekološka ograničenja

Vetroelektrane ne proizvode potpuno čistu električnu energiju. Nema hemijskih zagađenja tla, vode i vazduha, ali ima prigovora da visoki stubovi (na koje se montiraju vetroagregati) vizuelno kvare pejzaž. Takođe, ornitolozi strahuju da mogu da budu uzrok pomora ptica selica, ako im se lokacija nađe na putu njihovih seoba u jesen i proleće. Najozbiljnija zamerka vetroelektranama je da svojim radom podižu nivo buke, zbog čega treba da budu locirane ako je moguće dalje od naselja.

Nivo buke [dB] u funkciji rastojanja od vetrogeneratora [m], od 0 do 900 m: kriva opada od najveće vrednosti na samoj lokaciji, preko oko 40 dB na 300 m, ka oko 30 dB na 900 m.

Vidi se da je nivo buke na samoj lokaciji vetroagregata oko 70 dB, a na rastojanju od 300 m oko 40 dB [6].

Smatra se da je stacionarna vrednost nivoa buke od 30 dB prihvatljiva.

Iz predavanja: Novi i obnovljivi izvori električne energije: Karakteristike performansi vetroelektrana; Novi i obnovljivi izvori električne energije: Ekološka ograničenja.

17. Osnovne šeme priključenja vetrogeneratora na mrežu.

Pošto je brzina obrtanja turbine (10–40 ob/min) višestruko manja od uobičajenih brzina obrtanja generatora (1500 i 3000 ob/min), spoj između turbine i generatora izvodi se preko mehaničke spojnice. Savremeni vetrogeneratori su konstrukciono višepolni, niskobrzinski, naizmenični sinhroni, ili asinhroni generatori. Jednosmerni generatori danas se koriste samo za vetrogeneratore malih snaga (reda nekoliko desetina kW). Kada se opreme sa akumulatorskim baterijama za deponovanje električne energije, takvi generatori su pogodni za napajanje manjih, izolovanih potrošača i potrošačkih područja.

Naizmenični generatori mogu raditi sa konstantnom, ili promenljivom brzinom obrtanja.

Na Slici je prikazana opšta blok-šema vetrogeneratorskog postrojenja, priključenog na elektroenergetski sistem.

Blok-šema vetrogeneratorskog postrojenja: kinetička energija vetra se u vetroturbini pretvara u mehaničku, a u električnom generatoru u električnu energiju; vetroturbina preko mehaničke spojnice (menjačke kutije) pokreće električni generator, iza kog su elektronski pretvarači AC-DC-AC, kompenzator reaktivne snage i oprema za spregu sa mrežom prema elektroenergetskom sistemu; podsistem upravljanja vetrogeneratora prima merenja i signalizacije, merenje brzine vetra anemometrom i komande, a upravlja regulacijom ugla zakretanja elisa, regulacijom pobude i zaštitom.

Osnovne principijelne šeme priključenja na mrežu date su u Tabeli.

Asinhroni generatorSinhroni generator
a. Direktan priključak generatora na mrežu
Vetroturbina preko menjačke kutije (GB) pokreće asinhroni generator (AG), priključen direktno na mrežu 50 Hz.n=(1−s)F/pn = (1 - s)F/p; s=0⋯−0,08 {s = 0 \cdots -0{,}08}Potrošač reaktivne snagea1Vetroturbina preko menjačke kutije (GB) pokreće sinhroni generator (SG) sa sopstvenom pobudom, priključen direktno na mrežu 50 Hz.n=F/pn = F/pKontrolabilna reaktivna snagaa2
b. Priključak generatora na mrežu posredstvom AC-DC-AC pretvarača
Vetroturbina preko menjačke kutije (GB) pokreće asinhroni generator (AG), priključen na mrežu 50 Hz preko pretvarača AC-DC-AC sa jednosmernim (DC) međukolom.n=(0,8−1,2)F/pn = (0{,}8 - 1{,}2)F/pKontrolabilna reaktivna snagab1aVetroturbina preko menjačke kutije (GB) pokreće sinhroni generator (SG) sa sopstvenom pobudom, priključen na mrežu 50 Hz preko pretvarača AC-DC-AC sa jednosmernim (DC) međukolom.n=(0,5−1,2)F/pn = (0{,}5 - 1{,}2)F/pKontrolabilna reaktivna snagab1s

Iz predavanja: Novi i obnovljivi izvori električne energije: Osnovne šeme priključenja vetrogeneratora na mrežu.

18. Uklapanje vetroelektrana u elektroenergetski sistem (dozvoljeni stepen penetracije i farme vetrogeneratora).

Električna energija koju proizvode vetroelektrane po svojoj prirodi je stohastička i neupravljiva. Iako tehnika savremenih vetrogeneratora predstavlja prilično usavršenu tehnologiju, a osnovni energetski resurs (vetar) jedan od najrasprostranjenijih obnovljivih energenata, vetroelektrane nisu sposobne da samostalno, bez pomoći drugih vrsta izvora (ili akumulatora) električne energije, obezbede permanentno napajanje potrošača. Planiranje proizvodnje vetroelektrana moguće je samo na kratak rok, jer je predviđanje srednje brzine i trajanja vetra praktično nemoguće na dužim vremenskim horizontima. Zato je masovna primena vetroelektrana jedino moguća u elektroenergetskim sistemima gde druge, klasične elektrane (hidro i termo), moraju obezbediti operativnu i regulacionu rezervu, koju vetroelektrane ne mogu da obezbede.

U tom smislu definiše se dozvoljeni stepen penetracije kao količnik raspoložive snage (ili proizvedene energije) angažovanih vetroelektrana i ukupne snage (proizvedene električne energije) svih elektrana sistema, tj. on se može definisati na bazi snage i energije kao

pW=(PWT)a(Pgs)a;wW=(WWT)a(Wgs)a,p_W = \frac{(P_{WT})_a}{(P_{gs})_a}; \quad w_W = \frac{(W_{WT})_a}{(W_{gs})_a},

gde su: (PWT)a=∑i=1NWPWTi(P_{WT})_a = \sum_{i=1}^{N_W} P_{WTi}; (Pgs)a=∑i=1NPgi(P_{gs})_a = \sum_{i=1}^{N} P_{gi}, a NWN_W i NN broj vetroelektrana i ukupan broj elektrana u elektroenergetskom sistemu.

Smatra se da nije poželjno (ne i nemoguće) da stepen penetracije pređe 10 %.

Farme vetrogeneratora (vetroparkovi) i njihov priključak na elektroenergetsku mrežu

Jedinične snage savremenih vetrogeneratora su relativno male, pa se oni grade za niske nazivne napone (do 1000 V). Sasvim mali generatori, nazivnih snaga reda 100 kW, obično se direktno priključuju na niskonaponsku distributivnu mrežu 400/231 V. Ako se radi o vetrogeneratorima veće snage, kada ih je više, lociranih na bliskom geografskom prostoru, oni se udružuju u farmu vetrogeneratora (vetroparkove), koja se u elektroenergetskom sistemu tretira kao jedinstvena elektrana. Ako je instalisana snaga takvih elektrana između 100 kW i nekoliko MW, priključak se izvodi na srednjenaponsku distributivnu mrežu, nazivnog napona između 10 i 66 kV. Velike farme vetrogeneratora, snaga nekoliko desetina MW, preko odgovarajućih transformatora priključuju se direktno na prenosnu mrežu (čiji je nazivni napon iznad 90 kV). Nacionalnim i međunarodnim propisima reguliše se način priključka vetroelektrana na elektroenergetsku mrežu, koji zavisi od niza tehničkih uslova.

Iz predavanja: Novi i obnovljivi izvori električne energije: Uklapanje vetroelektrana u elektroenergetski sistem; Novi i obnovljivi izvori električne energije: Farme vetrogeneratora (vetroparkovi) i njihov priključak na elektroenergetsku mrežu.

19. Solarne elektrane (pojam, vrste i tehničke/konstrukcione karakteristike, kao i proračun odate snage – bez izlaganja detaljnih izraza).

Korišćenje energije sunca za proizvodnju električne energije može biti dvojako: korišćenjem solarnih koncentratora za proizvodnju vodene pare, koja pokreće parnu turbinu i sinhroni generator (klasična termoelektrana na paru), ili direktnim pretvaranjem solarne u električnu energiju, preko fotonaponskih pretvarača. Taj proces je osnov na kome počivaju solarne fotonaponske elektrane.

Osnov fotonaponskog efekta čini pobuđivanje i pokretanje elektrona, kada neki materijali apsorbuju dovoljno fotona svetlosti određene energije. Osnovna proizvodna jedinica fotonaponskog sistema je modul, koji se sastoji od tankog sloja fotoosetljivog poluprovodničkog materijala, sa čije obe strane su navučeni slojevi elektroprovodnog materijala. Stepen korisnog dejstva savremenih fotonaponskih modula je prilično nizak, odnosno 12–15 %.

Tipična naponsko-strujna karakteristika fotonaponskog modula prikazana je na Slici, za tri vrednosti specifične snage sunčevog zračenja. Maksimalna snaga modula ima se ako je radna tačka u blizini kolena V-I karakteristike. Treba napomenuti da izlazni napon opada sa porastom temperature, tako da je položaj tačke maksimalne izlazne snage promenljiv.

V-I karakteristike fotonaponskog modula za 1000, 700 i 500 W/m²: struja I [A] je skoro stalna (4 A, 1 A i niže) dok napon raste, pa posle kolena pada na nulu pri 21 V; strelica pokazuje pad napona pri porastu temperature.

U opštem slučaju, fotonaponski generatori su sastavljeni od nekoliko standardnih modula, prikazanih na Slici, povezanih u Nsg×NpgN_{sg} \times N_{pg} serijsko-rednih konfiguracija, gde su NsgN_{sg} i NpgN_{pg} ukupni brojevi serijskih i rednih modula, respektivno.

Zamenska šema fotonaponskog modula: strujni izvor, kome zračenje ℜ i temperatura modula T_m dolaze iz vremenske prognoze, paralelno sa diodama D1 i D2, zatim redna otpornost R_sm; na krajevima modula su struja I_m, napon V_m i snaga P_m^DC.

Strujno-naponska karakteristika fotonaponskog generatora je:

Ig=NpgIm;Vg=NsgVm,\begin{gathered} I_g = N_{pg}I_m; \\ V_g = N_{sg}V_m, \end{gathered}

gde su IgI_g (VgV_g) i ImI_m (VmV_m) struje (naponi) ukupnog fotonaponskog generatora ili pojedinačnog modula, respektivno.

Izlazna snaga fotonaponskog modula (PmDCP_m^{DC}) dobija se jednostavno množenjem izlazne struje (ImI_m) i napona na njegovim krajevima (VmV_m).

Može se pretpostaviti da je model invertora sličan modelu idealnog strujnog izvora, sa nultom DC komponentom i nultom ukupnom harmonijskom distorzijom izlaznog napona. Pod tim pretpostavkama, izlazna snaga iz DC/AC invertora je

PmAC=ηPmDC,P_m^{AC} = \eta P_m^{DC},

gde većina komercijalnih invertora imaju efikasnost (η\eta) između 90 i 96 %.

Iz predavanja: Novi i obnovljivi izvori električne energije: Solarne termo i fotonaponske elektrane.

Drugi deo

20. Pojam i vrste elektroenergetskih bilansa. Primer(i) popunjavanja površine ispod krive trajanja opterećenja.

Bilansna jednačina za snagu može se predstaviti kao:

Bilans snaga u jednoj tački: ulaze raspoloživi kapaciteti Pgr i dodati kapaciteti PgD, izlaze maksimalna potrošnja na pragu PpM i rezerva Prez, a isprekidano ulazi neisporučena snaga (redukcija) Pred.
Pgr+PgD≥PpM+Prez−PredP_{gr} + P_{gD} \ge P_{pM} + P_{rez} - P_{red}

Bilansna jednačina za energiju može se predstaviti kao:

Bilans energije u jednoj tački: ulaze proizvodnja postojećih elektrana Wg i proizvodnja dodatnih izvora WgD, izlazi potrošnja na pragu generatora Wp, a isprekidano ulazi neisporučena energija Wred.
Wg+WgD=Wp−WredW_g + W_{gD} = W_p - W_{red}

Bilans je sučeljavanje prihodne i rashodne strane fizičke veličine. U terminologiji elektroenergetskog bilansa sučeljavaju se proizvodnja i potrošnja električne energije i raspoložive snage.

Polazni osnov su dijagrami (krive trajanja) opterećenja.

Mora se obezbediti pouzdanost napajanja, zahtevana rezerva, kvalitet rada i slično.

Primer popunjavanja površine ispod krive trajanja opterećenja:

Kriva trajanja opterećenja Pp [MW] u funkciji vremena t [h], površina ispod nje podeljena na tri dela. Vršni deo: gasnoturbinske TE, akumulacione HE, reverzibilne HE. Srednji deo: TE na paru (čvrsta goriva), TE na paru (prirodni gas), TE na paru (tečna goriva). Osnovni (bazni) deo: protočne HE, TE-toplane, nuklearne TE, uvoz po dugoročnim ugovorima. TE – termoelektrane, HE – hidroelektrane.

Osnovni princip popunjavanja krive trajanja opterećenja je da se prvo angažuju proizvodne jedinice sa minimalnim proizvodnim troškovima:

1. Dno površine ispod krive trajanja popunjava uvoz po dugoročnim ugovorima (od ove vrednosti može se oduzeti izvoz po dugoročnim ugovorima, ili se jedinice za uvoz ne uzimaju u obzir pri bilansiranju).

2. U baznom delu krive trajanja angažuju se protočne HE, pošto nemaju regulacione mogućnosti, kao i biološki minimum akumulacionih HE.

3. Srednji deo ispod krive trajanja (najveći kod nas) popunjavaju klasične TE. Redosled angažovanja je od najjeftinije ka najskupljoj. Najlakši način je lista prioriteta po specifičnim troškovima.

4. Vršni deo ispod krive trajanja popunjava se akumulacionim HE. Osnovni cilj je da se što više iskoriste proizvodne mogućnosti.

5. Pumpno-akumulacione (reverzibilne) HE dolaze u sam vrh krive trajanja. Kod njih se mora voditi računa i o potrošnji u pumpnom režimu.

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta; Planiranje razvoja proizvodnih kapaciteta: Elektroenergetski bilansi; Planiranje razvoja proizvodnih kapaciteta: Angažovanje proizvodnih jedinica.

21. Pouzdanost proizvodnog podsistema, rezerva proizvodnih kapaciteta i faktor rezerve. Načini određivanja rezerve u karakterističnim tipovima elektroenergetskih sistema.

Osnovni cilj je da se odredi pouzdanost napajanja potrošača električnom energijom.

Osnovni kriterijum je rezerva proizvodnih kapaciteta, koja se definiše kao:

Rr=(Generatorska snaga u pogonu)−(Vrsˇna snaga)(Vrsˇna snaga)R_r = \frac{(\text{Generatorska snaga u pogonu}) - (\text{Vršna snaga})}{(\text{Vršna snaga})}

Nekada se daje u odnosu na: (Generatorsku snagu u pogonu)

Faktor rezerve u sistemu je:

ρM=Generatorska snaga u pogonuVrsˇna snaga\rho_M = \frac{\text{Generatorska snaga u pogonu}}{\text{Vršna snaga}}

Rezerva se koristi za: pokrivanje snage ispada, popravljanje grešaka u prognozi opterećenja, preraspodelu proizvodnje pri ispadima prenosnih kapaciteta, nedovoljnu pouzdanost interkonekcije ...

Načini određivanja rezerve

1. Dominantan udeo TE sa propisanim nivoom rezerve:

Propisuje se nivo rezerve u odnosu na godišnje vršno opterećenje (obično 15–25 %). Pristup je neosetljiv na intenzitet prinudnih ispada, veličinu agregata itd.

Kod HE mora se voditi računa i o rezervi u energiji.

2. Deterministički proračun rezerve na bazi jednostrukog ispada

Analizira se jednostruki ispad najveće generatorske jedinice, uz dodatak neke snage nezavisne od veličine ispada (Prez=PgM+P0P_{rez} = P_{gM} + P_0).

Ovaj metod ne vodi računa o intenzitetima prinudnih ispada.

3. Stohastički metod

Koristi se gde je dominantan udeo TE.

Bazira se na razmatranju verovatnoće simultanog ispada i modela dnevnih maksimuma.

Cilj je da se odredi broj dana u godini kada se očekuje manjak generatorske snage.

Ovaj indeks se naziva LOLP-indeks, ili indeks gubitka opterećenja.

LOLP (Loss of Load Probability) – verovatnoća gubitka opterećenja (prekida isporuke električne energije potrošačima).

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Pouzdanost proizvodnog podsistema; Planiranje razvoja proizvodnih kapaciteta: Načini određivanja rezerve.

22. LOLP indeks: pojam i način proračuna.

Bazira se na razmatranju verovatnoće simultanog ispada i modela dnevnih maksimuma.

Cilj je da se odredi broj dana u godini kada se očekuje manjak generatorske snage.

Ovaj indeks se naziva LOLP-indeks, ili indeks gubitka opterećenja.

LOLP (Loss of Load Probability) – verovatnoća gubitka opterećenja (prekida isporuke električne energije potrošačima).

Kriva trajanja dnevnih maksimuma Pg [MW] u funkciji dana [r.j.] do 1,0. PgrM je raspoloživa snaga, Prez potrebna rezerva iznad vrha krive, a ΔPgi ispad i-te generatorske jedinice. Prava PgrM − ΔPgi seče krivu, a Di je broj dana u kojima je kriva iznad nje; ta površina je osenčena.

Na prethodnoj slici DiD_i je broj dana u godini sa manjkom generatorskih kapaciteta tokom godine.

LOLP-indeks se sastoji iz dva dela:

  • Neraspoloživost generatorskih kapaciteta, okarakterisana intenzitetom prinudnih ispada (FOR – Forced Outage Rate).
  • Neraspoloživost generatorskih kapaciteta, okarakterisana intenzitetom planskih ispada (SOR – Scheduled Outage Rate).

Osnovna pretpostavka proračuna LOLP-indeksa je da je kriva trajanja deterministička veličina, dok su proizvodne jedinice međusobno nezavisne slučajne veličine.

Za ove pretpostavke stanje raspoloživosti definiše se samo na osnovu generatorskih jedinica u pogonu ili ispadu. Pri tome, obično se zanemaruju hidro-jedinice, jer su one pouzdanije od termičkih (po prirodi svog proizvodnog procesa).

Jedno stanje raspoloživosti određuju:

  1. veličina snage u ispadu i
  2. verovatnoća nastanka stanja.

Izraz za verovatnoću pojave stanja dobija se na osnovu metoda prebrojavanja stanja (prebrojati sva moguća stanja)

℘s=∏i=1N−Nraspqi∏i=1Nrasppi\wp_s = \prod_{i=1}^{\mathbb{N} - \mathbb{N}_{rasp}} q_i \prod_{i=1}^{\mathbb{N}_{rasp}} p_i

gde su:

  • N\mathbb{N} – ukupan broj mogućih komponenata u sistemu (određene jednim stanjem, koje može biti komponenta raspoloživa, ili komponenta neraspoloživa);
  • Nrasp\mathbb{N}_{rasp} – broj komponenata u stanju raspoloživosti (ispravnog rada);
  • N−Nrasp\mathbb{N} - \mathbb{N}_{rasp} – broj komponenata u stanju neraspoloživosti (ispada);
  • pip_i, qiq_i – verovatnoće pojave stanja raspoloživosti i stanja neraspoloživosti.

Neraspoloživost, određena intenzitetom prinudnih ispada (FOR – "Forced Outage Rate"), neke komponente izračunava se kao

qi=λiλi+μi=riri+mi=riTi=fiμi[r.j.],q_i = \frac{\lambda_i}{\lambda_i + \mu_i} = \frac{r_i}{r_i + m_i} = \frac{r_i}{T_i} = \frac{f_i}{\mu_i} \quad [\rj],

dok se raspoloživost izračunava kao

pi=μiλi+μi=miri+mi=miTi=fiλi[r.j.],p_i = \frac{\mu_i}{\lambda_i + \mu_i} = \frac{m_i}{r_i + m_i} = \frac{m_i}{T_i} = \frac{f_i}{\lambda_i} \quad [\rj],

gde su:

  • λi\lambda_i – intenzitet kvarova i-te komponente u [kvarova/god];
  • μi\mu_i – intenzitet obnavljanja i-te komponente u [obnavljanja/god];
  • ri=1/μir_i = 1/\mu_i – srednje vreme trajanja kvarova i-te komponente u [god];
  • mi=1/λim_i = 1/\lambda_i – srednje vreme do pojave kvara i-te komponente u [god];
  • Ti=ri+miT_i = r_i + m_i – srednje vreme između dva ispada i-te komponente u [god];
  • fi=1/(ri+mi)=1/Tif_i = 1/(r_i + m_i) = 1/T_i – srednja učestanost ispada i-te komponente u [ispada/god].

Godišnji LOLP-indeks izračunava se kao zbirna verovatnoća ukupnog očekivanog broja dana (DD) sa manjkom generatorskih kapaciteta (ΔPg\Delta P_g):

LOLP=℘(ΔPg)D=∑i℘[ΔPgi(rj)]Di,LOLP = \wp(\Delta P_g)D = \sum_i \wp[\Delta P_{gi}(r_j)]D_i,

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Načini određivanja rezerve.

23. Primer proračuna LOLP indeksa u sistemu sa tri proizvodne jedinice.

Primer:

Neka se raspolaže sa tri tipa agregata j = A, B i C, čiji je ukupan broj nAn_A, nBn_B i nCn_C. Njihova pojedinačna raspoloživost snaga je PgrAP_{grA}, PgrBP_{grB} i PgrCP_{grC} i verovatnoće prinudnih ispada qAq_A, qBq_B i qCq_C. Verovatnoće pojave ispada pojedinih tipova (gde ispada rAr_A, rBr_B i rCr_C agregata) su:

℘(rA)=(nArA)qArApAnA−rA;ΔPg(rA)=rAPgrA;\wp(r_A) = \binom{n_A}{r_A}q_A^{r_A}p_A^{n_A - r_A};\quad \Delta P_g(r_A) = r_AP_{grA};
℘(rB)=(nBrB)qBrBpBnB−rB;ΔPg(rB)=rBPgrB;\wp(r_B) = \binom{n_B}{r_B}q_B^{r_B}p_B^{n_B - r_B};\quad \Delta P_g(r_B) = r_BP_{grB};
℘(rC)=(nCrC)qCrCpCnC−rC;ΔPg(rC)=rCPgrC.\wp(r_C) = \binom{n_C}{r_C}q_C^{r_C}p_C^{n_C - r_C};\quad \Delta P_g(r_C) = r_CP_{grC}.

Verovatnoća simultanog ispada rAr_A agregata tipa A i rBr_B agregata tipa B je:

℘(rA+rB)=℘(rA)℘(rB);ΔPg(rA+rB)=rAPgrA+rBPgrB.\begin{gathered} \wp(r_A + r_B) = \wp(r_A)\wp(r_B); \\ \Delta P_g(r_A + r_B) = r_AP_{grA} + r_BP_{grB}. \end{gathered}

Verovatnoća simultanog ispada sva tri agregata je:

℘(rA+rB+rC)=℘(rA)℘(rB)℘(rC);ΔPg(rA+rB+rC)=rAPgrA+rBPgrB+rCPgrC.\begin{gathered} \wp(r_A + r_B + r_C) = \wp(r_A)\wp(r_B)\wp(r_C); \\ \Delta P_g(r_A + r_B + r_C) = r_AP_{grA} + r_BP_{grB} + r_CP_{grC}. \end{gathered}

Prethodni izrazi su verovatnoće pojave dva/tri nezavisna događaja. One daju jedno od mogućih stanja događaja, prikazano u tabeli (nA=1,nB=1n_A = 1, n_B = 1 i nC=1n_C = 1):

StanjeStatus proizvodnog agregataIspala snagaVerovatnoća pojave i-tog stanja
iABCΔPgi\Delta P_{gi}
1+++0(1−qA)(1−qB)(1−qC)(1 - q_A)(1 - q_B)(1 - q_C)
2++–PgrCP_{grC}(1−qA)(1−qB)qC(1 - q_A)(1 - q_B)q_C
3+–+PgrBP_{grB}(1−qA)qB(1−qC)(1 - q_A)q_B(1 - q_C)
4–++PgrAP_{grA}qA(1−qB)(1−qC)q_A(1 - q_B)(1 - q_C)
5+––PgrB+PgrCP_{grB} + P_{grC}(1−qA)qBqC(1 - q_A)q_Bq_C
6–+–PgrA+PgrCP_{grA} + P_{grC}qA(1−qB)qCq_A(1 - q_B)q_C
7––+PgrA+PgrBP_{grA} + P_{grB}qAqB(1−qC)q_Aq_B(1 - q_C)
8–––PgrA+PgrB+PgrCP_{grA} + P_{grB} + P_{grC}qAqBqCq_Aq_Bq_C
+ Agregat u pogonu ℘=p=1−q\wp = p = 1 - q– Agregat van pogona ℘=q\wp = qZbir: 1,00

Broj dana sa deficitom snage (DiD_i) očitava se sa godišnje krive trajanja opterećenja.

Na osnovu broja dana u godini sa deficitom proračunava se LOLP-indeks na već opisan način.

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Načini određivanja rezerve.

24. Osnovni i dopunski indeksi pouzdanosti koji se koriste u planiranju proizvodnih kapaciteta.

LOLP-indeks daje verovatnoću pojave deficita, ali ne daje informaciju o trajanju deficita.

Zato se koriste i drugi indeksi:

1. Očekivana neisporučena energija, odnosno LOE-indeks (LOE – Loss of Energy).

2. Indeks verovatnoće gubitaka energije, odnosno LOEP-indeks (LOEP – Loss of Energy Probability), kojim se definiše mera efekta deficita snage. Na godišnjem nivou definiše se kao:

(LOEP)g=℘⋅(LOE)=∑i=1Ng℘i⋅(LOE)i[r.j./god](\mathrm{LOEP})^{\mathrm{g}} = \wp \cdot (\mathrm{LOE}) = \sum_{\mathrm{i}=1}^{\mathbb{N}_{\mathrm{g}}} \wp_{\mathrm{i}} \cdot (\mathrm{LOE})_{\mathrm{i}} \quad [\text{r.j./god}]

3. Indeks učestanosti i trajanja, odnosno FD-indeks (FD – Frequency and Duration).

Pored stacionarnih verovatnoća (ne)raspoloživosti elemenata (pp i qq), za proračun (F&D)-indeksa neophodno je i poznavanje faktora napuštanja (tranzicije) s-tog stanja (τsk\tau_{sk}), koji određuju intenzitet kvarova (λk\lambda_k) i/ili obnavljanja (μk\mu_k) pojedinih (k-tih) prelaza iz s-tog stanja. Na osnovu svih mogućih (k-tih) prelaza vrši se proračun ukupnog faktora napuštanja s-tog stanja (∑kτsk\sum\limits_k \tau_{sk}), a onda i indeks učestanosti pojave s-tog stanja

fs=℘s∑kτskf_s = \wp_s \sum_k \tau_{sk}

i trajanja s-tog stanja

ds=1∑kτskd_s = \frac{1}{\sum\limits_k \tau_{sk}}

Indeks učestanosti i trajanja izračunava očekivano trajanje deficita generatorskih kapaciteta i očekivanu učestanost tog deficita tokom godine.

4. Srednje vreme trajanja kvarova i Srednje vreme trajanja stanja.

Na osnovu ovih indeksa izračunava se vreme između dve uzastopne pojave nekog stanja:

(Vreme ponavljanja stanja)=1(Ucˇestanost pojavljivanja stanja)(\text{Vreme ponavljanja stanja}) = \frac{1}{(\text{Učestanost pojavljivanja stanja})}
(Ocˇekivano srednje vreme trajanja stanja)=(Verovatnocˊa stanja)(Ucˇestanost pojavljivanja stanja)\begin{gathered} (\text{Očekivano srednje vreme trajanja stanja}) \\ = \frac{(\text{Verovatnoća stanja})}{(\text{Učestanost pojavljivanja stanja})} \end{gathered}

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Ostali indeksi pouzdanosti.

25. Deterministički izbor veličine dodatog proizvodnog agregata.

Najlakši način izbora je deterministički pristup: Prvo se nacrta dijagram vršnih snaga u dve susedne godine perioda planiranja kao prava, kao na slici:

Vršne snage Pg [MW] u godinama k i k + 1 kao dve prave. U godini k, od ose naviše: PpM^k, Prez^(k−1) i PgD^k; u godini k + 1: PpM^(k+1) i Prez^k, do istog vrha.

Sa prethodne slike, ukupna raspoloživa snaga u k-toj godini je:

Pgrk≥PpMk+Prezk−1+PgDk≥PpMk+1+Prezk,P_{gr}^k \ge P_{pM}^k + P_{rez}^{k-1} + P_{gD}^k \ge P_{pM}^{k+1} + P_{rez}^k,

odakle je:

PgDk≥(PpMk+1−PpMk)+(Prezk−Prezk−1)=nkPgjk;Pgjk≤βPpMk,\begin{gathered} \boxed{P_{gD}^k \ge \left(P_{pM}^{k+1} - P_{pM}^k\right) + \left(P_{rez}^k - P_{rez}^{k-1}\right) = n^kP_{gj}^k}; \\ P_{gj}^k \le \beta P_{pM}^k, \end{gathered}

gde koeficijent β\beta definiše relativnu snagu dodatog agregata.

Za praktičnu primenu uvode se neka uprošćenja:

1. Relativni porast opterećenja u dve susedne godine je konstantan:

PpMk+1PpMk=b,\frac{P_{pM}^{k+1}}{P_{pM}^k} = b,

gde je b pretpostavljena konstanta.

2. Ista pretpostavka kao prethodna važi i za rezerve:

PrezkPrezk−1=b.\frac{P_{rez}^k}{P_{rez}^{k-1}} = b.

3. Odnos rezerve u odnosu na snagu najvećeg agregata:

Prezk=aPgjk.P_{rez}^k = aP_{gj}^k.

Kada se zamene prethodne pretpostavke u jednačinu bilansa dobija se:

PgDk=PpMk(b−1)+Pgjka(1−1b)=(b−1)(PpMk+abPgjk)=nkPgjk.\begin{aligned} P_{gD}^k &= P_{pM}^k(b - 1) + P_{gj}^ka\left(1 - \frac{1}{b}\right) \\ &= (b - 1)\left(P_{pM}^k + \frac{a}{b}P_{gj}^k\right) = n^kP_{gj}^k. \end{aligned}

To znači da je snaga najvećeg agregata:

Pgjk=PpMknkb−1−ab=PpMkb−1nk−ab(b−1).\boxed{P_{gj}^k = \frac{P_{pM}^k}{\dfrac{n^k}{b - 1} - \dfrac{a}{b}} = P_{pM}^k\frac{b - 1}{n^k - \dfrac{a}{b}(b - 1)}}.

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Izbor veličine proizvodnih agregata.

26. Objasniti osnovne korake u studiji dogradnje proizvodnih kapaciteta.

Osnovni koraci u planiranju razvoja izvora su:

  • Analiza proizvodnih mogućnosti postojećih kapaciteta.
  • Određivanje nivoa pouzdanosti isporuke električne energije.
  • Proračun dodatnih kapaciteta: – Izgradnja novih kapaciteta da se zadovolji elektroenergetski bilans. – Potrebne rezerve (po veličini i strukturi).
  • Izrada plana gašenja i revitalizacije starih izvora (elektrana).
  • Izrada dinamičkog plana izgradnje.
  • Izbor lokacije.
  • Obezbeđenje finansijskih sredstava.
  • Potrebna proširenja prenosne mreže.

Ovo je najsloženija faza planiranja razvoja proizvodnih kapaciteta. Donosi se odluka o dogradnji ili izgradnji novih elektrana. To je najteži slučaj u praksi.

Po pravilu se postavlja kao optimizacioni problem:

Minimum investicionih i eksploatacionih troškova, uz zadovoljenje tehničko-ekonomskih kriterijuma.

Pri rešavanju prethodnog problema postoje dva tipa neizvesnih ulaznih podataka:

  • Buduće vrednosti prognozirane potrošnje (prognozirana kriva trajanja opterećenja). Modeluje se deterministički.
  • Prinudni ispadi termičkih jedinica (hidro jedinice su pouzdanije, tako da se njihovi ispadi ne razmatraju). Modeluje se probabilistički.

Studija dogradnje sprovodi se kroz sledeća dva koraka:

1. Popunjavanje površine ispod (determinističke) krive trajanja opterećenja. Kriterijum je minimum troškova sistema. Ovo je u stvari elektroenergetski bilans.

2. Proračun pokazatelja pouzdanosti, na osnovu kojih se zaključuje o adekvatnosti razvoja sistema.

Proračunavaju se:

  • Verovatnoća neisporučenog opterećenja (LOLP-indeks).
  • Očekivana neisporučena energija (EUE-indeks).

Ovi pokazatelji se proračunavaju na napred izloženi način, koji se skraćeno ponavlja:

Isti dijagram kao ranije, kriva trajanja dnevnih maksimuma Pg [MW] sa rezervom Prez, ispadom ΔPgi i raspoloživom snagom PgrM, a na apscisi i dani [r.j.] sa Di i 1,0 i vreme t [h] sa ti i T.

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta; Planiranje razvoja proizvodnih kapaciteta: Studija dogradnje proizvodnih kapaciteta.

27. Vrste troškova proizvodnih kapaciteta.

Investicioni troškovi:

  • Otplata zajma (kredita).
  • Kamate na pozajmljena sredstva.
  • Porezi, doprinosi, osiguranje ...
CIg=cIgCI,C_I^g = c_I^gC_I,

gde je cIgc_I^g godišnja stopa investicionih troškova.

Eksploatacioni troškovi:

CExg=∑HECOMg+∑TE(COMg+Cgg),C_{Ex}^g = \sum_{HE} C_{OM}^g + \sum_{TE}\left(C_{OM}^g + C_g^g\right),

gde je COMg=COMog+COMvgC_{OM}^g = C_{OM}^{og} + C_{OM}^{vg} (zbir fiksnih i varijabilnih troškova).

Ukupni troškovi u sistemu:

Cg=CIg+CExg=CFg+CVarg.C^g = C_I^g + C_{Ex}^g = C_F^g + C_{Var}^g.

Prosečni troškovi godišnje proizvedene električne energije:

cw=CgWg,c_w = \frac{C^g}{W^g},

gde je WgW^g očekivana godišnja proizvodnja na pragu elektrana.

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Troškovi proizvodnih agregata.

28. Optimizacioni postupci za izradu optimalnog vremenskog plana stavljanja u pogon.

Ovde se analizira i gašenje postojećih agregata (odsluženih, amortizovanih). Pošto se radi o vremenskom nizu, mora se izvršiti svođenje novčanih tokova, primenom nekog od ranije proučenih postupaka:

1. Metode ekvivalentnih vrednosti

  • Sadašnja ekvivalentna vrednost (PV)
  • Godišnja ekvivalentna vrednost (AV)
  • Buduća ekvivalentna vrednost (FV)

2. Profitne metode

  • Metod interne stope povraćaja (IRR)
  • Metod koristi i troškova (BC)

Očigledno je da vremenska dinamika stavljanja u pogon ima veliki uticaj.

Metod sadašnje ekvivalentne vrednosti

Sadašnja ekvivalentna vrednost fiksnih troškova je:

PVCF=∑k=1n∑j=1m[CIjg+COMjog(PVF)k]=∑j=1m[nCIjg+COMjog(PWF)],\begin{aligned} PVC_F &= \sum_{k=1}^{n}\sum_{j=1}^{m}\left[C_{Ij}^g + C_{OMj}^{og}(PVF)_k\right] \\ &= \sum_{j=1}^{m}\left[nC_{Ij}^g + C_{OMj}^{og}(PWF)\right], \end{aligned}

gde su:

  • nn – indeks godine;
  • jj – indeks agregata;
  • (PVF)k=(1+i)−k(PVF)_k = (1 + i)^{-k} – faktor svođenja na sadašnju ekvivalentnu vrednost;
  • (PWF)=(1+i)n−1i(1+i)n(PWF) = \dfrac{(1 + i)^n - 1}{i(1 + i)^n} – faktor svođenja uniformnog niza na sadašnju ekvivalentnu vrednost;

Ukupni troškovi su:

(PVC)=(PVCF)+(PVCVar).(PVC) = (PVC_F) + (PVC_{Var}).

Optimalno rešenje dobija se rešavanjem za h mogućih opcija:

min⁡h{(PVC)h}.\min_h\left\{(PVC)_h\right\}.

Metod poređenja koristi i troškova

Često se naziva i BC-metod (Benefit to Cost method).

Određuje se veličina snage za koju se umanjuje proizvodni kapacitet, kao i godišnja dobit (BgB^g). Ta dobit se poredi sa troškovima (CgC^g), kroz kriterijume Bg−Cg>0B^g - C^g > 0 ili Bg/Cg>1B^g/C^g > 1.

Efekti izgradnje novog izvora manifestuju se kroz isporuku električne energije i obezbeđenje potrebne snage potrošačima:

Bg=BkWg+BkWhg=aPEn+bWEg.B^g = B_{kW}^g + B_{kWh}^g = aP_{En} + bW_E^g.

Iz predavanja: Planiranje razvoja proizvodnih kapaciteta: Izrada vremenskog plana stavljanja u pogon; Planiranje razvoja proizvodnih kapaciteta: Metod sadašnje ekvivalentne vrednosti; Planiranje razvoja proizvodnih kapaciteta: Metod poređenja koristi i troškova.

29. Troškovi vazdušnih vodova (investicioni i eksploatacioni). Način proračuna i njihove aktualizacije.

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.

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.

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.

Postoje dve vrste troškova:

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

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}

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).

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}}.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Troškovi vazdušnih vodova; Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Eksploatacioni troškovi vazdušnih vodova.

30. Potrebni uslovi za optimalan prenos snage između dve tačke. Izvesti uslove za optimizaciju poprečnog preseka u formi Osnovnog Kelvin-ovog pravila.

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}

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.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Optimizacija prenosa između dve tačke; Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Optimizacija poprečnog preseka provodnika – Kelvin-ova pravila.

31. Potrebni uslovi za optimalan prenos snage između dve tačke. Izvesti uslove za optimizaciju poprečnog preseka u formi Modifikovanog Kelvin-ovog pravila.

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.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Optimizacija poprečnog preseka provodnika – Kelvin-ova pravila.

32. Opteretivost prenosnih vodova. Izvesti uslov za aktivno opterećenje prenosnog voda u formi detaljnijeg i maksimalno uprošćenog izraza.

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

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

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}}.

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

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

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.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Opteretivost prenosnih vodova.

33. Primena DC modela u planiranju prenosne mreže (pretpostavke primene uprošćenog DC modela tokova snaga). Promena tokova snaga u zavisnosti od promene injektiranja.

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],

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}

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},

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,
Δ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.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: DC model tokova snaga; Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Faktori preraspodele tokova snaga po granama pri promeni injektiranja.

34. Primena DC modela u planiranju prenosne mreže (primena DC modela i izvođenje faktora preraspodele tokova snaga po granama pri promeni injektiranja).

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

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.

Δ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.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Faktori preraspodele tokova snaga po granama pri promeni injektiranja; Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Izvođenje faktora osetljivosti.

35. Primena DC modela u planiranju prenosne mreže (izvođenje distribucionog faktora 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}.

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}}.

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)},

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Distribucioni faktori preraspodele tokova snaga po granama pri uključenju/isključenju grane.

36. Primena DC modela u planiranju prenosne mreže (osetljivost promene tokova snaga po granama pri simultanom uključenju/isključenju 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.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Osetljivost promene tokova snaga po granama pri simultanom ispadu grane i promeni injektiranja u čvorovima.

37. Primena metoda sukcesivne ekspanzije za planiranje prenosne mreže.

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},

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.

Č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\}.

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Metod sukcesivne ekspanzije.

38. Primena linearnog programiranja za planiranje prenosne mreže.

Egzaktnije rešenje je primena optimizacionih tehnika, od kojih je najveću primenu našlo linearno programiranje.

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.

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).

Iz predavanja: Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Primena linearnog programiranja za planiranje prenosnih mreža; Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Model postojeće mreže (označen indeksom P); Planiranje razvoja prenosnih mreža i studija dogradnje reaktivnih izvora: Model buduće proširene mreže (postojeći i novododati vodovi).

Teme

Planiranje · Potrošnja · Obnovljivi izvori

Folder predmeta Planiranje elektroenergetskih sistemaNa Drive-u: Predavanja · Slajdovi · Pitanja za usmeni ispit