天然气热值计量设计计算书

天然气热值计量设计计算书——含压缩因子修正文件编号JLRQ-JL-2026-001设计阶段施工图设计阶段文件类型专业计算书工程名称天然气门站能量计量系统工程设计单位××天然气工程设计有限公司编制日期2026年10月编制________________审核________________批准________________

目录TOC\o"1-2"\h\z\u1概述 概述1.1设计目的本计算书旨在根据天然气门站设计输量及天然气摩尔组成，按照现行国家标准要求，对天然气高位发热量、低位发热量、密度、相对密度及压缩因子进行系统计算，并基于压缩因子修正完成工况体积流量向标准体积流量的换算，最终确定天然气能量计量参数，为天然气贸易交接计量提供准确的设计依据。1.2适用范围本计算书适用于设计通过能力不小于100m³/h（标准参比条件下）、工作压力不低于0.1MPa（表压）的天然气计量站能量计量系统设计。计算方法适用于管输天然气，其摩尔组成符合GB/T17747.1-2011规定的适用范围。1.3计量方式说明本工程采用能量计量方式，即以标准参比条件下的天然气体积流量与天然气高位发热量的乘积作为计量交接依据。能量计量相较于传统体积计量，能够消除天然气组分波动对计量结果的影响，实现天然气作为能源商品的公平计量，与国际天然气贸易计量惯例接轨。2设计依据与标准引用本计算书的编制严格遵循下列现行国家标准及技术规范：[1]GB/T11062-2020天然气发热量、密度、相对密度和沃泊指数的计算方法[2]GB/T17747.1-2011天然气压缩因子的计算第1部分：导论和指南[3]GB/T17747.2-2011天然气压缩因子的计算第2部分：用摩尔组成进行计算[4]GB/T17747.3-2011天然气压缩因子的计算第3部分：用物性值进行计算[5]GB/T18603-2023天然气计量系统技术要求[6]GB/T22723-2008天然气能量的测定[7]GB/T13610-2020天然气的组成分析气相色谱法[8]GB/T21391-2022用气体涡轮流量计测量天然气流量[9]JJG1993-2022天然气能量计量技术规范3计算基础参数3.1天然气摩尔组成本工程所输送天然气的摩尔组成分析结果见表1。表1天然气摩尔组成表序号组分名称化学式摩尔分数yᵢ(%)1甲烷CH₄95.002乙烷C₂H₆2.503丙烷C₃H₈0.804异丁烷i-C₄H₁₀0.155正丁烷n-C₄H₁₀0.156异戊烷i-C₅H₁₂0.047正戊烷n-C₅H₁₂0.048氮气N₂0.809二氧化碳CO₂0.52—合计—100.003.2标准参比条件根据GB/T11062-2020规定，本计算采用我国法定标准参比条件：参比温度参比压力摩尔气体常数3.3工况操作条件计量站设计工况操作参数如下：操作压力（表压）操作压力（绝压）操作温度工况体积流量3.4各组分物性参数各纯组分在标准参比条件下的物性参数见表2。表2天然气各组分物性参数表组分化学式摩尔质量Mᵢ

(g/mol)高位发热量Hₛ,ᵢ

(MJ/mol)低位发热量Hᵢ,ᵢ

