應(yīng)用HYDRUS-1D模擬砂質(zhì)夾層土壤入滲特性①

范嚴(yán)偉，黃 寧，馬孝義，畢貴權(quán)，趙文舉

(1 蘭州理工大學(xué)能源與動力工程學(xué)院，蘭州 730050；2 西部災(zāi)害與環(huán)境力學(xué)教育部重點實驗室/蘭州大學(xué)土木工程與力學(xué)學(xué)院，蘭州730000；3 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室，陜西楊凌 712100)

?

應(yīng)用HYDRUS-1D模擬砂質(zhì)夾層土壤入滲特性①

范嚴(yán)偉1,2，黃 寧2，馬孝義3，畢貴權(quán)1，趙文舉1

(1 蘭州理工大學(xué)能源與動力工程學(xué)院，蘭州 730050；2 西部災(zāi)害與環(huán)境力學(xué)教育部重點實驗室/蘭州大學(xué)土木工程與力學(xué)學(xué)院，蘭州730000；3 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室，陜西楊凌 712100)

摘 要：依據(jù)非飽和土壤水分運動理論，采用HYDRUS-1D軟件，對砂質(zhì)夾層土壤入滲特性進(jìn)行數(shù)值模擬，分析各因素對砂質(zhì)夾層土壤入滲規(guī)律的影響。結(jié)果表明：砂質(zhì)夾層結(jié)構(gòu)對土壤入滲特性有較大影響，具有暫時的阻水和減滲作用；濕潤鋒穿過砂層上界面后，入滲過程變?yōu)榉€(wěn)滲階段，穩(wěn)滲率主要受砂層質(zhì)地、砂層埋深和壓力水頭影響，與土壤初始含水率和砂層厚度無關(guān)；砂質(zhì)夾層土壤剖面水分分布不連續(xù)，上層土壤基本飽和，砂層土壤未飽和，土壤剖面含水率主要受砂層質(zhì)地、砂層埋深和砂層厚度的影響。研究結(jié)果可為農(nóng)業(yè)水資源利用及工程防滲技術(shù)提供理論依據(jù)。

關(guān)鍵詞：砂質(zhì)夾層；影響因素；入滲特性；數(shù)值模擬；HYDRUS-1D

由于風(fēng)沙、水文和地質(zhì)等外力作用，田間土壤剖面普遍存在砂質(zhì)夾層結(jié)構(gòu)[1-2]。夾層結(jié)構(gòu)會改變土壤水分的入滲性能和分布狀況，對土壤入滲特性有重要影響[3-6]。與均質(zhì)土壤相比，砂質(zhì)夾層土壤的入滲特性主要受砂層質(zhì)地、埋深、厚度以及土壤初始含水率和壓力水頭等因素影響。許多學(xué)者采用室內(nèi)試驗方法，研究了砂層質(zhì)地、埋深和厚度對土壤入滲特性的影響[7-9]。隨著計算機技術(shù)的發(fā)展，數(shù)值模擬方法越來越多地應(yīng)用于土壤水分運動的研究中，成為定量分析土壤水分入滲特性的有效工具[10-12]。其中，美國農(nóng)業(yè)部鹽漬土實驗室開發(fā)的 HYDRUS-1D 軟件，可對一維飽和-非飽和土壤中水、熱及溶質(zhì)運移進(jìn)行數(shù)值模擬[13]。HYDRUS-1D 應(yīng)用廣泛，具有良好的適用性，可用來模擬層狀土壤水分入滲特性[14-16]。本文依據(jù)非飽和土壤水分運動理論，采用 HYDRUS-1D 軟件，對不同影響因素組合下的砂質(zhì)夾層土壤入滲特性進(jìn)行數(shù)值模擬，分析各影響因素條件下砂質(zhì)夾層土壤的入滲規(guī)律，為農(nóng)田水資源高效利用和工程防滲技術(shù)提供理論依據(jù)。

