孫 娜，李瑞平，苗慶豐，范雷雷，周利穎，李聰聰

（內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院，呼和浩特010018）

0 引言

農(nóng)田土壤水分入滲過程是影響畦灌性能的重要因素之一，是灌區(qū)進(jìn)行地面灌溉設(shè)計的重要基礎(chǔ)，同時也是農(nóng)田水文關(guān)系演變、土壤侵蝕和營養(yǎng)物質(zhì)遷移轉(zhuǎn)化等領(lǐng)域研究的熱點[1-3]。由于農(nóng)田土壤水力特征、土壤結(jié)構(gòu)及人類活動等因素的影響，不同區(qū)域位置的差異導(dǎo)致土壤水分入滲過程表現(xiàn)出一定的時空差異性和尺度依賴性，進(jìn)而給農(nóng)田畦灌系統(tǒng)的設(shè)計帶來困難[4,5]。因此，準(zhǔn)確認(rèn)識畦田不同位置土壤水分入滲特性的差異，選取適宜模型進(jìn)行水分入滲過程模擬，對農(nóng)田合理灌溉和節(jié)水增產(chǎn)具有重要意義。

國內(nèi)外專家學(xué)者基于野外試驗、室內(nèi)模擬試驗和模型模擬等手段，開展了大量關(guān)于土壤水分入滲過程的研究，并取得豐富的理論成果[6-8]。賈宏偉等[9]利用簡化的Phillip 公式和土壤轉(zhuǎn)換函數(shù)研究了入滲系數(shù)α、穩(wěn)滲率fc的空間變化規(guī)律；聶衛(wèi)波等[10]將標(biāo)定理論與地統(tǒng)計學(xué)相結(jié)合對Phillip 入滲公式進(jìn)行標(biāo)定，研究相關(guān)標(biāo)定因子空間變異特征；Rodríguez Vásquez等[11]利用Phillip 和Kostiakov 模型研究入滲參數(shù)的空間變異性與土壤容重、土壤密度、沙粒、黏粒和石灰肥料之間的關(guān)系；李麗梅等[12]分析指出土壤入滲性能參數(shù)中穩(wěn)滲率和累積入滲量具有強烈的空間變異性；劉繼龍等[13]利用多重分析得出穩(wěn)定入滲率的空間變異性主要受土壤容重、粗粉粒和黏粒含量的影響；何丹等[14]研究指出土壤容重、有機(jī)質(zhì)和初始含水率是影響耕作地入滲過程的主要變異源。

綜上所述，土壤水分入滲空間變異性的研究主要集中于基于數(shù)理統(tǒng)計研究入滲能力有關(guān)指標(biāo)的差異性，而針對黃河水灌溉含沙量較高和土壤理化性質(zhì)差異明顯的典型區(qū)域土壤入滲過程的差異性研究較少。因此，本文選取河套灌區(qū)典型試驗區(qū)，開展畦田土壤水分入滲問題研究，通過對比分析試驗區(qū)田塊內(nèi)畦首、畦中和畦尾的入滲過程，揭示畦田田塊內(nèi)不同空間位置土壤水分入滲過程、入滲率、累積入滲量和入滲系數(shù)的差異性，并運用不同入滲模型對入滲過程進(jìn)行擬合，選擇最適宜試驗區(qū)畦田灌溉的土壤入滲模型，為深入認(rèn)識河套灌區(qū)畦田土壤水分入滲空間變異狀況提供基礎(chǔ)資料，為河套灌區(qū)農(nóng)業(yè)農(nóng)田灌溉和耕作提供理論參考。

1 材料與方法

1.1 試驗區(qū)概況

