于小兵，陳 虹，吉中會，李陳靚

(1. 南京信息工程大學(xué) 氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，江蘇 南京210044；2. 南京信息工程大學(xué) 管理工程學(xué)院，江蘇 南京210044)

作為產(chǎn)糧大省，農(nóng)業(yè)在山東經(jīng)濟發(fā)展中占有舉足輕重的地位。據(jù)統(tǒng)計，2016年山東省發(fā)生了包括低溫冷凍、風(fēng)雹、洪澇、干旱在內(nèi)的多種自然災(zāi)害。其中，農(nóng)作物受災(zāi)面積達552.24 khm2，成災(zāi)面積為300.13 khm2，絕收面積達到62.12 khm2。直接經(jīng)濟損失70.39億元，農(nóng)業(yè)損失63.68億元，占比超過90%。可見，研究山東省農(nóng)業(yè)氣象災(zāi)害對農(nóng)作物產(chǎn)量的影響具有重要現(xiàn)實意義。

國外學(xué)者針對氣象災(zāi)害對農(nóng)作物產(chǎn)量的影響開展了大量研究[1-3]。國內(nèi)學(xué)者也對此進行了研究[4-7]。余會康等[8]分析了不同階段氣象變化對福建省糧食產(chǎn)量的影響。楊月鋒等[9]發(fā)現(xiàn)農(nóng)業(yè)現(xiàn)代化投入水平狀況、生產(chǎn)資料及糧食價格政策是福建省糧食產(chǎn)量變化的主要驅(qū)動因素。部分學(xué)者利用灰色關(guān)聯(lián)法對農(nóng)業(yè)氣象災(zāi)害評估進行研究[10-13]。溫丹蘋等[14]構(gòu)建了一種基于灰色關(guān)聯(lián)分析法的數(shù)理模型用以分析大豆產(chǎn)量和農(nóng)業(yè)氣象災(zāi)害變化特征。馬建勇等[15]采用灰色關(guān)聯(lián)法對東北地區(qū)1971-2009年間農(nóng)業(yè)氣象災(zāi)害的變化趨勢及其對糧食產(chǎn)量的影響開展了研究。

目前，已有大量關(guān)于農(nóng)業(yè)氣象災(zāi)害和糧食產(chǎn)量方面的研究，但主要集中討論全國或省(市、自治區(qū))氣象災(zāi)害的歷史變化趨勢及其對糧食總產(chǎn)量的影響，鮮有學(xué)者以相對災(zāi)損量指標為切入點對農(nóng)業(yè)氣象災(zāi)害與作物產(chǎn)量關(guān)系進行研究。本研究依據(jù)中國種植業(yè)管理司1950-2016年的農(nóng)業(yè)統(tǒng)計數(shù)據(jù)，利用灰色關(guān)聯(lián)分析法對農(nóng)作物總產(chǎn)量與氣象災(zāi)害受災(zāi)范圍及其單產(chǎn)與農(nóng)業(yè)氣象災(zāi)害強度之間的關(guān)系進行深入研究。

1 資料與研究方法

1.1 資料來源

本研究選取山東省農(nóng)作物的單產(chǎn)資料，旱災(zāi)、水災(zāi)、風(fēng)雹、低溫、臺風(fēng)的受災(zāi)面積、成災(zāi)面積和絕收面積。統(tǒng)計資料來源于《中國統(tǒng)計年鑒》及《中國農(nóng)業(yè)統(tǒng)計資料》和農(nóng)業(yè)部種植業(yè)管理司災(zāi)情數(shù)據(jù)庫，一共采用了67年的統(tǒng)計數(shù)據(jù)(1950-2016年)。

1.2 研究方法

1.2.1 受災(zāi)比率及相對災(zāi)損量

