程玲玲，夏 峰

(1．云南農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，云南 昆明 650201;2．云南省環(huán)境科學(xué)研究院，云南 昆明 650034)

研究表明，我國(guó)是世界上污水排放量最大的國(guó)家，也是污水排放速度最快的國(guó)家。當(dāng)前，作為保護(hù)水環(huán)境、防治水體污染的一個(gè)重要手段——水污染物總量控制已經(jīng)成為當(dāng)前的研究熱點(diǎn)之一。污染物總量如何公平、合理的分配，是污染總量控制的核心問(wèn)題。隨著科學(xué)技術(shù)的進(jìn)步與發(fā)展，總量分配的方法日漸增多，本文將著重介紹國(guó)內(nèi)外總量分配的發(fā)展?fàn)顩r，列舉國(guó)內(nèi)污染物總量分配的一般原則，并分析水污染總量分配兩種常用的具體方法。

1 背景

據(jù)《全國(guó)環(huán)境公報(bào)》顯示，全國(guó)廢水排放量呈逐年增長(zhǎng)趨勢(shì)。2005年，全國(guó)廢水排放量524.5億t;2006年，廢水排放量536.8億t;2007年，廢水排放量556.8億t;2008年比上年增加2.7%，全國(guó)廢水排放總量為571.7億t;2009年，全國(guó)廢水排放589.2億t，比上年增加3.0%。從2009年《中國(guó)環(huán)境狀況公報(bào)》看，在26個(gè)國(guó)控重點(diǎn)湖泊(水庫(kù))中，滿足Ⅱ類水質(zhì)的只有1個(gè)，Ⅲ類的5個(gè)，Ⅳ類的6個(gè)，Ⅴ類的5個(gè)，劣Ⅴ類的9個(gè)，占到了湖泊總數(shù)的34.6%。從這些數(shù)據(jù)不難看出，開展水污染物總量控制工作已經(jīng)刻不容緩。

2 總量分配的一般原則

目前，總量分配的方法很多，根據(jù)立足點(diǎn)和側(cè)重點(diǎn)的不同，分配的一般原則常見有以下兩個(gè)不同方向的闡述，一個(gè)方向是:等比例分配原則，費(fèi)用最小分配原則，按貢獻(xiàn)率削減排放量的分配原則，根據(jù)污染范圍和程度大小的分配原則，按污染物毒性大小承擔(dān)污染責(zé)任分擔(dān)率的原則，按企業(yè)污染治理先進(jìn)性考慮污染責(zé)任分擔(dān)率和削減率的分配原則;另一方向包括可持續(xù)性原則，公平性原則，效益性原則，技術(shù)可行性原則，方案可操作性原則，非經(jīng)濟(jì)要素標(biāo)準(zhǔn)原則，清潔生產(chǎn)原則，先易后難原則，不重復(fù)削減原則，重點(diǎn)控制原則，集中控制原則和公正原則。本人認(rèn)為，總量分配的一般原則可以列為以下幾條:

(1)等比例分配原則

等比例分配，指的是各污染源等比例分擔(dān)排放責(zé)任，以承認(rèn)污染源排污現(xiàn)狀為前提，把總量控制確定的允許排污總量等比例分配到源。

(2)重點(diǎn)控制原則

二八定律認(rèn)為一個(gè)關(guān)鍵的小的誘因、投入和努力，通?？梢缘玫酱蟮慕Y(jié)果、產(chǎn)出和酬勞。遵循二八定律，應(yīng)當(dāng)首先對(duì)重點(diǎn)排污單位進(jìn)行容量總量控制。

(3)先易后難原則

首先對(duì)濃度和行業(yè)總量未達(dá)標(biāo)企業(yè)進(jìn)行總量削減，在濃度達(dá)標(biāo)排放之后，再對(duì)總量負(fù)荷的削減量進(jìn)行分配。

(4)集中控制原則

首先考慮對(duì)位置臨近、污染物種類相同的污染源實(shí)行集中控制，然后再將排污量余量分配給其他的污染源。

(5)費(fèi)用最小分配原則

即經(jīng)濟(jì)優(yōu)化規(guī)劃分配原則，是一種以經(jīng)濟(jì)模型為基礎(chǔ)，以達(dá)到經(jīng)濟(jì)利益最大化為準(zhǔn)則，以治理費(fèi)用最小為目標(biāo)函數(shù)，以環(huán)境目標(biāo)值作為約束條件，使系統(tǒng)的污染治理投資費(fèi)用總和最小的原則。

