CN113050036B - 基于多谐振点传声器阵列的gis波束形成定位方法 - Google Patents

基于多谐振点传声器阵列的gis波束形成定位方法 Download PDF

Info

Publication number
CN113050036B
CN113050036B CN202110278865.4A CN202110278865A CN113050036B CN 113050036 B CN113050036 B CN 113050036B CN 202110278865 A CN202110278865 A CN 202110278865A CN 113050036 B CN113050036 B CN 113050036B
Authority
CN
China
Prior art keywords
microphone array
resonance point
gis
positioning
beam forming
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
CN202110278865.4A
Other languages
English (en)
Other versions
CN113050036A (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.)
Dongying Power Industry Bureau Of State Grid Shandong Electric Power Co
State Grid Corp of China SGCC
Original Assignee
Dongying Power Industry Bureau Of State Grid Shandong Electric Power Co
State Grid Corp of China SGCC
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 Dongying Power Industry Bureau Of State Grid Shandong Electric Power Co, State Grid Corp of China SGCC filed Critical Dongying Power Industry Bureau Of State Grid Shandong Electric Power Co
Priority to CN202110278865.4A priority Critical patent/CN113050036B/zh
Publication of CN113050036A publication Critical patent/CN113050036A/zh
Application granted granted Critical
Publication of CN113050036B publication Critical patent/CN113050036B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明公开了一种基于多谐振点传声器阵列的GIS波束形成定位方法,涉及波束形成定位领域,包括以下步骤:S1、采集传声器阵列接受到的声音信号;S2、根据声音信号,分别通过求根MUSIC算法获得多谐振点传声器阵列中不同频率范围的声场分布;S3、结合基于多谐振点传声器阵列改进的自适应加权融合算法进行多谐振点传声器阵列的不同频率范围声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源。本发明通过多谐振点传声器阵列接收到声音信号,综合运用多种多个频段的波束形成结果的融合,对故障点进行定位,不仅能确定噪声源的具体位置,提高定位精度,还能抑制干扰信号的波束输出,提高分辨率和抗干扰能力。

Description

