CN102777331A

CN102777331A - 风力发电机组风轮直径的确定方法

Info

Publication number: CN102777331A
Application number: CN2012102774506A
Authority: CN
Inventors: 刘卫
Original assignee: Guodian United Power Technology Co Ltd
Current assignee: Guodian United Power Technology Co Ltd
Priority date: 2012-08-06
Filing date: 2012-08-06
Publication date: 2012-11-14
Anticipated expiration: 2032-08-06
Also published as: CN102777331B

Abstract

本发明是有关于一种风力发电机组风轮直径的确定方法，包括以下步骤：以相同外部风资源条件下的已知风力发电机组作为参照机组；计算待设计机组的风轮直径，计算公式为D_B1＝D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2；通过迭代计算，对待设计机组风轮直径进一步修正，公式为D_B2-n=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A/(2M+D_B2-（n-1)))^1.5α。本发明可以1.5MW风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，从而为机组设计提供重要依据。

Description

风力发电机组风轮直径的确定方法

技术领域

本发明涉及风力发电技术领域，特别是涉及一种风力发电机组风轮直径的确定方法。

背景技术

对于风力发电机组设计者来说，确定风力发电机组的风轮直径是概念设计阶段非常核心的问题，一旦风轮直径设计选择失误，与市场需求偏差较大，则整个的设计工作都可能要推倒重来，因此如何正确选择和确定所设计的风力发电机组的风轮尺寸，或给出一个准确的尺寸范围，是一个非常重要和亟待解决的问题。

目前，国内主流机型为1.5MW风力发电机组，在国内外各类风资源环境中的应用非常普遍，对应各类风资源的风轮直径已经逐步趋于成熟化和系列化，但是2MW、2.5MW、3MW、5MW、6MW等机型开发时的风轮直径的选择一直是困扰设计师的一个问题，目前一些公司是根据国外厂家或其它公司的的叶片尺寸来直接应用，而由于不同叶片设计者设计叶片的参数、翼型等的不同，这种简单参照同类叶片的做法有较大的风险。

因此，鉴于1.5MW的风力发电机组可以基本涵盖我国各地的风资源条件的，如何以1.5MW风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，成为本领域针对风机直径尺寸选择方法的一个研究方向。

发明内容

本发明要解决的技术问题是提供一种风力发电机组风轮直径的确定方法，使其以1.5MW风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，从而克服现有的风力发电机组风轮直径的确定方法的不足。

为解决上述技术问题，本发明一种风力发电机组风轮直径的确定方法，包括以下步骤：

A.以相同外部风资源条件下的已知风力发电机组作为参照机组；

B.计算待设计机组的风轮直径，计算公式为

D_B1=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2，

其中，D_B1为计算得到的待设计机组的风轮直径，D_A为参照机组的风轮直径，P_EB为待设计机组的发电功率，P_EA为参照机组的发电功率，η_A为参照机组从气动输出功率到发电上网的总效率，η_B为待设计机组从气动输出功率到发电上网的总效率，C_PA为参照机组的叶片气动效率，C_PB为待设计机组的叶片气动效率；

C.通过迭代计算，对步骤B得到的待设计机组风轮直径进一步修正，公式为

D_B2-=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A)/(2M+D_B2-（n-1）))^1.5α，

其中，n为迭代次数，D_B2-n为进一步修正后的风轮直径，D_B2-（n-1)为上次迭代的风轮直径值，D_B2-0=D_B1，M为风轮最低点与地面的距离，α为风切变指数。

作为本发明的一种改进，所述的迭代次数随待设计机组的发电功率的增大而增大。

采用这样的设计后，本发明以1.5MW风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，从而可正确选择和确定出适应我国各地的风资源条件的风力发电机组的风轮尺寸，为机组设计提供重要依据。

具体实施方式

本发明风力发电机组风轮直径的确定方法，是以相同外部风资源条件下已知的模版风力发电机组作为参照机组（即机组A），来计算待设计机组（即机组B）的风轮直径D_B。

风力发电机组的气动功率是由以下公式给出：

P=1/2·π·ρ·(D/2)²·V³·C_p (1)

则有：

P_A＝1/2·π·ρ·(D_A/2)²·V_A ³·C_Pa (2)

P_B=1/2·π·ρ·(D_B/2)²·V_B ³·C_Pb (3)

P_A/P_B＝(1/2·π·ρ·(D_A/2)²·V_A ³·C_Pa)/(1/2·π·ρ·(D_B/2)²·V_B ³·C_Pb)

=(D_A ²·V_A ³·C_Pa)/(D_B ²·V_B ³·C_Pb) (4)

由于机组A和B的风机所处的外部风资源条件相同，则有V_B＝V_A，那么：

P_A/P_B=(D_A ²·C_PA)/(D_B ²·C_PB) （5）

D_A ²/D_B ²=(P_A·C_PB)/(P_B·C_PA) （6）

D_B=D_A·(P_B/P_A)^1/2·(C_PA/C_PB)^1/2 （7）

其中，D_B为待设计机组的风轮直径，D_A为参照机组的风轮直径，P_B为待设计机组的气动功率，P_A为参照机组的气动功率，C_PA为参照机组的叶片气动效率，C_PB为待设计机组的叶片气动效率。

如果对于同一设计者的产品，往往不同风轮直径的叶片气动效率也是基本一致的，即C_PA=C_PB，则等式（7）变为

D_B=D_A·（P_B/P_A）^1/2 （8）

但是，在实际设计中，不同功率等级的风机往往其轮毂中心位置和传动链的传动效率都不同，它们对风轮直径尺寸产生影响，因此，上述公式并不能直接得到可靠的结果，还需要对它们进行修正。

首先，从不同风机的传递效率进行风轮直径的修正

