摘要： 為了揭示黃河下游典型概化心灘繞流近壁面流場(chǎng)的水力特性，以下游馬峪溝斷面附近典型心灘作為研究對(duì)象，利用模型試驗(yàn)與數(shù)值模擬探究不同流量與水深條件下概化心灘繞流近壁面流場(chǎng)的水流流態(tài)、河床高程、水流挾沙力及床面阻力等特性，并驗(yàn)證平面二維水沙模型計(jì)算得到的流速分布與河床高程的正確性.結(jié)果表明：分流區(qū)斷面受到灘頭影響，流速變小，并且流速最大值出現(xiàn)在汊道.匯流區(qū)斷面在灘尾影響下發(fā)生邊界層分離，水槽中央流速達(dá)到最小，此后流速達(dá)到最大并趨于平均化；采用平面二維水沙模型計(jì)算得到的流速分布與河床高程的模擬值與試驗(yàn)值基本吻合，表明采用數(shù)值模擬探究心灘繞流近壁面流場(chǎng)水力特性是可行的.本研究將為揭示黃河下游典型沖積河段心灘形態(tài)特征調(diào)整過(guò)程的動(dòng)力學(xué)機(jī)理提供理論基礎(chǔ).

關(guān)鍵詞： 心灘繞流；水力特性；數(shù)值模擬；黃河下游；水沙模型

中圖分類(lèi)號(hào)： S277.9文獻(xiàn)標(biāo)志碼： A文章編號(hào)： 1674-8530（2024）11-1120-08

DOI：10.3969/j.issn.1674-8530.22.0304

張春晉，任小鳳，張敏，等.黃河下游典型概化心灘繞流近壁面流場(chǎng)水力特性[J].排灌機(jī)械工程學(xué)報(bào)，2024，42（11）：1120-1127.

ZHANG Chunjin，REN Xiaofeng， ZHANG Min，et al.Hydraulic characteristics of flow field near wall of typical generalized midchannel bar in Lower Yellow River[J].Journal of drainage and irrigation machinery engineering（JDIME），2024，42（11）：1120-1127.（in Chinese）

Hydraulic characteristics of flow field near wall of typical

generalized midchannel bar in Lower Yellow River

ZHANG Chunjin1，2， REN Xiaofeng3， ZHANG Min1*， FENG Qing1， ZHAO Chongxu1

（1. Key Laboratory of Lower Yellow River Channel and Estuary Regulation Sediment， MWR， Yellow River Institute of Hydraulic Research， Zhengzhou， Henan 450003， China; 2. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering， Hohai University， Nanjing， Jiangsu 210098， China; 3． Henan Yellow River Hydrology Reconnaissance amp; Design Institute Co.， Ltd.， Zhengzhou， Henan 450004， China）

Abstract： The typical midchannel bar near the Mayugou crosssection in the Lower Yellow River was taken as the research object， the hydraulic characteristics of the flow field near the wall of the midchannel， such as the flow pattern， the riverbed elevation， the sediment carrying capacity and the bed surface resistance were described by model experiment and numerical simulation under different flow rate and water depth conditions， and the correctness of the flow velocity and the riverbed elevation calculated by the planar twodimensional watersediment model was verified. The results show that the flow velocity of the crosssection in the diverting area decreases under the influence of the midchannel bar head， and the maximum flow velocity all appears in the branch. Under the influence of the midchannel bar tail， the boundary layer separation occurs in the crosssection of the confluence area， and the flow velocity in the center of the flume reaches the minimum， and then the flow velocity gradually reaches the maximum and tends to average. The simulated results of the flow velocity and the riverbed elevation calculated by the planar twodimensional watersediment model are basically consistent with the experimental results， which indicates that it is feasible to use the numerical simulation method to explore the hydraulic characteristics of the flow field near the wall of the midchannel bar. This paper will provide theoretical basis and technical support for revealing the dynamic mechanism of the adjustment process of the midchannel bar morphological characteristics in the typical alluvial reach of the Lower Yellow River.

Key words： flow around midchannel bar;hydraulic characteristics;numerical simulation;Lower Yellow River;water and sediment model

