姚 成,李致家,章玉霞

(河海大學(xué)水文水資源學(xué)院,江蘇南京 210098)

DEM分辨率對分布式水文模擬的影響

姚 成,李致家,章玉霞

(河海大學(xué)水文水資源學(xué)院,江蘇南京 210098)

以數(shù)字高程模型為基礎(chǔ),采用柵格型新安江模型進(jìn)行不同分辨率情況下的洪水模擬,分析了分辨率對分布式水文模擬的影響。研究結(jié)果表明:數(shù)字高程模型分辨率越高,流域數(shù)字化的精度相對越高,但模型的模擬結(jié)果并不一定越好;針對不同分辨率的數(shù)字高程模型,在重新率定匯流參數(shù)的情況下,柵格型新安江模型的模擬精度相當(dāng),均能取得較好的應(yīng)用效果;對于中尺度流域而言,采用30″分辨率的數(shù)字高程模型不僅能夠提高模型的運行時效,也能取得較高的模擬精度。

數(shù)字高程模型;柵格型新安江模型;水文模擬;分辨率;參數(shù)估計

流域水文模型作為一種模擬流域內(nèi)復(fù)雜水文現(xiàn)象的重要工具,一直是水文學(xué)研究的熱點課題之一。就目前水文模型發(fā)展趨勢而言,分布式水文模型由于其能夠更好地考慮氣候和下墊面因子的空間變異性,能夠更好地利用地理信息系統(tǒng)(GIS)與遙感(RS)技術(shù)等空間信息模擬流域的降雨徑流響應(yīng),現(xiàn)已成為水文模型發(fā)展的重要方向[1-2]。

在應(yīng)用分布式水文模型時,往往需要流域內(nèi)部詳細(xì)的水流流向、河網(wǎng)水系、流域形狀等地形地貌特征,需要對研究流域進(jìn)行數(shù)字化。隨著GIS的應(yīng)用與計算機能力的提高,利用數(shù)字高程模型(DEM)數(shù)字化流域和提取流域地貌特征的技術(shù)已得到了廣泛研究,DEM的引入大幅度促進(jìn)了分布式水文模型的進(jìn)一步發(fā)展[3]。然而,在DEM平臺上構(gòu)建分布式水文模型,DEM分辨率對流域地貌特征的提取結(jié)果以及產(chǎn)匯流模擬精度有著重要的影響[4-5]。理論上, DEM分辨率越高,所提取的流域地貌特征越接近于真實情況,但對基于DEM的分布式水文模型而言,高分辨率的DEM會使得模型的計算負(fù)擔(dān)過重,在很大程度上制約了模型的運行時效,這必將導(dǎo)致分布式水文模型在實際應(yīng)用中受到限制。為了更真實高效地描述現(xiàn)實世界中的降雨徑流形成機理,有必要進(jìn)一步開展DEM分辨率對分布式水文模擬影響的研究。

本文以皖南山區(qū)呈村流域為例,在3″、12″、21″與30″共4種DEM分辨率情況下,對研究流域進(jìn)行了數(shù)字化,分析了分辨率對地貌特征提取結(jié)果的影響,并在此基礎(chǔ)上,采用基于DEM的柵格型新安江模型[6-8]進(jìn)行呈村流域的洪水模擬,探討了不同DEM分辨率對模型參數(shù)以及模擬結(jié)果的影響。

1 柵格型新安江模型

1.1 模型原理

柵格型新安江模型以流域內(nèi)每個DEM柵格作為計算單元,并假設(shè)降雨和下墊面特征在每個單元內(nèi)空間分布均勻,模型只考慮各要素在不同單元之間的變異性。模型構(gòu)建以蓄滿產(chǎn)流原理為基礎(chǔ),其蒸散發(fā)、產(chǎn)流與分水源計算方法和新安江模型一致,采用概念性方法。因此,柵格型新安江模型屬于松散結(jié)構(gòu)分布式水文模型[1]。

在模擬洪水時,對于任意柵格單元,在考慮柵格間水量交換與河道排水網(wǎng)絡(luò)影響的基礎(chǔ)上,分別采用3層蒸散發(fā)模型與蓄滿產(chǎn)流原理對單元內(nèi)的蒸散發(fā)E與產(chǎn)流量R進(jìn)行計算,并采用自由水蓄水庫結(jié)構(gòu)對水源進(jìn)行劃分,包括地表徑流Rs、壤中流Ri以及地下徑流Rg3種水源,再根據(jù)柵格間的匯流演算次序,依次將各種水源演算至流域出口。對于單元間的匯流過程模擬,采用的是逐柵格馬斯京根匯流演算方法[6-9]。圖1為柵格單元內(nèi)產(chǎn)流與分水源示意圖,圖中P為降雨量;WM、SM分別為張力水、自由水蓄水容量;WUM、WLM分別為上層、下層張力水蓄水容量;S為實際的自由水含量;Kg、Ki分別為地下水徑流、壤中水徑流的出流系數(shù);θwp、θfc、θs分別為凋萎含水量、田間持水量與飽和含水量。

圖1 柵格單元內(nèi)產(chǎn)流與分水源示意圖

1.2 模型參數(shù)

柵格型新安江模型以新安江模型為基礎(chǔ),其參數(shù)與新安江模型參數(shù)基本一致,如圖1中的WM、SM、Kg與Ki等。不同的是,當(dāng)利用該模型進(jìn)行分布式水文模擬時,不能完全以集總的方式對模型參數(shù)賦值,而是要獲得參數(shù)的空間分布。因此,不能單一地根據(jù)觀測資料來率定參數(shù),而應(yīng)該在率定參數(shù)之前,先根據(jù)參數(shù)的物理意義,基于地形地貌、土壤、植被等下墊面條件,結(jié)合理論推導(dǎo),合理估計出參數(shù)的空間分布[10-12]?？陀^的估計方法不僅有利于甄選出敏感參數(shù),降低參數(shù)之間的相依性,簡化參數(shù)的率定工作,還可以保證模型參數(shù)在整個流域上的相容性,以更真實地模擬流域的水文響應(yīng)[13]。

以WM、SM、Kg與Ki4個參數(shù)為例,WM表示的是張力水蓄水容量,對應(yīng)的土層為“全土層”,即整個包氣帶[14-15],而張力水可以認(rèn)為介于凋萎含水量與田間持水量之間,則有:

式中:La為包氣帶土層厚度。SM表示的是自由水蓄水容量,對應(yīng)的土層為“表土層”,即腐殖質(zhì)土層[14-15],而自由水可以認(rèn)為介于田間持水量與飽和含水量之間,則有:

式中:Lh為腐殖質(zhì)土層厚度。

利用式(1)和式(2)計算每個柵格單元蓄水容量時,土壤水分常數(shù)θwp、θfc、θs可以根據(jù)每個單元對應(yīng)的土壤類型估算得到,土層厚度La、Lh可以根據(jù)每個單元對應(yīng)的地形指數(shù)及植被覆蓋估算得到。利用式(1)和式(2)不僅可以獲得WM與SM的空間分布,也可以統(tǒng)計得到研究流域張力水與自由水容量曲線。

兩個出流系數(shù)Kg與Ki屬于并聯(lián)參數(shù),其和代表自由水出流的快慢,與單元內(nèi)的土壤類型有關(guān),可以將θs與θfc作為衡量自由水出流快慢的指標(biāo);Kg/Ki代表地下水與壤中流的比,該比值可以通過θwp來反映[16],則有:

式中:moc與mr為兩個校正系數(shù),根據(jù)已有的研究成果,在估計參數(shù)空間分布時,mr≈1[13-16],而moc可以通過分水源中的結(jié)構(gòu)性約束(即Kg+Ki=0.7)試算得到。

當(dāng)采用逐柵格馬斯京根法進(jìn)行單元間匯流演算時,對于其參數(shù)流量比重因子xe與蓄量常數(shù)ke的確定,模型暫沒有考慮其空間分布不均的問題,而是利用流域內(nèi)柵格單元統(tǒng)一賦值的方法,需要根據(jù)流域的觀測資料進(jìn)行率定和檢驗。ke屬于敏感參數(shù),反映的是河道水流、地表徑流、壤中流與地下徑流在柵格單元內(nèi)的匯流時間,由于不同的徑流成分在匯流速度上存在著差異,因此在率定ke時應(yīng)予以區(qū)分。

2 研究流域概況

呈村流域位于皖南山區(qū),鄰近東南沿海,四季分明,氣候溫和,屬亞熱帶季風(fēng)氣候區(qū),流域控制面積290km2。流域內(nèi)植被良好,雨量充沛,多年平均降雨量約為2 100 mm,流域降水在年內(nèi)年際分配極不均勻,為典型的濕潤流域。呈村流域的植被類型主要包括常綠針葉林、落葉闊葉林、混合林、森林地、林地草原、牧草地與作物地,土壤類型主要為黏壤土。流域內(nèi)有10個雨量站,具有1986—1999年32場洪水過程的時段資料,其中,前22場洪水過程用于模型參數(shù)的率定,后10場洪水過程用于模型參數(shù)的檢驗。

3 結(jié)果分析

3.1 地貌特征提取結(jié)果

根據(jù)3″分辨率的DEM數(shù)據(jù)[17],利用最近鄰算法生成12″、21″與30″的DEM數(shù)據(jù),在這4種分辨率情況下,對呈村流域進(jìn)行了數(shù)字化,并提取了流域地貌特征。圖2為根據(jù)上述4種分辨率的DEM生成的流域邊界及水系,生成水系時,對應(yīng)圖2(a)～(d)所用的閾值(網(wǎng)格數(shù))分別為100、25、15和10。

圖2 根據(jù)不同分辨率DEM生成的呈村流域邊界與水系

由圖2可以看出,DEM分辨率越高,所生成的流域邊界及水系越精細(xì)。分析地貌特征結(jié)果后發(fā)現(xiàn),分辨率的變化對所提取的呈村流域面積影響相對較小,變化幅度在1 km2以內(nèi)。此外,隨著DEM分辨率的降低,流域平均高程及河道平均坡度有增大的趨勢,而最大匯流路徑長度有減小的趨勢,三者變化幅度分別在13 m、0.6％、12 km以內(nèi)。

3.2 洪水模擬結(jié)果

以流域數(shù)字化與地貌特征提取結(jié)果為基礎(chǔ),將柵格型新安江模型用于呈村流域的洪水模擬。考慮到分辨率對模型運算效率有較大影響,先針對30″DEM,結(jié)合土壤、植被等下墊面資料,對模型產(chǎn)流參數(shù)的空間分布進(jìn)行了估計,并在此基礎(chǔ)上,對模型匯流參數(shù)進(jìn)行了率定,其中,河道水流、地表徑流、壤中流和地下徑流對應(yīng)的ke率定結(jié)果分別為0.06、0.10、0.60和1.00,xe的率定結(jié)果均為0.3。然后在不改變模型參數(shù)的情況下對其他幾個不同分辨率DEM進(jìn)行了應(yīng)用,結(jié)果表明,完全根據(jù)30″DEM對應(yīng)的參數(shù)模擬呈村流域的洪水,對于3″、12″與21″DEM而言,峰值會減小,峰現(xiàn)時間也會滯后。分辨率越高,即柵格單元越小,洪水過程的推移坦化現(xiàn)象越明顯。分析其原因,可能主要是由于匯流參數(shù)ke取值未作改變所造成的。ke反映了不同水源在柵格單元上的匯流時間,隨著DEM分辨率的提高,計算單元尺度減小,相應(yīng)的匯流時間應(yīng)有所減小,因此需要對ke進(jìn)行重新率定。通過對ke重新率定后發(fā)現(xiàn),若以30″DEM對應(yīng)的ke率定結(jié)果為準(zhǔn),不同分辨率對應(yīng)的k′e≈ke/n(n為DEM分辨率變化的倍數(shù)),即可以根據(jù)計算單元的尺度對該參數(shù)進(jìn)行校正。表1為只率定ke而不改變其他參數(shù)情況下,不同分辨率DEM所對應(yīng)的柵格型新安江模型模擬結(jié)果的合格率、平均相對誤差水平以及確定性系數(shù)統(tǒng)計,各精度評定指標(biāo)均根據(jù)GB/T 22482—2008《水文情報預(yù)報規(guī)范》計算得到。由表1可以看出,通過對匯流參數(shù)ke進(jìn)行尺度校正,模型均能取得較高的模擬精度,表明流域匯流速度與總的匯流時間不受計算單元尺度的影響。

表1 不同分辨率DEM對應(yīng)的柵格型新安江模型模擬結(jié)果誤差統(tǒng)計

表1和圖2的應(yīng)用結(jié)果表明,雖然利用高分辨率DEM能夠獲得高精度的數(shù)字流域結(jié)果,但并不一定能夠獲得更高精度的分布式水文模擬結(jié)果。由表1可以看出,在3″、12″、21″與30″這4種DEM分辨率情況下,柵格型新安江模型的模擬精度相當(dāng),采用30″DEM不僅能夠提高模型的運行時效,也能取得較好的應(yīng)用效果。此外,模型在進(jìn)行柵格單元產(chǎn)流計算時,其參數(shù)可以根據(jù)每個單元內(nèi)的土壤類型、植被覆蓋等下墊面特征進(jìn)行提取,在分析DEM分辨率對模型影響時,所用的土壤及植被類型數(shù)據(jù)并沒有改變,因此提取的參數(shù)分布變化不大。以WM與SM兩個參數(shù)為例,圖3為不同分辨率DEM對應(yīng)的蓄水容量分布曲線(橫坐標(biāo)為不同的WM與SM值所占的流域面積比例),該曲線是利用式(1)和式(2)計算每個單元對應(yīng)的WM與SM后統(tǒng)計出的兩者在流域上總的分布情況。圖3表明不同分辨率DEM對應(yīng)的張力水和自由水蓄水容量曲線相似度較高,總體趨勢基本一致。由此可以看出,采用客觀合理的參數(shù)估計方法能夠在一定程度上降低參數(shù)對計算單元尺度的依賴性。

圖3 不同分辨率DEM對應(yīng)的蓄水容量分布曲線

4 結(jié) 語

隨著計算機、地理信息系統(tǒng)、遙感等新興技術(shù)以及水文學(xué)理論的進(jìn)一步發(fā)展,基于DEM的分布式水文模型成為水文模型新的發(fā)展方向。本文采用的柵格型新安江模型,即是以DEM為基礎(chǔ)構(gòu)建的松散結(jié)構(gòu)分布式模型,該模型可以較好地考慮降雨和地形地貌、土壤、植被等下墊面條件的空間變異性。

考慮到在實際應(yīng)用與研究中,DEM分辨率對分布式模型有著重要的影響,以皖南山區(qū)呈村流域為例,針對3″、12″、21″與30″這4種分辨率的DEM,探討了分辨率對柵格型新安江模型的影響。研究結(jié)果表明,在不同分辨率情況下進(jìn)行呈村流域的洪水模擬,通過對匯流參數(shù)ke進(jìn)行尺度校正,模型均可以取得較高的模擬精度,利用土壤、植被等下墊面特征估算模型參數(shù)有利于降低參數(shù)的尺度依賴性;此外, DEM分辨率越高,流域數(shù)字化的精度相對越高,但模型的應(yīng)用效果并不一定越好。一般而言,分辨率過高,會導(dǎo)致分布式模型的計算負(fù)擔(dān)過重,在一定程度上會降低其實用性。對于中尺度流域而言,采用30″DEM不僅能夠提高模型的運行時效,也能取得較好的應(yīng)用效果。

柵格型新安江模型在實際應(yīng)用時,只考慮了降雨與下墊面特征在不同計算單元之間的空間變異性,而在每個單元內(nèi),則認(rèn)為各個要素空間分布均勻,如何更好地考慮其在計算單元內(nèi)的空間變異性,還有待更深入的研究。

[1]芮孝芳,黃國如.分布式水文模型的現(xiàn)狀與未來[J].水利水電科技進(jìn)展,2004,24(2):55-58.(RUI Xiaofang, HUANG Guoru.The present and future of the distributed hydrological models[J].Advances in Science and Technology of Water Resources,2004,24(2):55-58.(in Chinese))

[2]熊立華,郭生練.分布式流域水文模型[M].北京:中國水利水電出版社,2004.

[3]芮孝芳,石朋.數(shù)字水文學(xué)的萌芽及前景[J].水利水電科技進(jìn)展,2004,24(6):55-58.(RUI Xiaofang,SHI Peng.Seeds and prospects of digital hydrology[J]. Advances in Science and Technology of Water Resources, 2004,24(6):55-58.(in Chinese))

[4]VIEUX B.Distributed hydrologic modeling using GIS [M].Dordrecht:Kluwer Academic Publishers,2001.

[5]BEVEN K J.Rainfall-runoff modelling:the primer[M]. West Sussex,UK:John Wiley＆Sons,2012.

[6]李致家,姚成,汪中華.基于柵格的新安江模型的構(gòu)建和應(yīng)用[J].河海大學(xué)學(xué)報:自然科學(xué)版,2007,35(2): 131-134.(LI Zhijia,YAO Cheng,WANG Zhonghua. Development and application of grid-based Xinanjiang model[J].Journal of Hohai University:Nature Sciences, 2007,35(2):131-134.(in Chinese))

