联合发射光谱与动力学计算的非热等离子体重整CO2-CH4机理

陈玉民1，卢钱程1，雷文涛1，岳文静1，张 兵1，黄文捷1，赵永椿2，张军营2

(1.中国矿业大学 低碳能源与动力工程学院，江苏 徐州 221116；2.华中科技大学 煤燃烧国家重点实验室，湖北 武汉 430074)

摘 要:非平衡等离子体重整CO2-CH4制合成气是实现2种温室气体资源化利用的新兴技术途径。综合采用发射光谱分析、反应动力学计算和连续质谱对常温常压介质阻挡放电(DBD)非热等离子体重整CO2-CH4反应途径进行探索。比功率SEI为52 J/cm3时，可得到CH4和CO2转化率最大分别为22.8%和9.4%。稳态气相产物连续质谱在线分析表明，C2H4为C2烃生成重要中间体，CH4比例越高C2烃生成量越大，且按C2H6<C2H4<C2H2顺序升高。发射光谱检测到·CH、·C2自由基和特征谱线，·CH相对强度随比功率增加明显降低，相对强度变化较小。通过动力学计算分析C2和C3产物生成消耗途径，发现89.2% C2H6来自复合反应CH3 CH3( M)C2H6( M)，78.9% C2H4主要由CH4 CHC2H4 H和C2H3 H( M)C2H4( M)反应生成，61.1% C2H2受电子碰撞反应影响形成，C2产物形成遵循路径：CH4→C2H6→C2H4→C2H2。电子碰撞离解CO2反应e CO2CO O-是CO重要形成途径，而反应和电子碰撞解离CH4是H2关键生成路径。研究表明，动力学模拟结合发射光谱可为合理揭示DBD重整CO2-CH4过程复杂反应机理提供有力工具。

关键词:非热等离子体；CO2-CH4重整；发射光谱；反应机理；化学动力学模型；等离子体化学

0 引 言

温室气体减排已成为重大科技挑战之一。开发高效、经济、洁净的碳捕集、利用和封存等技术迫在眉睫。温室气体CO2和CH4排放量巨大，干重整(DRM)可以将二者有效转化利用[1-3]。基于CH4和CO2制取合成气和高价值液态含氧化合物，对缓解能源危机、减少温室气体排放而导致的全球气候变暖等环境问题具有重要战略意义[4-6]。

由于CH4和CO2均为化学性质非常稳定的分子，其中C—O和C—H化学键能分别为532.2和439.3 kJ/mol，实现二者活化需高温、高压条件，过程能耗高[3]，且易导致催化剂积碳和烧结失活[7-8]。因此如何在温和件下实现二者活化是目前CO2-CH4转化利用最具挑战性的工作之一。

等离子体作为物质第4种状态，内部含有丰富的高活性带电粒子(电子、阳离子和阴离子)和电中性粒子(原子、分子和自由基)，化学活性极强[9-11]。大气压非热等离子体电子温度非常高(1～10 eV)，具有足够能量使CH4和CO2分子激发、离解、电离，而气体却可保持较低温度(接近室温)，从而有效避免传统热催化带来的高温问题。另外，高能电子也能扮演化学“催化剂”角色，实现温和条件下CH4和CO2活化[12-15]。研究表明不同类型非热等离子体源，如辉光放电[16-17]、介质阻挡放电(DBD)[18-21]、微波(MW)放电[22-24]、火花放电[25-27]和电弧放电[28-30]在非热等离子体辅助燃烧和燃料重整中的有效性。其中，介质阻挡放电等离子体由于具有装置简单、可操作性高、电子能量高等特点，在CO2-CH4转化利用中应用潜力巨大[31-33]。

非热等离子体辅助重整CO2-CH4在生产、消耗中间体和反应物时可能发生数百种反应，反应过程复杂。目前，主要通过试验诊断和数值模拟对CO2-CH4重整反应进行机理研究。DAI等[34]通过脉冲放电对CH4/CO2反应体系进行研究，发现反应气中加入定量CO2组分可有效促进CH4转化，但CO2添加量过大将抑制C2烃产物生成。ZHOU等[35]研究发现常压DBD放电下合成气中H2/CO体积比主要取决于进气中CH4/CO2体积比，并发现富CO2混合气能有效防止积碳。王晓玲等[36]研究表明脉冲特性(重复频率和脉宽)对DRM转化特性有显著影响。KOZLOV等[37]对常压DBD放电转化CH4/CO2过程进行机理分析和模拟研究。OZKAN等[38]阐述了介质阻挡放电转化CH4/CO2过程中间物质和高附加值产物最合理的生成消耗机理，共包含26步反应。WANG等[39]研究总结了较详细的CH4-CO2-O2-H2O等离子体化学反应机理，忽略了激发态粒子化学反应行为。SUN等[40-41]在模拟中考虑了激发态CO2、CH4反应行为，忽略了DBD脉冲非连续放电实际行为特征及其对转化路径的影响。此外，若将中间物质演变试验测量规律与机理仿真相结合，可更准确解析等离子体重整CO2-CH4机理。

笔者开展了介质阻挡放电非热等离子体重整CO2-CH4试验研究，并结合质谱仪稳态产物连续分析、发射光谱原位检测和瞬态Chemkin动力学计算，对DBD重整CO2-CH4反应机理进行探索。通过反应物和主要生成物连续质谱信号，分析了CH4-CO2重整产物变化规律；进气组分和比功率(SEI)对CH4、CO2转化率和产物CO、H2以及C2、C3烃生成效率影响。采用发射光谱原位检测CO2-CH4等离子体中自由基、粒子相对光谱强度随CH4比例和SEI的变化趋势及相应机理；分析了瞬态非热等离子体重整CO2-CH4动力学计算模型，分析了CO2-CH4反应路径和主要稳态产物生成途径。

1 试 验

试验在常压、常温下进行，试验系统如图1所示，主要由DBD反应器、交流激励电源、质量流量控制器、光谱仪、质谱仪(RGA)组成。DBD反应器采用圆盘形电极结构，石英玻璃作为阻挡介质，石英容器放电区域内部直径60 mm，放电间隙8 mm，石英厚度1 mm，放电区域体积22.608 cm3。混合气体组分为N2 CO2 CH4，其中CH4流量占比为x%(x=5.0、7.5、12.5、30.0、50.0、70.0)，N2流量为100 mL/min，混合气总流量为200 mL/min。反应器输入功率由交流激励电源(Coronablelab, CTP-2000K)控制，最大放电频率为100 kHz，总有效输入功率为90～130 W(功率间隔为10 W)。反应器出口处气体利用质谱仪(Tilon，SRD300M)在线测量主要产物CH4、CO2、CO、H2、C2烃和C3烃浓度。同时利用光谱仪(AVANTES，ULS2048L-USB2)测量CO2-CH4非热等离子体中自由基、反应物中间体特征谱带动态变化。

图1 试验装置示意
Fig.1 Schematic of the experimental setup

采用CH4、CO2转化率表征试验和模拟结果：

(1)

(2)

式中，Fin(CH4)、Fin(CO2)分别为CH4和CO2的进气流速，mL/min；Fout(CH4)、Fout(CO2)分别为CH4和CO2的出气流速，mL/min。

比功率(Specific Energy Input,SEI，J/cm3)定义为输入功率与气体流速比值，具体为

SEI=60P/F，

(3)

式中，P为输入功率，W；F为气体总流速，mL/min。

SEI为52 J/cm3、N2、CO2、CH4体积比为100∶50∶50稳定放电后电压电流波形图如图2(a)所示。可知此功率下放电电压为20 kV，放电峰值电流约20 mA，频率为9.3 kHz。计算得到单个连续脉冲时间约50 μs。由于单个微放电持续时间难以测量计算，后续模拟中设定为40 ns的矩形波。SEI为36、44和52 J/cm3时放电Lissajous图如图2(b)所示。通过图形面积法计算得到输入SEI为52 J/cm3时放电SEI为45.91 J/cm3，电源输出效率为88.3%。

图2 DBD反应器放电波形图
Fig.2 Diagram of discharge waveform of DBD reactor

2 化学动力学反应机理模型

基于CHEMKIN PRO 17.0中0D Plasma-PSR模型模拟了DBD放电激励下CO2-CH4等离子体反应过程，其为连续反应，但忽略气体扩散、流动时间。反应动力学机理模型包括电子碰撞反应和重粒子(原子、分子、自由基、离子和激发态物种)反应，包含92种物种，831个反应，详细物质见表1(v为振动量子数)。激发态物质包括N2、CH4、CO2和H2，这些分子吸收一定能量后，电子由低能级激发到高能级形成激发态。电子碰撞反应速率常数是基于电子碰撞截面数据，利用BOLSIG 软件[42]求解Boltzmann方程得到。电子-离子和重粒子重组反应阿伦尼乌斯速率系数取决于电子、气体温度，参数取自文献和NIST化学动力学数据库[43]。在动力学建模过程中忽略了电极附近电荷分离和鞘层形成。具体模拟工况见表2。

表1 详细动力学机理所包含物质

Table 1 Substances involved in detailed kinetic mechanism

表2 模拟工况条件

Table 2 Simulation conditions

3 结果与讨论

3.1 产物分布规律

CH4体积分数5%～70%和SEI为40～52 J/cm3时，非热等离子体重整CO2-CH4过程中产生气体组分与室温可凝结组分(GC-MS测量液相产物为CH2O、CH3OH、HCOOH和微量脂类)占比变化如图3所示。可知CH4体积分数为5.0%、7.5%和12.5%时，在不同SEI输入下反应后产物大部分为气体组分。表明CH4体积分数较低时，通过电子碰撞和复合只能生成H2、CO、C2H2、C2H4和C2H6等气体物质，且生成C2烃物质相对较少，无法继续由C2物质经重组反应生成高碳烃产物。不同SEI输入下，CH4体积分数由30%增至70%时，可凝结组分比例随CH4体积分数升高而增大。这是由于CH4体积分数较高时通过电子碰撞反应和自由基复合生成较多C2、C3烃，C2和C3碳氢化合物继续发生复合反应形成高碳烃产物。反应过程中并未明显生成固体碳，可能是由于CO2富余氧化了C物质。

图3 不同CH4比例和SEI下产生气态组分与可凝结组分占比
Fig.3 Proportions of gaseous and condensable components produced at different CH4 ratios and SEI

