基于多旋翼植保無人機(jī)的板栗飛防效果研究*

魏源，蘇強(qiáng)，曾揚(yáng)鵑，杜菁菁，齊永順，4，王東升，4

(1.河北科技師范學(xué)院園藝科技學(xué)院,河北秦皇島,066004;2.板栗產(chǎn)業(yè)技術(shù)教育部工程研究中心,河北秦皇島,066004;3.河北省特色園藝種質(zhì)挖掘與創(chuàng)新利用重點實驗室,河北秦皇島,066004;4.河北省板栗產(chǎn)業(yè)協(xié)同創(chuàng)新中心,河北秦皇島,066004)

0 引言

植保無人機(jī)由于其在作業(yè)效率、節(jié)水節(jié)時、適應(yīng)性、靈活性等多個方面具有優(yōu)勢,近年來得到了較快速的發(fā)展[1],在玉米、小麥等大田作物應(yīng)用較多[2],對水稻、棉花的病蟲防治有著非常顯著的效果,然而受地形地貌及作物自身屬性的影響,限制了植保無人機(jī)在山區(qū)的柑橘、茶葉、蘋果、板栗、桃樹等經(jīng)濟(jì)作物上的應(yīng)用,相關(guān)報道還較少[3]。

目前,由于板栗單產(chǎn)低、效益差、立地條件惡劣,使得板栗的生產(chǎn)多處于靠天吃飯的窘境,其栽培管理也多為粗放管理模式。導(dǎo)致栗園后期郁閉嚴(yán)重,產(chǎn)量低,效益差。據(jù)估計,板栗低產(chǎn)低效林占全國種植面積的80%左右,年均單產(chǎn)僅有70 kg左右。病蟲防治是板栗栽培管理中的重要環(huán)節(jié),當(dāng)前模式耗費人力巨大,管理用工費用逐漸增高,人工施藥對身體傷害極大,生態(tài)環(huán)境逐漸失衡等,極大制約著板栗產(chǎn)業(yè)可持續(xù)健康發(fā)展。通過植保無人機(jī)飛防大田作物已經(jīng)逐步成熟,然而其對于果樹的研究卻少之又少,果樹生長環(huán)境惡劣,且人力物力耗費巨大,如果可以通過試驗選定出針對果樹的最佳飛行參數(shù),并對其進(jìn)行深入研究是非常有必要的。本研究利用大疆植保無人機(jī)T16、T20、T30,在秦皇島市撫寧區(qū)進(jìn)行了多旋翼植保無人機(jī)飛防霧滴效果優(yōu)選試驗,制定12種不同的參數(shù)搭配組合,探索其在板栗飛防中的可行性及作業(yè)效果,為檢驗無人機(jī)在開心型板栗園噴灑農(nóng)藥和防治有害生物效果,制定板栗無人機(jī)飛行作業(yè)標(biāo)準(zhǔn)提供參考。

1 材料與方法

1.1 試驗材料

施藥機(jī)器：多旋翼植保無人機(jī)型號為大疆T16、T20、T30,利用扇形壓力噴頭,超低量噴霧飛防施藥。用藥量為119.40 L/hm2,T16最大噴頭流量4.8 L/min,T20最大噴頭流量6 L/min,T30最大噴頭流量8 L/min。

檢測用具：正品精準(zhǔn)科研數(shù)據(jù)分析專用水敏紙霧滴測試卡,萬向夾(尺寸8 cm)、自封袋(尺寸10 cm×15 cm)、標(biāo)簽貼紙(常規(guī)標(biāo)準(zhǔn)尺寸)以及希瑪高精度手持式風(fēng)力測試儀AS836。

1.2 試驗設(shè)計

該試驗于2021年10月27日進(jìn)行,試驗果園選在河北省秦皇島市撫寧區(qū)淼源堅果合作社板栗園內(nèi),位于N39°72′,E119°22′,地處燕山山脈東段丘陵地區(qū)與山前平原地帶。當(dāng)日平均氣溫為15 ℃,空氣濕度為62%,自然風(fēng)速為1.2 m/s。選擇3棵土壤肥力、果樹長勢相同,樹體指標(biāo)接近,土肥水管理條件一致的板栗樹,試驗品系為“撫寧15號”,株行距為3 m×4 m。樹體指標(biāo)如表1所示。本試驗共設(shè)12個處理,每個處理3次重復(fù),其中飛防作業(yè)根據(jù)大疆T16、T20、T30植保無人機(jī)的基本性能(技術(shù)參數(shù)),按照噴霧作業(yè)質(zhì)量要求,結(jié)合板栗園實際作業(yè)情況和工作經(jīng)驗,設(shè)計處理參數(shù)如表2所示。

