劉雯雯

(遼寧省沈陽市水文局，遼寧 沈陽 110094)

各大中城市隨著城市化進程的加快以及全球氣候的變化而面臨著峰高量大、短歷時、高強度的暴雨頻發(fā)以及洪水內(nèi)澇等問題的嚴峻調(diào)整。據(jù)不完全統(tǒng)計，我國自2008年以來約60%的城市均發(fā)生了不同程度的洪澇災害，在短期內(nèi)發(fā)生3次以上洪澇災害的城市約140個并以北京“7·21”特大暴雨最為嚴重，不僅嚴重影響著城市市政交通和社會經(jīng)濟的持續(xù)健康發(fā)展，而且對水資源系統(tǒng)平衡構(gòu)成嚴重威脅。海綿城市理論通過一系列的雨洪治理有效措施可顯著減少城市洪澇災害發(fā)生的損失和頻率，并使得城市防洪排澇壓力在很大程度上得到減緩。所以，依據(jù)我國城市基本特征和海綿城市新理念構(gòu)建城市雨洪模型，對于建設海綿城市提高雨洪模擬精度和雨洪管理水平具有十分重要的意義。

水文模型可通過物理方程和數(shù)學公式將較為復雜的水文演化過程進行簡要的描述，從而對較為復雜的水文過程實現(xiàn)抽象化和概念化處理，然而采用該方法經(jīng)常發(fā)生“失真”的情況并以此造成水文模型具有不確定性的特征。城市雨洪模型作為水文模型的重要分支在構(gòu)建過程中往往存在多種不確定性，如對水文水力特征描述失真造成的參數(shù)不確定性和雨洪過程認識不足造成的模型結(jié)構(gòu)不確定性等。模型參數(shù)數(shù)量在水文模型不斷改進、排水過程以及城市降雨徑流描述不斷完善的同時而顯著增大，并且模擬結(jié)果的不確定性隨不同參數(shù)之間復雜的相關性而逐漸增大。對雨洪模型的不確定性保持客觀、合理的認識對于提高雨洪模型精度以及加強模型參數(shù)的認識具有重要意義。

Beven等利用GLUE方法對后驗分布特征以及參數(shù)不確定性開展了深入的研究，并廣泛應用于水文模型中。本文以沈陽市為研究對象通過建立SWMM模型對研究區(qū)域雨洪過程進行了模擬預測，并對模型參數(shù)的不確定性利用GLUE方法進行更深層次的探討分析，以期為沈陽市城市雨洪模型的深入研究和發(fā)展提供一定的理論依據(jù)和技術支持。

1 實例應用與方法

1.1 構(gòu)建模型

沈陽市作為遼寧省的省會城市總面積約1.3萬km2，北依長白山麓、南鄰東遼半島屬于我國東北與渤海地區(qū)的關鍵性結(jié)合區(qū)域，是我國東北地區(qū)經(jīng)濟發(fā)展的核心區(qū)域；沈陽市屬于大陸性季風氣候區(qū)，年降水量約600～800mm，受環(huán)境影響降雨主要集中在夏季，年均氣溫為6.2～9.7℃，冬季寒冷漫長且降雪量較低，夏季短暫且暴雨較為集中；境內(nèi)大小支流約27條主要有繞陽河、遼河、渾河等支流。主城區(qū)商業(yè)較為發(fā)達其不透水面積占比約80%，城市主要以河道和管道作為排水通道。其水文站主要有渾南、沈河以及大東水文站，其水文觀測資料較為完善。

