Matlab编程天然气压缩因子计算模型.docx

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 模型适用条件 绝对压力：0MPa＜P＜12MPa 热力学温度：263K≤T≤338K 高位发热量：30MJ·m-3≤HS≤45 MJ·m-3 相对密度：0.55≤d≤0.80 天然气中各组分的摩尔分数应在以下范围内： CH4：0.7≤xCH4≤1.0 N2：0≤xN2≤0.20 CO2：0≤xCO2≤0.20 C2H6：0≤xC2H6≤0.10 C3H8：0≤x C3H8≤0.035 C4H10：0≤x C4H10≤0.015 C5H12：0≤x C5H12≤0.005 C6H14：0≤x C6H14≤0.001 C7H16：0≤x C7H16≤0.0005 C8H18和更高碳数烃类： C8H18：0≤x C8H18≤0.0005 H2：0≤x H2≤0.10 CO：0≤x CO≤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*，Uij，Kij，Gij* 。 2.3.2 待求量 压缩因子 Z 2.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方程方程，得到新的Z值，当（Z-Z0）的绝对值小于0.000001时，停止迭代，得到Z值。 3 程序代码 function [ Z ] = YSYZ( T,p,x) %计算天然气给定组分的压缩因子 % x为天然气组分，按照CO2 N2 H2 CO CH4 C2H6 C3H8 i-C4H10 n-C4H10 i-C5H12 n-C5H12 % n-C6H14 n-C7H16 n-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,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,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,0,0,0,0,0,0,0,0,0]; q=[0,0,0,0,0,0,1,0,0,0,