表1 樹體指標(biāo)Tab.1 Tree body index cm

表2 試驗設(shè)計Tab.2 Experiment design

1.3 試驗方案

為保證試驗的準(zhǔn)確性,先進(jìn)行噴幅預(yù)試驗,試驗田內(nèi)選定矩形區(qū)塊,設(shè)置一條長15 m,相距1 m平行的霧滴采集帶,從中心線橫向每隔1 m布置一根鋼管,隔離帶共15個點,將霧滴測試卡固定在每個點位。在無風(fēng)的情況下,操作植保無人機(jī)沿中間飛行,模擬噴霧條件。待預(yù)試驗結(jié)束后,將測試卡帶回實驗室,用霧滴計數(shù)卡計算出每張霧滴測試卡上的霧滴數(shù)(個/cm2),如果霧滴數(shù)高于20個/cm2,則為有效霧滴沉積密度,兩張最低有效霧滴沉積數(shù)卡片的距離,為測定植保無人機(jī)噴霧的航幅。

根據(jù)無人機(jī)起飛路徑,在距離起飛點30 m處,開始選定兩兩間距10 m的三棵試驗樹,此時飛機(jī)運行穩(wěn)定,速度高度已達(dá)到初定標(biāo)準(zhǔn)。選定樹冠中部南向和東向的主枝,沿主枝伸長方向,從內(nèi)到外分別水平放置三張測試卡。使用萬向夾將測試卡與葉片固定結(jié)實,確保正面朝上,保證與地面平行。操作無人機(jī)沿樹木中線飛行(用清水代替藥劑進(jìn)行噴霧),嚴(yán)格監(jiān)控記錄飛行時間、噴灑面積,填寫飛行記錄,同時使用風(fēng)力測試儀進(jìn)行風(fēng)力、風(fēng)溫測定。噴灑結(jié)束以后,按要求收取測試卡,進(jìn)行信息提取。

1.4 數(shù)據(jù)處理

將數(shù)據(jù)用六六山下霧滴軟件進(jìn)行霧滴分析,得出不同處理的DV1(取樣霧滴的體積按照從小到大進(jìn)行排序累積,當(dāng)累積值達(dá)到10%時所對應(yīng)的霧滴直徑)、VMD(取樣霧滴的體積按照從小到大進(jìn)行排序累積,當(dāng)累積值為50%時所對應(yīng)的霧滴直徑)、NMD(取樣霧滴的個數(shù)按照從小到大進(jìn)行順序累積,累積值為50%時所對應(yīng)的霧滴直徑)、霧滴譜寬度、霧滴密度、覆蓋率、農(nóng)藥沉積量。利用Excel 2010對數(shù)據(jù)進(jìn)行整理和相關(guān)性分析,利用SPSS 24.0進(jìn)行主成分分析,得到主成分分值,以所選主成分對應(yīng)的貢獻(xiàn)率占所提取的2個主成分的貢獻(xiàn)率總和的比例為權(quán)重,計算前2個主成分的分值與相應(yīng)權(quán)重的乘積累加,得到綜合分值F,分析研判各個指標(biāo)之間相關(guān)性。以此來評價不同組合處理的優(yōu)劣[4],選出最優(yōu)處理。

2 結(jié)果與分析

2.1 不同飛行參數(shù)下指標(biāo)效果分析

DV1值越大表明霧滴效果越好[5],不同處理下DV1數(shù)值差異顯著,最高的M10為97.47 μm,最低的M5,僅為M10的四分之一。T16無人機(jī)與T20無人機(jī)各處理之間無顯著差異,保持在24.35～40.59 μm之間,而T30無人機(jī)的DV1均值要比前者高出兩倍之多。當(dāng)無人機(jī)型號相同飛行參數(shù)不同的情況下DV1值均無顯著差異。通過對比不難發(fā)現(xiàn),當(dāng)無人機(jī)高度不變,隨著飛行速度增加,DV1值呈上升趨勢。

VMD值越大表明霧滴效果越好[6],其中M9的VMD值最高,達(dá)到229.93 μm,是最低的M5的2.2倍。對比不同機(jī)型可以發(fā)現(xiàn),T30與T16、T20無人機(jī)的VMD值差異顯著,且明顯優(yōu)于T16、T20。T30不同處理的VMD值均高于200 μm,而T16和T20試驗中最大的VMD值僅為166.71 μm。同一機(jī)型的不同處理間,VMD值雖有所變化,卻無顯著差異。此外,還發(fā)現(xiàn)使用T16和T20進(jìn)行試驗時,當(dāng)一個參數(shù)固定,VMD值會隨著另一個參數(shù)的升高而升高。然而在T30試驗中,當(dāng)一個參數(shù)固定,VMD值隨著另一個參數(shù)的升高而降低,表現(xiàn)了完全不同的變化規(guī)律。

不同處理下NMD數(shù)值的區(qū)間在28.53～143.05 μm之間,最高的為M10,最低的為M5,最高的M10是最低的M5的5倍。通過對比發(fā)現(xiàn),T30無人機(jī)在不同處理下無顯著差異,穩(wěn)定在136.08～143.05 μm之間,反觀T16、T20無人機(jī)在不同處理下則差異顯著,最小值為28.53 μm,最大值高達(dá)64.54 μm,NMD值越大表明霧滴效果越好,T30無人機(jī)要明顯優(yōu)于T16、T20無人機(jī)。在T20無人機(jī)進(jìn)行試驗的過程中,當(dāng)一個參數(shù)固定,NMD值會伴隨另一個參數(shù)呈現(xiàn)上升趨勢。

霧滴粒徑分布均勻度用霧滴譜寬度來表示,值越小表示均勻度越好[7],不同處理下霧滴譜寬度數(shù)值在0.90～1.69之間,最高的為M7,最低的為M9,M9處理霧滴分布最為均勻,效果最優(yōu),霧滴譜寬度為0.90,M7處理高達(dá)1.69。T30無人機(jī)試驗中霧滴譜寬度數(shù)值均穩(wěn)定在1.20以下,而T16無人機(jī)、T20無人機(jī)霧滴譜寬度數(shù)值介于1.21～1.69之間,T30植保無人機(jī)的霧滴粒徑分布均勻度要優(yōu)于T16植保無人機(jī),T16植保無人機(jī)要優(yōu)于T20植保無人機(jī),所有處理均無顯著性差異。當(dāng)T16無人機(jī)進(jìn)行試驗時,高度參數(shù)不變,霧滴譜寬度會隨著飛行速度的增加呈現(xiàn)上升趨勢。當(dāng)T30無人機(jī)進(jìn)行試驗時,無論設(shè)置哪一個參數(shù)為定值,霧滴譜寬度都會伴隨另一個參數(shù)呈現(xiàn)上升趨勢。

不同處理下霧滴密度數(shù)值在29.80～74.58個/cm2之間,最高的為M12,是最低的M5處理的2.5倍。T30無人機(jī)霧滴密度數(shù)值介于63.21～74.58個/cm2之間,T16無人機(jī)與T20無人機(jī)各處理之間無顯著性差異,T30無人機(jī)的霧滴密度均值要高于T20無人機(jī),T20無人機(jī)要優(yōu)于T16植保無人機(jī)。霧滴密度越高防治效果也就越好。當(dāng)T16無人機(jī)進(jìn)行試驗時,高度參數(shù)不變,隨著飛行速度增加霧滴密度呈現(xiàn)下降趨勢,當(dāng)T20無人機(jī)進(jìn)行試驗時,速度不變,隨著飛行高度增加霧滴密度呈現(xiàn)上升趨勢。

