发明名称 利用暂态主频分量的配网电缆单相接地故障测距方法
摘要 传统的线路双曲函数模型多用于工频稳态量的相量计算,本文利用拉氏变换将该模型进行了改进,结合配网电缆单相接地故障后的零模网络,提出了一种利用暂态主频量的单端测距算法,并分别通过频域法和时域法的求解,得到了故障距离。其中,频域法消除了过渡电阻的影响,时域法利用采样点的冗余性进行了结果优化。大量的EMTP仿真实验结果,验证了该方法的正确性,且不受过渡电阻、中性点运行方式、故障初始角等因素的影响,其测距的最大相对误差小于0.231%,平均测距误差小于20米,能够满足实际工程需求。
申请公布号 CN103941147B 申请公布日期 2016.08.17
申请号 CN201310652581.2 申请日期 2013.12.05
申请人 国家电网公司;重庆新世杰电气股份有限公司;国网重庆市电力公司璧山供电分公司 发明人 唐昆明;黄翰;冯彦;唐辰旭;胡建;杨伟;康丽红;陈洪波;孙小江
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 重庆为信知识产权代理事务所(普通合伙) 50216 代理人 余锦曦
主权项 一种利用暂态主频分量的配网电缆单相接地故障测距方法,包括有如下步骤:步骤(1)、采样故障电缆线路首端(M)的三相暂态电流信号i<sub>MA</sub>(t)、i<sub>MB</sub>(t)、i<sub>MC</sub>(t)和三相暂态电压信号u<sub>MA</sub>(t)、u<sub>MB</sub>(t)、u<sub>MC</sub>(t);步骤(2)、求取故障电缆线路首端电压的0模分量u<sub>M0</sub>(t)和电流的0模分量i<sub>M0</sub>(t):①、对步骤(1)中,采样得到的故障电缆线路的首端三相暂态电压信号序列值u<sub>MA</sub>(t)、u<sub>MB</sub>(t)、u<sub>MC</sub>(t)进行Karenbuaer相模变换,得到故障电缆线路首端电压的0模分量u<sub>M0</sub>(t),具体的变换矩阵如下:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mrow><mi>M</mi><mn>0</mn></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mi>M</mi><mn>1</mn></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mi>M</mi><mn>2</mn></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mrow><mi>M</mi><mi>A</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mi>M</mi><mi>C</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000964986660000011.GIF" wi="825" he="271" /></maths>②、对步骤(1)中,采样得到的故障电缆线路的首端三相暂态电流信号序列值i<sub>MA</sub>(t)、i<sub>MB</sub>(t)、i<sub>MC</sub>(t)进行Karenbuaer相模变换,得到故障电缆线路首端电流的0模分量i<sub>M0</sub>(t),具体的变换矩阵如下:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>i</mi><mrow><mi>M</mi><mn>0</mn></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>i</mi><mrow><mi>M</mi><mn>1</mn></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>i</mi><mrow><mi>M</mi><mn>2</mn></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>i</mi><mrow><mi>M</mi><mi>A</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>i</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>i</mi><mrow><mi>M</mi><mi>C</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000964986660000012.GIF" wi="796" he="271" /></maths>步骤(3)、采用滤波器对步骤(2)中得到的电压电流0模分量u<sub>M0</sub>(t)、i<sub>M0</sub>(t)进行滤波;步骤(4)、采用Karenbuaer相模变换,将三相系统解耦为0模、1模、2模系统,并对0模网络的时域数学表达式进行拉式变换,得到0模网络的象函数表达式如下:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>U</mi><mi>M</mi></msub><mo>(</mo><mi>s</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>M</mi></msub><mo>(</mo><mi>s</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>h</mi><mi>&gamma;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>l</mi></mrow></mtd><mtd><mrow><msub><mi>Z</mi><mi>C</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>s</mi><mi>h</mi><mi>&gamma;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>l</mi></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mn>1</mn><mrow><msub><mi>Z</mi><mi>C</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mi>s</mi><mi>h</mi><mi>&gamma;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>l</mi></mrow></mtd><mtd><mrow><mi>c</mi><mi>h</mi><mi>&gamma;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>l</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>U</mi><mi>N</mi></msub><mo>(</mo><mi>s</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>N</mi></msub><mo>(</mo><mi>s</mi><mo>)</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000964986660000021.GIF" wi="1355" he="293" /></maths>其中:U<sub>M</sub>(s)、I<sub>M</sub>(s)为线路首端的电压电流象函数、U<sub>N</sub>(s)、I<sub>N</sub>(s)为线路末端的电压电流象函数,l为线路全长,<img file="FDA0000964986660000022.GIF" wi="594" he="77" />为传播系数的运算式,<img file="FDA0000964986660000023.GIF" wi="641" he="77" />为波阻抗的运算式,其中的R<sub>0</sub>、L<sub>0</sub>、C<sub>0</sub>、G<sub>0</sub>分别为线路单位长度的电阻、电感、电容、电导,s为拉式变换中的复数;步骤(5)、利用步骤(3)中获取的暂态信号的0模分量序列,采用Prony算法,求得暂态0模电压电流信号的全频和主频时域表达式,并通过拉式变换得到其象函数表达式,其中全频象函数表达式如下:<img file="FDA0000964986660000024.GIF" wi="1629" he="439" />K为零模电压电流各个分量的个数、a<sub>MK</sub>为第K次分量的衰减时间常数、ω<sub>MK</sub>为第K次分量的频率、U<sub>M0K</sub>为零模电压第K次分量的幅值、I<sub>M0K</sub>为零模电流第K次分量的幅值、<img file="FDA0000964986660000031.GIF" wi="99" he="54" />为零模电压第K次分量的相角、δ<sub>MK</sub>为零模电流第K次分量的相角;主频象函数表达式如下:<img file="FDA0000964986660000032.GIF" wi="1244" he="359" />其中,a<sub>MZ</sub>为主频分量的衰减时间常数、ω<sub>MZ</sub>为主频分量的频率、U<sub>M0Z</sub>为零模电压主频分量的幅值、I<sub>M0Z</sub>为零模电流主频分量的幅值、<img file="FDA0000964986660000033.GIF" wi="94" he="54" />为零模电压主频分量的相角、δ<sub>MZ</sub>为零模电流主频分量的相角;步骤(6)、结合步骤(4)、步骤(5)中故障线路的零模等效网络及零模电压电流的象函数表达式,推导并得到单端主频的象函数测距方程;步骤(7)、将步骤(6)中的象函数测距方程,通过频域法或时域法求解该测距方程。
地址 100033 北京市西城区西长安街86号