劉 暢，陳俊英,2，蔡耀輝，李志軍,2，Leionid Gillerman

(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，陜西 楊凌 712100；2. 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室，陜西 楊凌 712100；3.本古里安大學(xué)Blaustein沙漠研究所，以色列 思德博克 84990)

0 引 言

土壤斥水性是指水分不能或很難濕潤土壤顆粒表面的物理現(xiàn)象，主要表現(xiàn)為水滴在斥水性土壤表面不能迅速入滲或鋪展的現(xiàn)象。世界多地發(fā)現(xiàn)了斥水土壤，其分布范圍廣泛，涉及土壤類型多。斥水土壤入滲能力低，很容易形成優(yōu)先流，導(dǎo)致土壤水分分布不均勻，使水中攜帶的溶質(zhì)更快地進(jìn)入地下水，降低土壤持水能力，形成沖刷，影響農(nóng)作物產(chǎn)量[1]。研究斥水土壤的土壤水分特征曲線，對斥水土壤水分運(yùn)動規(guī)律具有重要的理論意義。

土壤水分特征曲線(簡稱土-水曲線)是表示土壤水吸力和土壤含水率的關(guān)系曲線，反應(yīng)土壤水的能量和數(shù)量之間的關(guān)系[2]，是研究土壤水分運(yùn)動以及溶質(zhì)運(yùn)移的基礎(chǔ)，對評價土壤水分的有效性和持水性具有重要意義[3]。對于普通親水土壤以及含有混摻物親水土壤的土-水曲線，在國內(nèi)外已經(jīng)取得了很多的研究成果[4-6]，常用的土-水曲線擬合模型主要包括Brooks-Coery(BC)模型[7]、van Genuchten(VG)模型[8]、Dual-porosity(DP)模型[9]和Lognormal distribution(LND)模型[10]，一般都是基于實(shí)測吸力和土壤含水率數(shù)據(jù)，采用已有的數(shù)學(xué)模型進(jìn)行擬合，推求水力參數(shù)，評價土壤持水性能[11]。由于斥水土壤的水分運(yùn)動比親水土壤的緩慢，已有的親水土壤水力參數(shù)不能適用于一定條件下的斥水土壤，因此需要進(jìn)一步研究。斥水土壤的土-水曲線研究起步較晚，Batuters等[12]對人工配置的斥水砂土進(jìn)行入滲試驗(yàn)，研究發(fā)現(xiàn)斥水砂土的土-水曲線和斥水程度有關(guān)；Arey等[13]應(yīng)用毛管上升過程及相似理論，基于VG模型得出了斥水土壤的水力性質(zhì)模型，但文章指出模型的適用性還不明確；劉世賓等[14]測定了4種人工添加斥水劑配制土壤的土壤水分特征曲線，并用VG和BC模型進(jìn)行擬合和分析。國內(nèi)外對于斥水土壤水分特征曲線的研究雖然已經(jīng)取得了很多成果，但是對于不同天然斥水土壤，以及同類斥水土壤在不同斥水程度下的土-水曲線研究還不多見，因此本研究選用以色列3個質(zhì)地、2個斥水程度的6種斥水土壤，測定其吸力和土壤含水率，應(yīng)用RETC軟件中的VG、BC、DP和LND模型進(jìn)行擬合，并對擬合效果進(jìn)行評價，給出4個土壤水分特征曲線模型的適用性。

1 材料與方法

1.1 試驗(yàn)材料

研究所用土壤樣品取自以色列的基布茲Beery、基布茲Bitzaron和基布茲Magen三個地方不同質(zhì)地的表層土壤，共6種土樣，表1所示為各個土樣的編號代碼、土壤類型和相關(guān)參數(shù)。黏土L1和砂土M1的斥水等級為嚴(yán)重斥水，其余土壤樣品的斥水等級均為中等斥水，其中壤土B1和B2的斥水等級雖然相同，但B1的斥水持續(xù)時間比B2大。

表1 土壤采集地及土壤相關(guān)參數(shù)Tab.1 Soil spots and the related parameters

1.2 試驗(yàn)方法

將采集的土壤樣品自然風(fēng)干，去雜過2 mm標(biāo)準(zhǔn)孔篩篩子，按土壤采集地的實(shí)際容重(6種試驗(yàn)土樣容重均為1.30 g/cm3)裝入環(huán)刀。試驗(yàn)開始前將環(huán)刀置于蒸餾水中浸泡至飽和，試驗(yàn)結(jié)束后將環(huán)刀置于105 ℃烘箱內(nèi)干燥至恒重，以計算土壤含水率，各處理均重復(fù)3次，取其均值作為最終結(jié)果。

采用日本CR21GⅡ型高速恒溫冷凍離心機(jī)測定土壤水分曲線，測定時機(jī)內(nèi)恒溫4 ℃，將待測樣品放入離心機(jī)裝置中，加壓范圍為10～8 000 cm，隨著施加壓力增大，其平衡時間也增加，每次離心結(jié)束后，采用電子天平(ES-3002H型)稱重，獲得土壤質(zhì)量含水率，再將其轉(zhuǎn)換為體積含水率。

1.3 土壤水分特征曲線擬合模型

應(yīng)用美國國家鹽改中心(US Salinity Laboratory)提供的RETC軟件中不同的土壤水分特征曲線模型，擬合實(shí)測的試驗(yàn)數(shù)據(jù)，分析或預(yù)測斥水土壤的水力性質(zhì)。本文選擇求解土壤非飽和導(dǎo)水率K的Mualem模型，應(yīng)用RETC軟件中BC、VG、DP和LND模型對不同類型斥水土壤的實(shí)測土壤水分特征曲線進(jìn)行擬合，確定土壤水分特征曲線模型的參數(shù)，通過模型計算得出實(shí)測土壤水吸力對應(yīng)的含水率，并與實(shí)測值進(jìn)行對比。

(1)BC模型。

(1)

式中：Se是土壤的飽和度；θ是土壤的體積含水率；θr是土壤殘余體積含水率；θs是土壤飽和體積含水率；h是壓力水頭；λ是土壤孔隙尺寸的分布參數(shù)，與土-水曲線的斜率有關(guān)；α是進(jìn)氣值的倒數(shù)。

(2)VG模型。

(2)

式中：m、n是影響土壤水分特征曲線形態(tài)的經(jīng)驗(yàn)參數(shù)，m=1-1/n(n>1)，其余符號含義和式(1)相同。

(3)DP模型。

Se=w1[1+(α1h)n1]-m1+w2[1+(α2h)n2]-m2

(3)

式中：w1和w2分別是2個區(qū)域的權(quán)重因子；α1和α2分別是2個區(qū)域進(jìn)氣值的倒數(shù)，m1和m2、n1和n2分別是2個區(qū)域土壤孔隙尺寸的分布參數(shù)，其余符號含義與式(1)和(2)相同。

(4)LND模型。

(4)

式中：h0、σ等同于前述公式的α、n，分別是進(jìn)氣值的倒數(shù)和土壤孔隙尺寸的分布參數(shù)，其余符號含義與(1)和(2)式相同。

2 結(jié)果與分析

2.1 不同斥水土壤的土壤水分特征曲線

圖1為試驗(yàn)測得的不同斥水土壤水分特征曲線。由圖1可知，三種斥水土壤的水分特征曲線在形態(tài)上與親水土壤的較為一致，其變化趨勢也與親水土壤比較類似。當(dāng)吸力較低(S<1 000 cm)時，隨著吸力的增加，三種類型斥水土壤含水率下降速度較快，這主要是在低吸力時，土體通過大孔隙進(jìn)行排水，即使吸力變化不大，三種類型斥水土壤的含水率也會發(fā)生明顯改變，因此在低吸力時，斥水土壤的土-水曲線變化平緩；同時可以看出在相同吸力情況下，黏土L1、L2的含水率高于壤土和砂土，同時砂土的含水率又遠(yuǎn)低于黏土和壤土兩種類型土壤。當(dāng)吸力較高(s=1 000-8 000 cm)時，三種類型的斥水土壤含水率下降速度較慢，其土-水曲線呈陡直狀，在高吸力階段，隨著吸力的增加，斥水土壤只有較小的孔隙能保留水分，其土體對其吸持力較大，因此這一階段的土壤含水率表現(xiàn)出隨吸力增加無顯著變化；同時高吸力階段，在吸力相同的情況下，三種類型土壤含水率大小順序仍與低吸力階段相同，進(jìn)一步說明斥水土壤的土-水曲線與親水土壤的土-水曲線變化規(guī)律基本一致。

圖1 不同斥水土壤水分特征曲線Fig.1 Soil water characteristic curves of different repellent soils

從圖1還可以看出，對于同一種類型的土壤，在同一吸力條件下，土壤斥水程度等級越高，含水率越小，保水能力降低。這是因?yàn)槌馑寥烙信懦馑奶攸c(diǎn)，使得土壤對水的吸力減小，隨著外界吸力增大，會在土壤垂直剖面容易形成優(yōu)先流現(xiàn)象[15]，因此土壤斥水持續(xù)時間越大，土壤持水能力也越差，這與Ritsema等[16]關(guān)于斥水土壤水分運(yùn)移研究中的結(jié)果相同。

2.2 不同類型土壤水分曲線模型及評價分析

各斥水土樣的不同模型擬合統(tǒng)計特征值如表2所示。從表2可以看出，VG、BC和DP模型對三種類型的斥水土壤水分特征曲線的擬合度均較好，相關(guān)系數(shù)都在0.97以上，呈現(xiàn)極顯著相關(guān)關(guān)系(p<0.000 1)，殘差平方和小于0.005，所以VG、BC和DP模型均能較好的適用于斥水黏土、壤土和砂土。而LND模型在擬合黏土和壤土?xí)r相關(guān)系數(shù)均小于0.20(本文沒有列出)，因此LND模型不能應(yīng)用于斥水黏土和壤土，但對于斥水砂土M1和M2，LND模型擬合的相關(guān)系數(shù)接近1(分別為0.998和0.999)，殘差平方和幾乎為零，呈現(xiàn)極顯著相關(guān)關(guān)系(p<0.000 1)，說明LND模型能夠很好地適用于斥水砂土。

不同模型擬合斥水土壤的水分特征曲線效果如圖2所示。從圖2可以看出在高吸力到低吸力的過渡段，砂土M1、M2和壤土B2的VG模型擬合曲線與砂土M1、M2的BC模型擬合曲線均存在明顯的拐點(diǎn)，而黏土L1、L2和壤土B1的VG、BC模型擬合曲線與3種質(zhì)地土壤的DP模型擬合土-水曲線均相對比較平順，同時砂土M1、M2的LND模型擬合的土-水曲線與實(shí)測值擬合度較高。從圖2(a)可知，在高吸力階段，VG模型對斥水壤土和黏土的擬合值比實(shí)測值偏小，砂土的擬合值與實(shí)測值比較接近；而在低吸力階段，VG模型對3種質(zhì)地斥水土壤的擬合值效果均比較好。從圖2(b)可以看出，BC模型僅在高吸力段對黏土L1的擬合和實(shí)測數(shù)據(jù)存在一定的偏差，對其余土壤擬合效果均比較好。由圖2(c)可以看出，斥水砂土的DP模型擬合曲線與實(shí)測值差異很?。辉诘臀﹄A段，3種質(zhì)地土壤的DP模型擬合值與實(shí)測值非常接近；但在高吸力階段，黏土和壤土的DP模型擬合值比實(shí)測值稍微偏小。而由圖2(d)可知，斥水砂土的LND模型擬合值與實(shí)測值的差異也很小，進(jìn)一步說明LND模型能夠較好的適用斥水砂土。

表2 不同類型斥水土壤各模型擬合統(tǒng)計特征值Tab.2 Statistic eigenvalues of fitting models under different repellent soils

圖2 不同斥水土壤水分特征曲線模型擬合效果Fig.2 Comparison of models fitting effects on soil water characteristics curves under different repellent soils

3 結(jié) 語

