一种配电网单相接地故障行波波头的计算方法
【专利摘要】本发明公开了一种配电网单相接地故障行波波头的计算方法。线路数据采集器将故障前后各2个电压周波上传，对电压行波进行凯伦贝尔变换得到线模电压进行Hilbert变换，求导，并求取模值。对得到的模值去噪，取有效区间后得到行波波头到达时刻。本发明具有快速，准确的优点，能够准确计算出行波波头到达时刻，为确定故障位置创造条件，具有很好的应用前景。
【专利说明】
一种配电网单相接地故障行波波头的计算方法
技术领域
[0001] 本发明属于电力系统自动化技术领域，涉及配电网的故障诊断，是一种快速精确 的配电网单相接地故障行波波头的计算方法。
【背景技术】
[0002] 随着我国国民经济的发展，人们对供电可靠性的要求越高，与居民用电密切相关 的配电网运行可靠性便越高。这便要求对配电网的故障能够准确识别，迅速排除。然而，配 电网的单相接地故障定位问题长期以来一直没有得到很好的解决。因为配电网一般采用架 空线路一电缆线路的混合线路，且线路分支较多，结构复杂，使配电网在故障定位上要比输 电网困难，许多在输电网中已经成熟的技术在配电网中无法实现。在发生单相接地故障后， 由于故障电流比较微弱，对单相接地故障的定位更加困难。
[0003] 目前用于故障定位的阻抗法，s信号注入法都存在精度不高，自动化水平不足的问 题，而行波法虽然能解决以上问题，但是行波波头难以标定，使得行波法在应用上并没有得 到推广。在计算行波波头的方法中，小波法难以确定小波基函数，而Hilbert-huang法在高 精度的配电网定位中，由于数据量庞大，速度很慢，无法实用化。计算行波波头算法上的不 足，限制了单相接地故障行波定位技术的发展。
[0004] 配电网单相接地故障发生前以及发生后，系统线电压都为标准的正弦波，只有在 故障发生瞬间有一个短暂的行波过程，符合Hilbert变换的前提条件，可以对线电压直接进 行Hilbert变换，计算出波形奇异点，从而计算行波波头到达时刻。但是由于噪声的影响，直 接应用Hilbert变换容易计算出多个伪行波波头，必须进行有效的去噪处理。
【发明内容】
[0005] 本发明的目的在于对现有配电网行波定位的技术进行改进，提出一种配电网单相 接地故障行波波头的计算方法。
[0006] 本发明具体采用以下技术方案：
[0007] -种配电网单相接地故障行波波头的计算方法，其特征在于，所述计算方法包含 以下步骤：
[0008] 步骤1:在线路上安装故障数据采集器，采集故障发生后前后各两个周波的三相电 压数据；
[0009] 步骤2:对步骤1所得到的共四个周波的三相电压行波进行凯伦贝尔变换，得到线
模电压分量，其中，凯伦贝尔变换为： ，式中Ul，U2为线模电压，U〇 为零模电压，Ua，Ub，Uc为各相电压；
[0010]步骤3:选取步骤2得到的线模电压分量进行Hilbert变换，若配电网发生A相或B相 接地故障则选择线模电压m，若配电网发生C相接地故障则选择线模电压u2;
[0011] 步骤4:对步骤3得到的变换后的结果进行中心差分求导，然后取绝对值得到模值 序列Uh，其中，中心差分求导公式为
、,h(n+l)和h(n-l)分别为经过 Hilbert变换得到的线模电压数据在时刻n+1和时刻n-1的值，h'（η)为经过求导后在时刻η 的值；
[0012] 步骤5:对模值序列Uh中的第一个周波数据Xh进行判断，如果Xh(n)>0.5max{Uh}，则 将时刻η的模值视为脉冲干扰置零，得到去除干扰的第一个周波模值序列，其中，Xh(n)为X h 在时刻η的值；
[0013] 步骤6:将Uh中低于阈值Κ的模值视为噪声置零，其中Κ= 1.5max{Xh};
[0014] 步骤7:在步骤6置零后的结果中选取有效区间，区间内模最大值为行波波头到达 时刻。
[0015] 所述步骤(7)中的有效区间选取，以第一个模值大于0.5max {Uh}的点为起点，如果 在0.01个周波时长内至少有2个模值大于0.5max{Uh}，则标定为有效区间，否则，以下一个 模值大于〇. 5max{Uh}的点为起点继续判断直至选出有效区间。
[0016] 本发明与现有技术相比，具有以下优点：
[0017] 1、本发明采用Hilbert算法进行波头提取，加入了阈值去噪和有效空间内取值，能 够减少环境噪声以及小干扰的影响，定位准确。
[0018] 2、本发明采用算法简单易操作，与传统小波法和Hilbert-huang法提取波头相比， 运行时间更短，能够在故障发生后迅速定位出结果。
[0019] 3、本发明克服了传统emd分解所带来的模态混叠现象，以及数据拟合不准确带来 的边界效应，定位结果更具有适应性。
【附图说明】
[0020] 图1为实施配电网定位实验仿真结构图；
[0021] 图2为故障后线路数据采集器采集到的三相电压数据；
[0022] 图3为根据三相电压行波信号得到的线模分量；
[0023] 图4为经过变换后取模值的频率图；
[0024] 图5为经过阈值去噪后得到的频率图；
[0025] 图6为选取有效区间后的频率图；
[0026] 图7为本发明公开的配电网单相接地故障行波波头的计算方法流程示意图。
【具体实施方式】
[0027]下面结合说明书附图和具体实施例对本发明的技术方案作进一步详细介绍。
[0028]如图1所示，本申请列举的实施例是配电网单相接地故障测距实验，采用的是ATP-EMTP搭建线路模型进行实验。其中，电压源设计为理想电源，初始角为0°，来模拟三相无穷 大电源。变压器连结组别号为YD11，二次侧电压为10.5KV。变压器后为简单线路MN。单位长 度线路参数为：
[0029]正序参数：电阻心=0.21 Ω/km，电感Li = 1.2732mH/km，对地电容Ci = 0.00907yF/ km〇
[0030] 零序参数：电阻R〇 = 0.21 Ω/km，电感L〇 = 7.0028mh/km，对地电容C〇 = 0.00466yF/ km〇
[0031] 线路长度:MN长10km。
[0032] 下面以附图1为实施例，详细介绍配电网单相接地故障行波波头的计算方法，计算 方法的步骤如附图7所示：
[0033]步骤（1)，数据采集器采样率为1MHz，即每个点之间的间隔为lys。共采样4个周波 的三相电压数据一共80000ys，在第40000ys时系统发生C相接地故障，接地电阻为1 Ω，故障 一直持续1 s。采集到的电压数据如图2所示
[0034]步骤(2)，对步骤（1)所得到的各相电压行波进行凯伦贝尔变换，得到两个线电压 U1 矛口 U2,
[0035] 步骤(3)，由系统发生C相接地故障，所以选取线模电压分量1!2进行计算，如图3所 示为故障后行波线模分量。对步骤(2)得到的线模分量1! 2进行Hilbert变换。
[0036] 步骤(4)，然后按照公式
进行中心差分求导并求取模值得 到Uh，如图4所示。
[0037] 步骤(5)，对Uh的第一个周波数据，也就是区间（0,20000ys)的数据进行去脉冲干 扰，0 · 5max{Uh} = 327 · 6，将区间(0，20000ys)中大于 327 · 6 的值置零。
[0038] 步骤(6)，根据去除干扰后的区间（0，20000ys)选取阈值，K = 1 · 5max {Xh} = 111 · 3， 将全数据段Uh内低于阈值的模值视为噪声置零，得到去噪后模值序列如图5所示。
[0039]步骤（7)，由图5可得，区间（0，40002ys)全小于0 · 5max{Uh} = 327 · 6，而以第40003μ s为有效区间起始点，后面0.01个周波内模值有至少2个模值大于327.6，所以选取有效区间 为(40003ys，40203ys)，如图6所示，其中区间内最大值为第40004ys，标定线模分量波头到 达时刻为40004ys。
[0040]以上给出的实施例用以说明本发明和它的实际应用，并非对本发明作任何形式上 的限制，任何一个本专业的技术人员在不偏离本发明技术方案的范围内，依据以上技术和 方法作一定的修饰和变更当视为等同变化的等效实施例。
【主权项】
1. 一种配电网单相接地故障行波波头的计算方法，其特征在于:对故障电压行波进行 凯伦贝尔变换得到线模电压再进行Hubert变换，求导，并求取模值，对得到的模值去噪，取 有效区间后得到行波波头到达时刻。2. -种配电网单相接地故障行波波头的计算方法，其特征在于，所述计算方法包含W 下步骤： 步骤1:在线路上安装故障数据采集器，采集故障发生后前后各两个周波的Ξ相电压数 据； 步骤2:对步骤1所得到的共四个周波的Ξ相电压行波进行凯伦贝尔变换，得到线模电 压分量，其中，凯伦贝尔变换为，式中m，U2为线模电压，U0为零 模电压，Ua,Ub,Uc为各相电压； 步骤3:选取步骤2得到的线模电压分量进行化化ert变换，若配电网发生A相或B相接地 故障则选择线模电压m，若配电网发生C相接地故障则选择线模电压U2; 步骤4:对步骤3得到的变换后的结果进行中屯、差分求导，然后取绝对值得到模值序列 Uh，其中，中屯、差分求导公式为：，h(n+l)和h(n-l)分别为经过 HHbed变换得到的线模电压数据在时刻n+1和时刻n-1的值，h'（n)为经过求导后在时刻η 的值； 步骤5:对模值序列化中的第一个周波数据Xh进行判断，如果Xh(n)〉0.5max化h}，则将时 亥Ijn的模值视为脉冲干扰置零，得到去除干扰的第一个周波模值序列，其中，Xh(n)为Xh在时 亥。11的值； 步骤6:将化中低于阔值K的模值视为噪声置零，其中Κ= 1.5max{Xh}; 步骤7:在步骤6置零后的结果中选取有效区间，区间内模最大值为行波波头到达时刻。3. 根据权利要求2所述的配电网单相接地故障行波波头的计算方法，其特征在于： 所述步骤（7)中的有效区间选取，W第一个模值大于0.5max{Uh}的点为起点，如果在 0.01个周波时长内至少有2个模值大于0.5max{Uh}，则标定为有效区间，否则，W下一个模 值大于0. SmadUh}的点为起点继续判断直至选出有效区间。
【文档编号】G01R31/08GK106093705SQ201610424490
【公开日】2016年11月9日
【申请日】2016年6月15日 公开号201610424490.7, CN 106093705 A, CN 106093705A, CN 201610424490, CN-A-106093705, CN106093705 A, CN106093705A, CN201610424490, CN201610424490.7
【发明人】齐郑, 林健雄, 郑宇航, 庄舒仪
【申请人】华北电力大学