CN109342769B - 一种标定方法、流速测量方法及装置 - Google Patents

一种标定方法、流速测量方法及装置 Download PDF

Info

Publication number
CN109342769B
CN109342769B CN201811506650.8A CN201811506650A CN109342769B CN 109342769 B CN109342769 B CN 109342769B CN 201811506650 A CN201811506650 A CN 201811506650A CN 109342769 B CN109342769 B CN 109342769B
Authority
CN
China
Prior art keywords
formula
flow
coefficient
flow velocity
calibration
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
Application number
CN201811506650.8A
Other languages
English (en)
Other versions
CN109342769A (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.)
Laser Institute of Shandong Academy of Science
Original Assignee
Laser Institute of Shandong Academy of Science
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 Laser Institute of Shandong Academy of Science filed Critical Laser Institute of Shandong Academy of Science
Priority to CN201811506650.8A priority Critical patent/CN109342769B/zh
Publication of CN109342769A publication Critical patent/CN109342769A/zh
Application granted granted Critical
Publication of CN109342769B publication Critical patent/CN109342769B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P21/00Testing or calibrating of apparatus or devices covered by the preceding groups
    • G01P21/02Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
    • G01P21/025Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers for measuring speed of fluids; for measuring speed of bodies relative to fluids

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Details Of Flowmeters (AREA)

Abstract

本发明提供了一种标定方法、流速测量方法及装置,涉及光纤传感技术领域。所述方法包括:选定流速传感器置于确定的流场的参数;获取所述流速固定参数;最后根据所述流场的参数与所述固定参数,计算所述流速传感器的标定参数。通过获得标定参数并利用标定公式,避免了根据传热方程进行复杂微积分计算,减少了运算量,更利于工程利用;并且利用标定参数的方式,解决了常规方法、默认定律公式带来的误差问题。

Description

