冬小麥返青期地下害蟲危害空間分布型及其抽樣技術

李平

摘要：采用空間分布型檢驗、聚集強度指標檢驗和線性回歸方法研究了冬小麥返青期枯心苗空間分布規律及其抽樣方法。結果表明，冬小麥返青期枯心苗空間分布（即地下害蟲危害空間分布）呈聚集分布，聚集程度受種植環境影響較大。建立了冬小麥早春枯心苗最適抽樣模型及序貫抽樣模型。

關鍵詞：冬小麥;地下害蟲;空間分布型;理論抽樣模型;序貫抽樣模型

中圖分類號：S435.122 ? ?文獻標志碼：A ? ? 文章編號：1001-1463（2021）07-0009-05

doi：10.3969/j.issn.1001-1463.2021.07.003

Abstract：The spatial distribution and sampling of the withered wheat seedling in the period of returning to green was studied by using the spatial distribution test， the aggregation intensity index test and the linear regression methods. The results showed that the spatial distribution of the dead wheat seedlings（or the underground pest damages） was clustered. The degree of aggregation was greatly affected by the planting environment. The theoretical sampling model of early winter wheat dead seedlings and its sequential sampling model were established.

Key words：Winter wheat; Underground pests; Spatial distribution; Theoretical sampling model; Sequential sampling model

小麥是世界上播種面積最大、分布最廣泛的糧食作物，也是甘肅省第一大糧食作物。甘肅從東到西均有小麥種植，尤其是冬小麥年種植面積較大，約占全省小麥播種面積的2/3[1 - 2 ]。河西灌區是甘肅省冬小麥主要生產區域之一，冬小麥播種面積占全省的30%以上，冬小麥也是河西地區的傳統產業和主栽作物種類之一[3 - 4 ]。近年來，河西冬小麥生產正在向優質、特色、規模化、精深加工和產業現代化方向發展，部分市縣冬小麥產業已經在當地農村經濟發展和鄉村振興建設中起到了舉足輕重的作用。然而，隨著我國城鎮化的推進，耕地面積數量限制，資源環境約束加劇，國際糧食供求平衡形勢日趨嚴峻，糧食生產和糧食安全已經被擺到了極端重要的地位。因此，保障河西冬小麥生產安全對全省乃至全國糧食安全都具有十分重要的意義。調查發現，早春地下害蟲危害是制約河西冬小麥生產安全的早期重要因素之一。目前關于河西地區冬小麥早春地下害蟲危害空間分布規律及其預測預報的研究鮮有報道[5 - 7 ]。筆者在武威市涼州區研究了冬小麥返青期地下害蟲危害空間分布型及其抽樣方法，旨在為河西地區冬小麥地下害蟲的科學防治和預測預報提供參考。

1 ? 材料與方法

1.1 ? 調查地點和方法

調查地點位于武威市涼州區高壩鎮六壩村，當地平均海拔1 649 m，年均降水161 mm。土壤類型薄層灌漠土，耕層土壤有機質含量14.3 g/kg。指示冬小麥品種為寧冬6號，播種量450～525 kg/hm2。2021年3月10 — 11日選擇冬小麥返青期隨機調查了6塊田，用枯心苗數量反映田間總體地下害蟲危害空間分布，每塊田面積667～900 m2。每塊田即1個樣本田塊，每個樣本田塊均按棋盤式橫向選擇5個點，縱向選擇10個點，每個點為1個樣方，樣方面積0.25 m2。每個樣田調查50個樣方，統計每樣方枯心苗數量，制作χ 2頻次表。

1.2 ? 空間分布型檢驗

1.2.1 ? 聚集度指標檢驗 ? 采用擴散系數C、Cassie指標CA、Lloyd聚集指標M*/x、David&Moore叢生指數I以及聚集均數λ檢驗空間分布型。

1.2.2 ? 線性回歸檢驗 ? 以平均擁擠度M*與平均密度值做Iwao回歸M*=α+β。α為基本擴散指數，β為密度擴散系數。當α > 0，個體間相互吸引，分布的基本成分是個體群;當α=0，分布的基本成分是單個個體;當 ? ? ?α < 0，個體間相互排斥。當β=1時，隨機分布;當β < 1時，均勻分布;當β > 1時，聚集分布。以方差S2與平均密度取對數值后做Taylor回歸lg（S2）=lga+blg（）。當b=1時，空間分布為隨機分布;當b>1時，空間分布為聚集分布;當b趨近于0時，空間分布為均勻分布。

1.3 ? 理論抽樣模型和序貫抽樣模型

