王堂博,王廣建,孫萬(wàn)堂,牛鑫善,王 芳
?
綜述與展望
異丁烷催化脫氫反應(yīng)機(jī)理與失活動(dòng)力學(xué)研究進(jìn)展
王堂博,王廣建*,孫萬(wàn)堂,牛鑫善,王芳
(青島科技大學(xué)化工學(xué)院,山東 青島 266042)
摘要:異丁烯是化工行業(yè)重要的基礎(chǔ)原料,國(guó)內(nèi)外對(duì)異丁烯的需求量逐年遞增。僅靠石油催化裂解已無(wú)法滿足對(duì)異丁烯的需求,開(kāi)展異丁烷脫氫制異丁烯工藝的研究備受關(guān)注。綜述Cr系異丁烷脫氫催化劑的研究進(jìn)展,探討Cr系催化劑的活性中心以及發(fā)生在活性中心上的多種反應(yīng)機(jī)理和相應(yīng)的動(dòng)力學(xué)模型,詳述催化劑的失活機(jī)理,總結(jié)積炭的形成過(guò)程,指出Cr系異丁烷催化脫氫反應(yīng)和失活機(jī)理以及相關(guān)動(dòng)力學(xué)方面研究的不足,并對(duì)未來(lái)研究前景進(jìn)行展望。
關(guān)鍵詞:催化反應(yīng)工程;異丁烷;異丁烯;脫氫;動(dòng)力學(xué);反應(yīng)機(jī)理
CLC number:O643.12;TQ221.21+3Document code: AArticle ID: 1008-1143(2016)04-0001-06
異丁烯是重要的化工原料,主要用于制備甲基叔丁基醚、丁基橡膠、聚異丁烯和各種精細(xì)化學(xué)品[1-2]。異丁烯下游產(chǎn)品的需求量激增促進(jìn)了異丁烯產(chǎn)業(yè)的飛速發(fā)展,預(yù)計(jì)2016年異丁烯產(chǎn)量將達(dá)到29.78 Mt[3-4]。我國(guó)的異丁烯工業(yè)始于20世紀(jì)70年代,雖然近年來(lái)我國(guó)異丁烯產(chǎn)量有所增加,但遠(yuǎn)遠(yuǎn)不能滿足工業(yè)發(fā)展的需求。異丁烷脫氫轉(zhuǎn)化為異丁烯可緩解異丁烯供不應(yīng)求的矛盾,副產(chǎn)物氫氣也可用作加氫裝置的原料,是實(shí)現(xiàn)天然氣、煉廠氣資源優(yōu)化利用的重要途徑,未來(lái)應(yīng)用市場(chǎng)前景廣闊。
異丁烷脫氫技術(shù)關(guān)鍵在于催化劑的研制與選擇,催化劑決定著異丁烷脫氫工藝的流程和經(jīng)濟(jì)效益,Cr系催化劑具有成本低、催化活性好等優(yōu)點(diǎn),但穩(wěn)定性需要提高[5-7]。本文對(duì)異丁烷脫氫技術(shù)中Cr系催化劑的活性中心、反應(yīng)與失活機(jī)理和動(dòng)力學(xué)等方面的研究進(jìn)行綜述。
1Cr系催化劑的活性中心與反應(yīng)機(jī)理
1.1催化劑活性中心
異丁烷吸附主要發(fā)生在Cr—O位上,充分研究在異丁烷脫氫環(huán)境下Cr—O位的性質(zhì)對(duì)確定催化劑活性中心非常重要。研究表明,新制備的Cr2O3/Al2O3催化劑中存在大量的表面物種,包括Cr6+、Cr5+、Cr3+和Cr2+。各個(gè)物種的濃度與Cr負(fù)載量、催化劑制備過(guò)程中的焙燒溫度和選擇的載體有關(guān)。低Cr負(fù)載量時(shí),催化劑含有Cr6+和少量的Cr5+,此時(shí)Cr6+分為可溶性Cr6+和非可溶性Cr6+,Cr6+在載體上的密度隨著Cr負(fù)載量的增加而增加[8],最終可溶性的Cr6+穩(wěn)定在(5~10)個(gè)Cr·nm-2,非可溶性的Cr6+穩(wěn)定在2個(gè)Cr·nm-2。高Cr負(fù)載量時(shí),除了Cr6+、Cr5+、Cr3+和Cr2+外,還可以觀察到高度對(duì)稱的Cr2O3和嵌進(jìn)微晶相中的Cr2O3[9-10]。因?yàn)镃r2O3和SBA-15之間的相互作用較弱,相比于Al2O3載體,SBA-15上更容易形成高度對(duì)稱的Cr2O3和嵌進(jìn)微晶相中的Cr2O3。
部分研究者認(rèn)為,分散在催化劑表面的Cr6+是異丁烷脫氫反應(yīng)的主要活性中心,Kilicarslan S等[11]通過(guò)XPS分析觀察到在脫氫反應(yīng)完成后,Cr6+(2p3/2)在578 eV的峰消失,Cr3+(2p3/2)峰出現(xiàn)。Elzinga E J等[12]認(rèn)為脫氫反應(yīng)后,四面體結(jié)構(gòu)的Cr(Ⅵ)O4轉(zhuǎn)化為八面體結(jié)構(gòu)的Cr(Ⅲ)O6,八面體的Cr(Ⅲ)O6高度對(duì)稱,Cr在此結(jié)構(gòu)中沒(méi)有轉(zhuǎn)移電子的能力,隨著Cr(Ⅲ)O6的積累形成了Cr2O3結(jié)晶相,催化劑失去催化活性。也有學(xué)者持不同意見(jiàn),Weckhuysen B M等[13]通過(guò)XANES和紫外可見(jiàn)光譜觀察到在反應(yīng)初期Cr3+由Cr6+形成,并認(rèn)為Cr系催化劑脫氫活性主要?dú)w咎于Cr3+物種。催化劑表面存在3種類型的Cr3+物種:(1) 由Cr6+反應(yīng)后得到的氧化還原態(tài)的Cr3+;(2) 分散在催化劑表面的配位不飽和Cr3+,這兩類Cr3+物種都具有催化活性;(3) 高度對(duì)稱的Cr3+和微晶相中的Cr3+,微晶相中的Cr3+嵌進(jìn)載體Al2O3中被覆蓋,所以催化活性較低。具有催化活性的Cr3+物種通過(guò)吸附氣相中的異丁烷分子形成穩(wěn)定的對(duì)稱Cr3+物種。
1.2異丁烷脫氫反應(yīng)機(jī)理和動(dòng)力學(xué)
異丁烷脫氫是氣-固多相反應(yīng),動(dòng)力學(xué)模型復(fù)雜,總過(guò)程速率由反應(yīng)最慢的基元反應(yīng)決定,即速控步驟。多年來(lái),由于制備的催化劑物理化學(xué)性質(zhì)不同、反應(yīng)條件以及測(cè)試手段存在差異,有關(guān)異丁烷脫氫反應(yīng)的機(jī)理和速控步驟一直存在爭(zhēng)議。在上述活性位研究的基礎(chǔ)上,按照催化劑表面的活性中心數(shù)量以及氧離子是否參加反應(yīng),將反應(yīng)分為單一Cr—O活性位機(jī)理、協(xié)同Cr—O活性位機(jī)理和無(wú)氧協(xié)同活性位機(jī)理。
1.2.1單一Cr—O活性位機(jī)理
Burwell Jr R L等[14]提出異丁烷脫氫反應(yīng)是在一個(gè)單一Cr—O活性位上發(fā)生的,催化劑上異丁烷脫氫生成異丁烯分為以下3個(gè)步驟,如圖1所示。i-C4H9…Cr—O…H表示烷基和氫原子吸附在Cr—O活性位上,a為反應(yīng)中異丁烷的吸附過(guò)程;b為表面反應(yīng)過(guò)程。
圖 1 單一Cr—O活性位機(jī)理Figure 1 Reaction mechanism of singleCr—O active site
異丁烷在不飽和的Crn+中心吸附,然后異丁烷中的C—H鍵斷裂,O—H鍵和Cr—C鍵生成,最后異丁烯脫附,在催化劑表面生成H2,Crn+活性中心復(fù)原[15],認(rèn)為反應(yīng)中以烷烴的吸附過(guò)程(步驟1a)作為異丁烷脫氫過(guò)程中的速控步驟。
Weckhuysen B M等[16]提出了一個(gè)相似的異丁烷脫氫機(jī)理,只是將烷烴吸附與脫附后形成的烷基和氫原子分為兩個(gè)式子進(jìn)行表達(dá),認(rèn)為該反應(yīng)機(jī)理與Burwell Jr R L提出的機(jī)理相似。根據(jù)異丁烷脫氫反應(yīng)機(jī)理的結(jié)果建立了Cr系脫氫催化劑的動(dòng)力學(xué)模型。
Airaksinen S M K等[17]利用流化床反應(yīng)器,在常壓和(520~580) ℃條件下,對(duì)Cr系異丁烷催化劑的脫氫反應(yīng)動(dòng)力學(xué)進(jìn)行研究,分別以反應(yīng)過(guò)程中異丁烷的吸附和表面反應(yīng)作為異丁烷脫氫中的速控步驟,利用Langmuir-Hinshelwood原理建立動(dòng)力學(xué)方程,采用數(shù)學(xué)方法優(yōu)化各參數(shù),得到動(dòng)力學(xué)方程,如式(1)所示:
(1)
1.2.2協(xié)同Cr—O活性位機(jī)理
Carra S等[18]提出了協(xié)同Cr—O活性位機(jī)理,并假設(shè)活性位需要Cr和O離子同時(shí)存在,一些步驟發(fā)生在同一個(gè)Cr—O活性位上,同時(shí)用到Cr和O離子。烷烴在Cr—O鍵上發(fā)生吸附作用,吸附的烷基與相鄰的Cr—O鍵發(fā)生作用,最后烯烴和氫氣從表面脫附,如圖2所示。
圖 2 協(xié)同Cr—O活性位機(jī)理Figure 2 Reaction mechanism ofcollaborative Cr—O active site
Carra S以機(jī)理中的表面反應(yīng)(步驟2b)作為速控步驟,采用反應(yīng)動(dòng)力學(xué)方程(2)表示:
(2)
1.2.3無(wú)氧協(xié)同活性位機(jī)理
Kao J Y等[19]考察異丁烷在無(wú)載體Cr催化劑上的脫氫反應(yīng),提出無(wú)氧協(xié)同活性位機(jī)理,見(jiàn)圖3。
圖 3 無(wú)氧協(xié)同活性位機(jī)理Figure 3 Reaction mechanism of Cr active site
Kao J Y以無(wú)氧協(xié)同活性位機(jī)理中的表面反應(yīng)(步驟3b)為速控步驟,反應(yīng)動(dòng)力學(xué)方程見(jiàn)式(3):
(3)
Airaksinen S M K等[17]對(duì)上述三種動(dòng)力學(xué)模型分別進(jìn)行了殘差平方和SSR分析,結(jié)果顯示,模型(1)的擬合程度最好。單一Cr—O活性位機(jī)理中烷烴分子的吸附發(fā)生在獨(dú)立的Cr—O活性位上,協(xié)同Cr—O活性位機(jī)理和無(wú)氧協(xié)同活性位機(jī)理中假設(shè)的表面反應(yīng)(速控步驟)則需要多個(gè)相同的活性位同時(shí)存在。
Airaksinen S M K認(rèn)為,單一Cr—O活性位機(jī)理更加接近實(shí)際烷烴脫氫過(guò)程,所以上述模型的分析結(jié)果偏向于單一Cr—O活性位機(jī)理。
作者認(rèn)為異丁烷脫氫反應(yīng)機(jī)理可簡(jiǎn)化由圖4表述,烷烴在催化劑活性位上的吸附過(guò)程為反應(yīng)的速控步驟,今后的研究可以結(jié)合同位素標(biāo)記法,對(duì)低碳烷烴的脫氫機(jī)理進(jìn)行進(jìn)一步驗(yàn)證。
圖 4 異丁烷在CrOx催化劑上的脫氫反應(yīng)機(jī)理Figure 4 Reaction mechanisms for dehydrogenation of iso-butane on CrOxcatalyst
異丁烷脫氫生成異丁烯的反應(yīng)機(jī)理和反應(yīng)動(dòng)力學(xué)是異丁烷脫氫反應(yīng)理論研究的基礎(chǔ),對(duì)異丁烷轉(zhuǎn)化為異丁烯具有重要的指導(dǎo)意義。上述反應(yīng)動(dòng)力學(xué)方程是利用Langmuir-Hinshelwood原理分析反應(yīng)機(jī)理后得到,由于未考慮反應(yīng)中傳熱、傳質(zhì)等問(wèn)題,可能會(huì)造成偏差,導(dǎo)致動(dòng)力學(xué)模型不能客觀反映異丁烷脫氫的實(shí)際反應(yīng)歷程,有待于在今后的實(shí)驗(yàn)中進(jìn)一步修正和完善。
2催化劑的積炭失活機(jī)理
異丁烷制備異丁烯的失活過(guò)程研究目前主要集中在積炭反應(yīng)上,積炭的存在降低了催化劑活性,并且破壞載體的孔道結(jié)構(gòu),使得催化劑進(jìn)入非永久性失活狀態(tài)。
文獻(xiàn)[20]通過(guò)使用微量振蕩天平進(jìn)行實(shí)驗(yàn),結(jié)果表明,在反應(yīng)初期(10~40) min,催化劑的質(zhì)量隨著時(shí)間的延長(zhǎng)而增加,隨后保持穩(wěn)定。但是異丁烷轉(zhuǎn)化率隨著反應(yīng)時(shí)間的延長(zhǎng)依舊下降,反應(yīng)形成的積炭使催化劑中毒,逐漸失去活性。Cr系催化劑表面積炭機(jī)理比較復(fù)雜,如圖5所示,積炭通過(guò)活性中心深度脫氫、裂化、載體酸性位表面的聚合環(huán)化等過(guò)程形成。
圖 5 異丁烷脫氫過(guò)程中發(fā)生的裂化、深度脫氫和積炭反應(yīng)Figure 5 Proposed process of cracking,deep dehydrogenation and carbon depositionduring isobutane dehydrogenation reaction
Airaksinen S等[21]通過(guò)光譜、質(zhì)譜、色譜的聯(lián)用對(duì)積炭形成過(guò)程進(jìn)行研究,發(fā)現(xiàn)乙酸鹽的羧化物和脂肪族烴類在反應(yīng)初始階段形成,然后芳香烴開(kāi)始在催化劑上積累,最后催化劑失活。氫氣預(yù)處理、原料氣中加入氧化劑、改善載體酸性位、提高活性組分在催化劑表面的分散度可以有效抑制積炭,防止催化劑失活。通過(guò)分階段逐步升溫的方法進(jìn)行燒炭可以部分恢復(fù)催化劑的活性。但再生過(guò)程中燒炭所釋放的熱量導(dǎo)致活性Cr物種與Al2O3結(jié)合,形成一種穩(wěn)定的Cr2O3/Al2O3尖晶石結(jié)構(gòu),造成催化劑永久性失活。目前異丁烷脫氫制備異丁烯反應(yīng)失活過(guò)程的研究重點(diǎn)主要是積炭造成的非永久性失活。同時(shí),上述討論并未涉及碳母體問(wèn)題,無(wú)法對(duì)一些積炭現(xiàn)象做出合理解釋,因此,有必要研究碳母體的來(lái)源問(wèn)題。
目前,許多學(xué)者提出異丁烷脫氫反應(yīng)產(chǎn)生的積炭與反應(yīng)物系的組成和濃度有關(guān)。認(rèn)為催化劑顆粒上的炭分布是由表面向中心減少,但究竟與何因素相關(guān),至今仍沒(méi)有一致結(jié)論。Gascón J等[22]提出傳統(tǒng)經(jīng)典的積炭動(dòng)力學(xué)并沿用至今,單-多層積炭動(dòng)力學(xué)模型(MMCGM)認(rèn)為積炭的形成與濃度和溫度有關(guān)。催化劑表面的積炭濃度與時(shí)間關(guān)系有兩種,第一種是脫氫瞬間形成的瞬時(shí)表面積炭;第二種是多層積炭。積炭速率由單層瞬時(shí)積炭速率和多層積炭速率共同得到,動(dòng)力學(xué)方程式如下:
式中,t為反應(yīng)時(shí)間;cc代表積炭總濃度;cm和cM分別代表單層和多層積炭濃度;k1c和k2c分別代表單層和多層動(dòng)力學(xué)常數(shù)。
單層積炭量與分散在催化劑表面的活性位數(shù)量成正比,多層積炭伴隨著單層積炭的出現(xiàn)而出現(xiàn),多層積炭量與被單層積炭覆蓋的活性位的量成正比。
3結(jié)語(yǔ)
異丁烷制備異丁烯過(guò)程中的反應(yīng)與失活機(jī)理和反應(yīng)動(dòng)力學(xué)一直是該領(lǐng)域的焦點(diǎn)與熱點(diǎn),由于反應(yīng)過(guò)程的復(fù)雜性和多樣性,對(duì)其反應(yīng)機(jī)理的研究帶來(lái)巨大的挑戰(zhàn)。過(guò)去的研究大都集中在本征動(dòng)力學(xué)方面。因此,需要對(duì)Cr系催化劑脫氫反應(yīng)機(jī)理和反應(yīng)動(dòng)力學(xué)進(jìn)行更深一步的研究,尤其加強(qiáng)對(duì)宏觀動(dòng)力學(xué)方面的研究,為高性能催化劑的設(shè)計(jì)制備以及新工藝技術(shù)的開(kāi)發(fā)提供理論依據(jù)。同時(shí),研究新型催化劑與新反應(yīng)體系的動(dòng)力學(xué)機(jī)理,較全面進(jìn)行相應(yīng)宏觀動(dòng)力學(xué)模型研究是未來(lái)研究的方向。
參考文獻(xiàn):
[1]程正載,王洋,陳攀,等.混合碳四資源深加工綜合利用[J].精細(xì)石油化工進(jìn)展,2012,38(11):1245-1252.Cheng Zhengzai,Wang Yang,Chen Pan,et al.Deep processing and comprehensive utilization of mixed four carbon fraction[J].Advances in Fine Petrochemicals,2012,38(11):1245-1252.
[2]吳麗瓊,董炳利.丁烷催化脫氫催化劑及技術(shù)研究進(jìn)展[J].石油技術(shù)與應(yīng)用,2014,32(2):177-181.
Wu Liqiong,Dong Bingli.Research progress in catalysts and production technologies of butane catalytic dehydrogenation[J].Petrochemical Technology & Application,2014,32(2):177-181.
[3]張凌峰,劉亞錄,胡忠攀,等.丙烷脫氫制丙烯催化劑研究的進(jìn)展[J].石油學(xué)報(bào)(石油加工),2015,31(2):400-412.
Zhang Linfeng,Liu Yalu,Hu Zhongpan,et al.Advance in catalysts for propane dehydrogenation to propylene[J].Petroleum Processing and Petrochemicals,2015,31(2):400-412.
[4]崔丹丹,周紅軍.正丁烷脫氫異構(gòu)制異丁烯的研究進(jìn)展[J].工業(yè)催化,2014,22(12):901-904.
Cui Dandan,Zhou Hongjun. Research progress inn-butane dehydroisomerization to isobutene[J].Industrial Catalysis,2014,22(12):901-904.
[5]Ajayi B P,Jermy B R,Ogunronbi K E,et al.n-Butane dehydrogenation over mono and bimetallic MCM-41 catalysts under oxygen free atmosphere[J].Catalysis Today,2013,204:189-196.
[6]Inaba M,Murata K,Saito M,et al.Dehydrogenation and isomerization ofn-butane or isobutane over Cr catalysts supported on zeolites[J].Catalysis Letters,2002,84(3/4):273-279.
[7]Zhao H,Song H,Chou L.Synthesis and catalytic application in isobutane dehydrogenation of the mesoporous chromia/alumina catalysts based on a metal-organic framework[J].Microporous and Mesoporous Materials,2013,181:182-191.
[8]Shee D,Sayari A.Light alkane dehydrogenation over mesoporous Cr2O3/Al2O3catalysts[J].Applied Catalysis A:General,2010,389(1):155-164.
[9]Conley M P,Delley M F,Siddiqi G,et al.Polymerization of ethylene by silica-supported dinuclear CrⅢ sites through an initiation step involving C—H bond activation[J].Angewandte Chemie International Edition,2014,53(7): 1872-1876.
[10]Chakrabarti A,Wachs I E.The nature of surface CrOxsites on SiO2in different environments[J].Catalysis Letters,2015,145(4):985-994.
[11]Kilicarslan S,Dogan M,Dogu T.Cr incorporated MCM-41 type catalysts for isobutane dehydrogenation and deactivation mechanism[J].Industrial & Engineering Chemistry Research,2013,52(10):3674-3682.
[12]Elzinga E J,Cirmo A.Application of sequential extractions and X-ray absorption spectroscopy to determine the speciation of chromium in Northern New Jersey marsh soils developed in chromite ore processing residue(COPR)[J].Journal of Hazardous Materials,2010,183(1):145-154.
[13]Weckhuysen B M,Verberckmoes A A,Debaere J,et al.In situ UV-Vis diffuse reflectance spectros copy on line activity measurements of supported chromium oxide catalysts:relating isobutane dehydrogenation activity with Cr-speciation via experimental design[J].Journal of Molecular Catalysis A:Chemical,2000,151(1):115-131.
[14]Burwell Jr R L,Littlewood A B,Cardew M,et al.Reactions between hydrocarbons and deuterium on chromium oxide gel Ⅰ.General[J].Journal of the American Chemical Society,1960,82(24):6272-6280.
[15]Pasha F A,Cavallo L,Basset J M.Mechanism ofn-butane hydrogenolysis promoted by Ta-hydrides supported on silica[J].ACS Catalysis,2014,4(6):1868-1874.
[16]Weckhuysen B M,Schoonheydt R A.Alkane dehydrogenation over supported chromium oxide catalysts[J].Catalysis Today,1999,51(2): 223-232.
[17]Airaksinen S M K,Harlin M E,Krause A O I.Kinetic modeling of dehydrogenation of isobutane on chromia/alumina catalyst[J].Industrial & Engineering Chemistry Research,2002,41(23):5619-5626.
[18]Carra S,Forni L,Vintani C.Kinetics and mechanism in catalytic dehydrogenation ofn-butane over chromia-alumina[J].Journal of Catalysis,1967,9(2):154-165.
[19]Kao J Y,Piet-Lahanier H,Walter E,et al.Reversible reactions over non-metallic catalysts:a transient isotopic tracing of the isobutane-isobutene-hydrogen system over chromia[J].Journal of Catalysis,1992,133(2):383-396.
[20]Park H,Yun Y S,Kim T Y,et al.Kinetics of the dehydration of glycerol over acid catalysts with an investigation of deactivation mechanism by coke[J].Applied Catalysis B:Environmental,2015,176:1-10.
[21]Airaksinen S,Baares M A,Krause A I O.In situ characterisation of carbon-containing species formed on chromia/alumina during propane dehydrogenation[J].Journal of Catalysis,2005,230:507-513.
[22]Gascón J,Téllez C,Herguido J,et al.Propane dehydrogenation over a Cr2O3/Al2O3catalyst: transient kinetic modeling of propene and coke formation[J].Applied Catalysis A:General,2003,248(1):105-116.
Advance in catalytic reaction mechanisms and deactivation kinetics of isobutane dehydrogenation
WangTangbo,WangGuangjian*,SunWantang,NiuXinshan,WangFang
(Department of Chemical Engineering,Qingdao University of Science and Technology,Qingdao 266042,Shandong,China)
Abstract:Isobutene is an important raw material in petrolchemical industry,and the demand for isobutene is continuously increasing.The technologies for dehydrogenation of isobutane to isobutene have received much attention because the conventional petroleum cracking cannot satisfy the growing demand for isobutene.In this paper,the recent development in chromium-based catalysts for isobutane dehydrogenation was reviewed.The active center of chromium-based catalyst,a variety of reaction mechanisms and their corresponding dynamic model based on the active center were discussed.The catalyst deactivation mechanism and the formation process of carbon deposition were also summaried.The shortage of researches on the reaction and deactivation mechanisms of chromium-based catalysts for isobutane catalytic dehydrogenation and the related reaction kinetic study were pointed out.The future research prospects of the chromium-based catalysts for isobutane dehydrogenation are also outlined.
Key words:catalytic reaction engineering; isobutane; isobutene; dehydrogenation; kinetics; reaction mechanism
收稿日期:2015-09-24
基金項(xiàng)目:國(guó)家自然科學(xué)基金(21276130)資助項(xiàng)目;山東省自然科學(xué)基金(ZR2014BP009)資助項(xiàng)目
作者簡(jiǎn)介:王堂博,男,1990年生,在讀碩士研究生。
doi:10.3969/j.issn.1008-1143.2016.04.001 10.3969/j.issn.1008-1143.2016.04.001
中圖分類號(hào):O643.12;TQ221.21+3
文獻(xiàn)標(biāo)識(shí)碼:A
文章編號(hào):1008-1143(2016)04-0001-06
通訊聯(lián)系人:王廣建,男,博士,教授,主要從事環(huán)境凈化催化工程和催化新材料等研究。