本发明提供一种草原河流河道流量汇流计算方法,包括:根据极值选取原理,通过每个栅格高程数据获得每个栅格内水流的流向;通过高清遥感影像以及实地测量数据提取河流特征量;通过气温、降水、实际蒸散量和栅格面积,计算获取每个栅格的径流深时间序列;通过每个栅格面积和径流深计算河道的产流量、流速和河流水深;根据能量守恒定律,利用栅格内河流实际液体元流能量方程,分别计算栅格内的河流的沿程水头损失和弯道带来的局部水头损失,获得水流流出栅格的流速和水流流经栅格的实际所需时间;通过每个栅格的流向,获取每个格点流向流域泄流口所需经过的格点层数;根据径流深时间序列和河道的流量,获取当前格点在该时刻的流量。
1.一种草原河流河道流量汇流计算方法,其特征在于,包括以下步骤: 获取草原河流流域内每个栅格内水流的流向; 获取草原河流的河流特征量; 通过草原河流的气温、降水、实际蒸散量和栅格面积,计算每个栅格单位时间的径流深时间序列; 通过每个栅格面积和径流深,计算河道的产流量、流速和河流水深; 基于栅格内河流实际液体元流能量方程,根据河流特征量、河宽和河流水深,分别计算栅格内河流的沿程水头损失和弯道带来的局部水头损失,并确定水流流出栅格的流速和水流流经栅格的时间; 根据每个栅格内水流的流向,获取每个格点流向流域泄流口所经过的格点层数;并根据径流深时间序列、水流流出栅格的流速、水流流经栅格的时间和河道的产流量,获取当前格点在各个时刻的实际地表径流量; 所述栅格内的河流的沿程水头损失和弯道带来的局部水头损失的计算包括以下步骤: 根据1km弯曲半径当量的概念,将河流弯道的长度、角度和换算至同一量级,用于统一流域内河道的弯曲程度; 确定崩溃系数,即根据实时河流水深来判定洪峰过境时,是否会出现漫流的情况; 当漫流发生时,栅格河道将重置为无弯曲且拥有基础河长的状态;漫流结束后,河道随着地转偏向力等要素的影响,逐渐开始弯曲,即河长逐渐向实际河长恢复并出现弯曲河段; 构建栅格内河流实际液体元流能量方程: 式中,z1和z2为栅格入口和出口处的高程;p1和p2为栅格入口和出口处的气压;ρ为液体密度,g为重力加速度;v1和v2为栅格入口和出口处的流速;hw为总水头损失; 根据河流特征量以及河宽和河流水深,计算栅格内的河流的沿程水头损失和弯道带来的局部水头损失: hw=∑hf+∑hj n=(n0+n1+n2+n3+n4)×m5 式中,hw为损失的能量,包括水流的沿程水头损失hf和弯道处的局部水头损失hj;λ为沿程水头损失系数;LR为实际河长;R为水力半径;Re为雷诺数;ζ为局部水头损失系数;C为谢才系数;b为弯道河宽;r为河道的弯曲半径;n为河道摩阻系数,其中n0至n4分别为天然河道的基本粗糙度、不规则水面的影响、河道横截面形状和尺寸变化的影响、阻水物的影响和植物的影响;m5为河道曲折情况的影响; 所述当前格点在各个时刻的实际地表径流量的获取包括以下步骤: 通过每个栅格的流向,获取每个格点流向流域泄流口所需经过的格点层数j,设定流域格点行列号分别为m和n,则正在处理的格点层可表示为m(j),n(j);根据径流深时间序列和河道的流量,获取当前格点在t时刻的流量为Q(t)m(j),n(j);该时刻的径流量流向下一格点的时间(Δt)为: 式中,LR为河流长度;平均排放速度;当非整数的时候,将流量按照整数时间分割,令t时刻流出该栅格的流量为q(t)m(j),n(j),则: 其中,fix是一个向下四舍五入的函数。
2.根据权利要求1所述的草原河流河道流量汇流计算方法,其特征在于,所述草原河流流域内每个栅格内水流的流向获取,包括:根据极值选取原理,通过草原河流的流域边界和流域内每个栅格高程数据获得每个栅格内水流的流向。
3.根据权利要求1所述的草原河流河道流量汇流计算方法,其特征在于,所述草原河流的河流特征量的获取,包括:通过高清遥感影像以及实地测量数据提取草原河流的河流特征量。
4.根据权利要求1所述的草原河流河道流量汇流计算方法,其特征在于,还包括:当流域内存在洼地和闭塞湖地形时,每个栅格内水流的流向计算包括以下步骤: 确定流域出口,寻找每个栅格前往出口的流向路径; 当路径检测到进入死循环后,根据循环特征判断洼地形状,并寻找洼地的泄流口,进而向主流汇集,最终确定每个栅格内水流的流向。
5.根据权利要求3所述的草原河流河道流量汇流计算方法,其特征在于,所述河流特征量包括每个栅格内部的实际河长、平均河宽、河流弯曲角度及半径、河道糙率和坡降。
6.根据权利要求5所述的草原河流河道流量汇流计算方法,其特征在于,所述产流量Qsim、流速v、河宽WR和河流水深HR的计算公式如下: RD=QsimΔt/1000AG Qsim=AS×v=WR×HR×v 式中,AS为每个栅格的横截面积;AG为栅格面积;Δt为单位时间;RD为径流深时间序列。
7.根据权利要求1所述的草原河流河道流量汇流计算方法,其特征在于,还包括:当存在上游栅格的流量汇入时,计算该栅格获取的初始流量: 其中,Qsim(t)为由流量生成模块计算出的每个网格的流量,dir=1 to 7,表示1~7个上游收敛方向。
请联系平台
请联系平台