SEI为52 J/cm3时，测得CO、H2质谱信号随反应时间变化如图4所示。由图4可以看出，CH4体积分数为30%时，H2消耗比较明显，原因可能为H2经电离生成·H自由基与CH4解离生成的·CH、·CH2和·CH3 等自由基反应形成C2烃。然而，产物气中几乎未检测到O2。这是由于CO2解离形成O原子将进一步反应形成含氧化合物，只有小部分复合成O2。

图4 SEI为52 J/cm3时H2和CO生成随反应时间的变化
Fig.4 Change of H2 and CO generation with reaction time at SEI of 52 J/cm3

3.2 原料转化率

不同CH4比例下SEI对非热等离子体重整CO2-CH4过程中CH4、CO2转化率的影响如图5所示。在相同CH4体积分数下，随SEI增加，CH4和CO2转化率随之增大，CH4转化率在10.2%～22.8%，CO2转化率在0.6%～9.4%。

图5 非热等离子体重整CO2-CH4试验中CH4比例、SEI对CH4和CO2转化率的影响
Fig.5 Effect of CH4 ratio and SEI on CH4 and CO2conversion in non-thermal plasma reforming CO2-CH4 experiment

这是由于输入功率增加，放电区域电场强度增大，相应活性粒子数量及高能电子数量增加，且电子具有的能量得以提高，有利于CO2和CH4分解。同时，由于电子密度通常随输入功率增加而增加，从而转移更多能量解离CH4和CO2(如e CH4e CH3 H、e CH4e CH2 H2、e CH4e CH H2 H和e CO2e CO O等反应)。另外，电子与N2碰撞产生激发态(如N2(A)、N2(B)、N2(C)和N2(a′)等)也能促进一部分CO2和CH4分解。

在相同输入功率、不同CH4比例下，CH4转化率先增加后减小，且在不同SEI下，CH4比例由5%增至70%时，均在CH4体积分数30%时得到最大CH4转化率。CO2转化率随CH4比例和输入功率增加而增大，CH4体积分数为70%、SEI为52 J/cm3时达到最大转化率9.4%。这是由于CH4体积分数增加时CO2体积分数减少，电子碰撞解离和激发态N2碰撞激发促进了CO2分解。可见，不同CO2-CH4比例下，CO2和CH4参与等离子体重整反应中间反应过程有较大差异。另外，CH4体积分数较小时(CO2过量)，CO2分解产生O原子不仅有利于CH4分解，同时过量O原子还能与·H、C反应。

相同CH4比例和SEI下，CH4转化率始终大于CO2，这是由于高能电子能量达0.162和0.361 eV时，CH4可通过2个振动通道(CH4(v24)和CH4(v13))诱导断裂第1个C—H键，发生离解反应:e CH4e  CH4(v24)、e CH4(v24)e CH3 H和e CH4e  CH4(v13)、e CH4(v13)e CH3 H；而高能电子能量达到5.45 eV时，才能通过电子碰撞使CO2断裂1个CO键，发生解离反应e CO2e CO O生成CO和活性O基团[6]。由图5(b)可知，较高CH4体积分数会提高CO2转化率。

3.3 C2、C3烃生成规律

SEI为36、52 J/cm3时，CH4体积分数对C2和C3烃生成率的影响如图6所示。反应结束时C2H2、C2H6和C3H6质谱信号强度越大代表生成量越多。SEI由36 J/cm3增至52 J/cm3，C2和C3烃生成率逐渐增大，这是因为输入功率增加促进CH4解离产生更多·CH3、·CH2、·CH，进一步复合生成更多C2H2、C2H4、C2H6和C3H6等碳氢化合物。同理，CH4体积分数由5%增至70%、输入功率相同时，C2和C3烃生成率增加很明显。CH4体积分数为5.0%、7.5%和12.5%时，C2H2和C2H6生成率较低；归因于CH4 体积分数较小时，过量CO2会氧化C2烃，导致C2烃生成率降低。另外，CH4体积分数较小时并未生成C3H6，这是由于部分C2H4和C2H6被氧化成含氧化合物无法继续复合生成C3H6。

图6 CH4比例和SEI对生成物C2和C3烃生成率影响
Fig.6 Effect of CH4 ratio and SEI on the formation rate of C2 and C3 hydrocarbons

CH4体积分数为5.0%、7.5%和12.5%时，C2H2、C2H4和C2H6生成率为：C2H2<C2H6<C2H4；CH4体积分数为30%、50%和70%时，C2烃生成率为：C2H2<C2H4<C2H6。后续对其形成途径进行详细讨论。

3.4 CO2-CH4等离子体发射光谱

为进一步了解CO2-CH4在非热等离子体重整解离过程中反应机理，利用发射光谱对不同工况下DBD CO2-CH4等离子体发射光谱进行原位测量。CH4体积分数为50%、SEI为36和52 J/cm3时典型发射光谱(光谱强度均为相对强度，按336.06 nm最高强度值进行归一化处理)如图7所示。

图7 DBD CO2-CH4等离子体发射光谱
Fig.7 DBD CO2-CH4 plasma emission spectroscopy

