简介一份关于PSA变压吸附制氧设备设计的专业计算书由上海施耐德日盛气体设备有限公司编制适合从事气体分离、制氧设备设计及化工工艺计算的工程师与相关专业学生参考。文档以15Nm³/h产气量、93%纯度的制氧需求为例完整呈现从压缩空气耗气量、系统主管路管径、空压机排气量到空气净化组件选型、5A分子筛装填量、吸附塔直径与结构尺寸的设计计算过程。其中包含德国BF5A分子筛性能参数、基于氧气回收率的流量换算公式、空塔线速度校核等关键步骤并给出裕量选取与标准圆整原则。全文为单个doc格式文件共1个文档容量仅116KB便于直接编辑复用或按实际工况替换参数重新演算。截至目前已有289人学习适合作为PSA制氧工程计算模板或课程设计参考。1. PSA制氧设计计算书为什么先从15Nm³/h的案例看起这份由上海施耐德日盛气体设备完成的《PSA制氧设计计算书》开篇就把设计边界钉死产气量15Nm3/h、氧气纯度93%、工作压力0.5MPa、温度20℃分子筛选用德国BF5A型。变压吸附制氧有个特点设备外观看起来是几个罐和一堆阀门但真正的设计难度全在看不见的算量里——压缩空气耗气量决定空压机和净化组件5A分子筛装填量决定吸附塔直径和高度缓冲罐容积决定切换时的压力稳定性。看懂这份计算书的人也就掌握了中小型PSA制氧从工艺参数到设备尺寸的完整推导链条对气体设备设计、工程报价、现场设备调试都有直接帮助。下面按计算书的顺序把每个公式、系数来源和容易踩的坑逐一拆开。2. 压缩空气耗气量与主管路管径先把0.5MPa下的空气量算明白2.1 从氧气需求反推压缩空气流量 QaPSA制氧的原料是压缩空气空气中氧含量固定为20.9%但不可能全部回收。计算书先给出氧气产量、纯度和氧气回收率反推需要的压缩空气量Qa (Qn × C) / (20.9% × η)代入Qn15 Nm3/h、C93%、η35.2%得到Qa (15 × 0.93) / (0.209 × 0.352) ≈ 190 Nm3/h这个式子的物理意义很直白产品气里的氧量是15 × 0.93 13.95 Nm3/h空气中的氧含量比例是20.9%但整套吸附流程只能收回其中35.2%所以原料空气量要放大到约13.95 / (0.209 × 0.352)。回收率η不是固定值它与吸附压力、均压时间、分子筛性能、切换阀动作速度都有关系。如果厂家不能提供实测η前期估算可先取30%40%的区间偏向保守。还有一点要提醒C和20.9%都要写成小数写成93和20.9会让结果差出两个数量级我见过不止一次因为纯度漏除100而导致空压机选错的案例。2.2 用绝压和温度修正把标况流量换算到工况流量190Nm3/h是标态干空气流量但管道直径、吸附塔内流速必须用实际工况下的体积流量。理想气体状态方程做压力、温度修正Qa实 Qa × (PN / PS) × (TS / TN)PN0.101MPa是标态绝压PSPN气源压力0.1010.50.601MPa。当现场环境温度等于计算基准温度20℃时TS/TN1只保留压力修正Qa实 190 × 0.101 / 0.601 ≈ 31.93 m3/h这里最容易出错的是绝压和表压混用。气源压力0.5MPa是表压必须加上当地大气压0.101MPa直接拿0.5当分母Qa实会算成23.9m3/h后面所有跟流速有关的尺寸都会偏小。如果设备放在海拔1500米以上的地区PN要按当地大气压修正不能继续用0.101。温度修正也不是永远可以忽略气源温度比20℃高15℃以上时TS/TN会带来约5%的流量变化管线较长或空压机后冷却器效果不佳时要重新算。2.3 主管路管径 da 与气体流速选择有了工况流量管道内径用连续性方程计算da sqrt( (4 × Qa实 / 3600) / (π × υ) )Qa实31.93m3/h取气体流速υ8m/s代入得到da≈37.6mm圆整为DN40。气体流速是设计者可以主动选择的参数压缩空气总管推荐812m/s低压侧或短管可放到15m/s但流速越高压降越大长期运行电费会明显增加。下面这个表给出同一流量在不同流速下的理论内径流速 (m/s)681012理论内径 (mm)43.437.633.630.7推荐标准管径 DN50404032同样31.93m3/h取6m/s时DN50更稳妥但阀门和管件成本上升取12m/s时DN32即可但阻力大且阀门节流时噪音明显。选择8m/s折中是常见做法。计算书最终取DN40工程上没有争议。实际配管还要考虑长度、弯头数量和未来扩容一般会在理论管径上至少保持一档余量也就是说理论37.6mm按DN40选而不是按DN32。用代码把这步封装起来方便后续做多流速方案对比import math def pipe_diameter(q_actual_m3h, velocity8.0): # 工况流量从 m3/h 转成 m3/s q q_actual_m3h / 3600.0 # 截面积 流量 / 流速单位 m2 cross_area q / velocity # 直径 sqrt(4*A/pi)最后转成 mm d math.sqrt(4 * cross_area / math.pi) * 1000 return d qa_actual 190 * 0.101 / (0.101 0.5) d pipe_diameter(qa_actual, velocity8.0) print(d) # 输出: 37.6这个函数参数q_actual_m3h是工况流量velocity是设计流速返回单位是mm。函数里先把m3/h除3600变成m3/s再按截面积反算直径。当需要比较多种流速时只需要循环传入velocity即可。要注意的是代码里0.101和0.5分别是绝压和表压修改时要保持一致如果项目中有富余系数可以在q_actual_m3h上先乘1.1再送到函数里。2.4 空压机排气量折算不能直接用190Nm3/h选机空压机铭牌排气量通常是吸入状态下的湿空气体积而PSA设计基准是标态干空气。要把190Nm3/h换算到空压机吸入口状态计算书给出了一个工程公式V空 Qa × T0 / [264 × (P0 - ψ × Pb)]P01.01kgf/cm2T0293Kψ65%Pb20℃饱和水蒸汽分压0.0238kgf/cm2。代入V空 190×293 / [264×(1.01 - 0.65×0.0238)] ≈ 212.04 m3/h 3.54 m3/min这个264是单位换算经验常数把气压单位kgf/cm2、温度K、时间h整合成一套系数常温常压下可直接用。如果项目地大气压偏低或环境温度经常超过35℃建议把当地气象参数代入重新折算而不是沿用0.0238这个Pb。计算出3.54m3/min之后选空压机还要考虑分子筛再生、充压和阀门泄漏一般按10%15%服务系数放大选4.0m3/min左右。原计算书在这个公式的代号里出现过Qa实26.05m3/h的笔误复核应以Qa190Nm3/h和压力比0.101/0.601为准这不影响最终结论但提醒我们拿到任何设计计算书都要自己算一遍。2.5 回收率η变化时整套设计如何联动前面所有计算都把η固定为35.2%但实际运行时η可能因为均压时间、阀门内漏、分子筛老化而变化。假设其他条件不变只改ηηQa (Nm3/h)Qa实 (m3/h)理论管径 (mm)建议管径30%222.537.440.7DN5035.2%19031.937.6DN4040%166.928.035.2DN40可见η从35.2%降到30%压缩空气耗量增加17%主管路理论管径从37.6mm升到40.7mm标准管径就需要从DN40跳到DN50。对于已经建成的设备η下降最直接的后果是空压机排气量不够、吸附塔切换时压力下跌产品氧纯度跟着波动。因此设计阶段不要把η取得过于乐观如果分子筛已运行多年建议按30%重新校核一次空压机和管路裕量。3. 5A分子筛装填量与吸附塔直径从300kg到Φ550mm的关键校核3.1 理论装填量和实际装填量的差距在哪里分子筛装填量直接决定吸附塔的成本和处理能力。计算书使用产氧率K作为核心参数G分 Qn / K 15 / 50 0.3 (T) 300kgK50NL/kg的含义是在0.5MPa、20℃条件下每公斤5A分子筛每小时可产出50NL标升氧气。产品氧为15Nm3/h即15000NL/h所以理论装填量300kg。这里特别要注意K值不是常数值它随吸附压力升高而上升随温度升高而下降。计算书选用德国BF5A分子筛K值来自特定工况的实测数据如果换用国产5A或13X必须重新获取吸附等温线和产氧率数据。计算书在理论值上乘以1.1得到实际装填量330kg。这个10%的增量是工程经验的浓缩分子筛使用中会粉化粉体堵塞气流通道导致有效吸附容量下降压缩空气即使经过净化组件仍可能存在微量水分和油雾会使分子筛逐渐中毒吸附塔切换时气流冲刷壁面附近分子筛利用率偏低。如果原料气质量不稳定或现场维护条件较差这个富裕量可以提高到1.151.2但不要随意降低否则纯度衰减后很难补救。两塔交替运行时单塔装填量是总装填量的一半即165kg。后续的塔体尺寸计算都以单塔装填量为基础。3.2 单塔容积、堆密度单位与直径的初始估算单塔分子筛容积由质量除以堆密度得到。计算书给出的堆密度是670kg/L这个单位明显不对5A分子筛的堆密度通常在650700kg/m3范围内也就是0.650.70kg/L。如果真按670kg/L计算330kg分子筛只有0.49L体积吸附塔不可能装得下。正确换算V分 (330 / 2) / 670 165 / 670 ≈ 0.246 m3工程上取0.25m3作为单塔装填容积。这个密度单位错误其实是计算书打印时最常见的坑很多人在做设备复核时被带偏。统一的建议是所有物料计算都用kg/m3遇到kg/L一律先×1000。塔径需要满足空塔线速度的限制。立式PSA吸附塔的经验做法是空塔空气线速度ω≤0.045m/s高径比H:D≥2。空塔线速度是假定塔内没有分子筛、气体均匀通过整个塔截面时的线速度速度过高会导致气流在床层内偏流、压降大速度过低则塔径偏大、分子筛利用率下降。按0.045m/s取D sqrt(4 × Qa实 / (π × ω))Qa实31.93m3/h0.00887m3/s得到D≈0.501m圆整为Φ550mm。3.3 空塔线速度复核与高径比验证原计算书里的一处矛盾圆整后要回头校核线速度。Φ550mm塔截面积为π×0.55²/4≈0.2376m2实际线速度ω实 0.00887 / 0.2376 ≈ 0.0373 m/s低于0.045m/s说明选550mm是安全的。如果选Φ500mm截面积0.1963m2线速度约0.0452m/s已经略超上限在压力波动时容易扰动床层。所以Φ550mm是合理选择。接下来算装填高度。单塔容积0.25m3塔径Φ550mm可算出Ha 4V分 / (πD²) 4×0.25 / (π×0.55²) ≈ 1.05m高径比Ha:D≈1.91没有达到计算书要求的≥2。而原计算书写的是Ha1.58m高径比2.87反推其塔径约Φ450mm与前面圆整后的Φ550mm明显矛盾。这个矛盾在工程上很典型先按计算直径算出初步尺寸圆整标准直径后忘记重新迭代装填高度。正确做法是直接按标准直径Φ550mm反算装填高度再验证高径比若高径比不足1.91应把塔径改小一档到Φ500mm此时装填高度约1.27m高径比2.54同时空塔线速度0.0452m/s仅比额定值高一点配合更均匀的布料器也可以接受。所以实际设计这里要在两个方案里权衡Φ500mm满足高径比但线速度临界Φ550mm线速度余量大但高径比不足。很多制造厂最终会选Φ500mm并增加气流分布器。3.4 吸附塔内部结构尺寸与进出口管径吸附塔不是简单一个空罐子装满分子筛就完事。计算书把塔内纵向分成三个区域A区装3A分子筛用来吸附残余水分和二氧化碳保护主床层的5A分子筛B区装氧化铝干燥剂拦截水分并与5A分子筛分层装填C区是分子筛压紧附件通常用弹簧压板或气囊把床层压住防止气流冲击导致分子筛摩擦粉化。A区要在计算装填高度Ha基础上加0.05m修正然后依次加上h250mm、h3100mm左右的安装高度最终直边高H约1.81m。进出管径的确定逻辑是压缩空气入口管与空气净化组件后的主管路保持一致取d140mm因为PSA塔进口前没有其他节流装置氧气出口管先按产品气工况流量2.52m3/h计算理论内径约10.5mm考虑吸附塔升压时产品气回流和调阀需要一定余量取d212mm比空气进口小一档是合理的。这里注意氧气出口管径不能因为算出来细就随意缩小还要看塔内分布帽的开孔面积保证气流能均匀进入床层。3.5 用Python快速跑一台吸附塔的分子筛量分子筛装填量和塔径可以汇总成下面这个函数用于后续不同纯度、不同产量的方案对比def adsorber_calc(qn, qa_actual, K50, margin1.1, rho670): # qn 产品氧气流量 Nm3/hqa_actual 工况压缩空气流量 m3/h g_real qn * 1000 / K * margin # 总实际装填量 kg v_single g_real / 2 / rho # 单塔容积 m3 q_flow qa_actual / 3600.0 d_ads math.sqrt(4 * q_flow / (math.pi * 0.045)) * 1000 d_std math.ceil(d_ads / 50) * 50 # 标准塔径50mm一档 h_fill 4 * v_single / (math.pi * (d_std / 1000) ** 2) return {sieve_mass: g_real, sieve_volume: v_single, d_std: d_std, h_fill: h_fill} print(adsorber_calc(15, 31.93)) # 输出: {sieve_mass: 330.0, sieve_volume: 0.246, d_std: 550, h_fill: 1.05}这里g_real是总实际装填量v_single是单塔容积d_std按50mm一档向上取整得到标准塔径h_fill是用标准塔径反算的装填高度避免出现前面那种混用直径的情况。注意qn要传Nm3/hqa_actual要传工况实际流量两者单位不同不要混在一起。输出结果1.05m的装填高度已经暴露了高径比不足所以实际产品文档里的1.58m只能理解为制造图与计算书版本不一致复核时一定要以这个函数重算结果为准。4. 空气与氧气缓冲罐容积计算缓冲系数Ka怎么取才不掉链子4.1 为什么要两个缓冲罐Ka的物理含义吸附塔在切换瞬间进气阀和排气阀同时动作管线内流量会发生阶跃变化。空气缓冲罐的作用是把空压机连续输出的空气“缓存”起来在吸附塔充压时补充瞬时流量氧气缓冲罐则用来平滑产品气压力和纯度保证下游用户在切换期间不中断供气。两者容积都用缓冲系数Ka乘以对应流量求得但Ka的取值差别很大这不是随意定的而是由两侧的流量波动特性决定的。4.2 空气缓冲罐容积初算与结构反推空气缓冲罐初算容积V缓 Ka × Qa 0.011 × 190 2.09 m3这里的190是压缩空气耗气量单位Nm3/h不是氧气产量。Ka取0.011h相当于要给空压机供气侧约40秒的缓冲时间。有些设计人员想把空气缓冲罐做小一点把Ka降到0.008h结果吸附塔充压时进气压力掉压超过0.03MPa分子筛解吸不充分。因此Ka取值最好落在0.0100.015h区间并根据吸附塔容积和切换频率做微调。初取直边高1810mm标准椭圆封头估算直径D缓 sqrt(4×2.09 / (π×1.81)) ≈ 1.22m圆整为DN1200mm。然后反算直边高H缓 4×(V缓 - 2V封) / (πD²) 4×(2.09 - 2×0.225) / (π×1.2²) ≈ 1.45m圆整为1.5m。标准椭圆封头在DN1200下单个容积约0.225m3这部分数据来自GB/T 25198不能只用圆柱体积公式。反算完直边高后还要做最终校核V缓 π×1.2²×1.5/4 2×0.225 ≈ 2.15m3Ka实 2.15 / 190 ≈ 0.0113h与初取值基本吻合。整个流程的关键是先估直径-圆整-反算直边高-再校核Ka而不是直接用初算容积去画图。4.3 氧气缓冲罐的Ka为什么大一圈氧气缓冲罐初算容积V缓 Ka × Qn 0.031 × 15 0.465 m3Ka取0.0300.04h流量基准是产品氧15Nm3/h。计算书里有一行推导Ka(5.5/0.5)/(60×6)再乘绝压比(0.601/0.101)意思是吸附塔每小时切换6次每次切换产生约60秒的产气中断压力比越大需要的缓冲时间越长。虽然实际数字不一定精确对应但它点明了Ka的本质切换频率越高、吸附压力越高、产品气允许波动越小Ka就要取得越大。如果吸附塔采用3塔流程切换错开Ka可以适当取小2塔流程则必须取大值。按同样的结构算法初取直边高1810mm估算直径0.57m圆整为DN600mm。DN600椭圆封头容积约0.0352m3反算直边高约为1.4m最终总容积V缓 π×0.6²×1.4/4 2×0.0352 ≈ 0.466m3Ka实0.466/15≈0.031h符合预期。产品气流量越小罐体绝对容积越小但相对缓冲时间不能压缩这就是为什么15Nm3/h的小设备也要有一个直径600mm的氧气缓冲罐。4.4 两组缓冲罐参数对比与现场验证把两组计算放在一起看更清晰项目空气缓冲罐氧气缓冲罐缓冲系数Ka (h)0.0110.031流量基准压缩空气190 Nm3/h产品氧15 Nm3/h初算容积 (m3)2.090.465标准直径 DN (mm)1200600直边高 (mm)15001400实际总容积 (m3)2.150.466校核Ka (h)0.01130.031从表里能看到氧气罐的绝对容积只有空气罐的约1/4但它的Ka是空气罐的近3倍说明它承担的是高频小流量波动空气罐承担的是低频大流量波动。现场验证方法很简单在吸附塔切换的一瞬间看缓冲罐压力表空气罐压力下跌不超过0.03MPa、氧气罐压力波动不超过3%可以认为容积足够如果压力波动超限先检查切换阀开闭时间再考虑增加容积。还有一个容易犯的错Ka带时间量纲不要把它理解成“容积系数”。0.031×150.465m3这里的15是Nm3/h不是m3/h如果换成工况流量2.52m3/h去乘Ka就会变成0.184h罐体容积被放大约5倍。正确做法是设计和校核都统一用同一个流量基准并在计算书里写明取的是哪个流量。5. 把设计计算书写成可复用的Python小工具再避开四个工程坑5.1 一个覆盖主流程的粗算脚本前面每一步都可以单独手算但方案对比时会很慢。把主流程封装成一个函数输入氧气产量和纯度输出空气量、分子筛量、管径和塔径。下面是一个20行的版本def psa_design(qn, purity, eta0.352, p_gauge0.5, K50, rho_sieve670, omega0.045, velocity8.0): # qn 产品氧流量 Nm3/hpurity 氧纯度小数 qa qn * purity / (0.209 * eta) # 标态空气耗量 Nm3/h qa_act qa * 0.101 / (0.101 p_gauge) # 工况空气流量 m3/h d_pipe math.sqrt(4 * (qa_act / 3600) / (math.pi * velocity)) * 1000 g_ads qn * 1000 / K * 1.1 # 实际分子筛总装填量 kg v_sieve g_ads / 2 / rho_sieve # 单塔筛容积 m3 d_ads math.sqrt(4 * (qa_act / 3600) / (math.pi * omega)) * 1000 return { qa_标态: qa, qa_工况: qa_act, 主管内径mm: round(d_pipe), 分子筛装填量kg: round(g_ads), 单塔筛容积m3: round(v_sieve, 3), 吸附塔内径mm: math.ceil(d_ads / 50) * 50, } print(psa_design(15, 0.93)) # {qa_标态: 189.6, qa_工况: 31.9, 主管内径mm: 38, 分子筛装填量kg: 330, 单塔筛容积m3: 0.246, 吸附塔内径mm: 550}代码里的p_gauge必须传表压函数内部已经完成绝压修正qn是标况产品氧气流量purity用小数。返回的“主管内径mm”是理论值实际配管还要向上圆整到标准管径“吸附塔内径mm”已经做了50mm一档向上取整。当前参数下分子筛量330kg、塔径550mm和计算书主线一致。5.2 四个工程坑第一分子筛堆密度单位。计算书写的670kg/L是笔误真实单位是kg/m3任何质量除以密度算容积时都要先确认单位。第二吸附塔直径圆整后必须重新反算装填高度。原计算书Ha1.58m与Φ550mm塔径不匹配用标准直径反算出的1.05m才符合几何关系高径比不足时优先改塔径而不是改Ha。第三缓冲系数的流量基准不能混。空气罐乘压缩空气量氧气罐乘产品氧量两者差一个回收率折算。换用工况流量时Ka的数值必须同步重新设定不能只改一个乘数。第四K和η是工况实测值。0.5MPa、20℃下测得的50NL/kg和35.2%换到0.7MPa或40℃环境就不成立。方案对比可以用经验值正式设计一定要向分子筛厂家索取该温度压力下的吸附曲线并按吸附塔实际循环周期重新校核产氧率。遇到旧计算书里的矛盾不要急着否定整个方案。把几何尺寸独立复核一遍再拿现场半小时的压力和纯度记录对一下通常能判断是版本迭代还是抄写错误。对气体设备这类项目建议把几何尺寸、物料质量、阀门切换逻辑分开归档任何一次设计变更都重新跑一遍上面的粗算脚本再改图纸。本文还有配套的精品资源点击获取 PSA制氧设计计算书拆解：15Nm³/h案例从空气量到缓冲罐
简介一份关于PSA变压吸附制氧设备设计的专业计算书由上海施耐德日盛气体设备有限公司编制适合从事气体分离、制氧设备设计及化工工艺计算的工程师与相关专业学生参考。文档以15Nm³/h产气量、93%纯度的制氧需求为例完整呈现从压缩空气耗气量、系统主管路管径、空压机排气量到空气净化组件选型、5A分子筛装填量、吸附塔直径与结构尺寸的设计计算过程。其中包含德国BF5A分子筛性能参数、基于氧气回收率的流量换算公式、空塔线速度校核等关键步骤并给出裕量选取与标准圆整原则。全文为单个doc格式文件共1个文档容量仅116KB便于直接编辑复用或按实际工况替换参数重新演算。截至目前已有289人学习适合作为PSA制氧工程计算模板或课程设计参考。1. PSA制氧设计计算书为什么先从15Nm³/h的案例看起这份由上海施耐德日盛气体设备完成的《PSA制氧设计计算书》开篇就把设计边界钉死产气量15Nm3/h、氧气纯度93%、工作压力0.5MPa、温度20℃分子筛选用德国BF5A型。变压吸附制氧有个特点设备外观看起来是几个罐和一堆阀门但真正的设计难度全在看不见的算量里——压缩空气耗气量决定空压机和净化组件5A分子筛装填量决定吸附塔直径和高度缓冲罐容积决定切换时的压力稳定性。看懂这份计算书的人也就掌握了中小型PSA制氧从工艺参数到设备尺寸的完整推导链条对气体设备设计、工程报价、现场设备调试都有直接帮助。下面按计算书的顺序把每个公式、系数来源和容易踩的坑逐一拆开。2. 压缩空气耗气量与主管路管径先把0.5MPa下的空气量算明白2.1 从氧气需求反推压缩空气流量 QaPSA制氧的原料是压缩空气空气中氧含量固定为20.9%但不可能全部回收。计算书先给出氧气产量、纯度和氧气回收率反推需要的压缩空气量Qa (Qn × C) / (20.9% × η)代入Qn15 Nm3/h、C93%、η35.2%得到Qa (15 × 0.93) / (0.209 × 0.352) ≈ 190 Nm3/h这个式子的物理意义很直白产品气里的氧量是15 × 0.93 13.95 Nm3/h空气中的氧含量比例是20.9%但整套吸附流程只能收回其中35.2%所以原料空气量要放大到约13.95 / (0.209 × 0.352)。回收率η不是固定值它与吸附压力、均压时间、分子筛性能、切换阀动作速度都有关系。如果厂家不能提供实测η前期估算可先取30%40%的区间偏向保守。还有一点要提醒C和20.9%都要写成小数写成93和20.9会让结果差出两个数量级我见过不止一次因为纯度漏除100而导致空压机选错的案例。2.2 用绝压和温度修正把标况流量换算到工况流量190Nm3/h是标态干空气流量但管道直径、吸附塔内流速必须用实际工况下的体积流量。理想气体状态方程做压力、温度修正Qa实 Qa × (PN / PS) × (TS / TN)PN0.101MPa是标态绝压PSPN气源压力0.1010.50.601MPa。当现场环境温度等于计算基准温度20℃时TS/TN1只保留压力修正Qa实 190 × 0.101 / 0.601 ≈ 31.93 m3/h这里最容易出错的是绝压和表压混用。气源压力0.5MPa是表压必须加上当地大气压0.101MPa直接拿0.5当分母Qa实会算成23.9m3/h后面所有跟流速有关的尺寸都会偏小。如果设备放在海拔1500米以上的地区PN要按当地大气压修正不能继续用0.101。温度修正也不是永远可以忽略气源温度比20℃高15℃以上时TS/TN会带来约5%的流量变化管线较长或空压机后冷却器效果不佳时要重新算。2.3 主管路管径 da 与气体流速选择有了工况流量管道内径用连续性方程计算da sqrt( (4 × Qa实 / 3600) / (π × υ) )Qa实31.93m3/h取气体流速υ8m/s代入得到da≈37.6mm圆整为DN40。气体流速是设计者可以主动选择的参数压缩空气总管推荐812m/s低压侧或短管可放到15m/s但流速越高压降越大长期运行电费会明显增加。下面这个表给出同一流量在不同流速下的理论内径流速 (m/s)681012理论内径 (mm)43.437.633.630.7推荐标准管径 DN50404032同样31.93m3/h取6m/s时DN50更稳妥但阀门和管件成本上升取12m/s时DN32即可但阻力大且阀门节流时噪音明显。选择8m/s折中是常见做法。计算书最终取DN40工程上没有争议。实际配管还要考虑长度、弯头数量和未来扩容一般会在理论管径上至少保持一档余量也就是说理论37.6mm按DN40选而不是按DN32。用代码把这步封装起来方便后续做多流速方案对比import math def pipe_diameter(q_actual_m3h, velocity8.0): # 工况流量从 m3/h 转成 m3/s q q_actual_m3h / 3600.0 # 截面积 流量 / 流速单位 m2 cross_area q / velocity # 直径 sqrt(4*A/pi)最后转成 mm d math.sqrt(4 * cross_area / math.pi) * 1000 return d qa_actual 190 * 0.101 / (0.101 0.5) d pipe_diameter(qa_actual, velocity8.0) print(d) # 输出: 37.6这个函数参数q_actual_m3h是工况流量velocity是设计流速返回单位是mm。函数里先把m3/h除3600变成m3/s再按截面积反算直径。当需要比较多种流速时只需要循环传入velocity即可。要注意的是代码里0.101和0.5分别是绝压和表压修改时要保持一致如果项目中有富余系数可以在q_actual_m3h上先乘1.1再送到函数里。2.4 空压机排气量折算不能直接用190Nm3/h选机空压机铭牌排气量通常是吸入状态下的湿空气体积而PSA设计基准是标态干空气。要把190Nm3/h换算到空压机吸入口状态计算书给出了一个工程公式V空 Qa × T0 / [264 × (P0 - ψ × Pb)]P01.01kgf/cm2T0293Kψ65%Pb20℃饱和水蒸汽分压0.0238kgf/cm2。代入V空 190×293 / [264×(1.01 - 0.65×0.0238)] ≈ 212.04 m3/h 3.54 m3/min这个264是单位换算经验常数把气压单位kgf/cm2、温度K、时间h整合成一套系数常温常压下可直接用。如果项目地大气压偏低或环境温度经常超过35℃建议把当地气象参数代入重新折算而不是沿用0.0238这个Pb。计算出3.54m3/min之后选空压机还要考虑分子筛再生、充压和阀门泄漏一般按10%15%服务系数放大选4.0m3/min左右。原计算书在这个公式的代号里出现过Qa实26.05m3/h的笔误复核应以Qa190Nm3/h和压力比0.101/0.601为准这不影响最终结论但提醒我们拿到任何设计计算书都要自己算一遍。2.5 回收率η变化时整套设计如何联动前面所有计算都把η固定为35.2%但实际运行时η可能因为均压时间、阀门内漏、分子筛老化而变化。假设其他条件不变只改ηηQa (Nm3/h)Qa实 (m3/h)理论管径 (mm)建议管径30%222.537.440.7DN5035.2%19031.937.6DN4040%166.928.035.2DN40可见η从35.2%降到30%压缩空气耗量增加17%主管路理论管径从37.6mm升到40.7mm标准管径就需要从DN40跳到DN50。对于已经建成的设备η下降最直接的后果是空压机排气量不够、吸附塔切换时压力下跌产品氧纯度跟着波动。因此设计阶段不要把η取得过于乐观如果分子筛已运行多年建议按30%重新校核一次空压机和管路裕量。3. 5A分子筛装填量与吸附塔直径从300kg到Φ550mm的关键校核3.1 理论装填量和实际装填量的差距在哪里分子筛装填量直接决定吸附塔的成本和处理能力。计算书使用产氧率K作为核心参数G分 Qn / K 15 / 50 0.3 (T) 300kgK50NL/kg的含义是在0.5MPa、20℃条件下每公斤5A分子筛每小时可产出50NL标升氧气。产品氧为15Nm3/h即15000NL/h所以理论装填量300kg。这里特别要注意K值不是常数值它随吸附压力升高而上升随温度升高而下降。计算书选用德国BF5A分子筛K值来自特定工况的实测数据如果换用国产5A或13X必须重新获取吸附等温线和产氧率数据。计算书在理论值上乘以1.1得到实际装填量330kg。这个10%的增量是工程经验的浓缩分子筛使用中会粉化粉体堵塞气流通道导致有效吸附容量下降压缩空气即使经过净化组件仍可能存在微量水分和油雾会使分子筛逐渐中毒吸附塔切换时气流冲刷壁面附近分子筛利用率偏低。如果原料气质量不稳定或现场维护条件较差这个富裕量可以提高到1.151.2但不要随意降低否则纯度衰减后很难补救。两塔交替运行时单塔装填量是总装填量的一半即165kg。后续的塔体尺寸计算都以单塔装填量为基础。3.2 单塔容积、堆密度单位与直径的初始估算单塔分子筛容积由质量除以堆密度得到。计算书给出的堆密度是670kg/L这个单位明显不对5A分子筛的堆密度通常在650700kg/m3范围内也就是0.650.70kg/L。如果真按670kg/L计算330kg分子筛只有0.49L体积吸附塔不可能装得下。正确换算V分 (330 / 2) / 670 165 / 670 ≈ 0.246 m3工程上取0.25m3作为单塔装填容积。这个密度单位错误其实是计算书打印时最常见的坑很多人在做设备复核时被带偏。统一的建议是所有物料计算都用kg/m3遇到kg/L一律先×1000。塔径需要满足空塔线速度的限制。立式PSA吸附塔的经验做法是空塔空气线速度ω≤0.045m/s高径比H:D≥2。空塔线速度是假定塔内没有分子筛、气体均匀通过整个塔截面时的线速度速度过高会导致气流在床层内偏流、压降大速度过低则塔径偏大、分子筛利用率下降。按0.045m/s取D sqrt(4 × Qa实 / (π × ω))Qa实31.93m3/h0.00887m3/s得到D≈0.501m圆整为Φ550mm。3.3 空塔线速度复核与高径比验证原计算书里的一处矛盾圆整后要回头校核线速度。Φ550mm塔截面积为π×0.55²/4≈0.2376m2实际线速度ω实 0.00887 / 0.2376 ≈ 0.0373 m/s低于0.045m/s说明选550mm是安全的。如果选Φ500mm截面积0.1963m2线速度约0.0452m/s已经略超上限在压力波动时容易扰动床层。所以Φ550mm是合理选择。接下来算装填高度。单塔容积0.25m3塔径Φ550mm可算出Ha 4V分 / (πD²) 4×0.25 / (π×0.55²) ≈ 1.05m高径比Ha:D≈1.91没有达到计算书要求的≥2。而原计算书写的是Ha1.58m高径比2.87反推其塔径约Φ450mm与前面圆整后的Φ550mm明显矛盾。这个矛盾在工程上很典型先按计算直径算出初步尺寸圆整标准直径后忘记重新迭代装填高度。正确做法是直接按标准直径Φ550mm反算装填高度再验证高径比若高径比不足1.91应把塔径改小一档到Φ500mm此时装填高度约1.27m高径比2.54同时空塔线速度0.0452m/s仅比额定值高一点配合更均匀的布料器也可以接受。所以实际设计这里要在两个方案里权衡Φ500mm满足高径比但线速度临界Φ550mm线速度余量大但高径比不足。很多制造厂最终会选Φ500mm并增加气流分布器。3.4 吸附塔内部结构尺寸与进出口管径吸附塔不是简单一个空罐子装满分子筛就完事。计算书把塔内纵向分成三个区域A区装3A分子筛用来吸附残余水分和二氧化碳保护主床层的5A分子筛B区装氧化铝干燥剂拦截水分并与5A分子筛分层装填C区是分子筛压紧附件通常用弹簧压板或气囊把床层压住防止气流冲击导致分子筛摩擦粉化。A区要在计算装填高度Ha基础上加0.05m修正然后依次加上h250mm、h3100mm左右的安装高度最终直边高H约1.81m。进出管径的确定逻辑是压缩空气入口管与空气净化组件后的主管路保持一致取d140mm因为PSA塔进口前没有其他节流装置氧气出口管先按产品气工况流量2.52m3/h计算理论内径约10.5mm考虑吸附塔升压时产品气回流和调阀需要一定余量取d212mm比空气进口小一档是合理的。这里注意氧气出口管径不能因为算出来细就随意缩小还要看塔内分布帽的开孔面积保证气流能均匀进入床层。3.5 用Python快速跑一台吸附塔的分子筛量分子筛装填量和塔径可以汇总成下面这个函数用于后续不同纯度、不同产量的方案对比def adsorber_calc(qn, qa_actual, K50, margin1.1, rho670): # qn 产品氧气流量 Nm3/hqa_actual 工况压缩空气流量 m3/h g_real qn * 1000 / K * margin # 总实际装填量 kg v_single g_real / 2 / rho # 单塔容积 m3 q_flow qa_actual / 3600.0 d_ads math.sqrt(4 * q_flow / (math.pi * 0.045)) * 1000 d_std math.ceil(d_ads / 50) * 50 # 标准塔径50mm一档 h_fill 4 * v_single / (math.pi * (d_std / 1000) ** 2) return {sieve_mass: g_real, sieve_volume: v_single, d_std: d_std, h_fill: h_fill} print(adsorber_calc(15, 31.93)) # 输出: {sieve_mass: 330.0, sieve_volume: 0.246, d_std: 550, h_fill: 1.05}这里g_real是总实际装填量v_single是单塔容积d_std按50mm一档向上取整得到标准塔径h_fill是用标准塔径反算的装填高度避免出现前面那种混用直径的情况。注意qn要传Nm3/hqa_actual要传工况实际流量两者单位不同不要混在一起。输出结果1.05m的装填高度已经暴露了高径比不足所以实际产品文档里的1.58m只能理解为制造图与计算书版本不一致复核时一定要以这个函数重算结果为准。4. 空气与氧气缓冲罐容积计算缓冲系数Ka怎么取才不掉链子4.1 为什么要两个缓冲罐Ka的物理含义吸附塔在切换瞬间进气阀和排气阀同时动作管线内流量会发生阶跃变化。空气缓冲罐的作用是把空压机连续输出的空气“缓存”起来在吸附塔充压时补充瞬时流量氧气缓冲罐则用来平滑产品气压力和纯度保证下游用户在切换期间不中断供气。两者容积都用缓冲系数Ka乘以对应流量求得但Ka的取值差别很大这不是随意定的而是由两侧的流量波动特性决定的。4.2 空气缓冲罐容积初算与结构反推空气缓冲罐初算容积V缓 Ka × Qa 0.011 × 190 2.09 m3这里的190是压缩空气耗气量单位Nm3/h不是氧气产量。Ka取0.011h相当于要给空压机供气侧约40秒的缓冲时间。有些设计人员想把空气缓冲罐做小一点把Ka降到0.008h结果吸附塔充压时进气压力掉压超过0.03MPa分子筛解吸不充分。因此Ka取值最好落在0.0100.015h区间并根据吸附塔容积和切换频率做微调。初取直边高1810mm标准椭圆封头估算直径D缓 sqrt(4×2.09 / (π×1.81)) ≈ 1.22m圆整为DN1200mm。然后反算直边高H缓 4×(V缓 - 2V封) / (πD²) 4×(2.09 - 2×0.225) / (π×1.2²) ≈ 1.45m圆整为1.5m。标准椭圆封头在DN1200下单个容积约0.225m3这部分数据来自GB/T 25198不能只用圆柱体积公式。反算完直边高后还要做最终校核V缓 π×1.2²×1.5/4 2×0.225 ≈ 2.15m3Ka实 2.15 / 190 ≈ 0.0113h与初取值基本吻合。整个流程的关键是先估直径-圆整-反算直边高-再校核Ka而不是直接用初算容积去画图。4.3 氧气缓冲罐的Ka为什么大一圈氧气缓冲罐初算容积V缓 Ka × Qn 0.031 × 15 0.465 m3Ka取0.0300.04h流量基准是产品氧15Nm3/h。计算书里有一行推导Ka(5.5/0.5)/(60×6)再乘绝压比(0.601/0.101)意思是吸附塔每小时切换6次每次切换产生约60秒的产气中断压力比越大需要的缓冲时间越长。虽然实际数字不一定精确对应但它点明了Ka的本质切换频率越高、吸附压力越高、产品气允许波动越小Ka就要取得越大。如果吸附塔采用3塔流程切换错开Ka可以适当取小2塔流程则必须取大值。按同样的结构算法初取直边高1810mm估算直径0.57m圆整为DN600mm。DN600椭圆封头容积约0.0352m3反算直边高约为1.4m最终总容积V缓 π×0.6²×1.4/4 2×0.0352 ≈ 0.466m3Ka实0.466/15≈0.031h符合预期。产品气流量越小罐体绝对容积越小但相对缓冲时间不能压缩这就是为什么15Nm3/h的小设备也要有一个直径600mm的氧气缓冲罐。4.4 两组缓冲罐参数对比与现场验证把两组计算放在一起看更清晰项目空气缓冲罐氧气缓冲罐缓冲系数Ka (h)0.0110.031流量基准压缩空气190 Nm3/h产品氧15 Nm3/h初算容积 (m3)2.090.465标准直径 DN (mm)1200600直边高 (mm)15001400实际总容积 (m3)2.150.466校核Ka (h)0.01130.031从表里能看到氧气罐的绝对容积只有空气罐的约1/4但它的Ka是空气罐的近3倍说明它承担的是高频小流量波动空气罐承担的是低频大流量波动。现场验证方法很简单在吸附塔切换的一瞬间看缓冲罐压力表空气罐压力下跌不超过0.03MPa、氧气罐压力波动不超过3%可以认为容积足够如果压力波动超限先检查切换阀开闭时间再考虑增加容积。还有一个容易犯的错Ka带时间量纲不要把它理解成“容积系数”。0.031×150.465m3这里的15是Nm3/h不是m3/h如果换成工况流量2.52m3/h去乘Ka就会变成0.184h罐体容积被放大约5倍。正确做法是设计和校核都统一用同一个流量基准并在计算书里写明取的是哪个流量。5. 把设计计算书写成可复用的Python小工具再避开四个工程坑5.1 一个覆盖主流程的粗算脚本前面每一步都可以单独手算但方案对比时会很慢。把主流程封装成一个函数输入氧气产量和纯度输出空气量、分子筛量、管径和塔径。下面是一个20行的版本def psa_design(qn, purity, eta0.352, p_gauge0.5, K50, rho_sieve670, omega0.045, velocity8.0): # qn 产品氧流量 Nm3/hpurity 氧纯度小数 qa qn * purity / (0.209 * eta) # 标态空气耗量 Nm3/h qa_act qa * 0.101 / (0.101 p_gauge) # 工况空气流量 m3/h d_pipe math.sqrt(4 * (qa_act / 3600) / (math.pi * velocity)) * 1000 g_ads qn * 1000 / K * 1.1 # 实际分子筛总装填量 kg v_sieve g_ads / 2 / rho_sieve # 单塔筛容积 m3 d_ads math.sqrt(4 * (qa_act / 3600) / (math.pi * omega)) * 1000 return { qa_标态: qa, qa_工况: qa_act, 主管内径mm: round(d_pipe), 分子筛装填量kg: round(g_ads), 单塔筛容积m3: round(v_sieve, 3), 吸附塔内径mm: math.ceil(d_ads / 50) * 50, } print(psa_design(15, 0.93)) # {qa_标态: 189.6, qa_工况: 31.9, 主管内径mm: 38, 分子筛装填量kg: 330, 单塔筛容积m3: 0.246, 吸附塔内径mm: 550}代码里的p_gauge必须传表压函数内部已经完成绝压修正qn是标况产品氧气流量purity用小数。返回的“主管内径mm”是理论值实际配管还要向上圆整到标准管径“吸附塔内径mm”已经做了50mm一档向上取整。当前参数下分子筛量330kg、塔径550mm和计算书主线一致。5.2 四个工程坑第一分子筛堆密度单位。计算书写的670kg/L是笔误真实单位是kg/m3任何质量除以密度算容积时都要先确认单位。第二吸附塔直径圆整后必须重新反算装填高度。原计算书Ha1.58m与Φ550mm塔径不匹配用标准直径反算出的1.05m才符合几何关系高径比不足时优先改塔径而不是改Ha。第三缓冲系数的流量基准不能混。空气罐乘压缩空气量氧气罐乘产品氧量两者差一个回收率折算。换用工况流量时Ka的数值必须同步重新设定不能只改一个乘数。第四K和η是工况实测值。0.5MPa、20℃下测得的50NL/kg和35.2%换到0.7MPa或40℃环境就不成立。方案对比可以用经验值正式设计一定要向分子筛厂家索取该温度压力下的吸附曲线并按吸附塔实际循环周期重新校核产氧率。遇到旧计算书里的矛盾不要急着否定整个方案。把几何尺寸独立复核一遍再拿现场半小时的压力和纯度记录对一下通常能判断是版本迭代还是抄写错误。对气体设备这类项目建议把几何尺寸、物料质量、阀门切换逻辑分开归档任何一次设计变更都重新跑一遍上面的粗算脚本再改图纸。本文还有配套的精品资源点击获取