黑龍江省玉米低溫冷害天氣指數保險設計

余 蘭 ，鐘化然

（1.海倫市氣象局，黑龍江 海倫 152300；2.安達市氣象局，黑龍江 安達 151400）

1 引言

黑龍江省地處吉林玉米種植帶區域，以種植春玉米為主[1]。截至2018年，黑龍江省玉米總產量達3982萬噸，種植面積達6317千公頃 （國家統計局，2018）。玉米是一種極具代表性的喜溫作物，它在整個生育期對低溫都極為敏感[2]。黑龍江省地處我國中高緯度地區，玉米生長季熱量資源嚴重不足，低溫冷害頻繁發生。自1980年以來，黑龍江省的氣候明顯變暖，近50 a內氣溫上升幅度高達1.0℃[3]，原來的玉米種植區低溫冷害的發生頻率及強度也隨之下降[4，5]。然而，氣候變化并非是一成不變的，溫度年際變化頻繁發生，上述因素不僅在一定程度上抵消了積溫增加的效應，而且大大增加了低溫冷害發生的風險概率[1]。因此，低溫冷害仍是威脅玉米正常生長發育和糧食安全生產的主要農業氣象災害之一。

天氣指數保險是保險合同以指數為基礎，這種指數表示單一或者多個氣象要素對農作物的損害程度，如果指數達到一定水平且對農產品造成一定影響，投保人就能得到相應標準的賠償[6]。作為一種新型的保險模式，天氣指數保險具有客觀性、科學性、及時性等特性[7]。由于我國氣象災害頻發、重發，同時農業保險市場態勢較好，各地的天氣指數保險快速發展[7]。天氣指數保險既有利于規避農業生產風險，又有利于提高防災減災和災后恢復能力，自上世紀90年代以來國內外研究人員逐漸重視這一領域，同時進行了許多研究工作。然而，目前并未見到有關黑龍江省玉米低溫冷害天氣指數保險的文章。因此，本文依據黑龍江省玉米生長季低溫冷害造成的玉米產量損失和氣象要素之間存在的關系，建立玉米低溫冷害保險的賠付標準，并根據玉米低溫冷害風險分布情況，設計不同保險費，從而完成黑龍江省玉米低溫冷害天氣指數保險產品的設計。

2 資料與方法

2.1 研究資料

本研究選取了黑龍江省22個玉米農業氣象站（圖 1）1981－2016 年 5－9 月逐日平均氣溫數據、各站點經緯度數據，來源于黑龍江省氣象局。本文所選用的玉米產量數據，為氣象部門建有玉米農業氣象觀測站的22個縣 （市）1981－2016年的玉米單產數據，來源于黑龍江省統計局。

圖1 黑龍江省玉米農業氣象站點分布

利用Excel和SPSS軟件對玉米低溫冷害天氣指數和單產數據進行處理分析。

2.2 玉米低溫冷害天氣指數選取

天氣指數通常由一個或幾個氣象條件構成[7]。需具備客觀性、可驗證性及穩定性。指數可以利用已經頒布的農業氣象災害標準[7]，因此本研究利用最新的玉米低溫冷害國家標準[8]，選擇玉米生長季5－9月逐月平均氣溫之和的距平（△T）為玉米低溫冷害天氣指數。玉米低溫冷害強度等級[8]見表1。

玉米低溫冷害天氣指數的計算公式為

表1 黑龍江省玉米低溫冷害強度等級

式中，△T為玉米低溫冷害天氣指數 （℃），Tt為當年5－9月逐月平均氣溫之和（℃），T為5－9月逐月平均氣溫之和的多年平均值（℃）。

2.3 玉米單產減產率的計算

利用歷年的單產數據和步長為5 a的滑動平均法[9]，模擬趨勢產量Yt，進而分離出氣象產量，從而依據式（2）計算相對氣象產量Yw：

式中，Y為農作物實際產量[7]。若Yw＜0，則意味著氣象條件對作物生長不利，導致作物減產，視－Yw≥5%的年份為減產年[10]，－Yw即為減產率。

2.4 玉米低溫冷害天氣指數保險費的計算

依據張韻啟[11]的研究，保險公司的賠付率為保險公司期望賠償的金額與保險費的比值，農戶的損失為不發生賠付時需投入的保險費。將這兩個指標之和設為目標函數D，當D最小時，保險公司的賠付額和農戶的損失相對最小，就得出了最優保險費，則目標函數的計算公式為

