CN105737909B - 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 - Google Patents
基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 Download PDFInfo
- Publication number
- CN105737909B CN105737909B CN201610099983.8A CN201610099983A CN105737909B CN 105737909 B CN105737909 B CN 105737909B CN 201610099983 A CN201610099983 A CN 201610099983A CN 105737909 B CN105737909 B CN 105737909B
- Authority
- CN
- China
- Prior art keywords
- vortex
- transient impact
- interference
- impact vibration
- vibration interference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
- G01F1/20—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow
- G01F1/32—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow using swirl flowmeters
- G01F1/3209—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow using swirl flowmeters using Karman vortices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
- G01F1/20—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow
- G01F1/32—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow using swirl flowmeters
- G01F1/325—Means for detecting quantities used as proxy variables for swirl
- G01F1/3259—Means for detecting quantities used as proxy variables for swirl for detecting fluid pressure oscillations
- G01F1/3266—Means for detecting quantities used as proxy variables for swirl for detecting fluid pressure oscillations by sensing mechanical vibrations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
- G01F1/20—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow
- G01F1/32—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow using swirl flowmeters
- G01F1/325—Means for detecting quantities used as proxy variables for swirl
- G01F1/3287—Means for detecting quantities used as proxy variables for swirl circuits therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
- G01F1/20—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow
- G01F1/32—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow using swirl flowmeters
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明涉及流量检测领域,为一种以单片机(MCU)为核心、基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统。该方法包括:(1)查找包含瞬态冲击振动干扰的数据段,(2)配置卡尔曼滤波器,(3)分段进行卡尔曼滤波,(4)频域幅值谱分析。采用卡尔曼滤波器预测估计涡街流量信号,达到降低强瞬态冲击振动干扰能量和比重的目的。即使存在多个瞬态冲击振动干扰分量,且干扰能量大于涡街流量信号能量时,仍然能够消除强噪声的干扰,正确提取涡街流量频率,确保涡街流量计在复杂工况情况下的测量精度。
Description
技术领域
本发明属于流量检测技术领域,是基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统,具体涉及一种以单片机(MCU)为核心、分段卡尔曼滤波与周期图谱分析相结合的涡街流量计抗强瞬态冲击振动干扰的信号处理方法和系统。
背景技术
涡街流量计具有无机械可动部件、使用寿命长、测量范围大等优点,并且适用于多种介质如液体、气体和饱和蒸汽的测量,在流程工业中得到了广泛应用。但是,涡街流量计是一种基于流体振动工作原理的流量计,由漩涡发生体、涡街流量传感器和信号处理系统组成。涡街流量传感器对管道振动非常敏感,当振动干扰能量超过流量信号能量时,常规幅值占优的频谱分析方法无法排除该干扰,会影响到涡街流量计测量的准确性。因此,涡街流量计信号处理的难点在于如何从含有强振动干扰噪声的混合信号中提取出涡街流量信号。管道振动干扰可分为周期性正弦振动干扰和瞬态冲击振动干扰。
周期性正弦振动干扰主要是由应用现场的电机、气泵、水泵等设备的机械振动产生。每个振源分别产生单一、固定频率的振动干扰。当存在这种振动干扰时,涡街流量传感器输出多个不同频率的正弦信号叠加而成的平稳信号。针对周期性正弦振动干扰,国内外的学者进行了大量的研究,取得了较好的抗周期性正弦振动干扰的效果(L.GeraldSchlatter,Douglas William Barrett,F.John Waers,H.Lee Gilbert,and J.MarkElder.Signal processing method and apparatus for flowmeters.InternationalPatent,WO90/04230,Apr.19,1990;C.-L.Shao,K.-J.Xu,and M.Fang.Frequency-VarianceBased Antistrong Vibration Interference Method for Vortex Flow Sensor[J].IEEETrans.on Instrumentation and Measurement,2014,63(6):1566-1582;罗清林.抗强振动干扰数字涡街流量计信号处理方法研究与低功耗实现[D].合肥工业大学,2010;邢娟,张涛,郝松.管道振动对涡街流量计测量影响的试验研究[J].振动与冲击,2009,28(3):112-115)。
瞬态冲击振动干扰主要由管道敲击、水击冲击、气蚀冲击和脉动冲击等造成。管道系统在受到瞬态冲击力的作用时,涡街流量传感器输出信号的幅值突然增大,之后逐渐衰减至稳定状态。经频谱分析,瞬态冲击将产生多个干扰频率分量,甚至可能存在一个或多个干扰分量的能量超过流量信号的能量。针对涡街流量计抗瞬态冲击振动干扰,国内外的研究较少。其中,有些学者通过改进涡街流量传感器的结构来提升涡街流量计抗瞬态冲击振动干扰的能力(J.J.Miau,C.C.Hu,J.H.Chou.Response of a vortex flowmeter toimpulsive vibrations[J].Flow Measurement and Instrumentation,2000(11):41-49;潘岚,宋开臣,徐国梁.高抗干扰性能涡街流量计传感器的研究[J].中国计量学院学报,2005,16(4):268-270,278)。但是,当出现连续或较强瞬态冲击时,仍然会造成涡街流量计测量失误。针对机械冲击,有些学者提出盲源分离算法,对分离瞬态冲击振动干扰具有一定的效果(S.Hao,S.Jegelka,and A.Gretton.Fast Kernel-Based Independent ComponentAnalysis[J].IEEE Trans.on Signal Processing,2009,57(9):3498-3511;B.M.D.Vos,I.J.Taelman,and S.V.Huffel.Source separation from single-channel recordings by combining empirical-mode decomposition and independentcomponent analysis[J].IEEE Trans.on Biomedical Engineering,2010,57(9):2188-2196;.J.Antoni.Blind separation of vibration components:Principles anddemonstrations[J].Mechanical Systems and Signal Processing,2005,19(6):1166-1180)。但是,盲源分离算法运算量大,无法满足低功耗涡街流量计实时性的要求。中国发明专利公布了一种抗强瞬态冲击振动干扰的涡街流量计信号处理方法(徐科军,任保宏,邵春莉,方敏,舒张平.一种基于数据替换的涡街流量计抗低频强瞬态冲击振动的信号处理方法,申请发明专利,201510021818.6)。查找敲击振动起始点和替换数据段,然后通过用反向替换数据向左右两个方向不断镜像复制替换有敲击振动干扰的数据来消除强敲击振动干扰。
发明内容
本发明要解决的问题是:在存在强瞬态冲击振动干扰的时候,涡街流量传感器输出信号中含有多个干扰频率分量,绝大部分落在涡街流量信号的频率范围内,且干扰分量的能量大于流量信号的能量。目前,盲源分离算法无法用于低功耗涡街流量计的实时实现;带通滤波或自适应陷波不适用于该类强干扰的滤波;基于幅值占优的直接频谱分析也将产生错误的结果。因此,涡街流量计需要有效的抗强瞬态冲击振动干扰的数字信号处理方法。
本发明的技术方案是:对涡街流量传感器输出的一组数据进行分段监测,查找包含强瞬态冲击振动干扰的数据段;针对该数据段分段进行卡尔曼滤波,以降低强瞬态冲击振动干扰的能量和比重(比重为该数据段中瞬态冲击产生的最大干扰分量与涡街流量信号的幅值比),使经过滤波后涡街流量传感器输出的该组数据中流量信号的能量最大,即涡街流量信号在整体上能量占优;对滤波过后的整组数据进行频域幅值谱分析,提取幅值谱中最大峰值对应的频率即为涡街流量信号的频率。
具体的技术解决方案如下:
本发明对涡街流量传感器输出信号进行处理以降低强瞬态冲击振动干扰的能量和比重。首先,监测涡街流量传感器输出信号,分段比较后设置突变阈值,根据瞬态冲击振动干扰的幅值突变后振荡衰减的特性,查找包含瞬态冲击振动干扰的数据段;然后,配置用于预测和估计涡街流量信号的卡尔曼滤波器的变量和参数;并用此类卡尔曼滤波器对包含瞬态冲击振动干扰的数据段分段进行滤波,降低该数据段中强瞬态冲击振动干扰的能量和比重,使涡街流量信号在整体上能量占优;最后,通过频域幅值谱分析,提取涡街流量信号频率。基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法的具体处理步骤如下:
(1)查找包含瞬态冲击振动干扰的数据段
首先,采集涡街流量传感器输出信号,将2060点采样数据平均分段,根据每段数据的最大峰峰值计算突变阈值Vth;将信号的幅值与阈值进行比较,查找包含瞬态冲击振动干扰的起始数据段,并计数每次冲击数据段持续的段数N;查找瞬态冲击的起始点,并保存包含瞬态冲击振动干扰的数据段。
(2)配置卡尔曼滤波器
采用卡尔曼滤波器来预测和估计涡街流量信号。将涡街流量信号定义为状态变量Xk,将包含瞬态冲击的数据段定义为观测变量Yk,将瞬态冲击振动干扰信号定义为测量噪声vk。
该过程无控制参数,所以,控制输入变量Uk=0,系数B=0。
正常工况下涡街流量信号属于平稳正弦信号,可用上一时刻的涡街流量信号加上干扰噪声来预测估计当前状态的涡街流量信号,并用过程噪声来表示干扰噪声,则设置系数A=1,系数G=1,过程噪声wk的均值并可由涡街流量信号段数据的方差近似地计算过程噪声的方差Q。
待滤波数据段是由瞬态冲击振动干扰叠加涡街流量信号所组成,则设置系数H=1;因为强瞬态冲击振动干扰的能量一般大于流量信号的能量,为了降低振动干扰的比重,根据大量实验数据的验证,设置D的计算公式为:
式中,max(PPV)为20个峰峰值中的最大值,min(PPV)为20个峰峰值中的最小值;同时,测量噪声vk的均值并用包含瞬态冲击振动干扰数据段的方差近似地计算测量噪声方差R,其计算公式为:
式中,S2为方差,T为每段卡尔曼滤波包含的待滤波数据的点数,Yi为当前的观测变量值,为T个观测变量的均值。
(3)分段进行卡尔曼滤波
将每个包含瞬态冲击振动干扰的数据段均分成两段,分别计算每段测量噪声的方差,并进行卡尔曼滤波。包含瞬态冲击振动干扰数据段之外的数据保持不变。
(4)频域幅值谱分析
分段滤波有效地降低了强瞬态冲击振动干扰的能量和比重,使涡街流量信号在整体上能量占优。针对滤波后的采样数据,可直接通过频域幅值谱分析,查找幅值谱中最大峰值对应的频率,把它作为涡街流量信号的频率。
本发明的优点是:采用基于分段卡尔曼滤波的数字信号处理方法,将降低强瞬态冲击振动干扰能量和比重的目的巧妙地转换成预测和估计涡街流量信号,即使存在多个瞬态冲击振动干扰分量,且干扰分量的能量大于涡街流量信号的能量时,仍然能够消除强噪声干扰,正确提取涡街流量信号的频率,确保涡街流量计在复杂工况下的测量精度。
附图说明
图1是涡街流量计中信号处理系统硬件结构框图。
图2是涡街流量计中信号处理系统软件结构框图。
图3是涡街流量计中信号处理系统主监控程序流程图。
图4是涡街流量传感器输出信号的时域波形图。
图5是涡街流量传感器输出信号的频域幅值谱。
图6是基于卡尔曼滤波的抗瞬态冲击干扰算法的原理框图。
图7是查找包含瞬态冲击振动干扰数据段的示意图。
图8是分段进行卡尔曼滤波示意图。
图9是经过分段卡尔曼滤波后涡街流量传感器输出信号的时域波形图。
图10是经过分段卡尔曼滤波后涡街流量传感器输出信号的频域幅值谱。
具体实施方式
下面结合附图对本发明作进一步的说明。
涡街流量计中抗强瞬态冲击振动干扰的信号处理系统硬件结构框图如图1所示。本发明的系统硬件主要包括前向输入调理电路模块、数字信号处理模块、系统输出模块和电源转换电路模块。前向输入调理电路由电荷放大器、电压放大器、低通滤波器和电压跟随器组成。数字信号处理模块由MSP430F5418单片机、欠压监测电路、复位电路、铁电存储电路和外部看门狗电路组成。输出模块由通信电路和LCD显示组成。电源转换电路将24V的直流电源转换成3V的电压,供数字电路和调理电路使用。
涡街流量传感器选用压电式传感器,本发明的基本工作过程为:压电式传感器将涡街流量信号转换成电荷信号输出,输出信号先后经过电荷放大器、电压放大器、低通滤波器和电压跟随器的放大和滤波,然后进入MSP430F5418单片机的内部的ADC单元,经过ADC采样并转换成数字信号;MSP430F5418单片机采用基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法对该数字信号进行处理,得到涡街流量信号的频率,进而得到流体的流量;然后,通过输出模块将涡街信号的频率和流量值在LCD上实时显示。
涡街流量计中信号处理系统软件结构框图如图2所示。本发明的系统软件采用模块化设计,通过主监控程序来进行调控。主要包括主监控程序、初始化模块、中断模块、计算模块、看门狗模块、铁电存储模块和LCD显示模块。
涡街流量计中信号处理系统主监控程序流程图如图3所示。本发明的涡街流量计信号处理系统软件的工作过程为:系统上电,首先进行初始化,配置各个需要的模块;提取2060点ADC采样的最新数据;分段查找这组数据中是否存在强瞬态冲击振动干扰;若存在干扰,对包含瞬态冲击振动干扰的数据段分段进行卡尔曼滤波,降低瞬态冲击振动干扰的能量和比重,若不存在干扰,则忽略该步骤;计算处理过后的最新的2048点数据的幅值谱;提取当前涡街流量信号峰值对应的频率,根据涡街流量频率值计算瞬时流量;然后,刷新LCD,显示当前涡街流量的频率,完成一个循环;主监控程序返回,进入新一轮的循环。
本发明以图4所示的一组数据来说明抗强瞬态冲击振动干扰的具体方法。给定气体流量约57Hz时,ADC采样的涡街流量传感器输出的最新2060点数据的时域波形图如图4所示。此时,数据中存在强瞬态冲击振动干扰,其2048点频域幅值谱如图5所示。图5中最大峰值点(93.99Hz,0.1657V)是瞬态冲击振动的一个干扰分量的峰值点,其幅值远远大于气体流量信号对应的峰值点(57.37Hz,0.06677V)的幅值。此时,若采用基于幅值占优的幅值谱信号处理方法将提取93.99Hz的干扰分量作为涡街流量频率输出,不能得到正确的流量信息。
本发明的抗强瞬态冲击振动干扰算法原理框图如图6所示。首先,监测涡街流量传感器输出的信号,分段比较后设置突变阈值;根据瞬态冲击振动干扰的幅值突变后振荡衰减的特性,查找包含瞬态冲击振动干扰的数据段;然后,配置用于预测和估计涡街流量信号的卡尔曼滤波器的变量和参数;并用此类卡尔曼滤波器对包含瞬态冲击振动干扰的数据段分段进行滤波,降低该数据段中强瞬态冲击振动干扰的能量和比重,使涡街流量信号在整体上能量占优;最后,通过频域幅值谱分析,提取涡街流量信号的频率。下面,具体介绍基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法的处理流程:
(1)查找包含瞬态冲击振动干扰的数据段
突变阈值的设置如图4中双向箭头所示。首先将2060点采样数据均分成20段,每段包含103点数据。查找每段数据的最大值和最小值,然后计算两者差值(峰峰值)。设置20个峰峰值中最小的三个数值的均值为突变阈值Vth。从第一段开始,查找并保存最值(最值包括最大值和最小值)绝对值大于Vth的数据段。
通过大量的实验数据分析,存在瞬态冲击振动干扰的起始段数据主要有以下三大类特征。第一类:该段只有最大值的绝对值大于Vth,此最大值之后不存在相邻波谷,且下一段最大值和最小值的绝对值均超过突变阈值;第二类:该段只有最小值的绝对值大于Vth,此最小值之后不存在相邻波峰,且下一段最大值和最小值的绝对值均超过突变阈值;第三类:该段最大值和最小值的绝对值均超过突变阈值。
查找包含瞬态冲击振动干扰数据段的示意图如图7所示。根据上述三大类特征查找存在瞬态冲击振动干扰的起始数据段,并计数每次冲击数据持续的段数N。图7中所示的第8段即为瞬态冲击的起始数据段,且与其相邻的第9段、第10段和第11段都存在最值绝对值大于Vth的情况,则此4段数据属于1次瞬态冲击,即该次瞬态冲击的N=4。上述查找包含瞬态冲击振动干扰的数据段的方法基本上能够排除毛刺干扰,但还要区分是否是流量变大。综合考虑系统的响应速度和涡街流量信号提取的准确性,只有下面几个条件同时满足时才认为是流量变大。这几个条件分别是:2060点采样数据中只有1次幅值突变的情况;最后一段数据最值的绝对值大于Vth;N≥10,即至少有一半的数据属于流量变大后的数据。
确定存在强瞬态冲击振动干扰后,查找包含瞬态冲击振动干扰的起始段中第一个超过突变阈值的极值点作为第一个突变点,即图7中横坐标为792的点;以第一突变点左侧第一个过零点作为瞬态冲击起始点,即图7中横坐标785所对应的点。保存该起始点以及起始点之后的103*N点数据。该段数据即为包含瞬态冲击振动干扰的数据段,如图7中虚线框所包含的数据段[785,1197]。
(2)配置卡尔曼滤波器
卡尔曼滤波器的数学公式用状态空间的概念来描述,包含描述状态变量的状态方程(式(1))和描述观测变量的观测方程(式(2))。
Xk+1=AXk+BUk+Gwk (1)
Yk+1=HXk+1+Dvk+1 (2)
式中,Xk为k时刻的状态向量,Xk+1为k+1时刻的状态向量,Uk为k时刻的控制输入向量,Yk+1为k+1时刻的观测向量,A、B、D、G和H为已知的系数矩阵,详细的配置方法在下面通过结合具体的应用来介绍,wk和vk+1分别为过程噪声和测量噪声。同时,wk和vk+1是相互独立的高斯白噪声,其噪声均值分别为和协方差分别为Q和R。
本发明将降低强瞬态冲击振动干扰能量和比重的目的巧妙地转换成预测估计涡街流量信号。下面设置离散卡尔曼滤波器的变量和参数,以图4所示的涡街流量传感器输出信号为例来预测估计涡街流量信号。首先,将涡街流量信号定义为状态变量Xk,包含瞬态冲击振动干扰的数据段定义为观测变量Yk。瞬态冲击振动干扰信号定义为测量噪声vk。该过程无控制参数,所以,控制输入变量Uk=0,系数B=0。
正常工况下的涡街流量信号属于平稳正弦信号,可用上一时刻的涡街流量信号加上干扰噪声来预测估计当前状态的涡街流量信号,并用过程噪声来表示干扰噪声,则设置系数A=1,系数G=1,过程噪声wk的均值并可由涡街流量信号段数据的方差近似地计算过程噪声的方差Q,图4所示波形的Q约为0.05。
待滤波数据段是由瞬态冲击振动干扰叠加涡街流量信号所组成,则设置系数H=1;因为强瞬态冲击振动干扰的能量一般大于涡街流量信号的能量,为了降低振动干扰的比重,根据大量实验数据的验证,本发明按照式(3)设置D。无流量或接近流量下限时,涡街流量计受强瞬态冲击振动干扰的影响最大。为了保证卡尔曼滤波算法的有效性,可适当增加D的数值,图4所示波形的D=5.48。
式中,max(PPV)为20个峰峰值中的最大值,min(PPV)为20个峰峰值中的最小值;同时,测量噪声vk的均值并用包含瞬态冲击振动干扰数据段的方差近似地计算测量噪声方差R,其计算公式为
式中,S2为方差,T为每段卡尔曼滤波包含的待滤波数据的点数,Yi为当前的观测变量值,为T个观测变量的均值。
起始点的采样值作为初始状态变量X0|0;设置初始预测误差协方差P0|0=1。
(3)分段进行卡尔曼滤波
分段进行卡尔曼滤波示意图如图8所示。图4所示波形中包含3次瞬态冲击振动干扰,包含瞬态冲击振动干扰的数据段分别为[172,481]、[785,1197]和[1424,1733],如图8中双向箭头所示的范围。这3个数据段不仅包含瞬态冲击振动干扰,还包含涡街流量信号。因此,在对该类数据段进行卡尔曼滤波的同时,该段包含的涡街流量信号的能量也会有所衰减。然而,每个包含瞬态冲击振动干扰的数据段都包含幅值突变段和振荡衰减段。为了降低瞬态冲击振动干扰的能量,且尽可能地减小对涡街流量信号的衰减,通过大量实验数据的测试,本发明将每个包含瞬态冲击振动干扰的数据段均分成两段,分别进行卡尔曼滤波,分界线如图8中的红色短虚线所示。前半段近似为幅值突变段,后半段近似为振荡衰减段。这样,3段包含瞬态冲击振动干扰的待滤波数据段就分为6段。仅对这6段数据段进行卡尔曼滤波,而这6段之外的数据保持不变。针对这6段数据,分别计算其测量噪声方差R,并依次分别进行卡尔曼滤波。
经过分段进行卡尔曼滤波后,得到滤波后涡街流量传感器输出信号的时域波形图如图9所示。对比图4和图9滤波前后的时域波形图可以看出,分段卡尔曼滤波过程有效地降低了瞬态冲击振动干扰的能量。
(4)频域幅值谱分析
对图9所示的滤波后的涡街流量传感器输出数据进行频谱分析。取2060点数据中的后2048点数据,得到经过分段卡尔曼滤波后涡街流量传感器输出信号的频域幅值谱如图10所示。该幅值谱中气体流量信号的峰值点为(57.37Hz,0.05256V),瞬态冲击振动的最大干扰分量对应的峰值点为(93.99Hz,0.03219V)。通过图5和图10滤波前后涡街流量传感器输出信号的频域幅值谱对比可以看出,滤波前,由于受到强瞬态冲击振动的干扰,有多个超过流量信号幅值的干扰峰值;滤波后,干扰分量的能量大幅降低,涡街流量信号在整体上能量占优。此时,提取幅值谱中的最大峰值即为涡街流量信号对应的峰值。
正常工况情况下,在没有强瞬态冲击振动干扰的时候,涡街流量信号的能量占优,所以,直接查找幅值谱中的最大峰值即为涡街流量信号峰值;在有强瞬态冲击振动干扰的时候,通过图6所示的信号处理方法,先查找含有瞬态冲击振动干扰的数据段。然后配置卡尔曼滤波器变量和参数,分别对这些数据段分段进行卡尔曼滤波,降低强瞬态冲击振动干扰的能量和比重后,一般也能达到涡街流量信号在整体上能量占优的效果。此时,再提取幅值谱中的最大峰值即为涡街流量信号对应的峰值。
Claims (2)
1.基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法,包括前向输入调理电路模块、数字信号处理模块、系统输出模块、电源转换电路模块和基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法软件;压电式传感器将涡街流量信号转换成电荷信号输出,输出信号先后经过电荷放大器、电压放大器、低通滤波器和电压跟随器,然后进入MSP430F5418单片机的内部的ADC单元,经过ADC采样并转换成数字信号;MSP430F5418单片机采用基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法对该数字信号进行处理,得到涡街流量信号的频率,进而得到流体的流量;然后,通过输出模块将涡街信号的频率和流量值在LCD上实时显示;其特征在于:基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法包括查找包含瞬态冲击振动干扰的数据段,配置卡尔曼滤波器,分段进行卡尔曼滤波,频域幅值谱分析;采用卡尔曼滤波器预测估计涡街流量信号,达到降低强瞬态冲击振动干扰能量和比重的目的。
2.如权利要求1所述的基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法,其特征在于:对涡街流量传感器输出信号进行处理以降低强瞬态冲击振动干扰的能量和比重;首先,监测涡街流量传感器输出信号,分段比较后设置突变阈值,根据瞬态冲击振动干扰的幅值突变后振荡衰减的特性,查找包含瞬态冲击振动干扰的数据段;然后,配置用于预测和估计涡街流量信号的卡尔曼滤波器的变量和参数;并用此类卡尔曼滤波器对包含瞬态冲击振动干扰的数据段分段进行滤波,降低该数据段中强瞬态冲击振动干扰的能量和比重,使涡街流量信号在整体上能量占优;最后,通过频域幅值谱分析,提取涡街流量信号频率;基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法的具体处理步骤如下:
(1)查找包含瞬态冲击振动干扰的数据段
首先,采集涡街流量传感器输出信号,将2060点采样数据平均分段,根据每段数据的最大峰峰值计算突变阈值Vth;将信号的幅值与阈值进行比较,查找包含瞬态冲击振动干扰的起始数据段,并计数每次冲击数据段持续的段数N;查找瞬态冲击的起始点,并保存包含瞬态冲击振动干扰的数据段;
(2)配置卡尔曼滤波器
采用卡尔曼滤波器来预测和估计涡街流量信号;将涡街流量信号定义为状态变量Xk,将包含瞬态冲击的数据段定义为观测变量Yk,将瞬态冲击振动干扰信号定义为测量噪声vk;
该过程无控制参数,所以,控制输入变量Uk=0,系数B=0;
正常工况下涡街流量信号属于平稳正弦信号,可用上一时刻的涡街流量信号加上干扰噪声来预测估计当前状态的涡街流量信号,并用过程噪声来表示干扰噪声,则设置系数A=1,系数G=1,过程噪声wk的均值并可由涡街流量信号段数据的方差近似地计算过程噪声的方差Q;
待滤波数据段是由瞬态冲击振动干扰叠加涡街流量信号所组成,则设置系数H=1;因为强瞬态冲击振动干扰的能量一般大于流量信号的能量,为了降低振动干扰的比重,根据大量实验数据的验证,设置D的计算公式为:
式中,max(PPV)为20个峰峰值中的最大值,min(PPV)为20个峰峰值中的最小值;同时,测量噪声vk的均值并用包含瞬态冲击振动干扰数据段的方差近似地计算测量噪声方差R,其计算公式为:
式中,S2为方差,T为每段卡尔曼滤波包含的待滤波数据的点数,Yi为当前的观测变量值,为T个观测变量的均值;
(3)分段进行卡尔曼滤波
将每个包含瞬态冲击振动干扰的数据段均分成两段,分别计算每段测量噪声的方差,并进行卡尔曼滤波;包含瞬态冲击振动干扰数据段之外的数据保持不变;
(4)频域幅值谱分析
分段滤波有效地降低了强瞬态冲击振动干扰的能量和比重,使涡街流量信号在整体上能量占优;针对滤波后的采样数据,可直接通过频域幅值谱分析,查找幅值谱中最大峰值对应的频率,把它作为涡街流量信号的频率。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610099983.8A CN105737909B (zh) | 2016-02-23 | 2016-02-23 | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 |
PCT/CN2016/079350 WO2017143649A1 (zh) | 2016-02-23 | 2016-04-15 | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 |
US15/955,720 US10876869B2 (en) | 2016-02-23 | 2018-04-18 | Kalman filter based anti-transient-impact-vibration-interference signal processing method and system for vortex flowmeter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610099983.8A CN105737909B (zh) | 2016-02-23 | 2016-02-23 | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105737909A CN105737909A (zh) | 2016-07-06 |
CN105737909B true CN105737909B (zh) | 2018-07-13 |
Family
ID=56248414
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610099983.8A Active CN105737909B (zh) | 2016-02-23 | 2016-02-23 | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10876869B2 (zh) |
CN (1) | CN105737909B (zh) |
WO (1) | WO2017143649A1 (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105737909B (zh) | 2016-02-23 | 2018-07-13 | 合肥工业大学 | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 |
CN106123973B (zh) * | 2016-09-08 | 2019-02-01 | 合肥工业大学 | 基于互相关分析测量低雷诺数流量的涡街流量计 |
CN106871980B (zh) * | 2017-02-13 | 2019-03-26 | 合肥工业大学 | 基于相邻峰值最大差值的气体超声流量计间歇式激励和信号处理方法和系统 |
CN107063370B (zh) * | 2017-04-18 | 2019-05-07 | 中国工程物理研究院材料研究所 | 一种热分布式气体质量流量计的时间响应补偿方法 |
CN109974793B (zh) * | 2019-04-22 | 2020-08-04 | 合肥工业大学 | 一种电磁式涡街流量计测量含气导电液体流量的信号处理方法 |
FR3098297B1 (fr) * | 2019-07-03 | 2021-06-04 | Etat Francais Represente Par Le Delegue General De L’Armement | Dispositif de detection de l'approche d'un etat d'anneaux tourbillonaires, aerodyne a voilure tournante le comprenant, et procede associe |
CN110532509B (zh) * | 2019-09-29 | 2023-03-21 | 中国计量大学 | 一种基于不确定度分析的泵与风机性能预测方法 |
CN112747800B (zh) * | 2020-12-30 | 2024-04-05 | 金卡智能集团股份有限公司 | 一种用于计量仪表的流量滤波方法和计量仪表 |
CN113029258A (zh) * | 2021-02-05 | 2021-06-25 | 辽宁聚焦科技有限公司 | 基于涡街幅频特性抗干扰的涡街流量计数字信号处理系统 |
CN113959508A (zh) * | 2021-10-14 | 2022-01-21 | 安徽大学 | 基于目标函数最优估计的涡街传感器输出瞬态冲击干扰的识别方法 |
CN114166291B (zh) * | 2021-11-30 | 2024-09-10 | 合肥科迈捷智能传感技术有限公司 | 一种小口径涡街流量计抗振动干扰方法 |
CN114637960B (zh) * | 2022-03-31 | 2024-02-23 | 合肥工业大学 | 一种时频域数字涡轮流量传感器信号处理方法 |
CN115235554A (zh) * | 2022-07-12 | 2022-10-25 | 国网安徽省电力有限公司电力科学研究院 | 一种基于Kalman滤波的绝缘油取样流量检测方法及系统 |
CN115265723B (zh) * | 2022-09-28 | 2022-12-13 | 南通有来信息技术有限公司 | 基于数据处理的涡街流量计异常监测方法及系统 |
CN115307863B (zh) * | 2022-10-12 | 2022-12-09 | 中国空气动力研究与发展中心低速空气动力研究所 | 发动机进气模拟的稳流量进气控制方法、系统及存储介质 |
CN116595324B (zh) * | 2023-07-19 | 2023-11-10 | 北谷电子股份有限公司 | 信号瞬态冲击始点的提取方法 |
CN117909659B (zh) * | 2024-03-20 | 2024-05-14 | 中国空气动力研究与发展中心超高速空气动力研究所 | 变体飞行器风洞动态试验天平气动力信号处理方法及装置 |
CN118349791B (zh) * | 2024-06-18 | 2024-09-06 | 大连派立特电子科技发展有限公司 | 基于电参数分析的电源运行状态稳定性监测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1288149A (zh) * | 1999-09-14 | 2001-03-21 | 横河电机株式会社 | 旋涡流量计 |
CN101881638A (zh) * | 2009-05-12 | 2010-11-10 | 石家庄西比克仪表有限公司 | 基于mcu的涡街流量计信号数字处理系统 |
CN104729591A (zh) * | 2015-01-16 | 2015-06-24 | 合肥工业大学 | 一种基于数据替换的涡街流量计抗低频强瞬态冲击振动的信号处理方法 |
CN105181031A (zh) * | 2015-05-20 | 2015-12-23 | 杭州子午仪器有限公司 | 一种基于多传感器信息融合的涡街流量计 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU4519589A (en) | 1988-10-14 | 1990-05-01 | Engineering Measurements Company | Signal processing method and apparatus for flowmeters |
JP3216766B2 (ja) * | 1995-02-21 | 2001-10-09 | 富士電機株式会社 | カルマン渦流量計 |
JP2003149023A (ja) * | 2001-11-08 | 2003-05-21 | Yokogawa Electric Corp | 渦流量計 |
JP4405167B2 (ja) * | 2003-03-24 | 2010-01-27 | トキコテクノ株式会社 | 超音波式渦流量計 |
CN101476906B (zh) * | 2009-01-21 | 2010-08-18 | 合肥工业大学 | 抗强固定干扰的涡街流量计数字信号处理系统 |
CN105737909B (zh) | 2016-02-23 | 2018-07-13 | 合肥工业大学 | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 |
-
2016
- 2016-02-23 CN CN201610099983.8A patent/CN105737909B/zh active Active
- 2016-04-15 WO PCT/CN2016/079350 patent/WO2017143649A1/zh active Application Filing
-
2018
- 2018-04-18 US US15/955,720 patent/US10876869B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1288149A (zh) * | 1999-09-14 | 2001-03-21 | 横河电机株式会社 | 旋涡流量计 |
CN101881638A (zh) * | 2009-05-12 | 2010-11-10 | 石家庄西比克仪表有限公司 | 基于mcu的涡街流量计信号数字处理系统 |
CN104729591A (zh) * | 2015-01-16 | 2015-06-24 | 合肥工业大学 | 一种基于数据替换的涡街流量计抗低频强瞬态冲击振动的信号处理方法 |
CN105181031A (zh) * | 2015-05-20 | 2015-12-23 | 杭州子午仪器有限公司 | 一种基于多传感器信息融合的涡街流量计 |
Also Published As
Publication number | Publication date |
---|---|
CN105737909A (zh) | 2016-07-06 |
US20180231406A1 (en) | 2018-08-16 |
US10876869B2 (en) | 2020-12-29 |
WO2017143649A1 (zh) | 2017-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105737909B (zh) | 基于卡尔曼滤波的涡街流量计抗瞬态冲击干扰的信号处理方法和系统 | |
CN104729591B (zh) | 一种基于数据替换的涡街流量计抗低频强瞬态冲击振动的信号处理方法 | |
Lin et al. | Sequential spindle current-based tool condition monitoring with support vector classifier for milling process | |
Raj et al. | Early classification of bearing faults using morphological operators and fuzzy inference | |
Green et al. | Detection of Lagrangian coherent structures in three-dimensional turbulence | |
CN101701834B (zh) | 基于双传感器的抗强干扰的数字涡街流量计 | |
Frederich et al. | Modal analysis of complex turbulent flow | |
CN104217112B (zh) | 一种基于多类型信号的电力系统低频振荡分析方法 | |
CN105929022A (zh) | 钢丝绳脉冲电涡流无损探伤检测装置及检测方法 | |
CN106202977A (zh) | 一种基于盲源分离算法的低频振荡模式分析方法 | |
Xiao et al. | A multichannel fusion approach based on coupled hidden Markov models for rolling element bearing fault diagnosis | |
Zhu | Edge detection based on multi-structure elements morphology and image fusion | |
Deng et al. | Adaptive bandwidth Fourier decomposition method for multi-component signal processing | |
CN103493321A (zh) | 用于产生故障信号的方法 | |
Lv et al. | High-order synchroextracting chirplet transform for accurate instantaneous frequency estimation and its application in fault diagnosis of rotary machinery | |
CN109580224A (zh) | 滚动轴承故障实时监测方法 | |
CN104865508A (zh) | 基于数据分组量化的局部放电识别方法 | |
Saadawia et al. | Improved process monitoring and supervision based on a reliable multi-stage feature-based pattern recognition technique | |
CN107561359A (zh) | 一种基于布莱克曼窗的间谐波提取方法 | |
CN106596025A (zh) | 基于冲激响应的公路隧道悬挂风机基础稳定性检测方法及系统 | |
CN113537156B (zh) | 一种基于区间标准差结合频谱分析的振动数据清洗方法 | |
Shao et al. | Segmented Kalman filter based antistrong transient impact method for vortex flowmeter | |
CN112857806B (zh) | 一种基于移动窗口时域特征提取的轴承故障检测方法 | |
CN106124371A (zh) | 一种基于静电法的气固两相流细度测量装置及测量方法 | |
Zhou et al. | Multi-cavitation states identification of a sewage pump using CEEMDAN and BOA-SVM |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |