CN115077634B - 一种空气流量测量方法 - Google Patents

一种空气流量测量方法 Download PDF

Info

Publication number
CN115077634B
CN115077634B CN202210883457.6A CN202210883457A CN115077634B CN 115077634 B CN115077634 B CN 115077634B CN 202210883457 A CN202210883457 A CN 202210883457A CN 115077634 B CN115077634 B CN 115077634B
Authority
CN
China
Prior art keywords
value
gas flow
calculating
measurement
section
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
CN202210883457.6A
Other languages
English (en)
Other versions
CN115077634A (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.)
Hunan Minhang Automobile Technology Co ltd
Original Assignee
Hunan Minhang Automobile Technology 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 Hunan Minhang Automobile Technology Co ltd filed Critical Hunan Minhang Automobile Technology Co ltd
Priority to CN202210883457.6A priority Critical patent/CN115077634B/zh
Publication of CN115077634A publication Critical patent/CN115077634A/zh
Application granted granted Critical
Publication of CN115077634B publication Critical patent/CN115077634B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/05Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
    • G01F1/34Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明属于管道气流流量测量技术领域,具体涉及一种空气流量测量方法,测量精度高、抗干扰能力强,受外界环境的影响因素小,连续采样获得的数据在预估值与观测值之间进行方差计算,基于数据测量的高斯分布,建立防脉冲干扰平均滤波法递归模型,使测量结果更加接近真实值;并且可以适用于非稳定工况及稳定工况的数据测量,尤其是在非稳态测量中出现数据波动较大的情况下确保测量的精确度。

Description

