我国是化肥农药使用大国。根据联合国粮农组织相关数据，2016年我国化肥使用量为503.32kg/hm2，是世界平均使用量的3.58倍。同时，2014—2016年我国农药年均使用总量大致为170万t，高居世界第一[1,2]。化肥农药使用在中国粮食连续丰收增产过程中发挥了巨大作用，但由于使用方法不当等因素，滥用现象越发严重，对生态环境的危害日趋明显[3]。长期过量使用化肥易造成土壤板结、养分结构失调、肥力下降等问题[4]。同时，农田中施用的农药量仅有30%左右被作物利用，其余均扩散至土壤和大气中，土壤中农药残留量增加，造成土壤污染[5]。因此，农业部于2015年提出“到2020年农药化肥零增长行动”，对化肥农药的使用提出更高要求，化肥农药的混配施用成为解决这一问题的有效途径。

吡虫啉是1985年德国拜耳公司和日本特殊农药株式会共同研发的第一个商业化的新烟碱类农药，占整个新烟碱类农药市场的三分之一以上，中国是世界上吡虫啉最大的生产国与出口国。吡虫啉可以稳定存在于自然界中，在土壤表面具有约39d的光解半衰期，可长期残留于土壤[6,7]。近年来研究发现吡虫啉残留对多种非靶标生物包括蜜蜂、鸟类和哺乳动物等均有不同程度的危害[8,9,10]。因此研究吡虫啉在土壤中相关环境行为非常重要。吸附特性作为其主要环境行为之一，之前已有相关报道，包括在不同土壤中的吸附特性、生物炭对其在土壤中吸附特性的影响等[11,12]。尿素作为氮肥的主要形式，占据我国氮肥产量的60%以上。由于尿素利用率低，每年各种途径流失掉的尿素达到总使用量的一半左右，对生态环境造成严重影响，因此提高尿素利用率是达成化肥零增长的重要环节。

吡虫啉与肥料的混合施用目前已有研究。吡虫啉悬浮种衣剂与氨基酸液肥的互配在一定比例浓度下可有效地提高对棉花棉蚜的防治，同时对棉花的成苗率、株高、叶绿素等均无负面的影响[13]。吡虫啉与氮肥联合施用对番茄幼苗上的烟粉虱也具有良好的防治效果[14]。表明吡虫啉与肥料的混施是具有可行性的。

甘蔗是我国重要的工业制糖原料，2018年种植面积达到120万hm2，蔗糖产量944.5万t，占食糖总产的88%。甘蔗主要害虫包括甘蔗螟虫、甘蔗金龟子、甘蔗蚜虫等，吡虫啉均对其有较好的防治效果[15,16,17]。吡虫啉在我国粤西蔗区主要土壤类型中的吸附特性尚无研究，因此本研究选择粤西蔗区主要土壤类型红壤土进行研究，旨在通过研究尿素与吡虫啉混施对吡虫啉吸附特性的影响，为今后尿素与吡虫啉在粤西蔗区混配施用建立一定的理论基础。

1、材料与方法

1.1试验材料

吡虫啉标准品（纯度>98%）、乙腈（色谱纯）、尿素、氯化钙均为分析纯。供试土壤为红壤土，于2018年7月26日采自广东省湛江市的甘蔗农田表面土壤，取样深度为0～20cm，除杂、自然风干，过60目筛网备用，供试土壤为红壤土，基本性质：有机质32g/kg，有机碳18.6g/kg，阳离子交换量6.8cmol(+)/kg，机械组成：沙砾19.1%、粉砾13.4%、粘砾67.5%。

1.2试验仪器

IS-RSV1型立式恒温振荡器（美国精骐有限公司）、Agilent1260InfintyⅡ高效液相色谱仪（安捷伦科技有限公司）、H13-18K台式高速离心机（湖南可成仪器设备有限公司）、KQ-500B超声波清洗器（昆山市超声仪器有限公司）。

1.3试验方法

1.3.1吸附动力学试验

称取5g供试土壤于50mL锥形瓶，加入25mL浓度为1mg/L的吡虫啉溶液（背景溶液为0.01mol/LCaCl2），同时分别加入不同量尿素使尿素浓度达到0.8、4和8g/kg。不添加尿素为空白处理。将所有锥形瓶封闭瓶口后置于恒温振荡器于25±1℃、180r/min条件下振荡24h。分别于0、1、2、4、8、12、16、24h采样。将采集的样品置于离心机中以3500r/min离心7min，收集上清液过0.22µm水系微孔滤膜，测定吡虫啉含量。各处理3个重复。

1.3.2等温吸附试验

称取5g供试土壤于50mL锥形瓶，分别加入25mL浓度为0.5、1、2、10、20mg/L吡虫啉水溶液（背景溶液为0.01mol/LCaCl2溶液），同时分别加入不同量尿素使尿素浓度达到0.8、4和8g/kg，不添加尿素为空白处理。将所有锥形瓶封闭瓶口后置于恒温振荡器于25±1℃、180r/min条件下振荡24h，于24h采样。后续步骤同1.3.1。各处理重复3次。

1.3.3吸附热力学试验

设置摇床振荡温度为35±1℃、45±1℃，其他步骤同1.3.1操作进行。

1.4数据分析

1.4.1吸附量与吸附率的计算

根据我国《化学农药环境安全评价试验准则》吸附率的计算公式见式（1），吸附量计算公式见式（2）。

（公式）

式中：A-吸附率，%；M-农药总质量，mg;Cs-吸附平衡时土壤中吡虫啉的浓度，mg/kg;Ce-吸附平衡时溶液中吡虫啉的浓度，mg/L;C-吸附初始溶液中吡虫啉的浓度；V-溶液体积，mL;m-土壤质量，kg。

1.4.2吸附动力学模型

分别采用准一级吸附动力学方程（3）、准二级吸附动力学方程（4）和Elovich方程（5）对吸附动力学过程进行拟合分析。

（公式）

式中：Ct-土壤对吡虫啉在t时刻的吸附浓度，mg/kg;Cs-吸附平衡时土壤中吡虫啉浓度，mg/kg;t-吸附时间，h;K1-准一级吸附动力学方程速率常数；K2-准二级吸附动力学方程速率常数。

（公式）

式中：Ct-土壤对吡虫啉在t时刻的吸附浓度，mg/kg;t-吸附时间，h;a、b-常数。

1.4.3吸附等温线模型

分别采用Freundlich等温吸附模型（6）和Langmuir等温吸附模型（7）来描述吸附平衡时土壤中吡虫啉浓度和溶液中吡虫啉浓度的关系。

（公式）

式中：Cs-吸附平衡时土壤中吡虫啉浓度，mg/kg;Ce-吸附平衡时溶液中吡虫啉浓度，mg/L;Kf-Freundlich模型吸附常数；n-吸附强度常数；Qm-吸附平衡时土壤中吡虫啉最大吸附浓度，mg/kg;KL-Langmuir模型吸附亲和力常数。

1.4.4吸附热力学常数

吸附热力学过程通过式（8）、（9）、（10）计算得吉布斯自由能变ΔG、熵变ΔS和焓变ΔH进行分析。

（公式）

2、结果与分析

2.1尿素对吡虫啉吸附动力学的影响

尿素对吡虫啉吸附动力学过程的影响如图1所示。不同氮肥浓度下，吡虫啉吸附动力学过程保持一致。与CK相比，添加尿素后，吡虫啉吸附过程为3个阶段。0～1h，快速吸附阶段，尿素添加不影响快速吸附阶段的吸附能力，吸附量无明显差异；1～16h，慢速吸附阶段，吡虫啉在各处理下吸附量均缓慢增长，同时吸附量开始出现明显差异，随尿素浓度升高，吸附量逐渐降低；16～24h，吸附平衡阶段，吡虫啉在红壤中的吸附基本达到平衡。吸附平衡时，随尿素浓度（0.8、4、8g/kg）升高吡虫啉吸附量分别下降了6.00%、13.52%和20.41%，各处理之间具有明显差异。结果表明，尿素添加可以降低吡虫啉在红壤中的吸附量且随尿素浓度升高吸附量逐渐降低。

