土壤饱和导水率

1、引言
土壤饱和导水率是土壤重要的物理性质之一，它是计算土壤剖面中水的通量和设计灌溉、排水系统工程的一个重要土壤参数，也是水文模型中的重要参数，它的准确与否严重影响模型的精度。

下文介绍了确定饱和导水率的三类方法：按公式计算，实验室测定和田间现场测定，并对其研究现状进行分析，对同类研究有重要的参考价值。

饱和导水率由于土壤质地、容重、孔隙分布以及有机质含量等空间变量的影响空间变异强烈。

王小彬等[1]研究了容重及粒径大小对土壤持水性的影响，并对各种物料处理（或措施）的保水效果及其对土壤持水特征的影响进行了探讨。

研究结果表明，随着容重的增大，土壤的饱和导水率迅速下降；刘洪禄、杨培岭等[2]研究了波涌灌溉土壤表面密实层饱和导水率k与土壤机械组成、土壤容重、供水中断时间的定量关系。

研究结果表明，随着容重的增加，饱和导水率逐渐减小，但随着黏粒含量的增加，饱和导水率的变化率变小；吕贻忠等[3]针对鄂尔多斯沙地生物结皮进行调查，利用人工喷水模拟降雨分析结皮对土壤入渗性能的影响。

结果表明，3种土壤的饱和导水率随着土壤剖面深度的增加呈现出上土层高中间土层低、底土层又升高的趋势，扰动土与原状土的饱和导水率差异较大，达到显著水平，土壤容重、孔隙度、有机质含量、黏粒含量和全盐含量等均对土壤饱和导水率有一定的影响；Helalia认为有效孔隙率与土壤饱和导水率相关性明显。

单秀枝[4]通过测定并分析不同有机质含量的壤质土样的饱和导水率、水分特征曲线、水分扩散率及几个水分常数，研究结果表明，随着有机质含量的增加，土壤饱和导水率呈抛物线变化，当有机质含量为15 g/kg时，饱和导水率达到最大值。

汪志荣、张建丰等[5]根据不同温度条件下的入渗资料，分析了活塞（Green Ampt）公式在温度场中的适用性，认为
Green-Ampt公式适用于温度场影响下的土壤水分运动；Hopmans和Duley[6]研究了土壤温度对土壤特性的影响，结论表明，随着温度的增加，土壤饱和导水率增大。

邓西民等[7]在实验室对北京壤质黏土犁底层原状土柱进行模拟冻融处理，观测冻融对其容重、孔隙度、导水率的影响。

研究结
果表明，冻融处理后犁底层土壤饱和导水率提高1.4~7.7倍；Larson研究表明冻融会改变土壤结构、构造和其他物理形状，对土壤饱和导水率由增加的作用。

秦耀东等[8]对土壤中大孔隙流进行研究，并用一种简单的方法对土壤内的大孔隙和中小孔隙的饱和导水率进行分析，结果表明，造成土壤导水率较大空间变异的主要原因是土壤大孔隙分布的不均一性，一旦土样剔除大孔隙的影响，也就是在其质域范围内，土壤中小孔隙分布相对较为均一，因而其饱和导水率的变异性也就大大变小；陈风琴等研究了缙云山常绿阔叶林下土壤饱和导水率和土壤大孔隙的关系，结果表明，饱和导水率具有较大的空间变异性，变异系数达67%，其大小不仅取决于总孔隙度，更取决于能导水的大孔隙的数量和大小，饱和导水率的变化对大孔隙变化具有高度依赖性，且与半径大于0.1 cm的大孔隙体积有较好的相关性。

2 材料与方法
2.1 研究区概况
气候特征：邵阳市位于湖南省西南部，为半山半丘陵区，属中亚热带季风湿润气候，光照充足，水雨丰沛，四季分明，气候温和，夏少酷热，冬少严寒。

受地貌多样、高差悬殊影响，气候既有东、西部的地域差异，又有山地与丘平区的垂直差异，形成一定的小气候环境和立体气候效应。

境内年平均气温16. 1~17. 1C，无霜期272~304天，日照1347.3~1615.3小时，降水量1218.5~1473.5毫米；雨水大多集中在4~6月，易遇夏秋连旱。

地貌地势：邵阳市属江南丘陵大地形区。