(6)按污染物毒性大小承擔(dān)污染責(zé)任分擔(dān)率的原則

指在排污總量或排污責(zé)任分擔(dān)率的分配過(guò)程中，對(duì)毒性大、危害嚴(yán)重的危險(xiǎn)污染物應(yīng)提高其治污責(zé)任，加強(qiáng)其污染責(zé)任分配比例的原則。

(7)按貢獻(xiàn)率削減排放量的分配原則

即按各個(gè)污染源對(duì)總量控制區(qū)域內(nèi)水質(zhì)影響程度大小、污染物貢獻(xiàn)率的大小來(lái)削減污染負(fù)荷，對(duì)水質(zhì)影響大的污染源要多削減，反之則少削減。它體現(xiàn)了每個(gè)排污者平等共享水環(huán)境容量資源，同時(shí)也平等承擔(dān)超過(guò)其允許負(fù)荷量的責(zé)任。

(8)根據(jù)污染范圍和程度大小的分配原則

這種分配原則是考慮到在功能區(qū)或污染控制單元內(nèi)污染源的位置不同，有的污染源污染距離長(zhǎng)、面積大，有的污染源污染距離短、面積小，造成各污染源的污染影響的范圍和程度不同，從而將其作為污染責(zé)任分擔(dān)率的重要因素列入在內(nèi)。

(9)按企業(yè)污染治理先進(jìn)性考慮污染責(zé)任分擔(dān)率和削減率的分配原則

由于各企業(yè)污染治理技術(shù)的先進(jìn)性程度不同，造成各企業(yè)的污染排放總量及濃度對(duì)環(huán)境的影響不同，治污成本也存在差異，所以也將它作為污染責(zé)任分擔(dān)率和削減率的重要因素之一。

3 總量分配方法研究現(xiàn)狀

3.1 國(guó)外總量分配發(fā)展與現(xiàn)狀

“總量控制”這一概念由美國(guó)最早提出。1948年，美國(guó)國(guó)會(huì)制定了《聯(lián)邦水污染控制法》(簡(jiǎn)稱為《水污染法》)。1958年，日本開始實(shí)施《水質(zhì)保護(hù)法》、《工業(yè)污水限制法》等法案，在國(guó)家環(huán)境管理上經(jīng)歷了以“稀釋”、“架高”等為主要措施的早期限制時(shí)期和以濃度控制為核心的“單打一”治理時(shí)期。

美國(guó)也是最早關(guān)注非點(diǎn)源污染治理的國(guó)家。1972年修訂后的《聯(lián)邦水污染控制法》首次明確提出控制非點(diǎn)源污染。在此同時(shí)，美國(guó)開始在全國(guó)范圍內(nèi)實(shí)行水污染物排放許可證制度，即所謂的“排放抵消”政策和“泡泡政策”，隨后，為進(jìn)一步做好總量控制工作，美國(guó)提出了季節(jié)總量控制方法、“變量總量控制”，以及在污染源之間進(jìn)行污染負(fù)荷對(duì)換制度，其中包括“點(diǎn)源對(duì)換法”和“點(diǎn)源—非點(diǎn)源兌換法”。

總量控制相關(guān)法律制度出臺(tái)實(shí)施之后，國(guó)外研究學(xué)者在水環(huán)境總量分配的模型及分配方式方面展開了大量研究。如Thomannan and Sobel等，將污水處理費(fèi)用，即目標(biāo)函數(shù)線性化后，運(yùn)用線性規(guī)劃方法對(duì)確定性條件下優(yōu)化模型進(jìn)行分配計(jì)算;Converse，Ecker以確定水質(zhì)為約束條件，以區(qū)域污水處理費(fèi)用之和最小為目標(biāo)函數(shù)建立優(yōu)化模型，再運(yùn)用線性規(guī)劃、動(dòng)態(tài)規(guī)劃或整體規(guī)劃等優(yōu)化技術(shù)，實(shí)現(xiàn)污染物在排污口間的優(yōu)化分配;Lohani and Thanh在河流水質(zhì)規(guī)劃研究中建立了概率約束的規(guī)劃模型;Ellis采用嵌入概率約束條件的方式構(gòu)建了一個(gè)新的隨機(jī)水質(zhì)優(yōu)化模型，和以往水質(zhì)優(yōu)化模型的最大區(qū)別在于該模型中不僅河水流量是隨機(jī)變量，河段起始斷面BOD和DO、廢水中BOD和DO、好氧系數(shù)、復(fù)氧系數(shù)等也被視為隨機(jī)變量;Cardwell and Eills基于參數(shù)和模型的不確定性，對(duì)多點(diǎn)源的污染負(fù)荷分配進(jìn)行了研究，開發(fā)了隨機(jī)動(dòng)態(tài)規(guī)劃費(fèi)用最小模型;Li S.Y.等在考慮了河流斷面的橫向混合不均性基礎(chǔ)上，運(yùn)用優(yōu)化模型確定各排污口在給定水質(zhì)標(biāo)準(zhǔn)下的允許排放量;Catherine L K ling和Jinhuazhao基于污染物的不同區(qū)域特性，對(duì)不同分配方式的長(zhǎng)期效率進(jìn)行了分析;Peter Cramton和Suzi Kerr和其他一些研究者對(duì)拍賣方式和“祖父制”(grandfathering即根據(jù)歷史產(chǎn)量或排污量等因素?zé)o償分配)兩種主要分配方式的優(yōu)缺點(diǎn)進(jìn)行了比較，提出了各自的結(jié)論。

3.2 國(guó)內(nèi)總量分配發(fā)展與現(xiàn)狀

20世紀(jì)70年代末，我國(guó)開始進(jìn)行總量控制的研究工作，以1980年制定的第一松花江BOD總量控制標(biāo)準(zhǔn)為先導(dǎo)，進(jìn)行了最早的探索和實(shí)踐。

“六五”期間，上海開始在黃浦江上游水資源保護(hù)地區(qū)實(shí)行以水污染物排放總量控制為目的的排污許可證制度?！捌呶濉睍r(shí)期，上海頒布了《上海市黃浦江上游水源保護(hù)條例》，首次在地方立法中對(duì)總量控制制度進(jìn)行了規(guī)定。1988年3月，國(guó)家環(huán)保局下達(dá)了《關(guān)于以總量控制為核心的〈水污染物排放許可證管理暫行辦法〉和開展排放許可證試點(diǎn)工作的通知》，標(biāo)志著我國(guó)開始進(jìn)入通過(guò)總量控制強(qiáng)化水環(huán)境管理的新階段;1991年開始，國(guó)家環(huán)保局在16個(gè)城市進(jìn)行了排放污染物總量控制和許可證制度的試點(diǎn)工作，取得了一定的效果。

“九五”期間，全國(guó)人大通過(guò)的《國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展“九五”計(jì)劃和2010年遠(yuǎn)景目標(biāo)綱要》(1996)中，將污染物排放總量控制確立為我國(guó)環(huán)境保護(hù)的一項(xiàng)重大舉措。同一年，國(guó)家制定了《“九五”期間全國(guó)主要污染物排放總量控制計(jì)劃》，并逐年下達(dá)年度計(jì)劃，對(duì)二氧化硫、煙塵、粉塵、工業(yè)固體廢物、化學(xué)需氧量、石油類、氰化類、砷、汞、鉛、鎘、六價(jià)鉻12項(xiàng)污染物排放實(shí)行總量控制;“十五”期間，總量控制制度的管理范圍縮減為二氧化硫、塵 (煙塵和工業(yè)粉塵)、化學(xué)需氧量、氨氮、工業(yè)固體廢物5種污染物;“十一五”期間，總量控制制度的管理范圍進(jìn)一步縮減到二氧化硫、化學(xué)需氧量。2007年，我國(guó)二氧化硫和化學(xué)需氧量排放總量分別下降4.66%和3.14%，首次實(shí)現(xiàn)雙下降;在“十一五”化學(xué)需氧量和二氧化硫兩項(xiàng)主要污染物的基礎(chǔ)上，“十二五”期間國(guó)家將氨氮和氮氧化物納入總量控制體系，同時(shí)污染源普查口徑的農(nóng)業(yè)源也納入總量控制范圍。

從“六五”至今30多a的研究成果看，總量控制的研究主要是圍繞總量分配原則和分配方法展開的。污染物允許排放量分配則多是在目標(biāo)總量控制或水環(huán)境容量總量控制基礎(chǔ)上，基于公平合理原則或效率原則進(jìn)行的分配。如徐華君等評(píng)價(jià)了同等百分比削減分配方法和最小處理費(fèi)用法，探討了兼顧效益與公平的新分配思路;從效率和公平出發(fā)，吳亞瓊研究了初始排污權(quán)分配的一種協(xié)商仲裁制，對(duì)有關(guān)各方面的行為和結(jié)果以及結(jié)果的經(jīng)濟(jì)效率和公平性進(jìn)行了分析;方秦華等人提出了以經(jīng)濟(jì)總量為基礎(chǔ)按比例分配環(huán)境容量的模式;立足于效率原則，王亮等提出了行業(yè)水環(huán)境經(jīng)濟(jì)綜合指數(shù)的概念和確定方法，并建立了對(duì)目標(biāo)區(qū)域水污染物總量行業(yè)優(yōu)化分配模型;岳剛等采用定額達(dá)標(biāo)分配、治理費(fèi)用最小，及按現(xiàn)有排污總量進(jìn)行達(dá)標(biāo)分配的方法對(duì)水污染物進(jìn)行了總量分配;通過(guò)分析總量分配的關(guān)鍵技術(shù)與核心問(wèn)題，農(nóng)家等針對(duì)等比例分配方法、費(fèi)用最小分配方法、按貢獻(xiàn)率削減排放量分配方法以及數(shù)學(xué)規(guī)劃分配方法進(jìn)行了初步探討，并研究了總量分配技術(shù)中的公平性問(wèn)題，為入海污染物的總量分配提供科學(xué)參考。

3.3 總量分配的具體方法

3.3.1 層次分析法

層次分析法 (Analytichi－h(huán)yProcess，簡(jiǎn)稱AHP方法)是美國(guó)運(yùn)籌學(xué)家A.L.Saaty在20世紀(jì)70年代提出的，是一種對(duì)方案的多指標(biāo)系統(tǒng)進(jìn)行分析，將決策者對(duì)復(fù)雜系統(tǒng)的決策思維過(guò)程模型化、數(shù)量化、層次化、結(jié)構(gòu)化的決策方法。層次分析法運(yùn)用在流域水環(huán)境污染物總量分配可以分為以下幾個(gè)步驟:

(1)調(diào)查流域內(nèi)各分區(qū)經(jīng)濟(jì)、社會(huì)、環(huán)境和技術(shù)現(xiàn)狀。

(2)收集總量分配評(píng)價(jià)指標(biāo)、構(gòu)建指標(biāo)體系框架。

(3)確定層次結(jié)構(gòu)模型，設(shè)計(jì)專家咨詢方案，選擇專家并進(jìn)行專家咨詢。

(4)統(tǒng)計(jì)咨詢結(jié)果，計(jì)算分析層次。

(5)確定各單項(xiàng)評(píng)價(jià)指標(biāo)相對(duì)總目標(biāo)權(quán)重。

(6)確定各分區(qū)相對(duì)總目標(biāo)組合權(quán)重。

(7)確定排污總量分配方案。

層次分析法在我國(guó)已經(jīng)被廣泛研究應(yīng)用于污染物總量分配這一領(lǐng)域，如在建立定性與定量相結(jié)合的多指標(biāo)決策的排污總量分配層次結(jié)構(gòu)模型的基礎(chǔ)上，李如忠等將其運(yùn)用在安徽合肥市區(qū)域水污染物排放總量的分配;孫秀喜等人利用這種模型結(jié)合矩陣，對(duì)各分區(qū)排污總量進(jìn)行分?jǐn)?，結(jié)合實(shí)例進(jìn)行了計(jì)算分析;范平等以南方某市化學(xué)需氧量總量控制為例說(shuō)明了這種模型的具體運(yùn)用。2011年，李如忠在《基于Vague集的水污染負(fù)荷分配方法》中從水污染負(fù)荷分配涉及的社會(huì)、經(jīng)濟(jì)、環(huán)境、技術(shù)與管理等因素出發(fā)，篩選16個(gè)評(píng)價(jià)指標(biāo)，構(gòu)建了具有多層次結(jié)構(gòu)的水污染負(fù)荷分配評(píng)價(jià)指標(biāo)體系。

3.3.2 基尼系數(shù)修正法

基尼系數(shù) (Gini coefficient)是意大利經(jīng)濟(jì)學(xué)家基尼 (Corrado Gini，1884－1965)在1922年提出的，用于定量測(cè)定收入分配差異程度?；嵯禂?shù)是一個(gè)比例數(shù)值，在0和1之間，是國(guó)際上用來(lái)綜合考察居民內(nèi)部收入分配差異狀況的一個(gè)重要分析指標(biāo)。基尼系數(shù)修正法運(yùn)用在流域水環(huán)境污染物總量分配可以分為以下幾個(gè)步驟:

(1)劃分子流域的污染控制區(qū)，利用SWAT(陸域面源污染物輸移模擬)和EFDC(環(huán)境流體動(dòng)力學(xué))模型，核算并確定不同子流域的污染負(fù)荷及流域的水環(huán)境容量。

(2)確定初始總量分配方案。根據(jù)流域的現(xiàn)狀和未來(lái)發(fā)展模式以及各子流域的功能定位，確定各子流域的初始分配值。

(3)確定基本參數(shù) (主要包括各子流域人口或GDP、產(chǎn)值、水資源量等)，計(jì)算公平指數(shù)Gi。

(4)將Gi按由大到小的順序排列，將各子流域根據(jù)順序重新編號(hào)。繼而將對(duì)應(yīng)的參數(shù)按照編號(hào)重新整理出來(lái)。

(5)按重排的順序，依照相應(yīng)的公式，計(jì)算各子流域人口 (GDP)等指標(biāo)的累計(jì)百分比和污染負(fù)荷分配的累計(jì)百分比。

(6)繪制洛倫茲曲線和絕對(duì)公平曲線。

(7)計(jì)算出基尼系數(shù)。

(8)依據(jù)計(jì)算的基尼系數(shù)，評(píng)估分配的公平性。如果分配不是公平的，則從步驟 (2)開始重新調(diào)整，反之則分配過(guò)程結(jié)束。

在研究成果方面，學(xué)者王麗瓊從基尼系數(shù)的經(jīng)濟(jì)學(xué)內(nèi)涵出發(fā)，將基尼系數(shù)的概念和意義應(yīng)用于環(huán)境分析中，并以全國(guó)2006年各地區(qū)的水污染物總量分配為例，應(yīng)用基尼系數(shù)法分析了人口、國(guó)內(nèi)生產(chǎn)總值、水資源量指標(biāo)對(duì)水污染物總量分配的影響;樊華等在江西五大水系流域的水污染物總量分配方案的評(píng)估和調(diào)整中運(yùn)用了基尼系數(shù)法，使方案更加公平合理;趙娟等構(gòu)建重點(diǎn)城市污染物總量分配的基尼系數(shù)法，評(píng)價(jià)了全國(guó)重點(diǎn)城市污染物總量分配的合理性公平性。

另外，研究者盛虎參考美國(guó)TMDLs的分配方式 (在分配過(guò)程中預(yù)留一部分容量不作分配，以促進(jìn)水質(zhì)改善)，運(yùn)用層次分析法確定分配比例，運(yùn)用基尼系數(shù)法判別分配公平性，采取層次分析法與基尼系數(shù)法交換反饋的方法，研究了基于人口和經(jīng)濟(jì)的流域容量總量分配方案。

4 污染物總量分配方法中存在的問(wèn)題

基于公平性原則的總量分配方法的問(wèn)題在于實(shí)用與理論脫節(jié):要么停留在理論研究上，僅給出堅(jiān)持公平性原則需要遵循的原理和考慮的因素，而缺乏一套較完整的體系;要么給出了一套較完整的體系，卻過(guò)于理論化使得變量的選擇難以定量化，分配方案過(guò)于復(fù)雜而難以在實(shí)際中加以應(yīng)用，最終難以達(dá)到預(yù)期目的;或者片面強(qiáng)調(diào)某一方面的公平性，忽略了其他方面而致使真正的公平性并不能體現(xiàn)出來(lái)。

基于效率原則的總量分配方法，大多關(guān)注于如何使經(jīng)濟(jì)成本最小、治污費(fèi)用最小，忽視了總量分配過(guò)程中應(yīng)該考慮到的社會(huì)、環(huán)境、技術(shù)等因素。這種方法要達(dá)到經(jīng)濟(jì)成本最低，治污得力而花費(fèi)較少的企業(yè)承擔(dān)相應(yīng)較多的污染責(zé)任分擔(dān)率，分配的不合理性很容易引起企業(yè)之間的糾紛從而導(dǎo)致總量分配實(shí)施程序很難順利進(jìn)行。

5 結(jié)語(yǔ)

從“六五”的開始實(shí)行到“十二五”的規(guī)劃，污染物總量控制與分配問(wèn)題一直是我國(guó)環(huán)保工作中不可忽視且越來(lái)越關(guān)鍵的環(huán)節(jié)。本人認(rèn)為:在相對(duì)公平合理的基礎(chǔ)上，尋求環(huán)境、經(jīng)濟(jì)、技術(shù)系統(tǒng)整體最優(yōu)配置，結(jié)合多種污染物總量分配方法的優(yōu)點(diǎn)，同時(shí)把GIS、計(jì)算機(jī)等先進(jìn)科技手段融入到污染物總量控制中來(lái)，形成一套理論與實(shí)用較好結(jié)合的綜合體系是今后研究的重點(diǎn);對(duì)于非點(diǎn)源污染因排放特征復(fù)雜形成不易于把握的許多不確定性也需要深入分析。

［1］韓佳明．環(huán)境污染總量控制與環(huán)境監(jiān)測(cè)實(shí)用技術(shù)手冊(cè) ［M］．北京:中國(guó)環(huán)境科學(xué)出版社，2007．

［2］李如忠．區(qū)域水污染物排放總量分配方法研究 ［J］．環(huán)境工程，2002，20(6):61－63．

［3］田仁生，張治忠，鄒首民．關(guān)于污染物排放總量控制工作的若干思考和建議 ［J］．上海環(huán)境科學(xué)，2003，7:277－290．

［4］劉旭，秦南茜．我國(guó)污染防治法中的總量控制制度概述 ［J］．法制與社會(huì)，2011，(4):41－43．

［5］李如忠，錢家忠，汪家權(quán)．水污染物允許排放總量分配研究方法 ［J］．水力學(xué)報(bào)，2003，(5):112－115．

［6］王宏，孟凡宇，邢妍．水環(huán)境承載力與污染物總量分配研究進(jìn)展 ［J］．遼寧大學(xué)學(xué)報(bào) (自然科學(xué)版)，2009，6(4):364－367．

［7］宗永臣，張建新．水污染總量分配方法研究綜述 ［J］．內(nèi)江科技，2008，29(3):74－75．

［8］農(nóng)家，王金坑，陳克亮，等．入海污染物總量分配技術(shù)研究初探 ［J］．環(huán)境保護(hù)科學(xué)，2009，35(5):45－48．

［9］李如忠，舒琨．基于基尼系數(shù)的水污染負(fù)荷分配模糊優(yōu)化決策模型［J］．環(huán)境科學(xué)學(xué)報(bào)，2010，30(7):1518－1526．

［10］李如忠，舒琨．基于Vague集的水污染負(fù)荷分配方法 ［J］．水利學(xué)報(bào)，2011，42(2):127－135．

［11］鄧義祥，孟偉，鄭炳輝，等．基于響應(yīng)場(chǎng)的線性規(guī)劃方法在長(zhǎng)江口總量分配計(jì)算中的應(yīng)用［J］．環(huán)境科學(xué)研究，2009，22(9):995－1000．

［12］王乖虎，萬(wàn)繼偉，辛國(guó)興．區(qū)域污染源排污總量分配方案初探 ［J］．環(huán)境科學(xué)與管，2010，35(2):9－12．

［13］楊玉峰，傅國(guó)偉．區(qū)域差異與國(guó)家污染物排放總量分配［J］．環(huán)境科學(xué)學(xué)報(bào)，2001，21(2):129－133．

［14］林高松，李適宇，江峰．基于公平區(qū)間的污染物允許排放量分配方法［J］．水利學(xué)報(bào)，2006，37(1):52－56．

［15］方秦華，張珞平，王佩爾．象山港海域環(huán)境容量的二步分配法［J］．廈門大學(xué)學(xué)報(bào) (自然科學(xué)版)，2004，43(增刊):217－220．

［16］盛虎，李娜，郭懷成，等．流域容量總量分配及排污交易潛力分析［J］．環(huán)境科學(xué)學(xué)報(bào)，2010，30(3):655－663．

［17］范平，吳純德，張帆，等．南方某城市水污染物總量控制分配研究［J］．水資源保護(hù)，2009，25(4):20－39．

［18］孫秀喜，馮耀奇，丁和義．河道污染物總量分配模型的建立及分析方法研究 ［J］．地下水 ，2005，27(6):427－469．

［19］王紫雯，陳昌軍．湖泊水庫(kù)入流河道排污總量分配方法的探討 ［J］．浙江大學(xué)學(xué)報(bào)，2001，35(5):534－539．

