王 珍，郝鋒珍,2，李久生，栗巖峰

基于EPANET的再生水滴灌系統(tǒng)余氯分布模型構建

王 珍1，郝鋒珍1,2，李久生1※，栗巖峰1

（1. 中國水利水電科學研究院流域水循環(huán)模擬與調控國家重點實驗室，北京 100048；2. 山西農業(yè)大學城鄉(xiāng)建設學院，太谷 030801）

加氯處理是保證滴灌系統(tǒng)安全的重要措施之一，但是氯與灌溉水及滴灌管壁發(fā)生反應而引起的衰減會對余氯分布均勻性造成不利影響。該研究采用EPANET軟件構建了再生水滴灌系統(tǒng)水力性能及余氯運移、衰減和分布模型，基于試驗數(shù)據(jù)對模型水力和水質模擬參數(shù)進行率定和驗證，并應用該模型分析了加氯模式（濃度-時長：1.3 mg/L-185 min、3 mg/L-80 min、6 mg/L-40 min和8 mg/L-30 min）和毛管長度（10～150 m）對系統(tǒng)余氯分布特征的影響。結果表明，EPANET軟件可以較好地模擬滴灌系統(tǒng)水力性能及余氯分布特征，沿毛管方向余氯濃度模擬值與實測值一致性指數(shù)達到0.75以上。當毛管長度由10 m增加至150 m時，灌水單元平均余氯濃度與毛管末端余氯濃度隨毛管長度增加呈先升高而后降低趨勢，高濃度加氯處理平均余氯濃度和毛管末端余氯濃度隨毛管長度增加衰減速率相對較小。較短的加氯歷時（< 40 min）和過低的加氯濃度（如1.3 mg/L）均易造成余氯質量偏差率升高和均勻系數(shù)降低；建議使用加氯濃度3 mg/L且加氯時長80 min的加氯模式，該模式毛管極限鋪設長度達到66 m，余氯質量均勻系數(shù)超過90%。

氯；模擬；水力性能；水質；滴灌；EPANET

0 引 言

灌水器的堵塞防止和處理一直是再生水滴灌系統(tǒng)安全運行的核心問題[1]。通過向再生水滴灌系統(tǒng)加氯，可以有效抑制微生物生長，減緩灌水器堵塞形成[2]。Ravina等[3]研究指出余氯濃度控制在3～5 mg/L時可有效控制灌水器堵塞；Cararo等[4]則指出控制余氯濃度為0.5 mg/L對于控制灌水器堵塞更為有效；Li等[5-6]發(fā)現(xiàn)低濃度、長歷時的加氯模式更有利于減輕系統(tǒng)堵塞的發(fā)生。但是，因氯具有強氧化性，其在緩解灌水器堵塞的同時可能會對土壤環(huán)境及作物根系產生不利影響[7]。再生水滴灌適宜加氯模式的制定需綜合考慮加氯對系統(tǒng)性能和作物生長的共同影響[8-9]。Li等[8]發(fā)現(xiàn)控制余氯濃度為2～8 mg/L、加氯間隔為4～8周既能控制系統(tǒng)堵塞，又不會對番茄生長、產量和品質產生明顯影響；Hao等[9]發(fā)現(xiàn)加氯一定程度上抑制了土壤酶活性，較低的余氯濃度不會對玉米生長造成影響。上述研究為再生水滴灌系統(tǒng)加氯模式優(yōu)化提供了重要理論基礎。但是，加氯過程中余氯濃度會隨氯與灌溉水及管壁發(fā)生反應而衰減[10]，在灌水器流量不均勻性的疊加作用下，滴灌系統(tǒng)中可能出現(xiàn)部分灌水器加氯過量或加氯不足的情況，導致部分區(qū)域土壤環(huán)境惡化而部分區(qū)域堵塞控制效果不佳。定量研究再生水滴灌系統(tǒng)中余氯分布特征及其影響因素對進一步優(yōu)化加氯模式具有重要意義。

通過水力解析（如能量輪廓線法）或數(shù)值計算方法（如有限元法）對管網水力特性進行計算，并對灌水均勻性進行評價，是滴灌系統(tǒng)管網優(yōu)化的基礎[11-13]。但是，以往滴灌系統(tǒng)水力特性計算方法中較少考慮溶質遷移轉化過程及其在灌溉系統(tǒng)的運移分布特征，一定程度上限制了溶質（如肥料和氯）注入模式對系統(tǒng)性能的影響研究。近年來，基于EPANET構建的城鎮(zhèn)供水管網余氯分布模擬模型已在模擬不同管網布置形式、供水方式和加氯模式條件下余氯分布特征中得到廣泛應用[10,14-15]。與城鎮(zhèn)供水管網相比，滴灌系統(tǒng)出流節(jié)點（灌水器）更多，且節(jié)點流量受水力偏差影響明顯，相對較低的滴灌系統(tǒng)管網平均流速（約0.3 m/s）可能加大氯或肥料自加注點運移至不同位置灌水器的時間差異[16]，能否基于EPANET軟件實現(xiàn)滴灌系統(tǒng)中余氯濃度和肥料分布的模擬是一個值得關注的問題。本研究基于EPANET在供水管網水力和水質模擬中的優(yōu)勢，建立再生水滴灌系統(tǒng)余氯分布模型，并基于該模型評價滴灌系統(tǒng)參數(shù)及加氯模式對余氯分布特征的影響，以期為再生水灌溉條件下的系統(tǒng)設計及加氯模式優(yōu)選提供科學依據(jù)。

