基于主成分分析降水量指標預測棉鈴蟲百株卵量模型研究

李 胄，唐保善

(西北農林科技大學 農學院，陜西楊凌 712100)

基于主成分分析降水量指標預測棉鈴蟲百株卵量模型研究

李 胄，唐保善

(西北農林科技大學 農學院，陜西楊凌 712100)

棉鈴蟲是中國棉花主要害蟲之一，降水量對其種群數量具重要影響。采用連續10 a棉鈴蟲第2代和第3代的百株卵量與同期降水量指標資料，通過相關系數法篩選出10個預測第2代、第3代棉鈴蟲的降水量指標，采用主成分分析法建立第2代、第3代棉鈴蟲百株卵量的特征預測模型。檢驗結果表明，第2代棉鈴蟲的預測模型信度較好，歷史符合率為70%；而第3代棉鈴蟲的預測模型信度較差，歷史符合率為40%。由于入選因子均為前兆因子，因此，第2代棉鈴蟲的預測模型可用于區域棉鈴蟲發生程度的監測和預警研究。

棉鈴蟲；降水量指標；主成分分析；預測模型

棉鈴蟲[Helicoverpaarmigera(Hübner)]屬鱗翅目(Lepidoptera)夜蛾科(Noctuidae)，是多食性害蟲，寄主植物多達幾百種，其中栽培作物有60多種，栽培寄主植物中，農作物有12種，蔬菜30種，果樹林木7種，其他特種經濟作物11種[1]，棉花是受其危害最嚴重的經濟作物。棉鈴蟲在中國各個棉區常年都有發生，1992年和1993年中國棉鈴蟲空前爆發，對棉花的危害損失分別高達40億元和100億元[2]，嚴重威脅棉區農業和經濟發展。黃河流域棉區每年棉鈴蟲防治次數一般約25次，平均每3 d防治1次，若遇雨天或者所用農藥藥效差等情況，防治次數可達33次，是危害和損失最嚴重的棉區[2]。近年來，包括中國植保專家在內的廣大棉花科技工作者對棉鈴蟲的發生、測報以及防治進行了大量探索[1,3-8]，積累了豐富的資料和經驗，但關于降水量與棉鈴蟲百株卵量關系的研究報道較少，尤其是利用連續多年田間氣象因子中降水量的系統資料預測預報棉鈴蟲的百株卵量的研究報道更罕見。因此，本研究采用原陜西省棉花研究所鄭永善老師調查的連續10 a(1961-1970年)棉田第2代、第3代棉鈴蟲百株卵量原始資料，結合當地的氣象條件，采用主成分分析法建立棉鈴蟲百株卵量的預測模型，并對該模型進行回報檢驗，以期對棉鈴蟲發生發展的氣象條件尤其是降水量進行深入研究，為棉鈴蟲蟲災預警提供科學依據，防范并減輕棉鈴蟲的危害。

1 材料與方法

1.1 百株卵量調查

按照黑光燈誘集的第2代和第3代棉鈴蟲蛾量高峰日，確定調查卵量日期。每年以陜西省咸陽市涇陽縣斗口鎮為圓心，在半徑約5 km范圍內選取10塊有代表性的棉田調查百株卵量，每塊棉田面積≥0.133 hm2。每塊棉田確定并標記5個固定調查點，每點固定調查20株，所有調查均在1 d內完成。

1.2 數據獲取

1961-1970年的降水資料均來自陜西省咸陽市涇陽縣斗口氣象站(N34° 37′，E108° 53′)。

1.3 分析方法

降水因子以日、旬、月等為單位統計。采用線性傾向估計法分析各降水氣象要素資料，采用簡單相關系數和主成分分析法篩選影響棉鈴蟲的關鍵降水氣象因子，并運用多元回歸方法建立測報模型。

主成分分析法：首先將多個研究變量綜合成有限個互不關聯特征指標進行統計計算[9-11]，然后在保持總方差不變的前提下篩選方差貢獻最大的主成分作為研究依據。

設X=(x1，x2，…，xm)′為m維向量，則X的第1、2、…、m主成分定義為：

采用Excel 2003軟件計算相關系數，采用SPSS 19.0軟件進行主成分及多元線性回歸分析，采用Sigmaplot 12.0軟件制圖。

2 結果與分析

2.1 預報因子篩選

