CN104281780B - 线性水库滞留汇流及嵌套流域(多子流域)汇流方法 - Google Patents
线性水库滞留汇流及嵌套流域(多子流域)汇流方法 Download PDFInfo
- Publication number
- CN104281780B CN104281780B CN201410536584.4A CN201410536584A CN104281780B CN 104281780 B CN104281780 B CN 104281780B CN 201410536584 A CN201410536584 A CN 201410536584A CN 104281780 B CN104281780 B CN 104281780B
- Authority
- CN
- China
- Prior art keywords
- basin
- period
- runoff
- watershed
- nash
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种线性水库滞留汇流及嵌套流域(多子流域)汇流方法,收集研究流域降雨量、蒸发量和径流量的长系列数据,由降雨量、蒸发量和径流量应用水文模型计算产流过程,应用Nash线性水库方法,将流域内汇流路径视为<i>N</i>个串联的线性水库,上一个水库的出流为下一个水库的入流,得到Nash瞬时单位线,引入滞留系数()描述线性水库的蓄水能力,采用二项式近似模拟Nash瞬时单位线,对水文模型的汇流方法进行改进;本发明保留了瞬时单位线从时间和空间尺度上描述流域汇流过程的优势,同时避免了应用中复杂的参数问题,适用于在大尺度流域水文模拟,尤其在气候变化下流域水循环研究中具有重要的实用价值,该方法应用简便,且给出的模拟结果也更为客观科学、精度更高。
Description
技术领域
本发明涉及水文应用领域,具体涉及一种线性水库汇流方法。
背景技术
流域水文模型是模拟流域水文过程和认识流域水文规律的重要理论基础,对流域产汇流计算、洪水分析与预报以及水资源优化配置与调度等具有重要意义,一直是水文学和水资源等流域的重要基础研究课题之一,也是研究水文自然规律和解决水文实践问题的主要工具。
汇流过程是水文模型的重要组成部分,汇流模型的繁简取决于模型的时间分辨率和汇流滞时的相对关系。水文模型的时间尺度一般是日或者时,汇流过程对流域出口断面的流量尤其是峰值形成是非常重要的。对于尺度很小的流域,汇流过程相对比较快,对汇流过程的不同假设差别不大;对于大尺度流域,由于汇流的过程相对缓慢,采用不同的汇流方法可能影响较大。汇流过程一般包括坡面汇流和河道汇流。汇流演算模型可分为两类:水文学模型和水力学模型。水文学方法的控制方程一般只是时间的函数,水力学方法同时考虑了时间因素和空间因素。其中,马斯京根(Muskingum)及其改进方法(经典运动波方程的有限差分表示)、线性水库调蓄方法是典型的水文学汇流方法;水力学方法基本上都是基于对圣维南方程组的简化和近似,其中最主要的近似方法为扩散波和运动波。
等流时线是一种经典的流域汇流曲线,它从物理角度揭示了流域水文系统是一个有“忆滞”功能的系统。该方法经验性较强,一般假设流速在整个流域均匀分布,汇流时间主要由净雨到达流域出口断面的距离确定。该方法操作简单,在水文模型中广泛使用,但汇流流速是沿程变化的,且汇流路径并不是净雨到达流域出口断面的真实路径。
瞬时单位线汇流方法是把降雨产生的洪水过程概化为净雨经历多个相同的线性水库串联调节之后,在出口断面形成的流量过程。在实际应用时,瞬时单位线往往借助S曲线转换为时段单位线。但时段单位线将流域上呈现为分布状态的净雨集总成单一过程,对流域汇流过程的描述过于简化。
已有很多研究证明,气候变化对流域的水循环过程已经产生了影响,气候变化背景下流域的水文模拟对水文模型提出了新的挑战,亦是水文学及水资源专业领域研究的热点和难点问题。考虑到气候变化的尺度,气候变化对流域水文水资源影响的研究中,流域的尺度一般都比较大(一般为万平方公里以上),汇流过程对流域出口断面径流的模拟精度影响较大。上述传统方法要么公式本身比较复杂(如圣维南方程和瞬时单位线)在大尺度流域中难以直接应用,要么对流域汇流的描述过于简化缺乏客观表征,影响大尺度流域模拟的精度(等流时线法)。
发明内容
发明目的:本发明的目的在于针对现有技术的不足,提供一种线性水库滞留汇流及嵌套流域(多子流域)汇流方法,在保存瞬时单位线优势的同时可实现较方便、简单直接的应用,在大尺度流域水文模拟中具有重要的实用价值,可更好地为气候变化下流域水循环研究提供科学支撑。
技术方案:本发明提供了一种线性水库滞留汇流及嵌套流域(多子流域)汇流方法,包括以下步骤:
(1)收集研究流域降雨量、蒸发量和径流量的长系列数据,分析数据的代表性和一致性;
(2)选择研究流域的水文模型,由降雨量、蒸发量和径流量应用水文模型计算产流过程;
(3)应用Nash线性水库方法,将流域内汇流路径视为N个串联的线性水库,上一个水库的出流为下一个水库的入流,得到Nash瞬时单位线;
(4)由于Nash瞬时单位线难以直接应用,引入滞留系数α(0≤α≤1)描述线性水库的蓄水能力,采用二项式[(1-α)+α]N近似模拟Nash瞬时单位线,对水文模型的汇流方法进行改进,即:
U'(t)={(1-α)N,λ2(1-α)αN-1,...,λn(1-α)nαN-n,...,αN}
各项即为二项式[(1-α)+α]N展开项,各项系数1,λ2…λn组成杨辉三角(帕斯卡三角),例如当N=1时,汇流分配系数为1-α和α;当N=2,汇流分配系数为(1-α)2、2(1-α)α、α2,当α=0时,即为原汇流方法(等流时线法);
采用二项式[(1-α)+α]N方法近似模拟Nash瞬时单位线的构建,既解决了Nash瞬时单位线难以直接应用、往往需要转换为时段单位线的问题,又避免了时段单位线和等流时线法在时空上的均化效应(时段单位线将流域上呈现为分布状态的净雨过程集总成单一过程;等流时线法假定全流域流速相同,且只考虑了不同远近的来水的滞后现象)。需要说明的时,此改进只对有上游子流域的嵌套流域出口产生作用,对于无上游流域的子流域出口,由于不考虑河道汇流,其效果仍等同于应用原方法汇流;
(5)假定坡面汇流和河道汇流的速度分别为RV和CHV,则坡面汇流时间t0为:
t0=L(1)/CHV
其中,L(i)为第i级河道到达出口断面的最远距离,则河道的汇流时间为:
ti=t0+(L(i)-L(1))/RV
其中,ti(i=1,…,m)为第i级河道的汇流时间,则流域汇流时间为tm,图1给出了m=3时的汇流过程示意图;第i级河道径流的滞后时段为TCH(i)=ti/dt,其中dt为计算时段,当ti/dt为整数时最先到达出口的径流滞后时段为j0=INT(t0/dt),否则j0=INT(t0/dt)+1;当tm/dt为整数时,最远点径流滞后时段n=INT(tm/dt),否则n=INT(tm/dt)+1;如果汇流时间多于计算时段,即tm>dt,则在第k个时段产生的径流将在k+j0~k+n时段内陆续到达出口断面;
对多个子流域组成的嵌套流域情况,对上游某一子流域p,η表示当前子流域到达流域出口的时间延迟,则演进单位线可表示为二项式[(1-α)+α]j+η,j=j0+1,…,n,在第k个时段产生的径流将在k+j0+η~k+n+η时段内陆续到达出口断面,则在k+j0~k+j(j=j0+1,…,n)时段内,到达出口断面径流的比例A(j):
其中,ACH(i)为L(i)对应的流域面积占流域总面积的百分比,为二项式[(1-α)+α]j+η的第j+η+1个展开项;
假定第k时段流域的产流量Qk,第k+j时段内,到达出口的径流比例为ΔA(j)=A(j)-A(j-1),则第k+j段时k时段的产流在流域出口形成的流量为:
Qk+j=Qk+j×ΔA(j)
分别计算各时段产流量在流域出口形成的流量过程,在不考虑流域河槽调蓄情况下,将同时出现在出口的流量直接叠加,得整个降水的模拟流量过程。
有益效果:本发明保留了瞬时单位线从时间和空间尺度上描述流域汇流过程的优势,同时避免了应用中复杂的参数问题,适用于在大尺度流域水文模拟,尤其在气候变化下流域水循环研究中具有重要的实用价值,该方法应用简便,且给出的模拟结果也更为客观科学、精度更高。
附图说明
图1为三级河道汇流示意图(m=3);
图2为实施例中卢森堡Atter流域;
图3为Nash线性水库汇流方法;
图4中(a)为线性滞留汇流法一阶汇流示意图,(b)为二阶汇流示意图;
图5为U出口断面流量模拟情况。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
实施例:线性水库滞留汇流及嵌套流域(多子流域)汇流方法,其核心是引入滞留系数α(0≤α≤1)描述线性水库的蓄水能力,采用二项式[(1-α)+α]N方法近似模拟Nash瞬时单位线,对水文模型的汇流方法进行改进,其过程一般包括研究流域资料收集与分析,线性水库滞留汇流及嵌套流域(多子流域)汇流方法的建立,参数率定,以及嵌套流域的递阶推演,流域出口断面流量模拟评判等。以卢森堡Atter流域为例,具体步骤如下:
(1)以Bissen站为流域出口,分为9个嵌套的子流域,如图2所示,收集数据情况如下:每个子流域出口的流量数据为4年,4个雨量站的降雨量数据,1个蒸发观测站的蒸发量数据,流量数据、蒸发数据和降雨数据时段长均为1小时,数据系列长度为4年,根据Bissen流域水文气象条件选取水文模型TOPMODEL,计算流域产流。
(2)根据流域的DEM地形数据生成流域水系图,划分子流域,其汇流过程为:(2—4—9—1;3—7—8—9—1;4—9—1;5—8—9—1;6—9—1;7—8—9—1;8—9—1;9—1;其中:1—B(Bissen,下同,用子流域出口名字的首字母代称该子流域出口);2—H;3—C;4—S;5—N;6—P;7—E;8—R;9—U)。
(3)按照各个子流域断面分析降雨、流量的代表性和一致性,对不符合流域整体趋势的异常点进行剔除处理,尤其加强对流量资料的一致性分析。
(4)应用Nash线性水库方法,将流域内汇流路径视为N个串联的线性水库,上一个水库的出流为下一个水库的入流,通过一系列水库的演算,得到Nash瞬时单位线,即:
式中:U(t)为流域瞬时单位线,N为线性水库的个数,K为蓄水系数,t为时间步,Γ(N)为伽玛分布函数,应用该式联立入流方程Q(t)可求得流量过程线,如图3所示,a(1),a(2),…,a(N)表示一级级水库汇流过程,图b(1),b(2),…,b(N)表示出口断面流量过程线随汇流进程变化情况,即最终得到出口断面的流量过程线。
(5)引入滞留系数α(0≤α≤1)描述线性水库的蓄水能力,采用二项式[(1-α)+α]N近似模拟Nash瞬时单位线,图4中以一阶与二阶为例,对水文模型的汇流方法进行改进,即
U'(t)={(1-α)N,λ2(1-α)αN-1,...,λn(1-α)nαN-n,...,αN}
各项即为二项式[(1-α)+α]N展开项,各项系数1,λ2…λn组成杨辉三角(帕斯卡三角)。
(6)将参数α作为模型的汇流参数与模型其它参数一起应用选出的历史数据进行率定。
(7)引入确定性系数DC评判模型的模拟能力,见表1和图5:
表1各流域出口流量模拟情况
式中:y0(i)为实测值,yc(i)为预报值,是实测值的均值,n为资料序列长度。
由表1可以看出,本发明建立的方法对有上游子流域的情况有一定改善作用,且上游子流域越多(汇流路径越长)的流域出口断面的径流模拟精度的提高越显著(如总流域出口断面B)。图5中,横坐标表示时间序列(每小时/h),纵坐标主轴表示径流过程(mm/h),纵坐标副轴表示降雨过程(mm/h),虚线表示模拟值,实线表示实测值,从模拟过程和实测过程的接近程度可以看出,改进后的汇流方法使U出口断面的流量模拟效果有了一定的提高,从下图对较大流量的放大过程来看,从峰现时间差异和与实测值的逼近程度来看,α=0.1的模拟效果要优于α=0(原汇流方法)的模拟效果。需要注意的是,并不是参数α的值越大越好,表1中α的值与模拟效果有一个先增后减的过程,具体应用中需要率定以寻求最优值。由此表明,本发明可更有效地提高大尺度流域尤其是嵌套流域径流的模拟能力。
如上,尽管参照特定的优选实施例已经表示和表述了本发明,但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下,可对其在形式上和细节上作出各种变化。
Claims (1)
1.一种线性水库滞留汇流及嵌套流域汇流方法,其特征在于:包括以下步骤:
(1)收集研究流域降雨量、蒸发量和径流量的长系列数据,分析数据的代表性和一致性;
(2)选择研究流域的水文模型,由降雨量、蒸发量和径流量应用水文模型计算产流过程;
(3)应用Nash线性水库方法,将流域内汇流路径视为N个串联的线性水库,上一个水库的出流为下一个水库的入流,得到Nash瞬时单位线;
(4)引入滞留系数α(0≤α≤1)描述线性水库的蓄水能力,采用二项式[(1-α)+α]N近似模拟Nash瞬时单位线,对水文模型的汇流方法进行改进,即:
U'(t)={(1-α)N,λ2(1-α)αN-1,...,λn(1-α)nαN-n,...,αN}
各项即为二项式[(1-α)+α]N展开项,各项系数1,λ2…λn组成杨辉三角;
(5)假定坡面汇流和河道汇流的速度分别为RV和CHV,则坡面汇流时间t0为:
t0=L(1)/CHV
其中,L(i)为第i级河道到达出口断面的最远距离,则河道的汇流时间为:
ti=t0+(L(i)-L(1))/RV
其中,ti(i=1,…,m)为第i级河道的汇流时间,则流域汇流时间为tm;第i级河道径流的滞后时段为TCH(i)=ti/dt,其中dt为计算时段,当ti/dt为整数时最先到达出口的径流滞后时段为j0=INT(t0/dt),否则j0=INT(t0/dt)+1;当tm/dt为整数时,最远点径流滞后时段n=INT(tm/dt),否则n=INT(tm/dt)+1;如果汇流时间多于计算时段,即tm>dt,则在第k个时段产生的径流将在k+j0~k+n时段内陆续到达出口断面;
对多个子流域组成的嵌套流域情况,对上游某一子流域p,η表示当前子流域到达流域出口的时间延迟,则演进单位线可表示为二项式[(1-α)+α]j+η,j=j0+1,…,n,在第k个时段产生的径流将在k+j0+η~k+n+η时段内陆续到达出口断面,则在k+j0~k+j(j=j0+1,…,n)时段内,到达出口断面径流的比例A(j):
其中,ACH(i)为L(i)对应的流域面积占流域总面积的百分比,为二项式[(1-α)+α]j+η的第j+η+1个展开项;
假定第k时段流域的产流量Qk,第k+j时段内,到达出口的径流比例为ΔA(j)=A(j)-A(j-1),则第k+j时段中的k时段的产流在流域出口形成的流量为:
Qk+j=Qk×ΔA(j)
分别计算各时段产流量在流域出口形成的流量过程,在不考虑流域河槽调蓄情况下,将同时出现在出口的流量直接叠加,以得到整个降水的模拟流量过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410536584.4A CN104281780B (zh) | 2014-10-11 | 2014-10-11 | 线性水库滞留汇流及嵌套流域(多子流域)汇流方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410536584.4A CN104281780B (zh) | 2014-10-11 | 2014-10-11 | 线性水库滞留汇流及嵌套流域(多子流域)汇流方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104281780A CN104281780A (zh) | 2015-01-14 |
CN104281780B true CN104281780B (zh) | 2016-03-23 |
Family
ID=52256647
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410536584.4A Active CN104281780B (zh) | 2014-10-11 | 2014-10-11 | 线性水库滞留汇流及嵌套流域(多子流域)汇流方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104281780B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105160178B (zh) * | 2015-09-02 | 2016-08-24 | 国家电网公司 | 一种考虑降雨空间分布特征的水库流域子流域划分方法 |
CN105976103B (zh) * | 2016-05-03 | 2019-06-07 | 大连理工大学 | 一种基于动态蓄水容量的洪水预报方法 |
CN113379110B (zh) * | 2017-11-21 | 2022-08-05 | 中国水利水电科学研究院 | 一种中长期径流预报结果趋势检验方法 |
CN108446502B (zh) * | 2018-03-22 | 2019-01-29 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
CN108763829B (zh) * | 2018-06-22 | 2020-04-14 | 中国水利水电科学研究院 | 基于水动力相似的区域降雨模拟系统降雨过程设计方法 |
CN109033605B (zh) * | 2018-07-18 | 2020-10-23 | 华中科技大学 | 一种基于多阶段划分和多单位线选择的流域汇流模拟方法 |
CN109992868B (zh) * | 2019-03-25 | 2020-08-18 | 华中科技大学 | 一种基于异参离散广义Nash汇流模型的河道洪水预报方法 |
CN110543688B (zh) * | 2019-08-05 | 2022-10-14 | 河海大学 | 水库群调蓄下多阻断二维扩散波方程的水文模型汇流方法 |
CN110955924B (zh) * | 2019-09-11 | 2020-12-29 | 中国水利水电科学研究院 | 一种考虑淤地坝影响的子流域汇流模拟方法 |
CN112200386B (zh) * | 2020-10-30 | 2024-05-28 | 华中科技大学 | 一种基于等容瞬时单位线的水文预报方法及系统 |
CN112464493B (zh) * | 2020-12-10 | 2023-03-07 | 中国电建集团昆明勘测设计研究院有限公司 | 基于topmodel模型的改进模型、区域径流和洪水风险设计的方法 |
CN112651189B (zh) * | 2020-12-11 | 2023-03-10 | 中国电建集团昆明勘测设计研究院有限公司 | 一种基于自然子流域的通用流域水循环模拟计算方法 |
CN113919685B (zh) * | 2021-10-08 | 2023-04-18 | 中国水利水电科学研究院 | 一种优化生态流量分配的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034003B (zh) * | 2010-12-16 | 2012-11-28 | 南京大学 | 基于蓄水容量曲线和topmodel的流域水文模型的设计方法 |
CN102034002A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 一种高分辨率的全分布式水文模型topx的设计方法 |
JP5701919B2 (ja) * | 2013-03-14 | 2015-04-15 | 株式会社東芝 | 配水池水位設定装置、配水池水位設定方法及び配水池水位設定システム |
-
2014
- 2014-10-11 CN CN201410536584.4A patent/CN104281780B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104281780A (zh) | 2015-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104281780B (zh) | 线性水库滞留汇流及嵌套流域(多子流域)汇流方法 | |
CN102314554B (zh) | 一种基于陆气耦合的小流域洪水预报方法及系统 | |
CN109815305A (zh) | 一种无资料地区场次洪水径流过程反演的方法 | |
CN101290713A (zh) | 一种结合周相似特性的分形交通流预测方法 | |
CN109711095B (zh) | 一种基于水文模型获取河段区间入流的方法 | |
CN105046953A (zh) | 一种短时交通流组合预测方法 | |
CN104778503A (zh) | 一种考虑实时气象因素的公共楼宇空调短时基线负荷预测方法 | |
CN105320827A (zh) | 降雨径流预报系统 | |
CN109272186B (zh) | 一种水资源调度方法 | |
CN108876021A (zh) | 一种中长期径流预报方法及系统 | |
CN102155938A (zh) | 一种反推水库入库流量过程的测算方法 | |
CN105912836A (zh) | 一种纯遥感数据驱动的流域水循环模拟方法 | |
CN104268653A (zh) | 基于集束径流预报的梯级水库优化调度方法 | |
CN104182634B (zh) | 一种梯级水电站联合运行水位控制断面优选方法 | |
CN113065090A (zh) | 一种区域农田灌溉用水量分析计算方法 | |
CN104217260B (zh) | 一种风场邻近多台风电机测量风速缺损值的组合填充系统 | |
CN103745602A (zh) | 一种基于滑窗平均的交通流量预测方法 | |
CN108763829B (zh) | 基于水动力相似的区域降雨模拟系统降雨过程设计方法 | |
CN103870699B (zh) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 | |
CN108824349B (zh) | 基于水动力过程相似的模型流域单位线的构建方法 | |
CN104732110B (zh) | 基于混合weibull分布的森林生物量多尺度估测方法 | |
CN106446388A (zh) | 基于蒙特·卡洛算法的新安江模型参数优化方法 | |
Egawa et al. | A water flow forecasting for dam using neural networks and regression models | |
CN106407160B (zh) | 含多个风电场的电力系统概率潮流联合分布的计算方法 | |
Chunhui et al. | Natural runoff changes in the Yellow River Basin |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |