基于HYDRUS模型的暗管排水水鹽運(yùn)移模擬

石培君，劉洪光，何新林，李 紅，李開(kāi)明

(石河子大學(xué)水利建筑工程學(xué)院/現(xiàn)代節(jié)水灌溉兵團(tuán)重點(diǎn)實(shí)驗(yàn)室，新疆 石河子 832000)

土壤鹽漬化是當(dāng)前人類(lèi)面臨的一個(gè)全球性生態(tài)環(huán)境問(wèn)題[1-2]，尤其在干旱半干旱地區(qū)，水資源短缺和土壤鹽漬化已成為制約灌溉農(nóng)業(yè)可持續(xù)發(fā)展的關(guān)鍵因子，也是造成土地退化的重要原因[3-4]。新疆作為干旱半干旱農(nóng)業(yè)大區(qū)，其鹽漬化耕地面積占新疆總耕地面積的30.2%，已嚴(yán)重制約農(nóng)業(yè)的發(fā)展。膜下滴灌技術(shù)具有節(jié)水抑鹽、增溫保墑、調(diào)控土壤鹽分的作用，使作物根區(qū)處于脫鹽狀態(tài)，為作物生長(zhǎng)提供了相對(duì)良好的微環(huán)境，在干旱區(qū)和鹽堿化耕地具有應(yīng)用價(jià)值[5-7]。據(jù)統(tǒng)計(jì)，目前新疆膜下滴灌種植面積已達(dá)到3.3×106hm2，種植覆蓋率超過(guò)90%。許多學(xué)者對(duì)其水鹽運(yùn)移規(guī)律研究發(fā)現(xiàn)，膜下滴灌由于灌溉量小，淋溶作用微弱，一般不會(huì)產(chǎn)生深層滲透，很難利用灌溉水淋洗鹽分[8-9]，只有根區(qū)范圍處于脫鹽區(qū)，而其下土層處于積鹽區(qū)，其鹽分僅在土層中轉(zhuǎn)移而無(wú)法消除，在強(qiáng)蒸發(fā)條件下出現(xiàn)了次生鹽漬化威脅[10-11]，對(duì)干旱區(qū)綠洲生態(tài)環(huán)境的和諧發(fā)展產(chǎn)生較大影響。

暗管排水被認(rèn)為是地下水平衡和農(nóng)田水管理實(shí)踐中最適合的方法[12]，同時(shí)，也是防止和治理鹽堿化最直接、最有效的措施之一[13-14]。張金龍等[15]研究了暗管排水間距對(duì)濱海鹽土淋洗脫鹽效果的影響，得出排水間距對(duì)田面各點(diǎn)土壤脫鹽效果影響較大。張展羽等[16]在濱海鹽堿地模擬研究了不同埋深和間距的暗管排水系統(tǒng)的地下水埋深動(dòng)態(tài)變化，得出暗管技術(shù)參數(shù)對(duì)地下水埋深的動(dòng)態(tài)影響較大。王洪義等[17]在地下水埋深為1.5 m 的鹽堿地進(jìn)行不同埋深和間距的暗管排水試驗(yàn)，結(jié)果表明暗管間距5 m、埋深0.8 m為大慶地區(qū)蘇打鹽堿土降漬脫鹽的最佳工程布設(shè)參數(shù)。王振華等[18]在地下水埋深大于4 m 的鹽堿地上，研究了滴灌淋洗條件下不同暗管埋深和間距對(duì)土壤脫鹽效果的影響，結(jié)果表明暗管間距為15 m時(shí)土壤脫鹽效果較好。劉玉國(guó)等[19]研究了干旱區(qū)暗管排水對(duì)不同程度鹽堿化土壤脫鹽效果的影響，結(jié)果表明輕度和中度鹽漬化農(nóng)田鹽分由表聚型向脫鹽型轉(zhuǎn)變。李顯溦等[20]利用HYDRUS軟件對(duì)暗管排水的水鹽運(yùn)動(dòng)參數(shù)進(jìn)行了校驗(yàn)，結(jié)果表明模擬值與實(shí)測(cè)值吻合度較高，可以較好地描述暗管排水、排鹽過(guò)程中的土壤水鹽動(dòng)態(tài)。李亮等[21]利用HYDRUS模型對(duì)土壤水鹽的遷移進(jìn)行了模擬分析，結(jié)果表明模型對(duì)土壤含水率和含鹽率運(yùn)移的模擬具有較高精度，反映出鹽分積聚和水分運(yùn)移規(guī)律。

上述研究均表明暗管排水可有效降低地下水位和土壤鹽分含量，土壤脫鹽效果較好。但暗管排水與膜下滴灌相結(jié)合的水鹽運(yùn)移問(wèn)題，至今尚不清楚，阻礙了暗管排水技術(shù)在膜下滴灌條件下的推廣和應(yīng)用，而模型模擬是研究和解決該問(wèn)題的重要工具。因此，為探索膜下滴灌條件下暗管排水的水鹽分布特征及運(yùn)移規(guī)律，本文利用HYDRUS模型對(duì)暗管排水模型試驗(yàn)的水鹽運(yùn)移規(guī)律進(jìn)行模擬驗(yàn)證，分析膜下滴灌暗管排水條件下土壤水鹽分布特征和動(dòng)態(tài)變化規(guī)律，揭示膜下滴灌條件下暗管排水排鹽效果，為干旱區(qū)暗管排水技術(shù)的推廣提供理論支撐。

1 材料與方法

1.1 試驗(yàn)材料

本試驗(yàn)在新疆石河子大學(xué)水利建筑工程學(xué)院水利與土木工程實(shí)驗(yàn)中心進(jìn)行，試驗(yàn)用模型由混凝土砌筑而成，長(zhǎng)、寬、高分別為5 m、2 m、1.8 m。模型裝置分為土槽和蓄水池，長(zhǎng)、寬、高分別為5 m、1 m、1.8 m，底部以連通管模擬控制地下水位，其管徑為3.2 cm，管壁上開(kāi)5 mm圓孔，開(kāi)孔率為10%，包裹透水無(wú)紡布。埋設(shè)2根暗管，管徑為9 cm，埋深1 m，間距4 m，傾斜度為5‰，暗管上開(kāi)5 cm×0.2 cm矩形孔，開(kāi)孔率為6%。暗管外包裹一層透水無(wú)紡布，無(wú)紡布外鋪設(shè)3～5 cm細(xì)砂和粗砂，細(xì)砂平均粒徑為0.35～0.5 mm，粗砂平均粒徑為5～10 mm，細(xì)砂位于無(wú)紡布周?chē)?，粗砂鋪設(shè)于細(xì)砂上方。連通管和暗管均為PVC管。采用單翼迷宮式滴灌帶，滴頭流量為3.6 L·h-1，布管方式采用大田普遍使用的一膜三管鋪設(shè)方式，共鋪設(shè)6條滴灌帶。通過(guò)塑料水桶與水泵形成供水系統(tǒng)，利用水表控制灌溉水量，暗管排水模型剖面示意圖如圖1所示。

試驗(yàn)用土為新疆生產(chǎn)建設(shè)兵團(tuán)農(nóng)八師146團(tuán)鹽荒土，土壤初始含鹽量為14 g·kg-1，根據(jù)中國(guó)鹽土分類(lèi)標(biāo)準(zhǔn)，該土壤為重度鹽堿土。試驗(yàn)前將土風(fēng)干，過(guò)5 mm土篩除去雜質(zhì)，將土壤混和均勻形成均質(zhì)土壤，分層進(jìn)行裝填，每10 cm夯實(shí)一次，保證每次壓實(shí)度基本相同，裝填后土壤容重為1.48 g·cm-3，第一次灌水后略有沉降為8 cm，經(jīng)計(jì)算土壤容重達(dá)到了1.52 g·cm-3。試驗(yàn)用水為地下水，其礦化度為0.25 g·L-1，灌溉結(jié)束后地下水礦化度增加到2.7 g·L-1。

圖1 暗管排水模型剖面示意圖Fig.1 Diagram of model profile subsurface drainage pipe

