洪建權(quán)　鄒家榮　丁世洪　張志秀　賈忠華

摘要：QUAL2K水質(zhì)模型可以模擬污染物在排水通道中的降解過(guò)程，分析排水溝塘對(duì)農(nóng)業(yè)污染物的凈化效果。以江蘇省揚(yáng)州市江都區(qū)水稻灌區(qū)一典型排水系統(tǒng)為例，利用QUAL2K模型模擬農(nóng)田溝塘排水過(guò)程中氨氮和總磷兩種污染物的削減情況;在利用實(shí)地調(diào)查與監(jiān)測(cè)數(shù)據(jù)進(jìn)行模型檢驗(yàn)以后，根據(jù)實(shí)際水力聯(lián)系情況，預(yù)測(cè)分析了調(diào)整水力聯(lián)系及采用部分重點(diǎn)治理措施對(duì)改善整體水質(zhì)凈化效果的作用。結(jié)果顯示：① 研究區(qū)6個(gè)主要水質(zhì)監(jiān)測(cè)點(diǎn)模擬的平均相對(duì)誤差均在20%左右，相關(guān)系數(shù)R2和Nash-Sutcliffe系數(shù)NSE的計(jì)算值均大于0.5，驗(yàn)證了模型的可靠性，但是田間溝塘一些人為及自然等不確定因素對(duì)模型的精度具有一定的影響。② 3種治理措施的效果為截污>減小排水流量>調(diào)整溝塘流網(wǎng)系統(tǒng)，模擬顯示排水系統(tǒng)經(jīng)過(guò)少數(shù)可行的水力關(guān)系調(diào)整，就可以顯著提高整體水質(zhì)改善效果。QUAL2K模型可以靈活調(diào)整排水溝塘水力聯(lián)系，預(yù)演不同治理措施，為優(yōu)化農(nóng)田排水溝塘系統(tǒng)提供依據(jù)。

關(guān)鍵詞：農(nóng)田排水溝塘; QUAL2K模型; 水質(zhì)凈化; 農(nóng)業(yè)污染

中圖法分類號(hào)： S276.3

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2021.11.009

0引 言

隨著農(nóng)業(yè)生產(chǎn)過(guò)程中化肥農(nóng)藥的使用量不斷增加，農(nóng)田排水帶來(lái)的水環(huán)境惡化等問(wèn)題也越發(fā)突出，農(nóng)業(yè)非點(diǎn)源污染治理成為我國(guó)環(huán)境治理工作的一項(xiàng)重要內(nèi)容[1]。農(nóng)業(yè)非點(diǎn)源污染的產(chǎn)生一般是在降雨、灌溉的過(guò)程中，化肥和農(nóng)藥等污染物隨地表、地下徑流進(jìn)入農(nóng)田排水系統(tǒng)[2]。這不僅會(huì)造成農(nóng)田的養(yǎng)分流失嚴(yán)重，而且會(huì)導(dǎo)致接納水體污染物的超標(biāo)，引起水體富營(yíng)養(yǎng)化等一系列水環(huán)境污染問(wèn)題。農(nóng)業(yè)非點(diǎn)源污染具有污染源比較分散、涉及范圍很廣、控制難度較大等特點(diǎn)，是目前農(nóng)業(yè)水環(huán)境治理問(wèn)題的瓶頸[3]。農(nóng)業(yè)排水溝塘系統(tǒng)作為一種連接農(nóng)田與接納水體的天然緩沖帶，具有類似于濕地的功能，可以通過(guò)一系列生物化學(xué)過(guò)程去除農(nóng)田排水中的化肥、農(nóng)藥等污染物。處在排水通道上的溝塘具備去除污染物的便利條件，是一種天然的、不可替代的生態(tài)資源[4]。

研究發(fā)現(xiàn)農(nóng)田溝塘濕地系統(tǒng)對(duì)污染物的去除效果與溝塘在田間的分布特性密切相關(guān)，很多情況下，不利的分布和水力條件會(huì)影響溝塘系統(tǒng)的整體去污效率[5-6]。利用水質(zhì)模型不僅能夠模擬預(yù)測(cè)不同氣象和農(nóng)田灌排布置條件下復(fù)雜溝塘系統(tǒng)的去污能力，還能夠探求優(yōu)化溝塘水力聯(lián)系。其中，河流水質(zhì)模型可用來(lái)描述河流污染物的遷移和轉(zhuǎn)化過(guò)程及其它們之間的復(fù)雜聯(lián)系，預(yù)測(cè)河流水質(zhì)的變化情況[7]。由美國(guó)國(guó)家環(huán)保局推出的QUAL系列河流水質(zhì)模型經(jīng)歷了4個(gè)階段的開發(fā)研究，最新的、比較完善的QUAL2K水質(zhì)模型在國(guó)內(nèi)外的應(yīng)用比較廣泛，其功能全面、通用性強(qiáng)，能夠全面準(zhǔn)確地反映污染物在不同水體中的轉(zhuǎn)化規(guī)律[8-9]。QUAL2K模型在國(guó)內(nèi)主要將其用于模擬一維、穩(wěn)態(tài)的中小型河流水質(zhì)變化情況。如馬丁等[10]應(yīng)用QUAL2K模型對(duì)四川省某中小流域總磷的沿程變化進(jìn)行了模擬和預(yù)測(cè)，優(yōu)化了4種水污染治理措施，從而驗(yàn)證了QUAL2K模型在中小型河流水質(zhì)模擬中較高的模擬精度。劉瑾等[11]利用QUAL2K模型模擬渾河流域氨氮的變化情況，得到了氨氮在不同水期的變化規(guī)律，從而預(yù)測(cè)了該區(qū)域氨氮量。農(nóng)田溝塘排水系統(tǒng)類似于枝狀復(fù)雜河道水流匯集系統(tǒng)，但現(xiàn)有研究中對(duì)于應(yīng)用相關(guān)水質(zhì)模型模擬農(nóng)田排水污染物的變化規(guī)律研究較少。本文以江蘇省揚(yáng)州市江都區(qū)京杭大運(yùn)河?xùn)|側(cè)沿運(yùn)灌區(qū)的一個(gè)農(nóng)田排水溝塘系統(tǒng)為研究對(duì)象，選取污染物氨氮、總磷為控制因子，應(yīng)用QUAL2K模型模擬農(nóng)田排水過(guò)程中兩種因子在排水溝水系內(nèi)沿程的變化情況，在檢驗(yàn)?zāi)Ｐ瓦m用性的基礎(chǔ)上，探討優(yōu)化排水溝塘系統(tǒng)的水質(zhì)凈化能力，研究結(jié)果可為類似農(nóng)田溝塘系統(tǒng)污染物治理提供借鑒。

1材料與方法

1.1QUAL2K模型介紹

QUAL2K模型是美國(guó)國(guó)家環(huán)保局推出的一個(gè)一維穩(wěn)態(tài)綜合性河流水質(zhì)模型，適用于枝狀復(fù)雜河流網(wǎng)絡(luò)，允許沿河有多個(gè)取水口、排污口、支流分析入流點(diǎn)、面源負(fù)荷。該模型的水質(zhì)基本方程是一維平流-擴(kuò)散物質(zhì)遷移反應(yīng)方程[12-13]，考慮了平流彌散、稀釋、水質(zhì)組分的自反應(yīng)、水質(zhì)組分之間的相互作用以及組分外部源匯對(duì)組分濃度的影響[14-15]。

1.2研究區(qū)概況

研究區(qū)位于江蘇省揚(yáng)州市江都區(qū)京杭大運(yùn)河?xùn)|側(cè)的沿運(yùn)灌區(qū)（119°25′E、32°22′N），區(qū)內(nèi)地勢(shì)平坦，年平均溫度14.9 ℃，年降雨量約1 000 mm。研究區(qū)內(nèi)普遍實(shí)行稻麥輪作。圖1為研究區(qū)農(nóng)田與溝塘的分布及水力聯(lián)系。區(qū)內(nèi)農(nóng)田總面積為5.61 hm2，溝塘總面積為0.80 hm2，兩者的面積比為14.3%。研究區(qū)主要包括3個(gè)排水支路，但各支路控制的農(nóng)田面積差異較大，3條排水支路所占溝塘面積以及控制農(nóng)田面積如表1所列。支路1流程較長(zhǎng)，包括6個(gè)階段;其余2個(gè)支路的流程都很短，農(nóng)田排水直接排入附近的池塘和支溝。因此，本文選取支路1進(jìn)行QUAL2K模型的模擬研究。

根據(jù)實(shí)地調(diào)查監(jiān)測(cè)，基于QUAL2K模型的原理，對(duì)支路1中的13個(gè)排水溝單元進(jìn)行概化，其中排水溝分段的原則及水流匯入情況如下：

（1） 根據(jù)水質(zhì)模擬需要，首先將溝道劃分成一系列恒定的非均勻流溝段（同一溝段要求具有相同的水力參數(shù)和水質(zhì)特征），然后將劃分好的各個(gè)溝段再劃分為若干個(gè)等長(zhǎng)的計(jì)算單元[16-17]。

（2） 研究區(qū)模擬的排水溝主干段全長(zhǎng)500 m，分為5段，其中2條支流匯入，長(zhǎng)度分別為300 m和200 m，各劃分為3段和2段，共劃分10段。由于農(nóng)田排水溝排水距離短，排水范圍小，故在每個(gè)溝段設(shè)置2個(gè)計(jì)算單元;另外還有3處點(diǎn)源匯入。研究區(qū)溝段劃分及污染源匯入情況如圖3所示。

1.3QUAL2K模型參數(shù)率定

模型涉及的主要參數(shù)為研究區(qū)排水溝塘的水力參數(shù)和水質(zhì)參數(shù)，其中，實(shí)測(cè)輸入數(shù)據(jù)包括排水溝氨氮、總磷濃度均值，氣象數(shù)據(jù)以及排水溝流量均值等。本文將2017年研究區(qū)稻作期間水質(zhì)監(jiān)測(cè)數(shù)據(jù)用于模型的率定，2018年稻作期間水質(zhì)監(jiān)測(cè)數(shù)據(jù)用于模型的驗(yàn)證。

主要水質(zhì)參數(shù)包括有機(jī)氮水解系數(shù)khn、氨氮硝化系數(shù)kna、有機(jī)磷水解系數(shù)khp等，其他參數(shù)均采用模型推薦值。基于文獻(xiàn)報(bào)道的適宜參數(shù)取值范圍[18]，對(duì)上述參數(shù)進(jìn)行合理調(diào)試，直至得到滿意的擬合結(jié)果。主要水質(zhì)參數(shù)取值如表2所列。

QUAL2K模型提供了3種方法來(lái)計(jì)算各溝段水力學(xué)特征，分別是溢流堰法、流量系數(shù)法和曼寧系數(shù)法。本文采用了曼寧系數(shù)法，根據(jù)查閱相關(guān)資料，研究區(qū)較大的排水支溝取值0.04，排水農(nóng)溝取值0.02。率定后得到的水力學(xué)特征參數(shù)列于表3。

2結(jié)果與分析

2.1QUAL2K模型模擬及驗(yàn)證

基于上述水質(zhì)參數(shù)和水力參數(shù)，結(jié)合實(shí)際監(jiān)測(cè)數(shù)據(jù)進(jìn)行了模型的驗(yàn)證，2018年稻作期間排水溝平均流量為0.001 5 m3/s，圖4為排水溝流量過(guò)程圖。研究區(qū)稻作期間農(nóng)田排水溝污染物總磷、氨氮濃度實(shí)測(cè)值與QUAL2K模擬值的比較結(jié)果如表4～5所列。

數(shù)據(jù)顯示，研究區(qū)監(jiān)測(cè)點(diǎn)G2、G3、G4、G5、G6的QUAL2K總磷模擬值與實(shí)測(cè)值間平均相對(duì)誤差為6.72%;誤差最小值位于排水溝長(zhǎng)0.500 km的G6監(jiān)測(cè)點(diǎn)處，最小相對(duì)誤差為2.05%;誤差最大值位于排水溝長(zhǎng)0.405 km的G5監(jiān)測(cè)點(diǎn)處，最大相對(duì)誤差為13.94%。

研究區(qū)監(jiān)測(cè)點(diǎn)G2、G3、G4、G5、G6的QUAL2K氨氮模擬值與實(shí)測(cè)值的平均相對(duì)誤差為55.48%;誤差最小值位于排水溝長(zhǎng)0.500 km處的G6監(jiān)測(cè)點(diǎn)處，相對(duì)誤差為12.65%，誤差最大值位于排水溝長(zhǎng)0.405 km的G5監(jiān)測(cè)點(diǎn)處，相對(duì)誤差為185.15%。若考慮實(shí)驗(yàn)監(jiān)測(cè)過(guò)程中的過(guò)失誤差，去除最大相對(duì)誤差為185.15%的點(diǎn)，平均相對(duì)誤差為23.06%。

QUAL2K模型對(duì)平原河網(wǎng)地區(qū)農(nóng)田排水溝塘總磷、氨氮模擬結(jié)果的平均相對(duì)誤差均在20%左右，說(shuō)明QUAL2K模型能夠較好地對(duì)該研究區(qū)農(nóng)田排水溝塘水質(zhì)狀況進(jìn)行模擬。表6列出了氨氮、總磷模擬精度評(píng)價(jià)結(jié)果，表中模擬值與實(shí)測(cè)值之間的相關(guān)系數(shù)R2都在90%以上。模型驗(yàn)證期R2和NSE值大多大于0.5，R2值在0.9以上，模擬值和實(shí)測(cè)值擬合程度較好，表明QUAL2K模型能夠較好地模擬農(nóng)田溝塘排水過(guò)程中氨氮和總磷污染物的運(yùn)移削減規(guī)律。

2.2排水溝塘系統(tǒng)不同優(yōu)化措施提升水質(zhì)凈化效果模擬

為了探討提高研究區(qū)溝塘下游斷面污染物去除效果的優(yōu)化方案，根據(jù)以上模擬分析，結(jié)合排水系統(tǒng)現(xiàn)狀，針對(duì)排水溝塘水力條件提出以下3種情景設(shè)置進(jìn)行不同優(yōu)化方案的模擬（見表7）。

2.2.1調(diào)整水力聯(lián)系對(duì)溝塘水質(zhì)凈化效果的影響

南方平原地區(qū)排水溝塘交錯(cuò)分布，導(dǎo)致不同田塊的排水從多點(diǎn)、多級(jí)進(jìn)入/排出溝塘，幾乎不存在單一的“出口-入口”的關(guān)系[20-21]。目前的排水溝塘分布和水力聯(lián)系受到自然、歷史和農(nóng)業(yè)生產(chǎn)等多種因素的影響，分布上存在一定的隨機(jī)性[22]。溝塘與農(nóng)田面積匹配可能不合理，大面積農(nóng)田的附近可能分布著較小的溝塘，或是小面積農(nóng)田的附近分布著較大的溝塘，影響排水溝塘的水質(zhì)凈化功能。為了分析排水溝塘分布現(xiàn)狀對(duì)排水水質(zhì)凈化效果的影響，以研究區(qū)溝塘系統(tǒng)為例，在不改變現(xiàn)有溝塘總體分布的情況下，利用QUAL2K模型進(jìn)行排水系統(tǒng)流網(wǎng)優(yōu)化（見圖5），目的是在水力條件允許的情況下，發(fā)揮尺寸較大排水支溝的作用。

圖6～7為各情景下總磷、氨氮控制效果。

由圖6～7可知，情景1通過(guò)調(diào)整水力聯(lián)系，對(duì)排水系統(tǒng)流網(wǎng)進(jìn)行優(yōu)化后，在一定程度上改善了排水溝的水質(zhì)。下游斷面總磷質(zhì)量濃度從優(yōu)化前的146? μg/L降至136 μg/L，降幅為6.85%。下游斷面氨氮質(zhì)量濃度從優(yōu)化前的253 μg/L降至214 μg/L，降幅為15.4%，

可知優(yōu)化效果優(yōu)于總磷。從排水過(guò)程上來(lái)看，

排水主干溝前段部分沿程氨氮、總磷濃度下降較快，可見，較大支溝對(duì)于污染物的去除達(dá)到了較好的效果。

2.2.2排水溝入流流量變化對(duì)溝塘水質(zhì)凈化效果的影響

農(nóng)田排水溝塘在作物種植期間流量變化大，由于水稻種植需水量大，灌溉用水多，灌溉期排水流量大，排水速度過(guò)快，導(dǎo)致水力停留時(shí)間短，對(duì)氮磷污染物去除效果差。為保證排水溝對(duì)污染物的去除效果，可以通過(guò)控制田間排水技術(shù)，提高排水溝排水深度，減小排水流量，從而提高水力停留時(shí)間。在利用QUAL2K模型進(jìn)行模擬時(shí)，使排水溝流量保持在0.001 m3/s，以期達(dá)到類似的效果。此情景實(shí)施下，排水溝下游斷面總磷質(zhì)量濃度減至70 μg/L，同比下降52.1%，氨氮質(zhì)量濃度減至206 μg/L，同比下降18.6%（見圖6～7）。

2.2.3排水溝截污治理對(duì)溝塘水質(zhì)凈化效果的影響

研究區(qū)農(nóng)田排水集中分布在排水農(nóng)溝，經(jīng)統(tǒng)計(jì)，2018年研究區(qū)某些排水農(nóng)溝農(nóng)田排水氮磷入溝質(zhì)量濃度過(guò)高，排水溝取樣時(shí)明顯看出水樣顏色偏深，營(yíng)養(yǎng)物質(zhì)含量豐富，因此，未來(lái)可對(duì)營(yíng)養(yǎng)物質(zhì)嚴(yán)重聚集的排水溝進(jìn)行截污整治。在情景3實(shí)施后，排水溝水質(zhì)得到很大程度的提升，排水溝下游斷面總磷質(zhì)量濃度減少至57 μg/L，下降了約61%;氨氮質(zhì)量濃度減少至182 μg/L，下降了28.1%（見圖6～7）。

3討 論

研究結(jié)果表明：3種情景的排水溝塘水力條件優(yōu)化措施均可提高水質(zhì)凈化能力，其水質(zhì)凈化效果排序?yàn)榻匚?減小排水流量>調(diào)整溝塘流網(wǎng)系統(tǒng)。農(nóng)田排水溝塘系統(tǒng)，作為去除氮磷污染物的最直接的凈化場(chǎng)所，關(guān)于它的研究也不勝枚舉，國(guó)內(nèi)學(xué)者主要集中在溝塘底泥、水生植物、微生物方向，通過(guò)研究溝塘土壤吸附、微生物降解、水生植物吸收等作用實(shí)現(xiàn)對(duì)氮磷污染物的去除。針對(duì)溝塘現(xiàn)存問(wèn)題，現(xiàn)有研究關(guān)于溝塘分布特性以及溝塘與農(nóng)田之間的水力聯(lián)系對(duì)于污染物去除關(guān)注不多。通過(guò)QUAL2K模型針對(duì)溝塘存在問(wèn)題進(jìn)行溝塘潛力優(yōu)化，具有一定的創(chuàng)新性。但是對(duì)模型的精度評(píng)價(jià)可知NSE值并不是那么的理想，由于農(nóng)田溝塘的不確定因素，包括一些人為及自然因素，對(duì)模型的模擬精度影響較大。后期需要定期觀察水位的變化，加大排水水質(zhì)的監(jiān)測(cè)力度，注重排水流量過(guò)程監(jiān)測(cè)，為模擬水質(zhì)因子的日均沿程變化提供更加可靠的數(shù)據(jù)。

4結(jié) 論

本文以揚(yáng)州江都區(qū)京杭大運(yùn)河?xùn)|側(cè)的沿運(yùn)灌區(qū)為例，基于QUAL2K建立了農(nóng)田排水溝的水質(zhì)模型，對(duì)排水溝氨氮、總磷兩種污染物的去除效果進(jìn)行了模擬驗(yàn)證，并利用率定后的模型對(duì)研究區(qū)排水系統(tǒng)進(jìn)行優(yōu)化方案預(yù)測(cè)分析，結(jié)果表明：

（1） 應(yīng)用QUALZK模型對(duì)農(nóng)田排水溝主干段6個(gè)監(jiān)測(cè)點(diǎn)的水質(zhì)指標(biāo)模擬，除去個(gè)別監(jiān)測(cè)點(diǎn)，平均相對(duì)誤差均在20%左右，相關(guān)系數(shù)R2和Nash-Sutcliffe系數(shù)NSE大都大于0.50，表明模擬值與實(shí)測(cè)值相關(guān)性較好，符合模型精度要求，表明QUAL2K模型能夠較好地模擬農(nóng)田溝塘排水過(guò)程中污染物氨氮和總磷在排水溝中的運(yùn)移削減規(guī)律。

（2） 應(yīng)用QUAL2K模型分3種情景進(jìn)行模擬，對(duì)研究區(qū)溝塘水質(zhì)改善進(jìn)行了優(yōu)化分析，3種情景的實(shí)施效果為截污>減小排水流量>調(diào)整溝塘流網(wǎng)系統(tǒng)。

參考文獻(xiàn)：

[1]李秀芬，朱金兆，顧曉君，等.農(nóng)業(yè)面源污染現(xiàn)狀與防治進(jìn)展[J].中國(guó)人口·資源與環(huán)境，2010，20（4）：81-84.

[2]洪建權(quán)，馬勇，丁世洪，等.基于QUAL2K模型的農(nóng)田排水溝水質(zhì)模擬研究[J].江蘇水利，2020（4）：15-19.

[3]馬凡凡，邢素林，徐云連，等.巢湖圩區(qū)農(nóng)田排水溝渠氮磷時(shí)空分布及輸出特征[J].廣東農(nóng)業(yè)科學(xué)，2019，46（5）：67-73.

[4]賈忠華，陳誠(chéng)，羅紈，等.農(nóng)業(yè)排水溝塘系統(tǒng)污染物去除監(jiān)測(cè)區(qū)代表性分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2018，34（3）：110-117.

[5]孫少江，羅紈，鄒家榮，等.農(nóng)田排水系統(tǒng)分布特征對(duì)水質(zhì)凈化效果的影響研究[J].灌溉排水學(xué)報(bào)，2019，38（4）：122-128.

[6]羅紈，朱金城，賈忠華，等.排水溝塘分布特性及與農(nóng)田水力聯(lián)系對(duì)水質(zhì)凈化能力的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2017，33（10）：161-167.

[7]汪家權(quán)，陳眾，武君.河流水質(zhì)模型及其發(fā)展趨勢(shì)[J].安徽師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2004（6）：242-247.

[8]周華.河流綜合水質(zhì)模型QUAL2K應(yīng)用研究[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào)，2010，8（1）：71-75.

[9]張麗，吳金亮，楊國(guó)范，等.QUAL2K模型在蘇子河水質(zhì)模擬中的應(yīng)用[J].人民黃河，2013，35（12）：83-85.

