摘要:目的将氢气掺入天然气管道中会改变管道内气体的性质和流动状态，可能会影响标准孔板流量计计量精度,采用ANSYSYFluent对混氢天然气管道标准孔板流量计进行适应性研究。方法比较了不同混氢量的天然气对流出系数、可膨胀系数、相对密度系数、超压缩系数、流速及差压的影响。结果在303.15K.3MPa,混氢量为0%~30%的条件下,随着混氢量的增加,会导致差压上升;导致相对密度系数、可膨胀系数和超压缩系数下降;导致流速上升,使测量流量增加。结论由于氢气的发热量低于天然气,因此，针对混氢天然气，建议采用能量计量�；烨馓烊黄换岫员曜伎装辶髁考凭炔洗笥跋臁�
氢能是一-种绿色、低污染、可再生的燃料，被认为是最有前途的化石燃料替代品之一口。目前，利用可再生能源电解制氢,然后将氢气按照一定比例掺人天然气管道中进行输送是利用和运输氢能的有效途径[5。如IEAGHGR&.D项目掺人天然气管网中的氢气摩尔分数高达25%叫;AMeland项目掺人天然气管网中的氢气摩尔分数达到20%[1-8]。而掺氢天然气计量技术是掺氢天然气产规�；褪谐』闹匾�。标准孔板流量计由于其设计简单、成本低,仍然是石油与天然气行业中使用广泛的流量计。
由于氢气和甲烷物性差异巨大，在标况下其密度相差8倍以上[1],而密度是影响标准孔板流量计结果的重要因素[18]。当天然气中掺混氢气后，会导致其密度、黏度、比热容参数改变,进而影响标准孔板流量计计量精度。Dong等利用Fluent分析不同倾角孔板在测量天然气流量时对测量精度的影响;Jin等利用Fluent分析得到在测量液氢时不同孔板结构对流出系数和压力损失系数的影响;通过数值模拟技术得到在测量天然气流量时,流体相对密度变化值对测量值有较大的影响。
尽管前人已经做了很多研究,但目前对影响测量精度的研究主要集中在孔板结构的变化上,这将会增加流量计结构的复杂性，而且在实际的天然气管道中不易使用。此外,所研究的结论主要是对流出系数、差压等的影响,关于其他计量所需参数,如可膨胀系数、超压缩系数、相对密度系数的影响很少被研究。研究的介质主要是天然气或者液氢,关于混氢天然气的情况很少被研究。因此,本研究主要分析天然气管道中混人氢气后对标准孔板流量计测量精度的影响。
1标准孔板流量计工作原理
标准孔板流量计以能量守恒定律和流动连续性方程为基础,通过测量孔板前后产生的静压力差来衡量天然气流过节流装置的流量大小”。工况条件下的体积流量一般用流量计测量,然后换算成基本(标准)条件下的体积作为天然气贸易交接过程中的流量8]。GB/T21446-2008《用标准孔板流量计测量天然气流量》以293.15K.101.325kPa为条件,得到标准条件下天然气体积流量计算实用公式，如式(1)所示:
式中:qV。为标准条件下天然气体积流量,m³/s;Avn为体积流量系数，Avn。=3.1795X10-6;C为流出系数;E为渐进速度系数,E=1/(1-β)0.5;β为孔径比,β=d/D;d为孔板开孔直径,mm;D为测量管内径,mm;Fc为相对密度系数;ε为可膨胀系数;Fz为超压缩系数;Fr为流动温度系数;p1为孔板上游取压孔实测绝对压力,MPa;△p为孔板前后差压,MPa。
流出系数C的计算公式如式(2)~式(4)所示。
式中:ReD管径为雷诺数;L1为孔板上游端面到取压孔轴线的距离除以测量管内径得出的商;L2为孔板下游端面到取压孔轴线的距离除以测量管内径得出的商;M2为变量;A为变量。
2数值仿真模型建立及验证
2.1孔板结构
孔板结构示意图如图1所示。针对3种孔径比进行研究,孔板几何形状:孔板厚度为3.8mm,孔板开孔厚度为0.8mm,上游管径为150mm,孔板孔径分别为57mm、75mm、87mm,孔径比分别为0.38、0.50、0.58。本研究选择孔板上游直管段145D,下游直管段10D,以获得准确的模拟结果。
2.2计算网格划分
采用ANSYS建立了标准孔板流量计的三维模型,利用六面体网格对网格进行划分。在模拟中,整个几何形状被分为3个区域:上游、中心区域、下游。上游和下游区域使用较粗网格,中心区域采用更密的.网格,以获得压力梯度。墙附近的网格被细化,以满足标准墙功能的要求。管道模拟网格如图2所示。进行了网格尺寸独立性测试，用来数值模拟结果与网格尺寸和网格质量无关。以3MPa下氢气摩尔分数分别为0.0、0.4的CH-H2混合物为例,采用1267153、1893462、2637960、3439231个单元进行测试。网格数量从1893462增加到3439231时，网格数量对孔板前后的压力的影响已经很小了�？悸峭竦奈薰匦院图扑阈�，在以下模拟中采用2637960个单元的网格。
2.3控制方程
假设:实际流体在管道中做定常流运动;气质组分为甲烷和氢气混合物,且混合均匀;流体在管道内与外界无热量交换。因此,除了满足质量、动量和能量三大守恒方程外，还需满足气体状态方程。本研究使用SRK状态方程[21],如式(6)所示。
式中:p为压力，MPa;R为气体常数，8.314J/(mol·K);T为温度,K;V为摩尔体积,m³/mol;αe。为临界参数,是临界温度和临界压力的函数;α为引力函数,是对比温度和偏心因子的函数;b为斥力函数�；剐璺治黾淄楹颓馄诠艿乐衅宕使媛�，因此,开启组分输运模型,如式(7)所示:
式中:ρ为密度,kg/m³;ci为i组分的体积分数;t为时间,s;u为速度,m/s;Di为i组分的扩散系数,m²/s;Ri为单位时间、体积下产生i组分的质量，kg/(m³.s)。
针对天然气计量,还需结合湍流方程。K-εRNG模型在湍流模拟中得到了广泛的应用。与标准的kε模型相比，K-εRNG模型在表征具有强流线曲率、涡旋方面都有了显著的改进15]。因此,本研究选择kεRNG模型作为湍流方程。
2.4边界条件
选择3MPa压力边界进行计算。模拟的边界条件为:进口边界条件采用天然气压力,出口边界条件采用天然气出口流量。进口温度设置为303.15K,流体介质采用甲烷和氢气混合物,并由软件本身的数据库确定了其密度、黏度等参数。令x(CH4)和x(H2)分别为甲烷和氢气摩尔分数,边界条件设置见表1。
2.5有效性验证
基于流体相似原理,可利用Fluent计算在计量管内径为30mm,孔径比为0.42、0.59、0.65条件下水的流出系数，与实验值进行对比，对本研究模型有效性进行验证。验证结果如表2所列。
从表2可以看出,采用数值模拟方法计算出的流出系数与实验值吻合较好,偏差不超过-3.50%。
3结果与讨论
3.1混氢量对差压的影响
以孔板孔径比为0.38,x(H2)为0.00、0.10、0.20、0.30为例,Fluent仿真结果压力云图见图3�？拙侗任�0.38、0.50、0.58的标准孔板的差压随混氢量的变化如图4所示。
从图4可以看出,随着混氢量的增加,流过标准孔板的差压会逐步上升。从数值上看,孔径比越小,差压随混氢量的增加而上升的幅度越明显,这说明氢气对孔板的节流效应比较敏感
3.2混氢量对流速的影响
以孔板孔径比为0.38,x(H2)为0.00、0.10、0.20和0.30为例，Fluent仿真结果速度云图见图5。从图5可以看出,随着混氢量的增加,气流流过孔板后的速度更大。图6所示为混氢量与输送速度的关系图,从图中可看出,混氢量越高,流速越高。
因此,当天然气管道中掺入氢气后会导致流量增大。由于氢气的发热量小于甲烷,若仍然采用体积计量进行贸易交接,这将会对买方不利。若采用质量计量进行贸易交接,仍然不能合理体现掺氢天然气的实用价值,对供方不利。因此,针对混氢天然气,建议采用能量计量进行贸易交接。
3.3混氢量对流出系数的影响
采用式(2)计算得到不同混氢量下的流出系数,计算结果见图7。从图7可以看出:孔径比越大,流出系数越大;在混氢量小于0.3时,混氢量的变化几乎不会对流出系数产生影响。
3.4混氢量对相对密度系数的影响
相对密度系数变化与孔板结构无关,仅与组分的变化有关,图8所示为相对密度系数随混氢量的变化情况。从图8可看出,混氢量的增加会导致相对密度系数上升,这是由于氢气的摩尔质量远小于甲烷,混氢量的增加会导致其摩尔质量下降,进而导致相对密度系数上升。
3.5混氢量对可膨胀系数的影响
图9所示为可膨胀系数随混氢量的变化。从图9可以看出,随着混氢量的增加，会导致可膨胀系数下降,在低孔径比的情况下，其下降幅度要大于高孔径比,但整体下降幅度较小
3.6混氢量对超压缩系数的影响
超压缩系数是因天然气特性偏离理想气体定律而采用的修正系数，其与孔板结构无关。分析在303.15K,3MPa.5MPa和7MPa条件下的超压缩系数随混氢量的变化(见图10)。从图10可以看出，超压缩系数随混氢量的增加而下降,压力越大,下降幅度越大。
3.7混氢量对标准孔板流量计测量精度的影响
基于Fluent模拟结果,得到孔板前后压力、温度、黏度等参数,采用式(1)~式(5)计算得到的流量作为标准孔板流量计测量流量,以边界流量作为实际流量进行对比分析,分析结果见图11�；诒狙芯拷⒌募扑隳Ｐ偷玫奖曜伎装辶髁考频牟饬苛髁坑牍艿澜孛娴氖导柿髁恐涞牟饬课蟛�,其计算公式如式(8)所示。
式中:δ为测量误差,%;qbou为实际流量,m³/s;qea为测量流量(基于本研究建立的计算模型通过Fluent模拟计算得到的标准孔板流量计流量),m2/s。
从图11(a)可以看出随着混氢量的增加,标准孔板流量计测量流量也会显著增加。从图11(b)可以看出,标准孔板流量计计量精度几乎不受混氢量变化的影响。
4结论
采用数值模拟的方法,研究了标准孔板流量计应用于混氢天然气时的计量精度。研究了混氢量对差压、流速、流出系数、相对密度系数、可膨胀系数和超压缩系数的影响,可得到以下结论。
(1)在压力一定的情况下,混氢量的增加会导致体积流量测量的流量值增大。因此,针对混氢天然气,建议采用能量计量进行贸易交接。
(2)在压力一定的情况下,混氢量的增加会导致差压上升,导致相对密度系数、可膨胀系数和超压缩系数下降,而流出系数几乎不受氢气含量变化的影响。
(3)将氢气掺人天然气管网,在氢气摩尔分数小于30%的情况下,氢含量的变化不会对标准孔板流量计精度产生明显的影响。
本文来源于网络,如有侵权联系即删除！
