流域分布式水文模型的發(fā)展現(xiàn)狀與展望

陳 靜

（上海市政工程設(shè)計(jì)研究總院（集團(tuán)）有限公司，上海 200092）

1 引言

20世紀(jì)60年代出現(xiàn)的斯坦福（Stanford）模型是第一個(gè)真正意義上的流域水文模型，是計(jì)算機(jī)科學(xué)與水文學(xué)成功結(jié)合的產(chǎn)物，標(biāo)志著定量水文學(xué)在經(jīng)歷了漫長(zhǎng)而緩慢的發(fā)展之后開(kāi)始進(jìn)入快速發(fā)展階段[1]。在過(guò)去的半個(gè)多世紀(jì)中，水文模型的發(fā)展經(jīng)歷了從黑箱子、概念性模型向具有物理基礎(chǔ)的分布式水文模型，從集總式向分布式，從僅考慮降雨徑流模擬到耦合水質(zhì)、泥沙、生態(tài)系統(tǒng)及大氣系統(tǒng)等的發(fā)展歷程[2]。

根據(jù)對(duì)流域空間離散程度的不同描述方法，水文模型可分為集總式和分布式。其中集總式水文模型將流域視作一個(gè)整體，假定整個(gè)流域上降雨、流域下墊面和初始狀態(tài)在空間上均勻分布，不考慮三者的空間差異性；而分布式水文模型則充分考慮了降雨、流域下墊面和初始狀態(tài)在空間上的不均勻性，模型在求解描述水文循環(huán)過(guò)程中的偏微分方程時(shí)，采用離散化方法將流域分成多個(gè)單元，并考慮水分在離散單元之間的水量和動(dòng)量交換。由于這種假設(shè)與自然界下墊面的空間異質(zhì)性和降雨的時(shí)空不均勻性，而導(dǎo)致降雨徑流非線性的特征比較吻合，因此分布式水文模型對(duì)于水文循環(huán)過(guò)程的描述更加接近實(shí)際，是水文模型發(fā)展的必然趨勢(shì)[3-4]。1969年，F(xiàn)reeze和Harlan[5]第一次提出分布式物理模型的概念后，分布式水文模型的開(kāi)發(fā)和應(yīng)用開(kāi)始進(jìn)入快速發(fā)展階段。

2 分布式水文模型的研究現(xiàn)狀

歐美國(guó)家研究機(jī)構(gòu)在分布式水文模型方面起步較早，在這方面的研究一直走在世界前列。目前應(yīng)用較多的分布式水文模型也是以歐美國(guó)家或科研機(jī)構(gòu)開(kāi)發(fā)為主，主要包括歐洲科研機(jī)構(gòu)聯(lián)合開(kāi)發(fā)的SHE（System Hydrological European）模型、美國(guó)農(nóng)業(yè)部研發(fā)的SWAT（Soiland Water Assessment Tool）模型、英國(guó)水文學(xué)家開(kāi)發(fā)的IHDM（Insistute of Hydrology Distributed Model）模型、意大利水文學(xué)家開(kāi)發(fā)的TOPKAP（Topographic Kinematic Approximation and Integration）模型、美國(guó)科研機(jī)構(gòu)開(kāi)發(fā)的VIC（Variable Infiltration Capacity）模型等。

1986年,由歐洲多個(gè)國(guó)家水文學(xué)家聯(lián)合開(kāi)發(fā)的SHE水文模型是最早出現(xiàn)的分布式水文模型，也是應(yīng)用最廣泛的分布式水文模型之一，自出現(xiàn)以來(lái)在全世界不同國(guó)家和流域均有較為成功的應(yīng)用實(shí)踐。SHE模型是迄今為止唯一一個(gè)全部用連續(xù)性方程和動(dòng)力方程控制流域產(chǎn)匯流的流域水文模型，然而對(duì)于這樣一個(gè)包含復(fù)雜定解條件的偏微分方程組，得到解析解是相當(dāng)困難。SHE模型利用有限差分法進(jìn)行數(shù)值求解，在流域網(wǎng)格化方面SHE模型在水平方向上用網(wǎng)格表示，垂直方向上用層來(lái)表示。SHE模型的部分參數(shù)具有一定的物理意義，跟流域特征有關(guān)。數(shù)值化求解方法使得SHE模型在計(jì)算上具有較強(qiáng)的靈活性，而表征流域特性的參數(shù)亦使其具有一定的物理基礎(chǔ)，上述特點(diǎn)使SHE模型適用于多種資料條件，在歐洲和其他地區(qū)得到了廣泛應(yīng)用。

