匡翠萍， 鄭宇華， 顧 杰， 馬丹青

(1. 同濟(jì)大學(xué) 土木工程學(xué)院， 上海 200092； 2. 上海海洋大學(xué) 海洋科學(xué)學(xué)院， 上海 201306)

泥沙顆粒團(tuán)沉速

匡翠萍1， 鄭宇華2， 顧 杰2， 馬丹青2

(1. 同濟(jì)大學(xué) 土木工程學(xué)院， 上海 200092； 2. 上海海洋大學(xué) 海洋科學(xué)學(xué)院， 上海 201306)

針對(duì)現(xiàn)有的泥沙沉速公式大多適用于單顆粒泥沙沉降的情形，基于前人對(duì)單顆粒泥沙沉速規(guī)律的研究成果，利用試驗(yàn)測(cè)得不同泥沙顆粒團(tuán)在水中沉降時(shí)的沉速數(shù)據(jù)，引入泥沙顆粒團(tuán)特征粒徑及附加粒徑2個(gè)因子，對(duì)多個(gè)廣泛運(yùn)用的單顆粒泥沙沉速公式進(jìn)行了修正.用修正后的公式計(jì)算得到的泥沙顆粒團(tuán)沉速與用試驗(yàn)測(cè)得的泥沙顆粒團(tuán)沉速較接近,可推廣運(yùn)用到計(jì)算泥沙顆粒團(tuán)在水中的沉降速率.

泥沙； 顆粒團(tuán)； 沉速； 特征粒徑； 附加粒徑

在天然河道中，泥沙沉速的大小反映了泥沙顆粒在與水體相互作用時(shí)對(duì)水體阻力的抗拒能力[1].泥沙沉速是標(biāo)志泥沙運(yùn)動(dòng)特征的一個(gè)重要物理量.在泥沙輸運(yùn)過程中，泥沙顆粒不是以單顆粒形式沉降，而是以部分顆?；蛉款w粒的形式成群下沉.群體沉降與單顆粒沉降有較大不同，需額外考慮顆粒間的相互作用以及水體中顆粒的濃度對(duì)水流特性的影響.本文根據(jù)試驗(yàn)測(cè)得的泥沙顆粒團(tuán)在水中沉降時(shí)的沉速數(shù)據(jù)，研究不同粒徑組成及不同泥沙顆粒團(tuán)體積時(shí)泥沙顆粒團(tuán)的沉降特性，并基于前人研究所得的泥沙沉速公式，引入泥沙顆粒團(tuán)特征粒徑及附加粒徑因子，使泥沙單顆粒沉速公式同樣適用于泥沙顆粒團(tuán)的沉速計(jì)算，并使計(jì)算精度大幅提高.

1 泥沙沉速的研究進(jìn)展

泥沙顆粒沉速公式的研究最早是從單顆粒圓球在靜止流體中的受力分析出發(fā)，利用力的平衡原理，得到如下極限沉速公式[2]：

(1)

式中：ωd為泥沙顆粒(球體)的沉速；Cd為阻力系數(shù)；γs為泥沙容重；γ為流體(水)容重；g為重力加速度；d為泥沙顆粒直徑.

1851年，Stokes[3]忽略Navier-Stokes方程中所有的慣性項(xiàng)，用流函數(shù)法從理論上推導(dǎo)出了層流區(qū)單顆粒球體的阻力系數(shù)，得到了著名的Stokes球體沉速公式.Oseen和Goldstein[4]在Stokes的分析基礎(chǔ)上，做了進(jìn)一步的推導(dǎo)，給出了在層流沉降時(shí)圓球形固體顆粒的阻力系數(shù).竇國仁[5]采用Oseen型和Newton型阻力，按阻力迭加原則獲得繞流阻力，通過引進(jìn)與雷諾數(shù)有關(guān)的分離角提出關(guān)于圓球繞流在過渡區(qū)阻力的理論分析，拓寬了Stokes的理論，進(jìn)而導(dǎo)出單顆粒泥沙沉速公式.張瑞瑾等[1,6]根據(jù)阻力疊加原理，將 Stokes 型滯性阻力和 Newton型形狀阻力線性組合來表示繞流阻力，建立了一個(gè)適用于層流區(qū)、過渡區(qū)和紊流區(qū)的泥沙沉速公式.沙玉清[7]通過對(duì)沉速公式的數(shù)學(xué)轉(zhuǎn)換證明沉降規(guī)律可用無因次沉速ω*與無因次粒徑D*的函數(shù)來表達(dá)，進(jìn)而根據(jù)實(shí)測(cè)資料用曲線擬合的最小二乘法獲得了過渡區(qū)的顆粒沉速公式.