[10]馬丁，姚建.中小流域總磷污染模擬及控制對(duì)策研究[J].環(huán)境污染與防治，2018，40（6）：728-732.

[11]劉瑾，田博，項(xiàng)學(xué)敏，等.渾河流域沈陽(yáng)段氨氮污染源核算及水質(zhì)模擬分析[J].水電能源科學(xué)，2020，38（1）：40-43.

[12]方曉波.錢塘江流域水環(huán)境承載能力研究[D].杭州：浙江大學(xué)，2009.

[13]PARK S S，LEE Y S.A water quality modeling study of the Nakdong River，Korea[J].Ecological Modelling，2002，152（1）：65-75.

[14]唐偉.基于QUAL2K模型的水質(zhì)模擬研究：以武進(jìn)港小流域?yàn)槔齕D].南京：南京大學(xué)，2011.

[15]PARK S S，YONG S L.A water quality modeling study of the Nakdong River，Korea[J].Ecological Modelling，2002，152（1）：65-75.

[16]楊海林，楊順生.河流綜合水質(zhì)模型QUAL2E在河流水質(zhì)模擬中的應(yīng)用[J].云南環(huán)境科學(xué)，2003，22（2）：22-25.

[17]翟敏婷，辛卓航，韓建旭，等.河流水質(zhì)模擬及污染源歸因分析[J].中國(guó)環(huán)境科學(xué)，2019，39（8）：3457-3464.

[18]趙琰鑫，陳巖，吳悅穎.QUAL2K河流水質(zhì)模擬模型理論方法與應(yīng)用指南[M].北京：氣象出版社，2015.

[19]耿潤(rùn)哲，王曉燕，龐樹江，等.潮河流域非點(diǎn)源污染控制關(guān)鍵因子識(shí)別及分區(qū)[J].中國(guó)環(huán)境科學(xué)，2016，36（4）：1258-1267.

[20]GREGOIRE C，ELSAESSER D，HUGUENOT D，et al.Mitigation of agricultural nonpoint-source pesticide pollution in artificial wetland ecosystems[J].Environmental Chemistry Letters，2009，7（3）：205-231.

[21]李強(qiáng)坤，宋常吉，胡亞偉，等.模擬排水溝渠非點(diǎn)源溶質(zhì)氮遷移實(shí)驗(yàn)研究[J].環(huán)境科學(xué)，2016，37（2）：520-526.

[22]JIA Z，WU Z，LUO W，et al.The impact of improving irrigation efficiency on wetland distribution in an agricultural landscape in the upper reaches of the Yellow River in China[J].Agricultural Water Management，2013，121（4）：54-61.

（編輯：劉 媛）

Abstract：QUAL2K water quality model can simulate the degradation process of pollutants in the drainage channels，and analyze the purification effect of drainage ditches and ponds on agricultural pollutants.Taking a typical drainage system of a rice irrigation area in Jiangdu District of Yangzhou City，Jiangsu Province as an example，we simulated the reduction of two pollutants of ammonia nitrogen and total phosphorus in the process of farmland drainage in ditches and ponds by QUAL2K model.After using the field survey and monitoring data for model test，according to the actual hydraulic connection，the effect of adjusting hydraulic connection and some key treatment measures on improving the overall water quality purification were predicted and analyzed.The results showed that： ① The average relative errors of simulated values of water quality indicators at six main monitoring points in the study area were all about 20%，and the calculated values of R2 and Nash Sutcliffe coefficient （NSE） were greater than 0.5，verifying the reliability of the model，but some uncertain factors such as man-made and natural factors had a certain impact on the accuracy of the model in the ditches and ponds of field.② The effect of the three treatment measures was "sewage interception>drainage flow reduction>ditch pond network adjustment".QUAL2K model can flexibly adjust the hydraulic connection of the drainage ditches and ponds，simualte different treatment measures，and provide basis for optimizing farmland drainage ditches and ponds system.

Key words：farmland drainage ditches and ponds;QUAL2K model;water purification;agricultural pollution