常曉敏，王少麗*，管孝艷，黃佳盛，賈海峰，尤李俊

（1.中國(guó)水利水電科學(xué)研究院，北京 100038；2.國(guó)家節(jié)水灌溉北京工程技術(shù)研究中心，北京 100048；3.青海省水利水電科學(xué)研究院有限公司，西寧 810001）

0 引 言*

【研究意義】水環(huán)境污染主要來(lái)自點(diǎn)源和面源2 種不同類型的污染物排放。隨著點(diǎn)源污染的有效控制，面源污染逐漸成為水環(huán)境惡化的主要原因[1-2]。農(nóng)業(yè)面源污染特征解析及空間識(shí)別是面源污染治理的基礎(chǔ)，通過(guò)污染物解析確定其主要來(lái)源及各來(lái)源的貢獻(xiàn)率，對(duì)于制定流域水環(huán)境污染防治策略具有重要意義?！狙芯窟M(jìn)展】研究表明，農(nóng)業(yè)面源污染已成為水體污染的重要污染源[3]，如何更好地控制面源污染將是我國(guó)水環(huán)境保護(hù)面臨的主要問(wèn)題之一[4-5]。據(jù)統(tǒng)計(jì)，由氮、磷負(fù)荷造成的水體污染占據(jù)我國(guó)水體總污染的81%和93%[6]?！兜诙稳珖?guó)污染源普查公報(bào)》顯示，2017 年農(nóng)業(yè)排放源中的COD、NH3-N、TN、TP 排放量分別占全國(guó)水污染物排放總量的49.8%、22.4%、46.5%、67.2%，面源污染形勢(shì)依舊嚴(yán)峻[7]。定量評(píng)估面源污染排放量是污染特征解析的基礎(chǔ)，目前應(yīng)用較多的評(píng)估方法有輸出系數(shù)法[8]、排污系數(shù)法[9]、監(jiān)測(cè)法、模型法[10]?！厩腥朦c(diǎn)】以往針對(duì)面源污染特征解析的研究大多圍繞單一種植區(qū)、養(yǎng)殖區(qū)或特定降雨季節(jié)開展，關(guān)于典型區(qū)域農(nóng)業(yè)面源污染的時(shí)空異質(zhì)性、污染物排放特征的長(zhǎng)時(shí)間序列解析、污染物排放局部聚集特征尚不多見。黃河流域是中國(guó)重要的生態(tài)屏障，流域生態(tài)環(huán)境保護(hù)和高質(zhì)量發(fā)展已被列為國(guó)家重大戰(zhàn)略。湟水河作為黃河上游的重要一級(jí)支流和主要水源之一，其水質(zhì)直接決定下游地區(qū)的生態(tài)安全?！?017 年青海省環(huán)境狀況公報(bào)》指出，湟水河水質(zhì)總體呈輕度污染狀態(tài)[11]。由于流域存在種植、養(yǎng)殖和農(nóng)村生活疊加污染的問(wèn)題，污染來(lái)源復(fù)雜且分布范圍較廣，污染負(fù)荷具有較高的時(shí)空變異性[12]。同時(shí)，種植業(yè)與畜牧養(yǎng)殖業(yè)的發(fā)展規(guī)模在空間上存在一定錯(cuò)位，部分地區(qū)畜禽糞便量超過(guò)了農(nóng)田承載力?，F(xiàn)有的湟水流域水環(huán)境研究多集中在斷面水質(zhì)狀況、水環(huán)境容量、污染物削減與措施等方面，而關(guān)于流域面源污染特征的定量解析及污染源區(qū)的識(shí)別相對(duì)較少?！緮M解決的關(guān)鍵問(wèn)題】鑒于此，本研究以湟水河小峽橋斷面控制區(qū)為研究區(qū)，基于2021 年6 月發(fā)布的《排放源統(tǒng)計(jì)調(diào)查產(chǎn)排污核算方法和系數(shù)手冊(cè)》（以下簡(jiǎn)稱“《手冊(cè)》”）中的核算方法[13]，定量估算研究區(qū)內(nèi)不同單元TN、TP、COD、NH3-N 的排放量和排放強(qiáng)度，分析其時(shí)空變化特征及污染物來(lái)源，探討各類污染物的排放強(qiáng)度，揭示研究區(qū)農(nóng)業(yè)面源污染的聚集規(guī)律，確定重點(diǎn)防控區(qū)域，為湟水流域水環(huán)境治理提供理論參考。

1 研究區(qū)概況與數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

湟水干流是黃河上游最大的一級(jí)支流，全長(zhǎng)374 km，流域總面積為32 863 km2。湟水干流在青海省境內(nèi)的干流長(zhǎng)為335.4 km，省內(nèi)流域分布總面積為16 224.36 km2，約占青海省總面積的2.3%，干流控制流域內(nèi)的人口數(shù)約為296 萬(wàn)，占全省總?cè)丝跀?shù)的60%，流域內(nèi)耕地面積為29.39 萬(wàn)hm2，占全省耕地面積的49%。小峽橋斷面是湟水河干流的出境斷面，屬于青海省人口集中且經(jīng)濟(jì)發(fā)達(dá)區(qū)域[14]。本研究在全國(guó)控制單元、青海省“十三五”水環(huán)境控制單元?jiǎng)澐?、流域自然匯水邊界以及行政單元邊界的基礎(chǔ)上，確定小峽橋斷面控制區(qū)的分布，如圖1 所示。

圖1 湟水河小峽橋斷面控制區(qū)及子單元?jiǎng)澐諪ig.1 Location of Xiaoxia bridge section of Huangshui River and division results of sub control units

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

以鄉(xiāng)鎮(zhèn)為統(tǒng)計(jì)單元，研究區(qū)內(nèi)涉及的各個(gè)區(qū)縣的鄉(xiāng)鎮(zhèn)作物播種面積、作物種植類型、作物產(chǎn)量、化肥和農(nóng)藥施用量、畜禽養(yǎng)殖數(shù)量、農(nóng)村人口等數(shù)據(jù)來(lái)源于各區(qū)縣2010—2018 年國(guó)民經(jīng)濟(jì)及社會(huì)發(fā)展統(tǒng)計(jì)資料；生活源及農(nóng)業(yè)源的產(chǎn)、排污系數(shù)參考《手冊(cè)》[13]。

2 研究方法

2.1 污染物負(fù)荷量估算

2.1.1 農(nóng)村生活源

根據(jù)《手冊(cè)》中的生活源污染物核算方法對(duì)農(nóng)村生活源的污染物負(fù)荷量進(jìn)行估算。青海省各地區(qū)農(nóng)村生活污水的污染物產(chǎn)污強(qiáng)度、綜合去除率參考生活污染源產(chǎn)、排污系數(shù)手冊(cè)中的農(nóng)村生活污水污染物產(chǎn)生與排放系數(shù)[13]。青海省對(duì)生活污水進(jìn)行處理的行政村比例來(lái)自《中國(guó)城鄉(xiāng)建設(shè)統(tǒng)計(jì)年鑒》，2010—2018 年依次為0.6%、1.0%、1.0%、1.2%、1.2%、1.3%、6%、6%、7.56%。

2.1.2 種植業(yè)源

1）化肥

種植業(yè)污染物（NH3-N、TN、TP）排放量采用《手冊(cè)》中的產(chǎn)、排污系數(shù)法核算[13]。研究區(qū)內(nèi)鄉(xiāng)鎮(zhèn)化肥施用量數(shù)據(jù)來(lái)源于各縣統(tǒng)計(jì)局，化肥施用量均為折純量。青海省種植業(yè)的氮、磷排放系數(shù)參考農(nóng)業(yè)污染源產(chǎn)、排污系數(shù)手冊(cè)中的種植業(yè)氮、磷排放系數(shù)。由于缺少園地面積數(shù)據(jù)，且該類面積的比例較少，因此只考慮農(nóng)作物的總播種面積，TN、TP、NH3-N 流失系數(shù)分別為0.225、0.016、0.016 kg/hm2。

