易 彬,陳 璐

(1.華中科技大學(xué)土木與水利工程學(xué)院,湖北 武漢 430074；2.華中科技大學(xué)數(shù)字流域科學(xué)與技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430074)

流域匯流是地面徑流、壤中流、地下徑流匯集到流域出口斷面的徑流過(guò)程。凈雨滴通過(guò)不同介質(zhì)、不同路徑、不同粗糙面向流域出口斷面移動(dòng),實(shí)際匯流過(guò)程十分復(fù)雜。過(guò)去100 a間,大量學(xué)者從水動(dòng)力、水文、水滴運(yùn)動(dòng)以及系統(tǒng)視角觀察和分析流域匯流過(guò)程[1]。1932年,Sherman[2]在研究如何用降水量推求徑流量的過(guò)程中,提出了單位線的概念,目前該方法已成為使用最為廣泛的匯流計(jì)算方法。

單位線將匯流過(guò)程概化為線性時(shí)不變系統(tǒng),實(shí)際上,匯流特性不僅受控于流域地形地貌和下墊面特征,也與凈雨時(shí)空分布等因素有關(guān)[3]。基于流域地貌參數(shù),Rodríguez-Iturbe等[4]提出了地貌單位線方法,但該類方法無(wú)法考慮流域地形及下墊面特征的空間分布；為了考慮下墊面空間異質(zhì)性對(duì)匯流的影響,Maident等[5]提出了基于數(shù)字高程模型的分布式匯流單位線。分布式單位線計(jì)算方法將水質(zhì)點(diǎn)在流域中的滯留時(shí)間等價(jià)為流域瞬時(shí)單位線,其核心在于速度場(chǎng)的計(jì)算,常用的坡面流速計(jì)算公式有SCS流速公式、曼寧公式、均勻流方程、達(dá)西公式等[6],上述流速公式考慮了植被、坡度、土地利用等因素的影響[7-8],但忽略了由時(shí)變雨強(qiáng)、上游入流貢獻(xiàn)等因素引起的徑流匯流非線性問(wèn)題。

有學(xué)者通過(guò)建立坡面流速與時(shí)段凈雨強(qiáng)度、上游入流貢獻(xiàn)等因素的定量映射關(guān)系,提出了時(shí)變分布式單位線法[9-10]。孔凡哲等[11]基于改進(jìn)SCS流速計(jì)算公式提出了考慮降水強(qiáng)度的時(shí)變分布式單位線計(jì)算方法；Bunster等[12]考慮動(dòng)態(tài)的上游入流貢獻(xiàn),開發(fā)了時(shí)變動(dòng)態(tài)上游入流貢獻(xiàn)的分布式單位線,提高了洪峰預(yù)報(bào)精度。然而,上述方法均假定全流域水質(zhì)點(diǎn)存在靜態(tài)匯流路徑,即流域上均勻分布的凈雨滴通過(guò)唯一的匯流路徑向出口流動(dòng),事實(shí)上,以蓄滿產(chǎn)流為主的南方濕潤(rùn)地區(qū),在流域未達(dá)到蓄滿狀態(tài)之前,只有蓄滿區(qū)域產(chǎn)生凈雨,其匯流路徑應(yīng)隨流域蓄水狀態(tài)而動(dòng)態(tài)變化,當(dāng)流域蓄滿后,動(dòng)態(tài)匯流路徑轉(zhuǎn)變?yōu)殪o態(tài)匯流路徑。

為此,本文考慮下墊面空間分布異質(zhì)性和降水強(qiáng)度的時(shí)變特性,根據(jù)各個(gè)時(shí)段凈雨量級(jí)大小、下墊面土壤含水量狀態(tài)推求流域時(shí)變匯流路徑,嘗試根據(jù)流域產(chǎn)流情況將流域靜態(tài)匯流路徑轉(zhuǎn)換為動(dòng)態(tài)匯流路徑,推求同時(shí)考慮降水強(qiáng)度和下墊面土壤含水量時(shí)空分布特征的坡面流速計(jì)算公式,提出一種考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線。

1 研究方法

