鄒希云,莫松柏,張家望,張選民,鄒 胤
(湖南省益陽市氣象局,湖南 益陽 413000)
基于橋墩歷史水痕與氣象預(yù)報的中小橋水害預(yù)測
鄒希云,莫松柏,張家望,張選民,鄒 胤
(湖南省益陽市氣象局,湖南 益陽 413000)
該文將中小橋山洪水害作為研究對象,以人工神經(jīng)網(wǎng)絡(luò)中的BP網(wǎng)絡(luò)為工具,將橋墩歷史上的數(shù)次洪水水痕高程和產(chǎn)生當(dāng)次洪水的雨量記錄作為數(shù)據(jù)基礎(chǔ),提出了一種全新的流量計算和預(yù)報既有中小橋水害的方法。將“流域和橋渡”看成一個系統(tǒng),直接建立起輸入 (雨量)與輸出 (水位)的數(shù)學(xué)關(guān)系。因此,它不需要勘測流域面積、主河道長度及平均坡度等工作,可以節(jié)約大量的人力物力。多座橋梁的具體應(yīng)用表明,這種方法具有較高的準(zhǔn)確度,是一種有廣泛應(yīng)用前景的新方法。
氣象預(yù)報;橋墩水痕;人工神經(jīng)網(wǎng)絡(luò);流量計算
水毀橋梁是世界各國橋梁面臨的嚴重的自然災(zāi)害[1]。據(jù)統(tǒng)計,美國國家橋梁數(shù)據(jù)庫中有 58萬座橋梁,建在河流上的大約占 86%。1991年瑟赫爾 (Shirhole)和何特 (Holt)對美國已毀壞的 823座橋梁的毀壞原因進行分析,發(fā)現(xiàn)水毀橋梁占毀壞橋梁的 60%。同樣,在我國橋梁遭受洪水破壞也比較嚴重[1]。據(jù)有關(guān)部門不完全統(tǒng)計,僅 1977年9個省 (區(qū))就有水毀橋梁 900多座,水毀涵洞2 700多座,水毀損失約 5 400萬元;1981年四川發(fā)生嚴重的洪水災(zāi)害,全省公路和橋梁遭到了巨大的損失和破壞,共有 593條省市公路斷道阻車,沖毀和局部毀壞橋梁 457座,計 1 298m,沖毀涵洞4 647座,經(jīng)濟損失達 6 200多萬元。
對于大橋、特大橋,我國有專業(yè)部門對這些橋梁所跨過的河流作專業(yè)水文工作,這些橋梁發(fā)生水害的次數(shù)越來越少。但占 80%以上的中小橋卻沒有專門的水文觀測站,這些橋梁的水害已經(jīng)成了交通運輸?shù)耐怀鰡栴}。
引起中小橋梁水害的因素很多,有設(shè)計時對水文參數(shù)了解不夠的因素;有運用后橋梁環(huán)境 (如植被)發(fā)生變化的因素;有非法采砂等人為因素等等。如何對既有中小橋梁的抗洪能力作出簡單有效的評估;如何在山洪發(fā)生前對中小橋即將面對的狀況作出有效的預(yù)測,以便合理安排抗洪物資和抗洪隊伍,是人們十分關(guān)心的問題。
我國《鐵路橋梁檢定規(guī)范》[2]對鐵路橋梁水文檢定作出了原則上的規(guī)定。現(xiàn)在對中小橋梁檢定流量的估計還是采用由雨量到流量的方法,即根據(jù)橋梁所在地區(qū)的實際情況,采用地區(qū)經(jīng)驗公式法、推理公式法、單位線法、鐵路設(shè)計院法 (以下稱傳統(tǒng)方法)等方法[3],其中推理公式法又是各地區(qū)普遍采用的方法。
推理公式法認為,橋址河床斷面的流量與這個斷面的匯水面積、暴雨強度的關(guān)系是
式中,Q為通過橋址河床斷面的流量;αp為橋址上游匯水面積上的暴雨強度;A為橋址河床斷面匯水面積;K為比例系數(shù)(它與匯流方式、雨量損失、流域坡度、主河槽長度等因數(shù)有關(guān))。
利用推理公式法公式和其它傳統(tǒng)方法計算檢定流量,需要確定匯水面積、流域坡度、主河槽長度等,這些工作勞動強度大,收集時困難較多。
橋梁水文計算產(chǎn)生水位關(guān)系曲線,它確定了橋梁流量與水位的關(guān)系,如圖 1。
圖 1 金雞灣橋的水位流量關(guān)系曲線
令橋梁水位與流量的關(guān)系曲線為
雖然這是一條曲線,但在梁底附近的高水位區(qū)間,可以認為流量與水位是線性關(guān)系,如圖 1。將式(1)代入 (2),有
對于特定橋梁,式 (3)中的匯流面積 A是常數(shù);流域坡度、主河槽長度是常數(shù);與比例因子 K有關(guān)的匯流參數(shù)、損失參數(shù)也接近常數(shù)。因此,可以認為在高水位區(qū)間,洪水水位 H與暴雨強度αp接近線性關(guān)系。
橋梁在歷年洪水作用下,墩臺上不可避免的保留有各次洪水的記錄(水痕),由于混凝土中部分元素溶于水,混凝土橋墩的水痕更加明顯。
表 1 京廣線橋 1538歷史水痕實測與記錄對照表
表 1為京廣線橋 1538歷史水痕實測與記錄對照表。表 1表明:歷史洪水的當(dāng)時記錄與現(xiàn)時水痕實測值的誤差較小,利用實測水痕是可靠的。
人們最關(guān)心高水位區(qū)間的洪水,水位越高,水痕記錄的可靠性越好。自橋梁建成后,橋梁墩臺就無意識的收集了這些資料,并保留在墩臺上。
對于一些水痕不太明顯或交織不清的橋梁墩臺,可以使用模態(tài)模糊識別的方法處理,因不是本文目的,不在此贅述。
我國大部分地區(qū)建有雨量觀測站,氣象部門編輯有“××省地面氣候資料”,據(jù)此,不難確定中小流域的歷史暴雨強度、降雨前干旱天數(shù)等參數(shù)。除氣象部門建立的雨量站外,鐵路沿線還有許多鐵路建立的雨量站,對于橋梁而言這些是寶貴的資源,應(yīng)當(dāng)充分發(fā)揮它們的作用。
中小橋水位與暴雨肯定存在函數(shù)關(guān)系,但這種關(guān)系太復(fù)雜,很難使用普通的數(shù)學(xué)公式直接表示,本文使用人工神經(jīng)網(wǎng)絡(luò)中的 BP(Back Propagation)網(wǎng)絡(luò)[3]表示。
本文取網(wǎng)絡(luò)的拓撲結(jié)構(gòu)為 2-20-1。在這種網(wǎng)絡(luò)結(jié)構(gòu)模型下,取降雨前干旱天數(shù) d和日降雨量a為輸入信號。
式 (4)~(8)中,vji和wjk分別為隱層神經(jīng)元節(jié)點和輸出層神經(jīng)元節(jié)點的權(quán)向量;m取 20;L取 1;f為非線性轉(zhuǎn)移函數(shù),本文取 Sigmoid函數(shù)
學(xué)習(xí)初期取η=0.9。經(jīng)若干次學(xué)習(xí)后,神經(jīng)網(wǎng)絡(luò)得以穩(wěn)定。
獲得穩(wěn)定的神經(jīng)網(wǎng)絡(luò)后,利用預(yù)報暴雨強度,可以回代得到與預(yù)報暴雨強度對應(yīng)的橋梁水位。水位確定后不難確定相應(yīng)流量,并使用相應(yīng)公式求得沖刷深度等指標(biāo)。
圖 2是京廣線金雞灣橋的歷史水位圖,圖中虛線為計算機作模態(tài)識別后留下的中間結(jié)果。表 2是與圖 2對應(yīng)的雨量和水位紀(jì)錄。經(jīng)過 100 000次學(xué)習(xí)后,網(wǎng)絡(luò)達到穩(wěn)定狀態(tài)。金雞灣橋人工神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程見圖 3。
使用“鐵路橋涵水文檢算系統(tǒng)”[5]計算金雞灣橋的情況見圖 5。
表 2 金雞灣橋日最大暴雨與水位關(guān)系表
圖 4 使用“鐵路橋涵水文檢算系統(tǒng)”計算金雞灣橋
占卜橋的原設(shè)計流量 Qp=523m3/s;使用江西省單位線法,得其百年一遇的流量Qp1=521m3/s;使用本文方法,得其百年一遇流量
可見本文方法所得流量與江西省單位線法的結(jié)果比較接近。
滂河橋的原設(shè)計流量 Qp=910m3/s;使用山東省單位線法,得其百年一遇的流量Qp1=910.6m3/s;使用鐵四院法,得其百年一遇的流量 Qp2= 936.6m3/s;使用本文方法得百年一遇流量907.3m3/s;
可見本文方法所得流量與山東省單位線法的結(jié)果比較接近。
①本文提供的小流域流量計算方法,是從雨量直接到橋下水位。不需要確定流域參數(shù),不但減少了工作量,也減少了許多出錯的機會。
②BP網(wǎng)絡(luò)方法是一種成熟的數(shù)學(xué)工具。本文引用這種方法計算的結(jié)果比較穩(wěn)定、可靠。
③這種方法可以應(yīng)用到既有中小橋抗洪能力評估上。設(shè)計新橋需要使用傳統(tǒng)方法計算流量,因為新橋沒有歷史水痕可以利用。但既有中小橋上有水痕可以利用,加上氣象部門記錄的雨量參數(shù),就滿足了本文方法依靠的 2個基礎(chǔ)——雨量和水痕高程。在經(jīng)過幾十年的運用后,既有中小橋都有水痕高程可以利用,有的甚至已經(jīng)遇到了百年一遇的暴雨(如金雞灣橋),顯然利用這一水痕高程是十分可靠的。
④這種方法可以應(yīng)用到既有線中小橋水害預(yù)測上。隨著技術(shù)的進步,如多普勒氣象雷達站的建立、雨量站網(wǎng)的改善等,雨量專業(yè)預(yù)報的精度將會不斷提高。只要獲得未來雨量的預(yù)測值,就可以將它代入神經(jīng)網(wǎng)絡(luò)確定相應(yīng)的水位值,及時預(yù)報橋梁水害,這對于合理安排抗洪的人力、物力是十分有用的。
[1] 肖盛燮,凌天清,陳世民 .公路與橋梁抗洪分析[M].北京:人民交通出版社,1999.
[2] 中華人民共和國鐵道部.鐵路橋梁檢定規(guī)范[M].北京:中國鐵道出版社,2004.
[3] 韓力群 .人工神經(jīng)網(wǎng)絡(luò)理論、設(shè)計及應(yīng)用[M].北京:化學(xué)工業(yè)出版社,2002.
[4] 文雨松,張大偉 .鐵路橋涵可視化水文檢定[M].長沙:中南大學(xué)出版社,2003.
[5] 潘志祥,等 .湖南省大型水電廠流域?qū)I(yè)天氣預(yù)報方法研究[M].北京:氣象出版社,2000.
P457
B
2010-09-10
鄒希云 (1952-),男,高工,主要從事天氣預(yù)報工作。
1003-6598(2010)增刊 -0023-04