土壤水分动力学作业

1、土壤水分动力学作业用RETC软件中van Genuchten (1980)模型、Brooks and Corey (1964) 模型以及幂函数模型拟合下列三种土壤的水分特征曲线，用EXCEL作图比较水分特征曲线的拟合和实测结果。soil 1drywetsoil 2drywetsoil 3drywet(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.9581.00.222519.42507.10.24502

2、.0307.50.24550.0397.80.241348.01116.20.26431.4248.30.26405.6284.60.26787.1557.00.28374.4203.80.28307.6210.00.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.199.00.36113.670.90.36101.141.80.38190.182.30.3887

3、.353.40.3869.326.60.4162.867.40.464.938.90.446.116.50.42136.553.90.4245.026.30.4228.59.60.44110.041.10.4425.914.70.4414.34.41.土壤一1）脱湿过程a. VG函数模拟结果如下：r=0，s =0.5186，=0.006，n=1.6411，R2=0.9996，m=0.390634模拟结果和实测结果见表1： 表1 VG土壤1观测值与拟合值表(脱湿)NO P log-P WC-obs WC-fit WC-dev 1 .7144E+03 2.8539 .2000 .1971 .002

4、9 2 .5926E+03 2.7728 .2200 .2197 .0003 3 .5020E+03 2.7007 .2400 .2411 -.0011 4 .4314E+03 2.6349 .2600 .2617 -.0017 5 .3744E+03 2.5733 .2800 .2816 -.0016 6 .3269E+03 2.5144 .3000 .3012 -.0012 7 .2864E+03 2.4570 .3200 .3204 -.0004 8 .2509E+03 2.3995 .3400 .3396 .0004 9 .2191E+03 2.3406 .3600 .3589 .00

5、1110 .1901E+03 2.2790 .3800 .3784 .001611 .1628E+03 2.2117 .4000 .3986 .001412 .1365E+03 2.1351 .4200 .4196 .000413 .1100E+03 2.0414 .4400 .4420 -.0020RETC模拟结果图1： 图1 VG函数模拟土壤1水分特征曲线（脱湿）b.BC函数函数模拟结果如下：s=0 .5758439739 = 0. 01554 2608 n= 0.41167 R2=模拟结果和实测结果见表2:表2 BC土壤1观测值与拟合值表（脱湿）NO P log-P WC-obs WC-

6、fit WC-dev1 .7144E+03 2.8539 .2000 .2138 -.0138 2 .5926E+03 2.7728 .2200 .2309 -.0109 3 .5020E+03 2.7007 .2400 .2472 -.0072 4 .4314E+03 2.6349 .2600 .2631 -.0031 5 .3744E+03 2.5733 .2800 .2789 .0011 6 .3269E+03 2.5144 .3000 .2950 .0050 7 .2864E+03 2.4570 .3200 .3115 .0085 8 .2509E+03 2.3995 .3400 .3

7、289 .0111 9 .2191E+03 2.3406 .3600 .3478 .0122 10 .1901E+03 2.2790 .3800 .3687 .0113 11 .1628E+03 2.2117 .4000 .3930 .0070 12 .1365E+03 2.1351 .4200 .4226 -.0026 13 .1100E+03 2.0414 .4400 .4618 -.0218RETC模拟结果见图2： 图2 BC函数模拟土壤1水分特征曲线（脱湿）2）吸湿过程a. VG函数模拟结果如下:r= 0.09006 s=0.48969 =0.01399 n=1.64868 ,m=0.

8、393454, R2=0.9999模拟结果和实测结果见表3： 表3 VG土壤1观测值与拟合值表(吸湿)NO P log-P WC-obs WC-fit WC-dev 1 .5108E+03 2.7083 .2000 .2000 .0000 2 .3898E+03 2.5908 .2200 .2200 .0000 3 .3075E+03 2.4878 .2400 .2400 .0000 4 .2483E+03 2.3950 .2600 .2600 .0000 5 .2038E+03 2.3092 .2800 .2800 .0000 6 .1691E+03 2.2281 .3000 .3000 .

9、0000 7 .1413E+03 2.1501 .3200 .3200 .0000 8 .1184E+03 2.0734 .3400 .3400 .0000 9 .9900E+02 1.9956 .3600 .3600 .0000 10 .8230E+02 1.9154 .3800 .3800 .0000 11 .6740E+02 1.8287 .4000 .4000 .0000 12 .5390E+02 1.7316 .4200 .4200 .0000 13 .4110E+02 1.6138 .4400 .4400 .0000RETC模拟结果见图3： 图3 VG函数模拟土壤1水分特征曲线（吸

10、湿） b. BC函数模拟结果如下:s=0 .48969 = 0.01399 n= 1.64868模拟结果和实测结果见表4： 表4 BC土壤1观测值与拟合值表(吸湿)NO P log-P WC-obs WC-fit WC-dev 1 .5108E+03 2.7083 .2000 .2103 -.0103 2 .3898E+03 2.5908 .2200 .2287 -.0087 3 .3075E+03 2.4878 .2400 .2461 -.0061 4 .2483E+03 2.3950 .2600 .2630 -.0030 5 .2038E+03 2.3092 .2800 .2795 .00

11、05 6 .1691E+03 2.2281 .3000 .2962 .0038 7 .1413E+03 2.1501 .3200 .3131 .0069 8 .1184E+03 2.0734 .3400 .3308 .0092 9 .9900E+02 1.9956 .3600 .3496 .0104 10 .8230E+02 1.9154 .3800 .3702 .0098 11 .6740E+02 1.8287 .4000 .3939 .0061 12 .5390E+02 1.7316 .4200 .4221 -.0021 13 .4110E+02 1.6138 .4400 .4591 -.

12、0191RETC模拟结果见图4： 图4 BC函数模拟土壤1水分特征曲线（吸湿）3) .EXCEL作图比较水分特征曲线的拟合和实测结果a.VG函数模拟结果与实测值比较：b.BC函数模拟结果与实测值比较：c.幂函数模拟结果与实测值比较:2.土壤二1）脱湿过程a. VG函数模拟结果如下:r= 0.10614 s=0 .45999 = 0 .01282 n= 1.48868 m=0.328263965 R2=0.99994模拟结果和实测结果见表5： 表5 VG土壤2观测值与拟合值表(脱湿) NO P log-P WC-obs WC-fit WC-dev 1 .1165E+04 3.0663 .2000

