鄧新超

(深圳市水務(wù)規(guī)劃設(shè)計(jì)院股份有限公司，廣東 深圳 518022 )

隨著社會經(jīng)濟(jì)不斷發(fā)展，高速公路等建設(shè)項(xiàng)目逐漸增多，規(guī)模也逐漸增大，特別是在河道擬建橋墩后，使河道的行洪斷面、水流形態(tài)、沖淤平衡等發(fā)生一系列變化，進(jìn)而影響河道排放口排水，甚至還有可能影響河勢穩(wěn)定性。因此，需要通過Mike21FM數(shù)學(xué)模型結(jié)合SMS(Surface-Water Modeling System)軟件對工程區(qū)域進(jìn)行模擬計(jì)算，求出工程壅水高度及長度、流速流場變化及設(shè)計(jì)水位，形成三維可視化效果，為防洪評價(jià)提供可靠依據(jù)，從而避免因工程建設(shè)給河道防洪安全造成影響。

1 項(xiàng)目概況

機(jī)荷高速公路工程屬于G15沈海高速公路深圳段，位于深圳市中心區(qū)域，是粵港澳大灣區(qū)東西向的交通中軸線，也是深圳境內(nèi)通往機(jī)場的唯一一條東西走向的高速公路。起點(diǎn)位于龍崗區(qū)的荷坳立交，順接惠鹽高速；終點(diǎn)位于寶安區(qū)的鶴州立交，連接廣深高速，與深中通道順接，路線長約42km；該工程在道路樁號40+610—40+820段采用2跨型式通過九圍河，其中2座橋墩設(shè)置在河道行洪斷面內(nèi)，因此，需要對河道內(nèi)設(shè)置橋墩進(jìn)行數(shù)學(xué)模擬計(jì)算。

2 模型建立及計(jì)算

2.1 模型計(jì)算方法

計(jì)算采用丹麥丹華水利環(huán)境技術(shù)有限公司(DHI)的Mike21FM軟件，屬于Mike系列軟件中的一種，Mike系列軟件廣泛地應(yīng)用于潮汐、水流、水質(zhì)、熱流通、湖震、風(fēng)暴潮、密度流、波浪、泥沙侵蝕、船運(yùn)、防波堤布置等二、三維水力學(xué)現(xiàn)象的研究，Mike系列也因此被推薦為河流、河口、湖泊和海洋近岸水流的仿真模擬軟件。

Mike21FM是二維水動力模型，該模型水動力的控制方程為納維-斯托克斯(Navier-Stokes)方程，模型采用有限體積法將水動力控制的偏微分方程進(jìn)行離散并求解，模型可以采用不規(guī)則網(wǎng)格對地形進(jìn)行貼合，提升計(jì)算精度。

2.2 模型范圍及網(wǎng)格劃分

2.2.1模型范圍

本次模型范圍初步界定在擬建橋梁上下游河道順直、河勢較為穩(wěn)定的河段(樁號0+411.87—0+611.13)，首先將CAD地形數(shù)據(jù)轉(zhuǎn)化為Mike21FM模型可識別的地形數(shù)據(jù)，構(gòu)建模型計(jì)算區(qū)域(如圖1所示)，并通過網(wǎng)格將地形數(shù)據(jù)離散化求解。本次橋梁的壅水分析計(jì)算考慮了河道內(nèi)機(jī)荷高速現(xiàn)狀橋墩阻水作用，在設(shè)計(jì)標(biāo)準(zhǔn)洪水下通過網(wǎng)格將地形數(shù)據(jù)離散化進(jìn)而求解1#和2#橋梁建設(shè)前后的水面線高程差，得出1#和2#橋梁的壅水情況。

圖1 壅水分析計(jì)算模型范圍及主要計(jì)算因素(單位：m)

2.2.2網(wǎng)格劃分

在高程數(shù)據(jù)質(zhì)量及精度有限的情況下，生成的網(wǎng)格往往會出現(xiàn)局部異常，主要原因是地形數(shù)據(jù)缺失或高程異常所導(dǎo)致(如圖2所示)，本次計(jì)算模型區(qū)域采用SMS軟件對模擬區(qū)域地形進(jìn)行處理，橋墩通過拔高概化，最終形成非結(jié)構(gòu)三角網(wǎng)格，對2座擬建橋墩及5座現(xiàn)狀橋墩等主要構(gòu)筑物在模型中的網(wǎng)格進(jìn)行局部加密，如圖3所示。工程前后模型地形網(wǎng)格劃分如圖4所示。通過2種方法比較，聯(lián)合使用SMS軟件處理后生成的網(wǎng)格更加平滑、美觀，使得Mike21FM模擬計(jì)算結(jié)果更為精確、可靠。

圖2 未使用SMS軟件處理的地形三維效果圖

圖3 使用SMS軟件處理后的地形三維效果圖

圖4 工程前后模型地形Mike21 FM網(wǎng)格圖(單位：m)

2.3 計(jì)算工況

