甘肃省春玉米灾损风险评估

齐 月, 阳伏林, 姚玉璧, 王润元, 房世波

(1.中国气象局 兰州干旱气象研究所/甘肃省干旱气候变化与减灾重点实验室/中国气象局干旱气候变化与减灾重点实验室,兰州730020; 2.中国气象科学研究院 生态环境与农业气象研究所, 北京 100081)

全球变暖已成为不争的事实。伴随着气候变暖和降水变异的加剧，极端天气出现的频次在增加，强度也在加大[1-3]。干旱是最严重最为常见的自然灾害之一，其发生频率高、持续时间长、影响范围大，给国民经济特别是农业生产造成巨大的损失。据统计，我国农业受旱面积20世纪50年代年均为1 330万hm2，90年代年均0.243亿hm2，因旱损失粮食由年均43.5亿kg增加到195.7亿kg，给农业生产和人民生活带来严重影响[4-5]。甘肃省地处西北内陆地区，属温带大陆性气候。受地理因素影响，海洋暖湿气流难以在该地形成降水[6]，属于典型干旱半干旱区。由于干旱少雨导致农业产量低而不稳，干旱已经成为制约当地农业生产的主要农业气象灾害之一，并且近几十年全球变化使西北地区作物受干旱影响愈加严重[6]。

玉米是中国四大粮食作物之一，仅次于小麦。中国玉米种植面积和总产量仅次于美国，居世界第二[7-9]。中国现有的0.93亿hm2耕地中，玉米面积约为0.21亿hm2，占总耕地面积的1/5[10]。2017年甘肃省玉米种植面积约100万hm2，占全省粮食种植面积的35%，是甘肃省的主要粮食作物。在全球气候变化和干旱化趋势的背景下，干旱对甘肃省玉米生产的影响越来越严重，造成玉米大面积减产，给农业生产带来巨大损失。

目前，大部分气象灾害风险研究主要考虑4个要素，即致灾因子风险分析、承灾体易损性、灾情损失评估以及防灾减灾能力等进行综合评价，由于方法不同或权重不同，并且没有考虑到农业系统的核心——作物的损失状况，其评价结果差异较大[11-21]。在以往的研究中，多集中于某一县市的玉米风险进行研究，或者对玉米的生态适应性、资源利用进行分析，未见对甘肃省各县春玉米干旱灾害风险研究。本文根据风险理论，采用风险评估技术和方法，研究甘肃省71个市(县)不同年型、不同减产率下春玉米产量在干旱气候条件下的波动和减产的风险概率，并对甘肃省春玉米干旱灾害风险进行区划。

1 数据与方法

1.1 数据来源

本文的所用资料来自于甘肃省农村统计年鉴，包括1980—2011年甘肃省71个县(市)的春玉米种植面积和总产量资料；少数地区的资料可能有个别时段1～2 a资料缺测，统计时根据实际产量年数计算。根据蒲欣东等[22]的划分方法，将甘肃省划分为河西地区、陇中地区、陇南地区、陇东地区和甘南地区5个区域。

1.2 研究方法

1.2.1 作物气象产量分离 作物气象产量分离是目前研究气象因子与作物产量之间的关系最常用的研究方法。由于长时间序列的作物产量变化不仅与气象因子有关，也与科技进步、物质投入、环境、政策等有密切关系，其中科技进步水平对粮食单产的影响力最大。因此，一般将作物产量分解为趋势产量、气象产量和随机误差3部分[23-24]。趋势产量可看作是反映某一历史时期某一生产区域生产力发展水平的长周期产量分量，气象产量是以气象要素为主的短周期变化因子影响的产量分量。

Y=Yt+Yw+ε

(1)

式中：Y为实际产量(t)；Yt为趋势产量(t)；Yw为气象产量(t)；ε为其他随机产量(t)，又被称为“随机噪声”，可忽略不计。因此，一般研究中将实际作物产量直接分解为趋势产量和气象产量。即：

Y=Yt+Yw

(2)

本研究采用直线滑动平均法对趋势产量进行模拟[24]。将某一阶段内产量的时间序列的变化看作线性函数，成一直线。随着阶段的前进和推移，直线不断改变位置、后延滑动，从而反映产量的历史演变趋势。依次求取各阶段的线性回归函数，则经过某一时间点上各线性函数的均值，即为该时间点的趋势产量值。设某阶段的线性趋势方程为：

Yi(t)=ai+bit

(3)

式中：方程个数i=n+k+1；n为样本序列个数；k为时间步长，为了消除短周期波动的影响，滑动步长k取11；t为时间序号。本文选取1980—2011年的玉米单产资料，因此，共有20个趋势方程。计算每个方程在t点上的函数值Yi(t)，这样每个t点上分别有q个函数值，q的多少与n有关。当k≤n/2，则q=1,2,3,…,k,…,k,…,3,2,1；q连续为k的个数等于n-2(k-1)。当k≤n/2，则q=1,2,3,…,n-k+1,…,n-k+1,…,3,2,1；q连续为n-k+1的个数等于2k-n。本文单产序列共32 a，因此1980—1989年每个时间点t上的函数值个数q依次取1，2，…，10；1990—2001年这12 a每个时间点t上各有11个函数值，即q=11；2002—2010年每点上t依次有10，9，…，1个函数值。然后，再求出每个t点上q个函数的均值：

(4)

Yw=Y-Yt

(5)

由气象产量(Yw)和趋势产量(Yt)的比值可得到相对气象产量，见公式(6)。相对气象产量表明粮食波动的幅度，不受时间和地域的限制，具有可比性，因此能较好地反映气象因子对产量的影响。

(6)

式中：Yr为相对气象产量(%)。

当实际产量低于趋势产量时，相对气象产量为负值，称为减产率；反之，当实际产量高于趋势产量时，相对气象产量为正值，称为增产率。

1.2.2 相对气象产量正态分布型判别 相对气象产量与气象产量相比能更好地描述以气象要素为主的各种短期变动因子对产量序列的影响[25-26]。因为相对气象产量消除了生产力水平的影响，主要反映的是气象因子对产量的影响规律，而由于气候因子一般具有正态分布的特点，相对气象产量也应该具有正态分布的特征。为验证这一观点，本文采用偏度—峰度检验法对所需研究资料进行正态检验[24,27]。峰度和偏度是研究未知曲线性质的基本手段。理论上，正态分布的偏度：

(7)

正态分布的峰度：

(8)

表1 偏度-峰度检验临界值

经检验，研究区域有71个县市，其中9个县市的减产率未通过α=0.01水平的正态性检验，通过检验的占总体的88%，表明甘肃省玉米减产率基本服从正态分布。对于不符合正态分布的个别站点数据进行偏态分布正态化处理[26]。

1.2.3 灾害减产的风险评价指标 相对气象产量的负值主要是由于不利气象条件(灾害)所造成的减产[27]。在诸多灾害中，干旱是影响研究区最大的气象灾害，是造成春玉米减产最严重的自然灾害。大范围的春玉米相对气象产量波动的负值(减产率)主要是以干旱为主的一些因素造成的。因此，分析本研究区春玉米历年相对气象产量的波动，减产率和减产风险概率在很大程度上可以反映干旱对春玉米产量影响的风险程度。

(1) 减产率指标。相对气象产量是一个相对值，表明实际产量偏离趋势产量的波动幅值。实际单产低于当时趋势产量的百分率称为“减产率”，高于趋势产量的百分率称为“增产率”。对于某一相对气象产量{xi}，定义xi=0.0%为丰歉临界值，xi<0.0%的年份为歉收年份；xi=-5.0%为成灾临界值，xi<-5.0%的年份为成灾年份[22]。为了表示一个地区多年平均的干旱程度对产量的影响，采用历年平均减产率指标P。

(9)

式中：xi，n分别为不同年型对应的相对产量序列和年数。当n为歉收年份的年数时，得到歉年平均减产率；当n为成灾年份时，得到灾年平均减产率。

(2) 不同年型减产变异系数指标。变异系数表示减产的幅度偏离其平均值的程度，是反映研究区春玉米产量是否稳定的指标[25]。对于同样大小的平均减产率，年际间减产率差别大和小，其气象灾害的风险不同，波动性大，气象灾害的风险越大。减产年减产率的波动性，称为减产率变异系数。

(10)

式中：xi为减产率序列；P为平均减产率；n为对应xi的所有减产年。

变异系数是均方差和数学期望的比值，描述相对气象产量序列数值分散的程度如何，是高度集中在某个范围内，还是比较均匀地分散在整个变化范围内。某年型减产变异系数表明产量的分布状况，系数大，说明大丰收和大减产的年份多，产量稳定性差。

(3) 风险概率指标。风险概率指标是指根据分布函数曲线计算的春玉米单产不同增产率和减产率出现的概率[25]。正态分布函数为：

(11)

则随机变量x在区间(x1,x2)内发生的概率为：

(12)

本文按歉年和灾年两种年型对各风险评价指标进行分析。其中，歉年减产率为相对气象产量小于零，即歉收年份的减产率平均值；灾年减产率为相对气象产量小于5%(减产率大于5%)的成灾年份的减产率平均值[28]。

将甘肃省71个县(市)的春玉米平均减产率按照减产率高低分为5个级别：低值区(0～20%)、次低值区(20～30%)、中值区(30～40%)、次高值区(40～50%)、高值区(≥50%)。

1.2.4 综合风险区划方法 本文以甘肃省71个县市为样本，选取每个县(市)歉年平均减产率、灾年减产率变异系数、年减产率≥10%的风险概率3个指标作为风险指标，通过K-平均聚类算法对甘肃省春玉米进行减产风险区划，将其分为低、中、较高、高4类。

2 结果与分析

2.1 春玉米歉年平均减产率和变异系数分布规律

2.1.1 春玉米歉年平均减产率分布规律 从图1中可知，春玉米歉年平均减产率在该区的分布有明显的区域性：该区域以黄河为界，河西地区高于其他地区，大部分高值区和次高值区位于河西地区。该区域的北面高于南面；其中河西地区的12个县、陇中地区的24个县、陇东地区的9个县、陇南地区的7个县、甘南地区的1个县属于低值区，平均减产率在5.8%～20%的范围内；而除河西地区的古浪县、金昌县和民勤县，陇东地区的华池县和环县外，陇中地区其他9个县属于次低值区，平均减产率在20%～30%的范围内；中值区主要有3个县，全部位于陇中地区，平均减产率在30%～40%的范围内；次高值区有5个县，其中河西地区3个，平均减产率在40%～50%的范围内；平均减产率>50%的高值区有9个县，除了陇中地区的定西市、渭源市和漳县，陇东地区的静宁县，其他大部分位于河西地区。

图1春玉米歉年平均减产率分布

2.1.2 春玉米歉年减产率的变异系数分布规律 变异系数的大小反映了该地区春玉米产量的稳定状况，变异系数高表明该地区的玉米产量年际波动较大，玉米生产容易受外界条件的影响[24]。按变异系数的大小将其分成3个等级，分别为低值区(≤0.89)，中值区(≤1.26)，高值区(≤2.40)。

从图2可以看出，变异系数的分布没有明显的区域性和连续性，各个范围相间分布。高值区主要分布在甘肃省的中部地区，包括河西地区的酒泉市、张掖市、武威市、民勤县、金昌县、永昌县；陇中地区的皋兰县、临夏市、白银市、康乐县、和政县、临洮县、渭源县；陇东地区的正宁县、静宁县；陇南地区的徽县，大部分位于河西地区和陇中地区。低值区的变异系数值在0.89以下，主要分布在河西地区的阿克塞哈萨克族自治区、敦煌市、嘉峪关市、民乐县、肃北县；陇中地区的定西市、广河县、会宁县、靖远县、张家川县、漳县、陇西县；陇东地区的环县、泾川县、庄浪县；陇南地区的西和县，大部分位于河西地区和陇中地区。其他地区的变异系数为0.89～1.26，属于中值区。

图2春玉米歉年减产率的变异系数分布

2.2 春玉米灾年平均减产率和变异系数的分布规律

