1. 引言
农业生产中严重依赖除草剂来清除杂草 [1] ,除草剂的滥用使得大量除草剂及其代谢产物进入生态循环中,这些有毒物质会经过食物链逐级传递并不断蓄积,最终进入人体,对人类健康构成威胁。为了保护生态环境与维护人体健康,人们开始研究替代性杂草清除方法。
热除草是替代方法之一,其主要通过高温破坏杂草组织来抑制杂草生长。与其他热源(火焰、蒸汽、阳光等)相比,激光能够更好的聚集能量精确的作用于杂草的某一部位,为热除草提供了新的发展方向 [2] 。
Heisel T使用CO2激光器对三种杂草的茎部进行了辐照,研究表明激光作用于杂草茎部可对杂草生长产生抑制效果 [3] 。Solvejg K等人使用连续激光对三种杂草的顶端分生组织进行了辐照,实验结果表明通过将激光束更精确地指向作用位置,可以提高激光抑制杂草生长的效果 [4] 。C. Wöltjen和T. Rath等人分别使用了二极管激光器与CO2激光器,对两种杂激光草进行了辐照研究,发现激光的波长是影响激光处理杂草效果的主要因素,且处理不同生长阶段的杂草所需的激光能量也不同 [5] 。Marx Christian分别使用了连续激光器与脉冲激光器对杂草反枝苋的分生组织进行了辐照。研究结果表明,对杂草造成致命损伤所需的激光能量取决于能量耦合效率与激光的精度 [6] 。Yu‑Pin Lan使用了二极管激光器对薇甘菊茎部进行辐照,研究结果表明,杂草茎部的凹陷深度取决于激光的功率密度,且高功率密度的单光束激光对于抑制杂草的生长更加有效 [7] 。Guy Coleman等人对于激光清除一年生黑麦草进行了研究 [8] 。研究发现,杂草的生长阶段和杂草的大小对激光功效的影响很大。需要针对不同生长阶段的杂草采用不同的功率,并且要根据杂草植株的大小来调整激光光斑的大小。Ya Xiong等人为了证明激光除草的可行性,建造了一个机器人 [9] ,并且他们开发利用了多种算法来进行分析与研究。为了评估机器人的性能,在室内进行了一组除草实验。测试结果表明,对杂草进行640 ms的90 W、810 nm激光照射可有效控制杂草,该机器人具有高达97%的命中率。
上述研究主要对激光清除杂草的效果以及影响因素进行了研究,证实了激光清除杂草的可行性,但在相互作用过程中,缺少对激光与杂草温度与损伤程度的变化规律的研究。针对这一问题,本文建立了一种包含了杂草受激光辐照升温与杂草受热相变的理论模型,并以杂草马齿苋为研究对象,开展了模拟仿真研究,同时设计与开展了激光辐照杂草马齿苋茎部的实验,研究了杂草马齿苋在受到不同功率密度激光辐照下的杂草温度与损伤程度的变化规律。研究结果可为激光除草的参数选择提供一定的参考。
2. 模型
根据杂草茎部侧面形状,建立二维矩形模型,同时将激光视为一个位置固定的高斯热源作用于杂草茎部。杂草茎部受到激光辐照后,吸收激光能量温度升高,同时受激光辐照区域发生相变形成损伤区域,杂草的物理模型如图1所示。
设定如下:(1) 杂草的组成成分为水与纤维素 [10] ,其内部孔隙被液体充满,故将其视为饱和多孔介质;(2) 由于激光辐照杂草的时间很短,因此将内部水分视为不动流体。
3. 控制方程
(1) 能量方程:
(1)
其中:
(2)
(3)
上面两个参数分别为多孔介质恒定压力下的有效体积热容与有效导热系数。
式子,T为多孔介质温度(K),t为时间(s),Q为热源(W/m3),ε为孔隙率,ρ为密度(s为固体基质,f为液体) (kg/m³),
为恒压热容(s为固体基质,f为液体) (J/(kg·K)),λ为热导率(s为固体基质,f为液体) (W/(m·K))。
(2) 边界条件
① 热对流
(4)
式中,
为对流热通量(W/m2),h为传热系数(W/(m2·K)),T为表面温度(K),
为环境温度(K)。
② 热辐射
(5)
式中,
为表面向环境辐射的热通量(W/m2),
为表面发射率,
为Stefan-Boltzmann常数(5.67 × 10−8 W/(m2·K2))。
(3) Stefan条件
(6)
其中
(7)
(8)
上面两个参数分别为多孔介质的有效密度与多孔介质的有效潜热。
式中,
是界面处法向热通量的突变,在激光与杂草相互作用的过程中,
为激光通过杂草表面向杂草内部传递的能量(W/m2)。
为潜热(s为固体基质,f为液体) (J/kg),
为相变界面速度(m/s)。
(4) 激光热源
激光的光强分布在二维模型中为:
(9)
式中,Q为激光热源(W/m²),P为激光功率(W),
为激光焦点位置坐标,r为激光光斑半径(m)。
4. 实验研究
使用1064 nm连续激光辐照杂草茎部,研究激光功率密度大小(表1)对杂草马齿苋温度以及凹陷程度的影响,实验装置图如图2所示。整个实验系统由1064 nm连续激光器、分光镜、激光功率计、聚焦透镜、热像仪、数字延时发生器以及控制电脑组成。将热像仪与激光器连接到数字延时发生器,控制两者同时开启,使用激光功率计测量激光的功率,并使用热像仪记录杂草茎部中心区域的最高温度。实验结束后,对杂草的损伤程度进行测量。
5. 仿真与实验结果分析
(1) 凹陷程度
图3为激光辐照时间为2 s时,受到不同功率密度激光辐照后的杂草茎部的凹陷深度。从图3中可以
![](//html.hanspub.org/file/108-2571550x29_hanspub.png?20240603091103494)
Figure 3. The depth of the depression in the central area of the weed stem
图3. 杂草茎部中心区域的凹陷深度
![](//html.hanspub.org/file/108-2571550x30_hanspub.png?20240603091103494)
Figure 4. Stems of purslane irradiated by laser
图4. 受激光辐照后的马齿苋茎部
看出,实验测得的凹陷深度与仿真模拟得到的凹陷深度的变化趋势相同。整体变化趋势为杂草茎部的凹陷深度随着激光功率密度的增加而不断加大,当激光功率密度为33.310 W/mm2与39.720 W/mm2时,曲线上升速度略有提高,凹陷深度随着功率密度的变化近似线性增长。这说明在相同时间的激光辐照下,提升激光功率密度可有效的提升激光对杂草的破坏效果。
此外,实验测得的凹陷深度比仿真得到的凹陷深度更深。从图4可以看出,杂草茎部除了因受到激光辐照而直接损伤的区域外,受激光辐照的边缘区域也会发生轻微萎缩。仿真仅对激光与杂草相互作用的过程进行了研究,对植物组织的后续变化没有进行进一步的研究,因此实验测得的值比仿真得到的值更大。用实验所得深度减去仿真所得深度后,将所得差值除以实验值,计算得到误差为5%~10%,表明仿真结果与实验结果具有较好的一致性,该模型可为激光的参数选择提供一定的参考作用。
(2) 温度变化
![](//html.hanspub.org/file/108-2571550x31_hanspub.png?20240603091103494)
Figure 5. Variation of maximum temperature in the central region of weed stems with laser irradiation time
图5. 杂草茎部中心区域最高温度随激光辐照时间的变化
图5为在受到不同功率密度的连续激光辐照下,激光辐照时间为0~2 s的杂草茎部中心区域最高温度变化曲线图。所有曲线均可分为两个阶段,一为温度上升阶段,二为温度稳定阶段。当激光功率密度为47.761 W/mm2与39.720 W/mm2,辐照时间为0~0.5 s时,杂草茎部中心区域最高温度快速上升;辐照时间为0.5~1 s时,温度上升速度变缓;辐照时间为1~2 s时,温度趋于稳定,进入平台期。当激光功率密度为33.310 W/mm2与27.706 W/mm2,激光辐照时间为0~1.5 s时,杂草茎部中心区域最高温度开始逐渐上升;辐照时间为1.5~2 s时,温度趋于稳定,进入平台期。
根据上述曲线变化规律可知,杂草温度上升的速度与激光功率密度成正比,当激光功率密度提升了1.72倍时,杂草温度达到平台期所需时间缩短了3倍,说明高功率密度的激光可以更加快速地使杂草温度达到平台期。当所有温度变化曲线达到平台期后,杂草茎部温度所能达到的最高值相同,这说明一味地增加激光功率密度无法使杂草温度突破平台期。
从图5中还可以看出,仿真与实验所得变化曲线图总体趋势相同。在相同功率密度的激光辐照下,实验与仿真结果在杂草茎部中心区域最高温度进入平台期的时间上略有差异,实验比仿真达到平台期的时间更短。当杂草茎部中心区域最高温度趋于稳定时,仿真所得温度略高于实验测得的温度。将仿真温度与实验温度的差值除以实验值,得到其误差约为2%,说明仿真与实验结果吻合良好。该模型可用于预测杂草受激光辐照后,杂草茎部中心区域最高温度随时间的变化。
6. 结论
本文建立了激光与杂草相互作用的物理模型,基于该模型开展了数值模拟研究,并使用1064 nm连续激光器开展了不同功率密度的激光辐照杂草马齿苋茎部的损伤实验,获得了不同激光功率密度下,杂草中心区域最高温度和损伤程度的变化规律。杂草茎部的凹陷深度随着功率密度的变化近似线性增长,增加激光的功率密度可有效地提升激光对杂草的破坏效果。随着激光功率密度的增加,杂草中心区域的最高温度能更快的达到平台期,但无法使杂草温度突破平台期。
通过对比仿真与实验结果可知,二者吻合良好。表明该模型可用于预测激光辐照杂草茎部后,杂草茎部的凹陷深度与中心区域最高温度,为激光辐照杂草的结果预测提供了一定的参考作用。