主权项 |
一种基于广域测量信息的同调机群快速识别方法,设系统中有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>Δ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><</mo><msub><mi>ϵ</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>Δ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><</mo><msub><mi>ϵ</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>0.5时两发电机同调,同调矩阵c<sub>ij</sub>=1;否则,c<sub>ij</sub>=0;遍历所有发电机,更新同调矩阵C<sub>N×N</sub>,同调分群结束。 |