1 模型構建及試驗方法

EPANET是美國環(huán)境保護署（Environmental Protection Agency，EPA）開發(fā)的有壓管網水力和水質特性延時模擬軟件，可實現(xiàn)不同類型的供水系統(tǒng)水力及水質特性分析[17]，已在余氯衰減運移中得到了大量應用。

1.1 再生水滴灌加氯系統(tǒng)模型構建

基于EPANET軟件，本研究構建了滴灌系統(tǒng)典型灌水單元，包括水源1處、干管1根、支管1根、毛管若干條及灌水器若干個（圖1）。滴灌系統(tǒng)毛管數(shù)量、單條毛管上灌水器水量、毛管間距、灌水器間距和灌水器流量等技術參數(shù)可根據(jù)系統(tǒng)實際情況取值。本研究假定加氯點位于水源處，加入的氯通過干管—支管—毛管—灌水器后滴入土壤。

圖1 滴灌加氯系統(tǒng)組成示意圖

1.1.1 水力模擬

管網水力計算分析包括管網節(jié)點流量連續(xù)方程和管段能量方程的聯(lián)立求解。節(jié)點流量連續(xù)方程如式（1）所示：

式中Q為與節(jié)點相連的管段的流量，L/h；Q為節(jié)點流量，L/h；為節(jié)點編號；node為節(jié)點數(shù)量。管段能量方程如式（2）所示：

式中H和H分別為節(jié)點和處的水頭，m；h為管段水頭損失，m；為阻力系數(shù)；為流量指數(shù)。本研究中采用Hazen-Williams公式求解管段的水頭損失。EPANET采用Todini-Pilati梯度算法求解給定時間點管網水力狀態(tài)的流量連續(xù)性方程和水頭損失方程組[18]。滴灌系統(tǒng)水力性能模擬時，灌水器節(jié)點類型設置成“Emitters”，其節(jié)點出水量與節(jié)點壓力關系如式（3）所示：

=KP（3）

式中為灌水器實際流量，L/h；為流量系數(shù)；P為壓力水頭，m；為灌水器流態(tài)指數(shù)。

滴灌系統(tǒng)水力模擬時，所需的水力參數(shù)主要包括水源工作壓力、干管、支管和毛管阻力系數(shù)（R、R和R）、灌水器流量系數(shù)和流態(tài)指數(shù)及節(jié)點初始流量（Q，L/h）。其中，工作壓力可根據(jù)實際情況進行設置；灌水器流態(tài)指數(shù)和流量系數(shù)可由實測的壓力-流量關系曲線擬合得到；阻力系數(shù)根據(jù)管道材質設置初始值，并經實測的滴灌系統(tǒng)灌水器流量進行率定。

1.1.2 水質模擬

針對構建的滴灌系統(tǒng)，加氯條件下管網水質模擬的主要任務是在管網水力計算的基礎上，基于管網中溶質對流傳輸規(guī)律，預測余氯（水經過加氯消毒接觸一定時間后，水中所余留的以游離氯、化合氯或兩者并存形式存在的氯）[19]在管網內的時空變化。有效氯的傳輸衰減過程包括管道內物質的對流傳輸過程、物質動態(tài)反應過程及節(jié)點混合引起的物質濃度變化。管道內物質的對流傳輸方程見式（4）[17]：

式中C為管道中的余氯濃度，mg/L；u為管道中的流速，m/h；(C)為管道中反應物的反應速率，對于非反應物，(C)=0，mg/(L·h)；為總余氯衰減系數(shù)，1/h；為反應級數(shù)；為管道編號；pipe為管道數(shù)量。。

本研究忽略加入氯在水源中的混合過程，假設水源處的加氯濃度（有效氯，以Cl計，下同）[20]恒定，進而通過給定水源處加氯濃度的方式對余氯分布進行模擬。采用一級氯衰減模型來模擬滴灌系統(tǒng)內的余氯衰減，則=1。本研究中設定滴灌系統(tǒng)干管、支管和毛管采用PVC或PE管材，可忽略管壁腐蝕引起的氯衰減。因此，余氯衰減過程主要包括主體水余氯衰減和管壁衰減，二者之和構成管網余氯衰減[21]：

式中k為主體水衰減系數(shù)，1/h；k為管壁衰減系數(shù)，m/h；r為水力半徑，m；k為質量傳輸系數(shù)，模型中根據(jù)舍伍德數(shù)、管道長度、水力半徑、雷諾數(shù)等進行計算。

1.2 試驗材料與方法

1.2.1 滴灌系統(tǒng)水力性能試驗

試驗于2017年9月在國家節(jié)水灌溉北京工程技術研究中心大興試驗基地進行。在田間構建滴灌系統(tǒng)1套，滴灌系統(tǒng)包括：水源1處（4 m×4 m×2.5 m水池），用于提供灌溉用的二級處理再生水；鄰近水源處設置比例泵1臺（Mis Rite Model 2504，Tefen，以色列），用于向系統(tǒng)中注入有效氯溶液，使系統(tǒng)加氯濃度與設計值一致；精密壓力表1個（0.4%，0.40 MPa），用于測定干管首端工作壓力；干管（50 mm PVC管）1條，長度25 m，內徑46 mm；干管末端設置支管（40 mm PVC管）1條，長度36 m，內徑36 mm；支管末端設置內鑲貼片式滴灌帶1條（TYPHOON，16 mm，耐特菲姆），內徑15.6 mm，灌水器間距0.3 m，0.1 MPa條件下灌水器標稱流量1.6 L/h。試驗中，設置滴灌帶長度為100和50 m，分別記為滴灌系統(tǒng)A和B。

試驗開始前按照ASAE標準（S553）測定灌水器流量-壓力關系[22]，計算得到灌水器流量系數(shù)為0.495，流態(tài)指數(shù)為0.476。試驗中，利用閥門控制干管首部壓力為0.1 MPa。沿毛管方向，每隔2個灌水器設置1處灌水器流量監(jiān)測點，利用水桶承接20 min內灌水器出流，利用量筒測定承接水量，以此計算沿毛管方向灌水器流量。水力性能測試設置3個重復，每個重復測試時，干管和支管不變，對滴灌帶進行更換，將3次測定值進行平均得到不同觀測點處灌水器流量。

1.2.2 滴灌系統(tǒng)水質測定試驗

1）主體水衰減試驗

再生水選用二級處理再生水，取自北京市大興區(qū)黃村污水處理廠，水質指標如表1所示。主體水衰減由室內燒杯試驗確定。氯原材料選用10%有效氯（有效系數(shù)×含氯量）的次氯酸鈉，余氯濃度采用便攜式余氯計（ExStik CL200，Extech Instruments Corporation，美國）測定。為研究加氯濃度對再生水主體水衰減規(guī)律的影響，本研究控制水溫為20℃，參考常用加氯濃度[3-5]，設置加氯濃度為1.45、2.30、3.98、8.14 mg/L，研究加氯2 h內的氯衰減變化規(guī)律，分別于1、3、5、7、9、11、20、30、40、60、90、120 min測定余氯值并記錄。

表1 試驗用再生水水質指標

2）滴灌系統(tǒng)余氯濃度衰減試驗

在系統(tǒng)水力性能試驗的基礎上，通過比例泵向干管中注入恒定濃度的次氯酸鈉溶液，使加氯點有效氯濃度達到設計值。參照主體水衰減試驗，設置4個加氯濃度，分別為1.3、3.0、6.0和8.0 mg/L。為保證試驗環(huán)境溫度與主體水衰減試驗類似，所有試驗均在15:00－18:00進行，試驗過程中環(huán)境溫度范圍為20～25 ℃。對于滴灌系統(tǒng)，分別在干管末端和支管末端設置1處水質監(jiān)測點，沿毛管方向等間距布置10個水質監(jiān)測點，用以對余氯濃度分布進行測定。對于每組測試，加氯總時長為2 h。自加氯開始至加氯1 h期間，不進行取樣；加氯1 h后（假定此時系統(tǒng)不同位置余氯濃度不隨時間變化），按遠離加氯點方向依次在干管末端、支管末端和毛管上進行灌水器出流采樣，立即進行余氯測試（每次采樣及余氯測試約耗時4 min），得到各位置余氯濃度。與水力性能試驗對應，對于每個加氯濃度也設置3個重復，每個重復中更換滴灌帶后進行余氯濃度測試，將3次測定平均值作為監(jiān)測點處余氯濃度。

1.3 模型評價

本研究中，模型模擬效果評價分水力模擬評價和水質模擬評價兩部分組成，分別以灌水器流量及出流液中余氯濃度為對象進行評估，模擬值和觀測值的吻合程度采用均方根誤差（Root Mean Square Error，RMSE）和一致性指數(shù)（）進行評價[23]。RMSE最小值為0，越接近0模擬效果越好；值范圍為0~1，越接近1模擬效果越好。

2 模型率定與驗證結果與分析

2.1 水力參數(shù)

模型模擬時，利用系統(tǒng)A灌水器實測流量對參數(shù)進行率定，利用系統(tǒng)B灌水器實測流量對率定參數(shù)進行驗證。經率定和驗證，R、R和R分別為100、105和135。表2給出了灌水器流量模擬評價指標。灌水器流量模擬值和實測值平均流量基本一致，系統(tǒng)A、B均方根誤差（RMSE）分別為0.026和0.014 L/h，一致性指數(shù)（）分別為0.88和0.57，模擬值與實測值具有較好的一致性。綜上，利用EPANET軟件對滴灌系統(tǒng)水力性能能進行模擬具有較高的精度。

表2 灌水器流量模擬評價

2.2 水質參數(shù)

2.2.1 主體水衰減系數(shù)計算

不同加氯濃度條件下，再生水主體水衰減規(guī)律見圖2。不同加氯濃度條件下，余氯濃度衰減曲線的變化規(guī)律基本一致。

圖2 不同加氯濃度的余氯衰減規(guī)律

反應初期，余氯衰減迅速，大約30 min后衰減趨于平緩，這是因為該階段易于被氧化的有機物和無機物在反應初期占主導，氯的衰減速率快。當易于被氧化的物質消耗之后，不易與氯發(fā)生反應物質占主導，導致氯衰減速率相對較低。按照一級動力學式（7）對余氯衰減過程進行擬合[24]，得到加氯濃度為1.45、2.30、3.98和8.14 mg/L時的主體水衰減系數(shù)分別為0.415、0.301、0.222、0.132 1/h，2分別為0.791、0.750、0.955、0.810。

式中C為時刻的余氯濃度，mg/L；0為加氯濃度，mg/L。

利用冪函數(shù)擬合主體水衰減系數(shù)k與加氯濃度的關系（圖3），進而可基于擬合方程計算出系統(tǒng)試驗中加氯濃度為1.3、3、6和8 mg/L時，主體水衰減系數(shù)分別為0.446、0.257、0.163和0.136 1/h。

圖3 主體水衰減系數(shù)與加氯濃度關系

2.2.2 管壁衰減系數(shù)率定與驗證

參照EPANET技術手冊設定干管、支管及滴灌帶管壁的初始余氯衰減系數(shù)，分別基于系統(tǒng)A和B余氯測定結果對干管、支管和毛管管壁衰減系數(shù)（分別為wm、ws和wl）進行率定和驗證。圖4給出了滴灌系統(tǒng)A不同加氯濃度下的余氯模擬值與實測值比較結果。由圖可知，不同加氯濃度條件下，余氯濃度模擬值與實測值吻合程度較好。觀測點余氯濃度實測值及模擬值均隨距加氯點距離增大呈下降趨勢，這說明隨水流推進及運移時間增加，灌溉水中氯與灌溉水及管壁發(fā)生了明顯反應。表3給出了不同加氯濃度條件下余氯分布實測值和模擬值RMSE和。由表可知，參數(shù)率定后，不同加氯濃度條件模擬RMSE均相對較小（0.12～0.26 mg/L），余氯濃度實測值與模擬值差異不大；同時，一致性指數(shù)（）均大于0.75，模型模擬精度較高，這說明使用EPANET軟件進行不同規(guī)模滴灌系統(tǒng)余氯濃度的模擬是可行的。由表3還可看出，管壁衰減系數(shù)基本隨加氯濃度的升高而降低，這與Hallam等研究結果一致[21]，主要是因為加氯濃度低時，加入氯與灌溉水及管壁的反應比例較高造成的。

圖4 滴灌系統(tǒng)A不同加氯濃度下的不同位置余氯模擬值與實測值比較（率定）

表3 不同加氯濃度條件下管壁衰減系數(shù)及余氯模擬評價指標

3 模型應用與結果分析

3.1 模擬方案

模型應用考慮加氯模式（濃度-時長）和毛管長度對滴灌系統(tǒng)余氯濃度及余氯質量分布特征的影響。因水質會對余氯衰減過程產生影響，模擬評估時假定所用再生水水質與第2節(jié)中試驗條件一致。參照圖1滴灌系統(tǒng)，設定模擬模型中干管、支管長度分別為25和10 m，設定毛管數(shù)量為10條，以間距為1 m均勻分布于支管上，灌水器間距及模型水力及水質參數(shù)設置與第2節(jié)中率定結果一致。模擬時在保證加氯質量相同的條件下，設置1.3 mg/L-185 min（C1T4）、3 mg/L-80 min（C2T3）、6 mg/L-40 min（C3T2）和8 mg/L-30 min（C4T1）4種加氯模式。所有模式灌水時長均設定為4 h。為保證有效氯在加氯結束后仍部分存留于系統(tǒng)中，所有模式均先進行灌水，而后開始加氯，即對C1T4、C2T3、C3T2和C4T1模式分別在灌水35、160、200和210 min后開始加氯。毛管長度范圍設置為10~150 m，以10 m為間隔進行模擬，計算沿毛管方向水齡（自水源/加氯點運移至管網特定位置處所耗時間）、余氯濃度分布及平均濃度，評價加氯模式及毛管長度對余氯濃度分布的影響；計算灌水單元灌水器流量和灌水器出流余氯質量（∑（灌水器流量×余氯濃度×加氯時間））偏差率S和克里斯琴森均勻系數(shù)（Coefficient of Uniformity，CU）評價加氯模式和毛管長度的影響。

式中max和min分別是灌水器流量（L/h）或出流余氯質量（mg）最大值和最小值。

3.2 結果與分析

3.2.1加氯濃度及毛管長度對余氯濃度分布的影響

圖5給出了不同加氯濃度條件下典型毛管長度（10、80和150 m）水齡及余氯濃度沿毛管方向變化。由圖5可知，系統(tǒng)內不同位置節(jié)點水齡差異明顯，10、80和150 m滴灌系統(tǒng)沿毛管方向水齡變化范圍分別為5.8～12.8、1.2～16.9、1.1～27.3 min，水齡沿遠離進水口方向持續(xù)增加，在灌水單元后半段，因管段流速降低，灌水器節(jié)點水齡增加趨勢更為明顯。沿毛管方向水齡的變化決定了氯與灌溉水及管壁的反應時間差異[15]，是灌水器出流中余氯濃度的關鍵影響因素。隨著氯與水體及管壁的反應，不同加氯濃度條件下余氯濃度均沿遠離進水口方向呈現(xiàn)下降趨勢（圖5）。

圖5 不同加氯濃度典型毛管長度水齡及余氯濃度沿毛管方向變化

同一毛管長度，余氯濃度衰減速度隨加氯濃度的增大而減小。例如，毛管長度為80 m時，加氯濃度為1.3和6.0 mg/L時沿毛管方向余氯濃度最大值與最小值之比分別為2.0和1.4，這與試驗中低加氯濃度時較大的主體水衰減系數(shù)和管壁衰減系數(shù)有關。同一加氯濃度，毛管長度不同時，灌水單元中余氯濃度差異基本隨毛管長度的增加而增大。例如，對于3.0 mg/L加氯濃度，10、80和150 m毛管長度時余氯濃度最大值與最小值之差分別為0.39、1.00和1.53 mg/L，這是由灌溉水在長毛管滴灌單元中運移時間（水齡）及運移距離差異較大造成的。

圖6給出了不同加氯濃度條件下灌水單元平均余氯濃度和平均水齡、毛管末端余氯濃度和水齡隨毛管長度增加的變化。由圖可知，平均水齡和毛管末端水齡均隨毛管長度增加呈先降低而后升高的趨勢。毛管長度較短時，較大的平均水齡及毛管末端水齡與灌溉單元控制面積較小而導致的干、支管流速較小有關，這與圖5a中較大的毛管首端水齡一致。當毛管長度由10 m逐漸增加到40 m時，干、支管流速逐漸增大，平均水齡下降至最小值；隨著毛管長度的進一步增加，平均水齡和毛管末端水齡持續(xù)增加。

與平均水齡和毛管末端水齡變化趨勢相反，不同加氯濃度系統(tǒng)平均余氯濃度與毛管末端余氯濃度隨毛管長度增加呈先升高而后降低趨勢，水齡變化導致的有效氯與灌溉水及管壁反應時間的變化是引起這一變化的主要原因。不同加氯濃度條件下，在平均余氯濃度隨毛管長度增加至最大值后，隨毛管長度進一步增加，平均余氯濃度和毛管末端余氯濃度隨毛管長度增加衰減比例均隨加氯濃度的增加而減低。例如，對于加氯濃度為1.3 mg/L系統(tǒng)，當毛管長度由40 m逐漸增加到150 m時，平均余氯濃度與加氯濃度比值從82%逐漸降低至60%，而對于8 mg/L系統(tǒng)對應值則僅從91%降低至82%，這與圖5中不同加氯濃度處理沿毛管方向余氯衰減規(guī)律類似，也與低加氯濃度時較大的主體水衰減系數(shù)和管壁衰減系數(shù)有關，說明較高加氯濃度條件下平均余氯濃度對毛管長度變化的敏感性相對較弱。由圖6b還可看出，毛管末端余氯濃度隨毛管長度增加衰減速率大于平均余氯濃度。當加氯濃度為1.3、3、6和8 mg/L時，150 m毛管末端余氯濃度分別為平均余氯濃度的46%、61%、66%和70%。毛管末端余氯濃度、平均余氯濃度和加氯濃度的明顯差異給基于毛管末端余氯濃度控制系統(tǒng)加氯濃度、進而防止灌水器堵塞造成了一定困難。

圖6 不同加氯濃度灌水單元平均余氯濃度和平均水齡、毛管末端余氯濃度和水齡隨毛管長度增加的變化

3.2.2 加氯模式及毛管長度對余氯質量分布特征的影響

圖7給出了不同加氯模式及毛管長度條件下灌水單元內灌水器流量和余氯質量偏差率和均勻系數(shù)變化。灌水單元內灌水器流量偏差率隨毛管長度增加而增加，當毛管由10 m逐漸增加至150 m時，流量偏差率由0.01%增加到23%，以《微灌工程技術規(guī)范》（GB/T 50485-2009）[25]中規(guī)定灌水單元灌水器設計允許流量偏差率不超過20%為約束條件，僅考慮灌水均勻性情況下本研究所用滴灌帶極限鋪設長度為136 m，灌水器流量均勻系數(shù)對應值為95%。與灌水器流量偏差率呈現(xiàn)規(guī)律類似，灌水單元內余氯質量偏差率基本隨毛管長度的增加而升高，但所有加氯模式條件下余氯質量偏差率均明顯高于灌水器流量偏差率。例如，在毛管長度為10 m時，對于C1T4、C2T3、C3T2和C4T1加氯模式處理余氯質量偏差率分別達到22%、24%、36%和45%。灌水單元內余氯濃度差異是造成余氯質量偏差率較高的重要原因。除此之外，灌水單元內不同灌水器加氯歷時差異也是造成這一結果的重要原因。例如，由圖5a可知，毛管長度為10 m時，受管網中各灌水器節(jié)點水齡不同的影響，對于C4T1模式，自灌水結束前30 分鐘開始加氯后，毛管首端和毛管末端加氯歷時平均約為24.2和17.2 min。

圖7 不同加氯模式及毛管長度條件下灌水器流量和余氯質量偏差率和均勻系數(shù)

對比不同加氯模式余氯質量偏差率和均勻系數(shù)可知，C2T3加氯模式可以獲得相對較低的余氯質量偏差率和較高的余氯質量均勻系數(shù)。對于C3T2和C4T1加氯模式，較短的加氯時長導致的灌水器加氯歷時相對差異較大是導致高濃度加氯模式余氯質量偏差率較高的主要原因。例如，80 m毛管條件下，C3T2模式加氯歷時極差值為15 min，占到系統(tǒng)加氯時長的39%，遠高于C2T3模式加氯歷時極差占系統(tǒng)加氯時長比值（19%），這說明系統(tǒng)加氯過程中不宜設置過短的加氯時長，以保證灌水單元內不同位置灌水器加氯歷時相對差異較小。由圖7a還可看出，C1T4加氯模式余氯質量偏差率在不同毛管長度時基本高于C2T3模式，這與低濃度加氯條件下有效氯在管網中的明顯衰減有關。綜上可知，低濃度-長歷時加氯盡管可以降低灌水單元內灌水器加氯歷時相對差異，但是由于加氯濃度過小，有效氯在滴灌系統(tǒng)明顯衰減也會導致余氯質量偏差率的升高。

參照《微灌工程技術規(guī)范》（GB/T 50485-2009）中灌水器設計允許流量偏差率的規(guī)定，本研究假定余氯質量偏差率約束條件為40%，由此可推知C1T4、C2T3、C3T2和C4T1模式毛管極限鋪設長度分別為55、66、37和20 m，上述極限長度條件下對應余氯質量均勻系數(shù)分別為92%、93%、91%和92%，超過90%，說明C2T3模式在保障較高的余氯分布均勻性情況下可以明顯提高毛管鋪設長度。

4 結 論

本研究基于EPANET軟件，構建了再生水滴灌系統(tǒng)水力性能及余氯運移、衰減和分布模型，基于試驗數(shù)據(jù)對模型水力和水質模擬參數(shù)進行率定和驗證，并應用該模型分析了加氯模式（濃度-時長：1.3 mg/L-185 min、3 mg/L-80 min、6 mg/L-40 min和8 mg/L-30 min）和毛管長度（10～150 m）對系統(tǒng)余氯分布特征的影響，主要結論如下：

1）加氯濃度對氯衰減過程影響明顯，氯衰減速度隨隨加氯濃度的增加而降低。

2）基于EPANET軟件建立的再生水滴灌管網模型能較好模擬滴灌管網水力特征及系統(tǒng)余氯分布特征，沿毛管方向灌水器余氯濃度模擬值與實測值均方根誤差（RMSE）小于0.26 mg/L，一致性指數(shù)達到0.75以上。

3）灌水單元平均余氯濃度與毛管末端余氯濃度隨毛管長度增加呈先升高而后降低趨勢，高濃度加氯處理平均余氯濃度和毛管末端余氯濃度隨毛管長度增加衰減速率相對較小。

4）較短的加氯歷時（<40 min）和過低的加氯濃度（如1.3 mg/L）均易造成余氯質量偏差率升高和均勻系數(shù)降低；建議使用加氯濃度3 mg/L-加氯時長80 min的加氯模式，該模式毛管極限鋪設長度達到66 m，余氯質量均勻系數(shù)超過90%。

[1]Wang Zhen, Li Jiusheng, Li Yanfeng. Using reclaimed water for agricultural and landscape irrigation in China: A review[J]. Irrigation and Drainage, 2017, 66: 672-686.

[2]李云開，宋鵬，周博. 再生水滴灌系統(tǒng)灌水器堵塞的微生物學機理及控制方法研究[J]. 農業(yè)工程學報，2013，29(15)：98-107.