一种标定方法、流速测量方法及装置
技术领域
本发明涉及光纤传感技术领域,具体而言,涉及一种标定方法、流速测量方法及装置。
背景技术
目前,风速测量在矿山瓦斯抽采、巷道通风、采空区发火防控等监测领域占据着举足轻重的地位。为了获得准确的风速测量数据,就需要对风速传感器的一些参数进行标定,而现有的标定方法则是通过采用现有的多个复杂的常规公式计算流速传感器的一些设置参数,如经典King定律以及传热方程等,增加了计算的复杂性,并且因为多个公式的组合也导致误差的增大。
发明内容
为了克服上述现有技术中的不足,本发明提供一种标定方法及装置,用于改善上述问题。
为了实现上述目的,本发明实施例所提供的技术方案如下所示:
第一方面,本发明实施例提供一种标定方法,包括:获取所述流速传感器所放置的流场的参数;获取固定参数;根据所述流场的参数与所述固定参数,计算得到测风光栅波长漂移量Δλ以及光栅温度系数c,并结合流速标定公式计算得到所述流速传感器的无量纲比值β;
所述流场的参数包括:所述流速传感器处于热稳态状态时在预定流速u下产生的热量P;所述固定参数包括:所述流速传感器的圆管特征长度d、运动粘度v、流体导热系数k和所述流速传感器的探头面积S;
所述流速标定公式为:
Figure GDA0002663208620000021
其中,m1、m2以及m3分别为预先计算得到的所述流速传感器的第一标定系数、第二标定系数及第三标定系数,n为预先计算得到的定值。
结合上述第一方面提供的技术方案,在一些可能的实现方式中,所述标定方法还包括:根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure GDA0002663208620000022
对应的,通过以下步骤计算所述第一标定系数m1与第二标定系数m2:计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2;根据无量纲比值公式得到不同流速u1下的第一无量纲比值β1,所述无量纲比值公式为:
Figure GDA0002663208620000023
将不同流速u1下的所述第一无量纲比值β1与所述流速传感器在不同流速u1下的雷诺数Re拟合,得到不同流速u1下的第一拟合公式:
Figure GDA0002663208620000024
根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
结合上述第一方面提供的技术方案,在一些可能的实现方式中,所述标定方法还包括:计算所述流速传感器在不同流速u1下的换热系数h1:根据公式:
Figure GDA0002663208620000025
计算所述换热系数h1,其中,ΔT为所述流速传感器的初始温度与在不同流速u1下处于所述热稳态状态时的温度之间的温度变化量。
结合上述第一方面提供的技术方案,在一些可能的实现方式中,所述标定方法还包括:计算所述流速传感器在不同流速u1下的换热系数h2包括:根据公式
Figure GDA0002663208620000031
计算在不同流速u1下的努塞尔数Nu,其中,Pr为普朗特数;根据公式
Figure GDA0002663208620000032
计算得到所述换热系数h2
结合上述第一方面提供的技术方案,在一些可能的实现方式中,所述标定方法还包括:通过以下步骤计算n和所述第三标定系数m3:将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren;根据所述第二拟合公式,得到所述第三标定系数m3与n的数值。
结合上述第一方面提供的技术方案,在一些可能的实现方式中,所述标定方法还包括:通过计算公式:c=1-α计算所述光栅温度系数c,所述α为所述流速传感器的测风光栅和测温光栅温度系数比。
结合上述第一方面提供的技术方案,在一些可能的实现方式中,所述标定方法还包括:通过计算公式:Δλ=λ12计算所述测风光栅波长漂移量Δλ,所述λ1为所述流速传感器的测风光栅在流速u下对应的波长,所述λ2为所述流速传感器的温度补偿光栅在流速u下对应的波长。
第二方面,本发明实施例还提供一种流速测量方法,包括:获取所述流速传感器测得的热量P和测风光栅波长漂移量Δλ;获取固定参数,固定参数包括:所述流速传感器的光栅温度系数c、圆管特征长度d、运动粘度v、流体导热系数k、所述流速传感器的探头面积S、无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3以及定值n;根据流速标定公式:
Figure GDA0002663208620000033
计算得到流速u。
第三方面,本发明实施例还提供一种流速传感器标定装置,包括:获取单元及处理单元。
获取单元用于获取所述流速传感器所放置的流场的参数,所述流场的参数包括:所述流速传感器处于热稳态状态时在预定流速u下产生的热量P和测风光栅波长漂移量Δλ;并且还用于获取固定参数,所述固定参数包括:所述流速传感器的圆管特征长度d、运动粘度v、流体导热系数k和所述流速传感器的探头面积S;
处理单元用于根据所述流场的参数与所述固定参数,计算得到测风光栅波长漂移量Δλ以及光栅温度系数c,并结合流速标定公式计算得到所述流速传感器的无量纲比值β,所述流速标定公式为:
Figure GDA0002663208620000041
其中,m1、m2以及m3分别为预先计算得到的所述流速传感器的第一标定系数、第二标定系数及第三标定系数,n为预先计算得到的定值。
第四方面,本发明实施例还提供一种流速测量装置,包括测量单元、获取单元及处理单元。
测量单元,用于测得所述流速传感器产生的热量P以及测风光栅波长漂移量Δλ;
获取单元用于获取固定参数,所述固定参数包括:所述流速传感器的光栅温度系数c、圆管特征长度d、流体导热系数k、所述流速传感器的探头面积S、运动粘度v、无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3以及定值n;
处理单元用于根据所述流速传感器产生的热量P、所述测风光栅波长漂移量Δλ与所述固定参数,并结合流速标定公式计算得到所述流速传感器所放置的流场内的流速u,所述流速标定公式为:
Figure GDA0002663208620000042
本发明的有益效果是:
本发明通过提出一种新的流速标定公式,通过将需要标定的流速传感器放置于预先设定好了参数的流场中,根据获取到的参数(流场的参数和固定参数)和流速标定公式获得了所述流速传感器的无量纲比值β,即所述流速传感器的标定参数,完成了对所述流速传感器的标定。因为新的流速标定公式中需要的参数基本上是可控的(例如风速u下的热量)或者可以精确得到的参数(固定参数),避免了现有技术中使用常规方法、默认定律公式带来的复杂性及误差问题。
为使本发明的上述目的、特征和有点能更明显易懂,下文举出本发明实施例,并配合所附附图,作详细说明。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例1提供的一种标定方法的流程示意图;
图2为本发明实施例1提供的第一标定系数m1与第二标定系数m2的获取方法;
图3本发明实施例1提供的雷诺数与换热系数比关系图;
图4为本发明实施例1提供的雷诺数与努塞尔数标定关系图;
图5为本实施例2提供的一种流速测量方法的流程示意图;
图6为本发明实施例3提供的一种流速传感器标定装置的功能框图;
图7为本发明实施例4提供的一种流速测量装置的功能框图。
具体实施方式
本发明实施例中的技术方案为解决上述的技术问题,总体思路如下:
本发明通过预先设定一个确定的流场,将要标定的流速传感器放置于预先设定好的流场中,获取流速标定公式中除了无量纲比值外的各项参数,进而根据流速标定公式得出所述流速传感器的无量纲比值,便于所述流速传感器的工业应用,避免了传热方程进行微积分计算的复杂性,并且利用标定参数解决了常规方法、默认定律公式带来的误差问题。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
在本发明实施例的描述中,需要说明的是,指示方位或位置关系为基于附图所示的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,或者是本领域技术人员惯常理解的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
实施例1
请参考图1,图1为本发明实施例1提供的标定方法的流程示意图。本实施例提供的标定方法应用于流速传感器中,所述流速传感器包括风速传感器、液体流速传感器等。具体的,在本实施例中,以光纤热线风速传感器为例进行说明。
如图1所示,该标定方法包括:
步骤S101:获取所述流速传感器所放置的流场的参数;所述流场的参数包括:所述流速传感器处于热稳态状态时在预定流速u下产生的热量P;
步骤S102:获取固定参数;所述固定参数包括:所述流速传感器的圆管特征长度d、运动粘度v、流体导热系数k和所述流速传感器的探头面积S;
步骤S103:根据所述流场的参数与所述固定参数,计算得到测风光栅波长漂移量Δλ以及光栅温度系数c,并结合流速标定公式计算得到所述流速传感器的无量纲比值β,所述流速标定公式为:
Figure GDA0002663208620000071
需要说明的是,m1、m2以及m3分别为预先计算得到的所述流速传感器的第一标定系数、第二标定系数及第三标定系数,n为预先计算得到的定值。
通过先设定一个确定的流场,将所述流速传感器设置在所述确定的流场中从而获得一些测得的参数,进而根据所述的流速标定公式得到标定参数,解决了常规方法及默认定律公式带来的误差问题,并且减少了很多的运算量,避免了使用努赛尔公式等的复杂性,提高了传感器工作效率。
接下来将详细描述各个步骤的具体实施过程。
在步骤S101中,需要获取流速传感器处于热稳态状态时在预定流速u下产生的热量P,本领域技术人员知悉,当传感器置于流场中时,并且处于热稳态状态下,光纤热线产生的热量P等于流速传感器在该风速下所有散热量Q;所以在本实施例中,可以获取流速传感器在预定流速u下的所有散热量Q作为热量P。当然,在其它实施例中,还可以通过其它方式获取流速传感器在预定流速u下产生的热量P,例如在流速传感器上设置热量检测元件,本申请不作具体限定。
在步骤S102中获取的固定参数中,所述光栅温度系数c、圆管特征长度d以及探头面积S对于一个确定的传感器为定值,运动粘度v在一定温度范围内可以认为是一定值,在本实施例中,所述v与所述流体导热系数k为一个定值。
便于更清楚的描述本实施例中的标定方法的实施过程,在介绍步骤S103之前,下面先介绍步骤S103中使用到的流速标定公式中的第一标定系数m1、第二标定系数m2、第三标定系数m3和n的获取过程。
可选地,请参考图2,图2为本发明实施例1提供的第一标定系数m1与第二标定系数m2的获取方法。所述第一标定系数m1与第二标定系数m2的获取方法包括:
步骤S201:根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure GDA0002663208620000081
步骤S202:计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2
步骤S203:根据无量纲比值公式得到第一无量纲比值β1,所述无量纲比值公式为:
Figure GDA0002663208620000091
步骤S204:将所述第一无量纲比值β1与所述流速传感器在不同流速u1下的雷诺数Re拟合,得到第一拟合公式:
Figure GDA0002663208620000092
步骤S205:根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
需要说明的是,在本实施例中,所述流速u为一个定值,所述不同流速u1用于表达通过改变流速u统计一些参数在不同流速u下的数值,因而所述不同流速u1代表多个不同的流速u,为许多个流速u的集合,因而所述流速u为所述不同流速u1中的一个值。
请参考图3,图3为本发明实施例1提供的雷诺数与换热系数比关系图。从图3可以看出,对于本实施例1所述的光纤热线风速传感器,通过计算并记录不同风速u1下的第一无量纲比值β1与雷诺数Re,并将记录的第一无量纲比值β1与雷诺数Re的对应数值做成坐标图进行拟合,在相关系数R2=0.9986时,可以得到m1=0.0006,m2=-0.33,从而得到了本实施例中的光纤热线风速传感器的第一标定系数和第二标定系数。
通过将两个换热系数的无量纲比值进行拟合,避免了根据传热方程进行微积分计算的复杂性,提高了运算效率。
可选地,步骤S202包括:根据公式:
Figure GDA0002663208620000093
计算所述换热系数h1,其中,ΔT为所述流速传感器的初始温度与在所述风速u下热稳态状态时的温度之间的温度变化量。
可选地,步骤S202还包括:根据公式:
Figure GDA0002663208620000101
计算在不同流速u1下的努塞尔数Nu,其中Pr为普朗特数;之后根据公式
Figure GDA0002663208620000102
计算得到所述换热系数h2
需要说明的是,所述普朗特数Pr在本实施例中的变化可以忽略不计,均为一个定值。
所述计算得到定值n和第三标定系数m3的方法为:将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren;之后根据获得的第二拟合公式,得到第三标定系数m3与n的数值。
请参考图4,图4为本发明实施例1提供的雷诺数与努塞尔数标定关系图。从图4可以看出,对于本实施例1所述的光纤热线风速传感器,通过计算并记录不同风速u1下的雷诺数Re与努塞尔数Nu,将记录的雷诺数Re与努塞尔数Nu的对应数值做成坐标图进行最小二乘法拟合,在相关系数R2=0.9999时,可以得到m3=0.9138,n=0.4805,从而得到了本实施例中的光纤热线风速传感器的第三标定系数m3与n的值。
针对雷诺数Re和努塞尔数Nu,根据已测风速进行理论计算和最小二乘法拟合,简化了雷诺数与努赛尔数之间的关系;并且由于是针对King定律标定指数参数的不定值拟合,因而也使得公式更加准确。
可选地,步骤S103还包括通过计算公式:c=1-α计算所述光栅温度系数c,所述α为所述流速传感器的测风光栅和测温光栅温度系数比。对一个确定的光纤热线风速传感器,所述测风光栅和测温光栅温度系数的比值为一个定值,所以本发明实施例中的光栅温度系数c是一个定值。
可选地,步骤S103还包括通过计算公式:Δλ=λ12计算所述测风光栅波长漂移量Δλ,所述λ1为所述传感器测风光栅在流速u下对应的波长,所述λ2为所述传感器温度补偿光栅在流速u下对应的波长。
需要说明的是,在本实施例中,所述测风光栅波长漂移量Δλ的值通过光纤光栅解调仪通过上述计算公式得出,在其他实施例中,也可以采用如光谱仪等仪器,本申请不作具体限定。
实施例2
本发明实施例2提供一种流速测量方法。请参考图5,图5为本实施例2提供的流速测量方法的流程示意图。所述流速测量方法同样应用于流速传感器中,在本实施例2中同样采用光纤热线风速传感器。所述流速测量方法包括:
步骤S501:测得所述流速传感器产生的热量P以及测风光栅波长漂移量Δλ;
步骤S502:获取固定参数;
步骤S503:根据流速标定公式:
Figure GDA0002663208620000111
计算得到流速u。
其中,所述固定参数包括:所述流速传感器的光栅温度系数c、圆管特征长度d、运动粘度v、流体导热系数k、所述流速传感器的探头面积S、无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3以及定值n。需要说明的是,所述无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3以及定值n均通过实施例1获得。
需要说明的是,与实施例1相同,在步骤S501中,需要获取流速传感器处于热稳态状态时在预定流速u下产生的热量P,本领域技术人员知悉,当传感器置于流场中时,并且处于热稳态状态下,光纤热线产生的热量P等于流速传感器在该风速下所有散热量Q;所以在本实施例中,可以获取流速传感器在预定流速u下的所有散热量Q作为热量P。所述测风光栅波长漂移量Δλ的测量获取方式与实施例1相同,本实施例在此不再赘述。
除了步骤S501外,在步骤S502所述的固定参数中,所述光栅温度系数c、圆管特征长度d以及探头面积S对于一个确定的传感器为定值,运动粘度v在一定温度范围内可以认为是一定值,在本实施例中,所述v与所述流体导热系数k为一个定值。
本发明实施例2通过首先获取测量得到的参数以及流速传感器的固定参数,再仅通过流速标定公式,即获取得到了流速u,简化了流速u的计算获取过程,更利于工程应用。
实施例3
基于与实施例1的同一发明构思,本发明实施例还提供一种流速传感器标定装置,用于执行如图1所示的标定方法。请参考图6,图6为流速传感器标定装置的功能框图。与实施例1相同的是,所述流速传感器标定装置同样应用于流速传感器上,在本实施例中,所述流速传感器同样为光纤热线风速传感器。
所述流速传感器标定装置包括:获取单元11与处理单元12。获取单元11既用于获取所述流速传感器所放置的流场的参数,还用于获取固定参数。处理单元12用于根据所述流场的参数与所述固定参数,计算得到测风光栅波长漂移量Δλ以及光栅温度系数c,进而根据流速标定公式计算得到所述流速传感器的无量纲比值β。
具体地,在本实施例中所述流场的参数包括:所述流速传感器处于热稳态状态时在预定流速u下产生的热量P和测风光栅波长漂移量Δλ。
在本实施例中所述固定参数包括:所述流速传感器的圆管特征长度d、运动粘度v、流体导热系数k和所述流速传感器的探头面积S。
可选地,所述处理单元12还用于根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure GDA0002663208620000131
对应的,所述处理单元12还用于通过以下步骤计算所述第一标定系数m1与第二标定系数m2
计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2
根据无量纲比值公式得到不同流速u1下的第一无量纲比值β1,所述无量纲比值公式为:
Figure GDA0002663208620000132
将不同流速u1下的所述第一无量纲比值β1与所述流速传感器在的雷诺数Re拟合,得到不同流速u1下的第一拟合公式:
Figure GDA0002663208620000133
根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
可选地,所述处理单元12还用于根据公式:
Figure GDA0002663208620000134
计算所述换热系数h1
可选地,所述处理单元12还用于根据根据公式
Figure GDA0002663208620000135
计算在不同流速u1下的努塞尔数Nu;并且根据公式
Figure GDA0002663208620000136
所述处理单元12计算得到在不同风速u1所述换热系数h2
可选地,所述处理单元12还用于将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren;之后根据获得的第二拟合公式,得到第三标定系数m3与n的数值。
可选地,所述处理单元12还用于通过计算公式:c=1-α计算所述光栅温度系数c。
可选地,所述处理单元12还用于通过计算公式:Δλ=λ12计算所述测风光栅波长漂移量Δλ。
本实施例中的标定装置与实施例1所述的标定方法是基于同一构思下的发明,通过前述对标定方法及其各种变化形式的详细描述,本领域技术人员可以清楚的了解本实施例中标定装置的实施过程,所以为了说明书的简洁,在此不再赘述。
实施例4
基于与实施例2的同一发明构思,本发明实施例还提供一种流速测量装置,用于执行如图5所示的标定方法。请参考图7,图7为本发明实施例4提供的一种流速测量装置的功能框图。与实施例3相同的是,所述流速传感器标定装置同样应用于流速传感器上,在本实施例中,所述流速传感器同样为光纤热线风速传感器。
所述流速测量装置包括:测量单元21、获取单元22及处理单元20。所述测量单元21用于测得所述流速传感器产生的热量P以及测风光栅波长漂移量Δλ;所述获取单元22用于获取固定参数;所述处理单元20用于根据所述流速传感器产生的热量P、所述测风光栅波长漂移量Δλ与所述固定参数,根据实施例2中所述的流速标定公式计算得到所述流速传感器所放置的流场内的流速u。
本实施例中的流速测量装置与实施例2所述的流速测量方法是基于同一构思下的发明,通过前述对流速测量方法及其各种变化形式的详细描述,本领域技术人员可以清楚的了解本实施例中流速测量设备的实施过程,所以为了说明书的简洁,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种标定方法,其特征在于,应用于流速传感器,所述方法包括:
获取所述流速传感器所放置的流场的参数,所述流场的参数包括:所述流速传感器处于热稳态状态时在预定流速u下产生的热量P;
获取固定参数,所述固定参数包括:所述流速传感器的圆管特征长度d、运动粘度v、流体导热系数k和所述流速传感器的探头面积S;
根据所述流场的参数与所述固定参数,计算得到测风光栅波长漂移量Δλ以及光栅温度系数c,并结合流速标定公式计算得到所述流速传感器的无量纲比值β,所述流速标定公式为:
Figure FDA0002663208610000011
其中,m1、m2以及m3分别为预先计算得到的所述流速传感器的第一标定系数、第二标定系数及第三标定系数,n为预先计算得到的定值;
所述测风光栅波长漂移量Δλ、光栅温度系数c、第一标定系数m1、第二标定系数m2、第三标定系数m3以及n通过以下述方法计算得到;
根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure FDA0002663208610000012
计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2
根据公式:
Figure FDA0002663208610000013
计算所述换热系数h1,其中,ΔT为所述流速传感器的初始温度与在不同流速u1下热稳态状态时的温度之间的温度变化量;
根据公式
Figure FDA0002663208610000021
计算在不同流速u1下的努塞尔数Nu,其中,Pr为普朗特数;
根据公式
Figure FDA0002663208610000022
计算得到在不同风速u1所述换热系数h2
根据无量纲比值公式得到不同流速u1下的第一无量纲比值β1,所述无量纲比值公式为:
Figure FDA0002663208610000023
将不同流速u1下的所述第一无量纲比值β1与所述流速传感器在的雷诺数Re拟合,得到不同流速u1下的第一拟合公式:
Figure FDA0002663208610000024
根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren
根据所述第二拟合公式,得到所述第三标定系数m3与n的数值;
通过计算公式:c=1-α计算所述光栅温度系数c,所述α为所述流速传感器的测风光栅和测温光栅温度系数比;
通过计算公式:Δλ=λ12计算所述测风光栅波长漂移量Δλ,所述λ1为所述流速传感器的测风光栅在流速u下对应的波长,所述λ2为所述流速传感器的温度补偿光栅在流速u下对应的波长。
2.一种流速测量方法,其特征在于,应用于流速传感器,所述方法包括:
测得所述流速传感器产生的热量P以及测风光栅波长漂移量Δλ;
获取固定参数,所述固定参数包括:所述流速传感器的光栅温度系数c、圆管特征长度d、流体导热系数k、所述流速传感器的探头面积S、运动粘度v、无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3以及定值n;
根据流速标定公式:
Figure FDA0002663208610000031
计算得到流速u;
所述无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3、定值n通过以下述方法计算得到;
根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure FDA0002663208610000032
计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2
根据公式:
Figure FDA0002663208610000033
计算所述换热系数h1,其中,ΔT为所述流速传感器的初始温度与在不同流速u1下热稳态状态时的温度之间的温度变化量;
根据公式
Figure FDA0002663208610000034
计算在不同流速u1下的努塞尔数Nu,其中,Pr为普朗特数;
根据公式
Figure FDA0002663208610000035
计算得到在不同风速u1所述换热系数h2
根据无量纲比值公式得到不同流速u1下的第一无量纲比值β1,所述无量纲比值公式为:
Figure FDA0002663208610000036
将不同流速u1下的所述第一无量纲比值β1与所述流速传感器在的雷诺数Re拟合,得到不同流速u1下的第一拟合公式:
Figure FDA0002663208610000037
根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren
根据所述第二拟合公式,得到所述第三标定系数m3与n的数值。
3.一种流速传感器标定装置,其特征在于,所述标定装置包括:
获取单元,用于所述流速传感器所放置的流场的参数,所述流场的参数包括:所述流速传感器处于热稳态状态时在预定流速u下产生的热量P;
所述获取单元还用于获取固定参数,所述固定参数包括:所述流速传感器的圆管特征长度d、运动粘度v、流体导热系数k和所述流速传感器的探头面积S;
处理单元,用于根据所述流场的参数与所述固定参数,计算得到测风光栅波长漂移量Δλ以及光栅温度系数c,并结合流速标定公式计算得到所述流速传感器的无量纲比值β,所述流速标定公式为:
Figure FDA0002663208610000041
所述测风光栅波长漂移量Δλ、光栅温度系数c、第一标定系数m1、第二标定系数m2、第三标定系数m3以及n通过以下述方法计算得到;
根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure FDA0002663208610000042
计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2
根据公式:
Figure FDA0002663208610000043
计算所述换热系数h1,其中,ΔT为所述流速传感器的初始温度与在不同流速u1下热稳态状态时的温度之间的温度变化量;
根据公式
Figure FDA0002663208610000051
计算在不同流速u1下的努塞尔数Nu,其中,Pr为普朗特数;
根据公式
Figure FDA0002663208610000052
计算得到在不同风速u1所述换热系数h2
根据无量纲比值公式得到不同流速u1下的第一无量纲比值β1,所述无量纲比值公式为:
Figure FDA0002663208610000053
将不同流速u1下的所述第一无量纲比值β1与所述流速传感器在的雷诺数Re拟合,得到不同流速u1下的第一拟合公式:
Figure FDA0002663208610000054
根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren
根据所述第二拟合公式,得到所述第三标定系数m3与n的数值;
通过计算公式:c=1-α计算所述光栅温度系数c,所述α为所述流速传感器的测风光栅和测温光栅温度系数比;
通过计算公式:Δλ=λ12计算所述测风光栅波长漂移量Δλ,所述λ1为所述流速传感器的测风光栅在流速u下对应的波长,所述λ2为所述流速传感器的温度补偿光栅在流速u下对应的波长。
4.一种流速测量装置,其特征在于,应用于流速传感器,所述装置包括:
测量单元,用于测得所述流速传感器产生的热量P以及测风光栅波长漂移量Δλ;
获取单元,用于获取固定参数,所述固定参数包括:所述流速传感器的光栅温度系数c、圆管特征长度d、流体导热系数k、所述流速传感器的探头面积S、运动粘度v、无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3以及定值n;
处理单元,用于根据所述流速传感器产生的热量P、所述测风光栅波长漂移量Δλ与所述固定参数,并结合流速标定公式计算得到所述流速传感器所放置的流场内的流速u,所述流速标定公式为:
Figure FDA0002663208610000061
所述无量纲比值β、第一标定系数m1、第二标定系数m2、第三标定系数m3、定值n通过以下述方法计算得到;
根据雷诺数公式计算所述流速传感器在不同流速u1下的雷诺数Re,所述雷诺数公式为:
Figure FDA0002663208610000062
计算所述流速传感器在不同流速u1下的换热系数h1与换热系数h2
根据公式:
Figure FDA0002663208610000063
计算所述换热系数h1,其中,ΔT为所述流速传感器的初始温度与在不同流速u1下热稳态状态时的温度之间的温度变化量;
根据公式
Figure FDA0002663208610000064
计算在不同流速u1下的努塞尔数Nu,其中,Pr为普朗特数;
根据公式
Figure FDA0002663208610000065
计算得到在不同风速u1所述换热系数h2
根据无量纲比值公式得到不同流速u1下的第一无量纲比值β1,所述无量纲比值公式为:
Figure FDA0002663208610000066
将不同流速u1下的所述第一无量纲比值β1与所述流速传感器在的雷诺数Re拟合,得到不同流速u1下的第一拟合公式:
Figure FDA0002663208610000067
根据所述无量纲比值公式与所述第一拟合公式得到所述第一标定系数m1与第二标定系数m2
将在不同流速u1下的所述雷诺数Re与在不同流速u1下的所述努塞尔数Nu拟合得到第二拟合公式:Nu=m3Ren
根据所述第二拟合公式,得到所述第三标定系数m3与n的数值。
CN201811506650.8A 2018-12-10 2018-12-10 一种标定方法、流速测量方法及装置 Active CN109342769B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811506650.8A CN109342769B (zh) 2018-12-10 2018-12-10 一种标定方法、流速测量方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811506650.8A CN109342769B (zh) 2018-12-10 2018-12-10 一种标定方法、流速测量方法及装置

