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

常熟地区水稻土饱和导水率的间接方法研究基金项目： 研究工作得到所知识创新工程领域前沿项目（ISSASIP0719），支撑计划课题（2006BAD10A06-03），973计划课题（2005CB121103），863专题课题（2006AA10Z208）的资助。作者简介：张均华(1983-)，博士研究生，主要从事区域土壤水分溶质运移模拟研究。Email:*通讯作者：刘建立，Email:张均华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个样本数据。水稻土和潮土饱和导水率由ALFranklin 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 Source自变量个数Independent Variable输入变量Input variableCosby(1984)82Sand, Clay,Puckett(1985)9-101ClaySaxton(1986)9-103Sand, Clay, qSDanePuckett(1992)111ClayRawls(1992)123Sand, Clay, PorosityCampbellShiozawa(1994)9132Sand, ClayWosten(1997)144Silt, Clay, BD, OMWosten(1999)154BD, OM, Sand, ClayKC方程(1999)164Sand, Clay, qS, q33注：Sand土壤砂粒含量（%）；Silt土壤粉粒含量（%）； Clay土壤黏粒含量（%）；BD土壤容重(g/cm3)；OM土壤有机质（%）；Porosity土壤孔隙度（%）； q33 土壤水势为-33kPa时土壤水分含量（cm3/cm3）；qS 土壤饱和水分含量（cm3/cm3）。1.2.2 转换函数预测精度评价本文根据均方根误差、确定性系数以及Akaike信息准则评价土壤转换函数的预测精度。(1)均方根误差(Root mean square error, RMSE)，计算Ks预测值的平均精度： （1）(2)确定性系数(Coefficient of Determination,R2),反映预测值和实测值相关性的重要指标： （2）（3）Akaike信息准则（Akaikes information criterion, AIC）11,评价方法的实用性。一种方法中需要涉及的参数越多，其实用性就会下降，最好的方法应该有最小的AIC。 （3）公式(1)至（3）中，表示实测值，表示预测值，N表示样本数，np为涉及的参数数量。1.2.3饱和导水率的区域尺度扩展在对比ANN与PTFs基础上，利用预测Ks效果最好的方法预测常熟水稻土耕层土壤204个采样点的Ks。结合GIS的空间描述能力与函数的定量分析能力，进行点数据的区域尺度扩展，以期为区域尺度数值模拟的运行提供基础参数支持。图1为常熟市204个采样点的空间分布1 张福锁,朱兆良. 农田生态系统氮素行为与氮肥高效利用。图1 常熟市204个采样点空间分布Fig.1 The spatial distribution of Changshu in 204 samples2结果与分析2.1基于多元回归分析建立的土壤转换函数预测饱和导水率精度通过9种PTFS预测Ks对比实测值分析发现，各函数预测精度总体不高。AIC值平均2596.141，RMSE平均2.62（cm.d-1），R2平均0.49. 可以看出几乎所有函数模型均用到粒径分布数据，这些土壤物理性质比较容易获得，所以为预测Ks提供可能。但也有不少函数预测值与实测值差距较大，需要重新进行数据拟合，优化函数系数。2.2基于多元回归分析的土壤转换函数优化函数的可靠性通过对不同于构建函数时的独立序列的预测能力来评价17。本研究在兼顾各种质地属性的基础上，以随机选取的336个样本构建PTFs，剩余167个样本检验其可靠性。这样通过数据拟合得到的优化系数能更好的反映现有土壤饱和导水率。KC方程在公式中没有系数，不用进行优化设置。表2 调整系数后多元回归分析预测167个样本点饱和导水率与神经网络方法比较Table 2 Comparison of Predictive value between adjusted PTFs and ANN in 167 samples饱和导水率saturated hydraulic conductivity多元回归分析优化系数后比较Comparison of Ks by adjusted PTFs人工神经网络方法比较Comparison of Ks by ANN AICRMSE(cm.d-1)R2AICRMSE(cm.d-1)R2Cosby(1984)717.7740.6560.651592.7360.5870.794Puckett(1985)727.3660.6790.577618.8120.6620.708Saxton(1986)717.9010.6820.590603.0850.6240.736DanePuckett(1992)727.3660.6990.498618.8120.6620.708Rawls(1992)697.6720.6140.600587.2240.5950.727CampbellShiozawa(1994)717.7740.6770.641592.7360.5870.794Wosten(1997)702.1800.6150.495581.4370.5780.764Wosten(1999)1072.6591.8750.241610.2100.5760.746KC方程(1999)1195.9631.7290.598603.6430.6250.744表2为优化系数后PTFs生成的Ks预测值与ANN方法生成的Ks对比实测值分析。可以看出，PTFs调整系数后AIC比原来有很大程度的降低，线性相关性也有很大提高，Ks预测值与实测值之间差异性显著下降。Cosby(1984), Rawls(1992) 、CampbellShiozawa(1994)、Wosten(1997)这四个函数预测Ks效果较好。Wosten(1997)输入参数比较多，理论上参数越多，预测效果越理想。但随着输入参数的增多，采用土壤转换函数的主要优势（即实用性）也会减弱。所以本研究选择Cosby(1984)作为预测Ks的最佳函数。图1-a为Cosby(1984)饱和导水率实测值与预测值之间的对比分析。图2 Cosby(1984)多元回归函数与人工神经网络对饱和导水率的预测值与实测值比较Fig. 2 Comparison of Ks between ANN and Multiple regression analysis of Cosby(1984)通过与优化后PTFs预测能力比较，可以发现ANN方法总体上精度更高。其中以Cosby(1984)，Rawls(1992)，Wosten(1997)参数为输入变量的ANN预测效果比较好。Cosby(1984)输入参数少，在保证预测精度下能有较高的实用性，所以本研究选择Cosby(1984)的输入参数作为ANN方法的输入变量，即根据砂粒含量与黏粒含量预测Ks。图1-b为ANN方法根据砂粒含量与黏粒含量建立的Ks预测值对比实测值分析。2.3土壤基本性质对饱和导水率影响性分析通过相关性分析，可以识别出影响Ks的关键因素。表3为土壤基本性质与Ks相关性分析。可以看出Ks多与砂粒含量、孔隙度成正相关，而与粉粒、黏粒含量，-33Kpa时水分含量成负相关，这是因为大孔隙占饱和水运动的大多数。土壤质地和结构与导水率也有直接关系，砂质土壤通常比细质土壤具有更高的饱和导水率，同样，具有稳定团粒结构的土壤，比不稳定团粒结构的土壤，传导水分要快得多，后者在潮湿时结构被破坏，黏粒和粉砂粒阻塞较大孔隙的连接通道。表3转换函数所用参数与饱和导水率相关性分析Table 3 Correlation Analysis between parameters and Ks土壤基本性质soil physical properties土壤饱和导水率saturated hydraulic conductivitySand0.668*Silt-0.645*Clay-0.467*BD-0.183*Porosity0.428*OM-0.047q33-0.423*qS-0.091注：*表示置信度达95%，*表示置信度达到99%由表3还可以看出Ks与土壤粒径分布数据相关性较好，与容重，孔隙度等数据相关程度次之。所以本研究用Cosby(1984)预测Ks在实际应用中是可行的。图3为以Cosby(1984)的输入参数作为ANN输入变量预测常熟水稻土15个采样点Ks与实测值比较。RMSE为0.192，R2为0.937,可见以砂粒含量与黏粒含量预测Ks在较大尺度研究中能满足精度需要。图3 神经网络方法-Cosby(1984)预测常熟15个样本点饱和导水率Fig. 3 Predicted value of Ks in Changshu by ANN-Cosby(1984) in 15 samples2.4饱和导水率的区域尺度扩展图4 常熟市饱和导水率空间分布Fig. 4 Distribution of Saturated Hydraulic Conductivity in Changshu通过上述分析，可以发现ANN方法预测Ks效果总体上要优于PTFs。其中根据砂粒含量与黏粒含量预测的Ks与实测值最为接近。所以本研究以苏州常熟市范围内耕层土壤204个采样点的土壤粒径分布数据为基础，以Cosby(1984)的输入参数作为ANN输入变量，预测采样点Ks。通过对采样点Ks进行地统计学插值，将点源数据扩展到区域尺度。图4为常熟地区耕层Ks的空间分布状况。从图中可以看出粒径分布的空间状况，西部砂粒含量高于东部，黏粒含量与之相反，而Ks的空间分布规律与砂粒相似，这与表3相关分析中所揭示的结果相符合。3 结论土壤转换函数具有易获取输入数据，节省调查采样成本，计算简单等特点。本研究通过构建、校正PTFs，并与ANN方法比较，探求最大区域适用性的土壤转换函数。进一步的利用GIS的空间描述能力与函数的定量分析能力，对区域尤其是水稻土的Ks做出定量描述与评价，这对于该地区区域尺度数值模拟运行中基础参数的获得有非常重要的实际意义。通过分析发现PTFs与ANN中预测Ks效果最好的均为Cosby（1984）。相关性分析中，Ks与土壤粒径分布数据相关程度最高，所以本研究选择Cosby（1984）预测不同质地属性的Ks能满足精度需要。Cosby（1984）输入参数较易获得，这就为转换函数的实际应用提供了可能。而ANN对Ks的预测效果总体优于PTFS，这是由于ANN通过非线性迭代可以获得输入基本性质与输出Ks之间的最优关系。 参考文献1 Klute A, Dirken C. Hydraulic conductivity and diffusivity: Laboratory methods. In: Klute A(ed.) Methods of soil analysis. Part 1. Physical and mineralogical methods. 2nd edition. Soil Science Society of America, Madison, Wis. 1986, 687734.2黄元仿, 李韵珠. 土壤水力性质的估算土壤转换函数J. 土壤学报, 2002, 39(4): 517-523.3 van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soilsJ. Soil Sci. Soc. Am. J. 44, 892898.4朱安宁, 张佳宝, 陈效民等. 封丘地区土壤传递函数研究J. 土壤学报, 2003, 40(1): 53-58.5刘建立, 徐绍辉, 刘 慧. 估计土壤水分特征曲线的间接方法研究进展J. 水利学报. 2004, 2: 68-76.6 Minasny, B., A.B. McBratney. 2001. The Australian soil texture boomerang: a comparison of the Australian and USDA/FAO soil particle-size classification systems. Aust. J. Soil Res., 2001, 39, 14431451.7 Saxton, K.E. and W.J. Rawls. 2004. Soil water characteristic estimates by texture and organic matter for hydrologic solutions. Proc. , Soil Sci. Soc. Amer. International Conference, Seattle, WA Nov. 1-3, 2004.8 Cosby, B.J., G.M. Hornberger, R.B. Clapp, T.R. Ginn. 1984. A statistical exploration of the relationship of soil moisture characteristics to the physical properties of soils. Water Resources research, 20(6):682-690.9 Minasny, B., A.B. McBratney. 2000. Evaluation and development of hydraulic conductivity pedotransfer functions for Australian soil. Australian Journal of Soil Research, 38: 905-926.10 Puckett, W.E., J.H. Dane, B.F. Hajek. 1985. Physical and mineralogical data to determine soil hydraulic properties. SSSAJ, 49:831-836.11 Minasny, B., A.B. McBratney. 2000. Evaluation and development of hydraulic conductivity pedotransfer functions for Australian soil. Australian Journal of Soil Research, 38: 905-926.12 Rawls, W.J., L.R. Ahuja, D.L. Brakensiek. 1992. Estimating soil hydraulic properties from soils data. In: M. Th.van genuchten, F.J. Leij, and L.J. Lund (eds) Indirect Methods for Estimating the Hydrauli Properties of Unsaturated Soils. University of California, Riverside,329-340.13 Campbell, G.S., S. Shiozawa. 1992. Prediction of hydraulic properties of soils using PSD and bulk density data. pp. 317-328. In M. Th.van Genuchten et al. (eds) Proceedings of Inter. Workshop on Indirect Methods for Estimating the Hydraulic properties of Unsaturated Soils.14 Wosten, J.H.M. 1997. Pedotransfer functions to evaluate soil quality. In Gregorich E.G., Carter M.R. (eds), Soil Quality for Crop Production and Ecosystem Health. Development in Soils Science, Vol. 25, Elsevier, Amsterdam, pp. 221-245.15 Wosten, J.H.M., A. Lilly, A. Nemes, C. Le Bas. 1999. Development and use of a database of hydraulic properties of European soils. Geoderma, 90:169-185.16 Pachepsky, Y.A., D.J. Timlin, L.R. Ahuja. 1999. Estimating saturated soil hydraulic conductivity using water retention data and neural networks. Soil Science, 164(8): 552-560.17 Wosten J H M, Pachepsky Y A, Rawls W J. 2001. Pedo-transfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics. Journal of Hydrology . 251: 123150.Indirect methods to estimate saturated soil hydraulic conductivity in paddy soil of ChangshuZHANG Jun-hua1,2, LIU Jian-li1*, ZHU An-ning1(1. Institute of soil science, Chinese academy of sciences, Nanjing 210008,China; 2. Graduate University of Chinese Academy of Sciences, Beijing 100049,China)Abstract: As an alternative approach of direct experiment, using soil physic