心灘是河流地貌結(jié)構(gòu)的重要組成，其演變特點(diǎn)對(duì)于河流系統(tǒng)的完整性具有重要意義.心灘的存在使得河道由單股變?yōu)槎喙?，改變了水流的流?dòng)特性[1].近年來(lái)，許多專(zhuān)家對(duì)心灘繞流的水力特性開(kāi)展了大量的研究.于鵬杰等[2]與LOU等[3]以長(zhǎng)江下游典型多分汊河段東流水道為基礎(chǔ)，探究了不同來(lái)流條件下河道復(fù)雜的水流結(jié)構(gòu)特點(diǎn).XU[4]與LUCHI等[5]探討了游蕩性河道中心灘的發(fā)展過(guò)程，建立了河道心灘指標(biāo)與河道邊界條件、徑流量和輸沙量及水庫(kù)建設(shè)等控制變量之間的相互關(guān)系.余新明等[6]與金镠等[7]總結(jié)了心灘所在河段汊道的分流模式，探究了心灘的繞流結(jié)構(gòu)及其近底的輸沙特性.徐錫榮等[8]、張婕等[9]與KHAN等[10]分析了心灘演變特性對(duì)下游河勢(shì)的影響，闡明了當(dāng)心灘汊道呈現(xiàn)發(fā)展與衰退狀態(tài)時(shí)對(duì)下游分流分沙的影響，揭示了在心灘演變過(guò)程中水流的紊動(dòng)特性.已有研究成果主要闡述了天然河流分汊口與交匯口水力特性對(duì)河道平面形態(tài)的影響，揭示了河道調(diào)整過(guò)程的動(dòng)力學(xué)機(jī)理.

鑒于此，文中以黃河下游馬峪溝斷面附近心灘作為主要研究對(duì)象，對(duì)原型心灘模型概化，采用模型試驗(yàn)和數(shù)值模擬等研究方法分析不同流量和水深下心灘繞流近壁面流場(chǎng)水流流態(tài)、河床高程、水流挾沙力及床面阻力等水力特性，明確流量過(guò)程對(duì)心灘繞流水力特性的影響機(jī)制，并驗(yàn)證平面二維水沙模型的可靠性，為心灘附近流場(chǎng)特性研究提供一種切實(shí)可行的研究手段.

1模型試驗(yàn)

1.1模型比尺

文中采用了黃河下游馬峪溝斷面附近心灘作為研究對(duì)象，探究心灘繞流近壁面流場(chǎng)的水力特性，模型試驗(yàn)床沙級(jí)配與水流過(guò)程應(yīng)由下游床沙級(jí)配與水流過(guò)程所確定.馬峪溝斷面河寬約1 700 m，平均水深約5.0 m.根據(jù)重力相似準(zhǔn)則及試驗(yàn)場(chǎng)地，模型水平比尺設(shè)置為1∶850.由于原型河道寬深比可達(dá)340，若考慮為正態(tài)模型，則模型的水深偏小，無(wú)法滿(mǎn)足模型對(duì)水深需求.同時(shí)正態(tài)模型很難滿(mǎn)足水流阻力相似要求，在適應(yīng)水流運(yùn)動(dòng)與泥沙運(yùn)動(dòng)規(guī)律要求上靈活性小.為此，該模型設(shè)計(jì)為變態(tài)模型.考慮到模型水深滿(mǎn)足表面張力及測(cè)量需求，則模型垂向比尺設(shè)置為1∶42.針對(duì)模型幾何變率限制條件，模型相對(duì)保證率為0.8，表明模型80%以上過(guò)水?dāng)嗝媪鲌?chǎng)與原型相似.由于馬峪溝斷面附近心灘尺寸較大，外形不規(guī)則，為了便于分析心灘近壁面流場(chǎng)水力特性，對(duì)馬峪溝斷面附近心灘進(jìn)行概化處理.根據(jù)模型比尺，設(shè)置模型心灘長(zhǎng)度為2.00 m，寬度為1.00 m，高度為0.12 m.大流量時(shí)水流趨直存在水流流向與心灘中軸線(xiàn)夾角為0.因此將概化心灘中垂線(xiàn)與水流來(lái)流方向夾角設(shè)置為0.根據(jù)原型心灘高度，模型心灘沿長(zhǎng)度與寬度方向坡率分別設(shè)置為1∶8和1∶4.近3年，花園口以上河段比降均值為2.4×10-4.根據(jù)模型比尺，模型河床比降為4.8×10-3.基于鄭州熱電廠(chǎng)粉煤灰造價(jià)低、物理化學(xué)性能穩(wěn)定、選配加工方便，模型試驗(yàn)選用鄭州熱電廠(chǎng)粉煤灰作為模型沙樣.馬峪溝斷面在花園口水文站上游，距花園口水文站較近，選用花園口站懸沙與床沙中值粒徑代表原型懸沙與床沙級(jí)配，分別為0.025 mm和0.210 mm.

