摘要：卵石起動(dòng)流速是泥沙起動(dòng)流速的一種特殊形式。國內(nèi)許多專家學(xué)者對(duì)于既適用于散粒體又適用于黏性細(xì)顆粒的泥沙起動(dòng)流速公式進(jìn)行過大量研究，究竟哪一種公式更適合卵石起動(dòng)流速的計(jì)算與應(yīng)用尚無定論。為了探尋這些公式的合理性與實(shí)用性，依托水工模型試驗(yàn)水深測試值和理論值，對(duì)于一些常用的公式進(jìn)行了分析計(jì)算，又利用防沖槽沖坑深度反算卵石的起動(dòng)流速，并且將試驗(yàn)計(jì)算值與理論計(jì)算值以及沖深反算值進(jìn)行對(duì)比，間接地得出了長科院公式在卵石起動(dòng)流速計(jì)算與應(yīng)用方面更具有實(shí)用性的結(jié)論，以及該公式適用卵石粒徑的范圍。

關(guān)鍵詞：卵石; 起動(dòng)流速; 粒徑范圍; 水工模型試驗(yàn)

中圖法分類號(hào)： TV142

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

DOI：10.16232/j.cnki.1001-4179.2021.11.033

0引 言

現(xiàn)代地貌形態(tài)取決于內(nèi)、外營力的相互作用。內(nèi)營力指構(gòu)造作用所造成的地殼升降，從而構(gòu)成侵蝕、堆積、環(huán)境變化的宏觀地質(zhì)背景;外營力指太陽輻射、降水及微生物活動(dòng)等導(dǎo)致巖石風(fēng)化、剝蝕及輸移。由于侵蝕營力造成了水土流失，進(jìn)而產(chǎn)生了泥沙。卵石是較大顆粒的泥沙，是一種特殊的表現(xiàn)形式。本文對(duì)于自然界分類的卵石、漂石、頑石一律采用卵石的通俗說法，即均稱之為卵石（此處卵石粒徑≥15 mm）。

眾所周知，泥沙研究是水利專業(yè)學(xué)科研究的重要分支，其中的泥沙起動(dòng)流速研究又是泥沙研究的主要內(nèi)容。研究泥沙起動(dòng)問題時(shí)，往往從泥沙起動(dòng)流速入手[1]。早在200 多年前，就有人開始對(duì)泥沙起動(dòng)流速進(jìn)行研究。起動(dòng)流速是研究河流泥沙運(yùn)動(dòng)規(guī)律的重要參量，也是泥沙研究的難題之一[2]。推移質(zhì)運(yùn)動(dòng)是泥沙運(yùn)動(dòng)的一種重要形式[3]。而對(duì)于新疆眾多的山溪性多泥沙河流上修建的引水渠首，在較大的卵石防沙、排沙以及利用其進(jìn)行消能防沖方面，常常要對(duì)卵石這種推移質(zhì)進(jìn)行起動(dòng)流速的研究和計(jì)算。卵石的起動(dòng)有其特殊性，加強(qiáng)其起動(dòng)流速及起動(dòng)規(guī)律的研究，對(duì)工程泥沙及泥沙運(yùn)動(dòng)規(guī)律具有重要的價(jià)值。關(guān)于卵石起動(dòng)較早的一些研究，都是按一般泥沙運(yùn)動(dòng)概念進(jìn)行的。

一些研究以單顆粒泥沙的受力分析為出發(fā)點(diǎn)，考慮顆粒間黏結(jié)力，建立臨界起動(dòng)切應(yīng)力或者起動(dòng)流速公式[4]。較為典型的泥沙起動(dòng)研究分為兩個(gè)大的方面：“前蘇聯(lián)和中國”模式（以流速為主要參數(shù)）和“歐美”模式（以剪切力為主要參數(shù)）。20世紀(jì)30年代，Shields等分析了床面上泥沙顆粒的受力平衡，推導(dǎo)出無黏性均勻沙的起動(dòng)切應(yīng)力公式。Rouse根據(jù)Shields及他人的試驗(yàn)資料，繪測出著名的Shields曲線并推導(dǎo)出泥沙起動(dòng)的拖曳力公式。此公式曾一度被歐美工程師用作判別泥沙初始運(yùn)動(dòng)的唯一準(zhǔn)則[5]。還有學(xué)者通過研究得到了無黏性均勻沙的臨界起動(dòng)切應(yīng)力公式[6]，也有研究認(rèn)為只有沙粒切應(yīng)力τ才對(duì)推移質(zhì)泥沙運(yùn)動(dòng)起作用[7]。隨后，又有學(xué)者提出公式τ1=γhJ′，根據(jù)Shields曲線，在顆粒雷諾數(shù)Re>500時(shí)，τc=0.06，而通過大量試驗(yàn)，得到無因次起動(dòng)切應(yīng)力τc=0.047;對(duì)于非均勻沙，有研究提出粗糙紊流中泥沙50%起動(dòng)概率時(shí)相應(yīng)的無因次起動(dòng)切應(yīng)力τc50=0.046，而根據(jù)實(shí)測資料都得到紊流中τc50=0.045;通過系統(tǒng)分析80 a來泥沙起動(dòng)的研究成果和資料后得到，在沙卵石河流中，高雷諾數(shù)時(shí)，τc50在0.052～0.086之間，在低雷諾數(shù)時(shí)，τc50在0.030～0.073之間[8]。范玉通過研究表明對(duì)于散粒體粗顆粒泥沙，泥沙的起動(dòng)拖曳力只與泥沙粒徑有關(guān)，而和沙粒雷諾數(shù)無關(guān)[9]。

不過，基于流場和剪切力場之間的關(guān)系，更多的學(xué)者則是通過研究流速來判定泥沙的起動(dòng)。中國在泥沙起動(dòng)流速的理論研究上已達(dá)到了國際領(lǐng)先水平。目前僅針對(duì)非黏性泥沙提出的起動(dòng)流速公式已超過100種[10]，其中竇國仁公式、張瑞瑾公式、唐存本公式、沙玉清公式等都有很大的影響。

張瑞瑾利用實(shí)測資料對(duì)系數(shù)與指數(shù)率定得到既適用于散粒體又適用于黏性細(xì)顆粒泥沙的統(tǒng)一起動(dòng)流速公式，公式如下：

前3個(gè)公式雖然沒有交代適用粒徑范圍，但都明確了可以應(yīng)用于粗顆粒泥沙起動(dòng)流速計(jì)算，公式（4）則是明確了適用于卵石起動(dòng)流速計(jì)算。由此，本文依托水工模型試驗(yàn)水深測試值和理論值，對(duì)于這4個(gè)公式進(jìn)行了分析計(jì)算，又利用防沖槽沖坑深度反算卵石的起動(dòng)流速，并且將試驗(yàn)計(jì)算值進(jìn)行對(duì)比，得出更適用于卵石起動(dòng)流速計(jì)算的公式。

1模型設(shè)計(jì)和試驗(yàn)

1.1試驗(yàn)?zāi)康暮蛢?nèi)容

在具有一定泥沙組成的床面上，逐漸增加水流強(qiáng)度，直到使床面泥沙（簡稱床沙）由靜止轉(zhuǎn)入運(yùn)動(dòng)，這種現(xiàn)象稱為泥沙的起動(dòng)[12]。相應(yīng)的臨界水流條件稱為泥沙的起動(dòng)條件。泥沙的起動(dòng)條件是泥沙的基本水力特性之一，同時(shí)，它又是河床沖刷的臨界條件。確定泥沙起動(dòng)的水力條件是泥沙運(yùn)動(dòng)力學(xué)的重要內(nèi)容之一。吳彰松等[13]采用泥沙滾動(dòng)起動(dòng)模型開展了這方面的研究，劉蛟等[14]進(jìn)行了三維水流作用下散列體泥沙的起動(dòng)研究，李林林[15]也開展了卵礫石的起動(dòng)研究，張根廣等[16-17]研究了粗顆粒泥沙起動(dòng)機(jī)理和斜坡上泥沙起動(dòng)條件，毛寧[18]通過研究得到了塊石起動(dòng)流速公式，但以上研究與新疆多數(shù)修建于出山口較平緩開闊地帶引水渠首的卵礫石起動(dòng)規(guī)律有所不同。本文通過建立一個(gè)新疆山溪性多沙河流引水渠首閘后設(shè)置拋石防沖槽的標(biāo)準(zhǔn)斷面模型，在施放3種不同流量、防沖槽中拋填6種不同粒徑的卵石試驗(yàn)條件下，得出防沖槽正常的槽中水深試驗(yàn)值和理論值，為槽中卵石起動(dòng)流速的計(jì)算提供基礎(chǔ)參數(shù)。該模型試驗(yàn)主要用來研究防沖槽內(nèi)水沙特性[19-20]，本次對(duì)于起動(dòng)流速的分析只是采用了模型中的一部分試驗(yàn)數(shù)據(jù)，模型試驗(yàn)方案見表1。

