許弘 程素珍 尹艷青

(1.蘭州大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，甘肅 蘭州 730107；2.山東省水利科學(xué)研究院，山東 濟(jì)南 250013；3.北京市水務(wù)工程建設(shè)與管理事務(wù)中心，北京 100036)

我國(guó)是海洋大國(guó)，擁有漫長(zhǎng)的海岸線，山東省陸域岸線長(zhǎng)度3122km，已建大陸海堤1184km，海堤與岸線比值為38%，且現(xiàn)有海堤破損嚴(yán)重，因此山東省已經(jīng)明確提出要加大海堤建設(shè)，提高沿海重點(diǎn)地區(qū)的防風(fēng)暴潮能力，而潮水位是確定海堤規(guī)模的關(guān)鍵性參數(shù)?！逗５坦こ淘O(shè)計(jì)規(guī)范》(SL 435—2008)(以下簡(jiǎn)稱《規(guī)范》)提出在海岸地區(qū)設(shè)計(jì)潮水位頻率分析的線型可采用極值Ⅰ型分布曲線[1-2]，但《規(guī)范》僅列出了頻率、數(shù)據(jù)序列項(xiàng)數(shù)等相關(guān)系數(shù)的表格，通過(guò)查表求解相應(yīng)設(shè)計(jì)頻率的潮水位。查表法計(jì)算簡(jiǎn)單、應(yīng)用方便，但其方法的局限性和成果的缺陷明顯，在大中型工程中的應(yīng)用受到較大的限制。適線法在水文頻率分析中已經(jīng)普遍和成熟應(yīng)用，但主要基于皮爾遜Ⅲ型分布的應(yīng)用[3-9]，有關(guān)極限Ⅰ型分布的工程應(yīng)用很少，現(xiàn)有規(guī)范、相關(guān)文獻(xiàn)也缺少極限Ⅰ型分布適線法的計(jì)算程序、適線方法和成果合理性判斷標(biāo)準(zhǔn)，使該方法在潮水位頻率分析應(yīng)用中受到較大的限制。本文基于matlab計(jì)算軟件，結(jié)合山東省榮成市石島灣海堤設(shè)計(jì)，開(kāi)發(fā)了極值Ⅰ型分布求矩適線法程序，并與規(guī)范表格系數(shù)查算、理論曲線成果進(jìn)行比較，驗(yàn)證了該法的合理性、便利性，以期對(duì)類似工程提供借鑒。

1 極值Ⅰ型分布原理和適線求解方法

根據(jù)《規(guī)范》規(guī)定，設(shè)計(jì)潮水位計(jì)算應(yīng)采用頻率分析法，在海岸地區(qū)頻率曲線的線型可采用極值Ⅰ型分布。

潮水位觀測(cè)站歷年最高潮水位觀測(cè)值為h1，h2，…h(huán)n，觀測(cè)值呈極值Ⅰ型分布，則某一設(shè)計(jì)潮位hP對(duì)應(yīng)的設(shè)計(jì)頻率值P為

P=P(h≥hP)=1-F(hP)=1-e-e-α(hP-β)

(1)

σ——潮水位序列均方差。

潮水位序列均值和均方差采用矩法進(jìn)行估算，見(jiàn)式(2)、式(3)。

(2)

(3)

式中hi——序列第i年的年最高潮水位；

n——年最高潮水位序列項(xiàng)數(shù)。

對(duì)式(1)兩邊取兩次對(duì)數(shù)，把α、β化簡(jiǎn)代入得式(4)，得

(4)

令：λPn=-(0.45+0.78×ln(-ln(1-P)))，式(4)轉(zhuǎn)化為：

(5)

根據(jù)極值Ⅰ型分布理論對(duì)實(shí)測(cè)經(jīng)驗(yàn)頻率曲線進(jìn)行適線擬合須借助數(shù)學(xué)工具。在matlab語(yǔ)句中，有evcdf、evpdf、evinv等函數(shù)，分別表示極值Ⅰ型分布的累計(jì)函數(shù)、密度函數(shù)和逆函數(shù)，可以通過(guò)P=evcdf(-hP,-b,1/a)、hP=-evinv(k,-b,1/a)實(shí)現(xiàn)對(duì)給定設(shè)計(jì)頻率和相應(yīng)設(shè)計(jì)頻率的設(shè)計(jì)水位進(jìn)行交互計(jì)算實(shí)現(xiàn)，實(shí)現(xiàn)頻率適線法。

2 資料分析

2.1 基本資料

圖1 設(shè)計(jì)潮水位極值Ⅰ型分布適線

2.2 表格法計(jì)算設(shè)計(jì)潮水位

