Skip to content

Predavanje 07 i 08

Sigurnost elektroenergetskih sistema i rezerve proizvodnih kapaciteta

Predavanje prof. dr Andrije T. Sarića

Eksploatacija elektroenergetskih sistema, FTN

Sigurnost elektroenergetskih sistema

Sigurnost elektroenergetskih sistema (EES-a) označava sposobnost sistema da bez većih posledica podnese neke nepredviđene, ali verovatne poremećaje iz unapred formirane liste mogućih poremećaja. Znači, sigurnost se bavi problemima održavanja sistema u pogonu pri kvarovima.

Sistem je pogonski siguran ako smetnja sa liste: a) ne izaziva gubitak (oštećenje) elemenata, b) ne dovodi do poremećaja varijabli iznad dozvoljenih opsega (naponi, učestanost, tokovi snaga …).

Najčešće analizirani kvarovi: ispad generatora, ispad transformatora, ispad voda, nagla promena opterećenja …

Postoje tri moda sigurnosti:

  • Tranzijentna (mogu izazvati nestabilnost – gubitak sinhronizma). Brza kratkotrajna dinamika.
  • Dinamička (pojave duge dinamike – do 20 minuta). Vrši se analiza korena karakteristične jednačine.
  • Statička (posle iščezavanja prelazne pojave – stacionarno stanje).

Osnovne funkcije statičke sigurnosti su:

  1. Nadgledanje (nadzor, kontrola).
  2. Analiza poremećaja.
  3. Analiza korekcionih upravljačkih akcija i obezbeđenje sigurnosti.

Kontrolu sigurnosti sprovodi operator EES-a, praćenjem informacija o aktuelnom stanju (preko SCADA-e). Količina informacija koja ovde postoji je jako velika (koriste se procesni računari, sistem osnovnog (EMS) i dopunskog upravljanja).

Analiza mogućih poremećaja treba da omogući preduzimanje potrebnih mera za sprečavanje loših posledica. Cilj je sprečiti kaskadne poremećaje. Ovaj proces bi doveo do raspada sistema.

Obezbeđenje sigurnosti se sprovodi primenom korektivnih akcija koje treba sprovesti pre poremećaja, ako se utvrdi da je kvar kritičan.

Sve zajedno ovo je vrlo složen problem.

Važno je uočiti razliku između sigurnosti i pouzdanosti EES-a:

  • Sigurnost – determinističkog karaktera.
  • Pouzdanost – stohastički pokazatelj.

Sigurnost vrši ocenu performansi sistema na bazi tačno definisanih poremećaja, dok je drugi vezan za verovatnoću ispravnog rada komponenata i čitavog sistema.

Često se pouzdanost vezuje za planiranje EES-a, a za eksploataciju EES-a se vezuje sigurnost.

Radni režimi EES-a

Sistem može biti u više različitih režima. Prelaz je spontan (nekontrolisan), ili pod dejstvom poremećaja (kvarova), kao i usled korektivnih akcija.

Sistem karakterišu sledeći skupovi ograničenja:

  • Ograničenja potrošnje (jednačina bilansa aktivne i reaktivne snage):
G(x,u,z)=0.G(\boldsymbol{x}, \boldsymbol{u}, z) = 0.
  • Pogonska ograničenja (učestanost, naponi, tokovi snaga, PQ dijagrami):
H(x,u,z)≥0.H(\boldsymbol{x}, \boldsymbol{u}, z) \ge 0.
  • Ograničenja sigurnosti (prevodi potencijalno sistem u abnormalni režim):
S(x,u,z)≥0.S(\boldsymbol{x}, \boldsymbol{u}, z) \ge 0.

U prethodnim ograničenjima su:

  • x\boldsymbol{x} – vektor stanja (obično kompleksni naponi čvorova).
  • u\boldsymbol{u} – vektor upravljanja (obično snage generatora).
  • zz – vektor poremećaja (ispala snaga i slično).

