大壩防洪安全的評估和校核

1、大壩防洪平安的評估和校核摘要：洪水漫壩風(fēng)險(xiǎn)失事，是影響大壩平安的主要原因之一。本文尋求一種既能反映大壩防洪系統(tǒng)隨機(jī)性和模糊性，又能合理描繪大壩防洪才能的風(fēng)險(xiǎn)模型，以實(shí)現(xiàn)對大壩防洪才能的定量化，進(jìn)而為已建大壩和待建大壩的防洪平安評估和校核創(chuàng)造條件.關(guān)鍵詞：大壩防洪平安洪水漫壩失事隨機(jī)微分方程風(fēng)險(xiǎn)分析目前我國和世界上約三分之一的大壩失事，是洪水漫壩所造成的.因此，正確地對大壩防洪平安進(jìn)展評估和校核，具有非常重要的意義.洪水漫壩風(fēng)險(xiǎn)是和大壩洪水設(shè)計(jì)標(biāo)準(zhǔn)嚴(yán)密聯(lián)絡(luò)的.按我國現(xiàn)行的洪水設(shè)計(jì)標(biāo)準(zhǔn)對大壩的防洪平安進(jìn)展分析，從水文角度估算的理論漫壩風(fēng)險(xiǎn)率遠(yuǎn)大于實(shí)際漫壩失事率.這說明現(xiàn)有大壩通常具有一定的抗洪潛力

2、.這一抗洪潛力主要來源于兩個方面：由于水文、水力等隨機(jī)不確定性的影響，導(dǎo)致了設(shè)計(jì)者在調(diào)洪演算過程和泄洪建筑物設(shè)計(jì)規(guī)模、壩頂高程的決策中，留有一定的平安系數(shù)；由于工程、管理等模糊不確定性的影響，導(dǎo)致了洪水漫壩風(fēng)險(xiǎn)失事臨界限值的模糊化，常使洪水位略超壩頂高程而不發(fā)生失事事故。目前，國內(nèi)外對大壩防洪平安的分析主要從洪水設(shè)計(jì)標(biāo)準(zhǔn)的選擇出發(fā)，僅能考察大壩防洪的水文風(fēng)險(xiǎn)，帶有一定的片面性和局限性.諸多隨機(jī)不確定性和模糊不確定性因素均未能定量引入對大壩防洪平安的分析中，致使這一問題迄今未能解決?，F(xiàn)有大壩防洪平安分析1.1已建大壩防洪平安水準(zhǔn)評估通過分析國內(nèi)外大壩洪水漫壩風(fēng)險(xiǎn)失事的統(tǒng)計(jì)資料，可以理解已建大壩的

3、總體防洪平安水準(zhǔn)及趨勢。據(jù)九十年代初統(tǒng)計(jì)我國共建有大壩83000余座，其中土壩占90%以上，運(yùn)行多達(dá)3040年.表1我國各類大壩的漫頂失事率類型座數(shù)n漫壩失事座數(shù)np漫壩失事率Rp(N=3040年)漫壩失事率(預(yù)測)Rp(N=50年)大型3582中型2480110.46%0.66%小型8001011341.42%2.02%總計(jì)8284811471.38%1.97%在這些大壩的運(yùn)行期中，發(fā)生漫壩風(fēng)險(xiǎn)失事的共有1147座，約占失事總數(shù)的46.6%.各類大壩的漫頂失事座數(shù)及風(fēng)險(xiǎn)率見表1.這一統(tǒng)計(jì)說明，我國大中型壩的實(shí)際漫迎風(fēng)險(xiǎn)率Rp約在0.46%；小型壩的Rp偏高些，達(dá)1.42%；總的Rp可達(dá)1.3

4、8%.八十年代是我國大壩失事最少的年代.據(jù)文獻(xiàn)1的資料，這一時期我國年均失事大壩27座，其中漫壩失事的約占40.2%.據(jù)此，以35年運(yùn)行期推算，其Rp約為0.46%.這一風(fēng)險(xiǎn)率遠(yuǎn)小于前述3040年累計(jì)均值Rp=1.38%。文獻(xiàn)24統(tǒng)計(jì)了世界各國失事大壩的資料，約三分之一是洪水漫壩而導(dǎo)致失事的.20000座符合ILD標(biāo)準(zhǔn)的大壩和12500座土壩的統(tǒng)計(jì)資料說明，在50年左右的運(yùn)行期中，其Rp分別為0.2%和0.36%。.大壩防洪平安的水文分析傳統(tǒng)的大壩防洪平安分析認(rèn)為：漫壩風(fēng)險(xiǎn)主要來自超標(biāo)洪水.大壩在整個運(yùn)行期間抗御洪水的平安度SN，取決于設(shè)計(jì)洪水重現(xiàn)期Tr和運(yùn)行年限N，有：(1)N年內(nèi)的漫壩風(fēng)險(xiǎn)

5、率PFN那么為：(2)這一分析粗略認(rèn)定，大壩遭遇超標(biāo)洪水必然漫壩失事，而未遇超標(biāo)洪水那么一定平安.根據(jù)我國現(xiàn)行的大壩洪水設(shè)計(jì)標(biāo)準(zhǔn)5，由式(1)求得的各類各等級大壩在N=35年和50年內(nèi)的SN，見表2.表2我國的洪水設(shè)計(jì)標(biāo)準(zhǔn)和相應(yīng)的水文平安度壩型等級大()大()中小()小()土石壩洪水標(biāo)準(zhǔn)Tr(年)1000020001000500200運(yùn)行年限N(年)35503550355035503550平安度SN(%99.6599.5098.2697.5396.5695.1293.2390.4783.9177.83混凝土壩洪水標(biāo)準(zhǔn)Tr(年)50001000500200100運(yùn)行年限N(年)355035503

6、55035503550平安度SN(%)99.3099.0096.5695.1293.2390.4783.9177.8370.3460.50從表1和表2的分析比照中可見，盡管過去我國大壩的洪水設(shè)計(jì)標(biāo)準(zhǔn)較現(xiàn)行規(guī)定為低、設(shè)計(jì)洪水計(jì)算成果也可能偏小，但實(shí)際的Rp遠(yuǎn)小于式(2)估算的水文PFN.這說明現(xiàn)有大壩通常都具有一定的抗洪潛力，僅從水文角度考察大壩的防洪平安是不夠的.大壩的防洪平安受諸多不確定性因素的影響，主要可歸納為：水文、水力等隨機(jī)不確定性和工程、管理等模糊不確定性.只有綜合考察這些不確定性因素的作用，才能正確理解大壩的抗洪潛力，定量確定漫壩失事率PFN。轉(zhuǎn)貼于論文聯(lián)盟.ll.漫壩失事的隨機(jī)風(fēng)

7、險(xiǎn)分析2.1隨機(jī)不確定性的作用大壩的洪水漫迎風(fēng)險(xiǎn)率是與水庫的整個調(diào)洪過程聯(lián)絡(luò)著的.在整個調(diào)洪過程中，存在著許多人們難以意料和控制的隨機(jī)不確定性因素，如：入庫洪水過程Q(t)的水文條件，出庫泄洪過程q(h,)的水力條件、庫容和水位關(guān)系(h)的邊界條件、防洪起調(diào)水位H(t0)的初始條件等的隨機(jī)性.所有這些導(dǎo)致了不同時刻水庫蓄洪量(t)的隨機(jī)變化，這一變化又制約著庫水位H(t)的隨機(jī)消長.無疑地，H(t)與大壩的漫迎風(fēng)險(xiǎn)率Pf親密相關(guān)。傳統(tǒng)的調(diào)洪演算方法不能考察H(t)過程的隨機(jī)性.所以，設(shè)計(jì)者不得不采取偏保守的設(shè)計(jì)，在調(diào)洪演算過程和泄洪建筑物設(shè)計(jì)規(guī)模、壩頂高程的決策中，留有一定的平安系數(shù).由此可見

8、，從隨機(jī)水文、水力設(shè)計(jì)角度分析，在遭遇超標(biāo)洪水條件下，大壩仍可具一定的防洪平安度.當(dāng)然，亦不排除未遇超標(biāo)洪水時的漫壩風(fēng)險(xiǎn)。2.2隨機(jī)微分方程的建立對隨機(jī)變量H(t)的定量分析，是確定Pf的關(guān)鍵。調(diào)洪過程隨機(jī)微分方程的建立，為模擬和分析H(t)的隨機(jī)變化創(chuàng)造了條件。文獻(xiàn)6根據(jù)調(diào)洪過程中(t)具有iener過程特性的分析，建立了帶有隨機(jī)作用項(xiàng)、隨機(jī)系數(shù)和隨機(jī)初始條件的隨機(jī)微分方程：,(3)式中，Q(t)為均值洪水過程線，q(H,)為均值泄洪過程線G(H)=d(H/dH,H)為均值水位庫容曲線，H為均值庫水位，H0為t0時刻的初始隨機(jī)庫水位.dB(t)/dt為一正態(tài)白噪聲，該項(xiàng)的增加表示了隨機(jī)因素的

9、引入.這一隨機(jī)輸入是在調(diào)洪過程中實(shí)現(xiàn)的.Q(t),q(H,)和(H)等的綜合作用，導(dǎo)致了(t)圍繞其均值過程線(t)作隨機(jī)游走，其概率密度服從正態(tài)分布.在扣除了(t)的均值偏移(t)后，存在一無偏的iener過程B(t)，其均值EB(t)=0，方差DB(t)=2t.B(t)的概率密度f(B)為：(4)式中，2為過程的強(qiáng)度，它決定著(t)的離散程，據(jù)多個隨機(jī)過程結(jié)合分布的概念及對B(t)量綱分析的結(jié)果，從3個隨機(jī)過程標(biāo)準(zhǔn)差Q(t),q(t)和(t)可得：(5)取2的均值作為式(4)的過程強(qiáng)度。式(3)可簡化為典型的It方程形式：(6)就考察洪水漫壩風(fēng)險(xiǎn)率而言，人們關(guān)注的是這一It方程解過程的概率

10、密度，即在不同時刻庫水位的隨機(jī)分布狀況.在式(6)的情況下，文獻(xiàn)7推導(dǎo)了FkkerPlank方程的簡化表達(dá)式，以求解庫水位過程的概率密度f(h,t):(7)式(7)是一確定性的偏微分方程，可采用有限差分的計(jì)算方法求解，以解得調(diào)洪過程各個時刻的f(h,t)和相應(yīng)的H(t),H(t).同時，還可求得相應(yīng)的均值泄量q(t)：(8)以隨機(jī)數(shù)學(xué)模型考察H(t)的不確定性，能綜合多種隨機(jī)因素在整個過程中的影響，具有較高的可信程度。2.3漫壩失事的隨機(jī)風(fēng)險(xiǎn)表述對于大壩的洪水漫迎風(fēng)險(xiǎn)，?？梢?guī)定明確的極限標(biāo)志，即z=H-D0,(9)式中，D為隨機(jī)的壩頂高程，其不確定性可據(jù)工程情況分析而定.由式(9)確定的Pf為