受災(zāi)比率是農(nóng)作物受災(zāi)面積占該種農(nóng)作物播種面積的百分比，表示受氣象災(zāi)害影響的農(nóng)作物受災(zāi)范圍的大小; 災(zāi)損量是農(nóng)業(yè)氣象災(zāi)害致使作物的減產(chǎn)量，表示農(nóng)業(yè)氣象災(zāi)害對作物產(chǎn)量影響的強弱。因此，可采用受災(zāi)比率和災(zāi)損量來表征氣象災(zāi)害對某種農(nóng)作物受災(zāi)范圍和災(zāi)害強度的影響。受災(zāi)面積、成災(zāi)面積和絕收面積是災(zāi)情統(tǒng)計中最常用的指標。其中，因災(zāi)害致使減產(chǎn)10% 以上的種植面積為受災(zāi)面積; 因災(zāi)害致使減產(chǎn) 30% 以上的種植面積為成災(zāi)面積; 因災(zāi)害致使減產(chǎn)80% 以上的種植面積為絕收面積。減產(chǎn)分成法可對各類農(nóng)業(yè)氣象災(zāi)害造成的糧食減產(chǎn)量進行估算，故采用減產(chǎn)分成法對各類農(nóng)業(yè)氣象災(zāi)害造成的農(nóng)作物災(zāi)損量進行估計[14]，以1950-2016年發(fā)生的農(nóng)業(yè)氣象災(zāi)害為例，則災(zāi)損量可表達為：

Yij=(0.1×Aij+0.2×Bij+0.5×Cij)×Dj,i=1,2,3,4,5;j=1,2,…,67。

(1)

式中:i=1,2,3,4,5分別表示洪澇災(zāi)、干旱災(zāi)、風(fēng)雹災(zāi)、低溫災(zāi)和臺風(fēng)災(zāi)。j表示所研究的種植周期，j=1表示1950年；Yij為第i類農(nóng)業(yè)氣象災(zāi)害在第j年的造成的農(nóng)作物災(zāi)損量(kg / hm2)。Aij，Bij和Cij分別表示第i類農(nóng)業(yè)氣象災(zāi)害在第j年某種農(nóng)作物的受災(zāi)面積、成災(zāi)面積和絕收面積(hm2)。Dj表示第j年某種農(nóng)作物的單位面積產(chǎn)量(kg/ hm2)。

相對災(zāi)損量 表示不同年份不同農(nóng)業(yè)氣象災(zāi)害造成農(nóng)作物減產(chǎn)的相對損失程度，它避免了時空分布特征對農(nóng)作物產(chǎn)量的影響，具有可比較性。 表示第i類農(nóng)業(yè)氣象災(zāi)害在第j年的農(nóng)作物災(zāi)損量與該年農(nóng)作物總減產(chǎn)量的比值，即：

(2)

1.2.2 灰色關(guān)聯(lián)分析

灰色關(guān)聯(lián)分析是衡量因素之間關(guān)聯(lián)程度強弱的多因素統(tǒng)計分析方法。在系統(tǒng)發(fā)展過程中，關(guān)聯(lián)度越強，說明因素間關(guān)系越緊密; 反之，則越不緊密?；貧w分析、方差分析和主成分分析等方法要求樣本具有某類典型概率分布特征，而灰色關(guān)聯(lián)分析則不需要。鑒于本文所需數(shù)據(jù)不滿足典型概率分布特征，故在此運用灰色關(guān)聯(lián)法進行分析，具體步驟如下:

(1) 確定參考數(shù)列和比較數(shù)列。如以1950-2016年某種農(nóng)作物單產(chǎn)的時間序列數(shù)據(jù)為參考數(shù)列，以1950-2016年5類農(nóng)業(yè)氣象災(zāi)害相對災(zāi)損量的時間序列數(shù)據(jù)為比較數(shù)列。由于農(nóng)業(yè)氣象災(zāi)害相對災(zāi)損量與農(nóng)作物單產(chǎn)之間存在負相關(guān)關(guān)系，故將其相對災(zāi)損量數(shù)據(jù)倒數(shù)化后進行計算。假設(shè)原始數(shù)據(jù)為Xu=(Xu(1),Xu(2),…Xu(n))(u=0,1,2,3,4,5；n=1,2,…,15)。其中，參考數(shù)列為X0某種農(nóng)作物單產(chǎn)，比較數(shù)列分別為X1洪澇災(zāi)相對災(zāi)損量的倒數(shù)化，X2旱災(zāi)相對災(zāi)損量的倒數(shù)化，X3風(fēng)雹災(zāi)相對災(zāi)損量的倒數(shù)化，X4低溫災(zāi)相對災(zāi)損量的倒數(shù)化和X5臺風(fēng)災(zāi)相對災(zāi)損量的倒數(shù)化。

