一种基于GIS局部放电的精确定位方法与流程

本发明属于局部放电在线监测领域，尤其涉及一种基于gis局部放电的精确定位方法。

背景技术：

基于特高频的局部放电定位方法，其最基本的原理是采用电磁波到达的时间差关系。相应的计算方法主要包括搜索算法和迭代算法。牛顿迭代算法的核心是计算的收敛问题，相关的研究包括重庆大学的唐炬提出的taylor遗传定位算法，清华大学的高文胜等人提出的基于能量累积的空间搜索算法，华北电力大学唐志国提出的全局搜索迭代算法等。搜索算法(模式识别算法)是将被检测物划分为多个区域，将布置的传感器与各个区域之间的距离差作为判定模式。如果某区域对应的标准模式与判定模式相同，则该区域为故障区域。国内外在这个方面开展了很多研究，hasanrezamirzaei等人为修正搜索算法中的标准模式，提出了微粒群优化算法。李立学、张晓星、王辉等人通过采用bp算法、遗传算法、gk模糊分类算法，提取高位空间非线性特征，来对放电模型进行识别，结果表明具有良好的识别效果和鲁棒性。

在定位计算中，首波到达时刻或时间差会造成误差，可能会导致迭代算法无法局部收敛，而采用搜索算法得到的理论结果存在与实际的时间差不能相匹配的问题。这些问题都会导致定位不准确或失效。如何利用多个融合传感器的信息得出精确的定位结果，是特高频定位方法的核心问题。

技术实现要素：

本发明所要解决的技术问题是针对背景技术的不足提供了一种基于gis局部放电的精确定位方法，可以利用大量现有的投运检测设备，实现简单，具有较强的经济性和较好的实用价值。

本发明为解决上述技术问题采用以下技术方案

一种基于gis局部放电的精确定位方法，具体包含如下步骤；

步骤1，获取局放电磁波到达各个传感器的时间，并找到其中时间最小值对应的节点sm；

步骤2，确定节点sm到每个不相邻节点的最短路径，并根据扩展时差法将节点sm处的传感器作为扩展时差定位的一端，分别与其它k-1个传感器构成时差定位组合，计算出节点sm通过和每个相邻节点组成的最短路径中的故障局放位置，其中，k取整数且k≥1；

步骤3，根据w个相邻节点决定的路径计算得到w组测距结果，进而获取w×(k-1)个测距结果，并计算欧式测度距离矩阵e，其中w为正整数；

步骤4，找到欧式测度距离矩阵e中的最小值，以及其对应的行列来确定初步局放位置f'；

步骤5，根据初步局放位置f'计算局放发生的时间t0；

步骤6，融合多传感器计算出的局放位置，得到最终的局放距离，进而完成精确局放定位。

作为本发明一种基于gis局部放电的精确定位方法的进一步优选方案，在步骤2中，计算欧式测度距离矩阵e具体如下：利用euclidean距离测度法从w×(k-1)组测距结果与电磁波传播速度相除求得的电磁波传输时间与局放后传感器同步获取的电磁波到达时间得到两组向量的euclidean长度，通过两个向量的距离测度得到欧式距离测度矩阵e。

作为本发明一种基于gis局部放电的精确定位方法的进一步优选方案，欧式测度距离矩阵e具体如下：

其中，e(k-1)w表示对应的euclidean距离，d(k-1)w为传输时间矩阵，k为传感器的个数，tm为局放放电发生后各传感器获取的电磁波初始到达时间。

作为本发明一种基于gis局部放电的精确定位方法的进一步优选方案，局放发生的时间t0的具体计算如下：

其中，tmin为局放后各传感器获取的电磁波信号到达时间矩阵中的最小值，lbf’为初步局放位置f’到tmin对应的b号传感器的最短距离。

作为本发明一种基于gis局部放电的精确定位方法的进一步优选方案，在步骤6中，最终的局放距离具体计算如下：

其中，sact为距离所定位的区段首端的最终的局放距离，si为第i号传感器距离所定位的区段首端的距离，k为gis系统中放置的传感器总个数

有益效果

本发明采用以上技术方案与现有技术相比，具有以下技术效果：

1、本发明在同步获取局部放电产生的电磁波初始波头到达各传感器的时间后，利用扩展时差定位法，融合各传感器的波头信息及gis固有拓扑结构，引用euclidean距离得到初步测距结果，进而求得对应的局放发生时间，最后利用此局放时间进一步提高了测距精度；

2、本发明经仿真验证，相比于普通双端时差定位方法，该方法通过局放发生时间对测距结果的修正得到了较高的测距精度，通过人为模拟设置错误监测数据，验证了算法具有较好的容错性，可靠性较高，易于计算机实现；适应于复杂的大型电网，保护范围无死区，定位速度迅速准确，保护的有效范围为全网系统，不存在定位死区；

3、本发明可以利用大量现有的投运检测设备，实现简单，具有较强的经济性和较好的实用价值。

附图说明

图1是本发明的gis仿真模型；

图2是本发明的扩展时差定位法原理图；

图3是本发明的sm相邻传感器示意图；

图4是本发明的gis系统拓扑结构的无向图；