由图7可知，在DBD放电条件下，N2-CO2-CH4中产生了N2第2正带系第1负带系、·H、·CH、·C2、·CN、CO和等物质。其中N2第2正带系谱带强度较明显。检测到激发态离子谱带，该离子可能是由于电子直接碰撞激发导致和CO(c3∏→a3∏)等带系存在证明生成了主要产物CO。位于390、431 nm CH(B2∑→X2∏)、CH(A2Δ→X2∏)带系相对强度虽较小，但可反映等离子体重整过程中有高碳烃形成，而·CH 主要通过e CH4e CH3 H、e CH3e CH2 H、e CH2e CH H、e CH4e CH H2 H等反应产生。460～520 nm表明存在C2 Swan带系，主要对应C2产物，这些产物通过含C粒子(·CH、·CH2、·CH3、C2Hx及C3Hx等)之间复合反应及高能电子与这些粒子之间碰撞反应生成，该带系与积碳密切相关。试验中，其发射强度相对较弱，说明在DBD中积碳较少。·CN自由基可通过N2 CCN N、N2 C22CN等反应形成。由于激发态H易被其他粒子碰撞而淬灭，寿命极短，因此采用光谱仪难以检测到·H特征谱带。实际上CH4电子碰撞离解也会产生·CH3、·CH2等自由基，但本研究中并未检测到，这是因为CH3和CH2激发态发射波长在红外区。

不同SEI、CH4比例下CH和相对强度变化如图8所示。可知SEI和CH4比例对·CH(431 nm)相对强度有较大影响，而对影响较小。SEI为36 J/cm3时，·CH相对强度随CH4体积分数增加呈先增加后减小趋势，与CH4转化率变化趋势一致，CH4体积分数为5%和30%时其相对强度分别达最小值0.035和最大值0.064。SEI提高至52 J/cm3时，·CH相对强度明显变小，这是因为SEI增加提高了电子能量和密度，促进了·CH进一步转化。

图8 SEI对粒子发射光谱相对强度的影响
Fig.8 Effect of SEI on the relative intensity of CH and (370.4 nm) particle emission spectroscope

通过CO2-CH4等离子体发射光谱(图7)可知存在大量高振动激发态离子带系。这是由于电子碰撞将CO2激发到高电子或振动态，降低了CO2电子碰撞离解阈值能量(阈值能量由11 eV降至5.5 eV)，从而促进CO2离解形成CO、O2和物种。这也是粒子相对强度(图8(b))随CH4比例和SEI增加变化较小的原因。

3.5 动力学模拟分析

3.5.1 动力学模型验证

为验证动力学模型可靠性，对比了输入功率130 W、CH4体积分数30%、50%和70%工况下，试验和模拟CH4、CO2转化率，如图9所示，可知模拟结果与试验结果变化趋势一致，13.6 s计算时间后的模拟与试验结果误差<5%，详细对比见表3。

图9 不同CH4比例下CH4和CO2转化率试验与模拟结果对比
Fig.9 Comparison of experimental measurements and simulation results of CH4 and CO2 conversion in different CH4 ratios

动力学计算得到的CH4、CO2转化率偏高，主要原因为设置的放电功率过高、估算的单个微放电时间与实际有差异。

3.5.2 电子温度和电子质量分数

1次、10次微放电脉冲电子温度和电子质量分数随时间变化如图10所示。10-8 s时，施加40 ns矩形脉冲，功率沉积为1.2×106 W/cm3。电子温度最大约3.0 eV，发射光谱测量计算得到不同功率下电子温度Te为1.8～3.2 eV。10次脉冲放电时，电子温度和电子质量分数变化趋势相同。

图10 电子温度和电子质量分数随时间变化
Fig.10 Change of electron temperature and electronmass fraction with time

3.5.3 等离子体物种质量分数

1)CH4、CO2和H2、CO质量分数。CH4、CO2和H2、CO在1次、10次脉冲放电下质量分数随时间变化如图11所示。由图11(a)可知，1次脉冲下，CH4和CO2主要由电子碰撞解离、电子碰撞复合反应和自由基反应消耗，尤其是前2个阶段。对比10次脉冲放电，随反复放电CH4和CO2消耗更明显，H2和CO质量分数随脉冲次数增加也逐渐增大。

图11 CH4、CO2和H2、CO质量分数随时间变化
Fig.11 Change of CH4, CO2 and H2, CO mass fraction with time

2)稳态产物C2和C3烃质量分数。不同脉冲放电期间形成的重要稳态产物C2和C3烃(主要包括C2H2、C2H4、C2H6、C3H6和C3H8)质量分数随时间变化如图12所示。放电开始时，C2H6、C2H4和C2H2质量分数瞬时上升，而C3H6和C3H8质量分数上升出现延迟，这是因为C3分子只有在某些中间产物形成后才能产生。C2形成量排序为C2H6>C2H4>C2H2。

图12 稳态产物C2和C3质量分数随时间变化
Fig.12 Change of mass fraction of steady-state products C2 and C3 with time

3)自由基质量分数。不同脉冲次数下重要自由基质量分数随时间变化如图13所示。在单个脉冲期间，重要自由基为·CH3、·CH2、·H和·O，主要由CH4和CO2经电子碰撞解离直接形成。由图13(a)可知，放电诱发产生上述自由基质量分数达到极大值后由于自由基间的反应逐渐消耗而降低。但10个脉冲中，由于反复连续放电和初始反应物浓度较高，自由基质量分数下降速度变缓，甚至呈整体持续上升趋势。可见，准确模拟连续重复脉冲放电对合理描述DBD等离子体反应行为至关重要[37-39]。

图13 等离子体中形成的重要自由基质量分数随时间变化
Fig.13 Change of mass fraction of important radicals formed in the plasma with time