田間試驗區(qū)位于氣候、土壤和水鹽狀況較為典型的內(nèi)蒙古自治區(qū)巴彥淖爾市河套灌區(qū)雙河鎮(zhèn)進(jìn)步村六社酒莊試驗站（107°18′23.12″E，40°41′17.24″N），試驗區(qū)地貌多以平原為主，屬中溫帶大陸性半干旱氣候，降雨量少而蒸發(fā)較強，年均蒸發(fā)量2 032 mm，年均降水量僅為132 mm，降水量季節(jié)性差異特征明顯且年內(nèi)分配不均，雨季多集中在7-8月份，占全年降雨量的70%~80%，年平均日照時數(shù)為3 180 h，灌水生育期內(nèi)平均氣溫為20.9 ℃，無霜期約為126 d。試驗區(qū)土壤呈微堿性，土壤質(zhì)地類型分為2 層，上層0~20 cm 為粉土，土壤容重為1.42 g/cm3，田間持水量為22.35%；下層20~40 cm 為黏土，土壤容重為1.39 g/cm3，田間持水量為26.14%；試驗區(qū)土壤平均儲鹽量約為1 342 g/m3，為中度含鹽土壤。

1.2 試驗設(shè)計與數(shù)據(jù)處理

本文研究采用田間試驗與室內(nèi)試驗相結(jié)合的形式，于2017年4月對試驗田地進(jìn)行激光平地。為避免作物生長、施肥等因素的影響，分別于2017年10月和2018年10月在農(nóng)作物收割后進(jìn)行土樣采集和入滲試驗，并記錄采樣點周邊環(huán)境信息。綜合考慮試驗區(qū)地形和植被種植情況，選取長寬為45 m×21 m 的畦田田塊，將田塊按水流推進(jìn)方向分為畦首、畦中、畦尾，并按圖1所示布置入滲試驗點。參照圖1 中網(wǎng)格交點的位置，每個樣點重復(fù)取樣3次，使其具有代表性，能反映整塊畦田的土壤理化性質(zhì)分布格局。

圖1 試驗布置(單位：m)Fig.1 Test layout

土壤入滲過程采用雙環(huán)入滲儀進(jìn)行測定，雙環(huán)入滲儀內(nèi)環(huán)和外環(huán)直徑分別為26 cm 和60 cm，內(nèi)、外環(huán)高度均為25 cm。田間試驗時，首先清除表層土壤，將內(nèi)、外環(huán)埋入地面以下深度20 cm 處，然后關(guān)閉所有閥門，連接管線并給內(nèi)環(huán)供水，用馬氏瓶計量。入滲試驗采用當(dāng)?shù)攸S河水，測定含沙量為3%，試驗過程中水深控制在2 cm，內(nèi)外環(huán)水位齊平，外環(huán)水位采用水位平衡裝置控制，確保內(nèi)環(huán)屬于一維垂直入滲條件，減少積水水頭對水勢梯度的影響[15]。入滲水量分不同時段記錄，在試驗的0～15 min 每隔30 s 記錄一次；15～30 min 每隔1 min 記錄一次；30～60 min 每隔2 min 記錄一次；60～90 min每隔5 min記錄一次，90～120 min每隔10 min記錄一次。

土壤理化性質(zhì)測定指標(biāo)主要包括土壤質(zhì)地、土壤含水率、土壤含鹽量和土壤容重等。土壤含水率采用土鉆取土烘干法測定，土壤容重采用環(huán)刀法測量，土壤pH 和電導(dǎo)率（EC）按照水土比為5∶1 來配置樣品液，攪拌靜置后測定[16]；土壤樣品的機(jī)械組成采用MS 2000 型激光粒度儀測定，并按照國際制土壤粒徑分類標(biāo)準(zhǔn)把各粒級組含量定義成黏粒（＜0.002 mm）、粉粒（0.002～0.02 mm）和砂粒（0.02～2 mm）[17]。

本文采用GS+9.0 對土壤基本理化性質(zhì)的空間分布格局進(jìn)行研究分析其差異性特征；利用Orgin 2019軟件對土壤入滲過程進(jìn)行模型擬合及繪圖，選擇適宜本研究區(qū)的最優(yōu)入滲模型；采用SPSS 19.0 軟件分析入滲率與影響因素之間的相關(guān)性，探究土壤入滲空間變異性的原因。

1.3 研究方法

