CN106646121B - 一种配电网故障行波波头的辨识方法 - Google Patents
一种配电网故障行波波头的辨识方法 Download PDFInfo
- Publication number
- CN106646121B CN106646121B CN201611070835.XA CN201611070835A CN106646121B CN 106646121 B CN106646121 B CN 106646121B CN 201611070835 A CN201611070835 A CN 201611070835A CN 106646121 B CN106646121 B CN 106646121B
- Authority
- CN
- China
- Prior art keywords
- state
- signal
- fault
- traveling wave
- moment
- 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
- 238000012850 discrimination method Methods 0.000 title claims abstract description 10
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 8
- 230000001052 transient effect Effects 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract 1
- 238000013178 mathematical model Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
-
- 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/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
本发明公开了一种配电网故障行波波头的辨识方法。对故障电压行波进行凯伦贝尔变换得到故障电压线模分量,接着构建行波波头数学模型,根据该模型构造行波信号状态方程和观测方程,然后运用无迹卡尔曼滤波对故障电压行波去噪,并根据收敛的状态变量得到故障行波到达时刻。本发明能够在信噪比低的情况下依旧准确辨识行波波头,具有很好的工程应用前景。
Description
技术领域
本发明属于电力系统自动化技术领域,涉及配电网的故障诊断,是一种精确的配电网故障行波波头的辨识方法。
背景技术
随着我国国民经济的发展,人们对供电可靠性的要求越高,与居民用电密切相关的配电网运行可靠性便越高。这便要求对配电网的故障能够准确识别,迅速排除。然而,配电网的故障定位问题长期以来一直没有得到很好的解决。因为配电网一般采用架空线路—电缆线路的混合线路,且线路分支较多,结构复杂,使配电网在故障定位上要比输电网困难,许多在输电网中已经成熟的技术在配电网中无法实现。在发生故障后,由于故障电流比较微弱,对故障的定位更加困难。
行波测距技术最早应用在输电线路上,其测量的准确度高,反应迅速,并且不受线路结构不对称、互感器传变误差等因素的影响,因此得到了广泛的应用。在行波测距中,最为关键的技术是对故障发生后的暂态行波波头进行准确的识别。然而,由于暂态行波的频带在10k~100kHz不等,要完整采集故障行波,由采样定理可知,保护装置的采样率就必须非常高。高采样率会受到各种噪声的干扰,例如白噪声,脉冲噪声等,这就使得行波波头的识别非常困难。再加上在配电网中,线路电阻较大,行波能量衰减非常严重,使得行波的波头变缓,奇异性降低,在噪声的淹没下,波头的辨识就非常困难。极端情况下,系统在过零点或者经高阻接地的情况下,行波信号非常微弱,与噪声相比奇异性并未差异,利用传统的小波变换等方法难以获得准确的故障行波波头到达时刻。
无迹卡尔曼滤波是基于最小方差原理的滤波方法,能够在时域内滤除环境噪声的干扰,减少伪故障点的干扰。通过构建合适的状态方程,利用卡尔曼滤波在状态估计上的能力,可以对故障行波进行辨识,推算出行波到达检测装置的时刻。
发明内容
本发明对现有行波定位技术中的不足进行改进,提出一种配电网故障行波波头的辨识方法。
本发明采用以下技术方案:
一种配电网故障行波波头的辨识方法,其特征在于:所述方法包括以下步骤:
步骤1:对故障采集器得到的三相故障电压行波进行凯伦贝尔变换,得到线模电压分量,其中,凯伦贝尔变换为:式中u1,u2为线模电压,u0为零模电压,ua,ub,uc为各相电压;
步骤2:构建行波波头数学表达式,其表达式可以由以下函数表征:
其中,u(k)为k时刻故障线模电压信号的值,T是行波到达故障检测装置的时刻,当k<T时,故障行波还没有到达,系统正常,As表示基波信号的幅值,ωc为基波信号的角频率,为基波信号的初始相角;当k>T时,信号发生突变,在基波信号的基础上,叠加了暂态行波,其中Ae表示初始暂态行波波头的幅值,Ts为行波衰减时间常数。
步骤3:根据步骤2构建的行波信号表达式,选取状态变量构建离散化的行波信号状态方程和观测方程如下:
x(k)=x(k-1) k=1,2...
u(k)=H[x(k-1),k]+V(k) k=1,2...
式中,x(K)为k时刻系统状态变量,x(k‐1)为k‐1时刻系统状态变量,H[]表示故障电压信号u(k)与状态变量的函数关系,V(K)表示观测噪声。
步骤4:运用无迹卡尔曼滤波对故障行波电压进行处理,得到去噪的行波和状态变量T的收敛值。
4‐1、初始化状态变量均值:
其中,x0为初始值,为初始向量均值,P0是初始协方差矩阵。
4‐2、计算sigma采样点:
其中L为状态空间的维数,为6,则sigma点的个数为2L+1=13;为第k‐1时刻的状态变量均值;Xi(k-1)为第k‐1时刻的sigma点集;Sigma点对应的权重系数为:
其中,λ为比例系数,改变λ的大小可以调节sigma点与均值点之间的距离。
4‐3、计算状态变量均值的一步估计,将k‐1时刻的sigma点分别代入状态方程,得到k时刻的sigma点集以及其对应的均值
4‐4、计算协方差矩阵的一步估计:
4‐5、获得卡尔曼增益:
其中,是观测点集,是观测点集的均值,PU(k),U(k)是观测向量的自相关矩阵,PX(k),U(k)是观测向量与状态向量的互相关矩阵,Kk为卡尔曼增益。
4‐6、更新状态估计值和协方差矩阵:
式中u(k)为步骤1变换得到的故障行波线模分量在k时刻的值。
完成k时刻的状态变量和协方差矩阵的更新后,返回步骤4‐2,进行k+1时刻的估算。当k=N(N为故障信号的采样点数),状态变量和协方差矩阵停止更新,此时输出变换得到的电压信号U(k)和状态变量T的值,此时的U(k)为去噪后的行波信号,T为故障行波波头到达检测装置的时刻。
所述步骤4‐1协方差矩阵P0的取值如下:
其中0.0001<P11,P22...P66<0.1。
本发明与现有技术相比,具有以下优点:
本发明采用无迹卡尔曼滤波进行波头辨识,能够滤除行波传播过程中遇到的干扰信号,减少伪故障点的存在,精度更高。
本发明提出的行波波头辨识的方法,能够在行波信号微弱的情况下进行辨识,克服了原有算法在波头奇异性降低的情况下无法识别的问题,具有更强的适用性。
附图说明
图1为实施配电网定位实验仿真结构图;
图2为本发明公开的配电网故障行波波头辨识方法流程图;
图3为故障行波线模电压分量;
图4为经过无迹卡尔曼滤波后的行波线模电压分量。
具体实施方式
下面结合说明书附图和具体实施例对本发明的技术方案作进一步详细介绍。
如图1所示,本申请列举的实施例是配电网C相接地故障测距实验,采用的是ATP‐EMTP搭建线路模型进行实验。其中,电压源设计为理想电源,初始角为0°,来模拟三相无穷大电源。变压器连结组别号为YD11,二次侧电压为10.5KV。变压器后为简单线路T型线路。线路末端故障采集装置采样频率为1MHz,即每个采样点之间的时间间隔为1μs,故障发生在第2000μs,故障采集装置距离线路故障点19km。由于行波信号线模分量的传播速度基本不受频率、传播距离的影响,设定其传播速度为3×108m/s,则故障发生后行波信号线模分量传播到检测装置的时间为:
下面以附图1为实施例,详细介绍配电网故障行波波头的辨识方法,方法的步骤如附图2所示:
步骤(1),对故障采集器得到的三相故障电压行波进行凯伦贝尔变换,得到线模电压分量,由于系统发生C相接地故障,线模电压分量为ua,uc分别为A相,C相电压,如图3所示。
步骤(2),构建行波波头数学表达式,其表达式可以由以下函数表征:
其中,u(k)为k时刻故障线模电压信号的值,T是行波到达故障检测装置的时刻,当k<T时,故障行波还没有到达,系统正常,As表示基波信号的幅值,ωc为基波信号的角频率,为基波信号的初始相角;当k>T时,信号发生突变,在基波信号的基础上,叠加了暂态行波,其中Ae表示初始暂态行波波头的幅值,Ts为行波衰减时间常数。
步骤(3),根据步骤2构建的行波信号表达式,选取状态变量构建离散化的行波信号状态方程和观测方程如下:
x(k)=x(k-1) k=1,2...
u(k)=H[x(k-1),k]+V(k) k=1,2...
式中,x(K)为k时刻系统状态变量,x(k‐1)为k‐1时刻系统状态变量,H[]表示故障电压信号u(k)与状态变量的函数关系,V(K)表示观测噪声。
步骤(4):运用无迹卡尔曼滤波对故障行波电压进行处理,得到去噪的行波和状态变量T的收敛值。
4‐1、初始化状态变量均值:
其中,x0为初始值,根据仿真的10kV线路的实际情况,选取初始值为x0=(15000,314,0,10000,1800,500),
4‐2、计算sigma采样点:
其中L为状态空间的维数,为6,则sigma点的个数为2L+1=13;为第k‐1时刻的状态变量均值;Xi(k-1)为第k‐1时刻的sigma点集;Sigma点对应的权重系数为:
其中,选择λ=0.5来进行计算
4‐3、计算状态变量均值的一步估计,将k‐1时刻的sigma点分别代入状态方程,得到k时刻的sigma点集以及其对应的均值
4‐4、计算协方差矩阵的一步估计:
4‐5、获得卡尔曼增益:
其中,是观测点集,是观测点集的均值,PU(k),U(k)是观测向量的自相关矩阵,PX(k),U(k)是观测向量与状态向量的互相关矩阵,Kk为卡尔曼增益。
4‐6、更新状态估计值和协方差矩阵:
式中u(k)为步骤1变换得到的故障行波线模分量在k时刻的值。
完成k时刻的状态变量和协方差矩阵的更新后,返回步骤4‐2,进行k+1时刻的估算。当k=N(N为故障信号的采样点数,本实例中N=2500),状态变量和协方差矩阵停止更新,此时输出变换得到的去噪后的电压信号U(k)如图4所示,此时T=2062.134μs,与理论值T=2063μs相差不超过1μs,说明本发明具有很高的工程实用性。
以上给出的实施例用以说明本发明和它的实际应用,并非对本发明作任何形式上的限制,任何一个本专业的技术人员在不偏离本发明技术方案的范围内,依据以上技术和方法作一定的修饰和变更当视为等同变化的等效实施例。
Claims (2)
1.一种配电网故障行波波头的辨识方法,其特征在于:所述方法包括以下步骤:
步骤1:对故障采集器得到的三相故障电压行波进行凯伦贝尔变换,得到线模电压分量,其中,凯伦贝尔变换为:式中u1,u2为线模电压,u0为零模电压,ua,ub,uc为各相电压;
步骤2:构建行波波头数学表达式,其表达式可以由以下函数表征:
其中,u(k)为k时刻故障线模电压信号的值,T是行波到达故障检测装置的时刻,当k<T时,故障行波还没有到达,系统正常,As表示基波信号的幅值,ωc为基波信号的角频率,为基波信号的初始相角;当k>T时,信号发生突变,在基波信号的基础上,叠加了暂态行波,其中Ae表示初始暂态行波波头的幅值,Ts为行波衰减时间常数;
步骤3:根据步骤2构建的行波信号表达式,选取状态变量构建离散化的行波信号状态方程和观测方程如下:
x(k)=x(k-1) k=1,2...
u(k)=H[x(k-1),k]+V(k) k=1,2...
式中,x(K)为k时刻系统状态变量,x(k-1)为k-1时刻系统状态变量,H[]表示故障电压信号u(k)与状态变量的函数关系,V(K)表示观测噪声;
步骤4:运用无迹卡尔曼滤波对故障行波电压进行处理,得到去噪的行波和状态变量T的收敛值;
4-1、初始化状态变量均值:
其中,x0为初始值,为初始向量均值,P0是初始协方差矩阵;
4-2、计算sigma采样点:
其中L为状态空间的维数,为6,则sigma点的个数为2L+1=13;为第k-1时刻的状态变量均值;Xi(k-1)为第k-1时刻的sigma点集;Sigma点对应的权重系数为:
其中,λ为比例系数,改变λ的大小可以调节sigma点与均值点之间的距离;
4-3、计算状态变量均值的一步估计,将k-1时刻的sigma点分别代入状态方程,得到k时刻的sigma点集以及其对应的均值
4-4、计算协方差矩阵的一步估计:
4-5、获得卡尔曼增益:
其中,是观测点集,是观测点集的均值,PU(k),U(k)是观测向量的自相关矩阵,PX(k),U(k)是观测向量与状态向量的互相关矩阵,Kk为卡尔曼增益;
4-6、更新状态估计值和协方差矩阵:
式中u(k)为步骤1变换得到的故障行波线模分量在k时刻的值;
完成k时刻的状态变量和协方差矩阵的更新后,返回步骤4-2,进行k+1时刻的估算;当k=N,N为故障信号的采样点数,状态变量和协方差矩阵停止更新,此时输出变换得到的电压信号U(k)和状态变量T的值,此时的U(k)为去噪后的行波信号,T为故障行波波头到达检测装置的时刻。
2.根据权利要求1所述的一种配电网故障行波波头的辨识方法,其特征在于:所述步骤4-1协方差矩阵P0的取值如下:
其中0.0001<P11,P22...P66<0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611070835.XA CN106646121B (zh) | 2016-11-29 | 2016-11-29 | 一种配电网故障行波波头的辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611070835.XA CN106646121B (zh) | 2016-11-29 | 2016-11-29 | 一种配电网故障行波波头的辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646121A CN106646121A (zh) | 2017-05-10 |
CN106646121B true CN106646121B (zh) | 2019-01-22 |
Family
ID=58813973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611070835.XA Active CN106646121B (zh) | 2016-11-29 | 2016-11-29 | 一种配电网故障行波波头的辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646121B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103160A (zh) * | 2017-05-25 | 2017-08-29 | 长沙理工大学 | 基于贝叶斯滤波的微弱故障行波信号去噪和精确识别方法 |
CN108051705A (zh) * | 2017-12-26 | 2018-05-18 | 华北电力大学 | 配电网的单相接地故障定位方法和系统 |
CN109387733A (zh) * | 2018-12-21 | 2019-02-26 | 云南电网有限责任公司电力科学研究院 | 一种配电线路单相接地故障定位方法及系统 |
CN110703033B (zh) * | 2019-09-17 | 2021-10-15 | 国电南瑞科技股份有限公司 | 一种弱故障行波信号增强方法 |
CN112782532A (zh) * | 2020-10-19 | 2021-05-11 | 国网辽宁省电力有限公司 | 基于断路器合闸产生行波信号的配电网故障测距方法 |
CN112748306B (zh) * | 2020-12-09 | 2024-01-09 | 国网上海市电力公司 | 一种基于卡尔曼滤波对电缆故障进行定位的方法及系统 |
CN113884818B (zh) * | 2021-10-27 | 2024-04-09 | 国网江苏省电力有限公司徐州供电分公司 | 一种基于lstm的配电网故障行波到达时间精确估算方法 |
CN115965080A (zh) * | 2022-11-07 | 2023-04-14 | 河海大学 | 新能源发电单元运行状态识别方法、装置及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006022948A1 (de) * | 2006-05-17 | 2007-11-22 | Deuta-Werke Gmbh | Verfahren und Vorrichtung zur kostengünstigen Ermittlung der Geschwindigkeit aus einem Mikrowellen-Dopplersignal |
CN101291055A (zh) * | 2008-06-18 | 2008-10-22 | 昆明理工大学 | 一种输电线路故障行波初始波头到达时刻精确标定方法 |
CN102288869A (zh) * | 2011-05-10 | 2011-12-21 | 山东大学 | 一种输电线路单端行波故障测距方法 |
CN104965155A (zh) * | 2015-06-26 | 2015-10-07 | 国家电网公司 | 一种基于滤波特性的波头查找方法 |
CN106093705A (zh) * | 2016-06-15 | 2016-11-09 | 华北电力大学 | 一种配电网单相接地故障行波波头的计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10160746A (ja) * | 1996-11-29 | 1998-06-19 | Yokogawa Denshi Kiki Kk | 潮流計測装置 |
-
2016
- 2016-11-29 CN CN201611070835.XA patent/CN106646121B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006022948A1 (de) * | 2006-05-17 | 2007-11-22 | Deuta-Werke Gmbh | Verfahren und Vorrichtung zur kostengünstigen Ermittlung der Geschwindigkeit aus einem Mikrowellen-Dopplersignal |
CN101291055A (zh) * | 2008-06-18 | 2008-10-22 | 昆明理工大学 | 一种输电线路故障行波初始波头到达时刻精确标定方法 |
CN102288869A (zh) * | 2011-05-10 | 2011-12-21 | 山东大学 | 一种输电线路单端行波故障测距方法 |
CN104965155A (zh) * | 2015-06-26 | 2015-10-07 | 国家电网公司 | 一种基于滤波特性的波头查找方法 |
CN106093705A (zh) * | 2016-06-15 | 2016-11-09 | 华北电力大学 | 一种配电网单相接地故障行波波头的计算方法 |
Non-Patent Citations (2)
Title |
---|
基于卡尔曼滤波的行波波头检测算法研究;李传健;《电网与清洁能源》;20160125;第32卷(第1期);全文 |
无损卡尔曼滤波算法在目标跟踪中的应用;蔡琳;《自动化与仪器仪表》;20150725;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106646121A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646121B (zh) | 一种配电网故障行波波头的辨识方法 | |
CN107271768B (zh) | 一种最小二乘拟合动态频率测量方法 | |
He et al. | Natural frequency-based line fault location in HVDC lines | |
CN107621591B (zh) | 一种基于零模行波波速变化特性的输电线路迭代测距方法 | |
CN102944818B (zh) | 一种gis终端电缆故障在线测距方法 | |
CN101957443B (zh) | 声源定位方法 | |
CN103308804B (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN105699855B (zh) | 基于不受行波波速影响的单端行波故障测距计算方法和测距方法 | |
CN109001594A (zh) | 一种故障行波定位方法 | |
CN106093705B (zh) | 一种配电网单相接地故障行波波头的计算方法 | |
CN111222743B (zh) | 一种光纤感知事件的垂直偏移距离及威胁级别判断方法 | |
CN110285909B (zh) | 基于同步压缩变换的索承桥梁瞬时索力计算方法 | |
CN105223466B (zh) | 一种利用模极大值比的特高压直流输电线路单端测距方法 | |
CN104833851B (zh) | 基于分布式相关卡尔曼滤波的电力系统谐波估计方法 | |
CN105738760A (zh) | 一种频域法与行波法结合的高阻故障测距方法 | |
CN106228107B (zh) | 一种基于独立成分分析的超声导波断轨监测方法 | |
CN101982953B (zh) | 宽带无线通信信道频域多维参数化模型及建模方法 | |
CN104880648A (zh) | 架空线配电网单相接地故障的测距方法 | |
CN102967779A (zh) | 一种输电线路分布参数的辨识方法 | |
CN108181059A (zh) | 基于小波信号的多相流管道泄漏声波信号识别方法 | |
Wintoft | Study of the solar wind coupling to the time difference horizontal geomagnetic field | |
CN107179473B (zh) | 一种输电线路故障定位方法 | |
CN109324262A (zh) | 一种基于tt变换和波速优化的输电线路故障测距方法 | |
CN108845233A (zh) | 配电网架空线路单相接地双端检测定位方法 | |
CN105676169A (zh) | 一种脉冲声目标的定向方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |