氢能项目有个很有意思的现象:立项评审时大家都在聊电解槽、燃料电池和可再生电力,真到工艺计算阶段,最难啃的往往是被忽略的储运环节。常压氢气体积能量密度很低,即便70 MPa的Ⅳ型瓶,整车的储氢质量分数也就4%~6%;而液氢的体积密度能做到70 kg/m³以上,同样尺寸的储罐和槽车,液氢的运输容量高出压缩氢一个量级。于是氢液化循环就成了氢能产业链上绕不开的技术节点。Aspen Plus是我做这类工艺设计的主工具,这篇文章把氢液化循环从流程选型、物性方法、模型搭建、能耗核算到常见调试坑完整走一遍,正在做工艺包、写课题或者准备预可研的朋友可以直接拿来当参考。
1. 先算清楚账再动手:液氢在储运体系里的定位与两个物理门槛
1.1 压缩氢、液氢、低温压缩氢到底差多少
给一个直观对比表:
| 储运方式 | 体积储氢密度(kg/m³) | 系统质量储氢密度 | 典型能耗(占氢气LHV比例) | 主要应用场景 |
|---|---|---|---|---|
| 350~700 bar压缩氢 | 23~40 | 4%~7% | 约8%~12% | 加氢站、重卡、叉车 |
| 液氢 | 70左右 | 5%~9%(含储罐) | 30%~45% | 航天、长途海运、大规模储运 |
| 低温压缩氢 | 40~80 | 5%~9% | 介于两者之间 | 特殊场景、研究项目 |
| 金属氢化物/LOHC | 20~70 | 1%~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) | 说明 |
|---|---|---|---|
| 原料氢 | 300 | 1.1 | 纯氢,常压供应 |
| 压缩机出口 | 300 | 20 | 三级压缩加级间冷却后 |
| HX1氢出口 | 85 | 19.5 | 液氮预冷后,带0.5 bar压降 |
| 膨胀机入口 | 85 | 19.5 | 分流后去膨胀机 |
| 膨胀机出口 | 40~50 | 4.5 | 等熵效率0.85 |
| 主换热器出口(节流前) | 42 | 19 | 被返流冷却到设计值 |
| 节流后 | 20.5 | 1.3 | 两相,进分离器 |
| 液氢产品 | 20.4 | 1.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,就把设备投资和维护成本拉爆。