描述土壤水分入滲特征的指標(biāo)有初始入滲率、穩(wěn)定入滲率、90 min累積入滲量和入滲系數(shù)，土壤初始入滲率反映入滲初期單位時間內(nèi)地表單位面積土壤的入滲水量，穩(wěn)定入滲率反映入滲后期水流較為穩(wěn)定時的入滲速率，它等于或接近于飽和導(dǎo)水率，而土壤的入滲過程一般可在90 min 內(nèi)達(dá)到相對穩(wěn)定程度，因此，本文選取90 min 的累積入滲量作為反映土壤入滲能力的指標(biāo)[18,19]。分別采用具有代表性的Kostiakov 模型、Philip 模型、Horton 模型及通用經(jīng)驗?zāi)Ｐ蛯ν寥廊霛B過程進(jìn)行擬合分析[20]。

Kostiakov模型：

式中：I（t）為入滲速率，mm/min；t為入滲歷時，min；K、α為經(jīng)驗系數(shù)；f0為初始入滲率，mm/min。

Philip模型：

式中：f（t）為入滲速率，mm/min；t為入滲時間，min；S為土壤吸濕率，mm/min0.5；A為穩(wěn)定入滲率，mm/min。

Horton模型：

式中：f（t）為入滲速率，mm/min；t為入滲時間，min；f0為初始入滲率，mm/min；fc為穩(wěn)滲率，mm/min；k為模型參數(shù)，參數(shù)k決定f從f0減小到fc的速度。

通用經(jīng)驗?zāi)Ｐ停?/p>

式中：f（t）為入滲速率，mm/min；t為入滲時間，min；a、b為擬合的模型參數(shù)；n為經(jīng)驗系數(shù)。

2 結(jié)果與分析

2.1 土壤理化性質(zhì)特征分析

為了分析試驗區(qū)內(nèi)土壤理化性質(zhì)，直觀了解畦田內(nèi)土壤理化特性的空間分布狀況，現(xiàn)根據(jù)最優(yōu)半方差函數(shù)模型，選擇Kriging 插值法，利用GS+9.0 繪制了0~40 cm 耕作層土壤理化特性空間分布的平面等值線圖，如圖2所示。從圖2 可知，在整個試驗田塊范圍內(nèi)，土壤含水率呈現(xiàn)由田塊中心向四周逐漸增大的變化趨勢，含水率的變化范圍為10.98%~21.21%，含水率呈畦首＞畦尾＞畦中的分布特征，且含水率等值線圖比較密集，說明不同位置間含水率差異較明顯；試驗田塊內(nèi)土壤容重的變化范圍為1.50~1.63 g/cm3，呈條狀分布且空間范圍內(nèi)差異性較小，受大型機(jī)械壓實作用影響，土壤容重主要在畦首區(qū)域出現(xiàn)高值，而畦尾的土壤質(zhì)地疏松多孔，土壤容重主要出現(xiàn)低值；試驗田塊內(nèi)土壤pH 值為7.55~8.15，呈弱堿性，pH 高值區(qū)出現(xiàn)在試驗田塊中部區(qū)域，并向畦首和畦尾兩側(cè)依次呈現(xiàn)“高-低”交替的帶狀分布特征，空間分布差異性較明顯；試驗田塊內(nèi)土壤電導(dǎo)率呈斑塊狀分布，空間差異性較小，田塊內(nèi)3 個高值區(qū)其值分別為0.580、0.577 和0.478 mS/cm，并且分布于畦田兩側(cè)。土壤理化性質(zhì)特征受土壤質(zhì)地的分布狀況影響，試驗田塊內(nèi)土壤黏粒含量空間分布過渡比較平緩，土壤粉粒含量在畦中位置較高，土壤沙粒含量在畦首稍高，加之施肥、灌溉、翻耕等人為因素的影響，使得土壤理化性質(zhì)在空間分布上具有一定的差異性。

圖2 0~40 cm耕層土壤理化性質(zhì)等值線圖Fig.2 0~40 cm cultivated layer soil physical and chemical properties contour map

2.2 土壤入滲過程特征分析

