邵培寅，李亞龍，熊玉江，袁念念，蘇沛蘭

（1.太原理工大學(xué) 水利科學(xué)與工程學(xué)院，太原 030024；2.長江水利委員會長江科學(xué)院 農(nóng)業(yè)水利研究所，武漢 430010）

?農(nóng)田排水?

丘陵灌區(qū)排水循環(huán)灌溉模式節(jié)水減排效果研究

邵培寅1,2，李亞龍2，熊玉江2，袁念念2，蘇沛蘭1*

（1.太原理工大學(xué) 水利科學(xué)與工程學(xué)院，太原 030024；2.長江水利委員會長江科學(xué)院 農(nóng)業(yè)水利研究所，武漢 430010）

【目的】揭示丘陵灌區(qū)較為普遍的排水循環(huán)灌溉過程及其節(jié)水減排效果?！痉椒ā酷槍鄥^(qū)典型灌排單元，采用原位試驗監(jiān)測農(nóng)溝和斗溝中水量和氮磷質(zhì)量濃度變化，分析回歸率（回歸水量/抽灌水量）、負荷削減量、負荷削減率時空變化規(guī)律，并通過相關(guān)分析、逐步回歸方法求解負荷削減率的影響因素。【結(jié)果】監(jiān)測到抽灌-回歸過程24次，總回歸率為89.93%。生育期總磷、總氮、硝態(tài)氮和氨態(tài)氮負荷削減量分別為0.28、3.27、2.35、2.35 kg/hm2?？偭?、總氮、氨態(tài)氮負荷的削減發(fā)生在田間和農(nóng)溝，硝態(tài)氮負荷的削減發(fā)生在斗溝?？嘏糯胧┦沟媒邓笈潘摵芍械南鯌B(tài)氮和氨態(tài)氮達到了穩(wěn)定。4種負荷削減率都與其質(zhì)量濃度之比（回歸質(zhì)量濃度/抽灌質(zhì)量濃度）顯著相關(guān)，總磷、氨態(tài)氮削減率與回歸率顯著相關(guān)，總氮、硝態(tài)氮削減率與抽灌當(dāng)天及第2天降水量顯著相關(guān)，硝態(tài)氮削減率還與溫度顯著相關(guān)?！窘Y(jié)論】抽灌-回歸-控排灌溉模式在節(jié)水、方便農(nóng)戶灌溉的同時，提高了水肥利用率，降低了排水給下游水體帶來的風(fēng)險，在灌排單元尺度上節(jié)水和控污減排效果良好。

丘陵灌區(qū)；循環(huán)灌溉；氮磷負荷；水稻；控污減排

0 引 言

【研究意義】灌區(qū)灌溉和排水受人為干預(yù)較多[1]。丘陵灌區(qū)因塘堰、溝渠和水庫眾多，降水徑流、灌溉退水通過塘堰、溝渠等水體受到人為攔蓄進行重復(fù)利用現(xiàn)象十分普遍[2]。灌區(qū)排水再利用作為一種調(diào)蓄雨水、提高灌溉保證率和減少農(nóng)業(yè)面源污染的有效措施[3]，已被較多地區(qū)重視并加以采用[4-5]?！狙芯窟M展】根據(jù)排水再利用的對象，排水再利用可分為2種，一種是排水自流灌溉下級渠道（灌溉回歸水重復(fù)利用[6]），另一種是排水依靠水泵回灌本級渠道（循環(huán)灌溉[7]），相較而言，循環(huán)灌溉由于水肥重復(fù)利用，下游受到污染的風(fēng)險更小。Hama等[8-9]研究發(fā)現(xiàn)，循環(huán)灌溉能提供作物需氮量的8%～16%，Takeda等[10]匯總了8個灌溉期的監(jiān)測資料，發(fā)現(xiàn)循環(huán)灌溉的負荷削減量與田間水停留時間有關(guān)。Feng等[7]發(fā)現(xiàn)循環(huán)灌溉系統(tǒng)在暴雨期間稻田營養(yǎng)流失大幅減少。焦平金等[11]、于穎多等[12]評價了循環(huán)灌溉排水量估算模型，發(fā)現(xiàn)田面水和滲透水的磷素不受灌溉水中磷素的變化影響?！厩腥朦c】然而以上研究大多是針對流域尺度或田間尺度的節(jié)水效益及削減效果分析，涉及灌排單元尺度和循環(huán)灌溉各環(huán)節(jié)對于節(jié)水及削減效果的影響因素研究少見報道。另外，控制、存蓄排水是循環(huán)灌溉中不可缺少的一環(huán)。而關(guān)于溝渠控制排水的研究[13-15]大多是為調(diào)控灌排單元地下水位以增加作物對地下水的利用[16-17]、發(fā)揮溝渠濕地效應(yīng)以截污減排[18-19]，而非為了存蓄水分以供灌溉。涉及農(nóng)戶自發(fā)循環(huán)灌溉而促成的控制排水研究少見報道。【擬解決的關(guān)鍵問題】本研究針對丘陵灌區(qū)典型灌排單元進行原位觀測，將循環(huán)灌溉中的水流運動細化為抽灌、回歸和控排3個環(huán)節(jié)，具體分析抽灌-回歸-控排灌溉模式下各個環(huán)節(jié)水流和負荷的動態(tài)變化及其影響因素，綜合評價灌排單元的節(jié)水效果及控污減排機理，為丘陵灌區(qū)的排水再利用和非點源污染防控相結(jié)合進一步研究和推廣提供借鑒和參考。

1 材料與方法

1.1 研究地點概況

研究區(qū)位于漳河灌區(qū)三干渠南1.0 km處，灌區(qū)主要灌溉水源為10 km外的漳河水庫。灌區(qū)配有完備的干、支、斗、農(nóng)灌排渠系及各種小型水庫、塘堰，灌區(qū)地形主要為丘陵與平原結(jié)合區(qū)，其自然條件在南方丘陵地區(qū)具有一定的代表性。灌區(qū)中小型水庫300多座，塘堰約8.16萬座，星羅棋布，形成了“長藤結(jié)瓜”式灌溉系統(tǒng)。

漳河灌區(qū)氣候?qū)儆趤啛釒Ъ撅L(fēng)氣候，年平均氣溫16 ℃，多年平均降水量968.9 mm，降水集中在4—9月，此時期降水量約占全年降水量的70%，5—8月為密集降水期。主要農(nóng)作物為1年2季稻油（麥）輪作。試驗區(qū)的水稻田在5月下旬播種，9月上旬收割，水稻生長期近112 d。稻季農(nóng)民習(xí)慣采用底肥加追肥的方式，底肥在水稻移栽前整田施用，追肥在6月下旬施用。

1.2 試驗設(shè)計

整個循環(huán)灌溉過程可分為抽灌、回歸和控排3個過程。在水稻生長期，研究區(qū)內(nèi)農(nóng)戶通常用水泵將田塊周邊溝渠中的蓄水抽灌至田間（抽灌），田間水通過田埂側(cè)滲或降水后自流進入本級或下級斗溝（回歸），溝渠下游被堵住，用來存蓄雨后田間排水和田埂側(cè)向滲漏水（控排）。蓄水供平日抽灌使用，需要灌溉時再抽灌，如此循環(huán)往復(fù)。

選取循環(huán)灌溉典型灌排單元作為研究對象（圖1）。2相鄰田塊中間田埂寬15 cm，田塊1比田塊2低10 cm，田埂處開有缺口，田塊2田間水層超過5 cm時可自流排入田1。田塊1上方田埂寬15 cm，開有缺口，田塊1田間水層超過5 cm時可自流進入相鄰農(nóng)溝，農(nóng)溝與斗溝相通，長度分別為20 m和180 m，農(nóng)溝斷面為底寬0.5 m，深0.5 m的矩形。斗溝斷面為底寬0.7 m、坡比1∶1的梯形，斗溝深1 m。田塊總面積2 440 m2。保持下游涵洞沙袋高度不變，僅當(dāng)降水發(fā)生斗溝水位超過沙袋高度（約0.5 m）時才向下游排水。試驗期間為排除降水外其他水量輸入項，將農(nóng)溝上游和斗溝上游用麻袋堵住，并定期監(jiān)測斗溝上游段水位，水位較高時抽灌至下游堰塘，以防上游來水。周邊田塊相對灌排單元一側(cè)田埂鋪設(shè)防滲土工膜，防止周邊田塊田間水向灌排單元側(cè)滲。因此，整個灌排單元可看做封閉區(qū)域。本研究為原位試驗，對農(nóng)戶的田間管理和抽灌行為只做記錄，不做任何干涉。

