本文基于面板回归模型分析了宿迁水稻灾损率与暴雨洪涝灾害指数的相关性,建立了宿迁地区水稻暴雨洪涝灾损预测模型。以模型结果输出的水稻暴雨洪涝灾害损失率(灾损强度)和重灾乡镇比例(重灾面积)为依据,实现了宿迁地区水稻暴雨洪涝灾害的四级风险等级划分。通过统计分析历史灾情数据,得出各风险等级的灾害特征参数,作为应对水稻暴雨洪涝灾害风险的主要依据,用于指导生产实践。 By analyzing the relationship between the actual disaster loss rate and the index of storm flood disaster using panel regression model, prediction model of storm flood disaster of rice in Suqian area was designed. The risk rank of rice flood disaster was divided into 4 grades by the results of model. The disaster characteristic parameter was got by the statistical analysis of historical disaster data which can be used to reduce the flood disaster of rice.
沈伟,吴新胜,周航,王文清,张莹
宿迁市气象局,江苏 宿迁
收稿日期:2017年10月11日;录用日期:2017年10月23日;发布日期:2017年10月30日
本文基于面板回归模型分析了宿迁水稻灾损率与暴雨洪涝灾害指数的相关性,建立了宿迁地区水稻暴雨洪涝灾损预测模型。以模型结果输出的水稻暴雨洪涝灾害损失率(灾损强度)和重灾乡镇比例(重灾面积)为依据,实现了宿迁地区水稻暴雨洪涝灾害的四级风险等级划分。通过统计分析历史灾情数据,得出各风险等级的灾害特征参数,作为应对水稻暴雨洪涝灾害风险的主要依据,用于指导生产实践。
关键词 :水稻,暴雨洪涝,灾损,风险等级,风险特征
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
自然灾害风险是由自然灾害系统自身演变而导致未来损失的不确定性,对风险进行评估,就是借助于不确定性分析方法,定量的表达这种不确定性的大小,从而实现风险等级划分,并以此作为应对风险的重要依据 [
江苏省宿迁市地处淮河气候过渡带,天气气候变化异常,近30年气候资料分析结果表明,水稻种植季节(5~10月份)总降水量呈增加趋势,同时暴雨以上强度的降水日数明显增多,加之受淮、沂沭泗流域降水和泄洪的影响,暴雨洪涝灾害对水稻生产的危害加重。在这样的大背景下,开展水稻暴雨洪涝灾害风险等级划分研究,对布局粮食安全生产,指导防灾抗灾工程设施建设,科学合理的实施防灾抗灾对策措施,增强防灾抗灾能力,减少灾害损失,提高粮食种植的社会和经济效益有着十分重要的战略意义和实用价值。
本研究利用宿迁水稻暴雨洪涝灾损预报模型输出的灾损率预测结果和重灾发生乡镇比例两个因素,设定灾害等级划分阈值,实现对宿迁水稻种植区进行暴雨洪涝灾害风险等级划分的目的,并统计各风险等级的历史灾害参数特征,用以指导本地区水稻种植。
气象数据:2003~2013年宿迁境内各乡镇降水、大风及常年平均值,由宿迁市气象局提供,并经过严格的质量控制。
水稻产量数据:2003~2013年宿迁各乡镇水稻产量数据,来源于宿迁市及各县区的统计年鉴;水稻生育期数据由宿迁市农委提供。
水稻灾损数据:2008~2013年水稻灾害损失 ≥ 10%的重灾案例数据,来源于承担宿迁市农业保险业务的宿迁市平安保险公司。规定:当某乡镇出现暴雨洪涝灾害,有一个或以上自然村出现水稻灾损重灾案例时,视为该乡镇发生一次水稻暴雨洪涝灾害重灾案例。
沈伟等 [
图1. 宿迁水稻区划图
计算结果结合重灾案例发生时间分析了得到:水稻孕穗扬花期(7月21日至9月10日)和灌浆成熟期(9月10日至10月20日)是宿迁地区水稻种植受暴雨洪涝灾害威胁最大的阶段,而在分蘖拔节期及以前,暴雨洪涝灾害对宿迁水稻种植基本不会造成大的影响。
ξ = Y − Y t Y t ⋅ % (1)
董家山 [
I 2 = [ 0.4 ( L L ¯ + 1 2 ∑ R L ∑ R ¯ L ¯ + ∑ R ∑ R ¯ ) + 0.6 ( L s L s ¯ + 1 2 ∑ R s L s ∑ R s ¯ L s ¯ + ∑ R s ∑ R s ¯ ) ] ⋅ F max 8.0 (3)
(2)式和(3)式分别为孕穗扬花期和和灌浆成熟期的暴雨洪涝指数公式。 L 、 L s 、 ∑ R 、 ∑ R s 分别为孕穗扬花期(或灌浆成熟期)的降水日数、暴雨日数、总降水量、总暴雨量, L ¯ 、 L s ¯ 、 ∑ R ¯ 、 ∑ R s ¯ 为相应统计量的常年平均值, F max 为灌浆成熟期暴雨日的瞬时极大风速。
面板数据 [
y i t = μ + ∑ k = 1 n β k i x k i t + ν i t (4)
利用面板回归模型分析水稻灾损与暴雨洪涝指数的相互关系。(4)式 y i t 在文中代表某区域内乡镇i第t年水稻因暴雨洪涝造成的损失率, x k i t 表示乡镇i第t年第k个关键生育期内的暴雨洪涝灾害指数(k = 1或2,代表孕穗扬花期或灌浆成熟期), β k i 为的相关系数, ν i t 为随机变量。
Hausman检验 [
以水稻实际灾损率为研究对象,水稻孕穗扬花期和灌浆成熟期的暴雨洪涝灾害指数作为影响因子,代入模型的一般方程(4)式。通过Hausman检验,确定面板回归模型选用随机效应模型(图2,IV区的检验结果,其它各区图略),建立宿迁水稻暴雨洪涝灾害损失率预测模型(图3,IV区建模结果,其它各区图略),即各区水稻暴雨洪涝灾损灾损模型方程(表1)。将2008~2013年各乡镇的孕穗扬花期和灌浆成熟期暴雨洪涝指数带入表1中对应的模型方程,即可输出水稻暴雨洪涝灾害损失率的模型预测结果。
图2. IV区Hausman检验结果
图3. IV区建模结果
区 | 模型方程 | R |
---|---|---|
I | y = 0.648 − 0.376x1it − 0.545x2it | 0.6924 |
II | y = 1.090 − 0.428x1it − 0.509x2it | 0.7051 |
III | y = 2.176 − 1.565x1it − 0.345x2it | 0.7496 |
IV | y = 0.666 − 0.401x1it − 0.585x2it | 0.7526 |
V | y = 0.735 − 0.596x1it − 0.595x2it | 0.7711 |
VI | y = 0.702 − 0.583x1it − 0.607x2it | 0.7928 |
VII | y = 0.549 − 0.618x1it − 0.711x2it | 0.7843 |
VIII | y = 0.916 − 0.484x1it − 0.638x2it | 0.7165 |
IX | y = 0.694 − 0.575x1it − 0.498x2it | 0.7437 |
X | y = 0.821 − 0.846x1it − 0.375x2it | 0.8013 |
表1. 各区水稻暴雨洪涝灾害灾损模型方程
从表2可见,模型结果的平均绝对误差为0.07,平均相对误差为5%,说明总体上模拟效果较好。从分区情况来看,IX区、VII区、II区和VIII的效果明显好于其他各区,绝对误差均 ≤ 0.03,相对误差 ≤ 2%;IV区效果较差,相对误差达到28%。
项目 | I | II | III | IV | V | VI | VII | VIII | IX | X | 平均 |
---|---|---|---|---|---|---|---|---|---|---|---|
实际灾损率(−%) | 2.97 | 1.69 | 3.09 | 0.71 | 1.49 | 1.48 | 2.25 | 1.58 | 1.24 | 1.48 | 1.80 |
模型结果(−%) | 3.06 | 1.71 | 3.28 | 0.91 | 1.50 | 1.58 | 2.27 | 1.55 | 1.24 | 1.54 | 1.86 |
绝对误差 | 0.09 | 0.02 | 0.19 | 0.20 | 0.01 | 0.10 | 0.02 | 0.03 | 0.00 | 0.06 | 0.07 |
相对误差(%) | 3 | 1 | 6 | 28 | 1 | 7 | 1 | 2 | 0 | 4 | 5 |
表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%。
本研究进行暴雨洪涝灾害风险等级划分主要参考受灾程度和重灾面积两个参数,受灾程度用模型结果输出的暴雨洪涝灾害灾损率来表示,重灾面积用重灾乡镇比例来表示,即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可见,2008~2013年宿迁地区水稻暴雨洪涝灾害一级风险区位于沿淮河北岸和洪泽湖西岸;二级风险区位于泗洪县西部高岗地区和新沂河以北地区;三级风险区位于黄墩湖滞洪区及周边地区,四级风险区位于骆马湖以东、新沂河以南和京杭大运河沿岸及以北地区。
表5给出了根据2008~2013年水稻暴雨洪涝灾害情况统计得到的各个风险等级的灾害特征参数。从表中可见,一级风险区暴雨洪涝灾害造成的平均灾损率、重灾乡镇比例、受灾最重乡镇的平均灾损率和灾害最重年份的平均灾损率分别为−3.03%、0.61%、−5.94%和−7.46%,灾害最重年份灾情最重乡镇的平均灾损率高达到−23.27%;二级和三级风险区的各项灾害参数依次下降;四级风险区的平均灾损率、受灾最重乡镇的平均灾损率和灾害最重年份的平均灾损率均最低,分别为−1.30%、−2.58%和−3.21%,但重灾乡镇比例比三级略高,灾害最重年份灾情最重乡镇的平均灾损率比三级风险区偏高,我们认为四级风险区的乡镇数较多,相关参数更具有代表性,三级风险区的乡镇数较少,相关参数具有一定的偶然性,随着以后在应用中历史资料的进一步补充,这一缺陷将得到弥补。
项目 | I | II | III | IV | V | VI | VII | VIII | IX | X | 平均 | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
拟 | 漏报 | 2.4 | 0.0 | 0.0 | 0.0 | 3.7 | 1.0 | 6.8 | 0.0 | 7.7 | 0.0 | 2.7 |
合 | 空报 | 7.1 | 8.3 | 1.7 | 9.5 | 5.6 | 0.0 | 1.9 | 1.2 | 2.6 | 4.2 | 3.2 |
% | 准确 | 90.5 | 91.7 | 98.3 | 90.5 | 90.7 | 99.0 | 91.3 | 98.8 | 89.7 | 95.8 | 94.1 |
试 | 漏报 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 |
用 | 空报 | 14.3 | 0.0 | 0.0 | 28.6 | 0.0 | 0.0 | 0.0 | 7.1 | 0.0 | 25.0 | 4.8 |
% | 准确 | 85.7 | 100 | 100 | 71.4 | 100 | 100 | 100 | 92.9 | 100 | 75.0 | 95.2 |
表3. 2008~2013年灾损率 ≥ 10%的重灾乡镇拟合率及14年试用准确率统计
项目 | I | II | III | IV | V | VI | VII | VIII | IX | X | 平均 |
---|---|---|---|---|---|---|---|---|---|---|---|
模型灾损率(−%) | 3.06 | 1.71 | 3.28 | 0.91 | 1.50 | 1.58 | 2.27 | 1.55 | 1.24 | 1.54 | 1.86 |
重灾乡镇比例(%) | 0.54 | 0.57 | 0.68 | 0.23 | 0.27 | 0.23 | 0.37 | 0.30 | 0.40 | 0.35 | 0.39 |
YT (−) | 1.65 | 0.97 | 2.23 | 0.21 | 0.41 | 0.36 | 0.84 | 0.47 | 0.50 | 0.54 | 0.81 |
风险等级 | 一级 | 二级 | 一级 | 四级 | 四级 | 四级 | 二级 | 四级 | 四级 | 三级 | |
包含乡镇数(个) | 7 | 6 | 10 | 7 | 9 | 16 | 18 | 14 | 13 | 4 |
表4. 2008~2013年历史风险等级划分结果
图4. 2008~2013年宿迁水稻暴雨洪涝灾害风险区分布图
等级 | 实际平均 灾损率(%) | 重灾乡 镇比例 | 受灾最重乡镇 平均灾损率(%) | 最重灾年平 均灾损率(%) | 最重灾年最重乡 镇平均灾损率(%) | 含乡镇数 (个) |
---|---|---|---|---|---|---|
一级 | −3.03 | 0.61 | −5.94 | −7.46 | −23.27 | 17 |
二级 | −1.97 | 0.47 | −5.36 | −5.86 | −21.03 | 24 |
三级 | −1.48 | 0.35 | −3.65 | −4.34 | −9.33 | 4 |
四级 | −1.30 | 0.36 | −2.58 | −3.21 | −13.64 | 59 |
表5. 各风险等级历史灾害特征参数统计
项目 | I | II | III | IV | V | VI | VII | VIII | IX | X |
---|---|---|---|---|---|---|---|---|---|---|
YT | −0.0 | −0.0 | −1.41 | −0.0 | −0.0 | −0.0 | −0.13 | −0.0 | −0.0 | −0.0 |
风险等级 | 四级 | 四级 | 二级 | 四级 | 四级 | 四级 | 四级 | 四级 | 四级 | 四级 |
实际灾损率(%) | −0.8 | −0.5 | −2.4 | −0.6 | −0.7 | −0.7 | −1.1 | −0.6 | −0.6 | −0.2 |
实际重灾乡镇数,比例 | 0,0 | 0,0 | 4,0.57 | 0,0 | 0,0 | 0,0 | 1,0.06 | 0,0 | 0,0 | 0,0 |
表6. 2015年应用结果统计
用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%。
利用降水(暴雨)强度、降水(暴雨)日数和累计降水(暴雨)量参数自定义的暴雨洪涝灾害指数,能较好的描述暴雨洪涝灾害特征,与水稻暴雨洪涝灾害损失率的相关性较好。利用时间和界面序列二维面板数据建立随机效应回归预测模型,在一定程度上克服了资料年代短的缺陷,模型效果较好。
直接应用模型输出的预测结果进行水稻暴雨洪涝灾害风险等级划分,减少了灾情调查环节,节省了大量的时间和人力物力,经检验和使用效果较好,能够反映出宿迁地区水稻暴雨洪涝灾害的相关特征,具有较高的实用价值。从历史风险区的划分结果来看:一级风险区位于沿淮河北岸和洪泽湖西岸;二级风险区位于泗洪县西部高岗地区和新沂河以北地区;三级风险区位于黄墩湖滞洪区及周边地区,四级风险区位于骆马湖以东、新沂河以南和京杭大运河沿岸及以北地区,与宿迁地区水稻暴雨洪涝灾害的实际情况相符合。
2008~2013年四级风险区的重灾乡镇比例和灾害最重年份灾情最重乡镇的平均灾损率比三级风险区偏高,这可能是因为三级风险区的乡镇数较少,相关参数具有一定的偶然性造成的,随着历史资料的进一步补充,这一缺陷将得到弥补。
江苏省气象科研开放基金项目(项目编号:K201010)。
沈 伟,吴新胜,周 航,王文清,张 莹. 基于灾损预测模型的水稻洪灾风险等级划分The Risk Rank of Rice Flood Disaster Based on Damage Prediction Model[J]. 农业科学, 2017, 07(07): 530-538. http://dx.doi.org/10.12677/HJAS.2017.77069