劉倩，鄭洪濤，楊仁，陳曦

(哈爾濱工程大學(xué)動(dòng)力與能源工程學(xué)院，哈爾濱黑龍江150001)

化學(xué)回?zé)嵫h(huán)燃?xì)廨啓C(jī)(CRGT)作為一種先進(jìn)循環(huán)燃?xì)廨啓C(jī)，其熱效率由簡(jiǎn)單燃機(jī)的35.04%～40.1%增加到46.5% ～51.3%，NOx等污染物排放由幾百微降低到十幾微甚至1×10-6以下［1-2］?；瘜W(xué)回?zé)崞髯鳛镃RGT的核心部件，通過強(qiáng)吸熱的燃料蒸汽重整方式增加余熱利用深度。目前，工業(yè)用甲烷蒸汽重整反應(yīng)需要在高溫(700～900℃)高壓(2～3 MPa)的條件下在催化劑表面進(jìn)行，為了和CRGT渦輪排氣余熱溫度(350～500℃)相吻合，選擇一種新的催化方式顯得尤為重要。與傳統(tǒng)的熱化學(xué)方法相比，等離子體輔助甲烷蒸汽重整具有反應(yīng)速度快、反應(yīng)溫度低、裝置體積小、啟動(dòng)快等優(yōu)點(diǎn)［3-4］。因此作者提出了采用介質(zhì)阻擋放電(DBD)低溫等離子體催化燃料蒸汽重整技術(shù)用于化學(xué)回?zé)崞髟O(shè)計(jì)的新思路。由于DBD具有明顯的非平衡特性，為了提高目的產(chǎn)物的收率以提高能源利用率，需要對(duì)DBD催化蒸汽重整的反應(yīng)機(jī)理進(jìn)行研究。目前國(guó)內(nèi)外學(xué)者在數(shù)值研究方面有了一定進(jìn)展，主要建立了純流體加熱模型［5］、Fluid 模型［6］、0維替代模型［7］、湯生(Townsend)電離模型［8］等。上述學(xué)者多對(duì)DBD催化甲烷偶聯(lián)、部分氧化或二氧化碳重整等進(jìn)行了數(shù)值研究，而對(duì)DBD催化甲烷蒸汽重整的數(shù)值研究較少。

考慮到水蒸氣的來源問題，目前學(xué)者通過實(shí)驗(yàn)主要研究了水浴溫度(120℃左右)下等離子體催化甲烷蒸汽重整特性［9-10］。而在化學(xué)回?zé)嵫h(huán)燃?xì)廨啓C(jī)中，蒸汽通過蒸汽發(fā)生裝置產(chǎn)生，含量充足，因此研究高蒸汽流量下的等離子體催化甲烷蒸汽重整情況顯得十分必要。首先通過實(shí)驗(yàn)的方法研究較高溫度(350～500℃)時(shí)不同操作參數(shù)對(duì)重整過程中反應(yīng)物轉(zhuǎn)化及其產(chǎn)物選擇性的影響規(guī)律，為反應(yīng)機(jī)理的研究提供數(shù)據(jù)支持。然后發(fā)展和建立一套可用于描述DBD催化甲烷蒸汽重整的詳細(xì)反應(yīng)機(jī)理方程，并結(jié)合敏感性分析及其實(shí)驗(yàn)數(shù)據(jù)對(duì)其進(jìn)行修正和完善。再結(jié)合該動(dòng)力學(xué)機(jī)理模型和CHEMKIN中路徑分析模塊研究低溫常壓無催化劑時(shí)單獨(dú)等離子體催化甲烷蒸汽重整技術(shù)對(duì)甲烷轉(zhuǎn)化的動(dòng)力學(xué)影響。最后通過在放電區(qū)域添加Ni催化劑的實(shí)驗(yàn)即等離子協(xié)同催化實(shí)驗(yàn)對(duì)上述分析進(jìn)行了初步的驗(yàn)證。

1 實(shí)驗(yàn)裝置

實(shí)驗(yàn)裝置見圖1。反應(yīng)物為甲烷和水蒸氣，其中水蒸氣由蒸汽爐產(chǎn)生，同時(shí)在入口反應(yīng)物中加入體積分?jǐn)?shù)為10%的氮?dú)?，用于定量?jì)算出口處氣體總摩爾流量。甲烷和氮?dú)饬髁客ㄟ^質(zhì)量流量控制器(Brooks 5850E)控制，水流量通過平流泵(Beijing Xingda Co.，2PB00C)控制。氣體摻混均勻后進(jìn)入反應(yīng)器。DBD反應(yīng)器采用石英玻璃作為阻擋介質(zhì)，其放電區(qū)域長(zhǎng)130 mm，放電間隙為1 mm，總放電區(qū)域體積為12 cm3，見圖2。反應(yīng)器采用陶瓷加熱器進(jìn)行加熱，并利用溫度控制器控制其加熱溫度，保證溫度穩(wěn)定在設(shè)定值。由于熱慣性，實(shí)際溫度值與設(shè)定溫度值有5℃左右的偏差。反應(yīng)器輸入電壓和頻率是由低溫等離子體電源(CTP-2000K)控制，其放電電壓為4 000 V，頻率為41.93 kHz，總有效輸入功率為70 W。在每組實(shí)驗(yàn)的開始和結(jié)尾均記錄電壓和頻率值以確保實(shí)驗(yàn)的穩(wěn)定運(yùn)行。反應(yīng)器出口處氣體利用氣相色譜儀(Agilent GC 7890A)進(jìn)行定性和定量分析，其中CH4通過FID檢測(cè)器進(jìn)行分析，CO、CO2和 N2通過 TCD B檢測(cè)器進(jìn)行分析，H2通過TCD C檢測(cè)器進(jìn)行分析。為了保證色譜儀的正常工作，出口處未參加反應(yīng)的水需要通過冷凝器冷凝，氣體樣本在進(jìn)入色譜儀之前再次用變色硅膠進(jìn)行干燥，確保沒有水蒸氣進(jìn)入色譜儀。實(shí)驗(yàn)分別在不同的停留時(shí)間、水蒸氣/甲烷摩爾比和反應(yīng)溫度下進(jìn)行，各實(shí)驗(yàn)中壓力均為0.12 MPa。

圖1 介質(zhì)阻擋放電催化系統(tǒng)示意圖Fig.1 Diagram of dielectric barrier discharge-assisted system

利用可視化反應(yīng)器(見圖2(a))確定了等離子協(xié)同催化(discharge packed bed reactor，DPB)反應(yīng)器能夠正常工作的結(jié)構(gòu)形式為雙阻擋介質(zhì)結(jié)構(gòu)，即分別在內(nèi)、外電極的外側(cè)、內(nèi)側(cè)放置一環(huán)形阻擋介質(zhì)，見圖3。DPB反應(yīng)器中采用的催化劑為商用Ni質(zhì)量分?jǐn)?shù)為40%的 Ni/SiO2催化劑，顆粒大小為φ5 mm×2.5 mm圓柱的1/4(山東齊魯化學(xué)研究所供應(yīng))。

圖2 實(shí)驗(yàn)用反應(yīng)器Fig.2 Experimental reactor

圖3 定型DPB反應(yīng)器結(jié)構(gòu)示意圖Fig.3 Structural sketch of final DPB reactor

2 結(jié)果和討論

實(shí)驗(yàn)中采用甲烷轉(zhuǎn)化率χ及產(chǎn)物產(chǎn)量ζi作為反應(yīng)器性能評(píng)價(jià)標(biāo)準(zhǔn)，其定義如下:

2.1 動(dòng)力學(xué)模型建立

根據(jù)低溫等離子體化學(xué)理論，低溫等離子體化學(xué)催化過程主要包括電子碰撞反應(yīng)和自由基反應(yīng)。在折合電場(chǎng)作用下放電區(qū)域發(fā)生部分高能電子(1eV＜Te＜10eV)碰撞CH4和H2O分子發(fā)生激發(fā)、解離和電離等反應(yīng)，產(chǎn)生振動(dòng)活化組分、自由基、離子和電子等組分。其中自由基參與的反應(yīng)在甲烷轉(zhuǎn)化中占主導(dǎo)地位，離子在反應(yīng)中幾乎不起作用［11］。產(chǎn)生的自由基進(jìn)一步反應(yīng)生成合成氣(H2+CO)及低碳烴(C2/C3)。根據(jù)實(shí)驗(yàn)結(jié)果，在計(jì)算過程當(dāng)中包含了 C2、C3、C4、C5及其炭黑(C)。為了不使等離子體催化重整的化學(xué)反應(yīng)機(jī)理進(jìn)一步復(fù)雜化，機(jī)理的提出主要圍繞合成氣和低碳烴的產(chǎn)生過程進(jìn)行建立和優(yōu)化，而關(guān)于高碳烴的產(chǎn)生主要采用了2個(gè)反應(yīng)方程代替，并對(duì)其Arrhenius系數(shù)進(jìn)行部分修正。

等離子體催化過程中的電子碰撞反應(yīng)僅發(fā)生在高能電子群內(nèi)，即只有達(dá)到一定能量的電子才能有效碰撞反應(yīng)物分子，并以一定的概率發(fā)生化學(xué)反應(yīng)。采用全局平均電子能量替代局部高能電子能量，即建立電子碰撞反應(yīng)速率常數(shù)與平均電子能量的關(guān)系式。首先基于電子碰撞截面數(shù)據(jù)，利用BOLSIG+軟件［12］，求解Blotzmann方程，得到了反應(yīng)速率常數(shù)與平均電子能量之間的對(duì)應(yīng)關(guān)系;然后結(jié)合平均電子能量與電子溫度的關(guān)系式，通過ORIGIN進(jìn)行擬合，得到反應(yīng)速率常數(shù)關(guān)于電子溫度的Arrhenius形式關(guān)系式，而電子E也作為一種普通的氣相組分參與反應(yīng)［13］。自由基反應(yīng)的反應(yīng)速率常數(shù)通過 NIST Chemistry Webbook查詢得到。

低溫等離子體催化重整模型包括質(zhì)量守恒方程、物質(zhì)組分守恒方程及其氣體能量守恒方程和電子能量守恒方程。在計(jì)算過程中，假定自由基一旦產(chǎn)生便擴(kuò)散到整個(gè)區(qū)域，即為0維模型?；诘蜏氐入x子體催化理論知識(shí)及其實(shí)驗(yàn)過程中測(cè)定的溫度情況可知，等離子體催化過程中電子溫度較高，而中性組分和自由基的溫度基本保持不變，即與來流氣體溫度一致。因此在模型設(shè)置中假定氣體溫度保持不變。注意模型中沒有包括用于計(jì)算電勢(shì)的泊松方程，因此計(jì)算過程中沒有考慮電荷分布及其電極附近鞘層的形成等。選擇輸入功率P作為可以調(diào)節(jié)的參數(shù)，通過改變P保證甲烷出口濃度與實(shí)驗(yàn)數(shù)據(jù)接近，并根據(jù)相關(guān)文獻(xiàn)和實(shí)驗(yàn)數(shù)據(jù)，利用CHEMKINPRO軟件中的敏感性分析模塊，最終給出了包含175個(gè)反應(yīng)方程及39種組分的反應(yīng)方程組，使甲烷轉(zhuǎn)化及產(chǎn)物產(chǎn)量與實(shí)驗(yàn)結(jié)果基本吻合。反應(yīng)機(jī)理模型建立過程中采用壓力為1.2×10-6MPa，反應(yīng)溫度為400℃，停留時(shí)間為0.59 s，S/C 為2，輸入功率為70 W時(shí)的操作條件下的實(shí)驗(yàn)結(jié)果作為對(duì)照依據(jù)，該輸入功率一旦選定則保持不變，避免輸入功率變化對(duì)甲烷轉(zhuǎn)化率和產(chǎn)物分布的影響。

分子組分: ECH4H2O H2CO2COC2H6C2H4C2H2CH3OHCH3OCH2OC3H8C3H6C4H10C5H12，

自由基組分:CH3CH2CHCC2H5C2H3C2HHCOOH HO I*C3H7N*C3H7C4H9，

離子組分:CH4+CH3+H-CH2-OH-O-CH2+CH+C+。

2.2 模型的優(yōu)化和證實(shí)

圖4表明了實(shí)驗(yàn)研究中停留時(shí)間τ、水蒸氣/甲烷摩爾比(S/C)和反應(yīng)溫度T對(duì)甲烷轉(zhuǎn)化率和各產(chǎn)物產(chǎn)量的影響規(guī)律。為了證實(shí)得到的反應(yīng)機(jī)理方程，利用動(dòng)力學(xué)模型計(jì)算了不同操作條件下的甲烷轉(zhuǎn)化率和各產(chǎn)物組分的產(chǎn)量，如圖4所示。

圖4 不同操作條件對(duì)產(chǎn)物產(chǎn)量的影響規(guī)律Fig.4 The rules of different operating conditions on products yield

結(jié)果表明:該反應(yīng)機(jī)理能準(zhǔn)確預(yù)測(cè)不同操作條件下的甲烷轉(zhuǎn)化率的變化規(guī)律及其各產(chǎn)物產(chǎn)量的變化趨勢(shì)，其中C2H6的預(yù)測(cè)值與實(shí)驗(yàn)值偏差在4%以內(nèi)，但是該動(dòng)力學(xué)模型高估了H2的產(chǎn)量(偏差約為18% ～40%)、低估了CO(偏差約為60% ～80%)和C3H8的產(chǎn)量(偏差約為20% ～50%)。在實(shí)驗(yàn)過程中可以觀測(cè)到略顯黃色且?guī)в写碳ば詺馕兜挠偷?，這是在冷凝器中液化成油滴且隨著冷凝水流出的高碳烴。說明在實(shí)驗(yàn)過程中部分甲烷轉(zhuǎn)化為C4以上的高碳烴。由于在反應(yīng)機(jī)理中僅采用2個(gè)方程表明高碳烴的生成，因此對(duì)高碳烴產(chǎn)量的低估可能是高估H2產(chǎn)量的主要原因，從路徑圖5中也可以看出，大部分中間產(chǎn)物通過和H基的抽氫反應(yīng)生成了H2，這與文獻(xiàn)［14］分析結(jié)果一致。

圖5 反應(yīng)路徑圖Fig.5 Reaction path flux analysis

2.3 動(dòng)力學(xué)分析

圖4(a)表明隨著停留時(shí)間的增加，甲烷轉(zhuǎn)化率和產(chǎn)物產(chǎn)量均大幅度增加。這是因?yàn)楫?dāng)?shù)入x子體電源功率一定時(shí)，高能電子密度和能量基本保持不變，較長(zhǎng)的停留時(shí)間增加了電子碰撞反應(yīng)的概率，并且由于等離子體具有非平衡特性，高能電子碰撞CH4產(chǎn)生的 CH3之間發(fā)生偶聯(lián)反應(yīng)產(chǎn)生了 C2H6及其C3H8等，當(dāng)停留時(shí)間較長(zhǎng)時(shí)，C2、C3等組分又進(jìn)一步脫氫偶聯(lián)，生成了C4以上的高碳烴，使甲烷轉(zhuǎn)化率提高，具體見表1。

圖4(b)為S/C的影響，由于實(shí)驗(yàn)過程中為保證停留時(shí)間不變，S/C增加意味著甲烷的含量逐漸減小。從圖中可知隨著S/C增加，甲烷轉(zhuǎn)化率、H2和CO產(chǎn)量單調(diào)增加。從直方圖6(a)可知，影響甲烷轉(zhuǎn)化的主要反應(yīng)方程為R130、R129、R3及其R37，其中R130和R129分別為CH4與OH基和H基的反應(yīng)，而OH基和H基主要來自電子碰撞H2O解離反應(yīng)方程R12。圖7(a)表明了不同S/C對(duì)各主要反應(yīng)中的消耗速率的影響規(guī)律。從圖中可知，隨著S/C增加，R12中電子碰撞H2O的消耗速率大幅度增加，即電子碰撞H2O的概率增加，而R3中電子碰撞CH4的消耗速率略有下降;但由于OH基和H基的濃度增加，雖然甲烷含量降低，R129和R130的反應(yīng)速率基本保持不變。圖6(a)表明隨著S/C增加，CH4的總消耗速率是降低的，這與文獻(xiàn)［10］結(jié)果一致;但由于本實(shí)驗(yàn)中甲烷含量逐漸減小，因此相對(duì)消耗速率大幅度增加，CH4轉(zhuǎn)化率提高。

表1 不同停留時(shí)間對(duì)碳?xì)浠衔锂a(chǎn)量的影響Table 1 Effect of resident time on hydrocarbon yield %

圖4(c)表明了反應(yīng)溫度對(duì)甲烷轉(zhuǎn)化率和產(chǎn)物產(chǎn)量的影響規(guī)律。在數(shù)值模擬中，當(dāng)輸入功率保持不變即設(shè)定值為56 W時(shí)，甲烷轉(zhuǎn)化率呈線性增長(zhǎng)趨勢(shì)。根據(jù)圖7(b)可知，溫度對(duì)電子碰撞CH4和H2O(R3和R12)的反應(yīng)速率幾乎沒有影響，但隨著溫度增加，CH4和H基反應(yīng)速率增強(qiáng)(R129)，CH4的總消耗速率不斷增加。而實(shí)驗(yàn)研究表明隨著反應(yīng)溫度的增加，甲烷轉(zhuǎn)化率呈現(xiàn)先增加后減少的趨勢(shì)。這可能是由于溫度的增加導(dǎo)致石英玻璃管的電導(dǎo)率增加，歐姆損耗和介電損耗等增強(qiáng)，當(dāng)輸入功率不變時(shí)，作用于反應(yīng)氣體的功率減少［15］。因此在數(shù)值研究中，隨著反應(yīng)溫度的增加，相應(yīng)地降低輸入功率，對(duì)應(yīng)于反應(yīng)溫度為450℃和500℃時(shí)的功率分別調(diào)整為54 W和52 W，此時(shí)CH4轉(zhuǎn)化率結(jié)果與實(shí)驗(yàn)值吻合較好，這也證實(shí)了溫度對(duì)阻擋介質(zhì)材料屬性的影響。

圖6 CH4和H2O的主要的產(chǎn)生和消耗反應(yīng)的反應(yīng)速率Fig.6 Time-averaged reaction rates of the dominant reaction pathways for the production and consumption of CH4and H2O

圖7 S/C和反應(yīng)溫度對(duì)反應(yīng)物消耗速率的影響Fig.7 Consumption rate of reactants with increasing S/C，wall temperature

為了明確等離子體催化甲烷蒸汽重整過程中影響產(chǎn)物選擇性的主要反應(yīng)，結(jié)合優(yōu)化的反應(yīng)機(jī)理，利用CHEMKIN-PRO軟件中的路徑流分析模塊，進(jìn)行了反應(yīng)路徑流分析。操作條件是壓力為1.2×10-6MPa，反應(yīng)溫度為400 ℃，停留時(shí)間為0.59 s，S/C為2，輸入功率為70 W。不同操作條件下的反應(yīng)路徑基本一致，只是百分?jǐn)?shù)略有不同。甲烷主要通過電子碰撞反應(yīng)及其與OH基和H基結(jié)合產(chǎn)生CH3，而CH2、CH及C的產(chǎn)量較少。從實(shí)驗(yàn)結(jié)果可知，CO產(chǎn)量較低，僅為1%，低碳烴(C2、C3)產(chǎn)量在4%左右，另一部分甲烷轉(zhuǎn)化為高碳烴化合物。從路徑圖7分析可知，60%的CH3基結(jié)合生成了C2H6，并進(jìn)一步通過脫氫偶聯(lián)反應(yīng)轉(zhuǎn)化為高碳烴;僅有1.1%的CH3基和O基結(jié)合通過路徑CH3+O→CH2O→HCO→CO最終生成了CO，而CO幾乎沒有轉(zhuǎn)化為CO2，這與實(shí)驗(yàn)結(jié)果相吻合。

根據(jù)方程(1)、(2)可知，甲烷蒸汽產(chǎn)生合成氣的反應(yīng)為強(qiáng)吸熱反應(yīng)，而甲烷偶聯(lián)生成碳烴化合物的反應(yīng)則為弱吸熱反應(yīng)，不利于余熱回收。因此增加合成氣的產(chǎn)量和避免高碳烴的生成是提高余熱回收利用的主要手段，結(jié)合路徑分析圖，主要應(yīng)采取2個(gè)措施:1)抑制CH3基再結(jié)合生成C2H6;2)提高O基濃度，即提高CH3基與O基結(jié)合的反應(yīng)速率。通過H2O生成和消耗速率圖可知，通過電子碰撞反應(yīng)83%的反應(yīng)的H2O的生成了OH基和H基，僅有3%的反應(yīng)的H2O生成了O基，而且生成的OH基的大部分又和CH4再結(jié)合生成了H2O。因此抑制OH基與CH4的結(jié)合，同時(shí)促進(jìn)OH→O是提高O基濃度的關(guān)鍵。

由催化劑表面的甲烷蒸汽重整反應(yīng)機(jī)理可知［16］，通過選擇合適的催化劑，如 Ni催化劑，可以保證通過高能電子碰撞反應(yīng)產(chǎn)生的自由基及時(shí)吸附到催化劑表面，利用催化劑的選擇性特點(diǎn)，降低烷基與烷基的再結(jié)合概率，提高OH基向O基的轉(zhuǎn)化，從而提高合成氣的產(chǎn)量，降低高碳烴的產(chǎn)量，見式(3)～(10):

2.4 等離子協(xié)同催化實(shí)驗(yàn)驗(yàn)證

根據(jù)實(shí)驗(yàn)及數(shù)值研究結(jié)果可知，甲烷蒸汽重整過程中會(huì)產(chǎn)生高碳烴，該反應(yīng)過程近似絕熱甚至是放熱，不利于余熱的回收利用;且在化學(xué)回?zé)嵫h(huán)燃?xì)廨啓C(jī)中，化學(xué)回?zé)崞鳟a(chǎn)生的裂解氣將直接注入燃燒室燃燒，較高分子量的燃料不利于提高燃燒效率和拓寬熄火極限。因此，作者首次提出有效碳回收率的概念，用以表征裂解氣中碳的存在形式，作為衡量回?zé)崽匦院腿紵匦缘暮暧^評(píng)價(jià)指標(biāo)。

利用DPB反應(yīng)器，研究了壓力為1.2，反應(yīng)溫度為400℃，停留時(shí)間為0.39 s，S/C為2，輸入功率為70 W時(shí)的操作條件下的甲烷蒸汽重整情況，各重整特性參數(shù)見表2。從表中可知，在放電區(qū)域加入催化劑后，合成氣的產(chǎn)量有了顯著提高，C2及C3產(chǎn)量為0%，有效碳回收率達(dá)到了100%，即產(chǎn)物中沒有生成高碳烴等。因此該實(shí)驗(yàn)結(jié)果與動(dòng)力學(xué)預(yù)測(cè)結(jié)果相吻合，在接下來的工作中，作者將結(jié)合實(shí)驗(yàn)數(shù)據(jù)，在等離子體催化甲烷蒸汽重整反應(yīng)機(jī)理模型的基礎(chǔ)上進(jìn)一步研究等離子體與催化劑協(xié)同催化時(shí)的反應(yīng)機(jī)理及其重整特性。

表2 不同催化反應(yīng)器下重整結(jié)果對(duì)照表Table 2 Comparison of reforming results under different reactors

3 結(jié)論

1)實(shí)驗(yàn)研究發(fā)現(xiàn)DBD輔助甲烷蒸汽重整在停留時(shí)間為0.59 s時(shí)有較高的轉(zhuǎn)化率，但有效碳回收率較低，不同操作條件下CO產(chǎn)量為1%，低碳烴(C2、C3)產(chǎn)量為4%左右，剩余CH4轉(zhuǎn)化為高碳烴(C4以上);

2)發(fā)展并完善的動(dòng)力學(xué)機(jī)理可以準(zhǔn)確預(yù)測(cè)不同操作條件下DBD催化重整過程中的甲烷轉(zhuǎn)化率及其各產(chǎn)物產(chǎn)量的變化趨勢(shì)。其中CH3+CH3→C2H6和 OH+CH4→H2O+CH3是影響CO產(chǎn)量的主要限制反應(yīng);

3)等離子與催化劑協(xié)同催化實(shí)驗(yàn)中有效碳回收率達(dá)到100%，通過在放電區(qū)域添加Ni催化劑，可以提高合成氣的產(chǎn)量，避免高碳烴生成。在接下來的工作中將進(jìn)一步研究其動(dòng)力學(xué)機(jī)理。