基于多谐振点传声器阵列的GIS波束形成定位方法
技术领域
本发明涉及波束形成定位领域,尤其涉及一种基于多谐振点传声器阵列的GIS波束形成定位方法。
背景技术
近年来随着电力事业的发展,GIS振动噪声所带来的危害越来越多地受到人们的关注。GIS发生故障会产生振动,振动会产生噪声,噪声中蕴含GIS工作状态的重要信息。通过采集GIS的噪声信号,并对噪声信号进行时域、频域处理,可以得出GIS的工作状态,并且能够很好地反映GIS设备的故障。现阶段我国对GIS故障点的定位还处于起步阶段,且由于GIS内部结构复杂,有多种电力设备,理论计算很难准确分析噪声的产生及传播规律。
GIS的常规检修过程是采用人耳来分辨噪声,可靠性较低。现有的噪声研究还停留在测量频谱、声压级阶段,对采集到的噪声信号进行时域和频域上的处理,这些研究远不能满足深入了解GIS振动噪声传播机理。而采用本发明的定位技术来定位GIS的噪声源,能过分辨出GIS的噪声空间分布,类似于红外温度热成像技术,能够直观的发现GIS的故障点,简单实用,用于带电检测能够很好的发现GIS的潜在缺陷及故障,避免设备事故的发生,正逐渐成为一种重要的GIS状态分析诊断方法。
传统的GIS定位方法基于波束形成定位、到达时间差定位和高分辨率谱估计定位。在真实GIS上开展噪声源的定位,由于GIS内部复杂结构、振动噪声信号传播机理的复杂性以及环境噪声等因素均不利于传播时间值的准确测量,导致现有定位方法在时延估计方面存在一定的误差,容易造成定位精度变低,甚至出现定位失败的现象。因此,为了更好的实现GIS噪声源的定位,提出一种时延估计准确,定位精度高的GIS故障点定位方法,具有非常重要的意义。
发明内容
为解决上述技术问题,本发明提供一种分辨率高,抗干扰能力强,定位精度高,基于多谐振点传声器阵列的GIS波束形成定位方法。
本发明采用如下技术方案:
基于多谐振点传声器阵列的GIS的波束形成定位方法,所述方法通过多谐振点传声器阵列接收到声音信号,利用求根MUSIC算法获得不同频率范围的声场分布,结合基于多谐振点传声器阵列改进的自适应加权融合算法进行多谐振点传声器阵列的不同频率范围声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源。
基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于,包括以下步骤:
S1、采集传声器阵列接受到的声音信号;
S2、根据声音信号,分别通过求根MUSIC算法获得多谐振点传声器阵列中不同频率范围的声场分布;
S3、结合基于多谐振点传声器阵列改进的自适应加权融合算法进行多谐振点传声器阵列的不同频率范围声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源。
进一步的,所述步骤S2中通过求根MUSIC算法获得不同频率范围的声场分布,具体过程为:
(1)根据传声器阵列接受到的信号得到协方差矩阵的估计值R;
(2)对得到的协方差矩阵进行特征值分解,得到噪声子空间;
(3)构造求根MUSIC多项式,即
Figure BDA0002977581700000021
(5)对多项式求根,得到波达方向
Figure BDA0002977581700000022
进一步的,所述步骤S3中结合基于多谐振点传声器阵列改进的自适应加权融合算法进行多谐振点传声器阵列的不同频率范围声场叠加,通过主要成分分析法得到最终的噪声源定位,具体过程为:
(1)求不同频率范围的声场分布矩阵X的方差σi 2,根据多谐振点传声器阵列中传声器分组个数n,可以求得n个声场分布矩阵X的方差σi 2
(2)根据下面公式,求取各个声场分布矩阵的加权因子;
Figure BDA0002977581700000031
(3)不同频率范围声场根据其不同加权因子进行叠加,得到最终的融合声场;
(4)求最终融合声场的相关系数矩阵R;
(5)求解其特征值和特征向量;
(6)提取主成分;
(5)确定主成分,得到噪声源的定位。
进一步的,所述步骤5中具体的过程为,求解系数矩阵R的特征方程,解方程求出系数矩阵R的特征值和特征向量。
进一步的,所述步骤6中具体的过程为,将特征值由大到小排序,提取特征值大于1的m个主成分。
进一步的,在步骤S1之前,还包括进行GIS故障点精确定位硬件架设。
进一步的,所述进行GIS故障点精确定位硬件架设,包括传声器阵列、阵列支架、信号调理单元、多通道同步数据采集卡和终端计算机。
进一步的,所述传声器阵列由阵列支架固定,将各传声器采集到的信号与多通道同步数据采集卡相连,终端计算机与数据采集卡连接,用于接收和储存处理后的声音信号。
本发明具有以下有益效果:本发明的基于多谐振点传声器阵列的GIS波束形成定位方法,使用多谐振点传声器阵列,通过求根MUSIC算法获得多谐振点传声器阵列中不同频率范围的声场分布,结合基于多谐振点传声器阵列改进的自适应加权融合算法进行声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源,不仅能确定噪声源的具体位置,提高定位精度,还能抑制干扰信号的波束输出,提高分辨率和抗干扰能力。同时本发明的波束形成定位方法较传统宽带波束形成算法计算量低,计算速度快,充分发挥多谐振点传声器阵列的硬件优势舍弃传统波束形成算法中冗余计算的部分,大大提高计算速率,更加有效的利用噪声信号评估GIS的运行状态,确保GIS的安全稳定运行。
附图说明
下面结合附图对本发明作优选的说明:
图1为本发明的流程图;
图2为GIS噪声源定位的原理示意图;
图3为求根MUSIC算法的流程图;
图4为主成分分析法的流程图;
图5为自适应融合算法的流程图。
具体实施方式
本发明的目的是提供一种基于多谐振点传声器阵列的GIS波束形成定位方法,该方法通过通过求根MUSIC算法获得多谐振点传声器阵列中不同频率范围的声场分布,结合主成分分析法和基于多谐振点传声器阵列改进的自适应加权融合算法进行声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源,不仅可以提高噪声源的定位精度,还可提高GIS的噪声检测水准,同时本发明的波束形成定位方法较传统宽带波束形成算法计算量低,计算速度快,充分发挥多谐振点传声器阵列的硬件优势舍弃传统波束形成算法中冗余计算的部分,大大提高计算速率,更加有效的利用噪声信号评估GIS的运行状态。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,一种基于多谐振点传声器阵列的GIS波束形成定位方法,包括以下步骤:
步骤1,如图2所示,进行GIS故障点精确定位硬件架设,包括传声器阵列、阵列支架、信号调理单元、多通道同步数据采集卡和终端计算机,传声器阵列由阵列支架固定,将各传声器采集到的信号与多通道同步数据采集卡相连,终端计算机与数据采集卡连接,用于接收和储存处理后的声音信号,以声音信号作为定位触发信号,采用多种定位算法对声音信号作定位分析。
步骤2,如图3所示,对同一个待测噪声源,采用求根MUSIC算法,获得声场分布。
采用求根MUSIC算法求故障点方位角的具体过程如下:
(1)根据传声器阵列接受到的信号得到协方差矩阵的估计值R;
(2)对得到的协方差矩阵进行特征值分解,得到噪声子空间;
(3)构造求根MUSIC多项式,即
Figure BDA0002977581700000051
(4)对多项式求根,找到最接近单位圆内/外K个根,得到波达方向
Figure BDA0002977581700000052
步骤3,如图4和图5所示,结合基于多谐振点传声器阵列改进的自适应融合算法进行声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源。具体过程如下:
(1)求不同频率范围的声场分布矩阵X的方差σi 2,根据多谐振点传声器阵列中传声器分组个数n,可以求得n个声场分布矩阵X的方差σi 2
(2)根据下面公式,求取各个不同频率声场分布矩阵的加权因子。
Figure BDA0002977581700000053
(3)不同频率范围声场根据其不同加权因子进行叠加,得到最终的融合声场。
(4)求最终融合声场的相关系数矩阵R
(5)求解其特征值和特征向量;
(6)提取主成分;
(5)确定主成分,得到噪声源的定位。
本发明提供的一种基于多谐振点传声器阵列的GIS波束形成定位方法,使用多谐振点传声器阵列,通过求根MUSIC算法获得多谐振点传声器阵列中不同频率范围的声场分布,结合基于多谐振点传声器阵列改进的自适应加权融合算法进行声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源,不仅能确定噪声源的具体位置,提高定位精度,还能抑制干扰信号的波束输出,提高分辨率和抗干扰能力。同时本发明的波束形成定位方法较传统宽带波束形成算法计算量低,计算速度快,充分发挥多谐振点传声器阵列的硬件优势舍弃传统波束形成算法中冗余计算的部分,大大提高计算速率,更加有效的利用噪声信号评估GIS的运行状态,确保GIS的安全稳定运行。
以上仅为本发明的具体实施例,但本发明的技术特征并不局限于此。任何以本发明为基础,为解决基本相同的技术问题,实现基本相同的技术效果,所作出的简单变化、等同替换或者修饰等,皆涵盖于本发明的保护范围之中。