一种空气流量测量方法
技术领域
本发明属于管道气流流量测量技术领域,具体涉及一种空气流量测量方法。
背景技术
现有技术下,管道气流流量的测量往往直接采用热膜式、热线式和超声波式的空气流量计,这些流量计容易受到受灰尘、空气湿度和温度、传感器老化和传感器线性偏差的影响,抗干扰能力较弱,测量精度和重复一致性较差,并且在非稳定工况下适用性较差。
发明内容
发明要解决的技术问题在于克服现有技术中管道气流测量抗干扰能力弱、测量精度和重复一致性差并且在非稳定工况下适用性较差的缺陷,从而提供一种空气流量测量装置和方法。
一种空气流量测量方法,包括以下步骤:
步骤S1:根据气体的温度和湿度计算气体的密度;
步骤S2:根据待测管段第一分段和第二分段的两段的截面积和气压值,计算第二分段的气体流速计算值,并计算第二分段的气体流量估算值;
步骤S3:根据第二分段的气体流速测量值和第二分段的截面积,计算第二分段的气体流量测量值;
步骤S4:根据第一分段的和第二分段的两段的气压值的多个监测样本,分别计算气体流量估算值,并计算气体流量估算值的算术平均值、方差和标准差;
步骤S5:同理,根据气体流速测量值的多个监测样本,分别计算气体流量测量值,并计算气体流量测量值的算术平均值、方差和标准差;
步骤S6:根据气体流量估算值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量估算值,并计算气体流量估算值在置信区间内的算术平均值、方差和标准差;
步骤S7:同理,根据气体流量测量值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量测量值,并计算气体流量测量值在置信区间内的算术平均值、方差和标准差;
步骤S8:建立防脉冲干扰平均滤波法递归模型,引入迭代增益,并根据迭代增益、气体流量测量值的算术平均值和估算值的算术平均值,递归修正新的气体流量估算值的算术平均值。
进一步的,所述步骤S2:根据待测管段第一分段和第二分段的两段的截面积和气压值,计算第二分段的气体流速计算值,并计算第二分段的气体流量估算值,包括:
根据连续性方程和伯努利方程: 其中第一分段和第二分段的截面积分别记为A1、A2,气体流速分别记为V1、V2,气压值分别记为P1、P2,中心线高度分别记为h1、h2,气体密度记为ρ;
可得气体流速计算值其中Δp为两个分段的压力差,Δh为两个分段的高度差;
计算气体流量估算值Q=A2*V2
进一步的,所述步骤S4:根据第一分段的和第二分段的两段的气压值的多个监测样本,分别计算气体流量估算值,并计算气体流量估算值的算术平均值、方差和标准差,包括:
计算气体流量估算值的算术平均值
计算气体流量估算值的方差
计算气体流量估算值的标准差
其中Q1表示计算气体流量估算值的第一个样本,Qm表示计算气体流量估算值的第m个样本,m表示样本数量。
进一步的,所述步骤S6:根据气体流量估算值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量估算值,并计算气体流量估算值在置信区间内的算术平均值、方差和标准差,包括:
建立置信区间其中μ为气体流量估算值的算术平均值/>σ为气体流量估算值的标准差SQ,m代表样本数量,c为预设的标准分数,L为与标准分数对应的置信水平;
计算置信区间内样本流量估算值的算术平均值
计算置信区间内样本流量估算值的估算值方差
计算置信区间内样本流量估算值的估算值标准差
进一步的,所述步骤S8:建立防脉冲干扰平均滤波法递归模型,引入迭代增益,并根据迭代增益、气体流量测量值的算术平均值和估算值的算术平均值,递归修正新的气体流量估算值的算术平均值包括:
根据置信区间内样本流量估算值的估算值方差和测量值方差/>取值滤波增益值F,并且/>
递归修正新的气体流量估算值的算术平均值
进一步的,还包括以下步骤:根据测得的气体流量,标定空气流量计,并通过对比分析,得出空气流量计的精度与重复一致性。
进一步的,所述样本数量m取值为100~500个。
进一步的,所述标准分数c的取值为1.96,并对应置信水平L为95%。
一种空气流量测量装置,包括:
大气环境站,所述大气环境站安装在待测管段的入口处,用于监测入口气流的温度和湿度;
风速仪,所述风速仪安装在待测管段的出口处,用于监测出口气流的风速;
微压计,两个所述微压计分别连接待测管段两段不同内径的分段,用于监测两段不同内径的分段的气压;
数采设备前端,用于将大气环境站和风速仪采集的模拟量数据转换为数字量数据;
上位机,用于从数采设备前端获取转换的数字量数据,从微压计获取气压数据,并分析计算待测管道的气体流量。
进一步的,还包括相机,两台所述相机分别拍摄两个所述微压计,并将拍摄所得照片发送至所述上位机,所述上位机对所述照片进行图像识别获得所述微压计气压值读数。
有益效果:
1.本发明提供的一种空气流量测量方法测量精度高、抗干扰能力强,受外界环境的影响因素小,连续采样获得的数据在预估值与观测值之间进行方差计算,基于数据测量的高斯分布,建立防脉冲干扰平均滤波法递归模型,使测量结果更加接近真实值;并且可以适用于非稳定工况及稳定工况的数据测量,尤其是在非稳态测量中出现数据波动较大的情况下确保测量的精确度。
2.本发明提供的一种空气流量测量方法还可以对其他类型的流量测量装置进行标定,分析测量偏差和测量不确定度。
3.本发明提供的一种空气流量测量装置,通过微压计、风速仪、大气环境站、微压计、数采设备前端和上位机的配合设置,能够实现不同工况下进气流量标定,结构紧凑,抗干扰能力强,易于实验室和检测机构使用。相比于热膜式、热线式和超声波式的空气流量计,该流量计及标定装置不会受灰尘、空气湿度和温度、传感器老化和传感器线性偏差的影响,可最大程度降低测量误差,提高测量精度和重复一致性;能够精确计算管道中气体的流量,为后端用气设备或空气流量计提供准确的流量结果,同时能够适应不同管道的气流状态,能够计算得出非稳态工况和稳态工况下的气体流量。
4.本发明提供的一种空气流量测量装置,通过相机拍摄微压计,上位机对所述照片进行图像识别获得所述微压计气压值读数,在确保获得微压计高精度数据的同时,能够确保数据有效转换为数字信号,方便了后续计算,提高了流量测量效率。
附图说明
图1为本发明的装置整体结构示意图;
图2为本发明的方法流程示意图。
附图标记:1、待测管段;2、大气环境站;3、风速仪;4、微压计;5、数采设备前端;6、上位机;7、相机。
具体实施方式
实施例1
参照图1所示,本实施例提供了一种空气流量测量装置,包括:
大气环境站2,所述大气环境站2安装在待测管段1的入口处,用于监测入口气流的温度和湿度;
风速仪3,所述风速仪3安装在待测管段1的出口处,用于监测出口气流的风速;
微压计4,两个所述微压计4分别连接待测管段1两段不同内径的分段,用于监测两段不同内径的分段的气压;两段不同内径的分段设置开口,并通过橡胶软管与所述微压计4相连;所述微压计4为U型微压计最小刻度为1mm,确保气压的精确监测。
数采设备前端5,用于将大气环境站2和风速仪3采集的模拟量数据转换为数字量数据;
上位机6,用于从数采设备前端5获取转换的数字量数据,从微压计4获取气压数据,并分析计算待测管道的气体流量;所述上位机6为高性能的计算机,数采设备前端5通过CAN总线与上位机6连接,确保数据高速传输以及有效快速分析计算。
还包括相机7,两台所述相机7分别拍摄两个所述微压计4,并将拍摄所得照片发送至所述上位机6,所述上位机6对所述照片进行图像识别获得所述微压计4气压值读数;所述相机7为彩色相机,通过网线连接所述上位机6并传输拍摄的照片,确保微压计4的读数被有效读取。
所述相机7拍摄微压计4读数的频率为10Hz,风速仪3和大气环境站2通过数采设备前端5的采样频率为10Hz。
选择的待测管段1为平直的内壁光滑的金属管段,包括两段不同内径的分段,并且所述管段的中心线水平;避免因为管道的设置影响流量测量的准确性。
实施例2
参照图2所示,本实施例提供了一种空气流量测量方法,包括以下步骤:
步骤S1:根据气体的温度和湿度计算气体的密度
空气的密度ρ∈(T,Rh),通过查表法获得;
步骤S2:根据待测管段第一分段和第二分段的两段的截面积和气压值,计算第二分段的气体流速计算值,并计算第二分段的气体流量估算值:
根据连续性方程和伯努利方程: 其中第一分段和第二分段的截面积分别记为A1、A2,气体流速分别记为V1、V2,气压值分别记为P1、P2,中心线高度分别记为h1、h2,气体密度记为ρ;
可得气体流速计算值其中Δp为两个分段的压力差,所述压力差与大气压无关,是两个分段气压的相对值;Δh为两个分段的高度差,预设两个分段中心线为同一高度,因此Δh=0;
计算气体流量估算值Q=A2*V2
步骤S3:根据第二分段的气体流速测量值V2z和第二分段的截面积,计算第二分段的气体流量测量值:QZ=A2*V2z
步骤S4:根据第一分段的和第二分段的两段的气压值的多个监测样本,分别计算气体流量估算值,并计算气体流量估算值的算术平均值、方差和标准差;
采集20s的试验数据,采样频率为10Hz,样本数量为200;
计算气体流量估算值的算术平均值
计算气体流量估算值的方差
计算气体流量估算值的标准差
其中Q1表示计算气体流量估算值的第一个样本,Qm表示计算气体流量估算值的第m个样本,m表示样本数量。
步骤S5:同理,根据气体流速测量值的多个监测样本,分别计算气体流量测量值,并计算气体流量测量值的算术平均值、方差/>和标准差SQZ
步骤S6:根据气体流量估算值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量估算值,并计算气体流量估算值在置信区间内的算术平均值、方差和标准差;
建立置信区间其中μ为气体流量估算值的算术平均值/>σ为气体流量估算值的标准差SQ,m代表样本数量,c=1.96为预设的标准分数,L=95%为与标准分数对应的置信水平;
计算置信区间内样本流量估算值的算术平均值
计算置信区间内样本流量估算值的估算值方差
计算置信区间内样本流量估算值的估算值标准差
步骤S7:同理,根据气体流量测量值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量测量值,并计算气体流量测量值在置信区间内的算术平均值、方差和标准差;
步骤S8:建立防脉冲干扰平均滤波法递归模型,引入迭代增益F,并根据迭代增益、气体流量测量值的算术平均值和估算值的算术平均值,递归修正新的气体流量估算值的算术平均值。
根据置信区间内样本流量估算值的估算值方差和测量值方差/>取值滤波增益值F,并且/>
递归修正新的气体流量估算值的算术平均值
能够看出,如果估值太小,新的估计值会加上一个量/>如果/>估值太大,新的估计值会减去一个量/>从而确保新的估计值比现有值准确,多次递归修正气体流量估算值的算术平均值就会越来越准确。
并且防脉冲干扰平均滤波法计算的数值更倾向于测量或者估算结果更稳定一组的方差。同时,由于是递归估算结果,适应于稳态工况与非稳态工况。非稳态工况下的修正会依赖于估算值与/>在样本数量下的方差情况,方差越小的将会更加接近真实测量结果。
还包括以下步骤:根据测得的气体流量,标定空气流量计,并通过对比分析,得出空气流量计的精度与重复一致性。
将两种测量的流量值进行公示迭代计算出稳定工况下的流量结果,可以精确的标定热膜式、热线式和超声波式的空气流量计,并且通过对比分析,能够得出上述仪器测量的精度与重复一致性;
对于非稳态工况,管道内气流组织复杂,压差法的计算将会呈现较大的波动,95%置信区间内的数据将容易将脉冲情况下的数据剔除,保证数值选取的合理性,通过滤波迭代计算,测量结果更加趋近与测量值稳定的一方;并且由于采用迭代滤波计算,任何时刻的采集数据结果将是此刻前20秒的数据进行置信区间计算,提取2σ区间的数据进一步进行均值、方差和标准差计算,最后通过测量值与压差估算值之间的对比递归,获得流量的实验结果。
这种测量方法测量精度高、抗干扰能力强,受外界环境的影响因素小,连续采样获得的数据在预估值与观测值之间进行方差计算,基于数据测量的高斯分布,建立防脉冲干扰平均滤波法递归模型,使测量结果更加接近真实值;并且可以适用于非稳定工况及稳定工况的数据测量,尤其是在非稳态测量中出现数据波动较大的情况下确保测量的精确度。

