黄河包头河段泥沙模型设计(何文社,方铎,刘有录,张淑英)

热度209票 浏览77次 【共0条评论】【我要评论 时间:2000年7月01日 16:05

黄河包头河段泥沙模型设计水利论文W A-YgE V g

何文社1,方 铎1,刘有录2,张淑英2
oE)uE6V\*zN0(1.四川大学 高速水力学国家重点试验室;2.兰州铁道学院)

yL:N!T3p%ru0

摘要:针对黄河泥沙多且变化大的特点,采用泥沙运动相似条件的最新研究成果,完成了黄河包头河段动床模型沙选择和比尺设计。利用原型实测资料进行的验证试验,表明模型小河较好地复演了原型河势、水沙运动规律及河床冲淤变形等。

%J7b9ci#X._J0

关键词:泥沙模型;模型设计;黄河

2hh!oiD0

基金项目:国家自然科学基金委员会和中国水利部联合资助重大项目(59890200)及国家自然科学基金项目(49771055)
8|K8l UUFMH {+h0作者简介:何文社(1966-),男,四川大学博士生。
水利论文 d2P:xzR8` {!n

1 模拟河段河道概况水利论文;h)tW Zn3t-C

  黄河包头段位于著名的河套地区。包头画匠营子河段位于包头市新市区略偏东处,距新市区约12km,在昭君坟和磴口之间(昭君坟水文站与磴口水文站相距约53km)。试验河段包括包头黄河公路桥下游1800m,两桥间800m,铁路桥上游3200m,共计约6000m。对昭君坟水文站1934年以来的水文资料整理分析后,按中水3000m3/s的流量考虑,得到原型河段代表性资料如表1。

-r&_,\ A}:o Z7~j0

表1 原型河道水流要素特征值

Cm9|z oh ob-dz0

Characteristic values of river flow in prototype

Vd0VvK4s;By0

项 目取 值项 目取 值

河段长度L6000m河段比降i0.1~0.2‰
河道宽度B1000m悬移质中径d500.028mm
主河道平均宽度b400m河床床沙中径D500.18mm
平均水深H4.5m最大流量Q1%6200m3/s
平均流速V1.7m/s最小流量Qmin48m3/s
河床糙率nb0.015~0.020泥沙比重γs2.65~2.70
河滩糙率nt0.020~0.030泥沙沉速ω0.22cm/s

2 河流泥沙模型相似条件

$t/["y k9}+VI$ab vI0

  目前黄河动床模型相似律的研究较多[1~4],并结合黄河包头河段的特性,采用下列相似比尺关系作为本模型设计的依据。水利论文 C,V J} U:v[1h

  重力相似λVh1/2(1)
  阻力相似λnH2/3L1/2(2)

  由水流连续方程,可分别得出流量及时间比尺关系式

'g\d.D(n(C0
λQVλLλH(3)
λt1LV(4)

  悬移质相似水利论文Qzs}m@\#i

λωVHL)3/4水利论文}?&Iu#VF6z

(5)水利论文*zhJ,k1z+o.y wb.s+?

对于细沙河流,当悬沙粒径d<0.15mm时(原形悬沙均能满足这一些条件),悬沙沉速用滞留区公式计算,由此得悬沙粒径比尺水利论文G NJ3l]~u\q&L

λd=(λωλVγs-γ)1/2水利论文 n`9WT7AR d/@#Fe

(6)水利论文za+dc @ V$UI

  底沙起动及扬动相似λVVcVf(7)
  水流输沙率相似λSS*(8)
  河床冲淤变形相似λt2λoλ2LλH/λGs(9)

  对于悬沙模型,上式又可表示为

0R,V6AoQD y;@ ]P0

λt2λosλt1水利论文B}K){GO k @XP4O

(10)水利论文 I!rS#x`{[6L

  河型相似水利论文.q3}^vQX

[(γs-γ/γD50H)1/3/iB2/3]m≈[(γs-γ/γD50H)1/3/iB2/3]p水利论文 wV#G@S

(11)水利论文$zS3kur0N

3 几何比尺及模型沙选取水利论文xd.|6p:i:]h+k