1.2 試驗(yàn)設(shè)計(jì)及方法

試驗(yàn)前在蓄水池中注入水，水位控制在1.2～1.4 m范圍內(nèi)，保證下層土壤基本接近飽和狀態(tài)。試驗(yàn)開(kāi)始通過(guò)水泵持續(xù)灌水，至暗管開(kāi)始排水即停止灌水并記錄灌水量。為保證水分和鹽分在土壤內(nèi)重新分布且時(shí)空變化較小，每次灌水時(shí)間間隔基本相同，共進(jìn)行3次灌水淋洗，灌溉排水參數(shù)如表1所示。暗管排水結(jié)束后24 h內(nèi)在水平方向距離暗管0、40、80、120、160、200 cm處分別進(jìn)行取樣，縱向每隔10 cm取一個(gè)樣，取樣總深度為80 cm，取樣點(diǎn)如圖1所示。將所取土樣在105℃下烘干8 h，用烘干法測(cè)定其含水量，將烘干土樣研磨、過(guò)0.5 mm土壤篩，按土水比為1∶5配置成溶液，用玻璃棒攪拌、靜置使上層液澄清，用電導(dǎo)率儀測(cè)定其電導(dǎo)率，用干燥殘?jiān)ù_定土壤含鹽量與電導(dǎo)率之間的標(biāo)定關(guān)系式，如公式(1)所示。

y=0.0035EC+0.8972

(1)

式中,y為土壤含鹽量(%)；EC為電導(dǎo)率(μs·cm-1)。

1.3 模型建立

1.3.1 HYDRUS 模型介紹 HYDRUS軟件是一種可用來(lái)模擬水流和溶質(zhì)在非飽和多孔隙介質(zhì)中運(yùn)移的數(shù)值模型[22]，是利用土壤物理參數(shù)模擬水、熱及溶質(zhì)在非飽和土壤中的一維(二維)運(yùn)動(dòng)的有限元計(jì)算機(jī)模型[23]。該模型能夠較好地模擬水分、溶質(zhì)與能量在土壤中的分布以及時(shí)空變化和運(yùn)移規(guī)律。

1.3.2 模型的基本方程

土壤水分運(yùn)動(dòng)方程

(2)

式中,θ為土壤體積含水率(cm3·cm-3)；K(θ)為非飽和土壤導(dǎo)水率(cm·d-1)；t為時(shí)間(d)，z為土壤深度(cm)。

表1 灌溉排水參數(shù)

土壤水力參數(shù)van Genuchten模型

(3)

(4)

式中,Ks為土壤飽和導(dǎo)水率(cm·d-1)；θe為土壤相對(duì)飽和度(cm3·cm-3)；θr為土壤殘余體積含水率(cm3·cm-3)；θs為土壤飽和體積含水率(cm3·cm-3)；θ(h)為土壤體積含水率(cm3·cm-3)；h為負(fù)壓水頭(cm)；K(θ)為土壤非飽和導(dǎo)水率(cm·h-1)；n、m、α均為經(jīng)驗(yàn)參數(shù)。其中m=1-1/n，α是與土壤物理性質(zhì)有關(guān)的參數(shù)；l為經(jīng)驗(yàn)擬合參數(shù)，通常取平均值0.5。

鹽分運(yùn)移基本方程

(5)

式中,C為溶質(zhì)濃度(g·cm-3)；qi為水通量(cm·d-1)；Dij為擴(kuò)散度(cm2·d-1)；xi為空間坐標(biāo)(i=1，2)，x1=x，x2=z，D11=Dxx，D12=Dxz。

1.4 模型的初始條件和邊界條件

土壤水分上邊界采用零通量邊界，由于地膜的覆蓋忽略蒸發(fā)的影響。滴頭處采用第一類(lèi)邊界條件，輸入通過(guò)滴頭處的通量值，包括灌水量，其他邊界均為零通量邊界。由于在試驗(yàn)期間沒(méi)有降水，因此忽略降雨量的影響。上邊界初始含水率為0.025 cm3·cm-3，下邊界初始含水量為0.41 cm3·cm-3，上下邊界初始含鹽量均為14 g·kg-1，灌水過(guò)程中，上邊界為常壓，待灌水結(jié)束，上邊界為零通量。整個(gè)灌溉過(guò)程中，暗管邊界設(shè)為滲透邊界，將土壤含水率和土壤含鹽量的實(shí)測(cè)值及模擬值分別進(jìn)行對(duì)比，相應(yīng)調(diào)整土壤水力特性參數(shù)，使兩者充分接近為止，修正的模型參數(shù)如表2所示。修正后的基本參數(shù)不變，以每次灌水量和灌水時(shí)間的試驗(yàn)數(shù)據(jù)進(jìn)行模型的驗(yàn)證。模擬值和實(shí)測(cè)值的吻合程度采用決定系數(shù)R2和均方根誤差RMSE指標(biāo)進(jìn)行評(píng)價(jià)，計(jì)算公式如下:

(6)

式中,Si,Mi分別為模擬值和實(shí)測(cè)值；N為數(shù)值比較值，無(wú)量綱。

表2 土壤物理特性參數(shù)

注：θr為殘余土壤體積含水率；θs為飽和土壤體積含水率；Ks為飽和導(dǎo)水率；α、n為土壤水分特征曲線擬合參數(shù)；l為孔隙連通性參數(shù)。

Note:θris residual soil water content;θsis saturated soil volume water content;Ksis saturated water conductivity;αandnare fitting parameters of soil moisture characteristic curve;lis pore connectivity parameter.

1.5 模型模擬

為探索干旱區(qū)暗管排水條件下土壤水分和鹽分的動(dòng)態(tài)變化規(guī)律，將土槽簡(jiǎn)化為二維模型進(jìn)行計(jì)算模擬，計(jì)算區(qū)域?yàn)橐粋€(gè)高(垂直方向)180 cm、長(zhǎng)(水平方向)500 cm 的矩形區(qū)域。模型模擬0～80 cm深度范圍土壤水分和鹽分變化特征，模擬時(shí)間共計(jì)45 d，采用變時(shí)間步長(zhǎng)剖分方式，根據(jù)收斂迭代次數(shù)調(diào)整時(shí)間步長(zhǎng)。設(shè)定初始時(shí)間步長(zhǎng)為0.001 d，最小步長(zhǎng)為0.0001 d，最大步長(zhǎng)為5 d；土壤含水量容許偏差為0.001，壓力水頭容許偏差為1 cm。

2 結(jié)果分析

2.1 灌水后的水鹽分布特征及變化規(guī)律

圖2(a)，(b)，(c)分別為3次灌水后不同剖面處的土壤鹽分動(dòng)態(tài)變化特征。由圖2可知，每次灌水后0～30 cm土層內(nèi)鹽分下降速率較快，30～80 cm土層內(nèi)鹽分下降速率緩慢。各土層內(nèi)鹽分含量隨土壤深度的增加而增加，距離暗管不同間距處各土層鹽分含量變化趨勢(shì)一致，但相同剖面處鹽分含量存在差異。水平距離暗管200 cm處各土層內(nèi)鹽分含量最大，這是由于土壤鹽分隨水流在運(yùn)動(dòng)過(guò)程中受到對(duì)流彌散作用的影響，且縱向彌散作用大于橫向彌散作用，水平距離暗管間距越大，其剖面內(nèi)鹽分被灌溉水淋洗到暗管內(nèi)的越少?？傮w來(lái)看，每次灌水后鹽分含量在垂直方向上隨土壤深度的增加鹽分含量下降速率減緩；水平方向上距離暗管間距越大，各剖面處鹽分含量相對(duì)較大，土體脫鹽率較小。3次灌水淋洗過(guò)程中，各剖面處土壤脫鹽率從大到小依次為第2次灌后、第1次灌后、第3次灌后，第1次灌水后水平距離暗管不同間距處鹽分含量存在明顯差異，這是由于試驗(yàn)土樣經(jīng)篩分、分層裝填后存在大孔隙，灌水后形成孔隙流使鹽分得不到充分的溶解，只有少量鹽分隨水流流入暗管，因此，第1次灌水后鹽分含量隨土壤深度的變化趨勢(shì)較為平緩。由圖2(b)、(c)可知，第2次灌水和第3次灌水后土壤鹽分含量隨土壤深度的增加呈現(xiàn)先減小后增加的趨勢(shì)，距離暗管不同間距處鹽分含量變化規(guī)律不明顯，距離暗管200 cm處各土層內(nèi)鹽分含量最大，表層0～20 cm土層內(nèi)鹽分含量明顯高于20～40 cm土層鹽分含量，說(shuō)明地下水埋深較淺時(shí)強(qiáng)烈蒸發(fā)很容易使鹽分隨水分通過(guò)毛細(xì)管上升至表層，出現(xiàn)明顯的鹽分表聚現(xiàn)象。20～40 cm土層內(nèi)鹽分含量下降速率較快，鹽分含量下降至2 g·kg-1左右，40～60 cm土層內(nèi)鹽分呈現(xiàn)逐漸增大的趨勢(shì)，60～80 cm土層內(nèi)鹽分含量處于穩(wěn)定變化的趨勢(shì)，土壤鹽分下降速率變緩，脫鹽效果不明顯，這是由于靠近暗管周?chē)寥纼?nèi)的水分含量接近飽和，水勢(shì)梯度減弱，減少了與地下水之間的水量交換和鹽分淋洗量；隨著灌溉次數(shù)的增加，地下水中鹽分含量增加，溶質(zhì)濃度變大，暗管周?chē)馏w與地下水中的鹽分交換量也減小，說(shuō)明地下水礦化度大小與土壤脫鹽率以及暗管排鹽量有著密切的關(guān)系。第3次灌水后，0～80 cm土壤內(nèi)總鹽分含量下降速率較快，各剖面處鹽分的質(zhì)量分?jǐn)?shù)下降至8 g·kg-1以下，表明暗管排水技術(shù)在膜下滴灌條件下的土壤水鹽運(yùn)移具有較好的變化規(guī)律，對(duì)土壤鹽分的淋洗效果較好。

圖3為3次灌水后各剖面處土壤水分動(dòng)態(tài)變化特征。由圖3可知，水平距離暗管不同間距處各土層內(nèi)水分含量隨土壤深度的增加變化趨勢(shì)基本一致，均隨土壤深度的增加而增加。0～20 cm土層內(nèi)水分含量變化較大，這是由于表層土壤受到強(qiáng)蒸發(fā)作用的影響，另外，上層土壤水分在重力勢(shì)的作用下不斷向下運(yùn)移。垂直方向20～80 cm土層內(nèi)各剖面處的水分變化趨勢(shì)一致，隨著土壤深度的增加逐漸增加，水分含量差異不明顯，80 cm土層內(nèi)含水量幾乎接近飽和狀態(tài)。水平方向距離暗管不同間距處各剖面的水分含量存在明顯差異，且沒(méi)有明顯的變化規(guī)律，這可能是受覆膜邊界條件、蒸發(fā)作用、地下水的影響。由圖3(b)、(c)可以看出，水平距離暗管不同間距處各土層內(nèi)含水量隨土壤深度的增加而增加，但水平距離暗管0 cm處在60～80 cm土層內(nèi)呈減小趨勢(shì)，說(shuō)明暗管排水可以降低暗管上方土壤水分含量?？傮w來(lái)看，3次灌水后暗管周?chē)寥纼?nèi)含水量較高，基本接近飽和含水率，垂直方向上變化規(guī)律較好，均隨土壤深度的增加而增加，水平方向上距離暗管不同間距處各剖面含水量差異不大。

圖2 3次灌水后土壤不同剖面處的鹽分含量變化Fig.2 Changes of salt content in different profiles of soil after three times of irrigation

圖3 3次灌水后土壤不同剖面處的水分含量變化Fig.3 Changes of water content in different profiles of soil after three times of irrigation

2.2 模型驗(yàn)證

