专利名称：风力发电机组风轮直径的确定方法
技术领域：
本发明涉及风力发电技术领域，特别是涉及一种风力发电机组风轮直径的确定方法。
背景技术：
对于风力发电机组设计者来说，确定风力发电机组的风轮直径是概念设计阶段非常核心的问题，一旦风轮直径设计选择失误，与市场需求偏差较大，则整个的设计工作都可能要推倒重来，因此如何正确选择和确定所设计的风力发电机组的风轮尺寸，或给出一个准确的尺寸范围，是一个非常重要和亟待解决的问题。 目前，国内主流机型为I.5丽风力发电机组，在国内外各类风资源环境中的应用非常普遍，对应各类风资源的风轮直径已经逐步趋于成熟化和系列化，但是2丽、2. 5丽、3丽、5丽、6丽等机型开发时的风轮直径的选择一直是困扰设计师的一个问题，目前一些公司是根据国外厂家或其它公司的的叶片尺寸来直接应用，而由于不同叶片设计者设计叶片的参数、翼型等的不同，这种简单参照同类叶片的做法有较大的风险。因此，鉴于I. 5丽的风力发电机组可以基本涵盖我国各地的风资源条件的，如何以I. 5丽风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，成为本领域针对风机直径尺寸选择方法的一个研究方向。
发明内容
本发明要解决的技术问题是提供一种风力发电机组风轮直径的确定方法，使其以
I.5MW风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，从而克服现有的风力发电机组风轮直径的确定方法的不足。为解决上述技术问题，本发明一种风力发电机组风轮直径的确定方法，包括以下步骤A.以相同外部风资源条件下的已知风力发电机组作为参照机组；B.计算待设计机组的风轮直径，计算公式为Dbi=Da · (PEB/PEA)1/2 · ( ηΑ/ ηΒ) 1/2 · (CPA/CPB)1/2，其中，Dm为计算得到的待设计机组的风轮直径，Da为参照机组的风轮直径，Peb为待设计机组的发电功率，Pea为参照机组的发电功率，ΠΑ为参照机组从气动输出功率到发电上网的总效率，ΠΒ为待设计机组从气动输出功率到发电上网的总效率，CpaS参照机组的叶片气动效率，Cpb为待设计机组的叶片气动效率；C.通过迭代计算，对步骤B得到的待设计机组风轮直径进一步修正，公式为Db2-=Da · (ΡΕΒ/ΡΕΑ)1/2 · ( ΠΑ/ ΠΒ) 1/2 · (CPA/CPB)1/2 · ((2M+DA)/(2M+DB2_ (n_0))L5a,其中，n为迭代次数,DB2_n为进一步修正后的风轮直径,DB2_ (η_υ为上次迭代的风轮直径值，Db2^Dbi, M为风轮最低点与地面的距离，a为风切变指数。作为本发明的一种改进，所述的迭代次数随待设计机组的发电功率的增大而增大。采用这样的设计后，本发明以I. 5丽风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，从而可正确选择和确定出适应我国各地的风资源条件的风力发电机组的风轮尺寸，为机组设计提供重要依据。
具体实施例方式本发明风力发电机组风轮直径的确定方法，是以相同外部风资源条件下已知的模版风力发电机组作为参照机组(即机组A)，来计算待设计机组(即机组B)的风轮直径Db。风力发电机组的气动功率是由以下公式给出P=l/2 · Ji · P · (D/2)2 · V3 · Cp(I)则有 pA = 1/2 · Ji · P · (Da/2)2 · Va3 · Cpa (2)Pb= I/2 · Ji · P · (Db/2)2 · Vb3 · Cpb (3)PA/PB = (1/2 · Ji · p · (Da/2)2 · Va3 · CPa)/(l/2 · Ji · p · (DB/2)2 · Vb3 · Cpb)= (D/ · Va3 · CPa) / (Db2 · Vb3 · Cpb)(4)由于机组A和B的风机所处的外部风资源条件相同，则有Vb = Va，那么PA/PB= (Da2 · CPA) / (Db2 · Cpb)(5)Da2/Db2= (Pa · CPB) / (Pb · Cpa)(6)Db=Da · (PB/PA)1/2 · (CPA/CPB)1/2(7)其中，Db为待设计机组的风轮直径，Da为参照机组的风轮直径，Pb为待设计机组的气动功率，Pa为参照机组的气动功率，Cpa为参照机组的叶片气动效率，Cpb为待设计机组的叶片气动效率。如果对于同一设计者的产品，往往不同风轮直径的叶片气动效率也是基本一致的，即Cpa=Cpb，则等式(7)变为Db=Da · (PB/PA) 1/2 (8)但是，在实际设计中，不同功率等级的风机往往其轮毂中心位置和传动链的传动效率都不同，它们对风轮直径尺寸产生影响，因此，上述公式并不能直接得到可靠的结果，还需要对它们进行修正。首先，从不同风机的传递效率进行风轮直径的修正式(I)和式(2)的功率Pa和Pb，实际只是风力发电机组的气动功率，我们通常所说的风力发电机组的功率是指风力发电机组的上网功率，即发电功率，定为Pe，气动功率和电功率的关系为Pe=P · η(9)贝丨JP=Pe/η(10)n是风机从气动输出功率到发电上网的总得效率，在正常运行工况下，对于非直 n = nm · ne · nc对于直驱机组η = η e · η cη m为齿轮箱效率，η e为发电机效率，η c为变频器效率，在机组概念设计阶段，这些部件的效率是可以得到的，因此，对于风力发电机组A和B来说，^和ηΒ是可以得到的，由式(10)可知pa/pb=(pea/nA)/(PEB/nB) = (PEA/PEB) · (nB/nA)(H)设经过传输效率修正后的风轮直径为Dbi，则由式(11)代入等式(7 )，得出Dbi=Da · (ΡΕΒ/ΡΕΑ)1/2 · ( ΠΑ/ ΠΒ) 1/2 · (CPA/CPB)1/2 (12)其中，Peb为待设计机组的发电功率，Pea为参照机组的发电功率，nA为参照机组从气动输出功率到发电上网的总效率，jIb为待设计机组从气动输出功率到发电上网的总效率。之后，从风机轮毂高度导致的偏差对叶轮直径尺寸进行进一步修正。
在上述的推导中，我们假设机组A和B的风机轮毂中心高度是相同的，但是在实际中，不同风轮直径的轮毂高度基本是不一致的，必然受风切变的影响，因此需要对风轮直径Dbi进行修正。平均风速的风切变可以用赫尔曼的指数公式来表达即V⑶/V(Z2)= (Z1A 2) α(13)式中，V(Z1)和V(Z2)分别为高度21和&处的风速，α为风切变指数，一般条件下，指数α取O. 14。由此可见，由于不同轮毂高度的平均风速不一致，对风机的功率有影响，一般讲，风机风轮直径越大，轮毂的中心位置越高。假设机组B的风机从机组A的风机轮毂点Z1 (对应风速为\，风轮直径为Dbi)，升高的设计高度Z2,(对应风速为VB，风轮直径为DB2)，要求保持功率不变，则由等式(2)和(3)，得出(DB2/DB1)2=(VA/VB)3(14)一般讲在陆地，风机叶轮最低点距离地面为一定距离，设定为M米。Za = M+Da/2,B机组的轮毂高度为Zb = M+DB1/2,由等式(13)，可以得出 /\ = (ΖΑ/ΖΒ) α = ( (M+Da/2) / (M+DB1/2)) α= ((2M+Da)/(2M+DB1)) α(15)将式(15)代入式(14)则有(DB2/DB1)2二 ( (2M+Da) / (2M+DB1))3 a修正后的叶轮直径为Db2=Dbi · ((2M+Da)/(2M+DB1))L5a (16)将式(12)代入式(16)，则有Db2=Da · (PEB/PEA)1/2 · ( ΠΑ/ ΠΒ) 1/2 · (CPA/CPB)1/2 · ((2M+DA)/(2M+DB1))L5a (17)最后，进行迭代计算，得出准确风轮直径尺寸在这里，Dbi开始由(12)得出只是个初始值，然后需要对其进行迭代计算，即可得出较准确的叶轮直径。设由(12)计算得出初始值为Dbi，令Dici=Dbi,η为迭代次数并随待设计机组的发电功率的增大而增大，则，I次迭代为Dbh=Da · (ΡΕΒ/ΡΕΑ)1/2 · ( ΠΑ/ ΠΒ) 1/2 · (CPA/CPB)1/2 · ((2M+DA)/(2M+DB2_0))L5a2次迭代为DB2_2=Da · (PEB/PEA)1/2 · ( ΠA/ ΠB) 1/2 · (CPA/CPB)1/2 · ( (2M+Da) / (2M+DB2J )L5a
η次迭代为DB2_n=DA · (PEB/PEA)1/2 ·( nA/ nB)1/2 · (CPA/CPB)1/2 · ((2M+Da)/(2M+DB2_ ( _!)))L5a (18)。以下列举具体应用举例来进一步说明本发明的实施过程及其计算结果。第一步按照公式(12)，以I. 5MW风力发电机组的风轮尺寸为基础，S卩I. 5MW风轮直径为 Da，设 Cpa = CPB, nB=nA,我们可以初步得出的 2MW/2. 5MW/2. 5MW/3MW/3. 6MW/5MW/6MW/12MW系列风力发电机组的风轮尺寸，见表I。表I :
机组功率风轮直径I 风轮直径2 风轮直径3 风轮直径4 风轮直径5 (MW)(m)(ii)(ii)(ii){m)
1.5 70.0 77.0 82.0 86. 0 93.02.0 80.8 88.9 94.7 99.3 107.4
5.0127.8 140.6 149.7 136.0 169.8
12.0198.0 217.8 231.9 243.2 263.0第二步设a =0. 14, M=28m,根据公式(18)进行迭代修正计算，以风轮直径I为例，各次迭代修正见表2。表2
迭代次数"I Db2_0I Db2-1I Db2_2I Db2_3I 0β2-4I DB2_5 I
2.OMW80879Γ579 7 θΓθ θΓθ θΓθ
5. OMW127.8125.0125.4125. 3125. 3125. 3~
~12. OMW198.0188.0189. 5189. 3189. 3189. 3~从表2中可以看出，在保留小数点一位的情况下，对上述发电功率的风力发电机组风轮直径经过5次迭代，即可达到O. Im的精度。因此按照5次迭代，对各风轮直径进行迭代修正计算，得出最终风轮直径数据见表3。表3
权利要求
1.一种风力发电机组风轮直径的确定方法，其特征在于包括以下步骤 A.以相同外部风资源条件下的已知风力发电机组作为参照机组； B.计算待设计机组的风轮直径，计算公式为 Dbi=Da (PEB/PEA)1/2 ( nA/ nB)1/2 (CPA/CPB)1/2， 其中，Dbi为计算得到的待设计机组的风轮直径，Da为参照机组的风轮直径，Peb为待设计机组的发电功率，Pea为参照机组的发电功率，nA为参照机组从气动输出功率到发电上网的总效率，nB为待设计机组从气动输出功率到发电上网的总效率，Cpa为参照机组的叶片气动效率，Cpb为待设计机组的叶片气动效率； C.通过迭代计算，对步骤B得到的待设计机组风轮直径进一步修正，公式为 Db2-H=Da (PEB/PEA)1/2.(nA/nB)1/2. (CPA/CPB)1/2. ((2M+da)/(2M+dB2_^0))1-5°^ 其中，n为迭代次数，Db2为进一步修正后的风轮直径，DB2_ (n_D为上次迭代的风轮直径值，Dkmi=Dbi, M为风轮最低点与地面的距离，Ct为风切变指数。
2.根据权利要求I所述的风力发电机组风轮直径的确定方法，其特征在于所述的迭代次数随待设计机组的发电功率的增大而增大。
全文摘要
本发明是有关于一种风力发电机组风轮直径的确定方法，包括以下步骤以相同外部风资源条件下的已知风力发电机组作为参照机组；计算待设计机组的风轮直径，计算公式为DB1＝DA·(PEB/PEA)1/2·(ηA/ηB)1/2·(CPA/CPB)1/2；通过迭代计算，对待设计机组风轮直径进一步修正，公式为DB2-n=DA·(PEB/PEA)1/2·(ηA/ηB)1/2·(CPA/CPB)1/2·((2M+DA/(2M+DB2-(n-1)))1.5α。本发明可以1.5MW风力发电机组的风轮直径尺寸为基础，来确定其它风力发电机组的风轮直径尺寸，从而为机组设计提供重要依据。
文档编号F03D11/00GK102777331SQ201210277450
公开日2012年11月14日 申请日期2012年8月6日 优先权日2012年8月6日
发明者刘卫 申请人:国电联合动力技术有限公司