发明名称 一种双馈风电机组超速脱网临界时间的计算方法
摘要 本发明是一种双馈风电机组超速脱网临界时间的计算方法,其特点是:本文从建立DFIG受扰后机电暂态数学模型入手,分析了DFIG初始运行条件、扰动强度以及Crowbar保护动作等因素在机电暂态中的作用,推导了DFIG机组达到超速保护动作的条件,并且用所提出的方法通过对某风电场风电机群脱网实例的分析,根据所计算的各机组转子加速至超速保护动作定值的时间,正确识别了超速保护动作的风电机组。该方法计算简单,具有计算快速、评价精度高、应用价值高、在线应用前景大等优点。
申请公布号 CN104078999B 申请公布日期 2016.03.30
申请号 CN201410326344.1 申请日期 2014.07.09
申请人 国网吉林省电力有限公司;东北电力大学;国家电网公司 发明人 郑太一;穆钢;孙勇;严干贵;蔡宏毅;杨国新;王健;张继国;沈清坤
分类号 H02J3/38(2006.01)I;G06F19/00(2011.01)I 主分类号 H02J3/38(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种双馈风电机组超速脱网临界时间的计算方法,其特征在于,它包括以下步骤:1).建立DFIG机电暂态数学模型DFIG转子运动方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mrow><msub><mi>d&omega;</mi><mi>r</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><msub><mi>&omega;</mi><mi>N</mi></msub><msub><mi>T</mi><mi>j</mi></msub></mfrac><mrow><mo>(</mo><msubsup><mi>T</mi><mi>m</mi><mo>*</mo></msubsup><mo>-</mo><msubsup><mi>T</mi><mi>e</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000836293350000011.GIF" wi="717" he="134" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><msup><mi>T</mi><mo>*</mo></msup><mi>m</mi></msub><mo>=</mo><msub><msup><mi>P</mi><mo>*</mo></msup><mi>m</mi></msub><mo>/</mo><msub><msup><mi>&Omega;</mi><mo>*</mo></msup><mi>r</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><msup><mi>T</mi><mo>*</mo></msup><mi>e</mi></msub><mo>=</mo><msub><msup><mi>P</mi><mo>*</mo></msup><mrow><mi>e</mi><mi>s</mi></mrow></msub><mo>-</mo><msub><msup><mi>T</mi><mo>*</mo></msup><mrow><mi>e</mi><mi>r</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000836293350000012.GIF" wi="693" he="149" /></maths>式(1)和(2)中:上标*表示标幺值;ω<sub>r</sub>为转子电气角速度;Ω<sub>r</sub>为转子机械角速度;ω<sub>1</sub>为同步角速度;ω<sub>N</sub>为电网角频率;P<sub>m</sub>为机械功率;T<sub>j</sub>为惯性时间常数,T<sub>m</sub>为机械转矩;T<sub>e</sub>为电磁转矩,由同步转矩分量T<sub>es</sub>和异步转矩分量T<sub>er</sub>构成,若Crowbar保护动作,转子绕组经电阻R<sub>cb</sub>短路,DFIG进入异步电机运行状态,此时电磁转矩平均值T<sub>eav</sub>为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>e</mi><mi>a</mi><mi>v</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>3</mn><msubsup><mi>pU</mi><mn>1</mn><mn>2</mn></msubsup><mo>&times;</mo><mrow><mo>(</mo><msub><mi>R</mi><mi>r</mi></msub><mo>+</mo><msub><mi>R</mi><mrow><mi>c</mi><mi>b</mi></mrow></msub><mo>)</mo></mrow><mo>/</mo><mi>s</mi></mrow><mrow><mn>2</mn><msub><mi>&pi;f</mi><mi>N</mi></msub><mo>&lsqb;</mo><msup><mrow><mo>(</mo><msub><mi>R</mi><mi>s</mi></msub><mo>+</mo><mo>(</mo><msub><mi>R</mi><mi>r</mi></msub><mo>+</mo><msub><mi>R</mi><mrow><mi>c</mi><mi>b</mi></mrow></msub><mo>)</mo><mo>/</mo><mi>s</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>s</mi></msub><mo>+</mo><msub><mi>X</mi><mi>r</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000836293350000013.GIF" wi="925" he="163" /></maths>式(3)中:p为极对数,U<sub>1</sub>为定子相电压有效值;s=(ω<sub>1</sub>-ω<sub>r</sub>)/ω<sub>1</sub>为转差率;f<sub>N</sub>为电网频率;R<sub>s</sub>为定子绕组电阻;R<sub>r</sub>为转子绕组电阻;X<sub>s</sub>为定子绕组漏抗;X<sub>r</sub>为转子绕组漏抗;R<sub>cb</sub>为Crowbar保护电阻,通常取值为0.6p.u.~1.5p.u.,当电网发生不对称故障时,由于DFIG定子绕组中性点不接地,定子电压只存在正序分量和负序分量,无零序分量,故障导致的电压正序分量降低会使T<sub>eav</sub>减小,引起转子加速;电压的负序分量使电磁转矩中出现2倍频分量,但对T<sub>eav</sub>影响很小,分析不对称故障时,用定子电压正序分量有效值替代式(3)中U<sub>1</sub>,式(3)表明:电网扰动导致的电压跌落越深,电磁转矩平均值T<sub>eav</sub>下降幅度越大;2).Crowbar保护动作引起的转子转速变化正常运行的双馈风电机组可分为超同步和亚同步两种运行状态,超同步为发电机角速度大于旋转磁场角速度,亚同步即为发电机角速度小于旋转磁场角速度,当双馈风电机组处于超同步运行状态时,距超速保护定值ω<sub>rlim</sub>裕度较低,受电网故障扰动后更容易触发超速保护动作,由于所分析的过程时间较短,0&lt;分析的过程时间&lt;1s,再考虑到风力机的惯性较大,可假设T<sub>m</sub>在受扰后的暂态过程中保持不变,受扰前,T<sub>m</sub>与T<sub>e</sub>相等,转子转速保持不变,受扰后,定子电压正序分量降低,引起T<sub>eav</sub>发生跳变;3).转子加速至超速保护动作的时间Tcr由任何初态ω<sub>r0</sub>开始,Crowbar保护动作后转子转速随时间的变化可由以下积分式表示,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>&omega;</mi><mi>N</mi></msub><msub><mi>T</mi><mi>j</mi></msub></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mrow><mo>(</mo><msubsup><mi>T</mi><mi>m</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>T</mi><mrow><mi>e</mi><mi>a</mi><mi>v</mi></mrow></msub><msup><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>*</mo></msup><mo>)</mo></mrow><mi>d</mi><mi>u</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000836293350000021.GIF" wi="884" he="126" /></maths>当ω<sub>r</sub>达到超速保护定值ω<sub>rlim</sub>时,对应的时间即为转子加速至超速保护动作的时间T<sub>cr</sub>,于是有:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mrow><mi>r</mi><mi>lim</mi></mrow></msub><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>&omega;</mi><mi>N</mi></msub><msub><mi>T</mi><mi>j</mi></msub></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>T</mi><mrow><mi>c</mi><mi>r</mi></mrow></msub></msubsup><mrow><mo>(</mo><msubsup><mi>T</mi><mi>m</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>T</mi><mrow><mi>e</mi><mi>a</mi><mi>v</mi></mrow></msub><msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>*</mo></msup><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000836293350000022.GIF" wi="854" he="124" /></maths>式(5)表明:Crowbar动作后电磁转矩的平均值T<sub>eav</sub>越小,T<sub>cr</sub>越小;机组的初始转速ω<sub>r0</sub>越大,T<sub>cr</sub>也越小,考虑到在Crowbar保护动作的短暂过程中T<sub>eav</sub>(t)<sup>*</sup>变化不大,设为常数T<sub>eav</sub><sup>*</sup>,则转子加速至超速保护动作的时间T<sub>cr</sub>近似表达为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>c</mi><mi>r</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mrow><mi>r</mi><mi>lim</mi></mrow></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&omega;</mi><mi>N</mi></msub><mrow><mo>(</mo><msubsup><mi>T</mi><mi>m</mi><mo>*</mo></msubsup><mo>-</mo><msup><msub><mi>T</mi><mrow><mi>e</mi><mi>a</mi><mi>v</mi></mrow></msub><mo>*</mo></msup><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000836293350000023.GIF" wi="766" he="150" /></maths>
地址 130021 吉林省长春市朝阳区人民大街4629号