水利工程论文-黄河包头河段泥沙模型设计.doc

水利工程论文-黄河包头河段泥沙模型设计摘要：针对黄河泥沙多且变化大的特点,采用泥沙运动相似条件的最新研究成果,完成了黄河包头河段动床模型沙选择和比尺设计。利用原型实测资料进行的验证试验,表明模型小河较好地复演了原型河势、水沙运动规律及河床冲淤变形等。关键词：泥沙模型模型设计黄河1模拟河段河道概况黄河包头段位于著名的河套地区。包头画匠营子河段位于包头市新市区略偏东处，距新市区约12km，在昭君坟和磴口之间（昭君坟水文站与磴口水文站相距约53km）。试验河段包括包头黄河公路桥下游1800m，两桥间800m，铁路桥上游3200m，共计约6000m。对昭君坟水文站1934年以来的水文资料整理分析后，按中水3000m3/s的流量考虑，得到原型河段代表性资料如表1。表1原型河道水流要素特征值Characteristicvaluesofriverflowinprototype项目取值项目取值河段长度L6000m河段比降i0.10.2河道宽度B1000m悬移质中径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河流泥沙模型相似条件目前黄河动床模型相似律的研究较多14，并结合黄河包头河段的特性，采用下列相似比尺关系作为本模型设计的依据。重力相似V=h1/2(1)阻力相似n=H2/3/L1/2(2)由水流连续方程,可分别得出流量及时间比尺关系式Q=VLH(3)t1=L/V(4)悬移质相似=V(H/L)3/4(5)对于细沙河流,当悬沙粒径d0.15mm时(原形悬沙均能满足这一些条件),悬沙沉速用滞留区公式计算,由此得悬沙粒径比尺d=(V/s-)1/2(6)底沙起动及扬动相似V=Vc=Vf(7)水流输沙率相似S=S*(8)河床冲淤变形相似t2=o2LH/Gs(9)对于悬沙模型,上式又可表示为t2=o/st1(10)河型相似(s-/D50H)1/3/iB2/3m(s-/D50H)1/3/iB2/3p(11)3几何比尺及模型沙选取3.1几何比尺及模型变率根据试验河段平面范围及试验场地条件，选取水平比尺L=150。考虑到模型水流深度应满足表面张力及试验量测的要求，对模型几何变率前期的研究成果和模型沙特性等方面权衡后，取垂直比尺H=50，则几何变率Dt=3。采用张红武3提出的变态模型相对保证率的公式求得Dt5.6；利用张瑞瑾4提出的模型变态限制指标求得Dt6.9。验算结果表明本模型采用的变率在各家公式限制的范围之内,几何变态的影响有限，可以满足工程需要。3.2模型沙选取根据式(2)求得n=1.11(取R=H),因而模型糙率nm=0.0140.018,经比选后，决定采用经过专门加工处理的兰州煤厂精煤粉作为试验模型沙。该材料不仅比重小，而且密实度及凝聚力也较小(试验测得凝聚力C=0.018kg/m2,内摩擦角=35，容重s=1.401.42kg/m3，干容重00.620.71kg/m3，含水量为60%，D50=0.030.05mm)，因而能保证模型河床的活动性，以适应原型河床冲淤变化幅度较大的特点。原型沙容重sp=2.662.70kg/m3,干容重op=1.401.50kg/m3。模型设计时,取模型沙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(根据对原型及模型水流的温=0.8)。要求模型悬沙中径d50=0.028/0.78=0.036mm。对模型底沙，为保证模型小河的综合稳定性与原型相近，需按照河型相似条件进行选择，即将原型及模型的有关数据代入式(11)，得模型床沙中径D50=0.042mm。所对应的床沙粒径比尺D=0.18/0.042=4.29。为检查底沙的糙率，选用D50=0.040.06mm的煤粉作为模型沙进行预备试验，其动床河床糙率为0.0150.018，下限0.015为静平床糙率，上限0.018为河床出现沙纹后的动床糙率，故可满足阻力相似条件。为检查所选模型沙能否满足底沙起动及扬动相似条件，需分别确定模型沙及原型沙的起动流速和扬动流速。为此，采用煤粉进行了水槽试验，所得起动及扬动流速试验结果见表2。表2模型沙起动流速试验结果Theexperimentresultofthresholdvelocityofmodelsediment速度V(m/s)水深备注(cm)弱动普动扬动24.506.009.58sm=1.401.4245.507.5012.4066.508.5013.78d50=0.040.06mm87.209.2014.30为确定原型底沙的起动流速，点绘了黄河磴口、昭君坟、巴彦高勒及三湖河口水位站资料和引黄渠道沙土不冲流速与含沙量野外资料的关系(见图1)，以分析不冲流速与含沙量的关系，曲线与纵轴的交点，即为清水时的不冲流速。一般认为清水时的不冲流速等于起动流6，因而得到原型沙的起动流速Vcp=0.68m/s，相应水深范围内的模型沙起动流速Vcm=0.075-0.092m/s，可以得到Vc0.680.075-0.092=7.39-9.07下面采用沙玉清起动流速公式7计算天然河流床沙起动流速，该公式形式为图1不冲流速与含沙量关系RelationbetweennonscouringvelocityandsiltcontentVc=0.43d3/4+1.1(0.7-)4/d1/2h1/5(12a)式中为孔隙率，其稳定值约为0.4，d以mm计，Vc以m/s计。将d=0.18mm代入式（12a）得Vc=0.41h1/5(12b)当水深h=1-3.5m时，采用式(12b)及表2结果可求得相应的vc值列于表3。表3原型沙起动流速Vcp及起动流速比尺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所确定的起动流速比尺略小于流速比尺。根据以上情况综合分析可知，所选模型沙(包括底沙及悬沙中那部分与底沙经常交换的床沙质泥沙)基本能满足起动相似条件。黄河天然沙扬动流速一般为起动流速的1.541.75倍3，取原型底沙Vfp=1.65Vcp，则得到原型沙相应水深下的扬动流速为0.84m/s。由表1查得模型沙扬动流速Vfm=0.124-0.143m/s。从而Vf=6.77-5.87。所得扬动流速比尺与流速比尺较接近，表明所选模型沙还满足扬动相似条件(即vf=v)。从上述模型选沙结果来看，模型床沙中经常被水流冲起变为悬沙的那部分颗粒，其粒径接近模型沙的粒径，因而自然满足悬沙相似条件。3.3含沙量比尺及时间比尺式(9)是由河床变形方程式推导而来的，考虑了全沙变化对于河床变形的影响。由于原型实际上长期缺乏推移质泥沙观测资料，在试验中对推移质输沙进行定量模拟是相当困难的。不过，根据对原型情况分析后认为，由悬移质泥沙运动所引起的冲淤变化构成了该河段河床变形的主要部分，且推移质泥沙时常与悬移质泥沙交换，进入悬移质运动的行列中，因此本模型仅考虑悬移质输沙量变化的影响，并按照式(10)确定河床的冲淤变形时间比尺。事实上，通过河床验证试验确定模型进口加沙量时，自然反映了推移质运动对河床冲淤变形影响的相似性。另外，一般试验只模拟悬移质中的床沙质，但对于原型情况，悬移质中的冲泻质在河滩造床过程中起到了很大的作用。在试验中若相应扣除来沙量中的冲泻质，不仅使滩地淤积难以相似，而且还会使主槽变形产生偏离。因此本试验不再对悬移质中的床沙质和冲泻质加以区划。考虑到现有挟沙力公式不能同时应用于原型及模型，复杂河流的模型试验已不再直接利用由一些挟沙力公式直接导出的比尺关系式计算S*3。本试验分别确定原型和模型值，两者之比为含沙量比尺S*。经过试验河段附近实测资料验证后(见图2)，采用与原型资料颇为符合的黄科院张红武公式3，即S*=0.14(V3/ghlnh/6D50)0.6(13)图2水流挟沙力公式与原型资料对比Comparisonbetweenformulaandmeasureddata图3模型沙S*与V3/gh的关系RelationbetweenS*andV3/gh将水力泥沙因子代入上式，即V=1.7m/s,h=4.5m,R=4.5m,=0.22cm/s,D50=0.00018m,得S*p=5.26kg/m3利用模型沙进行挟沙力试验后，得到的结果如图3所示，点群遵循的关系可表示为S*m=0.18V3/gh(14)根据原型水力泥沙因子及相应的比尺，求得模型的Vm=24.05cm/s,Hm=9cm,=0.071cm/s，将其代入(14)式，可得S*m=4.0kg/m3。由于模型的挟沙能力小于水槽中的挟沙能力，参照以往黄河模型试验的经验，乘以折减系数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进一步改变水力及泥沙因子计算后可得到s=s*=1.5-2.0。将s及o代入(10)，得河床变形时间比尺t2=(1.08-1.44)t1=22.9-30.6。表明两个时间比尺相差较小，这对于非恒定动床模型试验很有益处，不仅能保证试验的可靠性，也便于试验操作。4模型验证对原型较为系统的两组同步水面线进行了施放,模型实测点与原型的偏差一般不超过2mm(相当于原型10cm),最大不超过3.2mm(相当于原型16cm)。在试验河段公路桥处有长期水位观测资料，将模型中相应的水位流量关系曲线与原型进行比较，可知两者颇为接近。对原型河段实测的5个流速断面进行验证，可以看到模型中流速沿河宽的分布与原型符合较好，表明满足了水流相似条件。按原型资料，对模型水沙过程进行概化。参照上述模型设计结果，分别选用不同的含沙量比尺及其相应的河床变形