发明名称 一种用于直驱永磁风电机变流器的多指标非线性控制方法
摘要 本发明公开了一种用于直驱永磁风电机变流器的多指标非线性控制方法,其包括以下步骤:1)确定风轮机、永磁同步发电机、电网侧、变流器直流环输入级,中间直流环节和输出级的电流、电磁转矩的数学模型,确定风电系统的约束方程;2)根据风电系统控制要求,建立非线性控制方法的微分代数系统;3)确定风电系统的非线性控制律。本发明针对直驱永磁风电系统脉冲宽度调制型变流器的非线性控制策略,控制占空比和调制角实现对脉冲宽度调制型变流器的控制,使输出电压的平滑性好和抗干扰性好,保持风电系统的小扰动稳定和改善风电系统的动态特性。即使电网电压不对称或电网电压跌落程度过大,也可以输出平滑的电压,保证电力和电气设备的正常运行。
申请公布号 CN103475009A 申请公布日期 2013.12.25
申请号 CN201310334453.3 申请日期 2013.08.02
申请人 李啸骢;王忠勇 发明人 李啸骢;王忠勇
分类号 H02J3/24(2006.01)I 主分类号 H02J3/24(2006.01)I
代理机构 福州君诚知识产权代理有限公司 35211 代理人 戴雨君
主权项 1.一种用于直驱永磁风电机变流器的多指标非线性控制方法,其特征在于:所述多指标非线性控制方法包括以下步骤: 1)确定风轮机、永磁同步发电机、电网侧、变流器直流环输入级,中间直流环节和输出级的电流、电磁转矩的数学模型,建立风电系统的约束方程; A1、风轮机产生的机械功率为风力机从空气中捕获的风能,风轮机的机械功率P<sub>m</sub>由式(I)所示:<img file="DEST_PATH_FDA0000392688110000011.GIF" wi="627" he="146" />其中λ=Rω/v<sub>w</sub>;     (I)式(I)中ρ为空气的密度,R为风轮机叶片半径,ν<sub>w</sub>为风速,C<sub>p</sub>为风能利用系数,是永磁体励磁磁链γ和桨距角β的非线性函数,β为叶片的桨距角,λ为风力机叶尖速比,ω为风力机的转速; B1、建立永磁同步发电机的数学模型: 分别以d、q轴建立坐标系,假设d、q坐标系以同步速度旋转,且q轴超前于d轴,将d轴定位于定子永磁体的磁链方向上,则在d、q轴坐标下的永磁同步发电机的数学模型为: <img file="DEST_PATH_FDA0000392688110000012.GIF" wi="867" he="172" /><img file="DEST_PATH_FDA0000392688110000013.GIF" wi="1518" he="191" />式(II)中:u<sub>1d</sub>、u<sub>1q</sub>分别为永磁同步发电机定子输出的d、q轴电压,i<sub>1d</sub>、i<sub>1q</sub>分别为永磁同步发电机定子输出的d、q轴电流;di<sub>1d</sub>、di<sub>1q</sub>分别表示i<sub>1d</sub>和i<sub>1q</sub>关于时间t的微分,L<sub>1d</sub>、L<sub>1q</sub>分别为永磁同步发电机定子直轴电感和交轴电感;γ为永磁体励磁磁链,不考虑温度影响时γ为一常数;n<sub>p</sub>为永磁同步发电机转子的极对数;ω<sub>m</sub>及R<sub>S</sub>分别为定子电角速度和相电阻; 根据电力电子技术理论可得到式(III)公式: u<sub>1d</sub>=u<sub>dc</sub>m<sub>1</sub>cosθ<sub>1</sub>,u<sub>1q</sub>=u<sub>dc</sub>m<sub>1</sub>sinθ<sub>1</sub>,u<sub>2d</sub>=u<sub>dc</sub>m<sub>2</sub>cosθ<sub>2</sub>,u<sub>2q</sub>=u<sub>dc</sub>m<sub>2</sub>sinθ<sub>2</sub>;    (III) 式(III)中,u<sub>dc</sub>为网侧直流侧电压;其中m<sub>1</sub>,θ<sub>1</sub>为变流器输入端正弦调制波的调制比和调制角;m<sub>2</sub>,θ<sub>2</sub>为变流器输出端正弦调制波的调制比和调制角。 当不考虑转子磁场的凸极效应并且电机气隙均匀时,L<sub>1d</sub>=L<sub>1q</sub>=L<sub>1</sub>,L<sub>1</sub>为定子电感;则式(II)转换为: <img file="DEST_PATH_FDA0000392688110000021.GIF" wi="966" he="162" /><img file="DEST_PATH_FDA0000392688110000022.GIF" wi="1574" he="188" />C1、建立电网侧的数学模型: 将电网电压综合矢量定向在d轴上,则d、q坐标下网侧模型为: <img file="DEST_PATH_FDA0000392688110000028.GIF" wi="1431" he="324" />式中:u<sub>2d</sub>、u<sub>2q</sub>分别为d、q轴电压,i<sub>2d</sub>、i<sub>2q</sub>分别为d、q轴电流;di<sub>2d</sub>、di<sub>2q</sub>分别表示i<sub>2d</sub>和i<sub>2q</sub>关于时间t的微分,e<sub>d</sub>、e<sub>q</sub>为电网d、q轴电压,ω<sub>f</sub>为电网角频率,L<sub>2</sub>、R分别为连接电感及等值电阻; 将式(III)代入式(V)后得: <img file="DEST_PATH_FDA0000392688110000029.GIF" wi="1522" he="328" />D1、变流器直流环输入和输出的电流为:设i<sub>dc1</sub>,i<sub>dc2</sub>分别为变流器直流环节输入和输出的电流,则<img file="DEST_PATH_FDA0000392688110000025.GIF" wi="1387" he="156" />式(VII)中,du<sub>dc</sub>表示u<sub>dc</sub>关于时间t的微分,i<sub>dc1</sub>,i<sub>dc2</sub>分别为变流器直流环节输入和输出的电流,C为变流器直流环节电容器的电容; 根据能量守恒,把式(III)代入式(VII)整理后得: <img file="DEST_PATH_FDA00003926881100000210.GIF" wi="1156" he="188" /><img file="DEST_PATH_FDA00003926881100000211.GIF" wi="1598" he="188" />E1、建立电磁转矩T<sub>e</sub>的数学方程:T<sub>e</sub>=1.5n<sub>p</sub>i<sub>1q</sub>[i<sub>1d</sub>(L<sub>1d</sub>-L<sub>1q</sub>)+γ];  (IX)由于L<sub>1d</sub>=L<sub>1q</sub>=L<sub>1</sub>,则T<sub>e</sub>=1.5n<sub>p</sub>i<sub>1q</sub>γ;   (X) 永磁同步发电机在d、q同步旋转坐标系下的d轴和<sub>q</sub>轴电压方程分别为 u<sub>1d</sub>-R<sub>s</sub>i<sub>1d</sub>+ω<sub>m</sub>L<sub>1</sub>i<sub>1q</sub>=0;u<sub>1q</sub>-R<sub>s</sub>i<sub>1q</sub>+ω<sub>m</sub>L<sub>1</sub>i<sub>1d</sub>+γ=0;  (XI)式中L<sub>1</sub>、R<sub>s</sub>及ω<sub>m</sub>分别为定子电感、相电阻、电角速度; F1、建立风电系统的约束方程: 结合式(I)至式(XI),可得风电系统的微分方程式(XII): <img file="DEST_PATH_FDA0000392688110000031.GIF" wi="957" he="166" /><img file="DEST_PATH_FDA0000392688110000032.GIF" wi="1112" he="173" /><img file="DEST_PATH_FDA0000392688110000033.GIF" wi="1153" he="179" /><img file="DEST_PATH_FDA0000392688110000034.GIF" wi="1132" he="187" /><img file="DEST_PATH_FDA0000392688110000035.GIF" wi="1659" he="163" /><img file="DEST_PATH_FDA0000392688110000036.GIF" wi="899" he="162" />G1、根据式(XII),得到风电系统的代数约束方程为式(XIII): T<sub>e</sub>-1.5n<sub>p</sub>i1<sub>q</sub>γ=0 u<sub>1d</sub>-R<sub>s</sub>i<sub>1d</sub>+ω<sub>m</sub>L<sub>1</sub>i<sub>1q</sub>=0;   (XIII) u<sub>1q</sub>-R<sub>s</sub>i<sub>1q</sub>+ω<sub>m</sub>L<sub>1</sub>i<sub>1d</sub>+γ=0 2)根据风电系统控制要求,建立非线性控制方法的微分代数系统: A2、微分代数系统的变换,多输入多输出非线性微分代数系统如式(XIV)所示: <img file="DEST_PATH_FDA0000392688110000041.GIF" wi="1543" he="551" />式(XIV)中x、y分别为多输出非线性微分代数系统的n和m维向量;f(x,y)∈R<sup>n</sup>,g<sub>i</sub>(x,y)∈R<sup>n</sup>,p(x,y)∈R<sup>m</sup>是对x和y都可微的非线性向量场;z=h(x,y)为光滑的向量场,u<sub>i</sub>为输入量;R<sup>n</sup>为n维变量空间,R<sup>m</sup>为m维变量空间; 定义1:M导数是一个标量函数,对某一标量函数h(x,y)做关于向量f(x,y)的M导数,即做如下运算: <img file="DEST_PATH_FDA0000392688110000042.GIF" wi="1270" he="227" />基于相同的原理,高阶M导数定义如下: <img file="DEST_PATH_FDA0000392688110000043.GIF" wi="1174" he="123" />定义2:当每个h<sub>i</sub>(x,y)对如式(XIV)所示的多输入多输出非线性微分代数系统的M导数<img file="DEST_PATH_FDA0000392688110000044.GIF" wi="493" he="85" />均为零,而<img file="DEST_PATH_FDA0000392688110000045.GIF" wi="220" he="85" />不为零,则多输入多输出非线性微分代数系统的相对阶为r<sub>i</sub>,多输入多输出非线性微分代数系统的总相对阶如式(XV)所示:<img file="DEST_PATH_FDA0000392688110000046.GIF" wi="1386" he="170" />当多输入多输出非线性微分代数系统总相对阶r小于系统维数n时,存在n-r个坐标映射; μ<sub>j</sub>(x,y)(j=1,2,…,n-r)满足M<sub>g</sub>μ<sub>j</sub>(x,y)=0(j=1,…,n-r); 如式(XIV)所示的多输入多输出非线性微分代数系统在坐标变换Φ(x,y)=[μ<sub>1</sub>,μ<sub>2</sub>,…μ<sub>r</sub>;η<sub>1</sub>,η<sub>2</sub>,…,η<sub>n-r</sub>]<sup>T</sup>时,表述如式(XVI)所示: <img file="DEST_PATH_FDA0000392688110000051.GIF" wi="1957" he="874" />B2、确定输出函数: 对输出量采用线性组合的形式,表达如式(XVII)所示: z<sub>i</sub>=h<sub>i</sub>(x,y)=αx+βy        (XVII) 式中未知参数α和β作为状态变量x和约束变量y的系数; C3、由线性控制律得到非线性控制律: 基于线性控制理论得到式(XVIII)所示的线性部分的反馈控制律: <img file="DEST_PATH_FDA0000392688110000052.GIF" wi="1486" he="168" />由式(XIX)所示的控制输入量u=B<sup>-1</sup>(x,y)(-A(x,y)+v)    (XIX) 式(XIX)中u=(u<sub>1</sub>,u<sub>2</sub>,...u<sub>m</sub>)<sup>T</sup>,v=(v<sub>1</sub>,v<sub>2</sub>,...v<sub>m</sub>)<sup>T</sup>,<img file="DEST_PATH_FDA0000392688110000053.GIF" wi="681" he="114" />得到多输入多输出非线性微分代数系统的非线性控制律如式(XX)所示: <img file="DEST_PATH_FDA0000392688110000054.GIF" wi="1635" he="327" />3)确定风电系统的非线性控制律: A3、式(XII)中风电系统的微分方程式的相关变量如式(XXI)及(XXII)所示: <img file="DEST_PATH_FDA0000392688110000061.GIF" wi="1299" he="117" /><img file="DEST_PATH_FDA0000392688110000062.GIF" wi="1555" he="361" /><img file="DEST_PATH_FDA0000392688110000063.GIF" wi="1676" he="629" />B3、确定风电系统的输出量:定子d轴电流i<sub>1d</sub>、输出电压d轴分量u<sub>2d</sub>、输出电压q轴分量u<sub>2q</sub>和网侧直流侧电压u<sub>dc</sub>; 将输出函数选为多个状态量的线性组合形式,因此风电系统的输出函数如式(XXIII)所示: h<sub>1</sub>(x,y)=β<sub>11</sub>T<sub>e</sub>+β<sub>12</sub>u<sub>1d</sub>h<sub>2</sub>(x,y)=α<sub>22</sub>i<sub>1q</sub>+α<sub>24</sub>i<sub>2d</sub>h<sub>3</sub>(x,y)=α<sub>34</sub>i<sub>2d</sub>+α<sub>33</sub>u<sub>dc</sub>        (XXIII) h<sub>4</sub>(x,y)=α<sub>44</sub>i<sub>2d</sub>+α<sub>45</sub>i<sub>2q</sub>+α<sub>43</sub>u<sub>dc</sub>C3、风电系统的非线性控制方法为: 根据输出函数h<sub>i</sub>(x,y),得到风电系统的总相对阶为<img file="DEST_PATH_FDA0000392688110000064.GIF" wi="444" he="147" />对以脉冲宽度调制型变流器为控制目标的非线性控制律,首先得到:<img file="DEST_PATH_FDA0000392688110000065.GIF" wi="388" he="249" /><img file="DEST_PATH_FDA0000392688110000066.GIF" wi="770" he="254" /><img file="DEST_PATH_FDA0000392688110000067.GIF" wi="941" he="378" /><img file="DEST_PATH_FDA0000392688110000071.GIF" wi="1241" he="146" /><img file="DEST_PATH_FDA0000392688110000072.GIF" wi="1513" he="343" />因r<sub>1</sub>=r<sub>2</sub>=r<sub>3</sub>=r<sub>4</sub>=1,则 <img file="DEST_PATH_FDA0000392688110000073.GIF" wi="2047" he="347" />将以上各式代入式(XX)后得控制输入如式(XXVI)所示: <img file="DEST_PATH_FDA0000392688110000074.GIF" wi="694" he="191" /><img file="DEST_PATH_FDA0000392688110000075.GIF" wi="716" he="192" /><img file="DEST_PATH_FDA0000392688110000076.GIF" wi="722" he="190" /><img file="DEST_PATH_FDA0000392688110000077.GIF" wi="1248" he="194" />其中 <img file="DEST_PATH_FDA0000392688110000078.GIF" wi="610" he="164" /><img file="DEST_PATH_FDA0000392688110000079.GIF" wi="568" he="152" /><img file="DEST_PATH_FDA00003926881100000710.GIF" wi="578" he="149" /><img file="DEST_PATH_FDA00003926881100000711.GIF" wi="601" he="161" />根据式(XXVI)即可对脉冲宽度调制型变流器中占空比m<sub>i</sub>和调制角θ<sub>i</sub>(i=1,2)进行控制,对风电系统的扰动进行调节。 
地址 530004 广西壮族自治区南宁市大学路100号(广西大学)