根據(jù)《規(guī)范》，對(duì)石島海陽(yáng)觀測(cè)站35年連續(xù)的年最高潮水位序列的均值、均方差及年頻率為0.1%、0.2%、0.5%、1%、2%、4%、5%、10%、25%、50%、75%、90%、95%、97%、99%、99.9%的最高潮水位分別按式(2)、式(3)和式(5)計(jì)算，系數(shù)λPn按《規(guī)范》表A.0.2查取，不同頻率下的設(shè)計(jì)潮水位見(jiàn)表1，不同頻率及相應(yīng)潮水位點(diǎn)繪在海森概率格紙上并連成曲線(見(jiàn)圖1)。

表1 不同頻率下設(shè)計(jì)潮水位計(jì)算成果 單位：m

2.3 理論計(jì)算設(shè)計(jì)潮水位

通過(guò)理論計(jì)算知，曲線形狀由均值、均方差兩個(gè)參數(shù)決定。調(diào)整均值僅使頻率曲線平行上下移動(dòng)，偏移經(jīng)驗(yàn)點(diǎn)，通過(guò)調(diào)整均值進(jìn)行適線法是不合適的(見(jiàn)圖2)。但當(dāng)驗(yàn)潮零點(diǎn)高程的標(biāo)準(zhǔn)不一致時(shí)，可以通過(guò)調(diào)整序列均值實(shí)現(xiàn)。

圖2 調(diào)整均值的極值Ⅰ型適線

設(shè)定理論頻率曲線P=42.97%為分界點(diǎn)，該點(diǎn)的設(shè)計(jì)潮水位等于序列均值，高于該頻率的設(shè)計(jì)潮水位均大于序列均值，小于該頻率的設(shè)計(jì)潮水位均小于序列均值。由于理論曲線以分界點(diǎn)旋轉(zhuǎn)，因此通過(guò)調(diào)整均方差，能夠適配出合理的曲線。

2.4 求矩經(jīng)驗(yàn)適線法計(jì)算設(shè)計(jì)潮水位

求矩經(jīng)驗(yàn)適線法是在一定的適線準(zhǔn)則下，求解與經(jīng)驗(yàn)點(diǎn)據(jù)擬合最優(yōu)的頻率曲線的統(tǒng)計(jì)方法，可避免因觀測(cè)序列短而產(chǎn)生計(jì)算結(jié)果可靠性差等問(wèn)題。采用矩法對(duì)觀測(cè)數(shù)據(jù)序列進(jìn)行統(tǒng)計(jì)，求出一組參數(shù)作為初值計(jì)算理論曲線，根據(jù)理論頻率曲線與經(jīng)驗(yàn)頻率點(diǎn)據(jù)的配合情況，通過(guò)經(jīng)驗(yàn)判斷調(diào)整參數(shù)，選定一條與經(jīng)驗(yàn)點(diǎn)據(jù)擬合良好的頻率曲線。經(jīng)驗(yàn)頻率曲線選擇時(shí)，盡可能選擇擬合點(diǎn)據(jù)好、結(jié)果相對(duì)保守的曲線作為經(jīng)驗(yàn)適線。適線要盡可能照顧點(diǎn)群的趨勢(shì)，使頻率曲線通過(guò)點(diǎn)群中心，點(diǎn)線配合不佳時(shí)，可重點(diǎn)考慮中上部點(diǎn)據(jù)，并應(yīng)分析經(jīng)驗(yàn)點(diǎn)據(jù)的精度，使曲線盡可能接近或穿過(guò)比較可靠的點(diǎn)據(jù)。

通過(guò)初估、適線和綜合對(duì)比分析，可以得到比較合理的、能夠滿足海堤工程設(shè)計(jì)要求的潮水位經(jīng)驗(yàn)頻率曲線。通過(guò)選定的經(jīng)驗(yàn)適線，根據(jù)極值Ⅰ型理論計(jì)算相應(yīng)設(shè)計(jì)頻率的潮水位。借助經(jīng)驗(yàn)頻率曲線予以延長(zhǎng)，推求相應(yīng)稀遇頻率的設(shè)計(jì)值。

適線法計(jì)算通過(guò)matlab軟件中的evinv、evcdf函數(shù)實(shí)現(xiàn)，在軟件環(huán)境下編制程序，能夠?qū)崿F(xiàn)數(shù)據(jù)的統(tǒng)計(jì)、分析、繪圖和曲線擬合同時(shí)完成，采用norminv函數(shù)生成海森概率格紙。根據(jù)理論曲線和經(jīng)驗(yàn)點(diǎn)據(jù)的適配情況，調(diào)整均方差值為0.9σ、1.1σ、1.2σ、1.4σ倍的初值，進(jìn)行適配，曲線見(jiàn)圖1，根據(jù)每條適配線計(jì)算年頻率為0.1%、0.2%、0.5%、1%、2%、4%、5%、10%、25%、50%、75%、90%、95%、97%、99%、99.9%的最高潮水位，結(jié)果見(jiàn)表1。

3 計(jì)算成果和方法的合理性選定

