入滲水頭對鹽堿土水分運移影響數值模擬研究

羅 朋

（山西省水利水電科學研究院，山西 太原 030002）

為緩解區域水危機，尤其是在淡水資源缺乏而咸水資源較豐富的鹽堿地分布區，充分利用區域內的水土資源，實現農業可持續發展，采用先進的灌水方法是未來農業發展的趨勢。鹽漬土是一種世界性的低產土壤，由于我國水資源短缺和土地資源的日益減少和退化，實行節水灌溉，提高土壤生產能力，改良和利用鹽漬化土壤已成為提高農業生產、合理利用水資源和土地資源的一個重要問題。近年來，國內外許多學者對鹽堿土的灌溉利用及改良方面做了大量的研究工作，特別是在鹽堿土壤中實施節水灌溉對土壤水鹽運動和作物生長的影響是許多學者十分關注的熱點問題，這就需要了解不同灌溉條件下的土壤水鹽運動規律。在不同節水灌溉方式對土壤水鹽運移的研究方面，國內外學者在噴灌、重力滴灌及地下滲灌做了一些研究工作[1-3]，關于不同水頭對土壤水鹽運動的影響，國內外學者大多研究了入滲水頭對土壤水分運動的影響[4-5]，對在鹽堿土壤中不同入滲水頭對土壤水分運動的影響的研究還較少。實際上，研究不同水頭對土壤水分變化的影響，對鹽堿土壤的水鹽運移特征有重要作用，它是鹽堿土實施重力滴灌和地下滲灌技術的一個重要理論依據。本文針對鹽堿土在不同入滲水頭條件下土壤水分運移特征進行分析，通過室內土柱試驗來揭示不同入滲水頭下鹽堿土垂直入滲對土壤中水分運移的影響進行模擬，為在節水灌溉條件下鹽堿地區的水分管理、鹽堿土壤的改良提供可靠的理論依據和現實指導。

1 試驗材料與方法

1.1 供試土壤

土樣取自甘肅省民勤縣雙茨鄉裕民村，取土深度40 cm。土樣的選取共分5個采樣點兩個地方，取回后土壤經過摻拌過2 mm篩子，使土樣均勻并經自然風干。采用烘干法測定風干土壤初始含水率為2.51%，采用重量法測得土壤初始含鹽量為1.634%。

1.2 試驗處理和方法

為研究鹽土在垂直入滲條件下，入滲水頭和初始容重對鹽土水分含量變化的影響，進行了一維垂直入滲試驗。試驗裝置包括試驗土柱、供水設備及其他測量裝置。入滲水采用無二氧化碳的蒸餾水。土柱試驗土壤容重分別為 1.35 g/cm3，1.45 g/cm3和 1.50 g/cm3，進行了在5個不同入滲水頭（2.5 cm，5 cm，10 cm，15 cm，20 cm）下的垂直入滲試驗。

2 結果與分析

土壤入滲特性是反映土壤入滲過程中濕潤鋒隨時間的變化規律以及土壤累積入滲量和入滲速率的變化特征。在積水入滲條件下，土壤中水分的運動過程受表層土壤與下層土壤的水勢梯度（基質勢梯度）及重力作用的影響，使水在土壤孔隙中做不穩定運動。

2.1 濕潤鋒與時間的擬合關系

對試驗處理下濕潤鋒的推進距離與入滲時間的實測數據關系進行分析，表明濕潤鋒的推進距離和入滲時間關系可用冪函數來描述：

式中：l——濕潤鋒推進距離，cm；

t——推進時間，min；

a，b——擬合參數。

對實測數據資料進行擬合[6]，擬合結果見表1。

表1 不同水頭作用下濕潤鋒與時間的關系擬合參數

由表1可以看出，不同容重以及不同水頭條件下的土壤濕潤鋒遷移與入滲時間的關系均呈顯著的冪函數相關。

在濕潤鋒遷移距離和入滲時間的關系式中，系數a為第一分鐘末的濕潤鋒推進距離，它是土壤質地、土壤初始含水量等因素控制條件下反映土壤濕潤鋒推進速度快慢的重要參數之一，也是反映入滲初期施加在土壤表面水勢梯度大小的重要因子。

在同一種土容重條件下，隨著入滲水頭的增加，冪函數中的系數a呈增加趨勢，而指數參數b則呈降低趨勢。為定量描述入滲水頭對參數a和b的影響，繪出不同水頭與參數a和b的關系，如圖1所示。

