段鳳儀，朱成立，王策，陳于，馬靚

灌溉優(yōu)先流對(duì)裂隙分布密度、灌水強(qiáng)度及土壤濕度的響應(yīng)機(jī)理

段鳳儀1，朱成立1，王策1*，陳于2，馬靚1

（1.河海大學(xué) 農(nóng)業(yè)科學(xué)與工程學(xué)院，南京 211100；2.江蘇省農(nóng)村水利科技發(fā)展中心，南京 210029）

【目的】研究土壤初始含水率（土壤濕度）、裂隙分布密度及灌水強(qiáng)度對(duì)裂隙優(yōu)先流強(qiáng)度的影響機(jī)理?！痉椒ā繕?gòu)建基于水量平衡原理的優(yōu)先流雙域滲透模型，并通過(guò)染色示蹤試驗(yàn)進(jìn)行物理模擬及模型驗(yàn)證，基于模型數(shù)值模擬了裂隙分布密度、初始含水率和灌水強(qiáng)度三維因素旋轉(zhuǎn)組合下的優(yōu)先流特征?！窘Y(jié)果】試驗(yàn)區(qū)的淺小型裂隙（裂隙平均寬度0.2 cm≤c≤1.0 cm，平均深度0.5 cm≤≤10 cm）與優(yōu)先流的發(fā)育無(wú)直接聯(lián)系，寬深型裂隙（裂隙平均寬度c≥2 cm，平均深度≥50 cm）是產(chǎn)生裂隙優(yōu)先流的主要原因；寬深型裂隙分布密度增大，基質(zhì)流深度、染色面積率以及灌水均勻度增加；模型對(duì)裂隙優(yōu)先流的模擬效果較好（2=0.951，＜0.01）；模擬結(jié)果顯示，裂隙深度很大程度上決定了優(yōu)先流的最大入滲深度；接近土壤飽和導(dǎo)水率的灌水強(qiáng)度無(wú)法觸發(fā)優(yōu)先流，隨著灌水強(qiáng)度增加，優(yōu)先流程度和最大入滲深度增加，基質(zhì)流深度和灌水均勻度降低；隨著土壤初始含水率和裂隙分布密度增大，基質(zhì)流深度和灌水均勻度增加，優(yōu)先流程度降低?！窘Y(jié)論】對(duì)于縮脹性較強(qiáng)、有機(jī)質(zhì)量豐富的農(nóng)田土壤，降低灌水強(qiáng)度、提高土壤初始含水率及裂隙分布密度能有效抑制裂隙優(yōu)先流并提高灌水質(zhì)量。

裂隙優(yōu)先流；雙域滲透模型；裂隙分布密度；土壤初始含水率；灌水強(qiáng)度

0 引言

【研究意義】土壤的空間異質(zhì)性導(dǎo)致水分繞過(guò)土壤基質(zhì)并沿著土壤中相互連通甚至貫穿的優(yōu)先通道集中流動(dòng)[1]，這種現(xiàn)象被稱為優(yōu)先流。農(nóng)田土壤干縮過(guò)程中產(chǎn)生的裂隙是誘發(fā)優(yōu)先流的常見因素[2]。由于干縮裂隙具有孔隙尺寸大、連通性強(qiáng)等特點(diǎn)，在高強(qiáng)度灌水或降水過(guò)程中，水分會(huì)沿裂隙通道向深層快速遷移，使灌溉初期優(yōu)先流成倍增加[3]。同時(shí)，入滲水流挾帶土壤中的養(yǎng)分或污染物沿著優(yōu)先通道直接并快速地下滲，造成土壤養(yǎng)分淋失和灌溉效率降低，并增加了污染地下水的風(fēng)險(xiǎn)[4-5]。因此，研究裂隙優(yōu)先流機(jī)理對(duì)降低水肥流失和地下水污染均具有實(shí)際意義。

【研究進(jìn)展】近年來(lái)，國(guó)內(nèi)外有關(guān)土壤優(yōu)先流的研究呈持續(xù)增長(zhǎng)的趨勢(shì)。Guo等[6]分析了190份有關(guān)優(yōu)先流的研究，發(fā)現(xiàn)大孔隙/裂隙通道特征、水分輸入特征（即灌水/降水特征）和初始土壤濕潤(rùn)條件是優(yōu)先流研究中較為關(guān)鍵的3個(gè)控制因素，但部分影響因素對(duì)裂隙優(yōu)先流的作用機(jī)制尚未存在定論。例如，Quisenberry等[7]通過(guò)田間試驗(yàn)發(fā)現(xiàn)，當(dāng)土壤初始含水率更高時(shí)，水分和溶質(zhì)遷移的深度更大，優(yōu)先流現(xiàn)象更明顯。而Hardie等[8]在農(nóng)場(chǎng)4個(gè)地點(diǎn)開展染色示蹤試驗(yàn)發(fā)現(xiàn)，與濕潤(rùn)土相比，干燥土中染色劑入滲深度和入滲率均成倍增加。其原因在于現(xiàn)場(chǎng)試驗(yàn)中初始含水率同時(shí)影響了土壤前期的導(dǎo)水特性及初始裂隙形態(tài)，復(fù)雜的土壤初始條件組合使得初始含水率對(duì)優(yōu)先流的影響機(jī)制難以得到準(zhǔn)確分析[9]。此外，Wells等[10]認(rèn)為，裂隙優(yōu)先流的入滲深度和入滲量與裂隙深度正相關(guān)。Luo等[11]研究發(fā)現(xiàn)，裂隙/大孔隙網(wǎng)絡(luò)連通性越好，土壤導(dǎo)水率越大，導(dǎo)水能力越強(qiáng)。而唐朝生等[12]研究表明，當(dāng)單位面積內(nèi)裂隙分布密度增大時(shí)，裂隙總長(zhǎng)度和網(wǎng)絡(luò)連通性增加，但裂隙的寬度和發(fā)育深度有所降低。裂隙特征的不同組合將影響優(yōu)先流的發(fā)育模式與強(qiáng)度。此外，降水/灌水強(qiáng)度對(duì)優(yōu)先流的激活效應(yīng)也受到了土壤初始狀態(tài)的強(qiáng)烈影響[13]?！厩腥朦c(diǎn)】國(guó)內(nèi)外關(guān)于灌水強(qiáng)度與土壤初始狀態(tài)對(duì)優(yōu)先流的研究中，主要以室外試驗(yàn)或模型研究為主，室外試驗(yàn)往往存在無(wú)關(guān)因素干擾，難以有效分析各因素對(duì)優(yōu)先流的影響機(jī)制；室內(nèi)試驗(yàn)或模型研究多為單因素分析，忽略了多因素組合的相互影響?！緮M解決的關(guān)鍵問(wèn)題】為此，本文基于裂隙優(yōu)先流雙域滲透模型，選取土壤初始含水率、裂隙分布密度及灌水強(qiáng)度3個(gè)關(guān)鍵因素開展多變量旋轉(zhuǎn)組合下的優(yōu)先流模擬。構(gòu)建模型基于表面入滲、層流運(yùn)移、裂隙邊壁吸滲與灌水強(qiáng)度相平衡的水量平衡原理[14]，并通過(guò)染色示蹤試驗(yàn)進(jìn)行驗(yàn)證。通過(guò)多變量組合模擬，建立各因素與裂隙優(yōu)先流的完整關(guān)系網(wǎng)絡(luò)，綜合分析各因素對(duì)裂隙優(yōu)先流的影響機(jī)制，為探索有效抑制裂隙優(yōu)先流的田間措施及灌水制度提供理論依據(jù)。

