肖 楠，葉 磊，吳 劍，姬榮彬，吳晨晨，周惠成

(大連理工大學水利工程學院，遼寧 大連 116024)

0 引 言

山丘區(qū)小流域因其復(fù)雜的地形條件和暴雨多發(fā)的氣候特征而極易造成多種類型的山洪災(zāi)害，已成為我國自然災(zāi)害的主要災(zāi)種之一。山洪具有突發(fā)性強、陡漲陡落、成峰較快的特點，預(yù)警預(yù)報難度非常大。因此，開展山丘區(qū)小流域水文模擬基礎(chǔ)研究，對于提高山洪災(zāi)害預(yù)警預(yù)報水平，降低人民生命財產(chǎn)損失，具有重大的現(xiàn)實意義。

流域水文模擬涉及水文氣象等多種輸入，依賴于對輸入、輸出信息進行解釋的人為判斷，多種因素導(dǎo)致水文模擬不確定性客觀存在，并始終影響防洪決策的正確性[1]。山丘區(qū)暴雨洪水受局部地形影響大、產(chǎn)匯流非線性顯著，更增加了山丘區(qū)小流域水文模擬的復(fù)雜性、不確定性。降雨作為控制流域水量平衡的主要輸入，是決定模型模擬結(jié)果好壞的關(guān)鍵因素[2]。尹雄銳等[3]通過總結(jié)多項研究成果指出:與模型結(jié)構(gòu)不確定性和模型參數(shù)不確定性相比，降雨輸入是最重要的不確定性因子；Michaud等[4]通過中小尺度流域不同類型水文模型的模擬對比表明降雨輸入精度對模擬結(jié)果影響很大；徐靜等[5]基于模糊集理論，應(yīng)用TOPMODEL對資水流域洪水過程進行模擬研究，結(jié)果表明雨量不確定性的傳播對洪水預(yù)報的影響處于主導(dǎo)地位；Chaubey等[6]發(fā)現(xiàn)降雨的空間分布和觀測精度對流量過程線影響很大，而且在模型參數(shù)率定時，降雨的空間變異性會導(dǎo)致參數(shù)的變異性。

目前國內(nèi)外關(guān)于降雨對水文模擬不確定性影響的研究主要集中在大中流域降雨空間分布的不均勻性和降雨強度的時間變異性對洪水過程特征的影響，而針對山丘區(qū)小流域降雨輸入誤差對水文模擬不確定影響方面的研究較少。不同于傳統(tǒng)大尺度流域水文預(yù)報較多關(guān)注洪水過程，山丘區(qū)小流域則更側(cè)重于峰值的預(yù)報，國內(nèi)外山洪災(zāi)害防治也以成災(zāi)流量或水位作為山洪預(yù)警指標?；诖?，本文選擇河南省欒川流域作為研究對象，通過設(shè)置不同的降雨輸入誤差方案，采用新安江水文模型研究不同降雨輸入誤差對洪峰模擬的影響，以期建立降雨輸入誤差與洪峰流量模擬結(jié)果之間的定量關(guān)系，為提高山洪預(yù)警預(yù)報水平提供實際支撐。

1 研究區(qū)概況

欒川縣位于河南省西南部，地勢起伏跌宕，地處東經(jīng)111°12′～111°02′，北緯33°39′～34°11′，流域面積343 km2。流域?qū)倥瘻貛Т箨懶约撅L氣候，氣候溫涼，雨量較多，冬長夏短。年平均氣溫12.1 ℃，年日照時數(shù)2 103 h，無霜期198 d，年平均降雨量784.7 mm，降雨年內(nèi)分配不均勻，主要集中在7-9月份。欒川流域如圖1所示。

圖1 欒川流域圖Fig.1 Location of Luanchuan basin

2 研究方法

2.1 模型簡介

本文采用新安江水文模型進行模擬計算。新安江水文模型是一個分散性的概念模型，已被廣泛應(yīng)用于我國半干旱半濕潤地區(qū)[7]。模型結(jié)構(gòu)具有三分特點，即分單元、分水源、分階段[8]。分單元指為了考慮降雨空間分布不均以及流域異質(zhì)性將計算流域分為多個單元流域進行產(chǎn)流計算；分水源是指將徑流分成地表、壤中、地下三種水源，它們的匯流速度不同，地表最快、地下最慢；分階段是指將匯流過程分為坡面匯流階段和河網(wǎng)匯流階段，在坡地，各種水源匯流速度不同，而在河網(wǎng)則無此差別。

2.2 數(shù)據(jù)資料可靠性分析

水文序列原始資料的可靠性直接影響水文模擬計算結(jié)果的精度。陳玲玲[9]等曾通過對降雨徑流系列建立相關(guān)關(guān)系來檢驗水文資料的可靠性，所以為對研究應(yīng)用的歷史雨洪資料進行可靠性分析，本文根據(jù)欒川流域多年降雨徑流系列，采用線性、二次多項式、對數(shù)、指數(shù)、乘冪五種不同曲線類型進行擬合，建立其相關(guān)關(guān)系式并計算相關(guān)系數(shù)，擬合結(jié)果見表1，其中x為年降雨量，y為年徑流量。