1.2模型設(shè)計(jì)和制作

模型比尺為1∶30。

模型范圍：以引水渠首閘后防沖槽為基準(zhǔn)進(jìn)行確定，設(shè)計(jì)為標(biāo)準(zhǔn)斷面模型[21]，原型寬度27 m（見圖1）。

1.3試驗(yàn)典型流量選取

全疆共有大小河流570條，其中年徑流量0.74億～2.39億m3的河流有百余條，這些河流上修建的引水渠首泄洪單寬流量為8.21～17.52 m2/s，而且這些渠首有很多都采用了拋石防沖槽的防沖結(jié)構(gòu)[22]。因此，結(jié)合試驗(yàn)場地條件，本次研究選取3個(gè)典型單寬流量（見表2）[23]，有其代表性與合理性。模型及試驗(yàn)如圖2～5所示。

2流速計(jì)算與對(duì)比分析

通過前述不同公式的表述可以看出，起動(dòng)流速主要與正常水深、卵石粒徑、容重、重力加速度、薄膜水厚度和孔隙率（這兩個(gè)參數(shù)主要為沙玉清公式應(yīng)用）等有關(guān)。本次沒有通過試驗(yàn)直接去測試卵石的起動(dòng)流速，而是通過與起動(dòng)流速有關(guān)的水深和沖刷坑深來計(jì)算起動(dòng)流速。起動(dòng)流速應(yīng)大于或等于某一個(gè)界限值，不沖流速應(yīng)小于這個(gè)界限值，這兩個(gè)界限值應(yīng)是同一個(gè)數(shù)值，小于不沖流速，卵石就不會(huì)起動(dòng);大于或者等于起動(dòng)流速，卵石就會(huì)起動(dòng)從而造成沖刷。將試驗(yàn)測試的水深值代入前述不同公式得出起動(dòng)流速試驗(yàn)計(jì)算值，又將防沖槽近似看作邊墻現(xiàn)澆混凝土、底部堆砌均一卵石的渠槽，據(jù)此計(jì)算得到槽中正常水深，也代入前述公式得到一個(gè)起動(dòng)流速理論計(jì)算值，再通過查閱相關(guān)不沖流速資料，在同等范圍條件下，將起動(dòng)流速與查詢的不沖流速進(jìn)行對(duì)比分析，對(duì)于起動(dòng)流速與不沖流速接近的數(shù)值確定為可信值，從而得出最適用于卵石起動(dòng)的起動(dòng)流速公式。

2.1起動(dòng)流速計(jì)算

2.1.1采用理論水深值的起動(dòng)流速計(jì)算

防沖槽近似看作渠槽的理論水深值計(jì)算見表3。

2.1.2采用試驗(yàn)水深值的起動(dòng)流速計(jì)算

通過拋填0.2～0.8 m不同粒徑卵石，分別施放3種流量，采用長科院、張瑞瑾、沙玉清、崗恰洛夫公式進(jìn)行起動(dòng)流速計(jì)算。發(fā)現(xiàn)以下規(guī)律：① 長科院公式和張瑞瑾公式的試驗(yàn)值與理論值起動(dòng)流速比較接近，反映出了理論與實(shí)際的一致性，沙玉清公式和崗恰洛夫公式在理論與實(shí)際的一致性方面偏差較大;② 各拋石粒徑下，卵石起動(dòng)流速與施放流量關(guān)系不大，隨著拋石粒徑的增大，卵石起動(dòng)流速也隨之增大;③ 拋填0.2～0.6 m卵石，各種工況條件下，隨著拋石粒徑的增大，卵石起動(dòng)流速也隨之增大，但是拋填0.8 m卵石條件下，起動(dòng)流速呈現(xiàn)明顯降低現(xiàn)象。計(jì)算統(tǒng)計(jì)結(jié)果見表4。

