胡葉立，米 鐵

（江漢大學(xué)工業(yè)煙塵污染控制湖北省重點(diǎn)實(shí)驗(yàn)室，武漢430056）

活性炭是一種孔隙結(jié)構(gòu)發(fā)達(dá)、比表面積大、吸附能力強(qiáng)的含碳材料，廣泛應(yīng)用于化工、環(huán)保、醫(yī)學(xué)、軍事等領(lǐng)域［1］.活性炭制備原料來(lái)源廣泛，幾乎所有含碳材料均可用于制備活性炭.常見的活性炭制備原料有植物類［2－6］、煤炭類［7］、石油類［8］等.制備原料不同，活性炭制備工藝條件、產(chǎn)品價(jià)格和吸附性能等也有很大的差別.因此，在選擇原料時(shí)，不僅要考慮活性炭的品位、制備工藝的難易、原料成本及來(lái)源是否豐富等因素，同時(shí)也要考慮到資源與環(huán)境問(wèn)題.

我國(guó)是中醫(yī)藥大國(guó)，中藥材多為植物類藥材，在中藥的加工過(guò)程中，產(chǎn)生大量中藥材廢渣.目前，國(guó)內(nèi)的大多數(shù)中藥材廠的大部分藥廢渣均作為廢棄物丟棄，主要處理方式為填埋和作為固體廢物任意排放.這種處理方式不僅造成嚴(yán)重的環(huán)境污染，而且還造成不必要的資源浪費(fèi).如能變廢為寶，既解決了浪費(fèi)和污染，也可帶來(lái)一定的經(jīng)濟(jì)效益.采用多種途徑回收利用中藥材廢渣將成為關(guān)注熱點(diǎn)［9］.將中藥材廢渣制成性能優(yōu)良活性炭，不僅是對(duì)其本身的一種環(huán)?；厥绽茫芡苿?dòng)其他的環(huán)保產(chǎn)業(yè)發(fā)展，是順應(yīng)社會(huì)可持續(xù)發(fā)展的正確途徑.

本文以化學(xué)（KOH）活化法制備活性炭，活性炭的制備方法主要有物理活化法和化學(xué)活化法等，而化學(xué)活化法具有活化溫度低，活化時(shí)間短，活性炭吸附性能好等優(yōu)點(diǎn).在化學(xué)活化中，常見的活化劑H3PO4、KOH、ZnCl2等，其中KOH活化所制備得到的活性炭具有孔隙結(jié)構(gòu)發(fā)達(dá)、比表面積大、吸附性能好等優(yōu)勢(shì)，已應(yīng)用于多種生物質(zhì)原料活化［10－12］，在活性炭制備過(guò)程中所面臨的一大挑戰(zhàn)是針對(duì)某一領(lǐng)域的應(yīng)用制備得到的活性炭應(yīng)具有特定的性能.活性炭最重要的性能就是它的吸附性能，而吸附性能較大程度上取決于制備條件［13］.因此，對(duì)于影響活性炭吸附性能因素的評(píng)價(jià)，合理的試驗(yàn)設(shè)計(jì)是相當(dāng)重要的.響應(yīng)曲面分析法（response surface methodology，RSM）被認(rèn)為是研究?jī)梢蛩鼗蚨鄠€(gè)因素之間交互作用的一種有效工具.RSM采用多元二次回歸方程來(lái)擬合因素和響應(yīng)值之間的函數(shù)關(guān)系，通過(guò)對(duì)回歸方程的分析來(lái)尋求最優(yōu)工藝參數(shù)，解決多變量問(wèn)題的一種統(tǒng)計(jì)方法［14］.相比傳統(tǒng)的正交設(shè)計(jì)等方法，RSM以最經(jīng)濟(jì)的方式、較少的實(shí)驗(yàn)次數(shù)和較短的時(shí)間對(duì)所選的實(shí)驗(yàn)參數(shù)進(jìn)行全面研究，得出正確結(jié)論，并從圖形方面分析尋求最優(yōu)實(shí)驗(yàn)考察因素值，其優(yōu)越性已為越來(lái)越多的實(shí)驗(yàn)工作者所關(guān)注，目前應(yīng)用到多個(gè)領(lǐng)域［15－16］，然而將RSM用于優(yōu)化中藥材廢渣基活性炭制備條件的相關(guān)研究在國(guó)內(nèi)外還尚未見報(bào)道.

