山東省玉米螟災害風險評估與區(qū)劃研究

張曉東, 左亞男, 季 民

（山東科技大學測繪與空間信息學院, 山東 青島 266000）

玉米螟（Pyrausta nubilalis）屬鱗翅目（Lepidoptera）螟蛾科（Pyralidae）, 又名玉米鉆心蟲, 是玉米上的重要害蟲, 也是世界性害蟲, 在山東省主要發(fā)生2～3 代, 玉米螟在玉米的各個生育時期都可以為害玉米植株的地上部分, 取食葉片、果穗、雄穗。其幼蟲孵化后群居于玉米心葉喇叭口處或嫩葉上取食, 葉片被幼蟲咬食后, 會降低其光合效率；雄穗被蛀, 常易折斷, 影響授粉；苞葉、花絲被蛀食, 會造成缺粒和秕粒；莖稈、穗柄、穗軸被蛀食后, 形成隧道, 破壞植株內(nèi)水分、養(yǎng)分的輸送, 使莖稈倒折率增加, 子粒產(chǎn)量下降。山東省是糧食大省, 也是玉米高產(chǎn)量省份, 玉米產(chǎn)量約占全國玉米總產(chǎn)量的10%, 為山東省主要種植作物。玉米產(chǎn)量每年因玉米螟都造成不同程度的損失, 因此, 對山東省玉米螟發(fā)生風險進行有效評估, 對玉米螟的防治有重要意義。