為了分析試驗田塊內(nèi)土壤入滲過程的變化特征，分別繪制不用空間位置下土壤累積入滲量和入滲率曲線如圖3所示。從圖3可知，試驗田塊內(nèi)土壤累積入滲量和入滲速率隨時間的變化過程明顯并且在田塊內(nèi)不同位置呈現(xiàn)一定特征。畦首、畦中和畦尾的土壤水分入滲隨時間的變化過程基本一致，累積入滲量隨著入滲時間的推移不斷增大，而入滲率隨時間的增加而減小；土壤入滲初期入滲率較大，隨后快速衰減并趨于一個相對穩(wěn)定的入滲率，衰減速度呈逐漸減小的變化特征。累積入滲量呈現(xiàn)畦首＞畦中＞畦尾的變化特征，入滲率則呈畦首＞畦尾＞畦中的變化特征，且在入滲開始后，畦首的初始入滲率下降速度較快，約為畦尾初始入滲率下降速度的2倍；畦首的土壤入滲率達(dá)到相對穩(wěn)定所需的時間遠(yuǎn)小于畦中和畦尾，在入滲開始55 min 后，畦中和畦尾的入滲率基本趨于穩(wěn)定，而畦首的土壤入滲率仍在不斷變化，直到78 min趨于穩(wěn)定。

圖3 不同空間位置的土壤入滲變化過程Fig.3 Change process of soil infiltration in different spatial positions

現(xiàn)對比試驗田塊內(nèi)畦首、畦中和畦尾不同位置土壤入滲特性參數(shù)，結(jié)果如表1所示。由表1 可知，初始入滲速率、90 min 累積入滲量、入滲系數(shù)和穩(wěn)定入滲率均為畦首最大，畦首的初始入滲率分別是畦中和畦尾的1.42 和1.81 倍，入滲系數(shù)和穩(wěn)定入滲率呈現(xiàn)畦首＞畦尾＞畦中的變化特征，累積入滲量與初始入滲率呈現(xiàn)畦首＞畦中＞畦尾的變化特征。經(jīng)方差分析顯示，畦首的初始入滲率、穩(wěn)定入滲率和入滲系數(shù)均與畦中和畦尾具有明顯差異性（P＜0.05），而畦中與畦尾之間的初始入滲率、穩(wěn)定入滲率和入滲系數(shù)無明顯差異性。

表1 不同位置土壤入滲特性參數(shù)Tab.1 Soil infiltration parameters at different locations

2.3 土壤入滲過程模型擬合

土壤水分入滲是一個復(fù)雜的過程，國內(nèi)外眾多專家學(xué)者建立了多種入滲模型來模擬入滲速率隨時間變化的過程。具體依據(jù)模型是否具有物理意義劃分為2 類：一類是物理模型，主要包括Horton 模型和Philip 模型；另一類是經(jīng)驗?zāi)Ｐ?，主要包括Kostiakov－liews 模型和通用經(jīng)驗?zāi)Ｐ蚚21]。本文為進(jìn)一步研究試驗田塊內(nèi)不同位置土壤水分入滲過程，選用以上4個模型對畦首、畦中和畦尾的土壤入滲過程進(jìn)行模擬，模擬結(jié)果如表2和圖4所示。

表2 土壤入滲模型擬合參數(shù)結(jié)果Tab.2 Parameter fit results of soil infiltration model

從表2 和圖4 可知，4 個模型擬合度存在一定差異，但決定系數(shù)R2均在0.9 以上，說明擬合效果較好，綜合考慮4 個模型的決定系數(shù)R2，并將模擬參數(shù)與實測數(shù)據(jù)對比分析可知，Horton 模型的R2為0.988～0.997，較其他3 個模型高，且參數(shù)反映的結(jié)果與實測值的偏差較小，表明Horton 模型更適宜于描述本試驗區(qū)田塊內(nèi)土壤入滲速率隨時間變化的過程。在Horton 模型中，初始入滲速率f0與實測值偏差較小，穩(wěn)定入滲率fc呈現(xiàn)為畦首＞畦尾＞畦中的變化特征，與試驗實測結(jié)果一致；經(jīng)驗常數(shù)k決定著土壤水分入滲過程中達(dá)到穩(wěn)滲速率的快慢程度，k值越大則越快達(dá)到穩(wěn)定入滲速率，從擬合結(jié)果看出，畦尾最先達(dá)到穩(wěn)滲狀態(tài)，畦中次之，畦首最慢，與實測結(jié)果一致。

