羅玉玲，張?zhí)鹛?，章玉霞，杜若?/p>

(1.河海大學水文水資源學院，江蘇 南京 210098；2.江蘇省徐州市水務局，江蘇 徐州 221000；3.南京河海大學出版社有限公司，江蘇 南京 210098)

無資料地區(qū)的洪水預報是目前水文工作者面臨的挑戰(zhàn)之一。2003年國際水文學會提出了水文資料匱乏流域的洪水預報10年計劃，期間涌現(xiàn)了大量相關(guān)研究成果[1- 2]。其中流域水文模型作為研究水文現(xiàn)象的有效工具是解決無資料流域水文預報的代表性方法之一，是無資料地區(qū)水文洪水預報的重要解決手段[3]。分布式水文模型因其對降雨和下墊面條件的空間變異性的考慮，可以更好地利用現(xiàn)代化技術(shù)條件，是目前水文模型的重要研究方向之一[4- 5]。匯流過程是洪水預報的重要環(huán)節(jié)，影響最終預報的精度，是水文模型研究的重點之一[6]。

目前的分布式水文模型研究，在匯流方法上，主要是基于圣維南方程組及其簡化算法的水力學方法如運動波、擴散波等，還有基于水量平衡和槽蓄方程的水文學方法，如線性水庫法、馬斯京根法[7- 8]。水力學方法具有較強的物理性，雖然可以有效反映匯流過程中的諸如河道斷面變化等的影響，但需要較多的河道實測資料，且因其計算復雜，目前尚沒有方法求得其精確解析解，在無資料地區(qū)運用受限。水文學方法是水力學的簡化，計算更加簡單，資料需求也相對較少，但在無資料地區(qū)的運用仍然受到一定限制[9]。為更好地利用下墊面條件，對無資料地區(qū)進行快速洪水預報，本文引入網(wǎng)格水滴匯流方法[10]，進行分布式水文模型的匯流演算，以更綜合地考慮下墊面條件和水文特征[11]。

網(wǎng)格新安江模型是分布式水文模型，其自創(chuàng)建以來就有著較為廣泛的應用[12- 17]。故本文在網(wǎng)格新安江基礎(chǔ)上，同時構(gòu)建采用網(wǎng)格水滴匯流網(wǎng)格新安江模型，將其應用于濕潤地區(qū)屯溪流域，與原有基于擴散波方法的網(wǎng)格新安江模型及基于馬斯京根法的網(wǎng)格新安江模型進行模擬精度對比[18]，驗證其可行性，分析網(wǎng)格新安江模型應用網(wǎng)格水滴匯流方法的合理性；并將基于網(wǎng)格水滴的網(wǎng)格新安江模型應用于半濕潤半干旱地區(qū)東灣流域，驗證其適用性。

1 模型與方法

1.1 網(wǎng)格水滴匯流方法

“網(wǎng)格水滴”是流域上將有降雨并發(fā)生產(chǎn)流的網(wǎng)格看作水滴，水滴受到重力和阻力的影響，沿坡度方向運動，其到達出口斷面的時間主要取決于匯流路徑及匯流速度[19]。

1.1.1坡面匯流

根據(jù)D8匯流方法確定網(wǎng)格水滴匯流的時間和路徑[20]，坡面匯流速度公式采用由美國水土保持局工程手冊提出的流速與地表坡度的函數(shù)關(guān)系式[21]：

(1)

式中，vs—坡面流水流速度，m/s；a—反應地表糙率的經(jīng)驗系數(shù)；Ss—沿水流方向的地表坡度。

1.1.2河道匯流

河道匯流采用根據(jù)曼寧公式演變的河道匯流速度公式[10]：

(2)

式中，vr—渠道流水流速度，m/s；α—反應糙率的參數(shù)；β—經(jīng)驗指數(shù)；A—網(wǎng)格匯水面積，m2；Sr—河道網(wǎng)格坡度。

1.1.3匯流時間矩陣

根據(jù)匯流路徑和匯流速度，將各個網(wǎng)格單元所需要的時間累加得到其到達流域出口所需要的時間，即為該網(wǎng)格的匯流時間。

(3)

所有網(wǎng)格的匯流時間組合一起，即為流域的匯流時間矩陣。

1.2 模型構(gòu)建

1.2.1基于擴散波的網(wǎng)格新安江模型[18]

基于擴散波的網(wǎng)格新安江模型(DW模型)是一個分布式水文模型，以DEM網(wǎng)格為計算單元，主要分為蒸散發(fā)、蓄滿產(chǎn)流、分水源、匯流4個計算模塊，其中匯流模塊采用一維擴散波匯流方法，依次將各種水源演算至流域出口，該方法物理基礎(chǔ)較好，可以利用下墊面信息估算參數(shù)[6]。

1.2.2基于逐網(wǎng)格馬斯京根法的網(wǎng)格新安江模型[18]

基于逐網(wǎng)格馬斯京根法的網(wǎng)格新安江模型(MK模型)其模塊區(qū)分與DW模型一致，其匯流方法采用的是逐網(wǎng)格馬斯京根法，依次將地表徑流、壤中流、地下徑流演算至流域出口。該方法具有計算簡便、實用性強的特點。

1.2.3基于網(wǎng)格水滴的網(wǎng)格新安江模型

基于網(wǎng)格水滴的網(wǎng)格新安江模型(CW模型)，產(chǎn)流以蓄滿產(chǎn)流原理為基礎(chǔ)，計算每個網(wǎng)格單元的產(chǎn)流量，采用自由水蓄水庫結(jié)構(gòu)，將其劃分為地表徑流、壤中流及地下徑流3種水源，再根據(jù)匯流演算規(guī)則演算至出口斷面[6]。采用的網(wǎng)格水滴模型進行匯流演算時假設柵格單元是由坡面及河道組成，坡面柵格分為地面徑流、地下徑流及壤中流，采用網(wǎng)格水滴匯流方法，根據(jù)網(wǎng)格水滴匯流速度、匯流路徑及時間矩陣，將徑流演算至流域出口。

2 研究流域和數(shù)據(jù)處理

2.1 研究流域

屯溪流域如圖1所示，位于安徽省，該流域靠近中國東南沿海地區(qū)，屬于亞熱帶季風氣候區(qū)、流域內(nèi)四季分明，氣候溫和，多年平均氣溫約為17℃。屯溪流域面積為2670km2，地勢西高東低，相對高差較大，水流方向自西向東。流域內(nèi)植被良好，雨量充沛，多年平均降雨為1800mm，降水在年內(nèi)、年際分配不均勻，汛期內(nèi)的降雨量占全年總降雨量60%；為典型的濕潤流域。

圖1 屯溪流域水系及水文站點圖

東灣流域(如圖2所示)伊河河源地區(qū)，流域面積2623km2。屬于大陸性季風氣候，是半濕潤地區(qū)。降雨量分布不均勻，流域內(nèi)多年平均降水量為773mm，日降水量一般在100mm以上，最大值可達600mm。暴雨主要集中在7—8月，占總降水量的60%，多年平均蒸發(fā)量為564mm，流域地勢西高東低，上游地區(qū)植被良好，主要為林地，下游河道附近裸地較多。

圖2 東灣流域水系及站點圖

2.2 數(shù)據(jù)

本文模型所需下墊面數(shù)據(jù)有DEM數(shù)字高程數(shù)據(jù)、土壤類型和土地利用數(shù)據(jù)。DEM數(shù)據(jù)通過地理空間數(shù)據(jù)云網(wǎng)站下載，采用90m分辨率數(shù)據(jù)。土壤類型數(shù)據(jù)為世界土壤數(shù)據(jù)庫HWSD中1∶100萬土壤柵格數(shù)據(jù)。土地利用數(shù)據(jù)為美國地質(zhì)勘探局網(wǎng)站下載，是分辨率為1km的矢量文件。

將研究流域數(shù)據(jù)包括降雨、蒸發(fā)及流量樞機們進行1h插值。選取屯溪流域內(nèi)2008—2016年23場洪水資料進行模擬，其中17場洪水進行率定，7場洪水進行驗證；選取東灣流域1996—2011年14場洪水資料進行模擬，其中10場洪水進行率定，4場洪水進行驗證。屯溪流域流量數(shù)據(jù)來自屯溪水文站，降雨資料來自月譚等17個雨量站(如圖1所示)；東灣流域流量數(shù)據(jù)來自東灣水文站，降雨資料來自譚頭等8個雨量站(如圖2所示)。

3 結(jié)果與討論

3.1 參數(shù)率定

在應用模型進行屯溪流域洪水模擬時，采用的柵格單元分辨率為1km×1km(30″×30″)，網(wǎng)格水滴匯流參數(shù)主要參數(shù)為a及α。其中2個流域的α率定為0.0017，屯溪流域的a值率定為0.84m/s，東灣流域的a值率定為0.32m/s。

坡面匯流公式及河道匯流公式與流域下墊面輸入極其相關(guān)，由于在處理DEM填洼時，對于洼地的處理為將其處理為與附近網(wǎng)格DEM值相當，如此會造成單元網(wǎng)格的坡度過小導致匯流速度過慢影響最終匯流，為解決此問題，采取措施為坡度小于0.1%時，網(wǎng)格采用流域坡面網(wǎng)格的平均坡度進行速度計算。

3.2 時間矩陣分析

根據(jù)CW模型，可以得到屯溪流域時間矩陣。如圖3所示，從屯溪流域的時間分布圖來看，匯流時間的分布為河道網(wǎng)格的匯流時間小于坡面匯流時間，且相對河道距離近的地方，匯流時間更短(如圖4所示)。造成如此匯流時間分布的主要原因是，河道的匯流速度比于坡面匯流速度大，水滴在河道匯流所花費時間更短，坡面水滴在匯流過程受到坡度及匯流路徑的影響大，部分坡度較小網(wǎng)格，匯流速度較慢，且因其遠離河道，坡面匯流路徑長，導致匯流時間相應變長。隨著坡面水滴流入河道，匯流速度也相應增大，故距離河道近的地方匯流時間相對較短。

圖3 屯溪流域匯流時間分布圖

圖4 匯流時間與水系分布圖

3.3 模擬結(jié)果分析對比

依據(jù)GB/T 22482—2008《水文情報預報規(guī)范》[22]，選取洪峰相對誤差、徑流深相對誤差、峰現(xiàn)時間誤差和確定性系數(shù)4個指標對模擬效果進行綜合評價。

3.3.1屯溪流域結(jié)果對比

根據(jù)屯溪水文站實測流量數(shù)據(jù)，采用DW模型、MK模型及CW模型3種不同匯流模型模擬屯溪流域洪水過程的合格率、確定性系數(shù)均值的統(tǒng)計結(jié)果見表1，圖5是3種模型洪水模擬結(jié)果與實測流量的對比圖。

表1 屯溪流域各模型模擬結(jié)果統(tǒng)計

率定期共有17場洪水，率定期的徑流深模擬，DW模型、MK模型和CW模型模擬相對誤差分別介于-27.6%～15.85%、23.16%～19.93%及-23.21%～18.27%之間。DW模型、MK模型和CW模型的洪峰相對誤差分別介于-32.5%～17.9%、-11.63%～27.27%和-16.18%～22.41%；DW模型的峰現(xiàn)時間誤差分別在-4～20h，-2～7h及-2～7h之間。DW模型、MK模型和CW模型的確定性系數(shù)分別介于0.81～0.98、0.81～0.97和0.8～0.97之間。

驗證期共有7場洪水，驗證期的的徑流深模擬，DW模型、MK模型及CW模型的相對誤差介于10.8%～15.2%、-10.43%～25.45%和-14.56%～18.35%；DW模型、MK模型和CW模型的洪峰相對誤差分別介于-17.3%～6%、-14.25%～0.58%和-17.82%～0.73%之間；3種匯流方法的時間合格率都為85.7%，誤差介于-2～4h、-2～4h及-3～4h之間。DW模型、MK模型和CW模型的確定性系數(shù)范圍在0.9～0.96、0.89～0.97及0.91～0.98之間

圖5 屯溪流域20080608號和20160420號洪水實測與模擬洪水過程線對比

