謝昕卓 陳偉 楊明強 王璐 楊菲菲 蔡新

摘要：

為模擬南京市三汊河河口閘附近水流流態(tài)，探究河口閘單孔開啟的可行性，采用Delft 3D軟件建立大范圍水深平均模型和局部三維數(shù)學(xué)模型，對三汊河河口附近水流流態(tài)進行數(shù)值模擬，研究最不利水位工況下的水流流態(tài)，并確定了閘孔不同運行模式對流速分布的影響范圍。進一步地，將河道底面流速與研究范圍內(nèi)泥沙起動速度進行比較，獲得了河口閘附近河道沖刷與淤積范圍的分布情況。研究成果可為研究河口閘單孔開啟提供依據(jù)。

關(guān) 鍵 詞：

水閘運行方式; 水流流態(tài); 流速分布; Delft 3D模型; 河道沖淤; 三汊河河口閘

中圖法分類號： TV131.4

文獻標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2021.07.026

0 引 言

三汊河河口閘位于外秦淮河入江口，是南京市秦淮河環(huán)境綜合整治工程的重要組成部分。河口閘采用雙孔護鏡門的閘門型式，在大閘門上疊加溢流活動小門，形式新穎獨特。河口閘工程投入運行以來嚴(yán)格按照相關(guān)規(guī)定，采用雙閘孔同時開啟的方式調(diào)度運行，并定期進行檢測，保證結(jié)構(gòu)安全，獲得了較多研究數(shù)據(jù)。

然而隨著城市化進程的加快，外秦淮河水環(huán)境需要進一步改善。河口閘需要增加常態(tài)引調(diào)水功能，進而保證外秦淮河武定門閘至三汊河入江口段水體水質(zhì)穩(wěn)定滿足Ⅴ類水及以上標(biāo)準(zhǔn)。基于外秦淮河新的功能需求，需要在現(xiàn)行條件下分析論證河口閘單孔開啟運行的可行性，達到節(jié)約能源與高效利用水資源的目的。

由于河口閘單孔運行方式超出了原工程設(shè)計工況范疇，未來水閘能否實施單孔開啟運行，需要通過科學(xué)有效的手段對該運行模式進行論證。為此需要進行河口閘單孔開啟工況下河道的流態(tài)研究，計算其水流流速，判斷其沖淤狀況，為判斷單孔開啟是否符合安全性要求提供依據(jù)。

對于此類問題，可以使用Delft 3D軟件建立數(shù)值模型的方法加以研究。國內(nèi)外學(xué)者采用此類方法進行流態(tài)研究，取得了較好的效果。唐洪武等[1]曾在2011年利用Delft 3D軟件建立數(shù)學(xué)模型，對三汊河河口閘雙孔開啟運行工況下閘下沖淤問題進行了系統(tǒng)研究。胡孜軍等[2]通過建立洪澤湖區(qū)二維水深平均模型與三河閘閘室局部三維全耦合計算模型，對洪澤湖區(qū)三河閘附近的過流流態(tài)進行了精細(xì)化模擬。雷曉玲[3]利用Delft 3D數(shù)學(xué)模型對三峽航道環(huán)保疏浚的水質(zhì)變化進行了數(shù)值模擬，并對三峽航道疏浚工作的成果進行了定量評估。左書華[4] 用Delft 3D軟件建立水動力模型，研究了鰲江口外平陽咀海域流場，計算結(jié)果與實際測量結(jié)果吻合良好，證明此模型可以為工程研究提供流場背景。廖庾強[5] 利用Delft 3D建立了柳河彰武新城段平面二維模型，得到了水動力特征、典型斷面形態(tài)變化、河道泥沙分布等結(jié)果，為柳河彰武新城段防洪工作提供了理論基礎(chǔ)。范翻平[6]曾用Delft 3D軟件建立鄱陽湖水動力模型，研究典型風(fēng)向、風(fēng)速對鄱陽湖流場的影響，為研究鄱陽湖的污染問題提供了理論支撐。趙明登等[7]將ArcGIS與Delft 3D結(jié)合，對渭河下游泛洪區(qū)洪水?dāng)?shù)值進行研究，實現(xiàn)了河道及潰堤洪水演進的實時動態(tài)演示。郎韻、陳偉等[8-9]曾利用Delft 3D建立了江-河-水閘大范圍水深平均數(shù)學(xué)模型，研究三汊河河口閘左右閘孔過流嚴(yán)重不均的問題，發(fā)現(xiàn)導(dǎo)致過閘水流不均的主要因素是河勢及河道深泓位置，并提出設(shè)置導(dǎo)流墻可以改善過閘水流不均，該模型取得了較好的模擬效果。

綜合來看，前人利用Delft 3D數(shù)值模擬方法對河道水動力展開研究已獲得滿足要求的成果，但前人的研究大都是在河口閘兩側(cè)閘孔全部開啟或關(guān)閉工況下開展的，對于河口閘單孔開啟工況下的水動力特性，還未獲得滿足要求的成果。本文在前人基礎(chǔ)上，運用Delft 3D建立數(shù)值模型的方法，對三汊河河口閘單孔開啟工況下的河道水動力情況進行計算。

1 工程概況

三汊河河口閘位于南京市外秦淮河入江口，上游距下關(guān)大橋約300 m，下游距三汊河口約150 m。三汊河河口閘是秦淮河環(huán)境整治工程的重要組成部分，其主要功能是非汛期關(guān)閘蓄水，抬高外秦淮河水位，同時形成親水景觀，改善城市水環(huán)境和城市形象;汛期來臨時則開閘放水，不影響外秦淮河行洪。

河口閘縱向軸線順?biāo)鞣较虿贾?，工程范圍?nèi)軸線長179 m。閘室采用整體式大跨度鋼筋混凝土塢式結(jié)構(gòu)，順?biāo)鞣较蜷L37 m，垂直水流方向長97 m。閘室共兩孔，單孔凈寬40.0 m，中閘墩寬4.0 m，邊閘墩寬4.5 m，單孔總寬度48.5 m。閘底板高程1.0 m（吳淞基面，下同），閘底板厚2.5 m，閘墩頂高程7.5 m?？紤]到與上下游河道平順連接和水流流態(tài)，閘上、閘下分別設(shè)漸變段。上游漸變段長50 m，包括長度為20 m的鋼筋混凝土護坦和長度為30 m的素混凝土護底，兩側(cè)為混凝土預(yù)制塊護坡和生態(tài)混凝土護坡;下游漸變段長92 m，包括長度為20 m的鋼筋混凝土消力池、長度為60 m的混凝土護底和長度為12 m的拋石防沖槽，兩側(cè)為混凝土方格護坡和生態(tài)混凝土護坡（見圖1）。

三汊河河口閘非汛期正常過水流量為30 m3/s，行洪流量為80 m3/s（關(guān)閘蓄水狀態(tài)）;汛期閘門開啟，行洪流量為600 m3/s。工程自2005年建成運用以來，總體運行情況良好，在提升外秦淮河的水環(huán)境狀況方面發(fā)揮了重要作用，已成為南京市標(biāo)志性水利工程。

2 數(shù)值模型的建立

為研究河口閘單孔開啟的流場流態(tài)，現(xiàn)建立大范圍水深平均數(shù)學(xué)模型，包括外秦淮河部分河段，以全面反映外秦淮河水動力條件的不同組合對樞紐河段流動特性的影響。

式中：u，v，ω分別為ξ，η，σ方向上的流速;Gξξ、Gηη為坐標(biāo)變換拉梅系數(shù);ζ為相對于模型參照平面的水深，下同。

通過連續(xù)方程可以計算垂向流速w：

式中：H為總水深，下同。

（2） 動量方程。

ξ和η方向的動量方程分別為

式中：ρ為水體密度，除了斜壓項，其他處忽略密度的變化;Pξ，Pη分別為ξ，η方向上的靜水壓力梯度;Fξ，F(xiàn)η分別為ξ，η方向上紊動動量通量;vV為垂向渦動黏性系數(shù);fu，fv為柯氏力系數(shù);Mξ，Mη分別為ξ，η方向上動量的源匯項。

（3） 邊界條件。

初始條件：

邊界條件：

出、入流開邊界上，給定水位、流速或流量過程。固壁邊界采用無滑移邊界條件。

時步長和空間步長滿足穩(wěn)定性條件：

2.2 計算區(qū)域的選取和數(shù)值模型的建立

為充分考慮長江水位流量對河口閘過流特性的影響，大范圍水深平均模型計算區(qū)域包括外秦淮河、河口閘工程及長江，并考慮長江中潛洲、江心洲所在位置對邊界選取的影響[9]。最終確定長江上游邊界取為河口上游約1 100 m處，長江下游邊界為河口下游900 m處，外秦淮河上游邊界根據(jù)實測地形資料的范圍確定為河口閘軸線以上約820 m處。水深平均模型計算區(qū)域及模型網(wǎng)格劃分如圖2所示。

其中，長江上游邊界使用流量邊界條件，長江下游邊界使用水位邊界條件，外秦淮河上游邊界使用流量邊界條件。

依據(jù)河口閘已有研究成果的綜合分析，針對河口閘可能的單孔開啟運行模式，最終擬定研究工況組合如表1所列。其中組合模式1為河口閘可能運行的最不利工況組合，組合模式2為組合模式1的補充。

由圖3可見，河口閘數(shù)值模擬結(jié)果與河工模型試驗結(jié)果基本吻合。最大絕對誤差僅約 0.1 m/s，出現(xiàn)在距上游 450 m 斷面的近左岸。驗證結(jié)果說明該模型的模擬結(jié)果基本可靠，可用于后期的研究和分析。

2.4 泥沙起動流速研究

目前對于泥沙沖淤的研究，多采用泥沙拖拽力和起動流速作為標(biāo)準(zhǔn)。鑒于流速場和剪力場之間存在著一定的關(guān)系，更多學(xué)者為方便使用，相對泥沙起動拖拽力而言，更熱衷于推求起動流速[10]。

本文將泥沙起動流速作為泥沙沖淤的判別標(biāo)準(zhǔn)，若河道底面流速大于泥沙起動流速，則判定為該位置河道沖刷。

目前就天然均勻沙起動公式來看，起動流速公式在形式上差別不大，一般與粒徑的1/3指數(shù)方、水深的1/6指數(shù)方成正比，但其他系數(shù)取值離散度較大[11]。

對于細(xì)顆粒泥沙的起動流速，中國學(xué)者展開了卓越的研究，其中竇國仁公式[12]、張瑞瑾公式[13]、唐存本公式[14]、沙玉清公式[15]等都有很大的影響。

對眾多公式進行比較后，從實驗驗證和工程應(yīng)用兩方面考慮，本文采用唐存本公式：

式中：m=4.7（hd）0.06，γs為泥沙容重，γ為水流容重，δ為淤積物的干容重，δ0為淤積物穩(wěn)定干容重，h為水深，d為泥沙粒徑。本文

取γs=2.65 g/cm，γ=1.00 g/cm，c=2.9 ×10-4 g/cm，ρ=1.02×10-3 g·s/cm3，取穩(wěn)定容重情況δ=δ0。

采用研究段的流速變化對河口處的河床穩(wěn)定性進行分析。通過模型計算得出觀測點的流速及水深情況，推導(dǎo)出對應(yīng)的泥沙起動粒徑，進而分析研究段泥沙的起動概率及變化范圍。

根據(jù)資料，秦淮河河段泥沙平均粒徑為0.18 mm[1]。河段的起動流速擬選用這一特征粒徑進行計算。

兩種組合模式下計算范圍內(nèi)河道水深及泥沙起動流速如表2所列。

3 預(yù)設(shè)水位流量下河口閘水流流態(tài)的計算

為直觀比較每種組合模式下不同工況水流流態(tài)差異，現(xiàn)將每種組合模式分為3種工況，選取典型斷面對斷面上流速分布進行分析比較。為了獲得河道底面流速，建立了局部三維數(shù)據(jù)模型，對河口閘附近河道沖淤狀況進行了分析。

3.1 3種工況下流速分布對比

為了對比分析，將組合模式1細(xì)分為工況1、工況2和工況3;將組合模式2細(xì)分為工況4、工況5和工況6，各個工況詳細(xì)信息如表3所列。

為直觀比較不同工況下外秦淮河水流流速分布，現(xiàn)選取多個典型斷面，對斷面流速橫向分布進行分析。選取斷面為：閘上游200 m、閘上游100 m、三汊河河口閘軸線處、閘下游80 m、閘下游100 m、河口處（閘下游約150 m）。

組合模式1各個典型斷面橫向流速分布曲線如圖4所示，

組合模式2各個典型斷面橫向流速分布曲線如圖5所示。

通過兩種組合模式下的斷面流速分布，可以得到如下規(guī)律性認(rèn)識：

單孔開啟對流速分布的影響主要出現(xiàn)在閘上100 m斷面至閘下150 m斷面范圍內(nèi)，在閘軸線斷面處流速差異達到最大，越遠離閘軸線流速分布差異越小。

3.2 研究范圍內(nèi)河道沖淤

前文所作分析皆基于大范圍水深平均模型，其優(yōu)勢在于可以分析較大范圍內(nèi)水流流態(tài)變化情況，但對于流速分析，不能提供充足數(shù)據(jù)。為提取沖淤分析所需的河道底面流速，也為了分析河道過水面積變化與流速分布的關(guān)系，需要重新選取計算區(qū)域，建立局部三維模型。

三汊河河口閘三維建模計算范圍的選取，既要能夠體現(xiàn)河口閘上下游流態(tài)又要保證水流順直。根據(jù)流速分布曲線，3種工況差異大約出現(xiàn)在河口閘上游200 m位置，三維建模上游邊界取到閘上游約300 m處;同時考慮到三維模擬主要關(guān)注閘孔開閉對閘下游附近區(qū)域的影響，當(dāng)下游匯入長江后，長江流量大，影響因素多，邊界不易確定，因此下游取到外秦淮河河口附近區(qū)域。這個范圍既包含了重點關(guān)注的區(qū)域，又方便了邊界條件的確定。局部三維水動力模型計算區(qū)域及網(wǎng)格劃分如圖6所示。

采用正交曲線網(wǎng)格對計算區(qū)域進行網(wǎng)格剖分，對河口閘樞紐局部進行了網(wǎng)格加密，垂向上沿水深均勻分為10層，其中垂向網(wǎng)格以閘下游100 m處斷面為例。網(wǎng)格尺寸在 3～35 m之間，垂向模型網(wǎng)格見圖7。

由圖8（a）～（b）可以看出：工況2，3條件下（外秦淮河下游邊界水位3.5 m，上游來流量60 m3/s），計算范圍內(nèi)閘上游250 m到閘墩范圍內(nèi)流速普遍大于0.57 m/s，高于河道泥沙起動流速0.53 m/s，河道會發(fā)生沖刷;閘軸線到河口范圍內(nèi)水流流速小于0.38 m/s，小于起動流速0.47 m/s，不會發(fā)生沖刷。

由圖8（c）～（d）可以看出，工況5、工況6條件下（秦淮河下游邊界水位4.5 m，上游來流量60 m3/s），計算范圍內(nèi)閘上游250 m到閘墩范圍內(nèi)流速普遍小于0.48 m/s，低于河道泥沙的起動流速0.53 m/s，河道不會發(fā)生沖刷，僅在閘上游200 m附近區(qū)域流速大于泥沙起動流速，會發(fā)生沖刷;閘軸線到河口范圍內(nèi)水流流速小于0.24 m/s，小于起動流速0.53 m/s，不會發(fā)生沖刷。

3.3 研究范圍內(nèi)流速變化原因分析

本節(jié)從過水面積變化的角度簡要分析研究范圍內(nèi)流速大小變化的原因。

圖9為研究范圍內(nèi)典型斷面處過水面積示意圖（以工況1為例）。

從圖9可以看出：在閘上200 m處，河道斷面基本呈“V”形，過水面積較小，因此流速較大。在閘上200 m到閘軸線范圍內(nèi)，河道斷面形狀逐漸變?yōu)椤癠”形，過水面積增大，流速略有減小。在閘墩附近，單孔開啟工況下，過水?dāng)嗝娴膶挾冗_到最小，流速也較大。在河口閘下游，由于河道加寬，流速逐漸減小。而到了河口處，由于河道底面高程抬高，過水面積急劇減小，流速再次增大。故在河口閘下游流速呈現(xiàn)先減小后增大的趨勢，防沖槽處的流速較小。

4 結(jié) 論

本文構(gòu)建了三汊河河口閘大范圍水深平均模型和局部三維模型，對河口閘各種預(yù)設(shè)單孔開啟模式進行了水動力計算，分析了各種工況下的水流流速特征，確定了閘孔不同運行模式對流速分布的影響范圍。并通過將河道床面流速與計算所得泥沙起動流速進行比較，對單孔開啟模式下的河床沖淤情況進行了分析，得到了以下結(jié)論。

（1） 單孔開啟對流速分布的影響出現(xiàn)在河口閘上游100 m處至河口閘下游150 m處，越靠近閘軸線斷面影響越大，并在閘軸線斷面處達到最大，越遠離閘軸線斷面影響越小。

