坡面降雨入渗产流规律的数值模拟研究论文.docVIP

　　坡面降雨入渗产流规律的数值模拟研究论文 摘要：本文采用运动波理论和两次改进后的Green-Ampt入渗模型建立了坡面降雨入渗产流的动力学模型，并得到了实验资料的良好验证。运用该模型研究了简单坡面上降雨入渗产流的动力学规律，分析了雨强、土壤初始含水量、渗透系数、坡面阻力，以及坡长、坡度等因素对坡面产流过程的影响规律，得出了一些有益的结论。 关键词：入渗 产流 坡面 动力学 1 概述 雨水降落在坡面上将产生雨水的聚集并形成坡面水流。坡面水流是土壤水蚀过程的主要动力，搞清产流的动力学特点是进一步研究侵蚀过程规律的基础。坡面水流不同于一般明渠流动，其水深极浅(一般只有几毫米)，沿程不断有质量源和动量源加入.freelpt积水入渗模型开始，相继有Horton(1940),Philip(1957)等模型出现，但G-A模型仍以其简单的形式，明晰的物理概念，良好的扩展性和可信的应用效果受到广泛重视，特别是经过Mein Larson(1973)和Chu(1978)的两次改进，使其可应用于不均匀降雨的入渗计算，更使它成为最有效和应用最广泛的模型。在国内，G-A模型尚未受到重视，Horton模型曾得到相当广泛的运用，但其参数的物理意义明显不如G-A模型明晰。也有研究者使用更基本的土壤水分运动微分方程，但所需的参数更加难于获取，计算也更为复杂。 本文工作旨在建立物理概念明晰的降雨入渗产流综合计算模式，并用以研究简单坡面的产流过程，分析各主要因素的影响和各主要因素的影响和各主要参量的变化规律。以期对坡面产流的动力学规律有清楚的认识。 2 计算模式 坡面流运动十分复杂，目前主要采用运动波理论、扩散波或完整圣维南方程进行描述。正如前文所述，运动波近似理论在大多数情况下可以很好地描述坡面流运动过程，且计算简单。因此本文仍采用一维运动波理论，即坡面流基本方程为 (1) 此处第二式直接使用了水力学中熟知的Chezy公式和Manning公式。其中，x为沿坡面向下的坐标，t为时间(s)，h为水深(m),q为单宽流量(m2/s)，p为降雨强度(m/s)，此处假设降雨方向垂直向下，i为入渗率(m/s)，S0为坡面坡度，S0=sinθ，θ为坡面倾角，n为Manning糙率系数。 土壤的入渗过程对坡面流的形成和流动过程影响很大，本文采用形式简单、物理概念明晰的G-A入渗模型，其计算方程为 i=dI/dt=K1+(θS-θi)S/I I=Kt+S(θS-θi)ln(1+I/S(θS-θi)(2) 其中 K为土壤饱和导水率(渗透系数)(m/s)，θS为土壤饱和含水率，即有效孔隙率(%)，θi为土壤初始含水率(%)，S为土壤吸力(m),I为累积入渗量(m)。 经典的Green-Ampt模型是干土积水入渗模型，其前提是在整个入渗过程中地表始终有积水。Mein Larson 1973年将其推广应用至降雨入渗的情况。设有稳定的雨强p，只有p大于土壤的入渗能力时，地表才能形成积水。而在降雨的初始阶段，全部降雨都渗入地下。由G-A模型知，入渗率是随累积入渗量的增加而减小的。设想当累积入渗量达到某一值时，i=p,此时开始积水，称此累积入渗量为Ip。因此由G-A模型入渗公式可以导出开始积水时的Ip值 Ip=(θS-θi)S/(p/K)-1 (3) 开始积水时间由tp=Ip/p给出。因此整个过程的入渗率可表示为 i=p t ≤tp i=K1+(θS-θi)S/I t tp (4)式中的I为积水开始后的累积入渗量(包含未积水时段的入渗量在内)。由于不是由t=0开始积水，I的计算须采用修正后的公式 Kt-(tp-tS)=I-S(θS-θi)ln1+I/S(θS-θi) t tp (5) tS表示假设由t=0开始积水，到入渗量I=Ip(或i=p)时所需时间，可理解为一个虚拟时间，可计算如下 KtS=Ip-S(θS-θi)ln1+Ip/S(θS-θi) (6) 改进的主要思想是将整个过程假设为从一开始就是积水入渗，这样该曲线在积水后部分相对于实际入渗曲线将向左平移tp-tS，将这条曲线向右平移tp-tS，再加上积水前的入渗强度等于降雨强度的关系，就得到真实的入渗过程。 但稳定降雨在实际应用中远不能满足需要，Chu(1978)将Mein Larson(1973)改进的G-A模型再作推广，使其可应用于变化的降雨过程。基本作法是，对每个计算时段将地表状态分为四种情况： 1.开始无积水，结束无积水 2.开始无积水，结束有积水 3.开始有积水，结束有积水 4.开始有积水，结束无积水 在每一时段开始，已知降雨总量与入渗总量，剩余总量。根据两个因子判断时段结束时是否有积水。 若时段开始无积水，使用因子cu,若时段开始有积水，使用因子cp，其表达式为 cu=P(tn)-R(