资源描述：
《matlab编程天然气压缩因子计算模型》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、1程序目的利用AGA8-92DC模型计算天然气的压缩因子，该程序主要应用于在输气和配气正常进行的压力P和温度T范围内的管输气的压缩因子计算2数学模型：AGA8-92DC模型2.1模型介绍此模型是已知气体详细的摩尔分数组成和相关压力、温度来计算气体压缩因子。输入变量包括绝对压力、热力学温度和摩尔组成。摩尔组成是以摩尔分数表示下列组分：CO2、N2、H2、CO、CH4、C2H6、C3H8、i-C4H10、n-C4H10、i-C5H12、n-C5H12、n-C6H14、n-C7H16、n-C8H18。2.2模型适用条件绝对压力
2、：0MPa＜P＜12MPa热力学温度：263K≤T≤338K高位发热量：30MJ·m-3≤HS≤45MJ·m-3相对密度：0.55≤d≤0.80天然气中各组分的摩尔分数应在以下范围内：CH4：0.7≤xCH4≤1.0N2：0≤xN2≤0.20CO2：0≤xCO2≤0.20C2H6：0≤xC2H6≤0.10C3H8：0≤xC3H8≤0.035C4H10：0≤xC4H10≤0.015C5H12：0≤xC5H12≤0.005C6H14：0≤xC6H14≤0.001C7H16：0≤xC7H16≤0.0005C8H18和更高碳数烃
3、类：C8H18：0≤xC8H18≤0.0005H2：0≤xH2≤0.10CO：0≤xCO≤0.03如果已知体积分数组成，则应将其换算成摩尔分数组成。所有摩尔分数大于0.00005的组分都不可忽略。2.3模型描述2.3.1已知条件绝对压力P、热力学温度T、组分数N；各组分的摩尔分数Xi，i=1～N；查附表1、2、3得到的以下数据：58种物质的状态方程参数an，bn，cn，kn，un，gn，qn，fn，sn，wn；14种识别组分的特征参数Mi，Ei，Ki，Gi，Qi，Fi，Si，Wi；14种识别组分的二元交互作用参数Eij*
4、，Uij，Kij，Gij*。2.3.2待求量压缩因子Z2.3.3计算步骤a)第二维利系数B的计算：二元参数Eij和Gij，由以下两式计算：b)计算系数Cn*，n=13～58用以下方程求解混合方程，计算混合物参数U，G，Q。c)计算混合物体积参数K；d)计算对比密度摩尔密度为：式中，P为绝对压力，Mpa；R为摩尔气体常数；T为热力学温度，K。对比密度ρr同摩尔密度ρm相关：e)利用AGA8-92DC方程，对压缩因子进行迭代计算迭代过程：给出Z0的初始值为1，先计算出ρm，将ρm、K和已知量带入AGA8-92DC方程方程，得
5、到新的Z值，当（Z-Z0）的绝对值小于0.000001时，停止迭代，得到Z值。3程序代码function[Z]=YSYZ(T,p,x)%计算天然气给定组分的压缩因子%x为天然气组分，按照CO2N2H2COCH4C2H6C3H8i-C4H10n-C4H10i-C5H12n-C5H12%n-C6H14n-C7H16n-C8H18的顺序输入%T为温度，单位为K%P为压力，单位为兆帕N=14;R=8.314;%状态参数值b=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,
6、2,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,5,5,5,5,5,6,6,7,7,8,8,8,9,9];k=[0,0,0,0,0,0,0,0,0,0,0,0,3,2,2,2,4,4,0,0,2,2,2,4,4,4,4,0,1,1,2,2,3,3,4,4,4,0,0,2,2,2,4,4,0,2,2,4,4,0,2,0,2,1,2,2,2,2];c=[0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,1,1,0
7、,0,1,1,1,1,1,0,1,1,1,1,0,1,0,1,1,1,1,1,1];g=[0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,1,0,1,0,0];f=[0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
8、0,0,0,0,0,0,0,0,0];q=[0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,1,0,0,0,0,0,1];s=[0,0,0,0,0,0,0,1,1