(1)斥水土壤的土壤水分特征水曲線，在低吸力階段隨著吸力增加，土壤含水率下降速度較快，土-水曲線呈平緩狀；在高吸力階段，土壤含水率下降速度較慢，土-水曲線呈陡直狀，這與親水土壤變化規(guī)律基本一致。

(2)對于同一種類型的土壤，在同一吸力條件下，土壤斥水程度等級越高，含水率越小，保水能力越低。

(3)對不同斥水土壤的土壤水分特征曲線模型計算值與實(shí)測值進(jìn)行回歸，由擬合特征值分析得出：斥水黏土和壤土的VG、BC和DP模型擬合值和實(shí)測值均呈極顯著相關(guān)關(guān)系，而LND模型則不適用斥水黏土和壤土；對于斥水砂土，VG、BC、LND和DP模型擬合值和實(shí)測值均呈極顯著相關(guān)關(guān)系，具有較好的適用性。

(4)從各模型擬合值與實(shí)測值的模擬效果來看，對于斥水砂土，VG和BC模型擬合的土-水曲線存在拐點(diǎn)，而LND和DP模型擬合的效果相對更好，因此對于斥水砂土推薦使用LND和DP模型；對于黏土和壤土，VG、BC和DP模型擬合效果差異不大，均可使用。

□

[1] 陳俊英，張智韜，汪志農(nóng)，等. 土壤斥水性影響因素及改良措施的研究進(jìn)展[J]. 農(nóng)業(yè)機(jī)械學(xué)報，2010，41(5)：84-89.

[2] 雷志棟，楊詩秀，謝傳森. 土壤水動力學(xué)[M]. 北京：清華大學(xué)出版社，1998.

[3] 邢旭光，趙文剛，馬孝義，等. 土壤水分特征曲線測定過程中土壤收縮特性研究[J]. 水利學(xué)報，2015，46(10)：1 181-1 188.

[4] 鄭 健，王 燕，蔡煥杰，等. 植物混摻土壤水分特征曲線及擬合模型分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報，2014，45(5)：107-112.

[5] 栗現(xiàn)文，周金龍，靳孟貴，等. 高礦化度土壤水分特征曲線及擬合模型適宜性[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(13)：135-141.

[6] 王志超，李仙岳，史海濱，等. 含殘膜土壤水分特征曲線模型構(gòu)建[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(14)：103-109.

[7] van Genuchten M Th. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980,44(44):892-898.

[8] Brooks R H, Corey A T. Hydraulic properties of porous media[R]. Hydrology Paper 3. Fort Collins: Colorado State University, 1964.

[9] Durner W. Hydraulic conductivity estimation for soils with heterogeneous pore structure[J]. Water Resources Research, 1994,30(2):211-223.

[10] Kosugi K. Lognormal distribution model for unsaturated soil hydraulic properties[J]. Water Resources Research, 1996,32(9):2 697-2 703.

[11] 張季如，胡 泳，余紅玲，等. 黏性土粒徑分布的多重分形特性及土-水特征曲線的預(yù)測研究[J]. 水利學(xué)報，2015，46(6)：650-657.

[12] Batuters T W J, Dicalo D A T, Steenhuis S, et al. Preferential flow in water repellent sands[J]. Soil Science Society of America Journal, 1998,62(5):1 185-1 190.

[13] Arye G, Tarchitzky J, Chen Y. Short-term reestablishment of water repellency after witting: Effect on Capillary Pressure-Saturation Relationship[J]. Soil Science Society of America Journal, 2007,71(3):692-702.

[14] 劉世賓，李 毅，Robert Horton. 斥水土壤的水力參數(shù)及水平吸滲規(guī)律[J]. 排灌機(jī)械工程學(xué)報，2013，31(11)：1 000-1 004.

[15] Diamantopoulos E, Durner W, Reszkowska A, et al. Effect of soil water repellency on soil hydraulic properties estimated under dynamic conditions[J]. Journal of Hydrology, 2013,486(4):175-186.

[16] Ritsema C J, Dekker L W. Preferential flow in water repellent sandy soils: principles and modeling implications[J]. Journal of Hydrology, 2000,s231-232(6):308-319.