天然气贸易计量能量计量设计计算书

×××天然气长输管道工程天然气贸易计量能量计量设计计算书——计量系统设计与能量计量计算——文件编号：TRQ-JS-2026-018设计阶段：施工图设计文件类型：设计计算书编制：______________________审核：______________________批准：______________________×××石油天然气工程设计有限公司二〇二六年十一月

目录TOC\z\o"1-3"\u\h1总则 11.1项目概况 11.2设计依据 11.3设计范围与计量界面 21.4主要设计参数 21.5符号与单位 32天然气基础物性参数计算 42.1计算方法与原理 42.2天然气组成数据 42.3摩尔质量与相对密度 52.4压缩因子 62.4.1参比条件压缩因子 62.4.2工况条件压缩因子 62.5高位发热量 72.6低位发热量 82.7密度与沃泊指数 82.8物性参数计算结果汇总 93计量系统设计与流量计算 103.1计量系统组成与等级 103.2流量计选型与计量管路 103.3标准参比条件流量换算 113.4超声流量计测量原理与流量计算 123.5流动适应性校核 124能量计量计算 144.1能量计量原理与数学模型 144.2发热量测定方式 144.3瞬时能量流量计算 154.4累积能量与贸易结算计算 154.5数据采集与处理 155测量不确定度评定 175.1评定依据与数学模型 175.2标准不确定度分量 175.3合成与扩展不确定度 185.4与标准限值比对 196安全核验 206.1安全泄放与超压保护 206.2管道强度与壁厚校核 216.3气质安全指标核验 226.4可燃气体检测报警与联锁 226.5防雷防静电与防爆 236.6计量与检测设备运行安全 237结果分析与工程指导意义 247.1主要计算结果汇总 247.2敏感性分析 247.3计量偏差与贸易风险分析 257.4工程指导意义 258结论与建议 27附录A符号与计量单位 28免责声明 29天然气贸易计量能量计量设计计算书PAGE271总则1.1项目概况本工程为×××天然气长输管道工程配套建设的贸易交接计量站（以下简称“本计量站”），站址位于陕西省延安市境内，承担上游气田来气经净化处理后，向下游城市燃气及工业用户进行分输交接计量的任务。按照我国天然气贸易由体积计量向能量计量转型的总体要求，本计量站按能量计量方式设计，计量系统同时具备标准参比条件下体积流量计量与单位体积高位发热量测定的能力，并以二者之积作为贸易交接的能量结算依据。本计算书依据现行国家标准与计量技术规范，对天然气组成物性参数、标准参比条件体积流量、发热量、能量计量结果、测量不确定度以及站场安全核验等内容进行系统计算，为计量系统设备选型、仪表精度配置、控制系统组态及贸易结算方式确定提供设计依据。1.2设计依据本计算书采用的主要技术标准、规范及计量技术文件见表1-1。标准均采用现行有效版本，凡引用文件中的条款通过本计算书的引用而构成为本计算书的条款。表1-1主要设计依据标准/规范编号标准名称本计算书适用内容GB/T18603—2023天然气计量系统技术要求系统分级、准确度、仪表配置GB/T22723—2024天然气能量的测定能量测定原理与数学模型GB/T11062—2020天然气发热量、密度、相对密度和沃泊指数的计算方法发热量、密度、沃泊指数计算GB/T17747.2—2011天然气压缩因子的计算第2部分：用摩尔组成进行计算压缩因子计算（AGA8—92DC法）GB/T13610—2020天然气的组成分析气相色谱法组分分析GB/T13609—2017天然气取样导则在线取样GB/T18604—2014用气体超声流量计测量天然气流量超声流量计流量计算GB17820—2018天然气商品天然气质量指标GB/T37124—2018进入天然气长输管道的气体质量要求管输气质要求GB50251—2015输气管道工程设计规范管道强度与站场设计GB150.1—2011压力容器第1部分：通用要求安全阀泄放能力计算GB3836.1—2021爆炸性环境第1部分：设备通用要求防爆电气设备选型JJF1993—2022天然气能量计量技术规范能量计量系统评定与验收JJF1059.1—2012测量不确定度评定与表示测量不确定度评定1.3设计范围与计量界面本计算书的设计范围为计量站内的贸易交接计量系统，主要包括：计量管路（含流量计及前后直管段）、压力与温度测量回路、在线气相色谱分析系统、流量计算机与能量计算单元、在线校准接口以及相应的取样、放空与安全保护设施。计量界面设定于计量站分输汇管下游的站界出站法兰处。交接界面上，上游方以交接点的标准参比条件下体积量及对应的高位发热量为结算基准，下游方以本站计量系统示值为接收基准；界面处的压力、温度测点由双方共同确认并加封缄。计量系统的运行参数、组分数据与能量累计值同步上传至站控系统及贸易结算数据平台。1.4主要设计参数本计量站的主要工艺设计参数见表1-2。标准参比条件按GB/T11062—2020的规定取20℃、101.325kPa（干基），燃烧参比条件与计量参比条件一致。表1-2主要设计参数序号参数名称符号/单位设计取值1天然气组成（摩尔分数）—见表2-12设计通过能力（标准参比条件）q_n/(m³·h⁻¹)120000（额定），144000（最大）3管道设计压力p/MPa10.04交接操作压力（绝对）p_f/MPa6.305交接操作温度t/℃10（设计计算取值）6标准参比压力p_n/kPa101.3257标准参比温度T_n/K293.15（20℃）8最低环境温度t_min/℃−209计量方式—超声流量计体积计量＋在线色谱发热量测定10能量计量系统等级—A级（GB/T18603—2023）11设计使用年限a201.5符号与单位本计算书采用的符号、名称、单位及取值约定遵循GB3100、GB3101及GB3102的规定，各符号与计量单位说明详见附录A。为避免歧义，特作如下约定：体积流量除特别注明外均指标准参比条件下的体积流量q_n；发热量除特别注明外均指计量参比条件下的单位体积高位发热量H_s；“高位”与“低位”分别对应燃烧产物中水为液态与气态的发热量。2天然气基础物性参数计算2.1计算方法与原理天然气作为多元混合气体，其发热量、密度、相对密度与沃泊指数等物性参数，按GB/T11062—2020采用已知气体摩尔组成的加权计算方法确定。该方法的基本思路为：先按理想气体模型对纯组分的基础物性参数进行摩尔分数加权求和，再引入混合气压缩因子对真实气体行为进行修正，最终得到真实气体在计量参比条件下的单位体积物性参数。真实气体偏离理想气体的程度由压缩因子表征。对于1mol气体，实际气体状态方程为： pVm=ZRT 式中：p为绝对压力，Pa；pVm=ZRT为摩尔体积，m³·mol⁻¹；Z为压缩因子；R为摩尔气体常数，取8.314462618J·mol⁻¹·K⁻¹；T理想气体状态下Z＝1，此时摩尔体积仅由温度与压力决定。计量参比条件（20℃、101.325kPa）下的理想气体摩尔体积按式（2）计算： Vm0=R代入数值可得Vm0=8.314462618×2.2天然气组成数据本计量站天然气来源为上游气田净化外输气，气源稳定，重烃含量低，属于典型的管输贫气。设计采用的气相色谱全组分分析数据见表2-1，各组分摩尔分数之和已归一化为100%。组分分析方法按GB/T13610—2020执行，取样按GB/T13609—2017执行。表2-1天然气组成分析数据（摩尔分数，干基）序号组分名称分子式摩尔分数x_i/%分析方法1甲烷CH₄93.20气相色谱法2乙烷C₂H₆2.60气相色谱法3丙烷C₃H₈0.55气相色谱法4异丁烷i-C₄H₁₀0.07气相色谱法5正丁烷n-C₄H₁₀0.08气相色谱法6异戊烷i-C₅H₁₂0.02气相色谱法7正戊烷n-C₅H₁₂0.02气相色谱法8正己烷（代表C₆⁺）n-C₆H₁₄0.02气相色谱法9氮气N₂2.44气相色谱法10二氧化碳CO₂1.00气相色谱法—合计—100.00—2.3摩尔质量与相对密度天然气的平均摩尔质量按各组分摩尔分数与纯组分摩尔质量的加权求和计算： M=i=1Nxi式中：M为混合气平均摩尔质量，kg·kmol⁻¹；M=i=1NxiMi为组分i的摩尔分数；Hm0=表2-2组分基础物性参数及加权计算（20℃、101.325kPa）组分M_i/