1 層狀土一維垂直入滲數(shù)學(xué)模型

1.1 基本方程

HYDRUS-1D 中，一維飽和-非飽和土壤水運動采用 Richards 方程進(jìn)行描述。其表達(dá)式為：

式中：z為垂向坐標(biāo)，規(guī)定z向下為正(cm)；C(h)為土壤比水容量(1/cm)；h在飽和區(qū)和非飽和區(qū)分別為壓力水頭和基質(zhì)勢(cm)；t為入滲時間(min)；K(h)為土壤非飽和導(dǎo)水率(cm/min)。

式(1)中涉及的土壤水分特征曲線θ(h)和土壤非飽和導(dǎo)水率K(h)采用van Genuchten-Mualem模型擬合[17]。即：

式中：θr為土壤殘余含水率(cm3/cm3)；θs為土壤飽和含水率(cm3/cm3)；α、n和m為土壤物理特性有關(guān)的擬合參數(shù)；m=1-1/n；Ks為土壤飽和導(dǎo)水率(cm/min)；l = 0.5。

1.2 定解條件

圖 1 為砂質(zhì)夾層土壤水分運動模擬計算簡圖。試驗開始時，土壤水分剖面為穩(wěn)定剖面，計算域內(nèi)各點土水勢相等。即：式中：0Ψ為土壤初始總水勢(cm)。

試驗中，土柱上邊界保持恒定水頭入滲，下邊界入滲水量未到達(dá)。即：

式中：h0為壓力水頭(cm)；L為土柱高度(比濕潤鋒所濕潤范圍大)(cm)。

1.3 數(shù)值求解

利用 HYDRUS-1D 軟件，對層狀土一維垂直入滲的數(shù)學(xué)模型進(jìn)行求解。模擬不同砂層質(zhì)地、初始含水率、砂層厚度、砂層埋深和壓力水頭條件下的砂質(zhì)夾層土壤入滲特性。求解過程中，采用Galerkin 有限元法對土壤剖面進(jìn)行空間離散，采用隱式差分格式進(jìn)行時間離散。模擬土層深度為 100 cm，時間步長為0.01 min，空間步長為 1 cm，模擬歷時 200 min。土壤質(zhì)地通過 van Genuchten-Mualem 模型中的參數(shù)來體現(xiàn)，土壤初始含水率、砂層埋深和砂層厚度通過初始條件來設(shè)定，而壓力水頭通過邊界條件來實現(xiàn)。

為保證研究土壤的廣泛性和成果的普適性。模擬中土壤的 van Genuchten-Mualem 模型水力特性參數(shù)分別取自參考文獻(xiàn)[6]、[13]和[15]，如表 1 所示。

圖 1 砂質(zhì)夾層土壤計算簡圖Fig.1 Calculation diagram of the soil with sand interlayer

表1　不同土質(zhì)van Genuchten-Mualem模型水力特性參數(shù)Table 1 Hydraulic parameters in van Genuchten-Mualem model of different soils

2 各因素對砂質(zhì)夾層土壤入滲特性影響分析

2.1 砂層質(zhì)地

為對比分析砂層質(zhì)地對土壤水分入滲特性的影響，在上層土壤初始含水率θ0= 0.165 cm3/cm3，砂層埋深Z = 30 cm，砂層厚度D = 15 cm，壓力水頭h0= 6 cm條件下，模擬得到不同砂層質(zhì)地和均質(zhì)壤土的入滲特性曲線及入滲結(jié)束時(t = 200 min)土壤剖面水分分布曲線(圖2)。

