黄河下游引黄灌溉对策研究——论提水灌溉的作用与效果(史红玲,蒋如琴)

热度303票 浏览68次 【共0条评论】【我要评论 时间:2000年3月01日 15:39

黄河下游引黄灌溉对策研究

(@5B$kaP9|M0

——论提水灌溉的作用与效果

1_ HJ2\ x6Y0

史红玲蒋如琴水利论文OWN u!Ol3_

中国水利水电科学研究院)

.c9YK,d6KL0

摘要:在新的形势下,对黄河下游引黄灌溉中泥沙处理和利用应提高到一个新的认识高度,变泥沙“包袱”为“资源”。泥沙长距离输送,使绝大部分泥沙进入田间,就可实现引黄灌溉的良性循环。除了以往有效的措施外,本文在总结灌区渠道不同部位设置提水泵站设施作用与效果的同时,通过建立灌区渠道悬移质冲刷计算模型,对其减淤清淤的作用与效果进行了理论分析研究。

.MQogd^uV m0

关键词:长距离输送;提水灌溉水利论文^+_,?C8W'^ T Y0p

作者简介:史红玲(1973-),女,中国水科院硕士。

-|$Nerv7u4|0

1 泥沙是制约引黄灌溉发展的重要因素

1h'I P"f ]g(}9~0

  黄河下游河南、山东两省的引黄灌区,自1985年至1990年底,累计引水2333亿m3,引沙38.65亿t。其中,约有33.22%的泥沙沉积在沉沙池内,35.32%淤积在灌溉渠系,22.9%被输入田间,8.56%进入排水河道。这就是说,除了被输入田间的约23%的泥沙成为资源,其余近80%的泥沙需要进行处理。这种长期采取的泥沙处理方式(沉沙池集中沉沙,渠系两岸堆积清淤泥沙),造成了引黄泥沙在平面上分布不合理的状况,给沿黄两岸的经济、社会和环境带来了严重的负面影响。引黄泥沙已成为制约引黄灌溉发展的主要因素。

YE(T)k c] Poq8}0Q*W0

  随着黄河整个流域经济的迅速发展,用水量不断增加,而黄河下游水资源却日趋贫乏。引黄灌区面临需水期无水可引的局面,不得不引用洪水期高含沙水流。但目前,大多数灌区渠首无可供高含沙水流沉沙的空间,以进行沉沙输水。面对这一形势,对泥沙的处理和利用应有新的认识,这一认识的基本点是变泥沙“包袱”为“ 资源”。从技术上讲就是要实现泥沙的长距离输送。

c.y'Jk7|]v B3k0

  实现泥沙长距离输送,关键在于提高渠道输送泥沙的能力。要减少输沙渠道的泥沙淤积,除了以往行之有效的措施外,可采用改输沙渠水流自流入沉沙池变为扬水入沉沙池的方式,这种措施必将增加渠道的水面比降,加大流速,增加渠道内水流输送泥沙的能力。如若在渠首设置提水泵站扬水入渠,同样可起到加大渠道水面比降的作用。

/]w2u,\ Xdz4{:b0

  利用提水泵站实现泥沙长距离输送的可行性已在部分灌区得到证实。例如,渠首设置提水灌溉设施的夹马口灌区和尊村灌区及渠尾(中)设置提水泵站设施的曹店灌区和麻湾灌区。这些灌区在不同部位设置了扬水设施来加大水面比降,提高水流挟沙能力,达到了少淤或不淤的目的。

NL2Br/PG*l3D0

  本文在参照典型灌区冲淤变化状况的基础上,建立了适用于灌溉渠道输沙计算的数学模型,就提水灌溉的效果与作用进行了定性的研究与分析。水利论文1^A4IB)g2x

2 渠道悬移质冲刷模型的建立水利论文+ec%@\.sy]

2.1 基本方程的建立及数值解法水利论文R/H*`EC*`;Z\

  基于灌溉渠道水流泥沙运动规律与明渠水流泥沙运动规律一致,模型采用非耦合饱和输沙一维模型的计算方法,计算中假定渠道发生冲刷淤积过程,计算时段内河床变形对水流条件影响不大,也不考虑水体中含沙量随时间的变化,则该方程组为

n]Q+kDL i!q[0

  水流连续方程

3vE*bg%Tdw d0

5gVfU6qZ0

(1)水利论文?5S&Ewk.Dm;f

  水流运动方程

;G(^%P2Xz:EbR"S0

%jMXJ zpD0

(2)

!K%A3c#{ En z0

  河床变形方程

!?j+wjgLL2S:x0

水利论文S;~)U4W rJg3aWm

(3)

;{$OMl m(\+F0

  水流挟沙力公式水利论文e}S"@VJ-b

S*=S*(J,u,d)水利论文1C g/N5z.p

(4)水利论文p.F8fNx{,T&ap{l

式中 Q为流量;A为过水面积;B为水面宽;Y为水深;u为流速;Y0为断面渠底平均高程;S*为水流挟沙力;G为断面总输沙率;ρ为泥沙干密度;g为重力加速度;n为河床糙率;x,t为距离及时间。水利论文D:ySGC"{2h X\

  上述简化后的基本方程,即为建立模型理论依据。采用Pressman的四点偏心差格式,其中时间偏导采用前差,空间偏导采用加权中心差,可得函数和导数离散关系式水利论文(X9hG-r `'nkHf5D(n

水利论文/~2^ Z"]!vt*Lj U3`-Pp/z

式中θ为加权系数0≤θ≤1,f表示任意函数;上标n为时间坐标序列号,下标j为断面序号。

h;p.`8p6E9m6Za|0

  按上述差分格式推演水流连续方程(1)和水流运动方程(2)中各项,对计算点的未知量按泰勒级数展开,忽略增量Δfj的高次项,可使多元高次方程化为线性方程组。为简化计算方法,考虑渠道断面形态为梯形,经整理,对式(1)和式(2)分别有

#e ZOb9e-vU8g(f#f0

水利论文6@"nE,} P,WN V2O

(5)水利论文n|-}l*K;x,yI6fy

2FFm3X7nk _}S0

(6)水利论文iNE\8b7vW5^O

    以上系数A1,A2,A3,A4,A5及B1,B2,B3,B4,B5均为时段初相邻断面变量值Zjn,Qjn,,,及进库水流条件的函数。

  上述两个方程中,共有四个未知数,对M个计算点共有2M个未知数,由差分方程可列出2M-2个类似于(5)、(6)的方程组。因此需要补充两个边界条件可以方程组闭合。

*{{p Ka~ Z0

  使用交错网格追赶法求解差分方程。为将Q与Y分别布置在两个不同断面上,在断面间再增加一个节点(以K表示),并假定增加节点断面流量与上一个断面相等。水利论文)Q2E!C4r$`M

2.2 初始条件及边界条件

WJK&WHD-D6H0

  初始条件:Q(x)t=0=Q0(x),h(x)t=0=h0(x)水利论文q$c}+v {c&Qnvr

  上边界条件:Q1=Q(t)水利论文;M0LmHFJ+|

  下边界条件:Ym=Y(t)水利论文$p+hRf'\&Ga

2.3 河床变形方程及水流挟沙能力

$g7bt]"fW{ P@0

  河床变形方程式(3),采用有限差分法进行离散,得到差分方程,如下

v"n W{\+V|H}a Mh0

(G2-G1)ΔT=ρΔxBΔY0水利论文 N"Ww?t!RO$w

(7)

2w"e \G qq&n0

  式中 G1、G2分别为进出口断面输沙率;ρ为渠床泥沙干密度;ΔY0为计算河段平均河床淤积或冲刷厚度,正值为冲,负值为淤。对于悬移质变化引进的冲淤计算,取G=QS。水利论文OWK B5r

  水流挟沙能力采用模型相似类比法计算渠道断面挟沙力,由下式确定水利论文M'kD"d$k2A'|6gvu

S*=λS*×S0水利论文+R?x.ncpo-jx

(8)

il0t#|P `}0

其中λS*J7/4λu3/2d2水利论文7D#Q#l.W TTX8p

式中 J、u、d分别为水面比降、流速、悬移质中值粒径;λS*为挟沙能力比尺;S*0为原始渠道平衡挟沙力。

.Wv;ORm0

3 提水灌溉作用与效果的分析研究

pl)]/|eu'M}$Lep&N0

3.1 提水灌溉的作用水利论文 VEC:bz

  以曹店灌区干渠为模拟计算对象,取代主要参数为:设计引水流量30m3/s,梯形断面,渠底宽度为14m,边坡系数为1∶2,渠底比降为1/7000,渠道水深1.373m,断面平均流速1.342m/s,过水面积22.36m2。计算了下游泵站水位由7.373m下降至6.873m,即下降水位50cm,相应分流比约50%条件下渠道沿程含沙量的变化,见图1、图2。

TO+jc3E7zd0D(d0

02t01.gif (5785 bytes)

%NT!L^ {D0

02t02.gif (6661 bytes)

aU9`8Q}A@&b0

图1 降水过程中含沙量的变化

R,[SmI+Lt!\ {0

图2 降水过程结束后含沙量的变化

{J&wSI/Qq0

Sediment concentration variation along the canal during pumping water at the end of it

'iSS,vb0

The sediment concentration variation after pumping ended水利论文WA8j"}8h:TB0d

  由图可看出,在降水过程中,沿程含沙量随降水幅度的增加而增长;降水过程结束后,沿程含沙量又趋于平缓。经计算分析,沿程水面比降和平均流速也存在着这样的规律。这是由于泵站提水,造成局部水位下降,泵站的上渠段水面比降加大,导致水流流速增加,该渠段发生冲刷,出口含沙量迅速增加;当降水结束时,断面最大含沙量达到150kg/m3;降水结束后,随着冲刷的不断发展,下游水深增大,水面比降和流速逐渐减小,加上冲刷过程中悬移质的粗化,水流挟沙能力减小,制约了冲刷的进一步发展。

?$a'nq3GRBw,@0

  因此,建立提水泵站设施,可以对渠道起到一定的拉沙作用,但其拉沙效果则与降水幅度、降水历时及入流条件、渠道床面形态等有关。

O9E4~/d%J%Y0

3.2 提水灌溉效果的单因素分析

u,H$_:\aM$ou^0

  利用建立的模型,进行了包括降水模式(包括降水方式、幅度、历时)、入流量、渠道断面形态、糙率等诸因素对提水灌溉效果的计算,见表1、表2、表3、表4。水利论文 q/zGy"N#a+B,k

表1 降水模式对冲刷效果的影响水利论文?{;I IJ0N"Ak

The effect of erosion in different duration of pumping

I{7EB)Ql,V0

引水模式降水历时(h)降低出口水位值(cm)延伸距离(km)平衡历时(h)冲刷深度(cm)冲刷量(万m3)

分段降水2.0305.53.615.700.770
集中降水0.2307.55.011.800.760
分段降水2.0508.05.237.462.072
集中降水0.2509.05.035.402.297

  表1说明,长时段的分段降水和短时段的集中降水对冲刷效果的影响差别不大,而降水幅度对冲刷效果的影响则很大。

*N7r6TJ#A6dF k0

表2 渠道入流量对冲刷效果的影响水利论文)ffM"q7A0v

The effect of erosion for different flow discharge水利论文?%Z f _Fz5D


河宽(m)流量(m3/s)挟沙力(kg/m3)延伸距离(km)平衡历时(h)冲刷深度(cm)冲刷量(万m3)

分段降水50cm入流量对冲刷的影响142023.9368.838.771.871
14303085.237.462.072
144035.159.54.437.162.442
145039.6911438.982.775

集中降水50cm入流量对冲刷的影响142023.9386.834.092.028
1430309534.952.297
144035.15104.435.42.422
145039.6910437.762.585

  显然,渠道入流量对降水冲刷效果的影响比较大。流量越大冲刷越深、拉沙距离越远、所需平衡历时越少、冲刷量也越多;即越易冲刷。

C(Nx*q/]1B"@,@ H0

表3 不同宽深比B/h条件下的冲刷效果水利论文)@u z9nYzt

The effect of erosion for different value of/h of canal水利论文G7L,}"hz"]ll


渠底宽(m)水深(m)水面宽(m)

宽深比水利论文!} }7ACq*UV

/h水利论文j,m,b7PU yQ.A~3g

过水面积水利论文Z` v~,A2i8Q'Ox

(m2)

`8zRR W IUeG0

挟沙力S*水利论文h b"{"q%G6u

(kg/m3)水利论文7?P+ws~b,H_*[ \

延伸距离水利论文L0e,P5T*ZD

(km)水利论文/q{k6g?9J z8Bb

平衡历时

3\1`B x l1Gi$\0

(h)

4C8Yw+?%\.Bv*[|0

冲刷深度水利论文9[0yV3[)Pf-Q0U

(cm)水利论文[?E1N$g3q

冲刷量

M^:a dR&o? B0

(万m3)

IGcOmKjOv0

分段降水50cm不同宽深比条件下的冲刷效果101.59615.322.4520.2136.8585.638.482.762
121.47416.912.7921.3130.0085.237.462.072
141.37318.583.1422.3626.4884.836.781.735
181.21222.043.8724.2822.8774.835.941.379

集中降水50cm不同宽深比条件下的冲刷效果101.59615.322.4520.2136.859.55.434.13.115
121.47416.912.7921.3130.009.5534.92.297
141.37318.583.1422.3626.4894.835.41.836
181.21222.043.8724.2822.878.54.837.11.454

  从表3中可看出,宽深比越小冲刷效果越好,且有随宽深比的减小而递增的趋势。这说明了当下游条件改变时,渠道越宽浅冲刷效果越差。相反,窄深断面形态的渠道相对而言更易于冲刷。水利论文(V#Dy+E#XA

表4 降低水位50cm糙率的影响水利论文b{,V0Vu&j])~T

The effect of erosion for different roughness coeff.水利论文 fUn:tt%N5IEFa_


降水模式糙率延伸距离(km)平衡历时(h)平衡冲深(cm)平衡冲刷(万m3)

分段降水0.0118.56.631.91.818
分段降水0.0206.53.219.880.999
分段降水0.0305.02.612.040.572
集中降水0.0119.07.831.221.842
集中降水0.0209.02.417.961.025
集中降水0.0308.01.410.240.579

  从表4可看出,糙率越小冲刷效果越显著。反之,糙率越大冲刷效果越差,而且回淤也越严重。水利论文i/eG T"K1H$T"_r

  各表说明水位下降幅度是影响冲刷效果的决定性因素,入流量、断面形态及糙率等对冲刷效果有直接影响。因此,降低水位是否能真正形成自下而上的溯源冲刷而不产生回淤,有赖于各种条件的组合。

4}%aC/Bvu&Z$?0

4 主要结论水利论文0O OoSek{x

  1.泥沙的处理和利用直接关系到黄河下游引黄灌区地区的经济、社会和环境质量,成为制约引黄发展的重要因素。水利论文5H Vm3Pfv s

  2.在新的形势下应改变单纯的把泥沙当作“包袱”、集中力量进行处理这一思想,把泥沙和水一样作为一种资源。当泥沙用于淤背固堤,沉沙造地和输沙入田时泥沙就成为资源。关键在于实现泥沙长距离输送入田。

Fl^kE1s0

  3.泥沙长距离输送入田,渠道比降起着举足轻重的作用,当今应通过人类的能动作用,设法增加其它动力,如提水灌溉设施,并通过优化调度运行,千方百计加大渠道的纵比降。水利论文*k:D V7o2`s MX&T'Jo

  4.典型灌区的实践经验表明渠道不同部位(如渠首、 渠中、 渠尾)设置提水灌溉设施并优化调度运行,将大大提高渠道的水面比降,增大渠道水流挟沙能力,起到减淤的效果。

UC[T1`d;BS0

  5.通过建立灌区渠道悬移质冲刷计算模型对提水灌溉作用和效果的理论分析研究表明,利用提水灌溉设施,可以通过降低其水位(或低水位运行)达到加大局部渠段的水面比降、流速和水流挟沙力的目的,在一定的条件下(主要是水位下降幅度),在渠道中可以产生自下而上的溯源冲刷,起到拉沙和清除淤积泥沙的作用。水利论文l|G'Zm e

参 考 文 献

yzW ] j$p~0

[1] 蒋如琴等。引黄渠系泥沙利用及对平原排水影响的研究。中国水利水电科学研究院,1996年12月。

w`,B4n \u'M.X0

[2] 吕秀贞,彭润泽。悬沙冲淤计算——隐式有限差分追赶法求解。泥沙研究,1985,(1).

Y7eK&i)uh(}0

[3] 武汉水利电力大学河流模拟研究室。黄河一维泥沙数学模型研究。1993年12月。

B"@$h VB*e0

[4] 张丽春,方红卫,府仁寿。一维非恒定非均匀泥沙数学模型研究。泥沙研究,1998,(3).

8V3S5h{\]0

[5] 李义天,尚全民。一维不恒定流泥沙数学模型研究。泥沙研究,1998,(1).[6] 王兆印,徐永年,苏晓波。挟沙水流的冲刷率及河床惯性的研究。泥沙研究,1998,(2).

3h oixwFH y0

[7] 蒋如琴,戴清。黄河中、下游引黄灌区渠道水流挟沙能力的研究。泥沙研究,1997,(2).水利论文*~5aLSa\.F

[8] 杨国录。四点时空偏心Preissmann格式的应用问题。泥沙研究,1991,(4).

:z"v?;B iO4h0

[9] 钱宁,万兆惠。泥沙运动力学。科学出版社,1983.

{N$d@.^0

[10] 曹祖德,王运洪。水动力泥沙数值模拟。天津大学出版社,1994年11月。水利论文Y~N0LZ;b

[11] 中国水利学会泥沙专业委员会主编。泥沙手册。中国环境出版社,1992.水利论文;I&Q#ug$i9r

"}*pb sfb/@GLP0[12]K.Mahmood and V.Yevievich Unsteady Flow in Open Channels 1975.

4TQ Xn`m n3h$S1e0水利论文|| gM.YT f

 

7E(F{:y1m0
TAG: 黄河 蒋如琴 下游 史红玲 提水灌溉
顶:40 踩:42
【已经有221人表态】
23票
极差
28票
很差
25票
较差
33票
稍差
34票
稍好
25票
较好
28票
很好
25票
极好
下一篇:黄河下游水旱灾害治理措施的研究(何富荣,胡春宏)
上一篇:黄河下游引黄灌溉中的泥沙处理利用(程秀文,尚红霞,陈发中,钱意颖)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放