1 實(shí)驗(yàn)材料與方法

1.1 實(shí)驗(yàn)原料

本實(shí)驗(yàn)選取的中藥材廢渣為蒸煮過(guò)后的貫眾，將其晾曬、磨碎、過(guò)篩（篩號(hào)20目，孔徑0.9 mm）及干燥后密封存放備用.

1.2 活性炭制備

稱取一定量的中藥材廢渣直接與固體氫氧化鉀按照浸漬比混合，加適量去離子水，攪拌均勻，浸漬24 h，在鼓風(fēng)干燥箱中105℃下干24 h，研磨，置于管式反應(yīng)爐中，在40 m L／min N2保護(hù)下，升溫至600～800℃，并恒溫活化0.5～2 h，在N2保護(hù)下，冷卻至室溫，依次用1 mol／L HCl溶液、煮沸的去離子水、冷去離子水洗滌，洗至溶液p H為7.0左右，干燥得到活性炭.

1.3 活性炭表征

活性炭的吸附性能采用碘吸附值和亞甲基藍(lán)吸附值進(jìn)行表征，活性炭碘吸附值按GB／T12496.8－1999測(cè)定，亞甲基藍(lán)吸附值按GB／T12496.10－1999測(cè)定；活性炭的表面形貌及其熱穩(wěn)定性分別采用掃描電鏡（SEM）、熱重（TGA）進(jìn)行測(cè)定.

1.4 實(shí)驗(yàn)設(shè)計(jì)

采用中心組合設(shè)計(jì)（CCD）響應(yīng)曲面分析方法，選取對(duì)活性炭碘吸附值、亞甲基藍(lán)吸附值有顯著影響的3個(gè)因素：浸漬比（x1）、活化溫度（x2）、活化時(shí)間（x3），分別以A、B、C表示，進(jìn)行三因素三水平共20個(gè)實(shí)驗(yàn)點(diǎn)的響應(yīng)面分析實(shí)驗(yàn).使用Design－Expert 7.1.6進(jìn)行數(shù)據(jù)分析，以碘吸附值（Y1）、亞甲基藍(lán)吸附值（Y2）為響應(yīng)值，求出二次回歸模型，進(jìn)而得到最佳制備工藝條件.試驗(yàn)因子及水平見表1.設(shè)該模型通過(guò)最小二乘法擬合的二次多項(xiàng)方程為：

表1 中心組合設(shè)計(jì)的各獨(dú)立變量以及其代碼水平值Tab.1 Independent variables and their coded levels for the CCD

其中，Y表示響應(yīng)值，在本實(shí)驗(yàn)中響應(yīng)值有碘吸附值（Y1）和亞甲基藍(lán)吸附值（Y2），b0為系數(shù)常數(shù)，bi為線性系數(shù)，bij為交互作用系數(shù)，bii為二次系數(shù)，xi、xj為因素變量值.

2 實(shí)驗(yàn)結(jié)果與討論

2.1 建立回歸模型

CCD用于建立活性炭制備因素對(duì)碘吸附值和亞甲基藍(lán)吸附值之間的相互影響關(guān)系.表2列舉了整個(gè)設(shè)計(jì)列表以及根據(jù)實(shí)驗(yàn)得到的響應(yīng)值，中心點(diǎn)15～20實(shí)驗(yàn)是用于檢驗(yàn)誤差以及數(shù)據(jù)可重復(fù)性.

