CN107270983B - 基于瞬时流量的数据搭建超声波流量模型的方法 - Google Patents
基于瞬时流量的数据搭建超声波流量模型的方法 Download PDFInfo
- Publication number
- CN107270983B CN107270983B CN201710540868.4A CN201710540868A CN107270983B CN 107270983 B CN107270983 B CN 107270983B CN 201710540868 A CN201710540868 A CN 201710540868A CN 107270983 B CN107270983 B CN 107270983B
- Authority
- CN
- China
- Prior art keywords
- flow
- data
- time difference
- instantaneous flow
- instantaneous
- 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
Images
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/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
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于瞬时流量的数据搭建超声波流量模型的方法。包括如下步骤:(1)产生稳定瞬时流量:利用瞬时流量产生装置产生需要的流量V;(2)时间差采样:在稳定的瞬时流量下,超声波计量装置采集一定数量的时间差数据t1....tn,所述时间差是指超声波在两点之间的逆流方向的传播时间和顺流方向的传播时间之差;(3)数据处理:对采集到的时间差数据进行数字滤波算法得到数据t;(4)数据模型:根据上述的瞬时流量V和时间差数据t通过公式A=V/t计算获得瞬时流量和时间差的流量系数A;(5)大样本数据模型:改变瞬时流量,重复步骤(1)、(2)、(3)、(4),从而得到多个流量系数A,通过数据拟合方法,得到时间差t与瞬时流量V之间的关系:V=A(t)。
Description
技术领域
本发明涉及一种基于瞬时流量的数据搭建超声波流量模型的方法,属于超声波流量计量技术领域。
背景技术
超声波流量计量是目前工业、民用流量计量中应用较广泛的一种流量计量方法。超声波流量计具有计量精度高、量程比高、使动流量低等特点,因此超声波流量计量装置越来越受到广泛的欢迎。目前的超声波流量模型是按照公式进行计算所得,公式中涉及到粘度、流量分布、管径大小等流体因素,不同的口径需要调整各种参数,然后根据参数获得模型公式。实验证明这种模型误差较大,因此需要大量的测试和修正,这种方法影响流量建模的效率。
发明内容
针对现有技术中存在的上述缺陷,本发明提供了一种简单、建模时间短、模型更加准确、超声波流量计量装置精度更高的基于瞬时流量的数据搭建超声波流量模型的方法。
本发明是通过如下技术方案来实现的:一种基于瞬时流量的数据搭建超声波流量模型的方法,其特征是:包括如下步骤:
(1)产生稳定瞬时流量:利用瞬时流量产生装置产生需要的流量V;
(2)时间差采样:在稳定的瞬时流量下,超声波流量计量装置采集一定数量的时间差数据t1....tn,所述时间差是指超声波在两点之间的逆流方向的传播时间和顺流方向的传播时间之差;
(3)数据处理:对采集到的时间差数据进行数字滤波算法得到数据t;
(4)数据模型:根据上述的瞬时流量V和时间差数据t通过公式A=V/t计算获得瞬时流量和时间差的流量系数A;
(5)大样本数据模型:改变瞬时流量,重复步骤(1)、(2)、(3)、(4),从而得到多个流量系数A,通过数据拟合方法,得到时间差t与瞬时流量V之间的关系:V=A(t)。
进一步的,步骤(1)中采集的流量V是在瞬时流量产生装置达到稳流状态之后进行采集。
进一步的,步骤(3)中,数据处理时采用均值法或者去掉最大值和最小值得均值的方法。
本发明的有益效果是:本发明的方法简单,完全摒弃了公式模型中粘度、流量分布、管径大小等因素对模型的影响,模型完全由测量数据获得,模型更加准确,超声波流量计量装置精度更高,且本发明的方法只需要产生短时间的稳定的瞬时流量,因此建模时间相比传统的质量法和体积法大大减小。
附图说明
图1是具体实施方式中的超声波流量计量装置的结构模型示意图;
图2是具体实施方式中的A与t的关系曲线图;
具体实施方式
下面通过非限定性的实施例对本发明作进一步的说明:
一种基于瞬时流量的数据搭建超声波流量模型的方法,其包括如下步骤:
步骤一:产生稳定瞬时流量:利用瞬时流量产生装置产生需要的流量V,且采集的流量V是在瞬时流量产生装置经过一定时间后达到稳流状态之后进行采集,在未达到稳流期间不计量或者计量无效。
步骤二:时间差采样:在稳定的瞬时流量下,超声波流量计量装置采集一定数量的时间差数据t1....tn,所述时间差是指超声波在两点之间的逆流方向的传播时间和顺流方向的传播时间之差。采样时,超声波流量计量装置接收到采样开始指令后开始采样,接收到停止采样后停止采样。本发明采用的超声波计量原理是时间差法,超声波在介质中的传播受到介质流动的影响,在逆流方向上的传播时间和顺流方向上的传播时间之差称为时间差,时间差的大小可以反映介质流速的大小。由附图1所示的超声波流量计量装置的结构模型中,A、B为换能器,Vm为介质的流速,D为管道的直径,φ为换能器A和换能器B相连接与管道轴心线的夹角,L为换能器A和换能器B之间的直线距离。
目前超声波流量计量模型按照以下公式计算所得:设超声波在介质中的传播速度为C,那么,超声波信号从换能器A传播到换能器B的时间即顺流时间为:超声波信号从换能器B传播到换能器A的时间即逆流时间为:则,时间差t=T逆-T顺,一般C远大于Vm,则时间差为 因此可得:上述公式为理论计算而得,实际流速Vm仍然受介质粘度、流体分布、管道结构等影响。
步骤三:数据处理:对采集到的时间差数据进行数字滤波算法得到数据t;数字滤波算法可以采用:算术平均值滤波、加权平均值滤波、滑动平均滤波算法等。以算术平均值滤波为例:
步骤四:数据模型:根据上述的瞬时流量V和时间差数据t通过公式A=V/t计算获得瞬时流量和时间差的流量系数A;
步骤五:大样本数据模型:改变瞬时流量,重复步骤(1)、(2)、(3)、(4),从而得到多个流量系数A,通过数据拟合方法,得到时间差t与瞬时流量V之间的关系:V=A(t)。
下面以DN20口径超声波水表为例,采集以下数据并计算流量系数A进行说明:
根据上表内的数据作出A与t的关系曲线图,如附图2所示。
根据A与t的关系曲线的变化趋势,通过数据拟合方法,得到时间差t与瞬时流量V之间的关系:V=A(t)。曲线拟合函数有很多,如线性函数、多项式函数、指数函数、三角函数。如附图2中的曲线,为了更加准确的进行曲线拟合,对上述曲线进行分段,分为三段:50L/h-110.4L/h、110.4L/h-361.7L/h、361.7L/h-1808.6L/h三段,根据每一段的走势,分别的模型函数为:
第一段50L/h-110.4L/h模型函数A=at+b
第二段110.4L/h-361.7L/h模型函数A=a(t+b)2
第三段361.7L/h-1808.6L/h模型函数A=a(t+b)2
由上述模型函数可知,每个函数只有两个参数,为了简化拟合函数参数的计算,把每一段的起始点和终止点的时间差代入上述函数中,分别求取函数参数a和b的值。
计算结果为:
第一段拟合函数:A=0.06267t+6.696
第二段拟合函数:A=-0.0004623(t-45)2+8.04
第三段拟合函数:A=-0.000012(t-215.7)2+8.39
根据上述的拟合函数计算误差如下表所示:
由以上误差数据可知,数据拟合误差小于1%,完全满足超声波水表对于精度的要求。
本实施例中的其他部分均为现有技术,在此不再赘述。
Claims (2)
1.一种基于瞬时流量的数据搭建超声波流量模型的方法,其特征是:包括如下步骤:
(1)产生稳定瞬时流量:利用瞬时流量产生装置产生需要的流量V;
(2)时间差采样:在稳定的瞬时流量下,超声波流量计量装置采集一定数量的时间差数据t1....tn,所述时间差是指超声波在两点之间的逆流方向的传播时间和顺流方向的传播时间之差;
(3)数据处理:对采集到的时间差数据进行数字滤波算法得到数据t;数据处理时采用算术平均值滤波、加权平均值滤波、滑动平均滤波算法;
(4)数据模型:根据上述的瞬时流量V和时间差数据t通过公式A=V/t计算获得瞬时流量和时间差的流量系数A;
(5)大样本数据模型:改变瞬时流量,重复步骤(1)、(2)、(3)、(4),从而得到多个流量系数A,通过数据拟合方法,得到时间差t与瞬时流量V之间的关系:V=A(t)。
2.根据权利要求1所述的基于瞬时流量的数据搭建超声波流量模型的方法,其特征是:步骤(1)中采集的流量V是在瞬时流量产生装置达到稳流状态之后进行采集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710540868.4A CN107270983B (zh) | 2017-07-05 | 2017-07-05 | 基于瞬时流量的数据搭建超声波流量模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710540868.4A CN107270983B (zh) | 2017-07-05 | 2017-07-05 | 基于瞬时流量的数据搭建超声波流量模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107270983A CN107270983A (zh) | 2017-10-20 |
CN107270983B true CN107270983B (zh) | 2020-04-07 |
Family
ID=60069933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710540868.4A Active CN107270983B (zh) | 2017-07-05 | 2017-07-05 | 基于瞬时流量的数据搭建超声波流量模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107270983B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110579251A (zh) * | 2019-09-17 | 2019-12-17 | 浙江清环智慧科技有限公司 | 流量测量方法、装置、***、电子设备及可读存储介质 |
CN112665666B (zh) * | 2020-12-18 | 2023-06-30 | 苏州东剑智能科技有限公司 | 一种流体计量表的计量方法 |
CN114235111B (zh) * | 2022-02-24 | 2022-07-15 | 青岛鼎信通讯股份有限公司 | 一种基于模型优化的超声波水表流量校准方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101886939A (zh) * | 2010-06-10 | 2010-11-17 | 宁波大学 | 一种时差法超声流量计静态漂移抑制模型及抑制方法 |
CN102435241A (zh) * | 2011-09-28 | 2012-05-02 | 中国计量科学研究院 | 基于全站仪的超声流量计换能器安装定位方法 |
CN105841762A (zh) * | 2016-03-17 | 2016-08-10 | 广州周立功单片机科技有限公司 | 超声波水表的流量计量方法和*** |
-
2017
- 2017-07-05 CN CN201710540868.4A patent/CN107270983B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101886939A (zh) * | 2010-06-10 | 2010-11-17 | 宁波大学 | 一种时差法超声流量计静态漂移抑制模型及抑制方法 |
CN102435241A (zh) * | 2011-09-28 | 2012-05-02 | 中国计量科学研究院 | 基于全站仪的超声流量计换能器安装定位方法 |
CN105841762A (zh) * | 2016-03-17 | 2016-08-10 | 广州周立功单片机科技有限公司 | 超声波水表的流量计量方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN107270983A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109977345B (zh) | 一种轴流泵叶顶间隙泄漏涡空化的数值模拟方法 | |
CN107270983B (zh) | 基于瞬时流量的数据搭建超声波流量模型的方法 | |
CN103353908B (zh) | 一种基于数值计算的管路阻力系数精确计算方法 | |
CN105758474B (zh) | 一种提高气体超声流量计测量精度的方法 | |
CN107014449B (zh) | 修正泵站流量测量结果的方法 | |
CN113483863B (zh) | 一种确定超声波燃气流量计微调系数的方法及*** | |
CN102087130A (zh) | 基于cfd技术的多声路超声流量计弯管安装声路的优化方法 | |
CN111400941A (zh) | 一种叶片泵内部回流及回流旋涡空化的数值预测方法 | |
CN107014451A (zh) | 基于广义回归神经网络推测超声波流量传感器系数的方法 | |
CN100348863C (zh) | 利用管路特性曲线测量泵特性曲线的方法 | |
CN114397475A (zh) | 一种适用于超声波水表的水流流速测量方法 | |
CN108763785B (zh) | 一种水轮机调节***实时仿真方法、装置和电子设备 | |
CN114091368A (zh) | 一种轴流式水轮机空化状态的识别方法 | |
CN103913203A (zh) | 超声波水表流量系数处理方法 | |
CN111259606B (zh) | 一种泵站引渠及前池含沙动床面流速分布确定方法和装置 | |
CN112304376B (zh) | 基于数据融合的超声波流量计流量测量方法 | |
CN106844964B (zh) | 一种压力管道输水***非恒定流模型的优化方法 | |
CN113514117B (zh) | 一种超声水表流道的设计方法 | |
Adamkowski et al. | The comparative analysis of the current-meter method and the pressure-time method used for discharge measurements in the Kaplan turbine penstocks | |
CN113408180B (zh) | 一种火箭发射场用涡街流量计流固耦合数值分析方法 | |
CN105181049B (zh) | 管道内流体流量的测量方法、测量***及辅助测量*** | |
CN111125867B (zh) | 基于混沌粒子群的化工生产管道实时瞬态模型的建立及计算方法 | |
CN114595596A (zh) | 一种水库淤积物密实沉降速率的模拟方法 | |
CN110388972B (zh) | 非等径超声波流量计串联检定装置及方法 | |
Krzemianowski et al. | Use of non-uniform rational b-splines for discharge calculation in the velocity area method |
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 |