CN107219501B - 一种基于gis局部放电的精确定位方法 - Google Patents
一种基于gis局部放电的精确定位方法 Download PDFInfo
- Publication number
- CN107219501B CN107219501B CN201710291719.9A CN201710291719A CN107219501B CN 107219501 B CN107219501 B CN 107219501B CN 201710291719 A CN201710291719 A CN 201710291719A CN 107219501 B CN107219501 B CN 107219501B
- Authority
- CN
- China
- Prior art keywords
- partial discharge
- distance
- sensor
- gis
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 6
- 238000000691 measurement method Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000010845 search algorithm Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- 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/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Testing Relating To Insulation (AREA)
- Locating Faults (AREA)
Abstract
本发明公开了一种基于GIS局部放电的精确定位方法,属于GIS局部放电监测的方法。该方法获取局放后各传感器同步获取电磁波到达时刻后,对利用扩展时差定位法得到的测距结果引入Euclidean测度距离矩阵得到初步局放位置,进而求得对应的局放发生时间。融合由多传感器信息计算出的局放位置,得到最终的局放距离。本发明设计合理,可以实时监测情况下利用多个测量点信息进行综合局部放电定位。本发明的定位方法可以应用于GIS局部放电精确定位中。
Description
技术领域
本发明属于局部放电在线监测领域,尤其涉及一种基于GIS局部放电的精确定位方法。
背景技术
基于特高频的局部放电定位方法,其最基本的原理是采用电磁波到达的时间差关系。相应的计算方法主要包括搜索算法和迭代算法。牛顿迭代算法的核心是计算的收敛问题,相关的研究包括重庆大学的唐炬提出的taylor遗传定位算法,清华大学的高文胜等人提出的基于能量累积的空间搜索算法,华北电力大学唐志国提出的全局搜索迭代算法等。搜索算法(模式识别算法)是将被检测物划分为多个区域,将布置的传感器与各个区域之间的距离差作为判定模式。如果某区域对应的标准模式与判定模式相同,则该区域为故障区域。国内外在这个方面开展了很多研究,Hasan Reza Mirzaei等人为修正搜索算法中的标准模式,提出了微粒群优化算法。李立学、张晓星、王辉等人通过采用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=[T1 T2 T3…TK]T,K个传感器分别对应节点为:S=[S1S2…SK]。
Step3.从K组时间数据TM=[T1 T2 T3…TK]T找到其中的最小值Tmin,并找到与其对应的节点Sm。
Step4.由节点Sm的相邻的w个节点SV=[Sv1 Sv2…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=[L1f L2f … LKf]T=c(TM-t0η) (6)
其中,
可进一步得到由各传感器计算出的局放距离矩阵:
S=[S1 S2 … 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%。
Claims (5)
1.一种基于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,融合多传感器计算出的局放位置,得到最终的局放距离,进而完成精确局放定位。
2.根据权利要求1所述的一种基于GIS局部放电的精确定位方法,其特征在于:在步骤3中,计算欧式测度距离矩阵E具体如下:利用Euclidean距离测度法从w×(K-1)组测距结果与电磁波传播速度相除求得的电磁波传输时间与局放后传感器同步获取的电磁波到达时间得到两组向量的Euclidean长度,通过两个向量的距离测度得到欧式距离测度矩阵E。
3.根据权利要求2所述的一种基于GIS局部放电的精确定位方法,其特征在于:在步骤3中,欧式测度距离矩阵E具体如下:
其中,E(k-1)w表示对应的Euclidean距离,D(K-1)w为传输时间矩阵,K为传感器的个数,TM为局放放电发生后各传感器获取的电磁波初始到达时间。
4.根据权利要求1所述的一种基于GIS局部放电的精确定位方法,其特征在于:在步骤5中,局放发生的时间t0的具体计算如下:
其中,Tmin为局放后各传感器获取的电磁波信号到达时间矩阵中的最小值,Lbf’为初步局放位置f’到Tmin对应的b号传感器的最短距离,c为电磁波传播速度。
5.根据权利要求1所述的一种基于GIS局部放电的精确定位方法,其特征在于:在步骤6中,最终的局放距离具体计算如下:
其中,Sact为距离所定位的区段首端的最终的局放距离,Si为第i号传感器距离所定位的区段首端的距离,K为GIS***中放置的传感器总个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710291719.9A CN107219501B (zh) | 2017-04-28 | 2017-04-28 | 一种基于gis局部放电的精确定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710291719.9A CN107219501B (zh) | 2017-04-28 | 2017-04-28 | 一种基于gis局部放电的精确定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107219501A CN107219501A (zh) | 2017-09-29 |
CN107219501B true CN107219501B (zh) | 2019-12-24 |
Family
ID=59943622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710291719.9A Active CN107219501B (zh) | 2017-04-28 | 2017-04-28 | 一种基于gis局部放电的精确定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107219501B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108279365B (zh) * | 2018-04-18 | 2024-05-07 | 国网江苏省电力有限公司电力科学研究院 | 一种gil管廊内空间电磁波监测定位装置、***及方法 |
CN109932627A (zh) * | 2019-04-30 | 2019-06-25 | 国网山东省电力公司烟台供电公司 | 一种gis耐压试验局部放电的定位方法及*** |
CN112363030A (zh) * | 2020-10-30 | 2021-02-12 | 上海均赫检测技术有限公司 | 一种运行状态下低压母线槽绝缘检测方法 |
CN112763868B (zh) * | 2020-12-26 | 2022-11-08 | 广东电网有限责任公司电力科学研究院 | 一种基于免疫粒子群的局部放电源定位方法和*** |
CN112986766B (zh) * | 2021-02-25 | 2023-09-15 | 西安西电开关电气有限公司 | 一种局部放电定位方法、装置、存储介质和设备 |
CN116184141B (zh) * | 2023-04-25 | 2023-08-29 | 南昌工程学院 | 一种气体绝缘设备故障诊断方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2574941A2 (en) * | 2011-09-29 | 2013-04-03 | Hyundai Heavy Industries Co., Ltd. | Method of estimating location of partial discharge using signal attenuation rate in GIS tank |
JP2015078882A (ja) * | 2013-10-16 | 2015-04-23 | 株式会社東芝 | 絶縁診断装置 |
CN105911437A (zh) * | 2016-04-11 | 2016-08-31 | 国网天津市电力公司 | 一种精确定位gis局部放电的方法 |
CN104714156B (zh) * | 2015-02-03 | 2017-08-04 | 华北电力大学 | 基于特高频放电在线监测数据的gis局部放电定位方法 |
-
2017
- 2017-04-28 CN CN201710291719.9A patent/CN107219501B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2574941A2 (en) * | 2011-09-29 | 2013-04-03 | Hyundai Heavy Industries Co., Ltd. | Method of estimating location of partial discharge using signal attenuation rate in GIS tank |
JP2015078882A (ja) * | 2013-10-16 | 2015-04-23 | 株式会社東芝 | 絶縁診断装置 |
CN104714156B (zh) * | 2015-02-03 | 2017-08-04 | 华北电力大学 | 基于特高频放电在线监测数据的gis局部放电定位方法 |
CN105911437A (zh) * | 2016-04-11 | 2016-08-31 | 国网天津市电力公司 | 一种精确定位gis局部放电的方法 |
Non-Patent Citations (2)
Title |
---|
Application of MV GIS partial discharge measurement and location;Min-Yen Chiu et al.;《2012 IEEE International Conference on Condition Monitoring and Diagnosis》;20120927;第462-465页 * |
GIS 中典型局部放电的频谱特征及传播特性;丁登伟 等;《高电压技术》;20141031;第40卷(第10期);第3243-3251页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107219501A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107219501B (zh) | 一种基于gis局部放电的精确定位方法 | |
CN107044852B (zh) | 不平整状态下全站仪测量方法 | |
CN103674015B (zh) | 一种无轨化定位导航方法及装置 | |
CN104507050B (zh) | 一种WiFi室内定位中概率型指纹匹配方法 | |
CN102740456B (zh) | 一种无线传感器网络的节点定位方法 | |
CN104080165B (zh) | 一种基于tdoa的室内无线传感器网络定位方法 | |
CN104143194A (zh) | 一种点云分割方法及装置 | |
CN108253976B (zh) | 一种充分借助车辆航向的三阶段在线地图匹配算法 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN109682336B (zh) | 用于车身精度检测的三坐标测量路径自动规划与优化方法 | |
CN103557792B (zh) | 一种空中加油锥套目标的视觉跟踪与位置测量方法 | |
CN106646412B (zh) | 一种用于隧道内的多雷达时间同步处理方法 | |
CN103954842B (zh) | 一种确定大型接地***接地电阻值的方法 | |
CN107621263A (zh) | 一种基于道路磁场特征的地磁定位方法 | |
CN107976182A (zh) | 一种多传感器融合建图***及其方法 | |
CN104684081A (zh) | 基于距离聚类遴选锚节点的无线传感器网络节点定位算法 | |
CN103278126A (zh) | 一种基于最小区域的零件球度误差评定方法 | |
CN105137412A (zh) | 一种2d激光雷达距离图像中线段特征精确拟合方法 | |
CN104155979A (zh) | 一种基于磁场对称性的磁导航智能车定位装置及方法 | |
CN106208050B (zh) | 一种基于pmu的电网支路静态参数检测辨识方法 | |
CN103294920B (zh) | 一种基于最小区域的零件圆柱度误差评定方法 | |
CN103344844A (zh) | 电网安全警示终端及其警示方法 | |
CN109239533A (zh) | 一种基于人工神经网络的特高压直流输电线路的故障定位方法 | |
CN112828883B (zh) | 一种未知环境下的机器人环境探索方法及*** | |
CN112858842A (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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 221116 Xuzhou Tongshan University Road, Xuzhou City, Jiangsu Province, Nanhu campus, China University of Mining and Technology Applicant after: China University of Mining & Technology Address before: 221116 Xuzhou University Road, Jiangsu, No. 1 Applicant before: China University of Mining & Technology |
|
GR01 | Patent grant | ||
GR01 | Patent grant |