用程序计算器计算天然气流量

　　1　天然气流量计算的数学模型

　　根据石油部SYL04-83标准,天然气流量计算公式为:

　　H格———差压百分开方格数;

　　P格———上游气流绝对静压百分开方格数;

　　Z———工作状态下的压缩系数;

　　T1———气流绝对温度,K;

　　t1———气流温度,℃;

　　β———孔径比;

　　d———计量孔板内径,mm;

　　D———计量管道外径,mm;

　　Hmax———流量计差压量程,MPa;

　　TC———天然气假临界温度,K;

　　PC———天然气假临界压力,MPa;

　　αo———特定流量系数;

　　Fr———常用流量下的雷诺系数;

　　rRE———计量管粗糙度系数;

　　bK———孔板锐利度系数;

　　G———天然气相对密度;

　　Ko———天然气流量计算K值。

　　2　程序编制

　　天然气流量计算程序采用BASIC语言编

　　程,在CASIO PB-120程序计算器上进行。程

　　序清单如下:

　　5　DEFM40

　　10　INPUT“n=”,X,“H=”,H,“P=”,P,“t=”,T:T=273.15+T

　　15　RESTORE(100+X)

　　20　READ A,F,R,B,C,G,I,J

　　25　RESTORE(100+10*X)

　　27　READ U,V,D

　　30　K=4.703 1/10 000*A*F*R*B*C↑2*(I*J/G)↑0.5

　　40　M=4.6/10 000*I*V↑3.72/U

　　45　Z=1-M*P↑2*T↑-3.72

　　50　L=(0.41+0.35*(C/D)↑4)/100*J/13 323/I

　　55　E=I-L*(H/P)↑2

　　60　Q=K*E*(1/Z/T)↑0.5*P*H

　　65　PRINT“Q”;STR $ (X);“=”;INT(Q+0.5):GOTO 10

　　101　DATA 0.626 8,1.0,1.001 4,1.0018,132.511,0.577 4,6.0,40 000

　　102　DATAαo,Fr,rRE,bk,d,G,Pmax,Hmax

　　………………………………………………

　　110　DATA 4.61,191.6,259.00

　　120　DATAPC,TC,D

　　………………………………………………

　　(1)程序中变量X代表各输气管线编号,“1”表示1#管线,“2”表示2#管线,依此类推。

　　(2)数组语句中的αo、Fr、rRE、bk、d、G、Pmax、Hmax、PC、TC、D分别用变量A、F、、R、　B、C、G、I、J、U、V、D表示,其数值由分管天然气计量的工程技术人员确定。