CN114485865B - 一种基于香农熵的超声波水表流量校准方法 - Google Patents
一种基于香农熵的超声波水表流量校准方法 Download PDFInfo
- Publication number
- CN114485865B CN114485865B CN202210392379.XA CN202210392379A CN114485865B CN 114485865 B CN114485865 B CN 114485865B CN 202210392379 A CN202210392379 A CN 202210392379A CN 114485865 B CN114485865 B CN 114485865B
- Authority
- CN
- China
- Prior art keywords
- flow
- weight
- follows
- segmentation
- updating
- 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
-
- 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
- Measuring Volume Flow (AREA)
Abstract
本发明涉及流量计量技术领域,公开了一种基于香农熵的超声波水表流量校准方法,包括以下步骤:收集水表测量数据作为数据集;确定初始流量分段点,为各流量分段点分配初始权重;构建权重与流量分段点的更新函数;利用梯度下降法进行二元线性拟合,获得分段校准系数;计算香农熵值,若香农熵值不满足条件,则更新权重与流量分段点,跳转上一步骤s4;若香农熵值满足条件则停止更新,利用最终的分段校准系数校准流量。本发明实现了流量分段点的自适应更新,校准过程不完全依赖人工经验,降低人为误差,自适应更新流量分段点;其次,重点关注流速和温度两大重要因素的影响,实现二维的自适应校准,从而提高超声波水表在多温度场景下的计量准确性。
Description
技术领域
本发明涉及流量校准技术领域,尤其涉及一种基于香农熵的超声波水表流量校准方法。
背景技术
超声波水表因计量精度高、量程比宽、压损小等优势,在民用、工业等领域得以较为广泛的应用,其基于校准系数来进行流量校准。当前,在流量校准过程中,流量分段点一般依据人工经验选取,选择过程反复迭代,消耗大量时间和人力。此外,超声波信号在水中传播速度受温度影响明显,若在校准过程中仅考虑流速因素影响,进行一维的数据拟合,则很难保证超声波水表在温度变化明显的工作场景下的计量精度。
发明内容
本发明针对现有技术存在的不足和缺陷,提供了一种基于香农熵的超声波水表流量校准方法,方法关注流速和温度两大重要因素的影响,实现流量分段点的自适应更新。
本发明的目的可以通过以下技术方案来实现:
一种基于香农熵的超声波水表流量校准方法,包括如下步骤:
s1,通过实验平台控制流量和温度,收集超声波水表的测量数据作为数据集;
s2,基于数据集,绘制流量特征曲线,确定初始流量分段点,为各流量分段点分配初始权重;
s3,构建权重与流量分段点的更新函数;
s4,基于测量数据与流量分段点,利用梯度下降法,进行二元线性拟合,获得分段校准系数;
s5,计算此时的流量分段情况对应的香农熵值,若香农熵值不满足预设条件,则基于s3中构建的更新函数更新权重,基于更新后的权重更新流量分段点,跳转步骤s4;若香农熵值满足预设条件,则停止更新,此时得到最终的分段校准系数,利用最终的分段校准系数校准流量。
优选地,所述步骤s1中通过实验平台控制的低区流量的波动应不超过±2.5%,高区应不超过±5%,一次测量的水温波动应不超过5℃。
优选地,所述测量温度由超声波水表的进出口处的温度计测量得出,上下游接收信号时间差由时差法计算得出。
式中,为第j次测量获取的实测流速均值,为第j次测量获取的上下游接收信号传播时间差的平均值,ΔT i,j 为在第j次测量时根据第i组上下游接收信号计算的时间差,N j 为第j次测量上下游接收信号传播时间差的总计算次数。
优选地,所述步骤s2中初始流量分段点包括层流/湍流分界点、流道切换点和流量特征曲线拐点,各流量分段点初始权重均默认等于1。
优选地,所述步骤s3中权重的更新原则如下:
其中,a 0,l,j ~a 3,l,j 为j次测量所属流量分段区间l对应的校准系数,由s4中二元线性拟合所得且初始值为零,T j 为测量温度;
若待更新权重对应流量分段点的前后两个相邻的流量分段区间内的相对误差RD均大于最大允许误差的二分之一,则将上半句所述的待更新权重更新为2;
若待更新权重对应流量分段点的前一个相邻的流量分段区间内相对误差大于最大允许误差的二分之一,后一个相邻流量区间内相对误差RD小于或等于最大允许误差的二分之一,则将上半句所述的待更新权重更新为-1;
若待更新权重对应流量分段点的前一个相邻的流量分段区间内相对误差小于或等于最大允许误差的二分之一,后一个相邻流量区间内相对误差RD大于最大允许误差的二分之一,则将上半句所述的待更新权重更新为1;
若待更新权重对应流量分段点的前后两个相邻的流量分段区间内的相对误差RD均小于或等于最大允许误差的二分之一,则将上半句所述的待更新权重更新为0。
优选地,所述步骤s3中流量分段点的更新原则如下:
若权重值等于-1,则认为该权重值对应的流量分段点需左移;
若权重值等于0,则认为该权重值对应的流量分段点需去除;
若权重值等于1,则认为该权重值对应的流量分段点需右移;
若权重值等于2,则认为该权重值对应的流量分段点需保留。
优选地,所述步骤s4中二元线性拟合公式如下:
式中,n为测量的总次数。
优选地,基于步骤s3中构建的更新函数更新s5中的流量分段点和权重,停止更新的条件为:香农熵值小于最大允许误差的1%;其中香农熵值的计算公式如下:
式中,H为香农熵值,m为分段区间总数,为第l个分段区间内的相对误差均值,x为第l个分段区间内的测量次数索引,RD lx 为第l个分段区间内第x次测量后计算得到的相对误差值,N l 为第l个分段区间内的总测量次数,MPE l 为第l个分段区间内的最大允许误差。
本发明的有益技术效果:实现了流量分段点的自适应更新,节省了时间和人力成本,减小了人工经验选取分段点导致的误差;方法关注流速和温度两大重要因素的影响,实现二维的自适应校准,从而提高了超声波水表在温度变化明显的工作场景下的计量准确性。
附图说明
图1为本发明的总体流程图,图2为本发明实施例中的流量特征曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不限定本发明。
实施例。
如图1所示,一种基于香农熵的超声波水表流量校准方法,包括如下步骤:
s1,通过实验平台控制流量和温度,收集超声波水表的测量数据作为数据集;其中通过实验平台控制的低区流量的波动应不超过±2.5%,高区应不超过±5%,一次测量的水温波动应不超过5℃。数据集包括多组实测平均流速、测量温度及上下游接收信号传播时间差。
其中测量温度由超声波水表的进出口处的温度计测量得出,上下游接收信号时间差由时差法计算得出。
s2,基于数据集,绘制流量特征曲线,确定初始流量分段点,为各流量分段点分配初始权重;其中流量特征曲线为单一温度下流量特征因子K随上下游接收信号传播时间差的平均值变化的曲线,其中流量特征因子K计算公式如下:
式中,ΔT i,j 为在第j次测量时根据第i组上下游接收信号计算的时间差,N j 为第j次测量上下游接收信号传播时间差的总计算次数。
其中初始流量分段点包括层流/湍流分界点、流道切换点和流量特征曲线拐点,各流量分段点初始权重均默认等于1。根据室温20℃下某型号DN15超声波水表的流量特征因子,设置初始流量分段点为5L/h、20L/h、90 L/h、150 L/h、400L/h、900 L/h、2500L/h。
s3,构建权重与流量分段点的更新函数。
其中权重的更新原则如下:
其中,a 0,l,j ~a 3,l,j 为j次测量所属流量分段区间l对应的校准系数,T j 为测量温度;
若待更新权重对应流量分段点的前后两个相邻的流量分段区间内的相对误差RD均大于最大允许误差的二分之一,则将上半句所述的待更新权重更新为2;
若待更新权重对应流量分段点的前一个相邻的流量分段区间内相对误差大于最大允许误差的二分之一,后一个相邻流量区间内相对误差RD小于或等于最大允许误差的二分之一,则将上半句所述的待更新权重更新为-1;
若待更新权重对应流量分段点的前一个相邻的流量分段区间内相对误差小于或等于最大允许误差的二分之一,后一个相邻流量区间内相对误差RD大于最大允许误差的二分之一,则将上半句所述的待更新权重更新为1;
若待更新权重对应流量分段点的前后两个相邻的流量分段区间内的相对误差RD均小于或等于最大允许误差的二分之一,则将上半句所述的待更新权重更新为0。
其中流量分段点的更新原则如下:
若权重值等于-1,则认为该权重值对应的流量分段点需左移;
若权重值等于0,则认为该权重值对应的流量分段点需去除;
若权重值等于1,则认为该权重值对应的流量分段点需右移;
若权重值等于2,则认为该权重值对应的流量分段点需保留。
s4,基于测量数据与流量分段点,利用梯度下降法,进行二元线性拟合,获得分段校准系数;其中二元线性拟合公式如下:
式中,n为测量的总次数。
s5,计算此时的流量分段情况对应的香农熵值,若香农熵值不满足预设条件,则基于s3中构建的更新函数更新权重,基于更新后的权重更新流量分段点,跳转步骤s4;若香农熵值满足预设条件,则停止更新,此时得到最终的分段校准系数,利用最终的分段校准系数校准流量。基于步骤s3中构建的更新函数更新s5中的流量分段点和权重,停止更新的条件为:香农熵值小于最大允许误差的1%;其中香农熵值的计算公式如下:
式中,H为香农熵值,m为分段区间总数,为第l个分段区间内的相对误差均值,x为第l个分段区间内的测量次数索引,RD lx 为第l个分段区间内第x次测量后计算得到的相对误差值,N l 为第l个分段区间内的总测量次数,MPE l 为第l个分段区间内的最大允许误差。
最终,根据上述计算获取的流量分段点为2.8L/h、3.8L/h、5.6L/h、10L/h、54.8L/h、151.7L/h、405.7L/h、979.1 L/h、2223.5L/h。利用分段区间内最终的校准系数计算流量,满足国标2级精度要求。
上述实施例是对本发明的具体实施方式的说明,而非对本发明的限制,有关技术领域的技术人员在不脱离本发明的精神和范围的情况下,还可做出各种变换和变化以得到相对应的等同的技术方案,因此所有等同的技术方案均应归入本发明的专利保护范围。
Claims (5)
1.一种基于香农熵的超声波水表流量校准方法,其特征在于,包括以下步骤:
s1,通过实验平台控制流量和温度,收集超声波水表的测量数据作为数据集;
s2,基于数据集,绘制流量特征曲线,确定初始流量分段点,为各流量分段点分配初始权重;其中流量特征曲线为单一温度下流量特征因子K随上下游接收信号传播时间差的平均值变化的曲线,其中流量特征因子K计算公式如下:
式中,为第j次测量获取的实测流速均值,为第j次测量获取的上下游接收信号传播时间差的平均值,ΔT i,j 为在第j次测量时根据第i组上下游接收信号计算的时间差,N j 为第j次测量上下游接收信号传播时间差的总计算次数;
s3,构建权重与流量分段点的更新函数;
权重的更新原则如下:
其中,a 0,l,j ~a 3,l,j 为j次测量所属流量分段区间l对应的校准系数,由s4中二元线性拟合所得且初始值为零,T j 为测量温度;
若待更新权重对应流量分段点的前后两个相邻的流量分段区间内的相对误差RD均大于最大允许误差的二分之一,则将上半句所述的待更新权重更新为2;
若待更新权重对应流量分段点的前一个相邻的流量分段区间内相对误差大于最大允许误差的二分之一,后一个相邻流量区间内相对误差RD小于或等于最大允许误差的二分之一,则将上半句所述的待更新权重更新为-1;
若待更新权重对应流量分段点的前一个相邻的流量分段区间内相对误差小于或等于最大允许误差的二分之一,后一个相邻流量区间内相对误差RD大于最大允许误差的二分之一,则将上半句所述的待更新权重更新为1;
若待更新权重对应流量分段点的前后两个相邻的流量分段区间内的相对误差RD均小于或等于最大允许误差的二分之一,则将上半句所述的待更新权重更新为0;
流量分段点的更新原则如下:
若权重值等于-1,则认为该权重值对应的流量分段点需左移;
若权重值等于0,则认为该权重值对应的流量分段点需去除;
若权重值等于1,则认为该权重值对应的流量分段点需右移;
若权重值等于2,则认为该权重值对应的流量分段点需保留;
s4,基于测量数据与流量分段点,利用梯度下降法,进行二元线性拟合,获得分段校准系数;二元线性拟合公式如下:
式中,n为测量的总次数
s5,计算此时的流量分段情况对应的香农熵值,若香农熵值不满足预设条件,则基于s3中构建的更新函数更新权重,基于更新后的权重更新流量分段点,跳转步骤s4;若香农熵值满足预设条件,则停止更新,此时得到最终的分段校准系数,利用最终的分段校准系数校准流量;
停止更新的条件为:香农熵值小于最大允许误差的1%;其中香农熵值的计算公式如下:
2.根据权利要求1所述的一种基于香农熵的超声波水表流量校准方法,其特征在于,所述步骤s1中通过实验平台控制的低区流量的波动应不超过±2.5%,高区应不超过±5%,一次测量的水温波动应不超过5℃。
4.根据权利要求3所述的一种基于香农熵的超声波水表流量校准方法,其特征在于,所述测量温度由超声波水表的进出口处的温度计测量得出,上下游接收信号时间差由时差法计算得出。
5.根据权利要求1所述的一种基于香农熵的超声波水表流量校准方法,其特征在于,所述步骤s2中初始流量分段点包括层流/湍流分界点、流道切换点和流量特征曲线拐点,各流量分段点初始权重均默认等于1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210392379.XA CN114485865B (zh) | 2022-04-15 | 2022-04-15 | 一种基于香农熵的超声波水表流量校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210392379.XA CN114485865B (zh) | 2022-04-15 | 2022-04-15 | 一种基于香农熵的超声波水表流量校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114485865A CN114485865A (zh) | 2022-05-13 |
CN114485865B true CN114485865B (zh) | 2022-09-27 |
Family
ID=81488909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210392379.XA Active CN114485865B (zh) | 2022-04-15 | 2022-04-15 | 一种基于香农熵的超声波水表流量校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114485865B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115773793B (zh) * | 2022-11-29 | 2023-07-21 | 青岛鼎信通讯科技有限公司 | 一种超声水表信号幅度动态调节方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114166318A (zh) * | 2022-02-14 | 2022-03-11 | 青岛鼎信通讯股份有限公司 | 一种基于深度学习的超声波水表流量数据校准方法 |
CN114166299A (zh) * | 2020-09-10 | 2022-03-11 | 克洛纳有限公司 | 用于运行超声流量测量仪的方法和超声流量测量仪 |
CN114235111A (zh) * | 2022-02-24 | 2022-03-25 | 青岛鼎信通讯股份有限公司 | 一种基于模型优化的超声波水表流量校准方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1825150A4 (en) * | 2004-12-02 | 2012-02-01 | Siemens Industry Inc | SYSTEM AND METHOD FOR CORRECTION OF PROFILE CALIBRATION FOR ULTRASONIC FLOW METERS |
-
2022
- 2022-04-15 CN CN202210392379.XA patent/CN114485865B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114166299A (zh) * | 2020-09-10 | 2022-03-11 | 克洛纳有限公司 | 用于运行超声流量测量仪的方法和超声流量测量仪 |
CN114166318A (zh) * | 2022-02-14 | 2022-03-11 | 青岛鼎信通讯股份有限公司 | 一种基于深度学习的超声波水表流量数据校准方法 |
CN114235111A (zh) * | 2022-02-24 | 2022-03-25 | 青岛鼎信通讯股份有限公司 | 一种基于模型优化的超声波水表流量校准方法 |
Non-Patent Citations (1)
Title |
---|
多通道超声波气体流量计的数学模型;鲍丽莎;《声学技术》;20100831;第29卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114485865A (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Provenzano et al. | Experimental analysis of local pressure losses for microirrigation laterals | |
CN114235111B (zh) | 一种基于模型优化的超声波水表流量校准方法 | |
Provenzano et al. | Simplified procedure to evaluate head losses in drip irrigation laterals | |
RU2478919C2 (ru) | Способ калибровки в условиях эксплуатации ультразвуковых расходомеров-счетчиков расхода и объема жидких однофазных сред | |
CN114485865B (zh) | 一种基于香农熵的超声波水表流量校准方法 | |
CN104390680B (zh) | 时差式超声波热量表流量误差修正方法 | |
CN106033000A (zh) | 一种利用雷达波测流仪推算流量的方法 | |
CN113483863B (zh) | 一种确定超声波燃气流量计微调系数的方法及系统 | |
CN107014449B (zh) | 修正泵站流量测量结果的方法 | |
CN112818495A (zh) | 一种新型管道压降测算算法参数的动态修正方法 | |
Yazdanshenasshad et al. | Neural-network-based error reduction in calibrating utility ultrasonic flow meters | |
DK56593D0 (da) | Stroemningsmaaler | |
CN117109709B (zh) | 一种超声水表校准方法 | |
CN116822180B (zh) | 一种模拟河渠输水过程的积分时滞模型参数辨识方法 | |
CN113268855B (zh) | 一种环形交叉口微观交通仿真模型的标定方法 | |
CN206223265U (zh) | 质量流量计的质量标定系统 | |
CN107270983A (zh) | 基于瞬时流量的数据搭建超声波流量模型的方法 | |
CN201637445U (zh) | 基于热线阵列和智能二维重建算法的流速场测量系统 | |
US20240077346A1 (en) | On-site calibration method of ultrasonic flowmeter based on acoustic velocity detection | |
CN114877975A (zh) | 一种基于管路入口雷诺数的超声水表流量误差修正方法 | |
CN114136387B (zh) | 基于支持向量机模型svm算法的多声道超声流量计误差补偿方法 | |
CN106844964B (zh) | 一种压力管道输水系统非恒定流模型的优化方法 | |
CN111932023B (zh) | 一种基于典型设计洪水过程线的小流域短期洪水预报方法 | |
RU2296959C1 (ru) | Способ градуировки объемных расходомеров теплосчетчика и устройство для его осуществления | |
Cockrell et al. | Production of shear profiles in a wind tunnel by cylindrical rods placed normal to the stream |
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 |