福州市閩江西河游泳場河段水流特性整治工程數學模擬研究

丁 堅，邵雨辰，郭 淮，王 霄

(1.河海大學 海岸災害及防護教育部重點實驗室，南京 210098；2.蘇交科集團股份有限公司，南京 210098；3.淮安市水利勘測設計研究院有限公司，淮安 223005)

福州市閩江西河游泳場河段水流特性整治工程數學模擬研究

丁 堅1，邵雨辰2，郭 淮3，王 霄3

(1.河海大學 海岸災害及防護教育部重點實驗室，南京 210098；2.蘇交科集團股份有限公司，南京 210098；3.淮安市水利勘測設計研究院有限公司，淮安 223005)

彎道環流、渦旋等不良流態往往是天然游泳場的安全的不利因素，不良流態的形成與河道的天然條件、人為活動兩方面因素相關。文章在分析閩江西河河段河勢、流態的基礎上，建立該河段的二維數學模型，針對不同的自然條件和整治方案，對該河段的水動力特性進行研究，在此基礎上，提出整治效果較優的河道整治方案。在天然河道游泳場整治過程中，除了考慮到對游泳場局部流場流態的優化調整之外，還要綜合考慮整治工程對河道通航、行洪的其他功能影響，進行綜合比選。

彎道環流；渦旋；河道整治

西河天然游泳場位于閩江福州北港洪山橋下游左岸(北岸)，長約1 000 m，目前來水清澈，沙質河床，具有開發為優良天然游泳場的基礎，但由于河勢、水流等因素制約了游泳場的發展，必須進行整治，以達到游泳場水流平順、流速適當、安全衛生等要求。為此需對該河段河床、水流進行全面分析，找出存在不足之處，對癥治理。

本文在分析閩江西河河段河勢、流態的基礎上，建立該河段的二維數學模型，針對不同的自然條件和整治方案，對該河段的水動力特性進行研究，在此基礎上，提出整治效果較優的河道整治方案。

1 河段概況

游泳場位于北港上段，主流自淮安分流后，進入北港，首經文山里，為窄深峽谷河段，河寬130～200 m，最窄處僅90 m，水流湍急，至洪山橋后，以下河道軸線走向左偏約27°，出舊洪山橋，河床明顯放寬至500m左右，來水流速減緩，至5 km以下蒼霞～解放大橋窄彎段時，河寬又明顯縮窄，在蒼霞河段處河寬僅200 m，加上解放大橋、江南橋阻水，使上游河段壅水，比降減小，以致舊洪山橋至解放大橋的水流進一步減緩，泥沙落淤，出現眾多邊灘、淺區及江心洲等。位于左岸的西河游泳場的凸岸邊灘也發育的比較大。通常河床的高程在-1.0～-3.0 m，但舊洪山橋在河心處有6個橋墩在歷年大洪水中被沖垮或拆除，為保福州市防洪安全，其余橋墩不能繼續垮塌，故多年來沿舊洪山橋軸線，在舊橋墩之間拋石護墩，沿橋軸線形成潛壩，有效地維持北港分流比穩定。由于中段沖掉六個橋墩，主流都要集中在該缺口下泄，將缺口下游河床沖深至-11.0～-14.0 m的深槽，其走向是沿著橋上段中軸線的方向，直逼橋下的右岸，故在20世紀90年代規劃防洪堤岸線時，確定了右岸的防洪堤(圖1)。

由于北港舊洪山橋以下河面放寬一倍以上，而河軸線又向左偏轉了27°，加上舊洪山橋軸線走向是斜交于河床斷面的，舊洪山橋起著潛壩的作用，來水溢過潛壩時，主流將趨向垂直于壩軸線方向出流，這種情況下過壩的水面比降最大，水往低處流，垂直于壩軸線下泄是客觀規律，因此，過舊洪山橋的水流動力軸線更向右偏，右岸為主流頂沖的凹岸，西河游泳場的左岸則成為凸岸，由洪枯水的水流動力軸線就很清楚顯示為彎道水流，出現螺旋流，表面水體向右偏，水面呈現橫比降，右高左低，使下層含沙量大的水體向左移，因此，將泥沙帶向左岸形成邊灘，這就是彎道水流造成凹岸沖刷，凸岸淤積的規律[1-7]，游泳場邊灘就是該流勢塑造的結果，沿邊灘上下及前沿環流都較為明顯，且直接威脅游泳的安全(圖2)。

圖1 閩江西河河勢圖Fig.1 River regime of Minjiang West River

2 數學模型的建立和驗證

2.1平面二維數學模型控制方程

(1)

(2)

(3)

2.2方程定解條件

圖2 舊洪山橋下游洪枯水水流動力軸線Fig.2 Flow dynamic axis at the downstream of old Hongshan bridge

為了使平面二維水流運動的基本方程在計算域Ω內有適定解，必須滿足定解條件，即在計算域上給出初始條件和邊界條件。

2.3模型計算區域及網格

根據研究內容的要求，數學模型邊界上起文山里，下至解放大橋，沿程長約11 km，水域面積4.2 km2，可參見圖1。進行數學模型計算時，主要采用實測水下地形，并在模型中采用插值的計算方法將地形賦值至各計算節點上。模型計算區域的離散采用三角形網格，并在工程區域進行局部加密，模型最大網格尺度為40 m，工程區最小網格尺度為10 m，可以滿足計算精度要求，模型生成水下地形見圖3。

圖3 模擬地形Fig.3 Simulation of the terrain

2.4模型參數設置及驗證

2.4.1 模型參數設置

模型采用的時間步長△t=20 s。計算水域的糙率是一個綜合影響因素，是數值計算中一個十分重要的參數，與水深、床面形態等因素有關，閩江河段淺灘深槽交錯復雜，因此根據水深條件，給出隨水深變化的糙率場。糙率范圍取0.02～0.05，根據驗證計算情況要求選取。

2.4.2 模型驗證

根據本次研究的要求與所掌握的實測資料，模型驗證計算時段：中水大潮期：9月18日11：00時～24：00時；中水小潮期：9月28日11：00時～23：00時。驗證位置選取參見圖1，分別為清水碼頭站、西河1#點、西河2#點、 在數學模型的驗證過程中，分別進行了潮位、潮流過程的驗證。潮位驗證圖見圖4～圖5,由圖可見，潮位模型計算值與實測值吻合較好。流速驗證圖見圖6，由圖可見,西河1#點、西河2#點的漲潮流速、落潮流及流向計算值與實測值吻合較好。由此說明建立的數學模型能較好地模擬西河游泳場河段的水流運動， 可用于工程方案的數值計算。

4-a清水碼頭站 4-b上游廠碼頭站圖4 大潮潮位驗證圖Fig.4 Verification of tidal level during spring tide

5-a清水碼頭站 5-b上游廠碼頭站圖5 小潮潮位驗證圖Fig.5 Verification of tidal level during neap tide

6-a 1#點大潮流速 6-b 1#點大潮流向

6-c 1#點小潮流速 6-d 1#點小潮流向

6-e 2#點大潮流速 6-f 2#點大潮流向

6-g 2#點小潮流速 6-h 2#點小潮流向

圖6 西河流速流向驗證圖
Fig.6 Verification of flow velocity and direction at West River

3 無工程措施時流場分析

圖7 西河游泳場河段漲急及落急流場圖Fig.7 Flow field at the times of flood maximum and ebb maximum in West River Natatorium reach

圖7為西河游泳場河段大、小潮漲急及落急流場圖(9月18日及9月28日)。由圖可見，漲急時西河游泳場河段水流較平順，僅在舊洪山橋左右兩側下游出現2個較小的漩渦；落急時西河游泳場河段水流特別紊亂，特別大潮落急在舊洪山橋左右兩側下游出現2個較大的漩渦，游泳場中段水域也存在漩渦，游泳場下游段深水區流速較大，所以大潮落急時期的流場是本次河道整治關注的重點。

4 工程措施及整治效果分析

4.1工程措施

方案1：首先必須對近1 000 m長的西河游泳場水域水下地形進行平整邊坡。圖8為西河游泳場水域布置圖，圖中虛線為游泳場水域的邊界線，岸邊虛線為5 m等深線，游泳場的深水虛線為0 m等深線，游泳場約有170 m寬(圖8)，對該區域內的河床進行削坡、填沙處理，平順河床。水閘泵房前水域保留其出口。邊坡平整后，游泳場水域等深線與河道右側水流動力軸線平行，將有利于游泳場水域水流的平順。

圖8 西河游泳場水域布置圖Fig.8 Water area layout of West River Natatorium

方案2： 在方案1平整游泳場邊坡的基礎上，對舊洪山大橋主泓兩側橋墩間水深挖至+2.0 m水深處，在舊洪山大橋下游深槽處拋一潛堤。由于加大了舊洪山大橋兩側的分流，增大了兩側的流速，主橋墩間流速將減小。這將從根本上消除舊洪山大橋下游兩側的漩渦，減弱彎道環流強度，歸順游泳場前沿的水流。

方案3：該方案分兩個區進行游泳場邊坡修正，西河泵站上游游泳場前沿水深從1.0 m等深線以3%坡度向岸修正，泵站下游游泳場前沿水深從0.0 m等深線以3%坡度向岸修正，游泳場寬度約160 m。泵站出口水域保留原狀，并在出口兩側邊坡進行護坡處理。在舊洪山大橋下游深槽處拋石平整至-8.0 m，拋石方約40 000 m3。在舊洪山橋左側下游沙灘頭部，沿等深線做弧形壩， 用于游泳場上游頭部固灘， 在游泳場下游尾部做一丁壩(頂部高程高出沙灘1.5 m, 頭部伸出游泳場0.0 m線以外10.0 m)，防止游泳場沙灘沖刷；在游泳場前沿填鵝卵石護底，防止水流沖刷。該方案同樣對舊洪山橋墩間不進行開挖維持原狀，保證北港分流比不變，對福州市區的防洪是有利的。

圖9 方案2潛堤位置圖Fig.9 Location of submerged bank of Plan 2

圖10 方案3工程布置圖Fig.10 Engineering layout of Plan 3

4.2整治效果分析

4.2.1 工程后流場變化

方案1游泳場邊坡修整后，游泳場原水閘出口外水域漩渦減弱，游泳場水域流場較邊坡修整前歸順。方案2實施后，加大了左右側分流及流速，減弱了主泓的分流及流速。使得流量在舊洪山大橋斷面分配更均勻，減小了主橋墩間的流速，減弱了彎道環流，使得游泳場前沿水域流場更歸順。該方案將有效地消除橋下兩側漩渦，減小彎道環流，歸順游泳場水域的水流流態。方案3分兩個區進行游泳場邊坡修正，游泳場邊坡修整后，與其他方案相比較，游泳場前沿水域流場比其他方案更歸順。方案3西河泵站出口下游游泳場寬度約160 m，游泳場前沿離航道比方案1和方案2要遠，減小船行波對游泳場沙灘的作用，也增大了游泳的安全性。該方案在主流深槽處拋石整平，減小河床底部水深的梯度變化，有利于平順主流區的流速流向。由于對舊洪山橋墩間不進行開挖，對閩江北港分流比無影響，對福州市區的防洪是有利的，但該方案對消除舊洪山橋下漩渦、減小彎道環流的效果不如方案2(圖11)。

圖11 各方案西河游泳場河段中水大潮落急流場圖Fig.11 Flow field of West River Natatorium reach in the median water discharge condition during spring tide at the times of ebb maximum in each plan

4.2.2 工程后流速變化

12-a 游泳場邊沿

12-b 游泳場中心圖12 游泳場河段流速分布Fig.12 Velocity distribution of natatorium reach

對河段工程前后流速進行提取比較(見圖12,圖中橫坐標1～11為游泳場上游頂端至下游頂端11個流速測點，相鄰測點距離約為100 m)，游泳場邊坡進行修整后(方案1)，與工程前比較，邊坡修整后斷面形態發生變化，西河泵站出口上游工程河段，流速沿斷面發生調整，左岸削坡區游泳流速增大，右岸動力、河中相對流速減小。泵站出口下游工程河段，上游過水斷面加大，使得來水加大，左岸場中流速減小，場邊沿、右岸動力、河中的流速均有所增大。方案2與方案1比較，舊洪山橋的阻水減弱明顯，橋下游的水流在西河泵站以上河段有所減小，西河泵站出口以下河段流速均有所增大。方案3與方案1相比，場沿點的流速相差不多，游泳場上游區水深要比前兩種方案要淺，場中點流速游泳場上游區要比方案1小一點，而游泳場下游區要比方案1比方案2大一點，該方案右岸點流速要比方案1和方案2略小，比工程前略大。

4.2.3 工程后水面線變化

西河游泳場河段治理工程(舊洪山大橋挖深、平整邊坡、潛堤、丁壩等)實施后，勢必會影響河流的過水斷面面積，對閩江北港及福州市區的防洪影響是工程方案確定前必須要評估的問題。通過模型模擬文山里一般豐水流量條件下(Q=4 000 m3s)，給出舊洪山大橋上500 m至橋下1 500 m河段工程方案后水面線對比圖(見圖13，橫坐標0點為舊洪山橋上游500 m處)。由圖可見，方案2實施后加大了舊洪山橋的過流能力，對水面線的改變較大，工程實施后舊洪山橋上水位下降約8 cm, 舊洪山橋下游水位升高約10 cm, 故加大了舊洪山橋下游福州市區的防洪壓力。對于方案1、3，由于沒有對舊洪山橋橋墩間水深進行開挖，河流的過流能力變化很小，對水面線的改變影響也很小，故對福州市區的防洪幾乎沒有影響。

圖13 工程區域河段一般豐水水面線對比圖Fig.13 Comparison of water surface profile of the general large flow in the project area

4.3推薦方案

本節提出的3種方案，對游泳場邊坡進行不同程度的平整，游泳場原水閘出口外水域漩渦減弱，游泳場水域流場較邊坡修整前歸順。方案2采用加大大橋兩側分流的方法，減小彎道環流及漩渦對游泳場產生的危害，效果較方案1和方案3要明顯，但方案2實際上增加了舊洪山大橋的過流能力，方案實施后大橋下游水面線會有較為明顯的抬高，對福州市區的防洪不利，故不采納。

由于西河泵站引水渠出口以上是淤積區，泵站出口以下為沖刷區，上游邊灘比下游邊灘要高，方案3采用游泳場上游外邊線定于+1 m高程，按3%坡度向岸邊填筑；下游段游泳場水域外邊線定于+0 m高程，向岸邊按3%坡度削坡，較方案1節省了工程量，且下游游泳場區離航道較遠，對游泳場的安全有利。此外，方案3采取了充填深坑的方法，抬高出舊洪山橋主流的水面高度，減弱跌水的射流動能，減小對右岸頂沖的強度，減小彎道環流[8]，保護右岸，較方案1為優，故本文推薦方案3。

5 結論

本文采用了數學模擬的方法，對西河游泳場河段的河勢、漲落急流場、斷面進行了分析，發現河段彎道環流及漩渦會對游泳場的安全造成威脅，是整治工程的重點。提出了3種整治工程方案，并對各方案實施后的流場流速進行了對比分析。方案2在舊洪山橋橋墩間進行開挖，將從根本上消除橋下的漩渦，減小彎道環流，歸順游泳場水域的水流流態，但存在防洪影響的問題。方案3對游泳場邊坡分區進行修整后，游泳場水域流場較邊坡修整前更歸順，下游場地離航道較遠，對游泳場的安全有利；在主流深槽處拋石整平，減小河床底部水深的梯度變化，有利于平順主流區的流速流向，整治效果較方案1為優。

