第38卷第2期2021年6月
Vol.38 No.2Jun* 2021苏州科技大学学报(自然科学版)
Journal of Suzhou University of Science and Technology (Natural Science Edition)doi : 10.12084/j.issn.2096-3289.2021.02.010
及恒功率控制
乔焰辉1，郝诗源2，郝万君3* *，尚友涛4，王 昊1，孙志辉1
［收稿日期］2020-04-23
［基金项目］国家自然科学基金资助项目(51477109；61703296)
［作者简介］乔焰辉(1967-)，女，河北抚宁人，硕士，高级实验师，研究方向：机电一体化，新能源发电测控系统。 * 通信作者：郝万君(1965—)，男，博士，教授，硕士生导师,E-mail :****************. o
(1.苏州科技大学机械工程学院，江苏苏州215009；2.丹麦科技大学电气工程系，哥本哈根灵比2800；3.苏州科技大学电子 与信息工程学院，江苏苏州215009；4，吉化集团有限公司包装制品厂，吉林吉林132021)
摘要：风能的随机不确定性以及风轮机自身特性使风电机组的发电功率、风辂转速等随风速变化而波动，降低了 风电机组发电效率和输出电能质量。在分析变桨变速双馈风力发电机组的机电特性的基础上，提出了积分型滑模控 制(ISMC)和PI 复合控制策略。考虑到风能的不确定性、模型参数的时变性等对系统性能的影响，提出了一种积分滑 模控制(ISMC)策略。为了进一步降低系统进入滑模稳态区域的抖动、提升系统响应的快速性，设计了带自适应修正 项的滑模趋近控制率。仿真结果表明，所提出的控制方法能快速地跟随风速的变化，有效地抑制风轮转矩持续波动 所引起的传动链振荡，既提升了低风速段的发电效率，又能保证了高风速段风机输出功率的平稳性。 关键词：风力发电系统；最大功率跟踪；输出功率稳定；积分滑模控制；自适应趋近控制率
中图分类号：TM614 文献标志码：A 文章编号：2096-3289(2021 )02-0061-07
随着风力发电机组向着大型化、高参数化方向发展，对风力发电机组风能利用率、输出电能品质和运
行 平稳性的要求也越来越高叫然而，风力发电机组是具有高复杂度、强非线性、结构柔性化等特征的机电集成 系统，其空气动力特性会随着自然环境的变化而动态改变，具有较大的不确定性。建模误差、参数时变和未 建模动态等使得风电机组的精确控制非常困难3〕。
针对上述问题，文献［4］根据风电场风速统计概率与Bladed 软件辨识参数，建立了风力发电机组变增益 PI 控制器参数整定与优化方法。文献［5］提出一种在全风速范围内同时采用变桨和变速控制来调整发电机输 出有功功率的滑模控制策略，在降低风电机组输出有功功率的波动方面具有良好的效果，但牺牲了风能的利 用效率。文献［6］为提升双馈风力发电系统宜流侧电压平稳性、最大风能跟踪等性能，提出一种基于指数趋近 律的滑模变结构控制策略。文献［7］在综合分析风速限功率控制特性基础上，提出一种主动变速和桨距角控 制相结合的新型限功率控制策略，充分利用了机组转动惯量，在一定程度上提高发电量。文献［8］推导出了最 优桨距角随叶尖速比偏离程度变化的规律，并提出一种额定风速以下的转矩一变桨协调控制策略。
笔者在分析风电机组的机械、电气特性基础上，根据风电系统在高、低风速段的不同控制目标要求，结合 滑模控制在处理不确定性和抗干扰的优势，给出ISMC 和ISMC+PI 控制两种策略。并通过仿真实验验证了方 法的可行性和有效性。
1风力发电系统模型
变速变桨风力发电系统由空气动力系统、传动系统、发电机系统、变桨/转速控制器、变桨执行器及参数 测量环节等部分组成，如图1所示。在自然风的推动下空气动力系统会持续旋转，将风能转化为机械能；传
62苏州科技大学学报(自然科学版)2021年
动链系统将低转速转变为高转速；发电机系统的作用是将高转速的机械能转化为交流电能，然后经过电力电子电路变换后输送到供电电网；变桨/转速控制器是风力发电机组的核心控制单元，它依据系统的预期控制目标和控制策略对系统的运行状态及参数进行调节与控制；参数测量和变桨执行器分别为检测环节和执行机构。
MS 功率测
转
速
测
量
图1风力发电机组的总体结构图
1.1空气动力系统模型
以一种具有柔性水平传动轴的气动风力机为例，运用空气动力学、转子动力学和电机学等方面理论建立该系统的连续动态模型％
当有效风速为”(m/s)的自然风沿轴向吹过风轮机时，风轮机实际捕获到的有效风功率为
(1)
式中,P为空气密度,4是风轮扫过的区域面积,C p(j8,A)为风功率的利用系数，其表达式为
Q8,A)=(0.44-0.0167j8)sin|]-0.00184(入-3)0(2)
Q.GB.A)是以叶片桨距角0和叶尖速比A为自变量的函数，且有k=Rw r lv o其中©为风轮转子的角速度,
R为风轮半径。
风轮气动力矩7；可以直接从公式(1)推导出来
T f匕-=¥pARCp(J3，入)牛(3)
(t)r Z人
1.2传动链系统模型
对于非直驱式风力发电机组而言，从风轮转子到发电机的动力传递是通过传动链系统来实现的，它由低速轴、齿轮变速箱和高速轴连接而成，其模型可用以下一阶微分方程来描述阿
(D k+K r)(o r+^-(i)g(4)
芈尸普厂(伶~+Kg)3g-Tg(5)气动力矩7；驱动低速轴以角速度©旋转，经过齿轮变速箱将角速度提升到高速其传递与转换后的转矩7；来带动风力发电机旋转并发电。厶和厶分别为低速轴侧(含风轮转子)和高速轴侧(含齿轮箱、发电机)的转动惯量,0为高、低转速轴间的扭转角,H k和D a分别表示传动系统的转矩硬度和阻尼系数，心、念分
第2期乔焰辉，等：基于积分滑模控制的风电机组最大功率跟踪及恒功率控制63
别为低、高速轴的黏性摩擦系数，他为齿轮箱转数比。
由(4),(5)式，消除中间变量0可得
JQ f T^K^i)r-NgJ^i)g-N^/D g—N^r g(6)若忽略功率损耗，发电机输出功率为
戸尸丁処(7)
2风电系统的控制目标
风电系统根据自然风速的不同，大体可划分为切入区、低风速区、过渡区和高风速区，不同风速运行区的控制目标略有不同㈣。其中低风速运行区的控制目标是实现风能的最大功率跟踪(MPPT)，即在保持桨距角固定(0。附近)情况下，通过改变发电机转矩来调节风轮转速©，以取得最佳的叶尖速比入甲，从而实现风力机捕捉风能的最大化。在高风速运行区，主要控制目标是将发电机功率P g和发电机转速3尽可能地稳定在额定参数附近。此时系统有两个控制器，分别通过调节桨距角给定值0云和发电机转矩给定值7；何来加以实现。
3积分滑模控制器和PID控制器设计
3.1风轮转速的积分滑模控制器设计
滑模控制(Sliding Mode Control,SMC)又称变结构控制，其基本原理是通过设计滑模切换函数和趋近控制率，使系统按着给定的趋近率从初始位置运动到滑模切换面，然后保持系统的运行轨迹沿着切换面滑动至原点叫
考虑(1)-(7)式描述的非线性风电系统，其控制目标是当存在不确定性和随机干扰的情况下，使系统的某个输出量y能够快速地跟踪给定期望值沧，同时使跟踪误差变化率也是收敛于零。
3.1.1积分滑模切换函数设计
根据系统的控制目标，设计滑模切换函数
S(i)=ce(t)+e(i)(8)其中,c必须满足Hurwitz条件，即c>0，其数值的大小决定了误差的指数收敛速度。
为进一步改善控制系统的跟踪性能和消除稳态误差，将积分作用加入滑模切换函数设计，则有
S⑴=ce(/)+d(t)+民J。e⑴dt(9)
其中民是积分滑模增益。
3.1.2系统的李雅普诺夫稳定性
根据控制系统的稳定性准则，可选取以下形式的Lyapunov函数
V=(1/2)S2(10)保证系统稳定的充分条件是Lyapunov函数的导数小于等于零，代入(9)式可得
7=SS=S(ce(t)+e(t)+A：ie(t))
满足系统稳定性的条件是Uwo(且仅当sp=o时，卩=0)，其等价条件为
S<0,for S>0
S=0,for S=0
S>0,for S<0
为了保证SS<0，需要设计滑模控制律u，它由等效控制(喝)和切换控制(“J两个部分组成
W=U eq+W sw
见，用于控制整个系统运动行为，即跟踪控制;“讯用于克服系统不确定性和各种干扰。(11)
(12)
(13)
64苏州科技大学学报(自然科学版)2021年
3.1.3等效控制律(忸)设计
对于风轮转速控制，其跟踪误差可定义为
%(/)=<0厂如(14)将上式代入(11)式，可得
V=SS=S(ce(t)+e(t)+kfi(i))=S(c(&>r-&>rref)+e(t)+kfi(i))(15)为求取等效控制量，令上式中S=o可得
c(d>r-d>rf)+e(t)=0(16)将(6)代入上式，最终可以得到
c(7；-“
加&+哼佝+心”+认)_加十(/)+血e(t)=0(17)
Jr
由上式得出等效控制律为
%eq=^,ief=&~-J Q厂Kfi)厂^~3厂厶0,"丿一玄(広)+(f)(18)
Tig TVg Tl r C T ig C T ig
3.1.4切换控制律(uQ设计
为了满足到达条件，并保证正常段运动的动态性能，采用指数和幕次趋近律的控制律设计方法，可以表示为
S=-T]S-k ISI^sgnS(19)
式中，sgn表不符号函数，表不为
sgnS=1,S>0
-1,S<0
(20)
其中,-"S为指数趋近项。在指数趋近过程中，趋近速度从一个较大值逐步趋向于零，收敛速度取决于系数4。由于单纯的指数趋近，运动点逼近切换面是一个渐近的过程，不能保证在有限时间到达，为
此增加了無次趋近率项EISFsgnS。其作用是使得S趋近于零时，趋近速度是k(s)\S\a，而不是零。通过调整该系数值，可以保证系统状态远离滑动模态(S较大)时，能以较大的速度趋近于滑动模态；当系统状态趋近滑动模态(S较小)时，保证较小的控制增益，以降低发电机转矩的抖振。因此，该系数的取值不仅要与风轮转速密切相关，还应该能够充分反映发电机转矩的波动情况。
当风机以参考风速稳态运行时，发电机转矩7；接近7；胪，即
T gtOSt=—(T^-K/Or^)
n g
(21)
当风机处于暂态运行时，发电机转矩偏差为
(22)
为了抑制发电机转矩的抖振以及可能由此激发的传动系统扭转振动和塔架共振，在幕次趋近律的系数〃中引入与发电机转矩的偏差&相关的机制，给出了一种新的自适应趋近律的参数定义
虹£)=%皿-@唤(23)咯、咯分别为系数虹e)的最大值和最小值。其设计的目的是随着系统趋近滑模面程度和发电机转矩的波动幅值来自适应地调整增益系数&@)，以实现更短的到达时间和更小转矩颤振的效果。
由此，切换控制律(“Q可选取为
%=--(f]S+k(s')ISI a sgn5)(24)
叫
第2期乔焰辉，等：基于积分滑模控制的风电机组最大功率跟踪及恒功率控制65最终可求得ISMC 控制量为
瞪nf=Ueq+%= —~厂^—(1)——<w r ,…/+
s (f) +
e(t)_-(“S+&(£)ISIPgllS )
(25)
Tbg TLg H/g C T l/g C Z ig C T Lg 3.2 (高风速区)发电机输出功率控制器设计
在额定风速区以上时，保证发电机输出功率稳定是风电系统的控制目标之一。此时以桨距控制为主，转 矩控制为辅。对于风速的随机快速变化，只有两者很好地结合才能使整机的性能达到最优问。
由(1),(2)式可知，风机系统输出功率与桨距角之间具有复杂的非线性关系，所以桨距角控制不适宜采 用上述积分滑模控制策略，所以直接采用PI 控制器。而转矩控制回路仍然采用ISMC 控制。
风机发电功率的跟踪误差定义为
epCt^Pg-Pg,^ (26)0=蛀(/)+仏］e(t)d/
(27)为了取得良好的控制性能，采用粒子寻优方法获得控制器参数饨、仏。4仿真结果及分析
为了验证所提出的控制策略，文中采用Matlab/Simulink 仿真平台建立5 MW 风力发电机组的全阶非线 性模型，其主要参数：额定功率为5 MW ；齿轮箱变比为97；风轮额定角速度为1.267 1 rad ・s"；风轮转动惯量 为5.915 7xl07 kg-m 2；发电机额定角速度为122.909 6 rad-s'1；发电机转动惯量为534.116 kg-n?；风轮直径为 129 m ；额定风速为20 m ・s"。
忽略自然风的塔影效应和剪切效应，在高、低两个风速段进行系统仿真实验四。低风速段的平均风速
为 7 m ・s"、湍流强度设为20%，如图2(a)所示，其中包括实际风速和估计风速2条曲线(此段的风速估计采用 了文献［13］给出的扩展卡尔曼滤波的风速估计方法)。高风速段的平均风速为20 m ・s=湍流强度设为12%, 如图3(a)所示。
ISMC 的控制器参数为:c =25，勺=0.52,a =0.65，民=0.008 9,%唤=0.95,%血=0.1 ；PI 的控制器参数为：為=0.005 0，仏=4.56x10"。将所提控制策略与应用广泛的PID 控制进行性能对比分析。低风速段下ISMC 与 PID 控制效果如图2(b)、(c)、(d)所示及见表1。
时间
(a)随机风速曲线(平均风速为7 m-s-1)—ISMC
20 40 60 80 100 120 140 160 180 200
时间/s (b)风机输出功率曲线
(d)发电机转速曲线
图2低风速段ISMC 与PI
的控制效果