Znaci u prethodnim ograničenjima definišu normalni sigurni režim.

Uobičajena praksa je da se sigurnost proverava za jednostruke ispade.

Mogući radni režimi EES-a dati su na donjoj slici:

Radni režimi EES-a i prelazi između njih. Normalni sigurni režim (G = 0, H ≥ 0, S ≥ 0) je sigurni nivo; sniženje nivoa rezerve i/ili povećana verovatnoća pojave poremećaja vodi ga u prethavarijski režim (G = 0, H ≥ 0, S ≱ 0), iz kog ga preventivno upravljanje vraća. Narušavanje pogonskih ograničenja vodi u poremećeni režim (G = 0, H ≱ 0), a razdvajanje sistema i/ili gubitak opterećenja u kritični režim (G ≠ 0, H ≱ 0). Korekciono (havarijsko) upravljanje vodi iz poremećenog u posthavarijski (G ≠ 0, H ≥ 0) ili u normalni režim, korekciono (obnavljajuće) iz kritičnog u posthavarijski, a korekciono (posthavarijsko) iz posthavarijskog u normalni ili prethavarijski. Prethavarijski i posthavarijski režim su nesigurni nivo, poremećeni i kritični ugroženi nivo; pune strelice su kontrolisani prelazi (pod dejstvom upravljanja), tačkaste spontani (pod dejstvom poremećaja).

Sa prethodne slike se vidi da postoje tri nivoa sigurnosti:

  • Sigurni nivo (Normalni sigurni režim).
  • Nesigurni nivo (Prethavarijski i Posthavarijski režim).
  • Ugroženi nivo (Poremećeni i Kritični režim).

Kriterijumi upravljanja i upravljačke akcije razlikuju se u pojedinim režimima:

  • Normalni sigurni

Kriterijum je ekonomičan rad sistema, uz zadovoljenje ograničenja sigurnosti.

Upravljačke akcije vezane za održavanje topologije, napona, učestanosti i razmene energije.

  • Prethavarijski

Kriterijum je održavanje sigurnosti.

Upravljačke akcije za povećanje sigurnosti, putem preventivnog delovanja na topologiju i raspodelu snaga, povećanje rezervi i slično.

  • Poremećen

Kriterijum je minimizacija vremena za otklanjanje narušavanja ograničenja.

Upravljačke akcije vezane za delovanje relejne zaštite i sl.

  • Kritičan

Kriterijum je minimizacija vremena dezintegracije sistema.

Upravljačke akcije su rasterećenje i razdvajanje sistema, resinhronizacija i uspostavljanje paralelnog rada.

  • Posthavarijski

Kriterijum je minimizacija neisporučene energije.

Upravljačke akcije su stavljanje u pogon i iskorišćenje hladne rezerve.

Vrste poremećaja

Svi poremećaji se dele u dve velike grupe:

  • Normalni poremećaji
  • Male promene opterećenja (deluju na spontane promene proizvodnje – primarna regulacija učestanosti)
  • Promene u strukturi vektora stanja i parametara
  • Promene u topologiji
  • Abnormalni poremećaji
  • Poremećaji u proizvodnji
  • Poremećaji u topologiji (promena statusa uključenosti)

Druga podela poremećaja je na interne i eksterne.

Listu poremećaja (ispada) obično čine:

  • Jednostruki ispadi (generator, transformator ili vod)
  • Verovatni dvostruki ispadi (generator-transformator, generator-vod)
  • Retki višestruki (kaskadni) ispadi

Analiza statičke sigurnosti

Jedan od mogućih pristupa je: proračunati tokove snaga za unapred zadatu listu ispada. Međutim, to zahteva dugo vreme za proračun. Zato se koristi DC model tokova snaga i faktori osetljivosti.

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 ii-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;
  • A={Aij}\boldsymbol{A} = \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 aℓia_\ell^i toka aktivne snage u ℓ\ell-toj grani usled promene injektiranja u i-tom čvoru je:

aℓi=ΔPℓΔPi=1Xℓg(Δθm−ΔθnΔPi),a_\ell^i = \frac{\Delta P_\ell}{\Delta P_i} = \frac{1}{X_\ell^g}\left(\frac{\Delta\theta_m - \Delta\theta_n}{\Delta P_i}\right),

gde je ℓ\ell-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 ℓ\ell-toj grani:

Pℓ(1)=Pℓ(0)+ΔPℓ=Pℓ(0)+aℓiΔPi.P_\ell^{(1)} = P_\ell^{(0)} + \Delta P_\ell = P_\ell^{(0)} + a_\ell^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:

Pℓ(1)=Pℓ(0)+∑i∈αpaℓiΔPi,P_\ell^{(1)} = P_\ell^{(0)} + \sum_{i \in \alpha_p}a_\ell^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:

aℓi={ΔPℓΔPi}=BℓgArTBr−1,\boxed{a_\ell^i = \left\{\frac{\Delta P_\ell}{\Delta P_i}\right\} = \boldsymbol{B}_\ell^g\boldsymbol{A_r}^{\mathrm{T}}\boldsymbol{B_r}^{-1}},

gde se može definisati matrica Aℓ=BℓgArTBr−1\boldsymbol{A}_\ell = \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)=(Br(0)−BmngekekT)−1Pr=θr(0)−(Br(0))−1ek(θm(0)−θn(0))1Bmng−ekT(Br(0))−1ek.\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} = \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^k\boldsymbol{P}_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:

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.

Sada je:

ΔPℓ=dℓkPk(0)+(aℓi+dℓkaki)ΔPi;ΔPℓ=dℓkPk(0)+(aℓi+dℓkaki)ΔPi=aℓi,kΔPi.\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)\Delta P_i = a_\ell^{i,k}\Delta P_i. \end{gathered}

Tok snage posle simultanog ispada je:

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

gde je, dakle, za simultani ispad generatora u i-tom čvoru (ΔPi=−Pgi\Delta P_i = -P_{gi}) 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.

Obezbeđenje sigurnosti

Glavno pitanje je šta treba preduzeti da se sistem povrati u siguran pogon?

Treba odgovoriti na pitanja šta TREBA, MOŽE i MORA da se uradi?

  • TREBA: Otkloniti potencijalne uzroke nesigurnosti akcijama preventivnog upravljanja.
  • MOŽE: Napraviti selekciju upravljačkih akcija i redosled primena (ekstremi su da se ne čini ništa, u nadi da se kritični kvar neće dogoditi, ili preduzeti sve potrebne akcije da se spreči nesigurnost).
  • MORA: Moraju se preduzeti potrebne akcije kada se kvar dogodi (korektivno upravljanje).

Moguće upravljačke akcije su:

1. Akcije na sredstva za proizvodnju ili potrošnju energije

  • Promene aktivnih i reaktivnih injektiranja generatora u cilju preraspodele tokova snaga.
  • Promene programa razmene energije sa susedima.
  • Rasterećenje potrošnje.

2. Akcije na prenosnu mrežu

  • Promena odnosa transformacije transformatora u cilju promene tokova reaktivnih snaga.
  • Promena topologije prenosne mreže.
  • Promene slike reaktivnih snaga u mreži (delovanje na regulatore pobude, nenominalne odnose transformacije, uključenje i isključenje statičkih izvora reaktivne snage i slično).

Rezerve proizvodnih kapaciteta

Rezerva je osnovni faktor koji određuje pouzdanost i sigurnost EES-a.

Ukupna rezerva (UR) proizvodnih kapaciteta je razlika između ukupne instalisane snage generatora (PgiP_{gi}) i stvarnog opterećenja EES-a (PpP_p):

UR=Pgi−Pp.UR = P_{gi} - P_p.

Ukupnu rezervu čine Remontna rezerva (ReR) i Raspoloživa rezerva proizvodnih kapaciteta (RPK):