Claims (7)

1.基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于,包括以下步骤:
S1、采集传声器阵列接受到的声音信号;
S2、根据声音信号,分别通过求根MUSIC算法获得多谐振点传声器阵列中不同频率范围的声场分布;
S3、结合基于多谐振点传声器阵列改进的自适应加权融合算法进行多谐振点传声器阵列的不同频率范围声场叠加,通过主要成分分析法得到最终的噪声源定位,找到最终的噪声源;
所述步骤S3中结合基于多谐振点传声器阵列改进的自适应加权融合算法进行多谐振点传声器阵列的不同频率范围声场叠加,通过主要成分分析法得到最终的噪声源定位,具体过程为:
(1)求不同频率范围的声场分布矩阵X的方差σi 2,根据多谐振点传声器阵列中传声器分组个数n,可以求得n个声场分布矩阵X的方差σi 2
(2)根据下面公式,求取各个不同频率声场分布矩阵的加权因子;
Figure FDA0003574434020000011
(3)不同频率范围声场根据其不同加权因子进行叠加,得到最终的融合声场;
(4)求最终融合声场的相关系数矩阵R;
(5)求解其特征值和特征向量;
(6)提取主成分;
(7)确定主成分,得到噪声源的定位。
2.根据权利要求1所述的基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于:所述步骤S2中通过求根MUSIC算法获得不同频率范围的声场分布,具体过程为:
(1)根据传声器阵列接受到的信号得到协方差矩阵的估计值R;
(2)对得到的协方差矩阵进行特征值分解,得到噪声子空间;
(3)构造求根MUSIC多项式,即
Figure FDA0003574434020000021
(4)对多项式求根,得到波达方向
Figure FDA0003574434020000022
3.根据权利要求1所述的基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于:所述步骤5中具体的过程为,求解系数矩阵R的特征方程,解方程求出系数矩阵R的特征值和特征向量。
4.根据权利要求1所述的基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于:所述步骤6中具体的过程为,将特征值由大到小排序,提取特征值大于1的m个主成分。
5.根据权利要求1所述的基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于:在步骤S1之前,还包括进行GIS故障点精确定位硬件架设。
6.根据权利要求5所述的基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于:所述进行GIS故障点精确定位硬件架设,包括传声器阵列、阵列支架、信号调理单元、多通道同步数据采集卡和终端计算机。
7.根据权利要求6所述的基于多谐振点传声器阵列的GIS波束形成定位方法,其特征在于:所述传声器阵列由阵列支架固定,将各传声器采集到的信号与多通道同步数据采集卡相连,终端计算机与数据采集卡连接,用于接收和储存处理后的声音信号。
CN202110278865.4A 2021-03-16 2021-03-16 基于多谐振点传声器阵列的gis波束形成定位方法 Active CN113050036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110278865.4A CN113050036B (zh) 2021-03-16 2021-03-16 基于多谐振点传声器阵列的gis波束形成定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110278865.4A CN113050036B (zh) 2021-03-16 2021-03-16 基于多谐振点传声器阵列的gis波束形成定位方法

