变化环境下的水库调度

变化环境下的水库调度——新挑战与新方法赵建世清华大学2012.7.21西安内容提要1.变化环境带来的理论挑战2.基于风险对冲规则的水库调度理论3.适应性动态调度框架4.密云水库案例分析5.结论水资源系统分析基本方法现代水资源系统分析理论与方法体系起源于1950年代的“哈佛水项目”（HarvardWaterProgram）水库调度是水资源系统分析的重要内容之一，其核心问题是确定“何时放水”以及“放多少水”（Loucks，2005）优化和模拟式水资源系统分析的两个基本工具，也是传统水库调度研究的基本手段。（Labadie,M.ASCE，2004）调度曲线基于历史水文信息和人类需求，通过优化和模拟研究，可以设计出实践中常用的调度曲线（Rulecurves），指导水库的调度管理工作。水库调度曲线代表了确定的年内的调度规则，其一经确立，一般不会随着未来的径流情势进行改变。由于曲线规则本身基于对历史水文信息的统计表达，因此很难与实时的水文预报信息相结合（Loucks，2005）。这些基本特性决定了其难以适应不断变化的环境。传统调度模型的一般性假定Stationarity（稳态性）径流概率分布随时间不变(Labadie,M.ASCE，2004)调度规则曲线（RuleCurve）过去曾经是打开未来之门的钥匙，但现在门锁已经换了！变化环境的挑战径流变化：非稳态性（Nonstationarity）密云水库入库径流1956-2009变化环境的挑战人类需求的变化密云水库:设计：防洪为主-灌溉为主(1960s-1980s)-城市供水(1990s-2008)-安全储备(2008-)重新思考水库调度问题的本质Humanshavechangedthewaytheworldworks.Nowtheyhavetochangethewaytheythinkaboutit,too.——TheEconomist,May26,2011水库:一个典型的人类-自然耦合系统（CoupledHuman-Naturesystem）人类系统自然系统经济学水文学交互作用变化环境下的水库调度耦合经济学原理和水文学成果进行适应性的动态调度Hydro-Economics内容提要1.变化环境带来的理论挑战2.基于风险对冲规则的水库调度理论3.适应性动态调度框架4.密云水库案例分析5.结论调度目标的经济特性与调度策略供水量效益值当调度目标为线性函数时（如缺水量最小），标准调度策略（SOP）是最优策略（DraperandLund，2004；YouandCai，2008）。标准调度策略（SOP）优先满足当前阶段用水需求，当前阶段需求完全满足后，水库开始蓄水以备未来使用，直至水库蓄满。调度目标的经济特性与调度策略供水量效益值风险对冲规则（HedgingRulePolicy,HRP）是Bower等人（1962）在哈佛水项目中提出的，其出发点是引用经济学理论来探讨水资源稀缺和不确定的条件下，水量在时间上的最有分配准则。DraperandLund[2004]：“在最优解处，存水的边际效益应该等于供水的边际效益”YouandCai[2008]：“当两个阶段的供水具有相同的边际效益时，两个阶段的总效益达到最大”供水量效益值边际效益：最后一滴水的价值；经济学中的边际效益递减原理；其数学表达为效用函数的一阶导数B‘(x).风险对冲理论模型DraperandLund（2004）提出了风险对冲规则的一般性经济学解释()()..00mMaxzBDCSstSDASSkDDd()()BDCSDS拉格朗日定理“在最优解处，存水的边际效益应该等于供水的边际效益”如何考虑物理约束的影响？如何考虑不确定性的影响？两阶段模型000000..)()(maxarg2121min1max1222111112211llxxKsKsSlxIsslxIStsxBxBei)81()71()61()51()41()31()21()11(KKT最优性条件如果存在常数和一组可行解能够满足上述方程，那么是两阶段水库调度模型的全局最优解.0,,,,,000)(0)(0)(0)(0)()(2121*22*11*1minmax*12*1*221*1*112122*2'211*1'1dcnnmmnndcemimdcmmnmnmxxsKKsIsSxISsxxBxB)106()96()86()76()66()56()46()36()26()16(0,,,,,2121dcnnmm),,(*1*2*1sxx),,(*1*2*1sxx最优性条件的经济学解释在满足所有物理约束的条件下，让当前阶段和未来阶段用水的边际效益尽量接近。2*2'21*1'1)()(ndncxBxB水文预报及其不确定性在理解水库调度经济学原理的基础上，如何耦合水文预报成果？如何考虑预报的不确定性？水文预报中的误差累计过程(Viner,2002))|,(PEPfI水文预报水资源管理Zhao,J.,X.Cai,andZ.Wang(2011),Optimalityconditionsforatwo-stagereservoiroperationproblem,WaterResourcesResearch,47,W08503,doi:10.1029/2010WR009971.考虑水文预报及其不确定性的对冲规则考虑一般性的水文预报及其不确定性,Zhaoet.al.[2011]定义了“风险调整后的边际效益”（Risk-adjustedmarginalbenefit,RAMB)严格的数学推导可以证明，耦合预报及其不确定性的最优性准则为：“在满足所有物理约束的条件下，让当前和未来阶段用水的风险调整后的边际效益尽量接近。”2'''')(5.0)()(iiiiiiixBxBxRAMB2*21*1)()(ndncxRAMBxRAMBmaxE[Bi(xi)i=1nå]s.tSini+(I1+e1)-x1-s1-l1=0s1+(I2+e2)-x2-s2-l2=0s2+(I3+e3)-x3-s3-l3=0......sn-1+(In+en)-xn-Send-ln=0xi³0li³0si³Kminsi£Kmax(1)1maxmin0()0()0()0()0,0,0iminimimicidiiiiimiiniiciciidiminicidiRAMBxssIxgsKKsRAMBi=¢Bi(xi)+0.5¢¢¢Bi(xi*)×si2KKTconditionsRule1:最理想的情况是各阶段的边际效益相同。即1iiRAMBRAMBRule2:如果某个阶段的流量非负约束被违反，则该阶段的最优放水量为0.11iiniiRAMBRAMBRAMBRule3:A如果某个阶段的最大或者最小库容约束被违反，则该阶段是一个RAMB的变化点。or1iciiRAMBRAMB1iidiRAMBRAMB最优准则:在满足物理约束的条件下，让各阶段的风险调整后的边际效益尽量接近。从两阶段模型到多阶段模型基于风险对冲规则的优化调度算法（Hedgingrulebasedalgorithm，HRBA)Step3(Rule3):If检查solution_1的各阶段放水量，如果某些阶段违反了最大或者最小库容约束，挑选出违反程度最大的阶段并把这个阶段的库容设置在边界上（最大或者最下库容），这样可以将原问题分解为两个子问题，然后重复上述步骤，直至得到最优解.Step1(Rule1):不考虑约束，寻找各阶段具有相同RAMB的理想解，定义为solution_0Step2(Rule2):检查solution_0各阶段的放水量，如果某些阶段出现负放水量，则将这些阶段放水量设置为0，然后重新计算其他阶段的防水量，得到solution_1。HRBA的效率评估算法变量离散总计算时间(秒)K=1K=2K=3K=4K=5HRBA1004.022.151.931.681.47100018.6214.1410.828.997.16DP10020.0219.6719.5719.4419.2510001968.691944.741923.761910.661902.61•入流:100个人工随机生成径流序列(每个序列100个阶段).•效用函数：32924iiixxx计算时间与变量数目为线性关系计算时间与变量数目为二次方关系内容提要1.变化环境带来的理论挑战2.基于风险对冲规则的水库调度理论3.适应性动态调度框架4.密云水库案例分析5.结论面向实时的动态决策过程Rollinghorizondecisionprocess（RHDP)面向实时的动态决策过程，是水库管理决策实际过程一种抽象，其核心是所有供决策使用的信息与管理实践完全相同。这一框架下的优化分析对实时调度管理具有直接的指导意义。预测的密云未来径流（预报年到2009）非稳态径流预报模型2009/(20091)tmeaniitIIt1-fii=14åBiæèçöø÷[Xt-Xt-1-m]=et年度径流预报模型:ARIMA(4,1,0)“4”代表自回归阶数；“1”代表线性趋势；“0”为移动平均阶数。observedpredicted内容提要1.变化环境带来的理论挑战2.基于风险对冲规则的水库调度理论3.适应性动态调度框架4.密云水库案例分析5.结论Ma,et.al,2010密云水库天然径流变化天然径流2003年之前采用实测数据；2003年（含）之后采用：可用水量+蒸发量（0.824）-补水量补水情况资料来源：北京水资源公报人类直接取水Ma,et.al,2010密云水库案例分析Bi(x)=2x3-114x2+1680x(1+r)i-1Kmin=4.37´108m3Kmax=39.64´108m3Sini=28.8´108m3Send=10.37´108m3调度策略定义SOP标准调度策略，尽量满足当前最大需水（10亿方/年）HR-1基于ARIMA(4,1,0)年水量预测模型和风险对冲的适应性动态调度框架HR-0基于完美预报的理想调度策略（或已知未来入流信息）AO1996-2009的实际调度策略效益函数效益函数的拟定：（1）最大总需水量为10亿方，且过原点。（2）将不同供水的优先级对应于不同边际效益生活用水：最高优先级，对应于缺水时最高的边际效益；工业用水：第二优先级，对应正常水量时的较高边际效益；基本河湖用水：第三优先级，对应较丰水时普通的边际效益；奢侈河湖用水：第三优先级，对应非常丰水时较低边际效益；Bi(x)=2x3-114x2+1680x(1+r)i-1模型拟合与预报情况（从起始年份到09年）预报未来总水量起始年份实际可用水量实际天然径流预报来流量20088.106.468.67200710.898.5212.00200615.4511.9415.80200520.9816.3517.26200426.1520.5016.58200329.4323.1324.69200231.0124.7120.37200135.7329.4332.05200037.2030.9047.90199939.3033.0089.64199849.9743.6778.05199755.5149.21101.06199667.8161.5179.90结果对比SOP:尽可能满足当前需求，总效益最低.AO:不是完全的SOP，隐含考虑了缺水的风险，总效益较SOP高.HR-0:完美预报信