一种基于超声回波信号相似度的动态阈值设置方法技术

本发明专利技术涉及一种基于超声回波信号相似度的动态阈值设置方法。本发明专利技术通过评估设定工况和实际工况下回波信号上升区域部分的相似度，进而根据最大相似度确定阈值电压，实现阈值的准确调整，确保了回波信号到达时刻点的准确定位，实现飞行时间的准确测量，提高了气体超声流量计的测量准确性与稳定性。

【技术实现步骤摘要】
一种基于超声回波信号相似度的动态阈值设置方法
本专利技术属于流量检测
，涉及一种基于超声回波信号相似度的动态阈值设置方法。
技术介绍
气体超声流量计具有可双向测量性、量程比高、重复性好、低压损等优势，已广泛应用于气体流量测量。气体超声流量计原理是根据测量声波在流体介质传播的顺逆流飞行的时间差计算得到截面的平均流速，其关键在于对顺逆流的飞行时间测量，飞行时间的测量准确性直接影响了气体超声流量计的计量精度。目前针对顺逆流飞行时间的测量方法，主要包括互相关法、曲线拟合法、双阈值法。双阈值法采用第一阈值检测和过零检测确认回波到达，得到回波渡越时间；该方法因原理简单，无需繁杂的计算量，硬件电路容易实现而被广泛应用。但是，双阈值法的缺点也很明显，其主要表现在该方法对回波信号的稳定性要求高。然而，超声回波信号在不同流速气体介质中衰减程度不同，回波信号幅值不稳定且回波特性易受环境因素影响，导致回波波形发生变化，这易导致阈值与回波信号产生误触发，得到错误的飞行时间到达点，导致回波渡越时间的误测量。
技术实现思路
本专利技术针对现有双阈值法的不足，提出了一种基于超声回波信号相似度的动态阈值设置方法。该方法通过评估设定工况和实际工况下回波信号上升区域部分的相似度，进而根据最大相似度确定阈值电压，实现阈值的准确调整，确保了回波信号到达时刻点的准确定位，实现飞行时间的准确测量。本专利技术的基于超声回波信号相似度的动态阈值设置方法具体步骤如下：步骤一：在设定工况下设置合理的阈值，设置的阈值幅值为回波的第2极大值与第3极大值的电压平均值，记为V[2,3]。并记录回波信号的第2,3,4极大值点的电压幅值，记为Vb2，Vb3，Vb4。步骤二：记录实际工况下回波信号的第1,2,3,4,5极大值点的电压幅值，分别记为Vc1，Vc2，Vc3，Vc4，Vc5。步骤三：对记录的实际工况下的极大值点电压幅值划分为3组，分别为a组：Vc1，Vc2，Vc3；b组：Vc2，Vc3，Vc4；c组：Vc3，Vc4，Vc5；然后分别计算每组数据的第一个值与Vb2的比例系数，如下：根据得到的比例系数对数据(Vb2,Vb3,Vb4)进行比例放缩处理，如下：Ven＝k1*Vbn,(n＝2,3,4)(4)Vfn＝k2*Vbn,(n＝2,3,4)(5)Vgn＝k3*Vbn,(n＝2,3,4)(6)可获得新的3组数据，分别为I组：Ve2，Ve3，Ve4；分别为II组：Vf2，Vf3，Vf4；分别为III组：Vg2，Vg3，Vg4；步骤四：分别计算a组与I组，b组与II组，c组与III组的欧式距离，用于评估相似程度，公式如下：步骤五：根据欧式距离评估的相似度进行准确的阈值调整。式(10)中，VThreshold代表实际工况下需要设置的阈值电压；V[1,2]代表实际工况下回波信号的第1极大值与第2极大值的电压平均值，V[2,3]与V[3,4]同理。min{}函数代表求取括号内数据的最小值，代表最高的相似度。本专利技术的有益效果在于：本专利技术能对由工况变化引起回波波形变化的回波进行准确地阈值调整设置，确保了飞行时间测量的准确性，提高了气体超声流量计的测量准确性与稳定性。附图说明图1为基于双阈值法的飞行时间测量原理图；图2为工况环境变化影响飞行时间测量示意图；图3为两回波信号部分区域之间相似度评估示意图；图4为设定工况下回波信号的阈值设置示意图；图5为根据相似度评估结果调整阈值幅值示意图；图6为基于超声回波信号相似度的动态阈值设置方法流程图；图7为本专利技术方法实施中使用的超声波流量计硬件电路框图。具体实施方式以下结合附图进一步说明本专利技术提出的一种基于超声回波信号相似度的动态阈值设置方法。图1为基于双阈值法的飞行时间测量原理图，图中第一阈值线用于阈值检测，第二阈值线用于过零检测。第一阈值与回波相交的黑色圆形点为特征点，根据特征点位置和第二阈值线与回波相交可得图中黑色的方形过零点，该点作为回波信号到达时刻点。由激励的产生做计时的起始，回波信号到达时刻点作为计时的终止，通过计时器件可以得到回波由发射到被接收的到达时间Ta，从而得到飞行时间Tf。到达时间Ta与飞行时间Tf相差n个周期时间Tn。Tf＝Ta-Tn(11)图2为工况环境变化影响飞行时间测量示意图。图中黑色实线为常压101kPa压力下的回波波形，黑色虚线为500kPa压力下的回波波形，可以看出随压力增大，回波波形发生明显变化。在常压101kPa下，根据阈值触发确定飞行时间到达点，得到的飞行时间为T1。而随着压力变化引起波形变化，阈值触发确定的飞行时间到达点较常压下往前推移一个周期，得到飞行时间为T2。飞行时间T1与T2的不同反应了：随工况环境变化，回波波形也发生了变化，回波波形的变化导致阈值的误触发影响了飞行时间的测量。且相比101kPa压力下的回波信号，500kPa压力下的回波信号最前端还有一个周期的小波凸起，这种情况使得基于回波极值点分布进行阈值调整的方法将不再适用。图3为两回波信号部分区域之间相似度评估示意图。本方法采用欧式距离评估2个回波上升区域的3个周期波形曲线之间的相似程度。其定义如下：式(12)中，I和J代表两个不同个体，(xi1,xi2,...,xin)与(xj1,xj2,...,xjn)为代表I、J个体的n维数值特征数据。结合图3，I、J个体是指黑色虚线框内回波信号的3个周期波形。其中的n维数值特征数据是指对黑色虚线框内提取的回波极大值，再经比例放缩处理后的数据。如图3中，星号标记处为回波信号的极大值。计算得到d(I,J)距离代表了两黑色虚线框内回波波形之间的相似程度，距离数值越小，相似度越高。若通过改变J个体所对应的回波上升区域的部分波形，分别与I个体进行欧式距离计算，可根据最小的欧式距离确认2个回波上升区域的部分波形之间的相互对应关系。图6为基于超声回波信号相似度的动态阈值设置方法流程图。本段同时结合图4与图5描述所提方法的流程步骤。详细方法步骤如下：步骤1：在设定工况下设置合理的阈值与回波特征数据获取。如图4虚线代表的阈值幅值为回波的第2极大值与第3极大值的电压平均值，记为V[2,3]，并记录回波信号的第2,3,4极大值点的电压幅值，作为对应周期波形的特征数据，记为Vb2，Vb3，Vb4。图4中，2,3,4标号台阶对应为回波信号的第2,3,4极大值幅值。步骤2：实际工况下特征数据获取。同步骤1方式记录实际工况下回波信号的第1,2,3,4,5极大值点的电压幅值，作为对应周期波形的特征数据，分别记为Vc1，Vc2，Vc3，Vc4，Vc5。步骤3：特征数据处理。对记录的实际工况下的极大值点电压幅值划分为3组，分别为a组：Vc1，Vc2，Vc3；b组：Vc2，Vc3，Vc4；c组：Vc3，V本文档来自技高网...