表2 活性炭制備試驗(yàn)設(shè)計(jì)列表Tab.2 Experimental design matrix for preparation of activated carbons

從表2中可以看出碘吸附值和亞甲基藍(lán)吸附值變化范圍分別為417.51 mg／g～957.22 mg／g和91.16 mg／g～394.22 mg／g.鑒于順序式模型偏差平方和，根據(jù)多項(xiàng)式最高階選擇模型，其附加形式應(yīng)有意義且該模型沒(méi)有其他表述形式.設(shè)計(jì)軟件所推薦碘吸附值和亞甲基藍(lán)吸附值模型為二次模型.通過(guò)多元回歸模型分析，使用二次多項(xiàng)式探究響應(yīng)值與3個(gè)變量之間的相互關(guān)系，如公式（1）所示.對(duì)于碘吸附值（Y1）和亞甲基藍(lán)吸附值（Y2）的最終二次回歸模型以因素代碼的形式表述于公式（2）和公式（3）.

單因素相關(guān)系數(shù)表示實(shí)際影響因素的影響程度，兩因素以及二階形式的相關(guān)系數(shù)分別表示兩因素相互影響程度和二次影響程度，系數(shù)前面為正值說(shuō)明是協(xié)同效應(yīng)，反之為拮抗效應(yīng).從公式（2）中可知浸漬比與活化溫度以及浸漬比與活化時(shí)間對(duì)碘吸附值相互影響均呈現(xiàn)出協(xié)同效應(yīng)，而活化溫度與活化時(shí)間對(duì)碘吸附值相互影響則表現(xiàn)出的是拮抗效應(yīng).同樣從公式（3）可以得到漬比與活化溫度、浸漬比與活化時(shí)間以及活化溫度與活化時(shí)間對(duì)亞甲基藍(lán)吸附值相互影響均呈現(xiàn)出協(xié)同效應(yīng).

一個(gè)多項(xiàng)式模型的建立通常是由相關(guān)系數(shù)R2評(píng)價(jià)其效果，相關(guān)系數(shù)R2越大，則表明二次模型越能滿足實(shí)驗(yàn)數(shù)據(jù)，即該模型擬合效果越好.本實(shí)驗(yàn)中碘吸附值和亞甲基藍(lán)吸附值的R2值分別為0.752和0.757，這說(shuō)明該模型能分別解釋75.2%的碘吸附值的變化和75.7%的亞甲基藍(lán)吸附值的變化［17］.這表明預(yù)測(cè)值是比較精確的，也比較接近實(shí)驗(yàn)值，同時(shí)也表明二次回歸模型能滿足實(shí)驗(yàn)數(shù)據(jù).

2.2 方差分析

模型的建立是否有意義以及穩(wěn)妥性還得進(jìn)一步通過(guò)方程分析進(jìn)行驗(yàn)證.表3、表4分別列出碘吸附值和亞甲基藍(lán)吸附值的方差分析結(jié)果.從表3可以看出，模型F值為3.37，Prob＞F差值為0.036，這表明該模型的建立是有意義的，同時(shí)也表明模型F值的存在僅有3.60%可能性是由于干擾造成的.當(dāng)Prob＞F的差值小于0.05時(shí)，則認(rèn)為該模型是有意義的，當(dāng)差值大于0.1則認(rèn)為該模型無(wú)意義.在此情形下，只有x1、x2對(duì)該模型來(lái)說(shuō)是有意義的，而其它形式對(duì)于響應(yīng)值均無(wú)顯著意義.同樣從表3可以得出這3個(gè)變量對(duì)碘吸附值影響顯著性水平依次為活化溫度＞浸漬比＞活化時(shí)間.

表3 碘吸附值響應(yīng)曲面二次回歸模型方差分析Tab.3 Analysis of variance（ANOVA）for response surface quadratic model for iodine desorption value of activated carbon

