中长期水文预报统计预报方法.ppt

中长期水文预报统计预报方法 * * 四、应用示例简介（多元回归） （2）以丹江口以上流域径流量为例 率定期（ 1984-2009 ）：年降雨总量误差10.8%，汛期各月平均误差约为15%~60%，年总量预报精度最高，枯季及汛期月份次之。 验证期（2010-2013）：年降雨总量误差13.5%，汛期各月平均误差约为14%~50%，年总量预报精度最高，枯季及汛期月份次之。 预报因子 环流特征量 提前年份 月份 x1 登陆台风 1 8 x2 北非副高脊线(20W-60E) 1 8 x3 北非副高北界(20W-60E) 1 12 x4 北半球极涡中心位置(JW) 1 6 x5 北非副高脊线(20W-60E) 1 6 x6 东亚槽位置(CW) 1 12 x7 亚洲经向环流指数(IM,60E-150E) 1 12 x8 东太平洋副高脊线(175W-115W) 1 6 x9 冷空气 1 4 x10 北美大西洋副高北界(110W-20W) 1 11 四、应用示例简介（多元回归） （2）以丹江口以上流域径流量为例 2010年 2011年 2012年 2013年 阶数 P=4 四、应用示例简介（自回归-AR(P)） （3）以丹江口以上流域4月径流量为例 率定期：1955-2009 验证期：2010-2013 AR(4) 方程 验证期误差统计 -9 29 38 18 -23 误差（%） 年份 合格？ 2010 合格 2011 合格 2012 不合格 2013 合格 2014 合格 四、应用示例简介（自回归-AR(P)） 预报结果示意图 （3）以丹江口以上流域4月径流量为例 四、应用示例简介（马尔科夫链） （4）以丹江口水库9月报10月入库径流量为例 率定期：1951-2008 验证期：2009-2013 为fi,j为第i状态经一步转移为第j状态的频数，转移概率为 状态划分 转移矩阵计算 四、应用示例简介（马尔科夫链） （4）以丹江口水库9月报10月入库径流量为例 预测年份 实际水平年 一步转移概率 枯水年 平水年 丰水年 2009 平水 0.024 0.878 0.098 2010 平水 0.024 0.878 0.098 2011 平水 0.000 0.700 0.300 2012 平水 0.024 0.878 0.098 2013 平水 0.024 0.878 0.098 预测水平年与实际情况相符 9月报10月径流预报成果表 四、应用示例简介（谐波分析） 则Xt的周期方程为： （5）以丹江口水库1月径流量为例 率定期：1951-2009 验证期：2010-2013 置信度： α=0.05 显著周期： 10、5年 周期成分1 周期成分2 采用上述方程，模拟1951~2009年1月径流量，其平均相对误差为24.4%；验证2010~2013年1月径流量，其平均相对误差为27.2%。 15.7% 419 362 2013 34.1% 495 369 2011 -28.5% 368 515 2012 -30.4% 293 421 2010 相对误差 模拟值 实测值 年份 相对误差 模拟值 实测值 年份 四、应用示例简介（谐波分析） 四、应用示例简介（随机森林） （6）丹江口逐月径流预报 预报因子：10项环流因子　　预报变量：逐月径流量 ··· 25棵决策树 率定期：1950~2009 验证期：1950~2009 读取训练数据 25 四、应用示例简介（随机森林） 验证期逐月预报结果： 四、应用示例简介（统计分析与水文模型耦合） （7）丹江口水库入库径流预报——多元回归预测+SWAT 预报月降雨过程 日降雨过程 相似性原理 SWAT 水文模型 预报日/月径流过程 模型参数率定与验证 率定期：1995-2005 验证期：2006-2011 多元线性回归模型 2012年月降雨预报 预报结果分析 选择预报因子 确定回归系数 模型构建与验证 最小距离原则: 月降雨统计预报模型（多元线性回归方程） 预报因子 环流特征量 提前年数 月份 x1 北美副高脊线(110W-60W) 1 11 x2 北美副高北界(110W-60W) 1 11 x3 太平洋区极涡面积指数(2区,150E-120W) 1 1 x4 南海副高北界(100E-120E) 2 12 x5 编号台风 1 5 x6 亚洲经向环流指数(IM,60E-150E) 1 5 x7 南海副高脊线(100E-120E) 1 2 x8 冷空气 1 3 x9 北美大西洋副高脊线(110W-20W) 2 11 x10 北半球极涡中心位置(JW) 1 7 预报因子选择 回归方程 最小二乘法计算回归系数 以汉江上游区 1 月份的月降雨为例 模型验证 月份 20