Iwao理論抽樣模型n=t2/D2[（α+1）/+ β-1]，n為最適抽樣數或理論抽樣數，為平均密度，D為允許誤差，t為置信度分布值，α、β同Iwao回歸模型參數。

Iwao序貫抽樣模型T（1、2）=nm0 ± t，加號計算可得枯心苗密度上限值T1，減號計算可得枯心苗密度下限值T2。n即抽樣數，m0為防治指標，t為置信度分布值，一般取95%置信區間即t=1.96;α、β同Iwao理論抽樣模型參數。田間調查時，若累計查得枯心苗數量大于上限值T1，表明地下害蟲發生量高于防治指標，需防治;若累計查得枯心苗數量低于下限值T2，表明地下害蟲發生量低于防治指標，不需防治;若累計查得枯心苗數量處于上下限值之間，應繼續取樣調查。

最大抽樣數模型Nmax=t2/{d2[（α+1）m0+（β-1）m02）]}，d即允許誤差D，m0、t、α、β同Iwao序貫抽樣模型參數[8 - 10 ]。當田間調查到最大抽樣數時，若累計查得枯心苗數量仍在上下限之間，則根據該點最靠近的邊界限值判斷是否需要防治。

采用Excel 2003和DPS17.10軟件處理試驗數據。

2 ? 結果與分析

2.1 ? 理論分布型檢驗

從表1可知，1、6號田χ2值均小于該自由度下負二項分布P0.05時的χ2值，表示1、6號田枯心苗實際分布與負二項分布模型顯著相符。2～5號田χ2值均小于該自由度下P-E分布P0.05時的χ2值，表示2～5號田枯心苗實際分布與P-E分布模型顯著相符。其中2、4號田χ2值均小于該自由度下泊松分布P0.05時的χ2值，表示2、4號田枯心苗實際分布也與泊松分布模型顯著相符，3號田χ2值也小于該自由度下奈曼分布、負二項分布P0.05時的χ2值，表示3號田枯心苗實際分布也與奈曼分布、負二項分布模型顯著相符;反之，上述各田χ2值大于該自由度下相應模型χ20.05值，則表示該田枯心苗實際分布與相應模型顯著不符。奈曼分布、P-E分布和負二項分布均為聚集分布，因此可得出1～6號田枯心苗空間分布（或地下害蟲危害空間分布）均顯著呈聚集分布。

2.2 ? 空間分布指標檢驗

由表2可見，1～6號田的擴散系數C > 1、Lloyd聚集指標M*/x > 1、Cassie指數CA > 0、David&Moore叢生指數I > 0，表示1～6號田枯心苗空間分布（或地下害蟲危害空間分布）呈聚集分布;上述樣田聚集均數λ < 2，表示冬小麥返青期枯心苗田間聚集是受環境因素決定。聚集均數λ和平均密度回歸模型為λ=0.937 2-0.144，R2=0.896 4，F=34.6 > F0.01，表示枯心苗聚集程度（或地下害蟲危害聚集程度）與平均密度極顯著正相關。

2.3 ? 線性回歸檢驗

2.3.1 ? Iwao回歸 ? 平均擁擠度M*和平均密度的值回歸模型顯著，方程式為M*=0.204 1+ 1.069，R2=0.685 6，F=8.7 > F0.05。式中基本擴散指數α=0.204 1 > 0，表示枯心苗個體間相互吸引，分布的基本成分是個體群;密度擴散系數β=1.06 9 > 1，表示枯心苗的空間分布型呈聚集分布。

2.3.2 ? Taylor回歸 ? 方差S2和平均密度的對數值的回歸模型極顯著，方程式為lg（S2）= ? 0.088 9+1.001 7 lg（ x ），R2=0.934 2，F=29.9 > ?F0.01。式中b=1.001 7 > 1，表示枯心苗發生分布呈聚集分布。

2.4 ? 理論抽樣模型與序貫抽樣模型

根據Iwao回歸模型和Iwao理論抽樣模型，一般取95%置信度（即t=1.96），可得出冬小麥田早期枯心苗最適抽樣模型n= ? ? 3.841 6/D2（1.204 1/+0.069）。一般允許誤差D可取0.1、0.2和0.3，可得出相應枯心苗密度（例如平均密度=1、2、3、4、5、6、7、8、9、10株/0.25 m2）的最適抽樣數（表3）。若密度相同，則抽樣數量隨著允許誤差增大而減少;若允許誤差相同，則抽樣數量隨密度增加而減少。

