20.01.2015 Views

Statistička metoda za procenu proizvodnje električne ... - Infoteh

Statistička metoda za procenu proizvodnje električne ... - Infoteh

Statistička metoda za procenu proizvodnje električne ... - Infoteh

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

INFOTEH-JAHORINA Vol. 11, March 2012.<br />

Statistička <strong>metoda</strong> <strong>za</strong> <strong>procenu</strong> <strong>proizvodnje</strong> električne<br />

energije iz fotonaponskog sistema<br />

Jovan Mikulović<br />

Elektrotehnički fakultet<br />

Univerzitet u Beogradu<br />

Beograd, Srbija<br />

mikulovic@etf.rs<br />

Željko Đurišić<br />

Elektrotehnički fakultet<br />

Univerzitet u Beogradu<br />

Beograd, Srbija<br />

djurisic@etf.rs<br />

Sadržaj—U ovom radu je prika<strong>za</strong>na statistička <strong>metoda</strong> <strong>za</strong><br />

<strong>procenu</strong> sezonske i godišnje <strong>proizvodnje</strong> električne energije iz<br />

fotonaponskog (PV) sistema. Na osnovu raspoloživih mernih<br />

podataka o srednjoj desetominutnoj horizontalnoj solarnoj<br />

iradijaciji i temperaturi ambijenta u toku jedne godine, dobijeni<br />

su dnevni dijagrami <strong>proizvodnje</strong> fotonaponskog sistema.<br />

Statističkom obradom dnevnih dijagrama <strong>proizvodnje</strong> i<strong>za</strong>brani<br />

su karakteristični sezonski dani, odnosno karakteristični sezonski<br />

dijagrami snage <strong>proizvodnje</strong> fotonaponskog sistema koji se mogu<br />

smatrati reprezentativnim <strong>za</strong> <strong>procenu</strong> proizvedene energije u<br />

posmatranom periodu vremena. Korišćenjem navedene metode<br />

izvršena je procena sezonske i godišnje <strong>proizvodnje</strong><br />

fotonaponskog sistema instalisane snage 495 kW.<br />

Ključne reči - fotonaponski sistem; dijagram <strong>proizvodnje</strong>;<br />

karakterisični dijagram; proizvedena energija; faktor kapaciteta<br />

I. UVOD<br />

Fotonaponski (PV) sistem predstavlja integrisan skup<br />

fotonaponskih modula i drugih komponenti, projektovan tako<br />

da primarnu solarnu energiju direktno pretvara u električnu<br />

energiju kojom se osigurava rad određenog broja potrošača<br />

jednosmerne (DC) i/ili naizmenične (AC) struje. Kod<br />

fotonaponskih sistema pove<strong>za</strong>nih na distributivnu mrežu,<br />

jednosmerna struja dobijena u solarnim modulima se pretvara<br />

u naizmeničnu struju pomoću invertora koji je pove<strong>za</strong>n na<br />

mrežu tako da se, osim napajanja potrošača, vrši i razmena<br />

energije sa mrežom [1].<br />

Projektovanje fotonaponskih sistema se najčešće vrši na<br />

osnovu njihove godišnje proizvedene energije, koja je takođe<br />

dobar parametar <strong>za</strong> praćenje dugoročnih performansi<br />

fotonaponskog sistema. Za predviđanje godišnje proizvedene<br />

energije fotonaponskog sistema neophodni su pouzdani<br />

modeli i metode s obzirom na stohastičku prirodu solarnog<br />

zračenja i na veliki broj uticajnih faktora (ambijentalni uslovi i<br />

performanse sistema) [2]. U literaturi su dati modeli <strong>za</strong><br />

predikciju snage i energije solarnog zračenja po jedinici<br />

površine (iradijacije i insolacije) [3-5]. Statističko<br />

modelovanje energije proizvedene fotonaponskim sistemom je<br />

dato u [6-8].<br />

U ovom radu je prika<strong>za</strong>na statistička <strong>metoda</strong> <strong>za</strong> <strong>procenu</strong><br />

sezonske i godišnje proizvedene energije fotonaponskog<br />

sistema koja se <strong>za</strong>sniva na karakterističnim sezonskim<br />

dijagramima snaga <strong>proizvodnje</strong>. Verifikacija metode je<br />

izvršena na primeru sezonske i godišnje <strong>proizvodnje</strong><br />

fotonaponskog sistema instalisane snage 495 kW.<br />

II. ANALIZA PROIZVODNJE FOTONAPONSKOG<br />

SISTEMA NA OSNOVU MERNIH PODATAKA<br />

Da bi se projektovali mrežno pove<strong>za</strong>ni PV sistemi moraju<br />

se dobro poznavati resursi solarne energije, karakteristike svih<br />

elemenata sistema i ambijentalni uslovi. Proračun resursa<br />

solarne energije se vrši na osnovu merenja i na osnovu<br />

proračuna insolacije na površini na kojoj se planira<br />

postavljanje panela. Polazna tačka <strong>za</strong> određivanje performansi<br />

sistema <strong>za</strong> poznatu ulaznu snagu zračenja je naznačena DC<br />

izlazna snaga fotonaponskog modula pri standardnim<br />

uslovima (iradijacija “jedno sunce“ tj. 1 kWh/m 2 , koeficijent<br />

vazdušne mase AM 1.5, tempeatura ćelije 25 0 C, moduli<br />

potpuno čisti). U realnim eksploatacionim uslovima snaga<br />

koju PV sistem predaje mreži P AC je manja od DC izlazne<br />

snage modula pri standardnim uslovima P DC (STC) zbog<br />

gubitaka odnosno efikasnosti konverzije:<br />

P = P ⋅η ⋅η ⋅η ⋅ η<br />

(1)<br />

AC DC (STC) Z N T inv<br />

gde η Z , η N i η T definišu smanjenje efikasnosti panela usled<br />

<strong>za</strong>prljanja, neuparenosti i povećanja temperature modula u<br />

odnosu na standardnu vrednost 25 0 C, a η inv definiše efikasnost<br />

invertora.<br />

Na efikasnost konverzije bitno utiču efikasnost invertora,<br />

smanjenje efikasnosti usled <strong>za</strong>prljanja modula neuparenost<br />

modula i razlike u ambijentalnim uslovima u odnosu na<br />

standardne. Pri iradijaciji od “jednog sunca“ rezultat ovih<br />

gubitaka može smanjiti izlaznu snagu <strong>za</strong> 20–40 % u odnosu na<br />

P DC (STC) . Efikasnost invertora varira u <strong>za</strong>visnosti od<br />

opterećenja. Mrežno pove<strong>za</strong>ni invertori imaju efikasnost preko<br />

90 %, osim pri vrlo malim opterećenjima. Neuparenost<br />

modula i<strong>za</strong>ziva smanjenje izlazne snage paralelno ve<strong>za</strong>nih<br />

modula zbog toga što strujno-naponske karakteristike modula<br />

nisu identične. Da bi efikasnost sistema bila što veća treba<br />

težiti da moduli koji formiraju panele budu što približnijih<br />

karakteristika.<br />

Važan faktor koji utiče na efikasnost PV sistema je i<br />

temperatura fotonaponskih modula. Da bi se mogla odrediti<br />

- 187 -


efikasnost modula pri različitim ambijentalnim uslovima<br />

potrebno je proračunati temperaturu modula. Na temperaturu<br />

modula dominantno utiče snaga zračenja Sunca i uslovi<br />

hlađenja, odnosno vetar. Za svaki modul proizvođač definiše<br />

temperaturu pri nominalnim uslovima eksploatacije (NOCT –<br />

Nominal Operation Cell Temperature). NOCT je temperatura<br />

modula pri ambijentalnoj temperaturi 20 0 C, solarnoj iradijaciji<br />

800 W/m 2 brzini vetra 1 m/s. Na osnovu podatka NOCT se<br />

može proceniti temperatura ćelije (modula ili panela):<br />

⎛ NOCT − 20 ⎞<br />

T = T + ⎜<br />

⎟⋅I<br />

⎝ 0.8 ⎠<br />

PV amb PV<br />

gde su: T PV temperatura solarnih modula, T amb temperatura<br />

ambijenta, I PV solarna iradijacija na površini modula.<br />

S obzirom da je pad efikasnosti fotonaponskih ćelija usled<br />

povećanja temperature solarne ćelije iznad standardne<br />

vrednosti (25 0 C) ΔP=-0,5 %/ 0 C, onda je snaga na DC<br />

priključcima sistema:<br />

( 1 0,005( ))<br />

DC(PTC) DC (STC) cell amb<br />

(2)<br />

P = P − T − T<br />

(3)<br />

Kada su poznati podaci <strong>za</strong> dnevne, mesečne ili godišnje<br />

prosečne iradijacije (kWh/m 2 ) na mestu analiziranog PV<br />

sistema, tada se može relativno jednostavno proceniti<br />

proizvodnja električne energije. Ako je proračunata dnevna<br />

insolacija <strong>za</strong> neku lokaciju npr. 5.6 kWh/m 2 dan, u pogledu<br />

<strong>proizvodnje</strong> električne energije PV sistema instaliranog na toj<br />

lokaciji može se posmatrati kao da 5.6 h/dan postoji<br />

iradijacija od “jednog sunca“ (1 kWh/m 2 ), ili 5,6 h “zenita<br />

sunca“. Energija na AC priključcima PV sistema pri<br />

standardnoj iradijaciji od “jednog sunca” se može izračunati<br />

tako što se izračunata iradijacija pomnoži sa brojem sati<br />

zenitnog sunca i na taj način se dobija ukupna očekivana<br />

energija <strong>proizvodnje</strong>:<br />

W I ⋅A ⋅ η (4)<br />

2 2<br />

[kWh/dan] =<br />

av<br />

[kWh/m /dan] [m ]<br />

gde su: W energija koju sistem proizvede u toku jednog dana,<br />

I av srednja dnevna insolacija na panel, A površina PV panela,<br />

η av prosečna efikasnost sistema u toku dana. Efikasnost<br />

sistema η av , pored smanjenja efikasnosti panela usled<br />

<strong>za</strong>prljanja, neuparenosti i povećanja temperature modula,<br />

uključuje efikasnost invertora η inv i efikasnost samih<br />

modula η m .<br />

Pri iradijaciji od “jednog sunca“, AC snaga PV sistema je:<br />

P ⋅A ⋅ η (5)<br />

2 2<br />

AC<br />

[kW] = 1 [kW/m ] [m ]<br />

1-sunce<br />

gde je η 1-sunce efikasnost sistema pri iradijaciji od “jednog<br />

sunca“. Kombinujući prethodne dve jednačine dobija se:<br />

W [kWh/dan] = P [kW]<br />

2<br />

I [kWh/m /dan] av<br />

AC<br />

⋅ ⋅ (6)<br />

2<br />

1 [kW/m ] η1-sunce<br />

η<br />

av<br />

Ako se pretpostavi da je prosečna efikasnost sistema po<br />

danu jednaka efikasnosti kada je iradijacija “jedno sunce“,<br />

tada je električna energija koja se predaje distributivnom<br />

sistemu:<br />

W I h<br />

2<br />

[kWh/dan]=<br />

av[kWh/m /dan] ⋅<br />

z[broj sati "zenita sunca"]<br />

Glavna pretpostavka u prethodnoj jednačini je da<br />

efikasnost sistema ostaje konstantna tokom dana. Ova<br />

pretpostavka je realna ako PV sistem ima uređaj <strong>za</strong><br />

optimi<strong>za</strong>ciju radne tačke, odnosno MPPT – Maximal Power<br />

Point Tracking, koji pri svim uslovima obezbeđuje da sistem<br />

radi u tački maksimalne snage. Pošto je snaga u tački<br />

maksimalne snage direktno proporcionalna iradijaciji,<br />

efikasnost sistema je obično konstantna. Uticaj promene<br />

temperature i vazdušne mase takođe imaju udela u efikasnosti<br />

konverzije ali su greške usled njihove dnevne varijacije<br />

relativno male. Efikasnost bi bila <strong>za</strong> nijansu iznad proseka<br />

ujutru, kada je hladnije i kada je manja iradijacija ali se ovi<br />

uticaji usrednjavanjem u dobroj meri kompenzuju. S obzirom<br />

da se maksimalna insolacija po pravilu javlja u periodima dana<br />

kada je i temperatura maksimalna, u anali<strong>za</strong>ma mesečne<br />

<strong>proizvodnje</strong> električne energije nekog PV sistema treba<br />

koristiti srednju vrednost maksimalnih dnevnih temperatura.<br />

Uobičajen i jednostavan način <strong>za</strong> prikazivanje efikasnosti<br />

<strong>proizvodnje</strong> energije bilo kog sistema <strong>za</strong> proizvodnju<br />

električne energije je prikazivanjem njegove AC energije i<br />

faktora kapaciteta. Faktor kapaciteta se obično definiše na<br />

godišnjem nivou, ali se može definisati i na dnevnom,<br />

nedeljnom i mesečnom nivou. Sledeća jednačina daje vezu<br />

između dnevne proizvedene električne energije i faktora<br />

kapaciteta na dnevnom nivou (CF – Capacity Factor):<br />

AC<br />

(7)<br />

W [kWh/dan] = P [kW] ⋅CF⋅ 24 [h] (8)<br />

Kombinovanjem prethodnih jednačina može se<br />

jednostavno definisati faktor kapaciteta <strong>za</strong> mrežno pove<strong>za</strong>ne<br />

PV sisteme:<br />

h [broj sati "zenita sunca"]<br />

24 h<br />

z<br />

CF = (9)<br />

III. STATISTIČKA METODA ZA ANALIZU<br />

DIJAGRAMA PROIZVODNJE FOTONAPONSKOG<br />

SISTEMA<br />

S obzirom na stohastičku prirodu solarnog zračenja na<br />

površini zemlje, anali<strong>za</strong> <strong>proizvodnje</strong> fotonaponskog panala se<br />

može izvršiti statističkim <strong>metoda</strong>ma. Statistička anali<strong>za</strong><br />

<strong>proizvodnje</strong> fotonaponskog panala se može izvršiti na osnovu<br />

karakterističnih sezonskih dijagrama snage <strong>proizvodnje</strong><br />

panela.<br />

Dnevni dijagrami snaga <strong>proizvodnje</strong> fotonaponskog sistema<br />

se mogu dobiti na osnovu izmerenih vrednosti iradijacije na<br />

- 188 -


fotonaponski panel I PV , površine panela A i ukupne efikasnosti<br />

fotonaponskog sistema η:<br />

P I ⋅A ⋅ η (10)<br />

2 2<br />

AC<br />

[kW] =<br />

PV<br />

[kW/m ] [m ]<br />

Uvažavajući smanjenje efikasnosti η Z , η N i η T usled<br />

<strong>za</strong>prljanja, neuparenosti i <strong>za</strong>grevanja modula, kao i efikasnost<br />

invertora η inv i efikasnost samih modula η m , efikasnost<br />

fotonaponskog sistema η može da se izrazi kao:<br />

η = η ⋅η ⋅η ⋅η ⋅ η<br />

(11)<br />

Z N T inv m<br />

Statističkom obradom dnevnih dijagrama snaga <strong>proizvodnje</strong><br />

može se odrediti karakterističan sezonski dan, odnosno<br />

karakteristični sezonski dijagrami snage <strong>proizvodnje</strong> sistema.<br />

Za sve dane u razmatranom periodu, iz dnevnih dijagrama<br />

snaga <strong>proizvodnje</strong> određuju se karakteristične veličine,<br />

odnosno karakteristični poka<strong>za</strong>telji. Karakteristične veličine<br />

su maksimalna dnevna snaga <strong>proizvodnje</strong> P max , srednja dnevna<br />

snaga <strong>proizvodnje</strong> P sr i srednja dnevna snaga <strong>proizvodnje</strong> na<br />

kvadrat i određuju na sledeći način [10]:<br />

Pmax = max<br />

i( Pi)<br />

i=1...m; (12)<br />

P<br />

P<br />

sr<br />

q<br />

1 m<br />

= ∑ P<br />

(13)<br />

m i = 1<br />

1 m<br />

m i = 1<br />

i<br />

= ∑ P<br />

(14)<br />

gde m predstavlja broj merenja insolacije u toku dana,<br />

odnosno broj vrednosti koje definišu dnevni dijagram snage<br />

<strong>proizvodnje</strong> fotonaponskog sistema.<br />

Maksimalne snage određuju termičko opterećenje vodova i<br />

transformatora. Srednje dnevne snage bitne su <strong>za</strong> ocenu<br />

padova napona, faktora snage i stepena kompen<strong>za</strong>cije<br />

reaktivnih snaga [10]. Srednje snage na kvadrat u korelaciji su<br />

sa gubicima snage i energije u distributivnoj mreži na koju se<br />

priključuje fotonaponski sistem.<br />

Za sve dane određuju se srednje vrednosti karakterističnih<br />

poka<strong>za</strong>telja (matematičko očekivanje) i verovatne donje i<br />

gornje granične vrednosti ovih poka<strong>za</strong>telja (granice intervala<br />

poverenja). Neka su X k vrednosti poka<strong>za</strong>telja X <strong>za</strong> dan "k".<br />

Matematičko očekivanje poka<strong>za</strong>telja X jednako je:<br />

2<br />

i<br />

1 n X<br />

k<br />

n k = 1<br />

X = ∑ (15)<br />

Za određivanje intervala poverenja korišćena je normalna<br />