由圖2可看出：砂質(zhì)夾層土壤入滲與均質(zhì)土入滲存在較大差異。入滲初期，砂質(zhì)夾層土壤的入滲規(guī)律和趨勢與均質(zhì)壤土相同，入滲速度快且迅速減小，累積入滲量和濕潤鋒運移距離隨時間成非線性增加；濕潤鋒到達(dá)砂層上界面時(入滲約52 min左右)，濕潤鋒稍有停滯，表現(xiàn)為阻水作用；濕潤鋒穿過砂層上界面后(入滲約64 min左右)，入滲率變?yōu)槌?shù)，入滲過程變?yōu)榉€(wěn)滲階段，穩(wěn)滲率小于相同時刻均質(zhì)壤土的瞬時入滲率，表現(xiàn)為減滲作用；隨著入滲時間的增加，穩(wěn)滲率逐漸大于均質(zhì)壤土的瞬時入滲率，表現(xiàn)為增滲作用。入滲結(jié)束時，上層壤土基本飽和，夾層砂土未飽和，下層壤土剖面水分分布的變化趨勢基本相同。

為定量分析各因素對砂質(zhì)夾層土壤入滲特性的影響，將穩(wěn)滲階段的累積入滲量采用線性關(guān)系表示[18]，即：

式中：I為累積入滲量(cm)；t1為濕潤鋒穿過砂層上界面的時間(min)；if為穩(wěn)滲率(cm/min)；I1為濕潤鋒穿過砂層上界面時的累積入滲量(cm)。

將模擬得到的不同砂層質(zhì)地下的入滲水量采用式(7)擬合，結(jié)果列于表2。

由表 2 可知：上述相關(guān)系數(shù)R2均大于 0.99，說明穩(wěn)滲階段的累積入滲量與入滲時間之間的關(guān)系均可以用線性關(guān)系表達(dá)。在相同上層土壤初始含水率、砂層埋深、砂層厚度和壓力水頭條件下，砂層質(zhì)地主要影響穩(wěn)滲率 if值。分析發(fā)現(xiàn)，穩(wěn)滲率 if與砂土 van Genuchten-Mualem 模型參數(shù)α的倒數(shù)成線性關(guān)系。對表 2 中的if和表 1 中的α擬合，得：

圖2　砂層質(zhì)地對砂質(zhì)夾層土壤入滲特性的影響Fig.2 Effect of sand layer texture on infiltration characteristics of the soil with sand interlayer

表2　不同砂層質(zhì)地下累積入滲量擬合參數(shù)Table 2 Fitting parameters of cumulative infiltration under different sand layer texture

2.2 土壤初始含水率

采用文獻(xiàn)[6]中的壤土和砂土，在D = 15 cm、Z = 30 cm、h0= 6 cm時，模擬得出不同上層土壤初始含水率條件下砂質(zhì)夾層土壤的入滲特性曲線及入滲結(jié)束時(t = 200 min)土壤剖面水分分布曲線，如圖3所示。

由圖3可看出：在其他條件相同時，土壤初始含水率對砂質(zhì)夾層土壤入滲特性影響較小。隨土壤初始含水率的增大，土壤累積入滲量略有減少，濕潤鋒運移距離稍有增大，主要是土壤初始含水率愈高，土水勢梯度小，入滲能力略有減小，而土壤易飽和，濕潤鋒運移有所加快。入滲結(jié)束時，隨土壤初始含水率的增大，砂土層含水率不飽和程度減小，下層壤土相同斷面處含水率稍有增大。

將模擬得到的不同初始含水率下的入滲水量采用式(7)擬合，結(jié)果列于表3。

由表3可知：土壤初始含水率主要影響 t1和I1值，而對穩(wěn)滲率 if無影響。初始含水率越大，濕潤鋒穿過砂層上界面的時間 t1越短，濕潤鋒穿過砂層上界面時的累積入滲量I1越小。分析發(fā)現(xiàn)，濕潤鋒穿過砂層上界面時的累積入滲量I1與土壤飽和差符合線性關(guān)系，對表3中的I1和θ0及表1中的θs擬合，得：

分析式(9)發(fā)現(xiàn)，I1與的斜率約等于砂層埋深Z。