Publications (2)

Publication Number Publication Date
CN109342769A CN109342769A (zh) 2019-02-15
CN109342769B true CN109342769B (zh) 2021-01-12

Family

ID=65304075

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811506650.8A Active CN109342769B (zh) 2018-12-10 2018-12-10 一种标定方法、流速测量方法及装置

Country Status (1)

Country Link
CN (1) CN109342769B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110398610B (zh) * 2019-08-29 2021-06-01 山东省科学院激光研究所 流速检测方法及光纤热线流速传感器探头
CN111610173B (zh) * 2020-05-27 2021-07-27 中国水利水电科学研究院 三维流体浓度场标定装置及标定方法
CN112380653B (zh) * 2020-11-17 2023-04-18 潍柴动力股份有限公司 换热器性能数据确定方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008113375A1 (de) * 2007-03-16 2008-09-25 Siemens Aktiengesellschaft Verfahren und vorrichtung zur messung von strömungsgeschwindigkeiten
CN104950133A (zh) * 2015-07-06 2015-09-30 浙江大学 微流体流速传感芯片、检测系统及检测方法
CN105699690A (zh) * 2016-01-28 2016-06-22 山东省科学院激光研究所 一种基于光纤光栅的涡街风速传感器及其组装方法
CN107462743A (zh) * 2017-08-08 2017-12-12 北京卫星环境工程研究所 适用于低气压下的风速标定装置及标定方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4678904A (en) * 1984-07-06 1987-07-07 Technology Dynamics, Inc. Optical measuring device using a spectral modulation sensor having an optically resonant structure
GB8520827D0 (en) * 1985-08-20 1985-09-25 York Ventures & Special Optica Fibre-optic sensing devices
US6813028B2 (en) * 2002-07-25 2004-11-02 Agilent Technologies, Inc. Calibration methodology and system for optical network analyzer

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008113375A1 (de) * 2007-03-16 2008-09-25 Siemens Aktiengesellschaft Verfahren und vorrichtung zur messung von strömungsgeschwindigkeiten
CN104950133A (zh) * 2015-07-06 2015-09-30 浙江大学 微流体流速传感芯片、检测系统及检测方法
CN105699690A (zh) * 2016-01-28 2016-06-22 山东省科学院激光研究所 一种基于光纤光栅的涡街风速传感器及其组装方法
CN107462743A (zh) * 2017-08-08 2017-12-12 北京卫星环境工程研究所 适用于低气压下的风速标定装置及标定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"A differentical pressure principle based mine wind velocity sensor using fiber bragg grating";王纪强 等;《16TH international conference on optical communication and networks》;20171231;全文 *
"High sensitivity hot-wire based wind velocity sensor using co-doped fiber and fiber bragg grating for use in mining applications";李振 等;《Journal of physics:IMEKO》;20181201;1-4页 *
"Novel negative pressure wave-based pipeline leak detection system using fiber bragg grating-based pressure sensors";王纪强 等;《journal of LT》;20171231;第35卷(第16期);3366-3371页 *
"一种海流传感器标定装置的研究与实验";陈迪志;《中国优秀硕士学位论文全文数据库》;20171231;I140-650 *

