大气压非平衡等离子体甲烷干法重整零维数值模拟*

钟旺燊 陈野力 钱沐杨† 刘三秋 张家良 王德真

1) (南昌大学物理系, 南昌 330031)

2) (大连理工大学物理学院, 大连 116024)

大气压非平衡等离子体由于其独特的非平衡特性, 可为甲烷和二氧化碳稳定温室气体分子活化和重整提供非热平衡和活化环境.本文采用了零维等离子体化学反应动力学模型, 考虑了详细的CH4/CO2 等离子体化学反应集, 重点研究了反应气体CH4/CO2 摩尔分数(5%—95%)对大气压非平衡等离子体甲烷干法重整制合成气和重要含氧化合物的影响.首先, 给出了进料气体不同体积比时电子密度和温度随时间的演化规律,结果表明初始甲烷摩尔分数的提高有利于获得较高的电子密度和电子温度.随后, 讨论了主要自由基和离子数密度在不同的甲烷摩尔分数下随着时间的变化规律, 并给出了反应气体的转化率、合成气体和重要含氧化合物的选择性.此外, 还明确了合成气和含氧化合物主要生成和损耗的化学反应路径, 发现甲基和羟基是合成含氧化合物的关键中间体.最后, 归纳总结给出了主要等离子体粒子之间的总体等离子体化学反应流程图.

1 引 言

全球能源需求的不断增长, 加速了化石燃料储量的枯竭, 同时还带来了气候变化的不利影响, 对人类生存环境构成了严重威胁, 引起了公众的广泛关注.在化石燃料占主导地位的世界, 天然气(主要成分为CH4)是继原油和煤炭之后的第三大能源, 越来越受到人们的重视.二氧化碳(CO2)是造成温室效应的主要因素, 也是碳的天然来源, 广泛存在于大气、沼气、油田伴生气中, 通过光合作用、燃料燃烧、有机分解等方式, 在自然生态平衡和现代社会运行中发挥着重要作用.CH4和CO2是两种主要的温室气体, 虽然CH4在大气中的浓度远低于CO2, 但其浓度增长率更快, 其温室效应也更具有长效性.另一方面, 合成气(CO 和H2)是环境友好的燃料和合成其他化学产品(如甲醇、乙醇,二甲醚等)的原料.目前, 工业上主要通过天然气的蒸汽重整反应(steam methane reforming, SMR)制备合成气[1]:

该反应是强烈地吸热和过分消耗甲烷的反应,容易造成宝贵的资源和能源的非高效利用.甲醇是结构最简单的饱和一元醇, 作为重要的基础有机化工原料, 主要用于塑料、精细化学品和石油化工等领域[2].甲醇同时也是一种新型的清洁能源, 得益于其良好的抗爆性和高辛烷值, 可以用作车用发动机燃料.目前工业上主要通过将煤转化成含有一氧化碳、二氧化碳和氢气的合成气, 以铜基催化剂高压(5—10 MPa)和高温(493—593 K)转化合成得到甲醇[3].此反应条件苛刻, 也会造成环境污染和消费大量的能源.因此, 利用CH4和CO2的重整反应制取合成气和高价值的液态含氧化合物(如CH3OH 和CH2O 等), 对缓解能源危机、减轻温室气体的排放而导致的全球气候变暖等环境问题具有重要的战略意义.在最近的国际催化会议上, 世界各地的化学家们又重提了甲烷和二氧化碳重组制合成气, 并且在这一研究领域做了大量的研究工作, 取得了一系列研究成果[4−8].一般通过甲烷干重整(dry reforming methane, DRM), 实现甲烷和二氧化碳重组制合成气.目前的研究工作主要集中于尝试各式催化剂和反应条件等, 对干法重整过程中反应气体转化率, 合成气选择性和产率的影响, 以此评估整体的催化性能和合成气形成机理.DRM 工艺最大的问题之一是催化剂材料处理有大量的碳沉积, 引起催化剂中毒, 这可能是DRM 过程中遇到的一个难于解决的大问题.此外, CH4和CO2都是化学性质非常稳定的分子, 其中的C—O和C—H 化学键能分别为532.2 kJ/mol 和439.3 kJ/mol[9].因此, 不管是通过甲烷的SMR 还是DRM反应制取合成气, 其活化条件都比较苛刻, 需要消耗大量的能量和反应废物, 整体的能效并不高.

大气压非平衡等离子体由于其独特的非平衡特性, 即电子温度一般为0—10 eV (1 eV ≈11605 K ≈ 96.32 kJ/mol), 而重粒子的温度仅为几百K, 电子温度远高于重粒子温度.高能电子的能量高于温室气体分子的化学键能, 能够扮演化学“催化剂”的角色, 碰撞解离温室气体分子, 从而实现常温常压下温室气体分子的活化.但不得不指出, 虽然大气压非平衡等离子体能够为CH4和CO2分子活化提供独特的非热平衡和活化环境, 但也存在反应物转化率低, 产物种类复杂, 目标产物选择性低等缺点.能量利用率和副产物是制约大气压非平衡等离子体催化工业化的两个最主要因素.这是因为在等离子体中的化学反应可控性较差, 导致最终目标产物的选择性和产率不够理想.目前, 学术界主流观点认为, 可通过大气压非平衡等离子和利于目标产物生成的化学催化剂组合, 实现大气压非平衡等离子体催化的大规模生产使用.

目前, 大气压非平衡等离子体作为一种新兴的放电技术, 作为“反应载体”用以实现二氧化碳和甲烷重组, 受到越来越多的科研人员关注, 主要通过实验诊断[10−17]和数值模拟[18−24]对甲烷放电和DRM 反应机理进行了详细的研究.例如2020 年,Zhang 等[10]通过发射光谱和化学气相色谱等表征手段, 研究了电子碰撞等离子体化学反应和温度热效应化学反应, 在温度可控的大气压介质阻挡反应腔中甲醇和二氧化碳重组制取合成气所扮演的角色.Brune 等[11]研究了介质阻挡放电(dielectric barrier discharge, DBD)反应器和催化剂对大气压甲烷DRM 过程的影响, 研究表明填充了多孔γ-Al2O3载体填充式DBD 反应器, 能够改变放电模式和化学活性, 从而实现CH4/CO2协同催化重组.此外, 催化剂为Cu 时, CH4的转化率高达90%,这是由于Cu 催化剂对于O 原子和CO 分子有更好的吸附和解吸附能力.Maqueo 等[12]研究了反应气体能量输入, 单脉冲能量和脉冲频率对大气压纳秒脉冲放电甲烷重组的影响, 实验结果论证了大气压纳秒脉冲反应器可作为燃料的重整腔室.脉冲频率是最为重要工作参数, 增加脉冲频率有助于提高反应气体的转化率和能量效益.Alawi 等[13]研究了掺杂气体N2和Ar 流量, 微波功率和CO2/CH4体积比, 对大气压低功率微波等离子体辅助DRM的影响.实验结果表明, 在N2和Ar 氛围内的DRM过程有着相同的反应性能.在相同的工作参数下,如恒定微波功率700 W, CO2/CH4体积比为2∶1,掺杂气体(N2或Ar), CO2, CH4的流量分别为1.5, 0.4 和0.2 L/min, Ar 氛围内的DRM 能产生更高的合成气体选择性和产率.王晓玲等[14]通过详细的电学和化学气相色谱, 对微秒脉冲电源驱动的DBD 甲烷干重整反应过程进行了诊断, 主要研究了脉冲DBD 放电特性及重复频率和脉宽对气体转化特性的影响.研究结果表明, 增大脉冲重复频率有利于提高CH4和CO2的转化率, 同时获得较高的合成气产率.Wu 等[15]利用耦合了磁场和切线流的旋转滑动弧等离子体, 实现了等离子体辅助的DRM, 主要研究了电压和CH4/CO2体积比对等离子体基本放电参数(电子激发温度, 电子密度和转动温度)及放电产物分布的影响.在论文综述中, 详细地介绍了DBD 等离子体DRM的研究现状、取得的研究进展以及遇到的研究瓶颈[16].现有的研究主要集中于DBD 反应结构, 如电极形貌、放电间隙、放电体积及填充的催化剂对DRM 过程的影响.研究结果表明, 反应气体的空间流动速率、能量密度以及反应气体的体积比对放电产物的分布起着至关重要的作用.