raspodela. Parametar σ n-1 je srednje kvadratno odstupanje:<br />

n<br />

σ<br />

n − 1<br />

n 1 k = 1<br />

( X X ) 2<br />

k<br />

= 1<br />

∑ − −<br />

(17)<br />

Karakteristični dnevni dijagram snage <strong>proizvodnje</strong> je<br />

dnevni dijagram čiji karakteristični poka<strong>za</strong>telji <strong>za</strong>dovoljavaju<br />

intervale poverenja (14) i najmanje odstupaju od svojih<br />

matematičkih očekivanja (13). Karakteristični dijagram se<br />

može smatrati reprezentativnim dijagramom u posmatranom<br />

periodu. Za analizu <strong>proizvodnje</strong> fotonaponskog sistema<br />

pogodno je analizirati četiri perioda u toku godine: prolećna<br />

sezona (mart, april, maj), letnja sezona (jun, jul, avgust),<br />

jesenja sezona (septembar, oktobar, novembar) i zimska<br />

sezona (decembar, januar, februar). Sezone su određene na<br />

osnovu mesečnih proizvodnji fotonaponskog sistema (zimsku<br />

sezonu čine tri meseca sa najmanjom proizvodnjom).<br />

IV. PROCENA SEZONSKE I GODIŠNJE<br />

PROIZVODNJE FOTONAPONSKOG SISTEMA<br />

Na osnovu raspoloživih mernih podataka o srednjoj<br />

desetominutnoj horizontalnoj iradijaciji [W/m 2 ] u toku 2009.<br />

godine <strong>za</strong> Beograd [9], proračunate su ukupne godišnje<br />

insolacije <strong>za</strong> različite nagibne i azimutne uglove PV modula.<br />

Izračunate insolacije su korigovane uvažavanjem uticaja<br />

temperature ambijenta prema izrazu (2) tako da su dobijene<br />

raspoložive energije modula. Najveća vrednost raspoložive<br />

energije iznosi 1316.1 kWh/m 2 i dobijena je pri nagibnom<br />

uglu modula od 32 0 i azimutnom uglu modula od -1 0 . Srednja<br />

dnevna raspoloživa energija pri nagibnom uglu modula od 32 0<br />

i azimutnom uglu modula od -1 0 iznosi 3.61 kWh/m 2 .<br />

Na osnovu desetominutnih vrednosti <strong>za</strong> horizontalnu<br />

iradijaciju izračunate su srednje dnevne insolacije na<br />

optimalno orijentisane fotonaponske module <strong>za</strong> svaki dan u<br />

toku godine i prika<strong>za</strong>ne su na slici 1. Vrednosti srednjih<br />

dnevnih insolacija su istovremeno i vrednosti broja sati<br />

“zenita sunca”. Srednja vrednost srednjih dnevnih insolacija<br />

na module u toku godine iznosi 4.01 kWh/m 2 .<br />

I av<br />

[ kWh / m 2 ]<br />

8<br />

7<br />

6<br />

5<br />

4<br />

3<br />

gde je n statistički uzorak, jednak broju dana u posmatranom<br />

periodu vremena.<br />

Sa verovatnoćom 0.95 veličina X se nalazi u intervalu:<br />

σ<br />

X X X<br />

n<br />

n−1 n−1<br />

− 1.96 < < + 1.96<br />

(16)<br />

σ<br />

n<br />

2<br />

1<br />

0<br />

0 50 100 150 200 250 300 350<br />

dan u godini<br />

Slika 1. Srednje dnevne insolacije (broj sati “zenita sunca”)<br />

- 189 -


Instalisana snaga fotonaponskog panela jednaka je<br />

proizvodu ukupnog broja solarnih modula koji čine panel i<br />

nominalne snage jednog modula. Ova DC snaga se ima pri<br />

standardnim test uslovima (STC). Za formiranje<br />

fotonaponskog panela i<strong>za</strong>brani su fotonaponski moduli<br />

proizvodjača SHARP model ND-220E1F, naznačene snage<br />

220 W, površine 1,642 m 2 i efikasnosti 13.4 %. Usvojeno je da<br />

se panel sastoji iz 2250 modula da bi istalisana snaga panela<br />

bila približno 500 kW. Tačna vrednost instalisane snage<br />

panela je P DC(STC) = 495 kW, a ukupna površina panela je<br />

A=3694.7 m 2 . Usvojeno je da su gubici usled <strong>za</strong>prljanja 4 %,<br />

gubici usled neuparenosti modula 3 % i da je efikasnost<br />