1 材料與方法

1.1 裂隙優(yōu)先流染色示蹤試驗(yàn)

1.1.1 供試土壤及其處理

試驗(yàn)區(qū)位于江蘇南京河海大學(xué)高效灌排與農(nóng)業(yè)水土環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室（31°86′N，118°60′E）。氣候?yàn)楸眮啛釒Ъ撅L(fēng)氣候，年均氣溫15.7 ℃，年均日照時(shí)間2 017.2 h，多年平均降水量1 200 mm，年均蒸發(fā)量900 mm。試驗(yàn)土壤取自試驗(yàn)區(qū)稻麥輪作旱田表層的黏壤土，將其過(guò)10 mm篩并人工去除根系雜草等雜質(zhì)后按1.25 g/cm3的體積質(zhì)量逐層（每層10 cm）回填至實(shí)驗(yàn)裝置內(nèi)，在每層回填后對(duì)土壤表面進(jìn)行磨毛處理，避免分層。實(shí)驗(yàn)裝置為用1.0 m×1.0 m×1.0 m的有機(jī)玻璃容器，其前側(cè)可拆卸，以便縱剖土壤，后側(cè)設(shè)置直徑1.0 cm、間隔10 cm的圓孔用于插入TDR（用于測(cè)量土壤含水率）的探針。填土完成后，采用模擬噴灌的方式使土壤濕潤(rùn)，隨后隔絕降水使其自然開裂。裂隙發(fā)育期間，定期觀測(cè)裂隙發(fā)育寬度與深度，待裂隙發(fā)育寬度與深度到達(dá)穩(wěn)定且土壤表層含水率降至凋萎系數(shù)后，裂隙培育結(jié)束。隨后利用相機(jī)記錄表面裂隙形態(tài)，采用彈性塑料桿測(cè)量裂縫深度。通過(guò)TDR測(cè)量土壤含水率，從而得到土壤初始含水率。

1.1.2 染色示蹤試驗(yàn)

裂隙培育完成后，采用碘化鉀-淀粉染色示蹤技術(shù)觀測(cè)該裂隙形態(tài)下的水流運(yùn)移。通過(guò)模擬降水器將含有碘化鉀示蹤劑（20 g/L）的灌溉水以60 mm/h[15]的灌水強(qiáng)度灌入開裂土壤中，灌水定額為100 mm[16]。待水分入滲重分布（12 h）后，每間隔4 cm對(duì)土壤進(jìn)行逐層縱剖，對(duì)每個(gè)剖面進(jìn)行平整，采用氣壓噴霧器對(duì)剖面均勻噴灑淀粉（50 g/L）和硝酸鐵（20 g/L）的混合液。由于灌溉水中加有碘化鉀示蹤劑，其中的碘離子會(huì)被三價(jià)鐵離子氧化為碘分子，碘分子與淀粉反應(yīng)將顯示出藍(lán)紫色，從而顯現(xiàn)出灌溉水的運(yùn)移區(qū)域。待20 min顯色穩(wěn)定后，對(duì)垂直剖面進(jìn)行拍照記錄。

1.1.3 染色圖像處理

將拍攝的剖面圖像通過(guò)Photoshop進(jìn)行扭曲矯正和剪裁，裁去底部和側(cè)邊受邊界效應(yīng)影響的區(qū)域。隨后利用Photoshop對(duì)染色區(qū)和未染色區(qū)進(jìn)行分割，并將染色區(qū)和未染色區(qū)域分別設(shè)置黑色和白色。最后通過(guò)MATLAB對(duì)圖像進(jìn)行二值化，將黑色和白色分別賦值為1和0，并將處理結(jié)果導(dǎo)出為數(shù)值化矩陣，用于后期計(jì)算。

1.1.4 裂隙優(yōu)先流評(píng)價(jià)參數(shù)

1）基質(zhì)流深度（cm），定義為土壤剖面表面至層間，每層染色率均不小于0.95 cm2/cm2時(shí)對(duì)應(yīng)的深度。層的染色率表達(dá)式為：

式中：為剖面像素圖中方向上的行數(shù)；dye()為剖面-Δ至+Δ層間的染色覆蓋率（cm2/cm2）；(,)為剖面圖上任意位置(,)處的賦值，未染色時(shí)(,)=0，染色時(shí)(,)=1；為剖面像素圖列數(shù)；為像素轉(zhuǎn)換為尺寸的轉(zhuǎn)換系數(shù)，這里取為1/20 cm/pxl。

2）優(yōu)先流指數(shù)pf（cm2/cm2），定義為優(yōu)先流區(qū)占總?cè)旧珔^(qū)dye面積比，計(jì)算式為：

式中：dye為總?cè)旧娣e（cm2），其值等于像素圖中染色像素總和乘以2；N為剖面像素圖對(duì)應(yīng)的列寬（cm），其值等于剖面像素圖列數(shù)乘以；其他指標(biāo)意義同上。

3）最大入滲深度M（cm），指染色溶液下滲的最大深度。剖面上任意點(diǎn)(,)的入滲深度1()可表示為：

式中：為剖面像素圖的行數(shù)；其他指標(biāo)意義同上。

4）灌水均勻度u（cm/cm），反映了染色區(qū)域分布的均勻度，是衡量灌水質(zhì)量的重要指標(biāo)，估算式為：

1.2 裂隙優(yōu)先流分析模型

1.2.1 模型原理

1）雙域滲透模型

灌溉水在重力作用下會(huì)沿優(yōu)先通道快速下滲產(chǎn)生優(yōu)先流，其流速遠(yuǎn)大于基質(zhì)流流速[17]。因此，可采用基質(zhì)流模式及快速遷移的優(yōu)先流模式分別對(duì)基質(zhì)域和裂隙域的水分運(yùn)移進(jìn)行模擬。本文采用基于水量平衡原理的裂隙優(yōu)先流雙域滲透模型[14]，基質(zhì)域的水分運(yùn)移仍視為毛細(xì)管驅(qū)動(dòng)流動(dòng)，采用Darcy-Richards 方程進(jìn)行描述；裂隙域中將優(yōu)先流簡(jiǎn)化為重力驅(qū)使的重力流，并假定為沿邊壁的層流運(yùn)動(dòng)，采用基于哈根-泊肅葉形式的層流方程模擬，同時(shí)考慮基質(zhì)域-裂隙域間水平吸滲，水平吸滲過(guò)程采用Green-Ampt模型進(jìn)行描述。由于層流參數(shù)復(fù)雜，本文對(duì)裂隙邊壁層流運(yùn)移過(guò)程進(jìn)行簡(jiǎn)化，通過(guò)水量平衡原理，即表面入滲、層流運(yùn)移、裂隙邊壁吸滲與灌水強(qiáng)度相平衡，迭代推求出層流運(yùn)移規(guī)律，以減少層流參數(shù)，提高模擬結(jié)果的準(zhǔn)確性。

2）水量平衡原理

將土壤收縮塊區(qū)簡(jiǎn)化為方形，裂隙網(wǎng)絡(luò)簡(jiǎn)化為縱橫交錯(cuò)且裂隙寬度與深度均勻的對(duì)稱網(wǎng)格。由于入滲曲線對(duì)稱分布，可將研究區(qū)對(duì)稱化（圖1）。

圖1 研究區(qū)概化模型與邊界條件