2.2.1 春玉米灾年平均减产率分布规律 从图3中可见，春玉米灾年平均减产率分布同样具有明显的区域性，其分布特征与歉年平均减产率大体一致。以黄河为界，河西地区高于其他地区，大部分高值区和次高值区位于河西地区。该区域的南面高于北面；其中河西地区的5个县、陇中地区5个县、陇南地区的3个县、陇东地区的7个县，属于低值区，平均减产率在5.8%～20%的范围内；而河西地区的4个县、陇中地区11个县、陇东地区的3个县、陇南地区的4个县属于次低值区，平均减产率在20%～30%的范围内；中值区主要有10个县，陇东地区环县、天水市、武山县、通渭县、永登县5个县，河西地区张掖市、民勤县2个，甘南地区的舟曲县以及陇中地区的临夏市和临夏县，平均减产率在30%～40%的范围内；次高值区有8个县，其中河西地区5个，平均减产率在40%～50%的范围内；平均减产率>50%的高值区有11个县，大部分位于河西地区和陇中地区。

2.2.2 春玉米灾年平均减产率的变异系数分布规律 根据春玉米灾年平均减产率变异系数的大小将其分成3个等级，分别为低值区(≤0.71)，中值区(≤0.97)，高值区(≤1.67)。

从图4可以看出，灾年平均减产率的变异系数的分布没有明显的规律性和区域性。高值区主要有白银市、皋兰县、古浪县、和政县、徽县、金昌县、静宁县、酒泉市、康乐县、临洮县、临夏市、民勤县、清水县、通渭县、渭源县、武威市、永昌县、榆中县和张掖市，大部分位于河西地区东部和陇中地区。中值区的变异系数在0.71～0.97，主要有河西地区的安西县、山丹县、肃南县、敦煌市，陇南地区的成县、崇信县、甘南地区的迭部县、舟曲县，陇中地区的甘谷县、积石山保安族东乡县、兰州市、岷县、秦安县，陇东地区的合水县、灵台县、宁县，陇南地区的礼县、文县，主要分布在变异系数高值区的周边。其他地区的变异系数分布在0.71以下，属于低值区。歉年减产率变异系数高的地区大多是歉年平均减产率大的地区。说明甘肃中部地区虽然减产率较小，但其产量波动的幅度要大于减产率大的地区。

图3春玉米灾年平均减产率分布

图4春玉米灾年减产率的变异系数分布

2.3 春玉米减产率≥5%的风险概率分布

减产率≥5%的风险概率是指偏离趋势产量负值的相对气象产量达到或者超过5%以上出现的概率[24]。从图5中可以看出，减产率≥5%的风险概率在本区分布具有明显的连续性和区域性。除河西地区的玉门市、金塔县、高台县、临泽县和武威市，及陇东和陇南地区的东部边缘区外，其他大部分地区属于次高值区和高值区，减产率≥5%的高风险区主要位于河西的大部分地区、陇中地区、甘南地区及陇东和陇南的中西部地区，整体上看，甘肃省的中西部地区减产率≥5%的风险高于东部边缘区域。

图5春玉米减产率≥5%的风险概率

2.4 春玉米减产率≥10%的风险概率分布

减产率≥10%的风险概率是指偏离趋势产量负值的相对气象产量达到或者超过10%以上出现的概率[24]。由图6可见，次高值区和高值区主要有河西地区的阿克塞哈萨克族自治县、敦煌市、安西县、古浪县、金昌县、民乐县、肃北县、嘉峪关市、肃南县，陇中地区的渭源县、景泰县、临夏市(县)、渭源县、漳县、张家川县，陇东地区的环县、静宁县、西峰市，陇南地区的西和县，大部分位于河西和陇中地区；低值区和次低值区主要有河西地区的金塔县、临泽县、武威市、玉门市，陇南地区的成县、康县、两当县、武都县、陇东地区的广河县，陇东地区的崇信县、合水县、华亭县、泾川县、灵台县、宁县、正宁县、庄浪县，陇中地区的清水县、永靖县、靖远县，大部分位于陇东地区；其他地区属于中值区。整体来看，甘肃中西部地区减产率≥10%的风险概率高于东部地区。

图6春玉米减产率≥10%的风险概率

2.5 春玉米减产率≥15%的风险概率分布