式（1）和式（2）的功率P_A和P_B，实际只是风力发电机组的气动功率，我们通常所说的风力发电机组的功率是指风力发电机组的上网功率，即发电功率，定为P_E，气动功率和电功率的关系为：

P_E=P·η （9）

则P=P_E/η （10）

η是风机从气动输出功率到发电上网的总得效率，在正常运行工况下，对于非直驱机组：η=ηm·ηe·ηc

对于直驱机组：η=ηe·ηc

ηm为齿轮箱效率，ηe为发电机效率，ηc为变频器效率，在机组概念设计阶段，这些部件的效率是可以得到的，因此，对于风力发电机组A和B来说，η_A和η_B是可以得到的，由式（10）可知：

P_A/P_B=(P_EA/η_A)/(P_EB/η_B)=(P_EA/P_EB)·(η_B/η_A) （11）

设经过传输效率修正后的风轮直径为D_B1，则由式（11）代入等式（7），得出：

D_B1=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2 （12）

其中，P_EB为待设计机组的发电功率，P_EA为参照机组的发电功率，η_A为参照机组从气动输出功率到发电上网的总效率，η_B为待设计机组从气动输出功率到发电上网的总效率。

之后，从风机轮毂高度导致的偏差对叶轮直径尺寸进行进一步修正。

在上述的推导中，我们假设机组A和B的风机轮毂中心高度是相同的，但是在实际中，不同风轮直径的轮毂高度基本是不一致的，必然受风切变的影响，因此需要对风轮直径D_B1进行修正。平均风速的风切变可以用赫尔曼的指数公式来表达：即：

V_(Z1)/V_(Z2)=(Z₁/Z ₂)^α （13）

式中，V_(Z1)和V_(Z2)分别为高度Z₁和Z₂处的风速，α为风切变指数，一般条件下，指数α取0.14。

由此可见，由于不同轮毂高度的平均风速不一致，对风机的功率有影响，一般讲，风机风轮直径越大，轮毂的中心位置越高。

假设机组B的风机从机组A的风机轮毂点Z₁(对应风速为V_A,风轮直径为D_B1),升高的设计高度Z₂,(对应风速为V_B,风轮直径为D_B2),要求保持功率不变，则由等式（2）和（3），得出：

(D_B2/D_B1)²=(V_A/V_B)³ （14）

一般讲在陆地，风机叶轮最低点距离地面为一定距离，设定为M米。Z_A＝M+D_A/2，B机组的轮毂高度为Z_B＝M+D_B1/2，由等式（13），可以得出

V_A/V_B＝(Z_A/Z_B)^α=((M+D_A/2)/(M+D_B1/2))^α

=((2M+D_A)/(2M+D_B1))^α （15）

将式（15）代入式（14）则有：

(D_B2/D_B1)²=((2M+D_A)/(2M+D_B1))^3α

修正后的叶轮直径为：

D_B2=D_B1·((2M+D_A)/(2M+D_B1))^1.5α （16）

将式（12）代入式（16），则有：

D_B2=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A)/(2M+D_B1))^1.5α（17）

最后，进行迭代计算，得出准确风轮直径尺寸

在这里，D_B1开始由（12）得出只是个初始值，然后需要对其进行迭代计算，即可得出较准确的叶轮直径。

设由（12）计算得出初始值为D_B1，令D_B2-0=D_B1，n为迭代次数并随待设计机组的发电功率的增大而增大，则,

1次迭代为：

D_B2-1=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A)/(2M+D_B2-0))^1.5α

2次迭代为：

D_B2-2=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A)/(2M+D_B2-1))^1.5α

n次迭代为：

D_B2-n=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A)/(2M+D_B2-（n-1）))^1.5α（18）。

以下列举具体应用举例来进一步说明本发明的实施过程及其计算结果。

第一步：按照公式（12），以1.5MW风力发电机组的风轮尺寸为基础，即1.5MW风轮直径为D_A，设C_PA＝C_PB，η_B=η_A，我们可以初步得出的2MW/2.5MW/2.5MW/3MW/3.6MW/5MW/6MW/12MW系列风力发电机组的风轮尺寸,见表1。

表1：

第二步：设α=0.14，M=28m，根据公式（18）进行迭代修正计算，

以风轮直径1为例，各次迭代修正见表2。

表2：

迭代次数	D_B2-0	D_B2-1	D_B2-2	D_B2-3	D_B2-4	D_B2-5
							2.0MW	80.8	79.5	79.7	79.6	79.6	79.6
5.0MW	127.8	125.0	125.4	125.3	125.3	125.3
							12.0MW	198.0	188.0	189.5	189.3	189.3	189.3

从表2中可以看出，在保留小数点一位的情况下，对上述发电功率的风力发电机组风轮直径经过5次迭代，即可达到0.1m的精度。因此按照5次迭代，对各风轮直径进行迭代修正计算，得出最终风轮直径数据见表3。

表3：

以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。

Claims

1.一种风力发电机组风轮直径的确定方法，其特征在于包括以下步骤：

B.计算待设计机组的风轮直径，计算公式为

D_B1=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2，

D_B2-n=D_A·(P_EB/P_EA)^1/2·（η_A/η_B）^1/2·(C_PA/C_PB)^1/2·((2M+D_A)/(2M+D_B2-（n-1）))^1.5α，

其中，n为迭代次数，D_B2为进一步修正后的风轮直径，D_B2-（n-1）为上次迭代的风轮直径值，D_B2-0=D_B1，M为风轮最低点与地面的距离，α为风切变指数。

2.根据权利要求1所述的风力发电机组风轮直径的确定方法，其特征在于所述的迭代次数随待设计机组的发电功率的增大而增大。