賀蘭山西麓綠洲土壤水分的入滲特性

鄭 策, 盧玉東, 郭建青, 李環(huán)環(huán)

(1.長(zhǎng)安大學(xué) 環(huán)境科學(xué)與工程學(xué)院, 陜西 西安 710054; 2.長(zhǎng)安大學(xué) 旱區(qū)地下水文與生態(tài)效應(yīng)教育部重點(diǎn)實(shí)驗(yàn)室, 陜西 西安 710054)

賀蘭山西麓綠洲土壤水分的入滲特性

鄭 策1,2, 盧玉東1,2, 郭建青1,2, 李環(huán)環(huán)1,2

(1.長(zhǎng)安大學(xué) 環(huán)境科學(xué)與工程學(xué)院, 陜西 西安 710054; 2.長(zhǎng)安大學(xué) 旱區(qū)地下水文與生態(tài)效應(yīng)教育部重點(diǎn)實(shí)驗(yàn)室, 陜西 西安 710054)

[目的] 研究賀蘭山西麓腰壩綠洲土壤水分入滲特性，為當(dāng)?shù)睾侠砝盟Y源和農(nóng)田灌溉提供理論依據(jù)。[方法] 在當(dāng)?shù)剡x取21個(gè)試驗(yàn)點(diǎn)進(jìn)行野外雙環(huán)入滲試驗(yàn)，利用不同入滲公式對(duì)實(shí)測(cè)結(jié)果進(jìn)行擬合，并通過(guò)多元分析及地統(tǒng)計(jì)學(xué)方法分析土壤水分入滲特性。[結(jié)果] Horton公式的擬合結(jié)果相關(guān)系數(shù)高達(dá)0.963，相對(duì)誤差與均方根誤差分別僅為4.09%與0.046 mm/min，擬合精度要優(yōu)于Philip公式與Kostiakov公式，且模型參數(shù)具有較好的物理意義。土壤質(zhì)地對(duì)入滲過(guò)程有顯著性影響，土壤中含砂量越高，入滲速率越大，達(dá)到穩(wěn)定所需時(shí)間越短，而土壤容重與初始含水率則分別對(duì)累積入滲量及初始入滲率有一定影響。穩(wěn)定入滲率最大值約為最小值的12倍，變異系數(shù)高達(dá)0.677，且由東北向西南呈現(xiàn)遞減的趨勢(shì)。[結(jié)論] Horton公式適用于描述腰壩綠洲土壤水分入滲過(guò)程，土壤性質(zhì)差異對(duì)入滲過(guò)程有較大影響，同時(shí)在空間上入滲率也存在較強(qiáng)變異性。

土壤水分; 雙環(huán)試驗(yàn); Horton公式; 入滲規(guī)律; 腰壩綠洲

入滲是指水分從土壤表面進(jìn)入土壤內(nèi)部形成土壤水的過(guò)程，同時(shí)也是降水、地表水、土壤水、地下水相互轉(zhuǎn)化的重要環(huán)節(jié)[1]。土壤水分入滲是一個(gè)復(fù)雜的過(guò)程，受到土壤質(zhì)地、容重、初始含水率及地下水位等多因素的影響，同時(shí)在時(shí)空分布上也存在著較多的不確定性。國(guó)內(nèi)外學(xué)者對(duì)水分入滲特性已進(jìn)行了大量的研究，如劉芝芹等[2]研究了金沙江流域土壤水分入滲特征，結(jié)果表明Kostiakov公式模擬雙環(huán)入滲法測(cè)得的入滲過(guò)程效果最佳；Duan等[3]模擬了不同入滲公式對(duì)草地的適用性，結(jié)果表明修正Kostiakov公式擬合效果最佳，而Philip公式擬合效果較差；宋小帥等[4]研究了典型森林類型的土壤水文效應(yīng)，結(jié)果表明Horton模型擬合結(jié)果最接近實(shí)測(cè)值；肖慶禮等[5]研究了黑河中游綠洲不同質(zhì)地土壤水分入滲情況，結(jié)果表明不同質(zhì)地土壤水分入滲與再分布存在明顯差異；武敏等[6]研究了不同質(zhì)地土壤水分入滲規(guī)律，結(jié)果表明壤砂土達(dá)到穩(wěn)定入滲所需時(shí)間少于砂壤土及粉砂壤土。

雖然以往的研究取得了相當(dāng)多的成果，但大多都是集中在水量相對(duì)充足的灌區(qū)之中，對(duì)干旱區(qū)的關(guān)注則相對(duì)較少。并且受土壤性質(zhì)及周?chē)h(huán)境的影響，土壤入滲規(guī)律往往具有一定的地域性。對(duì)于像腰壩綠洲這類地區(qū)，由于沒(méi)有地表水流，日常生活與農(nóng)業(yè)灌溉用水都要靠開(kāi)采地下水，連年的過(guò)量開(kāi)采導(dǎo)致當(dāng)?shù)氐叵滤徊粩嘞陆?，并因此引發(fā)了一系列水資源問(wèn)題。因此了解當(dāng)?shù)赝寥浪秩霛B過(guò)程，為其合理利用水資源提供指導(dǎo)成為亟待解決的問(wèn)題?；诖四康模疚臄M以野外雙環(huán)入滲試驗(yàn)為基礎(chǔ)，利用不同入滲公式擬合實(shí)測(cè)結(jié)果，確定適用于腰壩綠洲的最優(yōu)入滲公式。并通過(guò)多元分析及地統(tǒng)計(jì)學(xué)方法分析不同土壤參數(shù)對(duì)入滲過(guò)程的影響以及入滲特性的空間變異性，以期為腰壩綠洲合理利用水資源、指導(dǎo)當(dāng)?shù)剞r(nóng)田灌溉提供科學(xué)依據(jù)[7]。

1 材料與方法

1.1 研究區(qū)概況

腰壩綠洲位于38°25′—38°36′N(xiāo)，105°32′—105°39′E,是內(nèi)蒙古阿拉善高原的主要農(nóng)牧區(qū)之一，行政區(qū)劃隸屬內(nèi)蒙古自治區(qū)阿拉善左旗巴潤(rùn)別立鎮(zhèn)，面積約81.2 km2，地勢(shì)總體呈東高西低。試驗(yàn)區(qū)內(nèi)的氣候?qū)儆诘湫偷拇箨懶愿珊禋夂?，夏季酷熱而冬季寒冷，全年降雨稀少，氣候干旱。年平均降雨量?jī)H為198 mm，而平均蒸發(fā)量卻達(dá)到2 394 mm。腰壩綠洲位于賀蘭山西麓山前傾斜洪積平原水文地質(zhì)區(qū)，賦存有多層結(jié)構(gòu)的孔隙水及承壓水，地下水的埋深由東北向西南逐漸變淺[8]。

1.2 試驗(yàn)方法

本次野外試驗(yàn)采用定水頭雙環(huán)入滲方法[9]，雙環(huán)為同心環(huán)，內(nèi)、外環(huán)直徑分別為20，40 cm，環(huán)高40 cm，試驗(yàn)時(shí)將雙環(huán)插入地表下20 cm，并在環(huán)底鋪墊3～5 cm的小礫石作為緩沖層，防止表層的砂土被沖起。入滲時(shí)采用手動(dòng)加水的方法，保持內(nèi)外環(huán)的水頭始終在5 cm處不變。試驗(yàn)剛開(kāi)始時(shí)，每3 min記錄加入的水量，隨著單位時(shí)刻入滲量的減少，逐步調(diào)整為每5，10，20 min記錄。每個(gè)點(diǎn)均觀測(cè)150 min(試驗(yàn)點(diǎn)達(dá)到穩(wěn)定入滲的時(shí)間基本不超過(guò)120 min)，隨后計(jì)算每時(shí)刻的入滲速率。