地形地势的基本特点是：地形类型多样，山地、丘陵、岗地、平地、平原各类地貌兼有，以丘陵、山地为主，山地和丘陵约占全市面积的三分之二，大体是“七分山地两分田，一分水、路和庄园，东南、西南、西北三面环山，南岭山脉最西端之越城岭绵亘南境，雪峰山脉耸峙西、北，中、东部为衡邵丘陵盆地，顺势向中、东部倾斜，呈向东北敞口的筲箕形。

邵阳市为江南丘陵向云贵高原的过渡地带，西部雪峰山脉、系云贵高原的东缘，东、中部为衡邵丘陵盆地的西域。

市境北、西、南面高山环绕，中、东部丘陵起伏，平原镶嵌其中，呈由西南向东北倾斜的盆地地貌。

邵阳境内系江南丘陵向云贵高原过渡地带，
南岭山脉绵亘南境，雪峰山脉耸峙西、北，横邵丘陵盆地展布中、东部。

整个地势西南高而东北低，顺势向中、东部倾斜，呈东北向敞口的筲箕形。

最高峰为城步苗族自治县东部二宝顶，海拔2021米；最低处是邵东县崇山铺乡珍龙村测水岸边，海拔仅125米，地势比降为10.25%。

境内溪河密布，有5公里以上的大小河流595条，分属资江、沅江、湘江与西江四大水系。

资江干流两源透巡，支派纵横，自西南向东北呈“Y”字型流贯全境，流域面积遍及市辖9县3区。

巫水源出城步，横贯绥宁，西入沅江，为境内西南部的主要水道。

生物资源：邵阳位于南岭山脉、雪峰山脉与云贵高原余脉三大植物区系交会地带，是湖南四大林区之一。

1990年，全市林地面积为1588.89
万亩，活立木总蓄积量达2843万立方米，森林覆盖率为42.7%。

其中绥宁县活立木蓄积量1050万立方米，人均32. 94立方米，森林覆盖率为67.7%，均居全省各县之首，有“神奇的绿洲”之誉。

植物种类繁多，珍稀树种丰富。

邵阳市植物种类多达2826种，分属245科，792属，其中被子植物105科372属1659种，裸子植物8科33属67种，自然分布和引进栽培的木本植物115科409属1726种，用材林树种210种，以杉木、马尾松和阔叶用材林为大宗。

经济林树种432种，楠竹、油茶、油桐、漆树、板栗、乌桕、白蜡树、山苍子树等成片分布。

2.2 材料方法
本次试验共在邵阳市七里坪邵阳学院周边典型坡地(50m×360m)范围进行网格(10m×10m)取样,共布设8个采样点，每个采样点采集两个土壤样本。

采样点的基本情况如下:
一号点：经纬度27°11.494N，111°26.749E；高程为350m±4m 该坡坡向为NE60°；环刀编号1、2；土地类型：林地一号点表
层为枯枝落叶层，半分解腐殖质层，腐殖质层厚度为2cm，采样
深度为0—10cm，10cm以下为半风化母质层。

采样时间为：2012年4月1日15：00。

二号点：经纬度27°11.438N，111°26.687E；高程330±3m该坡坡向为SE150°；环刀编号3、4；该处土壤较疏松，为半阳坡，表
层为枯枝落叶层，半分解腐殖质层，腐殖质厚约1cm，土壤粗骨
性强，采样为半腐殖质层；采样时间为：2012年4月1日15：
18。

三号点：经纬度27°11.409N，111°26.652E；高程310±3m；坡向为SW210°环刀编号5、6；土壤类型为红壤，半腐殖质层厚度为
1cm，采样深度为0—8cm，采的是淋溶层植物根系较多，石头也
较多，粗骨性较强。

采样时间2012年4月1日15：40。

四号点：经纬度27°11.398N，111°26.578E；高程290 ±3m；环刀编号7、8；植被类型为菜地，采的是菜地耕作土，该点位于坡底，
三面环山，出口处走向为正南180°，距上瑞高速120m，表
层0—4cm根系较多，土质疏松，采样深度为7—13cm。

采样
时间为2012年4月1日16：00。

五号点：经纬度27°11.220N，111°26.620E；高程270±3m；坡向为正北；环刀编号9、10；植被类型为桔林，土壤中植物根系较多，
土壤较疏松，比菜地紧，表层为苔藓，没有枯枝落叶层，半分解
腐殖质层为1cm，采样深度为4—9cm，土壤粗骨性不明显，无
石头，根系较多。