Claims (4)

1.一种空气流量测量方法,其特征在于,包括以下步骤:
步骤S1:根据气体的温度和湿度计算气体的密度;
步骤S2:根据待测管段第一分段和第二分段的两段的截面积和气压值,计算第二分段的气体流速计算值,并计算第二分段的气体流量估算值;
步骤S3:根据第二分段的气体流速测量值和第二分段的截面积,计算第二分段的气体流量测量值;
步骤S4:根据第一分段的和第二分段的两段的气压值的多个监测样本,分别计算气体流量估算值,并计算气体流量估算值的算术平均值、方差和标准差;
步骤S5:同理,根据气体流速测量值的多个监测样本,分别计算气体流量测量值,并计算气体流量测量值的算术平均值、方差和标准差;
步骤S6:根据气体流量估算值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量估算值,并计算气体流量估算值在置信区间内的算术平均值、方差和标准差;
步骤S7:同理,根据气体流量测量值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量测量值,并计算气体流量测量值在置信区间内的算术平均值、方差和标准差;
步骤S8:建立防脉冲干扰平均滤波法递归模型,引入迭代增益,并根据迭代增益、气体流量测量值的算术平均值和估算值的算术平均值,递归修正新的气体流量估算值的算术平均值;
所述步骤S2:根据待测管段第一分段和第二分段的两段的截面积和气压值,计算第二分段的气体流速计算值,并计算第二分段的气体流量估算值,包括:
根据连续性方程和伯努利方程:A1*V1=A2*V2, 其中第一分段和第二分段的截面积分别记为A1、A2,气体流速分别记为V1、V2,气压值分别记为P1、P2,中心线高度分别记为h1、h2,气体密度记为ρ;
可得气体流速计算值其中Δp为两个分段的压力差,Δh为两个分段的高度差;
计算气体流量估算值Q=A2*V2
所述步骤S4:根据第一分段的和第二分段的两段的气压值的多个监测样本,分别计算气体流量估算值,并计算气体流量估算值的算术平均值、方差和标准差,包括:
计算气体流量估算值的算术平均值
计算气体流量估算值的方差
计算气体流量估算值的标准差
其中Q1表示计算气体流量估算值的第一个样本,Qm表示计算气体流量估算值的第m个样本,m表示样本数量;
所述步骤S6:根据气体流量估算值的算术平均值和标准差,通过2σ区间法,去除置信区间以外的气体流量估算值,并计算气体流量估算值在置信区间内的算术平均值、方差和标准差,包括:
建立置信区间其中μ为气体流量估算值的算术平均值/>σ为气体流量估算值的标准差SQ,m代表样本数量,c为预设的标准分数,L为与标准分数对应的置信水平;
计算置信区间内样本流量估算值的算术平均值
计算置信区间内样本流量估算值的估算值方差
计算置信区间内样本流量估算值的估算值标准差
所述步骤S8:建立防脉冲干扰平均滤波法递归模型,引入迭代增益,并根据迭代增益、气体流量测量值的算术平均值和估算值的算术平均值,递归修正新的气体流量估算值的算术平均值包括:
根据置信区间内样本流量估算值的估算值方差和测量值方差/>取值滤波增益值F,并且/>
递归修正新的气体流量估算值的算术平均值
2.根据权利要求1所述的一种空气流量测量方法,其特征在于,还包括以下步骤:根据测得的气体流量,标定空气流量计,并通过对比分析,得出空气流量计的精度与重复一致性。
3.根据权利要求1所述的一种空气流量测量方法,其特征在于,所述样本数量m取值为100~500个。
4.根据权利要求1所述的一种空气流量测量方法,其特征在于,所述标准分数c的取值为1.96,并对应置信水平L为95%。
CN202210883457.6A 2022-07-26 2022-07-26 一种空气流量测量方法 Active CN115077634B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210883457.6A CN115077634B (zh) 2022-07-26 2022-07-26 一种空气流量测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210883457.6A CN115077634B (zh) 2022-07-26 2022-07-26 一种空气流量测量方法

