富氢气流量测量与纯氢流量计算-江苏省苏科仪表有限公司

摘要：供需双方贸易交接计量以富氢气体中的纯氢量为准，由于测量的富氢气体中含有饱和水蒸气，随着温度的变化，其中的水蒸气与氢的体积分数均会发生变化，导致氢的平均相对分子质量不停变化，为测量带来困难。针对该问题，对比分析了涡街流量计方法、科氏力质量流量计方法、孔板流量计方法三种测量方案，最后选取孔板流量计结合热导式氢分析仪测量氢的体积分数，通过在流量演算器中进行密度补偿、干气部分计算、纯氢流量计算，实现了准确测量、数据可信。
　　某化工厂从附近的氯碱厂输入纯度为98％的氢气，用于生产氢体积分数φＨ2为70％～90％的富氢气体，再经1根内径为468ｍｍ的管道返送回氯碱厂，混合气中的水蒸气达到饱和状态，通过热导式氢分析仪测量混合气中的氢体积分数，测量范围为50％～100％。该分析仪测量出的是干气体中的氢体积分数，而工艺要求测量的是湿气体干部分流量，并与氢体积分数一起计算纯氢的体积流量，因此文中所述是一个变组分气体的湿气体干部分流量测量问题。
1测量方案的确定
1．1涡街流量计方法
1）该方法信号处理较为简单。涡街流量计是一种体积流量计，即输出信号仅与流过旋涡发生体处的流体平均流速成正比，而与物料的密度、黏度等无关。如果能用涡街流量计测量某种气体的体积流量，则信号经下面的处理后，就可得到纯氢体积流量。
ａ）将湿气体中的水汽量扣除。
B）将湿气体的干部分流量换算为标准状态流量。
ｃ）将测量到的干气体中氢体积分数值与湿气体的干部分体积流量相乘，即得到纯氢的体积流量。
由于测量对象具有管径大、流速低、静压低、流体密度小的特点，用涡街流量计测量是否可行，须经分析和具体验算才可确定。
2）已知条件。该测量对象与最低可测流速有关的参数包括：管道内径D20＝468ｍｍ，最高压力为70Kpa，常用压力为40Kpa，最低压力为20Kpa，干部分最大体积流量为1×104ｍ3/ｈ，干部分常用体积流量为7×103ｍ3/ｈ，最小体积流量为2．8×103ｍ3/ｈ，介质温度为21．0～63．5℃（常用50℃），常用工况流体密度ρ1＝0．4138KG/ｍ3。由于工况条件下的流体密度小于0．5KG/ｍ3，使用涡街流量计测量将受到流体密度的限制，因此根据上述数据计算得到常用工况条件下最小体积流量对应的流速为3．57ｍ/ｓ。
3）最低可测流速验算［1－2］。提供了大口径涡街流量计最低可测流速公式：

式中：ｖ———最低可测流速，ｍ/ｓ；ρ———流体密度，KG/ｍ3。
将ρ＝0．4138KG/ｍ3代入式（1）得到最低可测流速为8．73ｍ/ｓ，显然比最小体积流量对应的最低流速大得多。
4）缩小管径的可行性。缩径前后的管径与缩径流速的关系可用式（2）计算：

　　式中：D2———缩小后的管径，ｍｍ；D1———缩小前的管径，ｍｍ；ｖ2———缩小后的流速，ｍ/ｓ；ｖ1———缩小前的流速，ｍ/ｓ。
　　将D1＝468ｍｍ，ｖ1＝3．57ｍ/ｓ，ｖ2＝8．73ｍ/ｓ代入式（2），得D2＝299．3ｍｍ，即要选DＮ300涡街流量计才能将最小可测流量覆盖，如果混合气中的氢体积分数高于常用值或流体压力低于常用值，流体密度减小时将导致最小流量无法测量；而且管道流通截面积缩径后只有缩径前的41％，缩径导致流速增高，最大流量时的流速高达31ｍ/ｓ，远高于规程规定的安全流速［3］。因此，涡街流量计不适用富氢气体积流量的测量。
1．2科氏力质量流量计方法
　　用科氏力质量流量计直接测量变组分气体的质量流量，是个既省力又准确的方案，但要求被测气体静压在100Kpa以上，以便测量管内气体达到必要的流速。但在该案例中，静压最低时仅20Kpa，而且流速较高，也不符合安全流速的规定。
1．3孔板流量计方法
　　只要雷诺数大于5×103，采用该方法就能保证测量的精确度，不受下限流速的约束。但对于组分变化的气体，在温度、压力不变的情况下，其密度会随组分的变化而变化，从而产生误差，为此必须针对组分变化对测量的影响进行补偿。
2孔板流量计方法的实施

孔板流量计基本公式：

孔板流量计的基本公式如式（3）所示［4－6］：

　　式中：ｑm———质量流量，kg/s；Ｃ———流出系数；β———直径比，β＝D/D；ε1———节流件正端取压口平面上的可膨胀性系数；D———工作条件下节流件的开孔直径，ｍ；D———管道内径，ｍ；ΔP———差压，Pａ；ρ1———节流件正端取压口平面上的流体密度，KG/ｍ3。
当流量用体积流量表示时，式（3）转换为式（4）：