根據(jù)水量平衡原理，即灌水強(qiáng)度可轉(zhuǎn)換為土壤表面入滲率（Ⅰ）、裂隙邊壁水平吸滲率（Ⅱ）和裂隙邊壁層流通量（Ⅲ）3部分。因此，在任意時(shí)間段（，Δ）在水量平衡計(jì)算式為：

式中：sur為表層土壤入滲量（cm）；h為裂隙邊壁入滲量（cm）；Δlar為增加的層流量（cm2）；wf裂隙邊壁濕潤(rùn)段長(zhǎng)度（cm）；m為基質(zhì)域?qū)挾龋╟m）；0為灌水強(qiáng)度（cm/min）。

①表層土壤入滲量

土壤表層入滲量根據(jù)表層土壤入滲能力和灌水強(qiáng)度決定，假設(shè)灌水強(qiáng)度恒定且高于土壤表層飽和導(dǎo)水率s，即()=0，且0>s。若土壤表面初始狀態(tài)為完全干燥，則上表層入滲能力隨土壤含水率增高而逐漸減小。因此，不同階段的入滲率可采用Smith公式[18]描述為：

式中：0為經(jīng)驗(yàn)系數(shù)（min），其值等于漸近線橫坐標(biāo)；p為徑流發(fā)生時(shí)刻（min），即積水時(shí)刻或徑流發(fā)生時(shí)刻；和為入滲參數(shù)。

②裂隙邊壁入滲量

式中：h()為裂隙邊壁水平單寬入滲率（cm/min）。

當(dāng)0≤＜p時(shí)，模擬噴灌0全部通過(guò)土壤表面入滲，此時(shí)進(jìn)入裂隙的水流量為0，故裂隙邊壁入滲量也為0；當(dāng)≥p時(shí)，超出表層入滲能力的水分進(jìn)入裂隙內(nèi)，一部分沿裂隙邊壁層流下移，一部分沿裂隙邊壁水平入滲，裂隙邊壁入滲量由邊壁濕潤(rùn)段長(zhǎng)度wf和裂隙邊壁導(dǎo)水率決定，其裂隙邊壁濕潤(rùn)段長(zhǎng)度隨時(shí)間動(dòng)態(tài)變化，裂隙邊壁水平入滲率h()可根據(jù)Green-Ampt模型及土壤導(dǎo)水率函數(shù)獲得。

③裂隙內(nèi)增加的層流體積Δlar

Δlar指進(jìn)入裂隙的水流量，其層流厚度隨灌溉強(qiáng)度變化而變化，難以獲取，于是通過(guò)水量平衡原理將其簡(jiǎn)化，即先通過(guò)Hydrus 2D 軟件分段模擬，求得在時(shí)間步長(zhǎng)D（min）內(nèi)表層土壤入滲量sur以及裂隙邊壁入滲量h，后根據(jù)水量平衡原理建立等式，反復(fù)迭代推求出裂隙邊壁濕潤(rùn)段長(zhǎng)度wf的推移規(guī)律以及裂隙內(nèi)層流參數(shù)。其中裂隙濕潤(rùn)段長(zhǎng)度最大值不超過(guò)裂隙邊壁最大深度。

1.2.2 初始條件與邊界條件

根據(jù)模型原理可知，優(yōu)先流的發(fā)育可分為2個(gè)階段，各階段初始及邊界條件有所不同：

1）基質(zhì)流階段

當(dāng)0≤≤p時(shí)，水分沿土壤表面均勻入滲，水分上邊界（OA）為定通量邊界；左右邊界（OF、AC、DE）正負(fù)通量近似相等，下邊界（EF）通量可忽略不計(jì)，均設(shè)為零通量邊界，即：

2）裂隙優(yōu)先流階段

當(dāng)＞p時(shí)，土壤表面產(chǎn)生的徑流進(jìn)入裂隙內(nèi)，一部分沿裂隙邊壁層流下移，一部分水平吸滲至基質(zhì)域內(nèi)，形成優(yōu)先流。AB邊界在此階段為動(dòng)態(tài)邊界，于是有：

式中：()為土壤飽和導(dǎo)水率；為負(fù)壓水頭（cm）。

1.2.3 計(jì)算參數(shù)

土壤持水特征曲線和非飽和導(dǎo)水率分別基于VG模型和Mualem模型：

表1 試驗(yàn)土壤的物理-水力參數(shù)

1.2.4 模型驗(yàn)證方法

模型評(píng)價(jià)采用決定系數(shù)2和均方根誤差（root mean square of error，）驗(yàn)證模型有效性。一般認(rèn)為2越大，越小，模型擬合的效果越好。

2 結(jié)果與分析

2.1 試驗(yàn)結(jié)果與分析

從染色縱剖面中選取9個(gè)典型剖面進(jìn)行分析，由于裂隙隨機(jī)分布的特點(diǎn)，每個(gè)剖面可視為1組重復(fù)，土體表面裂隙網(wǎng)絡(luò)圖及典型剖面的染色分布如圖2所示，圖2（a）中S（=1, 2, …, 9）表示裂隙節(jié)點(diǎn)/端點(diǎn)。根據(jù)試驗(yàn)測(cè)定的裂隙形態(tài)數(shù)據(jù)，可將試驗(yàn)區(qū)裂隙分為3類：淺小型裂隙（裂隙平均寬度0.2 cm≤c≤1.0 cm，裂隙平均深度0.5 cm≤≤10 cm）、寬深型主裂隙（裂隙平均寬度c≥2 cm，裂隙平均深度≥50 cm）及介于二者之間的中等裂隙。圖中除3條呈“Y”字形連通分布的主裂隙（S1S4、S2S9、S3S6）以及1條最大深度達(dá)13 cm的中等裂隙（S7S8）外，其他裂隙均為表面淺小裂隙。由于淺小裂隙在灌水過(guò)程中會(huì)快速閉合，對(duì)水分入滲影響不明顯[19]，所以優(yōu)先流主要取決于寬深型裂隙。將圖2（a）和圖2（b）對(duì)應(yīng)分析可知，本試驗(yàn)的染色結(jié)果與該結(jié)論吻合，優(yōu)先流的產(chǎn)生位置均與主裂隙一一對(duì)應(yīng)：剖面1—剖面3的優(yōu)先流主要由主裂隙S5S6作用產(chǎn)生；剖面4—剖面6上呈一大一小的2部分優(yōu)先流分別由S1S2和S7S8共2條裂隙作用產(chǎn)生，而剖面7—剖面9經(jīng)過(guò)主裂隙 S3S5，并且鄰近裂隙S3S4、S2S9，其染色覆蓋率受到了S3S4、S2S9裂隙在縱軸方向上的水分水平運(yùn)移的影響，因此該優(yōu)先流分布模式受到以上3條主裂隙控制。

圖2 裂隙網(wǎng)絡(luò)及典型剖面染色分布

