发明名称 一种计及风电不确定性的暂态稳定风险在线评估方法
摘要 本发明涉及一种计及风电不确定性的暂态稳定风险在线评估方法。该方法包括下述步骤:1)确定风电功率超短期预测值及其置信区间估计;2)确定暂态稳定概率计算的基础条件;3)确定暂态稳定故障概率;4)判断是否发生暂态失稳;5)对暂态稳定故障后果严重度进行评估;6)对暂态失稳后果严重程度进行评估;7)计算暂态稳定风险评估指标;8)判断是否所有线路暂态风险评估指标都已计算;9)输出暂态稳定风险评估指标值;本发明的方法对输电线路乃至整个系统进行暂态稳定评估,很好地把系统暂态稳定的概率和造成的后果有机结合起来,调度运行人员依据得到的多种暂态稳定风险指标进行系统运行决策。
申请公布号 CN102842105B 申请公布日期 2015.04.01
申请号 CN201210245482.8 申请日期 2012.07.09
申请人 中国电力科学研究院;江苏省电力公司;国家电网公司 发明人 苏峰;严剑峰;于之虹;任玲玉;李海峰;罗建裕;李汇群
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 北京安博达知识产权代理有限公司 11271 代理人 徐国文
主权项 一种计及风电不确定性的暂态稳定风险在线评估方法,其特征在于,所述方法包括下述步骤:1)确定风电功率超短期预测值及其置信区间估计;2)确定暂态稳定概率计算的基础条件;3)确定暂态稳定故障概率;4)判断是否发生暂态失稳;5)对暂态稳定故障后果严重度进行评估;6)对暂态失稳后果严重程度进行评估;7)计算暂态稳定风险评估指标;8)判断是否所有线路暂态风险评估指标都已计算;9)输出暂态稳定风险评估指标值;所述步骤1)中,基于混沌时间序列模型,利用广域监测系统WAMS实时量测风电功率数据和风电功率历史数据,计算得到风电功率预测值P<sub>pred</sub>(t);根据功率预测误差概率分布计算得到风电功率预测值P<sub>pred</sub>(t)的置信区间;所述超短期指的是15分钟;所述混沌时间序列模型是基于相空间重构,利用加权一阶局域预测法得到风电功率超短期预测公式;采用风电场相对误差统计功率预测误差概率分布;所述风电场相对误差e<sub>t</sub>为t时刻的风电功率预测值P<sub>pred</sub>(t)与风电功率实测值P<sub>meas</sub>(t)的差,即用下述(19)式表示:e<sub>t</sub>=P<sub>pred</sub>(t)‑P<sub>meas</sub>(t)    (19);所述计算得到风电功率超短期预测值的置信区间包括下述步骤:a、判断风电功率预测值P<sub>pred</sub>(t)的所属功率区段;b、查找所述功率区段对应的误差概率密度曲线;c、确定功率概率密度;d、寻找累计概率大于或等于置信区间度1‑α的区间作为置信区间;所述步骤a中,所述功率区段包括等分的十个区段;所述步骤b中,所述误差概率密度曲线呈正态分布曲线,所述误差概率密度函数用下述(23)式表示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt><msub><mi>&sigma;</mi><mi>f</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>&mu;</mi><mi>f</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msup><msub><mi>&sigma;</mi><mi>f</mi></msub><mn>2</mn></msup></mrow></mfrac></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>23</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000593105460000021.GIF" wi="1289" he="204" /></maths>式中:μ<sub>f</sub>为随机变量的期望值;σ<sub>f</sub>为随机变量的标准差;所述步骤c中,根据(23)式确定各功率区段的功率概率密度;暂态稳定概率包括线路故障发生概率、线路故障类型概率和线路故障位置概率;若线路发生的故障均为永久性故障,且两条线路同时发生故障的概率不计,则线路L<sub>k</sub>上发生故障的概率用下述(24)式表示:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>L</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&lambda;</mi><mi>k</mi></msub><mi>t</mi></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000593105460000022.GIF" wi="1272" he="77" /></maths>式中:P<sub>r</sub>(L<sub>k</sub>)为线路发生的故障概率;λ<sub>k</sub>表示时间段内的线路L<sub>k</sub>的故障频率;t表示故障持续时间;m为自然数;所述故障类型按严重程度分为三相接地、两相接地、相间短路和单相接地;故障类型为C<sub>j</sub>的故障发生的频率为f<sub>j</sub>,则其发生的概率用下述(25)式表示:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>f</mi><mi>j</mi></msub><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></munderover><msub><mi>f</mi><mi>j</mi></msub></mrow></mfrac><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>4</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>25</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000593105460000023.GIF" wi="1198" he="192" /></maths>且有:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></munderover><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>26</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000593105460000024.GIF" wi="1073" he="139" /></maths>式中:P<sub>r</sub>(C<sub>j</sub>)为线路的故障类型概率;f<sub>j</sub>为线路故障发生的频率;若线路的总长度为L,故障点到线路首端的距离为L<sub>d</sub>,定义变量D:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><mfrac><msub><mi>L</mi><mi>d</mi></msub><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>27</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000593105460000025.GIF" wi="959" he="118" /></maths>所述线路故障位置概率用下述(28)式表示:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>d</mi></msub></munderover><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>D</mi><mi>h</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>28</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000593105460000026.GIF" wi="1046" he="142" /></maths>式中:P<sub>r</sub>(D<sub>h</sub>)表示线路上第h个故障点发生故障的概率;n<sub>d</sub>表示线路上离散故障点总数。
地址 100192 北京市海淀区清河小营东路15号