式中，B為保險金額，P為災害概率，M為保險費，L為玉米類保險投保人負擔比[11]。由于農業保險費不宜過高，假設約束條件為：0＜M＜100 (單位：元/666.7m2)。依據黑龍江省農業保險保費補貼管理辦法的相關規定，玉米類保險中央、省、市縣財政補貼分別為40%、25%和15%，投保人負擔20%，所以L為0.2。

2.5 玉米低溫冷害風險評估

通常可以利用氣象災害風險評估的結果來修訂天氣指數保險的保險費率。玉米低溫冷害風險大小既受天氣影響，又與種植區經濟發展水平、地理環境、抗災能力等因素有關[12]。因此，各種植區風險程度不同，應實施不同費率，本文依據不同種植區低溫冷害風險評估結果來修訂保險費率，以實現不同區域差別費率。致災因子危險性是本文玉米低溫冷害風險評估的主要考慮因素[13]。

致災因子危險性通常用來表示發生災害的可能性。該指數越大，表示發生玉米低溫冷害的可能性越大，從而對玉米產量的影響程度也就越大；反之亦然。本文用玉米低溫冷害頻率和強度乘積的絕對值來表示玉米低溫冷害致災因子危險性，

依據式（5）對致災因子危險性指數做無量綱化處理，進而得到黑龍江省玉米低溫冷害風險指數[12]。

式中，H′為無量綱化值，H為原始值，Hmin為最小值，Hmax為最大值。風險指數越大，表示玉米低溫冷害出現的風險越大，設定風險Ⅰ區為風險指數較高區域，該區費率浮動系數較大；其余為風險Ⅱ區。

3 結果與分析

3.1 玉米低溫冷害天氣指數與減產率的關系

將22個站點玉米低溫冷害明顯發生且導致玉米發生減產的樣本作為典型的低溫冷害樣本，利用SPSS軟件，分析黑龍江省玉米低溫冷害天氣指數與玉米單產減產率的相關性，并利用線性回歸分析法得到兩者的回歸方程

式中，X為減產率（%）。當玉米低溫冷害天氣指數減少1℃時，減產率增加5.85%。根據式（6）計算得到表2，即不同玉米低溫冷害天氣指數對應的玉米單產減產率。

表2 黑龍江省玉米低溫冷害天氣指數與減產率的關系

3.2 玉米低溫冷害風險評估

利用ArcGIS反距離權重法插值，采用自然斷點法將黑龍江省玉米低溫冷害風險等級分為3級，即高風險區、中等風險區和低風險區，詳見圖2。黑龍江省玉米低溫冷害高風險區位于西部的龍江，以及西南部的哈爾濱和雙城；中等風險區位于中部的勃利和集賢、東部的饒河，以及北部和西部大部分地區。設定以上地區為風險Ⅰ區，其余為風險Ⅱ區。根據經驗，設定風險Ⅰ區和風險Ⅱ區的費率浮動系數分別為 1.2和 1.0。

圖2 黑龍江省玉米低溫冷害風險空間分布

3.3 玉米低溫冷害天氣指數保險的賠付標準

按照目前玉米市場平均價格1.7元/kg推算，在正常情況下黑龍江省玉米畝產值可以達到850元，保障水平取70%，即保險金額B為595元/666.7m2。當玉米低溫冷害天氣指數 （△T）≤－5℃時起賠。將△T代入式（6）中即得減產率，賠償金額則為保險金額與減產率的乘積。將B=595元/666.7m2，災害概率P=0.08，L=0.2代入式（3）中，利用 LINGO 軟件得到保險費M=16元/666.7m2。不同地區保險費為保險費率浮動系數與M的乘積，則風險Ⅰ區和風險Ⅱ區的保險費分別為19.2元/666.7m2和16元/666.7m2。

4 結論和討論

本文依據中國氣象局 《農業保險氣象服務指南—天氣指數設計》的規定，利用最新的玉米低溫冷害國家標準，選擇5-9月逐月平均氣溫之和的距平作為黑龍江省玉米低溫冷害天氣指數，通過分析黑龍江省22個玉米農業氣象站的氣溫數據和玉米產量數據，利用SPSS軟件得到了黑龍江省玉米低溫冷害天氣指數與減產率的關系模型，從而構建了黑龍江省玉米低溫冷害天氣指數保險的賠付標準，并利用黑龍江省玉米低溫冷害風險評估結果修訂了保險費。

通過統計分析歷史數據，本文建立了黑龍江省玉米低溫冷害天氣指數模型，在理論上設計了黑龍江省玉米低溫冷害天氣指數保險產品，今后還有待開展試點應用來驗證和修訂有關參數。