郝 旭，張乃明，史 靜

（云南農(nóng)業(yè)大學(xué) 資源與環(huán)境學(xué)院，云南 昆明650201）

非點(diǎn)源污染的危害性隨點(diǎn)源污染控制能力的提高而逐漸表現(xiàn)出來(lái)，特別是當(dāng)點(diǎn)源污染控制水平達(dá)到一定程度后，非點(diǎn)源污染勢(shì)必成為水環(huán)境污染的主要來(lái)源［1］。從世界范圍來(lái)看，非點(diǎn)源污染已成為水環(huán)境污染的重要方式。美國(guó)60%的水環(huán)境污染起源于非點(diǎn)源［2］，據(jù)計(jì)算奧地利北部地區(qū)進(jìn)入水環(huán)境的非點(diǎn)源氮量比點(diǎn)源多，丹麥270條河流94%的氮負(fù)荷，52%的磷負(fù)荷是由非點(diǎn)源污染引起的［3］，荷蘭農(nóng)業(yè)非點(diǎn)源提供的總氮、總磷分別占水環(huán)境污染總量的60%和40%～50%［4］。在我國(guó)，長(zhǎng)江流域、云南洱海、北京密云水庫(kù)等許多水域，非點(diǎn)源污染比例均超過(guò)點(diǎn)源污染。在此背景下，進(jìn)行非點(diǎn)源污染的量化研究對(duì)水質(zhì)管理及飲用水源地保護(hù)顯得愈發(fā)重要。

應(yīng)用數(shù)學(xué)模型對(duì)非點(diǎn)源污染進(jìn)行模擬和估算是研究非點(diǎn)源污染的重要手段。已有的非點(diǎn)源模型大都是針對(duì)小區(qū)域的精細(xì)模擬，率定出的參數(shù)和結(jié)果帶有較大的區(qū)域性，很難適用于大中尺度流域的負(fù)荷估算。而輸出系數(shù)模型法所需參數(shù)少，操作簡(jiǎn)便易行，又具有一定的精度，且不考慮非點(diǎn)源污染發(fā)生和發(fā)展的復(fù)雜過(guò)程，因此，該模型在非點(diǎn)源污染負(fù)荷估算方面得到了一定的應(yīng)用［5－12］。

近幾年，西南地區(qū)水資源緊缺，該區(qū)連續(xù)3a大旱，加之人類活動(dòng)對(duì)昆明市云龍水庫(kù)匯水范圍內(nèi)的生態(tài)環(huán)境的干預(yù)逐步加重，導(dǎo)致水庫(kù)匯水區(qū)內(nèi)生態(tài)和水環(huán)境狀況日益惡化。但就目前情況而言，對(duì)云龍水庫(kù)庫(kù)區(qū)的基本狀況、基礎(chǔ)數(shù)據(jù)的調(diào)查與收集還遠(yuǎn)遠(yuǎn)不夠，對(duì)水庫(kù)水質(zhì)變化原因，污染來(lái)源與貢獻(xiàn)，水環(huán)境容量及庫(kù)區(qū)農(nóng)業(yè)面源污染控制技術(shù)等都缺乏系統(tǒng)研究。因此，本研究利用輸出系數(shù)模型，針對(duì)云龍水庫(kù)水源保護(hù)區(qū)的氮、磷非點(diǎn)源污染負(fù)荷進(jìn)行分析，為進(jìn)一步掌握研究區(qū)的氮磷負(fù)荷特征和合理調(diào)整農(nóng)業(yè)結(jié)構(gòu)提供科學(xué)依據(jù)。

1 研究方法

1.1 研究區(qū)概況

云龍水庫(kù)是昆明市最重要也是云南省最大的集中式飲用水源地。水庫(kù)位于金沙江水系二級(jí)支流掌鳩河中上游祿勸縣云龍鄉(xiāng)，距昆明市約86km，位于東經(jīng)102°16′—102°35′，北緯25°43′—26°07′，屬于北亞熱帶季風(fēng)氣候，多年平均氣溫10.2～13.8℃。云龍水庫(kù)控制徑流面積為745km2，流域呈矩形狀，地勢(shì)北高南低，流域平均海拔高程約為2 406m。水庫(kù)總庫(kù)容4.84×108m3，正常蓄水量3.97×108m3，主要入庫(kù)河流有石板河、老木河、水城河。云龍水庫(kù)水源區(qū)總?cè)丝诩s53 601人，經(jīng)濟(jì)以農(nóng)業(yè)為主。水源區(qū)實(shí)際耕地面積約2.27×104hm2，林地面積4.82×104hm2。全年牛、豬、羊、家禽等存欄數(shù)約21萬(wàn)頭。水源區(qū)經(jīng)濟(jì)以農(nóng)業(yè)為主，產(chǎn)業(yè)結(jié)構(gòu)單一。

1.2 非點(diǎn)源輸出模型及參數(shù)取值

輸出系數(shù)模型是利用半分布式途徑來(lái)計(jì)算某流域年均污染（總氮、總磷）總負(fù)荷的數(shù)學(xué)加權(quán)公式，其實(shí)質(zhì)上是一種集總模型。根據(jù)非點(diǎn)源污染負(fù)荷產(chǎn)生量和土地利用類型之間的聯(lián)系，利用流域相關(guān)數(shù)據(jù)直接建立流與非點(diǎn)源污染產(chǎn)出量的關(guān)系［13］，然后通過(guò)對(duì)不同來(lái)源污染負(fù)荷量求和，得到研究區(qū)域污染物總負(fù)荷。

Johnes等的輸出系數(shù)模型方程為:

式中:L——營(yíng)養(yǎng)物質(zhì)的總流失量；Ei——第i種營(yíng)養(yǎng)源的輸出系數(shù)；Ai——Ii第i類土地利用類型的面積或畜禽數(shù)量或人口數(shù)量；Ii——第i類營(yíng)養(yǎng)源的營(yíng)養(yǎng)物質(zhì)輸入量。

1.2.1 不同土地利用類型非點(diǎn)源負(fù)荷 不同土地利用類型產(chǎn)生的氮、磷非點(diǎn)源負(fù)荷計(jì)算式為:

式中:j——土地利用類型，共n種；Eij——第i種污染物在第j種土地利用類型單位面積的輸出系數(shù)〔kg/（km2·a）〕；Mij——第i種污染物在第j種土地利用類型中的營(yíng)養(yǎng)物輸入量（kg）；Aj——第j種土地利用類型的面積（km2）。

1.2.2 農(nóng)村生活非點(diǎn)源負(fù)荷 農(nóng)村生活污染源中，氮、磷非點(diǎn)源負(fù)荷的計(jì)算方法為:

式中:L（生活）i——第i種農(nóng)村生活污染物年輸出量（kg/a）；Di——第i種污染物的日輸出系數(shù)〔kg/（人·d）〕；H——區(qū)域內(nèi)人口數(shù)量（人）；M——污水處理過(guò)程中的機(jī)械去除營(yíng)養(yǎng)物質(zhì)系數(shù)；B——污水處理過(guò)程中的生物去除營(yíng)養(yǎng)物質(zhì)系數(shù)。

1.2.3 畜牧養(yǎng)殖非點(diǎn)源負(fù)荷 畜禽污染源中氮、磷非點(diǎn)源負(fù)荷計(jì)算方法為:

式 中:L（畜牧）i——第i種污染物年輸出量（kg/a）；Dij——第j種牲畜第i種污染物日輸出系數(shù)〔kg/（頭·d）〕；Hj——區(qū) 域 內(nèi) 第j種 牲 畜 的 數(shù) 量（頭）；M——污水處理過(guò)程中的機(jī)械去除營(yíng)養(yǎng)物質(zhì)系數(shù)；B——污水處理過(guò)程中的生物去除營(yíng)養(yǎng)物質(zhì)系數(shù)。1.2.4 輸出系數(shù)取值的確定 大量的試驗(yàn)研究［14］表明，影響輸出系數(shù)的因素有地形、氣候、水文、土地利用類型、植被覆蓋率等，都會(huì)對(duì)輸出系數(shù)產(chǎn)生較大的影響。輸出系數(shù)多采用單位時(shí)間、單位面積上的負(fù)荷量表示。一般確定非點(diǎn)源污染物輸出系數(shù)值可采用試驗(yàn)或調(diào)研的方法［15］。鑒于無(wú)云南地區(qū)的輸出系數(shù)試驗(yàn)數(shù)據(jù)，本研究計(jì)算所需的輸出系數(shù)值，參考國(guó)外不同土地利用方式下的氮和磷輸出系數(shù)的變化范圍的研究成果和參照國(guó)內(nèi)［1，11，16－23］相似 自然、氣候 條件下的其他地區(qū)研究結(jié)果并取其平均值作為云龍水庫(kù)徑流區(qū)的輸出系數(shù)值。不同土地利用類型的參考輸出系數(shù)值和本文的具體取值結(jié)果（表1）；畜禽糞尿和農(nóng)村生活部分的非點(diǎn)源污染物輸出系數(shù)值均參照國(guó)家環(huán)?？偩滞扑]的排泄系數(shù)和人口輸出系數(shù)［24］，并參考前人研究成果得到農(nóng)村生活及畜禽糞尿輸出系數(shù)（表2）。畜禽糞尿的氮磷輸出系數(shù)分別取為各自排泄系數(shù)的12%［24］。而農(nóng)村生活的輸出系數(shù)則反映了當(dāng)?shù)剞r(nóng)村人口對(duì)生活污水及垃圾的利用和處理水平。

估算非源污染負(fù)荷產(chǎn)生量所需的其它數(shù)據(jù)包括研究時(shí)段的土地利用或植被覆蓋數(shù)據(jù)、人口數(shù)量以及大牲畜、豬、羊和家禽數(shù)量，可通過(guò)查閱云南省昆明市祿勸統(tǒng)計(jì)年鑒、云南省數(shù)字鄉(xiāng)村網(wǎng)［25］以及文獻(xiàn)資料得到或?qū)嵉卣{(diào)查獲得。

表1 不同土地利用類型的輸出系數(shù) （kg/（hm2·a）

表2 農(nóng)村生活及畜禽糞便輸出系數(shù) kg/（只·a）

2 結(jié)果與分析

2.1 水源區(qū)各鄉(xiāng)鎮(zhèn)非點(diǎn)源氮、磷負(fù)荷產(chǎn)生量分析

根據(jù)云龍水庫(kù)自然和社會(huì)條件，應(yīng)用適用于大區(qū)域尺度非點(diǎn)源估算的輸出系數(shù)模型方法，估算云龍水庫(kù)徑流區(qū)氮磷非點(diǎn)源污染產(chǎn)生量、入庫(kù)量等。昆明市云龍水庫(kù)非點(diǎn)源污染物總氮和總磷負(fù)荷產(chǎn)生量計(jì)算結(jié)果詳見(jiàn)表3—4。

表3 2010年昆明市云龍水庫(kù)非點(diǎn)源污染物總氮負(fù)荷產(chǎn)生量

表4 2010年昆明市云龍水庫(kù)非點(diǎn)源污染物總磷負(fù)荷產(chǎn)生量

由表3—4可知，云龍水庫(kù)徑流區(qū)總氮和總磷的非點(diǎn)源負(fù)荷產(chǎn)生量分別為2 079.8和454.32t/a，且總氮負(fù)荷量是總磷負(fù)荷量的4.6倍。從空間分布上分析，撒營(yíng)盤(pán)鎮(zhèn)、馬鹿塘鄉(xiāng)、皎平渡鎮(zhèn)、團(tuán)街鎮(zhèn)、插甸鄉(xiāng)5個(gè)鄉(xiāng)鎮(zhèn)的非點(diǎn)源污染物總氮、總磷負(fù)荷量都較高。其中，撒營(yíng)盤(pán)鎮(zhèn)的最高，總氮、總磷年均負(fù)荷量達(dá)到了444.89和85.28t/a。因?yàn)檫@些鄉(xiāng)鎮(zhèn)的耕地面積、農(nóng)村人口、禽畜數(shù)量較多，均是產(chǎn)業(yè)結(jié)構(gòu)單一，經(jīng)濟(jì)以農(nóng)業(yè)為主，第二、三產(chǎn)業(yè)發(fā)展相對(duì)滯后的鄉(xiāng)鎮(zhèn)。相比之下，云龍鄉(xiāng)和茂山鎮(zhèn)的農(nóng)村人口和禽畜數(shù)量較少，因此非點(diǎn)源污染負(fù)荷產(chǎn)出總量較低。按其對(duì)水源區(qū)污染物總氮負(fù)荷產(chǎn)生總量貢獻(xiàn)率排序?yàn)?撒營(yíng)盤(pán)鎮(zhèn)＞插甸鄉(xiāng)＞皎平渡鎮(zhèn)＞團(tuán)街鄉(xiāng)＞馬鹿塘鄉(xiāng)＞云龍鄉(xiāng)＞茂山鎮(zhèn)。各鄉(xiāng)鎮(zhèn)對(duì)于總磷負(fù)荷量的有所變動(dòng)，其排列順序?yàn)?撒營(yíng)盤(pán)鎮(zhèn)＞插甸鄉(xiāng)＞團(tuán)街鄉(xiāng)＞皎平渡鎮(zhèn)＞馬鹿塘鄉(xiāng)＞云龍鄉(xiāng)＞茂山鎮(zhèn)。綜合分析可知，撒營(yíng)盤(pán)鎮(zhèn)的氮、磷污染產(chǎn)生量是水源庫(kù)區(qū)非點(diǎn)源污染物總氮和總磷負(fù)荷量的主要影響因素。

2.2 不同污染源類型對(duì)總氮和總磷負(fù)荷產(chǎn)生量的貢獻(xiàn)分析

由表3可以發(fā)現(xiàn)，農(nóng)村生活污染源對(duì)非點(diǎn)源污染物總氮負(fù)荷量貢獻(xiàn)率最高33.84%。不同污染源類型產(chǎn)生的總氮負(fù)荷量以及對(duì)總氮負(fù)荷量的貢獻(xiàn)率的大小順序?yàn)?農(nóng)村生活區(qū)＞畜禽糞便（其中，豬＞牛＞羊＞家禽）＞耕地＞林地。不同土地利用方式對(duì)流域非點(diǎn)源污染物總氮負(fù)荷產(chǎn)出量的影響，以農(nóng)村生活區(qū)污染最大，畜禽養(yǎng)殖業(yè)和耕地的次之，林地的最小。表4的數(shù)據(jù)顯示，污染源類型對(duì)總磷負(fù)荷貢獻(xiàn)率的影響與對(duì)總氮污染負(fù)荷的影響規(guī)律基本相似，但畜禽糞便污染對(duì)總磷負(fù)荷的貢獻(xiàn)明顯大于農(nóng)村生活區(qū)污染和耕地的貢獻(xiàn)。不同污染源類型對(duì)總磷負(fù)荷貢獻(xiàn)率的大小順序?yàn)?畜禽糞便（其中，豬＞家禽＞牛＞羊）＞農(nóng)村生活區(qū)＞耕地＞林地。這些表明，在各種土地利用類型中，目前水源的畜禽養(yǎng)殖和農(nóng)村居民生活用地對(duì)非點(diǎn)源污染物總氮和總磷的貢獻(xiàn)最大，因此為減少環(huán)境污染，農(nóng)村地區(qū)污水、生活垃圾及禽畜糞便的日常處理顯得非常重要。除此之外，大量畜禽養(yǎng)殖是造成該區(qū)非點(diǎn)源污染物總氮和總磷負(fù)荷量較高的直接原因。但不同鄉(xiāng)鎮(zhèn)畜禽的年存欄量和年出欄數(shù)量是不固定的，而且數(shù)量變動(dòng)隨機(jī)性大。因此，由于畜禽養(yǎng)殖業(yè)的發(fā)展引起環(huán)境質(zhì)量的不利影響需要引起足夠的重視，并且對(duì)各鄉(xiāng)鎮(zhèn)畜禽種類、年存欄量、年出欄數(shù)量定時(shí)統(tǒng)計(jì)上報(bào)，為研究工作提供基礎(chǔ)數(shù)據(jù)。

2.3 水源區(qū)入庫(kù)氮磷污染負(fù)荷

入庫(kù)污染物是影響水庫(kù)水質(zhì)的根本原因。根據(jù)《全國(guó)水環(huán)境容量核算指南》中提供的技術(shù)參數(shù)，結(jié)合云龍水庫(kù)的實(shí)際情況，分別確定各污染源的入庫(kù)系數(shù):耕地和林地的氮磷流失根據(jù)相關(guān)文獻(xiàn)分析可取入庫(kù)系數(shù)為0.1；而農(nóng)村生活和畜禽糞便產(chǎn)生的污染負(fù)荷大部分通過(guò)河道匯集后進(jìn)入水庫(kù)，取這類污染源的氮、磷入庫(kù)系數(shù)分別0.1和0.05。從而計(jì)算得到云龍水庫(kù)非點(diǎn)源總氮、總磷入庫(kù)量見(jiàn)圖1—2。

圖1 昆明市云龍水庫(kù)入庫(kù)總氮負(fù)荷

從圖1—2可以看出，由不同污染源產(chǎn)生的總氮和總磷的入庫(kù)負(fù)荷量分別為207.98和24.73t/a，總氮入庫(kù)總量是總磷入庫(kù)總量的8.4倍。入庫(kù)的總氮負(fù)荷的主要污染來(lái)源是農(nóng)村生活區(qū)污染和畜禽糞便污染，其入庫(kù)總氮負(fù)荷量分別為70.38和70.90t/a，占總?cè)霂?kù)總氮負(fù)荷量的33.84%～34.09%。水庫(kù)周邊村落是水庫(kù)水環(huán)境的重要污染源之一，周邊村落污水目前均未治理。村落污水非汛期常以氮磷污染物高濃度狀態(tài)的聚集于村落洼地及溝渠之中，汛期暴雨將其一起攜帶入庫(kù)，直接污染庫(kù)區(qū)水體。

圖2 昆明市云龍水庫(kù)入庫(kù)總磷負(fù)荷

根據(jù)不同污染源入庫(kù)的總氮負(fù)荷貢獻(xiàn)率的順序?yàn)?農(nóng)村生活＞畜禽糞便（其中大牲畜＞家禽）＞耕地＞林地。從圖2可看出，入庫(kù)總磷負(fù)荷的主要污染源為畜禽糞便污染和農(nóng)村生活污染，其入庫(kù)總磷負(fù)荷為12.50 和 8.21t/a，分 別 占 總 負(fù) 荷 的 50.54% 和33.21%，而林地產(chǎn)生的入庫(kù)總磷負(fù)荷為1.05t/a，占總?cè)霂?kù)量的百分比也最低，為4.23%。根據(jù)不同污染源的入庫(kù)總磷負(fù)荷貢獻(xiàn)率的順序?yàn)?畜禽糞便（其中，大牲畜＞家禽）＞農(nóng)村生活＞耕地＞林地。

3 結(jié)論

（1）云龍水庫(kù)徑流區(qū)氮磷負(fù)荷總的產(chǎn)生量分別為2 079.8和454.32t/a，按其對(duì)水源區(qū)污染物負(fù)荷產(chǎn)生總量貢獻(xiàn)率總氮排序?yàn)?撒營(yíng)盤(pán)鎮(zhèn)＞插甸鄉(xiāng)＞皎平渡鎮(zhèn)＞團(tuán)街鄉(xiāng)＞馬鹿塘鄉(xiāng)＞云龍鄉(xiāng)＞茂山鎮(zhèn)；總磷負(fù)荷貢獻(xiàn)排列順序?yàn)?撒營(yíng)盤(pán)鎮(zhèn)＞插甸鄉(xiāng)＞團(tuán)街鄉(xiāng)＞皎平渡鎮(zhèn)＞馬鹿塘鄉(xiāng)＞云龍鄉(xiāng)＞茂山鎮(zhèn)。

（2）不同污染源類型產(chǎn)生的總氮負(fù)荷及對(duì)總氮負(fù)荷的貢獻(xiàn)率大小順序?yàn)?農(nóng)村生活區(qū)＞畜禽糞便（其中，豬＞牛＞羊＞家禽）＞耕地＞林地。其中農(nóng)村生活污染對(duì)非點(diǎn)源污染物總氮負(fù)荷貢獻(xiàn)率最高33.84%。林地貢獻(xiàn)率最小為11.74%。不同污染源類型對(duì)總磷負(fù)荷貢獻(xiàn)率的大小順序?yàn)?畜禽糞便（其中，豬＞家禽＞牛＞羊）＞農(nóng)村生活區(qū)＞耕地＞林地。其中畜禽糞便污染對(duì)總磷負(fù)荷的貢獻(xiàn)明顯大于農(nóng)村生活區(qū)和耕地的貢獻(xiàn)。

（3）水庫(kù)徑流區(qū)不同污染來(lái)源產(chǎn)生的總氮和總磷的入庫(kù)負(fù)荷量分別為207.98和24.73t/a，總氮入庫(kù)總量是總磷入庫(kù)總量的8.4倍。入庫(kù)總氮負(fù)荷的主要污染來(lái)源是農(nóng)村生活區(qū)污染和畜禽糞便污染，其入庫(kù)總氮負(fù)荷分別為70.38和70.90t/a，占總?cè)霂?kù)總氮負(fù)的33.84%～34.09%。入庫(kù)總氮負(fù)荷貢獻(xiàn)率的順序?yàn)?農(nóng)村生活＞畜禽糞便（其中，大牲畜＞家禽）＞耕地 ＞林地。按入庫(kù)總磷負(fù)荷貢獻(xiàn)率的大小順序?yàn)?畜禽糞便（其中，豬＞家禽＞牛＞羊）＞農(nóng)村生活區(qū)＞耕地＞林地。

［1］邱涼，羅小勇，程紅光.長(zhǎng)江流域大尺度空間非點(diǎn)源污染負(fù)荷研究［J］.人民長(zhǎng)江，2011，42（18）:81.

［2］賀纏生，傅伯杰.非點(diǎn)源污染的管理及控制［J］.環(huán)境科學(xué)，1998，19（5）:15－18.

［3］Kronvang B.Diffusenutrient losses in Denmark［J］.Water Sci.Techno.，1996，133（1）:81－88.

［4］Boersp C M.Nutriente missions from agriculture in the Nether lands:Causes and remedies［J］.Water Sci.Techno.，1996，33（1）:183－191.

［5］ClesceriN L，Curran S J，Sedlak R I.Nutrient loads in Wiscons in lakes（Part 1）:Nitrogen and phosphorous export coefficients［J］.Water Resources Bulletin，1986，22（6）:983－989.

［6］Dickerh D，Haith D A.Loading functions for predicting nutrient losses from complex watersheds［J］.Water Resources Bulletin，1983，19（6）:951－959.

［7］Reckhow K H，Beaulac M N，Simpson J T.Modeling Phosphorus Loading and Lake Response Under Uncertainty:A Manual and Compilation of Export Coefficient［S］. U.S.EPA Report No.EPA－440/5－80－011，Washington D C:Office of Water Regulations，Criteria and Standards Division，U.S.Environmental Protection Agency，1980.

［8］Ierodiaconou D，Laurenson L，Leblanc M，et al.The consequences of land use change on nutrient exports:A regional scale assessment in South－west Victoria，Australia［J］.J.Environ Manage.，2005，74（4）:305－316.

［9］劉瑞民，楊志峰，丁曉雯，等.土地利用/覆蓋變化對(duì)長(zhǎng)江上游非點(diǎn)源污染影響研究［J］.環(huán)境科學(xué)，2006，27（12）:2407－2414.

［10］龍?zhí)煊澹撼５?，李繼承，等.基于SLURP模型和輸出系數(shù)法的三峽庫(kù)區(qū)非點(diǎn)源氮磷負(fù)荷預(yù)測(cè)［J］.環(huán)境科學(xué)學(xué)報(bào)，2008，28（3）:574－581.

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

［12］耿潤(rùn)哲，王曉燕，吳在興，等.北運(yùn)河下游不同土地利用非點(diǎn)源污染負(fù)荷估算［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2012，31（7）:1412－1420.

［13］Johnes P J，Heathwaite A L.Modeling the impact of land use change on water quality in agricultural catchments［J］.Hydrological Proeesses，1997，11（3）:269－286.

［14］莊永濤.渭河臨潼斷面以上流域非點(diǎn)源總氮負(fù)荷研究［D］.陜西 西安:西安理工大學(xué)，2002.

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

［16］梁常德，龍?zhí)煊?，李繼承，等.三峽庫(kù)區(qū)非點(diǎn)源氮磷負(fù)荷研究［J］.長(zhǎng)江流域資源與環(huán)境，2001，16（1）:26－30.

［17］段亮，段增強(qiáng)，夏四清.太湖旱地非點(diǎn)源污染定量化研究［J］.水土保持通報(bào)，2006，26（6）:40－43.

［18］李根，毛峰.我國(guó)水土流失型非點(diǎn)源污染負(fù)荷及其經(jīng)濟(jì)損失評(píng)估［J］.中國(guó)水土保持，2008（2）:9－11.

［19］Beaulac M N，Reckhow K H.An examination of land use nutrient export relationships［J］.Water Resources Bulletin，1982，18（6）:1013－1024.

［20］Frink C R.Estimating nutrient exports to estuaries［J］.Journal of Environmental Quality，1991，20（4）:717－724.

［21］陸建忠，陳曉玲，肖靖靖，等.改進(jìn)的輸出系數(shù)法在農(nóng)業(yè)污染源估算中的應(yīng)用［J］.華中師范大學(xué)學(xué)報(bào)，2012，46（3）:374.

［22］焦艷平，馬香玲，艾萬(wàn)祥，等.崗南、黃壁莊水庫(kù)上游流域非點(diǎn)源污染負(fù)荷估算［J］.南水北調(diào)與水利科技，2010，8（2）:96.

［23］曹彥龍，李崇明，郭勁松，等.重慶三峽庫(kù)區(qū)非點(diǎn)源污染來(lái)源分析及負(fù)荷計(jì)算［J］.重慶建筑大學(xué)，2007，29（4）:2.

［24］沈珍瑤，劉瑞民，葉閩，等.長(zhǎng)江上游非點(diǎn)源污染特征及其變化規(guī)律［M］.北京:科學(xué)出版社，2008:348－358.

［25］韓明明.云南省昆明祿勸彝族苗族自治縣［OL］.（2011－01－25）［2013－05－26］.云南:www.ynszxc.gov.cn，2009.