陳 靜
(上海市政工程設(shè)計(jì)研究總院(集團(tuán))有限公司,上海 200092)
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ā)展階段。
歐美國(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)性,可用于描述土壤蓄水能力的空間分布。
面對(duì)一般的水文問(wèn)題,集總式水文模型也可得到相對(duì)令人滿(mǎn)意的結(jié)果,而分布式水文模型的優(yōu)勢(shì)并不十分突出。但在面對(duì)復(fù)雜的水文問(wèn)題時(shí),如缺乏資料地區(qū)、空間異質(zhì)性明顯的流域,以及在研究氣候變化對(duì)水文水資源的影響、土地利用對(duì)水文過(guò)程的響應(yīng)等問(wèn)題時(shí),分布式水文模型就表現(xiàn)出明顯的優(yōu)勢(shì)和不可替代性。
在實(shí)際研究和應(yīng)用中流域水文資料不足是經(jīng)常遇到的情況,所謂的缺乏資料地區(qū)不僅有因區(qū)域測(cè)站建設(shè)不完善導(dǎo)致的資料缺乏,也有由于下墊面和氣候條件變化導(dǎo)致以往的水文資料不再適用而成為新的缺乏資料地區(qū)。經(jīng)驗(yàn)?zāi)P秃图偸剿哪P蛯?duì)于地區(qū)觀測(cè)資料有著很強(qiáng)的依賴(lài)性,因此面對(duì)資料缺乏地區(qū)的水文問(wèn)題常常無(wú)從下手。分布式水文模型通過(guò)求解水文過(guò)程偏微分方程的方法,弱化了對(duì)經(jīng)驗(yàn)參數(shù)的依賴(lài)性,具有很強(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ò)程的方法是可行的。
全球正在經(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)精度。
人類(lèi)活動(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)越被重視,尤其是近年來(lái)遙感技術(shù)、地理信息系統(tǒng)與分布式水文模型的耦合不斷深入,使得分布式水文模型能夠更加靈活地反映土地利用(覆被)變化的情景,模擬不同土地利用變化情境下的水文響應(yīng)[13]。
分布式水文模型的發(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)沙量模擬,獲得了令人滿(mǎ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ù)和治理方向。
在水資源利用程度不斷提高而水資源危機(jī)不斷加劇的背景下,水資源量與質(zhì)并重,管理更加精細(xì)化的要求使得傳統(tǒng)的水資源管理方式已不能滿(mǎn)足需要,分布式水文模型較好地適應(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]。
在過(guò)去的研究進(jìn)程中,分布式水文模型的研究和應(yīng)用獲得了較大的發(fā)展,但亦存在一些問(wèn)題亟待獲得突破和發(fā)展。
(1)能夠適應(yīng)不同尺度的新一代分布式水文模型研究和應(yīng)用亟待出現(xiàn)。尺度的問(wèn)題包含兩類(lèi),一類(lèi)是對(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)和方向之一;另一類(lèi)則是分布式水文模型與相關(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)性,一直以來(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)化。
分布式水文模型與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ī)理。