一种确定生物质焦-CO2气化修正随机孔模型（MRPM）参数的方法

《Progress in Reaction Kinetics and Mechanism》：A method for determining parameters of the modified random pore model (MRPM) for biomass char-CO2 gasification

　　

生物质焦-CO₂气化的动力学建模对气化技术发展、工艺分析及设备设计优化具有重要意义。与其他模型相比，修正随机孔模型（Modified Random Pore Model, MRPM）在描述生物质气化动力学机理方面展现出独特优势。然而，文献中尚未报道MRPM模型参数的确定方法，这增加了将MRPM应用于生物质气化研究的难度。本研究提出了一种系统性求解MRPM动力学参数的方法。选用两种性质截然不同的生物质焦，即咖啡渣焦（Coffee Grounds Char, CGC）和竹焦（Bamboo Char, BBC），对所提出的方法进行验证。结果表明，采用该方法构建的MRPM能够准确拟合CGC和BBC的r-X曲线和X-t曲线。本研究提出的MRPM参数求解方法被证实是有效且可靠的。

该论文发表于《Progress in Reaction Kinetics and Mechanism》，旨在解决生物质焦-CO₂气化动力学建模中的关键瓶颈问题。研究背景植根于气化技术优化需求：动力学建模是气化器设计与反应过程效率提升的基石，多种机理模型已被提出，包括容积反应模型（Volume-reaction Model, VM）、缩核模型（Shrinking Core Model, SCM）、综合模型（Integrated Model, ICM）、随机孔模型（Random Pore Model, RPM）及修正随机孔模型（MRPM）。由于煤与生物质原料在物理化学性质上存在显著差异，不同动力学模型的适用性也因材料特性变化而不同，不同生物质类型常需专门的动力学描述以实现精确建模。现有研究表明，MRPM在描述生物质气化行为方面具有明显优势——Zhang等人针对14种生物质建立了动力学模型，发现仅MRPM能够充分描述5种具有独特反应活性曲线样品的演化规律；Gupta等人对比多种模型后证实MRPM对AAEM富集生物质焦的CO₂气化动力学表征最为准确。生物质中高碱金属和碱土金属（Alkali and Alkaline Earth Metal, AAEM，即K、Ca、Na）含量引入了复杂的气化行为，增加了动力学建模挑战。然而，MRPM虽有最适潜力，其模型参数确定方法却长期缺失，限制了实际应用。

为此，研究人员开展了系统性研究，旨在建立一种简单、稳健且普遍适用的逐步参数确定方法，并最终证实该方法对两种典型生物质焦均有效可靠，具有重要的工程应用价值。

研究所用关键技术方法如下：采用固定床反应器在N₂气氛下1,203 K制备生物质焦，选取粒径0.063–0.075 mm样品，使用无侧壁平板氧化铝坩埚（Flat Alumina Disc, FAD）进行惰性载体实验；通过程序升温至目标温度后切换CO₂/N₂混合气进行等温气化实验，采用热重分析法获取转化率（X）与反应速率（r）数据；通过极限外推法确定初始反应速率以获得速率常数K，利用低转化率区间（X < 0.3或X < 0.5）RPM方程估算结构参数Ψ，最后通过代入中高转化率区反应速率解析求解修正参数c和p。

研究结果部分按以下结构展开：

**两种生物质焦气化行为差异**：基于CGC和BBC的r-X曲线轨迹特征，研究人员发现两者的r-X曲线变化趋势较为复杂，均呈现先升高后降低的趋势，且反应速率的峰值位置和变化趋势明显不同。CGC的低转化率区Ψ值稳定在约10.9 ± 0.3，而BBC在X < 0.3时Ψ约为8.3，表明两种生物质焦具有截然不同的孔隙结构特征和反应活性演化规律。

**速率常数K的确定**：研究人员通过对MRPM动力学函数取转化率趋近零的极限，推导出速率常数K对应于X = 0时的初始反应速率。为避免气体切换滞后、瞬态流动扰动及压力波动等因素影响，研究人员将X = 0.1–0.8范围内的反应速率用光滑函数拟合，再外推至X = 0处，从而获得准确的初始反应速率。不同温度和CO₂分压下CGC-FAD-1和BBC-FAD-1的K值被成功获取，如原文Table 2所列。

**结构参数Ψ的估算**：由于MRPM由常规RPM扩展而来，低转化率区（X < 0.3）稳定的Ψ代表不受AAEM催化作用影响的固有孔隙结构。在K已确定的基础上，Ψ成为RPM方程中的唯一未知参数，通过代入特定转化率下的实验反应速率即可计算。Table 3显示，CGC-FAD-1在X < 0.5时Ψ稳定在约10.9，而BBC-FAD-1在X < 0.3时约为8.3，但两者均在高转化率区急剧上升，证实传统RPM无法描述中高转化阶段行为。

**RPM拟合局限性验证**：研究人员进一步用RPM拟合r-X和X-t曲线，结果显示RPM对CGC-FAD-1在X > 0.5时拟合精度明显恶化，对BBC-FAD-1在X > 0.3时即显不足，与Ψ变化趋势高度一致，确认了RPM在描述生物质焦中高转化率行为时的固有局限性。

**修正参数c和p的确定与MRPM验证**：采用CGC-FAD-1的Ψ = 10.9和BBC-FAD-1的Ψ = 8.3作为固定值，通过将不少于两个中高转化率点的实验反应速率代入MRPM微分方程，解析求解得到c和p值（Table 4）。CGC的p约6.3–6.5、c约0.9–1.0；BBC的p约3.2–3.4、c约1.3–1.5。将这些参数代入MRPM方程后，模型与CGC实验数据在全转化率范围内均达成极佳吻合，对BBC在X < 0.9范围内也令人满意。研究人员指出，c和p数学上补偿了中高密度转化阶段传统RPM框架无法捕获的复杂结构变化，包括孔合并、碳基质坍塌以及灰分层引起的持续增加的扩散阻力等。

**研究结论**：研究人员提出了一种系统且可行的生物质焦-CO₂气化过程中MRPM动力学参数确定方法。以两种物理化学性质迥异的典型生物质焦为验证样本，速率常数K、结构参数Ψ及修正参数c和p被成功逐步求解。结果表明，通过该方法建立的MRPM能够准确拟合CGC和BBC的全范围反应速率-转化率（r-X）曲线和转化率-时间（X-t）曲线。与传统RPM无法准确描述中高转化率气化行为相比，参数优化后的MRPM拟合精度显著提高，对不同生物质焦具有更广泛的适用性。该方法针对生物质焦-CO₂气化体系建立，也适用于高灰分和高AAEM含量的生物质焦。

该参数求解程序简洁、可靠，无需复杂迭代算法或数值优化。对于本研究中考察的两种生物质焦（咖啡渣焦和竹焦）的CO₂气化动力学建模而言，它作为一个简单有效的框架——而非普遍适用的框架。虽然外推到其他焦-气反应和其他生物质类型需要进一步的实验验证，但该程序展现出明确的前景，在经过适当验证后，有望成为气化系统动力学机理研究、反应器模拟和工艺设计的有用工具。