13、 .2000 .0000 2 .7769E+03 2.8904 .2200 .2200 .0000 3 .5500E+03 2.7404 .2400 .2400 .0000 4 .4056E+03 2.6081 .2600 .2600 .0000 5 .3076E+03 2.4880 .2800 .2800 .0000 6 .2377E+03 2.3760 .3000 .3000 .0000 7 .1856E+03 2.2686 .3200 .3200 .0000 8 .1455E+03 2.1629 .3400 .3400 .0000 9 .1136E+03 2.0554 .3600 .36

14、00 .0000 10 .8730E+02 1.9410 .3800 .3800 .0000 11 .6490E+02 1.8122 .4000 .4000 .0000 12 .4500E+02 1.6532 .4200 .4200 .000013 .2590E+02 1.4133 .4400 .4400 .0000RETC模拟结果见图8：图8 VG函数模拟土壤2水分特征曲线（脱湿）b.BC函数模拟结果如下:s=0 .43000 = 0 .01909 n= 0.24355 R2=0.99523339模拟结果和实测结果见表6： 表6 BC土壤2观测值与拟合值表(脱湿) NO P log-P WC

15、-obs WC-fit WC-dev 1 .1165E+04 3.0663 .2000 .2020 -.0020 2 .7769E+03 2.8904 .2200 .2229 -.0029 3 .5500E+03 2.7404 .2400 .2425 -.0025 4 .4056E+03 2.6081 .2600 .2612 -.0012 5 .3076E+03 2.4880 .2800 .2794 .0006 6 .2377E+03 2.3760 .3000 .2975 .0025 7 .1856E+03 2.2686 .3200 .3160 .0040 8 .1455E+03 2.1629

16、 .3400 .3353 .0047 9 .1136E+03 2.0554 .3600 .3561 .0039 10 .8730E+02 1.9410 .3800 .3797 .0003 11 .6490E+02 1.8122 .4000 .4081 -.0081 12 .4500E+02 1.6532 .4200 .4300 -.0100 13 .2590E+02 1.4133 .4400 .4300 .0100RETC模拟结果见图9： 图9 BC函数模拟土壤2水分特征曲线（脱湿）2)吸湿过程a.VG函数模拟结果如下:r= 0.10106 s=0 .45998 = 0 .02291 n= 1

17、.42381 m=0.297659 R2=0.99994模拟结果和实测结果见表7: 表7 VG土壤2观测值与拟合值表(吸湿) NO P log-P WC-obs WC-fit WC-dev 1 .9044E+03 2.9564 .2000 .2000 .0000 2 .5810E+03 2.7642 .2200 .2200 .0000 3 .3978E+03 2.5997 .2400 .2400 .0000 4 .2846E+03 2.4542 .2600 .2600 .0000 5 .2100E+03 2.3222 .2800 .2800 .0000 6 .1583E+03 2.1995 .

18、3000 .3000 .0000 7 .1208E+03 2.0821 .3200 .3200 .0000 8 .9270E+02 1.9671 .3400 .3400 .0000 9 .7090E+02 1.8506 .3600 .3600 .0000 10 .5340E+02 1.7275 .3800 .3800 .0000 11 .3890E+02 1.5899 .4000 .4000 .0000 12 .2630E+02 1.4200 .4200 .4200 .0000 13 .1470E+02 1.1673 .4400 .4400 .0000RETC模拟结果见图10： 图10 VG函

