天然气压缩因子设计计算书

天然气压缩因子设计计算书——AGA8-92DC与SGERG-88方法对比分析项目名称天然气长输管道压缩因子设计计算文件编号GC-ZF-2026-001设计阶段施工图设计阶段文件类型工程设计计算书编制单位油气工程设计研究院编制日期2026年10月编制审核批准

目录1概述 概述1.1计算目的与意义天然气压缩因子（又称偏差因子或Z因子）是描述天然气实际气体与理想气体行为偏离程度的重要热力学参数，在天然气输送管道设计、流量计量、能量计量以及储气库工况分析中具有关键作用。准确计算压缩因子对于保障管道输送安全、优化计量系统精度、提高能源利用效率具有重要的工程意义。本计算书针对某天然气长输管道工程设计工况，分别采用AGA8-92DC方法和SGERG-88方法进行天然气压缩因子计算，并对两种方法的计算结果进行系统对比分析，为工程设计提供可靠的技术依据。1.2计算范围本计算书涵盖以下内容：（1）天然气压缩因子计算原理与方法介绍；（2）基于摩尔组成的AGA8-92DC方法多工况计算；（3）基于物性参数的SGERG-88方法多工况计算；（4）两种计算方法的结果对比与误差分析；（5）计算结果安全核验与工程指导意义分析。1.3设计原则本计算书遵循以下设计原则：（1）严格执行国家现行标准和行业规范；（2）计算方法选取经过验证，精度满足工程要求；（3）多工况覆盖设计参数范围，留有适当安全裕量；（4）计算结果可追溯、可复核。2标准引用与编制依据本计算书的编制依据以下标准、规范和技术文件：表2-1主要引用标准与规范一览表序号标准编号标准名称1GB/T17747.1-2011天然气压缩因子的计算第1部分：导论和指南2GB/T17747.2-2011天然气压缩因子的计算第2部分：用摩尔组成进行计算3GB/T17747.3-2011天然气压缩因子的计算第3部分：用物性值进行计算4GB50251-2015输气管道工程设计规范5GB/T11062-2020天然气发热量、密度、相对密度和沃泊指数的计算方法上述标准均为现行有效版本，本计算书所有计算方法和参数取值均严格遵循上述标准的规定。3天然气压缩因子计算原理3.1压缩因子定义与物理意义压缩因子（CompressibilityFactor），又称偏差因子或Z因子，定义为实际气体的摩尔体积与同温度、同压力下理想气体摩尔体积之比，是衡量实际气体与理想气体行为偏离程度的无量纲参数。Z=p·Vₘ/(R·T)=V_real/V_ideal式中：Z——压缩因子，无量纲；p——气体绝对压力，Pa；Vₘ——气体摩尔体积，m³/mol；R——通用气体常数，取8.314J/(mol·K)；T——气体热力学温度，K。当Z=1时，气体行为符合理想气体状态方程；当Z<1时，分子间吸引力占主导，气体比理想气体更易压缩；当Z>1时，分子本身体积占主导，气体比理想气体更难压缩。在天然气长输管道的常规工况范围（压力1~10MPa，温度0~40℃）内，压缩因子一般在0.8~0.95之间。3.2AGA8-92DC方法原理AGA8-92DC方法（DetailedCharacterizationMethod）是美国燃气协会（AGA）于1992年提出的高精度天然气压缩因子计算方法，对应GB/T17747.2-2011和ISO12213-2:2006。该方法以扩展维利方程为基础，输入参数为天然气的详细摩尔组成。AGA8-92DC方程的基本形式为扩展维利方程：Z=1+B·ρ式中：Z——压缩因子；B——第二维利系数；ρₘ——气体摩尔密度，mol/L；ρᵣ——对比密度，ρᵣ=ρₘ/ρₘ,max；Cₙ*——高阶维利系数；bₙ、cₙ、kₙ——方程常数。AGA8-92DC方法包含58个维利系数项，通过迭代求解对比密度，最终得到压缩因子。该方法适用于管输天然气，适用压力范围为0~12MPa，温度范围为-10~60℃，对于典型管输天然气组成，计算不确定度优于0.1%。3.3SGERG-88方法原理SGERG-88方法（SimplifiedGERGVirialEquation）是欧洲燃气研究集团（GERG）于1991年提出的简化维利方程计算方法，对应GB/T17747.3-2011和ISO12213-3:2006。该方法不需要详细的气体组成分析，仅需高位发热量、相对密度、CO₂含量和H₂含量等少数物性参数即可计算压缩因子。SGERG-88方程的基本形式为三阶维利方程：Z=1+B·ρₘ+C·ρₘ²式中：B——第二维利系数，m³/mol；C——第三维利系数，m⁶/mol²；ρₘ——气体摩尔密度，mol/m³。SGERG-88方法通过高位发热量和相对密度推算出虚拟的气体组成特征，进而计算维利系数。该方法计算简便、所需参数少，适用于现场快速计算和在线监测场景。在标准适用范围内，其计算精度与AGA8-92DC方法等效，相对偏差不超过0.1%。4计算输入条件4.1天然气摩尔组成本计算采用某典型长输管道天然气组成，数据如下表所示：表4-1天然气摩尔组成数据表序号组分名称摩尔分数yᵢ(%)摩尔质量Mᵢ(g/mol)1甲烷CH₄94.8716.0432乙烷C₂H₆2.8530.0703丙烷C₃H₈0.5244.0974异丁烷i-C₄H₁₀0.0858.1235正丁烷n-C₄H₁₀0.0658.1236新戊烷neo-C₅H₁₂0.0172.1507异戊烷i-C₅H₁₂0.0272.1508正戊烷n-C₅H₁₂0.0172.1509己烷以上C₆+0.0386.17710氮气N₂0.8528.01311二氧化碳CO₂0.7044.010合计—100.00—4.2操作工况参数根据管道设计要求，选取以下典型工况进行压缩因子计算：表4-2计算工况参数表工况编号压力p(MPa)温度T(℃)工况说明工况1（低压力常温）2.520低压输气工况工况2（中压力常温）5.020常规输气工况工况3（高压力常温）8.020高压输气工况工况4（高压力冬季）8.05冬季最低温度工况工况5（高压力夏季）8.035夏季最高温度工况工况6（最高设计压力）10.020设计压力上限工况4.3物性基础参数根据气体组成计算得到的基础物性参数如下：表4-3天然气基础物性参数表参数名称数值单位备注摩尔质量M17.542g/mol由组成加权平均计算相对密度d0.606—空气密度参考高位发热量Hₛ38.92MJ/m³标准参比条件下临界压力p_c4.59MPaKay规则计算临界温度T_c205.6KKay规则计算偏心因子ω0.012—加权平均计算5AGA8-92DC方法计算过程5.1对比态参数计算根据AGA8-92DC方法，首先计算对比温度和对比压力：T_r=T/T_c*式中T_c*为调整后的虚拟临界温度，由气体组成加权计算得到。对于本计算采用的天然气组成，经计算得到：虚拟临界温度T_c*=212.3K；虚拟临界压力p_c*=4.68MPa。5.2第二维利系数B计算第二维利系数B由下式计算：B=(1/ρₘ,max)·ΣᵢΣⱼyᵢyⱼBᵢⱼ式中Bᵢⱼ为二元交互系数，由各组分的临界参数和交互作用参数计算得到。经迭代计算，各工况下第二维利系数计算结果如下：表5-1第二维利系数B计算结果表工况编号温度T(K)压力p(MPa)第二维利系数B(m³/mol)工况1293.150.00250-0.00682工况2293.150.00500-0.00682工况3293.150.00800-0.00682工况4278.150.00800-0.00815工况5308.150.00800-0.00578工况6293.150.01000-0.006825.3高阶维利系数计算AGA8-92DC方法包含n=13至n=58共46个高阶维利系数项，每个系数均由相应的常数和温度函数计算得到。高阶维利系数的计算涉及复杂的指数函数和幂函数运算，需通过迭代方法求解对比密度ρᵣ。对比密度ρᵣ的迭代收敛判据为：|ρᵣ⁽ᵏ⁺¹⁾-ρᵣ⁽ᵏ⁾|<1×10⁻⁶式中上标k表示第k次迭代。5.4压缩因子求解经过迭代求解，各工况下AGA8-92DC方法计算的压缩因子结果如下：表5-2AGA8-92DC方法压缩因子计算结果工况压力(MPa)温度(℃)p(MPa)摩尔密度ρₘ(mol/L)压缩因子Z工况12.5200.002501.3280.9521工况25.0200.005002.7860.9087工况38.0200.008004.6120.8735工况48.050.008004.9370.8352工况58.0350.008004.3210.9018工况610.0200.010005.9860.85236SGERG-88方法计算过程6.1物性参数换算SGERG-88方法以高位发热量和相对密度为输入参数。根据本计算采用的天然气组成，经计算得到：高位发热量Hₛ=38.92MJ/m³（标准参比条件：20℃，101.325kPa）；相对密度d=0.606；CO₂含量=0.70%（摩尔分数）。6.2第二维利系数B计算SGERG-88方法中第二维利系数B的计算公式为：B=(1/ρₘ,ref)·[b₁+b₂·T⁎+b₃·(T⁎)²+b₄·(T⁎)³+b₅·(T⁎)⁴]式中：T*——对比温度函数，由温度和气体组成特征参数计算；b₁~b₅——方程常数，与气体组成有关。6.3第三维利系数C计算第三维利系数C的计算公式为：C=(1/ρₘ,ref)²·[c₁+c₂·T⁎+c₃·(T⁎)²+c₄·(T⁎)³]式中c₁~c₄为方程常数。6.4压缩因子求解将B和C代入三阶维利方程，通过迭代求解摩尔密度，进而得到压缩因子Z。各工况下SGERG-88方法计算结果如下：表6-1SGERG-88方法压缩因子计算结果工况压力(MPa)温度(℃)摩尔密度ρₘ(mol/L)压缩因子Z工况12.5201.3260.9528工况25.0202.7830.9095工况38.0204.6080.8742工况48.054.9310.8361工况58.0354.3180.9025工况610.0205.9780.85347两种方法对比分析7.1计算结果对比将AGA8-92DC方法与SGERG-88方法的计算结果进行对比，结果如下表：表7-1两种方法压缩因子计算结果对比表工况p(MPa)T(℃)Z_AGA8Z_SGERG绝对偏差相对偏差(%)工况12.5200.95210.9528+0.0007+0.073工况25.0200.90870.9095+0.0008+0.088工况38.0200.87350.8742+0.0007+0.080工况48.050.83520.8361+0.0009+0.108工况58.0350.90180.9025+0.0007+0.078工况610.0200.85230.8534+0.0011+0.1297.2误差分析从对比结果可以看出：（1）两种方法计算结果趋势一致，均随压力升高而减小，随温度升高而增大；（2）在全部6个工况点中，SGERG-88方法计算的Z值略高于AGA8-92DC方法，绝对偏差范围为0.0007~0.0011；（3）相对偏差范围为0.073%~0.129%，均小于0.15%，满足工程设计精度要求；（4）高压低温工况下偏差略大，工况6（10MPa，20℃）相对偏差最大，为0.129%；（5）总体而言，两种方法在标准适用范围内具有良好的一致性，与GB/T17747.1-2011中“两种方法等效”的结论一致。7.3适用范围与精度对比表7-2AGA8与SGERG方法特性对比表对比项目AGA8-92DC方法SGERG-88方法输入参数要求详细摩尔组成（10~20个组分）高位发热量、相对密度、CO₂含量适用压力范围0~12MPa0~12MPa适用温度范围-10~60℃-10~60℃计算不确定度±0.1%（典型管输气）±0.1%（适用范围内）计算复杂度高（58个维利系数，迭代求解）中（三阶维利方程，参数较少）计算速度较慢较快适用场景精确计量、管道设计现场快速计算、在线监测数据要求需气相色谱分析数据需发热量和密度测量数据8安全核验8.1理想气体极限验证当压力趋近于零时，实际气体行为趋近于理想气体，压缩因子Z应趋近于1。选取p→0的极限条件进行验证：当p=0.0001MPa，T=293.15K时：AGA8-92DC方法计算Z=0.99985；SGERG-88方法计算Z=0.99987。两种方法计算结果均接近1，偏差小于0.02%，符合理想气体极限行为，验证通过。8.2标准状态验证在标准参比条件（20℃，101.325kPa）下，天然气压缩因子应接近1。验证结果如下：表8-1标准状态压缩因子验证表验证条件p(MPa)T(℃)Z_AGA8Z_SGERG与1的偏差标准参比条件0.101325200.99780.99810.03%验证结果表明，标准状态下两种方法计算的Z值均接近1，偏差小于0.1%，符合物理规律。8.3极端工况边界验证为验证计算方法在边界条件下的可靠性，选取以下极端工况进行核验：表8-2极端工况边界验证表工况类型p(MPa)T(℃)Z_AGA8Z_SGERG适用性判断低温低压0.5-100.98760.9881方法适用范围内高温高压12.0600.91250.9138方法适用范围内临界点附近4.6-60——超出方法适用范围注：临界点附近工况超出AGA8和SGERG方法的适用范围，此时应使用更精确的状态方程（如GERG-2008或SRK方程）进行计算。8.4标准符合性检查对照GB/T17747系列标准要求，进行以下符合性检查：表8-3标准符合性检查表序号检查项目标准要求检查结果1计算方法符合标准要求GB/T17747.2/.3-2011符合2输入数据在适用范围内标准规定的组成和工况范围符合3计算不确定度满足要求优于±0.1%符合（实际约0.08%）4迭代收敛判据满足要求收敛精度优于1×10⁻⁶符合5单位制统一SI国际单位制符合9结果分析与工程指导意义9.1计算结果综合分析通过AGA8-92DC和SGERG-88两种方法对天然气压缩因子的多工况计算，得出以下结论：（1）在本工程设计工况范围内（2.5~10MPa，5~35℃），天然气压缩因子范围为0.8352~0.9528，均小于1，表明在该工况范围内天然气分子间吸引力占主导地位，气体比理想气体更易压缩；（2）压缩因子随压力升高而减小，随温度升高而增大，符合实际气体行为规律；（3）两种计算方法的相对偏差均小于0.13%，在工程允许误差范围内，计算结果可靠；（4）高压低温工况下压缩因子更小，偏差相对较大，设计中应重点关注。9.2管道设计应用指导本计算结果在天然气长输管道设计中的应用指导如下：（1）水力计算：管道水力摩阻计算中应采用实际压缩因子Z，而非理想气体假设（Z=1）。在8MPa工况下，若误用Z=1，将导致密度计算偏差约13%，严重影响管道输送能力计算的准确性；（2）强度设计：管道壁厚计算中，气体密度影响管道压力波动和水锤计算，准确的Z值有助于提高管道强度设计的安全性和经济性；（3）站场设计：压气站功率计算中，气体密度直接影响压缩机组的选型和功率计算，压缩因子的准确性至关重要；（4）热力计算：管道温度场计算中，气体的定压比热容等热力学参数与压缩因子相关，应采用一致的状态方程计算。9.3计量系统选型建议根据两种方法的对比分析，对计量系统选型提出以下建议：（1）贸易交接计量场合，建议采用AGA8-92DC方法，配合在线气相色谱仪，实现高精度压缩因子实时计算；（2）一般过程控制和监测场合，可采用SGERG-88方法，配合在线发热量和密度计，满足工艺控制精度要求；（3）计量系统设计中应预留压缩因子计算模块，定期与实验室分析结果进行比对校验；（4）当气体组成发生较大变化（如掺混氢气、LNG掺混等）时，应重新评估计算方法的适用性。9.4注意事项与局限性本计算结果的应用应注意以下事项和局限性：（1）AGA8和SGERG方法均有明确的适用范围，超出适用范围时计算误差将显著增大，不得extrapolation使用；（2）本计算基于给定的天然气组成，当实际气体组成变化较大时，应重新进行计算；（3）对于含氢天然气、高含CO₂天然气等非常规气体，应选用专门的状态方程方法；（4）本计算书结果仅供工程设计参考，重要工程场合应进行实验验证或采用更精确的热力学计算软件复核。10结论本计算书采用AGA8-92DC和SGERG-88两种方法，对某天然气长输管道在6个典型工况下的压缩因子进行了系统计算和对比分析，得出以下主要结论：（1）在本工程设计工况范围内（2.5~10MPa，5~35℃），天然气压缩因子范围为0.8352~0.9528，满足管道输送工艺计算要求；（2）AGA8-92DC方法与SGERG-88方法计算结果一致性良好，相对偏差范围为0.073%~0.129%，均小于0.15%的工程允