3.1 几何比尺及模型变率水利论文U(g*m'}n`@

  根据试验河段平面范围及试验场地条件,选取水平比尺λL=150。考虑到模型水流深度应满足表面张力及试验量测的要求,对模型几何变率前期的研究成果和模型沙特性等方面权衡后,取垂直比尺λH=50,则几何变率Dt=3。采用张红武[3]提出的变态模型相对保证率的公式求得Dt≤5.6;利用张瑞瑾[4]提出的模型变态限制指标求得Dt≤6.9。验算结果表明本模型采用的变率在各家公式限制的范围之内,几何变态的影响有限,可以满足工程需要。水利论文1f(w&IW}/G

3.2 模型沙选取

w/vj:}2G Q\,x0

  根据式(2)求得λn=1.11(取λRH),因而模型糙率nm=0.014~0.018,经比选后,决定采用经过专门加工处理的兰州煤厂精煤粉作为试验模型沙。该材料不仅比重小,而且密实度及凝聚力也较小(试验测得凝聚力C=0.018kg/m2,内摩擦角Φ=35°,容重γs=1.40~1.42kg/m3,干容重γ0=0.62~0.71kg/m3,含水量为60%,D50=0.03~0.05mm),因而能保证模型河床的活动性,以适应原型河床冲淤变化幅度较大的特点。

/a X$X~*fpc7M0

  原型沙容重γsp=2.66~2.70 kg/m3,干容重γop=1.40~1.50kg/m3。模型设计时,取模型沙γsm=1.41kg/m3, γom=0.67kg/m3;取原型沙γsp=2.68 kg/m3op=1.45 kg/m3。则水利论文eP$?\:W"Wu\w

λγs-γ=4.09水利论文/G2] jE-[J0Q

λγo=2.16水利论文e3CPs4d&v/j!@U&j6n3r q

  由式(5)得:λω=3.10;由式(6)得,λd=0.78(根据对原型及模型水流的温差分析,取λν=0.8)。

z$d!k{"?3Ml9YD0

  要求模型悬沙中径 d50=0.028/0.78=0.036mm。水利论文U+^u/dY6fKkb

  对模型底沙,为保证模型小河的综合稳定性与原型相近,需按照河型相似条件进行选择,即将原型及模型的有关数据代入式(11),得模型床沙中径D50=0.042mm。所对应的床沙粒径比尺λD=0.18/0.042=4.29。水利论文6^ n#b]6V `Db"\

  为检查底沙的糙率,选用D50=0.04~0.06mm的煤粉作为模型沙进行预备试验,其动床河床糙率为0.015~0.018,下限0.015为静平床糙率,上限0.018为河床出现沙纹后的动床糙率,故可满足阻力相似条件。

1J0y+z A3c2}0

  为检查所选模型沙能否满足底沙起动及扬动相似条件,需分别确定模型沙及原型沙的起动流速和扬动流速。为此,采用煤粉进行了水槽试验,所得起动及扬动流速试验结果见表2。

A"ys)P%TZ8i0

表2 模型沙起动流速试验结果水利论文!L&D){1n+`#? M

The experiment result of threshold velocity of model sediment

`x9Q'Mtf9^0

速度V(m/s)
水深
备注
(cm)弱动普动扬动

24.506.009.58γsm=1.40~1.42
45.507.5012.40
66.508.5013.78d50=0.04~0.06mm
87.209.2014.30

  为确定原型底沙的起动流速,点绘了黄河磴口、昭君坟、巴彦高勒及三湖河口水位站资料和引黄渠道沙土不冲流速与含沙量野外资料的关系(见图1),以分析不冲流速与含沙量的关系,曲线与纵轴的交点,即为清水时的不冲流速。一般认为清水时的不冲流速等于起动流速[6],因而得到原型沙的起动流速Vcp=0.68m/s,相应水深范围内的模型沙起动流速Vcm=0.075-0.092m/s,可以得到

λVc≈0.680.075-0.092=7.39-9.07

*R$g:P~"s'c0

  下面采用沙玉清起动流速公式[7]计算天然河流床沙起动流速,该公式形式为水利论文 X"EqIFQ6V

16t01.gif (5477 bytes)
图1 不冲流速与含沙量关系水利论文F6e+TH#A TBw
Relation between non scouring velocity水利论文^_%w%Z8B,f1s
and silt content

Vc=[0.43d3/4+1.1(0.7-ε)4/d]1/2h1/5

{)`.Bh oK8ZM0
(12a)

式中ε为孔隙率,其稳定值约为0.4,d以mm计,Vc以m/s计。将d=0.18mm代入式(12a)得水利论文kt6E`[vdn

Vc=0.41h1/5水利论文x]8M p A)w

(12b)

  当水深h=1-3.5m时,采用式(12b)及表2结果可求得相应的λvc值列于表3。水利论文1\"?(L4I)i#K

表3 原型沙起动流速Vcp及起动流速比尺λvc

;yer#E/y0

The values of threshold velocity Vcpand the velocity scale λvcof prototype sediment水利论文 g4OD9Hjgn


水深(m)1.01.52.03.03.5

Vcp(m/s)0.410.450.470.510.53

λvc6.836.626.276.005.96

