楊 松 黃劍峰 羅茂泉 馬澤慧 王 磊 吳玉琴

(云南農(nóng)業(yè)大學(xué)水利學(xué)院， 昆明 650201)

斥水性砂土水-氣形態(tài)及其對(duì)斥水-親水轉(zhuǎn)化的影響分析

楊 松 黃劍峰 羅茂泉 馬澤慧 王 磊 吳玉琴

(云南農(nóng)業(yè)大學(xué)水利學(xué)院， 昆明 650201)

處于地球表面的土壤，尤其是耕作層土壤，受外界因素影響其接觸角可能會(huì)發(fā)生改變。通過試驗(yàn)對(duì)不同接觸角的壤質(zhì)砂土中的孔隙水-氣形態(tài)分布狀況進(jìn)行研究，試驗(yàn)結(jié)果表明：砂土顆粒與孔隙水間的接觸角增大會(huì)使表層土壤中水封閉土層厚度減小，但接觸角增大到一定值后，水封閉層的厚度不再發(fā)生變化。隨著接觸角的增大，氣封閉層的厚度不斷減小。與連續(xù)固體表面不同，砂土顆粒的接觸角小于90°甚至降低至36°也會(huì)出現(xiàn)明顯的斥水現(xiàn)象，但隨著砂土中飽和度的增大，斥水現(xiàn)象會(huì)消失，砂土斥水與親水轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度與水封閉向雙開敞轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度基本一致，因此，砂土親水與斥水轉(zhuǎn)化時(shí)對(duì)應(yīng)的臨界含水率與孔隙水氣分布形態(tài)密切相關(guān)，通過理想模型對(duì)兩種不同水-氣形態(tài)下土壤的基質(zhì)吸力變化分析可以發(fā)現(xiàn)，產(chǎn)生這種現(xiàn)象是因?yàn)楫?dāng)土壤由水封閉變?yōu)殡p開敞時(shí)，水-氣交界面在液體側(cè)的曲率中心消失(接觸角小于90°)，氣-液界面引起的基質(zhì)吸力恒為“正”。

孔隙水-氣形態(tài)； 土壤斥水性； 接觸角； 顆粒模型； 孔隙模型

引言

隨著深度的增加，地表土壤層中的水分不斷增多，呈現(xiàn)出不同的形態(tài)。根據(jù)非飽和土壤中的水-氣形態(tài)可以分為3類[1]：水封閉非飽和土、雙開敞非飽和土和氣封閉非飽和土。當(dāng)土壤處于氣封閉狀態(tài)時(shí)，水-氣界面沒有與固體表面相交，此時(shí)氣相為封閉氣泡，只有液相連續(xù)，通常將其視為飽和土壤[2]；隨著土壤含水率不斷降低，氣相和液相分別連續(xù)，此時(shí)土壤為雙開敞狀態(tài)；如果土壤含水率進(jìn)一步降低，孔隙水被氣相隔離封閉，土壤進(jìn)入水封閉狀態(tài)。

斥水性土壤是指水分不能或者很難濕潤(rùn)土顆粒表面的土壤，即：土顆粒與孔隙水間的接觸角較大。處于地表的耕作層土壤，受外界因素影響其接觸角可能會(huì)發(fā)生改變，如：表層土壤含有機(jī)物、森林火災(zāi)、再生水灌溉等[3-6]。當(dāng)土顆粒與孔隙水的接觸角增大，土壤可能變?yōu)槌馑?。斥水性土壤由于降低了土壤的滲透性從而對(duì)耕作層土壤的土水特性有重要影響[7]。斥水性土壤增加了地表徑流，加快土壤侵蝕，還可能降低種子的發(fā)芽率和阻礙作物的生長(zhǎng)[8]。影響土壤斥水性的因素很多，含水率是其中一個(gè)重要因素[9]。當(dāng)含水率不斷增大，土壤可能由斥水變?yōu)橛H水，此時(shí)對(duì)應(yīng)的含水率為臨界含水率。其機(jī)理被解釋為：土壤中水分增加溶解了土壤中造成斥水的有機(jī)物成分[10]。然而，對(duì)于斥水成分穩(wěn)定的土壤也會(huì)出現(xiàn)臨界含水率的現(xiàn)象，其機(jī)理目前還不清楚。

通過理想顆粒模型或孔隙模型可以從微觀角度對(duì)土壤中的三相界面特性及毛細(xì)力展開研究，從而對(duì)土壤的土水特性有更深刻的認(rèn)識(shí)。欒茂田等[11]以理想球體顆粒模型為基礎(chǔ)，考慮張力吸力的影響，在雙開敞基礎(chǔ)上進(jìn)一步把非飽和土劃分為搭接雙開敞和不搭接雙開敞兩種狀態(tài)。張昭等[12-13]利用不等直徑的球體顆粒模型對(duì)非飽和土的水力特性及毛細(xì)力的演化規(guī)律進(jìn)行了深入探討。一般情況下，球體顆粒模型更適用于比表面積較小的壤質(zhì)砂土，采用片狀顆粒模型模擬黏土更為恰當(dāng)[14]。另外，還可以用理想孔隙模型對(duì)土壤的土水特性展開研究，很多土水特征曲線模型都是在圓柱形孔隙模型基礎(chǔ)上推導(dǎo)的[15-16]，有學(xué)者則通過建立正弦型孔隙模型對(duì)土壤中的水分遷移特性展開研究[17-18]。用于描述兩類模型的方程組都涉及到接觸角，因此兩類模型都可以反映土壤的斥水性，然而，兩類模型所對(duì)應(yīng)的水-氣形態(tài)是不同的。

本文對(duì)具有穩(wěn)定斥水性的壤質(zhì)砂土中孔隙水-氣形態(tài)展開研究，在試驗(yàn)基礎(chǔ)上分析斥水性砂土臨界含水率所對(duì)應(yīng)的孔隙水-氣形態(tài)，并據(jù)此分析理想顆粒模型和孔隙模型的適用性，最后通過理想模型對(duì)砂土斥水-親水隨含水率變化的原因進(jìn)行解釋。為進(jìn)一步認(rèn)清土壤斥水-親水的轉(zhuǎn)化機(jī)理提出一種思路。

1 試驗(yàn)材料與方法

1.1 試驗(yàn)材料

試驗(yàn)用土取自云南省昆明市團(tuán)結(jié)鄉(xiāng)某風(fēng)化土坡，土壤的基本物理性質(zhì)如表1所示，土樣經(jīng)過風(fēng)干碾碎后過2 mm篩，采用混合法改變砂土的表觀接觸角，在砂土中加入不同量的斥水劑(十八胺：C18H39N，白色蠟狀粉末)并混合均勻(由于十八胺不溶于水，所以配出的土樣具有穩(wěn)定的斥水性)，準(zhǔn)備5種風(fēng)干砂土試樣，十八胺質(zhì)量分?jǐn)?shù)分別為0、0.05%、0.1%、0.2%、0.3%，分別編號(hào)為S1、S2、S3、S4、S5。

表1 土樣的物理性質(zhì)Tab.1 Physical properties of test samples

1.2 接觸角測(cè)量試驗(yàn)

采用躺滴法對(duì)土顆粒接觸角進(jìn)行測(cè)量。待測(cè)土樣過0.5 mm篩，取長(zhǎng)度為4 cm左右的雙面膠帶貼到載玻片上，將篩好的土顆粒均勻覆蓋載玻片，并用200 g砝碼壓住載玻片，持續(xù)時(shí)間為2 min，移開砝碼后反復(fù)輕輕磕碰載玻片以去除多余的土顆粒，直到膠帶表面覆蓋有一層薄薄的土顆粒，且土顆粒分布均勻無明顯凸起。接觸角測(cè)量試驗(yàn)在JC2000型接觸角測(cè)量?jī)x上進(jìn)行，由于水滴能在S1樣制成的壓片上迅速攤開，因此認(rèn)為此時(shí)的接觸角為0°，其余4種試樣的接觸角測(cè)量結(jié)果如圖1所示，S2、S3、S4、S5接觸角分別為36°、57°、70°、77°。

圖1 S2、S3、S4、S5的接觸角Fig.1 Contact angles of samples S2, S3, S4 and S5

1.3 確定孔隙水-氣狀態(tài)的試驗(yàn)過程

文獻(xiàn)[1]給出了確定砂土孔隙水-氣形態(tài)的試驗(yàn)方法，本文在其基礎(chǔ)上進(jìn)行了適當(dāng)改進(jìn)，具體操作步驟如下：

(1)準(zhǔn)備5根長(zhǎng)度為1 m，內(nèi)徑為5 cm的透明玻璃管。玻璃管底部用濾紙封堵，將5種不同斥水度的砂土分別裝入不同玻璃管內(nèi)，每10 cm裝入314 g砂樣，振搗后使其干密度為1.6 g/cm3左右，直至裝滿玻璃管，如圖2所示。

圖2 試驗(yàn)裝置Fig.2 Test devices

(2)先把裝滿砂樣的玻璃管豎直放入盛滿蒸餾水的容器中，從玻璃管頂部注蒸餾水(由于混有十八胺砂土的接觸角較大，需要施加一定的初始水頭促使水下滲)，待整根玻璃管頂部有水漫出后停止注水，靜止10 min后再繼續(xù)注水，如此反復(fù)3次。把玻璃管從容器中移出，在玻璃管底部用充分吸水的海綿保濕，套上塑料膜防止水分蒸發(fā)，保持玻璃管底部的砂土處于飽和狀態(tài)。

(3)試樣在室內(nèi)靜置12 h，把1 m長(zhǎng)的玻璃管分成50小段，每段2 cm，每段取適當(dāng)土樣放入鋁盒干燥、稱量，測(cè)其含水率并計(jì)算相應(yīng)的飽和度。

1.4 滴水入滲試驗(yàn)

通過滴水入滲試驗(yàn)確定水滴入滲時(shí)間可以評(píng)價(jià)土壤的斥水性，其操作步驟如下：

(1)S2、S3、S4、S5每種土樣配制24個(gè)不同含水率，分別為：3.0%、3.5%、4.0%、4.5%、…、14.0%、14.5%，密封靜置24 h后制成干密度為1.6 g/cm3鋁盒試樣。

(2)如圖3所示，試樣表面保持平整，為延長(zhǎng)液滴的蒸發(fā)時(shí)間，試驗(yàn)時(shí)的相對(duì)濕度保持在90%以上。用注液器把體積為0.02 mL的蒸餾水滴到試樣表面，用秒表記錄液滴完全滲透進(jìn)試樣內(nèi)所用的時(shí)間。試樣表面不同部位進(jìn)行7次滴水入滲試驗(yàn)，將7滴水的入滲時(shí)間去掉最大值和最小值后求出的算術(shù)平均值為每個(gè)試樣的最終測(cè)定結(jié)果，取5 s作為斥水和不斥水的分界點(diǎn)[4]。

圖3 滴水入滲試驗(yàn)Fig.3 Test of water drop penetration

2 試驗(yàn)結(jié)果

2.1 水-氣形態(tài)隨土層深度變化

從地表面到潛水面，土層飽和度(含水率)的變化規(guī)律如圖4a所示[19]，在這個(gè)區(qū)域內(nèi)(包氣帶)，飽和度隨土層深度變化關(guān)系曲線可以分為3個(gè)階段，這3個(gè)階段分別對(duì)應(yīng)非飽和土的3種水-氣形態(tài)，在曲線上部，隨著土層深度增加，含水率基本沒有變化，此時(shí)的非飽和土處于水封閉狀態(tài)。當(dāng)?shù)竭_(dá)一定的土層深度后，土中的水-氣形態(tài)將由水封閉轉(zhuǎn)變?yōu)殡p開敞，在整個(gè)雙開敞區(qū)域內(nèi)，隨著土層深度的增加，土壤的飽和度不斷增大。到達(dá)曲線下部后，土壤進(jìn)入氣封閉狀態(tài)，趨于飽和，隨著深度的增加飽和度基本不變。