表4 亞甲基藍(lán)吸附值響應(yīng)曲面二次回歸模型方差分析Tab.4 Analysis of variance（ANOVA）for response surface quadratic model for methylene blue desorption value of activated carbon

從表4可以看出，模型F值為3.45，Prob＞F差值為0.0333，這表明該模型的建立是有意義的，同時(shí)也表明模型F值的存在僅有3.33%可能性是由于干擾造成的.此處x1x3、x22對(duì)響應(yīng)值有顯著性影響，對(duì)亞甲基藍(lán)吸附值影響顯著性水平依次為浸漬比＞活化溫度＞活化時(shí)間.碘吸附值和亞甲基藍(lán)吸附值的失擬項(xiàng)F值分別為13.95和629 588.34，其Prob＞F值分別為0.058、＜0.0001，這表明碘吸附值和亞甲基藍(lán)吸附值的失擬項(xiàng)均顯著，同時(shí)也說(shuō)明對(duì)于兩個(gè)模型的失擬F值僅有0.58%和0.01%的可能性因干擾發(fā)生，這也可能是因?yàn)榧僭O(shè)模型未考慮一些系統(tǒng)變量而導(dǎo)致失擬值不是很理想［17－18］.

基于數(shù)據(jù)結(jié)果分析，構(gòu)建的響應(yīng)曲面模型對(duì)于預(yù)測(cè)碘吸附值和亞甲基藍(lán)吸附值是可行的，也在研究變量范圍之內(nèi)，該結(jié)果進(jìn)一步由碘吸附值和亞甲基藍(lán)吸附值的預(yù)測(cè)值與實(shí)驗(yàn)值的關(guān)系進(jìn)行驗(yàn)證.圖1和圖2為響應(yīng)值的實(shí)驗(yàn)值與預(yù)測(cè)值的線性擬合圖，其中相關(guān)系數(shù)R2值分別為0.738和0.743，說(shuō)明預(yù)測(cè)值與實(shí)驗(yàn)值相對(duì)比較接近，同時(shí)也說(shuō)明響應(yīng)模型的建立能夠捕捉到活性炭制備變量與響應(yīng)值之間的相關(guān)性.

圖1 碘吸附值的預(yù)測(cè)值與實(shí)驗(yàn)值的關(guān)系Fig.1 Predicted and experimental values of iodine adsorption

圖2 亞甲基藍(lán)吸附值的預(yù)測(cè)值與實(shí)驗(yàn)值的關(guān)系Fig.2 Predicted and experimental values of methylene blue adsorption

2.3 碘吸附值

碘吸附值是活性炭吸附性能的一個(gè)重要指標(biāo)，活性炭制備條件對(duì)活性炭碘吸附值影響較大.從表3可以看出，浸漬比（x1）、活化溫度（x2）對(duì)響應(yīng)值有明顯的影響，而活化溫度F值11.12大于浸漬比F值5.4，這說(shuō)明活化溫度對(duì)響應(yīng)值的影響效果比浸漬比對(duì)響應(yīng)值的影響顯著.通過(guò)顯著模型描繪的響應(yīng)曲面是為了研究中間組分之間的相互關(guān)系.圖3中的三維響應(yīng)曲面圖呈現(xiàn)出的浸漬比與活化溫度及其交互作用對(duì)碘吸附值的影響，其中活化時(shí)間處于0水平（x3＝60 min）.可以看出，碘吸附值隨著浸漬比和活化溫度的增大而增大.浸漬比對(duì)活性炭孔結(jié)構(gòu)的形成以及活性炭比表面積的增大有顯著性的影響，這可能是由于在活化過(guò)程中原料中的碳與氫氧化鉀發(fā)生反應(yīng)，生成K2CO3，同時(shí)不斷產(chǎn)生CO2和CO氣體，而CO2和CO對(duì)活性炭孔結(jié)構(gòu)的形成有一定的促進(jìn)作用.除此之外，當(dāng)活化溫度高于700℃時(shí)，K2CO3開始分解，最終生成金屬鉀，鉀蒸氣擴(kuò)散嵌入不同碳層之間，導(dǎo)致活性炭材料劇烈膨脹，形成新的孔結(jié)構(gòu)，活性炭比表面積增大，活性炭吸附性能增強(qiáng)［18］.

