风力发电系统建模与仿真

风力发电系统建模与仿真摘要：风力发电作为一种清洁的可再主能源利用方式，近年来在世界范围内获得 了飞速的发展。本文基于风力机发电建立模型，主要完成了以下工作：（ 1） 基于风资源特点，建立了以风频、风速模型为基础的风力发电理论基础；（ 2） 运用叶素理论，建立了变桨距风力机机理模型；（ 3） 分析了变速恒频风力发电机的运行区域与变桨距控制的原理与方法，并给出了机组的仿真模型，为风力发电软件仿真奠定了基础；（ 4） 搭建了一套基于PSCAD/EMTDC 仿真软件的风力发电系统控制模型以及 完整的风力发电样例系统模型，并且已初步实现风力机特性模拟功能。关键词：风力发电；风频；风速；风力机；变桨距；建模与仿真1 风资源及风力发电的基本原理1.1 风资源概述（ 1） 风 能的基本情况山风的形成乃是空气流动的结果。风向和风速是两个描述风的重要参数。风向 是指风吹来的方向，如果风是从东方吹来就称为东风。风速是表示风移动的速度 即单 位时间内空气流动所经过的距离。风速是指某一高度连续 lOmin 所测得各瞬时风速的平均值。 一般以草地上空 10m高处的 lOmin 内风速的平均值为参考。风玫瑰图是一个给定地点一段时间内 的风向分布图。通过它可以得知当地的主导风向。风能的特点主要有：能量密度低、不稳定性、分布不均匀、可再生、须在有风地带、无污染、 分布广泛、 可分散利用、另外不须能源运输、可和其它能源相 互转换 等。（ 2） 风能资源的估算风能的大小实际就是气流流过的动能，因此可以推导出气流在单位时间内垂直流过单位截面积的风能，即风能密度，表示如下：co = 0.5/cn3（ 11 ）式中，3风能密度（ W/m2） , 是描述一个地方风能潜力的最方便最有价值的量；p空气密度(匕 / 用)；V 风速(加/s) O由于风速是一个随机性很大的量，必须通过一段时间的观测来了解它的平均 状况，一个地方风能潜力的多少要视该地常年平均风能密度的大小。因此需要求 出在一段时间内的平均风能密度，这个值可以将风能密度公式对时间积分后平均 来求得。有效风能密度还可根据下式求得Q = 0.5pv*P(y)dv( 1-2)式中，儿 - 启 动风速( mis ) ；v, 停 机风速( $) ；P(v) 有效风速范围内的条件概率分布密度函数 。平均风能密度则可用下式求得：o = P(v)clt( 1-3)1.2 风力发电的基本原理风能具有一定的动能，通过风轮机将风能转化为机械能，拖动发电机发电。 风力发电的原理是利用风带动风车叶片旋转，再通过增速器将旋转的速度提高来促使发电机发电的。依据H 前的风车技术，大约 3m/s 的微风速度便可以开始发 电。风力发电的原理说起来非常简单，最简单的风力发电机可山叶片和发电机两 部分构成如图 1-1 所示。空气流动的动能作用在叶轮上，将动能转换成机械能， 从而推动片叶旋转，如果将叶轮的转轴与发电机的转轴相连就会带动发电机发出 电来。八叶片风向发电机V.图1-1风力发电原理图1.3 风力发电的特点风力发电具有以下特点：可再生的洁净能源；建设周期短，装机规模灵活,可根据资金悄况决定一次装机规模，有一台资金就可以安装一台投产一台； 可靠性高，把现代高科技应用于风力发电机组使其发电可靠性大大提高,中、大型风力发电机组可靠性从 80年代的50%提高到了 98%,高于火 力发电且机组专命可达20年； 造价低，运行维护简单，实际占地面积小； 发电方式多样化，既可并网运行，也可以和其他能源如柴油发电、太阳能发电、水利发电机组形成互补系统，还可以独立运行；单机容M小2风能及风力机系统模型的建立2.2风频模型风速具有明显的随机性和间歇性。为了较精确地描述风速及其变化特性，引 入 风频分布的概念。风频分布就是风速的统汁概率分布，是衡M风能资源分布特性的重要指标，它反映了风电场莫个时段每一风速出现的概率，可以通过分析风电场实际测风的原始资料得到。根据风电场实际测风的结果，假设风速是以小时平均，按每小时正点前十分钟测取，那么在一年之内就有 N个测点，这样可得风电场实际的风频分布为：(2-1)式中F, 风速的实际分布频率；Nvi一年内风速出现的次数；N、 . 一年内总的测风点数，一般有Ny = 8760 o风电场风速符合威布尔分布：J (2-2)式中， 为风速 ( m/s), /(?) 为威布尔分布函数， A、 K 为威布尔尺度系 数(加/s)和形状系数。利用风电场测风的结果，对实际所得的风速数据进行统讣，得出年平均风速匕和风速频率分布片，并采用最小逼近法，% 2工 |/T =min 日(2-3)算出威布尔分布参数 A 、 K 的近似值。从而得到风速风频特性的数学模型, 进而得到风电场风能资源分布和评估、风力发电机组选型和发电量的预测以及风 电场 并网对系统的影响分析。2. 2 风速模型通常用四种成分的风速来模拟实际风速：基本风八阵风 R 、渐变风和随机风。(1)基本风基本风反映了风场平均风速的变化，风力发电机向电网输送功率的大小主要山基本风决定，它的测得山风电场测风所得的威布尔分布参数近似确定。一般认 为基 本风在一段时间内不随时间变化，可取常数。%=A? r 1+舟(2-4)(2)阵风图2J基本风随时间变化曲线 图阵风为描述风速突然变化的特性，可假设在该段时间内风速具有余弦特性。I Gmax 11 -cos2/r T g t T?r + TrTir t T 2r ( 2-6)T?R t T2R + Tr式中，/? max一渐变风的最大值R 渐变风开始时刻；T2R 一渐变风结束时刻Tk渐变风保持时间。图2? 3渐变风随时间变化曲线（4）随机风图随机风表示风速变化的随机特性：随机噪声风速。HI心=2 工Sv( H ; ) Avvl cos(vv/ + 0) i=is后)二咨严帆(2-7)w;. = （i - 0.5） Avv式中，叫一一笫i个分M的角频率；Aw 随机分M的离散间距；5在02兀间服从均匀概率密度的随机变M ;心地表粗糙系数，一般取 0.004;F 一 扰动范围（nr ）;“一一相对高度的平均风速（m/s ）；?（叫）一一风速随机分M分布谱密度（亦“），通过对其积分便可得短期 速数据。(5)合成风速 模拟实际作用在风力机上的风速为：%二人临+ 峡+ %, +人H(6)综合风速模型Wind Source Mean rES GUS 一 Ramp L-dNoise.图2-4综合风速模型 输入参数如下：(2-8)D 基本风： = 9m/s o 阵风：Vmax = 2m/ s , T ig = 3s , Tg =s 9 数M为 1。D 渐变风：Vmax = 2in/s , TR=4s ,TR =s ,数M 为 1.5. 随机风：Kn = 0.004 , F = 2000 沪,仿真结果如下:风逋模型输出图2-5综合风速模型仿真结果在前面我们已经讨论过，风是近似的服从威布尔分布，也就是说，近似的服从正态分布。如图2-5所示，在没有外力风速的情况下，山于受随机噪声风的影响，风速的曲线波动很大，在3s和4s时分别乂受到阵行风与渐变风的影响，波形也出现了相应的波动，其综合风速的最大值可达到15.96m./so所示说，用以 上的四个风的分量:在一定的程度上是可以大体的描述风的波形，但在一些细节上还需要进一步修正，所以它的使用范用是有限的，只是可以用在一些要求的精确程度不高的模型的仿真。2.3风力机建模与分析2. 3.1风力机能M转换过程风力机能M转换模型的功率及转矩计算公式是根据流体力学中气流的动能讣算公式，并结合贝兹理论得到的，详见资料4。风力机简化模型如下：几E风速图2? 6风力机简化模型风力机，风能的吸收和转换装置。传动装置主要包括轮毂、齿轮箱和传动轴,起连接和传动作用。发电机，能M转换装置。在变桨距风机中还应包括桨距角控制环节。能M转换过程是：风能f机械能f电能。山文献6得，风力机轴上的输出机械功率为：化=才欣2护Cp （入0）（2-9）式中，p空气密度（kg/心；R风机叶轮半径5 ）;2叶尖速比，定义为2 =一 ，其中为风力机叶轮转速（rad/s ） 9%怙为等效风速S/s）;P桨距角（deg）;Cp风能利用系数，是叶尖速比兄和叶片桨距角0的函数；对于给定的风力机系统，Cp的表达式是一定的。一种变桨距风力机的风能转化效率系数Cp (入 0)= 0? 222.5-0.4/7-5 e(2-10)%110.0357 2 + 0.08/7+1风力机获得转矩为(2-11)定义皿“)为转矩系数注:R%推出(2-12)对于给定的叶片桨距角AopiR 10,不同的叶尖速比所对应的(2-13)Cp值相差较大；对于给定的0,有且仅有一个固定的 入=心能使Cp达到最大值；;&0.50246810121416对、探骐罟隹比例R变化，才能保证风力机捕获的风能最大、效率最高。这是采用变速风电机组代替固定转速风电机组的初衷之一。图2? 7风机Cp-A特性曲线对于变桨距型风力发电机组，Cp特性可近似表示为(2-14)A-0.022A-2 L=0.5式中，？为叶片设计常数，一般取13。2. 3. 2风力机的稳定工作区空间曲面虽然能包含风力机运行的所有状态点，但是对于分析不太方便，所以在实际应用中多是取儿个离散的 0值，画出平面图的方法，如下图所示，取 6组”值，绘制如下：图24以电力曲线图2-图2-8风力机稳定工作区曲线图在7 （八,v）曲线中，以转矩T的最大值为顶点连成的一条线 AB,将曲线簇分 成 了两部分，其中右侧为稳定运行区域，左侧部分为不稳定工作区域。对比P （o）和几力）曲线，我们发现当 T达到最大时，P并没有达到最大，具体而 言，就是最大功率点对应的转速值要大于最大转矩点的转速值。如图 2-6所示，曲线 CD是山最大转矩点的连线而成的，曲线EF则是山最大功率点的连线而成。这 样一来，在PM曲线簇中，CD曲线和EF曲线之间的部分也是稳定区域。2. 3. 3基于叶素理论的风力机建模基于叶素理论的风轮建模是将叶片分为若干个微元，称为叶素，通过对叶素的受力分析求得微元转矩，再将所有微元转矩相加得到风力发电机风轮的输出转矩15-751，在风轮半径r处取一长度为的叶素，其弦长为/,节距角为0。图2? 9叶素微元受力分析图(2-16)(2-15)的作用。很小时，风轮(2-17)(2-18)如上图所示，来流方向的风速为h在半径为/?处的风轮机速度为？尸(?为风轮机角速度 八 气流相对于叶片的相对速度为W,则有:w = v-u叶素力 在相对速度为w的气流作用下，受到一个方向斜向上的气动力 将dF沿与相对速度w垂宜及水平方向可分解为升力 d Z4口阻力刃儿 当 可以近似的将叶素面积看成弦长与叶素长度的乘积，可得如下计算公式:dL = -pC flw2drdD = ApCdlw2dr气动力dF按垂宜和平行于旋转平面方向分解为法向力巧和切向力几, 转矩山切向力刃；产生，则有转矩微元 ：dT = r ?dF(i =)dL -sin I - dD ? cos I)令 = CdICA得总转矩汁算公式：RR rr =小 | dr =n ? J Q /W? G ? sin /(1矣? cot %式中，心一一轮毂半径；H风轮包含的桨叶个数；/ 倾斜角(桨距角P与攻角，之和)。上式为基于叶素理论的风轮模型函数，可以写为如下形式：T = fAu.PY 即风轮输出转矩为风速-风轮转速，桨距角0的函数。2. 3. 4基于PSCAD风力机模型与仿真基于PSCAD的风力机模型如下Beta ： Beta图2? 10风力机简化模型 基于上面的原理及理论公式，仿真结果如下：Title (for display on plots.meters.)GroupDisplay title on icon?Scale FactorUnit (for display on meters;)Save output during Multiple Run onWind Turbine MOD 2 Type7、一 广/Default Min / Max Limits p图2? 11风轮机机械转矩输出及参数设置图2J2风轮机机械转矩模拟仿真结果Vw12 010.08.0 604 0, 2.000h心 3 ，血 血 ado25。30 0IVbin : GraphsTmNoITopu | La st Run Only RTTitis (for display on plots(metersr.)GroupDis0 12Vtit拒 on icon?Scale FactorUnit (for display on meters: p,apkA,)S9ve output during Multiple Run onDefault Min / Max Limits W图2-13风轮机机械功率输出及参数设置PTNo二|ro|puLastRun Only jrIVbin : Graphs图2-14风轮机机械功率模拟仿真结果由上述两图曲线可知，风轮机的输出转矩和输出功率都是标幺值，则它们的曲线是完全一致的，在0、3s时变桨距控制系统在调节桨叶节距使转矩和功率输 出逐渐达到 稳定，山于乂突然受到在 3s与4s分别受到阵行风与渐变风的影响， 从而使波形在这两 个时间有的突变，之后继续达到稳定。3变桨距风力发电机组控制系统模型3. 1变桨距风力发电机组的运行状态变桨距风力发电机组根据变距系统所起的作用可分为三种运行状态，即风力 发电机组的起动状态（转速控制）、欠功率控制（不控制）和额定功率状态（功率控制）。由于变桨距系统的响应速度受到限制，对快速变化的风速，通过改变节距来控制输出功率的效果并不理想。因此，为了优化功率曲线，最新设计的变桨风力发电机组在进行功率控制的过程中，其功率反馈信号不再作为宜接控制桨叶节距的变变桨距系统山风速低频分M和发电机转速控制，风速的高频分M产生的机械能波动，通过迅速改变发电机的转速来进行平衡，即通过转子电流控制器对发电机转差率进行控制，当风速高于额定风速时，允许发电机转速升高，将瞬变的风能以风轮动能的形式储存起来；转速降低时，再将动能释放出来，使功率曲线达到理想的状态。3.2变桨距控制系统图3-1变桨距控制系统工作原理图在发电机并入电网时前，发电机转速山速度控制器 A根据发电机转速反馈信 号与给 定信号宜接控制；发电机并入电网后，速度控制B与功率控制器起作用。 功率控制器的任务主要是根据发电机转速给出相应的功率曲线，调整发电机转差率，并确定速度控制器B的速度给定。节距的给定参考值山控制器根据风力发电机组的运行状态给出。如图3-1所示，当风力发电机组并入电网前，山速度控制器 A给出；当风力发电机组并入 电网后山速度控制 B给出。(2)变距控制变距控制系统是一个随动系统，如图2-9所示。变距控制器是一个非线性比例控制器，它可以补偿比例阀的死带和极限。变距系统的执行机构是液压系统，节距控制器的输出信号经D/A转换后变成电压信号控制比例阀(或电液伺服阀),驱动液压缸活塞，推动变桨距机构，使桨叶节距角变化。活塞的位移反馈信号由位移传感器测经转换后输入比较器。图3-2变距控制结构框图(3)速度控制器A转速控制器A在风力发电机组进入待机状态或从待机状态重新起动时投入工作，如图2-10所示在这些过程中通过对节距角的控制，转速以一定的变化率上升。控制器也用于在同步速（50Hz时1500转/min）时的控制。当发电机转速在同步转速0min 内持续Is发电机将切入电网。控制器包含着常规的PD控制器和PI控制器，接着是节距角的非线性化环节,通过 非线性化处理，增益随节距角的增加而减小，以此补偿山于转子空气动力学产生的非线性，因为当功率不变时，转矩对节距角的比是随节距角的增加而增加的。当风力发电机组从待机状态进入运行状态时，变桨距系统先将桨叶节距角快速地转到45,风轮在空转状态进入同步转速。当转速从 0增加到1500min时，节距角给定 值从45线性的减小到5。这一过程不仅使转子具有高起动力矩，而且在风速快速地增大时能够快速起动。（ 4） 速度控制器B发电机切入电网后，速度控制系统 B 作用。如图 2-11 所示，速度控制器B 受发电机转速和风速的双重控制。在达到额定值前，速度给定值随功率给定值按 比例增加。额定的速度给定值是1569r/min, 相应的发电机转差率是 4%。如果风 速和功率输出一直低于额定值，发电机转差率将降低到 2%,节距控制将根据风 速调整到最佳状态，以优化叶尖速比。如果风速高于额定值，发电机转速通过改变节距来跟踪相应的速度给定值。 功率输出将稳定地保持在额定值上。从图中可知在风速信号输入端设有低通滤波 器，节距控制对瞬变风速并不响应。4风力发电控制系统的模拟仿真分析4.1 无穷大系统模型的建立图4-1风力发电机无穷大系统模型-HTQ 口00tN0-0004.2 风力发电机系统并网模拟仿真分析(1)发电机三相电东输钳及仿真结果分析Tifle (for display on plots,meterGroupDisplay title on icon?Scale FactorUnit (for display on meters: p,j.M -JEave output during Multiple Run onEbNo1,0Last F?un Only Default Min/Max Limits Lji图4-2异步发电机发电机三相电压输出及参数设 宜iVfain : Graphs旦册Fvein : GraphsU7Q5；O5Q7Q5Q0.05.010.015.020.025.030.0图4-3异步发电机发电机三相电压模拟仿真结果IVhn: Graphs6.0图4-4高压母线电压模拟仿真结果异步发电机三相电压模拟仿真结果分析如图4-3所示，在发电机并入电网前低压侧电压为0,在Is时发电机并入电网，电压0.50突然上升到接近于额定电压，发电机向电网输送功率。并网后的功率输出 逐渐趋于稳定,这时的电压曲线呈正弦曲线形状变化，向电网输送额定功率。(2)低压母线和高压母线的线电压输出及仿真结果分析2.00 -1it 11 t0.02.04.06,08.010.012.014.016.018,020.0Main : Graphsc Eab(121)T60 rI一I1一一tI一Ir100 -I I I I I I II I III II I I Illi I II t 0.0204.06.08010 012.014.016.018 020 0图4-5低压母线和高压母线的线电压输出及仿真结果在正常运行时，低压母线和高压母线电压均从0迅速上升并均达到各自的额定值,然后一宜保持稳定。低压母线电压稳定在0.69EU左右，高压母线电压稳 定在121 it V左(3)低压母线相电流输出及仿真结果分析并网前电流为0,在Is时断路器合闸并网，出现很大的冲击电流，其冲击电 流值达到11KA,最后开始衰减至0. 07KA,然后乂开始上升，最后趋于稳定，具电流最大稳定值为0. 64KAoMain : Graphs_ Eab(069)ZUO 1III111IF0.500.00Title (for display on plots.meters,)11bft 延rGroupDisplay title on icon?Scale FactorUnit (for display on meters; p.u.(kA ()Save output durinq Multiple Run onDefault Min! Max Limits p图4-6低压电流输出及参数设置Ivbin : Graphs-4.0-6.0?8.00.05.010.015.020.025.030.00.05.010.015.020.025.030.0Main : Graphs9.06.0300.0-3.00.05.010.015.020.025.030.0图4-7低压母线相电流输出及仿真结果4. 3变桨距控制系统模拟仿真分析通过控制桨距角的大小的改变就可以控制叶片吸收风功率的多少，桨距角的调节可以使发电机输出功率平稳。桨叶节迎角图4-8变桨距控制系统模拟仿真结果变桨距控制系统模仿真结果分析如下：风轮机启动时风力发电机组开始自动运行于风轮叶尖本来值 90,即桨矩 角初始值 为90度，在机组起动的过程中逐渐变小，这样叶片吸收风能逐渐增大 ,叶片的转速也逐 渐加快，最后在1.4s时桨矩角变为零，且保持不变，此时叶片 吸收风能达到了最大。5结论本文通过PSCAD/EMTDC电力系统模拟仿真软件，建立了变桨距风力发电机组控制系统模型，对加入控制系统的风力发电样例系统进行模拟仿真分析，验证了 控制系 统模型的可用性。风力发电系统控制策略是以风速的变化为依据，风能的最大利用效率为目的，为优化风力发电系统运行特性提出的控制方案。变桨距控制系统的设计主要采用 PI控制器，根据发电机有功功率输出和风轮机转速反馈来调节桨叶节距。通过风轮机桨距角控制系统对叶片桨距角进行控制，使风力发电机组的机械部分与发电机的电气部分配合，达到提高风能利用效率及改善供电质量:的LI的。利用风力发电样例系统来验证控制系统的可用性，并对各种仿真曲线进行分析，从而得出结论。根据风速模型的仿真曲线，分析风轮机和发电机各部分曲线的变化情况和整个系统的仿真曲线图。在并网以前电压的波形基本上是正弦形状的，转速基本上是稳定的。并网以后虽然受到了电网的干扰，但转速上升到额定转速后再没有多大变化；电流的波形虽然是正弦的，但整体的趋向也发生了相应的波动。变桨距控制系统在风力发电机组起动时，通过变距来获得足够的起动转矩；起动以后，当低于额定风速运行时的机组状态控制为转速，当高于额定风速运行时，通过调 整桨叶节距，改变气流对叶片的攻角，从而改变风力发电机组获得的空气动力转 矩，使功率输出保持稳定。额定风速之后的机组状态控制主要山桨距角调节实现。 得到的控制系统保持了风力发电机组运行的安全可靠性。参考文献1黄素逸.能源与节能技术.北京中国电力出版社.20042 叶杭冶 . 风力发电机组的控制技术. 机械工业出版社. 20023 宫靖远 . 风电场工程技术手册. 机械工业出版社. 20044 王承煦 张源 . 风力发电 . 中国电力出版社. 20025 李俊峰 . 风力 12 在中国 . 北京化学工业出版社. 20056 赵凤山 . 风力发电论文集. 金盾出版社. 20027 风电发展情况调研组， 2009年中国风电发展惜况调研报告20091, 28 姚兴佳，宋俊等，风力发电机组原理与应用，北京：机械工业出版社， 20089 高峰，风力发电机组建模与变桨距控制研究： 博士学位论文 , 北京：华北 电力大学， 200910 李东东，风力发电机组并网与仿真分析，水电能源科学， 2006年，第 24 卷第 1期. P： 1-311 何祚麻王亦楠. 我国新能源可持续发展的现实选择 . 中国三峡建设.2005. P： 25 2712 肖运启 . 双馈型风力发电机组励磁控制与优化运行研究： 博士学位论文 北京：华北电力大学， 2008.13 赵栋利，胡书举，赵斌，许洪华. 风力发电机、变流器及其低电压穿越概述 . 变频器世界， 2009, 2.14 张仰飞，袁越，陈小虎 . 风力机参数的可辨识分析 . 电力系统自动化， 2009, 33(6)： 86-89.15 S.M. B. Wilmshurst. Control strategies for Wind turbines. Wind Engineering. 1998 Vol. 12: P 236一一 24916 R. Chedid, F? Mrad and H? Basma? Intelligent Control for Wind Energy ConversionSystem. Wind Eng? 1998, Vol? 18:No. 1, 1 16