一种基于三级汇流结构的分布式侵蚀产沙模拟方法技术

本发明专利技术公开了一种基于三级汇流结构的分布式侵蚀产沙模拟方法，涉及土壤侵蚀、河道输沙模型与计算机数值计算技术领域。通过改进WEP‑L模型，将WEP‑L原有的“坡面‑河道”二级系统改进为“坡面‑沟道‑河道”三级系统，具体的，根据从坡面到流域出口的产输沙过程，按照土壤侵蚀发生的位置将水沙过程划分为坡面、沟道、河道三个环节，其中坡面侵蚀环节分为雨滴溅蚀、坡面径流侵蚀和坡面沟道产输沙过程；以各环节的侵蚀特点为依据，根据产输沙时空格局、承接传递过程等方面的相似性，分别以子流域和坡面等高带为计算单元，实现侵蚀产沙过程的多环节细化模拟，有效的反映了沟道侵蚀严重的特点，真实刻画了其侵蚀过程。

A Distributed Erosion and Sediment Yield Simulation Method Based on Three-stage Confluence Structure

【技术实现步骤摘要】
一种基于三级汇流结构的分布式侵蚀产沙模拟方法
本专利技术涉及土壤侵蚀、河道输沙模型与计算机数值计算
，尤其涉及一种基于三级汇流结构的分布式侵蚀产沙模拟方法。
技术介绍
在水土保持研究过程中，除了产流过程外，泥沙过程的同步模拟必不可少。然则，国内外水文模型在泥沙方面都采用半经验式的模拟方法，物理过程刻画薄弱，不能满足区域水土管理越来越精细化和精准化的需求。此外，已有的流域土壤侵蚀模型都只对坡面和河道两个侵蚀过程进行模拟，输沙过程计算时往往忽视了沟道过程的重要性。坡面侵蚀显然不能反映沟道的侵蚀特点，而沟道的比降一般都比河道比降大，河道产输沙也不能包含沟道产输沙过程。尤其，黄土高原沟道高度发育，沟道水流侵蚀量占总水力侵蚀量的比例接近50％。这一特点使沟道水力侵蚀过程在黄土高原土壤侵蚀中占据举足轻重的地位，也对水沙物理过程的模拟提出了更高的要求。目前，国内主流水文模型之一的大尺度流域水与能量转化过程模型(WaterandEnergytransferProcessesinLargeriverbasins，简称WEP-L)已很好地应用于黄土高原降水-径流研究，但该模型也仅考虑了“坡面-河道”两级汇流系统，并且不支持对流域侵蚀产沙等物理过程的模拟分析，模型不能直接应用于黄土高原水土保持分析与生态治理研究。
技术实现思路
本专利技术的目的在于提供一种基于三级汇流结构的分布式侵蚀产沙模拟方法，从而解决现有技术中存在的前述问题。为了实现上述目的，本专利技术采用的技术方案如下：一种基于三级汇流结构的分布式侵蚀产沙模拟方法，包括如下步骤：S1，地形数据提取：利用GIS平台和模拟区DEM数据，提取河网和坡度等地形数据；S2，子流域-等高带单元划分：利用GIS平台和S1中的河网提取结果，基于汇水区划分子流域；在每个子流域内，按照一定的高程间隔划分等高带；确定坡面-沟道-河道三级汇流结构；S3，基础数据库构建：收集模拟区内的水文、气象、土壤、土地利用、水土保持措施、人类活动取用水的模型输入数据；并按照模型需求和S2中的划分结果，对收集的模型输入数据进行空间或时间尺度上的处理，构建基础数据库；S4，计算坡面产输沙：利用S3所建数据库中各等高带的水文气象要素、土地利用条件、水保措施数据，以等高带为单元，按照上下承接的关系，依次计算各等高带的坡面沟间及沟道侵蚀量和输沙率；S5，计算沟道产输沙：基于S4中坡面产沙结果，以子流域为计算单元，进行沟道侵蚀、输沙模拟，计算沟道输沙率；S6，计算河道产输沙：基于S4中坡面、S5中沟道产输沙结果，以子流域为单元，进行河道侵蚀、输沙模拟；S7，侵蚀产沙量统计：坡面沟间及坡面沟道产沙由等高带汇总到子流域；坡面、沟道、河道各环节的产沙由子流域汇总到模拟区；最终，在模拟区内对所有环节的产沙进行汇总。优选地，S4包括如下步骤：S401，以等高带为计算单元，暴雨期以每小时降水为驱动，非暴雨期以日降水为驱动，计算坡面沟间雨滴溅蚀；S402，以产流计算所得坡面沟间径流为驱动，计算坡面沟间径流侵蚀；S403，每个等高带坡面沟间所产沙量进入当前等高带的坡面沟道，以坡面沟道内水流挟沙力为上限，判断坡面沟道冲淤情况，当上一等高带来沙和当前等高带来沙这两个沙源总的含沙量小于水流挟沙力时，则当前等高带坡面沟道发生冲刷；否则发生淤积；如此从上到下串联进行，实现产沙在坡面等高带间的传递。优选地，S401中，按照如下公式计算雨滴溅蚀：式中，SplMod为雨滴溅蚀侵蚀模数(kg/m2/s)；E为降雨侵蚀动能(J/m2)；I为雨强(mm/min)；α为坡面的坡度(°)；a1、b1、c1、a2、b2为经验参数。优选地，S402中，按照如下公式计算坡面径流侵蚀：SloMod＝a3·(q-b3)·(1+c3·α)式中，SloMod为坡面沟间径流侵蚀模数(kg/m2/s)；q为坡面沟间单宽流量(m2/s)；α为坡面的坡度(°)；a3、b3、c3为经验参数。优选地，S403中，按照如下公式计算水流挟沙力：式中，Tv为体积比挟沙力(m3/m3)；Cv为体积比含沙量(m3/m3)；v为水流流速(m/s)；ρs为土粒密度(kg/m3)；ρm为浑水密度(kg/m3)；g为重力加速度(N/kg)；R为水力半径(m)；ω为浑水中的泥沙沉速(m/s)；h为水深(m)；d50为泥沙的中值粒径(m)，粒径小于该值的泥沙占全部泥沙的50％；κ为浑水的卡门常数；κ0为清水的卡门常数；a4、b4、c4、d4、a5、b5为经验参数。优选地，S5具体为：以子流域为计算单元，沟道承接来自当前子流域坡面最低等高带的输沙，根据汇流结果进行沟道输沙计算。优选地，S6具体为：在子流域内，多条沟道中的输沙通过并联方式汇入河道，通过河道输沙模拟到达子流域出口的沙量；不同子流域出口的产输沙又基于上下游关系，通过串并联混合方式最终传递到模拟区出口。优选地，所述输沙计算方法中，对经典悬移质连续方程进行改进，添加侧向来沙项(qls)：式中，x为沿河方向到沟/河道入口的距离(m)；Cx为对应位置的水流含沙量(kg/m3)；Tx为对应位置的水流挟沙力(kg/m3)；α为恢复饱和系数,在一般水力因素条件下，平衡时恢复饱和系数在0.02～1.78之间，平均接近0.5；q为单宽流量(m2/s)，qls为单位长度沟道/河道侧向来沙输沙率(kg/s/m)。优选地，河段出口的含沙量按照如下公式进行计算：式中，C和C0分别为沟道或河道出口和入口的水流含沙量(kg/m3)；T和T0分别为沟道或河道出口和入口的水流挟沙力(kg/m3)；Cl为沟道/河道单位长度侧向来沙的含沙量(kg/m3)；L为河段长度(m)。本专利技术的有益效果是：本专利技术提供的基于三级汇流结构的分布式侵蚀产沙模拟方法，通过改进WEP-L模型，将WEP-L原有的“坡面-河道”二级系统改进为“坡面-沟道-河道”三级系统，并基于侵蚀产沙和输移机制增加流域侵蚀与产沙模块，具体的，根据从坡面到流域出口的产输沙过程，按照土壤侵蚀发生的位置将水沙过程划分为坡面、沟道、河道三个环节，其中坡面侵蚀环节分为雨滴溅蚀、坡面径流侵蚀和坡面沟道产输沙过程；以各环节的侵蚀特点为依据，根据产输沙时空格局、承接传递过程等方面的相似性，分别以子流域和坡面等高带为计算单元，实现侵蚀产沙过程的多环节细化模拟，有效的反映了沟道侵蚀严重的特点，真实刻画了其侵蚀过程。附图说明图1是本专利技术提供的基于三级汇流结构的分布式侵蚀产沙模拟方法流程示意图；图2是具体实施例中涉及到的模拟区示意图；图3是窟野河流域温家川水文站1968-1992年输沙率模拟结果示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，本专利技术提供了一种基于三级汇流结构的分布式侵蚀产沙模拟方法，包括如下步骤：S1，地形数据提取：利用GIS平台和模拟区DEM数据，提取河网和坡度等地形数据；S2，子流域-等高带单元划分：利用GIS平台和S1中的河网提取结果，基于汇水区划分子流域；在每个子流域内，按照一定的高程间隔划分等高带；确定坡面-沟道-河道三级汇流结构；S3，基础数据库构建：收集模拟区内的水文本文档来自技高网...

