许文达,张 涛,毕 英,李 宝,刘伟光
(1. 天津大学电气与自动化工程学院,天津 300072;2. 国家蒸汽流量计量站,烟台 264003)
蒸汽对涡街流量计计量特性影响分析
许文达1,2,张 涛1,毕 英2,李 宝2,刘伟光1
(1. 天津大学电气与自动化工程学院,天津 300072;2. 国家蒸汽流量计量站,烟台 264003)
通过涡街流量计分别在以水为介质的水流量计量标准装置、以空气为介质的音速喷嘴法气体流量标准装置和以蒸汽为介质的蒸汽流量标准装置进行比对实验,测试了涡街流量计仪表系数和其介质的关系,实验结果表明 3种介质条件下仪表系数存在规律性偏差,利用流体力学及卡门涡街原理,结合蒸汽独特的热力学性质,深入分析了蒸汽对涡街流量计计量特性造成计量性能影响的因素,对涡街流量计在 3种介质下仪表系数存在规律性偏差实验结论进行了解释和分析,得到了蒸汽和空气流体条件下的计算拟合公式,并提出了系数修正的计算方法.
蒸汽;涡街流量计;可压缩流体;旋涡发生体;等熵指数;仪表系数
从20世纪80年代起,涡街流量计在流量计量中应用广泛,特别是近几年能源计量不断得到强化,大量被应用于蒸汽流量的计量,由于其压损小、量程比大等特点,部分取代了传统差压式蒸汽流量计.由于缺乏蒸汽实流检测的手段,蒸汽对涡街流量计计量性能的影响,一直缺乏相关的实验数据和理论分析,人们对其蒸汽测量的研究还比较匮乏.目前,国内外对涡街流量计在蒸汽测量上的溯源,普遍采用空气或者水为介质的流量计量标准装置代替蒸汽计量标准装置进行检定或校准,随着人们对涡街流量计的深入研究,不少学者对此种方法提出了质疑[1-2],并且在实际蒸汽流量的计量中也碰到过由于检定使用计量标准装置的不同,出现的蒸汽计量纠纷.针对这些问题,笔者在所研制的蒸汽流量计量标准装置的基础上,利用水、空气、蒸汽 3种不同介质流量计量标准装置,开展了涡街流量计的计量特性实验测试,从涡街流量计漩涡产生机理出发,对不同流体介质影响因素进行研究,定量分析了各种流体介质对涡街流量计的计量特性造成影响的因素.
涡街流量计是根据“卡门涡街”原理研制成的流体振荡式流量测量仪表.在一定雷诺数范围内旋涡的分离频率与旋涡发生体的几何尺寸、管道的几何尺寸有关,旋涡的频率正比于流量,其结构原理如图 1所示.
图1 卡门涡街传感器Fig.1 Sensor of Carmen vortex
涡街流量计所测量的是旋涡发生体两侧的平均速度1u,而反映被测流量的是管道中的平均流速u,它们之间的关系式为
式中:f为旋涡频率,Hz;Sr为斯特劳哈尔数;u1为发生体两侧的平均流速,m/s;d为发生体迎流面的宽度,m;u为被测介质来流的平均流速,m/s;m为旋涡发生体两侧弓形面积与管道横截面面积之比,不可压缩流体中,由于流体密度ρ不变,由连续性方程[3]可得到1/m u u= .
由此可得体积流量为
式中K为涡街流量计的仪表系数,1/m3.
从式(3)可见,涡街流量计的仪表系数仅与其机械结构和斯特劳哈尔数相关,与来流流量无关. 而K又表征着涡街流量计的计量特性,在下文中,重点分析在各种不同流体介质条件下K的变化规律.
蒸汽对涡街流量计计量特性的影响主要有如下3个因素:首先从式(3)中可以看出,影响K值大小的因素主要有机械结构尺寸 D、m、d和斯特劳哈尔数Sr 4个变量.从涡街流量计基本原理上分析,在不同流体介质条件下,机械结构尺寸的变化主要是由于温度的变化带来的热胀冷缩引起,其次不同流体可压缩性的差异也会导致仪表系数的偏差.此外Sr受雷诺数的影响,而雷诺数又受黏度的影响,流体的不同带来黏度的不同,根据相关的研究黏度对涡街流量计仪表系数的影响可以忽略[4-6].下文主要分析前 2个因素对蒸汽流量计的影响.
2.1 温度的影响
对图1所示发生体,得到m的计算公式为
对流体流通面积而言,可以把旋涡发生体看作为矩形(宽为d,长为D,见图1),面积比为
将式(5)带入式(3)中得
相关的研究[7-9]表明,当涡街流量计发生体为图1所示形状时,在 d/D=0.28时,产生的旋涡强度和旋涡的稳定性最好,故取d/D=0.28代入式(6)中得
由金属材料的线性膨胀公式得
式中:tD为壳体温度为t的直径;20D 为20,℃时壳体的直径;α为材料线性膨胀系数.
将式(8)带入式(7)中得到由于温度变化而引起的仪表系数的变化为
发生体和壳体为不锈钢材质(1Crl8Ni9Ti),在20~300,℃时线性膨胀系数α为,将α和不同温度下的Kt带入式(9)中得到数据如表1所示.
表1 温度变化引起仪表系数的变化Tab.1 Instrument coefficient changes with temperature change
2.2 介质的可压缩性引起的影响
本文是基于水、空气、蒸汽 3种不同流体介质下开展的实验研究,由式(3)得
由式(10)中可以看出,发生体两侧的介质流动速度1u与K值成正比,对于不可压缩流体由连续性方程可推导得到1u与管道平均流速u的关系[10]为
但对于可压缩流体式(11)不再成立,流动过程遵循可压缩流体的伯努利方程为
流动连续性方程为
等熵过程方程为
式中:κ为可压缩流体的等熵指数;1p和2p为仪表上游和发生体两侧处的压力,Pa;1ρ和2ρ分别为仪表上游和发生体两侧处的介质的密度,kg/m3.由式(12)~式(14)整理得
式(15)描述了可压缩流体u和1u的关系,与描述不可压缩流体的关系式(11)不同,三者的关系还和被测介质的等熵指数κ、压力1p、密度1ρ相关.被测介质的流量是由管道平均流速u所反映出来,而涡街流量计所检测到的旋涡所表示的流速是旋涡发生体两侧的平均流速1u,以下分析给出了由于流体的压缩性所引起的1u的差异性,即采用不可压缩流体水作为介质对涡街流量计进行检定,而实际应用于可压缩流体所引起的计量偏差.
由于 m为旋涡发生体两侧弓形面积与管道横截面面积之比,按照涡街流量的通用设计[11],对三角柱形发生体,在时漩涡发生的频率信号最强,取根据m的计算式(4)得
水为不可压缩流体,可按照式(11)计算得到1wu .空气和蒸汽为可压缩流体,由式(15)可得到工况下的u1a、u1e,与按照不可压缩流体式(11)得到的 u1比较,可计算分析得到由于介质的可压缩性带来的计量偏差.
实验条件下空气、蒸汽、水为测试流体介质时,各测试条件参数如表2所示.
表2 实验测试条件各流体介质参数Tab.2 Fluid’s parameters under experimental conditions
将上述参数带入式(15)中,采用 matlab对公式进行分析计算,得到1au 、1eu随u变化的数值,并通过式(11)得到1u随u变化的数值,设定管道平均速度u,得到的数值如表3所示.
表3的后2列分别表示了空气、蒸汽由于可压缩性而引入的计量偏差,由计算数据可以得出,相同管道平均流速u的条件下,可压缩介质发生体两侧的流速均比不可压缩流体的流速1u大,并且随u的增大,引起的差值增大,1u的相对偏差也越来越大,在流速为50,m/s时空气由于可压缩性可引起1.54%的偏差,蒸汽由于可压缩性可引起 0.68%的偏差,因此在涡街流量计的使用中需要考虑可压缩介质引起的影响.
上述的分析中,只能体现出特定工况条件下,气体的可压缩性对仪表系数的影响,虽然式(15)全面描述了所有工况的压缩性影响,但该式无法直观、简化显示,为了直观地体现式(15)中各变量带来对可压缩性的影响程度,基于不同条件下的计算数据对式(15)进行简化拟合,根据表2中空气工况下计算得到的u和1u的数据,拟合曲线,对比各种拟合方法,指数拟合的拟合度最高.拟合曲线的形式为
表3 不同条件下发生体两侧平均速度比较Tab.3 Comparison of average velocity of generating body on both sides under different conditions
式中A、b、0y分别为与工况条件p、ρ、κ相关的系数.对于其他工况下也用同样的方法拟合,拟合度均不低于 0.998.在拟合过程中得到不同工况下的系数A、b、0y,分别拟合这3个系数与不同工况条件p、ρ、κ的多元函数得
由拟合公式及各因子系数可以看出,κ的增大会引起b和A的增大,但b增大会引起1/uu变小,A的增大会引起1/uu的变大,因此气体等熵指数κ不会引起仪表系数单方向的变化.
11/p ρ的增大会引起b的增大和A的减小,b增大会引起1/u u变小,A的减小会引起1/u u的减小,因此11/p ρ与引起仪表系数的变化成反比.在空气及蒸汽 2种介质条件下,空气11/p ρ要小于蒸汽,因此空气介质受气体可压缩性的影响比蒸汽大.
2.3 各因素综合影响
设温度、可压缩性两因素对涡街流量计仪表系数K的影响因子分别为tF、Fκ,以水为介质的仪表系数Kw为基准,由式(9)和式(16)得空气和蒸汽为介质的仪表系数分别为
仪表系数偏差为
在表2所列的测试条件下,水中的仪表系数不受3个因素的影响可认为是定值,黏度的影响可忽略,计算在不同流速下,综合因素影响,引起的仪表系数偏差和及空气和蒸汽之间的偏差如表4所示.
表4 分析计算得到工况条件下不同介质仪表系数的偏差Tab.4 Analysis and calculation of deviation in different medium instrument coefficients in working condition
3.1 实验装置
本研究所采用的实验装置分别为冷凝称重法蒸汽流量计量标准装置、负压法音速喷嘴气体流量计量标准装置、质量法水流量计量标准装置.其中冷凝称重法蒸汽流量计量标准装置是国内唯一采用蒸汽实流标定的流量标准装置,扩展不确定度为 0.1%(k= 2),检定流量范围为 0.01~10,t/h,其工作原理如图 2所示,以过热蒸汽作为检定介质,将通过蒸汽流量仪表的蒸汽进行冷凝称重的原理,对蒸汽流量仪表进行实流标定.气体装置为负压法音速喷嘴气体流量计量标准装置,扩展不确定度为 0.25%(k=2),检定流量范围为 2.5~10,000,m3/h.水装置为质量法水流量计量标准装置,扩展不确定度为 0.05%(k=2),检定流量范围为0.05~200,m3/h.
3.2 实验方案
将不同厂家、型号、口径的涡街流量计,按照涡街流量计国家检定规程要求[12],分别安装在蒸汽装置、空气装置上和水装置上进行测试,其中蒸汽装置在温度为 150,℃左右、压力在 0.35,MPa的条件下测试,检定流量依次从最小流量到最大流量的 20%、40%、60%、80%、100%,共计6个流量点,每个流量点重复检定3次.目前已开展了近100台涡街流量计的测试实验,剔除不合格的涡街流量计,总结这些涡街流量计实验结果,仪表系数偏差的方向上基本相似,偏差的大小具有差异,在一定的范围内.本文中仅对具有代表性的DN50和DN100 2台涡街流量计实验数据进行分析,比较在不同流体条件下,不同流量点仪表系数变化的情况.
图2 冷凝称重法蒸汽流量标准装置工作原理Fig.2 Working principle of the steam flow standard device using condensing and weighing method
3.3 实验数据
通过3套不同流体介质下的实验测试,3套装置下的仪表系数、仪表系数偏差以及通过分析计算得到的仪表系数偏差如表5所示,仪表系数和流量的曲线关系如图3和图4所示.由实验数据可以看出,涡街流量计在空气、水、蒸汽流体条件下仪表系数依次减小,空气及蒸汽流体介质各点仪表系数随流量增大而增大,与理论计算分析相一致.通过积累的实验数据看,空气比水为流体介质测量的 K值偏大 0.2%~1.0%,比蒸汽为流体介质测量的K偏大1.5%~3.0%,仪表系数K值偏差的大小因仪表厂家和规格型号的不同而变化,但偏差的正负趋势基本一致.
表5 各口径涡街流量计特性汇总Tab.5 Metrology character summary of vortex flowmeter with different diameters
图3 DN50口径不同流量条件下3种流体仪表系数比较Fig.3 Comparison of three fluid instrument coefficients under different fluxes in DN50
图4 DN100口径不同流量条件下3种流体仪表系数比较Fig.4 Comparison of three fluid instrument coefficients under different fluxes in DN100
(1) 在高温蒸汽测量时,需要考虑温度导致的线性膨胀引起的影响,否则会导致仪表系数偏大.造成的计量误差,可按照文章中推导公式进行温度补偿.
(2) 气体的压缩性对涡街流量计仪表系数的影响,受介质流速u、压力p、密度ρ以及等熵指数κ等相关参数的影响,与流速u成正比,与11/p ρ的变化成反比,气体等熵指数κ不会引起仪表系数单方向的变化.
(3) 不同介质的可压缩性会对涡街流量计仪表系数产生影响,空气介质下仪表系数最大,水次之,蒸汽最小,在同样的流速下蒸汽介质对仪表系数的影响要小于空气的影响,具体差别大小根据工况条件下分析计算可得到.
(4) 根据实验数据,不同介质下测量的涡街流量计仪表系数与理论计算分析相一致,但2个结果还有一定的偏差,可能是由于介质的其他物性参数引起的对仪表系数的影响,还需深入的理论分析,在实际使用中尽量采用实流标定获得仪表系数.
[1] 郑灿亭,孙晓峰. 蒸汽计量中存在的两大难题[J]. 石油化工自动化,2007(2):85-88.
Zheng Canting,Sun Xiaofeng. Two difficult problems existed in steam measurement[J]. Automation in Etrochemical Industry,2007(2):85-88(in Chinese).
[2] 郑灿亭. 蒸汽计量存在的问题及对策[J]. 仪表与计量技术,2006(1):31-33.
Zheng Canting. Problems and solutions for steam measurement[J]. Instrument and Metrological Technology,2006(1):31-33(in Chinese).
[3] 苏彦勋,梁国伟,盛 健. 流量计量与测试[M]. 北京:中国计量出版社,2007.
Su Yanxun,Liang Guowei,Sheng Jian. Flow Measurement and Instrument [M]. Beijing:China Metrology Publishing House,2007(in Chinese).
[4] 郑丹丹,唐兴华,张 涛. 上游单弯头和闸阀对涡街流量计测量性能影响[J]. 天津大学学报,2011,44(12):1124-1129.
Zheng Dandan,Tang Xinghua,Zhang Tao. Impact of upstream single bend and gate valve on the measurement performance of vortex flowmeters[J]. Jourmal of Tianjin University,2011,44(12):1124-1129(in Chinese).
[5] Goujon-Durand S. Technical note,linearity of the vortex meter as a function of fluid viscosity[J]. Flow Measurement Instrumentation,1995,6(3):235-238.
[6] Hans V,Poppen G,von Lavante E,et al. Vortex shedding flowmeters and ultrasound detection:Signal processing and influence of bluff body geometry [J]. Flow Measurement and Instrumentation,1998,9(2):79-82.
[7] Von Lavante E,Perpeet S,Hans V,et al. Optimization of acoustic signals in a vortex-shedding flowmeter using numerical simulation[J]. International Journal of Heat and Fluid Flow,1999,20:402-404.
[8] 孙宏军,刘 畅. 垂直上升管道内气液两相细泡流涡街特性的数值仿真[J]. 天津大学学报,2011,44(6):540-546.
Sun Hongjun,Liu Chang. Numerical simulaton of characteristic of gas-liquid two-phase vortex street in vertically upward pipeline[J]. Jourmal of Tianjin University,2011,44(6):540-546(in Chinese).
[9] Hans V,Windorfer H. Comparison of pressure and ultrasound measurements in vortex flow meters[J]. Measurement,2003,33(2):121-133.
[10] 王 池,王自和,张宝珠,等. 流量测量技术全书[M]. 北京:化学工业出版社,2012.
Wang Chi,Wang Zihe,Zhang Baozhu,et al. Flow Measurement Technology Book [M]. Beijing:Chemical Industry Press,2012(in Chinese).
[11] 蔡武昌. 涡街流量计实现干法校准——旋涡发生体的标准化[J]. 世界仪表与自动化,2004(12):43-45.
Cai Wuchang. Vortex street flowmeter realizes dry calibration:Vortex generating body standardization [J]. International Instrumentation and Automation,2004(12):43-45(in Chinese).
[12] 国家质量监督检验检疫总局. JJG 1029—2007 涡街流量计检定规程[S]. 北京:中国标准出版社,2007.
General Administration of Quality Supervision,Inspection and Quarantine of PRC. JJG 1029—2007 Verification Regulation of Vortex Flowmeter[S]. Beijing:China Standard Press,2007(in Chinese).
(责任编辑:孙立华)
Analysis of the Influence of Steam on the Measurement Characteristics of Vortex Flowmeter
Xu Wenda1,2,Zhang Tao1,Bi Ying2,Li Bao2,Liu Weiguang1
(1. School of Electrical Engineering and Automation,Tianjin University,Tianjin 300072,China;2. National Steam Flowrate Measurement Station,Yantai 264003,China)
The relations between media and vortex flowmeter instrument coefficient have been tested in a comparative trial by using water as the medium of water flow measurement standard apparatus,air as the medium of sonic nozzle gas flow standard device and steam as the medium of steam flow standard device respectively. Experimental results show that there exists a regularity of deviation in the instrument coefficient under the condition of three kinds of media. Based on the fluid mechanics and the Carmen vortex principle,considering to the unique thermodynamic properties of vapor,the influencing factors on the characteristics of steam measurement of vortex flowmeter were analyzed. Moreover,the test results that instrument coefficient presented regularity deviation under three kinds of media were explained and analyzed. Calculation fitting formula under the condition of air and steam fluid is obtained and the calculating method of coefficient modification is proposed.
steam;vortex flowmeter;compressible fluid;vortex generating body;isentropic index;instrument coefficient
TH701
A
0493-2137(2014)08-0683-06
10.11784/tdxbz201303042
2013-03-20;
2013-05-10.
国家质量技术监督检验检疫总局科技资助项目(2010QK221).
许文达(1978— ),男,博士研究生,xuwenda2008@126.com.
张 涛,zt50@tju.edu.cn.