摘要：為分析“方廳水院”工程對國家基本水文站金澤站水文監(jiān)測的影響程度，利用Mike21水動力模型，對不同水情下工程對該站水位、流量和流場的影響進行計算模擬。結(jié)果表明：工程實施前后水文站斷面水位最大變幅均在2 mm以內(nèi)，工程對水位觀測結(jié)果存在一定的影響，但幅度較?。凰恼緶y驗斷面流速分布發(fā)生變化，各水情條件下，流速變幅大于3%的點位占比16%～33%，改變了原有代表流速關(guān)系曲線；主流附近的點位流速略有加快，幅度不超過2%，改變了原有的代表流速關(guān)系曲線?！胺綇d水院”工程的實施對金澤水文站的水文監(jiān)測存在一定程度影響，提出了降低影響程度的補償措施，以保障水文資料一致性。

關(guān)鍵詞：水文監(jiān)測； Mike21模型； 水位； 流量； 流場； 金澤水文站

中圖法分類號：P332

文獻標志碼：A

DOI：10.15974/j.cnki.slsdkb.2024.09.005

文章編號：1006-0081（2024）09-0030-06

0 引 言

“方廳水院”是長三角生態(tài)綠色一體化發(fā)展示范區(qū)水鄉(xiāng)客廳的核心地標場所，工程塑造世界級水鄉(xiāng)人居文明典范，促進長三角地區(qū)文化交流與創(chuàng)新發(fā)展，集中展示一體化發(fā)展成就和江南水鄉(xiāng)文化，是列入《長三角一體化發(fā)展規(guī)劃“十四五”實施方案重大項目庫》中的核心項目。 根據(jù)《水文監(jiān)測環(huán)境和設(shè)施保護辦法》，“方廳水院”步行橋工程位于國家基本水文站金澤站水文監(jiān)測河段環(huán)境保護范圍內(nèi)，本文旨在利用平面二維水流數(shù)學模型Mike21分析 “方廳水院”步行橋工程建設(shè)對金澤水文站測驗設(shè)施、測站工作環(huán)境和水文資料整編方法等造成的影響。

橋梁工程的影響分析主要有經(jīng)驗公式法、數(shù)值計算法和水工模型法，隨著現(xiàn)代計算機技術(shù)與數(shù)值方法的不斷發(fā)展，水動力數(shù)值模型應(yīng)用前景廣闊［1］。其中，平面二維水流數(shù)學模型Mike21以垂線平均的水流因素作為研究對象，能夠較好地用于流體動力學、波浪、沉積物動力學、水質(zhì)和生態(tài)學的建模［2］，是目前國際上同類模型中功能比較齊全、應(yīng)用比較廣泛的模型，在對橋梁工程進行水位和流場模擬中已得到廣泛應(yīng)用［3］。目前，該模型在國內(nèi)的應(yīng)用發(fā)展也很迅速，在長江口綜合治理工程、南水北調(diào)工程、河北滄州防洪影響分析、昆山城市內(nèi)澇治理等實踐和研究中應(yīng)用廣泛［4-5］。

傳統(tǒng)Mike21水動力模型采用矩形網(wǎng)格計算模塊進行有限差分法，對海岸或防波堤等不規(guī)則邊界常處理成齒狀，計算結(jié)果不盡理想。矩形網(wǎng)格計算模塊難以進行小尺寸局部水工建筑物的繞流模擬。模型屬大范圍平面二維數(shù)學模型，不適用于模擬三維問題或急流等流動［6］。本文研究采用Mike21水動力模型三角形網(wǎng)格對建橋后計算區(qū)域進行剖分，對橋墩位置附近的網(wǎng)格局部加密來改善優(yōu)化模型，計算模擬金澤站在“方廳水院”工程運行期與施工期4種不同水情條件下的水位、流量、流場監(jiān)測的受影響程度。

1 測站及工程概況

1.1 測站概況

金澤水文站為國家基本水文站，為太湖流域黃浦江水系太浦河干流控制站，屬蘇浙滬省界控制站。金澤水文站位于太浦河下游，屬于黃浦江上游水系，距河口17 km，與“方廳水院”工程相對位置見圖1。該站監(jiān)測項目包括水位、流量、水質(zhì)、雨量和氣象。金澤站測驗河段較為順直，為感潮河段。受潮汐影響，每日水位呈現(xiàn)兩高兩低態(tài)勢，測驗斷面呈U形，兩岸為石駁岸，河床穩(wěn)定，過往船只較多。河槽河床質(zhì)為黏土，多年沖淤變化小。金澤站洪水來源于上游太湖洪水下泄及兩岸地區(qū)排澇，主要由降水產(chǎn)生，同時又受黃浦江潮水頂托影響，斷面流向為往復流，洪水具有洪峰歷時短、高水位持續(xù)時間長、漲潮快落潮慢等特點。河道漂浮物較少，一般以水葫蘆等水生植物為主。測站流量采用H-ADCP在線監(jiān)測，代表流速歷年關(guān)系較穩(wěn)定。

1.2 工程概況

“方廳水院”工程位于蘇浙滬兩省一市行政邊界交匯處，以長三角原點為中心，跨越太浦河，連通滬蘇浙三地。“方廳水院”采用四合院形制，四角分別設(shè)置了滬蘇浙皖的主題展館，在太浦河南岸以長三角原點為中心沿中軸線設(shè)置主會館。滬蘇浙皖4個主題館通過3座分別跨越太浦河及與元蕩漾聯(lián)通水系的步行橋相連?？绾訕蛄号c建筑之間無縫銜接，實現(xiàn)了跨域、跨河互連。此外，擬對水文站建筑外觀進行修繕翻新，以保證建筑與“方廳水院”的整體氛圍相融合。工程還將實施太浦河兩岸護岸拆除重建，并根據(jù)工程占用水域面積新開挖水域，起到水域補償?shù)淖饔谩?/p>

金澤水文站位于“方廳水院”東、西兩座步行橋之間，水文測驗斷面距離東步行橋約158 m，距離西步行橋約105 m，北步行橋位于現(xiàn)狀水文站場區(qū)北側(cè)約30 m?！胺綇d水院”步行橋工程位于金澤水文站水文監(jiān)測河段環(huán)境保護范圍20 km之內(nèi)，工程的建設(shè)可能對水文站測驗設(shè)施、測站工作環(huán)境和水文資料整編方法等存在直接影響。

2 數(shù)學模型分析

2.1 模型基本方程

采用丹麥水利研究所開發(fā)的Mike21水動力模型系統(tǒng)，分析工程對金澤站水文監(jiān)測的影響程度。模型利用ADI二階精度的有限差分法求解，基本方程考慮布辛涅司克（Boussinesq）近似和淺水假定，以及風應(yīng)力的影響，垂向積分的二維水動力學方程組［7-8］如下。

連續(xù)方程：

ξt+px+qy=S（1）

動量方程：

pt+xp2h+xpqh+ghξx+gp p2+q2C2h2－1ρwxhτxx+yhτxy-Ωq-fVVx+hρwxpa=0（2）

qt+yq2h+xpqh+ghξy+gq p2+q2C2h2－1ρwyhτyy+xhτxy+Ωq－fVVy+hρwypa=0（3）

式中：h為水深，m；ξ為水位，m；p，q為x，y方向的單寬流量，m2/s；C為謝才系數(shù)；f為風阻力系數(shù)；V，Vx，Vy為風速及其在x，y方向上的分量，m/s；Ω為科里奧利 （Coriolis）參數(shù)；pa為大氣壓，Pa；ρw為水的密度，kg/m3；τxx，τxy，τyy分別為剪切應(yīng)力分量。

2.2 模型建立

2.2.1 模型概化

計算范圍為“方廳水院”步行橋上游3.0 km、下游2.5 km河道范圍，采用三角形網(wǎng)格對建橋后計算區(qū)域進行剖分，共22 634個網(wǎng)格，11 869個節(jié)點。橋墩采用陸地的方式進行概化，新建東、西步行橋概化的橋墩形態(tài)和面積與實際尺寸一致，為2.5 m×8.6 m，共計兩組4個橋墩。施工圍堰采用陸地的方式進行概化，概化圍堰面積與實際施工圍堰平面尺寸一致。其中，東、西步行橋單個橋墩的圍堰尺寸為14 m×14 m，北步行橋的圍堰尺寸為10.2 m×17.25 m。水域補償包括北橋橋墩、東橋西南側(cè)、西橋南側(cè)3處，斷面補償包括東橋南側(cè)、西橋南側(cè)。

2.2.2 參數(shù)設(shè)置

（1） 邊界條件及初始條件。采用不同水情下的實測數(shù)據(jù)作為邊界條件，模擬計算工程前后水位、流場變化情況，邊界條件見表1。

（2） 糙率系數(shù)。采用曼寧糙率系數(shù)用于模型計算，糙率系數(shù)取值參考相關(guān)研究成果，不同區(qū)域取不同的值。

（3） 渦黏系數(shù)。橫向渦黏系數(shù)采用湍流公式（Smagorinsky）估算，相應(yīng)取值為0.28 m2/s，水平渦黏系數(shù)則采用對數(shù)公式估算。

（4） 時間步長。根據(jù)模型網(wǎng)格大小、水深條件動態(tài)調(diào)整模型計算時間步長，使庫朗數(shù) （ CFL數(shù)）小于0.8，滿足模型穩(wěn)定的要求，計算時步長為0.01～30 s。

（5） 干濕邊界處理。對計算區(qū)域內(nèi)灘地干濕過程采用網(wǎng)格凍結(jié)方法處理，當網(wǎng)格點水深小于0.005 m時，令該點為干點，灘地干出，不參與水動力計算；當網(wǎng)格點水深為0.005～0.05 m時，令該處流速為零，該點僅參與水流連續(xù)方程的計算，當該處水深大于0.1 m時，該網(wǎng)格點參與計算，水流上灘。

2.3 模型計算成果

根據(jù)數(shù)學模型計算分析，結(jié)合工程建設(shè)特點，分別對運行期和施工期的項目建設(shè)影響進行分析。選取模擬區(qū)域橋墩上下游100 m范圍內(nèi)12個斷面，84個代表點位分析不同水情條件下工程前后水位變化情況。

2.3.1 水位影響計算

（1） 運行期，橋墩上下游50 m以外水域范圍水位變化不足1 mm，水文站斷面水位變幅最大不足1 mm。因此，工程建設(shè)對水文站水位觀測結(jié)果存在一定影響，但幅度較小。計算結(jié)果見表2。

（2） 施工期，工程范圍內(nèi)水位變化情況見表3。各種水情條件下，施工期水文站斷面水位最大變幅均在2 mm以內(nèi)，對水位觀測結(jié)果存在一定的影響，但幅度較小。

2.3.2 流量影響計算

（1） 運行期，根據(jù)太浦河二維水流模型計算結(jié)果分析，各種水文條件下，工程實施后新建步行橋引起的水頭損失均較小，流量變化影響甚微，見表4。

（2） 施工期，圍堰施工期間水頭損失略大于運行期，詳情見表5。

2.3.3 流場影響計算

（1） 運行期流速變化的計算結(jié)果見表6。測驗斷面南、北兩端各35 m范圍受橋梁阻水影響，流速變化明顯，各種水文條件下，流速最大升、降幅在10%～15%；位于河道中央的點位流速略有加快，變化多在1%以內(nèi)。流向變化：各種水文條件下，測驗斷面各點位流向變化一般不超過1°，流向變化較小。

（2） 施工期圍堰施工后，水流在圍堰上游分流，在圍堰東側(cè)水域形成紊流。計算結(jié)果見表7。由于圍堰尺寸較大，施工期間金澤站測驗斷面各點位流速、流向變化相對偏大。對于流速變化，測驗斷面南、北兩端受圍堰阻水影響，流速減慢明顯，各種水文條件下的流速最大降幅在15%～40%；位于河道中央的點位流速加快，變化多在3%以內(nèi)。對于流向變化，各種水文條件下，測驗斷面除南北兩端外，其余點位流向變化一般不超過1°，流向變化較小。

3 工程建設(shè)對金澤水文站的影響分析

3.1 工程對觀測設(shè)施的影響分析

（1） 水位觀測設(shè)施。根據(jù)數(shù)模計算成果，金澤水文站水文測驗斷面南、北兩端受圍堰阻水影響，流速減慢明顯，位于河道中央的點位流速加快。圍堰施工過程中，近岸流速減緩，加之人為因素導致水體含沙量增大，有可能導致水文站水位井的進水口淤積堵塞。

（2） 雨量觀測設(shè)施。工程對人工雨量筒及自記雨量計等基礎(chǔ)設(shè)施不會產(chǎn)生明顯影響。

（3） 流量觀測設(shè)施。工程距水文站的測驗斷面較近，對流量觀測設(shè)施的影響僅限于施工期間。施工期間水面駐留的大型施工船舶及施工過程產(chǎn)生的振動，可能會對流量觀測設(shè)施產(chǎn)生不利影響，如H-ADCP安裝支架的穩(wěn)定性、儀器工作井的牢固性等。

3.2 工程對水文站測驗工作環(huán)境的影響分析

3.2.1 對水位觀測的影響分析

根據(jù)模擬計算成果，工程建設(shè)前后，橋墩上下游50 m以外水域范圍水位變化不足1 mm，水文站斷面水位變幅最大不足1 mm，工程建設(shè)對水文站水位觀測結(jié)果存在一定影響，但幅度較小?？紤]到東、西兩座橋在河道內(nèi)共布置4座主橋墩，占用過水面積較大，且自由浮動的柔性防撞墩進一步加大了阻水面積，橋墩的建設(shè)會造成上游水位壅高，下游水位降低。而金澤水位站處于東、西兩座橋之間，且太浦河為感潮河段，存在水位的周期性漲落過程，工程建設(shè)勢必改變水位站位置的水流流態(tài)，改變了水位站上下游的水面比降沿程分布規(guī)律，對水位站的水位觀測帶來明顯不利影響，導致金澤水位站的水位資料一致性受到影響。因此，開展水位的對比分析，研究工程實施對水位的影響程度十分必要。

3.2.2 對雨量觀測的影響分析

工程施工過程中產(chǎn)生的粉塵、振動等問題可能會影響雨量器的觀測精度。

3.2.3 對流量測驗的影響分析

由于河段上下游流場發(fā)生變化，致使水文站測驗斷面流速分布變化，并引起水流關(guān)系特性變化，改變了原有的代表流速關(guān)系曲線，使原有流量資料連續(xù)性遭到破壞。水文站測驗斷面中位于橋墩下游的點位流速減緩最大為0.042 m/s，最大降幅約為9.42%；主流附近的點位流速略有加快，加快幅度不超過2%。工程實施后，水文站需要開展定線比測，重新率定流量推算關(guān)系。

3.2.4 對河道地形的影響分析

工程的建設(shè)，特別主橋墩的建設(shè)，改變了工程附近河段的流態(tài)，改變了原來流速流向沿斷面的分布規(guī)律。根據(jù)計算成果，東、西橋橋墩對水流形成束縛，主流流速將有所加快，而近岸流速將有所減小。因此，工程建設(shè)可能導致河道近岸淤積、河道中央沖刷。受水流流態(tài)變化的影響，工程附近的河道地形將發(fā)生變化，以重新適應(yīng)新的河道邊界。開展工程前、工程施工期間、工程完工時、完工后不少于1 a的河道地形監(jiān)測十分必要，同時應(yīng)對水文測驗斷面地形加密測量。

3.3 對資料整編的影響分析

工程建設(shè)對水位資料整編方法無影響，但工程建成后，需要采取措施對工程前后的水位進行對比分析研究和水位整編專題論證，否則可能會影響水位資料的一致性。工程建設(shè)對雨量資料整編基本無影響；流量資料整編時，水文斷面附近的水流流態(tài)將發(fā)生改變，為了確保水文站的流量觀測精度，需要重新率定出穩(wěn)定的指標流速和實測斷面平均流速的關(guān)系線，確保金澤水文站流量資料的連續(xù)性、可靠性、一致性。

4 結(jié)論及建議