Publications (2)

Publication Number Publication Date
CN115077634A CN115077634A (zh) 2022-09-20
CN115077634B true CN115077634B (zh) 2023-11-03

Family

ID=83243283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210883457.6A Active CN115077634B (zh) 2022-07-26 2022-07-26 一种空气流量测量方法

Country Status (1)

Country Link
CN (1) CN115077634B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104034372A (zh) * 2014-06-12 2014-09-10 杭州冠一流体技术有限公司 一种密度修正型多参数瓦斯流量计及密度修正方法
CN108896438A (zh) * 2018-07-25 2018-11-27 重庆大学 孔底瓦斯含量测量装置和测量方法
WO2020005149A2 (en) * 2018-06-29 2020-01-02 Provtagaren Ab Method for digital flow measurement in pulsating flows
CN111426353A (zh) * 2020-04-08 2020-07-17 中国民用航空飞行学院 一种精确流量获取装置及方法
CN114298091A (zh) * 2021-12-13 2022-04-08 国网湖北省电力有限公司电力科学研究院 Sf6气体流量计量值修正方法、装置、设备及存储介质

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103939749B (zh) * 2014-04-28 2016-08-24 东北大学 基于大数据的输油管网泄漏智能自适应监控系统及方法
GB201700428D0 (en) * 2017-01-10 2017-02-22 Able Instr & Controls Ltd Apparatus and method for flare flow measurement

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104034372A (zh) * 2014-06-12 2014-09-10 杭州冠一流体技术有限公司 一种密度修正型多参数瓦斯流量计及密度修正方法
WO2020005149A2 (en) * 2018-06-29 2020-01-02 Provtagaren Ab Method for digital flow measurement in pulsating flows
CN108896438A (zh) * 2018-07-25 2018-11-27 重庆大学 孔底瓦斯含量测量装置和测量方法
CN111426353A (zh) * 2020-04-08 2020-07-17 中国民用航空飞行学院 一种精确流量获取装置及方法
CN114298091A (zh) * 2021-12-13 2022-04-08 国网湖北省电力有限公司电力科学研究院 Sf6气体流量计量值修正方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN115077634A (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
US5827963A (en) System and method for determining a density of a fluid
CN113916486B (zh) 一种基于多参数压力回归算法低速风洞速压测控方法
CN110081944B (zh) 一种基于实时压力变化的气体测量方法及所用装置
CN114964429A (zh) 一种超声波水表检定方法
CN102590557A (zh) 一种可变径负压式微风速标定装置
CN102538914A (zh) 一种具有校正功能的电子显示水表
CN109253781A (zh) 节流式压差流量传感器的校准方法及校准系统
CN114923547A (zh) 一种燃气表示值误差不确定度的自动评定装置及方法
CN115629027A (zh) 基于形变数据补偿处理的颗粒物检测装置和方法
CN115077634B (zh) 一种空气流量测量方法
CN211576248U (zh) 通风与空调专用风量风压测量测片
CN112229454A (zh) 测量烟用吸阻标准件流量压差曲线的系统和方法
CN110361077A (zh) 一种流量及风速传感器标校装置及标校方法
CN111272236B (zh) 气体层流流量计气体流量的计算方法及气体层流流量计
CN115389090A (zh) 基于参考端控制的压力测量系统及测量方法
CN211576247U (zh) 槽口型风量风压测量测片
CN210135981U (zh) 一种流量及风速传感器标校装置
RU86751U1 (ru) Измерительная аэродинамическая установка
CN1172164C (zh) 一种测量含尘气流流量的方法及其在线监测装置
CN210426717U (zh) 一种音速喷嘴检定装置
CN114046861A (zh) 一种提高差压式流量计动态流量测量准确度的方法
CN116202596B (zh) 伽马空管计数实时修正方法
CN219104080U (zh) 一种具有组分补偿功能的热式气体流量计及其系统
JP2000221065A (ja) 可変ベンチュリーの流量校正方法
CN112729433B (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