invertora 97.7 %.<br />

Da bi se uvažio uticaj temperature na efikasnost<br />

fotonaponske konverzije, izračunata je temperatura<br />

fotonaponskih modula na osnovu izmerene temperature<br />

ambijenta i poznate vrednosti <strong>za</strong> NOCT date od strane<br />

proizvođača (47.5 0 C). Temperature ambijenta i temperature<br />

modula u toku godine su prika<strong>za</strong>ne na slici 2.<br />

W [ kWh / dan ]<br />

3500<br />

3000<br />

2500<br />

2000<br />

1500<br />

1000<br />

500<br />

0<br />

0 50 100 150 200 250 300 350<br />

dan u godini<br />

Slika 3: Očekivane dnevne proizvedene energije fotonaponskog sistema<br />

35<br />

temperature [ 0 C ]<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

0<br />

T PV<br />

T amb<br />

CF [ % ]<br />

30<br />

25<br />

20<br />

15<br />

10<br />

5<br />

-10<br />

-20<br />

0 50 100 150 200 250 300 350<br />

dan u godini<br />

Slika 2: Temperatura fotonaponskih modula i izmerena temperatura ambijenta<br />

Procenjeno smanjenje efikasnosti fotonaponskog sistema<br />

usled povećanja temperature modula u odnosu na standardnu<br />

vrednost od 25 0 C na osnovu (2) i (3) je 6 %. Snaga koju PV<br />

sistem predaje mreži uz uvažavanje svih gubitaka odnosno<br />

efikasnosti konverzije je<br />

P = 495⋅0.96⋅0.97 ⋅0.94⋅ 0.977 kW=423 kW (18)<br />

AC<br />

Temperatura se menja u toku godine tako da je usvojeni<br />

procenat gubitaka usled nadtemperature relativan. Realan<br />

uticaj temperature je uključen prilikom proračuna <strong>proizvodnje</strong><br />

jer se poseduju merni podaci o vrednostima temperature <strong>za</strong><br />

celu godinu sa desetominutnom rezolucijom.<br />

Očekivane dnevne proizvedene energije solarnog sistema u<br />

toku godine su izračunate klasičnom metodom na osnovu<br />

snage koju PV sistem predaje mreži i broja sati “zenita sunca”.<br />

Ove energije su prika<strong>za</strong>ne na slici 3. Faktor kapaciteta<br />

fotonaponskog sistema je izračunat na osnovu (9) i prika<strong>za</strong>n je<br />

na slici 4.<br />

0<br />

0 50 100 150 200 250 300 350<br />

dan u godini<br />

Slika 4: Faktor kapaciteta fotonaponskog sistema<br />

Očekivane dnevne i ukupne proizvodene energije<br />

fotonaponskog sistema po sezonama i na godišnjem nivou su<br />

date u tabeli I.<br />

TABELA I. OČEKIVANE PROIZVEDENE ENERGIJE FOTONAPONSKOG SISTEMA I<br />

FAKTOR KAPACITETA<br />

Sezona<br />

Očekivana dnevna<br />

proizvedena energija<br />

[kWh]<br />

Očekivana energija<br />

proizvedena u sezoni<br />

[kWh]<br />

Faktor<br />

kapaciteta<br />

[%]<br />

prolećna 2036.6 187360 20.06<br />

letnja 2350.0 218550 23.14<br />

jesenja 1332.4 121250 13.12<br />

zimska 748.7 67384 7.37<br />

godišnje 1624.2 592830 15.99<br />

U tabeli I su takođe date vrednosti <strong>za</strong> faktor kapaciteta<br />

fotonaponskog sistema u posmatranim periodima vremena<br />

koje su izračunate na osnovu očekivanih proizvedenih energija<br />

i nazivne izlazne AC snage fotonaponskog sistema. Faktor<br />

- 190 -


kapaciteta fotonaponskog sistema na godišnjem nivou iznosi<br />

15.99%.<br />

Očekivane proizvedene energije fotonaponskog sistema su<br />

izračunate i na osnovu statističke metode koja se <strong>za</strong>sniva na<br />

karakterističnim sezonskim dijagramima snaga <strong>proizvodnje</strong>.<br />

Očekivane snage <strong>proizvodnje</strong> fotonaponskog sistema u toku<br />

godine su izračunate na osnovu insolacije na fotonaponski<br />

panel, površine panela i ukupne efikasnosti fotonaponskog<br />

sistema. U proračunu su korišćene srednje desetominutne<br />