4)离子质量分数。不同脉冲次数下主要离子随时间分布如图14 所示。形成的离子可分为2组：和在脉冲期间形成，但在余辉中寿命也相当长，而和主要在余辉中形成。在余辉中有明显消耗趋势，这是因为发生了反应；同时可由生成。对于主要由和产生，其消耗过程为和

图14 等离子体中形成的重要离子质量分数随时间变化
Fig.14 Change of mass fraction of important ions formed in the plasma with time

3.5.4 CO2-CH4反应路径

通过分析各物质生成消耗路径，得到模拟工况为SEI=52 J/cm3、CH4体积分数50%时反应路径如图15所示。可知CH4消耗主要通过离子、自由基反应，少部分直接由电子碰撞离解和电离。CH4离解主要生成·CH3、·CH2、·CH和·H等自由基，电离生成等离子。与CH4反应生成与C2H6重整主要生成与电子重组生成·CH3、·C2H3等基团；与电子重组生成·CH和·CH2(省略了·CH和·CH2)；可见电子碰撞CH4电离生成正离子最终都能转变为活性自由基。CH4与电子碰撞过程生成·CH3最多；大部分CH3通过反应H CH3( M)CH4( M)再形成CH4(CH4最主要生成来源)；89.2% CH3发生偶联反应生成C2H6，还有一部分·CH3与·C2H3复合生成C3H6。

图15 CO2-CH4反应路径分析
Fig.15 Analysis of CO2-CH4 reaction pathway

CO2消耗主要是电子碰撞激发形成激发态和离解成CO。含氧碳氢化合物CH3OH主要由·CH3与·OH反应生成，少量来源于·H与CH2OH、CH3O反应。发现含氧化合物生成消耗中都有·H自由基参与。另外，CH4和CO2部分消耗过程受到N2激发态影响，这说明加入N2能促进CO2和CH4转化。图15可了解H2、CO和C2、C3碳氢化合物及C2、C3中间物质的生成消耗路径以及对应贡献率。

3.5.5 C2和C3烃生成消耗途径

通过动力学计算分析了SEI=52 J/cm3、CH4体积分数50%模拟工况下，产物中C2H2、C2H4、C2H6和C3H6生成和消耗途径，结果如图16所示。可知由2个CH3结合生成C2H6反应R507对C2H6生成贡献率最大为89.2%，高能电子e与CH4碰撞离解生成CH3对C2H6生成促进作用最明显；H与C2H5基团结合反应R585对C2H6生成也有一定促进作用(贡献率为10.8%)。而与C2H6重组生成和CH3(R187)对C2H6消耗量最大，贡献率为90.3%；N2激发态N2(A)、N2(a′)能使部分C2H6分解。

图16 C2和C3烃生成消耗分析
Fig.16 Analysis of generation and consumption of C2 and C3 hydrocarbons

CH与CH4反应R475以及·H和·C2H3结合反应R621是C2H4主要生成途径，贡献率分别为40.2% 和38.7%；N2(A)、N2(a′)与C2H6反应R118、R146也会生成小部分C2H4和H2。而C2H4与·H反应R610生成C2H5及其与反应R188是C2H4消耗途径，贡献率为71.8%和28.2%。

电子e与和碰撞反应及C2H与CH4反应R480是C2H2主要生成途径，贡献率分别为32.5%、15.4%、11.2%、2.12%和5.89%；C2H2消耗途径主要是反应R646和反应R189，贡献率为57.9%和42.1%。

C3H6主要由·CH3和·C2H3结合反应形成；部分C3H6由·OH、·H、·CH2、·CH3自由基碰撞C3H7，使C3H7上1个C—H键断裂，·OH、·H、·CH2、·CH3自由基再夺取H基形成。·CH3、·CH2、·H则主要由电子与CH4碰撞解离产生。因此CH4离解产生·CH3自由基在随后C2H2、C2H4、C2H6和C3H6等形成中起主要作用。C2烃形成遵循CH4→CH3→C2H6→C2H4→C2H2转化途径，与试验结果吻合。

3.5.6 H2和CO生成消耗途径

计算得到H2和CO生成消耗路径如图17所示。由图17(a)可知，H2消耗反应中对H2消耗量很小，其中R207和R82反应主导了H2生成，生成贡献率分别54.6%和20.1%。另外，电子e与发生碰撞离解反应也生成部分H2；N2(a′)与CH4生成·CH2、C反应R139、R141对H2生成也有一定贡献，分别为9.5%、1.9%。由图17(b)可知，CO主要生成路径是高能电子e碰撞CO2分子发生离解反应R21(50.8%)；另一条路径为H碰撞激发态CO2(v1)分子使其断裂1个CO键生成CO和·OH，贡献率为26.1%。N2(A)和N2(B)与CO2反应R124和R102对CO生成也有一定促进作用。H2和CO生成主要受到电子碰撞离解反应的影响。

图17 H2和CO生成消耗分析
Fig.17 Analysis of generation and consumption of H2 and CO

4 结 论

1)考察不同输入功率、进气组分对等离子体重整效果的影响。提高输入功率可有效提高CH4、CO2转化率，促进CH4和CO2电子碰撞离解反应；随CH4体积分数增加，二者转化率先增加后减小，且在CH4体积分数为30%时获得最大CH4转化率22.8%。CH4体积分数较高时，会促进C2和C3烃生成，同时生成更多高碳烃等。C2烃生成量为：C2H6<C2H4<C2H2。

