土壤水分溶质动力学作业.doc

薛占琪 2017051941原始数据：soil 1drywetsoil 2drywetsoil 2drywet(cm/cm)h(cm)h(cm)(cm/cm)h(cm)h(cm)(cm/cm)h(cm)h(cm)0.2714.4510.60.21164.8904.40.25341.56618.70.22592.6389.80.22776.95810.222519.42507.10.24502307.50.24550397.80.2413481116.20.26431.4248.30.26405.6284.60.26787.15570.28374.4203.80.28307.62100.28489.4302.20.3326.9169.10.3237.7158.30.3318.4174.50.32286.4141.30.32185.6120.80.32213.8105.50.34250.9118.40.34145.592.70.34146.465.80.36219.1990.36113.670.90.36101.141.80.38190.182.30.3887.353.40.3869.326.60.4162.867.40.464.938.90.446.116.50.42136.553.90.424526.30.4228.59.60.4411041.10.4425.914.70.4414.34.41.用EXCEL文档绘制的3个土壤的土壤水分特征曲线的散点图如下：2.不同土壤的水分特征曲线做比较同：三种土壤的土壤水吸力均随土壤含水率的增大而减小异：在吸水和脱水过程中的水分特征曲线重合度越来越高，土壤1滞后现象最明显，土壤粘性较小，土壤3粘性最大3.分别用幂函数、BC函数和RETC软件拟合土壤水分特征曲线，再比较三个函数的拟合效果（1）土壤一土壤一脱湿过程（a）.VG函数模拟结果如下r=0，s=0.5186，=0.006，n=1.641，m=0.3906R2=0.99963433模拟结果和实测结果见表一表一 VG土壤一观测值和拟合值表（脱湿）RETC模拟结果如下图（b）.BC函数模拟结果如下 s=0.6685，=0.0223，n=0.4117，R2=0.98073343模拟结果和实测结果见表二表二 VG土壤一观测值和拟合值表（脱湿）RETC模拟结果如下图土壤一吸湿过程(a).VG函数模拟结果如下r=0.0899,s=0.4897，=0.014，n=1.6481，m=0.3932R2=0.99999994模拟结果和实测结果见表三表三 VG土壤一观测值和拟合值表（吸湿）RETC模拟结果如下图（b）.BC函数模拟结果如下 s=0.5411，=0.0413，n=0.3099，R2=0.55379471模拟结果和实测结果见表四表四 RETC土壤一观测值和拟合值表（吸湿）RETC模拟结果如下图（2）土壤二土壤二脱湿过程（a）.VG函数模拟结果如下r=0.1061，s=0.46，=0.0128，n=1.4887，m=0.3283R2=0.99999998模拟结果和实测结果见表五表五 VG土壤二观测值和拟合值表（脱湿）RETC模拟结果如下图（b）.BC函数模拟结果如下s=0.44，=0.0233，n=0.2302，R2=0.99253402模拟结果和实测结果见表六表六 BC土壤二观测值和拟合值表（脱湿）RETC模拟结果如下图土壤二吸湿过程(a).VG函数模拟结果如下r=0.1011，s=0.46，=0.0229，n=1.4238，m=0.2977R2=0.99999994模拟结果和实测结果见表七表七 VG土壤二观测值和拟合值表（吸湿）RETC模拟结果如下图（b）.BC函数模拟结果如下s=0.43，=0.0327，n=0.2236，R2=0.99553932模拟结果和实测结果见表八表八 RETC土壤二观测值和拟合值表（吸湿）RETC模拟结果如下图（3）土壤三土壤三脱湿过程（a）.VG函数模拟结果如下r=0.1122，s=0.46，=0.0284，n=1.2741，m=0.2151R2=0.99999998模拟结果和实测结果见表九表九 VG土壤三观测值和拟合值表（脱湿）RETC模拟结果如下图（b）.BC函数模拟结果如下s=0.43，=0.0336，n=0.1509，R2=0.99665362模拟结果和实测结果见表十表十 RETC土壤三观测值和拟合值表（脱湿）RETC模拟结果如下图土壤三吸湿过程(a).VG函数模拟结果如下r=0.1104,s=0.4597,=0.1065，n=1.2075，m=0.1718R2=0.99999919模拟结果和实测结果见表十一表十一 VG土壤三观测值和拟合值表（吸湿）RETC模拟结果如下图（b）.BC函数模拟结果如下s=0.43，=0.1008，n=0.1486，R2=0.99691566模拟结果和实测结果见表十二表十二 RETC土壤三观测值和拟合值表（吸湿）RETC模拟结果如下图4三种函数拟合效果比较根据土壤水分曲线的散点图画出趋势线，再根据趋势线得出幂函数公式，就可以得到幂函数曲线；再在EXCEL表中绘制三种函数模拟得到的曲线，结果分别如下：5.三种函数拟合效果、实测值与计算值比较观察三种模型中的土壤脱湿曲线与吸湿曲线，并比较R2值可知，模型的拟合程度均达到99%以上,而用VG函数模拟的曲线与实测结果拟合程度最高，