CN117109709A - 一种超声水表校准方法 - Google Patents
一种超声水表校准方法 Download PDFInfo
- Publication number
- CN117109709A CN117109709A CN202311361056.5A CN202311361056A CN117109709A CN 117109709 A CN117109709 A CN 117109709A CN 202311361056 A CN202311361056 A CN 202311361056A CN 117109709 A CN117109709 A CN 117109709A
- Authority
- CN
- China
- Prior art keywords
- flow
- water meter
- ultrasonic water
- point
- fitting
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000010586 diagram Methods 0.000 claims abstract description 12
- 230000000694 effects Effects 0.000 claims abstract description 11
- 238000007781 pre-processing Methods 0.000 claims abstract description 8
- 230000011218 segmentation Effects 0.000 claims description 12
- 238000011144 upstream manufacturing Methods 0.000 claims description 9
- 239000012530 fluid Substances 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000011088 calibration curve Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
- G01F25/10—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明涉及流量校准技术领域,公开了一种超声水表校准方法,包括以下步骤:控制流量连续变化,采集过程中的水表数据作为数据集;数据集预处理;绘制流量特征散点图;利用最小二乘法进行多项式拟合,获得流量特征曲线;计算瞬时流量;计算相对误差RD;根据RD评价流量特征曲线拟合效果,进而决定后续校准步骤。本发明通过实验平台控制流量不间断变化,快速覆盖整个流量区间,同时高测量频率保证流量变化过程中采集到足够多的瞬时数据,极大的提高了校准效率;根据误差自适应拟合校准曲线,保证校准精度。本方法在保证校准精度的前提下,改善了流量校准耗时耗力的问题,满足工厂快速生产的要求。
Description
技术领域
本发明涉及流量校准技术领域,尤其涉及一种超声水表校准方法。
背景技术
超声水表的校准是实现精准计量的前提。对此,目前的校表方式还是机械式的修正方法,流量测试数遍后修正误差,循环操作直到达到要求为止。此方法浪费时间的同时还增加了生产成本,既不适合超声水表的特性,又不能展示超声水表自身特点。而现有资料中还尚未发现其他可行性强的快速校表方法。
发明内容
本发明针对现有技术存在的不足和缺陷,提供了一种超声水表校准方法,在保证超声水表精度的前提下,实现超声水表的快速校准,满足工厂快速生产的要求。
本发明的目的可以通过以下技术方案来实现:
一种超声水表校准方法,包括以下步骤:
S1,通过实验平台,以不固定频率,控制流量在Q1~Q4间连续变化,采集过程中的超声水表数据作为数据集;
S2,对所述数据集进行预处理;
S3,基于预处理后的数据集绘制流量特征散点图;
S4,将流量特征散点图中待校准超声水表的上下游接收信号传播时间差Δt取对数后的结果作为横坐标x,利用最小二乘法进行多项式拟合,获得流量特征曲线K:
;
;
其中a、b、c是公式拟合系数;
S5,基于Δt与K计算瞬时流量vcalc;
S6,计算相对误差RD:
;
式中,vexp是数据集中超声水表采集的瞬时流量;
根据RD评价流量特征曲线拟合效果,以评价结果决定后续步骤:
情况1,若RD均小于或等于最大允许误差的二分之一,则认为流量特征曲线拟合效果满足要求,输出流量特征曲线,完成超声水表校准;
情况2,若存在RD大于最大允许误差的二分之一,且小于或等于最大允许误差,则
认为流量特征曲线拟合效果不满足要求,需将整个流量区间分为层流、湍流、过渡流三段
后,分段进行最小二乘法拟合,输出最终的流量特征曲线,完成超声水表校准;其中每一段
流量区间拟合公式均为:;
情况3,若存在RD大于最大允许误差,且小于最大允许误差的2倍,则认为流量特征曲线拟合错误,需执行如下流程:
以层流/湍流分界点与流量特征曲线拐点为分段点,将整个流量区间划分为若干子流量区间;
分段进行线性拟合;拟合公式均为:;
根据相对误差RD优化子流量区间分段:对线性拟合后的RD大于最大允许误差的二分之一的各子流量区间,分别在其中间位置增加一个分段点;
基于优化后的分段进行线性拟合,输出最终的流量特征曲线,完成超声水表校准;
其中拟合公式均为:;
情况4,若存在RD大于或等于最大允许误差的2倍,则认为超声水表异常,流程终止,需对该超声水表单独分析。
优选地,所述步骤S1中数据集包含待校准超声水表的上下游接收信号传播时间差以及标准表的瞬时流量;
其中标准表是指能够展示真实流量变化的高精度超声水表。
优选地,所述步骤S2中预处理的方法为:
以数据集点列中第3个点为当前点,计算当前点瞬时流量的特征因子k及前后各2
个点的五点均值;其中k=vexp/Δt;
判断是否大于2(max-min)/3:若是则从数据集中剔除当前点,否则以点列
中第4个点为当前点重复上述步骤;其中max、min分别是数据点列中所有点瞬时流量的特征
因子k中的最大值与最小值;
以此类推,直至以点列中的第n-2个点为当前点重复上述步骤完成,至此预处理完成;其中n为总点数。
优选地,所述步骤S3中流量特征散点图为反映单一温度下特征因子k随待校准超声水表上下游接收信号传播时间差Δt变化的散点图;其中k= vexp /Δt。
优选地,所述步骤S5中瞬时流量vcalc的计算公式如下:
。
优选地,所述步骤S6的情况2中将整个流量区间分为层流、湍流、过渡流三段的具体步骤如下:
确定层流、湍流、过渡流区间对应的雷诺数;
根据雷诺数与流量关系计算分段点处的流量值;其中ρ是流体密
度、Re是雷诺数、L是水力直径、μ是流体粘度系数、S是管段截面积。
本发明的有益技术效果:通过实验平台控制流量不间断变化,快速覆盖整个流量区间,同时高测量频率保证流量变化过程中采集到足够多的瞬时数据,极大的提高了校准效率;根据误差自适应拟合校准曲线,保证校准精度。本方法在保证校准精度的前提下,改善了流量校准耗时耗力的问题,满足工厂快速生产的要求。
附图说明
图1为本发明的总体流程图。
图2为本发明实施例中的流量特征曲线图。
图3为本发明实施例中的流量特征散点图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不限定本发明。
实施例:一种超声水表校准方法,参见图1,包括以下步骤:
S1,通过实验平台,以不固定频率,控制流量从Q1~Q4间连续变化,实验时长300s,其中层流区域流量变化速度是湍流区域的二倍,以每秒64组数据的速度采集超声波标准表的瞬时流量与待校表的瞬时上下游接收信号传播时间差作为数据集;
S2,对所述数据集进行预处理:假定采集到的数据集点列中最前2个点和最后2个
点是准确的,以数据集点列中第3个点为当前点,计算当前点瞬时流量的特征因子k及前后
各2个点的五点均值;其中k=vexp/Δt;
判断是否大于2(max-min)/3:若是则从数据集中剔除当前点,否则以点列
中第4个点为当前点重复上述步骤;其中max、min分别是数据点列中所有点瞬时流量的特征
因子k中的最大值与最小值;
以此类推,直至以点列中的第n-2个点为当前点重复上述步骤完成,至此预处理完成;其中n为总点数。
S3,基于预处理后的数据集绘制流量特征散点图,参见图3;其中流量特征散点图为反映单一温度下特征因子k随待校准超声水表上下游接收信号传播时间差Δt变化的散点图;其中k=vexp/Δt,vexp是超声水表采集的瞬时流量,由于Δt的单位是皮秒,因此图中纵坐标标注了特征因子的数量级为109,图2同理。
S4,之后将流量特征散点图中待校准超声水表的上下游接收信号传播时间差Δt取对数后的结果作为横坐标x,利用最小二乘法进行多项式拟合,获得流量特征曲线K,参见图2:
;
;
其中,a、b、c是公式拟合系数;
S5,基于Δt与K计算瞬时流量vcalc,计算公式如下:
;
式中,vcalc是计算的瞬时流量。
S6,之后计算相对误差RD:
;
式中,vexp是数据集中超声水表采集的瞬时流量;
根据RD评价流量特征曲线拟合效果,以评价结果决定后续步骤:
情况1,若RD均小于或等于最大允许误差的二分之一,则认为流量特征曲线拟合效果满足要求,输出流量特征曲线,完成超声水表校准;
情况2,若存在RD大于最大允许误差的二分之一,且小于或等于最大允许误差,则
认为流量特征曲线拟合效果不满足要求,需将整个流量区间分为层流、湍流、过渡流三段
后,分段进行最小二乘法拟合,输出最终的流量特征曲线,完成超声水表校准;其中每一段
流量区间拟合公式均为:;
其中将整个流量区间分为层流、湍流、过渡流三段,具体步骤如下:
层流与过渡流分界点对应雷诺数为2300,湍流与过渡流分界点对应的雷诺数为4000;
根据雷诺数与流量关系计算分段点处的流量;其中ρ是流体密度、
Re是雷诺数、L是水力直径、μ是流体粘度系数、S是管段截面积。
情况3,若存在RD大于最大允许误差,且小于最大允许误差的2倍,则认为流量特征曲线拟合错误,需执行如下流程:
以层流/湍流分界点与流量特征曲线拐点为分段点,将整个流量区间划分为若干子流量区间;
分段进行线性拟合;拟合公式均为:;
根据相对误差RD优化子流量区间分段:对线性拟合后的RD大于最大允许误差的二分之一的各子流量区间,分别在其中间位置增加一个分段点;
基于优化后的分段进行线性拟合,输出最终的流量特征曲线,完成超声水表校准;
其中拟合公式均为:;
情况4,若存在RD大于或等于最大允许误差的2倍,则认为超声水表异常,流程终止,需对该超声水表单独分析。
实施例中输出了最终的流量特征曲线,完成了超声水表的流量校准,满足国标2级表精度要求。
上述实施例是对本发明的具体实施方式的说明,而非对本发明的限制,有关技术领域的技术人员在不脱离本发明的精神和范围的情况下,还可做出各种变换和变化以得到相对应的等同的技术方案,因此所有等同的技术方案均应归入本发明的专利保护范围。
Claims (6)
1.一种超声水表校准方法,其特征在于,包括以下步骤:
S1,通过实验平台,以不固定频率,控制流量在Q1~Q4间连续变化,采集过程中的超声水表数据作为数据集;
S2,对所述数据集进行预处理;
S3,基于预处理后的数据集绘制流量特征散点图;
S4,将流量特征散点图中待校准超声水表的上下游接收信号传播时间差Δt取对数后的结果作为横坐标x,利用最小二乘法进行多项式拟合,获得流量特征曲线K:
;
;
其中a、b、c是公式拟合系数;
S5,基于Δt与K计算瞬时流量vcalc;
S6,计算相对误差RD:
;
式中,vexp是数据集中超声水表采集的瞬时流量;
根据RD评价流量特征曲线拟合效果,以评价结果决定后续步骤:
情况1,若RD均小于或等于最大允许误差的二分之一,则认为流量特征曲线拟合效果满足要求,输出流量特征曲线,完成超声水表校准;
情况2,若存在RD大于最大允许误差的二分之一,且小于或等于最大允许误差,则认为流量特征曲线拟合效果不满足要求,需将整个流量区间分为层流、湍流、过渡流三段后,分段进行最小二乘法拟合,输出最终的流量特征曲线,完成超声水表校准;其中每一段流量区间拟合公式均为:;
情况3,若存在RD大于最大允许误差,且小于最大允许误差的2倍,则认为流量特征曲线拟合错误,需执行如下流程:
以层流/湍流分界点与流量特征曲线拐点为分段点,将整个流量区间划分为若干子流量区间;
分段进行线性拟合;拟合公式均为:;
根据相对误差RD优化子流量区间分段:对线性拟合后的RD大于最大允许误差的二分之一的各子流量区间,分别在其中间位置增加一个分段点;
基于优化后的分段进行线性拟合,输出最终的流量特征曲线,完成超声水表校准;其中拟合公式均为:;
情况4,若存在RD大于或等于最大允许误差的2倍,则认为超声水表异常,流程终止,需对该超声水表单独分析。
2.根据权利要求1所述的一种超声水表校准方法,其特征在于,所述步骤S1中数据集包含待校准超声水表的上下游接收信号传播时间差以及标准表的瞬时流量;
其中标准表是指能够展示真实流量变化的高精度超声水表。
3.根据权利要求1所述的一种超声水表校准方法,其特征在于,所述步骤S2中预处理的方法为:
以数据集点列中第3个点为当前点,计算当前点瞬时流量的特征因子k及前后各2个点的五点均值;其中k=vexp/Δt;
判断是否大于2(max-min)/3:若是则从数据集中剔除当前点,否则以点列中第4个点为当前点重复上述步骤;其中max、min分别是数据点列中所有点瞬时流量的特征因子k中的最大值与最小值;
以此类推,直至以点列中的第n-2个点为当前点重复上述步骤完成,至此预处理完成;其中n为总点数。
4.根据权利要求1所述的一种超声水表校准方法,其特征在于,所述步骤S3中流量特征散点图为反映单一温度下特征因子k随待校准超声水表上下游接收信号传播时间差Δt变化的散点图;其中k=vexp/Δt。
5.根据权利要求1所述的一种超声水表校准方法,其特征在于,所述步骤S5中瞬时流量vcalc的计算公式如下:
。
6.根据权利要求1所述的一种超声水表校准方法,其特征在于,所述步骤S6的情况2中将整个流量区间分为层流、湍流、过渡流三段的具体步骤如下:
确定层流、湍流、过渡流区间对应的雷诺数;
根据雷诺数与流量关系计算分段点处的流量值;其中ρ是流体密度、Re是雷诺数、L是水力直径、μ是流体粘度系数、S是管段截面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311361056.5A CN117109709B (zh) | 2023-10-20 | 一种超声水表校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311361056.5A CN117109709B (zh) | 2023-10-20 | 一种超声水表校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117109709A true CN117109709A (zh) | 2023-11-24 |
CN117109709B CN117109709B (zh) | 2024-06-28 |
Family
ID=
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117367527A (zh) * | 2023-12-08 | 2024-01-09 | 青岛鼎信通讯科技有限公司 | 一种能提高超声水表可靠性的计量方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011226845A (ja) * | 2010-04-16 | 2011-11-10 | Azden Ltd | 超音波流量計 |
WO2013006090A1 (en) * | 2011-07-07 | 2013-01-10 | Belyaev Boris Mikhailovich | The calibration method, applied in operating conditions, for ultrasonic flow meters used for measuring volume and flow rate of single-phase liquid media |
US20140238148A1 (en) * | 2013-02-27 | 2014-08-28 | Daniel Measurement And Control, Inc. | Ultrasonic flow metering with laminar to turbulent transition flow control |
CN106441520A (zh) * | 2016-10-26 | 2017-02-22 | 齐鲁工业大学 | 一种小口径超声波水表精度校准方法及其校准系统 |
WO2020209806A1 (en) * | 2019-04-10 | 2020-10-15 | Pvte Company Limited | Calibration method for liquid flowmeter |
CN114235111A (zh) * | 2022-02-24 | 2022-03-25 | 青岛鼎信通讯股份有限公司 | 一种基于模型优化的超声波水表流量校准方法 |
CN114459573A (zh) * | 2022-01-29 | 2022-05-10 | 深圳友讯达科技股份有限公司 | 超声波燃气表的新型校准模型、校准方法及超声波燃气表 |
CN114964429A (zh) * | 2022-05-27 | 2022-08-30 | 安徽领水科技有限公司 | 一种超声波水表检定方法 |
CN115950495A (zh) * | 2023-02-09 | 2023-04-11 | 青岛鼎信通讯科技有限公司 | 一种超声水表测量频率调节方法 |
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011226845A (ja) * | 2010-04-16 | 2011-11-10 | Azden Ltd | 超音波流量計 |
WO2013006090A1 (en) * | 2011-07-07 | 2013-01-10 | Belyaev Boris Mikhailovich | The calibration method, applied in operating conditions, for ultrasonic flow meters used for measuring volume and flow rate of single-phase liquid media |
US20140238148A1 (en) * | 2013-02-27 | 2014-08-28 | Daniel Measurement And Control, Inc. | Ultrasonic flow metering with laminar to turbulent transition flow control |
CN106441520A (zh) * | 2016-10-26 | 2017-02-22 | 齐鲁工业大学 | 一种小口径超声波水表精度校准方法及其校准系统 |
WO2020209806A1 (en) * | 2019-04-10 | 2020-10-15 | Pvte Company Limited | Calibration method for liquid flowmeter |
CN114459573A (zh) * | 2022-01-29 | 2022-05-10 | 深圳友讯达科技股份有限公司 | 超声波燃气表的新型校准模型、校准方法及超声波燃气表 |
CN114235111A (zh) * | 2022-02-24 | 2022-03-25 | 青岛鼎信通讯股份有限公司 | 一种基于模型优化的超声波水表流量校准方法 |
CN114964429A (zh) * | 2022-05-27 | 2022-08-30 | 安徽领水科技有限公司 | 一种超声波水表检定方法 |
CN115950495A (zh) * | 2023-02-09 | 2023-04-11 | 青岛鼎信通讯科技有限公司 | 一种超声水表测量频率调节方法 |
Non-Patent Citations (1)
Title |
---|
金福一;: "流速仪检定系统中数据处理方法的研究", 水利信息化, no. 02, 25 April 2015 (2015-04-25), pages 59 - 62 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117367527A (zh) * | 2023-12-08 | 2024-01-09 | 青岛鼎信通讯科技有限公司 | 一种能提高超声水表可靠性的计量方法 |
CN117367527B (zh) * | 2023-12-08 | 2024-03-19 | 青岛鼎信通讯科技有限公司 | 一种能提高超声水表可靠性的计量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Durst et al. | LDA measurements in the near-wall region of a turbulent pipe flow | |
CN103808381B (zh) | 一种时差式超声波流量计的温度影响消除方法 | |
CN109240204B (zh) | 一种基于两步法的数控机床热误差建模方法 | |
CN114235111B (zh) | 一种基于模型优化的超声波水表流量校准方法 | |
CN104390680B (zh) | 时差式超声波热量表流量误差修正方法 | |
CN108119318B (zh) | 基于机组测风系统的风机技改优化效果评估方法及其系统 | |
CN114964429A (zh) | 一种超声波水表检定方法 | |
KR101852686B1 (ko) | 비행체의 유동속도 및 유동각을 도출하는 방법 | |
CN113947035A (zh) | 一种高超声速边界层转捩数据天地相关性方法 | |
CN115358026A (zh) | 一种基于多元线性回归与曲面拟合的五孔探针数据处理方法 | |
CN117109709B (zh) | 一种超声水表校准方法 | |
CN110751213B (zh) | 一种测风塔异常风速数据识别与补齐的方法 | |
CN117109709A (zh) | 一种超声水表校准方法 | |
Ndoye et al. | Parameterizable constant temperature anemometer: a new method for the analysis of velocity–temperature coupling in turbulent heat transfer | |
CN117249922B (zh) | 一种用于温度测试仪的温度校准方法及系统 | |
CN114485865B (zh) | 一种基于香农熵的超声波水表流量校准方法 | |
CN107818240A (zh) | 基于伴随同化的纵向离散系数反演方法 | |
CN219301749U (zh) | 一种机器视觉渠道流量测量设备标定装置 | |
Subasi et al. | Calibration of constant temperature anemometry with hot-film probes for low speed laminar water channel flows | |
CN112685954B (zh) | 一种汽车环境风洞风速和风机转速的预测方法及装置 | |
CN114414098A (zh) | 一种分布式光纤拉曼测温系统及其多点温度校正方法 | |
CN111504428B (zh) | 一种热式气体流量计快速标定方法 | |
CN116090340B (zh) | 一种基于数据分析的热工控制时间延迟估计方法 | |
CN116147724A (zh) | 一种适用于超声水表的计量方法 | |
Raffeiner et al. | Prediction of 2D aerodynamic flow fields on the example of a NACA0012 airfoil-A comparison of experimental data, CFD simulations and CNN predictions |
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 |