vrednosti <strong>za</strong> insolaciju i temperaturu panela. Očekivane snage<br />

<strong>proizvodnje</strong> fotonaponskog sistema su prika<strong>za</strong>ne na slici 5.<br />

P [ kW ]<br />

450<br />

400<br />

350<br />

300<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0 50 100 150 200 250 300 350<br />

dan u godini<br />

Slika 5: Očekivana snaga <strong>proizvodnje</strong> fotonaponskog sistema<br />

Statističkom metodom nađeni su karakteristični dijagrami<br />

snage <strong>proizvodnje</strong> fotonaponskog sistema <strong>za</strong> četiri sezone:<br />

prolećnu (mart, april, maj), letnju (jun, jul, avgust), jesenju<br />

(septembar, oktobar, novembar) i zimsku (decembar, januar,<br />

februar). Karakteristični dani koji su dobijeni su: u prolećnoj<br />

sezoni 27. mart, u letnjoj sezoni 28. avgust, u jesenjoj sezoni<br />

28. novembar i u zimskoj sezoni 20. decembar. Na slikama 6<br />

do 9 su prika<strong>za</strong>ne očekivane snage <strong>proizvodnje</strong> fotonaponskog<br />

sistema u karakterističnim danima po sezonama.<br />

P [ kW ]<br />

450<br />

400<br />

350<br />

300<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0 5 10 15 20<br />

sat<br />

Slika 6: Očekivana snaga <strong>proizvodnje</strong> fotonaponskog sistema u<br />

karakterističnom danu u prolećnoj sezoni<br />

P [ kW ]<br />

P [ kW ]<br />

P [ kW ]<br />

450<br />

400<br />

350<br />

300<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0 5 10 15 20<br />

sat<br />

Slika 7: Očekivana snaga <strong>proizvodnje</strong> fotonaponskog sistema u<br />

karakterističnom danu u letnjoj sezoni<br />

450<br />

400<br />

350<br />

300<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0 5 10 15 20<br />

sat<br />

Slika 8: Očekivana snaga <strong>proizvodnje</strong> fotonaponskog sistema u<br />

karakterističnom danu u jesenjoj sezoni<br />

450<br />

400<br />

350<br />

300<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0 5 10 15 20<br />

sat<br />

Slika 9: Očekivana snaga <strong>proizvodnje</strong> fotonaponskog sistema u<br />

karakterističnom danu u zimskoj sezoni<br />

- 191 -


Očekivane dnevne i ukupne proizvodene energije<br />

fotonaponskog sistema po sezonama i na godišnjem nivou<br />

dobijene na osnovu karakterističnih dijagrama su date u<br />

tabeli II. U tabeli su takođe date vrednosti faktora kapaciteta<br />

fotonaponskog sistema u pojedinim periodima. Faktor<br />

kapaciteta na godišnjem nivou iznosi 16.75%.<br />

TABELA II.<br />

Sezona<br />

OČEKIVANE PROIZVEDENE ENERGIJE FOTONAPONSKOG<br />

SISTEMA I FAKTOR KAPACITETA<br />

Očekivana dnevna<br />

proizvedena energija<br />

[kWh]<br />

Očekivana energija<br />

proizvedena u sezoni<br />

[kWh]<br />

Faktor<br />

kapaciteta<br />

[%]<br />

prolećna 2155.6 198320 21.23<br />

letnja 2367.4 220170 23.31<br />

jesenja 1440.9 131120 14.19<br />

zimska 789.7 71072 7.78<br />

godišnje 1688.4 620682 16.75<br />

Poređenjem vrednosti proizvednih energija i faktora<br />

kapaciteta iz tabela I i II <strong>za</strong>ključuje se da predložena<br />

statistička <strong>metoda</strong> <strong>za</strong> <strong>procenu</strong> <strong>proizvodnje</strong> fotonaponskog<br />

sistema u razmatranom slučaju odstupa od klasične metode<br />

proračuna 5.85 % u prolećnoj sezoni, 0.74 % u letnjoj sezoni,<br />

8.14 % u jesenjoj sezoni, 5.14 % u zimskoj sezoni i 4.7 % na<br />

godišnjem nivou. Povećanje tačnosti predložene statističke<br />

metode se može povećati ako bi se koristili merni podaci iz<br />

više godina. Na taj način bi predložena <strong>metoda</strong> bila u<br />

prednosti u odnosu na klasičnu metodu proračuna zbog<br />

korišćenja statističke obrade podataka u predikciji <strong>proizvodnje</strong><br />

