微咸水礦化度對重度鹽堿土壤入滲特征的影響

劉小媛，高佩玲，2，張晴雯，王乃江，趙連東

(1.山東理工大學 農業工程與食品科學學院， 山東 淄博 255000； 2.山東理工大學 資源與環境工程學院， 山東 淄博 255000;3.中國農業科學院農業環境與可持續發展研究所， 北京 100086)

我國淡水資源缺乏，時空分布不均，但淺層地下微咸水資源比較豐富。依據微咸水分類標準，我國一般認為含鹽量在2～5 g·L-1范圍內的水資源為微咸水[1]。據統計，中國地下微咸水資源約200億m3,其中可開采水量為130億m3,大部分存在于地下10～100 cm處，宜于開采利用。山東省處于我國黃淮海平原的中部，是我國嚴重缺水省份之一，水資源已經成為制約山東省社會經濟發展的“瓶頸”。山東省礦化度大于2 g·L-1的咸水資源主要分布在黃泛平原區和濱海平原區，總面積18 091 km2,其中濰坊、東營、濱州分布面積較大，科學合理開發利用微咸水資源，對于緩解淡水資源短缺，擴大農業水源，抗旱增產有著極其重要的作用。

近年來，國內外一些學者對微咸水灌溉與土壤水鹽運移及作物生長之間的關系做了相關研究。在微咸水灌溉對土壤理化性質的影響方面，蘇瑩等[2-5]研究表明，在一定礦化度范圍內，土壤入滲能力隨著微咸水礦化度的增加而增大；但超過一定的礦化度，隨著礦化度的增加，帶入土壤中鈉離子數量增加，土壤導水能力下降。任長江等[6]研究了不同初始氯化鈉含量對入滲特性的影響，認為土壤初始含鹽量越高，入滲濕潤鋒、累積入滲量和入滲率越小。畢遠杰等[7]得出入滲水頭對入滲系數、入滲指數、累積入滲量以及濕潤鋒推進距離都有顯著的影響，呈現出顯著性起伏與相對平穩的交替趨勢。針對微咸水灌溉對作物生長以及產量的影響，學者們主要研究了作物對鹽分的抵抗、吸收及其在作物體內的轉移[8-20]。研究表明：微咸水灌溉在帶入土壤水分的同時帶入鹽分，不僅造成鹽漬化的潛在危機，而且可能對植株生長造成直接影響。以上研究主要集中于對微咸水灌溉下非鹽堿土壤水鹽運移特征的研究，對于微咸水灌溉條件下重度鹽堿耕地土壤水鹽運移規律的研究尚未報道。

本文以山東省東營市重度鹽堿土為研究對象，分析微咸水灌溉條件下土壤水鹽運移規律。通過室內試驗，研究淡水及4種不同礦化度(2，3，4，5 g·L-1)微咸水對重度鹽堿土累積入滲量、濕潤鋒運移深度、土壤含水率、入滲歷時、入滲速率的影響；分析微咸水灌溉條件下，土壤水入滲模擬模型的適用性。以上研究將為重度鹽堿耕地微咸水合理高效利用，提供理論和數據支持。

1 材料與方法

1.1 供試土樣的基本性質

供試土壤取自山東省東營市墾利縣的棉花耕地，從表層至60 cm深度每隔20 cm分層取土，土樣經過風干、碾壓、篩分(2 mm篩)、均勻混合后制備成室內試驗土樣。對土壤的基本物理化學性質進行測定，風干土含水率為1.00%，飽和含水率為33.63%，EC5∶1為0.965 ms·cm-1，初始含鹽量為4.099 g·kg-1,按照我國華北平原土壤鹽堿化程度分級標準，供試土壤為重度鹽堿化土壤。利用Mastersizer3000型激光粒度儀測定土壤顆粒組成，顆粒組成見表1，且按照國際制土壤質地分類標準對供試土壤的質地進行劃分。在取擾動土的同時，利用200 mm3環刀分層(0～20，20～40，40～60 cm)取原狀土，在室內測定土壤容重為1.41 g·cm-3和田間持水率為27.85%。

表1 供試土壤顆粒組成及分類 Table 1 Granulometric composition and classification of the experimental soil

1.2 灌溉水量及水質

灌水量根據SL207-98《節水灌溉技術規范》[21]中的有關參數和試驗地的實際情況采用下式計算：