圖3 浸漬比與活化溫度對(duì)碘吸附值的響應(yīng)曲面圖Fig.3 The effects of impregation ratio and activation temperature on iodine adsorption value

圖4 活化溫度和活化時(shí)間對(duì)亞甲基藍(lán)吸附值的響應(yīng)曲面圖Fig.4 The effects of activation temperature and activation time on methylene blue adsorption value

2.4 亞甲基藍(lán)吸附值

亞甲基藍(lán)吸附值是活性炭吸附性能的又一重要指標(biāo)，活性炭制備條件對(duì)活性炭亞甲基藍(lán)吸附值影響也相對(duì)較大.圖4中的三維響應(yīng)曲面圖表示的是活化溫度和活化時(shí)間及其交互作用對(duì)亞甲基藍(lán)吸附值的影響，其中浸漬比處于0水平（x1＝2）.從圖4中可以看出，亞甲基藍(lán)吸附值開始隨著活化溫度和活化時(shí)間的增大而增大，然后當(dāng)活化溫度和活化時(shí)間進(jìn)一步增大時(shí)，亞甲基藍(lán)吸附值逐漸減小.這說(shuō)明這活化溫度和活化時(shí)間及其交互作用對(duì)亞甲基藍(lán)吸附值有一定的影響.隨著活化反應(yīng)的進(jìn)行，活化溫度升高，活化反應(yīng)速率加快，所形成的活性炭孔隙結(jié)構(gòu)發(fā)達(dá)，孔徑變大，活性炭亞甲基藍(lán)吸附性能增強(qiáng)，當(dāng)活化時(shí)間過(guò)長(zhǎng)，活化溫度過(guò)高，將導(dǎo)致活性炭孔徑進(jìn)一步增大，甚至造成原料碳燒蝕，孔塌陷，從而活性炭亞甲基藍(lán)吸附性能降低［19］.

2.5 響應(yīng)曲面優(yōu)化

本實(shí)驗(yàn)的主要研究目的是找出活性炭制備最佳工藝參數(shù)，以致于所制得的活性炭吸附性能最佳.表5給出了通過(guò)響應(yīng)曲面分析所得出的活性炭制備最優(yōu)工藝條件，以及在此條件下所得到的碘吸附值和亞甲基藍(lán)吸附值的預(yù)測(cè)值和實(shí)驗(yàn)值.最優(yōu)條件為浸漬比3、活化溫度744℃、活化時(shí)間75 min.在此條件下制備的活性炭碘吸附值和亞甲基藍(lán)吸附值為723.75 mg／g、350.82 mg／g，相對(duì)誤差分別為5.83%和0.17%，而相關(guān)文獻(xiàn)［13，17－18，20］中響應(yīng)值的實(shí)驗(yàn)值與預(yù)測(cè)值的相對(duì)誤差均在3%左右，有個(gè)別響應(yīng)值相對(duì)誤差超過(guò)5.6%.與這些文獻(xiàn)相比而言，本文碘吸附值的相對(duì)誤差稍微有點(diǎn)偏大，但整體的實(shí)驗(yàn)值與預(yù)測(cè)值基本一致.因此，基于響應(yīng)曲面法所得的最佳工藝參數(shù)準(zhǔn)確可靠，具有實(shí)用價(jià)值.

表5 模型最優(yōu)化驗(yàn)證Tab.5 Optimum model validation

2.6 活性炭的表面形貌

