PCI1714在提高水轮机空化噪声频谱分析精度的方法中的应用.doc

Application Story项目介绍：水轮机空化空蚀是造成机组运行不稳定的主要因素之一，其空化振动噪声常表现为非平稳性，采样非同步误差以及加窗截断必将会导致频谱能量泄漏、使频谱分析结果存在误差，进而影响机组空化诊断的准确性。本文在对空化噪声产生机理及其特征分析的基础上，兼顾工程中算法实现的可能性以及空化监测系统的分析实时性，通过选择Welch频谱分析方法、选取Hanning窗函数，并采用双峰谱线频谱校正算法等措施，以此提高空化噪声的频谱分析的准确性。文中给出了频谱幅值修正的算法，并列举了部分试验结果。案例名称：研华高速采集板卡CI-1714在电力行业提高水轮机空化噪声频谱分析精度的方法研究中的应用行业分类：电力地点：n 系统需求按照GB/T154691995对空化(Cavitation)的定义：空化是当流道中局部压力下降至临界压力(一般接近汽化压力)时，水中气核成长为气泡，气泡的聚积、流动、分裂、溃灭过程的总称。由于空化造成的过流部件材料损坏称为空蚀(Cavitation Erosion）。水轮机中的空化、空蚀是造成水轮机能量特性下降、运行不稳定的主要原因之一，其严重程度已经成为一般水电厂大修周期的决定性因素。由于空化、空蚀常常会引起水轮机内的噪声增强、振动加剧以及机组效率急速下降等现象，因此噪声法、振动法成为目前监测水轮机空化空蚀的主要手段1。由于空泡溃灭过程的持续时间极短，变化又非常剧烈，水轮机空化噪声呈现出脉冲性的振动声波，表现出非平稳性性。虽然小波变换可以根据信号不同的频率成分，自动调节采样密度，实现在时域和频域同时具有良好的局部化特性，从而可以很好的处理空化噪声的突变等情况，但数据经小波变换后得到的是小波系数，而不是信号的频谱值，缺乏直观性，在频谱分析方面并不具备较高的优势。在空化噪声的频谱分析中，基于傅立叶变换()的频谱分析方法仍然有着不可比拟的重要性。但是，由于空化信号的非平稳性，等间距自由采样必存在非同步采样误差，产生频谱能量泄漏，难以获得高精度的测量结果。尤其在水轮机开/停机过程中，频谱泄漏更为严重，因为水轮机在开/停机过程中的振动空化声响可以说是典型的带有调频调幅特性的瞬态信号，进行傅立叶频谱分析，获得的往往是整个开/停机过程中的振动声响的平均效应，这种情况往往会掩盖水轮机的开机或停机过程中的瞬态特性。因此，本文采用选取合适的窗函数、采用合理的频谱分析方法、并进行频谱校正等方法从不同的角度实现测量误差的最小化。n 系统描述 1.水轮机空化振动噪声特性及其监测原理 空化噪声的频谱取决于两个因素：(1)单位时间内溃灭的汽泡的数目和溃灭的速度，(2)溃灭发生时所引起冲击压的传播介质。由于气核的产生、成长以及溃灭的随机性，导致空化噪声的宽频特性，不仅包含20Hz20kHz的可闻听的声波，又包含20kHz以上的超声波信号，表现为非平稳性，在水轮机的开/停机过程中，这种非平稳性会表现的尤为突出。水轮机空化、空蚀的过程与水轮机的运行状态以及尾水管的压力脉动有着密切的关联。在空化、空蚀监测单元中，需要获得水轮机的工况状态信息。这些状态信息可以直接从水轮机状态监测与诊断单元中采用CAN通信方式获得2。整个监测单元的架构如图1所示。在本系统中，采用4只高频振动加速度传感器LC0151TA以及宽频超声波接收器SR15，分别安装在1#和9#导叶连杆拐臂套筒、4#导叶附近的靠近水导轴承的支持盖上、尾水管检修门等地方。所选用的研华采集卡PCI1712、PCI1714分别以150kHz/Chan和2.5MHz/Chan的采样频率采集振动加速度和超声波信号。采集单元采用研华加固的壁挂式工业控制计算机IPC6606，完成数据的采集、分析以及与其他单元的通信任务。2.频谱泄漏的原因以及改进措施21 频谱泄漏的原因简单地说，造成频谱泄漏的原因是采样频率与信号频率不同步，造成采样信号的相位在始端与终端不连续而引起的。在频谱分析中，如果采样值不能均匀的分布在一个信号周期内时，进行离散快速傅立叶变换()时，在频域内将会引起频谱泄露。从时域上看，在对非平稳的测量信号进行变换时，相当于对其进行周期延拓，由于采样非同步误差的存在，导致信号在边界点处不连续(或其导数在边界点不连续) 给信号带来了附加的高频分量。从频域看，在对于广域的信号进行频谱分析时，要求信号必须是有限长的离散信号，需要对被分析的信号进行加窗截断，此时在频域内则表现为信号频谱与所加窗函数频谱的卷积，信号的频谱将会受到所加窗频谱的调制。由于空化信号的非平稳性，导致窗函数的过零点与信号离散频点不统一，而导致信号分析的不连续，产生频谱能量泄漏。22 改进措施与方法221 频谱分析的方法选择在设备状态监测和故障诊断中，从频域中提取信号的特征是很常见的方式，常用的方法是计算信号的幅值谱或者功率谱密度(Power Spectrum DensityPSD)，以此提取故障特征。其中功率谱密度是沿频率轴分布的，不仅能反映出信号的频谱结构，而且能够反映出各频率成分所占的比重和数值，结合功率谱分布结构，就可以判别出各主要特征信号源的频率和能量分布情况，进而为基于频谱分析的状态监测和故障诊断提供较为科学的依据。经典谱估计主要有直接法和间接法两种：直接法是对有限个样本数据进行计算得到功率谱的方法，又称周期图法；间接法是先得到样本数据的自相关函数估计，然后进行计算得到功率谱的方法，由于计算复杂，工程中一 般不予采用。对于直接法的功率谱估计，当数据长度太大时，谱曲线起伏加剧，若太小，谱的分辨率较差，必须改进其算法。改进的算法有Bartlett法和Welch法两种。本设计采用Welch法，主要原因是Welch法在分段时，可使各段之间有一定采集点的重叠，相对而言增加了分段数，使基于DFT的频谱分析具有一定的遗传性，而Bartlett法没有，因而Welch法得到的谱估计要比Bartlett法的方差要小。实际上，Welch法对Bartlett法进行了两方面的修正：(1)选择适当的窗函数，并在周期图计算前直接加窗，得到的每一段的周期图为：为归一化因子，M为窗函数的宽度；(2)在分段时，可使各段之间有一定重叠，使方差减少。 222 窗函数的选择在对广域的信号进行频谱分析时，要求信号必须是有限长的离散信号，须对被分析的信号进行加窗截断。加适当窗可有效地改善时域截断处的不连续状态。但是截断所带来的频谱泄漏是本身固有的缺陷，是影响谱分析与估计精度的重要因素。频谱能量泄漏程度主要是由所选窗谱的形状决定的，工程中降低频谱能量泄漏的主要方法是选择适合所处理信号的窗函数。窗函数的主瓣宽度越窄，则频率分辨率越高，且由主瓣导致的频谱能量泄漏的范围越小；相对旁瓣幅度越小，则旁瓣导致的泄漏程度越小。因此，为防止频谱泄漏，窗函数的选择一般需要满足三个基本条件：(1)由于主瓣宽度影响信号频谱的展宽程度，因此尽可能地减少窗函数的主瓣宽度；(2)尽可能增加窗函数的主瓣高度与旁瓣高度之比，使得旁瓣迅速衰减，利于减少频谱的混叠效应；(3)当时，窗函数应满足，即窗函数应当是选频函数。但这三个要求有时是相互制约的，所有的旁瓣抑制技术均是以牺牲谱分辨率为代价的，因此在选择窗函数时，要尽可能兼顾这三个要求。工程中，几种常用的比较简洁窗函数有矩形窗、Hanning窗，Hamming窗和Blackman窗等，称为余弦窗。其一般表达式为： （），窗系数满足。当H1，就是矩形窗。H以及系数不同的值可以构成不同性能的窗。几种窗函数的最大旁瓣幅度和主瓣宽度的比较见表1所示。表中的N为窗的宽度。其中，矩形窗的主瓣最窄，但不幸的是旁瓣衰减较慢；用平滑变化的窗如Hanning、Hamming以及Blackman可以使旁瓣电平大幅度下降，可以使旁瓣的泄漏有较大幅度的下降，将泄漏主要控制在主瓣范围内，但其代价是主瓣宽度增加了一半以上。窗函数的主瓣宽度与N是成反比的，因此，频谱泄漏程度由窗函数的类型和宽度N共同确定，一般而言，窗函数的N越大，主瓣越窄，频谱能量泄漏的范围就越小。表1：几种常用窗的性能比较窗函数类型主瓣宽度（）旁瓣峰值幅度（dB）阻带最小衰减（dB）多项式项数H矩形窗213210Hanning窗432442Hamming窗441532Blackman窗657744从几种窗函数的特性来看，Blackman窗虽然旁瓣衰减大，但其主瓣相对较大，且其计算相对复杂，不利于对空化噪声实时在线监测与分析。Hamming窗的旁瓣衰减略比Hanning窗大，但随旁瓣增加，其衰减速度很慢，对于频谱分析不利。算法的复杂度与所选窗的复杂程度有关，若选用Hanning窗，不但计算量较小，并且可通过调节采样长度达到减少频谱泄漏的目的，且可满足信号分析的实时性。因此，在空化噪声的频谱分析中，选用Hanning窗作为窗函数。从时域上看，加Hanning窗实际上就是使周期采样信号的起始点与终止点的相位为0，在一定的程度上，可以克服由于空化噪声非平稳性的特性造成的采样信号相位在始端与终端不连续的现象。对于非平稳信号而言，任何窗函数均不能彻底解决主瓣处偏离零值点而导致的频谱泄漏而形成的误差。所以，选择较为合适的窗函数可以在一定的程度上使分析结果得到改善，若需进一步提高测量精度，还需借助频谱校正方法使之改善。 223 频谱校正算法 频谱校正是和加窗插值算法的发展，常规的频谱校正算法有：多点卷积幅值校正法3，能量重心校正法、峰值搜寻法45、比值公式法6以及相位差校正法7等。在实时监测系统中，频谱校正算法的处理必须考虑到由于其算法造成的额外的计算量与实时监测系统的实时性的矛盾。所以应采取尽可能简洁的校正算法减少运算时间。考虑到实际情况，本算法采用比值公式法。对于连续信号用采样频率进行自由采样得到的离散序列，再用长度为N的窗函数序列加权截断，得到一组新序列，其对应的频谱为： （1）式中： （2）式中为频谱分析的分辨率。由于水轮机的空化振动噪声的非平稳性，幅值频率很难正好位于离散谱线频点上，也即k0一般不为整数。设峰值点左右两侧的谱线分别为第k1和k2条谱线，这两条谱线也应该是峰值点附近幅值最大和次最大的谱线。显然，在离散频谱中找到这两根谱线，从而可以确定k1和k2。令这两条谱线幅值分别是和，且 （3），由于，可以引入一个辅助参数，显然，将式（3）进行变换可得：（4），令，假如所取的窗函数N较大时，式（4）一般可以简化为：，其反函数为：。当所选的窗函数为实系数时，其幅频响应是偶对称的，函数均是奇函数。在计算时，若过分追求旁瓣衰减特性而采用多项升余弦窗如Blackman窗，有可能无法求解，必须采用多项式逼近方法求解，计算复杂度较大，这也是采用简洁Hanning窗的原因之一。为了克服单峰谱线修正算法易受到频谱泄漏和噪声干扰影响的缺点，采用双峰谱线修正算法6。具体做法是：直接对k1和k2两根谱线幅值进行加权平均，权重的选择一般有：黄金分割法和幅值比例法。本次计算选择与各自的幅值成正比的幅值比例法，计算出实际的峰值点的幅值，从而可以使频谱的谱线以及幅值得到修正。根据上述修正计算的原理，信号频谱幅值修正的通用计算公式为：，式中的()为多项式偶次项系数。当采用不同的窗函数时得到所需的频谱幅值修正的公式：(1)采用Hanning窗时，其窗函数的离散函数为： （n0，1，N），其中，得到的简单实用的谱线幅值修正公式为： （5）(2)当采用矩形窗时：窗函数（n0，1，N），其中，谱线幅值的修正公式为： （6）3.对比实验31 原始空化噪声信号水轮机空化噪声根据其产生的机理是一连串的脉冲波形，而非连续波形，由于现场背景信号的干扰，致使空化噪声信号常常被淹没在较强的背景噪声之中，因此有必要在频谱分析之前对所采集的原始信号进行滤波处理。图2所示的为在水轮机尾水检修门实测经过去噪后的振动加速度空化噪声信号。图2：原始信号历时图32 窗函数对频谱分析的影响工程中改善频谱能量泄漏的主要方法是选择合适的频谱分析方法和适合所处理信号的窗函数，窗函数选择的不同，对谱线的光滑性有较大的影响。对如图2所示的信号分别选取矩形窗和Hanning窗，窗函数的宽度取，分析数据重叠率取50的窗函数长度，即，采用改进周期图法中方法进行频谱分析，所得到的频谱图见图3所示。为了对比分析方法对频谱泄漏的影响作用，列举了直接利用得到的频谱图，见图4所示。图4的谱峰较为模糊，谱线振荡较大，存在着较为严重的频谱能量泄漏问题。在图3中，为了使两谱线不重叠，特将采用Hanning窗的谱线提高了3dB。图3和图4比较，采用改进周期图法的方法时，不管采用何种窗函数截断，均比直接采用得到的谱线要光滑的多，谱线轮廓清晰。从图3可知，虽然矩形窗的主瓣宽度比Hanning窗的主瓣窄，但由于矩形窗旁瓣的衰减却比Hanning慢，导致采用矩形窗的谱峰的谱线相对而言比较模糊，而采用Hanning窗的谱峰则相对比较光滑和清晰，在一定程度上改善了频谱泄漏。 33 修正对比选择合适的谱分析方法和合适的窗函数，可以在一定的程度上降低频谱泄漏的程度，需进一步提高频谱分析精度，可以通过对其进行频谱修正来实现。图5为采用改进周期图法的频谱分析方法、窗函数选择Hanning窗、窗函数的宽度为，分析数据重叠率取50的窗函数程度，即时，分别采用式（5）的原理进行频谱修正和不进行频谱修正的频谱图。在图5中，为了使两谱线不重叠和易于辨认，同样特将采用经过修正过的谱线提高了3dB。从图5可以明显看出，空化噪声的频谱密度的主要集中在4kHz8kHz之间，在11kHz和16kHz附近也分别存在较强的谱峰能量。谱图不经