本次模擬計(jì)算在河道遭遇50a一遇設(shè)計(jì)標(biāo)準(zhǔn)洪水時(shí)建橋前后的壅水情況，模型設(shè)置2個(gè)開邊界，將九圍河上游設(shè)置為流量邊界，九圍河下游設(shè)置為水位邊界。根據(jù)設(shè)計(jì)洪水復(fù)核成果，九圍河上邊界50a一遇設(shè)計(jì)流量Q=84.32m3/s，同頻率下的九圍河下游設(shè)計(jì)洪水位H=28.13mm。

2.4 糙率選取

糙率(n值)是反映河床、岸壁形狀的不規(guī)則性和表面粗糙程度的一個(gè)系數(shù)，在水流運(yùn)動過程中，直接影響沿程能量損失的大小。根據(jù)河道特性、河道的地形地貌資料及岸壁情況綜合確定河道糙率。

河道糙率選擇主要參照《水力計(jì)算手冊》，綜合考慮河道斷面現(xiàn)狀及規(guī)劃情況，根據(jù)不同的材料選取不同的糙率。模型計(jì)算范圍河道深淺不一，邊灘和河床糙率各不相同，故模型在設(shè)定糙率時(shí)采用隨水深的變化值，根據(jù)每個(gè)網(wǎng)格的水深值設(shè)定不同的糙率值。根據(jù)相關(guān)研究，變動糙率的取值為：

n=2n1(h≤1.0m)

(1)

n=n1+n2/h(h≤1.0m)

(2)

式中，h—水深，m，取正值;n1—河床沙粒糙率系數(shù)，其取值范圍為0.01～0.025；n2—河床形態(tài)糙率系數(shù)，其取值范圍為0.005～0.01；n1和n2合稱糙率系數(shù)。

根據(jù)模擬計(jì)算涉河處河段50a一遇工況下水深成果，如圖5所示。該河道內(nèi)水深均大于1m，最深處水深達(dá)到5m，因此糙率計(jì)算方式采用式(2)取值，最后根據(jù)河道構(gòu)筑物特征及取值公式取綜合糙率取值為0.028。

圖5 模擬河道水深圖(單位：m)

2.5 行洪影響綜合分析

2.5.1過流能力影響分析

主要分析橋梁修建后對九圍河行洪區(qū)域整體過流能力影響。根據(jù)模型分析計(jì)算可知，工程后水位較工程前有所上升，但河道水位上漲高度微小，水位不漫過左右岸堤頂高程，不對河道過流造成影響，詳情見表1，如圖6所示。

2.5.2流速變化影響分析

根據(jù)模型計(jì)算擬建橋墩局部流速圖，所在河段因橋墩修建后斷面縮窄，造成橋墩與流水方向兩側(cè)流速增大，橋墩下游因橋墩阻水效應(yīng)而致流速減小，其中1#橋墩兩側(cè)流速為1.25m/s，2#橋墩兩側(cè)流速為1.05m/s，詳情如圖7所示，見表2。

表1 擬建1#和2#橋建橋前后九圍河水位變化情況(P=2%) 單位：m

表2 沖刷深度計(jì)算成果

由于本次擬建橋墩后將在橋墩處采用拋石進(jìn)行防護(hù)，表面采用砌石整平，防護(hù)措施能滿足橋墩處最大1.25m/s流速，因此，工程可滿足抗沖刷要求。

2.5.3壅水影響分析

通過模型分析計(jì)算工程前后水面線，可以得出橋墩的壅水高度及長度。根據(jù)上述計(jì)算成果可知，1#橋墩工程前水位為28.148m，工程后水位為28.202m；2#橋墩工程前水位為28.131m，工程后水位為28.139m，見表3。擬建1#、2#橋墩處壅水流場圖如圖8所示。由圖8可知，擬建1#、2#橋墩后產(chǎn)生的橋墩處最大壅水值為0.054m，橋墩擬建所造成的水位壅高對河道行洪影響甚微，且河道沿岸排放口不造成影響。

表3 擬建1#和2#橋墩建橋前后九圍河水位變化情況

圖6 工程后1#、2#橋墩處水位與河岸關(guān)系三維圖

圖7 擬建1#、2#橋墩局部流速圖

圖8 擬建1#、2#橋墩處壅水流場圖

3 結(jié)論

由以上對模型地形數(shù)據(jù)的處理、分析和計(jì)算可得出以下結(jié)論。通過SMS(Surface-Water Modeling System)軟件刪除地形數(shù)據(jù)中的異常高程點(diǎn)，并將擬建橋墩拔高網(wǎng)格概化。經(jīng)過處理后的模型網(wǎng)格精度更高，輪廓與實(shí)際更為貼近，以此可避免因網(wǎng)格地形凸變給模擬運(yùn)算造成影響。運(yùn)用Mike21FM模擬計(jì)算可以全面、準(zhǔn)確、直觀地反映橋墩處的流場變化、壅水高度及范圍變化等特性。根據(jù)Mike軟件自帶的Mike Animator Plus工具顯示三維模擬效果，可以直觀判斷工程前、后模擬設(shè)計(jì)水位與河道堤頂、排水管涵及橋墩等位置關(guān)系，為實(shí)際工程建設(shè)提供科學(xué)依據(jù)，有利于提升橋梁工程防洪評價(jià)分析計(jì)算水平。