現(xiàn)有時(shí)變分布式單位線存在2個(gè)假定：全流域蓄水量在降水結(jié)束前可以達(dá)到飽和狀態(tài)以及全流域的坡面、河網(wǎng)匯流路徑唯一。本文從流速計(jì)算和匯流路徑兩方面提出解決思路,一方面,根據(jù)土壤蓄水狀態(tài)提取流域動(dòng)態(tài)匯流路徑,另一方面,通過(guò)考慮時(shí)變土壤含水量改進(jìn)現(xiàn)有坡面流速計(jì)算公式,在此基礎(chǔ)上,結(jié)合改進(jìn)時(shí)變坡面流速計(jì)算公式和動(dòng)態(tài)匯流路徑推求考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線。

1.1 產(chǎn)流計(jì)算

SCS模型由美國(guó)農(nóng)業(yè)部水土保持局于1954年開發(fā),特點(diǎn)是結(jié)構(gòu)簡(jiǎn)單、參數(shù)較少,被廣泛用于小型集水區(qū)徑流預(yù)報(bào)[13],本文采用該模型進(jìn)行時(shí)段凈雨(Rt)計(jì)算。

1.2 動(dòng)態(tài)匯流路徑

本文將動(dòng)態(tài)匯流路徑定義為蓄滿柵格流向流域出口過(guò)程中所經(jīng)過(guò)的路徑,由于降水過(guò)程的持續(xù),流域蓄水狀態(tài)實(shí)時(shí)變化導(dǎo)致匯流路徑動(dòng)態(tài)變化,獲取任一時(shí)刻流域匯流路徑的關(guān)鍵在于確定流域當(dāng)前蓄水狀態(tài)分布情況,新安江模型常用帕累托分布表征土壤蓄水能力的空間分布不均勻性,其缺陷在于無(wú)法獲得流域蓄水能力在空間上的具體分布情況[14]。研究表明,地形指數(shù)與流域蓄水能力密切相關(guān)[15],本文采用地形指數(shù)(TI)表征流域土壤蓄水能力,根據(jù)地形指數(shù)計(jì)算柵格蓄水能力[16],表達(dá)式如下：

(1)

式中：WM,j為j柵格單元的最大蓄水容量；WM,min為全流域所有柵格中的最小蓄水容量；WM,max為全流域所有柵格中的最大蓄水容量；TI,min=min{TI,j,j=1,2,…,N}，N為流域柵格單元總數(shù)；TI,max=max{TI,j,j=1,2,…,N}；n為待率定參數(shù),反映蓄水容量不均勻分布的系數(shù),該值越大,流域蓄水容量的空間分布越不均勻。

(2)

式中：Fαt為流域蓄滿部分的面積；Pi為i時(shí)段降水量；Ri為i時(shí)段凈雨量；Wj,1為j柵格初始時(shí)段蓄水量；Aj為柵格單元面積。

采用動(dòng)態(tài)匯流路徑進(jìn)行流域匯流計(jì)算時(shí),由于僅蓄滿柵格發(fā)生產(chǎn)流,而SCS模型認(rèn)為全流域均勻產(chǎn)流,因此,t時(shí)段凈雨量(Rt)需重新分配至蓄滿面積,根據(jù)水量平衡原理：

(3)

(4)

1.3 基于動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線

分布式單位線方法已較為成熟,被廣泛用于中小流域匯流計(jì)算,具體流程詳見文獻(xiàn)[11],本文所提方法步驟如下：

(1) 根據(jù)流域DEMs將流域柵格化,確定柵格邊長(zhǎng),在此基礎(chǔ)上,采用D8算法確定柵格水流方向,并根據(jù)柵格水流方向提取該柵格至流域出口的匯流路徑集合Rroad,j(第j個(gè)柵格至流域出口路徑上的柵格單元集合),則全流域所有柵格匯流路徑集合可表示為

Rroad={Rroad,j|j=1,2,…,N}

(5)

由于該路徑集合是靜態(tài)的,然而實(shí)際匯流情況為只有蓄滿區(qū)域有匯流路徑,為此,結(jié)合分布式蓄水容量,推求不同蓄滿比例(αt)下蓄滿柵格的匯流路徑集合(Rroad,αt),滿足：

Rroad,αt={Rroad,j|Aj∈Fαt}

(6)

式中：Rroad,αt為流域蓄滿區(qū)域的匯流路徑集合,Rroad,αt?Rroad。隨著降水過(guò)程的持續(xù),αt不斷變化,因此,Rroad,αt構(gòu)成流域不同蓄水狀態(tài)下的動(dòng)態(tài)匯流路徑集合。