备注λv=7.7

  由表3所确定的起动流速比尺略小于流速比尺。根据以上情况综合分析可知,所选模型沙(包括底沙及悬沙中那部分与底沙经常交换的床沙质泥沙)基本能满足起动相似条件。水利论文og o+bc*y+q$A8P

  黄河天然沙扬动流速一般为起动流速的1.54~1.75倍[3],取原型底沙Vfp=1.65Vcp,则得到原型沙相应水深下的扬动流速为0.84m/s。由表1查得模型沙扬动流速Vfm=0.124-0.143m/s。从而λVf=6.77-5.87。所得扬动流速比尺与流速比尺较接近,表明所选模型沙还满足扬动相似条件(即λvfv)。水利论文'l.q0L*ez*fQ_

  从上述模型选沙结果来看,模型床沙中经常被水流冲起变为悬沙的那部分颗粒,其粒径接近模型沙的粒径,因而自然满足悬沙相似条件。

R"R9b9HD6b~7V0

3.3 含沙量比尺及时间比尺水利论文*nEU5cb},F)y

  式(9)是由河床变形方程式推导而来的,考虑了全沙变化对于河床变形的影响。由于原型实际上长期缺乏推移质泥沙观测资料,在试验中对推移质输沙进行定量模拟是相当困难的。不过,根据对原型情况分析后认为,由悬移质泥沙运动所引起的冲淤变化构成了该河段河床变形的主要部分,且推移质泥沙时常与悬移质泥沙交换,进入悬移质运动的行列中,因此本模型仅考虑悬移质输沙量变化的影响,并按照式(10)确定河床的冲淤变形时间比尺。事实上,通过河床验证试验确定模型进口加沙量时,自然反映了推移质运动对河床冲淤变形影响的相似性。另外,一般试验只模拟悬移质中的床沙质,但对于原型情况,悬移质中的冲泻质在河滩造床过程中起到了很大的作用。在试验中若相应扣除来沙量中的冲泻质,不仅使滩地淤积难以相似,而且还会使主槽变形产生偏离。因此本试验不再对悬移质中的床沙质和冲泻质加以区划。

)N1oS8M.^i0

  考虑到现有挟沙力公式不能同时应用于原型及模型,复杂河流的模型试验已不再直接利用由一些挟沙力公式直接导出的比尺关系式计算λS*[3]。本试验分别确定原型和模型值,两者之比为含沙量比尺λS*水利论文G/_@ c9D^!hL&J

  经过试验河段附近实测资料验证后(见图2),采用与原型资料颇为符合的黄科院张红武公式[3],即水利论文$Bm2d|.hO#B/@g!O

S*=0.14(V3/ghωlnh/6D50)0.6

"[Hu/o"V0

(13)水利论文4TF[;pK7T

16t02.gif (4859 bytes)16t03.gif (4572 bytes)
图2 水流挟沙力公式与原型资料对比
8]MhQ2PX7W0Comparison between formula and measured data
图3 模型沙S*与V3/ghω的关系
*C3vj8So0Relation between S*and V3/ghω

将水力泥沙因子代入上式 ,即V=1.7m/s,h=4.5m,R=4.5m,ω=0.22cm/s,D50=0.00018m,得

5@)l \uv'o0

S*p=5.26kg/m3水利论文:bF/y-T;Nl4V

利用模型沙进行挟沙力试验后,得到的结果如图3所示,点群遵循的关系可表示为水利论文,J%_ai-{ vvB

S*m=0.18V3/ghω水利论文Xt-m{"mn_JV|$y0['p

(14)水利论文^+E3w q1rX4v

  根据原型水力泥沙因子及相应的比尺,求得模型的Vm=24.05cm/s,Hm=9cm,ω=0.071cm/s,将其代入(14)式,可得S*m=4.0kg/m3水利论文VEiem

  由于模型的挟沙能力小于水槽中的挟沙能力,参照以往黄河模型试验的经验,乘以折减系数K1[3](一般取K1=0.7-0.8)得水利论文v F)qT4@3s7_3EN.Nd

λs*=(0.7-0.8)4.0=2.8-3.2kg/m3

b1L:YF'ac8U9gb0

则 λs*=5.26/2.8-3.2=1.64-1.88

Y }$otuq~^J^0

  进一步改变水力及泥沙因子计算后可得到λss*=1.5-2.0。将λs及λγo代入(10),得河床变形时间比尺λt2=(1.08-1.44)λt1= 22.9-30.6。表明两个时间比尺相差较小,这对于非恒定动床模型试验很有益处,不仅能保证试验的可靠性,也便于试验操作。水利论文re'z!a K

4 模型验证