分析各剖面染色覆蓋率隨深度的變化曲線（圖3）可知，開裂土體的前（剖面1—剖面3）、中（剖面7—剖面9）、后（剖面4—剖面6）部的覆蓋率減小趨勢(shì)呈明顯區(qū)別。剖面1—剖面3的染色覆蓋率在土層深度9 cm左右迅速減小，而剖面7—剖面9的染色覆蓋率從土層深度20 cm左右才逐漸減小。計(jì)算各剖面優(yōu)先流評(píng)價(jià)參數(shù)（表2）可知，隨著主裂隙的增多，基質(zhì)流深度、染色面積率以及灌水均勻度顯著增加。由此可知，各剖面主裂隙數(shù)量的不同是造成染色分布模式出現(xiàn)差異的主要原因。主裂隙數(shù)的增加會(huì)增大表面裂隙率及裂隙邊壁面積，導(dǎo)致入滲總表面積增加，土壤表面的入滲能力增強(qiáng)，基質(zhì)流階段延長(zhǎng)，基質(zhì)流深度因此增大；此外，主裂隙數(shù)的增加導(dǎo)致水分沿多個(gè)優(yōu)先通道下滲，使得滲漏量成倍增加，因此染色面積率增加；最后，主裂隙數(shù)的增加及裂隙間連通性的增加，有助于水分同時(shí)沿多個(gè)優(yōu)先通道下滲，避免了局部集中下滲，因此均勻度提高。此外，圖3顯示各剖面最大入滲深度沒(méi)有顯著差異，均處于60～70 cm之間，說(shuō)明最大入滲深度與裂隙數(shù)量分布差異沒(méi)有直接聯(lián)系，主要取決于裂隙最大深度。綜上所述，裂隙優(yōu)先流的形成不僅受到了裂隙形態(tài)（裂隙深度及寬度）的影響，還與裂隙分布特征（剖面上的主裂隙數(shù)）密切相關(guān)。其中，裂隙的分布特征主要通過(guò)影響裂隙間的連通性及土壤表面的導(dǎo)水能力影響優(yōu)先流的分布。

圖3 典型剖面染色覆蓋率隨土層深度變化

表2 各剖面優(yōu)先流評(píng)價(jià)參數(shù)

2.2 模型模擬與驗(yàn)證

采用Hydrus 2D對(duì)優(yōu)先流進(jìn)行分階段模擬。基于剖面1的1/2對(duì)稱區(qū)域進(jìn)行建模，模擬區(qū)水平尺寸為 40 cm，垂直尺寸為70 cm，半裂隙寬度為1 cm，裂隙深度為50 cm，各層土壤的物理-水力參數(shù)見表1。模擬灌溉強(qiáng)度0=0.10 cm/min[15]，灌溉總時(shí)長(zhǎng)為100 min[16]。上邊界條件采用Smith公式描述，其中的臨界時(shí)間p（又叫優(yōu)先流發(fā)生時(shí)刻或徑流點(diǎn)）通過(guò)數(shù)值模擬得到，基于土壤剖面形態(tài)與水力參數(shù)采用Hydrus2D模擬表面入滲，當(dāng)表面含水率達(dá)到表層土壤飽和含水率時(shí)視為優(yōu)先流發(fā)生，模擬可得p=14.4 min。之后采用SPSS 25.0擬合計(jì)算Smith公式參數(shù)、和0，其值分別為0.554、0.828、4.298 min。據(jù)此得到土壤上邊界條件，由此計(jì)算出進(jìn)入裂隙的流量，并根據(jù)水量平衡原理推算裂隙內(nèi)的層流參數(shù)及水平吸滲參數(shù)。隨后分階段對(duì)入滲過(guò)程進(jìn)行模擬，最大迭代數(shù)設(shè)定為10，不考慮入滲過(guò)程中裂隙的愈合情況，最后模擬得到該剖面優(yōu)先流的發(fā)育過(guò)程如圖4所示。

圖4 數(shù)值模擬優(yōu)先流發(fā)育過(guò)程

由圖4可知，裂隙優(yōu)先流的發(fā)育過(guò)程可分為：①基質(zhì)流階段（0≤≤14.4 min），水分垂直入滲，濕潤(rùn)峰均勻向下推進(jìn)，呈基質(zhì)流態(tài)，該階段持續(xù)時(shí)間受灌水強(qiáng)度和初始含水率控制；②裂隙內(nèi)層流推移階段（14.4 min≤≤75 min），當(dāng)土壤表面達(dá)到飽和含水率后，多余的水分進(jìn)入裂隙，并通過(guò)層流的方式向下推進(jìn)，同時(shí)水分通過(guò)裂隙邊壁向四周擴(kuò)散，呈現(xiàn)出漏斗狀的優(yōu)先流態(tài)。層流推移速度主要受灌水強(qiáng)度控制，其最大推移深度受裂隙深度限制，因此裂隙深度很大程度上決定了優(yōu)先流的最大入滲深度。③裂隙飽和入滲階段（75 min≤≤100 min），當(dāng)層流運(yùn)動(dòng)到達(dá)裂隙底端，裂隙內(nèi)逐漸積水飽和，之后裂隙底及邊壁在積水作用下開始飽和入滲，因此優(yōu)先流最大入滲深度往往略大于裂隙最大深度。該階段優(yōu)先流發(fā)育主要受裂隙邊壁飽和入滲率和裂隙內(nèi)壓力水頭的影響，優(yōu)先流發(fā)育緩慢。

根據(jù)染色剖面與模擬圖像，計(jì)算出染色覆蓋率的實(shí)測(cè)值與模擬值如圖5所示。經(jīng)計(jì)算，2=0.951（＜0.01），為0.095，說(shuō)明該模型能較好地模擬優(yōu)先流的分布特征。但模擬值在15～30 cm土層深度處的染色覆蓋率略小于實(shí)測(cè)值，其原因在于實(shí)際噴灌過(guò)程中水分直接落入裂隙內(nèi)并沿裂隙壁發(fā)生水平入滲造成實(shí)測(cè)結(jié)果與模擬結(jié)果存在一定差異。此外，由于裂隙深度測(cè)量不準(zhǔn)確，測(cè)量深度小于實(shí)際深度，因此最大入滲深度的模擬值偏小。整體來(lái)看，基質(zhì)流模擬的準(zhǔn)確性高于優(yōu)先流模擬部分。

圖5 裂隙優(yōu)先流模擬值與實(shí)測(cè)值（剖面1）對(duì)比

2.3 不同初始條件組合下的入滲模擬

為進(jìn)一步揭示不同土壤初始含水率、裂隙分布密度和灌水強(qiáng)度對(duì)優(yōu)先流程度和灌水質(zhì)量的影響，對(duì)裂隙優(yōu)先流雙域滲透模型進(jìn)行應(yīng)用。本文采用單位面積內(nèi)收縮塊區(qū)數(shù)量（）反映裂隙分布密度。選取了2組初始含水率0（0.2、0.35 cm3/cm3）、3組裂隙分布密度（2×2、4×4、6×6 cm/cm2）以及4組灌水強(qiáng)度0（0.10、0.06、0.03、0.02 cm/min）進(jìn)行旋轉(zhuǎn)組合設(shè)計(jì)，各組處理灌水定額均為6.0 cm?；谇捌谘芯砍晒芍?，裂隙面積率及裂隙深度與土壤含水率負(fù)相關(guān)[20]，根據(jù)其相關(guān)曲線，得到含水率為0.2 cm3/cm3和0.35 cm3/cm3時(shí)所對(duì)應(yīng)的裂隙面積率為9.5%和4.5%。同時(shí)，根據(jù)相同含水率下裂隙表面的裂隙率相同以及裂縫的平均寬度隨裂隙分布密度增大而減小的規(guī)律[12]確定裂隙分布密度下的裂隙寬度。裂隙深度按試驗(yàn)測(cè)定的寬深比=0.04 cm/cm確定。數(shù)值模擬不同初始條件和邊界條件的設(shè)定見表3。

表3 數(shù)值模擬初始條件和邊界條件的旋轉(zhuǎn)組合設(shè)計(jì)

根據(jù)模擬結(jié)果（圖6）可知，灌水強(qiáng)度0、初始土壤含水率0以及裂隙分布密度均將影響優(yōu)先流的發(fā)育程度。橫向?qū)Ρ认嗤跏己始傲严斗植济芏鹊膬?yōu)先流可知，隨著灌水強(qiáng)度的增加，基質(zhì)流深度降低，優(yōu)先流程度顯著增加。其原因在于灌水強(qiáng)度的增加使土壤表面更快達(dá)到飽和含水率，加快灌溉水由表層入滲轉(zhuǎn)為沿裂隙優(yōu)先下滲，基質(zhì)流階段縮短，基質(zhì)流深度降低。此外，灌水強(qiáng)度越大，層流推進(jìn)速度越快，相同時(shí)間內(nèi)最大入滲深度增大。當(dāng)灌水強(qiáng)度較小（如灌水強(qiáng)度0=0.02 cm/min）且趨近土壤飽和導(dǎo)水率s=0.018 6 cm/min時(shí)，灌溉水幾乎全部通過(guò)土壤基質(zhì)吸滲進(jìn)入土壤中，無(wú)優(yōu)先流現(xiàn)象。因此，接近土壤飽和導(dǎo)水率的灌水強(qiáng)度無(wú)法觸發(fā)優(yōu)先流?？v向?qū)Ρ认嗤嗨畯?qiáng)度及裂隙分布密度下的優(yōu)先流可知，隨著土壤含水率升高，基質(zhì)流深度增加，優(yōu)先流程度顯著降低。其原因在于土壤初始含水率的增加導(dǎo)致初始裂隙形態(tài)的深度和寬度降低，由于優(yōu)先流最大入滲深度很大程度取決于裂隙深度，因此最大入滲深度顯著降低；其次，高含水率下土壤基質(zhì)導(dǎo)水能力的提高和儲(chǔ)水能力的降低[21]，使得濕潤(rùn)峰推移速度加快，因此土壤基質(zhì)流的深度增大、水平入滲距離增加。最后，縱向?qū)Ρ认嗤始肮嗨畯?qiáng)度的優(yōu)先流可知，隨著裂隙分布密度的增加，優(yōu)先流程度顯著減弱。其主要原因在于隨著裂隙條數(shù)和總長(zhǎng)度增加，裂縫的平均深度和寬度相應(yīng)減小[12]，因此優(yōu)先流的最大入滲深度顯著降低；而裂隙條數(shù)和總長(zhǎng)度的增加會(huì)增大裂隙邊壁的總表面積，這不僅會(huì)增大表層土壤的入滲能力，使得基質(zhì)流深度增加，而且會(huì)使裂隙域與基質(zhì)域的水分交換量的占比增加，導(dǎo)致沿裂隙邊壁下滲的層流厚度和最大入滲深度降低；此外，隨著裂隙分布密度和網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)的增加，裂隙間的水平連通度提高[22]，增加了水分的橫向交換，超過(guò)土壤表面入滲能力的灌溉水能沿多個(gè)裂隙通道下滲，避免了集中下滲，提高了入滲的均勻性。綜上所述，不是任何一個(gè)連通、延長(zhǎng)的裂隙都能產(chǎn)生優(yōu)先流，它僅提供了傳導(dǎo)優(yōu)先流的一個(gè)基礎(chǔ)條件，灌水強(qiáng)度才是決定優(yōu)先流能否觸發(fā)和觸發(fā)快慢的重要條件；而土壤初始含水率反映了一種環(huán)境條件，影響了觸發(fā)優(yōu)先流的容易程度。

圖6 等灌水定額下初始含水率-灌水強(qiáng)度-裂隙密度組合設(shè)計(jì)的優(yōu)先流數(shù)值模擬結(jié)果

圖7分析了不同土壤初始含水率、裂隙分布密度和灌水強(qiáng)度數(shù)值模擬下的優(yōu)先流評(píng)價(jià)指標(biāo)，包括優(yōu)先流指數(shù)、最大入滲深度、基質(zhì)流深度和灌水均勻度4個(gè)參數(shù)。據(jù)此建立了各因素與優(yōu)先流程度的關(guān)系網(wǎng)絡(luò)。

由圖7可知，隨著灌水強(qiáng)度增大，基質(zhì)流深度和灌水均勻度降低，優(yōu)先流指數(shù)增大；隨著初始含水率降低，基質(zhì)流深度減小，最大入滲深度和優(yōu)先流指數(shù)增大；隨著裂隙分布密度增加，基質(zhì)流深度和灌水均勻度增大，優(yōu)先流程度減弱。當(dāng)土壤初始含水率0=0.20 cm3/cm3、裂隙分布密度為2×2 cm/cm2、灌水強(qiáng)度0=0.10 cm/min時(shí)，優(yōu)先流指數(shù)達(dá)到最大值，基質(zhì)流深度和灌水均勻度最小。針對(duì)低初始含水率時(shí)，最大入滲深度與灌水強(qiáng)度正相關(guān)；而初始含水率較高（0=0.35 cm3/cm3）時(shí)，最大入滲深度與灌水強(qiáng)度負(fù)相關(guān)，歸因于高含水率下土壤的儲(chǔ)水量和裂隙深度大大降低，此時(shí)較大的灌水強(qiáng)度可能會(huì)超過(guò)土壤表層和裂隙邊壁最大入滲能力的總和，大量水分在模擬灌溉結(jié)束時(shí)并沒(méi)有進(jìn)入土體中，入滲深度因此反而降低。綜上所述，灌水強(qiáng)度主要通過(guò)影響表層土壤達(dá)到飽和的時(shí)間以及裂隙內(nèi)的層流推移速度，分別影響基質(zhì)流深度和最大入滲深度，從而影響優(yōu)先流程度；土壤初始含水率主要通過(guò)影響裂隙深度和土壤導(dǎo)水性質(zhì)，分別影響最大入滲深度和基質(zhì)流深度，由此影響優(yōu)先流程度；而裂隙分布密度主要通過(guò)影響裂隙深度、土壤的導(dǎo)水性質(zhì)以及裂隙網(wǎng)絡(luò)連通性，分別影響最大入滲深度、基質(zhì)流深度以及灌水均勻度，從而對(duì)優(yōu)先流程度產(chǎn)生影響。

圖7 土壤初始含水率-灌水強(qiáng)度-裂隙密度旋轉(zhuǎn)設(shè)計(jì)模擬結(jié)果的優(yōu)先流評(píng)價(jià)參數(shù)

3 討論

3.1 裂隙特征對(duì)優(yōu)先流的影響

