CN110346004B - 一种双声道超声时差法的流量测量数据融合方法 - Google Patents

一种双声道超声时差法的流量测量数据融合方法 Download PDF

Info

Publication number
CN110346004B
CN110346004B CN201910756490.0A CN201910756490A CN110346004B CN 110346004 B CN110346004 B CN 110346004B CN 201910756490 A CN201910756490 A CN 201910756490A CN 110346004 B CN110346004 B CN 110346004B
Authority
CN
China
Prior art keywords
channel
flow
sound
time
data
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
Application number
CN201910756490.0A
Other languages
English (en)
Other versions
CN110346004A (zh
Inventor
蒋延付
王成李
卜勤超
章涛
赵伟国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Seck Intelligent Technology Co ltd
Original Assignee
Hangzhou Seck Intelligent Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Seck Intelligent Technology Co ltd filed Critical Hangzhou Seck Intelligent Technology Co ltd
Priority to CN201910756490.0A priority Critical patent/CN110346004B/zh
Publication of CN110346004A publication Critical patent/CN110346004A/zh
Application granted granted Critical
Publication of CN110346004B publication Critical patent/CN110346004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D1/00—Measuring arrangements giving results other than momentary value of variable, of general application
    • G01D1/02—Measuring arrangements giving results other than momentary value of variable, of general application giving mean values, e.g. root means square values
    • 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/66—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 measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
    • G01F1/667—Arrangements of transducers for ultrasonic flowmeters; Circuits for operating ultrasonic flowmeters
    • G01F1/668—Compensating or correcting for variations in velocity of sound

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Fluid Mechanics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明公开了一种双声道超声流量测量数据的融合方法。该发明首先采用格拉布斯准则对两个声道每秒测量的时差数据进行粗大误差剔除，然后分别将两个声道一秒内不包含粗大误差的时差数据进行流量计算并求平均流量，得到各个声道每秒的平均流量，最后采用改进的卡尔曼分布式融合方法将两个声道每秒的平均流量进行数据融合，得到双声道超声流量测量的平均流量。本发明的双声道超声流量测量数据的融合方法，能够减小测量过程中的随机误差，以及由于流场分布不均匀导致的流量测量误差，并且可以对双声路的两个声道工作状态进行判断，识别故障声道，剔除故障声道数据，有效地提高了双声道超声流量测量数据的准确性和可靠性。

Description

一种双声道超声时差法的流量测量数据融合方法
技术领域
本发明属于流量检测技术领域，涉及一种双声道超声时差法的流量测量数据融合方法。
背景技术
超声时差法的流量测量原理是通过检测测量管道中超声波信号顺逆流的渡越时间得到时间差，根据时差求得管道内的流体平均流量，该方法具有压损低、精度高、量程比宽和重复性好等优点。
基于时差法的超声波流量计通常按照声道分为单声道和多声道，小口径流量计受管道几何尺寸的限制一般设计为单声道。多声道超声波流量计应用于较大管径的测量，常用权重系数法、神经网络法等进行数据融合方法，对流场分布影响进行补偿。
对于中小管径的流量测量，常采用双声道的测量方法。然而，由于中小管径的超声波声程小，导致激励信号的干扰会严重影响回波信号，而且管径开孔的尺寸与管径之比较大，很容易在管道内产生扰流、漩涡，使回波信号发生幅值变化，流动噪声又会干扰回波信号，以上这些因素会干扰回波信号导致对顺逆流时间的测量产生误差或错误。如果采用两个声道权重方法进行数据融合，当某一声道受到干扰或产生故障时，系统仍按照设定的权重系数进行流量计算，会导致超声流量测量的数据偏差或错误，降低了超声流量测量的准确性和仪表的可靠性。
发明内容
本发明针对现有双声道超声流量测量数据的权重系数融合方法的不足，提出了一种具有故障自诊断功能的双声道超声流量测量数据的融合方法，该方法可以自动识别故障声道，剔除故障声道数据，能够提高超声流量测量数据的准确性和可靠性。
本发明的主要技术构思：首先，将各个声道每秒的时差数据根据格拉布斯准则进行粗大误差剔除；然后根据各个声道的时差与流量函数将这些时差数据转换为流量数据，并分别求出两个声道在一定时间内的平均流量；最后采用改进的卡尔曼分布式融合进行数据融合，融合过程中先进行双总体U检验，判断两个声道的工作状态，然后按照不同的方法计算得到双声道超声流量测量的平均流量。
本发明解决技术问题所采取的技术方案为：
本发明首先采用了格拉布斯准则对两个声道一定时间内测量的时差数据进行粗大误差剔除，提高两个声道时差数据的准确性。然后分别将两个声道不包含粗大误差的时差数据进行流量计算，分别得到各个声道平均流量，减少了时差测量过程中的随机误差。最后采用改进的卡尔曼分布式融合方法将两个声道平均流量进行数据融合，得到双声道超声流量测量的平均流量。改进的卡尔曼分布式融合方法在原有的卡尔曼分布式融合方法的基础上进行改进，首先将两个声道流量数据分别进行卡尔曼滤波，然后进行局部估计融合得到管道内平均流量，其中在局部估计融合过程中进行双总体U检验，判断两个声道的工作状态。若某一声道产生故障，则剔除该声道数据保证流量测量的准确性。
本发明的有益效果在于：本发明的双声道超声流量测量数据的融合方法，采用了格拉布斯准则对两个声道每秒测量的50个时差数据进行粗大误差剔除，降低了中小管径测量中激励信号对回波到达时间的测量影响，提高了两个声道时差数据的准确性。然后分别将两个声道一秒内不包含粗大误差的多个时差数据进行流量计算并求平均流量，减小了测量过程中的随机误差，提高了流量测量的动态测量性能。为了减少在管道内扰流、漩涡对回波信号幅值变化的影响及回波信号波形的干扰，采用改进的卡尔曼分布式融合方法，在局部估计融合过程中进行双总体U检验，判断两个声道的工作状态，进行相应的流量计算，从而融合得到管道内平均流量。因此，本发明提出的方法，不仅减少了由于流场分布不均匀导致的流量测量误差，并且可以对双声道的两个声道工作状态进行判断，识别故障声道，剔除故障声道数据，该方法提高了双声道超声流量测量数据的准确性和可靠性。
附图说明
图1为双声道超声流量测量的声道分布示意图；
图2为双声道超声流量测量的时差数据融合方法流程图；
图3为改进的卡尔曼分布式融合方法流程图。
具体实施方式
以下结合附图进一步说明本发明提出的一种双声道超声流量测量数据的融合方法。
参照图1，为双声道超声流量测量的声道分布示意图。换能器A和换能器B构成声道一，超声波传播路径长度为L1，方向的夹角为θ1，换能器C和换能器D构成声道二，超声波传播路径长度为L2，方向的夹角为θ2，两个声道为“X”型分布。由图1可分别求得声道一顺流时间TAB，如公式(1)；逆流时间TBA，如公式(2)；声道二顺流时间TCD，如公式(3)；逆流时间TDC，如公式(4)：
Figure BDA0002168878820000031
Figure BDA0002168878820000032
Figure BDA0002168878820000033
Figure BDA0002168878820000034
其中c为超声波在流场中的传播速度，V为流体的流速，则可以求得两个声道的时差ΔT1、ΔT2分别与流速的函数关系，如公式(5)、(6)：
Figure BDA0002168878820000035
Figure BDA0002168878820000036
由于c2＞＞(Vcosθ1)2和c2＞＞(Vcosθ2)2，则式(5)和(6)可简化为(7)和(8)：
Figure BDA0002168878820000037
Figure BDA0002168878820000038
则两个声道可以根据管道截面S和测量得到的时差ΔT1、ΔT2分别导出流量计算函数：
Figure BDA0002168878820000039
Figure BDA0002168878820000041
其中k1，k2分别为声道一和声道二的流量修正系数。
图2为双声道超声流量测量的时差数据融合方法流程图。超声流量测量系统设定每秒两个声道分别可以测量50次，即每秒两个声道可以分别得到50个时差数据。首先两个声道分别进行粗大误差剔除，然后根据式(9)和式(10)将两个声道不包含粗大误差的时差数据分别计算得到流量并取平均值，分别得到两个声道测量的平均流量
Figure BDA0002168878820000047
和
Figure BDA0002168878820000046
最后采用改进的卡尔曼分布式融合方法进行数据融合，得到每秒双声道超声流量测量的平均流量Q。根据样本大小选择合适的粗大误差剔除准则，由于两个声道的时差数据量大小均为50，因此本发明采用格拉布斯准则进行粗大误差剔除。
图3为改进的卡尔曼分布式融合方法流程图。首先两个声道分别进行卡尔曼滤波得到局部估计值，然后进行局部估计值融合。
在双声道超声流量测量中，系统的状态方程和观测方程为：
Figure BDA0002168878820000042
其中x(t)表示t时刻管道内的流量，q(t)表示t时刻被测系统的系统噪声，
Figure BDA0002168878820000043
表示t时刻第i个声道的平均流量，ri(t)表示t时刻第i声道的测量噪声。则可导出卡尔曼迭代方程组：
Figure BDA0002168878820000044
其中Pi(t)为t时刻第i声道的误差协方差，Ki(t+1)为t+1时刻的卡尔曼增益。
经过卡尔曼滤波后可以分别得到两个声道的估计量x1(t)和x2(t)，根据双总体U检验判断两个声道的工作作态选择对应的局部估计值融合策略。双总体U检验的检验统计量为：
Figure BDA0002168878820000051
其中
Figure BDA0002168878820000052
和
Figure BDA0002168878820000053
分别为样本一和样本二的平均值，
Figure BDA0002168878820000054
和
Figure BDA0002168878820000055
分别为样本一和样本二的方差，n1和n2分别为样本一和样本二的大小。则可以分别得到声道一t时刻不包含粗大误差的流量数组与其t-1时刻不包含粗大误差的流量数组的检验统计量U11(t)、声道二t时刻不包含粗大误差的流量数组与其t-1时刻不包含粗大误差的流量数组的检验统计量U22(t)、声道一和声道二t时刻不包含粗大误差的流量数组的检验统计量U12(t)。根据U11(t)、U22(t)和U12(t)可以得到双声道超声流量计测量状态如下表所示：
表1双声道超声流量计测量状态判定表
Figure BDA0002168878820000056
通过表1判断超声流量计测量状态，识别产生故障的声道，剔除该声道的测量数据以保证滤波的准确性，则局部值估计融合式为：
Figure BDA0002168878820000057
其中ε表示显著性水平α为0.01时双总体U检验的临界值。最终输出的Q(t)作为t时刻双声道超声流量计测量的平均流量。

Claims (2)

1.一种双声道超声时差法的流量测量数据融合方法，其特征在于：采用格拉布斯准则对两个声道每秒测量的时差数据进行粗大误差剔除，然后分别将两个声道内不包含粗大误差的时差数据进行流量计算并求平均流量，得到各个声道平均流量，最后采用改进的卡尔曼分布式融合方法将两个声道平均流量进行数据融合，得到双声道超声流量测量的平均流量；
所述的改进的卡尔曼分布式融合方法具体是：首先将两个声道流量数据分别进行卡尔曼滤波，然后进行局部估计融合得到管道内平均流量，其中在局部估计融合过程中进行双总体U检验，判断两个声道的工作状态，进而采用不同的数据融合。
2.根据权利要求1所述的一种双声道超声时差法的流量测量数据融合方法，其特征在于，设两个声道的估计量x1(t)和x2(t)，通过双总体U检验得到声道一t时刻不包含粗大误差的流量数组与其t-1时刻不包含粗大误差的流量数组的检验统计量U11(t)、声道二t时刻不包含粗大误差的流量数组与其t-1时刻不包含粗大误差的流量数组的检验统计量U22(t)、声道一和声道二t时刻不包含粗大误差的流量数组的检验统计量U12(t)；则数据融合方式为：
Figure FDA0002505636560000011
其中ε表示显著性水平α为0.01时双总体U检验的临界值，P1(t)为t时刻声道一的误差协方差，P2(t)为t时刻声道二的误差协方差，最终输出的Q(t)作为t时刻双声道超声流量计测量的平均流量。
CN201910756490.0A 2019-08-16 2019-08-16 一种双声道超声时差法的流量测量数据融合方法 Active CN110346004B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910756490.0A CN110346004B (zh) 2019-08-16 2019-08-16 一种双声道超声时差法的流量测量数据融合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910756490.0A CN110346004B (zh) 2019-08-16 2019-08-16 一种双声道超声时差法的流量测量数据融合方法

Publications (2)

Publication Number Publication Date
CN110346004A CN110346004A (zh) 2019-10-18
CN110346004B true CN110346004B (zh) 2020-08-21

Family

ID=68185233

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910756490.0A Active CN110346004B (zh) 2019-08-16 2019-08-16 一种双声道超声时差法的流量测量数据融合方法

Country Status (1)

Country Link
CN (1) CN110346004B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111426353B (zh) * 2020-04-08 2022-02-11 中国民用航空飞行学院 一种精确流量获取装置及方法
CN111982248B (zh) * 2020-08-14 2022-12-06 中国计量大学 一种改进的传感器动态数据融合方法
CN112697215B (zh) * 2020-11-19 2022-10-18 浙江工业大学 一种用于超声波水表数据滤波的卡尔曼滤波参数调试方法
CN112484798B (zh) * 2020-11-27 2024-03-19 中科水研(江西)科技股份有限公司 一种基于改进卡尔曼滤波器的时差超声波流量测量方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000042995A2 (en) * 1999-01-21 2000-07-27 Ortho-Mcneil Pharmaceutical, Inc. Use of anticonvulsant derivatives for treating transformed migraine
WO2003093499A2 (en) * 2002-05-03 2003-11-13 Pfizer Products Inc. Therapeutic use of selective pde10 inhibitors
CN101004352A (zh) * 2006-01-20 2007-07-25 深圳市建恒工业自控系统有限公司 多声道流量测量装置及其方法
CN101701828A (zh) * 2009-11-23 2010-05-05 常州达奇信息科技有限公司 基于立体视觉和信息融合的盲人自主导航方法
CN107367306A (zh) * 2017-08-31 2017-11-21 合肥工业大学 一种基于Kalman滤波的气体超声波流量计融合方法和系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003194604A (ja) * 2001-12-26 2003-07-09 Babcock Hitachi Kk 音波式ガス温度・流量計とガス中音波の伝播時間の検出方法
CN102981160B (zh) * 2012-11-08 2014-05-07 中国兵器科学研究院 一种确定空中目标航迹的方法及装置
CN203535444U (zh) * 2013-09-28 2014-04-09 西安航空学院 基于多传感器数据融合的污水处理远程监控系统
CN103870686B (zh) * 2014-03-07 2015-08-12 中国民航大学 一种基于信息融合的超声相控阵方法
CN105352535A (zh) * 2015-09-29 2016-02-24 河海大学 一种基于多传感器数据融合的测量方法
US10160448B2 (en) * 2016-11-08 2018-12-25 Ford Global Technologies, Llc Object tracking using sensor fusion within a probabilistic framework
CN209148563U (zh) * 2018-10-25 2019-07-23 广东华中科技大学工业技术研究院 一种双成像式透明构件缺陷检测装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000042995A2 (en) * 1999-01-21 2000-07-27 Ortho-Mcneil Pharmaceutical, Inc. Use of anticonvulsant derivatives for treating transformed migraine
WO2003093499A2 (en) * 2002-05-03 2003-11-13 Pfizer Products Inc. Therapeutic use of selective pde10 inhibitors
CN101004352A (zh) * 2006-01-20 2007-07-25 深圳市建恒工业自控系统有限公司 多声道流量测量装置及其方法
CN101701828A (zh) * 2009-11-23 2010-05-05 常州达奇信息科技有限公司 基于立体视觉和信息融合的盲人自主导航方法
CN107367306A (zh) * 2017-08-31 2017-11-21 合肥工业大学 一种基于Kalman滤波的气体超声波流量计融合方法和系统

Also Published As

Publication number Publication date
CN110346004A (zh) 2019-10-18

Similar Documents

Publication Publication Date Title
CN110346004B (zh) 一种双声道超声时差法的流量测量数据融合方法
RU2446393C2 (ru) Способ диагностики шероховатости трубопровода и ультразвуковой расходомер
CN103344288B (zh) 一种基于零点分析的时差式超声波流量计测量方法
CN101886939A (zh) 一种时差法超声流量计静态漂移抑制模型及抑制方法
CN102692627B (zh) 用于处理多普勒测量信号的仪器和方法
CN106643939A (zh) 用于超声波流量计计算超声波传播时间的方法
CN104089186B (zh) 一种基于组合滤波和动态阈值的管道压力异常诊断方法
CN105181793B (zh) 基于超声波双频信号测量两相流含气率的方法
JP2017502319A (ja) 自己チェック機能のある流量計および方法
CN201795819U (zh) 一种时差法超声流量计静态漂移抑制模型
CN111982248B (zh) 一种改进的传感器动态数据融合方法
CN112304376B (zh) 基于数据融合的超声波流量计流量测量方法
CN102095889B (zh) 三通道超声时差流速测量方法
CN117740090A (zh) 一种超声波的杂波消除方法、系统及流量计
CN117405187A (zh) 一种四声道气体超声波流量计的声道异常筛选方法
CN106678553B (zh) 一种泄漏动态压力波在管内气体中传播速度的计算方法
RU2376581C2 (ru) Способ и устройство для получения информации о распределении макроскопических частиц в жидкости по размеру
CN111964738A (zh) 一种基于超声回波信号相似度的动态阈值设置方法
CN115308434A (zh) 一种超声波测速方法及系统
CN111473827B (zh) V形声道零飘消除方法
CN103207011A (zh) 一种基于电路相移的超声精确时延检测方法
JP3604322B2 (ja) 超音波流量計
CN114674384A (zh) 超声波流量计的防错波检测方法、装置、设备及流量计
Li et al. A novel differential time-of-flight algorithm for high-precision ultrasonic gas flow measurement
CN112466330A (zh) 一种多途信道下噪声源的声源级估计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant