渾河上游清原縣非點(diǎn)源污染負(fù)荷及其空間分布特征

劉志波 魏建兵 程全國(guó)

摘要 以渾河上游清原縣為研究對(duì)象，采用輸出系數(shù)法和ArcGIS平臺(tái)對(duì)清原縣各鄉(xiāng)鎮(zhèn)和渾河主要支流非點(diǎn)源污染總氮（TN）、總磷（TP）負(fù)荷進(jìn)行估算和空間分布特征分析，探討從行政分區(qū)和河流集水區(qū)2個(gè)尺度控制非點(diǎn)源污染的方法與指導(dǎo)治理實(shí)踐。結(jié)果表明，2017年清原縣TN、TP的輸出負(fù)荷量分別達(dá)3 410.66、208.47 t，重點(diǎn)污染鄉(xiāng)鎮(zhèn)為清原鎮(zhèn)、英額門(mén)鎮(zhèn)、南口前鎮(zhèn)、南山城鎮(zhèn);綜合污染負(fù)荷較重的支流為長(zhǎng)山堡河、古城子河和南雜木河，TN、TP 的輸出負(fù)荷量分別為397.41、331.58、319.03 t，29.63、22.39、20.04 t;清原流域接近于主干流的區(qū)域和支流污染較重，TN、TP污染負(fù)荷空間分布有相似性。旱田是TN的主要污染來(lái)源，TP主要污染源為旱田、居民生活和養(yǎng)殖;非點(diǎn)源污染控制應(yīng)關(guān)注年內(nèi)雨季和春季融凍期徑流的監(jiān)測(cè)和治理，生活污水和垃圾的收集與處理，分散型養(yǎng)殖和小型加工企的環(huán)保管理。

關(guān)鍵詞 非點(diǎn)源污染;輸出系數(shù)法;氮磷負(fù)荷;空間分布;渾河上游清原縣

中圖分類(lèi)號(hào) X 522? 文獻(xiàn)標(biāo)識(shí)碼 A? 文章編號(hào) 0517-6611（2022）06-0064-05

doi：10.3969/j.issn.0517-6611.2022.06.014

開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Non-point Source Pollution Load and Its Spatial Distribution Characteristics of Qingyuan County in the Upper Reaches of Hunhe River

LIU Zhi-bo， WEI Jian-bing， CHENG Quan-guo

（Key Laboratory of Ecological Restoration of Regional Polluted Environment， Ministry of Education， School of Environment， Shenyang University， Shenyang， Liaoning 110044）

Abstract Taking Qingyuan County in the upper reaches of the Hunhe River as the research object，the output coefficient method and ArcGIS platform were used to estimate the total nitrogen （TN） and total phosphorus （TP） loads of non-point source pollution in the towns of Qingyuan County and the main tributaries of the Hun River and analyze the spatial distribution characteristics.The purpose of this paper was to explore the methods and practice of non-point source pollution control from two scales of administrative division and river catchment. The results showed that in 2017，the annual output load of TN and TP in Qingyuan County was 3 410.66 and 208.47 t，respectively， and the key pollution towns were Qingyuan Town， Yingemen Town， Nankouqian Town and Nanshan Town.The tributaries with heavy comprehensive pollution load were Changshanbao River， Guchengzi River and Nanzamu River， with the annual output loads of TN and TP were 397.41 t， 331.58 t， 319.03 t， 29.63 t，22.39 t and 20.04 t， respectively.The area and tributaries close to the main stream in Qingyuan watershed were heavily polluted， and the spatial distribution of TN and TP pollution load was similar. The main pollution sources of TN were dry-land， then TP emission mainly came from dry-land， residents life and breeding industry.The control of non-point source pollution should pay attention to the monitoring and treatment of runoff in rainy season and spring thawing period， the collection and treatment of domestic sewage and garbage， and the environmental management of decentralized aquaculture and small processing enterprises.

Key words Non-point source pollution;Output coefficient method;TN and TP load;Spatial distribution;Qingyuan County in upper reaches of Hunhe River

基金項(xiàng)目 國(guó)家自然科學(xué)基金項(xiàng)目（41771200，41171399）。

作者簡(jiǎn)介 劉志波（1996 —），男，遼寧沈陽(yáng)人，碩士研究生，研究方向：環(huán)境規(guī)劃與可持續(xù)發(fā)展。*通信作者，教授，博士，碩士生導(dǎo)師，從事景觀生態(tài)學(xué)、環(huán)境生態(tài)學(xué)教學(xué)和科研工作。

收稿日期 2021-07-05

河流上游是遼東地區(qū)維護(hù)水源涵養(yǎng)、土壤保持、水質(zhì)凈化和生物多樣性等生態(tài)服務(wù)功能的重點(diǎn)關(guān)注區(qū)。隨著人口的增長(zhǎng)和社會(huì)經(jīng)濟(jì)的快速發(fā)展，水資源缺乏和水污染問(wèn)題不斷凸顯，河流源頭區(qū)的生態(tài)建設(shè)和環(huán)境保護(hù)逐漸得到重視。在我國(guó)，經(jīng)過(guò)長(zhǎng)期的流域規(guī)劃與治理，特別是“十一五”以來(lái)國(guó)家水體污染控制與治理科技重大專(zhuān)項(xiàng)、農(nóng)村環(huán)境連片整治項(xiàng)目的大量實(shí)施，作為水源保護(hù)目標(biāo)的大部分河流上游點(diǎn)源污染的排放得到了良好的控制。河流、湖泊、水庫(kù)等水體的主要污染負(fù)荷也由點(diǎn)源污染[1]逐漸轉(zhuǎn)變?yōu)檗r(nóng)業(yè)非點(diǎn)源污染[2]，尤其是氮、磷等非點(diǎn)源污染已經(jīng)成為迫切需要解決的水環(huán)境問(wèn)題[3]。為更好地解決農(nóng)業(yè)非點(diǎn)源污染的問(wèn)題，采用水文水質(zhì)模型進(jìn)行非點(diǎn)源污染的時(shí)空模擬和負(fù)荷量化研究成為重要的手段。然而現(xiàn)有的 SWAT、HSP、GBNP 等分布式模型對(duì)數(shù)據(jù)資料要求高，參數(shù)眾多且率定困難，對(duì)于廣大偏遠(yuǎn)而缺乏監(jiān)測(cè)資料的流域和地區(qū)，分布式機(jī)理模型參數(shù)率定和精度驗(yàn)證等成為模型應(yīng)用的難題[4] 。因此，選取和使用對(duì)數(shù)據(jù)資料要求較低、簡(jiǎn)單實(shí)用的輸出系數(shù)模型就成為國(guó)內(nèi)外研究此類(lèi)問(wèn)題的有效解決途徑之一[5]。Johnes[6]提出的輸出系數(shù)法在綜合考慮污染物侵蝕、遷移、轉(zhuǎn)化的基礎(chǔ)上，利用土地利用狀況等資料即可估算流域污染物輸出量，該方法對(duì)缺乏長(zhǎng)期水量、水質(zhì)數(shù)據(jù)流域的非點(diǎn)源污染負(fù)荷估算有顯著優(yōu)勢(shì)。

渾河上游流域是國(guó)家九大飲用水源地之一大伙房水庫(kù)的集水區(qū)，水庫(kù)供給沈陽(yáng)、撫順、大連等7大城市2 300萬(wàn)人口的工業(yè)和生活用水，在遼寧省經(jīng)濟(jì)社會(huì)發(fā)展中具有舉足輕重的戰(zhàn)略地位。2014年9月被國(guó)務(wù)院批準(zhǔn)確定為國(guó)家重點(diǎn)保護(hù)和優(yōu)先修復(fù)的湖泊水庫(kù)之一，根據(jù)環(huán)境部門(mén)的監(jiān)測(cè)，COD出現(xiàn)下降，但N、P污染超標(biāo)問(wèn)題的解決依然面臨挑戰(zhàn)。通過(guò)實(shí)地考察，流域已無(wú)較大的工業(yè)點(diǎn)源，農(nóng)業(yè)非點(diǎn)源污染成為渾河上游流域的主要污染源。近年來(lái)，很多科研院所和高校在該地區(qū)開(kāi)展了有關(guān)非點(diǎn)源污染的大量研究，基本摸清了2條主要入庫(kù)河流水質(zhì)污染程度及其與土地利用的關(guān)系，但限于缺乏河流流量等水文水質(zhì)模型需要的監(jiān)測(cè)數(shù)據(jù)，流域非點(diǎn)污染負(fù)荷計(jì)算的工作一直較為缺乏。

筆者選取渾河上游清原縣為研究區(qū)域，基于輸出系數(shù)法，從鄉(xiāng)鎮(zhèn)分區(qū)和主要支流流域尺度進(jìn)行非點(diǎn)源污染負(fù)荷評(píng)估，并分析其空間分布特征，為基于行政分區(qū)管理和集水范圍治理方案制定提供科學(xué)依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

清原縣位于遼寧省東部（圖1），地理坐標(biāo)為124°20′06″～125°28′58″E、41°47′52″～42°28′25″N，包括14個(gè)鄉(xiāng)鎮(zhèn)，面積3 932.96 km2。截至2017年，清原縣總?cè)丝跒?2.2萬(wàn)，其中城鎮(zhèn)人口10.5萬(wàn)，農(nóng)業(yè)人口21.7萬(wàn)。產(chǎn)業(yè)結(jié)構(gòu)以種植業(yè)和養(yǎng)殖業(yè)為主，耕地面積為69 427.39 hm2，主要農(nóng)作物類(lèi)型為水稻、玉米和大豆，農(nóng)田管理方式和化肥施用管理較為粗放;主要養(yǎng)殖類(lèi)型有大牲畜（1.82 萬(wàn)頭）、豬（0.50萬(wàn)頭）、羊（3.88萬(wàn)只）和家禽（35.23萬(wàn)只），主要養(yǎng)殖方式尚以散養(yǎng)為主。

1.2 數(shù)據(jù)來(lái)源及處理

研究區(qū)各鄉(xiāng)鎮(zhèn)人口數(shù)據(jù)（農(nóng)業(yè)人口與非農(nóng)業(yè)人口）來(lái)源于《2017年清原縣統(tǒng)計(jì)年鑒》;土地利用數(shù)據(jù)以清原縣統(tǒng)計(jì)資料為基礎(chǔ)，結(jié)合清原縣ALSO衛(wèi)星遙感影像解譯和現(xiàn)場(chǎng)調(diào)查獲得;畜禽養(yǎng)殖情況以清原縣畜牧局 2017 年畜禽養(yǎng)殖調(diào)查表為依據(jù)，并篩選出重點(diǎn)養(yǎng)殖鄉(xiāng)鎮(zhèn)進(jìn)行現(xiàn)場(chǎng)調(diào)查核實(shí);水系數(shù)據(jù)和流域邊界采用ArcGIS 10.2軟件水文分析功能模塊由清原縣DEM（數(shù)字地形模型）提取而來(lái)。

1.3 研究方法

1.3.1 輸出系數(shù)法。

輸出系數(shù)法將非點(diǎn)源污染源分為城鎮(zhèn)用地、農(nóng)村生活、農(nóng)田、人口、牲畜等幾大類(lèi)，根據(jù)土地利用類(lèi)型的差異、居民非點(diǎn)源污染物的排放和處理狀況、牲畜的數(shù)量和分布，對(duì)不同的土地利用和不同牲畜采用不同的輸出系數(shù)[7-8]，流域的總負(fù)荷量計(jì)算公式如下：

L=ni=1（Ei×Ai×Ii）+P（1）

