摘要：為探究具有堰塞壩特征的大比降天然非均勻沙河道工程的泥沙顆粒起動流速，在前人基礎(chǔ)上重新定義隱蔽度的概念，綜合考慮泥沙顆粒的隱暴效應(yīng)，采用滾動起動模式對斜坡上泥沙顆粒進(jìn)行受力分析，建立斜坡非均勻無黏性泥沙顆粒的垂線平均流速公式.采用實(shí)測資料對該公式及前人總結(jié)出的公式進(jìn)行檢驗(yàn)，結(jié)果表明本文計算值與實(shí)測值吻合度更高，誤差更小.

關(guān)鍵詞：起動流速；滾動起動；無黏性泥沙；隱蔽度；暴露度

中圖分類號：TV142"" 文獻(xiàn)標(biāo)志碼：A" 文章編號：2095-6991（2025）01-0078-06

Study on Incipient Velocity Formulaof Non-uniform and Non-viscous Sediment Particles on Slope

LI Nan， TAN Long-jin， LIAO Hai-mei

（College of Civil Engineering， Guizhou University， Guiyang 550025， China）

Abstract：In order to explore the problem of sediment particle starting flow rate of natural non-uniform sand channel project with large specific drop characteristic of landslide dam， the rolling start mode was adopted to redefine the concept of hidden degree， and comprehensively considered the hidden degree and exposure degree in the upstream and downstream of sediment particles. Based on force analysis of the sediment particles on slope， the computing vertical mean velocity formula of non-uniform and non-viscous sediment particles on the slope were established. The results show that the calculated values in this paper are in better agreement with the measured values and the errors are smaller.

Key words：incipient velocity; rolling start mode; non-viscous sediment; hidden degree; exposure degree

0 引言

堰塞壩是由山崩、滑坡及泥石流等失穩(wěn)地質(zhì)體阻塞山谷、河道形成的堆積體^[1].壩體含有大量碎屑基質(zhì)和塊石，土石料空間分布不均勻性強(qiáng)，易發(fā)生滲流破壞.其中，中粗粒徑土石料顆粒的起動導(dǎo)致壩體加速失穩(wěn)潰決，將會給上下游及周邊地區(qū)帶來嚴(yán)重的生命及財產(chǎn)損失.潰決過程是漫壩水流對壩體土石料的沖刷侵蝕過程，與河道中的泥沙輸移有相似之處^[2-3].壩體土石料在堰塞湖下游河床發(fā)生堆積，促使下游段河道抬高.由于起動流速是研究河床演變的重要力學(xué)指標(biāo).因此，本文在泥沙研究理論基礎(chǔ)上對堰塞壩潰決過程中壩體內(nèi)中粗粒徑土石料顆粒起動條件進(jìn)行研究.

泥沙起動的水力指標(biāo)可用起動流速、起動切應(yīng)力和起動功率表示.其中，泥沙起動流速在水利工程的實(shí)際應(yīng)用中具有重要意義.早期學(xué)者將泥沙起動作為確定性問題進(jìn)行研究，后將概率論與泥沙力學(xué)相結(jié)合，考慮了起動過程中的隨機(jī)性^[4].竇國仁^[5]明確了河底泥沙的受力情況，并得出適用于沙石、粘土和黃土的泥沙起動流速公式.謝立全等^[6]、吳巖等^[7]考慮了水在孔隙中流動產(chǎn)生的滲透力，推導(dǎo)出岸坡上受到滲透力作用的泥沙起動流速公式.蔡蓉蓉等^[8]對文獻(xiàn)[9]中起動流速公式的摩擦系數(shù)函數(shù)進(jìn)行修正，克服原公式繁瑣的缺點(diǎn).

