草炭土及其有機質(zhì)組分對PCB138吸附動力學(xué)和熱力學(xué)研究

叢鑫，畢然，孫思坤

遼寧工程技術(shù)大學(xué)環(huán)境科學(xué)與工程學(xué)院，遼寧 阜新 123000

多氯聯(lián)苯（PCBs）是由209種同類物組成的一組氯代芳烴化合物，其理化性質(zhì)較為穩(wěn)定，在環(huán)境中具有長距離遷移的能力，目前PCBs廣泛分布在各環(huán)境介質(zhì)中（穆熙等，2019）。其中土壤中的PCBs難以分解，殘留時間高于在水體和大氣中，而且能夠通過地下水遷移至地表水中，對食物和飲用水安全造成重大風(fēng)險。由于吸附作用是污染物質(zhì)進入土壤環(huán)境中的第一個環(huán)境行為，因此吸附作用是影響污染物質(zhì)在環(huán)境中遷移轉(zhuǎn)化和最終歸宿的關(guān)鍵過程（崔兆杰等，2010）。有關(guān)PCBs在土壤中吸附的研究已經(jīng)取得了很多成果，一些學(xué)者探討了外源物質(zhì)的加入對PCBs在土壤中吸附行為的影響。如陳迪（2018）探討了一種天然可降解有機物木素磺酸鉀，分析其加入到土壤后可能對土壤吸附 PCBs影響的機理；陳濤等（2012）對不同類型表面活性劑的添加對PCBs在土壤中吸附行為的影響進行了分析。陳雄（2015）從生物學(xué)角度研究了附著劍菌對多氯聯(lián)苯的生物吸附和生物降解影響。王沛然（2015）對不同類型黏土礦物蒙脫石及腐殖質(zhì)-蒙脫石復(fù)合體對PCBs的吸附熱力學(xué)和分子動力學(xué)進行了模擬研究。

也有學(xué)者分析了有機質(zhì)對PCBs在土壤中殘留變化的影響（崔嵩等，2013）。但目前在這方面的研究多見于影響因素的相關(guān)性分析，對于不同有機質(zhì)組分對有機污染物在土壤吸附中所起作用的研究還缺少相關(guān)的數(shù)據(jù)。本研究以此為切入點，取有機質(zhì)含量豐富，北方地區(qū)分布較廣的草炭土為供試土壤，選取城市表層土壤中含量較高的六氯聯(lián)苯PCB138為研究對象，分析草炭土及其不同有機質(zhì)組分對PCB138的吸附特征，確定吸附動力學(xué)和熱力學(xué)模型，計算有關(guān)熱力學(xué)參數(shù)，分析其吸附機理，從而進一步明確多氯聯(lián)苯在土壤中的吸附行為，為PCBs的污染控制和生態(tài)風(fēng)險評價提供一些基礎(chǔ)實驗數(shù)據(jù)。

1 材料與方法

1.1 土壤樣品

供試的草炭土采自遼寧省阜新市彰武縣。土壤樣品采集后除去石塊和石礫，自然風(fēng)干，過2 mm篩后保存。供試土壤的理化性質(zhì)如表1所示。土壤有機質(zhì)測定采用燒失量法；土壤陽離子交換量采用乙酸鈉法；土壤粒徑分布由Horiba公司LA-300激光粒度儀測定。從表1可以看出供試草炭土有機質(zhì)質(zhì)量分數(shù)較大。

表1 供試土壤理化性質(zhì)Table 1 Physicochemical properties of the tested soils

1.2 胡敏酸和胡敏素提取實驗

稱取一定量土壤，按照土與鹽酸1:10的比例加入 0.05 mol·L-1鹽酸，攪拌，浸提，吸取上清液，清洗后再加入0.1 mol·L-1NaOH，間歇攪拌，浸提，上清液加鹽酸調(diào)至 pH為 1.0—1.5，使其產(chǎn)生褐色絮狀沉淀。靜置，離心（3000 r·min-1，15 min），得凝膠狀粗胡敏酸。向凝膠狀粗胡敏酸中加入 0.1 mol·L-1NaOH，間歇攪拌，浸提，上清液加鹽酸調(diào)至pH為1.0—1.5，使其產(chǎn)生褐色絮狀沉淀。靜置，離心，得細胡敏酸，冷凍干燥后研磨，保存?zhèn)溆谩?/p>