11、：在一定洪水頻率下，H(t)超過D事件發(fā)生的機(jī)率：Pf=PA/T=PfHtD=z0f(z)dz,(10)式中，A為H(t)大于D事件，T為一定頻率洪水發(fā)生的事件，PA/T為條件概率，f(z)為洪水超高z的概率密度函數(shù).只要確定了H和D的概型及相應(yīng)的統(tǒng)計(jì)特征參數(shù)，即可確定f(z)，進(jìn)而求得相應(yīng)的Pf值.圖1Pf,Pf和T的關(guān)系這樣，在不同的洪水重現(xiàn)期Ti所給出的洪水過程線Qi(t)條件下，可分別求得不同的Hi(t)和相應(yīng)的Pfi.由此，可繪得Pf-T關(guān)系曲線(見圖1).該曲線說明，大壩的隨機(jī)漫迎風(fēng)險(xiǎn)率Pf是隨T而變化的.在TTr的超標(biāo)洪水條件下，Pf不再認(rèn)定是1，而只是隨T的增加，逐漸趨向于1；

12、在TTr的條件下，也不再認(rèn)定絕無漫頂失事可能性，只是隨T的減小，漸趨于0；當(dāng)T=Tr時，Pf并未躍變，而是某一定值.顯然，由該曲線確定的大壩隨機(jī)風(fēng)險(xiǎn)率較傳統(tǒng)的水文風(fēng)險(xiǎn)方法更為合理。轉(zhuǎn)貼于論文聯(lián)盟.ll.3漫頂失事的模糊風(fēng)險(xiǎn)分析3.1模糊不確定性的作用洪水漫頂失事并沒有絕對清楚的界限，而具有著從量變到質(zhì)變的中介過渡過程.這一模糊的過渡過程是與大壩的工程和管理?xiàng)l件嚴(yán)密相關(guān)的，如：壩體的工程情況決定著大壩對洪水超越事件z的承受才能，混凝土壩、下游植被和防滲條件較好的土石壩，都可允許洪水位略超壩頂而不致失事破壞；在大洪水降臨時，采取搶險(xiǎn)加高壩頂?shù)却胧嗫擅獬赡馨l(fā)生的失事事故.所以，z與漫壩失事事件

