针对电磁流量计测量水煤浆时出现较大波动、甚至回零的问题,本文采集现场电磁流量计输出的水煤浆信号;在时域和频域对信号进行分析,找出了电磁流量计不能稳定测量水煤浆流量的原因;根据水煤浆信号特征,提出了基于励磁频率高次谐波分析的煤浆流量计信号处理方法;在基于DSP的电磁流量计变送器上实时实现该算法,进行现场验证。实验结果表明,测量结果较稳定,验证了所提出的算法的有效性。

1数据采集分析

1.1现场实验

针对电磁流量计测量水煤浆时出现较大波动,甚至回零这一问题,特去某煤化工企业甲醇分公司进行现场数据采集。该公司所使用的对置式四喷嘴气化有4个喷嘴,喷嘴管道口径为125mm,管中水煤浆流量基本稳定在19m3/h(流速约为0.48m/s)。每条喷嘴煤浆线上安装了3台电磁流量计,每台电磁流量计由传感器和变送器两部分组成。选择其中1条水煤浆管线上的1台电磁流量计进行数据采集,因为该台电磁流量计测量结果波动大,甚至出现回零的现象。将课题组研制的基于DSP的电磁流量变送器的信号线和励磁线接到该电磁流量传感器的电极和励磁线圈上,组合成完整的电磁流量计,进行水煤浆数据采集。使用的电磁流量变送器是以TI公司DSP芯片TMS320F28335为,采用高频励磁方案,其硬件主要包括励磁控制系统和信号采集处理系统,具体的模块有励磁驱动模块、信号调理采集模块、信号处理控制模块、人机接口模块、通信模块及电源管理模块。信号调理采集模块中的调理电路对仪表输出的信号进行放大和滤波,截止频率是2kHz,放大倍数约为230倍。通过NI公司USB-6216型号的数据采集卡进行数据采集,把调理电路的输出端连接到数据采集卡的一个差分输入端,并设置数据采集卡工作在差分的测量模式,设置采集卡的采样频率为10kHz。采集多组水煤浆信号数据,每组数据的时间长度为5min。

1.2数据分析

现场采集了25Hz方波励磁下的水煤浆信号,发现水煤浆信号的幅值非常大,甚至接近AD的量程上限,如图1所示。水煤浆信号主要由感应电动势信号和电极噪声组成。其中,感应电动势信号是由导电液体切割磁场产生的,其幅值和相同流量下介质为水的感应电动势幅值相同,仅约为数十毫伏。这是因为电磁流量计不受被测导电介质的温度、粘度、密度以及导电率的影响,只要经过水标定后,就可以用来测量其他导电液体的流量。电极噪声是水煤浆中的固体颗粒划过电极而引起的信号跳变,也称为浆液噪声,具有强非平稳性、随机性,频域具有近似1/f的特性。水煤浆信号中的浆液噪声幅值非常大,峰峰值可达数伏,远远高于与流量相关的感应电动势信号,如图2所示。这给流量信号的提取造成了极大的困难。

电磁流量计水煤浆信号整体图

采用方波励磁的电磁流量计,其传感器输出的与流量相关的感应电动势信号的波形也类似于方波。针对与流量相关的感应电动势信号f(t)的特点,可知其是由基波和奇次谐波叠加而成的。对于一个给定单峰值为Em的矩形波信号,其傅里叶展开为:


式中:g(t)表示浆液信号的幅值,特点为随机跳变的信号,波动比较大,f表示频率。浆液噪声在低频段幅值比较大,随着频率的增加,浆液噪声的幅值在减小。那么,传感器输出的信号s(t)形式为:

在传感器输出的信号中只有与流量相关的感应电动势信号才是有用信号,被用来计算流量。而提取感应电动势信号就需要包含频率等于fe,3fe,5fe,…等频率点的信号。但是,从水煤浆信号的频谱图可以看出,浆液噪声频带较宽,在频率点fe处的幅值较大,甚至将基波淹没,如图3所示。选择一组采集的水煤浆信号,把其等分成数段,利用MATLAB计算每段数据在基波处的幅值并提取保存在一个数组中,使用绘图工具画出来,如图4所示。可见,基波幅值在1~9mV波动,波动较大,而基波幅值在感应电动势信号中所占的比重又,所以,必然导致计算出的流量波动剧烈,出现测量不稳定的问题。从图3水煤浆信号的频谱图中还可以看出,随着频率的递增,水煤浆信号中的浆液噪声逐渐衰减,使高次谐波开始凸显。由式(1)可知,高次谐波的幅值也是与流量成线性关系的,因此,可以通过提取高次谐波计算流量,有效地避开浆液噪声的干扰,得到比较稳定的测量结果。

电磁流量计水煤浆信号频谱图

为了进一步研究水煤浆信号的特点,将其与纸浆信号进行对比。通过分析课题组采集的25Hz矩形波励磁下的纸浆信号发现,在同样流速下,测量介质为纸浆时,传感器输出信号经调理放大后能明显看到与流量相关的感应电动势信号,且其浆液干扰仅为数十毫伏,要远小于水煤浆信号中的浆液干扰,如图5所示。对图5所示的纸浆信号进行局部放大,得到如图6所示的信号。可见,纸浆信号中的浆液干扰持续的时间也远小于水煤浆信号中的浆液干扰,且频率较低。

电磁流量计纸浆信号整体图

在频域中对纸浆信号观察时发现,纸浆信号的浆液噪声频带在零频率点附近,距离流量信号基波频率点较远,对基波幅值和各奇次谐波幅值基本没有影响,纸浆信号在频域中的图形如图7所示。选择一组采集的纸浆信号,把其等分成数段,利用MATLAB计算每段数据在基波处的幅值并提取保存在一个数组中,使用绘图工具画出来,如图8所示。可见,基波幅值在4.7~4.95mV变化,波动较小。因此,提取到的与流量相关的感应电动势信号幅值会比较稳定。