发明名称 一种基于广域测量信息的同调机群快速识别方法
摘要 本发明涉及一种基于广域测量信息的同调机群快速识别方法,包括:记录分群数据窗内各发电机的功角信息和采样时间;采用差分法计算发电机转子角速度和转子角加速度;对发电机的转子角加速度进行三角函数拟合;对各电机的拟合后得到的三角函数求导,计算局部极值点进而确定各机组角点;判断两两机组的相应角点是否匹配;遍历所有发电机,更新匹配角点数矩阵;计算两两发电机的匹配角点数分别占两发电机角点的比例,判断两发电机是否同调;遍历所有发电机,更新同调矩阵,同调分群结束。本发明可实现电力系统同调机群的在线识别,为电力系统的动态等值和解列控制提供充裕的时间。
申请公布号 CN106487026A 申请公布日期 2017.03.08
申请号 CN201610727604.5 申请日期 2016.08.25
申请人 天津大学 发明人 张艳霞;李莹
分类号 H02J3/24(2006.01)I 主分类号 H02J3/24(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 程毓英
主权项 一种基于广域测量信息的同调机群快速识别方法,设系统中有N台发电机,分群信息存储在同调矩阵C<sub>N×N</sub>中,矩阵中元素c<sub>ij</sub>反映发电机i、j间的同调关系,当c<sub>ij</sub>=1时,发电机i、j同调;当c<sub>ij</sub>=0时,发电机i、j非同调;列向量Q<sub>1×N</sub>存储分群数据窗内的各发电机角点数,矩阵P<sub>N×N</sub>存储各发电机间的匹配角点数。该快速识别方法包括如下的步骤:1)同调矩阵C<sub>N×N</sub>和匹配角点数矩阵P<sub>N×N</sub>初始化;2)记录t时刻分群数据窗内各发电机的功角信息和采样时间,保存至功角矩阵D<sub>m×(N+1)</sub>,D<sub>m×(N+1)</sub>首列为采样时间,第i+1列为发电机i的功角信息;3)采用差分法计算发电机转子角速度和转子角加速度,分别存入角速度矩阵W<sub>(m‑1)×(N+1)</sub>、角加速度矩阵S<sub>(m‑2)×(N+1)</sub>中;4)设定拟合周期,对发电机的转子角加速度进行三角函数拟合;5)对各电机的拟合后得到的三角函数求导,计算局部极值点进而确定各机组角点,各发电机的角点对应时间暂存至矩阵T<sub>N×h</sub>,幅值暂存至矩阵E<sub>N×h</sub>,更新角点数列向量Q<sub>1×N</sub>;6)当角点数列向量值均不小于5个时,对T<sub>N×h</sub>和E<sub>N×h</sub>中各发电机角点的时间差和幅值差进行比较,若i与发电机j存在角点x、y满足如下公式<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&Delta;t</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mo>|</mo><msub><mi>t</mi><mrow><mi>i</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>|</mo><mo>&lt;</mo><msub><mi>&epsiv;</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;e</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mo>|</mo><msub><mi>e</mi><mrow><mi>i</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>e</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>|</mo><mo>&lt;</mo><msub><mi>&epsiv;</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001092438610000011.GIF" wi="365" he="135" /></maths>式中,t<sub>ix</sub>、t<sub>jy</sub>分别为发电机i、j的第x、y个角点对应的时间;e<sub>ix</sub>、e<sub>jy</sub>为发电机i、j的第x、y个角点对应的幅值;ε<sub>1</sub>为时间阈值;ε<sub>2</sub>为幅值阈值,则两机组的相应角点匹配,匹配角点数矩阵p<sub>ij</sub>、p<sub>ji</sub>加1;遍历所有发电机,更新匹配角点数矩阵P<sub>N×N</sub>;7)计算发电机i、j匹配角点数分别占两发电机角点的比例z<sub>i</sub>=p<sub>ij</sub>/q<sub>i</sub>、z<sub>j</sub>=p<sub>ji</sub>/q<sub>j</sub>,当(z<sub>i</sub>+z<sub>j</sub>)/2&gt;0.5时两发电机同调,同调矩阵c<sub>ij</sub>=1;否则,c<sub>ij</sub>=0;遍历所有发电机,更新同调矩阵C<sub>N×N</sub>,同调分群结束。
地址 300072 天津市南开区卫津路92号