水利工程論文-黃河包頭河段泥沙模型設(shè)計(jì).doc

水利工程論文-黃河包頭河段泥沙模型設(shè)計(jì)摘要：針對(duì)黃河泥沙多且變化大的特點(diǎn),采用泥沙運(yùn)動(dòng)相似條件的最新研究成果,完成了黃河包頭河段動(dòng)床模型沙選擇和比尺設(shè)計(jì)。利用原型實(shí)測(cè)資料進(jìn)行的驗(yàn)證試驗(yàn),表明模型小河較好地復(fù)演了原型河勢(shì)、水沙運(yùn)動(dòng)規(guī)律及河床沖淤變形等。關(guān)鍵詞：泥沙模型模型設(shè)計(jì)黃河1模擬河段河道概況黃河包頭段位于著名的河套地區(qū)。包頭畫匠營(yíng)子河段位于包頭市新市區(qū)略偏東處，距新市區(qū)約12km，在昭君墳和磴口之間（昭君墳水文站與磴口水文站相距約53km）。試驗(yàn)河段包括包頭黃河公路橋下游1800m，兩橋間800m，鐵路橋上游3200m，共計(jì)約6000m。對(duì)昭君墳水文站1934年以來(lái)的水文資料整理分析后，按中水3000m3/s的流量考慮，得到原型河段代表性資料如表1。表1原型河道水流要素特征值Characteristicvaluesofriverflowinprototype項(xiàng)目取值項(xiàng)目取值河段長(zhǎng)度L6000m河段比降i0.10.2河道寬度B1000m懸移質(zhì)中徑d500.028mm主河道平均寬度b400m河床床沙中徑D500.18mm平均水深H4.5m最大流量Q1%6200m3/s平均流速V1.7m/s最小流量Qmin48m3/s河床糙率nb0.0150.020泥沙比重s2.652.70河灘糙率nt0.0200.030泥沙沉速0.22cm/s2河流泥沙模型相似條件目前黃河動(dòng)床模型相似律的研究較多14，并結(jié)合黃河包頭河段的特性，采用下列相似比尺關(guān)系作為本模型設(shè)計(jì)的依據(jù)。重力相似V=h1/2(1)阻力相似n=H2/3/L1/2(2)由水流連續(xù)方程,可分別得出流量及時(shí)間比尺關(guān)系式Q=VLH(3)t1=L/V(4)懸移質(zhì)相似=V(H/L)3/4(5)對(duì)于細(xì)沙河流,當(dāng)懸沙粒徑d0.15mm時(shí)(原形懸沙均能滿足這一些條件),懸沙沉速用滯留區(qū)公式計(jì)算,由此得懸沙粒徑比尺d=(V/s-)1/2(6)底沙起動(dòng)及揚(yáng)動(dòng)相似V=Vc=Vf(7)水流輸沙率相似S=S*(8)河床沖淤變形相似t2=o2LH/Gs(9)對(duì)于懸沙模型,上式又可表示為t2=o/st1(10)河型相似(s-/D50H)1/3/iB2/3m(s-/D50H)1/3/iB2/3p(11)3幾何比尺及模型沙選取3.1幾何比尺及模型變率根據(jù)試驗(yàn)河段平面范圍及試驗(yàn)場(chǎng)地條件，選取水平比尺L=150?？紤]到模型水流深度應(yīng)滿足表面張力及試驗(yàn)量測(cè)的要求，對(duì)模型幾何變率前期的研究成果和模型沙特性等方面權(quán)衡后，取垂直比尺H=50，則幾何變率Dt=3。采用張紅武3提出的變態(tài)模型相對(duì)保證率的公式求得Dt5.6；利用張瑞瑾4提出的模型變態(tài)限制指標(biāo)求得Dt6.9。驗(yàn)算結(jié)果表明本模型采用的變率在各家公式限制的范圍之內(nèi),幾何變態(tài)的影響有限，可以滿足工程需要。3.2模型沙選取根據(jù)式(2)求得n=1.11(取R=H),因而模型糙率nm=0.0140.018,經(jīng)比選后，決定采用經(jīng)過(guò)專門加工處理的蘭州煤廠精煤粉作為試驗(yàn)?zāi)Ｐ蜕场Ｔ摬牧喜粌H比重小，而且密實(shí)度及凝聚力也較小(試驗(yàn)測(cè)得凝聚力C=0.018kg/m2,內(nèi)摩擦角=35，容重s=1.401.42kg/m3，干容重00.620.71kg/m3，含水量為60%，D50=0.030.05mm)，因而能保證模型河床的活動(dòng)性，以適應(yīng)原型河床沖淤變化幅度較大的特點(diǎn)。原型沙容重sp=2.662.70kg/m3,干容重op=1.401.50kg/m3。模型設(shè)計(jì)時(shí),取模型沙sm=1.41kg/m3,om=0.67kg/m3；取原型沙sp=2.68kg/m3,op=1.45kg/m3。則s-=4.09o=2.16由式(5)得：=3.10；由式(6)得,d=0.78(根據(jù)對(duì)原型及模型水流的溫=0.8)。要求模型懸沙中徑d50=0.028/0.78=0.036mm。對(duì)模型底沙，為保證模型小河的綜合穩(wěn)定性與原型相近，需按照河型相似條件進(jìn)行選擇，即將原型及模型的有關(guān)數(shù)據(jù)代入式(11)，得模型床沙中徑D50=0.042mm。所對(duì)應(yīng)的床沙粒徑比尺D=0.18/0.042=4.29。為檢查底沙的糙率，選用D50=0.040.06mm的煤粉作為模型沙進(jìn)行預(yù)備試驗(yàn)，其動(dòng)床河床糙率為0.0150.018，下限0.015為靜平床糙率，上限0.018為河床出現(xiàn)沙紋后的動(dòng)床糙率，故可滿足阻力相似條件。為檢查所選模型沙能否滿足底沙起動(dòng)及揚(yáng)動(dòng)相似條件，需分別確定模型沙及原型沙的起動(dòng)流速和揚(yáng)動(dòng)流速。為此，采用煤粉進(jìn)行了水槽試驗(yàn)，所得起動(dòng)及揚(yáng)動(dòng)流速試驗(yàn)結(jié)果見表2。表2模型沙起動(dòng)流速試驗(yàn)結(jié)果Theexperimentresultofthresholdvelocityofmodelsediment速度V(m/s)水深備注(cm)弱動(dòng)普動(dòng)揚(yáng)動(dòng)24.506.009.58sm=1.401.4245.507.5012.4066.508.5013.78d50=0.040.06mm87.209.2014.30為確定原型底沙的起動(dòng)流速，點(diǎn)繪了黃河磴口、昭君墳、巴彥高勒及三湖河口水位站資料和引黃渠道沙土不沖流速與含沙量野外資料的關(guān)系(見圖1)，以分析不沖流速與含沙量的關(guān)系，曲線與縱軸的交點(diǎn)，即為清水時(shí)的不沖流速。一般認(rèn)為清水時(shí)的不沖流速等于起動(dòng)流6，因而得到原型沙的起動(dòng)流速Vcp=0.68m/s，相應(yīng)水深范圍內(nèi)的模型沙起動(dòng)流速Vcm=0.075-0.092m/s，可以得到Vc0.680.075-0.092=7.39-9.07下面采用沙玉清起動(dòng)流速公式7計(jì)算天然河流床沙起動(dòng)流速，該公式形式為圖1不沖流速與含沙量關(guān)系RelationbetweennonscouringvelocityandsiltcontentVc=0.43d3/4+1.1(0.7-)4/d1/2h1/5(12a)式中為孔隙率，其穩(wěn)定值約為0.4，d以mm計(jì)，Vc以m/s計(jì)。將d=0.18mm代入式（12a）得Vc=0.41h1/5(12b)當(dāng)水深h=1-3.5m時(shí)，采用式(12b)及表2結(jié)果可求得相應(yīng)的vc值列于表3。表3原型沙起動(dòng)流速Vcp及起動(dòng)流速比尺vc值ThevaluesofthresholdvelocityVcpandthevelocityscalevcofprototypesediment水深(m)1.01.52.03.03.5Vcp(m/s)0.410.450.470.510.53vc6.836.626.276.005.96備注v=7.7由表3所確定的起動(dòng)流速比尺略小于流速比尺。根據(jù)以上情況綜合分析可知，所選模型沙(包括底沙及懸沙中那部分與底沙經(jīng)常交換的床沙質(zhì)泥沙)基本能滿足起動(dòng)相似條件。黃河天然沙揚(yáng)動(dòng)流速一般為起動(dòng)流速的1.541.75倍3，取原型底沙Vfp=1.65Vcp，則得到原型沙相應(yīng)水深下的揚(yáng)動(dòng)流速為0.84m/s。由表1查得模型沙揚(yáng)動(dòng)流速Vfm=0.124-0.143m/s。從而Vf=6.77-5.87。所得揚(yáng)動(dòng)流速比尺與流速比尺較接近，表明所選模型沙還滿足揚(yáng)動(dòng)相似條件(即vf=v)。從上述模型選沙結(jié)果來(lái)看，模型床沙中經(jīng)常被水流沖起變?yōu)閼疑车哪遣糠诸w粒，其粒徑接近模型沙的粒徑，因而自然滿足懸沙相似條件。3.3含沙量比尺及時(shí)間比尺式(9)是由河床變形方程式推導(dǎo)而來(lái)的，考慮了全沙變化對(duì)于河床變形的影響。由于原型實(shí)際上長(zhǎng)期缺乏推移質(zhì)泥沙觀測(cè)資料，在試驗(yàn)中對(duì)推移質(zhì)輸沙進(jìn)行定量模擬是相當(dāng)困難的。不過(guò)，根據(jù)對(duì)原型情況分析后認(rèn)為，由懸移質(zhì)泥沙運(yùn)動(dòng)所引起的沖淤變化構(gòu)成了該河段河床變形的主要部分，且推移質(zhì)泥沙時(shí)常與懸移質(zhì)泥沙交換，進(jìn)入懸移質(zhì)運(yùn)動(dòng)的行列中，因此本模型僅考慮懸移質(zhì)輸沙量變化的影響，并按照式(10)確定河床的沖淤變形時(shí)間比尺。事實(shí)上，通過(guò)河床驗(yàn)證試驗(yàn)確定模型進(jìn)口加沙量時(shí)，自然反映了推移質(zhì)運(yùn)動(dòng)對(duì)河床沖淤變形影響的相似性。另外，一般試驗(yàn)只模擬懸移質(zhì)中的床沙質(zhì)，但對(duì)于原型情況，懸移質(zhì)中的沖瀉質(zhì)在河灘造床過(guò)程中起到了很大的作用。在試驗(yàn)中若相應(yīng)扣除來(lái)沙量中的沖瀉質(zhì)，不僅使灘地淤積難以相似，而且還會(huì)使主槽變形產(chǎn)生偏離。因此本試驗(yàn)不再對(duì)懸移質(zhì)中的床沙質(zhì)和沖瀉質(zhì)加以區(qū)劃。考慮到現(xiàn)有挾沙力公式不能同時(shí)應(yīng)用于原型及模型，復(fù)雜河流的模型試驗(yàn)已不再直接利用由一些挾沙力公式直接導(dǎo)出的比尺關(guān)系式計(jì)算S*3。本試驗(yàn)分別確定原型和模型值，兩者之比為含沙量比尺S*。經(jīng)過(guò)試驗(yàn)河段附近實(shí)測(cè)資料驗(yàn)證后(見圖2)，采用與原型資料頗為符合的黃科院張紅武公式3，即S*=0.14(V3/ghlnh/6D50)0.6(13)圖2水流挾沙力公式與原型資料對(duì)比Comparisonbetweenformulaandmeasureddata圖3模型沙S*與V3/gh的關(guān)系RelationbetweenS*andV3/gh將水力泥沙因子代入上式，即V=1.7m/s,h=4.5m,R=4.5m,=0.22cm/s,D50=0.00018m,得S*p=5.26kg/m3利用模型沙進(jìn)行挾沙力試驗(yàn)后，得到的結(jié)果如圖3所示，點(diǎn)群遵循的關(guān)系可表示為S*m=0.18V3/gh(14)根據(jù)原型水力泥沙因子及相應(yīng)的比尺，求得模型的Vm=24.05cm/s,Hm=9cm,=0.071cm/s，將其代入(14)式，可得S*m=4.0kg/m3。由于模型的挾沙能力小于水槽中的挾沙能力，參照以往黃河模型試驗(yàn)的經(jīng)驗(yàn)，乘以折減系數(shù)K13(一般取K1=0.7-0.8）得s*=(0.7-0.8)4.0=2.8-3.2kg/m3則s*=5.26/2.8-3.2=1.64-1.88進(jìn)一步改變水力及泥沙因子計(jì)算后可得到s=s*=1.5-2.0。將s及o代入(10)，得河床變形時(shí)間比尺t2=(1.08-1.44)t1=22.9-30.6。表明兩個(gè)時(shí)間比尺相差較小，這對(duì)于非恒定動(dòng)床模型試驗(yàn)很有益處，不僅能保證試驗(yàn)的可靠性，也便于試驗(yàn)操作。4模型驗(yàn)證對(duì)原型較為系統(tǒng)的兩組同步水面線進(jìn)行了施放,模型實(shí)測(cè)點(diǎn)與原型的偏差一般不超過(guò)2mm(相當(dāng)于原型10cm),最大不超過(guò)3.2mm(相當(dāng)于原型16cm)。在試驗(yàn)河段公路橋處有長(zhǎng)期水位觀測(cè)資料，將模型中相應(yīng)的水位流量關(guān)系曲線與原型進(jìn)行比較，可知兩者頗為接近。對(duì)原型河段實(shí)測(cè)的5個(gè)流速斷面進(jìn)行驗(yàn)證，可以看到模型中流速沿河寬的分布與原型符合較好，表明滿足了水流相似條件。按原型資料，對(duì)模型水沙過(guò)程進(jìn)行概化。參照上述模型設(shè)計(jì)結(jié)果，分別選用不同的含沙量比尺及其相應(yīng)的河床變形