原《天然气化工—C1化学与化工》，经国家新闻出版署批复更名为《低碳化学与化工》。
门卓武，男，博士，教授级高级工程师，早期从事石油加氢工艺和煤直接液化油品加工技术开发和工业应用，从2004年起从事煤间接液化技术开发和工业应用，担任“十三五”、“十四五”国家重点研发计划“先进煤间接液化及产品加工成套技术开发”和“低二氧化碳选择性合成气直接制长链α-烯烃成套技术”项目负责人（首席科学家），开发的技术近5年累积创造11亿元的经济效益。近5年累计发表论文48篇，其中SCI收录19篇，EI收录9篇，其中前三和通讯作者18篇；与他人合著专著3部；授权发明专利57件，其中10件为国际专利。担任AJChE和《能源科技》编委、北京能源与环境学会京津冀专家委员会副主任委员、中国化工学会过程模拟及仿真专业委员会委员。
贺飞（1983—），博士研究生，高级工程师，研究方向为煤制油化工（煤间接液化、煤基材料和煤直接液化等），E-mail：
北京低碳清洁能源研究院煤间接液化研发技术团队成立于2017年，围绕费托合成催化剂、反应器和工艺开展基础研究、工程放大和工业应用。从催化剂基础活性相和结构出发，创新性提出了碳化铁活性相理论，开发出CNFT系列催化剂，获宁夏回族自治区2020年度科学技术进步一等奖。高质量完成国家重点研发计划“先进煤间接液化及产品加工成技术开发”项目，提升了煤间接液化行业技术水平，为促进煤炭清洁高效转化、保障国家能源安全战略提供了有力支撑。
1. 国能新疆化工有限公司，新疆 乌鲁木齐 831400；2. 北京低碳清洁能源研究院，北京 102211
摘 要 为揭示熔铁催化剂在H2气氛中的还原反应动力学行为，为还原反应器设计和放大提供基础数据，结合H2程序升温还原（H2-TPR）、气固相界面反应模型和谱图分析等方法研究了熔铁催化剂的本征还原反应动力学。熔铁催化剂还原受反应动力学控制，活化能为78.5~82.2 kJ/mol。提高还原温度能显著缩短还原时间，当还原温度不高于480 °C时，每提高10 °C，还原时间缩短大约31 min。通过4种还原反应动力学模型的比较，发现反应级数模型和几何收缩模型预测的活化能数值与实验值吻合较好（拟合度超过0.99），扩散模型和成核增长模型的预测结果较差。恒温H2-TPR还原根据结果得出，反应级数模型预测的完全还原所需要的时间和H2-TPR实验误差最小，综合看来，反应级数模型适合用于描述熔铁催化剂的还原过程。通过能量非均匀分布假设，发现活化能在还原过程中的变化较小，最大增幅不超过2.2%，说明熔铁催化剂活性组分的分布比较均匀，还原过程基本遵循能量均匀分布假设。
合成气直接转化制高的附加价值化学品（如 α-烯烃）是实现煤化工高端化、多元化与低碳化的重要方向之一。双功能催化剂和熔铁催化剂均可实现上述过程，前者存在转化率低、氢碳比高等问题，后者因具有甲烷选择性低、适用原料气范围宽、抗磨性能优越且廉价易得等优点受到了研究人员和工业界的重视。熔铁催化剂的主要成分是Fe 3O 4，Fe 3O 4被还原成 α-Fe后才具有反应活性，常用的还原气体是H 2，还原温度一般大于350 °C。研究熔铁催化剂的还原机理及反应动力学对将其进行工业应用具备极其重大意义 [ 1-2]。
目前，研究者大多通过热重（Thermogravimetric，TG）法来研究Fe 3O 4的还原机理和动力学，一般认为Fe 3O 4的还原路径与温度有关 [3]。PINEAU等 [4]采用热重法研究铁氧化物的还原时，发现当温度不高于420 °C时，Fe 3O 4直接还原为单质Fe；还原温度高于570 ℃时，Fe 3O 4首先变成维氏体（Wüstite），然后还原成单质Fe；当温度介于450 ℃和570 ℃时，维氏体和单质Fe共存。有研究者认为当还原温度不高于570 ℃时，无论是Fe 3O 4、FeO还是他们的混合物，都是一次性直接还原成单质Fe，只有当还原温度高于570 ℃时，才会出现Fe 3O 4 → FeO → Fe的还原过程 [ 5-6]。在还原反应速率控制步骤的研究中，王洪等 [7]采用热重法研究了水分压对铁基费托合成催化剂还原过程的影响，发现催化剂在还原温度较低时受内扩散控制，而在还原温度较高时则受晶相形成与生长模型控制，并且表观活化能随着水分压的增加而增大。TANG等 [8]用质谱（MS）研究了Fe 3O 4/ZrO 2和Fe 3O 4在H 2气氛中的还原反应动力学，发现Fe 3O 4还原为Fe遵循两个连续的步骤，且这两个步骤分别通过表观化学控制模型和扩散控制模型来描述。除了热重法，也有研究者采用程序升温还原（Temperature programmed reduction，TPR）来测量还原速率 [ 9-10]。和热重法相比，TPR灵敏度更高，能检测出只消耗10 -8 mol H 2的还原反应 [11]。杨锡尧 [ 12-13]详细推导了从TPR曲线获得反应活化能和指前因子的方程，并分析了升温速率和H 2消耗量对TPR结果的影响。FOLKE等 [14]利用TPR研究了助剂（Al、K和Ca）对Fe 1- xO还原过程的影响，发现助剂提高了催化剂的比表面积和孔隙率，有利于还原反应的进行，并能提升反应性能。MORA等 [15]采用H 2-TPR方法研究了氧空位缺陷载体对铁氧化物还原的推动作用，并用Avrami-Erofeyev成核模型描述了氧化铁还原反应机理。此外，用于描述还原反应动力学的数学模型也是研究重点，包括几何收缩模型（G model） [ 16-18]、成核增长模型（N model） [19]、扩散模型（D model） [8]和反应级数模型（R model） [18]等，KHAWAM等 [20]对以上诸多模型进行了详细地理论研究和比较。
目前，大多数TPR还原反应动力学研究只聚焦于反应活化能的求解，对反应动力学模型研究较少，而且大多数模型假设的活化能在还原过程中为常数，即认为催化剂表面是均匀的，但真实的催化剂表面很复杂，且某些特定部位对还原过程可能起到决定性作用，因此通过能量均匀分布假设得到的还原反应动力学模型的应用价值有限。本文利用H 2-TPR研究熔铁催化剂在H 2/Ar混合气条件下的还原反应动力学，讨论不同气固反应模型条件下的活化能以及模拟结果，并在恒温H 2-TPR还原实验中验证模型的准确性，最后考察能量非均匀分布假设对活化能的影响。
采用高温熔融法制备熔铁催化剂。精选磁铁矿、还原剂（纯Fe）以及Al 2O 3、KNO 3和CaCO 3等助催化剂按特殊的比例均匀混合后，加入电弧炉中熔融，水冷后的熔块经破碎和筛分至所需粒径范围后使用。催化剂的 n(Fe 2+)/ n(Fe 3+)通过调节磁铁矿粉与还原剂的比例控制，并通过乙二胺四乙酸（EDTA）容量法测定， n(Fe 2+)/ n(Fe 3+)为0.47。熔铁催化剂中各氧化物质量分数见表1。
H 2-TPR测试在Micromeritics Autochem II 2920化学吸附装置（美国）上进行。首先将约0.1 g样品在Ar气氛、300 ℃下预处理1 h，然后冷却至50 °C，再以一定流速向样品池中通入H 2/Ar混合气，同时以2 ℃/min将温度升至300 °C，之后以2 ℃/min的升温速率升温至850 °C，至此还原反应结束。实验过程中系统保持常压，采用热导检测器（TCD）对还原反应的H 2消耗量进行实时检测，以获得H 2-TPR数据。熔铁催化剂的H 2-TPR曲线，实验条件如下：催化剂0.12 g，还原气体10%H 2/Ar，流量90 mL/min。
熔铁催化剂在H 2气氛中的还原是一个典型的固相界面局部化学反应，也是一个无孔致密固体变成表面多孔结构的过程，还原反应式见式(1)。
还原速率与温度、气相组分分压以及还原度有关，其计算公式见式(2)。温度项 f 1( T)一般用Arrhenius方程表示，见式(3)。压力项 f 2( p)用各气相分压的幂函数形式表示，见式(4)。 f 3( x)为固相催化剂的反应模型：（1）成核增长模型（式(5)），常用于结晶过程，尤以三维Avrami-Erofeyev模型应用较广；（2）几何收缩模型（式(6)），假设反应主要发生在颗粒表面，比如在氨合成催化剂还原反应动力学中常用的球形催化剂收缩核模型 [23]；（3）扩散模型（式(7)），Jander三维模型适合转化率较低的还原反应，而基于球面径向扩散的稳态模型Ginstling-Brounshtein Model更适合熔铁催化剂的还原过程；（4）反应级数模型（式(8)），一般应用于均相动力学，但刘化章 [24]认为对于固体转化反应，反应级数 m并不是一个简单的数值。熔铁催化剂还原反应速率方程见式(9)。
式中， x为还原度（某一时刻被还原的量与可以被还原的总量之比），%； t为时间，min； p为某组分压力与系统总压之比，%； k 0为指前因子，min -1； R为气体常数，8.314 kJ/(mol·K) ； E为反应活化能，kJ/mol，一般认为在还原过程中保持不变，即遵从能量均匀分布假设 [21]； T为温度，K； n和 q为反应级数。
下文分别将成核增长模型、几何收缩模型、扩散模型和反应级数模型代入式(9)，可获得4类还原反应动力学模型。
舒季钊 [22]研究了还原气氛中水蒸气浓度对熔铁催化剂的还原速率的影响，发现水蒸气浓度增大会降低还原速率（水的氧化性），但此现象只在催化剂的还原度较高时才明显。由于实验中催化剂用量很小，和还原气量相比，还原反应过程中产生的水量极少，并且H 2也经过了脱水和脱氧，所以还原过程中的水分压可忽略 [5]。
为获得本征动力学数据，首先进行内扩散为主影响实验，确定催化剂粒径范围。选用4种粒径范围的熔铁催化剂，还原气量90 mL/min，还原气氛10%H 2/Ar，不同粒径范围的熔铁催化剂的还原度随温度的变化趋势见图2。
熔铁催化剂H 2气氛还原过程表现出了固相反应应有的动力学曲线特征：在还原初始阶段的还原速率较小，即所谓的诱导阶段，随着还原温度的升高，还原速率急速上升，进入主还原阶段，最后阶段的还原速率减小至0，整一个完整的过程的曲线呈S状。催化剂的粒径范围为25~150 μm时，其对还原速率的影响甚微，和文献 [ 25-27]研究结果相符，原因可能是熔铁催化剂的孔结构并不发达，所以粒径影响不明显。刘化章 [24]在研究合成氨催化剂纯H 2还原时，发现当催化剂粒径小于400 μm时，内扩散基本消除。WANG等 [26]在采用热重法研究铁氧化物H 2还原时发现，消除内扩散的催化剂临界粒径是240 μm。结合以上研究，本文采用粒径范围48~75 μm的熔铁催化剂。
内扩散影响排除后，进行外扩散为主影响实验，以确定最小还原气量。还原气氛为10%H 2/Ar，还原气量分别为30 mL/min、60 mL/min、90 mL/min和110 mL/min，对应H 2空速为1500~5500 mL/(g·h)。不同还原气量下的熔铁催化剂的还原度随温度的变化趋势见图3。还原气量对还原速率的影响较大，还原速率随着还原气量的增大而增大。当还原气量达到90 mL/min时，继续增大气量得到的还原度变化不大，此时外扩散的影响可以排除。
在H 2-TPR曲线中，峰温（ t m）也可拿来表征还原反应的难易程度 [17]， H 2-TPR峰温随气量的变化趋势见图4。随着气量的增大，H 2-TPR峰温从550 °C降至480 °C，说明大气量有利于还原反应的进行，气量的临界值为90 mL/min。超过临界值后，继续增大气量对还原反应影响不大。WANG等 [26]采用热重法研究铁氧化物的H 2还原时发现，消除外扩散的临界气量是88 mL/min。综上所述，本研究采用的气量为90 mL/min。
使用峰形法 [28]对H 2-TPR数据（催化剂粒径45~75 μm，气量90 mL/min）做多元化的分析，能够获得不一样的温度（300~850 °C）和时间下的还原速率（还原度为0.1%~99.8%），通过最小二乘法对数据来进行非线性拟合，得到动力学模型参数值，并对模型参数进行 F检验（ F-test）。取置信度（ α）为0.05，临界值（ F 0.05）为3.0，当模型 F值大于临界值的10倍，即动力学模型在置信度95%的区域内， F检验显著，在统计意义上无差异，不同还原反应动力学模型参数见表2。
对文献数据 [ 3-4,7,14,26,29]做多元化的分析，Fe 3O 4还原成 α-Fe的反应活化能一般为50~120 kJ/mol，而成核增长模型的活化能数值较小（37.5 kJ/mol），拟合度较差（ R 2为0.963），原因主要在于从模型假设来看，成核增长模型一般用来描述晶粒自身不断长大的类似结晶过程，要求形成的微细晶粒（ α-Fe）之间不发生粘结或长大，对温度和水蒸气浓度有限制。但是在H 2-TPR升温过程中，当还原温度达到Fe的Hutting温度（543 °C，0.3倍熔点）后，载体表面开始发生晶格表面质点的迁移，催化剂表面就会烧结，晶粒之间易发生聚并而长大。另外，还原生成的水蒸汽具有氧化性，也会促进小晶粒的长大，所以成核增长模型并不适用于熔铁催化剂的H 2还原过程，这也是模型拟合结果明显偏离实验值（图5(a)）的原因。同样，扩散模型的活化能数值太大（157.2 kJ/mol），拟合度最差（ R 2为0.916），是因为在排除内扩散和外扩散后得到的是本征动力学，扩散模型已经不适合描述熔铁催化剂的还原动力学，所以拟合结果不佳（图5(b)）。相比之下，几何收缩模型和反应级数模型的活化能数值（分别为82.2 kJ/mol和78.5 kJ/mol）和文献值比较吻合，模型拟合度较高（ R 2超过0.990），能较好地描述整个还原过程的速率变化。这也说明熔铁催化剂的还原过程是一个从表面开始，逐步向内核推进的过程，与合成氨催化剂还原动力学中常用的不变粒径收缩核模型类似 [26]。并且由于使用的催化剂的粒径较小，使得H 2在气膜和灰层的扩散速度较快，还原反应有很大的可能性成为速率控制步骤，所以几何收缩模型和反应级数模型均能很好地描述整个还原过程（图5(c)、图5(d)）。LIU等 [5]在研究合成氨催化剂还原时也认为几何收缩模型是应用最广和最重要的模型，而小颗粒催化剂还原更符合反应控制。
为考察温度对还原速率的影响，并验证反应动力学模型的准确性，进行了恒温H 2-TPR还原实验。首先将催化剂样品在Ar气氛下升温至目标温度，预处理1 h后切入还原气氛（10%H 2/Ar），当TCD信号稳定后停止实验。根据H 2-TPR曲线获得不一样的温度下催化剂被完全还原所需要的时间，并与几何收缩模型和反应级数模型的预测值比较，结果见图6。温度对还原速率的影响较大，尤其是还原温度不高于480 °C时，每提高10 °C，还原时间缩短大约31 min；当还原温度超过480 °C时，温度对还原速率的影响不明显，每提高30 °C，还原时间只缩短了约5 min。这可能是由于不一样的温度下的熔铁催化剂的还原机理各异造成的，当还原温度不超过480 °C时，Fe 3O 4直接被还原为金属Fe，当还原温度超过480 °C后，Fe 3O 4首先生成Wüstite中间相，该Wüstite中间相会在熔铁催化剂周围形成致密层，阻碍还原气体和产物气体的扩散，以此来降低了还原速率，此时继续提高还原温度对还原速率的推动作用不如低温时（480 °C以下）明显。总的来说，还原温度在380~510 °C时，反应级数模型计算得到的还原时间和H 2-TPR实验值吻合较好，而几何收缩模型在还原温度不高于450 °C时预测的还原时间和实验值差别较大。原因主要在于根据原位XRD对Fe 3O 4还原过程的分析表征根据结果得出，从前驱体到金属的还原过程，活性相发生了明显的结构变化 [ 30-31]，但根据传统的几何收缩理论模型，活性相 α-Fe的颗粒结构会继承前驱体氧化物，此假设和实际还原过程不符。综上，反应级数模型在预测还原时间上优于几何收缩模型，说明熔铁催化剂的还原过程主要受反应动力学控制。
▲ 图6 反应级数模型和几何收缩模型的完全还原时间的预测值与H2-TPR实验值比较
一般认为，活化能在还原过程中不随还原度变化，即活化能服从能量均匀分布假设。但真实的催化剂表面是不均匀的，在熔铁催化剂表面会有梯级组成的台阶，在表面的拐折位的反应速率比其他表面大得多，类似角位、边缘位和阶梯位的特定表面部位比一般晶面位在反应中更能起到决定性作用，所以催化剂表面的活性位的性质其实就是不一样的 [24,32]。严格来说，活化能受还原度的影响，有的活性位容易还原，在还原初期以很快的速率被还原（活化能较小），难还原的活性位随着温度的升高和时间的延长，在还原末期达到还原态（活化能较大）。研究者在CuO-Zn-Al 2O 3催化剂的H 2-TPR研究中引入了能量非均匀分布，建立了活化能与还原度的定量关系，包括线性和非线性，发现了类似的规律，即能量非均匀分布特征：有些部位给予较低的能量便可被还原，有些部位的还原则需要较高的能量 [ 33-34]。本文研究以反应级数模型为对象，假设活化能在还原过程中随还原度呈线)）、二次多项式（Quadratic Polynomial）函数变化（式(11)）、指数（Exponential）函数变化（式(12)）或幂律（Power）函数变化（式(13)）。
式中，下标L、Q、Ex和F分别代表呈线性函数变化、二次多项式函数变化、指数函数变化和幂律函数变化； a、 b和 c分别代表函数中的参数。
分别将式(10)~式(13)代入式(9)中，拟合后得到的模型参数见表3。
对不同的活化能模型假设，指前因子、反应级数 n和 m的数值很接近，还原初始阶段（还原度为0）的活化能差别很小，均为77.2~79.1 kJ/mol，随着还原反应的进行，活化能缓慢增大，到还原末期，活化能增大至78.3~80.5 kJ/mol。
根据表3的模型参数， 4种模型预测的活化能随还原度的变化趋势见图7，其中各模型的 RMSE均大于0.999， F 0.05检验均大于临界值的10倍。随着还原反应的进行，活化能均增加，但变化幅度较小，增幅最大的是指数函数模型，从78.3 kJ/mol增加到80.0 kJ/mol，增幅为2.2%。能够准确的看出，各模型的预测值差别较小，且均和H 2-TPR实验值接近，说明熔铁催化剂活性组分的分布比较均匀，其还原过程基本遵循能量均匀分布假设。
本文利用H 2-TPR和气固相界面反应模型研究了熔铁催化剂在H 2/Ar混合气条件下的本征还原反应动力学，得出如下结论。
（1）熔铁催化剂的还原是一个从表面开始，逐步向内核推进的过程，受反应动力学控制，活化能为78.5~82.2 kJ/mol。
（2）几何收缩模型和反应级数模型能较好地描述整个还原过程的速率变化过程。提高还原温度能显著缩短还原时间，在还原温度不高于480 °C时尤为明显。恒温H 2-TPR还原实验根据结果得出，还原温度为380~510 °C，反应级数模型预测的还原反应时间和实验值吻合程度优于几何收缩核模型。综合看来，反应级数模型更适合用来描述熔铁催化剂的还原过程。
（3）通过能量非均匀分布假如发现，活化能在整个还原过程中的变化较小，增幅不超过2.2%，说明熔铁催化剂活性组分的分布比较均匀，还原过程基本遵循能量均匀分布假设。
本研究获得的熔铁催化剂H 2气氛本征还原反应动力学模型将为还原反应器的设计和放大提供参考。