图1吡虫啉吸附动力学曲线

尿素对吡虫啉吸附动力学过程的影响拟合结果如图2(a、b、c）所示，拟合参数如表1所示。由表1可知：准一级动力学方程拟合相关系数随尿素浓度升高分别为0.684、0.871、0.862和0.863，拟合系数均小于0.900，表明添加尿素之后，准一级动力学方程不能较好描述吡虫啉在红壤中的吸附动力学过程。准二级动力学方程拟合相关系数随尿素浓度升高分别为0.904、0.937、0.931和0.924，可以较好地描述吡虫啉吸附动力学过程，吸附速率常数K2随尿素浓度的升高分别为0.591、0.923、0.813和0.955，吸附速率逐渐升高但各浓度与CK之间无显著差异，拟合吸附量Ce分别为2.005、1.900、1.875和1.729，吸附量逐渐降低，这与实际吸附容量结果一致；同时在高浓度尿素8g/kg处理条件下，吸附量与CK及低浓度尿素处理下均具有显著差异。上述结果表明，随尿素浓度升高，吡虫啉在红壤中的吸附速率逐渐上升，同时吸附量逐渐下降。Elovich方程拟合相关系数分别为0.925、0.974、0.966和0.949，也可以较好地描述吡虫啉吸附动力学过程。吸附常数b随尿素浓度升高分别为1.193、1.210、1.201和1.133，各处理之间无明显差异。吸附常数a随尿素浓度升高分别为0.288、0.265、0.225和0.200，吸附常数逐渐降低，同时高浓度尿素8g/kg处理条件下，与CK及低浓度尿素处理下均具有显著差异。上述结果表明，随尿素浓度升高，在快速吸附阶段，尿素对吡虫啉在土壤表面的物理吸附无显著影响。在慢速吸附阶段，吡虫啉在红壤中内部化学扩散过程受尿素浓度升高产生降低吸附的影响。

图2不同尿素浓度处理吡虫啉吸附动力学拟合曲线

表1动力学吸附模型拟合参数

2.2尿素对吡虫啉等温吸附特性的影响

吡虫啉与不同尿素浓度吸附等温过程拟合结果如图3(a、b、c）所示，拟合参数如表2所示。由表2可知：Freundlich和Langmuir拟合模型均可以较好地描述吸附等温过程，相关系数R2≥0.933。Freundlich拟合模型中，在25℃条件下，随尿素浓度升高吸附常数Kf分别为2.700、2.584、2.158和1.852,Kf逐渐下降，较高尿素浓度4g/kg和8g/kg处理下，与CK和低尿素浓度处理相比差异显著。表明高浓度尿素处理下，吡虫啉在红壤中的吸附能力逐渐降低。与CK相比，添加尿素后线型分子n增大且与CK有显著差异，表明添加尿素之后，吡虫啉在红壤中的吸附非线性程度逐渐升高。在35℃条件下，随尿素浓度升高吸附常数Kf分别为2.381、2.205、2.378、2.104,Kf下降但无显著影响，表明在35℃条件下氮肥浓度对吡虫啉的吸附无显著影响。与CK相比，线型分子n增大，但仅在4g/kg处理下差异显著。表明35℃条件下，尿素对吡虫啉的吸附线性程度无显著影响。在45℃处理下，结果与25℃处理下结果相似，吸附常数Kf逐渐降低，且在高尿素浓度8g/kg处理下差异显著。线性分子n逐渐增大，同样在8g/kg处理下差异显著。表明，随尿素浓度升高吡虫啉吸附在红壤中的吸附能力逐渐降低，非线性程度逐渐升高。上述结果表明，在25℃和45℃时，高浓度尿素处理显著降低吡虫啉在红壤中的吸附能力，同时提高吡虫啉吸附的非线性程度。Langmuir拟合模型中，25℃和35℃处理下，Langmuir吸附常数KL在不同处理中差异显著，45℃处理下，Langmuir吸附常数KL无显著差异。表明尿素对吡虫啉在红壤中单分子层吸附速率在较低温度时产生影响，在高温时无显著影响。在25℃和45℃处理条件下，拟合最大吸附容量Qm均随尿素浓度的升高逐渐降低且差异显著，在35℃处理条件下，尿素浓度对Qm无显著影响。表明在25℃和45℃条件下，吡虫啉吸附能力逐渐降低是由于土壤对吡虫啉吸附容量降低造成的，35℃条件下无明显影响，这一结果与Freundlich拟合结果一致。

图3不同尿素处理吡虫啉等温吸附特性拟合曲线

表2等温拟合模型特征值

2.3尿素对吡虫啉吸附热力学的影响

尿素对吡虫啉吸附热力学的影响如表3所示。红壤中添加尿素之后，吉布斯自由能ΔG数值无明显差异，表明添加尿素后，吡虫啉在红壤中的吸附仍然以物理吸附为主，尿素无显著影响，这与Elovich吸附动力学拟合结果一致，同时ΔG均小于0，表明吸附仍然是自发进行的。添加尿素之后，吸附焓变ΔH在高浓度尿素处理下升高，但ΔH均小于0，表明尿素添加后，吡虫啉在红壤中的吸附仍然为放热过程，氮肥的添加对放热过程产生影响但无显著差异。添加尿素后，吸附熵变ΔS均大于0且在高浓度尿素下逐渐升高，表明添加低浓度尿素后，对吡虫啉在红壤中的吸附混乱度无显著影响；添加高浓度尿素后，吡虫啉在红壤中的吸附混乱度升高。

表3等温拟合模型特征值

3、讨论

吡虫啉在土壤中的吸附动力学过程分为3个阶段，快速吸附、慢速吸附和平衡阶段。第一阶段是快速吸附阶段，吡虫啉与土壤进行接触产生振荡反应，吡虫啉很快就被吸附到土壤表面，吸附量快速增加；第二阶段是慢速吸附阶段，这一阶段主要是因为土壤中较易吸附的位点已经被占据，吡虫啉必须进入土壤颗粒内部不太容易吸附的位点，使得吸附速率开始降低；第三个阶段则是吸附平衡阶段。这一吸附过程与之前已进行的相关研究相吻合[18,19]。吸附动力学方程更加符合准二级动力学方程，准二级动力学方程包括物理吸附和化学吸附，说明吡虫啉的吸附不是单一方式的。Elovich方程是经验方程，假设吸附是一个先快后慢的过程，与本试验相符，也能较好地进行拟合。同时，随着尿素浓度的升高，吡虫啉在土壤的吸附量逐渐下降，这可能是因为吡虫啉分子中带有的极性集团使其较易与土壤中的有机质等形成氢键而发生吸附[20,21]，同时尿素进入土壤中后，小部分尿素分子通过氢键被土壤吸附，二者之间可能存在着竞争关系。并且尿素的添加可以改变土壤的酸碱性，在酸性土壤中加入尿素，土壤的pH都会在短期内随着尿素浓度的增大而急剧上升[22]，已有研究表明随着土壤pH的升高，吡虫啉在土壤中的吸附下降[23]。因此，添加尿素造成土壤pH的上升及尿素与吡虫啉之间的竞争吸附，可能是导致吡虫啉吸附量降低的原因。

本实验选择Freundlich和Langmuir两种等温吸附模型，两种模型均可以较好地进行拟合。Langmuir模型假设吸附是单层吸附，吸附是均匀分布的，Freundlich模型则是多分子层的吸附公式，用以说明不均匀的吸附过程。本实验数据同时符合两种吸附模型，这可能是由于土壤的结构复杂，包含多种离子和官能团，尿素的加入，土壤中加入更多不同分子和离子，可能导致土壤结构进一步复杂化，因此吸附过程既包括单层吸附也包括多层吸附[24]。在此之前也有过类似的报道，如Mg-Al和Zn-Al吸附去除磷酸盐时，便同时满足两种模型[25]。同时满足两种吸附模型还表明吸附过程当中物理吸附和化学吸附都是存在的。这与动力学试验得到的结果一致。Freundlich方程中吸附系数Kf反应土壤对农药的吸附强度，根据农药土壤吸附特性等级划分Kf<5属于难吸附，本试验中所有Kf均小于5（表2、表3），说明吡虫啉在红壤中均属于难吸附。Freundlich方程中的模型指数1/n表示农药在土壤中的吸附线性程度，当1/n<1时，表示吸附为非线性，同时数值越靠近1则非线性程度越低[26]。本试验当中所有模型指数均小于1，这种吸附特征曲线表明在低平衡浓度下，农药分子量较少，吸附位点较多，吸附是非线性；而高平衡浓度下，农药分子量增多，吸附位点减少，则主要是线性吸附[27,28]。当尿素浓度处于较高水平时，模型指数越靠近1，这可能是由于尿素浓度逐渐增高时，占据的吸附位点更多，即使农药在低平衡浓度时，吸附位点也减少，吸附的非线性程度也降低，因此尿素浓度较高时，吡虫啉在土壤中吸附的非线性程度降低。

所有处理中吉布斯自由能ΔG均为负值，表明吸附是一个自发过程。熵变ΔS均为正值，表明吸附过程是体系自由度增大，能量升高的过程[29]。无尿素添加的吡虫啉在土壤中的吸附，吸附熵ΔH>0,表明吸附过程是吸热反应，随着温度的升高，吡虫啉在土壤中的吸附逐渐升高，这与温度对吸附的影响一致，但是增加程度并不明显，这与他人的研究情况并不完全一致。之前的研究表明随着土壤温度的升高吡虫啉的吸附能力逐渐降低，BROZNIC等研究了克罗地亚沿海地区土壤当温度由20℃升高至40℃时，吡虫啉Freundlich常数由4.92降低至3.17[30]，这可能是由于不同土壤及其理化性质不同造成的结果。

4、结论

(1）吡虫啉在土壤中吸附动力学符合准二级动力学方程和Elovich方程，吸附过程是一个先快后慢的过程，添加尿素仍符合准二级动力学方程和Elovich方程，但随着尿素使用量的增大，吸附率和吸附量均出现下降，表明尿素加入降低了吡虫啉的吸附性。

(2）吡虫啉的吸附过程同时满足Freundlich和Langmuir方程，说明吸附过程同时存在单层吸附和多层吸附，尿素的加入降低了吸附的非线性程度。

(3）整个吸附过程ΔG<0、ΔH<0、ΔS>0，表明吸附是自发进行混乱度不断升高的放热反应。

参考文献：

[1]刘华军,孙淑惠,李超.环境约束下中国化肥利用效率的空间差异及分布动态演进[J].农业经济问题,2019(8):65-75.

[2]朱春雨,杨峻,张楠.全球主要国家近年农药使用量变化趋势分析[J].农药科学与管理,2017,38(4):13-19.

[3]谢邵文,杨芬,冯含笑,等.中国化肥农药施用总体特征及减施效果分析[J].环境污染与防治,2019,41(4):490-495.

[4]肖军,秦志伟,赵景波.农田土壤化肥污染及对策[J].环境保护科学,2005(5):36-38.

[5]高寒,陈娟,王沛芳,等.农药污染土壤的生物强化修复技术研究进展[J].土壤,2019,51(3):425-433.

[6]邴启征.木犀草素治疗吡虫啉诱导的鹌鹑肝脏损伤的效果观察[D].东北农业大学,2019:63.

[11]何丽芝,张小凯,吴慧明,等.生物质炭及老化过程对土壤吸附吡虫啉的影响[J].环境科学学报,2015,35(2):535-540.

[12]宣日成,王琪全,郑巍,等.吡虫啉在土壤中的吸附及作用机理研究[J].环境科学学报,2000(2):72-75.

[13]姜伟丽,付传翠,王丹,等.“新壮态”液肥对吡虫啉拌种防控棉蚜的增效作用及其对棉花生长的影响[J].中国棉花,2018,45(3):13-15.

王鑫,和锐敏,赖鉴添,邵丹青,何琴,付建涛,卢颖林,钟国华,安玉兴.粤西蔗区土壤中尿素对吡虫啉吸附特性的影响[J].中国糖料,2021(02):47-54.