2）作物秸稈

秸稈的面源污染物負(fù)荷排放量=秸稈污染物的產(chǎn)生量×秸稈面源污染物排放系數(shù)。其中，秸稈污染物的產(chǎn)生量=作物產(chǎn)量×秸稈與作物產(chǎn)量比值×秸稈產(chǎn)污系數(shù)。青海省主要農(nóng)作物包括小麥、玉米、豆類、薯類、油料作物。各類作物秸稈的TN、TP、COD產(chǎn)污系數(shù)參考該地區(qū)以往研究成果[15]，小麥、豆類、玉米、薯類和油料作物對(duì)應(yīng)的秸稈與作物產(chǎn)量比值分別為0.97、1.03、1.71、0.61 和3.0。

青海省秸稈資源利用主要以秸稈能源化、飼料化、秸稈還田、食用菌基料為主，未被利用的秸稈主要包括未收集與棄置亂堆的部分，其中未收集部分可假定為還田處理。本研究忽略能源化、飼料化、食用菌基料等利用方式的排污。作為肥料還田時(shí)的COD、TN、TP 排放系數(shù)分別為20%、10%和5%，棄置亂堆類秸稈的COD、TN、TP 排放系數(shù)均為50%[15]。秸稈利用現(xiàn)狀比例參考已有研究成果[15]，秸稈能源化、過(guò)腹還田、食用菌基料、秸稈還田、棄置亂堆量占秸稈總產(chǎn)量的比例分別為17%、38.8%、1.76%、27.64%、19.41%。

2.1.3 畜禽養(yǎng)殖業(yè)源

畜禽養(yǎng)殖業(yè)污染物的產(chǎn)生量和排放量均采用《手冊(cè)》中的產(chǎn)、排污核算法計(jì)算[13]。青海省規(guī)?；B(yǎng)殖場(chǎng)和養(yǎng)殖專業(yè)戶的排污系數(shù)參考農(nóng)業(yè)污染源產(chǎn)、排污系數(shù)手冊(cè)中的畜禽養(yǎng)殖排污系數(shù)。由于不同養(yǎng)殖處理糞便的方式、配套措施、飼養(yǎng)環(huán)境有所差別，因此排放系數(shù)的影響因素各異。本研究假設(shè)馬、驢、騾均為放養(yǎng)，其他畜禽養(yǎng)殖方式參考文獻(xiàn)[15]中的青海省畜禽養(yǎng)殖方式構(gòu)成比例?！妒謨?cè)》中只給出了規(guī)模養(yǎng)殖場(chǎng)和養(yǎng)殖戶的產(chǎn)、排污系數(shù)，未區(qū)分放養(yǎng)和散戶，本研究中散戶產(chǎn)、排污系數(shù)參考養(yǎng)殖專業(yè)戶的調(diào)查結(jié)果取值，放養(yǎng)產(chǎn)、排污系數(shù)參考規(guī)模養(yǎng)殖場(chǎng)和養(yǎng)殖戶的平均值。豬、肉牛、家禽的飼養(yǎng)量為當(dāng)年的出欄量，奶牛、羊、馬、驢、騾的飼養(yǎng)量為年末存欄量。馬、驢、騾的污染物產(chǎn)污系數(shù)及排放系數(shù)參考文獻(xiàn)[15]。

2.2 冷熱點(diǎn)分析

冷熱點(diǎn)分析是一種空間聚類方法，可展現(xiàn)指標(biāo)的高值和低值空間聚集分布規(guī)律，并彌補(bǔ)全局空間自相關(guān)和自然斷點(diǎn)分級(jí)對(duì)空間特征分析的不足，已應(yīng)用于生態(tài)系統(tǒng)服務(wù)價(jià)值空間異質(zhì)性及面源污染空間分布格局研究中。具體計(jì)算方法參考文獻(xiàn)[16]。

利用ArcGIS 冷熱點(diǎn)分析工具，在冷熱點(diǎn)分布圖屬性表中生成具有統(tǒng)計(jì)學(xué)意義的Z值和P值。按照Z(yǔ)值劃分出3 個(gè)熱點(diǎn)區(qū)、3 個(gè)冷點(diǎn)區(qū)、1 個(gè)不顯著區(qū)。劃分依據(jù)為一級(jí)冷點(diǎn)區(qū)（Z<-2.58），即99%置信區(qū)間；二級(jí)冷點(diǎn)區(qū)（-2.582.58），即 99%置信區(qū)間；二級(jí)熱點(diǎn)區(qū)（1.96

3 結(jié)果與分析

3.1 面源污染時(shí)間變化

2018 年，TN、TP、COD、NH3-N 的排放量分別為2 204.11、301.42、48 716.09、544.33 t/a，TN、COD、NH3-N 排放量比2010 年增加了5.3%、8.4%、10.4%，TP 排放量比2010 年減少了7.2%（表1）。由圖2 可知，2010—2018 年，來(lái)自種植業(yè)污染源的TN、TP、COD、NH3-N 排放量分別降低了18.39%、9.8%、22.7%、0.24%，主要原因是隨著綠色高效種植技術(shù)的推廣，區(qū)域內(nèi)種植業(yè)污染排放量顯著降低。來(lái)自畜禽養(yǎng)殖業(yè)污染源的TN、COD、NH3-N 排放量分別增加了9.65%、11.29%、15.91%，而TP 排放量降低了8.61%。其中，牛、羊養(yǎng)殖是畜禽養(yǎng)殖業(yè)污染中最主要的污染源，在TN、TP、COD、NH3-N的總排放量中，牛和羊養(yǎng)殖業(yè)的貢獻(xiàn)率分別為49.89%和15.50%、36.89%和23.95%、63.14%和9.46%、33.67%和17.30%。來(lái)自農(nóng)村生活污染源的TN、TP、NH3-N 排放量分別增加了3.42%、2.69%、3.08%，COD 排放量降低了0.35%。整體來(lái)看，流域面源污染與產(chǎn)業(yè)發(fā)展結(jié)構(gòu)密切相關(guān)，研究區(qū)TN、COD、NH3-N 排放量呈增加趨勢(shì)，TP 排放量減小，農(nóng)業(yè)面源污染排放量總體呈增長(zhǎng)趨勢(shì)，其中TN 排放量是TP 的7 倍，以往研究也指出TN 與TP 污染負(fù)荷量之比在5～10 之間[17]。

表1 2018 年研究區(qū)不同控制單元農(nóng)業(yè)面源污染排放量及貢獻(xiàn)率Table 1 Emissions and contribution rate of agricultural non-point source pollutants from different control units in the study area in 2018

圖2 研究區(qū)2010—2018 年農(nóng)業(yè)面源污染物TN、TP、COD、NH3-N 排放量變化Fig.2 Changes in emissions of agricultural non-point source pollutants TN, TP, COD and NH3-N in the study area from 2010 to 2018

3.2 面源污染空間分布

表1 為2018 年16 個(gè)子單元各類污染物排放量及貢獻(xiàn)率。從污染源結(jié)構(gòu)來(lái)看，排放貢獻(xiàn)率由高到低依次為：畜禽養(yǎng)殖業(yè)＞農(nóng)村生活＞種植業(yè)，畜禽養(yǎng)殖污染對(duì)TN、TP、COD、NH3-N 的排放貢獻(xiàn)率分別為73.54%、79.55%、78.20%、59.78%，農(nóng)村生活對(duì)TN、TP、COD、NH3-N 的排放貢獻(xiàn)率分別為17.65%、13.79%、21.54%、39.77%，種植業(yè)對(duì)TN、TP、COD、NH3-N 的排放貢獻(xiàn)率分別為8.81%、6.66%、0.27%、0.45%，研究區(qū)的首要污染源為畜禽養(yǎng)殖業(yè)。就面源污染指標(biāo)而言，COD 是研究區(qū)的首要污染物，占總排放量的比例為94.09%，TN、TP、NH3-N 占總排放量的比例分別為4.26%、0.60%、1.05%。這與以往研究[15]結(jié)論一致。鑒于此，研究區(qū)農(nóng)業(yè)面源污染治理應(yīng)將重點(diǎn)放在畜禽養(yǎng)殖業(yè)上，尤其應(yīng)重點(diǎn)關(guān)注COD 的減排，加強(qiáng)散戶畜禽廢污水中COD 的去除。

從空間分布來(lái)看，位于湟水干流下游的西鋼橋、橋頭橋控制區(qū)農(nóng)村人口多、農(nóng)業(yè)畜牧業(yè)發(fā)達(dá)，由此帶來(lái)的面源污染排放量較大。西鋼橋子單元污染排放總量最高，TN、TP、COD、NH3-N 排放量分別為558.61、77.91、11 725.96、119.51 t/a。因此，靠近西鋼橋區(qū)域的農(nóng)業(yè)面源污染不容忽視。西寧四區(qū)位于市區(qū)，經(jīng)濟(jì)發(fā)達(dá)，人口稠密，隨著城市化建設(shè)的加快，農(nóng)業(yè)種植面積及畜禽養(yǎng)殖量減少，農(nóng)村生活源成為該單元的主要污染源，對(duì)應(yīng)的TN、TP、COD、NH3-N 排放量占比分別為79.02%、73.59%、81.81%、88.33%。

在鄉(xiāng)鎮(zhèn)尺度上對(duì)2010 年和2018 年不同污染物排放強(qiáng)度進(jìn)行了對(duì)比，因4 類污染物排放強(qiáng)度存在差異，排放強(qiáng)度等級(jí)劃分也有所不同。由圖3 可知，TN、TP、COD、NH3-N 在2010 年的平均排放強(qiáng)度分別為282.85、44.68、6 037.39、64.84 kg/km2，2018 年的平均排放強(qiáng)度分別為289.74、41.26、6 418.91、66.35 kg/km2。整體來(lái)看，TN、COD、NH3-N 排放強(qiáng)度呈增加趨勢(shì)，TP 排放強(qiáng)度呈減小趨勢(shì)。

圖3 小峽橋斷面控制區(qū)TN、TP、COD 及NH3-N 排放強(qiáng)度分布Fig.3 Distribution of TN, TP, COD and NH3-N emission intensity of agricultural non-point source in Xiaoxiaqiao section control area

4 類污染物排放強(qiáng)度分布較為集中，在峽門橋和金灘斷面上游排放強(qiáng)度整體較小，進(jìn)入西寧城區(qū)的新寧橋斷面單元的排放強(qiáng)度整體較高。2010 年，TN 排放強(qiáng)度最高區(qū)域主要分布在大通縣良教鄉(xiāng)和石山鄉(xiāng)；2018 年新增了大通縣新莊鎮(zhèn)、塔爾鎮(zhèn)、多林鎮(zhèn)。2010 年，TP 排放強(qiáng)度最高區(qū)域主要分布于大通縣多林鎮(zhèn)、塔爾鎮(zhèn)、石山鄉(xiāng)以及湟中縣漢東鄉(xiāng)；2018 年新增了大通縣斜溝鄉(xiāng)和良教鄉(xiāng)。2010 年，COD 排放強(qiáng)度最高區(qū)域主要分布于大通縣良教鄉(xiāng)、石山鄉(xiāng)；2018 年新增了大通縣新莊鎮(zhèn)和西寧城中區(qū)。2010 年，NH3-N 排放強(qiáng)度最高的區(qū)域主要分布于大通縣城關(guān)鎮(zhèn)、塔爾鎮(zhèn)、良教鄉(xiāng)、石山鄉(xiāng)、長(zhǎng)寧鎮(zhèn)以及湟中縣漢東鄉(xiāng)和多巴鎮(zhèn)，2018 年新增了大通縣多林鎮(zhèn)。以上鄉(xiāng)鎮(zhèn)可作為研究區(qū)的污染重點(diǎn)治理區(qū)。

3.3 冷熱點(diǎn)格局分析

由于區(qū)域間污染物排放強(qiáng)度差異較大，排放強(qiáng)度可能呈一定的空間集聚特征，為了進(jìn)一步探究研究區(qū)面源污染集聚特征，利用冷熱點(diǎn)分析法對(duì)2010 年和2018 年的4 類污染物冷熱點(diǎn)空間格局進(jìn)行分析。由圖4 可知，不同污染物排放強(qiáng)度空間格局呈局部聚集性特征，TN、TP、COD 排放量呈較為明顯的空間溢出效應(yīng)，DT2、DT3、DT4 區(qū)域表現(xiàn)為空間連片的集聚性污染。TN、TP、COD 排放量熱點(diǎn)區(qū)主要分布在大通縣的橋頭橋（DT2、DT3）及新寧大通橋（DT4）區(qū)域，冷點(diǎn)區(qū)主要分布在海晏縣的金灘（HB）及互助土族自治縣的三其橋（HHZ）區(qū)域的部分鄉(xiāng)鎮(zhèn)。NH3-N 排放量熱點(diǎn)區(qū)主要分布在湟中縣的西鋼橋（HZ2）、新寧大通橋（DT4）、西寧城東區(qū)（XNN4）部分鄉(xiāng)鎮(zhèn)，冷點(diǎn)區(qū)分布范圍與TN、TP、COD 一致，但涉及鄉(xiāng)鎮(zhèn)面積較大。主要原因是大通縣牛羊養(yǎng)殖量較大，種植業(yè)較為發(fā)達(dá)；湟中縣的西鋼橋區(qū)域位于西寧人口聚集區(qū)，經(jīng)濟(jì)發(fā)達(dá)，污染排放量較多，互助縣種植業(yè)較為發(fā)達(dá)，近年來(lái)化肥零增長(zhǎng)政策驅(qū)動(dòng)、種植結(jié)構(gòu)調(diào)整、經(jīng)營(yíng)方式轉(zhuǎn)變使得種植業(yè)污染治理取得一定成效；海晏縣地廣人稀，單位面積的化肥施用量、畜禽養(yǎng)殖量和農(nóng)作物種植面積均較低，承載污染負(fù)荷能力相對(duì)較強(qiáng)，該區(qū)域污染程度差異性較小。其他地區(qū)污染物排放強(qiáng)度差異不大，未能形成明顯的聚集特征，為不顯著區(qū)。

圖4 小峽橋斷面控制區(qū)TN、TP、COD 及NH3-N 排放量的冷熱點(diǎn)空間分布Fig.4 Spatial distribution of TN, TP, COD and NH3-N hot and cold spots of agricultural non-point source pollution in Xiaoxiaqiao section control area

整體來(lái)看，區(qū)域污染物排放強(qiáng)度呈輕微增加趨勢(shì)，但局部地區(qū)污染聚集程度有所下降，熱點(diǎn)區(qū)面積減少，冷點(diǎn)區(qū)面積增加。盡管整體上污染物排放量和排放強(qiáng)度有所增加，但隨著生態(tài)文明建設(shè)的推進(jìn)，各地區(qū)污染物得到治理，局部地區(qū)污染程度的差異性有所降低，聚集程度減弱。

4 討 論

近年來(lái)，農(nóng)業(yè)面源污染防控已成為我國(guó)農(nóng)業(yè)工作的重點(diǎn)，由于受特定地域的自然條件、地形地貌、產(chǎn)業(yè)結(jié)構(gòu)等多因素的綜合影響，農(nóng)業(yè)面源污染產(chǎn)生和遷移具有不同的分布特征，污染關(guān)鍵源區(qū)也不同?；趨^(qū)域異質(zhì)性的農(nóng)業(yè)面源污染定量解析及空間格局研究更有利于區(qū)域污染治理策略的制定。