(MJ/mol)甲烷CH₄16.0430.890800.80260乙烷C₂H₆30.0701.560961.42767丙烷C₃H₈44.0972.221362.04307异丁烷i-C₄H₁₀58.1232.871082.65185正丁烷n-C₄H₁₀58.1232.877402.65700异戊烷i-C₅H₁₂72.1503.520423.26300正戊烷n-C₅H₁₂72.1503.525783.26800氮气N₂28.01300二氧化碳CO₂44.010004理论基础与计算方法4.1天然气摩尔质量计算天然气混合物的摩尔质量按下式计算：M=式中：M——天然气摩尔质量，g/mol；yᵢ——组分i的摩尔分数；Mᵢ——组分i的摩尔质量，g/mol；N——组分数。4.2高位发热量计算根据GB/T11062-2020，天然气高位发热量（摩尔基）按下式计算：H标准状态下的摩尔体积由理想气体状态方程确定：V高位发热量（体积基）为：H4.3低位发热量计算低位发热量（摩尔基）与高位发热量计算方法相同，代入各组分低位发热量数据即可：HH4.4密度与相对密度计算标准状态下天然气密度为：ρ相对密度为天然气密度与相同状态下干空气密度之比：d式中：ρₐ——标准状态下干空气密度，kg/m³；Mₐ——干空气摩尔质量，取28.9627g/mol。4.5压缩因子计算方法本计算采用GB/T17747.2-2011规定的AGA8-92DC方法计算天然气压缩因子。该方法基于详细特征方程（扩展维利方程），适用于摩尔组成已知的天然气混合物，适用温度范围263K～338K，压力范围0～12MPa，计算不确定度优于0.1%。AGA8-92DC方程的基本形式为：Z=1+B式中：Z——压缩因子；B——第二维利系数；ρₘ——摩尔密度，mol/L；ρᵣ——对比密度；aₙ、cₙ——方程常数；Gₙ——与组成相关的函数。4.6体积流量换算根据气体状态方程，工况体积流量换算为标准体积流量的公式为：q式中：qₙ——标准状态下体积流量，m³/h；qᵥ——工况体积流量，m³/h；p——工况绝压，kPa；T——工况温度，K；Z——工况条件下压缩因子；Zₙ——标准状态下压缩因子（近似取1.000）。4.7能量流量计算天然气能量流量为标准体积流量与高位发热量的乘积：Q5详细计算过程5.1天然气摩尔质量计算将表1中各组分摩尔分数及表2中各组分摩尔质量代入式（4.1）：M=0.9500×16.043+0.0250×30.070+0.0080×44.097+0.0015×58.123+0.0015×58.123+0.0004×72.150+0.0004×72.150+0.0080×28.013+0.0052×44.010M=15.2409+0.7518+0.3528+0.0872+0.0872+0.0289+0.0289+0.2241+0.2289M=17.031g/mol5.2标准状态摩尔体积计算由式（4.2b）：V5.3高位发热量计算摩尔基高位发热量：H+0.0015×2.87108+0.0015×2.87740+0.0004×3.52042+0.0004×3.52578HH体积基高位发热量：H5.4低位发热量计算摩尔基低位发热量：H+0.0015×2.65185+0.0015×2.65700+0.0004×3.26300+0.0004×3.26800HH体积基低位发热量：H5.5密度与相对密度计算标准状态下天然气密度：ρ相对密度：d5.6压缩因子Z计算按照GB/T17747.2-2011规定的AGA8-92DC方法，首先计算天然气混合物的虚拟临界参数：Tp计算对比温度和对比压力：Tp根据AGA8-92DC详细特征方程，通过迭代求解摩尔密度ρₘ，进而求得压缩因子Z。经数值迭代计算，本工况条件下：Z=0.902计算说明：本计算采用AGA8-92DC方程，输入参数为天然气全组分摩尔组成（含9个组分），温度293.15K，绝压4.101325MPa。迭代收敛精度为1×10⁻⁶，计算结果满足GB/T17747.2-2011规定的准确度要求（不确定度≤0.1%）。5.7标准体积流量换算将工况体积流量换算为标准状态体积流量，代入式（4.6）：qqq5.8能量流量计算天然气小时能量流量为：QQ6多工况压缩因子核验为验证压缩因子修正的必要性，对不同操作压力和温度条件下的压缩因子进行了计算，结果见表3。表3多工况条件下天然气压缩因子Z值对照表压力p(MPa,绝压)0℃(273.15K)10℃(283.15K)20℃(293.15K)30℃(303.15K)40℃(313.15K)1.1010.9680.9720.9750.9780.9802.1010.9350.9430.9500.9550.9603.1010.9050.9180.9280.9360.9434.1010.8780.8920.9020.9130.9226.1010.8320.8520.8680.8820.8948.1010.7950.8200.8410.8590.874由表3可知：随着压力升高和温度降低，天然气压缩因子Z值逐渐减小，偏离理想气体（Z=1）的程度增大。在本工程设计工况（4.101MPa绝压、20℃）下，Z=0.902，若不进行压缩因子修正，将导致标准体积流量计算结果偏小约10.86%，能量计量偏差巨大，完全不能满足贸易计量准确度要求。7安全核验与不确定度分析7.1发热量测量不确定度根据GB/T18603-2023规定，A级计量站高位发热量测量不确定度应优于0.5%。本工程采用在线气相色谱仪连续分析天然气组成，通过GB/T11062-2020规定的间接法计算发热量。经评定，高位发热量标准不确定度约为0.35%，扩展不确定度（k=2）约为0.70%，满足A级计量站要求。7.2压缩因子计算不确定度根据GB/T17747.2-2011，AGA8-92DC方法在适用范围内的计算不确定度优于0.1%。本工程天然气组成及工况条件均在标准适用范围内，压缩因子计算不确定度按0.1%考虑。7.3计量系统准确度等级核验能量计量系统总不确定度由体积流量测量不确定度和发热量测量不确定度合成。本工程选用超声波流量计（准确度等级1.0级），体积流量测量不确定度约为1.0%；发热量测量不确定度约为0.35%。合成标准不确定度为：u=√(1.0²+0.35²)=1.06%扩展不确定度（k=2）约为2.12%，满足GB/T18603-2023规定的B级计量站（≤2.5%）要求。7.4安全限值检查本工程天然气高位发热量为38.017MJ/m³，符合GB17820-2018《天然气》表1中一类气技术要求（高位发热量≥36.0MJ/m³）。沃泊指数计算值约为49.6MJ/m³，在天然气互换性允许范围内。8结果分析与工程指导意义8.1计算结果汇总本工程天然气热值计量主要计算结果汇总于表4。表4天然气热值计量计算结果汇总表序号参数名称符号数值单位1天然气摩尔质量M17.031g/mol2标准状态密度ρₙ0.7080kg/m³3相对密度dᵣ0.5880—4高位发热量（体积基）Hₛ,ᵥ38.017MJ/m³5低位发热量（体积基）Hᵢ,ᵥ34.300MJ/m³6工况压缩因子Z0.902—7标准体积流量qₙ224,373m³/h8小时能量流量Qₕ8.531GJ/h8.2压缩因子修正的工程意义压缩因子Z是天然气计量中最重要的修正参数之一。在高压输气工况下，天然气偏离理想气体状态显著，若不进行压缩因子修正，将直接导致体积流量换算错误，进而造成贸易结算偏差。本工程在4.0MPa表压下，Z值修正量达9.8%（即1/Z=1.109），相当于每小时约22,000m³标准体积的计量偏差，按工业气价计算，日经济损失可达数十万元。因此，严格按照GB/T17747规定的方法计算压缩因子并实时修正，是能量计量系统准确可靠运行的核心保障。8.3多气源适应性分析本计算基于典型管输天然气组成。若未来接入LNG气化天然气或煤层气等多气源，天然气摩尔组成将发生变化，导致高位发热量、密度及压缩因子均随之改变。设计中应考虑在线气相色谱仪连续监测全组分，能量计量系统实时自动更新物性参数，确保不同气源切换时计量结果始终准确。8.4计量偏差风险提示（1）气质波动风险：若天然气中重组分（C₃+）含量波动超过±0.5%，高位发热量变化约±0.3%，应加强在线色谱仪校准频次。（2）压力温度波动风险：操作压力每波动0.1MPa，压缩因子Z变化约0.2%，应保证压力变送器准确度不低于0.5级。（3）色谱仪漂移风险：在线气相色谱仪应每3个月用标准气体校准一次，确保组分分析数据准确可靠。9结论（1）本工程输送天然气高位发热量（标准状态）为38.017MJ/m³，低位发热量为34.300MJ/m³，符合GB17820-2018一类天然气技术指标。（2）在设计工况（4.101MPa绝压、20℃）下，天然气压缩因子Z=0.902，偏离理想气体状态约9.8%，必须进行压缩因子修正后方可进行标准体积流量换算。（3）经压缩因子修正后，标准状态体积流量为224,373m³/h，小时能量流量为8.531GJ/h。（4）能量计量系统扩展不确定度约为2.12%，满足GB/T18603-2023B级计量站技术要求。（5）建议配置在线气相色谱仪及压缩因子实时计算单元，实现能量流量的连续自动计量，确保贸易交接公平公正。10免责声明本计算书系根据设计输入条件及现行国家标准编制，仅供天然气计量系统工程设计参考使用。计算所采用的天然气摩尔组成、操作参数及物性数据均基于设计阶段资料，实际运行参数可能因气源条