基于地貌單位線的匯流模型在陳河流域的構(gòu)建與應(yīng)用

童冰星，李致家，溫婭惠，孫明坤

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

0 引 言

陳河流域位于北溫帶，屬大陸性季風氣候區(qū)，多年平均降水量700～900 mm；河流水量主要系雨水補給，局部暴雨是發(fā)生洪水的主要原因。流域地處秦嶺北麓，為秦嶺褶皺帶，上游大部分為高中山區(qū)，海拔高程700～3 500 m，山勢雄偉，地形起伏大，大部分區(qū)域的坡度在25°以上。這些因素導致陳河流域洪水陡漲陡落[1-5]，同時給流域的防洪工作帶來極大的挑戰(zhàn)；而經(jīng)常使用的線性水庫、滯后演算、馬斯京根法逐級匯流演算的方法難以對這種快速漲落的洪峰進行有效準確的模擬[6- 8]。

基于地貌擴散與水動力擴散的流域瞬時地貌單位線法是一種有一定物理基礎(chǔ)的流域匯流隨機模型[9-12]。它主要是在數(shù)字高程模型(DEM)的基礎(chǔ)上，根據(jù)流域匯流過程與地貌因子間的關(guān)系描述流域匯流的過程[13]。它所需要的參數(shù)少且對輸入條件響應(yīng)快，因而本文在陳河流域上以該單位線為基礎(chǔ)構(gòu)建匯流模型，并著重對洪水的洪峰部分進行模擬演算。

1 基于地貌擴散與水動力擴散的流域瞬時地貌單位線的原理

2002年石朋等提出的基于地貌擴散與水動力擴散的流域瞬時地貌單位線[14]中指出：地貌擴散作用一般取決于流域的大小、形狀等因素，主要反映流域中水系的疏密以及分布的形狀等因素對于流域匯流的影響。它可以由Surka于1968年提出的寬度函數(shù)[15-16]來描述。寬度函數(shù)所描述的是距流域出口斷面一定距離處的鏈數(shù)與距離之間的函數(shù)關(guān)系。其中的距離可以用拓撲長度來描述，也可以用幾何長度來表示。水動力擴散作用一般與流域上的流速分布有關(guān)系，產(chǎn)生的根源在于流域上的匯流速度分布不均[17-18]，水動力擴散對匯流的影響可以用河道響應(yīng)函數(shù)[19-21]來描述。基于寬度函數(shù)和河道響應(yīng)函數(shù)可以得到流域瞬時地貌單位線。由于本文中的陳河流域處于高中山區(qū)，坡度陡峻(見圖1)，河道中的洪水波一般屬于運動波；因此，本文采用結(jié)合運動波理論的流域瞬時地貌單位線

(1)

式中，n為水系河源數(shù)；μi為水系中內(nèi)鏈的平均長度；V為河道洪水波平均波速；t為時刻。

為求出1 h流域地貌單位線，應(yīng)將式(1)積分，得到流域的S曲線，然后用S曲線計算得到流域1 h單位線

(2)

u(Δt)=S(t)-S(t-Δt)

(3)

式中，Δt為單位線時段長，此處取1 h。

圖1 陳河流域坡度概率分布

2 陳河流域瞬時地貌單位線的提取

由式(2)可知，基于地貌擴散與水動力擴散的流域瞬時地貌單位線主要是由水系河源數(shù)、水系中內(nèi)鏈的平均長度以及河道洪水波的平均波速這三個參數(shù)確定的。

2.1 水系河源數(shù)和水系中內(nèi)鏈的平均長度的提取

水系河源數(shù)和水系中內(nèi)鏈的平均長度均為流域水系的特征值，因此可以從流域DEM數(shù)據(jù)中提取。本文利用地理信息系統(tǒng)軟件以90 m分辨率的DEM數(shù)據(jù)為基礎(chǔ)提取流域的數(shù)字水系，并將提取出的水系與陳河流域1∶55萬的Google地圖影像進行對比(見圖2)，確認兩者相吻合后從數(shù)字水系中提取出陳河流域的河源數(shù)和水系中內(nèi)鏈的平均長度(見表1)。

圖2 提取的水系與Google地圖對比

河源數(shù)目平均長度/m外鏈內(nèi)鏈鏈水系總長度/m河道洪水波平均波速/m·s2133664311382725693127

2.2 流域河道洪水波平均波速的計算

對于流域河道洪水波平均波速v，本文對其進行率定。率定時采用陳河流域2003年到2012年之間的13場洪水資料，對每一場洪水均以0.1m/s為步長，將v從1 m/s到10 m/s之間連續(xù)取90個值逐次計算，并對計算得到的洪峰相對誤差進行統(tǒng)計。

對計算結(jié)果進行統(tǒng)計分析，計算每取一個v值時洪峰相對誤差均值Δi，以及取每一個v值時的洪峰相對誤差平方均值Δj方。即

(4)

(5)

式中，j為v的值，以0.1為步長從1變化到10；i為每一場洪水的編號；n為洪水總的場數(shù)，取13。

由式(4)和式(5)計算得到洪峰相對誤差均值隨v變化曲線(圖3)和洪峰相對誤差平方均值隨v值變化曲線(圖4)。

圖3 洪峰相對誤差均值隨v變化

圖4 洪峰相對誤差平方均值隨v變化

由圖3可知：v從1 m/s變化到10 m/s時，計算的13場洪水的洪峰最初較低，相對于實測峰值整體偏小。隨著v的增大，計算出的洪峰逐漸升高，并有部分洪水計算出的洪峰高于實測。且在v=2.7 m/s時計算得到的各場洪水的洪峰相對誤差正負分布較為均勻，洪峰相對誤差均值較小。隨著v的進一步增大，計算出的洪峰普遍高于實測。結(jié)合圖4可知：v從1 m/s變化到10 m/s時，洪峰相對誤差平方均值首先隨著v的增大而減小，在v=2.7 m/s左右時相對誤差的平方均值最小，約為0.024。即各場洪水的洪峰誤差分布最為集中，變幅最小。隨著v的進一步增大洪峰相對誤差平方均值迅速增大。

表3 模擬演算結(jié)果

綜合上述分析，本文采用2.7 m/s作為陳河流域河道洪水波平均波速進行匯流計算較為合適。此時陳河流域瞬時地貌單位線的參數(shù)見表1。

采用表1中的參數(shù)值，結(jié)合式(2)及式(3)得到陳河流域瞬時地貌單位線u(1，t)(見圖5)。

圖5 陳河流域地貌單位線

3 基于地貌單位線匯流模型在陳河流域的應(yīng)用

本文采用三層蒸散發(fā)模型(主要參數(shù)見表2)計

算流域產(chǎn)生的徑流深，將其作為基于地貌單位線匯流模型(下簡稱“地貌單位線模型”)的輸入，對陳河流域的13場洪水的洪峰部分進行模擬演算；并使用同樣采用三層蒸散發(fā)模型進行產(chǎn)流計算的新安江模型，對陳河流域的13場洪水進行模擬演算，將兩個模型計算得到的結(jié)果進行對比分析(見表3)。

表2 三層蒸散發(fā)模型主要參數(shù)值

由表3可知：新安江模型的平均洪峰相對誤差為14.8%，地貌單位線模型的平均洪峰相對誤差為7.6%；平均洪峰現(xiàn)時間誤差為2.4 h，地貌單位線模型的平均峰現(xiàn)時間的誤差為1.7 h。由此可看出，地貌單位線模型的平均洪峰相對誤差和平均峰現(xiàn)時間的誤差均優(yōu)于新安江模型的計算結(jié)果。

基于表3分別對兩個模型的洪峰相對誤差的統(tǒng)計分析也表明，地貌單位線模型優(yōu)于新安江模型計算結(jié)果。

綜合上述分析認為，在陳河流域2003年到2012年之間的13場洪水中，以地貌單位線為基礎(chǔ)的流域匯流模型取得了理想的計算效果。

4 結(jié) 論

本文采用90 m分辨率的DEM數(shù)據(jù)提取了陳河流域水系，同時以0.1 m/s為步長率定分析得到適合于該流域的河道平均匯流速度為2.7 m/s。進而以此為基礎(chǔ)提取出陳河流域基于地貌擴散與水動力擴散的流域瞬時地貌單位線，構(gòu)建了以地貌單位線為基礎(chǔ)的流域匯流模型，并采用該模型與新安江模型對陳河流域上2003年到2012年間共13場洪水進行了流域匯流演算。計算結(jié)果的對比分析表明，在陳河流域上以地貌單位線為基礎(chǔ)的流域匯流模型的取得了理想的計算效果。

