发明名称 一种双馈式风力发电机转子暂态电流的计算方法
摘要 本发明提供一种双馈式风力发电机转子暂态电流的计算方法,包括:在同步速旋转坐标系中,建立电网电压骤升时双馈异步发电机的暂态数学模型;根据暂态数学模型,计算电网电压骤升时的定、转子暂态电流方程;对于设定的电压骤升投入机制,计算定子磁链和转子开路时的定子磁链方程以及从定子侧感应的转子磁链方程;根据电网电压骤升时的定、转子暂态电流方程,以及,定子磁链和转子开路时的定子磁链方程、从定子侧感应的转子磁链方程,得到电网电压对称骤升后,由于定子侧电气量的骤变而感生出的转子暂态电流表达式。本发明计算过程简单、结果准确,为采用双馈异步发电机的电力系统保护整定工作提供可靠的数据基础。
申请公布号 CN105703677A 申请公布日期 2016.06.22
申请号 CN201410687303.5 申请日期 2014.11.25
申请人 国家电网公司;国网新源张家口风光储示范电站有限公司;华北电力科学研究院有限责任公司;华北电力大学(保定) 发明人 刘少宇;常静;白恺;刘汉民;宋鹏;朱斯;刘京波
分类号 H02P9/10(2006.01)I;H02P21/14(2006.01)I 主分类号 H02P9/10(2006.01)I
代理机构 北京三友知识产权代理有限公司 11127 代理人 汤在彦
主权项 一种双馈式风力发电机转子暂态电流的计算方法,其特征在于,包括:步骤1,在同步速旋转坐标系中,建立电网电压骤升时双馈异步发电机的暂态数学模型,包括如下定、转子电压方程和定、转子磁链方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>u</mi><mi>s</mi></msub><mo>=</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>s</mi></msub><mo>+</mo><mfrac><mrow><mi>d</mi><msub><mi>&psi;</mi><mi>s</mi></msub></mrow><mi>dt</mi></mfrac></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>r</mi></msub><mo>=</mo><msub><mi>R</mi><mi>r</mi></msub><msub><mi>i</mi><mi>r</mi></msub><mo>+</mo><mfrac><mrow><mi>d</mi><msub><mi>&psi;</mi><mi>r</mi></msub></mrow><mi>dt</mi></mfrac><mo>-</mo><mi>j</mi><msub><mi>&omega;</mi><mi>r</mi></msub><msub><mi>&psi;</mi><mi>r</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000615907720000011.GIF" wi="555" he="309" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&psi;</mi><mi>s</mi></msub><mo>=</mo><msub><mi>L</mi><mi>s</mi></msub><msub><mi>i</mi><mi>s</mi></msub><mo>+</mo><msub><mi>L</mi><mi>m</mi></msub><msub><mi>i</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&psi;</mi><mi>r</mi></msub><mo>=</mo><msub><mi>L</mi><mi>m</mi></msub><msub><mi>i</mi><mi>s</mi></msub><mo>+</mo><msub><mi>L</mi><mi>r</mi></msub><msub><mi>i</mi><mi>r</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000615907720000012.GIF" wi="409" he="191" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>L</mi><mi>s</mi></msub><mo>=</mo><msub><mi>L</mi><mi>s&sigma;</mi></msub><mo>+</mo><msub><mi>L</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>r</mi></msub><mo>=</mo><msub><mi>L</mi><mi>r&sigma;</mi></msub><mo>+</mo><msub><mi>L</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000615907720000013.GIF" wi="361" he="182" /></maths>其中,u<sub>s</sub>、u<sub>r</sub>分别为定、转子电压矢量;R<sub>s</sub>、R<sub>r</sub>分别为定、转子电阻;i<sub>s</sub>、i<sub>r</sub>分别为定、转子电流矢量;ψ<sub>s</sub>、ψ<sub>r</sub>分别为定、转子磁链矢量;ω<sub>r</sub>为转子角速度;L<sub>sσ</sub>、L<sub>rσ</sub>分别为定、转子漏感,L<sub>m</sub>为定、转子互感,t为时间;步骤2,根据所述暂态数学模型,计算电网电压骤升时的定、转子暂态电流方程:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>i</mi><mi>s</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>L</mi><mi>s</mi></msub><mo>-</mo><mfrac><msub><mi>L</mi><mi>r</mi></msub><msub><mi>L</mi><mi>m</mi></msub></mfrac></mrow></mfrac><msub><mi>&psi;</mi><mi>s</mi></msub><mo>-</mo><mfrac><mn>1</mn><mrow><mfrac><mrow><msub><mi>L</mi><mi>s</mi></msub><msub><mi>L</mi><mi>r</mi></msub></mrow><msub><mi>L</mi><mi>m</mi></msub></mfrac><mo>-</mo><msub><mi>L</mi><mi>m</mi></msub></mrow></mfrac><msub><mi>&psi;</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>r</mi></msub><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mrow><mfrac><mrow><msub><mi>L</mi><mi>s</mi></msub><msub><mi>L</mi><mi>r</mi></msub></mrow><msub><mi>L</mi><mi>m</mi></msub></mfrac><mo>-</mo><msub><mi>L</mi><mi>m</mi></msub></mrow></mfrac><msub><mi>&psi;</mi><mi>s</mi></msub><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>L</mi><mi>r</mi></msub><mo>-</mo><mfrac><msup><msub><mi>L</mi><mi>m</mi></msub><mn>2</mn></msup><msub><mi>L</mi><mi>s</mi></msub></mfrac></mrow></mfrac><msub><mi>&psi;</mi><mi>r</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000615907720000014.GIF" wi="779" he="479" /></maths>步骤3,对于设定的电压骤升投入机制,计算定子磁链和转子开路时的定子磁链方程以及从定子侧感应的转子磁链方程:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&psi;</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>></mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mfrac><msub><mi>V</mi><mi>se</mi></msub><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><mi>t</mi></mrow></msup><mo>-</mo><mi>p</mi><mfrac><msub><mi>V</mi><mi>se</mi></msub><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><msub><mi>t</mi><mn>0</mn></msub></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>1</mn></msub></mrow></msup></mtd></mtr><mtr><mtd><msub><mi>&psi;</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>></mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>L</mi><mi>m</mi></msub><msub><mi>L</mi><mi>s</mi></msub></mfrac><mo>[</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mfrac><msub><mi>V</mi><mi>se</mi></msub><mrow><mi>j</mi><msub><mi>&psi;</mi><mi>s</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><mi>t</mi></mrow></msup><mo>-</mo><mi>p</mi><mfrac><msub><mi>V</mi><mi>se</mi></msub><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><msub><mi>t</mi><mn>0</mn></msub></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>]</mo></mtd></mtr><mtr><mtd><msub><mi>&tau;</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>L</mi><mi>s</mi></msub><msub><mi>R</mi><mi>s</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>&tau;</mi><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>L</mi><mi>r</mi></msub><msub><mi>R</mi><mi>r</mi></msub></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000615907720000021.GIF" wi="1228" he="605" /></maths>其中,所述设定的电压骤升投入机制为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>u</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>&lt;</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>V</mi><mi>se</mi></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><mi>t</mi></mrow></msup></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>></mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>V</mi><mi>se</mi></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><mi>t</mi></mrow></msup><mo>+</mo><mi>p</mi><msub><mi>V</mi><mi>se</mi></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><mi>t</mi></mrow></msup></mtd></mtr><mtr><mtd><mi>p</mi><mo>=</mo><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>u</mi><mi>s</mi></msub><mo>|</mo><mo>|</mo><mo>-</mo><msub><mi>V</mi><mi>se</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>V</mi><mi>se</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000615907720000022.GIF" wi="868" he="346" /></maths>其中,ψ<sub>r0</sub>为从定子侧感应的转子磁链;ω<sub>s</sub>为定子角速度;V<sub>se</sub>为正常运行时定子电压幅值;p为电网电压骤升度;τ<sub>1</sub>、τ<sub>2</sub>为衰减时间常数;t<sub>0</sub>为电网电压对称骤升时刻;步骤4,根据所述电网电压骤升时的定、转子暂态电流方程,以及,所述定子磁链和转子开路时的定子磁链方程、从定子侧感应的转子磁链方程方程,得到电网电压对称骤升后,由于定子侧电气量的骤变而感生出的转子暂态电流表达式:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>i</mi><mi>r</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mfrac><mrow><msub><mi>L</mi><mi>s</mi></msub><msub><mi>L</mi><mi>r</mi></msub></mrow><msub><mi>L</mi><mi>m</mi></msub></mfrac><mo>-</mo><msub><mi>L</mi><mi>m</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><mi>p</mi><mfrac><msub><mi>V</mi><mi>se</mi></msub><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&omega;</mi><mi>s</mi></msub><msub><mi>t</mi><mn>0</mn></msub></mrow></msup><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>1</mn></msub></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000615907720000023.GIF" wi="1580" he="285" /></maths>
地址 100031 北京市西城区西长安街86号