2. 国家柑橘产业技术体系机械化研究室,广东 广州 510642;
3. 广东省农情信息监测工程技术研究中心,广东 广州 510642;
4. 广东省山地果园机械创新工程技术研究中心,广东 广州 510642
2. Division of Citrus Machinery, China Agriculture Research System, Guangzhou 510642, China;
3. Guangdong Engineering Research Center for Moinitoring Agricultural Information, Guangzhou 510642, China;
4. Guangdong Engineering Technology Research Center for Creative Mountainous Orchard Machinery, Guangzhou 510642, China
近年来,随着我国农业生产的快速发展,对农业机械化的需求大幅增加,推进农机装备技术改进,研发适合国情、农民需要、先进适用的各类智能农机装备迫在眉睫。在农机装备智能领域,高效喷雾施药技术的一个主要标准是农药在靶标上的沉积率远高于其在非目标物或区域中的沉积率。这也是最大程度地提高农药有效利用率,同时消除农药对操作人员的危害以及降低环境污染的一个重要途径[1-3]。然而,我国在农业生产中对农药的有效利用率只有30%左右,流失量较高[4-5]。
农药雾化是农药以雾滴的形式分散到大气中,形成雾状分散体系的过程,其实质是喷雾液滴在外力的作用下实现比表面积的大幅度扩增[6-8]。将农药药液雾化为雾滴的过程依靠喷头来实现[9]。农药雾化程度直接影响农药漂移距离及沉积有效利用率。喷雾雾滴粒径分布是农药雾化程度的主要检测指标[10-11]。较大雾滴能够在较长时间内保持动量,到达靶标时间短,漂移较小,但雾滴过大又容易造成药液覆盖率降低、靶标附着性差和药液流失[12-14]。较小的雾滴由于其质量小,受到空气阻力影响大,通常会由于动量不足而无法到达靶标[15],但小雾滴有助于增大药液覆盖率、提升雾滴覆盖的均匀性和增大药液在作物冠层的穿透性[16]。
无论采用何种喷雾方式,植保机械的作业质量与雾滴直径大小都有直接关系[17]。如果选择的雾滴大小合适,可以用最小的药量,以最小的环境污染达到最大控制病虫害的目的。如果实际的雾滴比需要的雾滴大,所浪费的农药就会以雾滴直径三次方的速率增长[18]。张瑞瑞等[19]利用不同类型及浓度助剂分析比较空气诱导喷头IDK120-025和多量程平面喷头LU120-015的雾滴体积中径(Volume median diameter, VMD)及雾滴粒径相对分布跨度,结果表明不同类型及浓度的助剂对不同喷头效果不同,应根据所用喷头选择合适的助剂及配比浓度,降低农药雾滴漂移。王潇楠等[20]利用不同浓度抗蒸发助剂Agrospred 730、防漂移助剂Break-thru Vibrant、Silwet DRS-60、Greenwet 360分析比较了不同喷头雾滴漂移潜在指数,结果表明添加助剂浓度不同,雾滴体积中径不同,喷头潜在漂移指数不同。何玲等[21]利用不同表面张力的溶液分析了喷雾助剂对水稻冠层有效沉积率的影响。唐青等[22]利用不同气流条件对比分析了标准扇形喷头和空气诱导喷头的雾化特性。赵辉等[23]研究了不同喷雾液动态表面张力下,雾滴粒径的变化,结果表明动态表面张力改变对雾滴粒径的影响增大。
液体黏度影响喷头雾化性能[24],本文配制不同黏度的喷雾液,研究喷雾液黏度对农用喷头的喷雾雾滴粒径的影响,以期为农药喷施技术及机械性能的提高提供理论支撑与试验依据。
1 材料与方法 1.1 试验材料本文选择乳油、可湿性粉剂、悬浮剂和水分散粒剂4种常见剂型的农药作为研究对象,具体为:体积分数为15%的哒螨灵乳油、质量分数为10%的吡虫啉可湿性粉剂、质量浓度为240 g/L的螺螨酯悬浮剂和质量分数为70%的啶虫脒水分散粒剂。按照农药剂型的使用说明把4种农药分别配制成溶液,利用黏度计测试4种药液的黏度,获得不同剂型农药在使用配比下药液的黏度范围。据此,用水与甘油配制喷雾液代替农药进行大量的喷雾试验研究,避免用农药进行试验造成环境污染。同时利用甘油与水配制出不同黏度的喷雾液进行对比试验。
试验喷头为东莞市沙鸥喷雾系统有限公司生产的实心锥喷头,喷嘴直径为0.8 mm。
1.2 试验仪器试验所用的主要仪器有:丹麦DANTEC公司生产的多普勒动态粒子分析仪(Phase Doppler anemometry,简称PDA),粒子测量范围为1~10 000 μm,测量误差为1%;上海方瑞仪器有限公司生产的NDJ-8T数字旋转黏度计,测量范围为1~2 000 cp,测量误差为±1%,重复误差为±0.5%,仪器均满足相关参数测量要求。
试验环境温度对液体黏度的测试有一定程度的影响,在液体黏度测试过程中,环境温度是动态变化的。为此,本文在每次试验前,均多次测量环境温度,并使用NDJ-8T型数字旋转黏度计多次对同一个样品进行黏度测试,保证测量结果的稳定性和可信性。
1.3 雾滴参数定义本文主要对雾滴的3个指标参数进行了测量和分析,它们分别是:
算术平均直径(Arithmetic mean diameter,D10):指取样雾滴群的直径之和与雾滴群个数之和的比值,
${D_{10}} = \frac{1}{N}\sum\limits_{i = 1}^{{N_i}} {{n_i}{D_i}}; $ | (1) |
体积平均直径(Volume mean diameter,D30):指取样雾滴群平均体积所对应的直径,
${D_{{\rm{30}}}}{\rm{ = }}{\left[ {\frac{{\rm{1}}}{N}\sum\limits_{i = 1}^{{N_i}} {{n_i}D_i^3} } \right]^{\frac{1}{3}}};$ | (2) |
索尔特平均直径(Sauter mean diameter,D32):表明雾滴粒径对表面积的加权平均,
${D_{32}} = \frac{{\displaystyle\sum\limits_{i = 1}^{{N_i}} {{n_i}D_i^3} }}{{\displaystyle\sum\limits_{i = 1}^{{N_i}} {{n_i}D_i^2} }};$ | (3) |
在上述雾滴参数描述中,Di是指雾滴尺寸等级i的直径,Ni是指选择的雾滴尺寸等级i的个数,ni是指每个尺寸等级中的雾滴数,N是雾滴总数。
1.4 试验设计液体的黏度和温度有极密切的关系,通常温度升高,黏度降低。本文测试液体黏度是在环境温度为(25±0.5) ℃的前提下进行的。
1.4.1 4种常用于防治木虱的农药黏度范围测试为了揭示农药黏度对喷头雾滴参数的影响,根据农药使用说明的配比要求,配制了4种柑橘病虫害防治中常用的农药溶液,利用黏度计分别测量其黏度值,每隔2 min读取1次黏度值,共读取3次,取平均值作为最终黏度值,结果如表1所示。4种农药溶液黏度的平均最小值为1.25 cp,平均最大值为1.80 cp。
![]() |
表 1 4种常见农药的黏度 Table 1 Viscosities of four common pesticides |
由于甘油与水能够以任意比例互溶,且纯甘油黏度较大,本文采用纯甘油与水以一定的比例混合,配制4种不同质量分数的甘油溶液,结果如表2所示。随着甘油质量分数的增大,溶液的黏度也逐渐变大,黏度范围能够涵盖常用的农药溶液黏度。
![]() |
表 2 不同质量分数甘油溶液的黏度 Table 2 Viscosities of different mass fractions of glycerin solution |
喷雾液黏度对喷头雾化性能的影响试验包括:不同喷雾液黏度下对喷头轴心线方向上雾滴参数的测定,以及不同喷雾液黏度下对喷头径向(轴心方向高度一定时)雾滴参数的测定。喷雾试验时,喷头轴心方向垂直地面,喷雾压力为0.9 MPa。
前期试验中,为了能够设置合适的喷头喷雾区域采样点,在喷雾压力为1.0 MPa条件下,测量了喷头的雾锥角,结果为44.14°,由此,利用预试验测得的雾锥角正切值的结果,可以确定喷雾区域的边界位置,从而确定雾滴采样点的范围。
轴心方向上雾滴采样点设定:为研究不同黏度对喷头轴心方向上雾滴粒径的影响,建立以喷头处为原点、水平向左方向为x、垂直向里方向为y、垂直向下方向为z的笛卡尔坐标系,在喷头雾化区域轴心方向上,从喷头z=0 cm处开始,每隔2 cm利用PDA分别测量轴心方向上雾滴参数。
径向雾滴采样点设定:为研究黏度对喷头径向雾滴粒径的影响,在轴心方向上分别选取z=30 cm,z=40 cm处的水平面,测量所对应径向上的雾滴粒径;在z=30 cm处x选取−8~8 cm之间;在z=40 cm处,x选取−10~10 cm之间,其中径向(x方向)采样点间隔均为1 cm,利用PDA分别测量每个采样点上的雾滴参数。
1.4.4 数据处理运用统计分析软件SPSS16.0对轴心距离、液体黏度以及雾滴粒径的相关性进行统计分析。
2 结果与分析 2.1 轴心方向上雾滴粒径分布特性从喷头开始,每隔2 cm,测量喷头轴心方向上的雾滴参数,结果如表3所示。
![]() |
表 3 不同黏度下轴心方向上雾滴参数1) Table 3 Droplet parameters in the axis direction under different viscosities |
从表3可以看出,在轴心方向上距离喷头1~6 cm的区间内,雾滴粒径变小趋势明显,分析原因可能是:在一定的压力下,射流以连续液体的形式从喷头高速流出,受到外界气体的扰动作用,液体表面会产生不稳定的波动,当液体流速较大时,液体表面的波幅将快速变大,最终波峰被撕裂下来,形成较大雾滴。随着射出距离变大,这种表面波的波幅将逐渐扩大,使液体分裂成大颗粒液滴,液滴直径超过了一定的阈值,在周围空气切应力的作用下,它们将进一步碎裂成更细小的液滴,故在距离喷头1~6 cm的区间内呈现出变小的趋势。
Rayleigh理论导出了形成黏性射流最大不稳定性的比值,即有一个最小的扰动波长
${\lambda _{\min }} = {\text{π}} d,$ | (4) |
${\lambda _{{\rm{opt}}}} = \sqrt 2{\text{π}} d,$ | (5) |
对于黏性流体
${\lambda _{\min }} = {\text{π}} d,$ | (6) |
${\lambda _{{\rm{opt}}}} = \sqrt 2 {\text{π}} d{\left( {1 + \frac{{3{u_{_{\rm{l}}}}}}{{\sqrt {{\rho _{_{\rm{l}}}}{\sigma _{_{\rm{l}}}}d} }}} \right)^{{1 / 2}}},$ | (7) |
式中:
试验表明,在相同的喷雾条件下,对于不同黏度的喷雾液,液体密度
在轴心方向上距离喷头8~40 cm的区间内,不同轴心距离下的雾滴粒径呈现出从近喷头处到远喷头处逐渐变大的规律,其原因可能是:在喷头的轴心方向上,作用于雾滴的力主要是重力和空气阻力,当雾滴受到重力影响时,会加速下降,雾滴发生聚合碰撞,使雾滴粒径逐渐变大,形成了大雾滴。
2.2 轴心雾滴粒径的统计分析为了明确轴心距离、液体黏度与雾滴粒径之间的内在联系,更好地揭示影响雾滴粒径分布规律的主要原因,本文运用统计分析软件SPSS对轴心距离、液体黏度以及雾滴粒径进行相关性统计分析。
通过判断2个数据集是否在一条直线上来推断2个数据集的线性相关性,其计算公式为:
$r = \frac{{N\displaystyle\sum {{x_i}{y_i}} - \displaystyle\sum {{x_i}\displaystyle\sum {{y_i}} } }}{{\sqrt {N\displaystyle\sum {x_i^2 - {{\left(\displaystyle\sum {{x_i}} \right)}^2}} } \sqrt {N\displaystyle\sum {y_i^2 - {{\left(\displaystyle\sum {{y_i}} \right)}^2}} } }},$ | (8) |
式中,N为2个数据集样本个数,xi表示x数据集中某一样本(i)的样本值,yi表示y数据集中样本i的样本值,r值反映了2个变量线性相关性的强弱程度,r值越大说明两者相关性越强,据此,利用软件做出了液体黏度、轴心距离与雾滴粒径三者之间的Pearson相关系数(表4)。
![]() |
表 4 轴心雾滴粒径的Pearson相关分析 Table 4 Pearson correlation analysis of axial droplet size |
从表4可以看出,当轴心距离为4~40 cm时,也就是采用样本为4~40 cm区间内雾滴粒径数据,其轴心距离、液体黏度与雾滴粒径的Pearson相关系数(r)分别为0.531、0.795,可见液体黏度与雾滴粒径之间的相关性更强,接近于极强相关,显著性值较高(当显著性值小于0.05时认为显著),据此可得出,在距离喷嘴4~40 cm区间内,雾滴粒径的大小受到液体黏度的影响较大,同时轴心距离也是影响雾滴粒径变化的主要因素;当轴心距离为0~40 cm时,其轴心距离、液体黏度与雾滴粒径的Pearson相关系数分别为0.287、0.586,可见当计算范围包含0~4 cm区间时,液体黏度与雾滴粒径大小的相关性为中等程度相关,而轴心距离与雾滴粒径大小的相关性为弱相关,这说明雾滴粒径在区间0~4 cm与区间4~40 cm内变化趋势不统一,故而三者之间的相关性有所降低,结合表3不难发现,在雾滴粒径为0~4 cm区间内,喷雾液由喷嘴喷出后,在切应力的作用下撕裂成液膜,大雾滴变为小雾滴,这刚好与雾滴在4~40 cm区间内的变化相反,因此相关性有所降低。
综合考虑可以发现,雾滴从喷嘴喷出后液体黏度对雾滴粒径的影响大于轴心距离对雾滴粒径的影响。
2.3 径向雾滴粒径分布为了了解不同喷雾液黏度对喷头径向雾滴参数的影响(为便于分析,取z=30 cm处分析),绘制出如图1所示的雾滴体积平均直径D30随径向距离变化三维图。
![]() |
图 1 z=30 cm处径向雾滴体积平均直径(D30)变化
Fig. 1 Change of radial droplet D30 at
|
由图1可以看出,在相同喷雾液黏度下,雾滴粒径随着径向距离增大而增大,轴心位置始终保持着最小的雾滴粒径,径向左右两边雾滴粒径近似于对称分布;不同喷雾液黏度下,液体黏度越大,相同位置雾滴体积平均直径越大。
奥内佐格的研究表明,射流破碎后的液滴尺寸主要取决于喷口直径、液体的密度、黏性力和表面张力,并且整理出了无因次准数Z来进行描述,Z为射流破碎后的雾滴大小的无因次准数,它与液体的密度、黏性力和表面张力间存在一定的关系[25],如下式所示:
$Z = \frac{{{u_{_{\rm{l}}}}}}{{\sqrt {{\rho _{\rm l}}{\sigma _{\rm l}}{d_0}} }},$ | (9) |
式中:
在相同试验条件下,对于不同黏度的喷雾液,其液体密度、表面张力、射流直径变化不大,故液体的黏度是雾化过程中液体破碎的主要因素,黏性在射流表面扰动波增长的过程中起阻尼作用,黏性越小阻力越小,Z越小,故雾滴粒径越小。
2.4 径向雾滴粒径的统计分析为了更好地阐释雾滴粒径在径向上的分布规律,明确雾滴粒径在径向上与径向距离、液体黏度之间的内在制约因素,对径向距离、雾滴粒径以及液体黏度之间进行了Pearson相关分析,结果如表5所示。
![]() |
表 5 径向雾滴粒径的Pearson相关分析1) Table 5 Pearson correlation analysis of radial droplet size |
由表5可以看出,在径向上,雾滴粒径与径向距离的相关系数为0.932,具有极强相关性,且显著性值为0.000,小于0.05,故认为雾滴粒径与径向距离之间显著性高;雾滴粒径与液体黏度的相关系数为0.328,表现为弱相关,且显著性值为0.014,小于0.05,则认为雾滴粒径与液体黏度之间联系较径向距离弱。
由表4和表5可以发现,在喷头喷雾区域中,在轴心方向上雾滴粒径分布受液体黏度影响较大,二者相关系数为0.795,同时,在轴心方向上,雾滴粒径也会受到轴心距离的影响,相关系数为0.531,较液体黏度影响相对较弱;在喷头径向上,雾滴粒径与径向距离相关系数为0.932,表现出极强相关性,意味着径向位置不同时,雾滴粒径变化明显,由图1可以看出,当径向位置变大时,雾滴粒径增加明显;在径向上,雾滴粒径与液体黏度相关系数为0.328,表现为弱相关性。
3 结论本文研究了不同喷雾液黏度下,喷头的喷雾区域中垂直(轴心)与水平(径向)方向上喷雾雾滴粒径变化规律,得出以下结论:
1)喷雾液黏度越大,喷头所喷出雾滴的粒径也越大。随着黏度的增大,轴心方向上的雾滴粒径先变大后变小,喷雾液黏度越大,雾滴越难以雾化;在径向上,喷雾区域的边界处雾滴粒径大于内部雾滴粒径;
2)根据Pearson相关分析得出,轴心方向上,雾滴粒径的大小分布受液体黏度影响较大,相关系数为0.795;而在径向上,雾滴粒径大小分布受到径向距离影响较大,相关系数为0.932;
3)试验分析得出添加不同浓度的甘油溶液均可对雾滴粒径产生较大影响,因此应根据喷施的实际情况选择合适的浓度配比,适当提高喷雾液黏度,有助于达到降低农药雾滴漂移、增加雾滴穿透性的目的。
[1] |
许晏铭, 刘建河, 李星光. 液体物性对静电喷雾雾化性能的研究[J]. 排灌机械工程学报, 2018, 45(7): 793-798. ( ![]() |
[2] |
MANSOUR A, CHIGIER N. Air-blast atomization of non-Newtonian liquids[J]. Journal of Non-Newtonian Fluid Mechanics, 1995, 58(2): 161-194. ( ![]() |
[3] |
DAI X, XU Y, ZHENG J, et al. Analysis of the variability of pesticide concentration downstream of inline mixers for direct nozzle injection systems[J]. Biosystems Engineering, 2019, 180(13): 59-69. ( ![]() |
[4] |
范力更. 我国植保机械和施药技术的现状、问题及对策[J]. 农业开发与装备, 2018, 8(2): 15-21. DOI:10.3969/j.issn.1673-9205.2018.02.013 ( ![]() |
[5] |
任少伟, 杨传民, 陈国营, 等. 大豆蛋白液喷雾场粒度分布变化[J]. 包装工程, 2018, 39(17): 46-52. ( ![]() |
[6] |
吴亚垒, 祁力钧, 张亚, 等. 基于驻波与ZigBee实时监测雾滴蒸发系统设计与试验[J]. 农业工程学报, 2017, 33(17): 128-135. DOI:10.11975/j.issn.1002-6819.2017.17.017 ( ![]() |
[7] |
张东彦, 兰玉彬, 陈立平, 等. 中国农业航空施药技术研究进展与展望[J]. 农业机械学报, 2014, 45(10): 53-59. DOI:10.6041/j.issn.1000-1298.2014.10.009 ( ![]() |
[8] |
王沛, 祁力钧, 李慧, 等. 植物叶片表面结构对雾滴沉积的影响分析[J]. 农业机械学报, 2013, 44(10): 75-79. DOI:10.6041/j.issn.1000-1298.2013.10.013 ( ![]() |
[9] |
兰玉彬, 张海艳, 文晟, 等. 静电喷嘴雾化特性与沉积效果试验分析[J]. 农业机械学报, 2018, 49(4): 130-139. DOI:10.6041/j.issn.1000-1298.2018.04.015 ( ![]() |
[10] |
张慧春, GARY D, 郑加强, 等. 扇形喷头雾滴粒径分布风洞试验[J]. 农业机械学报, 2012, 43(6): 53-57. DOI:10.6041/j.issn.1000-1298.2012.06.010 ( ![]() |
[11] |
张京, 宋坚利, 何雄奎, 等. 扇形雾喷头雾化过程中雾滴运动特性[J]. 农业机械学报, 2011, 42(4): 66-69. ( ![]() |
[12] |
张慧春, 郑加强, 周宏平, 等. 农药喷施过程中雾滴沉积分布与脱靶飘移研究[J]. 农业机械学报, 2017, 48(8): 114-122. DOI:10.6041/j.issn.1000-1298.2017.08.012 ( ![]() |
[13] |
AUNG N Z, YANG Q, CHEN M, et al. CFD analysis of flow forces and energy loss characteristics in a flapper-nozzle pilot valve with different null clearances[J]. Energy Conversion and Management, 2014, 83(45): 284-295. ( ![]() |
[14] |
DE DOMENICO F, ROLLAND E O, HOCHGREB S. A generalised model for acoustic and entropic transfer function of nozzles with losses[J]. Journal of Sound and Vibration, 2019, 440(43): 212-230. ( ![]() |
[15] |
文晟, 兰玉彬, 张建桃, 等. 农用无人机超低容量旋流喷嘴的雾化特性分析与试验[J]. 农业工程学报, 2016, 32(20): 85-93. DOI:10.11975/j.issn.1002-6819.2016.20.011 ( ![]() |
[16] |
王果, 龚燕, 张晓, 等. 不同施药机具在玉米田间的雾滴沉积分布试验[J]. 农机化研究, 2017(6): 177-182. DOI:10.3969/j.issn.1003-188X.2017.06.035 ( ![]() |
[17] |
朱增强, 王冉冉. 烟雾机喷嘴雾化仿真研究[J]. 装备制造技术, 2018(7): 51-54. DOI:10.3969/j.issn.1672-545X.2018.07.017 ( ![]() |
[18] |
祁力钧, 傅泽田. 不同条件下喷雾分布试验研究[J]. 农业工程学报, 1999, 15(2): 113-117. ( ![]() |
[19] |
张瑞瑞, 张真, 徐刚, 等. 喷雾助剂类型及浓度对喷头雾化效果影响[J]. 农业工程学报, 2018, 34(20): 36-43. DOI:10.11975/j.issn.1002-6819.2018.20.005 ( ![]() |
[20] |
王潇楠, 何雄奎, 宋坚利, 等. 助剂类型及浓度对不同喷头雾滴飘移的影响[J]. 农业工程学报, 2015, 31(22): 49-55. DOI:10.11975/j.issn.1002-6819.2015.22.007 ( ![]() |
[21] |
何玲, 王国宾, 胡韬, 等. 喷雾助剂及施液量对植保无人机喷雾雾滴在水稻冠层沉积分布的影响[J]. 植物保护学报, 2017, 44(6): 1046-1052. ( ![]() |
[22] |
唐青, 陈立平, 张瑞瑞, 等. 高速气流条件下标准扇形喷头和空气诱导喷头雾化特性[J]. 农业工程学报, 2016, 32(22): 121-128. DOI:10.11975/j.issn.1002-6819.2016.22.017 ( ![]() |
[23] |
赵辉, 宋坚利, 曾爱军, 等. 喷雾液动态表面张力与雾滴粒径关系[J]. 农业机械学报, 2009, 40(8): 74-79. ( ![]() |
[24] |
肖丽萍, 刘木华, ZHU HEPING, 等. 喷嘴喷施不同生物农药雾滴特性研究[J]. 农业机械学报, 2018, 49(2): 100-106. DOI:10.6041/j.issn.1000-1298.2018.02.013 ( ![]() |
[25] |
LEE C H, REITZ R D. An experimental study of the effect of gas density on the distortion and breakup mechanism of drops in high speed gas stream[J]. International Journal of Multiphase Flow, 2000, 26(2): 229-244. DOI:10.1016/S0301-9322(99)00020-8 ( ![]() |