三种风速威布尔分布参数算法的比较概要.docVIP

三种风速威布尔分布参数算法的比较概要

三种风速威布尔分布参数算法的比较

徐卫民,孔新红,桂保玉

(江西省气象科学研究所,江西南昌330046

摘要:介绍计算威尔分布参数的累积分布函数拟合法、平均风速和标准差估计法和平均风速和最大风速估计等三种算法,并应用此算法计算了都阳气象站的风速威布尔分布参数。根据分布参数拟合了都阳县气象站的三种风速概率分布,将拟合的风速概率分布与同期的风速实际频率分布结果进行相关分析,依据相关系数判断拟合效果的好坏。通过比较得到了以下结论:平均风速和标准差估计法效果最好,累积分布函数拟合法次之,由于最大风速变化比较随机,平均风速和最大风速估计法效果波动最大,整体效果差。通过多年最大风速的平均数与平均风速计算,能减少最大风速抽样的随机性误差,结果更具代表性。关键词:风速;分布规律;威布尔;比较

0引言

近年来,我国并网运行的大中型风力发电厂建设逐渐纳入有计划、规范化发展的轨道。鄱阳湖风力发电站建设项目已经纳人江西省“十一五”规划重大建设项目中。为此,有必要开展风能分析及风电场设计等方面的研究工作。威布尔(Weibull分布双参数曲线,是一种形式简单且又能较好拟合实际风速分布的概率模型,只要给定了威布尔分布参数k和c,风速的分布形式便给定了,而毋需逐一查阅和统计所有的风速观测资料,可方便地求得平均风能密度、有效风能密度、风能可利用小时数,给实际使用带来许多方便

?

???

???=≤k

g

gc

Pυυυexp1((2

取对数整理后,有

([]

{}ckkpgglnln1lnln?=≤??υυυ(3

令([]{}

1lnlngPyυυ≤??=,gxυln=,ckaln?=,,kb=于是参数k和c可以由最小二乘法拟合bxay+=得到。具体做法如下:

将观测到的风速出现范围划分成n个风速间隔1~0υ,21~υυ,…,nnυυ~1?。统计每个间隔中风速观测值出现的频率1f,2f,…,nf和累积频率11fp=,212fpp+=,…,

nnnfpp+=?1取变换,并令

ln1υ=ix

([]11lnlnpyi??=

lncka?=(4kb=(5

因此,根据式(3及述风速累积频率观测资料,便可得到a、b的最小二乘估计值

∑∑??=

222i

i

i

ii

i

i

xxnyxxyxa

∑∑∑∑∑?+?=

22i

i

i

ii

i

xxny

xnyxb由式(4、式(5,得

?

??????=bacexpbk=

1.2根据平均风速?

υ和标准差vS估计威布尔分布参数

原理为:

(([]{}

11/122

?+Γ+Γ=????????μσ(6

从式(6可见μσ/是k的函数。因此当知道了分布的均值和方差,便可求解k。由于直接用

μσ/求解k比较困难,用上式的近似关系式求解k,即

086

.1??

??

?????=μσk(7

而由式(7有

1

1/(/k

c+Γ=μ(8

以平均风速_

υ估算μ,样本标准差VS估计σ即

∑=

iN

1

_

υυ∑?

?

?????=

2

_1SυυiN

式中:iυ为计算时段中每次的风速观测值;N为观测总次数。

由式(7和式(8便可求得k和c的估计值。1.3用平均风速和最大风速估计威布尔分布参数

设maxυ为时间T内观测到的10min平均最大风速,它出现的概率

(TcPk1

expmaxmax=

???????????????=≥υυυ(9

对式中(9作逆变换,得

((

T1/lnmax

+Γ=μ

υ(10

因此知道了maxυ和_υ,以_

υ作为μ的估计值,由式(10就能解出k。c则由式(8求得。

2计算结果[8]

收集了鄱阳县气象站1991-2000年的自记风资料,分别用三种方法计算每年及整10年的风速威布尔分布参数k、c,结果见表1。

表1计算结果

3效果检验[9]

3.1检验方法

利用鄱阳气象站1991-2000年逐小时风资料统计出每年及整10年中0-18m/s间,以1m/s为间隔的风速分布频率。

根据用3种方法估算出的威布尔分布参数k和c按公式计算出3种1991-2000年间每年及整10年中0-18.0m/s间,以1m/s为间隔的风速概率分布。

用拟合出的3种概率分布与同期实测风的分布频率进行相关分析,相关系数高