將土壤樣品放入容器中，加入質(zhì)量分數(shù)為40%的HF和2 mol·L-1HCl，體積比為1:1?；旌衔镌?0 ℃搖床下震蕩1 h，離心，棄去上清液，重復(fù)上述操作3次以上，洗滌，冷凍干燥得到胡敏素。

1.3 吸附試驗

吸附動力學(xué)實驗：稱取 1.00 g土壤樣品（過0.177 mm 篩）于錐形瓶中，加入 300 μg·L-1的PCB138溶液30 mL，將錐形瓶置于恒溫振蕩器上在室溫（25 ℃）條件下振蕩10、20、30、40、60、90、120、240、360、480 min-1，離心（4000 r·min-1，20 min），移取上清液10 mL于分液漏斗中，用正己烷萃取3次，收集有機相，過無水硫酸鈉柱去除水分，氮吹定容至1 mL，GC-MS檢測。

吸附熱力學(xué)實驗：稱取 1.00 g土壤樣品（過0.177 mm篩）于錐形瓶中，加入60、180、300、400、600、900 μg·L-1的 PCB138 溶液 30 mL，將錐形瓶置于恒溫振蕩器上，溫度分別控制在 (25±1)、(35±1)、(45±1) ℃振蕩 240 min-1，離心（4000 r·min-1，20 min），移取上清液10 mL于分液漏斗中，用正己烷萃取3次，收集有機相，過無水硫酸鈉柱去除水分，氮吹定容至1 mL，GC-MS檢測。

取 300、600、900 μg·g-13 種濃度 PCB138 標準樣品進行空白吸附實驗，實驗結(jié)果表明，回收率為72.8%—83.5%。

1.4 分析方法

采用氣相色譜質(zhì)譜分析方法（QP2010SE，日本島津公司）。色譜操作條件為：Rtx-5ms毛細管柱（30 m×0.25 mm×0.25 μm）；測試條件為：進樣口溫度 280 ℃；離子源溫度 230 ℃；程序升溫為50 ℃保持 1 min，然后以 20 ℃·min-1的速率升至280 ℃，保持3 min。不分流進樣；進樣量1 μL；載氣為高純氦氣（99.999%）。選擇離子掃描方式，外標法定量。

1.5 數(shù)據(jù)處理

吸附動力學(xué)采用準一級、準二級、Elovich和Weber-Morris動力學(xué)模型，計算公式如下：

式中：qe為吸附平衡時吸附量（μg·g-1）；qt為t時刻吸附量（μg·g-1）；k1為準一級吸附速率常數(shù)（min-1）；k2為準二級吸附速率常數(shù)[g·(μg·min)-1]；kp為顆粒內(nèi)擴散速率常數(shù)（mg·g-1·min-1/2）。

熱力學(xué)方程和熱力學(xué)常數(shù)計算公式為：

Freundlich方程：

Langmuir方程：

Temkin方程：

式中：ce為吸附平衡時液相溶質(zhì)濃度（μg·L-1）；kf和n分別為Freundlich方程的平衡吸附系數(shù)和等溫線的非線性參數(shù)；Qmax為最大吸附量（μg·L-1）；b為與吸附能相關(guān)的吸附常數(shù)（L·μg-1）；K取Freundlich模型的平衡吸附常數(shù)；ΔG0為標準吸附自由能變（kJ·mol-1）；ΔH0為標準吸附焓變（kJ·mol-1）；ΔS0為標準吸附熵變（kJ·mol-1·K-1）；R為摩爾氣體常數(shù)（8.314 J·mol-1·K-1）；T為熱力學(xué)溫度（K）。

1.6 紅外光譜實驗

稱量0.0015 g土壤樣品與0.15 g經(jīng)過預(yù)處理的KBr，在瑪瑙研缽中充分碾磨、壓成透明薄片，進行紅外光譜測試。所用儀器為德國布魯克 Bruker Tensor II傅里葉紅外光譜儀。

2 結(jié)果與討論

2.1 PCB138在草炭土及其組分中吸附動力學(xué)

