李大洋，梁忠民，周 艷

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

水文尺度問題在全球性問題如氣候變暖、水資源短缺等方面扮演著重要的角色，這是水文科學(xué)研究的重點和難點。水文特征展現(xiàn)出的規(guī)律對時空尺度變化的依賴性稱之為水文尺度效應(yīng)，國內(nèi)外學(xué)者對此進(jìn)行了廣泛而深入的研究。Vivoni等[1]在俄克拉何馬州東北部64 km2應(yīng)用分布式水文模型tRIBS研究了不同三角形地形分辨率對水文響應(yīng)的差異，并發(fā)現(xiàn)水文模型對各分辨率范圍內(nèi)的依賴強(qiáng)弱關(guān)系。吳險峰等[2]對100 m×100 m到1 000 m×1 000 m范圍內(nèi)6種水平分辨率的柵格DEM提取流域特征參數(shù)，發(fā)現(xiàn)DEM分辨率的改變對河網(wǎng)精度有影響，流域面積、長度等流域特征參數(shù)差別不大，但對坡度的影響較為強(qiáng)烈。林峰和陳桂芳[3]應(yīng)用分布式水文模型HEC-HMS對時間尺度效應(yīng)進(jìn)行了研究，討論不同模擬步長對模擬結(jié)果的影響程度，發(fā)現(xiàn)模擬步長增加降雨強(qiáng)度發(fā)生較為明顯的變化，洪峰流量與洪量均減小，峰現(xiàn)時間預(yù)報后延。葉許春和張奇[4]分析了不同網(wǎng)格劃分大小對分布式水文模型模擬結(jié)果的影響，并且分為靜態(tài)參數(shù)與變參數(shù)兩種情況；靜態(tài)參數(shù)法指不同網(wǎng)格均使用一套率定好的參數(shù)，而變參數(shù)法相對復(fù)雜，不同網(wǎng)格大小下水文模型重新率定參數(shù)計算結(jié)果，兩種方法比較發(fā)現(xiàn)網(wǎng)格精細(xì)化不一定能夠提高模擬效果，在一定輸入條件下分布式水文模型存在一個最佳的網(wǎng)格尺度。

本文選取淮河黃泥莊流域1983年～2009年間的21場洪水，對分布式水文模型MIKE SHE進(jìn)行模擬驗證，并研究了模型的尺度效應(yīng)。以1 000 m空間網(wǎng)格和1 h時間尺度確定的模型參數(shù)為基礎(chǔ)，分析了不同時空尺度變化對洪水各要素模擬結(jié)果的影響，為選擇合適的水文時空尺度提供一定的參考。

1 流域概況

本研究以史河上游黃泥莊流域作為研究區(qū)域。史河發(fā)源于大別山北脈梨樹尖，由白沙河、根河在黃泥莊上游匯合后，在梅山水庫區(qū)域與麻河、白水河、雙河匯合成史河。黃泥莊流域面積為805 km2，地處大別山區(qū)，平均海拔479 m，雨量豐沛，多年平均降雨量為1 405.3 mm，年平均徑流深738 mm。流域內(nèi)受人類活動影響相對較小，植被良好，土壤滲透性良好，有豐富的地下水。大別山區(qū)是淮河流域暴雨中心之一，對于淮河干流洪峰的形成起到重要作用，因此研究黃泥莊流域的洪水過程有著重要意義。黃泥莊水文站為黃泥莊流域的控制站，流域內(nèi)雨量站點較為豐富，為分布式模型計算提供了條件(見圖1)。

圖1 黃泥莊流域及其鄰近水文站點、水系、土壤質(zhì)地分布及土地利用分布

2 模型與方法

2.1 MIKE SHE模型簡介

MIKE SHE是丹麥水力研究所DHI開發(fā)的基于物理機(jī)制的分布式流域水文模型。MIKE SHE因其強(qiáng)大的水文模擬功能，友好的操作界面得到眾多水文學(xué)者的青睞，并被世界范圍內(nèi)的研究機(jī)構(gòu)、高校所使用。MIKE SHE可以對多種時空尺度、多種地理環(huán)境與氣候條件進(jìn)行模擬，核心模塊包括蒸散發(fā)模塊(ET)、非飽和帶水分模塊(USZ)、飽和地下水模塊(SZ)、坡面流(OL)與河道水流模塊(HD)。各個模塊提供了不同的模擬方法，使用較為靈活，可根據(jù)資料條件使用不同的方法。模型平面上將流域分成若干網(wǎng)格，垂向上將網(wǎng)格分成若干層，也可將流域劃分為若干子流域；蒸散發(fā)模塊有兩種可供選擇，分別是Kristensen-Jensen模型和雙層UZ/ET模型；非飽和帶主要采用帶源匯項的一維Richards方程；飽和帶水流采用三維Darcy方程或者線性水庫法；坡面流和河道水流模塊均采用Saint Venant方程組進(jìn)行差分求解。