UR=ReR+RPK.\mathrm{UR} = \mathrm{ReR} + \mathrm{RPK}.

Remontna rezerva (ReR) je onaj deo kapaciteta kojim se pokrivaju agregati u remontu, opravci i nezi.

Raspoloživa rezerva proizvodnih kapaciteta (RPK) služi za pokrivanje nepredviđenih ispada, a čine je Hladna rezerva (HR) i Operativna rezerva (OR):

RPK=HR+OR.\mathrm{RPK} = \mathrm{HR} + \mathrm{OR}.

Hladnu rezervu (HR) čine starije i skuplje termoelektrane. Mogu se staviti u pogon za 1 do 2 dana.

Operativna rezerva (OR) se neposredno koristi u eksploataciji, a čine je Rotirajuća operativna rezerva (RoR) i Nerotirajuća operativna rezerva (NRoR):

OR=RoR+NRoR.\mathrm{OR} = \mathrm{RoR} + \mathrm{NRoR}.

Rotirajuća operativna rezerva (RoR) je sinhronizovana na mrežu i može se startovati za 5 do 10 minuta, a čine je Regulaciona rezerva (RgR) i Rotirajuća havarijska rezerva (RoHR):

RoR=RgR+RoHR.\mathrm{RoR} = \mathrm{RgR} + \mathrm{RoHR}.

Regulaciona rezerva (RgR) služi za kompenzaciju grešaka između prognozirane i ostvarene potrošnje.

Rotirajuća havarijska rezerva (RoHR) služi za pokrivanje iznenadnih ispada.

Nerotirajuća operativna rezerva (NRoR) su kapaciteti koji nisu sinhronizovani na mrežu, ali se mogu brzo staviti u pogon, sinhronizovati i opteretiti. Čine je Hladna nerotirajuća operativna rezerva (HNRoR) i Topla nerotirajuća operativna rezerva (TNRoR):

NRoR=HNRoR+TNRoR.\mathrm{NRoR} = \mathrm{HNRoR} + \mathrm{TNRoR}.

Hladna nerotirajuća operativna rezerva (HNRoR) su brzo startujuće jedinice (uglavnom u hidroelektranama i gasnoturbinskim termoelektranama).

Topla nerotirajuća operativna rezerva (TNRoR) su agregati u parnim termoelektranama koje su van pogona, ali se održavaju u toplom stanju i spremni su za brzo pokretanje.

Na osnovu prethodno definisanih rezervi može se definisati i Ukupna havarijska rezerva (UHR) kao:

UHR=RoHR+NRoR=RoHR+HNRoR+TNRoR=OR−RgR.\boxed{\mathrm{UHR} = \mathrm{RoHR} + \mathrm{NRoR} = \mathrm{RoHR} + \mathrm{HNRoR} + \mathrm{TNRoR} = \mathrm{OR} - \mathrm{RgR}}.
Podela rezervi proizvodnih kapaciteta: ukupna rezerva UR deli se na remontnu rezervu ReR i ukupnu rezervu raspoloživih proizvodnih kapaciteta RPK; RPK na hladnu rezervu HR i operativnu rezervu OR; OR na rotirajuću RoR i nerotirajuću NRoR; RoR na rotirajuću regulacionu RgR i rotirajuću havarijsku RoHR; NRoR na hladnu nerotirajuću HNRoR (brzo startujuće jedinice) i toplu nerotirajuću TNRoR (TE održavane u toplom stanju); RoHR, HNRoR i TNRoR zajedno čine ukupnu havarijsku rezervu UHR.

Teme

Tokovi snaga · SCADA, EMS i DMS

Zadaci uz ovo predavanje

Predavanje u PDF-u, na Drive-u (otvara se u novom tabu)

Folder predmeta Eksploatacija elektroenergetskih sistemaNa Drive-u: Predavanja · Slajdovi · Statička sigurnost - MATLAB program · Pitanja za usmeni ispit