采样时间2012年4月1日16：25。

六号点：经纬度27°13.219N，111°25.408E；高程230±4m；环刀编号
11、12；白田洲中部偏东，植被类型为灌丛，顶级演替植物为构
树，下层为蒿草，腐殖质层厚度为2cm，并有蜗牛，千足虫，蚂
蚁，根系较少，采样深度分别为0—10cm，10cm—20cm；第一
个样（环刀11号）4—10cm 土壤黏性较大。

第二个样（环刀
12号）15—20cm表层与下层无区别，20cm内质地均一，没有
挖到沙粒石层。

采样时间2012年4月3日15：50。

七号点：经纬度27°13.232N，111°25.409E；高程210±3m；环刀编号
13、14与上点直线距离大概20m，白田洲河漫滩草地，植被类
型为杂草，表层为腐殖质层，下层有蚯蚓，土狗，土壤湿润，30cm
以上为黏土，以下为细砂，土壤较潮湿；环刀13采样深度4—10cm
环刀14采样深度14—20cm；采样时间2012年4月3日16：
00。

八号点：经纬度327°10.713N，111°25.92E；高程250±4m；土地类型为稻田表层秸秆较多，下层土壤中根系较多，有蜗牛，蚯蚓，昆
虫，腐殖质较厚，20cm以上为次表层，以下为泥底层，采样深度
分别为0—10cm，10—20cm，20—30cm；环刀15、16采样
深度2—8cm；采样时间2012年4月5日16:00。

2.3 室内分析
对上述8个采样点，每个采样点采集两个土壤样本。

用烘干法测定表层0～30cm土壤含水量。

每个采样点用环刀采原状土(为了减少表层枯枝落叶等杂物的影响,采样在表层5cm以下) ,重复2次,带回室内用定水头法测定土壤饱和导水率。

采样区山体走向为东西方向,整个坡面除距坡顶110m处较陡外,其他坡度为20°左右,坡上位110 m内为退耕20 a的苜蓿地,坡中位是退耕5a的杏树林地,坡下位为农田和少量杏树。

因采样前农田谷物已收割,退耕杏树林地种植较为稀疏,苜蓿地退耕年限较长,苜蓿的生产力很低,因而土壤表层含水状况均不考虑植被影响。

在采样期间,所处区域20d内无降雨记录,因而土壤含水量较低。

实验样本采集完毕后我们将其带回邵阳学院城市建设系土壤实验室进行进一步的分析。

首先在选定的实验地上，用环刀采取原状土，将垫有滤纸的低筛网盖好，并将环刀浸入事先准备装有水的容器中，注意水面不要超过换刀。

然后统一饱和12个小时。

其次在预定时间（12小时）将换刀取出，置于事先准备好的支架上，与此同时准备马氏瓶4-5个，根据实验要求在马氏瓶中装一定量的纯净水，并在马氏瓶出水口套好橡皮管。

将橡皮管一头放入换刀中并通入纯净水。

待重力水滴完后在环刀下部装上漏斗，漏斗下接一烧杯，待稳渗后并记录环刀水头高度。

待漏斗下面滴下滴一滴水时开始用秒表计时，每隔5秒更换漏斗下的烧杯（间隔时间短的，视渗漏快慢而定）并分别用事先准备好的量筒计算对应时间的渗出水量，并用温度计记录水温。

土壤饱和导水率系数在单位水压梯度下，通过垂直于水流方向的单位土
壤截面积的水流速度，又称土壤渗透系数。

在饱和水分的土壤中根据达西（H.Darcy ）定律得出土壤饱和导水率公式为：
)
1(t ⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⨯⨯⨯=h S L Q K 式（1）中：
K —饱和导水率（渗透系数），cm/s ；
Q —流量，渗透过的一定截面积S （2cm ）的水量，ml ； L —饱和土层厚度，渗透经过的距离，cm ； S —渗透管的横截面积，2cm ； t —渗透过水量Q 时所需要的时间，S ； h —水层厚度，水头（水位差），cm 。

饱和导水率（渗透系数）与土壤孔隙数量、土壤质地、结构、盐分含量、含水量、和温度等有关。

饱和导水率（渗透系数）K 的量纲为cm/s 或
mm/min 或cm/h 或m/d 。

在本次试验当中其中L （饱和土层厚度，渗透经过的距离，cm )取值为
5cm ；S （渗透管的横截面积，2cm ）取值为202cm ；本次试验的土样分析与
测定均是在邵阳学院城市建设系土壤实验室完成。