2.1.3采用試驗(yàn)沖刷坑深反算起動(dòng)流速

防沖槽沖坑最深處可認(rèn)為是沖刷最劇烈的位置，試驗(yàn)中分別測出了沖刷坑處的沖刷水深h（含沖坑深）。根據(jù)沖刷極限狀態(tài)法，有Q=BhU;另外，河床抗沖粗化層的計(jì)算可以通過水流連續(xù)公式q=Uh來計(jì)算。此兩式中Q為流量，B為寬度，h為沖刷達(dá)到極限狀態(tài)的水深即為沖刷水深，U為起動(dòng)流速，q=Q/B。因此，單寬流量q與沖刷水深h的比值即為起動(dòng)流速U。采用防沖槽沖坑試驗(yàn)值反算的起動(dòng)流速統(tǒng)計(jì)見表5。

從表5可以看出，在同一流量下，隨著卵石粒徑的增大，沖深反算值起動(dòng)流速總體呈現(xiàn)出上升趨勢，反映出了試驗(yàn)的客觀規(guī)律;由于試驗(yàn)過程中沖坑深與水深測量總會(huì)存在一些人為因素造成的誤差，使得該沖深反算值產(chǎn)生了異常的波動(dòng)趨勢，但就數(shù)值而言，還是有一定可靠度的。

2.2對(duì)比分析

將起動(dòng)流速的上述計(jì)算結(jié)果進(jìn)行對(duì)比分析，發(fā)現(xiàn)沖深反算值最小;其次是長科院公式的試驗(yàn)計(jì)算值和理論計(jì)算值（這兩個(gè)值比較接近），吻合程度最高;根據(jù)前述沙玉清公式和崗恰洛夫公式得到的試驗(yàn)計(jì)算值與卵石起動(dòng)流速吻合程度最差，張瑞瑾公式計(jì)算值相比長科院公式值較高，高出長科院公式值0.67～1.13 m。這也反映出現(xiàn)有公式在滲流作用明顯處，計(jì)算泥沙起動(dòng)流速時(shí)，會(huì)大于實(shí)際起動(dòng)流速[24]。將4個(gè)公式的試驗(yàn)計(jì)算值、長科院公式理論計(jì)算值與沖深反算值放在一起進(jìn)行對(duì)比（見表6）。

為了更清晰地對(duì)比，按照不同流量將表6中數(shù)值分別繪圖分析，如圖6所示。

防沖槽可近似看作渠槽，渠底卵石一般是松散堆積的，按照有關(guān)規(guī)范，可當(dāng)作非平面形卵石、砌筑來考慮，對(duì)于水力半徑等于2，粒徑為200，250，300 mm時(shí)，不沖流速分別為3.0，3.5，3.7 m/s[25];還有規(guī)范規(guī)定，干砌卵石渠道允許不沖流速為2.5～4 m/s[26]。從試驗(yàn)數(shù)據(jù)來看，當(dāng)流量為220 m3/s，卵石粒徑由200 mm逐漸增大至800 mm時(shí)，沖刷斷面處的最大流速由3.21 m/s逐漸增大至5.58 m/s;當(dāng)流量為300 m3/s，卵石粒徑由200 mm逐漸增大至800 mm時(shí)，沖刷斷面處的最大流速由4.16 m/s逐漸增大至8.15 m/s;當(dāng)流量為400 m3/s，卵石粒徑由200 mm逐漸增大至800 mm時(shí)，沖刷斷面處的最大流速由4.44 m/s逐漸增大至9.52 m/s。這些流速值是已經(jīng)發(fā)生沖刷時(shí)的數(shù)據(jù)，而起動(dòng)流速肯定會(huì)比沖刷時(shí)的流速值小很多才可能造成之后較大的沖刷。沙玉清公式明顯數(shù)值偏大，崗恰洛夫公式在本次試驗(yàn)中最大流量、最大卵石粒徑條件下數(shù)值反而下降很多，且多數(shù)數(shù)值基本大于沖深反算值的1.5～2倍，與現(xiàn)行規(guī)范也不相符，不適用于此類卵石起動(dòng)流速計(jì)算;張瑞瑾公式數(shù)值基本大于沖深反算值的1.0～2.6倍，而且200，300 mm粒徑的流速值也不滿足現(xiàn)行規(guī)范的要求，也不是理想的公式;只有長科院公式試驗(yàn)計(jì)算值與理論計(jì)算值較為接近，不但滿足現(xiàn)行規(guī)范的要求，而且除了800 mm粒徑外，沖坑深反算值均小于長科院公式值，并且也發(fā)現(xiàn)卵石粒徑為800 mm時(shí)，試驗(yàn)值規(guī)律性、代表性不強(qiáng)。