災害風險區(qū)劃是防災的重要工具, 在災害風險管理和防災減災方面具有重要的實際應用價值。賈改秀等[1］在探究甘肅省玉米病蟲害風險評估指標的基礎上, 構建了甘肅省玉米病蟲害風險區(qū)劃指標, 并進行了風險區(qū)劃和分區(qū)評述。袁福香等[2］整理了吉林省50 個縣市的主要農(nóng)作物產(chǎn)量面積資料和相關病蟲害資料, 根據(jù)災害發(fā)生的頻數(shù)確定病蟲害影響指數(shù), 分析各因子的影響程度, 并將各因子極差化處理, 計算出綜合風險防御指數(shù), 將吉林省主要農(nóng)作物病蟲害分成高、中、低風險防御區(qū)進行區(qū)劃評估。李陽等[3］利用彭陽縣的氣象、地形、土壤類型等數(shù)據(jù), 基于災害風險評估原理, 運用相關分析法、層次分析法, 結合GIS 技術對彭陽縣玉米干旱災害進行風險評估與區(qū)劃。Singh 等[4］在進行農(nóng)業(yè)區(qū)的風險區(qū)劃時, 從影響農(nóng)業(yè)區(qū)抗災能力的4 個方面考慮選取評價指標, 其多維度、多指標的討論過程為本研究評價指標的分析與選取提供了思路。

目前國內(nèi)外針對風險區(qū)劃[1-15］的研究已相對成熟, 但仍然缺少生物災害方向的區(qū)劃研究, 因此需對玉米螟災害的發(fā)生進行風險區(qū)劃, 以期為生物災害的防治及預測提供科學支撐。

對研究區(qū)域內(nèi)玉米螟的發(fā)生資料進行數(shù)據(jù)提取, 并將玉米螟發(fā)生程度分級國家農(nóng)業(yè)行業(yè)標準（NY/T 1611—2008）與山東省植物保護總站地方標準相結合, 提出了新的玉米螟發(fā)生程度分級準則, 并基于此準則進行發(fā)生程度定級, 統(tǒng)計發(fā)生程度頻數(shù)；以玉米螟發(fā)生頻數(shù)、玉米產(chǎn)量和玉米播種面積分別作為致災因子、承載體因子和孕災環(huán)境因子的影響指標, 基于參數(shù)估計法建立風險評估模型, 構建玉米螟災害風險評估指標體系, 進行了山東省玉米螟災害的綜合風險評估與區(qū)劃, 以期為生物災害的防治及預測提供科學支撐。

1 研究區(qū)與數(shù)據(jù)處理方法

1.1 研究區(qū)概況

山東省位于中國東部沿海、黃河下游, 地處34°25′—38°23′N, 114°36′—112°43′E, 東 西 長 約700 km, 南北寬約420 km, 陸地總面積約15.67萬km2, 是玉米生產(chǎn)和消費大省, 玉米面價格和總產(chǎn)量均居全國前列, 常年種植面積約2.7×106hm2。山東省屬溫暖帶季風氣候, 降水集中, 雨熱同季, 年平均氣溫為11～14 ℃, 年平均降水量為550～950 mm, 且集中于夏季, 光照時數(shù)年均為2 290～2 890 h, 自然條件非常適合玉米生長。

1.2 數(shù)據(jù)來源

所用數(shù)據(jù)包括山東省下屬16 個地市的玉米生產(chǎn)資料及玉米螟病蟲害發(fā)生資料。其中, 生產(chǎn)資料來源于《山東省統(tǒng)計年鑒》, 包括玉米產(chǎn)量及播種面積資料；病蟲害資料來自山東省農(nóng)業(yè)農(nóng)村廳發(fā)布的植保情況, 補充資料來源于各期刊文獻、新聞媒體及網(wǎng)站等發(fā)布的病蟲情報, 資料年限為2007—2019年。主要數(shù)據(jù)來源見表1。

表1 主要數(shù)據(jù)來源

1.3 數(shù)據(jù)預處理

1.3.1 發(fā)生程度分級準則 將收集到的數(shù)據(jù)資料進行整合, 由于數(shù)據(jù)資料包含被害株率、蟲株率、百株蟲量、百株卵量等多項指標, 且部分月份中已經(jīng)包含個別市的發(fā)生程度, 經(jīng)過比對, 發(fā)現(xiàn)其定級標準不符合中華人民共和國農(nóng)業(yè)行業(yè)標準——玉米螟測報技術規(guī)范（NY/T 1611—2008）, 因此, 在開展定級工作時, 通過查閱資料得知, ①山東省植物保護總站的玉米螟發(fā)生程度分級地方標準；②在農(nóng)作物蟲害方面, 被害株率、蟲株率、白穗率及卷葉率等被統(tǒng)稱為被害率, 其計算公式是統(tǒng)一的, 因此在定級工作中按照統(tǒng)一標準進行定級。為統(tǒng)一定級指標, 根據(jù)NY/T 1611—2008, 并與山東省植物保護總站農(nóng)作物病蟲發(fā)生程度分級地方標準相結合, 制定了以被害株率為主要劃分因子及百株蟲量為輔助因子的玉米螟發(fā)生程度分級準則（表2）。

表2 玉米螟發(fā)生程度分級

1.3.2 發(fā)生程度的定級 將上述數(shù)據(jù)根據(jù)表2 進行發(fā)生程度定級, 選用2009 年9 月的數(shù)據(jù)（資料來自《山東農(nóng)藥信息》）。首先, 將蟲株率按照被害株率的分級指標進行發(fā)生程度定級, 且當數(shù)據(jù)橫跨兩個級別時, 將較高級別定為其發(fā)生程度, 如2009 年9 月煙臺市玉米螟被害株率為10%～27%, 則其發(fā)生程度被定為2 級。其中沒有數(shù)據(jù)的地市定為未發(fā)生（棗莊市、臨沂市、日照市）, 其他指標按照表2 進行定級, 定級工作完成后按表3 所示統(tǒng)計玉米暝發(fā)生情況。

表3 2009 年9 月山東省各地市玉米螟發(fā)生情況統(tǒng)計

1.3.3 發(fā)生程度頻數(shù) 定級工作完成后, 分別統(tǒng)計發(fā)生程度為1、2、3 級的頻數(shù)V1、V2、V3, 得到2007—2019 年玉米螟發(fā)生頻數(shù)如表4 所示。

2 風險評估指標的確定及計算

2.1 風險評估指標的確定

風險評估指標體系的建立是進行病蟲害風險評估的關鍵環(huán)節(jié), 風險指標的選取應遵循科學性、實用性、系統(tǒng)性、人本性原則。依照上述原則, 通過查閱有關病蟲害研究文獻和歷史資料、危害程度、危害頻率, 結合病蟲害的特點, 提出病蟲害風險評估指標體系[1, 2, 4-9］, 包括致災因子、承載體因子及孕災環(huán)境因子。從這3 個方面闡述其對災害的影響, 選擇具有代表性且可以量化的影響因子作為風險區(qū)劃指標。

2.1.1 致災因子 致災因子是指造成糧食作物損失的各種病蟲害, 是重要的風險評估指標。本研究選用玉米螟2007—2019 年不同危害程度（輕發(fā)生、中發(fā)生、重發(fā)生）的發(fā)生頻數(shù)（表4）作為致災因子的影響指標。

表4 2007—2019 年山東省玉米螟發(fā)生頻數(shù)（V）（單位：次）

2.1.2 承載體因子 承載體因子是指當病蟲害發(fā)生時農(nóng)作物受病蟲害影響的因素, 是致災因子作用的實體。本研究中, 受玉米螟災害影響的承載體是玉米的產(chǎn)量。

2.1.3 孕災環(huán)境因子 孕災環(huán)境是指影響各種病蟲害生長、繁殖的外界環(huán)境條件。本研究所考慮的影響指標包括作物的種植面積, 還有導致病蟲害發(fā)生的氣象條件, 氣象條件的適宜與否致使病蟲害發(fā)生的程度不同, 此指標已包含在致災因子的影響中, 不再重復計算, 因此本研究只選用玉米播種面積作為病蟲害風險評估指標。

綜合以上分析, 所建立的玉米螟災害風險評估指標體系如圖1 所示。

圖1 玉米螟風險評估指標體系

2.2 風險評估指標的計算

2.2.1 玉米螟風險指數(shù) 病蟲害影響指數(shù)是衡量病蟲害對農(nóng)作物產(chǎn)量損失風險大小的指標。利用山東省16 個地市2007—2019 年玉米螟危害資料, 根據(jù)玉米螟歷年不同發(fā)生程度（1、2、3 級）的頻數(shù)（V1、V2、V3）（表3）, 利用加權平均得出玉米螟對產(chǎn)量的影響程度, 即病蟲害影響指數(shù), 也稱作玉米螟風險指數(shù)影響因子。重發(fā)生的產(chǎn)量損失大, 則權重系數(shù)大, 輕發(fā)生的產(chǎn)量損失小, 則權重系數(shù)小。風險指數(shù)（R）的計算公式如下。

式中,R是玉米螟風險指數(shù)；i是各地市（下同）。

根據(jù)表3 及式（1）計算出R, 根據(jù)各地市R, 利用反距離權重法（IDW）、掩膜提取及自然斷點分級法[11］將玉米螟影響指數(shù)分5 個等級進行區(qū)劃評估, 得到山東省玉米螟影響指數(shù)區(qū)劃（圖2）。由圖2 可知, 玉米螟災害發(fā)生的主要風險地集中在山東省的西部地區(qū), 該地區(qū)是山東省玉米主要種植區(qū)和高產(chǎn)區(qū), 同時也是玉米螟發(fā)生的高風險區(qū)。其中, 菏澤市、濱州市風險指數(shù)高, 為玉米螟發(fā)生的高風險地區(qū)；濟南市中西部為較高風險區(qū)。

圖2 玉米螟影響指數(shù)區(qū)劃

2.2.2 玉米產(chǎn)量影響因子 玉米螟發(fā)生區(qū)域產(chǎn)量越高, 損失越大, 則玉米螟危害風險越大, 因此本研究選擇玉米產(chǎn)量作為承載體因子的影響指標。產(chǎn)量影響因子的取值為各地市產(chǎn)量占全省總產(chǎn)量的比重, 產(chǎn)量值為近10 年平均值。玉米產(chǎn)量影響因子（K）的計算公式如下。

式中,Ky為各地市的玉米總產(chǎn)量；∑Ky為全省的玉米總產(chǎn)量。

根據(jù)式（2）計算出K, 并繪制玉米產(chǎn)量影響區(qū)劃（圖3）。由圖3可知, 山東省西部地區(qū)及濰坊市中部為玉米高產(chǎn)區(qū), 也是玉米螟發(fā)生高風險區(qū), 其中德州市產(chǎn)量最高；次高地區(qū)依次為菏澤市、聊城市及濰坊市。

圖3 玉米產(chǎn)量影響區(qū)劃

2.2.3 玉米播種面積影響因子 播種面積也是影響玉米螟發(fā)生的主要影響因子, 播種面積越大, 病蟲害越容易傳播和流行, 則受害風險性越大, 因此本研究選擇播種面積作為孕災環(huán)境因子的影響指標。播種面積計算公式為山東省各地市總播種面積占全省總播種面積的比重, 所用播種面積為近10 年的平均值。玉米播種面積影響因子（S）的計算公式如下。

式中,Sy為各地市的玉米播種總面積；∑Sy為全省的玉米播種總面積。

根據(jù)式（3）計算出S, 以S代表各地市的孕災環(huán)境因子繪制玉米播種面積影響區(qū)劃（圖4）。由圖4可知, 山東省玉米主要種植區(qū)分布在該省西部地區(qū)及濰坊市中部地區(qū), 也是山東省玉米高產(chǎn)區(qū), 同時也是玉米螟可能發(fā)生的高風險地區(qū)。其中高風險區(qū)位于德州市西北部及菏澤市中西部地區(qū)；聊城市、濰坊市為較高風險區(qū)。

圖4 玉米播種面積影響區(qū)劃

3 玉米螟綜合風險指數(shù)

為綜合評價山東省各地市玉米螟發(fā)生風險狀況, 為山東省生物災害的防治提供科學依據(jù), 對山東省玉米螟災害發(fā)生頻數(shù)進行綜合風險區(qū)劃[2, 4, 5, 7-9］。

3.1 數(shù)據(jù)歸一化

由于影響玉米螟災害風險的指標之間存在量級和量綱上的差異, 為了使不同的數(shù)據(jù)層能夠計算和比較, 需要將數(shù)據(jù)標準化處理。本研究采用極差化標準方法將單元數(shù)據(jù)轉(zhuǎn)換至[0, 1］范圍內(nèi)的值。設研究區(qū)域內(nèi)有m個評價指標, 有n個評價單元, 則有原始評價矩陣如下。

用極差標準化方法進行如下?lián)Q算。

式中,Xij是第i個指標第j個單元對應的原始數(shù)據(jù)；max｛Xij｝和min｛Xij｝分別是第i個指標下所有單元數(shù)據(jù)中的最大值和最小值；X′ij是標準化的單元數(shù)據(jù)。

3.2 綜合風險指數(shù)的計算

將極差化后的各風險指標相加, 得到包含原來幾個風險要素綜合影響的新指標, 稱為綜合風險指數(shù), 進而繪制綜合風險區(qū)劃, 進行區(qū)劃分析。綜合風險指數(shù)（G）的計算公式如下。

式中,R′i、K′i、S′i分別為致災因子、承載體因子、孕災環(huán)境因子極差化后的指標。

3.3 結果與分析

利用反距離權重法（IDW）、掩膜提取及自然斷點分級法將綜合風險指數(shù)分5 個等級進行區(qū)劃評估, 最后得出山東省玉米螟災害綜合風險區(qū)劃（圖5）。由圖5 可知, 綜合風險指數(shù)最大出現(xiàn)在菏澤市, 為0.897, 該地產(chǎn)量高, 種植面積大, 玉米螟發(fā)生頻率高, 是玉米螟為害造成損失的高風險區(qū)；最小出現(xiàn)在日照市, 是玉米螟發(fā)生最低風險區(qū)；其中德州市、濰坊市、濱州市都在0.53 以上, 為較高風險區(qū), 主要集中在山東省中西部主要產(chǎn)量區(qū)；日照市、威海市北部及東營市、淄博市、棗莊市等地都在0.20 以下, 為玉米螟發(fā)生低風險區(qū), 不易受到玉米螟為害；濟南市、濟寧市、臨沂市、泰安市、青島市都處于0.30～0.53, 為玉米螟發(fā)生中風險區(qū), 條件適宜就容易發(fā)生玉米螟為害, 或有局部暴發(fā)玉米螟的風險；全省其他地區(qū)都處于0.20～0.30, 是玉米螟發(fā)生較低風險區(qū), 比較不容易受到玉米螟為害（表5）。

圖5 山東省玉米螟綜合風險區(qū)劃

4 小結

1）對研究區(qū)域內(nèi)玉米螟災害的發(fā)生資料進行了搜集與提取, 并將中華人民共和國農(nóng)業(yè)行業(yè)標準——玉米螟測報技術規(guī)范（NY/T 1611—2008）與山東省植物保護總站的玉米螟發(fā)生程度分級地方標準相結合, 提出了新的玉米螟發(fā)生程度分級標準, 并在此標準的基礎上進行數(shù)據(jù)預處理。

2）綜合考慮孕災環(huán)境因子、致災因子及承載體因子三方面, 建立玉米螟災害風險評估指標體系, 對山東省玉米螟災害的發(fā)生風險展開分析與評價, 具體選取了玉米螟病蟲害影響指數(shù)、玉米產(chǎn)量、玉米播種面積3 項影響指標, 并根據(jù)影響指標將收集到的數(shù)據(jù)通過參數(shù)估計法風險評估模型和反距離權重法、掩膜提取法、自然斷點分級法進行處理, 繪制出研究區(qū)域內(nèi)各影響因子的區(qū)劃。

3）為綜合評價山東省各地玉米螟災害發(fā)生的風險狀況, 對山東省玉米螟進行綜合風險區(qū)劃。為消除不同量綱的影響, 將每一指標極差化處理, 然后將極差化后的風險指標相加, 得到綜合風險指數(shù), 根據(jù)綜合風險指數(shù)繪制出山東省玉米螟災害綜合風險區(qū)劃。

4）綜合風險區(qū)劃結果表明, 山東省玉米螟災害發(fā)生的重點風險區(qū)集中在該省中西部地區(qū)。其中, 菏澤市屬于玉米螟發(fā)生的高風險區(qū)；德州市、濰坊市、濱州市及聊城市屬于較高風險區(qū)；濟南市、濟寧市、臨沂市、泰安市及青島市為中風險區(qū)；煙臺市為玉米螟發(fā)生較低風險區(qū)；日照市、威海市、東營市、淄博市及棗莊市屬于低風險區(qū)。區(qū)劃結果對山東省防災減災事業(yè)及玉米螟災害的防治具有重要的參考價值。