然而，在將地貌單位線模型計算得到的洪水過程線與實測的流域出流過程線相比較時也發(fā)現(xiàn)：洪峰部分模擬的效果較好，但是漲洪與落洪部分常常過于快速，模擬效果并不是十分理想。這可能是由于將三層蒸散發(fā)計算出來的產(chǎn)流量全部作為單位線匯流演算的輸入，沒有針對不同的徑流進行徑流分割來分別作匯流演算，造成漲落的過程過快。因此，如何進行徑流劃分并分別對不同徑流采用相適合的方法進行匯流演算，是需要進一步研究的問題。

[1] 芮孝芳. 水文學研究進展[M]. 南京： 河海大學出版社, 2007．

[2] 張恭肅, 黃守信, 賀偉程. 小流域單位線的非線性分析[J]. 水利學報, 1981(3)： 1- 9.

[3] 李致家, 姜婷婷, 黃鵬年, 等. 降雨和地形地貌對水文模型模擬結(jié)果的影響分析[J]. 水科學進展, 2015, 26(4)： 473- 480.

[4] 李巧玲, 王榮克, 董小濤, 等. 涇河上游植被覆蓋動態(tài)及其與降雨徑流的關(guān)系[J]. 水力發(fā)電, 2015, 41(11)： 21- 33．

[5] 姚成, 孫如飛, 李致家, 等. 下墊面變化條件下合河流域設(shè)計洪水修訂[J]. 水力發(fā)電, 2015, 41(11)： 9- 13．

[6] 安東, 李致家, 闞光遠, 等. 數(shù)據(jù)驅(qū)動模型與概念性模型的應(yīng)用對比[J]. 水力發(fā)電, 2013, 39(12): 9- 12．

[7] 黃小祥, 姚成, 李致家, 等. 柵格新安江模型在天津于橋水庫流域上游的應(yīng)用[J]. 湖泊科學, 2016, 28(5)： 1134- 1140.

[8] 李致家, 包紅軍, 孔祥光, 等. 水文學與水力學相結(jié)合的南四湖洪水預(yù)報模型[J]. 湖泊科學, 2005, 17(4)： 299- 304.

[9] 陸桂華. 確定性方法推求地貌單位線[J]. 河海大學學報: 自然科學版, 1990, 18(6)： 80- 83.

[10] 芮孝芳. 由流路長度分布律和坡度分布律確定地貌單位線[J]. 水科學進展, 2003, 14(5)： 602- 606．

[11] 石朋, 芮孝芳. 由SSNs構(gòu)建流域地貌瞬時單位線的研究[J]. 水科學進展, 2005, 16(6)： 799- 803.

[12] 芮孝芳. 地貌瞬時單位線理論的若干評論[J]. 水科學進展, 1991, 2(2)： 195- 200．

[13] 胡健偉, 陸桂華, 吳志勇. 基于地理信息系統(tǒng)技術(shù)的GIUH通用公式的應(yīng)用[J]. 河海大學學報： 自然科學版, 2005, 33(3)： 269- 272.

[14] 芮孝芳, 石朋. 基于地貌擴散與水動力擴散的流域瞬時地貌單位線[J]. 水科學進展, 2002, 13(4)： 440- 444．

[15] SURKAN A J. Synthetic hydrographs: effects of network geometry[J]. Water Resources Research, 1968, 5(1): 112- 128．

[16] 包為民. 水文預(yù)報[M]. 北京： 中國水利水電出版社, 2009．

[17] 孔凡哲, 芮孝芳, 李燕. 基于空間分布流速場的單位線推求及應(yīng)用[J]. 河海大學學報：自然科學版, 2006, 34(5)： 485- 488．

[18] 孔凡哲, 芮孝芳. 一種地形指數(shù)計算方法在Topmodel洪水模擬計算中的應(yīng)用[J]. 水文, 2003, 23(3)： 16- 19.

[19] RINALDO A, MARANI A, RIGON R. Geomorphological dispersion[J]. Water Resources Research, 1991, 3(4): 512- 525．

[20] TROUTMAN B M, KARLINGER M R. Unit hydrograph approximations assuming linear flow through topologically random channel networks[J]. Water Resources Research, 1985, 21(5): 743- 754．

[21] RODRIGUEZ-ITURBE I, VALDES J B. The geomorphological structure of hydrologic response[J]. Water Resources Research, 1979, 15(6): 1409- 1420．