(2) 計(jì)算每個(gè)柵格逐時(shí)段的匯流速度(Vj,αt),獲得不同流域土壤含水量狀態(tài)下的流速分布場(chǎng)?，F(xiàn)有分布式單位線僅考慮降水時(shí)空分布不均勻性和下墊面空間分布異質(zhì)性對(duì)流域匯流速度的影響,其假定流域在降水脈沖結(jié)束之前可達(dá)到蓄滿狀態(tài),這與實(shí)際情況不符,忽略了時(shí)變土壤含水量造成下墊面狀態(tài)的改變,已有研究表明,坡面流速與植被覆蓋、雨強(qiáng)、坡度、土壤含水狀態(tài)有關(guān)[7-8,17]。為此,針對(duì)蓄滿單元,研究采用考慮雨強(qiáng)的SCS流速計(jì)算公式[11],針對(duì)處于匯流路徑的未蓄滿柵格單元,提出考慮降水強(qiáng)度和土壤含水量的時(shí)變流速計(jì)算公式：

(7)

式中：Vj,αt為當(dāng)流域蓄滿比例為αt時(shí)j柵格單元的匯流速度；k為流速系數(shù),依據(jù)下墊面植被、土壤類型進(jìn)行確定；βj為j柵格坡度；Ic為流域參考雨強(qiáng)；Wj,t為j柵格t時(shí)段蓄水量。

由式(7)可知,當(dāng)某一柵格屬于河道單元時(shí),流速固定為2 m/s[6]；當(dāng)該柵格處于未蓄滿區(qū)域且未處于蓄滿柵格匯流路徑上時(shí),柵格流速為0,不參與匯流；當(dāng)該柵格處于蓄滿區(qū)域時(shí),流速跟時(shí)變雨強(qiáng)有關(guān)；當(dāng)該柵格處于未蓄滿區(qū)域,但上游存在蓄滿區(qū)域時(shí),需要同時(shí)考慮雨強(qiáng)和土壤含水量對(duì)柵格流速的影響。

其中,未蓄滿柵格的土壤含水量表達(dá)式如下：

(8)

式中：Wj,t-1為j柵格t-1時(shí)段蓄水量；Qj,t為j柵格t時(shí)段的上游入流貢獻(xiàn)量;Pt為t時(shí)段降水量。

(9)

(10)

式中：Lj為柵格邊長(zhǎng)。

(11)

在此基礎(chǔ)上,將蓄滿柵格面積轉(zhuǎn)換為流域出口流量過(guò)程,得到任意時(shí)段的分布式單位線：

(12)

式中：Qm為第m個(gè)時(shí)段的單位線流量,m3/s；H為單位凈雨,mm。

1.4 參數(shù)率定與目標(biāo)函數(shù)

研究全面考慮洪水過(guò)程低流量誤差、高流量誤差、洪量誤差等多目標(biāo)因素,其中,高流量誤差強(qiáng)調(diào)洪峰預(yù)報(bào)誤差[18]。選擇Kling-Gupta系數(shù)(EKG)[19]、納什效率系數(shù)(ENS)[20]、納什效率系數(shù)對(duì)數(shù)形式及均方根誤差與標(biāo)準(zhǔn)差的比值(RSR)4個(gè)目標(biāo)參與優(yōu)化率定,優(yōu)化算法采用SCE-UA算法[21]。其中,EKG和ENS包含的平方誤差項(xiàng)強(qiáng)調(diào)了高流量的預(yù)報(bào)精度；RSR系數(shù)通過(guò)均方根誤差與標(biāo)準(zhǔn)差的比值量化了洪量誤差；ENS通過(guò)取對(duì)數(shù)壓縮變量尺度,強(qiáng)調(diào)了低流量預(yù)報(bào)精度。為全面考慮4個(gè)目標(biāo),建立了涵蓋上述指標(biāo)的綜合目標(biāo)函數(shù),依據(jù)Brunner等[22]研究成果,采用多目標(biāo)加權(quán)的方式對(duì)水文模型進(jìn)行優(yōu)化率定,表達(dá)式為

O=w1(1-ENS)+w2(1-EKG)+w3lg(1-ENS)+w4RSR

(13)

式中：w1、w2、w3、w4分別為4個(gè)目標(biāo)權(quán)重。本文重點(diǎn)關(guān)注匯流精度,洪峰準(zhǔn)確度較為重要,因此,ENS和EKG的比重相對(duì)較高,RSR和log(1-ENS)的比重相對(duì)較小。

2 例證研究

2.1 研究區(qū)域及數(shù)據(jù)

龍虎圩水和東石河是中國(guó)南方山區(qū)小河流,位于廣東省韓江流域,分別匯入程江和柚樹河,地處山區(qū)丘陵地帶,地形復(fù)雜,屬亞熱帶氣候,氣候高溫濕熱、暴雨頻繁,但時(shí)空分布不均,具有河床坡降陡、天然落差大、洪水易漲易退等特點(diǎn),降水徑流存在顯著非線性關(guān)系,容易引起山洪災(zāi)害。為此,研究選擇以上2個(gè)小流域?yàn)檠芯繉?duì)象,其中,龍虎圩水河長(zhǎng)20.8 km,流域面積130.8 km2,河流平均比降2.92‰；東石河全長(zhǎng)23.6 km,流域面積152.4 km2,河流平均比降3.56‰。流域DEMs分辨率為30 m×30 m,在此基礎(chǔ)上提取流域柵格坡度、流向等信息;流域植被覆蓋數(shù)據(jù)來(lái)自清華大學(xué)植被數(shù)據(jù)庫(kù)(data.ess.tsinghua.edu.cn)；用于參數(shù)率定和檢驗(yàn)的28場(chǎng)降水徑流數(shù)據(jù)由廣東省水文局梅州水文分局提供,其中,20場(chǎng)用于參數(shù)率定,8場(chǎng)用于結(jié)果檢驗(yàn)。

2.2 參數(shù)率定結(jié)果

選取2個(gè)流域近年各10場(chǎng)洪水進(jìn)行參數(shù)率定。模型率定分為產(chǎn)流參數(shù)率定、單位線參數(shù)率定、匯流參數(shù)率定3個(gè)階段。需要率定的參數(shù)有初損系數(shù)(η)、降雨前期流域特征參數(shù)(CN)、Ic、表征蓄水容量空間分布均勻性的系數(shù)(n)、WM，min及WM，max,其余參數(shù)(流域坡度、流速系數(shù))可根據(jù)流域下墊面信息確定。

(1) 根據(jù)歷史場(chǎng)次洪水過(guò)程推算流域凈雨量,以實(shí)際凈雨量與SCS模型預(yù)報(bào)凈雨量差值最小為依據(jù)率定SCS模型產(chǎn)流參數(shù)CN和η,實(shí)際凈雨量通過(guò)流域出口斷面流量過(guò)程反算得到。

(2) 采用試算法率定分布式單位線參數(shù)Ic。由于目標(biāo)流域場(chǎng)次洪水的降水強(qiáng)度多集中在5～15 mm/h,根據(jù)水量平衡原理重新分配后,蓄滿區(qū)雨強(qiáng)范圍多為20～60 mm/h,在參數(shù)Ic率定時(shí),先根據(jù)流域歷史平均降水強(qiáng)度假定幾組離散值,本文假定20 mm、40 mm、60 mm,然后在其中優(yōu)選,具體參見文獻(xiàn)[11]；此外,推求分布式單位線還需獲取流速系數(shù)、流域坡度等信息,其中,流速系數(shù)依據(jù)植被覆蓋情況進(jìn)行確定,不同植被類型對(duì)應(yīng)的流速系數(shù)見文獻(xiàn)[6],流域植被類型、地形坡度通過(guò)Arcgis工具提取。

(3) 根據(jù)時(shí)段凈雨及式(12)進(jìn)行匯流計(jì)算,該過(guò)程需要依據(jù)流域蓄滿比例選擇單位線,流域蓄滿比例的計(jì)算關(guān)鍵在于確定流域蓄水容量空間分布,因此,采用1.4節(jié)所提方法率定蓄水容量參數(shù)n、WM,max及WM,min,目標(biāo)函數(shù)權(quán)重確定參考文獻(xiàn)[22],w1—w4分別為0.5、0.25、0.15、0.1。根據(jù)率定后的流域CN值計(jì)算流域潛在蓄水能力(S),以該值作為流域平均蓄水容量,進(jìn)一步試算參數(shù)n和WM,max,研究發(fā)現(xiàn),在流域平均蓄水容量不變的前提下,參數(shù)WM,min不敏感,變化范圍較小,本文取值5 mm,參數(shù)n、WM,max對(duì)流域蓄水容量的空間分布有顯著影響,參數(shù)n取值越大,蓄水容量空間分布越不均勻,WM,max取值也越大。本文參數(shù)率定結(jié)果如表1。

表1 龍虎圩流域和東石流域參數(shù)取值方案Table 1Parameter schemes of the Longhuwei River basin and Dongshi River basin

表2 龍虎圩流域和東石流域場(chǎng)次洪水率定結(jié)果及方案選擇Table 2Calibration accuracy and schemes selection results in the Longhuwei River basin and Dongshi River basin

由表2可知,龍虎圩流域和東石流域率定期洪峰相對(duì)誤差在±20%以內(nèi),峰現(xiàn)時(shí)間誤差均在±6 h之間,2個(gè)流域的ENS均值分別為0.86和0.82,場(chǎng)次洪水合格率超80%。

2.3 動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線推求

根據(jù)2.2節(jié)確定的單位線參數(shù)可以推求考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線和現(xiàn)有時(shí)變分布式單位線,為方便工程實(shí)際應(yīng)用,結(jié)合相關(guān)研究[10-11,23],本文采用忽略上游入流貢獻(xiàn)的式(8)進(jìn)行未蓄滿柵格土壤含水量計(jì)算；同時(shí),為減少單位線條數(shù),將降水強(qiáng)度和流域蓄滿比例離散化,凈雨量和流域蓄滿比例因子離散區(qū)間見表3和表4。

表3 每個(gè)時(shí)段的雨強(qiáng)對(duì)應(yīng)的離散值Table 3Discrete value Is of excess rainfall intensity It/40 corresponding to each time interval

表4 每個(gè)時(shí)段的土壤含水量狀態(tài)對(duì)應(yīng)離散土壤含水狀態(tài)Table 4Discrete value αs of soil moisture content αt corresponding to each time interval

以龍虎圩流域?yàn)槔?獲得蓄水容量空間分布后,根據(jù)其蓄滿比例獲得具體蓄滿區(qū)域,不同蓄滿比例下(以αt取0.25、0.50、0.75為例),流域蓄滿區(qū)域分布如圖1所示,結(jié)合表3、表4和流域蓄滿面積可以推求考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線。以推求蓄滿比例為0.25時(shí)的單位線為例,即圖1(a)中綠色區(qū)域所產(chǎn)生的單位凈雨在流域出口形成的流量過(guò)程線(圖2中綠色的單位線),相應(yīng)地,可以推求不同蓄滿比例下的匯流單位線,現(xiàn)有時(shí)變分布式單位線結(jié)果如圖3所示。

圖1 龍虎圩流域不同蓄滿比例對(duì)應(yīng)的產(chǎn)流面積空間分布Fig.1 Runoff areas corresponding to different full storage ratio in the Longhuwei River basin

由圖2可知,在同一雨強(qiáng)下,考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線在流域蓄滿比例較小時(shí)單位線峰值也較小,單位線上漲歷時(shí)略有提前趨勢(shì),但不顯著；在同一蓄滿比例下,隨著降水強(qiáng)度的增加,流域匯流速度增加,匯流時(shí)間減少,單位線峰現(xiàn)時(shí)間提前,總歷時(shí)減小。本文所提考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線的物理意義為單位時(shí)段、單位凈雨均勻降落在蓄滿區(qū)域,在流域出口形成的流量過(guò)程線；圖3為現(xiàn)有時(shí)變分布式單位線,其物理意義為單位時(shí)段、單位凈雨均勻降落在全流域,在流域出口形成的流量過(guò)程線。兩者在降水及下墊面時(shí)空分布、匯流速度及匯流路徑等方面均存在顯著差異。

圖2 龍虎圩流域考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線Fig.2 Time-varying distributed unit hydrographs considering dynamic flow routing paths for the Longhuwei River basin

圖3 龍虎圩流域現(xiàn)有時(shí)變分布式單位線Fig.3 Current time-varying distributed unit hydrographs for the Longhuwei River basin

2.4 與現(xiàn)有分布式單位線對(duì)比分析

根據(jù)SCS模型產(chǎn)流計(jì)算結(jié)果,分別采用本文所提單位線和現(xiàn)有分布式單位線進(jìn)行模型檢驗(yàn),選取龍虎圩流域和東石流域各4場(chǎng)洪水進(jìn)行預(yù)報(bào),當(dāng)流域前期較為干燥時(shí)選用參數(shù)方案1,否則采用參數(shù)方案2,預(yù)報(bào)結(jié)果如表5。

