趙信峰　徐鵬　劉開磊　趙麗霞　徐十鋒　郟建

摘要：水文模型參數(shù)的選取通常依靠經(jīng)驗判斷或者依賴歷史庫中的不完備數(shù)據(jù)集進行自動優(yōu)選，所選參數(shù)并不一定能夠準確反映流域降雨徑流特點，更不足以反映不同洪水漲落階段洪水特征的變化。基于水文模型參數(shù)存在顯著不確定性的客觀事實，以隨機參數(shù)驅(qū)動水文模型，并結(jié)合數(shù)值模型實現(xiàn)概率預(yù)報。通過東灣流域36場洪水模擬試驗，揭示了水文參數(shù)不確定性對洪水預(yù)報結(jié)果的顯著影響，提出的概率預(yù)報算法能夠給出精確、可靠的預(yù)報結(jié)果，說明該算法能夠降低水文模型參數(shù)所帶來的洪水預(yù)報不確定性。

關(guān)鍵詞：洪水預(yù)報；新安江模型；參數(shù)不確定性；概率分布；概率預(yù)報

受參數(shù)不確定性影響，洪水預(yù)報模型往往難以達到足夠的精度，在實際洪水預(yù)報中基于傳統(tǒng)的洪水預(yù)報模型得到的預(yù)報結(jié)果，其不確定性程度較高，難以據(jù)此作出合適的防汛調(diào)度決策，在實際應(yīng)用中往往通過校正或概率預(yù)報的方式來降低洪水預(yù)報的不確定性。本研究從模型參數(shù)存在不確定性的客觀現(xiàn)實出發(fā)，分析參數(shù)的概率分布特征以生成隨機參數(shù)簇，驅(qū)動洪水預(yù)報模型產(chǎn)生初始預(yù)報解集，采用適合的方法對預(yù)報結(jié)果進行綜合，得到概率預(yù)報結(jié)果。

1模型介紹

1.1新安江模型及其參數(shù)

新安江模型是1983年由河海大學(xué)趙人俊教授帶領(lǐng)水文預(yù)報教研室研制，并逐步完善起來的流域產(chǎn)匯流模型。該模型在國內(nèi)水文預(yù)報中得到廣泛應(yīng)用，并在國際水文學(xué)研究中取得有價值的成果。

新安江模型是典型的概念性模型，由蒸散發(fā)、產(chǎn)流、分水源和匯流4個模塊組成。模型參數(shù)較多，一般需要根據(jù)優(yōu)化目標函數(shù)對參數(shù)的敏感性，將參數(shù)劃分為敏感、不敏感、區(qū)域敏感三類。對不敏感參數(shù)一般取經(jīng)驗值，不參與優(yōu)選；對于敏感或區(qū)域敏感參數(shù)，則需要充分考慮參數(shù)水文特性，采用客觀優(yōu)選或SCE-UA等自動優(yōu)化方法確定參數(shù)值。

河網(wǎng)蓄水消退系數(shù)CS常被歸為敏感參數(shù)，對洪水預(yù)報結(jié)果的影響較大，目前針對其水文特性及統(tǒng)計規(guī)律的研究較多，成果也較為豐富。根據(jù)李致家在沙埠流域?qū)S參數(shù)規(guī)律的研究成果可知，CS是時段長度和線性水庫的蓄泄系數(shù)的函數(shù)，反映流域匯流特性及線性水庫的時問尺度變化。陸曼皎嘗試通過蓄泄系數(shù)參數(shù)規(guī)律來間接推求CS值，其模擬試驗在皖南山區(qū)面積為100-3 000 km2的13個流域進行，結(jié)果表明地理因子公式推求得到CS的方法具有一定的可操作性，同時驗證了計算步長、時段內(nèi)人流分布可能帶來的參數(shù)不確定性。

1，2 BMA模型

BMA（Bayesian Model Averaging）是基于多預(yù)報序列的先驗信息進行模型綜合的數(shù)值模型。DUAN等在美國3個水文站進行模擬試驗，研究BMA方法的集合預(yù)報性能，其研究表明BMA對降低洪水預(yù)報結(jié)果的不確定性程度效果顯著。BMA模型依賴較長系列預(yù)報結(jié)果進行模型訓(xùn)練，進而根據(jù)各個時刻的初始（先驗）預(yù)報結(jié)果估計預(yù)報變量的后驗概率特征，生成概率預(yù)報以及均值、中位數(shù)形式的確定性預(yù)報結(jié)果。

在采用BMA模型進行集合預(yù)報時，無須關(guān)注各模型中哪個是最優(yōu)模型，也無須要求各模型均能夠提供高精度的預(yù)報結(jié)果，只須提供序列足夠長的資料以進行BMA模型訓(xùn)練。因此，在實際使用中，BMA模型能夠避免因?qū)ψ顑?yōu)模型的判斷不準而帶來的不確定性，避免產(chǎn)生較差的預(yù)報結(jié)果，同時能夠提供較為可靠的預(yù)報變量概率分布描述。

2基于參數(shù)不確定性的概率預(yù)報算法

實際上，BMA模型的先驗信息并不局限于多模型預(yù)報結(jié)果，只要給出多個時間序列的原始預(yù)報數(shù)據(jù)及相應(yīng)的實測序列，BMA就可以正常執(zhí)行運算?？紤]到在執(zhí)行實時洪水預(yù)報時，最優(yōu)參數(shù)并不能提前預(yù)知，因此可以依據(jù)參數(shù)的先驗概率分布特征，隨機給出某參數(shù)多個可能的值，以驅(qū)動水文模型產(chǎn)生相應(yīng)多個預(yù)報結(jié)果。在BMA框架下，把隨機參數(shù)所產(chǎn)生的多個預(yù)報結(jié)果視作BMA的集合預(yù)報成員，以基于以上各參數(shù)模擬歷史場次洪水的計算結(jié)果作為先驗信息，進行BMA模型訓(xùn)練；進而，以BMA模型綜合當前洪水的多個預(yù)報結(jié)果，得到在考慮參數(shù)不確定性的情境下，預(yù)報變量的概率分布的描述。其中，隨機生成的該參數(shù)的多個可能值被統(tǒng)一稱為參數(shù)簇。

根據(jù)以上思路，以新安江模型為例，考慮其參數(shù)CS的不確定性，采用如下步驟構(gòu)建基于參數(shù)不確定性的概率預(yù)報算法（簡稱PROP）。

（1）獲取參數(shù)的先驗概率分布。根據(jù)經(jīng)驗，選擇以新安江模型的參數(shù)CS為例，考察該參數(shù)在歷史各場洪水中的數(shù)值變化特征，考察各常見分布類型在描述CS的概率分布中的適用性。選出合適的分布類型之后，計算分布函數(shù)的參數(shù)，獲取CS的先驗概率分布。

（2）隨機生成參數(shù)簇。根據(jù)CS的先驗概率分布特征，隨機生成維度為Ⅳ的參數(shù)簇。

（3）構(gòu)建預(yù)報信息庫?；谝陨螻個參數(shù)，分別驅(qū)動新安江模型模擬所有場次的歷史洪水，計算得到各場洪水的次模模擬結(jié)果。在實時洪水預(yù)報中，步驟（2）中得到的參數(shù)簇可以在洪水預(yù)報之前生成，以降低運算量，保證實時性；這里提到的“歷史洪水”應(yīng)當變成“當前場次以前的歷史洪水”。

（4）訓(xùn)練BMA模型。根據(jù)成員數(shù)為N的歷史洪水預(yù)報結(jié)果的集合，訓(xùn)練BMA模型參數(shù)。受限于篇幅，此處不對相關(guān)技術(shù)細節(jié)作詳細描述，可在參考文獻[7]中查閱。

（5）生成預(yù)報變量后驗概率分布。設(shè)后驗分布的采樣數(shù)目為L，然后將當前的N個預(yù)報結(jié)果代入訓(xùn)練好的BMA模型中，基于蒙特卡羅采樣方法生成成員數(shù)為L的預(yù)報變量的解集。當L值足夠大時，該解集與預(yù)報變量的后驗概率分布相似，可以認為該解集的分布情況反映了預(yù)報變量的后驗概率分布特征，解集的均值可以視作預(yù)報變量的期望值。

在PROP算法中，模型參數(shù)的最優(yōu)值無須提前預(yù)知，因此該算法能夠避免洪水預(yù)報中不合理的參數(shù)對預(yù)報結(jié)果的負面影響。該算法僅依靠比較成熟、單一的新安江模型即可實現(xiàn)集合預(yù)報，無須引進其他模型，算法的實現(xiàn)簡便。該算法能夠提供變量的后驗概率分布及期望值等信息，比傳統(tǒng)水文模型的信息更豐富，能夠為防汛決策提供更多有價值的支持信息。

