CN111413649B - 基于近场宽带波束形成的大型电抗器故障检测方法及系统 - Google Patents

基于近场宽带波束形成的大型电抗器故障检测方法及系统 Download PDF

Info

Publication number
CN111413649B
CN111413649B CN202010327774.0A CN202010327774A CN111413649B CN 111413649 B CN111413649 B CN 111413649B CN 202010327774 A CN202010327774 A CN 202010327774A CN 111413649 B CN111413649 B CN 111413649B
Authority
CN
China
Prior art keywords
reactor
array
sound intensity
sub
sound
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
CN202010327774.0A
Other languages
English (en)
Other versions
CN111413649A (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.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
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 University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN202010327774.0A priority Critical patent/CN111413649B/zh
Publication of CN111413649A publication Critical patent/CN111413649A/zh
Application granted granted Critical
Publication of CN111413649B publication Critical patent/CN111413649B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups

Abstract

本发明公开了一种基于近场宽带波束形成的大型电抗器故障检测方法及系统，通过比较电抗器的声音强度或者能量图像与其正常工作时的声音强度或者能量图像，判断电抗器设备是否存在故障需要工作人员维护，进而实现了电抗器设备工作状态的在线检测。通过近场宽带波束形成的声音强度或者能量图像来判断电抗器工作状态，是一种在线检测方案，不仅降低了维护的成本，而且还提高了人们用电的质量，给大型电抗器设备的工作状态检测提供了一种新的方法。

Description

基于近场宽带波束形成的大型电抗器故障检测方法及系统
技术领域
本发明涉及大型电抗器故障检测领域，尤其涉及一种基于近场宽带波束形成的大型电抗器故障检测方法及系统。
背景技术
现如今，电力需求不断提升，用电质量问题也受到广泛的关注，电能的稳定供应与国民经济也有着密不可分的关系。大量实践证明，电抗器在出现严重事故之前，往往都有一定的潜伏性故障。电抗器内部的局部放电，局部过热，绕组变形、机械部件松动以及设备绝缘老化等潜伏性故障是随着时间积累引起的。电抗器运行时，由于机身与固件、零件之间或者零件自身之间的相互运动，设备就会发出声音，运行状态发生改变时，其发出的声音也随之改变。同时，电抗器内部的绕组和铁芯承担着电磁交换的重要功能，在高电压和强电磁的环境下可能发生不同的故障，从而运行的声音也会随之改变。
目前对电抗器的保护方法主要通过故障时电压、电流等电气参量进行继电保护，这种方法会降低用电质量以及提高供电公司的经济成本。并且相关潜伏性故障则难以检测，缺乏有效的在线检测方法、技术与装置。
发明内容
本发明的目的是提供一种基于近场宽带波束形成的大型电抗器故障检测方法及系统，能够检测到潜伏性故障，提前给相关人员预警来处理设备的状况，不仅提高人们的用电质量而且又能降低经济成本。
本发明的目的是通过以下技术方案实现的：
一种基于近场宽带波束形成的大型电抗器故障检测方法，包括：
在电抗器外部设置麦克风阵列；
利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描，将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像；
比较预设正常工作情况下与检测时的声音强度图像或声音能量图像，根据比较结果判断电抗器是否存在故障，并在判定发生故障时，发出预警。
一种基于近场宽带波束形成的大型电抗器故障检测系统，用于前述的方法，包括：
麦克风阵列，设置在电抗器外部；
采样处理装置，利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描，并获取麦克风阵列的接收信号；
检测处理装置，用于将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像；以及，将计算得到的电抗器的声音强度图像或声音能量图像，与预先获取的电抗器正常工作情况下的声音强度图像或声音能量图像对应的进行比较，根据比较结果判断电抗器是否存在故障；
非正常状态处理装置，用于在判定发生故障时，发出预警。
由上述本发明提供的技术方案可以看出，通过比较电抗器的声音强度或者能量图像与其正常工作时的声音强度或者能量图像，判断电抗器设备是否存在故障需要工作人员维护，进而实现了电抗器设备工作状态的在线检测。通过近场宽带波束形成的声音强度或者能量图像来判断电抗器工作状态，是一种在线检测方案，不仅降低了维护的成本，而且还提高了人们用电的质量，给大型电抗器设备的工作状态检测提供了一种新的方法。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。
图1为本发明实施例提供的一种基于近场宽带波束形成的大型电抗器故障检测方法的流程图；
图2为本发明实施例提供的平面麦克风阵列的信号模型示意图；
图3为本发明实施例提供的实现方法的一种系统构成示意图；
图中各标记对应的部件为：1-麦克风阵列；2-电抗器；3-采样处理装置；4-检测处理装置；5-非正常状态处理装置；6-目标；7-麦克风阵列的支架；8-麦克风阵元；9-参考阵元。
具体实施方式
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。
本发明实施例提供一种基于近场宽带波束形成的大型电抗器故障检测方法，如图1所示，其主要包括：
步骤1、在电抗器外部设置麦克风阵列。
本发明实施例中，所述麦克风阵列与电抗器表面平行，麦克风阵列平面与电抗器表面之间的距离为R，将M×N个麦克风按照X轴方向间隔d1(水平方向)，Y轴方向(垂直方向)间隔d2进行放置，最终形成[(M-1)×d1]×[(N-1)×d2]大小的均匀面阵，其中坐标为(0，0)的阵元作为参考阵元。
步骤2、利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描，将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像。
本发明实施例中，采用近场宽带波束形成方法对工作中的电抗器进行扫描，从而对扫描平面进行全方位的扫描；首先要建立近场宽带接收信号模型：如图2所示，以麦克风阵列的坐标原点阵元为参考阵元，X轴和Y轴作为阵列的基准线；假设有K个信号源
Figure BDA0002463833050000031
分别以二维角度
Figure BDA0002463833050000032
入射到麦克风阵列上，其中
Figure BDA0002463833050000033
θk分别表示第k个信号源入射的方位角和俯仰角，t为入射时刻；假设第m行第n列阵元的坐标为(xm,yn)，则第k个信号入射到该阵元上相对于第m行第n列阵元的距离dmnk表示为：
Figure BDA0002463833050000034
其中，rks＝R/cosθk表示第k个信号源到参考阵元的距离，R表示麦克风阵列平面与电抗器表面之间的距离，rmn表示第m行第n列阵元到参考阵元的距离，m＝1,2,...,M，n＝1,2,...,N；
以时间长度T为一个统计时段，在一个统计时段T内认为声源的位置基本不变。将统计时段T内的麦克风阵列接收信号划分为L个子段，每个子段包含I个时间采样点，相邻子段间的重叠率为γ，0≤γ＜1；
之后，对麦克风阵列接收到的每个子段宽带信号进行I点的离散傅里叶变换得到I个子带的窄带阵列信号，其中I不一定为2的整数次幂，若I为2的整数次幂会用快速算法节约该过程的运算时间。对于第i个子带，第k个信号源的导向矢量表示为
Figure BDA0002463833050000041
其中，
Figure BDA0002463833050000042
c表示声速大小；i＝1,2,...,I，fi表示第i个子带的频率；
对每个子段的阵列接收信号进行I点的离散傅里叶变换表示为：
Xl(fi)＝Al(fi)Sl(fi)+Nl(fi),i＝1,2,...,I，l＝1,2,...,L
其中，
Figure BDA0002463833050000043
表示阵列接收信号的第l个子段经I点离散傅里叶变换得到的第i个子带窄带阵列信号，x的下标数字为麦克风的序号，例如，
Figure BDA0002463833050000044
表示第1个麦克风接收信号的第l个子段经I点离散傅里叶变换得到的第i个子带窄带信号；
Figure BDA0002463833050000045
表示信号源的第l个子段经I点离散傅里叶变换得到的第i个子带窄带信号，s的下标数字为信号源的序号，例如，
Figure BDA0002463833050000046
为第1个信号源的第l个子段经I点离散傅里叶变换得到的第i个子带窄带信号；
Figure BDA0002463833050000047
表示第l个子段阵列接收噪声经I点离散傅里叶变换得到的第i个子带窄带阵列噪声，n的下标数字为麦克风的序号，例如，
Figure BDA0002463833050000048
表示第1个麦克风的第l个子段接收噪声经I点离散傅里叶变换得到的第i个子带窄带阵列噪声；A(fi)＝[a1(fi,r1s),a2(fi,r2s),…,aK(fi,rKs)]为阵列第l个子段第i个子带信号的导向矢量矩阵，这是基于认为在一个统计时段T内声源的位置基本不变得出的，a的下标数字为信号源的序号，例如，a1(fi,r1s)为第i个子带，第1个信号源的导向矢量矩阵。
本发明实施例设计了以下加权系数，并以此计算声音强度图像或声音能量图像。
在一个统计时段T内，对于第l个子段的阵列接收信号，在第i个子带的频率fi上，对于近场宽带波束形成而言，在X轴和Y轴两个方向上进行从左至右从上至下的幅度加权(或者在方位角俯仰角二位方向上进行从左至右从上至下的幅度加权)；对于在X轴和Y轴进行的每一个幅度加权，通过不同类型的窗函数进行幅度加权，加权矩阵为：
Figure BDA0002463833050000049
其中
Figure BDA00024638330500000410
为第i个子频率中X轴方向上窗函数的加权值，
Figure BDA00024638330500000411
为第i个子频率中Y轴方向上窗函数的加权值，
Figure BDA00024638330500000412
与
Figure BDA00024638330500000413
中的每一项表示一个幅度的窗函数的加权值；同时，令
Figure BDA00024638330500000414
其中vec表示矩阵或者向量的拉直运算。
对于第l个子段的阵列接收信号，在第i个子带的频率fi上，计算声音强度图像的公式为：
Figure BDA0002463833050000051
其中| |表示复数的取模，
Figure BDA0002463833050000052
表示
Figure BDA0002463833050000053
的共轭转置；在一个统计时段T内，在(x,y)坐标上的声音强度图像q(x,y)为各子段各子带频率的声音强度之和，其表达式为
Figure BDA0002463833050000054
在一个统计时段T内，阵列扫描电抗器的完整角度区域
Figure BDA0002463833050000055
获得的电抗器声音强度图像为：
Figure BDA0002463833050000056
在一个统计时段T内，在(x,y)坐标上的声音能量图像
Figure BDA0002463833050000057
为各子段各子带频率的声音强度之和，其表达式为
Figure BDA0002463833050000058
在一个统计时段T内，阵列扫描电抗器的完整角度区域
Figure BDA0002463833050000059
获得的电抗器声音能量图像为：
Figure BDA00024638330500000510
基于上述原理，通过步进扫描的方位角
Figure BDA00024638330500000511
和俯仰角Δθ的方式实时获取电抗器工作时的声音强度图像或者声音能量图像；其中，下标X,Y为完整角度区域
Figure BDA00024638330500000512
中麦克风阵列在X轴和Y轴获取所有信号的数量大小，X,Y的大小和步进选取的
Figure BDA00024638330500000513
和Δθ有关，即与完整角度区域
Figure BDA00024638330500000514
有关系。
步骤3、比较预设正常工作情况下与检测时的声音强度图像或声音能量图像，根据比较结果判断电抗器是否存在故障，并在判定发生故障时，发出预警。
根据先验知识，可以预先获取的电抗器正常工作情况下统计时段T内的声音强度图像，记为
Figure BDA00024638330500000515
计算得到统计时段T内的电抗器的声音强度图像记为
Figure BDA00024638330500000516
比较完整角度区域二者之差的绝对值：
Figure BDA00024638330500000517
其中|| ||表示矩阵的范数，若
Figure BDA00024638330500000518
则认为角度区域
Figure BDA00024638330500000519
出现故障，并发出预警(例如，通过灯光、声音等方式对工作人员进行预警)；其中，
Figure BDA00024638330500000520
为预先设置的完整角度区域
Figure BDA00024638330500000521
的声音强度图像与
Figure BDA00024638330500000522
的差异门限。
或者，根据先验知识，可以预先获取的电抗器正常工作情况下统计时段T内的声音能量图像，记为
Figure BDA0002463833050000061
计算得到统计时段T内的电抗器的声音强度图像记为
Figure BDA0002463833050000062
比较完整角度区域二者之差的绝对值：
Figure BDA0002463833050000063
若
Figure BDA0002463833050000064
则认为角度区域
Figure BDA0002463833050000065
出现故障，并发出预警；其中，
Figure BDA0002463833050000066
为预先设置的完整角度区域
Figure BDA0002463833050000067
的声音能量图像与
Figure BDA0002463833050000068
的差异门限。
本发明另一实施例还提供一种基于近场宽带波束形成的大型电抗器故障检测系统，用于实现前述的方法，如图3所示，其主要包括：
麦克风阵列1包含M×N个阵元8(参考阵元9)，麦克风阵列的支架为7，设置在电抗器2外部，接收电抗器2发出的声音信号6；
采样处理装置3，利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描，并获取麦克风阵列的接收信号；具体来说：利用麦克风阵列1采集电抗器2发出的声音信号6；在一个统计时段T内，通过分子段和离散傅里叶变换，将麦克风阵列1采集电抗器2声音信号6转换为各子段各子带信号；
检测处理装置4，用于将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像；具体来说，对于各子段各子带信号，阵列扫描电抗器的完整角度区域的每一组角度，用预先设计的加权矢量做近场宽带波束形成；在一个统计时段内，计算电抗器2的声音强度图像或声音能量图像；以及将计算得到的电抗器的声音强度图像或声音能量图像，与预先获取的电抗器2正常工作情况下的声音强度图像或声音能量图像对应的进行比较，根据比较结果判断电抗器2是否存在故障；
非正常状态处理装置5，用于在判定发生故障时，发出预警。
该系统的具体技术细节在之前的实施例中已经进行了详细的说明，故不再赘述。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将系统的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解，上述实施例的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM，U盘，移动硬盘等)中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。

Claims (5)

1.一种基于近场宽带波束形成的大型电抗器故障检测方法，其特征在于，包括：
在电抗器外部设置麦克风阵列；
利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描，将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像；
比较预设正常工作情况下与检测时的声音强度图像或声音能量图像，根据比较结果判断电抗器是否存在故障，并在判定发生故障时，发出预警；
其中，将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像的方式包括：
在一个统计时段T内，对于第l个子段的阵列接收信号，在第i个子带的频率fi上，在X轴和Y轴进行的每一个幅度加权，通过不同类型的窗函数进行幅度加权，加权矩阵为：
Figure FDA0002996103670000011
其中
Figure FDA0002996103670000012
为第i个子频率中X轴方向上窗函数的加权值，
Figure FDA0002996103670000013
为第i个子频率中Y轴方向上窗函数的加权值，
Figure FDA0002996103670000014
与
Figure FDA0002996103670000015
中的每一项表示一个幅度的窗函数的加权值；同时，令
Figure FDA0002996103670000016
其中vec表示矩阵或者向量的拉直运算；
对于第l个子段的阵列接收信号，在第i个子带的频率fi上，计算声音强度图像的公式为：
Figure FDA0002996103670000017
其中| |表示复数的取模，
Figure FDA0002996103670000018
表示
Figure FDA0002996103670000019
的共轭转置，Xl(fi)表示阵列接收信号的第l个子段经I点离散傅里叶变换得到的第i个子带窄带阵列信号；在一个统计时段T内，在(x,y)坐标上的声音强度图像q(x,y)为各子段各子带频率的声音强度之和，其表达式为
Figure FDA00029961036700000110
在一个统计时段T内，阵列扫描电抗器的完整角度区域
Figure FDA00029961036700000111
获得的电抗器声音强度图像为：
Figure FDA00029961036700000112
在一个统计时段T内，在(x,y)坐标上的声音能量图像
Figure FDA00029961036700000113
为各子段各子带频率的声音强度之和，其表达式为
Figure FDA0002996103670000021
在一个统计时段T内，阵列扫描电抗器的完整角度区域
Figure FDA0002996103670000022
获得的电抗器声音能量图像为：
Figure FDA0002996103670000023
其中，下标X,Y为完整角度区域
Figure FDA0002996103670000024
中麦克风阵列在X轴和Y轴获取所有信号的数量大小；
基于上述原理，通过步进扫描的方位角
Figure FDA0002996103670000025
和俯仰角Δθ的方式实时获取电抗器工作时的声音强度图像或者声音能量图像。
2.根据权利要求1所述的一种基于近场宽带波束形成的大型电抗器故障检测方法，其特征在于，所述在电抗器外部设置麦克风阵列，包括：
麦克风阵列平面与电抗器表面之间的距离为R，将M×N个麦克风按照X轴方向间隔d1，Y轴方向间隔d2进行放置，最终形成[(M-1)×d1]×[(N-1)×d2]大小的均匀面阵，其中坐标为(0，0)的阵元作为参考阵元。
3.根据权利要求2所述的一种基于近场宽带波束形成的大型电抗器故障检测方法，其特征在于，利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描时，先建立近场宽带接收信号模型：以麦克风阵列的坐标原点阵元为参考阵元，X轴和Y轴作为阵列的基准线；假设有K个信号源
Figure FDA0002996103670000026
分别以二维角度
Figure FDA0002996103670000027
入射到麦克风阵列上，其中
Figure FDA0002996103670000028
θk分别表示第k个信号源入射的方位角和俯仰角，t为入射时刻；假设第m行第n列阵元的坐标为(xm,yn)，则第k个信号入射到该阵元上相对于第m行第n列阵元的距离dmnk表示为：
Figure FDA0002996103670000029
其中，rks＝R/cosθk表示第k个信号源到参考阵元的距离，rmn表示第m行第n列阵元到参考阵元的距离，m＝1,2,...,M，n＝1,2,...,N；
以时间长度T为一个统计时段，将统计时段T内的麦克风阵列接收信号划分为L个子段，每个子段包含I个时间采样点，相邻子段间的重叠率为γ，0≤γ＜1；
之后，对麦克风阵列接收到的每个子段宽带信号进行I点的离散傅里叶变换得到I个子带的窄带阵列信号，对于第i个子带，第k个信号源的导向矢量表示为：
Figure FDA0002996103670000031
其中，
Figure FDA0002996103670000032
c表示声速大小；i＝1,2,...,I，fi表示第i个子带的频率；
对每个子段的阵列接收信号进行I点的离散傅里叶变换表示为：
Xl(fi)＝Al(fi)Sl(fi)+Nl(fi),i＝1,2,...,I，l＝1,2,...,L
其中，
Figure FDA0002996103670000033
表示阵列接收信号的第l个子段经I点离散傅里叶变换得到的第i个子带窄带阵列信号，x的下标数字为麦克风的序号；
Figure FDA0002996103670000034
表示源信号的第l个子段经I点离散傅里叶变换得到的第i个子带窄带信号，s的下标数字为信号源的序号；
Figure FDA0002996103670000035
表示第l个子段阵列接收噪声经I点离散傅里叶变换得到的第i个子带窄带阵列噪声，n的下标数字为麦克风的序号；A(fi)＝[a1(fi,r1s),a2(fi,r2s),…,aK(fi,rKs)]为阵列第l个子段第i个子带信号的导向矢量矩阵，a的下标数字为信号源的序号。
4.根据权利要求1所述的一种基于近场宽带波束形成的大型电抗器故障检测方法，其特征在于，所述比较预设正常工作情况下与检测时的声音强度图像或声音能量图像，根据比较结果判断电抗器是否存在故障，并在判定发生故障时，发出预警，包括：
根据先验知识，可以预先获取的电抗器正常工作情况下统计时段T内的声音强度图像，记为
Figure FDA0002996103670000036
计算得到统计时段T内的电抗器的声音强度图像记为
Figure FDA0002996103670000037
比较完整角度区域二者之差的绝对值：
Figure FDA0002996103670000038
其中||||表示矩阵的范数，若
Figure FDA0002996103670000039
则认为角度区域
Figure FDA00029961036700000310
出现故障，并发出预警；其中，
Figure FDA00029961036700000311
为预先设置的完整角度区域
Figure FDA00029961036700000312
的声音强度图像与
Figure FDA00029961036700000313
的差异门限；
或者，根据先验知识，可以预先获取的电抗器正常工作情况下统计时段T内的声音能量图像，记为
Figure FDA00029961036700000314
计算得到统计时段T内的电抗器的声音强度图像记为
Figure FDA00029961036700000315
比较完整角度区域二者之差的绝对值：
Figure FDA00029961036700000316
若
Figure FDA00029961036700000317
则认为角度区域
Figure FDA00029961036700000318
出现故障，并发出预警；其中，
Figure FDA00029961036700000319
为预先设置的完整角度区域
Figure FDA0002996103670000041
的声音能量图像与
Figure FDA0002996103670000042
的差异门限。
5.一种基于近场宽带波束形成的大型电抗器故障检测系统，其特征在于，用于实现权利要求1-4任一项所述的方法，包括：
麦克风阵列，设置在电抗器外部；
采样处理装置，利用麦克风阵列采用近场宽带波束形成方法对工作中的电抗器进行扫描，并获取麦克风阵列的接收信号；
检测处理装置，用于将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像；以及，将计算得到的电抗器的声音强度图像或声音能量图像，与预先获取的电抗器正常工作情况下的声音强度图像或声音能量图像对应的进行比较，根据比较结果判断电抗器是否存在故障；
非正常状态处理装置，用于在判定发生故障时，发出预警；
其中，将接收信号与预先设计的加权系数相结合，计算电抗器的声音强度图像或声音能量图像的方式包括：
在一个统计时段T内，对于第l个子段的阵列接收信号，在第i个子带的频率fi上，在X轴和Y轴进行的每一个幅度加权，通过不同类型的窗函数进行幅度加权，加权矩阵为：
Figure FDA0002996103670000043
其中
Figure FDA0002996103670000044
为第i个子频率中X轴方向上窗函数的加权值，
Figure FDA0002996103670000045
为第i个子频率中Y轴方向上窗函数的加权值，
Figure FDA0002996103670000046
与
Figure FDA0002996103670000047
中的每一项表示一个幅度的窗函数的加权值；同时，令
Figure FDA0002996103670000048
其中vec表示矩阵或者向量的拉直运算；
对于第l个子段的阵列接收信号，在第i个子带的频率fi上，计算声音强度图像的公式为：
Figure FDA0002996103670000049
其中| |表示复数的取模，
Figure FDA00029961036700000410
表示
Figure FDA00029961036700000411
的共轭转置，Xl(fi)表示阵列接收信号的第l个子段经I点离散傅里叶变换得到的第i个子带窄带阵列信号；在一个统计时段T内，在(x,y)坐标上的声音强度图像q(x,y)为各子段各子带频率的声音强度之和，其表达式为
Figure FDA00029961036700000412
在一个统计时段T内，阵列扫描电抗器的完整角度区域
Figure FDA00029961036700000413
获得的电抗器声音强度图像为：
Figure FDA0002996103670000051
在一个统计时段T内，在(x,y)坐标上的声音能量图像
Figure FDA0002996103670000052
为各子段各子带频率的声音强度之和，其表达式为
Figure FDA0002996103670000053
在一个统计时段T内，阵列扫描电抗器的完整角度区域
Figure FDA0002996103670000054
获得的电抗器声音能量图像为：
Figure FDA0002996103670000055
其中，下标X,Y为完整角度区域
Figure FDA0002996103670000056
中麦克风阵列在X轴和Y轴获取所有信号的数量大小；
基于上述原理，通过步进扫描的方位角
Figure FDA0002996103670000057
和俯仰角Δθ的方式实时获取电抗器工作时的声音强度图像或者声音能量图像。
CN202010327774.0A 2020-04-23 2020-04-23 基于近场宽带波束形成的大型电抗器故障检测方法及系统 Active CN111413649B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010327774.0A CN111413649B (zh) 2020-04-23 2020-04-23 基于近场宽带波束形成的大型电抗器故障检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010327774.0A CN111413649B (zh) 2020-04-23 2020-04-23 基于近场宽带波束形成的大型电抗器故障检测方法及系统

Publications (2)

Publication Number Publication Date
CN111413649A CN111413649A (zh) 2020-07-14
CN111413649B true CN111413649B (zh) 2021-07-06

Family

ID=71492087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010327774.0A Active CN111413649B (zh) 2020-04-23 2020-04-23 基于近场宽带波束形成的大型电抗器故障检测方法及系统

Country Status (1)

Country Link
CN (1) CN111413649B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112345855B (zh) * 2020-10-30 2021-10-08 浙江大学 一种并联电抗器运行状态检测系统及方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60100060A (ja) * 1983-11-04 1985-06-03 Electric Power Dev Co Ltd 部分放電検出装置
JPH0217465A (ja) * 1988-07-06 1990-01-22 Toyo Electric Mfg Co Ltd 電気機械巻線の絶縁破壊事前検知方法
CN100466061C (zh) * 2005-08-15 2009-03-04 华为技术有限公司 一种宽带波束形成方法和装置
US8503291B1 (en) * 2009-08-13 2013-08-06 Marvell International Ltd. Systems and methods for directing a beam towards a device in the presence of interference based on reciprocity
CN101907678B (zh) * 2010-07-12 2012-09-05 西安西电科大电力科技有限公司 电缆故障测试系统及故障点的确定方法
CN103149513B (zh) * 2013-02-26 2015-06-24 中国科学院电工研究所 一种变压器局部放电超声源重建的定位方法及装置
CN103438990B (zh) * 2013-09-10 2016-03-02 北京科技大学 一种高温环境下声场测量与分析装置及其分析方法
CN103969560B (zh) * 2014-05-28 2016-10-26 北达正视(北京)科技有限公司 一种用于局部放电故障检测的可视化超声波检测系统
CN105223544B (zh) * 2015-08-26 2018-01-12 南京信息工程大学 近场线性约束最小方差自适应加权频率不变波束形成方法
CN106023996B (zh) * 2016-06-12 2019-08-27 杭州电子科技大学 基于十字形声阵列宽带波束形成的声识别方法
CN106934421B (zh) * 2017-03-16 2020-11-06 山东大学 基于2dpca及svm的电力变压器故障检测方法
CN109001678B (zh) * 2017-06-06 2021-10-01 中国科学技术大学 一种基于三维麦克风阵列的雷声检测与定位方法
CN110987402B (zh) * 2019-12-20 2022-01-11 中国科学技术大学 基于声音能量的球磨机工作状态检测方法