图5是本发明的具体实施方案流程图。

具体实施方式

下面结合附图对本发明的技术方案做进一步的详细说明：

本发明的目的是提出一种gis实时监测下局部放电的精确定位方法，如图5所示，包含以下步骤：

step1.将gis各个节点编号，节点间的拓扑关系、线路长度绘制成表格，根据节点间的拓扑关系画出此gis系统的无向图。通过已有的优化配置方法得到优化后的节点布置方案。

step2.假设局放发生时，k个传感器同步检测到电磁波信号，各传感器可分别获取电磁波初始到达时间的数据：tm＝[t1t2t3…tk]^t，k个传感器分别对应节点为：s＝[s1s2…sk]。

step3.从k组时间数据tm＝[t1t2t3…tk]^t找到其中的最小值tmin，并找到与其对应的节点sm。

step4.由节点sm的相邻的w个节点sv＝[sv1sv2…svw]确定最短路径并根据扩展时差定位法得到w×(k-1)个测距结果，利用euclidean距离测度法从w×(k-1)组测距结果与电磁波传播速度相除求得的电磁波传输时间与局放后传感器同步获取的电磁波到达时间得到两组向量的euclidean长度，通过两个向量的距离测度得到欧式距离测度矩阵e，

其中，euclidean距离：

电磁波从各组测距结果中的局放位置到各传感器的传输时间dij＝[d1’d2’…dk’]可由最短路径求得，将w×(k-1)组传输时间用矩阵df表示如下：

step5.从euclidean距离测度矩阵e中找到其中的最小值，以及其对应的行列来确定初步局放位置f’，根据初步局放位置f’，计算局放发生的时间t0。

step6.融合多传感器计算出的局放位置完成精确局放定位。

所述的扩展时差定位法为：如图2，i和j为线路的始端和末端，电磁波透过i和j到达相邻传感器x和y。将x到y作为一条扩展的定位区段始端和末端，使用x和y的电磁波到达时刻计算局放发生位置的公式为：

其中：tx和ty分别为电磁波到达x和y的时间，c为电磁波的传播速度，lfx和lxy分别为局部放电点f到x的距离和x与y的距离。

在图3中，将节点sm处的传感器作为扩展时差定位的一端，分别与其它k-1个传感器构成时差定位组合。节点sm与节点sk处的传感器构成时差定位组合，sm、f、sk三者之间的最短距离为sm经f点到达sk，此距离中包含与节点sm直接相邻的一条线路。即sm、sk构成时差定位组合计算时，所采用的最短路径经过sm相邻的一条节点svr，根据w个相邻节点决定的路径由公式(4)计算得到w组测距结果，对于全部k-1个测距组合可以得到w×(k-1)个测距结果。

所述的融合多传感器计算出的局放位置完成的精确局放定位方法为：根据初步局放位置f’，由公式(5)计算出局放发生时间t0，

其中，tmin为局放后各传感器获取的电磁波信号到达时间矩阵中的最小值，lbf’为初步局放位置f’到tmin对应的b号传感器的最短距离。

如图1所示，再结合各传感器获取的电磁波到达时间，融合多传感器计算出的局放位置，得到局放位置f到各传感器的最短距离构成的矩阵为：

l＝[l1fl2f…lkf]^t＝c(tm-t0η)(6)

其中，

可进一步得到由各传感器计算出的局放距离矩阵：

s＝[s1s2…sk]^t(8)

融合各传感器信息，则局放距离最终确定为：

本发明的目的在于弥补现有技术的不足之处,提供一种设计合理、准确可靠的精确定位局部放电的方法。

本次算例使用500kvgis实验平台，将gis各个节点编号，节点间的拓扑关系、线路长度如表1所示。

表1

根据节点间的拓扑关系画出此gis系统的无向图，如图4所示。最终通过已有的优化配置方法求得需要配置传感器的节点为1,6,11,18,25。即为监测整个gis系统只需要在1,6,11,18,25节点处安装传感器。

假设在20ns时刻，距离8节点1000mm处，区段8-9上发生局部放电。通过各传感器能够监测到各个配置点的电压信号时域波形，对各传感器检测到的电磁波信号进行400mhz的低通滤波，滤波后信号中仅含有tem波。可以发现低通滤波后，各传感器检测的信号的幅值和波形的变化较小，可以轻易获取局部放电产生电磁波到达各传感器的初始时间tm，如表2所示：

表2

从传感器监测数据中可知，tmin＝t6。6号传感器分别与其他传感器构成时差定位组合，得到的初步测距结果如表3所示：

表3

表中，euclidean距离最小为6号传感器与11号传感器的组合，因此，以6号和11号得到的测距结果作为以后计算的基础，得到初步局放位置为区段8-9上距8号节点923.5mm，实际局放距离为1000mm,绝对误差为76.5mm，相对误差为7.65％。利用公式(5)求得局放发生时间t0为20.0283ns。

进一步利用公式(6)(8)得到校正后的局放距离如表4所示：

表4

最终，由公式(9)得到局放位置为区段8-9上距8节点962.7mm处，而实际局放位置为1000mm处，故绝对误差为37.3mm，相对误差为3.73％。