覃苑，胡海棠，淮賀舉，李存軍*，張巧玲，楊鐵利，王佳宇

1.遼寧科技大學(xué)

2.北京市農(nóng)林科學(xué)院信息技術(shù)研究中心

農(nóng)業(yè)面源污染（agricultural non-point source pollution，ANSP）指由于農(nóng)業(yè)生產(chǎn)中的氮磷營(yíng)養(yǎng)元素、農(nóng)藥及其他溶解態(tài)和固態(tài)污染物流失導(dǎo)致的污染。近年來(lái)，農(nóng)業(yè)面源污染對(duì)水體污染的貢獻(xiàn)越來(lái)越受到重視。2021年3月生態(tài)環(huán)境部和農(nóng)業(yè)農(nóng)村部聯(lián)合發(fā)布了《農(nóng)業(yè)面源污染治理與監(jiān)督指導(dǎo)實(shí)施方案（試行）》，明確了“十四五”至2035年農(nóng)業(yè)面源污染防治的總體要求、工作目標(biāo)，提出主要任務(wù)包括深入推進(jìn)農(nóng)業(yè)面源污染防治、完善農(nóng)業(yè)面源污染防治政策機(jī)制和加強(qiáng)農(nóng)業(yè)面源污染治理監(jiān)督管理。ANSP的污染負(fù)荷估算、關(guān)鍵源區(qū)源頭識(shí)別和治理決策是監(jiān)督管理的重要內(nèi)容。

ANSP在空間上分布分散且隱蔽，其形成受到區(qū)域的氣候、地形、土地利用類型和植被覆蓋等多方面因素的影響，導(dǎo)致ANSP在監(jiān)測(cè)、模擬和控制等方面面臨較大挑戰(zhàn)。GIS技術(shù)在空間數(shù)據(jù)管理、集成、分析方面有優(yōu)勢(shì)，已被廣泛應(yīng)用于ANSP的數(shù)據(jù)管理、空間分析、面源污染模型集成、信息系統(tǒng)開發(fā)和可視化等方面[1]。國(guó)內(nèi)外采用年度農(nóng)業(yè)面源污染模型（annualized agricultural nonpoint source pollution model，AnnAGNPS），水土評(píng)價(jià)模型（soil and water assessment tool，SWAT），輸出系數(shù)模型（export coefficient model，ECM）等模型和平均濃度法等方法進(jìn)行ANSP模擬與識(shí)別研究。如夏昊[2]結(jié)合AnnAGNPS模型，開發(fā)了基于GIS的面源信息系統(tǒng)，其擁有空間分析和可視化表達(dá)等功能，用于管理ANSP中不同來(lái)源的空間和屬性數(shù)據(jù)；喬衛(wèi)芳等[3]基于研究區(qū)域數(shù)字高程模型（digital elevation model，DEM）數(shù)據(jù)，將GIS與SWAT模型結(jié)合構(gòu)建了ANSP基礎(chǔ)信息庫(kù)，分析了研究區(qū)ANSP負(fù)荷的分布特征，并模擬了區(qū)域內(nèi)主要土地利用類型單位面積ANSP的負(fù)荷量。但AnnAGNPS和SWAT等模型所需參數(shù)多，計(jì)算復(fù)雜，在大區(qū)域的模擬上略顯乏力；生態(tài)系統(tǒng)服務(wù)和權(quán)衡的綜合評(píng)估（InVEST）模型中的營(yíng)養(yǎng)物傳輸率模型（NDR）所需參數(shù)少，機(jī)理清晰，適用于大區(qū)域的ANSP模擬研究[4]。

在進(jìn)行ANSP治理時(shí)，關(guān)鍵源區(qū)識(shí)別是面源污染管理的重要組成部分。研究表明，少數(shù)區(qū)域輸出的污染物通常是整個(gè)流域污染負(fù)荷的主要貢獻(xiàn)者，對(duì)水體環(huán)境質(zhì)量有著至關(guān)重要的作用[5]。《2015年中國(guó)環(huán)境狀況公報(bào)》顯示，海河、巢湖等流域的水質(zhì)狀況存在明顯的空間異質(zhì)性，表現(xiàn)為支流水質(zhì)在重度污染的情況下，其所屬的干流僅為輕度污染[6]。目前主要是針對(duì)面源污染負(fù)荷空間分布的異質(zhì)性特征進(jìn)行ANSP關(guān)鍵源區(qū)識(shí)別[7]，但是單純考慮污染負(fù)荷易忽略流域面積與產(chǎn)水量對(duì)污染負(fù)荷貢獻(xiàn)的影響，導(dǎo)致關(guān)鍵源區(qū)的識(shí)別準(zhǔn)確度下降，進(jìn)而影響治理效果。應(yīng)結(jié)合潛在氮、磷徑流濃度和入河負(fù)荷的空間分布進(jìn)行ANSP防控的關(guān)鍵源區(qū)識(shí)別與分析[8]。目前農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別涉及水文分析、農(nóng)業(yè)面源污染模型、柵格矢量運(yùn)算和統(tǒng)計(jì)，對(duì)于面源污染非專業(yè)人士和基層工作者存在跨平臺(tái)操作、數(shù)據(jù)類型多、資料匯總繁瑣、計(jì)算量大等問(wèn)題，亟需集成式的軟件工具平臺(tái)。

ANSP包括農(nóng)田種植業(yè)面源污染、畜禽養(yǎng)殖、農(nóng)村生活污染等，其中農(nóng)田種植業(yè)面源污染是農(nóng)業(yè)面源污染的重要組成部分，污染源相對(duì)分散且影響因素多。筆者針對(duì)在農(nóng)田種植業(yè)氮、磷化肥施用形成的面源污染關(guān)鍵源區(qū)識(shí)別過(guò)程中情景單一、數(shù)據(jù)量大、涉及環(huán)節(jié)多，存在計(jì)算復(fù)雜、操作繁瑣等問(wèn)題，將GIS與InVEST模型結(jié)合，研發(fā)農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別與治理模擬的集成計(jì)算機(jī)軟件系統(tǒng)，在多種情景下對(duì)關(guān)鍵源區(qū)進(jìn)行識(shí)別，并模擬治理后的效果，以期為農(nóng)田種植業(yè)面源污染治理提供有效的輔助手段。

1 ANSP關(guān)鍵源區(qū)識(shí)別軟件系統(tǒng)設(shè)計(jì)

基于GIS二次開發(fā)，集成水文網(wǎng)絡(luò)拓?fù)錁?gòu)建和面源污染關(guān)鍵源區(qū)識(shí)別技術(shù)，利用DEM數(shù)據(jù)提取河流節(jié)點(diǎn)、河段與子流域等水文地貌要素，整合InVEST產(chǎn)水量模型和營(yíng)養(yǎng)物傳輸率模型，結(jié)合基于入河污染負(fù)荷、潛在徑流污染物濃度和負(fù)荷與產(chǎn)水量比值3種情景的面源污染關(guān)鍵源區(qū)識(shí)別方法，以決策支持為重點(diǎn)，建立農(nóng)業(yè)面污染關(guān)鍵源區(qū)識(shí)別與治理模擬軟件系統(tǒng)。

1.1 系統(tǒng)目標(biāo)

通過(guò)研究ANSP的污染負(fù)荷情況與成因，運(yùn)用GIS技術(shù)，研制具有農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別與治理模擬的輔助決策軟件系統(tǒng)。對(duì)研究區(qū)域內(nèi)的水文網(wǎng)絡(luò)進(jìn)行拓?fù)錁?gòu)建、分析和專題信息可視化表達(dá)，快速高效得到研究區(qū)內(nèi)拓?fù)潢P(guān)系、子流域營(yíng)養(yǎng)元素分布規(guī)律和子流域潛在營(yíng)養(yǎng)元素徑流濃度等信息。在此基礎(chǔ)上，根據(jù)不同的方法分析面源污染關(guān)鍵源區(qū)空間位置，并模擬關(guān)鍵源區(qū)治理后的效果，為ANSP治理提供依據(jù)。

1.2 功能設(shè)計(jì)

以靈活的應(yīng)用集成接口構(gòu)建GIS軟件系統(tǒng)，對(duì)研究區(qū)內(nèi)的水網(wǎng)數(shù)據(jù)、DEM數(shù)據(jù)和土壤數(shù)據(jù)等多源異構(gòu)的資源環(huán)境數(shù)據(jù)進(jìn)行集成與調(diào)度，根據(jù)研究區(qū)已知的地理數(shù)據(jù)和資源環(huán)境數(shù)據(jù)，識(shí)別面源污染關(guān)鍵源區(qū)。系統(tǒng)按照實(shí)際功能分為6個(gè)模塊單元：1)拓?fù)潢P(guān)系分析單元，用于基于目標(biāo)區(qū)域的DEM數(shù)據(jù)構(gòu)建水文網(wǎng)絡(luò)拓?fù)潢P(guān)系；2)產(chǎn)水量分析單元，用于利用InVEST模型的產(chǎn)水量子模型，根據(jù)水文網(wǎng)絡(luò)拓?fù)潢P(guān)系，確定子流域產(chǎn)水量；3)入河負(fù)荷估算單元，根據(jù)地表營(yíng)養(yǎng)物負(fù)荷、傳輸率、次表層營(yíng)養(yǎng)物傳輸率和各柵格單元的入河負(fù)荷，計(jì)算子流域營(yíng)養(yǎng)物入河總負(fù)荷；4)潛在徑流濃度分析單元，用于根據(jù)子流域產(chǎn)水量，確定目標(biāo)區(qū)域內(nèi)各子流域的潛在氮、磷徑流濃度；5)關(guān)鍵源區(qū)識(shí)別單元，基于氮、磷入河負(fù)荷，潛在氮、磷徑流濃度，子流域氮、磷負(fù)荷與產(chǎn)水量比值3種情景，識(shí)別和確定目標(biāo)區(qū)域內(nèi)的面源污染關(guān)鍵源區(qū)；6)關(guān)鍵源區(qū)治理模擬單元，用于模擬對(duì)識(shí)別出的關(guān)鍵源區(qū)進(jìn)行治理后的結(jié)果，確定治理效果。

ANSP關(guān)鍵源區(qū)識(shí)別與治理模擬統(tǒng)一建模語(yǔ)言如圖1所示。由圖1可知，1種關(guān)鍵源區(qū)識(shí)別結(jié)果對(duì)應(yīng)1個(gè)關(guān)鍵源區(qū)識(shí)別模型，模擬治理結(jié)果依賴于關(guān)鍵源區(qū)識(shí)別結(jié)果；1個(gè)關(guān)鍵源區(qū)識(shí)別模型對(duì)應(yīng)若干農(nóng)田氮、磷入河負(fù)荷估算模型，子流域負(fù)荷與產(chǎn)水量比值估算模型與潛在氮、磷徑流濃度估算模型；農(nóng)田氮、磷入河負(fù)荷估算模型包含1個(gè)氮、磷流失負(fù)荷分配結(jié)果，1組地表營(yíng)養(yǎng)物傳輸率和1組次地表營(yíng)養(yǎng)物傳輸率；子流域負(fù)荷與產(chǎn)水量比值估算模型包含1個(gè)產(chǎn)水量模型，依賴于農(nóng)田氮、磷入河負(fù)荷估算模型結(jié)果；潛在氮、磷徑流濃度估算模型包含1個(gè)產(chǎn)水量模型，依賴于水網(wǎng)拓?fù)潢P(guān)系和農(nóng)田氮、磷入河負(fù)荷估算模型結(jié)果。ANSP關(guān)鍵源區(qū)識(shí)別與治理模擬軟件結(jié)構(gòu)如圖2所示。

圖1 ANSP關(guān)鍵源區(qū)識(shí)別與治理模擬統(tǒng)一建模語(yǔ)言Fig.1 Unified modelling language for the identification and control simulation of ANSP key source areas

圖2 ANSP關(guān)鍵源區(qū)識(shí)別與治理模擬軟件結(jié)構(gòu)Fig.2 Software structure diagram for identification and simulation of the treatment of ANSP key source areas

2 ANSP關(guān)鍵源區(qū)識(shí)別主要方法與模型

2.1 水文網(wǎng)絡(luò)拓?fù)錁?gòu)建

獲取研究區(qū)的DEM數(shù)據(jù)，基于DEM數(shù)據(jù)建立水文網(wǎng)絡(luò)拓?fù)潢P(guān)系。利用計(jì)算機(jī)編程提取河段以及子流域，徑流節(jié)點(diǎn)是基于幾何網(wǎng)絡(luò)的方法提取的，分為源頭節(jié)點(diǎn)、匯流節(jié)點(diǎn)、流域出口節(jié)點(diǎn)[9]。具體流程見圖3。

圖3 水文地貌要素提取Fig.3 Extraction of hydrological and geomorphic elements

流域水文網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的建立以流域匯流關(guān)系建立的基礎(chǔ)——河段為紐帶，對(duì)河網(wǎng)的徑流節(jié)點(diǎn)、河段、子流域進(jìn)行矢量要素的空間分析，構(gòu)建拓?fù)潢P(guān)系。河段屬性表里記錄對(duì)應(yīng)河段的起始節(jié)點(diǎn)和終止節(jié)點(diǎn)，以及對(duì)應(yīng)的子流域編碼；子流域?qū)傩员砝镉涗浥c其相對(duì)應(yīng)的河段編碼。定義節(jié)點(diǎn)的編碼為Point-ID，節(jié)點(diǎn)類型是NodeType，上下游節(jié)點(diǎn)分別為FROM_NODE、TO_NODE；河段編碼為 River_ID，長(zhǎng)度為River_len，匯入和匯出河段分別為Riv_in、Riv_out；子流域編碼為Basin_ID，每個(gè)子流域的匯入和匯出流域編碼分別為Basin_in、Basin_out。

2.2 流域氮、磷入河負(fù)荷計(jì)算方法

流域氮、磷入河負(fù)荷的計(jì)算分為2步，首先計(jì)算出流域氮、磷流失負(fù)荷，再將其與InVEST模型中的NDR模型進(jìn)行結(jié)合，計(jì)算得到流域氮、磷入河負(fù)荷。

2.2.1 氮、磷流失負(fù)荷模型

流域氮、磷流失負(fù)荷計(jì)算公式如下：

2.2.2 氮、磷入河負(fù)荷模型

通過(guò)模擬區(qū)域內(nèi)營(yíng)養(yǎng)物的來(lái)源和轉(zhuǎn)移過(guò)程，計(jì)算氮、磷入河系數(shù)，在考慮氮、磷流失負(fù)荷的情況下估算氮、磷入河負(fù)荷。首先按照式（2）將氮、磷流失負(fù)荷分配到表層和次表層：

最后，匯總計(jì)算氮、磷入河負(fù)荷，計(jì)算公式如下：

2.3 流域產(chǎn)水量計(jì)算方法

利用InVEST模型的產(chǎn)水量子模型計(jì)算流域產(chǎn)水量，計(jì)算公式如下：

為了反映流域多年平均產(chǎn)水量，降低年間的氣候差異和降水量差異造成不確定性的影響，潛在的蒸發(fā)散失量與降水量采用2001—2015年的平均值。

2.4 ANSP關(guān)鍵源區(qū)的識(shí)別

由于面源污染易受到地形、植被、氣候、施肥和土地利用類型等因素影響[12-15]，不同識(shí)別方法通常會(huì)呈現(xiàn)出不同程度的空間差異性。本研究采用3種不同情景識(shí)別關(guān)鍵源區(qū)。

2.4.1 基于入河負(fù)荷

該情景由分區(qū)統(tǒng)計(jì)獲得子流域氮、磷入河負(fù)荷，根據(jù)每個(gè)子流域氮、磷入河負(fù)荷大小進(jìn)行關(guān)鍵源區(qū)的識(shí)別[16-19]，計(jì)算公式如下：

2.4.2 基于潛在徑流濃度

該情景通過(guò)計(jì)算河段氮、磷入河通量與徑流量的比值，得到河段潛在徑流濃度，根據(jù)潛在徑流濃度的貢獻(xiàn)率大小進(jìn)行關(guān)鍵源區(qū)的識(shí)別[8]。

2.4.3 基于負(fù)荷與產(chǎn)水量比值

將子流域產(chǎn)水量與氮、磷入河負(fù)荷同時(shí)納入考慮，計(jì)算子流域產(chǎn)水量與氮、磷入河負(fù)荷的比值，根據(jù)該比值進(jìn)行關(guān)鍵源區(qū)的識(shí)別[20-22]，計(jì)算公式如下：

2.4.4 關(guān)鍵源區(qū)識(shí)別

根據(jù)式（8）～式（11）計(jì)算得到入河負(fù)荷、潛在徑流濃度和負(fù)荷與產(chǎn)水量比值，對(duì)子流域或河段對(duì)應(yīng)子流域進(jìn)行降序排列。按順序累加子流域入河負(fù)荷，每累加1次就通過(guò)子流域累計(jì)入河負(fù)荷與流域總負(fù)荷之比（Per）來(lái)判斷是否滿足設(shè)定閾值，若大于等于設(shè)定閾值則停止累加并將這些區(qū)域標(biāo)記為面源污染關(guān)鍵源區(qū)，流程如圖4所示。

圖4 關(guān)鍵源區(qū)識(shí)別流程Fig.4 Identification process for key source area

2.5 關(guān)鍵源區(qū)治理模擬

針對(duì)3種情境識(shí)別出的關(guān)鍵源區(qū)，設(shè)定相同的預(yù)期治理幅度，可直觀地觀察到在不同情景下對(duì)關(guān)鍵源區(qū)治理的效果。將3種情景下識(shí)別出的關(guān)鍵源區(qū)，按照與河流源頭的接近程度進(jìn)行排序，根據(jù)設(shè)定治理幅度依次降低這些區(qū)域的氮或磷負(fù)荷，再?gòu)纳嫌巍掠蔚捻樞蛞来胃铝饔蛩惺苡绊懙暮佣蔚牡蛄诐撛趶搅鳚舛取?/p>

3 農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別軟件應(yīng)用

3.1 研究區(qū)概況及特征

海河流域（112°E～120°E，30°N～43°N）地處我國(guó)華北地區(qū)，包括北京、天津、河南北部、山東北部、山西東部、內(nèi)蒙古和遼寧小部分區(qū)域與河北省大部分區(qū)域。流域地勢(shì)為西北高、東南低，包含灤河、北三河、大清河、永定河、子牙河、漳衛(wèi)南運(yùn)河、黑龍港運(yùn)東河、徒駭馬頰河與冀東沿海八大水系[23]。海河流域是我國(guó)三大糧食生產(chǎn)基地之一，流域土地利用類型主要為耕地，占流域總面積的48.9%。ANSP是海河流域主要污染源之一，也是海河水質(zhì)污染的主要原因之一[24]。

3.2 總氮關(guān)鍵源區(qū)識(shí)別

