CN112067953B - 一种提升变电站局部放电定向精度的方法及系统 - Google Patents

一种提升变电站局部放电定向精度的方法及系统 Download PDF

Info

Publication number
CN112067953B
CN112067953B CN202010944714.3A CN202010944714A CN112067953B CN 112067953 B CN112067953 B CN 112067953B CN 202010944714 A CN202010944714 A CN 202010944714A CN 112067953 B CN112067953 B CN 112067953B
Authority
CN
China
Prior art keywords
partial discharge
sensor
discharge signal
signal
representing
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
CN202010944714.3A
Other languages
English (en)
Other versions
CN112067953A (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202010944714.3A priority Critical patent/CN112067953B/zh
Publication of CN112067953A publication Critical patent/CN112067953A/zh
Application granted granted Critical
Publication of CN112067953B publication Critical patent/CN112067953B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Testing Relating To Insulation (AREA)

Abstract

本发明公开了一种提升变电站局部放电定向精度的方法,其包括步骤:(1)采用传感器阵列接收局部放电信号,并基于接收的局部放电信号构建局部放电信号的响应功率;(2)基于局部放电信号的响应功率构建确定局部放电信号的方向角θ的第一模型;(3)对局部放电信号阵列进行虚拟扩展,并基于第一模型得到虚拟扩展后的局部放电信号的定向模型,并基于定向模型对局部放电进行精确定向。此外,本发明还公开了一种提升变电站局部放电定向精度的系统。

Description

一种提升变电站局部放电定向精度的方法及系统
技术领域
本发明涉及一种提升定向精度的方法及系统,尤其涉及一种提升局部放电定向精度的方法及系统。
背景技术
众所周知,电力设备是电力系统的重要基础部分。电力设备的安全可靠运行对电力系统的稳定性和安全性具有重要意义。
在电力系统中,基于声波信号方向角(Direction of arrival,DOA)能够有效估计电力设备在线状态,实现对电力系统中电力设备的监测,其不仅可以为电力设备缺陷提供预警,而且还能进一步防止的故障扩散。
近年来,基于传感器阵列的局部放电定向方法受到了广泛的关注,众多研究人员开始针对该方法进行研究。该局部放电定向方法的性能与所用声传感器阵列的硬件参数密切相关,即阵列孔径,传感器数量,采样频率等。但是,通过硬件方法提高估计性能将始终导致额外的硬件和投资成本,这可能会限制其大规模的实际应用。
基于此,为了解决现有技术中存在的缺陷,期望获得一种提升变电站局部放电定向精度的方法,该方法无需硬件升级即可提高传感器阵列的性能,其可以有效提高在变电站中局部放电的定向精度。
发明内容
本发明的目的之一是提供一种提升变电站局部放电定向精度的方法,该提升变电站局部放电定向精度的方法无需硬件升级即可提高传感器阵列的性能,其通过将四阶累积量替换原始算法中二阶累积量,并利用最大值搜索求取局部放电方向角,从而有效提高在变电站现场中局部放电的定向精度。
根据上述发明目的,本发明提出一种提升变电站局部放电定向精度的方法,其包括步骤:
(1)采用传感器阵列接收局部放电信号,并基于接收的局部放电信号构建局部放电信号的响应功率;
(2)基于局部放电信号的响应功率构建确定局部放电信号的方向角θ的第一模型;
(3)对局部放电信号阵列进行虚拟扩展,并基于第一模型得到虚拟扩展后的局部放电信号的定向模型,并基于定向模型对局部放电进行精确定向。
在本发明所述的技术方案中,本发明提出了一种提升变电站局部放电定向精度的方法,该方法无需硬件升级即可提高传感器阵列的性能,其通过将四阶累积量替换原始算法中二阶累积量,并利用最大值搜索求取局部放电方向角,从而有效提高在变电站现场中局部放电的定向精度。
在本发明所述的提升变电站局部放电定向精度的方法中,设计的主要思想是将原来的信号低阶累积量替换为由泰勒展开建立的高阶累积量。需要说明的是,与低阶累积量相比,采用本发明所述方法所建立的高阶累积量包含在低阶累积量中被省略的高阶项,因此它们具有有关信号源的更详细的方向信息,可以提高声波信号方向角(DOA)估计性能。
进一步地,在本发明所述的提升变电站局部放电定向精度的方法中,在所述步骤(1)中:
将传感器阵列中的第m个传感器接收的局部放电信号记为ym(t),对其进行傅里叶变换后得为Ym(ω);
对局部放电信号Ym(ω)进行宽带波束形成处理,并将其频段划分为L个频率段,将每个频率段中各个传感器接收的局部放电信号加在一起以产生单个聚焦信号Z(ωl,θ):
Figure BDA0002674863220000021
其中,l=1,2,…,L;Yml)表示局部放电信号Ym(ω)划分出的第l个频率段;
Figure BDA0002674863220000022
表示与每个局部放电信号的第l个频率段对应的延时元件,τm表示第m个传感器信号与参考传感器的时延差值;Hml)表示第m个传感器对应的滤波器;M表示传感器阵列中的传感器的数量,m表示第m个传感器;
基于下式建立局部放电信号响应功率:
Figure BDA0002674863220000031
其中Z(ωl,θ)表示局部放电的单个聚焦信号;Z*l,θ)表示Z(ωl,θ)的共轭转置;
Figure BDA0002674863220000032
Figure BDA0002674863220000033
表示第n个传感器对应的滤波器的共轭转置;Wmnl)=E(Yml)·Yn *l)),其中E(·)表示数学期望,Yn *l)表示第n个传感器接收的局部放电信号傅里叶变换的共轭转置;第m个传感器与第n个传感器的组合滤波器Ψmn选为相变函数
Figure BDA0002674863220000034
Figure BDA0002674863220000035
表示自然底数e为底数的指数函数,j是复数的标识符号且有j2=-1;τn表示第n个传感器信号与参考传感器的时延差值。
进一步地,在本发明所述的提升变电站局部放电定向精度的方法中,在所述步骤(2)中,所述第一模型为:
θs=argmaxθP(θ)
其中,P(θ)表示局部放电信号的响应功率,argmax()表示函数取最大值时对应的θ。
进一步地,在本发明所述的提升变电站局部放电定向精度的方法中,在所述步骤(3)中,所述定向模型被构建为:
Figure BDA0002674863220000036
其中,
Figure BDA0002674863220000037
Yml)表示经过傅里叶变换后的第m个传感器接收的局部放电信号Ym(ω)划分出的第l个频率段;
Figure BDA0002674863220000038
表示经过傅里叶变换后的第n个传感器接收的局部放电信号Yn(ω)划分出的第l个频率段的共轭转置;cum4表示虚拟扩展后的四阶累积量;
Figure BDA0002674863220000039
表示自然底数e为底数的指数函数;ω表示傅里叶变换的角频率,τm表示第m个传感器信号与参考传感器的时延差值,τn表示第n个传感器信号与参考传感器的时延差值。
相应地,本发明的另一目的在于提供一种提升变电站局部放电定向精度的系统,该提升变电站局部放电定向精度的系统无需硬件升级即可提高传感器阵列的性能,其通过将四阶累积量替换原始算法中二阶累积量,并利用最大值搜索求取局部放电方向角,从而有效提高在变电站现场中局部放电的定向精度。
根据上述的发明目的,本发明提出了一种提升变电站局部放电定向精度的系统,其包括:
传感器阵列,其接收局部放电信号;
处理单元,其执行步骤:
(1)基于接收的局部放电信号构建局部放电信号的响应功率;
(2)基于局部放电信号的响应功率构建确定局部放电信号的方向角θ的第一模型;
(3)对局部放电信号阵列进行虚拟扩展,并基于第一模型得到虚拟扩展后的局部放电信号的定向模型,并基于定向模型对局部放电进行精确定向。
进一步地,在本发明所述的提升变电站局部放电定向精度的系统中,在所述步骤(1)中:
将传感器阵列中的第m个传感器接收的局部放电信号记为ym(t),对其进行傅里叶变换后得为Ym(ω);
对局部放电信号Ym(ω)进行宽带波束形成处理,并将其频段划分为L个频率段,将每个频率段中各个传感器接收的局部放电信号加在一起以产生单个聚焦信号Z(ωl,θ):
Figure BDA0002674863220000041
其中,l=1,2,…,L;Yml)表示局部放电信号Ym(ω)划分出的第l个频率段;
Figure BDA0002674863220000042
表示与每个局部放电信号的第l个频率段对应的延时元件;τm表示第m个传感器信号与参考传感器的时延差值;Hml)表示第m个传感器对应的滤波器;M表示传感器阵列中的传感器的数量,m表示第m个传感器;
基于下式建立局部放电信号响应功率:
Figure BDA0002674863220000043
其中Z(ωl,θ)表示局部放电的单个聚焦信号;Z*l,θ)表示Z(ωl,θ)的共轭转置;
Figure BDA0002674863220000051
Figure BDA0002674863220000052
表示表示第n个传感器对应的滤波器的共轭转置;Wmnl)=E(Yml)·Yn *l)),其中E(·)表示数学期望,Yn *l)表示第n个传感器接收的局部放电信号傅里叶变换的共轭转置;第m个传感器与第n个传感器的组合滤波器Ψmn选为相变函数
Figure BDA0002674863220000053
Figure BDA0002674863220000054
表示自然底数e为底数的指数函数,j是复数的标识符号且有j2=-1;τn表示第n个传感器信号与参考传感器的时延差值。
进一步地,在本发明所述的提升变电站局部放电定向精度的系统中,在所述步骤(2)中,所述第一模型为:
Figure BDA0002674863220000055
其中,
Figure BDA0002674863220000056
Yml)表示经过傅里叶变换后的第m个传感器接收的局部放电信号Ym(ω)划分出的第l个频率段;
Figure BDA0002674863220000057
表示经过傅里叶变换后的第n个传感器接收的局部放电信号Yn(ω)划分出的第l个频率段的共轭转置;cum4表示虚拟扩展后的四阶累积量;
Figure BDA0002674863220000058
表示自然底数e为底数的指数函数;ω表示傅里叶变换的角频率,τm表示第m个传感器信号与参考传感器的时延差值,τn表示第n个传感器信号与参考传感器的时延差值。
本发明所述的提升变电站局部放电定向精度的方法及系统相较于现有技术具有如下所述的优点以及有益效果:
综上所述可以看出,本发明所述的提升变电站局部放电定向精度的方法无需硬件升级即可提高传感器阵列的性能,其通过将四阶累积量替换原始算法中二阶累积量,并利用最大值搜索求取局部放电方向角,可以有效提高在变电站现场中局部放电的定向精度。
相应地,本发明所述的提升变电站局部放电定向精度的系统也同样具有上述的优点以及有益效果。
附图说明
图1为本发明所述的提升变电站局部放电定向精度的方法在一种实施方式下的步骤流程示意图。
图2示意性地显示了本发明所述的提升变电站局部放电定向精度的方法采用传感器阵列接收局部放电信号经宽带信号波束形成处理的示意图。
图3示意性地显示了在进行局部放电定向实验中,原始SRP算法和阵列虚拟拓展后的SRP方法在不同入射角情况下定向误差对比。
具体实施方式
下面将结合说明书附图和具体的实施例对本发明所述的提升变电站局部放电定向精度的方法及系统做进一步的解释和说明,然而该解释和说明并不对本发明的技术方案构成不当限定。
图1为本发明所述的提升变电站局部放电定向精度的方法在一种实施方式下的步骤流程示意图。
如图1所示,在本实施方式中,本发明所述的提升变电站局部放电定向精度的方法,可以包括如下步骤:
(1)采用传感器阵列接收局部放电信号,并基于接收的局部放电信号构建局部放电信号的响应功率;
(2)基于局部放电信号的响应功率构建确定局部放电信号的方向角θ的第一模型;
(3)对局部放电信号阵列进行虚拟扩展,并基于第一模型得到虚拟扩展后的局部放电信号的定向模型,并基于定向模型对局部放电进行精确定向。
需要说明的是,本发明所述的提升变电站局部放电定向精度的方法,设计的主要思想是将原来的信号低阶累积量替换为由泰勒展开建立的高阶累积量。与低阶累积量相比,采用本发明所述方法建立的高阶累积量包含在低阶累积量中被省略的高阶项,因此它们具有有关信号源的更详细的方向信息,从而可以有效提高声波信号方向角(DOA)估计性能。
其中,在本发明所述方法的步骤(1)中,需要采用传感器阵列接收局部放电信号,并构建局部放电信号的响应功率,其可以包括下述内容:
首先利用传感器阵列接收局部放电信号,将传感器阵列中的第m个传感器接收的局部放电信号记为ym(t),对其进行傅里叶变换后得为Ym(ω);
而后可以对局部放电信号Ym(ω)进行宽带波束形成处理,并将其频段划分为L个频率段,将每个频率段中各个传感器接收的局部放电信号加在一起以产生单个聚焦信号Z(ωl),并最终构建局部放电信号响应功率。
需要说明的是,假设由参考传感器采集获得的局部放电信号为y(t),其可以构建如公式(1)所示的信号模型:
y(t)=s(t)·h(ds,t)+v(t) (1)
其中,t表示的是采样点;s(t)是声源发出的信号记为信号源矢量;h(ds,t)表示的是阵列响应矢量,其与声源到阵列的位置和方向有关;v(t)代表的噪音。
在传感器阵列中,第m个传感器接收的局部放电信号可以记为ym(t),其与该参考传感器间的信号接收时间差为τm。因此,结合参考上述公式(1),其信号间关系有:
ym(t)=y(t–τm) (2)
相应地,对上述ym(t)利用傅里叶变换转换为频域可以得到Ym(ω):
Figure BDA0002674863220000071
其中,Y(ω)和Ym(ω)分别表示为y(t)和ym(t)的傅里叶变换,e-jωt表示自然底数e为底数的指数函数,ω表示傅里叶变换的角频率,t表示时间采样点;
Figure BDA0002674863220000072
表示自然底数e为底数的指数函数,τm表示第m个传感器信号与参考传感器的时延差值。
此外,在本发明中,由于局部放电信号为典型的宽带信号,不能直接应用窄带信号的各种算法,因此在本发明所述的步骤(1)中,需要对局部放电信号Ym(ω)进行宽带波束形成处理,如图2所示。
图2示意性地显示了本发明所述的提升变电站局部放电定向精度的方法采用传感器阵列接收局部放电信号经宽带信号波束形成处理的示意图。
需要说明的是,在进行宽带波束形成处理时,首先需要将频域宽带信号分成L个小段,这些小段称为频率段,在每个频率段中,信号可以视为窄带信号,这意味着可以使用窄带信号的方法。每个单元格的宽度是采样频率除以离散傅里叶变换(DFT)的长度L。
为此,传感器在接收到局部放电信号ym(t)之后,首先采用离散傅里叶变换(DFT)将时域信号传输到频域。在图2中,每个信号的第l个频率段首先与对应的延时元件e–jωτ相乘,以确保它们是时间对准的,然后,信号与滤波器H(ωl)相乘,这些滤波器通常在波束形成过程中进行调整,以增强所需的源信号,同时衰减其他信号。
将每个频率段中各个传感器接收的局部放电信号加在一起以产生单个聚焦信号Z(ωl,θ):
Figure BDA0002674863220000081
其中,l=1,2,…,L;Yml)表示局部放电信号Ym(ω)划分出的第l个频率段;
Figure BDA0002674863220000082
表示与每个局部放电信号的第l个频率段对应的延时元件;Hml)表示第m个传感器对应的滤波器;M表示传感器阵列中的传感器的数量,m表示第m个传感器。
根据Z(ωl)可以基于下式建立局部放电信号响应功率:
Figure BDA0002674863220000083
其中Z(ωl,θ)表示局部放电的单个聚焦信号;Z*l,θ)表示Z(ωl,θ)的共轭转置;
Figure BDA0002674863220000084
Figure BDA0002674863220000085
表示第n个传感器对应的滤波器的共轭转置;Wmnl)=E(Yml)·Yn *l)),其中E(·)表示数学期望,Yn *l)表示第n个传感器接收的局部放电信号傅里叶变换的共轭转置;第m个传感器与第n个传感器的组合滤波器Ψmn选为相变函数
Figure BDA0002674863220000086
Figure BDA0002674863220000087
表示自然底数e为底数的指数函数,j是复数的标识符号且有j2=-1;τn表示第n个传感器信号与参考传感器的时延差值。
因此,在本发明所述的步骤(2)中,基于局部放电信号的响应功率构建能够确定局部放电信号的方向角θ的第一模型:
θs=argmaxθP(θ) (6)
其中,P(θ)表示局部放电信号的响应功率,argmax()表示函数取最大值时对应的θ。
需要说明的是,在本发明所述的方法中,可以将原来的信号低阶累积量替换为由泰勒展开建立的高阶累积量。高阶累积量代表的是某一随机变量特征函数的泰勒展开的前k项的系数,因此对于高阶累积量来説,其阶数越高,内部所包含的信息则越多。
在上述原始的可控响应功率(Steered response power,SRP)定向方法公式(5)中,由于信号均值设为零,Wmnl)其实表示的就是信号的二阶累积量,若将其替换为四阶累积量,可以完成更高阶信号信息的利用,从而完成阵列的扩展。
将原始SRP算法公式(5)中信号(Yml),Ynl))间的二阶累积量Wmnl)=E(Yml)·Yn *l))写为:
Figure BDA0002674863220000091
其中,
Figure BDA0002674863220000092
表示的是噪声功率,其为一常数;ejωτ表示自然底数e为底数的指数函数,τ=τmn表示传感器m和传感器n接收信号的时间差值;Yml)表示的是经过傅里叶变换后的第m个传感器接收的局部放电信号Ym(ω)划分出的第l个频率段,Ym *l)表示的是Yml)的共轭转置。
将上述公式(7)可以替换为如下形式的四阶累积量:
Figure BDA0002674863220000093
其中,γ4,m表示的是Ym(ω)的自四阶累积量,同样为一常数;ejω2τ表示的是自然底数e为底数的指数函数,ω表示傅里叶变换的角频率,τ=τmn表示传感器m和传感器n接收信号的时间差值。由于γ4,m
Figure BDA0002674863220000094
均为常数,因此并不会影响本发明所述步骤(2)中的公式(6)所示的最大值搜索,因此替换后这两个常数不会影响最后的搜索结果。
比较上述公式(7)和公式(8)可知,扩展后的信号四阶累积量,两者的形式相同,其最大的区别在于二阶累积量中的时延值τ在四阶累积量中被提升了一倍至2τ。
将上述公式(8)所构建的高阶累积量代入的公式(6)中,则可得到扩展后的局部放电信号的定向模型,从而对局部放电进行精确定向。因此,在本发明所述的步骤(3)中,基于第一模型可以得到阵列虚拟扩展后的局部放电信号的定向模型:
Figure BDA0002674863220000101
其中,
Figure BDA0002674863220000102
Yml)表示表示经过傅里叶变换后的第m个传感器接收的局部放电信号Ym(ω)划分出的第l个频率段;
Figure BDA0002674863220000103
表示经过傅里叶变换后的第n个传感器接收的局部放电信号Yn(ω)划分出的第l个频率段的共轭转置;cum4表示虚拟扩展后的四阶累积量,具体函数形式如公式(8)所示;
Figure BDA0002674863220000104
表示自然底数e为底数的指数函数;ω表示傅里叶变换的角频率,τm表示第m个传感器信号与参考传感器的时延差值,τn同理。
为了更好地说明本发明所述的提升变电站局部放电定向精度的方法的应用情况,分别利用原始SRP算法以及阵列扩展后的SRP方法进行局部放电定向实验,以进行进一步说明。
在局部放电定向实验中,采用了一种提升变电站局部放电定向精度的系统,其可以用于实施本发明所述的提升变电站局部放电定向精度的方法。其中,该系统可以包括:传感器阵列和处理单元。
需要说明的是,系统中的传感器阵列能够用于接收局部放电信号;处理单元可以用于执行除接收局部放电信号以外的提升变电站局部放电定向精度方法中的其他所有操作。
表1中列出了在不同信噪比的情况下,基于原始SRP算法以及阵列拓展后SRP方法的定位试验结果。
表1.
Figure BDA0002674863220000105
Figure BDA0002674863220000111
由表1的定位试验结果可以看出,相较于原始的SRP算法,本发明所述的提升变电站局部放电定向精度的系统中在利用高阶累积量进行阵列扩展后。有利于提升局部放电的定向精度。定位试验结果表明,孔径虚拟扩展后,阵列方向定向精度提高了约40%,并且显示出一定的抗噪声能力,在高噪声水平的环境下,现场测试结果也验证了所提出方法的有效性。由此可见,本发明所述的提升变电站局部放电定向精度的方法及系统可能会促进基于传感器阵列的电力设备在线监测的应用。
图3示意性地显示了在进行局部放电定向实验中,原始SRP算法和阵列虚拟拓展后的SRP方法在不同入射角情况下定向误差对比。
在本实施方式中,分别利用原始SRP算法以及阵列拓展后的SRP方法进行局部放电试验,并不断变换局部放电对于传感器阵列的方向角,其结果如图3所示。
结合参阅图3可知,在各个方向上,局域拓展后的阵列在一定程度上能够提升局部放电源的定位精度。
综上所述可以看出,本发明所述的提升变电站局部放电定向精度的方法无需硬件升级即可提高传感器阵列的性能,其通过将四阶累积量替换原始算法中二阶累积量,并利用最大值搜索求取局部放电方向角,可以有效提高在变电站现场中局部放电的定向精度。
相应地,本发明所述的提升变电站局部放电定向精度的系统也同样具有上述的优点以及有益效果。
需要说明的是,本发明的保护范围中现有技术部分并不局限于本申请文件所给出的实施例,所有不与本发明的方案相矛盾的现有技术,包括但不局限于在先专利文献、在先公开出版物,在先公开使用等等,都可纳入本发明的保护范围。
此外,本案中各技术特征的组合方式并不限本案权利要求中所记载的组合方式或是具体实施例所记载的组合方式,本案记载的所有技术特征可以以任何方式进行自由组合或结合,除非相互之间产生矛盾。
还需要注意的是,以上所列举的实施例仅为本发明的具体实施例。显然本发明不局限于以上实施例,随之做出的类似变化或变形是本领域技术人员能从本发明公开的内容直接得出或者很容易便联想到的,均应属于本发明的保护范围。