泥沙沉速的試驗(yàn)分析從20世紀(jì)70年代末開始.王兆印等[8]利用沉降筒進(jìn)行泥沙懸浮液的沉降試驗(yàn)，推導(dǎo)出均勻粗顆粒在細(xì)顆粒稀懸浮液中的沉速公式.韓文亮等[9]利用自行設(shè)計(jì)的錐形沉降筒進(jìn)行非均勻沙動(dòng)水中沉降試驗(yàn)，提出了非均勻沙群體沉速公式.金文等[10]利用激光測(cè)量?jī)x器(PIV)對(duì)長江口泥沙水樣進(jìn)行沉降速度場(chǎng)測(cè)量，得到平均泥沙沉速.萬遠(yuǎn)揚(yáng)等[11]借助大型可溫控自動(dòng)攪拌沉降試驗(yàn)筒對(duì)長江口天然細(xì)顆粒泥沙的沉降機(jī)理展開了室內(nèi)試驗(yàn)，結(jié)果表明含沙量和粒徑均是影響細(xì)顆粒泥沙沉速的重要因子.

由于水流速度及水體含沙量的不同力學(xué)性質(zhì)，在研究泥沙顆粒團(tuán)沉速方面，基于單顆粒泥沙沉降公式擴(kuò)展應(yīng)用范圍往往需要進(jìn)行細(xì)致劃分，這樣細(xì)劃過程在實(shí)際應(yīng)用中會(huì)帶來一定的復(fù)雜性甚至存在較大的差距.在高濃度無黏性均勻沙顆粒團(tuán)沉速研究方面，被廣泛應(yīng)用的Richardson型公式[12]為：

(2)

式中：ωs，ω0分別為一定粒徑在懸液及在清水中的沉速；Sv為體積濃度；m為與雷諾數(shù)有關(guān)的待定系數(shù).基于式(2)，大多研究圍繞含沙量建立了類似的公式，僅m取值不同.例如Richardson[12]的試驗(yàn)研究顯示，層流區(qū)m取值為4.65；錢寧等[4]認(rèn)為m是一個(gè)與雷諾數(shù)有關(guān)的變數(shù)，隨著雷諾數(shù)的提高而減少，最小為2.25；王兆印等[13]的試驗(yàn)表明m為7.0～7.5.由此可見，m值變化很大，使得式(2)的應(yīng)用受到很大的限制.張耀哲[14]通過將非黏性均勻泥沙群體沉降運(yùn)動(dòng)概化為上升水流在懸浮顆粒層內(nèi)的運(yùn)動(dòng)，建立了層流區(qū)和紊流區(qū)均勻沙群體沉速公式，但泥沙下沉的重力特征及黏性作用與水流上升特性并不完全一致.封光寅等[15]在研究泥沙群體顆粒平均粒徑及平均沉速時(shí)認(rèn)為采用幾何平均值計(jì)算平均粒徑及平均沉速會(huì)帶來較大的誤差.費(fèi)祥俊[16]在研究非均勻沙群體沉速時(shí)認(rèn)為從非均勻沙中找出平均粒徑或中值粒徑作為代表粒徑來計(jì)算群體沉速往往難以符合實(shí)際.

目前，對(duì)于泥沙群體沉速的研究和試驗(yàn)方法已經(jīng)很多，計(jì)算公式也很多，但這些泥沙群體沉速公式因個(gè)別參數(shù)不易確定或者難以統(tǒng)一從而使應(yīng)用中存在一定的局限性.因此，結(jié)合現(xiàn)有的單顆粒泥沙沉速公式及試驗(yàn)研究成果，引入特征粒徑及附加粒徑2個(gè)因子，并推導(dǎo)出特征粒徑與附加粒徑的表達(dá)式，使現(xiàn)有的單顆粒泥沙沉速公式同樣適用于泥沙顆粒團(tuán)的沉速計(jì)算.

2 試驗(yàn)研究

2.1 試驗(yàn)設(shè)計(jì)

試驗(yàn)是在香港理工大學(xué)水力學(xué)實(shí)驗(yàn)室完成的[17].試驗(yàn)水槽長度為12.0 m，寬度為0.3 m，深度為0.5 m.水深為0.3 m.曼寧系數(shù)近似為0.01.通過一個(gè)高精度的計(jì)時(shí)器測(cè)量時(shí)間.投放裝置由1個(gè)小型電動(dòng)機(jī)和1個(gè)半球形碗組成，用于將裝在半球型碗中一定體積的泥沙樣品投放到水槽中.將半球形碗安裝在水槽中間水面上方的一根固定的軸上.試驗(yàn)一旦開始，由同步電子啟動(dòng)器觸發(fā)電動(dòng)機(jī)、計(jì)時(shí)器和高速數(shù)字?jǐn)z影機(jī)同時(shí)工作，電動(dòng)機(jī)帶動(dòng)軸旋轉(zhuǎn)180°，將碗里的泥沙垂直倒入水槽中，計(jì)時(shí)器開始計(jì)時(shí)，裝在水槽邊上的數(shù)字?jǐn)z影機(jī)開始連續(xù)拍攝泥沙顆粒團(tuán)在水中的運(yùn)動(dòng)情況.試驗(yàn)示意圖見圖1，圖中H0為水深，t為泥沙在水體中的沉降時(shí)間.

圖1 試驗(yàn)示意Fig.1 The schematic diagram of experiment

所用的泥沙取自香港海域海床，將泥沙置于實(shí)驗(yàn)室自然風(fēng)干，其中部分泥沙顆粒黏在一起.為使試驗(yàn)所用的泥沙為分散顆粒，將泥沙放于烘箱中烘干至70%，烘干后的泥沙置于實(shí)驗(yàn)室，冷卻后待用.由于沙樣的顆粒大小分布有些不均勻，用篩子進(jìn)行篩選，分別得到泥沙顆粒大小為0.060～0.118 cm，0.030～0.060 cm和0.015～0.030 cm的3組泥沙作為試驗(yàn)樣品.

試驗(yàn)時(shí)所用參數(shù)：初始體積為3～7 cm3，溫度T為25℃，水體運(yùn)動(dòng)黏滯系數(shù)ν為0.008 97 cm2·s-1，重力加速度g為980 cm·s-2.試驗(yàn)共分為9組, 其他主要試驗(yàn)參數(shù)及結(jié)果見表1.

表1 主要試驗(yàn)參數(shù)及結(jié)果Tab.1 The main experimental parameters and results

2.2 試驗(yàn)結(jié)果及分析

為了獲得準(zhǔn)確的試驗(yàn)數(shù)據(jù)，試驗(yàn)以0.25 s為時(shí)間間隔，從第1個(gè)0.25 s開始，量取泥沙顆粒團(tuán)在水中的擴(kuò)散寬度及泥沙顆粒團(tuán)垂向前鋒位置，由泥沙顆粒團(tuán)的擴(kuò)散寬度及前鋒位置(也就是泥沙顆粒團(tuán)的外包線)作一個(gè)橢圓，其中心可作為泥沙顆粒團(tuán)的中心位置，并以此計(jì)算泥沙顆粒團(tuán)的沉降速率.接著再量取第2個(gè)0.25 s的泥沙顆粒團(tuán)在水中的擴(kuò)散寬度及泥沙顆粒團(tuán)垂向前鋒位置，依次類推，直至泥沙顆粒團(tuán)前鋒碰到水槽底部為止.每組試驗(yàn)至少能獲得8個(gè)0.25 s的試驗(yàn)結(jié)果，對(duì)每一組所得是0.25 s間隔的泥沙顆粒團(tuán)沉降速率求平均，作為該試驗(yàn)組次的泥沙顆粒團(tuán)沉降速率.為了減少不確定因素對(duì)試驗(yàn)結(jié)果的影響，每組試驗(yàn)重復(fù)5次，再取其平均值作為最終的每組次泥沙顆粒團(tuán)的沉降速率.

表1第5列表示的是在不同試驗(yàn)條件下通過試驗(yàn)測(cè)得的泥沙顆粒團(tuán)在水中沉降時(shí)的沉速試驗(yàn)值.從表1可以看出，泥沙顆粒團(tuán)在水中沉降時(shí)，隨著泥沙顆粒團(tuán)顆粒粒徑范圍及初始體積的增加，泥沙顆粒團(tuán)的沉速增加.

3 泥沙顆粒團(tuán)沉速公式的研究

3.1 泥沙沉降規(guī)律及常用沉速計(jì)算公式

泥沙顆粒在靜水中下沉?xí)r，泥沙顆粒的不同運(yùn)動(dòng)狀態(tài)有不同的沉降阻力規(guī)律，運(yùn)動(dòng)狀態(tài)與沙粒雷諾數(shù)Re有關(guān).當(dāng)Re<0.25時(shí)，泥沙顆粒的下沉屬于層流狀態(tài)；當(dāng)Re>850時(shí)，泥沙顆粒的下沉屬于紊流狀態(tài)；當(dāng)0.25< Re <850時(shí)，泥沙顆粒的下沉處于過渡狀態(tài)[1].用試驗(yàn)測(cè)得的泥沙顆粒團(tuán)的沉速計(jì)算Re，計(jì)算結(jié)果見表1第6列，因此，在試驗(yàn)條件下泥沙顆粒團(tuán)沉降性質(zhì)處于過渡狀態(tài).

目前，泥沙沉速計(jì)算的公式往往在層流區(qū)、過渡區(qū)和紊流區(qū)具有不同形式.表2列舉了幾個(gè)常用的泥沙沉速ω計(jì)算公式[1,5-7].

表2 常用的泥沙顆粒沉速計(jì)算公式Tab.2 Commonly used calculation formulas of settling velocity of sediment particle

注：① 1938年; ② 1954年.

3.2 過渡區(qū)泥沙沉速的研究

表3 泥沙顆粒團(tuán)沉速試驗(yàn)值與各公式計(jì)算值的對(duì)比Tab.3 The comparison of experimental values of settling velocity of sediment particle cloud with the calculated values by formulas

由于上述公式中均未考慮泥沙顆粒團(tuán)初始體積對(duì)沉速的影響，因此在相同粒徑范圍下，不同初始體積的泥沙顆粒團(tuán)用各公式計(jì)算得到的沉速相同.由表3可知：在相同泥沙顆粒粒徑范圍下，用各公式計(jì)算得到的泥沙顆粒團(tuán)的沉速均遠(yuǎn)小于試驗(yàn)測(cè)得的泥沙顆粒團(tuán)的沉速.因此，考慮到泥沙顆粒之間的相互影響，需要對(duì)沉速公式進(jìn)行修正.

3.3 引入特征粒徑及附加粒徑對(duì)張瑞瑾泥沙沉速公式的修正

文獻(xiàn)[4]提到，如果一團(tuán)泥沙在水體中下沉，周圍相當(dāng)遠(yuǎn)的范圍內(nèi)都是清水，則它們的下沉速度要比單顆粒泥沙在無窮水體中下沉速率大.文獻(xiàn)[19]中也提到，用粒徑計(jì)法分析粒徑級(jí)配時(shí)，由于泥沙群團(tuán)沉降，使得沉速偏大，粒徑級(jí)配偏粗.因此，結(jié)合眾多學(xué)者對(duì)泥沙沉速的研究成果[1-22]，考慮到張瑞瑾的泥沙沉速公式為三區(qū)統(tǒng)一公式，對(duì)張瑞瑾公式中的粒徑項(xiàng)進(jìn)行修正，建立的表達(dá)式如下：

(3)

式中：D′為泥沙顆粒團(tuán)特征粒徑，cm，它代表泥沙顆粒團(tuán)沉降時(shí)的平均特性；d′為泥沙顆粒團(tuán)中泥沙顆粒的算術(shù)平均粒徑，cm；f(d′)為附加粒徑，cm，表明特征粒徑與平均粒徑的關(guān)系.

將D′代入表2中張瑞瑾公式，得到泥沙顆粒團(tuán)的沉速ω′公式如下：

(4)

將試驗(yàn)測(cè)得的前6組泥沙顆粒團(tuán)的沉速試驗(yàn)值代入式(4)，計(jì)算得到D′的值，再利用表達(dá)式(3)，計(jì)算得到附加粒徑f(d′)的值，計(jì)算結(jié)果見表4.

表4 泥沙顆粒團(tuán)特征粒徑及附加粒徑Tab.4 The characteristic partical size and additional partical size of sediment particle cloud