泥沙起動涉及復(fù)雜的水沙耦合作用，水流速度與泥沙顆粒相對位置都具有隨機(jī)性，泥沙顆粒間相對位置在河床中極其復(fù)雜.韓其為^[10]引入了暴露度以及相對暴露度概念來描述泥沙顆粒之間相對位置對泥沙起動的影響.何文社等^[11]提出附加質(zhì)量力概念，認(rèn)為其與泥沙顆粒在河床上的相對位置有直接關(guān)系.楊奉廣等^[12]通過經(jīng)典力學(xué)與概率論結(jié)合的辦法，對均勻沙的暴露度概率密度函數(shù)進(jìn)行修正，得出均勻沙大比降起動流速公式并應(yīng)用到唐家山堰塞湖下游河床，發(fā)現(xiàn)加大水深會使泥沙起動流速增大.WU^[13]、DUAN等^[14]在泥沙顆粒受力分析中考慮了上游顆粒的阻礙作用.馬菲等^[15]提出遮掩度概念，細(xì)化了泥沙顆粒在河床上的位置分布情況.白玉川等^[16]分析了三維維度下大粒徑泥沙顆粒的暴露度，考慮了局部繞流的影響.王愉樂等^[17]提出隱蔽角概念，體現(xiàn)出泥沙上游顆粒的隱蔽作用.麻卉港等^[18]通過試驗(yàn)研究結(jié)合暴露角分布概率，建立了考慮泥沙相對位置和水流條件的均勻沙起動全概率公式.

目前大多數(shù)學(xué)者在關(guān)于泥沙隱暴度研究上側(cè)重于下游泥沙顆粒的阻礙作用，而忽略了上游顆粒的隱蔽作用.相對暴露度的定義也側(cè)重于均勻沙顆粒之間的隱蔽關(guān)系.暴露度作為力矩平衡中的重要影響系數(shù)^[19]，直接關(guān)系泥沙顆粒受力分析的合理性，所以充分考慮泥沙顆粒在河床中的相對位置關(guān)系顯得極其重要.即使部分學(xué)者考慮到泥沙上游顆粒的隱蔽作用，但較多應(yīng)用于均勻沙或小比降非均勻沙.對于具有堰塞壩特征的大比降天然非均勻沙河道工程泥沙顆粒起動流速的研究甚少.初始水動力條件直接影響堰塞壩潰決時間長短和潰決洪水流量大小，明確壩體土石料顆粒的起動流速可為應(yīng)急搶險工作提供參考.為此，在泥沙研究的理論基礎(chǔ)上，充分考慮了上游顆粒的隱蔽作用，推導(dǎo)出滾動起動模式下斜坡非均勻土石料的起動流速公式，為改進(jìn)土石料顆粒起動流速及堰塞壩潰決過程關(guān)鍵水力參數(shù)的計算方法提供理論參考.

1 隱蔽度及暴露度分析

為體現(xiàn)泥沙顆粒在斜坡上位置排列的一般性，繪制了如圖1所示的泥沙顆粒位置分析圖.

圖中泥沙顆粒位于縱坡坡腳為β的斜坡上.根據(jù)韓其為相關(guān)定義可知，暴露度是泥沙顆粒最低點(diǎn)與相鄰下游顆粒接觸點(diǎn)在床面上的豎直距離，相對暴露度是暴露度與泥沙顆粒半徑的比值.周雙等^[20]定義暴露角θ為泥沙顆粒的圓心與下游顆粒接觸點(diǎn)之間的連線與床面垂線夾角.暴露度和相對暴露度并不能反映上游顆粒的隱蔽作用.實(shí)際上，即使兩個泥沙顆粒有相同暴露度和相對暴露度，若受到上游顆粒的隱蔽作用不同，兩者起動概率也可能不相同.

例如，如圖1所示的顆粒D1與顆粒D4為同一直徑，顆粒D2與顆粒D3為同一直徑，顆粒D5的直徑介于二者之間.在同一水流條件下，盡管顆粒D1與顆粒D4的暴露度和相對暴露度完全相同，由于顆粒D5對D4的隱蔽作用大于顆粒D2對D1的隱蔽作用，使得顆粒D1與D4的起動概率并不相同.因此分析泥沙顆粒之間相對位置應(yīng)綜合考慮下游顆粒的阻礙作用和上游顆粒的隱蔽作用.王愉樂等^[17]提出隱蔽角α，定義其為泥沙顆粒圓心與上游顆粒接觸點(diǎn)之間的連線與床面垂線夾角，并假設(shè)隱蔽角α與暴露角θ服從相同的分布規(guī)律.馬菲等^[15]提出了遮掩度，定義其為泥沙顆粒與上游顆粒接觸點(diǎn)和上游顆粒最低點(diǎn)之間的距離，[JP2]但遮掩度與泥沙顆粒粒徑相關(guān)性較弱.本文為使隱蔽度和所研究顆粒關(guān)聯(lián)性更強(qiáng)，重新定義隱蔽度及相對隱蔽度概念.隱蔽度1定義為泥沙顆粒最低點(diǎn)與上游泥沙顆粒接觸點(diǎn)在床面上的豎直距離.相對隱蔽度1定義為隱蔽度與泥沙顆粒半徑的比值.韓其為等^[10]提出的相對暴露度和本文提出的相對隱蔽度1分別如式（1）及式（2）所示：

Δ′=Δ/R，（1）

1′=1/R，（2）

式中，R為所研究泥沙顆粒半徑.

這里假設(shè)相對隱蔽度和相對暴露度服從相同的分布規(guī)律，結(jié)合韓其為對泥沙顆粒相對暴露度統(tǒng)計值可以得出在非均勻沙中，粗、中、細(xì)3種不同粒徑的泥沙顆粒相對隱蔽度的數(shù)學(xué)期望值：

當(dāng)D≥D/Δ′m=7.46D時為粗顆粒，其相對隱蔽度的數(shù)學(xué)期望值如式（3）：

1′=D/D；（3）

當(dāng)D/Δ′m＞D＞D時為中顆粒，其相對隱蔽度的數(shù)學(xué)期望值如式（4）：

1′=0.5D/D+Δ′m；（4）

當(dāng)D≤D時為細(xì)顆粒，其相對隱蔽度的數(shù)學(xué)期望值如式（5）：

1′=0.5D/D－Δ′m²/D/D－Δ′m；（5）

式中，D為泥沙顆粒直徑，D為泥沙顆粒平均直徑，Δ′m=0.134.

根據(jù)圖1中的幾何關(guān)系可以得到：

cosθ=R－Δ/R=1－Δ′ ，（6）

sinθ=1－cosθ²=2Δ′－Δ′²，（7）

cosα=R－1/R=1－1′，（8）

sinα=1－cosα²=21′－1′².（9）

式中，θ [0，90-β]；α [0，90-β].

2 泥沙受力分析

部分學(xué)者在研究泥沙顆粒隱暴效應(yīng)時多以水平面作為基準(zhǔn)面進(jìn)行探討，水平面的坡腳始終為零，可以理解為所設(shè)立的基準(zhǔn)面是固定不變的，而實(shí)際中斜坡坡度隨著泥沙輸移過程不斷變化，是一個動態(tài)值.堰塞壩的潰決過程中壩體土石料輸移強(qiáng)度較大，其基準(zhǔn)面隨潰決過程不斷發(fā)生變化.鑒于此，本文參考DEY^[21]對泥沙顆粒所處坡面的設(shè)定，結(jié)合堰塞壩實(shí)際潰壩過程中土石料的運(yùn)動規(guī)律，對土石料顆粒進(jìn)行如圖2所示的受力分析.

土石料顆粒位于橫坡坡腳為α1，縱坡坡腳為β的斜面上，斜面外法線與鉛垂方向的夾角為γ，各個角度值隨潰壩過程發(fā)生變化.γ與α1、β的關(guān)系如式（10）所示：

cosγ=cosβ²－sinα1².（10）

黏性細(xì)顆粒泥沙多以成片或成團(tuán)形式進(jìn)行起動，而本文建立的受力模型以單顆粒狀態(tài)起動，不適用于黏性細(xì)顆粒的起動形式，這里重點(diǎn)考察非黏性顆粒起動特性.黏性土石料所受力有：水下重力G、滲透力FS、拖曳力FD、上舉力FL、附加下壓力FM及黏著力Fμ.其中水下重力G、附加下壓力FM及黏著力Fμ方向鉛垂向下，滲透力FS、上舉力FL方向?yàn)樾泵嫱夥ň€方向，拖曳力FD方向與水流方向一致，各個物理量的值如式（11）～式（16）所示.對于土石料中的非黏性顆粒，認(rèn)為附加下壓力FM與黏著力FSymbolmA@的值為零.