由于試驗(yàn)區(qū)內(nèi)的地下潛水位由東北向西南不斷變淺，故沿著地下水流方向在試驗(yàn)區(qū)由東北向西南布置2條平行的試驗(yàn)線，共計(jì)21個(gè)試驗(yàn)點(diǎn)。在每個(gè)點(diǎn)進(jìn)行入滲試驗(yàn)之前，在其旁邊就近開(kāi)挖1 m深坑，從地表起每隔20 cm分別用100 cm3環(huán)刀與鋁盒取原狀土樣，每晚將當(dāng)日環(huán)刀所取土樣送往當(dāng)?shù)厮苷具M(jìn)行烘干試驗(yàn)測(cè)定其容重及初始含水量；野外試驗(yàn)結(jié)束后在室內(nèi)利用Bettersize 2000激光粒度分布儀對(duì)鋁盒所取土樣進(jìn)行顆粒分析試驗(yàn)，并對(duì)照國(guó)際制土壤質(zhì)地分類標(biāo)準(zhǔn)命名，試驗(yàn)點(diǎn)的土壤質(zhì)地共分為壤砂土、砂壤土與粉壤土3類；利用GPS定位儀確定各點(diǎn)的經(jīng)緯度。具體試驗(yàn)點(diǎn)參數(shù)如表1所示。

1.3 入滲公式擬合

不同入滲公式具有不同的側(cè)重點(diǎn)，選擇最適用于研究區(qū)的入滲公式是了解當(dāng)?shù)厝霛B過(guò)程的基礎(chǔ)[10]。本文選取在土壤入滲研究中最常用的Horton公式、Philip公式和Kostiakov公式作為研究對(duì)象，利用回歸分析方法分別對(duì)野外試驗(yàn)得到的實(shí)測(cè)數(shù)據(jù)進(jìn)行擬合，并確定各公式的入滲參數(shù)。

Horton公式：

ft=fc+(f0-fc)e-kt

(1)

Philip公式：

(2)

Kostiakov公式：

ft=nt-α

(3)

式中：ft——開(kāi)始入滲后t時(shí)刻的入滲速率(mm/min)；fc——穩(wěn)定后的入滲速率(mm/min)；f0——入滲初始時(shí)刻的入滲速率(mm/min)；t——時(shí)間(min)；k，A，n，α——相關(guān)入滲經(jīng)驗(yàn)參數(shù)。具體擬合方法參考文獻(xiàn)[11]。

表1 試驗(yàn)點(diǎn)土壤特性參數(shù)

2 結(jié)果與討論

2.1 入滲公式的適用性評(píng)價(jià)

分別利用Horton公式、Philip公式與Kostiakov公式對(duì)實(shí)測(cè)入滲數(shù)據(jù)進(jìn)行擬合，結(jié)果如圖1所示。因篇幅限制，3種質(zhì)地分別以砂壤土4號(hào)點(diǎn)、壤砂土5號(hào)點(diǎn)與粉壤土18號(hào)點(diǎn)為例。

由圖1的擬合結(jié)果可以明顯看出，Horton公式擬合效果最好，而Philip公式與Kostiakov公式擬合的結(jié)果與實(shí)測(cè)值相比則有一定差距。

為了具體評(píng)價(jià)入滲公式的適用性，選取相關(guān)系數(shù)(R)、相對(duì)誤差(RE)、均方根誤差(RMSE)作為評(píng)價(jià)指標(biāo)。不同入滲公式擬合結(jié)果如表2所示，其中RMSE通過(guò)下式計(jì)算[12]。

(4)

式中：fti——該試驗(yàn)點(diǎn)第i個(gè)時(shí)間點(diǎn)的入滲公式擬合值；fti′——該試驗(yàn)點(diǎn)第i個(gè)時(shí)間點(diǎn)的雙環(huán)試驗(yàn)實(shí)測(cè)值；n——該試驗(yàn)點(diǎn)入滲過(guò)程的記錄點(diǎn)數(shù)。

圖1 實(shí)測(cè)值與擬合值對(duì)比

統(tǒng)計(jì)值Horton公式RRE/%RMSE/(mm·min-1)Philip公式RRE/%RMSE/(mm·min-1)Kostiakov公式RRE/%RMSE/(mm·min-1)最大值0.9987.820.0970.9939.020.1210.98510.540.138最小值0.9022.110.0210.9042.670.0430.8821.980.032平均值0.9634.090.0460.9427.810.0540.9328.590.061

注：R表示相關(guān)系數(shù)；RE表示相對(duì)誤差； RMSE表示均方根誤差。

由表2結(jié)果可知，Horton公式擬合結(jié)果的相關(guān)系數(shù)R平均值最高，為0.963，Philip公式次之，Kostiakov公式值最小，表明Horton公式擬合效果最好。而Horton公式的相對(duì)誤差RE與均方根誤差RMSE的平均值分別為4.09%與0.046 mm/min，均為3個(gè)公式最小，表明Horton公式擬合誤差最小。

由Philip公式擬合得到穩(wěn)定入滲率fc值，其與實(shí)測(cè)的fc值對(duì)比如表3所示?？梢钥闯鲈跀M合得到的fc中，共有5個(gè)點(diǎn)為負(fù)值，不符合實(shí)際情況。其余為正值的點(diǎn)，平均相對(duì)誤差為15.26%，而其中11號(hào)點(diǎn)與15號(hào)點(diǎn)均超過(guò)了30%，與實(shí)測(cè)值相比有較大的差距。表明Philip公式并不適用于腰壩綠洲。

表3 Philip公式實(shí)測(cè)fc值與擬合fc值比較

而Kostiakov公式由于本身較為簡(jiǎn)單，導(dǎo)致擬合的結(jié)果在首尾兩端往往與實(shí)測(cè)值不符。由圖1也可以看出，在初始時(shí)刻，入滲速率往往偏大，在時(shí)間較長(zhǎng)時(shí)，入滲速率又小于穩(wěn)定入滲率，不符合實(shí)際情況。但在給定的時(shí)間長(zhǎng)度內(nèi)，該公式使用較為簡(jiǎn)便及精確。若需更廣泛的應(yīng)用，還需對(duì)公式作進(jìn)一步的改善。例如Lewis在Kostiakov公式的基礎(chǔ)上，引入了fc，提出了Kostiakov公式的三參數(shù)形式，適用范圍則更廣泛。

綜合來(lái)看，Horton公式擬合精度要優(yōu)于Philip公式與Kostiakov公式，式中的參數(shù)有良好的實(shí)際意義，可被用來(lái)描述當(dāng)?shù)赝寥浪秩霛B過(guò)程。

2.2 土壤參數(shù)對(duì)入滲過(guò)程的影響

土壤水分入滲過(guò)程受到土壤質(zhì)地、容重、初始含水率、地下水位及入滲時(shí)間等多因素綜合影響，在本次試驗(yàn)過(guò)程中，由于地下水位普遍大于10 m，同時(shí)所有試驗(yàn)點(diǎn)采用統(tǒng)一入滲時(shí)間，可忽略地下水位及入滲時(shí)間對(duì)試驗(yàn)的影響，因而主要影響因素為土壤參數(shù)。利用SPSS軟件分析試驗(yàn)數(shù)據(jù)可知，土壤質(zhì)地對(duì)入滲過(guò)程具有顯著性影響，其決定系數(shù)R2高達(dá)0.89，而土壤容重及初始含水率對(duì)入滲過(guò)程的影響則相對(duì)較小，且是基于同一土壤質(zhì)地基礎(chǔ)之上，下面具體討論各個(gè)參數(shù)對(duì)入滲過(guò)程的影響。

2.2.1 土壤質(zhì)地 為了定量研究不同質(zhì)地土壤對(duì)入滲過(guò)程的影響，分別統(tǒng)計(jì)了壤砂土、砂壤土、粉壤土的野外試驗(yàn)實(shí)測(cè)初始入滲率f0，穩(wěn)定入滲率fc及穩(wěn)定入滲時(shí)間t的分布，如表4所示。由于野外所有點(diǎn)觀測(cè)時(shí)間均為150 min，在作出入滲速率圖后，假定若10 min內(nèi)入滲速率相差小于0.01 mm/min，即認(rèn)為達(dá)到穩(wěn)定入滲。

回想自己“被選擇”的經(jīng)歷，梅先生卻有了另一番感悟：“相比較現(xiàn)代社會(huì)的選擇太多，誘惑太多，當(dāng)年的‘被選擇’也未必就是件壞事。被選擇了，卻省去了許多時(shí)間來(lái)做思前想后的考慮，迷茫也少些；被選擇了，也不會(huì)有太多的彎路讓人去走。‘鍥而不舍，金石可鏤’。任何一件事，你只要堅(jiān)持，就是專家。”

由表4中數(shù)據(jù)可以看出，不同質(zhì)地土壤入滲規(guī)律差異較大，以初始入滲率為例，最大的壤砂土入滲速率為2.218 mm/min,而最小的粉壤土僅為0.144 mm/min。總體上初始入滲率與穩(wěn)定入滲率均隨土壤中含砂量的降低而降低，說(shuō)明隨著土壤結(jié)構(gòu)由疏松變密實(shí)，其入滲能力也變差。

表4 不同質(zhì)地土壤入滲參數(shù)

注：f0為初始入滲率；fc為穩(wěn)定入滲率；t為穩(wěn)定入滲時(shí)間。下同。

穩(wěn)定入滲時(shí)間則隨土壤中含砂量的減少而增加，由圖1實(shí)測(cè)值曲線亦可以看出，壤砂土在試驗(yàn)初入滲速率快速下降，降幅較大，平均達(dá)到穩(wěn)定入滲僅需68 min；砂壤土入滲速率降幅變慢，平均達(dá)到穩(wěn)定入滲需79 min；而粉壤土入滲速率變化則較慢，變化幅度較小，平均達(dá)到穩(wěn)定入滲則需101 min。

2.2.2 容重 土壤容重是反應(yīng)土壤緊密程度的重要指標(biāo)，主要受土壤中孔隙分布的影響，而水分入滲的實(shí)質(zhì)即水分在孔隙中流動(dòng)并不斷深入的過(guò)程，因而容重對(duì)水分入滲過(guò)程必然存在影響[13]。按照土壤質(zhì)地分類，得到累積入滲量I與容重之間的關(guān)系如圖2所示。

圖2 累積入滲量與土壤容重關(guān)系

由圖2可以看出，壤砂土容重從1.62 g/cm3增加到1.66 g/cm3，累積入滲量由11.38 cm減小到6.95 cm；砂壤土容重從1.37 g/cm3增加到1.57 g/cm3，累積入滲量由7.58 cm減小到3.65 cm；粉壤土容重從1.38 g/cm3增加到1.44 g/cm3，累積入滲量由3.08 cm減小到1.11 cm。經(jīng)顯著性檢驗(yàn)，結(jié)果表明容重與累積入滲量之間存在顯著的線性相關(guān)關(guān)系(sig壤砂土=0.004<0.05,sig砂壤土=0.007<0.05,sig粉壤土=0.004<0.05)，并且3種情況下相關(guān)系數(shù)R均大于0.85，表明在同一質(zhì)地土壤中，累積入滲量隨容重增大呈減小的趨勢(shì)。分析此現(xiàn)象原因，在同一種質(zhì)地的土壤中，容重較小的孔隙較多，因此入滲能力較強(qiáng)，累積入滲量也較大。

