平面淹沒射流垂直沖刷泥沙試驗研究

顧磊　倪雁　黃佳麗　倪福生

摘要：射流在水利工程中作為沖沙手段已有較廣泛的應用，沖刷過程的射流沖坑特性與泥沙懸浮特性均是衡量沖刷效果的重要指標。通過改變射流速度、噴射靶距、泥沙粒徑，進行了平面射流沖刷沙床的系列試驗。結果表明：沖坑與泥沙的懸浮尺寸隨時間均呈現快速增長期、緩慢發展期和動態穩定期3個發展階段，在此過程中，泥沙懸浮尺寸會出現峰值，最終在某一數值上下波動;總體上沖坑深度與泥沙懸浮尺寸隨射流速度和噴射靶距的增大而基本呈線性增大，短時間（1 s）內靶距的增加對沖坑深度的影響較小;沖坑深度與泥沙懸浮尺寸隨粒徑的增大有所減小;沖刷參數可以用來綜合衡量射流的沖刷能力，隨著射流沖刷參數的增大，無量綱化的沖坑深度、懸沙高度和懸沙寬度均呈線性增長趨勢。

關鍵詞：沖刷參數;沖坑尺寸;泥沙懸浮尺寸;射流沖刷;泥沙;清淤工程

中圖分類號：TV142;TV697

文獻標志碼：A

doi：10. 3969/j .issn.1000- 1379.2019.06.009

在河流、航道、湖泊和水庫等處進行的清淤工程中，射流清淤方式日益得到重視。該方式的基本原理是利用射流沖擊沙床，增強水流局部紊動強度，使泥沙懸浮形成異重流，結合水流作用產生向下游的泥沙運動，總體上提高水流的挾沙能力，達到消除泥沙淤積的目的。近年來，射流沖刷方式已得到較廣泛的應用[1-2]，并獲得了顯著成效。

現有關于射流沖蝕的研究中，Rajaratnam N等[3]通過試驗測量了動態沖刷深度的相關規律：齊梅蘭等[4]基于泥沙起動理論和平面射流運動規律，建立了射流沖刷平衡深度計算式，并采用不同泥沙粒徑和射流參數的沖刷試驗數據，驗證了計算式的有效性：槐文信等[5]基于Fluent軟件對二維垂向射流沖刷沙床進行了數值模擬研究，提出了沖刷平衡時沖刷坑深度、寬度等特征長度的半經驗公式：張浩等[6]針對粗沙沙床條件下二維垂向淹沒射流沖坑深度進行了試驗研究，探索了射流速度對動態沙坑深度的影響：申振等[7]基于泥沙動力學和射流沖刷理論，在不同靶距下進行了針對粗沙和中沙沙床的淹沒射流沖刷試驗，分析了沖刷開始10 s內坑深的發展規律：顧磊等[8]探究了雙股射流沖刷時噴嘴間距對沖刷坑形狀、沖坑深度和沖坑截面積的影響：劉思源[9]基于兩相流模型和湍流模型，對不同粒徑的沙床進行了二維垂向淹沒射流沖刷的數值模擬，研究了粒徑和射流速度對沙床沖刷坑深的影響。可以看到，現有研究大多單獨將射流沖坑的相關參數作為衡量沖刷效果的關鍵指標。

實際上，射流沖坑僅反映了射流在沙床上形成的沖刷效果，射流清淤的目的不止于此，還需通過異重流和上游水流輸送的方式將泥沙挾帶至下游或泄沙口[10]。泥沙懸浮高度越高泥沙輸移距離越遠，懸浮范圍越廣輸沙量越大，而通過懸沙懸浮高度與寬度即可評估懸沙范圍。因此，泥沙懸浮特性同樣應作為衡量射流沖刷效果的重要指標。

本文將沖坑特性和泥沙懸浮特性作為射流沖刷效果的衡量指標，開展了平面淹沒射流沖刷泥沙的系列試驗，研究了關鍵參數對沖刷效果的影響，以期為射流清淤工程提供一定參考。

