一种电缆故障识别方法与流程

1.本发明涉及电缆故障识别技术领域，尤其涉及一种电缆故障识别方法。


背景技术：

2.电缆作为重要的基础性配套产业被誉为城市的“神经”和“血管”，为保证电缆工作时的安全可靠、安装更换方便以及城市美观，大多采用直埋敷设的方式。地下恶劣的环境、过负荷的使用以及自然老化都会导致电缆发生故障，且很难发现故障位置。因此，精确判别故障类型和定位故障点是至关重要的。
3.行波法作为目前常用的故障检测方法，它的工作原理是在电缆的一端发射激励信号，激励信号遇到故障点发生反射，通过入射信号和反射信号的时间差以及信号的传播速度来计算故障的位置。行波法常见的方法有：时域反射法、频域反射法、序列时域反射法、扩展频谱时域反射法等。tdr实现简单，且对电缆无损害，但是在测量过程中容易产生噪声信号的叠加，给反射信号提取带来困难，且测试信号随着传播距离的增大衰减会逐渐增大，很难识别反射信号，导致测量结果误差很大。随之兴起的fdr在tdr基础上有一定的改进，激励采用扫频信号，测试精度比tdr高，但是它只能检测高阻故障，且设备的体积较大成本较高。sstdr向电缆发送经余弦调制的伪随机码作为激励信号，该信号与电缆中的有效信号互不干扰，但是检测范围受码周期的限制，需要昂贵的伪随机码发生器，成本较高。stdr以其定位精度高、抗干扰能力强、实现简单、适用范围广而成为研究的重点。然而，传统stdr 只能识别基本的短路与开路故障，而不能识别诸如高阻、低阻等其它类型的故障，因此在精确的电缆故障诊断应用中受阻。


技术实现要素：

