曹慧群，趙 鑫

(1.長(zhǎng)江科學(xué)院 a.流域水環(huán)境研究所;b.流域水資源與生態(tài)環(huán)境科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室，武漢 430010;2.三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心，湖北宜昌 443002)

1 研究背景

在氣候變化與人類(lèi)活動(dòng)影響下，流域水沙過(guò)程及以其為載體的水環(huán)境演變過(guò)程均發(fā)生了很大變化。根據(jù)國(guó)務(wù)院最新出臺(tái)的《水污染防治行動(dòng)計(jì)劃》(國(guó)發(fā)〔2015〕17號(hào))，水環(huán)境保護(hù)是國(guó)家水安全保障的重中之重。為更有效地開(kāi)展流域水環(huán)境保護(hù)工作，有必要掌握新形勢(shì)下流域水環(huán)境演變過(guò)程。現(xiàn)有研究手段主要為野外監(jiān)測(cè)、物理模型試驗(yàn)和數(shù)值模擬。其中水環(huán)境數(shù)值模擬借助計(jì)算機(jī)技術(shù)，輸入電子化的水域地形條件和邊界條件，求解反映水動(dòng)力和水環(huán)境演變過(guò)程的理論方程，得到水動(dòng)力、泥沙及水環(huán)境要素的時(shí)空分布規(guī)律，可非常逼真地復(fù)演水環(huán)境系統(tǒng)中各要素的變化過(guò)程及其相互關(guān)系。具有不可替代的優(yōu)勢(shì):①通過(guò)大量方案模擬，可快速掌握水環(huán)境內(nèi)部因子的變化規(guī)律;②可模擬水體全范圍內(nèi)污染物濃度的時(shí)空變化過(guò)程，對(duì)水環(huán)境變化進(jìn)行定量描述［1－3］，同時(shí)可對(duì)未來(lái)變化進(jìn)行預(yù)測(cè);③耗資小、不受人力、物力的限制，研究周期短，可有效提高水環(huán)境管理工作的效率。因此，水環(huán)境數(shù)值模擬是水環(huán)境規(guī)劃、管理與保護(hù)工作中的重要研究手段。

目前國(guó)內(nèi)外針對(duì)流域水環(huán)境數(shù)值模擬開(kāi)展了大量研究工作，已取得了一定成果。但受到監(jiān)測(cè)數(shù)據(jù)有限、污染物輸移機(jī)制還未完全清晰等的影響，流域水環(huán)境模擬技術(shù)發(fā)展受到一定局限，未來(lái)仍有較大發(fā)展空間。本文總結(jié)了目前水環(huán)境數(shù)值模擬技術(shù)研究進(jìn)展，歸納了水環(huán)境模擬在流域水環(huán)境演變與調(diào)控研究、水利工程對(duì)流域水環(huán)境的影響研究、水污染防治工程設(shè)計(jì)、水資源保護(hù)規(guī)劃與管理等方面的應(yīng)用情況，并提出了流域水環(huán)境數(shù)值模擬技術(shù)的發(fā)展趨勢(shì)。

2 當(dāng)前主要水環(huán)境數(shù)值模擬技術(shù)

水環(huán)境數(shù)學(xué)模型的開(kāi)發(fā)，始于1925年Streeter和Phelps建立的水中溶解氧平衡模型(S-P模型)，此后逐漸由一維單因素河道穩(wěn)態(tài)模型擴(kuò)展至多維多因素動(dòng)態(tài)模型，同時(shí)考慮底泥釋放和大氣沉降的作用，并引入模糊數(shù)學(xué)、人工神經(jīng)網(wǎng)絡(luò)、3S等新技術(shù)［4］。水環(huán)境模擬技術(shù)主要由國(guó)外相關(guān)機(jī)構(gòu)研發(fā)，如美國(guó)環(huán)保局、美國(guó)陸軍工程兵團(tuán)、丹麥水動(dòng)力研究所等。

根據(jù)模擬對(duì)象的不同，將水環(huán)境模型分為3類(lèi):地表水環(huán)境模型、地下水環(huán)境模型、流域非點(diǎn)源污染模型。

2.1 地表水環(huán)境模型

地表水環(huán)境模型主要模擬水流、泥沙、污染物質(zhì)在河流、湖泊、河口等地表水體中的運(yùn)動(dòng)過(guò)程，代表性的模型有 EFDC，WASP，MIKE，Delft-3D，SMS，CE-QUAL-R1/W2/ICM 等［5－10］，具 體 見(jiàn) 表 1。 其 中EFDC模型是目前應(yīng)用最廣的三維模型，在三維水動(dòng)力計(jì)算方面具有突出優(yōu)勢(shì)，能夠較好適應(yīng)長(zhǎng)時(shí)間干濕交替情景，模型前期率定和后期可視化處理功能較強(qiáng);WASP模型在數(shù)值計(jì)算、參數(shù)率定上具有優(yōu)勢(shì)，但前期概化及網(wǎng)格化較繁瑣;MIKE系列模型可模擬多種水質(zhì)參數(shù)，適用于復(fù)雜水域，但模型源代碼不公開(kāi)，用戶(hù)拓展難度很大;Delft-3D模型嵌有水質(zhì)和生態(tài)過(guò)程庫(kù)，可拓展性較強(qiáng);SMS模型的優(yōu)點(diǎn)是具有強(qiáng)大的網(wǎng)格創(chuàng)建能力;CE-QUAL-R1/W2主要適用于湖庫(kù)的水質(zhì)模擬［11］。

表1 代表性水環(huán)境模型匯總Table 1 Representative models of water environment

2.2 地下水環(huán)境模型

地下水環(huán)境模型主要模擬地下水流與溶質(zhì)的運(yùn)移過(guò)程，代表性的模型有Visual MODFLOW，GMS等［12－13］，具體見(jiàn)表1。GMS為可視化三維地下水模擬系統(tǒng)，集成了 MODFLOW，MODPATH，MT3D，F(xiàn)EMWATER等近10個(gè)地下水模型，可進(jìn)行地下水的水流模擬、溶質(zhì)運(yùn)移模擬和反應(yīng)運(yùn)移模擬等;Visual MODFLOW同樣以MODFLOW模型為核心，特點(diǎn)是界面友好，可視化性能較強(qiáng)，使用方便。

2.3 流域非點(diǎn)源污染模型

流域非點(diǎn)源污染模型重點(diǎn)模擬陸面徑流產(chǎn)生以及污染物質(zhì)輸移過(guò)程，代表性的模型有 SWAT，HSPF，MIKE SHE，SWMM 等［14－17］，具體見(jiàn)表 1。SWAT模型是目前應(yīng)用最廣泛的流域非點(diǎn)源污染模型之一，適用于土壤類(lèi)型和土地利用方式復(fù)雜的流域，但應(yīng)用至我國(guó)流域需要進(jìn)行代碼轉(zhuǎn)換，工作量較大;HSPF模型以強(qiáng)大的水文模擬為基礎(chǔ)，模擬精度較高，但對(duì)輸入數(shù)據(jù)要求高，空間分辨率相對(duì)較低;MIKE SHE以完全物理分布式水文模型為基礎(chǔ)，通過(guò)多模塊配置可模擬多種要素，同樣對(duì)基礎(chǔ)資料要求高;SWMM模型適用于具有復(fù)雜下墊面城市的暴雨徑流和污染物運(yùn)動(dòng)過(guò)程。

3 流域水環(huán)境數(shù)值模擬技術(shù)應(yīng)用

水環(huán)境模擬通過(guò)對(duì)水體內(nèi)部發(fā)生的復(fù)雜污染過(guò)程定量化描述，識(shí)別主要污染物及其遷移途徑和轉(zhuǎn)化過(guò)程，具有較廣的應(yīng)用空間。目前，水環(huán)境模擬技術(shù)在基礎(chǔ)研究、工程應(yīng)用和規(guī)劃管理等方面均有應(yīng)用，其中基礎(chǔ)研究集中在流域水環(huán)境演變與調(diào)控研究方面，工程應(yīng)用集中在水利工程對(duì)流域水環(huán)境的影響、水污染防治工程設(shè)計(jì)方面，規(guī)劃管理集中在水資源保護(hù)規(guī)劃與管理方面。

3.1 流域水環(huán)境演變與調(diào)控研究

在全球氣候變化與人類(lèi)活動(dòng)影響下，流域水環(huán)境的演變發(fā)生復(fù)雜變化，可通過(guò)建立流域水環(huán)境數(shù)學(xué)模型，對(duì)水環(huán)境歷史演變進(jìn)行反演，通過(guò)定量的數(shù)值模擬結(jié)果探討變化的主要原因，并對(duì)未來(lái)變化趨勢(shì)進(jìn)行預(yù)測(cè)。全球氣候變化的影響研究，一般采用聯(lián)合國(guó)政府間氣候變化專(zhuān)門(mén)委員會(huì)發(fā)布的IPCC評(píng)估報(bào)告成果，將過(guò)去已經(jīng)發(fā)生的氣候變化和未來(lái)設(shè)計(jì)的氣候情景輸入模型，模擬分析流域水文循環(huán)過(guò)程及水量變化，重點(diǎn)分析洪澇與干旱等極端水文事件的發(fā)生頻次［18］，進(jìn)一步分析污染物遷移轉(zhuǎn)化和水體稀釋能力的變化，研究氣候變化對(duì)流域點(diǎn)源與非點(diǎn)源污染事件發(fā)生概率和量級(jí)的影響［19］;對(duì)于大規(guī)模人類(lèi)活動(dòng)的影響研究，可通過(guò)建立流域尺度水文與水環(huán)境數(shù)學(xué)模型，模擬分析土地利用方式改變、沿岸經(jīng)濟(jì)帶城市產(chǎn)業(yè)轉(zhuǎn)移等條件下，流域水環(huán)境的演變過(guò)程。在此基礎(chǔ)上，從流域宏觀層面，提出流域水沙調(diào)控對(duì)策，利用流域水環(huán)境模型模擬對(duì)比，篩選提出最優(yōu)方案，并分析流域水生態(tài)環(huán)境的響應(yīng)。

3.2 水利工程影響分析與水污染防治工程設(shè)計(jì)

重大水利工程(如水庫(kù)工程、水系連通工程、跨流域調(diào)水工程等)的運(yùn)行改變了江河湖庫(kù)的天然水文過(guò)程，對(duì)流域水環(huán)境產(chǎn)生較大的影響，一般可利用地表水環(huán)境數(shù)學(xué)模型開(kāi)展研究。水庫(kù)水環(huán)境數(shù)學(xué)模型一般考慮水動(dòng)力、泥沙和水質(zhì)等要素，入口邊界條件一般采用典型系列年的入庫(kù)水沙和水質(zhì)條件，出口邊界條件采用水庫(kù)調(diào)度的壩前水位過(guò)程。與自然河流相比，水庫(kù)水深變化幅度較大，且壩前水深較大，通常采用三維或垂向二維模型研究。根據(jù)水庫(kù)庫(kù)區(qū)水沙與水環(huán)境時(shí)空變化模擬結(jié)果，可進(jìn)一步開(kāi)展水庫(kù)的水質(zhì)達(dá)標(biāo)分析、水庫(kù)對(duì)水沙及污染物(營(yíng)養(yǎng)物)輸移過(guò)程的影響研究、水庫(kù)生態(tài)調(diào)度研究［20－21］等;河湖水系連通工程通過(guò)港渠等連通工程，將河流與湖泊連通，加快河湖之間的水體流動(dòng)，改善水環(huán)境質(zhì)量。但不同水體的匯合可能會(huì)帶來(lái)一定生態(tài)環(huán)境風(fēng)險(xiǎn)，可通過(guò)建立水動(dòng)力－水環(huán)境－水生態(tài)模型，研究水系連通對(duì)水環(huán)境改善作用的幅度以及產(chǎn)生生態(tài)風(fēng)險(xiǎn)的等級(jí)，并通過(guò)模擬比選，提出最優(yōu)的河湖水系連通方案［22］;此外，跨流域調(diào)水工程對(duì)調(diào)水區(qū)及受水區(qū)的水環(huán)境影響，亦可利用水環(huán)境數(shù)學(xué)模型開(kāi)展研究。

非點(diǎn)源污染控制是當(dāng)前水污染治理的難點(diǎn)，也是生態(tài)清潔小流域建設(shè)的重點(diǎn)。生態(tài)溝渠和人工濕地是非點(diǎn)源污染過(guò)程控制與末端治理的重要工程技術(shù)措施［23］，其工程設(shè)計(jì)需要借助流域非點(diǎn)源污染模型［24］。以小流域集水范圍為對(duì)象，建立流域水文與非點(diǎn)源污染數(shù)學(xué)模型，針對(duì)不同生態(tài)溝渠與人工濕地工程方案，重點(diǎn)考慮介質(zhì)特征、粒徑、水平布置與垂向布置、基質(zhì)與植物的類(lèi)型、植物種植密度與空間布置等，模擬比對(duì)提出最優(yōu)的工程配置與設(shè)計(jì)方案。

3.3 水資源保護(hù)規(guī)劃與管理

水環(huán)境模擬可定量確定污染負(fù)荷，為水功能區(qū)納污能力復(fù)核、水功能區(qū)達(dá)標(biāo)分析、入河排污口布局等水資源保護(hù)規(guī)劃工作提供決策依據(jù)。

水功能區(qū)納污能力是指滿(mǎn)足水功能區(qū)水質(zhì)目標(biāo)要求的污染物最大允許負(fù)荷量。水資源保護(hù)規(guī)劃需要提出現(xiàn)狀及規(guī)劃條件下的納污能力。數(shù)學(xué)模型計(jì)算法是納污能力計(jì)算的最有效方法之一。參照《水域納污能力計(jì)算規(guī)程》，根據(jù)河流污染物的混合程度，可采用零維、一維、二維模型計(jì)算納污能力，其設(shè)計(jì)流量一般采用最枯月均流量。隨著水環(huán)境模擬計(jì)算的發(fā)展，可進(jìn)一步考慮污染物的垂向擴(kuò)散、典型設(shè)計(jì)流量過(guò)程、河道地形變化等重要因素，開(kāi)發(fā)三維水環(huán)境數(shù)學(xué)模型，更高精度模擬計(jì)算河道的納污能力［25］，以滿(mǎn)足水環(huán)境精細(xì)化管理的需求;水功能區(qū)達(dá)標(biāo)分析根據(jù)水功能水質(zhì)目標(biāo)，判斷現(xiàn)狀水質(zhì)是否達(dá)標(biāo)，其中現(xiàn)狀水質(zhì)分布是水功能區(qū)達(dá)標(biāo)分析的基礎(chǔ)，需要利用有限的水質(zhì)監(jiān)測(cè)數(shù)據(jù)，結(jié)合數(shù)學(xué)模型計(jì)算分析得到;入河排污口的設(shè)置應(yīng)不影響水功能區(qū)的達(dá)標(biāo)，論證過(guò)程中可利用水環(huán)境數(shù)學(xué)模型，模擬計(jì)算不同排污口設(shè)置方案條件下，污染物的輸移過(guò)程、污染物濃度等值線分布及發(fā)展過(guò)程，分析污水排放對(duì)水功能區(qū)的影響范圍，從中比選出影響最小的排污口布局方案。

4 水環(huán)境數(shù)值模擬技術(shù)發(fā)展趨勢(shì)

目前水環(huán)境模擬技術(shù)存在以下不足:①模擬的過(guò)程相對(duì)單一，參數(shù)率定及模擬結(jié)果存在較大隨機(jī)性;②模擬因素考慮的不夠全面，模擬結(jié)果與實(shí)際情況有一定差距;③模型專(zhuān)業(yè)性較強(qiáng)，難以推廣應(yīng)用。隨著水環(huán)境監(jiān)測(cè)站網(wǎng)的完善、污染物在不同介質(zhì)中轉(zhuǎn)化機(jī)理與遷移過(guò)程研究的深入、以及計(jì)算機(jī)技術(shù)與地理信息技術(shù)的快速發(fā)展，為水環(huán)境數(shù)值模擬技術(shù)的進(jìn)一步發(fā)展提供可行性。

未來(lái)發(fā)展趨勢(shì)主要集中在流域水環(huán)境全過(guò)程模擬、水環(huán)境多因素耦合模擬、水環(huán)境模擬系統(tǒng)集成3方面。

4.1 流域全過(guò)程模擬

流域指由分水線所包圍的河流集水區(qū)，為水體自然流動(dòng)的邊界，具有自然整體性。對(duì)于流域水環(huán)境演變，水沙、營(yíng)養(yǎng)物與污染物首先從陸面侵蝕，再進(jìn)入河道對(duì)流與擴(kuò)散，并在地表水與地下水之間遷移轉(zhuǎn)化，是一個(gè)連續(xù)的輸移過(guò)程。隨著計(jì)算模擬技術(shù)與地理信息系統(tǒng)技術(shù)的發(fā)展，建立陸面、地表水與地下水水環(huán)境集成模型，客觀描述水環(huán)境全過(guò)程演變成為可能。此外，隨著水環(huán)境監(jiān)測(cè)站網(wǎng)的不斷完善，陸面、地表與地下水中的水文、水質(zhì)、水生生物的統(tǒng)一監(jiān)測(cè)體系逐漸建立，為流域全過(guò)程模擬以及參數(shù)準(zhǔn)確率定提供可行條件。目前研究重點(diǎn)是陸面－地表水環(huán)境耦合模擬、污染物在地表水和地下水中的遷移轉(zhuǎn)化過(guò)程模擬等［26］。研究難點(diǎn)是不同空間尺度之間的轉(zhuǎn)化方法，以及水沙、營(yíng)養(yǎng)物、污染物在不同空間尺度中的輸移機(jī)制及過(guò)程。

4.2 多因素耦合模擬

流域環(huán)境中主要包括水、泥沙、無(wú)機(jī)物和有機(jī)物、水生生物等物質(zhì)，這些物質(zhì)相互作用，存在復(fù)雜的耦合關(guān)系，如河流水體中的氮磷營(yíng)養(yǎng)物質(zhì)和重金屬物質(zhì)以溶于水的溶解態(tài)、吸附于泥沙上的顆粒態(tài)2種形態(tài)輸移;河流中營(yíng)養(yǎng)物質(zhì)能夠被浮游植物迅速吸收，然后通過(guò)浮游植物－浮游動(dòng)物－魚(yú)類(lèi)的食物鏈形式，影響水生生物的生長(zhǎng)和繁殖;河流中的有毒污染物，如環(huán)境激素和重金屬，會(huì)使水生生物遺傳變異、生理機(jī)能失常、機(jī)體變形甚至死亡。因此，隨著水環(huán)境綜合管理需求的凸顯，水沙－水環(huán)境－水生態(tài)多因素耦合模擬技術(shù)是未來(lái)的發(fā)展趨勢(shì)之一。

水體富營(yíng)養(yǎng)化是水動(dòng)力、營(yíng)養(yǎng)鹽、光照、水溫等因素綜合作用的結(jié)果，富營(yíng)養(yǎng)化模擬屬于一種典型的水環(huán)境多因素耦合模擬。目前，關(guān)于湖庫(kù)的富營(yíng)養(yǎng)化模擬已開(kāi)展了較多研究［27－28］。隨著河湖藻類(lèi)生長(zhǎng)機(jī)理以及水華發(fā)生機(jī)制的進(jìn)一步揭示，未來(lái)的趨勢(shì)是引入水動(dòng)力、水質(zhì)與富營(yíng)養(yǎng)化的關(guān)系，進(jìn)一步開(kāi)發(fā)三維富營(yíng)養(yǎng)化數(shù)學(xué)模型，更精準(zhǔn)地模擬營(yíng)養(yǎng)鹽和藻類(lèi)在水平和垂向的時(shí)空分布規(guī)律，預(yù)測(cè)水華發(fā)生的概率、區(qū)域及規(guī)模。

4.3 模擬系統(tǒng)集成

模型的前處理與后處理是水環(huán)境數(shù)值模擬工作的重要環(huán)節(jié)，其中前處理包括地形、水文、水質(zhì)數(shù)據(jù)的整理與輸入，后處理包括模擬結(jié)果的動(dòng)態(tài)展示和仿真演示。針對(duì)水環(huán)境高效管理的需求，有必要將模擬前后處理與數(shù)學(xué)模型集成于一體，建立水環(huán)境模擬系統(tǒng)，為相關(guān)的決策管理提供支撐。

根據(jù)側(cè)重點(diǎn)不同，可相應(yīng)開(kāi)發(fā)水質(zhì)信息化管理系統(tǒng)、水質(zhì)預(yù)警預(yù)報(bào)系統(tǒng)、突發(fā)性水污染事件應(yīng)急系統(tǒng)、流域/區(qū)域地下水信息及環(huán)境服務(wù)系統(tǒng)、水動(dòng)力與水環(huán)境演變?nèi)S動(dòng)態(tài)可視化演示系統(tǒng)、水環(huán)境實(shí)體景觀設(shè)計(jì)可視化演示系統(tǒng)等［29］。

5 結(jié)論

流域水環(huán)境數(shù)值模擬是水環(huán)境規(guī)劃、管理與保護(hù)工作中的重要研究手段。本文總結(jié)了目前流域水環(huán)境數(shù)值模擬技術(shù)及應(yīng)用研究進(jìn)展，并對(duì)水環(huán)境模擬技術(shù)的發(fā)展趨勢(shì)進(jìn)行了探討。

(1)根據(jù)模擬對(duì)象的不同，將水環(huán)境模型分為地表水環(huán)境模型、地下水環(huán)境模型、流域非點(diǎn)源污染模型3種類(lèi)型，并詳細(xì)列舉了各類(lèi)型的代表模擬軟件及其優(yōu)缺點(diǎn)。

(2)目前水環(huán)境模擬技術(shù)主要應(yīng)用在基礎(chǔ)研究、工程應(yīng)用和規(guī)劃管理等方面。①基礎(chǔ)研究:流域水環(huán)境演變與調(diào)控，主要包括全球氣候變化與人類(lèi)活動(dòng)影響下流域水環(huán)境的演變模擬，流域水沙調(diào)控及生態(tài)環(huán)境響應(yīng)模擬等。②工程應(yīng)用:水利工程影響分析與水污染防治工程設(shè)計(jì)，主要包括水庫(kù)、水系連通、跨流域調(diào)水工程等重大水利工程運(yùn)行對(duì)流域水環(huán)境影響模擬，非點(diǎn)源污染控制中生態(tài)溝渠和人工濕地方案設(shè)計(jì)等。③規(guī)劃管理:水資源保護(hù)規(guī)劃與管理。主要包括水功能區(qū)納污能力復(fù)核、水功能區(qū)達(dá)標(biāo)分析、入河排污口設(shè)置論證等。

(3)水環(huán)境數(shù)值模擬技術(shù)未來(lái)發(fā)展趨勢(shì)主要集中在3方面:一是陸面、地表水與地下水水環(huán)境集成的流域水環(huán)境全過(guò)程模擬;二是水沙－水環(huán)境－水生態(tài)多因素耦合模擬;三是將模擬前后處理與數(shù)學(xué)模型集成于一體的水環(huán)境模擬系統(tǒng)開(kāi)發(fā)。

［1］馮民權(quán)，鄭邦民，周孝德.水環(huán)境模擬與預(yù)測(cè)［M］.北京:科學(xué)出版社，2009.(FENG Min-quan，ZHENG Bang-min，ZHOU Xiao-de.Simulation and Prediction of Water Environment［M］.Beijing:Science Press，2009.(in Chinese))

［2］陸垂裕，肖偉華，趙 勇，等.復(fù)雜水環(huán)境系統(tǒng)數(shù)值模擬及其風(fēng)險(xiǎn)分析［M］.北京:中國(guó)水利水電出版社，2010.(LU Chui-yu，XIAO Wei-hua，ZHAO Yong，et al.Simulation and Risk Analysis of Complex Water Environment［M］.Beijing:China Water Power Press，2010.(in Chinese))

［3］彭澤洲，楊天行，梁秀娟，等.水環(huán)境數(shù)學(xué)模型及其應(yīng)用［M］.北京:化學(xué)工業(yè)出版社，2007.(PENG Zezhou，YANG Tian-xing，LIANG Xiu-juan，et al.Mathematical Model of Water Environment and Its Application［M］.Beijing:Chemical Industry Press，2007.(in Chinese))

［4］段 揚(yáng).基于EFDC的丹江口水庫(kù)水環(huán)境數(shù)值模擬分析［D］.北京:中國(guó)地質(zhì)大學(xué)(北京)，2014.(DUAN Yang.EFDC-Based Simulation on Water Environment in Danjiangkou Reservoir［D］.Beijing:China University of Geosciences(Beijing)，2014.(in Chinese))

［5］TetraTech，Inc.User’s Manual For Environmental Fluid Dynamics Code Hydro Version(EFDC-Hydro)Release 1.00［K］.USA:Atlanta，2002.

［6］TIM A，ROBERT B，JAMESL，et al.Water Quality A-nalysis Simulation Program(WASP)Version 6.0，Draft:User’s Manual［K］.USA:Atlanta，2001.

［7］DHI Water＆ Environment.MIKE 11 A Modeling System for Rivers and Channels User Guide［K］.Australia:DHI Water＆ Environment，2004.

［8］WL/Delft Hydraulics.Delft3D-WAQVersatile Water Quality Modelling in 1D，2D or 3D Systems Including Physical，(Bio)chemical and Biological Processes User Manual［K］.Netherlands:WL/Delft Hydraulics，2007.

［9］Environmental Modeling Research Laboratory of Brigham Young University.SMS Surface Water Modeling System Tutorials Version 8.1［K］.USA:Environmental Modeling Research Laboratory of Brigham Young University，2003.

［10］US Army Engineer Waterways Experiment Station.CEQUAL-R1:A Numerical One-dimensional Model of Reservoir Water Quality:User’s Manual［K］.Vicksburg，Mississippi:Environmental Laboratory，1986.

［11］張 昊，張代鈞.復(fù)雜水環(huán)境模擬研究與發(fā)展趨勢(shì)［J］.環(huán)境科學(xué)與管理，2010，(4):24－28.(ZHANG Hao，ZHANG Dai-jun.The Advance and Trend in the Simulation of Complicated Water Environment［J］.Environmental Science and Management，2010，(4):24－28.(in Chinese))

［12］鄭春苗，貝聶特.地下水污染物遷移模擬(第2版)［M］.北京:高等教育出版社，2009.(ZHENG Chunmiao，BENNEFF G D.Applied Contaminant Transport Modeling(Second Edition)［M］.Beijing:Higher Education Press，2009.(in Chinese))

［13］ANON D.Visual MODFLOW V.2.8.2 User’s Manual for Professional Applications in Three-Dimensional Groundwater Flow and Contaminant Transport Modeling［K］.Ontario:Waterloo Hydrogeologic Inc.，2000.

［14］夏 軍，翟曉燕，張永勇.水環(huán)境非點(diǎn)源污染模型研究進(jìn)展［J］.地理科學(xué)進(jìn)展，2012，(7):941－952.(XIA Jun，ZHAI Xiao-yan，ZHANG Yong-yong.Progress in the Research of Water Environmental Nonpoint Source Pollution Models［J］.Progress in Geography，2012，(7):941－952.(in Chinese))

［15］NEITSCH S L，ARNOLD J G，KINIRT J R.Soil and Water Assessment Tool Theoretical Documentation Version 2005［K］.Texas:Grassland，Soil and Water Research Laboratory，2005.

［16］BICKNELL B R，IMHOFF J C，KITTLE JR J L，et al.Hydrological Simulation Program-Fortran:User’s Manual for Version II［K］.Washington D C:National Exposure Research Laboratory，USEnvironmental Protection Agency，1997.

［17］ROSSMAN L A，SUPPLY W.Storm Water Management Model User’s Manual，Version 5.0［K］.Washington D C:National Risk Management Research Laboratory，Office of Research and Development，US Environmental Protection Agency，2010.

［18］曹麗娟，董文杰，張 勇.未來(lái)氣候變化對(duì)黃河和長(zhǎng)江流域極端徑流影響的預(yù)估研究［J］.大氣科學(xué)，2013，37(3):634－644.(CAO Li-juan，DONG Wenjie，ZHANG Yong.Estimation of the Effect of Climate Change on Extreme Streamflow over the Yellow River and Yangtze River Basins［J］.Chinese Journal of Atmospheric Sciences，2013，37(3):634－644.(in Chinese))

［19］夏星輝，吳 瓊，牟新利.全球氣候變化對(duì)地表水環(huán)境質(zhì)量影響研究進(jìn)展［J］.水科學(xué)進(jìn)展.2012，23(1):124－ 133.(XIA Xing-hui，WU Qiong，MOU Xin-li.Advances in Impacts of Climate Change on Surface Water Quality［J］.Advances in Water Science，2012，23(1):124－133.(in Chinese))