第2代棉鈴蟲卵在當地的發生高峰期為6月中下旬，第3代棉鈴蟲卵在當地的發生高峰期為7月中下旬，因此，選擇當年3～6月的22個降水量指標資料，分析其與棉鈴蟲百株卵量的相關關系，采用90%信度水平檢驗篩選出10個預報因子，在當年4～7月的22個降水量指標資料中篩選出10個預報因子，用于主成分分析，各降水量指標與棉鈴蟲第2代、第3代百株卵量的相關系數見表1。

由表1可知，入選的預報因子中，3月降水量、4月降水量、3月單日最大降水量、4月單日最大降水量、3月日降水量≥0.1 mm日數、3月日降水量≥10 mm日數、4月日降水量≥10 mm日數、4月最長連續降水量、6月上旬降水量與第2代棉鈴蟲的百株卵量呈正相關，說明這些降水量指標的升高有利于第2代棉鈴蟲的大發生，而4月最長無連續無降水日數與第2代棉鈴蟲百株卵量呈負相關，說明4月最長連續無降水日數的增加不利于第2代棉鈴蟲的大發生；4月降水量、6月降水量、4月日降水量≥0.1 mm日數、6月日降水量≥5.0 mm日數、6月日降水量≥10 mm日數、4月最長連續降水量、7月上旬降水量、4月單日最大降水量及6月單日最大降水量與第3代棉鈴蟲的百株卵量呈正相關，說明這些降水量指標的升高有利于第3代棉鈴蟲的大發生，而4月最長無連續無降水日數與第3代棉鈴蟲的百株卵量呈負相關，說明4月最長連續無降水日數的增加不利于第3代棉鈴蟲的大發生。

表1 陜西省咸陽市1961-1970年棉鈴蟲百株卵量與降水量指標的相關系數Table 1 The correlation between precipitation index and the amount of Helicoverpa armigeraper one hundred plants in Xianyang city, Shaanxi province from 1961 to 1970

2.2 主成分分析

按照第2代和第3代棉鈴蟲的入選預報因子主成分對總方差的解釋(表2)可知，第2代棉鈴蟲前4個主成分的累計貢獻率為92.767%，基本反映降水量指標的信息，從而可用4個彼此不相關的綜合指標分別綜合存在于原有的10個氣象指標的各類信息，從第2代棉鈴蟲百株卵量的前3個主成分的特征間量、因子載荷分析，第1主成分綜合3月降水量、4月降水量、3月單日最大降水量、3月日降水量≥10 mm日數、4月最長連續降水量、4月最長無連續無降水日數等6個降水量指標的變異信息，第2主成分綜合4月最長連續降水量、6月上旬降水量等2個降水量指標的變異信息，第3主成分反映4月單日最大降水量的變異信息，第4主成分反映3月日降水量≥0.1 mm日數的變異信息。

第3代棉鈴蟲前3個主成分的累計貢獻率為82.835%，基本反映原降水量指標的信息，從而可用3個彼此不相關的綜合指標分別綜合存在于原有的10個氣象指標的各類信息，從第3代棉鈴蟲百株卵量的前3個主成分的特征間量、因子載荷來看，第1主成分綜合4月降水量、6月降水量、6月日降水量≥5.0 mm日數、6月日降水量≥10 mm日數、7月上旬降水量、4月單日最大降水量、6月單日最大降水量等7個降水量指標的變異信息；第2主成分反映4月最長連續降水量的變異信息；第3主成分綜合4月日降水量≥0.1 mm日數、7月上旬降水量等2個降水量指標的變異信息。

表2 降水量指標主成分分析對總方差的解釋Table 2 The total variance explained by principal component analysis

2.3 預測模型建立

采用降水量特征指標第1、2、3、4主成分特征向量構造的新變量Z1、Z2、Z3、Z4作為預報因子，建立其與第2代棉鈴蟲百株卵量的多元回歸方程；用降水量特征指標第1、2、3主成分特征向量構造的新變量Z1、Z2、Z3作為預報因子，建立其與第3代棉鈴蟲百株卵量的多元回歸方程(表3)。其中，第2代棉鈴蟲的預測方程信度較好(P=0.046)，而第3代棉鈴蟲的預測方程信度較差(P=0.386)。由圖1可知，第2代棉鈴蟲模型的歷史符合率較高，為70%，而第3代棉鈴蟲模型的歷史符合率較低，為40%。說明，可用該模型來預測第2代棉鈴蟲百株卵量。

表3 棉鈴蟲百株卵量預測模型Table 3 Prediction model for number of boolworm eggs per 100 cotton plants

3 討 論