根據(jù)黃泥莊流域特點與資料現(xiàn)狀，蒸散發(fā)模塊選用Kristensen-Jensen模塊，非飽和帶模塊采用一維Richards方程近似求解，飽和地下水模塊選用線性水庫法計算，坡面流與河道匯流均為Saint Venant方程組近似求解，從而構(gòu)建了其水文模擬模型MIKE SHE(見圖2)。

圖2 黃泥莊流域水文模擬模型MIKE SHE流程

2.2 資料輸入與模型設(shè)置

一般而言，分布式模型的計算網(wǎng)格大小可設(shè)置為1 km，所以本文MIKE SHE的網(wǎng)格設(shè)置為1 km。時間步長設(shè)置為1 h。地理坐標(biāo)系選用WGS_84坐標(biāo)系，根據(jù)黃泥莊水文站位置中央子午線，使用3度投影分帶，投影坐標(biāo)系為WGS_1984_UTM_ZONE_50N，使用ARCGIS軟件對相關(guān)空間數(shù)據(jù)進(jìn)行統(tǒng)一處理。MIKE SHE需要的主要基礎(chǔ)資料如下：①網(wǎng)格大小為90 m的DEM數(shù)據(jù)(SRTMDEM，下載網(wǎng)址為http:∥www.gscloud.cn/)；②1∶10萬土地利用資料數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(下載網(wǎng)址為http:∥www.resdc.cn)； ③1∶100萬土壤數(shù)據(jù)來自世界土壤數(shù)據(jù)庫(HWSD)的中國數(shù)據(jù)集V1.1(下載網(wǎng)址為http:∥westdc.westgis.ac.cn)；④由于黃泥莊流域無國家氣象觀測站，所以使用臨近霍山站的氣象數(shù)據(jù)如日均氣溫、日最大/小氣溫、日平均風(fēng)速、日平均相對濕度、日照時數(shù)、日平均大氣壓等；⑤降雨數(shù)據(jù)以及流量數(shù)據(jù)來自雨量站以及水文站(見圖1)。

2.3 模型參數(shù)的確定

分布式模型與集總式模型一處明顯不同就是很多參數(shù)具有明確的物理意義。因此，本文模型參數(shù)采用物理意義及模型率定這兩種方式確定。

2.3.1通過物理意義確定參數(shù)

(1)參考蒸散發(fā)。它是指代表短綠色作物完全遮蔽地面，具有統(tǒng)一高度以及土壤中具有充足水分時的蒸散發(fā)率，它反映了可用于蒸發(fā)水的能量以及可用于運(yùn)輸水蒸氣的風(fēng)從地面進(jìn)入底層大氣的狀態(tài)，一般以苜蓿作為參照作物。MIKE SHE進(jìn)行場次洪水模擬需要的蒸散發(fā)資料為日參考蒸散發(fā)數(shù)據(jù)，其推薦使用FAO Penman-Monteith公式。本文采用霍山站1983年～2009年的氣象資料進(jìn)行參考蒸散發(fā)的計算(見圖3)。

圖3 1983年～2009年間日參考蒸散發(fā)過程

(2)葉面積指數(shù)與根系深度。葉面積指數(shù)LAI是1947年英國生態(tài)學(xué)家Watson[5]提出的用于反應(yīng)植物葉面數(shù)量變化的量。根據(jù)Watson的描述LAI定義為每單位地表面積上的綠葉面積總和或者單位面積上所有葉片向下投影后的面積總和，其為一個無量綱的量[6]。LAI一般在0～10之間，植被越濃密，LAI越大；反之，LAI越小。LAI與根系深度均為MIKE SHE計算實際蒸散發(fā)的參數(shù)，本文根據(jù)黃泥莊流域土地利用類型，并查詢MIKE SHE自帶的土地類型資料庫(MIKE SHE 2014版本)，從而確定LAI與根系深度的值，黃泥莊流域LAI與根系深度的分布圖(圖略)。

(3)曼寧系數(shù)(糙率系數(shù))。在水力學(xué)中，反映地表粗糙程度的量稱為糙率，用曼寧系數(shù)表示，即Manning n或者M(jìn)anning M(Manning n的倒數(shù))。根據(jù)黃泥莊流域土地利用類型，參考Engman等[7]推薦的不同地表條件下的糙率大小，得出曼寧系數(shù)Manning M流域分布(用于坡面流的計算)(圖略)。