由圖1可以看出，在三種土壤容重條件下，參數a與入滲水頭的關系呈對數函數相關，同一水頭條件下，容重對參數a的影響較小，說明初期的土壤入滲主要受土壤表層與干土層水勢梯度的影響，水頭的增加提高了水勢梯度，使濕潤鋒前進速率增加，表現為參數a增大。同時參數b與入滲水頭的關系呈線性負相關，但與土壤容重相比，入滲水頭對參數b的影響很小，土壤容重與參數b符合二次拋物線關系。

將參數a與入滲水頭的關系式和土壤容重與參數b的關系代入（1）式，就得到了考慮入滲水頭和土壤容重條件下濕潤鋒的推進距離和入滲時間的擬合關系，如式（2）所示：

2.2 累積入滲量與時間的擬合關系

圖1 γ=1.50入滲水頭、容重和參數a，b的關系

由試驗各處理下的累積入滲量與入滲時間的實測數據關系分析得出，不同土壤容重在不同入滲水頭條件下累積入滲量隨時間的變化曲線趨勢與濕潤鋒的推進距離和入滲時間關系曲線基本相同。分析表明，累積入滲量隨時間關系曲線也可用冪函數來描述，即常用的Kostiakov入滲經驗公式：

式中：I——累積入滲量，cm；

t——推進時間，min；

A，B——擬合參數。

在累積入滲量和時間的變化關系式中，系數A為第一分鐘末的累積入滲量，在土壤質地、土壤初始含水量等因素一定的情況下，它是反映土壤入滲能力大小的重要參數之一。隨著入滲時間的增加，進入土體的水分逐漸增加，被濕潤的土體越來越大，即累積入滲量隨時間逐漸增加。

表2 不同水頭作用下累積入滲量與時間擬合參數值

由表2可以看出，參數A和參數B都隨入滲水頭的增大而增大，且入滲擬合相關系數R均在0.99以上，呈顯著冪函數關系。入滲初期，土體結構沒有發生較大的變化，入滲過程主要靠基質勢梯度和重力梯度影響作用進行增滲。隨著入滲的進行，在入滲水頭壓力和重力的作用下，使土體結構發生變化，表層土體發生消散和分散，并受壓力而致密，從而減小土壤的空隙，使土壤的導水性和透水性相對減弱，也就是土壤水分輸運速度減慢，導致土壤入滲各個參數的整體減慢。

為定量描述入滲水頭對參數A和參數B的影響，繪出不同水頭與A和B的關系，如圖2所示。

由圖2可以看出，三個容重對參數A在入滲水頭的影響很小，這是因為A為入滲初始第一分鐘末的累積入滲量，說明初期的土壤入滲主要受土壤表層與干土層水勢梯度的影響，水頭的增加提高了水勢梯度，使濕潤鋒前進速率增加，表現為參數A增大；同一水頭條件下，隨著容重的增大，參數A也增大，但影響較小。同時可以看出，同一容重下參數B與入滲水頭呈線性正相關，隨著入滲水頭的增大，參數B增大；同一入滲水頭下，隨著容重的增大，參數B增大，且隨入滲水頭的增大各容重間參數B增大的差異變大。

3 結語

通過對鹽堿土垂直一維入滲試驗數據分析，研究了不同入滲水頭下濕潤鋒、累積入滲量與時間的關系，結果表明：一是隨著入滲水頭的增加，入滲速率逐漸減小，濕潤鋒和累積入滲量均可表示為時間的冪函數，且相關性很好，提出了考慮入滲水頭和土壤容重條件下濕潤鋒的推進距離和入滲時間的關系式；二是不同土壤容重在不同入滲水頭條件下，累積入滲量隨時間的變化也可用冪函數來描述，其中參數A受土壤容重和入滲水頭的影響較小，而參數B隨著入滲水頭的增加呈線性增加趨勢，隨著土壤容重的增加呈降低趨勢。

圖2 γ=1.50入滲水頭、容重和參數A和B的關系

[1] 齊學斌.節水灌溉的環境效應研究現狀及研究重點[J].農業工程學報,2000,16（2）:37-40.

[2] 馬東豪,王全九，來劍斌.膜下滴灌條件下灌水水質和流量對土壤鹽分分布影響的田間試驗研究[J].農業工程學報，2005，21（3）:42-46.

[3] 李毅，王文焰，王全九，等.非充分供水條件下滴灌入滲的水鹽運移特征研究[J].水土保持學報，2003，17（1）:1-4.

[4] 馬娟娟，孫西歡，等.入滲水頭對土壤入滲參數的影響[J].灌溉排水學報，2004，23（5）：53-55.

[5] 王芳，孫西歡，等.不同水頭作用下土壤垂直入滲特性試驗研究[J].山西水利，2005，（1）：75-76.

[6] 呂歲菊，李春光.土壤水—鹽運移規律數值模擬研究綜述[J].農業科學研究，2005，（1）：80-84.