張重英，王向陽

(1.蘭州石化職業(yè)技術(shù)學(xué)院土木工程學(xué)院，甘肅 蘭州 730060；2.甘肅省民航機場集團有限公司，甘肅 蘭州 730020)

洪水災(zāi)害是當今世界上損失最大的自然災(zāi)害。據(jù)聯(lián)合國統(tǒng)計，每年全世界60%的自然災(zāi)害是由于洪災(zāi)造成的。對于大型土木項目來說，防排洪系統(tǒng)至關(guān)重要。而且是一個復(fù)雜的開放系統(tǒng)[1- 2]，一方面，洪水自身具有隨機性、模糊性、灰色性、混沌性、分形分維性等種種不確定性；另一方面，防洪系統(tǒng)的排洪能力也具有不確定性，進一步加大了防洪系統(tǒng)的不確定性和風(fēng)險，更為主要的是，當今時代效益是建筑工程追求的目標，這也是在設(shè)計時要考慮的一個重要因素。傳統(tǒng)排洪溝的設(shè)計采用的是定值設(shè)計法[3]。即將排洪溝的排洪能力Qm看做是定值，只考慮實際洪水的概率分布，根據(jù)當?shù)氐谋┯曩Y料和設(shè)計基準期，由收集的洪水資料得出的不同重現(xiàn)期下的洪峰流量，只要設(shè)計的排洪溝的排水能力大于洪峰流量，就認為滿足要求。這樣計算并沒有實際反映結(jié)構(gòu)的可靠度，并且為保障排洪溝的安全性總是預(yù)留一定的安全水深高度，極有可能會造成資源的浪費。因此將可靠度的計算引入排洪溝的設(shè)計中，對其失效概率做出準確評價，得出最為合理的設(shè)計結(jié)果。

關(guān)于可靠度設(shè)計在防洪系統(tǒng)中的應(yīng)用，吳世偉等[4]和徐祖信等[5]分析了開敞式溢洪道水力設(shè)計中各種不確定因素，不僅考慮了洪水發(fā)生的概率分布，同時對溢洪道的斷面尺寸及施工造成的誤差進行統(tǒng)計，得出了溢洪道結(jié)構(gòu)的允許流量的概率分布(結(jié)構(gòu)抗力)。首次將結(jié)構(gòu)可靠度計算中廣泛采用的JC法用于泄洪風(fēng)險計算，并通過實例證明了JC法完全適合于泄洪風(fēng)險計算。但是由于當時條件限制，上述計算方法尚存在需要簡化改善的地方。本文以某機場的排洪溝設(shè)計為例，應(yīng)用簡化JC法和蒙特卡羅法對失效概率進行計算，并將計算結(jié)果進行對比分析，并在此基礎(chǔ)上提出新的計算思路。同時根據(jù)計算過程中遇到的水深取值問題進行探究分析，以求找到一種最為合理的排洪溝失效概率的計算方法。

1 排洪溝失效概率計算

1.1 排洪能力計算

采用可靠度的理論設(shè)計排洪溝，就要具體分析影響結(jié)構(gòu)可靠性的各個隨機變量，并求出它們的分布規(guī)律[7]。要計算出排洪能力的概率分布參數(shù)，需要獲得數(shù)據(jù)邊坡系數(shù)m=ctgθ、溝底寬b、溝深h、溝底坡度i、粗糙系數(shù)n的概率分布參數(shù)。

Qm=(bh+mh)5/3[b+2h(1+m2)0.5]-2/3i1/2/n

(1)

(2)

(3)

(4)

(5)

(6)

因此，可得如下結(jié)果：

(7)

1.2 設(shè)計洪水流量及其概率分布參數(shù)

通過計算得出設(shè)計洪水計算模式誤差和構(gòu)筑物防洪能力的均值和方差，計算排洪溝的失效概率，還需要設(shè)計洪水的概率分布參數(shù)。

收集了該地區(qū)氣象站從1951—1995年日雨量資料和1971—1995年自計雨量計資料(從10min到24h)，制定了該地區(qū)的暴雨公式。經(jīng)過分析計算獲得該地25年一遇暴雨雨力Sp=61.2(mm/h)，暴雨強度值系數(shù)0.761，所以25年一遇的暴雨公式為[7]：

a=61.2/t0.761(mm/h)

(8)

1.3 可靠度計算

我們已經(jīng)計算得出μ防洪能力=45.97m3/s，σ防洪能力=12.84m3/s。根據(jù)該機場的地理位置，查閱該地區(qū)的水文手冊，可以得到該機場排洪溝設(shè)計洪峰流量的計算模式誤差μ=0.72，σ=0.25。構(gòu)造排洪溝的防洪功能函數(shù)：

g(R，K，S)=R-KS

(9)

其中，防洪能力R服從正態(tài)分布N(45.97，12.84)，計算模式誤差K服從對數(shù)正態(tài)分布F(0.72，0.25)，實際洪水S服從極值Ⅰ型分布F(0.1233，17.82)，代入公式(9)，將實際洪水的分布進行當量正態(tài)化，采用JC法[8]得出可靠度指標β=1.846。因為實際洪水變量是以年最大值來分析的，所以該可靠度是每年的工程可靠度。參照可靠指標與失效概率的關(guān)系表(見表1)可知，該工程的年失效概率Pf=3.24×10-2。

為了檢驗結(jié)果的正確性，用蒙特卡羅法計算[9]，正如上面分析的三個變量：防洪能力、計算模式誤差、實際洪水。計算結(jié)果為年失效概率Pf=2.16×10-2。

在設(shè)計基準期為30年的運營過程中，總失事概率用下式計算：

Pf=1-(1-Pf′)N

(10)

式中，Pf—防洪工程在設(shè)計基準期內(nèi)總失事概率，Pf′—工程年失事概率，N—設(shè)計基準期的年數(shù)。所以用兩種方法計算得到該機場排洪溝在30年運營期間的總失事概率分別為Pf=1-(1-0.0324)30=0.628，Pf=1-(1-0.0216)30=0.481。

表1 可靠指標與失效概率的關(guān)系(×10-3)

用兩種算法求得的失效概率存在一定的誤差。蒙特卡羅法的結(jié)果可能因為迭代次數(shù)少的問題而存在一定的誤差，以及JC法在進行泰勒級數(shù)展開時略去二階及更高階項及當量正態(tài)化過程中會存在一定的誤差。兩種誤差的累積就造成了上面計算結(jié)果的誤差。

2 計算失效概率新方法

在上述方法中，是先求得了年失效概率，再用概率論的公式(10)計算得到基準期內(nèi)的失效概率。

可以直接考慮基準期內(nèi)的截排洪溝的排水能力，洪水流量的概率，來計算其失效概率。其中對于截排洪溝的排水能力，因為會有人員定期維護，并且排水能力主要與截排洪溝的截面積有關(guān)，在其發(fā)生坍塌時，會及時修復(fù)，所以，假定基準期內(nèi)的截排洪溝的排水能力R依然服從正態(tài)分布N(45.97，12.84)。計算模式誤差K與基準期的年限無關(guān)，仍服從對數(shù)正態(tài)分布F(0.72，0.25)。

主要考慮實際洪水S的分布變化。每年的實際洪水S服從極值Ⅰ型分布F(0.1364，17.82)，并且假定在基準期內(nèi)實際洪水相互獨立，于是按照統(tǒng)計學(xué)的原理，基準期內(nèi)實際洪水的概率分布函數(shù)為[10]

FsQT=[FsQτ]30

(11)

因為實際洪水S服從極值Ⅰ型分布，所以其概率密度函數(shù)為

FsQτ=exp{-exp[-αt(x-μt)]}

(12)

所以，基準期內(nèi)的實際洪水SQT也服從極值Ⅰ型分布，其概率分布中的參數(shù)為

(13)

由以上分析，經(jīng)過計算，基準期內(nèi)截排洪溝失效概率的三個相關(guān)隨機變量分布為：防洪能力R服從正態(tài)分布N(45.97，12.84)；計算模式誤差K服從對數(shù)正態(tài)分布F(0.72，0.25)；實際洪水S服從極值Ⅰ型分布F(0.1233，45.40)；

