发明名称 一种用于永磁直驱风电系统速度辨识策略的选择方法
摘要 本发明公开了一种用于永磁直驱风电系统速度辨识策略的选择方法,该方法为:建立系统的小信号线性化模型,分析FIA、EPLL、VPLL、SKO控制算法下系统的小信号稳定性和稳态特性;通过时域仿真和硬件实验,从实际应用的层面分析并验证了理论分析的结果,揭示了FIA、EPLL、VPLL、SKO算法下系统的动态性能、算法复杂度、参数依赖性及对系统的影响;最终得出选择用于永磁直驱风电系统速度辨识策略的判别方法。本发明对选择采用无速度传感器辨识方法的风电系统建立了评判标准,提高了永磁直驱风电系统的高性能水平,较好的解决了现有的技术缺乏针对永磁同步发电机的转子转速与磁链位置的辨识方法的动稳态特性、参数依赖性等性能的综合对比判别方法的问题。
申请公布号 CN103701388B 申请公布日期 2017.01.11
申请号 CN201310544616.0 申请日期 2013.11.06
申请人 国家电网公司;国网天津市电力公司 发明人 李隆基;刘宝成;项添春;唐庆华;郗晓光;满玉岩;王浩鸣
分类号 H02P21/00(2016.01)I;H02P21/18(2016.01)I 主分类号 H02P21/00(2016.01)I
代理机构 天津盛理知识产权代理有限公司 12209 代理人 王来佳
主权项 一种用于永磁直驱风电系统速度辨识策略的选择方法,其特征在于,该用于永磁直驱风电系统速度辨识策略的选择方法包括以下步骤:步骤一,建立系统的小信号线性化模型,分析FIA、EPLL、VPLL、SKO控制算法下系统的小信号稳定性和稳态特性,以比较不同算法稳态特性的差异;步骤二,通过时域仿真和硬件实验,从实际应用的层面分析并验证了分析的结果,揭示了FIA、EPLL、VPLL、SKO算法下系统的动态性能、算法复杂度、参数依赖性及对系统的影响;步骤三,最终得出选择用于永磁直驱风电系统速度辨识策略的判别方法;该用于永磁直驱风电系统速度辨识策略的选择方法选定闭环系统运行状态为相同的转速、相同的负载,且参数均调整至使得辨识转速响应相同;在步骤一中,建立线性化小信号模型的方法为:发电机、转速、转矩双闭环控制器及辨识模块的数学模型,分别建立发电机、转速、转矩的状态方程如下式示:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mi>g</mi></msub><mo>=</mo><msub><mi>F</mi><mi>g</mi></msub><mo>(</mo><msub><mi>x</mi><mi>g</mi></msub><mo>,</mo><msub><mi>u</mi><mi>g</mi></msub><mo>)</mo></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mi>c</mi></msub><mo>=</mo><msub><mi>F</mi><mi>c</mi></msub><mo>(</mo><msub><mi>x</mi><mi>c</mi></msub><mo>,</mo><msub><mi>u</mi><mi>c</mi></msub><mo>)</mo></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mi>o</mi></msub><mo>=</mo><msub><mi>F</mi><mi>o</mi></msub><mo>(</mo><msub><mi>x</mi><mi>o</mi></msub><mo>,</mo><msub><mi>u</mi><mi>o</mi></msub><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000934783890000011.GIF" wi="286" he="198" /></maths>其中,x<sub>g</sub>=[i<sub>d</sub> i<sub>q</sub> ω<sub>r</sub> θ<sub>r</sub>]<sup>T</sup>为发电机模型的状态变量,x<sub>c</sub>=[x<sub>ASR</sub> x<sub>dACR</sub> x<sub>qACR</sub>]<sup>T</sup>为控制器模型的状态变量,x<sub>o</sub>为辨识模块的状态变量;u<sub>g</sub>、u<sub>c</sub>、u<sub>o</sub>为三者输入量,若将给定转速指令为常数,分析发电机侧控制系统对输入原动转矩T<sub>L</sub>扰动的特性,状态方程为:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>F</mi><mo>(</mo><mi>x</mi><mo>,</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><mi>G</mi><mo>(</mo><mi>x</mi><mo>,</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000934783890000012.GIF" wi="242" he="125" /></maths>其中x为系统的状态变量,<maths num="0003"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub></mtd><mtd><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>r</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000934783890000013.GIF" wi="200" he="62" /></maths>为系统的输出量,建立小信号线性模型如下式所述:<maths num="0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>F</mi></mrow><mrow><mo>&part;</mo><mi>x</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mi>x</mi><mi>N</mi></msub><mo>,</mo><msub><mi>T</mi><mi>L</mi></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub></mrow></msub><mi>x</mi><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>F</mi></mrow><mrow><mo>&part;</mo><msub><mi>T</mi><mi>L</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mi>x</mi><mi>N</mi></msub><mo>,</mo><msub><mi>T</mi><mi>L</mi></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub></mrow></msub><msub><mi>T</mi><mi>L</mi></msub></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>G</mi></mrow><mrow><mo>&part;</mo><mi>x</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mi>x</mi><mi>N</mi></msub><mo>,</mo><msub><mi>T</mi><mi>L</mi></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub></mrow></msub><mi>x</mi><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>G</mi></mrow><mrow><mo>&part;</mo><msub><mi>T</mi><mi>L</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mi>x</mi><mi>N</mi></msub><mo>,</mo><msub><mi>T</mi><mi>L</mi></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub></mrow></msub><msub><mi>T</mi><mi>L</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000934783890000021.GIF" wi="668" he="286" /></maths>分析输入为原动转矩T<sub>L</sub>、输出为辨识转速<img file="FDA0000934783890000022.GIF" wi="42" he="53" />和磁链位置<img file="FDA0000934783890000023.GIF" wi="34" he="61" />的系统幅频特性,三者均采用标幺值表示且基值分别为T<sub>LN</sub>、<img file="FDA0000934783890000024.GIF" wi="62" he="55" />及2π;在步骤一中,FIA算法辨识策略的建立,纯积分器中,PI调节器用于生成转子磁链幅值的补偿量,可写作:<maths num="0005"><math><![CDATA[<mrow><msub><mi>&psi;</mi><mrow><mi>c</mi><mi>o</mi><mi>m</mi><mi>p</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>k</mi><mi>p</mi></msub><mo>+</mo><mfrac><msub><mi>k</mi><mi>i</mi></msub><mi>s</mi></mfrac><mo>)</mo></mrow><mfrac><mrow><msub><mi>&psi;</mi><mi>&alpha;</mi></msub><msub><mi>u</mi><mi>&alpha;</mi></msub><mo>+</mo><msub><mi>&psi;</mi><mi>&beta;</mi></msub><msub><mi>u</mi><mi>&beta;</mi></msub></mrow><mrow><mo>|</mo><mi>&psi;</mi><mo>|</mo></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000934783890000025.GIF" wi="523" he="116" /></maths>在步骤一中,EPLL法辨识策略的建立,PI调节器与反正切函数可抑制环外扰动,u<sub>abc</sub>、i<sub>abc</sub>,磁链位置误差描述为:<maths num="0006"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>&theta;</mi><mo>=</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msub><mover><mi>e</mi><mo>^</mo></mover><mi>d</mi></msub><msub><mover><mi>e</mi><mo>^</mo></mover><mi>q</mi></msub></mfrac><mo>)</mo></mrow><mo>=</mo><mi>&alpha;</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msub><mover><mi>u</mi><mo>^</mo></mover><mi>d</mi></msub><mo>+</mo><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub><msub><mi>L</mi><mi>d</mi></msub><msub><mover><mi>i</mi><mo>^</mo></mover><mi>d</mi></msub></mrow><mrow><msub><mover><mi>u</mi><mo>^</mo></mover><mi>q</mi></msub><mo>-</mo><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub><msub><mi>L</mi><mi>q</mi></msub><msub><mover><mi>i</mi><mo>^</mo></mover><mi>q</mi></msub></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000934783890000026.GIF" wi="687" he="127" /></maths>EPLL只有一对PI参数,在EPLL的空载试验中,仅需测试辨识转子转速的阶跃响应;在步骤一中,VPLL法辨识策略的建立,VPLL直接将发动机端电压变换至d‑q坐标系,为得到转子磁链位置,需由端电压矢量位置求得的定子磁链位置再进行功率角的补偿,d轴电压指令参考值可以表述为:<maths num="0007"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>d</mi><mo>*</mo></msubsup><mo>=</mo><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub><msub><mi>L</mi><mi>q</mi></msub><msubsup><mi>i</mi><mi>q</mi><mo>*</mo></msubsup></mrow>]]></math><img file="FDA0000934783890000027.GIF" wi="190" he="62" /></maths>其中,<img file="FDA0000934783890000028.GIF" wi="44" he="51" />为定子电压矢量的旋转速度,<img file="FDA0000934783890000029.GIF" wi="35" he="62" />为定子磁链位置,δ为定子磁链位置与转子磁链位置的角度差,<img file="FDA00009347838900000210.GIF" wi="43" he="52" />为辨识转子转速,<img file="FDA00009347838900000211.GIF" wi="39" he="61" />为转子磁链位置;在步骤一中,SKO法辨识策略的建立,SKO法包括一组降维的状态变量与常量增益矩阵;在步骤二中,通过MATLAB/Simulink时域仿真与硬件实验的方法为:在时域仿真中,建立永磁直驱风电系统的Simulink模型,网侧采用电网电压定向的矢量控制,辨识转速与磁链位置的关系不一致,由于转速与磁链位置的积分关系,即使辨识转速有较大误差波动,但若转速误差曲线与误差零轴所包围的面积小,磁链位置误差便会小;硬件实验在15kW永磁直驱风电系统平台上进行,风力机由带变频器调速的异步电机模拟,FIA、EPLL、VPLL、SKO四种方法分别作系统辨识部分后,稳态运行时的电机转速为67.7rpm,负载转矩为97.8Nm,FIA法、EPLL法、VPLL法及SKO法四种辨识方法下发电机三相电流的频谱分析,每种方法控制下,最大的谐波含量为偶次谐波;辨识的磁链位置误差越大,控制下的电流谐波THD越大,磁链位置误差决定了系统的运行特性,在矢量控制中,磁链位置信息被用作坐标变换,获得更平稳的转矩应使得电流谐波THD尽量小,即需减小辨识磁链位置的误差;通过比较定子电感发生±25%变化时的运行特性,对比参数依赖性,仿真结果得出由弱至强依次为SKO、EPLL、VPLL、FIA;在步骤三中,FIA法由于算法复杂性限制实际应用,VPLL法的参数较难调试,控制效果差于FIA法、EPLL法,SKO法控制效果最差,EPLL法的有效性且易于实际数字控制器实现,适用于兆瓦级永磁直驱风电系统无速度传感器控制策略的转子转速、磁链位置辨识策略。
地址 100031 北京市西城区西长安街86号
您可能感兴趣的专利