3試驗流域及數(shù)據(jù)介紹

3.1流域概況

本研究所選試驗流域為東灣流域（見圖1），位于東經(jīng)111°-112°、北緯33.0°-34.5°之間的伊河河源地區(qū)，流域面積為2 856 km2。流域地勢西高東低，上游林地面積大，屬大陸性季風氣候區(qū)。降水量的年內(nèi)分布極不均勻，每年7-9月降水量占年降水量的一半以上。年降水量隨高程增加而遞增，山地為多雨區(qū)，河谷及附近丘陵為少雨區(qū)。降水年際變化較大，最大年降水量為最小年降水量的2倍左右。

3.2流域數(shù)據(jù)

本研究選擇東灣流域1962-2011年36場洪水用于數(shù)值模擬試驗。模型參數(shù)率定中，采用先優(yōu)化日模參數(shù)，再模擬出各場次洪水初始時刻的流域土壤飽和程度，然后排列出敏感次模參數(shù)的順序。由于本研究所關(guān)注的是參數(shù)的不確定性程度及基于參數(shù)不確定性的相關(guān)理論方法，每場洪水的最優(yōu)參數(shù)值不同，這就需要知道各場洪水中的最優(yōu)參數(shù)值，因此每場洪水都需要優(yōu)化得到一個最優(yōu)參數(shù)值。

根據(jù)《水文情報預(yù)報規(guī)范》（GB/T 22482-2008），本研究中涉及的洪水預(yù)報精度評定，均采用NSE（Nash-Sutcliffe Efficiency Coefficient）指標作為評價依據(jù)。以每場洪水的最優(yōu)參數(shù)為依據(jù)，分別執(zhí)行預(yù)報并計算NSE指標值。統(tǒng)計結(jié)果顯示，所選各場洪水均達到丙級以上精度，其中乙級以上32場，甲級以上11場，說明新安江模型在東灣流域適用性好，本研究的成果對于推進新安江模型在該流域及相似流域的應(yīng)用具有一定的參考價值。

4參數(shù)CS的先驗概率分布

本研究采用SCE-UA算法。分別對每場洪水進行參數(shù)率定，以獲取CS在每場洪水的最優(yōu)值。計算率定得到36個CS參數(shù)值的累計概率值，并應(yīng)用常見的Beta、Gamma、Laplace、Logistic、P-Ⅲ、Weibull（兩參數(shù)）、對數(shù)正態(tài)、正態(tài)、指數(shù)分布共9種分布函數(shù)對CS值的統(tǒng)計概率分布點進行擬合。擬合情況見圖2，圖中橫坐標表示直接統(tǒng)計最優(yōu)CS值的頻率；縱坐標表示相應(yīng)CS值在不同分布類型中的頻率估計值；NSE為擬合確定性系數(shù)，其數(shù)值范圍為（-∞，1]，數(shù)值越大，則擬合度越高，相應(yīng)的分布函數(shù)越適合于描述最優(yōu)CS值的概率分布特征。

從圖2可以看出，所選9個分布函數(shù)的確定性系數(shù)值均在0.9以上，各分布函數(shù)對最優(yōu)CS值的擬合度均較高。Beta、P-Ⅲ相對于其他分布函數(shù)表現(xiàn)更優(yōu)，然而P-Ⅲ分布對較大CS值的擬合效果較差，因而可以認為Beta分布是這9個分布函數(shù)中最好的，適合于描述不同場次洪水中最優(yōu)CS值的概率分布特征。在確定參數(shù)CS所服從分布類型之后，記錄其分布函數(shù)的各項參數(shù)值，隨機生成50組（L=50）不同的CS參數(shù)。

5原始預(yù)報集合的獲得

將隨機生成的50組參數(shù)分別代入新安江模型，獨立地對所選歷史洪水進行模擬計算，可以獲得集合預(yù)報成員數(shù)為50的預(yù)報值序列，然后結(jié)合對應(yīng)的實測數(shù)據(jù)對BMA模型進行訓(xùn)練。針對每場洪水，利用訓(xùn)練好的BMA模型綜合50個原始預(yù)報結(jié)果進行集合預(yù)報運算。圖3、圖4分別展示了實測一原始預(yù)報一PROP均值預(yù)報流量過程比較、PROP概率預(yù)報流量過程。限于篇幅，在圖3、圖4中將36場洪水彼此首尾相連，繪制在同一張圖上。由于圖中相鄰兩時刻的時間跨度可能比較大，因此該洪水過程并不代表1962-2011年完整的流量過程，僅用于說明各集合預(yù)報成員及PROP在這36場洪水中的表現(xiàn)。

