渠道滲漏HYDRUS模擬驗證及影響因素分析

付 強，李 玥，李天霄，崔 嵩，劉 東（東北農業(yè)大學水利與土木工程學院，哈爾濱 150030）

渠道滲漏HYDRUS模擬驗證及影響因素分析

付 強，李 玥，李天霄，崔 嵩，劉 東
（東北農業(yè)大學水利與土木工程學院，哈爾濱 150030）

渠道滲漏影響灌溉水利用效率，分析影響滲漏的主要因素，評估各影響因素對渠道滲漏的影響程度，可以為提高灌溉水利用效率提供理論參考。為了得到渠道滲漏的主要影響因素，進行了渠道滲漏室內模擬試驗，觀測不同邊坡系數、入滲水頭、渠道底寬下的累積入滲量和濕潤鋒推進過程，建立試驗條件下的土壤水分運動模型，選用HYDRUS-2D軟件求解所建模型，并將求解得到的模擬值與實測值進行比較，結果表明，模型模擬R2均高于0.9、均方根誤差為0.43～4.99 L，檢驗了所建模型可以用于模擬渠道滲漏的土壤水分入滲過程。模擬表明，邊坡系數對濕潤鋒運移距離和累積入滲量的影響很小，渠道底寬對渠道土壤水分入滲過程的影響較大，通過雙因素方差分析，渠道底寬對累積入滲量和濕潤鋒運移距離的影響極顯著（P＜0.01），而入滲水頭對水平距離的影響顯著（P＜0.05），渠道底寬和入滲水頭之間不存在顯著的交互影響（P＞0.05）。研究結果可以為提高灌溉水利用效率、對灌溉渠道防滲提供技術參考。

渠道；入滲；數值分析；HYDRUS-2D；濕潤鋒；滲漏

0 引 言

灌溉水損失的一大部分來自渠道滲漏，渠道水滲漏會降低灌區(qū)的灌溉水利用效率，造成水資源的浪費。由于通過大田試驗研究渠道滲漏土壤入滲規(guī)律比較困難，室內模擬試驗開始受到關注：孫美等進行了夾砂層狀土的渠道滲漏室內模擬試驗，發(fā)現(xiàn)夾砂層會使?jié)駶欎h出現(xiàn)不連續(xù)的現(xiàn)象，具有一定的阻水減滲效果[1]；張銳等開展了室內土箱模擬試驗，發(fā)現(xiàn)土壤體積質量越大，入滲水頭越高，灌溉需水量越少[2]；Zhang等[3]實施了不同土壤類型、初始含水量和壟溝尺寸下的6組室內試驗，觀測了累積入滲量、濕潤鋒分布的變化，為設計壟溝種植系統(tǒng)提供了參考；Phogat等[4]利用土箱試驗驗證了二維水平衡模型，并分析了渠床高程對滲流和地下水位上升的影響；冀榮華等[5]在室內模擬了負壓灌溉下的土壤水分運動，得到了灌水器半徑對土壤入滲有顯著影響的結論；聶衛(wèi)波等[6]通過模擬溝灌入滲，建立了溝灌土分運動數值模型，并以此為基礎，推導出溝灌條計算累積入滲量的方法。數值模擬方面，有限元計算機模型HYDRUS具有豐富的功能[7-12]：Bufon等[13]用HYDRUS-2D計算了滴灌條件的棉花含水量，發(fā)現(xiàn)計算結果與測量數據擬合程度較高；Kandelous等[14]用HYDRUS-2D模擬了地下滴灌系統(tǒng)的室內和大田試驗，并用均方根誤差對試驗結果和模擬結果進行了評價，認為其對應關系非常好；張勇勇[15]建立了壟溝灌溉下的土壤水動力學模型，用HYDRUS-2D軟件進行求解，得到了不同試驗條件下的土壤含水量、累積入滲量和濕潤鋒運移距離并驗證了設立模型的合理性。以上研究成果表明，用HYDRUS-2D模擬土壤水分入滲過程具有較高的可靠性，但并未考慮渠道邊坡系數對土壤入滲的影響及各因素之間的交互特性；邊坡系數作為渠道的重要參數，其變化對渠道入滲的影響非常重要，而各因素之間的交互影響對渠道斷面的優(yōu)化設計也起著決定性的作用。本文開展渠道滲漏室內模擬試驗，建立試驗條件下的土壤水分運動方程，通過試驗和模擬方法分析邊坡系數、入滲水頭和渠道底寬及其交互作用對濕潤鋒運移距離和累積入滲量的影響，研究渠道土壤水分運動的主要影響因素，以期為減少渠道水滲漏，提高灌溉水利用效率提供理論支持。

1 材料與方法

1.1 室內試驗

1.1.1 試驗設備

試驗設備由模擬土箱和供水設備2部分構成，模擬土箱規(guī)格為60 cm×20 cm×110 cm，由有機玻璃制成，為了防止氣阻，土箱側面和正面均留有排氣孔；為了讀數方便，土箱后側玻璃板貼有刻度條。試驗由馬氏瓶供水，馬氏瓶直徑20 cm，高80 cm[16-17]。試驗土箱內按照試驗設計要求，設置恰當的渠道斷面尺寸，試驗裝置如圖1。

圖1 室內模擬試驗裝置Fig. 1 Indoor simulation experimental setup

1.1.2 供試土樣

試驗土樣于2016年7月取自哈爾濱市東北農業(yè)大學水利學院試驗田（126°45′32″E、45°44′41″N）[18]。取樣深度為10～30 cm，土壤質地為壤土，顆粒比較均勻，試驗前，土樣風干破碎，過2 mm篩，配置初始質量含水率7%，按設計干容重1.4 g/cm3分層（10 cm）裝入試驗土箱中，層間打毛[19]。

1.1.3 試驗設計

室內試驗于2016年夏季在東北農業(yè)大學水利學院水工廳內進行。試驗采用三因素完全試驗設計，設計的渠道底寬分別為5、10、15 cm；入滲水頭分別為3、6、9 cm；邊坡系數分別為1、1.2、1.5，共27個處理，各3次重復，需進行81次試驗?？紤]到工作量較大，而且試驗目的是為了模型合理性分析，實際進行了15個處理試驗。每次試驗持續(xù)時間為480 min。

試驗開始時首先按照設計的渠道斷形狀裝填土樣，隨后調節(jié)馬氏瓶進氣管低端至設計入滲水頭高度，快速向土箱內加水至設計水位，同時打開其止水夾，再對進氣管低端高度進行微調。試驗開始后每隔5 min記錄馬氏瓶讀數，每隔10 min在土箱貼有刻度條的一側描繪濕潤峰曲線；試驗進行一段時間后根據濕潤鋒運移距離和入滲速度適當延長讀數、描繪時間間隔。試驗結束后用分辨率為1334×750像素的相機拍下繪制的濕潤峰曲線，整理數據。

1.2 模型建立

1.2.1 渠道水分運動模型

設計試驗條件下的渠道水分入滲過程可以簡化為二維非飽和水平及垂直入滲土壤水分運動問題，采用Richards方程進行描述[20]。

式中θ為土壤體積含水率，cm3/cm3；t為入滲時間，min；h為基質勢，cm；x、z分別水平坐標和垂直距離，cm；K(h)為土壤非飽和導水率，cm/min。

上述Richards方程中提及的非飽和導水率、土壤基質勢及含水率的關系采用VG-M模型進行描述[21-22]。

式中θr為殘余含水率，cm3/cm3；θs為飽和含水率，cm3/；KS為飽和導水率，cm/min；m、n、α為經驗系數，m=1-1/n。

根據離心機實測的土壤水分特征曲線，環(huán)刀法實測的飽和導水率等數據，采用MATLAB2012a擬合得到相關模型參數[23-26]：飽和含水率0.438 cm3/cm3；殘余含水率0.042 cm3/cm3；α為 0.029；n為1.274；飽和導水率0.049 cm3/cm3。

1.2.2 定解條件

所建模型的初始和邊界條件如圖2所示。

圖2 模型求解的邊界條件Fig. 2 Boundary conditions for model solving

初始條件：在設計的試驗條件下，土壤水分剖面為穩(wěn)定剖面。因此，模型的初始條件為計算區(qū)域內各點具有相同的水土勢。

式中H0為初始水土勢，cm；Ω為計算區(qū)域（即圖2中的陰影部分）；模擬區(qū)域最低點的z=0。