2)发射光谱检测到CO2-CH4等离子体中有·CH、·C2自由基和物种。随SEI增加，CH相对强度明显下降，这是由于·CH继续复合生成大量C2、C3碳氢化合物，而相对强度变化不大。

3)构建的动力学机理可准确预测非热等离子体重整CO2-CH4反应行为。由反应路径分析可知，CH4经电子碰撞电离生成的正离子最终均转变为活性自由基。其中，89.2% CH3发生偶联反应生成C2H6；CO2消耗主要是电子碰撞激发形成激发态和离解产生CO。CH4离解产生·CH3自由基在C2H2、C2H4、C2H6和C3H6等形成中起重要作用。C2碳氢化合物形成遵循转化途径CH4→CH3→C2H6→C2H4→C2H2，与试验结果相互验证。

参考文献(References):

[1] QUÉRÉ C L, ANDREW R M, FRIEDLINGSTEIN P, et al. Globalcarbon budget 2017[J]. Earth System Science Data, 2018, 10(1): 405-448.

[2] WANG L, YI Y H, WU C F, et al. One-step reforming of CO2 and CH4 into high-value liquid chemicals and fuels at room temperature by plasma-driven [J]. Angewandte Chemie International Edition, 2017, 56(44): 13679-13683.

[3] PAKHARE D, SPIVEY J. A review of dry (CO2) reforming of methane over noble metal catalysts [J]. Chemical Society Reviews, 2014, 43(22): 7813-7837.

[4] JU Y G, SUN W T. Plasma assisted combustion: Dynamics and chemistry [J]. Progress in Energy and Combustion Science, 2015, 48: 21-83.

[5] JU Y G, LEFKOWITZ J K, REUTER C B, et al. Plasma assisted low temperature combustion [J]. Plasma Chemistry and Plasma Processing, 2015, 36(1): 85-105.

[6] CHUNG Wei-Chieh, CHNAG Moo-Been. Review of and plasma on dry reforming of CH4 and possible synergistic effects [J]. Renewable and Sustainable Energy Reviews, 2016, 62: 13-31.

[7] 陈玉民，卢钱程，任燕燕，等. Ni-xFe/Mayenite促进CO2-CH4重整的稳定性调控实验研究[J]. 燃烧科学与技术，2022，28 (3): 229-238.

CHEN Yumin, LU Qiancheng, REN Yanyan, et al. Experimental investigations on stability tailoring of Ni-xFe/Mayenite for CO2-CH4 reforming[J]. Journal of Combustion Science and Technology, 2022, 28(3): 229-238.

[8] 吴兴亮，吕凌辉，马清祥，等. 甲烷二氧化碳重整镍基催化剂的研究进展[J]. 洁净煤技术，2021，27 (3): 129-137.

WU Xingliang, LYU Linghui, MA Qingxiang, et al. Research progress of nickel-based catalysts for carbon dioxide reforming of methane[J]. Clean Coal Technology, 2021, 27(3): 129-137.

[9] 朱丽华，张悦，田瑶瑶，等. 电源参数和气体组分对低温等离子体转化煤层甲烷的影响[J]. 洁净煤技术，2022，28 (5): 81-86.

ZHU Lihua, ZHANG Yue, TIAN Yaoyao, et al. Effects of power supply parameters and gas composition on transformation of coal bed methane by low-temperature plasma[J]. Clean Coal Technology, 2022, 28(5): 81-86.

[10] 戴栋，宁文军，邵涛. 大气压低温等离子体的研究现状与发展趋势[J]. 电工技术学报, 2017, 32(20): 1-9.

DAI Dong, NING Wenjun, SHAO Tao. Research status and development trend of atmospheric pressure low-temperature plasma[J]. Journal of Electrotechnical Technology, 2017, 32(20): 1-9.

[11] 邵涛，章程，王瑞雪. 大气压脉冲气体放电与等离子体应用[J]. 高电压技术, 2016, 42(3): 685-705.

SHAO Tao, ZHANG Cheng, WANG Ruixue. Atmospheric pressure pulsed gas discharge and plasma application [J]. High Voltage Technology, 2016, 42(3): 685-705.

[12] GUO X, FANG G, LI G, et al. Direct,nonoxidative conversion of methane to ethylene,aromatics,and hydrogen [J]. Science, 2014, 344: 616-619.

[13] LIU H, WIERZBICKI D, DEBEK R, et al. La-promoted Ni-hydrotalcite-derived catalysts for dry reforming of methane at low temperatures [J]. Fuel, 2016, 182: 8-16.

[14] LOENDERS B,ENGELMANN Y, BOGAERTS A. Plasma-catalytic partial oxidation of methane on Pt(111): A microkinetic study on the role of different plasma species [J]. The Journal of Physical Chemistry C, 2021, 125(5): 2966-2983.

[15] ABBAS H F, DAUD W W. Hydrogen production by methane decomposition:A review [J]. International Journal of Hydrogen Energy, 2010, 35(3): 1160-1190.

[16] LI D, LI X, BAI M, et al. CO2 reforming of CH4 by atmospheric pressure glow discharge plasma: A high conversion ability [J]. International Journal of Hydrogen Energy, 2009, 34(1): 308-313.