代入公式g(R，K，S)=R-KS，將實際洪水的分布進行當量正態(tài)化，采用JC法得出可靠度指標β=0.681，參照可靠指標與失效概率的關(guān)系表可知，該工程的在基準期內(nèi)的失效概率為Pf=0.248。

為了檢驗結(jié)果的正確性，用蒙特卡羅法編程計算，結(jié)果為Pf=0.218，兩種算法的計算結(jié)果基本一致，相差14%，說明結(jié)果具有一定的可靠性。同時，計算結(jié)果也說明了，對于相同分布的各個參數(shù)，當其取值不同時，JC法和蒙特卡羅法對于失效概率的計算的誤差會發(fā)生變化，這主要是由于參數(shù)發(fā)生變化時，相同抽樣次數(shù)抽到結(jié)構(gòu)失效的次數(shù)會變化，結(jié)果的精確度就會變化，這也就是我們采用重要抽樣法來提高模擬精度的原因。

對于失效概率的新方法與傳統(tǒng)的方法相比，結(jié)果相差很大。以JC法為例，結(jié)果見表2。

表2 兩種方法計算結(jié)果的差異

對于排洪溝在基準期內(nèi)的失效概率，新方法的計算結(jié)果更符合實際情況。因為要計算整個基準期內(nèi)的失效概率，就要以基準期內(nèi)的各個影響排洪溝失效概率的因素為研究對象來分析，而不是先計算每年的失效概率，再求出整個基準期內(nèi)的失效概率。同時從計算結(jié)果可以看出新方法的計算結(jié)果也符合實際情況，這樣的失效概率也是被人們所接受的。

3 計算水深對失效概率影響

在計算截排洪溝的排水能力時，用的是計算水深，而實際的截排洪溝發(fā)生失效時的水深為截排洪溝的總水深[11]，那么計算水深是否能滿足實際防洪的需要呢？為了弄清這個問題，就要研究水深與排洪能力和截排洪溝失效概率的關(guān)系。

分別選取計算水深如下，計算出截排洪溝的輸水能力和失效概率(實際洪水和計算模式誤差都不變)。計算結(jié)果如圖1所示。

圖1 計算水深對失效概率的影響

從圖1可以看出，隨著計算水深的增大，截排洪溝的排水能力基本呈線性增長，而失效概率卻下降很快，說明計算水深對結(jié)果(特別是失效概率)影響很大。同時考慮當計算水深變?yōu)闇仙顣r，排洪溝會有一定的淤積，當淤積過多時會有專人清理，所以溝深的均值及均方差都會發(fā)生變化，并且當溝深變化時，溝底寬度也會相應(yīng)變化。通過計算發(fā)現(xiàn)，當淤積造成的溝深值的均方差由0.15m增長到0.45m時，失效概率基本呈線性增長，說明計算水深的離散程度對截排洪溝的失效概率有一定影響，但沒有水深均值影響明顯。計算結(jié)果如圖2所示。

圖2 計算水深的離散程度對失效概率的影響

通過以上分析可以看出，排洪溝計算水深是影響其失效概率的一個主要因素，因此，當排洪系統(tǒng)的可靠度不能滿足要求時，首先考慮增大排洪溝的溝深。并且在日常維護中，要及時檢查排洪溝是否有淤積、溝壁是否發(fā)生破壞等直接影響水深的因素。

4 結(jié)語

(1)文章利用工程上常用的JC法和蒙特卡洛法對排洪溝的失效概率進行了計算對比和誤差分析，將計算年失效概率改為直接計算基準期內(nèi)總體失效概率，經(jīng)計算對比發(fā)現(xiàn)新方法的計算結(jié)果更容易被人接受。

(2)根據(jù)計算過程中遇到的水深取值問題進行分析，結(jié)果表明，排洪溝計算水深是影響其失效概率的一個主要因素，滿足計算的需要，最終確定較為合理的排洪溝失效概率計算方法。

(3)文章依然存在需要改進的地方，最主要的問題是新方法的計算結(jié)果雖然較為理想，但與傳統(tǒng)方法的計算結(jié)果之間差異較大，其中的原理依然有待研究，這也是下一步的研究重點。