1.一种变压器绕组的机械状态检测方法,其特征在于,包括如下步骤:步骤S1,在变压器箱壁布置振动加速度传感器,以预定的采样频率采集变压器振动信号;步骤S2,通过对振动信号傅里叶谱进行划分,构造自适应小波函数,得到具有紧支撑频率的窄带子信号,将原始振动信号分解成多个固有模式函数之和;步骤S3,对振动信号的固有模式函数进行Hilbert变换得到振动信号的时频谱图;步骤S4,对振动信号时频谱图提取局部二值特征,构造振动信号特征向量;包括步骤S4a、S4b、S4c、S4d、S4e:步骤S4a、由于变压器振动信号主要频率在500Hz以下,所以关注时频谱图500Hz范围内;步骤S4b、将时频谱时间轴均匀划分为m段,频率轴均匀划分为n段,得到m乘以n个网格,将每个网格内的时频幅值相加,记为h kj ,时频谱网格记为H(k,j);步骤S4c、不考虑H(k,j)边界元素,对每个非边界元素,以其为中心,取3x3大小窗口,以窗口中心点的幅值为阈值,对该中心点8个邻域幅值做二值化处理,即窗口内其它位置的幅值分别与窗口中心的幅值进行比较,当大于等于中心幅值时,其对应位置赋值为1,否则赋值为0;步骤S4d、对领域赋值加权求和,作为中心元素的LBP值:式中,H c 为所关注的时频谱网格元素的幅值,H i 为所关注元素邻域的幅值;步骤S4e、统计时频网格的LBP特征,即统计每个LBP值的元素的个数,将统计结果按LBP值大小从小到大排列,形成一个16维的向量V LBP =[LBP 0 ,LBP 2 ,…,LBP 15 ],其中LBP i 表示在时频网格中,LBP值为i的网格的个数,对V LBP 归一化,形成16维的特征向量V;步骤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],并按递减规律排列(M 1 ≥M 2 ≥...M M );步骤S2d、定义阈值为M M +α(M 1 -M M ),比较极大值点幅值是否大于阈值,若大于,保留,否则删除,所保留的极值点个数记为N;步骤S2e、分割傅里叶谱,划分边界为相邻极值点的中点,且规定0和π为第1条和最后一条划分边界,则每段频谱可表示为:Λ n =[ω n-1 ,ω n ],n=1,2,…,N式中,Λ n 为第n段频谱,∪为并集符号,ω n-1 ,ω n 为第n段频谱的边界;步骤S2f、构造经验小波函数 和经验尺度函数 如下式中:β(x)=x 4 (35-84x+70x 2 -20x 3 )式中,β为过渡函数,γ为频段系数;步骤S2g、定义经验小波变换,构造细节系数为经验小波函数与信号的内积,构造近似系数为经验尺度函数与信号的内积:式中, 代表细节系数, 代表近似系数,x(t)代表原信号,F -1 代表逆傅立叶变换;步骤S2h、根据经验小波变换构造固有模式函数:经过上述步骤,原始振动信号x(t)被分解为若干个固有模式函数之和:
3.如权利要求2所述的方法,其特征在于,所述步骤S5进一步包括:根据下述公式计算所述基准特征向量和待测信号的特征向量的相似度d:其中,V 1 表示基准振动信号的特征向量,V 2 表示待检测振动信号的特征向量。
4.根据权利要求3所述的方法,其特征在于,在所述步骤S6中,若所述相似度d大于或等于预定阈值,则判定所述变压器绕组的机械状态为正常;若所述相似度d小于预定阈值,则判定所述变压器绕组的机械状态为异常。
5.根据权利要求4所述的方法,其特征在于,所述预定阈值为0.8。