趙利安，王鐵力

（1.遼寧工程技術(shù)大學(xué)礦業(yè)學(xué)院，遼寧阜新 123000；2.中煤科工集團(tuán)武漢設(shè)計(jì)研究院有限公司管道輸煤所，武漢 430070）

伴有滑動(dòng)床的水砂漿體傾斜管道流動(dòng)研究

趙利安1，＊，王鐵力2

（1.遼寧工程技術(shù)大學(xué)礦業(yè)學(xué)院，遼寧阜新 123000；2.中煤科工集團(tuán)武漢設(shè)計(jì)研究院有限公司管道輸煤所，武漢 430070）

為進(jìn)一步研究伴有滑動(dòng)床的傾斜漿體管道砂顆粒輸送規(guī)律，應(yīng)用力平衡理論，研究了懸浮層、滑動(dòng)床層速度以及滑動(dòng)床厚度求解方法，由此給出了傾斜管道淤積臨界速度迭代求解方法。同時(shí)，應(yīng)用懸浮層剪應(yīng)力線(xiàn)性分布假設(shè)，并借助普朗特經(jīng)驗(yàn)公式，給出了滑動(dòng)床上方速度分布模型。粗砂管道水力輸送實(shí)驗(yàn)數(shù)據(jù)表明：（1）淤積臨界流速與顆粒粒徑和輸送濃度負(fù)相關(guān)，輸送濃度不變時(shí)，淤積臨界流速與管道傾角大小正相關(guān)。與模型計(jì)算結(jié)果對(duì)比表明淤積臨界速度計(jì)算值和實(shí)測(cè)值最大偏差為13%；（2）隨著管道傾角增加，最大速度點(diǎn)位置有向上偏移的趨勢(shì)。與計(jì)算結(jié)果對(duì)比表明上層速度分布理論模型計(jì)算值和實(shí)測(cè)值偏差不超過(guò)10%。結(jié)果表明力平衡理論和剪應(yīng)力線(xiàn)性分布假設(shè)能夠較好地用于伴有滑動(dòng)床的傾斜粗砂漿體管道淤積臨界速度和滑動(dòng)床上方速度分布預(yù)測(cè)研究。

沉降性漿體；淤積臨界流速；滑動(dòng)床；擴(kuò)散方程；速度分布

0 引 言

隨著漿體管道輸送技術(shù)的發(fā)展，越來(lái)越多的行業(yè)和部門(mén)要用到物料的管道輸送技術(shù)，如水力采煤、海洋采礦、地面金屬礦管道輸送、水庫(kù)河流清淤和地面管道輸煤等。在特定條件下，如地勢(shì)起伏不平，或者管道容易受到外界干擾時(shí)，就要涉及到傾斜管道漿體輸送技術(shù)。傾斜向上管道輸送中，淤積臨界流速和速度分布的研究對(duì)于傾斜漿體管道輸送的參數(shù)優(yōu)化、堵塞的防止以及輸送成本的降低都有重要意義。

近5年來(lái)，國(guó)內(nèi)外漿體管道輸送技術(shù)中，淤積臨界速度的研究主要集中在水平管道方面［1－15］。這些研究均著眼于水平管道流動(dòng)，未能考慮管道角度的因素。此外，上述關(guān)于淤積臨界流速的研究基本上是從單個(gè)顆粒受力考慮的，而不是從群體顆粒角度考慮。為了便于研究，本文中的淤積臨界速度指的是管道底部固體顆粒全部開(kāi)始滑動(dòng)時(shí)的漿體平均流速?；瑒?dòng)底床流動(dòng)是管道中顆粒由靜止轉(zhuǎn)入運(yùn)動(dòng)的過(guò)渡形式之一，與淤積臨界流速有密切的關(guān)聯(lián)，其在傾斜管道試驗(yàn)中被多次觀(guān)測(cè)到［16－17］。因而，要確定固體顆粒管道輸送的臨界淤積條件，一般需要對(duì)滑動(dòng)底床的運(yùn)動(dòng)進(jìn)行研究（實(shí)際管道輸送中，出現(xiàn)層移運(yùn)動(dòng)的可能性較?。Ｆ駷橹?，伴有滑動(dòng)床的傾斜管道漿體流動(dòng)方面的研究主要集中于阻力損失和床上方濃度分布方面［18－20］。而滑動(dòng)床上方速度分布的研究也未見(jiàn)公開(kāi)發(fā)表?？梢?jiàn)，從滑動(dòng)床運(yùn)動(dòng)角度考慮傾斜漿體管道淤積臨界速度和速度分布的研究國(guó)內(nèi)外尚無(wú)記載。而粗顆?；瑒?dòng)床輸送形式在工程中的應(yīng)用越來(lái)越多，因此，有必要開(kāi)展這方面的研究。

1 滑動(dòng)床厚度求解模型的推導(dǎo)

在傾斜管道中建立坐標(biāo)系，以沿著管道軸線(xiàn)傾斜向上方向?yàn)閤軸，以垂直于管道軸線(xiàn)向上方向?yàn)閥軸，垂直主流方向，距離管道底部y處顆粒體積濃度為Sv。眾多研究表明，在固液兩相流中，固體顆粒的濃度是隨高度的變化而變化的，流動(dòng)參數(shù)只在y方向發(fā)生變化［21－22］。根據(jù)紊流擴(kuò)散理論，可以得到［20］：

考慮邊界條件并對(duì)式（1）進(jìn)行積分，結(jié)果為：

對(duì)式（2）積分，可以求得懸浮層的平均濃度公式為：

式中：SV1為懸浮層（上層）的平均體積濃度；SVB為下層（滑動(dòng)層）的平均體積濃度；A1為管道中懸浮層所占的面積；εs為固體顆粒的擴(kuò)散系數(shù)，其計(jì)算方法見(jiàn)有關(guān)文獻(xiàn)［20］。

設(shè)傾斜管道內(nèi)漿體的平均體積濃度為SV；下層（滑動(dòng)層）的平均體積濃度為SV2（等于SVB）；下層斷面積為A2；上層的平均流速為U1；下層的平均流速為U2；P1和P2為長(zhǎng)度為l的管段兩邊的壓力，如圖1所示。

圖1 漿體傾斜管道流動(dòng)示意圖Fig.1 Schematic diagram of slurry pipe flow in inclined pipe

根據(jù)顆粒和液體質(zhì)量守恒，可以得到：

設(shè)所研究的漿體管道長(zhǎng)度為l，垂高為h；管道上部?jī)?nèi)壁與漿體之間的剪應(yīng)力為τw；滑動(dòng)床上部漿體與滑動(dòng)床之間的剪應(yīng)力為τ0；管道下部滑動(dòng)床內(nèi)部流體與管壁之間的剪應(yīng)力為τ2。對(duì)于滑動(dòng)床上部漿體，當(dāng)其處于均勻流動(dòng)時(shí)，根據(jù)力的平衡條件，如圖1所示，有：

式中：x1為圓弧ACB長(zhǎng)度；x0為弦AB的長(zhǎng)度，P為長(zhǎng)度為l為管段兩邊的壓力差（即P1－P2）。根據(jù)圖1，x1和x0可以分別表示為：x1＝D（π－θ），x0＝D sinθ。其中D為管道直徑，θ為滑動(dòng)床中心角。

式（6）可以變?yōu)椋?/p>

ρ1、τw和im的計(jì)算見(jiàn)文獻(xiàn)［23］。

式（7）中剪應(yīng)力τ0計(jì)算公式為：

式中：f12為懸浮層和滑動(dòng)床間界面的阻力系數(shù)，此處采用C.A.Shook和M.C.Roco提出的公式［24］：

式中：d為顆粒平均粒徑。

將式（3）、（4）、（5）和（7）聯(lián)合組成方程組求解，可以求解出滑動(dòng)床厚度hb，以及U1、U2和SV1。

為了驗(yàn)證通過(guò)式（3）、（4）、（5）和（7）聯(lián)合組成方程組求解hb、U1、U2和SV1迭代方法的收斂性，進(jìn)行了顆粒平均粒徑為1.71mm漿體的水平管道和傾角為30°傾斜管道的試算，輸送濃度CVd為0.21，漿體平均速度為5.5m／s。利用Matlab軟件編制了迭代法求解非線(xiàn)性方程組的程序，設(shè)置求精度為0.001，迭代20次后，所求解量與迭代次數(shù)n的關(guān)系如圖2所示。

圖2 迭代次數(shù)n與各變量的關(guān)系Fig.2 Relation between the iteration number and the variables

從圖2各變量隨迭代次數(shù)變化的曲線(xiàn)可以看出，計(jì)算結(jié)果收斂。這說(shuō)明，通過(guò)式（3）、（4）、（5）和（7）聯(lián)合組成方程組求解hb、U1、U2和SV1是可行的。

2 滑動(dòng)床上方速度分布模型的推導(dǎo)

