发明名称 汽轮发电机转子匝间短路故障位置及短路匝数的判定方法
摘要 汽轮发电机转子匝间短路故障位置及短路匝数的判定方法,属于电力系统主设备继电保护领域,其特征在于,首先利用基于多回路模型的发电机转子匝间故障分析技术对故障进行计算,择优选择检测对象,利用特征谐波分量的有效值及相角确定转子故障类型及故障位置判据;然后基于计算结果拟合得到不同槽内绕组短路时定子分支故障特征谐波电流与短路匝数的函数关系;应用时只需从发电机机端采集定子一个分支的电流,通过快速傅里叶分析滤取分支电流的故障特征谐波分量的有效值及相角,根据故障类型和故障位置判据可确定故障位置,再通过求解定子分支故障特征谐波电流与短路匝数的函数关系所构成的方程即可确定短路匝数。
申请公布号 CN102087329B 申请公布日期 2013.05.22
申请号 CN201010554528.5 申请日期 2010.11.19
申请人 清华大学 发明人 郝亮亮;孙宇光;王祥珩
分类号 G01R31/06(2006.01)I 主分类号 G01R31/06(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 楼艮基
主权项 1.汽轮发电机转子匝间短路故障位置及短路匝数的判定方法,其特征在于,所述方法是在计算机中依次按以下步骤进行的:步骤(1),计算转子不同槽发生的不同匝数短路故障时定子支路电流的各次谐波的有效值和相角,步骤如下:步骤(1.1),向计算机输入的定子和转子所有回路电压方程的矩阵形式可表示为:<img file="FSB00000951506000011.GIF" wi="2130" he="868" /><img file="FSB00000951506000012.GIF" wi="1738" he="807" />其中:L<sub>1</sub>,…,L<sub>Q</sub>,…,L<sub>N</sub>是各定子支路自感,Q=1,2,…N,Q是定子支路序号,N是定子支路总数,L<sub>1d</sub>,…,L<sub>gd</sub>,…,L<sub>ld</sub>是各阻尼回路自感,gd=1d,2d,…ld,gd是阻尼回路序号,ld是阻尼回路总数,L<sub>f</sub>是励磁绕组回路自感,L<sub>fkL</sub>是励磁绕组故障附加回路自感,M<sub>Q,1</sub>,…,M<sub>Q,N</sub>是第Q条定子支路与其它各定子支路之间的互感,M<sub>Q,1d</sub>,…,M<sub>Q,gd</sub>,…,M<sub>Q,ld</sub>是定子第Q条支路与各阻尼回路gd之间的互感,gd=1d,2d,…,ld,M<sub>1,f</sub>,…,M<sub>Q,f</sub>,…,M<sub>N,f</sub>是各定子支路Q与励磁绕组回路f之间的互感,M<sub>1,fkL</sub>,…,M<sub>Q,fkL</sub>,…,M<sub>N,fkL</sub>是各定子支路Q与励磁绕组故障附加回路fkL之间的互感,M<sub>gd,1</sub>,…,M<sub>gd,Q</sub>,…,M<sub>gd,N</sub>是第gd个阻尼回路与各定子支路Q之间的互感,M<sub>gd,1d</sub>,…,M<sub>gd,ld</sub>是第gd个阻尼回路与其它各阻尼回路之间的互感,M<sub>1d,f</sub>,…,M<sub>gd,f</sub>,…,M<sub>ld,f</sub>是各阻尼回路gd与励磁绕组回路f之间的互感,M<sub>1d,fkL</sub>,…,M<sub>gd,fkL</sub>,…,M<sub>ld,fkL</sub>是各阻尼回路gd与励磁绕组故障附加回路fkL之间的互感,M<sub>f,1</sub>,…,M<sub>f,Q</sub>,…,M<sub>f,N</sub>是励磁绕组回路f与各定子支路Q之间的互感,M<sub>f,1d</sub>,…,M<sub>f,gd</sub>…,M<sub>f,ld</sub>是励磁绕组回路f与各阻尼回路gd之间的互感,M<sub>f,fkL</sub>是励磁绕组回路f与励磁绕组故障附加回路fkL之间的互感,M<sub>f,fkL</sub>=M<sub>fkL,f</sub>,M<sub>fkL,1</sub>,…,M<sub>fkL,Q</sub>,…,M<sub>fkL,N</sub>是励磁绕组故障附加回路fkL与各定子支路Q之间的互感,M<sub>fkL,1d</sub>,…,M<sub>fkL,gd</sub>…,M<sub>fkL,ld</sub>是励磁绕组故障附加回路fkL与各阻尼回路gd之间的互感,i<sub>1</sub>,…,i<sub>Q</sub>,…,i<sub>N</sub>是各定子支路Q电流,i<sub>1d</sub>,…,i<sub>gd</sub>,…,i<sub>ld</sub>是各阻尼回路gd电流,i<sub>f</sub>、i<sub>fkL</sub>分别是励磁绕组回路f和励磁绕组故障附加回路fkL的电流,i<sub>A</sub>,i<sub>B</sub>,i<sub>C</sub>是定子三相负载电流,r<sub>T</sub>、L<sub>T</sub>、u<sub>A</sub>′、u<sub>B</sub>′、u<sub>C</sub>′分别为折算到发电机一侧的变压器的电阻、电感和电网各相电压,U<sub>1</sub>,…,U<sub>Q</sub>,…,U<sub>N</sub>是各定子支路电压,U<sub>A</sub>,U<sub>B</sub>,U<sub>C</sub>是定子各相电压,r<sub>Q</sub>是定子支路电阻,r<sub>1d</sub>,…,r<sub>gd</sub>,…,r<sub>ld</sub>是各阻尼回路电阻,E<sub>ZF</sub>、R<sub>ZF</sub>分别为励磁系统电源的电动势和内电阻,r<sub>f</sub>为励磁绕组回路电阻,r<sub>fk</sub>为励磁绕组故障附加回路电阻,将上述矩阵形式的方程简记为:[U]=p{[L]·[I]}+[R]·[I]+[B′]其中:[U]是以列形式表示的定子支路电压是由U<sub>1</sub>,…U<sub>1</sub>,…,U<sub>N</sub>,各阻尼回路电压U<sub>1d</sub>,…,U<sub>gd</sub>,…,U<sub>ld</sub>=0,励磁绕组电压E<sub>ZF</sub>,励磁绕组故障附加回路电压U<sub>fkL</sub>=0以及定子侧负载电压U<sub>A</sub>,U<sub>B</sub>,U<sub>C</sub>共同组成的序列,[L]是矩阵形式表示的所有定子支路和转子各回路的自感或互感,[R]是矩阵形式表示的所有定子支路和转子各回路的电阻,[I]是矩阵形式表示的所有定子支路和转子各回路的电流,[B′]是电网三相电压的列形式,p是微分算子,由于上式中转子电压方程是回路电压,而定子电压方程是支路电压,为了处理方便,定转子电压方程都采用回路电压,这样就需要重新处理定子回路,将定子支路电压方程转换为定子回路电压方程,从而得到以定子回路电流和转子各回路电流为状态变量的同步发电机状态方程:p[I′]=[A]·[I′]+[B]其中:[A]=-[L′]<sup>-1</sup>·[R′]·[I′];[B]=[L′]<sup>-1</sup>·[U′]-[L′]<sup>-1</sup>·[H]·[B′][L′]=[H]·[L]·[H<sup>T</sup> ];[R′]=p[L′]+[H]·[R]·[H<sup>T</sup> ],[I′]=[H<sup>T</sup>]<sup>-1</sup>·[I],[I′]包括定子和转子回路电流,[I]包括定子支路电流、转子回路电流,[U′]=[H]·[U],[U′]包括定子和转子回路电压,[U]包括定子支路电压、转子回路电压,[H]为定子支路对回路的变换阵,步骤(1.2),计算定子和转子各回路的参数,步骤(1.3),采用四阶龙格-库塔法数值方法对步骤(1,1)的微分方程组进行求解,即可求得转子不同槽数m下发生的不同匝数短路故障时定子回路电流的二次和四次谐波环流的有效值I<sub>m</sub>和二次和四次谐波环流的相角α<sub>m</sub>,m=1,2,3,…,用I<sub>m,μ2</sub>,I<sub>m,μ4</sub>,α<sub>m,μ2</sub>,α<sub>m,</sub><sub>μ4</sub>表示,步骤(1.4),不同槽数m下I<sub>m,μ2</sub>/I<sub>m,μ4</sub>之值和α<sub>m,μ2</sub>-α<sub>m,μ4</sub>的值分别表示为K<sub>m</sub>和S<sub>m</sub>,对于发生在第m槽的短路故障,不论短路匝数N为多少,K<sub>m</sub>和S<sub>m</sub>为常数;步骤(2):拟合不同槽内绕组短路时定子分支故障特征谐波电流与短路匝数的函数关系,步骤如下:基于步骤1的计算结果分别得到转子各槽发生故障时定子分支二、四次谐波电流与短路匝数的关系;通过一次多项式拟合将转子第m槽发生故障时引起的定子支路电流的二、四次谐波I<sub>μ2</sub>、I<sub>μ4</sub>分别表示为以短路匝数N为自变量的函数:I<sub>m,μ2</sub>=fm<sub>,μ2</sub>(N)、I<sub>m,μ4</sub>=f<sub>m,μ4</sub>(N),m是槽的序号,对于A1553实验电机,以第5槽短路时,2次谐波电流I<sub>5,μ2</sub>与短路匝数N的关系:I<sub>5,μ2</sub>=f<sub>5,μ2</sub>(N)=0.0121+0.03088N,步骤(3),判定故障类型及位置,步骤如下:步骤(3.1),在发电机机端采集定子一个分支的电流,利用傅立叶滤波算法计算出二次和四次谐波分量的有效值I′<sub>μ2</sub>、I′<sub>μ4</sub>与相角α′<sub>μ2</sub>、α′<sub>μ4</sub>,求出K′=I′<sub>μ2</sub>/I′<sub>μ4</sub>、S′=α′<sub>μ2</sub>-α′<sub>μ4</sub>,步骤(3.2),若<img file="FSB00000951506000031.GIF" wi="255" he="140" />或<img file="FSB00000951506000032.GIF" wi="262" he="154" />时,则为同槽内部故障,转步骤(3.4),否则转步骤(3.3),步骤(3.3),若满足以下两种情况之一:若<img file="FSB00000951506000041.GIF" wi="392" he="166" />且K′介于K<sub>m</sub>和K<sub>m+1</sub>之间,或若<img file="FSB00000951506000042.GIF" wi="392" he="130" />且S′介于S<sub>m</sub>和S<sub>m+1</sub>之间,则为第m槽和第m+1槽之间的短路故障,步骤(3.4)若判断为第m槽内的短路故障,利用求解步骤(2)中得到的函数I<sub>μ2</sub>=f<sub>m,μ2</sub>(N)或I<sub>μ4</sub>=f<sub>m,μ4</sub>(N),列写方程:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>I</mi><mrow><mi>&mu;</mi><mn>2</mn></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>f</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>N</mi><mo>&le;</mo><msub><mi>N</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>上式中N<sub>m</sub>为发电机励磁绕组第m槽匝数,通过求解上述方程,计算出短路匝数,若判断为第m槽和第m+1槽之间的故障,需分别对二、四次谐波列出方程,并加上两槽导体匝数的约束条件,即:<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>I</mi><mrow><mi>&mu;</mi><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msqrt><msubsup><mi>f</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>f</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>N</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>f</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>f</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><msub><mi>N</mi><mn>2</mn></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mi>&alpha;</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>&alpha;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>&mu;</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></msqrt></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><mrow><mi>&mu;</mi><mn>2</mn></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msqrt><msubsup><mi>f</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>4</mn></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>f</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>&mu;</mi><mn>4</mn></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>N</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>f</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>4</mn></mrow></msub><mrow><mo>(</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>f</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>&mu;</mi><mn>4</mn></mrow></msub><mrow><mo>(</mo><msub><mi>N</mi><mn>2</mn></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mi>&alpha;</mi><mrow><mi>m</mi><mo>,</mo><mi>&mu;</mi><mn>4</mn></mrow></msub><mo>-</mo><msub><mi>&alpha;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>&mu;</mi><mn>4</mn></mrow></msub><mo>)</mo></mrow></msqrt></mtd></mtr><mtr><mtd><msub><mi>N</mi><mn>1</mn></msub><mo>&le;</mo><msub><mi>N</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>N</mi><mn>2</mn></msub><mo>&le;</mo><msub><mi>N</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced></math>]]></maths>上式中N<sub>1</sub>、N<sub>2</sub>分别为两槽的短路匝数,解该方程即可计算出两槽的短路匝数,若计算出的短路匝数不足1,或方程无正值解,则说明发电机无转子匝间短路故障。
地址 100084 北京市海淀区清华园1号