Also Published As

Publication number Publication date
CN111413649A (zh) 2020-07-14

Similar Documents

Publication Publication Date Title
JP3100757B2 (ja) 監視診断装置
CN112526451B (zh) 基于麦克风阵列成像的压缩波束形成及系统
CN112180329B (zh) 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法
CN104614653A (zh) 基于天线阵列的局放检测装置的多局放点定位和分辨方法
CN111413649B (zh) 基于近场宽带波束形成的大型电抗器故障检测方法及系统
CN106066468A (zh) 一种基于声压、振速互谱法的矢量阵左右舷分辨方法
CN203037774U (zh) 一种应用在变压器及电抗器上的噪声源定位装置
CN111582355B (zh) 基于最优宽带波束形成的大型电抗器故障检测方法及系统
CN110687506A (zh) 一种基于矢量传声器阵列的低频噪声源定位装置及方法
CN105785349A (zh) 一种相控阵三维声学摄像声呐的噪声去除方法
CN113064129A (zh) 一种高频地波雷达海流合成方法
WO2023207047A1 (zh) 一种高分辨率谱估计声阵列成像方法
CN102546049B (zh) 天线指标计算方法及装置
CN107992801B (zh) 一种振动设备的故障诊断方法及系统
CN111812581A (zh) 基于原子范数的球面阵列声源波达方向估计方法
CN115542091A (zh) 一种无人机局放检测系统及方法
CN103048641A (zh) 具有联合约束优化形式的矢量阵稳健聚焦处理方法
CN112098938A (zh) 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法
CN116818321A (zh) 基于超声相控阵的轴承故障检测方法、装置及存储介质
CN115453300A (zh) 一种基于声传感器阵列的局部放电定位系统及方法
CN111965628B (zh) 一种垂直出水航行体瞬时波浪参数估计方法
CN112612027B (zh) 一种浅海环境下利用声能量起伏的海洋内波监测方法
CN114706084A (zh) 一种基于yolo算法的水声被动跟踪多目标的方法
Han et al. Distributed and parallel subarray beamforming for underwater real-time 3-D acoustic imaging
CN117572435B (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