土壤水分溶质动力学课后复习

1、薛占琪 2017051941原始数据:soil 1drywetsoil 2drywetsoil 2drywet0 (cm3/cm3)h(cm)h(cm)0 (cm3/cm3)h(cm)h(cm)0 (cm3/cm3)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.2830

2、7.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.368887.353.40.3869.326.60.4162.867.40.464.938.90.42136.553.90.424526.30.4411041.10.4425

3、.914.70.4414.34.41.用EXCE文档绘制的3个土壤的土壤水分特征曲线的散点图如下:土壤水分特征曲线1田（力吸水壤土800700600500300100000.40.51400土壤含水率（cm3/cm ）soil 1(dry) dsoil 1(wet)土壤水分特征曲线2)me(力吸水壤土120010008006004000.40.5土壤含水率（cm3/cm）2000:soil 2(dry) J soil 2(wet)土壤水分特征曲线37000 十60005000(力吸水壤土4000300020001000十,0.40.5土壤含

4、水率(cm3/cm )soil 3(dry)广、soil 3(wet)2. 不同土壤的水分特征曲线做比较同：三种土壤的土壤水吸力均随土壤含水率的增大而减小异：在吸水和脱水过程中的水分特征曲线重合度越来越高，土壤1滞后现象最明显，土壤粘性较小，土壤3粘性最大3. 分别用幕函数、BC函数和RETC软件拟合土壤水分特征曲线，再比较三个函数的拟合效果(1) 土壤 土壤一脱湿过程(a) .VG函数模拟结果如下?= 0, ?= 0.5186, a= 0.006 , n = 1.641 , m = 0.3906? = 0.99963433模拟结果和实测结果见表一表一 VG 土壤一观测值和拟合值表(脱湿)HO

5、Plog-PWC-obaWC-fitHC-dev1.了144E十032 S539.2 000.1971.002 92592石E击032 7723E20O,2197 00033.5020E+032.7007.00.2411OOll勺.23丄坦E+032.6349.2 SOiiO.2 617-”00175.3744E十032.5733.2 50110.2ai 石-.0015石.3259E+032 - 5144-3000.3012-.001272S后4E+032.勺570 32010 32040004.2509E+CI33.3 事95.340110.339后-0004g.21-91E+032.340

6、i.3 SCiiiO.3Sfi5.OOll10.1301E+032.220-3900,3754,OOl石11162SE+032.2117.000.3 98 石.001415.dS占SE十032.1351.4二勤6,000413110QE牛032 0114.4400 4420-.0020RETC模拟结果如下图Hydraulic Properties:logh vs Theta0 +-Quo010 20_30.4OS 0_&Water Content M(b) .BC函数模拟结果如下?= 0.6685 , a= 0.0223 ,n = 0.4117, ? = 0.98073343模拟结果和实测结果

7、见表二NOPlog-EWOcb3TflC-fit闻匚一dev1.7144+033.9535.2000.213S-.01332.59比E+032.7723.2200.2 JOS.010旧3.5020E+032.7007.2400-.0072T.呼31E+032,53今9,2 石00,2631-,003152.5733.2 000.27S9.00116.326SE+032.51弓-3000-2950.00507.286E+032.且570.3200-3115.0085Q.2S0I9E+032.3995.3400-3239.01119.2191E+032.30石.3600 3号 76.012210.

8、1301E+032.2790.3800.36S7.0113111石28E+032.2117.40003D30.007012.1365E+032.1351.4200.422-.003石13 liaiOE+032.0414,4400-4618-.0213土壤一观测值和拟合值表（脱湿）表二VGRETC莫拟结果如下图Hydraulic Properties: log h vs. ThetaWater Content -土壤一吸湿过程(a).VG函数模拟结果如下?= 0.0899,?= 0.4897, a= 0.014, n = 1.6481 , m = 0.3932? = 0.99999994模拟结果

9、和实测结果见表三表三VG 土壤一观测值和拟合值表(吸湿)MOPlog-PwC-obsMC-fltWC-dev1510SE+032,7081*2Q002000*00002.3S9aE+032.S90e.22 00.3200.00003.3075+032.4S7S 2900 2900*00004”2注S3E+032.3950.2600.2600.00005.2033E+032,3032.2300.2300,0000石,1691E+032.2231,3000.3000.00007-1413E+032 .1501.3200.3200OQOO3.11S4E+032.0734400.000051.3955

10、.SMO.000010.3230E-021.915今.3800* 3Q00.000011-740i-021.327,4000,4000.000012-SSaOE+021,7316.4500.4200.000013-4110Z-021.0000RETC莫拟结果如下图Hydraulic Properties:log h vs. ThetaWater Content -1(b) .BC函数模拟结果如下?= 0.5411，a= 0.0413，n = 0.3099，? = 0.55379471模拟结果和实测结果见表四表四RETC 土壤一观测值和拟合值表（吸湿）HOFlog-PWC-obaWC-fitWC

11、-dev1.510E4032 JOSl.2000.2103-.01032-3898E+032.5 909-2200.2237-.008732.4S7S-2400-.00614-2453Z-032.3950-26 00*2630-.00305.203SE+O32.3 092.2SQ0.2795.00056.1691E+032.2251 3000.2962.00337iqi3E+O32.1501-32 DO.3131 0 五 93.1134Z+032*0734.34M.3303.00529 9900E*021.9356-3600.3436.010410.3230E+021*915弓.3300.37

12、02.10 OSS11-740E+021.8287-4000.3939.113.539CE0：1.7316-.002113MllOE-02:花132处00.4531-.0191RETC莫拟结果如下图Hydraulic Properties: log h vs. Theta巨邑-PST gun胡Mejd-爲空(2) 土壤二土壤二脱湿过程(a) .VG函数模拟结果如下?= 0.1061 , ?= 0.46, a= 0.0128, n = 1.4887, m = 0.3283? = 0.99999998模拟结果和实测结果见表五表五VG 土壤二观测值和拟合值表(脱湿)HOPlog-?HC-otsKC-

13、fiEtfC-dev1,iiesE+M3.0663 2 00C.254一,0石5今2.7769E+032.e?04.2 300.2朋5-.04653.5500E+032.7Q04.2-IOC.26S7-.02374.405石E+032 6OS1.2 600.2727-412 7.3076E+032.SSC.2 500.2*757.00056.2377E+032.3760.3000.2S11.008971S56E+032.26S6.3：0G.305 5.01123,1吟55E+032.1629.300.333 4.G06S9.1136E+032.0554.3600.3总30-.003010.S7

14、30E+Q219吒10.3S0.3917-.011711.6409E+043.eo&s.弓 000.2647.135312G500E+021.6552.42 00.4245-.004513.2 59I0E+O21.133.010.42S6.0104RETC莫拟结果如下图Hydraulic Properties: log h vs. ThetaWater Content 卜(b) .BC函数模拟结果如下?= 0.44, a= 0.0233 , n = 0.2302 , ? = 0.99253402模拟结果和实测结果见表六NOPlog-FWC-Ot3wc-fitWC-dev1.1155E+043.

15、0663.2 682-.0322.7759E+032.3904*2200*2700-.05003.3500E+032.7404.24002T25-.032岂4.eoS6E+a32 60S1*2eoo* 27 石9-.01石95.3076E+032.4SS0-zSCO.2530-.00306.2377E + CI32.3760*3000.2916,003 些7.1SS6E + 032.26?fi.3200,3042.0155o.1455E + 032,1629.3400.3225,01759.1136E+032.0554.3600*35.101.009910.Q73CIE + 021.9410.

16、3000.3947-.014711.6403E+0T3.3063.2663.1337丄2.4500E+021.5532.4200.4500-.010013.2530E+021.4133.0100土壤二观测值和拟合值表（脱湿）表六BCRETC莫拟结果如下图Hydraulic Properties: log h vs. Theta1_0 +0.0C_10.20-30.40_5Water Content -土壤二吸湿过程(a).VG函数模拟结果如下?= 0.1011 , ?= 0.46, a= 0.0229, n = 1.4238, m = 0.2977? = 0.99999994模拟结果和实测结果

17、见表七表七VG 土壤二观测值和拟合值表(吸湿)NOPlog-PWC-ohsWC-fTtWC-dev12.3564,2 050.2652-.06522SSIOE+032.72.22 00.2663-.04BS3.3973E+032.5997.220-02S2q2e今6E+032.4542.2600.2726-,012&5.OOE+032.3222.2 550.2794.000若6153+口32.-93?.3Q0Q.00347.12 0i8E+032.OS21.32 0110.30S2.Oils3.9270E+021.367L3g.3331.00599.7090E+021.S306.3 600.3

18、631-0O3110.S540E井051.72?S.3SOiO.3531-.012111-3809E+03 .SSOS.4000-2ES2.13弓S12.2630E+021.4200.4200.4250.0050131吩7QE+口m73.a 二05RETC模拟结果如下图Hydraulic Poperties:log h vs. Theta-u旦 pmeH gun SSL )汙一Waler Content -(b) .BC函数模拟结果如下?= 0.43, a= 0.0327 , n = 0.2236 , ? = 0.99553932模拟结果和实测结果见表八NOPlog-F闻C-ofc3枇-fit

19、KC-dev1.9044E+032.9554.2000.27S-.06732.581GE+032,76气2.22002石97.0弓973,357SE+032.5997.2725-.0325弓2.4542 2 eg-.0135.SlOGE+032.32222SOO2S30-.003061533E+032.1995,3000,2319.00517,120aE+032.0521.3200,304.C154e9270E+021.9671-3900.3230.01709,7D9GE+021.8506.3600,3506.009410，534C+021.7275.3943-.01311,3203E+03.

20、5202 4000.133912丄4200-4200.4300-.010013*1470E+021.1673.4400*4300-0100RETC 土壤二观测值和拟合值表（吸湿）表八RETC模拟结果如下图Hydraulic Poperties:log h vs. ThetaWater Content -1(3) 土壤三土壤三脱湿过程(a) .VG函数模拟结果如下?= 0.1122 , ?= 0.46, a= 0.0284, n = 1.2741 , m = 0.2151? = 0.99999998模拟结果和实测结果见表九表九VG 土壤三观测值和拟合值表(脱湿)MOFlog-FWC-ot3WCf

21、itWC-dev1.53ZE+03.7277.2 000.2000.OOOiOn.5519iE + 043.4013.2200士 00.OOOlO3,133E+03.12 沪.2400 OOOO.7S71E + 032.SSSO.2 6010.2 eoc 00 005.4394E+032.G37.2S0nQ.2300.000063ia弓E+032.5030.3000 3000.00007.2138E+032.3300.32 00.320110.com.1464E+032.1655.300.340ii0.00003-lOllE+032 343.3600-3600-00 0010.6930E+02

22、1.5407.3S0nQ.300.000011,马tlOE+Q21.6637弓000.oom12,2S50E+021.4548.42 00.420.000013.1430E+021.1533.44 0ii0.4400.0000RETC莫拟结果如下图Hydraulic Poperties: log h vs- ThetaWatsr ContAnt -(b) .BC函数模拟结果如下?= 0.43, a= 0.0336 , n = 0.1509 , ? = 0.99665362模拟结果和实测结果见表十表十RETC 土壤三观测值和拟合值表(脱湿)KOPlog-FWC-0fc3低-litWC-dev1.

23、5342E+043-7277 2000.1965.00352.2513E+。理3”013,2200 2201-,00013.134SE+03,1297.2400.2419-.0i019些.7e71E+032-3960.200.如4-.0i02 45.4E94E+Q32.6597.2300 2315-.0019.31S4E+D32-5030.3000-30DS-OQOS7.2=35E+032.33DQ.3200313豆,OQOS8.16E+032-1655,3400-3382.00189.lOllE+032-00483石00-357S.002 410 6930E+021-S40T.3300,37

24、55*001 弓11,q610E+021.6637.4000.4026-.002615.5S50E+021-454S.42C0.4300-.OlCO13.1430E+021.1553*4400.4300*0100RETC模拟结果如下图Hydraulic Properties: log h vs. ThetaWater Contenl |土壤三吸湿过程(a).VG函数模拟结果如下=0.1718?= 0.1104,?= 0.4597, a= 0.1065 , n= 1.2075 , m? = 0.99999919模拟结果和实测结果见表十表一 VG 土壤三观测值和拟合值表(吸湿)NOPiog-PWC

25、-obsWC-fitWC-dev1.55二9E十023.32032OQ0.tOiQO,00002.2507E-043.35A3.2300,5 200.00003.llLfiE亠043.0477.2400.2400.0000住.753*00305.3022E-032.4S03.2500.2300.00006.174SE+032.410 3000_3000.000072,02334叮00e1.S1S3.3400,3400.00009.4130E+031.6212,3500-3600.000010.Z66OE-O21-423* 3500.3500400011.1650E-021.2175.400 D.4001-.0031129石OOEVOl.982 3 200弓 198,000213.aiooE-01.643 5 4400.401-.OOOlRETC莫拟结果如下图Hydraulic Properties: log h vs. ThetaWater Content -(b) .BC函数模拟结果如下?= 0.43, a = 0.1008 , n = 0.1486 , ? = 0.99691566模拟结果和实测结果见表十二表十二RETC 土壤三观测值和拟合值表(吸湿)OfiSfiEvedfluted dat