水利工程論文-洛惠渠高含沙水流的特性分析.doc

水利工程論文-洛惠渠高含沙水流的特性分析摘要：利用洛惠渠近二十年的實測資料，從泥沙顆粒沉速、流變特性等基本概念著手分析洛惠渠高含沙水流的特性。指出水流中含沙量的增大和細(xì)顆粒含量的增加，一方面使流體粘度增加；另一方面使流體容重增大，從而導(dǎo)致顆粒的沉速大幅度降低，甚至形成不分選泥漿。文章為進(jìn)一步研究黃河等多沙河流的高含沙水流輸沙提供了科學(xué)依據(jù)。關(guān)鍵詞：洛惠渠高含沙沉速粘性賓漢體1灌區(qū)概況人民引洛灌區(qū)位于陜西關(guān)中東部，是河源泥沙干擾比較突出、泥沙含量較大的灌區(qū)。灌區(qū)平面圖見圖1，總干渠全長21.5km，沿洛河?xùn)|岸溝壑區(qū)經(jīng)五號隧洞到達(dá)義井，再分為東、中、西三條干渠，呈扇形分布。圖1人民引洛渠灌區(qū)平面圖SketchofRenminYinLuoIrrigationDistrict灌區(qū)水源為北洛河，它是黃河多沙支流之一，發(fā)源于陜北定邊縣白宇山區(qū)。北洛河多年平均平均流量為25.7m3/s；洪水期79月，最大洪峰流量為4420m3/s；枯水期56月，最小流量為1.2m3/s；一般流量在10m3/s左右。氵狀頭站多年平均水、沙及各級含沙量(以重量百分比計)情況見表1。表1北洛河氵狀頭站多年平均水、沙情況統(tǒng)計表(19521974)AnnualrunoffandsedimentloadatZhuangtoustation年徑流量年沙量最大含沙量(%)各級含沙量出現(xiàn)歷時(小時)總量(億m3)79月(%)總量(萬t)79月(%)1525%2535%35%小計8.144595009065.0120.470.6159.3350.3據(jù)前十五年的資料統(tǒng)計，一次洪峰持續(xù)高含沙的情況如下：含沙量大于40%、歷時超過四天的有三年，超過三天的有九年；洪峰含沙量大于50%的幾乎每年都出現(xiàn)，其中持續(xù)兩天以上的有兩年。北洛河洪水一般多出現(xiàn)于78月份，最早于6月中、下旬出現(xiàn)。洪峰和沙峰基本相適應(yīng)，沙峰稍遲于洪峰，洪峰峰型陡瘦，沙峰陡漲緩落。人民引洛渠直接引用北洛河水，水沙關(guān)系與北洛河水變化過程基本相同。從1974年至1986年實測資料分析，引洛渠汛期含沙量與懸沙粒徑有一定關(guān)系，即：含沙量越大，粒徑越粗(見圖2)；從汛期各次洪峰相比較，由于來水來沙區(qū)不同，懸沙粒徑也有所不同。一般來說，每年第一次洪峰的懸沙粒徑較粗。據(jù)19631970年資料統(tǒng)計；北洛河懸沙平均粒徑為0.03990.0437mm，中值粒徑為0.03530.0420mm。2泥沙顆粒在高含沙水流中的特性顆粒沉速是泥沙運動中的一個重要參數(shù)，如沉淀落淤速度、挾沙能力等都與這一參數(shù)密切相關(guān)。圖2東干義井?dāng)嗝婧榉鍟r段Sd50關(guān)系圖Relationbetweensiltconcentrationandd50atYijingsectioninpeaktime已有的研究表明含沙量對泥沙顆粒的沉速有著一定的影響。粘土漿液的含沙量達(dá)到一定程度以后，具有非牛頓流體的特性。國內(nèi)外學(xué)者多采用賓漢方程來近似描述其流變特性1，2。2.1洛惠渠高含沙水流的流變特性1985年費祥刻授3以愛因斯坦在球形顆粒的稀懸液條件下，推得的相對粘度理論公式液/0=(1+2.5Sv)為基礎(chǔ)，分析求得不同粒度組成下高濃度懸液的相對粘度表達(dá)式，并做了比較系統(tǒng)的粘性試驗進(jìn)行驗證，最后得出高濃度泥沙懸液相對粘度的公式液/0=1-KSv/Svm-2.5(1)式中系數(shù)K是對固體濃度的修正系數(shù)，即考慮了細(xì)顆粒表面的薄膜水，又考慮了粒間封閉水對濃度的影響。Sv和Svm分別表示固體體積比濃度及其極限濃度，液和0則表示泥沙懸液及同溫度的清水的粘滯系數(shù)。當(dāng)懸液中固體濃度較高時，已屬非牛頓流體，可用賓漢模型來描述其剪切力與切變率關(guān)系。賓漢極限剪切力的試驗結(jié)果可用下列經(jīng)驗式表達(dá)4B=9.810-2exp(B+1.5)(2)式中BN/m2,系數(shù)B=8.45，=Sv-Sv0/Svm這里Sv0為由牛頓體轉(zhuǎn)變到非牛頓體的臨界濃度。公式(1)和(2)的特點是除了濃度Sv作為影響懸液流變參數(shù)的因素以外，又加入了反映顆粒組成特性的Svm值，使流變參數(shù)的應(yīng)用范圍明顯擴(kuò)大。根據(jù)采集黃河中下游有關(guān)水文站的懸沙沙樣進(jìn)行的流變試驗3，得出如下關(guān)系式(3)Svm=0.92-0.2lgpi/di(4)K=1+2.0（Sv/Svm）0.3（1-Sv/Svm）4(5)洛惠渠引水是引黃河之水，上述公式可以應(yīng)用于此。2.2泥沙顆粒在賓漢流體中的受力分析下面，我們分析一下作用于流體中下沉顆粒上的力。當(dāng)泥沙顆粒的水下有效重量與作用于顆粒的阻力平衡時，它便達(dá)到其最終沉速(s-)d3/6=CDd2/42/2g(6)式中s為泥沙顆粒比重，s=sg；為流體比重。知道了阻力系數(shù)如何變化，就可以推得沉速如何變化。作用于泥沙顆粒的阻力等于作用于泥沙顆粒表面的切應(yīng)力和壓應(yīng)力沿運動方向的合力。泥沙顆粒在牛頓流體中作層流沉降時，作用于其上的總阻力(即式(6)等號右邊的項)為3d，其中三分之一是壓應(yīng)力的合力，三分之二是切應(yīng)力的合力。若CDD2/42/2g=3d(7)CD=24/Re(8)Re=d/v(9)泥沙顆粒在做紊流沉降時，切應(yīng)力的作用便退居次席。在充分紊動區(qū)，阻力系數(shù)C趨于一個常數(shù)；在層流沉降與紊流沉降之間，則有一個過渡區(qū)?？偟膩碚f，阻力系數(shù)CD是雷諾數(shù)的函數(shù)。泥沙顆粒在賓漢流體中沉降時，不僅有相當(dāng)于牛頓流體動力粘滯系數(shù)的剛度系數(shù)，賓漢切應(yīng)力B也會引起水流阻礙泥沙顆粒的運動。參照安斯列的概念，下面討論作用于賓漢體中層流沉降顆粒上的力5。可以假定由剛度系數(shù)(相當(dāng)于)引起的力仍為3d，由B引起的力由兩部分組成：1.設(shè)作用于球體表面的切應(yīng)力為B，則它們沿垂直方向的合力為1/42d26；2.根據(jù)塑性流滑動線理論，作用于球面的壓應(yīng)力沿垂直方向的合力為5/82d2B。所有這些力的合力亦即總阻力，應(yīng)為3d+7/82d2B。如果同樣用一般阻力表達(dá)形式，可得CDd2/42/2g=3d+7/8X)2d2B(10)若認(rèn)為在層流區(qū)CD仍與雷諾數(shù)Re1成反比CD=24/Re1(11)則由上式可以得出修正雷諾數(shù)Re1的形式為Re1=2/d+7/24B(12)對于牛頓流體，B=0這時Re1即轉(zhuǎn)化成式(9)形成的Re。對比泥沙顆粒在牛頓體中的沉降情況，泥沙顆粒在賓漢流體中沉降時，阻力系數(shù)CD也應(yīng)是修正后的雷諾數(shù)Re1的函數(shù)。將上述方程式(6)和(10)聯(lián)和求解(s-)d3/6=3d+7/82d2B上式為清水比重。如果泥沙顆粒在高含沙水流中運動，那么就用渾水比重液來代替清水比重，即(s-液)d3/6=3d液+7/82d2B液=(s-液)d2/6-7/8dB/3(13)式中液是泥沙顆粒在渾水中的沉速。根據(jù)流體的流變特性可知，當(dāng)含沙濃度SvSv0時，液體是牛頓體，B=0,剛度系數(shù)(=)，式(13)可化為液=(s-液)d2/18液(14)式中液根據(jù)方程式(1)求得，含沙量越大，液越大，且大于相同溫度下清水的粘滯系數(shù)0。通過上式可以看出泥沙顆粒在高清不同，根據(jù)洛惠渠資料可以點出關(guān)系曲線，如圖3所示。從圖3可以看出，相同的粒徑組成在高含沙水流中的平均沉速遠(yuǎn)遠(yuǎn)小于在清水中的平均沉速，有的只相當(dāng)于清水中沉速的1/7。究其原因：一方面使泥沙顆粒受的上浮力增大；另一方面是顆粒在高含沙水流中的沉速減小。由于泥沙顆粒在高含沙水流中的沉速減小，圖3泥沙顆粒在清水和渾水中的沉速Fallvelocityinclearwaterandsiltywater如果采用Sv3/gh關(guān)系作為挾沙能力的基礎(chǔ)S=k(v3/gh)m,那么在流速v和水深h一定的情況下，挾沙能力就增大，從而使得泥沙淤積量減小。這更有利于泥沙的遠(yuǎn)距離輸送。當(dāng)SvSv0時，高含沙水流成為賓漢體，切應(yīng)力B、液和剛度(相當(dāng)于粘性)三種因素共同影響造成泥沙顆粒在賓漢體中的沉速與在清水中的沉速不一樣。而當(dāng)高含沙水流一旦成*.齪禾澹營瑽的切應(yīng)力便可以托住一般細(xì)小顆粒，使其不下沉。如果令式(13)等于零，如下式(s-液)d2max/6-7/8dmaxB=0dmax=21B/4(s-液)便可以得到相應(yīng)含沙量情況下能使泥沙顆粒不下沉(即=0)的最大粒徑dmax。根據(jù)洛惠渠的資料，計算得到dmax如表2。從表2可以看出dmax是很大的。從另一個角度，我們也可以認(rèn)為當(dāng)高含沙水流形成賓漢體后，挾沙能力大大增加。這也是高含沙水流產(chǎn)生揭河底現(xiàn)象的一個重要原因。表2高含沙(賓漢體)水流不沉的最大粒徑值Themaximumsizeofnosettinginhyperconcentratedflow測次時間D平均(mm)Sv液(kg/cm3)B(N/m2)dmax(mm)S71974.7.29.7:450.04590.3290571.5429438.64328412.0452S87.29.18:480.0440.3362261.55477411.2790515.87396S97.30.7:220.02520.3218871.53111320.8247128.73945S107.30.17:300.03220.3188681.52613233.5617846.12892S117.31.19:450.03440.2347171.3872834.7864535.908415S138.1.18:100.0330.2177361.3592645.0369056.09231S31977.6.27.8:560.04640.3037741.5012265.1795816.977082S71978.7.28.17:000.03270.2237741.3692264.734135.767407S11980.6.18.11:000.04230.307171.506839.15633612.38951S51974.7.29.7:300.04990.3362261.5547749.01856512.692593結(jié)語通過對洛惠渠水沙資料的分析，可以看出，含沙量越大，泥沙粒徑越粗，流體容重就越大，粘性系數(shù)也增大。從而導(dǎo)致沉速大大降低，流體挾沙能力增強(qiáng)，河道淤積量減少。參考文獻(xiàn)1錢寧，萬兆惠。高含沙水流研究述評。水利學(xué)報，1985，(5).2ByronBird,R,DaiG.C.andBarbaraJ.Yarusso,TheRheologyandFlowofViscoplasticMaterials.ReviewsinChemicalEngineering,Vol.I,No.1,1982.3費祥俊。黃河中下游含沙水流粘度的計算模型。泥沙研究，1991，(2).4FeiXiangjun,YangMeiqing,ThePhysic