(4)土壤水文常數(shù)與水力傳導(dǎo)度。土壤水分是水量平衡計算的重要組成部分，是植物獲取水分的主要來源[8]。土壤中水分受到的作用力有分子力、毛管力和重力，根據(jù)受力的不同，可以將土壤水分分為吸濕水、薄膜水、毛管水、重力水等四類，這四類土壤水分均有一個臨界值，土壤水文常數(shù)就是用來描述土壤水分閾值效應(yīng)的參量。MIKE SHE中計算非飽和帶土壤水分流動需要用到的土壤水文常數(shù)為凋萎系數(shù)、田間持水量與飽和含水量以及水力學(xué)參數(shù)飽和水力傳導(dǎo)度。參考美國農(nóng)業(yè)部土壤質(zhì)地三角形及其飽和水力傳導(dǎo)度圖，根據(jù)黃泥莊流域土壤類型(見圖1b)可以查到各土壤水文常數(shù)(容積含水率表示)與飽和水力傳導(dǎo)度(見表1)。

表1 黃泥莊流域不同土壤類型土壤水文常數(shù)與飽和水力傳導(dǎo)度

2.3.2通過模型率定確定參數(shù)

MIKE提供了適用于各模塊參數(shù)率定的工具AUTOCAL，并有兩種參數(shù)率定方法SCE-UA[9]與PSE[10]可供選擇。本文使用SCE-UA算法進(jìn)行參數(shù)自動率定。SCE-UA是一種全局優(yōu)化算法，融合了一些算法的優(yōu)勢，在一定程度上能夠解決水文模型高維參數(shù)全局搜索問題[11]。確定性系數(shù)DC可用于評價水文模型模擬與實測流量過程的吻合程度，大小在-∞到1之間，值越大表明模擬與實測過程越接近，等于1表明模擬值與實測值相同[12]。本文以1-DC作為SCE-UA的目標(biāo)函數(shù)。即

(1)

表2為SCE-UA各參數(shù)的取值范圍以及參數(shù)意義，表3給出了本文需要率定的模型參數(shù)。因此，SCE-UA各參數(shù)可取n=8，m=17，q=9，β=17；p取值較為靈活，與率定參數(shù)個數(shù)無關(guān)，值越大計算精度越高但計算時間越長，Duan[13]推薦p取2對于一般水文模型即可滿足精度要求。以黃泥莊水文站1983年～2009年間21場場次洪水的前15場作為率定期，后6場作為驗證期。參數(shù)率定結(jié)果如表3所示。

表2 SCE-UA參數(shù)

表3 黃泥莊流域MIKE SHE/11參數(shù)率定結(jié)果

表4 黃泥莊水文站1983年～2009年間場次洪水模擬結(jié)果

3 模型模擬效果及尺度效應(yīng)分析

3.1 模型率定與驗證

GB/T 22482—2008《水文情報預(yù)報規(guī)范》推薦采用絕對誤差、相對誤差、確定性系數(shù)DC三種指標(biāo)對洪峰流量、洪峰水位、洪峰出現(xiàn)時間、洪水過程等洪水要素的模擬結(jié)果進(jìn)行評價。本文以洪峰預(yù)報相對誤差、洪峰出現(xiàn)時間(簡稱“峰現(xiàn)時間”)誤差、確定性系數(shù)作為模型計算精度評定指標(biāo)。洪峰預(yù)報許可誤差為20%，峰現(xiàn)時間誤差為3 h，根據(jù)經(jīng)驗可將確定性系數(shù)0.7以上作為許可范圍。模擬洪水場次滿足以上許可條件稱之為合格場次。

黃泥莊水文站1983年～2009年間的21場洪水模擬結(jié)果見表4。率定期15場洪水有10場合格，平均合格率67%，驗證期6場洪水其中5場合格，平均合格率83%，全部洪水場次合格率為71%，達(dá)到乙級評價標(biāo)準(zhǔn)。由驗證期各場次洪水模擬流量與實測流量過程(見圖4)可以看出，黃泥莊水文站模擬流量過程線與實測流量過程線吻合度較高。

3.2 尺度效應(yīng)分析

在上述率定好的參數(shù)的基礎(chǔ)上，分別研究5種網(wǎng)格尺度以及5種時間尺度變化對分布式水文模型MIKE SHE模擬洪水各要素如洪峰、洪量、峰現(xiàn)時間等的影響。為節(jié)省篇幅，本文以20080622號洪水為例進(jìn)行說明。