2.2.3 初始含水率 土壤初始含水率是影響水分入滲過(guò)程的重要因子，對(duì)濕潤(rùn)鋒的運(yùn)移速率及穩(wěn)定入滲時(shí)間都有較大影響。分別選取本次試驗(yàn)3種不同質(zhì)地土壤中初始含水率最大與最小的樣品，統(tǒng)計(jì)相關(guān)土壤入滲參數(shù)(如表5所示)。

表5 不同初始含水率土壤入滲參數(shù)

注：I為累積入滲量。

由表5中3組對(duì)比結(jié)果可以看出，初始含水率較低的土壤，初始入滲速率遠(yuǎn)大于初始含水率較高土壤，平均增幅高達(dá)36%，表明其受初始含水率影響較大；同時(shí)初始含水率較低的土壤達(dá)到穩(wěn)定入滲所需時(shí)間也相對(duì)較長(zhǎng)，但增幅僅為11%，表明其受初始含水率影響相對(duì)較小。因?yàn)樵谌霛B初期，土壤水分處于非飽和狀態(tài)，土壤含水率越低則土水勢(shì)越低，此時(shí)土壤對(duì)水分子吸引力較大，因而初始含水率越低的土壤初始入滲率則越大，達(dá)到穩(wěn)定入滲所需的時(shí)間也越久。而穩(wěn)定入滲率與累積入滲量則并未有表現(xiàn)出任何規(guī)律，主要是因?yàn)楸敬卧囼?yàn)之前所有試驗(yàn)點(diǎn)均已處于長(zhǎng)期未灌溉狀態(tài)，初始含水率普遍較低且差距不大，并且試驗(yàn)時(shí)并未刻意去控制其他變量，因而初始含水率對(duì)入滲過(guò)程的影響主要表現(xiàn)在影響初始入滲率上，而對(duì)最終穩(wěn)定含水率及累積入滲量則影響不大。

2.3 入滲率的空間變異性分析

腰壩綠洲處于山前傾斜洪積平原的中下部地區(qū)，沉積物由不同時(shí)期洪積扇重疊堆積而成，由山前向西沉積物顆粒逐漸變細(xì)。

根據(jù)顆粒分析試驗(yàn)結(jié)果可以看出，土壤中砂含量由東向西有減少的趨勢(shì)，為了研究試驗(yàn)區(qū)內(nèi)入滲率的空間分布情況，利用地統(tǒng)計(jì)學(xué)方法對(duì)入滲率數(shù)據(jù)進(jìn)行分析，將土壤入滲參數(shù)視為相互獨(dú)立隨機(jī)變量，分別計(jì)算初始入滲率f0與穩(wěn)定入滲率fc樣本的均值、方差、標(biāo)準(zhǔn)差、變異系數(shù)及偏態(tài)系數(shù)，得到相應(yīng)的結(jié)果如表6所示。

表6 入滲率的統(tǒng)計(jì)特征值

由表6中數(shù)據(jù)可以看出初始入滲率與穩(wěn)定入滲率均表現(xiàn)出較大的空間變異性。其最大值分別約為最小值的16與12倍；變異系數(shù)分別為0.547與0.677，屬于中等強(qiáng)度變異；偏態(tài)系數(shù)分別為1.612與4.459，均表現(xiàn)出正偏趨勢(shì)。為了更加直觀描述入滲率在空間上的分布情況，利用Kriging內(nèi)插值法繪制穩(wěn)定入滲率在研究區(qū)內(nèi)的等值線圖(如圖3所示)。

圖3 研究區(qū)穩(wěn)定入滲率fc空間分布

