李小斐，宋坤莉，趙東風，于 瀾，李 石

(中國石油大學(華東) 化學工程學院，山東 青島 266580)

揮發(fā)性有機物(VOCs)治理是近年來改善中國環(huán)境空氣質量的重要內容[1-3]，苯既是VOCs的重要組分，也是石油化工企業(yè)VOCs排放的典型特征污染物。苯不僅自身具有較高的毒性，同時也是環(huán)境空氣中PM2.5和O3的潛在前體物質。石油化工企業(yè)的苯排放主要來自苯儲罐，是石油化工企業(yè)VOCs排放的重點管控對象之一，目前常用的治理技術是利用活性炭吸附對苯蒸氣進行吸附回收[4-7]。雖然活性炭吸附技術具有運行成本低、機械強度高、操作簡單等優(yōu)點，但也存在不足，一方面目前對于活性炭的選擇缺乏理論上的指導，往往過于追求具有超高比表面積的活性炭作為吸附劑，而缺乏對活性炭吸附苯過程微觀內在機理的考慮，例如活性炭表面官能團對吸附過程的影響、苯分子在孔結構中的吸附位置與構型等，造成了活性炭材料選擇上的盲目性和選擇因素的單一性；另一方面，由于苯在活性炭上的吸附過程是放熱反應，吸附過程中會放出大量的熱，易出現“飛溫現象”，不僅會造成活性炭吸附能力的大幅下降，還會帶來較大的安全風險，容易引起吸附罐內無火焰燃燒，即陰燃現象，而目前對于吸附熱的研究較少。

活性炭吸附苯過程的微觀機理和能量變化情況難以通過實驗手段實現，而分子模擬技術可以有效解決這些問題。筆者在構建具有氫基、羥基、羰基和羧基官能團活性炭分子模型的基礎上，利用巨正則蒙特卡羅(Grand Canonical Monte Carlo,GCMC)方法，探索了在303.15 K下活性炭分子模型中不同表面官能團對苯分子吸附過程中吸附等溫線、等量吸附熱、相互作用能和吸附構型等的微觀影響和本質機理，從微觀層面觀察了苯分子在活性炭中的吸附過程，并可直觀地觀察到苯分子的吸附狀態(tài)和吸附位。不僅有效彌補了實驗技術的不足和基礎理論分析的缺陷[6-11]，而且對于活性炭類吸附劑的選擇以及具有更高吸附性能的活性炭研發(fā)均具有重要的指導意義[12-14]。同時對于石油化工企業(yè)其他輕質油品儲罐(如汽油儲罐)揮發(fā)的VOCs吸附治理，提供一定的借鑒和參考。

1 活性炭模型的構建和表征

1.1 帶有不同官能團活性炭模型的構建

首先利用MS軟件的“Volume visualization”模塊構建了組成活性炭微觀結構的基本單元，基本結構單元中包括純炭結構和含有氫基、羥基、羰基、羧基的炭結構，如圖1所示。

Blue ball—Carbon atom；Red ball—Oxygen atom；White ball—Hydrogen atom圖1 活性炭微觀基本結構單元Fig.1 The activated carbon microstructural unit(a)Pure carbon structure；(b)Carbon structure containing hydrogen；(c)Carbon structure containing hydroxyl；(d)Carbon structure with carbonyl；(e)Carbon structure with carboxyl

通過MS中的“Amorphous cell tools”模塊完成活性炭微觀結構模型的構建。將純炭基本結構單元(見圖1(a))填充到“Build”模塊構建的空晶胞中，形成基本活性炭結構模型(如圖2(a)所示)，記為AC。為了分析吸附過程中官能團的影響，同步構建了表面帶有不同官能團的活性炭分子模型(如圖2(b)～(e)所示)，并且不同模型中氫基、羥基、羰基和羧基數量相同，均為模型總原子數的15%，分別記為AC-H、AC-OH、AC-O和AC-COOH。

為了更加清晰地看到所構建活性炭分子模型的內部結構，創(chuàng)建了模型的等位面，如圖2(f)～(j)所示。為了保證構建模型與真實活性炭結構的相符性，模型中元素比例、炭密度等結構參數根據文獻[6,8,12]中報道的真實數據確定，同時利用“Forcite”模塊對模型結構進行優(yōu)化，以保證構建的模型具有最低能量和最優(yōu)構型。

1.2 活性炭模型的表征

活性炭模型的表征主要包括活性炭密度、徑向分布函數(Radial distribution function,RDF)、表面積、孔隙率、孔體積和孔徑范圍等?；钚蕴棵芏瓤梢灾苯訌腗S軟件的“Properties”模塊獲得?；钚蕴吭娱g的RDF曲線通過“Forcite”模塊分析，如圖3所示。由圖3可以看出，AC模型的第1個峰值出現在約0.14 nm位置，而增加官能團后模型的第1個峰值出現在0.1 nm附近，說明官能團的引入使得活性炭的孔徑結構減小，RDF曲線峰值與大多數實驗數據非常吻合[15-16]，表明了構建模型與真實活性炭結構的相似性。

r—Distance from reference particle圖3 不同模型活性炭原子間的徑向分布函數(RDF)曲線Fig.3 Radial distribution function (RDF)curves between carbon atoms of different models

活性炭分子模型的表面積、孔隙率、孔體積和孔徑范圍通過“Atom volumes &Surfaces”模塊來分析。表面積(S，m2/g)、孔體積(Vp，cm3/g)和孔隙率(P，%)分別由式(1)、式(2)和式(3)計算。

(1)

(2)

(3)

式(1)～式(3)中：Sp為每個單元的比表面積，nm2；ρ為密度，g/cm3；V0為晶胞體積，nm3；VF為自由體積，nm3。

構建模型的孔徑范圍以分子探針直徑對孔分數進行微分獲得，如圖4所示。由圖4可以看出，不同活性炭結構模型孔徑范圍相似，主要為0.30～0.85 nm，而且官能團引入后使模型結構的孔徑尺寸有一定程度的縮小，考慮到苯的吸附過程主要在微孔中進行，因此研究中忽略了中、大孔尺寸對吸附過程的影響。

Pd—Distribution probability;Dp—Pore diameter圖4 不同活性炭模型結構的孔徑范圍Fig.4 Pore diameter scope of different model structures

將活性炭分子模型的表征結果與文獻[17-20]報道的活性炭結構特征參數進行了比較，如表1所示。由表1可以看出，筆者所構建的活性炭分子模型參數與實驗結果非常相似，進一步說明了所構建模型的合理性。

表1 活性炭結構模型的表征參數與實驗結果的比較Table 1 Comparison between the characterization parameters of the activated carbon structure model and the experimental results

2 巨正則蒙特卡羅(GCMC)模擬

采用GCMC方法研究了不同模型活性炭對苯分子的吸附行為和性能。GCMC程序包括粒子的轉變、旋轉、破壞和產生，以確保體相和孔相之間潛在的化學平衡，同時在三維方向(x，y，z)設置周期性邊界條件。在模擬開始時，執(zhí)行1×104步平衡操作以將構建模型結構的能量降低到合理范圍，然后進行1×106步計算操作。使用Metropolis算法使吸附苯分子的結構模型系統具有最低能量，力場選擇基于Ewald加和與Atom方法計算的“COMPASS”[21-22]。通過MS軟件中的“Sorption”模塊，對苯在不同活性炭模型結構中吸附等溫線、吸附熱、相互作用能和吸附構型等進行了詳細分析，活性炭結構模型和苯的勢能參數取值如表2所示。

表2 吸附物和吸附質的蘭納.瓊斯(LJ)勢能參數和部分原子電荷[8,23]Table 2 Lennard-Jones (LJ)potential energy parameters and partial atomic charges of adsorbents and adsorbents[8,23]

3 活性炭表面不同官能團對苯吸附的影響分析

3.1 活性炭模型中氫基、羥基、羰基和羧基對苯吸附等溫線的影響

不同活性炭模型在303.15 K下對苯的吸附等溫線如圖5所示。由圖5可以看出，活性炭對苯的吸附容量在0.1 kPa之前呈線性增加，斜率越大說明吸附速率越快。與AC相比，添加氫基和羥基后的AC-H和AC-OH對苯的吸附速率減??；而添加羰基和羧基后，AC-O和AC-COOH對苯的吸附速率增大，且AC-COOH對苯的吸附速率最大。隨著壓力的增加，吸附逐漸達到平衡，不同活性炭模型對苯的吸附量由大到小順序為AC-COOH、AC-O、AC、AC-OH、AC-H。由此可見，氫基和羥基官能團的引入并不利于活性炭對苯分子的吸附，說明氫基和羥基對苯分子具有一定的排斥作用。而羰基和羧基的引入則不同，盡管活性炭的孔徑有一定減小(如圖4所示)，但活性炭對苯分子的吸附速率和吸附量均得到了一定的提高，說明氧的引入為活性炭內部結構提供了更多的吸附位點，并且苯分子是極性分子，易于形成氫鍵作用。因此，氧原子的密度越大，活性炭對苯分子的吸附量越高，并且羧基官能團的影響要大于羰基官能團的影響。羥基官能團雖然也含有氧原子，但由于羥基是堿性官能團，降低了活性炭表面的極性，也降低了活性炭對苯分子的吸附能力。

p—Pressure;Ac—Adsorptive capacity圖5 苯在不同活性炭模型結構中的吸附等溫線Fig.5 Adsorption isotherm of benzene in different activated carbon model structures T=303.15 K

3.2 活性炭模型中氫基、羥基、羰基和羧基對苯等量吸附熱的影響

等量吸附熱是衡量活性炭吸附性能的另一個關鍵參數，它表征了吸附過程中的熱效應。不同活性炭分子模型吸附苯過程的等量吸附熱變化如圖6所示。由圖6可以看出：在壓力小于0.1 kPa的范圍，不同活性炭模型對苯的等量吸附熱明顯減少；然后隨著壓力的增加，等量吸附熱緩慢增加，這與文獻[23-28]報道的結果相一致。與AC相比，AC-O和AC-COOH對苯的等量吸附熱增大，而且羧基官能團對等量吸附熱的影響大于羰基官能團，而AC-H和AC-OH對苯的等量吸附熱減小。說明含氧官能團的引入增強了活性炭與苯分子之間的結合能力，分析原因是羰基和羧基作為AC-O和AC-COOH模型結構中的高能吸附位點，在低壓下即可吸附較多的苯分子，靜電作用起主導作用。隨著壓力的增加，模型結構中的低能吸附位點開始吸附填充苯分子，炭原子和苯分子之間的相互作用占主導地位。而羥基和氫基的引入作用相反，降低了活性炭與苯分子之間的作用力強度，吸附熱較低。

p—Pressure圖6 苯在不同活性炭模型結構中的等量吸附熱曲線Fig.6 Isosteric heat curves of benzene in different activated carbon model structures

圖7為苯在不同活性炭模型結構中的平均等量吸附熱。由圖7可以看出，在壓力0～1 kPa范圍內活性炭對苯的平均吸附熱由大到小順序為AC-COOH、AC-O、AC、AC-OH、AC-H，這與苯在不同活性炭模型結構中的吸附等溫線變化(見圖5)相一致。因此，在實際應用中，表面官能團的組成與性質也應是活性炭選擇上重要的考慮因素。

圖7 壓力0～1 kPa條件下苯在不同活性炭模型結構中的平均等量吸附熱Fig.7 Average isosteric heat of benzene in different activated carbon model structures at 0-1 kPa

3.3 活性炭模型中氫基、羥基、羰基和羧基對苯吸附作用能的影響

圖8為不同活性炭模型結構中苯的吸附作用能的比較曲線。由圖8可以看出，不同活性炭模型在苯吸附過程中均存在分子間作用能、靜電作用能和范德華作用能，并且分子間作用能在整個吸附過程中占主導地位，而靜電作用能比范德華作用能更小。從圖8(a)可以看出,苯和活性炭分子間相互作用能是恒定的，不隨壓力增加而變化，而且由于含氧官能團的引入，有利于模型結構與苯分子間作用能的增強，不同活性炭模型與苯的分子間相互作用能由大到小順序為AC-COOH、AC-O、AC-OH、AC、AC-H。從圖8(b)和(c)可以發(fā)現，靜電作用能和范德華作用能都呈現周期性下降的趨勢，并且變化范圍相似。說明在吸附過程中，苯首先吸附在作用能較高的高能位點，然后吸附在作用能較低的低能位點，這與前面的吸附等溫線和等量吸附熱分析結果相似。

STEP—Number of Grand Canonical Monte Carlo steps圖8 不同活性炭模型結構中苯的吸附作用能曲線的比較Fig.8 Comparison of adsorption energy curves of benzene in different activated carbon model structures(a)Intermolecular energy；(b)Electrostatic energy；(c)Van der Waals energy

3.4 活性炭模型中氫基、羥基、羰基和羧基對苯吸附構型的影響

在303.15 K條件下，分析了苯在不同活性炭分子模型中的吸附過程。圖9為AC、AC-H、AC-OH、AC-O和AC-COOH模型吸附苯分子后的構型圖。由圖9可以清晰地看到苯分子在各個模型結構中的吸附位置和數量，帶有羧基的活性炭模型中吸附的苯分子最多，帶有羰基的活性炭模型次之，這也說明了含氧官能團的引入有利于活性炭對苯分子的吸附，與前述吸附等溫線反映的吸附量結果相一致。

4 結 論

(1)利用MS軟件，構建表面具有氫基、羥基、羰基、羧基官能團等的活性炭模型，并利用GCMC方法考察了活性炭模型對苯吸附的微觀過程和機理。

(2)酸性含氧官能團的引入為活性炭結構引入高能吸附位點，使活性炭與苯分子之間的吸附強度增加，有利于活性炭對苯吸附量和吸附速率的提高，且羧基官能團的效果要優(yōu)于羰基官能團，吸附熱也最大。而羥基自身的堿性性質降低了活性炭表面的極性，使活性炭對苯的吸附性能下降；氫基官能團使炭密度降低，沒有引入更多的吸附位點，使該結構對苯的吸附能力最低。

(3)從活性炭對苯的吸附作用能來看，苯首先吸附在高能位點，然后吸附在作用能較低的低能位點，吸附過程中分子間作用能占主導地位，而靜電作用能比范德華作用能更小。

Red—Adsorbed benzene molecules圖9 不同活性炭模型結構吸附苯的構型圖Fig.9 Configuration diagram of benzene adsorption by different activated carbon model structures(a)AC；(b)AC-H；(c)AC-OH；(d)AC-O；(e)AC-COOH