在数值模拟方面, 前人也做了很多有意义的研究工作.赵曰峰等[18]建立了二维大气压甲烷针-板放电流体模型, 并进行了相关实验验证, 研究得到了主要粒子浓度的时空分布和反应路径.Slaets 等[19]采用零维(zero-dimensional, 0D)化学反应动力学模型研究了滑动弧等离子管中, 掺入少量的氮气和氧气对大气压甲烷的干重整的影响.研究结果表明, 氮气有助于提高二氧化碳的转化率, 但是部分输入能量会耗散在氮气分子转动激发上.此外, 氧气有助于提高甲烷的转化率, 降低能量损耗, 但二氧化碳的转化率基本保持不变.王伟宗等[20]总结出一套详细完备的等离子体化学反应集, 可实现CH4, CO2, N2, O2和H2O 等反应气体任意组合所涉及的等离子体化学反应模拟, 从而构建出这一研究领域中, 大气压DBD 等离子体更为宽广的等离子体化学图像.Snoeck 等[21]采用了零维化学反应动力学模型, 研究了等离子体放电脉冲及其余辉中的等离子体行为和化学反应机理, 计算给出了反应气体的转化率、高价值含氧化合物的选择性和产率, 以及DBD 等离子体的能量效益.Liu 等[22]对等离子体催化CH4/CO2转化成三种不同的含氧化合物(CH3OH, HCHO 和CH3COOH) 潜在的反应机理, 从ns-ms 和nm-mm 的时空尺度进行了详细的总结.比利时安特卫普大学Bogaerts 教授课题组[23,24]分别采用0D 和1D 等离子体化学反应动力学模型, 模拟了大气压CH4/CO2DBD 等离子体复杂的等离子体化学转化反应机理, 计算给出了DBD 活性粒子数密度, 放电产物的转化率、产率和反应气体的转化率, 及整体反应流程图等.

综上所述, 目前大气压非平衡等离子体DRM的研究主要通过实验诊断(发射光谱、电学特性、化学气相色谱等)和流体模型数值模拟, 研究DBD反应器工作参数, 如CO2/CH4体积比, 反应停留时间, 催化剂类型和填充方式, 掺杂气体, 电源工作参数, 电极结构等对等离子体DRM 过程中产物分布的影响, 从而分析出DRM 过程中潜在的反应机理.前人的研究结果表明, 初始甲烷摩尔分数是等离子体DRM 过程的一个很重要的影响参数.本文采用零维等离子体化学反应动力学模型, 研究了初始甲烷摩尔分数对等离子体DRM 过程的影响,重点关注合成气(CO 和H2)以及高价值的含氧化合物(CH2O 和CH3OH)主导的生成和损耗反应路径, 反应气体的转化率和重要产物的选择性.最后总结归纳出了主要等离子体粒子之间的整体化学反应流程图, 直观地展示出等离子体粒子之间的化学转化关系.

2 模型描述

本文利用0D 等离子体化学反应动力学模型[25−29], 考虑了详细的CH4和CO2等离子体化学反应集, 以描述大气压非平衡等离子体DRM 潜在的化学反应机理.该模型主要用于求解包括非平衡等离子体中各类物质数密度随时间变化在内的复杂气相等离子体化学反应机理.模型中的等离子体被视为“等离子体化学反应器”, 从而忽略等离子体粒子的输运, 认为粒子具有均匀的数密度.该模型主要由三部分组成, 即等离子体化学反应集模块、玻尔兹曼求解器模块和电路模块.在本文中, 只使用了前面两个模块.更具体地, 该模型由包含化学反应集与用户自定义等离子体放电特性仿真设置的输入文件, 求解电子碰撞反应速率系数的玻尔兹曼求解器, 和求解微分方程的DVODE_F90 组成[30].基于粒子的产生和损耗反应过程, 粒子数密度随时间的演化可由等离子体化学反应计算如下:

其中aij(1)和aij(2)分别表示物种i在反应j左侧和右侧的化学计量系数;ni是物种i在反应左侧的数密度;kj是反应j的速率系数.重粒子之间碰撞反应的速率系数kj通过标准Arrhenius 公式k=ATBexp(–C/T)计算得到, 其中T是背景气体的温度, 本文设定为常数300 K.电子碰撞反应的速率系数kj是通过将LXCAT 网站上获取的电子碰撞截面数据导入到内置的玻尔兹曼求解器模块中计算得出.电子温度取决于焦耳加热和电子碰撞反应, 可以用下面的能量平衡方程来计算:

其中ne是电子密度;Te是电子温度;j和E分别是等离子体中的电流密度和电场;vmi是物质i的电子动量传递碰撞频率;me是电子质量;Mi和Ti分别是物种i的质量和温度.kl是电子撞击反应l的速率系数;Nl是参与反应的重物质的数密度;Δεl是电子能量的变化.(3)式右边的第一项是由施加的功率密度脉冲引起的能量增益, 第二项和第三项分别表示由弹性和非弹性碰撞引起的能量损失.上述微分方程使用DVODE 进行积分, 焦耳热用功率密度表示.在该模型中, 等离子体的总体积设置为0.08 cm3, 与等离子体接触的反应器壁的表面积为1.68 cm2, 等离子体的横向扩散长度为0.0637 cm, 反应器壁的温度为300 K, 初始电子温度为0.025 eV.

本模型中考虑的等离子体化学反应集, 主要参考了比利时安特卫普大学Bogaerts 教授课题组[21,23,24]相关的大气压非平衡CH4/CO2等离子体零维模型数值模拟工作.在此基础上, 增加了15 个C 和H相关的中性粒子碰撞反应[29,31].最终本模型所考虑的CH4/CO2等离子体化学反应集, 包含82 种等离子体粒子和885 个等离子体化学反应(60 个电子碰撞反应, 334 个离子反应和491 个中性粒子反应).

3 仿真结果及讨论

3.1 电子时间演化特性

对于大多数反应气体, 大气压DBD 通常处在所谓的丝状状态, 包含大量的独立微放电丝状通道, 其持续时间为纳秒量级, 通常为1—100 ns[21].在微放电丝状通道中, 电子的能量通过电子直接碰撞的形式, 用于反应气体分子的激发、电离和解离等反应过程, 从而产生大量的电子、离子、激发态粒子和中性自由基等.这也就是这些丝状放电对于流体模型描述化学反应至关重要的原因.在本文的0D 模型中, 并不能刻画描述丝状放电的空间演化过程, 但可以描述微放电脉冲及其余辉的等离子体反应过程.因此, 施加到等离子体反应器的为三角脉冲, 其持续时间为30 ns(约为一个丝状放电持续时间), 图1 虚线是脉冲功率的波形图, 其峰值功率密度为6.0 × 105W/cm3.0 ns 时脉冲功率开始线性增加, 15 ns 时达到峰值功率, 在随后的15 ns时间内线性下降至0, 因此三角形脉冲功率的持续时间为30 ns.一个脉冲周期中余辉期间为60 ns,即脉冲的周期为90 ns.图1 实线给出的是甲烷摩尔分数为10%, 30%, 50%, 70%和90%时, 在五个周期内电子温度随时间的演化曲线, 从图1 中可以清楚地看出, 电子温度随着脉冲功率增加而线性快速上升, 然后上升速率略微变缓, 在脉冲功率未达到最大值时, 电子温度已达到峰值约为2.0 eV, 这主要是由于电子快速响应电场加速获得的能量.随后, 在脉冲功率达到最大值之前, 电子温度又快速线性下降至0 左右, 并在余辉中也保持在0 附近.我们也注意到, 在参考文献[21]中, 也观察到类似的电子温度随时间的演化特性.数值模拟结果表明, 甲烷的摩尔分数越大, 电子温度的增加速率越快, 电子温度的峰值也越高.例如, 甲烷摩尔分数为10%, 30%, 50%, 70%和90%时, 电子温度的峰值分别为1.85, 1.98, 2.18, 2.38 和2.55 eV.电子密度随时间演化特性与电子温度不尽相同.如图2所示, 随着脉冲功率的线性增加并未达到峰值前,电子密度也线性增加至峰值(1014—1015cm–3), 随后快速小幅度地下降, 再随着时间缓慢下降.不同甲烷的摩尔分数时, 电子密度随时间的整体演化规律相类似, 但较高的甲烷含量放电产生的电子密度更高.De Bie 等[24]的研究也发现, 在CH4/CO2反应气体等离子体体系中, 较高的甲烷占比能使得电子密度维持在一个更高的水平.

图1 甲烷摩尔分数为10%, 30%, 50%, 70%和90%时的电子温度随时间演化规律Fig.1.Electron temperature as a function of time for methane mole fractions of 10%, 30%, 50%, 70% and 90%.

图2 甲烷摩尔分数为10%, 30%, 50%, 70%和90%时的电子密度随时间变化趋势Fig.2.Electron density as a function of time for methane mole fractions of 10%, 30%, 50%, 70% and 90%.

3.2 粒子时间演化特性

图3 展示了不同甲烷含量时主要自由基数密度随时间的变化趋势以及自由基值周期平均值随甲烷摩尔分数的变化规律.如图3 所示, 大气压非平衡等离子体DRM 的主要自由基为H, CH2, CH3,O, OH, C, CH3O 和CH2OH 等, 其中H 原子的密度最高维持在1015—1016cm–3数量级, 其次是CH2,CH3等.随着甲烷含量的增加, H, CH2, CH3自由基密度会随着增加, 但由于CO2含量在减小, 因而缺乏足够的O 原子参与反应, 导致最终生成的含氧化合物的密度会降低, 如图3(d)所示.例如甲烷摩尔分数为10%时, H 原子的数密度达到了1015cm–3的数量级, O 原子的数密度达到1014cm–3的数量级, 而当甲烷摩尔分数增加到90%时, O 原子的数密度显著降低为1013cm–3的数量级, H 原子的数密度小幅度增加接近1016cm–3的数量级.这主要是由于O 原子和H 原子只能分别来源于进料气体CO2和CH4, CH4含量的增加(同时CO2的含量在减小), 解离产生的O 原子更少, 从而生成的含氧化合物密度会降低.从图3(a)—(c)中可看出,CH3自由基的数密度随时间的变化呈现出较为明显周期性波动, 上下波动范围可达 1 个数量级.值得一提的是, 图3 中所展示的只是仿真起始一段脉冲周期内(2 µs, 约22 个周期)的中性自由基数密度随时间的演化趋势, 这段起始时间内CH3自由基的密度比CH2密度略低, 但经过一段反应时间后, CH3自由基的密度会高于CH2自由基.一方面, CH3和CH2自由基的主要生成反应路径均是电子和CH4分子的直接碰撞解离反应, 即E+CH4→E + CH3+ H 和E+ CH4→E + CH2+ H2,但CH3自由基主要生成反应的反应速率比CH2自由基主要生成反应的反应速率至少高了1 个数量级(数据未给出).另一方面, 本模型中三角形脉冲周期内余辉持续时间仅为60 ns, 导致起始一段周期内CH3自由基的密度会略低于CH2自由基, 且CH3自由基随时间的变化呈现出较为明显的周期性波动, 而CH2自由基未能及时呈现周期性波动.

图3 (a) 10%, (b) 50%和(c) 90%甲烷摩尔分数时主要自由基的数密度随时间变化趋势, 以及(d) 主要自由基的周期平均值随甲烷摩尔分数的变化Fig.3.The number densities of main radicals as a function of time for methane mole fractions of (a) 10%, (b) 50%, (c) 90%, and(d) time averaged number densities of main radicals as a function of initial CH4 fraction.

图4 展示了不同甲烷含量时主要离子的数密度随时间的变化规律以及离子周期平均值随甲烷摩尔分数的变化规律.如图4 所示, 主要离子为, OH–, H–, O–,等, 其中阳离子的密度最高维持在1014cm–3数量级, 其次是OH–, H–等阴离子.值得一提的是, 在甲烷含量较低时, H–离子密度整体低于1012cm–3, 且周期性波动的振幅较大, 而甲烷含量较高时, H–离子密度整体接近1012cm–3, 且周期性波动的振幅较小.这些主要离子随时间的演化趋势均呈现幅度大小不一的周期性波动, 幅度较大的是(幅度跨越8 个数量级), 其次是O–和, 较小的是和OH–.周期性波动幅度较大的离子主要来源于CH4和CO2分子的电子碰撞电离反应.对于阴离子O–也呈现幅度较大的周期性涨落, 是因为其主要的生成反应路径为E+ CO2→CO + O–.因此, 在本模型中,是最高浓度的阳离子, OH–是最高浓度的阴离子, 这与Bai 等[32]的结果相一致.此外, 如图4(d)所示, 当甲烷含量增加时, 除了H–, 其他的离子均呈现下降趋势.相比于主要自由基数密度, 主要离子的数密度整体要低2 至3 个数量级, 这表明离子在等离子体化学中的作用并不大.对于含氧的阴离子OH–, O–, 和的浓度下降可能是由于起始CH4含量的增加, 从而CO2含量的较少所导致的.图5 展示了不同甲烷含量时主要分子的数密度随时间的变化规律以及分子周期平均值随甲烷摩尔分数的变化规律.大气压非平衡等离子体DRM 主要 分 子 为CO, H2, H2O, CH2O, O2, CH3OH 等,其中合成气是主要的分子产物, 而重要环氧化合物CH2O 和CH3OH 的浓度较低.如图5(d)所示,随着CH4含量的增加, 仅H2密度有较为明显的增加, 而其他分子的密度均出现不同程度的下降, 其中O2受到的影响最大, 从1014cm–3下降到109cm–3数量级.这主要是由于CO2含量的减小, 解离产生的O 和O2均会很明显地减小, 因而缺乏足够的O 原子参与其他含氧化合物.

图4 (a) 10%, (b) 50%和(c) 90%甲烷摩尔分数时主要离子的数密度随时间的变化趋势, 以及(d)主要离子的周期平均值随甲烷摩尔分数的变化Fig.4.The number densities of main ions as a function of time for methane mole fractions of (a) 10%, (b) 50%, (c) 90%, and(d) time averaged densities of main ions as a function of initial CH4 fraction.

图5 (a) 10%, (b) 50%和(c) 90%甲烷摩尔分数时主要分子数密度随时间的变化规律, 及(d)主要分子的周期平均值随甲烷摩尔分数的变化Fig.5.The number densities of main molecules as a function of time for methane mole fractions of (a) 10%, (b) 50%, (c) 90%, and(d) time averaged densities of main molecules as a function of initial CH4 fraction.

3.3 反应气体转化率和主要产物选择性

反应气体的转化率X, 和主要产物选择性S的计算公式分别为:

图6(a)展示了甲烷和二氧化碳的转化率随甲烷摩尔分数的变化趋势.显然, CO2转化率具有更大的变化范围, 在CH4的摩尔分数超过70%时CO2的转化率就会高于CH4的转化率, 且随着CH4的摩尔分数的进一步增加, CO2转化率远高于CH4转化率.具体地, 当甲烷的摩擦分数从5%增加到95%时, 甲烷的转化率从8.0%缓慢下降至3.0%左右, 而二氧化碳的转化率从0.69%逐渐上升至27%.我们也注意到, Zhang 等[10]对于大气压等离子体DRM 过程中电子碰撞化学反应和背景温度有关的热化学反应所扮演的角色, 开展了详细的实验研究.当DBD 反应器温度控制在293 K 且功率固定为10 W 时, 实验结果表明随着甲烷的摩尔分数从20%增加到66.7%, 甲烷和二氧化碳的转化率分别维持在7.5%和3.3%左右.且在这个范围内甲烷含量的变化, 对反应气体二氧化碳和甲烷的转化率影响并不大.在本模型中, 我们的等离子体气体温度固定为300 K, 如图6(a)所示, 当甲烷的摩尔分数从20%增加到66.7%时, 甲烷和二氧化碳的转化率均不高且变化不太大, 但甲烷的转化率比二氧化碳转化率更高, 而且整体上和实验数据符合的较好.反应研究详细了合成气(H2和CO)是DRM 过程中形成的主要产物, 主要关注了两种最重要的含氧化合物CH3OH 和CH2O, 图6(b)中给出了这两种含氧化合物的选择性随着CH4摩尔分数的变化规律.如图6(b)所示, 合成气的选择性明显高于含氧化合物的选择性, CH2O 的选择性整体低于5%, 而CH3OH 的选择性更低(低于1%).随着CH4摩尔分数从5%增加到95%, 只有H2的选择性在整体上表现为从13.0%增加到24.6%,而CO 的选择性从58.9%下降至9.7%.涉及到的主要物理机理为, 氢原子来自于甲烷, 高甲烷含量有利于H2的形成, 且不同的CH4摩尔分数时, H2的生成都主要来源于电子和CH4的碰撞解离反应(如图10 所示).CH2O 的选择性相比于CH3OH 更高, CH2O 和CH3OH 的选择性均随着CH4的摩尔分数增加呈现为缓慢降低的趋势.这主要是因为,如后文图8 和9 所示, CH2和CO2以及CH3和O之间的反应是主要的CH2O 生成路径, 而CH3和OH 的三体复合反应是CH3OH 的主要生产路径.另一方面, 如前文图3 所示, 随着甲烷摩尔分数的增加, CH2, CH3, OH 的密度并没有显著升高.因此, 最终两种重要的含氧化合物CH2O 和CH3OH的选择性均随着CH4摩尔分数的增加呈现为缓慢降低的趋势.

图6 (a)进料气体的转化率和(b)合成气和重要含氧化合物的选择性随着甲烷摩尔分数的变化趋势Fig.6.Time-averaged (a) conversion, (b) selectivity as a function of initial CH4 mole fraction.

3.4 主导反应路径

为了进一步探究合成气和重要含氧化合物的生成和损耗机理, 图7—图10 依次给出了CH2O,CH3OH, CO, H2主要的生成和损耗反应的时间平均反应速率随甲烷摩尔分数(10%, 50%, 90%)的柱状图.如图7 所示, 当甲烷摩尔分数较低时, CH2和CO2的反应是主要的CH2O 生成反应路径, 而当甲烷摩尔分数增加时, CH3和O 之间的反应变成了主要的CH2O 生成反应路径.CH2O 的主要损耗反应路径是CH2O 分别和O 和OH 生成CHO的反应.如图8 所示, CH3OH 主要的生成化学反应路径是CH3, OH 和CO2的三体复合反应(CH3+OH+CO2→CH3OH+CO2), 这是因为这两种自由基的高浓度和活性, 各自孤立的电子恰好形成一对共价键.值得注意的是, 随着CH4摩尔分数增加为90%时, CH3, OH 和CH4的三体复合反应速率比OH–和CH3之间的反应速率更大, 成为CH3OH分子生成的次要反应.这是因为甲烷摩尔分数从10%增加到90%时, 一方面OH–密度下降了至少一个量级(如图4(c)), 另一方面是CH4的含量远大于CO2.CH3OH 的主要损耗反应路径是CH3OH分子分别与O 和OH 反应生成CH2OH 的反应.此外, O 原子可分别从CH3OH 的甲基(-CH3)或者羟基(-OH)中夺取H 原子分别生成CH2OH 和CH3O.数值模拟结果表明, O 和CH3OH 反应生成CH2OH 的反应速率远大于O 和CH3OH 反应生成CH3O 的反应速率, 这表明CH3OH 中的甲基比羟基更容易失去氢原子.

如图9 所示, CO 主要来源于CO2的电子碰撞反应, 该电子碰撞反应在所有的CO 生成反应中反应速率最大.当甲烷含量较低为10%时,E+ CO2→O–+ CO 的反应速率比其他CO 生成反应大了近三个数量级.而当甲烷摩尔分数逐渐升高时,E+CO2→O–+CO 的反应速率逐渐减小, 而其他的CO 生成反应的反应速率逐渐增加, 如C+O2→O +CO,E+ CO2+→O + CO, 但E+ CO2→O–+CO始终是CO 的主导生成反应路径.电子和CO 碰撞解离和电离反应, 是CO 的主要消耗反应路径.图10给出的是H2的主要生成和损耗反应路径.结果表明,CH4的电子碰撞解离反应(E+ CH4→E + CH2+H2,E+ CH4→E + CH + H + H2)是H2分子最重要的生成反应路径.而电子和H2碰撞电离反应则是H2主要的损耗反应路径.如图9 和图10 所示, H2主要生成反应的反应速率远大于CO 的生成反应速率, 而且H2主要损耗反应的反应速率却小于CO 的主要损耗反应速率, 这也就是生成物中H2的密度高于CO 的原因(如图5 所示).

图7 CH2O 的主要生成和损耗反应的时间平均反应速率随甲烷摩尔分数的变化柱状图 (a) 10%, (b) 50%, (c) 90%Fig.7.Time-averaged reaction rates of the dominant reaction pathways for the production and consumption of CH2O as a function of methane mole fraction: (a) 10%, (b) 50%,(c) 90%.

图8 CH3OH 的主要生成和损耗反应的时间平均反应速率随甲烷摩尔分数的变化柱状图 (a) 10%, (b) 50%, (c)90%Fig.8.Time-averaged reaction rates of the dominant reaction pathways for the production and consumption of CH3OH as a function of methane mole fraction: (a) 10%,(b) 50%, (c) 90%.

图9 CO 的主要生成和损耗反应的时间平均反应速率随甲烷摩尔分数的变化柱状图 (a) 10%, (b) 50%, (c) 90%Fig.9.Time-averaged reaction rates of the dominant reaction pathways for the production and consumption of CO as a function of methane mole fraction (a) 10%, (b) 50%,(c) 90%.