fotonaponskog sistema.<br />

V. ZAKLJUČAK<br />

S obzirom na stohastičku prirodu solarnog zračenja i na<br />

veliki broj uticajnih faktora (ambijentalni uslovi i performanse<br />

sistema), <strong>za</strong> <strong>procenu</strong> <strong>proizvodnje</strong> energije iz fotonaponskog<br />

sistema je pogodno koristiti statističke metode. U radu je<br />

prika<strong>za</strong>na <strong>metoda</strong> koja se <strong>za</strong>sniva na statističkoj obradi<br />

dnevnih dijagrama <strong>proizvodnje</strong> fotonaponskog sistema pri<br />

čemu se određuju karakteristični sezonski dijagrami snage<br />

<strong>proizvodnje</strong> koji se mogu smatrati reprezentativnim <strong>za</strong><br />

<strong>procenu</strong> proizvedene energije u posmatranim periodima<br />

vremena. Korišćenjem predložene metode izvršena je procena<br />

sezonske i godišnje <strong>proizvodnje</strong> fotonaponskog sistema<br />

instalisane snage 495 kW. Vrednosti proizvedenih energija i<br />

faktora kapaciteta na godišnjem nivou izračunate statističkom<br />

metodom odstupaju oko 5 % u odnosu na odgovarajuće<br />

vrednosti koje se dobijaju klasičnom metodom proračuna.<br />

LITERATURA<br />

[1] Gilbert M. Masters, Renewable and Efficient Electric Power Systems,<br />

Stanford University<br />

[2] D. L. King, W. E. Boyson, J. A. Kratochvil, “Analysis of Factors<br />

Influencing the Annual Energy Production of Photovoltaics Systems”,<br />

Photovoltaic Specialists Conference, 2002.<br />

[3] M. J. Ahmat, G. N. Tiawari, “Solar radiation models – review”,<br />

International Journal of Energy and Environment, Volume 1, Issue 3,<br />

2010, pp. 513-532.<br />

[4] S. Kaplanis, E. Kaplani, “A model to predict expected mean and<br />

stochastic hourly global solar radiation I(h, n j ) values”, Renewable<br />

Energy 32, 2007, pp. 1414-1425.<br />

[5] P. Bacher, “Short-term Solar Power Forecasting”, Kongens Lyngbay<br />

2008, IMM-2008-13.<br />

[6] O. Perpinan, E. Lorenzo, M. A. Castro, “On the calculation of energy<br />

produced by PV grid-connected system”, Progress in Photovoltaics<br />

Research and Applications (2007), Volume: 15, Issue: 3, Publisher:<br />

Wiley Online Library, Pages: 265-274.<br />

[7] M. Brabec, E. Pelikan, P. Krč, K. Eben, P. Musilek, “Statistical<br />

Modeling of Energy Production by Photovoltaic Farms”, IEEE Electrical<br />

Power and Energy Conference, 2010.<br />

[8] O. Perpinan, “Statistical analysis of performance and simulation of twoaxis<br />

tracking PV system”, Solar Energy, 83:11(2074-2085)<br />

[9] Ž. Đurišić, J. Mikulović, “Prethodna studija isplativosti gradnje<br />

fotonaponske elektrane na mikrolokaciji u Obrenovcu”, Januar 2012.<br />

[10] J. Nahman, D. Salamon, Z. Stojković, J. Mikulović, “Rationali<strong>za</strong>tion of<br />

Operation of An Industrial Network” Electric Power System Research,<br />

Volume 78, Issue 10, October 2008, Pages 1664-1671.<br />

Abstract— This paper presents a method for estimation of<br />

seasonal and annual electrical energy production from<br />

photovoltaic (PV) system. Based on the measured ten minute<br />

average horizontal solar radiation and ambiental temperature<br />

data during a year, dayly production diagrams for PV system are<br />

obtained. Dayly production diagrams are statistically processed<br />

in order to choose characteristic days, e.g. characteristic sesional<br />

power production diagrams which are representative for energy<br />

production estimate. By using the proposed method, seasonal and<br />

annual energy production from PV system of 495 kW rated<br />

power has been estimated.<br />

Key words- photovoltaic system; production diagram; energy<br />

production, capacity factor<br />

STATISTICAL METHOD FOR ESTIMATION OF<br />

ELECTRICAL ENERGY PRODUCTION FROM PV<br />

SYSTEM<br />

Jovan Mikulović<br />

Željko Đurišić<br />

- 192 -

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!