通過對第2代、第3代棉鈴蟲百株卵量與前期降水量指標的相關關系分析，并通過主成分分析法建立棉鈴蟲百株卵量預測模型，檢驗結果表明，第2代棉鈴蟲百株卵量預測方程信度較好，歷史符合率為70%，而第3代棉鈴蟲百株卵量預測方程信度較差，歷史符合率僅為40%。常曉娜等[12]研究表明，夏季各月的降水量是決定各代棉鈴蟲發生程度的關鍵決因素，此期間的降雨天數與第4代棉鈴蟲的發生關系最為密切；褚茗莉等[13]調查研究也表明，1992-1997年6月下旬至

圖1 棉鈴蟲百株卵量實測值與模型預報值對比Fig.1 Comparison of measured value and predicted value of the amountof Helicoverpa armigera per one hundred plants

8月中旬累計降雨量達200 mm左右時棉鈴蟲大發生，而1994-1996年同期累計降雨量平均都在400 mm以上時棉鈴蟲輕發生；李永昌等[14]研究認為，不但包括降雨在內的溫濕度顯著影響棉鈴蟲的發生，且風雨強度以及降雨量的分布對棉鈴蟲的發生程度也有明顯影響。本研究基本與上述研究結果一致。由于所選降水量指標均為前兆因子，第2代棉鈴蟲百株卵量預測模型可用于研究區域棉鈴蟲發生程度的監測、預警與研究。而從第3代棉鈴蟲百株卵量預測模型的檢驗效果來看，某些年份的預測結果還存在較大偏差，這可能是由于初選因子的局限性，由于資料有限，該模型有待進一步完善。

本研究雖然分析降水量與棉鈴蟲發生的關系，得出一些結論，氣候資料較為系統，代表性也較強，研究結果對棉花蟲害預測預報工作有一定的參考價值。但不足之處是沒有第1代棉鈴蟲百株卵量的資料，因此未能結合田間的蟲口基數預測第2代棉鈴蟲的百株卵量。

Reference：

[1]王厚振,華堯楠,牟吉元.棉鈴蟲預測預報與綜合治理[M].北京:中國農業出版社,1999:15,179-187.

WANG H ZH,HUA Y N,MU J Y.Prediction and Comprehensive Control for Cotton Bollworm [M].Beijing:Chinese Agriculture Press,1999:15,179-187(in Chinese).

[2]毛樹春.中國棉花可持續發展研究[M].北京:中國農業出版社,1999:149.

MAO SH CH.Studying for Sustained Development of Chinese Cotton[M].Beijing:Chinese Agriculture Press,1999:149(in Chinese).

[3]天門縣病蟲測報站,華中農學院農學系植保專業73級工農兵學員.棉鈴蟲的發生情況及發生期預測[J].華農科技,1976(1):144-149.

Tianmen County Surveying and Announcing Station on Plant Disease and Insect Pests,Workers-peasants-soldier Students Begin Studying from 1973 at Plant Protection Specialities in Agriculture Department of Huazhong Agriculture College.Forecasting occurring state and period of bollworm[J].JournalofHuazhongUniversity,1976(1):144-149(in Chinese).

[4]戈 峰,劉向輝,丁巖欽,等.華北棉區各代棉鈴蟲生命表及南北棉鈴蟲發生特征研究[J].應用生態學報,2003,4(2):242-245.

GE F ,LIU X H,DING Y Q,etal.Life-table ofHelicoverpaarmigerain northern China and characters of population development in southern and northern China[J].ChineseJournalofAppliedEcology,2003,4(2):242-245(in Chinese with English abstract).

[5]劉曉平,孫興全,吳樹芳.棉鈴蟲的發生規律及其綜合防治方法[J].安徽農學通報,2008,14(24):107-108.

LIU X P,SUN X Q,WU SH F,etal.Occurrence regularity and comprehensive control method of cotton bollworm [J].AnhuiAgriculturalScienceBulletin,2008,14(24):107-108(in Chinese).

[6]郭松景,李世民,卓喜牛,等.玉米田棉鈴蟲的發生危害特點及分布型研究[J].河南農業科學,2004(11):45-47.

GUO S J,LI SH M,ZHUO X N,etal.Studies on the occurrence and damage characteristics and its distribution type of cotton [J].JournalofHenanAgriculture,2004(11):45-47(in Chinese).

[7]郭培宗,苗連河,呂國輝,等.豫北地區棉鈴蟲發生規律研究[J].中國棉花,1998,25(5):15-17.

