白 輝 高 偉 陳 巖# 吳悅穎
(1.環(huán)境保護(hù)部環(huán)境規(guī)劃院,北京 100012;2.云南大學(xué)生態(tài)學(xué)與環(huán)境學(xué)院,云南 昆明 650091)
在大氣環(huán)境受到嚴(yán)重挑戰(zhàn)的今天,水環(huán)境污染問(wèn)題同樣不能輕視,其仍是目前我國(guó)面臨的主要環(huán)境問(wèn)題之一。造成目前這些環(huán)境問(wèn)題的根本原因是部分有限的地理空間上承載了過(guò)量的能源消費(fèi)、工業(yè)生產(chǎn)及其他各種人類(lèi)活動(dòng)等,社會(huì)經(jīng)濟(jì)承載現(xiàn)狀已達(dá)到或遠(yuǎn)超環(huán)境承載能力。面對(duì)如此嚴(yán)峻的形勢(shì),黨中央和國(guó)務(wù)院高度重視,相繼出臺(tái)相關(guān)文件要求摸清水資源、水環(huán)境承載力現(xiàn)狀,建立監(jiān)測(cè)預(yù)警機(jī)制。國(guó)務(wù)院印發(fā)的《水污染防治行動(dòng)計(jì)劃》提出要充分考慮水資源、水環(huán)境承載能力,以水定城、以水定地、以水定人、以水定產(chǎn);建立水資源、水環(huán)境承載能力監(jiān)測(cè)評(píng)價(jià)體系,實(shí)行承載能力監(jiān)測(cè)預(yù)警,已超過(guò)承載能力的地區(qū)要實(shí)施水污染物削減方案,加快調(diào)整發(fā)展規(guī)劃和產(chǎn)業(yè)結(jié)構(gòu);到2020年,組織完成市、縣域水資源、水環(huán)境承載能力現(xiàn)狀評(píng)價(jià)[1]。水環(huán)境承載力的計(jì)算及管理是區(qū)域水環(huán)境保護(hù)規(guī)劃工作的一個(gè)重要基礎(chǔ),通過(guò)研究社會(huì)經(jīng)濟(jì)發(fā)展的“環(huán)境承載上限”,以環(huán)境容量為約束,科學(xué)評(píng)估水環(huán)境承載力,不僅可為區(qū)域項(xiàng)目環(huán)境影響評(píng)價(jià)和環(huán)境基礎(chǔ)設(shè)施建設(shè)提供參考數(shù)據(jù),為地區(qū)經(jīng)濟(jì)發(fā)展模式、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與轉(zhuǎn)型、優(yōu)化空間布局等提供環(huán)境保護(hù)基礎(chǔ)數(shù)據(jù)支持,還可為流域、區(qū)域間實(shí)施污染物排放總量控制、排污許可證制度和排污權(quán)交易等現(xiàn)代環(huán)境管理制度提供重要的科學(xué)依據(jù)。因此,從環(huán)境質(zhì)量管理的實(shí)際需求出發(fā),亟需開(kāi)展水環(huán)境承載力評(píng)價(jià)研究,為社會(huì)經(jīng)濟(jì)發(fā)展和環(huán)境管理提供決策依據(jù)與科學(xué)支撐。
近年來(lái),我國(guó)學(xué)者在水環(huán)境承載力的理論與實(shí)踐方面積極探索并取得了很大進(jìn)展。目前,在學(xué)術(shù)界,水環(huán)境承載力的概念尚未達(dá)成統(tǒng)一,尚無(wú)公認(rèn)、統(tǒng)一的界定,與水資源承載力、生態(tài)承載力和水生態(tài)承載力等概念存在較大的交叉,使得承載力評(píng)估缺乏固定標(biāo)準(zhǔn)[2-3]。目前,關(guān)于水環(huán)境承載力的概念主要分為3類(lèi):第1類(lèi)主要闡述水環(huán)境能夠支撐的社會(huì)經(jīng)濟(jì)和人口的規(guī)模,如賈振邦[4]在1995年給出的定義——“在一定的自然環(huán)境條件和特定的社會(huì)經(jīng)濟(jì)發(fā)展模式下,區(qū)域水環(huán)境對(duì)社會(huì)發(fā)展的支撐能力”;第2類(lèi)主要闡述水環(huán)境能夠接納的污水和污染物的含量,如前水利部部長(zhǎng)汪恕誠(chéng)[5]在2003年給出的定義——“一定的水域,在其水體能夠繼續(xù)使用并能保持良好的生態(tài)系統(tǒng)的前提下,所能容納污水和污染物的最大能力”;第3類(lèi)主要闡述水環(huán)境能夠接納的污水和污染物的含量、能夠支撐的社會(huì)經(jīng)濟(jì)和人口的規(guī)模,如崔寧等[6]在2010年給出的定義——“在某一水域(包括自然環(huán)境條件和社會(huì)經(jīng)濟(jì)發(fā)展模式),其水環(huán)境系統(tǒng)功能在持續(xù)正常發(fā)揮的前提下,所能夠容納污染物的最大限量及對(duì)人類(lèi)社會(huì)、經(jīng)濟(jì)活動(dòng)可持續(xù)發(fā)展的最大支撐能力”。
在水環(huán)境承載力評(píng)價(jià)方法研究方面,學(xué)者們嘗試采取多種方法對(duì)水環(huán)境承載力進(jìn)行定性與定量的研究[7]。應(yīng)用于水環(huán)境承載力評(píng)價(jià)的方法主要有綜合指標(biāo)體系評(píng)價(jià)法、多目標(biāo)決策優(yōu)化法、系統(tǒng)動(dòng)力學(xué)方法、人工神經(jīng)網(wǎng)絡(luò)法等[8-20]。這些研究大多處于探索階段,學(xué)術(shù)界未形成公認(rèn)的量化評(píng)價(jià)方法。
結(jié)合水環(huán)境承載力的基本理念和本研究的實(shí)際需求,本研究認(rèn)為水環(huán)境承載力是指在一定的水文條件和水質(zhì)目標(biāo)約束下,區(qū)域水體可支撐的一定結(jié)構(gòu)、布局和技術(shù)水平下經(jīng)濟(jì)社會(huì)發(fā)展規(guī)模。為了科學(xué)評(píng)價(jià)水環(huán)境承載力,本研究在現(xiàn)有水環(huán)境承載力研究的基礎(chǔ)上,提出了一種基于環(huán)境容量的水環(huán)境承載力評(píng)價(jià)方法,并且結(jié)合目前水環(huán)境管理需求,提出了一種根據(jù)水環(huán)境承載力結(jié)果進(jìn)行總量控制的方法,以期解決目前水環(huán)境管理面臨的問(wèn)題。
水環(huán)境承載力評(píng)價(jià)的主要目的是判斷區(qū)域經(jīng)濟(jì)社會(huì)活動(dòng)對(duì)水環(huán)境的壓力是否超過(guò)承載力水平。因此,水環(huán)境承載力評(píng)價(jià)的基本內(nèi)容包括壓力計(jì)算、承載力計(jì)算和承載力評(píng)估3個(gè)方面。在本研究中,重點(diǎn)關(guān)注的是人類(lèi)活動(dòng)對(duì)水環(huán)境質(zhì)量的壓力,特別是水體中的COD和氨氮兩項(xiàng)指標(biāo)。在水環(huán)境承載力評(píng)估的計(jì)算中,一定水文條件和水質(zhì)目標(biāo)約束可表達(dá)為某種特征污染物的水環(huán)境容量;一定結(jié)構(gòu)、布局和技術(shù)水平下經(jīng)濟(jì)社會(huì)發(fā)展規(guī)模難以用單一指標(biāo)直接表述,但考慮在水環(huán)境承載力系統(tǒng)中,經(jīng)濟(jì)社會(huì)發(fā)展是通過(guò)污染物排放作用于水環(huán)境系統(tǒng)的,因此可用污染物進(jìn)入水體的量(入河(湖)量)間接表達(dá)經(jīng)濟(jì)社會(huì)的壓力強(qiáng)度。綜上,本研究的水環(huán)境承載力評(píng)估涉及3個(gè)方面內(nèi)容:水環(huán)境容量計(jì)算、污染物入河(湖)量計(jì)算和基于環(huán)境容量的水環(huán)境承載力評(píng)估。
水環(huán)境容量和污染物入河(湖)量計(jì)算方法已較成熟,本研究中不對(duì)其計(jì)算方法進(jìn)行描述,只對(duì)基于其計(jì)算結(jié)果的水環(huán)境承載力評(píng)價(jià)方法進(jìn)行研究。
水環(huán)境承載力評(píng)估主要包括3個(gè)步驟:(1)計(jì)算基于承載率的單要素(污染物)的水環(huán)境承載力;(2)計(jì)算綜合的水環(huán)境承載力;(3)根據(jù)水環(huán)境承載力的等級(jí)劃分結(jié)果,確定水環(huán)境承載力的狀態(tài)。本研究主要針對(duì)水環(huán)境中的污染物COD和氨氮進(jìn)行水環(huán)境承載力的評(píng)估。
由于水污染物種類(lèi)較多,表征水環(huán)境質(zhì)量的指標(biāo)較多,每個(gè)指標(biāo)對(duì)應(yīng)的環(huán)境容量和排放量存在差異,因此有必要針對(duì)單個(gè)指標(biāo)進(jìn)行單要素的評(píng)估。本研究采用承載率指標(biāo)表征各污染物的承載力狀況。承載率是指區(qū)域或流域環(huán)境承載量(各環(huán)境要素指標(biāo)的現(xiàn)狀值)與該區(qū)域或流域環(huán)境承載量閾值(各環(huán)境要素指標(biāo)上限值)的比值。應(yīng)用承載率指標(biāo)進(jìn)行評(píng)價(jià),可以清晰地看出某地區(qū)或流域環(huán)境發(fā)展現(xiàn)狀與理想值的差距,評(píng)價(jià)其環(huán)境承載的壓力現(xiàn)狀。承載率的計(jì)算公式如下:
(1)
式中:CI為承載率;PD為污染物入河(湖)量,t/a;EC為水環(huán)境容量,t/a。
承載率的取值為大于0的實(shí)數(shù),在進(jìn)行指標(biāo)解釋時(shí)不便使用,為此本研究選擇極差法把各項(xiàng)指標(biāo)進(jìn)行歸一化處理。歸一化的結(jié)果為所有的指標(biāo)均在0~100取值,數(shù)值越大表示指標(biāo)所反映的承載狀態(tài)越好。
(2)
式中:Zi為第i個(gè)污染物指標(biāo)的歸一化承載率;CImin、CImax、CIvalue分別為承載率的最小值、最大值和當(dāng)前值。
承載率無(wú)法反映水環(huán)境系統(tǒng)的整體承載水平,水環(huán)境承載力評(píng)估要求將不同指標(biāo)綜合起來(lái),獲得綜合評(píng)估值。本研究借助綜合評(píng)估模型將各個(gè)指標(biāo)的承載率數(shù)值進(jìn)行加權(quán)綜合,獲得研究區(qū)的綜合水環(huán)境承載力,公式如下:
(3)
式中:ECIr為第r個(gè)行政區(qū)的水環(huán)境承載力指數(shù);n為污染物指標(biāo)總數(shù);ωi為第i個(gè)污染物指標(biāo)的權(quán)重。
指標(biāo)賦權(quán)是根據(jù)各個(gè)指標(biāo)重要性分別賦予相應(yīng)權(quán)重的過(guò)程。目前,指標(biāo)賦權(quán)方法有主觀賦權(quán)和客觀賦權(quán)兩種方法,其中主觀賦權(quán)法以層次分析法為主,該方法以專(zhuān)家打分為基礎(chǔ),根據(jù)打分結(jié)果對(duì)指標(biāo)進(jìn)行賦權(quán);客觀賦權(quán)法以熵權(quán)法為典型代表,主要根據(jù)指標(biāo)自身的熵值進(jìn)行賦權(quán),以反映指標(biāo)信息量的重要性。兩種賦權(quán)法各有利弊,但共同的特點(diǎn)是權(quán)重一旦確定,即為恒定值,在整個(gè)評(píng)價(jià)過(guò)程中與評(píng)價(jià)指標(biāo)的取值無(wú)關(guān)[21-22]。在現(xiàn)實(shí)中,同一系統(tǒng)的不同指標(biāo)往往具有互補(bǔ)性,在可持續(xù)發(fā)展理論中被表述為“弱可持續(xù)性”,即系統(tǒng)的不同屬性之間在一定程度上是可以相互補(bǔ)償?shù)?。此外,木桶原理也表明,系統(tǒng)的狀態(tài)常受制于處于最差狀態(tài)的變量。因此,系統(tǒng)中各個(gè)指標(biāo)的權(quán)重應(yīng)當(dāng)是與指標(biāo)自身取值相關(guān)的,對(duì)于取值較差的指標(biāo),應(yīng)當(dāng)賦予更高的權(quán)重。為體現(xiàn)不同指標(biāo)取值變化對(duì)評(píng)價(jià)結(jié)果的影響,本研究使用變權(quán)法對(duì)各個(gè)指標(biāo)進(jìn)行賦權(quán)。該方法以等權(quán)法為基礎(chǔ),考慮指標(biāo)與評(píng)價(jià)標(biāo)準(zhǔn)的偏離度,偏離度越大,其權(quán)重越高,體現(xiàn)了極端變化指標(biāo)對(duì)評(píng)價(jià)結(jié)果的重要影響。變權(quán)法的計(jì)算公式為:
(4)
式中:Zbest為評(píng)價(jià)標(biāo)準(zhǔn)的最優(yōu)值,在本研究中Zbest=1。
目前,采用指數(shù)評(píng)價(jià)的模型一般將指數(shù)等分為不同等級(jí),從而便于評(píng)價(jià)結(jié)果的分析。本研究將水環(huán)境承載力指數(shù)按20組距劃分為5等級(jí),按優(yōu)承載、可承載、輕超載、中超載、重超載命名,并分配5種不同顏色(見(jiàn)表1)。優(yōu)承載說(shuō)明地區(qū)的水環(huán)境面臨的經(jīng)濟(jì)社會(huì)壓力較小,尚有較大的開(kāi)發(fā)空間;重超載表示研究區(qū)當(dāng)前的經(jīng)濟(jì)社會(huì)發(fā)展已嚴(yán)重超出水環(huán)境的支撐能力,其他等級(jí)按線性內(nèi)插方式設(shè)定。
表1 水環(huán)境承載力等級(jí)標(biāo)準(zhǔn)
污染物總量控制是從管理上實(shí)現(xiàn)提出的環(huán)境目標(biāo),但與水質(zhì)和環(huán)境改善之間缺乏直接的聯(lián)系。在區(qū)域污染物排放總量目標(biāo)確定的情況下,如何實(shí)現(xiàn)區(qū)域污染物排放的最優(yōu)配置是總量控制中必須首先解決的重要問(wèn)題。目前,總量分配常用的方法有等比例分配法、貢獻(xiàn)率分配法、基尼系數(shù)法和博弈法等[23-24]。這些方法均從經(jīng)濟(jì)社會(huì)發(fā)展程度或污染貢獻(xiàn)的角度進(jìn)行總量分配,難以有效兼顧環(huán)境和人類(lèi)活動(dòng)的綜合因素,本研究提出一種基于區(qū)域水環(huán)境承載力的總量分配方法。該方法以各區(qū)域的水環(huán)境承載力為基礎(chǔ)進(jìn)行排放量的分配,在水環(huán)境承載力較高的區(qū)域,應(yīng)分配較高允許排放量;反之,在水環(huán)境承載力較低的區(qū)域,應(yīng)減少允許排放量,降低區(qū)域人類(lèi)活動(dòng)的壓力。在該方法體系中,總量分配是通過(guò)二次分配實(shí)現(xiàn)的。在第1層次分配時(shí),根據(jù)研究區(qū)總的污染物總量控制目標(biāo),計(jì)算總削減率,采用等比例削減原則將削減目標(biāo)分配到各行政區(qū),該步驟主要體現(xiàn)分配的公平性原則;在此基礎(chǔ)上,從管理、工程和結(jié)構(gòu)減排3個(gè)角度計(jì)算各行政區(qū)的最大削減能力,進(jìn)而估算最小的可能排放量,以此作為等比例分配的約束,該步驟主要考慮分配結(jié)果的可行性;經(jīng)過(guò)最大削減能力的校正,可能會(huì)出現(xiàn)所有行政區(qū)的排放量加和高于允許排放量的情況,此時(shí)需要進(jìn)一步增加部分區(qū)域的削減量,即二次分配的問(wèn)題。在本研究中,二次分配的主要考慮因素是水環(huán)境承載能力,以現(xiàn)狀評(píng)價(jià)的水環(huán)境承載力為基數(shù),對(duì)存在削減余量的區(qū)域進(jìn)一步分配削減指標(biāo)。基于水環(huán)境承載力的污染物總量分配技術(shù)具體見(jiàn)圖1。
總量分配的具體計(jì)算過(guò)程分為兩個(gè)部分,即第1層次分配和第2層次分配。第1層次分配主要考慮公平性和可行性問(wèn)題,計(jì)算公式如下:
FPrj=Max(EPrj,APrj-PPrj)
(5)
(6)
AWCrjm
(7)
式中:r、j、k、m分別為研究區(qū)內(nèi)的行政區(qū)、污染物種類(lèi)、行業(yè)類(lèi)型、可關(guān)閉行業(yè)類(lèi)型的序號(hào);FPrj、EPrj、APrj、PPrj分別為第r個(gè)行政區(qū)第j種污染物第1層次分配的排放量、等比例削減原則下分配的排放量、實(shí)際排放量、削減最大潛力,mg;TCPj、TAPj分別為區(qū)域第j種污染物實(shí)際排放總量和總量控制目標(biāo),mg;WTVrk為第r個(gè)行政區(qū)第k種行業(yè)的污水排放量,L;WRMArk、WRMIrk分別為第r個(gè)行政區(qū)第k種行業(yè)的污水最大、實(shí)際處理率,%;PCIrjk、PCOrjk分別為第r個(gè)行政區(qū)第j種污染物第k種行業(yè)的污水處理設(shè)備入水、出水質(zhì)量濃度,mg/L;SDRrk、ADRrk分別為第r個(gè)行政區(qū)第k種行業(yè)的污水最大、實(shí)際達(dá)標(biāo)排放率,%;SDCrjk、ADCrjk分別為第r個(gè)行政區(qū)第j種污染物第k種行業(yè)的污水達(dá)標(biāo)、實(shí)際排放質(zhì)量濃度,mg/L;WTVrm為第r個(gè)行政區(qū)第m種可關(guān)閉行業(yè)的污水排放量,L;AWCrjm為第r個(gè)行政區(qū)第j種污染物第m種可關(guān)閉行業(yè)的污水排放質(zhì)量濃度,mg/L。
圖1 基于水環(huán)境承載力的污染物總量分配技術(shù)示意圖Fig.1 The waste load allocation technology roadmap based on water environmental carrying capacity
第2層次分配是在第1層次分配的基礎(chǔ)上進(jìn)行的,主要考慮各地區(qū)的水環(huán)境承載力水平,計(jì)算公式如下:
(8)
式中:OPrj為第r個(gè)行政區(qū)第j種污染物最終分配的排放量,mg。
本研究提出的基于環(huán)境容量的水環(huán)境承載力評(píng)價(jià)方法,方法科學(xué)合理,可操作性強(qiáng),比較符合現(xiàn)實(shí)需求,具有一定的實(shí)踐意義,能對(duì)全國(guó)水環(huán)境承載力現(xiàn)狀評(píng)價(jià)起一定的指導(dǎo)作用?;谒h(huán)境承載力的總量控制方法是對(duì)目前總量控制制度的改進(jìn),能將現(xiàn)有總量控制與地方實(shí)際水環(huán)境承載力現(xiàn)狀很好地結(jié)合起來(lái),能為“十三五”總量控制管理提供科學(xué)技術(shù)方法和思路。
本研究從環(huán)境管理的實(shí)際需求出發(fā),重點(diǎn)關(guān)注的是水環(huán)境承載力中人類(lèi)活動(dòng)對(duì)水環(huán)境質(zhì)量的壓力,而對(duì)水環(huán)境承載力中的人口和社會(huì)經(jīng)濟(jì)規(guī)模的研究將會(huì)在后期的優(yōu)化調(diào)控研究中進(jìn)行,以總的水環(huán)境承載力為約束,兼顧不同區(qū)域環(huán)境容量利用情況,實(shí)現(xiàn)區(qū)域的最大承載能力和水環(huán)境目標(biāo),以達(dá)到不同區(qū)域之間水環(huán)境承載能力的優(yōu)化配置。
[1] 國(guó)務(wù)院.水污染防治行動(dòng)計(jì)劃[M].北京:人民出版社,2015.
[2] 侯麗敏,岳強(qiáng),王彤.我國(guó)水環(huán)境承載力研究進(jìn)展與展望[J].環(huán)境保護(hù)科學(xué),2015,41(4):104-108.
[3] 郭懷成,唐劍武.城市水環(huán)境與社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展對(duì)策研究[J].環(huán)境科學(xué)學(xué)報(bào),1995,15(3):363-369.
[4] 賈振邦.本溪市水環(huán)境承載力及指標(biāo)體系[J].環(huán)境保護(hù)科學(xué),1995,21(2):8-12.
[5] 汪恕誠(chéng).水環(huán)境承載能力分析與調(diào)控[J].水利發(fā)展研究,2002(2):2-6.
[6] 崔寧,梁冬梅.水環(huán)境承載力評(píng)述[J].山西水利,2010,21(12):113-114.
[7] 李川.水環(huán)境承載力量化方法的研究進(jìn)展[J].環(huán)境科學(xué)與管理,2008,33(8):66-70.
[8] 來(lái)雪慧,王小文,徐杰峰,等.基于向量模法的陜南地區(qū)水環(huán)境承載力評(píng)價(jià)[J].水土保持通報(bào),2010,30(2):56-59.
[9] 耿雅妮.基于向量模法的西安市水環(huán)境承載力研究[J].中國(guó)農(nóng)學(xué)通報(bào),2013,29(11):168-172.
[10] 孟繁宇,樊慶鋅,李穎.哈爾濱水環(huán)境承載力評(píng)價(jià)指標(biāo)體系的構(gòu)建[J].黑龍江水利科技,2012,40(3):1-4.
[11] 高偉.湖州市水環(huán)境承載力評(píng)價(jià)[J].水利科技與經(jīng)濟(jì),2009,15(11):947-949.
[12] 江明峰,馬太玲,孫晶.呼和浩特市水環(huán)境承載力綜合評(píng)價(jià)[J].干旱區(qū)資源與環(huán)境,2010,24(9):60-63.
[13] 王莉芳,陳春雪.濟(jì)南市水環(huán)境承載力評(píng)價(jià)研究[J].環(huán)境科學(xué)與技術(shù),2011,34(5):199-202.
[14] 曾現(xiàn)進(jìn),李天宏,溫曉玲.基于AHP和向量模法的宜昌市水環(huán)境承載力研究[J].環(huán)境科學(xué)與技術(shù),2013,36(6):200-205.
[15] 李美榮,鄭欽玉,劉娟,等.基于AHP法的重慶市水環(huán)境承載力研究[J].水利科技與經(jīng)濟(jì),2012,18(5):1-6.
[16] 朱銀銀,柴磊.基于系統(tǒng)動(dòng)力學(xué)的某市水環(huán)境承載力動(dòng)態(tài)變化研究[J].水利科技與經(jīng)濟(jì),2010,16(9):1039-1041.
[17] 邢有凱,余紅,肖楊,等.基于向量模法的北京市水環(huán)境承載力評(píng)價(jià)[J].水資源保護(hù),2008,24(4):1-4.
[18] 張軍.基于向量模法的陜西省水環(huán)境承載力評(píng)價(jià)[J].河南科學(xué),2012,30(11):1634-1637.
[19] 王儉,孫鐵珩,李培軍,等.基于人工神經(jīng)網(wǎng)絡(luò)的區(qū)域水環(huán)境承載力評(píng)價(jià)模型及其應(yīng)用[J].生態(tài)學(xué)雜志,2007,26(1):139-144.
[20] 趙衛(wèi),劉景雙,蘇偉,等.遼寧省遼河流域水環(huán)境承載力的多目標(biāo)規(guī)劃研究[J].中國(guó)環(huán)境科學(xué),2008,28(1):73-77.
[21] 李如忠.基于指標(biāo)體系的區(qū)域水環(huán)境動(dòng)態(tài)承載力評(píng)價(jià)研究[J].中國(guó)農(nóng)村水利水電,2006(9):42-46.
[22] 張祥娟,李新,李秀霞,等.用層次分析法建立京杭大運(yùn)河蘇州高新區(qū)段水環(huán)境承載力指標(biāo)體系研究[J].環(huán)境保護(hù)與循環(huán)經(jīng)濟(jì),2011(2):42-44.
[23] 王秀蘭,李紅亮.河北省水環(huán)境承載力及污染物總量控制方案研究[J].海河水利,2004(4):31-33.
[24] 趙然航,曹升樂(lè),高輝國(guó).城市水環(huán)境承載力與可持續(xù)發(fā)展策略研究[J].山東大學(xué)學(xué)報(bào),2005,35(2):90-94.