(2) 處理原始數(shù)據(jù)Xu。為消除原始數(shù)據(jù)在量綱和數(shù)量級上的差異，采用均值化方法處理原始數(shù)據(jù)，即：

(3)

(4)

(4) 計算關(guān)聯(lián)系數(shù)ρ0u。Δ(min)和Δ(max)分別為全部u個比較數(shù)列在各期絕對差值中的最小者和最大者。ρ0u為比較數(shù)列與參考數(shù)列在t時期的關(guān)聯(lián)系數(shù)。

(5)

式中：ε為分辨系數(shù)，用以削弱Δ(max)過大而致使關(guān)聯(lián)系數(shù)失真的影響，文中ε= 0.2。

(5) 關(guān)聯(lián)度ρ0u的計算與比較。

(6)

式中：ρ0u表示比較數(shù)列u與參考數(shù)列的關(guān)聯(lián)度。關(guān)聯(lián)度越大，兩者關(guān)系越緊密。根據(jù)關(guān)聯(lián)度的大小進行排序，確定比較數(shù)列對參考數(shù)列的相關(guān)程度。

2 山東省主要農(nóng)作物與農(nóng)業(yè)氣象災(zāi)害分析

2.1 山東省小麥產(chǎn)量與農(nóng)業(yè)氣象災(zāi)害分析

1950-2016年，山東省小麥單產(chǎn)和小麥總產(chǎn)整體均呈現(xiàn)上升的趨勢(圖1)。其中: 第1階段(1950-1959年) ，山東省小麥單產(chǎn)和總產(chǎn)均較為平穩(wěn)；第2階段(1960-1969年) ，山東省小麥單產(chǎn)和總產(chǎn)均呈先減少后增加的趨勢，并在1969年達到歷史最低，小麥單產(chǎn)為25.81 kg/hm2，總產(chǎn)為131.5萬t；第3階段(1970-1979年)、第4階段(1980-1989年)、第5階段(1990-1999年)和第6階段(2000-2009年) 四個階段均呈現(xiàn)波動上升趨勢。第7階段(2010-2016年)，小麥單產(chǎn)和總產(chǎn)均為平穩(wěn)上升趨勢，均在2015年達到歷史最高。

圖1 山東省1950-2016年小麥單產(chǎn)及總產(chǎn)量變化

在受災(zāi)范圍的關(guān)聯(lián)分析中，以1950-2016年山東省小麥總產(chǎn)量為參考數(shù)列，以 1950-2016年山東省小麥5類農(nóng)業(yè)氣象災(zāi)害受災(zāi)比率為比較數(shù)列，計算山東省小麥總產(chǎn)量和5類農(nóng)業(yè)氣象災(zāi)害受災(zāi)比率的關(guān)聯(lián)度(表1)。

結(jié)果表明，第1、2、3階段的臺風(fēng)受災(zāi)范圍與山東省小麥總產(chǎn)量的關(guān)聯(lián)度最高，其次是風(fēng)雹，臺風(fēng)是山東省1950-1979年影響小麥總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。第4階段各類氣象災(zāi)害的受災(zāi)范圍與小麥總產(chǎn)的關(guān)聯(lián)度由高到低排序為:旱災(zāi)、洪澇、風(fēng)雹、低溫和臺風(fēng)，旱災(zāi)是山東省1980-1989年影響小麥總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。第5、6階段洪澇、旱災(zāi)、風(fēng)雹、低溫和臺風(fēng)的受災(zāi)范圍與山東省小麥總產(chǎn)的關(guān)聯(lián)度由高到低排序均為: 風(fēng)雹、洪澇、旱災(zāi)、低溫和臺風(fēng)。第7階段中，洪澇是山東2010-2016年影響小麥總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。

