1. 引言
胡萝卜是双子叶植物纲,伞形科,胡萝卜属的植物。为春、冬两季的主要蔬菜之一,享有“小人参”和“金笋”之美誉 [1]。近年来,由于农业种植产业结构调整和退耕还林后经济作物的发展,胡萝卜成为农民的主要种植产业和经济来源之一,种植面积超过1万公顷。但由于胡萝卜茎叶部分接近地面,根部深入土壤,容易引发微管蚜、茴香凤蝶和黑腐病、菌核病、软腐病等病虫害。据调查,通常发病率为10%~30%,严重的可以达到50%以上 [2]。敌敌畏(Dichlorvos)又名O,O-二甲基-O-(2,2-二氯乙烯基)磷酸酯,分子式为C4H7Cl2O4P,原药为无色至浅棕色液体,有挥发性,有芳香味,闪点75℃,熔点−60℃,密度1.42 g/cm3,沸点74℃,蒸气压1.7 mPa (25℃),水中溶解度10 g/L (室温),是一种有机磷杀虫剂,具有胃毒、触杀和熏蒸作用,残效期较短,并且残毒低,对害虫击倒力快而强。广泛应用于防治果蔬、蔬菜和茶叶等害虫。目前,敌敌畏在我国小麦、玉米、水稻、十字花科蔬菜等作物上的残留分析方法有GC、GC-MS等 [3] [4] [5] [6]。
随着生活品质的不断提升,各级政府及广大民众对食品安全问题越来越重视。食品特别是蔬菜中的农药残留问题成为国内外专家学者关注的重问题点之一。对于农产品中农药残留或农药降解规律的研究也取得了许多成果,这些研究工作主要体现在两个方面:第一,从定性方面出发研究农残产生的机理及危害,提出各种建议等,如 [2] [3] [4] [5] 等等。第二,从定量方面出发探索农药降解与残留的规律,依据不同的侧重点建立相应的数学模型,给出农药降解量或残留值随时间的演化方程,用以估计未来时刻的农药降解量或残留值,并给以安全建议。如文献 [7] 等。正如大家所熟知,农药的降解是一个复杂的物理化学过程,不同品种的蔬菜、不同类别的农药其降解机理不尽相同,并且大量研究表明农药的降解与环境有着密切相关 [7]。影响农药降解的环境因素主要是日照时长、温度、平均降雨量。1) 日照时长:日照时间的长短不仅对数植物的生长周期是一个影响因素因素,对于农药的降解速率也是一个非常重要的影响因素。光照在某些化学反应中能促进或抑制化学反应的过程,也是某些化学反应的必要条件。2) 温度:各种植物的生长、发育都要求有一定的温度条件,植物的生长和繁殖要在一定的温度范围内进行,在此温度范围的两端是最低和最高温度。低于最低温度或高于最高温度都会引起植物体死亡。最低与最高温度之间有一最适温度,在最适温度范围内植物生长繁殖得最好。同一植物在温度不同的条件下生长情况会有所不同,且温度的高低也会影响农药在植物体内的降解速度。3) 降雨量:水不仅是人类赖以生存的重要资源也是植物体的重要组成部分、更是植物生命活动不可缺少的重要环境因子之一。在植物体内,含水量可达70%至80%,—些水生植物的含水量可高达90%以上,可见水和植物的关系非常重要。水不单是植物体的重要组成部分,水也影响着体内酶的生化反应、植物的光合作用、有机物的代谢过程等。然而降雨量的多少会影响到土壤的PH值,干湿程度等,这些因素都会直接影响胡萝卜中水分的含量,胡萝卜中水分减少或增多会直接使得农药在植物中含量不同。因此降雨量对敌敌畏在胡萝卜中的残留值影响是一个非常重要的因素。
本文主要考虑环境因素对农药降解影响问题。通过敌敌畏在胡萝卜和土壤中的消解动态和最终残留试验,建立了不同区域敌敌畏在胡萝卜上的残留动态模型,获得了敌敌畏的残留演化过程,为有机磷农药在蔬菜上的合理使用提供科学参考。
2. 模型的设定与实证分析
2.1. 变量的选取及数据来源
本文选取解释变量:长沙市、石家庄市的平均日照时长(X1)、温度(X2)、降雨量(X3)。被解释变量:Y (敌敌畏在胡萝卜中的残留量)。解释变量的数据来源于2016年各省《统计年鉴》以及中国气象数据网,被解释变量数据来源于由农业部农药所审批,湖南农业大学农业环保研究所承担敌敌畏在胡萝卜上的最终残留水平的补点试验。
2.2. 模型的构建
为研究上述变量对敌敌畏在胡萝卜上农药残留的影响程度,建立如下回归模型:
采用OLS (ordinary least square)方法估计模型,利用SPSS分析各因素对敌敌畏在胡萝卜上的降解过程和残留影响。分别作长沙市和石家庄市Y (敌敌畏在胡萝卜中的残留量)对长沙市和石家庄市平均日照时长(X1)、温度(X2)、降雨量(X3)的曲线回归(包括线性、二次、立方、复合函数),选取拟合效果最好的曲线回归。长沙市各变量的R2结果如表1所示,石家庄市各变量的R2结果如表2所示。
![](Images/Table_Tmp.jpg)
Table 1. R2 value of each variable in Changsha
表1. 长沙市各变量R2值
![](Images/Table_Tmp.jpg)
Table 2. R2 value of each variable in Shijiazhuang
表2. 石家庄市各变量R2值
由上表可知,分别作长沙市和石家庄市Y (敌敌畏在胡萝卜中的残留量)对平均日照时长(X1)、温度(X2)、降雨量(X3)的曲线回归中拟合效果最好的均是三次函数。且平均日照时长(X1)和温度(X2)跟Y (敌敌畏在胡萝卜中的残留量)在无论什么函数形式的拟合度都很差,说明日照时长和温度对敌敌畏在胡萝卜中的残留量相关性很小,故只选取降雨量(X3)作为自变量。所以根据试验所得数据,分别对长沙市、石家庄市两地做曲线回归分析。长沙市的结果如表3所示,石家庄市的结果如表4所示。
![](Images/Table_Tmp.jpg)
Table 3. Changsha’s curve regression result
表3. 长沙曲线回归结果
![](Images/Table_Tmp.jpg)
Table 4. Shijiazhuang’s curve regression result
表4. 石家庄曲线回归结果
长沙市降雨量对农药残留量影响的拟合图如图1所示,石家庄市降雨量对农药残留量影响的拟合图如图2所示:
![](//html.hanspub.org/file/1-2621042x10_hanspub.png)
Figure 1. Fitting map of the impact of rainfall on pesticide residues in Changsha
图1. 长沙市降雨量对农药残留量影响的拟合图
![](//html.hanspub.org/file/1-2621042x11_hanspub.png)
Figure 2. Fitting map of the impact of rainfall on pesticide residues in Shijiazhuang
图2. 石家庄市降雨量对农药残留量影响的拟合图
所以长沙市降雨量与敌敌畏在胡萝卜中的残留量的函数模型为:
所以长沙市降雨量与敌敌畏在胡萝卜中的残留量的函数模型为:
3. 结果与分析
本文从农药残留产生的机理入手,分析了降雨量、温度、日照时长等环境因素对敌敌畏在胡萝卜中的残留量产生的影响。结果表明:不同区域的平均温度和平均日照时长对敌敌畏在胡萝卜中的残留量的影响很小,而不同区域的降雨量对敌敌畏在胡萝卜中的残留量的影响显著。在此基础上建立了降雨量与敌敌畏在胡萝卜中的残留量关系的非线性回归模型,为农药在蔬菜中的残留量估计及安全合理使用提供理论依据。