1 試驗裝置及方法

本試驗所用裝置見圖1。所用水槽由透明有機玻璃板構成，長1 300 mm，寬300 mm，高1 000 mm，兩側900 mm處設有溢流擋板。平面噴嘴出口寬度為2.5mm，采用精加工工藝保證出口寬度誤差在+0.005 mm內，平面噴嘴長邊尺寸為290 mm，遠大于噴嘴寬度，以保證射流的二維平面特征。噴嘴由夾具固定，松動夾具即可使噴嘴上下活動以獲得不同的噴射靶距。水流動力由一臺離心水泵提供，水泵由變頻電機驅動，調節電機頻率即可改變水泵轉速，由水泵出口流出的水流經由一個調節閥門后進人平面噴嘴，改變水泵轉速或調節閥門開度均可獲得不同的射流速度。射流速度由電磁流量計測得流量后計算得到，該流量計精度可達0.2%。整個裝置位于一個大水池內，當試驗水槽內水由兩側溢流擋板溢出時，由大水池內集水坑收集溢出水并與水泵人口通過管道相連，從而實現水的循環使用。

試驗時可通過透明玻璃觀察平面噴嘴的動態沖刷過程，在其中一面有機玻璃上貼有透明的標準坐標紙，以此作為參照，可以讀出沖刷時的泥沙懸浮高度和懸浮寬度;在透明玻璃外側放置直尺，以此作為參照，可以讀出射流沖坑深度。在正對該面前方，放置一頻率最高可達100 fps（即每秒拍攝100幀圖像）的攝像機，對沖刷過程進行實時拍攝。

本研究主要對射流沖刷動態沖坑以及懸浮泥沙的發展進行分析，重點關注沖坑在剛開始時短時間內的發展過程，并對不同沙床條件、不同射流靶距、不同射流速度下的動態沖坑以及懸浮泥沙進行對比分析。當沖坑和懸浮泥沙的變化均趨于穩定時結束試驗。通過視頻處理軟件對拍攝的視頻進行逐幀分析，讀出各時刻的沖坑深度、泥沙懸浮高度和懸浮寬度。

為了很好地模擬實際疏浚工程施工過程中復雜的水下河床條件，分別選取中值粒徑為0.5、1.8 mm的兩種典型沙床條件進行射流試驗，通過不同沙質沙床的對比試驗，對實際疏浚施工中水下射流沖刷有更深入的了解。

射流沖刷試驗的主要影響因素為沙床要素和射流要素。為了更深入地研究射流的影響要素和沖刷機理，選取不同射流要素的試驗結果進行比較，試驗工況見表1。分別進行同一噴射靶距下不同射流速度和同一射流速度下不同噴射靶距的試驗，受相機拍攝范圍和試驗空間的限制，射流速度較大時無法同時拍攝到完整的沖坑和懸沙形態，因此在結果分析時選擇合適的工況分別進行沖坑和懸沙數據的提取。

2 試驗結果與分析

2.1 沖刷形貌

通過各工況下的試驗視頻，觀察沖坑形貌以及懸沙形貌的發展過程，發現沖坑形貌以及懸沙形貌隨時間的發展規律基本相同。

圖2為射流速度為3 m/s、射流靶距為150 mm條件下粗沙河床沖坑的發展過程，可以看出沖坑基本呈初始、發展和平衡3個變化階段。在初始階段，床面沙粒在射流作用下起動，主要以推移質方式運動，伴隨著射流流線的偏轉而被全部帶出沖坑并在兩側迅速堆積形成沙丘;在發展階段，射流與更多的沙粒發生作用，射流中軸線處沙粒運動最為劇烈，一部分沙粒以推移質運動方式沿沙坡上爬，左右沙丘不斷向兩側移動，沖坑寬度越來越大，沖坑中明顯出現兩個上下擺動的旋渦[11]，另一部分沙粒則以懸移質形式在沖坑中隨旋渦旋轉，沖坑形態逐漸清晰;平衡階段，大多數泥沙顆粒已經被甩出沖坑，隨旋渦一起旋轉的泥沙顆粒減少，直至基本恒定數量的沙粒隨旋渦運動，此時沖刷速率已不明顯，沖坑形貌基本不再變化，沖坑趨于穩定。

