多介质和多类型气体流量计数学模型的应用

吴富伟，董文润，黄益昌

(1.上海理工大学光电信息与计算机工程学院，上海 200093；2.天信仪表集团有限公司，浙江 温州 325800)

0 引言

气体流量计量是气体输送工程中重要的测量环节，它是企业进行贸易交接、经济分析以及降低运行成本的主要因素，直接影响企业的经济效益与用户利益[1]。气体流量计是一种用于计量气体体积的装置，在民用、商业和工业场合都得到了广泛的应用。根据被测管道气体的预期流量范围、被测气体的类型以及环境因素等，需要选用不同类型的气体流量计。

目前，用于气体贸易计量的流量计按工作原理分类，主要有速度式、容积式和差压式等类型。速度式流量计包括涡轮、涡街、超声、旋进旋涡等流量计，容积式流量计包括腰轮流量计等，一般具备输出脉冲信号与二次装置(流量显示仪表)连接；差压式流量计或一次装置，如孔板、V锥、喷嘴和文丘里管等，则不同于上述输出脉冲信号的方式，通常是输出模拟信号与二次装置连接。

由于当前计算机技术的迅速发展，流量显示仪表(下称体积修正仪)所采用的单片机运算能力越来越强大。根据目前的技术现状，体积修正仪可以同时实现多种介质、多种输入信号(脉冲信号、模拟信号等)的运算处理，即同一台体积修正仪能适用于不同介质、不同工作原理的流量装置。因此，有必要对多种流量计的数学模型进行研究。

1 脉冲式流量计数学模型

脉冲输出式流量计的气体流速通常与脉冲频率成正比。按照信号发生方式的不同，信号分为高频脉冲和低频脉冲两种信号形式。涡轮流量计通过计算转子的转速来测量流量。叶片与气流速度成正比[2]，涡轮流量计的高频脉冲信号通常取自叶轮或叶轮轴上的光栅盘，腰轮流量计的高频脉冲信号取自腰轮转子，则操作条件下的瞬时流量为：

(1)

式中：qv为操作条件下的瞬时流量，m3/h；f为脉冲频率，Hz；K为仪表系数，m-3。

仪表系数K除了与涡轮导程、叶片数、叶片宽度和螺旋升角等结构因素有关外，还与流体黏性、轴承阻尼和润滑油黏度等有关[3]。低频脉冲信号通过一套机械传动机构及磁耦合连接件传送到机械式计数器，操作条件下的瞬时流量为：

(2)

式中：n为间隔时间内测得的脉冲数，无量纲；tn为n个脉冲所占用的时间，s；K为仪表系数，m-3。

贸易结算时，需要将操作条件下的气体体积转换为标准状态下的气体体积。天然气在标准状态下的瞬时流量为：

(3)

式中：qn为标准条件下的瞬时流量，m3/h；Zn、Z分别为操作和标准条件下的压缩因子、无量纲；p、T分别为操作条件下的气体压力和温度，kPa、K；pn、Tn分别为标准条件下的气体压力(101.325 kPa)和温度(293.15 K)。

空气、氧气和氮气等湿气体在标准状态下的瞬时流量为：

(4)

式中：φ、φn分别为操作和标准条件下的相对湿度，%；psmax为温度为T时水蒸气的最大压力，kPa；pnsmax为温度为Tn时水蒸气的最大压力，kPa；

2 差压式流量计数学模型

差压式流量计是基于流体流动的节流原理，利用流体流经节流装置时产生的上下游压力差实现流量测量的。孔板流量计的节流装置由标准孔板、取压装置和上下游直管段组成。标准孔板是一块由机械加工形成的圆形穿孔薄板，它的节流孔圆柱面与孔板上游端面垂直，边缘尖锐[3-4]。其原理是气体通过节流孔时产生压降，通过测量压力差来计算气体流量。由于没有运动部件，标准孔板性能稳定可靠且使用期限长[5]。对于孔板、喷嘴和文丘里管节流装置，气体在操作条件下的体积流量计算公式为：

(5)

式中：C为流出系数，即通过节流装置的实际流量值与理论流量值的比值，无量纲；β为节流孔或喉部的直径与上游测量管内径的比值，无量纲；ε为可膨胀系数，用于修正当气体通过节流件时密度的变化引起的流量变化，无量纲；d为节流孔或喉部直径，m；Δp为在上下游取压孔位置上取得的静压力之差，Pa；ρ1为操作条件下气体在上游取压孔处的密度，kg/m3。

V锥流量计是基于文丘里管测量原理，以一个同轴安装在测量管内的尖圆锥体作为节流件的差压式流量测量装置。它将节流布局从中心孔节流改为环状节流，在流体流经内锥体时会在其周围形成边界层，并疏导流体离开锥体尾部边缘，从而减少节流件磨损。对于V锥节流装置，气体在操作条件下的体积流量计算公式为：

(6)

接下来分析各种不同类型节流装置的流出系数C、可膨胀系数ε和气体密度ρ1的求解过程。C为通过节流装置的实际流量与理论流量之间的系数。利用不可压缩流体对节流装置的校准结果表明，在一定的安装条件下，C仅与雷诺数有关[6]。只要节流装置几何相似，在雷诺数相同的情况下C是相等的。孔板流量计C的计算公式为：

(7)

(8)

当D<71.12 mm时，C=C1+C2；当D≥71.12 mm时，C=C1。

式中：ReD为管径雷诺数，无量纲；L1为孔板上游端面到上游取压孔轴线的距离与测量管内径的比值，无量纲；L2为孔板下游端面到下游取压孔轴线的距离与测量管内径的比值，无量纲。

