中圖分類號：X524；X821文獻(xiàn)標(biāo)志碼：A

Investigation and Risk Assessment of Non-point Source Pollution in a Small Hilly Watershed

A Case Study in Qiaodian Reservoir in Jinan City

SUN Diantai ¹ ，WANG Xin2， JIANG Xin2， WU Yuanchao ³ ， LIU Yuyu ¹ （1.School of Water Conservancyand Environment，University of Jinan，Jinan 25oO22，Shandong，China; 2.Water Resources Research Institute of Shandong Province， Jinan 25O014，Shandong，China; 3.Service Center for Jinan Municipal Water Conservancy Project，Jinan 25O099，Shandong，China）

Abstract：To grasp thecurrent situationof watershed non-point source polution inreservoir-type water sources in hily areas，and explore theriskofregional pollutant discharge and countermeasures forpreventionandcontrol，QiaodianReservoir n Jinan Citywas taken asthe studyarea，andtheoutputcoeficient method combined with theequivalent polution loadintensitywas used to estimateand analyze the agricultural non-point source pollution discharge，discharge intensity anddistributioncharacteristicsofthestudyarea in2O22，andclassfytherisk levelof pollution discharge.Theresults showthattheemissionintensityofagricultural non-point sourcepolutioninthe studyareaasa wholeisonthelow side， andthe highly poluted areas are mainly located in the upper reaches of the basin.The main sources of pollution in the study areaare livestock and poultry breeding sources，and the polution emissions from domestic sources and planting sourcesare relativelysmall，the dominant areas ofbreeding sources account forahigh percentage andare mostlydistributed in theupperreachesof the basin.The polution risk inthestudyareashows atrendof increasing from the downstream （the southwestern region） to the upper reaches （the northeastern region）.

Keywords：non-point source polltion；output factor method；equal standard polutionload intensity；smallhilly watershed;Qiaodian Reservoir

近年來，水庫逐漸成為城鎮(zhèn)生產(chǎn)、生活用水的重要來源，水庫型水源地日益成為經(jīng)濟(jì)社會可持續(xù)發(fā)展的重要保障[1]。與地下水水源地相比，地表水水源地具有開放性的特點，周圍自然環(huán)境和水體水質(zhì)易受人類活動的影響[2]。水庫型水源地流域中的鄉(xiāng)村養(yǎng)殖和農(nóng)業(yè)特色種植在促進(jìn)地方經(jīng)濟(jì)的發(fā)展的同時，也帶來了水源污染風(fēng)險[3-4]，此外，水庫周邊的農(nóng)村生活污水也會對水環(huán)境健康造成一定威脅[5]。近年來，我國對供水水庫水環(huán)境越來越重視，點源污染已逐步得到控制，面源污染成為眾多小流域水質(zhì)惡化的重要因素[6]

流域污染排放量的估算方法繁多，常用的方法主要有水文模型法、輸出系數(shù)法[7]、平均濃度法等[8]。在近年的研究中，水文模型的應(yīng)用較為廣泛，相關(guān)模型眾多，如土壤和水評估工具（SWAT）模型、水文模擬程序（HSPF）模型等[9]。由于大多數(shù)小流域存在基礎(chǔ)資料不足、面源污染監(jiān)測基礎(chǔ)資料匱乏[10]的問題，因此建立的水文模型精度不能保證。平均濃度法在實際應(yīng)用時容易受到環(huán)境變化和采樣誤差的影響。相比之下，輸出系數(shù)法避開了面源污染發(fā)生的復(fù)雜過程，所需參數(shù)少，操作簡便，適用性廣且結(jié)果較為準(zhǔn)確[]

位于省市鋼城區(qū)東南部的喬店水庫流域的污染來源較廣，種類多，污染現(xiàn)狀不明且基礎(chǔ)資料較少，因此，本文中選擇輸出系數(shù)法估算喬店水庫面源污染排放量，利用等標(biāo)污染負(fù)荷強度結(jié)合聚類分析法分析不同污染排放量和主要污染源?？紤]地形坡度影響，采用秩和比法對估算得到的喬店水庫流域的污染排放強度進(jìn)行等級劃分。開展供水水庫流域的污染負(fù)荷評估以及等級劃分，估算污染排放量，提出合理的評價結(jié)果，可以進(jìn)一步評估水源地流域內(nèi)的面源污染風(fēng)險，為污染的治理與防控提供依據(jù)和參考，提出有針對性的污染控制策略。

材料與方法

1.1 研究區(qū)概況

喬店水庫位于市鋼城區(qū)東南部的辛莊街道，地理位置為東經(jīng) 117^°50^′-118^°0^′ 、北緯 36^°9^′ 一36^°18^′ ，是市萊蕪、鋼城區(qū)的主要水源地，屬國家重要地表水水源地[12]，水庫總庫容為2799 萬 m³ ，總控制流域面積為 26km² 。全流域均為山丘區(qū)，地勢由南、北、東三面向西南方傾斜。多年平均降水量為750mm 左右，降水量的 70%～80% 集中在汛期，洪水暴漲暴落，水源較充足。本文中根據(jù)匯水原理及農(nóng)業(yè)管理現(xiàn)狀將整個流域分為9個區(qū)，如圖1所示（資料來源：國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站，經(jīng)ArcGIS10.2軟件數(shù)字化處理。）

1. 2 現(xiàn)場調(diào)研及數(shù)據(jù)搜集

在研究區(qū)開展流域調(diào)查，共調(diào)查23個村莊，通過問卷調(diào)查、口頭訪問、實地考察等方式，收集各個村落居住點的人口、排污方式、畜禽養(yǎng)殖類型和數(shù)量、種植類型和面積及耕作方式等相關(guān)資料。

數(shù)字高程模型（DEM）來源于中國科學(xué)院計算機網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺（http：//www.gscloud.cn），以此進(jìn)行子流域的劃分和坡度計算。根據(jù)高程數(shù)據(jù)結(jié)合匯水原理，將喬店水庫以上流域劃分為9個子流域，以子流域為基本單元進(jìn)行估算分析。

1.3 研究方法

1.3.1 輸出系數(shù)法

輸出系數(shù)法是根據(jù)土地利用類型與污染物輸出之間的關(guān)系而建立一種方法[13]，能夠充分利用相對容易得到的土地利用狀況等資料直接建立土地利用與受納水體非點源污染負(fù)荷的關(guān)系。

研究區(qū)流域內(nèi)不存在工業(yè)活動，廢氣排放較少，對雨水水質(zhì)影響不大，因此本文中排除降雨污染源和工業(yè)污染源的影響，僅計算流域農(nóng)業(yè)源與生活源的污染物排放量

Q_i=Y_i+Z_i+S_i，

式中： Q_i 為面源污染物 i 的年排放量； Y_i 為養(yǎng)殖源中污染物 i 的年排放量； Z_i 為種植源中污染物 i 的年排放量； S_i 為生活源中污染物 i 的年排放量。

式中： q_j 為第 j 類養(yǎng)殖的數(shù)量（出欄量）; f_ij 為第 j 類養(yǎng)殖的污染物 i 的排放系數(shù)； n 為養(yǎng)殖種類的數(shù)量。

式中： A_g 為農(nóng)作物總播種面積； e_ki 為第 k 類農(nóng)作種植的污染物 i 的流失系數(shù)； A_y 為園地作物總播種面積； e_yi 為第 y 類園地種植的污染物 i 的流失系數(shù)； σ_o 為種植種類的數(shù)量； q_i 為調(diào)查年度用于種植業(yè)的含氮（含磷）化肥使用面積折純用量； q₀ 為2017年度用于種植業(yè)的含氮（含磷）化肥使用面積折純用量。

S_i=rT×365，

式中： r 為常住人口； T 為產(chǎn)污濃度系數(shù)。

由于各子流域的面積不同，即使污染排放量相近，對環(huán)境的影響及污染程度也不同，因此將計算的污染排放量根據(jù)面積折算成排放強度，計算公式為

式中： I_i 為污染物 i 的排放強度； A 為某流域的面積。

選取化學(xué)需氧量（COD）、總磷（TP）、總氮（TN）3個指標(biāo)，在輸出系數(shù)中，農(nóng)業(yè)用地系數(shù)根據(jù)實地調(diào)研的作物管理方法，結(jié)合《排放源統(tǒng)計調(diào)查產(chǎn)排污核算方法和系數(shù)手冊》[14]（以下簡稱《手冊》）中給出的系數(shù)計算，氮肥、復(fù)合肥折純用量參考化肥種類相關(guān)資料以及2018年《省統(tǒng)計年鑒》[15]。畜禽養(yǎng)殖采用《手冊》以及相關(guān)文獻(xiàn)資料給出的產(chǎn)排污系數(shù)和折算標(biāo)準(zhǔn)[16]。農(nóng)村生活部分的面源污染物輸出系數(shù)采用國家環(huán)?？偩滞扑]的人口輸出系數(shù)及《手冊》中的省綜合去除率。開展水庫流域的基礎(chǔ)情況調(diào)查，估算2022年研究區(qū)各污染物排放水平[17]

1.3.2 等標(biāo)污染負(fù)荷強度

為了統(tǒng)一對比各指標(biāo)，將排放強度轉(zhuǎn)化為等標(biāo)污染負(fù)荷強度，即單位面積所排放污染物稀釋至環(huán)境控制標(biāo)準(zhǔn)所需的水量，以消除不同指標(biāo)的不同標(biāo)準(zhǔn)造成的影響[18] 。

式中： M_i 為污染物 i 的等標(biāo)負(fù)荷強度； ρ₀ 為污染物i 含量標(biāo)準(zhǔn)值（取《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中Ⅲ類水標(biāo)準(zhǔn)），以質(zhì)量濃度計； M_a 為污染源 α_a 的等標(biāo)負(fù)荷強度； M 為區(qū)域的等標(biāo)負(fù)荷強度； L 為污染物種類的數(shù)量； ?_m 為污染源種類的數(shù)量。

計算各區(qū)域的等標(biāo)污染負(fù)荷強，以每個子流域的等標(biāo)污染負(fù)荷強度為基本單元，對比省平均等標(biāo)污染負(fù)荷強度和流域平均強度開展分析討論[19] 。

1.3.3 聚類分析

聚類是數(shù)據(jù)挖掘、模式識別等研究方向的重要研究內(nèi)容之一，在識別數(shù)據(jù)的內(nèi)在結(jié)構(gòu)方面具有極其重要的作用[20]。本文中選擇 K 均值聚類算法，對各污染物的等標(biāo)污染負(fù)荷數(shù)據(jù)進(jìn)行聚類分析。采用系統(tǒng)聚類方法劃分不同區(qū)域的污染物主導(dǎo)因素，分析不同地區(qū)的污染類型[21]?；谖廴竞怂憬Y(jié)果，將ArcGIS軟件中的分組分析工具用于 K 均值聚類，系統(tǒng)聚類分析使用統(tǒng)計產(chǎn)品與服務(wù)解決方案軟件（SPSS）。

1.3.4 秩和比

秩和比（RSR）法具有易于理解和操作、克服主觀定權(quán)缺陷、綜合能力強等特點[22]，適用于有異常值或數(shù)值為0的情況，能夠反映所有評價指標(biāo)的信息及綜合水平。本文中選擇加權(quán)秩和比法劃分流域的污染排放等級，直觀體現(xiàn)流域各區(qū)域農(nóng)業(yè)面源污染排放強度。加權(quán)秩和比法依據(jù)各評價指標(biāo)的權(quán)重值，結(jié)合秩和比 R_sr 值得到加權(quán)秩和比（WRSR）值[23]，計算公式為

式中： R_sr，x 為第 x 個指標(biāo)的秩和比； R_b 為第 b 個被評價單位的秩； c 為被評價單位數(shù)； d 為指標(biāo)數(shù)。

當(dāng)各指標(biāo)權(quán)重不同時，計算WRSR值，

式中： W_rsr 為加權(quán)秩和比； w_j 為權(quán)重系數(shù)； R_cx 為第c 個被評價單位的第 x 個指標(biāo)的秩。

輸出系數(shù)法并沒有考慮地形因素對污染量排放的影響。相關(guān)研究[24]表明，降雨強度、下墊面因素大致相同時，坡度和土壤中氮、磷元素流失能力具有較大的正相關(guān)性。通過ArcGIS軟件計算每個子流域的平均坡度，根據(jù)污染物進(jìn)入水體的影響因素，在避免太高相關(guān)性的前提下，選擇TN和TP的等標(biāo)污染負(fù)荷強度、子流域平均坡度3個指標(biāo)，以特征值法主觀賦權(quán)，在計算 R_sr，x 值時根據(jù)不同權(quán)重加權(quán)，以 W_rsr 值分檔排序，進(jìn)而評估污染排放風(fēng)險等級。秩和比指標(biāo)的權(quán)重參數(shù)見表1。

2 結(jié)果與討論

2.1 污染排放量

污染物排放強度計算結(jié)果表明，2022年喬店水庫流域3類污染源共排放 cosΩ 、TN和TP污染物338.2t ，排放量分別為 317.4，18.1，2.7t ，分別占總排放量的 93.84%.5.36% 和 0.80% 。各污染源排放量由大到小的排序為養(yǎng)殖源、生活源、種植源，排放量（占比）分別為 258.4t（76.4%） 、 75.0t（22.2%） 人4.8t（1.4%） 。

各子流域的污染物排放量及流域的平均值如圖2所示。由圖可知：COD、TN和TP的排放量最大的區(qū)域均為4號子流域，分別占比 29.2% 、 25.4% 和 27.5% 。1、2、4、5號子流域COD排放量大于流域的平均值，排放量最小值出現(xiàn)于6號子流域。2、4、5、7號子流域TN排放量大于流域的平均值，排放量最小值位于9號子流域。2、4、5號子流域的TP排放量大于流域的平均值，排放量最小值位于8號流域。綜上分析，喬店水庫流域COD、TN以及TP的高排放區(qū)域相似度較高，出現(xiàn)于2、4、5號子流域，分布于流域3條支流的上游，可能與子流域面積較大、農(nóng)村居民生活較為密集和養(yǎng)殖規(guī)模較大有關(guān)。

2.2 等標(biāo)污染負(fù)荷

2.2.1等標(biāo)污染負(fù)荷排放強度

根據(jù)式（6）將計算得到的排放量轉(zhuǎn)化為等標(biāo)污染負(fù)荷強度，研究區(qū)各子流域等標(biāo)污染負(fù)荷強度如圖3所示。由圖可以看出，COD排放量較大，等標(biāo)污染負(fù)荷強度平均值約為省的平均值的 86% ，其中，1、2、4、5號子流域COD排放強度大于流域的平均值，最大值出現(xiàn)于4號子流域，最小值出現(xiàn)于6號子流域。 TN 、TP的等標(biāo)污染負(fù)荷強度較小，流域的平均值均遠(yuǎn)小于省的平均值，僅分別為省的平均值的 6.9% 和 5.5% 。 ΔTN 、TP的最大等標(biāo)污染負(fù)荷強度分別出現(xiàn)于4、2號子流域。各類污染物的高負(fù)荷區(qū)域均出現(xiàn)于流域上游畜禽養(yǎng)殖密集地區(qū)，說明養(yǎng)殖業(yè)對流域污染負(fù)荷影響較大；污染物最低負(fù)荷區(qū)域分別為6、8、9號子流域，均分布于流域下游及庫區(qū)附近，說明水庫管理范圍內(nèi)的環(huán)境保護(hù)措施較有成效。

喬店水庫流域各類污染物及污染源的等標(biāo)污染負(fù)荷占比如表2所示。污染物等標(biāo)污染負(fù)荷強度由大到小的順序為TN、COD、TP，在總量中占比分別為 39.94% ！ 31.93% 和 28.13% ，各污染物排放強度較為均衡。就污染源而言，等標(biāo)污染負(fù)荷強度由大到小的順序為養(yǎng)殖源、生活源、種植源，在總量的占比分別為 60.75% 、 23.40% 和 15.85% ，養(yǎng)殖源和生活源為主要污染源，貢獻(xiàn)了 73.97% 的COD和69.13% 的 TP 。種植源貢獻(xiàn)率較低，主要原因是研究區(qū)位于山丘區(qū)，耕地面積小，農(nóng)藥、化肥投入也較少，生活源和養(yǎng)殖源的污染占比相對較高。

COD、TN、TP等標(biāo)污染負(fù)荷強度最大區(qū)域的等標(biāo)污染負(fù)荷強度分別為最低區(qū)域的7.4、5.1、10.0倍，變異系數(shù)分別為0.65、0.46、0.66，說明不同地區(qū)的人類活動及自然差異對污染排放強度具有較大影響。分析流域污染負(fù)荷區(qū)域特征，識別區(qū)域主要污染貢獻(xiàn)源，可以因地制宜地提出相應(yīng)的管理措施，對流域面源污染防控有重大意義。

2.2.2 聚類分析

將研究區(qū)各污染物排放量劃分為4個污染等級，各子流域劃分為4個污染主導(dǎo)類型，結(jié)果分別如圖4、5所示。由圖4可知，COD、TN、TP等標(biāo)污染負(fù)荷強度最大值均在流域上游東北部2、4號子流域，污染負(fù)荷強度自下游到上游由小到大分布，這是由上、下游環(huán)境管理措施及產(chǎn)業(yè)結(jié)構(gòu)不同所致：上游地區(qū)較為發(fā)達(dá)的養(yǎng)殖業(yè)導(dǎo)致污染物排放量較大；中游地區(qū)為較為密集的居住區(qū)，因此污染負(fù)荷強度為中等。等標(biāo)污染負(fù)荷強度最小值均在庫區(qū)周邊，說明水庫直接匯流區(qū)域污染物排放較少。

由圖5可以看出，流域上游污染源多以養(yǎng)殖源為主，中游污染源主要為生活、種植源以及混合型，下游污染源則主要為養(yǎng)殖、生活源。從面積占比來看，養(yǎng)殖源主導(dǎo)型區(qū)域居多，且越靠近上游，養(yǎng)殖源污染占比越多，不同污染源主導(dǎo)型區(qū)域應(yīng)采取不同的污染治理措施，如研究區(qū)養(yǎng)殖業(yè)污染問題突出，治理重點應(yīng)放在畜禽養(yǎng)殖方面。

（資料來源：國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站，經(jīng)ArcGIS10.2軟件數(shù)字化處理。）

（資料來源：國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站，經(jīng)ArcGIS10.2軟件數(shù)字化處理。）

2.3 污染風(fēng)險等級劃分

根據(jù)各指標(biāo)的等標(biāo)污染負(fù)荷強度值計算各子流域的WRSR值，按數(shù)值大小排序、區(qū)間估計，共分為3個檔，分別對應(yīng)低污染（ （0～0.222） ）、中污染（ （gt;0.222～0.758） 和高污染（ gt;0.758～1.000） 。喬店水庫流域各區(qū)域WRSR值及污染等級分區(qū)分別如表3、圖6所示。由表3可知：WRSR值較大的1、2、4、5號子流域，環(huán)形分布于種植面積較多、養(yǎng)殖規(guī)模較大的上游地區(qū)，最大值出現(xiàn)于4號子流域；WRSR值中等的3、7號子流域，分布于人口較為密集、兼有少量養(yǎng)殖活動的中游區(qū)域；WRSR值較小的6、8、9號子流域分布于下游庫區(qū)附近。由圖6可知：流域污染程度自西南向東北遞增，2、4號子流域為高污染地區(qū)，該區(qū)域污染負(fù)荷較為嚴(yán)重，污染發(fā)生風(fēng)險最高，應(yīng)作為重點治理區(qū)域進(jìn)行防控。1、3、5、6、7、8號子流域為中污染地區(qū)，污染發(fā)生風(fēng)險偏高，須要針對區(qū)域排放特征開展防控。9號子流域為低污染地區(qū)，是養(yǎng)殖種植活動較少的庫區(qū)和自然保護(hù)區(qū)，該區(qū)域污染排放較少，污染風(fēng)險較低。污染等級由下游到上游呈遞增趨勢，說明喬店水庫流域的環(huán)境保護(hù)措施效果顯著。

2.4 污染防控建議

喬店水庫流域中、上游存在的畜禽養(yǎng)殖污染排放問題，主要原因為流域養(yǎng)殖業(yè)以小規(guī)模養(yǎng)殖戶為主，沒有形成規(guī)?；酿B(yǎng)殖場，管理較為粗放，大量畜禽糞便未經(jīng)處理就直接還田。對于養(yǎng)殖污染最為嚴(yán)重的流域東北部的2、4號子流域，應(yīng)考慮改進(jìn)管理方式，進(jìn)行畜禽糞便垃圾的收集工作，推廣糞便還田新工藝，調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，減小養(yǎng)殖規(guī)模，降低污染物的產(chǎn)生與排放。

生活源排放主要存在于3號子流域及庫區(qū)附近，主要問題是居民環(huán)保意識不足，仍存在污水和生活垃圾排入河道及在河道內(nèi)洗衣服的現(xiàn)象。對此應(yīng)開展水環(huán)境保護(hù)教育宣傳，建立獎懲機制，提高流域中已建成的污水收集設(shè)施的利用率，防止中游人口密集區(qū)域大量的生活源污染物直接進(jìn)入河道。也可于下游入庫區(qū)域修建矮壩、人工濕地等設(shè)施，阻擋與削減水體污染物

種植源排放問題在流域較為普遍，原因是流域土壤含沙量高，且山丘區(qū)坡度大，土壤營養(yǎng)物易流失。應(yīng)針對當(dāng)?shù)赝寥罈l件，選擇更為合適的農(nóng)藥、化肥種類，適時適量施用，同時提高秸稈、地膜的回收率，調(diào)整種植結(jié)構(gòu)，選擇需肥量更少且當(dāng)?shù)剡m宜種植的農(nóng)林作物等

綜合而言，相對于省污染物平均等標(biāo)污染負(fù)荷強度（COD、TN、TP的平均等標(biāo)污染負(fù)荷強度分別為 189917，2910657，2590938m³/km²） 和省內(nèi)相近大流域南四湖流域2013年污染物排放水平最小值（TN、TP的等標(biāo)污染負(fù)荷強度分別為2 156 000、 780 500m³/km² ）[25]，喬店水庫流域等標(biāo)污染負(fù)荷強度明顯較低（COD、TN、TP的等標(biāo)污染負(fù)荷強度分別為 ，說明在氣候、地理位置、種植結(jié)構(gòu)相似的情況下，由于單位面積農(nóng)業(yè)規(guī)模和管理方式上的差異，喬店水庫流域農(nóng)業(yè)面源污染物排放量較少，污染排放現(xiàn)象不嚴(yán)重。

從實地調(diào)研中發(fā)現(xiàn)，流域內(nèi)部分區(qū)域仍存在水體富營養(yǎng)化問題，說明即使污染排放量相對較少，開展各類污染源的防控治理工作依然十分重要。對于喬店水庫的污染防控研究，后續(xù)的工作應(yīng)重點放在以下2個方面：一是研究內(nèi)源污染；二是比較流域范圍內(nèi)下墊面污染流失能力。

3結(jié)論

本文中通過等標(biāo)污染負(fù)荷法對喬店水庫流域污染物排放水平進(jìn)行估算，得到以下結(jié)論：

1）喬店水庫流域2022年COD、 ΔTN 、TP排放量估算值分別為 317.4，18.1，2.7t，COD 排放量最大，占比約 93.84% ，是流域的主要污染排放物。TP排放量最小，占比為 0.80%

2）流域內(nèi)的等標(biāo)污染負(fù)荷強度由大到小的污染物順序為 ΔTN 、 cosΩ 、TP，主要污染物為TN和COD，占比分別為 39.94%.31.93% 。COD 排放強度最大值遠(yuǎn)大于省的平均值，負(fù)荷強度平均值略小于省的平均值。 TN 、TP等標(biāo)污染負(fù)荷強度的平均值和最大值均遠(yuǎn)小于省的平均值。各個子流域、各類污染物的等標(biāo)污染負(fù)荷強度存在明顯不同，不同區(qū)域相同污染物的污染負(fù)荷強度比值達(dá)到10.0，變異系數(shù)最大達(dá)到0.65，說明不同地區(qū)的人類活動及自然差異對污染排放強度具有較大影響。

3）研究區(qū)污染物排放量由大到小的污染源排序為養(yǎng)殖源、生活源、種植源，占比分別為 60.75% ！23.40% 和 15.85% 。污染物高排放區(qū)域分布于流域上游，且多為養(yǎng)殖源主導(dǎo)。流域下游多為生活源養(yǎng)殖源污染混合主導(dǎo)，下游及庫區(qū)附近的污染負(fù)荷強度較小，水庫管理區(qū)范圍內(nèi)的環(huán)境治理措施效果顯著。

4）喬店水庫流域上游的2、4號子流域污染風(fēng)險等級最高，9號子流域污染風(fēng)險最低，研究區(qū)污染風(fēng)險自西南到東北呈現(xiàn)遞增趨勢。

參考文獻(xiàn)：

［1］錢旭，蔡梅，王元元，等．浙江省銅山源水庫水體偏堿污染成因研究[J]．人民長江，2023，54（8）：67.

[2] 黃乾，薛雁，高印軍．山丘區(qū)水庫型水源地保護(hù)措施及建議[J]．水利，2019（3）：2.

[3] 史沖，王浩宇，王慧亮，等．基于SWAT的鲇魚山水庫流域氮磷面源污染時空變化研究[J]．華北水利水電大學(xué)學(xué)報（自然科學(xué)版），2023，44（2）：39.

[4] XUYX，MUT，YUAN Z，etal.Spatial patterns in pollution dis-charges from livestock and poultry farm and the linkage betweenmanure nutrients load and the carrying capacity of croplands in166006.

[5]LU MY， YANG S S，YU X L，et al. Decision support frameworkadapted to local conditions to select technologies for rural domesticsewage treatment in the Yangtze River Economic Belt[J]. Journalof Cleaner Production，2023，426：139067.

［6］可欣，于維坤，尹煒，等．小流域面源污染特征及其控制對策[J]．環(huán)境科學(xué)與技術(shù)，2009，32（7）：202.

［7］李建勛，李明萬，劉曉雨，等．基于輸出系數(shù)法的陜西農(nóng)業(yè)面源污染負(fù)荷研究[J]．生態(tài)與農(nóng)村環(huán)境學(xué)報，2023，39（4）：523.

［8］李懷恩．估算非點源污染負(fù)荷的平均濃度法及其應(yīng)用[J]．環(huán)境科學(xué)學(xué)報，2000（4）：398.

[9]ASADI S， MOUSAVI S， LOPEZ B A， et al. Analyzing hydrologi-calalterationandenvironmentalflows inahighlyanthropizedagri-cultural river basin system using SWAT+ ，WEAP and IAHRIS[J].Journal ofHydrology：Regional Studies，2024，52：101738.

［10］陶雙駿，邵光成，蘇江霖，等．小流域面源污染風(fēng)險評估研究：基于多分類有序離散選擇模型[J]．農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2017，36（7）：1293.

[11］薛利紅，楊林章．面源污染物輸出系數(shù)模型的研究進(jìn)展［J].生態(tài)學(xué)雜志，2009，28（4）：756.

［12］劉新剛，張兆迅，池京立．喬店水庫水源地保護(hù)存在的問題及對策[J]．中國水土保持，2019（8）：66.

[13］黃永剛，付玲玲，胡筱敏．基于河流斷面監(jiān)測資料的非點源負(fù)荷估算輸出系數(shù)法的研究和應(yīng)用［J]．水力發(fā)電學(xué)報，2012，31（5）：159.

[14］生態(tài)環(huán)境部．排放源統(tǒng)計調(diào)查產(chǎn)排污核算方法和系數(shù)手冊[EB/OL].（2021-06-21）[2024-03-20] https：//www.mee.gov.cn/xxgk2018/xxgk/xxgk01/202106/t20210618_839512. html.

［15］省統(tǒng)計局，國家統(tǒng)計局調(diào)查總隊．統(tǒng)計年鑒.北京：中國統(tǒng)計出版社，2018.

［16］李丹陽，元傳仁，衛(wèi)亞楠，等．中國北方地區(qū)羊養(yǎng)殖業(yè)產(chǎn)污系數(shù)測算[J]．農(nóng)業(yè)工程學(xué)報，2021，37（6）：225.

［17］彭舜磊，李鵬，王梓臣，等．河南省農(nóng)業(yè)面源污染負(fù)荷估算及區(qū)劃[J]．水土保持研究，2018，25（1）：226.

［18］丘麗清，吳根義，李想，等．基于等標(biāo)污染負(fù)荷強度的海南省農(nóng)村面源污染特征研究［J]．農(nóng)業(yè)資源與環(huán)境學(xué)報，2023，40（1） ： 57.

［19］于志曉．省農(nóng)業(yè)面源污染 TN、TP、COD 時空特征及防治區(qū)劃[D]．臨沂：臨沂大學(xué)，2023：18-19.

[20］孫吉貴，劉杰，趙連宇．聚類算法研究[J]．軟件學(xué)報，2008（1）： 48.

［21］蒙小俊，張建東，王秋利，等．漢江流域安康段面源污染特征分析[J]．湖南農(nóng)業(yè)科學(xué)，2023（10）：41.

[22］田鳳調(diào)．秩和比法在醫(yī)院統(tǒng)計中的應(yīng)用[J]．中國醫(yī)院統(tǒng)計，1994（1） ： 42.

［23］付學(xué)謙，陳皓勇．基于加權(quán)秩和比法的電能質(zhì)量綜合評估[J]．電力自動化設(shè)備，2015，35（1）：129.

［24］王麗，王力，王全九．不同坡度坡耕地土壤氮磷的流失與遷移過程[J]．水土保持學(xué)報，2015，29（2）：71.

［25］孫笑笑．南四湖流域非點源污染負(fù)荷估算及其對土地利用變化的響應(yīng)研究[D]．濟(jì)寧：曲阜師范大學(xué)，2018：27-29.

（責(zé)任編輯：于海琴）