表5 本文方法和現(xiàn)有方法預(yù)報(bào)結(jié)果對(duì)比Table 5Comparisons of the results predicted by the proposed and current methods

由表5可知,采用考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線進(jìn)行匯流計(jì)算后,龍虎圩流域和東石流域場(chǎng)次洪水平均納什效率系數(shù)在0.85以上,洪峰相對(duì)誤差均值約為8%,平均峰現(xiàn)時(shí)間為1.65 h,預(yù)報(bào)精度達(dá)到工程應(yīng)用乙級(jí)及以上精度,且優(yōu)于現(xiàn)有時(shí)變分布式單位線預(yù)報(bào)結(jié)果,其中龍虎圩流域預(yù)報(bào)效果優(yōu)于東石流域。以龍虎圩流域20030517號(hào)、20120527號(hào)場(chǎng)次洪水為例展開分析,采用本文所提動(dòng)態(tài)匯流單位線進(jìn)行徑流預(yù)報(bào)時(shí),這2場(chǎng)洪水預(yù)報(bào)精度有較大提升,2場(chǎng)洪水預(yù)報(bào)徑流過(guò)程如圖4所示；龍虎圩流域20060601號(hào)、東石流域20161021號(hào)、20190609號(hào)場(chǎng)次洪水存在精度提升不明顯或各指標(biāo)沒(méi)有同步提升現(xiàn)象,原因如下：龍虎圩流域20060601號(hào)洪水前期土壤含水量達(dá)106 mm,流域接近蓄滿狀態(tài),動(dòng)態(tài)流徑的時(shí)變分布式單位線和現(xiàn)有方法推求的單位線差距較??；東石流域20161021號(hào)、20190609號(hào)洪水產(chǎn)流參數(shù)不具有普適性。

圖4 場(chǎng)次洪水預(yù)報(bào)流量與實(shí)測(cè)流量對(duì)比Fig.4 Comparisons of the predicted and observed flow for flood events 20030517 and 20120528

由于采用相同的產(chǎn)流參數(shù),采用2種單位線進(jìn)行預(yù)報(bào)的產(chǎn)流量是一樣的,因此,預(yù)報(bào)洪水過(guò)程誤差主要由單位線的差異造成。由圖4可知,采用現(xiàn)有時(shí)變分布式單位線進(jìn)行匯流計(jì)算時(shí),20030517號(hào)洪水和20120527號(hào)洪水預(yù)報(bào)誤差均高于所提方法,其中,20030517號(hào)洪水洪峰相對(duì)誤差超過(guò)14%,ENS為0.88；20120527號(hào)洪水峰現(xiàn)時(shí)間誤差為3 h,ENS為0.87。采用所提單位線進(jìn)行預(yù)報(bào)后,2場(chǎng)洪水過(guò)程的ENS分別提高至0.94和0.93。結(jié)果表明,考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線可以更好地表征流域的匯流過(guò)程。

2.5 流域匯流時(shí)間規(guī)律解析

分布式單位線的核心為柵格匯流時(shí)間的推求,為解析動(dòng)態(tài)匯流路徑對(duì)流域匯流時(shí)間的影響機(jī)制,下面深入分析流域匯流時(shí)間的動(dòng)態(tài)變化規(guī)律。采用動(dòng)態(tài)匯流路徑后,有一部分蓄滿柵格由于流經(jīng)未蓄滿區(qū)域,導(dǎo)致匯流時(shí)間發(fā)生變化,為比較流域不同蓄滿比例下土壤含水量對(duì)柵格匯流速度的影響,離散凈雨強(qiáng)度設(shè)定為最大值Is=2.0,通過(guò)提取流域最先達(dá)到蓄滿狀態(tài)的25%流域面積(F0%—F25%)的柵格匯流時(shí)間,分析該部分面積上柵格在不同流域蓄滿比例下的匯流時(shí)間變化,并與采用現(xiàn)有時(shí)變流速計(jì)算公式的匯流時(shí)間進(jìn)行比較,以龍虎圩流域?yàn)槔?流域最先蓄滿的25%面積上的柵格匯流時(shí)間分布如圖5。

圖5 龍虎圩流域25%蓄滿柵格匯流時(shí)間分布Fig.5 Flow routing time distribution of the 25% full storage grid in the Longhuwei River basin