表1 欒川流域年降雨徑流擬合結(jié)果Tab.1 Fitting results between yearly rainfall andrunoff in Luanchuan Basin

由表1擬合結(jié)果可知，欒川流域年降雨量與徑流量相關(guān)系數(shù)很高，具有很強的相關(guān)性，表明水文資料系列的可靠性較高，可為后續(xù)開展降雨輸入對洪峰模擬不確定性影響研究提供較好的數(shù)據(jù)支持。

2.3 降雨輸入誤差方案設(shè)計

為研究降雨輸入對山丘區(qū)小流域洪峰模擬的不確定性，本文通過對降雨量的放大縮小來設(shè)計不同降雨輸入誤差方案。欒川流域共7個雨量站，經(jīng)抽站法計算得到當4個雨量站觀測數(shù)據(jù)有效時，數(shù)據(jù)允許誤差為20%的保證率為81%，一定程度上可以滿足水文模擬的需求，而在實際情況中3個雨量站同時出現(xiàn)故障的概率很小。因此，在方案設(shè)計時取±20%為降雨輸入誤差的上限，加設(shè)誤差為±5%、±10%四種方案。方案設(shè)置具體步驟如下：①采用泰森多邊形法計算流域面雨量，作為降雨“真值”；②對降雨“真值”序列按±5%、±10%、±20%的誤差進行放大或縮小，生成6種不同誤差條件下的面雨量序列；③將各種方案輸入模型進行模擬，為消除模型參數(shù)的影響，每種方案均采用由降雨“真值”率定所得參數(shù)模擬。

2.4 目標函數(shù)選取與參數(shù)率定

在大流域的水文模擬研究中，為了能夠全面反映流域水文過程中不同的水文變化特性，一般要設(shè)置多個目標函數(shù)對水文模型參數(shù)進行率定[10]。由于山丘區(qū)河道坡降較大，河床狹窄，行洪區(qū)小，調(diào)蓄能力差，導(dǎo)致其洪水具有陡漲陡落、歷時較短、峰型尖瘦的特點[11]，洪水預(yù)報過程中在洪峰處易發(fā)生顯著突變，洪峰預(yù)報誤差較大。山區(qū)小流域災(zāi)害的發(fā)生取決于洪峰流量是否超過成災(zāi)流量，洪峰流量大小也直接影響山洪災(zāi)害破壞性的大小，準確預(yù)報洪峰流量對于提高山洪預(yù)警預(yù)報水平至關(guān)重要。因此，為提高洪峰預(yù)報精度，參考洪水預(yù)報精度評定標準，將洪峰流量相對誤差小于20%的場次視為合格，以洪峰流量合格場次越多越好作為目標函數(shù)進行參數(shù)率定，峰值相對誤差作為評定指標，計算如式(1)所示。與多目標參數(shù)率定相比，本文參數(shù)率定部分更加突出了山洪過程的洪峰特征，根據(jù)預(yù)報結(jié)果能夠更加準確的做出災(zāi)害防治決策。

(1)

式中：Qsim、Qobs分別為場次洪峰流量的模擬值與實測值, m3/s；Qrel為場次洪峰流量相對誤差,%。

根據(jù)歷史雨洪資料，在欒川流域共選取了13場典型洪水。經(jīng)率定驗證，洪峰流量合格率可達93.33%，新安江模型參數(shù)如表2所示。由表2可知，參數(shù)率定值均在經(jīng)驗取值范圍內(nèi)，滿足物理意義，可適用于該流域水文模擬。

表2 欒川流域新安江模型參數(shù)Tab.2 The parameters of the Xinanjiang modelin the Luanchuan Basin

3 結(jié)果分析

3.1 降雨輸入對洪峰模擬結(jié)果的影響分析

采用率定后的新安江水文模型參數(shù)，按照不同的降雨誤差輸入方案分別進行模擬計算。本文以實測的降雨數(shù)值作為“真值”，認為其降雨誤差為“0”，以“真值”輸入得出的洪峰流量作為“標準值”，將各方案下的洪峰模擬結(jié)果與“標準值”相比較，得到洪峰模擬偏差。表3給出了欒川流域不同降雨輸入誤差的洪峰模擬偏差結(jié)果。

由表3可知，降雨輸入對山丘區(qū)小流域水文模擬影響較大，以“20030828”號洪水為例，當降雨輸入誤差為5%時，洪峰模擬偏差就已經(jīng)接近11%；當降雨輸入誤差增大到10%時，洪峰模擬偏差增大到22%左右；降雨輸入誤差繼續(xù)增大到20%，相應(yīng)的洪峰模擬偏差達到41%左右。

6種降雨輸入誤差方案下所有場次的洪峰模擬平均偏差統(tǒng)計結(jié)果見表4和圖2。其中“方案0”為利用“真值”的模擬結(jié)果，認為其降雨輸入誤差為0，洪峰平均偏差為0。

從模擬統(tǒng)計結(jié)果中可以看出隨著降雨輸入誤差的增加，洪峰模擬偏差隨之增加，降雨輸入誤差與洪峰模擬偏差之間呈正相關(guān)，表明洪峰流量與降雨密切相關(guān)，降雨輸入誤差對洪峰模擬偏差有顯著影響。

表3 欒川流域不同降雨輸入誤差的洪峰模擬偏差結(jié)果統(tǒng)計表Tab.3 Flood peak simulation deviation results versusdifferent rainfall errors in Luanchuan basin

表4 欒川流域不同降雨輸入誤差的洪峰模擬平均偏差統(tǒng)計結(jié)果Tab.4 Statistical results of simulated average deviation of floodpeak versus different rainfall errors in Luanchuan basin

圖2 降雨輸入誤差與洪峰模擬偏差關(guān)系圖Fig.2 Relation between rainfall input errors and simulated flood peak deviation

3.2 考慮前期土壤干濕程度的模擬結(jié)果分析

山洪洪峰流量與土壤干濕程度密切相關(guān)，相同降雨輸入誤差對于前期土壤干濕程度不同的場次洪水，影響程度不同。為進一步探究不同前期土壤干濕條件下降雨輸入誤差與洪峰模擬偏差之間的規(guī)律，統(tǒng)計所有場次洪水的前期土壤干濕程度(表5)，在此基礎(chǔ)上給出欒川流域不同降雨輸入誤差隨著前期土壤干濕程度變化的洪峰模擬偏差結(jié)果(圖3)，圖中k為擬合直線的斜率。

表5 各場次洪水前期土壤干濕程度Tab.5 Soil water contents at the early stage of the flood

注：經(jīng)參數(shù)率定得本流域WM=169 mm。

注：橫座標表示前期土壤由干到濕。圖3 考慮前期土壤干濕程度的模擬結(jié)果對比圖Fig.3 Comparison of simulation results for soil water contents at the early stage of the flood

由圖3可以看出，隨著土壤由干到濕的變化，洪峰模擬偏差變化較小，表明相比于降雨輸入誤差，前期土壤干濕程度對山丘區(qū)小流域洪峰模擬的影響相對較小。對于相同的降雨輸入誤差，前期土壤濕度越大，洪峰模擬偏差越小。

降雨輸入誤差變化范圍為-20%～+20%，隨著模型降雨輸入的相對增加，趨勢線斜率絕對值也隨之增大，即模型降雨輸入越大，前期土壤干濕程度對小流域洪峰模擬影響越敏感。從-20%～-5%降雨誤差圖可以看出，隨著前期土壤由干到濕變化，洪峰偏差絕對值雖呈逐漸減小趨勢，但減小趨勢不明顯，而當降雨誤差為+5%～+20%時，隨著前期土壤由干到濕變化，洪峰偏差的絕對值明顯變小。總體來看，對于選定的13場洪水，當降雨輸入誤差一定時，次洪前期土壤干濕程度越接近土壤最大蓄水容量WM，則洪峰模擬偏差越小。洪峰模擬偏差與前期土壤干濕程度呈一穩(wěn)定變化趨勢，且該趨勢隨著降雨輸入的增加變得顯著。當前期土壤較濕潤時，若降雨較為充沛，流域土壤相對較容易蓄滿，符合新安江水文模型蓄滿產(chǎn)流規(guī)律。

4 結(jié) 語

降雨作為水文模型輸入的關(guān)鍵數(shù)據(jù)，是水文模擬中重要的不確定性影響因子。本文以河南省欒川流域為研究對象，通過設(shè)置不同的降雨誤差改變新安江水文模型輸入，開展降雨輸入對山丘區(qū)小流域洪峰模擬不確定性影響研究。結(jié)論如下。

(1)洪峰流量與降雨密切相關(guān),降雨輸入誤差對山丘區(qū)小流域洪峰模擬的不確定性影響顯著,降雨輸入誤差的增大會導(dǎo)致洪峰模擬偏差增大，欒川流域洪峰模擬偏差約為降雨輸入誤差的2倍。

(2)對于同一降雨輸入誤差，隨著前期土壤由干到濕變化，洪峰模擬偏差呈減小趨勢，且降雨輸入誤差增大方案(+5%、+10%、+20%)的趨勢相比減小方案(-5%、-10%、-20%)更為顯著。

受山丘區(qū)局部地形、下墊面條件以及實測資料匱乏等因素影響，山丘區(qū)小流域水文模擬難度較大。本文研究成果可為山丘區(qū)小流域水文模擬結(jié)果修正及后期防洪決策提供參考依據(jù)，為提升山洪預(yù)警預(yù)報水平提供基礎(chǔ)支撐。