基于Gibbs自由能最小化的甘油空气–蒸汽气化制合成气热力学平衡建模与无碳约束效率优化

《Biomass》：Thermodynamic Equilibrium Modeling and Optimization of Glycerol Air–Steam Gasification for Syngas Production

　　

生物柴油生产会产生约10%的粗甘油副产物，其过量供应压低了价值；粗甘油空气–蒸汽气化仍研究不足。本研究通过Cantera中的Gibbs最小化对该路线进行平衡建模，并将固体碳作为凝聚相。模型依据已发表平衡结果进行校验，并以27个独立实验点进行验证（T ≥ 500 °C时H₂摩尔分数的均方根误差（RMSE）为0.027）。研究评估温度以及蒸汽/甘油比（S/G）和空气/甘油比（A/G）对组成和热值的影响，并在两种外部热（Q）核算约定下优化考虑外部热与蒸汽焓的效率。效率在≈433 °C处呈现内部极大值（91.3%，接近自热），此时固体碳生成显著（1.02 mol/mol），无碳边界在T* = 685.1 °C处与S/G = 1相交。约束优化将实际最优点置于T ≈ 686 °C、S/G = 1、A/G = 0.5（η = 88.6%，C(s) = 0，下热值（LHV） = 14.6 MJ/kg）。粗甘油杂质调节碳生成：水使其降低（10 wt%时为0.77对1.02 mol/mol），而甲醇和非甘油有机物（MONG）使其增加（5 wt%时分别为1.07和1.26 mol/mol）；对于80/10/5/5 wt%粗甘油，以粗甘油为基准的效率约为90%。这些结果界定了粗甘油增值利用的无碳条件。

研究背景与问题：生物柴油需求受能源安全、石油储量下降和温室气体减排推动而增长，但每生产生物柴油约副产10%甘油，粗甘油供应与燃料需求 stoichiometrically 绑定，过剩供应导致其出口价格在2022—2024年明显下降。粗甘油含未反应甲醇、工艺水、催化剂中和无机盐及MONG等杂质，直接应用受限，纯化成本高。气化可将粗甘油转化为以H₂和CO为主的合成气，但单独针对粗甘油空气–蒸汽气化的平衡研究有限。既有甘油氧化蒸汽重整（OSR）和自热重整（ATR）平衡分析多以H₂产率为目标，常通过高S/G或排除固体碳规避碳生成；当目标转为燃料气能量效率并计入外部热和蒸汽焓时，低S/G、低A/G区域更有吸引力，而该区域正是固体碳热力学有利区。因此，研究人员提出：在计入蒸汽焓和外部热的效率目标下，能量最优条件会移向低S/G和低A/G，即固体碳热力学有利区，效率最大化与无碳操作存在冲突。

研究内容与意义：研究人员建立甘油空气–蒸汽气化的过程级平衡模型，以Cantera中Gibbs能最小化求解，包含固体碳凝聚相，并与已发表平衡结果和实验数据比较；开展温度、S/G和A/G对合成气组成、热值和当量比影响的参数分析；在两种Q核算约定下优化含外部热和蒸汽焓的效率；将最优解与无碳边界和自热条件对照，并进行无碳约束优化及粗甘油杂质影响量化。研究结论表明，无碳约束前效率最优位于≈433 °C、S/G = 1、A/G = 0.5，效率η_A = 91.33%，但固体碳约1.018 mol/mol，是效率上界而非实际操作点；无碳约束后实际最优在T ≈ 686 °C、S/G = 1、A/G = 0.5，η_A = 88.62%、C(s) = 0、LHV = 14.60 MJ/kg。其意义在于界定了粗甘油无碳气化的热力学条件，并明确效率最大化与无碳操作之间的权衡。

关键技术与方法：研究采用Cantera 3.2.0进行Gibbs能最小化平衡计算，气相候选物种为CO₂、CO、H₂、CH₄、H₂O、N₂和O₂，气体热化学性质取自GRI-Mech 3.0，甘油热化学取自Burcat数据库；以纯凝聚石墨相评估固体碳及气–固多相平衡。模型用已发表甘油蒸汽重整和OSR平衡结果校验，并用27个实验点验证，实验数据来源于Zhang 2007、Kotnala 2021和Chen 2011等三个独立研究的四个子集（21个气体组成和6个H₂产率）。优化以效率为目标，采用SLSQP序列二次规划，多起点策略排除局部最优；并开展无碳边界、无碳约束优化、粗甘油替代物和热力学敏感性分析。

