莊 巖,蔣斌波,王靖岱,陽永榮
(浙江大學(xué) 化學(xué)工程聯(lián)合國家重點實驗室,浙江 杭州310027)
ZnO催化對苯二甲酸脫羧制苯的反應(yīng)歷程
莊 巖,蔣斌波,王靖岱,陽永榮
(浙江大學(xué) 化學(xué)工程聯(lián)合國家重點實驗室,浙江 杭州310027)
通過比較不同金屬氧化物催化對苯二甲酸脫羧制苯的活性及反應(yīng)產(chǎn)物分布情況,優(yōu)選ZnO作為脫羧催化劑。采用FT-IR、XRD和Py-GC/MS表征了催化對苯二甲酸脫羧反應(yīng)前后的ZnO催化劑及對苯二甲酸鋅配合物。同時還對對苯二甲酸鋅配合物進(jìn)行了TG-DTA及DTG分析,以探討ZnO催化對苯二甲酸脫羧反應(yīng)機(jī)理以及催化劑積炭的原因。結(jié)果表明,當(dāng)反應(yīng)溫度為550℃、質(zhì)量空速為0.48 h-1時,ZnO催化對苯二甲酸脫羧的轉(zhuǎn)化率達(dá)到100%,產(chǎn)物中苯的質(zhì)量分?jǐn)?shù)達(dá)40%。ZnO催化對苯二甲酸脫羧反應(yīng)分兩步進(jìn)行,首先對苯二甲酸化學(xué)吸附在ZnO表面形成對苯二甲酸鋅配合物,然后再熱分解生成苯、CO2等。對苯二甲酸鋅配合物的熱分解溫度在410~530℃,在N2中熱分解時發(fā)生的深度脫氫反應(yīng)導(dǎo)致積炭的形成。
對苯二甲酸(TPA);氧化鋅(ZnO);脫羧;苯
精對苯二甲酸(PTA)生產(chǎn)過程中,氧化單元排放的固體廢棄物中含有大量的對苯二甲酸(TPA)等芳香酸,及少量的有色雜質(zhì)和鈷、錳金屬催化劑等,其含量因生產(chǎn)工藝的不同而改變[1],導(dǎo)致其回收利用比較困難。通常采用焚燒、掩埋和生物降解等方法處理PTA殘渣,但容易造成資源浪費和二次污染。目前,有研究者將其應(yīng)用于生產(chǎn)增塑劑、不飽和聚酯樹脂和聚酯涂料等,但其中的有色雜質(zhì)和苯甲酸嚴(yán)重影響了產(chǎn)品的色澤和質(zhì)量[2]。
筆者首次提出PTA殘渣催化脫羧的處理方法,即將其中的芳香酸通過催化脫羧反應(yīng)制備苯、甲苯等芳烴。該方法得到的產(chǎn)物簡單,易分離精制,可以實現(xiàn)資源有效利用和環(huán)境保護(hù)。然而,用于該反應(yīng)的多相催化劑卻不多,目前僅有某些堿土金屬氧化物[3]和過渡金屬氧化物[4-8]可以催化芳香酸脫羧。Masuda等[6]研究了FeOOH催化對苯二甲酸脫羧,反應(yīng)產(chǎn)物主要是二氧化碳、二苯甲酮、苯甲酸等,苯的含量較少。Yoshioka等[7-8]考察了對苯二甲酸在CaO上的催化脫羧行為,苯是主要產(chǎn)物,但CaO在反應(yīng)后生成CaCO3,降低了催化脫羧活性,且積炭量在20%以上。現(xiàn)有的脫羧催化劑存在苯收率低、副產(chǎn)物多以及積炭量偏高的問題。
筆者考察了PTA殘渣中主要組分對苯二甲酸(TPA)在不同金屬氧化物催化下的脫羧行為,并采用FT-IR、XRD、Py-GC/MS等方法表征催化劑,討論了TPA的脫羧反應(yīng)機(jī)理。
1.1 原料
Al(NO)3·9H2O、Zn(NO)2·6H2O、Mg(NO)2·6H2O、Ce(NO)3·6H2O、吡啶、對苯二甲酸(TPA),分析純,國藥集團(tuán)化學(xué)試劑有限公司產(chǎn)品。
1.2 金屬氧化物催化劑的制備
配制一定濃度的硝酸鹽溶液,于60℃恒溫下逐滴加入氨水調(diào)節(jié)溶液pH值,至出現(xiàn)糊狀沉淀。過濾,用純凈水洗滌至中性。經(jīng)110℃干燥,600℃焙燒4 h,研磨篩分,得所需粒徑的金屬氧化物催化劑顆粒。
1.3 對苯二甲酸鋅配合物的制備
稱取1.5 gTPA和4 g Zn(NO)2·6 H2O溶于15 mL吡啶中,在110℃下加熱回流4 h。生成的白色沉淀物經(jīng)過濾、乙醇洗滌后, 60℃真空干燥至恒重, 280℃焙燒2 h,制得白色對苯二甲酸鋅配合物(TP-Zn)。
1.4 催化劑的表征
采用島津公司XRD-6000型X射線衍射分析儀對樣品進(jìn)行XRD分析,管電壓40 kV,管電流80 mA,CuKα射線,連續(xù)掃描,2θ掃描范圍5°~90°。采用Nicolet 公司5700型紅外光譜儀對樣品進(jìn)行FT-IR分析,KBr壓片,波數(shù)范圍400~4000 cm-1,分辨率4 cm-1,掃描次數(shù)32次。采用METTLER公司 TGA/SDTA型熱分析儀對TP-Zn樣品進(jìn)行TG分析,升溫速率15℃/min。采用Thermo Fisher Scientific公司CDS 5200/DSQ II型裂解色質(zhì)聯(lián)用儀對TP-Zn樣品進(jìn)行Pyrolysis-GC-MS分析。稱取一定量TP-Zn樣品置于微型石英裂解管中,裂解溫度700℃,保持10 s。色譜柱為DB-WAX (30 m×0.25 mm×0.25 μm),載氣He,分流比50, 色譜柱程序升溫,從50℃(5 min)→20℃/min→260℃(15 min)。質(zhì)譜分析采用EI離子源,電子能量70 eV,接口溫度250℃,相對分子質(zhì)量掃描范圍25~500。
1.5 催化劑評價
采用固定床反應(yīng)器(管直徑25 mm×2.5 mm、管長1000 mm)評價催化劑催化脫羧性能,實驗裝置如圖1所示。催化劑粒徑40~80目,催化劑質(zhì)量10 g,TPA總進(jìn)料量4.8 g。將質(zhì)量濃度為0.08 g/mL的TPA吡啶溶液,以一定的流量(MHSV=0.48 h-1)泵入固定床反應(yīng)器中,在N2中于一定的溫度下進(jìn)行脫羧反應(yīng)。生成物經(jīng)冷卻,氣、液分離后分別收集。反應(yīng)結(jié)束后,在N2吹掃下將反應(yīng)器溫度降至室溫,取出催化劑進(jìn)行表征。采用在線氣相色譜儀(GC-TCD)分析氣體產(chǎn)物組成,色譜柱為自制聚合物填充柱,柱溫70℃,橋電流120 mA。采用離線氣相色譜儀(GC-FID)分析液體產(chǎn)物組成,DB-5HT毛細(xì)管色譜柱(30 m×0.32 mm×0.10 μm),色譜柱箱程序升溫,從50℃(5 min)→5℃/min→70℃(1 min)→25℃/min→250℃(3 min)。
圖1 管式固定床脫羧反應(yīng)實驗裝置
2.1 不同金屬氧化物的催化脫羧性能
為了消除吡啶對TPA脫羧反應(yīng)的影響,以玻璃珠為惰性物質(zhì)代替催化劑進(jìn)行空白實驗。4種金屬氧化物催化TPA脫羧反應(yīng)的產(chǎn)物分布示于圖2。由圖2可以看出,在500℃時,TPA的熱穩(wěn)定性較高,空白實驗的轉(zhuǎn)化率僅為13.3%。而在ZnO、MgO、Al2O3、CeO2催化下,TPA的轉(zhuǎn)化率有不同程度的提高,其中ZnO的催化脫羧活性和苯選擇性最高,產(chǎn)物中苯的質(zhì)量分?jǐn)?shù)達(dá)到35%。加之價廉易得,ZnO將具有廣闊的應(yīng)用前景。
圖2 500℃時4種金屬氧化物催化TPA脫羧反應(yīng)的產(chǎn)物分布
圖3為TPA脫羧反應(yīng)空白實驗中TPA轉(zhuǎn)化率和苯收率隨反應(yīng)溫度的變化。由圖3可知,在無金屬氧化物催化劑時,TPA脫羧反應(yīng)的TPA轉(zhuǎn)化率和苯收率都隨著反應(yīng)溫度的升高而不斷增加。在TPA熱分解過程中,吡啶中的N原子可能與羧酸中的羰基結(jié)合形成中間過渡態(tài)結(jié)構(gòu)[9-10],從而起到催化脫羧作用。因此,TPA的熱分解可認(rèn)為是溫度和吡啶共同作用的結(jié)果。
圖3 TPA脫羧反應(yīng)空白實驗中反應(yīng)溫度對TPA轉(zhuǎn)化率和苯產(chǎn)率的影響
2.2 反應(yīng)溫度對ZnO催化TPA脫羧反應(yīng)的影響
理論上TPA在完全轉(zhuǎn)化時會生成質(zhì)量分?jǐn)?shù)為47%的苯和53%的CO2。圖4為反應(yīng)溫度對ZnO催化TPA脫羧反應(yīng)產(chǎn)物分布的影響。由圖4可知,隨著反應(yīng)溫度的升高,TPA脫羧反應(yīng)的轉(zhuǎn)化率不斷增加,產(chǎn)物中苯的質(zhì)量分?jǐn)?shù)也逐漸增加,由450℃的19.34%增至550℃的41.69%;CO2的質(zhì)量分?jǐn)?shù)也呈增加的趨勢,由450℃的28.72%增至550℃的52.36%;而殘余TPA質(zhì)量分?jǐn)?shù)逐漸降低,由450℃的25.74%降至550℃的1.85%。反應(yīng)產(chǎn)物中均含有苯甲酸,其含量隨反應(yīng)溫度的升高先增加后降低。因此,在550℃下,ZnO催化TPA脫羧反應(yīng)的活性最高,且副產(chǎn)物含量少。
隨著反應(yīng)溫度的升高,TPA脫羧產(chǎn)物中苯、苯甲酸含量先同時增加,隨后苯甲酸含量逐漸降低,而苯含量卻繼續(xù)增加。由此可以說明,苯甲酸是TPA脫羧反應(yīng)的中間生成物,它可以進(jìn)一步脫羧生成苯??梢哉J(rèn)為,TPA的催化脫羧反應(yīng)同時存在串聯(lián)反應(yīng)和一步反應(yīng)。串聯(lián)反應(yīng)即TPA首先脫去1個羧基生成苯甲酸,然后苯甲酸繼續(xù)催化脫羧生成苯和CO2等產(chǎn)物;一步反應(yīng)即TPA直接催化脫羧生成苯和CO2等產(chǎn)物。
圖4 反應(yīng)溫度對ZnO催化TPA脫羧反應(yīng)產(chǎn)物分布的影響
2.3 對苯二甲酸催化脫羧反應(yīng)機(jī)理
對脫羧反應(yīng)機(jī)理的深入研究,有助于解決催化劑的積炭問題和減少副反應(yīng)的發(fā)生,并對催化劑配方的優(yōu)選提供指導(dǎo)。一般認(rèn)為,金屬氧化物催化乙酸脫羰基制丙酮[11-12]反應(yīng)過程中,乙酸首先吸附在金屬氧化物表面生成羧酸鹽中間物,然后經(jīng)脫水、烯酮化等歷程生成丙酮。金屬氧化物催化TPA脫羧反應(yīng)過程中也有可能生成羧酸鹽。因此,筆者采用FT-IR、XRD、Py-GC/MS表征了催化TPA脫羧反應(yīng)前后的ZnO催化劑和對苯二甲酸鋅配合物TP-Zn。
圖5為TPA、TP-Zn和450℃催化TPA脫羧反應(yīng)后ZnO的FT-IR譜。由圖5可看出,TPA在1678 cm-1處有C=O振動引起的吸收峰和在2500~3000 cm-1范圍內(nèi)由O—H振動引起的較寬的吸收峰。TP-Zn在以上位置無明顯的吸收峰,而在1591 cm-1和1386 cm-1處出現(xiàn)了2個新的吸收峰[13]。羧酸鹽的羧基—COO-具有多電子π鍵體系,2個C—O振動頻率相近,有強(qiáng)烈的振動偶合作用,結(jié)果導(dǎo)致1678 cm-1處C=O 的振動峰消失,在1591 cm-1和1386 cm-1處出現(xiàn)2個新的吸收峰[13]。以上結(jié)果表明,—COOH上的氫解離后以羧酸根的形式存在。TP-Zn在3200~3600 cm-1范圍內(nèi)較寬的吸收峰[14]則是由于其表面存在吸附水而引起。450℃催化TPA脫羧反應(yīng)后的ZnO的紅外吸收峰位置與TP-Zn的特征峰位置相同,表明TPA可能在ZnO表面被化學(xué)吸附生成了TP-Zn。
圖5 TPA、TP-Zn和450℃催化TPA脫羧反應(yīng)后ZnO的FT-IR譜
圖6為TP-Zn和450℃催化TPA脫羧反應(yīng)前后ZnO的XRD譜。由圖6可以看出,450℃催化TPA脫羧反應(yīng)后的ZnO與催化反應(yīng)前ZnO相比較,在2θ為10.1°、16.4°、18.0°、23.5°、30.1°處出現(xiàn)了新的衍射峰。這些衍射峰的位置與TP-Zn的特征衍射峰相同。此結(jié)果進(jìn)一步表明,TPA在脫羧反應(yīng)過程中在ZnO表面生成了TP-Zn。圖7為不同溫度催化TPA脫羧反應(yīng)后ZnO與新鮮ZnO的XRD譜。對比圖6和圖7可以看出,500℃反應(yīng)后ZnO的XRD譜中并沒有TP-Zn的特征峰,其衍射峰的位置與純ZnO一致,表明在500℃反應(yīng)時ZnO表面的TP-Zn分解完全,且沒有新的晶型出現(xiàn)。
圖6 TP-Zn和450℃催化TPA脫羧反應(yīng)前后ZnO的XRD譜
圖7 不同溫度催化TPA脫羧反應(yīng)后ZnO與新鮮ZnO的XRD譜
圖8為TP-Zn的Pyrolysis-GC-MS總離子流圖。利用標(biāo)準(zhǔn)質(zhì)譜庫檢索的信息對主要的13個峰所對應(yīng)的裂解產(chǎn)物進(jìn)行定性分析,所得定性、定量結(jié)果列于表1。從圖8和表1可見,TP-Zn裂解產(chǎn)物中苯的含量占絕對優(yōu)勢,其次是吡啶、二苯甲酮、聯(lián)苯等。吡啶是TP-Zn樣品制備過程中被吸附的溶劑。ZnO催化TPA脫羧得到的產(chǎn)物中主要是苯和CO2,此外還含有少量的二苯甲酮、聯(lián)苯等,與TP-Zn裂解產(chǎn)物相同。因此,TP-Zn是對苯二甲酸在ZnO催化下脫羧反應(yīng)過程中的中間生成物。
由FT-IR、XRD和Pyrolysis-GC-MS分析結(jié)果可以推出TPA的脫羧反應(yīng)機(jī)理,如圖9所示。即TPA首先化學(xué)吸附在ZnO表面生成TP-Zn,然后TP-Zn在高溫下熱分解生成苯、CO2等產(chǎn)物。
圖8 TP-Zn的Pyrolysis-GC-MS總離子流色譜
表1 TP-Zn的Pyrolysis-GC-MS分析結(jié)果
Table 1 Analysis results from Pyrolysis-GC-MS of TP-Zn
PeakNo 1)Retentiontime/minCompoundPercentoftotalarea/%13 50Benzene66 4324 45Pyridine15 5139 19Benzaldehyde0 38412 98Benzoicacid0 25516 13Biphenyl1 03616 83Diphenylmethane0 67718 77Fluorene0 31819 38Benzophenone11 10920 609?Fluorenone0 831021 16Phenanthrene0 101123 049,10?Anthracenedione0 581223 32Triphenylmethane0 101324 649?Phenylfluorene0 40
1) Shown in Fig.8
圖9 TPA在ZnO催化下的脫羧反應(yīng)機(jī)理
2.4 ZnO催化TPA脫羧反應(yīng)積炭形成的原因
TP-Zn的熱分解是TPA 脫羧反應(yīng)過程中的重要環(huán)節(jié),其熱分解行為與積炭的形成密切相關(guān)。圖10為TP-Zn在空氣和N2中的TG-DTA曲線,圖11為相應(yīng)的DTG曲線。由圖10可發(fā)現(xiàn),TP-Zn在空氣和N2中熱分解的趨勢有所區(qū)別。在空氣中,TG曲線有3個失重階段,DTA曲線存在2個較弱的吸熱峰和1個強(qiáng)放熱峰。DTA曲線在203.92℃處的弱吸熱峰對應(yīng)的失重,歸屬于樣品骨架中束縛水及溶劑的脫除;526.03℃處的強(qiáng)放熱峰對應(yīng)的失重則是有機(jī)配體在空氣中燃燒反應(yīng)的結(jié)果。在N2中,TG曲線有4個失重階段,DTA曲線存在3個較弱的吸熱峰和2個弱放熱峰。DTA曲線在195.3℃處的弱吸熱峰對應(yīng)的失重,歸屬于樣品骨架中束縛水及溶劑吡啶的脫除;530.94℃處的弱放熱峰對應(yīng)的失重則是TP-Zn熱分解反應(yīng)的結(jié)果。TP-Zn在空氣和N2中的失重量不同。TP-Zn與空氣中的O2反應(yīng),最終生成ZnO,530℃后質(zhì)量達(dá)到穩(wěn)定,失重量67.3%;程序升溫至530℃時,TP-Zn在N2中的失重量59.9%,明顯低于前者。TP-Zn裂解產(chǎn)物中含有少量的芴、9,10-蒽醌、菲等(見表1),可見,在N2中,TP-Zn同時進(jìn)行了熱分解反應(yīng)和深度脫氫反應(yīng),最終生成C/ZnO復(fù)合產(chǎn)物。TP-Zn深度脫氫反應(yīng)生成的焦炭與ZnO可能發(fā)生氧化還原反應(yīng),導(dǎo)致TP-Zn在N2中程序升溫至1000℃時的失重量(68.7%)高于其在空氣中的失重量(67.4%)。
圖10 TP-Zn在空氣和N2中的TG-DTA曲線
由圖11可以看出,TP-Zn在300℃時,溶劑和水已基本完全脫除;TP-Zn在空氣和N2中的熱解溫度均在410~530℃之間;隨著溫度的升高,TP-Zn的熱分解速率加快,且在500℃左右失重速率達(dá)到最大值。以ZnO為催化劑在較低溫度,如450℃下催化TPA脫羧反應(yīng)時, 其中間產(chǎn)物TP-Zn未分解完全,因此而導(dǎo)致ZnO催化劑的表觀積炭量偏高(見圖4)。
圖11 TP-Zn在空氣和N2中的DTG曲線
(1) ZnO具有較好的催化TPA脫羧制苯性能,脫羧反應(yīng)主要產(chǎn)物是苯和CO2。在反應(yīng)溫度550℃、常壓和質(zhì)量空速0.48 h-1條件下,TPA轉(zhuǎn)化率高達(dá)100%,產(chǎn)物中苯的質(zhì)量分?jǐn)?shù)可達(dá)40%,積炭量約為2%。
(2) ZnO催化TPA的脫羧反應(yīng)歷程分為兩步,首先TPA在ZnO表面被化學(xué)吸附生成對苯二甲酸鋅配合物,然后繼續(xù)熱分解生成苯、二氧化碳等。
(3) 對苯二甲酸鋅配合物在空氣和N2中的熱分解溫度均在410~530℃之間,其DTG峰對應(yīng)的溫度為500℃。其在N2中熱分解時,還能發(fā)生深度脫氫反應(yīng),導(dǎo)致催化劑積炭。
[1] DIESSEL K, MODIC R, STRUSS F. Recovery and reuse of heavy-metal oxidation catalyst from the Witten DMT process:US, 4410449[P].1983-10-18.
[2] ROFFIA P, CALINI P, MOTTA L, et al. Byproduct identification in the terephthalic acid production process and possible mechanisms of their formation[J]. Industrial & Engineering Chemistry Product Research and Development, 1984, 23(4):629-634.
[3] ZHANG A, MA Q,WANG K, et al. Naphthenic acid removal from crude oil through catalytic decarboxylation on magnesium oxide[J]. Applied Catalysis A:General, 2006, 303(1):103-109.
[4] DE LANGE M W, VAN OMMEN J G, LEFFERTS L. Deoxygenation of benzoic acid on metal oxides 1 The selective pathway to benzaldehyde[J]. Applied Catalysis A:General, 2001, 220(1-2):41-49.
[5] DE LANGE M W, VAN OMMEN J G, LEFFERTS L. Deoxygenation of benzoic acid on metal oxides 2 Formation of byproducts[J]. Applied Catalysis A:General, 2002, 231(1):17-26.
[6] MASUDA T, MIWA Y, HASHIMOTO K, et al. Recovery of oil from waste poly(ethylene terephthalate) without producing any sublimate materials[J]. Polymer Degradation and Stability, 1998, 61(2):217-224.
[7] KUMAGAI S, GRAUSE G, KAMEDA T, et al. Decomposition of gaseous terephthalic acid in the presence of CaO[J]. Industrial & Engineering Chemistry Research, 2011, 50(4):1831-1836.
[8] GRAUSE G, HANDA T, KAMEDA T, et al. Effect of temperature management on the hydrolytic degradation of PET in a calcium oxide filled tube reactor[J]. Chemical Engineering Journal, 2011, 166(2):523-528.
[9] ARTOK L, SCHOBERT H H. Reaction of carboxylic acids under coal liquefaction conditions 1 Under nitrogen atmosphere[J]. Journal of Analytical and Applied Pyrolysis, 2000, 54(1):215-233.
[10] CLARK L W. The mechanism of the decomposition of trichloroacetic acid in aromatic amines[J]. The Journal of Physical Chemistry, 1959, 63(1): 99-101.
[11] RAJADURAI S. Pathways for carboxylic acid decomposition on transition metal oxides[J]. Catalysis Reviews, 1994, 36(3): 385-403.
[12] VOHS J M, BARTEAU M A. Reaction pathways and intermediates in the decomposition of acetic and propionic acids on the polar surfaces of zinc oxide[J]. Surface Science, 1988, 201(3):481-502.
[14] BRZYSKA W, WOODKIEWICZ W. Thermal decomposition of copper (II) benzenedicarboxylates[J]. Journal of Thermal Analysis and Calorimetry, 1988, 34(5):1207-1215.
Catalytic Decarboxylation Mechanism of Terephthalic Acid to Benzene Over ZnO Catalyst
ZHUANG Yan, JIANG Binbo, WANG Jingdai, YANG Yongrong
(StateKeyLaboratoryofChemicalEngineering,ZhejiangUniversity,Hangzhou310027,China)
Several metal oxides catalysts were prepared for catalytic decarboxylation of terephthalic acid(TPA), and then ZnO was selected after comparison of their catalytic activities and reaction products distribution of TPA decarboxylation. FT-IR, XRD and Py-GC/MS were used to characterize the ZnO samples before and after catalytic decarboxylation and zinc terephthalate (TP-Zn), which was also analyzed by TG-DTA and DTG, to study the mechanism of the decarboxylation and coke formation. The results showed that conversion of TPA decarboxylation over ZnO was about 100% and the content of benzene in products was over 40% at 550℃ and MHSV of 0.48 h-1. TP-Zn intermediate was identified on ZnO after its catalyzing decarboxylation. The TPA decarboxylation over ZnO involved two steps, in which the first step was chemical adsorption of TPA on ZnO to form TP-Zn and the second step was thermal decomposition of TP-Zn to produce benzene and CO2. The decomposition temperature of TP-Zn was between 410℃ and 530℃ in N2atmosphere. The decomposition companied with deep dehydrogenation of TP-Zn resulted in carbon deposition on ZnO surface.
terephthalic acid(TPA);zinc oxide(ZnO);decarboxylation;benzene
2013-12-17
國家自然科學(xué)基金項目(21176208)和國家高技術(shù)研究發(fā)展計劃項目(2012AA030304)資助
莊巖,男,碩士研究生,從事工業(yè)催化研究;E-mail:zhuangzi7295@126.com
蔣斌波,男,副教授,博士,從事多相流反應(yīng)工程、多相流檢測與信息處理研究;Tel:0571-87952254; E-mail:jiangbb@zju.edu.cn
1001-8719(2015)03-0698-07
TQ032
A
10.3969/j.issn.1001-8719.2015.03.013