EBⅫ型分布曲線(xiàn)在北方干旱地區(qū)應(yīng)用

張 楠

(遼寧省營(yíng)口市水利勘測(cè)建筑設(shè)計(jì)院，遼寧 營(yíng)口 115000)

變化環(huán)境下，由于全球氣候變化及流域人類(lèi)活動(dòng)加劇的影響，我國(guó)防汛抗洪壓力逐年增大[1]。隨著氣候變化、氣溫普遍上升，全國(guó)范圍內(nèi)的極端氣象現(xiàn)象時(shí)常發(fā)生，部分局部地區(qū)暴雨、洪水等水文事件頻繁出現(xiàn)，往往造成嚴(yán)重災(zāi)害及重大經(jīng)濟(jì)損失[2]。受氣候變化及人類(lèi)活動(dòng)的影響，部分地區(qū)水文要素統(tǒng)計(jì)特性發(fā)生變化，以往的水文頻率分析成果在這些地區(qū)已存在不適用問(wèn)題。在實(shí)際應(yīng)用過(guò)程中發(fā)現(xiàn)，我國(guó)規(guī)范中推薦使用的P-Ⅲ型頻率曲線(xiàn)線(xiàn)型[3-6]對(duì)北方部分地區(qū)水文變量的擬合效果較差，其在北方干旱半干旱地區(qū)的適用性有待進(jìn)一步評(píng)估，同時(shí)也需要探索針對(duì)變化環(huán)境下適用于北方干旱半干旱地區(qū)水文頻率分布曲線(xiàn)線(xiàn)型。同時(shí)，隨著經(jīng)濟(jì)社會(huì)發(fā)展，人類(lèi)活動(dòng)改變自然環(huán)境的力度逐漸加大，在高強(qiáng)度人類(lèi)活動(dòng)影響下，天然流域下墊面狀況發(fā)生了劇烈變化，造成流域產(chǎn)匯流情況發(fā)生改變，因此需探究北方干旱半干旱地區(qū)流域降雨徑流水文非線(xiàn)性特征及其模擬，為北方干旱半干旱地區(qū)設(shè)計(jì)洪水研究提供依據(jù)。

1 研究方法

本文主要采用P-Ⅲ型頻率曲線(xiàn)、EBⅫ型分布、正態(tài)分布、對(duì)數(shù)正態(tài)分布、廣義正態(tài)分布5種分布曲線(xiàn)對(duì)北方干旱地區(qū)水文頻率曲線(xiàn)的影響進(jìn)行分析，主要介紹EBⅫ型分布曲線(xiàn)的計(jì)算原理，其他曲線(xiàn)計(jì)算可詳見(jiàn)參考文獻(xiàn)[7-10]。Burr系列分布來(lái)源于如下微分方程的解：

(1)

式中，g(x)—能使F(X)當(dāng)-∞

FEBⅫ(x)=1-{1+(x/b)c}-β(b,c,β>0)

(2)

EBⅫ型分布對(duì)BurrⅫ型分布進(jìn)行擴(kuò)展，令且，定義EBⅫ型分布函數(shù)和概率密度函數(shù)分別為：

(3)

和

(4)

式中，λ、c、k—分別為EBⅫ型分布尺度參數(shù)、形狀參數(shù)和不等式參數(shù)。當(dāng)k≤0時(shí)，0≤x≤∞；k>0時(shí),0≤x≤λ/k1/c。由于在洪水頻率分析中，洪水變量不會(huì)小于零，因此，應(yīng)用于洪水變量分析中的概率密度曲線(xiàn)左端不應(yīng)該出現(xiàn)負(fù)值。當(dāng)把EBⅫ型分布應(yīng)用于洪水頻率分析時(shí)，參數(shù)k應(yīng)該滿(mǎn)足約束k≤0。根據(jù)EBⅫ型分布概率密度函數(shù)形式，以下給出了參數(shù)λ=1.0時(shí)，參數(shù)c、k不同組合下EBⅫ分布6種典型的概率密度形狀曲線(xiàn)，如圖1所示。

圖1 EBⅫ分布典型概率密度形狀曲線(xiàn)

2 研究結(jié)果

2.1 P-Ⅲ分布參數(shù)對(duì)頻率曲線(xiàn)的影響

為分析P-Ⅲ分布曲線(xiàn)對(duì)參數(shù)EX對(duì)頻率曲線(xiàn)的影響。選定Cv=1.0、Cs=2.0，分析了參數(shù)EX分別為100、200、300、500以及700情況下P-Ⅲ分布的頻率曲線(xiàn)，為分析P-Ⅲ分布曲線(xiàn)對(duì)偏態(tài)系數(shù)Cs對(duì)頻率曲線(xiàn)的影響。選定EX=500、Cv=1.0，分析了了參數(shù)Cs分別為0.1、1.0、2.0、4.0、10.0情況下P-Ⅲ分布的頻率曲線(xiàn)，分析結(jié)果如圖2所示。

圖2 P-Ⅲ分布曲線(xiàn)對(duì)參數(shù)的影響結(jié)果

圖2為參數(shù)分別為0.5、1.0、1.5、2.0、3.0情況下P-Ⅲ分布的頻率曲線(xiàn)，從圖2(a)中可以得到：隨著增大，頻率曲線(xiàn)顯得越來(lái)越陡。從圖2(b)中可以得到：隨著逐漸變大，頻率曲線(xiàn)的上尾部分變得越來(lái)越陡，而下尾部分逐漸變得平緩。

2.2 EBⅫ分布參數(shù)對(duì)頻率曲線(xiàn)的影響