研究结果：
3.1 模拟模型：通过基准工况（650 °C、S/G = 2、A/G = 1）的Gibbs平衡计算，得到产物气LHV为11.95 MJ/kg，并给出相分离器出口组成。
3.2 数值模型的实验验证：通过已发表平衡结果校验和27个实验点比较，发现标准反应焓与热中性O₂/甘油比在0.7%内一致，无氧化剂及稀释条件下最大H₂产率在5%内重现；T ≥ 500 °C时H₂摩尔分数RMSE为0.027、平均相对偏差3.8%，CH₄ RMSE为0.012，CO和CO₂ RMSE分别为0.048和0.058，反映CO/CO₂分布受水煤气变换催化活性影响；低于500 °C时实验偏离平衡，H₂平均相对偏差13.2%。
3.3 过程参数分析：通过温度及S/G、A/G扫描，发现H₂摩尔分数随温度升高而增加并在≈750 °C达最大（y_H₂ = 0.521），CH₄显著下降，CO持续增加；S/G升高提高H₂分数但降低合成气LHV，A/G升高增强氧化并因N₂稀释降低H₂分数和LHV。
3.4 效率优化与最优点评述：通过在两种Q约定下最大化效率，发现效率随S/G和A/G降低而升高，最优位于下界S/G = 1、A/G = 0.5；规范约定下内部温度极大值在≈433 °C（η_A = 91.33%，Q≈0），同点η_B = 88.52%，变体B最大值88.94%位于400 °C下界；该点固体碳约1.0 mol/mol，故为效率上界。
3.5 固体碳生成与无碳边界：通过含石墨的多相平衡，发现无约束最优433 °C处固体碳为1.018 mol/mol（约进料碳34%）；A/G = 0.5时避免固体碳所需最小S/G在600 °C约2.3、700 °C约0.8，边界在T* = 685.1 °C穿过S/G = 1；空气加入缩小碳生成区，边界对固体碳表示敏感。
3.6 无碳约束优化：通过以C(s) ≤ 0.005 mol/mol为约束最大化效率，得到T ≈ 686 °C、S/G = 1、A/G = 0.5，η_A = 88.62%、η_B = 86.42%、C(s) = 0、LHV = 14.60 MJ/kg；该点外部热仍为正，O₂/G = 0.105、ER = 0.03，若实现自热需A/G ≈ 5.3（O₂/G ≈ 1.11、ER ≈ 0.32），超出探索范围。
3.7 粗甘油杂质影响：通过水、甲醇和MONG替代物计算，发现在无约束效率最大点，10 wt%水使有效S/G升至1.57、LHV降至13.37 MJ/kg、固体碳由1.018降至0.770 mol/mol，效率89.4%；5 wt%甲醇提高LHV至13.88 MJ/kg并使碳增至1.074 mol/mol；5 wt%MONG使碳增至1.261 mol/mol；80/10/5/5 wt%粗甘油有效S/G为1.64、LHV为14.51 MJ/kg、粗基效率约90.4%、碳生成1.094 mol/mol。

讨论与结论：原文总结认为，该平衡模型经已发表平衡结果校验和27个实验点验证，在T ≥ 500 °C可较好重现H₂分数，但斜率与截距检验统计功效有限；参数分析表明提高S/G通过重整和水煤气变换提高H₂但稀释降低热值，提高A/G增强氧化并降低H₂和能量含量。无碳约束前效率最优在低S/G、低A/G及≈433 °C处，接近自热但固体碳显著，仅为效率上界；无碳操作在A/G = 0.5时需600 °C下S/G ≈ 2.3或近700 °C下S/G ≈ 0.8，约束最优为T ≈ 686 °C、S/G = 1、A/G = 0.5（η = 88.62%、C(s) = 0、LHV = 14.60 MJ/kg），该条件自热需A/G ≈ 5.3。粗甘油中水降低碳生成，甲醇和MONG增加碳生成，典型80/10/5/5 wt%粗甘油粗基效率约90%。原文指出局限包括：粗甘油以固定组成替代物建模；平衡方法代表热力学极限，真实反应器可能偏离，且甘油/脂肪共气化实验显示H₂和CO分数被高估，故组成和效率为上限；验证数据多对应S/G = 9–12，而最优S/G = 1附近实验数据仍少。未来应实验验证无碳条件、整合能量回收策略，并开展该路线的技术经济和生命周期评价。