M=H(θmax-θ0)γ

(1)

式中，H為土壤計劃濕潤層深度，由于在地下水位較高的鹽堿化地區，計劃濕潤層深度不宜大于60 cm，故H取60 cm；θmax為土壤計劃濕潤層所允許的最大含水率(占干土重)，一般為田間持水率，取27.85%；θ0為土壤計劃濕潤層初始含水率(占干土重)，取1.00%；γ為土壤容重，取1.41 g·cm-3。由公式(1)計算得到一次實驗的灌溉水量為22.71 cm。

室內試驗所采用的不同礦化度的微咸水是利用化學試劑配制而成，根據試驗水質需要，利用研究區淺層地下水的化學組成進行配水。不同礦化度微咸水所需各種可溶性鹽的量見表2。同時使用蒸餾水作為淡水(0g·L-1)對照。

1.3 試驗設備

試驗設備主要包括試驗土柱和供水設備。試驗土柱采用有機玻璃材料制成，內徑為8cm，高為90cm，土柱側面10cm以下，每隔5cm開一圓形取樣口，其直徑為15mm，便于取土進行土壤含水率分析；利用馬氏瓶供水，其內徑為8cm，高為50cm，供水水頭控制在2cm左右。在試驗前按照當地平均容重分16層(每層5cm厚)均勻填裝。試驗裝置如圖1所示。

表2 不同礦化度水質化學組成Table 2 The chemical composition of different saline water

注：1.輸水軟管;2.試驗土柱;3.馬氏瓶;4.進氣管;5.進氣口;6.桌面;7.取土口

Note：1,waterhose; 2,soilcolumn; 3,markovbottle; 4,airinletvent; 5,airinlet; 6,desk; 7,soilextractionport

圖1 一維垂直積水入滲試驗裝置圖

Fig.1Experimentdevicediagramofonedimensionverticalinfiltrationtest

1.4 試驗過程

將初始含水量為1%的供試土樣按土壤容重1.41g·cm-3分16層(每層5cm厚)均勻裝入土柱，每層裝入質量為358.00g土壤后,利用壓實器壓實土壤，使其達到規定高度5cm,填裝完畢后，在土表放置與土柱截面積相同的濾紙，以防止灌水時對表土的沖刷。利用馬氏瓶供水，進行垂直入滲試驗。

在試驗過程中，按照先密后疏的時間間隔定時觀察、記錄土柱濕潤鋒遷移深度和馬氏瓶水位。當灌水定額入滲結束后，立即從土表至濕潤鋒處每隔5cm沿垂向提取土樣，用烘干法測定土壤含水量。

1.5 土壤入滲過程模型簡介

土壤水分入滲指降水或灌溉水由地表進入土壤的過程，屬于土壤水分運動的一部分，并與地表產流、降雨后土壤水再分配、農田水分最優調控、土壤侵蝕、養分隨水分的遷移、農業面源污染等問題密切相關。

土壤入滲模型能夠表示入滲速率隨時間變化規律，從而對土壤入滲過程作出定性或定量的評價。國內外一些學者對此進行了研究，并提出了許多入滲模型，目前影響較大的有Green-Ampt模型，Kostiakov模型，Kostiakov-Leiws模型,philip模型和Horton模型等[22-26]?；趯崪y數據，利用以下三種模型進行分析。

(1)Kostiakov入滲模型：

I=Ktα

(2)

式中，I為累積入滲量(cm)；t為入滲時間(min)；K,α為入滲參數。

i(t)=Kt-α

(3)

式中，i(t)為t時刻的入滲率(cm·min-1)；K為入滲系數，表示第一時間末的土壤入滲率(cm·min-1)；t為入滲時間(min)。

Kostiakov入滲方程是一個典型的經驗公式，模型參數較容易通過實測資料或試驗獲取，因而在田間灌溉，短歷時入滲研究和實踐中得到普遍應用。

(2)Philip入滲模型：

I=1/2St-0.5+A

(3)

式中，I為入滲率(cm·min-1);S為土壤吸濕率(cm·min-0.5)；t為入滲時間(min)；A為常數(cm·min-1)，通常代表穩定入滲率或穩定下滲強度。

(3) Green-Ampt入滲模型：

i(t)=Ks(Zf+Sf+H)/Zf

(4)

式中，i(t)為t時刻的入滲率(cm·min-1)；Ks為土壤表征飽和導水率(cm·min-1)；H為土壤表面積水深度(cm)；Sf為濕潤鋒面吸力(cm)；Zf為濕潤鋒深度(cm)。

該模型研究的是初始干燥土壤在薄層積水條件下入滲問題。基本假定：入滲時存在明確水平濕潤鋒面，同時具有固定不變的吸力Sf，土壤含水率分布呈現階梯狀，濕潤區為飽和含水量，濕潤鋒前為初始含水率。

本文將采用這些模型，對微咸水灌溉條件下重度鹽堿土入滲過程進行模擬，分析其適用性。

2 結果與分析

2.1 微咸水礦化度對累積入滲量的影響

累積入滲量指在一定時間內通過單位面積的總水量。不同礦化度微咸水灌溉下，累積入滲量隨入滲時間的變化如圖2。5個處理的累積入滲量均隨時間的延長而增加，但增長的速度逐漸變緩，灌水定額入滲結束時入滲歷時大小順序為：5 g·L-1<3 g·L-1<4 g·L-1<2 g·L-1<淡水。在相同入滲歷時內，隨著入滲水礦化度的增加，累積入滲量也隨之增加。造成這種變化規律的原因是微咸水改善了土壤結構，增強了土壤的導水能力，使累積入滲量也隨之增加。但與淡水相比，2，3，4，5 g·L-1微咸水累積入滲量隨入滲時間的變化差異較小，曲線基本重合，主要原因可能是研究土樣是重度鹽堿土壤，離子含量豐富，且為砂質壤土，大孔隙多，入滲能力較強，故入滲水礦化度增加對其影響相對較小。

圖2 不同礦化度微咸水累積入滲量隨時間的變化

Fig.2 Relationship between cumulative infiltration and time under different mineralization degrees of infiltration water

由圖2可知，不同礦化度累積入滲量隨入滲時間的增加呈現冪函數增加趨勢，符合Kostiakov模型對累積入滲量I與入滲時間t之間關系的表述。因此，利用Kostiakov模型對試驗數據進行擬合，擬合結果如表3所示，其中R2為決定系數。

由表3可知，采用Kostiakov模型模擬不同礦化度微咸水灌溉下累積入滲量與時間關系時，決定系數均大于0.99，冪函數關系成立。在不同礦化度微咸水灌溉條件下，入滲指數α變化較小，在0.498～0.514之間。但淡水與微咸水的入滲指數α差異明顯，淡水低于微咸水。隨著入滲水礦化度的增加，累積入滲量也隨之增加?？傮w而言，在不同入滲水礦化度下，入滲系數K變化幅度不大。

表3 不同礦化度下累積入滲量與時間t擬合結果 Table 3 The fitting results between cumulative infiltration and time

2.2 微咸水礦化度對濕潤鋒運移深度的影響

濕潤鋒指在濕潤帶的末端，土壤含水量突變，與干層土有明顯的界面。不同礦化度微咸水入滲過程中濕潤鋒隨入滲時間的變化過程見圖3。不同礦化度微咸水在入滲初期，濕潤鋒運移深度增加較快，當入滲持續一段時間后，濕潤鋒的增加速率逐漸變小。在相同入滲歷時內，隨著入滲水礦化度的增加，濕潤鋒運移深度也相應增大。這表明：在微咸水入滲過程中，在一定濃度范圍內，隨著微咸水礦化度的增加，土壤中的鹽分濃度不斷增大，土壤顆粒之間的排斥力降低，進而促進顆粒間絮凝，改變了土壤團粒結構，增加了土壤顆粒的團聚性，使土壤中大孔隙增加，提高了水分在土壤中的入滲能力。

圖3 不同礦化度微咸水濕潤鋒隨時間的變化

Fig.3 Relationship between wetting front depth and time under different mineralization degree of infiltration water

2.3 濕潤鋒與累積入滲量的關系

由圖3可知，濕潤鋒與時間亦呈冪函數變化趨勢關系,對不同礦化度微咸水濕潤鋒與時間進行了擬合，得到冪函數系數A，β，擬合結果見表4。

由表4可知，微咸水灌溉下，t的指數非常小，可近似認為利用微咸水灌溉時，濕潤鋒與累積入滲量成線性關系，即H=CI。

