CN118243186A - 一种超声波流量计计量方法及系统 - Google Patents
一种超声波流量计计量方法及系统 Download PDFInfo
- Publication number
- CN118243186A CN118243186A CN202410449035.7A CN202410449035A CN118243186A CN 118243186 A CN118243186 A CN 118243186A CN 202410449035 A CN202410449035 A CN 202410449035A CN 118243186 A CN118243186 A CN 118243186A
- Authority
- CN
- China
- Prior art keywords
- flight
- difference
- filtering
- time
- filtering result
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000001914 filtration Methods 0.000 claims abstract description 152
- 239000000523 sample Substances 0.000 claims description 32
- 238000012545 processing Methods 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 10
- 239000012530 fluid Substances 0.000 description 4
- 230000005284 excitation Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F15/00—Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
- G01F15/06—Indicating or recording devices
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明提出了一种超声波流量计计量方法及系统。该方法为:计算各飞行时间样本的飞行时间差;当飞行时间样本数量满足样本数预设值时,将所有飞行时间差进行中位数滤波,得到中位数样本;当中位数样本数量满足中位数样本数预设值时,将所有中位数样本进行平均滤波,得到当前飞行时间差滤波结果;当当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值不在设定范围内时,调整滤波参数,重新对所有飞行时间差进行平均滤波;当两者差值在设定范围内时,基于当前飞行时间差滤波结果进行流量计量。该超声波流量计计量方法通过软件滤波的方式使在水流不稳定的情况下仍能准确地对流速、流量进行计量,提升了计量精度。
Description
技术领域
本发明涉及超声波流量计领域,具体涉及一种超声波流量计计量方法及系统。
背景技术
超声波流量计是通过检测流体流动对超声束或超声脉冲的作用以测量流量的仪表,可广泛应用于医疗、石油、化工、冶金、电力、给排水等领域。
外夹式或者管段式超声波流量计是以“速度差法”为原理来对对流量进行计量。但在计量时往往会存在一些干扰因素,如流体中存在气泡,水流不稳等情况,从而导致对流量的计量不准确。
在现有技术中,如CN116202581A所记载的一种结合FPGA技术和高性能算法的超声波流量计对气泡问题进行了优化,但对水流不稳问题仍是束手无策。目前解决水流不稳的问题通常是采用整流器的方式去解决水流不稳定造成的时间差波动问题,而在一些微小管径下,整流器无法安装,这样就无法测量微小管径水流不稳定的情况。
发明内容
为了克服上述现有技术中存在的缺陷,本发明的目的是提供一种可有效应对水流不稳问题的超声波流量计计量方法及系统。
为了实现本发明的上述目的,本发明提供了一种超声波流量计计量方法,包括以下步骤:
超声波流量计持续采集飞行时间样本,计算各飞行时间样本中上行飞行时间和下行飞行时间的飞行时间差;
当飞行时间样本数量满足样本数预设值时,将所有飞行时间差进行中位数滤波处理,得到中位数样本;
当中位数样本数量满足中位数样本数预设值时,将所有中位数样本进行平均滤波处理,得到当前飞行时间差滤波结果;
将当前飞行时间差滤波结果与上次飞行时间差滤波结果进行比较,当两者差值不在设定范围内时,调整滤波参数,重新对所有飞行时间差进行滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内;当两者差值在设定范围内时,基于当前飞行时间差滤波结果进行流量计量。
该超声波流量计计量方法通过软件滤波的方式使在水流不稳定的情况下仍能准确地对流速、流量进行计量,提升了计量精度。
可选的,所述平均滤波采用滑动平均滤波。
在超声波流量计计量方法的一种可选方案中,当当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值不在设定范围内时,调整平均滤波参数,重新对所有中位数样本进行平均滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内。
该可选方案通过反复调整平均滤波参数,提高当前飞行时间差滤波结果的准确度,提升了计量的精度。
在超声波流量计计量方法的一种可选方案中,若调整平均滤波参数的次数达到设定次数,但当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值仍不在设定范围内时,则调整中位数滤波参数,重新对所有飞行时间差进行滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内。该可选方案进一步提升了计量精度。
在超声波流量计计量方法的一种可选方案中,记录当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值不在设定范围内的持续时间,如果该持续时间小于设定时间,则调整平均滤波参数,重新进行平均滤波处理;否则,基于当前飞行时间差滤波结果进行流量计量。
在该可选方案中,当持续时间小于设定时间,说明流速虽然发生变化,但仅是短时间内的波动不稳,因此调整平均滤波参数,重新进行平均滤波处理,当持续时间不小于设定时间,说明流速发生变化,且在一定时间内趋于稳定了,因此就可以当前飞行时间差滤波结果进行流量计量。因此该可选方案可有效确定流速的变化是否趋于稳定,进一步的提高了计量的准确性。
在该超声波流量计计量方法的一种可选方案中,基于当前飞行时间差滤波结果进行流量计量时,先将当前飞行时间差滤波结果进行卡尔曼滤波处理,再基于该滤波结果进行流量计量。
在该超声波流量计计量方法的一种可选方案中,若当前飞行时间差滤波结果与上次飞行时间差滤波结果满足公式:
Δt2×(1-p)<Δt1<Δt2×(1+p),判定当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内,其中Δt1为当前飞行时间差滤波结果,Δt2为上次飞行时间差滤波结果,p为计算精度值。
本发明还提供了一种超声波流量计计量系统,包括飞行时间样本获取模块、存储模块和处理模块;
所述飞行时间样本获取模块与处理模块连接,将获取的飞行时间样本发送对处理模块,所述处理模块与存储模块连接,所述存储模块中存放有至少一个可执行指令,所述可执行指令使所述处理模块基于所述飞行时间样本执行如上述的超声波流量计计量方法对应的操作。
在该超声波流量计计量系统的一种可选方案中,所述飞行时间样本获取模块包括至少一对用于采集飞行时间样本的收发一体式超声波探头。
本发明的有益效果是:本发明通过中位值滤波去除异常数据,通过滑动平均滤波降低数据的波动,还通过与上次飞行时间差滤波结果的差异性处理来调节滑动平均滤波的参数,滑动平均滤波后再经卡尔曼滤波处理,基于该滤波结果进行流量计量,克服了水流不稳而导致的计量不准的问题,提高了计量精度。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是实施例一的原理流程图;
图2是实施例一的原理框图;
图3是超声波流量计原理框图;
图4是实施例二的原理流程图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
实施例一
如图1和图2所示,本发明提供了一种超声波流量计计量方法。具体如下:
超声波流量计采集飞行时间样本,飞行时间样本包括同一采集时刻的上行飞行时间和下行飞行时间,计算各飞行时间样本的飞行时间差。
具体地,在同一时刻记录超声波流量计采集到的上行飞行时间t1和下行飞行时间t2。如图3所示,由第一收发一体式超声波探头C发送激励信号,由第二收发一体式超声波探头B接收信号得到上行飞行时间t1,由第二收发一体式超声波探头B发送激励信号,第一收发一体式超声波探头C接收信号得到下行飞行时间t2,由下行飞行时间t2减去上行飞行时间t1得到Δt。
其中,l为第二收发一体式超声波探头B到第一收发一体式超声波探头C的距离;C0为声速;v为流体流速。
由上式得,
持续采集飞行时间样本,当采集的飞行时间样本数量N满足样本数预设值n时,n为正整数,至少大于等于2,将所有飞行时间样本的飞行时间差Δt进行中位数滤波处理,得到中位数样本,记录中位数柆本数量,记为N1。
当中位数样本数量N1小于中位数样本数预设值n1时,则继续采集飞行时间样本,直至中位数样本数量N1不小于中位数样本数预设值n1;当中位数样本数量N1不小于中位数样本数预设值n1时,将所有中位数样本进行平均滤波处理,得到当前飞行时间差滤波结果Δt1。本实施例中,平均滤波采用滑动平均滤波,即假设中位数样本数预设值n1为10,则参与平均滤波处理的中位数样本形成一个数组arry,每得到一个中位数样本就将这个值赋值到arry[0],上一个中位数样本则滑动到arry[1],当中位数样本大于10时,先进入数组arry的中位数样本则溢出,不参与平均滤波。
将当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2进行比较,当两者差值不在设定范围内时,调整滤波参数,重新对所有飞行时间差进行滤波处理,直至当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2的差值在设定范围内;当两者差值在设定范围内时,基于当前飞行时间差滤波结果Δt1进行流量计量,并记录当前飞行时间差滤波结果Δt1以作下轮计算时使用。
本实施例中,若当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2满足公式Δt2×(1-p)<Δt1<Δt2×(1+p),则判定当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2的差值在设定范围内,基于当前飞行时间差滤波结果Δt1进行流量计量其中,p为计算精度值。若Δt1<Δt2×(1-p)或Δt1>Δt2×(1+p),则判定当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2的差值不在设定范围内,需调整平均滤波参数,重新对所有中位数样本进行平均滤波处理,直至当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2的差值在设定范围内。
在调整滤波参数时,先调整平均滤波参数,重新对所有中位数样本进行平均滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内。
若调整平均滤波参数的次数达到设定次数,但当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值仍不在设定范围内时,则调整中位数滤波参数加大丢弃飞行时间差Δt的个数,一般成2的次方增长,然后按上述步骤重新对所有飞行时间差进行滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内。在调整了中位数滤波参数后,平均滤波参数优选但不限于恢复到调整前。
在基于当前飞行时间差滤波结果Δt1进行流量计量时,可先将当前飞行时间差滤波结果Δt1进行卡尔曼滤波处理,再基于该滤波结果进行流量计量,输出最终结果。
具体计算时,流体流速单位时间的流量其中,r为流道半径,π为圆周率。
实施例二
本实施例还提供了实施例一的一种可选方案:如图4所示,当当前飞行时间差滤波结果Δt1与上次飞行时间差滤波结果Δt2的差值不在设定范围内时,开始记录该情况的持续时间,如果该持续时间小于设定时间T,则调整平均滤波参数,重新对所有中位数样本进行平均滤波处理,并按上述实施例的步骤执行;如果该持续时间不小于设定时间T,则基于当前飞行时间差滤波结果Δt1进行流量计量。
记录当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值不在设定范围内的持续时间,如果该持续时间大于设定时间,则调整平均滤波参数,重新进行平均滤波处理;否则,基于当前飞行时间差滤波结果进行流量计量。
实施例三
本发明还提出了一种超声波流量计计量系统的实施例。该实施例中的超声波流量计计量系统包括飞行时间样本获取模块、存储模块和处理模块。
飞行时间样本获取模块与处理模块连接,将获取的飞行时间样本发送对处理模块。飞行时间样本获取模块可以是至少一对用于采集飞行时间样本的收发一体式超声波探头,如实施例一中所述的第一收发一体式超声波探头C和第二收发一体式超声波探头B,也可以是其他形式的数据接收模块。
处理模块与存储模块连接,存储模块中存放有至少一个可执行指令,可执行指令使所述处理模块基于飞行时间样本执行如实施例一或实施例二所述的超声波流量计计量方法对应的操作。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。
Claims (9)
1.一种超声波流量计计量方法,其特征在于,包括以下步骤:
超声波流量计持续采集飞行时间样本,计算各飞行时间样本中上行飞行时间和下行飞行时间的飞行时间差;
当飞行时间样本数量满足样本数预设值时,将所有飞行时间差进行中位数滤波处理,得到中位数样本;
当中位数样本数量满足中位数样本数预设值时,将所有中位数样本进行平均滤波处理,得到当前飞行时间差滤波结果;
将当前飞行时间差滤波结果与上次飞行时间差滤波结果进行比较,当两者差值不在设定范围内时,调整滤波参数,重新对所有飞行时间差进行滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内;当两者差值在设定范围内时,基于当前飞行时间差滤波结果进行流量计量。
2.根据权利要求1所述的超声波流量计计量方法,其特征在于,当当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值不在设定范围内时,调整平均滤波参数,重新对所有中位数样本进行平均滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内。
3.根据权利要求2所述的超声波流量计计量方法,其特征在于,若调整平均滤波参数的次数达到设定次数,但当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值仍不在设定范围内时,则调整中位数滤波参数,重新对所有飞行时间差进行滤波处理,直至当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内。
4.根据权利要求1所述的超声波流量计计量方法,其特征在于,记录当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值不在设定范围内的持续时间,如果该持续时间小于设定时间,则调整平均滤波参数,重新进行平均滤波处理;否则,基于当前飞行时间差滤波结果进行流量计量。
5.根据权利要求1所述的超声波流量计计量方法,其特征在于,基于当前飞行时间差滤波结果进行流量计量时,先将当前飞行时间差滤波结果进行卡尔曼滤波处理,再基于该滤波结果进行流量计量。
6.根据权利要求1所述的超声波流量计计量方法,其特征在于,若当前飞行时间差滤波结果与上次飞行时间差滤波结果满足公式:
Δt2×(1-p)<Δt1<Δt2×(1+p),判定当前飞行时间差滤波结果与上次飞行时间差滤波结果的差值在设定范围内,其中Δt1为当前飞行时间差滤波结果,Δt2为上次飞行时间差滤波结果,p为计算精度值。
7.根据权利要求1所述的超声波流量计计量方法,其特征在于,所述平均滤波采用滑动平均滤波。
8.一种超声波流量计计量系统,其特征在于,包括飞行时间样本获取模块、存储模块和处理模块;
所述飞行时间样本获取模块与处理模块连接,将获取的飞行时间样本发送对处理模块,所述处理模块与存储模块连接,所述存储模块中存放有至少一个可执行指令,所述可执行指令使所述处理模块基于所述飞行时间样本执行如权利要求1至7任一项所述的超声波流量计计量方法对应的操作。
9.根据权利要求8所述的超声波流量计计量系统,其特征在于,所述飞行时间样本获取模块包括至少一对用于采集飞行时间样本的收发一体式超声波探头。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410449035.7A CN118243186A (zh) | 2024-04-15 | 2024-04-15 | 一种超声波流量计计量方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410449035.7A CN118243186A (zh) | 2024-04-15 | 2024-04-15 | 一种超声波流量计计量方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN118243186A true CN118243186A (zh) | 2024-06-25 |
Family
ID=91557022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410449035.7A Pending CN118243186A (zh) | 2024-04-15 | 2024-04-15 | 一种超声波流量计计量方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118243186A (zh) |
-
2024
- 2024-04-15 CN CN202410449035.7A patent/CN118243186A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102519542B (zh) | 压力液位传感器 | |
US5827963A (en) | System and method for determining a density of a fluid | |
CN114001804B (zh) | 一种基于时间差法超声波计量装置的校准方法及系统 | |
CN112267996B (zh) | 一种液压泵的流量脉动测试装置 | |
CN113029263A (zh) | 基于样本方差剔除误差的燃气损失流量测定方法及系统 | |
CN115824331A (zh) | 一种适用于超声水表的低功耗测量方法 | |
CN109716074A (zh) | 超声波流量计及使用部分流量测量的方法 | |
CN102538914A (zh) | 一种具有校正功能的电子显示水表 | |
CN115950495A (zh) | 一种超声水表测量频率调节方法 | |
CN112665666B (zh) | 一种流体计量表的计量方法 | |
CN108458759A (zh) | 测量超声流量计信号延迟时间、流量的方法与处理装置 | |
CN114777887A (zh) | 一种超声波水表的零点漂移消除方法 | |
JPH07280603A (ja) | 機械の異常判定方法 | |
CN118243186A (zh) | 一种超声波流量计计量方法及系统 | |
CN112504410B (zh) | 一种校正射流水表流量计算的方法、射流水表以及校正补偿系数的计量方法 | |
CN210894193U (zh) | 一种发电厂水汽系统在线pH表的准确性检验系统 | |
CN106021855B (zh) | 一种反应堆周期计算方法 | |
CN113188615B (zh) | 一种基于动态偏移补偿的电磁流量计系统及测量方法 | |
KR20170035555A (ko) | 초음파 유량계의 신호처리시스템 | |
CN115388965A (zh) | 一种气体超声流量传感器的测量算法 | |
CN107796492B (zh) | 一种超声波燃气表在线式校准方法 | |
CN106124371A (zh) | 一种基于静电法的气固两相流细度测量装置及测量方法 | |
CN111174842B (zh) | 用于电子式计量表的随机数采样方法 | |
CN109298376B (zh) | 一种基于标准电能表组的电能量值传递的方法及系统 | |
CN117367527B (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 |