根據Iwao序貫抽樣模型，若選擇枯心苗防治指標為4株/0.25 m2，即m0=4;取95%置信區間即t=1.96，可得出其序貫抽樣方程T（1、2）=4n±4.713。例如取調查樣方數n=10、20、30、40、50、60、70、80、90、100時，可得到相應枯心苗數量的序貫抽樣表（表4）。在田間調查中，若累計查得枯心苗數量大于表中上限值T1，即地下害蟲風險高于防治指標，需要開展防治;若累計查得枯心苗數量小于表中下限值T2，即地下害蟲風險低于防治指標，不需防治;若累計查得枯心苗數量處于T1和T2之間，仍需繼續取樣調查。

根據最大抽樣模型，取95%置信度即t=1.96，一般允許誤差d可取0.1、0.2或0.3。例如d=0.1時，可得出Nmax≈2 274，即在防治指標4±0.1株/0.25 m2時田間調查的最大抽樣數是2 274個;d=0.2時，可得出Nmax≈569，即在防治指標4±0.2株/0.25 m2時田間調查的最大抽樣數是569個;d=0.3時，可得出Nmax=253，即在防治指標4±0.3株/0.25 m2時田間調查的最大抽樣數是253個。實際應用中，在一定允許誤差內調查到最大抽樣數時，若累計查得枯心苗數量仍在T1和T2之間，則根據該數值最靠近的邊界限值來決定是否開展防治。

3 ? 結論與討論

田間調查和相關檢驗表明，甘肅金武地區冬小麥返青期枯心苗空間分布（或地下害蟲危害空間分布）型呈聚集分布的結論，該結論與大麥田地下害蟲危害空間分布規律一致[11 - 13 ]。冬小麥早春枯心苗田間聚集分布（或地下害蟲危害空間分布）受種植環境的影響較大，且聚集程度與密度極顯著正相關。建立了冬小麥返青期枯心苗最適抽樣模型n=3.8416/D2（1.204 1/x+0.069）和序貫抽樣模型T （1、2）=4n±4.723。

本研究建立的序貫抽樣方法可供河西地區冬小麥種植大戶或基層專業化統防統治組織參考。在實際應用中，可根據序貫抽樣表適時組織和開展冬小麥枯心苗調查，決定是否開展冬小麥田地下害蟲防治。

參考文獻：

[1] 楊長剛，楊文雄，王世紅. ?甘肅省小麥產業發展對策[J]. ?中國種業，2017（11）：1-5.

[2] 李貴喜，白靈軍. ?甘肅隴東冬小麥安全生產現狀及對策[J]. ?甘肅農業科技，2017（4）：68-72.

[3] 倪勝利，李興茂. ?冬小麥新品種隴鑒110選育報告[J]. ?甘肅農業科技，2020（4）：11-14.

[4] 馬小麗. ?小麥地下害蟲的識別與防治[J]. ?河南農業，2014（3）：38.

[5] 張 ? 菁. ?彬縣小麥地下害蟲為害特點及其綜合防控措施[J]. ?陜西農業科學，2015，61（5）：70-71.

[6] 楊榮明. ?小麥地下害蟲發生特點及防治方法[J]. ?農家致富，2015（6）：34-35.

[7] 劉效華，劉建華，虎夢霞. ?甘肅河西地區冬小麥光合生產潛力分析[J]. ?甘肅農業科技，2010（9）：9-12.

[8] 李 ? 平，戴 ? 偉. ?藜科雜草在洋蔥育苗田的空間分布型及其抽樣技術[J]. ?甘肅農業科技，2021，52（4）：49-52.

[9] 李 ? 平. ?洋蔥根腐病在育苗初期的空間分布型及抽樣技術[J]. ?甘肅農業科技，2021，52（5）：26-29.

[10] 李 ? 平，戴 ? 偉. ?潛葉蠅幼蟲在二月蘭的田間空間分布型及其抽樣[J]. ?甘肅農業科技，2021，52（5）：53-57

[11] 劉凡杰，張清旺，王恩明. ?冬小麥主要地下害蟲分布特點與防治方法[J]. ?河北農業科技，2008，9（17）：26.

[12] 李耀發，黨志紅，安靜杰，等. ?河北省主要作物田地下害蟲種類及其分布[J]. ?中國農學通報，2018，34（28）：114-119.

[13] 王文娟，唐小琴，孟長軍. ?裸大麥田蠐螬的空間分布型及抽樣技術研究[J]. ?江蘇農業科學，2012，40（8）：105-107.

（本文責編：楊 ? 杰）