2.3 砂層厚度

采用文獻(xiàn)[6]中的壤土和砂土，在θ0= 0.165 cm3/cm3、Z= 30 cm、h0= 6 cm時，模擬得出均質(zhì)壤土和不同砂層厚度條件下砂質(zhì)夾層土壤的入滲特性曲線及入滲結(jié)束時(t = 200 min)土壤剖面水分分布曲線，如圖4所示。

由圖4可看出：在其他條件相同時，砂層厚度對砂質(zhì)夾層土柱入滲特性影響微弱，達(dá)到穩(wěn)滲狀態(tài)時，穩(wěn)滲率相同，砂層厚度對累積入滲量無影響，減滲效果相同。隨砂層厚度的增大，濕潤鋒穿過砂層所需時間增加，穿過砂層后，運移距離基本相等。入滲結(jié)束時，砂層土壤剖面水分分布的變化趨勢相同，隨砂層厚度的增大，砂層土壤不飽和程度有所增大。

2.4 砂層埋深

采用文獻(xiàn)[6]中的壤土和砂土，在0θ= 0.165 cm3/cm3、D =15 cm、h0= 6 cm時，模擬得出均質(zhì)壤土和不同砂層埋深條件下砂質(zhì)夾層土壤的入滲特性曲線及入滲結(jié)束時(t = 200 min)土壤剖面水分分布曲線，如圖5所示。

將模擬得到的不同砂層埋深下的入滲水量采用式(7)擬合，結(jié)果列于表4。

圖3　土壤初始含水率對砂質(zhì)夾層土壤入滲特性的影響Fig.3 Effect of initial water content on infiltration characteristics of the soil with sand interlayer

表3　不同初始含水率下累積入滲量擬合參數(shù)Table 3 Fitting parameters of cumulative infiltration under different initial water content

由圖5及表4可知：在其他條件相同時，砂層埋深對砂質(zhì)夾層土壤入滲特性影響較大，砂層埋深越小，到達(dá)穩(wěn)滲階段的時間越短，減滲作用越早發(fā)生，但穩(wěn)滲率越大，導(dǎo)致先減滲后增滲現(xiàn)象發(fā)生。隨砂層埋深增大，濕潤鋒運移距離同樣存在先減小后增大的趨勢。砂層埋深對砂層土壤含水率分布有一定影響，砂層埋深較小時，砂層含水率隨深度逐漸增加，砂層埋深較大時，砂層含水率隨深度逐漸減少，且含水率不飽和程度增大。

分析發(fā)現(xiàn)，穩(wěn)滲率if與砂層埋深Z的倒數(shù)成線性關(guān)系。對表 4 中的if和Z擬合，得：

圖4　砂層厚度對砂質(zhì)夾層土壤入滲特性的影響Fig.4 Effect of sand layer thickness on infiltration characteristics of the soil with sand interlayer

圖5　砂層埋深對砂質(zhì)夾層土壤入滲特性的影響Fig.5 Effect of sand layer burial depth on infiltration characteristics of the soil with sand interlayer

濕潤鋒穿過砂層上界面時的累積入滲量I1與砂層埋深Z符合線性關(guān)系，對表 4 中的I1和Z擬合，得：

分析式(11)發(fā)現(xiàn)，I1與Z的斜率約等于土壤飽和差

綜合式(9)和式(11)，濕潤鋒穿過砂層上界面時的累積入滲量I1可表示為：

表4　不同砂層埋深下累積入滲量擬合參數(shù)Table 4 Fitting parameters of cumulative infiltration under different sand layer burial depth

2.5 壓力水頭

采用文獻(xiàn)[6]中的壤土和砂土，在0θ= 0.165 cm3/ cm3、D = 15 cm、Z = 30 cm時，模擬得出均質(zhì)壤土和不同壓力水頭條件下砂質(zhì)夾層土壤的入滲特性曲線及入滲結(jié)束時(t = 200 min)土壤剖面水分分布曲線，如圖6所示。