由圖3結(jié)果可以明顯看出，總體上穩(wěn)定入滲率東部較西部高，北部較南部高，即沿試驗(yàn)線布設(shè)方向入滲率呈逐步減小趨勢(shì)，具有較強(qiáng)的空間變異性。中西部個(gè)別點(diǎn)入滲率較大，主要是由于該點(diǎn)處土壤含砂量較高，分析其原因。中部地區(qū)主要是受耕作影響，個(gè)別農(nóng)戶為了耕作需要引入砂土；而西部地區(qū)如20號(hào)點(diǎn)則是由于位于腰壩綠洲最西邊，緊鄰騰格里沙漠，土地沙化現(xiàn)象嚴(yán)重。綜合來(lái)看，入滲率的空間分布主要受洪積扇堆積物影響，呈現(xiàn)出由東北向西南逐步減小的趨勢(shì)，但在局部地區(qū)會(huì)受農(nóng)業(yè)耕作及周邊環(huán)境影響。

3 結(jié) 論

(1) 通過(guò)擬合結(jié)果可以看出，Horton公式擬合結(jié)果的相關(guān)系數(shù)平均值高達(dá)0.963，相對(duì)誤差與均方根誤差分別僅為4.09%與0.046 mm/min，表明Horton公式的擬合精度要好于Philip公式與Kostiakov公式。且Horton公式的參數(shù)有較好的物理意義，符合實(shí)際情況，說(shuō)明Horton公式可被用于描述腰壩綠洲土壤水分入滲過(guò)程。

(2) 通過(guò)多元回歸及地統(tǒng)計(jì)學(xué)分析可以發(fā)現(xiàn)，土壤水分入滲過(guò)程受土壤性質(zhì)差異影響較大，其中受土壤質(zhì)地差異影響尤為顯著；同時(shí)入滲率存在較強(qiáng)空間變異性，穩(wěn)定入滲率最大值為最小值的12倍，變異系數(shù)高達(dá)0.677，空間分布主要受洪積扇堆積物分布影響，呈現(xiàn)出由東北向西南遞減的趨勢(shì)。

[1] 邵明安,王全九,黃明斌.土壤物理學(xué)[M].北京:高等教育出版社,2006:126-128.

[2] 劉芝芹,郎南軍,彭明俊,等.金沙江流域典型森林土壤水分入滲特征試驗(yàn)研究[J].水土保持通報(bào),2014,34(2):43-47.

[3] Duan Runbin, Fedler C B, Borrelli J. Field evaluation of infiltration models in lawn soil[J]. Irrigation Science,2011,29(5):379-389.

[4] 宋小帥,康峰峰,韓海榮,等.遼河源典型森林類型的土壤水文效應(yīng)[J].水土保持通報(bào),2015,35(2):101-105.

[5] 肖慶禮,黃明斌,邵明安,等.黑河中游綠洲不同質(zhì)地土壤水分的入滲與再分布[J].農(nóng)業(yè)工程學(xué)報(bào),2014,30(2):124-131.