本研究發(fā)現(xiàn)，淺小裂隙對(duì)增強(qiáng)優(yōu)先流無(wú)明顯作用；裂隙深度越深，水分最大入滲深度越大，優(yōu)先流程度越高。這一結(jié)論與Wells等[10]的試驗(yàn)結(jié)果吻合。其主要原因在于裂隙內(nèi)由重力勢(shì)主導(dǎo)的層流運(yùn)動(dòng)常高出基質(zhì)流速幾個(gè)數(shù)量級(jí)，在灌溉水有充足補(bǔ)充的情況下，水分會(huì)以極快的速度沿裂隙邊壁層流推移，且層流推移的最大深度取決于裂隙深度，因此裂隙越深，最大入滲深度越大，優(yōu)先流越明顯。但本試驗(yàn)及模型均采用均質(zhì)土，土壤原有孔隙結(jié)構(gòu)被破壞，實(shí)際情況下優(yōu)先流形成的深層濕潤(rùn)鋒可能會(huì)與土壤中其他類型的孔隙流道（蟲洞、根通道）相連接，優(yōu)先流更為顯著。因此，具有一定延伸深度或縱向連通性的裂隙/大孔隙網(wǎng)絡(luò)是優(yōu)先流傳導(dǎo)的前提條件，裂隙的延伸深度通過(guò)影響水分的最大入滲深度與總?cè)霛B量成為了影響優(yōu)先流程度的關(guān)鍵因素。而裂隙發(fā)育寬度會(huì)影響灌水過(guò)程中裂隙的愈合速度，是影響優(yōu)先流發(fā)育的另一重要因素。此外，優(yōu)先流程度與裂隙分布特征也有密切聯(lián)系，裂隙分布特征包括裂隙連通性及裂隙分布密度。對(duì)于一些具有較強(qiáng)的縱向連通性和垂直連續(xù)性的裂隙或孔隙網(wǎng)絡(luò)，即使在低強(qiáng)度降雨下，也會(huì)發(fā)生明顯的優(yōu)先流[23]。裂隙/大孔隙網(wǎng)絡(luò)的連通性越強(qiáng)，土壤導(dǎo)水能力越強(qiáng)[11]，因此，裂隙連通性一般會(huì)促進(jìn)優(yōu)先流發(fā)育。而本文數(shù)值模擬發(fā)現(xiàn)，隨著裂隙分布密度的增加，裂隙間的橫向連通性增強(qiáng)，避免了土壤水的局部集中下滲，提升了入滲的均勻性，優(yōu)先流程度減弱。因此，連通度對(duì)優(yōu)先流的影響效果主要取決于裂隙間連通度的增加是抑制還是促進(jìn)了水分的集中下滲。此外，根據(jù)土壤開裂規(guī)律及模擬結(jié)果發(fā)現(xiàn)，若土體收縮塊區(qū)減小，盡管裂隙分布密度增大，但裂隙發(fā)育的深度和寬度降低[12]，優(yōu)先流程度和最大入滲深度因此大幅降低，這為抑制土壤裂隙優(yōu)先流提供了新思路。綜上所述，裂隙網(wǎng)絡(luò)是優(yōu)先流產(chǎn)生的前提條件，裂隙的幾何形態(tài)、彎曲度以及網(wǎng)絡(luò)連通度等特征通過(guò)影響水流在土壤中流動(dòng)的路徑與速率，直接影響了優(yōu)先流的空間結(jié)構(gòu)與發(fā)育程度。因此，抑制裂隙優(yōu)先流的根本方法是抑制裂隙發(fā)育。除了通過(guò)增施有機(jī)肥、植被覆蓋等方法抑制土壤開裂外，在土壤中添加秸稈等植物殘?bào)w或纖維能通過(guò)改變土壤的塑性和破碎強(qiáng)度，縮小收縮塊區(qū)，從而有效抑制裂隙寬度和深度的開展[24]，但其對(duì)優(yōu)先流的抑制效果，有待在試驗(yàn)與生產(chǎn)實(shí)際中驗(yàn)證。

3.2 灌水特征對(duì)優(yōu)先流的影響

灌水或降水過(guò)程中，當(dāng)灌水強(qiáng)度超過(guò)周圍土壤基質(zhì)的入滲能力，或在土壤表層含水率達(dá)到一個(gè)閾值時(shí)，便會(huì)觸發(fā)土壤水的優(yōu)先流動(dòng)[1]。Edwards等[25]發(fā)現(xiàn)，灌水強(qiáng)度越大，土壤水入滲量越大，優(yōu)先流發(fā)生時(shí)間越短。本文數(shù)值模擬結(jié)果與該結(jié)果一致。其原因之一在于灌水強(qiáng)度影響了優(yōu)先流的觸發(fā)時(shí)間。當(dāng)灌水強(qiáng)度低于土壤表面的飽和導(dǎo)水率時(shí)，灌溉水全部以基質(zhì)流的形式均勻下滲，難以觸發(fā)優(yōu)先流，而較高的灌水強(qiáng)度會(huì)使土壤表面更快達(dá)到飽和，加快表層均勻入滲轉(zhuǎn)為沿裂隙優(yōu)先下滲，縮短優(yōu)先流觸發(fā)時(shí)間，增強(qiáng)了優(yōu)先流程度。其次，高強(qiáng)度灌溉提高了裂隙通道內(nèi)的水流通量，改變了土壤水分運(yùn)動(dòng)時(shí)以毛管勢(shì)能為主導(dǎo)的能量驅(qū)動(dòng)方式。當(dāng)灌水強(qiáng)度較低時(shí)，進(jìn)入優(yōu)先通道內(nèi)的水分較少，優(yōu)先流只能以薄膜的狀態(tài)沿著孔壁和團(tuán)聚體接觸點(diǎn)緩慢遷移[1]，此時(shí)毛管勢(shì)為主要的驅(qū)動(dòng)力，限制了水分快速流動(dòng)。只有當(dāng)水膜達(dá)到一定厚度（高達(dá)約100 μm）[26]，水膜才能在以重力勢(shì)能為主導(dǎo)的驅(qū)動(dòng)方式下實(shí)現(xiàn)快速流動(dòng)。較高的灌水強(qiáng)度能幫助土壤表層建立更多的正壓勢(shì)，提高進(jìn)入優(yōu)先通道的水流通量，增加水膜厚度，改變水分驅(qū)動(dòng)方式，推動(dòng)水膜實(shí)現(xiàn)連續(xù)快速的優(yōu)先流動(dòng)。因此，灌水強(qiáng)度是決定優(yōu)先流能否觸發(fā)和觸發(fā)快慢的重要條件。此外，Gjettermann等[27]研究表明，灌溉強(qiáng)度越高，在土壤表面35 cm以下區(qū)域，被染色的優(yōu)先流通道越多，這表明較高的灌溉強(qiáng)度會(huì)激活更多的優(yōu)先流通道，促進(jìn)優(yōu)先流的形成。此外，灌水持續(xù)時(shí)間、灌溉時(shí)間節(jié)點(diǎn)及灌溉方式等也是影響優(yōu)先流程度的重要外部條件[28-29]。因此，因地制宜地選擇灌溉制度或灌溉方式是抑制優(yōu)先流、提高灌水質(zhì)量的關(guān)鍵。例如，對(duì)于優(yōu)先通道延伸深度較大且地下水污染敏感性較強(qiáng)的地區(qū)，可通過(guò)低強(qiáng)度持續(xù)噴灌，以防止水分深層滲漏造成污染[30]。

3.3 土壤初始含水率對(duì)裂隙優(yōu)先流的影響

