一种耦合长、中、短期径流预报信息的水库优化调度方法技术

本发明专利技术提供一种耦合长、中、短期径流预报信息的水库优化调度方法,该方法包括以下步骤：1)分析径流自身随机性及其预报不确定性，并进行描述；2)构建耦合长、中、短期径流预报信息的水库优化调度模型，确定其目标函数、约束条件及递推方程；3)结合递推方程对所建立的水库优化调度模型进行求解。本发明专利技术的效果和益处是将长、中、短三种预见期径流预报信息耦合至水电调度中，对不确定性描述更为合理，对预报径流的处理方式更为有效，有利于提高水库调度对不确定信息的使用效率，保证水库发电效益和系统稳定性。

【技术实现步骤摘要】

本专利技术属于水电调度运行
，涉及一种耦合长、中、短期径流预报信息的水库优化调度方法。
技术介绍
近年来，随着水文预报技术的发展，预报、调度逐渐走向融合。耦合降雨和径流等综合预报信息的水库调度有利于延长预见期、提高水能资源利用率，而预报信息受预报精度、下垫面等诸多因素影响，存在着较大的不确定性，故考虑预报信息的水库调度决策受预报不确定性的影响大。如何合理地将预报信息与水库调度进行耦合以减少预报不确定性的影响，成为研究者及调度人员亟待解决的问题。随机动态规划(SDP)能够将水库蓄水、径流等状态及不确定性考虑到算法中，与水库多阶段决策的特点具有较高的耦合性。近年来，越来越多的研究开始尝试将径流预报信息及其不确定性耦合到SDP模型中。其中单一预见期的径流预报信息耦合使用方式的研究较多，考虑不同预见期的径流预报信息的耦合使用方式的研究也相继开展起来，如徐炜等(XuW.,ZhangC.,PengY.,etal.WaterResourcesResearch[J].2014,50(12):9267-9286.)使用短期和中期径流预报信息作为水文状态变量，建立了短、中期径流预报套接的贝叶斯随机动态规划模型(TS-BSDP)；唐国磊(唐国磊.考虑径流预报及其不确定性的水电站水库调度研究.大连理工大学,2009.)构建了考虑中、长期径流预报信息的随机动态规划模型。SDP模型具有使用的水文状态变量越多越能描述入流的随机特性的特点(Tejada-GuibertJ.A.,JohnsonS.A.,StedingerJ.R.Waterresourcesresearch[J].1995,31(10):2571-2579)，当前研究中仅使用了中、长期或短、中期径流预报信息，尚未将长、中、短期径流预报信息同时作为水文状态变量耦合至随机动态规划模型中。如何有效地处理长、中、短期径流预报信息，并耦合至水库调度中，亟需做进一步的研究。
技术实现思路
针对现有技术的不足，本专利技术提供一种耦合长、中、短期径流预报信息的水库优化调度方法。本专利技术的技术方案具体如下：一种耦合长、中、短期径流预报信息的水库优化调度方法，具体包括以下步骤：第一步，分析、描述径流及径流预报的不确定性；所述的径流及径流预报的不确定性针对长、中、短期径流预报信息。将欧洲气象中心(ECMWF)提供的未来1-10天集合降雨预报信息平均值输入到新安江模型得到10天径流预报信息，并划分为1～5天(短期)和6～10天(中期)两段，长期径流预报信息是指利用相关系数法从74项大气环流因子选取相关性较高的因子并结合ANN模型得到的6～10天所对应时段至汛末总来水量的预报值。1.1径流自身随机性由于径流表现出一定的随机性，常把径流作为随机过程来处理。假定水库入库流量具有简单马尔科夫特性，采用径流状态转移概率P[Qt|Qt-1]描述径流随机性，反映径流本身的不确定性，其中，P[Qt|Qt-1]表示当时段t-1流量为Qt-1条件下时段t发生Qt的概率，Q表示流量。1.2径流预报的不确定性反映了径流预报的不确定性，可用预报流量Ft与实测值间Qt的统计概率P[Ft|Qt](称为似然概率)表示。1.3后验的径流状态转移概率由于径流状态转移概率P[Qt|Qt-1]是在未获得任何径流预报信息之前已经获得，常将其称为先验的径流状态转移概率。一旦获得新的预报信息，便可由统计方法或贝叶斯定理得到后验的径流状态转移概率。当时段t-1径流量为Qt-1和t时段预报流量为Ft时，t时段径流实际发生在Qt的后验的径流状态转移概率可综合先验的径流状态转移概率和似然概率由贝叶斯定理求得，如下所示：P[Qt|Qt-1,Ft]=P[Ft|Qt]×P[Qt|Qt-1]ΣQtP[Ft|Qt]×[Qt|Qt-1]---(1)]]>当时段t-1径流量为Qt-1、t时段中期预报流量表示为Ft、t时段长期预报流量为FLt，t时段径流实际发生在Qt的后验概率表示为P[Qt|Qt-1,Ft,FLt]，可通过统计方法求得。1.4径流的可预测概率当时段t-1径流量为Qt-1和t时段预报流量为Ft时，径流的可预测概率为P[Ft|Qt-1]，可结合已知的先验概率P[Qt|Qt-1]和似然概率P[Ft|Qt]，通过全概率公式计算求得P[Ft|Qt-1]=ΣQt+1P[Ft|Qt]×P[Qt|Qt-1]---(2)]]>当时段t-1径流量为Qt-1、t时段中期预报流量表示为Ft、t时段至汛末总的长期预报流量表示为FLt，径流的可预测概率为P[Ft|Qt-1,FLt]·P[FLt|Qt-1]，可通过统计方法求得。第二步，构建耦合长、中、短期径流预报信息的水库优化调度模型(SML-BSDP模型)，确定其目标函数、约束条件及递推方程；2.1目标函数水库发电优化调度模型以调度期在满足发电保证出力的条件下，总发电量最大为目标函数，计算公式如下：Max[Σt=1TE[B(Zt,Ht,Zt+1)]]---(3)]]>B(Zt,Ht,Zt+1)＝[b(Zt,Ht,Zt+1)-α{Max(e-b(Zt,Ht,Zt+1),0)本文档来自技高网...

【技术保护点】
一种耦合长、中、短期径流预报信息的水库优化调度方法，其特征在于，包括以下步骤：第一步，针对长、中、短期径流预报信息，分析、描述径流及径流预报的不确定性；1.1)采用径流状态转移概率P[Qt|Qt‑1]描述自身径流随机性，反映径流本身的不确定性，其中，P[Qt|Qt‑1]表示当时段t‑1流量为Qt‑1条件下时段t发生Qt的概率，Q表示流量；所述的流状态转移概率P[Qt|Qt‑1]为先验的径流状态转移概率；1.2)用预报流量Ft与实测值间Qt的似然概率P[Ft|Qt]表示径流预报的不确定性；1.3)将获得的新的预报信息，用统计方法或贝叶斯定理得到后验的径流状态转移概率；当时段t‑1径流量为Qt‑1和t时段预报流量为Ft时，t时段径流实际发生在Qt的后验的径流状态转移概率通过综合先验的径流状态转移概率和似然概率由贝叶斯定理求得，如下所示：P[Qt|Qt-1,Ft]=P[Ft|Qt]×P[Qt|Qt-1]ΣQtP[Ft|Qt]×P[Qt|Qt-1]---(1)]]>当时段t‑1径流量为Qt‑1、t时段中期预报流量表示为Ft、t时段长期预报流量为FLt，t时段径流实际发生在Qt的后验概率表示为P[Qt|Qt‑1,Ft,FLt]，通过统计方法求得；1.4)径流的可预测概率当时段t‑1径流量为Qt‑1和t时段预报流量为Ft时，径流的可预测概率P[Ft|Qt‑1]为：P[Ft|Qt-1]=ΣQt+1P[Ft|Qt]×P[Qt|Qt-1]---(2)]]>当时段t‑1径流量为Qt‑1、t时段中期预报流量表示为Ft、t时段至汛末总的长期预报流量表示为FLt，径流的可预测概率为P[Ft|Qt‑1,FLt]·P[FLt|Qt‑1]，通过统计方法求得；第二步，构建耦合长、中、短期径流预报信息的水库优化调度模型，确定其目标函数、约束条件及递推方程；2.1)目标函数水库发电优化调度模型以调度期在满足发电保证出力的条件下，总发电量最大为目标函数，计算公式如下：Max[Σi=1TE[B(Zt,Ht,Zt+1)]]---(3)]]>B(Zt,Ht,Zt+1)＝[b(Zt,Ht,Zt+1)‑α{Max(e‑b(Zt,Ht,Zt+1),0)}β]·Δt       (4)式中，Zt+1为时段t的最优决策末水位；Zt为时段t的初水位；Ht为时段t的水文状态变量；B(·)是时段t的发电量效益函数，即t时段惩罚后的发电量；b(·)是t时段水电站的发电出力，由耗水率、水头等计算获得；α和β是惩罚系数，由水电站发电保证率确定；e为水库水电站发电保证出力；Δt是计算时段长度；E[B(·)]为时段t在给定Zt和Ht下的期望发电量；T为调度期的总时段数；2.2)约束条件在水库发电优化调度模型的优化过程中，水库的蓄水、泄流及发电出力需要满足约束条件；2.3)递推方程水库发电优化调度模型的递推方程在非汛期考虑短、中期径流预报信息，以Qt、Ft+1为水文状态变量，递推方程为：foptnp5(Zt,Qt,Ft+1)=Max{Bt(Zt,Qt,Zt+1)+ΣQt+1PQt+1|Qt,Ft+1·ΣFt+1PFt+2|Qt+1·foptn-1(Zt+1,Qt+1,Ft+2)}foptna5(Zt,Qt,Ft+1)=Max{PQt+1|Qt,Ft+1·ΣFt+2PFt+2|Qt+1·foptn-1(Zt+1,Qt+1,Ft+2)}foptnt10(Zt,Qt,Ft+1)=QtQt+Ft+1foptnp5(Zt,Qt,Ft+1)+Ft+1Qt+Ft+1foptna5(Zt,Qt,Ft+1)---(10)]]>水库发电优化调度模型的递推方程在汛期考虑长、中、短期径流预报信息，以Qt、Ft+1和长期预报信息FLt+1为水文状态变量，基于这三种变量的信息更新下时段流量Qt+1的发生概率，递推方程为：在汛期和非汛期的过渡期，当t时刻和t+1时刻有长期预报信息，而t+2时刻无长期预报信息时，递推方程为：foptnp5(Zt,Qt,Ft+1,FLt+1)=Max{Bt(Zt,Qt,Zt+1)+ΣQt+1PQt+1|Qt,Ft+1,FLt+1·ΣFt+2PFt+2|Qt+1·foptn-1(Zt+1,Qt+1,Ft+2)}foptna5(Zt,Qt,Ft+1,FLt+1...

【技术特征摘要】
1.一种耦合长、中、短期径流预报信息的水库优化调度方法，其特征在
于，包括以下步骤：
第一步，针对长、中、短期径流预报信息，分析、描述径流及径流预
报的不确定性；
1.1)采用径流状态转移概率P[Qt|Qt-1]描述自身径流随机性，反映径流
本身的不确定性，其中，P[Qt|Qt-1]表示当时段t-1流量为Qt-1条件下时段t
发生Qt的概率，Q表示流量；所述的流状态转移概率P[Qt|Qt-1]为先验的径
流状态转移概率；
1.2)用预报流量Ft与实测值间Qt的似然概率P[Ft|Qt]表示径流预报的不
确定性；
1.3)将获得的新的预报信息，用统计方法或贝叶斯定理得到后验的径
流状态转移概率；
当时段t-1径流量为Qt-1和t时段预报流量为Ft时，t时段径流实际发
生在Qt的后验的径流状态转移概率通过综合先验的径流状态转移概率和似
然概率由贝叶斯定理求得，如下所示：
P[Qt|Qt-1,Ft]=P[Ft|Qt]×P[Qt|Qt-1]ΣQtP[Ft|Qt]×P[Qt|Qt-1]---(1)]]>当时段t-1径流量为Qt-1、t时段中期...

【专利技术属性】
技术研发人员：彭勇，张小丽，张弛，王本德，唐榕，黄朔，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21