圖3為射流速度為4 m/s、噴射靶距為30 mm時粗沙的懸浮發展過程。可以看到，當射流剛剛到達床面時（0 s），泥沙開始起動，被兩側回流挾帶而懸浮，此時兩側懸起形態基本對稱。隨著時間推移（0.4-5 s），懸浮泥沙逐漸上升、擴展，懸浮泥沙逐漸增多，沙坑逐漸增大，在沖坑內出現明顯旋渦，兩側泥沙懸浮形態呈非對稱性，有時左側懸沙增多，有時右側懸沙增多。某一時刻（5 s）后，沖坑基本達到平衡，懸沙范圍基本動態穩定，只是懸沙在渦流挾帶下在兩側交替增多，左右側懸沙峰值呈現一定的周期性。為進一步獲取沖坑深度和泥沙懸浮尺寸的詳細發展情況，測量各時刻的沖坑深度、泥沙懸浮高度和懸浮寬度，其隨時間的變化見圖4、圖5。

圖4為射流速度為3 m/s、射流靶距為150 mm工況下沖坑深度隨時間的變化情況，其他工況與其基本相似。其中圖4（a）為10 s時間內沖坑深度的發展過程，可以看出在前3s坑深急劇變大，6s后坑深已基本無變化.3-6 s坑深變化緩慢，整體呈指數增長趨勢。在一些沖淤工程中，移動噴嘴與沙床上某點的作用時間較短，由1 s時間內沖坑深度的發展情況（見圖4（b））可知，1 s內沖坑深度呈線性增長，增長斜率較大，可見該時段射流與泥沙處于快速交互作用階段。

圖5以射流速度為4 m/s、噴射靶距為30 mm工況為例顯示了泥沙懸浮尺寸的變化過程。由圖5可知，在沖刷開始1 s內，與沖坑深度相似，泥沙懸浮高度與懸浮寬度均呈線性快速增長。1-4 s內，其數值盡管有所波動，但整體上仍然呈增長趨勢。在4s左右泥沙懸浮高度與懸浮寬度均出現最大值。4s后泥沙懸浮尺寸有所下降，最后達到動態平衡并在某一數值上下波動，波動頻率在Is以內，這可能是紊動射流的非定常性所決定的。

綜上所述，無論是沖坑尺寸還是泥沙懸浮尺寸，均經歷了快速增長期、緩慢發展區和動態穩定期。在射流作用初期，沖刷尺寸均呈線性快速增長，泥沙懸浮尺寸會出現最大值。動態平衡時，沖坑尺寸基本穩定，而泥沙懸浮尺寸則由極大值有所下降，并因射流動力的非定常性而呈現周期波動。

2.2 射流速度的影響

分別將沖坑深度和泥沙懸浮尺寸隨射流速度的變化繪制在圖6、圖7中。圖中各數值點的噴射靶距、泥沙粒徑等條件均保持一致，僅改變射流速度，取值為各工況下的最大沖坑深度和泥沙懸浮尺寸。為反映作用時間的影響，分別選取10 s（見圖6（a））和1 s（見圖6（b））時間內的沖坑深度值。由圖6可知，無論是10s時段還是1 s時段，沖坑深度均隨射流速度線性增大，且兩個時段內該增大的趨勢線基本平行。可見增大射流速度對加深沖坑效果顯著，而且無論作用時間長短，其最大沖深的增長趨勢都由射流對床面作用初期的沖擊效果決定。發展初期的沖刷速率直接決定了所能達到的最大沖坑深度。

由圖7可知，泥沙懸浮尺寸隨射流速度的變化與沖坑深度的相似，均隨射流速度的增大而線性增大。可見射流速度對泥沙懸浮尺寸的影響也非常顯著。

