CN108761279B - 一种用于局部放电源定位的极坐标网格搜索方法及系统 - Google Patents

一种用于局部放电源定位的极坐标网格搜索方法及系统 Download PDF

Info

Publication number
CN108761279B
CN108761279B CN201810247930.5A CN201810247930A CN108761279B CN 108761279 B CN108761279 B CN 108761279B CN 201810247930 A CN201810247930 A CN 201810247930A CN 108761279 B CN108761279 B CN 108761279B
Authority
CN
China
Prior art keywords
grid
value
discharge source
partial discharge
sensor
Prior art date
Application number
CN201810247930.5A
Other languages
English (en)
Other versions
CN108761279A (zh
Inventor
薛峰
魏东亮
王植
李通
阙伟平
刘珂
刘勤锋
何俊达
林伯琴
张雅洁
吴秀全
胡岳
黄洪逸
Original Assignee
广东电网有限责任公司东莞供电局
上海交通大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 广东电网有限责任公司东莞供电局, 上海交通大学 filed Critical 广东电网有限责任公司东莞供电局
Priority to CN201810247930.5A priority Critical patent/CN108761279B/zh
Publication of CN108761279A publication Critical patent/CN108761279A/zh
Application granted granted Critical
Publication of CN108761279B publication Critical patent/CN108761279B/zh

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/12Testing 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
    • G01R31/1227Testing 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 of components, parts or materials

Abstract

本发明提供了一种用于局部放电源定位的极坐标网格搜索方法及系统,包括:利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形的网格;计算特征值数组(A1,A2,…An);查找特征值数组中最小值Am,最小值Am所在的网格即局放点所在的网格。判断最小值Am所在的网格的大小是否符合精度要求,若判断结果为是则作为最终解,若判断结果为否则将最小值Am所在的网格作为下一个待搜索的空间返回再次进行网格划分。本发明提供一种利用特高频时延信号在角度上准确定位局部放电源的方法;算法计算量小,时间性能优越,相较既有算法而言有较大提升;简单易行,对外部设备环境依赖低,检测准确度高。

Description

一种用于局部放电源定位的极坐标网格搜索方法及系统

技术领域

本发明涉及电力技术领域,具体地,涉及一种用于局部放电源定位的极坐标网格搜索方法及系统。

背景技术

局部放电是引起电力设备绝缘故障的重要原因,局部放电检测作为绝缘监测的重要内容,能够在发生绝缘介质击穿前及时排查到设备绝缘介质中存在的缺陷与劣化。高效的局放定位,有利于制定的检修方案,提高检修效率,减少停电检修时间与防止破坏性故障发生,对维护电力系统的稳定性有重要意义。

局部放电发生时,会伴随出现相应的声、光、热、电现象。基于这些特征量,人们发展出许多相应的检测局部放电的方法。如脉冲电流法、声测法、光测法以及化学检测法等非电量检测有着不易受电磁干扰影响、无需考虑电容等长处。但非电量检测存在灵敏度低,难以进行详细诊断,因此常作为辅助手段。

传统的局部放电检测方法的频率较低,随着技术发展,近年来推出了一个较为新颖的特高频检测方法。运动的电场与磁场会相互转化。电磁波由二者相互转化,并对外传播而产生。变电站中电设备一般在绝缘介质内部发生局部放电时,其放电范围较小,击穿时间很短,局部放电脉冲的上升沿十分陡,由此便会产生快速变化的电场,进而激发频率很高的电磁波,并可被设备外部特高频传感器所接收。故可利用传感器监测电气设备发生局放所产生的UHF信号,通过分析,获取局放定位信息,实现设备绝缘故障的诊断。

通过建立特高频传感器天线阵列,现有的分布式传感器特高频定位算法利用局部放电源产生的特高频信号到达四个特高频传感器的时间差建立空间坐标方程,通过求解出局放源的坐标,实现定位局放源。其中牛顿迭代法在存在误差时,迭代解的误差很大,甚至可以达到数千米的等级,且与牛顿迭代法本身性质有关,难以通过优化参数改变。网格搜索方法存在计算量大,方向角存在误差等缺点。

发明内容

针对现有技术中的缺陷,本发明的目的是提供一种用于局部放电源定位的极坐标网格搜索方法及系统。

根据本发明提供的一种用于局部放电源定位的极坐标网格搜索方法,包括:

网格划分步骤:利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形的网格,并分别标记为1,2,3…n,网格中心点记为Pk(rkk),k=1,2,3…n;

特征值数组计算步骤:特征值数组(A1,A2,…An)的计算公式包括:

其中,c是光速,τi是假设局部放电源位于(r,θ),其放电信号到达第i个传感器与到达第一个传感器的时间差,i是传感器序号,ri是第i个传感器的极坐标距离值,fi(R,θ)是局放源到第i个传感器与到第一个传感器的理论距离差与依照时延计算的距离差的误差值,Ak值为到三个传感器误差值的集合指标作为特征值反应定位误差大小情况;

定位步骤:查找特征值数组中最小值Am,最小值Am所在的网格即局放点所在的网格。

较佳的,还包括:

判断循环步骤:判断最小值Am所在的网格的大小是否符合精度要求,若判断结果为是则作为最终解,若判断结果为否则将最小值Am所在的网格作为下一个待搜索的空间返回网格划分步骤。

较佳的,在循环时新的r值为原半径除以缩小倍数,新的θ值为原角度除以角度缩小倍数。

较佳的,所述四个传感器成矩形放置。

根据本发明提供的一种用于局部放电源定位的极坐标网格搜索系统,包括:

网格划分模块:利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形的网格,并分别标记为1,2,3…n,网格中心点记为Pk(rkk),k=1,2,3…n;

特征值数组计算模块:特征值数组(A1,A2,…An)的计算公式包括:

其中,c是光速,τi是假设局部放电源位于(r,θ),其放电信号到达第i个传感器与到达第一个传感器的时间差,i是传感器序号,ri是第i个传感器的极坐标距离值,fi(R,θ)是局放源到第i个传感器与到第一个传感器的理论距离差与依照时延计算的距离差的误差值,Ak值为到三个传感器误差值的集合指标作为特征值反应定位误差大小情况;

定位模块:查找特征值数组中最小值Am,最小值Am所在的网格即局放点所在的网格。

较佳的,还包括:

判断循环模块:判断最小值Am所在的网格的大小是否符合精度要求,若判断结果为是则作为最终解,若判断结果为否则将最小值Am所在的网格作为下一个待搜索的空间返回再次进行网格划分。

较佳的,在循环时新的r值为原半径除以缩小倍数,新的θ值为原角度除以角度缩小倍数。

较佳的,所述四个传感器成矩形放置。

与现有技术相比,本发明具有如下的有益效果:

1、提供一种利用特高频时延信号在角度上准确定位局部放电源的方法;

2、算法计算量小,时间性能优越,相较既有算法而言有较大提升;

3、简单易行,对外部设备环境依赖低,检测准确度高。

附图说明

通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:

图1为本发明的流程图;

图2为本发明实施例的传感器放置示意图;

图3为本发明的模型示意图;

图4为本发明实施例的现场定位模型;

图5为本发明实施例的现场定位结果。

具体实施方式

下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。

如图1所示,本发明提供的一种用于局部放电源定位的极坐标网格搜索方法,包括:

网格划分步骤,利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形网格,并分别标记为1,2,3…n,网格中心点记为Pk(rkk),k=1,2,3…n;r与θ初始值为外加设定,在循环时新的r值为原搜索半径除以缩小倍数,新的θ值为原有搜索网格角度除以角度缩小倍数。

特征值数组(A1,A2,…An)计算步骤:如下式(1)(2)所列,其中c是光速,τi是传感器2,3,4与传感器1的局部放电信号到达时间差,ri是第i个传感器的极坐标距离值(即到原点的距离),fi(R,θ)是局放源到第i个传感器与到第一个传感器的理论距离差与依照时延计算的距离差的误差值,Ak值为到三个传感器误差值的集合指标作为特征值反应定位误差大小情况。

定位步骤,通过查找特征值数组中最小值Am,其所在的网格即局放点所在的网格点,并将此网格作为下一个待搜索的空间。

判断循环步骤,判断网格大小是否符合精度要求,若是则将此网格中心点作为最终解返回,否则循环搜索,进入下一步新的网格划分。

误差分析证明:

如附图2所示,将四个传感器呈矩形放置,1,2,3,4为同一水平面内呈矩形放置的四个传感器位置,1,2之间的距离为2r1,2,3之间的距离为2r2。并以矩形中心为原点建立如图坐标系。局放源坐标为(x,y,z),假设从局放源发出信号到1,2,3,4号传感器接收到信号所用的时间为t1,t2,t3,t4,2,3,4号传感器接收到信号的时间点相对于1号传感器滞后的时间为t21,t31,t41。则可建立方程组如下:

求解方程组(1)可分别得到

在一定的时延误差内,t21/t41变化小,同理加上包含距离信息的t1,(t21+2t1)/(t41+2t1)接近1。至于距离r,时延变化可能导致t21+t41-t31减小,放大误差。因此时延误差对距离定位影响较大,而方向角对时延误差不敏感。

在上述一种用于局部放电源定位的极坐标网格搜索方法的基础上,本发明还提供一种用于局部放电源定位的极坐标网格搜索系统,包括:

网格划分模块:利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形的网格,并分别标记为1,2,3…n,网格中心点记为Pk(rkk),k=1,2,3…n;

特征值数组计算模块:特征值数组(A1,A2,…An)的计算公式包括:

其中,c是光速,τi是假设局部放电源位于(r,θ),其放电信号到达第i个传感器与到达第一个传感器的时间差,i是传感器序号,ri是第i个传感器的极坐标距离值,fi(R,θ)是局放源到第i个传感器与到第一个传感器的理论距离差与依照时延计算的距离差的误差值,Ak值为到三个传感器误差值的集合指标作为特征值反应定位误差大小情况;

定位模块:查找特征值数组中最小值Am,最小值Am所在的网格即局放点所在的网格。

判断循环模块:判断最小值Am所在的网格的大小是否符合精度要求,若判断结果为是则作为最终解,若判断结果为否则将最小值Am所在的网格作为下一个待搜索的空间返回再次进行网格划分。

实施例1:局部放电源的定位仿真

首先设置如下参数:

(1)如图3,有四个特高频传感器分别放置于

以及

(2)一个局部放电源可能任意出现在1000个半径为60米的圆形范围内,圆心与特高频传感器中心重合;

(3)算法中所使用的时延是在真实时延的基础上外加一个服从σ~U(-ε00)分布的随机误差,其中ε0=1ns。

在直角坐标系下,为了达到搜索精度为1米的收敛条件,如果对一个网格的搜索作为一次计算,随着搜索条件设置,计算量在675到1800之间可以将超过90%的局放源定位到2度的误差范围内。而在极坐标下,通过设置初始搜索半径为60米,搜索半径缩小倍数为2,初始分割角度为40度,搜索角度缩小倍数为5,可以在计算量为69的情况下精准定位。

实施例2:局部放电源的现场定位

在中国南方电网某220Kv变电站中利用此算法定位局放源。如图4所示,天线阵列布置于中控室屋顶,坐标分布为S1(-7.8,3.75,1.2),S_2(7.8,3.75,1.2),S3(7.8,-3.75,-3.15)以及S_4(7.8,-3.75,-3.15)阵列中心点为o(0,0,0),绝缘子上面的沿面放电作为局部放电源。十六组时延信号组成的数据集合被获取用于算法计算。极坐标系下收敛角度为2度,直角坐标系下收敛边长为2米。设置极坐标系下初始搜索半径为30米,初始搜索角度为40度,半径缩小倍数为2,角度缩小倍数为5,直角坐标系下初始搜索边长为6米,边长缩小倍数为5。结果如图5所示,在悬殊的计算量下,大部分的角度定位精度都能高至2度以内。

本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以,本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构;也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。

Claims (8)

1.一种用于局部放电源定位的极坐标网格搜索方法,其特征在于,包括:
网格划分步骤:利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形的网格,并分别标记为1,2,3…n,网格中心点记为Pk(rkk),k=1,2,3…n;
特征值数组计算步骤:特征值数组(A1,A2,…An)的计算公式包括:
其中,c是光速,τi是假设局部放电源位于(r,θ),其放电信号到达第i个传感器与到达第一个传感器的时间差,i是传感器序号,ri是第i个传感器的极坐标距离值,fi(R,θ)是局放源到第i个传感器与到第一个传感器的理论距离差与依照时延计算的距离差的误差值,Ak值为到三个传感器误差值的集合指标作为特征值反应定位误差大小情况;
定位步骤:查找特征值数组中最小值Am,最小值Am所在的网格即局放点所在的网格。
2.根据权利要求1所述的用于局部放电源定位的极坐标网格搜索方法,其特征在于,还包括:
判断循环步骤:判断最小值Am所在的网格的大小是否符合精度要求,若判断结果为是则作为最终解,若判断结果为否,则将最小值Am所在的网格作为下一个待搜索的空间,并自定义按照一定比例扩大搜索网格或不扩大以原有网格作为下一个搜索空间,返回网格划分步骤。
3.根据权利要求2所述的用于局部放电源定位的极坐标网格搜索方法,其特征在于,在循环时新的r值为原半径除以缩小倍数,新的θ值为原角度除以角度缩小倍数。
4.根据权利要求1所述的用于局部放电源定位的极坐标网格搜索方法,其特征在于,传感器的数量为4个,呈矩形放置。
5.一种用于局部放电源定位的极坐标网格搜索系统,其特征在于,包括:
网格划分模块:利用极坐标系,将搜索空间划分为半径为r,角度为θ,数量为nr的一个个扇形的网格,并分别标记为1,2,3…n,网格中心点记为Pk(rkk),k=1,2,3…n;
特征值数组计算模块:特征值数组(A1,A2,…An)的计算公式包括:
其中,c是光速,τi是假设局部放电源位于(r,θ),其放电信号到达第i个传感器与到达第一个传感器的时间差,i是传感器序号,ri是第i个传感器的极坐标距离值,fi(R,θ)是局放源到第i个传感器与到第一个传感器的理论距离差与依照时延计算的距离差的误差值,Ak值为到三个传感器误差值的集合指标作为特征值反应定位误差大小情况;
定位模块:查找特征值数组中最小值Am,最小值Am所在的网格即局放点所在的网格。
6.根据权利要求5所述的用于局部放电源定位的极坐标网格搜索系统,其特征在于,还包括:
判断循环模块:判断最小值Am所在的网格的大小是否符合精度要求,若判断结果为是则作为最终解,若判断结果为否则将最小值Am所在的网格作为下一个待搜索的空间返回再次进行网格划分。
7.根据权利要求6所述的用于局部放电源定位的极坐标网格搜索系统,其特征在于,在循环时新的r值为原半径除以缩小倍数,新的θ值为原角度除以角度缩小倍数。
8.根据权利要求5所述的用于局部放电源定位的极坐标网格搜索系统,其特征在于,传感器的数量为4个,呈矩形放置。
CN201810247930.5A 2018-03-23 2018-03-23 一种用于局部放电源定位的极坐标网格搜索方法及系统 CN108761279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810247930.5A CN108761279B (zh) 2018-03-23 2018-03-23 一种用于局部放电源定位的极坐标网格搜索方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810247930.5A CN108761279B (zh) 2018-03-23 2018-03-23 一种用于局部放电源定位的极坐标网格搜索方法及系统

Publications (2)

Publication Number Publication Date
CN108761279A CN108761279A (zh) 2018-11-06
CN108761279B true CN108761279B (zh) 2020-06-05

Family

ID=63980187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810247930.5A CN108761279B (zh) 2018-03-23 2018-03-23 一种用于局部放电源定位的极坐标网格搜索方法及系统

Country Status (1)

Country Link
CN (1) CN108761279B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU995111A1 (ru) * 1980-10-08 1983-02-07 Куйбышевский институт инженеров железнодорожного транспорта Нагл дное пособие по математике
EP1227340A1 (fr) * 2001-01-30 2002-07-31 Meteo-France Procédé de détection radar
CN102445641A (zh) * 2011-11-01 2012-05-09 上海交通大学 移动式变电站电力设备局部放电检测装置和定位方法
CN103913680A (zh) * 2014-03-17 2014-07-09 上海交通大学 基于时延误差随机模拟与统计分析的局部放电定位方法
CN104237750A (zh) * 2014-09-05 2014-12-24 中国西电电气股份有限公司 一种gis绝缘缺陷局部放电故障图谱绘制方法
CN105866636A (zh) * 2016-03-23 2016-08-17 国家电网公司 基于时差定位的变电站定位方法
CN107505548A (zh) * 2017-08-29 2017-12-22 华北电力大学(保定) 一种基于柔性阵列传感器的局放超声定位方法
CN107589352A (zh) * 2017-09-05 2018-01-16 浙江群力电气有限公司 一种局部放电源的定位方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH696646A5 (de) * 2003-01-15 2007-08-31 Alstom Technology Ltd Verfahren zur Analyse und/oder Überwachung des Teilentladungsverhaltens eines elektrischen Betriebsmittels.

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU995111A1 (ru) * 1980-10-08 1983-02-07 Куйбышевский институт инженеров железнодорожного транспорта Нагл дное пособие по математике
EP1227340A1 (fr) * 2001-01-30 2002-07-31 Meteo-France Procédé de détection radar
CN102445641A (zh) * 2011-11-01 2012-05-09 上海交通大学 移动式变电站电力设备局部放电检测装置和定位方法
CN103913680A (zh) * 2014-03-17 2014-07-09 上海交通大学 基于时延误差随机模拟与统计分析的局部放电定位方法
CN104237750A (zh) * 2014-09-05 2014-12-24 中国西电电气股份有限公司 一种gis绝缘缺陷局部放电故障图谱绘制方法
CN105866636A (zh) * 2016-03-23 2016-08-17 国家电网公司 基于时差定位的变电站定位方法
CN107505548A (zh) * 2017-08-29 2017-12-22 华北电力大学(保定) 一种基于柔性阵列传感器的局放超声定位方法
CN107589352A (zh) * 2017-09-05 2018-01-16 浙江群力电气有限公司 一种局部放电源的定位方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GIS体外UHF法检测局部放电的研究与应用;钱江;《科技情报开发与经济》;20111231;第21卷(第34期);第198-201页 *
变压器绕组对局部放电特高频定位方法的影响;郑书生 等;《高电压技术》;20130228;第39卷(第2期);第324-329页 *

Also Published As

Publication number Publication date
CN108761279A (zh) 2018-11-06

Similar Documents

Publication Publication Date Title
CN102788926B (zh) 小电流接地系统单相接地故障区段定位方法
CN102305891B (zh) 一种电力系统低频振荡在线监测方法
Viberg et al. Analysis of the combined effects of finite samples and model errors on array processing performance
Wan et al. PD source diagnosis and localization in industrial high-voltage insulation system via multimodal joint sparse representation
CN103558512B (zh) 一种基于矩阵运算的配电网10kV馈线故障定位方法
Demarest et al. FDTD modeling of scatterers in stratified media
CN103176077B (zh) 一种数控成品电路板在环境综合作用下的可靠性快速测评方法
Zhang et al. Improved DOA estimation algorithm for co‐prime linear arrays using root‐MUSIC algorithm
CN102735999B (zh) Gis局部放电在线监测系统及其故障模式识别方法
Hou et al. Localization algorithm for the PD source in substation based on L-shaped antenna array signal processing
EP1611446A2 (en) Systems and methods for time corrected lightning detection
JP2016539333A (ja) ターゲットの位置決め方法及びシステム
US20160011292A1 (en) System and Method for Determining Angle
CN106096184B (zh) 一种基于声强法的变压器多点声源模型的噪声计算方法及系统
CN103064073B (zh) 一种基于频率捷变改变雷达目标特性的方法
CN102866334A (zh) 一种车载式变电站局部放电定位系统及其定位方法
US20180274520A1 (en) Method for evaluating power characteristics of wind turbines, apparatus and storage medium
CN203930058U (zh) 一种合成孔径雷达面目标回波信号模拟器
WO2013059636A1 (en) Time of arrival based wireless positioning system
CN102111876A (zh) 一种选择用于定位的参考标签的方法及装置
CN103954940B (zh) 雷达网基于交叉定位点聚类的集中式与分布式压制干扰鉴别方法
CN106324340B (zh) 一种同步相量和频率测量动态性能的方法
CN106154037B (zh) 一种基于校验的同步相量自适应计算方法
CN103116688A (zh) 用于机载航电系统的多源异类传感器目标航迹关联方法
CN105116435A (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