G=π6ρs－ρgD3，（11）

FS=iρgπD361－e，（12）

FD=18CDπρD²u²，（13）

FL=18CLπρD²u²，（14）

FM=3πK2ρgHR3－tδ1δ1－t，（15）

Fμ=32πq0Rδ303－tδ11t²－tδ²1， （16）

式中：ρs和ρ分別為泥沙及水的密度；g為重力加速度，取9.81 m/s²；D為所研究泥沙顆粒的直徑；i為滲透比降；e為孔隙率；CD為水流拖曳力系數(shù)，取0.4；CL為水流上舉力系數(shù)，取0.1；u為泥沙起動時底部瞬時流速，K2=2.25810^-3；t為土石料之間的縫隙；δ1為薄膜水的厚度，取值為δ1=410^-7；δ0為一個水分子的厚度，取值為δ0=310^-10；q0=1.3109 kg/m².

3 泥沙顆粒起動流速公式建立

鑒于堰塞壩潰決過程中大部分顆粒為滾動起動^[22]，少部分顆粒為滑動起動和躍起起動，故本文采用滾動起動模式，根據(jù)圖2泥沙顆粒受力圖，由力矩平衡方程可得：

FDLD+FLLL+FSLS+GLG1+FMLM1+FμLμ1=GLG2+FMLM2+FμLμ2，LG2=LM2=Lμ2=cosγcosβRsinθ.（20）

部分學(xué)者遵循韓其為對拖曳力力臂LD的取值，取LD=a+Rcosθ；a為拖曳力到顆粒中心的垂直距離，當(dāng)水流流速為拋物線分布時，a=R/3.實(shí)際上，由于水流速度與泥沙顆粒相對位置的隨機(jī)變化，拖曳力在泥沙顆粒上的作用點(diǎn)也不斷變化.本文采取王愉樂做法，通過對水流積分得到拖曳力作用點(diǎn)的位置，并假設(shè)所研究泥沙顆粒與上游顆粒接觸點(diǎn)位于河床基準(zhǔn)面上，該接觸點(diǎn)的水流速度為零.

在明渠中，時均流速U與摩阻流速u的關(guān)系如式（21）所示：

U=5.75lg30.2yxKSu*，（21）

式中：y為底部流速作用點(diǎn)的位置；x為關(guān)于黏性底層影響的修正系數(shù)，對于無黏性泥沙顆粒，x=1；K[KG-*3]S 為床面粗糙度，K[KG-*3]S=2D^[15，23].

yD=∫^H105.75u*lg30.2yxKS·ydy/

∫^H105.75u*lg30.2yxKSdy，H1=2R－=R+Rcosα，（23）

式中，H1為泥沙顆粒在河床上的突出高度.

由圖1可知：

yD+=a+R.（24）

將式（24）的a代入韓其為的拖曳力力臂公式中可得：

LD=yD－Rcosα+Rcosθ.（25）

將各力（對于非黏性顆粒，認(rèn)為附加下壓力FM與黏著力Fμ的值為零）及力臂代入式（17）可推求得到非黏性顆粒的瞬時起動流速u如式（26）所示：

u=φi，α，β，Δ′，′·ωD.（26）

φi，α，β，Δ′，′=ρs－ρρcosγcosβ2Δ′－Δ′²－sinβ1－Δ′－i1－e2Δ′－Δ′²185－95SymbolQC@′+342Δ′－Δ′²，（27）

ωD=gD4CL.（28）

泥沙顆粒的瞬時起動流速是指作用在泥沙顆粒底部的水流速度，而為了實(shí)用性多換算為垂線平均流速（沿水深的平均流速）.垂線平均流速U與摩阻流速uSymbolj@@的關(guān)系如式（29）所示：

Uu=6.5HD14+lgHD=φ.（29）

由式（21）和式（29）可得：

U=uφ=0.174φU.（30）

