博士□基地班硕士□硕博连读研究生□兽医硕士专业学位□学术型硕士☑工程硕士专业学位□农业推广硕士专业学位□全日制专业学位硕士□同等学力在职申请学位□中职教师攻读硕士学位□高校教师攻读硕士学位□风景园林硕士专业学位□西北农林科技大学研究生课程结课论文封面（课程名称：土壤水分溶质动力学）学位课☑选修课☐研究生年级、姓名 2 vccccccccccccc研究生学号 XXXXXXXXXXX所在学院（系、部） XXXXXXXXXXXXXXX学院专业学科农业工程任课教师姓名 XXXXXXXX考试日期考试成绩评卷教师签字处土壤入渗实验报告一、垂直入渗实验1、实验目的测定土壤的垂直入渗特征曲线，掌握测定方法。
了解土壤一维入渗特性，确定入渗条件下土壤累积入渗量曲线以及入渗速率数学表达式，用不同的入渗经验公式描述入渗速率并绘制相应的图表。
2、实验要求（1）土柱圆筒高约29cm ，内径10cm 。
控制装土容重为1.43g/cm 。
垂直入渗过程中，进水端的水位由马氏瓶控制。
入渗过程中，观测不同时间的累积入渗量。
（2）根据实验数据在方格纸上点绘入渗过程线（速度~入渗时间），确定饱和入渗速度k 值。
（3）根据实验数据在双对数纸上点绘入渗曲线，确定α及k 值，写出该种土壤的入渗公式。
（4）略述土壤入渗过程，入渗性强弱，分析原因。
3、实验原理（1）实验利用马氏瓶供水并维持稳定水压；（2）对于均质土的入渗强度，已有若干计算公式，菲利普根据严格的数学推导，求得解析解为：f 21i t2i +=S式中，i ——t 时刻的入渗速率；S ——与土壤初始含水率有关的特性常数，成为吸水率； f i ——稳定入渗率，即饱和土壤渗透系数。
在非饱和土壤入渗初期，S 起主要作用，所以菲利普公式可以改写为：21t2i S =考斯加可夫根据野外实测资料的分析，发现入渗强度与时间之间成指数关系，其形式为：-a 1t i i =式中，1i ——第一个单位时间的入渗强度；a ——反应土壤性质与入渗初始时土壤含水率的经验常数；累积入渗量：)1(0110101)(αααα---=-===⎰⎰t i t i dt t i dt t i I tt4、实验仪器马氏瓶、玻璃土柱、天平、滤纸、秒表、烧杯5、实验方法和步骤（1）测量玻璃土柱的内径R =10cm ，控制土壤容重为1.43g/cm ，计算2cm 高土柱所需的土量g V m 8.2192)210(4.12=⨯⨯=⋅=πγ； （2）分层装土：在玻璃柱底部放入一片滤纸，然后装土。
称土219.8g ，控制每次装土高度为2cm ，压实，装下一层时表面打毛。
（3）调节马氏瓶，使其进气孔高出土柱面2cm ，形成水头差；（4）用烧杯迅速向玻璃柱中加水，立即打开供水阀，同时打开秒表计时，三者要求同时进行，然后开始实验；（5）每隔一定时间测定一组实验数据（记录马氏瓶读数、时间）； （6）处理与分析实验数据。
6、实验数据处理（1）测定的实验数据及数据处理见表1（2）累积入渗量、入渗速率与时间的关系曲线分别见图1、图2。
图1 累积入渗量I与入渗时间t的关系曲线图2 入渗速率i与入渗时间t的关系曲线由图可以看出整体结果还是比较符合的，但是相关性不太高。
7、对三种经验公式的拟合效果进行对比分析（1）土壤垂直累计入渗量散点图及各经验公式拟合曲线①使用Philip入渗公式进行拟合如下图3：图3 philip入渗公式拟合图②使用考斯加可夫入渗公式进行拟合如下图4：图4 考斯加可夫入渗公式拟合图③使用Green-Ampt入渗公式进行拟合如下图5：图5 Green-Ampt入渗公式拟合图（2）对土壤垂直入渗三种经验公式拟合效果进行对比分别用三种经验公式拟合的回归系数见下表2表2 不同拟合公式回归系数对比表拟合效果对比图如下图6由上述图表可知：（1）累积入渗量随时间的增加而增大，但增加的速率越来越小，最终趋于稳定，入渗速率随时间的增加而减小，最后趋于稳定，其值为稳定入渗率；（2）采用考Philip入渗公式拟合效果最好，其次是采用斯加可夫公式拟合，而采用Green-Ampt入渗公式拟合的效果也比较好，但是相对于前两种模型来说误差还是比较大的。
二、水平入渗实验1、实验目的在熟练掌握水平土柱吸渗法测定非饱和土壤水扩散率原理的基础上，了解土壤水分水平入渗的过程、入渗特性，确定入渗条件下湿润锋x 和时间t 之间的关系，了解入渗条件下土壤累计入渗量曲线以及数学表达式，计算土壤的入渗速率以及数学表达式，同时得到土壤水扩散率)(θD 的关系，并绘制相应的图表。
2、实验要求水平土柱（长30cm ），是由直径5cm 、厚度为2cm 的单环组装形成的，控制装土容重为1.43g/cm 。
水平入渗过程中，进水端的水位由马氏瓶控制。
入渗过程中，观测不同时间的累积入渗量以及湿润锋的距离。
实验结束后，用烘干法分层测定土壤重量含水率计算体积含水率。
3、实验仪器马氏瓶（3cm×5cm ）、玻璃土柱、天平、滤纸、秒表、烧杯。
4、实验原理水平土柱入渗法测非饱和土壤水扩散率)(θD 要求土柱的土壤质地均一,且初始含水量分布均匀,在进水端边界含水量稳定不变,忽略重力作用的条件下,该土壤水属于一维流问题,其微分方程及定解条件为:⎪⎪⎪⎩⎪⎪⎪⎨⎧====⎥⎦⎤⎢⎣⎡∂∂∂∂=∂∂)3(0t 0x ),()2(0t 0),()1()(s 0＞，，＞θθθθθθθt x x t x x D x t 式中：t 为入渗时间（min ）；x 为水平入渗距离（cm ）；θ为体积含水量（cm³/cm³）；0θ，s θ分别为土壤的初始含水量和饱和含水量。
对方程(1)施以Boltzmann 变换可得： ⎰⋅⋅-=θθθλθλθ021)(d d d D带入Boltzmann 变换参数2/1-=xt λ，上式可变为：θλθλθ∆⋅∆∆-=∑)(21)(D式中：λ为相邻两点θ对应的λ平均值；θλ∆∆为相邻两点λ和θ增量比值θλ∆∆的平均值。
Philip 水平入渗公式：21t S I ⋅=或者2121-⋅=t S i5、实验方法和步骤（1）称土：控制土壤容重为1.43g/cm ，计算4cm 高土柱所需的土量：g V m 9.1094)25(4.12=⨯⨯=⋅=πγ（2）装土柱：在水平实验土槽底部垫上滤纸，然后将实验用土按设计容1.43g/cm 的标准分层装入水平土槽中，为保证土的均匀性，我们将土按4cm 高度分层装入。
（3）装水：在马氏瓶中装入一定量的水，将下部进气阀和出水阀关闭。
（4）连接仪器：用橡皮输水管将马氏瓶的出水口与水平土槽进水口相连，然后打开马氏瓶顶部的加水孔的橡皮塞和出水阀，同时将水平槽的排气孔打开，给水平土槽下部的水室进行排气和充水，保证水能够均匀的入渗。
（5）待水室充满水后，立即将马氏瓶的加水孔和水平土槽的排气孔密封，打开马氏瓶下部的进气阀，将水平土柱放平，让水平土柱中心轴与马氏瓶的进气阀相平，这样才能保证水平入渗在无压条件下进行，同时，打开秒表开始计时，并记下马氏瓶上的刻度数。
（6）进行连续观测，每当湿润锋前进2cm 时，记下时间和马氏瓶上的刻度数，达到稳定入渗时，停止实验，然后打开水平土槽，将其中的土按2cm 长度分层装入事先准备好的铝盒中，然后称重，并放入烘箱进行烘干、称重。
6、实验数据处理与分析（1）计算累积入渗量I，结果见表1；绘制湿润锋x与时间t的关系曲线，见图1；再绘制累积入渗量I与时间t的关系曲线，见图2。
图1 湿润锋x与时间t的关系曲线图2 累积入渗量I 与时间t 的关系曲线（2）根据5.0ST 总总=I 反解5.0/T I 总总=S ，计算入渗率0.5S/2t i = ，结果见表2，然后绘制入渗速率i 与时间的关系曲线，见图3。
图3 入渗速率i 与时间t 的关系曲线（3）计算土壤体积含水率，计算结果见表3。
（4）计算 5.0x/t 总=λ，其中总t 为总累计入渗时间，可得出不同含水率θ对应的λ值。
结果见表4，然后绘制二者关系图，见图4。
图4 含水率θ与λ的关系曲线(5)根据θ和λ的值，计算土壤扩散率)D(θ(θD，结果见表5，然后绘制θ和)的关系曲线，见图5。
图5 扩散率)(θD 与体积含水率θ的关系曲线（6）经过对实验数据的有关处理，可以得到如下发现：①通过拟合得到累计入渗量I 与时间t 之间满足函数关系：5816.0909.6x y =相关系数2R=0.9971。
②同样的经拟合得到入渗率i和时间t之间满足函数关系式：5.0-y=9515.5x相关系数2R=1。
③在得到水平土柱吸渗法测定非饱和土壤水扩散率)D及累计入渗量I和(θ入渗率i的基础上,可以确定水平入渗条件下湿润锋x和时间t之间满足函数关系：4327.0y=1931.2x相关系数R2=0.9979。
从上面有关曲线的拟合相关关系来看，决定系数R2值都在0.997以上，说明该实验数据结果比较可靠。