圖 1為 PCB138在土壤中吸附動力學(xué)曲線。從吸附動力學(xué)曲線可以看出，PCB138在草炭土及其各組分中吸附量隨著時間的增加呈現(xiàn)增大趨勢，吸附量在0—60 min內(nèi)增長較快，PCB138吸附量為 6.85—7.82 μg·g-1，其吸附量為 480 min 吸附量的83%—94%。在60—480 min區(qū)間，隨著時間的延長，PCB138在草炭土中吸附量上升趨勢逐漸變緩，在 240 min后吸附基本達到平衡，此時PCB138 的吸附量值為 8.11—8.57 μg·g-1。實驗表明草炭土及其有機質(zhì)組分對 PCB138的吸附是一個快速吸附過程。這與宋善軍（2010）進行的五氯聯(lián)苯 PCB118吸附動力學(xué)實驗結(jié)果相一致。畢新慧等（2001）也認為PCBs同系物在土壤上吸附數(shù)小時的結(jié)果與15 d基本上沒有顯著差異。從圖1可以看出吸附前期吸附量增長較快，這可能是因為吸附主要是在顆粒表面進行，這種物理性吸附作用和分配作用通常在相當(dāng)短的時間內(nèi)完成（Chen et al.，2012）。隨著反應(yīng)的進行，表面較易吸附的點位已被占據(jù)，具有疏水性的分子將克服顆粒的水分子層到微孔隙，然后通過基質(zhì)擴散進入到顆粒內(nèi)部孔隙而發(fā)生吸附，此階段吸附速率逐漸減慢（Karapanagioti et al.，2000）。

圖1 PCB138在土壤組分中吸附動力學(xué)曲線Fig. 1 Adsorption kinetic curves of PCB138 on soil organic matter fractions

采用準一級動力學(xué)模型、準二級動力學(xué)模型、Elovich動力學(xué)模型和 Weber-Morris模型分別對PCB138在草炭土及其組分中的吸附動力學(xué)過程進行分析，擬合結(jié)果如圖2和表2所示。由圖2和表2可以看出，準一級和準二級動力學(xué)方程擬合結(jié)果要略好于Elovich動力學(xué)模型和Weber-Morris模型，其可決系數(shù)r2在 0.907—0.994之間。圖 2顯示W(wǎng)eber-Morris模型擬合直線不經(jīng)過原點，這表明顆粒內(nèi)擴散并不是唯一的速率控制步驟，吸附過程受到其它階段的共同控制（?zcan et al.，2005）。黃木華（2017）認為 Elovich模型適用于反應(yīng)過程中活化能變化較大的吸附過程。從圖2可知準一級動力學(xué)方程能較好的描述吸附反應(yīng)的初始階段，但隨著吸附反應(yīng)的進行，吸附數(shù)據(jù)逐漸偏離擬合曲線，因此PCB138在草炭土上的吸附用準二級動力學(xué)方程描述更為準確。這可能是由于準一級動力學(xué)方程一般適用于過程比較單一的吸附動力學(xué)；準二級動力學(xué)模型包含吸附的所有過程，如液膜擴散，表面吸附和顆粒內(nèi)部擴散等；在準二級動力學(xué)方程中，PCB138在草炭土原土、胡敏酸及胡敏素中的吸附速率常數(shù)k1值分別為0.014、0.011、0.005，其吸附速率由快到慢依次為原土、胡敏酸和胡敏素。污染物質(zhì)在胡敏素中的吸附速率值最小，這可能是由于胡敏素是慢吸附的主要區(qū)域，短時間內(nèi)有機污染物很難到達并占據(jù)胡敏素上的吸附位（Pan et al.，2006）。

表2 PCB138在土壤組分中吸附動力學(xué)參數(shù)Table 2 Adsorption kinetic parameters of PCB138 on soil organic matter fractions

2.2 PCB138在草炭土及其組分中吸附熱力學(xué)

圖3為不同溫度條件下，不同濃度PCB138在草炭土及其組分中的吸附量。從圖3可知土壤及其有機質(zhì)組分對PCB138的吸附能力略有差別，由強到弱依次為：胡敏酸＞原土壤＞胡敏素。這與袁敏等（2012）的研究結(jié)果相一致，在其研究中發(fā)現(xiàn)胡敏酸對有機污染物質(zhì)的吸附容量更大。圖3也表明，隨著污染物質(zhì)PCB138濃度的增大，土壤及其有機質(zhì)組分對其吸附量也隨之增大。司雄元等（2012）在對幾種土壤及黏土礦物對多氯聯(lián)苯吸附特性的研究中也發(fā)現(xiàn)隨著PCB77和PCB18濃度增加，土壤對PCBs吸附量也隨之增大。

圖2 4種動力學(xué)方程擬合曲線Fig. 2 Fitting curves of four different kinetic equations

