南京大学学报(自然科学), 2024, 60(1): 130-140 doi: 10.13232/j.cnki.jnju.2024.01.013

场地有机污染物吸附行为多参数线性自由能模型研究

刘昆1,2, 南晨曦1, 孔令冉1, 刘慧婷3, 付翯云1, 瞿晓磊,1

1.污染控制与资源化研究国家重点实验室,南京大学环境学院,南京,210023

2.江苏省环境工程技术有限公司,南京,210023

3.辽宁省废水治理技术重点实验室,沈阳理工大学环境与化学工程学院,沈阳,110159

Poly⁃parameter linear free energy relationship models for organic pollutant sorption prediction at contaminated sites

Liu Kun1,2, Nan Chenxi1, Kong Lingran1, Liu Huiting3, Fu Heyun1, Qu Xiaolei,1

1.State Key Laboratory of Pollution Control and Resource Reuse, School of Environment, Nanjing University, Nanjing, 210023, China

2.Jiangsu Environmental Engineering Technology Co. , Ltd. , Nanjing, 210023, China

3.Key Laboratory for Technology of Wastewater Treatment in Liaoning Province, School of Environmental and Chemical Engineering, Shenyang Ligong University, Shenyang, 110159, China

通讯作者: E⁃mail:xiaoleiqu@nju.edu.cn

收稿日期: 2023-09-18  

基金资助: 国家重点研发计划.  2019YFC1804201.  2020YFC1807002
污染控制与资源化研究国家重点实验室开放基金.  PCRRF22034
国家自然科学基金青年基金.  22206132

Received: 2023-09-18  

摘要

准确评估污染物在污染场地土壤中的吸附分配过程是污染场地安全管理和再利用开发的重要环节.以我国典型电子垃圾拆解场地为研究对象,解析了场地土壤中有机质的吸附特征,建立了适用于场地有机污染物分配评估的有机碳标化分配系数(KOC)预测模型,揭示了有机污染物在场地土壤中有机质上的吸附机制.该场地土壤有机质对有机污染物具有较强的吸附能力.构建了预测场地土壤中有机质吸附有机污染物的KOC的多参数线性自由能模型(pp⁃LFER).pp⁃LFER模型相较于常用的单参数线性自由能模型(sp⁃LFER)展现出更好的场地吸附数据的预测能力(R2=0.919).同时,采用其他文献报道的pp⁃LFER模型对该场地污染物KOC进行预测,发现预测偏差较高(RMSE>1.12),这表明pp⁃LFER模型的预测效果受场地土壤中有机质的性质影响较大,其跨区域应用性仍待提升.进一步结合pp⁃LFER模型中的分子结构描述符和各项系数解析了土壤中有机质的吸附机制,发现疏水作用和极化作用是吸附过程的重要作用力,氢键作用显著影响极性化合物的吸附过程.本研究基于实际污染场地土壤构建了高精度吸附预测模型,为场地污染风险评估和修复利用提供了更准确的技术手段.

关键词: 场地污染 ; 土壤有机质 ; 有机污染物 ; 有机碳标化分配系数 ; 预测模型

Abstract

Accurate assessment of the sorption process of pollutants is an important step for the safe management and development of contaminated sites. We used a representative e⁃waste dismantling site in China as a research object to analyse the sorption properties of soil organic matter,developed a prediction model for the partition coefficient (KOC) of organic carbon,and revealed the sorption mechanism of organic pollutants on the soil organic matter in the site. The results show that the site soil organic matter has a strong sorption capacity for organic pollutants. Poly parameter linear free energy relationships(pp⁃LFER) show a better fit to the site sorption data thansingle parameter linear free energy relationships (sp⁃LFER) (R2=0.919). However,the validation of the established pp⁃LFER for soil organic matter from different areas shows a significant deviation (RMSE>1.12). The performance of pp⁃LFERs is affected by the soil organic matter propreties. It may not be able to accurately predict the sorption behavior of soils in different sites. The sorption mechanism was further analysed through molecular structure descriptors and the physicochemical significance of the coefficients. It was found that hydrophobic effect and polarisation are important forces in the sorption process,with hydrogen bonding having a significant effect on strongly polar compounds. This study provides the high⁃precision sorption prediction model for the actual contaminated site. It provides the managers with reliable model options for risk assessment.

Keywords: site contamination ; soil organic matter ; organic pollutants ; organic carbon normalized partition coefficient ; prediction model

