基于变压器故障噪声监测诊断技术的音频特征提取方法研究

摘

要：变压器在运行过程中可能会出现不同类型的故障，导致变压器发出异常声响。通过采集变压器故障运行时的异常音频，并提取音频特征量进行分析处理，就可判断出变压器的故障类型。此种变压器故障诊断方法的关键点在于采集变压器运行中的故障音频信号并提取其特征量，尤其是从包含环境噪声的混合音频中准确提取有效的故障特征量是该方法的技术核心。现通过研究Hilbert-Huang变换在音频分析中的应用，探讨了变压器故障音频特征分离和提取的有效方法。关键词：变压器故障诊断；音频特征；Hilbert-Huang变换0

引言变压器是电力系统中的重要设备，其运行稳定性直接影响供电质量，因此变压器运行状态监测及故障诊断已成为电网维护的一项重要工作。目前，传统的变压器监测一般需要使用接触式测量仪器进行巡检，人力成本高，效率低。基于此，本文提出了一种以音频分析为基础的新检测手段：根据变压器运行时声音的变化来反映变压器运行状态的变化，即在变压器发生故障时提取音频特征信息，使用经验模态分解（EMD）方法将信号分解成若干IMF特征分量并进行Hilbert变换，得到相应的Hilbert谱，最后通过科学算法和软件分析将异常运行时的音频特征量与正常运行时的特征量做对比分析，从而实现变压器故障诊断功能。在利用音频特征监测变压器状态时，对音频信号的处理是其核心步骤，主要包括音频特征的提取及故障特征的识别与判断。由于变压器型号种类多、结构复杂、工作环境多样，因此在实际采样过程中采集到的音频除了含有变压器自身固有机械振动产生的正常音频以外，还有许多环境噪声干扰混杂于其中，需要对采集的音频信号做降噪处理，从中分离并提取有效的音频特征。本文就如何分离并提取音频特征进行探讨分析。

1

Hilbert-Huang变换简介Hilbert-Huang变换是由黄锷等人提出的一种新的信号分析法。该分析方法是在经验模态分解（EMD）方法的基础上，结合Hilbert谱的概念和分析改进Hilbert谱之后得出的方法，被命名为Hilbert-HuangTransform（简称HHT），即希尔伯特黄变换。在实际工程中对于变电站设备故障噪声音频特征的提取，一般多采用梅尔频率倒谱系数（MFCC）方法。该方法是一种频域分析方法，它模拟了人耳非线性听觉特性，并且基于梅尔刻度对频带进行等距划分，使得间隔频带更加接近于人类的听觉，故其被广泛运用于各类声音信号的处理中。但这种提取方法在加性噪声的情况下稳定性和精度并不高，所以还需要进一步改进，目前采用较多的改进方案是将其与隐马尔可夫模型结合运用。Hilbert-Huang变换是一种时域分析方法。首先利用经验模态分解（EMD）方法将一个混合信号分解成若干IMF分量，通过Hilbert变换得到各个IMF的Hilbert谱，将每个IMF分量的Hilbert谱分解成独立的子向量，最后汇总所有IMF的Hilbert谱，合成原始信号的Hilbert谱，从而达到分离混合信号的目的。同经典的傅里叶变换信号分析方法相比，EMD方法克服了用傅里叶变换方法提取信号频谱时必须要利用信号的全部时域信息这一缺点，具有良好的时频分辨率和自适应性，能够完美地重构原始信号。Hilbert-Huang变换适用于对非线性、非平稳信号的分析处理，具有完全自适应性，可以在时间和频率上同时达到很高的精度，可用于突变信号的分析。这种方法非常适用于变电站背景噪声强的混合音频处理，能够从复杂的背景噪声中分离出有效的音频信号。

2

Hilbert-Huang变换的应用2.1

经验模态分解（EMD）Hilbert-Huang变换首先要进行的就是通过EMD方法将信号分解，这个信号可以是平稳的、线性的，也可以是不平稳的、非线性的，EMD的主要作用是将信号分解成若干个IMF，每一个IMF必须满足两个条件：（1）信号的极值点（极大值或极小值）数目和过零点数目相等或最多相差一个；（2）由局部极大值构成的上包络线和由局部极小值构成的下包络线的平均值为零。其过程如下：首先获得给定的信号的所有极点，用三次样条插值函数求出上下包络线，上下包络线的均值记为，与的差记为，则：由于通常情况下此时的还不是一个IMF分量序列，这就需要对重复上述操作，直到过零点的数目和极点的数目相等或仅相差一个，且各个瞬时平均值都等于零，从而得出第一个IMF分量。然后用减去，得到残余信号，再用重复上述操作，得到第二个IMF分量、第三个IMF分量……，最后得到第个IMF分量和残余信号，满足给定的终止条件时结束（足够小或者是单调函数）。最终获得的表达式为：2.2

Hilbert谱和Hilbert边际谱对于每一个分解出来的IMF分量，做Hilbert变换：生成的解析信号为：由此得到瞬时振幅和瞬时相位的函数：进一步可求出瞬时频率：将所有分量的Hilbert谱汇总，就得到了：这个表达式省略了残余函数，它可能是一个单调函数，也可能是一个常数。表示取实部，展开式称为Hilbert谱，记为：由此对上式进行积分就可以进一步得到Hilbert边际谱：3

仿真实验实验采用两组信号源，其中一组为多周期混合信号线性叠加的模拟信号函数：另外一组为一段实地采样音频，为混合音频信号，该音频信号同时包含随机环境噪声以及具有固定周期和规律的特定信号，对两组信号分别进行Hilbert-Huang变换。首先对模拟函数信号进行分解，分解结果如图1所示。其中，图1（a）为函数的原始信号，图1（b）为函数原始信号的频谱。经过Hilbert-Huang变换后，通过EMD方法将函数信号分解，按频率由高到低排列得：IMF1为分解得到的特征分量；IMF2为分解得到的分量；IMF3为分解得到的分量；IMF4为分解后的残余项。函数IMF分量如图2所示，图2（a）（b）（c）（d）为各IMF分量波形图，图2（e）（f）（g）（h）为IMF分量对应的频谱。以IMF2特征分量为例，对函数信号进行Hilbert-Huang变换后，得到能够反映IMF2分量特征的物理量。IMF2的特征波形及其频谱、包络线、瞬时频率、调制信号和调制信号的频谱如图3所示。在完成对多周期混合信号线性叠加模拟信号的分析之后，再对含有大量随机噪声的实地采样音频信号进行处理，根据原始信号绘制波形图及其对应的频谱，如图4所示。

经过Hilbert-Huang变换的EMD分解后，混合音频被成功分解成15个IMF分量，其中前4个频率最高，所含信息量相对较大的4个IMF分量的波形图如图5（a）（b）（c）（d）所示，其对应的频谱如图5（e）（f）（g）（h）所示。对分解求得的15组IMF分量分别进行Hilbert-Huang变换的EMD分解后，便可提取出每一个IMF分量的特征波形及其频谱、包络线、瞬时频率、调制信号和调制信号的频谱等可以反映此分量特征信息的特征量，通过将这些特征量与变压器正常运行时的音频特征量进行对比分析，即可及时判断变压器的运行状态或者故障类型。

4

结语Hilbert-Huang变换在不平稳、非线性信号的分析处理上体现出来的优越性，使得这种方法在地震波分析、机械故障诊断、语音信号处理等多