研發(fā)的農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別軟件主界面見圖5。以TN為例，使用該軟件對(duì)海河流域的農(nóng)業(yè)面源污染關(guān)鍵源區(qū)進(jìn)行分析。流域ANSP的面積與TN、TP的關(guān)系通常是呈對(duì)數(shù)上升趨勢(shì)，往往不到10%的面積就貢獻(xiàn)了30%以上的氮、磷污染負(fù)荷[25]，據(jù)此將軟件識(shí)別幅度設(shè)定為30%，分別基于入河污染負(fù)荷、潛在徑流污染物濃度、負(fù)荷與產(chǎn)水量比值3種情景進(jìn)行TN面污染關(guān)鍵源區(qū)識(shí)別，結(jié)果如圖6、圖7所示。

圖5 農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別軟件用戶界面Fig.5 User interface of the key source area identification software for ANSP

由圖6、圖7可知，TN入河負(fù)荷較高的區(qū)域主要分布在灤河、北三河、璋衛(wèi)河水系與黑龍港運(yùn)東河下游，子牙河、大清河水系中部以及徒駭馬頰河全域，呈現(xiàn)自西北向東南方向增長(zhǎng)的分布特征。潛在TN徑流濃度升高的區(qū)域主要分布中部—東南部地區(qū)的黑龍港運(yùn)東河、徒駭馬頰河與大清河水系。TN負(fù)荷與產(chǎn)水量比值高的區(qū)域集中分布在中部—南部地區(qū)的徒駭馬頰河上游、大清河水系及黑龍港運(yùn)東河。基于入河負(fù)荷情景下識(shí)別出的關(guān)鍵源區(qū)較為分散，除了永定河與子牙河水系外皆有分布?；跐撛趶搅鳚舛惹榫跋碌年P(guān)鍵源區(qū)分布較為集中，多為中部—南部分布，其余在東南部也有分布?；诹饔蜇?fù)荷產(chǎn)水比值情景下識(shí)別出的關(guān)鍵源區(qū)分布非常集中，以中部—南部分布為主，識(shí)別出的區(qū)域與基于潛在徑流濃度情景下的區(qū)域重合度較高。

圖6 3種情景下的TN污染區(qū)域分布情況Fig.6 Regional distribution of total nitrogen pollution under three scenarios

圖7 3種情景下識(shí)別的關(guān)鍵源區(qū)Fig.7 Key source areas identified by three scenarios

3.3 治理模擬結(jié)果與分析

在識(shí)別出關(guān)鍵源區(qū)的基礎(chǔ)上，對(duì)關(guān)鍵源區(qū)進(jìn)行模擬治理。為了更好地反映治理效果，將識(shí)別出的關(guān)鍵源區(qū)TN、TP入河負(fù)荷清零，并在考慮水網(wǎng)的拓?fù)潢P(guān)系下，對(duì)河段TN、TP潛在徑流濃度進(jìn)行重新統(tǒng)計(jì)。模擬結(jié)果如圖8所示。由圖8可知，基于入河負(fù)荷情景下，模擬治理的效果欠佳，流域中部水系的潛在徑流濃度下降不明顯；基于潛在徑流濃度、負(fù)荷與產(chǎn)水量比值2種情景下，治理效果相近，治理區(qū)域的潛在徑流濃度均有較為明顯的下降，但對(duì)東南部和南部的部分水系治理效果不明顯。

圖8 原始徑流濃度與3種情景下治理效果對(duì)比Fig.8 Comparison chart of original runoff concentration and that of three methods after treatment

計(jì)算了3種情景識(shí)別出的面污染關(guān)鍵源區(qū)模擬治理后TN、TP負(fù)荷的下降率，結(jié)果見表1。