“方廳水院”工程對金澤站水位、流量及流場存在一定影響，且施工期的影響要大于工程運行期，特別是圍堰施工期間，影響較為突出。對于水位監(jiān)測，施工前及整個施工期應(yīng)在河段適當位置設(shè)立臨時水位站，對臨時水位站與金澤水位站的水位同步比測，便于對工程前后水位資料進行銜接處理。在施工期及運行期適當增加對水位觀測設(shè)施的維護保養(yǎng)頻次，對水位井實施清淤、加固工程，并開展水位整編專題論證。對于流量、流場監(jiān)測，施工前應(yīng)布設(shè)臨時水文測驗斷面，開展臨時水文測驗斷面流量比測工作；施工期間開展臨時斷面與原斷面同步觀測，并收集金澤站下游約 5.2 km 處的東蔡水文站的同步流量成果，為修正現(xiàn)狀斷面的流量成果提供參考依據(jù)；運行期，在原址對流量測驗基礎(chǔ)設(shè)施進行加固改造。另外，為了研究全斷面覆蓋的測驗方式，建議在現(xiàn)狀斷面的右岸增設(shè)一套H-ADCP系統(tǒng)，收集垂向潮流資料，同時考慮到東、西兩座橋共4座橋墩的建設(shè)。河道過流斷面面積減小，對水流形成束縛，建議在主流兩側(cè)的安全位置布設(shè)2條坐底ADCP觀測系統(tǒng)，收集垂向潮流資料，與現(xiàn)狀左岸H-ADCP構(gòu)建新的流量觀測系統(tǒng)［9］，維持金澤水文站的流量測驗精度。“方廳水院”工程項目完成后，經(jīng)多方比測，若金澤站依然無法率定出穩(wěn)定的代表流速關(guān)系線，建議在新的河道邊界條件下，構(gòu)建流量調(diào)算模型，反演現(xiàn)狀斷面流量，消除橋墩對測驗資料的影響。

參考文獻：

［1］ 任梅芳，徐宗學，蘇廣新.基于二維水動力模型與經(jīng)驗公式的橋梁壅水計算及其對比分析［J］.水力發(fā)電學報，2017，36（5）：78-87.

［2］ 胡麗琴，王進，陳雙飛.基于MIKE21模型對橋梁施工期洪水影響分析［J］.水利科技與經(jīng)濟，2023，29（10）：92-97.

［3］ 杜麥，陳小威.基于MIKE21的城市橋梁防洪影響分析研究［J］.中國防汛抗旱，2023，33（10）：73-77.

［4］ 許婷.丹麥MIKE21模型概述及應(yīng)用實例［J］.水利科技與經(jīng)濟，2010，16（8）：867-869.

［5］ 侯俊，陳誠，鄭玉磊，等．基于Mike Flood模型的平原水網(wǎng)地區(qū)內(nèi)澇治理效果模擬研究——以昆山市淀山湖鎮(zhèn)為例［J/OL］．水利水電科技進展.https：∥link.cnki.net/urlid/32.1439.TV.20240522.0913.002

［6］ 袁雄燕，徐德龍.丹麥MIKE21模型在橋渡壅水計算中的應(yīng)用研究［J］.人民長江，2006，37（4）：31-32，52.

［7］ 彭興，馬祺瑞，果鵬.雅魯藏布江江熱橋防洪影響評價分析［J］.水利水電快報，2023，44（增2）：5-10．

［8］ 周雁潭，劉立軍，李鵬輝，等.基于多層級目標耦合的城市內(nèi)澇防治研究［J］.人民長江，2024，55（5）：15-22，32.

［9］ 秦海斌，李偉.鳳凰頸排灌站自動流量比測率定分析［J］.陜西水利，2023（8）：46-49.

（編輯：李 慧）

Impact analysis of Fangting Water Institute Project on hydrometry

SONG Chengwei

（Shanghai Hydrological Survey Team of Qingpu District，Qingpu District Water Bureau，Shanghai 201700，China）

Abstract：

To analyze the impact of Fangting Water Institute Project on the hydrological monitoring of the national basic Jinze Hydrological Station，the advanced Mike21 hydrodynamic model was used to simulate the impact of the project on water level，flow and velocity field under different flow conditions.The results showed that the maximum variation of the cross-section water level of the hydrological station before and after the implementation of the project was within 2 mm，which had a relatively small impact on the water level observation results.The distribution of flow velocity in the observation section of the hydrological station was changed，among which，the velocity of 16% to 33% of the measuring points changed with a margin larger than 3%.The flow velocity near the mainstream was slightly increased，with an increase of no more than 2%，changing the original representative flow velocity relationship curve.Therefore，the project will have an inevitable impact on the hydrological monitoring of Jinze Hydrological Station，and further compensation measures was proposed to reduce the impact and ensure the consistency of hydrological data.

Key words：

hydrometry； Mike21 model； water level； flow； velocity field； Jinze hydrological station