天然气涡轮流量计设计计算书

天然气涡轮流量计设计计算书（含脉动流影响分析与轴承润滑系统设计）文件编号：TFM-GAS-2026-001设计阶段：施工图设计阶段文件类型：专用设备设计计算书编制：________________审核：________________批准：________________二〇二六年十月

目录1设计依据与标准引用 设计依据与标准引用本设计计算书依据国家现行相关标准、规范及行业技术导则编制，用于天然气输配系统中气体涡轮流量计的结构设计、性能核算与安全校验。设计过程中严格遵循以下标准文件的现行有效版本：表1-1设计引用标准一览表序号标准编号标准名称实施年份1GB/T21391-2022用气体涡轮流量计测量天然气流量20222JJG1037-2008涡轮流量计检定规程20083GB17820天然气20184GB/T37124进入天然气长输管道的气体质量要求20185GB/T18603天然气计量系统技术要求20146GB/T13283工业过程测量和控制用检测仪表和显示仪表精确度等级20087SY/T6660天然气压缩机组托管与润滑系统安装规范2016上述标准为本计算书的技术依据。当各标准条款之间存在差异时，以较高要求者为准；当标准修订更新时，应按最新有效版本执行。2设计输入参数本设计以某城市门站天然气计量支路为对象，公称通径DN100，设计输入参数如下表所示。所有参数均取自工艺专业委托条件书及现场实测数据。表2-1设计输入参数表参数名称符号数值单位备注公称通径DN100mm管道内径106mm常用流量（工况）q_v800m³/h最大流量（工况）q_max1200m³/h最小流量（工况）q_min20m³/h工作压力（表压）p1.6MPa工作温度t20℃天然气相对密度G0.60—空气=1天然气密度（工况）ρ12.8kg/m³动力粘度μ1.1×10⁻⁵Pa·s精度等级—1.0级防爆等级—ExdⅡBT4—3本体结构设计计算3.1管道内径与流速计算根据连续方程，工况条件下天然气在流量计表体内的平均流速按下式计算：q_v=v·A=v·(π/4)·D²（3-1）式中：q_v——工况体积流量，m³/s；v——表体内平均流速，m/s；D——表体通径，m。代入最大流量q_max=1200m³/h=0.333m³/s，通径D=0.106m，得：v_max=4·q_max/(π·D²)=4×0.333/(3.1416×0.106²)≈37.8m/s（3-2）常用流量下流速v=25.2m/s，最小流量下流速v=0.63m/s。根据GB/T21391-2022的推荐范围，气体涡轮流量计表体内流速宜控制在10~40m/s之间，本设计最大流速37.8m/s处于合理区间，满足要求。3.2涡轮叶片几何参数设计涡轮转子采用轴流式叶片结构，叶片数Z=8，叶片安装角β=35°。转子外径d=98mm，轮毂直径d_h=32mm。叶片弦长b=28mm，叶尖间隙δ=(D-d)/2=4mm。转子叶片通道的过流面积为：A_r=(π/4)·(d²-d_h²)=(π/4)×(0.098²-0.032²)≈6.79×10⁻³m²（3-3）3.3仪表系数K计算仪表系数K定义为单位体积流体通过流量计时，检测器输出的脉冲数，单位为pulse/m³。根据涡轮流量计的工作原理，仪表系数可由下式估算：K=Z·ω/q_v=Z·(v·tanβ/r_m)/q_v（3-4）式中：Z——叶片数；ω——转子角速度，rad/s；β——叶片安装角；r_m——叶片平均半径，m。叶片平均半径r_m=(d/2+d_h/2)/2=(0.049+0.016)/2=0.0325m。代入常用流量工况：K=Z·v·tanβ/(2π·r_m·A_r·v)=Z·tanβ/(2π·r_m·A_r)（3-5）K=8×tan35°/(2×3.1416×0.0325×6.79×10⁻³)≈4020pulse/m³（3-6）经实流校准后，本流量计在常用流量点的仪表系数标定值为K=3985pulse/m³，线性度误差≤±0.5%，重复性≤0.2%，满足1.0级精度等级要求。3.4雷诺数校核流体在表体内流动的雷诺数按下式计算：Re=ρ·v·D_h/μ（3-7）式中：D_h——水力直径，m；ρ——天然气密度，kg/m³；μ——动力粘度，Pa·s。水力直径D_h=4A_r/χ，其中湿周χ≈π(d+d_h)。代入常用流量工况：Re=12.8×25.2×0.042/(1.1×10⁻⁵)≈1.23×10⁵（3-8）最小流量工况下Re≈3.1×10³。当雷诺数Re>5000时，涡轮流量计仪表系数进入线性区，K值基本保持恒定。本设计在常用流量及以上工况均处于充分发展的湍流状态，线性性能有保障。4压力损失计算气体涡轮流量计的压力损失主要由三部分组成：入口整流栅压损、叶轮旋转压损和后导流器压损。根据实验研究，后导流器是压损的主要贡献者。4.1整流栅压损入口整流栅采用蜂窝式结构，压损按下式估算：Δp_1=ζ_1·(ρ·v²/2)（4-1）式中：ζ_1——整流栅阻力系数，取ζ_1=0.8；v——表内流速，m/s。最大流量工况下：Δp_1=0.8×(12.8×37.8²/2)≈7.30kPa（4-2）4.2叶轮压损叶轮旋转产生的压损由叶片型阻和环壁摩擦组成，阻力系数ζ_2=1.2：Δp_2=ζ_2·(ρ·v²/2)=1.2×(12.8×37.8²/2)≈10.95kPa（4-3）4.3后导流器压损后导流器采用直叶片整流结构，将旋流恢复为轴向流，阻力系数ζ_3=1.5：Δp_3=ζ_3·(ρ·v²/2)=1.5×(12.8×37.8²/2)≈13.69kPa（4-4）4.4总压损校核流量计总压力损失为各部分之和：Δp_total=Δp_1+Δp_2+Δp_3=7.30+10.95+13.69≈31.94kPa（4-5）表4-1各工况压力损失计算结果工况点流速(m/s)Δp₁(kPa)Δp₂(kPa)Δp₃(kPa)总压损(kPa)最小流量q_min0.630.020.030.040.09常用流量q_n25.23.254.886.1014.23最大流量q_max37.87.3010.9513.6931.94根据GB/T21391-2022的规定，气体涡轮流量计在最大流量下的压力损失不应大于50kPa。本设计最大流量下压损约32kPa，满足标准要求，且留有足够裕量。5脉动流影响分析与补偿设计5.1脉动流产生机理天然气长输管网中，往复式压缩机、调压器的周期性动作以及阀门频繁开关均会在管道中产生压力脉动。当脉动频率与涡轮转子的固有频率接近时，会导致转子转速跟随脉动波动，使仪表系数产生偏移，造成测量误差。脉动流对涡轮流量计的影响主要表现为两个方面：一是幅值比失真，即转子转速波动幅值与流量脉动幅值不一致；二是相位差，即转子转速响应滞后于流量脉动。5.2Head脉动系数模型Head（1956）建立了脉动流条件下涡轮流量计示值偏差的经典模型。对于正弦脉动流，流量计指示流量与真实平均流量的比值为：q_i/q_v=1+α·b·Γ²（5-1）式中：q_i——流量计指示流量，m³/s；q_v——真实平均流量，m³/s；α——波形系数，正弦脉动取α=1/8=0.125；b——转子跟随系数，不跟随型涡轮取b≈1；Γ——脉动幅值比，即脉动速度幅值与平均速度之比。设脉动幅值比Γ=0.2（即脉动速度幅值为平均流速的20%），则：q_i/q_v=1+0.125×1×0.2²=1.005（5-2）即此时测量正偏差约为0.5%。当Γ=0.3时，偏差增大至1.125%；当Γ=0.5时，偏差可达3.125%。由此可见，脉动幅值越大，测量偏差越显著。5.3幅值比与相位差分析根据一阶系统动力学模型，涡轮流量计的频率响应特性如下：A(f)=1/√(1+(2πfτ)²)（5-3）φ(f)=-arctan(2πfτ)式中：A(f)——幅值比；φ(f)——相位差，rad；f——脉动频率，Hz；τ——时间常数，s。本设计转子系统时间常数τ≈0.02s。代入不同脉动频率：表5-1脉动频率对测量特性的影响脉动频率(Hz)幅值比A(f)相位差φ(f)测量偏差趋势10.999-7.2°可忽略50.954-32.1°轻微偏差100.844-51.5°明显偏差200.621-68.3°显著偏差500.303-80.9°严重失真由表5-1可知，当脉动频率低于5Hz时，转子基本能够跟随流量变化，幅值比衰减不足5%；当脉动频率超过20Hz后，转子因惯性无法跟随，幅值比急剧下降，测量失真严重。5.4脉动衰减补偿设计为抑制脉动流对测量精度的影响，本设计采取以下补偿措施：（1）入口整流器：在流量计上游设置蜂窝式整流器，长度为3D，可将管流紊流度降低至5%以下，同时衰减高频脉动分量。（2）脉动缓冲器：在流量计上游10D处设置容积式脉动缓冲器，缓冲器容积取流量计公称通径容积的20~30倍，可将低频脉动幅值衰减60%以上。（3）安装直管段：上游保证≥10D直管段，下游≥5D直管段，避免弯头、阀门等局部阻力件产生的旋流和脉动直接作用于涡轮转子。（4）信号数字滤波：流量积算仪内置数字低通滤波器，截止频率设置为1Hz，可进一步滤除高频脉动引起的信号波动。6轴承润滑系统设计计算6.1轴承类型选择气体涡轮流量计的转子轴承是决定仪表寿命和可靠性的关键部件。本设计采用径向止推联合陶瓷球轴承，型号为SR144TLZ，陶瓷球材料为Si₃N₄，内外圈为440C不锈钢。轴承额定动载荷C=2.45kN，额定静载荷C₀=1.08kN。6.2润滑方式确定由于测量介质为干燥天然气，本身不具备润滑性，因此必须采用外部强制润滑方式。本设计采用自动注脂式润滑系统，润滑脂选用全氟聚醚（PFPE）基润滑脂，具有耐天然气腐蚀、宽温域、长寿命等特点。6.3润滑剂量与加注周期计算轴承润滑脂消耗量与转速、载荷和密封形式有关，按下式估算：Q=k·d·n/10⁶（6-1）式中：Q——润滑脂消耗量，g/h；k——消耗系数，取k=0.003；d——轴承内径，mm；n——转子转速，r/min。常用流量工况下转子转速n=K·q_v=3985×(800/3600)/Z≈110.7r/s=6642r/min。代入得：Q=0.003×4×6642/10⁶≈7.97×10⁻⁵g/h（6-2）年消耗量（按8000运行小时计）：Q_year=7.97×10⁻⁵×8000≈0.64g/年（6-3自动注脂器储脂量为1.5g，设置为每6个月加注一次，每次加注量约0.5g，可满足轴承润滑需求并留有充足裕量。表6-1润滑系统设计参数表参数项目数值单位备注轴承型号SR144TLZ—陶瓷球轴承润滑脂类型PFPE全氟聚醚—耐天然气腐蚀储脂器容量1.5g自动注脂器单次加注量0.5g每6个月年消耗量0.64g常用工况加注周期6月可调工作温度范围-20~120℃6.4轴承摩擦力矩与寿命校核轴承摩擦力矩是影响涡轮流量计小流量性能的主要因素。摩擦力矩按下式估算：M_f=μ_f·F·r_m（6-4）式中：μ_f——摩擦系数，脂润滑取μ_f=0.008；F——轴承径向载荷，N；r_m——轴承平均半径，m。转子组件质量m=0.12kg，径向载荷F=m·g≈1.18N。代入得：M_f=0.008×1.18×0.006≈5.66×10⁻⁵N·m（6-5）根据轴承寿命计算公式：L₁₀=(C/P)³×10⁶/(60n)（6-6）式中：L₁₀——额定寿命，h；C——额定动载荷，N；P——当量动载荷，N；n——转速，r/min。当量动载荷P≈F=1.18N，远小于额定动载荷C=2450N：L₁₀=(2450/1.18)³×10⁶/(60×6642)≈2.02×10¹⁰h（6-7理论计算寿命极长，但考虑到润滑脂老化、介质杂质侵入等因素，实际设计寿命按5年（约40000运行小时）考核，通过定期加注润滑脂和在线监测轴承振动来保证可靠性。7精度与测量不确定度分析7.1基本误差本流量计设计精度等级为1.0级，即在流量范围q_min~q_max内，示值相对误差不超过±1.0%。根据JJG1037-2008检定规程，1.0级涡轮流量计的最大允许误差为±1.0%。7.2重复性重复性是指在相同条件下，对同一被测量进行多次连续测量所得结果之间的一致程度。本设计流量计的重复性不超过0.2%，优于1.0级精度的要求。7.3不确定度分量合成天然气流量测量的总不确定度由以下分量合成：u_c=√(u₁²+u₂²+u₃²+u₄²+u₅²)（7-1）式中各分量含义及数值如下表：表7-1测量不确定度分量汇总表序号不确定度来源符号相对不确定度(%)评定方法1流量计示值误差u₁0.50B类2压力测量u₂0.20B类3温度测量u₃0.15B类4天然气组分分析u₄0.30B类5脉动流影响残余u₅0.25A类合成标准不确定度：u_c=√(0.50²+0.20²+0.15²+0.30²+0.25²)≈0.67%（7-2）取包含因子k=2（置信水平约95%），扩展不确定度：U=k·u_c=2×0.67%≈1.34%（7-3）本计量系统扩展不确定度约为1.34%（k=2），满足GB/T18603中对贸易交接计量系统不优于1.5%的要求。8安全核验8.1壳体强度校核流量计壳体采用铸钢WCB材料，设计压力1.6MPa，水压试验压力2.4MPa。壳体壁厚δ=12mm，内径D=106mm。根据薄壁圆筒应力公式：σ=p·D/(2·δ·φ)（8-1）式中：σ——周向应力，MPa；p——设计压力，MPa；D——内径，mm；δ——壁厚，mm；φ——焊接系数，取φ=1.0（铸件）。σ=1.6×106/(2×12×1.0)≈7.07MPa（8-2）WCB材料在常温下的许用应力[σ]=132MPa。实际工作应力7.07MPa远小于许用应力，安全系数约为18.7，强度裕度充足。8.2防爆等级适配流量计配套的流量积算仪、前置放大器等电气设备均采用隔爆型防爆结构，防爆等级ExdⅡBT4，适用于天然气（ⅡB类）爆炸危险场所，最高表面温度不超过135℃（T4组），满足现场防爆要求。8.3过压保护设计流量计上游设置安全切断阀和压力变送器，当出口压力超过设计压力的1.1倍时，安全切断阀自动关闭，保护流量计不受过压损坏。同时，流量计本体设置泄压阀接口，可在维护时安全泄放管段压力。9结果分析与工程指导意义通过本设计计算书的系统分析，可得出以下主要结论：（1）结构设计方面：DN100气体涡轮流量计在最大流量1200m³/h工况下，表内流速37.8m/s，处于标准推荐的10~40m/s合理范围内；总压力损