(kg·kmol⁻¹)H°_m,i(高)/

(kJ·mol⁻¹)H°_m,i(低)/

(kJ·mol⁻¹)√b_ix_i·M_i/

(kg·kmol⁻¹)CH₄16.043891.09802.650.2088114.9521C₂H₆30.0701561.411428.740.299000.7818C₃H₈44.0972220.132043.230.358890.2425i-C₄H₁₀58.1232869.382648.260.412680.0407n-C₄H₁₀58.1232878.572657.450.422230.0465i-C₅H₁₂72.1503530.243264.890.465620.0144n-C₅H₁₂72.1503537.173271.830.484250.0144n-C₆H₁₄86.1774196.583888.600.533480.0172N₂28.0135000.131530.6835CO₂44.010000.432560.4401合计————17.2333注：表中M_i、H°_m,i、√b_i等基础数据引自GB/T11062—2020；n-C₆H₁₄的低位摩尔发热量按烷烃同系物高低位差值规律外推确定。由表2-2可得混合气平均摩尔质量M＝17.2333kg·kmol⁻¹。理想气体相对密度按混合气摩尔质量与干燥空气摩尔质量之比计算，空气摩尔质量按GB/T11062—2020取M_air＝28.9626kg·kmol⁻¹： d0=MM计算得d0 d=d0Zair取空气在参比条件下的压缩因子Z_air＝0.99958，代入得d=0.5950×0.999582.4压缩因子2.4.1参比条件压缩因子计量参比条件下（20℃、101.325kPa）的混合气压缩因子按GB/T11062—2020规定的求和因子法计算： Zmix=1-i=1N式中：b_i为组分i的求和因子，其平方根见表2-2；Z_mix为混合气在参比条件下的压缩因子。逐项计算得∑xibZ该值为参比条件下由理想气体体积发热量换算为真实气体体积发热量的修正系数。2.4.2工况条件压缩因子在工作条件下（p_f＝6.30MPa、t＝10℃，即T_f＝283.15K），压缩因子按GB/T17747.2—2011（AGA8—92DC法）计算。为便于设计阶段复核，本节同时给出基于对应态原理（Pitzer—Curl三参数关联式）的近似计算，两者结果偏差应在±0.5%以内，工程应用中以AGA8—92DC法计算结果为准。对应态法采用混合气假临界参数与偏心因子，先计算对比参数： Tr=TT pr=pp本天然气组成按摩尔分数加权得假临界温度T_c＝194.61K、假临界压力p_c＝4.600MPa、偏心因子ω＝0.0178。代入工作条件得T_r＝1.4550、p_r＝1.3694。第二维里系数的Pitzer关联式为： B0=0.083-0.422T B1=0.139-0.172T B=RTcp计算得B⁰＝−0.14858、B¹＝0.10340，第二维里系数B＝−5.162×10⁻⁵m³·mol⁻¹。压缩因子按截断式维里方程计算： Z=1+BpRT （代入数值可得Zf=1+-5.162×10-5×6.30×1068.3145×283.15=0.8620。采用AGA8—92DC法复核，结果为2.5高位发热量天然气的高位发热量（即总发热量）按各组分高位发热量的摩尔分数加权求和计算。理想气体摩尔高位发热量为： Hm0=i=1逐项加权求和得H_m⁰＝889.03kJ·mol⁻¹。理想气体体积高位发热量按理想气体摩尔体积换算： Hv0=i=1计算得Hv Hs=Hv代入得Hs=36.9580.95471=38.712MJ⋅m-3。2.6低位发热量低位发热量按燃烧产物中水为气态时的发热量计算，用于能量利用效率分析及下游用户用能核算。理想气体摩尔低位发热量为： Hl,m0=i=1逐项加权求和得H_l,m⁰＝801.74kJ·mol⁻¹。理想气体体积低位发热量为801.74×10 Hl=Hl,v计算得Hl=33.3300.95471=34.911MJ⋅m-3。高位与低位发热量之差为3.801MJ·m⁻³，相当于每立方米天然气燃烧生成水蒸气的凝结潜热，该差值与按燃烧反应计量系数与水的汽化潜热（25℃2.7密度与沃泊指数理想气体密度按平均摩尔质量与理想气体摩尔体积计算： ρ0=MV计算得ρ0 ρ=MpZRT （代入数值可得ρ=17.2333×10沃泊指数是反映燃气互换性的重要指标，按高位发热量与相对密度平方根之比计算： Ws=Hs计算得Ws=38.7120.6230=49.05MJ⋅2.8物性参数计算结果汇总天然气在计量参比条件（20℃、101.325kPa）下的主要物性参数计算结果汇总于表2-3。表2-3天然气主要物性参数计算结果物性参数符号单位计算结果平均摩尔质量Mkg·kmol⁻¹17.2333理想气体摩尔体积V_m⁰m³·kmol⁻¹24.0545参比条件压缩因子Z_mix—0.95471工况条件压缩因子（6.30MPa、10℃）Z_f—0.8620真实气体相对密度d—0.6230真实气体密度ρkg·m⁻³0.7504高位发热量（摩尔基）H_m⁰kJ·mol⁻¹889.03高位发热量（体积基）H_sMJ·m⁻³38.712低位发热量（体积基）H_lMJ·m⁻³34.911高位沃泊指数W_sMJ·m⁻³49.053计量系统设计与流量计算3.1计量系统组成与等级依据GB/T18603—2023，按设计通过能力将天然气计量站贸易计量系统划分为A、B、C三级。本计量站设计通过能力为120000m³·h⁻¹（标准参比条件），属于大型贸易交接计量站，按A级系统设计。A级系统应采用流量计、温度与压力测量装置及在线气相色谱分析仪的完整配置，具备实时组分分析与发热量在线测定的能力，并对计量设备准确度、可靠性与校准提出更高要求。计量系统由计量管路与流量计、压力测量回路、温度测量回路、在线气相色谱分析系统、流量计算机与能量计算单元、在线校准接口及数据采集与传输装置等组成，主要构成见表3-1。表3-1能量计量系统主要组成序号子系统/设备主要技术要求数量1气体超声流量计DN250，准确度等级0.5级，双向计量1台2压力变送器0.1级，量程0～10MPa2台（冗余）3铂电阻温度计A级，Pt100，允许误差±0.5K2支（冗余）4在线气相色谱分析仪分析C₁～C₆⁺、N₂、CO₂等组分1套5流量计算机/能量计算单元支持实时压缩因子与发热量计算1套6在线实流校准接口DN250，具备双向校准条件1套7数据采集与传输装置与站控系统及结算平台通信1套3.2流量计选型与计量管路天然气贸易交接计量流量计的选型应综合考虑量程比、压力损失、准确度、直管段要求及在线校准条件。气体超声流量计具有量程比宽（可达1∶50以上）、无阻碍测量、压力损失小、可双向计量且易于实流校准等优点，为A级计量系统的首选流量计类型，故本计量站选用气体超声流量计，口径DN250，准确度等级0.5级。计量管路内径按流量计上游直管段内径确定，取管道规格Φ273×12.7，内径D_i＝273.1−2×12.7＝247.7mm，流通面积按式（21）计算： A=πDi2计算得A=π×0.247724=0.04819m3.3标准参比条件流量换算气体超声流量计直接测量的是工作条件下的实际体积流量q_f，而贸易交接以标准参比条件下的体积流量q_n为基准，二者按气体状态方程换算： qn=qf式中：p_f、T_f分别为工作条件下的绝对压力与热力学温度；p_n、T_n分别为标准参比压力与温度；Z_f、Z_n分别为工作条件与参比条件下的压缩因子。按表1-2及第2章计算结果代入各参数，额定流量工况下：q_n＝120000m³·h⁻¹，pfpn=由此可得额定工况下的实际体积流量qf=12000062.176×1.0353×1.1076=1683.1m3⋅表3-2流量换算计算结果计算项目符号单位额定工况最大工况标准参比条件体积流量q_nm³·h⁻¹120000144000工作条件体积流量q_fm³·h⁻¹1683.12019.7工作条件体积流量q_fm³·s⁻¹0.46750.5610管道内径D_im0.24770.2477流通面积Am²0.048190.04819工作条件轴向平均流速vm·s⁻¹9.7011.64工作条件密度ρ_fkg·m⁻³53.5053.50质量流量q_mkg·h⁻¹900461080553.4超声流量计测量原理与流量计算气体超声流量计采用传播时间差法测量气体流速。流量计在管道上布置若干对超声换能器，声波在顺流与逆流方向的传播时间之差与气体轴向流速成正比。单一声道的轴向流速按下式计算： vj=Lj式中：L_j为第j声道换能器之间的声程长度，m；θ_j为声道与管道轴线的夹角，（°）；t_1j、t_2j分别为顺流与逆流传播时间，s。流量计的输出流量为各声道流速按权重系数的加权组合： qf=Aj=1n式中w_j为第j声道权重系数，由流量计实流校准确定，Σw_j＝1。由于超声流量计直接输出工作条件下体积流量，实际应用时需结合压力、温度及压缩因子按式（22）换算为参比条件下的体积流量，换算在流量计算机内在线完成。3.5流动适应性校核为保证超声流量计在有效范围内工作，需对工作条件下雷诺数进行校核。雷诺数按式（25）计算： Re=ρfvD取工作条件下天然气动力黏度μ＝1.2×10⁻⁵Pa·s，额定工况下Re=53.50×9.70×0.24771.2×10-5=1.07×10同时校核流量计的量程比适用范围：本计量站最小流量一般不低于额定流量的10%，对应q_n,min＝12000m³·h⁻¹，换算后工作条件流速约0.97m·s⁻¹，仍高于超声流量计可靠测量的下限流速，量程比满足1∶10以上的工艺要求。管道压力损失按沿程阻力公式校核： Δp=λLDiρ λ=0.11ΔDi0.25取管道粗糙度Δ＝0.05mm，计算得λ＝0.0129；计量管段长度约15m，额定工况下沿程压力损失约0.4kPa，占操作压力的0.006%，可忽略不计，表明超声流量计无阻碍测量带来的压力损失极小，节能效益显著。4能量计量计算4.1能量计量原理与数学模型天然气能量计量是通过对气体体积（或质量）流量与单位体积（或质量）发热量两个独立参数的测量，进而合成得到能量值的计量方式。按GB/T22723—2024《天然气能量的测定》的规定，能量测定采用以时间变化为基础的间接测定法，即在某一计量周期内，能量值为气体流量与对应发热量的乘积对时间的积分： E=t1t2式中：E为计量周期内的天然气能量，MJ；q_n(t)为标准参比条件下的体积流量，m³·h⁻¹；H_s(t)为单位体积高位发热量，MJ·m⁻³。当计量周期内流量与发热量变化平缓时，可用周期起止时刻的平均值简化计算；对于本计量站这类连续平稳输气的贸易交接场合，采用分段累积的离散形式更为实用： E=k=1KQn,k式中：Q_n,k为第k时段的标准参比条件体积量，m³；H_s,k为第k时段的高位发热量，MJ·m⁻³；K为计量周期内的时段总数。该式为能量计量的基本计算模型，也是流量计算机与结算平台的核心算法。4.2发热量测定方式单位体积高位发热量的获取方式依据GB/T18603—2023的规定分为在线测定、离线测定与赋值三类。本计量站为A级系统，采用在线测定方式：由在线气相色谱分析仪连续分析天然气全组分，经归一化处理后，按GB/T11062—2020的规定实时计算单位体积高位发热量H_s，并将结果同步送至流量计算机参与能量计算。在线气相色谱分析仪的采样周期不大于5min，测量结果经滑动平均后作为当前时段的发热量计算值，以抑制组分分析随机波动的影响。发热量测定的相对扩展不确定度应不大于0.5%（k＝2），与本计算书第5章不确定度评定结果（0.50%）一致。当在线色谱分析仪故障或检修时，按GB/T22723—2024的规定采用经确认的代表性发热量进行赋值计算，并在结算中注明赋值方式及有效期。4.3瞬时能量流量计算瞬时能量流量为单位时间内的能量值，按标准参比条件下的体积流量与单位体积高位发热量之积计算： e=qnHs额定工况下（q_n＝120000m³·h⁻¹，H_s＝38.712MJ·m⁻³）：e=120000×38.712=4645440MJ⋅h-1，即464.5GJ·h⁻¹。按当日连续运行24h计，日输送能量为11149GJ；最大工况下瞬时能量流量可达4.4累积能量与贸易结算计算贸易交接通常以月为结算周期。计量周期内的累积能量为各时段体积量与对应发热量乘积的累加值。当采用周期平均发热量时，平均高位发热量按下式确定： Hs=k=1K累积能量则可由周期累积体积量与平均发热量之积直接计算： E=QnHs以某月（30d）连续运行结算为例：该月标准参比条件累积体积量Q_n＝8.50×10⁷m³，加权平均高位发热量H̄_s＝38.65MJ·m⁻³，则该月交接能量为：E=8.50×按设计通过能力测算，本计量站年标准参比条件输气量约9.6×10⁸m³（按年运行8000h计），对应年交接能量约3.72×10⁷GJ。与传统体积计量方式相比，能量计量在同等体积量下能够反映天然气实际热值的差异，结算结果更公平、更能反映天然气的使用价值。4.5数据采集与处理能量计量系统采用流量计算机作为核心计算单元，实时接收超声流量计、压力变送器、温度变送器及在线气相色谱分析仪的测量信号，按式（22）完成标准参比条件体积流量换算，按式（30）完成瞬时能量流量计算，并按式（29）进行累积。所有原始测量数据与计算结果同步上传站控系统与贸易结算平台，存储周期不小于10a。为满足量值溯源与数据可追溯要求，压力、温度、流量与组分分析通道均应具备实时数据质量标识（有效性标志），当任一测量通道数据无效时，自动切换至冗余通道或按规定的替代值策略处理，并在数据记录中标注替代原因与持续时间。5测量不确定度评定5.1评定依据与数学模型能量计量系统的测量不确定度按JJF1059.1—2012《测量不确定度评定与表示》和GB/T18603—2023的规定评定，量值不确定度评定参照GB/T35186—2017《天然气计量系统》。能量值由体积流量与发热量两个相对独立的参数合成，其相对合成标准不确定度按平方和开方的方式计算： uEE2=式中：u_E/E为能量测量的相对合成标准不确定度；u_Q/Q为体积流量测量的相对标准不确定度；u_H/H为高位发热量测量的相对标准不确定度。其中体积流量测量的相对标准不确定度又由流量计、压力、温度与压缩因子等分量合成： uQQ2=式中各分量的相对标准不确定度由对应仪表的允许误差或校准结果按均匀分布折算。对于以引用误差表示的仪表，其相对标准不确定度按式（35）计算： upp=Δ uTT=Δ式中：Δ_P为压力测量的最大允许误差，MPa；Δ_T为温度测量的绝对误差，K；分母中的√3为均匀分布的包含因子。5.2标准不确定度分量（1）流量计分量。气体超声流量计准确度等级0.5级，经在线实流校准并考虑安装条件影响，取相对标准不确定度u_q/q＝0.20%。（2）压力分量。压力变送器准确度0.1级，量程上限10MPa，工作压力6.30MPa，最大允许误差ΔP=0.1%（3）温度分量。铂电阻温度计允许误差±0.5K，工作温度T_f＝283.15K，故uT（4）压缩因子分量。按GB/T17747.2—2011的AGA8—92DC法计算，其预期不确定度约0.1%，取u_Z/Z＝0.10%。（5）发热量分量。在线气相色谱组分分析引入的相对标准不确定度约0.22%，GB/T11062—2020基础数据的相对标准不确定度约0.10%，按发热量的组分敏感度合成： uHH=1计算得uHH=表5-1标准不确定度分量汇总序号不确定度来源评定依据相对标准不确定度/%1流量计及安装影响0.5级，实流校准0.202压力测量0.1级，量程10MPa0.093温度测量±0.5K0.104压缩因子AGA8—92DC法0.105高位发热量在线色谱＋GB/T110620.255.3合成与扩展不确定度将表5-1中流量计、压力、温度与压缩因子四个分量合成为体积流量的相对标准不确定度：u进一步将体积流量与发热量分量合成为能量测量的相对合成标准不确定度：u合成标准不确定度按式（40）取包含因子k＝2，得到能量计量的相对扩展不确定度： U=kuc （ Urel=kuc,rel计算得UE,rel5.4与标准限值比对依据GB/T18603—2023对A级计量系统的准确度要求，能量计量系统能量测量结果的相对扩展不确定度应不大于1.0%（k＝2）。本计量站能量计量系统的相对扩展不确定度为0.72%，优于标准限值，满足贸易交接计量的准确度要求，具备按能量方式计量结算的技术条件。作为对比，若发热量采用赋值法（赋值不确定度按0.5%计），则uEE=0.262+0.502%6安全核验6.1安全泄放与超压保护计量站内工艺管道与设备应设置可靠的压力泄放设施，防止因下游阀门误关闭、调压系统故障或气源超压导致管道超压。安全泄放量按最不利工况确定：当下游截断阀关闭时，上游最大来气量需全部泄放。按最大工况标准参比条件体积流量144000m³·h⁻¹、参比条件密度0.7504kg·m⁻³计算，所需泄放质量流量为：W安全阀泄放能力按临界流动条件计算（依据GB150.1—2011附录B及API520临界流方程）： W=KdAp式中：W为安全阀额定泄放量，kg·h⁻¹；K_d为泄放系数，全启式安全阀取0.975；A为最小泄放面积，m²；p₁为泄放压力（绝对），Pa；κ为气体绝热指数，取1.30；M为摩尔质量，kg·kmol⁻¹；Z为泄放工况下的压缩因子；T₁为泄放温度，K；R为摩尔气体常数。安全阀整定压力取1.1倍操作压力，即p_s＝1.1×6.30＝6.93MPa；泄放压力取1.1倍整定压力，即p₁＝1.1×6.93＝7.62MPa；泄放温度取操作温度283.15K；泄放工况压缩因子按2.4.2节方法计算得Z＝0.833。计算得临界流系数项RT=κMZRT1⋅2选用DN100全启式安全阀，喉径d₀＝65mm，最小泄放面积A＝π×65²/4＝3318mm²＝3.318×10⁻³m²，则W=5.291×107×3.318×10-3=175560kg⋅h-1表6-1安全阀泄放能力核算计算参数符号取值/结果最大工况标准体积流量q_n,max144000m³·h⁻¹所需泄放质量流量W_req108055kg·h⁻¹整定压力p_s6.93MPa泄放压力（绝对）p₁7.62MPa泄放温度T₁283.15K绝热指数κ1.30泄放工况压缩因子Z0.833安全阀公称通径/喉径DN/d₀DN100/65mm最小泄放面积A3318mm²安全阀额定泄放量W175560kg·h⁻¹裕量系数W/W_req1.636.2管道强度与壁厚校核计量站内工艺管道按GB50251—2015的规定进行直管段壁厚计算： δ=pD2σs式中：δ为管道计算壁厚，mm；p为设计压力，MPa；D为管道外径，mm；σ_s为管材最小屈服强度，MPa；F为强度设计系数；φ为焊缝系数；t为温度折减系数；C为腐蚀裕量，mm。计量站内管道采用L360无缝钢管，σ_s＝360MPa，外径D＝273.1mm，设计压力p＝10.0MPa；按站场管道取强度设计系数F＝0.4，无缝钢管焊缝系数φ＝1.0，操作温度低于120℃时温度折减系数t＝1.0，腐蚀裕量C＝1.0mm。代入式（43）：δ=10.0×选用标准壁厚12.7mm（Φ273×12.7），实际壁厚大于计算壁厚，裕量充足。对选用壁厚进行应力校核： σ=pD-2δ2δ 计算得σ=10.0×273.1-2×12.72×12.7=97.52MPa；许用应力σ=σs表6-2管道壁厚与强度校核计算参数符号取值/结果设计压力p10.0MPa管道外径D273.1mm管材最小屈服强度σ_s360MPa（L360）强度设计系数F0.4计算壁厚δ_calc10.48mm选用壁厚δ12.7mm实际环向应力σ97.52MPa许用应力[σ]144.0MPa安全裕量系数[σ]/σ1.486.3气质安全指标核验本计量站交接的天然气气质应符合GB17820—2018《天然气》一类气的质量要求，并满足GB/T37124—2018《进入天然气长输管道的气体质量要求》。按上游气田提供的设计气质数据，对关键指标进行核验，结果见表6-3。水分与烃类露点的控制应保证在交接压力下低于最低环境温度5℃以上，以防止在计量管路中形成水合物或凝析液，影响流量测量准确度。表6-3天然气气质指标核验核验项目标准要求设计气质结论高位发热量≥34.0MJ·m⁻³38.712MJ·m⁻³合格总硫（以硫计）≤20mg·m⁻³≤6mg·m⁻³合格硫化氢≤6mg·m⁻³≤2mg·m⁻³合格二氧化碳（摩尔分数）≤3.0%1.00%合格水露点（6.3MPa）低于环境温度5℃≤−10℃合格烃露点不高于最低环境温度≤−5℃合格6.4可燃气体检测报警与联锁计量站属爆炸危险环境，应设置固定式可燃气体检测报警系统。在流量计区、色谱分析小屋、阀门组及法兰连接等易泄漏部位设置可燃气体探测器，探测器采用催化燃烧式或红外式原理，响应时间不大于30s。报警设定值为：一级报警25%LEL，二级报警50%LEL；甲烷爆炸下限LEL取4.4%（体积分数），对应一级报警浓度约1.1%。二级报警应联动站场紧急切断系统，自动关闭进站与出站截断阀，切断气源并启动放空。检测报警系统应独立于过程控制系统，具备声光报警、故障自诊断及数据记录功能，并按GB/T17626系列标准满足电磁兼容要求。6.5防雷防静电与防爆计量站的防爆电气设备选型应符合GB3836.1—2021及相应各部分的防爆要求，爆炸危险区域内的仪表、接线箱、照明及电气设备应采用隔爆型（Exd）或本质安全型（Exia），防爆等级不低于ExdⅡBT4。站内管道、设备及仪表外壳均应可靠接地，防静电接地电阻不大于100Ω，防雷接地电阻不大于10Ω，并采取等电位联结措施。超声流量计换能器接线、色谱分析仪载气管路及仪表信号电缆的敷设应符合防爆与电磁兼容要求，信号电缆采用屏蔽电缆并单端接地。6.6计量与检测设备运行安全在线气相色谱分析仪及其分析小屋应设置可燃气体检测、强制通风与超温联锁保护；标准气瓶应固定并远离热源，载气管路应设置超压泄放。超声流量计换能器应能承受管道最高工作压力并具备在线更换条件。流量计算机、压力与温度变送器的供电与信号回路应设置过压保护与浪涌保护装置，保证在雷暴天气下的运行安全。计量系统整体应按GB/T18603—2023的要求进行投产试运与验收，并在运行中执行周期检定/校准与实流比对，确保量值溯源与计量准确度长期稳定。7结果分析与工程指导意义7.1主要计算结果汇总本计算书的主要计算结果汇总于表7-1。由计算结果可见，各物性参数之间相互印证、结果自洽，计量系统配置与不确定度评定结果均满足现行标准要求。表7-1主要计算结果汇总项目符号结果备注天然气平均摩尔质量M17.2333kg·kmol⁻¹GB/T11062参比条件压缩因子Z_mix0.95471求和因子法工况条件压缩因子Z_f0.8620AGA8—92DC高位发热量H_s38.712MJ·m⁻³20℃、101.325kPa低位发热量H_l34.911MJ·m⁻³20℃、101.325kPa真实相对密度d0.6230—高位沃泊指数W_s49.05MJ·m⁻³互换性指标额定瞬时能量流量e464.5GJ·h⁻¹q_n＝120000m³·h⁻¹年交接能量E_a3.72×10⁷GJ按8000h/a能量计量相对扩展不确定度U_E,rel0.72%（k＝2）优于1.0%限值7.2敏感性分析为识别影响能量计量准确度的关键因素，对各输入参量进行敏感性分析。敏感度系数定义为输出量相对变化与输入量相对变化之比： S=ΔEEΔ依据式（22）、式（30）的数学关系，标准参比条件体积流量与能量值对压力、温度、压缩因子与发热量的敏感度系数见表7-2。需要特别说明的是：标准参比条件体积流量与工作条件压力近似成正比（q_n∝p_f），与工作条件热力学温度近似成反比（q_n∝1/T_f），与工作条件压缩因子成反比（q_n∝1/Z_f），与发热量成正比。表7-2各输入参量对能量计量的敏感度输入参量相对变化能量值相对变化主要影响工作压力p_f＋1%约＋1.00%压力测量与取压工作温度T_f＋1%约−1.00%温度测量与安装压缩因子Z_f＋1%约−1.00%组分分析与状态方程高位发热量H_s＋1%约＋1.00%色谱分析与组分体积流量q_f＋1%约＋1.00%流量计与校准敏感性分析表明，压力、温度、压缩因子与发热量对能量计量结果均具有接近1∶1的传递关系，任一参量1%的偏差均会导致约1%的能量结算偏差。因此，计量系统必须同步保证压力、温度测量的准确度与稳定性，并保持压缩因子与发热量计算的正确性；其中温度测量采用冗余配置、压力测点靠近流量计、色谱分析仪定期校准等措施是控制计量偏差的关键。7.3计量偏差与贸易风险分析按本计算书的年交接能量3.72×10⁷GJ、以单位体积气价2.85元·m⁻³（对应标准参比条件体积量9.6×10⁸m³）估算，本计量站年贸易额约27.4亿元。在能量计量系统相对扩展不确定度0.72%的条件下，理论上可能出现的最大计量偏差对应的经济量约为：ΔE=0.72%×27.4=0.197（亿元），约合这一量级说明，计量系统的准确度直接关系到贸易双方的重大经济利益。在体积计量方式下，由于未考虑发热量差异，当上游气源热值波动或下游掺混不同气源时，折算偏差可能远大于计量系统本身的不确定度，进而引发结算争议。采用能量计量方式后，发热量参与结算，不同品质天然气的价值得以公平体现，可从机制上减少计量争议。7.4工程指导意义（1）关于计量方式。本计算结果表明，对于本计量站这类设计通过能力大、气源组分相对稳定的A级贸易交接计量站，采用“超声流量计体积计量＋在线色谱发热量测定”的能量计量方案，能量计量系统的相对扩展不确定度可控制在1.0%以内，技术上是成熟可靠的，经济上也有明显的公平性收益。（2）关于仪表精度配置。由敏感性分析可知，压力与温度测量对能量计量的影响与流量计同等重要。设计中应避免“重流量计、轻温度压力”的倾向，压力变送器准确度不宜低于0.1级，温度测量允许误差不宜大于±0.5K，并应保证取压点与测温点靠近流量计。（3）关于发热量测定方式。对比计算表明，赋值法会使能量计量系统的不确定度增大至1.13%，超出A级系统限值。因此，本计量站必须配置在线气相色谱分析仪，实时测定天然气发热量；赋值计算仅作为色谱仪故障时的临时措施，并应严格限定使用条件与时间。（4）关于量值溯源与运行管理。计量系统应建立完整的量值溯源链，流量计定期实流校准、温度与压力仪表定期检定、色谱分析仪定期用标准气校