張 潔，梁忠民，胡義明，王 軍，李彬權(quán)

(河海大學(xué)水文水資源學(xué)院，南京 210098)

洪水資料的一致性是現(xiàn)行水文頻率分析計(jì)算的基本前提[1]。近年來(lái)，人類(lèi)活動(dòng)和氣候變化使得洪水形成的下墊面和氣候條件發(fā)生了重大改變，水文極值系列不再具有一致性，那么一致性條件下的水文頻率分析方法在變化環(huán)境下也就不再適用[2-7]。因此，變化環(huán)境下非一致性水文頻率分析方法的研究成為廣大學(xué)者的研究重點(diǎn)，其中，基于變參數(shù)概率分布模型途徑的非一致性水文頻率分析是目前研究最為廣泛的方法之一[8-10]。其通過(guò)分析不同變參數(shù)概率分布模型的擬合效果，選取擬合最優(yōu)的模型作為最終使用模型。在此基礎(chǔ)上，采用期望等待時(shí)間法、期望發(fā)生次數(shù)法、設(shè)計(jì)壽命水平法[9]或等可靠度法[10]推求變化環(huán)境下特定重現(xiàn)期對(duì)應(yīng)的非一致性水文設(shè)計(jì)值。水文設(shè)計(jì)值推求的可靠性受諸多因素影響，如模型類(lèi)型、參數(shù)估計(jì)等。胡義明等研究了變參數(shù)概率分布函數(shù)模型中參數(shù)估計(jì)不確定性、模型輸入不確定性以及2者耦合對(duì)水文設(shè)計(jì)值估計(jì)的影響[9]。變化環(huán)境下非平穩(wěn)的水文極值系列理論上的變參數(shù)概率分布模型無(wú)法得知，這使得依據(jù)擬合最優(yōu)原則獲得的最優(yōu)模型也未必能精準(zhǔn)地代表非平穩(wěn)極值系列的概率分布特征，導(dǎo)致依據(jù)單一最優(yōu)擬合模型計(jì)算的結(jié)果不可避免地存在不確定性。

為此，本文選取赤池信息準(zhǔn)則(AIC)評(píng)估了多個(gè)變參數(shù)概率分布模型的擬合效果，提出將擬合效果較好的若干種變參數(shù)概率分布函數(shù)模型進(jìn)行綜合，以減小模型選擇對(duì)水文設(shè)計(jì)值產(chǎn)生的不確定性影響，增強(qiáng)水文頻率分析計(jì)算得到的水文設(shè)計(jì)值的可靠性。

1 方法原理

1.1 “等可靠度”法基本原理

下面介紹“等可靠度”法。假定工程設(shè)計(jì)壽命期為L(zhǎng)，且n年的歷史洪峰資料x(chóng)1，x2，…,xn存在減少趨勢(shì)；采用“等可靠度”法推求重現(xiàn)期為T(mén)的設(shè)計(jì)洪水的基本流程可表述如下[9]。

首先采用分布參數(shù)隨時(shí)間等協(xié)變量的變化而變化的變參數(shù)PE3概率分布函數(shù)描述該洪峰系列。假定尺度參數(shù)β和形狀參數(shù)γ為不隨時(shí)間變化的常數(shù)，而位置參數(shù)μ隨時(shí)間線性變化，即：

μ(t)=μ0+μ1t

(1)

β=β0

(2)

γ=γ0

(3)

此時(shí)，變參數(shù)的PE3概率密度函數(shù)可表示為：

(4)

變參數(shù)PE3概率分布函數(shù)為Ft[μ(t),β,γ],t=1,2,…,n，其中的參數(shù)μ0、μ1、β和γ可由n年的歷史洪峰數(shù)據(jù)x1，x2，…,xn估計(jì)得到，進(jìn)而推求工程壽命期內(nèi)(n+1～n+L)逐年的變參數(shù)PE3概率分布函數(shù)Ft[μ(t),β,γ],t=n+1,n+2,…,n+L。則變化環(huán)境下水文設(shè)計(jì)值XT,NS的水文可靠度為：

(5)

(6)

(7)