13、A是相關(guān)而不一樣的，它可以表示為一個模糊子集A(z).顯然，A(z)的模糊性也直接影響著Pf。論文聯(lián)盟.LL.編輯。經(jīng)典數(shù)學(xué)無法研究和處理模糊現(xiàn)象.所以，在大壩防洪平安的分析中，只能忽略模糊不確定性的影響，簡單地以準(zhǔn)確數(shù)學(xué)方法處理，即認(rèn)為z發(fā)生時，A一定發(fā)生，從而造成Pf計(jì)算結(jié)果的片面性和局限性.由此可見，從工程和管理的角度出發(fā)，大壩的實(shí)際漫頂失事率Pf并不會象理論上估計(jì)的那么大.3.2模糊隸屬函數(shù)確實(shí)定隸屬函數(shù)是描繪模糊性的關(guān)鍵.為了考察模糊失事率Pf，必須首先分析特征因子z對A的隸屬函數(shù)A(z).它既要表達(dá)工程設(shè)計(jì)、管理人員的主觀意識對A的斷定和信度，又要受客觀條件的制約，成為客觀性的一種

14、量度.顯然，z越大，A的可能性就越大，對“洪水超高z大這一模糊概念，可采用升半正態(tài)的模糊分布8來描繪：(11)亦可以升半梯形分布粗略給出這一隸屬函數(shù)(12)式中，Z為超高臨界限值，需依實(shí)際工程和管理的情況，經(jīng)歷地確定.壩體允許漫溢條件好的、抗洪搶險(xiǎn)條件好的，都可采用較大的Z值.3.3拓展的漫壩失事模糊風(fēng)險(xiǎn)表述為了研究和處理漫壩失事過程中的模糊現(xiàn)象，必須將前述Pf定義拓展為：“發(fā)生H超過D事件，并在某種程度上造成大壩失事破壞的機(jī)率.即認(rèn)為，z僅是A的一個模糊子集A(z)，其模糊集合記為z0(13)在“漫壩風(fēng)險(xiǎn)失事率定義模糊拓展的根底上，根據(jù)模糊事件的概率計(jì)算方法，對式(10)修正，得到隨機(jī)模糊風(fēng)

15、險(xiǎn)率Pf：(14)假設(shè)考慮功能函數(shù)z=H-D服的正態(tài)分布N(z,z)，選取A(z)為常用的升半正態(tài)分布，即式(11)，那么有：(15)式(14)和(15)兼顧了漫壩失事的隨機(jī)性和模糊性，能較為全面、合理反映漫壩失事的實(shí)際風(fēng)險(xiǎn).圖1中Pf-T關(guān)系曲線說明，考慮了漫壩失事的模糊性，大壩的整個防洪平安度將有所增加.4大壩防洪平安分析方法的討論在水工領(lǐng)域，大壩的防洪風(fēng)險(xiǎn)度PFf通常指：在一定的壩高和泄洪建筑物設(shè)計(jì)規(guī)模下，大壩在運(yùn)行期N年內(nèi)，發(fā)生漫壩失事的機(jī)率.前述Pf,Pf僅為一定洪水發(fā)生條件下的失事率.在Pf-T圖上，根據(jù)計(jì)算精度需將T分作假設(shè)干區(qū)間，即可導(dǎo)出大壩一年運(yùn)行期PF1為：,i=1,2n(

16、16)其中，Ti可為非等值，但應(yīng)滿足：i=1,2n(17)這樣，在整個N年運(yùn)行期中PFN那么為：(18)顯然，大壩的防洪風(fēng)險(xiǎn)度PFN取決于N和PF1，而PF1那么由Pf-T曲線的位置所決定.該曲線主要由設(shè)計(jì)控制點(diǎn)Tr和相應(yīng)的Pf所制約.當(dāng)然，整條曲線的陡竣程度，即在假設(shè)干非設(shè)計(jì)控制點(diǎn)Ti上的Pfi，亦會對PF1產(chǎn)生一定的影響。選擇較大的Tr，將有利于降低PF1；同樣，在Tr點(diǎn)上Pf的減小，也將有利于降低PF1.在選定Tr的情況下，通過調(diào)洪演算所確定的壩高和泄洪建筑物必然要使大壩承受漫頂失事風(fēng)險(xiǎn).設(shè)計(jì)的保守程度將導(dǎo)致不同的Pf和Pf，從而使整條Pf-T曲線的位置變化.式(2)僅從水文角度分析，得

17、出的PFf其控制因素僅為Tr和N，顯得非常粗略.在此根底上，即可對詳細(xì)的大壩進(jìn)展防洪平安評估和校核.對已建大壩，可根據(jù)原設(shè)計(jì)給出的水文、水力等資料，補(bǔ)充分析隨機(jī)不確定性.同時，結(jié)合大壩建成后的實(shí)際，進(jìn)展模糊不確定性分析.由式(14)、(18)等計(jì)算確定PFf.只要對Pf所涉及的隨機(jī)模糊不確定性參數(shù)選擇得當(dāng)，計(jì)算并綜合統(tǒng)計(jì)各詳細(xì)大壩的PFN，總體受騙接近前述實(shí)際統(tǒng)計(jì)所得的Rp.當(dāng)某一大壩的PFN較Rp為大，那么說明該壩的防洪平安度達(dá)不到現(xiàn)有該類壩的一般水準(zhǔn)，當(dāng)采取措施予以改良.對新建大壩，可在現(xiàn)行洪水設(shè)計(jì)標(biāo)準(zhǔn)的根底上開展工作.即認(rèn)為，現(xiàn)行標(biāo)準(zhǔn)確定的Tr根本合理.根據(jù)工程實(shí)際，分析各種隨機(jī)和模糊不

18、確定性，計(jì)算確定PFN，使其滿足適當(dāng)?shù)姆篮槠桨捕纫?通過調(diào)整Pf，將使各類新建大壩具有較為統(tǒng)一的平安尺度.這樣，將使設(shè)計(jì)和管理中留有的平安富裕有據(jù)可依，從而做到平安和經(jīng)濟(jì)的統(tǒng)一.5算例和討論根據(jù)文獻(xiàn)7給出的某大壩情況作一算例，以說明本文建議的防洪平安分析方法.圖2H,Q、q和t的關(guān)系該壩為大()型混凝土壩，按Tr=1000a標(biāo)準(zhǔn)設(shè)計(jì).設(shè)計(jì)洪水過程線繪于圖2中，Qax=12400/s.T=20a時，Qax=6300/s;T=2000a時，Qax=16200/s.按放大典型洪水過程線方法，可分別給出其洪水過程線.該壩壩頂高程D=19.5(以堰頂高程為基準(zhǔn))，其泄洪建筑物有：主壩5孔、副壩1孔溢流

19、堰，另設(shè)1底孔.起調(diào)水位h0=12.5，當(dāng)入庫流量Q3000/s時，不泄洪，q=0；當(dāng)3000/sQ6300/s時，q=3000/s；當(dāng)Q6300/s時，洪設(shè)施敞泄。據(jù)此調(diào)洪規(guī)那么，采用前述隨機(jī)調(diào)洪演算方法，求得千年一遇洪水條件下的f(h,t)分布繪于圖3中.顯然，隨著過程的進(jìn)展，不同時刻的庫水位不再是定值，而是具有某種分布的隨機(jī)變量.據(jù)f(h,t)分布，可以求出相應(yīng)的H(t)和H(t)，以及q(t).不同情況下的H(t)，與傳統(tǒng)方法計(jì)算結(jié)果的比擬亦顯示在圖2中.圖3f(h,t)分布在這一計(jì)算過程中，取Q(t)=Q(t)/Q(t)=0.1，q=0.01，=0.05，以洪峰時刻相應(yīng)的值為計(jì)算根據(jù)

20、，按式(5)求得B(t)的=920000s.取初值H0為離散型的分布，即(19同時，還進(jìn)展了T=20a,2000a等情況下的隨機(jī)調(diào)洪計(jì)算.D符合正態(tài)分布，D=19.5,D=0.01.在前述調(diào)洪方案下，當(dāng)Tr=1000a洪水發(fā)生時，洪水漫壩的最大隨機(jī)風(fēng)險(xiǎn)出如今t=60h時，相應(yīng)Hax=18.65,Hax=1.18，由此求得Pfax=24.21%.此外，還計(jì)算求得T=20a,2000a等點(diǎn)的Pfax(見圖1).考慮該壩能具1以下的漫壩水頭承受才能或搶險(xiǎn)加高可能性，即取Z=1，那么以隨機(jī)模糊風(fēng)險(xiǎn)失事模型，計(jì)算得千年一遇洪水時的Pf=14.36%及相應(yīng)的Pf-T曲線.按式(16)積分求得PF1=0.0

21、0032相應(yīng)PF35=1.11%、PF50=1.59%。這兩數(shù)值遠(yuǎn)小于僅從水文角度考慮的PF1=0.001、PF35=3.44%和PF50=4.88%.在此根底上進(jìn)展的該壩防洪平安評估說明，考慮隨機(jī)模糊因素的作用后，該壩的PFN仍偏高，超過了國內(nèi)現(xiàn)有大中型壩Rp=0.46%和0.66%的水準(zhǔn).假設(shè)根據(jù)該壩設(shè)計(jì)洪水過程線的特點(diǎn)，改變調(diào)洪方案，將敞泄的控制限值Q降低至5500/s，敞泄的時間相應(yīng)提早到t=27h。這一方案下，ax=18.34,Hax=0.73，相應(yīng)的Pf和Pf大幅度降低，分別為6.83%和2.39%.同時，整條Pf-T曲線的位置降低，積分法求得PF1=0.000069，相應(yīng)的PF3

22、5=0.24%、PF50=0.34%。改變調(diào)洪方案后，該壩的防洪平安度將可超過我國現(xiàn)有大中型壩的水準(zhǔn).因此，當(dāng)預(yù)報(bào)有較大洪水出現(xiàn)時，該庫應(yīng)提早敞泄，并相應(yīng)做好抗洪搶險(xiǎn)準(zhǔn)備，以確保大壩的防洪平安.6結(jié)論洪水漫壩是大壩風(fēng)險(xiǎn)失事的主要原因之一.通過對國內(nèi)外失事大壩資料的統(tǒng)計(jì)分析，評估了現(xiàn)有大壩的總體防洪平安水準(zhǔn).在闡述大壩防洪系統(tǒng)隨機(jī)不確定性和模糊不確定性的根底上，建立了漫壩失事的隨機(jī)模糊風(fēng)險(xiǎn)模型.這一模型能較為合理地描繪大壩的防洪才能，定量給出漫壩失事的風(fēng)險(xiǎn)率.這樣，就為已建和待建大壩的防洪平安評估和校核創(chuàng)造了條件.參考文獻(xiàn)1張秀玲，文明宣。我國水庫失事的統(tǒng)計(jì)分析及平安對策討論.水利管理論文集.水利部水管司，1992.3.2NatinalAadeyPress.SafetyfExi