Li Yunkai, Song Peng, Zhou Bo. Microbiology mechanism and controlling methods for emitter clogging in the reclaimed water drip irrigation system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(15): 98-107. (in Chinese with English abstract)

[3]Ravina I, Paz E, Sofer Z, et al. Control of emitter clogging in drip irrigation with stored treated municipal sewage effluent[J]. Agricultural Water Management, 1997, 33(2/3): 127-137.

[4]Cararo D, Botrel T, Hills D, et al. Analysis of clogging in drip emitters during wastewater irrigation[J]. Applied Engineering in Agriculture, 2006, 22(2): 251-257.

[5]Li Jiusheng, Chen Lei, Li Yanfeng, et al. Effects of chlorination schemes on clogging in drip emitters during application of sewage effluent[J]. Applied Engineering in Agriculture, 2010, 26(4): 565-578.

[6]Song Peng, Li Yunkai, Zhou Bo, et al. Controlling mechanism of chlorination on emitter bio-clogging for drip irrigation using reclaimed water[J]. Agricultural Water Management, 2017, 184: 36-45.

[7]Song Peng, Feng Gi, Brooks J, et al. Environmental risk of chlorine-controlled clogging in drip irrigation system using reclaimed water: The perspective of soil health[J]. Journal of Cleaner Production, 2019, 232: 1452-1464.

[8]Li Jiusheng, Li Yanfeng, Zhang Hang. Tomato yield and quality and emitter clogging as effected by chlorination schemes of drip irrigation systems applying sewage effluent[J]. Journal of Integrative Agricultural, 2012, 11(10): 1744-1754.

[9]Hao Fengzhen, Li Jiusheng, Wang Zhen, et al. Influence of chlorine injection on soil enzyme activities and maize growth under drip irrigation with secondary sewage effluent[J]. Irrigation Science, 2018, 36(6): 363-379.

[10]Maier S H, Powell R S, Woodward C A. Calibration and comparison of chlorine decay models for a test water distribution system[J]. Water Research, 2000, 34(8): 2301-2309.

[11]Kang Yuehu, Nishiyama S. Finite element method analysis of microirrigation system pressure distribution[J]. Transactions of the Japanese Society of Irrigation, Drainage and Reclamation Engineering, 1994, 169: 19-26.

[12]Baiamonte G. Advances in designing drip irrigation laterals[J]. Agricultural Water Management, 2018, 199: 157-174.

[13]李援農，朱鋒. 基于最大控制面積和最低費用的微灌小區(qū)管網優(yōu)化[J]. 農業(yè)工程學報，2015，31(23)：80-87.