（2） 組合模式1工況下，河口閘上游段河道發(fā)生沖刷;組合模式2工況下，河口閘上游200 m處斷面發(fā)生沖刷。其余河道位置河道床面流速小于泥沙起動流速，不會發(fā)生沖刷。

（3） 通過典型斷面過水面積與計算區(qū)域內(nèi)河道流速變化對比，刻畫了計算區(qū)域內(nèi)河道過水面積與流速的關(guān)系，對河道流速變化的原因進行了分析，認(rèn)為河道過水面積的變化是導(dǎo)致流速變化的原因。

參考文獻：

[1] 河海大學(xué).秦淮河三汊河口閘下河口泥沙沖淤研究[R].南京：河海大學(xué)，2011.

[2] 胡孜軍，王玲玲，成高峰，等.三河閘泄流流態(tài)與進流糾偏數(shù)值研究[J].水利學(xué)報，2011，42（11）：1349-1354.

[3] 雷曉玲，袁廷，楊程，等.基于Delft3D模型的三峽航道環(huán)保疏浚水質(zhì)數(shù)值模擬研究[J].工業(yè)安全與環(huán)保，2016，42（3）：65-68.

[4] 左書華.Delft3D在鰲江口外平陽咀海域流場模擬中的應(yīng)用[J].水文，2007（6）：55-58.

[5] 廖庚強.基于Delft3D的柳河水動力與泥沙數(shù)值模擬研究[D].北京：清華大學(xué)，2013.

[6] 范翻平.基于Delft3D模型的鄱陽湖水動力模擬研究[D].南昌：江西師范大學(xué)，2010.

[7] 趙明登，曾威，唐小春，等.基于Arcgis與Delft 3D的渭河下游洪泛區(qū)洪水?dāng)?shù)值模擬[J].中國農(nóng)村水利水電，2012（10）：145-147，151.

[8] 郎韻，臧英平，唐洪武，等.三汊河河口閘泄流特性與整流措施[J].水利水電科技進展，2014，34（2）：11-15.

[9] 陳偉，臧英平，沃玉報.三汊河河口閘過閘水流流態(tài)改善措施數(shù)值分析[J].人民長江，2014，45（13）：14-16，21.

[10] 張紅武.泥沙起動流速的統(tǒng)一公式[J].水利學(xué)報，2012，43（12）：1387-1396.

[11] 韓其為.泥沙起動規(guī)律及起動流速[J].泥沙研究，1982（2）：11-26.

[12] 竇國仁.論泥沙起動流速[J].水利學(xué)報，1960（4）：44-60.

[13] 張瑞瑾，謝鑒衡，陳文彪.河流動力學(xué)[M].北京：中國工業(yè)出版社，1961.

[14] 唐存本.泥沙起動規(guī)律[J].水利學(xué)報，1963（2）：1-12.

[15] 沙玉清.泥沙運動學(xué)引論[M].北京：中國工業(yè)出版社，1965.

（編輯：胡旭東）

引用本文：

謝昕卓，陳偉，楊明強，等.南京市三汊河河口閘單孔開啟水流流態(tài)模擬研究

[J].人民長江，2021，52（7）：153-159，167.

Simulation on flow pattern of Sanchahe estuary sluice under single hole opening in Nanjing City

XIE Xinzhuo1，CHEN Wei2，YANG Mingqiang2，WANG Lu2，YANG Feifei2，CAI Xin1

（1.Collage of Mechanics and Materials，Hohai University，Nanjing 211100，China; 2.Nanjing Sanchahe River Gate Management Office，Nanjing 210036，China）

Abstract：

In order to simulate the flow pattern near the Sanchahe estuary sluice in Nanjing City and explore the feasibility of opening a single hole of the sluice，Delft 3D software was used to establish a large-scale water depth average model and a local three-dimensional mathematical model.Then the flow pattern near the Sanchahe estuary sluice was numerically simulated，the flow pattern under the most unfavorable water level condition was studied，and the influence range of different operation modes of the sluice on the velocity distribution was determined.Furthermore，by comparing the flow velocity at the bottom of the river with the sediment incipient velocity in the study area，the distribution of the deposition-erosion range of the river near the estuary gate was obtained.The research results can provide a basis for the study of single hole opening of the estuary sluice.

Key words：

operation mode of sluice;flow pattern;velocity distribution;Delft 3D model;river deposition-erosion;Sanchahe estuary gate