3 结果与分析
3.1 时间间隔和质量差的测定
将实验样本带回邵阳学院城市建设系土壤实验室，采用渗桶法测得质量差数据如下表：
表1-1 时间差与质量差测量数据
1
2
3 4 时间间隔（min) 质量差（g) 时间间隔
（min) 质量差（g
时间间隔
（min) 质量差
（g) 时间间
隔（min)
质量差
（g) 10 2.09 5 11.21 5 19.86 5 79.79 10 1.95 5 10.98 5 20.17 5 79.39 10
2.07
5
10.76
5
20.01
5
79.68
10 1.99 5 10.9 5 19.45 5 79.72 10 2.07 5 10.45 5 19.52 5 80.27 10 2.09 5 10.53 5 19.93 5 80.59 10 1.98 5 10.28 5 19.67 5 80.85 10 1.97 5 10.27 5 19.90 5 81.2 10 1.88 5 10.22 5 19.44 5 81.81 10 1.89 5 10.33 5 19.41 5 82.25 10 1.82 5 10.34 5 19.40 5 80.88
表1-2 时间差与质量差测量数据
5 6 7 8
时间间隔（min) 质量差
（g)
时间间隔
（min)
质量差
（g)
时间间隔
（min)
质量差
（g)
时间间隔
（min)
质量差
（g)
5 47.09 5 4.33 5 80.80 5 62.77 5 47.21 5 4.25 5 74.22 5 62.25 5 46.97 5 4.33 5 75.02 5 61.03 5 47.08 5 4.17 5 78.87 5 66.22 5 47.10 5 4.27 5 77.88 5 65.79 5 46.97 5 4.32 5 77.79 5 66.45 5 46.99 5 4.27 5 78.01 5 65.70 5 45.88 5 4.33 5 78.11 5 65.48 5 47.88 5 4.34 5 78.09 5 65.52 5 46.99 5 4.3
6 5 78.08 5 65.33 5 46.98 5 4.34 5 78.0
7 5 64.55
表1-3 时间差与质量差测量数据
9 10 11 12
时间间隔（min) 质量差
（g)
时间间隔
（min)
质量差
（g)
时间间隔
（min)
质量差
（g)
时间间
隔（min)
质量差
（g)
5 187.55 5 132.1
6 5 2.9
7 5 4.15 5 185.4 5 128.5
8 5 2.98 5 4.12 5 183.5
9 5 126.36 5 3.05 5 4.02 5 182.66 5 124.16 5 2.95 5 3.79 5 181.23 5 122.19 5 2.86 5 3.74 5 180.3 5 143.57 5 3.03 5 3.71 5 180.11 5 145.55 5 3.19 5 3.5 5 180.13 5 143.44 5 3.2 5 3.5 5 180.09 5 145.88 5 3.17 5 3.48 5 181.01 5 146.45 5 3.19 5 4.15 5 180.98 5 146.13 5 3.18 5 4.12
表1-4 时间质量差与质量差测量数据
13 14 15 16
时间间隔（min) 质量差
（g)
时间间隔
（min)
质量差
（g)
时间间隔
（min)
质量差
（g)
时间间隔
（min)
质量差
（g)
10 1.57 5 20.24 5 32.51 5 47.12 10 1.55 5 20.02 5 31.34 5 45.69 10 1.53 5 19.46 5 30.59 5 44.49 10 1.53 5 18.46 5 30.45 5 43.17 10 1.55 5 18.5 5 28.4 5 42.42 10 1.55 5 18.26 5 28.18 5 41.41 10 1.33 5 17.96 5 27.94 5 40.58 10 1.35 5 18.05 5 27.55 5 39.45
10 1.34 5 17.61 5 27.48 5 38.75
10 1.4 5 17.52 5 27.14 5 37.95
10 1.39 5 17.72 5 27.02 5 37.21
注：其中1号土壤样本和13号土壤样本因其土质较黏着，里面粘土较多，因此测量时间采取10分钟每次。

其他土样统一5分钟每次。

3.2 土壤容重和饱和导水率的计算
通过采用渗桶法测得八个点共十三个样本数据后进行整合得出质量差
如表1-1到1-4所示。

质量差数据测定完毕后，将环刀取下，放入烘箱中105℃烘干8个小时，等土壤冷却至室温后，将样本过秤测得1-13号样本的土壤容重。

为保证实验的精确性，在饱和导水率的计算过程中要使用达到稳定的数据进行计算（单位时间内渗出水量相等为止），所以在表1-1到表1-4中要选择相对稳定的数据进行土壤饱和导水率的计算。

其中1号样本测得数据相对稳定，因此全部采用并带入公式（1）计算出饱和导水率；2号样本采取的数据为11.21-10.28之间的数据；3号样本采取的数据为19.86-19.67之间的数据；4号样本采取的数据为79.79-80.85之间的数据；5号样本数据基本稳定，所以全部采用；6号样本数据基本稳定，因此全部采用；7号样本数据基本稳定，因此全部采用；8号样本采取的数据为66.22-65.33之间的数据；9号样本采用的是数据为187.55-180.30之间的数据；10号样本采用的数据为132.16-122.19之间的数据；11号样本所有数据基本稳定，因此全部采用；12号样本采用的数据为4.15-3.48之间的数据；13号样本所有数据基本稳定，因此全部采用；14号样本采用的数据为18.46-17.72之间的数据；15号样本采用数据为28.40-27.02之间的数据；16号样本数据基本稳定，因此全部采用。

数据采集完毕后将稳定数据带入公式（1）中进行土壤饱和导水率的计算，算得数据如下表：
表2各样本土壤容重和土壤饱和导水率计算表
样点编号高程（m) 土壤容重(g/cm3) 土壤饱和导水率(cm/s)
一号点1
350m
1.3662 0.2542
2 1.4141 0.5961
二号点3
330m
1.2221 0.5851 4 1.2309
2.9374
三号点5
310m
1.3489 1.6793 6 1.4024 0.1642
四号点7
290m
1.0778 3.2881 8 1.0778 3.2891
五号点
9
270m
1.0837 7.0561 10 1.1195 3.9590
六号点11
230m
1.2359 0.1250 12 1.2008 0.0001
七号点13
210m
1.0462 0.0292 14 0.9783 0.7504
八号点15
250m
0.8788 0.9223 16 0.8875 0.3674
测得样本土壤饱和导水率后，将一号点到八号点中包含样本的容重和土壤饱和导水率求平均值获得一号点到八号点土壤容重和土壤饱和导水率，如下表：
表3 土壤容重和土壤饱和导水率
高程(m) 土壤容重(g/cm3) 土壤饱和导水(cm/s) 一号点 350m 1.3901 0.4252
二号点 330m 1.2265 1.7612
三号点 310m 1.3757 0.9218
四号点 290m 1.0778 3.2885
五号点 270m 1.1016 5.5076
六号点 230m 1.2184 0.0626
七号点 210m 1.0123 0.3898
八号点 250m 0.8832 0.6440
3.3 相关性分析
测得以上数据后我们对一至八号点的高程与土壤饱和导水率和土壤容重与饱和导水率进行一次相关系数的比较分析。

首先进行一至八号点高程与土壤饱和导水率相关性分析，如表4-1和图1：
表4-1 高程和土壤饱和导水率记录表
高程(m) 土壤饱和导水率(cm/s)
350 0.4252
330 1.7612
310 0.9218
290 3.2885
270 5.5076
230 0.0626
210 0.3898
250 0.6449
图1 高程和土壤饱和导水率对比分析曲线图
由表1-4和图1的对比分析中我们看出低山林地不同海拔高度表层土壤
高程与土壤饱和导水率的相关关系较为明显。

其关系式为：
57.13.3789177.351729.00005.079.6101.423456+-+-+--+--=x x x x x x y
注：关系式中x 为高程，y 为低山林地不同海拔高度表层土壤饱和导水率；
计算得出：R 2=0.9395；由此可以看出高程与低山林地不同海拔高度表层
土壤饱和导水率成明显的线性正相关关系。

同时我们也对一至八号点的土壤容重和土壤饱和导水率进行了相关系
数的分析，具体内容如表4-2和图2：
表4-2 土壤容重和饱和导水率
土壤容重(g/cm 3)
土壤饱和导水率(cm/s) 1.3902 0.4252 1.2265
1.7612 1.3757
0.9218 1.0778
3.2885 1.1016
5.5076 1.2184
0.0626 1.0123
0.3898 0.8832 0.6449
图2土壤容重和饱和导水率对比分析曲线图
由表4-2和图2的相关性对比分析我们可以看出土壤容重与低山林地表层土壤饱和导水率无明显的线性相关性。

其关系式为:
+
y2-
33
x
.
=x
-
.
328
124
39
501
74
.
计算得出;R2=0.2792;式中x为土壤容重，y为低山林地不同海拔高度表层土壤饱和导水率。

4 结论与讨论：
我们对邵阳市邵阳学院七里坪校区附近的低山林地不同还把高度表层土壤饱和导水率进行测定。

期间在邵阳市七里坪邵阳学院周边典型坡地(50 m×360m)范围进行网格(10m ×10m)取样,共布设8个采样点，每个采样点采集两个土壤样本，带回邵阳学院城市建设系土壤实验室进行饱和导水率的测定分析。

本次试验我们通过渗透筒法，测量了邵阳学院七里坪校区附近低山林地不同海拔高度表层土壤饱和导水率，结果表明低山林地不同海拔高度表层土壤饱和导水率的变化范围为0.06255-5.5075815cm/s,其中海拔高度为270m的五号点植被类型为桔林，土壤中植物根系较多，粗骨性不明显，所含石头等颗粒较少，土质较疏松，因此五号点的土壤保和导水率较高。

而高程为230m的六号点由于土壤植物根系较少，粘土多土壤粘性
较大，因此导致土壤饱和导水率偏低。

通过对高程和土壤饱和导水率，土壤容重和土壤饱和导水率之间的线性相关性分析我们了解到，高程和低山林地不同海拔高度表层土壤饱和导水率成线性正相关关系，而土壤容重与低山林地不同海拔高度表层土壤饱和导水率无明显的线性相关性。