［20］樊華，陳然．基尼系數(shù)法在水污染物總量分配中的應(yīng)用［J］．安徽農(nóng)業(yè)科學(xué)，2009，37(24):11678－11680，11685．

［21］王麗瓊，基于公平性的水污染物總量分配基尼系數(shù)分析［J］．生態(tài)環(huán)境，2008，17(5):1796－1801．

［22］黃輝柏，羅敏，毛燁，等．江西省九江市水環(huán)境容量測(cè)算和總量分配研究 ［J］．江西能源，2007，(4):67－70．

［23］郭森，韓保新，楊靜．納污海域水環(huán)境容量計(jì)算與總量分配方法研究 －以欽州灣為例 ［J］．環(huán)境科學(xué)與技術(shù)，2008(29):19－22．

［24］趙娟，羅宏，薛婕，等．全國(guó)重點(diǎn)城市工業(yè)廢水污染物排放總量分配［J］．環(huán)境保護(hù)科學(xué)，2010，36(2):8－10．

［25］胡康萍，許振成．水體污染物允許排放總量分配方法研究［J］．中國(guó)環(huán)境科學(xué)，1991，(6):447－452．

［26］錢駿，肖杰，蔣夏．四川省水污染物總量分配測(cè)算研究［J］．西華大學(xué)學(xué)報(bào) (自然科學(xué)版)，2008，27(6):16－18．

［27］吳亞瓊，趙勇，吳相林，等．污染物排放總量分配的機(jī)制設(shè)計(jì)方法研究 ［J］．管理工程學(xué)報(bào)，2004，(4):65－68．

［28］徐華君，徐百福．污染物允許排放總量的公平協(xié)調(diào)思路與方法 ［J］．新疆大學(xué)學(xué)報(bào)，1996，13(3):86－89．

［29］岳剛，呂焰．撫順市水污染物總量分配探討［J］．遼寧城鄉(xiāng)環(huán)境科技，2007，27(4):41－59．

［30］趙華林，黃小贈(zèng)，張震宇．日本水污染物總量控制技術(shù)政策及對(duì)小城鎮(zhèn)分散型污水處理的思考 ［J］．環(huán)境保護(hù)，2009，425(15):70－72．

［31］馮金鵬，吳洪壽，趙帆．水環(huán)境污染總量控制回顧、現(xiàn)狀及發(fā)展探討［J］．南水北調(diào)與水利科技，2004， (1):1672－1683．

［32］李摯萍．中國(guó)排污許可制度立法研究——兼談環(huán)境保護(hù)基本制度之間協(xié)調(diào)［A］．環(huán)境法治與建設(shè)和諧社會(huì)——2007年全國(guó)環(huán)境資源法學(xué)研討會(huì)論文集［B］．

［33］封金利，楊維，施爽，等．水污染物總量分配方法研究［J］．環(huán)境保護(hù)與循環(huán)經(jīng)濟(jì)，2010，6(30):34－37．

［34］ Ellis．J．H．Stochastic water quality optimization using in bedded chance constrains［J］．Water Resources Research，1987，(2):2227－2238．

［35］ Cardwell．H．a(chǎn)nd Ells．H．Stochastic dynamic programming models for water quality management［J］．Water Resources Research，1993，29(4):803－813．

［36］ Thomann R V，Sobel M．S．Estuarine．Water Quality Management and Forecasting［J］．Journal of Sanitary Engineering Division，1964，89(5):9－36．

［37］ K ling C L，Zhao J H．On the long－run efficiency of auctioned vs．freepermits ［J］．Economics Letters，2000，69(2):235－238．

［38］ Cramton P，Kerr S．Trade able carbon permit auctions－How and why to auction not grandfather［J］．Energy Policy，2002，30(4):333－345．

［39］ Jensen Jesper，Rasmussen Tobias．Allocation of CO2emission permits:A general equilibrium analysis of policy instruments［J］．Journal of Environmental Economics＆ Management，2000，40(2):111－ 136．

［40］Jesper J，Tobias R．Allocation of CO2emission permits:a general equilibrium Analysis of policy instruments［J］．Journal of EnvironmentalEconomics＆Management， 2000， 40 (2):111－136．

［41］ US EPA．Handbook for Developing Watershed TMDLs［R］．Washington D．C．:United States Environmental Protection A-gency，Office of Water(4503F)，2008:1－10．