流域产沙分布式模型地研究.pdfVIP

．66． 第六届全国泥沙基本理论研究学术讨论会论文集 流域产沙分布式模型的研究 曹文洪 祁 伟 郭庆超 (中国水利水电科学研究院泥沙研究所 北京 100044) 摘要：为反映因下垫面因子空间分布不均和人类活动对流域侵蚀产沙的影响，开发和建立了基于场次暴 雨的流域侵蚀产沙分布式数学模型。该模型能够模拟出流域在不同水土保持措施(或不同土地利用类型) 下的径流和侵蚀产沙的时空过程，从而加强了模型在水土保持措施制定中的应用，为配置流域内水土保 持措施和检测流域管理提供技术支撑。 关键词：流域分布式模型降雨侵蚀产沙 流域侵蚀与产沙是一个极其复杂的系统，很早就受到了国内外学者的重视，并相应 地开发出不同类型的产沙模型。但由于受到诸多自然因素的制约，同时又受到人类活动 的影响，各因素之间存在着错综复杂的相互作用。 目前国内的流域产沙模型多为集总式模型(LumpedModel)，即把影响过程的各种不 同的参数进行均一化处理，进而对流域水文泥沙过程的空间特性实行平均化模拟，其模 型结果不包含流域水文泥沙过程空间特性的具体信息。这类模型所采用的实际上是一种 经验方法，精确度较低，通用性不高。 分布式模型(DistributedModel)将流域细化为多个连续的单元，结构上可反映流域内水文 泥沙的时空变化过程，还能对流域内任一单元进行模拟计算，是当前流域水文模型(降雨径流 模型)的研究热点，相关研究文献也颇多，然而分布式模型在侵蚀产沙方面的研究文献还较少 见。随着我国西部大开发战略的实施，生态环境建设的力度不断加强，也迫切需要发展侵蚀 产沙分布式模型，以便对流域不同措施及其配置的减沙治理效果进行预测和评价。 有鉴于此，开发基于场次暴雨的流域侵蚀产沙分布式数学模型，可以模拟流域在不 同水土保持措施(不同土地利用类型)下的径流和侵蚀产沙的时空过程，从而加强模型在水 保措施制定中的应用及检测流域管理措施对径流输沙过程产生的影响，进而能为配置流 域内水土保持措施和检测流域管理提供技术支撑和便利。 1模型结构 1．1降雨径流子模型 1．1．1有效降雨计算 有效降雨等于天然实际降雨扣除降雨损失所剩下的部分，而降雨损失一般包括植物 基金项目：科技部公益专项资助(项目编号：2004DIB4J169)。 作者简介：曹文洪，男，1963年出生，所长，教授级高级工程师。 二、流域侵蚀及产沙 ‘67· 截留、填洼、蒸发蒸腾等损失的雨量。我国西部地区以超渗产流为主，一般情况下一次 暴雨历时不长，降雨损失主要是植物截留损失和土壤入渗，蒸发蒸腾和填洼损失雨量一 般很小，可以近似忽略。故有效降雨为： le=I—f—J (1) 式中：厶为有效降雨强度；，为实际降雨强度；f为土壤入渗率；J为植物截留强度。 植物截留过程与土壤的渗透过程很相似，即当降雨被植物林冠截留后，首先会产生 初始截留；当植物林冠截留饱和后，仍具有一定的截留能力，称稳定截留强度。它们分 别相当于Horton土壤下渗方程中的初始下渗率和稳定下渗率。因此，采用Horton方程的 形式来描述植物截留的物理过程： Jr=J。+(Jo—J。)e“‘ (2) 式中：Z为降雨过程中任意时刻的植物截留强度，mm／h；以为植物林冠稳定截留强度， mm／h；Jo为植物林冠初始截留强度，mm／h；。c为植物林冠特性系数，％，取值与降雨强 度和植物种类(林分)有关；f为降雨历时。 研究表明，初始截留强度厶与降雨强度，和郁闭度彳直接相关，且有aoAI；而稳 定截留强度上与降雨强度，和林冠特性系数仅有关，可由试验观测结果计算得到[11。 1．1．2土壤入渗计算 土壤人渗计算采用物理概念明晰的Green．Ampt入渗方程： 厂=ddFf=K[1+(B—oi)sF／F】 (3) F=舱坶c良吲tn(·+丽F] c4， 式中：F为土壤累积入渗量