Claims (4)

1.一种提升变电站局部放电定向精度的方法,其特征在于,包括步骤:
(1)采用传感器阵列接收局部放电信号,并基于接收的局部放电信号构建局部放电信号的响应功率;
(2)基于局部放电信号的响应功率构建确定局部放电信号的方向角θ的第一模型;其中所述第一模型为:
θs=argmaxθP(θ)
其中,P(θ)表示局部放电信号的响应功率,argmaxθ()表示函数取最大值时对应的θ;
(3)对局部放电信号阵列进行虚拟扩展,并基于第一模型得到虚拟扩展后的局部放电信号的定向模型,并基于定向模型对局部放电进行精确定向,其中所述定向模型被构建为:
Figure FDA0003301729520000011
其中,
Figure FDA0003301729520000012
Yml)表示经过傅里叶变换后的第m个传感器接收的局部放电信号Ym(ω)划分出的第l个频率段;
Figure FDA0003301729520000013
表示经过傅里叶变换后的第n个传感器接收的局部放电信号Yn(ω)划分出的第l个频率段的共轭转置;cum4表示虚拟扩展后的四阶累积量;
Figure FDA0003301729520000014
表示自然底数e为底数的指数函数;ω表示傅里叶变换的角频率,τm表示第m个传感器信号与参考传感器的时延差值,τn表示第n个传感器信号与参考传感器的时延差值;
其中所述定向模型基于如下形式的四阶累积量获得:
Figure FDA0003301729520000015
其中,γ4,m表示的是Ym(ω)的自四阶累积量,同样为一常数;ejω2τ表示的是自然底数e为底数的指数函数,ω表示傅里叶变换的角频率,τ=τmn表示第m个传感器和第n个传感器接收信号的时间差值。
2.如权利要求1所述的提升变电站局部放电定向精度的方法,其特征在于,在所述步骤(1)中:
将传感器阵列中的第m个传感器接收的局部放电信号记为ym(t),对其进行傅里叶变换后得为Ym(ω);
对局部放电信号Ym(ω)进行宽带波束形成处理,并将其频段划分为L个频率段,将每个频率段中各个传感器接收的局部放电信号加在一起以产生单个聚焦信号Z(ωl,θ):
Figure FDA0003301729520000021
其中,l=1,2,…,L;Yml)表示局部放电信号Ym(ω)划分出的第l个频率段;
Figure FDA0003301729520000022
表示与每个局部放电信号的第l个频率段对应的延时元件,τm表示第m个传感器信号与参考传感器的时延差值;Hml)表示第m个传感器对应的滤波器;M表示传感器阵列中的传感器的数量,m表示第m个传感器;
基于下式建立局部放电信号响应功率:
Figure FDA0003301729520000023
其中Z(ωl,θ)表示局部放电的单个聚焦信号;Z*l,θ)表示Z(ωl,θ)的共轭转置;
Figure FDA0003301729520000024
Figure FDA0003301729520000025
表示第n个传感器对应的滤波器的共轭转置;Wmnl)=E(Yml)·Yn *l)),其中E(·)表示数学期望,Yn *l)表示第n个传感器接收的局部放电信号傅里叶变换的共轭转置;第m个传感器与第n个传感器的组合滤波器Ψmn选为相变函数
Figure FDA0003301729520000026
Figure FDA0003301729520000027
表示自然底数e为底数的指数函数,j是复数的标识符号且有j2=-1;τn表示第n个传感器信号与参考传感器的时延差值。
3.一种提升变电站局部放电定向精度的系统,其特征在于,包括:
传感器阵列,其接收局部放电信号;
处理单元,其执行步骤:
(1)基于接收的局部放电信号构建局部放电信号的响应功率;
(2)基于局部放电信号的响应功率构建确定局部放电信号的方向角θ的第一模型;其中所述第一模型为:
θs=argmaxθP(θ)
其中,P(θ)表示局部放电信号的响应功率,argmaxθ()表示函数取最大值时对应的θ;
(3)对局部放电信号阵列进行虚拟扩展,并基于第一模型得到虚拟扩展后的局部放电信号的定向模型,并基于定向模型对局部放电进行精确定向;其中所述定向模型被构建为:
Figure FDA0003301729520000031
其中,
Figure FDA0003301729520000032
Yml)表示经过傅里叶变换后的第m个传感器接收的局部放电信号Ym(ω)划分出的第l个频率段;
Figure FDA0003301729520000033
表示经过傅里叶变换后的第n个传感器接收的局部放电信号Yn(ω)划分出的第l个频率段的共轭转置;cum4表示虚拟扩展后的四阶累积量;
Figure FDA0003301729520000034
表示自然底数e为底数的指数函数;ω表示傅里叶变换的角频率,τm表示第m个传感器信号与参考传感器的时延差值,τn表示第n个传感器信号与参考传感器的时延差值;
其中所述定向模型基于如下形式的四阶累积量获得:
Figure FDA0003301729520000035
其中,γ4,m表示的是Ym(ω)的自四阶累积量,同样为一常数;ejω2τ表示的是自然底数e为底数的指数函数,ω表示傅里叶变换的角频率,τ=τmn表示第m个传感器和第n个传感器接收信号的时间差值。
4.如权利要求3所述的提升变电站局部放电定向精度的系统,其特征在于,在所述步骤(1)中:
将传感器阵列中的第m个传感器接收的局部放电信号记为ym(t),对其进行傅里叶变换后得为Ym(ω);
对局部放电信号Ym(ω)进行宽带波束形成处理,并将其频段划分为L个频率段,将每个频率段中各个传感器接收的局部放电信号加在一起以产生单个聚焦信号Z(ωl,θ):
Figure FDA0003301729520000036
其中,l=1,2,…,L;Yml)表示局部放电信号Ym(ω)划分出的第l个频率段;
Figure FDA0003301729520000041
表示与每个局部放电信号的第l个频率段对应的延时元件,τm表示第m个传感器信号与参考传感器的时延差值;Hml)表示第m个传感器对应的滤波器;M表示传感器阵列中的传感器的数量,m表示第m个传感器;
基于下式建立局部放电信号响应功率:
Figure FDA0003301729520000042
其中Z(ωl,θ)表示局部放电的单个聚焦信号;Z*l,θ)表示Z(ωl,θ)的共轭转置;
Figure FDA0003301729520000043
Figure FDA0003301729520000044
表示第n个传感器对应的滤波器的共轭转置;Wmnl)=E(Yml)·Yn *l)),其中E(·)表示数学期望,Yn *l)表示第n个传感器接收的局部放电信号傅里叶变换的共轭转置;第m个传感器与第n个传感器的组合滤波器Ψmn选为相变函数
Figure FDA0003301729520000045
Figure FDA0003301729520000046
表示自然底数e为底数的指数函数,j是复数的标识符号且有j2=-1;τn表示第n个传感器信号与参考传感器的时延差值。
CN202010944714.3A 2020-09-10 2020-09-10 一种提升变电站局部放电定向精度的方法及系统 Active CN112067953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010944714.3A CN112067953B (zh) 2020-09-10 2020-09-10 一种提升变电站局部放电定向精度的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010944714.3A CN112067953B (zh) 2020-09-10 2020-09-10 一种提升变电站局部放电定向精度的方法及系统

