□文/鞠思祺
在氣候變化背景下,主要農業氣象災害發生的頻率和危害程度以及影響范圍均有增加趨勢。為了提升農業氣象災害風險評估的精度和效率,我們對主要農業氣象災害風險評估技術和其應用的研究進行了深入分析。基于貝葉斯網絡技術,我們將氣象參數、農業產量指標、農作物品質指標、農產品價格與氣象影響因子之間的關系建立了多元非線性回歸模型,并對全國主要農區的主要農業氣象災害風險進行了評估。研究結果表明,我國主要農區的主要農產品價格與氣象風險因子之間存在顯著相關性。
對于農業氣象災害風險評估的相關研究,通常以每年的氣象條件為基礎,再結合當地實際情況對其進行劃分。標準主要是根據受災面積、影響范圍等來判定風險等級。為了便于計算,本研究采用貝葉斯網絡技術進行農業氣象災害風險評估。其中,參數設置為:損失函數(i=1/2,i=0);氣象條件與作物產量指標之間的關系用二元非線性回歸模型表示;損失函數用線性模型表達。其中a為自變量,b為因變量。
選取影響因子是為了進行更好地統計分析,因此首先要對指標進行篩選,主要的指標包括:氣候因素、農業氣象災害、社會經濟因素、地理位置、自然地理環境以及水文情況。確定了指標的篩選方案,具體如下:作物產量、品質和價格對氣象風險的影響因子篩選分析。通過分析中國農作物氣象災害風險評估的相關研究成果,可以發現中國主要農作物氣象災害對農產品價格的影響顯著,并且農作物氣象災害和農產品價格之間呈現正相關關系。對于農作物氣象風險的預測預報,傳統的基于單一預測模型,如BP 神經網絡、灰色GM(1,1)模型等的預測結果都存在較大誤差。本文采用基于貝葉斯網絡分析(Bayesian Network Analysis)的預測方法,通過構建一系列多元非線性回歸模型,對中國主要農區的主要農產品價格與氣象風險因子之間進行建模和預測[1]。
農業氣象災害風險評估模型建立完成后,利用其模型,根據中國北方農業區劃數據和各地主要農產品產量信息,對全國范圍內的主要農作物產量和價格與氣象影響因子之間進行相關性分析,結果表明氣象影響因子對全國主要農區的農業產量影響較大。在不同氣象災害風險等級下,不同農產品的價格與氣象影響因子之間存在顯著相關關系。其中:玉米、小麥、水稻三種農產品的價格與氣象參數間相關關系最密切,其次是大豆和其他農產品;由于我國糧食產量占全國的比重較大,所以在農業氣象災害風險評估中首先考慮主要農作物價格與氣象參數之間的相關關系。為了減少因變量對自變量的過度依賴而產生的過擬合問題,可將自變量對數轉換為對數形式;如果有多個自變量時,則采用中心差分等方法進行處理。
除了上述方式,還可以通過建立作物產量模型來研究不同氣象條件對作物生長的影響,即確定影響各個品種產量的氣象因子及其閾值,然后根據上述指標設計產量模型,從而得到每個品種不同生長期產量的模擬值。將該數據作為產量預估參數,用于預測作物生長狀況。建立了生長模型之后,通過計算得到各個品種各發育期產量值及其分布情況。根據已有品種生長模型對不同氣象條件的敏感性、適應性和可調控性分析,建立品種不同發育期的種植制度,確定其最佳生產方案。將各品種種植制度與氣象因子相結合就得到區域氣候模式下作物生長的模擬結果。將不同生長期的結果按面積進行匯總分類并與作物最佳生長方案進行比較分析,得到區域氣候模式下最優生產方案。在氣候模式下按照區域氣候模型中農作物種植品種以及所需氣象因子閾值建立作物生長模型,得到每個品種農作物適宜的產量模擬值。將各品種適宜生產的農業氣象條件輸入至計算機中進行計算并生成相應的作物生長模型。將上述計算結果與實際數據相比較進行檢驗分析后得出最優種植方案,并確定未來該季節所需管理措施并制定生產計劃[2]。最后對各作物模型修正后用于模擬農業氣象災害發生時所造成農作物產量損失情況。
基于綜合作物生長模型,該系統由數據采集、數據管理與分析等子系統構成。作物產量數據管理:利用歷史氣象資料,結合農作物生長周期理論、作物發育模型與氣候模式等手段,建立了農作物生長模型。災害風險評估:以不同的風險等級為基準,在該等級的基礎上,計算各因素組合帶來的風險值,并以此作為對主要農業氣象災害造成損失程度判斷的標準。模型參數設置:該系統通過收集歷史氣象資料以及作物生長發育、生理生化過程及環境條件等資料,利用多元線性回歸與逐步回歸等數學方法,構建了不同農作物(棉花、水稻、玉米)的產量模擬模型以計算各主要農作物在不同產量水平下所對應的自然災害及農業氣象災害致災因子。模型驗證與模擬:利用建立好的作物生長模型對各個因素組合后所得到的產量水平進行模擬分析,并對模擬結果進行評估、檢驗與反饋改進,最終達到提高產量水平的目的[3]。
根據上述模型和技術體系,在基于作物生長模型的農業氣象災害風險評估系統的基礎上,設計了綜合作物產量風險評估系統。軟件系統主要由數據采集、數據處理、結果顯示三個部分組成,該軟件利用Visual Studio 2008 開發工具,在VB6.0 環境下開發和實現。基于VB6.0開發的綜合作物產量風險評估系統提供了查詢、統計、計算等數據處理功能,可以查詢到全國各地、各等級的氣象災害信息、歷年農作物種植面積以及當年棉花播種面積;計算出當年棉花單產;對不同等級農作物單產進行比較分析;輸入歷史氣象災害信息,計算相應區域的氣象災害影響程度;支持多用戶同時在線查詢等。