彭興 馬祺瑞 果鵬
摘要:為研究雅魯藏布江流域江熱橋?qū)拥婪篮榈挠绊?,基于MIKE21軟件建立該研究區(qū)域的二維水動(dòng)力數(shù)值模型,并對(duì)該工程的行洪水位和流場(chǎng)影響進(jìn)行了模擬計(jì)算分析。結(jié)果表明:在10 a一遇和50 a一遇洪水條件下,工程建設(shè)引起水位最大壅高值分別為7.8 cm和9.8 cm,引起流速增加的最大值分別為0.11 m/s和0.12 m/s,工程建設(shè)不會(huì)對(duì)河道防洪和河勢(shì)造成明顯的不利影響。此外,對(duì)橋墩沖刷、岸坡沖刷、岸坡邊坡穩(wěn)定進(jìn)行了計(jì)算分析,計(jì)算結(jié)果滿足規(guī)范要求。研究結(jié)果可為江熱橋工程項(xiàng)目提供防洪影響評(píng)價(jià)方面的技術(shù)支撐。
關(guān)鍵詞:防洪影響評(píng)價(jià); 二維水動(dòng)力模型; 橋梁工程; 江熱橋; 雅魯藏布江
中圖法分類號(hào):TV214
文獻(xiàn)標(biāo)志碼:A
DOI:10.15974/j.cnki.slsdkb.2023.S2.002
文章編號(hào):1006-0081(2023)S2-0005-06
0 引 言
橋梁工程跨越河流架設(shè)于兩岸之間,在施工及運(yùn)行過程中,往往會(huì)導(dǎo)致河流形態(tài)變化,破壞河道局部沖淤平衡,有的甚至?xí)?duì)區(qū)域河段行洪及河勢(shì)產(chǎn)生較大影響。因此,跨河橋梁工程的防洪影響評(píng)價(jià)歷來是水行政主管部門關(guān)注的重點(diǎn)對(duì)象,同時(shí),做好跨河橋梁工程的防洪影響分析及評(píng)價(jià)是確保區(qū)域防洪安全的重要保障。
本文以西藏自治區(qū)山南市浪卡子縣江熱橋?yàn)槔?,從多方面開展分析論證,將橋梁對(duì)流域局部壅水、流速變化、岸坡穩(wěn)定的影響進(jìn)行了系統(tǒng)分析,可為相似項(xiàng)目研究提供參考,同時(shí)也可為實(shí)際橋梁設(shè)計(jì)建設(shè)及工程附近區(qū)域防洪減災(zāi)提供參考。
1 工程概況
山南市浪卡子縣江熱橋位于雅魯藏布江干流中游索朗嘎咕峽谷段,左、右岸行政區(qū)劃分屬拉薩市尼木縣、山南市浪卡子縣,工程位置見圖1。
江熱橋按照四級(jí)公路標(biāo)準(zhǔn)設(shè)計(jì),橋面寬度為7.5 m(行車道)+2×1.25 m(人行道、護(hù)欄),共10 m,設(shè)計(jì)車速為30 km/h,橋型布置見圖2。工程從左岸往右岸布置了0號(hào)橋臺(tái)、1號(hào)橋墩、2號(hào)橋墩和3號(hào)橋臺(tái),其中1號(hào)橋墩、2號(hào)橋墩位于河道內(nèi)。
2 洪水影響分析計(jì)算
2.1 二維水動(dòng)力模型
采用丹麥水力研究所(DHI)開發(fā)的MIKE21 FM建立工程區(qū)域河道平面二維水流數(shù)學(xué)模型,計(jì)算分析10 a一遇設(shè)計(jì)洪水和50 a一遇設(shè)計(jì)洪水條件下,擬建工程對(duì)工程區(qū)域河段防洪設(shè)計(jì)洪水位和流場(chǎng)的影響。二維水流模型基本方程如下[1-3]:
水流連續(xù)方程:
水流運(yùn)動(dòng)方程:
式中:z為水位;p、q為x、y方向單寬流量;h為水深;s為源、匯項(xiàng);six、siy為源、匯項(xiàng)在x、y方向的分量;c為謝才阻力系數(shù);Ω為科氏力;E為渦黏擴(kuò)散系數(shù),根據(jù)Smagorinsky公式確定。
式中:u、v為x、y方向垂線平均流速;
Δ為網(wǎng)格間距;CS為計(jì)算參數(shù),0.25<CS<1.0。
2.2 邊界條件參數(shù)取值
2.2.1 研究范圍及網(wǎng)格布置
模型計(jì)算范圍為橋址以上2.0 km至橋址以下2.2 km、全長(zhǎng)4.2 km的河段,見圖3。采用無結(jié)構(gòu)三角形混合網(wǎng)格剖分計(jì)算區(qū)域,網(wǎng)格單元數(shù)為45 669個(gè),特別對(duì)工程局部區(qū)域進(jìn)行了加密,工程附近網(wǎng)格尺度為2 m。
2.2.2 計(jì)算工況與邊界條件
選取以下2種計(jì)算工況開展計(jì)算分析。
工況1:橋址10 a一遇設(shè)計(jì)洪水。相應(yīng)洪水流量為4 680 m3/s、設(shè)計(jì)洪水位為3 642.21 m,分別作為工況1的模型上下邊界條件。
工況2:橋址50 a一遇設(shè)計(jì)洪水。相應(yīng)洪水流量為6 220 m3/s、設(shè)計(jì)洪水位為3 644.01 m,以此作為工況2的模型上下邊界條件。
在計(jì)算過程中,網(wǎng)格區(qū)域內(nèi)的部分單元在漲水時(shí)被“淹沒”,在落水時(shí)則“干出”。為正確反映這些單元的干濕變化,模型中采用了如下的動(dòng)邊界模擬技術(shù):選定一臨界水深(hmin一般取為0.01 m),當(dāng)某時(shí)刻某濕單元的實(shí)際水深(水位減去河底高程)小于hmin時(shí),認(rèn)為該單元“干出”,在本步時(shí)末退出模型計(jì)算;干單元的水位值由附近濕單元的水位插值得到,當(dāng)某時(shí)刻某干單元實(shí)際水深大于臨界水深且有流入的動(dòng)量通量時(shí),則讓其變濕,并恢復(fù)該單元的模型計(jì)算。
對(duì)于橋墩阻水處理,則通過計(jì)算現(xiàn)有橋墩柱體上減少的拖曳力來反映橋墩對(duì)水流的阻力,而每個(gè)柱體上減少的拖曳力在MIKE21 FM中通過動(dòng)量公式轉(zhuǎn)化為切應(yīng)力來計(jì)算。
2.2.3 模型率定與驗(yàn)證
由于缺少實(shí)測(cè)資料,經(jīng)多次試算來確定河段的糙率系數(shù)。計(jì)算結(jié)果顯示:當(dāng)河段綜合糙率系數(shù)為0.048,進(jìn)口流量為4 680 m3/s時(shí),橋址計(jì)算水位為3 652.521 m;再采用相同的糙率系數(shù),進(jìn)口流量為6 220 m3/s(50 a一遇)時(shí),橋址計(jì)算水位為3 654.959 m。由此可見,橋址計(jì)算水位與設(shè)計(jì)洪水位的絕對(duì)誤差均小于0.1 cm,由此確定河段的綜合糙率為0.048。
2.3 10 a一遇設(shè)計(jì)洪水計(jì)算結(jié)果分析
2.3.1 阻水面積計(jì)算
根據(jù)橋梁相關(guān)設(shè)計(jì)圖紙,將加入橋梁后的斷面與原始斷面進(jìn)行對(duì)比,計(jì)算大橋建成前后斷面過水面積與橋墩水下面積,得到橋墩阻水面積占河道過水?dāng)嗝婵偯娣e的比例。結(jié)果顯示:在10 a一遇設(shè)計(jì)洪水位工況下,原河道斷面過水面積1 377.15 m2,橋墩占用面積32.11 m2(表1),橋墩阻水面積占河道過水?dāng)嗝婵偯娣e的2.33%。
2.3.2 工程對(duì)水位的影響分析
在江熱橋上游與該橋距離為5 m的橫斷面內(nèi),等間距布置6個(gè)評(píng)價(jià)點(diǎn),采用前述構(gòu)建的二維模型分析擬建工程修建前后的壅水影響[2]。不同評(píng)價(jià)點(diǎn)工修建程前后水位變化統(tǒng)計(jì)結(jié)果見表2,可以看出,工程后水位最大壅高為7.8 cm。
工程修建前后橋梁附近局部壅水見圖4,可以看出,水位變幅0.5 cm以上的影響范圍為橋址上游62.4 m至橋址下游21.9 m;水位變幅0.1 cm以上的影響范圍為橋址上游86.1 m至橋址下游 26.5 m。
2.3.3 工程對(duì)流速的影響分析
10 a一遇洪水條件下,不同評(píng)價(jià)點(diǎn)工程修建前后流速變化統(tǒng)計(jì)結(jié)果見表3,可以看出,工程修建后流速基本沒有變化,最大流速變化0.11 m/s。
10 a一遇洪水條件下,工程修建前后流速等值線分布和流場(chǎng)分布見圖5和圖6,可以看出,工程修建前后流速大小和流場(chǎng)變化不大。
2.4 50 a一遇設(shè)計(jì)洪水計(jì)算結(jié)果分析
2.4.1 阻水面積計(jì)算
根據(jù)橋梁相關(guān)設(shè)計(jì)圖量算,在橋址50 a一遇設(shè)計(jì)洪水位(3 654.96 m)工況下,原河道斷面過水面積1 587.75 m2,橋墩占用面積40.41 m2,橋墩阻水面積占河道過水?dāng)嗝婵偯娣e的比例為2.50%。
2.4.2 工程對(duì)水位的影響分析
不同評(píng)價(jià)點(diǎn)工程修建前后水位變化統(tǒng)計(jì)結(jié)果見表4,可以看出,工程修建后水位最大壅高為98 cm。
工程修建后橋梁附近局部壅水見圖7,可以看出,水位變幅0.5 cm以上的影響范圍為橋址上游 71.2 m至橋址下游17.6 m;水位變幅0.1 cm以上的影響范圍為橋址上游103.2 m至橋址下游 37.6 m。
2.4.3 工程對(duì)流速的影響分析
50 a一遇洪水條件下,不同評(píng)價(jià)點(diǎn)工程修建前后流速變化統(tǒng)計(jì)結(jié)果見表5,可以看出,工程修建后流速基本沒有變化,最大流速變化0.12 m/s。
50 a一遇洪水條件下,工程修建前后流速等值線分布和流場(chǎng)分布見圖8和圖9,可以看出,工程修建前后流速大小和流場(chǎng)變化不大。
3 橋墩沖刷分析計(jì)算
根據(jù)JTGC 30-2015《公路工程水文勘測(cè)設(shè)計(jì)規(guī)范》,建橋后橋梁墩臺(tái)附近的河床沖刷由自然演變沖刷、一般沖刷和局部沖刷組成[4-7]。根據(jù)河道演變分析成果,工程河段形勢(shì)總體穩(wěn)定,河床沖淤變化較小。本次計(jì)算主要分析一般沖刷和局部沖刷兩種沖刷形式。根據(jù)地質(zhì)鉆孔揭示情況,按照非黏性土河槽一般沖刷和非黏性土橋墩局部沖刷計(jì)算。
3.1 10 a一遇設(shè)計(jì)洪水橋墩沖刷分析
10 a一遇設(shè)計(jì)洪水條件下橋址天然設(shè)計(jì)洪水位為3 652.52 m。在此水位下,擬建工程后橋址所在位置的斷面平均流速為3.48 m/s,1號(hào)與2號(hào)橋墩受到水流沖刷影響,0號(hào)與3號(hào)橋臺(tái)沒有受到水流淹沒,因此開展1號(hào)與2號(hào)橋墩的沖刷分析計(jì)算。
10 a一遇設(shè)計(jì)洪水條件下,一般沖刷后的最大水深約為27.46 m,工程修建后橋址所在斷面河槽部分一般沖刷深度為1.92 m。
橋墩局部沖刷受水流流速、水深、橋墩結(jié)構(gòu)及土質(zhì)因素的影響。橋墩結(jié)構(gòu)決定墩形系數(shù)及阻水面積。擬建橋梁橋墩均為直徑統(tǒng)一的圓柱體結(jié)構(gòu),故其墩形系數(shù)取1.00。橋址因缺少床沙級(jí)配資料,故選用河段下游羊村水文站實(shí)測(cè)床沙級(jí)配近似代替,取其中值粒徑38 mm。經(jīng)計(jì)算,1號(hào)和2號(hào)橋墩的局部沖刷深度分別為5.02 m和7.22 m。由一般沖刷和局部沖刷計(jì)算結(jié)果可知,1號(hào)和2號(hào)橋墩的最大沖深分別為6.94 m和9.14 m。1號(hào)和2號(hào)橋墩分別埋入河底以下18.72 m和13.77 m,即最大沖深以下11.78 m和4.63 m,滿足橋梁基底埋深安全值要求。
3.2 50 a一遇設(shè)計(jì)洪水橋墩沖刷分析
經(jīng)計(jì)算,50 a一遇設(shè)計(jì)洪水條件下,工程修建后橋址所在斷面河槽部分一般沖刷后的最大水深約為30.04 m,一般沖刷深度為2.08 m;50 a一遇設(shè)計(jì)洪水條件下1號(hào)和2號(hào)橋墩的局部沖刷深度分別為6.14 m和8.83 m,最大沖深分別為8.22 m、10.91 m。1號(hào)和2號(hào)橋墩分別埋入河底以下18.72 m和13.77 m,即最大沖深以下10.50 m和2.86 m,滿足JTGC 30-2015《公路工程水文勘測(cè)設(shè)計(jì)規(guī)范》橋梁基底埋深安全值要求。
4 岸坡沖刷分析計(jì)算
橋梁施工完成后,橋墩附近水位將略微壅高,流速將有所增加,引起左岸岸坡沖刷加大[8-9],需分析工程前后左岸岸坡坡腳處的沖刷變化情況。岸坡沖刷深度分析采用GB 50286-2013《堤防工程設(shè)計(jì)規(guī)范》附錄D.2.2中所列公式予以計(jì)算。
工程前后左岸岸坡坡腳處沖刷深度計(jì)算結(jié)果見表6。
由表6可知,在10 a一遇設(shè)計(jì)洪水位情況下,江熱橋建設(shè)后左岸岸坡坡腳沖刷深度由0.80 m增加至0.86 m;在50 a一遇設(shè)計(jì)洪水位情況下,江熱
橋建設(shè)后左岸岸坡坡腳沖刷深度由1.27 m增加至1.33 m;江熱橋建設(shè)對(duì)左岸岸坡沖刷無明顯不利影響。
5 工程河段岸坡穩(wěn)定計(jì)算分析
5.1 正常運(yùn)用條件
工況1:橋梁施工前現(xiàn)狀岸坡,河道水位為10 a一遇洪水位3 652.52 m。工況2:橋梁施工前現(xiàn)狀岸坡,河道水位為50 a一遇洪水位3 654.96 m。
工況3:橋梁施工后岸坡,河道水位為10 a一遇洪水位3 652.52 m。工況4:橋梁施工后岸坡,河道水位為50 a一遇洪水位3 654.96 m。
5.2 非常運(yùn)用條件
根據(jù)GB 18306-2015《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖(1∶400萬)》,江熱橋工程50 a超越概率10%地震動(dòng)峰值加速度0.20g,對(duì)應(yīng)地震基本烈度Ⅷ度,非常運(yùn)用條件考慮以下工況。
工況5:橋梁施工前現(xiàn)狀岸坡,河道水位為10 a 一遇洪水位3 652.52 m,遭遇Ⅷ度地震。
工況6:橋梁施工前現(xiàn)狀岸坡,河道水位為50 a一遇洪水位3 654.96 m,遭遇Ⅷ度地震。
工況7:橋梁施工后岸坡,河道水位為10 a一遇洪水位3 652.52 m,遭遇Ⅷ度地震。
工況8:橋梁施工后岸坡,河道水位為50 a一遇洪水位3 654.96 m,遭遇Ⅷ度地震。
荷載按公路Ⅰ級(jí)荷載考慮,作用于G318路面。
工程河段位于高山峽谷區(qū),河道防洪標(biāo)準(zhǔn)為10 a一遇,邊坡為5級(jí)。根據(jù)SL 386-2007《水利水電工程邊坡設(shè)計(jì)規(guī)范》,5級(jí)邊坡在正常運(yùn)用條件下的抗滑穩(wěn)定安全系數(shù)可取1.10~1.05,在遭遇地震的非常運(yùn)用條件下的抗滑穩(wěn)定安全系數(shù)可取 1.05~1.00。本次計(jì)算偏安全考慮,抗滑穩(wěn)定最小安全系數(shù)取規(guī)定范圍內(nèi)的大值。
抗滑穩(wěn)定采用理正邊坡穩(wěn)定分析軟件進(jìn)行計(jì)算,該程序采用最優(yōu)化原理,搜索最小安全系數(shù)和相應(yīng)臨界滑面。本次計(jì)算方法采用簡(jiǎn)化畢肖普法[10-11],不同工況下抗滑穩(wěn)定安全系數(shù)計(jì)算結(jié)果見表7。計(jì)算結(jié)果表明,各工況岸坡抗滑穩(wěn)定安全系數(shù)大于規(guī)范要求允許值,滿足規(guī)范要求,工程實(shí)施對(duì)左岸岸坡抗滑穩(wěn)定無不利影響。但0號(hào)橋臺(tái)開挖應(yīng)注意采取插打鋼板樁等措施做好防護(hù),確保施工期G318國(guó)道通行安全。
6 結(jié) 論
基于工程水文地質(zhì)、設(shè)計(jì)圖紙等基本資料,對(duì)江熱橋工程臨近河道的洪水影響、橋墩沖刷、岸坡沖刷以及岸坡穩(wěn)定進(jìn)行了分析計(jì)算,綜合評(píng)價(jià)了擬建項(xiàng)目對(duì)防洪的影響,得到的主要結(jié)論如下。
(1)工程位于高山峽谷河段,兩岸無防洪工程,無防洪保護(hù)對(duì)象,工程位置10 a一遇洪水位和50 a一遇洪水位疊加工程建設(shè)產(chǎn)生的最大壅高值后,仍然低于全橋底板最低高程(3 658.955 m)和G318國(guó)道高程(3 661.80 m),工程建設(shè)不會(huì)對(duì)河道防洪造成不利影響。
(2)在10 a一遇和50 a一遇洪水條件下,工程的興建引起流速增加的最大值分別為0.11 m/s和0.12 m/s,工程對(duì)流速的影響局限于工程附近,對(duì)計(jì)算河段整體流場(chǎng)影響較小,不會(huì)引起較大范圍的河床沖淤變化。工程建設(shè)不會(huì)對(duì)河勢(shì)造成明顯的不利影響。
(3)10 a一遇洪水情況下,工程修建后1號(hào)和2號(hào)橋墩最大沖深分別為6.94 m和9.14 m;50 a一遇洪水情況下,1號(hào)和2號(hào)橋墩最大沖刷深度分別為8.22 m和10.91 m。橋墩最大沖深小于橋墩埋深,滿足JTGC 30-2015《公路工程水文勘測(cè)設(shè)計(jì)規(guī)范》橋梁基底埋深安全值要求。
(4)在10 a一遇洪水情況下,工程修建后橋址所在左岸岸坡坡腳沖刷深度由0.80 m增加至0.86 m;在50 a一遇洪水情況下,左岸岸坡坡腳沖刷深度由1.27 m增加至1.33 m;增加幅度為0.06 m,因此工程建設(shè)對(duì)左岸岸坡沖刷無明顯不利影響。
(5)在10 a一遇設(shè)計(jì)洪水和50 a一遇設(shè)計(jì)洪水位情況下,正常運(yùn)用條件工程前后岸坡抗滑穩(wěn)定安全系數(shù)為1.31~1.49,非常運(yùn)用條件下工程前后岸坡抗滑穩(wěn)定安全系數(shù)為1.15~1.35,各工況岸坡抗滑穩(wěn)定安全系數(shù)大于等于規(guī)范要求最大值,滿足規(guī)范要求,工程建設(shè)對(duì)左岸岸坡穩(wěn)定無不利影響。
參考文獻(xiàn):
[1] 許婷.丹麥MIKE21模型概述及應(yīng)用實(shí)例[J].水利科技與經(jīng)濟(jì),2010,16(8):867-869.
[2] 袁雄燕,徐德龍.丹麥MIKE21模型在橋渡壅水計(jì)算中的應(yīng)用研究[J].人民長(zhǎng)江,2006(4):31-32,52.
[3] 郭鳳清,屈寒飛,曾輝,等.基于MIKE21 FM模型的蓄洪區(qū)洪水演進(jìn)數(shù)值模擬[J].水電能源科學(xué),2013,31(5):34-37.
[4] 張胡,閆杰超,陳凱華.典型橋墩局部沖刷及防護(hù)特性數(shù)值模擬研究[J].水利水電快報(bào),2018,39(4):45-47,52.
[5] 張佰戰(zhàn),李付軍.橋墩局部沖刷計(jì)算研究[J].中國(guó)鐵道科學(xué),2004(2):49-52.
[6] 李奇,王義剛,謝銳才.橋墩局部沖刷公式研究進(jìn)展[J].水利水電科技進(jìn)展,2009,29(2):85-88,94.
[7] 王玲玲,俞佳亮,朱海,等.橋墩局部沖刷研究進(jìn)展[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2019,47(3):223-229.
[8] 葉威,章光,胡少華,等.河流沖刷作用下堤岸穩(wěn)定性演化分析[J].水利水電快報(bào),2019,40(1):42-47.
[9] 謝立全,于玉貞.滲流作用下的岸坡泥沙起動(dòng)條件[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2006(9):1534-1537.
[10] 楊曉松,高濤.瑞典條分法與簡(jiǎn)化畢肖普法在海堤抗滑穩(wěn)定分析中的比較[J].水運(yùn)工程,2017(2):27-32.
[11] 王新奇,孫勝利,趙洪嶺.邊坡穩(wěn)定分析方法評(píng)述與應(yīng)用[J].水利技術(shù)監(jiān)督,2010,18(1):7-10.