Also Published As

Publication number Publication date
CN109342769A (zh) 2019-02-15

Similar Documents

Publication Publication Date Title
CN109342769B (zh) 一种标定方法、流速测量方法及装置
CN107407590B (zh) 具有用于流体成分补偿的mems热式流量传感器
US7636640B2 (en) Multi-gas flow device
EP3032230B1 (en) Flow meter and a method of calibration
Mozley Predicting dynamics of concentric pipe heat exchangers
US9671266B2 (en) Method for thermally determining mass flow of a gaseous medium and thermal mass flow meter
CN103234662A (zh) 一种温度自动检测的补偿方法及温度自动检测系统
CN109100051B (zh) 温度传感器的动态响应的温度修正方法及装置
US20050109100A1 (en) Thermal mass flowmeter apparatus and method with temperature correction
CN109580033A (zh) 一种混凝土坝分布式光纤测温数据误差补偿方法
US20230258506A1 (en) Thermometer having a diagnostic function
CN105286812B (zh) 一种体温测量方法和装置
JP5695671B2 (ja) 流体の流速測定システム
Jamroz Relationship between dynamic coefficients of two temperature sensors under nonstationary flow conditions
CN104913862A (zh) 基于圆箔热流计的对流热流测量方法
RU2566641C2 (ru) Способ учета тепловой энергии, отдаваемой отопительным прибором
CN114577367B (zh) 一种光纤温度传感器标定的方法、装置和计算机设备
US20240068859A1 (en) Systems and methods for detecting a cryogenic fluid level in a container
RU2549256C1 (ru) Способ измерения параметров потоков жидкостей и газов
JP2018105789A (ja) 流量測定装置、設定装置、及び設定プログラム
JP4907959B2 (ja) フローセンサ用補正ユニット、流体判別装置、及び、流量計測装置
Schellekens et al. Design and results of a new interference refractometer based on a commercially available laserinterferometer
Simmons Continuous temperature‐logging equipment
SE449035B (sv) Sett och anordning for bestemning av en metkropps absoluta temperatur
CN114787599A (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
GR01 Patent grant
GR01 Patent grant