漿體管道輸送時(shí)，上層漿體紊流流動(dòng)時(shí)，剪應(yīng)力一般呈現(xiàn)k型分布，如圖3所示。剪應(yīng)力τ0和上部剪應(yīng)力τw可以通過(guò)計(jì)算得到。根據(jù)圖3可知，漿體懸浮層流動(dòng)時(shí)，存在剪應(yīng)力為零的點(diǎn)，由對(duì)稱(chēng)性可知，零剪力點(diǎn)一定在斷面的中垂線(xiàn)上，且是最大流速u(mài)max的所在位置［25］。假定滑動(dòng)床上方垂向剪應(yīng)力呈線(xiàn)性分布，零剪應(yīng)力點(diǎn)在滑動(dòng)床面以上y1處。則根據(jù)圖3所示，可以得到：

圖3 伴有滑動(dòng)床的傾斜管道速度分布和剪應(yīng)力分布示意圖Fig.3 Schematic diagram of velocity distribution and shear stress distribution of inclined slurry pipe with sliding bed

最大速度可用以下各式求解［26］：

上面3式中，ρ為清水密度，u為清水的動(dòng)力粘性系數(shù)，u＊為剪切速度，計(jì)算式為：

求得最大速度umax后，利用任意點(diǎn)剪切應(yīng)力τr與普朗特半經(jīng)驗(yàn)公式得到的剪應(yīng)力相等，可以得到速度分布計(jì)算公式為：

式中：r為距離最大速度點(diǎn)的垂直距離；u＊1為當(dāng)?shù)丶羟兴俣?，k為卡門(mén)常數(shù)，取值為0.4；u＊1的計(jì)算公式為［26］：

式中：Sv為任意點(diǎn)體積濃度，S為固體顆粒比重，ρ為清水密度。Sv可以通過(guò)式（2）確定。

3 傾斜管道中臨界速度求解模型推導(dǎo)

以?xún)A斜管道中滑動(dòng)床為研究對(duì)象，如圖1所示，當(dāng)滑動(dòng)床處于受力平衡狀態(tài)時(shí)，設(shè)τ為滑動(dòng)床上方處于跳動(dòng)的顆粒對(duì)于滑動(dòng)床的碰撞力沿著管道軸方向的分量，N為碰撞力垂直于管軸的分量。使滑動(dòng)床保持受力平衡時(shí)的速度就是管道顆粒淤積的臨界速度。此時(shí)，使滑動(dòng)床移動(dòng)的動(dòng)力Fd為：

阻止滑動(dòng)床移動(dòng)的阻力Fz為：

式中：x2為圓弧AEB長(zhǎng)度；fs為滑動(dòng)床與管道壁間的摩擦系數(shù)；G為滑動(dòng)床水下的有效重力，其表達(dá)形式為：

通常，認(rèn)為τ和N很小，因而τ和N在此可以忽略。再考慮到τ2是滑動(dòng)床內(nèi)部流體與管壁之間的剪應(yīng)力，可認(rèn)為近似與上部漿體與管壁的剪應(yīng)力τw相等，則式（17）和（18）可以分別簡(jiǎn)化為：

因而淤積臨界速度具體求解方法步驟為：

（1）令漿體初始的平均速度Um（0）足夠大，使管道底部無(wú)靜止顆粒層存在；

（2）k為用于計(jì)算次數(shù)的整數(shù)，開(kāi)始k＝1；

（3）U（k）＝U（k－1）－ΔU，ΔU為一微小的速度增量；

mmmm

（4）將式（3）、（4）、（5）和（7）組成方程組，聯(lián)合求解出hb、U1、U2和SV1；

（5）利用式（21）和（22）分別計(jì)算滑動(dòng)床所受的動(dòng)力Fd和阻力Fz；

（6）若Fd＞Fz，則k＝k＋1，返回步驟（3），重復(fù)（3）～（6）各步，否則計(jì)算結(jié)束；

（7）最后得到的Um（k－1）即為淤積臨界流速度Udl。

4 實(shí)驗(yàn)研究

4.1 實(shí)驗(yàn)裝置總體描述

試驗(yàn)總體管路布置如圖4所示，管路總體長(zhǎng)度40m左右。為了減少測(cè)量誤差，試驗(yàn)中采用了等動(dòng)力探頭及樣品旁路管路裝置，如圖4右側(cè)所示［27］。為了能夠改變測(cè)量管路的傾角，采用電機(jī)7、滑輪8，將測(cè)試管段吊起，從而改變管路傾角。

