土壤水分特征曲线的测定及经验模型对比

　　【摘要】土壤水的基质势或土壤水吸力是土壤含水率的函数，它们之间的关系曲线称为土壤水分特征曲线。该曲线反映了土壤水的能量与数量关系，是反映土壤水分运动基本特征的曲线〔1〕。它是表示土壤基本水力特征的重要指标，对研究土壤水滞留与运移有十分重要的作用。
　　【关键词】土壤水分特征曲线压力膜仪经验模型参数拟合
　　1研究意义
　　土壤水分运动是陆地水循环的重要组成部分，是地表水与地下水相互作用的纽带。是降雨产流计算、农田灌溉与排水设计、地下水补给计算、土壤植物水分定量关系预测的基础〔2〕。
　　土壤水分运动3个参数中以预测非饱和导水率最为困难，土壤水分特征曲线则最容易得到，准确性也最好，方法较多，且通过水分特征曲线模型可以推求其他2个参数，因此，水分特征曲线的获取对预测土壤水分运动参数至关重要。
　　2水分特征曲线测试方法
　　（1）直接方法。分实验室法和田间方法两种方式。实验室内测定主要有张力计法、砂性漏斗法、压力膜法、离心机法和热电偶温度计测定等。田间原位测定大都用张力计法。
　　（2）经验公式法。经验公式法中比较常用的有：BrooksCorey（1964）模型，vanGenuechten（1980）模型、GardnerRusso（1988）模型等。
　　（3）间接推求法。可以分为3类：土壤转换函数方法、物理经验方法、分形几何方法。土壤转换函数就是利用已有的土壤基本性质（如粒径分布、容重、有机质含量等）通过某种算法构建起来的预测吸力与水分含量之间关系的函数〔3〕。
　　3水分特征曲线的影响因素
　　（1）土壤质地和结构：相同的含水量下，质地越细，水吸力就愈大，曲线愈陡；反之质地越粗，吸力就越小，曲线愈平缓。（2）温度：在同一吸力条件下，温度升高，土壤持水量减少，温度低时，其持水能力增强；或者，在同一含水量条件下，温度高时，吸力较低，而温度降低时，则吸力升高。（3）滞后现象：土壤水分特征曲线的滞后作用对任何质地的土壤均存在，吸水和脱水过程，负压与含水率曲线是不同的。滞后影响的程度不同，土质越轻，滞后的影响越大。目前对滞后现象的解释存在三种理论，即瓶颈理论〔1〕、接触角理论〔4〕和弯月面延迟形成理论。
　　4水分特征曲线模型
　　本文借助于RETC软件，选择具有代表性的BC模型、VG模型，对土壤水分特征曲线中的参数进行拟合，比较拟合效果。
　　BrooksCorey模型形式简单，便于推求描述土壤水分运动模型和确定土壤水分运动参数，对具有较窄孔径分布的均质和各向同性的粗质地样本效果较理想，而对于细质地土壤和未扰动的原状土通常精度较差〔56〕。
　　vanGenuchten提出的经验公式不仅能够表征整个压力水头范围内的水分特征数据，还可方便地利用统计孔径分布模型来估计水力传导率，因此在土壤水研究中比较流行〔7〕。适用土壤质地范围比较宽，同时可以使饱和土壤吸力为0，符合吸湿过程中土壤吸力变化特点，模型对粗质地的土壤拟合效果最好。
　　5研究材料
　　实验样品取自某水文地质环境地质试验场。样品编号分别为1、2、4、7、12、14号（见表一）。
　　6实验过程
　　（1）试验开始前首先检查压力室的密封性能，如无漏气现象，则压力室可以使用。（2）对陶土板、土壤样品进行饱和。排除陶土板微细孔中的空气。（3）将饱和好的陶土板擦去表面多余的水，放入压力室内。并测量饱和土样的重量，记录数据。顶盖安装好并上紧螺丝，确保压力室不漏气。（4）给压力室通入气体以增加压力，样品室的压力由压力表控制，样品在压力的作用下开始释水，并通过压力室出水孔排出室外，当达到平衡时，土的基质吸力就等于施加的气压力。（5）当48小时不再出水时即可认为达到平衡。在达到平衡后，称量每级吸力下集水瓶的重量，以便测定其含水量的变化。依次施加的吸力分别为0。02、0。05、0。10、0。20、0。32、0。54、0。84、1。20、1。66、2。24、2。94（单位：Bar）。（6）在施加最高一级基质吸力达到稳定后取出土样、烘干称重，测定应于最高吸力下的含水量。绘制基质吸力与含水量关系曲线即水分特征曲线。
　　7结果分析
　　使用RETC软件拟合数据时，选取BrooksCorey模型、vanGenuchten模型拟合参数，并进行对比，得到拟合结果（见表2）。
　　从图1、图2、图3中可看出，1、2号土样的两种模型拟合曲线与实际曲线接近，上下波动不大。12号粉土的VG、BC模型曲线与实际想接近，但与1、2号拟合结果相比，其误差最大。
　　三种土样的同种模型拟合结果，其精确度大小不同，由R2值对比可看出：VG模型拟合效果粘土（0。99857）好于好于粉土（0。99804）亚粘土（0。99779）；BC模型拟合效果粘土（0。99799）好于粉土（0。99429）好于亚粘土（0。98995）。
　　3种土壤的BC模型拟合水分特征曲线结果可以看出：在负压大于1。66Bar（基质势1500cm）时，BC模型拟合结果相对于实测值出现偏移，实测点不能完全落在拟合曲线上。即对于同一种土壤，BC模型在不同的基质势区间拟合的效果都是不一样的，VG模型在负压增大时，曲线偏移并不明显。
　　总体上比较三种模型，VG模型拟合效果好于BC模型的拟合效果。
　　实验条件环境和操作等方面会对实验结果的准确性产生影响。同时考虑测定温度对土壤水分特征曲线的影响对实验结果加以校正。
　　参考文献：
　　〔1〕张蔚榛。土壤水动力学〔M〕。北京：中国水利水电出版社，1996。9。
　　〔2〕王全九，邵明安，郑纪勇。土壤中水分运动与溶质迁移〔M〕。北京：中国水利水电出版社，2007。4。
　　〔3〕vanGenuchtenMTh，LeijFJ，LundLJ（eds。）。Proc。Int。WorkshoponIndirectMethodsforEstimatingtheHydraulicProeprtiesofUnsaturatedoSils〔M〕，UniversityofClaifonria，Rivesride，CA，1992。
　　〔4〕邵明安，王全九，黄明斌。土壤物理学〔M〕。高等教育出版社，2006。11。
　　〔5〕LeijF，RussellWB，LeschSM，ClosedformexpressionsforwaterretentionandconductivitydataLJJ，GroundWater，1997，35（5）：848858。
　　〔6〕vanGenuchtenMTh，LeijFJ，YatesSR，TheRETCcodeforquantifyingthehydraulicfunctionsofunsaturatedsoils〔R〕。USEPARepotr600291065。U。S。EnvironmenatlProtectionAgency，Ada，Oklahoma，1991。
　　〔7〕刘建立，徐绍辉，刘慧。估计土壤水分特征曲线的间接方法研究进展〔J〕。水利学报，2004，2，3：26875。