[17] GHORBANZADEH A M, MATIN N S. Methane conversion to hy-drogen and higher hydrocarbons by double pulsed glow discharge [J]. Plasma Chemistry and Plasma Processing, 2005, 25(1): 19-29.

[18] TU X, WHITEHEAD J C. Plasma-catalytic dry reforming of met-hane in an atmospheric dielectric barrier discharge: Understandi-ng the synergistic effect at low temperature [J]. Applied B:Environmental， 2012, 125: 439-448.

[19] KHOJA A H, TAHIR M, AMIN N A S. Dry reforming of meth-ane using different dielectric materials and DBD plasma reactor configurations [J]. Energy Conversion and Management, 2017, 144: 262-274.

[20] MEI D H, ZHANG P, DUAN G H, et al. CH4 reforming with CO2 using a nanosecond pulsed dielectric barrier discharge plasma [J]. Journal of CO2 Utilization, 2022, 62: 102073.

[21] ZHENG X, TAN S, DONG L, et al. LaNiO3@SiO2 core-shell nano-particles for the dry reforming of CH4 in the dielectric barrier discharge plasma [J]. International Journal of Hydrogen Energy, 2014, 39(22): 11360-11367.

[22] FIDALGO B, DOMINGUZ A, Pis J J, et al. Microwave-assisted dry reforming of methane [J]. International Journal of Hydrogen Energy, 2008, 33(16): 4337-4344.

[23] PHAM T P, RO K S, CHEN L, et al. Microwave-assisted dry reforming of methane for syngas production:A review [J]. Chemistry Letters, 2020, 18(6): 1987-2019.

[24] MINEA T, BEKEROM Dirk C M, PEETERS Floran J J, et al. Non-oxidative methane coupling to C2 hydrocarbons in a microwave plasma reactor [J]. Plasma Processes and Polymers, 2018, 15(11): 1800087.

[25] SCAPINELLO M, DELIKONSTAINTIS E, STEFANIDIS G D. Direct methane-to-ethylene conversion in a nanosecond pulsed discharge [J]. Fuel, 2018, 222: 705-710.

[26] GAO Y, ZHANG S, SUN H, et al. Highly efficient conversion of methane using microsecond and nanosecond pulsed spark discharges [J]. Applied Energy, 2018, 226: 534-545.

[27] DELIKONSTANTIS E, SCAPINELLO M, STEFANIDIS G D. Low energy cost conversion of methane to ethylene in a hybrid plasma-catalytic reactor system [J]. Fuel Processing Technology,2018, 176: 33-42.

[28] TU X, WHITEHEAD J C. Plasma dry reforming of methane in an atmospheric pressure AC gliding arc discharge:Co-generation of syngas and carbon nanomaterials [J]. International Journal of Hydrogen Energy, 2014, 39(18): 9658-9669.

[29] LIAO M, WANG Y, WU H, et al. Study of non-thermal DC arc plasma of CH4/Ar at atmospheric pressure using optical emission spectroscopy and mass spectrometry [J]. Plasma Science and Technology, 2015, 17(9): 743-748.

[30] ZHANG Hao, WANG Weizong, LI Xiaodong, et al. Plasma activation of methane for hydrogen production in a N2 rotating gliding arc warm plasma: A chemical kinetics study [J]. Chemical Engineering Journal, 2018, 345: 67-78.

[31] MARTINI L M, DILECCE G, GUELLA G, et al. Oxidation of CH4 by CO2 in a dielectric barrier discharge [J]. Chemical Physics Letters, 2014, 593: 55-60.

[32] 丁玉柱, 夏胜国，王琼芳，等. 介质阻挡放电对甲烷离解及燃烧火焰的影响[J]. 中国电机工程学报，2011，31 (31): 204-210.

DING Yuzhu, XIA Shengguo, WANG Qiongfang, et al. Effects of dielectric barrier discharge on methane dissociation and combustion flame [J]. Chinese Journal of Electrical Engineering, 2011, 31(31): 204-210.

[33] RAY D, REDDY P M K, CHALLAPALLI S. Glass beads packed DBD-plasma assisted dry reforming of methane [J]. Topics in Catalysis， 2017, 60(12/14): 869-878.

[34] DAI Bin, ZHANG X L, GONG W M, et al. Study on the methane coupling under pulse corona plasma by using CO2 as oxidant [J]. Plasma Science & Technology, 2000, 2(6): 577-580.

[35] ZHOU L M, XUE B, KOGELSCHATZ U, et al. Nonequilibrium plasma reforming of greenhouse gases to synthesis gas [J]. Energy & Fuels, 1998, 12(6): 1191-1199.

[36] 王晓玲，高远，张帅，等. 脉冲参数对介质阻挡放电等离子体CH4干重整特性影响的实验[J]. 电工技术学报，2019，34 (6): 1329-1337.

WANG Xiaoling, GAO Yuan, ZHANG Shuai, et al. Experiment on the effect of pulse parameters on the CH4 dry reforming characteristics of dielectric barrier discharge plasma[J]. Chinese Journal of Electrotechnical Technology, 2019, 34(6): 1329-1337.

[37] KOZLOV K V, MICHEL P, WAGNER H E. Synthesis of organic compounds from mixture of methane with carbon dioxide in dielectric-barrier discharge at atmospheric pressure [J]. Plasmas and Polymers, 2000, 5(3/4): 129-150.