圖5展示了圖1(a)中綠色區(qū)域的柵格匯流時(shí)間,αt=1.00的橙色區(qū)域代表現(xiàn)有時(shí)變流速計(jì)算公式的結(jié)果,橫坐標(biāo)代表柵格,按照蓄水容量從小到大排序,以圖5F18%處對(duì)應(yīng)的柵格為例,其在蓄滿狀態(tài)為0.25、0.50、0.75、1.00時(shí)對(duì)應(yīng)的匯流時(shí)間分別為9 h、7 h、6 h、5 h,可以看出隨著流域蓄滿比例的增加,蓄滿柵格的匯流時(shí)間逐漸減少,這是因?yàn)殡S著降水的持續(xù),處于蓄滿路徑的未蓄滿柵格單元也逐漸趨于飽和,未蓄滿柵格單元逐漸減少,由式(7)可知,柵格蓄滿后流速達(dá)到最大值,因此,隨著流域蓄滿比例的增加,柵格匯流時(shí)間逐漸縮短,直至流域達(dá)到全蓄滿狀態(tài),此時(shí)動(dòng)態(tài)匯流路徑轉(zhuǎn)化為靜態(tài)匯流路徑,流域匯流時(shí)間趨于穩(wěn)定。

圖6為龍虎圩流域不同蓄滿比例下考慮動(dòng)態(tài)匯流路徑時(shí)計(jì)算的匯流時(shí)間和靜態(tài)匯流路徑計(jì)算的匯流時(shí)間箱線圖,結(jié)果表明,現(xiàn)有時(shí)變流速公式會(huì)高估坡面匯流速度,匯流時(shí)間不隨蓄滿比例而改變,考慮土壤含水量對(duì)流域匯流速度的影響后,隨著流域逐漸達(dá)到全部蓄滿狀態(tài),即αt趨于1.00時(shí),2種方法匯流時(shí)間趨于一致,匯流時(shí)間分布更加集中；此外,隨著流域蓄滿面積的增加,土壤含水量對(duì)匯流時(shí)間的影響有減小趨勢(shì),所提匯流速度計(jì)算方程對(duì)前期較為干旱的情況影響較大。

圖6 本文所提方法和現(xiàn)有方法不同蓄滿比例下流域匯流時(shí)間箱線圖Fig.6 Box plots of flow routing time under different storage ratio of the proposed and current methods

3 結(jié) 論

(1) 針對(duì)現(xiàn)有模型中的蓄水容量曲線只能定性描述流域蓄水容量空間不均勻分布的缺陷,引入地形指數(shù)刻畫流域柵格蓄水容量的空間不均勻分布,定量表征每一個(gè)柵格蓄水容量,結(jié)合產(chǎn)流計(jì)算結(jié)果,進(jìn)而從柵格尺度劃分流域蓄滿柵格和未蓄滿柵格,可為柵格流速場(chǎng)計(jì)算提供更高精度輸入場(chǎng)。

(2) 提出了考慮時(shí)變降水強(qiáng)度和土壤含水量的匯流速度計(jì)算公式,通過(guò)劃分4種柵格狀態(tài)(蓄滿柵格、未蓄滿但處于匯流路徑柵格、未蓄滿且不處于匯流路徑柵格、河道柵格),更加精確地刻畫了流速場(chǎng)的空間分布規(guī)律。

(3) 建立了考慮動(dòng)態(tài)匯流路徑的時(shí)變分布式單位線,針對(duì)傳統(tǒng)單位線假定全流域均勻產(chǎn)流的不合理假定,本文以南方山區(qū)龍虎圩流域和東石流域?yàn)檠芯繉?duì)象,基于蓄滿區(qū)域更易產(chǎn)流的理論基礎(chǔ),推求了流域不同蓄滿比例下的動(dòng)態(tài)匯流路徑,完善了單位線推求的理論基礎(chǔ)。

(4) 結(jié)合本文所提改進(jìn)時(shí)變流速計(jì)算公式和動(dòng)態(tài)匯流路徑理論,計(jì)算了不同蓄滿比例下柵格匯流時(shí)間分布場(chǎng),具體而言,通過(guò)改進(jìn)流速公式提高了匯流時(shí)間計(jì)算的精度,通過(guò)劃分蓄滿區(qū)域和未蓄滿區(qū)域,將匯流方法的匯流路徑從靜態(tài)匯流路徑發(fā)展為動(dòng)態(tài)匯流路徑,應(yīng)用結(jié)果表明,所提方法應(yīng)用精度較高,提高了分布式單位線匯流時(shí)間計(jì)算的準(zhǔn)確性。