CN114235111B - 一种基于模型优化的超声波水表流量校准方法 - Google Patents
一种基于模型优化的超声波水表流量校准方法 Download PDFInfo
- Publication number
- CN114235111B CN114235111B CN202210168812.1A CN202210168812A CN114235111B CN 114235111 B CN114235111 B CN 114235111B CN 202210168812 A CN202210168812 A CN 202210168812A CN 114235111 B CN114235111 B CN 114235111B
- Authority
- CN
- China
- Prior art keywords
- population
- individuals
- mathematical model
- water meter
- test
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明涉及智慧水务技术领域,公开了一种基于模型优化的超声波水表流量校准方法,包括以下步骤:s1,控制流量,获得标准组和测试组数据;s2,基于标准组数据对流量特征曲线作非线性拟合,计算初始数学模型;s3,基于测试组数据评价数学模型;s4,基于测试组数据优化模型参数;s5,基于数学模型,通过差分进化算法,获得单只水表最优的分段校准系数;s6,重复s3至s5,完成所有水表的流量校准。本发明直接使用原始数据求解校准系数,简单高效,避免因人工选择分段点不合理引起精度误差,从源头修正校准系数,为每支水表优化独立的模型,消除加工或安装误差带来的影响,将精度控制在要求范围内,实现同批次或型号水表的高度一致性。
Description
技术领域
本发明涉及智慧水务技术领域,尤其涉及一种基于模型优化的超声波水表流量校准方法。
背景技术
超声波水表因计量精度高、量程比宽、压损小等优势,在民用、工业等领域得以广泛应用,其通过计算上下游超声波信号的传播时间差计算流速,基于校准系数,获得管道内平均流速。当前,在流量校准过程中,一般参考流体力学的理论公式,依据人工经验,进行流量分段校准,过程反复迭代,且很难保证最优解。此外,考虑时间成本和经济成本,在对同型号多支水表的校准系数进行二次修正时,利用的测量数据往往较少,通过传统二次修正方法很难保证水表在全流量区间内的精度。
发明内容
本发明针对现有技术存在的不足和缺陷,提供了一种基于模型优化的超声波水表流量校准方法。该方法直接使用原始测量数据求解校准系数,避免了因人工选择分段点不合理而引起的精度误差。此外,从源头修正校准系数,在初始数学模型基础上,为每支水表优化独立的数学模型参数,消除因加工或安装误差带来的影响,将精度控制在要求范围内。
本发明的目的可以通过以下技术方案来实现:
一种基于模型优化的超声波水表流量校准方法,其特征在于,包括以下步骤:
s1:使用一只已校准的超声波水表在全流量区间内测量到的流量数据作为标准组,以所有待校准的超声波水表在全流量区间内测量到的流量数据作为测试组;
s2:基于标准组数据,绘制流量特征曲线,对流量特征曲线进行非线性拟合,获得初始数学模型;
s3:基于测试组数据计算相对误差RD,根据相对误差RD评价数学模型,以评价结果决定是否跳过步骤s4;
s4:基于测试组数据,通过差分进化算法,对数学模型参数进行优化;
s5:基于数学模型,通过差分进化算法,获得单只待校准超声波水表最优的分段校准系数;
s6:重复步骤s3至步骤s5,完成所有待校准超声波水表的流量校准。
进一步地,所述上下游接收信号传播时间差是以上下游接收信号作为输入,通过时差法计算得到的。
进一步地,所述步骤s2中流量特征曲线为校准系数随上下游接收信号传播时间差的平均值变化的曲线,其中校准系数K计算公式如下:
式中,为在第j次测量时第i组的上下游接收信号传播时间差,Nj为第j次测量时的上下游接收信号传播时间差总数,为第j次测量的已校准的超声波水表获得的标准组的平均流速,为第j次测量的上下游接收信号传播时间差的平均值。
进一步地,所述步骤s3中相对误差RD的计算公式如下:
数学模型评价原则如下:
若RD l 均小于或等于最大允许误差的二分之一,则认为数学模型与测试组数据吻合,跳过步骤s4;
若存在RD l 大于最大允许误差的二分之一,且小于或等于最大允许误差,则认为数学模型与测试组数据不吻合,需进入步骤s4进行数学模型参数的优化;
若存在RD l 大于最大允许误差,则认为数学模型不满足要求,需脱离本流程,对该超声波水表单独分析。
进一步地,所述步骤s4具体步骤如下:
s4.2:设定差分进化算法对数学模型的搜索范围,初始化种群个体数目、最大迭代次数、交叉因子和变异因子;
s4.3:在数学模型参数θ的搜索范围内随机生成种群作为父代种群,迭代次数加1,根据目标函数和约束条件,计算父代种群中个体的适应度;
s4.4:基于交叉因子与变异因子,在父代种群中随机选择个体进行交叉和变异,生成试验种群;
s4.5:计算试验种群中个体的适应度,将父代种群与试验种群中个体的适应度进行比较,根据比较结果生成子代种群;
s4.6:进行子代种群中个体的非支配排序;
s4.7:计算子代种群中个体的拥挤距离,剔除拥挤距离相对较小的个体以保持子代种群个体数目与初始种群个体数目一致;
s4.8:判断现有迭代次数是否到达最大迭代次数要求,若满足,则完成数学模型参数优化,结束流程;否则跳转至步骤s4.3。
进一步地,所述步骤s5具体步骤如下:
s5.2:设定差分进化算法对参数γ的搜索范围,初始化种群个体数目、最大迭代次数、交叉因子和变异因子;
s5.3:在参数γ的搜索范围内随机生成种群作为父代种群,迭代次数加1,根据目标函数和约束条件,计算父代种群中个体的适应度;
s5.4:基于交叉因子与变异因子,在父代种群中随机选择个体进行交叉和变异,生成试验种群;
s5.5:计算试验种群中个体的适应度,将父代种群与试验种群中个体的适应度进行比较,根据比较结果生成子代种群;
s5.6:进行子代种群中个体的非支配排序;
s5.7:计算子代种群中个体的拥挤距离,剔除拥挤距离相对较小的个体以保持子代种群个体数目与初始种群个体数目一致;
s5.8:判断现有迭代次数是否到达最大迭代次数要求,若满足,则完成参数优化,结束流程;否则跳转至步骤s5.3。
本发明的有益技术效果:直接使用原始测量数据求解校准系数,简单高效,避免了因人工选择分段点不合理而引起的精度误差。此外,从源头修正校准系数,在初始数学模型基础上,为每支水表优化独立的数学模型参数,消除因加工或安装误差带来的影响,将精度控制在要求范围内,从而实现同批次或型号水表的高度一致性。
附图说明
图1为本发明的总体流程图。
图2为本发明实施例中超声波水表的流量特征曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不限定本发明。
实施例:
如图1所示,一种基于模型优化的超声波水表流量校准方法,包括以下步骤:
s1:调节流量控制装置,待流量稳定后,使用一只已校准的超声波水表在全流量区间内测量到的流量数据作为标准组,以所有待校准的超声波水表在全流量区间内测量到的流量数据作为测试组。数据集包括每支超声波水表测得的平均流速及上下游接收信号传播时间差的平均值,其中平均流速的计算公式如下:
所述上下游接收信号传播时间差是以上下游接收信号作为输入,通过时差法计算得到的。
s2:基于标准组数据,绘制流量特征曲线,对流量特征曲线进行非线性拟合,获得初始数学模型。流量特征曲线为校准系数随上下游接收信号传播时间差的平均值的变化曲线,其中校准系数K计算公式如下:
式中,为在第j次测量时第i组的上下游接收信号传播时间差,Nj为第j次测量时的上下游接收信号传播时间差总数,为第j次测量的已校准的超声波水表获得的标准组的平均流速,为第j次测量的上下游接收信号传播时间差的平均值。
如图2所示为某型号DN15超声波水表的流量特征曲线,校准系数随上下游接收信号传播时间差的平均值的变化而改变。基于图2数据进行最小二乘拟合后,获得初始数学模型的表达式如下:
式中,a、b、c、d均为数学模型中的可调参数,取值与超声波水表型号相关。参数 a的取值范围为[0.02016, 0.02107],参数b的取值范围为[-0.0001136, 0.0003671],参数c的取值范围为[-0.009227, -0.008075],参数d的取值范围为[-0.2465, -0.1726]。
s3:基于测试组数据计算相对误差RD,根据相对误差RD评价数学模型,以评价结果决定是否跳过步骤s4。相对误差RD的计算公式如下:
数学模型评价原则如下:
若RD l 均小于或等于最大允许误差的二分之一,则认为数学模型与测试组数据吻合,跳过步骤s4;
若存在RD l 大于最大允许误差的二分之一,且小于或等于最大允许误差,则认为数学模型与测试组数据不吻合,需进入步骤s4进行数学模型参数的优化;
若存在RD l 大于最大允许误差,则认为数学模型不满足要求,需脱离本流程,对该超声波水表单独分析。
实施例中,某型号DN15超声波水表为2级表,设定最大允许误差即为国标要求的2级精度。
s4:基于测试组数据,通过差分进化算法,对数学模型参数进行优化。具体步骤如下:
s4.2:设定差分进化算法对数学模型的搜索范围,初始化种群个体数目、最大迭代次数、交叉因子和变异因子;
s4.3:在数学模型参数θ的搜索范围内随机生成种群作为父代种群,迭代次数加1,根据目标函数和约束条件,计算父代种群中个体的适应度;
s4.4:基于交叉因子与变异因子,在父代种群中随机选择个体进行交叉和变异,生成试验种群;
s4.5:计算试验种群中个体的适应度,将父代种群与试验种群中个体的适应度进行比较,根据比较结果生成子代种群;
s4.6:进行子代种群中个体的非支配排序;
s4.7:计算子代种群中个体的拥挤距离,剔除拥挤距离相对较小的个体以保持子代种群个体数目与初始种群个体数目一致;
s4.8:判断现有迭代次数是否到达最大迭代次数要求,若满足,则完成数学模型参数优化,结束流程;否则跳转至步骤s4.3。
s5.2:设定差分进化算法对参数γ的搜索范围,初始化种群个体数目、最大迭代次数、交叉因子和变异因子;
s5.3:在参数γ的搜索范围内随机生成种群作为父代种群,迭代次数加1,根据目标函数和约束条件,计算父代种群中个体的适应度;
s5.4:基于交叉因子与变异因子,在父代种群中随机选择个体进行交叉和变异,生成试验种群;
s5.5:计算试验种群中个体的适应度,将父代种群与试验种群中个体的适应度进行比较,根据比较结果生成子代种群;
s5.6:进行子代种群中个体的非支配排序;
s5.7:计算子代种群中个体的拥挤距离,剔除拥挤距离相对较小的个体以保持子代种群个体数目与初始种群个体数目一致;
s5.8:判断现有迭代次数是否到达最大迭代次数要求,若满足,则完成参数优化,结束流程;否则跳转至步骤s5.3。
s6:重复步骤s3至步骤s5,完成所有待校准超声波水表的流量校准。
上述实施例是对本发明的具体实施方式的说明,而非对本发明的限制,有关技术领域的技术人员在不脱离本发明的精神和范围的情况下,还可做出各种变换和变化以得到相对应的等同的技术方案,因此所有等同的技术方案均应归入本发明的专利保护范围。
Claims (3)
1.一种基于模型优化的超声波水表流量校准方法,其特征在于,包括以下步骤:
s1:使用一只已校准的超声波水表在全流量区间内测量到的流量数据作为标准组,以所有待校准的超声波水表在全流量区间内测量到的流量数据作为测试组;最终得到的数据集包括每支超声波水表测得的平均流速及上下游接收信号传播时间差的平均值,其中平均流速的计算公式如下:
s2:基于标准组数据,绘制流量特征曲线,对流量特征曲线进行非线性拟合,获得初始数学模型;其中流量特征曲线为校准系数随上下游接收信号传播时间差的平均值变化的曲线,其中校准系数K计算公式如下:
式中,为在第j次测量时第i组的上下游接收信号传播时间差,Nj为第j次测量时的上下游接收信号传播时间差总数,为第j次测量的已校准的超声波水表获得的标准组的平均流速,为第j次测量的上下游接收信号传播时间差的平均值;
s3:基于测试组数据计算相对误差RD,根据相对误差RD评价数学模型,以评价结果决定是否跳过步骤s4;其中相对误差RD的计算公式如下:
数学模型评价原则如下:
若RD l 均小于或等于最大允许误差的二分之一,则认为数学模型与测试组数据吻合,跳过步骤s4;
若存在RD l 大于最大允许误差的二分之一,且小于或等于最大允许误差,则认为数学模型与测试组数据不吻合,需进入步骤s4进行数学模型参数的优化;
若存在RD l 大于最大允许误差,则认为数学模型不满足要求,需脱离本流程,对该超声波水表单独分析;
s4:基于测试组数据,通过差分进化算法,对数学模型参数进行优化;具体步骤如下:
s4.2:设定差分进化算法对数学模型的搜索范围,初始化种群个体数目、最大迭代次数、交叉因子和变异因子;
s4.3:在数学模型参数θ的搜索范围内随机生成种群作为父代种群,迭代次数加1,根据目标函数和约束条件,计算父代种群中个体的适应度;
s4.4:基于交叉因子与变异因子,在父代种群中随机选择个体进行交叉和变异,生成试验种群;
s4.5:计算试验种群中个体的适应度,将父代种群与试验种群中个体的适应度进行比较,根据比较结果生成子代种群;
s4.6:进行子代种群中个体的非支配排序;
s4.7:计算子代种群中个体的拥挤距离,剔除拥挤距离相对较小的个体以保持子代种群个体数目与初始种群个体数目一致;
s4.8:判断现有迭代次数是否到达最大迭代次数要求,若满足,则完成数学模型参数优化,结束流程;否则跳转至步骤s4.3;
s5:基于数学模型,通过差分进化算法,获得单只待校准超声波水表最优的分段校准系数;具体步骤如下:
s5.2:设定差分进化算法对参数γ的搜索范围,初始化种群个体数目、最大迭代次数、交叉因子和变异因子;
s5.3:在参数γ的搜索范围内随机生成种群作为父代种群,迭代次数加1,根据目标函数和约束条件,计算父代种群中个体的适应度;
s5.4:基于交叉因子与变异因子,在父代种群中随机选择个体进行交叉和变异,生成试验种群;
s5.5:计算试验种群中个体的适应度,将父代种群与试验种群中个体的适应度进行比较,根据比较结果生成子代种群;
s5.6:进行子代种群中个体的非支配排序;
s5.7:计算子代种群中个体的拥挤距离,剔除拥挤距离相对较小的个体以保持子代种群个体数目与初始种群个体数目一致;
s5.8:判断现有迭代次数是否到达最大迭代次数要求,若满足,则完成参数优化,结束流程;否则跳转至步骤s5.3;
s6:重复步骤s3至步骤s5,完成所有待校准超声波水表的流量校准。
2.根据权利要求1所述的一种基于模型优化的超声波水表流量校准方法,其特征在于,所述上下游接收信号传播时间差是以上下游接收信号作为输入,通过时差法计算得到的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210168812.1A CN114235111B (zh) | 2022-02-24 | 2022-02-24 | 一种基于模型优化的超声波水表流量校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210168812.1A CN114235111B (zh) | 2022-02-24 | 2022-02-24 | 一种基于模型优化的超声波水表流量校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114235111A CN114235111A (zh) | 2022-03-25 |
CN114235111B true CN114235111B (zh) | 2022-07-15 |
Family
ID=80747867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210168812.1A Active CN114235111B (zh) | 2022-02-24 | 2022-02-24 | 一种基于模型优化的超声波水表流量校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114235111B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114485865B (zh) * | 2022-04-15 | 2022-09-27 | 青岛鼎信通讯股份有限公司 | 一种基于香农熵的超声波水表流量校准方法 |
CN115950495B (zh) * | 2023-02-09 | 2023-06-09 | 青岛鼎信通讯科技有限公司 | 一种超声水表测量频率调节方法 |
CN116793462B (zh) * | 2023-08-21 | 2023-11-03 | 成都千嘉科技股份有限公司 | 基于超声波燃气表的标定方法和标定装置 |
CN117109709B (zh) * | 2023-10-20 | 2024-06-28 | 青岛鼎信通讯科技有限公司 | 一种超声水表校准方法 |
CN117268512B (zh) * | 2023-11-23 | 2024-02-09 | 青岛鼎信通讯科技有限公司 | 一种适用于超声水表的一致性优化方法 |
CN117367527B (zh) * | 2023-12-08 | 2024-03-19 | 青岛鼎信通讯科技有限公司 | 一种能提高超声水表可靠性的计量方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101048645A (zh) * | 2004-08-13 | 2007-10-03 | 恩特格里公司 | 用于流量设备的校准的***和方法 |
CN102435265A (zh) * | 2011-09-19 | 2012-05-02 | 苏州聚元微电子有限公司 | 一种时差法超声波流量计精度修正的方法 |
CN105181040A (zh) * | 2015-09-25 | 2015-12-23 | 辽宁聚焦科技有限公司 | 一种差压式流量计的数字化标定及优化方法 |
CN107270983A (zh) * | 2017-07-05 | 2017-10-20 | 青岛积成电子股份有限公司 | 基于瞬时流量的数据搭建超声波流量模型的方法 |
CN109724657A (zh) * | 2018-12-29 | 2019-05-07 | 杭州莱宸科技有限公司 | 基于改进的Elman神经网络的水表流量计量方法和*** |
CN113033093A (zh) * | 2021-03-24 | 2021-06-25 | 南京航空航天大学 | 一种基于仿真模型的***设计参数多目标优化方法 |
-
2022
- 2022-02-24 CN CN202210168812.1A patent/CN114235111B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101048645A (zh) * | 2004-08-13 | 2007-10-03 | 恩特格里公司 | 用于流量设备的校准的***和方法 |
CN102435265A (zh) * | 2011-09-19 | 2012-05-02 | 苏州聚元微电子有限公司 | 一种时差法超声波流量计精度修正的方法 |
CN105181040A (zh) * | 2015-09-25 | 2015-12-23 | 辽宁聚焦科技有限公司 | 一种差压式流量计的数字化标定及优化方法 |
CN107270983A (zh) * | 2017-07-05 | 2017-10-20 | 青岛积成电子股份有限公司 | 基于瞬时流量的数据搭建超声波流量模型的方法 |
CN109724657A (zh) * | 2018-12-29 | 2019-05-07 | 杭州莱宸科技有限公司 | 基于改进的Elman神经网络的水表流量计量方法和*** |
CN113033093A (zh) * | 2021-03-24 | 2021-06-25 | 南京航空航天大学 | 一种基于仿真模型的***设计参数多目标优化方法 |
Non-Patent Citations (4)
Title |
---|
A Novel Mathematical Model for Transit-time Ultrasonic Flow Measurement;Lei Kang;《2019 IEEE International Ultrasonics Symposium (IUS)》;20191209;全文 * |
Neural-network-based error reduction in calibrating utility ultrasonic flow meters;B. Yazdanshenasshad;《Flow Measurement and Instrumentation》;20181011;全文 * |
Ultrasonic flow detection and research on regression analysis calibration method;Ming-liang Liu;《2012 International Conference on Systems and Informatics (ICSAI2012)》;20120625;全文 * |
基于回波峰值拟合的气体超声流量计信号处理;田雷;《电子测量与仪器学报》;20170731;第31卷(第7期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114235111A (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114235111B (zh) | 一种基于模型优化的超声波水表流量校准方法 | |
CN109001136B (zh) | 一种基于紫外—可见光吸收光谱的cod在线监测方法 | |
CN111291937A (zh) | 基于支持向量分类与gru神经网络联合的处理污水水质预测方法 | |
CN112926265A (zh) | 基于遗传算法优化神经网络的大气多孔探针测量校准方法 | |
CN101473196A (zh) | 测量装置及其使用方法 | |
CN110346005B (zh) | 基于深度学习的科里奥利质量流量计数字信号处理方法 | |
CN103177289B (zh) | 一种噪声不确定复杂非线性动态***的建模方法 | |
CN114964429A (zh) | 一种超声波水表检定方法 | |
CN108362911A (zh) | 一种基于最小二乘法和插值法的流速仪标定方法 | |
CN114487976B (zh) | Mcm电子式互感器校验仪溯源不确定度评定方法及*** | |
CN112818495A (zh) | 一种新型管道压降测算算法参数的动态修正方法 | |
CN116147724B (zh) | 一种适用于超声水表的计量方法 | |
CN112925984A (zh) | 一种基于gcn推荐的样本密度聚合方法 | |
CN114485865B (zh) | 一种基于香农熵的超声波水表流量校准方法 | |
CN111351862A (zh) | 一种超声测量较准方法及测厚方法 | |
CN114485863A (zh) | 超声波水表的流量误差校正方法、***、计算机及介质 | |
CN109101759A (zh) | 一种基于正逆响应面法的参数识别方法 | |
CN111914487B (zh) | 一种基于对抗神经网络的无资料地区水文参数率定方法 | |
CN115034140A (zh) | 一种基于关键控制因子的地表水水质变化趋势预测方法 | |
CN114924489A (zh) | 一种适用于流程工业预测控制的模型自主学习方法 | |
Yi et al. | An improved intelligent calibration method for vortex flowmeter | |
CN113945729B (zh) | 一种基于渠道垂直切面的平均流速计算方法 | |
CN111199774B (zh) | 一种三坐标设备位置校准方法 | |
CN110705186A (zh) | 通过rbf粒子群优化算法即时在线仪表校验和诊断方法 | |
CN116304968B (zh) | 一种基于bp神经网络的超声水表流量数据融合方法和装置 |
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 |