試驗(yàn)得到的曲線與示意圖基本一致，但土顆粒的接觸角不同，土中水-氣形態(tài)轉(zhuǎn)化的分界點(diǎn)位置可能會(huì)發(fā)生改變。如圖4b～4f所示，親水土樣S1中，水封閉狀態(tài)轉(zhuǎn)化為雙開敞的分界點(diǎn)A位于深度為28 cm左右的土層中，對(duì)應(yīng)的飽和度為36%，當(dāng)土顆粒的接觸角增大，S2、S3樣的分界點(diǎn)A深度為24 cm和20 cm左右，相應(yīng)的飽和度減小為32%和27%。土顆粒的接觸角進(jìn)一步增大，S4、S5樣中，水-氣形態(tài)轉(zhuǎn)化的分界點(diǎn)A位置則維持在深度為20 cm的土層中，飽和度也保持在20%左右。接觸角的增大會(huì)使水封閉土層厚度減小，但接觸角增大到一定值后，水封閉土層的厚度就基本保持不變。分界點(diǎn)A處的飽和度隨接觸角的增大而減小，接觸角增大到一定值后，分界點(diǎn)A處的飽和度不受接觸角的影響。

當(dāng)試樣中孔隙水-氣形態(tài)由雙開敞轉(zhuǎn)化為氣封閉時(shí)，5個(gè)試樣在分界點(diǎn)B處的飽和度基本都保持在90%～92%，相應(yīng)的S1、S2、S3、S4、S5試樣轉(zhuǎn)化點(diǎn)所在的土層深度分別為82、84、88、94、96 cm。由此可見，接觸角對(duì)試樣雙開敞狀態(tài)轉(zhuǎn)變?yōu)闅夥忾]狀態(tài)時(shí)的飽和度影響不大，但隨著接觸角的增大，氣封閉層的厚度不斷減小。

圖4 飽和度隨土層深度的變化曲線Fig.4 Changing curves of saturation with soil depth

2.2 滴水入滲時(shí)間隨飽和度變化

當(dāng)水滴滴到親水土樣S1表面時(shí)，無論土樣的飽和度如何變化，水滴都會(huì)迅速入滲，入滲時(shí)間趨近于零。其余土樣(S2～S5)，由于接觸角較大，水滴落到土樣表面后并不會(huì)立刻入滲，其入滲時(shí)間與飽和度的關(guān)系如圖5所示，土壤處于風(fēng)干狀態(tài)時(shí)(Srlt;5%)，水滴在土樣表面直至蒸發(fā)結(jié)束后也不入滲，土壤表現(xiàn)出了劇烈的斥水性。隨著飽和度的增大，4種土樣表面水滴的入滲時(shí)間不斷減小。S2樣在飽和度增加到28%左右，水滴入滲時(shí)間縮短到5 s以內(nèi)；S3、S4、S5樣在飽和度為20%～25%之間時(shí)，水滴入滲時(shí)間縮短到5 s以內(nèi)。對(duì)比前面的孔隙水-氣形態(tài)確定試驗(yàn)可以發(fā)現(xiàn)，土樣斥水與親水轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度與其水封閉向雙開敞轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度基本一致。

圖5 滴水入滲時(shí)間T隨飽和度Sr變化曲線Fig.5 Changing curves of water penetration time with saturation

3 討論

3.1 土壤的斥水性與接觸角

通常情況下認(rèn)為，固體表面的接觸角大于90°為斥水，小于90°為親水。然而，通過滴水入滲試驗(yàn)可以看出，土樣S2、S3、S4、S5的接觸角小于90°仍然出現(xiàn)了明顯的斥水現(xiàn)象，特別是S2，其接觸角僅為36°，斥水現(xiàn)象依然很明顯。其實(shí)，土壤的斥水和親水并不是以90°為分界點(diǎn)，SHIRTCLIFFE等[20]從能量的角度證明了等直徑球體顆粒親水與斥水轉(zhuǎn)化時(shí)對(duì)應(yīng)的接觸角小于90°。LIKOS等[21]和YANG等[22]通過等直徑球體顆粒模型計(jì)算的結(jié)果表明：球體顆粒間可能在接觸角小于40°時(shí)就發(fā)生斥水現(xiàn)象，這與本文試驗(yàn)結(jié)果是一致的。除接觸角外，土顆粒的幾何形態(tài)對(duì)土壤親水和斥水性質(zhì)也有重要影響，在斥水性砂土中加入黏土顆?？梢悦黠@改善土壤的斥水性[23]。

