流域TMDL計(jì)劃中的關(guān)鍵技術(shù)參考模板

1、流域TMDL計(jì)劃中的關(guān)鍵技術(shù)徐宗學(xué) 徐華山 吳曉猛摘要: 闡述了流域TMDL 計(jì)劃中的關(guān)鍵技術(shù)問題，包括流域水環(huán)境特征識(shí)別、水污染現(xiàn)狀分析、污染負(fù)荷估算、水文模型和點(diǎn)源污染物削減分配模型構(gòu)建、非點(diǎn)源污染削減措施制定、點(diǎn)源和非點(diǎn)源削減措施實(shí)施效果評(píng)估等，分析了這些技術(shù)的應(yīng)用情況，指出了我國(guó)制定和實(shí)施流域TMDL 計(jì)劃時(shí)應(yīng)注意的相關(guān)問題。關(guān)鍵詞: 最大日負(fù)荷總量; 流域; 污染負(fù)荷分配; 水文水質(zhì)模型; 關(guān)鍵技術(shù); 綜述中圖分類號(hào): X321 文獻(xiàn)標(biāo)志碼: A 文章編號(hào): 1006-7647(2014) 01-0008-06Key technologies of total maximum dai

2、ly load(TMDL) plan for river basinXU Zongxue1，2，XU Huashan1，2，3，WUXiaomeng3Abstract: The key technologies of the total maximum daily load(TMDL) plan for river basin are discussedThose key technologies include feature recognition of water environment，analysis of water pollution status，estimation of p

3、ollution load，construction of hydrological model and model of point source pollution reduction and allocation，formulation of nonpoint source pollution reduction measures，and evaluation of the implementation effects of point source and non-point source pollution reduction measures，the applications of

4、 those key technologies were ananlyzed，and the relevant problems，which should be noted in the formulation and implementation of the TMDL plan for river basin are pointed outKey words: total maximum daily load; river basin; pollution load allocation; hydrological and water quality model; key technolo

5、gy; review最大日負(fù)荷總量(total maximum daily load，TMDL) 計(jì)劃是國(guó)際上最具代表性的流域總量控制技術(shù)體系1。1985 年，美國(guó)在修訂清潔水法時(shí)提出了TMDL 計(jì)劃，制定并頒布了TMDL 實(shí)施導(dǎo)則，1997 年制定了TMDL 技術(shù)指南，2008 年出版了流域TMDL 開發(fā)手冊(cè)草稿，明確了制定流域TMDL 計(jì)劃的步驟: 分析流域特征。收集并分析流域水體數(shù)據(jù)，描述流域狀況和受損狀態(tài)，獲取水體基本信息和流域受損因素，確定河流水質(zhì)標(biāo)準(zhǔn)及TMDL 目標(biāo)，識(shí)別潛在污染源。計(jì)算污染負(fù)荷容量。建立污染負(fù)荷和水質(zhì)關(guān)聯(lián)模型，評(píng)估現(xiàn)有污染負(fù)荷，計(jì)算最大容許污染負(fù)荷總量。確定污染負(fù)

6、荷分配。兼顧地理特征、時(shí)間序列和污染源，選擇最合適、最可行的污染負(fù)荷分配方案。提交TMDL 報(bào)告。起草TMDL報(bào)告，準(zhǔn)備TMDL 文檔，編制TMDL 報(bào)告。流域TMDL 計(jì)劃的關(guān)鍵技術(shù)主要包括流域地表水水質(zhì)受損水體和潛在污染源識(shí)別，選擇與構(gòu)建污染負(fù)荷和水質(zhì)關(guān)聯(lián)技術(shù)模型，計(jì)算最大容許污染負(fù)荷總量，確定點(diǎn)源污染負(fù)荷和非點(diǎn)源污染負(fù)荷分配方案，構(gòu)建公平、合理、可行的點(diǎn)源污染負(fù)荷削減分配模型，優(yōu)化非點(diǎn)源污染控制措施。目前，美國(guó)很多州都對(duì)其行政區(qū)域內(nèi)的水質(zhì)受限水體制定并實(shí)施了TMDL 計(jì)劃，這對(duì)受限水體的水質(zhì)改善發(fā)揮了明顯作用2。我國(guó)關(guān)于流域TMDL 計(jì)劃的研究起步較晚，加上我國(guó)國(guó)情與歐美國(guó)家差別較大，在

7、制定實(shí)施流域TMDL 計(jì)劃時(shí)上述關(guān)鍵技術(shù)也有別于歐美國(guó)家。2 / 141 水質(zhì)受損水體識(shí)別流域水質(zhì)受損水體識(shí)別是制定流域TMDL 計(jì)劃的前提和基礎(chǔ)，一般需確定流域水體水質(zhì)空間差異，識(shí)別污染源并分析污染源與受損河段之間的關(guān)系，同時(shí)通過時(shí)間序列分析評(píng)估受損水體、潛在污染源或其他水體受損條件隨時(shí)間變化的趨勢(shì)，分析水質(zhì)參數(shù)與水文條件之間或者不同水質(zhì)參數(shù)之間的關(guān)系。流域TMDL 計(jì)劃可能涉及不同的子流域、受損河段、特征污染物、可能污染源，這些都會(huì)影響流域水體水質(zhì)。因此，詳細(xì)識(shí)別受損水體對(duì)實(shí)施流域TMDL 計(jì)劃、最大限度利用水資源是非常重要的1。我國(guó)水質(zhì)受損水體識(shí)別多采用水質(zhì)評(píng)價(jià)法，包括單因子指數(shù)法、綜合

8、指數(shù)法、模糊數(shù)學(xué)法等，這些方法對(duì)于分析流域尺度長(zhǎng)時(shí)間序列的水質(zhì)變化略顯不足。多元統(tǒng)計(jì)分析方法包括系統(tǒng)聚類分析法、季節(jié)性Kendall 法、因子-主成分分析法和多元線性回歸分析法，在流域地表水水質(zhì)時(shí)空變化、污染源定性識(shí)別及流域地表水水質(zhì)監(jiān)測(cè)斷面布設(shè)合理性評(píng)估等方面均得到廣泛的應(yīng)用3。國(guó)內(nèi)研究者采用系統(tǒng)聚類分析法對(duì)我國(guó)陜西秦嶺南坡金水河、三峽庫(kù)區(qū)香溪河、漢江流域堵河、漢江上游、湘江4-9等流域地表水水質(zhì)進(jìn)行評(píng)估，分析其時(shí)空變化規(guī)律; 國(guó)外學(xué)者也采用該法分析了西班牙皮蘇埃加河、馬來群島的河口、日本富士河、土耳其柏瑞馬茲河、尼泊爾巴格馬提河、阿根廷蘇琪亞河、韓國(guó)洛東江以及印度的恒河、默哈納迪河、髙緹河

9、、婆羅門河等流域地表水水質(zhì)時(shí)空演變特征10-20。季節(jié)性Kendall 法常用于流域地表水水質(zhì)非參數(shù)趨勢(shì)分析21-28; 因子-主成分分析法可對(duì)河流污染源進(jìn)行定性識(shí)別11，14，17-18，29-33; 多元線性回歸分析法常與主成分分析法一起用于計(jì)算不同污染源對(duì)水體不同指標(biāo)的貢獻(xiàn)率5，17，29，31，34-35。2 水文水質(zhì)模型選擇水文水質(zhì)模型是分析環(huán)境問題的重要工具，是輔助決策支持系統(tǒng)的有力支撐，協(xié)助管理者更好地理解流域系統(tǒng)的動(dòng)態(tài)變化，比如人為污染源和土地管理活動(dòng)污染源如何輸入并影響受納水體以及他們輸入受納水體的條件，如何通過改變這些條件來改變污染物輸入負(fù)荷等，對(duì)所有這些復(fù)雜響應(yīng)關(guān)系的描述

10、主要依靠水文水質(zhì)模型來實(shí)現(xiàn)。TMDL 計(jì)劃的目標(biāo)形式類似于我國(guó)流域總量控制，其核心思想是對(duì)污染負(fù)荷與水質(zhì)的動(dòng)態(tài)響應(yīng)關(guān)系進(jìn)行定量描述與表征，進(jìn)行目標(biāo)可控的管理情景分析，實(shí)現(xiàn)成本效益最大化。水文水質(zhì)模型是制定實(shí)施TMDL 計(jì)劃的關(guān)鍵技術(shù)之一，在核定水環(huán)境容量、估算污染負(fù)荷輸入、評(píng)估受納水體水質(zhì)、確定污染負(fù)荷分配等方面都得到廣泛應(yīng)用。TMDL 計(jì)劃要求模型必須具備量化污染負(fù)荷變化的能力，因此具有物理機(jī)制的模型最合適。如何正確合理地選擇和使用模型是整個(gè)TMDL 計(jì)劃實(shí)施成敗的關(guān)鍵。Depinto等36從哲學(xué)角度研究了選擇模型的基本原則，指出不可能存在一種模型對(duì)所有的TMDL 計(jì)劃都適合，地表水環(huán)境管理

11、的相關(guān)議題使期望的模型復(fù)雜化，但現(xiàn)實(shí)中可利用的數(shù)據(jù)卻又限制了模型的復(fù)雜化，因此，在實(shí)際研究中只能選擇符合管理目標(biāo)的簡(jiǎn)易模型。選擇模型時(shí)要考慮的因素如下: 管理目標(biāo)，包括關(guān)注的污染物、合適的時(shí)空尺度、污染負(fù)荷類型和控制措施、描述水平的精度和準(zhǔn)確度、利益相關(guān)者的接受水平、與其他流域規(guī)劃的相容性等; 特定研究區(qū)特征，包括流域或受納水體尺度和類型、污染物輸移的復(fù)雜性、重要水文過程; 可利用資源，包括數(shù)據(jù)、時(shí)間、預(yù)算、計(jì)算機(jī)需求、工作人員技能等。適合TMDL 計(jì)劃的模型一般具備以下特征:基于流域尺度; 可以連續(xù)模擬; 考慮點(diǎn)源和非點(diǎn)源污染;考慮地表和地下水的貢獻(xiàn); 考慮污染物的溶解態(tài)和顆粒態(tài); 具備表征

12、流域和陸地土地利用特征對(duì)污染物的影響能力。美國(guó)國(guó)家環(huán)保總局(USEPA) 考慮到模型對(duì)于TMDL 計(jì)劃的重要性，兼顧州與州之間水質(zhì)標(biāo)準(zhǔn)的差異、數(shù)據(jù)的可利用性以及使用者的偏好等因素，提供了幾十種備選模型，但是它們當(dāng)中沒有一個(gè)能同時(shí)適應(yīng)所有的流域。應(yīng)用于TMDL 計(jì)劃的模型較多，如AGNPS /AnnAGNPS、GWLF、HSPF、LSPC(HSPF 的升級(jí)版本) 、SWAT和SWMM等。Borah等37將TMDL 計(jì)劃中使用的模型分為負(fù)荷模型、受納水體模型及流域模型。GWLF、PLOAD等為負(fù)荷模型，用來估算從流域出口進(jìn)入到受納水體的沉積物或化學(xué)物質(zhì)的污染負(fù)荷。AQUATOX、BATHTUB、C

13、E-QUAL-W2和QUAL2E等為受納水體模型，用來分析受納水體(如溪流、水壩蓄水、湖泊和河口) 接受污染負(fù)荷后的水量和水質(zhì)變化。AGNPS、AnnAGNPS、HSPF、LSPC 和SWAT等為流域模型，包含負(fù)荷模型和受納水體模型。流域模型又可以進(jìn)一步劃分為兩類，一類是長(zhǎng)系列連續(xù)模擬模型，另外一類是暴雨事件模型。AnnAGNPS、HSPF、LSPC 和SWAT等是連續(xù)模擬模型，用來分析長(zhǎng)序列氣候變化、水文變化、管理措施變化對(duì)流域水量和水質(zhì)的影響。AGNPS 是暴雨事件模型，常用來預(yù)測(cè)單場(chǎng)降雨事件(尤其是極端暴雨事件) 過程中和降雨后洪水、大量泥沙以及營(yíng)養(yǎng)元素對(duì)水質(zhì)和水量的影響，并分析管理措施

14、的效果。暴雨事件模型常被用來分析評(píng)估TMDL 計(jì)劃中最佳管理措施(bestmanagement practices，BMPs) 的實(shí)施效果。另外還有一些模型也具備連續(xù)模擬和單場(chǎng)暴雨事件模擬的能力，如CASC2D/GSSHA、MIKE SHE等。美國(guó)環(huán)境保護(hù)局38總結(jié)了上述這些模型的特征、優(yōu)勢(shì)、局限性、應(yīng)用歷史、時(shí)空尺度、模擬能力、假設(shè)假定及支持的區(qū)域等，并對(duì)模型的界面能力進(jìn)行了詳細(xì)的描述和總結(jié)。TMDL 計(jì)劃另外一個(gè)關(guān)鍵技術(shù)是非點(diǎn)源污染負(fù)荷的定量化，流域農(nóng)業(yè)非點(diǎn)源污染控制技術(shù)多借助分布式水文模型對(duì)流域非點(diǎn)源污染負(fù)荷進(jìn)行估算，然后設(shè)置不同的最佳管理措施并對(duì)比其控制效果、環(huán)境效益和所需費(fèi)用，為管理

15、部門決策提供科學(xué)依據(jù)和技術(shù)支撐。SWAT 模型是USEPA 實(shí)施TMDL計(jì)劃的首選模型，被廣泛應(yīng)用于流域非點(diǎn)源污染關(guān)鍵區(qū)識(shí)別、非點(diǎn)源污染負(fù)荷定量化和BMPs 實(shí)施效果模擬評(píng)估等方面。3 點(diǎn)源污染負(fù)荷分配制定和實(shí)施流域TMDL 計(jì)劃的核心問題是如何科學(xué)合理地分配流域水體污染物削減量，即根據(jù)排污地點(diǎn)、數(shù)量和方式，結(jié)合污染源排污總量削減技術(shù)和優(yōu)先削減污染物順序、經(jīng)濟(jì)可行性等因素，在各控制區(qū)域分配環(huán)境資源，其實(shí)質(zhì)是確定各排污方利用環(huán)境資源的權(quán)利和削減污染物的義務(wù)，即利益的分配和矛盾的協(xié)調(diào)。污染負(fù)荷分配實(shí)際上是利益最大化和污染最小化的矛盾統(tǒng)一體，涉及水體污染物的削減技術(shù)以及經(jīng)濟(jì)的可行性、負(fù)荷分配的公平性

16、，污染負(fù)荷分配合理與否是TMDL 計(jì)劃能否成功實(shí)施的關(guān)鍵39。關(guān)于污染物總量的削減與分配研究在國(guó)外開展比較早，美國(guó)是最早實(shí)施污染物總量控制的國(guó)家，研究制定了污染負(fù)荷分配削減技術(shù)指南40; onald等41在1985 年比較評(píng)價(jià)了8 種比較流行的污染負(fù)荷分析方法; 傳統(tǒng)的TMDL 計(jì)劃在分配容許污染負(fù)荷時(shí)不考慮經(jīng)濟(jì)成本，而Zaidi等42研究了在分配污染負(fù)荷時(shí)如何將經(jīng)濟(jì)分析納入其中。我國(guó)對(duì)水體污染物分配的研究主要以水環(huán)境容量或目標(biāo)總量控制為基礎(chǔ)，并考慮經(jīng)濟(jì)優(yōu)化和公平合理的原則40，基尼系數(shù)法是我國(guó)最常用的一種分配方法?，F(xiàn)行的基尼系數(shù)法分為單指標(biāo)基尼系數(shù)法和綜合指標(biāo)基尼系數(shù)法，這些分配方法均存在一

17、定的不足，如單項(xiàng)社會(huì)經(jīng)濟(jì)指標(biāo)的選擇不能充分反映各區(qū)域的社會(huì)、經(jīng)濟(jì)及資源狀況; 約束條件的設(shè)置過于苛刻，不能體現(xiàn)各單項(xiàng)指標(biāo)重要性差異等。因此，在優(yōu)選相關(guān)指標(biāo)的基礎(chǔ)上，歸一化處理相關(guān)數(shù)據(jù)，引入權(quán)重系數(shù)，構(gòu)建多維水環(huán)境基尼系數(shù)分配模型，并設(shè)置不同的約束條件，優(yōu)化求解分配方案是今后一段時(shí)期點(diǎn)源污染負(fù)荷削減與分配的重要工作之一。4 非點(diǎn)源污染負(fù)荷BMPs優(yōu)化非點(diǎn)源污染的控制可以從兩個(gè)方面著手: 一是控制污染源，即從源頭上控制污染物進(jìn)入水體; 二是控制污染物遷移轉(zhuǎn)化途徑。由于認(rèn)識(shí)問題的深淺不同，過去人們往往僅依靠一種措施來控制復(fù)雜的非點(diǎn)源污染問題，隨著對(duì)非點(diǎn)源污染問題的深入研究，逐漸認(rèn)識(shí)到依靠單一措施控制

18、非點(diǎn)源污染的效果并不理想，需要根據(jù)特定流域的水文氣候因素和地形要素制定適合該流域的管理措施43-44。由于非點(diǎn)源污染具有分布廣、隨機(jī)性強(qiáng)、潛伏期長(zhǎng)、難以估算等特點(diǎn)，致使非點(diǎn)源污染的控制難度大、成本高45，各國(guó)在實(shí)踐中提出過多種調(diào)控方法，涉及多個(gè)領(lǐng)域，歸納起來有以下4 種系統(tǒng)管理措施:農(nóng)田系統(tǒng)管理措施，包括養(yǎng)分管理和家畜管理;生態(tài)系統(tǒng)管理措施，包括河流、湖泊、海灣的地表水管理; 農(nóng)作物系統(tǒng)管理措施，包括農(nóng)作物管理和土壤管理; 水文系統(tǒng)管理措施，包括水資源管理、土地使用管理和景觀管理。高效、經(jīng)濟(jì)、科學(xué)的最佳管理措施BMPs 在非點(diǎn)源污染控制中得到廣泛應(yīng)用46-47。在歐洲、美國(guó)和其他一些國(guó)家，BM

19、Ps 優(yōu)化常常借助于SWAT 模型，許多研究專注于BMPs 的環(huán)境效益評(píng)價(jià)48-52，也結(jié)合經(jīng)濟(jì)學(xué)理論進(jìn)行BMPs 的費(fèi)用-效益分析研究53-56。Qiu等57利用SWAT 模型對(duì)農(nóng)業(yè)流域的河岸緩沖帶進(jìn)行了經(jīng)濟(jì)評(píng)價(jià)，發(fā)現(xiàn)研究區(qū)的河岸緩沖帶能有效消除農(nóng)田帶來的農(nóng)藥，強(qiáng)烈建議農(nóng)民發(fā)展和保留臨近水體的河岸緩沖帶; Panagopoulos等58在2007 年提出了BMPs 費(fèi)用-效益分析計(jì)算方法，并在2011 年采用SWAT 模型分析了不同的BMPs 對(duì)流域非點(diǎn)源負(fù)荷削減的效果，進(jìn)行了效益對(duì)比分析; abotyagov等59指出了不同種類的BMPs費(fèi)用，為研究者選擇經(jīng)濟(jì)合理的BMPs 提供了參考;

20、Maringanti等60發(fā)展了一種多目標(biāo)優(yōu)化方法來選擇和實(shí)施BMPs，成功地考慮了經(jīng)濟(jì)和生態(tài)效應(yīng)。采用SWAT 模型輸出基準(zhǔn)情景下HU 水平的污染負(fù)荷，同時(shí)將SWAT 模型不同的運(yùn)行結(jié)果、經(jīng)濟(jì)數(shù)據(jù)和等位基因序列輸入BMPs 優(yōu)化模型。在優(yōu)化過程中，首先在特定土地利用類型給出的等位基因序列中搜尋特定的管理措施，然后估算所選取的最佳管理措施費(fèi)用，最后在流域尺度上加權(quán)平均HU 水平上的污染負(fù)荷和總費(fèi)用。優(yōu)化過程是使流域總污染負(fù)荷和總費(fèi)用這兩個(gè)目標(biāo)函數(shù)最小，結(jié)合SWAT 模型和優(yōu)化模型采用多目標(biāo)遺傳算法求解。5 結(jié)語與展望我國(guó)流域TMDL 計(jì)劃的相關(guān)研究起步較晚，已有成果多停留在對(duì)流域TMDL 計(jì)劃

21、關(guān)鍵技術(shù)的介紹以及對(duì)技術(shù)體系的理論探討上，尚未開展具體的案例研究。同時(shí)，我國(guó)河流受到人類活動(dòng)強(qiáng)烈干擾，TMDL 技術(shù)體系中的點(diǎn)源污染負(fù)荷削減與分配、非點(diǎn)源污染控制措施優(yōu)化、水文模型構(gòu)建等關(guān)鍵技術(shù)都有別于歐美國(guó)家，使得該技術(shù)體系在我國(guó)的應(yīng)用受到限制。因此，我國(guó)流域TMDL 計(jì)劃的制定和實(shí)施要根據(jù)我國(guó)國(guó)情區(qū)別對(duì)待，概括起來需要注意以下問題:a流域水質(zhì)時(shí)空演變分析需要借助流域管理機(jī)構(gòu)多年的水質(zhì)監(jiān)測(cè)數(shù)據(jù)，采用多元統(tǒng)計(jì)方法分析水質(zhì)在流域尺度上的空間差異和多年變化趨勢(shì)，結(jié)合污染源調(diào)查，分析引起水體水質(zhì)受損的原因，為制定流域TMDL 計(jì)劃提供基礎(chǔ)數(shù)據(jù)。b根據(jù)流域的不同特點(diǎn)選擇并修訂模型，如在選擇SWAT 模

22、型作為流域TMDL 計(jì)劃的模型時(shí)，對(duì)于我國(guó)部分流域閘壩較多、取用水頻繁等實(shí)際情況，對(duì)相應(yīng)模塊進(jìn)行修改是可行的，也是必要的。c現(xiàn)階段我國(guó)點(diǎn)源污染負(fù)荷依然是造成流域水體水質(zhì)受損的主要因素，基于我國(guó)國(guó)情建立一套科學(xué)、公平、各排污方均能接受且能順利推行實(shí)施的分配體系，是我國(guó)流域TMDL 計(jì)劃制定實(shí)施的重點(diǎn)，也是難點(diǎn)所在。d根據(jù)流域特點(diǎn)并配合必要的經(jīng)濟(jì)效益分析，制定適合流域非點(diǎn)源污染負(fù)荷的管理措施并付諸實(shí)施，是將來制定和實(shí)施流域TMDL 計(jì)劃必須考慮的關(guān)鍵問題之一。參考文獻(xiàn):1張楠基于不確定性的流域TMDL 及其安全余量研究D北京: 北京師范大學(xué)，20092柯強(qiáng)，趙靜，王少平，等最大日負(fù)荷總量(TMDL

23、) 技術(shù)在農(nóng)業(yè)面源污染控制與管理中的應(yīng)用與發(fā)展趨勢(shì)J生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2009，25(1) : 85-91(KE Qiang，ZHAO Jing，WANG Shaoping，et alApplication of total maximum daily load(TMDL) in control of agricultural nonpoint source pollution and its developmental trendJJournal of Ecology and ural Environment，2009，25(1) : 85-91(in Chinese) )3周豐，郭懷成，

24、黃凱，等基于多元統(tǒng)計(jì)方法的河流水質(zhì)空間分析J水科學(xué)進(jìn)展，2007，18(4 ) : 544-551(ZHOU Feng，GUO Huaicheng，HUANG Kai，et alMultivariate statistical technique for spatial variation in river water qualityJAdvances in Water Science，2007，18(4) ，544-551(in Chinese) )4卜紅梅，劉文治，張全發(fā)多元統(tǒng)計(jì)方法在金水河水質(zhì)時(shí)空變化分析中的應(yīng)用J資源科學(xué)，2009，31(3) :429-434(BU Hongmei，L

25、IU Wenzhi，ZHANG Quanfa，et alApplication of multiple statistical analysis to spatialtemporal variations of water quality of the Jinshui iver Jesources Science，2009，31(3 ) ，429-434(in Chinese) )5 BU Hongmei，TAN Xiang，LI Siyue，et alWater quality assessment of the Jinshuiiver(China ) using multivariate

26、statistical techniques JEnvironmental Earth Sciences，2010，60(8) : 1631-16396葉麟，黎道豐，唐濤，等香溪河水質(zhì)空間分布特性研究J應(yīng)用生態(tài)學(xué)報(bào)，2003，14(11 ) : 1959-1962(YE Lin，LI Daofeng，TANG Tao，et alSpatial distribution of water quality in Xiangxiiver of China JChinese Journal of Applied Ecology，2003，14(11 ) : 1959-1962(in Chinese)

27、)7顧勝，李思悅，張全發(fā)漢江堵河流域地表水質(zhì)時(shí)空變化特征J長(zhǎng)江流域資源與環(huán)境，2009，18(1) : 41-46 (GU Sheng，LI Siyue，ZHANG QuanfaSpatital-temporal distribution of surface water quality in the Du iver Basin Jesources and Environment in the Yangtze Basin，2009，18(1) ，41-46(in Chinese) )8LI Siyue，ZHANG QuanfaSpatial characterization of dissol

28、ved trace elements and heavy metals in the upper Han iver(China ) using multivariate statistical techniquesJJournal of Hazardous Materials，2010，176 (1 /2 /3) : 579-5889ZHANG Qi，LI Zhongwu，ZENG Guangming，et alAssessment of surface water quality using multivariate statistical techniques in red soil hi

29、lly region: a case study of Xiangjiang watershed of China JEnvironmental Monitoring and Assessment，2009，152: 123-13110VEGA M，PADO ，BAADO E，et alAssessment of seasonal and polluting effects on the quality of river water by exploratory data analysisJWater esearch，1998，32 (12) : 3581-359211ALKAKHI A F

30、M，AHMAD A，EASA A MAssessment of surface water quality of selected estuaries of Malaysia: multivariate statistical techniquesJEnvironmentalist，2009，29(3) : 255-26212SHESTHA S，KAZAMA FAssessment of surface water quality using multivariate statistical techniques: a case study of the Fuji iver Basin of

31、JapanJEnvironmental Modelling Software，2007，22(4) : 464-47513 VAOL M，SEN BAssessment of surface water quality using multivariate statistical techniques: a case study of Behrimaz Stream of TurkeyJEnvironmental Monitoring and Assessment，2009，159: 543-55314KANNEL P ，LEE S，KANEL S ，et alChemometric appl

32、ication in classification and assessment of monitoring locations of an urban river systemJAnalytica Chimica Acta，2007，582(2) : 390-39915 SINGH K P，MALIK A，MOHAN D，et alMultivariate statistical techniques for the evaluation of spatial and temporal variations in water quality of Gomti iver (India) : a

33、 case studyJWater esearch，2004，38 (18) : 3980-399216 SUNDAAY S K，PANDA U C，NAYAK B B，et alMultivariate statistical techniques for the evaluation of spatial and temporal variations in water quality of the Mahanadi river-estuarine system(India ) : a case study JEnvironmental Geochemistry and Health，20

34、06，28 (4) : 317-33017SINGH K P，MALIK A，SINHA SWater quality assessment and apportionment of pollution sources of Gomti iver(India ) using multivariate statistical techniques: a case studyJAnalytica Chimica Acta，2005，538(1 /2) : 355-37418 SUNDAAY S KApplication of multivariate statistical techniques

35、in hydrogeochmical studies: a case study of Brahmani-Koel iver(India) JEnviron Monit Assess，2010，164: 297-31019ALBETO W D，DEL PILA D M，VALEIA A M，et alPattern recognition techniques for the evaluation of spatial and temporal variations in water quality: a case study of Suquia iver Basin(Cordoba-Arge

36、ntina) JWater esearch，2001，35(12) : 2881-289420HAN S，KIM E，KIM SThe water quality management in the Nakdong iver Watershed using multivariate statistical techniquesJKsce Journal of Civil Engineering，2009，13(2) : 97-10521ZHANG J，BUKE P，IICANIN N，et alLong-term water quality trends in the Lake Okeecho

37、bee Watershed of FloridaJCritical eviews in Environmental Science and Technology，2011，41(Sup1) : 548-57522LEE H W，BHANG K J，PAK S SEffective visualization for the spatiotemporal trend analysis of the water quality in the Nakdong iver of KoreaJEcological Informatics，2010，5(4) : 281-29223YIN X，HU G H，

38、SOC I CTrend analysis of surface water quality in Hunan ProvinceMLos Alamitos: Ieee Computer Society，200924HIONO Y，WATANABE I，NONAKA KTrends in water quality around an intensive tea-growing area in Shizuoka，JapanJSoil Science and Plant Nutrition，2009，55(6) : 783-79225 SUN Dong，LI Junfeng，SUN FeifeiS

39、easonal Kendall test and its application in water quality trend analysisM Liverpool: World Acad Union-World Acad Press，200826JOHNSON H O，GUPTA S C，VECCHIA A V，et al Assessment of water quality trends in the Minnesota iver using Non-Parametric and parametric methodsJJournal of Environmental Quality，2

40、009，38(3) : 1018-103027BOUZA-DEANO ，TENEO-ODIGUEZ M，F(xiàn)ENANDEAESPINOSA A JTrend study and assessment of surface water quality in the Ebro iver(Spain) JJournal of Hydrology，2008，361(3 /4) : 227-23928CHANG HSpatial analysis of water quality trends in the Han iver Basin of South KoreaJWater esearch，2008，

41、42(13) : 3285-330429SIMEONOV V，STATIS J A，SAMAA C，et alAssessment of the surface water quality in Northern Greece JWater esearch，2003，37(17) : 4119-412430KAZI T G，AAIN M B，JAMALI M K，et alAssessment of water quality of polluted lake using multivariate statistical techniques: a case studyJEcotoxicolo

42、gy and Environmental Safety，2009，72(2) : 301-30931PEKEY H，KAAKAS D，BAKOGLU MSource apportionment of trace metals in surface waters of a polluted stream using multivariate statistical analysesJ Marine Pollution Bulletin，2004，49(9 /10) : 809-81832OMO-IABO O O，OLOBANIYI S B，ODUYEMLI K，et alSurface and

43、groundwater water quality assessment using multivariate analytical methods: a case study of the Western Niger Delta of NigeriaJPhysics and Chemistry of the Earth，2008，33: 666-67333AZZELLINO A，SALVETTI ，VISMAA ，et alCombined use of the EPA-Qual2E simulation model and factor analysis to assess the sou

44、rce apportionment of point and non point loads of nutrients to surface watersJScience of the Total Environment，2006，371(1 /2 /3 ) : 214-22234ZHOU Feng，HUANG G H，GUO Huaicheng，et alSpatiotemporal patterns and source apportionment of coastal water pollution in eastern Hong Kong JWater esearch，2007，41(

45、15) : 3429-343935SU Sjiliang，ZHI Junjun，LOU Liping，et alSpatiotemporal patterns and source apportionment of pollution in Qiantang iver(China) using neural-based modeling and multivariate statistical techniques JPhysics and Chemistry of the Earth，2011，36(9 /10 /11) : 379-38636DEPINTO J V，F(xiàn)EEDMAN P L，

46、DILKS D M，et alModels quantify the total maximum daily load processJJournal of Environmental Engineering-Asce，2004，130 (6) : 703-71337BOAH D K，YAGOW G，SALEH A，et alSediment and nutrient modeling for TMDL development and implementationJTransactions of the Asabe，2006，49 (4) : 967-98638美國(guó)環(huán)境保護(hù)局美國(guó)TMDL 計(jì)劃

47、管理模型實(shí)施實(shí)踐M王東，趙越，徐敏，等，譯北京: 中國(guó)環(huán)境科學(xué)出版社，201239劉贛明最大日負(fù)荷總量(TMDL) 模式下的污染負(fù)荷分配D廣州: 中山大學(xué)，200540張杏杏，彭曉春，賀濤，等水污染物總量分配方法研究進(jìn)展J環(huán)境科學(xué)與管理，2011，36(4) : 3-7(ZHANG Xingxing，PEN Xiaochun，HE Tao，et alAdvances in the research of total pollutant load allocat ion for surface water JEnvironmental Science and Management，2011，36

48、 (4) : 3-7(in Chinese) )41ONALD C，IENE KAn evaluation of eight waste-load allocation methodsJWateresources Buletin，1985，21 (5) : 833-83942 ZAIDI A Z，DEMONSABET S M，El-FAHAN How to include economic analysis in TMDL allocationJ Journal of Water esources Planning and Management-Asce，2008，134(3) : 214-2

49、2343PANAGOPOULOS Y，MAKOPOULOS C，MIMIKOU Meducing surface water pollution through the assessment of the cost-effectiveness of BMPs at different spatial scales JJournal of Environmental Management，2011，92 (10) : 2823-283544AO N S，EASTON Z M，SCHNEIDEMAN E M，et al Modeling watershed-scale effectiveness

50、of agricultural best management practices to reduce phosphorus loadingJ Journal of Environmental Management，2009，90(3 ) : 1385-139545周慧平，高超，朱曉東關(guān)鍵源區(qū)識(shí)別: 農(nóng)業(yè)非點(diǎn)源污染控制的方法J生態(tài)學(xué)報(bào)，2005，25(12) : 3368-3374(ZHOU Huiping，GAO Chao，ZHU XiaodongIdentification of critical source areas: an efficient way for agricultura

51、l nonpoint source pollution controlJActa Ecologica Sinica，2005，25(12) : 3368-3374(in Chinese) )46王曉燕非點(diǎn)源污染與管理M北京: 海洋出版社，200347MCKISSOCK G，JEFFEIES C，DACY B JAn assessment of drainage best management practices in ScotlandJJ Chart Inst Water Environ Manage，1999，13 (1) : 47-5148AABI M，GOVINDAAJU S，HANTU

52、SH M M，et alole of watershed subdivision on modeling the effectiveness of best management practices with swatJ Journal of the American Water esources Association，2006，42(2) : 513-52849BALUND I，KIKKALA T，MALVE O，et alAssessing SWAT model performance in the evaluation of management actions for the imp

53、lementation of the water framework directive in a finnish catchment JEnvironmental Modelling Software，2007，22(5) : 719-72450SAHU M，GU Modeling the effects of riparian buffer zone and contour strips on stream water quality JEcological Engineering，2009，35(8) : 1167-117751SANTHI C，SINIVASAN ，ANOLD J，et alA modeling approach to evaluate the impacts of water quality management plans implemented in a watershed in Texas JEnvironmental Modelling Software，2006，21(8) : 1141-115752WHITE M J，ANOLD J GDevelopment of