一种变压器绕组的机械状态检测方法技术

本发明专利技术公开了一种变压器绕组的机械状态检测方法，包括如下步骤：步骤S1，在变压器箱壁布置振动加速度传感器，以预定的采样频率采集变压器振动信号；步骤S2，将原始振动信号分解成多个固有模式函数之和；步骤S3，对振动信号的固有模式函数进行Hilbert变换得到振动信号的时频谱图；步骤S4，对振动信号时频谱图提取局部二值特征，构造振动信号特征向量；步骤S5,获得基准特征向量以及待测振动信号的特征向量，并计算两者的相似度；步骤S6，将所述相似度与一预定阈值进行比较，确定所述变压器绕组机械状态。实施本发明专利技术，可以通过在线监测变压器突发短路时变压器箱壁振动信号对绕组的工作状态进行判别。

【技术实现步骤摘要】
一种变压器绕组的机械状态检测方法
本专利技术涉及电力设备状态检测方法，尤其是涉及一种变压器绕组的机械状态检测方法。
技术介绍
变压器是电力系统最重要的设备之一，它的安全可靠运行影响着整个电力系统，因此，研究其状态在线监测方法是智能电网建设的重要课题，具有重要的经济效益和社会效益。统计数据表明，变压器的绕组是发生故障最多的部件之一，变压器出口短路事故时有发生，由冲击电流引起的强大电磁力作用在绕组上，引起绕组强烈振动，造成绕组机械结构损坏。传统的短路电抗法和频响分析法灵敏度不高，仅在变压器绕组发生严重变形时能够给出准确的判断。目前，针对变压器绕组故障的带电检测已有较多研究，如红外检测、局放检测等，这些方法在变压器在线监测方面各有优势，受到较多关注。但红外检测和局放检测主要针对电接触缺陷和绝缘缺陷，无法对机械故障进行有效诊断。近年来，振动法检测变压器绕组机械状态逐渐受到关注，一方面，对振动信号的采集简单方便易操作，非常适用于在线监测，广泛应用于大型机械故障诊断技术中；另一方面，振动信号包含大量能反映变压器机械状态的信息，因此振动分析法灵敏度高，能检测出变压器绕组的松动、微小变形等前期故障，及时对变压器进行检修。其原理本质，是因为变压器绕组不同状态时，其机械状态发生改变，必然会导致振动特性的不同，因此可以根据振动信号的特征，判断变压器绕组是否发生故障。振动传感器非接触式信号采集方式不会对其正常工作造成影响，因此具有极好的应用前景。应用振动法对绕组状态检测，其重点和难点在于振动信号的分析方法以及状态判据，目前尚无统一标准。变压器短路冲击所产生的振动信号，具有强时变和非平稳性，信号特征不易提取，时域和频域特征难以表征信号更细致变化，无法灵敏反应变压器状态。近年来，非平稳信号分析方法快速发展，兼顾时域和频域特征的时频谱分析方法在振动信号处理领域取得了很多成果，如Hilbert-Huang变换、小波分析法以及一些改进算法都成功运用到变压器振动信号处理上，但这些算法也存在自适应性不够、时频分解结果不够准确等缺点。
技术实现思路
本专利技术所要解决的技术问题在于提供一种变压器绕组状态诊断方法，通过在线监测变压器突发短路时变压器箱壁振动信号对绕组的工作状态进行判别。为了解决上述技术问题，本专利技术所采用的技术方案为，提供一种变压器绕组的机械状态检测方法，包括如下步骤：步骤S1，在变压器箱壁布置振动加速度传感器，以预定的彩样频率采集变压器振动信号；步骤S2，通过对振动信号傅里叶谱进行划分，构造自适应小波函数，得到具有紧支撑频率的窄带子信号，将原始振动信号分解成若干个固有模式函数之和；步骤S3，对振动信号的固有模式函数进行Hilbert变换得到振动信号的时频谱图；步骤S4，对振动信号时频谱图提取局部二值特征，构造振动信号特征向量；步骤S5,根据上述步骤以绕组正常状态下的振动信号计算特征向量，记为基准特征向量；传感器安装位置不变，采集变压器突发短路时振动信号，根据上述步骤计算特征向量，作为待测振动信号的特征向量，根据所述基准特征向量和待测信号的特征向量，计算两者的相似度；步骤S6，将所述相似度与一预定阈值进行比较，确定所述变压器绕组机械状态。优选地，步骤S2具体包括如下步骤：步骤S2a、使用傅里叶变换对振动信号x(t)处理，得到信号的频域表示X(w)，并将频率范围归一化到[0,π]；步骤S2b、检测频域内的所有极大值点:对频谱函数取导，得到一次导数X'(ω)；对一次导数取符号函数，当X'(ω)>0，X'(ω)＝1，当X'(ω)<0，X'(ω)＝-1，当X'(ω)＝0，X'(ω)＝0；对符号函数取导，得到频谱函数二次导X”(ω)，若X”(ω)＝-2，则X(w+1)为极大值点。步骤S2c、记为频谱函数范围内检测到的极大值点幅值，将其归一化到[0,1]，并按递减规律排列(M1≥M2≥...MM)；步骤S2d、定义阈值为MM+α(M1-MM)，比较极大值点幅值是否大于阈值，若大于，保留，否则删除，所保留的极值点个数记为N；步骤S2e、分割傅里叶谱，划分边界为相邻极值点的中点，且规定0和π为第1条和最后一条划分边界，则每段频谱可表示为:Λn＝[ωn-1,ωn],n＝1,2,L,N式中，Λn为第n段频谱，U为并集符号，ωn-1,ωn为第n段频谱的边界；步骤S2f、构造经验小波函数和经验尺度函数如下式中：β(x)＝x4(35-84x+70x2-20x3)式中，β为过渡函数，γ为频段系数；步骤S2g、定义经验小波变换，构造细节系数为经验小波函数与信号的内积，构造近似系数为经验尺度函数与信号的内积：式中，代表细节系数，代表近似系数，x(t)代表原信号，F-1代表逆傅立叶变换；步骤S2h、根据经验小波变换构造固有模式函数：经过上述步骤，原始振动信号x(t)被分解为若干个固有模式函数之和：优选地，所述步骤S4进一步包括如下步骤：步骤S4a、由于变压器振动信号主要频率在500Hz以下，所以关注时频谱图500Hz范围内；步骤S4b、将时频谱时间轴均匀划分为m段，频率轴均匀划分为n段，得到m乘以n个网格，将每个网格内的时频幅值相加，记为hkj，时频谱网格记为H(k,j)；步骤S4c、不考虑H(k,j)边界元素，对每个非边界元素，以其为中心，取3x3大小窗口，以窗口中心点的幅值为阈值，对该中心点8个邻域幅值做二值化处理，即窗口内其它位置的幅值分别与窗口中心的幅值进行比较，当大于等于中心幅值时，其对应位置赋值为1，否则赋值为0；步骤S4d、对领域赋值加权求和，作为中心元素的LBP值:式中，Hc为所关注的时频谱网格元素的幅值，Hi为所关注元素邻域的幅值；步骤S4e、统计时频网格的LBP特征，即统计每个LBP值的元素的个数，将统计结果按LBP值大小从小到大排列，形成一个16维的向量VLBP＝[LBP0,LBP2,…,LBP15]，其中LBPi表示在时频网格中，LBP值为i的网格的个数，对VLBP归一化，形成16维的特征向量V。优选地，所述步骤S5进一步包括：根据下述公式计算所述基准特征向量和待测信号的特征向量的相似度d：其中，V1表示基准振动信号的特征向量，V2表示待检测振动信号的特征向量。优选地，在所述步骤S6中，若所述相似度d大于或等于预定阈值，则判定所述变压器绕组的机械状态为正常；若所述相似度d小于预定阈值，则判定所述变压器绕组的机械状态为异常。优选地，所述预定阈值为0.8。实施本专利技术实施例，具有如下的有益效果：本专利技术提出一种结合经验小波分解和局部二值编码对变压器振动信号进行时频谱特征提取，在信号处理的自适应性和准确性方面做出了改进，提高了绕组机械状态检测结果的可靠性。本专利技术的实施例，通过对变压器突发短路时的振动信号分析，可及时、有效检测绕组机械状态。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。图1为本专利技术提供的一种变压器绕组的机械状态检测方法；图2为本专利技术一个实施例中第一次突发短路时的振动波形；图3为本专利技术一个实施本文档来自技高网...