判斷模型床沙級(jí)配是否滿(mǎn)足泥沙起動(dòng)流速與揚(yáng)動(dòng)流速相似律依據(jù)是泥沙起動(dòng)流速比尺、揚(yáng)動(dòng)流速比尺及流速比尺達(dá)到接近.對(duì)于輕壤土，其級(jí)配與黃河下游的床沙級(jí)配較為接近.根據(jù)壤土及黏土不沖流速的表達(dá)式，當(dāng)水深在3～5 m時(shí)，原型床沙的起動(dòng)流速為0.84～0.91 m/s.當(dāng)鄭州熱電廠(chǎng)粉煤灰的中值粒徑為0.060 mm時(shí)，模型沙樣的起動(dòng)流速在0.11～0.14 m/s，泥沙起動(dòng)流速比尺為6.01～8.32，與流速比尺6.45接近，表明床沙滿(mǎn)足泥沙起動(dòng)相似條件.在天然河流部分懸沙與床沙存在交換現(xiàn)象，需復(fù)核懸沙落淤至床面后泥沙起動(dòng)流速是否符合相似條件.根據(jù)原型懸沙中值粒徑和水流黏性系數(shù)，當(dāng)水深在3～5 m時(shí)，原型懸沙落淤后起動(dòng)流速為1.06～1.16 m/s.對(duì)于已有模型沙樣，模型懸沙落淤至床面后起動(dòng)流速在0.14～0.16 m/s，泥沙起動(dòng)流速比尺為6.64～8.27，與流速比尺6.45接近，表明懸沙落淤床面滿(mǎn)足泥沙起動(dòng)流速相似.根據(jù)水槽試驗(yàn)結(jié)果[11]，與原型接近的天然沙揚(yáng)動(dòng)流速為起動(dòng)流速的1.5～1.7倍.當(dāng)水深為3～5 m時(shí)，原型床沙起動(dòng)流速為0.84～1.16 m/s，床沙揚(yáng)動(dòng)流速為1.38～1.91 m/s.由于模型相應(yīng)床沙揚(yáng)動(dòng)流速為0.20～0.22 m/s，則泥沙揚(yáng)動(dòng)比尺為6.32～9.57，與流速比尺6.45接近，表明模型所選床沙近似滿(mǎn)足泥沙揚(yáng)動(dòng)相似條件.因此，模型泥沙采用中值粒徑為0.060 mm的鄭州熱電廠(chǎng)粉煤灰作為床沙沙樣較為合理[11].

1.2試驗(yàn)系統(tǒng)

試驗(yàn)系統(tǒng)由電磁流量計(jì)、離心泵、進(jìn)水池、消能柵、矩形水槽、尾門(mén)、沉沙池、量測(cè)設(shè)備及其輔助裝置等組成，其中矩形水槽采用鋼化玻璃制作，長(zhǎng)度18.5 m，寬度2.0 m，高度1.0 m，坡降4.8×10-3.在槽底平鋪長(zhǎng)度15.0 m、寬度2.0 m、厚度0.3 m粉煤灰作為測(cè)試段，采用粉煤灰制作概化心灘，如圖1所示.流量采用矩形薄壁堰計(jì)算，尾門(mén)由翻板門(mén)結(jié)合小水閥控制，水面線(xiàn)由槽頂活動(dòng)測(cè)針采集，水位通過(guò)水位測(cè)針讀取，河床沖淤形變采用三維地形儀測(cè)量.流速采用剖面流速儀采集，采樣頻率50 Hz，測(cè)量時(shí)間60 s.

1.3試驗(yàn)方案

試驗(yàn)布置13個(gè)測(cè)試斷面，每個(gè)斷面布置17個(gè)測(cè)點(diǎn)，測(cè)點(diǎn)間距為0.125 m.鑒于黃河干支流水庫(kù)運(yùn)用后，年來(lái)沙量呈顯著下降，故在本試驗(yàn)不考慮來(lái)沙情況.模型試驗(yàn)共開(kāi)展了6組清水試驗(yàn)，并且每組試驗(yàn)持續(xù)5 h.模型試驗(yàn)的進(jìn)口邊界控制流量，而出口邊界控制水位，試驗(yàn)工況見(jiàn)表1，表中Qin為進(jìn)口流量，Zout為出口水位，hout為出口水深，X為濕周，R為水力半徑，hy為心灘淹沒(méi)深度.試驗(yàn)結(jié)束以后，測(cè)量不同流量與水深下心灘繞流近壁面流場(chǎng)水流流態(tài)、河床高程、水流挾沙力及床面阻力等水力特性.

