基于能量变化率的气体超声波流量计信号处理方法
- 格式:pdf
- 大小:423.55 KB
- 文档页数:7
第36卷第9期 2015年9月
仪 器 仪 表 学 报
Chinese Journal of Scientific Instrument V01.36 No.9
Sep.2015
基于能量变化率的气体超声波流量计 信号处理方法
沈子文 ,徐科军 ,方敏 ,朱文姣 ,汪伟 (1.合肥工业大学电气与自动化工程学院合肥230009; 2.工业自动化安徽省工程技术研究中心合肥230009)
摘要:针对气体超声流量计信号处理中不易找到稳定的测量特征点的问题,提出了一种基于能量变化率的过零检测信 号处理方法,以及最优阈值参数的选取方法。基于能量变化率的过零检测信号处理方法首先求取能量变化率曲线上的 关键阈值点,然后回溯到原始波形信号上,找到对应的过零点,并以此作为信号的稳定特征点,进而求得信号的传播时间 及气体流量。在以FPGA和DSP为核心的气体超声流量计信号处理系统上实时实现了算法,并进行标定实验,实验结果 验证了方法的有效性。 关键词:超声波;传播时间;能量变化率;过零检测;数字信号处理 中图分类号:TH814.92 文献标识码:A 国家标准学科分类代码:470.40
Rate of energy based signal processing method for ultrasonic gas flowmeter
Shen Ziwen ,Xu Kejun 一,Fang Min ,Zhu Wenjiao ,Wang Wei ( .School of Electrical and Automation Engineering,Hefei University of Technology,Hefei 230009,China; 2.Engineering Technology Research Center of Industrial Automation,Anhui Province,Hefei 230009,China)
Abstract:Aiming at the problem that it is difficult to find stable feature point in the signal processing of ultrasonic gas flowmeter, a zero—crossing detection signal processing method based on the rate of energy and a method for determining the optimum threshold parameter are proposed.Firstly,this zero—crossing detection signal processing method obtains the key threshold point on the rate of energy cuive,then traces back to the original waveform signal and obtains the correlated zero—crossing point,which is taken as the stable feature point of the signa1.Finally,the propagation time of the ultrasonic signal is calculated and the gas flow rate is meas— ured accurately.This method was implemented in real time in the signal processing system of the ultrasonic gas flowmeter with DSP and FPGA as the dual—cores,and calibration experiment was conducted.Gas flow calibration experiment resuhs verify the effective— ness of the proposed method. Keywords:ultrasonic;propagation time;rate of energy;zero-crossing detection;digital signal processing
1 引 言 气体超声流量计以其测量精度较高、量程比大、 无压损等特点,在天然气和大口径气体流量测量领域
收稿日期:2015-06 Received Date:2015-06
具有独特的优势和广阔的发展前景 。气体超声波 流量计基于传播时间差的测量原理,需要在超声回波 信号上找到一个稳定的特征点来确定传播时间。但 是,超声信号在气体中幅值衰减严重,导致信号强度 较弱,信噪比较低。而且,由于气体回波信号是类似 幅值调制的正弦波,在有噪声干扰的情况下,常见处 第9期 沈子文等:基于能量变化率的气体超声波流量计信号处理方法 2139 理过程中采用的特征点(如最大最小值点、极大极小 值点、斜率最大值点等)都容易波动,具有不确定性, 影响流量的准确计算。 国内外对气体超声流量计均有相关研究。国外 起步较早,大多数采用互相关和过零检测的方法,计 算超声信号从发射到接收的传播时间 。Brassier P等采用互相关算法计算传播时间和时间差。但 是,由于互相关算法的计算量很大,对系统的处理能 力要求较高,实时性难以保证。美国的Daniel公司 申报的关于计算信号传播时间的方法与系统的专利 介绍了一种相同的信号在时间轴上的两个不同位置 的时间差的测量方法 。具体来说,该方法首先采 用多种可选的滤波方法去除噪声干扰,采用信号叠 加等方法来增强原始信号,然后计算信号的能量,再 利用多种可选的窗函数来对信号的能量作平均处 理,接着求得能量在时间轴上的变化率,在变化率的 曲线上,根据一个确定的百分比找到临界点,回溯到 原始的信号波形上进行最终的时间计算。但是,专 利没有披露详细的技术细节,例如,没有给出处理过 程中的多个关键参数。根据专利的思路,利用MAT. LAB对该方法进行了研究。针对专利中较多参数难 以确定的问题,利用穷举的方法,选择了一组相对最 优的计算参数。在对实际信号进行处理过程中,发 现处理结果会随机地出现错误,处理效果并不理想, 并不能满足高精度的要求。 国内的气体超声流量计研究起步较晚,大多采用模 拟过零检测的方法计算信号的传播时间,求得气体流 量 。但是,与数字信号处理方法相比,模拟信号处理方 法在稳定性、抗干扰性以及信号处理能力等方面存在较 大的劣势。国内部分高校对数字信号处理方法也有研 究,例如互相关算法 、小波分析方法 ” 。但是,仅仅 停留在仿真阶段,没有实时实现。 针对气体超声流量计中数字信号处理方法中难 以找到稳定特征点的问题,本文提出了处理气体超声 流量信号这类幅值调制的周期信号的处理方法。首 先在原始信号中提取出较长一段单调递增或者递减 的信号特征波形;然后,通过阈值选择特征波形上的 某一特征点,再回溯到原始的幅值调制的周期信号 上,找到相应的过零点。同时,提出与之相对应的最 优阈值选取方法,以获得最大容错区间,提高阈值的 可靠性。并且在以DSP和FPGA为双核心的硬件平 台上实时实现该算法,进行气体流量标定实验,取得 了较好的实验结果。 2基于能量变化率的过零检测算法及最优 阈值参数的选取 气体超声流量计由流量计表体、电子组件及微处 理器系统、超声换能器等组成。超声换能器通常沿管 壁安装,由一个超声换能器发射的超声脉冲被另一个超 声换能器所接收,反之亦然。当管道中的气体流速不为 零时,沿气流方向顺流传播的信号将加快速度,而逆流传 播的信号将减慢,据此得到气体流速与传播时间的关 系为:
= L :to-T ̄
式中: 是气体沿着声道的平均流速,,J是声道长度, 0是声道与管道轴线的夹角, 和 分别是超声脉冲 逆流、顺流传播时间 ”J。可见,逆流时间 和顺流时 间 的准确计算是气体超声流量计精确计量的关键。 发射换能器开始激励到接收换能器接收到超声回波信号 的时间为 ,示意图如图1所示,即为式(1)中的 和
。
0 O
一0 —0 0 l60 240 时问,
图1 超声信号激励及接收波形图 Fig.1 The excitation and echo waveforms of ultrasonic signal
超声回波信号的起始位置没有明确的特征,而且在 有噪声干扰情况下难以准确定位。如图1所示,如果能 够在超声回波信号上找到一个稳定的特征点,通过特征 点计算得到时间 。只要在 和 之间存在一个恒定 不变的时间差 ,就能够通过式(2)来得到准确的 。 所以,在接收到的超声信号上,找到一个稳定的特征点, 是气体超声流量信号处理的关键。
=T— (2)
;喜 2140 仪器仪表学报 第3 6卷 研究超声回波信号的波形发现,超声回波信号是类 似幅值经过调制的正弦波,其幅值逐渐增大,然后再逐渐 减小,直至衰减为零。信号的能量也会经历一个从无到 有,再从有到无的变化过程。虽然随着流量的不同,传播 时间 会发生变化,但是,超声回波信号的波形相对比 较稳定,其整体的变化趋势不发生改变,即 保持恒 定值。 为了寻找回波信号中的一个稳定特征点,就需要准 确地找到信号中一个固定的周期区间。由于超声信号与 正弦信号有关,所以,以正弦信号为例,如图2所示。信 号由于干扰的影响,会在不同的采样周期中得到图2中 标记的波动1、波动2、波动3信号。再根据信号的特征, 寻找一个在较长的区间内单调的特征曲线,而且特性曲 线的波形只与信号特性有关,与干扰无关。假设图2中 的斜线段为满足要求的特征曲线,虽然信号由于噪声干 扰会有波动,但是与干扰无关的特性曲线并不会变化。 只要信号的波动不会使图中正弦信号的两个负斜率过零 点确定的波动冗余变为零,那么就可以在纵坐标方向得 到一个阈值的选取区问。根据阈值选取区间确定的阈 值,就能够在每次采样周期中都准确地找到信号中一个 固定位置的正弦周期区间。而且,波动冗余和斜线段的 斜率越大,获得阈值的选取区间就越大,就越有利于提高 计算方法的抗干扰能力。 时间 图2方法示意图 Fig.2 Schematic diagram of the method 气体超声波信号是类似于正弦调制的信号,信号的 特点并不会因为气体介质和管道直径的变化而改变。按 照上述方法,只要通过对信号的处理,得到一段与信号特 性相关的稳定的单调曲线,就能够通过上述方法确定一 个固定的周期区间,然后通过求取过零点的方法,确定最 终的传播时间,就能够使该方法确定的传播时间误差在 可接受范围内。 基于超声信号幅值变化的特点,通过能量变化率 的方法,找到了一条满足上述要求的特征曲线。首先 对超声回波信号的幅值进行平方处理,得到超声回波 信号的能量,如图3中实线所示;然后,求取能量的包 络线,如图3中虚线所示。由于系统计算资源有限, 为了减少计算量,采用最简单的线性插值的算法求取 包络。具体步骤为:找到能量信号的两个相邻的峰值 点( 。,E。)、( ,E ),根据式(3)即可得到插值点 ( ,E )的坐标。 — ,