[发明专利]一种基于聚合水库蓄放水模拟的洪水预报方法有效
申请号: | 201610286165.9 | 申请日: | 2016-05-03 |
公开(公告)号: | CN105760703B | 公开(公告)日: | 2018-04-10 |
发明(设计)人: | 彭勇;张弛;周惠成;孙新国;张小丽;李昱 | 申请(专利权)人: | 大连理工大学 |
主分类号: | G06F19/00 | 分类号: | G06F19/00 |
代理公司: | 大连理工大学专利中心21200 | 代理人: | 梅洪玉,潘迅 |
地址: | 11602*** | 国省代码: | 辽宁;21 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 聚合 水库 放水 模拟 洪水 预报 方法 | ||
1.一种基于聚合水库蓄放水模拟的洪水预报方法,其特征在于,包括以下步骤:
第一步,利用天然期场次洪水信息率定TOPMODEL模型参数;所述的天然期为未建水利工程或水利工程很少的时期;
第二步,将流域内众多水利工程聚合成一个聚合水库,计算聚合水库的库容V,确定聚合水库的初始蓄满率δ0:
δ0=Bχ0τ+(1-B)ρ0(1)
其中,P7为洪水发生前7天总降雨量,χ0为洪水过程中某一时刻土壤饱和程度,参数B、Pmax、Pmin、τ都通过优化得到,0<δ0<1;
第三步,根据水利工程的拦蓄、放流规律制定聚合水库的蓄放水模拟图,确定水利工程影响量;所述的聚合水库的蓄放水模拟图能够表示一场洪水的水库蓄满率的变化过程、时段净雨放流率随水库蓄满率的变化过程和时段净雨拦蓄率随水库蓄满率的变化过程;
所述的时段净雨放流率随水库蓄满率的变化过程对应的曲线为k1=[(A-δ)/A]m;所述的时段净雨拦蓄率随水库蓄满率的变化过程对应的曲线为k2=[(δ-A)/(1-A)]n;水利工程影响量为:
(1)当δt<A时,时段拦蓄量△R1,t=k1,t Rt;
(2)当δt>A时,时段放流量△R2,t=k2,t Rt;
其中,参数A、m、n通过优化得到;k1,t和k2,t由聚合水库蓄放模拟图得到,k1,t为放水率,k2,t为拦蓄率;Rt为时段产流量;△R1,t、△R2,t为时段水利工程影响量;
第四步,将第一步率定的TOPMODEL模型参数作为模拟无水利工程影响或水利工程影响较小情况下的流域产流情况的基础参数;根据实测场次洪水信息,结合考虑水利工程影响的洪水计算过程,率定出洪水发生前聚合水库蓄水量参数B、Pmax、Pmin、τ和聚合水库的蓄放水模拟图参数A、m、n,进行洪水预报;
所述的考虑水利工程影响的洪水计算为:
4.1)由第一步率定的TOPMODEL模型参数得到天然预报结果Rt;
4.2)根据第二步得到的聚合水库的初始蓄满率δ0和库容V,计算聚合水库的初始库容v0=δ0V;
4.3)计算各时段水利工程影响拦蓄量△R1,t或放流量△R2,t;
初始化时段t=1,根据δ0和第三步的聚合水库蓄放模拟图得到时段水利工程影响拦蓄量△R1,1或放流量△R2,1,并计算聚合水库的该时段末库容v1及水库的该时段蓄满率δ1=v1/V,其中v1=v0+△R1,1或v1=v0-△R2,1;
令时段t=t+1,根据δt-1和第三步的聚合水库蓄放模拟图得到时段径流放水率k1,t或拦蓄率k2,t,从而得到时段水利工程影响拦蓄量△R1,t或放流量△R2,t,并计算聚合水库的时段末库容vt及水库的时段蓄满率δt=vt/V,其中vt=vt-1+△R1,t或vt=vt-1-△R2,t;
4.4)根据步骤4.3)所得的水利工程影响拦蓄量△R1,t或放流量△R2,t校正Rt,得到水利工程影响下的时段产流量R’t;
4.5)根据步骤4.4)得到的时段产流量R’t,得到水利工程影响下的场次洪水模拟结果Q’t。
2.根据权利要求1所述的一种基于聚合水库蓄放水模拟的洪水预报方法,其特征在于,所述的第二步计算聚合水库库容V步骤为:计算流域中大、中小型水库及塘坝水利工程各自聚合时的库容;将流域中大、中小水库及塘坝聚合成一个聚合水库,以大、中小型水库及塘坝各自聚合时的库容之和作为该聚合水库的聚合库容;大、中小型水库及塘坝水利工程各自聚合时的库容,分别为:
(1)大型水库:大型水库的兴利库容为聚合时的库容;
(2)中小型水库:分析水库现有运行资料得出实际运行中库容占总库容的比例,并以该比例与中小型水库总库容的乘积作为聚合时的库容;
(3)塘坝:塘坝的总库容作为聚合时的库容。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于大连理工大学,未经大连理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201610286165.9/1.html,转载请声明来源钻瓜专利网。
- 上一篇:软件篡改识别方法和系统
- 下一篇:一种海表盐度反演方法和装置
- 同类专利
- 专利分类
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用