发明名称 一种电力系统小干扰振荡稳定性分解式模态分析方法
摘要 本发明公开了一种电力系统小干扰振荡稳定性分解式模态分析方法,将常规的计算残差Ri的模态分析进行分解计算,从而实现物理意义清晰的模态分析计算结果,其首先对电力系统稳态数据和电力系统静态数据进行测量后,然后分别对通道因子、重构系数和参与因子进行计算后,得到残差Ri的模态分析进行分解计算结果,基于该结果,可清晰揭示控制器影响阻尼的内部机理,实现对附加阻尼控制器的整定,以有效抑制电力系统低频振荡。本发明的电力系统小干扰振荡稳定性分解式模态分析方法,能够清晰详细地给出指定控制器或装置是如何影响指定的电力系统机电振荡模态,从而为电力系统小干扰振荡稳定性分析和控制设计提供物理意义清晰的指导。
申请公布号 CN101917003A 申请公布日期 2010.12.15
申请号 CN201010240356.4 申请日期 2010.07.29
申请人 东南大学 发明人 杜文娟;王海风;陈中
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种电力系统小干扰振荡稳定性分解式模态分析方法,其特征在于:所述分析方法包括如下步骤:(1)通过数据采集和监控系统SCADA系统、能量管理系统EMS获得电力系统稳态数据:发电机机端电压、机端有功、母线有功和母线无功;(2)输入电力系统稳态数据:发电机机端电压、机端有功、母线有功和母线无功;(3)输入电力系统静态数据:电力网络拓扑数据、线路阻抗导纳数据、变压器阻抗变比数据;(4)发电机组固有数据:发电机内部电抗数据、励磁系统数据;(5)利用潮流计算工具和初值计算工具,得到含指定控制器或装置的开环系统的线性化矩阵: <mrow> <mi>&Delta;</mi> <mover> <mi>X</mi> <mo>&CenterDot;</mo> </mover> <mo>=</mo> <mi>A&Delta;X</mi> <mo>+</mo> <mi>B&Delta;u</mi> </mrow>Δy=CΔX以及: <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>&Delta;</mi> <mover> <mi>&delta;</mi> <mo>&CenterDot;</mo> </mover> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;</mi> <mover> <mi>&omega;</mi> <mo>&CenterDot;</mo> </mover> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;</mi> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>&omega;</mi> <mn>0</mn> </msub> <mi>I</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>A</mi> <mn>22</mn> </msub> </mtd> <mtd> <msub> <mi>A</mi> <mn>23</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mn>31</mn> </msub> </mtd> <mtd> <msub> <mi>A</mi> <mn>32</mn> </msub> </mtd> <mtd> <msub> <mi>A</mi> <mn>33</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>&Delta;&delta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&omega;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;Z</mi> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mn>3</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>&Delta;u</mi> </mrow>其中,X为系统状态变量,A为系统线性化矩阵,B为控制矩阵,C为输出矩阵,u为控制器输出信号,δ为发电机功角状态变量向量,ω为发电机转速状态变量向量,Z为系统其它状态变量向量,Δ为线性化算子,变量加点为微分算子,ω0I为对角阵,ω0为额定角速度;(6)计算指定控制器或装置向各台发电机提供的电磁转矩的通道因子Bj(λi),j=1,2,…N,即B(λi)的第j个分量,其中,N为系统中发电机的台数,λi= ξi±jωi为系统中指定分析的某个机电的第i个振荡模态,根据步骤(5)中得到:B(s)=B2+A23(sI A33) 1B3;(7)计算指定控制器或装置的输入反馈信号用各台发电机转速状态变量重构的重构系数γj(λi),j=1,2,…N: <mrow> <msub> <mi>&gamma;</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>C</mi> <msub> <mi>v</mi> <mi>i</mi> </msub> </mrow> <msub> <mi>v</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>j</mi> </mrow> </msub> </mfrac> </mrow>其中vi为相应于λi的右特征向量,vi2j为vi中对应于Δωj,j=1,2,…N的分量,j为第j台发电机,系统的输出方程为:Δy=CX(8)计算系统中各台发电机相应于指定振荡模态的参与因子Sij: <mrow> <msub> <mi>S</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mo>&PartialD;</mo> <msub> <mi>&lambda;</mi> <mi>i</mi> </msub> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mrow> <mo>(</mo> <mo>=</mo> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>j</mi> </mrow> </msub> <msub> <mi>v</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mi>N</mi> </mrow>其中,DjΔωj,j=1,2,…N为各台发电机上获得的额外阻尼转矩,wi为相应于λi的左特征向量,wi2j为wi中对应于Δωj,j=1,2,…N的分量;(9)通过对指定机电振荡模态的分解式模态分析方法计算残差: <mrow> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>S</mi> <mi>ij</mi> </msub> <msub> <mi>B</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>&gamma;</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow>(10)通过输出步骤(9)中得到的残差,安排指定控制器或装置在电力系统中的位置。
地址 210009 江苏省南京市江宁开发区东南大学路2号