氢能项目有个很有意思的现象立项评审时大家都在聊电解槽、燃料电池和可再生电力真到工艺计算阶段最难啃的往往是被忽略的储运环节。常压氢气体积能量密度很低即便70 MPa的Ⅳ型瓶整车的储氢质量分数也就4%~6%而液氢的体积密度能做到70 kg/m³以上同样尺寸的储罐和槽车液氢的运输容量高出压缩氢一个量级。于是氢液化循环就成了氢能产业链上绕不开的技术节点。Aspen Plus是我做这类工艺设计的主工具这篇文章把氢液化循环从流程选型、物性方法、模型搭建、能耗核算到常见调试坑完整走一遍正在做工艺包、写课题或者准备预可研的朋友可以直接拿来当参考。1. 先算清楚账再动手液氢在储运体系里的定位与两个物理门槛1.1 压缩氢、液氢、低温压缩氢到底差多少给一个直观对比表储运方式体积储氢密度(kg/m³)系统质量储氢密度典型能耗(占氢气LHV比例)主要应用场景350~700 bar压缩氢23~404%~7%约8%~12%加氢站、重卡、叉车液氢70左右5%~9%(含储罐)30%~45%航天、长途海运、大规模储运低温压缩氢40~805%~9%介于两者之间特殊场景、研究项目金属氢化物/LOHC20~701%~5%放氢需额外热量小众场景液氢的最大优势就是体积密度。同样是20吨运力的液氢槽车如果改用700 bar压缩氢车体自重和设备成本会显著增加单位运输成本立刻上去。这也是为什么长距离跨国运输和航天发射这类要求高密度存储的场景液氢几乎是唯一解。液氢的代价在能耗和低温。氢气要从常温一路冷却到沸点20 K附近再冷凝这个过程理论最小功耗大约在3.3 kWh/kg上下实际工业装置普遍在10~15 kWh/kg。按氢气低热值33.3 kWh/kg算液化就要吃掉氢自身能量的30%~45%。这个数字很扎眼所以工艺设计的目标一直很明确把液化能耗压下来压不下来液氢就只能在少数刚需场景里用。1.2 20 K低温与正仲氢转化液化绕不开的两道坎氢液化难不只是因为温度低。第一个坎是温区跨度大从常温到20 K冷量品位完全不同单靠一台制冷机没法高效覆盖第二个坎是正仲氢转化。常温下氢是75%正氢、25%仲氢的混合物而20 K液氢温度下的平衡态仲氢含量超过99%。正氢向仲氢转化是放热过程转化热约527 kJ/kg——这个数比液氢本身的汽化潜热446 kJ/kg还大。如果液化出来的是正仲比例不平衡的“正常氢”储罐里的氢会一边自发转化一边放热导致液氢持续蒸发几天内就能白白损失两三成的液体。所以工业液化装置必须在流程里加装催化剂转化器通常在80 K、40 K、20 K左右分几段把正氢转化为仲氢把转化热随冷量一起排掉。这两道坎决定了氢液化流程绝不是简单把气体压缩冷却就能完成的它需要一个多温位、多设备的深冷系统。而这种系统的物料平衡、能量平衡和冷量匹配靠手算几乎不可能Aspen Plus的价值就在这里把循环画出来把物性算准就能在电脑上先跑几十个方案再谈建装置。2. 循环是怎么一步步变复杂的反转温度、膨胀机与预冷策略2.1 焦汤节流的反转温度决定了氢不能常温直接液化先讲最底层的原理。气体节流膨胀焦汤效应能不能降温取决于它的焦汤系数。大多数气体在常温下节流会降温但氢有个特殊的反转温度约202 K。高于这个温度氢节流反而升温只有先把氢预冷到202 K以下节流膨胀才有制冷效果。这意味着最经典的Linde-Hampson循环——气体压缩、冷却、节流、分离——如果直接用常温氢做不仅液化不了节流后温度还会上升。必须先用外部冷源把氢从常温预冷到80 K左右再进节流阀才能得到液氢。Linde-Hampson循环结构最简单全部冷量靠焦汤效应产生节流过程本身不可逆损失极大液化单位产品的能耗非常高。纯理论算下来这种循环用于氢液化时单位能耗可以到40 kWh/kg以上基本没有工程吸引力它更大的意义是教学和验证性质。2.2 Claude循环的核心用膨胀机把冷量“挣”出来Claude循环在节流阀前面并联了一台透平膨胀机。高压氢气分流一部分进膨胀机气体对外做功出膨胀机后温度大幅下降这股低温气体进入换热器为主流氢提供冷量另一部分高压氢继续冷却后进节流阀液化。膨胀机的引入带来两个好处。第一冷量直接从工质的内能“买”出来而不是全靠节流的焦汤效应效率高得多第二膨胀机输出轴功可以回收抵消一部分压缩机功耗。一个设计合理的Claude循环液化单位能耗能降到12~20 kWh/kg比Linde-Hampson低一个档次。再把Claude循环的膨胀级数增加变成两级、三级膨胀就是Collins循环每增加一级膨胀冷量供给的温位匹配更细不可逆损失更小能耗进一步下降。2.3 液氮预冷和混合制冷剂预冷什么时候该上Claude循环通常在常温到80 K这段再加一道预冷。最常用的是液氮预冷液氮常压在77 K左右相变能把高压氢从常温直接冷却到90 K以下把压缩机功耗和膨胀机负荷都降下来。液氮来源广泛装置简单中小型氢液化装置很喜欢用。缺点是需要液氮供应链且液氮预冷只覆盖高温段80 K以下还是要靠膨胀机和节流。规模更大的装置比如LNG冷能可利用或氢气源就在大型化工园区的场合也可以用混合制冷剂预冷。混合制冷剂循环用多组分工质在不同温位相变效率比液氮高省掉了液氮的采购和储存成本但流程复杂、组分配比调起来很费时间。用Aspen Plus做筛选时我一般按下面这张表分类比较循环类型预冷方式膨胀配置典型SEC(kWh/kg)复杂程度适用规模Linde-Hampson必须有外部预冷仅节流30以上低实验室、原理演示简单Claude可不预冷1级透平节流15~25中中小型预冷Claude液氮或混合制冷剂1级透平节流8~12中高现代工业装置Collins液氮或混合制冷剂2~3级透平节流6~9高大型集中液化表里的SEC是粗略范围具体数值会受规模、效率和冷源条件影响但它能帮你快速判断如果只要5 t/d的小装置上Collins纯属给自己找麻烦如果是50 t/d以上的大型基地简单Claude预冷又不够看。3. 物性方法挑不好20 K附近全是坑3.1 Aspen Plus常用状态方程在深冷段的真实表现物性方法是流程模拟的根。氢的临界点很低临界温度33.2 K临界压力1.3 MPa在常温高压和中低温区间氢的气相行为其实很接近理想气体但一旦进入20 K附近的液态区分子间作用力、量子效应和正仲氢状态都会让常规状态方程表现变差。Aspen Plus里最常见的几种物性方法在氢液化场景下的表现需要甄别物性方法优点问题我的用法PR Peng-Robinson参数全、收敛稳、气相计算准液氢密度和焓值在低温段偏差偏大流程快速筛选、预收敛PR-BM修正了Boston-Mathias alpha函数外推行为好一点深冷液相区仍不够精确基算例首选SRK与PR类似同样的低温液相偏差问题对比校验用LK-PLOCK面向低温气体处理的Lee-Kesler类方法组分参数覆盖有限设计阶段主力GERG-2008参考级精度尤其混合气Aspen Plus常规版本不直接提供需要外部集成最终校验基准为什么PR在20 K附近会有偏差PR的alpha函数是根据临界点和偏心因子拟合的氢的临界参数很低低于临界温度很多时PR对饱和液相密度的预测误差可以到百分之几甚至更大这个误差放在冷量平衡里换热器和膨胀机功率都会跟着偏。我做过一次对比用PR和用REFPROP查表数据算同一个80 K到20 K换热的焓差相差接近5%对功耗核算来说这个量级已经不能忽略了。3.2 我的物性校验套路PR-BM筛流程LK-PLOCK定设计REFPROP查底我的做法是分阶段用物性方法。方案比选阶段几十个流程要跑稳定性和速度优先用PR-BM跑趋势完全没问题到了有希望的候选方案换成LK-PLOCK重新算一遍SEC和换热器负荷都会变一点最后把关键节点的温度、压力、焓值、熵值导出和NIST REFPROP的氢物性表或者文献里的饱和性质数据比对确认偏差在可接受范围。有个细节容易忽略Aspen Plus组件库里的H2是把正氢和仲氢当作同一组分处理的热力学性质是基于平衡混合物或正常氢拟合的。流程模拟本身不受太大影响但涉及正仲氢转化热的计算Aspen默认算不出来需要你在转化段单独加一个热负荷来表示转化放热。后面第7节会讲具体怎么补。提示物性方法不是一次定死就不动的。同一个流程方案比选、初步设计、详细设计这三个阶段用的物性方法完全可以不同关键是要知道每个阶段在做什么精度的决策。4. 预冷Claude循环的建模过程从空白文件到收敛出数4.1 流程图先搭骨架模块清单与连接关系我用Aspen Plus搭预冷Claude循环模块清单基本固定原料氢温度300 K、压力1.1 bar近似常压流量按目标产能定比如先设100 kg/h。压缩机群用MCompr或三个Compr串联中间加水冷器目标压力20~30 bar。氢的绝热指数约1.41单级直接从1.1 bar压到25 bar理想出口温度要奔着700 K去材料、润滑油和安全性都撑不住所以必须多级压缩带级间冷却一般按每级压比2.5~2.8配三级。预冷换热器HX1高压氢被液氮冷却到80~90 K。液氮侧可以简化为77 K饱和液氮进、气态氮出。分流器预冷后的高压氢分两股一股去膨胀机一股去主换热器。膨胀机EXP用Compr块选Isentropic等熵效率0.82~0.88出口压力4~5 bar。主换热器HX2剩余高压氢和返流的低压低温氢换热高压氢被冷却到40~50 K。节流阀高压氢从20~30 bar节流到1.3 bar等焓膨胀。气液分离器Flash2节流后的两相流体在这里分离底部出液氢产品顶部低温氢气返流经过HX2、HX1复热然后循环回压缩机入口或放空。模块连接关系就是一张“高压回路返流回路”的循环图压缩机出口→HX1→分流→膨胀机支路和HX2→节流阀→分离器分离器气相→HX2热侧→HX1热侧→压缩机入口。膨胀机出口也接回HX2热侧与返流气汇合。4.2 参数设置里的门道效率、夹点与压力比几个关键参数给出来直接抄压缩机等熵效率0.78~0.85。氢压缩机的实测效率不会太高别用0.9去乐观估计否则算出来的功耗和实际电费对不上。级间冷却温度300 K也就是冷却到室温水冷或空冷都按这个温度设定。每级压力比2.5~2.8。20 bar目标时三级压缩每级约2.7左右。透平膨胀机等熵效率0.82~0.88出口压力设置成4~5 bar。出口压力不能太低否则膨胀机出口温度过低后面返流气体与主流的温度匹配会出问题。换热器最小换热温差冷端3 K热端5 K。氢液化系统温差压得很低因为温差就是冷量损失但Aspen Plus里换热器如果设温差过小收敛会非常痛苦先按3~5 K算后面再逐步收紧。节流阀用Valve块只设压降。这就是等焓节流的正确模拟方式千万别在阀上设热负荷。分离器压力1.3 bar左右。实际装置不会让分离器处于常压稍微带点正压有利于克服下游阻力也能减少空气渗入风险。进膨胀机的分流比例和液氮预冷量需要配合设计这个用设计规格来做见下一节。4.3 收敛顺序和设计规格让流程按你要的结果走循环流程最大的问题是返流换热器构成的回路。Aspen Plus会自己选撕裂流但初值不好就会有“流程拓扑没问题、数值震荡不收敛”的尴尬。我的收敛顺序是先用较宽松的换热器出口温度赋初值或者把回路断开试一次收敛后再把回路接上然后用设计规格逐步收紧目标。两条设计规格必设改变量膨胀机分流比目标节流阀进口温度40~45 K。这个温度直接决定节流后液体产率温度越低液氢产率越高但不能低到让主换热器出现温度交叉。改变量液氮用量目标HX1出口氢温度80~90 K。液氮用量直接决定预冷负荷也是能耗大头。4.4 一个示意基算例的流股温度压力表下边给一个100 kg/h原料气、采用液氮预冷Claude循环的示意算例流股参数供对量级用具体数值会因物性方法和效率设定略有差别位置温度(K)压力(bar)说明原料氢3001.1纯氢常压供应压缩机出口30020三级压缩加级间冷却后HX1氢出口8519.5液氮预冷后带0.5 bar压降膨胀机入口8519.5分流后去膨胀机膨胀机出口40~504.5等熵效率0.85主换热器出口(节流前)4219被返流冷却到设计值节流后20.51.3两相进分离器液氢产品20.41.3饱和液相返流氢气20.5→接近常温1.3从分离器气相复热液氢产量取决于节流前温度、压力和返流气温度一般100 kg/h原料氢的流程液氢产量大概在60~90 kg/h剩下的返流气循环再液化。这个“一次液化率”不追求高让返流气循环就行总能耗才是关键指标。5. 能耗怎么核算SEC、可用能与Aspen里的计算器块5.1 理想功、实际功耗与效率指标氢液化的理论最小功不是拍脑袋猜的它由始末状态决定从环境态比如300 K、1 bar的氢气到目标态20 K、1.3 bar的饱和液氢最小功等于焓变化减去环境温度乘熵变化也就是 W_min ΔH - T0·ΔS。按常用的环境参考态算氢液化的理想功耗大致3.3 kWh/kg。这个数很有用它是所有对比的基准。实际装置的单位液化能耗SEC 总输入电功率 / 液氢产量工业范围10~15 kWh/kg。把理想功除以实际SEC得到无量纲的效率分数FOM通常在0.25~0.4之间。换句话说现有装置比理想过程多花了2.5到4倍的能量这就是工艺优化要啃的空间。5.2 损失分布告诉你下一步该优化谁能耗高能量去哪了用可用能分析exergy analysis也叫㶲分析能看到每一项设备的不可逆损失。对预冷Claude循环做一个典型损失分布设备/环节典型可用能损失占比原因压缩机与级间冷却30%~35%压缩功转化为热冷却器排走换热器20%~25%有限温差、压降节流阀15%~20%等焓膨胀不可逆膨胀机5%~10%等熵效率不足、入口温差损失液氮预冷环节10%~15%液氮生产与传热温差这个分布说明三件事压缩机是最大的能耗黑洞但它不可绕过换热器的有限温差是被低估的优化点节流阀的损失可以通过提高节流前压力和温度匹配来压缩。如果你把节流阀替换成液力透平湿膨胀机那15%~20%的损失能回收一部分这就是现在液氢装置都在研究“液体膨胀机替代节流阀”的原因。在Aspen Plus里做可用能分析有两种路径新版软件自带Exergy分析工具或者自己在计算器块里按 e (h - h0) - T0(s - s0) 对每个流股求可用能值然后用Mixer和加法器汇总出入平衡。两者都能用后者对流程变量的控制更灵活。5.3 灵敏度分析改哪个参数最见效模型收敛后马上跑灵敏度分析。我最常扫的三个变量膨胀机分流比、压缩机出口压力、液氮用量。趋势基本是分流比增大冷量多了节流前温度降低、液氢产率上升但膨胀机处理量过大时换热器内冷热流匹配变差返流温度过低SEC先降后升存在最优区间。压缩机出口压力升高焦汤效应增强但压缩功耗上升更快通常也存在一个最优压力点。液氮用量加大预冷温度更低主换热器和膨胀机负荷减小但液氮本身的制备能耗要计入总账这个最优值藏得比较深最容易漏算。灵敏度分析结果直接用Aspen Plus的Sensitivity模块跑出来一条U型曲线找到谷底就是该工况下的较优设计点。6. 调试实录从“不收敛”到“结果可信”的完整链路6.1 收敛问题的排查顺序氢液化流程第一次跑不收敛概率最高的几个原因按排查顺序列一下看错误报告里的物料和能量失衡标志。氢的分子量小很小的流股误差都会在能量平衡里放大。换热器温度交叉。Aspen Plus报“Temperature cross”十有八九是换热器出口温度给得不符合热力学比如冷端温差设成了负数。先放开出口温度让换热器自由算再把夹点设置加回去。撕裂流初值太离谱。返流气温度给300 K去换20 K的冷流第一次迭代就会爆。手动把撕裂流的初值改成大致合理的温度压力比如返流气温度先给80 K再让流程迭代。把阻尼因子调小。Aspen的默认收敛参数对强非线性深冷系统经常太激进把TOLERANCE适当放大阻尼设0.5~0.7先求“能收敛”再求“收敛准”。6.2 换热器夹点与多股流冷箱的处理换热器是氢液化模拟里最大的麻烦制造者。单股流用HeatX就够了返流气体超过两股的地方必须用MHeatX。MHeatX收敛有一个经典技巧把每个流股的出口温度指定成设计值Design模式先算通等全流程收敛后再切回Rating模式验算面积和UA。夹点分析要单独做。在Aspen里把换热器冷热复合曲线导出来看最小温差出现在哪个温位。氢液化系统最小温差通常出现在低温端也就是20 K附近因为那个温位的比热和潜热变化很剧烈最容易夹死。如果低温端接近零温差说明高压流和低压返流的流量比没调好优先检查膨胀机分流比而不是硬去调换热器面积。6.3 物性之外的软件坑字体、安装包与数据来源最后聊两个和流程模拟无关、但几乎每个Aspen用户都会撞上的事。一个是SHX字体问题。Aspen Plus在打开某些带绘图或报告模板的文件时会提示缺少SHX字体文件界面图形区和打印输出直接乱码或者报错。解决办法就是从可信渠道把对应的SHX字体文件放到Aspen安装目录的Fonts文件夹里。这个坑跟模拟水平没关系纯属环境配置但卡住很多人。网上搜“aspen plus cad shx字体下载”能找到一些帖子注意甄别来源别随便下载来路不明的压缩包。另一个是软件获取方式。网上找个网盘链接下载所谓“绿色版”“破解版”Aspen Plus是很危险的做法——且不说法律风险这类安装包经常被植入后门工艺数据泄露了损失无法估量。AspenTech官方提供试用授权和学生版本起步完全够用项目上需要长期用该买授权就买授权。别为了省一点软件成本把工艺包的核心数据置于风险里。7. 模型进阶与工程落地从基算例走向可用的工艺包7.1 把正仲氢转化补进模型能耗更接近现实前面提过Aspen默认不区分正仲氢。但工程上这个转化热会影响能耗和储运安全。实际操作是在换热流程的温位节点上加一个负的热负荷或者用很小的Heater块设置负Duty来代表转化放热这个热负荷由冷箱内的冷流带走。段落设置一般两级到三级80 K一级大部分转化热在这里释放、40 K一级、20 K一级。转化器的装填量可以按转化动力学估算稳态模拟里先按平衡转化率算总放热量再分配到各温位。加上转化热之后SEC会比纯Aspen默认结果高5%~8%这是个很真实的修正写可研报告时建议保留。7.2 装置规模、设备选型与跨软件校验模拟收敛只是第一步。下一步把关键设备参数从Aspen导出来压缩机功率和出口温度、冷箱换热面积和UA值、膨胀机功率和转速量级、分离器尺寸。冷箱面积的估算尤其重要20 K温区的换热器通常用铝制板翅式或绕管式单位面积造价极高面积差10%就是几百万元的成本差。我习惯还用Aspen HYSYS把同一流程再搭一遍对比关键节点的流量和温度。HYSYS在低温气体处理和LNG换热器方面积累的经验参数更多两套软件算出来的SEC如果能对到5%以内这个模型才敢拿出去给别人评审。注意两套软件的物性方法必须设成一致否则对比的是物性差异而不是流程差异。7.3 再往前一步混合制冷剂、氦循环与动态模拟模型框架搭好以后扩展方向就很清晰了。想摆脱液氮供应链就把预冷段改成混合制冷剂循环在Aspen里用换热器复合曲线优化MR组分配比想挑战更低的单位能耗可以在膨胀机后面再接氦气布雷顿循环作为低温冷源做装置设计和安全分析时把稳态模型导入HYSYS Dynamics或者Unisim做动态泄压、储罐滚动周期和BOG处理模拟能提前发现很多稳态看不到的问题。我个人经验是不要把流程模拟当成终点。Aspen模型最大的作用是帮你在纸面上把方案筛掉80%剩下20%值得花时间和钱去做小试、中试和详细设计。每一次优化都要回到能耗、投资和可靠性三个维度去权衡——氢液化本身已经够贵了不要为了追求账面能耗降低1 kWh/kg就把设备投资和维护成本拉爆。