1.4模型率定

文中還需率定概化模型，進(jìn)一步說(shuō)明概化模型設(shè)計(jì)的可靠性.率定結(jié)果如下：① 河段進(jìn)出口流速分布較為平順，可滿(mǎn)足模型設(shè)計(jì)進(jìn)出口條件；② 試驗(yàn)平均流速與實(shí)測(cè)流速較為接近，可研究實(shí)際河道水動(dòng)力條件；③ 模型試驗(yàn)灘頭壅水淤高，灘尾出現(xiàn)逆坡，該現(xiàn)象與下游實(shí)測(cè)結(jié)果吻合；④ 模型水流流態(tài)與原型一致，即水流流向心灘時(shí)受心灘壅阻，比降減小，流速降低，在心灘上游出現(xiàn)負(fù)比降，形成壅水.汊道水流在灘尾形成交匯，使流速減小，形成回流；⑤ 根據(jù)流速實(shí)測(cè)資料可計(jì)算分流比，該分流比與模型設(shè)計(jì)分流比較為吻合；⑥ 原型和模型床沙起動(dòng)流速比尺、揚(yáng)動(dòng)流速比尺均與流速比尺接近.綜上，該模型可更好地研究黃河下游心灘對(duì)河道水流結(jié)構(gòu)的影響.

2試驗(yàn)結(jié)果

2.1水流流態(tài)

圖2為心灘繞流的水流流態(tài)，由圖可知：分流區(qū)為心灘上游，分汊區(qū)為存在心灘的河道，而匯流區(qū)為心灘下游.由于心灘存在，水流受到阻力增大，迫使上游來(lái)流在心灘頭部點(diǎn)N形成汊道.水流流向心灘時(shí)受心灘壅阻，比降減小，流速降低，在心灘上游出現(xiàn)負(fù)比降，形成壅水.隨后部分水流繞過(guò)心灘點(diǎn)A和點(diǎn)B從汊道流向下游，水流在經(jīng)過(guò)灘體最寬位置后兩汊水流交匯，使得在灘頂尾部存在橫流，水深越淺，橫流越明顯.由于汊道水流流向不同，使得水流在交匯過(guò)程中耗散大量水流動(dòng)能，引起水流流速減小，形成回流區(qū)或緩流區(qū).

圖3為心灘繞流過(guò)程水面線(xiàn)H1分布，由圖可知：不同方案下心灘繞流過(guò)程水面線(xiàn)分布基本一致，水槽中間沿程水面線(xiàn)分布趨勢(shì)與汊道存在明顯不同，且中間水面線(xiàn)略低于汊道；灘頂前區(qū)水面線(xiàn)急劇下降，最低點(diǎn)位于汊道展寬8#斷面.灘頂后區(qū)至灘尾水面呈逆坡，在灘尾處水面線(xiàn)達(dá)到最大，且心灘下游河段水面線(xiàn)又回歸正常，表明水流所受阻力在汊道明顯大于單一河段；當(dāng)水深一定，水面線(xiàn)沿程趨勢(shì)基本一致，均呈先降低后升高，并隨流量增加，水面線(xiàn)分布將越高.當(dāng)流量一定，水面線(xiàn)沿程分布基本一致，并且隨著水深的增加，比降也隨之減小.

由圖4可知：當(dāng)水深一定時(shí)，隨流量增大各測(cè)點(diǎn)流速增大.同時(shí)當(dāng)流量一定時(shí)，隨水深增大各測(cè)點(diǎn)流速減小.從流速橫向分布看，5#斷面從左至右流速呈先減小后增大趨勢(shì)，在距左岸1.0 m附近受灘頭影響，流速變小，最大流速出現(xiàn)在汊道.9#斷面流速在橫向分布上呈先變小后變大趨勢(shì)，在匯流區(qū)流速達(dá)到最小，心灘完全被淹沒(méi)時(shí)斷面流速變化較小，心灘剛好被淹沒(méi)時(shí)斷面流速變化較大，這是由于在水深較小時(shí)水流流態(tài)相對(duì)紊亂；心灘高程起伏對(duì)于流速的影響較大，當(dāng)距離灘頭位置較遠(yuǎn)時(shí)，其水流的脈動(dòng)強(qiáng)度除個(gè)別測(cè)點(diǎn)波動(dòng)以外，整體沿著橫向分布均勻.進(jìn)入汊道以后，心灘的存在使得地形起伏較大，灘面的水深較小而流速較大，因此灘面水流的脈動(dòng)強(qiáng)度也相對(duì)較??；由于9#斷面處兩股水流發(fā)生交匯，水流紊動(dòng)強(qiáng)烈，使得水槽中間水流的脈動(dòng)強(qiáng)度較大.而在汊道出口位置水深相對(duì)較大，水流較為平穩(wěn)，水流脈動(dòng)強(qiáng)度也較小，因此表現(xiàn)出了“中間高和兩側(cè)低”的變化趨勢(shì).

2.2河床高程

圖5為河床高程H2分布，由圖可知：由于5#斷面距離分汊區(qū)較遠(yuǎn)，使得主流分流不明顯，斷面中間沖刷較深，而7#斷面位置主流已發(fā)生分流，因此沖刷體現(xiàn)在汊道；心灘灘體除尾部外均受到?jīng)_刷，高程降低，灘體降低幅度從灘頭到灘尾逐漸減小.而灘尾為水流回流區(qū)，流速較小，故從灘體沖刷挾帶泥沙在此落淤，造成灘體高程增加.雖然心灘尾部為淤積區(qū)，但汊道沖刷的泥沙不會(huì)在此淤積，淤積泥沙大多來(lái)自于灘頂；因分流區(qū)汊道較沒(méi)有分汊前流速明顯增大，使得挾沙能力變強(qiáng)，故汊道河床沖刷變形強(qiáng)烈，且出現(xiàn)深槽，高程降低；在相同流量條件下，水深越小，心灘淹沒(méi)程度越低，水流流速越大，挾沙能力越強(qiáng)，灘體及汊道沖淤變形越劇烈.在相同水深條件下，流量越大，灘體及汊道沖淤變形同樣越劇烈.

11#橫斷面因汊道水流在此交匯，流速減緩，泥沙在此沉積落淤，使得斷面局部高程增加，而斷面兩側(cè)也在汊道水流沖刷下繼續(xù)降低高程而形成深槽；在灘體淹沒(méi)程度相同時(shí)，流量越大，汊道深槽越明顯，深槽高程越低，而中間淤積區(qū)淤積高程反而越高.當(dāng)流量相同時(shí)，淹沒(méi)程度越低，河床變形越明顯，中間淤積厚度越高，汊道沖刷深度越大.由于是清水試驗(yàn)，心灘變化趨勢(shì)表現(xiàn)為長(zhǎng)度增大，而高度減小.心灘長(zhǎng)度增加是因?yàn)樾臑┪膊恳騾R流引起泥沙在匯流區(qū)沉積落淤，灘尾向下游延伸.灘體長(zhǎng)度與灘尾匯流區(qū)位置有關(guān)，而灘體長(zhǎng)度則與沖刷后心灘高度成反比.

2.3水流挾沙力

圖6為水流挾沙力F分布，由圖可知：進(jìn)入汊道前水流挾沙力較小，各斷面水流挾沙力相差不大.進(jìn)入汊道后水流挾沙力變大，在心灘最寬處達(dá)到最大，經(jīng)過(guò)灘體最寬處，汊道水流挾沙力變小.最后水流通過(guò)汊道在9#斷面匯合，水流挾沙力又趨于平穩(wěn)；當(dāng)水深相同時(shí)，隨著流量的增大，則水流挾沙力增大.而當(dāng)流量相同時(shí)，隨著水深的減小，水流挾沙力則表現(xiàn)為增大.

2.4床面阻力

圖7為床面阻力沿程分布，圖中σ為床面切應(yīng)力，由圖可知：分汊區(qū)床面阻力明顯大于單一河段床面阻力，從灘頭5#斷面開(kāi)始增大到灘體最寬處7#斷面達(dá)到最大，經(jīng)過(guò)7#斷面后到灘尾9#斷面，床面阻力又逐漸減??；同一水深下，隨著流量增大床面阻力越大.方案5床面阻力最大，其次為方案3，最小為方案1；同一流量下，隨著水深減小床面阻力將增大.方案2床面阻力最大，其次為方案3，最小為方案4.

2.5結(jié)果討論