邊界條件：試驗過程中EFQ保持恒定入滲水頭h0(cm)；上邊界QG和GC覆蓋塑料薄膜，可以忽略蒸發(fā)，視為零通量邊界；左邊界EA為對稱軸，水平通量為0；右邊界CB視為零通量面；下邊界AB設置了排水孔，而且在試驗過程中濕潤鋒未到達土槽底部，不影響土壤水分入滲過程，保持恒定初始條件，設為自由排水邊界。

綜上，其邊界條件可表述為

式中t0為試驗延續(xù)時間，min。

1.2.3 數值計算

為了求解上述偏微分方程，采用HYDRUS-2D軟件對試驗設計方案下的濕潤鋒運移水平、豎直距離和累積入滲量進行數值模擬[27-29]，并將所得結果與試驗結果進行對比分析，采用均方根誤差（root mean square error，RMSE）和決定系數（R2）[30]檢驗所建立的土壤水分運動方程的合理性。

2 結果與分析

2.1 模型驗證

2.1.1 累積入滲量對比

從15組試驗中隨機選取9組試驗組合研究其累積入滲量實測值和模擬值動態(tài)變化（圖3）。結果表明，模擬值與實測值變化趨勢一致，均表現(xiàn)出初期入滲較快，隨入滲時間的延長入滲速率逐漸降低，并趨向穩(wěn)定的狀態(tài)。各次試驗累積入滲量的模擬值和實測值的RMSE和 R2如表1所示。由表1可知，RMSE的值在0.43～4.99 L之間變化，R2值在0.91～0.97之間變化（均高于0.9），說明實測值與模擬值之間的對應關系較好，所建立的土壤水分運動方程合理。

表1 累積入滲量實測值和模擬值的統(tǒng)計分析Table 1 Statistical analysis of measured and simulated cumulative infiltration

圖3 各試驗條件下累積入滲量的模擬值與實測值對比Fig. 3 Comparison of simulated and measured cumulative infiltration under various experimental conditions

2.1.2 濕潤鋒運移距離對比

以4種組合試驗為例（圖4），對比480 min內濕潤鋒運移距離模擬值和實測值，結果表明，各處理條件下濕潤鋒運移距離的實測值與模擬值基本吻合，均表現(xiàn)出初期運移速度較快，隨試驗時間的延長漸漸減緩，并趨于穩(wěn)定。

綜上，實測數據與模擬結果吻合較好，可以認為本文建立的土壤水動力學模型是可靠的，選取的相關參數值是適當的，用于模擬試驗條件下的累積入滲量和濕潤鋒運移過程是可行的。

2.2 基于HYDRUS模擬的累積入滲量影響因素分析

為了分析渠道底寬w、入滲水頭h0、邊坡系數f對累積入滲量的影響，用HYDRUS軟件模擬了27組渠道滲漏試驗。以圖5為例表明累積入滲量隨三者的變化。由圖可知，累積入滲量隨渠道底寬和入滲水頭變化較大，但邊坡系數f對累積入滲量的影響很小。

圖4 各試驗條件下濕潤峰的模擬值與實測值對比Fig. 4 Comparison of simulated and measured wetting front under various experimental conditions

圖5 各試驗條件下的累積入滲量對比Fig. 5 Comparison of cumulative infiltration under different experimental conditions

對w和h0下渠道480 min累積入滲量做雙因素方差分析，結果見表2。由表2可以看出，w對累積入滲量影響的F值為59.46，遠大于臨界值4.78，說明渠道底寬對入滲量的影響極顯著（P＜0.01）。當f=1、h0=3 cm時，若w由5 cm增加到15 cm時，累積入滲量由7 688 mL增加到37 027 mL，增加了382%。表2還表明h0對累積入滲量影響的F值為2.612，小于臨界值4.778，說明其對累積入滲量的影響不顯著（P＞0.05）。二者之間交互影響F值為0.30，說明w和h0的交互作用不顯著（P＞0.05）。這是因為，隨著f、h0和w的增加，濕周增加，入滲界面受水面積增加，因而累積入滲量增加。而f對濕周的影響很小，因此其對渠道入滲的影響不作考慮。且由表2可以看出，濕周隨w的增幅比隨h0的增幅大，因此w對累積入滲量的影響比h0顯著。

2.3 基于HYDRUS模擬的濕潤鋒運移距離影響因素分析

不同水頭下底寬、邊坡系數對濕潤鋒運移影響類似，故以水頭為9 cm表明不同底寬和邊坡系數對濕潤鋒運移的影響，如圖6所示?？梢钥闯?，邊坡系數對濕潤鋒運移距離幾乎沒有影響，而渠道底寬、入滲水頭與濕潤鋒運移距離呈正相關。