Li Yuannong, Zhu Feng. Optimization of basic irrigation controlling area network for micro-irrigation based on maximum control area and minimum cost models[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(23): 80-87. (in Chinese with English abstract)

[14]孫傅，陳吉寧，曾思育. 基于EPANET-MSX的多組分給水管網水質模型的開發(fā)與應用[J]. 環(huán)境科學，2008，29(12)：3360-3367.

Sun Fu, Chen Jining, Zeng Siyu. Development and application of a multi-species water quality model for water distribution systems with EPANET-MSX[J]. Environmental Science, 2008, 29(12): 3360-3367. (in Chinese with English abstract)

[15]Mohapatra S, Sargaonkar A, Labhasetwar P. Distribution network assessment using EPANET for intermittent and continuous water supply[J]. Water Resources Management, 2014, 28(11): 3745-3759.

[16]Lamm F R, Ayars J E, Nakayama F S. Microirrigation for Crop Production: Design, Operation and Management[M]. Amsterdam: Elsevier, 2007.

[17]Rossman L A. EPANET 2 User’s Manual[M]. Cincinnati: US Environmental Protection Agency, 2000.

[18]Todini E, Pilati S. A gradient method for the analysis of pipe networks[C]//International Conference on Computer Applications for Water Supply and Distribution. UK, Leicester Polytechnic, 1987: 1-20.

[19]環(huán)境保護部. HJ 586-2010：水質游離氯和總氯的測定[S]. 北京：中國環(huán)境科學出版社，2010.

[20]中華人民共和國國家質量監(jiān)督檢驗檢疫總局，中國國家標準化管理委員會. GB 19106-2013：次氯酸鈉[S]. 北京：中國標準出版社，2013.

[21]Hallam N B, West J R, Forster C F, et al. The decay of chlorine associated with the pipe wall in water distribution systems[J]. Water Research, 2002, 36(14): 3479-3488.

[22]ASAE. S553: Collapsible emitting hose (drip tape): Specifications and performance testing[S]. St Joseph, MI: ASABE, 2003.

[23]Wang Zhen, Li Jiusheng, Li Yanfeng. Simulation of nitrate leaching under varying drip system uniformities and precipitation patterns during the growing season of maize in the North China Plain[J]. Agricultural Water Management, 2014, 142: 19-28.

[24]Charles N Haas, Sankaram B Karra. Kinetics of microbial inactivation by chlorine: I review of results in demand-free systems[J]. Water Research, 1984, 18(11): 1443-1449.

[25]中華人民共和國住房和城鄉(xiāng)建設部，中華人民共和國國家質量監(jiān)督檢驗檢疫總局. GB/T 50485-2009 微灌工程技術規(guī)范[S]. 北京：中國計劃出版社，2009.

Modeling of the residual chlorine distribution in drip irrigation system applying reclaimed water based on EPANET

Wang Zhen1, Hao Fengzhen1,2, Li Jiusheng1※, Li Yanfeng1

(1.,100048,;2.030801,)

Chlorination is an important practice to ensure the safety of drip irrigation system with reclaimed water application. However, the decaying process caused by chlorine reacting with irrigation water and drip irrigation pipeline will adversely affect the uniformity of residual chlorine distribution. The objective of this study was to evaluate the effect of different chlorination modes on the residual chlorine distribution in drip irrigation system using a model, which was constructed using the EPANET software package and could simulate the hydraulic performance and residual chlorine transport, attenuation and distribution in a drip irrigation system. The hydraulic and water quality parameters for the model were calibrated and validated based on the experimental data. The distribution of residual chlorine under 4 chlorination modes (chlorination concentration injected– chlorination duration of 1.3 mg/L-185 min (C1T4), 3 mg/L-80 min (C2T3), 6 mg/L-40 min (C3T2), and 8 mg/L-30 min (C4T1)) and 15 lateral lengths (10-150 m at 10-m intervals) were analyzed. The results showed that the hydraulic performance and residual chlorine distribution in the drip irrigation system were well simulated by the constructed model using the EPANET software. The index of agreement between the simulated and measured residual chlorine concentration of emitter discharge along lateral is above 0.75. The age of water discharged from emitters increased with the distance from the inlet of the lateral, while it ranged from 5.8-12.8, 1.2-16.9, and 1.1-27.3 min for the system with lateral length of 10, 80, and 150 m, respectively. The residual chlorination concentration decreased gradually from the inlet of the lateral due to the reaction of chlorine with water and pipe wall. The decay rate of the residual chlorine concentration decreased with the chlorine concentration injected. The relatively higher decay rate was attributed to the higher bulk reaction coefficient and wall reaction coefficient in the drip irrigation system with low chlorination concentration. Under a given injected chlorination concentration, the system with longer length of lateral produced higher variation of the residual chlorine concentration because of the relatively higher variation of water age. The average residual chlorine concentration of the drip irrigation system increased first and then decreased with increasing length of the lateral. A relatively lower decay rate of the average residual chlorine concentration and the residual chlorine concentration at the end of the lateral was observed for drip irrigation system with high injected chlorination concentration. For example, the ratio of average residual chlorine concentration to chlorination concentration injected decreased from 82% to 60% and from 91% to 82% when the lateral increased from 40 to 150 m for the system with 1.3 and 8.0 mg/L of chlorination concentration, respectively. Much higher deviation of the residual chlorine mass was observed than the water applied for all the chlorination modes. The deviations of the residual chlorine mass for the system with 10-m length of lateral were 22%, 24%, 36%, and 45% for the C1T4, C2T3, C3T2 and C4T1 modes, respectively. The C2T3 mode generally produced the lowest deviation and the highest Christiansen uniformity coefficient (CU) of residual chlorine mass while the length of lateral ranged from 10 to 150 m. A relatively shorter duration of chlorination (smaller than 40 min) and too low injected chlorination concentration (1.3 mg/L) probably produced higher deviation of the residual chlorine mass and thus decreased the coefficient of uniformity of the residual chlorine mass. The injected chlorination concentration of 3 mg/L with chlorination duration of 80 min was recommended, while the maximum lateral length could extend to 66 m with the coefficient of uniformity of residual chlorine mass exceeding 90%.

chlorine; simulation; hydraulic performance; water quality; drip irrigation; EPANET

王珍，郝鋒珍，李久生，等. 基于EPANET的再生水滴灌系統(tǒng)余氯分布模型構建[J]. 農業(yè)工程學報，2020，36(10)：99-106.doi：10.11975/j.issn.1002-6819.2020.10.012 http://www.tcsae.org

Wang Zhen, Hao Fengzhen, Li Jiusheng, et al. Modeling of the residual chlorine distribution in drip irrigation system applying reclaimed water based on EPANET[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(10): 99-106. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.10.012 http://www.tcsae.org

2019-09-15

2020-04-10

國家重點研發(fā)計劃課題（2016YFC0403103）；國家自然科學基金項目（51790531）；流域水循環(huán)模擬與調控國家重點實驗室自由探索課題（SKL2018TS06）

王珍，高級工程師，主要從事灌溉原理及技術研究。Email：wangzhen-686@163.com

李久生，研究員，博士生導師，主要從事灌溉原理與技術研究。Email：lijs@iwhr.com

10.11975/j.issn.1002-6819.2020.10.012

S275.6

A

1002-6819(2020)-10-0099-08