Skip to content

Predavanje 04

Gubici električne energije

Predavanje prof. dr Andrije T. Sarića

Distributivne i industrijske mreže, TF Čačak

Proračun gubitaka energije u elektroenergetskim sistemima (EES-ima), kao i distributivnim mrežama je vrlo važan, pošto ti gubici mogu biti desetak, pa i više procenata. Ovde nisu uključeni netehnički gubici (posledica grešaka mernih uređaja, nejednovremeno očitavanje brojila, subjektivne greške, neregistrovana potrošnja i drugi).

Gubici utiču na bilansiranje i eksploataciju EES-a.

Osnovni problemi pri proračunu gubitaka su:

  • Tačnost izračunatih gubitaka.
  • Razdvajanje tehničkih i komercijalnih gubitaka.
  • Raspodela gubitaka po elementima.
  • Određivanje mesta sa najvećim gubicima.
  • Analiza uzroka gubitaka.
  • Izbor mera za smanjenje gubitaka i drugi.

Generalno, gubici snage i energije dele se na:

  • Gubici koji zavise od napona.
  • Gubici koji zavise od struje.

Pošto se napon menja u uskim granicama oko nominalnog (±10\pm 10 %), onda se gubici energije koji zavise od napona mogu smatrati konstantnim. Izračunavaju se za nominalnu vrednost napona (jednostavan proračun). Približno su stalni (konstantni).

Određivanje dela gubitaka koji zavise od struje je vrlo složeno. Ovi gubici zavise od načina opterećivanja pojedinih elemenata (dijagrama opterećenja). Ovo je vrlo teško odrediti u praksi.

Zato se u praksi izbegava pristup preko dijagrama opterećenja, već se razvijaju uprošćeni postupci.

Postupci proračuna se dele na: determinističke i probabilističke.

Najčešće korišćene determinističke metode su: metode koje koriste ekvivalentno vreme trajanja maksimalnih gubitaka snage, metode bazirane na srednje-kvadratnoj struji, metode ekvivalentne otpornosti i druge, i one će biti predmet proučavanja u ovom kursu.

Pored prethodne podele na stalne i promenljive gubitke, prema elementima u kojima nastaju dele se na:

  • Gubitke u vodovima.
  • Gubitke u transformatorima.

Proračun gubitaka u jednom elementu

Posmatra se element između i-tog i j-tog čvora. Gubici aktivne snage u tom elementu su:

ΔPgij=3RijIij2.\Delta P_{gij} = 3 R_{ij} I_{ij}^2 .

Ako se uvedu pretpostavke konstantnog napona (Ui=Const.U_i = \text{Const.}) i faktora snage (cos⁡φij=Const.\cos\varphi_{ij} = \text{Const.}), uz Pij=3 UijIijcos⁡φijP_{ij} = \sqrt{3}\, U_{ij} I_{ij} \cos\varphi_{ij}, gubici snage su:

ΔPgij=3RijPij23cos⁡2φijUi2=cijPij2,\Delta P_{gij} = 3 R_{ij} \frac{P_{ij}^2}{3 \cos^2\varphi_{ij} U_i^2} = c_{ij} P_{ij}^2 ,

gde je cij=Rijcos⁡2φijUi2c_{ij} = \dfrac{R_{ij}}{\cos^2\varphi_{ij} U_i^2} konstanta koja vezuje gubitke i kvadrat toka aktivne snage.

Gubici energije u posmatranom elementu u vremenu T, koje je podeljeno na K intervala trajanja Δtk\Delta t_k (Pij,kP_{ij,k}), i konstantnog opterećenja su:

ΔWgij=cij∑k=1KPij,k2Δtk=cijT∑k=1KPij,k2℘(Δtk),\Delta W_{gij} = c_{ij} \sum_{k=1}^{K} P_{ij,k}^2 \Delta t_k = c_{ij} T \sum_{k=1}^{K} P_{ij,k}^2 \wp(\Delta t_k) ,

gde je ℘(Δtk)=ΔtkT\wp(\Delta t_k) = \dfrac{\Delta t_k}{T} verovatnoća pojavljivanja intervala Δtk\Delta t_k.

Opterećenje u k-tom vremenskom intervalu (Pij,kP_{ij,k}) je slučajna promenljiva koja podleže određenim pravilima iz teorije verovatnoće:

  • Matematičko očekivanje (srednja vrednost):