時均起動底速U與瞬時起動底速u的關(guān)系涉及到起動標(biāo)準(zhǔn).起動概率參考DOU^[24]中動及普動所對應(yīng)的起動概率ε1=0.0228、ε2=0.1590作為起動標(biāo)準(zhǔn)，由李林林等^[25]得到的時均起動底速U與瞬時起動底速u的關(guān)系，并參考韓其為給出的扁度函數(shù)f（λ）取值表，本文取f（λ）=1.59后可得：

U=0.914u，ε1=0.0228，

U=1.130u，ε2=0.1590.（31）

最終可得垂線平均流速U如式（32）所示：

U=kεuHD14+lgHD，（32）

式中：k（ε）和起動概率相對應(yīng)，k（0.0228）=1.034，k（0.1590）=1.278，H為水深.

4 起動流速公式驗(yàn)證及敏感性分析

針對垂線平均流速公式（32），選擇中動起動標(biāo)準(zhǔn)，采用文獻(xiàn)[10]中武漢水利電力大學(xué)在循環(huán)變坡水槽中的若干種下限固定在0.05 mm、上限變動范圍為25～12 mm，且級配連續(xù)的非均勻沙分級起動實(shí)測數(shù)據(jù)進(jìn)行驗(yàn)證.其水深變化范圍4.4～15.2 cm，比降變化范圍0.04～2.由比降可獲知縱坡坡度β，而滲透比降i和橫坡坡度α1在文獻(xiàn)[10]中不作為研究變量，始終為0.將所研究泥沙顆粒的直徑與該組泥沙顆粒的平均直徑代入式（3）、（4）和（5）中，即可求取相對隱蔽度1的數(shù)學(xué)期望值.最后通過式（32）得到垂線平均流速的計算值，并與文獻(xiàn)[15]和文獻(xiàn)[10]中的數(shù)據(jù)進(jìn)行對比，結(jié)果如表1和圖3所示.

從表1可以看出，本文大多數(shù)計算值與實(shí)測值差異不超過5%，有極少數(shù)達(dá)到11.從圖3可以看出，本文計算值與實(shí)測值較為接近.通過計算3個公式的可決系數(shù)R²和平均絕對誤差，得到如表2所列的結(jié)果，發(fā)現(xiàn)本文公式擬合效果較好.

垂線平均流速對縱坡坡度和滲透比降的敏感性分析結(jié)果如圖4所示.

在平均粒徑為0.736 mm的粒徑組中，對于粒徑為0.55 mm的泥沙顆粒，在水深為6.6 cm、滲透比降變化范圍為0.1～0.6、縱坡坡度變化范圍為0的工況下，由垂

線平均流速公式（32）得到泥沙起動流速變化情況.該圖表明，滲透比降和縱坡坡度與泥沙起動流速都呈負(fù)相關(guān)關(guān)系.滲透比降越大，泥沙起動流速越?。豢v坡坡度越大，泥沙起動流速越小.這與吳巖所得到的結(jié)果一致.對于具有一定坡度的堰塞壩潰壩土石料而言，滲透比降越大，土石料的可蝕性

就越高，堰塞壩泄流槽破裂幾何形狀也越大，從而加劇壩體垮塌，破壞過程所經(jīng)歷的時間就越短.

5 結(jié)語

本文建立了考慮隱暴效應(yīng)的斜坡非均勻無黏性泥沙顆粒的垂線平均流速公式.通過武漢水利電力大學(xué)測驗(yàn)資料對其進(jìn)行驗(yàn)證，結(jié)果表明公式計算值精度較高；對公式進(jìn)行敏感性分析發(fā)現(xiàn)，滲透比降和縱坡坡度與泥沙起動流速都呈負(fù)相關(guān)關(guān)系，泥沙起動流速隨著滲透比降和縱坡坡度的增大而減小.但不均勻性強(qiáng)且存在一定河床坡度的泥沙起動流速實(shí)測值較少，公式是否適用于具有堰塞壩土石料特征的大比降非均勻泥沙起動流速的計算還需進(jìn)一步檢驗(yàn).

參考文獻(xiàn)：

[1] COSTA J E，SCHUSTER R L.The Formation and Failure of Natural Dams[J].Geological Society of America Bulletin，1988，100（7）：1054-1068.

