天然气计量温压补偿设计计算书

第50页共50页天然气计量温压补偿设计计算书——××工业园区天然气分输计量站差压式孔板流量计温压补偿系统文件编号：WX-RQ-2026-JS-008设计阶段：施工图设计文件类型：设计计算书专业类别：仪表自控／油气储运版本：A版册数：共1册第1册编制：＿＿＿＿＿＿日期：＿＿＿＿审核：＿＿＿＿＿＿日期：＿＿＿＿批准：＿＿＿＿＿＿日期：＿＿＿＿××石油天然气工程设计有限公司二〇二六年十月

目录TOC\z\o"1-3"\u\h1总则 41.1编制目的 41.2适用范围 41.3编制依据 51.4设计原则 61.5计算内容与流程 61.6主要符号 72设计基础条件 92.1站场概况与工艺流程 92.2天然气组成 92.3参比条件与工况条件 102.4设计流量范围 102.5计量管与节流装置基本参数 112.6计量仪表配置及精度 113天然气物性参数计算 133.1摩尔质量与相对密度 133.2标准参比条件下的压缩因子与密度 143.3工作条件下的压缩因子与密度 153.4高位发热量、低位发热量与沃泊指数 153.5等熵指数与动力黏度 163.6物性参数计算成果汇总 174温压补偿原理与数学模型 184.1温度、压力对天然气计量的影响机理 184.2理想气体偏差与超压缩因子 184.3温压补偿的数学模型 194.4差压式流量计的温压补偿模型 194.5节流件热膨胀修正 214.6补偿误差传递与灵敏度分析 214.7温压补偿实现方式 225节流装置设计计算 235.1计算流程与迭代方法 235.2开孔直径与直径比计算 245.3流出系数计算 245.4可膨胀性系数计算 255.5雷诺数校核 265.6各流量点特性计算 265.7永久压力损失与输气能耗 275.8节流装置设计计算成果汇总 286温压补偿系统配置与参数整定 296.1节流装置与取压系统 296.2压力测量回路 296.3温度测量回路 296.4差压测量回路 306.5流量计算机与补偿算法 306.6补偿参数整定 316.7补偿回路故障处理策略 316.8能量计量接口 327测量不确定度评定 337.1评定依据与数学模型 337.2各不确定度分量 337.3合成与扩展不确定度 347.4评定结果分析 358温压补偿效果与适应性分析 378.1补偿前后计量偏差对比 378.2补偿后残余偏差 388.3量程比与低流量适应性 388.4计量方式适应性比选 398.5经济性分析 408.6气候与工况适应性 419安全核验 429.1管道强度与压力等级核验 429.2节流装置与引压系统安全核验 439.3仪表防爆与电气安全 439.4可燃气体检测与安全联锁 449.5水合物与冻堵风险核验 449.6超压保护与放散系统 459.7安全核验结论 4510结果分析与工程指导意义 4610.1主要计算结果汇总 4610.2计算结果综合分析 4610.3对工程设计的指导意义 4710.4对采购与施工的指导意义 4810.5对运行维护与计量管理的指导意义 4811结论 50附录A主要计算数据汇总 52免责声明 541总则1.1编制目的××工业园区天然气分输计量站（以下简称“本站”）是园区中压燃气管网的起点站，承担园区12家工业用户及3座CNG加气站的天然气供应与贸易交接计量任务。站内计量系统采用标准孔板差压式流量计（以下简称“孔板流量计”），设计最大通过能力8000m³/h（标准参比条件），工作压力0.40MPa（表压）。天然气是可压缩流体，其体积随温度、压力的变化可达每开尔文0.34%、每千帕0.2%量级。本站地处内陆地区，冬季最低气温可达−10℃，夏季地表管网温度可达35℃，昼夜温差与季节温差叠加，若计量系统不进行实时温度、压力补偿（以下简称“温压补偿”），计量结果将随环境与工况波动产生显著的系统性偏差，直接影响供用双方的贸易结算公平性。为系统完成本站温压补偿的方案设计，编制本计算书，主要目的如下：（1）确定本站天然气计量系统温压补偿的数学模型、计算公式与取值原则；（2）完成标准孔板节流装置的全套几何设计计算，并校核其适用条件；（3）完成温度、压力、差压测量回路及流量计算机（流量积算仪）的配置计算与参数整定；（4）评定温压补偿后体积流量测量系统的合成标准不确定度与扩展不确定度；（5）开展温压补偿效果定量分析、安全核验与经济性分析，形成工程实施指导意见。1.2适用范围本计算书适用于本站新建天然气贸易交接计量系统，并适用于与本站气质、压力、流量条件相近的同类天然气计量系统的设计计算参考。计算所依据的适用条件如下。（1）介质条件：单相、清洁、无游离水与固体颗粒的管输天然气，气质符合GB17820—2018《天然气》二类气要求。（2）节流装置条件：取压方式为法兰取压的标准孔板节流装置，孔板开孔直径d≥12.5mm，测量管内径50mm≤D≤1000mm，直径比0.1≤β≤0.75。（3）流动条件：气流处于充分发展的紊流状态，管道雷诺数Re_D≥5000。（4）工况条件：工作压力0.25~0.60MPa（表压），工作温度0~35℃，设计点压力0.40MPa（表压），设计点温度20℃。（5）系统条件：计量系统按GB/T18603—2023《天然气计量系统技术要求》中区域性分输站（B级）计量系统的技术要求配置。超出上述条件时，本计算书的结论需重新核算后方可引用。1.3编制依据本计算书引用的标准、规范及技术文件见表1-1。表1-1设计计算主要依据的标准、规范及技术文件序号标准／文件编号标准／文件名称1GB/T21446—2025用标准孔板流量计测量天然气流量2GB/T18603—2023天然气计量系统技术要求3SY/T6143—2004用标准孔板流量计测量天然气流量4GB/T35186—2025天然气计量系统性能评价5GB/T17747.1—2011～GB/T17747.3—2011天然气压缩因子的计算6GB/T11062—2020天然气发热量、密度、相对密度和沃泊指数的计算方法7GB17820—2018天然气8GB/T13610—2020天然气的组成分析气相色谱法9JJG640—2016差压式流量计检定规程10JJF1059.1—2012测量不确定度评定与表示11JJF1993—2022天然气能量计量技术规范12GB50028—2006（2020年版）城镇燃气设计规范13GB50251—2015输气管道工程设计规范14GB55009—2021燃气工程项目规范15GB/T50493—2019石油化工可燃气体和有毒气体检测报警设计标准16GB50183—2004石油天然气工程设计防火规范17GB50093—2013自动化仪表工程施工及质量验收规范18SH/T3005—2016石油化工自动化仪表选型设计规范19GB/T50770—2013石油化工安全仪表系统设计规范20GB3836.1爆炸性环境第1部分：设备通用要求21GB/T18604用气体超声流量计测量天然气流量22GB/T21391用气体涡轮流量计测量天然气流量注：表中未注明日期的引用文件，其最新版本（包括所有的修改单）适用于本计算书；凡引用文件中提出的要求，均以现行有效版本为准。1.4设计原则（1）合法合规原则。温压补偿方式、节流装置设计、仪表配置与不确定度水平均满足GB/T21446—2025、GB/T18603—2023及国家计量法规的规定，可直接用于施工图设计与贸易交接计量。（2）准确可靠原则。补偿模型覆盖温度、压力、压缩因子与节流件热膨胀四个修正环节，流量测量系统扩展不确定度满足B级计量站的技术要求，并为后续升级为A级能量计量预留接口。（3）宽范围补偿原则。温压补偿的有效范围覆盖全年环境温度与管网压力波动区间，在最小流量点仍能保持不低于10%的差压信号利用率。（4）故障安全原则。温度、压力补偿回路任一路故障时，流量计算机给出明确的报警并按预置的安全策略运行，避免因补偿失效造成计量结果大幅偏离。（5）经济合理原则。在满足准确度要求的前提下，控制节流装置的永久压力损失，降低输气能耗；仪表选型兼顾性能与全生命周期成本。1.5计算内容与流程本计算书的计算内容与流程见表1-2。表1-2温压补偿设计计算内容与流程步骤计算内容主要依据输出结果1基础条件整理、气质与工况参数确定工艺专业条件、P&ID第2章2天然气物性参数计算（摩尔质量、密度、压缩因子、发热量等）GB/T11062、GB/T17747第3章3温压补偿数学模型建立与误差传递分析GB/T21446、GB/T18603第4章4节流装置几何设计计算（孔板孔径、直径比、流出系数、可膨胀性系数等）GB/T21446、ISO5167-2第5章5温压补偿仪表配置与参数整定SH/T3005、GB/T18603第6章6测量不确定度评定JJF1059.1、GB/T21446第7章7温压补偿效果与经济性分析——第8章8安全核验GB50028、GB50251等第9章9结果分析与工程指导意义——第10章1.6主要符号本计算书使用的主要符号及其含义、单位见表1-3，符号含义在全文中保持一致。表1-3主要符号及含义符号含义单位q标准参比条件下的体积流量m³/hq工作条件下的体积流量m³/hq质量流量kg/sp孔板上游侧取压口处绝对压力Pap孔板下游侧取压口处绝对压力Pap标准参比条件压力，取101.325kPaPaT孔板上游侧气流热力学温度KT标准参比条件温度，取293.15KKΔp孔板上下游侧差压Paβ直径比，β=无量纲d工作条件下孔板开孔直径mD工作条件下上游测量管内径mC流出系数无量纲ε可膨胀性系数无量纲R以管内径计算的雷诺数无量纲ρ工作条件下孔板上游侧天然气密度kg/m³ρ标准参比条件下天然气密度kg/m³Z工作条件下天然气压缩因子无量纲Z标准参比条件下天然气压缩因子无量纲B第二维里系数m³/molκ等熵指数无量纲μ动力黏度Pa·sM摩尔质量kg/kmolG相对密度（对空气）无量纲H高位体积发热量MJ/m³W沃泊指数MJ/m³E能量流量GJ/hu标准不确定度——U扩展不确定度（k=2）——2设计基础条件2.1站场概况与工艺流程本站位于××工业园区北侧，来气自上游次高压管网以0.60~0.80MPa（表压）进入站内，经调压计量橇降压、计量后以0.40MPa（表压）向园区中压管网外输。站内计量橇采用“一开一备”双路布置，单路设计通过能力8000m³/h（标准参比条件）。单路计量系统的主要工艺构成如下：上游汇管→进口切断球阀→过滤器（过滤精度20μm）→调压器（监控调压＋工作调压）→计量直管段→标准孔板节流装置→下游汇管→出口切断球阀。孔板上游、下游按GB/T21446—2025的要求设置足够长度的直管段，取压方式采用法兰取压，差压由高静压差压变送器测量，气流压力与温度分别由压力变送器和温度变送器测量，三路信号送入流量计算机完成温压补偿计算、流量积算与数据远传。2.2天然气组成本站天然气组成按上游气源提供的在线气相色谱分析数据（依据GB/T13610—2020）确定，摩尔分数经归一化处理后见表2-1。表2-1天然气组成（摩尔分数）序号组分分子式摩尔分数y摩尔质量Mi1甲烷CH₄0.925016.0432乙烷C₂H₆0.040030.0703丙烷C₃H₈0.015044.0974异丁烷i-C₄H₁₀0.003058.1235正丁烷n-C₄H₁₀0.003058.1236氮气N₂0.010028.0137二氧化碳CO₂0.004044.010—合计—1.0000—由表2-1可见，本站天然气以甲烷为主（92.50%），C₂及以上重烃组分合计6.10%，非烃组分（N₂+CO₂）合计1.40%，属典型管输富气。气质参数满足GB17820—2018对二类气的要求：高位发热量不小于31.4MJ/m³，CO₂摩尔分数不大于3.0%，总硫不大于100mg/m³。2.3参比条件与工况条件温度、压力补偿计算涉及两类状态：标准参比条件（用于流量计量与贸易结算）与工作条件（现场实际状态），两者取值及依据见表2-2。表2-2参比条件与工作条件类别参数符号取值依据／说明标准参比条件参比压力p101.325kPaGB/T21446—2025、GB/T18603—2023标准参比条件参比温度T293.15K（20℃）同上标准参比条件参比湿度—干气同上燃烧参比条件燃烧温度t20℃GB/T11062—2020工作条件（设计点）工作压力（表压）p0.40MPa工艺专业提供工作条件（设计点）工作绝压p0.501325MPap工作条件（设计点）工作温度t20℃工艺专业提供工作条件（设计点）热力学温度T293.15KT工作条件（波动范围）压力波动p0.25～0.60MPa管网运行压力区间工作条件（波动范围）温度波动t0～35℃全年环境与输气温度区间大气条件当地大气压p101.325kPa按标准大气压计2.4设计流量范围按园区用气负荷特性与调峰要求，本站单路计量系统的设计流量范围见表2-3。流量量程比按最大流量与最小流量之比计。表2-3设计流量范围流量工况符号标准参比条件下体积流量qn/占最大流量比例最大流量q8000100%常用流量q500062.5%最小流量q250031.3%量程比—3.2∶1—差压式流量计的差压与流量平方成正比，量程比3.2∶1对应最小流量点差压为最大差压的9.8%，接近10%的工程下限，满足GB/T21446—2025对差压信号利用率的要求。若后续用气量进一步下降，应通过更换孔板或减小差压变送器量程的方式保证计量精度。2.5计量管与节流装置基本参数计量管与节流装置的设计基础参数见表2-4。表2-4计量管与节流装置基本参数参数名称符号取值说明管道规格—φ219×720#无缝钢管，GB/T8163测量管内径（20℃）D205.0mm按实测内径取值管道绝对粗糙度上限K≤0.001内壁无锈蚀、无结垢取压方式—法兰取压取压口距孔板端面25.4mm差压上限Δ16kPa差压变送器量程上限孔板材质—06Cr19Ni10耐蚀不锈钢，线膨胀系数16.6×10管道线膨胀系数λ11.6×10碳钢上游直管段长度L≥20D4100mm下游直管段长度L≥10D2050mm孔板两侧取压口间距—2×25.4mm法兰取压标准结构2.6计量仪表配置及精度计量系统主要仪表的测量范围与准确度等级见表2-5，仪表准确度等级的选择以满足B级计量站的技术要求并留有裕量为原则。表2-5计量系统主要仪表配置及精度序号仪表名称型号／规格测量范围准确度等级1标准孔板节流装置法兰取压、d20Δp≤16kPa几何尺寸符合GB/T214462差压变送器高静压智能型0～16kPa±0.075%FS3压力变送器智能型0～1.0MPa±0.1%FS4温度变送器（Pt100）铠装铂电阻−50～150℃±0.3℃（含变送器）5流量计算机多通道流量积算仪—计算误差优于0.05%6在线气相色谱仪用于组分与发热量分析—满足GB/T136107可燃气体探测器催化燃烧式／红外式0～100%LEL±3%LEL3天然气物性参数计算天然气物性参数是温压补偿计算的基础，本节按GB/T11062—2020与GB/T17747.1—2011～GB/T17747.3—2011给出的方法，逐项计算本站天然气的摩尔质量、相对密度、密度、压缩因子、发热量、等熵指数与动力黏度。3.1摩尔质量与相对密度天然气为多组分混合物，其摩尔质量按各组分摩尔分数加权计算： M=i=1Nyi式中：M为天然气摩尔质量，kg/kmol；yi为组分i的摩尔分数；Mi为组分i的摩尔质量，kg/kmol；N代入表2-1数据，各组分对摩尔质量的贡献见表3-1。表3-1各组分对摩尔质量的贡献序号组分yMiyi1CH₄0.925016.04314.83982C₂H₆0.040030.0701.20283C₃H₈0.015044.0970.66154i-C₄H₁₀0.003058.1230.17445n-C₄H₁₀0.003058.1230.17446N₂0.010028.0130.28017CO₂0.004044.0100.1760—合计1.0000—17.5090计算得天然气摩尔质量M=17.509kg/kmol。相对密度为天然气密度与同温同压干空气密度之比，按摩尔质量之比计算： G=MMair 式中：G为相对密度；Mair为干空气摩尔质量，取28.9626kg/kmol（GB/T11062—2020代入得： G=17.50928.9626=0.6045 相对密度0.6045对应GB/T11062—2020中典型管输富气的取值范围，与气质组成相符。3.2标准参比条件下的压缩因子与密度天然气在计量压力下偏离理想气体行为，压缩因子Z是温压补偿中除压力、温度之外的第三个关键修正量。本站标准参比条件压力仅为101.325kPa，对比压力极低，本节采用第二维里系数法计算压缩因子，该方法在低压比条件下与GB/T17747.2—2011（AGA8详实法）的一致性优于5×10压缩因子的维里方程（截取第二维里项）为： Z=1+BTpRT式中：Z为压缩因子；BT为第二维里系数，m³/mol；p为绝对压力，Pa；T为热力学温度，K；R为通用气体常数，取8.314462J/(mol·K)混合物第二维里系数按二次混合规则计算： BT=i=1N纯组分及交叉第二维里系数采用Pitzer-Curl关联式计算： Bij=RT B0=0.083-0.422T B1=0.139-0.172T式中：Tc,ij、pc,ij、ωij分别为交叉临界温度、交叉临界压力与交叉偏心因子，按Tsonopoulos混合规则由纯组分临界参数确定；按上述各式计算，标准参比条件（293.15K、101.325kPa）下天然气的第二维里系数为B=-52.758cm³/mol，压缩因子为： Zn=1+-52.758×10标准参比条件下的天然气密度按实际气体状态方程计算： ρn=pn代入得： ρn=101325×17.5090.99781×8.314462×293.153.3工作条件下的压缩因子与密度工作条件（0.501325MPa、293.15K）下天然气对比压力约为0.11，第二维里截断法仍然适用。计算得工作条件下第二维里系数为B=-52.758cm³/mol（温度相同），压缩因子为： Z1=1+-52.758×10工作条件下的天然气密度为： ρ1=p1结果表明：工作条件下压缩因子为0.98915，偏离理想气体约1.09%；若在流量计算中忽略压缩因子修正（取Z1=1），将引入约3.4高位发热量、低位发热量与沃泊指数按GB/T11062—2020，天然气高位体积发热量由各组分理想气体体积发热量按摩尔分数加权求和： Hs=i=1N式中：Hs为高位体积发热量，MJ/m³；Hs,i为组分i的理想气体体积发热量（燃烧参比条件20℃、计量参比条件20℃、各组分取值及加权结果见表3-2。表3-2高位体积发热量计算序号组分yHs,iyi1CH₄0.925037.04434.2662C₂H₆0.040064.9102.5963C₃H₈0.015092.2901.3844i-C₄H₁₀0.0030121.0600.3635n-C₄H₁₀0.0030121.4200.3646N₂0.0100007CO₂0.004000—合计1.0000—38.974计算得高位体积发热量Hs=38.974 Hs,m=Hs低位体积发热量按各组分的低位／高位热值比加权修正，计算得Hl=34.864沃泊指数表征燃气互换性与能量输送能力，按下式计算： W=HsG=能量流量按体积流量与高位发热量的乘积计算： E=qnHs本站最大流量点与常用流量点的能量流量分别为311.79GJ/h与194.87GJ/h；按年运行8000h、平均负荷率60%估算，年交接能量约为1.50×1063.5等熵指数与动力黏度等熵指数κ用于可膨胀性系数计算，由各组分比热容按摩尔分数加权后求得： κ=CpCv计算得混合气定压比热容Cp=37.20J/(mol·K)，定容比热容Cv=28.89J/(mol·K)动力黏度μ影响雷诺数与流出系数，按工况条件取μ=1.08×10-53.6物性参数计算成果汇总天然气物性参数计算结果汇总见表3-3。表3-3天然气物性参数计算成果汇总序号参数名称符号单位数值1摩尔质量Mkg/kmol17.5092相对密度G—0.60453标准参比条件密度ρkg/m³0.72954工作条件密度ρkg/m³3.64085标准参比条件压缩因子Z—0.997816工作条件压缩因子Z—0.989157高位体积发热量HMJ/m³38.9748低位体积发热量HMJ/m³34.8649质量高位发热量HMJ/kg53.4310沃泊指数WMJ/m³50.1311等熵指数κ—1.28812动力黏度μPa·s1.08×4温压补偿原理与数学模型4.1温度、压力对天然气计量的影响机理天然气是可压缩流体，其体积随温度、压力显著变化。在工作压力接近常压、温度接近常温的工况下，气体的体积流量可近似按理想气体状态方程描述： qn=qv由该式可见，温度升高使气体膨胀、实际密度减小，若按固定参数换算将高估标准体积流量；压力升高使气体压缩、密度增大，若按固定参数换算将低估标准体积流量。对理想气体，温压补偿的误差满足如下与压力、温度成正比的关系： δqnqn对本站工作状态（p1=501325Pa、T1=293.15K）作灵敏度分析：温度每偏差1K，未补偿的标准体积流量偏差约0.34%；压力每偏差1kPa，未补偿的标准体积流量偏差约0.20%。本站全年温度波动区间0～35℃、压力波动区间4.2理想气体偏差与超压缩因子实际气体分子间存在相互作用力且分子本身占有体积，其行为偏离理想气体，这一偏离用压缩因子Z表征。在天然气计量中定义超压缩因子Fpv Fpv=1ZFpv大于1表示气体较理想气体更难压缩。本站工作条件下Fpv=压缩因子的计算方法按GB/T17747.1—2011～GB/T17747.3—2011分为两类：详实法（AGA8详实法）需完整摩尔组分，适用于贸易交接计量；简化法（SGERG-88等）只需相对密度、发热量及非烃组分含量，适用于趋势监控。本站配置在线气相色谱仪，具备完整组分数据，温压补偿计算中压缩因子按GB/T17747.2—2011详实法计算，并在低压比区间采用第二维里截断法进行快速核算，两者偏差不大于5×104.3温压补偿的数学模型天然气标准体积流量的温压补偿数学模型为： qn=qv定义温压补偿系数（又称工况转换系数）Kc Kc=p1则 qn=Kc式中：qv为流量计实测的工作条件体积流量，m³/h；p1、T1为工作条件绝对压力与热力学温度；pn、Tn为标准参比条件压力与温度；该模型包含四个修正环节：压力修正、温度修正、压缩因子修正以及节流件热膨胀修正（见4.5节）。前三项通过式（4-4）统一体现，第四项在孔板开孔直径与管内径的实时计算中体现。按设计点参数计算本站在标准参比条件的补偿系数： Kc=501325×293.15×0.99781101325×293.15×0.98915即1m³工作条件下的天然气对应4.989m³标准参比条件下的天然气（设计点）。4.4差压式流量计的温压补偿模型孔板流量计的流量计算以能量方程与连续性方程为基础，其质量流量方程（GB/T21446—2025）为： qm=C1-式中：qm为质量流量，kg/s；C为流出系数；ε为可膨胀性系数；β=dD为直径比；d为工作条件下孔板开孔直径，m；D为工作条件下上游测量管内径，m；Δp为差压，Pa；ρ1工作条件下天然气密度按下式计算： ρ1=p1标准参比条件下的体积流量为： qn=qm将密度表达式代入质量流量方程可见，压力与温度通过密度ρ1 qm∝p1即差压式流量计的压力、温度误差对流量结果的影响被“开平方”，其灵敏度约为工况体积流量型仪表的一半。差压式与工况体积流量式两类仪表的温压补偿模型对比见表4-1。表4-1两类流量计温压补偿模型对比比较项目差压式（孔板）流量计工况体积流量型（涡轮、超声等）测量量差压Δp工作条件体积流量q流量方程qqv压力灵敏度0.5（∂qn1.0温度灵敏度−0.5−1.0密度影响以平方根方式进入仅用于体积换算补偿精度要求相对宽松相对严格由上表可见，在相同的压力、温度测量准确度条件下，差压式流量计的温压补偿结果更为稳健，这是本站选择孔板流量计的重要技术依据之一。4.5节流件热膨胀修正孔板开孔直径与测量管内径随介质温度变化，当工作温度偏离孔板标定温度（通常为20℃）时，需按材料线膨胀系数进行修正： dt=d20 Dt=D20式中：dt、Dt分别为工作温度t1下的孔板开孔直径与管内径；d20、D20为20℃下的相应尺寸；λd、λD为孔板（不锈钢）与管道（碳钢）的线膨胀系数，分别取温度由20℃升高至35℃时，孔板开孔直径的相对变化为： Δdd=16.6×10-6即0.025%，对流量结果的影响约为0.05%（孔径以二次方进入流量方程）。温度由20℃降低至0℃时，影响方向相反、量值相近。热膨胀修正量虽小，但在贸易交接计量中不可忽略，本站流量计算机按实时温度自动完成该项修正。4.6补偿误差传递与灵敏度分析温压补偿的误差由各输入量的测量误差按误差传递规律合成。对式（4-4）取对数后微分，可得标准体积流量的相对误差传递关系： δqnqn对差压式流量计，考虑密度以平方根方式进入流量方程，其综合误差传递关系为： δqnqn对工况体积流量型流量计，压力、温度的弹性系数（灵敏度系数）S定义如下： Sp=∂q对差压式流量计，相应弹性系数为Sp=0.5、本站温压补偿系统的灵敏度分析结果见表4-2。该表表明：温度测量的1K误差在差压式系统中引起0.17%的流量偏差，压力测量1kPa误差引起约0.10%的流量偏差，因此温度、压力变送器的准确度等级必须与其对流量结果的贡献相匹配。表4-2温压补偿系统灵敏度分析输入量单位偏差差压式流量偏差工况体积流量式偏差说明压力p1kPa0.10%0.20%工作压力501.325kPa压力p1%0.50%1.00%相对偏差温度T1K0.17%0.34%工作温度293.15K温度T1%0.50%1.00%相对偏差压缩因子Z0.1%0.05%0.10%详实法计算精度4.7温压补偿实现方式本站采用“实时全补偿”方式，由流量计算机按采样周期读取工作压力、温度、差压与组分数据，逐周期完成下列补偿计算：（1）按实测温度修正孔板开孔直径与管内径（式4-12、式4-13）；（2）按实测压力、温度与在线组分计算工作条件下压缩因子与密度（式3-3～式3-13）；（3）按实测差压、密度与几何参数计算质量流量（式4-8）；（4）按标准参比条件下的密度换算标准体积流量（式4-10），并累计交接气量与能量流量。当在线色谱仪故障或组分数据中断时，流量计算机自动切换为按固定组分（按最近一次有效分析值锁定）计算压缩因子，同时发出报警并经通信接口上传，避免因组分缺失导致计量中断。5节流装置设计计算5.1计算流程与迭代方法标准孔板的设计计算不能由流量方程直接解出开孔直径。原因在于流量方程中的流出系数C与可膨胀性系数ε均与直径比β、雷诺数ReD相关，而ReD又取决于流量与管内径，因而孔径（1）假定初始直径比β0，由d=βD（2）由管内径计算雷诺数Re（3）按Reader-Harris/Gallagher公式（式5-5～式5-7）计算流出系数C，按式5-10计算可膨胀性系数ε；（4）由流量方程反算新的开孔直径d，并更新β；（5）比较新的β与上一次迭代值，若差值大于收敛容差则重复步骤（2）～（4）。迭代收敛判据取βk设计条件为：最大标准体积流量qn,max=8000m³/h，最大差压Δpmax=16kPa，标准参比条件密度ρn=0.7295kg/m³，工作条件密度ρ1=3.6408kg/m³，测量管内径D=0.205m，工作压力对应最大流量的质量流量为： qm=qn,由图5-1的迭代流程，反算开孔直径的迭代公式为： d=qm1-β5.2开孔直径与直径比计算按5.1节迭代方法计算，迭代过程见表5-1。计算表明，经5次迭代后直径比变化量已小于1×10表5-1孔板开孔直径迭代计算过程迭代序号假定直径比β计算开孔直径d/mm雷诺数R流出系数C可膨胀性系数ε回代所得β'β'-β0（初值）0.500000102.500—————10.500000102.5009322780.603210.990780.482831.72×20.48282798.9809322780.602740.990860.484041.22×30.48404399.2299322780.602770.990850.483968.25×40.48396099.2129322780.602770.990850.483975.61×50.48396599.2139322780.602770.99085—1×10迭代收敛后，孔板在20℃下的开孔直径取整为： d20=99.21mm （相应的直径比为： β=d20D20直径比0.48396位于标准推荐的0.4～0.6区间内，既可保证计量准确度，又能有效控制永久压力损失，符合GB/T21446—2025与SY/T6143—2004的建议。5.3流出系数计算流出系数按ISO5167-2:2022与GB/T21446—2025采用的Reader-Harris/Gallagher公式计算： C=0.5961+0.0261β2-0.216β +0.043+0.080e-10L -0.031M2'-0.8M式中： A=19000βReD M2'=2L法兰取压时，L1=L2=L2'=25.4将β=0.48396、ReD=932278、L1=L2=0.12390代入，中间量计算结果为：A=0.02484，M2'=0.48021，M2流出系数各项计算见表5-2。表5-2流出系数分项计算序号计算项表达式数值1基准项0.59610.5961002直径比平方项0.02610.0061133直径比八次方项-0.216−0.0006504雷诺数项0.0005210.0003295组合项Ⅰ0.0188+0.0063A0.0015276组合项Ⅱ0.043+0.0800.0008407取压口位置项-0.031−0.001487—合计C0.602772计算得流出系数C=0.60277。该值对应法兰取压标准孔板在β≈0.484、ReD≈9.3×1055.4可膨胀性系数计算天然气流经孔板时因压力降低而膨胀、密度减小，需引入可膨胀性系数修正。按ISO5167-2:2022给出的公式： ε=1-0.351+0.256β4+0.93式中：p2为孔板下游侧取压口处绝对压力，p2=p代入Δpmax=16kPa、p1=501325Pa，得p2p1=0.968085；代入 ε=1-0.36756×1-0.975128=0.99085 （计算得最大差压点可膨胀性系数ε=0.99085。由于本站工作压力相对差压较高（Δpmaxp1=3.19%），可膨胀性系数接近5.5雷诺数校核雷诺数是判断流动是否满足标准适用条件的关键参数，以管内径计算的雷诺数按下式核算： ReD=4代入qm=1.6211kg/s、D=0.205m、μ=1.08×10-5Pa·s校核结论：ReD=9.3228×105≫5000，满足GB/T21446—2025与ISO5167-2:2022对雷诺数下限的要求；同时β=0.48396处于0.1≤β≤0.75的适用区间，D=205mm处于50≤D≤1000mm5.6各流量点特性计算按孔板口径确定后，各流量点对应的差压与流动参数按下式换算（差压与流量的平方成正比）： Δp=Δpmaxqn工作条件下的体积流量与管内平均流速分别为： qv=qn v1=qv各流量点的计算结果汇总于表5-3。表5-3各流量点特性参数计算流量工况qnΔp/Pa差压占量程比qv管内流速v1RCε最大流量800016000100.0%1602.913.499.32×0.60280.9909常用流量5000625039.1%1001.88.435.83×0.60310.9964最小流量25001562.59.8%500.94.222.91×0.60390.9991由表5-3可见：管内流速在最大流量点达到13.49m/s，处于中压燃气管道的常用流速范围（10～20m/s）之内；最小流量点差压为量程的9.8%，接近10%的工程下限，差压信号仍具有足够的利用率；孔板开孔处流速为57.60m/s，未超过标准对气体流速的限制。5.7永久压力损失与输气能耗孔板节流装置的永久压力损失反映被测流体经过节流件后无法恢复的那部分压力，其常用计算式为： δp=1-β2Δp代入β=0.48396、Δpmax δpmax=1-0.234217即12.25kPa，占工作压力（501.325kPa绝压，401.325kPa表压）的3.05%（以表压计约3.05%）。最大流量点的压损功率按流量与压力损失的乘积计算： P=qvδp代入qv=1602.9m³/h0.4453m³/s，得P=5.46kW。按年运行8000h、平均负荷率0.6、电价0.7元/(kW·h)估算，年折算输气能耗费用约为1.83万元需要说明的是，差压式流量计的永久压力损失与其准确度是一对矛盾：减小β可提高差压信号利用率，但压损随之增大；增大β可降低压损，但会降低信噪比。本站选定的β=0.484是在两者之间取平衡的结果。若后续对能耗有更高要求，可考虑改用压力损失更小的超声流量计或涡轮流量计（其压损约为孔板的1/5～1/3），但需同步提高温压补偿精度，相关比选见8.4节。5.8节流装置设计计算成果汇总节流装置设计与校核成果汇总见表5-4。表5-4节流装置设计计算成果汇总序号项目符号单位计算结果标准要求／限值判定1上游测量管内径Dmm205.050～1000满足2孔板开孔直径dmm99.21≥12.5满足3直径比β—0.483960.1～0.75满足4取压方式——法兰取压法兰／角接取压满足5最大差压ΔPa16000由设计选定—6雷诺数R—9.32×≥5000满足7流出系数C—0.60277——8可膨胀性系数ε—0.99085——9永久压力损失δpPa12253宜小于工作压力5%满足10上游直管段长度Lmm≥4100≥20D满足11下游直管段长度Lmm≥2050≥10D满足12最小流量点差压比—%9.8≥10（工程建议）基本满足6温压补偿系统配置与参数整定6.1节流装置与取压系统节流装置采用法兰取压标准孔板，孔板材质为06Cr19Ni10不锈钢，取压口中心距孔板上游端面25.4mm。孔板安装于计量直管段中部，上游直管段长度不小于20D（4100mm），下游直管段长度不小于10D（2050mm），孔板上游2D范围内管道内壁无焊缝、无凸出物。引压管路采用φ14×2不锈钢管，取压口设于管道上部（气体介质），引压管以不小于1∶12的坡度向取压口倾斜，使冷凝液能够自行回流至主管道；引压管最高点设排气阀，最低点设排污阀。差压变送器与引压管整体保温，冬季投入电伴热，防止冷凝液结冰造成测量失真。6.2压力测量回路压力变送器量程按工作压力波动范围的1.3～1.5倍选取。本站工作压力范围为0.25～0.60MPa（表压），取量程0～1.0MPa（表压），精度等级±0.1%FS。压力变送器接液部分材质为316L，过程连接为M20×1.5外螺纹，防护等级不低于IP65，防爆等级不低于ExdⅡBT4Gb。压力测量对流量结果的相对不确定度贡献按下式估算： uqnqn代入up1=0.1%×1.0MPa1.0kPa、p1=501.325kPa，得相对贡献为6.3温度测量回路温度测量采用A级铂电阻（Pt100）配智能温度变送器，插入式安装于孔板下游侧、距孔板端面不小于5D处。测温元件套管材质为316L，插入深度不小于管道内径的1/3，且不小于75mm。温度变送器量程−50～150℃，精度±0.3℃（含变送器）。温度测量对流量结果的相对不确定度贡献按下式估算： uqnqn代入uT1=0.3K、T1=293.15K，得相对贡献为6.4差压测量回路差压变送器为高静压智能型，量程0～16kPa，精度±0.075%FS，额定静压4MPa，单向过载能力不低于4MPa。差压变送器安装在孔板下方，两根引压管对称布置、长度相等，避免附加差压。差压测量对流量结果的相对不确定度贡献为： uqnqn以常用流量点（Δp=6250Pa）为例，uΔp=0.075%×16000=12Pa，相对贡献为0.5×126250=0.096%6.5流量计算机与补偿算法流量计算机为多通道流量积算仪，按1s计算周期执行温压补偿与流量积算，主要功能配置如下。（1）参比条件设置：参比压力101.325kPa，参比温度20℃，与GB/T21446—2025及贸易合同约定一致。（2）物性计算模块：内置GB/T17747.2—2011（AGA8详实法）压缩因子计算程序，按在线气相色谱仪提供的摩尔组分实时计算Z1与Zn；组分数据刷新周期（3）热膨胀修正模块：按实测温度自动修正孔板开孔直径与管内径（式4-12、式4-13）。（4）流量计算模块：按式4-8、式4-10计算质量流量与标准体积流量，并按式3-17计算能量流量。（5）积算与记录模块：按小时、日、月累计标况体积流量与能量，掉电保持时间不少于10年；具备90d的分钟级历史数据存储能力。（6）通信模块：支持ModbusRTU/TCP与RS-485，向站控系统上传瞬时流量、累计流量、温度、压力、差压、压缩因子及报警状态。6.6补偿参数整定温压补偿相关参数的整定值见表6-1。表6-1温压补偿参数整定表序号参数符号／说明整定值1参比压力p101.325kPa2参比温度T20℃3管道内径（20℃）D205.0mm4孔板开孔直径（20℃）d99.21mm5管道线膨胀系数λ11.6×106孔板线膨胀系数λ16.6×107差压量程Δp0～16kPa8压力量程p0～1.0MPa（表压）9温度量程t−50～150℃10等熵指数κ1.288（按组分自动计算）11动力黏度μ1.08×10-512压缩因子计算方法—GB/T17747.2—2011详实法13计算周期—1s14组分刷新周期—5min6.7补偿回路故障处理策略温度、压力任一补偿回路失效都将导致计量结果严重失真。为使系统在故障状态下仍保持“可计量、可追溯”，本站流量计算机设置如下故障处理策略，见表6-2。表6-2温压补偿回路故障处理策略序号故障情形系统响应计量处理方式1压力信号超量程或断线声光报警、上传故障位切换至上一有效压力值，并标记该时段数据2温度信号超量程或断线声光报警、上传故障位切换至上一有效温度值，并标记该时段数据3差压信号异常声光报警停止积算，等待恢复4在线色谱仪故障报警按最近一次有效组分计算压缩因子5仪表电源失电报警、记录流量计算机由UPS供电，维持计量与记录6通信中断报警本地继续积算与存储，恢复后补传数据上述策略的核心思想是：任何单点故障都不应使计量值发生阶跃性突变。故障期间的数据应被明确标记，作为后续贸易结算协商的依据，避免以失真数据直接用于结算。6.8能量计量接口按GB/T18603—2023与JJF1993—2022的要求，本站计量系统预留能量计量功能：由在线气相色谱仪提供组分与高位发热量数据，流量计算机按式3-17计算瞬时与累计能量流量，并同时输出体积流量与能量流量两套数据。当前阶段以体积流量作为交接结算依据，能量数据用于用户能耗分析与后续计量方式升级的技术储备。7测量不确定度评定7.1评定依据与数学模型测量不确定度按JJF1059.1—2012《测量不确定度评定与表示》规定的方法评定，各输入量不确定度的取值依据GB/T21446—2025与相关仪表的技术条件。标准参比条件下体积流量的数学模型可表示为各输入量的函数： qn=fd,D,C,ε,Δp,ρ其中工作条件密度ρ1 ρ1=p1由于各输入量之间互不相关，合成标准不确定度按方和根法合成，标准体积流量的相对合成标准不确定度为： uqnqn工作条件密度的相对标准不确定度为： uρ1ρ1扩展不确定度取包含因子k=2： U=kuqn （7.2各不确定度分量各输入量的标准不确定度及其来源见表7-1。表7-1各不确定度分量评定序号分量不确定度来源量值相对标准不确定度对流量结果的影响1u孔板孔径加工公差与测量（±0.05mm）0.05mm0.05%二次方进入，贡献0.10%2u流出系数公式不确定度（GB/T21446）—0.50%直接进入，0.50%3u可膨胀性系数公式不确定度—0.10%直接进入，0.10%4u差压变送器±0.075%FS及安装影响12Pa0.19%（常用点）0.075%（最大点）平方根进入，0.096%5u压力变送器±0.1%FS1.0kPa0.20%经密度平方根进入6u温度变送器±0.3℃0.3K0.10%经密度平方根进入7u压缩因子（AGA8详实法）—0.10%经密度平方根进入8u组分分析（在线色谱）—0.05%经密度平方根进入9u标准参比条件密度（含发热量）—0.05%直接进入由表7-1中第5～8项合成工作条件密度的相对标准不确定度： uρ1ρ17.3合成与扩展不确定度以常用流量点（qn=5000m³/h，Δp=6250分量数值／%平方／%²2u0.100.0100u0.500.2500u0.100.01000.5u0.0960.00920.5u0.1250.0156u0.050.0025合计—0.2973由表计算得常用流量点标准体积流量的相对合成标准不确定度： urelqn=扩展不确定度为： Urelqn=2×0.545按同样方法计算不同流量点的相对扩展不确定度，结果见表7-2。表7-2不同流量点的测量不确定度流量工况qnΔp/PauuUrelqn最大流量8000160000.075%0.538%1.08%常用流量500062500.192%0.545%1.09%最小流量25001562.50.768%0.660%1.32%低于量程点1500562.52.133%1.194%2.39%7.4评定结果分析（1）结果水平。本站温压补偿后的体积流量测量系统扩展相对不确定度为1.08%～1.32%（k=2），在常用流量点约为1.09%，满足GB/T18603—2023对区域性分输站（B级）计量系统的准确度要求，可用于贸易交接计量。（2）主要误差源。由表7-1可见，流出系数C的不确定度（0.50%）在合成不确定度中占绝对主导地位，其平方项占总方差的84.1%；其次为压力与温度测量（经密度平方根传递）。这说明：提高系统准确度的关键是提高流出系数的可信度，而非盲目提高变送器精度等级。（3）校准优化的潜力。若对本站孔板流量计实施在线实流校准（标准表法），流出系数的相对不确定度可由0.50%降至0.30%左右，此时合成相对标准不确定度降为： urelqn= Urelqn=2×0.371即扩展不确定度由1.09%降至0.74%，提升幅度达32%，可满足A级交接站的技术要求。建议在投运后按JJG640—2016的规定开展首次检定，并视条件配置在线实流校准装置。（4）不确定度的工况依赖。由表7-2可见，随着流量下降，差压信号减弱，uΔpΔp迅速增大，扩展不确定度由最大流量点的1.08%上升至1500m³/h时的8温压补偿效果与适应性分析8.1补偿前后计量偏差对比为定量说明温压补偿的必要性，本节按“以设计点参数为固定输入”的方式计算不实施温压补偿时的计量偏差。偏差定义为： δ=qn,disp-q式中：qn,disp为按固定设计参数计算所得的示值；qn,act对差压式流量计，忽略压缩因子在固定参数条件下的差异，偏差主要来自密度偏差： δΔp=ρ1,des对工况体积流量型流量计（涡轮、超声等），偏差为： δv=p1,des各典型工况下的计算见表8-1。表8-1不实施温压补偿时的计量偏差序号实际工况（表压、温度）实际密度/(kg·m⁻³)与设计点密度偏差差压式偏差工况体积流量式偏差10.40MPa、20℃（设计点）3.640800020.40MPa、30℃3.5163−3.42%+1.75%+3.41%30.40MPa、5℃3.8456+5.62%−2.70%−5.12%40.30MPa、20℃2.9082−20.13%+11.89%+24.92%50.50MPa、20℃4.3766+20.21%−8.79%−16.63%60.30MPa、30℃2.8094−22.84%+13.84%+29.18%由表8-1可以得出三点重要结论：（1）温度影响不可忽略。仅温度由20℃升至30℃，不补偿的差压式系统即产生+1.75%的计量偏差，相当于每1K约0.17%；工况体积流量型仪表的偏差可达+3.41%。按本站年交接气量约2400万m³、气价2.5元/m³测算，1.75%的偏差对应约105万元的年度计量差异，远大于温压补偿装置的投入。（2）压力影响更为严重。压力由0.40MPa降至0.30MPa时，不补偿的差压式系统偏差高达+11.89%，工况体积流量型仪表偏差达+24.92%。中压管网压力随负荷波动属常态，若压力回路失效而未及时发现，将造成重大计量纠纷。（3）差压式流量计的偏差被“压缩”。由于密度以平方根方式进入流量方程，同样工况偏差下，差压式系统的计量偏差约为工况体积流量型的1/2（如表8-1中+1.75%与+3.41%、+11.89%与+24.92%），体现出差压式流量计在温压补偿方面的相对稳健性。8.2补偿后残余偏差实施实时温压补偿后，由温度、压力、压缩因子及热膨胀修正引入的残余偏差被压缩至测量不确定度所限定的范围内。以常用流量点为例，补偿后残余偏差的各分项贡献见表8-2。表8-2补偿后残余偏差的分项贡献（常用流量点）序号分项贡献方式相对偏差1温度测量误差（±0.3℃）经密度平方根传递0.051%2压力测量误差（±1.0kPa）经密度平方根传递0.100%3压缩因子计算误差经密度平方根传递0.050%4组分分析误差经密度平方根传递0.025%5热膨胀修正残差（温度±0.3℃）孔径二次方传递0.010%—补偿环节合计方和根合成0.122%补偿环节引入的残余偏差合计仅0.122%，与流出系数不确定度（0.50%）相比很小，说明本站温压补偿方案的仪表配置是合理且有余量的，整体不确定度由节流装置本身（流出系数）而非补偿环节决定。8.3量程比与低流量适应性差压式流量计的根本特点是差压与流量平方成正比，由此导致低流量区的不确定度显著上升，见表8-3。表8-3低流量工况下的适应性分析序号qn占最大流量比Δp/Pa占差压量程比Urel（k=2评价18000100%16000100%1.08%优良2500062.5%625039.1%1.09%优良3250031.3%1562.59.8%1.32%可接受4150018.8%562.53.5%2.39%不满足贸易交接要求由表8-3可见，本站计量系统的有效计量下限约为2500m³/h。针对可能出现的小流量工况，建议采取下列措施之一：（1）减小差压量程。将差压变送器量程由0～16kPa调整为0～4kPa（切换式双量程），在1500m³/h工况下uΔpΔp由2.133%降至0.533%，Urel（2）更换孔板口径。按低流量工况重新设计孔径（增大β使差压升高），但需同步校核压损与流速。（3）改用宽量程仪表。在低流量时段切换至并联的超声流量计或涡轮流量计支路计量。8.4计量方式适应性比选为便于后续改造决策，本节对差压式（孔板）、气体超声、气体涡轮三种计量方式在相同工况条件下的适应性进行比选，见表8-4。表8-4天然气计量方式比选比较项目标准孔板（差压式）气体超声流量计气体涡轮流量计主要标准GB/T21446—2025GB/T18604GB/T21391典型准确度1.0%～1.5%（未校准）0.5%～1.0%0.5%～1.0%量程比3∶1（受差压限制）20∶1以上10∶1～20∶1永久压力损失大（约Δp的0.77倍）极小中等温压补偿要求需补偿（灵敏度0.5）需补偿（灵敏度1.0）需补偿（灵敏度1.0）对直管段要求高（上游20D）高（上游10D～20D）中等对气质要求宽松对脏污、液滴敏感对颗粒敏感初始投资低高中等维护量低（孔板需定期检查）中（探头需清洗标定）中（轴承需维护）适用场合中小流量、压损不敏感大口径、高价值交接中高压力、洁净气比选结论：本站流量规模（≤8000m³/h）、中压工况与投资规模条件，采用标准孔板差压式流量计配实时温压补偿在技术与经济两方面均合理；若未来园区用气量跃升、或对压力损失与量程比提出更高要求，可优先考虑改造为气体超声流量计，并对配套温压补偿仪表同步提高精度等级（压力变送器±0.075%FS、温度变送器±0.2℃）。8.5经济性分析温压补偿系统的经济性可从“投入—避免损失”两个角度评价，见表8-5。表8-5温压补偿经济性分析序号项目计算依据数值1年交接气量（按常用流量折算）5000m³/h×8000h×0.62400万m³2温度波动导致的未补偿偏差全年5～30℃−2.70%～+1.75%3仅温度一项对应的年度计量差异按平均1.0%、气价2.5元/m³约60万元4压力波动导致的未补偿偏差0.30～0.50MPa−8.79%～+11.89%5温压补偿装置投资变送器、流量计算机、安装调试约10万元6孔板永久压力损失能耗费用年运行8000h、负荷率0.6约1.83万元/年7投资回收期第3项相对第5项不足1个季度分析表明，温压补偿装置的投入相对于其避免的计量偏差损失可以忽略不计，具有极高的经济性。同时应注意到，孔板的永久压力损失（1.83万元/年）是差压式计量方式的固有代价，在流量规模扩大后应纳入改造评估。8.6气候与工况适应性（1）温度适应性。本站全年环境与输气温度区间0～35℃，处于温压补偿模型与仪表量程的中部，补偿计算无外推风险；冬季极端低温（−10℃）时，通过引压管保温与电伴热保证压力、差压回路正常工作。（2）压力适应性。工作压力区间0.25～0.60MPa（表压），对应绝压0.35～0.70MPa。压力变送器量程0～1.0MPa时，最低工作点信号利用率为35%，变送器工作在有效量程区间内。（3）组分适应性。在线色谱仪按5min周期更新组分，流量计算机按实测组分实时计算压缩因子与密度，可适应季节调峰期的气质波动（如冬季掺混LNG后C₂以上组分升高）。（4）提压运行适应性。若上游提压至0.60MPa（表压），绝压701.325kPa，工作条件下压缩因子由0.98915降至约0.985，密度增至5.02kg/m³，差压不变时质量流量约增大17.5%，最大流量将超出设计值，需同步核算孔板过载能力与差压变送器量程；由于温压补偿模型对压力变化的适应是连续的，流量计算本身无需修改参数，仅需校核硬件裕量。9安全核验温压补偿设计计算不仅应满足计量准确度要求，还必须通过强度、密封、防爆、检测报警等方面的安全核验，确保计量系统在各种工况下均不构成安全隐患。本章按现行标准逐项核验。9.1管道强度与压力等级核验计量管道的计算壁厚按GB50028—2006（2020年版）的规定计算： δ=pD2σφ式中：δ为管道计算壁厚，mm；p为设计压力，MPa；D为管道外径，mm；σ为管材许用应力，MPa；φ为焊缝系数（无缝钢管取1.0）；C为腐蚀裕量，mm。代入p=0.8MPa、D=219mm、σ=137MPa（20#无缝钢管）、φ=1.0、C=0.5mm δ=0.8×2192×137×1.0+0.5=1.139mm 同时按GB50251—2015输气管道强度公式校核： δ=pD2σs代入管材最低屈服强度σs=245MPa、设计系数F=0.4、温度折减系数t=1.0，得δ=1.394两种方法计算结果均远小于选用管材的实际壁厚7mm，壁厚裕量达5.02～6.15倍，管道强度满足要求。管道、法兰、阀门及仪表按PN1.6（1.6MPa）压力等级配置，为设计压力0.8MPa的2倍，压力等级满足要求。表9-1管道强度与压力等级核验序号核验项目计算值选用值／限值判定1计算壁厚（GB50028—2006）1.139mm实选7mm满足2计算壁厚（GB50251—2015）1.394mm实选7mm满足3工作压力下计算壁厚0.900mm实选7mm满足4管道压力等级0.8MPa（设计）PN1.6满足5壁厚裕量倍数—≥5倍满足9.2节流装置与引压系统安全核验节流装置与引压系统的安全核验项目见表9-2。表9-2节流装置与引压系统安全核验序号核验项目设计要求设计取值判定1孔板厚度与强度按GB/T21446—2025确定，承受Δp不产生永久变形按标准选型满足2孔板夹持器密封金属缠绕垫片，耐压≥PN1.6不锈钢缠绕垫满足3引压管规格与耐压不锈钢管，耐压远高于工作压力φ14×2，316L满足4引压管坡度不小于1∶12，坡向取压口按1∶12施工满足5引压管防冻保温＋电伴热全程保温、伴热满足6排污与排气设排污阀、排气阀各1只／回路满足7差压变送器静压额定静压≥1.5倍工作压力4MPa（工作0.8MPa）满足8差压变送器单向过载≥额定静压4MPa满足9直管段长度上游≥20D、下游≥10D4100mm／2050mm满足9.3仪表防爆与电气安全本站计量橇处于天然气泄漏可能形成的爆炸性气体环境（2区），仪表防爆与电气安全核验见表9-3。表9-3仪表防爆与电气安全核验序号核验项目设计要求设计取值判定1防爆型式与等级不低于ExdⅡBT4Gb或ExiaⅡCT4ExdⅡBT4Gb／ExiaⅡCT4满足2仪表防护等级不低于IP65IP65～IP67满足3保护接地电阻≤4Ω≤4Ω满足4防静电接地电阻≤100Ω≤100Ω满足5信号电缆阻燃、屏蔽、铠装阻燃屏蔽铠装电缆满足6供电可靠性仪表由UPS供电，后备时间≥30minUPS供电满足7电磁兼容满足GB/T17626系列要求通过EMC试验满足9.4可燃气体检测与安全联锁按GB/T50493—2019的要求，在计量橇孔板节流装置、引压管路接头、法兰连接处等易泄漏部位设置可燃气体探测器，探测器安装高度按天然气密度小于空气的特点设于泄漏点上方0.5～1.0m范围内。探测器报警设定与联锁逻辑如下：（1）一级报警：可燃气体浓度达到25%LEL，站控系统声光报警，值班人员现场检查；（2）二级报警：可燃气体浓度达到50%LEL，联锁关闭计量橇进出口切断阀，同时启动强制通风与放散；（3）安全联锁回路按GB/T50770—2013的要求配置，安全完整性等级不低于SIL1，联锁逻辑独立于过程控制系统。9.5水合物与冻堵风险核验天然气经孔板节流后压力降低、温度随之下降，温降按焦耳-汤姆逊效应估算： ΔT=μJTΔp 式中：μJT为焦耳-汤姆逊系数，对本站天然气取0.45K/MPa；Δp为节流压差，MPa代入最大差压Δp=0.016MPa： ΔT=0.45×0.016=0.0072K （9-5）即最大差压点节流温降仅0.007K，可忽略不计。结合本站天然气含水量符合GB17820—2018规定的露点要求，孔板节流部位不存在水合物生成与冻堵风险。引压管路的冻堵风险主要来自环境低温和冷凝液滞留，已通过保温、电伴热与坡度敷设措施予以消除。9.6超压保护与放散系统（1）超压保护。计量橇进口设安全切断阀，设定压力0.60MPa（表压）；出口设安全放散阀，设定压力0.65MPa（表压），放散管引至站场放散总管，放散管高度不低于周围建（构）筑物2m且不小于10m。（2）超压校核。上游来气最高压力0.80MPa（表压），低于管道设计压力0.8MPa与设备压力等级PN1.6，且在安全切断阀动作压力之上，超压保护逻辑成立。9.7安全核验结论综合本节各项核验，本站天然气计量温压补偿系统的安全核验结论见表9-4。全部核验项目均满足现行标准要求，设计方案在安全方面可行。表9-4安全核验结论汇总序号核验类别核验结论1管道强度与压力等级壁厚裕量≥5倍，压力等级PN1.6，满足要求2节流装置与引压系统密封、耐压、防冻、排污措施齐全，满足要求3仪表防爆与电气安全防爆等级ExdⅡBT4Gb，接地电阻合格，满足要求4可燃气体检测与联锁两级报警＋联锁切断，满足GB/T50493—20195水合物与冻堵风险节流温降0.007K，无水合物风险6超压保护与放散切断阀0.60MPa、放散阀0.65MPa，满足要求10结果分析与工程指导意义10.1主要计算结果汇总本计算书的主要计算结果汇总见表10-1。表10-1温压补偿设计计算主要成果汇总序号计算项目符号单位计算结果1天然气摩尔质量Mkg/kmol17.5092相对密度G—0.60453标准参比条件密度ρkg/m³0.72954工作条件密度ρkg/m³3.64085工作条件压缩因子Z—0.989156高位体积发热量HMJ/m³38.9747沃泊指数WMJ/m³50.138孔板开孔直径dmm99.219直径比β—0.4839610流出系数C—0.6027711可膨胀性系数ε—0.9908512最大流量点雷诺数R—9.32×13永久压力损失δpPa1225314温压补偿系数（设计点）K—4.98915扩展相对不确定度（常用点）U%1.09（k=2）16年折算压损能耗费用—元/年1833010.2计算结果综合分析（1）物性参数合理。计算所得摩尔质量17.509kg/kmol、相对密度0.6045、标况密度0.7295kg/m³、高位发热量38.974MJ/m³，与气质组成（甲烷92.5%、C₂以上6.1%）相符，各项参数均落在GB17820—2018二类气与GB/T11062—2020典型管输富气的取值区间内。（2）节流装置设计合理。直径比0.48396位于0.4～0.6的推荐区间，最大流量点雷诺数9.32×105远高于5000的下限，永久压力损失占工作压力约3.05%，各项几何与流动参数全部落在GB/T21446—2025与（3）温压补偿精度满足要求。补偿后体积流量测量系统扩展相对不确定度在常用流量点为1.09%（k=2），在最大流量点为1.08%，满足B级计量站的准确度要求；补偿环节自身引入的残余偏差仅0.122%，不构成不确定度的主导因素。（4）温压补偿的必要性得到量化证明。不实施温压补偿时，仅温度由20℃升高至30℃即可造成+1.75%（差压式）的计量偏差，压力由0.40MPa降至0.30MPa时偏差达+11.89%，远超贸易交接允许的误差范围。（5）存在明确的优化空间。流出系数不确定度（0.50%）占总方差的84.1%，是系统准确度的制约因素；实施在线实流校准可使扩展不确定度由1.09%降至0.74%，达到A级交接站水平。10.3对工程设计的指导意义（1）口径与量程的匹配性设计。本设计表明：孔板直径比的选择是“准确度—压力损失—量程比”三者的折中。当设计流量范围较宽时，应优先保证常用流量点处于差压量程的30%～70%区间，最小流量点不低于10%。本站最小流量点差压占量程9.8%，已接近下限，说明量程比不宜进一步扩大。（2）差压量程的确定原则。差压上限应结合孔板孔径迭代确定，不应先定孔径再配变送器。本设计中若将差压上限取为25kPa，则迭代所得孔径为76.2mm、直径比约0.372，虽然压损略有下降，但远离推荐区间、孔径加工与直管段要求更严，说明16kPa的取值是合理的。（3）补偿环节的优先级排序。由不确定度评定可知，提高系统准确度应优先改善流出系数的可信度（实流校准），其次才是提高变送器精度等级。将差压变送器精度由0.075%FS提高到0.05%FS，扩展不确定度仅由1.09%降至约1.07%，投入产出比很低。（4）为后续计量方式升级预置条件。本设计在管道、直管段、通信与电源等方面为后续改造为气体超声流量计预留了条件，包括直管段长度（上游20D满足超声流量计要求）、仪表供电与通信接口、以及流量计算机的多通道输入能力。10.4对采购与施工的指导意义（1）采购技术要求。孔板开孔直径应按d20=99.21mm加工，公差控制在±0.05mm以内，孔板入口直角边必须锐利无损；差压变送器须为高静压型，额定静压不低于4MPa；压力变送器量程0～1.0MPa、精度±0.1%FS；温度变送器采用A级Pt100，精度（2）安装施工要点。孔板安装方向必须与流向一致（锐边朝上游）；上游20D、下游10D范围内不得设置任何阻力件与取压口；引压管两侧对称敷设、长度相等，坡度1∶12坡向取压口；引压管与孔板夹持器连接处不得有泄漏，安装完成后应进行严密性试验。（3）仪表接线与接地。信号电缆屏蔽层应在流量计算机侧单端接地，避免形成地环流；防爆区域内的接线应符合GB3836.1的要求，隔爆面不得损伤，本安回路应使用蓝色标识的专用电缆。10.5对运行维护与计量管理的指导意义（1）运行监视重点。日常运行应重点监视三项参数：孔板上游压力是否处于0.25～0.60MPa（表压）设计区间；管内流速是否超过20m/s（最大流量点设计值13.49m/s）；差压占量程比是