气体限流孔板的工艺计算
- 格式:pdf
- 大小:1.03 MB
- 文档页数:6
CHEMICALENGINEERINGDESIGN化工设计2021,31(2)
气体限流孔板的工艺计算
孔令伟 中国寰球工程有限公司北京分公司 北京 100012
摘要 限流孔板作为石油化工行业常用的限流降压设施,在工程设计中得到广泛应用。对气体而言,限流孔板孔径的计算方法有很多种。本文介绍了设计当中较为常用的三种亚临界气体流动限流孔板孔径计算方法以及五种临界气体流动限流孔板孔径计算方法。通过对比分析,指出每种计算方法的特点和易用性,并给出使用建议。通过两个实例计算发现,对临界流动气体而言,几种方法的计算结果偏差较大,针对这种情况,本文提出了建议的解决方法。关键词 孔板 限流 降压体 气体 亚临界 临界
孔令伟:高级工程师。2006年毕业于中国海洋大学应用化学专业获硕士学位。主要从事油气集输与处理方面的设计与研究工作。联系电话:18618127116,E-mail:konglingwei-hqc@cnpccomcn。 限流孔板作为节流元件,具有结构简单、易加
工、制造成本低、安装方便等优点,在满足工艺要
求的前提下,使用限流孔板代替调节阀来限定流量
或降低压力,可大幅降低投资和操作维修费用。根
据工艺要求,限流孔板可以用来:①限定流量;
②降低压力;③同时限流降压。
气体流过孔板时,其流动状态分为两种:亚临
界流动和临界流动。亚临界流动主要用于管道降压
目的,临界流动主要用于限制流量目的。实际上,
不管是单纯的限定流量或是降低压力,单级孔板即
可满足要求。如果对限定流量和降低下游压力均有
要求,则需要采用多级限流孔板来实现。通常情况
下,对限制流量和降低压力均有要求情况的非常少
见。本文分别针对亚临界流动和临界流动情况下的
限流孔板的计算进行讨论。
有关气体限流孔板孔径的确定,相关的计算公
式很多。比如,《工艺系统工程设计技术规定》[1]
HG/T20570-95、《炼油装置工艺管道安装设计手
册(下册)》[2]、《佩里工程师手册》[3]以及ISO
5167[4]等,各工程公司和设计院所依据的计算方
法也不完全相同。笔者在工作过程中发现,一些设
计人员在设计过程中对限流孔板的原理及公式理解
不透彻,会导致错误使用某些公式。本文选取了几
种较为常用的计算方法进行对比分析,并给出使用
建议。需要指出的是,本文未选取ISO5167中的方法进行讨论,因为ISO5167主要为孔板流量计
规范,涉及过多参数且有严格的使用条件限制,比
如要求孔板下游与上游压力比要大于等于075,
这就大大限制了其公式的适用范围。
1 孔径计算方法
11 亚临界流动
当流体流过一定孔径的孔板时,孔板前后就会
存在一定压差,它与流体流量(流速)的大小有
关,流量越大,压差越大;当一定压力、流量的流
体流过某一孔板时,通过改变孔径的大小则可以控
制孔板下游压力的大小。正是基于此原理,限流孔
板可起到管道降压作用。
(1)方法一
国内最常用的方法是依据《工艺系统工程设
计技术规定》HG/T20570-95[1][5-7]给出的单板
计算公式,其计算孔板孔径公式如下:
W=4378×C×d20×
P1M()ZTkk()-1P4P()
12k-P4P()
1k+1[]槡k(1)
式中,W为流体的质量流量,kg/h;
C为孔板流量
系数;d0为孔板孔径,m;D为管道内径,m;P1为孔板前压力,Pa;P4为孔板后压力,Pa;M为
相对分子质量;Z为压缩系数;T为孔板前流体温
度,K;k为绝热指数。012021,31(2)孔令伟 气体限流孔板的工艺计算
由公式(1)可知,欲求出限流孔板孔径d0,
需先确定孔板流量系数C,而孔板流量系数C又由
雷诺数Re和孔径比β(d0/D)共同决定,故需要
先假定C值,经迭代计算求出限流孔板孔径d0。
然后根据孔径比β(d0/D)值和Re值由“C-Re
-d0/D关系图”查得孔板流量系数C′,根据C′与
C的偏差大小,判断计算出的d0是否有效,否则
按查出的C′值作为输入值,重新计算,直到满足
要求为止。
由计算过程可知,应用式(1)计算限流孔板
的孔径既有C和d0的迭代计算,同时又有查图过
程,较为不方便,不利于大量计算。
(2)方法二
文献[8]中介绍了一种亚临界气体流动限流孔
板的计算方法:
q=Y×C×A2×ΔPρ槡1(2)
式中,q为流体的实际体积流量(按孔板上游状态
计),m3/s;Y为气体膨胀因子,可根据压比△P/
P1及孔径比β(d0/D)从以下图1、图2中查取;
C为孔板流量系数,意义同公式(1);A为孔径截
面积,m2;△P为孔板前后压差,P1-P4,MPa;
ρ1为孔板前气体密度,kg/m3。
由公式(2)可知,欲求限流孔板孔径d0,需
先确定孔板流量系数C及膨胀因子Y。而孔板流量
系数C与Re和孔径比β(d0/D)值有关,膨胀因
子Y与β(d0/D)值有关,故需要先假定d0值,
查图得到孔板流量系数C及膨胀因子Y,经迭代计
算求出限流孔板孔径d0′。然后根据孔径比β(d0/
D)值和Re值由“C-Re-d0/D关系图”孔板流
量系数C′,同时,根据图1或图2查出膨胀因子
Y′值。根据C′与C的偏差大小、Y′与Y的偏差大
小以及d0′与d0的偏差大小,判断计算出的d0′是
否有效,否则计算出的d0′值作为输入值,重新计
算,直到满足要求为止。
由计算过程可知,应用公式(2)计算限流孔
板的孔径既有C、Y以及d0的迭代计算,同时又
包含查图过程,极为不方便。
(3)方法三
对亚临界气体流动,有工程公司采用一种来自
于佩里工程师手册第10
章[3]
中介绍的方法:图1 膨胀因子(k≈13)
图2 膨胀因子(k≈14)
d0=4WπKYCc/2gc(P1-P4)槡槡ρ(3)
K=05982+01941β-04981β2+13679β3
(4)
Y=1-1-γk(041+035β4)(5)γ=P2P1(6)
P2=P1-P1-P41-β2(7)
式中,d0为限流孔板孔径,m;W为气体质量流
量,kg/s;gc为重力加速度,=98,m/s2;P1为
孔板前压力,kg/cm2;P2为孔板中心处压力,kg/
cm2;P4为孔板后压力,kg/cm2;ρ1为孔板前气体
密度,kg/m3;
K为排出系数;Y为气体膨胀系数;
Cc为转换常数,=100;β为孔径比,d0/D;k为
绝热指数。11 CHEMICALENGINEERINGDESIGN化工设计2021,31(2)
由公式(3)~(7)可知,欲求限流孔板孔径
d0,需先确定孔板流量系数C和膨胀系数K,而孔
板流量系数C及膨胀系数K均与压比β(d0/D)
有关,故需要先假定β值,经迭代计算求出限流
孔板孔径d0,进而得到β′值,然后根据β与β′的
偏差大小判断计算出孔板的孔径d0是否可以接受,
否则按计算出的β′值作为新输入值,重新计算,
直到满足要求为止。
由计算过程可知,应用式(3)~(7)计算限
流孔板的孔径仅包含迭代计算,不含查图过程,利
用Excel表即可实现自动计算。
12 临界流动
当孔板后的压力降至低于一定数值(临界压
力Pc),流体通过孔板缩孔处的流速达到音速,不
论如何降低出口压力,只要孔板上游的压力保持一
定,流量将维持一定的数值而不再增加。基于气体
临界流动的性质,限流孔板可用作限流作用。
实际上,根据国外相关设计资料,以及作者参
与的国外工程设计,气体管线上的孔板几乎全部采
用的是单板设计。并且在实际的操作中,单板的效
果很好,检维修也十分快捷方便,这点在文献[6]
中也有提及。
(1)方法一
此方法从亚临界流体流动公式(1)演化得
来,用气体临界压力Pc替换孔板后压力P4:
W=4378Cd20×
P1M()ZTkk()-1PcP()
12k-PcP()
1k+1[]槡k(8)
式中,Pc为临界限流压力,Pa;其他参数意义同
公式(1)。
其计算过程同公式(1)完全相同,计算限流
孔板的孔径既有C值的迭代计算,同时又包含查
图过程,较为不方便。
(2)方法二
文献[8]中介绍了一种临界流体的计算方法,
其公式(9)从公式(2)演化而来,用气体临界
压力Pc替换孔板后压力P4:
q=Y×C×A2P1(1-yc)槡ρ(9)
式中,γc为临界压比,Pc/P1;其他参数意义同公式(2)。
此方法计算过程同公式(2)完全相同,应用
公式(9)计算限流孔板的孔径既包含C、Y以及
d0的迭代计算,同时又有查图过程,很不方便。
(3)方法三
有工程公司采用了一种临界流体的计算方
法[9]:
d0=4WπCP1/gckM()ZRT2k()+1k+1k槡槡-1(10)
C=08484-0099P4P()
1+01167P4P()
12-0630P4P()
13
(11)式中,d0为限流孔板孔径,ft;W为气体质量流
量,lb/s;gc为重力加速度,=3217,ft/s2;P1为孔板前压力,lb/ft2;P4为孔板后压力,lb/ft2;
C为流量系数,无量纲;k为绝热指数;Z为压缩
因子;R为1546,ft·lb/(lb-mol)·F;T为孔板
前流体温度,°R;其他参数意义同公式(3)。
由式(10)、(11)可知,采用此方法计算气
体临界流动状态下限流孔板无迭代计算,也无需查
图,因此非常简单。需要注意的是,其计算过程需
要采用英制单位。
(4)方法四
第四种计算方法来自《炼油装置工艺管道安
装设计手册(下册)》[2]
d20=W1252×α×ερ1×Δ槡P(12)
式中,d0为孔板孔径,mm;W为流体的质量流量,kg/h;α为流量系数;ε为膨胀系数;ρ1为孔
板前气体密度,kg/m3;△P为孔板前后的压力降,kgf/cm2。
公式(12)中涉及到的参数,比如流量系数α和膨胀系数ε需要查图或查表,因此便捷性差一
些。另外,实际上式中的△P仅与孔板前压力P1有
关,与孔板后压力无关,因此公式(12)的表达不
清晰。
(5)方法五
文献[
5][10]中介绍了一种临界气体流动条件下
限流孔板孔径的计算公式:
d0=2812Wφ×1P1×ρ槡槡1(13)21