方差分析表明（表2），w對水平距離影響的F值為59.55，垂直距離為7.42，均大于臨界值4.78，而且水平距離的F值遠大于垂直距離，說明w對濕潤鋒運移距離的影響顯著（P＜0.01），且對水平距離的影響特別顯著（P＜0.01），比如在t=480 min，f=1、h0=3 cm的條件下，當w由5 cm增加到15 cm時，濕潤鋒運移水平、垂直距離分別增加了55.89%和47.82%。這是因為在f和h0不變的情況下，隨著w的增加，入滲界面濕周增加，濕潤鋒運移距離有所增加，而w的增加直接導致垂向入滲面積增加，因而對水平距離的影響更為顯著。入滲水頭對水平距離影響的F值為19.14，垂直距離為3.0。說明h0對水平距離的影響顯著（P＜0.05），如在w=5 cm的條件下，

當h0由3 cm增加到9 cm時，濕潤鋒運移水平增加了33.45%。在w和f不變的條件下，隨著h0增加，入滲界面壓力勢和受水面積增加，而h0的增加對垂向入滲面的影響較大，因而對水平距離的影響比較顯著。w、h0對濕潤鋒運移水平、垂直距離的交互影響F值分別為0.15和0.004，說明w、h0對濕潤鋒運移距離不存在顯著的交互作用（P＜0.05）。

圖6 各試驗條件下的濕潤鋒運移距離對比Fig. 6 Comparison of wetting front migration distance under different experimental conditions

2.4 入滲速率變化

所有試驗條件下的入滲速率均展現(xiàn)出隨時間的延長漸漸減小并趨向穩(wěn)定的征象；在試驗開始初期，入滲速率下降得較快，隨著時間的延長，土壤水土勢梯度漸漸減小，入滲速率也隨之減小并趨向一個定值。以f =1.5為例（圖7），在w和f相同的條件下，隨h0的增加，入滲速率增加。其原因為，h0增加時，入滲界面處壓力勢增大導致入滲速率增加；隨入滲時間的推移，土壤濕潤區(qū)域變大，水土勢不斷減小，并趨于穩(wěn)定，因此入滲速率漸漸趨于穩(wěn)定。在相同h0和f條件下，入滲速率隨渠道底寬的增加而增加，且增幅較大。其原因為，渠道底寬增加，入滲界面濕周增加，受水面積變大，使得入滲速率增大。在w和h0相同的條件下，由于f對入滲斷面濕周的影響很小，因此f對入滲速率幾乎沒有影響。

圖7 邊坡系數為1.5時不同底寬和入滲水頭下的入滲率對比Fig. 7 Comparison of infiltration rates under different canal base width and infiltration head with slope coefficient of 1.5

3 結 論

選用基于Richards方程的有限元模擬軟件HYDRUS-2D進行數值模擬得到的模擬值和實測值的的擬合程度較高（R2均高于0.9），可以契合地模擬試驗條件下的渠道滲漏土壤水分入滲規(guī)律。

模擬結果均表明，邊坡系數對土壤水分入滲過程的影響較為微弱；累積入滲量、濕潤鋒運移距離及入滲速率均隨入滲水頭增加而增加，且入滲水頭對濕潤鋒運移水平距離的影響顯著（P＜0.05）；隨渠道底寬的增加，濕潤鋒運移距離、累積入滲量以及入滲率均表現(xiàn)出增加的趨勢，且渠道底寬增加對累積入滲量和運移距離的影響極顯著（P＜0.01）。當邊坡系數為1、入滲水頭為3 cm時，若底寬由5 cm增加到15 cm時，480 min累積入滲量增加了382%，濕潤鋒運移水平、垂直距離分別增加了55.89%和47.82%。通過雙因素方差分析，發(fā)現(xiàn)渠道底寬和入滲水頭之間的沒有顯著的交互作用。

本文只針對各影響因素對渠道滲漏的影響做了定性分析，如何確定各試驗條件下累積入滲量和濕潤鋒運移距離隨影響因素變化的定量關系，尋求最優(yōu)經濟效益下的渠道斷面設計，將作為下一步研究的重點。