可膨胀系数是考虑到流体的可压缩性而使用的系数，其值取决于压力比和等熵指数。将可膨胀系数应用于理论方程，以调节气体在通过孔板时密度的减小引起的流量变化[7]。孔板流量计的可膨胀系数按下式计算：

(9)

式中：τ为孔板下游气流与上游气流绝对静压力的比值，无量纲；к为等熵指数，无量纲。

喷嘴节流装置根据内部结构的不同，可分为ISA1932喷嘴、长径喷嘴和文丘里喷嘴。喷嘴在管道内的部分是圆形的，由圆弧廓形的收缩部分和圆筒形喉部组成[8]。ISA1932喷嘴的流出系数由下式计算：

C=0.99-0.226 2β4.1-(0.001 75β2-

(10)

长径喷嘴由四分之一椭圆状收缩入口和圆筒喉部组成，其圆筒喉部比ISA1932喷嘴要长，长度为喉部直径的0.6。其流出系数计算式为：

(11)

文丘里喷嘴与以上两种喷嘴有很大的区别。它的上游部分结构与ISA1932喷嘴相同，但它增加了一个扩散段，这种结构的压力损失会更低。其流出系数按下式计算：

C=0.985 8-0.196β4.5

(12)

喷嘴流量计的可膨胀系数按下式计算：

(13)

文丘里管由入口圆筒段、圆锥收缩段、圆筒喉部和圆锥扩散段依次连接而成[8-9]，它与V锥流量计的流出系数直接由厂家标定。V锥流量计的可膨胀系数计算公式为[10]：

(14)

利用干空气的物质的量、干空气在标准条件下的压缩因子和天然气相对密度，可推导出操作条件下天然气在上游取压孔处的密度为：

(15)

式中：p1为操作条件下天然气在上游取压孔处的静压力，MPa；Gr为天然气相对密度，无量纲。

当压缩因子作为参数输入或已知天然气详细组分物质的量分数的情况下，天然气相对密度由下式计算：

(16)

式中：xj为天然气各组分的物质的量分数，无量纲；Gij为天然气各组分的理想相对密度，无量纲。

空气、氮气和氧气等湿性气体的密度应该由湿气体的干气部分和蒸汽部分相加得到：

(17)

式中：ρn为标准条件下的气体密度，kg/m3；ρsmax为操作条件下蒸汽的最大密度，kg/m3。

蒸汽根据状态可分为饱和蒸汽和过热蒸汽。当水分子的蒸发与凝结处于动平衡状态时，空间中蒸汽分子的密度保持不变，此时的状态称为饱和状态。当蒸汽达到饱和状态后如果继续加热则温度会上升，成为过热蒸汽。饱和蒸汽的压力和温度呈线性关系，因此只要测得温度或压力，即可通过查表得到蒸汽密度；而过热蒸汽则同时需要测得压力和温度，通过查表获得蒸汽密度。

天然气在标准条件下的体积流量可由式(3)计算得到，氮气、氧气等湿气体在标准条件下的体积流量可由式(4)计算得到，蒸汽的质量流量和热量流量分别由操作条件下的体积流量乘以蒸汽密度和热值得到，蒸汽热值与密度通过查表得到。

3 数学模型的计算处理

脉冲输出式流量计高频和低频脉冲信号的单片机处理方式略有不同。对于高频信号，只需测得间隔时间内的脉冲数即可计算流量值，而低频脉冲信号需要精确测量间隔时间内测得的脉冲数所占用的时间。因此，利用定时器计算32.768 kHz晶振输出信号的累计时间，使用外部中断测量低频脉冲数，从而求出流量值。两个脉冲信号的间隔时间测量流程如图1所示。

图1 间隔时间测量流程图

差压式流量计的数学模型计算公式有很多的幂运算。这样的运算对运行频率相对较低的单片机来说，运算负荷是相当大的。由于在流量计的计量现场一般没有外电源，体积修正仪的内置锂电池模块通常需要连续工作几年以上，因此需要尽可能快速地完成数学模型计算工作，以进入休眠状态，从而降低功耗。对于一些数学运算，在不影响精度的情况下，可以采用快速的计算方法。平方根的运算可以运用牛顿迭代法来近似求解，算法如下[11]。

float InvSqrt(float x)

{

float xhalf = 0.5f * x;

int i = *(int*)&x;

i = 0x5f375a86 - (i >> 1);

x = *(float*)&i;

x = x * (1.5f - xhalf * x * x);

return x;

}

如需更高的精度，则可以重复运行第7行的牛顿迭代式。运算量较大的公式可以根据变量的取值范围，预先建立数据表格的方法。使用时，根据输入变量的值查询表格中的数据并进行线性插值计算。在建立数据表格时，需要对线性插值的计算值与公式实际计算值进行误差比较，适当调整变量步长以保证计算精度。

如C只与β和ReD有关，故可根据β和ReD的取值范围建立数据表格，需要时再通过查询值进行精确插值计算。线性插值计算图如图2所示。

图2 线性插值计算图

为了计算C，需要先查询表格中相邻四点的流出系数值C1、C2、C3和C4，再利用线性插值计算C：

(18)

查表计算流程如图3所示。

图3 查表计算流程图

4 结束语

本文通过对不同工作原理的流量计用于多种介质的数学模型进行分析及研究，使同一台体积修正仪能与不同工作原理流量计配套使用。速度式、容积式和差压式使用同一台体积修正仪，并实现了多种介质的准确测量。体积修正仪实现了通用气体流量计量的功能。