從表4的計(jì)算結(jié)果可知，在不同泥沙顆粒團(tuán)粒徑組成及不同泥沙顆粒團(tuán)體積下，計(jì)算得到的泥沙顆粒團(tuán)的特征粒徑值及附加粒徑值均不同.由于泥沙顆粒中值粒徑的變化趨勢(shì)與水流紊動(dòng)強(qiáng)度類似，即粒徑自水面向下呈遞增趨勢(shì)，沉速呈指數(shù)形式下降[9]，因此為了研究泥沙顆粒粒徑與泥沙團(tuán)沉速的關(guān)系，設(shè)附加粒徑方程為f(d′)=ad′nx，其中，a為附加粒徑與平均粒徑的正比系數(shù)，x為附加粒徑與泥沙顆粒個(gè)數(shù)所呈指數(shù)關(guān)系的指數(shù)，n為單位體積中泥沙顆粒的個(gè)數(shù).

用試驗(yàn)所得數(shù)據(jù)的前6組代入方程計(jì)算分析，可推出附加粒徑的表達(dá)式如下：

(5)

用試驗(yàn)測(cè)得的后3組泥沙顆粒團(tuán)的沉速試驗(yàn)值對(duì)所推導(dǎo)的泥沙沉速公式(4)進(jìn)行驗(yàn)證，驗(yàn)證結(jié)果見表5.

表5 泥沙顆粒團(tuán)沉速試驗(yàn)值與式(4)計(jì)算值的對(duì)比

Tab.5 The comparison of experimental values of settling velocity with the calculated values by formula (4) for sediment particle cloud

組號(hào)ωT/(cm·s-1)ωC/cm·s-1η/%712．6712．23-3．47810．6711．305．9099．6710．003．41

從表5可知，用修正后的泥沙沉速公式(4)計(jì)算得到的泥沙顆粒團(tuán)的沉速與試驗(yàn)測(cè)得的泥沙顆粒團(tuán)的沉速進(jìn)行比較，最大誤差為5.90%，小于10%，計(jì)算精度較高，即張瑞瑾的泥沙沉速公式可寫成如下形式：

(6)

4 不同公式與附加粒徑的關(guān)系

4.1 附加粒徑適用性分析

進(jìn)一步研究發(fā)現(xiàn)，D′可適于所有泥沙顆粒團(tuán)的沉速計(jì)算，將其代入表2中過渡區(qū)的經(jīng)驗(yàn)公式進(jìn)行計(jì)算，結(jié)果見表6.

表6 泥沙顆粒團(tuán)沉速試驗(yàn)值與特征粒徑計(jì)算沉速值的對(duì)比Tab.6 The comparison of experimental values of settling velocity of sediment particle cloud with the calculated values by characteristic particle size

表6中，4，5，6這3組試驗(yàn)數(shù)據(jù)所有公式中計(jì)算值與試驗(yàn)值的誤差均偏高.泥沙粒徑范圍分別為0.060～0.118 cm的1，2，3組和0.015～0.030 cm的7，8，9組的試驗(yàn)結(jié)果均體現(xiàn)了泥沙顆粒團(tuán)在沉降過程中隨著初始體積的增加沉速也相應(yīng)增加的特征，而泥沙粒徑范圍為0.030～0.060 cm的4，5，6組中卻出現(xiàn)了隨初始體積的增加泥沙沉速反而減小的情況，說明4，5，6組試驗(yàn)結(jié)果存在一定的誤差.

除去4，5，6組計(jì)算結(jié)果，用沙玉清、武漢水利電力學(xué)院、張瑞瑾和Van Rijn的公式計(jì)算的結(jié)果精度都較高，而竇國仁和岡恰洛夫公式的計(jì)算結(jié)果與試驗(yàn)值相比存在較大誤差，需要對(duì)f(d)再進(jìn)行修正.

4.2 對(duì)岡恰洛夫公式的分析

早期的岡洽洛夫公式①僅表達(dá)ω和d呈線性關(guān)系，計(jì)算精度不高，經(jīng)反復(fù)計(jì)算分析后，附加粒徑的修正表達(dá)式如下：

(7)

用修正后的附加粒徑計(jì)算式計(jì)算D′，并代入岡恰洛夫過渡區(qū)的2個(gè)公式中，計(jì)算結(jié)果見表7.

表7 泥沙顆粒團(tuán)沉速試驗(yàn)值與用岡恰洛夫公式

修正后的計(jì)算值的對(duì)比

Tab.7 The comparison of experimental values of settling velocity of sediment particle cloud with the calculated values by modified Concharov’s formulas

組號(hào)ωT/(cm·s-1)岡恰洛夫公式①岡恰洛夫公式②ωC/(cm·s-1)η/%ωC/(cm·s-1)η/%114．0014．845．6613．88-0．86213．6713．18-3．7213．01-5．07312．6712．740．5511．92-6．29413．3311．98-11．2711．92-11．83513．6710．71-27．6410．73-27．40611．679．14-27．689．26-26．03712．6713．506．1512．730．47810．6711．083．7011．124．0599．679．65-0．219．11-6．15

由表7可知，除去4，5，6這3組情況，用修正后的附加粒徑f(d′)計(jì)算得到的泥沙沉速精度較高，岡恰洛夫公式①和②的沉速計(jì)算值與試驗(yàn)值的相對(duì)誤差的絕對(duì)值均小于10 %.

4.3 對(duì)竇國仁公式的分析

竇國仁在Stokes球體沉降公式的基礎(chǔ)上引入分離角的概念，且認(rèn)為公式中的系數(shù)不應(yīng)該是常數(shù)，而是雷諾數(shù)的函數(shù)，這使得公式的計(jì)算更具合理性，但同時(shí)也增加了公式使用的復(fù)雜性.研究發(fā)現(xiàn)，在不同的粒徑范圍內(nèi)，附加粒徑的修正值會(huì)隨著實(shí)際泥沙顆粒粒徑的變化而變化，無法建立一個(gè)各參數(shù)相同的統(tǒng)一表達(dá)式，因此需要根據(jù)不同的粒徑范圍對(duì)附加粒徑進(jìn)行分別修正.經(jīng)反復(fù)計(jì)算分析后，所得結(jié)果見表8.

表8 泥沙顆粒團(tuán)沉速試驗(yàn)值與用竇國仁公式修正后

計(jì)算值的對(duì)比及修正表達(dá)式

Tab.8 The comparison of calculated values by modified Dou Guoren’s formula with experimental values of settling velocity of sediment particle cloud and corresponding modified expressions

組號(hào)ωT/(cm·s-1)修正后竇國仁公式ωCη/%修正表達(dá)式114．0015．007．16213．6713．48-1．37312．6711．54-8．96413．3314．186．35513．6712．75-6．71611．6710．89-6．70712．6711．96-5．60810．6710．781．0599．679．22-4．61f1(d′)=0．006d′n0．722f2(d′)=0．006d′n0．650f3(d′)=0．006d′n0．592

由表8可知，由于每一組附加粒徑表達(dá)式中都選用一個(gè)不同的指數(shù)，公式修正后的各組計(jì)算值與試驗(yàn)值比較，相對(duì)誤差的絕對(duì)值均小于10%，計(jì)算精度大幅提高.對(duì)每一組附加粒徑表達(dá)式中的指數(shù)進(jìn)行進(jìn)一步分析可知，隨著粒徑范圍的增大指數(shù)逐級(jí)遞增，據(jù)此可以建立指數(shù)x與粒徑d′的關(guān)系：

(8)

則應(yīng)用竇國仁公式計(jì)算泥沙顆粒團(tuán)沉速時(shí)，其附加粒徑項(xiàng)為：

(9)

上述對(duì)各公式附加粒徑的修正說明不同的附加粒徑表達(dá)式對(duì)應(yīng)的各公式中單顆粒泥沙沉速與顆粒團(tuán)泥沙沉速的關(guān)系不同.雖然本次試驗(yàn)粒徑范圍在0.030～0.118 cm之間，但由于沙玉清、武漢水利電力學(xué)院、張瑞瑾、Van Rijn和岡恰洛夫的附加粒徑表達(dá)式中x與n有關(guān)，因此，理論上上述修正結(jié)果能應(yīng)用于大多粒徑范圍的計(jì)算，且指數(shù)x取值范圍為0.469～0.486；而應(yīng)用竇國仁公式其附加粒徑表達(dá)式中由于建立了x與d′的關(guān)系，在粒徑0.030～0.118 cm范圍之外是否適用，還有待進(jìn)一步的試驗(yàn)驗(yàn)證.

5 結(jié)語

利用試驗(yàn)測(cè)得的泥沙顆粒團(tuán)沉速數(shù)據(jù)對(duì)常用泥沙沉速公式進(jìn)行了研究，通過引入泥沙顆粒團(tuán)的特征粒徑及附加粒徑，推導(dǎo)出特征粒徑與附加粒徑的計(jì)算表達(dá)式，并對(duì)本文所列的泥沙沉速公式中泥沙的粒徑進(jìn)行了修正.用本文特征粒徑代入各泥沙沉速公式中計(jì)算，得到的泥沙顆粒團(tuán)的沉速與試驗(yàn)測(cè)得的泥沙顆粒團(tuán)的沉速基本一致，計(jì)算精度較高.

[1] 張瑞瑾,謝鑒衡,陳文彪. 河流動(dòng)力學(xué)[M]. 武漢:武漢大學(xué)出版社, 2007.

ZHANG Ruijin, XIE Jianheng, CHEN Wenbiao. Fluvial dynamics[M]. Wuhan: Wuhan University Press, 2007.

[2] 趙德招,陳立,周銀軍,等. 單顆粒泥沙沉速公式的對(duì)比研究[J].人民黃河, 2009, 31(1): 36.

ZHAO Dezhao,CHEN Li,ZHOU Yinjun,etal. The study on comparison of single sediment settling velocity formulas[J]. Yellow River, 2009, 31(1): 36.

[3] Stokes G G. On the effect of the internal friction of fluids on the motion of pendulums[J]. Transactions of Cambridge Philosophical Society, 1851, 9(2): 8.

[4] 錢寧,萬兆惠. 泥沙運(yùn)動(dòng)力學(xué)[M].北京:科學(xué)出版社, 2003.

QIAN Ning, WAN Zhaohui. Mechanics of sediment transport[M]. Beijing: Science Press, 2003.

[5] 竇國仁. 泥沙運(yùn)動(dòng)理論[M]. 南京:南京水利科學(xué)研究所, 1963.

DOU Guoren. Theory of sediment transport[M]. Nanjing: Nanjing Hydraulic Research Institute, 1963.

[6] 張瑞瑾. 河流泥沙動(dòng)力學(xué)[M]. 北京:中國水利水電出版社, 1998.

ZHANG Ruijin. River sediment dynamics[M]. Beijing: China Water & Power Press, 1998.

[7] 沙玉清. 泥沙運(yùn)動(dòng)學(xué)引論[M]. 北京:中國工業(yè)出版社, 1965.

SHA Yuqing. Introduction of sediment dynamics[M]. Beijing: China Industry Press, 1965.

[8] 王兆印,錢寧.高濃度泥沙懸浮液物理特性的實(shí)驗(yàn)研究[J].水利學(xué)報(bào),1984 (4):1.

WANG Zhaoyin,QIAN Ning.Experimental study on the physical properties of sediment suspensions with hyperconcentration[J].Journal of Hydraulic Engineering,1984 (4):1.

[9] 韓文亮,郜國明,惠遇甲.動(dòng)水中非均勻沙沉降規(guī)律的試驗(yàn)研究[J].泥沙研究,1998 (3):62.

HAN Wenliang, HAO Guoming, HUI Yujia. Experimental study on settlement law of non uniform sediment in dynamic water[J].Journal of Sediment Research, 1998(3):62.

[10] 金文，王增道. PIV直接測(cè)量泥沙沉速試驗(yàn)研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A), 2005, 20(1):19.

JIN Wen, WANG Zengdao.Direct settling velocity detection of sediment using PIV[J]. Chinese Journal of Hydrodynamics(A),2005, 20(1):19.

[11] 萬遠(yuǎn)揚(yáng),吳華林,顧峰峰.長江口細(xì)顆粒泥沙沉降速度室內(nèi)試驗(yàn)研究[J].人民長江,2014,45(1):98.

WAN Yuanyang,WU Hualin,GU Fengfeng. Experimental study on settling velocity of fine sediment in Changjiang River Estuary[J].Yangtze River, 2014, 45(1):98.

[12] Richardson J F. Sedimentation and fluidisation of fine particles and two component mixtures of solids[J]. Transactions of the Institution of Chemical Engineers, 1961, 39(5): 348.

[13] 王兆印，宋振琪. 非均勻非恒定流中泥沙運(yùn)動(dòng)規(guī)律初探[J]. 水利學(xué)報(bào), 1997(6): 2.

