土壤水力参数 matlab,一种由土壤剖面含水率估算土壤水力参数的方法

1研究背景在非饱和带中,负压和水力传导度均是土壤含水率的非线性函数,vanGenuchten建立的土壤负压和水力传导度与土壤含水率的关系式在土壤水分模拟中得到广泛应用,为表达方便,称vanGenuchten表达式中的参数为土壤水力参数。目前,确定这些参数的方法有直接拟合法、间接确定法及其他方法。直接拟合法采用试验得到的土壤含水率负压数据和土壤含水率水力传导度数据,通过数学优化方法求参。在国外,无论是试验测量还是数学优化方面的研究都已比较成熟,并开发了相关软件,如美国盐土实验室开发的RETC[1]软件。有学者提出在水力传导度数据缺失情况下推求完整土壤水力参数的方法[2]。目前,国内研究主要集中于优化方法选择及软件拟合结果分析。李春友[3]、廖林仙[4]、许小健[5]等分别将单纯形调优法、遗传算法、差异演化算法用于拟合土壤水力参数,并阐述了各数学方法的特点;范严伟等[6]运用软件RETC和MATLAB对土壤水力参数进行了拟合和对比分析。以上方法直接通过试验数据拟合参数,精度较高。但非饱和水力传导度试验难度很大,数据可信度很差[2]。由于土壤负压测试手段的问题,常常较难得到高精度的测量结果。间接确定法通过土壤颗粒级配和容重等特性来间接估计土壤水力参数。虽然间接确定法得到的参数准确度不如试验数据拟合法,但其所需的土壤信息容易获取,在实践中得到广泛应用。美国盐土实验室开发的ROSETTA[7]软件通过输入土壤特性参数,以土壤传递函数(pedotransferfunction)来估计土壤水力参数。刘建立等将间接确定法推求土壤水力参数划分为土壤转换函数法、物理经验法、分形理论[8]3种方法,并从这3种角度出发,总结了近20年来国内外土壤水分运动参数的间接推求方法。其他方法通过土壤水分运动规律反推参数。瞬时剖面法[9]利用含水率和土壤基质吸力的时空分布反求非饱和水力传导度。邵明安[10]等采用入渗试验的吸湿率和湿润锋运移速度等数据来推求土壤水力参数,并选取6种土样做了验证。土壤水分是农田水分管理、区域水文条件研究和流域水量平衡计算的重要参量[11]。土壤含水率数据作为直观反映土壤水分的指标,得到了广泛的应用。本文根据瞬时剖面法的思想,在入渗或蒸发条件下,利用土壤含水率的测量数据,通过非饱和土壤渗流的达西定律,推求饱和水力传导度和VanGenuchten模型中的非饱和土壤水力参数。2土壤水力参数优化方程在描述非饱和土壤水力性质的模型中,应用最广泛的是VanGenuchten模型。VanGenuchten模型[12]将土壤水分特征曲线描述如下:Se=1+(h)n-m(1)模型将非饱和水力传导度[13]描述为:K=KsSe1/21()1-Se1/mm2(2)式中:Se为土壤有效含水率[],Se=-rs-r;为土壤体积含水率[L3L-3],r为残余体积含水率[L3L-3],s为饱和体积含水率[L3L-3];h为土壤基质吸力[L],h0;K为非饱和水力传导度[LT-1];Ks为饱和水力传导度[LT-1];[L-1]、n[]和m[]均为经验常数,其中m=1-1/n,0