不同處理下霧滴覆蓋率在4%～21.33%之間,最高的為M9,是最低處理M2的5.3倍。霧滴覆蓋率越高霧滴效果越好。T30植保無人機(jī)的霧滴覆蓋率均值要優(yōu)于T20、T16植保無人機(jī),T20、T30無人機(jī)各處理間無顯著性差異。當(dāng)T16無人機(jī)試驗高度不變時,霧滴覆蓋率會伴隨著飛行速度的增加呈下降趨勢;然而當(dāng)速度參數(shù)不變時,霧滴覆蓋率又隨著飛行高度增加呈現(xiàn)上升趨勢。T20無人機(jī)高度不變,隨著飛行速度增加霧滴覆蓋率呈現(xiàn)上升趨勢;速度不變,隨著飛行高度增加霧滴覆蓋率呈現(xiàn)下降趨勢。

在葉面上沉積藥量越多,噴灑效果就越好,流失在環(huán)境中的藥量就越少。不同處理下農(nóng)藥沉積量在0.15～1.24 μL/cm2之間,農(nóng)藥沉積量最高的M9,是最低的M2的8倍。T30無人機(jī)各處理霧滴農(nóng)藥沉積量均值為0.825 μL/cm2,T16與T20無人機(jī)各處理農(nóng)藥沉積量均值都低于0.58 μL/cm2,T16無人機(jī)在進(jìn)行試驗時,高度參數(shù)不變,隨著飛行速度增加農(nóng)藥沉積量呈現(xiàn)下降趨勢;當(dāng)T20無人機(jī)速度不變,隨著飛行高度增加農(nóng)藥沉積量呈現(xiàn)下降趨勢。

農(nóng)藥通過無人機(jī)噴灑的霧滴譜寬度、覆蓋率、農(nóng)藥沉積量、覆蓋密度、覆蓋的均勻性都會對藥效發(fā)揮有顯著影響[8]。是判定植保無人機(jī)最適合組合的重要因素,不同飛行參數(shù)下各指標(biāo)差異性分析如表3所示。

表3 各指標(biāo)差異性分析Tab.3 Differential analysis of each indicator

2.2 12個不同處理的主成分分析

利用SPSS24.0對7個試驗指標(biāo)進(jìn)行主成分分析,分別計算了7個指標(biāo)的特征向量、特征值及貢獻(xiàn)率(表4)和成分矩陣(表5)。

表4 七個指標(biāo)的特征向量、特征值及貢獻(xiàn)率Tab.4 Eigenvectors,eigenvalues,and contribution rates of the seven indicators

表5 七個指標(biāo)的成分矩陣Tab.5 Component matrix of the seven indicators

前2個主成分(特征值>1)的累計貢獻(xiàn)率達(dá)到89.087%,即前2個主成分可以代表所測所有指標(biāo)89.087%的信息量,解釋了絕大部分的原始信息。由各特征向量和貢獻(xiàn)率可以得出,第1主成分的貢獻(xiàn)率最大,達(dá)到了71.716%,特征值為5.02,主要由DV1、VMD、NMD、霧滴譜寬度、霧滴密度等因子決定;第2主成分的貢獻(xiàn)率為17.371%,特征值為1.216,主要由農(nóng)藥沉積量因子決定。各指標(biāo)間的相關(guān)性分析可以了解其內(nèi)在聯(lián)系,為篩選最佳處理提供依據(jù)。

通過表6可以發(fā)現(xiàn),DV1和NMD呈極顯著正相關(guān)(R=0.997,P<0.01),霧滴譜寬度和覆蓋率呈極顯著負(fù)相關(guān)(R=0.735,P<0.01),霧滴密度和NMD呈極顯著正相關(guān)(R=0.802,P<0.01),霧滴覆蓋率和農(nóng)藥沉積量呈極顯著正相關(guān)(R=0.848,P<0.01)。

表6 各指標(biāo)之間相關(guān)性分析Tab.6 Correlation analysis among the various indicators

2.3 七個不同指標(biāo)的綜合評價

由表7得分系數(shù)矩陣可以計算出2個主成分的得分,最終算出綜合得分F。從綜合指標(biāo)上看,綜合分值F越高,通過指標(biāo)表現(xiàn)出來的處理越理想。不同處理的主成分得分及綜合評價的分析結(jié)果如表8所示。

表7 成分得分系數(shù)矩陣Tab.7 Component score coefficient matrix

表8 12個處理的綜合得分以及排名Tab.8 Composite scores for the 12 treatments as well as rankings

