风力发电系统建模与仿真.docx

风力发电系统建模与仿真摘要： 风力发电作为一种清洁的可再生能源利用方式， 近年来在世界范围内获得 了飞速的发展。本文基于风力机发电建立模型，主要完成了以下工作：(1)基于风资源特点，建立了以风频、风速模型为基础的风力发电理论基 础；(2)运用叶素理论，建立了变桨距风力机机理模型；(3)分析了变速恒频风力发电机的运行区域与变桨距控制的原理与方法， 并给出了机组的仿真模型，为风力发电软件仿真奠定了基础；(4)搭建了一套基于 PSCAD/EMTDC 仿真软件的风力发电系统控制模型以及 完整的风力发电样例系统模型，并且已初步实现风力机特性模拟功能。关键词：风力发电；风频；风速；风力机；变桨距；建模与仿真1 风资源及风力发电的基本原理1.1 风资源概述(1)风能的基本情况1风的形成乃是空气流动的结果。 风向和风速是两个描述风的重要参数。 风向 是指风吹来的方向， 如果风是从东方吹来就称为东风。 风速是表示风移动的速度 即单位时间内空气流动所经过的距离。风速是指某一高度连续 10min 所测得各瞬时风速的平均值。一般以草地上空 10m 高处的 10min 内风速的平均值为参考。 风玫瑰图是一个给定地点一段时间内 的风向分布图。通过它可以得知当地的主导风向。风能的特点主要有：能量密度低、不稳定性、分布不均匀、可再生、须在有 风地带、无污染、分布广泛、可分散利用、另外不须能源运输、可和其它能源相 互转换等。(2)风能资源的估算风能的大小实际就是气流流过的动能， 因此可以推导出气流在单位时间内垂 直流过单位截面积的风能，即风能密度，表示如下： = 0.5pv3 (1- 1)式中, 风能密度(W / m2)， 是描述一个地方风能潜力的最方便最有价值的量；p 空气密度( kg / m3)；v 风速( m / s )。由于风速是一个随机性很大的量， 必须通过一段时间的观测来了解它的平均 状况， 一个地方风能潜力的多少要视该地常年平均风能密度的大小。 因此需要求 出在一段时间内的平均风能密度， 这个值可以将风能密度公式对时间积分后平均来求得。有效风能密度还可根据下式求得 = j v2 0.5pv3 P(v)dvv1式中，(1-2)v 启动风速( m / s )；1v 停机风速( m / s )；2P(v) 有效风速范围内的条件概率分布密度函数 2。平均风能密度则可用下式求得： = 1 j 0.5pv3P(v)dt (1-3)T1.2 风力发电的基本原理风能具有一定的动能，通过风轮机将风能转化为机械能，拖动发电机发电。 风力发电的原理是利用风带动风车叶片旋转， 再通过增速器将旋转的速度提高来 促使发电机发电的。依据目前的风车技术，大约 3m/s 的微风速度便可以开始发 电。 风力发电的原理说起来非常简单， 最简单的风力发电机可由叶片和发电机两 部分构成如图 1- 1 所示。空气流动的动能作用在叶轮上，将动能转换成机械能， 从而推动片叶旋转， 如果将叶轮的转轴与发电机的转轴相连就会带动发电机发出电来。1.3 风力发电的特点风力发电具有以下特点： 可再生的洁净能源； 建设周期短，装机规模灵活，可根据资金情况决定一次装机规模，有一 台资金就可以安装一台投产一台； 可靠性高， 把现代高科技应用于风力发电机组使其发电可靠性大大提高， 中、大型风力发电机组可靠性从 80 年代的 50%提高到了 98%，高于火 力发电且机组寿命可达 20 年； 造价低，运行维护简单，实际占地面积小； 发电方式多样化，既可并网运行，也可以和其他能源如柴油发电、太阳 能发电、水利发电机组形成互补系统，还可以独立运行； 单机容量小2 风能及风力机系统模型的建立2.2 风频模型风速具有明显的随机性和间歇性。 为了较精确地描述风速及其变化特性， 引 入风频分布的概念。 风频分布就是风速的统计概率分布， 是衡量风能资源分布特 性的重要指标， 它反映了风电场某个时段每一风速出现的概率， 可以通过分析风 电场实际测风的原始资料得到。根据风电场实际测风的结果， 假设风速是以小时平均， 按每小时正点前十分 钟测取，那么在一年之内就有 N 个测点，这样可得风电场实际的风频分布为：F iNvi(2- 1)Ny式中F 风速 v 的实际分布频率；i wiN 一年内风速 v 出现的次数；vi wiNy 一年内总的测风点数，一般有 Ny = 8760。风电场风速符合威布尔分布：f (v)= K (| V )|K -1 e-(|(A(V)|KA ( A )式中， v 为风速( m /s )， f (.)为威布尔分布函数，(2-2)A、 K 为威布尔尺度系数( m /s )和形状系数。利用风电场测风的结果， 对实际所得的风速数据进行统计， 得出年平均风速V 和风速频率分布 F ，并采用最小逼近法，P iNxy f - F 2 = mini ii=1(2-3)算出威布尔分布参数 A、 K 的近似值。从而得到风速风频特性的数学模型， 进而得到风电场风能资源分布和评估、风力发电机组选型和发电量的预测以及风 电场并网对系统的影响分析。2.2 风速模型通常用四种成分的风速来模拟实际风速：基本风v 、阵风v 、渐变风和随wb wg机风。(1) 基本风 vwb基本风反映了风场平均风速的变化， 风力发电机向电网输送功率的大小主要 由基本风决定， 它的测得由风电场测风所得的威布尔分布参数近似确定。 一般认 为基本风在一段时间内不随时间变化，可取常数。vwb = A . T 1 + K(1) (2-4)vwbt图 2- 1 基本风随时间变化曲线图(2) 阵风 vwg阵风为描述风速突然变化的特性，可假设在该段时间内风速具有余弦特性。(2-5)vwg = G 2(max) (|(t - T1GTG)| T1G 其( t)他()T1时 TG式中，G 阵风幅值( m /s )；m a xT 阵风周期( s )；GT 阵风开始时刻( s )。1G图 2-2 阵风随时间变化曲线图(3) 渐变风 vwr渐变风用以描述风场稳态能量随时间缓慢变化的过程，以风速由小变大为 例，渐进风可用下式模型：v = 1R R0t-Twr T2R ax(R) maxt T + T1R 2R RT t T1R 2RT t T + T2R 2R R(2-6)式中，R 渐变风的最大值；m a xT 渐变风开始时刻； T 渐变风结束时刻；1R 2RT 渐变风保持时间。R图 2-3 渐变风随时间变化曲线图(4) 随机风随机风表示风速变化的随机特性：随机噪声风速。v = 2xn S (w )Aw2(1) cos(w t +v )wn V i i ii=12K F 2 wN iS (w )= 2 1 + (|( A(F)(w)i )|2 3(4)V iw = (i - 0.5)Awi式中，(2-7)w 第 i 个分量的角频率；iAw 随机分量的离散间距；v 在 0 2 间服从均匀概率密度的随机变量； iK 地表粗糙系数，一般取 0.004；NF 扰动范围( m 2 )；A 相对高度的平均风速( m / s )；S (w ) 风速随机分量分布谱密度( m2 / s )，通过对其积分便可得短期V i风速数据。(5) 合成风速模拟实际作用在风力机上的风速为：(2-8)v = v + v + v + vw wb wg wr wn(6)综合风速模型Wind SourceMeanES Gust Vw Ramp VwNoise图 2-4 综合风速模型输入参数如下： 基本风： v = 9m /s。wb 阵风： v = 2m /s， T = 3s， T = 1s ，数量为 1。max IG G 渐变风： v = 2m /s， T = 4s， T = 1s ，数量为 1.5. max IR R 随机风： K = 0.004， F = 2000m2， n = 50。 N仿真结果如下：图 2-5 综合风速模型仿真结果在前面我们已经讨论过， 风是近似的服从威布尔分布， 也就是说， 近似的服 从正态分布。如图 2-5 所示,在没有外力风速的情况下，由于受随机噪声风的影 响， 风速的曲线波动很大， 在 3s 和 4s 时分别又受到阵行风与渐变风的影响， 波 形也出现了相应的波动，其综合风速的最大值可达到 15.96m/s。所示说，用以 上的四个风的分量在一定的程度上是可以大体的描述风的波形， 但在一些细节上 还需要进一步修正， 所以它的使用范围是有限的， 只是可以用在一些要求的精确 程度不高的模型的仿真。2.3 风力机建模与分析2.3.1 风力机能量转换过程风力机能量转换模型的功率及转矩计算公式是根据流体力学中气流的动能 计算公式，并结合贝兹理论得到的，详见资料4。风力机简化模型如下：风速vw风力机bE( .)发电机Og传动装置桨距角控制T aeT m图 2-6 风力机简化模型风力机， 风能的吸收和转换装置。 传动装置主要包括轮毂、 齿轮箱和传动轴， 起连接和传动作用。 发电机， 能量转换装置。 在变桨距风机中还应包括桨距角控 制环节。能量转换过程是：风能 机械能 电能。由文献6得，风力机轴上的输出机械功率为：P = 1 pR 2 v3 C (入, b) w 2 P式中，(2-9)p 空气密度(kg / m3 )； R 风机叶轮半径( m )；入 叶尖速比，定义为入 = ROtur，其中 O 为风力机叶轮转速( rad / s )， v tureqv 为等效风速( m /s )；eqb 桨距角( deg )；C 风能利用系数，是叶尖速比 入 和叶片桨距角 b 的函数；对于给定的P风力机系统， C 的表达式是一定的。P一种变桨距风力机的风能转化效率系数：CP (入, b)= 0.22(|(16(1)6 - 0.4b - 5)|e-16(2.)51 1 0.0356 = 入 + 0.08b - b 3 + 1风力机获得转矩为：T = Pww Otur(2- 10)(2- 11)定义 C (入, b)为转矩系数，TC (入, b)= C (入, b)PT 入注：由入 = Roturveq推出o = 入optvtur R eq(2- 12)(2- 13) 对于给定的叶片桨距角 b ，不同的叶尖速比所对应的C 值相差较大； P 对于给定的 b ，有且仅有一个固定的入 = 入 能使 C 达到最大值； opt P入 在风速不断变化的情况下,要保持 入 、o 必须随着风速按照 opt 的比例 opt tur R变化，才能保证风力机捕获的风能最大、效率最高。这是采用变速风电机组代替固定转速风电机组的初衷之一。图 2-7 风机C - 入 特性曲线P对于变桨距型风力发电机组， C 特性可近似表示为：PCP = 0.5(|( R入(C)f - 0.022 b - 2)|e-0.255 R 入(C)f(2- 14)式中， C 为叶片设计常数，一般取 13。f2.3.2 风力机的稳定工作区空间曲面虽然能包含风力机运行的所有状态点， 但是对于分析不太方便， 所以在实际应用中多是取几个离散的 b 值，画出平面图的方法，如下图所示，取6组 b 值，绘制如下：图 2-8 风力机稳定工作区曲线图在T(o, v)曲线中，以转矩T 的最大值为顶点连成的一条线 AB，将曲线簇分成了两部分，其中右侧为稳定运行区域，左侧部分为不稳定工作区域。 6对比P(o , v)和T(o, v)曲线， 我们发现当T 达到最大时， P 并没有达到最大， 具体而言，就是最大功率点对应的转速值要大于最大转矩点的转速值。如图 2-6 所示，曲线 CD 是由最大转矩点的连线而成的，曲线 EF 则是由最大功率点的连线而成。这 样一来，在 P(o , v)曲线簇中， CD 曲线和 EF 曲线之间的部分也是稳定区域。2.3.3 基于叶素理论的风力机建模基于叶素理论的风轮建模是将叶片分为若干个微元， 称为叶素， 通过对叶素的受力分析求得微元转矩， 再将所有微元转矩相加得到风力发电机风轮的输出转 矩 5,75 ，在风轮半径 r 处取一长度为 dr 的叶素，其弦长为 l ，节距角为 b 。( 1| dL = 1(2) pCl |dD = pC lw2 drl 2 d图 2-9 叶素微元受力分析图如上图所示，来流方向的风速为 v ，在半径为 r 处的风轮机速度为 u = O . r( O 为风轮机角速度)，气流相对于叶片的相对速度为 w ，则有：tw = v - ut(2- 15)叶素 dr 在相对速度为 w 的气流作用下，受到一个方向斜向上的气动力 dF 的作用。将 dF 沿与相对速度 w 垂直及水平方向可分解为升力 dL 和阻力 dD ，当 dr 很小时，可以近似的将叶素面积看成弦长与叶素长度的乘积，可得如下计算 公式：气动力 dF 按垂直和平行于旋转平面方向分解为法向力 dF 和切向力 dF，风a轮转矩 dT 由切向力 dF 产生，则有转矩微元：udT = r . dF = r(dL . sin I - dD . cos I)u令 c = C / C ，得总转矩计算公式：d lT = n . jdT = n . 12 prlw2 Cl . sin I(1 - c . cot I)drr0式中，r 轮毂半径；0n 风轮包含的桨叶个数；I 倾斜角(桨距角b 与攻角 i 之和)。u(2- 17)(2- 18)上式为基于叶素理论的风轮模型函数，可以写为如下形式： T = f(v, u, b )，即风轮输出转矩为风速 v ，风轮转速 u ，桨距角 b 的函数。2.3.4 基于 PSCAD 风力机模型与仿真基于 PSCAD 的风力机模型如下：图 2- 10 风力机简化模型基于上面的原理及理论公式，仿真结果如下：图 2- 11 风轮机机械转矩输出及参数设置图 2- 12 风轮机机械转矩模拟仿真结果图 2- 13 风轮机机械功率输出及参数设置图 2- 14 风轮机机械功率模拟仿真结果由上述两图曲线可知， 风轮机的输出转矩和输出功率都是标幺值， 则它们的 曲线是完全一致的，在 03s 时变桨距控制系统在调节桨叶节距使转矩和功率输 出逐渐达到稳定，由于又突然受到在 3s 与 4s 分别受到阵行风与渐变风的影响，从而使波形在这两个时间有的突变，之后继续达到稳定。3 变桨距风力发电机组控制系统模型3.1 变桨距风力发电机组的运行状态变桨距风力发电机组根据变距系统所起的作用可分为三种运行状态， 即风力 发电机组的起动状态(转速控制)、欠功率控制(不控制)和额定功率状态(功 率控制)。由于变桨距系统的响应速度受到限制 ,对快速变化的风速 ,通过改变节距来 控制输出功率的效果并不理想。 因此， 为了优化功率曲线， 最新设计的变桨风力 发电机组在进行功率控制的过程中， 其功率反馈信号不再作为直接控制桨叶节距 的变量。 变桨距系统由风速低频分量和发电机转速控制， 风速的高频分量产生的 机械能波动， 通过迅速改变发电机的转速来进行平衡， 即通过转子电流控制器对 发电机转差率进行控制， 当风速高于额定风速时， 允许发电机转速升高， 将瞬变的风能以风轮动能的形式储存起来； 转速降低时， 再将动能释放出来， 使功率曲线达到理想的状态。3.2 变桨距控制系统(1)变桨距控制系统工作原理图 3- 1 变桨距控制系统工作原理图在发电机并入电网时前， 发电机转速由速度控制器 A 根据发电机转速反馈信 号与给定信号直接控制；发电机并入电网后，速度控制 B 与功率控制器起作用。 功率控制器的任务主要是根据发电机转速给出相应的功率曲线， 调整发电机转差 率，并确定速度控制器 B 的速度给定。节距的给定参考值由控制器根据风力发电机组的运行状态给出。如图 3- 1 所示， 当风力发电机组并入电网前， 由速度控制器 A 给出； 当风力发电机组并入 电网后由速度控制 B 给出。(2)变距控制变距控制系统是一个随动系统， 如图 2-9 所示。 变距控制器是一个非线性比 例控制器，它可以补偿比例阀的死带和极限。变距系统的执行机构是液压系统， 节距控制器的输出信号经 D/A 转换后变成电压信号控制比例阀(或电液伺服阀)， 驱动液压缸活塞， 推动变桨距机构， 使桨叶节距角变化。 活塞的位移反馈信号由 位移传感器测量，经转换后输入比较器。图 3-2 变距控制结构框图(3)速度控制器 A转速控制器 A 在风力发电机组进入待机状态或从待机状态重新起动时投入 工作，如图 2- 10 所示在这些过程中通过对节距角的控制，转速以一定的变化率 上升。控制器也用于在同步速(50Hz 时 1500 转/min)时的控制。当发电机转速在同步转速10rmin 内持续 1s 发电机将切入电网。图 3-3 速度控制器 A控制器包含着常规的 PD 控制器和 PI 控制器，接着是节距角的非线性化环节， 通过非线性化处理， 增益随节距角的增加而减小， 以此补偿由于转子空气动力学 产生的非线性， 因为当功率不变时， 转矩对节距角的比是随节距角的增加而增加 的。当风力发电机组从待机状态进入运行状态时， 变桨距系统先将桨叶节距角快 速地转到 45，风轮在空转状态进入同步转速。当转速从 0 增加到 1500 rmin 时， 节距角给定值从 45线性的减小到 5。这一过程不仅使转子具有高起动力矩，而且在风速快速地增大时能够快速起动。(4)速度控制器 B发电机切入电网后，速度控制系统 B 作用。如图 2- 11 所示，速度控制器 B 受发电机转速和风速的双重控制。 在达到额定值前， 速度给定值随功率给定值按 比例增加。额定的速度给定值是 1569r/min,相应的发电机转差率是 4%。如果风 速和功率输出一直低于额定值，发电机转差率将降低到 2%，节距控制将根据风 速调整到最佳状态，以优化叶尖速比。如果风速高于额定值，发电机转速通过改变节距来跟踪相应的速度给定值。 功率输出将稳定地保持在额定值上。从图中可知在风速信号输入端设有低通滤波器，节距控制对瞬变风速并不响应。4 风力发电控制系统的模拟仿真分析4.1 无穷大系统模型的建立图 4- 1 风力发电机无穷大系统模型4.2 风力发电机系统并网模拟仿真分析(1) 发电机三相电压输出及仿真结果分析EbEcEd图 4-2 异步发电机发电机三相电压输出及参数设置图 4-3 异步发电机发电机三相电压模拟仿真结果图 4-4 高压母线电压模拟仿真结果异步发电机三相电压模拟仿真结果分析如图 4-3 所示,在发电机并入电网前低压侧电压为 0,在 1s 时发电机并入电 网, 电压突然上升到接近于额定电压,发电机向电网输送功率。并网后的功率输出 逐渐趋于稳定，这时的电压曲线呈正弦曲线形状变化，向电网输送额定功率。(2)低压母线和高压母线的线电压输出及仿真结果分析Main : GraphsEab(0.69)V K2.001.501.000.500.00-0.50-1.00-1.50-2.000.02.0 4.06.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0 .Main : GraphsEab(121)16014012010080V K60402000.02.0 4.06.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0 .图 4-5 低压母线和高压母线的线电压输出及仿真结果在正常运行时,低压母线和高压母线电压均从 0 迅速上升并均达到各自的额 定值，然后一直保持稳定。低压母线电压稳定在 0.69 kV 左右，高压母线电压稳 定在 121 kV 左右。(3) 低压母线相电流输出及仿真结果分析并网前电流为 0,在 1s 时断路器合闸并网，出现很大的冲击电流，其冲击电 流值达到 11KA， 最后开始衰减至 0.07KA， 然后又开始上升， 最后趋于稳定,其电 流最大稳定值为 0.64KA。IaIbIc图 4-6 低压电流输出及参数设置图 4-7 低压母线相电流输出及仿真结果4.3 变桨距控制系统模拟仿真分析通过控制桨距角的大小的改变就可以控制叶片吸收风功率的多少， 桨距角的调节可以使发电机输出功率平稳。图 4-8 变桨距控制系统模拟仿真结果变桨距控制系统模仿真结果分析如下：风轮机启动时风力发电机组开始自动运行于风轮叶尖本来值 90，即桨矩 角初始值为 90 度， 在机组起动的过程中逐渐变小， 这样叶片吸收风能逐渐增大， 叶片的转速也逐渐加快，最后在 1.4s 时桨矩角变为零，且保持不变，此时叶片吸收风能达到了最大。5 结论本文通过 PSCAD/EMTDC 电力系统模拟仿真软件， 建立了变桨距风力发电机组 控制系统模型， 对加入控制系统的风力发电样例系统进行模拟仿真分析， 验证了 控制系统模型的可用性。 风力发电系统控制策略是以风速的变化为依据， 风能的 最大利用效率为目的，为优化风力发电系统运行特性提出的控制方案。变桨距控制系统的设计主要采用 PI 控制器，根据发电机有功功率输出和风 轮机转速反馈来调节桨叶节距。通过风轮机桨距角控制系统对叶片桨距角进行控 制， 使风力发电机组的机械部分与发电机的电气部分配合， 达到提高风能利用效 率及改善供电质量的目的。 利用风力发电样例系统来验证控制系统的可用性， 并 对各种仿真曲线进行分析，从而得出结论。根据风速模型的仿真曲线， 分析风轮机和发电机各部分曲线的变化情况和整 个系统的仿真曲线图。 在并网以前电压的波形基本上是正弦形状的， 转速基本上 是稳定的。 并网以后虽然受到了电网的干扰， 但转速上升到额定转速后再没有多大变化； 电流的波形虽然是正弦的， 但整体的趋向也发生了相应的波动。 变桨距 控制系统在风力发电机组起动时，通过变距来获得足够的起动转矩；起动以后， 当低于额定风速运行时的机组状态控制为转速, 当高于额定风速运行时，通过调 整桨叶节距， 改变气流对叶片的攻角， 从而改变风力发电机组获得的空气动力转 矩，使功率输出保持稳定。额定风速之后的机组状态控制主要由桨距角调节实现。得到的控制系统保持了风力发电机组运行的安全可靠性。参考文献1 黄素逸. 能源与节能技术. 北京中国电力出版社. 2004 2 叶杭冶. 风力发电机组的控制技术. 机械工业出版社. 2002 3 宫靖远. 风电场工程技术手册. 机械工业出版社. 2004 4 王承煦 张源. 风力发电. 中国电力出版社. 20025 李俊峰. 风力 12 在中国. 北京化学工业出版社. 20056 赵凤山. 风力发电论文集. 金盾出版社. 20027 风电发展情况调研组， 2009 年中国风电发展情况调研报告 20091， 28 姚兴佳，宋俊等，风力发电机组原理与应用，北京：机械工业出版社， 2008 9 高峰，风力发电机组建模与变桨距控制研究： 博士学位论文，北京：华北电力大学， 200910 李东东， 风力发电机组并网与仿真分析， 水电能源科学， 2006 年，第 24 卷第 1 期.P： 1-311 何祚庥 王亦楠. 我国新能源可持续发展的现实选择. 中国三峡建设. 2005.P： 25-2712 肖运启双馈型风力发电机组励磁控制与优化运行研究： 博士学位论文 北京：华北电力大学， 200813 赵栋利，胡书举，赵斌，许洪华风力发电机、变流器及其低电压穿越概 述变频器世界， 2009， 214 张仰飞， 袁越， 陈小虎 风力机参数的可辨识分析 电力系统自动化， 2009， 33(6)： 86-8915 S.M.B.Wilmshurst. Control strategies for Wind turbines. WindEngineering. 1998 Vol. 12: P 236-24916 R. Chedid, F. Mrad and M. Basma. Intelligent Control for Wind EnergyConversion System. Wind Eng. 1998,Vol. 18:No. 1,1 16