CN106019080B - 一种基于沿线能量突变的同塔双回直流线路单端行波测距方法 - Google Patents
一种基于沿线能量突变的同塔双回直流线路单端行波测距方法 Download PDFInfo
- Publication number
- CN106019080B CN106019080B CN201610334759.2A CN201610334759A CN106019080B CN 106019080 B CN106019080 B CN 106019080B CN 201610334759 A CN201610334759 A CN 201610334759A CN 106019080 B CN106019080 B CN 106019080B
- Authority
- CN
- China
- Prior art keywords
- line
- voltage
- traveling wave
- fault
- along
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Locating Faults (AREA)
Abstract
本发明涉及一种基于沿线能量突变的同塔双回直流线路单端行波测距方法,属于电力系统继电保护技术领域。当同塔双回直流线路发生接地故障时,首先,对量测端获取电压求取故障分量,并采用相模变换求取线模电压分量;其次,应用贝杰龙传输方程和量测端线模电压分量计算沿线电压和电流行波变化量;再次求沿线电压行波和沿线电流行波变化量乘积的三次方,并于观测时窗内进行积分来构造测距函数;结合测距函数的沿线突变规律及故障电压行波的小波变换模极大值点极性来实现故障测距。本发明针对同塔双回直流线路进行故障定位,其原理简单,毋需标定故障行波波头,且不受故障瞬时性、故障过渡电阻变化等因素的影响,测距结果准确可靠。
Description
技术领域
本发明涉及一种基于沿线能量突变的同塔双回直流线路单端行波测距方法,属于电力系统继电保护技术领域。
背景技术
直流线路是同塔双回直流输电系统的重要组成部分,其电压等级高,输电距离远,在输送大容量电能的过程中途经区域环境恶劣,极易引发故障。对直流线路的故障点进行精确定位技术有助于快速、准确地确定故障位置,减轻巡线盲目性、加速恢复送电。直流线路故障定位大多采用行波原理。目前的行波测距方法大多是基于故障行波时域特征并于时间轴上对行波进行观测、刻画和波头标定,以及故障距离的计算。其中,时域单端行波测距需要在行波波头标定和波头识别的可靠性、测距分析的自动化方面做进一步研究;时域双端行波测距由于利用故障线路两侧的初始行波波到时差,其初始行波标定可靠性和准确性易得以保证,且毋需对故障点反射波进行辨识,但双端行波测距对线路两端时钟精确同步要求较高。因此,急需提出一种新的故障测距方法,不受故障行波波头的有效辨识和测距时钟精确同步对故障定位准确性的影响。
发明内容
本发明要解决的技术问题是克服传统时域行波测距要求故障行波有效辨识和测距时钟精确同步的局限性,提出一种基于沿线能量突变的同塔双回直流线路单端行波测距方法。
本发明的技术方案是:一种基于沿线能量突变的同塔双回直流线路单端行波测距方法,当同塔双回直流线路发生接地故障时,首先,对量测端获取电压求取故障分量,并采用相模变换求取线模电压分量;其次,应用贝杰龙传输方程和量测端线模电压分量计算沿线电压和电流行波变化量;再次求沿线电压行波和沿线电流行波变化量乘积的三次方,并于观测时窗内进行积分来构造测距函数;结合测距函数的沿线突变规律及故障电压行波的小波变换模极大值点极性来实现故障测距。
具体步骤为:
第一步、当同塔双回直流线路发生接地故障时,在采样率1MHz下,对量测端电压行波信号进行采样,用量测端所获故障极线电压减去故障前的极线电压得到极线电压故障分量:
式(1)中,下标ⅠP、ⅠN、ⅡP和ⅡN分别表示Ⅰ回线正极、Ⅰ回线负极、Ⅱ回线正极和Ⅱ回线负极;下标R表示整流侧量测端M;下标|0|表示故障前电气量;
第二步、对步骤(1)的极线电压突变量作解耦变换,得到线模电压故障分量:
式(2)中,S-1为电压解耦矩阵;u′R,0(t)为零模电压分量;u′R,1(t)、u′R,2(t)和u′R,3(t)为解耦变换得到的三个线模电压分量;
第三步、选取步骤(2)得到的线模分量电压,并利用式(3)和式(4)计算同塔双回直流线路的沿线电压行波和沿线电流行波变化量:
式(3)和式(4)中,下标s表示线模分量,取s=1,2,或3;rs,Zc,s,vs分别为s模行波的电阻、波阻抗和波速;x为距量测端的距离;uR(x,Δt)是t时刻距离R端x处的s模沿线电压行波变化量;iR(x,Δt)是t时刻距离R端x处的s模沿线电流行波变化量;
第四步、采用式(5),将步骤(3)得到的量测端s模沿线电压行波变化量和沿线电流行波变化量相乘并求其三次方,最后于时窗长度l/vs进行积分来构造测距函数;
式(5)中,l为故障极线的长度;t0为故障初始行波到达量测端R的时刻;
根据式(5)可以得到测距函数沿线的突变点解集x=[x1,x2,……xn];
第五步、故障位置的获取:
对突变点解集x=[x1,x2,……xn]中,寻找对应距离之和等于故障线长的两个突变点,即满足判据:
x+x*=l x,x*∈[x1,x2,……xn] (6)
满足式(6)的x和x*均含有故障位置的信息,对整流侧量测端R所获故障电压行波进行小波变换,并求取小波变换模极大值;
根据x对应的故障位置确定第2个故障行波:
若第二个模极大值点的极性与第一个模极大值点的极性相反,则判断故障点位于半线长之内,且故障点为离开量测端M的距离为x处;
若第二个模极大值点的极性与第一个模极大值点的极性相同,则判断故障点位于半线长之外,且故障点为离开量测端M的距离为l-x处。
本发明的有益效果是:针对同塔双回直流线路进行故障定位,其原理简单,毋需标定故障行波波波头,且不受故障瞬时性、故障过渡电阻变化等因素的影响,测距结果准确可靠。
附图说明
图1是本发明实施例1、实施例2的同塔双回直流线路结构图;
图2是本发明实施例1半线长之内故障下量测端M电压行波;
图3是本发明实施例1半线长之内故障下量测端M电压行波模量;
图4是本发明实施例1半线长之内测距函数fM(x)的突变分布结果;
图5是本发明实施例1半线长之内故障电压行波小波变换模极大值结果;
图6是本发明实施例2半线长之外故障下量测端M电压行波行波;
图7是本发明实施例2半线长之外故障下量测端M电压行波模量;
图8是本发明实施例2半线长之外测距函数fM(x)的突变分布结果;
图9是本发明实施例2半线长之外故障电压行波小波变换模极大值结果。
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步说明。
一种基于沿线能量突变的同塔双回直流线路单端行波测距方法,当同塔双回直流线路发生接地故障时,首先,对量测端获取电压求取故障分量,并采用相模变换求取线模电压分量;其次,应用贝杰龙传输方程和量测端线模电压分量计算沿线电压和电流行波变化量;再次求沿线电压行波和沿线电流行波变化量乘积的三次方,并于观测时窗内进行积分来构造测距函数;结合测距函数的沿线突变规律及故障电压行波的小波变换模极大值点极性来实现故障测距。
具体步骤为:
第一步、当同塔双回直流线路发生接地故障时,在采样率1MHz下,对量测端电压行波信号进行采样,用量测端所获故障极线电压减去故障前的极线电压得到极线电压故障分量:
式(1)中,下标ⅠP、ⅠN、ⅡP和ⅡN分别表示Ⅰ回线正极、Ⅰ回线负极、Ⅱ回线正极和Ⅱ回线负极;下标R表示整流侧量测端M;下标|0|表示故障前电气量;
第二步、对步骤(1)的极线电压突变量作解耦变换,得到线模电压故障分量:
式(2)中,S-1为电压解耦矩阵;u′R,0(t)为零模电压分量;u′R,1(t)、u′R,2(t)和u′R,3(t)为解耦变换得到的三个线模电压分量;
第三步、选取步骤(2)得到的线模分量电压,并利用式(3)和式(4)计算同塔双回直流线路的沿线电压行波和沿线电流行波变化量:
式(3)和式(4)中,下标s表示线模分量,取s=1,2,或3;rs,Zc,s,vs分别为s模行波的电阻、波阻抗和波速;x为距量测端的距离;uR(x,Δt)是t时刻距离R端x处的s模沿线电压行波变化量;iR(x,Δt)是t时刻距离R端x处的s模沿线电流行波变化量;
第四步、采用式(5),将步骤(3)得到的量测端s模沿线电压行波变化量和沿线电流行波变化量相乘并求其三次方,最后于时窗长度l/vs进行积分来构造测距函数;
式(5)中,l为故障极线的长度;t0为故障初始行波到达量测端R的时刻;
根据式(5)可以得到测距函数沿线的突变点解集x=[x1,x2,……xn];
第五步、故障位置的获取:
对突变点解集x=[x1,x2,……xn]中,寻找对应距离之和等于故障线长的两个突变点,即满足判据:
x+x*=l x,x*∈[x1,x2,……xn] (6)
满足式(6)的x和x*均含有故障位置的信息,对整流侧量测端R所获故障电压行波进行小波变换,并求取小波变换模极大值;
根据x对应的故障位置确定第2个故障行波:
若第二个模极大值点的极性与第一个模极大值点的极性相反,则判断故障点位于半线长之内,且故障点为离开量测端M的距离为x处;
若第二个模极大值点的极性与第一个模极大值点的极性相同,则判断故障点位于半线长之外,且故障点为离开量测端M的距离为l-x处。
实施例1:采用如图1所示的±500kV同塔双回直流输电系统。线路全长为1286km,采用四分裂导线,线路两侧装有0.3H的平波电抗器,每极单12脉动阀组接线,额定输送功率为6400MW,额定电流为3200A。采样率设为1MHz。
根据步骤一,得到整流侧量测端M的各极线电压故障分量波形,如图2所示;根据步骤二,采用式(2)对量测端M的极线电压突变量进行电压解耦变换,得到量测端M的电压线模分量,如图3所示;根据步骤三,采用式(3)和式(4)计算沿线电压和电流变化量;根据步骤四,采用式(5)构造测距函数fM(x),其沿线能量分布结果如图4所示。由图4可知,测距函数沿线的突变点解集x=[598 683 1195];根据步骤五,得到x1+x2=l,可知突变点x1和x2含有故障位置的信息。对量测端M的故障电压行波进行小波变换模极大值分析,如图5所示,根据x1=598,得到第二个模极大值点的极性与第一个模极大值点的极性相反,故判断故障发生在半线长内,且故障距离开量测端M598km。
实施例2:同样采用如图1所示的±500kV同塔双回直流输电系统。线路全长为1286km,采用四分裂导线,线路两侧装有0.3H的平波电抗器,每极单12脉动阀组接线,额定输送功率为6400MW,额定电流为3200A。采样率设为1MHZ。假设整流侧电气量可测,若Ⅰ回线正极ⅠP的半线长之外距离整流侧量测端M900km处发生金属性接地故障。
根据步骤一,得到整流侧量测端M的各极线电压故障分量波形,如图6所示;据步骤二,采用式(2)对量测端M的极线电压突变量进行电压解耦变换,得到量测端M的电压线模分量,如图7所示;根据步骤三,采用式(3)和式(4)计算沿线电压和电流变化量;根据步骤四,采用式(5)构造测距函数fM(x),其沿线能量分布结果如图8所示。由图8可知,测距函数沿线的突变点解集x=[384 770 889];根据步骤五,得到x1+x3=l,可知突变点x1和x3含有故障位置的信息。对量测端M的故障电压行波进行小波变换模极大值分析,如图9所示,根据x1=384,得到第二个模极大值点的极性与第一个模极大值点的极性相,故判断故障发生在半线长内,且故障距离开量测端M:l-384=902km。
以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (1)
1.一种基于沿线能量突变的同塔双回直流线路单端行波测距方法,其特征在于:当同塔双回直流线路发生接地故障时,首先,对量测端获取电压求取故障分量,并采用相模变换求取线模电压分量;其次,应用贝杰龙传输方程和量测端线模电压分量计算沿线电压和电流行波变化量;再次求沿线电压行波和沿线电流行波变化量乘积的三次方,并于观测时窗内进行积分来构造测距函数;结合测距函数的沿线突变规律及故障电压行波的小波变换模极大值点极性来实现故障测距;
具体步骤为:
第一步、当同塔双回直流线路发生接地故障时,在采样率1MHz下,对量测端电压行波信号进行采样,用量测端所获故障极线电压减去故障前的极线电压得到极线电压故障分量:
式(1)中,下标ⅠP、ⅠN、ⅡP和ⅡN分别表示Ⅰ回线正极、Ⅰ回线负极、Ⅱ回线正极和Ⅱ回线负极;下标R表示整流侧量测端M;下标|0|表示故障前电气量;
第二步、对步骤(1)的极线电压突变量作解耦变换,得到线模电压故障分量:
式(2)中,S-1为电压解耦矩阵;u′R,0(t)为零模电压分量;u′R,1(t)、u′R,2(t)和u′R,3(t)为解耦变换得到的三个线模电压分量;
第三步、选取步骤(2)得到的线模分量电压,并利用式(3)和式(4)计算同塔双回直流线路的沿线电压行波和沿线电流行波变化量:
式(3)和式(4)中,下标s表示线模分量,取s=1,2,或3;rs,Zc,s,vs分别为s模行波的电阻、波阻抗和波速;x为距量测端的距离;uR(x,Δt)是t时刻距离R端x处的s模沿线电压行波变化量;iR(x,Δt)是t时刻距离R端x处的s模沿线电流行波变化量;
第四步、采用式(5),将步骤(3)得到的量测端s模沿线电压行波变化量和沿线电流行波变化量相乘并求其三次方,最后于时窗长度l/vs进行积分来构造测距函数;
式(5)中,l为故障极线的长度;t0为故障初始行波到达量测端R的时刻;
根据式(5)可以得到测距函数沿线的突变点解集x=[x1,x2,……xn];
第五步、故障位置的获取:
对突变点解集x=[x1,x2,……xn]中,寻找对应距离之和等于故障线长的两个突变点,即满足判据:
x+x*=l x,x*∈[x1,x2,……xn] (6)
满足式(6)的x和x*均含有故障位置的信息,对整流侧量测端R所获故障电压行波进行小波变换,并求取小波变换模极大值;
根据x对应的故障位置确定第2个故障行波:
若第二个模极大值点的极性与第一个模极大值点的极性相反,则判断故障点位于半线长之内,且故障点为离开量测端M的距离为x处;
若第二个模极大值点的极性与第一个模极大值点的极性相同,则判断故障点位于半线长之外,且故障点为离开量测端M的距离为l-x处。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610334759.2A CN106019080B (zh) | 2016-05-19 | 2016-05-19 | 一种基于沿线能量突变的同塔双回直流线路单端行波测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610334759.2A CN106019080B (zh) | 2016-05-19 | 2016-05-19 | 一种基于沿线能量突变的同塔双回直流线路单端行波测距方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019080A CN106019080A (zh) | 2016-10-12 |
CN106019080B true CN106019080B (zh) | 2019-04-12 |
Family
ID=57095493
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610334759.2A Active CN106019080B (zh) | 2016-05-19 | 2016-05-19 | 一种基于沿线能量突变的同塔双回直流线路单端行波测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019080B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108808635B (zh) * | 2018-04-19 | 2020-06-19 | 华北电力大学(保定) | 一种高压直流输电线路单极故障隔离方法 |
CN109119977B (zh) * | 2018-09-20 | 2019-07-23 | 山东大学 | 基于单端电压的多端柔性直流电网直流线路快速保护方法及系统 |
CN113945806B (zh) * | 2021-11-15 | 2024-04-09 | 昆明理工大学 | 一种柔性直流输电线路单端故障测距方法及系统 |
CN115015686A (zh) * | 2022-01-26 | 2022-09-06 | 昆明理工大学 | 一种lcc-vsc混合高压直流输电线路故障测距方法及系统 |
CN115902530A (zh) * | 2023-03-10 | 2023-04-04 | 昆明理工大学 | 一种接地极线路故障测距方法及系统 |
CN115963358B (zh) * | 2023-03-13 | 2023-08-04 | 昆明理工大学 | 一种混合三端柔性直流输电线路故障测距方法及系统 |
CN117192292B (zh) * | 2023-11-07 | 2024-02-06 | 昆明理工大学 | 一种雷击接地极线路故障测距方法及系统 |
CN117805555B (zh) * | 2024-02-28 | 2024-05-03 | 昆明理工大学 | 双向推演的四端柔性直流输电线路故障测距方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590694A (zh) * | 2012-02-21 | 2012-07-18 | 昆明理工大学 | 一种基于集中参数t模型的同杆双回输电线路区内外故障测后模拟识别方法 |
CN103412240A (zh) * | 2013-07-24 | 2013-11-27 | 昆明理工大学 | 一种不依赖波头辨识的同塔双回输电线路单端行波故障测距方法 |
CN103954885A (zh) * | 2014-05-20 | 2014-07-30 | 华北电力大学 | 基于分布参数的双回线故障单端定位系统及定位方法 |
CN104655988A (zh) * | 2015-02-28 | 2015-05-27 | 华南理工大学 | 一种基于单回信息的同塔双回直流线路行波故障定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030085715A1 (en) * | 2001-08-15 | 2003-05-08 | David Lubkeman | System and method for locating a fault on ungrounded and high-impedance grounded power systems |
-
2016
- 2016-05-19 CN CN201610334759.2A patent/CN106019080B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590694A (zh) * | 2012-02-21 | 2012-07-18 | 昆明理工大学 | 一种基于集中参数t模型的同杆双回输电线路区内外故障测后模拟识别方法 |
CN103412240A (zh) * | 2013-07-24 | 2013-11-27 | 昆明理工大学 | 一种不依赖波头辨识的同塔双回输电线路单端行波故障测距方法 |
CN103954885A (zh) * | 2014-05-20 | 2014-07-30 | 华北电力大学 | 基于分布参数的双回线故障单端定位系统及定位方法 |
CN104655988A (zh) * | 2015-02-28 | 2015-05-27 | 华南理工大学 | 一种基于单回信息的同塔双回直流线路行波故障定位方法 |
Non-Patent Citations (1)
Title |
---|
耦合双回线路电弧故障测距的新相模变换方法;束洪春 等;《高电压技术》;20090331;第35卷(第3期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106019080A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106019080B (zh) | 一种基于沿线能量突变的同塔双回直流线路单端行波测距方法 | |
CN106019079B (zh) | 一种同塔双回直流线路新型双端故障测距方法 | |
Suonan et al. | A novel fault-location method for HVDC transmission lines | |
CN104898021B (zh) | 一种基于k‑means聚类分析的配电网故障选线方法 | |
CN106093708B (zh) | 一种波速修正的同塔双回直流输电线路双端故障测距方法 | |
CN102135571B (zh) | 超高压/特高压多回路输电线路零序阻抗抗干扰测量方法 | |
CN101762774B (zh) | 基于遗传算法参数辨识高压直流输电线路故障测距方法 | |
CN109541392A (zh) | 一种适用于柔性直流输电系统的单端故障测距方法 | |
CN105353268A (zh) | 一种用于输电线路分布式行波故障判断及定位方法 | |
CN107359603B (zh) | 基于单端量的高压直流线路多段式距离保护方法及系统 | |
CN105548802B (zh) | 一种基于故障行波沿线分布特性的t接线路三端不同步故障测距方法 | |
CN103630814A (zh) | 高压电缆在交叉互联下绝缘介质损耗角趋势在线监测方法 | |
CN103399209A (zh) | 一种特高压双极直流输电线路工频参数测量方法 | |
CN101923137A (zh) | 小电流接地系统单相接地故障选线方法 | |
CN106501612B (zh) | 基于非均匀传输线的双回部分同塔线路零序参数测量方法 | |
CN102590704A (zh) | 一种基于贝杰龙模型的同杆双回输电线路区内外故障识别的测后模拟方法 | |
CN103941153A (zh) | 一种基于波形相似性的k-NN算法的多出线辐射网故障测距方法 | |
CN105842582B (zh) | 基于emtr的柔性直流线路故障测距方法 | |
CN106093705A (zh) | 一种配电网单相接地故障行波波头的计算方法 | |
CN106443540B (zh) | 一种基于仿真数据的行波单端测距装置测试方法 | |
CN103605003B (zh) | 超/特高压同塔四回输电线路零序电容高精度测量方法 | |
CN107179473B (zh) | 一种输电线路故障定位方法 | |
CN108845233A (zh) | 配电网架空线路单相接地双端检测定位方法 | |
CN107121617B (zh) | 一种采用滤波器支路电流的直流输电线路k近邻测距方法 | |
CN105652156A (zh) | 特高压交流输电线路单相接地电压相位突变测距方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |