甘 甜，談廣鳴

(武漢大學(xué)水利水電學(xué)院，武漢 430072)

自然界中的泥沙是流體運動中受水流、風(fēng)力、波浪、冰川及重力作用移動后沉積下來的固體顆粒碎屑，而黏性泥沙在多沙河流中占有很大一部分比重，同時也大量存在于水庫、港灣、河口以及淤泥質(zhì)海岸中[1,2]。由于黏性泥沙顆粒細(xì)，重量輕等特性，其表面存在一定的物理化學(xué)作用使其在起動和沖刷特性上不同于粗顆粒沙，尤其是流速小時，黏性泥沙沉積在河床表面，經(jīng)過長時間的淤積固結(jié)，顆粒間水分排出，形成具有不同團聚結(jié)構(gòu)的微團。固結(jié)后微團的密度、黏結(jié)力、抗剪強度都發(fā)生了變化，增強了河床的抗沖性，表現(xiàn)出“淤積容易沖刷難”等特點[3]。泥沙起動問題是泥沙動力學(xué)中的基本問題，也是生產(chǎn)實踐中常涉及的問題，黏性泥沙的起動，與散體沙有著本質(zhì)的區(qū)別，黏性泥沙受黏結(jié)力的影響，固結(jié)后起動為微團形式。目前我國對于細(xì)顆粒泥沙受力及粗顆粒泥沙流速統(tǒng)一規(guī)律理論等的研究在國際上處于領(lǐng)先水平[4]，但由于黏性泥沙起動沖刷運動的復(fù)雜性，還有很多環(huán)節(jié)有待進一步研究。

天然粗顆粒泥沙動起動試驗一般采用明渠水槽和環(huán)形水槽作為試驗裝置，而對于固結(jié)歷時較長的黏性泥沙，顆粒間黏結(jié)力較強，臨界起動條件和相應(yīng)的摩阻流速變化很大，一般的明渠水槽控制實驗條件比較困難，也難以滿足較強的臨界起動條件，因此王軍、舒彩文、呂平等[5-7]均采用封閉矩形管道進行了黏性泥沙起動沖刷等系列試驗，雖然壓力管道水流條件滿足固結(jié)黏性泥沙起動的臨界水流條件，但管道內(nèi)流速分布與壓強分布均與明渠流不同，其結(jié)果是否能在天然河道明渠流中推廣，仍需仔細(xì)探討。

1 壓力管流與明渠流中固結(jié)黏性泥沙起動公式

土力學(xué)中黏性土邊坡滑動破壞與河流、水庫中黏性泥沙起動沖刷有一定的相似性，因此，為探究有壓管流與明渠流中黏性泥沙起動的異同，可借鑒土力學(xué)中黏性土邊坡穩(wěn)定性分析的研究成果[8]，針對2種模式下黏性泥沙的起動模式進行受力分析，運用黏性微團起動瞬間服從力矩平衡建立方程，推導(dǎo)黏性泥沙起動流速和起動切應(yīng)力公式。

1.1 壓力管流中固結(jié)黏性泥沙起動受力分析

為建立固結(jié)黏性泥沙起動流速及切應(yīng)力公式，假定受力微團為黏性泥沙顆粒組成的黏性微團，形狀為圓球體。假定床面淤積物為均勻分布，微團側(cè)面抗剪力忽略不計。起動模式為以0點為中心轉(zhuǎn)動起動(見圖1)。

由圖1可知，微團起動主要受重力、水流上舉力、拖曳力、黏結(jié)力、薄膜水附加壓力等，具體受力分析如下。

(1)微團脫離體的有效重力Ws?？紤]形狀影響：

Ws=a1(ρ0-ρ)gd3

(1)

式中：a1為微團的形狀的系數(shù)，與微團的幾何形狀有關(guān)；ρ0為微團干密度；ρ為水的密度；d為微團脫離體的直徑。

(2)水流拖曳力FD。水流流經(jīng)河床時，由于泥沙表面粗糙不平，水流和淤積物表面接觸后產(chǎn)生摩擦力以及由于微團頂部曲線分離造成形狀阻力，2者之和即為拖曳力，方向沿起動坡面，其表達式為：

(2)

式中：a2為面積系數(shù)，與微團迎水面幾何因素有關(guān)；CD為沿斜坡面阻力系數(shù)，與床面周圍的繞流流態(tài)有關(guān)；ρ為水的密度；u為管流中作用于淤積物表面的流速。

(3)水流上舉力FL。水流流動時，淤積物表面流速為水流運動速度，底部為微團與床面間滲透水的流動速度，比頂部流速小很多。根據(jù)伯努里定律，流速大小差異造成了壓力差，進而產(chǎn)生上舉力FL，方向垂直微團朝上。其表達式如下：

(3)

式中：a3為面積系數(shù)，與淤積物表面形狀有關(guān)；CL為上舉力系數(shù)。

黏性顆粒間的黏結(jié)力十分復(fù)雜，它一方面決定于泥沙的礦物組成、液體的粒子性質(zhì)和濃度、液體的溫度和pH值等，另一方面又與微團所形成的骨架有關(guān)，前者直接影響泥沙顆粒表面的薄膜水性質(zhì)和顆粒間的吸引力，后者則關(guān)系到整個系統(tǒng)的穩(wěn)定性。

(4)黏結(jié)力N。考慮細(xì)顆粒泥沙之間存在Van der Waals力，則黏結(jié)力為阻礙黏性泥沙起動的作用力。采用唐存本[9]的研究成果：

(4)

式中：ξ為黏結(jié)力系數(shù)；ρ′s為淤積物穩(wěn)定干密度(參考干密度)；ρ′/ρ′s為淤積物的相對干密度。

綜上，黏結(jié)力與水深無關(guān)，它宏觀地反映了微團脫離體在受拉時由于薄膜水接觸而引起的拉應(yīng)力。

對于淤積固結(jié)的黏性顆粒，與松散顆粒不同，靜水壓力無法通過顆粒之間的液體傳遞，使顆粒面上靜水壓力的平衡遭到破壞，產(chǎn)生薄膜水附加壓力。

(5)薄膜水附加壓力G與管道壓力P。薄膜水具有過渡性，薄膜水附加壓力是由于薄膜水不完全符合帕斯卡定理引起的。采用竇國仁[10,11]的研究成果：

G=ρgHωk

(6)

即：

(7)

式中：δ為顆粒薄膜水接觸厚度，取為0.21×10-7m；H為微團表面中心點處水深。

在壓力管道中，由于薄膜水附加壓力單向傳遞水壓力，故當(dāng)微團與床面淤積物接觸時，受到向下的力，在管道試驗中，假設(shè)薄膜水附加壓力傳遞管道水壓力，可簡化表示為：

(8)

(6)黏性微團正應(yīng)力Fn。微團下部河床施加給微團的支持力，方向垂直微團向上。當(dāng)微團脫離下部河床時Fn=0。

1.2 壓力管流中固結(jié)黏性泥沙起動公式推導(dǎo)

微團起動瞬間，假設(shè)其沿滑動面滑移脫離河床，水流的拖曳力和上舉力是微團起動的主要動力，而微團重力、黏結(jié)力、水壓力和薄膜水附加壓力則阻礙微團的起動，滿足繞0點的力矩平衡：

WsLWs+FDLFD+PLp+FnLFn+NLN+FLLFL+GLG=0

(9)

即：

(10)

將式(1)～式(8)代入式(10)得：

(11)

化簡得：

(12)

淤積固結(jié)條件下的黏性微團，阻礙其起動的主要原因為黏結(jié)力與薄膜水附加壓力，重力的影響很小[2]，故忽略式(12)中的重力項，進一步簡化公式得：

(13)

(14)

管流內(nèi)流速符合對數(shù)分布公式[12]：

(15)

式中：u*為摩阻流速；uy為距底部y的流速；Ks為床面粗糙度；χ為校正參數(shù)。

(16)

用微團底部作用流速作為固結(jié)黏性微團起動流速參考，從式(14)可以看出，影響固結(jié)黏性泥沙起動主要與水壓力和泥沙的相對干密度有關(guān)。

1.3 明渠流中黏性泥沙起動公式

對比文章上一節(jié)壓力管道中固結(jié)黏性泥沙起動受力分析，對黏性微團在明渠中的受力進行分析，得出起動流速公式。

由式(1)～式(4)可知，明渠中黏性泥沙起動所受重力、拖曳力、上舉力、黏結(jié)力均與壓力管流中相同，只有薄膜水附加壓力項不同。

根據(jù)竇國仁的研究成果，明渠流薄膜水附加壓力公式見式(7)。

假設(shè)在微團起動一瞬間，微團繞0點滑動起動，則滿足繞0點力矩平衡，同理忽略重力項，簡化得：

(18)

同樣的，在明渠中，微團底部作用流速與斷面平均流速為對數(shù)函數(shù)關(guān)系。從公式(18)可以看出，固結(jié)的黏性泥沙在明渠中的起動流速主要與水深與泥沙的相對干密度有關(guān)。

1.4 壓力管流與明渠流中黏性泥沙起動對比分析

竇國仁在研究細(xì)顆粒黏性泥沙起動流速時，考慮了黏結(jié)力與薄膜水附加壓力的影響，利用對數(shù)流速公式將床面作用流速轉(zhuǎn)化為斷面平均流速，得出公式(19):

(19)

對比1.2，1.3中推導(dǎo)出的公式可以看出，無論是壓力管流中的黏性泥沙起動流速公式(16)還是明渠中的黏性泥沙起動公式(18)，與竇國仁推導(dǎo)的公式形式相似，均與黏性微團的相對干密度和水深的1/2次方有關(guān)，故推導(dǎo)的流速公式合理。

對比式(15)與式(18)可以看出，影響?zhàn)ば晕F起動的實質(zhì)為壓力項與黏結(jié)力項共同作用的結(jié)果，2種情況下推導(dǎo)的底部流速公式形式相似，都是水深(壓力水頭)、微團粒徑、相對干密度的函數(shù)。其系數(shù)項略有不同，但并非影響?zhàn)ば阅嗌称饎拥闹饕蛩?。故從作用在微團表面的流速來看，2者差異不大，作用流速公式可簡寫為：

(20)

2 管流與明渠流對照試驗

為進一步驗證管道中黏性泥沙起動的結(jié)果是否能推廣到明渠中，需加入對照試驗。由于明渠水槽的水流條件弱，在明渠中進行淤積固結(jié)黏性泥沙起動試驗較為困難，洪大林[13]的實驗數(shù)據(jù)選取5組抗沖性較差的沙樣進行同組對照試驗。

試驗選擇5組抗沖性能較差的土樣，同時在矩形壓力管道與明渠水槽中進行起動試驗。測量明渠水槽中垂線流速分布，假設(shè)明渠壁面粗糙，根據(jù)明渠摩阻流速計算公式[式(15)]，由不同水深處流速比值求解Ks,進而解出u*，進而求出起動切應(yīng)力τ。表1為洪大林通過計算得出的明渠與壓力管道黏性泥沙的起動切應(yīng)力。將2者點繪至圖2，可見，數(shù)據(jù)點分布在45°線附近，即u*M≈u*G，τM≈τG。

表1 明渠與有壓管道中黏性泥沙起動切應(yīng)力 N/m2

圖2 黏性泥沙明渠與管流中起動切應(yīng)力對比Fig.2 The comparison of cohesive sediment incipient shear stress in open channel flow and pipe flow

同組對照的試驗結(jié)果同樣可以證明黏性泥沙在壓力管道中起動的結(jié)果可以推廣到一般明渠中使用。

3 結(jié) 語

通過對固結(jié)黏性泥沙起動的受力分析，推導(dǎo)出黏性泥沙起動公式并結(jié)合明渠黏性泥沙起動的對照試驗，得出以下結(jié)論。

(2)通過選取5組抗沖性較差的黏性土分別在明渠和管流中進行同組對照試驗，通過起動時不同深度水深，結(jié)合公式計算出明渠與管流中的起動切應(yīng)力，通過對比發(fā)現(xiàn)2者差異不大，故在小水深條件下，用管道進行試驗時，管道壓力項代替了明渠大水深條件下的靜水壓強和薄膜水附加壓力項，試驗結(jié)果基本可以代表沙樣在明渠水槽中的試驗結(jié)果。

文章中黏性泥沙起動公式中流速指標(biāo)為底部作用流速，若用斷面平均流速作為對比參數(shù)時仍需進行換算。

□

[1] 錢 寧，萬兆惠.泥沙運動力學(xué)[M].北京：科學(xué)出版社，2003.

[2] 張瑞瑾，謝鑒衡.河流泥沙動力學(xué)[J].北京：水利電力出版社，1998.

[3] 韓其為.水庫淤積[M].北京：科學(xué)出版社，2003.

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

[5] 王 軍.淤積固結(jié)條件下黏性泥沙起動與沖刷問題初步研究[D].武漢：武漢大學(xué)，2007.

[6] 呂 平，談廣鳴，王 軍.黏性泥沙淤后起動流速試驗研究[J].中國農(nóng)村水利水電，2008,(2)：56-58.

[7] 舒彩文，王 軍，談廣鳴.干密度對黏性淤積物起動和沖刷影響研究[J].武漢大學(xué)學(xué)報(工學(xué)版)，2007,(1)：25-28.

[8] 張?zhí)m丁.黏性泥沙起動流速的探討[J].水動力學(xué)研究與進展，2000，15(1)：82-88.

[9] 唐存本.泥沙起動規(guī)律[J].水利學(xué)報，1963,(2)：1-12.

[10] 張紅武.泥沙起動流速的統(tǒng)一公式[J].水利學(xué)報，2012，43(12)：1 378-2 396.

[11] 竇國仁.再論泥沙起動流速[J].泥沙研究，1999,(6)：1-8.

[12] 談廣鳴，舒彩文，陳一鳴.黏性泥沙淤積固結(jié)特性[M].北京：中國水利水電出版社,2014.

[13] 洪大林.黏性原狀土沖刷特性研究[D].南京：河海大學(xué)，2005.

[14] 劉 潔，劉潔玉，白玉川.黏性泥沙流變特性及其臨界起動研究[J].泥沙研究，2015,(6)：59-64.