一种噪声背景下准确估计超声流量计传播时间方法
【专利说明】一种噪声背景下准确估计超声流量计传播时间方法 所属技术领域
[0001] 本发明属于流量测量技术领域，涉及超声流量计。
【背景技术】
[0002] 超声流量计作为新兴的流量仪表近十几年发展迅速，在大型水利水电项目以及天 然气贸易结算领域已有广泛应用，其中时差法超声流量计作为计量仪表应用最多。精确确 定超声换能器接收信号的起始位置，得到精确的超声波传播时间是时差法超声流量计能够 准确测量流量的关键。目前，传统测量方法中普遍采用的是双门限电平法（简称双门限 法），首先通过非过零比较，确定使用哪个周期的信号作为计时起始位置，再通过过零比较 确定超声波传播时间。但这种方法的缺点是当有噪声存在时，过零比较得到的时间会有较 大波动，严重影响了流量计的精度，而且当噪声大到一定程度时，非过零比较无法准确确定 信号位置，从而导致错波发生，此时由于误差太大，流量计已不能使用。
[0003] 为了在有噪声时仍能获得准确的超声波传播时间，研宄者提出基于超声脉冲接收 信号模型对超声波传播时间进行估计的方法。目前超声脉冲接收信号的经验模型主要有两 种，高斯模型和混合指数模型。高斯模型用于描述比较对称的宽带高频信号。周方等学者 将此模型与高斯一牛顿法、模拟退火法、蚁群算法等相结合，证明了蚁群算法是一种有效的 回波参数估计方法，具有较高精度[1];针对混合指数模型，姚振静等人将此模型与卡尔曼 滤波及其优化算法进行结合求取超声波传播时间，用来提高超声测距精度[2]。虽然已有学 者基于超声脉冲接收信号模型提出了一些超声波传播时间估计算法，但研宄仅针对算法本 身，应用也仅限于超声测距领域。对于超声流量测量来说，时间测量精度要求远高于超声测 距，已有方法不能直接用来解决噪声背景下超声流量计的准确测时问题。
[0004] [1]周方，张小凤，张光斌，超声回波参数的蚁群算法估计，陕西师范大学学报（自 然科学版），2012,40 (2)，35 ~40
[0005] [2]Zhen-JingYao,Qing-HaoMeng,MingZeng,Improvementintheaccuracyof estimatingthetime-〇f-flightinanultrasonicrangingsystemusingmultiple square-rootunscentedKalmanfilters,ReviewofScientificInstruments，2010,81， 104901-1 ~104901-7
【发明内容】
[0006] 本发明为解决噪声背景下超声流量计中超声波传播时间准确测量问题，基于超声 脉冲接收信号混合指数模型，提供一种超声波传播时间准确估计方法，达到降低测量误差， 提高超声流量计测量精度的目的。
[0007] 本发明的技术方案如下：
[0008] -种噪声背景下准确估计超声流量计传播时间方法，包括下列步骤：：
[0009] 步骤一：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型如 下：
[0010] Am (x) =A(kts)sin[2 3ifc (kts-t) + 0 ] (1)
[0011] 其中，
[0012] 式中，Am(X)为超声脉冲接收离散模型信号，k为采样点数序号（k= 1,2，…N)，N 为采样点数，ts为采样时间间隔，f。为超声换能器的中心频率，0为初始相位角，A。为接收 信号幅值，T和m是超声换能器的两个特性参数，t为超声波传播时间，U(kts-T)为单位 阶跃信号，当采样频率、采样点数、超声换能器及其发射频率确定后，模型中的k、N、ts、f。和 9即为定值，随着被测流量改变，模型信号Am(X)中的和t会相应变化，其中，参数 向量x=[心mTt]，对x估计得越准确，模型信号越逼近实验获得的真实接收信号，利用 下面的遗传一蚁群算法对x中这四个参数进行估计，最终t的最优估计值即为超声波传播 时间；
[0013] 步骤二：根据实验采样获得的超声接收信号与模型信号Am(X)确定目标函数 f(x):基于最小二乘思想确定目标函数，将4中各采样点电压与Am(X)中采样点电压之差的 平方和作为目标函数f(x)，当目标函数f(x)取得最小值时，x为最优估计值；
[0014] 步骤三：基于目标函数f(x)构造适应度函数eval(x)，目标函数值越小，适应度函 数值就越大：
[0016] 其中，a是使得a_f(Xi)为正的任意实数，n为种群数目，n值设定越大，收敛精度 越高，收敛速度越慢，角标i代表第i个种群个体，在参数边界内随机产生；
[0017] 步骤四：设定x中四个参数边界，根据适应度函数值确定每个个体被复制的次数： 首先确定每个个体的预复制次数，它等于相应个体的适应度函数值与种群数目乘积的整数 部分，设R为总的预复制次数与种群数目之差，为了保持种群数目不变，从最优个体到第R 个优秀个体，每个个体的复制数目加1，即适应度函数值较大的个体，被复制的概率越大；
[0018] 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个 体分为一组，选定一个固定的判断概率〇. 9,并随机生成一个在0到1之间的实数P。，如果 P。小于判断概率〇. 9则进行交叉，否则不交叉，为获得最优个体，每次交叉、变异之后都用上 一代的最优个体代替下一代的最差个体；步骤二至步骤五为遗传算法，通过迭代获得优化 解作为后续蚁群算法的信息素初始条件；
[0019] 步骤六：利用上述的遗传算法得到的优化解初始化蚁群信息素；
[0020]步骤七：蚂蚁移动步长a初始化，根据实验获得的超声接收信号确定参数向量X 的变化范围；
[0021] 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参 数向量X中的每个参数单独进行，首先对Ac!进行，令
，如果/?) 小于/(41)，则
，如果/(〇小于
;其中，；^为第t代第s只蚂蚁的参数向量 值，xf1为第t_l代迭代之后按照信息素更新第s只蚂蚁最好的参数向量值，(1\在[_a，a] 之间随机产生；同理，对T、m和t三个参数进行优化，得到新的优化解；
[0022] 步骤九：信息素挥发和更新：按照下列规则设定信息素挥发系数：当迭代次数较 小时，信息素挥发系数设定较大，有利于全局寻优；当迭代次数较大时，信息素挥发系数设 定较小，有利于提高局部收敛精度、加快收敛速度；重新计算信息素，蚂蚁重新分配，返回步 骤七进行迭代计算，迭代结束判据设为：在算法迭代t代后进行判断，如果第t代的最优目 标函数值与第t-j(j= 1，2, 3, 4, 5)代的最优目标函数值之差的绝对值均小于1(T7,则迭代 结束，此时参数向量x为最优解，其中的t即为超声波传播时间最优估计值。
[0023] 本发明将遗传算法和蚁群算法相结合，利用遗传算法的全局寻优能力使得在噪声 背景下仍能找到全局最优解，且遗传算法收敛速度比蚁群算法快。同时结合蚁群算法的高 度局部寻优能力，提高参数估计的精确度。
【具体实施方式】
[0024] 以下将以实验获得的气体超声接收信号为目标对象，结合技术方案中的步骤方 法，给出各步的参数设置和实施方法：
[0025] 步骤一 ：DN100单声道气体超声流量计的超声换能器中心频率为120KHZ，实验采 用泰克TBS1102示波器采集超声换能器接收到的原始信号，时基选择lOus/div，整个示波 器窗口采样时间为l〇us/div*10格=lOOus，在2500个数据点存储深度下，相邻采样点的时 间间隔为lOOus/2. 5k= 40ns，则超声脉冲接收信号混合指数模型[2]中（式（1) (2))的相 关参数设置为：采样点序号k= 1，2,…2500,采样时间间隔ts= 40ns，换能器中心频率 = 120KHz，初始相位角0=0rad。其余参数AQ、T、m和t为待估计参数，定义为参数向量x= [A0mTt] 0
[0026] 步骤二：根据实验采样获得的超声接收信号与模型信号Am(X)确定目标函数 f(x) 〇
[0028] 其中，是示波器采样获得的真实超声接收信号，Am(X)是基于超声脉冲接收信号 混合指数模型重构的信号，x= [A。mTt]为待估计参数向量，采样点数N= 2500,角标 k代表在第k个采样点处的信号。
[0029] 步骤三：基于目标函数f(x)构造适应度函数eval(x)。式⑶中，设定a= 100， 种群数目n= 20,种群个体由Matlab中rand函数随机生成。
[0030] 步骤四：通过分析示波器采样超声接收信号特点，将x中四个参数边界分别设为 A0G[0? 025,0. 045]V，mG[3.4,3. 8]，TG[12. 3, 12.7]，tG[384, 439]us。根据适应度 函数值确定每个个体被复制的次数乂：
[0031] Mj=abs(eval(xj) *n)
[0033] Mi是第i个种群个体的预复制数目，R是种群数目与总的预复制数目之差。为了 保持种群数目不变，从最优个体到第R个优秀个体，每个个体的复制数目加1。
[0034] 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个 体分为一组，选定一个固定的判断概率〇. 9,采用Matlab中rand函数随机生成一个