面源污染定量估算是污染特征解析的基礎(chǔ)，關(guān)于農(nóng)業(yè)面源污染估算方法較多，為規(guī)范排放源產(chǎn)、排量核算方法，統(tǒng)一產(chǎn)、排污系數(shù)，2021 年6 月生態(tài)環(huán)境部編制的《手冊(cè)》包含了農(nóng)業(yè)源產(chǎn)、排污核算方法和系數(shù)。本文采用該手冊(cè)核算法對(duì)研究區(qū)不同污染源排放量進(jìn)行了核算，核算結(jié)果與第二次全國(guó)污染源普查[18]、2015 年《中國(guó)環(huán)境年鑒》數(shù)據(jù)[19]、黃河流域等相關(guān)研究結(jié)果一致[20]?！肚嗪Ｊ〉诙稳珖?guó)污染普查》指出[21]，全省畜禽養(yǎng)殖業(yè)各類污染物排放量占農(nóng)業(yè)源排放量的比例高達(dá)90%，其中牛羊是畜禽養(yǎng)殖污染中最重要的污染源，這與張?zhí)冫惖萚22]研究結(jié)論一致。

畜禽養(yǎng)殖是湟水流域的首要污染源，因地制宜提出相應(yīng)的防控措施十分必要。當(dāng)前青海省畜禽養(yǎng)殖污染物處理的難點(diǎn)主要集中在區(qū)域性大量小規(guī)模的散養(yǎng)戶。已有數(shù)據(jù)表明，青海省畜禽養(yǎng)殖業(yè)廢水污染物削減率普遍較低，尤其是規(guī)模以下畜禽養(yǎng)殖戶的廢水污染物削減率僅為75.22%～86.69%，低于全省平均水平[21]。因此，應(yīng)加強(qiáng)規(guī)模以下畜禽養(yǎng)殖戶的糞污收集、處理、資源化利用水平，建議推廣“分散收集、集中處理”的模式。具體可采用固液分離處理技術(shù)降低畜禽排污水COD，其中絮凝分離技術(shù)方法簡(jiǎn)單、成本低廉，在青海省較為適用[23]。在未來(lái)發(fā)展中，農(nóng)牧融合、平衡種植養(yǎng)殖模式有望成為熱點(diǎn)，農(nóng)村散養(yǎng)的農(nóng)戶根據(jù)周圍土壤與種植業(yè)的實(shí)際承載量合理布局畜禽的養(yǎng)殖規(guī)模，配套糞污處理設(shè)施，讓畜禽的糞便和養(yǎng)殖污水經(jīng)過(guò)無(wú)害化處理后再次還田利用，真正實(shí)現(xiàn)種養(yǎng)良性生態(tài)循環(huán)[24]。

5 結(jié) 論

1）從時(shí)序變化來(lái)看，2010—2018 年研究區(qū)來(lái)自種植業(yè)污染源的排放量減小，來(lái)自畜禽養(yǎng)殖源和農(nóng)村生活源的排放量增加；TN、COD、NH3-N 排放量及排放強(qiáng)度總體呈增加趨勢(shì)，TP 減少。

2）從空間變化來(lái)看，4 類污染物排放強(qiáng)度分布比較集中，在峽門橋和金灘斷面上游單元整體較小，北川河支流從峽門橋斷面開始進(jìn)入到新寧橋斷面單元整體較高；從污染源結(jié)構(gòu)來(lái)看，排放貢獻(xiàn)率為畜禽養(yǎng)殖業(yè)源＞農(nóng)村生活源＞種植業(yè)源；就面源污染指標(biāo)而言，COD 是研究區(qū)的首要污染物。

3）從冷熱點(diǎn)空間格局來(lái)看，4 類污染物冷熱點(diǎn)空間格局呈局部集聚性分布特征，區(qū)域污染物排放強(qiáng)度整體呈輕微增加趨勢(shì)，但局部地區(qū)污染程度的聚集程度有所減弱。

（作者聲明本文無(wú)實(shí)際或潛在的利益沖突）