Publications (2)

Publication Number Publication Date
CN113050036A CN113050036A (zh) 2021-06-29
CN113050036B true CN113050036B (zh) 2022-05-17

Family

ID=76512720

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110278865.4A Active CN113050036B (zh) 2021-03-16 2021-03-16 基于多谐振点传声器阵列的gis波束形成定位方法

Country Status (1)

Country Link
CN (1) CN113050036B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110687506A (zh) * 2019-10-11 2020-01-14 国网陕西省电力公司电力科学研究院 一种基于矢量传声器阵列的低频噪声源定位装置及方法
CN111175698A (zh) * 2020-01-18 2020-05-19 国网山东省电力公司菏泽供电公司 一种基于声振结合的变压器噪声源定位方法、系统和装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201510220D0 (en) * 2015-06-11 2015-07-29 Renishaw Plc Additive manufacturing apparatus and method
CN106054130B (zh) * 2016-06-06 2019-01-22 南京工程学院 一种基于music算法的室内定位方法及装置
CN109901110A (zh) * 2019-02-25 2019-06-18 西安电子科技大学 基于主成分分析的支持向量机近场声源定位方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110687506A (zh) * 2019-10-11 2020-01-14 国网陕西省电力公司电力科学研究院 一种基于矢量传声器阵列的低频噪声源定位装置及方法
CN111175698A (zh) * 2020-01-18 2020-05-19 国网山东省电力公司菏泽供电公司 一种基于声振结合的变压器噪声源定位方法、系统和装置

Also Published As

Publication number Publication date
CN113050036A (zh) 2021-06-29

Similar Documents

Publication Publication Date Title
CN109375171B (zh) 一种基于正交匹配追踪算法的声源定位方法
CN111024209B (zh) 一种适用于矢量水听器的线谱检测方法
CN111175698B (zh) 一种基于声振结合的变压器噪声源定位方法、系统和装置
CN110837001B (zh) 一种电力系统中谐波和间谐波的分析方法与装置
CN109188362B (zh) 一种麦克风阵列声源定位信号处理方法
CN109683134B (zh) 一种面向旋转声源的高分辨率定位方法
CN108957403B (zh) 一种基于广义互相关的高斯拟合包络时延估计方法及系统
CN110687528B (zh) 自适应波束形成器生成方法及系统
Zheng et al. A RSSI-AOA-based UHF partial discharge localization method using MUSIC algorithm
CN110082741A (zh) 一种基于伪数据重构的超分辨波达角估计算法
Lan et al. Partial discharge location of power cables based on an improved phase difference method
CN110940933B (zh) 一种用于测量陡脉冲上升沿起始时刻的综合计算方法
CN109991519B (zh) 基于神经网络和无线传感阵列的局部放电测向方法及系统
CN113050036B (zh) 基于多谐振点传声器阵列的gis波束形成定位方法
CN109752633B (zh) 一种对变电站局部放电信号进行定位的方法及系统
CN107861115A (zh) 一种基于瞬时自相关矩阵稀疏分解的othr机动目标参数估计方法
CN115497501A (zh) 基于sw-music的变压器故障声纹定位方法及系统
Li et al. A new iterative near-field coherent subspace method for rub-impact fault localization using AE technique
CN114201991A (zh) 一种基于超声传感器阵列的局放信号检测方法和系统
CN115291054A (zh) 一种基于小波图像卷积的非梯度定位方法
CN109581422A (zh) 一种基于子空间误差修正的线性扫频干扰抑制方法及系统
Zhang et al. A Low Frequency Source of Transformer Localization Method Using the Combination of CLEAN algorithm based on Spatial Source Coherence and Multi-signal Classification Algorithm
CN114858271B (zh) 一种声音探测用的阵列放大方法
Liu et al. Localization of Multi-source UHF Signals Based on Random Matrix Theory
CN113721184B (zh) 基于改进music算法的近场信号源定位方法

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