[7]李致家,姚成,章玉霞,等.柵格型新安江模型的研究[J].水力發(fā)電學(xué)報,2009,28(2):25-33.(LI Zhijia, YAO Cheng,ZHANG Yuxia,et al.Study on gird-based Xinanjiangmodel[J].JournalofHydroelectric Engineering,2009,28(2):25-33.(in Chinese))

[8]YAO Cheng,LI Zhijia,BAO Hongjun,et al.Application ofadevelopedGrid-XinanjiangmodeltoChinese watersheds for flood forecasting purpose[J].Journal of Hydrologic Engineering,2009,14(9):923-93.

[9]姚成,章玉霞,李致家.擴散波與馬斯京根法在柵格匯流演算中的應(yīng)用比較[J].河海大學(xué)學(xué)報:自然科學(xué)版, 2013,41(1):6-10.(YAO Cheng,ZHANG Yuxia,LI Zhijia.Application and comparison of cell-to-cell diffusion wave and Muskingum routing methods[J].Journal of Hohai University:Nature Sciences,2013,41(1):6-10.(in Chinese))

[10]芮孝芳,凌哲,劉寧寧,等.新安江模型的起源及對其進(jìn)一步發(fā)展的建議[J].水利水電科技進(jìn)展,2012,32 (4):1-5.(RUI Xiaofang,LIN Zhe,LIU Ningning,et al. Origin of Xinanjiang model and its further development [J].Advances in Science and Technology of Water Resources,2012,32(4):1-5.(in Chinese))[11]姚成,紀(jì)益秋,李致家,等.柵格型新安江模型的參數(shù)估計及應(yīng)用[J].河海大學(xué)學(xué)報:自然科學(xué)版,2012,40 (1):42-47.(YAO Cheng,JI Yiqiu,LI Zhijia,et al. Parameter estimation and application of grid-based Xin' anjiang model[J].Journal of Hohai University:Natural Sciences,2012,40(1):42-47.(in Chinese))

[12]YAO Cheng,LI Zhijia,YU Zhongbo,et al.A priori parameterestimatesforadistributed,grid-based Xinanjiang model using geographically based information [J].Journal of Hydrology,2012,468/469:47-62.

[13]ANDERSON R,KOREN V,REED S.Using SSURGO data to improve Sacramento Model a priori parameter estimates [J].Journal of Hydrology,2006,320:103-106.

[14]趙人俊.流域水文模擬[M].北京:水利電力出版社, 1984.

[15]趙人俊,王佩蘭.新安江模型參數(shù)的分析[J].水文, 1988(6):2-9.(ZHAO Renjun,WANG Peilan.Analysis of the parameters of the Xinanjiang model[J].Journal of China Hydrology,1988(6):2-9.(in Chinese))

[16]KOREN V,SMITH M,WANG D,et al.Use of soil property data in the derivation of conceptual rainfall-runoff model parameters[C]//Proceedings of the 15th Conference on Hydrology.Long Beach:AMS,2000:103-106.

[17]中國科學(xué)院計算機網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)服務(wù)平臺.SRTM 90米分辨率數(shù)字高程原始數(shù)據(jù)[EB/OL]. [2012-09-28].http://datamirror.csdb.cn.

Effects of DEM resolution on distributed hydrologic simulation

//YAO Cheng,LI Zhijia,ZHANG Yuxia(College of Hydrology and Water Resources,Hohai University,Nanjing210098,China)

Based on the digital elevation model(DEM),the Grid-Xin'anjiang model was employed to simulate flood events at different DEM resolutions to assess the effects of resolution on distributed hydrologic simulations.The results indicate that the simulated results may not be improved by using higher resolution DEM data while the accuracy of the watershed delineation can be enhanced.The Grid-Xin'anjiang model is able to produce comparable and promising results for different DEM resolutions with a recalibration of the routing parameters.For medium-sized watersheds,improved computational efficiency and simulation accuracy can be achieved by using the 30 arc-second DEM.

digital elevation model;grid-Xin'anjiang model;hydrologic simulation;resolution;parameter estimation

10.3880/j.issn.10067647.2013.05.003

P333.2

A

10067647(2013)05001104

20121128 編輯:熊水斌)

國家自然科學(xué)基金(41101017,41130639,51179045);公益性行業(yè)(氣象)科研專項(GYHY201006037)

姚成(1982—),男,江蘇揚州人,講師,博士,主要從事水文預(yù)報研究。E-mail:yaocheng@hhu.edu.cn