美國(guó)農(nóng)業(yè)部于1988年開(kāi)發(fā)的SWAT模型是目前應(yīng)用最廣泛的分布式水文模型之一，是一個(gè)具有很強(qiáng)物理基礎(chǔ)的長(zhǎng)時(shí)段流域分布式水文模型，適用于下墊面和管理?xiàng)l件差異性明顯的復(fù)雜大流域[6]。SWAT在模型架構(gòu)上采用模塊化設(shè)計(jì)思路，模型功能強(qiáng)大，不僅能夠模擬降雨徑流過(guò)程，而且在水土流失、非點(diǎn)源污染和水資源管理等問(wèn)題上也得到了很好的應(yīng)用。

分布式水文模型在國(guó)內(nèi)發(fā)展較晚，國(guó)內(nèi)對(duì)于流域水文模型的探索是應(yīng)用與開(kāi)發(fā)相結(jié)合的模式，一方面積極引入國(guó)外相對(duì)成熟的流域水文模型，嘗試將其應(yīng)用于國(guó)內(nèi)不同流域的水文模擬、洪水預(yù)報(bào)、水資源管理等水文過(guò)程，同時(shí)驗(yàn)證分布式水文模型對(duì)流域的適應(yīng)性，獲得符合流域?qū)嶋H的參數(shù)，服務(wù)于流域的研究和應(yīng)用實(shí)踐；另一方面也積極開(kāi)發(fā)新的流域水文模型。

張建云[7]建立了一種參數(shù)網(wǎng)格化的分布式徑流模型，并在華北、江淮流域的水資源動(dòng)態(tài)模擬評(píng)估中得到較好應(yīng)用；李蘭等[8]提出了一種分布式水文模型，將數(shù)學(xué)物理反問(wèn)題與洪水預(yù)報(bào)相結(jié)合，給出了包含流域產(chǎn)匯流和水庫(kù)洪水演進(jìn)的動(dòng)態(tài)分布預(yù)報(bào)耦合模型，不僅可以用于分析地區(qū)的降雨徑流規(guī)律，也可以用于洪水預(yù)報(bào)，并在多個(gè)水庫(kù)流域得到應(yīng)用；熊立華等[9]提出了一個(gè)基于數(shù)字高程模型（DEM）的分布式水文模型（DDRM），將土壤蓄水能力作為反映流域土壤特征的參數(shù)，同時(shí)假設(shè)各點(diǎn)地形指數(shù)與土壤蓄水能力存在相關(guān)性，可用于描述土壤蓄水能力的空間分布。

3 分布式水文模型的應(yīng)用

面對(duì)一般的水文問(wèn)題，集總式水文模型也可得到相對(duì)令人滿意的結(jié)果，而分布式水文模型的優(yōu)勢(shì)并不十分突出。但在面對(duì)復(fù)雜的水文問(wèn)題時(shí)，如缺乏資料地區(qū)、空間異質(zhì)性明顯的流域，以及在研究氣候變化對(duì)水文水資源的影響、土地利用對(duì)水文過(guò)程的響應(yīng)等問(wèn)題時(shí)，分布式水文模型就表現(xiàn)出明顯的優(yōu)勢(shì)和不可替代性。

3.1 分布式水文模型在缺少資料地區(qū)水文模擬中的應(yīng)用

在實(shí)際研究和應(yīng)用中流域水文資料不足是經(jīng)常遇到的情況，所謂的缺乏資料地區(qū)不僅有因區(qū)域測(cè)站建設(shè)不完善導(dǎo)致的資料缺乏，也有由于下墊面和氣候條件變化導(dǎo)致以往的水文資料不再適用而成為新的缺乏資料地區(qū)。經(jīng)驗(yàn)?zāi)Ｐ秃图偸剿哪Ｐ蛯?duì)于地區(qū)觀測(cè)資料有著很強(qiáng)的依賴性，因此面對(duì)資料缺乏地區(qū)的水文問(wèn)題常常無(wú)從下手。分布式水文模型通過(guò)求解水文過(guò)程偏微分方程的方法，弱化了對(duì)經(jīng)驗(yàn)參數(shù)的依賴性，具有很強(qiáng)的普適性，因而在缺少觀測(cè)資料地區(qū)的水文模擬中得到了廣泛應(yīng)用，也成為缺乏資料地區(qū)模擬降雨徑流過(guò)程、進(jìn)行水文預(yù)報(bào)的一種有效的工具。黃粵等[10]選取具有物理基礎(chǔ)的分布式水文模型MIKE SHE來(lái)模擬塔里木河主要源區(qū)之一——開(kāi)都河流域的徑流量變化過(guò)程，在觀測(cè)缺乏的條件下,研究分布式水文模型在缺乏資料地區(qū)的日徑流模擬中的適用性，結(jié)果表明模擬結(jié)果與實(shí)測(cè)結(jié)果高度相關(guān)，在資料相對(duì)缺乏的客觀條件下，應(yīng)用分布式水文模型來(lái)模擬地區(qū)水文過(guò)程的方法是可行的。

3.2 分布式水文模型在研究水文循環(huán)對(duì)氣候變化響應(yīng)機(jī)理中的應(yīng)用

全球正在經(jīng)歷的氣候變化直接影響了降雨的時(shí)空分布，局部極端天氣發(fā)生的頻率和降雨強(qiáng)度的明顯增強(qiáng)，地區(qū)降雨徑流特性和水資源量也在不斷變化，因此研究氣候變化對(duì)水文水資源的影響也成為近年來(lái)水文學(xué)的熱點(diǎn)之一，不少學(xué)者將大氣環(huán)流模型與流域水文模型耦合而形成陸-氣耦合模式（GCMs），該模式成為研究水文循環(huán)對(duì)氣候變化響應(yīng)機(jī)理的重要方法。分布式水文模型對(duì)陸面的空間差異性有著更為真實(shí)的描述和表征，這一特點(diǎn)使其在陸-氣耦合模式中尤被重視。彭艷[11]針對(duì)三峽庫(kù)區(qū)流域特點(diǎn)，利用分布式水文模型與數(shù)值天氣預(yù)報(bào)模式進(jìn)行耦合，構(gòu)建了三峽庫(kù)區(qū)陸-氣耦合洪水預(yù)報(bào)模型，并開(kāi)發(fā)了陸-氣耦合的洪水預(yù)報(bào)系統(tǒng)，有效延長(zhǎng)了入庫(kù)洪水預(yù)報(bào)預(yù)見(jiàn)期，提高了洪水預(yù)報(bào)精度。

3.3 分布式水文模型在流域土地利用（覆被）變化對(duì)水文過(guò)程影響研究中的應(yīng)用

人類活動(dòng)對(duì)水文過(guò)程的影響主要表現(xiàn)為改變土地利用方式及陸面覆被情況，這種改變直接導(dǎo)致蒸散發(fā)、截留、填洼、下滲等因素的改變，從而改變流域（區(qū)域）產(chǎn)匯流的過(guò)程及流量，繼而導(dǎo)致流域（區(qū)域）不同時(shí)間和空間尺度下的水文循環(huán)的變化。土地利用（覆被）變化對(duì)水文過(guò)程影響的研究，對(duì)于揭示區(qū)域及全球尺度水文循環(huán)規(guī)律、相互影響機(jī)理、生態(tài)安全格局等有著十分重要的意義[12]。而分布式水文模型的特點(diǎn)使其在這一類研究中越來(lái)越被重視，尤其是近年來(lái)遙感技術(shù)、地理信息系統(tǒng)與分布式水文模型的耦合不斷深入，使得分布式水文模型能夠更加靈活地反映土地利用（覆被）變化的情景，模擬不同土地利用變化情境下的水文響應(yīng)[13]。

3.4 分布式水文模型在物質(zhì)輸移和水質(zhì)模擬中的廣泛應(yīng)用

分布式水文模型的發(fā)展從最初的單一模擬降雨徑流逐漸向耦合水質(zhì)、泥沙等物質(zhì)運(yùn)移的方向發(fā)展，國(guó)外學(xué)者較早就對(duì)分布式模型在物質(zhì)運(yùn)移方面的應(yīng)用進(jìn)行了探索，利用模型對(duì)流域的產(chǎn)流、產(chǎn)沙、非點(diǎn)源污染等過(guò)程進(jìn)行了大量的模擬研究。Schomberg等[14]應(yīng)用分布式水文模型在美國(guó)的不同州進(jìn)行了產(chǎn)流量、產(chǎn)沙量模擬，獲得了令人滿意的模擬結(jié)果。國(guó)內(nèi)對(duì)于分布式水文模型的研究和應(yīng)用起步相對(duì)較晚，但近年來(lái)在應(yīng)用分布式水文模型分析產(chǎn)流、產(chǎn)沙量及非點(diǎn)源污染負(fù)荷等方面也有較多的應(yīng)用實(shí)踐。解志林[15]利用SWAT模型研究了阜陽(yáng)市沙潁河流域非點(diǎn)源總氮的時(shí)空分布特征，不僅分析確認(rèn)了農(nóng)業(yè)面源污染物的重點(diǎn)貢獻(xiàn)區(qū)域，同時(shí)研究了降雨、下滲、地表徑流等各個(gè)水文環(huán)節(jié)總氮進(jìn)入水體的變化過(guò)程；陳祥義[16]選取三峽庫(kù)區(qū)典型流域龍河流域作為研究區(qū)域，利用分布式水文模型SWAT對(duì)研究區(qū)內(nèi)的非點(diǎn)源污染進(jìn)行模擬，找到污染物排放關(guān)鍵源區(qū)，同時(shí)探討了非點(diǎn)源污染防控最佳管理措施，提高非點(diǎn)源污染治理的效率；楊寶林等[17]利用分布式水文模型模擬分析了浙江省蓮塘口流域的關(guān)鍵污染源區(qū)，為流域后續(xù)的面源污染控制和治理提供科學(xué)依據(jù)和治理方向。

3.5 分布式水文模型在水資源管理中的應(yīng)用

在水資源利用程度不斷提高而水資源危機(jī)不斷加劇的背景下，水資源量與質(zhì)并重，管理更加精細(xì)化的要求使得傳統(tǒng)的水資源管理方式已不能滿足需要，分布式水文模型較好地適應(yīng)了這一發(fā)展要求。分布式水文模型已經(jīng)在流域水資源綜合管理、農(nóng)業(yè)灌溉和城市取用水研究、地表水與地下水評(píng)價(jià)與計(jì)算等領(lǐng)域發(fā)揮了重大作用。王忠靜等[18]利用分布式水文模型應(yīng)用于變化環(huán)境下的水資源評(píng)價(jià)，以永定河山區(qū)為例進(jìn)行了不同時(shí)期下墊面的水資源評(píng)價(jià)，驗(yàn)證了分布式水文模型用于水資源評(píng)價(jià)的有效性。在南水北調(diào)工程中，分析人員應(yīng)用分布式水文模型進(jìn)行情景分析，研究水文系統(tǒng)的動(dòng)態(tài)變化，對(duì)地下水資源進(jìn)行評(píng)估與計(jì)算，從而制定地下水恢復(fù)的方案[19]。

4 存在的一些問(wèn)題

在過(guò)去的研究進(jìn)程中，分布式水文模型的研究和應(yīng)用獲得了較大的發(fā)展，但亦存在一些問(wèn)題亟待獲得突破和發(fā)展。

（1）能夠適應(yīng)不同尺度的新一代分布式水文模型研究和應(yīng)用亟待出現(xiàn)。尺度的問(wèn)題包含兩類，一類是對(duì)不同時(shí)間和空間尺度的適應(yīng)性。不同尺度的水文循環(huán)機(jī)理、水文規(guī)律和特征均有不同，所使用的水文模型的結(jié)構(gòu)也就不盡相同。開(kāi)發(fā)一種模型結(jié)構(gòu)能夠適應(yīng)于不同尺度的分布式水文模型，使得模型對(duì)于不同尺度的水文問(wèn)題都有著較好的表達(dá)，是水文模型進(jìn)一步發(fā)展的主要難點(diǎn)和方向之一；另一類則是分布式水文模型與相關(guān)模型耦合時(shí)的尺度協(xié)調(diào)問(wèn)題。例如，水文模型與生態(tài)系統(tǒng)模型的耦合，水文模型與大氣環(huán)流模型的耦合等，均存在時(shí)間和空間尺度的協(xié)調(diào)問(wèn)題，如何實(shí)現(xiàn)模型之間尺度的轉(zhuǎn)換是耦合模型進(jìn)一步發(fā)展的重點(diǎn)和難點(diǎn)之一。

（2）遙感與地理信息系統(tǒng)、分布式水文模型的耦合與集成有待進(jìn)一步發(fā)展。分布式水文模型的特點(diǎn)決定了其對(duì)地理信息有著較高的依賴性，一直以來(lái)受數(shù)據(jù)來(lái)源及數(shù)據(jù)密度的制約，并沒(méi)有出現(xiàn)真正意義上的分布式水文模型。地理信息系統(tǒng)（GIS）和遙感技術(shù)（RS）的發(fā)展為分布式水文模型解決了數(shù)據(jù)源及密度的阻礙，RS技術(shù)強(qiáng)大的動(dòng)態(tài)數(shù)據(jù)采集以及GIS數(shù)據(jù)管理和處理技術(shù)在分布式水文模型中的交互式應(yīng)用，使得真正意義上分布式水文模型的出現(xiàn)成為可能。一般來(lái)說(shuō)分布式水文模型與地理信息系統(tǒng)的耦合有4種形式[20]，模型嵌入GIS平臺(tái)、GIS功能集成到模型、模型和GIS的松散連接、模型和GIS的緊密連接。松散連接是目前分布式水文模型與GIS集成研究中最普遍的方法，通過(guò)數(shù)據(jù)文件或者共享數(shù)據(jù)庫(kù)的方式借助數(shù)據(jù)接口程序來(lái)實(shí)現(xiàn)二者的結(jié)合，接口程序相對(duì)簡(jiǎn)單，運(yùn)行效率高，模型的穩(wěn)定性較好，但是缺點(diǎn)在于易用性較差，經(jīng)常會(huì)出現(xiàn)數(shù)據(jù)交換錯(cuò)誤，從而影響到模型的精度及效果。通過(guò)地理信息系統(tǒng)的宏語(yǔ)言或二次平發(fā)，將分布式水文模型與地理信息系統(tǒng)緊密連接，實(shí)現(xiàn)二者的完全集成，是模型耦合集成的難點(diǎn)和發(fā)展方向。

（3）分布式水文模型在參數(shù)率定及運(yùn)行方面有待進(jìn)一步優(yōu)化。分布式水文模型的參數(shù)眾多，多參數(shù)率定一直以來(lái)都是模型的難點(diǎn)，自動(dòng)率定成為一種必然的發(fā)展趨勢(shì)。利用粒子群算法、神經(jīng)網(wǎng)絡(luò)及人工智能等方法應(yīng)用于模型的多參數(shù)優(yōu)化，實(shí)現(xiàn)多參數(shù)的自動(dòng)率定，是模型提升計(jì)算效率的重要手段。另外，在運(yùn)行方面現(xiàn)有的分布式水文模型多為串行計(jì)算方式，但隨著分布式水文模型在大流域、高密度數(shù)據(jù)源條件下的應(yīng)用，串行計(jì)算在運(yùn)行方面會(huì)遇到較大的瓶頸，并行計(jì)算的研究和應(yīng)用能夠幫助分布式水文模型實(shí)現(xiàn)在運(yùn)行效率上的優(yōu)化。

5 發(fā)展展望

分布式水文模型與RS、GIS的緊密耦合和完全集成是未來(lái)研究的熱點(diǎn)之一，在開(kāi)發(fā)性體系架構(gòu)的GIS平臺(tái)基礎(chǔ)上，通過(guò)GIS提供的宏語(yǔ)言或者二次開(kāi)發(fā)高級(jí)語(yǔ)言，將水文模型和GIS緊密連接，實(shí)現(xiàn)分布式模型和GIS的完全集成。完全集成的分布式水文模型可實(shí)現(xiàn)對(duì)水文動(dòng)態(tài)的實(shí)時(shí)模擬及演示，而且穩(wěn)定性和易用性均較好，是分布式水文模型發(fā)展的一個(gè)方向。

分布式水文模型在多學(xué)科的交叉融合中發(fā)揮更大作用。分布式水文模型與大氣環(huán)流模型耦合發(fā)展起來(lái)的陸-氣模式也是學(xué)科交叉的一個(gè)熱點(diǎn)和難點(diǎn)。新一代陸-氣耦合模型的出現(xiàn)，將促進(jìn)人們對(duì)水文循環(huán)與氣候演變交互影響機(jī)理的進(jìn)一步認(rèn)識(shí)。而在生態(tài)水文學(xué)領(lǐng)域，隨著分布式水文模型與相關(guān)模型在尺度匹配和參數(shù)傳遞等方面的突破，模型之間的耦合從松散向緊密耦合，植被生長(zhǎng)過(guò)程與水文循環(huán)之間的作用與反饋模擬將得到進(jìn)一步發(fā)展，從而不斷提高生物動(dòng)力與水文循環(huán)相互作用機(jī)理。