土壤初始含水率通過(guò)影響水分的入滲和傳導(dǎo)，成為影響優(yōu)先流發(fā)育的另一關(guān)鍵因素。本研究顯示，較低的土壤含水率會(huì)使優(yōu)先流更明顯。其原因在于，對(duì)于縮脹性較強(qiáng)的土壤，如本試驗(yàn)中的黏壤土，含水率的降低會(huì)加劇土壤裂隙的收縮[20]，裂隙深度和體積的增加將引起水分的大量滲漏[19]。此外，土壤的疏水性也是造成低含水率下優(yōu)先流加劇的關(guān)鍵原因[31]。疏水性常發(fā)生在沙質(zhì)土壤或有機(jī)質(zhì)量豐富的土壤中，與土壤濕度成反比[32]。因此，干燥土壤的表面常會(huì)因?yàn)槭杷源嬖谳^大的入水勢(shì)，導(dǎo)致灌溉水難以從土壤基質(zhì)入滲而被迫流向大孔或裂隙通道內(nèi)產(chǎn)生優(yōu)先流。由于孔隙通道內(nèi)的生物活動(dòng)和有機(jī)質(zhì)降解，很多裂隙或大孔隙內(nèi)壁也存在一定的疏水性[9]，這將阻止孔隙通道內(nèi)的水分與基質(zhì)水分交換，進(jìn)一步加劇優(yōu)先流的發(fā)育。此外，低含水率下的低導(dǎo)水率也是迫使水分進(jìn)入裂隙產(chǎn)生優(yōu)先流的原因之一。然而，較高的初始含水率會(huì)促進(jìn)優(yōu)先流的發(fā)育，這種情況通常發(fā)生在以壤土或粉土為主的穩(wěn)定大孔隙或生物孔隙（而不是收縮裂縫）中[8]。較高的土壤初始含水率降低了土壤基質(zhì)吸力和土壤基質(zhì)中的儲(chǔ)水空間，減少了從大孔隙到土壤基質(zhì)的側(cè)向滲透[33]，推動(dòng)水分向深層土壤中移動(dòng)，因此優(yōu)先流的總滲流量和穿透深度增加。此外，濕潤(rùn)土壤的孔隙活躍度更高，大孔間的連通性更強(qiáng)[34]，有助于建立更高效的流動(dòng)網(wǎng)絡(luò)，增強(qiáng)優(yōu)先流動(dòng)。綜上所述，由于土壤基質(zhì)的導(dǎo)水特性、孔隙特征及疏水性的強(qiáng)弱程度隨初始含水率的變化趨勢(shì)受土壤質(zhì)地結(jié)構(gòu)的影響較大，土壤初始含水率對(duì)優(yōu)先流的影響效果難以進(jìn)行普遍性歸納。

4 結(jié)論

1）淺小型裂隙（裂隙平均寬度0.2 cm≤c≤1.0 cm，平均深度0.5 cm≤≤10 cm）與優(yōu)先流發(fā)育沒(méi)有直接聯(lián)系，寬深型裂隙（裂隙平均寬度c≥2 cm，平均深度≥50 cm）是產(chǎn)生裂隙優(yōu)先流的主要原因；寬深型裂隙分布密度增大，基質(zhì)流深度、染色面積率以及灌水均勻度顯著增加。

2）優(yōu)先流雙域滲透模型對(duì)裂隙優(yōu)先流的模擬效果較好，典型剖面染色覆蓋率的模擬值與實(shí)測(cè)值極顯著相關(guān)（決定系數(shù)2=0.951，均方根誤差=0.095）。

3）接近土壤飽和導(dǎo)水率的灌水強(qiáng)度無(wú)法觸發(fā)優(yōu)先流，灌水強(qiáng)度提高，基質(zhì)流深度和灌水均勻度顯著降低，優(yōu)先流程度增加；初始含水率和裂隙分布密度增大，基質(zhì)流深度和灌水均勻度增大，優(yōu)先流程度減弱。

[1] NIMMO J R. The processes of preferential flow in the unsaturated zone[J]. Soil Science Society of America Journal, 2021, 85(1): 1-27.

[2] DADFAR H, ALLAIRE S E, DE JONG R, et al. Development of a method for estimating the likelihood of crack flow in Canadian agricultural soils at the landscape scale[J]. Canadian Journal of Soil Science, 2010, 90(1): 129-149.

[3] NOVáK V, ?IM?UNEK J, VAN GENUCHTEN M T. Infiltration of water into soil with cracks[J]. Journal of Irrigation and Drainage Engineering, 2000, 126(1): 41-47.

[4] LIU C W, CHENG S W, YU W S, et al. Water infiltration rate in cracked paddy soil[J]. Geoderma, 2003, 117(1/2): 169-181.

[5] LI J H, ZHANG L M. Geometric parameters and REV of a crack network in soil[J]. Computers and Geotechnics, 2010, 37(4): 466-475.

[6] GUO L, LIN H. Addressing two bottlenecks to advance the understanding of preferential flow in soils[J]. Advances in Agronomy, 2018, 147: 61-117.

[7] QUISENBERRY V L, PHILLIPS R E. Percolation of surface-applied water in the field[J]. Soil Science Society of America Journal, 1976, 40(4): 484-489.

[8] HARDIE M A, COTCHING W E, DOYLE R B, et al. Effect of antecedent soil moisture on preferential flow in a texture-contrast soil[J]. Journal of Hydrology, 2011, 398(3/4): 191-201.

[9] JARVIS N J. A review of non-equilibrium water flow and solute transport in soil macropores: Principles, controlling factors and consequences for water quality[J]. European Journal of Soil Science, 2020, 71(3): 279-302.

[10] WELLS R R, DICARLO D A, STEENHUIS T S, et al. Infiltration and surface geometry features of a swelling soil following successive simulated rainstorms[J]. Soil Science Society of America Journal, 2003, 67(5): 1 344-1 351.

[11] LUO L F, LIN H, LI S C. Quantification of 3-D soil macropore networks in different soil types and land uses using computed tomography[J]. Journal of Hydrology, 2010, 393(1/2): 53-64.

[12] 唐朝生, 施斌, 劉春, 等. 影響?zhàn)ば酝帘砻娓煽s裂縫結(jié)構(gòu)形態(tài)的因素及定量分析[J]. 水利學(xué)報(bào), 2007, 38(10): 1 186-1 193.

TANG Chaosheng, SHI Bin, LIU Chun, et al. Factors affecting the surface cracking in clay due to drying shrinkage[J]. Journal of Hydraulic Engineering, 2007, 38(10): 1 186-1 193.

[13] DEMAND D, BLUME T, WEILER M. Spatio-temporal relevance and controls of preferential flow at the landscape scale[J]. Hydrology and Earth System Sciences, 2019, 23: 4 869-4 889.

