基于隱蔽度的泥沙起動研究

基于隱蔽度的泥沙起動研究

1散體非均勻沙擊穿方法的改進沙砂巖的遷移研究已有200多年的歷史。在研究初期，主要集中在散砂上，如著名的“六方法則”和不規(guī)則的薩哈林水流量。這一時期對泥沙起動的處理方法一般分為3個步驟:①把泥沙顆粒概化為球體,忽略粒間作用力和粘結(jié)力;②通過對泥沙顆粒受力分析,建立滾動起動、滑動起動或跳躍起動方程,并推求泥沙近底作用起動流速;③采取某一水流流速分布來建立近底作用流速與垂線平均流速的關(guān)系,或者把近底作用流速轉(zhuǎn)化成無因次臨界起動切應(yīng)力。散體均勻沙的起動研究總體上還是比較成熟的,但是對散體非均勻沙的起動研究仍存在很多困難,至今起動公式仍沒有統(tǒng)一的形式。一般來講,散體非均勻沙起動的研究方法主要體現(xiàn)在4個方面:①取寬級配非均勻沙中某一級配泥沙做代表,直接代入散體均勻沙的起動(流速、切應(yīng)力、功率等,下同)公式;②引入一個附加力(阻力、動力及質(zhì)量力),建立泥沙起動平衡方程,推求出具有一定物理意義的起動公式(秦榮昱、張啟衛(wèi)等);③通過研究顆粒的絕對暴露度,建立被“增大”或“減小”后的等效粒徑,并代入起動平衡方程推求泥沙起動公式(劉興年等);④研究顆粒的相對暴露度,分析泥沙顆粒的起動力臂和阻力臂,進而推求出以相對暴露度為參數(shù)的起動公式(韓其為等)。從上述中間的2種研究方法中,我們不難得到寬級配非均勻沙的起動特性:粗顆粒較之同粒徑的均勻沙易于起動,細顆粒較之同粒徑的均勻沙難于起動,中等粒徑泥沙顆粒與同粒徑的均勻沙起動條件相當(dāng)。本文擬從3個方面來研究泥沙顆粒相對隱蔽度在散體沙起動中的應(yīng)用:①借鑒韓其為的研究思想,建立以相對隱蔽度為參數(shù)的散體沙近底作用起動流速公式;②結(jié)合非均勻沙顆粒相對隱蔽度的研究成果,數(shù)學(xué)上論證寬級配非均勻散體沙的起動特性;③分析起動標準及相對隱蔽度分布特性,數(shù)值上建立散體沙起動臨界切應(yīng)力與相對隱蔽度在不同起動標準下的對應(yīng)關(guān)系。以下所述隱蔽度若無特殊說明均表示泥沙顆粒相對隱蔽度。2泥沙粒徑的暴露度如圖1所示,參考前人的研究方法,將顆粒概化為二維“圓餅”狀。泥沙顆粒隱蔽度為床面上所研究表層泥沙顆粒的最低點到與下游顆粒接觸點之間沿床面的垂向距離(隱蔽長度)與所研究顆粒粒徑的比值,記為η。并定義1-η為泥沙顆粒的暴露度,記為ξ。研究表明:對于均勻沙來講,統(tǒng)計意義上隱蔽度的取值區(qū)間為[0.067,0.5];對于非均勻沙講,任意級配粒徑顆粒隱蔽度的期望的計算式為E(ηi)=nΣj=1ηi,jΡi,j=1i-1Σj=1C2ΝjC1Νi+C3Νi+nΣj=i+1C2ΝjC1Νi{i-1Σj=1C1ΝiC2Νj?12[12-√D2i+2DiDj2(Di+Dj)+Dj(Di+Dj)]+C3Νi12(0.067+0.5)+nΣj=i+1C1ΝiC2Νj12[12-√D2i+2DiDj2(Di+Dj)+12]}。(1)式中:ηi,j為顆粒Di在顆粒Dj影響的綜合隱蔽度;Pi,j為對應(yīng)組合的概率;Ni為顆粒Di在表層床沙中的粒數(shù)。3沙量方程3.1泥沙系數(shù)的建立泥沙起動流速(切應(yīng)力)是一個瞬時量,所研究的泥沙顆粒是表層床沙,也就是說所建立的泥沙起動公式僅僅適合“第一層”床沙。對于床面表層的散體沙,可以忽略粒間作用力和粘結(jié)力。水下泥沙顆粒所受力一般為水流拖拽力、水流升力(脈動上舉力)、水下重力。若忽略滲流梯度力的影響,可把泥沙受力、泥沙起動力矩和阻力矩的相關(guān)表達式列成表1。3.2cd3d2iu2ib2iidiidiidi型如圖2所示,根據(jù)表1中的表達式,我們可以建立粒徑為Di的泥沙顆粒滾動起動模型的力矩平衡方程Md=Mf,即CDαDD2iρu2ib2kξiξiDi=[α(γs-γ)D3i-CLαLD2iρu2ib2]√ηi-η2iDi。(2)其中,ξi+ηi=1,可以解得u2ib=2α[CDαDkξi√1-ηiηi+CLαL](ρs-ρρ)gDi。(3)3.3對特征高度的計算理論上講,可以從2個角度來解決該系數(shù)的取值問題:①經(jīng)驗計算,即借鑒韓其為對特征高度的經(jīng)驗取值,等價于本文利用的水流推力有效力臂,可以求出該系數(shù);②數(shù)值計算,即利用微積分計算出泥沙顆粒所受的水流力矩及推力,二者的比值為水流推力的有效力臂。3.3.1/2-+1/6根據(jù)韓其為的研究思想,假定水流推力作用點距球心上方1/6直徑處,這樣就得到kξ(1-η)=(1/2-η+1/6),(4)解得kξ=1-1/3(1-η)。代入均勻沙隱蔽度的期望0.2835,這樣得到均勻沙水流推力力臂修正系數(shù)kξ=0.5348。3.3.2顆粒-下游-非單次讓數(shù)算例2對于推移質(zhì)泥沙顆粒來講,河床上泥沙顆粒起動后,仍將“下落”(滾落、滑落、沉落)到河床的“凹槽”處,等待下次起動。當(dāng)所研究顆粒與下游顆粒發(fā)生碰撞(或向下游顆粒翻越)時,可以忽略下游顆粒對其隱蔽性的影響,該模型可近似概化為光滑平板上單顆粒起動(見圖3)。當(dāng)所研究顆?！跋侣洹钡较掠巍鞍疾邸碧帟r,下次起動時,要考慮下游顆粒對其隱蔽作用的影響(見圖4)。于是要分2種情況來研究水流推力力臂系數(shù):(1)泥沙顆?！懊娣e元”積分該模型考慮了顆粒間碰撞力的作用,此時可以忽略下游顆粒的影響,該模型可概化為光滑平板上單顆粒起動。根據(jù)對泥沙顆粒的“面積元”進行積分來求解其力矩及受力可得kξ=∫2R0yCDρu2b?2√2yR-y2dy∫2R02RCDρu2b?2√2yR-y2dy=∫2R0y1+2m√2yR-y2dy∫2R02Ry2m√2yR-y2dy。(5)這里的流速采用指數(shù)型流速分布公式,指數(shù)m取6,可以計算出kξ≈0.550。(2)動態(tài)力臂修正系數(shù)該模型考慮了顆粒間隱蔽度的作用,忽略水流在泥沙隱蔽長度范圍內(nèi)水流速度(孔隙流)的影響,同理可以計算出水流推力力臂系數(shù),kξ=∫(1-η)D0y1+2m√D24-[y-(12-η)D2]dy∫(1-η)D0(1-η)Dy2m√D24-[y-(12-η)D2]dy。(6)代入均勻沙隱蔽度的期望0.2835及流速分布指數(shù)m=6,可以計算出kξ=0.5091。耦合上述2種模型,可取泥沙顆粒水流推力力臂系數(shù)為式(5)與式(6)的平均值0.530?？梢园l(fā)現(xiàn)這2個角度得出的結(jié)論是基本一致的。當(dāng)推力力臂系數(shù)越大泥沙越容易起動,為提高泥沙起動條件的預(yù)測精度,可以確定水流推力力臂系數(shù)為0.530。對于寬級配非均勻沙來講,由于任意級配顆粒的隱蔽度受床沙級配特性的影響,無法在數(shù)值上構(gòu)建出任意級配粒徑顆粒的水流推力力臂修正系數(shù)。但可以定性分析:由于較大顆粒的隱蔽度期望較小,依照上述計算方法,對應(yīng)的水流推力力臂系數(shù)要略大于0.53;小顆粒對應(yīng)的水流推力力臂系數(shù)要略小于0.53。為統(tǒng)一研究寬級配非均勻沙的分級起動,可假定寬級配非均勻沙水流推力力臂系數(shù)的平均值與均勻沙水流推力系數(shù)相等。4偽完全隱蔽狀態(tài)研究泥沙隱蔽度的目的是為更好地解決泥沙起動流問題,以往人們對泥沙顆粒隱蔽度的分析都是在統(tǒng)計意義上進行的,即泥沙顆粒隱蔽度要小于0.5。事實上,床面泥沙確實存在一些極端情況,在分析文中所推求的泥沙顆粒近底作用起動流速公式之前,首先給出3點說明:(1)統(tǒng)計意義上來講,泥沙起動為泥沙輸移初始階段,河床泥沙在水流的作用下自動組織成穩(wěn)定的穩(wěn)定的床面,即表層泥沙趨向于起動,“次表層”泥沙處于穩(wěn)定狀態(tài)(瞬時意義)。(2)對于散體非均勻沙來講,大部分表層泥沙顆粒的隱蔽度都不大于0.5。當(dāng)某個顆粒的隱蔽度趨近于零時,則稱該顆粒處于“完全暴露”狀態(tài)(見圖5(a))。而對于極個別的表層沙顆粒來講,顆粒隱蔽度也有可能大于0.5,此時水流升力(脈動上舉力)起主要作用,該泥沙仍可以跳躍起動,則稱該顆粒處于“偽完全隱蔽”狀態(tài)(見圖5(b))。由于泥沙顆粒以這種起動(非輸移)形式的概率很小(約占3%),在研究泥沙顆粒隱蔽度時,往往忽略瞬時“偽完全隱蔽”狀態(tài)的影響。當(dāng)上游顆粒在所研究顆粒起動之前“下落”(包括沉落、滑落和滾落)到該顆粒上游的“空隙”處(見圖5(b))時,該顆粒則會因為被上下游2個顆?！翱ㄗ　倍荒芩矔r起動,則稱該顆粒處于“完全隱蔽”狀態(tài)(見圖5(c))。由于此狀態(tài)下“粒間作用力”不能忽略,所建立的平衡方程(2)不成立。為修正該方程,在此給出2個假定:其一,由于水流推力力矩作用不大,我們約定此時所研究顆粒隱蔽度趨于1(暴露度長度趨于零,則水流推力力臂近似為零);其二,把“粒間作用力”N一部分與水流升力(脈動上舉力)抵消,也即方程(2)中的上舉力系數(shù)趨近于零,另一部分附加到顆粒自身水下重力上,用顆粒體積系數(shù)α′修正。由于“粒間作用力”N大大增加了泥沙顆粒起動的阻力矩,使得該顆粒成為真正的不能瞬時起動的“次表層”泥沙。事實上,王興奎(1993)通過對顆粒臨界起動脈動分析,提出脈動上舉力不僅要大于顆粒水下重力,還要持續(xù)一段時間ΔT,泥沙顆粒才能起動。由于“空隙”的存在,假定上游顆粒趨向于會在Δt時間內(nèi)“下落”(包括沉落、滑落和滾落)“空隙”處,泥沙顆粒隱蔽度在大于0.5時是否仍能起動,取決于Δt與ΔT的大小關(guān)系:若上游顆?！跋侣洹睍r間Δt小于水流脈動上舉力持續(xù)時間ΔT,該顆粒在起動前會被上下游2個顆?！翱ㄗ　?即該顆粒處于“完全隱蔽”狀態(tài);否則當(dāng)“下落”時間Δt大于水流脈動上舉力持續(xù)時間ΔT時,這說明在上游顆?！跋侣洹敝?該顆粒早已被水流脈動上舉力揚起向下游運動。(3)對于散體均勻沙來講,由于顆粒間一般不會形成足夠大的“空隙”,無論是在理論意義還是現(xiàn)實意義上,床面表層泥沙的隱蔽度都不會大于0.5,也就是說:均勻沙顆粒理論上不存在處于“偽完全隱蔽”狀態(tài),其主要以滾動起動形式為主,泥沙起動平衡方程(2)完全適合床面表層的任何散體均勻沙顆粒。從動態(tài)角度上講,某些均勻沙顆?？梢源嬖谔幱凇巴耆┞丁睜顟B(tài)(發(fā)生碰撞或翻越下游顆粒時)。4.1cddk1-1-cll本文已經(jīng)推求出了泥沙近底作用起動流速公式(考慮函數(shù)的連續(xù)性,去掉了下標i):u2b=2α[CDαDkξ√1-ηη+CLαL](ρs-ρρ)gD。(7)式中,CD,CL,αD,αL,α,kξ為對應(yīng)的一些常數(shù)。下面分3種情況對式(7)進行極限分析。4.1.1微動滾落泥沙顆粒的初始組織結(jié)構(gòu)如圖6所示,由于此時所研究顆粒的隱蔽度趨于零,不難計算出limηi→0u2b=limηi→02α[CDαDkξ√1-ηη+CLαL]?(ρs-ρρ)gDi=0。(8)這說明當(dāng)某個顆粒處于“完全暴露”狀態(tài)時,微小的流速的水流都可以使該泥沙顆粒起動。也就是說:泥沙輸移初始階段之前,極個別“完全暴露”泥沙顆粒會通過微動滾落到“低凹”處(具有一定的隱蔽度)來自動組織成較穩(wěn)定的床面。這表明把泥沙顆粒的“極個別動”作為泥沙起動標準并不合適。4.1.2種類型的對比以“偽完全隱蔽”狀態(tài)的極端“極限跳躍”為例(完全不考慮水流推力力矩的作用),如圖7所示,個別表層泥沙顆粒的隱蔽度趨近1時,水流推力不能有效地促使該顆粒向下游轉(zhuǎn)動,但是由于上游顆粒并不能給予該顆粒有效的“壓力”,水流升力(脈動上舉力)仍能促使該顆粒以跳躍方式起動,此時上舉力系數(shù)為某一常數(shù)。若把隱蔽度趨于1代入式(7)可以得到limη→1u2b=limη→12α[CDαDkξ√1-ηη+CLαL]?(ρs-ρρ)gDi=2αCLαL(ρs-ρρ)gDi。(9)式(9)也就是被廣泛采用的泥沙躍移起動流速公式?？梢园l(fā)現(xiàn):該式的適用范圍是處于“偽完全隱蔽”狀態(tài)的泥沙顆粒(忽略了水流推力作用的極端情況),并且式(9)的躍移起動流速較之式(7)的滾動起動流速為大。事實上,這與Choi(2000)的試驗結(jié)論(滾動起動的拖拽力最小,跳躍起動的拖拽力最大)也是一致的。但是泥沙顆粒以這種形式起動的概率很小(充分紊流區(qū),約占3%)。4.1.3不同表層的泥沙顆粒如圖8所示(η=1),受“粒間作用力”N的影響,認為上舉力系數(shù)CL趨近與零。把此時的隱蔽度和上舉力系數(shù)代入式(7)得到limη=1CL→0u2b=limη→1CL→02α′[CDαDkξ√1-ηη+CLαL]?(ρs-ρρ)gDi=∞。(10)這說明“次表層”泥沙顆粒不會起動(除非在某一時刻能成為表層泥沙顆粒)。參考M.S.Yalin的觀點(在相鄰顆粒沒有他移以前,該顆粒不可能進入運動狀態(tài)),可認為“完全隱蔽”的泥沙顆粒處于“自鎖”狀態(tài)。綜上所述,我們可以得到一個論點:所研究泥沙顆粒在不同時間不同地點的隱蔽度變化可以直接影響該泥沙顆粒的運動狀態(tài),這也就能很好地解釋床面上某些顆粒在不同情況下為何出現(xiàn)“容易移動”、“難以起動”和“走走停停”的現(xiàn)象。4.2顆粒的擊穿特性根據(jù)泥沙起動流速的物理意義,可以認為u2b是一個關(guān)于隱蔽度的單值函數(shù),對η求導(dǎo)數(shù)后得d(ub2)dη=αCDαDkξ/ηη-η2[CDαDkξ1-ηη+CLαL]2?(ρs-ρρ)gD>0。(11)這表明ub2是一個關(guān)于隱蔽度的增函數(shù),也即,泥沙顆粒的起動流速將隨著該顆粒隱蔽度的增大而增大。由式(1)可知:對于寬級配散體非均勻沙來講,粗顆粒的隱蔽度較之同等粒徑下的均勻沙的為小;細顆粒的隱蔽度較之同等粒徑下的均勻沙的為大;中等顆粒的隱蔽度與同等粒徑下的均勻沙的相近。也即,粗顆粒促使較細顆粒的隱蔽度增大;細顆粒促使粗顆粒的隱蔽度減小;中等顆粒同時受到粗、細顆粒的影響,“增大”與“減小”效應(yīng)被相互抵消。若把式(1)代入到式(3)中,就可以得到結(jié)論:對于寬級配散體非均勻沙來講,粗顆粒的起動流速較之同等粒徑下的均勻沙的為小(易于起動);細顆粒的起動流速較之同等粒徑下的均勻沙的為大(難于起動);中等粒徑泥沙顆粒的起動流速與同等粒徑下的均勻沙的相當(dāng)。這就從理論上證明了寬級配散體非均勻沙的起動特性。另外,本文也為文獻中對“大顆?！薄ⅰ爸械阮w?！焙汀靶☆w粒”概念的定義提供了論據(jù)。5無因次擊穿切應(yīng)力線0.067褚君達對國內(nèi)外散體均勻沙起動條件的研究進行了匯總,為了方便與其他學(xué)者的研究成果進行比較,他取近底流速的特征高度為泥沙顆粒粒徑,即認為近底作用流速與摩阻流速的轉(zhuǎn)化關(guān)系為ub=8.5u*。若把該式代入式(7)可得u*2=2α72.25[CDαDkξ1-ηη+CLαL](ρs-ρρ)gD。(12)把CD=0.4,CL=0.1,kξ=0.53,αD=π/4,α=π/6代入式(12)得起動切應(yīng)力為θc=ρu*2(γs-γ)D=1/372.25[0.0531-ηη+0.025]。(13)文獻已分析出均勻沙隱蔽度的取值區(qū)間為[0.067,0.5]。對于床面上的散體沙來講,總有某些顆粒的隱蔽度很小,略大于0.067;也總存在某些顆粒的隱蔽度很大,略小于0.5;也有一些顆粒隱蔽度介于期望0.2835左右。把泥沙顆粒隱蔽度從0.067到0.5變化代入式(13),可以得到無因次起動切應(yīng)力(shields數(shù)),即θci=ρu*2(γs-γ)D=1/372.25[0.0531-ηiηi+0.025]∈[0.0277,0.0571]。(14)這也說明:對于泥沙顆粒隱蔽度不同的取值,可以得到對應(yīng)的泥沙起動條件。該區(qū)間基本涵蓋了褚君達對眾多學(xué)者的統(tǒng)計值,與錢寧的統(tǒng)計區(qū)間[0.03～0.06]也是一致的。由于泥沙顆粒在不同的情況下完全有可能呈現(xiàn)出不同的隱蔽度,也就對應(yīng)不同的起動條件,因此眾多學(xué)者的研究成果都具有一定的合理性。另外,根據(jù)克雷默(H.Kramer)對泥沙起動的判別標準:①弱動—床面上這里或那里有屈指可數(shù)的細顆粒處于運動狀態(tài);②中動—床面上有中等大小顆粒在運動,運動強度已無法計數(shù);③普動—各種大小的顆粒已投入運動,并持續(xù)地普及床面各處。床面上表層泥沙在沖刷前都被各自的“次表層”泥沙“支撐”,表層顆粒間互不影響,此時表層顆粒的隱蔽度較小,甚至某些顆粒的隱蔽度較小,很容易在較小水流作用下起動;隨著水流強度的增大