該試驗(yàn)屬于引水樞紐下游泥沙沖刷的水力學(xué)研究范疇，雖然與河流動(dòng)力學(xué)研究范疇有所差異，但是下游泥沙沖刷與拋石的起動(dòng)流速有著直接的關(guān)系，也就因此與河流動(dòng)力學(xué)范疇有了緊密聯(lián)系，而且防沖槽內(nèi)拋石的起動(dòng)與沖刷比之天然河道卵石起動(dòng)情況要復(fù)雜的多，水流條件也劇烈很多，開展此項(xiàng)研究具有一定意義。通過上述分析可以初步認(rèn)為：長科院公式適用于卵石粒徑d≤600 mm的起動(dòng)流速計(jì)算，d>600 mm粒徑卵石的起動(dòng)流速計(jì)算還需要去除顆粒之間黏結(jié)力等因素開展進(jìn)一步研究。

3結(jié) 論

（1） 經(jīng)過驗(yàn)證，文中提到的可用于散粒體或粗顆粒泥沙、卵石起動(dòng)流速的4個(gè)經(jīng)驗(yàn)公式中，長科院公式理論值與實(shí)際值接近，用于卵石起動(dòng)流速計(jì)算最為可靠和實(shí)用。

（2） 卵石起動(dòng)流速應(yīng)大于或等于界限值，卵石不沖流速應(yīng)小于界限值，這兩個(gè)界限值應(yīng)是同一個(gè)數(shù)值;小于不沖流速，卵石就不會(huì)起動(dòng);大于或者等于起動(dòng)流速，卵石就會(huì)起動(dòng)從而造成沖刷。

（3）長科院公式適用于卵石粒徑d≤600 mm的起動(dòng)流速計(jì)算。

參考文獻(xiàn)：

[1]蔡蓉蓉，張羅號(hào)，張紅武.李保如泥沙起動(dòng)流速公式修正研究[J].水利學(xué)報(bào)，2019，50（5）：547-554.

[2]馬志敏，鄒先堅(jiān)，趙小紅，等.基于B超成像的泥沙起動(dòng)流速測量方法[J].水科學(xué)進(jìn)展，2015，26（1）：59-59.

[3]呂升奇，郭美嬌，袁禮偉.淹沒柔性植物對(duì)明渠推移質(zhì)運(yùn)動(dòng)影響的實(shí)驗(yàn)研究[J].人民長江，2016，47（2）：67-71.

[4]劉潔，劉潔玉，白玉川.黏性泥沙流變特性及其臨界起動(dòng)的研究[J].泥沙研究，2015（6）：59-64.

[5]閆紅杰，史雅君，周萍，等.河床泥沙起動(dòng)流速公式修正及可靠性驗(yàn)證[J].科學(xué)技術(shù)與工程，2018，18（17）：297-301.

[6]李林林，張根廣.滲流作用下岸坡泥沙起動(dòng)流速研究[J].人民長江，2017，48（23）：84-88.

[7]黃才安，奚斌.推移質(zhì)輸沙率公式的統(tǒng)一形式[J].水利水運(yùn)科學(xué)研究，2000（2）：72-78.

[8]湯潔.壅水條件下水流結(jié)構(gòu)及床面泥沙起動(dòng)研究[D].武漢：長江科學(xué)院，2016.

[9]范玉.水沙兩相流理論及引水工程管道輸沙問題的研究[D].天津：天津大學(xué)，2014.

[10]朱正雷，邵宇陽，趙曉晴.基于FLUENT的便攜式泥沙起動(dòng)測量裝置模擬效率與精度[J].水電能源科學(xué)，2019，37（6）：156-156.