Publications (2)

Publication Number Publication Date
CN112067953A CN112067953A (zh) 2020-12-11
CN112067953B true CN112067953B (zh) 2022-02-25

Family

ID=73663519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010944714.3A Active CN112067953B (zh) 2020-09-10 2020-09-10 一种提升变电站局部放电定向精度的方法及系统

Country Status (1)

Country Link
CN (1) CN112067953B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113050038A (zh) * 2021-03-23 2021-06-29 上海交通大学 一种基于虚拟阵列扩展的变电站声源定位方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702000A (zh) * 2009-11-26 2010-05-05 华北电力大学(保定) 变压器局部放电的相控超声定位方法及系统
CN102253127A (zh) * 2011-07-01 2011-11-23 西安交通大学 用于变压器局部放电定位的l形超声阵列传感器及其方法
CN103995221A (zh) * 2014-05-29 2014-08-20 西安交通大学 用于变压器局部放电定位的超声波接收阵列的定向方法
CN108037410A (zh) * 2017-11-24 2018-05-15 华北电力大学 一种基于可控响应功率的局部放电超声波定位方法与装置
CN108181557A (zh) * 2017-12-29 2018-06-19 上海交通大学 一种确定特高频局部放电信号方位的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702000A (zh) * 2009-11-26 2010-05-05 华北电力大学(保定) 变压器局部放电的相控超声定位方法及系统
CN102253127A (zh) * 2011-07-01 2011-11-23 西安交通大学 用于变压器局部放电定位的l形超声阵列传感器及其方法
CN103995221A (zh) * 2014-05-29 2014-08-20 西安交通大学 用于变压器局部放电定位的超声波接收阵列的定向方法
CN108037410A (zh) * 2017-11-24 2018-05-15 华北电力大学 一种基于可控响应功率的局部放电超声波定位方法与装置
CN108181557A (zh) * 2017-12-29 2018-06-19 上海交通大学 一种确定特高频局部放电信号方位的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于虚拟阵列扩展的特高频多源局部放电DOA估计方法;周南等;《高电压技术》;20190531;第45卷(第5期);第1551-1558页 *
用于变压器中局部放电定向的超声阵列RSS改进算法;王国利等;《高电压技术》;20190831;第45卷(第8期);第2509-2514页 *