【技术保护点】
1.一种基于三级汇流结构的分布式侵蚀产沙模拟方法，其特征在于，包括如下步骤：S1，地形数据提取：利用GIS平台和模拟区DEM数据，提取河网和坡度等地形数据；S2，子流域‑等高带单元划分：利用GIS平台和S1中的河网提取结果，基于汇水区划分子流域；在每个子流域内，按照一定的高程间隔划分等高带；确定坡面‑沟道‑河道三级汇流结构；S3，基础数据库构建：收集模拟区内的水文、气象、土壤、土地利用、水土保持措施、人类活动取用水的模型输入数据；并按照模型需求和S2中的划分结果，对收集的模型输入数据进行空间或时间尺度上的处理，构建基础数据库；S4，计算坡面产输沙：利用S3所建数据库中各等高带的水文气象要素、土地利用条件、水保措施数据，以等高带为单元，按照上下承接的关系，依次计算各等高带的坡面沟间及沟道侵蚀量和输沙率；S5，计算沟道产输沙：基于S4中坡面产沙结果，以子流域为计算单元，进行沟道侵蚀、输沙模拟，计算沟道输沙率；S6，计算河道产输沙：基于S4中坡面、S5中沟道产输沙结果，以子流域为单元，进行河道侵蚀、输沙模拟；S7，侵蚀产沙量统计：坡面沟间及坡面沟道产沙由等高带汇总到子流域；坡面、沟道、河道各环节的产沙由子流域汇总到模拟区；最终，在模拟区内对所有环节的产沙进行汇总。...

