基于MATLAB的风力发电系统仿真研究
摘要：本文介绍了风力发电机组的结构组成及原理，并建立了风力发电系统风速的数学模型、传动系统模型、发电机的数学模型, 并用MATLAB软件对风速模型进行了仿真, 结果证明了这些模型的正确性和有效性,说明了风力发电系统的仿真在对风力发电系统分析中的重要作用。
关键词：风力发电;MATLAB仿真; 动态模型; 风力发电机组
绪论
近几年来，风力发电机组单机容量和风电场建设规模都日益扩大，成为电网电源中的重要组成部分。风力的随机性和间歇性以及机组运行时的对无功的需求都会影响电力系统稳定运行。所以，在风电场建设前，需要论证分析风电场接入电网的可行性和确定允许接入的容量水平。作为分析的基础，需要建立正确的风电机组和风电场的数学模型。另外，针对新型风力发电机组，也需要根据其特性建立适当的数学模型，并应用于电力系统中，分析它的运行结果。因此，关于风力发电的课题研究是非常有必要的，对我国的能源结构调整将起到重要的推动作用。
1风力发电机结构组成原理
风力发电机组通常亦被称为风能转换系统。典型的并网型风力发电机组主要包括起支撑作用的塔架、风能的吸收和转换装置—风轮机（叶片、轮毂及其控制器）、起连接作用的传动机构—传动轴、齿轮箱、能量转换装置—发电机及其它风机运行控制系统—偏航系统和制动系统等。风力发电过程是：自然风吹转叶轮，带动轮毂转动，将风能转变为机械能，然后通过传动机构将机械能送至发电机转子，带动着转子旋转发电，实现由机械能向电能的转换，最后风电场将电能通过区域变电站注入电网。其能量转换过程是：风能→机械能→电能。
2 风力发电系统对并网运行的影响
2.1 风力发电机并网过程对电网的冲击影响
异步电机作为发电机运行时，没有独立的励磁装置，并网前发电机本身没有电压，因此并网时必然伴随一个过渡过程。异步发电机并网时的冲击电流的大小，与并网时网络电压的大小、发电机的暂态电抗以及并网时的滑差有关。滑差越大则交流暂态衰减时间就越长，并网时冲击电流有效值也就越大。目前可以通过加装软起动装置和风机非同期并网来削弱冲击电流，但同时给电网带来一定的谐波污染。
2.2 对电能质量的影响
风资源的不确定性和风电机组本身的运行特性使风电机组的输出功率是波动的，可能影响电网
的电能质量，包括对电网频率的影响和电网电压的影响，对电压的影响如电压偏差、电压波动和闪变、谐波以及周期性电压脉动等。
2.3 对电压稳定性的影响
对于电网稳定性的主要威胁一方面是风速的波动性和随机性引起风电场出力随时间变化且难以准确预测，而导致风力发电接入系统时潜在安全隐患。另一方面是弱电网中风电注入功率过高引起的电压稳定性降低。
2.4 对现有电网保护装置的影响
与常规配电网保护不同，通过风电场与电力系统联络线的潮流有时是双向的。风力发电机组在有风期间都是和电网相连的，当风速在起动风速附近变化时，为防止风电机组频繁投切对接触器的损害，允许风电机组短时电动机运行，此时会改变联络线的潮流方向，继电保护装置应充分考虑到这种运行方式。其次，并网运行的异步发电机没有独立的励磁机构，在电网发生短路故障时由于机端电压显著降低，异步发电机仅能提供短暂的冲击短路电流。
总之，风电场故障电流主要是公用电网电源提供的，风电场保护的技术困难是怎样根据有限的故障电流来识别故障的发生，使保护装置快速而准确的动作。
3 风力发电机组的数学模型 3.1 风速模型
风速是风力机的原动力，它的模型相对于风电机组比较独立。在电力系统稳态研究中，为了较精确地描述风的随机性和间歇性的特点，本研究中沿用国内外使用较多的风力四分量模型，各分量分别为基本风V A 、阵风V B 、渐变风V C 和随机风V D 。
由于风力机感受到的风速主要是轮毂高度H 处的风速V W ，风速从测风高度H 0到风力机轮毂高度H 必须进行修正。这在风速数据的处理和分析过程中是应该考虑的因素。修正公式为：
00w w H v v H α
⎡⎤
=⎢⎥⎣⎦
(3.1)
(其中α为高度修正系数，一般工程应用取1/7) 下面分别介绍各风力分量的计算公式。 （1）基本风
可以由风电场测风数据获得的威布尔分布参数近似确定，由威布尔分布的数学期望值可得：
11A V A K ⎛
⎫=∙Γ+ ⎪⎝⎭
(3.2)
V A 基本风速(m/s)，A 和K 是威布尔分布的尺度参数和形状参数，Γ(1+1/K)
表示伽马函数。 （2）阵风
描述风速突然变化的特性一般用阵风来表示。
()()()
111100G B S G G G G G t T V V T t T T t T T <⎧⎪
=≤≤+⎨⎪
>+⎩ (3.3)
其中1(max /2){1cos[2(/)(/)]}s G G G V G t T T T π=--
V B ，T 1G ，T G ，maxG 分别为阵风风速(m/s)，起动时间(s)，周期(s)和最大值(m/s)。 （3）渐变风
对风速的渐变特性可以用渐变风成分来表示。
()()()()
1112220max 0R R R C R R R R R t T V T t T V R T t T T T T T γ<⎧⎪
≤<⎪=⎨≤<+⎪⎪≥+⎩
(3.4) 其中212max [1(/)/()]R R R V R t T T T γ=--
式中V C ，maxR ，T 1R ，T 2R ，T R 分别为渐变风风速(m/s)、最大值(m/s )、起动时间(s)、终止时间(s)和保持时间(s)。 （4） 随机风
风速的随机性一般用随机噪声风分量来表示。
1
2
1
2[()]cos()N
D V i i i i V S ωωωϕ==∆-∑ （3.5）
其中：()()24
2321221/i N i
V i i i K F S F ωω
ωωπωμπ⎧⎛⎫=-∙∆ ⎪⎪⎝⎭⎪
⎨=⎪⎪⎡⎤
+⎣⎦⎩
式中：i ϕ指0～2π之间均匀分布的随机变量；K N 指地表粗糙系数；F 指扰动范围(2
m )；μ指相对高度的平均风速(m/s)；N 指频谱取样点数，i ω指各个频率段的频率。
综合上述四种风速成分，模拟实际作用在风力机上的风速为：