[6] 武敏,馮紹元,孫春燕,等.北京市大興區(qū)典型土壤水分入滲規(guī)律田間試驗(yàn)研究[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2009,14(4):98-102.

[7] Franzluebbers A J. Water infiltration and soil structure related to organic matter and its stratification with depth[J]. Soil Tillage Research, 2002,66(2):197-205.

[8] 姜凌,李佩成,胡安炎,等.內(nèi)蒙古阿拉善腰壩綠洲地下水水化學(xué)特征[J].干旱區(qū)資源與環(huán)境,2009,23(11):105-110.

[9] 朱良軍,張光輝,任宗萍.4種土壤入滲測(cè)定方法的比較[J].水土保持通報(bào),2012,32(6):163-168.

[10] 周鑫洋,王仰仁,劉群昌.四個(gè)入滲公式對(duì)農(nóng)田潮土的適應(yīng)性研究[J].灌溉排水學(xué)報(bào),2015,34(8):96-100.

[11] 徐海芳,劉恩民,鄭麗萍.農(nóng)田土壤表層入滲指數(shù)的空間相關(guān)性研究[J].中國(guó)農(nóng)村水利水電,2004(11):11-13.

[12] 聶衛(wèi)波,武世亮,馬孝義,等.農(nóng)田土壤入滲特性研究[J].干旱地區(qū)農(nóng)業(yè)研究,2013,31(4):31-37.

[13] 趙文舉,郁文,范嚴(yán)偉,等.滴灌條件下容重對(duì)壓砂土壤水分入滲規(guī)律的影響[J].水土保持通報(bào),2016,36(1):169-173.

SoilWaterInfiltrationCharacteristicsinOasisonWestSideofHelanMountains

ZHENG Ce1,2, LU Yudong1,2, GUO Jianqing1,2, LI Huanhuan1,2

(1.SchoolofEnvironmentScienceandEngineering,Chang’anUniversity,Xi’an,Shaanxi710054,China; 2.KeyLaboratoryofSubsurfaceHydrologyandEcologicalEffectsinAridRegion,MinistryofEducation,Chang’anUniversity,Xi’an,Shaanxi710054,China)

[Objective] The aim of this study was to analyze the characteristics of soil water infiltration in Yaoba Oasis on the west side of Helan Mountains, and to provide theoretical basis for utilizing water resource efficiently and determining a reasonable way for field irrigation. [Methods] The double-ring infiltration experiments from 21 test spots were carried out. The three different infiltration formulas were used to fit the experimental data. And the soil water infiltration characteristics were analyzed using multivariate analysis and geostatistics methods. [Results] The correlation coefficient of fitting results of Horton formula was as high as 0.963, and the relative error and root mean square error were only 4.09% and 0.046 mm/min, respectively, indicating that the fitting accuracy of Horton formula was better than those of Philip formula and Kostiakov formula. Soil texture had a significant impact on infiltration process, great proportion of sand components in the soil had a faster infiltration rate and had relatively shorter time required to achieve stable infiltration. While the soil bulk density and initial water content had a certain influence on cumulative infiltration volume and initial infiltration rate, respectively. The maximum of stable infiltration rate was about 12 times than the minimum one and the coefficient of variation was up to 0.677. The spatial distribution plot showed a decreasing trend of infiltration rate from the northeast to the southwest of the oasis. [Conclusion] Horton formula could be used to describe the process of soil water infiltration in Yaoba Oasis. The difference of soil properties had a great influence on infiltration process, and the infiltration rate presented a strong spatial variability.

soilwater;double-ringtest;Hortonformula;infiltrationlaw;YaobaOasis

A

1000-288X(2017)05-0146-06

S152.7

文獻(xiàn)參數(shù)： 鄭策, 盧玉東, 郭建青, 等.賀蘭山西麓綠洲土壤水分的入滲特性[J].水土保持通報(bào)，2017,37(5)：146-151.

10.13961/j.cnki.stbctb.2017.05.025； Zheng Ce, Lu Yudong, Guo Jianqing, et al. Soil water infiltration characteristics in oasis on west side of Helan Mountains[J]. Bulletin of Soil and Water Conservation, 2017,37(5)：146-151.DOI:10.13961/j.cnki.stbctb.2017.05.025

2017-03-25

2017-04-19

國(guó)家自然科學(xué)基金面上項(xiàng)目“灌區(qū)包氣帶氮素生物地球化學(xué)動(dòng)態(tài)對(duì)水分的響應(yīng)過(guò)程研究”(41273104); 阿拉善SEE生態(tài)協(xié)會(huì)資助項(xiàng)目(220029150182)

鄭策(1993—),男(漢族),陜西省西安市人,博士研究生,研究方向?yàn)楣喔热霛B。E-mail:2015129041@chd.edu.cn。

盧玉東(1969—),男(漢族),河北省唐山市人,教授,博士生導(dǎo)師,主要從事水資源利用方面的研究。E-mail:luyudongphd@163.com。