圖4 試驗(yàn)管路布置及等動(dòng)力測(cè)量裝置布置圖Fig.4 Layout of test pipeline and isokinetic measuring pipe

4.2 采樣單元原理

根據(jù)伯努利方程，當(dāng)圖4中水頭損失h1等于h2時(shí)，管道中的漿體流速與采樣頭12中的流速相同，此時(shí)測(cè)量比較準(zhǔn)確。等動(dòng)力學(xué)采樣探頭12可以旋轉(zhuǎn)，因此可以在管道橫斷面內(nèi)任一點(diǎn)取樣。當(dāng)水頭損失h1等于h2時(shí)，閥門(mén)15和16之間的漿體樣本與管道內(nèi)漿體誤差最小。此時(shí)關(guān)閉閥門(mén)15和16，通過(guò)排泄閥17獲得一定體積的漿體樣本。對(duì)獲取的漿體樣本進(jìn)行處理，可以得到干燥的固體顆粒體積及顆粒級(jí)配，進(jìn)而得到漿體的各個(gè)測(cè)量點(diǎn)當(dāng)?shù)貪舛?。通過(guò)流量取樣閥4排出的漿體樣本，可以得到輸送濃度CVd的值。漿體各個(gè)測(cè)量點(diǎn)流動(dòng)速度可通過(guò)流量計(jì)18得到。在透明有機(jī)玻璃管段，采用JVC GZ－N1PAC高清閃存攝像機(jī)拍攝顆粒在管道內(nèi)的滑動(dòng)床運(yùn)動(dòng)狀態(tài)圖像。來(lái)自于漿體罐2的漿體通過(guò)漿體泵1加壓，進(jìn)入測(cè)試管段10。

試驗(yàn)中采用砂顆粒平均粒徑分別為0.4、0.6和1.71mm，顆粒的級(jí)配情況如圖5所示，3種顆粒均級(jí)配均為窄級(jí)配分布。粗砂顆粒密度為2650 kg／m3，形狀為近似球形。

圖5 顆粒級(jí)配曲線(xiàn)Fig.5 Map of particle size distribution

測(cè)試段管道內(nèi)徑為150mm，管道傾角分別為0°、10°、16°、20°、25°和30°。漿體濃度變化范圍為0.10～0.27。測(cè)試中將管道中心點(diǎn)O上、下各劃分為4個(gè)條帶，從而需要測(cè)定9個(gè)點(diǎn)的濃度值（假定各條帶內(nèi)顆粒濃度均勻）。透明有機(jī)玻璃測(cè)試管段的直徑為150mm，長(zhǎng)度為460mm，用來(lái)觀(guān)測(cè)不同角度及流動(dòng)狀態(tài)。管道內(nèi)漿體的平均體積濃度SV通過(guò)對(duì)濃度分布進(jìn)行積分而得到。

試驗(yàn)中，通過(guò)連續(xù)調(diào)整漿體平均速度Vm，觀(guān)測(cè)和記錄管道底部顆粒運(yùn)行狀態(tài)。當(dāng)管道底部滑動(dòng)床從連續(xù)運(yùn)動(dòng)狀態(tài)即將進(jìn)入靜止?fàn)顟B(tài)時(shí)，記錄此時(shí)對(duì)應(yīng)的平均速度Um，即為淤積臨界速度值。經(jīng)過(guò)試驗(yàn)確定，顆?；瑒?dòng)床內(nèi)濃度SVB近似為一常數(shù)，約0.5左右，固定床濃度大致為0.6左右。管道與粗砂顆粒間的滑動(dòng)摩擦系數(shù)fs為0.55左右。

4.3 圖像采集與分析

用高速攝影機(jī)測(cè)量漿體管道中不同運(yùn)動(dòng)狀態(tài)時(shí)，顆粒床的厚度和運(yùn)動(dòng)速度等信息。拍攝的圖像蘊(yùn)含顆粒運(yùn)動(dòng)信息。試驗(yàn)時(shí)，通過(guò)逐漸增大漿體流速，同時(shí)拍攝顆粒群體管道中運(yùn)行狀態(tài)。攝像機(jī)最高拍攝速度為6×105幀／s，鏡頭分辨率為2350pixel× 900pixel。攝像機(jī)通過(guò)usb接口將圖像一幅一幅地傳輸?shù)诫娔X中進(jìn)行信息分析，從而確定各個(gè)流速下是否存在滑動(dòng)床，以及床高度和速度等信息。

4.4 不同工況條件

水平管道中的淤積臨界流速度測(cè)試試驗(yàn)工況如表1所示；傾斜管道中測(cè)試不同傾角下的淤積臨界流速度以及淤積時(shí)床角度和管道傾角關(guān)系的試驗(yàn)工況如表2所示。

表1 水平管道流動(dòng)實(shí)驗(yàn)工況Table 1 Experimental conditions for horizontal pipe flow

表2 傾斜管道流動(dòng)實(shí)驗(yàn)工況Table 2 Experimental conditions for inclined pipe flow

試驗(yàn)中發(fā)現(xiàn)，漿體速度逐漸減低，當(dāng)形成滑動(dòng)床時(shí)，進(jìn)一步降低平均速度時(shí)，滑動(dòng)床的移動(dòng)變得不穩(wěn)定。從不穩(wěn)定滑動(dòng)床移動(dòng)，變動(dòng)到固定床形式流動(dòng)狀態(tài)的速度范圍大致為0.35m／s。

輸送濃度在0.20～0.22范圍變動(dòng)時(shí)，通過(guò)圖6所示的調(diào)節(jié)管道傾角的裝置使管道傾角逐漸增大時(shí)，傾角對(duì)滑動(dòng)床高度的影響微乎其微。

圖6 調(diào)節(jié)管道傾角的裝置Fig.6 Device for regulating inclination angle of pipeline

5 實(shí)驗(yàn)數(shù)據(jù)分析及模型驗(yàn)證

5.1 滑動(dòng)床上方速度分布模型的驗(yàn)證

速度分布的計(jì)算和實(shí)測(cè)值如圖7和8所示，試驗(yàn)工況為表1和2中第3組所列。

計(jì)算發(fā)現(xiàn)，上層速度分布理論和實(shí)測(cè)度誤差不超過(guò)10%，這說(shuō)明前文中垂向剪應(yīng)力呈線(xiàn)性分布的假設(shè)未導(dǎo)致過(guò)大誤差，假設(shè)成立。

圖7 水平管道速度分布計(jì)算值和實(shí)測(cè)值對(duì)比Fig.7 Comparison of calculated and measured values of velocity distribution in horizontal pipe

圖8 傾斜管道速度分布計(jì)算值和實(shí)測(cè)值對(duì)比Fig.8 Comparison of calculated and measured values of velocity distribution in inclined pipe

圖9 淤積臨界流速與顆粒粒徑和輸送濃度的關(guān)系Fig.9 Relationship between critical deposit velocity and particle transport concentration

圖10 淤積臨界流速和管道傾角關(guān)系Fig.10 The relationship between the critical deposit velocity and the inclination of the pipeline

計(jì)算中發(fā)現(xiàn)，隨著管道傾角的增大，零剪應(yīng)力（即最大速度點(diǎn)）有向上偏移的趨勢(shì)。這是由于在實(shí)驗(yàn)中，輸送濃度不變條件情況下，上下層速度差U1－U2的增量大于速度U1的增量，導(dǎo)致界面剪切應(yīng)力隨著管道傾角增大而增大，從圖3可知，必然導(dǎo)致y1增大，進(jìn)而零剪應(yīng)力（即最大速度點(diǎn)）向上偏移。

5.2 臨界淤積速度計(jì)算模型的驗(yàn)證

分析實(shí)驗(yàn)數(shù)據(jù)后發(fā)現(xiàn)，得到的淤積臨界流速Udl隨著顆粒平均粒徑的d增大而減小，如圖9所示。同時(shí)，圖9中試驗(yàn)數(shù)據(jù)還表明，輸送濃度CVd越大，則淤積臨界流速Udl的值越小。

圖10顯示了輸送濃度CVd不發(fā)生較大波動(dòng)（CVd＝0.21～0.22）條件下，管道傾角對(duì)淤積臨界流速Udl的影響。從圖中可以看出，隨著管道傾角的增大，3種級(jí)配顆粒的漿體淤積臨界流速Udl的值也會(huì)相應(yīng)增大。在管道角度超過(guò)25°以后，淤積臨界流速增加的趨勢(shì)有所減緩。

試驗(yàn)中還發(fā)現(xiàn)，其他條件相同時(shí)，粗顆粒滑動(dòng)床的速度比細(xì)顆?；瑒?dòng)床快。這種現(xiàn)象可以從理論上進(jìn)行解釋。由于作用在滑動(dòng)床床面的剪應(yīng)力τ0（顆粒床滑動(dòng)驅(qū)動(dòng)力）與上下層界面阻力系數(shù)f12正相關(guān)（見(jiàn)式（9）），而f12與滑動(dòng)床表面顆粒的粒徑正相關(guān)（見(jiàn)式（9））。由于粗顆粒床面的顆粒較大，導(dǎo)致界面阻力系數(shù)f12較大，從而粗顆粒床面的剪應(yīng)力τ0較細(xì)顆粒床面大，表現(xiàn)為粗顆?；瑒?dòng)床具有較快的滑動(dòng)速度。

在維持輸送濃度不變條件下，計(jì)算了臨界淤積狀態(tài)時(shí)滑動(dòng)床角度θ的值，發(fā)現(xiàn)對(duì)于不同的管道傾角，其變化幅度為20.12°～21.40°。且顆粒越粗，則該角度越小。管道傾角越大，該值越小。

當(dāng)輸送濃度CVd不變時(shí)，滑動(dòng)床角度隨著管道傾角增大而呈現(xiàn)輕微減小趨勢(shì)，如圖11所示。由此可以推斷，管道傾角對(duì)于滑動(dòng)床厚度的影響輕微。從而可以推斷，從水平管道開(kāi)始，隨著管道傾角的增大，為了維持輸送濃度不變，必須增大漿體速度Um，增大的Um主要用于平衡滑動(dòng)床越來(lái)越大的重力下滑力分量，在此過(guò)程中，也稍稍減小了滑動(dòng)床的厚度，如圖11所示。

圖12為計(jì)算的淤積臨界速度和實(shí)測(cè)的淤積臨界速度的對(duì)比，計(jì)算值和實(shí)測(cè)值最大偏差為13%。可見(jiàn)，本文提出的淤積臨界速度計(jì)算模型，用來(lái)預(yù)測(cè)顆粒粒徑0.4、0.6和1.71mm的砂漿體淤積臨界速度是較為可靠的。本文提出的方法對(duì)于其它條件下漿體淤積臨界速度的預(yù)測(cè)，模型是否適用，還需要進(jìn)一步檢驗(yàn)。

圖11 管道傾角φ變化對(duì)滑動(dòng)床角度θ的影響Fig.11 The influence of the change of the pipe angleφon angle of the sliding bedθ

圖12 淤積臨界流速Udl計(jì)算值和實(shí)測(cè)值對(duì)比Fig.12 Comparison of calculated value and measured value of critical deposit velocity Udl

6 結(jié) 論

（1）在水平管道擴(kuò)散方程的基礎(chǔ)上，提出了傾斜管道懸浮層顆粒濃度分布、滑動(dòng)床厚的計(jì)算模型。

（2）通過(guò)分析滑動(dòng)床受力平衡關(guān)系，給出了傾斜管道淤積臨界速度的計(jì)算模型。

（3）通過(guò)垂向剪應(yīng)力呈線(xiàn)性分布假設(shè)，給出了滑動(dòng)床上方懸浮層速度分布模型。

（4）通過(guò)分析顆粒粒徑為0.4、0.6和1.71mm，輸送濃度在0.10～0.27的3種粗砂漿體在直徑為150mm的試驗(yàn)管路中的流動(dòng)試驗(yàn)數(shù)據(jù)，發(fā)現(xiàn)上層速度分布模型理論計(jì)算值和實(shí)測(cè)度誤差不超過(guò)10%。淤積臨界流速與平均粒徑和輸送濃度負(fù)相關(guān)。當(dāng)輸送濃度不變時(shí)，管道傾角對(duì)于滑動(dòng)床厚度的影響甚微。淤積臨界速度的實(shí)驗(yàn)值模型計(jì)算值的對(duì)比表明，最大偏差不超過(guò)13%。

［1］張士林.漿體管道淤積流速數(shù)值計(jì)算方法研究［J］.水力采煤與管道運(yùn)輸，2011，（4）：6－9.Zhang S L.Research on numerical calculation method of slurry pipeline deposition velocity［J］.Hydraulic Coal Mining ＆Pipeline Transportation，2011，（4）：6－9.

［2］郗夏楠，孫西歡，李永業(yè)，等.管道水力輸送臨界流速研究進(jìn)展綜述［J］.山西水利，2012，（3）：36－38.Xi X N，Sun X H，Li Y Y，et al.Review of research progress on critical flow velocity of pipeline hydraulic transport［J］.Shanxi Water Resources，2012，（3）：36－38.

