孫宇飛
(遼寧省水利水電勘測設計研究院,遼寧 沈陽 110006)
跨河建筑物修建對大遼河河道行洪安全的影響
孫宇飛
(遼寧省水利水電勘測設計研究院,遼寧 沈陽 110006)
大遼河系指太子河入渾河匯流口三岔河起,至河口止的一段河流。該河段水深較大,且受感潮影響,水流形態較為復雜。本文以大遼河鮑家至賞軍河段為例,分析該河段水力特性,對比跨河建筑物修建前后河段水力特性的變化,分析跨河建筑物的修建對該河段行洪安全的影響,為今后類似河段跨河建筑物修建及管理提供經驗借鑒。
大遼河;跨河建筑物;行洪安全
大遼河系指太子河入渾河匯流口三岔河起,至河口止的一段河流。1958年春外遼河堵截后,僅宣泄渾河、太子河洪水。大遼河先后流經海城市、盤山縣、大洼縣、大石橋市、營口市區等市縣,于營口市郊注入渤海。河道全長96km。
大遼河地處遼河下游沖積平原,曲折多彎,而且多急彎和“對頭彎”,歷史上記載大遼河有多處自然裁彎的河段,所以舊河道和河叉較多。整個大遼河的流向是向西南,由三岔河經下坎子、石佛、賞軍、田莊臺、水源,至榮興農場附近,流向折向東南,經鴨島進入營口市。
大遼河系感潮河段,河口系不規則半日周潮,潮流速、潮流量隨潮水向上游推進而減小,同一次漲潮營口站潮流量最大可達5400m3/s、田莊臺站1975m3/s、三岔河站僅462m3/s。
本次分析選擇大遼河鮑家至賞軍河段,該河段位于大遼河三岔河至田莊臺段,具有洪沖枯淤的特性,河道蜿蜒曲折,灘槽分明,是較典型的平原彎曲型河道,造床流量約1000~1200m3/s。該河段既受河水控制又受潮水影響,使彎道凹岸不斷上提下挫,造成塌岸,形成多處險工。河岸土質多為均勻分布的亞粘土和沙壤土。
通過對該河段進行二維水力計算,分析該河段水力特性,在二維水力計算模型中通過增加阻水建筑物的方式分析跨河建筑物修建前后河段水力特性的變化。
2.1 計算范圍
本次水力計算范圍為自上游三岔河至下游大遼河河口的大遼河全河段,計算河長96km。大遼河兩岸均有連續堤防,1999年大遼河防洪整治工程完成后,大遼河防洪標準可達50年一遇,本次計算以堤防為兩岸邊界范圍。
2.2 地形處理
本次計算地形采用2005年1∶10000平面圖和2011年實測橫斷面資料,計算范圍以外地形均做抬高處理,使其在計算中不過流。為反映模型內計算區域的地形特點,對模型中的公路、鐵路等進行了高程點加密,以確保網格剖分后高程的連續。
2.3 邊界條件
本次計算三岔河組合流量沿用《大遼河防洪工程初步設計》成果,典型年選取1960年。模型出口潮位參照《大遼河防洪工程初步設計》成果,選取歷年實測最高潮位3.2m為設計潮位。
2.4 計算參數
(1)計算時段及時間步長
模型采用非恒定流方式計算,計算時長選為153h,時間步長為30s。
(2)計算區域的曼寧值(糙率)
計算區域的曼寧值反映區域內不同地物對水流的阻力作用,采用和地形網格相對應的糙率網格。
本次水力計算參照《大遼河防洪工程初步設計》、《1995年洪水遼河、渾河、大遼河三河水面線驗證》及《大遼河行洪能力分析與防汛特征水位確定》,結合現狀河道的地形、 地貌、河槽組成、水流條件等特性,并考慮近年來河段地貌、植被變化及清障等因素對糙率的影響,最終確定河道糙率。主槽糙率為0.012~0.025,灘地糙率為0.06~0.13。
(3)模型干濕條件
模型干濕度是為了模型在計算二維洪水演進時判斷水流在網格間傳遞的兩個參數,低于“干水深”(dryingdepth)洪水不再演進,高于“干水深”低于“洪水深”(floodingdepth)洪水演進但不參與計算,高于“洪水深”洪水演進且參與計算,此參數越大模型越穩定,根據現場查勘對灘面組成的分析,本次計算“干水深”選用0.005m,“洪水深”選用0.05m。
2.5 計算方案
工程所在河段防洪標準為50年一遇,本次對河段50年及100年一遇洪水頻率進行計算,以分析相應洪水頻率下該河段水力特性。水力計算成果見圖1~圖4。

圖1 天然情況河段5O年一遇水位計算成果圖

圖2 天然情況河段1OO年一遇水位計算成果圖

圖3 天然情況河段5O年一遇流速計算成果圖

圖4 天然情況河段1OO年一遇流速計算成果圖
本河段屬蜿蜒型河道,灘地占河道行洪寬度較大,考慮本河段河道特性及跨河建筑物工程建設特點,本次阻水建筑物主要添加于河道兩側灘地之中,分別于計算河段右岸灘地及其下游左岸灘地添加阻水建筑物,右岸灘地添加阻水建筑物6座,左岸灘地添加阻水建筑物2座,阻水建筑物概化為矩形,尺寸為10m×10m,將阻水建筑物概化后添加于網格地形中使其整體不過流,概化后河段模型地形見圖5。