[11]水利部水利水電規(guī)劃設(shè)計(jì)總院.堤防工程設(shè)計(jì)規(guī)范：GB 50286-2013[S].北京：中國計(jì)劃出版社，2013.

[12]周雙，謝飛，李福育，等.負(fù)坡上均勻散粒體泥沙起動(dòng)流速研究[J].中國水運(yùn)，2016，16（12）：225-225.

[13]吳彰松，張根廣，寇潭，等.基于概率統(tǒng)計(jì)理論的泥沙起動(dòng)流速研究[J].長江科學(xué)院院報(bào)，2017，34（7）：7-11.

[14]劉蛟，嚴(yán)軍，劉鐵，等.三維水流作用下散列體泥沙的起動(dòng)研究[J].泥沙研究，2015（1）：7-13.

[15]李林林.基于統(tǒng)計(jì)理論的推移質(zhì)起動(dòng)概率及輸移規(guī)律[D].楊凌：西北農(nóng)林科技大學(xué)，2018.

[16]張根廣，周雙，邢茹，等.基于相對(duì)暴露度的無黏性均勻泥沙起動(dòng)流速公式[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2016，24（4）：687-697.

[17]張根廣，王新雷，吳彰松，等.泄水建筑物下游沖刷坑二維數(shù)值模擬[J].泥沙研究，2017，42（3）：42-47.

[18]毛寧.論泥沙礫石的起動(dòng)流速[J].長江科學(xué)院院報(bào)，2011，28（1）：6-11.

[19]李江峰，尹輝，耿凡坤，等.拋填不同粒徑卵石對(duì)閘后防沖槽的影響研究水工模型試驗(yàn)研究報(bào)告[R].烏魯木齊：新疆水利水電科學(xué)研究院，2018.

[20]李江峰，李娟，尹輝，等.基于標(biāo)準(zhǔn)斷面的防沖槽水沙特性研究水工模型試驗(yàn)研究報(bào)告[R].烏魯木齊：新疆水利水電科學(xué)研究院，2019.

[21]李江峰，尹輝，戚印鑫，等.拋石粒徑對(duì)防沖槽結(jié)構(gòu)和沖刷的影響研究[J].水利與建筑工程學(xué)報(bào)，2019，17（3）：115-120.

[22]戚印鑫，張磊，孫娟，等.白楊河阿克蘇渠首除險(xiǎn)加固工程可行性研究報(bào)告[R].烏魯木齊：新疆水利水電科學(xué)研究院，2018.

[23]李江峰，尹輝，戚印鑫，等.渠首閘后防沖槽的拋石粒徑與消能試驗(yàn)研究[J].水利與建筑工程學(xué)報(bào)，2019，17（6）：117-181.

[24]閆紅杰，史雅君，周萍，等.河床泥沙起動(dòng)流速公式修正及可靠性驗(yàn)證[J].科學(xué)技術(shù)與工程，2018，18（17）：297-301.

[25]中華人民共和國水利部.渠道防滲工程技術(shù)規(guī)范：GB/T 50600-2010[S].北京：中國計(jì)劃出版社，2011.

[26]中華人民共和國水利部.灌溉與排水工程設(shè)計(jì)標(biāo)準(zhǔn)：GB 50288-2018[S].北京：中國計(jì)劃出版社，2018.

（編輯：謝玲嫻）

Abstract：Pebble starting velocity is a special form of sediment starting velocity.Many experts and scholars have conducted a large number of studies on the starting velocity formula which is applicable to both loose particles and viscous fine particles，but there is still no final conclusion on which formula is more suitable for the pebble starting velocity.In order to explore the rationality and practicability of these formulas，we relied on the water depth test values of hydraulic model test and theoretical value to analyze and calculate some commonly used formulas，and used the anti-flushing groove depth to inversely deduce the starting velocity of the pebble.Then we compared the experimental calculation values，theoretical calculation values and the inversed deduced values of the scouring depth.The results showed that the formula of Changjiang River Scientific Research Institute was more practical in the calculation and application of the pebble starting velocity，and the range of the particle size of the pebble applicable to the formula was obtained.

Key words：pebble;starting velocity;range of the particle size;hydraulic model test