圖4 黃泥莊水文站驗證期各洪水場次實測流量與模擬流量過程

3.2.1空間尺度效應(yīng)

選取大小為250 、500 、1 000 、2 000 、3 000 m的5種網(wǎng)格尺度，時間步長固定為1 h，以20080622號洪水為例進(jìn)行模擬。從20080622號場次洪水各網(wǎng)格尺度下的模擬洪水過程(見圖5)可以看出，洪水過程線隨著網(wǎng)格尺度的增加而變高，然而當(dāng)網(wǎng)格達(dá)到3 000 m時洪水過程線開始變的低矮；網(wǎng)格尺度在1 000 m以內(nèi)時洪水過程線隨網(wǎng)格尺度變動幅度較大，網(wǎng)格大于1 000 m時變化較小。從各洪水要素來看(見表5)，峰現(xiàn)時間未發(fā)生變化，當(dāng)網(wǎng)格≤2 000 m時，洪峰與洪量隨著網(wǎng)格尺度的增加而增加，如網(wǎng)格為250 m洪峰、洪量(徑流深)分別為915 m3/s和27.48 mm；網(wǎng)格為2 000 m時洪峰、洪量為1 900 m3/s和50.35 mm；網(wǎng)格尺度為3 000 m時的洪峰流量與洪量介于網(wǎng)格1000 m與500 m之間，且更靠近網(wǎng)格尺度為1 000 m時的洪流流量與洪量。從模擬精度來看(見表5)，網(wǎng)格≥1 000 m時洪峰誤差、洪量誤差、峰現(xiàn)時間誤差均符合精度要求，而網(wǎng)格<1 000 m時除峰現(xiàn)時間誤差外其余指標(biāo)均不符合精度要求。

圖5 黃泥莊流域20080622號場次洪水各網(wǎng)格尺度下模擬流量與實測流量過程

表5 黃泥莊流域20080622號場次洪水在各網(wǎng)格尺度下的洪水要素

網(wǎng)格大小/m洪峰/m3·s-1洪峰誤差/%洪量/mm洪量誤差/%峰現(xiàn)時間誤差/h250915-48.60 27.48 -48.16 05001380-22.47 38.73 -26.93 01 00018101.69 48.44 -8.61 02 00019006.74 50.35 -5.01 03 0001740-2.25 48.30 -8.88 0

3.2.2時間尺度效應(yīng)

分別將MIKE SHE模擬的時間步長設(shè)置為0.1、0.5、1、2、3 h，網(wǎng)格大小均為1 000 m，并將1 h時間間隔的各個雨量站的降雨數(shù)據(jù)與參考蒸散發(fā)數(shù)據(jù)處理成相應(yīng)時段，以20080622號洪水為例進(jìn)行模擬。如圖6所示，洪水過程線隨著時間步長的增加而變的平坦。從洪水要素來看(見表6)，洪峰與洪量均隨著時間步長的增加而減小，且洪峰減小幅度較為明顯，峰現(xiàn)時間未呈現(xiàn)趨勢性變化。從計算精度來看(見表6)，除3 h步長所模擬的洪峰誤差較差外，其余均符合精度要求。

圖6 黃泥莊流域20080622號場次洪水各時間步長下的模擬流量與實測流量過程

表6 黃泥莊流域20080622號場次洪水各時間步長下的洪水要素

4 主要結(jié)論

本文使用分布式水文模型MIKE SHE對淮河黃泥莊流域1983年～2009年間的場次洪水進(jìn)行了模擬與驗證，并研究了模型的時空尺度效應(yīng)，主要結(jié)論如下：

(1)21場洪水有15場模擬結(jié)果合格，合格率為71%，達(dá)到乙級評價標(biāo)準(zhǔn)；這表明MIKE SHE模型對研究流域的洪水模擬預(yù)報具有一定的適用性。

(2)洪水模擬結(jié)果對網(wǎng)格與步長的變化較為敏感且主要體現(xiàn)在洪峰和洪量上。在一定網(wǎng)格尺度范圍內(nèi)(250～2 000 m)，洪峰與洪量隨網(wǎng)格尺度增加而增大，但網(wǎng)格達(dá)到3 000 m時，洪峰與洪量開始變??；洪峰與洪量隨時間步長(0.1 ～3 h)的增加而增大。受場次洪水模擬時段以及流域特征影響，網(wǎng)格與時間步長變化均未對峰現(xiàn)時間產(chǎn)生明顯影響。