发明名称 暂态安全域最短半径估计方法
摘要 本发明公开了一种暂态安全域最短半径估计方法,涉及电力系统安全分析技术领域,包括:S1:计算电力系统暂态稳定性量化指标;S2:计算稳定性量化指标分别对发电机节点注入向量的灵敏度和负荷节点注入向量的灵敏度;S3:找出当前运行处节点注入向量的最危险增长方向;S4:在最危险增长方向上增加负荷及发电出力增长步长,得到增加后的节点注入向量;S5:在增加后的节点注入向量处重复步骤S1,计算稳定量化指标,并判断系统是否稳定,如果系统稳定,重复步骤S2至步骤S5,直至系统失去稳定,得到电力系统的负荷裕度、发电裕度。本发明实现了对电力系统的在线稳定量化分析,且计算复杂度低,不存在数值收敛性问题。
申请公布号 CN102664405B 申请公布日期 2014.04.02
申请号 CN201210135066.2 申请日期 2012.04.28
申请人 清华大学 发明人 何光宇;李嘉;刘铠诚;常乃超
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 王莹
主权项 1.一种暂态安全域最短半径估计方法,其特征在于,包括以下步骤:S1:计算电力系统暂态稳定性量化指标;S2:计算所述暂态稳定性量化指标分别对发电机节点注入向量的灵敏度和负荷节点注入向量的灵敏度;S3:找出当前运行状态下节点注入向量的最危险增长方向,所述最危险增长方向为发电机节点注入向量的灵敏度和负荷节点注入向量的灵敏度的最大值元素的方向;S4:在所述最危险增长方向上增加负荷及发电出力增长步长,得到增加后的节点注入向量;S5:在所述增加后的节点注入向量处重复步骤S1,计算所述暂态稳定性量化指标,并判断系统是否稳定,如果系统稳定,重复步骤S2至步骤S5,直至系统失去稳定,得到电力系统的负荷裕度、发电裕度;所述步骤S1中选取基于单机能量函数的混合方法SHM提供的电力系统暂态稳定性量化指标:<![CDATA[<math><mrow><mi>&eta;</mi><mo>=</mo><msubsup><mi>V</mi><mrow><mi>PE</mi><mo>,</mo><mi>i</mi></mrow><mi>max</mi></msubsup><mo>-</mo><msub><mi>V</mi><mrow><mi>PE</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>cl</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mrow><mi>KE</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>cl</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,t<sub>cl</sub>为故障切除时刻,V<sub>PE,i</sub>(t<sub>cl</sub>)、V<sub>KE,i</sub>(t<sub>cl</sub>)分别为故障切除时刻的发电机i的势能、动能,<img file="FDA0000395427240000012.GIF" wi="105" he="76" />为持续故障下发电机i的势能最大值,所述持续故障是指在时刻t<sub>ins</sub>,t<sub>ins</sub>&gt;t<sub>cl</sub>重新插入故障,发电机i的势能所能达到的最大值,t<sub>ins</sub>满足:KE<sub>total</sub>(t<sub>ins</sub>)=min KE<sub>total</sub>(t),t&gt;t<sub>cl</sub>其中KE<sub>total</sub>(t)为时刻t系统总的动能,<img file="FDA0000395427240000013.GIF" wi="535" he="132" />其中,M<sub>i</sub>为发电机i的转动惯量常数,<img file="FDA0000395427240000014.GIF" wi="57" he="67" />为发电机i相对于系统惯性中心的角速度,对于故障F,注入空间中的暂态安全域定义为:Ω={P|<sub>η</sub>(P)≥ε,<sub>g</sub>(x,P)≤0,P∈R<sub>n</sub>}P为节点注入向量,P=[P<sub>1</sub>,P<sub>2</sub>,...,P<sub>n</sub>]<sup>T</sup>,其中P<sub>i</sub>表示各节点的注入功率,i=1,2,...,n,η(P)为故障F下的暂态稳定性量化指标,g(x,P)≤0为运行约束不等式,ε为门槛值,R<sub>n</sub>表示n维实数向量;若前系统运行状态下的节点注入为P<sup>0</sup>,<img file="FDA0000395427240000021.GIF" wi="456" he="79" />相应的针对故障F的SHM稳定裕度指标为:η<sup>0</sup>=η(P<sup>0</sup>)其中<![CDATA[<math><mrow><msup><mi>P</mi><mn>0</mn></msup><mo>=</mo><msup><mrow><mo>&lsqb;</mo><msubsup><mi>P</mi><mi>G</mi><mn>0</mn></msubsup><mo>,</mo><msubsup><mi>P</mi><mi>L</mi><mn>0</mn></msubsup><mo>&rsqb;</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><![CDATA[<math><mrow><msubsup><mi>P</mi><mi>G</mi><mn>0</mn></msubsup><mo>=</mo><msup><mrow><mo>&lsqb;</mo><msubsup><mi>P</mi><mrow><mi>G</mi><mo>,</mo><mn>1</mn></mrow><mn>0</mn></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mi>G</mi><mo>,</mo><mn>2</mn></mrow><mn>0</mn></msubsup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msubsup><mi>P</mi><mrow><mi>G</mi><mo>,</mo><mi>m</mi></mrow><mn>0</mn></msubsup><mo>&rsqb;</mo></mrow><mi>T</mi></msup></mrow></math>]]></maths>为发电机节点注入向量,<img file="FDA0000395427240000024.GIF" wi="474" he="91" />为负荷节点注入向量,v+m=n。
地址 100084 北京市海淀区清华园北京100084-82信箱
您可能感兴趣的专利