表4 濕潤鋒與累積入滲量的關系 Table 4 Relationship between wetting front of soil and cumulative infiltration volume

對濕潤鋒與累積入滲量關系進行擬合，見表5。隨著入滲水礦化度的增加，擬合系數C也隨之增大，即累積入滲量隨著礦化度的增加而增大。擬合結果與試驗所得結果一致。

表5 濕潤鋒與累積入滲量H-I的擬合結果 Table 5 The fitting results between wetting front depth and cumulative infiltration

2.4 微咸水礦化度對土壤含水率分布特征的影響

土壤含水率指近地表層的含水率。不同礦化度微咸水灌溉條件下，實測土壤含水率與土層深度間變化關系見圖4。由圖4可知，由于入滲過程中土壤表面積水的存在，土壤表層含水率基本為飽和含水率，隨著入滲深度的增加，土壤質量含水率呈陡峭的下降趨勢，在濕潤鋒處達到最小值。隨著微咸水礦化度的增加，同一土層的含水率呈現增大的趨勢。礦化度為4 g·L-1的微咸水入滲結束后土壤含水率最大，大于2，3，5 g·L-1,均遠遠大于淡水。這主要是由于隨著微咸水礦化度的增加，由入滲水帶入土壤中的Na+數量隨之增加，Na+使土壤顆粒收縮，膠體顆粒的分散和膨脹，破壞土壤團聚體結構，導致土壤導水能力降低，從而導致土壤含水率減少。

2.5 入滲率模型分析

土壤入滲率又稱土壤滲透速率，其在任意時刻的大小i(t)與此時代表處的土壤水分通量q(0,t)相等，反映了土壤的入滲性能。

目前的土壤入滲模型大多針對淡水入滲過程建立起來的，這些模型是否適合描述微咸水入滲過程需要深入研究。因此本研究分別用philip模型，Green-Ampt模型，Kostiakov模型對不同礦化度下微咸水入滲試驗數據進行擬合，比較三種模型的適用性。采用三種模型擬合后的入滲參數、相關性以及顯著性見表6。

圖4 不同礦化度下土壤含水率隨土層深度的變化

Fig.4 Moisture content of soil profile under different mineralization degree of infiltration water

從表6可以看出，采用philip，Green-Ampt，Kostiakov模型對不同礦化度下微咸水入滲試驗進行模擬，決定系數都在0.94以上。其中，Philip模型，Green-Ampt模型在對淡水入滲過程擬合時，得到的穩定入滲率A，土壤表征飽和導水率Ks,濕潤峰處的土壤水吸力Sf為負數，與實際情況不符。Kostiatov模型中，入滲指數α隨入滲水礦化度的增加而降低，入滲系數k隨入滲水礦化度的增加呈現增加的趨勢，說明微咸水礦化度越高，入滲率衰減越慢。與Green-Ampt模型和Philip模型相比，Kostiakov公式的擬合結果更符合于入滲的真實過程。

3 結論與討論

山東省東營市重度鹽堿土壤在華北地區非常典型，本文對微咸水一維垂直入滲試驗中微咸水礦化度對重度鹽堿土壤入滲特征的影響進行了研究。通過研究得到以下結論：

1) 累積入滲量和濕潤鋒運移深度隨微咸水礦化度的增加呈現增大的趨勢，但是增加幅度在減小。不同礦化度微咸水濕潤鋒與累積入滲量呈線性關系，累積入滲量、濕潤鋒運移深度與入滲時間均呈冪函數關系。

表6 三種模型擬合分析結果 Table 6 The analysis results of three models

注 Note：*0.01

2) 不同礦化度微咸水土壤含水率隨著土層深度的增加呈現減小的趨勢。隨著入滲的進行，濕潤鋒不斷前移，含水率的分布曲線由比較陡直逐漸變為相對平緩。當微咸水礦化度增加時，同一土層的含水率亦呈現增大的趨勢。

3) 采用Kostiakov模型，Philip模型和Green-Ampt模型模擬了微咸水在重度鹽堿土中的入滲過程。研究結果表明，Kostiakov模型可以很好地描述重度鹽堿土壤入滲率與入滲時間的關系。