根據(jù)3種模型的模擬結(jié)果及圖5三種模型模擬洪水過程與實測流量對比圖，表明3種匯流方法在屯溪流域洪水過程的模擬精度相當，從徑流深合格率、洪峰合格率和平均確定性系數(shù)上看均達到甲級預報精度，采用這3種方法進行網(wǎng)格新安江的匯流計算，均可以取得較好的應用效果，但20080608號洪水，3種匯流方法的峰現(xiàn)時間都相對提前，MK方法的峰現(xiàn)時間較為接近實際峰現(xiàn)時間，DW方法峰現(xiàn)時間較為提前一些，CW方法居于中間，3種匯流方法精度沒有較大區(qū)別，都有較好的模擬效果。3種匯流方法都能很好的反應洪水的走勢，對于降雨的反應也較好，能夠很好的應用于濕潤地區(qū)。在物理意義上，DW模型的參數(shù)具有較為明確的物理意義，對于下墊面特征信息有著良好的考慮；MK模型是水文學方法，對下墊面情況的考慮較少，CW模型物理意義較為明確，且對下墊面特征信息考慮較多；在計算效率方面，DW模型計算一次所花費時間比MK模型及CW模型長，計算效率較低，MK模型的計算效率與CW模型的計算效率相當，均具有較高的計算效率。隨著數(shù)據(jù)信息的增多，如分辨率的提高，DW模型與CW模型及K模型計算時間的差異增大。DW模型物理意義較為明確，適用性較為廣泛，在平原地區(qū)也可以較好應用，對下墊面條件對匯流影響考慮較為充分，在無資料地區(qū)應用具有一定的適用性，但是難以對無資料地區(qū)洪水進行快速預報；MK模型計算簡便，運算效率較高，實用性強，但較為依賴降雨徑流的觀測資料，且在無資料地區(qū)的適用性較弱；CW模型的參數(shù)物理意義較為明確，計算簡便，具有較高的運算效率，且其輸入項可以根據(jù)下墊面特征提取，參數(shù)估算較為簡單，故對于無資料地區(qū)的快速洪水預報具有一定的適用性。

3.3.2東灣流域結(jié)果分析

根據(jù)東灣流域流量數(shù)據(jù)，采用CW模型對東灣流域洪水過程進行模擬，驗證網(wǎng)格水滴匯流模型在半濕潤半干旱地區(qū)的適用性，CW模型在東灣流域的模擬結(jié)果誤差統(tǒng)計見表2。

從表2可以看出，東灣流域用于率定的10場洪水種，7場洪水的徑流深誤差小于20%，率定期平均徑流深合格率為70%；模擬的平均洪峰合格率為90%；平均峰現(xiàn)時間合格率為90%；平均確定性系數(shù)為0.70；達到了乙級預報精度。率定期的平均徑流深合格率為75%，平均洪峰合格率為75%，平均峰現(xiàn)時間合格率為100%，平均確定性系數(shù)為0.70；圖6為東灣流域19980814號洪水的過程，率定期的10場洪水對匯流參數(shù)進行了很好的率定，使得率定期和驗證期次洪模擬達到了較好的精度，結(jié)果表明，網(wǎng)格水滴匯流模型在半濕潤半干旱地區(qū)的平均模擬精度可以達到乙級標準，即網(wǎng)格水滴模型的匯流方法在半濕潤半干旱地區(qū)具有一定的適用性。

表2 東灣流域CW模型模擬結(jié)果誤差統(tǒng)計

圖6 東灣流域19980814號洪水

4 結(jié)語

引入網(wǎng)格水滴匯流(CW模型)的網(wǎng)格新安江模型，將其應用于濕潤地區(qū)屯溪流域及半濕潤半干旱流域東灣流域。研究結(jié)果表明，與采用擴散波方法(DW模型)的網(wǎng)格新安江模型及采用逐網(wǎng)格馬斯京根方法(MK模型)的網(wǎng)格新安江模型相比，CW模型在計算精度相當?shù)那闆r下，可以兼具計算效率、實用性及較為明確的物理意義的特點；CW模型在濕潤半濕潤地區(qū)東灣流域的適用性良好，為解決無資料地區(qū)洪水的快速預報問題提供了一個較好的解決思路。