• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于貝葉斯理論的洪水頻率不確定性分析

    2016-03-23 05:30:37張冬冬劉少華
    中國農(nóng)村水利水電 2016年4期
    關(guān)鍵詞:后驗置信區(qū)間參數(shù)估計

    張冬冬,魯 帆,劉少華

    (1.長江水利委員會水文局,武漢 430010;2. 中國水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國家重點實驗室,北京 100038 )

    洪水頻率分析主要是利用已有的實測洪水?dāng)?shù)據(jù)構(gòu)造洪水分布函數(shù),通過分布函數(shù)計算一定重現(xiàn)期下的洪水設(shè)計值。洪水頻率分析中的難點在于極值事件的頻率計算,由于極值洪水相對于一般的洪水樣本,發(fā)生頻率較小,而在實際計算的過程中,極值洪水以及可以利用的歷史特大洪水觀測數(shù)據(jù)往往較為缺乏,影響分布函數(shù)參數(shù)的估計和高分位數(shù)計算,從而導(dǎo)致洪水頻率分析中存在較大的不確定性。如何定量的評估洪水頻率計算中線型參數(shù)以及推求設(shè)計值的不確定性,成為洪水頻率分析中的熱點問題。

    本文選取大渡河流域典型站點洪水資料作為分析對象,選取廣義極值分布作為洪水分布的線型,利用貝葉斯原理,通過先驗分布和似然函數(shù)構(gòu)造后驗分布;利用MCMC模擬方法,根據(jù)抽樣值構(gòu)造一定頻率下的洪水設(shè)計值,并得到其相應(yīng)的置信區(qū)間,以定量描述洪水頻率分析的不確定性,并與傳統(tǒng)方法的置信區(qū)間進(jìn)行對比,從而對該方法的有效性和優(yōu)良性進(jìn)行評估。

    1 貝葉斯洪水頻率分析的基本原理

    貝葉斯估計的理論基礎(chǔ)是貝葉斯原理,主要是綜合考慮未知總體分布和樣本分布中的參數(shù)先驗信息,根據(jù)信息求出樣本的后驗分布,再根據(jù)后驗分布推斷未知總體分布的參數(shù)[1]。貝葉斯估計與傳統(tǒng)的極大似然估計法、矩法和權(quán)函數(shù)法等參數(shù)估計方法的主要區(qū)別在于:貝葉斯估計將未知參數(shù)視為隨機(jī)變量,在利用樣本信息推斷總體分布參數(shù)之前,需要規(guī)定樣本參數(shù)的先驗分布,先驗分布的選取通常是利用已有的歷史洪水資料或利用決策者的經(jīng)驗進(jìn)行設(shè)定。在洪水頻率計算中,能夠運(yùn)用一些有效的信息(歷史洪水等),表示成先驗分布的形式,通過最大化的利用有效信息,可以獲得較其他傳統(tǒng)參數(shù)估計方法更精確的推斷。與傳統(tǒng)估計方法相比,貝葉斯估計方法優(yōu)勢在于對模型不確定性的評估,貝葉斯方法利用大量隨機(jī)樣本模擬,概率分布的形式輸出模型的參數(shù)以及一定重現(xiàn)期下的設(shè)計值,概率分布相比與單一估計值包含了額外的信息,可以反映洪水頻率分析中的不確定性。近年來,許多學(xué)者嘗試在水文頻率分析研究中采用貝葉斯的估計方法。Wood等[2]最早提出了基于貝葉斯理論的水文模型參數(shù)估計方法的技術(shù)框架,Kuczera[3]根據(jù)貝葉斯理論同時探索了應(yīng)用分布參數(shù)不確定性問題,主要是利用重要抽樣方法分析了皮爾遜III型和對數(shù)皮爾遜III型曲線的分布參數(shù)。W Femandes等[4]認(rèn)為PMF的不確定性和先驗分布線型有關(guān),針對洪水頻率線型多數(shù)沒有考慮上下界問題,他采用貝葉斯方法和具有上界的分布函數(shù)相結(jié)合,用以福爾瑟姆水庫(美國)為例進(jìn)行洪水頻率分析。后驗統(tǒng)計推斷中分母積分的計算是應(yīng)用貝葉斯方法的一個難點。為了簡化后驗分布中的積分計算,可以通過選擇合理的先驗分布族的方法來實現(xiàn),但某些情況下,數(shù)值積分的方法無法計算參數(shù)較多且結(jié)構(gòu)復(fù)雜的后驗分布。馬爾科夫鏈蒙特卡羅(MCMC)方法能夠解決貝葉斯后驗分布高維積分的問題[5]。

    假設(shè)洪水序列以及分布線型是已知的情況下,利用貝葉斯估計方法進(jìn)行洪水頻率分析的基本框架如圖1所示。

    圖1 貝葉斯參數(shù)估計方法的基本原理Fig.1 Basic principle of bayesian parameters estimation

    (1)先驗分布。依據(jù)一定的原則確定總體分布參數(shù)的先驗分布。在先驗分布的選擇上,目前沒有一個統(tǒng)一的規(guī)范。但是有一些推薦的原則,例如,在沒有歷史數(shù)據(jù)以及經(jīng)驗可以作為參考的情況下,可以根據(jù)參數(shù)的物理意義,采用取值范圍內(nèi)的均勻分布作為總體分布參數(shù)的先驗分布。如果認(rèn)為參數(shù)在取值范圍內(nèi)更傾向取較小的值,推薦應(yīng)用方差較大的正態(tài)分布,相反,可以應(yīng)用方差較小的正態(tài)分布。

    (2)似然函數(shù)。以θ表示總體分布參數(shù)向量,可以用π(θ)表示其先驗分布的密度。由于Xi是相互獨立的,所以樣本x=(X1,X2,…,Xn)的似然函數(shù)f(x|θ)可用下式計算:

    (1)

    (3)后驗分布。由貝葉斯原理,后驗分布可以由先驗密度π(θ)以及式(1)確定的似然函數(shù)計算得到,如下式:

    (2)

    (4)總體分布的參數(shù)估計。根據(jù)式(2)求得總體分布參數(shù)θ的后驗分布,貝葉斯統(tǒng)計推斷與傳統(tǒng)的參數(shù)估計相比,最大的不同是貝葉斯給出的總體分布參數(shù)θ的分布函數(shù)?;诜植己瘮?shù)可以對θ的進(jìn)行統(tǒng)計推斷,一般情況下,θ的估計值有兩種選取方法,第一種方法是將后驗分布的50%分位數(shù)作為 的估計值,第二種方法是將使得后驗密度最大的θ值作為其估計值。

    (5)洪水頻率設(shè)計值計算。假設(shè)已知洪水觀測值x的情況下,θ的后驗分布可以用f(θ|x)來表示,設(shè)定 表示需要計算的一定重現(xiàn)期下的洪水設(shè)計值,那么在已知x的條件下,z的密度函數(shù)可以表示為:

    (3)

    由式(3)得:

    (4)

    Pr(Z≤z|x)可以理解為,已知洪水資料條件下,某一重現(xiàn)期洪水頻率設(shè)計值z的分布函數(shù)。

    如果假設(shè)洪水重現(xiàn)期為m年一遇,即Pr(Z≤z|x)=1-1/m,通過求解方程即可求得設(shè)計值。由于在輸入的時候使用了樣本參數(shù)的先驗分布以及似然函數(shù),因此貝葉斯方法可以將模型參數(shù)的不確定性考慮進(jìn)去,得出的設(shè)計值也是概率分布的形式,體現(xiàn)了由于參數(shù)不確定性引起的頻率設(shè)計的不確定性。

    在求解方程的過程中,由于后驗概率密度較為復(fù)雜,需要借助特定的抽樣模擬方法估計式(4)后驗分布的估計值。

    2 MCMC模擬

    MCMC方法[6]的基本流程如下。

    (1)根據(jù)參數(shù)的后驗分布,通過MCMC方法產(chǎn)生服從式(2)要求的一階馬爾科夫鏈的模擬序列θ0,θ1,θ2,…,其中θ0為任意初始值,θi+1只依賴于當(dāng)前的θi,與之前的序列θ0,θ1,…,θi-1無關(guān),即θi+1由條件分布q(·|θi)產(chǎn)生。

    (2)設(shè)定轉(zhuǎn)移核。為了保證馬爾科夫鏈?zhǔn)菚r間齊次的,即任意參數(shù)經(jīng)過一步迭代以后的邊緣分布與i無關(guān),這里定義轉(zhuǎn)移核為p(θ,θ′)=q(θi+1=θ′|θi=θ)。并使得后驗分布f(θ|x)為平穩(wěn)分布,即要求:

    ?θ′∈Θ

    (5)

    (3)通過MCMC方法產(chǎn)生隨機(jī)樣本。設(shè)定任意一個初始的θ0,由(1)和(2)中定義的馬爾克夫鏈進(jìn)行抽樣模擬,根據(jù)訓(xùn)練步長產(chǎn)生隨機(jī)序列θ1,θ2,…,θn。文獻(xiàn)[7]指出,由于初始設(shè)定θ0的邊緣分布并不一定是f(θ|x),需要一定的訓(xùn)練時期使得序列達(dá)到平穩(wěn),因此在利用MCMC模擬的時候,通常去掉非平穩(wěn)序列的θ1,θ2,…,θk,用剩下的θk+1,θk+2,…,θn作為后驗分布f(θ|x)的抽樣即可。

    (4)利用(3)中產(chǎn)生的序列進(jìn)行蒙特卡洛積分,本文選擇模擬樣本的50%分位數(shù)作為后驗分布的均值,而后驗分布密度由模擬序列的頻率直方圖生成。

    在使用MCMC方法進(jìn)行模擬的時候,馬爾科夫鏈的穩(wěn)定性依賴于轉(zhuǎn)移核的選取。目前,常用于轉(zhuǎn)移核構(gòu)造的方法有兩種,分別是Metropolis-Hastings算法[7]和Gibbs抽樣[8]。Metropolis-Hastings算法(簡稱M-H算法)是最早提出是在1953年,經(jīng)過后人改進(jìn)成為目前應(yīng)用較為廣泛的MCMC方法。

    Metropolis-Hastings算法基本步驟如下。

    ①確定參數(shù)的初值θ0以及轉(zhuǎn)移核q(·|θi)。

    ②進(jìn)行迭代,從q(·|θi)抽取一個推薦值θ*。

    ③計算相應(yīng)的接受概率αi。

    (6)

    ④以概率αi接受θ*為下一個θi+1,即:

    (7)

    式中:μ是均勻分布在一定區(qū)間內(nèi)的隨機(jī)數(shù),其取值范圍在0~1之間。

    (5)重復(fù)①~④的步驟n次,去除前面k個不穩(wěn)定序列后,認(rèn)為θk+1,θk+2,…為滿足抽樣要求的平穩(wěn)序列,根據(jù)這些平穩(wěn)序列進(jìn)行后驗分布的各種統(tǒng)計推斷(后驗均值及各個分位數(shù)等)。從極值理論角度來講,通過式(7)給出的概率進(jìn)行優(yōu)選可以使得序列達(dá)到平穩(wěn)狀態(tài),并且保證模擬序列具備所求的邊緣分布的統(tǒng)計特征。

    3 GEV模型分布參數(shù)與洪水的估計

    作為搭建極值變量的一個重要分布,GEV分布在水文、氣象、環(huán)境、保險和金融等領(lǐng)域得到廣泛的應(yīng)用[9, 10]。對于GEV分布,可以采用極大似然估計法和矩估計等對其參數(shù)進(jìn)行估計,其中極大似然估計方法是目前采用較多的方法[11]。盡管極大似然估計法可以計算極值變量分位數(shù)的置信區(qū)間,以此表達(dá)相應(yīng)估計值的不確定性,然而,極大似然估計方法要求待估計的樣本具有漸近正態(tài)性的特征,對于樣本長度要求較高,當(dāng)樣本系列較小時,極大似然估計方法不能表達(dá)分位數(shù)估計的中所帶來不確定性。Martins和Stedinger[12]指出,序列長度對極大似然估計有一定的影響,較小的樣本會造成極大似然估計方法估計出的形狀參數(shù)偏小,從而影響對于總體分位數(shù)估計的準(zhǔn)確性。

    GEV分布函數(shù)可以表示為:

    (8)

    式中:μ、σ、ξ分別表示GEV分布模型的位置參數(shù)(Location parameter)、尺度參數(shù)(Scale parameter)和形狀參數(shù)(Shape parameter),且同時滿足條件為μ∈R,σ>0,ξ∈R,1+ξ(x-μ)σ>0。

    分別選取6個站點的年最大洪峰流量系列建立GEV分布模型。為方便計算和表達(dá),洪峰流量系列的單位統(tǒng)一為103m3/s,以下類同。令φ=logσ,以保證σ>0。先驗密度參數(shù)確定如下:選擇先驗密度函數(shù)為π(μ,φ,ξ)=πμ(μ)πφ(φ)πξ(ξ)。其中πμ(·)、πφ(·)和πξ(·)為正態(tài)密度函數(shù),其統(tǒng)計特征是均值為零,方差分別為υμ、υφ和υξ,同時保證參數(shù)μ、φ和ξ是相互獨立的。根據(jù)文獻(xiàn)[11]的推薦,初步設(shè)定υμ=υφ=104,υξ=100,從而幫助密度函數(shù)變得更加平滑。

    參考上文中M-H算法的步驟[13],利用概率αi對向量(μ,σ,ξ)的每個分量進(jìn)行篩選,同時定義每個分量在坐標(biāo)軸上隨機(jī)移動的步長為:μ*=μ+εμ,φ*=φ+εφ,ξ*=ξ+εξ,其中,εμ、εφ、εξ是均值為0,方差分別為ωμ、ωφ、ωξ的正態(tài)變量。ωμ、ωφ、ωξ的數(shù)值根據(jù)試驗確定。

    4 大渡河流域典型站點分析

    選取大渡河流域6個典型站點年最大洪峰系列進(jìn)行分析,分別為丹巴、大金、猴子巖、瀘定、瀑布溝和雙江口,6個站點覆蓋流域大部分地區(qū),具有一定代表性,將各個站點年最大洪峰流量系列作為洪水頻率計算的樣本數(shù)據(jù)。

    4.1 參數(shù)估計

    設(shè)定每個參數(shù)的訓(xùn)練次數(shù)為20 000次,通過MCMC模擬產(chǎn)生丹巴站GEV分布的3個參數(shù)序列如圖2所示。其中,位置參數(shù)和形狀參數(shù)為原序列,而尺度參數(shù)是經(jīng)過σi=eφi轉(zhuǎn)化以后得到的模擬序列。由于初始值的邊緣分布并不一定滿足后驗分布要求,因此去掉前面不平穩(wěn)的序列之后,可以認(rèn)為將其余比較穩(wěn)定的模擬值當(dāng)做符合要求的后驗分布的觀測值。

    由于方差ωμ、ωφ、ωξ的設(shè)定具有一定的主觀性,為了進(jìn)一步明晰不同的ωμ、ωφ、ωξ對參數(shù)收斂速度的影響,設(shè)置了2組不同的數(shù)值進(jìn)行對比分析如圖2所示,在圖2(a)的MCMC模擬中,設(shè)置ωμ、ωφ、ωξ的取值分別為0.02、0.01、0.1,從圖2(a)中可以看出,經(jīng)過2500次迭代后序列趨于穩(wěn)定;而圖2(b)的MCMC模擬中,設(shè)置ωμ、ωφ、ωξ的取值分別為0.1、0.05、0.5,為圖2(a)的5倍,從圖2(b)中可以看出,經(jīng)過1000次迭代后序列趨于穩(wěn)定。通過對比可以看出,設(shè)置較大的方差,可以一定程度上加快MCMC序列的收斂速度,同時,對比兩個情景下參數(shù)的收斂值可以發(fā)現(xiàn),兩個情景下的參數(shù)收斂值基本上是一致的,這說明選取ωμ、ωφ、ωξ的大小僅僅對于模型收斂速度產(chǎn)生影響,而并不影響對模型參數(shù)的估計。

    按照上述的方法,選取ωμ=0.1、ωφ=0.05、ωξ=0.5的值得出其余5個站點的MCMC序列如圖3所示。

    根據(jù)20 000組參數(shù),得到每個站點GEV各個參數(shù)后驗分布統(tǒng)計特征如表1所示。相比傳統(tǒng)的參數(shù)估計方法,貝葉斯分布不僅給出參數(shù)的估計值,同時也給出參數(shù)置信區(qū)間,通過置信區(qū)間表示參數(shù)估計的不確定性。

    4.2 擬合優(yōu)度檢驗

    (1)圖像分析法。圖像分析法一種較為直觀的利用圖形來描述擬合的優(yōu)劣程度的方法。主要是通過點繪理論聯(lián)合概率值和經(jīng)驗聯(lián)合概率值,如果得到的點距較均勻地分布在45°線附近,則說明建立的概率分布模型是合理的[14]。圖4給出了大渡河流域各個站點最大洪峰流量系列GEV模型經(jīng)驗頻率與理論頻率擬合效果圖,從圖4可以看出,樣本系列基本都落在了經(jīng)驗頻率與理論頻率所在的45°線上,說明通過MCMC序列估計出來的GEV分布模型與實際樣本系列擬合程度較好。

    圖2 丹巴站Bayes參數(shù)估計的MCMC模擬圖Fig.2 Bayes MCMC simulation for parameters at Danba station

    圖3 各站Bayes參數(shù)估計的MCMC模擬圖 Fig.3 Bayes MCMC simulation for parameters of GEV model

    圖4 大渡河流域各個站點GEV模型經(jīng)驗頻率與理論頻率擬合效果圖Fig.4 The fitting effect chart for empirical frequency and theoretical frequency of GEV model at different stations in Dadu River basin

    站點參數(shù)置信下限期望值置信上限位置參數(shù)3.2253.3963.568丹巴尺度參數(shù)0.3910.4980.644形狀參數(shù)-0.0900.0570.282位置參數(shù)1.8482.0642.277大金尺度參數(shù)0.4700.6120.802形狀參數(shù)-0.262-0.0560.203位置參數(shù)2.8603.0343.218猴子巖尺度參數(shù)0.4350.5430.682形狀參數(shù)-0.1260.0120.212位置參數(shù)3.1073.3273.556瀘定尺度參數(shù)0.5440.6870.868形狀參數(shù)-0.157-0.0240.154位置參數(shù)4.3744.5634.787瀑布溝尺度參數(shù)0.4890.6310.829形狀參數(shù)-0.0210.1700.400位置參數(shù)2.0302.2202.419雙江口尺度參數(shù)0.4380.5670.727形狀參數(shù)-0.0910.0690.284

    (2)K-S檢驗及OSL、AIC準(zhǔn)則。

    離差平方和最小準(zhǔn)則(OSL):

    (9)

    OSL值越小,說明GEV模型擬合得越好。

    AIC信息準(zhǔn)則:

    (10)

    AIC=nln(MSE)+2k

    (11)

    式中:Fi表示經(jīng)驗頻率;Ci表示理論頻率;n為分布樣本的個數(shù);k為模型參數(shù)的個數(shù),在GEV模型中k=3。

    與OSL相同,AIC值越小,說明GEV模型擬合得越好。

    采用OSL和AIC、Kolmogrov-Smirnow(簡稱K-S)法[15]3種指標(biāo),對傳統(tǒng)的兩種參數(shù)估計方法和貝葉斯MCMC方法得到模型擬合效果進(jìn)行了對比分析,評定結(jié)果如表2所示。3種參數(shù)估計方法均通過K-S方法的檢驗(顯著水平5%),說明3種方法均適用于估計GEV模型的參數(shù)。對比各個參數(shù)估計方法的OSL和AIC可以看出,對于丹巴、大金、猴子巖、瀑布溝和瀘定站點,由貝葉斯估計參數(shù)得出的OSL和AIC的值是最小的,也就是說線型擬合的效果最好,而對于雙江口來說,矩估計得出的OSL和AIC值最小,但是矩估計的OSL和AIC值與貝葉斯估計相差較小,可以認(rèn)為矩估計和貝葉斯估計擬合出效果近似相同,極大似然估計相對比其他兩種估計方法,擬合效果相對較差。與傳統(tǒng)的參數(shù)估計方法相比,貝葉斯估計方法計算的OSL和AIC值總體上更小,這說明GEV模型擬合的效果更好,貝葉斯估計方法具有一定的優(yōu)越性。

    表2 各個站點擬合優(yōu)度檢驗結(jié)果Tab.2 Results of goodness of fit test

    4.3 洪水頻率不確定性分析

    估計一定重現(xiàn)期下的洪水設(shè)計值是水文頻率分析的主要目的之一,廣義極值分布的p(0

    xp=μ-σ[1-(-logp)-ξ]/ξ

    (12)

    當(dāng)(μi,σi,ξi)時,即Gumbel分布的分位數(shù)為:

    xp=μ-σlog(-logp)

    (13)

    (14)

    (15)

    貝葉斯方法將向量(μi,σi,ξi)的每個模擬值代入式(16)就得到相應(yīng)的1/(1-p)年重現(xiàn)期的設(shè)計洪峰流量值的后驗分布樣本及置信區(qū)間。

    (16)

    以丹巴站為例,根據(jù)式(16)得出丹巴站年最大洪峰流量的重現(xiàn)水平圖如圖5所示,可以看出隨著重現(xiàn)期的增加,設(shè)計洪峰流量的置信區(qū)間也在增大,洪水頻率的不確定性也在增加。圖6給出了丹巴站各個典型重現(xiàn)期下的設(shè)計洪峰流量的后驗密度估計圖。根據(jù)后驗密度估計,計算出流域各站點各典型重現(xiàn)期下的設(shè)計洪峰流量值及95%置信度下的置信區(qū)間,并與傳統(tǒng)Delta進(jìn)行比較,結(jié)果如表3所示。由表3可以看出,根據(jù)貝葉斯MCMC方法得到的各個站點各個重現(xiàn)期的設(shè)計值總是小于其置信下限和置信上限的平均值,而傳統(tǒng)Delta方法得出的各個站點各個重現(xiàn)期的設(shè)計值是等于其置信下限和置信上限的平均值的,實際中,由于大的洪水資料比較有限,在估計某個洪水頻率設(shè)計值上限時往往不確定性要大于下限的不確定性,由此可以看出根據(jù)貝葉斯MCMC方法得到的置信區(qū)間往往更接近實際的情況。同時,從表3可以看到,重現(xiàn)期越長,重現(xiàn)水平越大,置信區(qū)間的寬度越寬,代表不確定性越大,與Delta方法相比,貝葉斯MCMC方法得到的置信區(qū)間的寬度相對要小一些,說明該方法在評估洪水頻率不確定性時要比Delta方法得到的不確定性要小。綜上所述,貝葉斯MCMC方法在評估洪水頻率較傳統(tǒng)的Delta方法具有一定的優(yōu)越性。

    圖5 丹巴站Bayes重現(xiàn)水平圖Fig.5 Return level diagram of GEV model at Danba station

    5 結(jié) 語

    以GEV模型為洪水頻率分析線型,應(yīng)用貝葉斯MCMC算法對大渡河流域洪水頻率分析不確定性進(jìn)行研究,研究結(jié)果表明,Metropolis-Hastings抽樣的貝葉斯MCMC方法是洪水頻率參數(shù)估計的一個有效方法,與傳統(tǒng)的極大似然估計和矩估計相比,貝葉斯MCMC估計方法的擬合效果略優(yōu)于另外兩種方法。貝葉斯MCMC方法通過利用參數(shù)的先驗信息得到參數(shù)以及設(shè)計值的后驗密度函數(shù),通過參數(shù)和設(shè)計值概率分布的形式表達(dá)洪水頻率中的不確定性信息。同時,將貝葉斯MCMC方法得出的置信區(qū)間與傳統(tǒng)Delta方法得出的置信區(qū)間進(jìn)行比較,發(fā)現(xiàn)貝葉斯MCMC方法得到的置信區(qū)間的寬度相對要小 一些,置信上限和與估計值之間距離大于置信下限與估計值之間距離,這種不對等性與實際更加接近,更能精準(zhǔn)地估計洪水頻率的置信區(qū)間。

    圖6 丹巴站各個典型重現(xiàn)期設(shè)計洪峰流量后驗概率密度圖Fig.6 Posterior probability density map of flood peak for different return periods at Danba station

    m3/s

    對于洪水頻率不確定性研究,本文僅僅采用GEV一種線型進(jìn)行研究,沒有考慮線型選擇對于洪水頻率分析不確定性造成的影響,未來工作需要進(jìn)一步將國內(nèi)較為常見的P-Ⅲ、Lognormal等線型加入,綜合分析洪水頻率分析中各個環(huán)節(jié)的不確定性。

    [1] 朱慧明,韓玉啟. 貝葉斯多元統(tǒng)計推斷理論[M]. 上海:科學(xué)出版社, 2006.

    [2] Wood E F, Rodríguez Iturbe I. Bayesian inference and decision making for extreme hydrologic events[J]. Water Resources Research, 1975,11(4):533-542.

    [3] Claps P, Laio F. Can continuous streamflow data support flood frequency analysis? An alternative to the partial duration series approach[J]. Water Resources Research, 2003,39(8):1 216-1 225.

    [4] Fernandes W, Naghettini M, Loschi R. A Bayesian approach for estimating extreme flood probabilities with upper-bounded distribution functions[J]. Stochastic Environmental Research and Risk Assessment, 2010,24(8):1 127-1 143.

    [5] 朱 嵩,毛根海,劉國華,等. 改進(jìn)的 MCMC 方法及其應(yīng)用[J]. 水利學(xué)報,2009,(8):1 019-1 023.

    [6] 梁忠民,戴 榮,李彬權(quán). 基于貝葉斯理論的水文不確定性分析研究進(jìn)展[J]. 水科學(xué)進(jìn)展,2010,21(2):274-281.

    [7] 陳 平,徐若曦. Metropolis-Hastings 自適應(yīng)算法及其應(yīng)用[J]. 系統(tǒng)工程理論與實踐,2008,28(1):100-108.

    [8] 張志華,屈 斐. 軟件可靠性模型的 Bayes 推斷及 Gibbs 算法[J]. 高校應(yīng)用數(shù)學(xué)學(xué)報(A 輯),2005,20(4):401-408.

    [9] Coles S G, Tawn J A. A Bayesian analysis of extreme rainfall data[J]. Applied Statistics, 1996:463-478.

    [10] Walshaw D. Modelling extreme wind speeds in regions prone to hurricanes[J]. Journal of the Royal Statistical Society (Series C Applied Statistics), 2000,49(1):51-62.

    [11] Huang W, Xu S, Nnaji S. Evaluation of GEV model for frequency analysis of annual maximum water levels in the coast of United States[J]. Ocean Engineering, 2008,35(11):1 132-1 147.

    [12] Reis Jr D S, Stedinger J R. Bayesian MCMC flood frequency analysis with historical information[J]. Journal of Hydrology, 2005,313(1):97-116.

    [13] 魯 帆,嚴(yán)登華. 基于廣義極值分布和 Metropolis-Hastings 抽樣算法的貝葉斯 MCMC 洪水頻率分析方法[J]. 水利學(xué)報,2013,44(8):942-949.

    [14] 曾惠芳,朱慧明,李素芳,等. 基于 MH 算法的貝葉斯分位自回歸模型[J]. 湖南大學(xué)學(xué)報(自然科學(xué)版),2010,(2):88-92.

    [15] 張冬冬,魯 帆,嚴(yán)登華,等. 基于 Archimedean Copula 函數(shù)的洪水多要素聯(lián)合概率分布研究[J]. 中國農(nóng)村水利水電,2015,(1):68-74.

    猜你喜歡
    后驗置信區(qū)間參數(shù)估計
    定數(shù)截尾場合三參數(shù)pareto分布參數(shù)的最優(yōu)置信區(qū)間
    基于新型DFrFT的LFM信號參數(shù)估計算法
    p-范分布中參數(shù)的置信區(qū)間
    多個偏正態(tài)總體共同位置參數(shù)的Bootstrap置信區(qū)間
    基于對偶理論的橢圓變分不等式的后驗誤差分析(英)
    貝葉斯統(tǒng)計中單參數(shù)后驗分布的精確計算方法
    列車定位中置信區(qū)間的確定方法
    一種基于最大后驗框架的聚類分析多基線干涉SAR高度重建算法
    Logistic回歸模型的幾乎無偏兩參數(shù)估計
    基于向前方程的平穩(wěn)分布參數(shù)估計
    天堂动漫精品| xxxhd国产人妻xxx| 啦啦啦在线免费观看视频4| 国产精品国产高清国产av| 久久久久久人人人人人| 久久精品国产清高在天天线| 一级,二级,三级黄色视频| 欧美黄色淫秽网站| 精品久久久精品久久久| 男人舔女人的私密视频| 亚洲视频免费观看视频| 久久久久久久精品吃奶| 亚洲精品国产一区二区精华液| 亚洲国产欧美网| 精品国产亚洲在线| 国产xxxxx性猛交| av在线天堂中文字幕 | 日韩欧美在线二视频| 美国免费a级毛片| 动漫黄色视频在线观看| 99精品欧美一区二区三区四区| 亚洲精华国产精华精| 欧美一级毛片孕妇| 欧美乱色亚洲激情| 又紧又爽又黄一区二区| 亚洲男人的天堂狠狠| 亚洲五月婷婷丁香| 很黄的视频免费| 成人亚洲精品av一区二区 | 夜夜看夜夜爽夜夜摸 | 大陆偷拍与自拍| 成人精品一区二区免费| 色综合欧美亚洲国产小说| 国产欧美日韩综合在线一区二区| 美女 人体艺术 gogo| 少妇的丰满在线观看| 精品日产1卡2卡| 精品国内亚洲2022精品成人| 久久青草综合色| 国产精品电影一区二区三区| 在线播放国产精品三级| xxxhd国产人妻xxx| 亚洲 欧美 日韩 在线 免费| 亚洲成人国产一区在线观看| 桃红色精品国产亚洲av| 久久热在线av| 精品一区二区三区视频在线观看免费 | 成人三级黄色视频| 一夜夜www| 夜夜夜夜夜久久久久| 中文字幕高清在线视频| 好男人电影高清在线观看| 超色免费av| 国产精品亚洲av一区麻豆| 久久国产亚洲av麻豆专区| 妹子高潮喷水视频| 国产一区二区在线av高清观看| 一边摸一边抽搐一进一小说| 国产单亲对白刺激| 亚洲自偷自拍图片 自拍| 99国产极品粉嫩在线观看| 美女高潮到喷水免费观看| 成人黄色视频免费在线看| 亚洲情色 制服丝袜| 亚洲欧美日韩另类电影网站| 国产欧美日韩一区二区三| 中亚洲国语对白在线视频| 精品久久久久久成人av| 国产亚洲av高清不卡| 这个男人来自地球电影免费观看| 成人18禁高潮啪啪吃奶动态图| 激情在线观看视频在线高清| 亚洲视频免费观看视频| 在线观看免费日韩欧美大片| 亚洲精品av麻豆狂野| 女人被躁到高潮嗷嗷叫费观| 精品久久久久久久毛片微露脸| 午夜老司机福利片| netflix在线观看网站| 久久久久亚洲av毛片大全| 成人亚洲精品一区在线观看| 无人区码免费观看不卡| 亚洲自偷自拍图片 自拍| 成人三级黄色视频| 免费看十八禁软件| 国产激情欧美一区二区| 欧美日本亚洲视频在线播放| 久久中文字幕一级| 久久伊人香网站| 九色亚洲精品在线播放| 亚洲精品一二三| 久久精品亚洲av国产电影网| 男男h啪啪无遮挡| 成人国产一区最新在线观看| 国产精品自产拍在线观看55亚洲| 性少妇av在线| 色老头精品视频在线观看| 精品国产一区二区三区四区第35| 国产黄色免费在线视频| 国产99久久九九免费精品| 欧美黑人精品巨大| 真人一进一出gif抽搐免费| 精品一区二区三区av网在线观看| 国产99白浆流出| www.999成人在线观看| 亚洲一区二区三区欧美精品| 精品熟女少妇八av免费久了| 欧美中文综合在线视频| 久久精品国产综合久久久| 久久香蕉精品热| 999久久久精品免费观看国产| 亚洲免费av在线视频| 日日夜夜操网爽| 亚洲av片天天在线观看| 看黄色毛片网站| 757午夜福利合集在线观看| 涩涩av久久男人的天堂| 国产一区二区三区综合在线观看| 中文字幕另类日韩欧美亚洲嫩草| 色婷婷av一区二区三区视频| 天堂中文最新版在线下载| 午夜免费激情av| 一进一出抽搐动态| 在线播放国产精品三级| 精品久久久精品久久久| 国产99白浆流出| 在线观看一区二区三区| 天天添夜夜摸| 午夜免费激情av| 久久久国产成人精品二区 | 男女下面进入的视频免费午夜 | 99香蕉大伊视频| 在线观看日韩欧美| 成人三级黄色视频| 欧美黑人精品巨大| 亚洲精品美女久久久久99蜜臀| 亚洲欧美一区二区三区黑人| 免费观看人在逋| 18禁裸乳无遮挡免费网站照片 | 国产精品久久久久成人av| 亚洲精品国产区一区二| av视频免费观看在线观看| 男女之事视频高清在线观看| 亚洲人成电影免费在线| 亚洲欧美日韩高清在线视频| 国产精品久久久av美女十八| 真人一进一出gif抽搐免费| 热re99久久国产66热| 天天躁狠狠躁夜夜躁狠狠躁| 久久九九热精品免费| 夜夜看夜夜爽夜夜摸 | 国产片内射在线| 激情视频va一区二区三区| 国产三级黄色录像| 日韩欧美在线二视频| 国产一区二区三区视频了| 亚洲av电影在线进入| 精品日产1卡2卡| 成人特级黄色片久久久久久久| 国产一区二区在线av高清观看| 十分钟在线观看高清视频www| 老熟妇乱子伦视频在线观看| 亚洲熟妇中文字幕五十中出 | 国产高清激情床上av| 欧美乱妇无乱码| 亚洲av成人一区二区三| 国产精品乱码一区二三区的特点 | 免费高清视频大片| 男女高潮啪啪啪动态图| 国产精品综合久久久久久久免费 | www.www免费av| 国产精品电影一区二区三区| 国产一区在线观看成人免费| 久久精品91无色码中文字幕| 精品国产国语对白av| av欧美777| 青草久久国产| 国产欧美日韩一区二区三区在线| 国内毛片毛片毛片毛片毛片| 国产成人啪精品午夜网站| 亚洲精品成人av观看孕妇| 久久国产精品男人的天堂亚洲| 久久亚洲真实| 精品国产国语对白av| 日韩欧美国产一区二区入口| 男插女下体视频免费在线播放| 成人精品一区二区免费| 午夜福利在线在线| 最后的刺客免费高清国语| 国产综合懂色| 午夜免费成人在线视频| 亚洲欧美日韩卡通动漫| 色综合婷婷激情| 2021天堂中文幕一二区在线观| 中亚洲国语对白在线视频| 观看免费一级毛片| 精品久久久久久久久亚洲 | 国内精品久久久久久久电影| 国产午夜精品论理片| 黄色女人牲交| 欧美3d第一页| 欧美日韩综合久久久久久 | 在线a可以看的网站| 国产欧美日韩精品亚洲av| 精品人妻熟女av久视频| 又紧又爽又黄一区二区| 哪里可以看免费的av片| 精品99又大又爽又粗少妇毛片 | 精品久久久久久成人av| 欧美在线一区亚洲| 嫩草影院入口| 老司机午夜十八禁免费视频| 亚洲精品乱码久久久v下载方式| 99久久精品国产亚洲精品| 欧美国产日韩亚洲一区| 一本综合久久免费| 成人毛片a级毛片在线播放| 色哟哟哟哟哟哟| av在线蜜桃| 国产私拍福利视频在线观看| 亚洲精品影视一区二区三区av| 成年女人毛片免费观看观看9| 色5月婷婷丁香| 国产白丝娇喘喷水9色精品| 美女被艹到高潮喷水动态| 成人国产综合亚洲| 悠悠久久av| 老熟妇仑乱视频hdxx| 男人舔女人下体高潮全视频| 一区二区三区四区激情视频 | 别揉我奶头 嗯啊视频| 日日摸夜夜添夜夜添小说| 国产一区二区三区视频了| 草草在线视频免费看| 琪琪午夜伦伦电影理论片6080| 97热精品久久久久久| 国产精品久久久久久久久免 | 特级一级黄色大片| 国产精品久久久久久亚洲av鲁大| 免费电影在线观看免费观看| 亚洲aⅴ乱码一区二区在线播放| 一级av片app| a级一级毛片免费在线观看| 国产白丝娇喘喷水9色精品| 九九久久精品国产亚洲av麻豆| 午夜亚洲福利在线播放| 成人美女网站在线观看视频| av在线观看视频网站免费| 午夜免费激情av| 国内精品久久久久精免费| 男女下面进入的视频免费午夜| 18禁黄网站禁片午夜丰满| 99国产精品一区二区三区| 99在线人妻在线中文字幕| 亚洲成a人片在线一区二区| 999久久久精品免费观看国产| a级毛片a级免费在线| 直男gayav资源| 日韩欧美 国产精品| 久久天躁狠狠躁夜夜2o2o| 欧美最黄视频在线播放免费| 亚州av有码| 老司机福利观看| 国产精品久久久久久久电影| 99国产综合亚洲精品| 免费黄网站久久成人精品 | 亚洲最大成人手机在线| 欧美+日韩+精品| 久久这里只有精品中国| 亚洲av不卡在线观看| 午夜两性在线视频| 国产精品久久视频播放| 国产精品电影一区二区三区| 无遮挡黄片免费观看| 国产精品爽爽va在线观看网站| 国产激情偷乱视频一区二区| 深夜a级毛片| 国产淫片久久久久久久久 | 我的女老师完整版在线观看| 俄罗斯特黄特色一大片| 久久精品国产亚洲av天美| 亚洲av五月六月丁香网| 久久这里只有精品中国| 一夜夜www| 特大巨黑吊av在线直播| 老司机深夜福利视频在线观看| 久久精品国产自在天天线| 在线观看美女被高潮喷水网站 | 欧美xxxx性猛交bbbb| 欧美极品一区二区三区四区| 波多野结衣巨乳人妻| 国语自产精品视频在线第100页| 我的女老师完整版在线观看| 亚洲精品影视一区二区三区av| 色综合欧美亚洲国产小说| 国语自产精品视频在线第100页| 久久亚洲精品不卡| 99久国产av精品| 在线观看免费视频日本深夜| 69av精品久久久久久| 亚洲成av人片在线播放无| 欧美午夜高清在线| 国产精品永久免费网站| 日韩人妻高清精品专区| 1024手机看黄色片| 日韩精品中文字幕看吧| 我要看日韩黄色一级片| 十八禁国产超污无遮挡网站| 免费电影在线观看免费观看| 亚洲天堂国产精品一区在线| 搡老岳熟女国产| 午夜精品久久久久久毛片777| 99热精品在线国产| 午夜影院日韩av| 国产亚洲欧美在线一区二区| 亚洲人与动物交配视频| 亚洲av中文字字幕乱码综合| 亚洲乱码一区二区免费版| 99在线人妻在线中文字幕| 国产美女午夜福利| 亚洲成av人片在线播放无| 能在线免费观看的黄片| 嫩草影院入口| 亚洲av不卡在线观看| 亚洲成人免费电影在线观看| 我的女老师完整版在线观看| 国产真实乱freesex| 脱女人内裤的视频| 国产精品久久久久久亚洲av鲁大| 久久久久久国产a免费观看| 亚洲成a人片在线一区二区| 久久香蕉精品热| 欧美zozozo另类| av专区在线播放| 久久久久久久久久黄片| 内地一区二区视频在线| 婷婷六月久久综合丁香| 99热这里只有精品一区| 亚洲真实伦在线观看| 国产成人av教育| 精品不卡国产一区二区三区| 性色av乱码一区二区三区2| 老司机福利观看| 国产av一区在线观看免费| 国产中年淑女户外野战色| 亚洲一区高清亚洲精品| 特级一级黄色大片| 自拍偷自拍亚洲精品老妇| av视频在线观看入口| 亚洲久久久久久中文字幕| 搞女人的毛片| 国产主播在线观看一区二区| 97热精品久久久久久| 亚洲三级黄色毛片| 成人毛片a级毛片在线播放| 一个人免费在线观看电影| 精品99又大又爽又粗少妇毛片 | 久久久久久久久中文| 国产免费av片在线观看野外av| 亚洲美女搞黄在线观看 | 两性午夜刺激爽爽歪歪视频在线观看| 一夜夜www| 脱女人内裤的视频| 亚洲av日韩精品久久久久久密| 99久久精品一区二区三区| 美女黄网站色视频| 久久久久久久亚洲中文字幕 | 直男gayav资源| 亚洲片人在线观看| 亚洲成a人片在线一区二区| 国产男靠女视频免费网站| 午夜精品在线福利| 天天一区二区日本电影三级| 九九久久精品国产亚洲av麻豆| 身体一侧抽搐| 五月玫瑰六月丁香| 欧美xxxx性猛交bbbb| 中文资源天堂在线| 村上凉子中文字幕在线| а√天堂www在线а√下载| 精品午夜福利视频在线观看一区| 国产亚洲精品久久久久久毛片| www.999成人在线观看| 老司机午夜福利在线观看视频| 国产黄片美女视频| 欧美最新免费一区二区三区 | 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 在线国产一区二区在线| 黄色女人牲交| 麻豆成人午夜福利视频| 欧美日本视频| 国产精品免费一区二区三区在线| 国产视频内射| 中文亚洲av片在线观看爽| 91麻豆精品激情在线观看国产| 18禁在线播放成人免费| 观看免费一级毛片| 日韩免费av在线播放| 国产黄色小视频在线观看| 内地一区二区视频在线| 小说图片视频综合网站| 午夜影院日韩av| 极品教师在线免费播放| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲美女久久久| 欧美日本视频| 亚洲国产精品sss在线观看| 日日干狠狠操夜夜爽| 又黄又爽又刺激的免费视频.| 国内精品美女久久久久久| 99久久无色码亚洲精品果冻| 亚洲欧美日韩高清在线视频| 搡老岳熟女国产| 一边摸一边抽搐一进一小说| 91狼人影院| 久久精品久久久久久噜噜老黄 | 老女人水多毛片| 一级黄色大片毛片| or卡值多少钱| 99热这里只有精品一区| 成人av一区二区三区在线看| www.999成人在线观看| 少妇丰满av| 在线观看66精品国产| 婷婷色综合大香蕉| 女人十人毛片免费观看3o分钟| 国产精品电影一区二区三区| 国产野战对白在线观看| 久久久国产成人免费| 一个人看视频在线观看www免费| 99热只有精品国产| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 欧美午夜高清在线| 国产视频内射| 九九在线视频观看精品| 老熟妇仑乱视频hdxx| 亚洲 欧美 日韩 在线 免费| 69av精品久久久久久| 欧美日本视频| 亚洲自偷自拍三级| 在线播放国产精品三级| 色噜噜av男人的天堂激情| 成人鲁丝片一二三区免费| 欧美日韩黄片免| 亚洲精品影视一区二区三区av| 久久精品国产亚洲av涩爱 | 99国产综合亚洲精品| 国产成人a区在线观看| 少妇人妻一区二区三区视频| 国产淫片久久久久久久久 | 国产私拍福利视频在线观看| 一区二区三区免费毛片| 老司机深夜福利视频在线观看| 成人特级av手机在线观看| 美女被艹到高潮喷水动态| 美女黄网站色视频| 色综合欧美亚洲国产小说| 悠悠久久av| 69av精品久久久久久| 好看av亚洲va欧美ⅴa在| 亚洲精品粉嫩美女一区| 国产单亲对白刺激| 一级黄色大片毛片| 欧美日韩中文字幕国产精品一区二区三区| 熟女人妻精品中文字幕| 日本成人三级电影网站| 最近视频中文字幕2019在线8| 天堂网av新在线| av中文乱码字幕在线| 简卡轻食公司| 一级黄色大片毛片| 国产高清有码在线观看视频| 日日夜夜操网爽| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 一区二区三区高清视频在线| 成人午夜高清在线视频| 国产成人a区在线观看| 国产欧美日韩一区二区三| 国产成人aa在线观看| 国产一区二区三区视频了| 日韩欧美精品免费久久 | 日本精品一区二区三区蜜桃| 国内精品久久久久精免费| 乱人视频在线观看| 国产黄片美女视频| 高清毛片免费观看视频网站| 91午夜精品亚洲一区二区三区 | 级片在线观看| 日韩欧美三级三区| 在线播放无遮挡| 国内毛片毛片毛片毛片毛片| 九九在线视频观看精品| 给我免费播放毛片高清在线观看| 搞女人的毛片| 亚洲性夜色夜夜综合| 国产成人a区在线观看| 欧美在线黄色| 国产成人av教育| 人妻制服诱惑在线中文字幕| 黄色日韩在线| www日本黄色视频网| 国产日本99.免费观看| 一区二区三区激情视频| 黄色女人牲交| 亚洲七黄色美女视频| 国产亚洲精品综合一区在线观看| 亚洲精品成人久久久久久| 亚洲欧美清纯卡通| 成人国产综合亚洲| 久久精品综合一区二区三区| 久久人妻av系列| 亚洲精品日韩av片在线观看| 国产日本99.免费观看| 亚洲av二区三区四区| 免费大片18禁| 又黄又爽又刺激的免费视频.| 日韩欧美一区二区三区在线观看| 一级黄片播放器| 婷婷精品国产亚洲av在线| 国产真实乱freesex| 熟女电影av网| 国产亚洲精品久久久com| 伊人久久精品亚洲午夜| 国产精品女同一区二区软件 | 哪里可以看免费的av片| 欧美xxxx性猛交bbbb| 亚洲在线自拍视频| 国内少妇人妻偷人精品xxx网站| а√天堂www在线а√下载| 日韩国内少妇激情av| 亚洲七黄色美女视频| 99精品在免费线老司机午夜| 婷婷色综合大香蕉| x7x7x7水蜜桃| 亚洲av成人不卡在线观看播放网| 国产淫片久久久久久久久 | 禁无遮挡网站| 免费一级毛片在线播放高清视频| 床上黄色一级片| 国产精品伦人一区二区| 国产av在哪里看| 国产精品久久视频播放| 免费av不卡在线播放| 91在线精品国自产拍蜜月| 婷婷亚洲欧美| 欧美性感艳星| 国内精品久久久久精免费| 男女床上黄色一级片免费看| 99热这里只有是精品50| 久久精品国产自在天天线| 老司机午夜福利在线观看视频| 成年女人看的毛片在线观看| 男女之事视频高清在线观看| 最近最新中文字幕大全电影3| 永久网站在线| 国产高清有码在线观看视频| 久久草成人影院| 欧美日韩黄片免| 毛片一级片免费看久久久久 | 亚洲最大成人手机在线| 国产成+人综合+亚洲专区| 国产91精品成人一区二区三区| 精品无人区乱码1区二区| 亚洲,欧美,日韩| 亚州av有码| 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄 | 婷婷丁香在线五月| 国产亚洲精品久久久久久毛片| 亚洲av一区综合| 日本一本二区三区精品| 深夜a级毛片| a在线观看视频网站| 成人亚洲精品av一区二区| 国产乱人视频| 中亚洲国语对白在线视频| 久久久久久久久中文| 亚洲美女黄片视频| 亚洲av免费高清在线观看| 色哟哟哟哟哟哟| 欧美成人一区二区免费高清观看| av福利片在线观看| 听说在线观看完整版免费高清| 国产色爽女视频免费观看| 亚洲中文字幕一区二区三区有码在线看| 超碰av人人做人人爽久久| 日韩欧美在线乱码| 制服丝袜大香蕉在线| 精品久久久久久久久久免费视频| 91麻豆精品激情在线观看国产| 国产中年淑女户外野战色| 色哟哟哟哟哟哟| 欧美成人a在线观看| .国产精品久久| 亚洲国产精品成人综合色| 亚洲人成网站在线播放欧美日韩| 97热精品久久久久久| 我的老师免费观看完整版| 亚洲精品久久国产高清桃花| 亚洲经典国产精华液单 | 一卡2卡三卡四卡精品乱码亚洲| 久久久久久九九精品二区国产| 麻豆国产av国片精品| 伊人久久精品亚洲午夜| 欧美在线一区亚洲| 欧美不卡视频在线免费观看| 性色avwww在线观看| 夜夜看夜夜爽夜夜摸| 性色avwww在线观看| 久久热精品热| 99在线视频只有这里精品首页| 此物有八面人人有两片|