實際上，無論是沖坑深度還是泥沙懸浮尺寸，其形成的本質為泥沙的起動和懸浮，這取決于作用于床面泥沙的流速或切應力以及反沖水流速度。當其他條件恒定時，射流速度正是上述條件的綜合反映。射流速度增大時，到達床面的流速與床面切應力均有所增大，更多的泥沙被挾帶懸浮，而反沖水流的速度與旋渦能量也有所增大，就會為懸浮泥沙拋離沖坑或保持懸浮提供更多動力，因此沖坑與泥沙懸浮尺寸都會有所增大。

2.3 噴射靶距的影響

圖8和圖9分別為沖坑尺寸和泥沙懸浮尺寸隨噴射靶距的變化，其中圖8分別選取了10 s時段內和Is時段內的最大沖坑深度。可以看出，在本文試驗條件下，較長時段內的沖坑深度、泥沙懸浮高度和懸浮寬度均隨噴射靶距的增大而線性增大。由射流的水下衰減規律可知，隨著截面與射流出口距離的增加，射流中心速度逐漸減小[12]，那么靶距減小時到達床面的水流動力應有所增大，但試驗結果卻顯示，小靶距下的沖坑深度和泥沙懸浮尺寸均比大靶距下小得多，原因可能是射流不僅為泥沙起動提供了動力，而且方向向下的射流對泥沙的向上運動還具有一定的抑制作用，尤其是小靶距時射流還未充分發展即與泥沙發生作用，核心區的強剪切大大阻礙了中線附近泥沙的起動。此時，從相應的試驗視頻也可以看到，整個沖刷流場明顯受到較強約束，其紊動非常平緩。這種抑制作用會隨噴射靶距的增加而衰減，其衰減速度比射流動力衰減的速度快得多。這一原因還未得到證實，需測量更多詳細信息，深入研究兩相流場進行驗證。

此外，從圖8（b）可知，當射流作用時間較短時，噴射靶距對沖坑深度的影響較小，在噴射靶距為150 mm時沖坑深度略有升高，但幅度非常小。

可以肯定的是，當噴射靶距增大到足夠大時，必然會出現床面射流速度小于泥沙起動速度的情況，此時無泥沙懸浮，沖坑和懸沙尺寸都將為零。因此，必然存在一個噴射靶距使得沖刷尺寸達到最大。清淤工程中應設法找到該最佳噴射靶距。當噴嘴移動沖刷時，由于地形的波動，因此噴射靶距會有所變化，但這種變化在較短作用時間內對沖刷效果的影響并不顯著。

2.4 泥沙粒徑的影響

選擇中值粒徑為0.5 mm的中沙和1.8 mm的粗沙兩種典型泥沙進行分析，其在相同射流條件下達到動態平衡時的沖刷形態見圖10。由圖10可知，中沙的沖坑深度、泥沙懸浮高度和寬度分別為為148、190、410 mm.粗沙的沖坑深度、泥沙懸浮高度和寬度分別為98、190、240 mm。可見當粒徑增大時，沖坑深度與泥沙懸浮寬度均有所減小。這由泥沙的起動原理[13]可以很好地解釋：當粒徑增大時，顆粒慣性增大，其起動條件增加，懸浮所需動力也相應增大，對射流能量的消耗也就有所增大。顯然，當面對粗顆粒泥沙時，為保證相同的沖刷結果，需提供更強的射流動力。

2.5 隨沖刷參數的變化

Rajaratnam N等提出了圓柱噴嘴的沖刷參數E。，并獲得了較為準確的經驗公式[14]：式中：U0為射流速度;d為噴嘴直徑;h為噴射靶距;g為重力加速度;D為泥沙粒徑;△p為泥沙密度與水的密度之差;p為水的密度。

參數E采用無量綱數反映了射流在床面處的強度與泥沙的水下慣性之比，沖坑特性與泥沙懸浮特性均與這兩方面因素相關。本文仿照Rajaratnam N等的做法，結合平面射流的衰減規律，對圓柱噴嘴的沖刷參數稍作變化，見式（2）： 式中：2b為平面噴嘴的寬度。

