常熟地区水稻土饱和导水率间接方法的研究.doc

常熟地区水稻土饱和导水率的间接方法研究基金 基金项目： 研究工作得到所知识创新工程领域前沿项目（ISSASIP0719），支撑计划课题（2006BAD10A06-03），973计划课题（2005CB121103），863专题课题（2006AA10Z208）的资助。 作者简介：张均华(1983-)，博士研究生，主要从事区域土壤水分溶质运移模拟研究。Email:zhangjh@issas.ac.cn *通讯作者：刘建立，Email:jlliu@issas.ac.cn 张均华1,2 刘建立1* 朱安宁1 （1. 中国科学院南京土壤研究所，南京210008；2. 中国科学院研究生院，北京100049） 摘要：作为直接试验的一种替代方法，利用土壤基本物理性质通过土壤转换函数预测饱和导水率简便易行，成本低廉，并且预测精度能满足实际研究的需要。本研究利用目前得到较多应用的9种基于多元回归分析建立的转换函数来构建、校正预测土壤饱和导水率的经验公式，并与人工神经网络方法相比较。结果表明，人工神经网络总体预测效果要优于基于多元回归分析建立的转换函数，并且Cosby(1984)在输入参数较少的基础上预测饱和导水率精度最高。本文以Cosby(1984)预测常熟水稻土壤饱和导水率，进一步利用GIS的空间描述能力与函数的定量分析能力，得到区域尺度饱和导水率的分布状况，为该地区区域尺度数值模拟的运行提供基础参数支持。 关键词：土壤饱和导水率；多元回归函数；人工神经网络；GIS 中图分类号：S152.7 文献标识码：A 饱和水力传导率（Ks）是重要的土壤水力学参数，是定量分析土壤水分入渗、地表产流等过程的基础。土壤特别是表层土壤饱和导水率空间变异性非常强烈，需要采集大量样本或进行田间现场测定才能得到一个较大空间尺度上的饱和导水率空间分布[1]。但是通常认为室内填装样本的测定结果无法表征田间原始状态。田间原位测定土壤饱和传导率的方法主要有Guelph渗透仪、流速渗透仪、圆盘渗透仪和双管法。由于时间，人力成本等问题，原位测定仅能用于有限的点，无法满足较大空间尺度（如区域上）的研究。 已有的研究表明，土壤粒径组成，容重，有机质含量等基本土壤性质与饱和导水率之间存在一定关系，可以通过常规土壤调查中积累的大量基本土壤性质资料间接估算或预测区域尺度饱和导水率的空间分布[2-5]。并且在许多实际问题中，间接方法的预测结果也能提供足够精度的Ks预测值。但到目前为止，在长江中下游典型水稻土上的应用很少有报道。本文即由此出发，以江苏常熟市为研究区域，利用目前得到较多应用的基于多元回归分析的土壤转换函数(Pedo-transfer functions，简称PTFs)来构建、校正预测土壤饱和导水率的经验公式，并与人工神经网络方法(ANN)相比较。进一步利用GIS的空间描述能力与函数的定量分析能力，对区域尺度的饱和导水率做出定量描述与评价，为区域尺度数值模拟的运行提供基础参数支持。 1材料与方法 1.1数据来源 土壤转换函数方法普遍存在明显的地域局限性。为克服转换函数的区域适用性问题，建立PTFs时除采用45个常熟水稻土样本外，还利用地跨河南封丘、长垣、延津、原阳四县的天然文岩渠流域潮土样本169个，以及随机选取UNSODA数据库[6]中各种质地类型的土壤样本289个，共计503个样本数据。水稻土和潮土饱和导水率由A．L．Franklin Engineers生产的圆盘渗透仪(CSIRO disc permeameter)现场测定。土壤基本性质包括粒径组成、容重、有机质含量等，常规方法室内测定。土壤特征水分含量数据(如饱和含水量和-33 kPa含水量)根据土壤转换函数预测[7]。 1.2分析方法 1.2.1土壤转换函数 选用上世纪八十年代以来得到较多应用的9种PTFs，从503个样本数据中随机选取336个构建PTFs。通过与实测值比较，重新优化函数系数。剩余167个样本检验函数的预测能力，以期得到一个饱和导水率预测值与实测值最为接近的PTFs。在对函数系数进行优化处理时采用了自行编制的基于改进的Levenberg-Marquard算法的FORTRAN程序。ANN方法通过非线性迭代可以获得输入基本性质与输出饱和导水率之间（理论上）的最优关系。本研究应用无内部反馈的误差反向传播前向神经网络（BP网络）对随机选取的336个样本数据进行训练，剩余167个样本数据进行验证。进而对ANN与PTFs进行比较。研究表明，Ks多服从正态分布，因此本文在对Ks实测值进行对数转换的基础上进行Ks的预测、比较。 表1 用于本研究的预测饱和导水率的土壤转换函数 Table 1 The Pedo-transfer functions in this study 函数来源 Function S