由美國環(huán)保部門首先研發(fā)的主要用于管道和地表徑流計算的SWMM模型，是通過對連續(xù)或單一降雨條件下的水體環(huán)境狀況和水量進行預測分析，對各個河道或管道中水質(zhì)、水質(zhì)狀況、水量進行預測，對各個子流域在任意時刻不同時間步長的水體環(huán)境以及水量變化過程進行時時跟蹤模擬實現(xiàn)城市暴雨洪水的降雨徑流設計與管理的重要手段和方法。對各子流域中透水、有滯和無滯蓄不透水面積3部分所產(chǎn)生的污染負荷、徑流和降雨可利用SWMM模型的徑流模塊部分進行綜合處理；而模型的匯流計算模塊主要是對水泵、蓄水、渠道等水利調(diào)節(jié)措施進行水量傳輸?shù)哪M處理，其中動力波、恒定流以及運動波理論方法是用于計算河道匯流的常用方法。在解決城市暴雨洪水過程模擬、排水系統(tǒng)、雨洪調(diào)蓄處理規(guī)劃設計分析及其他相關領域，SWMM模型已得到了廣泛的推廣和應用。

研究區(qū)域的數(shù)字化下墊面處理是SWMM模型進行雨洪預測分析的重要基礎和前提條件，依據(jù)其基本內(nèi)涵可將下墊面數(shù)字化劃分3大模塊即子流域、排水系統(tǒng)的信息概化與提取、子流域邊界條件的劃分。管網(wǎng)系統(tǒng)通過節(jié)點連接而組成的體系即為城市排水系統(tǒng)，可分為河渠和地下管道2種類型?？紤]到管網(wǎng)資料的可獲取性，研究區(qū)域內(nèi)排水管網(wǎng)主要沿道路布設，遵循就近原則，按照城市排水檢查井節(jié)點對各個匯水子流域的地表進行設置。結(jié)合研究區(qū)域?qū)嶋H狀況對SWMM模型進行概化處理，各子流域共采用13個河道管網(wǎng)進行控制。因此本研究概化條件中河道、道路排水通道以及總排水通道分別為10、10、21條。結(jié)合相關研究資料，對參數(shù)的不確定性開展了詳細的分析探討。

雨洪模擬數(shù)據(jù)采用沈陽市2015—2016年的5場暴雨洪水水文相關資料，各場次暴雨參數(shù)見表1。依據(jù)研究區(qū)域的雨量站提取降雨數(shù)據(jù)資料，而對分洪流量以及降雨徑流資料分別由各區(qū)段內(nèi)的水文站和河道閘室提供。

表1 沈陽市2015—2016年各場次暴雨特征

1.2 GLUE方法

模擬結(jié)果主要與一系列參數(shù)的選取密切相關而與單一參數(shù)的不確定性相關不大，此條件是利用GLUE法的核心理論。其基本理論是對參數(shù)取值范圍進行預先設定，然后在取值范圍內(nèi)進行隨機取樣并形成參數(shù)值組合輸入模型進行運算。以似然目標函數(shù)為判定依據(jù)，通過對似然函數(shù)值及其權重進行求解最終可得到相應組合條件下各參數(shù)的似然值。似然值邊界條件應首先選取某一值并以此作為邊界值進行判定，低于此值則似然值被賦為0，若高于該值則參數(shù)組合全部被重新歸一化處理，然后在某置信條件下按處理后值的大小進行不確定性區(qū)間的預測和求解。結(jié)合GLUE基本理論和方法，可將其程序應用劃分為以下幾個過程。

步驟一：定義似然目標函數(shù)。設定似然目標函數(shù)用以表征模擬與實測之間的吻合程度，若系統(tǒng)與預測表現(xiàn)出不相似時則對函數(shù)賦值為0，若二者保持良好的相似性且不斷增大則函數(shù)值應滿足線性增大的變化趨勢。

步驟二：確定先驗分布形式即參數(shù)的取值范圍。利用參數(shù)的物理特性確定取值范圍，并且考慮到先驗分布特征在一般情況下較難確定，通常情況可選取符合均勻分布條件的參數(shù)作為先驗分布特征。

步驟三：依據(jù)相關準則進行不確定性分析。利用參數(shù)組模擬狀況確定有效參數(shù)組并描繪似然值與模型參數(shù)的散點圖，以此對不確定性進行定性分析。目標函數(shù)受似然分布影響不明顯時，其似然度則表現(xiàn)出與初始狀態(tài)存在明顯差異的特征，在此條件下參數(shù)具有顯著的敏感性特性。

