Rayleigh raspodjela brzine vjetra

Rayleigh raspodjela brzine vjetra je statistička metoda s kojom se služimo za matematičko opisivanje mjerenih podataka brzine vjetra. To je ustvari poseban slučaj Weibullove raspodjele.

Weibullova raspodjela

Weibullova raspodjela je dvoparametarska i varijacije u brzini vjetra karakterizirane su dvjema funkcijama: funkcijom gustoće razdiobe vjerojatnosti i kumulativnom funkcijom raspodjele.^{^[1]}

Funkcija gustoće razdiobe vjerojatnosti

Ona označava vjerojatnost da je brzina vjetra upravo jednaka vrijednosti brzine v i prikazana je sljedećom jednadžbom:

[math]\displaystyle{ f(v)=\frac{k}{c}\left(\frac{v}{c}\right)^{k-1}e^{-\left(\frac{v}{c}\right)^k} }[/math]

[math]\displaystyle{ c }[/math] - faktor skale, [m/s]
[math]\displaystyle{ k }[/math] - faktor oblika (bezdimenzijski faktor)
[math]\displaystyle{ v }[/math] - brzina vjetra, [m/s]

Kumulativna funkcija raspodjele

Ona daje vjerojatnost da je brzina vjetra manja ili jednaka brzini v i prikazana je jednadžbom:

[math]\displaystyle{ F(v)=\int_{0}^{v}f(v)\,dv }[/math]

Uvrštavanjem [math]\displaystyle{ f(v) }[/math] u gornju formulu dobijemo

[math]\displaystyle{ F(v)=1-e^{-\left(\frac{v}{c}\right)^k} }[/math]

Faktor oblika k i njegov utjecaj

Male vrijednosti faktora oblika k označavaju brzine vjetrova koje jako variraju i zato nam funkcija gustoće razdiobe vjerojatnosti za taj slučaj daje manju vjerojatnost pojavljivanja određene brzine vjetra v, dok su konstantne vrijednosti brzine vjetra karakterizirane većim vrijednostima faktora. To znači da će vjerojatnost pojavljivanja brzine vjetra biti negdje oko ukupne srednje brzine vjetra, odnosno sa manjim vrijednostima faktora oblika k krivulja nam postaje položitija, a sa većim vrijednostima strmija.^{^[2]}

Rayleighova raspodjela je ustvari specijalni slučaj Weibullove raspodjele u kojem faktor oblika k poprima vrijednost 2, te je ona jednoparametarska. Tada nam funkcija gustoće razdiobe i kumulativna funkcija glase:

[math]\displaystyle{ f(v)=\frac{2v}{c^2}e^{-\left(\frac{v}{c}\right)^2} }[/math]

[math]\displaystyle{ F(v)=1-e^{-\left(\frac{v}{c}\right)^2} }[/math]

Razlog primjene

Poznato je da brzina vjetra konstantno varira. Da bi bili u mogućnosti predvidjeti proizvodnju električne energije iz vjetroturbine, moramo znati učestalost i intenzitet vjetra. To se može postići na način da brzinu vjetra mjerimo anemometrom i računamo srednju vrijednost brzine vjetra u vremenskim intervalima od deset minuta. Ti podaci mogu biti sortirani u raspone brzine vjetra od 1 m/s^{^[3]}. Ako bi takvo mjerenje provodili duži period (recimo na godišnjoj razini) dobili bi histogram brzine vjetra i vidjeli bi da je Rayleighova raspodjela brzine vjetra često dobra aproksimacija tih izmjerenih vrijednosti^{^[4]}.

funkcija gustoće Rayleighove raspodjele

Faktor skale c određuje se iz izraza:^{^[5]}

[math]\displaystyle{ c=\frac{\bar{v}}{\Gamma\left(1+\frac{1}{k}\right)} }[/math]

[math]\displaystyle{ \bar{v}=\frac{1}{N}\sum_{i=1}^\N v_i }[/math] - srednja brzina vjetra, [m/s]
[math]\displaystyle{ \Gamma\left(1+\frac{1}{k}\right) }[/math] - gama funkcija

Za Rayleighovu raspodjelu faktor skale c iznosi:

[math]\displaystyle{ c=\frac{2}{\sqrt{\Pi}}\bar{v} }[/math]

Dodatne mogućnosti

Pomoću funkcije gustoće raspodjele vjerojatnosti moguće je odrediti vjerojatnost da su brzine vjetra iznad neke brzine v₁. To je vrlo korisno jer na taj način možemo znati isplati li nam se investicija u vjetroelektranu. U Hrvatskoj, investicija u vjetroelektranu se isplati ako je prosječna godišnja brzina vjetra veća od 5,8 [m/s]^{^[6]}.

[math]\displaystyle{ p(v \ge v_1)=\int_{v_1}^\infty f(v)\,dv = e^{-\left(\frac{v}{c}\right)^2} }[/math]

[math]\displaystyle{ p(v \ge v_1) }[/math] - vjerojatnost da je brzina vjetra v veća od neke minimalne brzine v₁, ispod koje nam se ne isplati izgradnja vjetroelektrane

Moguće je još odrediti vjerojatnost da se brzina vjetra nalazi između neke dvije brzine v₁ i v₂. To je pogodno za slučaj kada želimo isplativost investicije u vjetroelektranu, ali ne želimo da je izložena brzinama većim od brzine v₂ koja bi mogla dovesti do oštećenja lopatica turbine:

[math]\displaystyle{ p(v_1\lt v\lt v_2)=F(v_2)-F(v_1)=\int_{v_1}^{v_2} f(v)\,dv = e^{-\left(\frac{v_1}{c}\right)^2} - e^{-\left(\frac{v_2}{c}\right)^2} }[/math]

Izvori

[1] ttp://www.intechopen.com/books/wind-energy-management/weibull-distribution-for-estimating-the-parameters

[2] ttp://data.naturalcapitalproject.org/invest-releases/documentation/2_5_2/wind_energy.html

[3] ttp://wind-data.ch/tools/weibull.php?lng=en]

[4] ttp://www.wind-power-program.com/wind_statistics.htm]

[5] ttp://www.intechopen.com/books/wind-energy-management/weibull-distribution-for-estimating-the-parameters

[6] ttp://project-rec.eu/obnovljivi-izvori-energije/energija-vjetra/

[1]

[2]

[3]

[4]

[5]

[6]