【技术保护点】
1.一种变压器绕组的机械状态检测方法，其特征在于，包括如下步骤：步骤S1，在变压器箱壁布置振动加速度传感器，以预定的采样频率采集变压器振动信号；步骤S2，通过对振动信号傅里叶谱进行划分，构造自适应小波函数，得到具有紧支撑频率的窄带子信号，将原始振动信号分解成多个固有模式函数之和；步骤S3，对振动信号的固有模式函数进行Hilbert变换得到振动信号的时频谱图；步骤S4，对振动信号时频谱图提取局部二值特征，构造振动信号特征向量；步骤S5,根据上述步骤以绕组正常状态下的振动信号的基准特征向量，以及获得变压器突发短路时振动信号对应的待测振动信号的特征向量，并计算所述基准特征向量和待测信号的特征向量相似度；步骤S6，将所述相似度与一预定阈值进行比较，确定所述变压器绕组机械状态。

【技术特征摘要】
1.一种变压器绕组的机械状态检测方法，其特征在于，包括如下步骤：步骤S1，在变压器箱壁布置振动加速度传感器，以预定的采样频率采集变压器振动信号；步骤S2，通过对振动信号傅里叶谱进行划分，构造自适应小波函数，得到具有紧支撑频率的窄带子信号，将原始振动信号分解成多个固有模式函数之和；步骤S3，对振动信号的固有模式函数进行Hilbert变换得到振动信号的时频谱图；步骤S4，对振动信号时频谱图提取局部二值特征，构造振动信号特征向量；步骤S5,根据上述步骤以绕组正常状态下的振动信号的基准特征向量，以及获得变压器突发短路时振动信号对应的待测振动信号的特征向量，并计算所述基准特征向量和待测信号的特征向量相似度；步骤S6，将所述相似度与一预定阈值进行比较，确定所述变压器绕组机械状态。2.如权利要求1所述的方法，其特征在于，步骤S2具体包括如下步骤：步骤S2a、使用傅里叶变换对振动信号x(t)处理，得到信号的频域表示X(w)，并将频率范围归一化到[0,π]；步骤S2b、检测频域内的所有极大值点:对频谱函数取导，得到一次导数X'(ω)；对一次导数取符号函数，当X'(ω)>0，X'(ω)＝1，当X'(ω)<0，X'(ω)＝-1，当X'(ω)＝0，X'(ω)＝0；对符号函数取导，得到频谱函数二次导X”(ω)，若X”(ω)＝-2，则X(w+1)为极大值点。步骤S2c、记为频谱函数范围内检测到的极大值点幅值，将其归一化到[0,1]，并按递减规律排列(M1≥M2≥...MM)；步骤S2d、定义阈值为MM+α(M1-MM)，比较极大值点幅值是否大于阈值，若大于，保留，否则删除，所保留的极值点个数记为N；步骤S2e、分割傅里叶谱，划分边界为相邻极值点的中点，且规定0和π为第1条和最后一条划分边界，则每段频谱可表示为:Λn＝[ωn-1,ωn],n＝1,2,L,N式中，Λn为第n段频谱，U为并集符号，ωn-1,ωn为第n段频谱的边界；步骤S2f、构造经验小波函数和经验尺度函数如下式中：β(x)＝x4(35-84x+70x2-20x3)式中，β...

【专利技术属性】
技术研发人员：李勋，艾精文，黄荣辉，李舒适，王丰华，
申请(专利权)人：深圳供电局有限公司，上海交通大学，
类型：发明
国别省市：广东,44