［20］陳明千，劉中峰，李 然，等.大型水庫(kù)三維水質(zhì)模型研究［J］.水利水電科技進(jìn)展，2010，30(2):5－9.(CHEN Ming-qian，LIU Zhong-feng，LI Ran，et al.Three-dimensional Model for Water Quality of Large Reservoirs［J］.Advances in Science and Technology of Water Resources，2010，30(2):5－9.(in Chinese))

［21］曹慧群，李青云，陳 野.三峽水庫(kù)運(yùn)行對(duì)長(zhǎng)江干流磷傳輸?shù)挠绊懠皩?duì)策［J］.人民長(zhǎng)江，2014，45(15):14－17.(CAO Hui-qun，LI Qing-yun，CHEN Ye.Impacts of Three Gorges Reservoir Operation on Phosphorus Transportation in Mainstream of Yangtze River and Countermeasure［J］.Yangtze River，2014，45(15):14－17.(in Chinese))

［22］田 勇.湖泊三維水動(dòng)力水質(zhì)模型研究與應(yīng)用［D］.武漢:華中科技大學(xué)，2012.(TIAN Yong.Development and Application of a Three Dimensional Hydrodynamic and Water Quality Lake Model［D］.Wuhan:Huazhong University of Science and Technology，2012.(in Chinese))

［23］王振華，李青云，黃 茁.生態(tài)清潔小流域建設(shè)研究現(xiàn)狀及展望［J］.人民長(zhǎng)江，2011，42(增2):115－118.(WANG Zheng-hua，LI Qing-yun，HUANG Zhuo.Research Status and Prospect on Construction of Small Eco-clean Watersheds ［J］. Yangtze River， 2011，42(Sup.2):115－118.(in Chinese))

［24］葉芝菡，段淑懷，吳敬東.流域水文模型在生態(tài)清潔小流域規(guī)劃中的應(yīng)用［J］.中國(guó)水土保持，2007，(9):12－ 13.(YE Zhi-han，DUAN Shu-huai，WU Jingdong.Application of Hydrological Model at Small Ecoclean Watersheds Planning［J］.Soil and Water Conservation in China，2007，(9):12－13.(in Chinese))

［25］趙 鑫，黃 茁，李青云，等.我國(guó)現(xiàn)行水域納污能力計(jì)算方法的思考［J］.中國(guó)水利，2012，(1):29－32.(ZHAO Xin，HUANG Zhuo，LI Qing-yun，et al.Consideration on Current Calculating Method of Pollution Carrying Capacity in China［J］.China Water Resources，2012，(1):29－32.(in Chinese))

［26］凌敏華，陳 喜，程勤波，等.地表水與地下水耦合模型研究進(jìn)展［J］.水利水電科技進(jìn)展，2010，30(4):79－84.(LING Min-hua，CHEN Xi，CHENG Qin-bo，et al.Advances in Coupled Surface Water and Groundwater Models［J］.Advances in Science and Technology of Water Resources，2010，30(4):79－ 84.(in Chinese))

［27］LI Y P，ACHARYA K，YU Z B.Modeling Impacts of Yangtze River Water Transfer on Water Ages in Lake Taihu，China［J］.Ecological Engineering，2011，37(2):325－334.

［28］唐天均，楊 晟，尹魁浩，等.基于EFDC模型的深圳水庫(kù)富營(yíng)養(yǎng)化模擬［J］.湖泊科學(xué)，2014，26(3):393－400.(TANG Tian-jun，YANG Sheng，YIN Kuihao，et al.Simulation of Eutrophication in Shenzhen Reservoir Based on EFDC Model［J］.Journal of Lake Sciences，2014，26(3):393－400.(in Chinese))

［29］陳 進(jìn)，李青云.長(zhǎng)江流域水環(huán)境綜合治理的技術(shù)支撐體系探討［J］.人民長(zhǎng)江，2011，42(2):94－97.(CHEN Jin，LI Qing-yun.Discussion on Technical Support System for Comprehensive Water Environment Harnessing in Yangtze River Basin［J］.Yangtze River，2011，42(2):94－97.(in Chinese))