图10 H2 的主要生成和损耗反应的时间平均反应速率随甲烷摩尔分数的变化柱状图 (a) 10%, (b) 50%, (c) 90%Fig.10.Time-averaged reaction rates of the dominant reaction pathways for the production and consumption of H2 as a function of methane mole fraction (a) 10%, (b) 50%,(c) 90%.

图11 总结给出了CH4/CO2摩尔分数比为1∶1 时大气压非平衡等离子体DRM(CH4/CO2=1∶1)生成合成气和CH2O, CH3OH 等含氧化合物的总体反应流程图.为了流程图的清晰和紧凑起见, 仅考虑时间平均的净反应速率超过1018cm3/s的反应路径, 且图中箭头的宽度与时间平均的净反应速率成线性正比.电子和甲烷分子的碰撞解离反应生成CH3, CH2, H2, 其中生成CH3自由基的净反应速率最大.形成的甲基中, 一部分CH3自由基将重新结合成高级烃(C2H5, C2H6和C3H8等), 一部分参与反应形成含氧化合物如CH2O 和CH3OH等, 还有一部分反应生成H2.含氧化合物中的O来源于CO2分子, 主要通过电子和CO2分子的碰撞解离形成CO, O, O2等, 其中绝大部分形成了CO, CO 进一步生成O 原子, 因此O 原子主要是由CO2解离形成的.O 原子进一步转移到OH 等粒子, 并与CH3反应参与形成CH2O 与CH3OH等含氧化合物.此外, CO2也可和甲烷解离产物CH2反应形成CH2O.如图11 所示, 反应气体CH4/CO2的电子碰撞解离反应生成CO 和H2的净反应速率明显高于其他含氧化合物, 此外甲基和羟基对于含氧化合物CH2O 与CH3OH 的生成至关重要.总体反应流程图进一步揭示了大气压非平衡等离子体DRM 的反应机理, 为大气压非平衡等离子体DRM 提供了理论参考.

图11 CH4/CO2 摩尔分数比为1:1 的大气压非平衡等离子体DRM 反应总体流程图.箭头的粗细与时间平均的净反应速率成线性正比Fig.11.Schematic overview of the dominant reaction pathways for the conversion of CH4 and CO2 into representative higher oxygenates and syngas in atmospheric non-equilibrium plasma for a 1:1 CH4/CO2 gas mixture.The thickness of the arrows is linearly proportional to time-averaged rate of net reaction.

4 结 论

本文以温室气体CH4/CO2为反应气体构建了大气压非平衡等离子体DRM 零维等离子体化学反应动力学模型, 重点研究了进料气体CH4/CO2摩尔分数对DRM 制合成气(CO, H2)和含氧化合物(CH2O, CH3OH 等)的影响.数值模拟结果表明, 电子温度和电子密度随施加的三角形脉冲功率密度呈现周期性变化, 且较高的甲烷摩尔分数有利于放电产生更高的电子密度和电子温度.放电产生的主要自由基和分子为H, CO, CH2, CH3, O, H2等, 随着甲烷摩尔分数的增加, H, H2, CH3的密度会随着增加, 但由于二氧化碳摩尔分数的减小, 因而缺乏足够的氧原子参与反应, 导致含氧化合物的密度会降低.主要的离子为, OH–等, 随着甲烷摩尔分数的增加, 除H–阴离子外, 其他离子密度均会降低.二氧化碳的转化率随着甲烷摩尔分数的增加而增加, 但甲烷的转化率对甲烷含量的变化并不太敏感, 在甲烷摩尔分数较低的时候甲烷的转化率较高.含氧化合物的选择性均较低(低于5%),甲烷摩尔分数从5%增加至95%时, 只有H2的选择性在整体上表现为从13.0%增加到24.6%, 而CO 的选择性从58.9%下降至9.7%.此外, 甲基和羟基的三体复合反应是甲醇的主要生产反应, 电子和甲烷分子的碰撞解离是氢气的主导生成反应, 电子和CO2分子的碰撞反应(E+ CO2→O–+ CO)是CO 的主要生成反应.对于CH2O 分子, 当甲烷摩尔分数较低时, CH2和CO2的反应是主要的生成反应, 但当甲烷摩尔分数较高时, CH3和O 的反应逐渐占据了主导位置.最后, 归纳总结给出了大气压非平衡等离子体DRM 整体的反应流程图, 结果表明甲基和羟基对于含氧化合物CH2O 与CH3OH等的生成至关重要.