[14] 王策, 張展羽, 陳曉安, 等. 基于水量平衡原理的裂隙優(yōu)先流雙域滲透模型及其應(yīng)用[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2021, 52(10): 314-326, 348.

WANG Ce, ZHANG Zhanyu, CHEN Xiaoan, et al. Dual-permeability model for crack preferential flow based on principle of water volume balance and its application[J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(10): 314-326, 348.

[15] 王康, 張仁鐸, 繳錫云. 多孔介質(zhì)中非均勻流動(dòng)特性的染色示蹤試驗(yàn)研究[J]. 水科學(xué)進(jìn)展, 2007, 18(5): 662-667.

WANG Kang, ZHANG Renduo, JIAO Xiyun. Characterizing heterogeneity of water flow and solute transport in the porous media using dye tracer[J]. Advances in Water Science, 2007, 18(5): 662-667.

[16] 朱磊, 尤今, 陳玖泓. 裂隙網(wǎng)絡(luò)對(duì)坡面流及土壤水分入滲影響[J]. 灌溉排水學(xué)報(bào), 2017, 36(12): 95-100.

ZHU Lei, YOU Jin, CHEN Jiuhong. The influence of fracture network on overland flow and water infiltration in soil[J]. Journal of Irrigation and Drainage, 2017, 36(12): 95-100.

[17] COPPOLA A, COMEGNA A, DRAGONETTI G, et al. Simulated preferential water flow and solute transport in shrinking soils[J]. Vadose Zone Journal, 2015, 14(9): 1-22.

[18] 雷志棟, 楊詩(shī)秀, 謝森傳. 土壤水動(dòng)力學(xué)[M]. 北京: 清華大學(xué)出版社, 1988.

LEI Zhidong, YANG Shixiu, XIE Senchuan. Soil hydrodynamics[M]. Beijing: Tsinghua University Press, 1988.

[19] 張展羽, 陳于, 孔莉莉, 等. 土壤干縮裂縫幾何特征對(duì)入滲的影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015, 46(10): 192-197.

ZHANG Zhanyu, CHEN Yu, KONG Lili, et al. Impacts of desiccation crack geometric characteristics on infiltration in soil[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(10): 192-197.

[20] WANG C, FENG G X, ZHANG Z Y, et al. Geometrical and statistical analysis of dynamic crack morphology in shrink-swell soils with addition of maize roots or salinity (NaCl)[J]. Soil and Tillage Research, 2021, 212: 105 057.

[21] 張俊, 牛文全, 張琳琳, 等. 初始含水率對(duì)微潤(rùn)灌溉線源入滲特征的影響[J]. 排灌機(jī)械工程學(xué)報(bào), 2014, 32(1): 72-79.

ZHANG Jun, NIU Wenquan, ZHANG Linlin, et al. Effects of soil initial water content on line-source infiltration characteristic in moistube irrigation[J]. Journal of Drainage and Irrigation Machinery Engineering, 2014, 32(1): 72-79.

[22] SMET S, BECKERS E, PLOUGONVEN E, et al. Can the pore scale geometry explain soil sample scale hydrodynamic properties?[J]. Frontiers in Environmental Science, 2018, 6: 20.

[23] SAMMARTINO S, MICHEL E, CAPOWIEZ Y. A novel method to visualize and characterize preferential flow in undisturbed soil cores by using multislice helical CT[J]. Vadose Zone Journal, 2012, 11(1): 1-13.

[24] TANG C S, SHI B, CUI Y J, et al. Desiccation cracking behavior of polypropylene fiber-reinforced clayey soil[J]. Canadian Geotechnical Journal, 2012, 49(9): 1 088-1 101.

[25] EDWARDS W M, SHIPITALO M J, OWENS L B, et al. Rainfall intensity affects transport of water and chemicals through macropores in no-till soil[J]. Soil Science Society of America Journal, 1992, 56(1): 52-58.

[26] TOKUNAGA T K, WAN J M. Approximate boundaries between different flow regimes in fractured rocks[J]. Water Resources Research, 2001, 37(8): 2 103-2 111.

[27] GJETTERMANN B, NIELSEN K L, PETERSEN C T, et al. Preferential flow in sandy loam soils as affected by irrigation intensity[J]. Soil Technology, 1997, 11(2): 139-152.

[28] HOPKINS I, GALL H, LIN H. Natural and anthropogenic controls on the frequency of preferential flow occurrence in a wastewater spray irrigation field[J]. Agricultural Water Management, 2016, 178: 248-257.

[29] QI W, ZHANG Z Y, WANG C, et al. Crack closure and flow regimes in cracked clay loam subjected to different irrigation methods[J]. Geoderma, 2020, 358: 113 978.

[30] CHEN C C, ROSEBERG R J, SELKER J S. Using microsprinkler irrigation to reduce leaching in a shrink/swell clay soil[J]. Agricultural Water Management, 2002, 54(2): 159-171.

[31] T?UMER K, STOFFREGEN H, WESSOLEK G. Seasonal dynamics of preferential flow in a water repellent soil[J]. Vadose Zone Journal, 2006, 5(1): 405-411.

[32] BISDOM E B A, DEKKER L W, SCHOUTE J F T. Water repellency of sieve fractions from sandy soils and relationships with organic material and soil structure[J]. Geoderma, 1993, 56(1/2/3/4): 105-118.

[33] BEVEN K, GERMANN P. Macropores and water flow in soils[J]. Water Resources Research, 1982, 18(5): 1 311-1 325.

[34] WEILER M, NAEF F. Simulating surface and subsurface initiation of macropore flow[J]. Journal of Hydrology, 2003, 273(1/2/3/4): 139-154.

Effects of Cracking Patterns and Antecedent Soil Water Content on Preferential Flow of Irrigated Water

DUAN Fengyi1, ZHU Chengli1, WANG Ce1*, CHEN Yu2, MA Liang1

(1. College of Agricultural Science and Engineering, Hohai University, Nanjing 211100, China;2.Jiangsu Rural Water Conservancy Science and Technology Development Center, Nanjing 210029, China)

【Objective】Macropores formed by physical and biological processes provide fast pathways for rainfall and irrigation water to bypass soil matrix. It has important implications not only for hydrological processes but also for biogeochemical processes. In this paper, we investigated the effect of cracking patterns and antecedent soil moisture content on preferential flow of irrigated water.【Method】Our study was based on numerical simulation and experiment. We first validated a dual-permeability model against a dye tracing experiment, and then applied it to elucidate how crack patterns, antecedent soil moisture content and irrigation intensity combined to modulate the formation and development of preferential flow of irrigated water in soil.【Result】Comparison with experimental results showed that the model is accurate for predicting preferential flow, with2=0.951 and<0.01. Shallow small cracks with aperture in the range of 0.2～1.0 cm and in the depth of 0.5～10 cm had little impact on preferential flow, and large cracks with aperture more than 2 cm and expanding to the depth beyond 50 cm were the mechanisms underlying the formation and development of preferential flow. An increase in crack density enhanced both infiltration depth in soil matrix and soil water uniformity. Irrigation intensity less than saturated soil hydraulic conductivity did not trigger preferential flow. Increasing irrigation intensity increased infiltration depth in soil matrix and soil water uniformity in the soil. An increase in antecedent soil water content and cracking density increased water flow in soil matrix and soil water uniformity because of the reduced preferential flow.【Conclusion】In the fields with soil shrinkages and rich organic matters, reducing irrigation intensity and irrigating before soil water content drops below a certain level can reduce preferential flow and improve irrigation uniformity consequently.

crack preferential flow; dual-permeability models; fracture distribution density; antecedent water content; irrigation intensity

1672 - 3317（2023）02 - 0032 - 10

S152.7

A

10.13522/j.cnki.ggps.2022263

段鳳儀, 朱成立, 王策, 等. 灌溉優(yōu)先流對(duì)裂隙分布密度、灌水強(qiáng)度及土壤濕度的響應(yīng)機(jī)理[J]. 灌溉排水學(xué)報(bào), 2023, 42(2): 32-41.

DUAN Fengyi, ZHU Chengli, WANG Ce, et al.Effects of Cracking Patterns and Antecedent Soil Water Content on Preferential Flow of Irrigated Water[J]. Journal of Irrigation and Drainage, 2023, 42(2): 32-41.

2022-05-11

國(guó)家自然科學(xué)基金青年項(xiàng)目（52109053）；江蘇省自然科學(xué)基金青年基金項(xiàng)目（BK20200523）；江西省水利科學(xué)院開放研究基金項(xiàng)目（2021SKTR03）；國(guó)家科技獎(jiǎng)后備培育項(xiàng)目（20212AEI91011）；江西省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（20203BBGL73226）

段鳳儀（1997-），女。碩士研究生，主要從事農(nóng)田裂隙及裂隙優(yōu)先流理論研究。E-mail: dfy@hhu.edu.cn

王策（1989-），男。教授，主要從事農(nóng)田裂隙及其優(yōu)先流理論、高效灌排理論研究。E-mail: wangce@hhu.edu.cn

責(zé)任編輯：白芳芳