CN108152687B - 一种利用超声反演的电力变压器局部放电定位方法 - Google Patents

一种利用超声反演的电力变压器局部放电定位方法 Download PDF

Info

Publication number
CN108152687B
CN108152687B CN201711420829.7A CN201711420829A CN108152687B CN 108152687 B CN108152687 B CN 108152687B CN 201711420829 A CN201711420829 A CN 201711420829A CN 108152687 B CN108152687 B CN 108152687B
Authority
CN
China
Prior art keywords
power transformer
partial discharge
sensor array
ultrasonic
array element
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
Application number
CN201711420829.7A
Other languages
English (en)
Other versions
CN108152687A (zh
Inventor
赵荣普
赵威
王欣
万子逸
杨进
韦瑞峰
陈永琴
胡昌斌
陈井锐
李梦滔
张国志
吴辰阳
陈欣
段永生
周涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming Power Supply Bureau of Yunnan Power Grid Co Ltd
Original Assignee
Kunming Power Supply Bureau of Yunnan Power Grid Co Ltd
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 Kunming Power Supply Bureau of Yunnan Power Grid Co Ltd filed Critical Kunming Power Supply Bureau of Yunnan Power Grid Co Ltd
Priority to CN201711420829.7A priority Critical patent/CN108152687B/zh
Publication of CN108152687A publication Critical patent/CN108152687A/zh
Application granted granted Critical
Publication of CN108152687B publication Critical patent/CN108152687B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/1209Testing 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 using acoustic measurements

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Testing Relating To Insulation (AREA)

Abstract

本发明公开了一种利用时间反演的电力变压器局部放电定位方法,首先对测量的电力变压器局部放电数据进行时间反演处理,然后依据时间反演超声波信号确定电力变压器局部放电超声的波达方向,再沿局部放电超声的波达方向搜索确定电力变压器局部放电位置,这样不仅可以极大减小搜索范围、大大缩短电力变压器局部放电位置搜索所需时间,可以提高电力变压器局部放电位置的测量效率和定位精度。

Description

一种利用超声反演的电力变压器局部放电定位方法
技术领域
本发明属于电通信诊断技术领域,涉及电力变压器运行中的局部放电诊断技术,尤其涉及一种利用超声反演的电力变压器局部电定位方法。
背景技术
电力变压器是用来将某一数值的交流电压(电流)变成频率相同的另一种或几种数值不同的电压(电流)的设备,其是保证电力系统正常可靠运行的重要组成部分。
随着经济的发展、人民生活水平的提高,用电量不断增加,需要大量的电力变压器投入电力系统的运行。但是,由于制作工艺以及长期处在不同的运行工况下的影响,电力变压器往往会在局部区域产生绝缘缺陷,如果不对电力变压器的局部缺陷进行处理,电力变压器局部缺陷在强电场的作用下就会迅速发展,产生局部放电现象,最终导致电力变压器绝缘的失效,从而使电力变压器的维修更加困难,甚至直接报废,导致大量的人力、财力浪费。此外,由于电力变压器结构复杂,如果无法预先知道电力变压器局部放电的位置,将会在电压变压器检修时浪费大量的时间进行局部放电位置的定位,从而大大延长检修周期,降低供电可靠性。
现有比较成熟的电力变压器局部放电定位方法,如(1)申请号为200910175252.7的专利申请文件公开的变压器局部放电的相控超声定位方法及系统,采用置于变压器上两个不同位置的相控超声传感器接收局放源发出的超声宽频信号,然后应用宽带聚焦算法、窄带测向算法和交叉算法对接收的超声宽频信号进行处理确定局部放电的具体位置,由于该方法需要多次测量利用不同的方向线进行定位,增加了现场人员的工作量;(2)申请号为201310060568.8的专利申请文件公开的变压器局部放电超声源重建的定位方法及装置,提供了一种利用时间反转技术结合变压器传播模型计算局部放电源位置及估计放电量大小的方法,由记录的声压信号计算反演声场,实现声场的聚集和声源的重建,再根据聚焦声场的位置和强度确定放电点位置和估算放电量大小,但该方法需要准确获取电力变压器的传播模型,由于电力变压器结构复杂,因此难以准确地建立电力变压器的传播模型,因此这种方法难以实现对电力变压器局部放电源位置的精准定位,同时,对于不同的电力变压器,其传播模型也相差很大,这种方法缺乏很好的通用性。因此开发有效的电力变压器局部放电定位技术意义重大。
发明内容
针对目前缺少现电力变压器局部放电位置准确有效的定位技术,本发明旨在提供一种利用超声反演的电力变压器局部电定位方法,以实现对电力变压器局部放电进行定向和定位。
本发明的思路,首先利用时间反演对电力变压器局部放电信号进行时间反演处理,然后在笛卡尔空间坐标系下,根据不同方位角和不同仰角θ下各个传感器阵元相对于第一个传感器阵元的时间差ti对时间反演超声波信号si(T-t)进行时延处理,并依据时延处理后的超声波信号得到辐射能量最大的电力变压器局部放电沿方向实现电力变压器局部放电超声信号的定向;再在确定的方向上构建搜索点,根据搜索点到各个传感器阵元的距离hi得到搜索点到各个传感器阵元相对于到第一个传感器阵元的时间差mi,并利用该时间差mi再次对时间反演超声波信号si(T-t)进行时延处理,并依据时延处理后的超声波信号得到方向上能量极值点,能量极值点对应的搜索点即为电力变压器的局部放电位置。
本发明提供的利用超声反演的电力变压器局部放电定位方法,其技术方案的构成包括以下步骤:
(1)局部放电数据测量
采用设置于电力变压器外壁表面的超声传感器阵列检测电力变压器发生局部放电时产生的超声波信号si(t),以超声传感器阵列几何中心为原点,建立笛卡尔空间坐标系,并获取超声传感器阵列各个传感器空间位置Pi(xi,yi,zi),其中i为超声传感器阵列中传感器阵元序号,t检测时间,x、y、z为各传感器的横坐标、纵坐标和高度坐标;
(2)局部放电定向,包括以下分步骤:
(21)将超声传感器阵列各个传感器阵元检测得到的超声波信号si(t)进行时间反演,得到时间反演超声波信号si(T-t);
(22)将建立的笛卡尔空间坐标系进行相位划分,其中将方位角在0到2π范围内划分为N1等份,将仰角θ在0到π/2范围内划分为N2等份,得到不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di
(23)根据不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di计算各个传感器阵元相对于第一个传感器阵元的时间差ti
其中v为超声在变压器油中的传播速度;
(24)根据得到的各个传感器阵元相对于第一个传感器阵元的时间差ti,对步骤(21)所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si′(t)为:
其中,FFT为快速傅里叶变换,IFFT为快速傅里叶反变换,e为指数函数,ω为角频率;
(25)根据步骤(24)所得时延处理的超声波信号si′(t)得到电力变压器局部放电沿方向的辐射能量:
其中,n为超声传感器阵列阵元的个数;
遍历N1等份方位角和N2等份仰角θ,得到一个N1×N2的能量矩阵;选取能量矩阵中最大点对应的方位角和仰角θmax,即为电力变压器局部放电超声的波达方向,即实现了电力变压器局部放电超声信号的定向;
(3)局部放电定位,包括以下分步骤:
(31)根据得到的方位角和仰角θmax构建空间搜索点坐标,其为由坐标原点开始,方位角为仰角为θmax方向的直线上的点o(x,y,z);
(32)搜索点到各个传感器阵元的距离hi
(33)根据搜索点到各个传感器阵元的距离hi计算搜索点至各个传感器阵元相对于搜索点至第一个传感器阵元的时间差mi
(34)根据得到的搜索点至各个传感器阵元相对于第一个传感器阵元的时间差mi,对步骤(21)所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si″(t)为:
(35)根据步骤(34)所得时延处理的超声波信号si″(t)得到电力变压器局部放电在搜索点o(x,y,z)的能量:
遍历由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点,并计算各个搜索点的能量,选取能量最大对应的搜索点omax(xmax,ymax,zmax)即为电力变压器的局部放电位置。
上述利用超声反演的电力变压器局部放电定位方法,所述超声传感器阵列是由设置在电力变压器外壁表面的若干超声传感器组成,超声传感器在电力变压器外壁表面的位置为随机设计。超声传感器阵列可以只设置在一侧的外壁表面,其排列形式可以是随机的,也可以是按照一定规律排列成规则阵列,例如由若干超声传感器在电力变压器外壁表面排列构成l×k矩阵,l、k为非零正整数,或者由若干超声传感器在电力变压器外壁表面排列构成圆形阵列。
上述利用超声反演的电力变压器局部放电定位方法,所述步骤(2)的目的在于根据电力变压器局部放电信号在各个方向的辐射能量,找到辐射能量最大的方向。为此,需要先将建立的笛卡尔空间坐标系按照方位角和仰角θ进行相位划分,将各个传感器阵元相对于坐标原点的距离采用方位角和仰角θ进行表示;在优选的实施方式中,将方位角在0到2π范围内划分为N1等份,100≤N1≤2000;将仰角θ在0到π/2范围内划分为N2等份,100≤N2≤500。然后依据不同方位角和不同仰角θ下各个传感器阵元相对于第一个传感器阵元的时间差ti,对步骤(21)所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si′(t),再根据所得时延处理的超声波信号si′(t)得到电力变压器局部放电沿方向的辐射能量,能量最大点对应的方位角和仰角θmax,即为电力变压器局部放电超声的波达方向,即实现了电力变压器局部放电超声信号的定向。
上述利用超声反演的电力变压器局部放电定位方法,所述步骤(3)的目的在于重新沿着方位角和仰角θmax方向搜索,其辐射能量最大的搜索点即为电力变压器局部放电位置。为此,需要先根据得到的方位角和仰角θmax构建空间搜索坐标,并获得搜索点到各个传感器阵元的距离。然后计算各个传感器阵元相对于第一个传感器阵元的时间差mi,重新对步骤(21)所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si″(t),再根据所得时延处理的超声波信号si″(t)得到电力变压器局部放电在搜索点的辐射能量,能量最大点对应的搜索点,即为电力变压器局部放电位置。所述步骤(35)中,对由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点进行遍历时可以采用等间距方式选取搜索点位置;或者先采用第一间距遍历由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点,得到最大能量对应的搜索点,然后采用小于第一间距的第二间距遍历最大能量对应的搜索点前后0.1米范围内的搜索点,其能量极大值点即为电力变压器的局部放电位置。
与现有技术相比,本发明提供的利用超声反演的电力变压器局部电定位方法,具有如下十分突出的优点和有益技术效果:
1.本发明基于先定向后定位的发明思路,首先利用求解电力变压器局部放电超声信号能量极值的方法确定电力变压器局部放电超声的波达方向,沿局部放电超声的波达方向再次利用求解电力变压器局部放电超声信号能量极值的方法确定电力变压器局部放电位置,这样不仅可以极大减小搜索范围、大大缩短电力变压器局部放电位置搜索所需时间,而且无需利用多个方位线交点确定局部放电源位置,减小了方位线不相交时存在的误差,从而有效提高电力变压器局部放电位置的测量效率和定位精度;
2.本发明对检测到的超声波信号进行时间反演处理,能够有效弥补局部放电点到各个传感器的时间延迟;
3.本发明利用超声传感器阵列各个传感器阵元信号的相互叠加,有效提高了低信噪比条件下的定向和定位精度;
4.本发明由于需要对整个电力变压器进行搜索定位,因此受量化误差及采样率影响较小,在量化误差较大或采样率较低时仍能实现电力变压器局部放电的精确定位;
5.本发明相比于传统的电力变压器局部放电定位方法,仅需进行一次超声波信号测试即可实现电力变压器局部放电定向及定位分析,极大降低了电力变压器局部放电定位难度,且具有很好的通用性,适于在本领域内推广使用。
附图说明
图1为本发明利用超声反演的电力变压器局部放电定位方法流程示意图。
图2为本发明利用超声反演的电力变压器局部放电定位方法中局部放电定向步骤的流程示意图。
图3为本发明利用超声反演的电力变压器局部放电定位方法中局部放电定位步骤的流程示意图。
图4为本发明实施例1利用超声反演的电力变压器局部放电定位中得到的电力变压器局部放电定向谱图。
图5为本发明实施例1利用超声反演的电力变压器局部放电定位中得到的电力变压器局部放电定位图。
图6为本发明实施例2利用超声反演的电力变压器局部放电定位中得到的电力变压器局部放电定向谱图。
图7为本发明实施例2利用超声反演的电力变压器局部放电定位中得到的电力变压器局部放电定位图。
具体实施方式
以下将结合附图给出本发明实施例,并通过实施例对本发明的技术方案进行进一步的清楚、完整说明。显然,所述实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明内容,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例,都属于本发明所保护的范围。
实施例1
本实施例针对的电力变压器物理尺寸为4m×4m×4m,局部放电点位于(1m,1m,2m)处,利用超声反演的电力变压器局部电定位过程如图1至图3所示,包括以下步骤:
S1、局部放电数据测量
本实施例采用的超声传感器阵列为由16个传感器组成的4×4矩阵,相邻两个超声传感器之间的距离为5mm,将超声传感器阵列表贴于被测电力变压器某一侧壁的外表面中心处,检测电力变压器发生局部放电时产生的超声波信号si(t),以超声传感器阵列几何中心为原点,建立笛卡尔空间坐标系,并获取超声传感器阵列各个传感器空间位置Pi(xi,yi,zi),其中i为超声传感器阵列中传感器阵元序号,i为1-16,t检测时间,x、y、z为各个超声传感器的横坐标、纵坐标和高度坐标。
S2、局部放电定向,包括以下分步骤:
S21、将超声传感器阵列各个传感器阵元检测得到的超声波信号si(t)进行时间反演t→T-t处理,得到时间反演超声波信号si(T-t);
S22、将建立的笛卡尔空间坐标系进行相位划分,其中将方位角在0到2π范围内划分为201等份,将仰角θ在0到π/2范围内划分为251等份,得到不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di
S23、将得到的不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di与第一个传感器阵元相对于坐标原点的距离d1相减,得到各个传感器阵元相对于第一个传感器阵元的距离差,进而得到各个传感器阵元相对于第一个传感器阵元的时间差ti
其中v为超声在变压器油中的传播速度,取值为1500m/s;
S24、由时间反演性质可知,时间反演处理后的信号能够弥补信号的时延时间,根据得到的各个传感器阵元相对于第一个传感器阵元的时间差ti,对步骤S21所得经时间反演处理后的超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si′(t)为:
其中,FFT为快速傅里叶变换,IFFT为快速傅里叶反变换,e为指数函数,ω为角频率;
S25、根据步骤S24所得时延处理的超声波信号si′(t)得到电力变压器局部放电沿方向的辐射能量:
其中,n为超声传感器阵列阵元的个数,n=16;
遍历201等份方位角和251等份仰角θ,得到一个201×251的能量矩阵,由时间反演性质可知,当遍历的方位角和仰角θ恰好为实际局部放电的方向时,此时阵列信号将在该方向上实现能量聚焦,因此选取能量矩阵中最大点的方位角和仰角θmax即为电力变压器局部放电超声的波达方向,即实现了电力变压器局部放电超声信号的定向。
由图4可知,当方位角为45度、仰角θmax为54.72度时能量达到最大值,该方位恰好为真实的电力变压器局部放电方位。
S3、局部放电定位,包括以下分步骤:
S31、根据得到的方位角和仰角θmax构建空间搜索点坐标,其为由坐标原点开始,方位角为仰角为θmax方向的直线上的点o(x,y,z);
S32、搜索点到各个传感器阵元的距离hi
S33、将得到的搜索点与各个传感器阵元间的距离hi与搜索点与第一个传感器阵元间的距离h1相减,得到搜索点为o(x,y,z)时,搜索点至各个传感器阵元相对于第一个传感器阵元的距离差,进而得到搜索点至各个传感器阵元相对于第一个传感器阵元的时间差mi
S34、根据得到的搜索点至各个传感器阵元相对于第一个传感器阵元的时间差mi,对步骤S21所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si″(t)为:
S35、根据步骤S34所得时延处理的超声波信号si″(t)得到电力变压器局部放电在搜索点o(x,y,z)的能量:
遍历由坐标原点开始,方位角为为45度、仰角为θmax为54.72度方向的直线上的所有点,并计算各个遍历点的能量。由时间反演性质可知,当遍历位置恰好为电力变压器实际局部放电位置时,能量将达到最大;选取能量最大的遍历点omax(xmax,ymax,zmax)即为电力变压器的局部放电位置;本实施例采用第一间距0.05米,在由坐标原点开始,方位角为45度、仰角θmax为54.72度的直线选取搜索点,重复步骤S32-S35,得到电力变压器位置-能量谱图;进一步针对获得的电力变压器局部放电位置-能量谱图中极大位置前后0.1米范围内的线段,采用第二间距0.001米,在所得线段上选取搜索点,重复步骤S32-S35,得到更加详细的电力变压器局部位置-能量谱图,则局部电力变压器局部位置-能量谱图中的能量极大值点即为电力变压器的局部放电位置。
由图5可知,当遍历点为(0.998m,0.998m,1.995m)时能量达到最大,该位置与实际局部放电点的位置非常接近。
实施例2
本实施例针对的电力变压器物理尺寸为4m×4m×4m,局部放电点位于(-1m,-0.5m,1m)处,利用超声反演的电力变压器局部电定位过程如图1至图3所示,包括以下步骤:
S1、局部放电数据测量
本实施例采用的超声传感器阵列为由16个传感器组成的4×4矩阵,相邻两个超声传感器之间的距离为5mm,将超声传感器阵列表贴于被测电力变压器某一侧壁的外表面中心处,检测电力变压器发生局部放电时产生的超声波信号si(t),以超声传感器阵列几何中心为原点,建立笛卡尔空间坐标系,并获取超声传感器阵列各个传感器空间位置Pi(xi,yi,zi),其中i为超声传感器阵列中传感器阵元序号,i为1-16,t检测时间,x、y、z为各个超声传感器的横坐标、纵坐标和高度坐标。
S2、局部放电定向,包括以下分步骤:
S21、将超声传感器阵列各个传感器阵元检测得到的超声波信号si(t)进行时间反演t→T-t处理,得到时间反演超声波信号si(T-t);
S22、将建立的笛卡尔空间坐标系进行相位划分,其中将方位角在0到2π范围内划分为1001等份,将仰角θ在0到π/2范围内划分为251等份,得到不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di
S23、将得到的不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di与第一个传感器阵元相对于坐标原点的距离d1相减,得到各个传感器阵元相对于第一个传感器阵元的距离差,进而得到各个传感器阵元相对于第一个传感器阵元的时间差ti
其中v为超声在变压器油中的传播速度,取值为1500m/s;
S24、由时间反演性质可知,时间反演处理后的信号能够弥补信号的时延时间,根据得到的各个传感器阵元相对于第一个传感器阵元的时间差ti,对步骤S21所得经时间反演处理后的超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si′(t)为:
其中,FFT为快速傅里叶变换,IFFT为快速傅里叶反变换,e为指数函数,ω为角频率;
S25、根据步骤S24所得时延处理的超声波信号si′(t)得到电力变压器局部放电沿方向的辐射能量:
其中,n为超声传感器阵列阵元的个数,n=16;
遍历1001等份方位角和251等份仰角θ,得到一个1001×251的能量矩阵,由时间反演性质可知,当遍历的方位角和仰角θ恰好为实际局部放电的方向时,此时阵列信号将在该方向上实现能量聚焦,因此选取能量矩阵中最大点的方位角和仰角θmax即为电力变压器局部放电超声的波达方向,即实现了电力变压器局部放电超声信号的定向。
由图6可知,当方位角为206.64度、仰角θmax为41.76度时能量达到最大值,该方位恰好为真实的电力变压器局部放电方位。
S3、局部放电定位,包括以下分步骤:
S31、根据得到的方位角为206.64和仰角θmax为41.76度构建空间搜索点坐标,其为由坐标原点开始,方位角为仰角为θmax方向的直线上的点o(x,y,z);
S32、搜索点到各个传感器阵元的距离hi
S33、将得到的搜索点与各个传感器阵元间的距离hi与搜索点与第一个传感器阵元间的距离h1相减,得到搜索点为o(x,y,z)时,搜索点至各个传感器阵元相对于搜索点至第一个传感器阵元的距离差,进而得到搜索点至各个传感器阵元相对于搜索点至第一个传感器阵元的时间差mi
S34、根据得到的搜索点至各个传感器阵元相对于搜索点至第一个传感器阵元的时间差mi,对步骤S21所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si″(t)为:
S35、根据步骤S34所得时延处理的超声波信号si″(t)得到电力变压器局部放电在搜索点o(x,y,z)的能量:
遍历由坐标原点开始,方位角为为206.64度、仰角为θmax为41.76度方向的直线上的所有点,并计算各个遍历点的能量。由时间反演性质可知,当遍历位置恰好为电力变压器实际局部放电位置时,能量将达到最大;选取能量最大的遍历点omax(xmax,ymax,zmax)即为电力变压器的局部放电位置;本实施例采用第一间距0.05米,在由坐标原点开始,方位角为206.64度、仰角θmax为41.76度的直线选取搜索点,重复步骤S32-S35,得到电力变压器位置-能量谱图;进一步针对获得的电力变压器局部放电位置-能量谱图中极大位置前后0.1米范围内的线段,采用第二间距0.001米,在所得线段上选取搜索点,重复步骤S32-S35,得到更加详细的电力变压器局部位置-能量谱图,则局部电力变压器局部位置-能量谱图中的能量极大值点即为电力变压器的局部放电位置。
由图7可知,当遍历点为(-0.998m,-0.501m,0.997m)时能量达到最大,该位置与实际局部放电点的位置非常接近。

Claims (7)

1.一种利用超声反演的电力变压器局部放电定位方法,其特征在于包括以下步骤:
(1)局部放电数据测量
采用设置于电力变压器外壁表面的超声传感器阵列检测电力变压器发生局部放电时产生的超声波信号si(t),以超声传感器阵列几何中心为原点,建立笛卡尔空间坐标系,并获取超声传感器阵列各个传感器空间位置Pi(xi,yi,zi),其中i为超声传感器阵列中传感器阵元序号,t为检测时间,x、y、z为各传感器的横坐标、纵坐标和高度坐标;
(2)局部放电定向,包括以下分步骤:
(21)将超声传感器阵列各个传感器阵元检测得到的超声波信号si(t)进行时间反演,得到时间反演超声波信号si(T-t);
(22)将建立的笛卡尔空间坐标系进行相位划分,其中将方位角在0到2π范围内划分为N1等份,将仰角θ在0到π/2范围内划分为N2等份,得到不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di
(23)根据不同方位角和不同仰角θ下各个传感器阵元相对于坐标原点的距离di计算各个传感器阵元相对于第一个传感器阵元的时间差ti
其中v为超声在变压器油中的传播速度;
(24)根据得到的各个传感器阵元相对于第一个传感器阵元的时间差ti,对步骤(21)所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si′(t)为:
其中,FFT为快速傅里叶变换,IFFT为快速傅里叶反变换,e为指数函数,ω为角频率;
(25)根据步骤(24)所得时延处理的超声波信号si′(t)得到电力变压器局部放电沿方向的辐射能量:
其中,n为超声传感器阵列阵元的个数;
遍历N1等份方位角和N2等份仰角θ,得到一个N1×N2的能量矩阵;选取能量矩阵中最大点对应的方位角和仰角θmax,即为电力变压器局部放电超声的波达方向,即实现了电力变压器局部放电超声信号的定向;
(3)局部放电定位,包括以下分步骤:
(31)根据得到的方位角和仰角θmax构建空间搜索点坐标,其为由坐标原点开始,方位角为仰角为θmax方向的直线上的点o(x,y,z);
(32)搜索点到各个传感器阵元的距离hi
(33)根据搜索点到各个传感器阵元的距离hi计算从搜索点到各个传感器阵元与从搜索点到第一个传感器阵元的时间差mi
(34)根据得到的从搜索点到各个传感器阵元相对于到第一个传感器阵元的时间差mi,对步骤(21)所得时间反演超声波信号si(T-t)进行时延处理,得到时延处理后的超声波信号si″(t)为:
(35)根据步骤(34)所得时延处理的超声波信号si″(t)得到电力变压器局部放电在搜索点o(x,y,z)的能量:
遍历由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点,并计算各个搜索点的能量,选取能量最大对应的搜索点omax(xmax,ymax,zmax)即为电力变压器的局部放电位置。
2.根据权利要求1所述利用超声反演的电力变压器局部放电定位方法,其特征在于所述超声传感器阵列是由设置在电力变压器外壁表面的若干超声传感器组成,超声传感器在电力变压器外壁表面的位置为随机设计。
3.根据权利要求2所述利用超声反演的电力变压器局部放电定位方法,其特征在于若干超声传感器在电力变压器外壁表面排列构成l×k矩阵,l、k为非零正整数。
4.根据权利要求2所述利用超声反演的电力变压器局部放电定位方法,其特征在于若干超声传感器在电力变压器外壁表面排列构成圆形阵列。
5.根据权利要求1至4任一权利要求所述利用超声反演的电力变压器局部放电定位方法,其特征在于所述步骤(22)中,100≤N1≤2000,100≤N2≤500。
6.根据权利要求5所述利用超声反演的电力变压器局部放电定位方法,其特征在于所述步骤(35)中,对由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点进行遍历时采用等间距方式选取搜索点位置。
7.根据权利要求5所述利用超声反演的电力变压器局部放电定位方法,其特征在于所述步骤(35)中,对由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点进行遍历时,采用第一间距遍历由坐标原点开始,方位角为仰角为θmax方向的直线上的搜索点,得到最大能量对应的搜索点;然后采用小于第一间距的第二间距遍历最大能量对应的搜索点前后0.1米范围内的搜索点,其能量极大值点即为电力变压器的局部放电位置。
CN201711420829.7A 2017-12-25 2017-12-25 一种利用超声反演的电力变压器局部放电定位方法 Active CN108152687B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711420829.7A CN108152687B (zh) 2017-12-25 2017-12-25 一种利用超声反演的电力变压器局部放电定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711420829.7A CN108152687B (zh) 2017-12-25 2017-12-25 一种利用超声反演的电力变压器局部放电定位方法

Publications (2)

Publication Number Publication Date
CN108152687A CN108152687A (zh) 2018-06-12
CN108152687B true CN108152687B (zh) 2019-12-31

Family

ID=62465807

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711420829.7A Active CN108152687B (zh) 2017-12-25 2017-12-25 一种利用超声反演的电力变压器局部放电定位方法

Country Status (1)

Country Link
CN (1) CN108152687B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021152377A1 (en) * 2020-01-29 2021-08-05 Ecole Polytechnique Federale De Lausanne (Epfl) Partial discharge localization using time reversal: application to power transformers and gas-insulated substations

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308280B (zh) * 2019-12-11 2022-02-01 云南电网有限责任公司临沧供电局 非接触式超声波检测局放噪声与放电的方法
CN112083297B (zh) * 2020-08-20 2021-08-13 西安交通大学 一种基于声压的抗干扰的变压器局部放电定位方法及装置
CN112466627B (zh) * 2020-11-11 2021-12-24 西安西电变压器有限责任公司 一种变压器油箱的设计方法、系统及变压器油箱
CN117388644B (zh) * 2023-06-27 2024-04-16 国网宁夏电力有限公司电力科学研究院 一种电力设备局部放电定位方法与系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6340890B1 (en) * 1994-12-05 2002-01-22 Abb Research, Ltd. Method and device for locating partial discharges in an electric high-voltage apparatus
CN101702000A (zh) * 2009-11-26 2010-05-05 华北电力大学(保定) 变压器局部放电的相控超声定位方法及系统
CN102707208A (zh) * 2012-06-08 2012-10-03 华北电力大学(保定) 变压器局部放电定位系统及其定位方法
CN103149513A (zh) * 2013-02-26 2013-06-12 中国科学院电工研究所 一种变压器局部放电超声源重建的定位方法及装置
CN104155586A (zh) * 2014-08-21 2014-11-19 华北电力大学(保定) 变压器局部放电超声矢量阵列定位装置及其制造方法
CN104880651A (zh) * 2015-05-28 2015-09-02 国网山西省电力公司电力科学研究院 一种变压器内部局部放电测定装置及其定位方法
CN105974281A (zh) * 2016-05-09 2016-09-28 中国电力科学研究院 一种用于变压器局部放电定位的检测系统和检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6340890B1 (en) * 1994-12-05 2002-01-22 Abb Research, Ltd. Method and device for locating partial discharges in an electric high-voltage apparatus
CN101702000A (zh) * 2009-11-26 2010-05-05 华北电力大学(保定) 变压器局部放电的相控超声定位方法及系统
CN102707208A (zh) * 2012-06-08 2012-10-03 华北电力大学(保定) 变压器局部放电定位系统及其定位方法
CN103149513A (zh) * 2013-02-26 2013-06-12 中国科学院电工研究所 一种变压器局部放电超声源重建的定位方法及装置
CN104155586A (zh) * 2014-08-21 2014-11-19 华北电力大学(保定) 变压器局部放电超声矢量阵列定位装置及其制造方法
CN104880651A (zh) * 2015-05-28 2015-09-02 国网山西省电力公司电力科学研究院 一种变压器内部局部放电测定装置及其定位方法
CN105974281A (zh) * 2016-05-09 2016-09-28 中国电力科学研究院 一种用于变压器局部放电定位的检测系统和检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PD positioning method and experimental study based on circular-shaped ultrasonic array sensors;Li Yanqing 等;《2014 International Conference on Power System Technology》;20141022;第1-6页 *
局部放电超声阵列传感器的声学性能评价及其稀疏设计;程述一;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20141215;C042-52 *
电力变压器智能化改造技术应用研究;赵荣普;《云南电业》;20120331(第3期);第37-38页 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021152377A1 (en) * 2020-01-29 2021-08-05 Ecole Polytechnique Federale De Lausanne (Epfl) Partial discharge localization using time reversal: application to power transformers and gas-insulated substations

Also Published As

Publication number Publication date
CN108152687A (zh) 2018-06-12

Similar Documents

Publication Publication Date Title
CN108152687B (zh) 一种利用超声反演的电力变压器局部放电定位方法
CN108873083B (zh) 一种人工场源频率域电磁视电阻率测量方法
Zhang et al. Efficient immersion imaging of components with nonplanar surfaces
CN108169731A (zh) 基于单近场校正源的拖曳线列阵阵形估计方法和装置
CN105866636B (zh) 基于时差定位的变电站定位方法
CN109696480B (zh) 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法
CN107748129A (zh) 一种用于多孔介质绝缘板材渗水的检测系统及其检测方法
CN107505548A (zh) 一种基于柔性阵列传感器的局放超声定位方法
CN108169511A (zh) 三维空间来风的风速测量系统及方法
CN108169327A (zh) 基于激励波束成型和加权图像融合的music腐蚀监测方法
CN109374748A (zh) 声发射源时反定位成像方法
CN109188534B (zh) 一种基于主动电场原理的水下金属形状探测方法及装置
CN103529428A (zh) 一种基于空间十元阵的点声源被动声定位方法
CN101865361B (zh) 双基阵数据融合电站锅炉四管泄漏被动声测定位方法
Liu et al. Damage imaging of a U-shaped boom using improved periodic permanent magnet electromagnetic acoustic transducers
CN103869096B (zh) 一种超声波风速仪量程拓宽方法
CN111239560B (zh) 基于多传感器的局部放电定位方法
CN113514833A (zh) 一种基于海浪图像的海面任意点波向的反演方法
CN110440896B (zh) 一种超声波测量系统及测量方法
CN112114046A (zh) 一种基于时间反转的铝板损伤混合成像定位方法
CN108508444A (zh) 一种采动影响区输电线路弧垂计算方法及系统
Mengjie et al. Study on gas leakage localization method based on ultrasonic sensor area array
CN113126029B (zh) 适用于深海可靠声路径环境的多传感器脉冲声源定位方法
CN204719179U (zh) 一种变压器内部局部放电测定装置
CN101782644A (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