4.本发明的目的在于提供一种电缆故障识别方法，解决背景技术中提到的技术问题。
5.一种电缆故障识别方法，所述方法包括如下步骤：
6.步骤1：设置及电缆参数模型，根据参数模型分析出电缆故障的种类；
7.步骤2：生成两路相同检测信号，一路作为参考信号，另一路作为测试信号注入待测电缆，测试信号通过t型头注入待测故障电缆；
8.步骤3：当电缆发生故障时，故障点处的阻抗与电缆的特征阻抗不匹配，激励信号在电缆内传输时遇到阻抗不匹配的点时发生反射；
9.步骤4：将采集到的参考信号与反射信号进行互相关运算分析，根据反射信号的延迟特性和幅值特性，来识别故障类型和定位故障位置。
10.进一步地，所述步骤1的具体过程为：设l0、r0、c0、g0分别为单位长度长线的分布电感、分布电阻、分布电容和分布电导，根据基尔霍夫电流定律和电压定律，距离电缆初始端为x处的电压电流方程为：
[0011][0012]
根据初始条件和边界条件求出电压电流方程的解，得到距离电缆x处入射信号的反射系数为：
[0013][0014]
其中，z0为电缆故障点的阻抗，zl为电缆的特征阻抗，分析反射情况的类型为，电缆正常，无故障：z0＝zl，г＝0，无反射存在，低阻故障，z0<10z， z为电缆的波阻抗，不超过40ω，，г在(
‑
1，0)之间，信号发生部分反射，且入射信号与反射信号的极性相反，其中z0＝0，г＝
‑
1为短路故障，电缆开路， z0＝∞，г＝1，信号发生全反射，且入射信号与反射信号的极性相同，高阻故障，z0>10z，г在(0，1)之间，且入射信号与反射信号的极性相同。
[0015]
进一步地，所述步骤2中生成检测信号具体过程为：通过若干个线性移位寄存器，并通过移位时钟脉冲在线性移位寄存器进行移动生成检测信号，生成算式为:
[0016][0017]
c
i
表示反馈系数，i为正整数,x
i
表示元素的相应位置。
[0018]
对检测信号进行功率谱分析，分析算式为：
[0019][0020]
表示不同谱线，a表示序列电平大小，
[0021]
其中，n为检测信号的长度，fc/n为检测信号的谱线间距，当时钟频率 f
c
不变，信号长度n增大时，谱线之间的间距减小，强度下降，但不会改变在指定频率内的总功率，在相同点数下的检测信号自相关曲线随着阶数的增加，其自相关更好、更容易辨识，。
[0022]
进一步地，所述步骤2中，有控制信号发生器产生的检测信号经过功率分配器后被复制成两路相同的信号，其中一路信号作为参考信号s
ref
送入示波器的1端口，另一路通过t型连接头时分成两部分，一部分作为激励信号进入电缆，另一部分作为基础信号s0送入示波器的2端口。
[0023]
进一步地，所述步骤3中，激励信号进入到电缆内一旦遇到故障点，就会返回到发射端口，反射信号与激励信号之间存在延迟，延迟时间为t型连接头到故障点的往返时间τ，反射信号在返回到t型连接头后会分成两路，一路反射信号s’ret
经过功率分配器后被示波器端口1采集，另一部分反射信号s
ret
直接被示波器端口2采集。
[0024]
进一步地，所述步骤4中，pc机采集示波器两个端口的数据进行互相关运算，经过分析计算出信号往返时间τ，通过算式计算出故障位置。
[0025]
进一步地，互相关运算与分析计算的过程为:
[0026]
假设控制信号发生器产生的信号为c(t)，则示波器接收到的各个信号为： s
ref
＝k1c(t
‑
τ1)，s0＝k2c(t
‑
τ2)，s
ret
＝k3c(t
‑
τ3)，s’ret
＝k4c(t
‑
τ4)，其中τ1、τ2、τ3、τ4分为参考信号、基础信号、反射信号1、反射信号2从功率分配器到示波器的传播时间，k1、k2、k3、k4为相对于信号发生器输出的衰减系数，将组合信号s
ref
+s’ret
与s0+s
re
t进行相关运算可得：
[0027][0028]
上式(5)包含三项，其中第一项为参考信号自相关峰值，第二项为反射信号2互相关的峰值，第三项为反射信号1互相关的峰值，即衰减系数k4>k3，假设信号从输出端到t型连接头的传输时间为τ5，则τ3‑
τ1＝τ，τ4‑
τ2＝τ+2τ5，得出反射信号1的传播时间对于定位故障点更精确。
[0029]
进一步地，计算出故障位置的具体过程为:
[0030]
假设入射信号遇到故障返回至入射点的时间为t，且入射信号在电缆内的传播速度为v，则故障距离l可由下式表示：
[0031][0032]
然而当故障点与测距点小于固定值时，测试信号传播时间短，发射信号与反射信号混叠在一起，致使无法测距，因此，被测电缆的测试距离范围通常设定为：
[0033][0034]
其中，t
s
＝1/f
s
为一个码元持续时间，f
s
为码元生成速率，n为伪随机序列的长度，c为光速。
[0035]
进一步地，步骤4中还包括建立基于互相关峰值的电缆负载阻抗识别模型，具体过程为，将电缆加载不同的负载阻抗，得到不同负载阻抗下反射信号1与参考信号的互相关信号，互相关信号包括参考信号的自相关峰值、反射信号2与基础信号的互相关峰值和反射信号1与参考信号的互相关峰值，互相关峰值进行处理：
[0036][0037]
其中，r
z
为某负载下的互相关峰值，min(r
xy
)、max(r
xy
)分别为不同负载下互相关峰值的最大值和最小值，反射系数г的取值范围在[
‑
1,1]，
[0038]
其次，将式(8)代入反射系数公式(2)中，得到互相关峰值rz对应的负载阻抗的z0估算公式：
[0039]
[0040]
其中，z
l
为电缆的特征阻抗，随负载电阻由小变大，对应的互相关峰值由负数逐渐变为正数，利用式(9)计算得到负载阻抗的估算值，并与真实的负载阻抗做相关性分析，得到负载阻抗的估算值与真实值的相关度为98.80％。
[0041]
本发明采用了上述技术方案，本发明具有以下技术效果：
[0042]
本发明判断故障点处的负载大小，进而对该故障类型进行判定，因此，可以推断该方法能够实现电缆故障类型的精确识别，可用于精确的电缆故障类型识别和故障点定位，通过对不同长度同轴电缆进行故障实验测试，结果表明，该方法能够实现故障点的精确定位，通过对同一长度同轴电缆反射信号互相关峰值分析表明，其与不同负载之间具有很强的相关性，可实现对开路、短路等故障精确识别。
附图说明
[0043]
图1为本发明原理示意图。
[0044]
图2为本发明功率谱分析中的时域波形图。
[0045]
图3为本发明功率谱分析中的功率谱密度图。
[0046]
图4为本发明功率谱分析中的自相关曲线图。
[0047]
图5为本发明开路故障与短路故障的峰值对比图。
[0048]
图6为本发明不同负载对应的峰值图。
具体实施方式
[0049]
为使本发明的目的、技术方案及优点更加清楚明白，举出优选实施例，对本发明进一步详细说明。然而，需要说明的是，说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解，即便没有这些特定的细节也可以实现本发明的这些方面。
[0050]
一种电缆故障识别方法，如图1所示，所述方法包括如下步骤：
[0051]
步骤1：设置及电缆参数模型，根据参数模型分析出电缆故障的种类。
[0052]
设l0、r0、c0、g0分别为单位长度长线的分布电感、分布电阻、分布电容和分布电导，根据基尔霍夫电流定律和电压定律，距离电缆初始端为x处的电压电流方程为：
[0053][0054]
根据初始条件和边界条件求出电压电流方程的解，得到距离电缆x处入射信号的反射系数为：
[0055][0056]
其中，z0为电缆故障点的阻抗，zl为电缆的特征阻抗，分析反射情况的类型为，电缆正常，无故障：z0＝zl，г＝0，无反射存在，低阻故障，z0<10z， z为电缆的波阻抗，不超过40ω，，г在(
‑
1，0)之间，信号发生部分反射，且入射信号与反射信号的极性相反，其中z0＝0，г＝
‑
1为短路故障，电缆开路， z0＝∞，г＝1，信号发生全反射，且入射信号与反射信号
的极性相同，高阻故障，z0>10z，г在(0，1)之间，且入射信号与反射信号的极性相同。
[0057]
步骤2：生成两路相同检测信号，一路作为参考信号，另一路作为测试信号注入待测电缆，测试信号通过t型头注入待测故障电缆。
[0058]
通过若干个线性移位寄存器，并通过移位时钟脉冲在线性移位寄存器进行移动生成检测信号，生成算式为:
[0059][0060]
对检测信号进行功率谱分析，分析算式为：
[0061][0062]
其中，n为检测信号的长度，fc/n为检测信号的谱线间距，当时钟频率 fc不变，信号长度n增大时，谱线之间的间距减小，强度下降，但不会改变在指定频率内的总功率，在相同点数下的检测信号自相关曲线随着阶数的增加，其自相关更好、更容易辨识。
[0063]
其中，具体的分析图如图2
‑
4所示，因此，本文采用十阶检测信号作为电缆故障检测的激励信号。
[0064]
控制信号发生器产生的检测信号经过功率分配器后被复制成两路相同的信号，其中一路信号作为参考信号s
ref
送入示波器的1端口，另一路通过t 型连接头时分成两部分，一部分作为激励信号进入电缆，另一部分作为基础信号s0送入示波器的2端口。
[0065]
步骤3：当电缆发生故障时，故障点处的阻抗与电缆的特征阻抗不匹配，激励信号在电缆内传输时遇到阻抗不匹配的点时发生反射。激励信号进入到电缆内一旦遇到故障点，就会返回到发射端口，反射信号与激励信号之间存在延迟，延迟时间为t型连接头到故障点的往返时间τ，反射信号在返回到t 型连接头后会分成两路，一路反射信号s’ret
经过功率分配器后被示波器端口1 采集，另一部分反射信号s
ret
直接被示波器端口2采集。
[0066]
步骤4：将采集到的参考信号与反射信号进行互相关运算分析，根据反射信号的延迟特性和幅值特性，来识别故障类型和定位故障位置。pc机采集示波器两个端口的数据进行互相关运算，经过分析计算出信号往返时间τ，通过算式计算出故障位置。
[0067]
互相关运算与分析计算的过程为:
[0068]
假设控制信号发生器产生的信号为c(t)，则示波器接收到的各个信号为： s
ref
＝k1c(t
‑
τ1)，s0＝k2c(t
‑
τ2)，s
ret
＝k3c(t
‑
τ3)，s’ret
＝k4c(t
‑
τ4)，其中τ1、τ2、τ3、τ4分为参考信号、基础信号、反射信号1、反射信号2从功率分配器到示波器的传播时间，k1、k2、k3、k4为相对于信号发生器输出的衰减系数，将组合信号s
ref
+s’ret
与s0+s
re
t进行相关运算可得：
[0069][0070]
上式(5)包含三项，其中第一项为参考信号自相关峰值，第二项为反射信号2互相关的峰值，第三项为反射信号1互相关的峰值，即衰减系数k4>k3，假设信号从输出端到t型连接头的传输时间为τ5，则τ3‑
τ1＝τ，τ4‑
τ2＝τ+2τ5，得出反射信号1的传播时间对于定位故障点更精确。
[0071]
计算出故障位置的具体过程为:
[0072]
假设入射信号遇到故障返回至入射点的时间为t，且入射信号在电缆内的传播速度为v，则故障距离l可由下式表示：
[0073][0074]
然而当故障点与测距点小于固定值时，测试信号传播时间短，发射信号与反射信号混叠在一起，致使无法测距，因此，被测电缆的测试距离范围通常设定为：
[0075][0076]
其中，t
s
＝1/f
s
为一个码元持续时间，f
s
为码元生成速率，n为伪随机序列的长度，c为光速。
[0077]
步骤4中还包括建立基于互相关峰值的电缆负载阻抗识别模型，具体过程为，将电缆加载不同的负载阻抗，得到不同负载阻抗下反射信号1与参考信号的互相关信号，互相关信号包括参考信号的自相关峰值、反射信号2与基础信号的互相关峰值和反射信号1与参考信号的互相关峰值，互相关峰值进行处理：
[0078][0079]
其中，r
z
为某负载下的互相关峰值，min(r
xy
)、max(r
xy
)分别为不同负载下互相关峰值的最大值和最小值，反射系数г的取值范围在[
‑
1,1]，
[0080]
其次，将式(8)代入反射系数公式(2)中，得到互相关峰值rz对应的负载阻抗的z0估算公式：
[0081][0082]
其中，z
l
为电缆的特征阻抗，随负载电阻由小变大，对应的互相关峰值由负数逐渐变为正数，利用式(9)计算得到负载阻抗的估算值，并与真实的负载阻抗做相关性分析，得到负载阻抗的估算值与真实值的相关度为98.80％。
[0083]
在实验室利用对syv50
‑3‑
1同轴电缆在50m、100m、150m、200m处分别进行开路故障和短路故障的模拟测试。实验中使用的测试信号为十阶的检测信号，其码元生成率为30mhz，测试点数为12500。
[0084]
电缆故障的识别和定位，图5为电缆故障的识别和定位结果，其中图5 (a)为同轴
电缆在50m、100m、150m、200m处的开路故障测试结果。可以看出反射信号互相关峰值与参考信号自相关的峰值极性相同，因此由式(2) 判断出电缆是开路故障。下表1实际开路故障位置与测量故障位置的比较。
[0085]
图5(b)为同轴电缆在50m、100m、150m、200m处的短路故障测试结果。可以看出反射信号互相关峰值与参考信号自相关的峰值极性相反，因此由式(2)判断出电缆是短路故障。下表2实际短路故障位置与测量故障位置的比较。
[0086]
此外，由于衰减系数k4>k3，因此反射信号2互相关峰值高于反射信号1 互相关峰值，即图5表示的各信号其左侧互相关峰值大于其右侧峰值。
[0087]
表1实际开路故障位置与测量故障位置的比较
[0088][0089]
表2实际短路故障位置与测量故障位置的比较
[0090][0091]
通过在200m同轴电缆的末端添加不同的负载进行电缆故障的模拟测试，得到如图6所示的不同阻抗下的实测响应，可以看出随负载由小变大，对应的互相关峰值也由负数逐渐变为正数。
[0092]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。