《ENERGY CONVERSION AND MANAGEMENT》：Dual-stage Sabatier reactor with differential thermal management and membrane-based methane upgrading for high-purity methane production
编辑推荐：
本研究提出了一种采用差分化热管理策略的双级Sabatier反应器，并通过实验与数值分析，连同基于膜的甲烷提质工艺对其性能进行了评估。为克服单级反应器运行中存在的热不稳定性和平衡限制问题，研究人员将反应过程在空间上划分为动力学有利区域和热力学有利区域。所提出的热
本研究提出了一种采用差分化热管理策略的双级Sabatier反应器，并通过实验与数值分析，连同基于膜的甲烷提质工艺对其性能进行了评估。为克服单级反应器运行中存在的热不稳定性和平衡限制问题，研究人员将反应过程在空间上划分为动力学有利区域和热力学有利区域。所提出的热管理策略缓和了过度温升，稳定了反应区，并抑制了一氧化碳（CO）的生成。在化学计量比条件（H?/CO? = 4.0）下，双级反应器实现了89.8 ± 2.4%的甲烷收率（三次平行实验的均值±标准差），而在H?/CO? = 5.0条件下运行时，出口CO?浓度降至报告阈值以下（<1 vol%），甲烷收率提高至94.0 ± 1.6%。该反应器还表现出改善的热稳定性，第一级最高温度被限制在约560 °C，第二级最高温度低于320 °C。研究人员对聚砜中空纤维膜的气体渗透特性进行了实验评估，并据此开发了与组件流动布置一致的一维逆流膜模型。在不拟合混合气数据的情况下，该模型再现了实验研究的进料组成范围（H?/CO? = 3.0–5.0）和工作压力范围（1–3 bar表压）内的测量切割率、甲烷纯度和甲烷回收率。使用经实验评估的模型进行工艺模拟预测表明，在两级渗透物串联配置中，甲烷纯度和回收率可以同时达到至少95%，同时保持进料侧压力为1 bar表压，并将渗透物压力降至约0.20 bar绝压。在定义的过程边界内，化学转化效率和系统效率分别被确定为74.3 ± 2.0%和69.8 ± 1.9%；后者包括反应器加热、进料压缩以及膜相关的真空/级间压缩和循环压缩。这些发现为涉及原位资源利用（ISRU）和地面CO?甲烷化应用的基于Sabatier的系统提供了有用的设计指导。
**双级Sabatier反应器与膜法甲烷提质耦合工艺的集成设计与性能评价——论文解读**
**1. 研究背景与问题提出**
随着人类空间探索活动从近地轨道向深空拓展，原位资源利用（In-Situ Resource Utilization, ISRU）已成为保障任务可持续性和经济可行性的关键技术。将火星大气中丰富的CO?（约95%）与当地可获取的水（H?O）通过电解和CO?加氢反应转化为甲烷（CH?）和氧气（O?），被认为是技术成熟度和可行性最高的路径之一。该系统的核心在于Sabatier反应器与质子交换膜水电解器（PEMWE）的集成。然而，Sabatier反应（CO? + 4H? → CH? + 2H?O，ΔH°???K = ?165 kJ·mol?1）是强放热反应，给反应器设计和运行带来了严峻的热管理挑战。在实际条件下，催化剂床层内可能形成超过600 °C的局部热点，导致催化剂烧结和失活；同时，高温会使热力学平衡向反应物方向移动，限制甲烷收率，并促进逆水煤气变换（Reverse Water-Gas Shift, RWGS）反应生成CO，降低产物纯度。传统空间基Sabatier系统（如国际空间站生命支持系统）主要目标是CO?去除和水回收，而非高纯度甲烷生产。相比之下，火星ISRU系统对甲烷纯度、转化效率和热控可靠性提出了更高要求。已有的单级绝热反应器研究普遍存在热积累导致的转化性能受限问题，温度漂移接近600 °C，且未反应气体的存在需要额外的循环或分离单元。尽管已有研究探索了多级反应器构型、结构化催化剂和强化换热设计等方案，但仍缺乏对反应器热行为、甲烷生产性能、下游甲烷提质和整体能量效率的综合实验与建模评估。
**2. 研究目标与主要技术方法**
针对上述问题，本研究旨在通过系统实验与数值建模相结合，为高纯度甲烷生产系统建立实用的设计指南。研究设计并实验评估了一种采用差分化热管理的双级Sabatier反应器，并研究了基于聚砜中空纤维膜的下游甲烷提质工艺。主要关键技术方法包括：
（1）**双级反应器与差分化热管理策略**：将反应过程空间分离为动力学有利区（第一级，高温快速转化）和热力学有利区（第二级，低温促进平衡转化并抑制CO生成）。采用选择性隔热设计——催化剂床层区域暴露以被动散热，非反应区（进出口及连接管线）保温，无需主动冷却回路或工程化传热基底。反应器为标称1英寸不锈钢管（内径2.29 cm），每级填充17 cm床层（66.6 g）的0.5 wt% Ru/Al?O?催化剂（3.2 mm颗粒，Sigma-Aldrich），两级物理分离。整个实验活动共进行138次运行（每次3小时反应运行），单级与双级构型使用同一物理床层，未更换或重新填装催化剂。
（2）**膜分离实验与建模**：采用聚砜中空纤维膜组件（MCH-1512A, Airrane），有效膜面积1.25 m2。实验评估了CH?、CO?、H?的单组分气体渗透特性，并测定了与反应器出口组成相匹配的混合气（H?/CO? = 3.0–5.0）在1–3 bar表压下的分离性能。建立了与组件流动布置一致的一维逆流膜模型，并使用文献报告的纯气理想选择性估算H?和CO?渗透率，未对混合气数据进行拟合。此外还进行了反应器-膜直接耦合实验，将反应器实际出口气体（经冷凝和气液分离后）直接通入膜组件。
（3）**系统能量效率分析**：基于低热值（Lower Heating Value, LHV）定义了化学转化效率（ηchem）和系统效率（ηsys），后者在明确的过程边界内计入了反应器加热功率、进料压缩功以及膜相关的真空/级间压缩和循环压缩功。
**3. 主要研究结果**
**3.1 热控制特性与稳定性分析**
单级绝热反应器呈现出热点从床层中部向入口迁移的现象，平均迁移速度在10–60 min内为7.5 cm/h，且峰值温度在运行最后1.5 h仍以+1.7 °C/h的速率上升，表明3 h内未达到热稳态。相比之下，双级反应器两级热点在最初10 min后位置保持稳定（漂移小于0.5 cm），残余温度漂移分别为+0.2和?0.5 °C/h。第一级最大轴向温度梯度在整个运行期间保持约110 °C/cm。180 min时的能量平衡显示，双级系统中反应热106.0 W加上加热器输入33.2 W，仅有3.4 W转化为工艺流股的显热增益，约98%的热量通过非隔热床层区段被动排出。单级构型仅约52%的反应热通过残余热损失排出，其余被下游显热带走或在床层中积累。测得收率与各级温度下的化学平衡一致：单级收率（67.6 ± 3.2%）对应约530 °C的平衡收率；双级第一级在477 °C（最下游测温点）下达到77.7%的平衡收率，第二级在305 °C下将总收率推进至89.8 ± 2.4%，填补了约65%的剩余差距。在催化剂耐久性方面，单级反应器峰值温度达到596.6 °C且仍在上升，超过了此前记录的同一催化剂失效事件（586 °C床层最高温度）[10]；而双级反应器第一级最高温度被限制在约560 °C以下。启动性能方面，双级反应器（300 °C启动温度）的点火延迟为2.4 min（出口CH?浓度超过10 vol%），短于单级反应器（220 °C启动）的约5 min。
**3.2 反应器构型对出口气体组成和流动特性的影响**
单级绝热反应器出口含有较高浓度的未反应H?和CO?，而双级反应器的CH?浓度显著提高，且出口组成在180 min运行期间几乎保持恒定，表明其快速达到稳态并维持稳定性能。双级反应器表现出的甲烷流量更高，H?和CO?流量更低，反映出更充分的反应物利用。稳态比较（三次平行实验均值±标准差，120–180 min时均）定量确认了各组分差异。
**3.3 转化率、收率与副产物生成的定量分析**
双级反应器在H?转换率和CH?收率上均优于单级构型。CO作为副产物的分析表明，单级绝热反应器因热量积累导致高温区域扩展，为RWGS反应创造了条件，出口CO浓度在3 h内逐渐上升至约1.8 vol%。而双级反应器的CO浓度始终低于报告阈值（<1 vol%），这归因于第一级温度被限制在约560 °C、第二级低于320 °C的受控温度条件，不利于RWGS反应而有利于甲烷生成。
**3.4 H?/CO?进料比对甲烷化性能的影响**
在H?/CO? = 3.0（缺氢条件）下，H?为限制反应物，CO?转化不完全，出口CO?浓度为21.2 ± 0.2 vol%，甲烷收率仅为73.4 ± 2.3%。在化学计量比（4.0）下，性能改善，但出口仍残留3.9 ± 0.1 vol%的CO?，说明仅靠化学计量配比无法实现完全转化。在H?/CO? = 5.0（富氢条件）下，出口CO?浓度降至报告阈值以下（<1 vol%），甲烷收率提高至94.0 ± 1.6%。单级反应器遵循相同趋势但水平一致较低（3.0、4.0、5.0进料比下分别为59.0 ± 3.3%、67.6 ± 3.2%和73.8 ± 3.0%），双级反应器在各进料比下具有14–22个百分点的优势。从工艺设计角度看，H?/CO? = 5.0运行将出口气体从CH?/H?/CO?三元混合物简化为近似CH?/H?二元混合物，有利于下游甲烷纯化。
**3.5 富氢条件下的碳平衡考量**
在H?/CO? = 5.0条件下，所有三次平行运行的出口CO?浓度均低于报告阈值，而甲烷收率为94.0 ± 1.6%，表观碳平衡存在约6个百分点的缺口。误差传播分析表明，仅仪器规格不确定度即为约±3.4个百分点；允许出口CO?取至1 vol%报告阈值可将碳回收提高约2个百分点，使缺口从约6个百分点减至约4个百分点。催化剂表面含碳物种（如吸附态CO和甲酸盐物种）的暂时滞留可能是物理层面的贡献因素，但本研究中未对这些表面物种进行定量，该机制仍是合理但未经证实的次要贡献。
**3.6 CH?和H?O生产速率的定量分析**
将H?/CO?比从3.0提高至5.0，两种构型的甲烷和水生产速率均增加。在所有进料条件下，双级反应器的CH?和H?O生产速率均一致高于单级反应器。在H?/CO? = 5.0时，甲烷生产速率达到约54.1 L·h?1，水生产速率达到约86 mL·h?1。测得的生产速率与Sabatier反应化学计量关系（每摩尔甲烷生成两摩尔水）一致，H?O/CH?摩尔比约为2。
**3.7 膜系统的气体渗透特性**
H?和CO?的切割率显著高于CH?，这归因于聚合物膜的扩散控制传输机制——较小分子（H?、CO?）具有更高扩散性而更易渗透。CH?渗透率远低于H?和CO?，有利于膜分离过程中的甲烷截留。CH?渗透率在参考压力1 bar表压下测得约11 GPU（1 GPU = 10?? cm3(STP)·cm?2·s?1·cmHg?1），在所研究的压力范围内变化约±10%。假设H?/CH?和CO?/CH?理想选择性分别为56.0和22.4[48]，H?和CO?渗透率分别估计为约616 GPU和246.4 GPU。随着操作压力升高，跨膜压降（ΔP）降低，原因在于渗透通量增加减小了截留物流速，从而降低了模块内的气体速度和摩擦损失。
**3.8 两级膜工艺的数值设计与膜模型的实验验证**
建立了中空纤维膜级的一维逆流模型——沿纤维长度积分进料侧和渗透侧组分平衡，局部渗透通量由Ji = Qi·(Ph·xi ? Pl·yi)给出。模型使用基准渗透率组，未对混合气数据进行拟合。在1–3 bar表压、常压渗透侧条件下，逆流模型再现了实验测量的切割率、甲烷纯度和甲烷回收率，测量点几乎全部落在逆流和错流预测之间，定量验证了模型在实验范围内（三种反应器出口组成）的适用性。由于单级膜受甲烷纯度与回收率之间的固有权衡限制，研究采用了双级渗透物串联构型。在进料侧1 bar表压（约2.0 atm绝压）、渗透侧0.2 atm的设计点，模型预测可同时实现95.0%的甲烷纯度和98.3%的甲烷回收率（化学计量比进料），第一级和第二级膜面积在六倍缩放实验流量基础上分别为0.655 m2和0.155 m2（实际反应器规模下为0.109 m2和0.026 m2）。对缺氢和富氢出口组成，相同设计程序得出95.0%纯度下回收率分别为96.5%和98.5%。直接耦合实验表明，在反应器实际产物流量下，模块的切割率为0.81–0.94，测量的甲烷纯度（95.5–99.6%）均落在逆流/错流预测带内；H?/CO? = 5.0时的测量回收率（34.4%）与逆流预测（34.3%）吻合至0.1个百分点以内。敏感性分析显示，设计性能对假设的H?/CH?选择性比对CO?/CH?选择性更敏感；使用文献中十一组聚砜选择性对进行补充分析，固定面积甲烷纯度为89.8–98.1%，而重新定径可在每种情况下恢复95.0%纯度目标。在选定设计点，逆流与错流两种公式所需膜面积差异小于15%。
**3.9 实验混合气分离性能**
表观（混合气）选择性αMG定义为(yA/yB)/(xA/xB)。在1 bar表压下，表观H?/CH?选择性（5.81、4.59、8.17，分别对应H?/CO? = 3.0、4.0、5.0）远低于模拟中使用的理想选择性56，并随压力显著增加（2 bar表压下H?/CO? = 4.0和5.0时为26.3和34.6）。这一看似表明固有选择性丧失的行为实际上被膜模型本身使用理想渗透率所预测：在有限切割率下，快渗透组分的驱动力因其在渗透侧积累和在进料侧耗尽而崩溃。模型在相同条件下预测的表观选择性（1 bar表压下4.8、4.1、9.3；2 bar表压下20.4和38.4）复现了测量值的量级、跨进料比的排序以及随压力的陡峭依赖性。直接耦合实验在高切割率（0.81–0.94）下测得表观选择性升至30.8–37.8，比相同压力比下低切割率值高一个数量级，进一步支持了表观选择性强烈依赖操作点、不能直接解释为固有膜选择性的结论。
**3.10 能量性能分析**
在化学计量比条件（H?/CO? = 4.0）下，化学转化效率（ηchem）为74.3 ± 2.0%，仪器规格不确定度估计约±2.9个百分点。在定义的系统边界内（包括反应器加热功率33.2 W（实测）、进料压缩功6.13 W、渗透侧真空/级间压缩2.71 W、富氢第二级渗透物循环压缩至反应器压力2.27 W，均采用理想等熵模型），系统效率（ηsys）为69.8 ± 1.9%，仪器规格不确定度估计约±2.7个百分点。ηchem与ηsys之差约4.5个百分点，表明所包含的辅助功率需求相对于氢气化学能输入较小。若以膜产物中截留的甲烷（回收率98.3%）为有用产物，对应的产物基准效率为68.6%。所有压缩项假设理想等熵行为；若考虑实际小型压缩机50–80%的等熵效率，辅助功率将相应增加，ηsys最多降低约1个百分点。
**3.11 与文献报道的Sabatier反应器设计比较及放大部署意义**
与文献中代表性设计的对比表明，本研究的双级反应器在1 bar表压、未稀释进料条件下实现了89.8 ± 2.4%的甲烷收率（H?/CO? = 4.0），CO低于报告阈值，峰值床温约560 °C，且无需主动冷却回路或工程化传热基底。2.4 min的点火延迟虽短于文献报告的结构化基底ISRU反应器的约20 min冷启动时间，但两者的启动定义和初始条件不同，应谨慎比较。跨研究的直接定量排名需谨慎，因为操作条件存在异质性：压力范围约1–10 bar，[32]的转化率在约67 vol% Ar稀释下测得，催化剂装填量跨近三个数量级（0.3 g至250 g），且效率指标的系统边界不同。在放大方面，直接增大反应器管径会降低表面积与体积比，可能削弱在本规模下验证的被动散热行为；因此，平行管束编号放大比简单扩大直径更合适，但需实验验证管间热相互作用。对于火星部署，实验在地面条件下进行，未复现火星热环境；稀薄大气限制自然对流，需专门的导热和辐射散热设计。对于地面Power-to-Gas应用，预测的95% CH?产物展示了在1 bar表压进料下实现实质性甲烷提质的潜力。
**4. 讨论与结论总结**
在设计双级Sabatier反应器时，将反应过程划分为动力学有利区和热力学有利区并采用差分化热管理策略，成功解决了单级绝热操作中热点迁移、温度持续上升和CO副产物生成等问题。双级反应器实现了固定反应区、较低时间温度漂移，且在富氢进料条件下可获得94.0 ± 1.6%的甲烷收率并将出口CO?降至报告阈值以下。需要注意的是，低表观选择性在有限切割率和高渗透侧背压下是模型本身可预测的行为，不应被误解为固有膜选择性的下降，这为正确理解膜分离性能提供了理论依据。能量分析表明辅助功率需求仅使系统效率比化学转化效率低约4.5个百分点，有利于能量受限场景。研究结论如下：第一，双级反应器表现出比单级构型明显更好的时空热稳定性和更低的CO生成量；第二，进料比选择强烈影响反应器性能和下游分离要求；第三，经实验验证的逆流膜模型预测在0.2 atm渗透压、1 bar表压进料下可同时达到至少95%的甲烷纯度和回收率；第四，在所定义系统边界内，化学转化效率和系统效率分别为74.3 ± 2.0%和69.8 ± 1.9%。总体而言，所提出的双级反应器与膜法甲烷提质策略为同时提高甲烷收率和纯度并量化系统级能量性能提供了集成方法，为ISRU应用和其他CO?甲烷化系统的双级Sabatier反应器及下游提质工艺提供了实用设计指导。实际部署仍需通过编号放大、应用特定热管理、相关压力环境下膜操作、冷凝液处理和瞬态运行等方面的验证。
