CN105222839A - 超声波流量计仪表系数的非线性修正方法 - Google Patents

超声波流量计仪表系数的非线性修正方法 Download PDF

Info

Publication number
CN105222839A
CN105222839A CN201510519962.2A CN201510519962A CN105222839A CN 105222839 A CN105222839 A CN 105222839A CN 201510519962 A CN201510519962 A CN 201510519962A CN 105222839 A CN105222839 A CN 105222839A
Authority
CN
China
Prior art keywords
max
range section
instrument coefficient
flow
preset relation
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
Application number
CN201510519962.2A
Other languages
English (en)
Other versions
CN105222839B (zh
Inventor
杨鸣
金定飞
徐永鹏
鲁法利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TANCY INSTRUMENT GROUP CO Ltd
Zhejiang Tancy Instrument Science & Technology Co ltd
Goldcard Smart Group Co Ltd
Original Assignee
Ningbo Pengsheng Technology Development Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ningbo Pengsheng Technology Development Co Ltd filed Critical Ningbo Pengsheng Technology Development Co Ltd
Priority to CN201510519962.2A priority Critical patent/CN105222839B/zh
Publication of CN105222839A publication Critical patent/CN105222839A/zh
Application granted granted Critical
Publication of CN105222839B publication Critical patent/CN105222839B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Volume Flow (AREA)

Abstract

本发明提供一种超声波流量计仪表系数的非线性修正方法。该方法包括:选取超声波流量计所在的管道内待测流体中的多个第一测量点的流量的理论值;测量流体在各第一测量点的流量的实际值和各第一测量点的流量的实际值对应的仪表系数的实际数值;根据各第一测量点选取各量程段;在各量程段内,根据各第一测量点的流量的实际值和各第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各量程段内的仪表系数的表达式;根据各量程段内的仪表系数的表达式计算各量程段内的仪表系数的计算值。采用本发明的超声波流量计仪表系数的计算方法减小了仪表系数的误差,进而使超声波流量计仪表系数计算结果更准确,超声波流量计测量更精确。

Description

超声波流量计仪表系数的非线性修正方法
技术领域
本发明涉及超声波技术领域,尤其涉及一种超声波流量计仪表系数的非线性修正方法。
背景技术
超声波流量计是近年来发展最为迅速的流量计之一,它是一种利用超声波信号在流体中传播时流体的流速信息来测量流体流量的测量技术,它具有非接触式测量、测量精度高、测量范围宽、安装维护方便等一些突出优点,常被应用到煤气、天然气、油气等气体的流量测量中。
时差法是采用超声波进行流量测量的常用方法之一,其原理是根据超声波信号顺流传播时间和逆流传播时间之差来计算流速,从而计算出流量。超声波流量计的仪表系数反应了超声气体流量计的仪表特性,其相对于平均仪表系数值的偏差程度表示了基本误差、线性度等决定仪表是否合格的关键指标。现有技术的超声波大多是将流速分布修正系数作为仪表系数,减小仪表系数的偏差,通过仪表系数把线平均速度折算成面积平均流速,减小超声波流量计的测量误差。
但是现有的流速分布修正系数是通过一些理论修正公式计算得来,这些公式是在各自作了某些假设的条件下导出,再用实验方法确定其中一些参数为常数,最后求得流速分布修正系数。而流速分布修正系数与流动剖面有关,流速分布修正系数为雷诺系数的函数,是实时变化的,按照前述的理论修正公式计算得到流速分布修正系数,会导致流速分布修正系数误差较大,即仪表系数误差较大。
发明内容
本发明提供一种超声波流量计仪表系数的非线性修正方法用以解决由于超声波流量计的仪表系数是实时变化的,而按现有理论公式计算得来的仪表系数导致误差较大的问题。
本发明提供一种超声波流量计仪表系数的非线性修正方法,所述方法包括:
选取所述超声波流量计所在的管道内待测流体中的多个第一测量点的流量的理论值;
测量所述流体在各所述第一测量点的流量的实际值和各所述第一测量点的流量的实际值对应的仪表系数的实际数值;
根据各所述第一测量点选取多个量程段;
在各所述量程段内,根据各所述第一测量点的流量的实际值和各所述第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各所述量程段内的仪表系数的表达式;
根据各所述量程段内的仪表系数的表达式计算各所述量程段内的所述仪表系数的计算值,并根据所述仪表系数的计算值实时动态的修正所述仪表系数。
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,各所述第一测量点的流量的理论值包括Qmin、Qt、0.1Qmax、0.1875Qmax、0.25Qmax、0.4Qmax、0.7Qmax和Qmax;其中Qmin为最小流量,Qt为分界流量,Qmax为最大流量。
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,在各所述量程段内,根据各所述第一测量点的流量的实际值和各所述第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各所述量程段内的仪表系数的表达式,具体包括:
选取各所述量程段对应的所述预设关系式;
对于每一个所述量程段,根据所述预设关系式中的常数的数量,在所述量程段内选取与所述常数数量相同的所述第一测量点的流量的实际值和对应的所述仪表系数的实际值;
将所述第一测量点的流量的实际值和对应的所述仪表系数的实际值代入所述预设关系式,求得所述预设关系式中的常数的数值;
将所述预设关系式中的所述常数的数值代入所述预设关系式,得到所述超声波流量计的所述仪表系数的表达式。
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,选取各所述量程段对应的所述预设关系式,具体包括:
当第一量程段为Qmin~0.1Qmax时:对应的所述预设关系式为: K = 10 ( a q + b ) + c ;
当第二量程段为0.1Qmax~0.25Qmax时,对应的所述预设关系式为:K=aq2+bq+c;
当第三量程段为0.25Qmax~0.4Qmax时,对应的所述预设关系式为:K=aq+b;
当第四量程段为0.4Qmax~0.7Qmax时,对应的所述预设关系式为:K=aq+b;以及
当第五量程段为0.7Qmax~Qmax时,对应的所述预设关系式为:K=aq+b。
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,对于每一个所述量程段,根据所述预设关系式中的常数的数量,在所述量程段内选取与所述常数数量相同的所述第一测量点的流量的实际值和对应的所述仪表系数的实际值,具体包括:
在所述第一量程段内选取所述Qmin、所述Qt和所述0.1Qmax三个所述测量点的第一测量数据;
在所述第二量程段内选取所述0.1Qmax、所述0.1875Qmax、所述0.25Qmax三个所述测量点的第二测量数据;
在所述第三量程段内选取所述0.25Qmax和所述0.4Qmax两个所述测量点的第三测量数据;
在所述第四量程段内选取所述0.4Qmax和所述0.7Qmax两个所述测量点的第四测量数据;以及
在所述第五量程段内选取所述0.7Qmax和所述Qmax两个所述测量点的第五测量数据。
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,将所述第一测量点的流量的实际值和对应的所述仪表系数的实际值代入所述预设关系式,求得所述预设关系式中的常数的数值,具体包括:
利用所述第一测量数据结合牛顿迭代法算出在所述第一量程段内的所述预设关系式中的常数a、b和c;
利用所述第二测量数据结合二次方程解法算出在所述第二量程段内的所述预设关系式中的常数a、b和c;
利用所述第三测量数据结合直线方程解法算出在所述第三量程段内的所述预设关系式中的常数a和b;
利用所述第四测量数据结合直线方程解法算出在所述第四量程段内的所述预设关系式中的常数a和b;以及
利用所述第五测量数据结合直线方程解法算出在所述第五量程段内的所述预设关系式中的常数a和b。
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,将所述预设关系式中的所述常数的数值代入所述预设关系式,得到所述超声波流量计的所述仪表系数的表达式,具体包括:
K = 10 ( a q + b ) + c Q min ≤ q ≤ 0.1 Q max aq 2 + b q + c 0.1 Q max ≤ q ≤ 0.25 Q max a q + b 0.25 Q max ≤ q ≤ 0.4 Q max a q + b 0.4 Q max ≤ q ≤ 0.7 Q max a q + b 0.7 Q max ≤ q ≤ Q max .
进一步地,上述所述超声波流量计仪表系数的非线性修正方法中,根据各所述量程段内的所述仪表系数的表达式计算各所述量程段内的所述仪表系数的计算值具体包括:
在所述各量程段内选取多个第二测量点的流量的理论值,各所述第二测量点包括各所述第一测量点;
将各所述第二测量点的流量的理论值代入各所述量程段内的所述仪表系数的表达式;
计算出各所述量程段内的各所述第二测量点对应的所述仪表系数的计算值。
本发明的超声波流量计仪表系数的非线性修正方法,通过选取超声波流量计的多个第一测量点的流量值,在不同量程内建立仪表系数的表达式,实现了对超声波流量计仪表系数的实时动态的修正。采用本发明的超声波流量计仪表系数的非线性修正方法减小了仪表系数的误差,进而使超声波流量计的仪表系数计算结果更准确,超声波流量计测量更精确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的超声波流量计仪表系数的非线性修正方法的实施例流程图;
图2为DN50型的超声波流量计的仪表系数K—流量Q的曲线与实测流量点分布图;
图3为DN50型的超声波流量计的流量仪表系数K相对示值误差图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明的超声波流量计仪表系数的非线性修正方法的实施例流程图。如图1所示,本实施例的超声波流量计仪表系数的非线性修正方法,具体可以包括如下步骤:
101、选取超声波流量计所在的管道内待测流体中的多个第一测量点的流量的理论值;
根据国际标准ISOCD17089、美国标准A.G.A.ReportNo.9和中国标准GB/T18604-2014提出的贸易计量的超声波气体流量计需进行实流校准要求,选取多个具有代表性的测量点。例如,各第一测量点的流量的理论值包括Qmin、Qt、0.1Qmax、0.1875Qmax、0.25Qmax、0.4Qmax、0.7Qmax和Qmax,其中Qmin为最小流量,Qt为分界流量,Qmax为最大流量。由于每种型号的超声波流量计的规格不同,超声波流量计的Qmin、和Qmax根据超声波流量计的出厂值确定,在气体流量中,Qt对应的流速应不大于3m/s,所以这个值是不一定的,需要计算得到,在使用之前就计算好,作为理论值使用。其它流量点的值优选为0.1Qmax、0.1875Qmax、0.25Qmax、0.4Qmax和0.7Qmax,可以适当的添加一个或者多个。
102、测量流体在各第一测量点的流量的实际值和各第一测量点的流量的实际值分别对应的仪表系数的实际数值;
选取的多个测量点为超声波流量计的理论值,实际使用中会存在偏差,测量出超声波气体流量计的实际流量值和多个测量点分别对应的仪表系的实际数值,以便于后面构造仪表系数和流量的关系。
103、根据各第一测量点选取多个量程段;
根据选取的多个第一测量点的流量的理论值选取两个测量点的流量的理论值作为一个量程段,两个测量点的流量的理论值之间可以包括其它测量点的理论流量值也可以不包括,相邻两个量程段中,后一个量程段的起点值为前一个量程段的终点值。
104、在各量程段内,根据各第一测量点的流量的实际值和各第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各量程段内的仪表系数的表达式;
将多个第一测量点的流量的实际值和各第一测量点的流量的实际值分别对应的仪表系数的实际数值代入预设关系式,求得各量程段内的仪表系数的表达式。
105、根据各量程段内的仪表系数的表达式计算各量程段内的仪表系数的计算值,并根据仪表系数的计算值实时动态的修正所述仪表系。
根据“步骤104”得到的仪表系数的表达式,计算超声波流量计在该流量计量程内的任一流量值对应的仪表系数,并根据仪表系数的计算值实时动态的修正所述仪表系,使超声波流量的仪表系数计算结果更准确,超声波流量计测量更精确。
本实施例的超声波流量计仪表系数的非线性修正方法,通过选取超声波流量计的多个第一测量点的流量值,在不同量程内建立仪表系数的表达式,实现了对超声波流量计仪表系数的实时动态的修正计算。采用本实施例的超声波流量计仪表系数的非线性修正方法减小了仪表系数的误差,进而使超声波流量的仪表系数计算结果更准确,超声波流量计测量更精确。
进一步地,在上述实施例中图1所示实施例的超声波流量计仪表系数的非线性修正方法的“步骤104”的基础上进一步更加详细地说明本发明的技术方案。“步骤104”在各量程段内,根据各第一测量点的流量的实际值和各第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各量程段内的仪表系数的表达式具体包括:
(1)选取各量程段对应的预设关系式;
将多个第一流量点分成五段,即五个量程段,在各量程段内分别用不同的曲线标定K。例如,选取各量程段对应的预设关系式,具体包括:
当第一量程段为Qmin~0.1Qmax时:对应的设关系式为:
当第二量程段为0.1Qmax~0.25Qmax时,对应的预设关系式为:K=aq2+bq+c;
当第三量程段为0.25Qmax~0.4Qmax时,对应的预设关系式为:K=aq+b;
当第四量程段为0.4Qmax~0.7Qmax时,对应的预设关系式为:K=aq+b;以及
当第五量程段为0.7Qmax~Qmax时,对应的预设关系式为:K=aq+b。
各量程段的预设关系式中a、b或者c分别为各预设关系式中的常数,但各常数的数值随着不同型号的超声波流量计而不同,因此a、b或者c的具体数值需根据不同型号的超声波流量计得到。其中K为仪表系数;q为流量。
需要说明的是,每个预设关系中a、b或者c在所有的关系中是独立的常数,任意两个或多个预设关系中的常数a、b或者c没有关系,需要根据实际情况计算a、b或者c,其结果可以相同也可以不同。
(2)对于每一个量程段,根据预设关系式中的常数的数量,在量程段内选取与常数数量相同的第一测量点的流量的实际值和对应的仪表系数的实际值;
例如,因为第一量程段内的预设公式中含有a、b和c三个常数,所以在第一量程段内选取Qmin、Qt和0.1Qmax三个测量点的第一测量数据;
同上,在第二量程段内选取0.1Qmax、0.1875Qmax、0.25Qmax三个测量点的第二测量数据;
同上,在第三量程段内选取0.25Qmax和0.4Qmax两个测量点的第三测量数据;
同上,在第四量程段内选取0.4Qmax和0.7Qmax两个测量点的第四测量数据;以及
同上,在第五量程段内选取0.7Qmax和Qmax两个测量点的第五测量数据。
(3)将第一测量点的流量的实际值和对应的仪表系数的实际值代入预设关系式,求得预设关系式中的常数的数值;
例如,对于五个量程段,具体采用如下方式计算:
利用第一测量数据结合牛顿迭代法可算出在第一量程段内的常数a、b和c;
利用第二测量数据结合二次方程解法可算出在第二量程段内的常数a和b和c;
利用第三测量数据结合直线方程解法可算出在第三量程段内的常数a和b;
利用第四测量数据结合直线方程解法可算出在第四量程段内常的数a和b;以及
利用第五测量数据结合直线方程解法可算出在第五量程段内的常数a和b。
(4)将预设关系式中的常数的数值代入预设关系式,得到超声波流量计的仪表系数的表达式。
经步骤(1)—步骤(3)所得到的仪表系数的表达式为
K = 10 ( a q + b ) + c Q min ≤ q ≤ 0.1 Q max aq 2 + b q + c 0.1 Q max ≤ q ≤ 0.25 Q max a q + b 0.25 Q max ≤ q ≤ 0.4 Q max a q + b 0.4 Q max ≤ q ≤ 0.7 Q max a q + b 0.7 Q max ≤ q ≤ Q max .
此仪表系数的表达式中a、b或者c,具体的为各型号的超声波流量计对应的具体数值。
进一步地,在上述实施例中图1所示实施例的超声波流量计仪表系数的非线性修正方法的“步骤105”的基础上进一步更加详细地说明本发明的技术方案。“步骤105”,根据各量程段内的仪表系数的表达式计算各量程段内的仪表系数的计算值,具体包括:
(11)在各量程段内选取多个第二测量点的流量的理论值,各所述第二测量点包括各所述第一测量点;
(12)将多个第二测量点的流量的理论值代入各量程段内的仪表系数的表达式;
(13)计算出各量程段内的各所述第二测量点对应的仪表系数的计算值。
在得到一种型号的超声波流量计的仪表系数的表达式后,可以选择该超声波流量计总量程内的任一流量点的理论值,确定好量程段,代入对应的表达式,即可得到该超声波流量计的仪表系数的计算值。
为了验证由此方法计算得来的超声波流量计仪表系数是否精确,在各量程段内选取多个第二测量点的流量的理论值,尽可能的多选,并且包括第一测量点的理论值。
测量出各量程段内选取多个第二测量点的流量的实际值和对应的仪表系数的实际值;并标记在仪表系数K与流量Q的坐标系中或者记录在表格内。
将计算出各量程段内的仪表系数的计算值用曲线平滑连接,得到仪表系数K与流量Q的坐标系中仪表系数K与流量Q的曲线或者与测量的多个第二测量点的流量的实际值和对应的仪表系数的实际值记录在同一表格内。根据仪表系数K与流量Q的曲线或者表格内数据比较实测值和计算值的大小,可以计算出偏差结果。
本实施例的超声波流量计仪表系数的非线性修正方法,通过选取超声波流量计的多个第一测量点的流量值作,在不同量程内建立仪表系数的表达式,实现了对超声波流量计仪表系数的实时动态的修正。采用本实施例的超声波流量计仪表系数的非线性修正方法减小了仪表系数的误差,进而使超声波流量计的仪表系数计算结果更准确,超声波流量计测量更精确。
进一步地,上述实施例的超声波流量计仪表系数的非线性修正方法,具体以DN50型的超声波流量计本发明的技术方案进一步更加详细地进行说明。
第一步:确定流量Qmin,Qt和Qmax的大小为Qmin=3m3/h,Qt=8m3/h,Qmax=160m3/h。
第二步:测量流量点Qmin、Qt、0.1Qmax、0.1875Qmax、0.25Qmax、0.4Qmax、0.7Qmax、Qmax的流量与K系数。见表1,表1为DN50型的超声波流量计的流量与K系数的实测值。
表1
第三步:将流量分成五个量程段分别用不同的公式拟合。
在第一量程段Qmin~0.1Qmax内,利用3、8和16三个流量点的数据结合牛顿迭代法计算得到该量程段内的公式。
在在第二量程段0.1Qmax~0.25Qmax内,利用16、30和40三个流量点的数据,利用方程组计算得到该量程段内的公式。
在第三量程段0.25Qmax~0.4Qmax内,利用40和64两个流量点数据,利用方程组计算得到该量程段内的公式。
在第四量程段0.4Qmax~0.7Qmax内,利用64和112两个流量点数据,利用方程组计算得到该量程段内的公式。
在第五量程段0.7Qmax~Qmax内,利用112和160两个流量点数据,利用方程组计算得到该量程段内的公式。
整个流量段的K系数可以表示为:
K = 10 ( 0.6583 q + 4.3198 ) + 45800.64 3 m 3 / h ≤ q ≤ 16 m 3 / h 1.34 q 2 - 127.109 q + 70445.746 16 m 3 / h ≤ q ≤ 40 m 3 / h - 17.576 q + 68204.406 40 m 3 / h ≤ q ≤ 64 m 3 / h - 9.721 q + 67696.235 64 m 3 / h ≤ q ≤ 112 m 3 / h - 3.789 q + 67031.84 112 m 3 / h ≤ q ≤ 160 m 3 / h .
使用同一型号的超声波流量计DN50测从Qmin~Qmax的任意一组数据进行K系数拟合误差验证。具体实验结果详见表2、表3、图2和图3。
表2为DN50型的超声波流量计的仪表系数K—流量Q的拟合计算值与实测数据误差第一分析表。
表2
表3为DN50型的超声波流量计的仪表系数K—流量Q的拟合计算值与实测数据误差第二分析表。
表3
由表2和表3中的相对示值误差可以看出,采用上述实施例的K与q的关系式拟合计算得到的K与K实测值误差很小,在整个量程内不会超过±0.5%。
图2为DN50型的超声波流量计的仪表系数K—流量Q的曲线与实测流量点分布图。如图2所示,仪表系数的实测分布点基本位于采用上述实施例的K与q的关系式拟合得到的仪表系数在整个量程内的曲线上。表明采用上述实施例的K与q的关系式拟合计算得到的K与K实测值误差很小。
图3为DN50型的超声波流量计的流量仪表系数K相对示值误差图。如图3所示,其中相对示值误差基本分布在±0.4%之间,采用上述实施例的K与q的关系式拟合计算得到的K与K实测值误差很小。
通常情况下,在超声波流量计的最小流量值和分界流量值之间为±1.4%,在超声波流量计的分界流量值和最大流量值之间为±1.0%。由此可见,本实施例的超声波流量计仪表系数的非线性修正方法减小了仪表系数的误差。
本实施例的超声波流量计仪表系数的非线性修正方法,通过选取超声波流量计的多个代表性的测量点的流量值作为基准值,在不同量程内建立仪表系数与流量的关系式,得到多个不同的关系式,实现了对超声波流量计仪表系数的实时状态的拟合计算。采用本实施例的超声波流量计仪表系数的非线性修正方法减小了仪表系数的误差,进而使超声波流量计的仪表系数计算结果更准确,超声波流量计测量更精确。
需要说明的是,上述实施例采用DN50型的超声波流量计为例对本发明的超声波流量计仪表系数的非线性修正方法进行详细说明的,其它型号的超声波流量计仪表系数的非线性修正方法同样可以采取本发明的超声波流量计仪表系数的非线性修正方法进行计算。其原理和计算步骤如上,在此不再赘述。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种超声波流量计仪表系数的非线性修正方法,其特征在于,所述方法包括:
选取所述超声波流量计所在的管道内待测流体中的多个第一测量点的流量的理论值;
测量所述流体在各所述第一测量点的流量的实际值和各所述第一测量点的流量的实际值对应的仪表系数的实际数值;
根据各所述第一测量点选取多个量程段;
在各所述量程段内,根据各所述第一测量点的流量的实际值和各所述第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各所述量程段内的仪表系数的表达式;
根据各所述量程段内的仪表系数的表达式计算各所述量程段内的所述仪表系数的计算值,并根据所述仪表系数的计算值实时动态的修正所述仪表系数。
2.根据权利要求1所述的方法,其特征在于,各所述第一测量点的流量的理论值包括Qmin、Qt、0.1Qmax、0.1875Qmax、0.25Qmax、0.4Qmax、0.7Qmax和Qmax;其中Qmin为最小流量,Qt为分界流量,Qmax为最大流量。
3.根据权利要求1所述的方法,其特征在于,在各所述量程段内,根据各所述第一测量点的流量的实际值和各所述第一测量点的流量的实际值对应的仪表系数的实际数值以及预设关系式,构造各所述量程段内的仪表系数的表达式,具体包括:
选取各所述量程段对应的所述预设关系式;
对于每一个所述量程段,根据所述预设关系式中的常数的数量,在所述量程段内选取与所述常数数量相同的所述第一测量点的流量的实际值和对应的所述仪表系数的实际值;
将所述第一测量点的流量的实际值和对应的所述仪表系数的实际值代入所述预设关系式,求得所述预设关系式中的常数的数值;
将所述预设关系式中的所述常数的数值代入所述预设关系式,得到所述超声波流量计的所述仪表系数的表达式。
4.根据权利要求3所述的方法,其特征在于,选取各所述量程段对应的所述预设关系式,具体包括:
当第一量程段为Qmin~0.1Qmax时:对应的所述预设关系式为: K = 10 ( a q + b ) + c ;
当第二量程段为0.1Qmax~0.25Qmax时,对应的所述预设关系式为:K=aq2+bq+c;
当第三量程段为0.25Qmax~0.4Qmax时,对应的所述预设关系式为:K=aq+b;
当第四量程段为0.4Qmax~0.7Qmax时,对应的所述预设关系式为:K=aq+b;以及
当第五量程段为0.7Qmax~Qmax时,对应的所述预设关系式为:K=aq+b。
5.根据权利要求4所述的方法,其特征在于,对于每一个所述量程段,根据所述预设关系式中的常数的数量,在所述量程段内选取与所述常数数量相同的所述第一测量点的流量的实际值和对应的所述仪表系数的实际值,具体包括:
在所述第一量程段内选取所述Qmin、所述Qt和所述0.1Qmax三个所述测量点的第一测量数据;
在所述第二量程段内选取所述0.1Qmax、所述0.1875Qmax、所述0.25Qmax三个所述测量点的第二测量数据;
在所述第三量程段内选取所述0.25Qmax和所述0.4Qmax两个所述测量点的第三测量数据;
在所述第四量程段内选取所述0.4Qmax和所述0.7Qmax两个所述测量点的第四测量数据;以及
在所述第五量程段内选取所述0.7Qmax和所述Qmax两个所述测量点的第五测量数据。
6.根据权利要求5所述的方法,其特征在于,将所述第一测量点的流量的实际值和对应的所述仪表系数的实际值代入所述预设关系式,求得所述预设关系式中的常数的数值,具体包括:
利用所述第一测量数据结合牛顿迭代法算出在所述第一量程段内的所述预设关系式中的常数a、b和c;
利用所述第二测量数据结合二次方程解法算出在所述第二量程段内的所述预设关系式中的常数a、b和c;
利用所述第三测量数据结合直线方程解法算出在所述第三量程段内的所述预设关系式中的常数a和b;
利用所述第四测量数据结合直线方程解法算出在所述第四量程段内的所述预设关系式中的常数a和b;以及
利用所述第五测量数据结合直线方程解法算出在所述第五量程段内的所述预设关系式中的常数a和b。
7.根据权利要求6所述的方法,其特征在于,将所述预设关系式中的所述常数的数值代入所述预设关系式,得到所述超声波流量计的所述仪表系数的表达式,具体包括:
K = 10 ( a q + b ) + c Q min ≤ q ≤ 0.1 Q max aq 2 + b q + c 0.1 Q max ≤ q ≤ 0.25 Q max a q + b 0.25 Q max ≤ q ≤ 0.4 Q max a q + b 0.4 Q max ≤ q ≤ 0.7 Q max a q + b 0.7 Q max ≤ q ≤ Q max .
8.根据权利要求1所述的方法,其特征在于,根据各所述量程段内的所述仪表系数的表达式计算各所述量程段内的所述仪表系数的计算值具体包括:
在所述各量程段内选取多个第二测量点的流量的理论值,各所述第二测量点包括各所述第一测量点;
将各所述第二测量点的流量的理论值代入各所述量程段内的所述仪表系数的表达式;
计算出各所述量程段内的各所述第二测量点对应的所述仪表系数的计算值。
CN201510519962.2A 2015-08-21 2015-08-21 超声波流量计仪表系数的非线性修正方法 Active CN105222839B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510519962.2A CN105222839B (zh) 2015-08-21 2015-08-21 超声波流量计仪表系数的非线性修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510519962.2A CN105222839B (zh) 2015-08-21 2015-08-21 超声波流量计仪表系数的非线性修正方法

Publications (2)

Publication Number Publication Date
CN105222839A true CN105222839A (zh) 2016-01-06
CN105222839B CN105222839B (zh) 2018-04-24

Family

ID=54991921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510519962.2A Active CN105222839B (zh) 2015-08-21 2015-08-21 超声波流量计仪表系数的非线性修正方法

Country Status (1)

Country Link
CN (1) CN105222839B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105571661A (zh) * 2016-02-23 2016-05-11 合肥工业大学 基于示值误差拟合的电磁流量计仪表特征系数的计算方法
CN105606146A (zh) * 2016-03-10 2016-05-25 江苏泰斯特电子设备制造有限公司 一种非线性修正技术
CN105758474A (zh) * 2016-05-30 2016-07-13 天津大学 一种提高气体超声流量计测量精度的方法
CN106679773A (zh) * 2017-02-07 2017-05-17 合肥工业大学 基于分段线性的电磁流量计仪表特征系数计算方法
CN108072428A (zh) * 2017-12-05 2018-05-25 成都秦川物联网科技股份有限公司 误差修正方法及燃气表
CN109932026A (zh) * 2017-12-19 2019-06-25 丹尼尔测量和控制公司 多流体校准
CN114485863A (zh) * 2021-12-24 2022-05-13 广东艾科技术股份有限公司 超声波水表的流量误差校正方法、系统、计算机及介质
CN115655423A (zh) * 2022-12-22 2023-01-31 成都千嘉科技股份有限公司 一种超声波燃气表的校准方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1330764A (zh) * 1998-10-26 2002-01-09 松下电器产业株式会社 设定流率系数的方法及利用该方法的流率计
US20050209793A1 (en) * 2004-02-27 2005-09-22 Fuji Electric Systems Co., Ltd. Doppler ultrasonic flowmeter, and processor and method thereof with quantization error correction
CN102435265A (zh) * 2011-09-19 2012-05-02 苏州聚元微电子有限公司 一种时差法超声波流量计精度修正的方法
CN102538913A (zh) * 2011-12-26 2012-07-04 宁波水表股份有限公司 单声道超声水表流量测量特性校正方法
CN102879061A (zh) * 2012-09-19 2013-01-16 宁波水表股份有限公司 一种基于拟合方程的水表误差校正方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1330764A (zh) * 1998-10-26 2002-01-09 松下电器产业株式会社 设定流率系数的方法及利用该方法的流率计
US20050209793A1 (en) * 2004-02-27 2005-09-22 Fuji Electric Systems Co., Ltd. Doppler ultrasonic flowmeter, and processor and method thereof with quantization error correction
CN102435265A (zh) * 2011-09-19 2012-05-02 苏州聚元微电子有限公司 一种时差法超声波流量计精度修正的方法
CN102538913A (zh) * 2011-12-26 2012-07-04 宁波水表股份有限公司 单声道超声水表流量测量特性校正方法
CN102879061A (zh) * 2012-09-19 2013-01-16 宁波水表股份有限公司 一种基于拟合方程的水表误差校正方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105571661A (zh) * 2016-02-23 2016-05-11 合肥工业大学 基于示值误差拟合的电磁流量计仪表特征系数的计算方法
CN105571661B (zh) * 2016-02-23 2018-08-24 合肥工业大学 基于示值误差拟合的电磁流量计仪表特征系数的计算方法
CN105606146A (zh) * 2016-03-10 2016-05-25 江苏泰斯特电子设备制造有限公司 一种非线性修正技术
CN105758474A (zh) * 2016-05-30 2016-07-13 天津大学 一种提高气体超声流量计测量精度的方法
CN105758474B (zh) * 2016-05-30 2019-01-08 天津大学 一种提高气体超声流量计测量精度的方法
CN106679773A (zh) * 2017-02-07 2017-05-17 合肥工业大学 基于分段线性的电磁流量计仪表特征系数计算方法
CN108072428A (zh) * 2017-12-05 2018-05-25 成都秦川物联网科技股份有限公司 误差修正方法及燃气表
CN109932026A (zh) * 2017-12-19 2019-06-25 丹尼尔测量和控制公司 多流体校准
CN114485863A (zh) * 2021-12-24 2022-05-13 广东艾科技术股份有限公司 超声波水表的流量误差校正方法、系统、计算机及介质
CN115655423A (zh) * 2022-12-22 2023-01-31 成都千嘉科技股份有限公司 一种超声波燃气表的校准方法
CN115655423B (zh) * 2022-12-22 2023-03-21 成都千嘉科技股份有限公司 一种超声波燃气表的校准方法

Also Published As

Publication number Publication date
CN105222839B (zh) 2018-04-24

Similar Documents

Publication Publication Date Title
CN105222839A (zh) 超声波流量计仪表系数的非线性修正方法
CN103808381B (zh) 一种时差式超声波流量计的温度影响消除方法
CN104040299B (zh) 流量计的管线内验证方法
US9816847B2 (en) Ultrasonic flowmeter and method for measuring flow
RU2478919C2 (ru) Способ калибровки в условиях эксплуатации ультразвуковых расходомеров-счетчиков расхода и объема жидких однофазных сред
CN105043511A (zh) 基于雷诺数的涡街流量计系数修正方法
EP2909587B1 (en) Ultrasonic flow metering system with an upstream pressure transducer
JP5942085B2 (ja) 流量補正係数設定方法とこれを用いた流量計測装置
EP4220098A1 (en) Flow metering system condition-based monitoring and failure to predictive mode
CN105603955B (zh) 一种基于实测数据的剩余沉降计算方法及应用
CN102435265A (zh) 一种时差法超声波流量计精度修正的方法
EA039878B1 (ru) Способ калибровки пустотной фракции
CN102538913B (zh) 单声道超声水表流量测量特性校正方法
CN104806889A (zh) 天然气长输管道管输损耗监测方法及装置
Provenzano et al. Assessing a local losses evaluation procedure for low-pressure lay-flat drip laterals
CN105571661A (zh) 基于示值误差拟合的电磁流量计仪表特征系数的计算方法
Rhinehart et al. A power law approach to orifice flow rate calibration
DE19605652A1 (de) Verfahren zur Kalibrierung eines Durchflußmessers
EP3889555B1 (en) Flow meter calibration method
US10809107B2 (en) Multi-fluid calibration
Oriji et al. Application of metering process in oil and gas production in Niger delta fields
Ribeiro et al. Measurement uncertainty of turbine flow meter calibration used in conformity assessment for water management
KR101585175B1 (ko) 병렬 운전 펌프의 소모 전력 절감 운용 장치 및 그 방법
RU2533329C1 (ru) Установка для поверки и калибровки счетчиков, расходомеров и расходомеров-счетчиков газа
Henry et al. Towards assessing online uncertainty for three-phase flow metering in the oil and gas industry

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20161215

Address after: Jinqiao Hangzhou economic and Technological Development Zone Hangzhou city Zhejiang province 310000 Street No. 158 workshop room 1-401

Applicant after: Zhejiang Tianxin Ultrasonic Technology Co.,Ltd.

Address before: Hongtang town Jiangbei District Ningbo City, Zhejiang province 315038 City Qiu

Applicant before: NINGBO PENGSHENG TECHNOLOGY DEVELOPMENT Co.,Ltd.

GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 310018 workshop 1, 158 Jinqiao street, Qiantang New District, Hangzhou City, Zhejiang Province - Room 401

Patentee after: Tianxin instrument group (Hangzhou) Co.,Ltd.

Address before: 310000 room 1-401, workshop 158, Jinqiao street, Hangzhou Economic and Technological Development Zone, Hangzhou City, Zhejiang Province

Patentee before: Zhejiang Tianxin Ultrasonic Technology Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220729

Address after: No. 198, Hualian Road, Cangnan Industrial Park, Wenzhou City, Zhejiang Province 325800

Patentee after: Tancy Instrument Group Co.,Ltd.

Address before: 310018 workshop 1, 158 Jinqiao street, Qiantang New District, Hangzhou City, Zhejiang Province - Room 401

Patentee before: Tianxin instrument group (Hangzhou) Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231117

Address after: No. 198, Hualian Road, Cangnan Industrial Park, Wenzhou City, Zhejiang Province 325800

Patentee after: Tancy Instrument Group Co.,Ltd.

Patentee after: ZHEJIANG TANCY INSTRUMENT SCIENCE & TECHNOLOGY Co.,Ltd.

Address before: No. 198, Hualian Road, Cangnan Industrial Park, Wenzhou City, Zhejiang Province 325800

Patentee before: Tancy Instrument Group Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240426

Address after: No. 198, Hualian Road, Cangnan Industrial Park, Wenzhou City, Zhejiang Province 325800

Patentee after: Tancy Instrument Group Co.,Ltd.

Country or region after: China

Patentee after: Goldcard Smart Group Co.,Ltd.

Patentee after: ZHEJIANG TANCY INSTRUMENT SCIENCE & TECHNOLOGY Co.,Ltd.

Address before: No. 198, Hualian Road, Cangnan Industrial Park, Wenzhou City, Zhejiang Province 325800

Patentee before: Tancy Instrument Group Co.,Ltd.

Country or region before: China

Patentee before: ZHEJIANG TANCY INSTRUMENT SCIENCE & TECHNOLOGY Co.,Ltd.