风力发电系统建模与仿真

《风力发电系统建模与仿真》由会员分享，可在线阅读，更多相关《风力发电系统建模与仿真（25页珍藏版）》请在金锄头文库上搜索。

1、风力发电系统建模与仿真风力发电系统建模与仿真摘要：风力发电作为一种清洁的可再生能源利用方式，近年来在世界范围内获得了飞速的发展。本文基于风力机发电建立模型，主要完成了以下工作：（1）基于风资源特点，建立了以风频、风速模型为基础的风力发电理论基Tali础；（2）运用叶素理论，建立了变桨距风力机机理模型；（3）分析了变速恒频风力发电机的运行区域与变桨距控制的原理与方法，并给出了机组的仿真模型，为风力发电软件仿真奠定了基础；（4）搭建了一套基于PSCAD/EMTDC仿真软件的风力发电系统控制模型以及完整的风力发电样例系统模型，并且已初步实现风力机特性模拟功能。关键词：风力发电；风频；风速；风力机；变

2、桨距；建模与仿真1风资源及风力发电的基本原理1.1风资源概述（1）风能的基本情况1风的形成乃是空气流动的结果。风向和风速是两个描述风的重要参数。风 向是指风吹来的方向，如果风是从东方吹来就称为东风。风速是表示风移动的 速度即单位时间内空气流动所经过的距离。风速是指某一高度连续10min所测得各瞬时风速的平均值。一般以草地上 空10m高处的10min内风速的平均值为参考。风玫瑰图是一个给定地点一段时 间内的风向分布图。通过它可以得知当地的主导风向。风能的特点主要有：能量密度低、不稳定性、分布不均匀、可再生、须在 有风地带、无污染、分布广泛、可分散利用、另外不须能源运输、可和其它能 源相互转换等。

3、（2）风能资源的估算风能的大小实际就是气流流过的动能，因此可以推导出气流在单位时间内 垂直流过单位截面积的风能，即风能密度，表示如下：= 0.5pv3（1-1）i=J式中e-风能密度（W Im2 ）是描述一个地方风能潜力 的最方便|有价值的，量；p空气密度(kg / m 3)；v风速(m / s )。由于风速是一个随机性很大的量，必须通过一段时间的观测来了解它的平 均状况，一个地方风能潜力的多少要视该地常年平均风能密度的大小。因此需 要求出在一段时间内的平均风能密度，这个值可以将风能密度公式对时间积分 后平均来求得。有效风能密度还可根据下式求得(1-2)(1-3)=j v2 0.5pv 3 P

4、(v)dvvi式中，V 1启动风速(m/s );V2停机风速(m/s )；P(v)有效风速范围内的条件概率分布密度函数。平均风能密度则可用下式求得：1=j 0.5 pv3P(v)dt1.2风力发电的基本原理风能具有一定的动能，通过风轮机将风能转化为机械能，拖动发电机发电。 风力发电的原理是利用风带动风车叶片旋转，再通过增速器将旋转的速度提高 来促使发电机发电的。依据目前的风车技术，大约3m/s的微风速度便可以开始 发电。风力发电的原理说起来非常简单，最简单的风力发电机可由叶片和发电 机两部分构成如图1-1所示。空气流动的动能作用在叶轮上，将动能转换成机 械能，从而推动片叶旋转，如果将叶轮的转轴

5、与发电机的转轴相连就会带动发 电机发出电来。图IT风力发电原理图1.3风力发电的特点风力发电具有以下特点： 可再生的洁净能源； 建设周期短，装机规模灵活，可根据资金情况决定一次装机规模，有一 台资金就可以安装一台投产一台； 可靠性高，把现代高科技应用于风力发电机组使其发电可靠性大大提高，中、大型风力发电机组可靠性从80年代的50%提高到了 98%，高于火力发电且机组寿命可达20年； 造价低，运行维护简单，实际占地面积小； 发电方式多样化，既可并网运行，也可以和其他能源如柴油发电、太阳 能发电、水利发电机组形成互补系统，还可以独立运行； 单机容量小2风能及风力机系统模型的建立2.2风频模型风速具

6、有明显的随机性和间歇性。为了较精确地描述风速及其变化特性， 引入风频分布的概念。风频分布就是风速的统计概率分布，是衡量风能资源分 布特性的重要指标，它反映了风电场某个时段每一风速出现的概率，可以通过 分析风电场实际测风的原始资料得到。根据风电场实际测风的结果，假设风速是以小时平均，按每小时正点前十 分钟测取，那么在一年之内就有N个测点，这样可得风电场实际的风频分布为:F 土（2-1）i N式中F风速：,的实际分布频率;Nvi年内风速V仅出现的次数;Ny风电场风速符合威布尔分布:年内总的测风点数,一般有N广8760。K-1e(2-2)式中，v为风速（m/s），f（）为威布尔分布函数，A、K为威布

7、尔尺度系 数（m/s）和形状系数。利用风电场测风的结果，对实际所得的风速数据进行统计，得出年平均风速Vp和风速频率分布F，并采用最小逼近法，y |fi - F |2 = min算出威布尔分布参数A、K的近似值。从而得到风速风频特性的数学模型， 进而得到风电场风能资源分布和评估、风力发电机组选型和发电量的预测以及 风电场并网对系统的影响分析。2.2风速模型通常用四种成分的风速来模拟实际风速：基本风V仍、阵风V“g、渐变风 和随机风。（1）基本风Vwb基本风反映了风场平均风速的变化，风力发电机向电网输送功率的大小主 要由基本风决定，它的测得由风电场测风所得的威布尔分布参数近似确定。一 般认为基本风