采用式（2）計算了各工況的沖刷參數，并用噴射靶距對沖坑深度、泥沙懸浮高度和懸浮寬度的最大值進行無量綱化，獲得的無量綱沖刷尺寸與沖刷參數的關系見圖11。

由圖11可知，隨著射流沖刷參數E的增大，無量綱化的沖坑深度、泥沙懸浮高度和懸浮寬度均呈線性增長趨勢。線性擬合關系式分別為

3 結語

分別在中值粒徑為0.5、1.8 mm的沙床上，改變射流速度和噴射靶距，進行了平面淹沒射流的垂直沖刷系列試驗，觀測了沖坑、懸沙的形態和尺寸發展，在本文試驗條件范圍內結論如下。

（1）沖刷主要經歷快速增長期、緩慢發展期和動態平衡期3個階段，其中快速增長期的沖刷速率直接決定了整體沖刷趨勢。沖坑深度基本呈指數增長趨勢，而泥沙懸浮尺寸會出現最大值，然后有所下降并維持在某一數值上下波動。因此，工程中若能在快速增長期適當延長射流的作用時間，則有利于提高沖沙量。

（2）長時段沖刷的沖坑深度和懸沙尺寸隨射流速度和沖刷靶距的增大呈線性增大，隨泥沙粒徑的增大而有所減小。當射流作用時間較短時，噴射靶距的變化與長時段內的作用規律不同，當噴射靶距變化時，沖坑尺寸變化不明顯。工程中應選擇合適的噴射靶距，而對于運動噴嘴，沖刷效果對噴射靶距的敏感性則會大大減弱。

（3）平面噴嘴的沖刷參數可以綜合反映射流相對泥沙的沖刷強度，采用該參數能夠較好地預測平面射流沖刷的沖坑尺寸和泥沙懸浮尺寸。

參考文獻：

[1] 張原鋒，張隆榮，射流沖刷試驗及其在黃河潼關河段的應用[J].水利水電技術，2004，35（ 12）：19-22.

[2] 陸宏圻.楊勇，何培杰，等，小浪底水庫射流沖吸式清淤設備研究[J].人民黃河，2011，33（4）：15-16.

[3]RAJARATNAM N， MAZUREK K A.Erosion of Sand by Cir-cular Impinging Water Jets with Small Tailwater[J].Joumalof Hydraulic Engineering， 2003， 129（3）：225-229.

[4]齊梅蘭，府仁壽，陳稚聰，射流沖刷平衡深度研究[J].水動力學研究與進展，2005，20（3）：368-372.

[5] 槐文信，王增武，錢忠東，等，二維垂向射流沙質河床沖刷的數值模擬[J].中國科學：技術科學，2012（1）：72-81.

[6] 張浩，倪福生，顧磊，某靶距下射流垂直沖刷坑深的試驗研究[J].科學技術與工程，2014，14（2）：238-240.

[7] 申振，徐立群，張浩，等，射流作用下不同沙床動態沖坑變化規律的試驗研究[J].水資源與水工程學報，2015（2）：182-185.

[8] 顧磊，倪福生，徐立群，等，垂直射流沖刷沙床的噴嘴間距試驗研究[J].泥沙研究，2017，42（3）：54-58.

[9]劉思源，泥沙粒徑及射流速度對沙床沖刷的影響[J].人民長江，2017，48（ 23）：79-83.

[10] 范家驊，異重流與泥沙工程試驗與設計[M].北京：中國水利水電出版社，2011：320-355.

[11]鄒志利，海岸動力學[M].北京：人民交通出版社，2009：176-180.

[12]董志勇，射流力學[M].北京：科學出版社，2005：15-23.

[13]童思陳，許光祥，河灣岸坡泥沙起動流速研究[J].水利學報，2008，39（ 11）：1167-1173.

[14]ADERIBICBE 0 0， RAJARATNAM N.Erosion of Loose Bedsby Submerged Circular Impinging Vertical Turbulent Jets[J].Journal of Hydraulic Research， 1996， 34（1）：19-33.