圖4 土壤入滲模型擬合曲線Fig.4 Fit diagram of soil infiltration model

2.4 土壤入滲空間變異性影響因素分析

為探究影響土壤入滲特性空間變異性的因素，現(xiàn)對反映土壤入滲特性的參數(shù)初始入滲率、穩(wěn)定入滲率和累積入滲量與土壤理化性質(zhì)指標(biāo)進(jìn)行相關(guān)性分析，結(jié)果如表3所示。由表3 可知，畦田土壤入滲特性與土壤理化性質(zhì)之間具有一定聯(lián)系。畦首、畦中和畦尾不同位置與土壤理化性質(zhì)之間相關(guān)性具有一定差異，其中土壤容重和含水率與不同位置土壤的入滲特性之間具有較高的相關(guān)性，說明土壤入滲過程明顯受到土壤容重和含水率的影響，土壤質(zhì)地和結(jié)構(gòu)發(fā)生變化，土壤質(zhì)地由重變輕，其入滲能力則逐漸增大[22]。初始入滲率在畦首位置與沙粒具有極顯著關(guān)系，主要由于畦首表層土壤含沙量較大，孔隙大，保水性能較差，導(dǎo)致其滲透速率較快；累積入滲量與畦尾的黏粒含量具有極顯著相關(guān)，黏粒含量越高，土壤孔隙較小，土壤通氣透水性能受到限制，從而導(dǎo)致入滲能力下降，累積入滲量越??；穩(wěn)定入滲率在畦中與土壤電導(dǎo)率呈極顯著負(fù)相關(guān)，隨電導(dǎo)率增加而減小，說明電導(dǎo)率高的土壤鹽分含量相對較多，而鹽分對土壤入滲特性具有一定的阻礙作用[23,24]。

表3 入滲特性參數(shù)與土壤特性相關(guān)關(guān)系Tab.3 Correlation between infiltration parameters and soil properties

3 結(jié)論

（1）試驗區(qū)田塊內(nèi)不同位置土壤理化性質(zhì)具有一定差異性，含水率呈畦首＞畦尾＞畦中的分布特征；土壤呈弱堿性，pH 在畦中較高且向畦首和畦尾兩側(cè)依次呈現(xiàn)“高-低”交替的帶狀分布特征；土壤電導(dǎo)率呈斑塊狀分布且空間差異性較?。煌寥廊葜爻蕳l狀分布且空間范圍內(nèi)差異性較小。試驗田塊內(nèi)不同位置土壤入滲特性與土壤理化性質(zhì)之間具有一定聯(lián)系，土壤容重和含水率對不同位置土壤的入滲特性影響明顯。

（2）試驗田塊內(nèi)不同位置土壤入滲特性具有一定差異性，畦首位置的初始入滲速率、90 min累積入滲量、入滲系數(shù)和穩(wěn)定入滲率最大，且初始入滲率分別是畦中和畦尾的1.42和1.81倍，入滲系數(shù)和穩(wěn)定入滲率呈現(xiàn)畦首＞畦尾＞畦中的變化特征，累積入滲量與初始入滲率呈現(xiàn)畦首＞畦中＞畦尾的變化特征。

（3）通過選用4個模型對試驗田塊內(nèi)不同位置土壤入滲過程進(jìn)行模擬可知，4個模型擬合度存在一定差異，但決定系數(shù)R2均在0.9 以上，擬合效果較好，其中Horton 模型模擬參數(shù)反映的結(jié)果與實測值的偏差較小，更適宜于描述本試驗區(qū)田塊內(nèi)土壤的入滲過程。