[1] 孫美，毛曉敏，陳劍，等. 夾砂層狀土條件下渠道滲漏的室內試驗和數值模擬[J]. 農業(yè)工程學報，2010，26(8)：33-38. Sun Mei, Mao Xiaomin, Chan Jian, et al. Laboratory experiment and simulation on canal seepage through sand interlayer[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(8): 33-38. (in Chinese with English abstract)

[2] 張銳，閆超，劉博，等. 入滲水頭及土壤體積質量對壟作溝灌土壤水分入滲和分布的影響[J]. 灌溉排水學報，2016，35(6)：37-41. Zhang Rui, Yan Chao, Liu Bo, et al. Effects of infiltration heads and soil volume-quality on soil water infiltratrion and distribution under ridge tillage furrow irrigation[J]. Journal of Irrigation and Drainage, 2016, 35(6): 37-41. (in Chinese with English abstract)

[3] Zhang Yongyong, Wu Pute, Zhao Xining, et al. Simulation of soil water dynamics for uncropped ridges and furrows under irrigation conditions[J]. Canadian Journal of Soil Science, 2013, 93(1): 85-98.

[4] Phogat V, Malik R S, Kumar S. Modelling the effect of canal bed elevation on seepage and water table rise in a sand box filled with loamy soil[J]. Irrigation Science, 2009, 27(3): 191 -200.

[5] 冀榮華，王婷婷，祁力鈞，等. 基于HYDRUS-2D的負壓灌溉土壤水分入滲數值模擬[J]. 農業(yè)機械學報，2015，46(4)：113-119. Ji Ronghua, Wang Tingting, Qi Lijun, et al. Numerical simulation of soil moisture infiltration under negative pressure irrigation based on HYDRUS-2D[J]. Transactions of the Chinese Society of Agricultural Machinery, 2015, 46(4): 113-119. (in Chinese with English abstract)

[6] 聶衛(wèi)波，馬孝義，王術禮. 溝灌土壤水分運動數值模擬與入滲模型[J]. 水科學進展，2009，20(5)：668-676. Nie Weibo, Ma Xiaoyi, Wang Shuli. Infiltration model andnumerical simulation of the soil water movement in furrow irrigation[J]. Advances in Water Science, 2009, 20(5): 668-676. (in Chinese with English abstract)

[7] Deb S K, Shukla M K, Mexal J G, et al. Evaluation of spatial and temporal root water uptake patterns of a flood-irrigated pecan tree using the HYDRUS (2D/3D) model[J]. Journal of Irrigation and Drainage Engineering, 2013, 139(8): 599-611.

[8] Han Ming, Zhao Chengyi, Feng Gary, et al. Evaluating the effects of mulch and irrigation amount on soil water distribution and root zone water balance using HYDRUS-2D[J]. Water, 2015, 7(6): 2622-2640.

[9] Saso J K, Parkin G W, Drury C F, et al. Chloride leaching in two Ontario soils: Measurement and prediction using HYDRUS-1D[J]. Canadian Journal of Soil Science, 2016, 92(2): 285-296.

[10] 李耀剛，王文娥，胡笑濤. 基于HYDRUS-3D的涌泉根灌土壤入滲數值模擬[J]. 排灌機械工程學報，2013，31(6)：546-552. Li Yaogang, Wang Wen′e, Hu Xiaotao. Numerical simulation of soil water infiltration under bubbled root irrigation based on HYDRUS-3D[J]. Journal of Drainage and Irrigation Machinery Engineering, 2013, 31(6): 546-552. (in Chinese with English abstract)

[11] 潘紅霞，付恒陽，賀屹. 基于HYDRUS-2D的地下滴灌下水分運移數值模擬研究[J]. 灌溉排水學報，2015，34(3)：52—55. Pan Hongxia, Fu Hengyang, He Yi. Numerica simulation of water movement under subsurface drip irrigation using HYDRUS-2D[J]. Journal of Irrigation and Drainage, 2015, 34(3): 52-55. (in Chinese with English abstract)

[12] El-Nesr M N, Alazba A A, ?im?nek J. HYDRUS simulations of the effects of dual-drip subsurface irrigation and a physical barrier on water movement and solute transport in soils[J]. Irrigation Science, 2014, 32(2): 111-125.

[13] Bufon V B, Lascano R J, Bednarz C, et al. Soil water content on drip irrigated cotton: Comparison of measured and simulated values obtained with the Hydrus 2D model[J]. Irrigation Science, 2012, 30(4): 259-273.

[14] 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): 1070-1076.

