发明名称 一种基于谐振滑模的DFIG系统控制方法
摘要 本发明公开了一种基于谐振滑模的DFIG系统控制方法,其采用直接转矩/功率控制而能够直接对电机电磁转矩和无功功率、网侧变流器输出的有功功率和无功功率进行控制,且消除了复杂的负序及各次谐波分量提取以及转子电流参考值计算步骤,故而不会引入分解延时。基于谐振的滑模控制可消除电机电磁转矩、无功功率及网侧变流器有功功率、无功功率中由于电网电压谐波所引入的负面影响,从而达到稳定的输出,较小的稳态误差以及较好的动态响应特性,进而达到预期控制效果。本发明无需进行复杂的负序及各次谐波分量提取,也无需进行电网电压相位的估计和复杂的电流参考值计算,能够大大减少控制延时,增强系统的快速性及稳定性。
申请公布号 CN103281025B 申请公布日期 2015.04.22
申请号 CN201310186591.1 申请日期 2013.05.17
申请人 浙江大学 发明人 年珩;全宇
分类号 H02P21/00(2006.01)I 主分类号 H02P21/00(2006.01)I
代理机构 杭州天勤知识产权代理有限公司 33224 代理人 周丽娟
主权项 一种基于谐振滑模的DFIG系统控制方法,包括如下步骤:对于DFIG机侧变流器控制:A1.采集DFIG的三相定子电压、三相定子电流、三相转子电流、转速以及转子位置角,根据转子位置角通过坐标变换确定DFIG的三相定子电压、三相定子电流、三相转子电流以及定子磁链在定子静止α‑β坐标系中的分量,进而计算出DFIG的输出有功功率P<sub>s</sub>、输出无功功率Q<sub>s</sub>以及电磁转矩T<sub>e</sub>;A2.使给定的电磁转矩参考量T<sub>eref</sub>和无功功率参考量Q<sub>sref</sub>分别减去DFIG的电磁转矩T<sub>e</sub>和输出无功功率Q<sub>s</sub>,得到电磁转矩误差量ΔT<sub>e</sub>和无功功率误差量ΔQ<sub>s</sub>;根据以下算式分别对电磁转矩误差量ΔT<sub>e</sub>和无功功率误差量ΔQ<sub>s</sub>进行PIR调节,得到电磁转矩滑模值S<sub>Te</sub>和无功功率滑模值S<sub>Qs</sub>;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>S</mi><mi>Te</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>K</mi><mrow><mi>sp</mi><mn>1</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>K</mi><mrow><mi>si</mi><mn>1</mn></mrow></msub><mi>s</mi></mfrac><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mrow><mi>sr</mi><mn>1</mn></mrow></msub><mi>s</mi></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mn>6</mn><mi>&omega;</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>T</mi><mi>e</mi></msub></mrow>]]></math><img file="FDA0000665048760000011.GIF" wi="887" he="195" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>S</mi><mi>Qs</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>K</mi><mrow><mi>sp</mi><mn>2</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>K</mi><mrow><mi>si</mi><mn>2</mn></mrow></msub><mi>s</mi></mfrac><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mrow><mi>sr</mi><mn>2</mn></mrow></msub><mi>s</mi></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mn>69</mn><mi>&omega;</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>Q</mi><mi>s</mi></msub></mrow>]]></math><img file="FDA0000665048760000012.GIF" wi="918" he="188" /></maths>其中:s为拉普拉斯算子,K<sub>sp1</sub>和K<sub>sp2</sub>均为给定的比例系数,K<sub>si1</sub>和K<sub>si2</sub>均为给定的积分系数,K<sub>sr1</sub>和K<sub>sr2</sub>均为给定的谐振系数,ω=2πf,f=50Hz;A3.根据所述的电磁转矩滑模值S<sub>Te</sub>和无功功率滑模值S<sub>Qs</sub>通过以下滑模控制算法,计算得到转子电压指令在定子静止α‑β坐标系中的分量;U<sub>rα</sub>=(U<sub>sα</sub>·C‑ψ<sub>sβ</sub>·D)/(U<sub>sα</sub>ψ<sub>sα</sub>‑U<sub>sβ</sub>ψ<sub>sβ</sub>)U<sub>rβ</sub>=(U<sub>sβ</sub>·C‑ψ<sub>sα</sub>·D)/(U<sub>sα</sub>ψ<sub>sα</sub>‑U<sub>sβ</sub>ψ<sub>sβ</sub>)C=σL<sub>m</sub>[K<sub>si1</sub>ΔT<sub>e</sub>+K<sub>Te</sub>sgn(S<sub>Te</sub>)‑pQ<sub>s</sub>‑1.5(A·ψ<sub>sα</sub>‑B·ψ<sub>sβ</sub>)]/1.5pD=‑σL<sub>m</sub>[K<sub>si2</sub>ΔQ<sub>s</sub>+K<sub>Qs</sub>sgn(S<sub>Qs</sub>)‑ωP<sub>s</sub>+1.5(A·U<sub>sβ</sub>‑B·U<sub>sα</sub>)]/1.5pA=[‑L<sub>r</sub>(U<sub>sα</sub>+ω<sub>r</sub>ψ<sub>sβ</sub>)/L<sub>m</sub>‑ω<sub>r</sub>σL<sub>m</sub>I<sub>sβ</sub>]/σL<sub>m</sub>B=[‑L<sub>r</sub>(U<sub>sβ</sub>‑ω<sub>r</sub>ψ<sub>sα</sub>)/L<sub>m</sub>+ω<sub>r</sub>σL<sub>m</sub>I<sub>sα</sub>]/σL<sub>m</sub><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>sgn</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>Te</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Te</mi></msub><mo>></mo><mn>0.05</mn></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>Te</mi></msub><mo>/</mo><mn>0.05</mn><mo>,</mo></mtd><mtd><mo>-</mo><mn>0.05</mn><mo>&le;</mo><msub><mi>S</mi><mi>Te</mi></msub><mo>&le;</mo><mn>0.05</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Te</mi></msub><mo>&lt;</mo><mo>-</mo><mn>0.05</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000665048760000021.GIF" wi="1170" he="300" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>sgn</mi><mrow><mo>(</mo><msub><mi>s</mi><mi>Qs</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Qs</mi></msub><mo>></mo><mn>0.05</mn></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>Qs</mi></msub><mo>/</mo><mn>0.05</mn><mo>,</mo></mtd><mtd><mo>-</mo><mn>0.05</mn><mo>&le;</mo><msub><mi>S</mi><mi>Qs</mi></msub><mo>&le;</mo><mn>0.05</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Qs</mi></msub><mo>&lt;</mo><mo>-</mo><mn>0.05</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000665048760000022.GIF" wi="1187" he="312" /></maths>其中:U<sub>rα</sub>和U<sub>rβ</sub>分别为转子电压指令在定子静止α‑β坐标系中的α轴分量和β轴分量,p为DFIG的极对数,Ψ<sub>sα</sub>和Ψ<sub>sβ</sub>分别为定子磁链在定子静止α‑β坐标系中的α轴分量和β轴分量,I<sub>sα</sub>和I<sub>sβ</sub>分别为三相定子电流在定子静止α‑β坐标系中的α轴分量和β轴分量,U<sub>sα</sub>和U<sub>sβ</sub>分别为三相定子电压在定子静止α‑β坐标系中的α轴分量和β轴分量,σ为DFIG的漏感系数且σ=1‑(L<sub>s</sub>L<sub>r</sub>/L<sub>m</sub><sup>2</sup>),L<sub>s</sub>和L<sub>r</sub>分别为DFIG的定子电感和转子电感,L<sub>m</sub>为DFIG的定转子互感,ω<sub>r</sub>为DFIG的转速,K<sub>Τe</sub>和K<sub>Qs</sub>均为给定的收敛控制系数;A4.对转子电压指令在定子静止α‑β坐标系中的分量进行Park变换,得到转子电压指令在转子静止α‑β坐标系中的分量;进而根据转子电压指令在转子静止α‑β坐标系中的分量通过SVPWM技术构造得到一组PWM信号以对DFIG机侧变流器进行控制;对于DFIG网侧变流器控制:B1.采集DFIG网侧变流器的三相进线电流和直流母线电压V<sub>dc</sub>,通过坐标变换确定三相进线电流在定子静止α‑β坐标系中的分量,进而计算出DFIG网侧变流器的输出有功功率P<sub>g</sub>和输出无功功率Q<sub>g</sub>;B2.使给定的直流母线电压参考量V<sub>dcref</sub>减去直流母线电压V<sub>dc</sub>,得到直流母线电压误差量ΔV<sub>dc</sub>;对直流母线电压误差量ΔV<sub>dc</sub>进行PI调节,得到直流侧有功功率P<sub>d</sub>;B3.提取定子电压基波、定子电流基波和定子电流谐波分别在定子静止α‑β坐标系中的分量;进而根据实际控制目标和直流侧有功功率P<sub>d</sub>,计算出在该控制目标下DFIG网侧变流器的输出有功功率参考值P<sub>gref</sub>和输出无功功率参考值Q<sub>gref</sub>;B4.使输出有功功率参考值P<sub>gref</sub>和输出无功功率参考值Q<sub>gref</sub>分别减去输出有功功率P<sub>g</sub>和输出无功功率Q<sub>g</sub>,得到有功功率误差量ΔP<sub>g</sub>和无功功率误差量ΔQ<sub>g</sub>;根据以下算式分别对有功功率误差量ΔP<sub>g</sub>和无功功率误差量ΔQ<sub>g</sub>进行PIR调节,得到有功功率滑模值S<sub>Pg</sub>和无功功率滑模值S<sub>Qg</sub>;<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>S</mi><mi>Pg</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>K</mi><mrow><mi>gp</mi><mn>1</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>K</mi><mrow><mi>gi</mi><mn>1</mn></mrow></msub><mi>s</mi></mfrac><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mrow><mi>gr</mi><mn>1</mn></mrow></msub><mi>s</mi></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mn>6</mn><mi>&omega;</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>P</mi><mi>g</mi></msub></mrow>]]></math><img file="FDA0000665048760000031.GIF" wi="939" he="204" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>S</mi><mi>Qg</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>K</mi><mrow><mi>gp</mi><mn>2</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>K</mi><mrow><mi>gi</mi><mn>2</mn></mrow></msub><mi>s</mi></mfrac><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mrow><mi>gr</mi><mn>2</mn></mrow></msub><mi>s</mi></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mn>6</mn><mi>&omega;</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>Q</mi><mi>g</mi></msub></mrow>]]></math><img file="FDA0000665048760000032.GIF" wi="968" he="204" /></maths>其中:K<sub>gp1</sub>和K<sub>gp2</sub>均为给定的比例系数,K<sub>gi1</sub>和K<sub>gi2</sub>均为给定的积分系数,K<sub>gr1</sub>和K<sub>gr2</sub>均为给定的谐振系数;B5.根据所述的有功功率滑模值S<sub>Pg</sub>和无功功率滑模值S<sub>Qg</sub>通过以下滑模控制算法,计算得到网侧电压指令在定子静止α‑β坐标系中的分量;U<sub>gα</sub>=ρ(EU<sub>sα</sub>‑FU<sub>sβ</sub>) U<sub>gβ</sub>=ρ(EU<sub>sβ</sub>+FU<sub>sβ</sub>)<img file="FDA0000665048760000035.GIF" wi="481" he="115" /><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><msub><mrow><mn>1.5</mn><mi>L</mi></mrow><mi>g</mi></msub><mo>[</mo><msub><mi>K</mi><mrow><mi>gi</mi><mn>1</mn></mrow></msub><mi>&Delta;</mi><msub><mi>P</mi><mi>g</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Pg</mi></msub><mi>sgn</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>Pg</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mi>g</mi></msub><msub><mi>P</mi><mi>g</mi></msub><mo>/</mo><msub><mi>L</mi><mi>g</mi></msub><mo>+</mo><mi>&omega;</mi><msub><mi>Q</mi><mi>g</mi></msub><mo>]</mo><mo>+</mo><mfrac><mn>1</mn><mi>&rho;</mi></mfrac></mrow>]]></math><img file="FDA0000665048760000034.GIF" wi="1519" he="179" /></maths>F=‑1.5L<sub>g</sub>[K<sub>gi2</sub>ΔQ<sub>g</sub>+K<sub>Qg</sub>sgn(S<sub>Qg</sub>)+R<sub>g</sub>Q<sub>g</sub>/L<sub>g</sub>‑ωP<sub>g</sub>]<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>sgn</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>Pg</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Pg</mi></msub><mo>></mo><mn>0.05</mn></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>Pg</mi></msub><mo>/</mo><mn>0.05</mn><mo>,</mo></mtd><mtd><mo>-</mo><mn>0.05</mn><mo>&le;</mo><msub><mi>S</mi><mi>Pg</mi></msub><mo>&le;</mo><mn>0.05</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Pg</mi></msub><mo>&lt;</mo><mo>-</mo><mn>0.05</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000665048760000033.GIF" wi="1206" he="309" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>sgn</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>Qg</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Qg</mi></msub><mo>></mo><mn>0.05</mn></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>Qg</mi></msub><mo>/</mo><mn>0.05</mn><mo>,</mo></mtd><mtd><mo>-</mo><mn>0.05</mn><mo>&le;</mo><msub><mi>S</mi><mi>Qg</mi></msub><mo>&le;</mo><mn>0.05</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>S</mi><mi>Qg</mi></msub><mo>&lt;</mo><mo>-</mo><mn>0.05</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000665048760000041.GIF" wi="1206" he="309" /></maths>其中:U<sub>gα</sub>和U<sub>gβ</sub>分别为网侧电压指令在定子静止α‑β坐标系中的α轴分量和β轴分量,K<sub>Pg</sub>和K<sub>Qg</sub>均为给定的收敛控制系数,L<sub>g</sub>为DFIG网侧变流器的进线电感,R<sub>g</sub>为DFIG网侧变流器进线电感的内阻;进而根据网侧电压指令在定子静止α‑β坐标系中的分量通过SVPWM技术构造得到一组PWM信号以对DFIG网侧变流器进行控制。
地址 310027 浙江省杭州市西湖区浙大路38号