基于多重不確定性因素的施工導(dǎo)流風(fēng)險(xiǎn)分析

1、基于多重不確定性因素的施工導(dǎo)流風(fēng)險(xiǎn)分析第25卷第2期2007年4月青海大學(xué)(自然科學(xué)版)JournalofQinshaiUniversity(NatureScience)V01.25No.2Apr.2007基于多重不確定性因素的施工導(dǎo)流風(fēng)險(xiǎn)分析李假設(shè)東,胡志根2,徐森泉2,蒲建平3(1.青海大學(xué),青海西寧810016;2.武漢大學(xué),湖北武漢430072;3.成都勘測(cè)設(shè)計(jì)院,四川成都610072)摘要:通過(guò)分析施工洪水過(guò)程中的洪峰流量,洪量和歷時(shí)的不確定性及導(dǎo)流建筑物的水力參數(shù)的不確定性,確定施工洪水過(guò)程和導(dǎo)流建筑物泄流能力的模擬計(jì)算模型.運(yùn)用MonteCarlo方法模擬計(jì)算施工導(dǎo)流系統(tǒng)渡訊各時(shí)

2、段的狀態(tài)參數(shù),得到上游圍堰堰前水位的分布函數(shù),在此基礎(chǔ)上確定導(dǎo)流系統(tǒng)的風(fēng)險(xiǎn).關(guān)鍵詞:施工導(dǎo)流;風(fēng)險(xiǎn)分析;MonteCarlo方法中圖分類號(hào):TV854文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):10068996(20Cr7)o2005704礎(chǔ)skanalysisforconstructiondiversionbasedonthemultiuncertaintyLIRuodong1,HUZhigen2,XUSenqua,PUJianpiIl(1.QinghaiUniversity,Xining810016,China;2.SchoolofWaterReseureeandHydropowerEngineering,W

3、uhanUniversity,Wuhan430072,China;3.ChengduSurveyandDesignInstitute,Chengdu610072,China)Abstract:Basedonlyinguncertaintiesofconstructiondiversionsystem,thispaperpresentsthesimulationmodelofconstructionfloodanddischargecapacityofdiversionstructure,andthemodelincludesuncertaintiessuch/ISpeakflow,floodv

4、olumes,flooddurationandhydraulicparametersofdischargeconstruction.3"heMonteCarlomethodisapFdjedtosimulatethestateparametersoftheconstructiondiversionsystemineveryperiodoftime.Asaresult,arandomdistributionfunctionoftheupstreamwaterlevelisobtained,thentheriskofthediversionsystemisdetermined.Keywo

5、rds:constructiondiversion;riskanalysis;MonteCarlomethod施工導(dǎo)流的風(fēng)險(xiǎn)分析在施工導(dǎo)流工程中占有重要的地位.由于施工導(dǎo)流系統(tǒng)的復(fù)雜性,不確定性因素有很多,如水文的不確定性,水力的不確定性和其它數(shù)據(jù)的不確定性等.國(guó)內(nèi)外不少學(xué)者在施工導(dǎo)流風(fēng)險(xiǎn)分析方面做了很多探討,提出一些導(dǎo)流風(fēng)險(xiǎn)模型1.如采用概率組合的方法間接估算施工導(dǎo)流風(fēng)險(xiǎn),模型中只考慮了洪水過(guò)程的不確定性,在實(shí)際工程中流量資料有限的情況下難以得到理想的結(jié)果l2;基于隨機(jī)微分方程分析調(diào)洪過(guò)程庫(kù)水位的隨機(jī)變化規(guī)律,求解庫(kù)水位的概率密度分布l3j3,對(duì)于施工導(dǎo)流工程,微分方程的求解存在困難.近來(lái)又

6、有學(xué)者提出利用Monte.Carlo方法,通過(guò)統(tǒng)計(jì)分析模型確定動(dòng)態(tài)風(fēng)險(xiǎn)l4J.但是在模擬施工洪水過(guò)程中,只考慮了洪峰流量的不確定性,而沒(méi)有洪水過(guò)程的洪量和歷時(shí)的不確定性,以及導(dǎo)流建筑物的各種水力參數(shù)的不確定性.為此,本文基于多重不確定性因素,運(yùn)用Monte.Carlo方法,模擬施工洪水過(guò)程和導(dǎo)流建筑物泄流能力,綜合分析施工導(dǎo)流系統(tǒng)的風(fēng)險(xiǎn).1施工洪水過(guò)程的不確定性模擬計(jì)算洪水過(guò)程變化受眾多因素影響,極其復(fù)雜且有明顯的隨機(jī)性.主要表現(xiàn)在:施工洪水過(guò)程洪峰流量的不確定性;施工洪水過(guò)程洪量的不確定性;施工洪水過(guò)程歷時(shí)的不確定性.收稿日期:20061114作者簡(jiǎn)介:李假設(shè)東(19曉一),男,青海西寧人,