總體而言，風(fēng)雹、旱災(zāi)和洪澇對山東省小麥產(chǎn)量危害較大。風(fēng)雹在7個階段中的關(guān)聯(lián)度均位列前茅，是影響山東省小麥總產(chǎn)的關(guān)鍵農(nóng)業(yè)氣象災(zāi)害。在研究周期內(nèi)，風(fēng)雹、旱災(zāi)、洪澇三者的關(guān)聯(lián)度分別為0.907、0.882和0.878，故從受災(zāi)范圍角度看，風(fēng)雹、旱災(zāi)和洪澇對山東省小麥生產(chǎn)具有較大影響。

在災(zāi)害強度分析中，以山東省小麥5類農(nóng)業(yè)氣象災(zāi)害的相對災(zāi)損量為比較數(shù)列。以1950-2016年山東省小麥單產(chǎn)為參考數(shù)列，計算山東省小麥單產(chǎn)與5類農(nóng)業(yè)氣象災(zāi)害的關(guān)聯(lián)度(表2)。

結(jié)果表明，第1階段洪澇災(zāi)、旱災(zāi)、風(fēng)雹、低溫和臺風(fēng)的災(zāi)害強度與山東省小麥單產(chǎn)的關(guān)聯(lián)度由高到低為: 低溫、臺風(fēng)、洪澇、旱災(zāi)和風(fēng)雹; 第2階段5類農(nóng)業(yè)氣象災(zāi)害強度與小麥單產(chǎn)關(guān)聯(lián)度由高到低為: 洪澇、臺風(fēng)、低溫、旱災(zāi)和風(fēng)雹；第3、4、5、6、7階段5類農(nóng)業(yè)氣象災(zāi)害強度與山東省小麥單產(chǎn)的關(guān)聯(lián)度排序中，關(guān)聯(lián)度最高的均為旱災(zāi)。在總的研究周期內(nèi)，旱災(zāi)對山東省小麥單產(chǎn)的影響程度最大，是第三階段至第七階段的關(guān)鍵農(nóng)業(yè)氣象災(zāi)害，洪澇和風(fēng)雹次之，而臺風(fēng)和低溫在該階段對小麥單產(chǎn)影響較小。

表1 小麥總產(chǎn)與受災(zāi)比率的灰色關(guān)聯(lián)度分析

表2 小麥單產(chǎn)與相對災(zāi)損量的灰色關(guān)聯(lián)度分析

表3 玉米總產(chǎn)與受災(zāi)比率的灰色關(guān)聯(lián)分析

表4 玉米單產(chǎn)與相對災(zāi)損量的灰色關(guān)聯(lián)分析

2.2 山東省玉米產(chǎn)量與農(nóng)業(yè)氣象災(zāi)害分析

1950-2016年，山東省玉米單產(chǎn)和玉米總產(chǎn)整體上均呈現(xiàn)增加的趨勢(圖2)。第1階段，玉米單產(chǎn)和總產(chǎn)均較為平穩(wěn)，無明顯波動。第2階段的山東省玉米單產(chǎn)和總產(chǎn)在1960年達到歷史最低，玉米單產(chǎn)52.43 kg/hm2，玉米總產(chǎn)為74.9萬t。第3階段至第6階段這四個階段均呈現(xiàn)波動上升，其中第5階段和第6階段波動較為劇烈。第7階段，山東省玉米單產(chǎn)在2012年達到歷史最高440.57 kg/hm2，而總產(chǎn)呈現(xiàn)緩慢上升趨勢，在2016年達到歷史最高2 064.95萬t。1960年山東省玉米單產(chǎn)約為研究期內(nèi)玉米單產(chǎn)峰值的11.9%，而1960年山東省玉米總產(chǎn)量約為研究期內(nèi)玉米總產(chǎn)量峰值的3.63%。

圖2 1950-2016年玉米單產(chǎn)及總產(chǎn)量變化