隨著體系反應(yīng)溫度的升高，草炭土及其組分對PCB138的吸附量呈現(xiàn)增加趨勢，這表明實驗條件下草炭土及其組分對PCB138的吸附過程為吸熱反應(yīng)。25、35、45 ℃時草炭土原土對PCB138吸附量分別為 1.60—17.4，1.67—17.4、1.70—17.7 μg·g-1，吸附量有所增加，但增幅較小，吸附量約增加0.26%—4.36%。PCB138在胡敏酸和胡敏素中吸附量隨著溫度的升高也有所增加，胡敏酸由25 ℃的1.70—17.7 μg·g-1增長到 45 ℃的 1.73—17.9 μg·g-1。胡敏素由 25 ℃的 1.40—17.1 μg·g-1增長到 45 ℃的 1.56—17.3 μg·g-1。

在 25、35、45 ℃ 3種溫度條件下，PCB138在草炭土及其組分中的吸附等溫線如圖4所示，采用Freundlich、Langmuir和Temkin 3種等溫吸附模型對吸附數(shù)據(jù)的擬合結(jié)果見表3。從表3中可以看出，F(xiàn)reundlich等溫吸附模型擬合結(jié)果的可決系數(shù)在0.92—0.99，擬合效果最好，Langmuir等溫吸附模型的擬合效果最差。Temkin等溫吸附模型的擬合效果居中。這可能是因為 Langmuir吸附模型適用于吸附質(zhì)只是單層的吸附在吸附劑表面，吸附劑表面吸附點位均勻的吸附反應(yīng)，而Freudlich吸附模型認為吸附既存在單層吸附，也存在多層吸附。非線性的 Temkin吸附也是一種不均勻的表面吸附，其內(nèi)部也發(fā)生多層吸附。研究表明土壤對有機物的吸附，不單純是表面吸附，還包括在土壤中的分配過程（展惠英等，2004）。

表3 PCB138在土壤不同組分中吸附參數(shù)及可決系數(shù)Table 3 Adsorption parameters of PCB138 on soil organic fractions and coefficient of determination of equation

圖3 不同溫度條件下不同濃度PCB138在草炭土及其組分中的吸附量Fig. 3 Adsorption capacity of PCB138 having different concentration on soil organic matter fractions at different temperatures

圖4 PCB138在草炭土及其組分中的吸附等溫線Fig. 4 Adsorption isotherms of PCB138 on soil organic fractions at different temperatures

本研究選用 Freundlich等溫吸附模型來描述PCB138在草炭土及其組分中的吸附過程。Freundlich模型中Kf是與吸附容量有關(guān)的常數(shù)。Kf值越大則吸附劑的吸附能力越強。由表3可知，隨著溫度的升高Kf值的大小略有變化。為反映土壤有機質(zhì)結(jié)構(gòu)特征對PCBs吸附的影響，選擇100 μg·L-1的質(zhì)量濃度計算有機碳歸一化分配系數(shù)Koc（Koc=(qe/ce)/foc），其值列在表3中。由表3可知3種土壤組分中，胡敏素的Koc值相對較小，這與林秀梅等（2006）在對菲在土壤有機質(zhì)中的分配行為進行研究時的結(jié)論相一致，在其實驗的 24 h內(nèi)，F(xiàn)reundlich方程中擬合得到胡敏素組分中的lgKoc值要低于胡敏酸的lgKoc值，但在吸附48 h后，胡敏素中菲含量在720 h內(nèi)在連續(xù)不斷地增加。這也表明胡敏素是慢吸附過程的主要區(qū)域。n值反映吸附的非線性程度，也可以反應(yīng)吸附點位能量分布，從表4可以看出草炭土及其組分對PCB138的等溫吸附曲線有一定程度的非線性。Freundlich方程擬合得到的n值均小于且偏離1，表明吸附等溫線的非線性及吸附點位的異質(zhì)性，有學(xué)者認為這可能與土壤的非均質(zhì)性有關(guān)（Zhang et al.，2012）。

表4 PCB138在草炭土中吸附的熱力學(xué)參數(shù)數(shù)值Table 4 Thermodynamic parameters of PCB138 on peat soil

等溫吸附數(shù)據(jù)顯示，PCB138在炭土原土及胡敏酸中吸附能力略有波動，但基本表現(xiàn)為在胡敏酸中的吸附能力要大于在原土中的吸附能力。這可能是因為胡敏酸是有機質(zhì)中含量較高，較穩(wěn)定的組分，含有多種活性官能團，具有較高的吸附能力，可與許多有機無機物質(zhì)發(fā)生相互作用（蔣煜峰等，2005）。此外土壤的胡敏酸等可溶性有機質(zhì)還可以通過對疏水性污染物的增溶作用來提高吸附量（劉靜等，2010）。