本研究對于華北地區利用微咸水進行灌溉具有一定的指導意義。重度鹽堿耕地土壤在不同微咸水灌溉模式下的水鹽運移規律，將在后續研究中進一步分析。

[1] 張啟海，周玉香.微咸水灌溉發展的基礎與措施探討[J].中國農村水利水電,1998,(10):12-13.

[2] 吳忠東，王全九，蘇 瑩，等.不同礦化度微咸水對土壤入滲特征的影響[J].人民黃河,2005,27(12):49-50.

[3] 劉婷珊，高艷明，李建設.不同礦化度微咸水入滲下土壤水鹽運移特征研究[J].灌溉排水學報,2014,33(3):68-72.

[4] 郭太龍，遲道才，王全九，等.入滲水礦化度對土壤水鹽運移影響的試驗研究[J].農業工程學報,2005,21(Z1)：84-87.

[5] 蘇 瑩.微咸水礦化度對土壤水鹽運移的影響研究[J].楊凌職業技術學院學報,2008,6(7):1-4,17.

[6] 任長江，白 丹，周 文，等.土壤初始含鹽量對水分入滲特性的影響[J].干旱區研究,2014,31(2):222-225.

[7] 畢遠杰，孫西歡，馬娟娟，等.不同入滲水頭作用下一維垂直入滲特性[J].山西水利,2006，(2)：75-78.

[8] 馬文君，程琴娟，李良濤，等.微咸水灌溉下土壤水鹽動態對作物產量的影響[J].農業工程學報，2010,26(1):73-80.

[9] 王克全,何新林,王振華，等.不同灌水處理對滴灌春小麥生長及產量的影響研究[J].節水灌溉,2010,(9):41-42.

[10] 王全九，單魚洋.微咸水灌溉與土壤水鹽調控研究進展[J].農業機械學報,2015,46(12):117-126.

[11] 栗 濤,王全九.淡水入滲模型對微咸水入滲過程的適用性[J].排灌機械工程報,2014,32(1):80-85.

[12] 吳忠東,王全九.不同初始含水率條件下的微咸水入滲試驗[J].農業機械學報,2010,41(增刊):53-58.

[13] Dudely L M, Ben-Gal A, Shani U. Influence of plant,soil,and water on the leaching fraction[J]. Vadose Zone Journal, 2008,7(2):420-425.

[14] 樊麗琴，楊建國，尚紅鶯，等.淋洗水質和水量對寧夏龜裂堿土水鹽運移的影響[J].水土保持學報，2015,29(6):258-262.

[15] 陳書飛，何新林，汪宗飛，等.微咸水滴灌研究進展[J].節水灌溉，2010,(2)：6-9.

[16] 畢遠杰，王全九，雪 靜.微咸水造墑對油葵生長及土壤鹽分分布的影響[J].干旱區資源與環境，2010,24(3)：183-189.

[17] 邵建榮，張鳳華，董 艷，等.干旱區微咸水滴灌條件下典型土壤鹽堿化影響因素研究[J].干旱地區農業研究,2015,33(6):216-221.

[18] 胡文明.微咸水灌溉對作物生長影響的試驗研究[J].灌溉排水學報,2007,26(1)：86-88.

[19] 寇偉峰，劉兆普，陳銘達，等.不同濃度海水對油葵幼苗光合作用和葉綠素熒光特性的影響[J].西北植物學報,2006,26(1):73-77.

[20] 來劍斌.土壤水分運特征及其參數確定[M].西安:西安理工大學,2003.

[21] 水利部農村水利司，水利部農田灌溉研究所.SL207-98.節水灌溉技術規范[S].北京：中國標準出版社，1998.

[22] Or D, Silva H R. Prediction of surface irrigation advance using soil intake properties[J]. Irrigation Science, 1996,16(4):159-167.

[23] Horton R E. An Approach toward a physical interpretation of infiltration-capacity[J]. Soil Sci，1940，5(3):399-417.

[24] Green W, Ampt G. Studies on soil physics. The flow of air and water through soils[J]. J Agric Sci, 1911,4(1):1-24.

[25] Philip J. The theory of infiltration :Theinfilotration equation and its solution[J]. Soil Science, 1957,83(5):345-357.

[26] Smith RE. The infiltration envelope:Results from a theoretical infiltrometer[J]. Journal of Hydrology, 1972,17:1-22.