式中，L 為營(yíng)養(yǎng)物（污染物）的輸出負(fù)荷（t/a）; Ei為第i 種營(yíng)養(yǎng)源（污染源）的輸出系數(shù)[t/（km2·a）或kg/（人·a）]; Ai為第i 類(lèi)土地利用類(lèi)型的面積（hm2）或第 i 種牲畜的數(shù)量、人口數(shù)量; Ii為第i 種營(yíng)養(yǎng)源的營(yíng)養(yǎng)物（污染物）的輸入量（t/a）; P 為降雨輸入的營(yíng)養(yǎng)物（污染物）的數(shù)量（t/a）。對(duì)于各種土地利用類(lèi)型來(lái)說(shuō)， Ei 表示土地利用中所流失的營(yíng)養(yǎng)物（污染物）的比例;對(duì)于牲畜而言，它表示牲畜排泄物直接進(jìn)入受納水體的比例，中間應(yīng)考慮人類(lèi)收集和儲(chǔ)存糞肥過(guò)程中的氨揮發(fā);對(duì)于人類(lèi)而言，它反映當(dāng)?shù)厝巳簩?duì)含磷去污劑的使用狀況、飲食和生活，對(duì)于農(nóng)村人口來(lái)說(shuō)，由于農(nóng)村沒(méi)有配套的污水處理設(shè)施且農(nóng)村人口分布分散，該系數(shù)主要反映農(nóng)村生活污水中的營(yíng)養(yǎng)物（污染物）進(jìn)入受納水體的比例[9]。

確定合理的輸出系數(shù)是輸出系數(shù)模型的關(guān)鍵，確定輸出系數(shù)有3種常用方法，即文獻(xiàn)分析法、試驗(yàn)?zāi)M法和水文統(tǒng)計(jì)法[10]。該研究采用文獻(xiàn)分析法[4，11-15]結(jié)合前期研究支流水質(zhì)采樣監(jiān)測(cè)數(shù)據(jù)綜合得到不同類(lèi)型污染源的TN、TP輸出系數(shù)（表1）。

1.3.2 支流集水區(qū)TN、TP輸出負(fù)荷計(jì)算及空間分布圖繪制。

使用河流集水區(qū)邊界裁切縣域TN、TP輸出負(fù)荷，跨多個(gè)鄉(xiāng)鎮(zhèn)的集水區(qū)采用面積加權(quán)求平均值的方法確定河流集水TN、TP輸出負(fù)荷，其計(jì)算公式如下：

Ei=EjAj+…+EnAnn（2）

式中，Ei為第i個(gè)河流集水區(qū)的負(fù)荷，Ej為第j個(gè)鄉(xiāng)鎮(zhèn)的負(fù)荷，Aj為第j個(gè)河流集水區(qū)在第j個(gè)鄉(xiāng)鎮(zhèn)內(nèi)的面積，n 為河流集水區(qū)跨越的鄉(xiāng)鎮(zhèn)數(shù)量。使用ArcGIS 10.2矢量數(shù)據(jù)屬性表輸入功能賦值支流集水區(qū)負(fù)荷并繪制空間分布圖。

2 結(jié)果與分析

2.1 行政分區(qū)角度TN、TP輸出負(fù)荷估算

根據(jù)輸出系數(shù)模型計(jì)算得到研究區(qū) 2017 年 TN和 TP 輸出負(fù)荷量分別為3 410.66、208.47 t，不同鄉(xiāng)鎮(zhèn)不同類(lèi)型污染源輸出量及其占比情況見(jiàn)表2。從表2可以看出，清原縣經(jīng)濟(jì)和人口中心清源鎮(zhèn)TN和TP的負(fù)荷均為最高，分別為399.17、29.83 t，經(jīng)濟(jì)較為落后與人口密度較小的敖家堡鄉(xiāng)TN和TP負(fù)荷最低，分別為149.31、7.70 t。

2.2 行政分區(qū)角度TN、TP污染負(fù)荷空間格局與構(gòu)成

清原縣產(chǎn)業(yè)以農(nóng)林業(yè)發(fā)展為主，城鎮(zhèn)區(qū)配套的污水管網(wǎng)措施不健全，農(nóng)村地區(qū)居民產(chǎn)生的生活污水和畜禽養(yǎng)殖的污水大多是直接排放或通過(guò)簡(jiǎn)易化糞池和沼氣池處理后排放;產(chǎn)生的沼液和沼渣被施入農(nóng)田，作為農(nóng)作物生產(chǎn)或苗木生產(chǎn)的基肥，與施入農(nóng)田和林下中草藥而未被植物吸收利用的化肥一起，在降雨徑流和融雪徑流的沖刷下逐漸累積和匯入水體，成為地表水污染的主要來(lái)源。

從圖2可看出，TN、TP 輸出負(fù)荷的空間分布格局基本一致，最高負(fù)荷量出現(xiàn)在中部的清原鎮(zhèn)，而位于清源鎮(zhèn)南北兩側(cè)的敖家堡、大蘇河鎮(zhèn)和枸乃甸污染負(fù)荷則相對(duì)較低。TN 輸出負(fù)荷較高的鄉(xiāng)鎮(zhèn)為南口前、英額門(mén)、清原鎮(zhèn)、南山城，共占清原縣鄉(xiāng)鎮(zhèn)總污染負(fù)荷的 38.01%;TP 輸出負(fù)荷較高的鄉(xiāng)鎮(zhèn)為清原鎮(zhèn)、南山城和夏家堡，共占全縣污染負(fù)荷的34.91%，主要是因?yàn)檫@些鄉(xiāng)鎮(zhèn)耕地面積廣、人口密集，受到農(nóng)業(yè)面源和農(nóng)村生活污染的影響較大。而敖家堡、枸乃甸人口密度小，主要以林地、草地為主。污染物輸出負(fù)荷量較高的鄉(xiāng)鎮(zhèn)對(duì)地表水環(huán)境威脅性大。在各污染源分類(lèi)中，TN主要來(lái)源于農(nóng)田污染，其污染負(fù)荷排放量占總體的71.82%，農(nóng)田中尤以旱田污染為重;TP主要來(lái)源于農(nóng)田污染，占比達(dá)53.19%，農(nóng)村生活和畜禽養(yǎng)殖污染貢獻(xiàn)率也較大，占比達(dá)到46.81%，畜禽養(yǎng)殖尤以豬和家禽的TP負(fù)荷貢獻(xiàn)率較大。因此減少TP污染負(fù)荷不僅要注意農(nóng)村生活污水收集與治理，也要關(guān)注旱田磷肥的合理施用和豬、家禽養(yǎng)殖污染排放的有效控制。

2.3 主要支流集水區(qū)角度TN、TP輸出負(fù)荷

除以行政分區(qū)為單元來(lái)加強(qiáng)環(huán)境的管理外，目前地表水環(huán)境的保護(hù)多采取“河長(zhǎng)制”的管理方式，以整條河流或某河段為單元，以當(dāng)?shù)氐男姓鞴茴I(lǐng)導(dǎo)和相關(guān)技術(shù)人員組成河流環(huán)境保護(hù)委員會(huì)，加強(qiáng)河流湖泊水環(huán)境、水生態(tài)和水景觀的巡視管理和制定治理規(guī)劃。因此，提取出研究區(qū)流域邊界和主要河流，估算河流集水區(qū)的污染負(fù)荷，對(duì)于明確重點(diǎn)治理河流和確定治理優(yōu)先順序具有指導(dǎo)價(jià)值。圖3為渾河上游清原流域主要支流水系及其集水區(qū)TN、TP負(fù)荷分級(jí)。

從TN負(fù)荷（圖3a）來(lái)看，各支流年負(fù)荷在150.00～397.41 t，高負(fù)荷支流為下玉城河、英額河、柳木橋河、長(zhǎng)春屯河、長(zhǎng)山堡河、古城子河和南雜木河，年負(fù)荷在319.03～397.41 t，較高負(fù)荷的支流為北腰堡河、十八嶺河、向陽(yáng)河和南口前河，年負(fù)荷在295.45～292.26 t，這些支流應(yīng)該優(yōu)先考慮治理。

TP負(fù)荷空間分布（圖3b）與TN趨勢(shì)大體相似，各支流年負(fù)荷在9.40～29.63 t，高負(fù)荷支流為長(zhǎng)山堡河、古城子河和南雜木河，年負(fù)荷在20.40～29.63 t，較高負(fù)荷的支流為北腰堡河、十八嶺河、向陽(yáng)河、南口前河、斗虎屯河、下玉城河、英額河、柳木橋河和長(zhǎng)春屯河，年負(fù)荷在15.00～17.69 t，這些支流應(yīng)該優(yōu)先考慮治理。從TN、TP綜合考慮的角度，長(zhǎng)山堡河、古城子河和南雜木河是非點(diǎn)源污染TN、TP控制較優(yōu)先的支流，TN、TP 的輸出負(fù)荷量分別為397.41、331.58、319.03，29.63、22.39、20.04 t。 結(jié)合野外調(diào)查分析，長(zhǎng)山堡河、古城子河集水區(qū)域農(nóng)業(yè)種植和養(yǎng)殖密度大，南雜木河集水區(qū)域居民密度大和小型加工業(yè)企業(yè)多，而這些區(qū)域污染物排放管理和集中處理設(shè)施薄弱。

3 討論

基于輸出系數(shù)模型和借助GIS技術(shù)對(duì)區(qū)域或流域尺度非點(diǎn)源污染TN、TP負(fù)荷進(jìn)行評(píng)估，并分析其空間分布特征與治理優(yōu)先順序。渾河上游清原縣主要以農(nóng)業(yè)生產(chǎn)為主，水環(huán)境重點(diǎn)的治理方向?yàn)榉屈c(diǎn)源污染。農(nóng)村生活污水、農(nóng)田地表徑流是研究區(qū)重點(diǎn)的治理方向，這個(gè)結(jié)論與張靜等[16]對(duì)渾河流域非點(diǎn)源污染負(fù)荷估算與分析結(jié)論類(lèi)似。研究區(qū)畜禽養(yǎng)殖以散養(yǎng)為主，其排放的畜禽糞便、有機(jī)物殘?bào)w等對(duì)地表水環(huán)境造成污染，這與趙倩等[17]對(duì)渾河上游大蘇河鄉(xiāng)面源污染負(fù)荷及現(xiàn)狀評(píng)價(jià)結(jié)果相符。在鄉(xiāng)鎮(zhèn)中心建立配套的污水管網(wǎng)作為試點(diǎn)是治理資金投入的重點(diǎn)，農(nóng)村地區(qū)污水管網(wǎng)未覆蓋的地方污水至少應(yīng)該簡(jiǎn)單處理后再排放，這與城鎮(zhèn)化程度和經(jīng)濟(jì)發(fā)展程度密切相關(guān)，同時(shí)也說(shuō)明該地區(qū)經(jīng)濟(jì)發(fā)展程度尚不足以支撐環(huán)境污染治理投入，需要關(guān)注經(jīng)濟(jì)發(fā)展與環(huán)境保護(hù)的協(xié)調(diào)問(wèn)題。盡早將農(nóng)村污水建設(shè)成網(wǎng)和進(jìn)一步加強(qiáng)農(nóng)田施肥管理仍是今后研究區(qū)水環(huán)境治理的重點(diǎn)和難點(diǎn)。

TN、TP負(fù)荷空間分布趨勢(shì)相似，這與李鑄衡等[18]對(duì)土地利用情景下渾河—太子河流域非點(diǎn)源污染模擬結(jié)論一致。根據(jù)各鄉(xiāng)鎮(zhèn)污染物TN、TP空間分布可看出水環(huán)境治理的重點(diǎn)應(yīng)放在清原鎮(zhèn)、英額門(mén)鎮(zhèn)、南口前鎮(zhèn);從主要支流TN、TP負(fù)荷空間分布格局來(lái)看，長(zhǎng)山堡河、古城子河和南雜木河是非點(diǎn)源污染控制較優(yōu)先的支流。

該研究確定污染源輸出系數(shù)采用了文獻(xiàn)分析法，為了使模型估算的準(zhǔn)確性更高，未來(lái)的研究應(yīng)考慮加入模擬試驗(yàn)和實(shí)地監(jiān)測(cè)的方法確定土地利用、農(nóng)村生活和畜禽養(yǎng)殖區(qū)域的輸出系數(shù)。對(duì)比胡正等[7，19-20]基于改進(jìn)的輸出系數(shù)模型對(duì)我國(guó)南方地區(qū)非點(diǎn)源污染估算評(píng)價(jià)，該研究采用的輸出系數(shù)模型尚未考慮降雨所帶來(lái)的營(yíng)養(yǎng)物質(zhì)對(duì)研究數(shù)據(jù)的影響，據(jù)該地區(qū)以往降雨測(cè)定表明，雨水中以氮沉降為主的營(yíng)養(yǎng)物質(zhì)濃度較低。因此，針對(duì)河流上游以森林覆蓋為主的流域，該研究選用的輸出系數(shù)模型對(duì)于非點(diǎn)源污染負(fù)荷的估算用于指導(dǎo)管理實(shí)踐是適用的，GIS技術(shù)也能夠?yàn)榉治鑫廴矩?fù)荷空間分布及確定優(yōu)先治理區(qū)域提供有力的支持。

4 結(jié)論

（1）渾河上游清原縣2017年非點(diǎn)源污染 TN、TP 的輸出負(fù)荷量分別為3 410.66、208.47 t，TN、TP 污染負(fù)荷的空間分布較為相似，重點(diǎn)污染鄉(xiāng)鎮(zhèn)為清原鎮(zhèn)、英額門(mén)鎮(zhèn)、南口前鎮(zhèn)、南山城鎮(zhèn)，其次為灣甸子鄉(xiāng)、土口子鄉(xiāng)、枸乃甸鄉(xiāng)。

（2）渾河上游清原縣2017年TN負(fù)荷主要來(lái)源為農(nóng)田面源，TP負(fù)荷中農(nóng)村生活和畜禽養(yǎng)殖污染貢獻(xiàn)率也較大，綜合污染負(fù)荷較重的支流為長(zhǎng)山堡河、古城子河和南雜木河，TN、TP 的輸出負(fù)荷量分別為397.41、331.58、319.03，29.63、22.39、20.04 t。

（3）渾河上游旱田是TN 的主要污染來(lái)源，TP主要污染源為旱田、居民生活和養(yǎng)殖，應(yīng)關(guān)注年內(nèi)雨季和春季融凍期徑流的監(jiān)測(cè)和治理，生活污水和垃圾的收集與處理，分散型養(yǎng)殖和小型加工企的環(huán)保管理。

參考文獻(xiàn)

[1] 中華人民共和國(guó)環(huán)境保護(hù)部，中華人民共和國(guó)國(guó)家統(tǒng)計(jì)局，中華人民共和國(guó)農(nóng)業(yè)部.第一次全國(guó)污染源普查公報(bào)[N].人民日?qǐng)?bào)，2010-02-10（016）.

[2] 第二次全國(guó)污染源普查公報(bào)[J].環(huán)境保護(hù)，2020，48（18）：8-10.

[3] 國(guó)務(wù)院印發(fā)《“十三五”生態(tài)環(huán)境保護(hù)規(guī)劃》? 提出12項(xiàng)約束性指標(biāo)? 到2020年生態(tài)環(huán)境質(zhì)量總體改善[J].環(huán)境經(jīng)濟(jì)，2016（ZA）：8-9.