7、副教授.研究方向:施工導(dǎo)流.58青海大學(xué)第25卷1.1洪峰流量和洪量的不確定性分析洪水過(guò)程受多種復(fù)雜因素影響,這樣的隨機(jī)變量總體分布是未知的,洪峰流量和洪量是洪水過(guò)程不確定性中的兩個(gè)重要因素,經(jīng)多年實(shí)踐證明,在我國(guó)一般假定為P一型是合理的.施工導(dǎo)流風(fēng)險(xiǎn)分析中,假定洪峰流量和洪量服從P一型分布,其密度函數(shù)為:)=(一e(1)式中,口,和口0為P一11I型分布的形狀,刻度和位置參數(shù);/-,(口)為口的伽瑪參數(shù).并且,口=4/2=2/At9coc,口o=.(12co/c,);c|為PHI型分布的離差系數(shù);co為P一型分布的離勢(shì)系數(shù);.為P一型分布的均值.1.2洪水過(guò)程歷時(shí)的不確定性分析當(dāng)前我國(guó)的大型

8、水利水電工程主要集中在西南部,大局部工程不具備獲得多年的實(shí)測(cè)洪水過(guò)程資料.由于洪水過(guò)程歷時(shí)也存在著不確定性,現(xiàn)就某工程歷年實(shí)測(cè)洪水過(guò)程資料(選用其中14年)來(lái)初步確定洪水過(guò)程歷時(shí)的密度函數(shù)分布.由于水利水電工程的實(shí)測(cè)洪水過(guò)程資料一般較少,所以采用w檢驗(yàn)法來(lái)檢驗(yàn)J(表1).表1洪水歷時(shí)分析k(I)(+1一I)(+1一I)一(I)akk(I)(+1一I)(n+1一I)一(I)a4328849l26181260.116224107203100.299295106009o0.096534396972580.251010510570.600.o77844596662070.21

9、5ll152O569490.059854626361740.185712522564420.O42464626241620.160025374656241590.137214537552150.0084通過(guò)統(tǒng)計(jì)檢驗(yàn)?zāi)Ｐ陀?jì)算得W=O.9872,Wo.05=O.924,所以認(rèn)為洪水過(guò)程歷時(shí)服從正態(tài)公布,其密度函數(shù)為:f(x)=1e(2)式中u為歷時(shí)的數(shù)學(xué)期望,d2為歷時(shí)的方差.1.3施工洪水過(guò)程的綜合模擬計(jì)算綜合考慮施工洪水過(guò)程中洪峰流量,洪量和歷時(shí)的不確定性,運(yùn)用MonteCarlo方法對(duì)洪峰流量,洪量和歷時(shí)進(jìn)行隨機(jī)模擬,采用分時(shí)段同頻率放大典型洪水過(guò)程線,結(jié)合歷時(shí)調(diào)整

10、洪水過(guò)程線的方法,模擬出一系列供導(dǎo)流風(fēng)險(xiǎn)分析計(jì)算的施工洪水過(guò)程.2導(dǎo)流建筑物泄流能力的不確定性模擬計(jì)算導(dǎo)流系統(tǒng)泄流能力的不確定性和泄流建筑物水力參數(shù)的不確定性有關(guān),如糙率,過(guò)水?dāng)嗝婷娣e,濕周,底坡等.導(dǎo)流建筑物的施工,運(yùn)行過(guò)程中存在誤差等是主要導(dǎo)致水力參數(shù)不確定性的原因.2.1泄流能力和水力參數(shù)的關(guān)系滿寧方程廣泛應(yīng)用于隧洞,明渠等,對(duì)于導(dǎo)流建筑物,其泄流能力可采用修正因子入對(duì)滿寧方程進(jìn)行修正,那么泄流能力Q為J:?生?A5/3?P一2/3?$1/2(3)式中/7,為滿寧糙率系數(shù);A為過(guò)水?dāng)嗝婷娣e;P為濕周;S為底坡.水力參數(shù)n,A,P,S等具有不確定性,如果各水力參數(shù)的密度函數(shù)為,根據(jù)泄流能力

11、和水力參數(shù)的關(guān)系,理論上就可以得出泄流能力的密度函數(shù).2.2水力參數(shù)的不確定性分析現(xiàn)在國(guó)內(nèi)外缺乏大量資料來(lái)準(zhǔn)確確定這些水力參數(shù)的隨機(jī)分布概率,所以常對(duì)它們的分布進(jìn)行假設(shè).對(duì)于泄流建筑物為隧洞時(shí),由施工或者試驗(yàn)引起過(guò)水?dāng)嗝婷娣e,濕周和底坡不確定性的因素有很多,假定過(guò)水?dāng)嗝婷娣e,濕周和底坡服從正態(tài)分布;糙率根據(jù)一些工程經(jīng)驗(yàn)和學(xué)者的分析探討,假定服從三角形分布.2.3泄流能力的不確定性分析模擬根據(jù)上面假設(shè)的各水力參數(shù)的隨機(jī)分布,直接確定泄流能力的密第2期李假設(shè)東等:基于多重不確定性因素的施工導(dǎo)流風(fēng)險(xiǎn)分析59度函數(shù)比擬困難.分析具體工程導(dǎo)流洞的泄流能力時(shí),可以運(yùn)用MonteCarlo方法模擬各水力參數(shù)

12、,通過(guò)假設(shè)檢驗(yàn)分析泄流能力的隨機(jī)分布.圖1是模擬某工程導(dǎo)流洞泄流能力統(tǒng)計(jì)得到的統(tǒng)計(jì)密度函數(shù)曲線和擬合的三角形分布密度函數(shù)曲線.通過(guò)對(duì)圖1的統(tǒng)計(jì)密度函數(shù)進(jìn)行檢驗(yàn),證明三角形分,n,布優(yōu)于正態(tài)分布,認(rèn)為導(dǎo)流洞的泄流能力較為接近三角形0腫25分布,其密度函數(shù)為:)=口b6<c(4)其他式中口為泄流能力下限;b為平均泄流能力;C為泄流能力上限;口,b,C參數(shù)通過(guò)模擬泄流能力的統(tǒng)計(jì)資料來(lái)確定.71049710987Q(m)圖1泄流能力的密度函數(shù)曲線由于水利水電工程施工導(dǎo)流系統(tǒng)的復(fù)雜性,必然存在其他的不確定性因素,如上游庫(kù)區(qū)現(xiàn)場(chǎng)測(cè)量和繪制的誤差,上游運(yùn)行多年的泥沙淤積和坍塌等.限于現(xiàn)階段這方

13、面資料的搜集缺乏和研究不深人,導(dǎo)流風(fēng)險(xiǎn)計(jì)算暫不考慮這局部不確定性因素.3基于多重不確定性因素的導(dǎo)流風(fēng)險(xiǎn)計(jì)算施工導(dǎo)流的風(fēng)險(xiǎn)是根據(jù)設(shè)計(jì)資料,綜合考慮施工洪水過(guò)程和導(dǎo)流建筑物泄流能力等不確定性因素的影響,模擬上游洪水過(guò)程和導(dǎo)流建筑物泄流能力,通過(guò)施工洪水的調(diào)洪演算,統(tǒng)計(jì)分析計(jì)算上游圍堰堰前水位超過(guò)上游圍堰設(shè)計(jì)擋水位的概率,即R=尸()式中為模擬上游圍堰堰前水位;爿加為上游圍堰設(shè)計(jì)擋水位.確定施工導(dǎo)流風(fēng)險(xiǎn),要系統(tǒng)考慮河道來(lái)流洪水過(guò)程,導(dǎo)流建筑物泄流能力及樞紐的其他特征,通過(guò)系統(tǒng)模擬的方法來(lái)實(shí)現(xiàn).綜合考慮多重不確定性因素,運(yùn)用MonteCarlo方法,模擬施工洪水過(guò)程及導(dǎo)流建筑物泄流能力,通過(guò)調(diào)洪演算和