[2] 趙天龍，馬廷森，付長靜，等.考慮顆粒暴露角的寬級配泥沙起動流速公式研究[J].中國農(nóng)村水利水電，2023，No.483（1）：8-16.

[3] 趙天龍，付長靜，馬廷森，等.漫壩水流條件下寬級配土石料沖刷特性試驗(yàn)研究[J].人民長江，2023，54（1）：212-219.

[4] 錢寧，萬兆惠.泥沙運(yùn)動力學(xué)[M].北京：科學(xué)出版社，1983：240-242.

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

[6] 謝立全，于玉貞.滲流作用下的岸坡泥沙起動條件[J].清華大學(xué)學(xué)報（自然科學(xué)版），2006（9）：1534-1537.

[7] 吳巖，白玉川.滲流作用下的岸坡非均勻沙起動[J].水利水電技術(shù)，2014，45（1）：136-138.

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

[9] 李保如.泥沙起動流速的計算方法[J].泥沙研究，1959（1）：71-77.

[10] 韓其為，何明民.泥沙起動規(guī)律及起動流速[M].北京：科學(xué)出版社，1999.

[11] 何文社，方鐸，楊具瑞，等.泥沙起動流速研究[J].水利學(xué)報，2002（10）：51-56.

[12] 楊奉廣，劉興年，黃爾，等.唐家山堰塞湖下游河床泥沙起動流速研究[J].四川大學(xué)學(xué)報（工程科學(xué)版），2009，41（3）：84-89.

[13] WU F C，CHOU Y J.Rolling and lifting probabilities for sediment entrainment[J].Journal of Hydraulic Engineering，2003，129（2）：110-119.

[14] DUAN J G，SCOTT S.Selective bed-load transport in Las Vegas Wash，a gravel-bed stream[J].Journal of Hydrology，2007，342（3-4）：320-330.

[15] 馬菲，韓其為，李大鳴.非均勻沙分組起動流速[J].天津大學(xué)學(xué)報，2010，43（11）：977-980.

[16] 白玉川，王鑫，曹永港.雙向暴露度影響下的非均勻大粒徑泥沙起動[J].中國科學(xué)（技術(shù)科學(xué)），2013，43（9）：1010-1019.

[17] 王愉樂，張根廣，周雙，等.均勻無黏性單顆粒泥沙起動切應(yīng)力的探討[J].泥沙研究，2020，45（1）：1-6.

[18] 麻卉港，張根廣，周雙，等.均勻沙起動概率試驗(yàn)研究[J].泥沙研究，2022，47（4）：8-14.

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

[20] 周雙，張根廣，梁宗祥，等.斜坡上均勻散粒體泥沙的起動流速研究[J].泥沙研究，2015（4）：7-13.

[21] DEY S.Threshold of Sediment motion on combined transverse and longitudinal sloping beds[J].Journal of Hydraulic Research，2003，41（4）：405-415.

[22] 鄧明楓，陳寧生，鄧虎，等.堰塞壩漫頂潰決過程與機(jī)理試驗(yàn)研究[J].成都理工大學(xué)學(xué)報（自然科學(xué)版），2011，38（3）：359-365.

[23] 白玉川，陳有華，韓其為.泥沙顆粒躍移運(yùn)動機(jī)理[J].天津大學(xué)學(xué)報，2012，45（3）：196-201.

[24] DOU G R.Incipient motion of sediment under currents[J].China Ocean Engineering，2000（4）：391-406.

[25] 李林林，張根廣，吳彰松，等.正負(fù)坡上均勻散粒體泥沙起動流速的研究[J].泥沙研究，2016（5）：54-59.

［責(zé)任編輯：李嵐 杜佳］

基金項(xiàng)目：國家自然科學(xué)青年基金項(xiàng)目（42107189）

作者簡介：李楠（1998-），男，貴州銅仁人，在讀碩士，研究方向?yàn)檠呷麎螡Q過程研究.E-mail：1762550404@qq.com.

*通信作者：廖海梅（1988-），女，廣西梧州人，副教授，研究方向?yàn)樗姽こ碳盀?zāi)害研究.E-mail：291408646@qq.com.