根據(jù)構(gòu)建的數(shù)值模型，輸入土壤水力特性參數(shù)，同時(shí)設(shè)置時(shí)間、空間離散化處理參數(shù)后運(yùn)行計(jì)算，并進(jìn)行模型可靠性驗(yàn)證分析。通過(guò)比較第1次灌水、第3次灌水后距離暗管不同間距處各土層內(nèi)含水量和含鹽量的實(shí)測(cè)值與模擬值，兩者差異較小且總體變化趨勢(shì)一致，如圖4、5、6所示。第1次灌水后距離暗管不同間距處土壤鹽分模擬值與實(shí)測(cè)值的RMSE分別為0.364和0.632，R2分別為0.995和0.992，距離暗管不同間距處土壤水分模擬值與實(shí)測(cè)值的RMSE分別為1.024和1.324，R2分別為0.906和0.921；第3次灌水后距離暗管不同間距處土壤鹽分模擬值與實(shí)測(cè)值的RMSE分別為0.505和0.197,R2分別為0.996和0.994。由此看出，水平距離暗管不同間距處土壤鹽分和水分實(shí)測(cè)值與模擬值的RMSE均較小，R2均較大，表明兩者之間具有較好的一致性。從圖6可以看出，盡管水分變化在某些觀測(cè)點(diǎn)上存在一定的誤差，但總體規(guī)律基本一致。

圖4(a)和(b)表示第1次灌水后水平距離暗管0 cm和200 cm處各剖面鹽分含量實(shí)測(cè)值和模擬值的變化趨勢(shì)。從決定系數(shù)和均方根來(lái)看，模擬值和實(shí)測(cè)值吻合較好，第1次灌水結(jié)束后，水平距離暗管0 cm處0～30 cm土層內(nèi)鹽分變化較小，鹽分含量下降較大，實(shí)測(cè)值和模擬值存在較大的誤差，是由于灌水過(guò)程中強(qiáng)烈蒸發(fā)使鹽分隨毛管水上升至表面，造成實(shí)測(cè)值偏大，40～80 cm土層內(nèi)鹽分迅速增加，這是由于在灌溉水流的擴(kuò)散作用下，鹽分通過(guò)對(duì)流、彌散作用不斷向下遷移。水平距離暗管200 cm處0～40 cm土層內(nèi)鹽分含量呈線性增長(zhǎng)趨勢(shì)，實(shí)測(cè)值與模擬值吻合度較好，50～80 cm土層內(nèi)鹽分含量相對(duì)保持穩(wěn)定狀態(tài)。圖5(a)和(b)顯示第3次灌水后距離暗管不同剖面處鹽分含量實(shí)測(cè)值和模擬值數(shù)據(jù)，決定系數(shù)和均方根均顯示模擬值和實(shí)測(cè)值吻合很好，水平距離暗管200 cm處0～20 cm土層內(nèi)出現(xiàn)鹽分表聚現(xiàn)象，是由于水平距離暗管200 cm處沒(méi)有覆蓋地膜，蒸發(fā)作用使鹽分通過(guò)毛細(xì)管上升至表面，20～60 cm土層內(nèi)鹽分含量增加速率較快。水平距離暗管0 cm處60～80 cm土層內(nèi)鹽分含量呈下降趨勢(shì)，而水平距離暗管200 cm處60～80 cm鹽分含量保持不變。圖4和圖5相比較，可以看出經(jīng)過(guò)3次灌水后鹽分總含量明顯減小，說(shuō)明膜下滴灌與暗管排水相結(jié)合的技術(shù)可以有效地改良鹽堿地。

圖6(a)和(b)為第1次灌水后水平距離暗管0 cm和200 cm處各剖面土壤水分含量實(shí)測(cè)值和模擬值變化，由圖6可以看出，土壤水分的運(yùn)動(dòng)具有較好的規(guī)律，隨著土層深度的增加呈現(xiàn)逐漸增加的趨勢(shì)，水平距離暗管0 cm和200 cm處表層20 cm范圍內(nèi)土層含水量相差較大，這是由于水平距離暗管200 cm處沒(méi)有覆膜，受蒸發(fā)影響大，且暗管對(duì)水分的運(yùn)動(dòng)影響相對(duì)較小。從相關(guān)系數(shù)和均方根來(lái)看，實(shí)測(cè)值和模擬值吻合好，是因?yàn)橥寥篮Y分后，分層裝填，比較符合達(dá)西定律及土壤水鹽運(yùn)動(dòng)方程的均勻連續(xù)性假設(shè)。

2.3 模型應(yīng)用

利用實(shí)測(cè)值和模擬值修正后的參數(shù)進(jìn)行更長(zhǎng)時(shí)間序列的模擬計(jì)算，共進(jìn)行6次灌水的模擬，總時(shí)長(zhǎng)為90 d。圖7(a)和(b)表示距離暗管不同位置處各土層內(nèi)鹽分含量變化曲線。圖7(a)顯示0～40 cm土層內(nèi)鹽分含量<2 g·kg-1，達(dá)到了非鹽化土水平，40～80 cm土層內(nèi)鹽分含量隨著深度的增加呈對(duì)數(shù)型增長(zhǎng)趨勢(shì)，整體鹽分質(zhì)量分?jǐn)?shù)<8 g·kg-1，達(dá)到中度鹽化水平。圖7(b)顯示距離暗管不同間距處鹽分含量呈現(xiàn)波動(dòng)變化，這是由于覆膜造成邊界條件的變化，無(wú)膜覆蓋處土壤水分既受重力作用的影響，又受到大氣環(huán)境的影響，而覆膜位置處土壤水分受大氣環(huán)境的影響較小，運(yùn)移主要依靠重力勢(shì)和基質(zhì)勢(shì)，因此土壤中的可溶性鹽在對(duì)流、彌散和水流運(yùn)移作用下呈現(xiàn)波動(dòng)變化。在水平距離暗管間距150 cm范圍內(nèi)，20～40 cm深度土層內(nèi)鹽分含量趨于1 g·kg-1左右，60 cm深度土壤內(nèi)鹽分含量趨于2～4 g·kg-1變化，80 cm深度土層內(nèi)鹽分含量在6～8 g·kg-1范圍內(nèi)變化，各土層內(nèi)鹽分含量在水平距離暗管200 cm處(兩暗管中間位置)相對(duì)較大。

圖4 第1次灌水后土壤鹽分變化Fig.4 Changes of soil salt content after the first irrigation

圖5 第3次灌水后土壤鹽分變化Fig.5 Changes of soil salt content after the third irrigation

圖6 第1次灌水后土壤水分變化Fig.6 Changes of soil moisture after the first irrigation

圖7 模擬不同剖面出鹽分變化曲線Fig.7 Simulation curve of salt change at different profiles

3 結(jié) 論

1)通過(guò)膜下滴灌與暗管排水相結(jié)合的模型試驗(yàn)，經(jīng)過(guò)3次灌水淋洗后，土槽內(nèi)總鹽分含量降低較快，表層0～20 cm土壤鹽分含量下降速率最大，達(dá)到了非鹽化土水平，膜下滴灌與暗管排水技術(shù)的結(jié)合可以有效地降低土壤鹽分含量。

2)利用HYDRUS模型對(duì)膜下滴灌暗管排水條件下的水鹽運(yùn)移規(guī)律進(jìn)行了模擬，土壤水分含量和鹽分含量的實(shí)測(cè)值與模擬值具有很好的一致性，模型能夠較好地模擬土壤水鹽分布特征以及隨時(shí)間和土壤深度變化的規(guī)律，這將為干旱區(qū)膜下滴灌暗管排水條件下的水鹽運(yùn)移理論研究提供幫助。

3)基于修正的模型參數(shù)進(jìn)行HYDRUS模型的應(yīng)用，通過(guò)時(shí)長(zhǎng)為90 d的灌水模擬試驗(yàn)，發(fā)現(xiàn)表層40 cm范圍內(nèi)鹽分含量下降至2 g·kg-1，已達(dá)到非鹽化土水平。暗管排水條件下灌水排鹽可以有效地保證作物的正常生長(zhǎng)，膜下滴灌與暗管排水相結(jié)合的技術(shù)可高效地改良鹽堿地。