從表8可以看出,綜合分值排名前5位的分別為M9、M10、M12、M11、M3,說明這幾個處理綜合表現(xiàn)理想;M8、M7、M6、M4、M2得分居中,說明這幾個處理的綜合表現(xiàn)一般;M1、M5評分較低,綜合表現(xiàn)較差。

3 討論

植保無人機(jī)體積較小,能夠適應(yīng)各種場地,在農(nóng)業(yè)病蟲害防治工作中具有安全、高效的優(yōu)勢,能夠提升實際應(yīng)用的便利度,與傳統(tǒng)農(nóng)業(yè)生產(chǎn)相比,使用植保無人機(jī)可以節(jié)約農(nóng)藥和水以及減少污染。目前我國植保無人機(jī)主要機(jī)型有大疆T16、大疆T20、大疆T30、天鷹-3、極飛P20、極飛V40等,其中大疆T16、T20、T30是農(nóng)戶們應(yīng)用最多、性價比最高、總體反饋最好的三款無人機(jī),因此對這三款無人機(jī)進(jìn)行試驗具有切實的可行性和參考性。

近年來關(guān)于無人機(jī)飛行參數(shù)對霧滴沉積效果影響的研究較多[9]。高圓圓[10]與張聰合等[11]研究發(fā)現(xiàn)植保無人機(jī)飛防作業(yè)質(zhì)量受作業(yè)技術(shù)參數(shù)的影響很大,無人機(jī)作業(yè)技術(shù)參數(shù)包括飛行高度、飛行速度、噴幅寬度等,其中飛行高度、飛行速度對霧滴沉積有顯著性影響。合理的飛行參數(shù)組合對于植保無人機(jī)的高效利用是至關(guān)重要的。考慮到植保無人機(jī)的施藥范圍、施藥區(qū)域等問題,在小麥、玉米、高粱等大田作物中一般采用超低容量噴霧技術(shù)進(jìn)行噴施[12-13]。高軍等[14]認(rèn)為無人機(jī)飛行高度不宜過低,若飛行高度較低,旋翼風(fēng)場之外的霧滴受力情況較為復(fù)雜,受到環(huán)境風(fēng)場、旋翼風(fēng)場的邊界效應(yīng)、環(huán)境溫度等多重因素的影響,一部分霧滴在還沒有落到葉片上就已經(jīng)消散,另一部分霧滴由于復(fù)雜風(fēng)力的作用會發(fā)生飄移[15-17],霧滴均勻性變差,藥害風(fēng)險增加,藥液大量沉積地面從而造成浪費,飛行高度一般為1.5～2 m。然而,果樹作為一種高大樹種,一般種植在山坡、丘陵地帶,生長環(huán)境艱苦,人工作業(yè)較為困難,采用植保無人機(jī)開展防治工作可提高工作效率。王聰?shù)萚18]認(rèn)為不同沉降距離下的霧滴覆蓋密度存在非常大的差異,當(dāng)沉降距離從1.5 m增加到3 m時,霧滴覆蓋密度減幅超過65%;顯然無人機(jī)飛防果樹的作業(yè)高度直接影響到了霧滴的沉降效果,飛行高度一般設(shè)定在2 m左右。薛新宇等[19]通過篩選N-3型無人機(jī)施藥參數(shù)研究發(fā)現(xiàn),無人機(jī)在3 m高度以下施藥的效果最好。余文勝等[20]認(rèn)為運用較低速度、較大流量可以更好地使用植保無人機(jī)均勻噴灑至茶樹葉面,當(dāng)速度較快時,霧滴無法有效沉降在樹體內(nèi)部,繼而影響樹葉的沉積與吸收,飛行速度為1 m/s、1.5 m/s時可以讓樹葉更好地吸收藥液。對于開心型板栗樹,冠層厚度較大,當(dāng)飛行高度達(dá)到3 m時,霧滴沉降距離顯著上升,受飄離影響,對于樹冠內(nèi)部葉面影響較大。植保無人機(jī)的噴幅有限,當(dāng)植保無人機(jī)飛行高度低于2 m,又會影響到樹冠外圍樹葉的噴灑效果,所以高度范圍為2～2.5 m。綜上所述,在對板栗樹進(jìn)行飛防效果研究時,飛行速度設(shè)置為1 m/s、1.5 m/s,飛行高度設(shè)置為2 m/s、2.5 m/s。