WANG Zhaoyin, SONG Zhenqi. A preliminary study on the movement of sediment in the non uniform unsteady flow[J]. Journal of Hydraulic Engineering, 1997(6): 2.

[14] 張耀哲.非黏性均勻沙群體沉降運(yùn)動(dòng)沉速公式的研究[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)：自然科學(xué)版,2006,34(4)：117.

ZHANG Yaozhe. Study on theoretic formula of sedimentation velocity in the sediment group settlement zone[J]. Journal of Northwest A & F University: Natural Science Edition,2006,34(4)：117.

[15] 封光寅,章厚玉,張孝軍,等.泥沙群體顆粒平均粒徑及平均沉速計(jì)算方法的修正[J].南水北調(diào)與水利科技,2003,6(1):36.

FENG Guangyin, ZHANG Houyu, ZHANG Xiaojun,etal. Silt calculation method & formula revision of average particle size and setting velocity[J].South-to-North Water Transfers and Water Science & Technology,2003,6(1):36.

[16] 費(fèi)俊祥.泥沙的群體沉降——兩種典型情況下非均勻沙沉速計(jì)算[J].泥沙研究,1992(3):11.

FEI Xiangjun. Settling of particle in suspension—Two typical cases in calculating the settling velocity of nonuniform sediment[J].Journal of Sediment research,1992(3):11.

[17] Gu J. Dispersion of unsorted particle cloud in ambient cross flow[D]. Hong Kong: The Hong Kong Polytechnic University, 2004.

[18] Gu Jie, Li Chiwai, Yang Hong,etal. Numerical study on the mixing of unsorted sediment particles discharged into a cross-flow by Multiphase Particle-In-Cell(MP-PIC) method[J]. International Journal of Sediment Research, 2007,22(2):150.

[19] 蘇杭麗,盧永生. 粒徑計(jì)顆粒級(jí)配分析的泥沙群團(tuán)沉降規(guī)律[J].河海大學(xué)學(xué)報(bào), 1999, 27(1):34.

SU Hangli, LU Yongsheng. A theory of silt group settling for particle size analysisby particle size apparatus[J].Journal of Hohai University, 1999,27(1): 34.

[20] Cheng N S. Simplified settling velocity formula for sediment particle[J].Journal of Hydraulic Engineering, 1997,132(2):149.

[21] 邵學(xué)軍,王興奎. 河流動(dòng)力學(xué)概論[M].北京:清華大學(xué)出版社,2005.

SHAO Xujun, WANG Xingkui. Introduction of fluvial dynamics[M]. Beijing: Tsinghua University Press, 2005.

[22] 顧杰,黃靜,李志偉. 瞬時(shí)排放泥沙顆粒團(tuán)在橫流中對(duì)流擴(kuò)散特性的試驗(yàn)研究[J].水科學(xué)進(jìn)展,2008, 19(4):483.

GU Jie, HUANG Jing, LI Zhiwei. Experimental study on dispersion of instantaneous discharging sediment cloud in cross-flow[J].Advances in Water Science, 2008, 19(4):483.

Settling Velocity of Sediment Particle Clouds

KUANGCuiping1，ZHENGYuhua2,GUJie2,MADanqing2

(1. College of Civil Engineering, Tongji University, Shanghai 200092, China; 2. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China)

The most of previous sediment settling velocity formulas were suitable for single sediment particle. Based on the previous research results on single sediment particle settling velocity laws and the experimental data of sediment particle cloud during the settlement process in water body, the formulas for calculating the settling velocity of sediment particle cloud were proposed by introducing two factors, i.e., the characteristic particle size of sediment particle cloud and additional particle size, to correct the existing single sediment particle settling velocity formulas. The settling velocity values of sediment particle cloud calculated by the corrected formulas were very close to the experimentally measured settling velocity values of sediment particle cloud, which indicates the corrected formula can be applied to calculate the settling velocity of sediment particle clouds.

sediment; particle cloud; settling velocity; characteristic particle size; additional particle size

2016-03-11

國家“九七三”重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2012CB957704)

匡翠萍(1966—)，女，教授，博士生導(dǎo)師，工學(xué)博士，主要研究方向?yàn)楹涌诤０豆こ蹋瓻-mail：cpkuang@#edu.cn

鄭宇華(1993—)，女，碩士生，主要研究方向?yàn)楹恿骱０秳?dòng)力學(xué)．E-mail: yuvaz@sina.com

TV142+.1

A