发明名称 一种基于暂态量特征的超/特高压输电线路故障测距方法
摘要 本发明涉及一种基于暂态量特征的超/特高压输电线路故障测距方法。现有的故障测距精确度不高。本发明首先采样故障系统保护安装处的电流和母线电压并计算故障分量系统中的电流和电压。其次计算系统等效电感并消去工频分量。然后通过计算各次暂态频率分量的幅频特性找到各次主频频率值。最后根据各次主频计算故障距离。本发明频率提取误差小、故障测距精确度高。
申请公布号 CN102508114A 申请公布日期 2012.06.20
申请号 CN201110329463.9 申请日期 2011.10.26
申请人 浙江大学 发明人 王慧芳;高效海;何奔腾
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.一种基于暂态量特征的超/特高压输电线路故障测距方法,其特征在于该方法包括以下步骤:步骤1.采样故障系统保护安装处的电流i(n)和母线电压u(n);步骤2.计算故障分量系统中的电流i<sub>g</sub>(n)和电压u<sub>g</sub>(n);具体是:正常运行系统的电流为负荷电流为i<sub>f</sub>(n),电压为额定电压u<sub>e</sub>(n),根据叠加原理,可以用故障系统中采集到的电压电流减去正常运行系统的电压电流,获得故障分量系统中的电压电流,即u<sub>g</sub>(n)=u(n)-u<sub>e</sub>(n),i<sub>g</sub>(n)=i(n)-i<sub>f</sub>(n);步骤3.计算系统等效电感L<sub>N</sub>;具体是:根据故障分量系统有:<img file="FDA0000102275120000011.GIF" wi="424" he="111" />因此通过任意两组不相关的采样数据,就可以求得系统的等效电阻R<sub>N</sub>和电感L<sub>N</sub>;为了降低误差,利用i<sub>g</sub>(n)、u<sub>g</sub>(n)采样数据,通过最小二乘法拟合,得到L<sub>N</sub>值;步骤4.为突出故障电流各次暂态主频频率分量的特性,消去工频分量;具体是:对电流i<sub>g</sub>(n)进行积分,变成I(n)序列;然后再对I(n)序列进行积分,变成II(n)序列,最后通过以下迭代公式<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>II</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>i</mi><mi>g</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>I</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>i</mi><mi>g</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>s</mi></msub><mo>+</mo><mi>I</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>II</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>s</mi></msub><mo>+</mo><mi>II</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>X</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup></mfrac><mi>i</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mi>II</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>可以计算到消去工频分量影响的X(n)序列,其中T<sub>s</sub>为采样周期,ω<sub>0</sub>为工频频率;步骤5.为减少直流分量的影响,对X(n)序列求导得到序列Y(n);求导公式为:<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>Y</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mi>X</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>Y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>X</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mi>X</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msub><mi>T</mi><mi>s</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>步骤6.对Y(n)序列进行离散傅里叶变换,通过计算各次暂态频率分量的幅频特性找到各次主频频率值;傅里叶变换公式为:<img file="FDA0000102275120000022.GIF" wi="407" he="120" />其中N为数据长度,<img file="FDA0000102275120000023.GIF" wi="230" he="93" />幅频特性为:|i(k)|<sup>2</sup>=R[i(k)]<sup>2</sup>+I[i(k)]<sup>2</sup>,k对应的频率为<img file="FDA0000102275120000024.GIF" wi="223" he="121" />而ω=2πf;幅频特性最大的频率值为一次主频率ω<sub>1</sub>,其次的为二次主频频率ω<sub>2</sub>,以此类推可得m次主频ω<sub>m</sub>;步骤7.根据各次主频计算故障距离;根据公式<img file="FDA0000102275120000025.GIF" wi="555" he="108" />代入各次主频频率,即可计算出故障距离。
地址 310027 浙江省杭州市西湖区浙大路38号