為分析EBⅫ分布的不等式參數(shù)k對(duì)頻率曲線(xiàn)的影響。選定c=2.0、λ=200，分析其分布參數(shù)對(duì)頻率曲線(xiàn)的影響，此外，為分析EBⅫ分布的尺度參數(shù)對(duì)頻率曲線(xiàn)的影響。選定c=2.0、k=-0.50，分析其分布參數(shù)對(duì)頻率曲線(xiàn)的影響，如圖3所示。

圖3 EBⅫ分布曲線(xiàn)對(duì)參數(shù)的影響結(jié)果

從圖3(a)中可以得到：在k<0的情況下，隨著k的減小，頻率曲線(xiàn)上端顯得越來(lái)越陡；不等式參數(shù)k對(duì)頻率曲線(xiàn)的上尾部分影響顯著，而頻率曲線(xiàn)的下尾部分(p>50%)幾乎不受參數(shù)k值變化的影響。從圖3(b)中可以得到：對(duì)于給定頻率，λ越大對(duì)應(yīng)的變量值越大；隨著λ的逐漸增大，頻率曲線(xiàn)整體變得越來(lái)越陡。

2.3 正態(tài)分布分布參數(shù)對(duì)頻率曲線(xiàn)的影響

為分析正態(tài)分布參數(shù)的均值EX對(duì)頻率曲線(xiàn)的影響。選定Cv=1.0，分析參數(shù)EX分別為200、400、600、800、1000情況下正態(tài)分布的頻率曲線(xiàn)，此外為分析正態(tài)分布參數(shù)的變差系數(shù)Cv對(duì)頻率曲線(xiàn)的影響，分析了參數(shù)Cv分別為0.5、1、1.5、2、3情況下正態(tài)分布的頻率曲線(xiàn)，分析結(jié)果如圖4所示。

圖4 正態(tài)分布分布曲線(xiàn)對(duì)參數(shù)的影響結(jié)果

從圖4中(a)可以得到：隨著EX逐漸變大，頻率曲線(xiàn)上尾部分抬升，頻率曲線(xiàn)逐漸變陡。從圖4(b)中可以得到：隨著Cv逐漸變大，頻率曲線(xiàn)逐漸變陡。

2.4 對(duì)數(shù)正態(tài)分布參數(shù)對(duì)頻率曲線(xiàn)的影響

為分析對(duì)數(shù)正態(tài)分布均值EX對(duì)頻率曲線(xiàn)的影響。選定Cv=1.0，分析了參數(shù)EX分別為50、100、200、500、1000情況下對(duì)數(shù)正態(tài)分布的頻率曲線(xiàn)，為分析對(duì)數(shù)正態(tài)分布變差系數(shù)Cv對(duì)頻率曲線(xiàn)的影響。選定EX=500，分析了參數(shù)Cv分別為0.2、0.4、0.5、1.0、2.0情況下正態(tài)分布的頻率曲線(xiàn)，分析結(jié)果如圖5所示。

圖5 對(duì)數(shù)正態(tài)分布分布曲線(xiàn)對(duì)參數(shù)的影響結(jié)果

從圖5(a)中可以得到：隨著EX逐漸變大，頻率曲線(xiàn)整體上升，頻率曲線(xiàn)逐漸變陡。從圖5(b)中可以得到：隨著Cv逐漸變大，頻率曲線(xiàn)上端抬升且逐漸變陡，而下端緩慢下降且趨于平緩。

2.5 廣義正態(tài)分布參數(shù)對(duì)頻率曲線(xiàn)的影響

為分析廣義正態(tài)分布的位置參數(shù)μ對(duì)頻率曲線(xiàn)的影響。選定σ=200、k=-1.0，分析了參數(shù)μ分別為100、200、500、1000、1500情況下GND分布的頻率曲線(xiàn)，為分析廣義正態(tài)分布的的形狀參數(shù)k對(duì)頻率曲線(xiàn)的影響。選定μ=200、σ=200，分析了參數(shù)k分別為-0.5、-0.8、-1.0、-1.5、-2.0情況下GND分布的頻率曲線(xiàn)，分析結(jié)果如圖6所示。

圖6 廣義正態(tài)分布分布曲線(xiàn)對(duì)參數(shù)的影響結(jié)果

從圖6(a)中可以得到：隨著σ值大的逐漸增大，頻率曲線(xiàn)上端上升，下端下降，且頻率曲線(xiàn)逐漸變陡。從圖6(b)中可以得到：隨著k值負(fù)的越來(lái)越大，頻率曲線(xiàn)上端逐漸變陡，而下端具有細(xì)微程度的抬升且變得平緩。

3 結(jié)語(yǔ)

(1)我國(guó)規(guī)范推薦的P-Ⅲ型頻率曲線(xiàn)在遼寧省部分地區(qū)或部分站點(diǎn)的適用性確實(shí)較差；同時(shí)，通過(guò)對(duì)5種頻率分布線(xiàn)型適用性的綜合研究比較發(fā)現(xiàn)，整體上EBⅫ型頻率曲線(xiàn)在北方干旱半干旱地區(qū)具有更好的適用性。在北方干旱半干旱水文頻率分析中可考慮使用EBⅫ型分布作為P-Ⅲ型分布的驗(yàn)證，對(duì)P-Ⅲ型分布擬合結(jié)果進(jìn)行復(fù)核。

(2)在P-Ⅲ型分布難以得到滿(mǎn)意的擬合效果情況下，如果EBⅫ型分布有更好的擬合效果，可以使用EBⅫ型分布進(jìn)行替代。