圖5 阻水建筑物概化后模型地形圖(紅色節點為添加阻水建筑物位置)
添加阻水建筑物后對該河段進行二維水力計算,并與相同水力條件下天然河道水力計算成果進行對比,成果見圖6~圖9。

圖6 工程修建前后5O年一遇水位差值對比圖

圖7 工程修建前后1OO年一遇水位差值對比圖

圖8 工程修建前后5O年一遇流速差值對比圖

圖9 工程修建前后1OO年一遇流速差值對比圖
工程修建前后河段水力特性差值對比分析見表1,由分析成果可見,新建阻水建筑物侵占河道行洪面積,會導致建筑物迎水側水位壅高,水位壅高值隨著距建筑物的距離的增加而逐漸減小,直至與現狀情況一致。建筑物背水側水位會降低,水位降低值隨著距建筑物距離的增加降低而逐漸減小。建筑物迎、背水側由于建筑物阻水流速會降低,相應的建筑物兩側流速會增大,流速變化值隨著距建筑物距離增加而逐漸減小。
50年一遇洪水流量下,左岸灘地阻水建筑物位置天然情況下河道水位為4.75m,工程修建后水位為4.75m,右岸灘地阻水建筑物位置天然情況下河道水位為 5.65m,工程修建后河道水位為5.65m,壅水實際影響值為mm級,阻水建筑物最大壅水高度均小于0.01m,壅水最大影響范圍為150m,阻水建筑物的建設對河道水位影響不大,不會影響河道行洪。
100年一遇洪水流量下,左岸灘地阻水建筑物位置天然情況下河道左灘水位為5.86m,工程修建后水位為5.86m,右岸灘地阻水建筑物位置天然情況下河道水位為 6.74m,工程修建后水位為6.74m,壅水實際影響值仍為mm級,阻水建筑物最大壅水高度均小于0.01m,壅水最大影響范圍為200m,工程的修建對河道水位影響不大,壅水影響范圍同樣較小,不會影響河道行洪及堤防防洪安全。
阻水建筑物工程修建前后,河道流速大小變化不大,流速分布基本不變。建筑物附近由于工程阻擋,水流形態較天然情況復雜紊亂,建筑物迎流側和背流側由于阻水作用流速減小,建筑物兩側增加過流流速增大,但流速變化值均不大。50年一遇洪水條件下,左岸灘地阻水建筑物位置附近主槽最大流速2.2m/s,灘地最大流速0.61m/s,建筑物附近流速變化最大值0.08m/s。右岸灘地阻水建筑物位置附近主槽最大流速1.03m/s,灘地最大流速0.51m/s,建筑物附近流速變化最大值0.06m/s。垂直于流線方向流速最大影響范圍為200m。
100年一遇洪水流量下,左岸灘地阻水建筑物位置附近主槽最大流速2.45m/s,灘地最大流速0.73m/s,建筑物附近流速變化最大值0.13m/s。右岸灘地阻水建筑物位置附近主槽最大流速1.75m/s,灘地最大流速0.62m/s,建筑物附近流速變化最大值0.11m/s。垂直于流線方向流速最大影響范圍為350m。

表1 工程修建前后計算成果對比表
經以上分析可知,由于大遼河為沖積平原河流,河道比降較緩且灘地較寬,只要阻水建筑物修建于灘地之上,工程占河道行洪斷面較小且并未擠占河道主槽,阻水建筑物對河道水位及流速影響均相對較小,對該河段行洪能力影響較小,不會對河段行洪安全造成明顯不利影響。大遼河為蜿蜒型河段,本文所選河段河道水力特性與大遼河其他河段特性均較為相似,阻水建筑物對河道水力特性的影響對我國北方其他蜿蜒型中下游河段有相對普遍的參考價值。然而,不同跨河建筑物工程型式及特點不盡相同,模型中對相應跨河建筑物概化方式及計算方法仍有待探討,需進一步予以完善。
[1]李煒,徐孝平.水力學[M].武漢:武漢水利電力大學出版社,1999.
[2]謝鑒衡.河床演變及整治[M].北京:中國水利水電出版社,1997.
[3]宋廣堯.水力學與橋涵水文[M].北京:中國鐵道出版社,1987.
[4]張學齡.橋涵水文(第二版)[M].北京:人民交通出版社,1986.
[5]俞高明.橋涵水力水文[M].北京:人民交通出版社,2002.
[6]太子河(含大遼河)防洪治理工程可行性研究報告[R].遼寧省水利水電勘測設計研究院,2015.
[7]盤錦220千伏鶴鄉變-遼河變電力線路工程跨越大遼河防洪評價報告[R].遼寧省水利水電勘測設計研究院,2015.
[8]楊揚,韓鳳霞,劉軍梅,等.山區河道水面線和管涵漫水橋壅水計算探討[J].水利規劃與設計,2013(01):50-53.
[9]熊繼斌,李江波,石昊,等.三維數學模型在跨河工程防洪評價中的應用[J].水利規劃與設計,2013(10):26-30.
[10]丁立國,雷炎,任聃,等.遼蒙界老哈河橋壅水影響分析[J].水利建設與管理,2007(08):21-23.
[11]駱文廣,陸晶,賴昊,等.近距多座并行橋梁墩前壅水研究[J].四川大學學報(工程科學版),2015(04):9-12.
TV131.4
B
1672-2469(2017)02-0090-04
10.3969/j.issn.1672-2469.2017.02.029
2016-05-19
孫宇飛(1986年—),男,工程師。