式中：ｑＶ———体积流量，ｍ3/ｓ。
2．2混合气密度计算公式
1）标准状态下湿气体干部分的密度ρgn用式（5）计算［7－8］：

　　式中：φｉ———干部分各组分的体积分数，％；ρGｎｉ———干部分各组分的密度，KG/ｍ3．
　　由于氢体积分数较大，该案例中只测量湿气体干部分中的氢体积分数，即认为除氢之外的几种其他气体的体积分数恒定不变，则式（5）简化成式（6）：

　　式中：φＨ2———氢体积分数，％；ρＨ2———标准状态下的氢气密度，KG/ｍ3；ρＮｎ———标准状态下氢之外其他气体的平均密度，KG/ｍ3。
2）工作状态下湿气体干部分密度ρG的计算如式（7）所示［7－8］：

　　式中：P1———工作状态下差压装置正端取压口绝对压力，Pａ；1———湿气体相对湿度，取1＝100％；Pｓ1ｍａｘ———工作状态下饱和水蒸气的分压，Pａ；Pｎ———标准状态下气体绝对压力，Pａ；Tｎ，Tｆ———标准状态和工作状态下气体的热力学温度，K；Zｎ，Zｆ———标准状态和工作状态下气体的压缩系数，取Zｎ＝Zｆ＝1。其中Pｓ1ｍａｘ由表1查得.。

3）工作状态下湿气体密度ρ1的计算如式（8）所示：

　　式中：ｑＶ1———工作状态下混合气干部分体积流量，ｍ3/ｈ；ｑＶｎ———标准状态下混合气干部分体积流量，通常由测量任务书给出，包括刻度最大流量ｑＶｎｍａｘ，常用流量ｑＶｎｃｏｍ和最小流量ｑＶｎｍｉｎ，ｍ3/ｈ。
2．4差压装置的计算
　　求得ｑＶ1和ρ1后，可按GB/T2624—2006设计计算差压装置。
3温度压力及ε1的补偿
3．1温度压力补偿氢体积分数的补偿公式
　　气体温度、压力偏离设计值时，氢体积分数也会偏离设计值，此时湿气体干部分流量可用式（11）进行补偿［10］：

式中：带“′”的物理量为实际使用工况条件下的参数，不带“′”的物理量为设计工况所对应的参数，在差压装置计算书中可以获得。
式（9）中，由于工况变化，ρ1已经从式（8）所表示的值变成式（12）所表示的值：

3．2流量示值的可膨胀性系数补偿
　　为了提高湿气体干部分流量测量的精确度，必须进行气体可膨胀性系数（ｅｘPａｎｓｉBｉｌｉTｙ）的补偿，并按GB/T2624—2006中的公式Ｃ＝ｆ（β，RｅD）进行雷诺数补偿，具体操作方法详见参考文献［10－11］。
3．3补偿公式的实现
　　式（15）中有一半的物理量是已知的，只需按仪表说明书的要求填入菜单，其余的则是测量值或计算值，具体见表2所列。

3．4结构上的处理
　　由于混合气中的水蒸气体积分数总是达到饱和程度，为防止冷凝水在孔板前、三阀组内、差压变送器高低压室内的聚积对测量的精确度产生影响，该项目仪表采用偏心孔板，并将差压变送器布置在水平管道上方，而且弃用针形阀三阀组，改用不锈钢球阀，从而杜绝了积水。
　　为便于清理差压信号导压管内可能生成的污物，差压装置采用径距取压方法。
4纯氢流量的计算
　　孔板流量计经密度补偿测得的湿气体干部分体积流量，换算为标准状态下的体积流量后，乘上干气体中的氢体积分数，就得到纯氢的体积流量。
5结束语
1）文中计算是变组分气体流量测量问题，组分的变化源自混合气中氢体积分数的变化和气体温度变化引起的水汽体积分数的变化。介质温度冬季最低达21℃，夏季最高达63．5℃，水汽体积分数相差9倍之多，因此混合气密度变化很大。
2）该案例中采用热导式氢分析仪测量气体干部分的氢体积分数，并根据气体中饱和水蒸气压力与饱和水蒸气密度表，扣除水汽体积分数，计算得到纯氢体积流量。
3）该案例中采用流量演算器实现各项计算，由于使用的是汇编语言，所以杜绝了病毒的侵扰，运行安全可靠。
4）该案例中的系统在东北某化工企业投运6ａ以来，运行可靠，未发生过故障。测量结果与产品产量一起计算得到的氢单耗不受季节影响，因此作为贸易交接手段得到了双方认可，从而解决了经营管理中的难题。

以上内容源于网络，如有侵权联系即删除!


江苏省苏科仪表有限公司是一家专业提供涡街流量计、涡轮流量计、电磁流量计的企业，公司将以优质的服务优惠的价格，服务新老客户。
版权所有：江苏省苏科仪表有限公司				技术支持│易品网络
温度仪表事业部压力仪表事业部流量仪表事业部校验仪表事业部显示仪表事业部变送器仪表事业部电线电缆事业部