［1］CARCSACI C，F(xiàn)ACCHINI B，HARVEY S.Modular approach to analysis of chemically recuperated gas turbine cycles［J］.Energy Conversion Management，1998;39(16/18):1693-1703.

［2］PAN Fumin，ZHENG Hongtao，LIU Qingzhen，et al.Design and performance calculations of chemically recuperated gas turbine on ship［J］.Proc IMechE Part A:J Power and Energy，2013，227(8):908-918.

［3］HAMMER T，KAPPES T，BALDAUF M.Plasma catalytic hybrid processes:gas discharge initiation and plasma activation of catalytic processes［J］.Catalysis Today，2004，89:5-14.

［4］顏士鑫，李曉東，鐘犁，等.滑動(dòng)弧等離子體協(xié)助甲烷蒸汽重整制氫［J］.太陽能學(xué)報(bào)，2011，32(5):766-770.YAN Shixin，LI Xiaodong，ZHONG Li，et al.Gliding arc discharge plasma assisted reforming of methane into hydrogen with oxygen and water vapor［J］.Acta Energiae Solaris Sinica，2011，32(5):766-770.

［5］RUSU I，CORMIER J M.On a possible mechanism of the methane steam reforming in a gliding arc reactor［J］.Chemical Engineering Journal，2003，91:23-31.

［6］PETROVI D，MARTENS T，Van DIJK J，et al.Modeling of a dielectric barrier discharge used as a flowing chemical reactor［J］.Journal of Physics:Conference Series，2008，133:1-8.

［7］NAIR S A，NOZAKI T，OKAZAKI K.Methane oxidative conversion pathways in a dielectric barrier discharge reactor-Investigation of gas phase mechanism［J］.Chemical Engineering Journal，2007，132:85-95.

［8］李潔，王彧婕，龍華麗，等.強(qiáng)電場(chǎng)中甲烷活化的湯生模型及反應(yīng)動(dòng)力學(xué)分析［J］.化學(xué)工程，2007，35(8):25-28.LI Jie，WANG Yujie，LONG Hihua，et al.Townsend ionization model and kinetic analysis of methane activated in strong electric field［J］.Chemical Engineering，2007，35(8):25-28.

［9］SUGASAWA M，TERASAWA T，F(xiàn)UTAMURA S.Effects of initial water content on steam reforming of aliphatic hydrocarbons with nonthermal plasma［J］.Journal of Electrostatics，2010，68:212-217.

［10］ZHANG X，WANG B W，LIU Y W，et al.Conversion of methane by steam reforming using dielectric-barrier discharge［J］.Chinese Journal of Chemical Engineering，2009，17(4):625-629.

［11］HIRAOKA K，AOYAMA K，MORISE K.A study of reaction mechanisms of methane in a radio-frequency glow discharge plasma using radical and ion scavengers［J］.Canadina Journal of Chemistry，1985，63:2899-2905.

［12］HAGELAAR G J M，PITCHFORD L C.Solving the Boltzmann equation to obtain electron transport coefficients and rate coefficients for fluid models［J］.Plasma Sci Sources and Tech，2005，14:722.

［13］KEE R J，RUPLEY F M.，MEEKS E，et al.CHEMKINIII:a Fortran chemical kinetics package for the analysis of gas-phase chemical and plasma kinetics.SAND96-8216，1996 CHEMKIN-PRO［R］.San Diego:Reaction Design，2008.

［14］SUN W T，UDDI M，WON S H，et al.Kinetic effects of non-equilibrium plasma-assisted methane oxidation on diffusion flame extinction limits［J］.Combustion and Flame，2012，159(1):221-229.

［15］HAMMER T H，KAPPES T H，BALDAUF M.Plasma catalytic hybrid processes:gas discharge initiation and plasma activation of catalytic processes［J］.Catalysis Today，2004，89:5-14.

［16］JANARDHANAN V M，DEUTSCHMANN O.CFD analysis of a solid oxide fuel cell with internal reforming:coupled interactions of transport，heterogeneous catalysis and electrochemical processes［J］.J Power Sources，2006，162(2):1192-1202.