由表1可知，在研究區(qū)的數(shù)據(jù)與條件下，基于入河負(fù)荷情景識(shí)別出的面污染關(guān)鍵源區(qū)的TN治理效果弱于基于潛在徑流濃度和基于負(fù)荷產(chǎn)水量比值情景的治理效果。由于海河流域西北地區(qū)入河負(fù)荷高，但林地和草地面積較大，攔截和稀釋了大量的氮、磷污染物，而中南部地區(qū)耕地面積占比高，且采取了一年兩熟的種植模式，化肥的施用量大且耗水嚴(yán)重，導(dǎo)致地表徑流量偏低，所以3種情景下識(shí)別出的關(guān)鍵源區(qū)有著較顯著的空間差異性。由于基于潛在徑流濃度情景反映的是河段與其上游輸入污染物的匯集和徑流過(guò)程共同作用的結(jié)果，而子流域負(fù)荷與產(chǎn)水量比值反映子流域產(chǎn)水量與氮、磷污染負(fù)荷比值，二者更能反映出ANSP對(duì)水質(zhì)的影響程度。

表1 3種情景下關(guān)鍵源區(qū)模擬治理效果Table 1 Simulation of governance effects in key source areas under three scenarios

對(duì)海河流域的治理模擬結(jié)果表明，僅針對(duì)面源污染入河負(fù)荷識(shí)別出的氮和磷污染關(guān)鍵源區(qū)在治理的結(jié)果上并不理想，基于潛在徑流濃度和基于負(fù)荷與產(chǎn)水量比值情景下的治理效果優(yōu)于前者，且結(jié)果接近。這是因?yàn)槊嬖次廴臼艿降匦?、植被和土地利用類型等因素影響較大，單一情景下識(shí)別效果并不一定能準(zhǔn)確反映面源污染對(duì)水質(zhì)的影響程度?？紤]到不同地區(qū)地理環(huán)境的復(fù)雜性，單一識(shí)別方法具有一定的局限性，在實(shí)際使用中，ANSP防控關(guān)鍵源區(qū)的識(shí)別應(yīng)將入河負(fù)荷、潛在徑流濃度和負(fù)荷產(chǎn)水量比值的空間分布情況結(jié)合分析，以期得到更好的治理效果，支持農(nóng)業(yè)面源污染治理科學(xué)決策。

4 結(jié)論

(1)針對(duì)農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別過(guò)程中情景單一、數(shù)據(jù)量大、涉及環(huán)節(jié)多，以及模型復(fù)雜、操作繁瑣和治理效果不明確等問(wèn)題，設(shè)計(jì)開發(fā)了多情景分析的農(nóng)業(yè)面源污染關(guān)鍵源區(qū)識(shí)別軟件。軟件具備拓?fù)潢P(guān)系分析、產(chǎn)水量分析、入河負(fù)荷估算、徑流濃度分析、關(guān)鍵源區(qū)識(shí)別和關(guān)鍵源區(qū)模擬治理等功能，可模擬3種常見情景下農(nóng)業(yè)面源污染分布的特點(diǎn)，為農(nóng)業(yè)面源污染治理決策提供有效的輔助手段。

(2)以海河流域?yàn)槔M(jìn)行了農(nóng)業(yè)面源污染關(guān)鍵源區(qū)的識(shí)別和治理模擬應(yīng)用，關(guān)鍵源區(qū)識(shí)別結(jié)果顯示，流域TN負(fù)荷呈現(xiàn)自西北向東南方向增長(zhǎng)的分布特征，中部—東南部地區(qū)的潛在TN徑流濃度較高，中部—南部地區(qū)的TN負(fù)荷與產(chǎn)水量比值明顯高于其他地區(qū)。模擬治理的結(jié)果表明，基于入河負(fù)荷情景對(duì)氮或磷的面源污染關(guān)鍵源區(qū)治理效果較不理想，而基于潛在徑流濃度、負(fù)荷與產(chǎn)水量比值情景的治理效果較好，且后者略優(yōu)于前者。