基于数字流域模型的多沙粗沙区侵蚀产沙计算.pdfVIP

中国科学E辑：技术科学 2007年第37卷第3期：446～454 《中国科学杂志社 http：／／www．scichina．tom S(：Ir=NcEIN 基于数字流域模型的多沙粗沙区侵蚀 产沙计算 刘家宏术王光谦李铁键 薛 海 贺 莉 (清华大学水沙科学与水利水电工程国家重点实验室，北京100084) 摘要 黄河难以治理的症结是水少沙多，水沙不平衡．根据多年治黄实践，黄 河中游多沙粗沙区的治理是黄土高原水土流失治理的重点，模拟计算多沙粗沙区 的侵蚀产沙量意义重大．利用数字流域模型框架建立了黄土高原多沙粗沙区产流 产沙数学模型，模型在岔巴沟小流域进行参数率定．应用率定后的模型计算了黄 产沙量，得到了多沙粗沙区的汛期径流深分布图和产沙模数分布图．计算结果能 够较好地解释黄河流域的“少水多沙”现象，为多沙粗沙区的水土保持工作提供了 科学依据． 关键词 数字流域模型 多沙粗沙区产沙计算 黄土高原多沙粗沙区主要分布在黄河流域的河口镇至龙门区间及泾河、北洛河上游地区 的黄土丘陵沟壑区(如图1)，总面积7．86x104km2【1】．该区域地形破碎，沟壑纵横，水土流失严 重，土壤侵蚀产沙机理复杂，是黄河泥沙及其粗泥沙的主要来源区．模拟计算这一区域的侵蚀 产沙量及其空间分布特征，对于研究多沙粗沙区的产沙规律、进行淤地坝规划以及评价水保 和SWAT模型等，这3个模型都是基于缓坡条件提出来的，对于黄土丘陵沟壑区这样的陡坡条 件并不适用．国内早期的流域产沙模型多为集总式模型【2】，即通过对历史资料的分析研究建立 产沙量和降雨量、降雨强度和植被因子等的经验关系式[3”61；后来一些研究者将人工神经网络 技术应用于流域产沙预报，建立了一些模型[7,81．集总式模型和人工神经网络模型均不能反映 流域产沙的空间分布特性；同时基于小流域野外观测数据或坡面实验数据的集总式经验产沙 模型在向大流域或整个多沙粗沙区推广时，普遍存在尺度效应问题．本文在文献【9]建立的数 elevation model)数据，建立了有一定物理基 字流域模型框架下，基于高分辨率的DEM(digital 础的黄土高原多沙粗沙区产流产沙分布式数学模型，实现了大范围侵蚀产沙过程的精细模拟． 收稿日期：2005．07．18；接受日期：2006．09—06 国家自然科学创新研究群体基金资助项目(批准号 4E—mail：Liujiahong@tsinghua．edu．ca 第3期 刘家宏等：基于数字流域模型的多沙粗沙区侵蚀产沙计算447 模型以“坡面”为基本单元，其面积尽量逼近真实的天然坡面面积，从而消除了模型参数的尺度 效应．模型模拟了多沙粗沙区汛期的侵蚀产沙过程，不仅得到了各年汛期的侵蚀总量，而且较 好地再现了该区侵蚀产沙的空间分布特征． 河 图1 黄土高原多沙粗沙区的范围及主要支流分布 1多沙粗沙区侵蚀产沙数学模型 黄河数字流域模型系统具有模型分层分级的优势，可以构建大区域的整体模型，将数字 流域模型框架和侵蚀产沙模型相结合就可以研究黄土高原多沙粗沙区的产沙情况．本文对黄 土高原多沙粗沙区进行两级划分：第一级分为85个小流域单元，每个单元小流域的面积平均 约为1000 图2多沙粗沙区划分及河网图 中国科学E辑技术科学 第37卷 0．37 km2．然后在每个“坡面”对应的“坡面”上进行产沙计算，最后统计得到了黄土高原多沙粗 沙区的产沙总量．以每个坡面的平均产沙模数为样本进行插值，就得到了整个多沙粗沙区的 产沙模数分布图． 在两级划分的基础上，先进行坡面产流计算，然后用坡面产沙公式进行坡面产