［3］陳琴瑞，李甲.漿體管道輸送臨界流速經(jīng)驗(yàn)公式適宜性分析［J］.中國(guó)礦山工程，2015，4（6）：73－75.Chen Q R，Li J.Suitability analysis of empirical formulas for critical velocity in slurry pipeline transportation［J］.China Mine Engineering，2015，4（6）：73－75.

［4］王新民，李天正，張欽禮.基于GA－ELM漿體管道輸送臨界流速預(yù)測(cè)模型研究［J］.中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2015，11（8）：101－105.Wang X M，Li T Z，Zhang Q L.Study on prediction model of critical flow velocity in slurry pipeline transportation based on GA－ELM［J］.Journal of Safety Science and Technology，2015，11（8）：101－105.

［5］秦德慶，曹斌，夏建新.不同顆粒物料管道水力輸送不淤臨界流速的確定［J］.礦冶工程，2014，34（1）：9－11.Qin D Q，Cao B，Xia J X.Study on non－depositing velocity of different particle materials in pipeline by hydraulic transportation［J］.Mining and Metallurgical Engineering，2014，34（1）：9－11.

［6］曠水泉.大石河鐵尾礦高濃度管道輸送臨界流速研究［J］.水力采煤與管道運(yùn)輸，2015，（2）：1－5.Kuang S Q.Study on critical flow velocity of high concentration pipeline transportation of iron tailings in Dashihe［J］.Hydraulic Coal Mining ＆Pipeline Transportation，2015，（2）：1－5.

［7］Souza P T C，Moraes J D，Slatter P T，et al.Modelling the critical velocity for heterogeneous flow of mineral slurries［J］.In International Journal of Multiphase Flow，2014，65：31－37.

［8］Sanni E S，Olawale S A，Adefila S S.Theoretical study of sand entrainment and deposits in horizontal oil transport［J］.Science Research，2015，3（6）：314－323.

［9］Rice H P，F(xiàn)airweather M，Peakall J，et al.Constraints on the functional form of the critical deposition velocity in solid－liquid pipe flow at low solid volume fractions［J］.Chemical Engineering，2015，（126）：759－770.

［10］Boulanger J，Joan B.Sand suspension deposition in horizontal low－concentration slurry pipe flows［J］.Granular Matter，2016，18（2）：1－10.

［11］Harbottle D，F(xiàn)airweather M，Biggs S.The minimum transport velocity of colloidal silica suspensions［J］.Chemical Engineering Science，2011，66（11）：2309－2316.

［12］Rice H.Transport and deposition behaviour of model slurries in closed pipe flow［D］.Leeds：University of Leeds，2013：112－122.

［13］KotzéM，Sutherland A，KotzéR.A system to estimate coarse particle velocities at the pipe wall in settling slurry flow［J］.Flow Measurement and Instrumentation，2013，（32）：63－70.

［14］Mario R R，Eduardo S A.Two－layer model for horizontal pipe flow of newtonian and non－newtonian settling dense slurries［J］.Ind Eng Chem Res，2012，51（20）：7095－7103.

［15］Abd Al Aziz A I，Mohamed H I.A study of the factors affecting transporting solid－liquid suspension through pipelines［J］.Open Journal of Fluid Dynamics，2013，（3）：152－162.

［16］Vaclav M，Jan K.Stationary－and sliding beds in pipe flows of settling slurry［C］／／The Proceedings of the 15th International Freight Pipeline Society Symposium.Prague：Taylor ＆Francis，2014：231－237.

［17］Matouseka V，Krupickaa J，Penika V.Distribution of mediumto－coarse glass beads in slurry pipe flow：evaluation of measured concentration profiles［J］.Particulate Science and Technology，2014，32（2）：186－196.

［18］Matouseka V，Krupickaa J，Penika V.Distribution of mediumto－coarse glass beads in slurry pipe flow：evaluation of measured concentration profiles［J］.Particulate Science and Technology，2014，32（2）：186－196.

［19］Sape A M.An overview of theories describing head losses in slurry transport：a tribute to some of the early researchers［J］.Pipeline and Riser Technology，2013，4（03）：233－241.

［20］趙利安.大顆粒漿體管內(nèi)流動(dòng)規(guī)律研究［D］.阜新：遼寧工程技術(shù)大學(xué)，2011：98－109.Zhao Li－an.Study on flow law of large particle slurry in pipeline［D］.Fuxin：Liaoning Technical University，2011：98－109.