参考文献:
[1] 王小彬,叶元林,王安明,等.容重及粒径大小对土壤持水性的影响研究[J].灌溉排水学报, 2003, 22(3): 15-18.
[2] 刘洪禄、杨培岭等.波涌灌溉土壤表面密实层饱和导水率k与土壤机械组成、土壤容重、供水中断时间的定量关系的研究[J].灌溉排水学报，2001，18(3)：19-22
[3] 单秀枝.土壤有机质与水动力学的研究[M].北京:科学出版社,2002.
[4] 吴华山,陈效民,叶民标,等.太湖地区主要水稻土的饱和导水率及其影响因素研究[J].灌溉排水学报,2006, 25(2): 46-49.
[5] Hopmans和Duley.土壤温度对土壤特性的影响研究分析第三版[M].北京:中国农业出版社,1999，(4):40-47.
[6] 邓西民.对北京壤质黏土犁底层原状土柱进行模拟冻融的处理研究[J].水利学报,2002,(2):3640-3647.
[7]秦耀东等.土壤中大孔隙流进行研究.南京农业大学学报,1994,17(4)∶134-137
[7]Kammar R S, Rizvi H A, Ahmed M,etal.Measurement of field-saturated hydraulic conductivityby usingGuelph and Velocity permeameters.Transactions of the ASAE, 1989,32(6):1885-1890
[8] 华孟、王坚.土壤物理学.北京农业大学出版社,1993,284~285,82~84
[9]Rrynoids W D et al.Institu measurement of field-saturated hydraulic conductivity sorptivity and the aypa-ramete using the Guelph.permeamenter. Soil Sci, 1985, 140(4)∶292-302
[10]樊艺峰，土壤饱和导水率研究现状分析，2001，18(3):11-13
[11] 许一飞,雷廷武.降水入渗过程的模型研究[J].北京农业工程大学学报,1993,13(4):39-46.
[12] 许迪,蔡林根.冬小麦—夏玉米种植模式下的农田水量平衡模拟及入渗补给规律分析[J].水利学报,1997,12: .
[13] 周择福,洪玲霞.不同林地土壤水分入渗和入渗模拟的研究[J].林业科学,1997,33(1):9-17.
[14]郝振纯.黄土地区降雨入渗模型初探[J].水科学进展,1994,5(3):186-193.
[15] 郭维东.坐水播种时耕层土壤水分入渗的二维数值模拟[J].农业工程学报,2001,17(2):24-27.
[16] 何园球,等.红壤丘岗区人工林土壤水分、养分流失动态研究[J].水土保持学报,2002,16(4):91-93.
[17] 李鸿杰,杜历.土壤水、盐、入渗变异特性及其相互关系的空间序列分析[J].土壤学报,1993,30(1):60-68.
[18] 刘贤赵,康绍忠,等.黄土区坡地降雨人渗产流中的滞后机制及其模型研究[J].农业工程学报,1999,15(4):95-99.
[19] 田瑾,项静恬,等.经济、环境等非线性系统的预测和调控[M].北京:中国统计出版社,2001，30(1):55-58
毕业设计(论文)附录
目录
1. 原始数据 (17)
2. 文献综述 (18)
3. 外文原文 (22)
4. 外文译文 (28)
附录1:原始数据
饱和导谁率原始记录数据
1 2 3 4
水头(cm) 0.2
水头
(cm)
0.9
水头
(cm)
1.7
水头
(cm)
1.3
温度(℃)20
温度
(℃)
20
温度
(℃)
122.21
温度
(℃)
20
土重(g) 136.62 土重(g) 141.41 土重(g) 122.21 土重(g) 123.09 时间重量(g) 时间重量(g) 时间重量(g) 时间重(g) 0＇00 96.78 0＇00 102.64 0＇00 102.83 0＇00 155.87 10＇00 98.87 5＇00 113.85 5＇00 123.32 5＇00 230.39 20＇00 100.82 10＇00 124.83 10＇00 143.18 10＇00 304.38 30＇00 102.89 15＇00 135.59 15＇00 163.35 15＇00 377.58 40＇00 104.88 20＇00 146.49 20＇00 183.36 20＇00 450.36 50＇00 106.95 25＇00 156.94 25＇00 202.81 25＇00 523.17 60＇00 109.04 30＇00 167.47 30＇00 222.33 30＇00 595.42 70＇00 111.02 35＇00 177.75 35＇00 242.26 35＇00 667.42 80＇00 111.01 40＇00 177.74 40＇00 261.93 40＇00 738.89 90＇00 110.12 45＇00 177.77 45＇00 281.83 45＇00 96.83 100＇00 110.88 50＇00 177.56 50＇00 301.27 50＇00 809.78
5 6 7 8
水头(cm) 1.4
水头
(cm)
1.3
水头
(cm)
0.7
水头
(cm)
1.0
温度(℃)20
温度
(℃)
19
温度
(℃)
19
温度
(℃)
20
土重(g) 134.89
土重
(g)
140.24
土重
(g)
107.78
土重
(g)
103.0
时间重量(g) 时间重量(g) 时间重量(g) 时间重量(g) 0＇00 102.59 0＇00 96.76 0＇00 155.94 0＇00 101.64 5＇00 149.29 5＇00 100.96 5＇00 224.90 5＇00 101.85 10＇00 196.38 10＇00 105.29 10＇00 305.70 10＇00 102.83 15＇00 243.59 15＇00 109.54 15＇00 379.92 15＇00 103.59 20＇00 290.56 20＇00 113.87 20＇00 454.94 20＇00 102.49 25＇00 337.64 25＇00 118.04 25＇00 533.81 25＇00 103.94 30＇00 384.74 30＇00 122.31 30＇00 533.01 30＇00 103.47 35＇00 431.71 35＇00 126.63 35＇00 532.77 35＇00 102.75 40＇00 432.01 40＇00 130.90 40＇00 532.12 40＇00 103.74 45＇00 431.88 45＇00 130.33 45＇00 532.09 45＇00 103.77 50＇00 431.78 50＇00 130.43 50＇00 532.98 50＇00 102.56
附录2:文献综述
低山林地不同海拔高度表层土壤饱和导水率文献综述
1 引言：
土壤饱和导水率是土壤重要的物理性质之一，它是计算土壤剖面中水的通量和设计灌溉、排水系统工程的一个重要土壤参数，也是水文模型中的重要参数，它的准确与否严重影响模型的精度。

下文介绍了确定饱和导水率的三类方法：按公式计算，实验室测定和田间现场测定，并对其研究现状进行分析，对同类研究有重要的参考价值。

饱和导水率由于土壤质地、容重、孔隙分布以及有机质含量等空间变量的影响空间变异强烈。

