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

薛占琪 2017051941原始数据：soil 1drywetsoil 2drywetsoil 2drywetθ(cm3/cm3)h(cm)h(cm)θ(cm3/cm3)h(cm)h(cm)θ(cm3/cm3)h(cm)h(cm)0.2714.4510.60.21164.8904.40.25341.56618.70.22592.6389.80.22776.95810.222519.42507.105080.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函数模拟结果如下模拟结果和实测结果见表一表一 VG土壤一观测值和拟合值表（脱湿）RETC模拟结果如下图（b）.BC函数模拟结果如下模拟结果和实测结果见表二表二 VG土壤一观测值和拟合值表（脱湿）RETC模拟结果如下图土壤一吸湿过程(a).VG函数模拟结果如下模拟结果和实测结果见表三表三 VG土壤一观测值和拟合值表（吸湿）RETC模拟结果如下图（b）.BC函数模拟结果如下模拟结果和实测结果见表四表四 RETC土壤一观测值和拟合值表（吸湿）RETC模拟结果如下图（2）土壤二土壤二脱湿过程（a）.VG函数模拟结果如下模拟结果和实测结果见表五表五 VG土壤二观测值和拟合值表（脱湿）RETC模拟结果如下图（b）.BC函数模拟结果如下模拟结果和实测结果见表六表六 BC土壤二观测值和拟合值表（脱湿）RETC模拟结果如下图土壤二吸湿过程(a).VG函数模拟结果如下模拟结果和实测结果见表七表七 VG土壤二观测值和拟合值表（吸湿）RETC模拟结果如下图（b）.BC函数模拟结果如下模拟结果和实测结果见表八表八 RETC土壤二观测值和拟合值表（吸湿）RETC模拟结果如下图（3）土壤三土壤三脱湿过程（a）.VG函数模拟结果如下模拟结果和实测结果见表九表九 VG土壤三观测值和拟合值表（脱湿）RETC模拟结果如下图（b）.BC函数模拟结果如下模拟结果和实测结果见表十表十 RETC土壤三观测值和拟合值表（脱湿）RETC模拟结果如下图土壤三吸湿过程(a).VG函数模拟结果如下模拟结果和实测结果见表十一表十一 VG土壤三观测值和拟合值表（吸湿）RETC模拟结果如下图（b）.BC函数模拟结果如下模拟结果和实测结果见表十二表十二 RETC土壤三观测值和拟合值表（吸湿）RETC模拟结果如下图4．三种函数拟合效果比较根据土壤水分曲线的散点图画出趋势线，再根据趋势线得出幂函数公式，就可以得到幂函数曲线；再在EXCEL表中绘制三种函数模拟得到的曲线，结果分别如下：5.三种函数拟合效果、实测值与计算值比较观察三种模型中的土壤脱湿曲线与吸湿曲线，并比较R2值可知，模型的拟合程度均达到99%以上,而用VG函数模拟的曲线与实测结果拟合程度最高，效果最佳，BC函数次之，幂函数最差。同时VG函数能更好的模拟曲线上高含水量低土壤吸力的部分，有明显的进气吸力,因为BC函数和幂函数在数学表达式上本身就类似，所以BC函数和幂函数则更加接近，且模拟区间偏向于曲线上低含水量高土壤吸力的部分。