[4] 劉增進(jìn)，張關(guān)超，楊育紅，等.河南省農(nóng)業(yè)非點(diǎn)源污染負(fù)荷估算及空間分布研究[J].灌溉排水學(xué)報(bào)，2016，35（11）：1-6.

[5] WORRALL F，BURT T P.The impact of land-use change on water quality at the catchment scale：The use of export coefficient and structural models[J].Journal of hydrology，1999，221（1/2）：75-90.

[6] JOHNES P J.Evaluation and management of the impact of land use change on the nitrogen and phosphorus load delivered to surface waters：The export coefficient modelling approach[J].Journal of hydrology，1996，183（3/4）：323-349.

[7] 胡正，敖天其，李孟芮，等.改進(jìn)的輸出系數(shù)模型在缺資料地區(qū)面源綜合評(píng)價(jià)[J].灌溉排水學(xué)報(bào)，2019，38（2）：108-114.

[8] 張辰，陸建忠，陳曉玲.基于輸出系數(shù)模型的云南洱海流域農(nóng)業(yè)非點(diǎn)源污染研究[J].華中師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，51（1）：108-114.

[9] 劉亞瓊，楊玉林，李法虎.基于輸出系數(shù)模型的北京地區(qū)農(nóng)業(yè)面源污染負(fù)荷估算[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27（7）：7-12.

[10] 張洪波，李俊，黎小東，等.缺資料地區(qū)農(nóng)村面源污染評(píng)估方法研究[J].四川大學(xué)學(xué)報(bào)（工程科學(xué)版），2013，45（6）：58-66.

[11] 胡富昶，敖天其，胡正，等.改進(jìn)的輸出系數(shù)模型在射洪縣的非點(diǎn)源污染應(yīng)用研究[J].中國(guó)農(nóng)村水利水電，2019（6）：78-82，92.

[12] 楊雯，敖天其，王文章，等.基于輸出系數(shù)模型的瓊江流域（安居段）農(nóng)村非點(diǎn)源污染負(fù)荷評(píng)估[J].環(huán)境工程，2018，36（10）：140-144.

[13] 程靜，賈天下，歐陽(yáng)威.基于STELLA和輸出系數(shù)法的流域非點(diǎn)源負(fù)荷預(yù)測(cè)及污染控制措施[J].水資源保護(hù)，2017，33（3）：74-81.

[14] 蔡明，李懷恩，莊詠濤，等.改進(jìn)的輸出系數(shù)法在流域非點(diǎn)源污染負(fù)荷估算中的應(yīng)用[J].水利學(xué)報(bào)，2004，35（7）：40-45.

[15] 鄭堯文.渾河清原流域地表水氮素污染特征及其源匯研究[D].沈陽(yáng)：沈陽(yáng)大學(xué)，2019.

[16] 張靜，何俊仕，周飛，等.渾河流域非點(diǎn)源污染負(fù)荷估算與分析[J].南水北調(diào)與水利科技，2011，9（6）：69-73.

[17] 趙倩，馬建，問(wèn)青春，等.渾河上游大蘇河鄉(xiāng)農(nóng)業(yè)非點(diǎn)源污染負(fù)荷及現(xiàn)狀評(píng)價(jià)[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2010，26（2）：126-131.

[18] 李鑄衡，劉淼，李春林，等.土地利用變化情景下渾河-太子河流域的非點(diǎn)源污染模擬[J].應(yīng)用生態(tài)學(xué)報(bào)，2016，27（9）：2891-2898.

[19] 陳亞榮，阮秋明，韓鳳翔，等.基于改進(jìn)輸出系數(shù)法的長(zhǎng)江流域面源污染負(fù)荷估算[J].測(cè)繪地理信息，2017，42（1）：96-99，104.

[20] 段揚(yáng)，蔣洪強(qiáng)，吳文俊，等.基于改進(jìn)輸出系數(shù)模型的非點(diǎn)源污染負(fù)荷估算：以嫩江流域?yàn)槔齕J].環(huán)境保護(hù)科學(xué)，2020，46（4）：48-55.