【技术特征摘要】
1.一种基于三级汇流结构的分布式侵蚀产沙模拟方法，其特征在于，包括如下步骤：S1，地形数据提取：利用GIS平台和模拟区DEM数据，提取河网和坡度等地形数据；S2，子流域-等高带单元划分：利用GIS平台和S1中的河网提取结果，基于汇水区划分子流域；在每个子流域内，按照一定的高程间隔划分等高带；确定坡面-沟道-河道三级汇流结构；S3，基础数据库构建：收集模拟区内的水文、气象、土壤、土地利用、水土保持措施、人类活动取用水的模型输入数据；并按照模型需求和S2中的划分结果，对收集的模型输入数据进行空间或时间尺度上的处理，构建基础数据库；S4，计算坡面产输沙：利用S3所建数据库中各等高带的水文气象要素、土地利用条件、水保措施数据，以等高带为单元，按照上下承接的关系，依次计算各等高带的坡面沟间及沟道侵蚀量和输沙率；S5，计算沟道产输沙：基于S4中坡面产沙结果，以子流域为计算单元，进行沟道侵蚀、输沙模拟，计算沟道输沙率；S6，计算河道产输沙：基于S4中坡面、S5中沟道产输沙结果，以子流域为单元，进行河道侵蚀、输沙模拟；S7，侵蚀产沙量统计：坡面沟间及坡面沟道产沙由等高带汇总到子流域；坡面、沟道、河道各环节的产沙由子流域汇总到模拟区；最终，在模拟区内对所有环节的产沙进行汇总。2.根据权利要求1所述的基于三级汇流结构的分布式侵蚀产沙模拟方法，其特征在于，S4包括如下步骤：S401，以等高带为计算单元，暴雨期以每小时降水为驱动，非暴雨期以日降水为驱动，计算坡面沟间雨滴溅蚀；S402，以产流计算所得坡面沟间径流为驱动，计算坡面沟间径流侵蚀；S403，每个等高带坡面沟间所产沙量进入当前等高带的坡面沟道，以坡面沟道内水流挟沙力为上限，判断坡面沟道冲淤情况，当上一等高带来沙和当前等高带来沙这两个沙源总的含沙量小于水流挟沙力时，则当前等高带坡面沟道发生冲刷；否则发生淤积；如此从上到下串联进行，实现产沙在坡面等高带间的传递。3.根据权利要求2所述的基于三级汇流结构的分布式侵蚀产沙模拟方法，其特征在于，S401中，按照如下公式计算雨滴溅蚀：式中，SplMod为雨滴溅蚀侵蚀模数(kg/m2/s)；E为降雨侵蚀动能(J/m2)；I为雨强(mm/min)；α为坡面的坡度(°)；a1、b1、c1、a2、b2为经验参数。4.根据权利要求2所述的基于三级汇流结构的分布式侵蚀产沙模拟方法，其特征在于，S40...

【专利技术属性】
技术研发人员：周祖昊，蔡静雅，刘佳嘉，徐东坡，龚家国，严子奇，王浩，贾仰文，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京,11