减产率≥15%的风险概率是指偏离趋势产量负值的相对气象产量达到或者超过15%以上出现的概率[24]。从图7中可以看出，减产率≥15%的风险概率整体偏低，大部分地区属于低值区和中值区。高值区和次高值区分布在河西地区的阿克塞哈萨克族自治县、安西县、古浪县、嘉峪关市、金昌县、肃北县、肃南县、民乐县，陇中地区的景泰县、漳县，陇东地区的静宁县、渭源县、西峰市，陇南地区的西和县；其他地区主要以中低值区为主。整体来看，河西地区减产率≥15%高于其他地区。

图7春玉米减产率≥15%的风险概率

2.6 春玉米产量灾害风险区划

从图8可以看出，甘肃省春玉米单产减产风险区呈整体分散，小范围聚集连片的特点。较高风险区域较多，分布在37个县(市)，其次为高风险区域，中风险区域分布最少。甘肃省东部和河西地区中部为春玉米减产高风险区，陇中、陇东大部分地区为春玉米减产较高风险区，而中、低风险区主要分布在河西地区的西部和南部、陇中地区的东部。

甘肃省春玉米种植模式主要是春种秋收，影响春玉米生长发育的主要气象灾害包括苗期的霜冻害、生长中后期的干旱(主要为夏、秋旱)。春玉米生长的中后期需水量较大，占全生育期的50%以上，此时遇到干旱的发生，使春玉米光合作用降低，产量明显下降。甘肃省中西部河西地区属于灌溉农业，虽然各类减产率均属于高风险区，但是可通过灌溉缓解干旱对春玉米生长发育的影响。甘肃省其他地区属于雨养农业，干旱是导致该地区减产的主要因素。

图8春玉米产量灾害风险区划

3 讨 论

以往的研究大多考虑风险研究的4个要素即致灾因子风险分析、承灾体易损性、灾情损失评估以及防灾减灾能力等进行综合评价，由于权重的不同，并且没有考虑到农业系统核心——作物灾损状况，其评价结果差异较大。基于以风险理论为基础，采用风险评估技术方法，运用甘肃省各县的春玉米实际的产量资料对春玉米干旱风险进行区划，避免权重的影响，使得评价结果更符合实际。

本文采用实际产量偏离趋势产量负值的百分比主要是由于受气象条件影响造成的春玉米产量，干旱是甘肃省春玉米受到的最普遍最严重的自然灾害，该地区春玉米历年灾年减产率基本上是干旱造成的。河西地区为典型的干旱区，中部地区为春玉米减产的高风险区，变异系数也较大，但是河西地区的西部和南部属于低风险区域，由于河西地区属于灌溉农业，干旱发生时降水减少，通过灌溉来缓解干旱导致的减产，使其干旱风险降低，这与白月明[6]、姚小英[29]等的结果相一致。陇中地区是典型的半干旱区，以雨养农业为主，干旱直接导致作物减产，与上述的结果相一致。本文只是对甘肃省春玉米干旱风险做了一个尝试下的探讨，还需要进行深入的分析，根据不同的趋势产量分析法进一步进行风险区划分析。

4 结 论

(1) 春玉米歉年和灾年的平均减产率分布均具有明显的区域性，大致以黄河为界，河西地区高于其他地区。大部分的高值区和次高值区位于河西地区；中值区位于陇中地区；低值区和次低值区大部分位于陇中、陇东和河西地区。整体是北面高南面低。

(2) 春玉米单产的歉年和灾年减产率的变异系数分布没有明显的区域性和连续性，高值区和低值区主要分布在河西和陇中地区，说明河西地区和陇中地区的年平均减产率波动较大，春玉米容易受到干旱的影响。

(3) 减产率≥5%的风险概率均具有明显的连续性和区域性，高风险区主要分布在甘肃省的中西部地区；减产率≥10%的高风险区主要分布在河西和陇中地区，主要位于甘肃省的中西部地区；减产率≥15%的风险概率整体偏低，大部分地区属于低值区和中值区，高值区、次高值区主要分布在河西地区。

(4) 甘肃省春玉米减产风险区呈整体分散，小范围聚集连片的特征。河西地区的玉门市、金塔县、临泽县、张掖市、武威市，陇中地区的靖远县、永靖县以及陇东地区的东部各县春玉米减产属于高风险区，陇中、陇东大部分地区为春玉米减产较高风险区，其他县市为中、低风险区。