19、数模拟土壤2水分特征曲线（吸湿）b.BC函数模拟结果如下:s=0 .44000 = 0 .04035 n=0 .21198模拟结果和实测结果见表8： 表8 BC土壤2观测值与拟合值表(吸湿)NO P log-P WC-obs WC-fit WC-dev 1 .9044E+03 2.9564 .2000 .2053 -.0053 2 .5810E+03 2.7642 .2200 .2254 -.0054 3 .3978E+03 2.5997 .2400 .2443 -.0043 4 .2846E+03 2.4542 .2600 .2623 -.0023 5 .2100E+03 2.3222 .2

20、800 .2797 .0003 6 .1583E+03 2.1995 .3000 .2970 .0030 7 .1208E+03 2.0821 .3200 .3145 .0055 8 .9270E+02 1.9671 .3400 .3327 .0073 9 .7090E+02 1.8506 .3600 .3521 .0079 10 .5340E+02 1.7275 .3800 .3739 .0061 11 .3890E+02 1.5899 .4000 .3999 .0001 12 .2630E+02 1.4200 .4200 .4345 -.0145 13 .1470E+02 1.1673 .

21、4400 .4400 .0000RETC模拟结果见图11图11 BC函数模拟土壤2水分特征曲线（吸湿）3) EXCEL作图比较水分特征曲线的拟合和实测结果a.VG函数模拟结果与实测值比较：b.VG函数模拟结果与实测值比较：c.幂函数模拟结果与实测值比较：3.土壤三1）脱湿过程a.VG函数模拟结果如下:r= 0.11223 s=0.46003 = 0.02844 n= 1.27406 m=0.21297025 R2= 0.99999998模拟结果和实测结果见表9: 表9 VG土壤3观测值与拟合值表(脱湿)NO P log-P WC-obs WC-fit WC-dev 1 .5342E+04 3.

22、7277 .2000 .2000 .0000 2 .2519E+04 3.4013 .2200 .2200 .0000 3 .1348E+04 3.1297 .2400 .2400 .0000 4 .7871E+03 2.8960 .2600 .2600 .0000 5 .4894E+03 2.6897 .2800 .2800 .0000 6 .3184E+03 2.5030 .3000 .3000 .0000 7 .2138E+03 2.3300 .3200 .3200 .0000 8 .1464E+03 2.1655 .3400 .3400 .0000 9 .1011E+03 2.0048

23、 .3600 .3600 .0000 10 .6930E+02 1.8407 .3800 .3800 .0000 11 .4610E+02 1.6637 .4000 .4000 .0000 12 .2850E+02 1.4548 .4200 .4200 .0000 13 .1430E+02 1.1553 .4400 .4400 .0000RETC模拟结果见图15： 图15 VG函数模拟土壤3水分特征曲线（脱湿）b.BC函数模拟结果如下： s=0 .44000 = 0 .04238 n= 0.14688 R2= 0.99665362模拟结果和实测结果见表10： 表10 BC土壤3观测值与拟合值表

24、(脱湿)NO P log-P WC-obs WC-fit WC-dev 1 .5342E+04 3.7277 .2000 .1984 .0016 2 .2519E+04 3.4013 .2200 .2216 -.0016 3 .1348E+04 3.1297 .2400 .2429 -.0029 4 .7871E+03 2.8960 .2600 .2629 -.0029 5 .4894E+03 2.6897 .2800 .2819 -.0019 6 .3184E+03 2.5030 .3000 .3002 -.0002 7 .2138E+03 2.3300 .3200 .3183 .0017

25、8 .1464E+03 2.1655 .3400 .3365 .0035 9 .1011E+03 2.0048 .3600 .3553 .0047 10 .6930E+02 1.8407 .3800 .3756 .0044 11 .4610E+02 1.6637 .4000 .3988 .0012 12 .2850E+02 1.4548 .4200 .4280 -.0080 13 .1430E+02 1.1553 .4400 .4400 .0000RETC模拟结果见图16： 图16 BC函数模拟土壤3水分特征曲线（脱湿）2)吸湿过程a.VG函数模拟结果如下：r= 0.11016 s=0.459

26、82 = 0.10694 n= 1.20705 R2=0.99999902模拟结果和实测结果见表11 表11 VG土壤3观测值与拟合值表(吸湿)NO P log-P WC-obs WC-fit WC-dev 1 .6619E+04 3.8208 .2000 .2000 .0000 2 .2507E+04 3.3992 .2200 .2200 .0000 3 .1116E+04 3.0477 .2400 .2400 .0000 4 .5570E+03 2.7459 .2600 .2600 .0000 5 .3022E+03 2.4803 .2800 .2800 .0000 6 .1745E+03 2.2418 .3000 .3000 .0000 7 .1055E+03 2.0233 .3200 .3200 .0000 8 .6580E+02 1.8182 .3400 .3400 .0000 9 .4180E+02 1.6212 .3600 .3600 .0000 10 .2650E+02 1.4232 .3800 .3801 -.0001 11 .1650E+02 1.2175 .4000 .4000 .0000 12 .9600E+01 .9823 .4200 .4198 .0002 13 .4400E+01 .6435 .4400 .4401 -.000