根據(jù)“等可靠度”法，變化環(huán)境的非一致性條件下不再具有重現(xiàn)期的概念，但非一致性條件下計(jì)算的水文設(shè)計(jì)值所具有的可靠度至少應(yīng)該達(dá)到一致性條件下計(jì)算的設(shè)計(jì)值的可靠度。令：

RNS=RS

(8)

求解式(8)，即可得到變化環(huán)境下非一致性條件的水文設(shè)計(jì)值XT,NS。

1.2 考慮模型選擇不確定性的設(shè)計(jì)值計(jì)算

采用變參數(shù)概率分布模型進(jìn)行變化環(huán)境下非一致性水文頻率分析時(shí)，由于人類(lèi)活動(dòng)和氣候變化對(duì)水文極值系列分布規(guī)律產(chǎn)生了很大影響，因此將分布函數(shù)中的參數(shù)設(shè)定為隨時(shí)間等協(xié)變量而改變的變量，而不再采用常數(shù)。然而，概率分布中參數(shù)隨協(xié)變量變化的準(zhǔn)確關(guān)系是未知的，現(xiàn)階段擬合效果最好的模型并不能保證在未來(lái)也始終優(yōu)于其他模型，隨著時(shí)間的外延，水文設(shè)計(jì)值估計(jì)的可靠性越來(lái)越低，而不確定性隨之增大。為降低模型選擇導(dǎo)致的不確定性，增強(qiáng)水文頻率分析計(jì)算結(jié)果的可靠性，可綜合考慮多個(gè)擬合效果較好的變參數(shù)概率分布模型，其中，各模型的權(quán)重可由其擬合程度確定。下面闡述推求設(shè)計(jì)重現(xiàn)期為T(mén)，工程設(shè)計(jì)壽命為L(zhǎng)的水文設(shè)計(jì)值的具體過(guò)程。

(1)以不同的概率分布函數(shù)(如PE3概率分布函數(shù)、廣義極值分布函數(shù)等)為基礎(chǔ)，構(gòu)建不同的變參數(shù)概率分布函數(shù)模型來(lái)擬合變化環(huán)境下非平穩(wěn)的水文極值系列。假定概率分布函數(shù)中位置參數(shù)、尺度參數(shù)、形狀參數(shù)等參數(shù)隨時(shí)間這一協(xié)變量的變化而變化，采用AIC指標(biāo)評(píng)估各模型擬合效果，并根據(jù)擬合效果選出擬合效果最優(yōu)的前m個(gè)模型，記為Mi，i=1,2，…，m。

(2)基于“等可靠度”法計(jì)算步驟(1)中第i個(gè)模型Mi對(duì)應(yīng)的水文設(shè)計(jì)值X(Mi)，i=1,2，…，m：

(9)

(3)采用AIC指標(biāo)對(duì)擬合效果進(jìn)行評(píng)估。AIC值越小的模型擬合效果越好，其所占權(quán)重也應(yīng)該越大。根據(jù)評(píng)估結(jié)果，計(jì)算模型Mi的權(quán)重wi，i=1,2，…，m；

(10)

式中：AICi為第i個(gè)模型Mi的AIC值。

(4)將m個(gè)模型計(jì)算的水文設(shè)計(jì)值X(Mi)加權(quán)平均，獲得綜合水文設(shè)計(jì)值X，進(jìn)而降低模型選擇的不確定性影響。即：

(11)

值得注意的是，公式(10)僅適用于經(jīng)AIC評(píng)估后2個(gè)或2個(gè)以上模型的AIC值接近且明顯小于其他模型的情況。權(quán)重的計(jì)算僅限于擬合效果較好的模型，而不能應(yīng)用于建立的全部模型。當(dāng)某個(gè)模型的AIC值明顯小于其他模型時(shí)，可直接采用此最優(yōu)模型進(jìn)行水文頻率分析計(jì)算。AIC值作為評(píng)估模型擬合效果的指標(biāo)，具有對(duì)各模型篩選比較的功能，AIC值明顯較高的模型擬合效果差，應(yīng)該及時(shí)去除，否則不能發(fā)揮AIC指標(biāo)的功能，將大大增加水文設(shè)計(jì)值的誤差。公式(10)的作用是在AIC指標(biāo)排除較差模型的基礎(chǔ)上，將篩選得到的優(yōu)秀模型進(jìn)行綜合，以進(jìn)一步增強(qiáng)水文頻率分析計(jì)算的可靠性。

2 實(shí)例研究

本文研究對(duì)象為黃龍灘站1956-2014年共59 a的7 d洪量歷史觀測(cè)資料，黃龍灘站所在流域的位置、水系和水庫(kù)分布等情況[4]見(jiàn)圖1。圖2為7 d洪量系列的時(shí)間序列圖。由圖2可知，7 d洪量系列明顯地隨時(shí)間而減小。

圖1 黃龍灘所在流域基本情況Fig.1 Basic information map of Huanglongtan river basin

圖2 7 d洪量時(shí)間序列Fig.2 Time series of 7 day flood volume

選定時(shí)間作為協(xié)變量，通過(guò)假設(shè)PE3和GEV分布函數(shù)的位置參數(shù)或尺度參數(shù)并非常數(shù)而是隨協(xié)變量變化的變量，建立下面4個(gè)不同的變參數(shù)概率分布模型。每個(gè)模型中參數(shù)與協(xié)變量的關(guān)系如下。

(1)模型1：PE3分布函數(shù)中的尺度參數(shù)β和形狀參數(shù)γ為不隨時(shí)間變化的常數(shù)，而位置參數(shù)μ是隨時(shí)間線性變化的變量，記為PE3-Loc，即：

μ(t)=μ0+μ1t,β(t)=β0,γ(t)=γ0

(12)

(2)模型2：PE3分布函數(shù)中的形狀參數(shù)γ為不隨時(shí)間變化的常數(shù)，而位置參數(shù)μ和尺度參數(shù)β是隨時(shí)間線性變化的變量，記為PE3-Loc-Scl，即：

μ(t)=μ0+μ1t,β(t)=exp(β0+β1t),γ(t)=γ0

(13)

(3)模型3：GEV分布函數(shù)中的尺度參數(shù)α和形狀參數(shù)γ為不隨時(shí)間變化的常數(shù)，而位置參數(shù)ξ是隨時(shí)間線性變化的變量，記為GEV-Loc，即：

ξ(t)=ξ0+ξ1t,α(t)=α0,γ(t)=γ0

(14)

(4)模型4：GEV分布函數(shù)中的形狀參數(shù)γ為不隨時(shí)間變化的常數(shù)，而位置參數(shù)ξ和尺度參數(shù)α是隨時(shí)間線性變化的變量，記為GEV-Loc-Scl，即：

ξ(t)=ξ0+ξ1t,α(t)=exp(α0+α1t),γ(t)=γ0

(15)

采用貝葉斯方法對(duì)模型參數(shù)進(jìn)行估計(jì)，模型的估計(jì)參數(shù)為最大后驗(yàn)概率值對(duì)應(yīng)的參數(shù)。平行運(yùn)行5條鏈，在經(jīng)過(guò)1 000次采樣后，各個(gè)模型的采樣樣本已收斂。每條鏈上除去用于預(yù)熱的前面900個(gè)樣本后余下100個(gè)參數(shù)，總共余下500個(gè)參數(shù)組，而參數(shù)的最大后驗(yàn)估計(jì)就是得到最大參數(shù)后驗(yàn)密度值的參數(shù)。計(jì)算各個(gè)模型的AIC值，以此作為指標(biāo)評(píng)估各模型的擬合效果，見(jiàn)表1。

表1 不同模型擬合的AIC指標(biāo)值Tab.1 AIC for assessing the model performance

由表1可知，AIC指標(biāo)值最小的模型是模型2即PE3-Loc-Scl模型，由此得到此模型具有最好的擬合效果。模型4即GEV-Loc-Scl次之，且2模型的AIC值較為接近。直接采用PE3-Loc-Scl模型進(jìn)行7 d洪量設(shè)計(jì)值計(jì)算可能存在誤差，故綜合考慮PE3-Loc-Scl和GEV-Loc-Scl 2種模型，以降低模型選擇對(duì)水文設(shè)計(jì)值產(chǎn)生的不確定性影響，增強(qiáng)水文頻率分析計(jì)算結(jié)果的可靠性。