[38] OZKAN A, DUFER T, et al. CO2-CH4 conversion and syngas formation at atmospheric pressure using a multi-electrode dielectric barrier discharge [J]. Journal of CO2 Utilization, 2015, 9: 74-81.

[39] WANG W Z, SNOECKX R, ZHANG X M, et al. Modeling plasma-based CO2 and CH4 conversion in mixtures with N2, O2, and H2O: The bigger plasma chemistry picture [J]. Journal of Physical Chemistry C, 2018, 122(16): 8704-8723.

[40] SUN J T, CHEN Q, GUO Y W, et al. Quantitative behavior of vibrational excitation in AC plasma assisted dry reforming of methane [J]. Journal of Energy Chemistry, 2020, 46: 133-143.

[41] SUN J T, CHEN Q. Kinetic roles of vibrational excitation in RF plasma assisted methane pyrolysis [J]. Journal of Energy Chemistry, 2019, 39: 188-197.

[42] HAGELAAR G J M, PITCHFORD L C. Solving the boltzmann equation to obtain electron transport coefficients and rate coefficients for fluid models [J]. Plasma Source Science and Technology, 2005, 14(4): 722-733.

[43] National Institute of Standards and Technology. NIST chemical kinetics database [EB/OL].[2022-06-28]. https://kinetics.nist.gov/kinetics/index.jsp.

[44] WANG, X, GAO S, ZHANG H, et al. Nanosecond pulsed pla-sma assisted dry reforming of CH4: The effect of plasma operating parameters [J]. Applied Energy, 2019, 243: 132-144.

[45] GOUJARD V, TATIBOUET J M, BATIOT-DUPEYRAT C. Use of a non-thermal plasma for the production of synthesis gas from biogas [J]. Applied A: General, 2009, 353: 228-235.

Mechanism exploration of non-thermal plasma reforming of CO2-CH4 bycombining emission spectroscopy and kinetics calculations

CHEN Yumin1,LU Qiancheng1,LEI Wentao1,YUE Wenjing1,ZHANG Bing1,HUANG Wenjie1,ZHAO Yongchun2,ZHANG Junying2

(1.School of Low-Carbon Energy and Power Engineering,China University of Mining and Technology,Xuzhou 221116,China;2.State Key Laboratory of Coal Combustion,Huazhong University of Science and Technology,Wuhan 430074,China)

Abstract:CO2-CH4 reforming by non-equilibrium plasma to syngas is an important emerging technology pathway to realize the resource utilization for the two greenhouse gases. The reaction pathway of CO2-CH4 reforming by non-thermal plasma in a dielectric barrier discharge (DBD) reactor at ambient temperature and pressure was thoroughly investigated by combing emission spectrometry analysis, reaction kinetics simulation and continuous mass spectrometry measurements. A maximum conversion of CH4 and CO2 of 25.8% and 9.6% respectively, can be achieved under a specific energy input (SEI) of 52 J/cm3. The continuous mass spectrometry online analysis of the stable gaseous products shows that C2H4 acts as a key intermediate in the formation of C2 hydrocarbons. The higher the proportion of CH4, the greater the generation of C2 hydrocarbons, and it increases in the order of C2H6<C2H4<C2H2.·CH, ·C2 radicals and characteristic spectral lines of CO and are detected in the emission spectrum. The relative intensity of ·CH decreases significantly with the increase of specific power, and the relative intensity of changes little. The generation and consumption pathways of C2 and C3 products are analyzed by kinetic calculation, and it is found that 89.2% of C2H6 comes from the combination reactionCH3 CH3( M)C2H6( M), 78.9% of C2H4 is mainly formed byCH4 CHC2H4 H and C2H3 H( M)C2H4( M) reaction, 61.1% of C2H2 is formed by electron collision reaction, and the formation of C2 product follows the path: CH4→C2H6→C2H4→C2H2. The reaction of electron collision dissociation e CO2CO O-of CO2 is an important formation pathway of CO, while the H2 and electron collision dissociation of CH4 are the key generation pathway of H2. Integration of emission spectroscopy and kinetic calculation provides a useful tool to uncover the reaction mechanism of plasma-enhanced CO2-CH4 reforming process.

Key words:non-thermal plasma;CO2-CH4 reforming;emission spectroscopy;reaction mechanism;chemical kinetics model;plasma chemistry

移动阅读

收稿日期:2022-06-30;责任编辑:白娅娜

DOI:10.13226/j.issn.1006-6772.HK22063003

基金项目：国家自然科学基金资助项目(51606216,52176144)；中央高校基本科研业务费专项资金资助项目(2020ZDPYMS23)；国家重大科研仪器研制资助项目(51827808)

作者简介:陈玉民(1986—)，男，湖南衡阳人，副教授，博士。E-mail:yuminchen@cumt.edu.cn

引用格式:陈玉民，卢钱程，雷文涛，等.联合发射光谱与动力学计算的非热等离子体重整CO2-CH4机理[J].洁净煤技术，2022,28(10):54-67.

CHEN Yumin,LU Qiancheng,LEI Wentao，et al.Mechanism exploration of non-thermal plasma reforming of CO2-CH4 by combining emission spectroscopy and kinetics calculations[J].Clean Coal Technology，2022,28(10):54-67.

中图分类号:TQ53；TK114

文献标志码:A

文章编号:1006-6772(2022)10-0054-14