［21］林雪松，陳殿強(qiáng)，何峰，等.固液兩相流中顆粒濃度分布統(tǒng)一公式的研究［J］.水資源與水工程學(xué)報(bào)，2013，24（4）：82－84.Lin X S，Chen D Q，He F，et al.Study on unified formula of particle concentration distribution in two－phase flow of liquidsolid［J］.Journal of Water Resources ＆Water Engineering，2013，24（4）：82－84.

［22］傅旭東，王光謙.低濃度固液兩相流中的粗顆粒濃度分布［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，42（10）：1361－1364.Fu X D，Wang G Q.Concentration distribution of coarse particles in dilute solid／liquid two－phase floes［J］.Journal of Tsinghua University（Science ＆Technology），2002，42（10）：1361－1364.

［23］趙利安，孟慶華.伴有滑動(dòng)床的漿體傾斜管道摩阻損失研究［J］.礦業(yè)研究與開(kāi)發(fā)，2008，28（2）：3－6.Zhao L A，Meng Q H.Research on friction loss of slurry with sliding bed in inclined pipeline［J］.Mining Research and Development，2008，28（2）：3－6.

［24］Shook C A，Roco M C.Slurry flow：principles and practice［M］.Oxford：Butterworth－Heinemann，1991：123－124.

［25］張士林，許振良，邵龍?zhí)?弓形斷面管道的紊流速度分布［J］.遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2006，25（1）：63－65.Zhang S L，Xu Z L，Shao L T.Water velocity in pipe of arching section［J］.Journal of Liaoning Technical University，2006，25（1）：63－65.

［26］Hisham A N，Shook C A，Nabil E.Isokinetic probe sampling from slurry pipelines［J］.The Canadian Journal of Chemical Engineering，1984，62（2）：179－185.

［27］趙利安，姜威.粗顆粒高濃度流體管道流動(dòng)濃度分布研究［J］.泥沙研究，2016，（1）：37－41.Zhao L A，Jiang W.Study on concentration distribution of pipe flow with high concentration coarse particles［J］.Journal of Sediment Research，2016，（1）：37－41.

Study on the flow of water sand slurry with sliding bed in inclined pipe

Zhao Li’an1，＊，Wang Tieli2
（College of Mining，Liaoning Technical University，F(xiàn)uxin Liaoning 123000，China；Department of Pipeline Coal Transporting，Wuhan Design and Research Institute of China Coal Technology and Engineering Group，Wuhan 430064，China）

In order to further study sand particles’transport rules in the inclined slurry pipe associated with the sliding bed，force balance theories are used to study the method for evaluating the suspended layer and sliding bed velocities as well as the sliding bed thickness，and an iterative method for computing the deposit critical velocity is given.Meanwhile，the shear force linear distribution assumption is used and with the aid of prandt empirical formula，a velocity distribution model above the sliding bed is given.Experimental data of coarse sand pipe hydraulic transportation show that the deposit critical velocity is negatively correlated with the particle size and concentration，and it is positively correlated with the pipe dip angle when the transportation concentration is kept constant.The comparison between the calculated results and the measured results show that the maximum deviation of the calculated value from the measured value is 13%.The experimental data also shows that with the increase of the angle of the pipe，the maximum velocity point position has an upward trend.The comparison between the calculated and measured results shows that the deviation of the theoretical velocity distribution model and the measured value is not more than 10%.Results show that the force balance theory and the linear hypothesis of the shear stress distribution apply well in the research of deposition critical velocity and velocity distribution above sliding bed in the inclined coarse sand slurry pipe accompanied with the sliding bed.

settling slurry；deposit critical velocity；sliding bed；diffusion equation；velocity distribution

TV142

A

（編輯：楊 娟）

1672－9897（2016）06－0043－07

10.11729／syltlx20160073

2016－04－27；

2016－08－18

中國(guó)煤炭科工集團(tuán)科技創(chuàng)新基金項(xiàng)目（2014QN040）

＊通信作者E－mail：anlizhao1＠163.com

Zhao L A，Wang T L.Study on the flow of water sand slurry with sliding bed in inclined pipe.Journal of Experiments in Fluid Mechanics，2016，30（6）：43－49.趙利安，王鐵力.伴有滑動(dòng)床的水砂漿體傾斜管道流動(dòng)研究.實(shí)驗(yàn)流體力學(xué)，2016，30（6）：43－49.

趙利安（1973－），男，陜西長(zhǎng)安人，博士，講師。研究方向：漿體輸送理論研究。通信地址：遼寧省阜新市中華路47＃遼寧工程技術(shù)大學(xué)礦業(yè)學(xué)院采礦系（123000）。E－mail：an lizhao1＠163.com