已有水流挾沙力公式都可反映水流挾沙力由床面阻力、流速及水深決定.若利用已有公式分析各因素與水流挾沙力的關(guān)系似有不妥，應(yīng)從水流挾沙力本質(zhì)出發(fā)分析各因素與水流挾沙力關(guān)系.床面阻力指河床對(duì)水流運(yùn)動(dòng)阻礙作用，并以熱能形式消耗.床面阻力對(duì)能量消耗體現(xiàn)在水流紊動(dòng)能，部分紊動(dòng)能又轉(zhuǎn)化為泥沙懸浮功.當(dāng)水深、流速及比降相同，床面阻力越大，轉(zhuǎn)化為泥沙懸浮功越多，越能促進(jìn)泥沙懸浮，提高水流挾沙力.若上述因素變化較大，則床面阻力與水流挾沙力很難獲得定性關(guān)系.本模型試驗(yàn)用于分析清水狀態(tài)下心灘對(duì)河道水流結(jié)構(gòu)影響，為此以馬峪溝斷面附近典型心灘作為研究對(duì)象，并對(duì)原型心灘模型概化.模型概化是一種有效研究手段，不僅可簡(jiǎn)化研究對(duì)象使研究結(jié)論更具普適性，還可對(duì)原型河道起指導(dǎo)作用.需強(qiáng)調(diào)概化模型僅采用黃河下游心灘附近河道邊界條件，至于河道邊界形態(tài)并未完全考慮.概化模型根據(jù)流速比尺換算得到流速與實(shí)際流速較為接近，因此該模型可研究實(shí)際河道水流結(jié)構(gòu).余新明等[6]開(kāi)展了分汊河道水沙特征研究，其結(jié)論是進(jìn)口無(wú)泥沙補(bǔ)給，分流區(qū)及上下游會(huì)形成沖刷深槽；在匯流區(qū)特別是灘尾，由于其能量損失使得泥沙在該處落淤，上述現(xiàn)象與文中結(jié)論一致.就心灘存在對(duì)水流結(jié)構(gòu)影響而言，本研究能反映原型河道真實(shí)流態(tài)，對(duì)原型水流結(jié)構(gòu)探究具有指導(dǎo)意義.

3數(shù)值模擬

3.1數(shù)學(xué)模型

平面二維水沙模型的定解條件包括初始條件與邊界條件[12].幾何模型邊界條件涉及進(jìn)、出口及水槽兩側(cè)，其中進(jìn)口采取流量邊界，出口采取水位邊界，水槽兩側(cè)采取無(wú)滑移邊界.在幾何模型中，床沙中值粒徑0.060 mm，床沙鋪設(shè)厚度0.3 m.為保證模型計(jì)算效率與求解精度，幾何模型流體域采用非結(jié)構(gòu)化四面體網(wǎng)格加密.由0.02 m與0.03 m這2種體網(wǎng)格計(jì)算的幾何模型出口斷面平均流速的相對(duì)誤差不超過(guò)1%.一般定義相鄰2種體網(wǎng)格尺寸劃分的幾何模型流體域，當(dāng)出口斷面平均流速的相對(duì)誤差小于1%時(shí)，則認(rèn)為體網(wǎng)格尺寸對(duì)于幾何模型出口斷面平均流速的影響程度可忽略不計(jì).因此，文中將幾何模型流體域的體網(wǎng)格尺寸設(shè)置為0.02 m，可有效地滿(mǎn)足網(wǎng)格無(wú)關(guān)性檢驗(yàn)的要求.

為了更好地捕捉近壁面的黏性底層，在兩側(cè)固壁采取了局部邊界層加密.文中根據(jù)尼庫(kù)拉茲經(jīng)驗(yàn)公式計(jì)算了不同雷諾數(shù)條件下沿程阻力系數(shù)，并求解水槽壁面的黏性底層厚度.根據(jù)黏性底層厚度計(jì)算值0.27 mm，第1層網(wǎng)格厚度設(shè)置為0.15 mm（y+=30），網(wǎng)格層間比例因子設(shè)置為1.1，且邊界層區(qū)域的網(wǎng)格層數(shù)設(shè)置為11.采用ADI有限差分法離散基本控制方程組，水流連續(xù)性方程中水位對(duì)于時(shí)間偏導(dǎo)數(shù)項(xiàng)采用向前差分格式，水流運(yùn)動(dòng)方程中對(duì)流項(xiàng)則采用一階迎風(fēng)格式，其余各項(xiàng)均采用中心差分格式[13-16].

3.2結(jié)果驗(yàn)證