因此，在天然河道游泳場整治過程中，除了考慮到對游泳場局部流場流態的調整之外，還要綜合考慮整治工程對河道通航、行洪的其他功能影響，并進行綜合比選，才能得出最優整治方案。

[1]張瑋，倪兵，陳乾陽.長江澄通河段通州西水道整治工程對分流比影響研究[J].水道港口，2013,34(1):39-44. ZHANG W,NI B,CHEN Q Y. Research on influence of regulation works on diversion ratio of Tongzhousha west channel in Chengtong reach,Yangtze River[J].Journal of Waterway and Harbor,2013,34(1):39-44.

[2]陳曉云.長江太倉—南通河段河道演變特性與航道治理思路分析[J].水運工程，2012(4)：102-110. CHEN X Y. Regulation scheme of Taicang-Nantong deep-water fairway in Yangtze River[J]. Port & Waterway Engineering, 2012(4)：102-110.

[3]仲志余，王永忠.論長江澄通河段的綜合治理與開發[J].人民長江，2009(6)：1-5. ZHONG Z Y, WANG Y Z. Discussion on comprehensive regulation and development of Chengtong reach of the Yangtze River[J]. Yangtze River, 2009(6)：1-5.

[4]張麗.彎道水流二維數值模擬及岸線整治的水動力分析方法研究[J].水道港口，2012,33(2):136-141. ZHANG L. 2-D numerical simulation of bend flow and research of hydrodynamic method on renovating bank[J].Journal of Waterway and Harbor，2012,33(2):136-141.

[5]唐金芝.河流橫比降試驗及改正方法[J].湖南水利，1994(6):22-24.

[6]唐宏濤，吳偉.河流育道水流特點與沖刷深底探討[J].浙江水利科技，2004(1):7-9.

[7]宋志堯.計算彎道環流的通用公式[J].泥沙研究，2003,8(4):19-23. SONG Z Y. The Universal Formula on Transverse Circulation in Channel Bends[J]. Journal of Sediment Research，2003,8(4):19-23.

[8]張麗萍,楊達源,朱大奎.長江渝宜河段深槽分布特征及成因[J].南京大學學報:自然學版，2001,37(3):311-316. ZHANG L P, YANG D Y,ZHU D K. Formation Causes and Distributing Characteristics of Deep Troughs in the Chongqing-Yichang Reaches of the Yangtze River[J]. Journal of Nanjing University:Natural Sciences,2001,37(3):311-316.

Flow characteristics and regulation engineering mathematics simulation research of Minjiang River near West River Natatorium in Fuzhou

DINGJian1,SHAOYu-chen2,GUOHuai3,WANGXiao3

(1.KeyLaboratoryofCoastalDisasterandDefence,MinistryofEducation,HohaiUniversity,Nanjing210098,China;2.JSTIGROUP,Nanjing210098,China; 3.Huai′anWaterConservancySurveyingandDesignInstituteCo.,Ltd.,Huai′an223005,China)

Adverse flow conditions such as circulation flow and vortex are often considered as the disadvantageous factors affecting the safety of the natural bathing beach, and the formation of adverse flow relates to both the natural channel condition and human activities. In this paper, the hydrodynamic characteristics of the river and the optimized channel regulation scheme were studied by setting up a two-dimensional mathematical model that based on the analysis of Minjiang River morphology and flow condition. Research results indicate that in the process of natatorium river regulation, in addition to considering the local flow field of natatorium flow optimization adjustment, other functions of the channel such as navigation, flood discharge must be considered in order to make a comprehensive comparison.

circulation flow; vortex; river regulation

TV 133;O 242.1

：A

：1005-8443(2017)04-0374-06

2016-10-10；

：2017-02-16

丁堅(1960-)，男，江蘇省南京人，高級工程師，主要從事港口與航道工程專業研究。

Biography：DING Jian(1960-)，male，senior engineer.