分別將PE3-Loc-Scl和GEV-Loc-Scl的 AIC指標(biāo)值的倒數(shù)占2種模型AIC指標(biāo)值倒數(shù)之和的比重作為其權(quán)重，見(jiàn)表2。

表2 PE3-Loc-Scl和GEV-Loc-Scl的權(quán)重計(jì)算Tab.2 The weights of PE3-Loc-Scl and GEV-Loc-Scl

計(jì)算重現(xiàn)期為100 a時(shí)，PE3-Loc-Scl模型、GEV-Loc-Scl模型和綜合模型計(jì)算的7 d洪量設(shè)計(jì)值隨時(shí)間的變化趨勢(shì)，見(jiàn)圖3。

圖3 PE3-Loc-Scl模型、GEV-Loc-Scl模型和綜合模型100 d一遇7 d洪量設(shè)計(jì)值的時(shí)間變化規(guī)律Fig.3 The change with time of design values of one-hundred-year 7 day flood volume for the PE3-Loc-Scl model、the GEV-Loc-Scl model and the comprehensive model

根據(jù)圖3，PE3-Loc-Scl模型、GEV-Loc-Scl模型和綜合模型計(jì)算的重現(xiàn)期為100 a時(shí)的7 d洪量設(shè)計(jì)值均隨時(shí)間呈減少趨勢(shì)，與7 d洪量實(shí)測(cè)資料的趨勢(shì)相一致。PE3-Loc-Scl模型和GEV-Loc-Scl模型的100 a一遇7 d洪量設(shè)計(jì)值存在差異，2者加權(quán)平均后的綜合模型得到的綜合設(shè)計(jì)值可能具有更高的可靠度。在1956-2014年總共59 a的7 d洪量實(shí)測(cè)資料中，有0次超過(guò)綜合模型100 a一遇7 d洪量設(shè)計(jì)值，表明綜合模型的擬合效果較好。

基于PE3-Loc-Scl模型、GEV-Loc-Scl模型和綜合模型，采用“等可靠度”法計(jì)算設(shè)計(jì)標(biāo)準(zhǔn)為100 a一遇和20 a一遇時(shí)，40、60、80、100 a的工程設(shè)計(jì)壽命對(duì)應(yīng)的7 d洪量設(shè)計(jì)值，結(jié)果見(jiàn)圖4。盡管PE3-Loc-Scl模型和GEV-Loc-Scl模型的擬合效果接近，但根據(jù)2者計(jì)算的7 d洪量設(shè)計(jì)值卻存在較大差異。由于無(wú)法知曉其理論分布，因此，2個(gè)模型都未必能精準(zhǔn)地代表該洪量系列的概率分布特征，故對(duì)2個(gè)模型的計(jì)算結(jié)果進(jìn)行加權(quán)綜合獲得綜合設(shè)計(jì)值，以盡可能地降低模型選取對(duì)設(shè)計(jì)值的影響。

圖4 不同工程設(shè)計(jì)壽命期及重現(xiàn)期下3種模型推求的7 d洪量期望估計(jì)Fig.4 Design value for 7 day flood volume of three models with different return periods under the conditions of different design engineering life

3 結(jié) 論

(1)采用AIC指標(biāo)評(píng)估多個(gè)模型對(duì)7 d洪量系列的擬合效果，并根據(jù)AIC評(píng)估結(jié)果計(jì)算較優(yōu)模型的權(quán)重，綜合不同變參數(shù)概率分布模型的計(jì)算結(jié)果。

(2)較優(yōu)模型具有相近的擬合效果，但推求的設(shè)計(jì)值卻表現(xiàn)出較大差異性。通過(guò)綜合多個(gè)較優(yōu)模型的計(jì)算結(jié)果，一定程度上增強(qiáng)了計(jì)算結(jié)果的可靠性。

(3)當(dāng)采用其他方法(如期望等待時(shí)間法等)推求非一致性條件下水文設(shè)計(jì)值時(shí)，同樣可以采用類(lèi)似方法計(jì)算綜合水文設(shè)計(jì)值以減小模型選擇的不確定性。