滴水入滲試驗(yàn)結(jié)果表明：隨著含水率的增大，土壤的斥水性會(huì)減小甚至消失。一般情況下，天然土壤隨著深度的增加，含水率是不斷增大的，直到進(jìn)入飽和區(qū)。因此，孔隙水與土顆粒的接觸角主要對(duì)水封閉和雙開敞狀態(tài)下的非飽和土產(chǎn)生影響，特別是處于水封閉狀態(tài)下的土層，可能存在基質(zhì)吸力消失的情況。

3.2 斥水性土壤中的保水性

斥水性土壤帶來負(fù)面影響的同時(shí)，也有其有益的一面。如：土壤斥水性可以增強(qiáng)土壤團(tuán)聚體的穩(wěn)定性，保持土壤中的有機(jī)碳含量[24-25]。很顯然，當(dāng)土壤處于水封閉狀態(tài)時(shí)含水率要低于雙開敞狀態(tài)。通過前面的試驗(yàn)可知，親水性土壤表層水封閉土層的厚度要明顯大于斥水性土壤。這表明，雖然斥水性土壤阻止了地表水的入滲，但同時(shí)也能很好地保持土壤中的原有含水率，減少土壤水分的蒸發(fā)量，這也可以視為斥水性土壤的有益面。

3.3 顆粒模型與孔隙模型

各種理想模型為土壤中基質(zhì)吸力計(jì)算及水分遷移分析提供了方便，然而各模型的適用范圍至今也沒有明確的限定。就模型的幾何形態(tài)而言，球體顆粒模型(圖6a)更適合比表面積較小的砂土或粉土，片狀顆粒模型(圖6b)則更符合黏土顆粒的特性，由于模型中土顆粒由液橋連接，所以兩種模型應(yīng)該對(duì)應(yīng)于水封閉狀態(tài)?？紫赌Ｐ?圖6c、6d)適合于各種類型的土壤，然而其幾何形態(tài)對(duì)應(yīng)于雙開敞狀態(tài)更為恰當(dāng)。

圖6 理想顆粒模型與孔隙模型Fig.6 Ideal particle models and pore models

現(xiàn)有文獻(xiàn)中并沒有嚴(yán)格按照孔隙水-氣形態(tài)來選擇模型，通過理想球體顆粒模型和孔隙模型建立的土水特征曲線都能得到試驗(yàn)的驗(yàn)證。當(dāng)然，前提是假設(shè)土顆粒與孔隙水接觸角很小或者直接假設(shè)為零。當(dāng)接觸角較大時(shí)，通過前面的試驗(yàn)可知：土壤在水封閉狀態(tài)為斥水，基質(zhì)吸力為零，而隨著含水率的增大，土壤進(jìn)入雙開敞狀態(tài)，此時(shí)基質(zhì)吸力又會(huì)出現(xiàn)。當(dāng)土顆粒與孔隙水的接觸角較大時(shí)，需要根據(jù)不同的水-氣形態(tài)選擇相適應(yīng)的模型。

3.4 斥水性土壤的臨界含水率

由土水特性曲線可知，含水率越小土壤的基質(zhì)吸力越大。然而，當(dāng)土壤斥水性增大時(shí)，即使含水率很小，基質(zhì)吸力也可能很小甚至消失。即：土壤中不存在負(fù)的孔隙水壓力。因此，可以通過計(jì)算基質(zhì)吸力來反映“吸引”與“排斥”[26]。本文試驗(yàn)用土為比表面積較小的壤質(zhì)砂土，當(dāng)其處于水封閉狀態(tài)時(shí)用顆粒模型描述更為恰當(dāng)，取土顆粒間的液橋作為研究對(duì)象，如圖7a所示，由Young-Laplace方程可知，基質(zhì)吸力φ為

(1)

式中γ——液體的表面張力

r1、r2——任意一對(duì)相互垂直的平面在液面上所截取曲線的曲率半徑

當(dāng)曲率中心在液體側(cè)時(shí)r1、r2為“負(fù)”，反之為“正”。設(shè)φ為“正”時(shí)表示吸引、親水，φ為“負(fù)”時(shí)表示排斥、斥水。接觸角越小，則r1越小，φ越大，土壤的親水性就越好。當(dāng)接觸角增大后，r1迅速增大，而r2仍然維持在較小值，此時(shí)φ為負(fù)，土壤是斥水的，且接觸角越大斥水性越強(qiáng)，但基質(zhì)吸力由正變負(fù)所對(duì)應(yīng)的接觸角顯然小于90°。

圖7 不同水-氣狀態(tài)下的孔隙水Fig.7 Pore water with different pore water-air configurations

土壤中的含水率增大，土中水-氣形態(tài)將由水封閉變?yōu)殡p開敞，此時(shí)應(yīng)當(dāng)選用毛細(xì)管模型，雙開敞中的孔隙水形態(tài)如圖7b所示?；|(zhì)吸力φ為

(2)

式中r3——毛細(xì)管半徑

此時(shí)水-氣交界面嚴(yán)格遵循表面力學(xué)中的定義，接觸角小于90°親水，接觸角大于90°斥水，因?yàn)楫?dāng)接觸角小于90°時(shí)，彎液面的任意曲率中心都是在空氣一側(cè)，φ為“正”。當(dāng)接觸角大于90°，彎液面的任意曲率中心都是在液體一側(cè)，因此φ為“負(fù)”。

由此可見，當(dāng)斥水性土壤由水封閉變?yōu)殡p開敞時(shí)，土壤隨之由斥水變?yōu)橛H水，這是因?yàn)樵诖诉^程中水-氣交界面在液體側(cè)的曲率中心消失(接觸角小于90°)，氣-液界面引起的基質(zhì)吸力恒為“正”。當(dāng)然，這種解釋只適用于固-液接觸角不變的情況。

4 結(jié)論

(1)非飽和砂土中土顆粒與孔隙水的接觸角增大會(huì)使包氣帶中水封閉土層厚度減小，但接觸角增大到一定值后，水封閉土層的厚度就基本保持不變。接觸角對(duì)試樣雙開敞狀態(tài)轉(zhuǎn)變?yōu)闅夥忾]狀態(tài)時(shí)的飽和度影響不大，但隨著接觸角的增大，氣封閉層的厚度不斷減小。

(2)非飽和砂土中接觸角小于90°甚至降低至36°也會(huì)出現(xiàn)明顯的斥水現(xiàn)象，但隨著砂土中飽和度的增大，砂土的斥水性消失，斥水與親水轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度與水封閉向雙開敞轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度基本相同。

(3)采用理想模型對(duì)土壤的微觀特性進(jìn)行研究時(shí)，如果非飽和砂土中接觸角較大，需要根據(jù)土壤中不同的水-氣形態(tài)選擇相應(yīng)的理想模型。

(4)土壤由斥水向親水轉(zhuǎn)化時(shí)對(duì)應(yīng)的飽和度與水封閉向雙開敞對(duì)應(yīng)的飽和度基本一致，通過理想模型對(duì)兩種不同狀態(tài)下土壤的基質(zhì)吸力變化分析可以發(fā)現(xiàn)，產(chǎn)生這種現(xiàn)象是因?yàn)楫?dāng)土壤由水封閉變?yōu)殡p開敞時(shí)，水-氣交界面在液體側(cè)的曲率中心消失(接觸角小于90°)，氣-液界面引起的基質(zhì)吸力恒為“正”。

1 俞培基,陳愈炯.非飽和土的水-氣形態(tài)及其與力學(xué)性質(zhì)的關(guān)系[J].水利學(xué)報(bào),1965(1):16-23.

YU Peiji, CHEN Yujiong. Water-air patterns and mechanics properties of unsaturated soil[J]. Journal of Hydraulic Engineering, 1965(1):16-23.

2 殷宗澤. 土工原理[M]. 北京：中國(guó)水利水電出版社,2007.

3 CZACHOR H, DOERR S H, LICHNER L. Water retention of repellent and subcritical repellent soils: new insights from model and experimental investigations[J]. Journal of Hydrology, 2010, 380(1):104-111.

4 巨娟麗, 李毅, 宋紅陽,等. CaCl2溶液灌溉對(duì)土壤水鹽與斥水性分布的影響[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2014, 45(10):159-166. http:∥www.j-csam.org/ch/reader/view_abstract.aspx?file_no=20141025amp;flag=1amp;journal_id=jcsam. DOI：10.6041/j.issn.1000-1298. 2014.10.025.

JU Juanli, LI Yi, SONG Hongyang, et al. Effects of CaCl2solution irrigation on distributions of soil water, salt and water repellency[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(10):159-166.(in Chinese)

5 WOCHE S K, GOEBEL M O, KIKHAM M B,et al. Contact angle of soils as affected by depty, texture, and land management[J]. European Journal of Soil Science, 2005, 56(2):239-251.

6 REGALADO C M, RITTER A. Characterizing water dependent soil repellency with minimal parameter requirement[J]. Soil Science Society of America,2005,69(6)：1955-1966.

7 VOGELMANN E S, REICHERT J M, PREVEDELL J, et al. Threshold water content beyond which hydrophobic soils become hydrophilic: the role of soil texture and organic matter content[J]. Geoderma, 2013,209-210:177-187.

8 BLACKWELL P S. Management of water repellency in Australia, and risks associated with preferential flow, pesticide concentration and leaching[J]. Journal of Hydrology, 2000, 231-232(231):384-395.

9 陳俊英, 吳普特, 張智韜,等.土壤斥水性對(duì)含水率的響應(yīng)模型研究[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2012,43(1):63-67. http:∥www.j-csam.org/ch/reader/view_abstract.aspx?file_no=20120113amp;flag=1amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2012.01.013.

CHEN Junying, WU Pute, ZHANG Zhitao, et al. Response models for soil water repellency and soil moisture[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(1):63-67.(in Chinese)

10 TAUMER K, STOFFREGEN H, WESSOLET G. Determination of repellency distribution using soil organic matter and water content[J]. Geoderma, 2005, 125(1-2):107-115.

11 欒茂田, 李順群, 楊慶. 非飽和土的基質(zhì)吸力和張力吸力[J]. 巖土工程學(xué)報(bào), 2006, 28(7):863-868.

LUAN Maotian, LI Shunqun, YANG Qing. Matric suction and tension suction of unsaturated soils[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(7):863-868.(in Chinese)

12 張昭，劉奉銀，張國(guó)平，等. 不等徑濕顆粒與液橋相互作用的微觀水力特性[J]. 水利學(xué)報(bào)，2013，44(7)：810-817.

ZHANG Zhao，LIU Fengyin，ZHANG Guoping，et al. Microscopic hydraulic behavior from the interactions between uneven-sized wet particles and liquid bridge[J]. Journal of Hydraulic Engineering，2013，44(7)：810-817. (in Chinese)

13 張昭,劉奉銀,齊吉琳，等. 粗顆粒間液橋毛細(xì)力演化規(guī)律的動(dòng)態(tài)計(jì)算方法[J]. 巖土力學(xué), 2016,37(8): 2263-2270.

ZHANG Zhao，LIU Fengyin，QI Jilin, et al. A dynamic calculation method for evolution law of capillarity forces of liquid bridge between coarse particles[J]. Rock and Soil Mechanics, 2016,37(8): 2263-2270. (in Chinese)

14 楊松, 吳珺華, 董紅艷,等. 砂土和黏土的顆粒差異對(duì)土壤斥水性的影響[J]. 土壤學(xué)報(bào), 2016, 53(2):421-426.

YANG Song, WU Junhua, DONG Hongyan, et al. Effects of particle differences between sands and clay on soil water repellency[J]. Acta Pedologica Sinica, 2016,53(2): 421-426. (in Chinese)

15 FREDLUND D G, XING A. Equations for the soil-water characteristic curve[J]. Canadian Geotechnical Journal, 1994, 31(4): 521-532.

16 RUSSELL A R, BUZZI O. A fractal basis for soil-water characteristics curves with hydraulic hysteresis[J]. Géotechnique, 2012, 62(3):269-274.

17 CZACHOR H. Modelling the effect of pore structure and wetting angles on capillary rise in soils having different wettabilities[J]. Journal of Hydrology, 2006, 328(328):604-613.

18 GUARRACINO L, ROTTING T, CARRERA J. A fractal model to describe the evolution of multiphase flow properties during mineral dissolution[J]. Advances in Water Resources,2013, 67(4):78-86.

19 許兆義. 包氣帶的水分分布特征及給水度的探討[J]. 長(zhǎng)春地質(zhì)學(xué)院學(xué)報(bào), 1986(3): 89-94.

XU Zhaoyi. Feature of the moisture disturbution in aeration zone and specific yield[J]. Journal of Changchun University of Earth Sciences, 1986(3): 89-94. (in Chinese)

20 SHIRTCLIFFE N J，MCHALE G，NEWTON M I，et al. Critical conditions for the wetting of soils[J]. Applied Physics Letters，2006，89(9)：094101.

21 LIKOS W J, LU N. Hysteresis of capillary stress in unsaturated granular soil[J]. Journal of Engineering Mechanics, 2004, 130(6):646-655.

22 YANG S, LU T H. Study of soil-water characteristic curve using microscopic spherical particle model[J]. Pedosphere, 2012, 22(1):103-111.

23 LEELAMANIE D A L, KARUBE J, YOSHIDA A. Clay effects on the contact angle and water drop penetration time of model soils[J]. Soil Science amp; Plant Nutrition, 2010, 56(3):371-375.

24 MATAIX-SOLERA J, DOERR S H. Hydrophobicity and aggregate stability in calcareous topsoils from fire-affected pine forests in southeastern Spain[J]. Geoderma, 2004, 118(1-2):77-88.

25 SPACCINI R, PICCOLO A, CONTE P, et al. Increased soil organic carbon sequestration through hydrophobic protection by humic substances.[J]. Soil Biology amp; Biochemistry, 2003, 34(12):1839-1851.

26 LIKOS W J, LU N. Hysteresis of capillary stress in unsaturated granular soil[J]. Journal of Engineering Mechanics, ASCE, 2004, 130(6):646-655.

PoreWater-AirConfigurationsinWaterRepellentSandySoilandItsEffectsonTransformationofHydrophilicitytoHydrophobicity

YANG Song HUANG Jianfeng LUO Maoquan MA Zehui WANG Lei WU Yuqin

(CollegeofWaterConservancy,YunnanAgriculturalUniversity,Kunming650201,China)

The earth surface soil, especially cultivated soil, contact angles of soil may be changed by external factors. Hydrophilic soil will translate into hydrophobic soil because of the bigger contact angle. Repellent soil will bring a series of agricultural and environmental problems. To develop the study to the pore water-air shape distribution in different contact angle soils, some tests were carried out. The experimental results showed that the increase of contact angles between soil particles and pore water can make the thickness of the closed-water soil layer reduced. When the contact angle was increased to a certain value, the closed-water soil layer thickness would remain unchanged. With the increase of the contact angle, the thickness of closed-air layer was reduced. It was different from the continuous solid surface, the contact angel of soil granule was less than 90°, even reduced to 36°, it would also appear the water-repellent phenomenon obviously, but with the increase of the saturation degree in the sand, the water-repellent phenomenon would disappear. The corresponding saturation of water repellent and hydrophilic conversion was basically consistent with the corresponding saturation of closed-water system and bi-opened system transformation. Therefore, the critical water content of hydrophilicity to hydrophobicity was closely related to pore water-air configurations. The matric suction in different pore water-air configurations was analyzed with ideal models, the causes of this phenomenon was the disappearance of water-air surface centers of curvature in water side (contact angle was less than 90°), the matric suction caused by water-air surface was greater than zero. The results suggested that it was not enough to focus on chemical property of soil water repellency, the mechanics of interface was also an important aspect of repellent soil study.

pore water-air configurations; soil water repellency; contact angle; particle models; pore models

10.6041/j.issn.1000-1298.2017.11.030

S152.7

A

1000-1298(2017)11-0247-06

2017-02-10

2017-05-19

國(guó)家自然科學(xué)基金項(xiàng)目(41662021、51209182)

楊松(1982—)，男，副教授，博士，主要從事農(nóng)業(yè)水土工程研究，E-mail： yscliff007@126.com