表格法是《規(guī)范》推薦采用的設(shè)計(jì)潮水位計(jì)算方法，根據(jù)極值Ⅰ型分布，編寫了P-n-λPn參數(shù)表，計(jì)算簡(jiǎn)單，通過(guò)手算即可完成。極值Ⅰ型分布為全樣本理論曲線，為避免觀測(cè)系列較短對(duì)擬合計(jì)算結(jié)果的影響，采用觀測(cè)值序列項(xiàng)數(shù)多少調(diào)整λPn值，減小因系列短而產(chǎn)生的誤差。由圖1知，表格法計(jì)算成果和1.15倍均方差初值擬合曲線一致，說(shuō)明表格法由通過(guò)調(diào)整均方差初值計(jì)算生產(chǎn)，是適線中的一條經(jīng)驗(yàn)頻率曲線成果。

表格法把復(fù)雜的計(jì)算過(guò)程簡(jiǎn)化為查表計(jì)算，并把樣本序列長(zhǎng)度的影響考慮在內(nèi)，其計(jì)算成果能夠滿足小型工程的需要，但規(guī)范表格是在特定序列項(xiàng)數(shù)n、頻率P下查算相應(yīng)λPn值，對(duì)于表格中沒(méi)有列出的n、P對(duì)應(yīng)的λPn值，只能通過(guò)內(nèi)插求解，存在一定的誤差，同時(shí)由于成果是根據(jù)固定適配線確定，其合理性、與經(jīng)驗(yàn)點(diǎn)據(jù)的適配性存在較大的局限。

表格法計(jì)算的設(shè)計(jì)潮水位受序列長(zhǎng)度的影響較大，觀測(cè)序列的長(zhǎng)度不同，其λPn值不同，時(shí)間越短，其系數(shù)越大。為驗(yàn)證計(jì)算成果受時(shí)間序列長(zhǎng)短的影響程度，選取該工程20年、25年、30年、35年最大潮水位系列資料，同時(shí)計(jì)算20年一遇設(shè)計(jì)潮水位，其計(jì)算成果見(jiàn)表2，根據(jù)計(jì)算成果繪制的序列項(xiàng)數(shù)n與潮水位關(guān)系曲線見(jiàn)圖3。從圖3可知，序列長(zhǎng)度為25年的計(jì)算設(shè)計(jì)潮水位值最大，說(shuō)明表格法并不能保證序列長(zhǎng)度越小選定的成果越保守，這是該方法的局限性。

表2 表格法計(jì)算20年一遇潮水位成果 單位：m

圖3 序列項(xiàng)數(shù)n與潮水位關(guān)系曲線

適線法根據(jù)極值Ⅰ型分布和觀測(cè)資料統(tǒng)計(jì)參數(shù)進(jìn)行理論分析和經(jīng)驗(yàn)判斷，通過(guò)調(diào)整參數(shù)，適配出與經(jīng)驗(yàn)頻率點(diǎn)據(jù)擬合最好的曲線，得到滿足工程設(shè)計(jì)要求的不同設(shè)計(jì)頻率下的潮水位設(shè)計(jì)值，計(jì)算成果包含了表格法。它是基于matlab軟件環(huán)境編制的程序，能夠完成數(shù)據(jù)的統(tǒng)計(jì)分析，多條適配曲線擬合，圖表自動(dòng)生成，計(jì)算速度快，功能強(qiáng)大，能夠把經(jīng)驗(yàn)點(diǎn)據(jù)、規(guī)范表格計(jì)算值、適配曲線同時(shí)繪制在同一張海森概率格紙上，并把成果匯總于同一表格中，便于設(shè)計(jì)人員比較和分析，確定適合設(shè)計(jì)要求的設(shè)計(jì)潮水位計(jì)算值。

4 結(jié) 論

本文采用極值Ⅰ型分布，結(jié)合工程設(shè)計(jì)潮水位實(shí)例計(jì)算，得出以下結(jié)論：

a.規(guī)范表格法簡(jiǎn)單依靠參數(shù)表查算完成，過(guò)程機(jī)械，結(jié)果受資料序列長(zhǎng)度影響，且不能確定與經(jīng)驗(yàn)頻率點(diǎn)據(jù)的擬合程度，具有較大的局限性。

b.表格法計(jì)算成果與特定倍數(shù)均方差適線成果一致，能夠通過(guò)適線法計(jì)算完成。

c.適線法基于matlab軟件的極值Ⅰ型分布函數(shù)完成，計(jì)算理論符合《規(guī)范》要求，計(jì)算功能強(qiáng)大，能夠把經(jīng)驗(yàn)點(diǎn)據(jù)、《規(guī)范》表格計(jì)算值、適配曲線同時(shí)繪制在同一張海森概率格紙上，并把成果匯總于同一個(gè)表格中，便于設(shè)計(jì)人員比較和分析，確定適合設(shè)計(jì)要求的設(shè)計(jì)潮水位計(jì)算值。