由圖6可看出：在其他條件相同時，壓力水頭對砂質(zhì)夾層土壤入滲特性影響較小，隨壓力水頭增大，土壤累積入滲量、穩(wěn)滲率、濕潤鋒運移距離和下層壤土相同斷面處含水率稍有增大，主要是由于土壤水分入滲主要由基質(zhì)勢、重力勢和壓力勢作用，壓力水頭大時，導(dǎo)致入滲下界面處壓力勢大，入滲稍有增大。壓力水頭越大，砂土層含水率不飽和程度越小。

將模擬得到的不同壓力水頭下的入滲水量采用式(7)擬合，結(jié)果列于表5。

圖6　壓力水頭對砂質(zhì)夾層土壤入滲特性的影響Fig.6 Effect of pressure head on infiltration characteristics of the soil column with sand interlayer

表5　不同壓力水頭下累積入滲量擬合參數(shù)Table 5 Fitting parameters of cumulative infiltration under different pressure head

由表 5 可知：壓力水頭主要影響t1和if值，對濕潤鋒穿過砂層上界面時的累積入滲量I1無影響。壓力水頭越大，濕潤鋒穿過砂層上界面的時間t1越短，穩(wěn)滲率if越大。分析發(fā)現(xiàn)，穩(wěn)滲率if與壓力水頭h0成線性關(guān)系。對表 5 中的if和h0擬合，得：

3 結(jié)論

1)砂質(zhì)夾層結(jié)構(gòu)對土壤入滲特性有較大影響，具有暫時的阻水和減滲作用。

2)濕潤鋒穿過砂層上界面時，上層土壤基本飽和，該時刻的累積入滲量等于砂層埋深與土壤飽和差的乘積。

3)濕潤鋒穿過砂層上界面后，入滲過程變?yōu)榉€(wěn)滲階段，該階段的累積入滲量隨時間成線性變化，穩(wěn)滲率受砂層質(zhì)地、砂層埋深和壓力水頭影響，與土壤初始含水率和砂層厚度無關(guān)。

4)砂質(zhì)夾層土壤剖面水分分布不連續(xù)，砂層土壤未飽和，土壤剖面含水率主要受砂層質(zhì)地、砂層埋深和砂層厚度的影響。

5)在農(nóng)業(yè)水資源利用和工程防滲技術(shù)中，應(yīng)重點考慮不同砂層質(zhì)地下的合理埋深，以提高其阻水減滲效果。

參考文獻(xiàn)：

[1]王文焰,王全九,沈冰,等.甘肅秦王川地區(qū)雙層土壤結(jié)構(gòu)的入滲特性[J].土壤侵蝕與水土保持學(xué)報,1998,4(2):36-40

[2]王雄師.疏勒河項目區(qū)土體構(gòu)型分類及對鹽堿地改良的影響[J].甘肅水利水電技術(shù),1998,34(3):66-69

[3]Hill D E,Parlange J Y.Wetting front instability in layered soils[J].Soil Science Society of America Journal,1972,36(5):697-702

[4]王文焰,張建豐,汪志榮,等.砂層在黃土中的阻水性及減滲性的研究[J].農(nóng)業(yè)工程學(xué)報,1995,11(1):104-110

[5]Yang H,Rahardjo H,Leong E C,et al.A study of infiltration on three sand capillary barriers[J].Canadian Geotechnical Journal,2004,41(4):629-643

[6]王春穎,毛曉敏,趙兵.層狀夾砂土柱室內(nèi)積水入滲試驗及模擬[J].農(nóng)業(yè)工程學(xué)報,2010,26(11):61-67

[7]王文焰,張建豐,汪志榮,等.黃土中砂層對入滲特性的影響[J].巖土工程學(xué)報,1995,17(5):33-41