14、統(tǒng)計(jì)分析模型得到導(dǎo)流設(shè)計(jì)水位的風(fēng)險(xiǎn)率.基于多重不確定性因素的導(dǎo)流風(fēng)險(xiǎn)系統(tǒng)模擬流程見(jiàn)圖2.4結(jié)語(yǔ)圖2導(dǎo)流風(fēng)險(xiǎn)系統(tǒng)模擬流程考慮導(dǎo)流系統(tǒng)中多重不確定性因素,本文在考慮施工洪水過(guò)程的不確定性中,除考慮了洪峰流量和洪量,把洪水歷時(shí)也考慮進(jìn)去;通過(guò)分析各水力參數(shù)的不確定性,確定泄流能力的分布函數(shù).運(yùn)用MonteCarlo方法,模擬施工洪水過(guò)程和導(dǎo)流建筑物的泄流能力,并通過(guò)調(diào)洪演算和統(tǒng)計(jì)分析模型確定施工導(dǎo)流的風(fēng)險(xiǎn)率.分析了基于多重不確定性因素的導(dǎo)流風(fēng)險(xiǎn)計(jì)算模型的可行性和有效性.由于水利水電施工導(dǎo)流影響因素很多,隨機(jī)因素不容易確定,文中提出的導(dǎo)流風(fēng)險(xiǎn)分析方法還需要不斷地改良.參考文獻(xiàn):1肖煥雄,韓采燕.施工導(dǎo)

15、流系統(tǒng)超標(biāo)洪水風(fēng)險(xiǎn)率模型研究J.水利,1993,24(11):76832王卓甫.考慮洪水過(guò)程不確定的施工導(dǎo)流風(fēng)險(xiǎn)分析計(jì)算J.水利,1998,24(4):3337.告島二(二(二二魚(yú)口旦,0,址青海大學(xué)第25卷3姜樹(shù)海.隨機(jī)微分方程在泄洪風(fēng)險(xiǎn)分析中的應(yīng)用J.水利,1994,25(3):19.4胡志根,劉全,賀昌海,等.基于Monte-Carlo方法的土石圍堰擋水導(dǎo)流風(fēng)險(xiǎn)分析J水科學(xué)進(jìn)展,2002,22(9):6346385王燕生.工程水文學(xué)M.北京:水利電力出版社,1992.4145.6長(zhǎng)江水利委員會(huì).水利水電工程設(shè)計(jì)洪水計(jì)算手冊(cè)M.北京:水利電力出版社,1995.3236.7肖煥雄,韓采燕.施

16、工導(dǎo)流標(biāo)準(zhǔn)與方案優(yōu)選M.武漢:湖北科學(xué)技術(shù)出版社,1996.182o.(責(zé)任編輯王寶通)(上接第51頁(yè))C1試件中可見(jiàn)卷席狀的水化產(chǎn)物層層疊起.C2試件中的水化產(chǎn)物形成密不可分的結(jié)構(gòu),幾乎連接在一起.cl試件(28d.3.00KX)c2試件(28d.3.00KX)圖6不I司配比砂漿水化28d的SEM圖4結(jié)論(1)水泥顆粒的粒徑分布與最緊密堆積理論的要求尚有很大差距,因此摻人比水泥顆粒更細(xì)的礦物摻合料,可以起到改善膠凝材料顆粒級(jí)配的作用.(2)通過(guò)試驗(yàn)證明,不同細(xì)度的礦物摻合料復(fù)合摻配時(shí),它們的顆?？梢韵嗷パa(bǔ)充,使膠凝材料顆粒的粒徑分布更接近緊密堆積狀態(tài)的要求,從而提高水泥漿體的密實(shí)度.(3)砂漿中膠凝材料的顆粒級(jí)配在接近Andreasen顆粒最緊密堆積理論要求時(shí),砂漿強(qiáng)度較高,水泥凝膠體的微觀結(jié)構(gòu)也比擬致密,這時(shí)可以適當(dāng)增加礦物摻合料替代水泥的百分率,而不會(huì)明顯影響強(qiáng)度.參考文獻(xiàn):1蒲心誠(chéng),王勇威.高效活性礦物摻料與混凝土的高性能化J.混凝土,2002,(2):36.2李瀅,楊靜.膠凝材料顆粒級(jí)配對(duì)水泥凝膠體結(jié)構(gòu)及強(qiáng)度的影響