8、在一段时间内不随时间变化，可取常数。11Vwb -r1 + 京（2-4）Vwbt图2-1基本风随时间变化曲线图（2）阵风vwg阵风为描述风速突然变化的特性，可假设在该段时间内风速具有余弦特性。V = wgGmax21 - cos 2 兀(T* Tg0Tg t Tig + Tg其他时间(2-5)式中，G -max阵风幅值（m/s ）;Tg 阵风周期（s）；T1G图2-2阵风随时间变化曲线图（3）渐变风vwr渐变风用以描述风场稳态能量随时间缓慢变化的过程，以风速由小变大为 例，渐进风可用下式模型：式中，vwrt - T1RT2 R-，RRmaxRmaxt T2R + TrTr t T2RTr t

9、Tr + Tr(2-6)R渐变风的最大值;maxT1R 渐变风开始时刻；T2 R 渐变风结束时刻;渐变风保持时间。图2-3渐变风随时间变化曲线图（4）随机风随机风表示风速变化的随机特性：随机噪声风速。v =翌 k （w kw cos（w t + 中）i=1S (w )=V i(2-7)w, = （i- 0.5） Aw式中，w,第i个分量的角频率；Aw 随机分量的离散间距；中一一在02兀间服从均匀概率密度的随机变量; iKN地表粗糙系数，一般取0.004；F动范围（m 2 ）；相对高度的平均风速（m/s ）;七（.）风速随机分量分布谱密度（m 2 /s ），通过对其积分便可得短期 风速数据。（5

10、）合成风速模拟实际作用在风力机上的风速为：（2-8）（6）综合风速模型w=Vwb + Vwg +，+，图2-4综合风速模型输入参数如下： 基本风： v = 9m / s。wb 阵风：V = 2m / s，Tg = 3s，气=1s，数量为1。 渐变风：v = 2m / s， Tr = 4s， Tr = 1s，数量为 1.5. 随机风： 匕=0.004，F = 2000m2，n = 50。仿真结果如下：风速模型输出050505050 2522017.偶12.10.75. Vw:.俳-年lullMikU 闻 hj ：山 h：I 北航 HL0.05.010015.020.025.030.0图2-5综合

11、风速模型仿真结果在前面我们已经讨论过，风是近似的服从威布尔分布，也就是说，近似的 服从正态分布。如图2-5所示，在没有外力风速的情况下，由于受随机噪声风的 影响，风速的曲线波动很大，在3s和4s时分别又受到阵行风与渐变风的影响， 波形也出现了相应的波动，其综合风速的最大值可达到15.96m/s。所示说，用 以上的四个风的分量在一定的程度上是可以大体的描述风的波形，但在一些细 节上还需要进一步修正，所以它的使用范围是有限的，只是可以用在一些要求的精确程度不高的模型的仿真。2.3风力机建模与分析2.3.1风力机能量转换过程风力机能量转换模型的功率及转矩计算公式是根据流体力学中气流的动能 计算公式，

12、并结合贝兹理论得到的，详见资料4。风力机简化模型如下：风七TmE图2-6风力机简化模型风力机，风能的吸收和转换装置。传动装置主要包括轮毂、齿轮箱和传动 轴，起连接和传动作用。发电机，能量转换装置。在变桨距风机中还应包括桨 距角控制环节。能量转换过程是：风能一机械能一电能。由文献6得，风力机轴上的输出机械功率为：P = - pnR2v3C （k, P）（2-9）w 2P式中，P空气密度（kg /m3 ）;R 风机叶轮半径（m ）;k 叶尖速比，定义为k =捋十，其中出为风力机叶轮转速（rad/s ）， eqV为等效风速（m/s ）；P桨距角（deg ）;cp 风能利用系数，是叶尖速比k和叶片桨距

13、角p的函数；对于给定的风力机系统，cp的表达式是一定的。一种变桨距风力机的风能转化效率系数：C (k, p)= 0.22 p哄-0.4p- 5、I 8）(2-10)1 _10.035=8 人 + 0.08。 P 3 +1风力机获得转矩为:w Otur(2-11)定义ct Q,。）为转矩系数，C 炽,p）= CP。）T人注：iRo由人=_产推出eqCO= opt vtur R eq(2-12)(2-13) 对于给定的叶片桨距角。，不同的叶尖速比所对应的Cp值相差较大; 对于给定的。，有且仅有一个固定的x = xop能使Cp达到最大值； 在风速不断变化的情况下,要保持x列、o 必须随着风速按照孑的

14、比例 变化，才能保证风力机捕获的风能最大、效率最高。这是采用变速风电机组代替固定转速风电机组的初衷之一。图2-7风机 入特性曲线P对于变桨距型风力发电机组，cp特性可近似表示为：C = 0.5IRCL - 0.022P - 2：-0-255筌（2-14）P I人J式中，cf为叶片设计常数，一般取1-3。2.3.2风力机的稳定工作区空间曲面虽然能包含风力机运行的所有状态点，但是对于分析不太方便， 所以在实际应用中多是取几个离散的P值，画出平面图的方法，如下图所示，取6组P值，绘制如下：图2-5以的网曲线图25代瓯状）曲线图2-8风力机稳定工作区曲线图在T侦,v）曲线中，以转矩T的最大值为顶点连成的一条线AB，将曲线簇分 成了两部分，其中右侧为稳定运行区域，左侧部分为不稳定工作区域。对比 P效v）和T侦,v）曲线，我们发现当T达到最大时，P并没有达到最大，具体而言， 就是最大功率点对应的转速值要大于最大转矩点的转速值。如图2-6所示，曲线 CD是由最大转矩点的连线而成的，曲线EF则是由最大功率