土壤饱和导水率ks是单位土水势梯度下的土壤水通量，在水势用水势头表示时，饱和导水率的单位与通量单位相同，都是cm/s，或m/d。

饱和导水率综合反映了多孔介质对某种流体在其中流动阻碍作用，因此ks值一方面取决于孔隙介质的基质特征，同时，也和流体的某些物理性质，如粘度和密度有关。

曾有人将ks分成两个因子，土壤内透水率k和反映流体流动性的参数f，这样有：ks=kf（1）这种将ks分解为两个因子的方法在理论上是可能的，也有一定实意义，但由于我们主要兴趣只是研究水这种流体在土壤介质中的运动，而且在绝大多数情况下不去专门考虑温度对水的物理性质的影响，所以认为水的物理性质接近于恒量，因此不将ks分解为两个因素来研究。

2土壤饱和导水率的确定方法
确定饱和导水率的方法大致可分为以下三类：按公式计算、实验室测定和田间现场测定。

2.1公式计算确定饱和导水率的公式都是经验性的，因为影响饱和导水率的因素很复杂，许多试图发现普遍可应用的函数关系迄今为止得到的结果都令人失望，这些公式只能在极有限的条件下应用而无普遍意义，因为我们对饱和导水率的研究不应集中在这方面，这里不作详细介绍。

2.2实验室测定实验室测定饱和导水率的仪器类似Darcy的实验仪器，测得水头损失ΔH和流量Q后，如果实验中的土壤服从线性定律，则可求得到饱和导水率ks。

k=-Qω·荦H（2）其中：ω为土壤横断面面积，荦H=ΔHΔZ为水力势梯度。

实际测定应在几个不同的水力势梯度下进行，为此实验设备需略加改
进，使水头损失在一定范围内变化，具体做法是改进口土层的水头。

用这种方法测定，由于测定过程中进口土层的水头都是不变的，因此这种方法称为定水头法。

另一种实验室较常用的测定方法称为变水头法，这种方法测定颗粒较细的土壤的饱和导水率较为合适。

无论是定水头还是变水头，如何使测定土样有足够的代表性是应用这些方法测定必须考虑的一个问题。

2.3田间现场测定。

双环法双环法是使用直径不同的两个圆环（钢板或硬塑料）；降雨模拟法目前，降雨模拟器的种类很多，大都为手提式的，因此田间应用相当方便。

在做试验前，应先平整试验地块，然后插入铁制框架（1 m2左右）以分割开各重复处理区域。

准备好后，装上降雨模拟器，以特定降雨强度均匀喷洒。

各处理地块分别在一定时间间隔内测量径流量，在达到稳定径流量后，即可停止试验。

从降雨强度中减去径流量即为土壤入渗率，然后按照双环入渗法分析计算，即可求得饱和导水率。

此方法虽然达到稳定流所需时间较短，但若加上前期准备工作，则完成每次试验所需时间和双环入渗法大体相当。

以上两种方法都是按照一维垂直入渗解进行计算的，因此为了获得更为精确的饱和导水率值，必须努力使试验条件满足假定条件，避免水分侧向运动。

圭尔夫渗透仪法。

圭尔夫渗透仪是根据三维稳定入渗原理研制而成，主要由供水与量测系统、支架和入渗部件组成。

供水部分按照马氏瓶原理，由同心双管组成。

内管为进气管，外管为供水管，用以测定维持某一水位时所需的稳定流流量；入渗部件指与土壤施测面直接联系的部分，水通过该部分形成入渗。

测孔深达欲测土层，即可测定该层饱和导水率。

圭尔夫入渗仪法圭尔夫入渗仪法用于测定土壤表层饱和导水率，是圭尔夫渗透仪的一种特例，即测孔中水深为零的情况（实际测定时，一般需有2~3 mm深水位，以保持充分湿润）。

此种方法的试验数据处理及参数求解与圭尔夫渗透仪法相同。

圭尔夫渗透仪法及入渗仪法具有操作简单、试验所需时间较短、用水量少等优点。

以上介绍的测定方法只能测定个别点的土壤饱和导水率，对大面积农田土层的平均饱和导水率的测定，一般要在现场挖测试坑或打测试井，然后用注水或抽水的方法测定其流量、水势等数值，再根据达西定律求得平均饱和导水率，这种方法在水文地质、水利工程部门用得较多。

3土壤饱和导水率研究现状。