圖8為方案3的流速分布與河床高程模擬值與試驗(yàn)值對(duì)比，從圖中可以看出：模擬值與試驗(yàn)值吻合，誤差不超過(guò)6%；隨著灘頂水深增加，旋槳式流速儀對(duì)水流干擾越小，模擬值與試驗(yàn)值吻合越理想；5#斷面水流在灘頭影響下流速降低，而灘尾9#斷面由于邊界層分離導(dǎo)致流速降低，引起灘尾形成局部回流；在灘頭5#斷面，由于灘頭存在改變河道原有水流結(jié)構(gòu)，引起心灘附近流速重新分配，增大水流紊動(dòng)，引起河床沖刷.而11#斷面在灘尾回流影響下泥沙在該區(qū)域落淤，導(dǎo)致心灘向下游逐步延長(zhǎng).

4結(jié)論

1） 心灘繞流過(guò)程中灘頂前部水面線(xiàn)急劇下降，最低點(diǎn)位于展寬斷面，灘頂后部至灘尾區(qū)域水面線(xiàn)逐漸升高，在灘尾水面線(xiàn)達(dá)到最大，且心灘下游河段水面線(xiàn)又回歸正常.

2） 在分流區(qū)，從左至右流速呈先減小后增大，最大流速出現(xiàn)在汊道.匯流區(qū)受灘尾影響，流速達(dá)到最小，此后流速開(kāi)始逐漸增大.心灘的灘體除尾部外均受到?jīng)_刷，高程降低，且灘體降低幅度從灘頭到灘尾逐漸減小.而灘尾為水流匯流區(qū)，流速較小，灘體沖刷挾帶泥沙在此落淤，造成灘體高程增加.

3） 采用平面二維水沙模型計(jì)算得到的流速分布與河床高程模擬值與試驗(yàn)值基本吻合，表明采用數(shù)值模擬探究心灘繞流近壁面流場(chǎng)水力特性是可行的，將為下游心灘研究提出一種切實(shí)可行的方法.

參考文獻(xiàn)（References）

[1]KHAN M A， SHARMA N. Investigation of coherent flow turbulence in the proximity of midchannel bar [J]. KSCE journal of civil engineering， 2019， 23（12）： 5098-5108.

[2]于鵬杰，文杰，談廣鳴，等.多分汊河型水流結(jié)構(gòu)概化模型試驗(yàn)研究[J].武漢大學(xué)學(xué)報(bào)（工學(xué)版）， 2020， 53（5）： 385-393.

YU Pengjie， WEN Jie， TAN Guangming， et al. Experimental study on generalized model of braided river flow structure [J]. Engineering journal of Wuhan University， 2020， 53（5）： 385-393. （in Chinese）

[3]LOU Y Y， MEI X F， DAI Z J. et al. Evolution of the midchannel bars in the middle and lower reaches of the Changjiang （Yangtze） River from 1989 to 2014 based on the Landsat satellite images： impact of the Three Gorges Dam [J]. Environmental earth sciences， 2018， 77： 1-18.

[4]XU J X. Evolution of midchannel bars in a braided river and complex response to reservoir construction： an example from the middle Hanjiang River， China [J]. Earth surface processes amp; landforms，1997， 22： 953-965.

[5]LUCHI R， ZOLEZZI G， TUBINO M. Modelling midchannel bars in meandering channels [J]. Earth surface processes amp; landforms，2010， 35： 902-917.

[6]余新明， 談廣鳴， 張悅， 等. 分汊河道水沙輸移特征試驗(yàn)[J].武漢大學(xué)學(xué)報(bào)（工學(xué)版）， 2007，40（4）： 9-12.

YU Xinming，TAN Guangming，ZHANG Yue，et al.Experimental on transportation characteristics of flow and sediment in distributaries channel [J]. Engineering journal of Wuhan University， 2007，40（4）： 9-12. （in Chinese）

[7]金镠， 樓飛. 長(zhǎng)江口深水航道分流比與汊道穩(wěn)定性評(píng)估[J].水運(yùn)工程， 2021（2）： 53-57.

JIN Liu， LOU Fei. Assessment of water diversion ratio and channel stability in the Yangtze estuary deepwater channel [J]. Port amp; waterway engineering， 2021（2）： 53-57. （in Chinese）

[8]徐錫榮， 鐘凱， 白金霞.長(zhǎng)江小黃洲演變對(duì)下游汊道分流特性的影響[J].河海大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2014， 42（3）： 211-216.

XU Xirong， ZHONG Kai， BAI Jinxia. Influence of evolution of Xiaohuangzhou bifurcated reach on flowdividing characteristics of downstream bifurcated reaches of Yangtze River[J]. Journal of Hohai University（natural sciences）， 2014， 42（3）： 211-216. （in Chinese）

