发明名称 一种考虑不确定性影响的电力系统故障分析方法
摘要 本发明公开了一种考虑不确定性影响的电力系统故障分析方法,用于对各线路开断故障引起系统过负荷的可能性进行静态分析,属于电力系统分析与调度技术领域。本发明针对现有电力系统静态安全分析技术中所采用的故障行为指标无法适用于存在负荷波动和间歇性能源出力等不确定性因素的情况,提出了一种全新的故障行为指标,该指标将电网运行时不确定因素对支路功率的影响和支路的负载率结合起来,从而可有效地表征电网运行存在不确定因素情况下的线路开断故障引起系统过负荷的可能性,并基于该指标进行进一步地系统安全性分析、调度。
申请公布号 CN104240039B 申请公布日期 2017.04.26
申请号 CN201410506194.2 申请日期 2014.09.26
申请人 东南大学;国家电网公司;中国电力科学研究院;国网宁夏电力公司 发明人 高山;陈跃彬;姚建国;杨胜春;耿天翔;摆世彬;王珂;曾丹
分类号 G06Q10/06(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/06(2012.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种考虑不确定性影响的电力系统故障分析方法,用于对各线路开断故障引起系统过负荷的可能性进行静态分析,其特征在于,用故障行为指标PI<sub>k</sub>来表征所述电力系统中的第k条支路出现开断故障引起系统过负荷的可能性,故障行为指标值越大,引起系统过负荷的可能性越大,故障行为指标按照下式得到:<maths num="0001"><math><![CDATA[<mrow><msub><mi>PI</mi><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>w</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mo>|</mo><mfrac><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mi>l</mi></mrow></msub><msub><mi>P</mi><mrow><mi>l</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mfrac><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mo>|</mo><msub><mi>P</mi><mrow><mi>k</mi><mi>l</mi></mrow></msub><mo>|</mo><mo>&gt;</mo><msub><mi>P</mi><mrow><mi>l</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001177157160000011.GIF" wi="628" he="151" /></maths>式中,<img file="FDA0001177157160000012.GIF" wi="58" he="70" />表示所述第k条支路断开后,所述电力系统中负载率大于w的m条其它支路中的第l条支路的有功功率期望值;P<sub>lmax</sub>表示所述第l条支路的传输容量;F(|P<sub>kl</sub>|>P<sub>lmax</sub>)表示所述第l条支路的越线概率;w为预设的负载率阈值,其取值范围为[0.5,1];所述<img file="FDA0001177157160000013.GIF" wi="82" he="68" />F(|P<sub>kl</sub>|>P<sub>lmax</sub>)按照以下方法得到:步骤1、用普通直流潮流计算方法算出所述电力系统在正常运行情况下的潮流分布,得到当前运行点上的节点阻抗矩阵X、关联矩阵A、支路功率与节点电压相角的关系矩阵T、支路功率与节点注入功率的关系矩阵H;步骤2、根据发电机组出力和负荷功率的期望值和方差,求出节点注入有功功率随机变量的一阶半不变量P<sup>(1)</sup>和二阶半不变量P<sup>(2)</sup>:设节点i的发电机注入有功功率的期望值和方差分别为u<sub>1</sub>和<img file="FDA0001177157160000014.GIF" wi="83" he="68" />负荷注入有功功率的期望值和方差分别为u<sub>2</sub>和<img file="FDA0001177157160000015.GIF" wi="83" he="69" />则节点i的注入有功功率的一阶和二阶半不变量分别为P<sub>i</sub><sup>(1)</sup>=u<sub>1</sub>+u<sub>2</sub>,<img file="FDA0001177157160000016.GIF" wi="315" he="65" />步骤3、利用基于补偿法的支路开断模拟方法获取所述第k条支路断开后,其他支路有功功率与节点注入有功功率的关系矩阵H';然后由<img file="FDA0001177157160000017.GIF" wi="230" he="70" />其中P为节点注入有功功率期望值,得到所述第k条支路断开后,所述电力系统中其他支路的有功功率期望值矩阵<img file="FDA0001177157160000018.GIF" wi="89" he="71" />选出其中负载率大于w的支路,从H'提取出负载率大于w的m条支路的有功功率与节点注入有功功率的关系矩阵,组成新的关系矩阵H<sub>k</sub>;步骤4、由节点注入有功功率的一阶半不变量和二阶半不变量,根据公式<img file="FDA0001177157160000019.GIF" wi="322" he="71" />v=1,2求出第k条支路开断后,其他支路中负载率大于w的m条支路的有功功率P<sub>k</sub>的一阶半不变量<img file="FDA00011771571600000110.GIF" wi="76" he="63" />和二阶半不变量<img file="FDA00011771571600000111.GIF" wi="102" he="70" />进而得到负载率大于w的m条支路的有功功率P<sub>k</sub>的期望值和方差,其中<img file="FDA0001177157160000021.GIF" wi="450" he="63" />l=1,2,…,m;步骤5、根据所得到的支路功率P<sub>kl</sub>的期望值和方差,l=1,2,…,m,求出其概率密度函数<img file="FDA0001177157160000022.GIF" wi="556" he="165" />从而可得其越线概率:<maths num="0002"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mo>|</mo><msub><mi>P</mi><mrow><mi>k</mi><mi>l</mi></mrow></msub><mo>|</mo><mo>&gt;</mo><msub><mi>P</mi><mrow><mi>l</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mrow><msub><mi>&sigma;</mi><mrow><mi>k</mi><mi>l</mi></mrow></msub><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt></mrow></mfrac><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>l</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow><msub><mi>P</mi><mrow><mi>l</mi><mi>max</mi></mrow></msub></msubsup><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>k</mi><mi>l</mi></mrow></msub><mo>-</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mi>l</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mrow><mi>k</mi><mi>l</mi></mrow><mn>2</mn></msubsup></mrow></mfrac></mrow></msup><msub><mi>dP</mi><mrow><mi>k</mi><mi>l</mi></mrow></msub><mo>.</mo></mrow>]]></math><img file="FDA0001177157160000023.GIF" wi="980" he="159" /></maths>
地址 210096 江苏省南京市四牌楼2号