PDF (986KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘昆, 南晨曦, 孔令冉, 刘慧婷, 付翯云, 瞿晓磊. 场地有机污染物吸附行为多参数线性自由能模型研究. 南京大学学报(自然科学)[J], 2024, 60(1): 130-140 doi:10.13232/j.cnki.jnju.2024.01.013

Liu Kun, Nan Chenxi, Kong Lingran, Liu Huiting, Fu Heyun, Qu Xiaolei. Poly⁃parameter linear free energy relationship models for organic pollutant sorption prediction at contaminated sites. Journal of nanjing University[J], 2024, 60(1): 130-140 doi:10.13232/j.cnki.jnju.2024.01.013

随着工业化进程的持续推进,电子、化工、石油等行业的不断发展,场地土壤有机污染问题日益突出1.近年来,我国相继出台并实施了《土壤污染防治行动计划》和《中华人民共和国土壤污染防治法》,场地污染的调查、监管和修复工作受到了高度关注.场地有机污染物的主要赋存相为土壤和地下水,土壤有机质(Soil Organic Matter,SOM)对有机污染物的吸附作用较强2,显著影响有机污染物的迁移归趋.场地污染具有隐蔽性强和迁移规律复杂等特点3,准确评估场地有机污染的吸附分配行为是实现污染场地土壤安全利用开发的重要环节之一.有机碳标化分配系数(Organic Carbon Normalized Partition Coefficient,KOC)是有机污染物在土壤环境中迁移、归趋和风险评估的重要参数之一4-5,可通过吸附分配实验测定.但由于场地土壤有机质性质各异、有机污染物种类众多,不同污染物的KOC差异显著,因此,KOC的实验测定方法难以满足实际应用需求,开发准确、可靠的KOC预测模型对场地有机污染风险评估具有重要意义.

目前,常用于KOC预测的模型是线性自由能模型(Linear Free Energy Relationships,LFERs),包括单参数线性自由能模型(Single Parameter LFER,sp⁃LFER)和多参数线性自由能模型(Poly Parameter LFER,pp⁃LFER)46-10.线性自由能模型通过将有机化合物结构性质作为参数,实现KOC的预测.sp⁃LFER以正辛醇⁃水分配系数(Octanol⁃Water Partition Coefficient,KOW)或溶解度(Solubility,S)作为量化有机污染物性质的参数411,由于仅需一个有机化合物常规参数便可达到KOC预测的目的,sp⁃LFER具有很强的实用性.然而,sp⁃LFER主要适用于以疏水分配为主要作用机制的吸附行为预测,很难准确地预测吸附机制较为复杂的土壤吸附212.pp⁃LFER是在Abraham模型基础上发展起来的13-14,利用分子结构描述符(也称Abraham溶质参数描述符)量化有机污染物与吸附相之间的相互作用.通常认为五个分子结构描述符就可以全面地描述有机污染物与吸附相之间的相互作用,因此,pp⁃LFER为预测土壤有机质的吸附行为提供了模型选择71015-16,但不同区域性质差异较大的土壤有机质的pp⁃LFER公式系数差异显著17.因而,构建有针对性的场地土壤有机质pp⁃LFER,是实现污染场地KOC精准预测的保障.

电子垃圾拆解场地是我国受历史影响形成的一种典型污染场地,该类场地中的电子垃圾通常包含各类金属、含阻燃剂塑料制品、电气设备产品外壳等18,若处置不当,极易造成严重的有机污染,给人们带来显著的健康风险.本研究的主要目的是聚焦我国电子垃圾拆解场地,构建基于pp⁃LFER的准确的KOC预测模型,提升场地有机污染风险的评估精度.研究围绕浙江省台州地区某典型电子垃圾拆解场地开展,利用批吸附实验法获取场地土壤有机质KOC,结合Abraham溶质参数描述符,通过多元线性回归法建立pp⁃LFER;进一步对pp⁃LFER的分子结构参数进行贡献率分析,解析有机污染物在土壤有机质上的吸附机制,为场地污染风险评估和修复与利用提供更准确的技术手段.

1 材料与方法

1.1 样品采集

本研究随机采集了浙江省台州市某电子垃圾拆解区(28°32'0.96'' N,121°21'53.15'' E)的七个土壤样品,深度为0~20 cm的表层土壤,该地区是我国典型的电子废弃物拆解地19.表层土壤样品采集后,经自然风干并挑出植物残体和砖粒石子等异物后,研碎过60目筛后即获得用于实验的土壤样品,将其放置于广口瓶中备用.

1.2 实验材料

2,2′,4,4′⁃四溴联苯醚(BDE⁃47,100.0%)购自美国AccuStandard公司;1,3,5⁃三氯苯(99%)和芘(98%)购自美国Sigma⁃Aldrich公司;蒽(≥99%)购自瑞典Fluka公司;阿特拉津(97%)购自上海麦克林生化科技有限公司;对硝基苯酚(99%)购于Acros organics公司;氯化钠(分析纯)购自南京化学试剂有限公司;叠氮化钠(高级纯)购自上海睿铂赛生物科技公司;正己烷(色谱纯)和甲醇(色谱纯)购自美国Tedia公司.实验中所用的去离子水由ELGA Labwater系统(PURELAB Ultra,ELGA LabWater Global Operations,英国)制备(电阻率为18.2 MΩ·cm,25 °C).

1.3 批吸附实验

取15 mg土壤于配有聚四氟乙烯瓶盖的22 mL EPA瓶中,加入适量的背景液预润湿土壤样品,本实验背景液为含0.02 mol·L-1氯化钠和200 mg·L-1叠氮化钠的混合溶液,其中叠氮化钠是微生物生长抑制剂.预润湿后向EPA瓶中加入不同浓度梯度的有机污染物储备液,加入阿特拉津初始溶液的浓度范围为3.1×10-4~2.8×10-3 mmol·L-1,加入蒽初始溶液的浓度范围为1.7×10-4~1.2×10-3 mmol·L-1,加入芘初始溶液的浓度范围为1.8×10-4~1.4×10-3 mmol·L-1,加入对硝基苯酚初始溶液的浓度范围为5.5×10-4~4.2×10-3 mmol·L-1,加入1,3,5⁃三氯苯初始溶液的浓度范围为4.4×10-4~3.3×10-3 mmol·L-1,加入BDE⁃47初始溶液的浓度范围为3.9×10-6~2.9×10-5 mmol·L-1.加入甲醇的总体积小于溶液总体积的0.1%,防止其发生溶剂效应.用背景液充满EPA瓶,将EPA瓶放置在旋转混合器(转速为6 r·min-1)中,室温条件下避光翻转震荡3 d,以达到吸附平衡.所有样品在达到平衡状态后取下,静置约6 h,以3800 r·min-1转速离心15 min,用玻璃注射器抽取上清液.阿特拉津、蒽、芘和对硝基苯酚采用配备光电二极管阵列检测器的高效液相色谱仪测定,色谱柱选择4.6×150 mm Zorbax Eclipse XDB⁃C18,1,3,5⁃三氯苯和BDE⁃47采用配备ECD检测器和DB⁃WAX色谱柱的气相色谱测定.上清液中有机污染物的具体测定方法见表1.有机污染物的浓度根据标准曲线计算,标准曲线的绘制方法与前述做法一致,只是标准溶液中不含土壤样品.标准曲线至少包含六个浓度点,覆盖整个实验浓度范围,线性相关性R20.99.在所有吸附实验中每个样品点均设置两组平行样.

表1   有机化合物的检测方法

Table 1  Analytical methods for the organic compounds

有机化合物色谱方法流动相(vv检测波长(nm)
阿特拉津液相色谱水∶甲醇(40∶60)222
液相色谱水∶甲醇(10∶90)254
液相色谱水∶甲醇(10∶90)270
对硝基苯酚液相色谱水∶乙腈(40∶60)220
有机化合物色谱方法升温条件
1,3,5⁃三氯苯气相色谱进样口和检测器温度分别设为270 °C和280 °C,炉温起始温度设为100 °C,以10 °C·min-1的速度升温至200 °C,再以30 °C·min-1的速度升温至250 °C,保持3 min.
BDE⁃47气相色谱

进样口和检测器温度分别设为290 °C和310 °C,柱温初始温度设为80 °C,保持1 min,再以

20 °C·min-1的速度升温至280 °C,保持5 min,最后以2 °C·min-1的速度升温至300 °C,保持2 min.

新窗口打开| 下载CSV


1.4 吸附模型构建

lgKOC的计算方法见式(1)和式(2)20

Kd=qe/Ce
lgKOC=lgKd/fOC

其中,qe是有机化合物在土壤有机质上的平衡浓度,Ce是有机化合物在水相的平衡浓度,fOC是土壤有机碳含量.

本文采用的sp⁃LFER表达式如下:

lgKOC=b+algKOW

其中,KOW为有机化合物疏水性定量参数,a为线性回归系数,b为常数.

本研究采用的pp⁃LFER见式(4),常用于两个凝聚相之间的分配系数建模9表达:

lgKOC=c+eE+sS+aA+bB+vV

其中,E为化合物的过量摩尔折射率,S为化合物的极化率/偶极性参数,A为化合物的氢键酸度,B为化合物的氢键碱度,V为化合物的McGowan摩尔体积;esabv为多元线性回归系数,表示有机化合物与有机和水两相间的相互作用差异,c为常数.

有机化合物的分子结构描述符获取:数据参考UFZ⁃LSER数据集网站21,有机化合物的lgKOW和分子结构描述符见表2.

表2   有机化合物的lgKOW和分子结构描述符

Table 2  lgKOW and molecular structure descriptors of organic compounds

lgKOWESABV
阿特拉津2.651.221.290.171.011.62
4.682.291.3400.281.45
5.132.811.7100.281.58
对硝基苯酚1.911.071.720.820.260.95
1,3,5⁃三氯苯4.190.980.73001.08
BDE⁃476.802.321.4500.342.08

新窗口打开| 下载CSV


模型拟合方法:利用SPSS统计软件对数据进行多元线性回归分析,以lgKOC为因变量,以分子结构描述符作为自变量建立pp⁃LFER,以KOW为自变量建立sp⁃LFER.

2 结果与讨论

2.1 场地土壤吸附行为

本研究涉及的有机污染物包括非极性和极性有机化合物,蒽、芘、1,3,5⁃三氯苯和BDE⁃47为非极性有机化合物,阿特拉津和对硝基苯酚为极性有机化合物.表2给出了六种有机化合物的分子结构描述符数值,分子结构描述符范围分别是E:0.98~2.81,S:0.73~1.72,A:0~0.82,B:0~1.01,V:0.95~2.08.图1总结了七个场地土壤样品吸附有机化合物的吸附等温线.表3总结了土壤有机碳含量和有机碳标化分配系数.该场地土壤对阿特拉津、蒽、芘、对硝基苯酚、1,3,5⁃三氯苯和BDE⁃47的吸附KOC平均值分别为11228,286686,310849,61528,39914和796206 L·kg-1·C-1.本研究涉及的非极性有机化合物吸附量基本呈现出随有机化合物疏水性(lgKOW)的增加而增大的趋势.土壤有机质对多环芳烃等非极性有机化合物的吸附主要受到疏水作用控制1222-23.值得注意的是,本研究场地的土壤有机质对三氯苯的吸附能力为105 L·kg-1·C-1数量级,根据报道三氯苯在辽中南地区、湖南和广东地区土壤有机质上的吸附能力均为104 L·kg-1·C-1数量级24-25.本场地土壤有机质对多环芳烃的吸附能力达到了106 L·kg-1·C-1数量级,接近煤炭源有机质所能达到的吸附能力23.国内典型土壤类型—草甸棕壤对多环芳烃的吸附能力大约在105 L·kg-1·C-1数量级26.土壤有机质对非极性有机污染物的亲和性与有机质结构聚合度密切相关,其中聚合度可以用极性指数(O+N)/C表征,有机质结构聚合度越大,硬碳类有机质的含量越高,吸附能力增强27.草甸棕壤多出现在中国东北地区和内蒙古地区,与辽中南地区和广东地区土壤有机质极性指数相似且均远小于褐煤有机质24-25,说明草甸棕壤、辽中南地区和广东地区土壤有机质极性更强,疏水性较弱,与非极性有机污染物亲和性较弱.本研究所用土壤有机质对非极性有机化合物很强的吸附能力,可能与其结构聚合度大、硬质碳含量高有关228-29.研究发现场地土壤有机质对极性有机污染物表现出相似的强吸附现象,对阿特拉津的吸附能力高于国内东北地区黑土、冻融黑土和北京地区砂质壤土30-32,对对硝基苯酚的吸附能力接近表面活性剂改性的膨润土33.土壤有机质对极性有机化合物的吸附过程由极性(如氢键)相互作用主控,羧基和羟基官能团是形成氢键的关键结构34.该场地土壤有机质对极性化合物表现出的强吸附能力说明其具有丰富的极性官能团.

图1

图1   土壤有机质对非极性有机污染物(a,右上角图示BDE⁃47吸附等温线)和极性有机污染物(b)的吸附等温线

Fig.1   Soil organic matter sorption isotherms for nonpolar organic pollutants (a) with an insert figure of BDE⁃47 and polar organic pollutants (b)


表3   土壤有机碳含量及土壤有机质吸附有机污染物的有机碳标化分配系数(KOC,单位:L·kg-1·C-1

Table 3  Soil organic carbon content and organic carbon normalized partition coefficients (KOC,unit: L·kg-1·C-1) for organic compounds sorption to soil organic matter

土壤编号有机碳含量KOC(阿特拉津)KOC(蒽)KOC(芘)KOC(对硝基苯酚)KOC (1,3,5⁃三氯苯)KOC (BDE⁃47)
训练集
1号土壤2.45%53831216191764982966975509578096
2号土壤1.68%57151396371606134211071285566239
6号土壤0.75%1276445289845924010649970469550808
测试集
3号土壤1.96%601221428921231638561106171541700
4号土壤1.50%90991807171860435142910139616595
5号土壤1.30%86503006083194005915616943860994
7号土壤0.78%3097459703566183310327624434859014

新窗口打开| 下载CSV


2.2 预测模型建立与验证

基于该场地七种表层土壤对六种有机污染物构建的吸附数据集,随机选择三种土壤有机质对六种有机污染物的吸附数据构建训练集,结合目前常用于KOC预测的线性自由能模型,利用多元线性回归方法,以lgKOC为因变量,pp⁃LFER以分子结构描述符为自变量,sp⁃LFER以lgKOW为自变量,构建用于场地土壤有机质吸附的预测模型.pp⁃LFER和sp⁃LFER拟合结果如式(5)和式(6)所示:

lgKOC=3.863(±0.436)+1.096(±0.677)E-1.716(±1.419)S+2.341(±1.578)A-1.270(±0.424)B+1.088(±0.259)V
n=18,R2=0.919,adjusted R2=0.885
lgKOC=0.316±0.058lgKOW+3.640(±0.264)
n=18,R2=0.647,adjusted R2=0.625

拟合结果显示,pp⁃LFER对训练集数据拟合R2=0.919,sp⁃LFER的R2=0.647.利用训练集构建模型一般认为R2>0.6是模型对数据具有较好拟合能力的表现35-36.尽管pp⁃LFER和sp⁃LFER对训练集数据拟合均达到了R2>0.6,但pp⁃LFER显著优于sp⁃LFER.这主要由于sp⁃LFER本身构成简单,仅具备描述有机化合物疏水性的指标,其主要考虑吸附过程中的疏水分配作用.但训练集包含非极性和极性有机化合物.极性有机化合物的吸附受极性(氢键)相互作用显著影响,无法在sp⁃LFER模型中体现出来.因此,吸附机制适用更为广泛的pp⁃LFER相较于sp⁃LFER展现了更好的拟合能力,对训练集的描述性更为精准223.

进一步利用剩余的四种土壤对六种有机化合物吸附数据构建验证集,检验本研究得到的两种吸附预测模型的预测精度.通常对于lgKOC的预测认为均方根误差(Root Mean Squared Error,RMSE)小于1便可接受,若小于0.5则认为具有较好的预测性10.图2展示了pp⁃LFER和sp⁃LFER的预测精度,虽然两个模型的预测精度均在可接受范围内,即RMSE小于1,但pp⁃LFER更为准确.这说明在预测有机污染物种类多样、机制较为复杂的吸附行为时,pp⁃LFER更加准确.

图2

图2   验证集中四种土壤有机质对六种有机化合物吸附实验值与pp⁃LFER(a,式(5)))和sp⁃LFER(b,式(6))预测值比较,虚线表示1∶1的线

Fig.2   Comparison between experimental lgKOC and predictions made by sp⁃LFER (a,eqs (5)) and pp⁃LFER models (b,eqs (6)) for organic pollutants sorption to soil organic matter in the validation set,the dashed line represents the 1∶1 line


2.3 模型比较与通用性分析

基于本研究场地构建的验证集,选择文献中已建立的土壤有机质sp⁃LFER(式(7)26)和pp⁃LFER(式(8)7,(9)10,(10)20),比较不同场地间预测模型的准确性与通用性.

lgKOC=1.10 lgKOW-0.68
lgKOC=0.14+1.10E-0.72S+0.15A-1.98B+2.28V
lgKOC=0.25+0.81E-0.16S+0.20A-2.34B+2.11V
lgKOC=-1.04+0.31E+1.27S-0.10A-3.94B+3.71V

根据图3中的比较结果不难发现,文献中针对其他区域土壤有机质建立的sp⁃LFER和pp⁃LFER的预测准确性较差,RMSE均大于1.式(7)中的sp⁃LFER是基于七种多环芳烃有机污染物构建的,吸附机制相对简单,但是本研究所选择的有机化合物性质具有多样性,吸附机制更为复杂,如前所述并不符合sp⁃LFER适用的机制条件,造成了sp⁃LFER预测的严重偏差.另一方面,由于线性自由能模型仅考虑了有机化合物性质,并未很好地考虑有机质性质,使得模型通用性较差37,一旦有机质性质具有较大区别,便会出现较大的预测偏差9.本场地土壤有机质与其他地区土壤性质差异显著,因此其他地区建立得到的pp⁃LFER并未在本场地展现良好的预测性.pp⁃LFER的这种地区不通用性或模型依赖标化的现象在文献中也有所报道.之前有研究围绕波罗的海的波的尼亚湾和另一指定区域针对性地构建了两个高精度pp⁃LFER.然而由于波罗的海水体溶解性有机质的异质性较强,pp⁃LFER并未在波罗的海其他区域的预测中展现很好的预测效果37.Bronner and Goss10通过柱实验和批吸附实验收集了泥炭土有机质对农药等有机化合物的吸附数据,尽管利用该吸附数据得到了高精度pp⁃LFER,但利用此次实验得到的新数据验证Niederer et al38建立的pp⁃LFER时,模型偏差均明显较大,甚至其中的两个模型对KOC预测偏差超过一个数量级.总之,通过以上的模型验证与比较发现,pp⁃LFER是一种可以应对吸附机制复杂的预测模型,相较于sp⁃LFER更适合用于描述污染场地有机化合物的迁移规律.但是受限于模型未很好地考虑有机质性质,其并不具备区域通用性.为实现不同区域或污染场地土壤有机质吸附性能的准确预测,pp⁃LFER的标化工作有待开展.

图3

图3   土壤有机质吸附有机化合物的lgKOC实测值与sp⁃LFER(a)和pp⁃LFER(b~d)模型的lgKOC预测值比较

Fig.3   Comparison between experimental lgKOC and predictions made by sp⁃LFER (a) and pp⁃LFER (b~d) models for organic pollutions sorption to soil organic matter


2.4 吸附机理解析与分子间作用力贡献率比较

鉴于pp⁃LFER分子结构描述符和各项系数具有特定的物理化学意义,对其进行分析和比较可以有效地解析土壤有机质吸附过程的主要机制.本文总结了文献中土壤/底泥/泥炭有机质pp⁃LFER各项系数(表4).通过比较不难发现,本研究基于土壤有机质拟合得到的pp⁃LFER系数的大小不同于其他区域,说明土壤有机质与有机化合物间的相互作用强度具有差异性,针对本研究区域土壤有机质构建特异性pp⁃LFER很有必要.但本研究得到的系数与文献中以土壤/底泥有机质为吸附相的正负趋势一致,eav为正值,sb为负值,系数为正表明有机化合物与有机相之间的作用力更强,系数为负则表明有机化合物与水相之间的作用力更强.进一步分析本研究得到的pp⁃LFER,以解析土壤有机质吸附有机化合物的主要作用力.e为正值说明有机化合物更易与土壤有机质发生色散和诱导力的作用,表现出溶质分子更易于从水相扩散到土壤有机相的趋势.v为正值说明溶质分子趋向有机相的过程更易通过空穴作用分配其中,这是由于水分子具有排列高度有序的结构特点,凝聚性较强1040.溶质分子若要打破水分子的氢键作用形成空穴,需要消耗较大的自由能39,且溶质分子越大,需要消耗的能量越大,因此更倾向于通过空穴作用分配到有机相.S代表化合物的极化/偶极作用,地质吸附相的偶极性/极化率往往比水小41,有机化合物更易与水相发生极化/偶极作用,因此s一般为负.本研究s绝对值相对文献值更大,表明了本研究区域土壤有机质与水更大的极化率/偶极性差异,与其他区域土壤有机质存在较强异质性.AB分别表示化合物的氢键酸度和碱度,是pp⁃LFER中的特异性作用参数.这一对相互作用是土壤有机质与极性化合物相互作用的专属作用力,通常与土壤有机质极性官能团紧密相关42.本研究a为正值说明土壤有机质接受电子的能力大于水相,提供更强的氢键碱作用43a在系数中大小相对更突出,说明固液两个分配相的氢键碱作用性质相差较大.然而土壤有机质的氢键碱作用仅对高酸度化合物有意义7,虽然阿特拉津和对硝基苯酚均被定义为极性有机化合物,但阿特拉津的氢键酸度(0.17)远小于对硝基苯酚(0.82),因此a值相对突出主要反映了土壤有机质与对硝基苯酚之间的相互作用.b为负值表明水相可提供更强的氢键酸作用,具有更强的提供电子的能力.通过进一步对不同种类有机化合物各项作用力贡献率的比较(图4),证明了氢键作用仅对极性有机化合物的KOC预测有明显贡献,极化作用、空穴作用和色散作用是实现土壤有机质KOC准确预测的重要作用力.

表4   不同pp⁃LFER分子结构描述符系数汇总

Table 4  Summary of the coefficients of different pp⁃LFER molecular structure descriptors

esabv参考文献
土壤/底泥有机质0.740-0.31-2.272.09Poople et al[39]
土壤/底泥有机质1.10-0.720.15-1.982.28Nguyen et al[7]
土壤/底泥有机质1.08-0.830.28-1.852.55Nguyen et al[7]
土壤/底泥有机质0.89-0.580-1.991.29Kipka et al[15]
土壤/底泥有机质1.20-0.08-0.19-1.811.16Kipka et al[15]
泥炭0.311.27-0.10-3.943.71Endo et al[20]
泥炭0.430.190.02-3.833.51Endo et al[20]
泥炭0.81-0.61-0.21-3.442.99Bronner et al[10]
土壤有机质1.10-1.722.34-1.271.09本研究

新窗口打开| 下载CSV


图4

图4   不同种类污染物的各项作用力贡献率比较

Fig.4   Comparison of the contribution of each force for different types of pollutants


3 结论

(1)本研究场地土壤有机质吸附有机污染物的能力强,与其他地区土壤有机质存在显著区别.与辽中南、湖南和广东地区土壤有机质相比,对非极性有机化合物的吸附能力强一个数量级,与草甸棕壤和东北黑土相比,对极性有机化合物的吸附能力强一个数量级.

(2)基于分子结构描述符建立的pp⁃LFER可以准确描述机制复杂的吸附过程,模型对训练集数据的拟合能力远高于常用的sp⁃LFER.然而,pp⁃LFER在使用时需要标化过程,地区间公式通用性不强.结合实际场地土壤有机质性质,构建针对各地区特点的pp⁃LFER,是实现污染场地吸附准确预测的重要环节.

(3)通过比较和分析分子结构描述符各项系数,证明极化作用、空穴作用和色散作用是KOC准确预测的关键因子,氢键作用仅对强极性有机化合物的吸附过程产生显著影响.

(4)该模型基于我国典型电子垃圾拆解污染场地建立,KOC的准确预测有助于更好地理解污染场地有机化合物的迁移规律和归趋,为我国有机污染物场地的管理提供了有效的模型.

参考文献

Zhao C, Dong Y, Feng Y P,et al.

Thermal desorption for remediation of contaminated soil:A review

Chemosphere,2019,221841-855.

[本文引用: 1]

Cornelissen G, Gustafsson O, Bucheli T D,et al.

Extensive sorption of organic compounds to black carbon,coal,and kerogen in sediments and soils:Mechanisms and consequences for distribution,bioaccumulation,and biodegradation

Environmental Science & Technology,2005,39(18):6881-6895.

[本文引用: 4]

陈晶中,陈杰,谢学俭.

土壤污染及其环境效应

土壤,2003,35(4):298-303.

[本文引用: 1]

Chen J Z, Chen J, Xie X J,et al.

Soil pollution and its environmental impact

Soil,2003,35(4):298-303.

[本文引用: 1]

Chiou C T, Peters L J, Freed V H.

A physical concept of soil⁃water equilibria for nonionic organic compounds

Science,1979,206(4420):831-832.

[本文引用: 3]

Schwarzenbach R P, Gschwend P M, Imboden D M. Environmental organic chemistry. 3rd ed. New JerseyWiley,20171024.

[本文引用: 1]

Goss K U, Schwarzenbach R P.

Linear free energy relationships used to evaluate equilibrium partitioning of organic compounds

Environmental Science & Technology,2001,35(1):1-9.

[本文引用: 1]

Nguyen T H, Goss K U, Ball W P.

Polyparameter linear free energy relationships for estimating the equilibrium partition of organic compounds between water and the natural organic matter in soils and sediments

Environmental Science & Technology,2005,39(4):913-924.

[本文引用: 5]

Zhu D Q, Pignatello J J.

A concentration⁃dependent multi⁃term linear free energy relationship for sorption of organic compounds to soils based on the hexadecane dilute⁃solution reference state

Environmental Science & Technology,2005,39(22):8817-8828.

Endo S, Goss K U.

Applications of polyparameter linear free energy relationships in environmental chemistry

Environmental Science & Technology,2014,48(21):12477-12491.

[本文引用: 2]

Bronner G, Goss K U.

Predicting sorption of pesticides and other multifunctional organic chemicals to soil organic carbon

Environmental Science & Technology,2011,45(4):1313-1319.

[本文引用: 7]

Chiou C T, Porter P E, Schmedding D W.

Partition equilibriums of nonionic organic compounds between soil organic matter and water

Environmental Science & Technology,1983,17(4):227-231.

[本文引用: 1]

Xing B S, Pignatello J J.

Dual⁃mode sorption of low⁃polarity compounds in glassy poly(vinyl chloride) and soil organic matter

Environmental Science & Technology,1997,31(3):792-799.

[本文引用: 2]

Abraham M H.

Hydrogen bonding. 31. Construction of a scale of solute effective or summation hydrogen⁃bond basicity

Journal of Physical Organic Chemistry,1993,6(12):660-684.

[本文引用: 1]

Abraham M H, Poole C F, Poole S K.

Classification of stationary phases and other materials by gas chromatography

Journal of Chromatography A,1999,842(1-2):79-114.

[本文引用: 1]

Kipka U, Di Toro D M.

A linear solvation energy relationship model of organic chemical partitioning to dissolved organic carbon

Environmental Toxicology and Chemistry,2011,30(9):2023-2029.

[本文引用: 3]

Vitale C M, Di Guardo A.

Predicting dissolved organic carbon partition and distribution coefficients of neutral and ionizable organic chemicals

Science of the Total Environment,2019,6581056-1063.

[本文引用: 1]

刘昆,付翯云,瞿晓磊.

天然有机质疏水性定量方法及吸附预测模型

环境化学,2023,42(1):71-79.

[本文引用: 1]

Liu K, Fu H Y, Qu X L.

The hydrophobicity quantification methods and sorption models for natural organic matter

Environmental Chemistry,2023,42(1):71-79.

[本文引用: 1]

矫旭东,陆泗进,王业耀.

电子垃圾拆解场地污染特征及土壤监测评估的研究概况

土壤通报,2014,45(5):1265-1272.

[本文引用: 1]

Jiao X D, Lu S J, Wang Y Y,et al.

Pollution characteristics of an e⁃waste recycling area and its soil monitoring and assessment

Chinese Journal of Soil Science,2014,45(5):1265-1272.

[本文引用: 1]

邓绍坡,骆永明,宋静.

电子废弃物拆解地PM10中多氯联苯、镉和铜含量调查及人体健康风险评估

环境科学研究,2010,23(6):733-740.

[本文引用: 1]

Deng S P, Luo Y M, Song J,et al.

Concentrations of polychlorinated biphenyls,cadmium and copper in inhalable particulate (PM10) from an e⁃waste recycling area and human health risk assessment

Research of Environmental Sciences,2010,23(6):733-740.

[本文引用: 1]

Endo S, Grathwohl P, Haderlein S B,et al.

LFERs for soil organic carbon⁃water distribution coefficients (KOC) at environmentally relevant sorbate concen⁃trations

Environmental Science & Technology,2009,43(9):3094-3100.

[本文引用: 4]

Ulrich N, Endo S, Brown T N,et al.

UFZ⁃LSER database v 3

.2.1. Leipzig,GermanyHelmholtz Centre for Environmental Research⁃UFZ,2017.

[本文引用: 1]

Xia G S, Ball W P.

Adsorption⁃partitioning uptake of nine low⁃polarity organic chemicals on a natural sorbent

Environmental Science & Technology,1999,33(2):262-269.

[本文引用: 1]

Liu K, Fu H Y, Zhu D Q,et al.

Prediction of apolar compound sorption to aquatic natural organic matter accounting for natural organic matter hydrophobicity using aqueous two⁃phase systems

Environmental Science & Technology,2019,53(14):8127-8135.

[本文引用: 3]

石昕. 土壤不同组分对极性和非极性苯系化合物的吸附行为研究. 博士学位论文. 南京南京大学,2010.

[本文引用: 2]

Shi X. Investigating roles of different soil components in sorption of polar and nonpolar aromatic compounds. Ph.D. Dissertation. NanjingNanjing University,2010.

[本文引用: 2]

舒月红,黄小仁,贾晓珊.

1,2,4⁃三氯苯/四氯乙烯和

1,2,4,5⁃四氯苯在土壤/沉积物上的竞争吸附. 土壤学报,2009,46(4):714-720.

[本文引用: 2]

Shu Y H, Huang X R, Jia X S.

Competitive sorption between 1,2,4⁃Trichlorobenzene/Tetra⁃chloroethene and

1,2,4,5–Tetrachlorobenzene on soils/sediments. Acta Pedologica Sinica,2009,46(4):714-720.

[本文引用: 2]

何耀武,区自清,孙铁珩.

多环芳烃类化合物在土壤上的吸附

应用生态学报,1995,6(4):423-427.

[本文引用: 2]

He Y W, Ou Z Q, Sun T H.

Adsorption of poly⁃cyclic aromatic hydrocarbons on soils

Chinese Journal of Applied Ecology,1995,6(4):423-427.

[本文引用: 2]

Wang X L, Sato T, Xing B S.

Sorption and displacement of pyrene in soils and sediments

Environmental Science & Technology,2005,39(22):8712-8718.

[本文引用: 1]

Ghosh U, Zimmerman J R, Luthy R G.

PCB and PAH speciation among particle types in contaminated harbor sediments and effects on PAH bioavailability

Environmental Science & Technology,2003,37(10):2209-2217.

[本文引用: 1]

Song J Z, Peng P A, Huang W L.

Black carbon and kerogen in soils and sediments

1. Quantification and characterization. Environmental Science & Technology,2002,36(18):3960-3967.

[本文引用: 1]

霍丽娟,任理,毛萌.

阿特拉津及其代谢物在砂质壤土中的吸附

中国环境科学,2018,38(1):254-262.

[本文引用: 1]

Huo L J, Ren L, Mao M,et al.

Sorption of atrazine and its metabolites on a sandy loam soil

China Environmental Science,2018,38(1):254-262.

[本文引用: 1]

屈梦雄. 阿特拉津在东北冻融土壤中的吸附降解行为研究. 硕士学位论文. 大连大连理工大学,2014.

Qu M X. Adsorption and degradation behavior of atrazine in freezing and thawing soils of northeast China. Master Dissertation. DalianDalian Univer⁃sity of Technology,2014.

李金梅. 阿特拉津在狼尾草根际与非根际土壤的吸附特性及机制. 硕士学位论文. 哈尔滨东北农业大学,2018.

[本文引用: 1]

Li J M. Adsorption characteristics and mechanism of atrazine in rhizosphere and non⁃rhizosphere soil of pennisetum. Master Dissertation. HarbinNortheast Agricultural University,2018.

[本文引用: 1]

沈学优,卢瑛莹,朱利中.

对⁃硝基苯酚在水/有机膨润土界面的吸附行为——热力学特征及机理

中国环境科学,2003,23(4):367-370.

[本文引用: 1]

Shen X Y, Lu Y Y, Zhu L Z.

Sorption behavior of p⁃nitrophenol on the boundary of water and organobentonite:The thermodynamic character and mechanism

China Environmental Science,2003,23(4):367-370.

[本文引用: 1]

Wang X L, Guo X Y, Yang Y,et al.

Sorption mechanisms of phenanthrene,lindane,and atrazine with various humic acid fractions from a single soil sample

Environmental Science & Technology,2011,45(6):2124-2130.

[本文引用: 1]

Liu H J, Wei K Y, Yu Y S,et al.

Predicting adsorption coefficients of VOCs using polyparameter linear free energy relationship based on the evaluation of dispersive and specific interactions

Environmental Pollution,2019,255(Pt 1):113224.

[本文引用: 1]

Golbraikh A, Shen M, Xiao Z Y,et al.

Rational selection of training and test sets for the development of validated QSAR models

Journal of Computer⁃Aided Molecular Design,2003,17(2):241-253.

[本文引用: 1]

Ripszam M, Paczkowska J, Figueira J,et al.

Dissolved organic carbon quality and sorption of organic pollutants in the Baltic Sea in light of future climate change

Environmental Science & Technology,2015,49(3):1445-1452.

[本文引用: 2]

Niederer CGoss K USchwarzenbach R P.

Sorption equilibrium of a wide spectrum of organic vapors in leonardite humic acid:Modeling of experimental data

Environmental Science & Technology,200640(17):5374-5379.

[本文引用: 1]

Poole S K, Poole C F.

Chromatographic models for the sorption of neutral organic compounds by soil from water and air

Journal of Chromatography A,1999,845(1-2):381-400.

[本文引用: 2]

Vitha M, Carr P W.

The chemical interpretation and practice of linear solvation energy relationships in chromatography

Journal of Chromatography A,2006,1126(1-2):143-194.

[本文引用: 1]

Kile D E, Chiou C T, Zhou H D,et al.

Partition of nonpolar organic pollutants from water to soil and sediment organic matters

Environmental Science & Technology,1995,29(5):1401-1406.

[本文引用: 1]

Spurlock F C, Biggar J W.

Thermodynamics of organic chemical partition in soils. 2. Nonlinear partition of substituted phenylureas from aqueous solution

Environmental Science & Technology,1994,28(6):996-1002.

[本文引用: 1]

张良静,顾慧燕,毕二平.

含特殊官能团有机物分配系数预测的多参数线性自由能关系(PP⁃LFERs)研究

环境化学,2018,37(5):1045-1053.

[本文引用: 1]

Zhang L J, Gu H Y, Bi E P.

Polyparameter linear free energy relationships (PP⁃LFERs) for predicting partition coefficients of organic compounds with specific functional groups

Environmental Chemistry,2018,37(5):1045-1053.

[本文引用: 1]

/