含砾石土壤降雨入渗过程模拟 Simulation of rainfall infiltration in stony soil.pdfVIP

第2l卷第6期 水科学进展 V01．21，No．6 2010年11月 ADVANCESINWATERSCIENCE Nov．，2010 含砾石土壤降雨入渗过程模拟 朱元骏1”，邵明安1’2 (1．西北农林科技大学水土保持研究所，陕西杨凌712100；2．中国科学院水利部水土保持研究所，陕西杨凌712100) 摘要：在室内模拟降雨试验的基础上，利用多重相互作用连续体(MINC)组件，对含砾石土壤降雨入渗过程进行了 数值模拟，以期深入理解含砾石土壤中的水土过程。研究结果表明，土壤中砾石的存在能够导致土壤水分饱和度增 加，使得土壤总人渗量以及降雨入渗深度降低。含砾石土壤中存在一个陡的湿润锋面，在湿润锋面上土壤水分饱和 度急剧下降。MINC组件可以很好地模拟含砾石土壤降雨入渗过程，其对含砾石土壤湿润锋和土壤含水率变化的模 拟值与实测值吻合较好。降雨期间土壤水分饱和度沿土壤剖面的变化可用来解释含砾石土壤降雨入渗深度降低的 原因。 关键词：含砾石土壤；降雨入渗；多重相互作用连续体；数值模拟；土壤水分饱和度 中图分类号：S152．7 文献标志码：A 文章编号：1001-6791(2010)06-0779-09 降雨入渗是陆地表层水分进入土壤的起点，它不仅决定进入土壤中水分的量…，还强烈影响植物根系 吸水旧1、土壤溶质迁移以及陆地表面蒸腾蒸发过程¨…，是地表水文循环中的一个重要环节。对降雨入渗过 程进行数值模拟有助于提高陆地水文循环模型精度，还可以定量评估不同土壤类型以及土地利用方式对人渗 的影响”o。基于均值介质的达西定律及其一系列改进方程可以很好地描述土壤中水分通量与土壤水势梯度 的关系，是土壤水分运动的基本方程旧o。在假定入渗时土壤中存在一个明显的干湿交界面条件下，Philip和 水分运动过程¨…。这些方法和模型的提出，加深了人们对土壤入渗过程的理解，为基于过程的水文模型提 供了理论基础和数值计算方法。然而，相关研究较少考虑土壤中出现的镶嵌体，如砾石、钙质层风化形成的 结核等。这类土壤的分布具有一定的地带性，反映了自然成土和沉积过程中土壤质地和类型的变化，例如在 黄土高原北部地区成土过程中出现的钙质结核¨…，以及河流沉积地貌区出现的夹砾石土壤。这类土壤既不 同于均质土壤，与层状土壤也有一定差异，降雨人渗过程也因此而不同。 土壤中镶嵌体对降雨入渗的影响是多方面而且复杂的。国内外研究表明，其对入渗的影响与镶嵌体类 型、尺寸、含量和位置等有关…4“。在特定的试验条件下，镶嵌体可能促进入渗；然而大部分条件下镶嵌 体抑制了土壤水分向深层入渗。由于对土壤中镶嵌体的分布很难进行描述，对入渗过程的研究大多数是建立 在模拟试验的基础上，不仅研究结果变异很大，而且存在尺度问题¨6。18J。加之对此类土壤中水势、水分的 测定存在一定难度，导致相关研究结果在实际中很难应用。本研究在室内模拟试验的基础上，结合PetraSim 软件中的多重相互作用连续体(MultipleInteracting 过程中湿润锋以及含水率进行了数值模拟和验证，以期深入理解该类型土壤中发生的水土过程，为相关水文 模型参数的选择提供可靠的依据。 收稿日期：2010-03．15 (10502．Z10-1) 作者简介：朱元骏(1978一)，男，湖北竹山人，博士，助理研究员，主要从事土壤水分运动及坡面水土过程研究。 E·mail：zhuyuanjun@foxmail．corn 万方数据 780 水科学进展 第2l卷 1材料与方法 1．1 试验材料 模拟试验在黄土高原土壤侵蚀与旱地农业国家重点实验室降雨大厅进行。试验用土填入一铁制柱中，规 格1 1 cm m(长)×0．5m(宽)xm(深)。柱体一侧中央镶嵌20cm宽带刻度有机玻璃，有机玻璃中央每隔10 开直径3 收集口。试验土壤取自杨凌农田表层0—20cm土壤，过 2him筛，自然风干，土壤质地为粉质粘土(砂粒 1