论文部分内容阅读
摘要 为提高农药利用率，研究山地果园管道恒压喷雾系统中喷雾压力和孔径对雾滴粒径参数的影响，利用激光粒度仪，测量3种孔径空心圆锥雾喷头在8种压力下的6种雾滴粒径参数数据，分析各粒径级雾滴的分布情况，给出雾滴粒径大小及其随压力和孔径的变化趋势，进行了雾滴参数的多元线性回归，建立基于压力的雾滴参数模型。结果表明：粒径大于100 μm的雾滴（小于0.3%）和粒径为0～20 μm的雾滴（小于0.5%）可忽略不计，雾滴粒径主要分布于20～40 μm（79.5%～92.8%）；雾滴均较细小，全部为气溶胶；孔径越小，压力越大，雾滴越细小；6种雾滴粒径参数与孔径和压力均有良好的二元线性关系（R2分别为0.928、0.956、0.949、0.949、0.889和0.815）；6种雾滴粒径参数均随压力呈二次多项式变化规律，决定系数R2均达到了0.894以上。研究结果对山地果园管道恒压喷雾中喷头的选型、喷雾压力的调整及喷雾效果的优化有重要参考意义。
关键词 山地果园； 管道恒压喷雾； 雾滴粒径参数； 喷雾压力； 孔径
中图分类号： S 491
文献标识码： A
雾滴粒径参数是反映雾滴大小的主要指标[12]，直接影响雾滴的飘移[3]、沉积[47]和药效的发挥[8]，从而影响农药的利用率[9]和喷雾效果[1011]。国内外学者采用各种方法对雾滴粒径进行测量[12]，用图像处理技术对雾滴进行显微识别[13]，利用高速摄像机结合数字图像处理技术对雾滴粒径分布和雾滴的运动进行分析[1415]，利用激光粒度仪测量雾滴谱[16]，同时对影响喷头雾滴粒径参数的各种因素进行了分析[1718]。但这些研究多为测量在较低压力喷雾时喷头的雾滴，大多集中在单因素对雾滴粒径的影响分析，而对适合果园需求的高压力喷雾时喷头的雾滴大小研究较少，主要影响因素与雾滴大小之间的相关关系研究也较少。
我国南方果园多为山地果园，多采用低矮密植型种植模式，移动式喷雾机械很难进入[1920]，为此笔者所在团队研究了果园管道恒压喷雾控制系统。由于管道喷雾无需机具在果园中移动，具有省力、省工、喷雾效果好等优点，适合山地果园使用[21]。管道喷雾依靠液泵形成的较高压力，将药液在喷头处进行雾化并施向靶标，但管道内药液的压力随作业人数的不同会产生波动，同时，药液在管道中流动时会出现压力损失，导致整个管网系统药液压力不稳定且各点压力各不相同。管道恒压控制装置可根据管道中压力的实际值与设定值间的误差及误差变化趋势，在线调整模糊PID的参数，经带有变速积分、微分先行优化算子的增量式PID算法计算，获得控制量以控制管道中药液的压力[22]。因此压力成为影响管道喷雾系统中喷头雾滴参数的最重要的因素之一。除压力外，在喷头类型选定后，喷头孔径也是影响雾滴粒径参数的重要因素。然而到目前为止，还未见有文献对管道恒压喷雾系统中应如何设置喷雾压力及如何选择喷头的孔径进行研究。
本文通过改变压力和孔径，利用激光粒度仪采集雾滴粒径参数数据，研究压力和孔径对果园用空心圆锥雾喷头雾滴大小的影响，掌握雾滴大小随压力和孔径的变化规律，可以根据喷施对象的所需雾滴大小设置管道恒压喷雾控制装置的压力，选择合适的喷头类型，从而提高喷雾效果。
1 描述雾滴粒径的参数
（1）VMD：体积中值粒径（volume median diameter，VMD）。指取样雾滴的体积按雾滴大小顺序进行累计，其累计值为取样雾滴体积总和的50%所对应的雾滴直径。根据雾滴群VMD的大小将雾滴分为4类：①粗雾：VMD≥400 μm；②细雾：100 μm≤VMD<400 μm；③弥雾：50 μm≤VMD<100 μm；④气溶胶：VMD<50 μm。
（2）VAD：体积平均粒径（volume average diameter，VAD）。指取样雾滴平均体积所对应的直径。
（3）S/V：体积比表面积。指单位体积雾滴的表面积。
（4）SMD：表面积平均粒径（superficial area median diameter，SMD）。指粒径对表面积的加权平均。
（5）NAD：数量平均粒径（number average diameter，NAD）。指取样雾滴群的直径之和与雾滴群个数之和的比值。
（6）NMD：数量中值粒径（number median diameter，NMD）。指取样雾滴的个数按雾滴大小顺序进行累计，其累计值为取样雾滴个数总和的50%所对应的雾滴直径。
2 数据获取
雾滴粒径参数数据通过激光粒度仪获取。激光粒度仪通过测量颗粒群的散射谱，散射谱的强度与被测颗粒群的大小有关，并经透镜再次汇聚后，被位于透镜后焦面上的光电阵列探测器所接收，转换成电信号后经放大和A/D转换经通讯口送入计算机，进行反演运算和数据处理后，即可给出被测颗粒群的雾滴粒径参数。
2.1 试验设备
2.1.1 喷雾性能综合试验平台
喷雾性能综合试验平台主要由机械结构、电气控制和数据采集处理3部分组成。可提供本试验中所需雾滴并测量雾滴粒径参数。
2.1.2 激光粒度仪
Winner318B激光粒度仪可测定雾滴粒径参数。该仪器采用信息光学原理，通过测量颗粒群的散射谱，来分析雾滴粒径大小。其准确性误差<3%；重复性误差<3%；绿色激光器，波长532 nm，功率30 mW；量程为15～711 μm。
2.1.3 拉丝模孔径测量仪
拉丝模孔径测量仪可测量喷头孔径。测量仪将光学与数码、软件相结合，在软件中对孔径进行测量，可得出圆孔的直径。测量仪的分辨率为0.1 μm；重复精度为±0.5 μm；测量范围为0.05～30 mm；变倍方式为连续变倍，最大倍率为800倍。 2.2 试验对象
对Lechler TR80系列3种型号TR8001C、TR8002C和TR8003C喷头的雾滴粒径参数进行测量，喷头内部结构图如图1所示。该系列喷头压力范围为0.3～2.0 MPa，是一种带陶瓷芯的空心圆锥雾喷头，高耐磨，耐化学药剂，有极细雾滴谱，适用于压力要求较高的山地果园管道喷雾系统。
该系统在果园埋下塑料管网，管网进液端配有药池、恒压喷药控制系统，埋设的管道每隔一段距离就安装上升至地面的立式管，每个立管在打药时再与软管、喷枪相连接，使用电动机带动药泵加压，把药液输送到果树旁，人工逐株打药，适合多人同时打药，如图2所示。恒压喷药控制系统由变频器、变频电动机、药泵、压力变送器、恒压控制箱组成，可使出水口压力保持恒定，出水口压力可根据本文研究成果进行合理设定。
2.3 试验方法
试验用水为清水；用喷雾性能综合试验平台调节压力在0.70～1.40 MPa间变化，每次步进0.10 MPa；用拉丝模孔径测量仪对喷头孔径进行测量，测得TR8001C的喷嘴孔径为1.0 mm，TR8002C为1.4 mm，TR8003C为1.8 mm；用激光粒度仪对喷头中轴线上距喷嘴15 cm处的雾滴颗粒群的散射谱进行测量，得到雾滴能谱数据，再将能谱数据转换为雾滴粒径参数数据。
试验前用标准物质（直径为25 μm）对激光粒度仪进行标定。每次测量前先进行背景测试，背景测试累计10次后再进行样品测试以去除背景光的影响。每条样品测试数据为某个时刻样品的雾滴参数。本试验在每种工况下对样品连续测量100次得到100条测试数据，100条测试数据取平均值后作为每种工况下样品的测试数据。
3 数据处理与分析
首先将试验获得的能谱数据转换为雾滴粒径数据，然后分析各粒径级雾滴的分布情况，分析雾滴粒径大小，给出雾滴参数随压力和孔径的变化趋势，并对雾滴参数进行多元线性回归，最后建立基于压力的雾滴参数模型。
能谱数据向雾滴粒径数据转换时产生的计算误差称拟合误差。本试验拟合误差均小于0.012，符合试验要求。
本试验利用Origin 9.0进行绘图及曲线拟合，利用SPSS 16.0对试验数据进行多元线性回归分析，利用Excel 2007进行数据处理。
3.1 各粒径级雾滴的分布情况
根据3种喷头8种压力下雾滴的粒径分布数据，绘制各粒径级雾滴百分比的分布如图3所示。由于粒径D≥100 μm的雾滴很少，均少于0.3%，可忽略不计，所以仅绘制0～100 μm范围内雾滴的分布情况。可以看出：
（1） 雾滴均较细小。粒径为0～20 μm的雾滴比例均少于0.5%，可以忽略不计；主要是粒径为20～40 μm的雾滴，其比例均多于79.5%；粒径为40～60 μm的雾滴较少，其比例均少于10.6%；粒径为60～80 μm的雾滴较少，其比例均少于7.4%；粒径为80～100 μm的雾滴较少，其比例均少于2.1%。
（2） 按雾滴所占百分比从大到小排列为20～40 μm> 40～60 μm>60～80 μm>80～100 μm。
（3） 压力越大，孔径越小，雾滴粒径分布越集中，雾滴越细小。
3.2 雾滴粒径参数随压力和孔径的变化
为详细了解雾滴粒径参数的变化趋势，分别绘制出6种雾滴粒径参数VMD、VAD、S/V、SMD、NAD和NMD随压力和孔径变化的曲线，如图4所示。
从图4可以看出：
（1） 压力在0.70～1.40 MPa变化时，TR8001C、TR8002C和TR8003C的VMD变化范围分别为30.6～31.9 μm、30.9～32.0 μm和31.0～32.317 μm。说明在试验压力范围内，雾滴均较细小且全部为气溶胶。
（2） 孔径不变、压力越大，VMD、VAD、SMD、NAD和NMD越小，S/V越大。说明压力越大，雾滴的粒径越小。原因是压力越大，液体流出喷孔的初速度就越大，喷孔中液体扰动程度及喷出喷孔后所受到的空气介质阻尼也越大，雾滴更易破碎。
（3） 压力不变，孔径越小，VMD、VAD、SMD、NAD和NMD越小，S/V越大，即雾滴的粒径越小；原因是喷孔直径越小，液体喷出喷孔时会受到更大的节流，喷孔内扰动增加，雾滴更易破碎。
（4） VMD和VAD随压力和孔径的变化趋势较NAD和NMD更明显。说明：侧重反映体积的粒径参数比侧重反映数量的粒径参数随压力和孔径的变化趋势更明显。
3.3 基于压力和孔径的多元回归统计模型
雾滴的粒径大小与压力和孔径密切相关。从图4可以看出，6种雾滴粒径参数的变化趋势与压力和孔径的变化趋势基本一致，具有很强的相关性。其中，压力取值范围为0.80～1.40 MPa，步进0.10 MPa；孔径取值范围为1.0～1.8 mm，步进0.4 mm。分析结果如表1～3所示，可以看出：
（1）每种雾滴参数多元回归统计模型如表1所示，其中：y为各种雾滴粒径参数，x1为喷头孔径（单位为mm），x2为喷雾压力（单位为MPa）。模型的复相关系数R、决定系数R2、调整决定系数R2C均较高，其中VMD、VAD、S/V、SMD和NAD的决定系数R2分别为0.928、0.956、0.949、0.949、0.889和0.815。说明回归方程拟合优度较高，符合拟合要求。
（2） 由F检验可知，所有雾滴参数均有显著值P=0.000<0.01，表明VMD、VAD、S/V、SMD、NAD及NMD与孔径和压力之间均存在极显著的线性回归关系（见表2）。
（3） 由t检验可知，VMD、VAD、S/V、SMD及NAD的回归系数b0、b1和b2均有显著值P=0.000<0.01，表明方程的回归系数b0项（常数项）、b1项（孔径系数）和b2项（压力系数）影响均极显著；NMD的回归系数b0、b2均有显著值P=0.000<0.01，b1有显著值P=0.032<0.05，表明方程的回归系数b0项和b2项影响均极显著，回归系数b1项影响显著（见表3）。 综上，回归通过了F检验和t检验，且拟合优度较高，说明本文建立的雾滴粒径参数与压力和孔径的多元回归模型是准确的。
3.4 雾滴粒径参数与压力的关系
经曲线回归依次分析VMD、VAD、S/V、SMD、NAD及NMD随压力的变化规律。在进行多项式回归时，同时利用线性函数、指数函数、对数函数和幂函数进行了回归分析，这些函数的相关系数均比多项式回归的相关系数要小，故采用多项式表示各雾滴参数随压力的变化规律（见表4）。以喷雾压力（单位MPa）为横坐标，雾滴参数为纵坐标，得3种型号喷头的6种雾滴参数与喷雾压力（0.80～1.40 MPa）的拟合曲线如图5所示。从表4和图5可以看出，雾滴参数随压力的变化规律均能用二次多项式表示，且回归均具有较高的决定系数，R2均达到了0.894以上。
4 结论
Lechler TR8001C、TR8002C和TR8003C型喷头在0.70～1.40 MPa压力下，粒径D≥100 μm的雾滴均少于0.3%，粒径为0～20 μm的雾滴均少于0.5%，可忽略不计，按雾滴所占百分比从大到小排列20～40 μm>40～60 μm>60～80 μm>80～100 μm；主要是粒径为20～40 μm的雾滴（79.5%～92.8%）；VMD为30.6～32.3 μm，VAD为31.3～36.9 μm，雾滴均较细小，全部为气溶胶；压力越大、孔径越小，雾滴粒径越小。
VMD、VAD、S/V、SMD、NAD及NMD与孔径和压力均有较好的二元线性关系，建立了这6个雾滴参数的二元线性模型，模型决定系数R2分别为0.928、0.956、0.949、0.949、0.889和0.815。
VMD、VAD、S/V、SMD、NAD及NMD均随压力呈二次多项式变化规律，建立了这6个雾滴参数的二次多项式模型，决定系数R2均达到了0.894以上。
不同测试点雾滴粒径参数的变化规律有待进一步研究；其他因素对雾滴粒径的影响有待进一步研究。
参考文献
[1] Fife J P，Ozkan H E，Derksen R C.Viability of a biological pest control agent through hydraulic nozzles[J]. Transactions of the ASAE，2005，48（1）：4554.
[2] Dorr G J，Hewitt A J，Adkins S W，et al. A comparison of initial spray characteristics produced by agricultural nozzles[J]. Crop Protection，2013，53（11）：109117.
[3] 吕晓兰，傅锡敏，宋坚利，等. 喷雾技术参数对雾滴飘移特性的影响[J]. 农业机械学报，2011，42（1）：5963.
[4] 袁会珠，王忠群，孙瑞红，等. 喷洒部件及喷雾助剂对担架式喷雾机在桃园喷雾中的雾滴沉积分布的影响[J]. 植物保护，2010，36（1）：106109.
[5] 贾卫东，胡化超，陈龙，等. 风幕式静电喷杆喷雾喷头雾化与雾滴沉积性能试验[J]. 农业工程学报，2015，31（7）：5359.
[6] 吕晓兰，傅锡敏，吴萍，等. 喷雾技术参数对雾滴沉积分布影响试验[J]. 农业机械学报，2011，42（6）：7075.
[7] 杨希娃，代美灵，宋坚利，等. 雾滴大小、叶片表面特性与倾角对农药沉积量的影响[J]. 农业工程学报，2012，28（3）：7073.
[8] 徐德进，顾中言，徐广春，等. 雾滴密度及大小对氯虫苯甲酰胺防治稻纵卷叶螟效果的影响[J].中国农业科学，2012，45（4）：666674.
[9] 袁会珠，杨代斌，闫晓静，等.农药有效利用率与喷雾技术优化[J]. 植物保护，2011，37（5）：1420.
[10]陈丹，任广伟，王秀芳，等. 4种喷雾器在茶树上喷雾效果比较[J]. 植物保护，2011，37（5）：110114.
[11]李飞，王相晶，吴青君，等. 三种药剂喷雾和灌根施药方式对西花蓟马的残留毒力[J]. 植物保护，2013，39（3）：173177.
[12]王双双，何雄奎，宋坚利，等. 农用喷头雾化粒径测试方法比较及分布函数拟合[J]. 农业工程学报，2014，30（20）：3442.
[13]Krishnan P，Evans T，Ballal K.Scanning electron microscopic studies of new and used fan nozzles for agricultural sprayers [J]. Applied Engineering in Agriculture，2004，20（2）：133137.
[14]毛益进，王秀，马伟. 农药喷洒雾滴粒径分布数值分析方法[J]. 农业工程学报，2009，25（S2）：7882.
[15]史春建，邱白晶，汤伯敏，等. 基于高速图像的雾滴尺寸分布统计与运动分析[J]. 农业机械学报，2006，37（5）：6366.
[16]Santangelo P E.Characterization of highpressure watermist sprays： experimental analysis of droplet size and dispersion [J]. Experimental Thermal and Fluid Science，2010，34（8）：13531366.
[17]胡桂琴，许林云，周宏平，等.影响空心圆锥雾喷头雾滴粒径的多因素分析[J].南京林业大学学报：自然科学版，2014，38（2）：133136.
[18]Santangelo P E.Characterization of highpressure watermist sprays： experimental analysis of droplet size and dispersion [J]. Experimental Thermal and Fluid Science，2010，34（8）：13531366.
[19]何雄奎. 改变我国植保机械和施药技术严重落后的现状[J]. 农业工程学报，2004，20（1）：1315.
[20]洪添胜，杨洲，宋淑然，等. 柑橘生产机械化研究[J]. 农业机械学报，2010，41（12）：105110.
[21]代秋芳，洪添胜，宋淑然，等. 基于PSTN和PLMN的果园泵房智能监控系统[J]. 农业机械学报，2012，43（4）：162167.
[22]宋淑然，阮耀灿，洪添胜，等. 果园管道喷雾系统药液压力的自整定模糊PID控制[J]. 农业工程学报，2011，27（6）：157161.
（责任编辑：田 喆）
关键词 山地果园； 管道恒压喷雾； 雾滴粒径参数； 喷雾压力； 孔径
中图分类号： S 491
文献标识码： A
雾滴粒径参数是反映雾滴大小的主要指标[12]，直接影响雾滴的飘移[3]、沉积[47]和药效的发挥[8]，从而影响农药的利用率[9]和喷雾效果[1011]。国内外学者采用各种方法对雾滴粒径进行测量[12]，用图像处理技术对雾滴进行显微识别[13]，利用高速摄像机结合数字图像处理技术对雾滴粒径分布和雾滴的运动进行分析[1415]，利用激光粒度仪测量雾滴谱[16]，同时对影响喷头雾滴粒径参数的各种因素进行了分析[1718]。但这些研究多为测量在较低压力喷雾时喷头的雾滴，大多集中在单因素对雾滴粒径的影响分析，而对适合果园需求的高压力喷雾时喷头的雾滴大小研究较少，主要影响因素与雾滴大小之间的相关关系研究也较少。
我国南方果园多为山地果园，多采用低矮密植型种植模式，移动式喷雾机械很难进入[1920]，为此笔者所在团队研究了果园管道恒压喷雾控制系统。由于管道喷雾无需机具在果园中移动，具有省力、省工、喷雾效果好等优点，适合山地果园使用[21]。管道喷雾依靠液泵形成的较高压力，将药液在喷头处进行雾化并施向靶标，但管道内药液的压力随作业人数的不同会产生波动，同时，药液在管道中流动时会出现压力损失，导致整个管网系统药液压力不稳定且各点压力各不相同。管道恒压控制装置可根据管道中压力的实际值与设定值间的误差及误差变化趋势，在线调整模糊PID的参数，经带有变速积分、微分先行优化算子的增量式PID算法计算，获得控制量以控制管道中药液的压力[22]。因此压力成为影响管道喷雾系统中喷头雾滴参数的最重要的因素之一。除压力外，在喷头类型选定后，喷头孔径也是影响雾滴粒径参数的重要因素。然而到目前为止，还未见有文献对管道恒压喷雾系统中应如何设置喷雾压力及如何选择喷头的孔径进行研究。
本文通过改变压力和孔径，利用激光粒度仪采集雾滴粒径参数数据，研究压力和孔径对果园用空心圆锥雾喷头雾滴大小的影响，掌握雾滴大小随压力和孔径的变化规律，可以根据喷施对象的所需雾滴大小设置管道恒压喷雾控制装置的压力，选择合适的喷头类型，从而提高喷雾效果。
1 描述雾滴粒径的参数
（1）VMD：体积中值粒径（volume median diameter，VMD）。指取样雾滴的体积按雾滴大小顺序进行累计，其累计值为取样雾滴体积总和的50%所对应的雾滴直径。根据雾滴群VMD的大小将雾滴分为4类：①粗雾：VMD≥400 μm；②细雾：100 μm≤VMD<400 μm；③弥雾：50 μm≤VMD<100 μm；④气溶胶：VMD<50 μm。
（2）VAD：体积平均粒径（volume average diameter，VAD）。指取样雾滴平均体积所对应的直径。
（3）S/V：体积比表面积。指单位体积雾滴的表面积。
（4）SMD：表面积平均粒径（superficial area median diameter，SMD）。指粒径对表面积的加权平均。
（5）NAD：数量平均粒径（number average diameter，NAD）。指取样雾滴群的直径之和与雾滴群个数之和的比值。
（6）NMD：数量中值粒径（number median diameter，NMD）。指取样雾滴的个数按雾滴大小顺序进行累计，其累计值为取样雾滴个数总和的50%所对应的雾滴直径。
2 数据获取
雾滴粒径参数数据通过激光粒度仪获取。激光粒度仪通过测量颗粒群的散射谱，散射谱的强度与被测颗粒群的大小有关，并经透镜再次汇聚后，被位于透镜后焦面上的光电阵列探测器所接收，转换成电信号后经放大和A/D转换经通讯口送入计算机，进行反演运算和数据处理后，即可给出被测颗粒群的雾滴粒径参数。
2.1 试验设备
2.1.1 喷雾性能综合试验平台
喷雾性能综合试验平台主要由机械结构、电气控制和数据采集处理3部分组成。可提供本试验中所需雾滴并测量雾滴粒径参数。
2.1.2 激光粒度仪
Winner318B激光粒度仪可测定雾滴粒径参数。该仪器采用信息光学原理，通过测量颗粒群的散射谱，来分析雾滴粒径大小。其准确性误差<3%；重复性误差<3%；绿色激光器，波长532 nm，功率30 mW；量程为15～711 μm。
2.1.3 拉丝模孔径测量仪
拉丝模孔径测量仪可测量喷头孔径。测量仪将光学与数码、软件相结合，在软件中对孔径进行测量，可得出圆孔的直径。测量仪的分辨率为0.1 μm；重复精度为±0.5 μm；测量范围为0.05～30 mm；变倍方式为连续变倍，最大倍率为800倍。 2.2 试验对象
对Lechler TR80系列3种型号TR8001C、TR8002C和TR8003C喷头的雾滴粒径参数进行测量，喷头内部结构图如图1所示。该系列喷头压力范围为0.3～2.0 MPa，是一种带陶瓷芯的空心圆锥雾喷头，高耐磨，耐化学药剂，有极细雾滴谱，适用于压力要求较高的山地果园管道喷雾系统。
该系统在果园埋下塑料管网，管网进液端配有药池、恒压喷药控制系统，埋设的管道每隔一段距离就安装上升至地面的立式管，每个立管在打药时再与软管、喷枪相连接，使用电动机带动药泵加压，把药液输送到果树旁，人工逐株打药，适合多人同时打药，如图2所示。恒压喷药控制系统由变频器、变频电动机、药泵、压力变送器、恒压控制箱组成，可使出水口压力保持恒定，出水口压力可根据本文研究成果进行合理设定。
2.3 试验方法
试验用水为清水；用喷雾性能综合试验平台调节压力在0.70～1.40 MPa间变化，每次步进0.10 MPa；用拉丝模孔径测量仪对喷头孔径进行测量，测得TR8001C的喷嘴孔径为1.0 mm，TR8002C为1.4 mm，TR8003C为1.8 mm；用激光粒度仪对喷头中轴线上距喷嘴15 cm处的雾滴颗粒群的散射谱进行测量，得到雾滴能谱数据，再将能谱数据转换为雾滴粒径参数数据。
试验前用标准物质（直径为25 μm）对激光粒度仪进行标定。每次测量前先进行背景测试，背景测试累计10次后再进行样品测试以去除背景光的影响。每条样品测试数据为某个时刻样品的雾滴参数。本试验在每种工况下对样品连续测量100次得到100条测试数据，100条测试数据取平均值后作为每种工况下样品的测试数据。
3 数据处理与分析
首先将试验获得的能谱数据转换为雾滴粒径数据，然后分析各粒径级雾滴的分布情况，分析雾滴粒径大小，给出雾滴参数随压力和孔径的变化趋势，并对雾滴参数进行多元线性回归，最后建立基于压力的雾滴参数模型。
能谱数据向雾滴粒径数据转换时产生的计算误差称拟合误差。本试验拟合误差均小于0.012，符合试验要求。
本试验利用Origin 9.0进行绘图及曲线拟合，利用SPSS 16.0对试验数据进行多元线性回归分析，利用Excel 2007进行数据处理。
3.1 各粒径级雾滴的分布情况
根据3种喷头8种压力下雾滴的粒径分布数据，绘制各粒径级雾滴百分比的分布如图3所示。由于粒径D≥100 μm的雾滴很少，均少于0.3%，可忽略不计，所以仅绘制0～100 μm范围内雾滴的分布情况。可以看出：
（1） 雾滴均较细小。粒径为0～20 μm的雾滴比例均少于0.5%，可以忽略不计；主要是粒径为20～40 μm的雾滴，其比例均多于79.5%；粒径为40～60 μm的雾滴较少，其比例均少于10.6%；粒径为60～80 μm的雾滴较少，其比例均少于7.4%；粒径为80～100 μm的雾滴较少，其比例均少于2.1%。
（2） 按雾滴所占百分比从大到小排列为20～40 μm> 40～60 μm>60～80 μm>80～100 μm。
（3） 压力越大，孔径越小，雾滴粒径分布越集中，雾滴越细小。
3.2 雾滴粒径参数随压力和孔径的变化
为详细了解雾滴粒径参数的变化趋势，分别绘制出6种雾滴粒径参数VMD、VAD、S/V、SMD、NAD和NMD随压力和孔径变化的曲线，如图4所示。
从图4可以看出：
（1） 压力在0.70～1.40 MPa变化时，TR8001C、TR8002C和TR8003C的VMD变化范围分别为30.6～31.9 μm、30.9～32.0 μm和31.0～32.317 μm。说明在试验压力范围内，雾滴均较细小且全部为气溶胶。
（2） 孔径不变、压力越大，VMD、VAD、SMD、NAD和NMD越小，S/V越大。说明压力越大，雾滴的粒径越小。原因是压力越大，液体流出喷孔的初速度就越大，喷孔中液体扰动程度及喷出喷孔后所受到的空气介质阻尼也越大，雾滴更易破碎。
（3） 压力不变，孔径越小，VMD、VAD、SMD、NAD和NMD越小，S/V越大，即雾滴的粒径越小；原因是喷孔直径越小，液体喷出喷孔时会受到更大的节流，喷孔内扰动增加，雾滴更易破碎。
（4） VMD和VAD随压力和孔径的变化趋势较NAD和NMD更明显。说明：侧重反映体积的粒径参数比侧重反映数量的粒径参数随压力和孔径的变化趋势更明显。
3.3 基于压力和孔径的多元回归统计模型
雾滴的粒径大小与压力和孔径密切相关。从图4可以看出，6种雾滴粒径参数的变化趋势与压力和孔径的变化趋势基本一致，具有很强的相关性。其中，压力取值范围为0.80～1.40 MPa，步进0.10 MPa；孔径取值范围为1.0～1.8 mm，步进0.4 mm。分析结果如表1～3所示，可以看出：
（1）每种雾滴参数多元回归统计模型如表1所示，其中：y为各种雾滴粒径参数，x1为喷头孔径（单位为mm），x2为喷雾压力（单位为MPa）。模型的复相关系数R、决定系数R2、调整决定系数R2C均较高，其中VMD、VAD、S/V、SMD和NAD的决定系数R2分别为0.928、0.956、0.949、0.949、0.889和0.815。说明回归方程拟合优度较高，符合拟合要求。
（2） 由F检验可知，所有雾滴参数均有显著值P=0.000<0.01，表明VMD、VAD、S/V、SMD、NAD及NMD与孔径和压力之间均存在极显著的线性回归关系（见表2）。
（3） 由t检验可知，VMD、VAD、S/V、SMD及NAD的回归系数b0、b1和b2均有显著值P=0.000<0.01，表明方程的回归系数b0项（常数项）、b1项（孔径系数）和b2项（压力系数）影响均极显著；NMD的回归系数b0、b2均有显著值P=0.000<0.01，b1有显著值P=0.032<0.05，表明方程的回归系数b0项和b2项影响均极显著，回归系数b1项影响显著（见表3）。 综上，回归通过了F检验和t检验，且拟合优度较高，说明本文建立的雾滴粒径参数与压力和孔径的多元回归模型是准确的。
3.4 雾滴粒径参数与压力的关系
经曲线回归依次分析VMD、VAD、S/V、SMD、NAD及NMD随压力的变化规律。在进行多项式回归时，同时利用线性函数、指数函数、对数函数和幂函数进行了回归分析，这些函数的相关系数均比多项式回归的相关系数要小，故采用多项式表示各雾滴参数随压力的变化规律（见表4）。以喷雾压力（单位MPa）为横坐标，雾滴参数为纵坐标，得3种型号喷头的6种雾滴参数与喷雾压力（0.80～1.40 MPa）的拟合曲线如图5所示。从表4和图5可以看出，雾滴参数随压力的变化规律均能用二次多项式表示，且回归均具有较高的决定系数，R2均达到了0.894以上。
4 结论
Lechler TR8001C、TR8002C和TR8003C型喷头在0.70～1.40 MPa压力下，粒径D≥100 μm的雾滴均少于0.3%，粒径为0～20 μm的雾滴均少于0.5%，可忽略不计，按雾滴所占百分比从大到小排列20～40 μm>40～60 μm>60～80 μm>80～100 μm；主要是粒径为20～40 μm的雾滴（79.5%～92.8%）；VMD为30.6～32.3 μm，VAD为31.3～36.9 μm，雾滴均较细小，全部为气溶胶；压力越大、孔径越小，雾滴粒径越小。
VMD、VAD、S/V、SMD、NAD及NMD与孔径和压力均有较好的二元线性关系，建立了这6个雾滴参数的二元线性模型，模型决定系数R2分别为0.928、0.956、0.949、0.949、0.889和0.815。
VMD、VAD、S/V、SMD、NAD及NMD均随压力呈二次多项式变化规律，建立了这6个雾滴参数的二次多项式模型，决定系数R2均达到了0.894以上。
不同测试点雾滴粒径参数的变化规律有待进一步研究；其他因素对雾滴粒径的影响有待进一步研究。
参考文献
[1] Fife J P，Ozkan H E，Derksen R C.Viability of a biological pest control agent through hydraulic nozzles[J]. Transactions of the ASAE，2005，48（1）：4554.
[2] Dorr G J，Hewitt A J，Adkins S W，et al. A comparison of initial spray characteristics produced by agricultural nozzles[J]. Crop Protection，2013，53（11）：109117.
[3] 吕晓兰，傅锡敏，宋坚利，等. 喷雾技术参数对雾滴飘移特性的影响[J]. 农业机械学报，2011，42（1）：5963.
[4] 袁会珠，王忠群，孙瑞红，等. 喷洒部件及喷雾助剂对担架式喷雾机在桃园喷雾中的雾滴沉积分布的影响[J]. 植物保护，2010，36（1）：106109.
[5] 贾卫东，胡化超，陈龙，等. 风幕式静电喷杆喷雾喷头雾化与雾滴沉积性能试验[J]. 农业工程学报，2015，31（7）：5359.
[6] 吕晓兰，傅锡敏，吴萍，等. 喷雾技术参数对雾滴沉积分布影响试验[J]. 农业机械学报，2011，42（6）：7075.
[7] 杨希娃，代美灵，宋坚利，等. 雾滴大小、叶片表面特性与倾角对农药沉积量的影响[J]. 农业工程学报，2012，28（3）：7073.
[8] 徐德进，顾中言，徐广春，等. 雾滴密度及大小对氯虫苯甲酰胺防治稻纵卷叶螟效果的影响[J].中国农业科学，2012，45（4）：666674.
[9] 袁会珠，杨代斌，闫晓静，等.农药有效利用率与喷雾技术优化[J]. 植物保护，2011，37（5）：1420.
[10]陈丹，任广伟，王秀芳，等. 4种喷雾器在茶树上喷雾效果比较[J]. 植物保护，2011，37（5）：110114.
[11]李飞，王相晶，吴青君，等. 三种药剂喷雾和灌根施药方式对西花蓟马的残留毒力[J]. 植物保护，2013，39（3）：173177.
[12]王双双，何雄奎，宋坚利，等. 农用喷头雾化粒径测试方法比较及分布函数拟合[J]. 农业工程学报，2014，30（20）：3442.
[13]Krishnan P，Evans T，Ballal K.Scanning electron microscopic studies of new and used fan nozzles for agricultural sprayers [J]. Applied Engineering in Agriculture，2004，20（2）：133137.
[14]毛益进，王秀，马伟. 农药喷洒雾滴粒径分布数值分析方法[J]. 农业工程学报，2009，25（S2）：7882.
[15]史春建，邱白晶，汤伯敏，等. 基于高速图像的雾滴尺寸分布统计与运动分析[J]. 农业机械学报，2006，37（5）：6366.
[16]Santangelo P E.Characterization of highpressure watermist sprays： experimental analysis of droplet size and dispersion [J]. Experimental Thermal and Fluid Science，2010，34（8）：13531366.
[17]胡桂琴，许林云，周宏平，等.影响空心圆锥雾喷头雾滴粒径的多因素分析[J].南京林业大学学报：自然科学版，2014，38（2）：133136.
[18]Santangelo P E.Characterization of highpressure watermist sprays： experimental analysis of droplet size and dispersion [J]. Experimental Thermal and Fluid Science，2010，34（8）：13531366.
[19]何雄奎. 改变我国植保机械和施药技术严重落后的现状[J]. 农业工程学报，2004，20（1）：1315.
[20]洪添胜，杨洲，宋淑然，等. 柑橘生产机械化研究[J]. 农业机械学报，2010，41（12）：105110.
[21]代秋芳，洪添胜，宋淑然，等. 基于PSTN和PLMN的果园泵房智能监控系统[J]. 农业机械学报，2012，43（4）：162167.
[22]宋淑然，阮耀灿，洪添胜，等. 果园管道喷雾系统药液压力的自整定模糊PID控制[J]. 农业工程学报，2011，27（6）：157161.
（责任编辑：田 喆）