1.3 試驗觀測與測試方法

1.3.1 水量觀測與分析

本研究對抽灌-回歸-控排灌溉模式中的抽灌水量、回歸水量、排水量和降水量進行監(jiān)測。

1）抽灌水量為斗溝抽灌至田間的水量。農(nóng)戶每次抽灌需要4～5 h。根據(jù)斗溝每日水位波動及斗溝梯形斷面推求各次抽灌水量，生育期斗溝水位變化及典型抽灌過程曲線見圖2，人工抽灌前后斗溝庫容差即為抽灌水量。

圖1 典型灌排單元及控排抽灌現(xiàn)象Fig.1 Schematic diagram of typical irrigation and drainage unit and phenomenon of pumping irrigation and drainage controlled

圖2 生育期斗溝水位變化及典型抽灌過程曲線Fig.2 Variation of water level in bucket ditch during growth period and typical pumping irrigation process curve

2）回歸水量為抽灌田間水層超過田埂排水口高度自流排出田間的水量與田埂側(cè)滲水量之和，回歸水先匯入農(nóng)溝再排入斗溝，回歸水量即為農(nóng)溝排水量。農(nóng)溝出口設(shè)有巴歇爾槽，使用HOBO水位計監(jiān)測每小時水位，對應(yīng)巴歇爾槽公式[20]求得回歸水量。

3）排水量為斗溝水位高出涵洞沙袋，通過涵洞向下游排水量。在斗溝下游涵洞出口處設(shè)置HOBO水位計監(jiān)測每小時水位，對應(yīng)圓形斷面自由出流水位流量關(guān)系曲線[18]，求得排水量。

4）降水量。通過研究區(qū)內(nèi)設(shè)置的HOBO RG3翻斗式自記雨量記錄儀記錄。

1.3.2 水質(zhì)觀測與分析

在農(nóng)溝和斗溝上下游分別設(shè)置2個取樣點，每次取樣對2種溝別各取上下游混合樣。取樣頻率為3 d 1次，遇降水、施肥則在事件發(fā)生第1、第3、第5、第7天加取。另外在生育期取3次降水水樣。2種溝別及雨水樣總計取樣35次，143個。每次取樣200 mL，原狀水樣送實驗室分析。水樣分析指標包括總磷（TP）、總氮（TN）、氨態(tài)氮（NH4+-N）和硝態(tài)氮（NO3+-N）。分析方法按照《水和廢水監(jiān)測分析方法》[21]中相關(guān)規(guī)定進行。

2 結(jié)果與分析

本研究循環(huán)灌溉模式發(fā)生在返青期、分蘗期、拔節(jié)孕穗期和抽穗開花期，其他生育期抽灌發(fā)生較少。經(jīng)過篩選，計算分析6月8日—8月8日的抽灌-回歸-控排事件。

2.1 排水量、抽灌量和回歸水量變化

試驗區(qū)降水集中發(fā)生于6月27日—7月19日，合計降水量為120.8 mm（圖4）。連續(xù)降水期間，田間水層較高，不需要灌溉，抽灌事件發(fā)生較少。斗溝蓄水超過下游涵洞沙袋高度外排，因此有連續(xù)排水事件發(fā)生，排水量也較大。而在其他時間，降水較少且農(nóng)戶多采用濕潤灌溉（間歇淹水），田間水層保持的較低（5 cm）。因此抽灌事件發(fā)生較多，而由于控排措施的存在，灌排單元排水量較少?；貧w過程一般在抽灌過程結(jié)束后1～2 d內(nèi)完成，一般視田間濕潤情況而定。試驗期間，共計監(jiān)測排水3 430.17 m3，抽灌回歸過程24次，平均每次抽灌水量77 m3。這與Hama等[9]在日本比瓦湖面流域85.3%的結(jié)果相比稍大，原因是本研究的灌排單元尺度與流域尺度相比，從稻田到蓄水點回歸距離較短，蒸發(fā)損失更少，且本研究生育期降水量更大。

2.2 灌排單元內(nèi)氮磷濃度變化特征研究

為分析循環(huán)灌溉回歸過程從田間到農(nóng)溝再到斗溝氮磷污染物質(zhì)量濃度的變化特征，選用標準差S和標準差率