步驟四：確定概率密度分布。設定臨界值并運行程序，對高于預先設定值的所有似然值進行歸一化計算，并因此得到0～1范圍分布的值，對其進行加權求解得到似然判據(jù)值，利用權重系數(shù)分布條件對空間分布狀況進行求解和分析。

步驟五：預報不確定性范圍?？衫盟迫恢档拇笮∨判蚪Y(jié)果預報一定置信水平的不確定性范圍，通常條件下不確定性范圍的上、下界限為累計似然分布的5%和95%。

表2 SWMM模型各參數(shù)取值范圍

1.3 參數(shù)的確定

水文模型眾多參數(shù)的取值變化可對模擬結(jié)果產(chǎn)生不同程度的影響，模擬結(jié)果的置信水平、不確定性與取值范圍的大小具有顯著相關性。依據(jù)不同參數(shù)的內(nèi)涵，有些參數(shù)可通過實際測量求解而部分參數(shù)只能依據(jù)經(jīng)驗取值無法根據(jù)已知情況確定，因此參數(shù)的取值可對評價結(jié)果的不確定性造成顯著影響。本研究結(jié)合相關文獻以及SWMM模型文獻選取了11個參數(shù)進行不確定性分析。結(jié)合研究區(qū)域特征對參數(shù)取值，并依據(jù)實際物理意義將參數(shù)設定在合理的范圍內(nèi)，本文參考SWMM模型水文參數(shù)取值范圍相關文獻和實測數(shù)據(jù)確定了各參數(shù)的不確定性取值范圍，見表2。表中Imperv參數(shù)物理意義為不透水率，Zero-imperv、N-imperv、N-perv、N-conduit、N-river、D-imperv、D-perv、MaxRate、MinRate以及Decay其物理意義分別為無洼蓄水不透水區(qū)百分比、不透水區(qū)曼寧系數(shù)、透水區(qū)曼寧系數(shù)、管道曼寧系數(shù)、河道曼寧系數(shù)、不透水區(qū)洼蓄量、透水區(qū)洼蓄水量、最大和最小入滲率、衰減系數(shù)。

1.4 不確定性區(qū)間評價指標

針對不確定性區(qū)間估計問題，水文學家給出了一系列優(yōu)良性區(qū)間估計方法，主要是對寬度和覆蓋率等不確定性區(qū)間參數(shù)進行研究。熊立華等在此基礎上針對不確定性區(qū)間的對稱性原則給出了評價指標。本文結(jié)合相關問題分別選取了評價偏移幅度(D)、覆蓋率(CR)、平均相對寬度(RB)以及平均對稱度(S)4個評價指標對不確定性、有效性進行分析，各指標計算公式如下：

(1)

(2)

(3)

(4)

式中，n—流量過程序列點個數(shù)；Qs,upper,Qs,lower—分別代表模擬徑流點估計置信區(qū)間的上下邊界值；Qo、nQin—分別代表實測流量和實測值落在置信區(qū)間的個數(shù)。

在一定的置信水平條件下，預測不確定性區(qū)間的優(yōu)良性隨區(qū)間的評價相對寬度呈顯著的正相關性即寬度越窄則代表優(yōu)良性越好；預測區(qū)間包含實測流量的比例通常采用覆蓋率表示為百分比形式，模擬值接近于實測值的可能性由覆蓋率的高低表征，若覆蓋率趨近于100%則代表其越可能貼近實測值，并且模型的模擬效果越優(yōu)。S值在理想條件或絕對對稱的條件下為0，若S值小于0.5則代表預測流量過程線在平均水平上位于不確定性范圍之內(nèi)，若S值大于等于0.5則表示預測過程線在平均水平上位于預測不確定性區(qū)間之上或之下，此時S值越大則表示不確定區(qū)間性能越差?？紤]到高度非線性的水文現(xiàn)象，觀測流量過程線與不確定性預測區(qū)間的中心線通常是不重合的，因此為了量化這一偏離程度采用平均偏移幅度指標D進行表征，原則上預測區(qū)間的對稱性隨平均偏移幅度的降低而逐漸變好，其不確定性區(qū)間的性能也就越優(yōu)良。

