风切变与湍流度对风含功率的影响-20141210.docx

风切变与湍流度对风含功率的影响罗绍卓1，陈卓1，原野1（1.北京天源科创风电技术有限责任公司，100176）摘 要：以风力机前方自由来流风为研究对象，通过在叶轮扫风面内风速矢量在时间上积分，结合指数切变方程及能量定理，考察纵向风切变及湍流度对风含功率的影响，建立了计算模型及影响系数的表达式。研究表明，在同时考虑风切变及湍流时，在通常风切变指数范围内（=0.10.2），切变独立影响系数kS约为0.99。湍流度与风切变共同对风含功率产生影响，湍流度的影响随切变指数增大而增大，轮毂高度处风速超过15m/s时，影响系数低于0.01。关键词：风力发电；风含功率；影响系数；风力机中图分类号：TK513.5文献标识代码：A0 引 言随着风电市场的发展及风电服务项目的深化，对已投入运行的风力机进行功率提升成为风电投资商越来越重视的项目。研究风力机前方来流风的风含功率是有助于制定功率提升方案。风场地形的复杂化，使得风剪切、塔影效应和湍流对风力机的影响越来越显著12。风资源最显著的特征是其不均匀性及随机性，主要表现为风资源在空间和时间上的分布不均3。风切变是指在大气边界层中，由于受地球表面的作用，平均风速随高度变化的现象4，是风资源在空间上分布不均的典型特征之一。湍流指的是由于地形差异及气温变化的热效应，引起的风速在一定时间内的风速波动，表征了风资源在时间上的分布不均特性。图1 积分区域Fig.1 Integral domain风含功率指的是不受阻挡流过叶轮扫风面积上的自由来流风所蕴含的功率。在不考虑空气密度的变化时，由式（1）计算：Pwind=2T-T2T2V2udsdt (1)式中，Pwind风含功率；空气密度；T积分时长； V合成风速；u风速纵向分量。如图1所示建立积分坐标系，其中z轴垂直于地面，在半径为R的圆域内对式（1）进行积分，得式（2）。Pwind=2T-T2T20R02V2urdrddt (2)风电行业常用下式简单计算风含功率。Pwind=uhub3R22 (3)式中，uhub轮毂高度处平均风速。1 风切变对稳态风的风含功率影响对于均匀定常流动，仅有纵向分量u。考虑指数形式的纵向风切变5：uz=uhubzzhub (4)式中， zhub轮毂高度；风切变指数。采用极坐标系，令z=zhub+rsin，把式（4）代入式（2），得：Pwind=uhub320R02zhub+rsinzhub3rdrd(5) 定义风切变影响系数为：ks=0R02zhub+rsinzhub3rdrdR2 (6)取uhub=322m/s，zhub=75m, R=43.5m，=0.10.2，算得ks如图2所示。从图2中可见，在上述典型的风力机设计参数范围内，ks约等于0.99，并在=0.150.175间达到极小值。图2 不同风切变指数及平均风速下的ksFig.2 ks for different wind shear exponents and averaged wind speeds2 湍流度对风含功率影响空间上某点的风速矢量可以分解为三分量，分别为纵向u，横向v，及竖向。时间轴上，每一个分量由平均风速及脉动风速组成，通常取定纵向为主风向（垂直于叶轮平面），三个风速分量由下式表示：u=u+uv=v = (7)式中，u纵向分量平均风速；u纵向分量脉动风速；v横向分量脉动风速；垂直分量脉动风速。三个脉动风速分量可以用Kaimal谱（单边谱）描述5：Skf=k24Lkuhub1+6fLkuhub53 (8)式中，k速度分量索引号，1为纵向，2为横向，3为竖向；f频率；k三个脉动速度的标准差；Lk积分尺度因子。式（8）满足以下关系：k2=0Skfdf (9)此外还有下式成立6：1=Iu (10)当轮毂高度超过60米（IEC 64100-1，2005），上述标准差及积分因子可以根据表1进行计算。表1 不同速度分量的标准差及积分尺度因子Table1 Standard deviation and integral scale factor for three components of velocity vector索引号k标准差k积分尺度因子Lk11340.220.81113.430.5127.72根据表1可将式（8）改写为：Skf=k2124Lkuhub1+6fLkuhub53 (11)式中，k系数，k=1时为1，2时为0.8，3时为0.5。假定，功率谱在整个叶轮旋转平面内一致，在考虑的时间范围内，平面内各点的平均风速不变，把式（7）代入式（2），假设任意两个脉动风速分量互不相关，并忽略脉动风速的三次乘积项，得Pwind=20R02u3+ulimT+1T-T2T23u2+v2+2dTdrd (12)根据Parseval等式7，利用式（9）及式（11），可将式（12）改写为：Pwind=20R02u3drd+20R02u12312+22+32drd (13)利用式（10），定义湍流度风切变影响系数：kIS=Ihub2312+22+32R20R02zhub+rsinzhubrdrd(14)结合式（6）、式（14），式（13）可改写为：Pwind=u3R22 ks+kIS (15) 当0.1,0.2时，k=ks+kI0.99+kI (16)文献6给出不同轮毂高度处平均风速与湍流度的关系：Ihub=I1515+auhub1+auhub (17)其中，I15轮毂高度平均风速为15m/s时的参考湍流强度，取0.18或0.16；a 系数，当I15=0.18时，取2，当I15=0.16时，取3。将式（17）代入式（14）得kIS=312+22+32R215+auhub1+auhub2 I1520R02zhub+rsinzhubrdrd(18)取uhub=322，zhub=75, R=43.5，=0.1，0.15，0.2，I15=0.18，0.16时，kIS变化如图3所示。图中可见，kIS随着剪切系数增大而增大。相同风切变水平下，kIS随风速增大而减小，且很快收敛。轮毂高度平均风速超过15m/s时，kIS集中在0.01以下；平均风速低于5m/s且湍流度较高时，kIS将超过0.04。图3 不同湍流强度及平均风速下的kIsFig.3 kIs for different turbulent intensities and averaged wind speeds3 结论为了提高风力发电机功率提升量化分析准确性，本文对风速矢量在时间及叶轮平面内进行积分考察风切变及湍流度对风含功率的影响。本文利用风力机设计标准、规范中的风切变模型、Kaimal风功率谱模型，通过Parseval等式建立了风切变及湍流度对风含功率的影响函数计算模型，结合规范中湍流度设计要求对典型风切变及切入切出风速范围内对影响系数进行了计算，得出以下结论：1）影响系数可以分解为风切变影响系数ks以及湍流度风切变影响系数kIs两项。2）ks在设计参数变化范围内越等于0.99，是构成影响系数的最主要成分，可看作是直流分量。3）kIs决定着实际风含功率偏离常用简化计算式，即式（3），的程度。4）采用风力机规范中的湍流度设计方法计算结果表明，kIS随着剪切系数增大而增大。相同风切变水平下，kIS随风速增大而减小，且很快收敛。轮毂高度平均风速超过15m/s时，kIS集中在0.01以下；平均风速低于5m/s且湍流度较高时，kIS将超过0.04参考文献1 Dolan D S L,Lehn P W.Simulaion model of wind turbine 3p torque oscillations due to wind shear and tower shadowJ. IEEE Transactions on Energy Conversion，2006，21(3)：717-724.2 Eggers A J，Digumarthi R，Chaney K. Wind shear and turbulence effects on rotor fatigue and loads controlJ.ASME Journal of Solar Energy Engineering，2003，125(4)：402-4.3 Tony Burton，Nick Jenkins，David Sharpe，Ervin Bossanyi. 风能技术M.武鑫译. 北京：科学出版社，2014.14 贺德馨，陈坤，张亮亮，等.风工程与工业空气动力学M.北京：国防工业出版社，2006.5 IEC61400-1，Wind turbines-Part1：Design requirmentsS. Switzerland，International Electrotechnical Commission，2005.6 GL2010，Guidline for Certification of Wind TurbinesS. Ge-rmany，Germanischer Lloyd Industrial Services Gmbh， 2010.7 张元林. 积分变换M. 北京：高等教育出版社，2011.EFFECT OF WIND SHEAR AND TURBULENCE INTENSITY ON WIND POWERLuo Shaozhuo1, Chen zhuo1, Yuan ye1(1.CHINA TIANYUAN NEW ENERGY TECHNOLOGY CO.,LTD, Beijing 100176, China)Abstract: To analyze the effect of changes of vertical wind shear and the turbulence intensity on wind power, the free stream in front of the wind turbine is chosen as the research object, the computational model and the effect coefficient are derived via integral of the wind speed vectors on the surface of the rotor over time combined with exponential power law and the energy theorem. The research shows that while the wind shear index () is in the range of 0.10.2, the independent wind shear coefficient kS would be around 0.99. The turbulence intensity and the wind shear would affect the wind power at the same time; the effect of change in turbulence intensity would get more signif