The invention provides a method and device for predicting runoff and sediment yield in a multi spatial scale watershed. The method includes: dividing the river basin to be predicted into m third-level ditches according to DEM topographic data, remote sensing images and river network of the river basin to be predicted; calculating the confluence flow of each second-level ditch by using Naxi instantaneous unit line according to the section outlet discharge of each third-level ditch; and then calculating the continuous flow of Muskingando reach. Based on the sediment yield of each three-level gully, the sediment yield of each three-level gully is determined according to the sediment yield of each three-level gully landform unit, the sediment yield of the gully landform unit and the sediment yield of the gully landform unit. Then, the sediment yield of each three-level gully is determined by the sediment transport ratio. The model coupling calculation can simulate and predict the multi-spatial scale watershed, thereby improving the accuracy and popularization practicability of the simulation and prediction of soil erosion and sediment yield.
【技术实现步骤摘要】
一种多空间尺度流域产流产沙预测方法及装置
本专利技术涉及分布式产流产沙预测术领域,尤其涉及一种多空间尺度流域产流产沙预测方法及装置。
技术介绍
当前流域洪水泥沙预报模型主要采用单一种类模型,结合降雨观测数据与下垫面数据,进行产流产沙模拟与预测预报,然而由于水文现象在不同时空尺度上呈现出性质迥异的变化特征,不同尺度上水文变量和参数的空间变异性很大,而且现有不同空间尺度上水沙循环的机理、流域水沙循环全要素过程的动态耦合模拟的关键问题尚未解决,因此,单一种类模型在进行流域水沙过程模拟过程时,存在流域范围受限、难以大范围推广的问题。
技术实现思路
为解决现有技术中存在的上述问题,本专利技术提供一种多空间尺度流域产流产沙预测方法及装置。一方面,本专利技术提供一种多空间尺度流域产流产沙预测方法,该方法包括:步骤1、根据待预测流域的DEM地形数据,采用D8算法计算每个DEM栅格流向、DEM栅格坡度和DEM栅格汇流累积量,依据Strahler河网分级方法,确定待预测流域的河网;步骤2、根据待预测流域的DEM地形数据、遥感影像和河网,将待预测流域划分为m个三级沟;步骤3、根据DEM地形 ...
【技术保护点】
1.一种多空间尺度流域产流产沙预测方法,其特征在于,包括:步骤1、根据待预测流域的DEM地形数据,采用D8算法计算每个DEM栅格流向、DEM栅格坡度和DEM栅格汇流累积量,依据Strahler河网分级方法,确定待预测流域的河网;步骤2、根据待预测流域的DEM地形数据、遥感影像和河网,将待预测流域划分为m个三级沟;步骤3、根据DEM地形数据的分辨率,对每个三级沟进行空间离散,将每个三级沟划分为多个网格单元;步骤4、根据每个网格单元中的水流从网格中心点至三级沟出口的网格汇流时间,分别采用等流时线法和线性水库法计算每个三级沟网格单元的地表汇流和地下汇流,确定每个三级沟的断面出口产 ...
【技术特征摘要】
1.一种多空间尺度流域产流产沙预测方法,其特征在于,包括:步骤1、根据待预测流域的DEM地形数据,采用D8算法计算每个DEM栅格流向、DEM栅格坡度和DEM栅格汇流累积量,依据Strahler河网分级方法,确定待预测流域的河网;步骤2、根据待预测流域的DEM地形数据、遥感影像和河网,将待预测流域划分为m个三级沟;步骤3、根据DEM地形数据的分辨率,对每个三级沟进行空间离散,将每个三级沟划分为多个网格单元;步骤4、根据每个网格单元中的水流从网格中心点至三级沟出口的网格汇流时间,分别采用等流时线法和线性水库法计算每个三级沟网格单元的地表汇流和地下汇流,确定每个三级沟的断面出口产流量;步骤5、根据每个三级沟的断面出口产流量,采用纳希瞬时单位线计算每个二级沟的沟网汇流量;步骤6、根据每个二级沟的沟网汇流量,采用马斯京根多河段连续流量演算模型计算待预测流域的断面出口产流量;步骤7、根据每个三级沟的格网坡度和遥感影像,将每个三级沟划分为梁峁坡地貌单元、沟谷坡地貌单元和沟槽地貌单元;步骤8、根据每个三级沟中梁峁坡地貌单元的产沙量、沟谷坡地貌单元的产沙量和沟槽地貌单元的产沙量,确定每个三级沟的产沙量;步骤9、根据每个三级沟的产沙量,采用泥沙输移比确定待预测流域的断面出口产沙量。2.根据权利要求1所述的方法,其特征在于,所述步骤4中根据式(1)计算网格汇流时间:其中,tc表示网格汇流时间;Kv表示速度常数;L表示网格中心点至三级沟出口的水流路径长;S表示从网格中心点至三级沟出口之间河道的坡度。3.根据权利要求1所述的方法,其特征在于,所述步骤4中分别根据式(2)、式(3)和式(4)计算每个三级沟网格单元的地表汇流、地下汇流和每个三级沟的断面出口产流量:Qj,3(t)=Qu,i+Qg,i(4)其中,Qu,i表示第i时段的地表汇流量,Qg,i表示第i时段的地下汇流量,Qj,3(t)表示第j个三级沟的断面出口量,Qg,i-1表示为第i-1时段地下汇流量,hu,i表示第i时段地面净雨量,fn为每两条等流时线之间的面积,hg,i-1表示第i-1时段地下净雨量,F表示待预测流域的总面积,△t表示时段长度,Kg表示地下水库的蓄泄系数,j=1,2…….,m,t为时段序号,1,2…,30。4.根据权利要求3所述的方法,其特征在于,所述步骤5中根据式(5)计算每个二级沟的沟网汇流量:其中,Qp,2(t)为第p个二级沟的沟网汇流量,t为时段序号,1,2…..,30;n为线性水库个数;k为线性水库的蓄泄系数;D为微分算子d/dt。5.根据权利要求4所述的方法,其特征在于,所述步骤6中根据式(6)和式(7)计算待预测流域的断面出口产流量:Ql=C0Ql-1,2+C1Ql-2,2+C2Ql-1(6)其中:C0+C1+C2=1,Q(t)表示待预测流域的断面出口产流量,Ql为第l个河段的出流量,Ql-1...
【专利技术属性】
技术研发人员:王玲玲,左仲国,姚文艺,秦奋,韩志刚,孔祥兵,杨吉山,王昌高,杨二,王志慧,李勉,张攀,倪用鑫,鲍宏喆,郭凯,侯欣欣,
申请(专利权)人:黄河水利委员会黄河水利科学研究院,河南大学,
类型:发明
国别省市:河南,41
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。