在受災(zāi)范圍的關(guān)聯(lián)分析中，以1950-2016年山東省玉米總產(chǎn)量為參考數(shù)列，以 1950-2016年山東省玉米5類農(nóng)業(yè)氣象災(zāi)害受災(zāi)比率為比較數(shù)列，計算山東省玉米總產(chǎn)量和5類農(nóng)業(yè)氣象災(zāi)害受災(zāi)比率的關(guān)聯(lián)度(表3)。結(jié)果表明，第1、2階段臺風(fēng)災(zāi)害與山東省玉米總產(chǎn)的關(guān)聯(lián)度最高，其次是風(fēng)雹。臺風(fēng)是山東省1950-1979年影響玉米總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。第3、4階段受災(zāi)范圍與山東省玉米總產(chǎn)的關(guān)聯(lián)度最高分別為洪澇和旱災(zāi)。第5、6階段風(fēng)雹災(zāi)的受災(zāi)范圍與山東省玉米總產(chǎn)的關(guān)聯(lián)度最高，風(fēng)雹是山東省1990-2009年影響玉米總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。第7階段關(guān)聯(lián)度由高到低排序為：低溫、洪澇、風(fēng)雹、旱災(zāi)和臺風(fēng)。從受災(zāi)范圍看，在研究周期內(nèi)，風(fēng)雹是關(guān)鍵農(nóng)業(yè)氣象災(zāi)害，低溫和旱災(zāi)次之，三者關(guān)聯(lián)度分別為0.881、0.864和0.864。

在災(zāi)害強度分析中，以山東省玉米5類農(nóng)業(yè)氣象災(zāi)害的相對災(zāi)損量為比較數(shù)列。以1950-2016年山東省玉米單產(chǎn)為參考數(shù)列，計算山東省玉米單產(chǎn)與5類農(nóng)業(yè)氣象災(zāi)害的關(guān)聯(lián)度(表4)。

結(jié)果表明，第1階段洪澇、旱災(zāi)、風(fēng)雹、低溫和臺風(fēng)的災(zāi)害強度與山東省玉米單產(chǎn)的關(guān)聯(lián)度由高到低排序為: 低溫、臺風(fēng)、洪澇、旱災(zāi)和風(fēng)雹; 第2階段5類農(nóng)業(yè)氣象災(zāi)害強度與山東省玉米單產(chǎn)的關(guān)聯(lián)度由高到低排序為: 洪澇、臺風(fēng)、低溫、旱災(zāi)和風(fēng)雹；第3～7階段5類農(nóng)業(yè)氣象災(zāi)害強度與山東省玉米單產(chǎn)的關(guān)聯(lián)度排序中，旱災(zāi)均為關(guān)聯(lián)度最高者，洪澇和風(fēng)雹次之。在研究周期內(nèi)，5類農(nóng)業(yè)氣象災(zāi)害強度與山東省玉米單產(chǎn)的關(guān)聯(lián)度排序由高到低為: 旱災(zāi)、洪澇、風(fēng)雹、臺風(fēng)和低溫。從災(zāi)害強度看，總體來說，旱災(zāi)對山東省玉米單產(chǎn)的危害最大，洪澇和風(fēng)雹次之，而臺風(fēng)和低溫對山東省玉米單產(chǎn)影響較小。

2.3 山東省棉花產(chǎn)量與農(nóng)業(yè)氣象災(zāi)害分析

1950-2016年，山東省棉花單產(chǎn)和棉花總產(chǎn)整體上均呈現(xiàn)上升的趨勢(圖3)。其中: 第1階段、第2階段和第3階段，山東省棉花單產(chǎn)和總產(chǎn)均較為平穩(wěn)。第4階段和第5階段均先劇烈上升后劇烈下降，且棉花總產(chǎn)量在1984年達到歷史最高172.5萬t后開始劇烈下降，棉花單產(chǎn)在劇烈下降后有所上升；第6階段和第7階段，山東省棉花單產(chǎn)呈波動上升趨勢，在2016年達歷史最高78.57 kg/hm2，而棉花總產(chǎn)則先上升后下降。1961年山東省棉花單產(chǎn)約為研究期內(nèi)棉花單產(chǎn)峰值的6.75%，而1961年山東省棉花總產(chǎn)量約為研究期內(nèi)棉花總產(chǎn)量峰值的1.86% 。

圖3 1950-2016年棉花單產(chǎn)和總產(chǎn)量的變化

表6 棉花單產(chǎn)與相對災(zāi)損量的灰色關(guān)聯(lián)分析

在受災(zāi)范圍的關(guān)聯(lián)分析中，以1950-2016年山東省棉花總產(chǎn)量為參考數(shù)列，以1950-2016年山東省棉花5類農(nóng)業(yè)氣象災(zāi)害受災(zāi)比率為比較數(shù)列，計算山東省棉花總產(chǎn)量和5類農(nóng)業(yè)氣象災(zāi)害受災(zāi)比率的關(guān)聯(lián)度(表5)。結(jié)果表明，第1、2、3階段臺風(fēng)的受災(zāi)范圍與山東省棉花總產(chǎn)的關(guān)聯(lián)度最高，臺風(fēng)是山東省1950-1979年影響棉花總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。第4階段農(nóng)業(yè)氣象災(zāi)害的受災(zāi)范圍與棉花總產(chǎn)的關(guān)聯(lián)度由高到低排序為:旱災(zāi)、風(fēng)雹、洪澇、低溫和臺風(fēng)，旱災(zāi)是山東省1980-1989年影響棉花總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。第5～7階段的關(guān)聯(lián)度由高到低排序均為: 風(fēng)雹、洪澇、旱災(zāi)、低溫和臺風(fēng)，風(fēng)雹是山東省1990-2016年影響棉花總產(chǎn)的主要農(nóng)業(yè)氣象災(zāi)害。但總體來看，風(fēng)雹、旱災(zāi)和洪澇是影響山東棉花作物總產(chǎn)的關(guān)鍵農(nóng)業(yè)氣象災(zāi)害。

從受災(zāi)范圍看，風(fēng)雹、旱災(zāi)和洪澇對山東省棉花產(chǎn)量危害較大。風(fēng)雹是關(guān)鍵農(nóng)業(yè)氣象災(zāi)害，旱災(zāi)和洪澇次之，三者關(guān)聯(lián)度分別為0.900、0.875和0.869。在災(zāi)害強度分析中，以山東省棉花5類農(nóng)業(yè)氣象災(zāi)害的相對災(zāi)損量為比較數(shù)列。以1950-2016年山東省棉花單產(chǎn)為參考數(shù)列，計算山東省棉花單產(chǎn)與5類農(nóng)業(yè)氣象災(zāi)害的關(guān)聯(lián)度(表6)。結(jié)果表明，第1階段洪澇、旱災(zāi)、風(fēng)雹、低溫和臺風(fēng)的災(zāi)害強度與山東省棉花單產(chǎn)的關(guān)聯(lián)度由高到低排序為: 低溫、洪澇、臺風(fēng)、旱災(zāi)和風(fēng)雹; 第2階段的關(guān)聯(lián)度由高到低排序為: 洪澇、臺風(fēng)、低溫、旱災(zāi)和風(fēng)雹；第3～6階段5類農(nóng)業(yè)氣象災(zāi)害強度與山東省棉花單產(chǎn)的關(guān)聯(lián)度排序中，旱災(zāi)均為關(guān)聯(lián)度最高者，洪澇和風(fēng)雹次之。在研究周期內(nèi)，5類農(nóng)業(yè)氣象災(zāi)害強度與山東省棉花單產(chǎn)的關(guān)聯(lián)度排序由高到低為: 旱災(zāi)、洪澇災(zāi)、風(fēng)雹災(zāi)、臺風(fēng)災(zāi)和低溫災(zāi)。

3 結(jié)論

本文采用灰色關(guān)聯(lián)法對山東省主要農(nóng)作物產(chǎn)量與氣象災(zāi)害受災(zāi)范圍及災(zāi)害強度間的關(guān)系進行研究。從整個研究周期來看，在受災(zāi)范圍方面，風(fēng)雹、旱災(zāi)是影響山東省三種農(nóng)作物總產(chǎn)量的主要農(nóng)業(yè)氣象災(zāi)害。從受災(zāi)強度看，農(nóng)作物單產(chǎn)主要受到旱災(zāi)和洪澇的影響。旱災(zāi)無論從受災(zāi)范圍還是災(zāi)害強度上對農(nóng)作物的生產(chǎn)影響都較大。低溫和臺風(fēng)則無論從受災(zāi)范圍還是受災(zāi)強度方面，對農(nóng)作物的生產(chǎn)影響都較小。因此，種植具有耐干旱特質(zhì)的農(nóng)作物品種對山東省農(nóng)戶而言是一個較好的選擇。