Also Published As

Publication number Publication date
CN112067953A (zh) 2020-12-11

Similar Documents

Publication Publication Date Title
CN101852826B (zh) 一种电力系统的谐波分析方法及其装置
JP2846125B2 (ja) 超音波測定信号の適応的最適化方法
CN111222088B (zh) 一种改进的平顶自卷积窗加权电力谐波幅值估计方法
CN111798869B (zh) 一种基于双麦克风阵列的声源定位方法
CN112067953B (zh) 一种提升变电站局部放电定向精度的方法及系统
CN109889186A (zh) 一种基于多级滤波器组的宽带波束形成方法
CN109343058B (zh) 基于混合算法的正交非线性调频信号生成方法及装置
CN110705041A (zh) 一种基于easi的线性结构工作模态参数识别方法
CN109507643A (zh) 旁瓣相消中展宽零陷的方法
CN109074814A (zh) 一种噪声检测方法及终端设备
CN110047507A (zh) 一种声源识别方法及装置
JP6441769B2 (ja) クラスタリング装置、クラスタリング方法及びクラスタリングプログラム
CN113299307A (zh) 麦克风阵列信号处理方法、系统、计算机设备及存储介质
US11482239B2 (en) Joint source localization and separation method for acoustic sources
CN102598126A (zh) 信息处理装置及其附属设备、信息处理系统及其控制方法和控制程序
JP6266178B2 (ja) 方探装置
CN109061597B (zh) 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法
Lu et al. Online blind source separation method with adaptive step size in both time-invariant and time-varying cases
CN108507743B (zh) 一种基于lmd的冲击响应谱时域信号合成的方法及装置
Washizawa et al. A flexible method for envelope estimation in empirical mode decomposition
CN110063751B (zh) 电源系统及其控制方法
CN114184838A (zh) 基于sn互卷积窗的电力系统谐波检测方法、系统及介质
Tao et al. A novel approach to fast DOA estimation of multiple spatial narrowband signals
Ando Time–frequency representation with variant array of frequency-domain Prony estimators
Fan et al. Matching pursuit based on nonparametric waveform estimation

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