[9]張婕， 王平義， 胡杰龍， 等. 分汊河道近底水流特性試驗(yàn)研究[J]. 科學(xué)技術(shù)與工程， 2021， 21（8）： 3297-3303.

ZHANG Jie， WANG Pingyi， HU Jielong， et al. Experimental study on the characteristics of nearbottom flow in bifurcated channel [J]. Science technology and engineering， 2021， 21（8）： 3297-3303. （in Chinese）

[10]KHAN M A， SHARMA N， PU J H， et al. Midchannel braidbarinduced turbulent bursts： analysis using octant events approach [J]. Water， 2022， 14（3）：14030450.

[11]劉宇，辜玉慧，劉炎鑫，等.涇河流域水文氣象演變及其變化歸因[J].水利水電技術(shù)，2023，54（10）：34-48.

LIU Yu，GU Yuhui，LIU Yanxin，et al.Hydrometeorological evolution and its change attribution in Jinghe River Basin[J].Water resources and hydropower engineering，2023，54（10）：34-48.（in Chinese）

[12]李志祥，徐輝，馮建剛，等.虹吸式出水管虹吸形成過(guò)程的壁面壓力及脈動(dòng)特性[J].水利水電科技進(jìn)展，2023，43（1）：43-49.

LI Zhixiang，XU Hui，F(xiàn)ENG Jiangang，et al.Wall pressure and fluctuating characteristics of a siphon outlet pipe during siphoning formation process[J].Advances in science and technology of water resources，2023，43（1）：43-49.（in Chinese）

[13]袁淑霞，裴龍瑤，張宇鵬，等.非牛頓流體灘涂鋼板樁圍堰泄漏量分析[J].水利水電科技進(jìn)展，2023，43（4）：67-72.

YUAN Shuxia，PEI Longyao，ZHANG Yupeng，et al. Analysis on leakage of steel sheet cofferdam in mud flat of nonNewtonian flow[J].Advances in science and technology of water resources，2023，43（4）：67-72.（in Chinese）

[14]張雷，曹尊毅，王金亮，等，運(yùn)行工況對(duì)多沙河流水輪機(jī)導(dǎo)葉區(qū)磨損影響研究[J].水利水電技術(shù)，2022，53（4）：148-156.

ZHANG Lei， CAO Zunyi， WANG Jinliang， et al. Study on influence of operating conditions on turbine vane wear in sandy river[J]. Water resources and hydropower engineering，2022，53（4）：148-156.（in Chinese）

[15]趙占超，侯素珍，賴(lài)瑞勛，等，黃河內(nèi)蒙古河段水沙及河床演變特性變化分析[J].水利水電技術(shù)，2022，53（2）：109-120.

ZHAO Zhanchao， HOU Suzhen， LAl Ruixun， et al， Analysis on changes of watersediment and riverbed evolution characteristics in Inner Mongolia reach of Yellow River[J]. Water resources and hydropower engineering， 2022，53（2）： 109-120.（in Chinese）

[16]張春晉， 孫西歡， 李永業(yè). 龍?zhí)ь^式泄洪洞安全監(jiān)測(cè)指標(biāo)分析[J]. 排灌機(jī)械工程學(xué)報(bào)， 2022， 40（11）： 1112-1119.

ZHANG Chunjin， SUN Xihuan， LI Yongye. Analysis of safety monitoring indicators of dragonhead spillway tunnel[J]. Journal of drainage and irrigation machinery engineering， 2022， 40（11）： 1112-1119.（in Chinese）

（責(zé)任編輯談國(guó)鵬）

收稿日期： 2022-12-12； 修回日期： 2023-04-30； 網(wǎng)絡(luò)出版時(shí)間： 2024-11-08

網(wǎng)絡(luò)出版地址： https：//link.cnki.net/urlid/32.1814.TH.20241108.0950.020

基金項(xiàng)目： 黃科院發(fā)展基金資助項(xiàng)目（202006）；黃河水科學(xué)研究聯(lián)合基金資助項(xiàng)目（U2243220）；水利部泥沙科學(xué)與北方河流治理重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金資助項(xiàng)目（IWHR-SEDI-2022-07）

第一作者簡(jiǎn)介： 張春晉（1989—），男，山西平遙人，高級(jí)工程師，博士后（zhangchunjintyut@163.com），主要從事河床演變及河道整治研究.

通信作者簡(jiǎn)介： 張敏（1979—），女，陜西大荔人，正高級(jí)工程師，博士（zmiii@163.com），主要從事河床演變及河道整治研究.