m{Pij}=1T∑k=1KPij,kΔtk.m\{P_{ij}\} = \frac{1}{T} \sum_{k=1}^{K} P_{ij,k} \Delta t_k .
  • Odstupanje:
ΔPij,k=Pij,k−m{Pij}.\Delta P_{ij,k} = P_{ij,k} - m\{P_{ij}\} .
  • Disperzija:
σ2{Pij}=1T∑k=1K(ΔPij,k)2Δtk.\sigma^2\{P_{ij}\} = \frac{1}{T} \sum_{k=1}^{K} \left( \Delta P_{ij,k} \right)^2 \Delta t_k .

To znači da je:

Pij,k=m{Pij}+ΔPij,k.P_{ij,k} = m\{P_{ij}\} + \Delta P_{ij,k} .

Ako se to zameni u prethodni izraz za gubitke aktivne energije dobija se:

ΔWgij=cij∑k=1KPij,k2Δtk=cijT∑k=1K(m{Pij}+ΔPij,k)2ΔtkT=cijT∑k=1K(m2{Pij}+2m{Pij}ΔPij,k+ΔPij,k2)ΔtkT\begin{aligned} \Delta W_{gij} &= c_{ij} \sum_{k=1}^{K} P_{ij,k}^2 \Delta t_k = c_{ij} T \sum_{k=1}^{K} \Big( m\{P_{ij}\} + \Delta P_{ij,k} \Big)^2 \frac{\Delta t_k}{T} \\ &= c_{ij} T \sum_{k=1}^{K} \Big( m^2\{P_{ij}\} + 2 m\{P_{ij}\} \Delta P_{ij,k} + \Delta P_{ij,k}^2 \Big) \frac{\Delta t_k}{T} \end{aligned}

Po definiciji matematičkog očekivanja i odstupanja je:

∑k=1KΔPij,kΔtk=0,\sum_{k=1}^{K} \Delta P_{ij,k} \Delta t_k = 0 ,

pa je:

ΔWgij=cijT(m2{Pij}+σ2{Pij}).\boxed{\Delta W_{gij} = c_{ij} T \Big( m^2\{P_{ij}\} + \sigma^2\{P_{ij}\} \Big)} .

Ovo je vrlo važna relacija, pošto omogućava proračun gubitaka energije u periodu T, ako je poznato matematičko očekivanje (srednja vrednost) i disperzija opterećenja u posmatranom periodu.

Gubici energije se mogu napisati i kao:

ΔWgij=cijT(Pijsr)2,\Delta W_{gij} = c_{ij} T \left( P_{ij}^{sr} \right)^2 ,

odakle je

(Pijsr)2=m{(Pijsr)2}=m2{Pij}+σ2{Pij}.\left( P_{ij}^{sr} \right)^2 = m\left\{ \left( P_{ij}^{sr} \right)^2 \right\} = m^2\{P_{ij}\} + \sigma^2\{P_{ij}\} .

Faktor opterećenja može se izraziti kao

m=PijsrPijmax=WijT⋅Pijmax=TMT.m = \frac{P_{ij}^{sr}}{P_{ij}^{max}} = \frac{W_{ij}}{T \cdot P_{ij}^{max}} = \frac{T^M}{T} .

Iz prethodnog izraza je

Pijsr=m{Pijsr}=PijmaxTMT.P_{ij}^{sr} = m\{P_{ij}^{sr}\} = P_{ij}^{max} \frac{T^M}{T} .

Drugi način proračuna gubitaka (rađen ranije) je preko maksimalne aktivne snage (PijmaxP_{ij}^{max}) i ekvivalentnog vremena trajanja maksimalnih gubitaka (τ\tau):

ΔWgij=ΔPgijmaxτg=cij(Pijmax)2τg.\Delta W_{gij} = \Delta P_{gij}^{max} \tau_g = c_{ij} \left( P_{ij}^{max} \right)^2 \tau_g .

Kombinujući prethodna dva izraza, dobija se:

ΔWgij=cijT(Pijmax)2(TMT)2+cijTσ2{Pij}=cij(Pijmax)2[T(TMT)2+Tσ2{Pij}(Pijmax)2]\begin{aligned} \Delta W_{gij} &= c_{ij} T \left( P_{ij}^{max} \right)^2 \left( \dfrac{T^M}{T} \right)^2 + c_{ij} T \sigma^2\{P_{ij}\} \\ &= c_{ij} \left( P_{ij}^{max} \right)^2 \left[ T \left( \dfrac{T^M}{T} \right)^2 + T \frac{\sigma^2\{P_{ij}\}}{\left( P_{ij}^{max} \right)^2} \right] \end{aligned}

U prethodnom izrazu je TMT^M ekvivalentno vreme trajanja maksimalne snage.

Ekvivalentno vreme trajanja maksimalnih gubitaka poređenjem prethodna dva izraza je:

τg=T(TMT)2+Tσ2{Pij}(Pijmax)2,\boxed{\tau_g = T \left( \dfrac{T^M}{T} \right)^2 + T \frac{\sigma^2\{P_{ij}\}}{\left( P_{ij}^{max} \right)^2}} ,

odnosno:

τg=T(m2+σ2{Pij}(Pijmax)2)=Tμ,\boxed{\tau_g = T \left( m^2 + \frac{\sigma^2\{P_{ij}\}}{\left( P_{ij}^{max} \right)^2} \right) = T \mu} ,

gde je m faktor opterećenja, a μ\mu faktor gubitaka.

U prethodnim izrazima član σ2{Pij}\sigma^2\{P_{ij}\} uvažava promenljivost dijagrama opterećenja.

Determinističke metode proračuna godišnjih gubitaka energije

Metod integracije

Ovaj metod zahteva poznavanje dijagrama opterećenja, ili krive trajanja opterećenja.

Proračun se vrši preko izraza:

ΔWgij=3Rij∫0TIij2(t) dt≈cij∫0TPij2(t) dt.\Delta W_{gij} = 3 R_{ij} \int_0^T I_{ij}^2(t)\, \mathrm{d}t \approx c_{ij} \int_0^T P_{ij}^2(t)\, \mathrm{d}t .

Drugi način je diskretizacija dijagrama opterećenja (na primer u koracima Δtk=1\Delta t_k = 1 h), čime se prethodni integral svodi na sumu:

ΔWgij=cij∑k=1KPij,k2⋅Δtk.\Delta W_{gij} = c_{ij} \sum_{k=1}^{K} P_{ij,k}^2 \cdot \Delta t_k .

Metod srednje-kvadratne struje

Polazi se od poznatog integrala:

ΔWgij=3Rij∫0TIij2(t) dt=3Rij(Iijsk)2T,\Delta W_{gij} = 3 R_{ij} \int_0^T I_{ij}^2(t)\, \mathrm{d}t = 3 R_{ij} \left( I_{ij}^{sk} \right)^2 T ,

gde je srednje-kvadratna struja za period T:

Iijsk=1T∫0TIij2(t) dt=∑k=1KIij,k2T.I_{ij}^{sk} = \sqrt{\frac{1}{T} \int_0^T I_{ij}^2(t)\, \mathrm{d}t} = \sqrt{\frac{\sum\limits_{k=1}^{K} I_{ij,k}^2}{T}} .

Pošto je dominantno opterećenje aktivno, onda se kao i ranije za pretpostavke Vi=Const.V_i = \text{Const.} i cos⁡φij=Const.\cos\varphi_{ij} = \text{Const.} dobija:

ΔWgij=cij∫0TPij2(t) dt=cij(Pijsk)2T.\Delta W_{gij} = c_{ij} \int_0^T P_{ij}^2(t)\, \mathrm{d}t = c_{ij} \left( P_{ij}^{sk} \right)^2 T .

Metod srednje struje

Prethodni postupak zahteva poznavanje dijagrama opterećenja, što je teško u praksi. Zato se koristi veza srednje-kvadratne i srednje struje:

Iijsk=k0Iijsr,I_{ij}^{sk} = k_0 I_{ij}^{sr} ,

koja pod pretpostavkom Vi=Const.V_i = \text{Const.} i cos⁡φij=Const.\cos\varphi_{ij} = \text{Const.} važi i za aktivne snage:

Pijsk=k0Pijsr,P_{ij}^{sk} = k_0 P_{ij}^{sr} ,

gde je k0k_0 koeficijent oblika, koji povezuje srednje-kvadratnu i srednju snagu (struju), koji se obično definiše kao:

k0=m{Pij2}m2{Pij}.k_0 = \sqrt{\frac{m\{P_{ij}^2\}}{m^2\{P_{ij}\}}} .

Uvažavajući prethodno, godišnji gubici energije su:

ΔWgij=cijk02(Pijsr)2T.\boxed{\Delta W_{gij} = c_{ij} k_0^2 \left( P_{ij}^{sr} \right)^2 T} .

Koeficijent oblika k0k_0 izračunava se na osnovu dijagrama opterećenja. Međutim, k0k_0 se menja u relativno uskom opsegu. Minimalna vrednost k0k_0 je 1, kada se opterećenje ne menja tokom vremena. Za maksimalnu vrednost koeficijent k0k_0 treba sprovesti detaljniju analizu. Polazeći od napred izvedenog izraza:

m{Pij2}=m2{Pij}+σ2{Pij},m\{P_{ij}^2\} = m^2\{P_{ij}\} + \sigma^2\{P_{ij}\} ,

i zamenom u prethodni izraz za k0k_0 dobija se:

k0=1+σ2{Pij}m2{Pij}.\boxed{k_0 = \sqrt{1 + \frac{\sigma^2\{P_{ij}\}}{m^2\{P_{ij}\}}}} .

Iz prethodnog izraza zaključuje se da maksimalna vrednost koeficijenta oblika k0k_0 odgovara maksimalnom odnosu disperzije i matematičkog očekivanja (srednje vrednosti).

Metod ekvivalentnog vremena trajanja gubitaka

Najčešće korišćeni metod za procenu gubitaka električne energije, pošto je jednostavan za primenu. Za njega nije neophodno znati dijagram opterećenja, već samo pokazatelj dijagrama opterećenja koji se naziva ekvivalentno vreme trajanja maksimalnih gubitaka (izvedeno ranije):

τg=∫0TPij2(t) dt(Pijmax)2.\boxed{\tau_g = \frac{\int_0^T P_{ij}^2(t)\, \mathrm{d}t}{\left( P_{ij}^{max} \right)^2}} .

Međutim, problem nastaje što nije lako dati vezu između τ\tau i faktora opterećenja (takođe ranije izveden):

m=PijsrPijmax=TMT.m = \frac{P_{ij}^{sr}}{P_{ij}^{max}} = \frac{T^M}{T} .

Najčešće korišćena veza je:

τg=0,17TM+0,83(TM)2T=0,17mT+0,83m2T.\boxed{\tau_g = 0{,}17 T^M + 0{,}83 \frac{\left( T^M \right)^2}{T} = 0{,}17 m T + 0{,}83 m^2 T} .

Proračun gubitaka energije u distributivnim mrežama

Najčešće se koristi metod ekvivalentne otpornosti, zbog brzine proračuna i malog broja potrebnih informacija. Ideja je da se realna mreža zameni ekvivalentnom koja ima jednake gubitke, kao i realna mreža.

U izvođenju izraza za ekvivalentnu otpornost posmatra se distributivna mreža sa N grana, od čega je NtN_t transformatora i NvN_v vodova.

Gubici u gvožđu malo se menjaju sa promenom opterećenja, tako da su ukupni gubici u gvožđu čitave mreže:

ΔPFe=∑ℓ=1NtΔPFeℓ.\Delta P_{Fe} = \sum_{\ell=1}^{N_t} \Delta P_{Fe\ell} .

Gubici snage u vodovima i bakru transformatora zavise od opterećenja i čine varijabilne (promenljive) gubitke. Mogu se napisati kao:

ΔPg var=∑ℓ=1NvΔPgℓ+∑ℓ=1NtΔPCuℓ=∑ℓ=1Nv3RℓIℓ2+∑ℓ=1NtPCunℓ(IℓIℓn)2.\Delta P_{g\,var} = \sum_{\ell=1}^{N_v} \Delta P_{g\ell} + \sum_{\ell=1}^{N_t} \Delta P_{Cu\ell} = \sum_{\ell=1}^{N_v} 3 R_\ell I_\ell^2 + \sum_{\ell=1}^{N_t} P_{Cun\ell} \left( \frac{I_\ell}{I_{\ell n}} \right)^2 .

Preko ℓ\ell-te grane vrši se snabdevanje energijom nℓn_\ell čvorova, tako da je struja ℓ\ell-te grane:

Iℓ=∑i=1nℓIi.I_\ell = \sum_{i=1}^{n_\ell} I_i .

Uvažavanjem koeficijenata učešća KiK_i pojedinih čvorova u struji napojne deonice, struja ℓ\ell-te grane može se izračunati kao:

Iℓ=Ind∑i=1nℓKi,I_\ell = I_{nd} \sum_{i=1}^{n_\ell} K_i ,

gde su:

  • IndI_{nd} – struja napojne deonice;
  • Ki=IiIndK_i = \dfrac{I_i}{I_{nd}} – faktor učešća.

Sada se gubici mogu prikazati kao:

ΔPg var=Ind2[∑ℓ=1Nv3Rℓ(∑i=1nℓKi)2+∑ℓ=1NtPCunℓIℓn2(∑i=1nℓKi)2],\Delta P_{g\,var} = I_{nd}^2 \left[ \sum_{\ell=1}^{N_v} 3 R_\ell \left( \sum_{i=1}^{n_\ell} K_i \right)^2 + \sum_{\ell=1}^{N_t} \frac{P_{Cun\ell}}{I_{\ell n}^2} \left( \sum_{i=1}^{n_\ell} K_i \right)^2 \right] ,

odnosno:

ΔPg var=Ind2Rek,\Delta P_{g\,var} = I_{nd}^2 R_{ek} ,

gde ekvivalentnu otpornost RekR_{ek} predstavlja izraz u zagradi u prethodnom izrazu.

Uvođenje prethodne ekvivalentne otpornosti počiva na činjenici da se ni jedna veličina u zagradi ne menja sa promenom opterećenja u mreži.

Proračun gubitaka izloženim metodom u distributivnoj mreži sprovodi se kroz sledeće korake:

1. Za jedan k-ti (karakterističan) režim mreže vrši se proračun promenljivih gubitaka ΔPg vark\Delta P_{g\,var}^k, gubitaka u gvožđu transformatora ΔPFe\Delta P_{Fe} (ne zavise od režima) i struje napojne deonice IndkI_{nd}^k.

2. Gubici u gvožđu transformatora za proizvoljno opterećenu mrežu su konstantni.

3. Proračun ekvivalentne otpornosti na osnovu karakterističnog režima:

Rek=ΔPg vark(Indk)2.R_{ek} = \frac{\Delta P_{g\,var}^k}{\left( I_{nd}^k \right)^2} .

Ova otpornost ostaje konstantna za sve režime.

4. Promenljivi gubici snage za proizvoljno opterećenje računaju se iz izraza:

ΔPg varn=Rek(Indn)2.\Delta P_{g\,var}^n = R_{ek} \left( I_{nd}^n \right)^2 .

Važno je pitanje kako izabrati karakterističan režim za proračun ekvivalentne otpornosti? Najčešće je to režim maksimalnog opterećenja, pošto se onda ide na stranu sigurnosti.

U praksi je dosta teško raspolagati i sa faktorima učešća deonica (KiK_i) u struji napojne deonice. U tom slučaju uvodi se još jedno uprošćenje: faktor učešća KiK_i je srazmeran odnosu naznačene snage transformatora u tom čvoru i ukupne snage svih transformatora koje taj vod napaja.

Tada se ekvivalentna otpornost izračunava kao:

Rek=∑ℓ=1Nv3Rℓ(SℓSΣn)2+∑ℓ=1NtPCunℓVjn2SΣn2.R_{ek} = \sum_{\ell=1}^{N_v} 3 R_\ell \left( \frac{S_\ell}{S_{\Sigma n}} \right)^2 + \sum_{\ell=1}^{N_t} \frac{P_{Cun\ell} V_{jn}^2}{S_{\Sigma n}^2} .

Da bi se primenio ovaj pristup neophodno je znati samo podatke o elementima i opterećenju napojne stanice. Dakle, postupak proračuna je jednostavan.

Međutim, problem sa ovim pristupom je što se pretpostavlja da je opterećenje transformatora srazmerno njegovoj nominalnoj snazi. Ovo kod novougrađenih transformatora nije uvek tačno (oni su najčešće malo opterećeni, kako bi se ostavila mogućnost porasta opterećenja u budućnosti), pa tada postupak može znatno da greši.

Teme

Gubici · Distributivne mreže