從圖3可以看到，基于隨機參數(shù)所構(gòu)建的各集合預(yù)報成員，模擬得到的流量過程的總體趨勢與實測的相差不大，但是仍然可見大量的原始預(yù)報結(jié)果偏離實測值。各集合預(yù)報成員所選CS參數(shù)值不同，是各自預(yù)報結(jié)果之間產(chǎn)生顯著差異的唯一原因。這進一步印證了參數(shù)不確定性的存在會導(dǎo)致預(yù)報結(jié)果顯著發(fā)散的現(xiàn)象，也強調(diào)了參數(shù)不確定性是水文預(yù)報中必須考慮的一個重要因素。PROP算法均值預(yù)報結(jié)果與實測流量結(jié)果的總體變化趨勢一致，說明基于PROP方法進行洪水預(yù)報時，其確定性預(yù)報結(jié)果更為穩(wěn)定、可靠，且無須事先進行參數(shù)率定，因而可以避免參數(shù)不確定性導(dǎo)致的計算精度降低。

6概率預(yù)報結(jié)果的獲得與分析

概率預(yù)報時常采用90%置信區(qū)問來表征預(yù)報變量的概率分布情況，置信區(qū)間中的實測點占全部實測點的比例（即覆蓋率CR）被認為能夠反映概率預(yù)報結(jié)果的可靠性。CR的取值范圍為[0，1]，CR值越大，實測點落在置信區(qū)問內(nèi)的比例越大，依據(jù)概率預(yù)報結(jié)果漏報的可能性越小。

從圖4可以看出，PROP所提供的90%置信區(qū)間基本上能夠?qū)崪y值包含在內(nèi)，CR值達到93.0%。PROP算法的洪峰預(yù)報結(jié)果可靠性高，所選36場洪水中有31場的實測洪峰落在概率預(yù)報的置信區(qū)間之內(nèi)，例如在1975年大洪水中，實測洪峰流量為4 200 m3/s，對應(yīng)的PROP預(yù)報90%上、下限分別為4 430 ma/s與490 m3/s。綜合PROP的概率預(yù)報結(jié)果對整體與洪峰附近洪水過程的匹配程度，可以認為該模型所獲得的概率預(yù)報結(jié)果較為可靠。

7結(jié)論

本研究基于水文模型的參數(shù)存在顯著不確定性的客觀事實，以隨機生成的參數(shù)驅(qū)動水文模型，并結(jié)合數(shù)值模型構(gòu)建PROP算法實現(xiàn)集合預(yù)報。通過東灣流域36場洪水模擬試驗，揭示了水文模型參數(shù)不確定性對洪水預(yù)報結(jié)果的顯著影響，并驗證了PROP所提供的確定性及概率預(yù)報結(jié)果的精確性、可靠性，證明PROP能夠降低水文模型參數(shù)所帶來的洪水預(yù)報不確定性。

在實際洪水預(yù)報中，參數(shù)的優(yōu)選往往依靠經(jīng)驗判斷或者依賴歷史庫中的不完備數(shù)據(jù)集進行自動優(yōu)化，然而由于洪水特征無法準確預(yù)知，甚至在一場洪水不同階段所要求參數(shù)的值也存在較大差異，因此所選參數(shù)并不一定適合當前洪水的預(yù)報。PROP算法為考慮參數(shù)不確定性的條件下實現(xiàn)準確的洪水預(yù)報問題，提供了一個可靠的解決方案。該算法強化了對參數(shù)概率分布特征的描述，弱化了對求解最優(yōu)參數(shù)值的要求，降低了參數(shù)不確定性導(dǎo)致較差預(yù)報結(jié)果的可能性：依賴現(xiàn)有較為成熟的BMA模型，為洪水預(yù)報工作提供更為豐富、可靠的預(yù)報信息，對于完善并提高現(xiàn)有的洪水預(yù)報技術(shù)具有參考價值。

然而也需要看到，本研究所采用的參數(shù)僅僅針對單一的參數(shù)CS，實際洪水預(yù)報中往往有多個參數(shù)存在明顯不確定性，如何準確描述多參數(shù)的聯(lián)合概率分布，以及如何在洪水預(yù)報中同時考慮多參數(shù)不確定性的影響，將是一個有價值的研究方向。