【技术保护点】
1.一种基于超声回波信号相似度的动态阈值设置方法，其特征在于该方法包括以下步骤：/n步骤一：在设定工况下设置合理的阈值，设置的阈值幅值为回波的第2极大值与第3极大值的电压平均值，记为V

【技术特征摘要】
1.一种基于超声回波信号相似度的动态阈值设置方法，其特征在于该方法包括以下步骤：
步骤一：在设定工况下设置合理的阈值，设置的阈值幅值为回波的第2极大值与第3极大值的电压平均值，记为V[2,3]；并记录回波信号的第2,3,4极大值点的电压幅值，记为Vb2，Vb3，Vb4；
步骤二：记录实际工况下回波信号的第1,2,3,4,5极大值点的电压幅值，分别记为Vc1，Vc2，Vc3，Vc4，Vc5；
步骤三：对记录的实际工况下的极大值点电压幅值划分为3组，分别为a组：Vc1，Vc2，Vc3；b组：Vc2，Vc3，Vc4；c组：Vc3，Vc4，Vc5；然后分别计算每组数据的第一个值与Vb2的比例系数，如下：









根据得到的比例系数对数据(Vb2,Vb3,Vb4)进行比例放缩处理，如下：
Ven＝k1*Vbn(4)

【专利技术属性】
技术研发人员：马也驰，赵伟国，章圣意，
申请(专利权)人：中国计量大学，
类型：发明
国别省市：浙江;33