2.3 PCB138在草炭土中吸附熱力學(xué)參數(shù)

一些學(xué)者認為表征熱力學(xué)過程的吸附焓變等參數(shù)可以作為研究污染物吸附機理的重要手段（張勁強，2008）。應(yīng)用吉布斯方程計算在 25—45 ℃范圍內(nèi)草炭土對PCB138吸附的熱力學(xué)參數(shù)數(shù)值，結(jié)果見表4。表4中數(shù)據(jù)顯示，實驗所得的熱力學(xué)參數(shù) ΔG0、ΔH0和 ΔS0值均＞0。這表明 PCB138在草炭土上的吸附為非自發(fā)進行的吸熱反應(yīng)，吸附過程中固液界面的無序性和體系混亂度增加。

ΔG0值可以在一定程度上反映吸附過程推動力的大小，其絕對值越大，表明吸附的推動力越大（張曉蕾等，2012）。由表4可知，隨著溫度的升高，ΔG0值逐漸減小，這可能是由于吸附的驅(qū)動力占用了高能吸附位點而有所降低（陳守益等，2018）。Tan et al.（2015）認為這種數(shù)值的波動表明在低濃度的吸附過程中伴隨有熵的補償。ΔH0的數(shù)值可以用來說明吸附機制，當(dāng)ΔH0在2—40 kJ·mol-1時，氫鍵可能在吸附過程中起主導(dǎo)作用（張濤等，2014）。在一些研究中氫鍵主要被用于描述PAHs，芳香族碳氫化合物和疏水區(qū)域較大的農(nóng)藥被強吸附到腐殖質(zhì)土壤的原因（謝國紅等，2006；Spark et al.，2002）。本研究計算得到的ΔH0大于0，表明吸附反應(yīng)是吸熱反應(yīng)，這與吸附量隨著溫度升高而增加的實驗結(jié)果相一致。

2.4 紅外光譜分析

草炭土及實驗所提取的草炭土有機質(zhì)組分的紅外光譜如圖5所示。從圖5可以看出草炭土及其組分在官能團區(qū)3428、3125、1400 cm-1等頻率處，指紋區(qū)中的1040 cm-1和886 cm-1等頻率處有明顯吸收。一般認為，3500—3100 cm-1的寬峰為O-H伸縮振動譜帶。1400 cm-1附近有明顯的脂肪族CH3和CH2中C-H的振動吸收峰。1040 cm-1處吸收峰多為碳水化合物或多糖結(jié)構(gòu)中的 C-O伸縮振動及無機物的Si-O伸縮振動吸收峰。900—650 cm-1的是芳環(huán)上 C-H的面外變形振動位置（包桂奇等，2013）。

圖5 草炭土及其組分的紅外光譜Fig. 5 FTIR spectra of the peat soil and its organic matter fractions

由圖4可以看出草炭土原土與胡敏酸的吸收峰比較相似，這表明其官能團的組成情況也較為相似，實驗提取的胡敏素在譜圖上的吸收峰較少，且吸收峰面積也較小，這可能是由于胡敏素中芳香性組分含量較低而脂肪族組分含量較高所致（Gunasekara et al.，2003）。

3 結(jié)論

（1）六氯聯(lián)苯PCB138在草炭土原土、胡敏酸和胡敏素中吸附是一個快速過程，240 min可以達到吸附平衡。吸附速率表現(xiàn)為在原土中最快，k1值為0.014，在胡敏素中最慢，k1值為0.005。準一級、準二級、Elovich和 Weber-Morris模型分析結(jié)果表明準二級動力學(xué)方程能更好的描述其吸附過程。

（2）在實驗條件下，增大PCB138濃度，提高反應(yīng)體系溫度都能增大草炭土原土壤、胡敏酸和胡敏素對 PCB138吸附量。土壤及其組分對PCB138的吸附能力表現(xiàn)為胡敏酸的吸附能力最強。與Langmuir和Temkin等溫吸附模型相比，F(xiàn)reundlich模型能更好的描述PCB138的吸附過程。

（3）等溫吸附實驗計算得到的ΔG0、ΔH0和ΔS0值均＞0，這表明PCB138在草炭土上的吸附為非自發(fā)進行的吸熱反應(yīng)，吸附過程中體系混亂度增加。