[8]曲晨曉,王煒.土壤剖面中砂質(zhì)夾層的儲水作用及機理研究[J].華中農(nóng)業(yè)大學(xué)學(xué)報,1997,16(5):349-356

[9]邱玥,魏新平,廖華勝,等.夾砂層土壤水分入滲試驗研究[J].水資源與水工程學(xué)報,2009,20(1):120-123

[10]Kandelous M M,?im?nek J.Numerical simulations of water movement in a subsurface drip irrigation system under field and laboratory conditions using HYDRUS-2D[J].Agricultural Water Management,2010,97(7):1 070-1 076

[11]虎膽·吐馬爾白,吳爭光,蘇里坦,等.棉花膜下滴灌土壤水鹽運移規(guī)律數(shù)值模擬[J].土壤,2012,44(4):665-670[12]余根堅,黃介生,高占義.基于HYDRUS 模型不同灌水模式下土壤水鹽運移模擬[J].水利學(xué)報,2013,44(7):826-834

[13]?im?nek J,van Genuchten M T,?ejna M.Development and applications of the HYDRUS and STANMOD software packages and related codes[J].Vadose Zone Journal,2008,7(2):587-600

[14]毛曉敏,尚松浩.計算層狀土穩(wěn)定入滲率的飽和層最小通量法[J].水利學(xué)報,2010,41(7):810-817

[15]任利東,黃明斌,樊軍.不同類型層狀土壤持水能力的研究[J].農(nóng)業(yè)工程學(xué)報,2013,29(19):105-111

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

[17]Van Genuchten M T.A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J].Soil Science Society of America Journal,1980,44(5):892-898

[18]李毅,任鑫,Horton R.不同質(zhì)地和夾層位置對層狀土入滲規(guī)律的影響[J].排灌機械工程學(xué)報,2012,30(4):485-490

Simulation of Infiltration Characteristics in Soil with Sand Interlayer Using HYDRUS-1D

FAN Yanwei1,2,HUANG Ning2,MA Xiaoyi3,BI Guiquan1,ZHAO Wenju1
(1 College of Energy & Power Engineering,Lanzhou University of Technology,Lanzhou 730050,China; 2 Key Laboratory of Mechanics on Disaster and Environment in Western China,The Ministry of Education of China/ School of Civil Engineering and Mechanics,Lanzhou University,Lanzhou 730000,China; 3 Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas,Northwest A & F University,Yangling,Shaanxi 712100,China)

Abstract:Based on the theory for water movement in non-saturated soil,the HYDRUS-1D was applied to simulate the infiltration characteristics of the soil with sand interlayer.The impacts of various factors on the infiltration of soil with sand interlayer was analyzed.The results showed that:sand interlayer structure has great effect on soil infiltration characteristics,with a temporary water-blocking and infiltration-reducing.The process of infiltration reached the steady infiltration stage when the wetting front passed through the upper interface of sand layer.The steady infiltration rate was affected by the sand layer texture,sand layer burial depth and pressure head,and not affected by the initial soil water content and sand layer thickness.Soil water distribution in profile of the soil with sand interlayer was not continuous,the soil of upper layer was basically saturated,and soil of sand layer was unsaturated.Soil profile moisture was mainly affected by sand layer texture,sand layer burial depth and sand layer thickness.The results obtained in present study will provide theoretical basis for agricultural water resource utilization and engineering seepage control technique.

Key words:Sand interlayer; Influencing factors; Infiltration characteristics; Numerical simulation; HYDRUS-1D

作者簡介：范嚴(yán)偉(1982—)，男，山東聊城人，博士研究生，講師，主要從事農(nóng)業(yè)水土工程研究。E-mail:fanyanwei24@163.com

基金項目：①國家自然科學(xué)基金項目(51409137、51269008)和甘肅省自然科學(xué)基金項目(145RJYA293)資助。

DOI:10.13758/j.cnki.tr.2016.01.029

中圖分類號：S152.7