1. 引言
自然灾害风险是由自然灾害系统自身演变而导致未来损失的不确定性,对风险进行评估,就是借助于不确定性分析方法,定量的表达这种不确定性的大小,从而实现风险等级划分,并以此作为应对风险的重要依据 [1] 。目前全球各种自然灾害造成的损失中,暴雨洪涝灾害(含热带气旋)占60%,干旱占15%,地震占15%,其它占10% [2] 。在气候变暖的背景下,全球多数地区极端强降水事件有增多的趋势 [3] - [14] ,给水稻生产带来巨大的影响 [15] [16] [17] 。
江苏省宿迁市地处淮河气候过渡带,天气气候变化异常,近30年气候资料分析结果表明,水稻种植季节(5~10月份)总降水量呈增加趋势,同时暴雨以上强度的降水日数明显增多,加之受淮、沂沭泗流域降水和泄洪的影响,暴雨洪涝灾害对水稻生产的危害加重。在这样的大背景下,开展水稻暴雨洪涝灾害风险等级划分研究,对布局粮食安全生产,指导防灾抗灾工程设施建设,科学合理的实施防灾抗灾对策措施,增强防灾抗灾能力,减少灾害损失,提高粮食种植的社会和经济效益有着十分重要的战略意义和实用价值。
本研究利用宿迁水稻暴雨洪涝灾损预报模型输出的灾损率预测结果和重灾发生乡镇比例两个因素,设定灾害等级划分阈值,实现对宿迁水稻种植区进行暴雨洪涝灾害风险等级划分的目的,并统计各风险等级的历史灾害参数特征,用以指导本地区水稻种植。
2. 资料来源及研究方法
2.1. 资料及来源
气象数据:2003~2013年宿迁境内各乡镇降水、大风及常年平均值,由宿迁市气象局提供,并经过严格的质量控制。
水稻产量数据:2003~2013年宿迁各乡镇水稻产量数据,来源于宿迁市及各县区的统计年鉴;水稻生育期数据由宿迁市农委提供。
水稻灾损数据:2008~2013年水稻灾害损失 ≥ 10%的重灾案例数据,来源于承担宿迁市农业保险业务的宿迁市平安保险公司。规定:当某乡镇出现暴雨洪涝灾害,有一个或以上自然村出现水稻灾损重灾案例时,视为该乡镇发生一次水稻暴雨洪涝灾害重灾案例。
2.2. 研究方法
2.2.1. 水稻灾损率计算
沈伟等 [18] 在宿迁地区水稻暴雨洪涝灾害灾损分析研究中,基于孕灾环境的敏感性、承载体易损性和暴露度等因素将宿迁水稻种植区划分为10个区域(见图1),并利用公式(1)分区计算各乡镇的水稻暴雨洪涝灾害实际灾损率,其中Y为实际产量,
为采用最优曲线拟合法模拟各乡镇水稻产量趋势。将
Figure 1. Rice cropping regionalization in Suqian City
图1. 宿迁水稻区划图
计算结果结合重灾案例发生时间分析了得到:水稻孕穗扬花期(7月21日至9月10日)和灌浆成熟期(9月10日至10月20日)是宿迁地区水稻种植受暴雨洪涝灾害威胁最大的阶段,而在分蘖拔节期及以前,暴雨洪涝灾害对宿迁水稻种植基本不会造成大的影响。
(1)
2.2.2. 水稻暴雨洪涝灾害指数
董家山 [18] 在沈伟 [19] 的基础上,为了定量的描述水稻孕穗扬花期和灌浆成熟期暴雨洪涝对水稻的影响,选取了雨日、暴雨日数、总雨量、暴雨量这4个因子,并在灌浆成熟期考虑暴雨日的大风影响,定义了水稻暴雨洪涝指数
和
。
(2)
(3)
(2)式和(3)式分别为孕穗扬花期和和灌浆成熟期的暴雨洪涝指数公式。
、
、
、
分别为孕穗扬花期(或灌浆成熟期)的降水日数、暴雨日数、总降水量、总暴雨量,
、
、
、
为相应统计量的常年平均值,
为灌浆成熟期暴雨日的瞬时极大风速。
2.2.3. 面板回归模型
面板数据 [20] 具有截面和时间序列两个维度,可以更好的控制样本差异,减少变量间的多重共线性。面板回归模型的一般形式如下:
(4)
利用面板回归模型分析水稻灾损与暴雨洪涝指数的相互关系。(4)式
在文中代表某区域内乡镇i第t年水稻因暴雨洪涝造成的损失率,
表示乡镇i第t年第k个关键生育期内的暴雨洪涝灾害指数(k = 1或2,代表孕穗扬花期或灌浆成熟期),
为的相关系数,
为随机变量。
Hausman检验 [21] [22] 用于检验暴雨洪涝指数
和误差项
之间是否存在相关关系,当
与
相关时,乡镇之间的差异用不同的截距来反映,即固定效应模型;反之,选用随机效应模型。文中面板回归分析的工具为EViews8。
3. 水稻暴雨洪涝灾损模型建立与效果检验
3.1. 水稻暴雨洪涝灾损模型建立
以水稻实际灾损率为研究对象,水稻孕穗扬花期和灌浆成熟期的暴雨洪涝灾害指数作为影响因子,代入模型的一般方程(4)式。通过Hausman检验,确定面板回归模型选用随机效应模型(图2,IV区的检验结果,其它各区图略),建立宿迁水稻暴雨洪涝灾害损失率预测模型(图3,IV区建模结果,其它各区图略),即各区水稻暴雨洪涝灾损灾损模型方程(表1)。将2008~2013年各乡镇的孕穗扬花期和灌浆成熟期暴雨洪涝指数带入表1中对应的模型方程,即可输出水稻暴雨洪涝灾害损失率的模型预测结果。
Figure 2. Hausman test in region IV
图2. IV区Hausman检验结果
Figure 3. The modeling results in region IV
图3. IV区建模结果
Table 1. Disaster loss model caused by rainstorm and flood disasters in each rice planting region
表1. 各区水稻暴雨洪涝灾害灾损模型方程
3.2. 模型效果检验
从表2可见,模型结果的平均绝对误差为0.07,平均相对误差为5%,说明总体上模拟效果较好。从分区情况来看,IX区、VII区、II区和VIII的效果明显好于其他各区,绝对误差均 ≤ 0.03,相对误差 ≤ 2%;IV区效果较差,相对误差达到28%。
Table 2. Errors of the model results from 2008 to 2013
表2. 2008~2013年模型结果误差统计
模型对于2008~2013年各区水稻暴雨洪涝灾害灾损率 ≥ 10%的重灾乡镇拟合结果检验结果见表3,历史拟合准确率除IX区为89.7%以外,其他9个区均>90%,平均拟合准确率为94.1% (漏报率为2.7%,空报率为3.2%)。
从表3可见2014年试用效果较好,重灾乡镇平均准确率为95.2%,没有出现漏报,空报率为4.8%,IV区试用效果最差,准确率只有71.4%,其次是X区,准确率也仅有75.0%。
4. 水稻暴雨洪涝灾害风险等级划分及应用
4.1. 阈值确定与风险等级划分
本研究进行暴雨洪涝灾害风险等级划分主要参考受灾程度和重灾面积两个参数,受灾程度用模型结果输出的暴雨洪涝灾害灾损率来表示,重灾面积用重灾乡镇比例来表示,即YT = 各区模型计算出的平均灾损率*各区重灾乡镇比例。重灾乡镇比例为2008~2013年出现灾损率 ≥ 10%的重灾乡镇个数/各区包含的乡镇数*年数。
从表4可见,YT的平均值为−0.81,最大值为−0.21,最小值为−2.23,负值越大表示受灾越严重。以YT的平均值、(YT平均值 + YT最小值)/2、(YT平均值 + YT最大值)/2为三个阈值,即Y1 = −(2.23 + 0.81)/2 = −1.52、Y2 = −0.81、Y3 = −(0.81 + 0.21)/2 = −0.51,将宿迁地区水稻种植区划分为四级风险等级。一级风险:YT ≤ Y1;二级风险:Y1 < YT ≤ Y2;三级风险:Y2 < YT ≤ Y3;四级风险:YT > Y3。2008~2013年水稻暴雨洪涝灾害历史风险等级划分结果和区域分布情况见表4和图4。
4.2. 各风险等级历史灾害分布特征统计
从图4可见,2008~2013年宿迁地区水稻暴雨洪涝灾害一级风险区位于沿淮河北岸和洪泽湖西岸;二级风险区位于泗洪县西部高岗地区和新沂河以北地区;三级风险区位于黄墩湖滞洪区及周边地区,四级风险区位于骆马湖以东、新沂河以南和京杭大运河沿岸及以北地区。
表5给出了根据2008~2013年水稻暴雨洪涝灾害情况统计得到的各个风险等级的灾害特征参数。从表中可见,一级风险区暴雨洪涝灾害造成的平均灾损率、重灾乡镇比例、受灾最重乡镇的平均灾损率和灾害最重年份的平均灾损率分别为−3.03%、0.61%、−5.94%和−7.46%,灾害最重年份灾情最重乡镇的平均灾损率高达到−23.27%;二级和三级风险区的各项灾害参数依次下降;四级风险区的平均灾损率、受灾最重乡镇的平均灾损率和灾害最重年份的平均灾损率均最低,分别为−1.30%、−2.58%和−3.21%,但重灾乡镇比例比三级略高,灾害最重年份灾情最重乡镇的平均灾损率比三级风险区偏高,我们认为四级风险区的乡镇数较多,相关参数更具有代表性,三级风险区的乡镇数较少,相关参数具有一定的偶然性,随着以后在应用中历史资料的进一步补充,这一缺陷将得到弥补。
Table 3. Fitting rate of towns whose rice loss rate is over 10% from 2008 to 2013 and trial accuracy of 14 years
表3. 2008~2013年灾损率 ≥ 10%的重灾乡镇拟合率及14年试用准确率统计
Table 4. The results of historical risk grades from 2008 to 2013
表4. 2008~2013年历史风险等级划分结果
Figure 4. Risk grades layout of rainstorm and flood disasters in Suqian city from 2008 to 2013
图4. 2008~2013年宿迁水稻暴雨洪涝灾害风险区分布图
Table 5. The historical disaster characteristic parameter of each risk grade
表5. 各风险等级历史灾害特征参数统计
Table 6. Application results in 2015
表6. 2015年应用结果统计
4.3. 水稻暴雨洪涝灾害风险等级划分结果的应用
用2015年各乡镇的I1和I2实时数据分别导入表1中对应区的模型方程,模型输出各乡镇的灾损率预测结果,将模型结果输出的各乡镇灾损率*各区重灾乡镇比例即得到各乡镇的YT值,统计得到2015年各区YT平均值,并根据各区YT阈值进行风险等级划分,划分结果见表6。
从表6可见,2015年宿迁大部分地区风险等级较低,多为四级风险,只有III区为二级风险,对照表5可就以得出各区2015年水稻暴雨洪涝灾害风险参数特征。
实况是当年宿迁地区降水比常年偏少,没有出现大范围的暴雨洪涝灾害,只有水稻抽穗扬花期(8月19~20日)洪泽湖西部出现局地性大暴雨,导致III区水稻暴雨洪涝灾害风险等级上升,造成III区平均灾损率为−2.4%,洪泽湖岸边的龙集镇平均灾损率达到−4.2%,孙园乡平均灾损率−3.1%。
5. 结论与讨论
利用降水(暴雨)强度、降水(暴雨)日数和累计降水(暴雨)量参数自定义的暴雨洪涝灾害指数,能较好的描述暴雨洪涝灾害特征,与水稻暴雨洪涝灾害损失率的相关性较好。利用时间和界面序列二维面板数据建立随机效应回归预测模型,在一定程度上克服了资料年代短的缺陷,模型效果较好。
直接应用模型输出的预测结果进行水稻暴雨洪涝灾害风险等级划分,减少了灾情调查环节,节省了大量的时间和人力物力,经检验和使用效果较好,能够反映出宿迁地区水稻暴雨洪涝灾害的相关特征,具有较高的实用价值。从历史风险区的划分结果来看:一级风险区位于沿淮河北岸和洪泽湖西岸;二级风险区位于泗洪县西部高岗地区和新沂河以北地区;三级风险区位于黄墩湖滞洪区及周边地区,四级风险区位于骆马湖以东、新沂河以南和京杭大运河沿岸及以北地区,与宿迁地区水稻暴雨洪涝灾害的实际情况相符合。
2008~2013年四级风险区的重灾乡镇比例和灾害最重年份灾情最重乡镇的平均灾损率比三级风险区偏高,这可能是因为三级风险区的乡镇数较少,相关参数具有一定的偶然性造成的,随着历史资料的进一步补充,这一缺陷将得到弥补。
基金项目
江苏省气象科研开放基金项目(项目编号:K201010)。