储氢容器日蒸发率与真空多层绝热性能衰减计算书

储氢容器日蒸发率与真空多层绝热性能衰减计算书第1页储氢容器日蒸发率与真空多层绝热性能衰减计算书（工程计算与安全核验版）文件编号：H2-TC-2026-001设计阶段：详细工程设计文件类型：绝热性能计算书编制：___________日期：___________校核：___________日期：___________审核：___________日期：___________标准化审查：___________日期：___________批准：___________日期：___________2026年10月

目录TOC\o"1-2"\h\z\u1概述 11.1计算目的 11.2适用范围 11.3引用标准 12设计输入参数 22.1储罐几何参数 22.2液氢介质物性参数 22.3环境设计条件 33高真空多层绝热漏热计算 33.1漏热机理分析 33.2辐射换热计算 43.3固体导热计算 43.4残余气体导热计算 53.5支撑结构与管路漏热 54静态日蒸发率计算 64.1基本计算公式 64.2额定工况日蒸发率计算 64.3不同容积规格日蒸发率对比 74.4蒸发率温度敏感性分析 75真空多层绝热性能衰减计算 85.1真空度衰减机理 85.2材料放气率计算 85.3吸附剂吸气能力与寿命 95.4真空度与漏热的定量关系 95.5服役年限-性能衰减预测 106安全核验与裕度分析 106.1标准符合性核验 106.2真空失效安全分析 116.3运输工况动态蒸发率修正 117结果分析与工程指导意义 127.1主要计算结论 127.2工程设计优化建议 127.3在役检测与维护建议 138免责声明 13第1页1概述1.1计算目的本计算书针对采用高真空多层绝热（HighVacuumMulti-LayerInsulation,HV-MLI）技术的低温液态储氢容器，进行系统的日蒸发率（Boil-OffRate,BOR）理论计算，并建立真空多层绝热性能随时间衰减的预测模型。通过多维度漏热分析、真空度衰减计算、吸附剂寿命评估及安全裕度核验，为储氢容器的工程设计、制造验收、在役检测及安全运维提供完整的技术依据。计算结果将用于验证储氢容器在设计寿命周期内的绝热性能是否满足GB/T40060-2021《液氢贮存和运输技术要求》中关于日蒸发率≤0.3%/d的强制要求，并给出不同服役年限下的绝热性能演化规律及真空维护周期建议。1.2适用范围本计算书适用于公称容积5~200m³、设计压力0.2~1.6MPa、工作温度20.3K（液氢沸点）至313K（环境温度）范围内的立式/卧式固定式液氢储氢容器，以及40ft液氢罐式集装箱。计算模型基于稳态传热理论与真空放气动力学理论，对采用镀铝聚酯薄膜/玻璃纤维纸交替组合的高真空多层绝热结构进行计算。1.3引用标准表1-1引用标准一览表标准编号标准名称GB/T40060-2021液氢贮存和运输技术要求GB/T18443.5-2010真空绝热深冷设备性能试验方法第5部分：静态蒸发率测量GB/T18443.2-2010真空绝热深冷设备性能试验方法第2部分：真空度测量GB/T18442.3-2011固定式真空绝热深冷压力容器第3部分：设计GB/T34542.1-2017氢气储存输送系统第1部分：通用要求GB/T29729-2022氢系统安全的基本要求GB150.3-2011压力容器第3部分：设计NB/T10356-2019真空绝热深冷压力容器高真空多层绝热技术条件2设计输入参数2.1储罐几何参数本计算以某典型100m³立式液氢储氢容器为基准算例进行详细计算，其他容积规格可按相同方法进行参数化计算。储罐采用同心双壁圆筒结构，主要几何参数如下表所示。表2-1储罐几何参数表序号参数名称符号单位数值备注1公称几何容积Vgm³100设计任务书2内筒内径Dim2.80优化设计3内筒筒体直段长度Lim15.00长径比约5.4:14内筒壁厚δimm1206Cr19Ni10，强度计算确定5椭圆封头曲面深度hbm0.70标准椭圆封头a/b=26绝热层厚度δinsm0.30高真空多层绝热夹层7外筒内径Dom3.44Do=Di+2δi+2δins8多层绝热总层数N层60镀铝箔+玻璃纤维纸交替9设计充装率φ—0.85GB/T40060规定≤90%10有效换热面积（含封头）Am²148.5计算值2.2液氢介质物性参数液氢在标准大气压（0.101325MPa）下的饱和温度为20.3K，对应物性参数依据美国国家标准与技术研究院（NIST）热力学数据库及《低温工程热物性手册》选取，具体参数如下表所示。表2-2液氢（仲氢）热物性参数表（20.3K,0.101325MPa）序号参数名称符号单位数值备注1饱和温度TLH2K20.27标准沸点2饱和液体密度ρlkg/m³70.79液态仲氢密度3饱和气体密度ρgkg/m³1.338气态仲氢密度4气化潜热rkJ/kg446.0标准沸点下汽化潜热5液态定压比热容cplkJ/(kg·K)9.67饱和液态6气态定压比热容cpgkJ/(kg·K)12.15饱和气态7液体动力粘度μlμPa·s13.2—8体积膨胀系数βlK⁻¹0.016—2.3环境设计条件表2-3环境设计条件参数表序号参数名称符号单位数值备注1环境温度（夏季极端）TaK313.1540℃，极端工况2环境温度（年平均）Ta,avgK288.1515℃，年平均3环境温度（标准测试）T0K293.1520℃，GB/T18443基准4大气压力PakPa101.325标准大气压5外壳表面发射率εo—0.85碳钢外表面涂漆6设计寿命Y年15固定式储氢容器3高真空多层绝热漏热计算3.1漏热机理分析高真空多层绝热（MLI）是目前液氢温区性能最优的绝热形式，其总漏热量由三部分组成：辐射换热、固体导热（通过间隔材料）和残余气体导热。当夹层真空度优于1×10⁻²Pa时，残余气体导热被有效抑制，总漏热主要由辐射换热和固体导热构成。根据Lockheed逐层计算模型，MLI的表观导热系数满足叠加原理：λeff=λr+λs+λg(3-1)式中：λeff为MLI表观导热系数，W/(m·K)；λr为辐射换热等效导热系数，W/(m·K)；λs为固体导热等效导热系数，W/(m·K)；λg为残余气体导热系数，W/(m·K)。3.2辐射换热计算多层绝热中相邻两层反射屏之间的辐射换热采用Stefan-Boltzmann定律计算。对于N层反射屏，冷壁温度Tc=20.3K，热壁温度Th=313K，考虑反射屏的有效发射率εeff，单位面积辐射热流密度计算公式如下：qr=σ·(Th⁴−Tc⁴)/[(N+1)·(2/ε−1)](3-2)式中：σ为Stefan-Boltzmann常数，σ=5.67×10⁻⁸W/(m²·K⁴)；ε为反射屏表面发射率，双面镀铝聚酯薄膜ε=0.04；N为绝热层总层数。代入参数计算：Th=313K，Tc=20.3K，N=60层，ε=0.04。由于Tc⁴=20.3⁴≈1.7×10⁵K⁴，远小于Th⁴=313⁴≈9.59×10⁹K⁴，冷壁辐射项可近似忽略。计算得：qr=5.67×10⁻⁸×9.59×10⁹/[61×(2/0.04−1)]=542.8/2989≈0.182W/m²。辐射等效导热系数：λr=qr·δins/(Th−Tc)=0.182×0.30/292.85≈1.86×10⁻⁴W/(m·K)。3.3固体导热计算多层绝热中，固体导热主要通过间隔材料（玻璃纤维纸）以及反射屏之间的接触点进行。层密度n（层/cm）是影响固体导热的关键参数——层数过少辐射抑制不足，层数过多层间接触压力增大导致固体导热升高，存在最优层密度使总漏热最小。对于液氢温区地面储罐，工程最优层密度约为20~25层/cm。本设计采用n=20层/cm（即300mm绝热层60层）。固体导热热流密度采用经验公式：qs=Cs·f·ks·(Th+Tc)/2·(Th−Tc)/δ(3-3)式中：Cs为经验系数，取0.008；f为层间接触系数，取0.35；ks为间隔材料导热系数，玻璃纤维纸在平均温度下ks≈0.012W/(m·K)；δ为单层间距，δ=δins/N=0.005m。代入计算：qs=0.008×0.35×0.012×166.65×292.85/0.005≈0.328W/m²。固体导热等效导热系数：λs=qs·δins/(Th−Tc)=0.328×0.30/292.85≈3.36×10⁻⁴W/(m·K)。3.4残余气体导热计算夹层残余气体导热与真空度密切相关。当气体分子平均自由程λmfp大于夹层间距d时（Knudsen数Kn＞1，即自由分子流区），气体导热与压力成正比；当真空度较差进入连续介质区时，气体导热急剧增大。残余气体导热计算公式为：qg=α0·P·(γ+1)/(γ−1)·√(R/(2πMT))·(Th−Tc)(3-4)式中：α0为热适应系数，氢气对铝/玻璃表面α0≈0.4；P为夹层残余气体压力，Pa；γ为气体比热比，氢气γ=1.4；R为通用气体常数，R=8.314J/(mol·K)；M为气体摩尔质量，氢气M=0.002kg/mol；T为气体平均温度，K。设计封结真空度P≤5×10⁻⁴Pa（冷态），代入计算得：qg≈0.007W/m²，残余气体导热在总漏热中占比＜1.5%，可忽略。真空恶化时的气体导热贡献见第5章性能衰减计算。3.5支撑结构与管路漏热除绝热层本体漏热外，支撑结构（8根G10环氧玻璃钢吊杆/支撑管）和工艺管路（充装管、排气管、仪表管）构成的固体热桥也是重要漏热路径。采用加长热桥长度、减小截面积、选用低导热材料的设计方案。各路径漏热计算汇总如下表所示。表3-1各漏热路径漏热量汇总表序号漏热路径热流密度(W/m²)漏热量(W)占比(%)1MLI绝热层（辐射+固体+气体）0.51776.849.22内筒支撑结构（G10玻璃钢）—46.429.73工艺管路系统—25.016.04测量仪表引线—3.62.35接管颈与法兰热桥—4.32.8合计总漏热量Qtotal—156.1100.0总漏热面积A=148.5m²，最终总漏热量Qtotal≈156.1W，作为日蒸发率计算的输入。需要说明的是，支撑和管路漏热属于集中热桥，按总换热面积平均的热流密度不适用，故在表中单独列项计算。4静态日蒸发率计算4.1基本计算公式静态日蒸发率（BOR）定义为储氢容器在额定充满率、热平衡稳定状态下，每24小时内因外界热量侵入导致液氢蒸发损失的质量百分比。根据能量守恒定律，侵入总热量全部用于液氢的汽化：BOR=(Qtotal×86400)/(ρl·Veff·r)×100%(4-1)式中：Qtotal为总漏热量，W（即J/s）；86400为一天的秒数，s/d；ρl为液氢密度，kg/m³；Veff为有效充装容积，m³，Veff=Vg×φ=100×0.85=85m³；r为液氢气化潜热，J/kg。4.2额定工况日蒸发率计算代入额定工况参数（40℃极端高温）计算：（1）每日总侵入热量：Qday=156.1×86400=1.349×10⁷J/d=13487kJ/d；（2）有效装载液氢质量：m_LH2=ρl×Veff=70.79×85=6017.2kg；（3）总汽化潜热容量：H_total=m_LH2×r=6017.2×446.0=2,683,671kJ；（4）质量日蒸发率（40℃）：αm=Qday/H_total×100%≈0.503%/d。上述计算为环境温度40℃极端工况下的被动储存蒸发率。按GB/T18443.5标准，换算至标准环境温度20℃（293.15K）及年平均温度15℃，并考虑配置主动制冷机制冷量抵消部分漏热，修正结果如下。对于多层绝热容器，温度修正公式为：α0=α'·(h1/h0)·[(293.15⁴−T2⁴)/(293.15⁴−T1⁴)]·(0.7T0/T2+0.3T1/T2)(4-2)年平均环境温度15℃（288K）工况下，无制冷机被动储存日蒸发率BOR≈0.42%/d；若配置1台50W@20KGM制冷机，可抵消约32%的漏热，有效日蒸发率降至约0.28%/d，满足GB/T29729-2022规定的≤0.3%/d安全限值要求。大型固定式液氢球罐通常配置主动制冷以实现近零蒸发储存（ZBO）。4.3不同容积规格日蒸发率对比储氢容器的容积效应（体积/表面积比）对日蒸发率有显著影响。容积越大，单位容积对应的漏热面积越小，日蒸发率越低。不同容积规格的日蒸发率计算结果如下表所示（均为40℃极端工况、无制冷机被动储存）。表4-1不同容积液氢储罐日蒸发率对比公称容积(m³)内径(m)换热面积(m²)总漏热(W)极端BOR(%/d)年均BOR(%/d)备注52.226.526.51.871.34无制冷102.633.633.61.170.84无制冷40ft罐箱≈4268.568.50.840.60无制冷1002.8148.5156.10.5030.36可选配制冷机2003.2238.0235.20.3580.26建议配制冷机1000(球罐)12.4554.0446.80.0960.069配制冷机近零蒸发2000(球罐)15.6880.0662.00.0420.030配制冷机近零蒸发4.4蒸发率温度敏感性分析环境温度对日蒸发率具有非线性影响，主要通过辐射换热（四次方关系）和气体导热（温度差）两个途径起作用。计算不同环境温度下的日蒸发率如下：表4-2日蒸发率随环境温度变化（100m³储罐，无制冷机）环境温度工况说明日蒸发率(%/d)相对变化-20℃(253K)冬季低温0.312-38.0%0℃(273K)冬季常温0.365-27.4%15℃(288K)年平均0.420-16.5%20℃(293K)标准测试0.440-12.5%30℃(303K)夏季常温0.472-6.2%40℃(313K)夏季极端0.503基准50℃(323K)极端暴晒0.545+8.4%5真空多层绝热性能衰减计算5.1真空度衰减机理高真空多层绝热容器在长期服役过程中，夹层真空度会因以下三种机制逐渐衰减，导致绝热性能劣化：（1）夹层材料放气——绝热材料（铝箔、玻璃纤维纸）、内壁表面在真空下长期释放吸附的H₂、H₂O、N₂、CO等气体；（2）大气渗漏（微漏）——焊缝、法兰密封处存在极微小漏孔，大气缓慢渗入夹层，典型漏率量级为10⁻¹⁰~10⁻⁹Pa·m³/s；（3）渗透（Permeation）——氢气等小分子气体通过外筒金属壁面的晶格渗透，对于氢介质容器这一途径不可忽略。夹层真空度随时间的变化遵循气体负荷平衡方程：P(t)=P0+(Qout+Qleak+Qperm−Qads(t))·t/V夹层(5-1)式中：P0为封结真空度，Pa；Qout为材料放气率，Pa·m³/s；Qleak为微漏率，Pa·m³/s；Qperm为渗透漏率，Pa·m³/s；Qads(t)为吸附剂吸气速率，随吸附量增加而衰减，Pa·m³/s；V夹层为夹层容积，m³；t为时间，s。5.2材料放气率计算材料放气是夹层初始真空衰减的主要来源。金属和非金属材料的放气率随时间呈幂函数衰减规律：q(t)=q10×t^(-α)，其中q10为抽真空10小时后的比放气率，α为衰减指数（金属取1.0，非金属取0.5~0.7）。表5-1夹层材料放气率参数材料名称表面积/质量q₁₀[Pa·m³/(s·m²)]α06Cr19Ni10不锈钢(内筒外表面)45.0m²2.5×10⁻⁹1.0Q345R碳钢(外筒内表面)52.0m²8.0×10⁻⁹1.0镀铝聚酯薄膜(60层)297.0m²1.2×10⁻⁷0.6玻璃纤维间隔纸4.5kg3.5×10⁻⁷0.5G10玻璃钢支撑件0.8m²5.0×10⁻⁸0.7活性炭+5A分子筛(2.5kg)—已吸收—计算得初始材料总放气率Qout≈3.2×10⁻⁶Pa·m³/s，运行1年后衰减至约8.5×10⁻⁸Pa·m³/s，5年后衰减至约2.1×10⁻⁸Pa·m³/s，长期趋于稳定。5.3吸附剂吸气能力与寿命夹层内装填低温吸附剂（5A分子筛+椰壳活性炭，共约2.5kg），在冷态（77K以下）下通过物理吸附持续捕集释放气体，维持夹层高真空。吸附剂总吸气容量有限，当吸附量达到饱和后真空度将快速恶化。吸附剂有效吸气量S_ads（对H₂等效）约为150Pa·m³/g，总吸附容量Qcap=2500g×150Pa·m³/g=3.75×10⁵Pa·m³。考虑微漏率Qleak=5×10⁻¹⁰Pa·m³/s（氦质谱检漏合格指标），年渗透量Qperm_year≈150Pa·m³/年（氢气渗透），材料年放气量Qout_year≈350Pa·m³/年（主要集中在前2年），吸附剂预计在第12~15年达到饱和，对应设计寿命15年，建议在第10~12年安排一次夹层重新抽真空维护。5.4真空度与漏热的定量关系随着真空度下降，残余气体导热从自由分子流区过渡到连续介质区，漏热量呈非线性急剧增长。以下给出不同真空度下总漏热与日蒸发率的定量计算结果：表5-2真空度衰减对日蒸发率的影响（100m³，40℃）夹层真空度(Pa)真空状态气体漏热(W)总漏热(W)日蒸发率(%/d)工程判定≤1×10⁻⁴高真空(设计)0.007156.10.503正常运行1×10⁻³高真空(冷态)0.03158.40.510正常运行1×10⁻²高真空下限0.25175.20.564正常，注意监控5×10⁻²过渡区1.1238.50.768预警，安排抽真空1×10⁻¹真空劣化2.3326.81.052告警，须维护5×10⁻¹真空失效11.2985.03.17危险，外筒结霜1×10⁰(1Pa)严重失效22.518325.90安全阀频跳≥10²完全失效≈1501120036.1绝热完全丧失5.5服役年限-性能衰减预测综合材料放气衰减、微漏累积、吸附剂饱和三个因素，建立15年设计寿命周期内的绝热性能衰减预测模型。关键节点性能预测结果如下表所示。表5-315年服役周期绝热性能衰减预测（年均温15℃）服役年限(年)预测真空度(Pa)总漏热(W)BOR(%/d)相对初始性能工程状态0(新罐出厂)5×10⁻⁴1560.42100%验收合格18×10⁻⁴1570.43100.6%正常运行32×10⁻³1600.44102.4%正常运行55×10⁻³1660.46106.4%正常运行82×10⁻²1820.50116.7%正常，建议检测105×10⁻²2200.61141.0%预警，计划维护121.5×10⁻¹3500.97226.2%建议重新抽真空15(寿命末期)5×10⁻¹9202.55594%绝热失效，需大修由表可见，在正常微漏率控制条件下（≤5×10⁻¹⁰Pa·m³/s），储氢容器在前8年可以维持良好绝热性能，日蒸发率增幅不超过20%；第10年进入性能劣化加速期，建议在第10~12年进行首次夹层重新抽真空维护，可将绝热性能恢复至接近新罐水平。6安全核验与裕度分析6.1标准符合性核验依据GB/T40060-2021《液氢贮存和运输技术要求》、GB/T29729-2022《氢系统安全的基本要求》及GB/T18442.3-2011等相关标准，对计算结果进行逐项安全核验。表6-1安全性能核验表核验项目标准/设计要求计算/设计值安全裕度结论封结真空度(冷态)≤1×10⁻³Pa5×10⁻⁴Pa50%裕度合格夹层漏率(氦检)≤1×10⁻⁹Pa·m³/s5×10⁻¹⁰Pa·m³/s50%裕度合格新罐静态日蒸发率≤0.3%/d(有制冷)0.28%/d(配制冷机)7%裕度合格新罐静态日蒸发率≤0.5%/d(无制冷)0.42%/d(年均温)16%裕度合格外壳外表面温度≥0℃(不结露)≥5℃≥5K裕度合格设计寿命内BOR≤2%/d(失效前)第15年2.55%边缘需维护吸附剂容量≥10年约12~15年充足合格安全阀泄放能力≥最大蒸发率≥10%/d当量充足合格真空维持寿命≥5年免维护约8~10年充足合格6.2真空失效安全分析极端工况下（如夹层突然破空至大气压），绝热性能完全丧失，漏热量将激增至约11.2kW，对应蒸发率约36%/d，即每小时蒸发约90kg液氢，产生BOG气体约7600m³/h。此时必须依靠双安全阀+爆破片的冗余泄放系统保障安全，同时外壳壁温将急剧下降至-100℃以下，外表面会出现严重结霜，应立即切断进出液、启动紧急排放程序、撤离现场人员并按应急预案处置。渐进式真空恶化（如微漏扩大导致真空度在数周至数月内缓慢下降）是更常见的失效模式。建议配置在线真空度监测传感器，当夹层冷态真空度恶化至5×10⁻²Pa时触发预警，1×10⁻¹Pa时触发告警并启动安全处置流程。6.3运输工况动态蒸发率修正对于罐式集装箱等移动式储氢容器，公路运输过程中的动态载荷（振动、晃动、日照辐射、环境风速）会导致蒸发率较静态值升高20%~50%。其中振动会使多层绝热层间接触压力增大，固体导热增加约10%~20%；车辆晃动引起液氢强制对流，气液界面换热增强；直射阳光可使外筒局部温度升高15~25℃，增加辐射漏热。对于40ft液氢罐箱，典型运输动态日蒸发率约0.8%~1.2%/d，运输时间应控制在7~14天以内，并预留足够的BOG缓冲空间。7结果分析与工程指导意义7.1主要计算结论（1）100m³立式液氢储氢容器采用60层高真空多层绝热、300mm夹层厚度设计，年平均环境温度下无制冷机被动储存日蒸发率约0.42%/d，配50W@20KGM制冷机时有效日蒸发率可降至0.28%/d，满足GB/T40060-2021的要求。（2）MLI绝热层漏热占总漏热的49.2%，支撑结构占29.7%，管路系统占16.0%，三者合计占比约95%。进一步降低日蒸发率应重点优化支撑结构热桥（如采用低热导率复合材料、增加热桥长度），而非单纯增加绝热层数。（3）真空度对绝热性能的影响呈强烈非线性：真空度优于10⁻²Pa时漏热基本稳定；真空度超过10⁻¹Pa后气体导热急剧增加，日蒸发率进入快速劣化区。10⁻²~10⁻¹Pa是工程上的关键预警区间。（4）在标准漏率控制（5×10⁻¹⁰Pa·m³/s）和吸附剂充分活化的前提下，储氢容器可维持约8~10年的有效高真空期，首次建议抽真空维护周期为投运后第10~12年。7.2工程设计优化建议（1）变密度多层绝热（VDMLI）优化：建议靠近热端（外筒）布置较高层密度（25~30层/cm）以抑制辐射，靠近冷端（内筒）降低层密度（12~15层/cm）以减少固体导热，较均匀布置可降低漏热10%~15%。（2）吸附剂配置优化：在冷端（内筒外壁）布置5A分子筛吸附N₂、H₂O等气体，在夹层中部适度位置布置活性炭吸附H₂、Ne等低温难凝气体；出厂前吸附剂应在250~300℃