气体压缩因子Z的在线修正与单片机实现

    一、 简述
    通常人们认为气体分子是弹性的、不占体积的质点，分子间除了相互碰撞瞬间外没有任何作用力，是理想气体，遵循状态方程PV=RT。实际上符合这些条件的理想气体是不存在的。实验表明，只有在低压、高温下，实际气体才可以近似看作理想气体。

    由于实际气体与理想气体的差异，使得在气体流量计量中的度和可靠性难以评价，特别是西气东送这样的低温、高压管道流量的计量，所以实际气体不能用理想气体状态方程来进行处理。

    例如：已知容器体积为0.030m3，里面有质量为3.7kg、温度为243.2K的CO，已知压力的实测值为7.09?105Pa，用理想气体状态方程计算摩尔体积为：

    Vm=2.625m3/kmol，

    摩尔体积常数为：

    Rm=8.314?103 J/（kmol.K）

    其压力为：

    P=  =7.71?105  Pa[1]，

    压力的计算值与实测值的偏差为：d=8.74%。

    可见，当压力增大时，用理想气体状态方程计算出来的结果与实际结果已有较大的偏离。

    二、 压缩因子Z定义
    
    在高压低温下，任何气体对理想状态方程都会出现明显的偏差，而且压力愈高，温度愈低，偏离就愈多。实际气体的这种偏离，采用PV与RT的比值来说明，这个比值称为压缩因子，以符号Z表示，定义为：  = 

    式中，V—实际气体在（P，T）下的比容，

    V0 —代表在相同的（P，T）状态下理想气体的比容。

    压缩因子Z是相同的压力和温度下，实际气体与理想气体的比容比值。

    显然，理想气体的Z=1，实际气体的Z一般不等于1。Z值偏离1的大小，反映了实际气体对理想气体性质的偏离程度。

    1949年的R-K方程，根据  ，推出压缩因子Z的对比态方程形式[2]：

    其中：  ，    

   同样是上例，查表可得CO的Pc=35×105Pa，Tc =133K，用R-K方程计算得到P=6.97 ?105Pa，与实际测量值的偏差为：d=-1.69% ，可见，引入压缩因子后，可使气体测量的度大大提高。

    三、 压缩因子Z的确定
   
    如何地确定压缩因子Z，以获得满意的测量精度，自然就成了流量计量的主要问题。前人在此已做了大量的工作，根据实验数据整理成多种二参数的（Pr，Tr ）压缩因子图，其中，以Nelson-Obert提供的通用压缩因子图（简称N-O）图zui为准确。

    通过压缩因子图，根据Pr和Tr的值，就可以确定压缩因子Z。

    也可以根据压缩因子Z的对比态形式：      

    只要知道Pc和Tc，就可以通过迭代法用上式计算任意温度与压力下各种气体的压缩因子Z及相应的容积V 。开始迭代时，可先设Z=1， 确定出h值，代入Z方程中可得到一个新的Z值，然后再求出相应的新的h值和Z值，如此迭代到h值和Z值的变化小于允许误差为止。

    四、 压缩因子Z的单片机的实现

    单片机技术特别适合用在各种、小型、省电、可移动的场合。以单片机为内核，对流量进行测量的智能化仪表也随着单片机技术的发展而得到了广泛的应用。对气态物质进行在线测量时，由于气体的可压缩性，对测量体积的动态修正就成了必然的要求。引入压缩因子Z后，由于单片机字长的有限性、处理速度的限制、运算指令的限制，如何实时、准确的确定压缩因子Z，成为单片机算法处理的至关重要的问题。

    为了提高气体流量计量的精度，我们采用了美国TI公司的产品MSP430系列单片机[3]实现了气体流量的在线温度压力补偿、压缩因子Z修正的高精度气体流量计量。

    该系统由温度传感器、压力传感器、流量传感器、电池电源管理电路、系统时钟、数据远程传输接口、大容量可修改存储器、键盘、超低温LCD显示器等部分组成。见图2所示。

    温度、压力数据的采集采用14位的A/D；流量数据的采集采用硬件脉冲计数和软件测量相结合的方法，大大的提高了数据采集的精度，下面重点介绍有关流量采集的数据处理方法。

   1、被测气体的体积流量总量与信号脉冲数间的换算公式：V1=  m3

    式中： V1—气体工况下的体积流量（m3）；

    K—流量计仪表系数（1/m3）；

    f—流量信号脉冲数。

    为了在比较宽的量程比范围内实现等精度流量计量，考虑到流量传感器的非线性，我们在f—V1数据关系之间采用了非线性数据插值表处理，大大提高了原始测量数据的精度。

    2、被测气体的体积流量总量换算为标准状态下体积流量总量的公式：      （m3）

    式中： P0—标准大气压（P0=101325Pa）；      

    P1—气体工况压力（Pa）；

    T0—标准状态下的温度（293K）；     

    T1—气体工况温度（K）；

    V0—标准状态下的体积流量（m3）；   

    V1—气体工况下的体积流量（m3）；

    Z—气体压缩系数。

    3、压缩系数Z的计算方法：

    Pr = P/Pm

    式中：P—实际测量压力（压力）；

    Pm—临界压力（通过键盘输入）。  

    Tr = T/Tm

    式中： T—实际测量的温度，单位为K；

    Tm—气体介质的临界温度（通过键盘输入），单位为K。

    利用 Z=g（Tr，Pr，T，P） 函数关系可以得到压缩系数。

    （1）Z值的*种确定法（插值法）：

     关于压缩因子Z的在线数据修正，在精度要求比较低的场合，可以采用通用的N-O图进行处理。具体作法是将不同Tr时的Z1=g1（Pr）关系制成一张数据表，根据测量的P值求得Pr，利用Z1=g1（Pr）关系求得Z1。再对不同Tr时的Z进行插值得到Z，即可求得气体的流量V0。

    （2）Z值的第二种确定法（迭代法）：

     由R-K方程推出的压缩因子Z的对比态方程形式，通过迭代法，可以获得满意的迭代精度。但是，这样处理也存在两个方面的问题。一是由于单片机的处理速度的限制，迭代运算耗费了单片机的大量资源，而使得运行处理速度变得非常的缓慢，甚至有崩溃的可能。二是R-K方程本身仍然是对真实气体的一个近似描述的状态方程，测量的精度*依赖于R-K方程的精度。

    （3）Z值的第三种确定法（数学模型解算法）：

    由于单片机的处理速度限制、R-K方程以及所测量煤气、天燃气的适应性、精度的限制。我们根据天燃气、煤气的特点和*化的数学方法，在所研究的的对比温度Tr范围内，我们建立了一个Z=g（Tr，Pr，T，P） 数学模型[4]，[5]，[6]。通过在线实时压缩因子Z修正，可以获得高精度的流量计量。实践证明，简单易行，获得了满意的测量精度。

    4、工况下的瞬时体积流量与累积体积流量间的换算公式：            

     （m3/h）                 

    式中： Q1—气体工况下的瞬时体积流量（m3/h）；

    t—时间（h）。      

    5、标准状态下的瞬时体积流量与累积体积流量间的换算公式：   （m3/h）                

    式中： Q0—标准状态下的瞬时体积流量（m3/h）。

    五、 结论

    实践证明该单片机系统集温度、压力、流量测量，温度与压力补偿、压缩因子Z在线数学模型修正、信息存储与管理、网络数据传输、远程检测、收费和管理于一体的智能数字化流量计量系统。创造性地采用磁电隔离技术测量流量，有效降低始动压力（20Pa）；采用非线性脉冲—流量插值，提高量程比（接近1:15）；并可实时进行温度和压力补偿，通过优化的数学模型算法解算进行压缩因子Z修正，从而实现了气体流量的高精度计量（优于±0.5%）。可广泛应用于煤气、天然气、沼气等非严重腐蚀气体的流量计量和控制。