2 結(jié)果與分析

采用隨機采樣法隨機提取1萬組參數(shù)運行模型，其中似然判據(jù)選取確定性系數(shù)并以0.7作為判據(jù)臨界值，由此可對各場次暴雨有效參數(shù)組進行模擬分析。研究表明，不透水區(qū)域曼寧系數(shù)和不透水區(qū)域洼蓄量存在一定范圍的變化但變化幅度不大，而河道與管道曼寧系數(shù)具有較為顯著的高峰范圍，而其影響范圍較低，其他各參數(shù)在相應區(qū)間內(nèi)的變化對似然值的影響較低但存在較為顯著的不確定性；結(jié)合參數(shù)后驗概率分布可知，除河道與管道曼寧系數(shù)之外其他各參數(shù)均接近于均勻分布，此研究結(jié)果與相關研究保持良好的一致性。

在進行參數(shù)不確定分析時似然判據(jù)臨界值取為0.70，當參數(shù)組似然值低于該臨界值時為0，若高于該臨界值則被重新歸一化。對歸一化處理后的似然值進行由大到小的排列并求解某置信水平下模型預報的不確定范圍。本研究利用90%置信區(qū)間分析模擬流量過程，并模擬分析研究區(qū)域在相應期間內(nèi)的各場次暴雨徑流過程。擬合結(jié)果顯示，暴雨流量過程擬合曲線具有良好的吻合度，各場次洪水峰現(xiàn)時間模擬誤差均在2h以內(nèi)，而個別場次的峰值模擬效果較差如20150724場次，模型模擬結(jié)果與峰現(xiàn)時間模擬效果較好的場次為20160718場次和20150625場次。依據(jù)最優(yōu)模擬過程和各場次50%分位點可以看出，最優(yōu)模擬效果與50%分位點流量模擬效果大致相同，甚至可出現(xiàn)峰值處的模擬效果優(yōu)于最優(yōu)模擬過程的現(xiàn)象，因此采用本模型可為峰現(xiàn)預報提供一定的依據(jù)和參考。

本文對沈陽市2015—2016年各場次區(qū)間估計的評價指標值見表3。由表3可知，不確定性模擬區(qū)間寬度在6場暴雨中均處于較低水平，而D值普遍較大。對比分析CR與RB值可知，二者具有明顯的線性相關性，此研究結(jié)果與熊立華等相關結(jié)論保持一致。結(jié)合4項指標值可知，20150721場次與20160622場次的模擬結(jié)果相對較好均符合不確定性分析要求，而20150625場次雖然具有較低的RB和CR值，但具有相對較好的D值；20150724場次以及20160718場次具有相對較好的CR值和D值，其不確定性區(qū)間的性能整體較低。

表3 不確定性區(qū)間評價結(jié)果

3 結(jié)論

本文在詳細分析雨洪預測模型應用現(xiàn)狀的基礎上，通過分析探討GLUE方法理論構(gòu)建了雨洪預測模型，并以沈陽市2015—2016年的雨洪為例利用模型進行預測分析，對不確定性取值范圍進行了簡要概述，結(jié)論如下。

(1)90%置信度水平的GLUE不確定性范圍分析方法雖然不能完全對5場洪水的實測流量全過程進行完全模擬，但是對洪水流量過程、峰值以及峰現(xiàn)之間模擬得很好，并能夠給出較高精度的不確定性范圍。

(2)依據(jù)評價指標的具體內(nèi)涵可以將11個參數(shù)劃分為不敏感和敏感2大類參數(shù)。本研究成果有利于提高雨洪模擬的效率和準確性，可為沈陽市城市雨洪模型的深入研究和發(fā)展提供一定的理論依據(jù)和技術支持。

(3)影響沈陽市洪水模擬結(jié)果的主要因素為河道與城市管道的曼寧系數(shù)，所以研究區(qū)域雨洪過程受排水通道的糙率影響較為顯著。為有效防止洪澇災害的發(fā)生并確保城市排水的通暢，應加強對城市排水管道的日常疏通與維護工作。