[15] 張勇勇. 壟溝灌溉土壤水分入滲模擬研究[D]. 北京：中國科學院大學，2013. Zhang Yongyong. Simulation Research on Soil Water Infiltration under Ridge-furrow Irrigation[D]. Beijing: University of Chinese Academy of Science. (in Chinese with English abstract)

[16] 張銳. 壟作溝灌土壤水分入滲規(guī)律的試驗研究[D]. 哈爾濱：東北農業(yè)大學，2016. Zhang Rui. Research on the Infiltration Law of Ridge Tillage Furrow Irrigation Soil Water[D]. Harbin: Northeast Agricultural University. (in Chinese with English abstract)

[17] 孫美，毛曉敏，陳劍. 渠道滲漏室內試驗土壤水分運動數值模擬[C]. 紀念中國農業(yè)工程學會成立30周年暨中國農業(yè)工程學會2009年學術年會，2009.

[18] 付強，侯仁杰，王子龍，等. 凍融期積雪覆蓋下土壤水熱交互效應[J]. 農業(yè)工程學報2015，31(15)：101-107. Fu Qiang, Hou Renjie, Wang Zilong, et al. Soil moisture thermal interaction effects under snow cover during freezing and thawing period[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(15): 101-107. (in Chinese with English abstract)

[19] 聶衛(wèi)波. 畦溝灌溉水流運動模型與數值模擬研究[D]. 楊凌：西北農林科技大學，2009. Nie Weibo. Research on Water Flow Model and Numerical Simulation for Border and Furrow Irrigation[D]. Yangling: Northwest Agriculture and Forestry University. (in Chinese with English abstract)

[20] 雷志棟，楊詩秀，謝森傳. 土壤水動力學[M]. 北京：北京清華大學出版社，1988.

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

[22] Mualem Y. A new model for predicting the hydraulic conductivity of unsaturated porous media[J]. Water Resources Research, 1976, 12(3): 513-522.

[23] 彭建平，邵愛軍. 用Matlab確定土壤水分特征曲線參數[J].土壤，2007，39(3)：433-438. Peng Jianping, Shao Aijun. Determination of parameters of soil water characteristic carve by Matlab[J]. Soils, 2007, 39(3): 433-438. (in Chinese with English abstract)

[24] 楊改強，霍麗娟，楊國義，等. 利用MATLAB擬合vanGenuchten方程參數的研究[J]. 土壤，2010，42(2)：268 -274. Yang Gaiqiang, Huo Lijuan, Yang Guoyi, et al. Research on fitting van Genuchten equation parameter with MATLAB Software[J]. Soils, 2010, 42(2): 268-274. (in Chinese with English abstract)

[25] 付強，蔣睿奇，王子龍，等. 基于改進螢火蟲算法的土壤水分特征曲線參數優(yōu)化[J]. 農業(yè)工程學報，2015，31(11)：117-122. Fu Qiang, Jiang Ruiqi, Wang Zilong, et al. Optimization of soil water characteristic curves parameters by modified firefly algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(11): 117-122. (in Chinese with English abstract)

[26] 劉洪波，張江輝，虎膽·吐馬爾白，等. 土壤水分特征曲線VG模型參數求解對比分析[J]. 新疆農業(yè)大學學報，2011，34(5)：437-441. Liu Hongbo, Zhang Jianghui, Hudan Tuma'erbai, et al. Contrast analysis on soil water characteristic curve of VG model parameters[J]. Journal of Xinjiang Agricultural University, 2011, 34(5): 437-441. (in Chinese with English abstract)

[27] Iqbal S, Guber A K, Khan H Z. Estimating nitrogen leaching losses after compost application in furrow irrigated soils of Pakistan using HYDRUS-2D software[J]. Agricultural Water Management, 2016, 168: 85-95.

[28] 李睿冉. HYDRUS-2D模型在渠道滲漏數值模擬中的應用[J]. 節(jié)水灌溉，2012(11)：71-74.

[29] 毛曉敏，姚立強，馮紹元，等. 層狀土條件下混凝土襯砌渠道滲漏及土壤水分分布的數值模擬[J]. 水利學報，2011，42(8)：949-955. Mao Xiaomin, Yao Liqiang, Feng Shaoyuan, et al. Numerical simulation on canal seepage and soil water distribution for concrete lining canal with layered soil structure[J]. Journal of Hydraulic Engineering, 2011, 42(8): 949-955. (in Chinese with English abstract)

[30] Singh D K, Rajput T B S, Singh D K, et al. Simulation of soil wetting pattern with subsurface drip irrigation from line source[J]. Agricultural Water Management, 2006, 83(12): 130-134.

HYDRUS simulation and verification of canal leakage and its influencing factors analysis

Fu Qiang, Li Yue, Li Tianxiao, Cui Song, Liu Dong
(School of Water Conservancy & Civil Engineering, Northeast Agricultural University, Harbin 150030, China)

Canal leakage is the main factor affecting the efficiency of irrigation water. In this paper, the main factors of canal leakage (referring to infiltration) were analyzed and numerical analysis was carried out. In order to obtain the main factors that affected canal leakage, the indoor leakage simulation experiment was carried out. The experimental site was located in Northeast Agricultural University, Harbin(126°45′32″E、45°44′41″N). The test case was made of Plexiglas material, with the size 60 cm ×20 cm×110 cm (length×width×height). Mariotte flask was used to ensure water supply and its diameter was 20 cm and the height was 80 cm. A total of 15 treatments with 3 replicates were conducted. During the experiment, the cumulative infiltration and wet front migration distance of each experiment were recorded. In addition, the soil water characteristic curve of the soil samples was tested and the relevant parameters in the VG model were fitted. The parameters of the fitting and the saturated hydraulic conductivity measured by the ring knife method were input into the HYDRUS-2D software to solve the soil water movement model under the experimental conditions, and the simulate values were obtained. The error analysis of the simulated and measured values was carried out. The root mean square error and the determination coefficient of each experiment were calculated. The results showed that the model was reasonable with the root mean square error of 0.43-4.99 L and the R2higher than 0.9, and the model could simulate the soil water movement process under the experimental conditions. Based on the measured and simulated values, the influence of the 3 factors on the canal leakage and the influence of these factors on the infiltration process were analyzed. The influence of the slope coefficient on the infiltration process was weak and its influence on the infiltration process could be ignored. Two-way analysis of variance was performed in order to analyze the influence of canal base width and infiltration head and its interaction on cumulative infiltration and wetting front migration distance. The results showed that the statistic value of the canal base width was 55.62 for the cumulative infiltration, 59.46 for the horizontal distance of wet front migration, and 7.42 for the vertical distance. All the 3 statistic values were greater than the critical value of 4.78 under the significance level of 0.01. It indicated that the canal base width had a significant impact on the cumulative infiltration and the wetting front migration distance. The statistic value of the infiltration head at the cumulative infiltration was 2.612, the statistic value of the horizontal distance of wet front migration was 19.14, and the vertical distance was 3.00. The infiltration head had a significant effect on the horizontal distance. The statistic values of the effects of canal width and infiltration head on the cumulative infiltration and the horizontal and vertical distance were 0.30, 0.15 and 0.004, respectively. They were less than the critical value of 3.47, indicating that there was no significant interaction between the canal width and infiltration head (P＞0.05). The results of the study can provide technical support for improving the utilization efficiency of irrigation water from the aspect of canal seepage control.

canals; infiltration; numerical analysis; HYDRUS-2D; wetting front; leakage

10.11975/j.issn.1002-6819.2017.16.015

TV732.6; S152.7+2

A

1002-6819(2017)-16-0112-07

付 強，李 玥，李天霄，崔 嵩，劉 東. 渠道滲漏HYDRUS模擬驗證及影響因素分析[J]. 農業(yè)工程學報，2017，33(16)：112-118.

10.11975/j.issn.1002-6819.2017.16.015 http://www.tcsae.org

Fu Qiang, Li Yue, Li Tianxiao, Cui Song, Liu Dong. HYDRUS simulation and verification of canal leakage and its influencing factors analysis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(16): 112-118. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.16.015 http://www.tcsae.org

2017-02-24

2017-06-10

國家自然科學基金（51479032、51579044）；國家重點研發(fā)計劃（2017YFC0406002）

付 強，男，遼寧錦州人，教授，博士生導師，主要從事農業(yè)水土資源系統(tǒng)分析、凍融土壤水熱作用機理等方面研究。哈爾濱 東北農業(yè)大學水利與土木工程學院，150030。Email：fuqiang0629@126.com。