活性炭的內(nèi)部結(jié)構(gòu)決定了其吸附性能，通過(guò)掃描電鏡圖像（SEM）分析給出活性炭的表面形貌可對(duì)活性炭的內(nèi)部結(jié)構(gòu)進(jìn)行定性分析.圖5（a）、圖5（b）分別給出了活性炭前軀體（500℃下，N2氛圍中炭化30 min的炭化料）表面形貌以及在最優(yōu)工藝條件（浸漬比3、活化溫度744℃、活化時(shí)間75 min）下制備得到的活性炭的表面形貌.從圖5（a）、圖5（b）中均可清晰看出，樣品表面光滑，形狀不規(guī)則，主要呈管狀結(jié)構(gòu).圖5（a）中樣品表面未出現(xiàn)明顯孔隙結(jié)構(gòu)，而圖5（b）中樣品表面形成了大量的排列一致的孔結(jié)構(gòu).結(jié)果表明在活化過(guò)程中，KOH有效地促進(jìn)了活性炭孔結(jié)構(gòu)的形成，因而形成了比表面積大、孔隙結(jié)構(gòu)發(fā)達(dá)的活性炭.類似結(jié)果也在一些相關(guān)文獻(xiàn)［21］中報(bào)道過(guò).

圖5 活性炭前軀體（a）和活性炭（b）SEM圖Fig.5 SEM images：（a）Precursor，（b）Activated carbon（magnifications：1200X）

2.7 熱重分析

為了解中藥材廢渣基活性炭的熱穩(wěn)定性，將最優(yōu)條件下活化所制備的活性炭在N2保護(hù)下進(jìn)行熱重分析，TG曲線如圖6所示.從圖6中可以看出，在300℃以下，活性炭樣品重量有輕微下降，在300～800℃，樣品重量損失加快，但幅度不大，整體上活性炭樣品損失量較小，由此可以說(shuō)明活性炭的熱穩(wěn)定性較好.

圖6 活性炭樣品的熱穩(wěn)定性分析Fig.6 The thermal stability analysis of activated carbon sample

3 結(jié)論

1）采用CCD響應(yīng)曲面分析法優(yōu)化了KOH活化中藥材廢渣基活性炭制備條件.

2）根據(jù)響應(yīng)曲面分析，活化溫度對(duì)碘吸附值的影響最為顯著，而對(duì)于亞甲基藍(lán)吸附值，浸漬比對(duì)其影響最為顯著.

3）根據(jù)響應(yīng)曲面分析所得到的最優(yōu)工藝條件為浸漬比3、活化溫度744℃、活化時(shí)間75 min，在此條件下制備的活性炭碘吸附值和亞甲基藍(lán)吸附值分別為723.75 mg／g、350.82 mg／g，與理論預(yù)測(cè)值幾乎一致，說(shuō)明基于響應(yīng)曲面法所得的最佳工藝參數(shù)準(zhǔn)確可靠，具有一定的實(shí)用價(jià)值.

4）通過(guò)SEM、熱重分析可知最優(yōu)條件下所制備的活性炭具有孔隙結(jié)構(gòu)發(fā)達(dá)、熱穩(wěn)定性高、吸附性能較好等特點(diǎn).

［1］ 孫 康，蔣劍春.國(guó)內(nèi)外活性炭的研究進(jìn)展及發(fā)展趨勢(shì)［J］.林產(chǎn)化學(xué)與工業(yè)，2009，29（6）：98－104.

［2］ Issa I Salame，Teresa J Bandosz.Comparison of the surface features of two wood－based activated carbon［J］.Ind Eng Chem Res，2000，39（2）：301－306.

［3］ 李園園，陳少華，張召基，等.KOH活化絲瓜絡(luò)制備高比表面積活性炭［J］.化工進(jìn)展，2012，31（6）：1274－1279.

［4］ Sait Yorgun，Naile Vural，Hakan Demiral.Preparation of high－surface area activated carbons from Paulownia wood by ZnCl2activation［J］.Microporous and Mesoporous Materials，2009，122：189－194.