常見的航空噴施作業(yè)質(zhì)量性能指標(biāo)主要有單位面積霧滴沉積量、霧滴譜寬度、霧滴覆蓋密度及均勻度[21]。余文勝等[20]采用體積中值直徑(VMD)或數(shù)量中值直徑(NMD)等指標(biāo)對霧滴沉積效果進(jìn)行研究。本文基于不同機(jī)型、不同速度、不同高度組合下的VD1、VMD、NMD、農(nóng)藥沉積量、霧滴密度、霧滴譜寬度、覆蓋率七個性能指標(biāo)數(shù)據(jù),利用六六山下霧滴軟件進(jìn)行霧滴分析[22-24],研究發(fā)現(xiàn)其總體趨勢基本一致,所選指標(biāo)可以有效反映植保無人機(jī)的噴霧質(zhì)量,最后結(jié)合主成分分析法進(jìn)行分析,發(fā)現(xiàn)T30植保無人機(jī)的飛防效果顯著優(yōu)于T16、T20植保無人機(jī),T30植保無人機(jī)最理想的飛防參數(shù)為：速度1 m/s,高度2 m,噴幅7 m。

已有多項研究表明植保無人機(jī)噴霧相對于常規(guī)噴霧,工作更加高效、防治費用成本更低,具有良好的應(yīng)用前景[25-27]。伴隨著科技的進(jìn)步,多旋翼無人機(jī)技術(shù)的蓬勃發(fā)展引發(fā)了植保無人機(jī)產(chǎn)業(yè)迅速成熟,目前無人機(jī)研發(fā)機(jī)構(gòu)正在著手精準(zhǔn)噴藥方式提高、先進(jìn)傳感器開發(fā)、變量噴藥控制方式改進(jìn)、霧滴沉降規(guī)律研究、藥物漂移預(yù)防等。無人機(jī)還可以運用于蘋果、葡萄、山楂、櫻桃等一切生存環(huán)境困難的植物,如果可以早期研究甚至可以馬上運用到實用當(dāng)中去,對于農(nóng)藥的節(jié)省,人力物力的大量節(jié)約,資源的利用開發(fā)都有很大的貢獻(xiàn),在不久的將來山區(qū)、牧區(qū),均可以噴灑,甚至還可以涉及到人類無法到達(dá)的區(qū)域,這些都是在日前的文獻(xiàn)中少有記載,未來也將會彌補(bǔ)我國在農(nóng)業(yè)無人機(jī)領(lǐng)域的空白局面。

4 結(jié)論

本文使用大疆T16、T20、T30植保無人機(jī)對板栗樹進(jìn)行噴灑試驗,通過設(shè)置不同的飛行參數(shù),研究了12種不同飛行組合對霧滴密度及覆蓋密度等七個指標(biāo)的影響,并對數(shù)據(jù)進(jìn)行處理分析,提出了無人機(jī)板栗噴灑作業(yè)參數(shù)的優(yōu)化建議,結(jié)論如下。

1) 通過對比數(shù)據(jù)發(fā)現(xiàn),在飛行參數(shù)相同的情況下,大疆T30植保無人機(jī)各項指標(biāo)要優(yōu)于T16、T20植保無人機(jī),T30無人機(jī)的DV1值約為T16、T20的3.1倍,VMD值約為1.56倍,NMD值約為3倍,霧滴密度、覆蓋率、農(nóng)藥沉積量三個指標(biāo)的倍數(shù)分別為1.5倍、1.92倍、2.1倍。

2) 通過SPSS主成分分析法發(fā)現(xiàn),大疆T30植保無人機(jī)在飛行速度為1 m/s、飛行高度為2 m時,總得分最高,分值為1.268,分值最低的為-1.12,飛行參數(shù)是速度為1 m/s、高度為2 m的T20植保無人機(jī),通過對比不同飛行參數(shù)下植保無人機(jī)在開心型板栗園的飛防效果,達(dá)到防治有害生物的效果,建議選擇T30無人機(jī),飛行參數(shù)為速度為1 m/s、高度為2 m。