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