［5］ 黃明堦，陳衛(wèi)群，陳燕丹，等.草酸鉀活化法制備榴蓮殼活性炭及其表征［J］.環(huán)境工程學(xué)報(bào)，2012，6（10）：3730－3734.

［6］ Tham Y J，Latif P A，AbdullahB A M，et al.Performances of toluene removal by activated carbon derived from durian shell［J］.Bioresource Technology，2011，102（2）：724－728.

［7］ 梁大明.中國(guó)煤質(zhì)活性炭［M］.北京：化學(xué)工業(yè)出版社，2008.

［8］ 鄧梅根，王仁清，馮義紅.石油焦氧化改性制備活性炭及其電容性能研究［J］.功能材料，2013，44（13）：1960－1963.

［9］ 米 鐵，胡葉立，余新明，等.中藥材廢渣的資源化利用的研究［J］.工業(yè)安全與環(huán)保，2013，39（8）：26－28.

［10］ 毛 宇，馬承愚，范藝苑，等.廢棄辣椒秸稈高比表面積活性炭的制備及表征［J］.應(yīng)用化工，2012，41（5）：776－770.

［11］ 陳 永，周柳江，洪玉珍，等.椰殼纖維基高比表面積中孔活性炭的制備［J］.新型炭材料，2010，25（2）：151－155.

［12］ Zaira Zaman Chowdhury，Sharifuddin M Z，Rashid A K，et al.Preparation，characterization and adsorption performance of the KOH－activated carbons derived from kenaf fiber for lead（II）removal from waste water［J］.Scientific Research and Essays，2011，6（29）：6185－6196.

［13］ Tan I A W，Ahmad A L，Hameed B H.Optimization of preparation conditions for activated carbons from coconut husk using response surface methodology［J］.Chemical Engineering Journal，2008，137：462－470.

［14］ Alam M Z，Muyibi S A，Toramae J.Statistical optimization of adsorption processes for removal of 2，4－dichlorophenol by activated carbon derived from oil palm empty fruit bunches［J］.Environ Sci，2007，19：674－677.

［15］ 李家元，吳彥瑜，周少奇.響應(yīng)曲面法優(yōu)化絮凝處理木薯淀粉廢水［J］.環(huán)境工程學(xué)報(bào)，2010，7（4）：1555－1560.

［16］ 謝明勇，萬(wàn) 成，聶少平，等.響應(yīng)面法優(yōu)化低聚果糖液脫色工藝［J］.食品科學(xué)，2011，32（2）：9－12.

［17］ Mohd Azmier Ahmad，Rasyidah Alrozi.Optimization of preparation conditions for mangosteen peel－based activated carbons for the removal of Remazol Brilliant Blue R using response surface methodology［J］.Chemical Engineering Journal，2010，165：883－890.

［18］ Ahmad A A，Hameed B H.Effect of preparation conditions of activated carbon from bamboo waste for real textile wastewater［J］.Journal of Hazardous Materials，2010，273：487－493.

［18］ Lozano－Castello D，Lillo－Rodenas M A，Cazorla－Amoros D，et al.Preparation of activated carbons from spanish anthracite by KOH［J］.Carbon，2001，39：741－749.

［19］ Lua A C，Yang T.Characteristics of activated carbon prepared from pistachio－nut shell by zinc chloride activation under nitrogen and vacuum conditions［J］.Journal of Colloid and Interface Science，2005，290（2）：505－513.

［20］ Salman J M.Optimization of preparation conditions for activated carbon from palm oil fron－ds using response surface methodology on removal of pesticides from aqueous solution［J］.Arabian Journal of Chemistry，2014，7：101－108.

［21］ Hameed B H，Tan I A W，Ahmad A L.Preparation of oil palm empty fruit bunch－based activated carbon for removal of 2，4，6－trichlorophenol：optimization using response surface methodology［J］.Hazard Mater，2009，164：1316－1324.