GUO P Z,MIAO L H,Lü G H,etal.Study on the occurrence regularity of cotton bollworm in the north of Henan [J].ChinaCotton,1998,25(5):15-17(in Chinese).

[8]王孝法.新疆墾區棉鈴蟲發生規律及防治措施[J].中國棉花,1999,26(7):43.

WANG X F.Occuring law and preventing and curing measure on boolworm in Xinjiang reclamation area[J].ChinaCotton,1999,26(7):43(in Chinese).

[9]柯惠新,黃京華,沈 浩.調查研究中的統計分析法[M].北京:北京廣播學院出版社,1992:476.

KE H X,HUANG J H,SHEN H.Statistical Analysis in Investigation and Research [M].Beijing:Beijing Broadcasting Institute Press,1992:476(in Chinese).

[10]肖炳光,張燕春,盧秀萍,等.烤煙品種主成分分析和聚類分析[J].種子,2000(2):27-29.

XIAO B G,ZHANG Y CH,LU X P,etal.The principal component analysis and cluster analysis of the flue-cured tobacco varieties[J].Seeds,2000(2):27-29 (in Chinese with English abstract).

[11]張玉榮,張秀華,周顯青,等.主成分分析法綜合評價大米的食味品質[J].河南工業大學學報,2008,29(5):1-5.

ZHANG Y R,ZHANG X H,ZHOU X Q,etal.The comprehensive evaluation of rice eating qualities by principal component analysis[J].JournalofHenanUniversityofTechnology,2008,29(5):1-5 (in Chinese with English abstract).

[12]常曉娜,高慧璟,陳法軍,等.環境濕度和降雨對昆蟲的影響[J].生態學雜志,2008,27(4):619-625.

CHANG X N,GAO H J,CHEN F J,etal.Effects of environmental moisture and precipitation on insects:A review[J].ChineseJournalofEcology,2008,27(4):619-625 (in Chinese with English abstract).

[13]褚茗莉,許國慶,田本志,等.遼寧省棉鈴蟲蟲情劃區的調查研究[J].遼寧農業科學,2000(1):15-18.

CHU M L,XU G Q,TIAN B ZH,etal.Contour survey of cotton bollworm in Liaoning provine[J].LiaoningAgriculturalScience,2000(1):15-18 (in Chinese with English abstract).

[14]李永昌,吳臨平,郭斌華,等.氣象要素對棉鈴蟲發生規律的影響[J].山西農業科學,1997,25(4):72-75.

LI Y CH,WU L P,GUO B H,etal.The effect of climite onHelicoverpaarmigeraoccurence[J].JournalofShanxiAgriculturalSciences,1997,25(4):72-75 (in Chinese with English abstract).

PrincipalComponentAnalysisForecastModelofHelicoverpaarmigeraBasedonPrecipitationIndex

LI Zhou and TANG Baoshan

(College of Agronomy, Northwest A&F University, Yangling Shaanxi 712100, China)

Helicoverpaarmigera(Hübner) is one of the main pests in cotton grown in China, meteorological factors have very important implications with its population dynamics. Using observational data ofH.armigeraand meteorological factors from 1961 to 1970, this paper adopted correlation coefficient method to select ten precipitation index and as forecasting index of the amount of second and third generationH.armigeraper one hundred plants, and using principal component analysis method to establish predictive models. The results show the model of second generation accuracy is 70%, and the third generation is 40%. Forecast based the model has strong timeliness and favorable effect, therefore, it can be used for monitoring and early warning of occurrence degree ofH.armigera.

Helicoverpaarmigera; Precipitation index; Principal component analysis; Prediction model

2017-08-07

2017-09-15

State 863 Projects (No. 2002AA241021).

LI Zhou, male，associate professor. Research area: cotton breeding and major pest controlling. E-mail: lz19580219@163.com

S435.132

A

1004-1389(2017)10-1554-05

日期：2017-10-18

網絡出版地址：http://kns.cnki.net/kcms/detail/61.1220.S.20171018.1733.040.html

2017-08-07

2017-09-15

國家863計劃(2002AA241021)。

李 胄，男，副教授，從事棉花育種及主要病蟲害防治研究。E-mail: lz19580219@163.com

唐保善，男，助理研究員，從事棉花植保與育種研究。E-mail: Tbs19620802@163.com

CorrespondingauthorTANG Baoshan, male, research assistant. Research area: cotton breeding and protecting. E-mail:Tbs19620802@163.com

(責任編輯：郭柏壽Responsibleeditor:GUOBaishou)