微流控实验与模拟揭示的地下氢甲烷化的传输与微生物控制

《International Journal of Hydrogen Energy》：Transport and microbial controls on subsurface hydrogen methanation revealed by microfluidic experiments and simulation

　　

地质甲烷化（Geomethanation），即在地质构造内通过产甲烷菌（methanogens）的微生物活动将氢气和二氧化碳原位转化为甲烷，仍是一项同时实现二氧化碳利用、氢气储存和天然气生产的有前景技术。尽管对该过程感兴趣，但其转化效率及在地下空间中的竞争力仍存在不确定性，特别是在氢气储存和碳封存需求日益增长的情况下。研究人员需要结合相关储层的详细知识与实验室测试，以评估地下环境对氢甲烷化的潜力。多个因素，如底物可用性以及盐水中必需营养物质的供应和存在，控制着微生物生长和转化效率。研究人员开发了一种耦合实验-数值方法，结合新型高分辨率微流控平台与直接数值模拟，以在孔隙尺度解析传质和甲烷产生动力学。工作流程包括微模型定殖、底物厌氧引入以及通过色谱测量进行气体分析。研究人员使用微流控反应器装置研究生物量分布和气体动力学，并校准数值模拟结果，这对气体转化和采收具有重要影响。(1) 在底物气体注入期间，微生物聚集从菌落主导的生活方式转变为浮游生活方式。(2) 观察到的菌落解体和细胞向气-液界面的迁移可能支持底物摄取和转化效率。实验产生的甲烷演化速率（methane evolution rate）峰值约为0.35 mmol/L·h。(3) 无量纲分析揭示了孔隙网络内从分子扩散到平流混合的不同传输机制。这些发现表明，微生物活性、界面传质和平流营养供应的平衡主导了非饱和流环境中的地质甲烷化过程。

**论文解读：微流控实验与模拟揭示地下氢甲烷化的传输与微生物控制**

**研究背景与问题**
在能源转型背景下，欧洲委员会等机构积极推动可再生能源技术以减少对进口化石燃料的依赖。氢气（H₂）作为可再生低碳气体，其地下储存（如多孔地质构造）被视为平衡季节性供需失衡的长期储能方案。然而，地下储存中微生物诱导的甲烷化过程（即地质甲烷化，Geomethanation）——通过产甲烷菌将H₂和二氧化碳（CO₂）原位转化为甲烷（CH₄）——虽具二氧化碳利用、氢气储存和天然气生产潜力，但其转化效率及与氢气储存和碳封存的空间竞争问题仍不明确。现有研究表明，微生物过程在枯竭气藏和深部盐水层中可影响气体组成和储存性能，但多孔介质传输与微生物活动的复杂相互作用导致氢气损失和甲烷生成的不确定性，尤其是非饱和条件下气-液两相共存时，底物可及性受气-液传质控制，而饱和条件下研究已不适用。为此，研究人员开发了耦合微流控实验与直接数值模拟的方法，旨在解析孔隙尺度下的传质、微生物动态和甲烷产生动力学，并评估传输限制对反应器性能的影响。该研究发表在《International Journal of Hydrogen Energy》。

**主要关键技术方法**
研究人员采用微流控芯片（硼硅酸盐玻璃）模拟多孔岩心结构，孔隙度0.57，渗透率2.53 D（达西），平均孔径170 μm。实验使用来自奥地利Molasse盆地Haller组现场试验点的混合微生物群落，在合成盐水（synthetic brine, SB）中培养，并于40°C、2 barg下注入40 vol.% H₂和10 vol.% CO₂的混合气。关键方法包括：(1) 高分辨率显微成像与Ilastik工具进行图像分割，量化气、液、生物量相饱和度；(2) 基于斯托克斯-布林克曼（Stokes-Brinkman）框架的直接数值模拟（GeoDict软件），结合粒子追踪统计提取界面摩尔通量；(3) 无量纲分析：通过Péclet数（Pe）、Damk?hler数（Da_I和Da_II）、Sherwood数（Sh）和Thiele模量（Φ）表征传输-反应机制。微GC（micro-gas chromatography）用于定量甲烷产量。

**研究结果**

**3.1 非饱和流实验**
通过注入接种物使微模型定殖后，观察到生物量积累导致渗透率下降，并在气体注入后形成准稳定通道。气体突破导致微生物菌落解体，分散细胞向气-液界面迁移（趋化性），菌落生活方式转为浮游生活方式。图像分割显示气体饱和度在35%-55%间波动，生物量最高占孔隙体积17%。甲烷演化速率（MER）在约320-350 h时达到0.30-0.35 mmol/L·h。实验表明，气体-液体界面临近区域生物量密度增加，支持底物摄取和转化效率提升。

**3.2 无量纲分析**
粒子追踪模拟显示：盐水中溶解的H₂和CO₂的Pe≈1-10，D/D₀≈2-10，处于分子扩散向对流混合过渡区；而生物量中底物的Pe≈10^-5-100，D/D₀≈100，以扩散为主。Sherwood-Thiele图（Sh-Φ）表明，H₂的Sh较高（10^-3-10¹），受外部传质限制较轻；CO₂的Sh低一至两个数量级，因溶解度高、扩散系数低而受边界层控制更严重。在Pe-Da^-1图上，大部分生物量和液相数据落在反应限制扩散区（Pe<10，Da_I^-1~10），而气相数据落在平流限制反应区（Pe>10）。研究确定Pe～10和Da^-1～10为最优操作平衡点，此时平流开始主导扩散且反应时间尺度略大于传输时间尺度。

**3.3 反应器性能与敏感性**
三个气体循环周期中，MER从低值持续上升至0.35 mmol/L·h。模型预测与微GC测量吻合良好，H₂从40 vol%降至<5 vol%，CO₂从10 vol%降至<2 vol%，CH₄达2-3 vol%。混合群落中的平行电子-碳汇（如发酵作用）使CH₄产率低于理想4:1化学计量比。内部有效性因子（η）与Φ的关系表明，H₂和CO₂的内部扩散限制弱（Φ?1时η≈100%），整体性能受外部传质（Sh）控制。温度升至80°C可提高扩散系数和η，而孔隙曲折度降低有效扩散系数（D_eff），降低η。

**结论与展望**
研究结论：通过发展的反应性传输模型，结合无量纲分析，研究人员在微流控芯片上成功描述了循环H₂和CO₂条件下的反应性传输动力学。主要结论包括：（1）非饱和流条件下，微生物活动、界面传质和平流营养供应的平衡主导地质甲烷化过程；（2）菌落解体和细胞向气-液界面迁移增强底物摄取和转化效率；（3）无量纲映射（Pe-Da^-1平面）结合Thiele模量和有效性因子，揭示了内外部梯度联系，最优操作点位于Pe～10、Da^-1～10。研究意义：该框架有助于识别运行点以最大化转化率和有效性，并为从孔隙尺度到更大反应器的放大提供机理基础。未来方向包括在岩心样品中深化Pe和Da表征、研究岩心中气体分散和动态流体构型、开发尺度依赖反应器评级方案以及解决现场应用工程考虑。