• <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ù)估計
    91字幕亚洲| 天天操日日干夜夜撸| 亚洲一区二区三区欧美精品| 日本一区二区免费在线视频| 2018国产大陆天天弄谢| 国产1区2区3区精品| 一边摸一边抽搐一进一出视频| 蜜桃国产av成人99| 久久久精品免费免费高清| 大香蕉久久网| 十八禁人妻一区二区| 操美女的视频在线观看| 久久热在线av| 夫妻午夜视频| 老司机影院成人| 12—13女人毛片做爰片一| 新久久久久国产一级毛片| 精品一区二区三卡| www.999成人在线观看| 午夜福利一区二区在线看| av有码第一页| 麻豆国产av国片精品| 老司机午夜福利在线观看视频 | 亚洲天堂av无毛| 精品久久蜜臀av无| 精品卡一卡二卡四卡免费| 秋霞在线观看毛片| 美女中出高潮动态图| 日本黄色日本黄色录像| 国产成人影院久久av| 99热网站在线观看| 亚洲精品av麻豆狂野| 九色亚洲精品在线播放| 国产亚洲av片在线观看秒播厂| 精品乱码久久久久久99久播| 色精品久久人妻99蜜桃| 大码成人一级视频| 丰满少妇做爰视频| 每晚都被弄得嗷嗷叫到高潮| 欧美老熟妇乱子伦牲交| av网站免费在线观看视频| 人人妻人人澡人人看| 久久国产精品人妻蜜桃| www日本在线高清视频| 久久久久久人人人人人| 人妻一区二区av| 天天操日日干夜夜撸| 亚洲伊人久久精品综合| 一区二区三区四区激情视频| 99热全是精品| 两个人免费观看高清视频| 国产一区二区三区在线臀色熟女 | 国精品久久久久久国模美| 精品福利观看| 日韩 欧美 亚洲 中文字幕| 成人免费观看视频高清| 美女主播在线视频| 青青草视频在线视频观看| 免费在线观看视频国产中文字幕亚洲 | 亚洲少妇的诱惑av| 亚洲精品美女久久久久99蜜臀| 大型av网站在线播放| 女警被强在线播放| 久久久久久久精品精品| 成年美女黄网站色视频大全免费| 两人在一起打扑克的视频| 黄色a级毛片大全视频| 下体分泌物呈黄色| 欧美在线黄色| 免费观看人在逋| 别揉我奶头~嗯~啊~动态视频 | 91成人精品电影| 国产欧美日韩精品亚洲av| 最近最新免费中文字幕在线| 久久久精品94久久精品| 下体分泌物呈黄色| 一个人免费看片子| 丝袜美腿诱惑在线| 99热国产这里只有精品6| 久久影院123| 欧美精品啪啪一区二区三区 | 国产伦人伦偷精品视频| 久久久精品国产亚洲av高清涩受| 亚洲国产精品一区三区| 秋霞在线观看毛片| 另类亚洲欧美激情| 久久久久久久久久久久大奶| 777米奇影视久久| 国产免费一区二区三区四区乱码| 国产在线一区二区三区精| 久久国产精品男人的天堂亚洲| 亚洲国产日韩一区二区| 在线观看免费日韩欧美大片| 丰满迷人的少妇在线观看| 亚洲一区中文字幕在线| 亚洲国产看品久久| 国产伦人伦偷精品视频| 婷婷丁香在线五月| 欧美激情久久久久久爽电影 | 精品高清国产在线一区| 久久久久久久国产电影| 大香蕉久久网| 国产国语露脸激情在线看| 91老司机精品| 美女视频免费永久观看网站| 99re6热这里在线精品视频| 视频在线观看一区二区三区| 亚洲精品国产一区二区精华液| 夜夜骑夜夜射夜夜干| 三上悠亚av全集在线观看| 国产日韩欧美在线精品| 亚洲av日韩在线播放| 亚洲免费av在线视频| 日日摸夜夜添夜夜添小说| 精品国产超薄肉色丝袜足j| 亚洲成av片中文字幕在线观看| 99香蕉大伊视频| 男男h啪啪无遮挡| 丝瓜视频免费看黄片| 少妇裸体淫交视频免费看高清 | 欧美 日韩 精品 国产| 色94色欧美一区二区| 十分钟在线观看高清视频www| 五月天丁香电影| 80岁老熟妇乱子伦牲交| 免费看十八禁软件| 久久久久久久久免费视频了| 成人av一区二区三区在线看 | 精品高清国产在线一区| 久久国产精品影院| 亚洲视频免费观看视频| 久久九九热精品免费| av一本久久久久| 9热在线视频观看99| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 国产国语露脸激情在线看| 亚洲 欧美一区二区三区| 激情视频va一区二区三区| 黄片小视频在线播放| 亚洲av国产av综合av卡| 9色porny在线观看| 国产成人欧美| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| 国产激情久久老熟女| 91九色精品人成在线观看| 亚洲av成人不卡在线观看播放网 | 一级,二级,三级黄色视频| 日本猛色少妇xxxxx猛交久久| 久久久久网色| 日本五十路高清| 新久久久久国产一级毛片| 国精品久久久久久国模美| 一级黄色大片毛片| 国产亚洲精品第一综合不卡| 黄片小视频在线播放| 中文字幕最新亚洲高清| 欧美乱码精品一区二区三区| 欧美日韩av久久| 免费黄频网站在线观看国产| 亚洲视频免费观看视频| 国产有黄有色有爽视频| 亚洲天堂av无毛| 亚洲精品国产色婷婷电影| 精品国产国语对白av| 老司机在亚洲福利影院| 免费在线观看影片大全网站| 国产精品久久久av美女十八| 日本av免费视频播放| 亚洲免费av在线视频| 国产亚洲欧美精品永久| 一级毛片女人18水好多| 日韩欧美国产一区二区入口| 久久精品国产亚洲av高清一级| 中文精品一卡2卡3卡4更新| 在线 av 中文字幕| 国产xxxxx性猛交| a级毛片在线看网站| 久久精品久久久久久噜噜老黄| 欧美国产精品一级二级三级| 亚洲国产欧美日韩在线播放| 天天添夜夜摸| 久久久国产一区二区| 欧美 亚洲 国产 日韩一| 久久国产精品男人的天堂亚洲| 高清黄色对白视频在线免费看| 国产激情久久老熟女| 久热爱精品视频在线9| 国产精品国产av在线观看| 久久久久国内视频| 久久久久久久久免费视频了| 三上悠亚av全集在线观看| 国产成人精品在线电影| 女人精品久久久久毛片| 国产精品一区二区精品视频观看| 亚洲精品一二三| 亚洲五月婷婷丁香| 黄色怎么调成土黄色| 久久人人97超碰香蕉20202| 久久国产精品大桥未久av| 日韩大片免费观看网站| 一级黄色大片毛片| 性高湖久久久久久久久免费观看| 爱豆传媒免费全集在线观看| 下体分泌物呈黄色| 91精品三级在线观看| 一进一出抽搐动态| 五月天丁香电影| 国产精品久久久人人做人人爽| 汤姆久久久久久久影院中文字幕| 日日摸夜夜添夜夜添小说| 新久久久久国产一级毛片| 黄频高清免费视频| 99热国产这里只有精品6| 午夜免费鲁丝| 在线观看人妻少妇| 中国美女看黄片| 免费久久久久久久精品成人欧美视频| 狂野欧美激情性bbbbbb| 亚洲精品国产精品久久久不卡| 免费在线观看黄色视频的| 精品卡一卡二卡四卡免费| 搡老岳熟女国产| 老鸭窝网址在线观看| 国产人伦9x9x在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲av日韩在线播放| 欧美国产精品一级二级三级| 亚洲精品国产区一区二| 免费观看人在逋| 老司机福利观看| 亚洲欧美一区二区三区久久| 波多野结衣av一区二区av| 亚洲精华国产精华精| 国产日韩一区二区三区精品不卡| 亚洲精品中文字幕在线视频| 黄色视频,在线免费观看| 中文字幕色久视频| 久久久精品区二区三区| 99精国产麻豆久久婷婷| 国产日韩一区二区三区精品不卡| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 男女国产视频网站| 99国产极品粉嫩在线观看| 亚洲欧美一区二区三区黑人| 天堂8中文在线网| www.自偷自拍.com| 男女下面插进去视频免费观看| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 亚洲精品中文字幕一二三四区 | 免费看十八禁软件| 久久国产亚洲av麻豆专区| 久久久国产欧美日韩av| 欧美亚洲 丝袜 人妻 在线| 啦啦啦在线免费观看视频4| 国产野战对白在线观看| 大香蕉久久成人网| 国产精品久久久久成人av| 国产欧美日韩一区二区三 | 99久久国产精品久久久| 国产精品一区二区精品视频观看| 两个人免费观看高清视频| a在线观看视频网站| 在线观看免费高清a一片| 少妇精品久久久久久久| 99精品久久久久人妻精品| 一本综合久久免费| 男女床上黄色一级片免费看| 最黄视频免费看| 中文字幕精品免费在线观看视频| kizo精华| 日韩制服骚丝袜av| 黄色视频在线播放观看不卡| 欧美日本中文国产一区发布| 国产97色在线日韩免费| 九色亚洲精品在线播放| 国产av精品麻豆| 国产精品一区二区在线不卡| 99久久人妻综合| 亚洲av欧美aⅴ国产| 脱女人内裤的视频| 18禁观看日本| 久久99一区二区三区| 国产精品一区二区免费欧美 | 精品少妇黑人巨大在线播放| 日本wwww免费看| 男人爽女人下面视频在线观看| 99精品欧美一区二区三区四区| 桃花免费在线播放| 一区二区三区激情视频| 日本a在线网址| 久久热在线av| 日韩 欧美 亚洲 中文字幕| 久久久国产成人免费| 欧美性长视频在线观看| 成人国产一区最新在线观看| av网站免费在线观看视频| 啦啦啦在线免费观看视频4| 飞空精品影院首页| 成在线人永久免费视频| 我的亚洲天堂| 美女午夜性视频免费| 美女主播在线视频| 97精品久久久久久久久久精品| videosex国产| 国产精品欧美亚洲77777| 91字幕亚洲| 电影成人av| 精品高清国产在线一区| 亚洲视频免费观看视频| 国产成人精品无人区| 中文字幕高清在线视频| 人妻一区二区av| 久久久久久亚洲精品国产蜜桃av| 热re99久久精品国产66热6| 中文字幕人妻丝袜制服| 人人澡人人妻人| 欧美精品高潮呻吟av久久| 两性夫妻黄色片| 亚洲va日本ⅴa欧美va伊人久久 | 美女大奶头黄色视频| 一级毛片女人18水好多| 另类亚洲欧美激情| 大码成人一级视频| 悠悠久久av| 男女下面插进去视频免费观看| 美女扒开内裤让男人捅视频| 午夜两性在线视频| 叶爱在线成人免费视频播放| 国产精品一区二区免费欧美 | 精品一区二区三区四区五区乱码| 亚洲熟女毛片儿| 黑丝袜美女国产一区| 一级毛片电影观看| 天天躁狠狠躁夜夜躁狠狠躁| av福利片在线| 韩国高清视频一区二区三区| 亚洲七黄色美女视频| 久久综合国产亚洲精品| 男女无遮挡免费网站观看| 一区在线观看完整版| 欧美性长视频在线观看| 老司机影院毛片| 一级毛片电影观看| 国产xxxxx性猛交| 中文字幕人妻熟女乱码| 少妇 在线观看| 大陆偷拍与自拍| 久久久久精品国产欧美久久久 | 亚洲美女黄色视频免费看| 久久人人爽人人片av| 久久久国产欧美日韩av| 日本黄色日本黄色录像| 久久青草综合色| 一二三四社区在线视频社区8| 在线av久久热| 亚洲欧美色中文字幕在线| 久久久久久亚洲精品国产蜜桃av| 国产在线观看jvid| 亚洲精品国产一区二区精华液| 美女午夜性视频免费| 欧美一级毛片孕妇| 欧美黑人精品巨大| 久久久久精品国产欧美久久久 | bbb黄色大片| 国产高清视频在线播放一区 | videos熟女内射| av片东京热男人的天堂| 波多野结衣av一区二区av| 在线av久久热| 美女中出高潮动态图| 国产亚洲欧美精品永久| 国产av一区二区精品久久| 丰满人妻熟妇乱又伦精品不卡| 精品人妻在线不人妻| 午夜成年电影在线免费观看| 另类亚洲欧美激情| 久久影院123| 亚洲av美国av| 精品国内亚洲2022精品成人 | 亚洲专区中文字幕在线| 精品一区二区三区四区五区乱码| 无遮挡黄片免费观看| 中文字幕色久视频| 一区二区日韩欧美中文字幕| 午夜久久久在线观看| 久久精品国产亚洲av香蕉五月 | 成人国产av品久久久| 他把我摸到了高潮在线观看 | 建设人人有责人人尽责人人享有的| 在线观看免费日韩欧美大片| 久久久久久久大尺度免费视频| 久久香蕉激情| 国产精品影院久久| 成年女人毛片免费观看观看9 | 欧美性长视频在线观看| 亚洲国产欧美一区二区综合| 久久人人爽人人片av| 久久精品亚洲av国产电影网| 在线亚洲精品国产二区图片欧美| 亚洲国产av新网站| 国产在线观看jvid| 丰满人妻熟妇乱又伦精品不卡| 1024香蕉在线观看| 看免费av毛片| 黄色视频,在线免费观看| 一本大道久久a久久精品| 在线看a的网站| 男人爽女人下面视频在线观看| 国产精品 国内视频| 国产日韩欧美视频二区| 别揉我奶头~嗯~啊~动态视频 | 美女扒开内裤让男人捅视频| 宅男免费午夜| 免费一级毛片在线播放高清视频 | 国产麻豆69| 在线观看人妻少妇| 法律面前人人平等表现在哪些方面 | av福利片在线| 亚洲成人免费av在线播放| 亚洲精品乱久久久久久| 国产av一区二区精品久久| 国产日韩一区二区三区精品不卡| 国产精品九九99| 在线观看舔阴道视频| 黄频高清免费视频| 少妇 在线观看| 成年人黄色毛片网站| 三上悠亚av全集在线观看| 性少妇av在线| 啦啦啦免费观看视频1| 亚洲avbb在线观看| 国产91精品成人一区二区三区 | 妹子高潮喷水视频| 中国美女看黄片| 老司机影院成人| 在线观看一区二区三区激情| 男人舔女人的私密视频| 亚洲一区二区三区欧美精品| 午夜91福利影院| 日韩,欧美,国产一区二区三区| 狂野欧美激情性xxxx| 97在线人人人人妻| 免费不卡黄色视频| 自线自在国产av| 最黄视频免费看| 亚洲欧美激情在线| 久久人妻熟女aⅴ| 精品高清国产在线一区| 午夜福利,免费看| www.精华液| 精品亚洲成a人片在线观看| 别揉我奶头~嗯~啊~动态视频 | 久久国产精品大桥未久av| 一区二区三区四区激情视频| 人妻一区二区av| 国产片内射在线| 国产淫语在线视频| 国产精品亚洲av一区麻豆| 亚洲全国av大片| 男人操女人黄网站| 啦啦啦啦在线视频资源| 丝袜人妻中文字幕| 欧美日韩av久久| 亚洲七黄色美女视频| 一本久久精品| 精品久久久久久久毛片微露脸 | 久久精品国产a三级三级三级| 手机成人av网站| 波多野结衣一区麻豆| 午夜视频精品福利| 老司机午夜福利在线观看视频 | 亚洲,欧美精品.| 美女福利国产在线| 在线看a的网站| 亚洲成人免费电影在线观看| 99久久国产精品久久久| 国产精品香港三级国产av潘金莲| 人人妻人人添人人爽欧美一区卜| 免费高清在线观看日韩| 国产精品1区2区在线观看. | 亚洲男人天堂网一区| 亚洲色图综合在线观看| 国产亚洲精品久久久久5区| 亚洲欧美精品综合一区二区三区| 99精品久久久久人妻精品| 午夜福利视频在线观看免费| 国产高清videossex| 大香蕉久久网| 亚洲国产精品成人久久小说| 亚洲天堂av无毛| av在线老鸭窝| 欧美久久黑人一区二区| 俄罗斯特黄特色一大片| 欧美成人午夜精品| 少妇人妻久久综合中文| 日韩欧美一区视频在线观看| 久久青草综合色| 黄频高清免费视频| 亚洲精品在线美女| 成人国产av品久久久| 天天添夜夜摸| 久久国产精品影院| 国产极品粉嫩免费观看在线| 国产淫语在线视频| 不卡av一区二区三区| 在线观看免费视频网站a站| 99热网站在线观看| 丁香六月天网| 午夜激情久久久久久久| 日韩精品免费视频一区二区三区| 久久久精品94久久精品| 欧美变态另类bdsm刘玥| 午夜久久久在线观看| 欧美变态另类bdsm刘玥| 欧美精品av麻豆av| www.999成人在线观看| 亚洲av电影在线观看一区二区三区| 精品人妻熟女毛片av久久网站| 黄色毛片三级朝国网站| 亚洲性夜色夜夜综合| 伦理电影免费视频| 精品一品国产午夜福利视频| 精品一区在线观看国产| 亚洲av美国av| 久久国产亚洲av麻豆专区| 叶爱在线成人免费视频播放| 又黄又粗又硬又大视频| 大香蕉久久网| 精品久久久精品久久久| 久久精品国产亚洲av高清一级| 下体分泌物呈黄色| 99国产精品一区二区蜜桃av | 一本色道久久久久久精品综合| 精品人妻1区二区| 99热国产这里只有精品6| 老司机靠b影院| 99热国产这里只有精品6| 狠狠精品人妻久久久久久综合| 国产熟女午夜一区二区三区| 久久久久久久国产电影| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av成人不卡在线观看播放网 | 亚洲第一av免费看| 国产高清videossex| 久久久水蜜桃国产精品网| 亚洲av片天天在线观看| 国产精品自产拍在线观看55亚洲 | 国产有黄有色有爽视频| 夜夜骑夜夜射夜夜干| 亚洲av欧美aⅴ国产| 亚洲欧美日韩另类电影网站| 欧美+亚洲+日韩+国产| 免费看十八禁软件| 韩国高清视频一区二区三区| 黄片大片在线免费观看| 亚洲成人免费电影在线观看| 久久性视频一级片| 飞空精品影院首页| 人人妻人人澡人人看| 狠狠精品人妻久久久久久综合| 日本一区二区免费在线视频| 日韩欧美国产一区二区入口| 中文字幕人妻熟女乱码| 久久精品国产a三级三级三级| 99热网站在线观看| 欧美97在线视频| 男女高潮啪啪啪动态图| 亚洲第一青青草原| 免费在线观看日本一区| a 毛片基地| 亚洲欧美精品自产自拍| 亚洲av片天天在线观看| 亚洲欧美一区二区三区黑人| 又黄又粗又硬又大视频| 亚洲av日韩在线播放| 水蜜桃什么品种好| 免费看十八禁软件| 国产精品久久久久久精品古装| 亚洲av电影在线进入| 十分钟在线观看高清视频www| 欧美av亚洲av综合av国产av| av网站在线播放免费| 两个人免费观看高清视频| 国产高清国产精品国产三级| 精品少妇黑人巨大在线播放| 脱女人内裤的视频| 香蕉丝袜av| 免费在线观看视频国产中文字幕亚洲 | 99国产精品一区二区蜜桃av | 18禁裸乳无遮挡动漫免费视频| 一区福利在线观看| 丰满迷人的少妇在线观看| 日韩三级视频一区二区三区| 国产亚洲精品一区二区www | 午夜福利视频在线观看免费| 丰满人妻熟妇乱又伦精品不卡| 男女床上黄色一级片免费看| 久久精品国产综合久久久| av网站免费在线观看视频| 国产无遮挡羞羞视频在线观看| bbb黄色大片| 亚洲成人免费av在线播放| 亚洲avbb在线观看| 一本一本久久a久久精品综合妖精| 在线观看免费高清a一片| 亚洲精品中文字幕一二三四区 | 伊人久久大香线蕉亚洲五| 少妇人妻久久综合中文| 国产主播在线观看一区二区| 青春草亚洲视频在线观看| 国产在线免费精品|