9U,k(Ee \wD0

对原型较为系统的两组同步水面线进行了施放,模型实测点与原型的偏差一般不超过2mm (相当于原型10cm),最大不超过3.2mm(相当于原型16cm)。在试验河段公路桥处有长期水位观测资料,将模型中相应的水位流量关系曲线与原型进行比较,可知两者颇为接近。

H7HA:x f @U2ay0

  对原型河段实测的5个流速断面进行验证,可以看到模型中流速沿河宽的分布与原型符合较好,表明满足了水流相似条件。

7zp!A} H X2T g&c0

  按原型资料,对模型水沙过程进行概化。参照上述模型设计结果,分别选用不同的含沙量比尺及其相应的河床变形时间比尺施放水沙过程。先后开展了4次验证试验,结果表明当含沙量比尺λS=2,相应的河床变形时间比尺λt2=23时,模型与原型相似良好。试验值不仅真实地复演了原型河势变化,而且由洪后实测的断面地形资料来看,模型河床冲淤变形与原型也颇为一致。水利论文0a#w%~ W6W VR%F

5 结论水利论文*\B UHZ6I+Z

  黄河系著名的多沙河流,原型河段又位于河套地区,含沙量大,且河床冲淤变化剧烈。为更真实地模拟原型天然河道中水流运动、泥沙运动以及河床演变的规律,在现有模型相似率的基础上,开展了大量的工作,采用多方案比选,对黄河包头河段模型设计进行了探索。通过该模型的设计与验证,可得出如下几点认识:水利论文q4zY"ij0])|r%Ui

  1.在确定模型悬沙粗度时,采用式(5)进行比尺计算,从理论上同时兼顾了时均流速输移及紊动扩散对悬沙相似的影响,特别是采用昭君坟、三湖河口、巴彦高勒及磴口等水文站大量资料,对来沙量及当地河段水流挟沙力的相互关系进行分析后,认为时均流速输移及紊动扩散对本河段的悬沙质都起着相当重要的作用。

"u@kKCjGxO a Z0

  2.在其它学者研究成果的基础上专门进行了模型沙特性试验。选取经过加工处理后的精煤粉做模型沙,具有比重适中、加工及分选方便的优点。不仅为本模型的设计提供了可靠的依据,且为模型验证及以后的生产试验提供了有利条件。水利论文t ng AiV;n'd

  3.模型不需人为加糙即能满足水面线相似要求,这对于模型制作和试验,都是很有意义的。水利论文B&j[9k%I!AEY)v

  4.模型设计以分别确定的原型及模型水流挟沙力两者之比作为含沙量比尺,这一作法相对于一般处理与实际更为接近。模型设计给出了λs及λs*的取值范围,减小了试验的盲目性。通过验证试验确定的河床冲淤变形时间比尺与水流运动时间比尺颇为接近,避免了常遇到的两个时间比尺相差过多所引起的一系列问题,这对于开展非恒定流河段动床模型具有重要的意义。

T mC-l'Q:od DO!c0

参 考 文 献

diq+SD"{"E0

[1] 张红武。论动床变态河工模型的相似律。黄科所科学研究论文集(第二集).河南科学技术出版社,1990.水利论文T'eCy_#yu,t/FX*L9]

[2] 李保如。我国河流泥沙物理模型的设计方法。水动力学研究与进展,1991(增刊).

T5Xf;zg2Jl0

[3] 张红武,江恩惠等。黄河高含沙洪水模型的相似律。河南科学技术出版社,1994.水利论文/F^"}O4]Z

[4] 屈孟浩。黄河动床模型试验相似原理及设计方法。黄科所科学研究论文集(第二集).河南科学技术出版社,1990.

V9}a6N1g/lLV0

[5] 张瑞瑾。论河道水流比尺模型变态问题。第二次河流泥沙国际学术讨论会论文集。南京。水利水电出版社,1983.

P$[y'lh9LW S Y0

[6] 徐正凡等。水力计算手册。水利出版社,1980.水利论文?@ v+c shqo

[7] 沙玉清。泥沙运动学引论。中国工业出版社,1965.水利论文$`&|L3}#N6r!u q!^4Q5h

 水利论文 uE*YO9L%~cO/K

TAG: 方铎 河段 黄河 包头 何文社
顶:22 踩:25
【已经有153人表态】
25票
极差
14票
很差
23票
较差
19票
稍差
18票
稍好
20票
较好
17票
很好
17票
极好
下一篇:水流含沙量对磨蚀的影响(邓军,杨永全,沈焕荣,许唯临,董建伟)
上一篇:论植物工程在西部环境生态建设中的地位和作用(拾兵,班久次仁,曹叔尤)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放